A Semi-Analytical Geomechanical Approach for Forecasting Production Performance in Multifractured Composite Systems

Stress-dependence of reservoir matrix and fractures can strongly affect the performance of multifractured horizontal wells (MFHWs) completed in unconventional hydrocarbon reservoirs. In order to model fluid flow in unconventional reservoirs exhibiting this stress-dependence, most traditional reservoir flow simulators, and many simulators described in published work, use conventional reservoir fluid flow model formulations. These formulations typically neglect the influence of the rate of change of volumetric strain of the reservoir matrix and fractures, even though reservoir stress and pressure change significantly during the course of production. As a result, the effect of matrix and fracture deformation on production is neglected, which can lead to errors in predicting production performance in most stress-sensitive reservoirs. To address this problem, some studies have proposed the use of porosity and transmissibility multipliers to model stress-sensitive reservoirs. However, in order to apply this approach, multipliers must be estimated from laboratory experiments, or used as a history-match parameter, possibly resulting in large errors in well performance predictions. Alternatively, fully-coupled, fully numerical geomechanical simulation can be performed, but these methods are computationally costly, and models are difficult to setup. This paper presents a new fully-coupled, two-way analytical modeling approach that can be used to simulate fluid flow in stress-sensitive unconventional reservoirs produced through MFHWs. The model couples poroelastic geomechanics theory with fluid flow formulations. The two-way coupled fluid flow-geomechanical analytical model is applied simultaneously to both the matrix and fracture regions. In the proposed algorithm, a porosity-compressibility coupling parameter for the two physical models is setup to update the stress- and pressure-dependent fracture/matrix properties iteratively, which are later used as input data for the fracture-matrix reservoir fluid flow model at each iteration step. The analytical approach developed for the fully-coupled, two-way analytical model, using the enhanced fracture region conceptual model, is validated by comparing the results with numerical simulation. Predictions using the fully-coupled enhanced fracture region model are then compared with the same enhanced fracture region model but with the conventional pressure-dependent modeling approach implemented. A sensitivity study performed by comparing the new fully-coupled model predictions with and without geomechanics effects accounted for reveals that, without geomechanics effects, production performance in stress-sensitive reservoirs might be overestimated. The study also demonstrates that use of the conventional stress-dependent modeling approach may cause production performance to be underestimated. Therefore, the proposed fully-coupled, two-way analytical model can be useful for practical engineering purposes.