Assessing 3D Bedding Plane Effects on Well Performance in a Field-Scale Unconventional Reservoir Model Using EDFM

Abstract Bedding-plane slip effects during hydraulic fracturing have recently gained interest in unconventional plays due to their influence in hydraulic fracture growth in vertical and horizontal directions. However, most of the current workflows cannot fully model field-scale sub-horizontal orientation of bedding planes because of complications with gridding techniques, or due to simplifications related to the use of 2D models. These challenges have motivated the assessment of 3D bedding plane interactions on well performance using the embedded discrete fracture model (EDFM) technology for field case scenarios. An efficient hydraulic fracture propagation model is used to model hydraulic fracture growth in the presence of bedding layers. The model captures shear slippage at the bedding layer interfaces and corrects the calculated stress intensity factor to account for height containment. A hydraulic fracture model, constrained by geomechanical information, is built in a corner point grid. Resulting hydraulic fracture geometries and identified bedding layer fractures are transferred to EDFM by using a 3D bedding plane generator, which places sub-horizontal polygons across the well trajectory, honoring its orientation and geometry. To locate the spatial position of bedding layers, geostatistical constrains, core analysis and petrophysical interpretations – including well image logs – can be taken into account. Lastly, a reservoir simulation model is built to evaluate the effects of bedding planes on well performance. 3D effects of bedding planes in a shale gas reservoir were captured in a field case scenario using numerical models. Higher contribution to production was observed in the results of this study. The main reasons are larger fracture lengths generated along the pay zone caused by bedding plane influence in the fracture propagation process and shear slippage along bedding plane fractures, which create a larger effective conductive surface area. When modeling bedding planes, computational efficiency is substantial due to the EDFM method, preserving spatial orientation and geometry of each bedding plane. Direct assessment of bedding plane properties is provided, which highlights the importance of capturing their interactions with hydraulic fracture growth and well performance. A seamless integration of bedding plane models can be achieved in an efficient workflow that provides key lessons for future fracture design and well spacing optimization.