An Improved Grid Generation Approach for Discrete Fracture Network Modelling Using Line Fracture Concept for Two-Phase Flow Simulation
- Keke Wang (SouthWest Petroleum University) | Xiaolong Peng (SouthWest Petroleum University) | Zhimin Du (SouthWest Petroleum University) | Manicher Haghighi (Zhenjiang You) | Lu Yu (Notheast Sichuan Gas production plant, Southwest Branch company, Sinopec Corporation)
- Document ID
- Society of Petroleum Engineers
- SPE Asia Pacific Oil & Gas Conference and Exhibition, 25-27 October, Perth, Australia
- Publication Date
- Document Type
- Conference Paper
- 2016. Society of Petroleum Engineers
- 5 Reservoir Desciption & Dynamics, 5.5 Reservoir Simulation, 5.8 Unconventional and Complex Reservoirs, 4.1.2 Separation and Treating, 4 Facilities Design, Construction and Operation, 5.8.6 Naturally Fractured Reservoir, 4.1 Processing Systems and Design
- discrete fracture network model, line fracture, MPFA method, two-phase flow, unstructured quadrilateral/hexahedral mesh
- 1 in the last 30 days
- 117 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 9.50|
|SPE Non-Member Price:||USD 28.00|
In traditional Discrete Fracture Network (DFN) method for fractured porous media simulation, fractures are modelled in different scales and orientations using single computational grid. However, a large number of small grids near both ends of fracture are generated, which substantially reduce the computation efficiency. In this paper, a new line fracture approach is presented to avoid generation of small grids. Combined with Multi-Point Flux Approximation (MPFA) algorithm, the simulation efficiency is highly improved without losing accuracy.
We present a modified mesh generation algorithm in which the volume of large fractures is set zero by neglecting the fracture width. The fracture aperture (very short line) is treated as single node during mesh generation using the paving method. Therefore, the number of small grids around the fracture is significantly decreased. In this paper, a quadrilateral unstructured grid system is generated using the proposed paving method. Then, by taking the median dual of quadrilateral mesh as the control volume, all fracture properties such as the aperture and permeability are accounted for as independent parameters in discrete fractures. The governing equations of DFN model are derived and solved using MPFA algorithm for oil-water two-phase flow.
A new discrete fracture network model using line fracture approach is developed. To show the efficiency and accuracy of the new model, the proposed method is applied to a case study of two separated and two crossed fractures, in the computation domain with the scale of 50 by 50 meters. Compared with the standard DFN model including the width of fractures, the computation efficiency of the new model is highly improved (by 2.5 times), due to the large reduction of grid numbers (20%). Good agreement in water and oil saturation distributions resulted from the two models validates the newly developed method.
Simulation of naturally fractured reservoirs has been a challenge for decades. Discrete fracture network modelling is a promising approach. However, without an efficient grid generation method, DFN model is not feasible and practical due to the large amount of computational time. The novelty of the present work is the improved model using line fracture approach. The new grid generation algorithm is coupled with multi-point flux approximation approach for flow simulation in fracture network system. High efficiency and accuracy of the developed method is readily applied to fractured reservoir simulation.
Sponsored by NSFC (Natural Science Fund of China - "the numerical simulation of natural fractured reservoir based on the unstructured hexahedral grid system and concurrent computation") (No. 51474179)
|File Size||5 MB||Number of Pages||19|
Castellini, A. Flow Based Grids for Reservoir Simulation MS thesis, Stanford University, Stanford, CA USA.Courrioux, G., Nullans, S., Guillen, A., Boissonnat, J.D., Repusseau, P., Renaud, X. and Thibaut, M. 3D volumetric modelling of Cadomian terranes (Northern Brittany, France): an automatic method using Voronoi diagrams. Tectonophysics, 2001, 331, 181–196.
Diskin, B., Thomas, J., Nielsen, E., White, J. and Nishikawa, H.  Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations Part I: Viscous Fluxes. 47th AIAA Aerospace Sciences Meeting including The New Horizons Forum and Aerospace Exposition. 5-8 January Orlando Florida.
Sun, J., Schechter, D. S. 2015a. Optimization-Based Unstructured Meshing Algorithms for Simulation of Hydraulically and Naturally Fractured Reservoirs with Variable Distribution of Fracture Aperture, Spacing, Length and Strike. SPE Reservoir Evaluation & Engineering 18 (04):463–480. SPE-170703-PA.
Merland, R., Lévy, B. and Caumon G. Building PEBI Grids Conforming to 3D Geological Features Using Centroidal Voronoi Tessellations. In: Marschallinger, R. and Zolb, R. (Eds.) Mathematical Geosciences at the Crossroads of Theory and Practice – Proceedings of IAMG conference. Slazburg, 2011, 254–265.