Research Papers: Oil/Gas Reservoirs

Modeling of Hydraulic Fracture Propagation in Shale Gas Reservoirs: A Three-Dimensional, Two-Phase Model

[+] Author and Article Information
Chong Hyun Ahn

Petroleum and Natural Gas Engineering,
Department of Energy and Mineral Engineering
and EMS Energy Institute,
The Pennsylvania State University,
3S Laboratory for Petroleum Research,
202 Hosler Building,
University Park, PA 16802
e-mail: cza5010@psu.edu

Robert Dilmore

U.S. Department of Energy,
National Energy Technology Laboratory,
626 Cochrans Mill Road,
P.O. Box 10940,
Pittsburgh, PA 15236-0940
e-mail: robert.dilmore@netl.doe.gov

John Yilin Wang

Petroleum and Natural Gas Engineering,
Department of Energy and Mineral Engineering
and EMS Energy Institute,
The Pennsylvania State University,
3S Laboratory for Petroleum Research,
202 Hosler Building,
University Park, PA 16802
e-mail: john.wang@psu.edu

Contributed by the Petroleum Division of ASME for publication in the JOURNAL OF ENERGY RESOURCES TECHNOLOGY. Manuscript received January 1, 2016; final manuscript received May 21, 2016; published online June 27, 2016. Assoc. Editor: Egidio Marotta.This material is declared a work of the U.S. Government and is not subject to copyright protection in the United States. Approved for public release; distribution is unlimited.

J. Energy Resour. Technol 139(1), 012903 (Jun 27, 2016) (13 pages) Paper No: JERT-16-1001; doi: 10.1115/1.4033856 History: Received January 01, 2016; Revised May 21, 2016

A three-dimensional, two-phase, dual-continuum hydraulic fracture (HF) propagation simulator was developed and implemented. This paper presents a detailed method for efficient and effective modeling of the fluid flow within fracture and matrix as well as fluid leakoff, fracture height growth, and the fracture network propagation. Both a method for solving the system of coupled equations, and a verification of the developed model are presented herein.

Copyright © 2017 by ASME
Your Session has timed out. Please sign back in to continue.


Ahn, C. H. , Chang, O. C. , Dilmore, R. , and Wang, J. Y. , 2014, “ A Hydraulic Fracture Network Propagation Model in Shale Gas Reservoirs: Parametric Studies to Enhance the Effectiveness of Stimulation,” SPE Unconventional Resources Technology Conference, Denver, CO, Aug. 25–27, Paper No. SPE 1922580-MS.
Cipolla, C. L. , Warpinski, N. R. , Mayerhofer, M. J. , and Lolon, E. P. , 2010, “ The Relationship Between Fracture Complexity, Reservoir Properties, and Fracture-Treatment Design,” SPE Prod. Oper., 25(4), pp. 438–452. [CrossRef]
Li, M. , and Lior, N. , 2015, “ Analysis of Hydraulic Fracturing and Reservoir Performance in Enhanced Geothermal Systems,” ASME J. Energy Resour. Technol., 137(4), p. 041203. [CrossRef]
Lorenz, J. C. , and Warpinski, N. R. , 1996, “ Natural Fracture Characteristics and Effects,” Leading Edge, 15(8), pp. 909–911. [CrossRef]
Osholake, T. , Wang, J. Y. , and Ertekin, T. , 2012, “ Factors Affecting Hydraulically Fractured Well Performance in the Marcellus Shale Gas Reservoirs,” ASME J. Energy Resour. Technol., 135(1), p. 013402. [CrossRef]
Rahman, M. , and Rahman, M. , 2011, “ Optimizing Hydraulic Fracture to Manage Sand Production by Predicting Critical Drawdown Pressure in Gas Well,” ASME J. Energy Resour. Technol., 134(1), p. 013101. [CrossRef]
Teufel, L. W. , and Clark, J. A. , 1984, “ Hydraulic Fracture Propagation in Layered Rock: Experimental Studies of Fracture Containment,” SPE J., 24(1), pp. 19–32. [CrossRef]
Blanton, T. L. , 1982, “ Experimental Study of Interaction Between Hydraulically Induced and Pre-Existing Fractures,” SPE Unconventional Gas Recovery Symposium, Pittsburgh, PA, May 16–18, Paper No. SPE 10847.
Renshaw, C. E. , and Pollard, D. D. , 1998, “ An Experimetally Verified Criterion for Propagation Across Unbounded Frictional Interfaces in Brittle Linear Elastic Materials,” Int. J. Rock Mech. Min. Sci. Geomech., 32(3), pp. 237–249. [CrossRef]
Warpinski, N. R. , and Teufel, L. W. , 1987, “ Influence of Geologic Discontinuities on Hydrauli Fracture Propagation,” J. Pet. Technol., 39(2), pp. 209–220. [CrossRef]
Ahn, C. H. , Dilmore, R. , and Wang, J. Y. , 2014, “ Development of Innovative and Efficient Hydraulic Fracturing Numerical Simulation Model and Parametric Studies in Unconventional Naturally Fractured Reservoirs,” J. Unconv. Oil Gas Resour., 8, pp. 25–45. [CrossRef]
Hofmann, H. , Babadagli, T. , and Zimmermann, G. , 2014, “ Numerical Simulation of Complex Fracture Network Development by Hydraulic Fracturing in Naturally Fractured Ultratight Formations,” ASME J. Energy Resour. Technol., 136(4), p. 042905. [CrossRef]
Keshavarzi, R. , and Mohammadi, S. , 2012, “ A New Approach for Numerical Modeling of Hydraulic Fracture Propagation in Naturally Fractured Reservoirs,” SPE/EAGE European Unconventional Resources Conference & Exhibition—From Potential to Production, Vienna, Austria, Mar. 20–22, Paper No. SPE-152509-MS.
Meyer, R. B. , and Bazan, L. W. , 2011, “ A Discrete Fracture Network Model for Hydraulically Induced Fractures—Theory, Parametric and Case Studies,” SPE Hydraulic Fracturing Technology Conference, The Woodlands, TX, Jan. 24–26, Paper No. SPE 140514.
Salehi, S. , and Nygaard, R. , 2014, “ Full Fluid–Solid Cohesive Finite-Element Model to Simulate Near Wellbore Fractures,” ASME J. Energy Resour. Technol., 137(1), p. 012903. [CrossRef]
Taleghani, A. D. , and Olson, J. E. , 2011, “ Numerical Modeling of Multistranded-Hydraulic-Fracture Propagation: Accounting for the Interaction Between Induced and Natural Fractures,” SPE J., 16(3), pp. 575–581. [CrossRef]
Taleghani, A. D. , and Olson, J. E. , 2013, “ How Natural Fractures Could Affect Hydraulic-Fracture Geometry,” SPE J., 19(1), pp. 161–171. [CrossRef]
Weng, X. , Kresse, O. , Cohen, C. , Wu, R. , and Gu, H. , 2011, “ Modeling of Hydraulic Fracture Network Propagation in a Naturally Fractured Formation,” SPE Prod. Oper., 26(4), pp. 368–380. [CrossRef]
Jahromi, M. Z. , Wang, J. Y. , and Ertekin, T. , 2013, “ Development of a Three-Dimensional Three-Phase Fully Coupled Numerical Simulator for Modeling Hydraulic Fracture Propagation in Tight Gas Reservoirs,” SPE Hydraulic Fracturing Technology Conference (SPE HFTC), Woodlands, TX, Feb. 4–6, Paper No. SPE-163850-MS.
Zhang, X. , Jeffery, R. G. , and Thiercelin, M. , 2007, “ Deflection and Propagation of Fruid-Driven Fractures at Frictional Bedding Interfaces: A Numerical Investigation,” J. Struct. Geol., 29(3), pp. 396–410. [CrossRef]
Zhang, X. , Jeffery, R. G. , and Thiercelin, M. , 2007, “ Effects of Frictional Geological Discontinuities on Hydraulic Fracture Propagation,” SPE Hydraulic Fracture Technology Conference (SPE HFTC), College Station, TX, Jan. 29–31, Paper No. SPE 106111.
Stoner, M. A. , 1969, “ Steady-State Analysis of Gas Production, Transmission and Distribution Systems,” Fall Meeting of the Society of Petroleum Engineers of AIME, Denver, CO, Sept. 28–Oct. 1, Paper No. SPE 2554.
Zhang, X. , and Jeffery, R. G. , 2008, “ Reinitiation or Termination of Fluid Driven Fractures at Frictional Bedding Interfaces,” J. Geophys. Res., 113(B8), p. B08416. [CrossRef]
Zhou, D. , Zheng, P. , Peng, J. , and He, P. , 2015, “ Induced Stress and Interaction of Fractures During Hydraulic Fracturing in Shale Formation,” ASME J. Energy Resour. Technol., 137(6), p. 062902. [CrossRef]
Gu, H. , Weng, X. , Lund, J. B. , Mack, M. G. , Ganguly, U. , and Suarez-Rivera, R. , 2011, “ Hydraulic Fracture Crossing Natural Fracture at Nonorthogonal Angles: A Criterion, Its Validation and Applications,” SPE Hydraulic Fracture Technology Conference (SPE HFTC), The Woodlands, TX, Jan. 24–26, Paper No. SPE-139984-MS.
Poltluri, N. , Zhu, D. , and Hill, A. D. , 2005, “ The Effect of Natural Fractures on Hydraulic Fracture Propagation,” SPE European Formation Damage Conference, Sheveningen, The Netherlands, May 25–27, Paper No. SPE 94568.
Howard, G. C. , and Fast, C. R. , 1970, Hydraulic Fracturing (SPE Monograph Series, Vol. 2), Society of Petroleum Engineers, Richardson, TX.
Perkins, T. K. , and Kern, L. R. , 1961, “ Widths of Hydraulic Fractures,” J. Pet. Technol., 13(9), pp. 934–949. [CrossRef]
Geertsma, J. , and de Klerk, F. , 1969, “ A Rapid Method of Predicting Width and Extent of Hydraulically Induced Fractures,” J. Pet. Technol., 21(12), pp. 1571–1581. [CrossRef]
Carter, R. D. , 1957, “ Derivation of the General Equation for Estimation the Extent of the Fracture Area,” Drilling and Production Practice, G. C. Howard and C. R. Fast, eds., American Petroleum Institute, New York, pp. 261–268.
Settari, A. , 1985, “ A New General Model of Fluid Loss in Hydraulic Fracturing,” SPE J., 25(4), pp. 491–501. [CrossRef]
Abousleiman, Y. N. , 1991, “ A Poroelastic PKN Model With Pressure Dependent Leakoff and Formation Permeability Determination,” Ph.D. dissertation, University of Delaware, Newark, DE.
Baree, R. D. , and Mukherjee, H. , 1996, “ Determination of Pressure Dependent Leakoff and Its Effect on Fracture,” SPE Annual Technical Conference and Exhibition, Denver, CO, Oct. 6–9, Paper No. SPE-36424-MS.
Carslaw, H. , and Jaeger, J. C. , 1956, Conduction of Heat in Solids, 2nd ed., Oxford University Press, Oxford, UK.
Fan, Y. , and Economides, M. J. , 1995, “ Fracture Dimensions in Frac & Pack Stimulation,” SPE J., 1(4), p. SPE-30469-PA.
Warren, J. E. , and Root, P. J. , 1963, “ The Behavior of Naturally Fractured Reservoirs,” SPE J., 3(3), pp. 245–255. [CrossRef]
Kazemi, H. , Merrill, L. S. , Porterfield, K. L. , and Zeman, P. R. , 1976, “ Numerical Simulation of Water-Oil Flow in Naturally Fracture Reservoirs,” SPE J., 16(6), pp. 317–326. [CrossRef]
Thomas, L. K. , Dixon, T. N. , and Pierson, R. G. , 1983, “ Fractured Reservoir Simulation,” SPE J., 23(1), pp. 42–54. [CrossRef]
Zhang, X. , Du, C. , Deimbacher, F. , Crick, M. , and Harikesavanallur, A. , 2009, “ Sensitivity Studies of Horizontal Wells With Hydraulic Fractures in Shale Gas Reservoirs,” International Petroleum Technology Conference, Doha, Qatar, Dec. 7–9, Paper No. IPTC-13338-MS.
Saidi, A. M. , 1983, “ Simulation of Naturally Fractured Reservoirs,” SPE Reservoir Simulation Symposium, San Francisco, CA, Nov. 15–18, Paper No. SPE-12270-MS.
Gilman, J. R. , 1986, “ An Efficient Finite-Difference Method for Simulating Phase Segregation in the Matrix Blocks in Double-Porosity Reservoirs,” SPE Reservoir Eng., 1(4), p. SPE-12271-PA.
Wu, Y.-S. , and Pruess, K. , 1988, “ A Multiple-Porosity Method for Simulation of Naturally Fractured Petroleum Reservoirs,” SPE Reservoir Eng., 3(1), pp. 327–336. [CrossRef]
Beckner, B. L. , Chan, H. M. , McDonald, A. E. , Wooten, S. O. , and Jones, T. A. , 1991, “ Simulating Naturally Fractured Reservoirs Using a Subdomain Method,” SPE Symposium on Reservoir Simulation, Anaheim, CA, Feb. 17–20, Paper No. SPE-21241-MS.
Lamb, H. , ed., 1932, Hydrodynamics, 6th ed., Dover Publications, New York.
England, A. H. , and Green, A. E. , 1963, “ Some Two-Dimensional Punch and Crack Problems in Classical Elasticity,” Math. Proc. Cambridge Philos. Soc., 59(2), pp. 489–500. [CrossRef]
Warpinski, N. R. , and Smith, M. B. , 1989, “ Rock Mechanics and Fracture Geometry,” Recent Advances in Hydraulic Fracturing (Monograph Series), J. L. Gidley , ed., SPE, Richardson, TX.
Sneddon, I. N. , 1946, “ The Distribution of Stress in the Neighbourhood of a Crack in an Elastic Solid,” Proc. R. Soc. London A, 187(1009), pp. 2229–2260. [CrossRef]
Olson, J. E. , 2004, “ Predicting Fracture Swarms—The Influence of Subcritical Crack Growth and the Crack-Tip Process Zone on Joint Spacing in Rock,” The Initiation, Propagation, and Arrest of Joints and Other Fractures Geological Society Special Publication, Vol. 231, T. Engelder , and J. W. Cosgrove , eds., Geological Society of London, London, pp. 73–87.
Pollard, D. D. , and Segall, P. , 1987, “ Theoretical Displacements and Stresses Near Fractures in Rock: With Applications to Faults, Joints, Veins, Dikes, and Solution Surfaces,” Fracture Mechanics of Rock (Academy Press Geology Series), B. K. Atkinson , ed., Elsevier, London, pp. 277–349.
Hossain, M. M. , and Rahman, M. K. , 2008, “ Numerical Simulation of Complex Fracture Growth During Tight Reservoir Stimulation by Hydraulic Fracturing,” J. Pet. Sci. Eng., 60(2), pp. 86–104. [CrossRef]
Gidley, J. L. , Holditch, S. A. , Nierode, D. E. , and Veatch, R. W., Jr. , 1989, Recent Advances in Hydraulic Fracturing (SPE Monograph Series, Vol. 12), Society of Petroleum Engineers, Richardson, TX.
Valko, P. , and Economides, M. J. , 1997, Hydraulic Fracture Mechanics, Wiley, New York.
Jacot, R. H. , Bazan, L. W. , and Meyer, B. R. , 2010, “ Technology Integration—A Methodology to Enhance Production and Maximize Economics in Horizontal Marcellus Shale Wells,” SPE Annual Technical Conference and Exhibition, Florence, Italy, Sept. 19–22, Paper No. SPE 135262-MS.


Grahic Jump Location
Fig. 1

Visual representation of height growth model

Grahic Jump Location
Fig. 2

(a) Length (ft) versus time (min) and (b) width (in.) versus time (min). The grid sizes of 0.1, 1, 10, and 100 ft are studied.

Grahic Jump Location
Fig. 3

Schematic of gridding system in both fracture and matrix domain

Grahic Jump Location
Fig. 4

Fracture domain discretization along x-direction

Grahic Jump Location
Fig. 5

Schematic flowchart of the three-dimensional, two-phase numerical solution

Grahic Jump Location
Fig. 6

Comparison results showing favorable comparison between PKN and model developed herein for fracture height evolution through time (a), and fracture width versus time (b)

Grahic Jump Location
Fig. 7

Visual representation of three-layer height growth model input values

Grahic Jump Location
Fig. 8

Height growth comparison result for seven-layer case. Solid lines indicate the result from the model developed herein while square symbols represent the results from the reference model.

Grahic Jump Location
Fig. 9

Vertical stress profile for each layer

Grahic Jump Location
Fig. 10

Height comparison between the model developed herein and DFN. (Left: height versus network length from the model developed herein and the black and grey to white color bar indicates the width of a fracture in ft; right: height versus width of the fracture at the perforation (primary fracture) from DFN. The contour map shows the width of a primary fracture as length increases.)

Grahic Jump Location
Fig. 11

Final fracture network comparison (top view). Top network is from Jacot's result and the bottom is from the model developed herein. The horizontal well is located along the zero in x-axis and perforation is placed at (0,0).

Grahic Jump Location
Fig. 12

Final fracture network characteristics (height, length, and aperture) from the model developed herein. This representation is for only half of a fracture network.



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In