Imaging Seismic Deformation Induced by Hydraulic Fracture Complexity

2009 ◽  
Vol 12 (01) ◽  
pp. 48-52 ◽  
Author(s):  
Shawn C. Maxwell ◽  
Charles Waltman ◽  
Norman R. Warpinski ◽  
Michael J. Mayerhofer ◽  
Neda Boroumand

Summary Microseismic mapping is extensively used in the Barnett Shale to map hydraulic fracture complexity associated with interactions of the stimulation with pre-existing fractures (fracs). Previous studies have indicated a fair correlation between the well performance and extent of the seismically active volume. However, in addition to this measure of the extent of the stimulated fracture network, the characteristics of this fracture network is also expected to impact the well performance. In particular, the fracture spacing is believed to be an important factor controlling the potential gas flow. In this paper, we use the density of the total seismic moment release (a robust measure of the microseism strength) as an indication of the seismic deformation that may correlate to the fracture density. The study uses a set of microseismic maps of hydraulic fracture stimulations, including cases in which the stimulated reservoir volume measured by the extent of the seismically active region poorly correlated with the well performance. Incorporating the seismic moment density to assess the fracture density with the network extent, an improved correlation with the well performance was observed. Introduction Microseismic mapping of hydraulic fracture stimulations has become a common technique to map the fracture growth and geometry (Warpinski et al. 2004; Fisher et al. 2002; Maxwell et al. 2002; Fisher et al. 2004; Rutledge et al. 2004; Shapiro et al. 2004; Chambers et al. 2008; Lu et al. 2008; Warpinski et al. 2005). Microseismic images provide details of the fracture azimuth, height, length, and complexity resulting from interaction with pre-existing fratures. The resulting images can be used to calibrate numerical simulations of the fracture growth, allowing more confident modeling of other stimulations in the field, and a better identification of the stimulated region that may ultimately be drained by the well. Arguably, the Barnett Shale is the field that has had the most fracs mapped over the last several years. Microseismic mapping in the Barnett Shale has repeatedly demonstrated extreme fracture complexity resulting from interaction between the injection and a pre-existing fracture network (Fisher et al. 2002; Maxwell et al. 2002; Fisher et al. 2004; Rutledge et al. 2004; Shapiro et al. 2004; Chambers et al. 2008; Lu et al. 2008; Warpinski et al. 2005; Mayerhofer et al. 2006). Even between neighboring wells, the geometry of the stimulated fracture network shows a high degree of variability caused by localized differences in the fracture network (Fisher et al. 2002). The Barnett Shale has a low-intrinsic matrix permeability, and the permeability enhancement associated with the fracture stimulation results in permeable fracture networks sufficient for economic gas recovery in the field. Previous studies have shown a correlation between the volume of the reservoir stimulated as measured by the volume of the reservoir that emits microseisms during the stimulation, and the production ultimately realized from the well (Fisher et al. 2002; Fisher et al. 2004; Mayerhofer et al. 2006). The correlation is attributed to larger fracture networks being stimulated in wells in which a large microseismically active volume of the reservoir has been realized, resulting in more permeable fracture pathways connected to the well and therefore a higher potential for gas flow to the well. Recently, many operators in the Barnett Shale have attempted horizontal completions, which have allowed large volumes of the reservoir to be stimulated with large fracture networks. Many of these completions use perforated, cemented liners, and the microseismic images allow for indentification of improved perforation staging to maximize the stimulated reservoir volume (SRV) (Fisher et al. 2004). Many of the Barnett Shale stimulations are water fracs in which large volumes of water are injected at a high rate (Mayerhofer et al. 1997). One possible mechanism for the success of waterfracs is that increased fluid pressure in natural fractures induced shear failure, resulting in fracture dilation associated with mismatched surfaces on opposite sides of the fracture. Within this conceptual framework, the microseismic events correspond to the actual fracture movement. The earlier investigations of the SRV measured the total volume of the microseismically active region. However, this measure of the stimulated volume does not take into account the properties of the fracture network, which has also been indicated to impact well performance (Mayerhofer et al. 2006). Furthermore, the permeability enhancement of the fracture may be related to deformation associated with fracturing. Beyond the basic hypocentral locations of the microseisms used to calculate the SRV, additional seismic signal characteristics allow investigation of the source of the mechanical deformation resulting in the microseisms. In particular, the seismic moment (Aki and Richards 1980), a robust measure of the strength of an earthquake or microearthquake can be used to quantify the seismic deformation (Maxwell et al. 2003). In this paper, we examine several published microseismic projects in the Barnett Shale formation for correlation between the production and seismic-deformation attributes. In the next section, we describe seismic moments and the calculation of seismic deformation. We illustrate how a seismic moment can be used to remove a recording bias present in most microseismic monitoring applications and the importance for calculating the seismic deformation. Finally, we present the comparison between production, seismic deformation, and SRV for several published datasets.

Author(s):  
Yunsuk Hwang ◽  
Jiajing Lin ◽  
David Schechter ◽  
Ding Zhu

Multiple hydraulic fracture treatments in reservoirs with natural fractures create complex fracture networks. Predicting well performance in such a complex fracture network system is an extreme challenge. The statistical nature of natural fracture networks changes the flow characteristics from that of a single linear fracture. Simply using single linear fracture models for individual fractures, and then summing the flow from each fracture as the total flow rate for the network could introduce significant error. In this paper we present a semi-analytical model by a source method to estimate well performance in a complex fracture network system. The method simulates complex fracture systems in a more reasonable approach. The natural fracture system we used is fractal discrete fracture network model. We then added multiple dominating hydraulic fractures to the natural fracture system. Each of the hydraulic fractures is connected to the horizontal wellbore, and some of the natural fractures are connected to the hydraulic fractures through the network description. Each fracture, natural or hydraulically induced, is treated as a series of slab sources. The analytical solution of superposed slab sources provides the base of the approach, and the overall flow from each fracture and the effect between the fractures are modeled by applying the superposition principle to all of the fractures. The fluid inside the natural fractures flows into the hydraulic fractures, and the fluid of the hydraulic fracture from both the reservoir and the natural fractures flows to the wellbore. This paper also shows that non-Darcy flow effects have an impact on the performance of fractured horizontal wells. In hydraulic fracture calculation, non-Darcy flow can be treated as the reduction of permeability in the fracture to a considerably smaller effective permeability. The reduction is about 2% to 20%, due to non-Darcy flow that can result in a low rate. The semi-analytical solution presented can be used to efficiently calculate the flow rate of multistage-fractured wells. Examples are used to illustrate the application of the model to evaluate well performance in reservoirs that contain complex fracture networks.


2021 ◽  
pp. 014459872098153
Author(s):  
Yanzhi Hu ◽  
Xiao Li ◽  
Zhaobin Zhang ◽  
Jianming He ◽  
Guanfang Li

Hydraulic fracturing is one of the most important technologies for shale gas production. Complex hydraulic fracture networks can be stimulated in shale reservoirs due to the existence of numerous natural fractures. The prediction of the complex fracture network remains a difficult and challenging problem. This paper presents a fully coupled hydromechanical model for complex hydraulic fracture network propagation based on the discontinuous deformation analysis (DDA) method. In the proposed model, the fracture propagation and rock mass deformation are simulated under the framework of DDA, and the fluid flow within fractures is simulated using lubrication theory. In particular, the natural fracture network is considered by using the discrete fracture network (DFN) model. The proposed model is widely verified against several analytical and experimental results. All the numerical results show good agreement. Then, this model is applied to field-scale modeling of hydraulic fracturing in naturally fractured shale reservoirs. The simulation results show that the proposed model can capture the evolution process of complex hydraulic fracture networks. This work offers a feasible numerical tool for investigating hydraulic fracturing processes, which may be useful for optimizing the fracturing design of shale gas reservoirs.


Energies ◽  
2020 ◽  
Vol 13 (3) ◽  
pp. 633
Author(s):  
Guangzhi Yang ◽  
Shicheng Zhang ◽  
Jia Wang ◽  
Ning Li ◽  
Xinfang Ma ◽  
...  

Exploring engineering methods for increasing fracture network complexity is important for the development of unconventional oil and gas reservoirs. In this study, we conducted a series of fracturing experiments on naturally fractured volcanic samples. An injection method, multiple flow pulses, is proposed to increase fracture complexity. The results show that fluid leaked into the natural fracture network (NFN) when the injection rate was low (0.2 mL/min); hydraulic-fracture-dominant fracture geometry was created with an injection rate of 2 and 5 mL/min. Under the 2 mL/min-injection scheme with 3 pulses, the injection pressure during the intermittent stage was low (<5 MPa), resulting in a limited increase in fracture complexity. When the number of the flow pulses increased to 5, the pressure drop rate in the fourth and fifth intermittent stage significantly increased, indicating an increase in the aperture of natural fractures (NFs) and in the fluid leak-off effect. Under the 5 mL/min injection scheme containing 5 pulses, besides the enhanced fluid leak-off, a sharp injection pressure drop was observed, indicating the activation of NFs. The complexity and the aperture of the ultimate fracture network further increased. The injection method, multiple flow pulses, can be used to create complex fracture networks effectively.


Materials ◽  
2019 ◽  
Vol 12 (15) ◽  
pp. 2387 ◽  
Author(s):  
Bo Zhang ◽  
Yong Li ◽  
Nicholas Fantuzzi ◽  
Yuan Zhao ◽  
Yan-Bao Liu ◽  
...  

Coal contains a large number of fractures, whose characteristics are difficult to describe in detail, while their spatial distribution patterns may follow some macroscopic statistical laws. In this paper, several fracture geometric parameters (FGPs) were used to describe a fracture, and the coal seam was represented by a two-dimensional stochastic fracture network (SFN) which was generated and processed through a series of methods in MATLAB. Then, the processed SFN image was able to be imported into COMSOL Multiphysics and converted to a computational domain through the image function. In this way, the influences of different FGPs and their distribution patterns on the permeability of the coal seam were studied, and a finite element model to investigate gas flow properties in the coal seam was carried out. The results show that the permeability of the coal seam increased with the rising of fracture density, length, aperture, and with the decrease of the angle between the fracture orientation and the gas pressure gradient. It has also been found that large-sized fractures have a more significant contribution to coal reservoir permeability. Additionally, a numerical simulation of CBM extraction was carried out to show the potential of the proposed approach in the application of tackling practical engineering problems. According to the results, not only the connectivity of fractures but also variations of gas pressure and velocity can be displayed explicitly, which is consistent well with the actual situation.


Fractals ◽  
2019 ◽  
Vol 27 (04) ◽  
pp. 1950057 ◽  
Author(s):  
TONGJUN MIAO ◽  
SUJUN CHENG ◽  
AIMIN CHEN ◽  
YAN XU ◽  
GUANG YANG ◽  
...  

Fractures with power law length distributions abound in nature such as carbonate oil and gas reservoirs, sandstone, hot dry rocks, etc. The fluid transport properties and morphology characterization of fracture networks have fascinated numerous researchers to investigate for several decades. In this work, the analytical models for fracture density and permeability are extended from fractal fracture network to general fracture network with power law length distributions. It is found that the fracture density is related to the power law exponents [Formula: see text] and the area porosity [Formula: see text] of fracture network. Then, a permeability model for the fracture length distribution with general power law exponent [Formula: see text] and the power law exponent [Formula: see text] for fracture length versus aperture is proposed based on the well-known cubic law in individual fracture. The analytical expression for permeability of fractured networks is found to be a function of power law exponents [Formula: see text], area porosity [Formula: see text] of fracture network, and the micro-structural parameters (maximum fracture length [Formula: see text], fracture azimuth [Formula: see text] and fracture dip angle [Formula: see text]). The present model may shed light on the mechanism of seepage in fracture networks with power law length distributions.


Processes ◽  
2018 ◽  
Vol 6 (8) ◽  
pp. 113 ◽  
Author(s):  
Shen Wang ◽  
Huamin Li ◽  
Dongyin Li

To investigate the mechanism of hydraulic fracture propagation in coal seams with discontinuous natural fractures, an innovative finite element meshing scheme for modeling hydraulic fracturing was proposed. Hydraulic fracture propagation and interaction with discontinuous natural fracture networks in coal seams were modeled based on the cohesive element method. The hydraulic fracture network characteristics, the growth process of the secondary hydraulic fractures, the pore pressure distribution and the variation of bottomhole pressure were analyzed. The improved cohesive element method, which considers the leak-off and seepage behaviors of fracturing liquid, is capable of modeling hydraulic fracturing in naturally fractured formations. The results indicate that under high stress difference conditions, the hydraulic fracture network is spindle-shaped, and shows a multi-level branch structure. The ratio of secondary fracture total length to main fracture total length was 2.11~3.62, suggesting that the secondary fractures are an important part of the hydraulic fracture network in coal seams. In deep coal seams, the break pressure of discontinuous natural fractures mainly depends on the in-situ stress field and the direction of natural fractures. The mechanism of hydraulic fracture propagation in deep coal seams is significantly different from that in hard and tight rock layers.


2018 ◽  
Vol 140 (3) ◽  
Author(s):  
Oliver Chang ◽  
Michael Kinzel ◽  
Robert Dilmore ◽  
John Yilin Wang

Horizontal drilling with successful multistage hydraulic fracture treatments is the most widely applied and effective method to enable economic development of hydrocarbon-bearing shale reservoirs. Once fracture networks are established, they must be propped open to maintain pathways for fluid migration through the production phase. As such, the design and application of effective and efficient proppant treatment is considered a key step to successfully develop the targeted resource. Unfortunately, the available literature and simulation tools to describe proppant transport in complex fracture networks are inadequate, and some of the fundamental mechanisms of proppant transport are poorly understood. The present study provides a critical review of relevant published literature to identify important mechanisms of particle transport and related governing equations. Based on that review, a mathematical model was developed to quantitatively predict the transport behavior of proppant particles in model fracture networks. Aspects of this mathematical model are compared against computational fluid dynamic (CFD) simulation, and implications of this work are discussed.


Lithosphere ◽  
2021 ◽  
Vol 2021 (Special 1) ◽  
Author(s):  
Jiazheng Liu ◽  
Xiaotong Liu ◽  
Hongzhang Zhu ◽  
Xiaofei Ma ◽  
Yuxue Zhang ◽  
...  

Abstract The gridless analytical and semianalytical methodologies can provide credible solutions for describing the well performance of the fracture networks in a homogeneous reservoir. Reservoir heterogeneity, however, is common in unconventional reservoirs, and the productivity can vary significantly along the horizontal wells drilled for producing such reservoirs. It is oversimplified to treat the entire reservoir matrix as homogeneous if there are regions with extremely nonuniform properties in the reservoir. However, the existing analytical and semianalytical methods can only model simple cases involving matrix heterogeneity, such as composite, layered, or compartmentalized reservoirs. A semianalytical methodology, which can model fracture networks in heterogeneous reservoirs, is still absent; in this study, we propose a decomposed fracture network model to fill this gap. We discretize a fractured reservoir into matrix blocks that are bounded by the fractures and/or the reservoir boundary and upscale the local properties to these blocks; therefore, a heterogeneous reservoir can be represented with these blocks that have nonuniform properties. To obtain a general flow equation to characterize the transient flow in the blocks that may exhibit different geometries, we approximate the contours of pressure with the contours of the depth of investigation (DOI) in each block. Additionally, the borders of each matrix block represent the fractures in the reservoir; thus, we can characterize the configurations of complex fracture networks by assembling all the borders of the matrix blocks. This proposed model is validated against a commercial software (Eclipse) on a multistage hydraulic fracture model and a fracture network model; both a homogeneous case and a heterogeneous case are examined in each of these two models. For the heterogeneous case, we assign different permeabilities to the matrix blocks in an attempt to characterize the reservoir heterogeneity. The calculation results demonstrate that our new model can accurately simulate the well performance even when there is a high degree of permeability heterogeneity in the reservoir. Besides, if there are high-permeability regions existing in the fractured reservoir, a BDF may be observed in the early production period, and formation linear flow may be indistinguishable in the early production period because of the influence of reservoir heterogeneity.


Solid Earth ◽  
2020 ◽  
Vol 11 (5) ◽  
pp. 1731-1746
Author(s):  
Romesh Palamakumbura ◽  
Maarten Krabbendam ◽  
Katie Whitbread ◽  
Christian Arnhardt

Abstract. Understanding the impact of fracture networks on rock mass properties is an essential part of a wide range of applications in geosciences from understanding permeability of groundwater aquifers and hydrocarbon reservoirs to erodibility properties and slope stability of rock masses for geotechnical engineering. However, gathering high-quality, oriented-fracture datasets in the field can be difficult and time-consuming, for example, due to constraints on field work time or access (e.g. cliffs). Therefore, a method for obtaining accurate, quantitative fracture data from photographs is a significant benefit. In this paper we describe a method for generating a series of digital fracture traces in a geographic information system (GIS) environment, in which spatial analysis of a fracture network can be carried out. The method is not meant to replace the gathering of data in the field but to be used in conjunction with it, and it is well suited when field work time is limited or when the section cannot be accessed directly. The basis of the method is the generation of the vector dataset (shapefile) of a fracture network from a georeferenced photograph of an outcrop in a GIS environment. From that shapefile, key parameters such as fracture density and orientation can be calculated. Furthermore, in the GIS environment more complex spatial calculations and graphical plots can be carried out such as heat maps of fracture density. Advantages and limitations compared to other fracture network capture methods are discussed.


Solid Earth ◽  
2020 ◽  
Vol 11 (6) ◽  
pp. 2119-2140
Author(s):  
Billy James Andrews ◽  
Zoe Kai Shipton ◽  
Richard Lord ◽  
Lucy McKay

Abstract. Fault architecture and fracture network evolution (and resulting bulk hydraulic properties) are highly dependent on the mechanical properties of the rocks at the time the structures developed. This paper investigates the role of mechanical layering and pre-existing structures on the evolution of strike–slip faults and fracture networks. Detailed mapping of exceptionally well exposed fluvial–deltaic lithologies at Spireslack Surface Coal Mine, Scotland, reveals two phases of faulting with an initial sinistral and later dextral sense of shear with ongoing pre-faulting, syn-faulting, and post-faulting joint sets. We find fault zone internal structure depends on whether the fault is self-juxtaposing or cuts multiple lithologies, the presence of shale layers that promote bed-rotation and fault-core lens formation, and the orientation of joints and coal cleats at the time of faulting. During ongoing deformation, cementation of fractures is concentrated where the fracture network is most connected. This leads to the counter-intuitive result that the highest-fracture-density part of the network often has the lowest open fracture connectivity. To evaluate the final bulk hydraulic properties of a deformed rock mass, it is crucial to appreciate the relative timing of deformation events, concurrent or subsequent cementation, and the interlinked effects on overall network connectivity.


Sign in / Sign up

Export Citation Format

Share Document