OPTIMIZING MULTISTAGE HYDRAULIC FRACTURING DESIGN BASED ON THREE-DIMENSIONAL (3D) CONTINUUM DAMAGE MECHANICS

• US 20170169137A1
• Filed: 08/28/2014
• Published: 06/15/2017
• Est. Priority Date: 08/28/2014
• Status: Active Grant

• Current Assignee: Landmark Graphics Corporation (Halliburton Company)
• Original Assignee: Landmark Graphics Corporation (Halliburton Company)
• Inventors: Shen, Xinpu, Standifird, William B., Evans, Scot
• Granted Patent: US 9,836,561 B2
• Time in Patent Office: Days
• CPC Class Codes: E21B 41/00 Equipment or details not covered by groups E21B15/00 - E21B40/00
• CPC Class Codes: E21B 43/26 by forming crevices or fractures
• CPC Class Codes: E21B 49/00 Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
• CPC Class Codes: E21B 7/046 horizontal drilling drilling with mechanical conveying means E21B7/003
• CPC Class Codes: G01V 1/282 Application of seismic models, synthetic seismograms
• CPC Class Codes: G01V 20/00 Geomodelling in general
• CPC Class Codes: G01V 2210/646 Fractures
• CPC Class Codes: G06F 2113/08 Fluids
• CPC Class Codes: G06F 30/20 Design optimisation, verification or simulation optimisation, verification or simulation of circuit designs G06F30/30
• CPC Class Codes: G06F 30/23 using finite element methods [FEM] or finite difference methods [FDM]

Abstract

System and methods for optimizing multistage hydraulic fracturing design based on three-dimensional (3D) continuum damage mechanics are provided. A 3D global model of a field of hydrocarbon reservoirs is generated. A finite element analysis of the 3D global model is performed to calculate values of one or more mechanical variables for the field. A smaller-scale 3D sub-model of a selected portion of the field is generated from the 3D global model based on the calculated values of the mechanical variables. The hydraulic fracturing effects of fluid injection stimulation within the targeted reservoir formation are simulated by applying one or more numerical damage models to the sub-model. Optimal design parameters for multistage hydraulic fracturing of the targeted formation with a plurality of horizontal wells are determined based on numerical results of the simulation.

Claims

1. A computer-implemented method for optimizing multistage hydraulic fracturing design based on three-dimensional (3D) continuum damage mechanics, the method comprising:

generating a 3D global model of a field of hydrocarbon reservoirs;

calculating values of one or more mechanical variables for the field based on a finite element analysis of the generated 3D global model;

generating a smaller-scale 3D sub-model of a selected portion of the field based on the 3D global model and the calculated values of the one or more mechanical variables, the selected portion corresponding to a formation of at least one of the field'"'"'s hydrocarbon reservoirs that is targeted for fluid injection stimulation;

applying one or more numerical damage models to the generated 3D sub-model to simulate hydraulic fracturing effects of the fluid injection stimulation within the targeted reservoir formation; and

determining one or more optimal design parameters for multistage hydraulic fracturing of the targeted reservoir formation with a plurality of horizontal wells, based on numerical results of the simulation.

2. The method of claim 1, wherein the hydraulic fracturing effects are simulated in vertical, transverse, and horizontal directions within a 3D volume of the 3D sub-model.

3. The method of claim 1, wherein the one or more numerical damage models applied to the 3D sub-model include a plasticity-based continuum damage model applied to the 3D sub-model to stimulate the hydraulic fracturing effects at a material level of the targeted reservoir formation and a coupled poro-elastoplastic finite element model applied to the 3D sub-model to simulate the hydraulic fracturing effects at a structural level of the targeted reservoir formation as a whole.

4. The method of claim 1, wherein the one or more mechanical variables include a geo-stress distribution, a pore pressure distribution, and a displacement distribution across one or more external surfaces of the 3D global model.

5. The method of claim 4, wherein the one or more external surfaces include a bottom surface and one or more lateral surfaces of the 3D global model.

6. The method of claim 4, wherein the values of the one or more mechanical variables calculated using the 3D global model are applied to the 3D sub-model as a set of boundary and initial conditions for simulating the hydraulic fracturing effects of the fluid injection stimulation within the targeted reservoir formation.

7. The method of claim 6, wherein the numerical results of the simulation include continuum damage distribution values that are calculated based on the set of boundary and initial conditions applied to the 3D sub-model, the continuum damage distribution values representing the hydraulic fracturing effects of the fluid injection stimulation within the targeted reservoir formation relative to a fluid injection point.

8. The method of claim 1, wherein the plurality of horizontal wells are drilled using a zipper frac technique involving two or more parallel horizontal wells, the multistage hydraulic fracturing includes performing a multistage stimulation along each of the two or more parallel horizontal wells, and the one or more optimal design parameters include an optimal stage interval between successive stages of the multistage stimulation in addition to an optimal well spacing between the two or more parallel horizontal wells.

9. The method of claim 8, wherein the 3D global model is generated based on subsurface data captured for the field, the subsurface data including microseismic data obtained from stimulation of one or more offset wells drilled for the field.

10. The method of claim 9, further comprising:

determining a critical value for a synthetic continuum damage variable of the targeted reservoir formation based on information derived from the microseismic data, the synthetic continuum damage variable representing the hydraulic fracturing effects resulting from compressive and tensile stresses under the fluid injection stimulation, the critical value corresponding to a value of the synthetic continuum damage variable that is sufficient to form an effective fracture within a zone of fracture propagation under the fluid injection stimulation, and the effective fracture having an opening of at least a predetermined size for holding a propping agent used in the fluid injection stimulation; and

determining the optimal stage interval and the optimal well spacing based on the critical value determined for the synthetic continuum damage variable.

11. The method of claim 10, wherein the optimal stage interval is determined based on a distance of fracture propagation along a first path between a fluid injection point and a point on the first path corresponding to half of the critical value of the synthetic continuum damage variable, the first path extends in an axial direction relative to a trajectory of the two or more parallel horizontal wells, and the point on the first path represents a location between adjacent stimulation stages corresponding to a same one of the two or more parallel horizontal wells where a cumulative effect of overlapping fracture propagation zones formed by the adjacent stimulation stages allows the synthetic continuum damage variable to reach the critical value and thereby increase a total area of the effective fracture.

12. The method of claim 11, wherein the optimal well spacing is determined based on a distance of fracture propagation along a second path between the fluid injection point and a point on the second path corresponding to half of the critical value, the second path extends in a perpendicular direction relative to the first path and the trajectory of the two or more parallel horizontal wells, and the point on the second path represents a location between at least two different ones of the two or more parallel horizontal wells where the cumulative effect of overlapping fracture propagation zones formed by stimulation stages corresponding to the at least two different parallel horizontal wells further increases the total area of the effective fracture.

13. A system for optimizing multistage hydraulic fracturing design based on three-dimensional (3D) continuum damage mechanics, the system comprising:

at least one processor; and

a memory including instructions stored therein, which when executed by the processor, cause the processor to perform functions including functions to;

generate a 3D global model of a field of hydrocarbon reservoirs;

calculate values of one or more mechanical variables for the field based on a finite element analysis of the generated 3D global model, the one or more mechanical variables including a geo-stress distribution, a pore pressure distribution, and a displacement distribution across one or more external surfaces of the 3D global model;

generate a smaller-scale 3D sub-model of a selected portion of the field based on the 3D global model and the calculated values of the one or more mechanical variables, the selected portion corresponding to a formation of at least one of the field'"'"'s hydrocarbon reservoirs that is targeted for fluid injection stimulation;

apply one or more numerical damage models to the generated 3D sub-model to simulate hydraulic fracturing effects of the fluid injection stimulation within the targeted reservoir formation, wherein the hydraulic fracturing effects are simulated in vertical, transverse, and horizontal directions within a 3D volume of the 3D sub-model; and

determine one or more optimal design parameters for multistage hydraulic fracturing of the targeted reservoir formation with a plurality of horizontal wells, based on numerical results of the simulation.

14. The system of claim 13, wherein the one or more numerical damage models applied to the 3D sub-model include a plasticity-based continuum damage model applied to the 3D sub-model to stimulate the hydraulic fracturing effects at a material level of the targeted reservoir formation and a coupled poro-elastoplastic finite element model applied to the 3D sub-model to simulate the hydraulic fracturing effects at a structural level of the targeted reservoir formation as a whole.

15. The system of claim 13, wherein the values of the one or more mechanical variables calculated using the 3D global model are applied to the 3D sub-model as a set of boundary and initial conditions for simulating the hydraulic fracturing effects of the fluid injection stimulation within the targeted reservoir formation.

16. The system of claim 15, wherein the numerical results of the simulation include continuum damage distribution values that are calculated based on the set of boundary and initial conditions applied to the 3D sub-model, the continuum damage distribution values representing the hydraulic fracturing effects of the fluid injection stimulation within the targeted reservoir formation relative to a fluid injection point.

17. The system of claim 13, wherein the plurality of horizontal wells are drilled using a zipper frac technique involving two or more parallel horizontal wells, the multistage hydraulic fracturing includes performing a multistage stimulation along each of the two or more parallel horizontal wells, and the one or more optimal design parameters include an optimal stage interval between successive stages of the multistage stimulation in addition to an optimal well spacing between the two or more parallel horizontal wells.

18. The system of claim 17, wherein the functions performed by the processor further include functions to:

determine a critical value for a synthetic continuum damage variable of the targeted reservoir formation based on information derived from the microseismic data, the synthetic continuum damage variable representing the hydraulic fracturing effects resulting from compressive and tensile stresses under the fluid injection stimulation, the critical value corresponding to a value of the synthetic continuum damage variable that is sufficient to form an effective fracture within a zone of fracture propagation under the fluid injection stimulation, and the effective fracture having an opening of at least a predetermined size for holding a propping agent used in the fluid injection stimulation; and

determine the optimal stage interval and the optimal well spacing based on the critical value determined for the synthetic continuum damage variable.

19. The system of claim 18, wherein:

the optimal stage interval is determined based on a distance of fracture propagation along a first path between a fluid injection point and a point on the first path corresponding to half of the critical value of the synthetic continuum damage variable, the first path extending in an axial direction relative to a trajectory of the two or more parallel horizontal wells, the point on the first path representing a location between adjacent stimulation stages corresponding to a same one of the two or more parallel horizontal wells where a cumulative effect of overlapping fracture propagation zones formed by the adjacent stimulation stages allows the synthetic continuum damage variable to reach the critical value and thereby increase a total area of the effective fracture; and

the optimal well spacing is determined based on a distance of fracture propagation along a second path between the fluid injection point and a point on the second path corresponding to half of the critical value, the second path extending in a perpendicular direction relative to the first path and the trajectory of the two or more parallel horizontal wells, and the point on the second path representing a location between at least two different ones of the two or more parallel horizontal wells where the cumulative effect of overlapping fracture propagation zones formed by stimulation stages corresponding to the at least two different parallel horizontal wells further increases the total area of the effective fracture.

20. A computer-readable storage medium having instructions stored therein, which when executed by a processor cause the processor to execute functions, including functions to:

generate a 3D global model of a field of hydrocarbon reservoirs;

calculate values of one or more mechanical variables for the field based on a finite element analysis of the generated 3D global model;

generate a smaller-scale 3D sub-model of a selected portion of the field based on the 3D global model and the calculated values of the one or more mechanical variables, the selected portion corresponding to a formation of at least one of the field'"'"'s hydrocarbon reservoirs that is targeted for fluid injection stimulation;

apply one or more numerical damage models to the generated 3D sub-model to simulate hydraulic fracturing effects of the fluid injection stimulation within the targeted reservoir formation; and

determine one or more optimal design parameters for multistage hydraulic fracturing of the targeted reservoir formation with a plurality of horizontal wells, based on numerical results of the simulation.