{'fraction_non_alphanumeric': 0.035269709543568464, 'fraction_numerical': 0.0005186721991701245, 'mean_word_length': 5.473154362416107, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
An Efficient Numerical Hybrid Model for Multiphase Flow in Deformable Fractured-Shale Reservoirs


 After hydraulic fracturing, a shale reservoir usually has multiscale fractures and becomes more stress-sensitive. In this work, an adaptive hybrid model is proposed to simulate hydromechanical coupling processes in such fractured-shale reservoirs during the production period (i.e., the hydraulic-fracturing process is not considered and cannot be simulated). In our hybrid model, the single-porosity model is applied in the region outside the stimulated reservoir volume (SRV), and the matrix and natural/induced fractures in the SRV region are modeled using a double-porosity model that can accurately simulate the matrix/fracture fluid exchange during the entire transient period. Meanwhile, the fluid flow in hydraulic fractures is modeled explicitly with the embedded-discrete-fracture model (EDFM), and a stabilized extended-finite-element-method (XFEM) formulation using the polynomial-pressure-projection (PPP) technique is applied to simulate mechanical processes. The developed stabilized XFEM formulation can avoid the displacement oscillation on hydraulic-fracture interfaces. Then a modified fixed-stress sequential-implicit method is applied to solve the hybrid model, in which mixed-space discretization [i.e., finite-volume method (FVM) for flow process and stabilized XFEM for geomechanics] is used. The robustness of the proposed model is demonstrated through several numerical examples. In conclusion, several key factors for gas exploitation are investigated, such as adsorption, Klinkenberg effect, capillary pressure, and fracture deformation. In this study, all the numerical examples are 2D, and the gravity effect is neglected in these simulations. In addition, we assume there is no oil phase in the shale reservoirs, thus the gas/water two-phase model is used to simulate the flow in these reservoirs.