Efficient Modeling of Unconventional Well Performance with Millions of Natural and Hydraulic Fractures Using Embedded Discrete Fracture Model EDFM

The complexity of dynamic modeling for naturally fractured reservoirs has increased in recent years to incorporate more data and physics, as well as to handle advanced completion designs and development scenarios. While these complex models can provide more insight to difficult problems, they come with higher computational costs. Such a limitation prohibits an asset team from working with a large number of well/fracture scenarios that correctly represent geological uncertainty. This study presents a powerful non-intrusive Embedded Discrete Fracture Model (EDFM) method to efficiently handle millions of natural and hydraulic fractures with hundreds of horizontal wells, which has never been modeled in the literature. Specifically, we built a 3D geological model using a black oil reservoir simulator with 100 square miles in the horizontal area and 11 layers of 165 ft thickness. The total number of matrix cells without considering fractures is over 3 million. In total, 400 horizontal wells with well length of 6000 ft were modeled in two target layers. Each layer contains 200 wells. Each well has 112 hydraulic fractures with cluster spacing of 50 ft. The total number of hydraulic fractures is 44,800. In addition, we generated three cases with 10K, 100K and 1 million 3D natural fractures with dip angle from 70 to 90 degrees. For the case with 1 million natural fractures, the total number of cells is over 42 million. Well performance for the field example, with and without natural fractures, was investigated. This work adds significant value to the well and fracture spacing optimization process during field development planning. The non-intrusive EDFM method has been proven to be an efficient fracture modeling tool for simulating million-level complex hydraulic/natural fractures, which significantly improves accuracy and reduces computational time.

Discrete Fracture Model for Hydro-Mechanical Coupling in Fractured Reservoirs

The process of coupled flow and mechanics occurs in various environmental and energy applications, including conventional and unconventional fractured reservoirs. This work establishes a new formulation for modeling hydro-mechanical coupling in fractured reservoirs. The discrete-fracture model (DFM), in which the porous matrix and fractures are represented explicitly in the form of unstructured grid, has been widely used to describe fluid flow in fractured formations. In this work, we extend the DFM approach for modeling coupled flow-mechanics process, in which flow problems are solved using the multipoint flux approximation (MPFA) method, and mechanics problems are solved using the multipoint stress approximation (MPSA) method. The coupled flow-mechanics problems share the same computational grid to avoid projection issues and allow for convenient exchange between them. We model the fracture mechanical behavior as a two-surface contact problem. The resulting coupled system of nonlinear equations is solved in a fully-implicit manner. The accuracy and generality of the numerical implementation are accessed using cases with analytical solutions, which shows an excellent match. We then apply the methodology to more complex cases to demonstrate its general applicability. We also investigate the geomechanical influence on fracture permeability change using 2D rock fractures. This work introduces a novel formulation for modeling the coupled flow-mechanics process in fractured reservoirs, and can be readily implemented in reservoir characterization workflow.

Impact of Complex Fracture Networks on Rate Transient Behavior of Wells in Unconventional Reservoirs Based on Embedded Discrete Fracture Model

It has recently been demonstrated that complex fracture networks (CFN) especially activated natural fractures (ANF) play an important role in unconventional reservoir development. However, traditional rate transient analysis (RTA) methods barely investigate the impact of CFN or ANF. Furthermore, the influence of CFN on flow regime is still ambiguous. Failure to consider these effects could lead to misdiagnosis of flow regimes and underestimation of original oil in place (OOIP). A novel numerical RTA method is therefore presented herein to improve the quality of reserves assessment. A new methodology is introduced. Propagating hydraulic fractures (HF) can generate different stress perturbations to allow natural fractures (NF) to fail, forming various ANF pattern. An embedded discrete fracture model (EDFM) of ANF is stochastically generated instead of local grid refinement (LGR) method to overcome the time-intensive computation time. These models are coupled with reservoir models using non-neighboring connections (NNCs). Results show that except for simplified models used in previous studies subjected to traditional concept of stimulated reservoir volume (SRV), in our study, the ANF region has been discussed to emphasis the impact of NF on simulation results. Henceforth, ANF could be only concentrated around the near-wellbore region, and it may also cover the whole simulation area. Obvious distinctions could be viewed for different kinds of ANF on diagnostic plots. Instead of SRV-dominated flow mentioned in previous studies, ANF-dominated flow developed in this work is shown to be more reasonable. Also, new flow regimes such as interference flow inside and outside activated natural fracture flow region (ANFR) are found. In summary, better evaluation of reservoir properties and reserves assessment such as OOIP are achieved based on our proposed model compared with conventional models. The novel RTA method considering CFN presented herein is an easy-to-apply numerical RTA technique that can be applied for reservoir and fracture characterization as well as OOIP assessment.

Numerical Simulation of Gas-Water Two-Phase Flow in Deep Shale Gas Reservoir Development Based on Mixed Fracture Modeling

In this paper, the development of a three-dimensional, two-phase fluid flow model (Modified Embedded Discrete Fracture Model) to study flow performances of a fractured horizontal well in deep-marine shale gas is presented. Deep-marine shale gas resources account for nearly 80% in China, which is the decisive resource basis for large-scale shale gas production. The dynamic characteristics of deep shale gas reservoirs are quite different and more complex. This paper uses the embedded discrete fracture model to simulate artificial fractures (main fractures and secondary fractures) and the dual-media model to simulate the mixed fractured media of natural fractures and considers the flow characteristics of partitions (artificial fractures, natural fractures, and matrix). Gas desorption is considered in the matrix. Different degrees of stress sensitivity are considered for natural and artificial fractures. Aiming at accurately simulating the whole production history of horizontal well fracturing, especially the dynamic changes of postfracturing flowback, a postfracturing fluid initialization method based on fracturing construction parameters (fracturing fluid volume and pump stop pressure) is established. The flow of gas and water in the early stage after fracturing is simulated, and the regional phase permeability and capillary force curves are introduced to simulate the process of flowback and production of horizontal wells after fracturing. The influence of early fracture closure on the gas-water flow is characterized by stress sensitivity. A deep shale gas reservoir of Sinopec was selected for the case study. The simulation results show it necessary to consider the effects of fractures and stress sensitivity in the matrix when considering the dynamic change of production during the flowback and production stages. The findings of this study can help for better understanding of the fracture distribution characteristics of shale gas, shale gas production principle, and well EUR prediction, which provide a theoretical basis for the effective development of shale gas horizontal well groups.

Numerical Study of Reactive Flow in Fractured Carbonate Rock

Acidizing technology is an effective reformation method of oil and gas reservoirs. It can also remove the reservoir pollution near wellbore zones and enhance the fluid transmissibility. The optimal injection rate of acid is one of the key factors to reduce cost and improve the effect of acidizing. Therefore, the key issue is to find the optimal injection rate during acid corrosion in fractured carbonate rock. In this work, a novel reactive flow mathematical model based on two-scale model and discrete fracture model is established for fractured carbonate reservoirs. The matrix and fracture are described by a two-scale model and a discrete fracture model, respectively. Firstly, the two-scale model for matrix is combined with the discrete fracture model. Then, an efficient numerical scheme based on the finite element method is implemented to solve the corresponding dimensionless equations. Finally, several important aspects, such as the influence of the injection rate of acid on the dissolution patterns, the influence of fracture aperture and fracture orientations on the dissolution structure, the breakthrough volume of injected acid, and the dynamic change of fracture aperture during acidizing, are analyzed. The numerical simulation results show that there is an optimal injection rate in fractured carbonate rock. However, the fractures do not have an impact on the optimal acid injection rate, they only have an impact on the dissolution structure.

Study on Smart Numerical Simulation Mechanism in Complex Reservoir Using Embedded Discrete Fracture Model

The embedded discrete fracture model (EDFM) method is used to study the mechanism model of the reservoir with a complex fracture system, and the sensitivity analysis of the development effect is carried out. The results show that: 1) the EDFM method can directly characterize the fracture morphology and fracture attributes in the reservoir; 2) Through comparative analysis of reservoir type, capillary force, and fracture development degree, it is found that different reservoir types have great differences in development effect, fracture opening has little influence on the reservoir, capillary force and fracture density have a great influence on the reservoir development effect, and with the increase of fracture density, the reservoir development effect increases significantly.

Unified Reservoir And Seismic Simulation With Explicit Representation Of Fractures And Faults

In the context of remote sensing, the vast disparity in characteristic scales between seismic deformation (e.g. milliseconds) and transient flow (e.g. hours) allows a "two-model paradigm" for geophysics and reservoir simulation. In the context of flow-induced geohazard risk mitigation and micro-seismic data integration, this paradigm breaks down. Under micro-seismic deformation, events occur with high-frequency, and over sustained duration during which the rock-fluid coupling is significant. In risk mitigation scenarios, the onset of seismic deformation is directly tied to quasi-static coupling periods. This work develops an approach to reservoir simulation modeling that allows simultaneous resolution of transient (inertial) poromechanics and multiphase fluid flow in the presence of fracture. A mixed discretization scheme combining the extended finite element method (XFEM) and the embedded discrete fracture model (EDFM) is extended using a second-order implicit Newmark time integration scheme for the inertial mechanics. A Lagrange multiplier method is developed to model pressure-dependent contact traction in fractures. The contact constraints are adapted to accommodate fracture opening. Slip-weakening fracture friction models are incorporated. Finally, a time-step controller is proposed to combine local discretization error with contact traction and slip-rate control along the fractures. This strategy allows automatic adaptation to resolve quasi-static, inter-seismic triggering, and co-seismic spontaneous rupture periods within one model. The model is verified to simulate complete induced earthquake sequences, including inter-seismic and dynamic rupture phases. The performance of the adaptive model is illustrated for cases with various set-ups of production and injection periods in a fractured reservoir with explicit fracture representation.

Constructing Dual-Porosity Models from High-Resolution Discrete-Fracture Models Using Deep Neural Networks

Detailed geological description of fractured reservoirs is typically characterized by the discrete-fracture model (DFM), in which the rock matrix and fractures are explicitly represented in the form of unstructured grids. Its high computation cost makes it infeasible for field-scale applications. Traditional flow-based and static-based methods used to upscale detailed geological DFM to reservoir simulation model suffer from, to some extent, high computation cost and low accuracy, respectively. In this paper, we present a novel deep learning-based upscaling method as an alternative to traditional methods. This work aims to build an image-to-value model based on convolutional neural network to model the nonlinear mapping between the high-resolution image of detailed DFM as input and the upscaled reservoir simulation model as output. The reservoir simulation model (herein refers to the dual-porosity model) includes the predicted fracture-fracture transmissibility linking two adjacent grid blocks and fracture-matrix transmissibility within each coarse block. The proposed upscaling workflow comprises the train-validation samples generation, convolutional neural network training-validating process, and model evaluation. We apply a two-point flux approximation (TPFA) scheme based on embedded discrete-fracture model (EDFM) to generate the datasets. We perform trial-error analysis on the coupling training-validating process to update the ratio of train-validation samples, optimize the learning rate and the network architecture. This process is applied until the trained model obtains an accuracy above 90 % for both train-validation samples. We then demonstrate its performance with the two-phase reference solutions obtained from the fine model in terms of water saturation profile and oil recovery versus PVI. Results show that the DL-based approach provides a good match with the reference solutions for both water saturation distribution and oil recovery curve. This work manifests the value of the DL-based method for the upscaling of detailed DFM to the dual-porosity model and can be extended to construct generalized dual-porosity, dual-permeability models or include more complex physics, such as capillary and gravity effects.

Efficient Modeling of Hydraulic Fractures with Diamond Boundary in Shale Gas Reservoirs

The effects of hydraulic fractures with complex boundaries on the well performance of a shale gas well, considering a more realistic corner point geological model, have rarely been studied previously. In this study, the nonintrusive embedded discrete fracture model (EDFM) method was employed to investigate the effects of different boundary shapes of hydraulic fracture, coupling a network of sophisticated natural fractures discrete fracture network (DFN) on well performance. First, by implementing the powerful EDFM technology, concepts of two categories (rectangle and diamond) of hydraulic fracture with different boundaries were designed. Next, the geometric equations defining vertices of multiple rectangular- or diamond-shaped hydraulic fractures in arbitrary coordinate systems were derived. Subsequently, the horizontal well with multistaged hydraulic fractures and sophistically oriented 3D natural fractures was inputted into the reservoir model to perform history matching. After history matching, the results were further analyzed to compare the production forecast from the two categories. The results show that 20-year cumulative gas productions for rectangle- and diamond-shaped fractures are approximately 1.237×108 m3 and 1.486×108 m3, respectively. In other words, the diamond category can produce 20.1% more gas than the rectangle category. For cumulative water production, the diamond category produces 3.8×104 m3, as against the 3.0×104 m3 produced by the rectangle category (or 26.7% more). This implies that the diamond-shaped fractures have the potential to reach the far field region of the reservoir away from the wellbore. This means that more intersections with natural fractures DFN can be achieved, and more drainage area is unlocked. The visualization of pressure distributions and drainage volume was easily shown, and these results further confirm that the extent of fluid drainage by the diamond fracture is larger compared to that by the rectangular fracture given the same total surface area.