IP Library Granted Patent US 12,509,975
Granted Patent B2
US 12,509,975 · App. 17/537,901 · Granted Dec 30, 2025

Systems and methods for hydraulic fracture and reservoir simulation

Inventor: Mark W. McClure (Palo Alto, CA)
Assignee: ResFrac Corporation
E21B43/26E21B47/0025E21B49/008E21B2200/20
View Patent ↗
Loading inventors, assignments & file history…
Monitor This Case
Get email alerts when status or documents change.
Order Certified Copies
Most orders are placed with the USPTO same day — all within 24 business hours.
Order via The Patent Place →
Pre-filled with this patent's details
Quick Facts
Patent No.
US 12,509,975
App. No.
17/537,901
Granted
Dec 30, 2025
Kind
B2
Abstract

Provided herein are systems and methods for modeling and simulating reservoir, wellbore, and hydraulic fracturing. The systems and methods provided herein may facilitate well life cycle simulation by integrating a three-dimensional model representative of hydraulic fracturing and fluid flow in a wellbore and reservoir. The systems and methods may couple fluid flow in the wellbore and reservoir during injection and extraction with propagation of fractures through subsurface materials during fluid injection. Integrated three-dimensional reservoir, wellbore, and hydraulic fracture simulation may be useful for the design of hydraulic fracture treatments and prediction of future reservoir production.

Claims (25)

1 . A method for determining hydraulic fracture treatment of a production well, comprising:

(a) receiving one or more input parameters for a production well, wherein at least one of the one or more input parameters are based on data from an integrated three-dimensional (3-D) model of a calibration well;

(b) simultaneously simulating crack propagation and fluid transport through simulated cracks in the production well by performing an integrated 3-D simulation of hydraulic fracturing and fluid flow in a wellbore and reservoir of the production well based on the one or more input parameters for the production well, wherein performing the integrated 3-D simulation comprises:

i) creating an integrated 3-D model representative of the hydraulic fracturing and the fluid flow in the wellbore and the reservoir of the production well;

ii) adaptively treating a component of a matrix element, a fracture element, and a wellbore element of the integrated 3-D model as an explicit component or an implicit component with respect to a time step to simultaneously solve the matrix element, the fracture element, and the wellbore element during the crack propagation and fluid transport; and

(c) outputting the hydraulic fracture treatment of the production well based on the integrated 3-D simulation of hydraulic fracturing and fluid flow in the wellbore and reservoir of the production well.

2 . The method of claim 1 , wherein the component of the matrix element, the fracture element, and the wellbore element of the integrated 3-D model is adaptively treated as the explicit component or the implicit component based on numerical instability of the component.

3 . The method of claim 1 , wherein the one or more input parameters comprise one or more hydraulic fracture treatment conditions of the production well.

4 . The method of claim 3 , wherein the one or more hydraulic fracture treatment conditions of the production well include (i) spacing of perforation clusters, (ii) spacing between wells, (iii) amount of proppant injected into a perforation cluster, (iv) injection rate, (v) injection volume, (vi) length of each stage along the production well, (vii) type of proppant injected, (viii) type of fluid injected, (ix) sequencing of fluid and proppant injection during a stage, or (x) sequencing of injection stages.

5 . The method of claim 3 , wherein the one or more hydraulic fracture treatment conditions of the production well are provided in an input file.

6 . The method of claim 5 , wherein the input file is generated by a user with the assistance of a graphical user interface.

7 . The method of claim 1 , wherein the calibration well is disposed adjacent to or in the same geological formation as the production well.

8 . The method of claim 1 , wherein the integrated 3-D model of the calibration well is created by adaptively treating a component of the integrated 3-D model of the calibration well as explicit or implicit with respect to a time step.

9 . The method of claim 1 , wherein the one or more input parameters are provided in an input file.

10 . The method of claim 1 , wherein the one or more input parameters comprise a model parameter modified based on the data from the integrated 3-D model of a calibration well.

11 . The method of claim 10 , wherein the model parameter is modified based on a comparison between the model parameter and a production data from the integrated 3-D model of a calibration well.

12 . The method of claim 11 , wherein the production data includes production rate, production pressure, injection pressure during fracturing, or fracture length.

13 . The method of claim 1 , wherein the integrated 3-D model representative of the hydraulic fracturing and the fluid flow in the wellbore and the reservoir of the production well provides a sensitivity analysis for the geological data.

14 . The method of claim 13 , wherein the geological data includes permeability or fracture conductivity.

15 . The method of claim 1 , wherein the integrated 3-D model representative of the hydraulic fracturing and the fluid flow in the wellbore and the reservoir of the production well assesses which of the one or more input parameters have the largest impact on performance of the integrated 3-D simulation.

16 . The method of claim 1 , further comprising displaying to a user via a graphical user interface one or more output properties representing the response and state of the production well at a given time.

17 . The method of claim 16 , wherein the one or more output properties are selected from the group consisting of fluid pressure, temperature, fluid saturation, molar composition, fluid phase density, fluid phase viscosity, proppant volume fraction, and fracture aperture.

18 . The method of claim 1 , wherein the one or more input parameters comprise production boundary conditions.

19 . The method of claim 1 , further comprising solving the explicit component to obtain explicit variables, using the explicit component and the explicit variables to solve the implicit component, and using numerical differentiation to approximate derivatives of the implicit component to obtain implicit variables.

20 . The method of claim 19 , further comprising evaluating the implicit component after substitution of the explicit component into the implicit component.

Assignments (2)
SECURITY INTEREST Recorded Aug 27, 2026
From: RESFRAC CORPORATION
To: BANC OF CALIFORNIA
Reel/Frame 075802/0534 →
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Oct 3, 2022
From: MCCLURE, MARK W.
To: RESFRAC CORPORATION
Reel/Frame 061289/0960 →
Continuity (5)
Continuation In Part 16991727 · Aug 12, 2020
Continuation In Part 16359498 · Mar 20, 2019
Continuation In Part PCTUS2019023216 · Mar 20, 2019
Provisional Application 62646150 · Mar 21, 2018
Related Publication 20220127940A1 · Apr 28, 2022
References Cited (62)
US 9540911B2 · Stone et al. · 2017 [cited by applicant]
US 10954757B2 · Shetty et al. · 2021 [cited by applicant]
US 11220889B2 · McClure · 2022 [cited by applicant]
US 11225855B2 · McClure · 2022 [cited by applicant]
US 20040199329A1 · Stone · 2004 [cited by applicant]
US 20090164189A1 · Bourbiaux · 2009 [cited by applicant]
US 20110024125A1 · Wallace et al. · 2011 [cited by applicant]
US 20110257944A1 · Du · 2011 [cited by examiner]
US 20130073268A1 · Abacioglu et al. · 2013 [cited by applicant]
US 20130238304A1 · Glinsky · 2013 [cited by applicant]
US 20130346035A1 · Madasu et al. · 2013 [cited by applicant]
US 20140151033A1 · Xu · 2014 [cited by applicant]
US 20140222403A1 · Lepage et al. · 2014 [cited by applicant]
US 20140297235A1 · Arora et al. · 2014 [cited by applicant]
US 20140372089A1 · Weng et al. · 2014 [cited by applicant]
US 20160010443A1 · Xu · 2016 [cited by applicant]
US 20160123118A1 · Downton et al. · 2016 [cited by applicant]
US 20160283651A1 · Knight et al. · 2016 [cited by applicant]
US 20170002652A1 · Kampfer et al. · 2017 [cited by applicant]
US 20170030166A1 · Shen et al. · 2017 [cited by applicant]
US 20170275970A1 · Crawford et al. · 2017 [cited by applicant]
US 20170315266A1 · Myers et al. · 2017 [cited by applicant]
US 20170328179A1 · Dykstra et al. · 2017 [cited by applicant]
US 20180010429A1 · Willberg et al. · 2018 [cited by applicant]
US 20180058212A1 · Ali et al. · 2018 [cited by applicant]
US 20180252102A1 · Anschutz et al. · 2018 [cited by applicant]
US 20190292884A1 · Mcclure · 2019 [cited by applicant]
US 20200032623A1 · McClure · 2020 [cited by applicant]
US 20200371262A1 · McClure · 2020 [cited by applicant]
CA 2240580A1 · 1999 [cited by examiner]
CN 102575510A · 2012 [cited by applicant]
CN 105452599A · 2016 [cited by applicant]
CN 106844909A · 2017 [cited by applicant]
EP 3101225A1 · 2016 [cited by applicant]
FR 3000579A1 · 2014 [cited by applicant]
FR 3036210A1 · 2016 [cited by applicant]
WO WO2011081544A1 · 2011 [cited by applicant]
WO WO2016032489A1 · 2016 [cited by applicant]
WO WO2016183043A1 · 2016 [cited by applicant]
WO WO2017059152A1 · 2017 [cited by applicant]
WO WO2017078671A1 · 2017 [cited by applicant]
WO WO2017078674A1 · 2017 [cited by applicant]
WO WO2017082862A1 · 2017 [cited by applicant]
WO WO2017216594A1 · 2017 [cited by applicant]
WO WO2018033524A1 · 2018 [cited by applicant]
WO WO2018165352A1 · 2018 [cited by examiner]
WO WO2019183252A1 · 2019 [cited by applicant]
Oballa_1990 (Adaptive-Implicit Method in Thermal Simulation, SPE Reservoir Engineering, Nov. 1990). (Year: 1990). [cited by examiner]
Mcclure et al. A three-dimensional reservoir, wellbore, andhydraulic fracturing simulator that is compositional and thermal, tracksproppant and water solute transport, includes non-darcy and non-newtonian flow, and hand… [cited by applicant]
Mcclure et al. ResFrac technical writeup. pp. 1-54. (Latestversion—2019). [cited by applicant]
PCT/US2019/023216 International Search Report dated Jul. 11, 2019. [cited by applicant]
US 20040199329 Al, a compressed version, Oct. 2004. [cited by applicant]
US 20110024125 A1, a compressed version, Feb. 2011. [cited by applicant]
U.S. Appl. No. 16/359,498 Notice of Allowance dated Sep. 27, 2021. [cited by applicant]
U.S. Appl. No. 16/359,498 Office Action dated Jun. 30, 2021. [cited by applicant]
U.S. Appl. No. 16/507,638 Notice of Allowance dated Oct. 20, 2021. [cited by applicant]
U.S. Appl. No. 16/507,638 Office Action dated Jul. 2, 2021. [cited by applicant]
U.S. Appl. No. 16/991,727 Office Action dated Sep. 14, 2021. [cited by applicant]
U.S. Appl. No. 16/359,498 Notice of Allowance dated Nov. 10, 2021. [cited by applicant]
U.S. Appl. No. 16/359,498 Notice of Allowance dated Nov. 23, 2021. [cited by applicant]
U.S. Appl. No. 16/507,638 Notice of Allowance dated Nov. 10, 2021. [cited by applicant]
EP19770871.2 Extended European Search Report dated Nov. 26, 2021. [cited by applicant]