A 3D peridynamic simulation of hydraulic fracture process in a heterogeneous medium

Abstract A non-local state-based peridynamic method (PD) is employed to simulate the initiation and propagation of the hydraulic fracturing. This study presents a novel 3D simulation of the hydraulic fracturing initiation and propagation in a heterogeneous geological medium by a linear-viscoelastic PD model. The interaction between induced fractures and a pre-existing fracture, as well as the influence of the angle of approach and differential horizontal stress on the fracture propagation behavior was investigated through the simulation. In the lower approach angles and differential horizontal stress, the pre-exiting fracture was dilated by the hydraulic fracture; in higher approach angles with higher differential horizontal stress, the hydraulic fracture crossed the pre-existing fracture. The accuracy of the 3D model was validated by experimental hydraulic fracturing data published based on the laboratory tests. Finally, the influences of the three different injection rates on fracture initiation and propagation, induced fracture surface and shear dilation were investigated. The simulation showed that the maximum dilatation and fractured zone occurred at the injection rate of 0.61 m 3 /min. The 0.61 m 3 /min injection rate caused the highest complete damage (0.9–1) with 5.24% of the total number of atoms. As a result, the peridynamic approach presents promising results in predicting fracture propagation and damage area.

[1]  Philippe H. Geubelle,et al.  Non-ordinary state-based peridynamic analysis of stationary crack problems , 2014 .

[2]  Robert Lipton,et al.  Multiscale Dynamics of Heterogeneous Media in the Peridynamic Formulation , 2012 .

[3]  Mofazzal Hossain,et al.  A shear‐dilation‐based model for evaluation of hydraulically stimulated naturally fractured reservoirs , 2002 .

[4]  Ove Stephansson,et al.  Fundamentals of discrete element methods for rock engineering , 2007 .

[5]  S. Silling,et al.  Peridynamic modeling of membranes and fibers , 2004 .

[6]  Correlative XRM and FIB-SEM for Non/Organic Pore Network Characterization in Woodford Shale Rock Matrix , 2015 .

[7]  R. Lehoucq,et al.  Peridynamic Theory of Solid Mechanics , 2010 .

[8]  M. Ortiz,et al.  FINITE-DEFORMATION IRREVERSIBLE COHESIVE ELEMENTS FOR THREE-DIMENSIONAL CRACK-PROPAGATION ANALYSIS , 1999 .

[9]  Kun Zhou,et al.  Analysis and Approximation of Nonlocal Diffusion Problems with Volume Constraints , 2012, SIAM Rev..

[10]  Failure prediction in fully coupled thermal and deformational fields with peridynamics , 2012, 2012 IEEE 62nd Electronic Components and Technology Conference.

[11]  Shanyong Wang,et al.  2D-numerical analysis of hydraulic fracturing in heterogeneous geo-materials , 2009 .

[12]  T. L. Blanton.  Propagation of Hydraulically and Dynamically Induced Fractures in Naturally Fractured Reservoirs , 1986 .

[13]  R. Lehoucq,et al.  Peridynamics for multiscale materials modeling , 2008 .

[14]  Qiang Du,et al.  A New Approach for a Nonlocal, Nonlinear Conservation Law , 2012, SIAM J. Appl. Math..

[15]  Jianqiao Ye,et al.  Dynamic simulation of crack initiation and propagation in cross-ply laminates by DEM , 2011 .

[16]  Yan Jin,et al.  Analysis of fracture propagation behavior and fracture geometry using a tri-axial fracturing system in naturally fractured reservoirs , 2008 .

[17]  Bahattin Kilic,et al.  Peridynamic theory for progressive failure prediction in homogeneous and heterogeneous materials , 2008 .

[18]  N. Sakhavand.  Parallel simulation of reinforced concrete sructures using peridynamics , 2011 .

[19]  Florin Bobaru,et al.  A peridynamic formulation for transient heat conduction in bodies with evolving discontinuities , 2012, J. Comput. Phys..

[20]  Stewart Andrew Silling,et al.  Dynamic fracture modeling with a meshfree peridynamic code , 2003 .

[21]  S. Silling,et al.  A meshfree method based on the peridynamic model of solid mechanics , 2005 .

[22]  Daniel Z. Turner,et al.  A Non-Local Model for Fluid-Structure Interaction with Applications in Hydraulic Fracturing , 2012, ArXiv.

[23]  F. Donze,et al.  Discrete element modelling of concrete submitted to dynamic loading at high strain rates , 2004 .

[24]  S. Silling,et al.  Peridynamic States and Constitutive Modeling , 2007 .

[25]  Bhanu Kiran Tuniki.  Peridynamic constitutive model for concrete , 2012 .

[26]  Hossein Honarvar Gheitanbaf,et al.  Parallel simulation of particle dynamics with application to micropolar peridynamic lattice modeling of reinforced concrete Structures , 2013 .

[27]  Mostafa Sharifzadeh,et al.  Triaxial creep tests and back analysis of time-dependent behavior of Siah Bisheh cavern by 3-Dimensional Distinct Element Method , 2011 .

[28]  S. P. Hunt,et al.  Modelling the Kaiser effect and deformation rate analysis in sandstone using the discrete element method , 2003 .

[29]  Mukul M. Sharma,et al.  A peridynamic formulation of pressure driven convective fluid transport in porous media , 2014, J. Comput. Phys..

[30]  P. Cundall,et al.  A discrete numerical model for granular assemblies , 1979 .

[31]  Colin Thornton,et al.  Numerical studies of uniaxial powder compaction process by 3D DEM , 2004 .

[32]  K. Shahriar,et al.  Experimental creep tests and prediction of long-term creep behavior of grouting material , 2014, Arabian Journal of Geosciences.

[33]  S. Silling.  Reformulation of Elasticity Theory for Discontinuities and Long-Range Forces , 2000 .