Preventing frac hits and well interferences with fast marching simulation using embedded discrete fracture models constrained by poroelastic geomechanical modelling of enhanced permeability

First Break ◽  
2019 ◽  
Vol 37 (4) ◽  
pp. 49-54
Author(s):  
A. Ouenes ◽  
A. Bachir ◽  
R. Smaoui ◽  
C. Hammerquist ◽  
M. Paryani
SPE Journal ◽  
2019 ◽  
Vol 24 (06) ◽  
pp. 2590-2608 ◽  
Author(s):  
Xu Xue ◽  
Changdong Yang ◽  
Tsubasa Onishi ◽  
Michael J. King ◽  
Akhil Datta–Gupta

Summary Recently, fast–marching–method (FMM) –based flow simulation has shown great promise for rapid modeling of unconventional oil and gas reservoirs. Currently, the application of FMM–based simulation has been limited to using tartan grids to model the hydraulic fractures (HFs). The use of tartan grids adversely impacts the computational efficiency, particularly for field–scale applications with hundreds of HFs. Our purpose in this paper is to extend FMM–based simulation to incorporate local grid refinements (LGRs) and an embedded discrete fracture model (EDFM) to simulate HFs with natural fractures, and to validate the accuracy and efficiency of the methodologies. The FMM–based simulation is extended to LGRs and EDFM. This requires novel gridding through the introduction of triangles (2D) and tetrahedrons (2.5D) to link the local and global domain and solution of the Eikonal equation in unstructured grids to compute the diffusive time of flight (DTOF). The FMM–based flow simulation reduces a 3D simulation to an equivalent 1D simulation using the DTOF as a spatial coordinate. The 1D simulation can be carried out using a standard finite–difference method, thus leading to orders of magnitude of savings in computation time compared with full 3D simulation for high–resolution models. First, we validate the accuracy and computational efficiency of the FMM–based simulation with LGRs by comparing them with tartan grids. The results show good agreement and the FMM–based simulation with LGRs shows significant improvement in computational efficiency. Then, we apply the FMM–based simulation with LGRs to the case of a multistage–hydraulic–fractured horizontal well with multiphase flow, to demonstrate the practical feasibility of our proposed approach. After that, we investigate various discretization schemes for the transition between the local and global domain in the FMM–based flow simulation. The results are used to identify optimal gridding schemes to maintain accuracy while improving computational efficiency. Finally, we demonstrate the workflow of the FMM–based simulation with EDFM, including grid generation; comparison with FMM with unstructured grid; and validation of the results. The FMM with EDFM can simulate arbitrary fracture patterns without simplification and shows good accuracy and efficiency. This is the first study to apply the FMM–based flow simulation with LGRs and EDFM. The three main contributions of the proposed methodology are (i) unique mesh–generation schemes to link fracture and matrix flow domains, (ii) DTOF calculations in locally refined grids, and (iii) sensitivity studies to identify optimal discretization schemes for the FMM–based simulation.


Geofluids ◽  
2020 ◽  
Vol 2020 ◽  
pp. 1-12
Author(s):  
Tao Chen

Equivalent fracture models are widely used for simulations of groundwater exploitation, geothermal reservoir production, and solute transport in groundwater systems. Equivalent permeability has a great impact on such processes. In this study, equivalent permeability distributions are investigated based on a state-of-the-art numerical upscaling method (i.e., the multiple boundary method) for fractured porous rocks. An ensemble of discrete fracture models is created based on power law length distributions. The equivalent permeability is upscaled from discrete fracture models based on the multiple boundary method. The results show that the statistical distributions of equivalent permeability tensor components are highly related to fracture geometry and differ from each other. For the histograms of the equivalent permeability, the shapes of k x x and k y y change from a power law-like distribution to a lognormal-like distribution when the fracture length and the number of fractures increase. For the off-diagonal component k x y , it is a normal-like distribution and its range expands when the fracture length and the number of fractures increase. The mean of diagonal equivalent permeability tensor components increases linearly with the fracture density. The analysis helps in generating stochastic equivalent permeability models in fractured porous rocks and reduces uncertainties when applying equivalent fracture models.


Geofluids ◽  
2020 ◽  
Vol 2020 ◽  
pp. 1-20
Author(s):  
Aaditya Khanal ◽  
Ruud Weijermars

The present study compares flow paths in reservoirs with natural fractures, solved with Complex Analysis Methods (CAM), to those solved with Embedded Discrete Fracture Models (EDFM). One aim is to define scaling rules for the strength (flux) of the discrete natural fractures used in CAM models, which was previously theoretically defined based on the expected flow distortion. A major hurdle for quantitative benchmarks of CAM with EDFM results is that each of the two methods accounts for natural fractures with different assumptions and input parameters. For example, EDFM scales the permeability of the natural fractures based on a cubic equation, while CAM uses a flux strength. The results from CAM and EDFM are used to scale the flux strength of the natural fractures and improve the equivalent permeability contrast estimation for CAM. The permeability contrast for CAM is calculated from the ratio of the enhanced velocity inside natural fractures to the unperturbed matrix fluid velocity. A significant advantage of flow and pressure models based on CAM is the high resolution without complex gridding. Particle tracking results are presented for fractures with different hydraulic conductivity ranging from highly permeable to impervious.


1995 ◽  
Vol 412 ◽  
Author(s):  
John J. Nitao ◽  
Thomas A. Buscheck

AbstractAn in situ heater test was performed at G-Tunnel, Nevada Nuclear Test Site, to investigate the thermal-hydrological response of unsaturated, fractured volcanic tuff under conditions similar to those at Yucca Mountain. The NUFT flow and transport code was used to model the test using discrete-fracture and equivalent-continuum approaches. Nonequilibrium fracture flow and thermal buoyant gas-phase convection were found to be the likely causes for observed lack of condensate imbibition into the matrix. The potential repository at Yucca Mountain was also modeled. Disequilibrium fracture flow is predicted to occur for less than a hundred years after emplacement followed by a period of fracture-matrix equilibrium, during which the equivalent-continuum and discrete-fracture models give almost identical results.


2020 ◽  
Vol 8 ◽  
Author(s):  
Zhi-dong Yang ◽  
Yong Wang ◽  
Xu-yang Zhang ◽  
Ming Qin ◽  
Shao-wen Su ◽  
...  

Tight oil is a kind of unconventional oil and gas resource with great development potential. Due to the unconventional characteristics of low porosity and low permeability in tight oil reservoirs, single wells generally have no natural productivity, and industrial development is usually conducted in combination with horizontal wells and hydraulic fracturing techniques. To capture the flow behavior affected by fractures with complex geometry and interaction, we adopted embedded discrete fracture models (EDFMs) to simulate the development of fractured reservoirs. Compared with the traditional discrete fracture models (DFMs), the embedded discrete fracture models (EDFMs) can not only accurately represent the fracture geometry but also do not generate a large number of refine grids around fractures and intersections of fractures, which shows the high computational efficiency. To be more consistent with the real characteristic of the reservoir and reflect the advantage of EDFMs on modeling complex fractures, in this work, the hydraulic fractures are set as oval shape, and we adopted 3-dimensional oil–gas two-phase model considering capillary forces and gravity effects. We developed an EDFM simulator, which is verified by using the fine grid method (FGM). Finally, we simulated and studied the development of tight oil without and with random natural fractures (NFs). In our simulation, the pressure varies widely from the beginning to the end of the development. In real situation, tight oil reservoirs have high initial pressure and adopt step-down bottom hole pressure development strategy where the bottom hole pressure of the last stage is below the bubble point pressure and the free gas appears in the reservoir. Modeling studies indicate that the geometry of fracture has a great influence on the pressure and saturation profiles in the area near the fractures, and dissolved gas flooding contributes to the development of tight oil, and NFs can significantly improve production, while the effect of the stress sensitivity coefficient of NFs on production is more significant in the later stage of production with lower reservoir pressure.


Sign in / Sign up

Export Citation Format

Share Document