Application of Lightning Breakdown Simulation in Inversion of Induced Fracture Network Morphology in Stimulated Reservoirs

2021 ◽  
Author(s):  
Zhao Hui ◽  
Sheng Guanglong ◽  
Huang Luoyi ◽  
Zhong Xun ◽  
Fu Jingang ◽  
...  

Abstract Accurately characterizing fracture network morphology is necessary for flow simulation and fracturing evaluation. The complex natural fractures and reservoir heterogeneity in unconventional reservoirs make the induced fracture network resulting from hydraulic fracturing more difficult to describe. Existing fracture propagation simulation and fracture network inversion methods cannot accurately match actual fracture network morphology. Considering the lightning breakdown similar as fracture propagation, a new efficient approach for inversion of fracture network morphology is proposed. Based on the dielectric breakdown model (DBM) for lightning breakdown simulation and similarity principle, an induced fracture propagation algorithm integrating reservoir in-situ stress, rock mechanical parameters, and stress shadow effect is proposed. The fractal index and random function are coupled to quantitatively characterize the probability distribution of induced fracture propagation path. At the same time, a matching rate function is proposed to quantitatively evaluate the fitting between fracture network morphology and the micro seismic data. Combined with automatic history matching method, the actual fracture network morphology can be inverted with the matching rate as objective function. The proposed approach is applied to fracture network simulation of mult-fractured horizontal wells of shale oil reservoir in China, and the fracture networks from inversion fit well with the micro seismic data. A simulation of 94 fractures in the 32 section of Well X2 shows that the well propagates more obvious branch fractures. The single-wing fracture network communicates approximately 200m horizontally and approximately 10m vertically. In single fracture flow simulation, it is necessary to consider the influence of complex fracture network morphology, but when simulating fluid flow for a single well or even a reservoir, only the main fracture needs to be considered. This paper proposes an induced fracture propagation algorithm that integrates reservoir in-situ stress, rock mechanical parameters, and stress shadowing effects. This algorithm greatly improves the calculation efficiency on the premise of ensuring the accuracy of induced fracture network morphology. The approach in this paper provides a theoretical basis for flow simulation of stimulated reservoirs and optimization of fracture networks.

Energies ◽  
2020 ◽  
Vol 13 (18) ◽  
pp. 4718
Author(s):  
Song Wang ◽  
Jian Zhou ◽  
Luqing Zhang ◽  
Zhenhua Han

Hydraulic fracturing is a key technical means for stimulating tight and low permeability reservoirs to improve the production, which is widely employed in the development of unconventional energy resources, including shale gas, shale oil, gas hydrate, and dry hot rock. Although significant progress has been made in the simulation of fracturing a single well using two-dimensional Particle Flow Code (PFC2D), the understanding of the multi-well hydraulic fracturing characteristics is still limited. Exploring the mechanisms of fluid-driven fracture initiation, propagation and interaction under multi-well fracturing conditions is of great theoretical significance for creating complex fracture networks in the reservoir. In this study, a series of two-well fracturing simulations by a modified fluid-mechanical coupling algorithm were conducted to systematically investigate the effects of injection sequence and well spacing on breakdown pressure, fracture propagation and stress shadow. The results show that both injection sequence and well spacing make little difference on breakdown pressure but have huge impacts on fracture propagation pressure. Especially under hydrostatic pressure conditions, simultaneous injection and small well spacing increase the pore pressure between two injection wells and reduce the effective stress of rock to achieve lower fracture propagation pressure. The injection sequence can change the propagation direction of hydraulic fractures. When the in-situ stress is hydrostatic pressure, simultaneous injection compels the fractures to deflect and tend to propagate horizontally, which promotes the formation of complex fracture networks between two injection wells. When the maximum in-situ stress is in the horizontal direction, asynchronous injection is more conducive to the parallel propagation of multiple hydraulic fractures. Nevertheless, excessively small or large well spacing reduces the number of fracture branches in fracture networks. In addition, the stress shadow effect is found to be sensitive to both injection sequence and well spacing.


Geofluids ◽  
2019 ◽  
Vol 2019 ◽  
pp. 1-14 ◽  
Author(s):  
Chuanyin Jiang ◽  
Xiaoguang Wang ◽  
Zhixue Sun ◽  
Qinghua Lei

We investigated the effect of in situ stresses on fluid flow in a natural fracture network. The fracture network model is based on an actual critically connected (i.e., close to the percolation threshold) fracture pattern mapped from a field outcrop. We derive stress-dependent fracture aperture fields using a hybrid finite-discrete element method. We analyze the changes of aperture distribution and fluid flow field with variations of in situ stress orientation and magnitude. Our simulations show that an isotropic stress loading tends to reduce fracture apertures and suppress fluid flow, resulting in a decrease of equivalent permeability of the fractured rock. Anisotropic stresses may cause a significant amount of sliding of fracture walls accompanied with shear-induced dilation along some preferentially oriented fractures, resulting in enhanced flow heterogeneity and channelization. When the differential stress is further elevated, fracture propagation becomes prevailing and creates some new flow paths via linking preexisting natural fractures, which attempts to increase the bulk permeability but attenuates the flow channelization. Comparing to the shear-induced dilation effect, it appears that the propagation of new cracks leads to a more prominent permeability enhancement for the natural fracture system. The results have particularly important implications for predicting the hydraulic responses of fractured rocks to in situ stress fields and may provide useful guidance for the strategy design of geofluid production from naturally fractured reservoirs.


1982 ◽  
Vol 22 (03) ◽  
pp. 341-349 ◽  
Author(s):  
H.A.M. van Eekelen

Abstract One of the main problems in hydraulic fracturing technology is the prediction of fracture height. In particular, the question of what constitutes a barrier to vertical fracture propagation is crucial to the success of field operations. An analysis of hydraulic fracture containment effects has been performed. The main conclusion is that in most cases the fracture will penetrate into the layers adjoining the pay zone, the depth of penetration being determined by the differences in stiffness and in horizontal in-situ stress between the pay zone and the adjoining layers. For the case of a stiffness contrast, an estimate of the penetration depth is given. Introduction Current design procedures for hydraulic fracturing of oil and gas reservoirs are based predominantly on the fracturing theories of Perkins and Kern, Nordgren, and Geertsma and de Klerk. In the model proposed by Perkins and Kern, and improved by Nordgren, the formation stiffness is concentrated in vertical planes perpendicular to the direction of fracture propagation, The fracture cross section in these planes is assumed elliptical, and the stiffness of the formation in the horizontal plane is neglected. In the model proposed by Geertsma and de Klerk, the stiffness of the formation is concentrated in the horizontal plane. The fracture cross section in the vertical plane is assumed rectangular, and the stiffness in the vertical plane is neglected. In both models, the fluid pressure is assumed a function of the distance from the borehole, independent of the transverse coordinates. The theory by Perkins and Kern is more appropriate for long fractures (L/H >1, where L and H are length and height of the fracture), whereas the model by Geertsma and de Klerk is applicable for short fractures, L/H less than 1. The main shortcoming of these fracture-design procedures is that they assume a constant, preassigned fracture height. H. The value of H has a strong influence on the result, for fracture length, fracture width, and proppant transport. Usually, the estimated fracture height is based on assumed "barrier action" of rock layers above and below the pay zone. This situation is rather unsatisfactory. Moreover, if these layers do not contain the fracture, large volumes of fracturing fluid may be lost in fracturing unproductive strata, and communication with unwanted formations may be opened up. Whether an adjacent formation will act as a fracture barrier may depend on a number of factors: differences in in-situ stress, elastic properties, fracture toughness, ductility, and permeability; and the bonding at the interface. We analyze these factors with respect to their relative influence on fracture containment. Differences in in-situ stress and differences in elastic properties affect the global or overall stress field around the fracture, and, hence, the three-dimensional shape of the fracture. This shape, together with the horizontal and vertical fracture propagation rates, determines the fluid pressure distribution in the fracture, which in turn affects the stress field around the fracture. Consequently, the elastic stress field, the fluid pressure field, and the fracture propagation pattern are intimately coupled, which makes the fracture propagation problem a complicated one. Whether at a certain point of the fracture edge the fracture will propagate is determined by the intensity of the stress concentration at that point. This stress concentration depends on the global stress distribution in and around the fracture, but it also is affected directly by local ductility, permeability, and elastic modulus in the tip region. SPEJ P. 341^


2020 ◽  
Vol 54 ◽  
pp. 149-156
Author(s):  
Ajay K. Sahu ◽  
Ankur Roy

Abstract. It is well known that fracture networks display self-similarity in many cases and the connectivity and flow behavior of such networks are influenced by their respective fractal dimensions. In the past, the concept of lacunarity, a parameter that quantifies spatial clustering, has been implemented by one of the authors in order to demonstrate that a set of seven nested natural fracture maps belonging to a single fractal system, but of different visual appearances, have different clustering attributes. Any scale-dependency in the clustering of fractures will also likely have significant implications for flow processes that depend on fracture connectivity. It is therefore important to address the question as to whether the fractal dimension alone serves as a reasonable proxy for the connectivity of a fractal-fracture network and hence, its flow response or, if it is the lacunarity, a measure of scale-dependent clustering, that may be used instead. The present study attempts to address this issue by exploring possible relationships between the fractal dimension, lacunarity and connectivity of fractal-fracture networks. It also endeavors to study the relationship between lacunarity and fluid flow in such fractal-fracture networks. A set of deterministic fractal-fracture models generated at different iterations and, that have the same theoretical fractal dimension are used for this purpose. The results indicate that such deterministic synthetic fractal-fracture networks with the same theoretical fractal dimension have differences in their connectivity and that the latter is fairly correlated with lacunarity. Additionally, the flow simulation results imply that lacunarity influences flow patterns in fracture networks. Therefore, it may be concluded that at least in synthetic fractal-fracture networks, rather than fractal dimension, it is the lacunarity or scale-dependent clustering attribute that controls the connectivity and hence the flow behavior.


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.


Geofluids ◽  
2019 ◽  
Vol 2019 ◽  
pp. 1-15 ◽  
Author(s):  
Yongxiang Zheng ◽  
Jianjun Liu ◽  
Yun Lei

The formation of the fracture network in shale hydraulic fracturing is the key to the successful development of shale gas. In order to analyze the mechanism of hydraulic fracturing fracture propagation in cemented fractured formations, a numerical simulation about fracture behavior in cemented joints was conducted based firstly on the block discrete element. And the critical pressure of three fracture propagation modes under the intersection of hydraulic fracturing fracture and closed natural fracture is derived, and the parameter analysis is carried out by univariate analysis and the response surface method (RSM). The results show that at a low intersecting angle, hydraulic fractures will turn and move forward at the same time, forming intersecting fractures. At medium angles, the cracks only turn. At high angles, the crack will expand directly forward without turning. In conclusion, low-angle intersecting fractures are more likely to form complex fracture networks, followed by medium-angle intersecting fractures, and high-angle intersecting fractures have more difficulty in forming fracture networks. The research results have important theoretical guiding significance for the hydraulic fracturing design.


2006 ◽  
Vol 324-325 ◽  
pp. 383-386 ◽  
Author(s):  
Zhi Long Lian ◽  
Xiu Xi Wang ◽  
Heng An Wu ◽  
Bing Xue ◽  
J. Zhang ◽  
...  

Numerical simulation of hydraulic fracturing propagations in the permeable reservoirs was carried out with the finite element analysis software (ABAQUS). A model of coupling the stress equilibrium and fluid continuity equations was proposed and implemented. The nonuniform of sink pore pressure on the fracture surfaces which changes associated with the propagation of fracture was described by a self-developed subroutine through the FLOW in ABAQUS. Samples under different conditions were conducted for studying the rules of the propagation of hydraulic fracturing. The results show that the permeability at the fracture tip is more serious than any other places of the fracture face. The model also illustrates that the fracture geometry is mainly determined by the minimal in-situ stress. The model can be used to simulate the effects of hydraulic fracturing pressures and injection rates on fracture propagation. The results are of much significance for the design of hydraulic fracturing treatments.


Sign in / Sign up

Export Citation Format

Share Document