Mathematical Modeling and Numerical Simulation of Water‐Rock Interaction in Shale Under Fracturing‐Fluid Flowback Conditions

Abstract

Water‐rock interaction cannot be ignored for shale reservoirs with high‐salinity formation brine and complex rock composition, and stimulated by massive slick‐water fracturing treatment. However, there have been few studies on the flowback model fully coupled with different effects of water‐rock interaction. This paper presents the development of a coupled hydro‐chemical‐mechanical model for modeling water‐rock interaction in fractured shale during the post‐fracturing flowback period. The model considers distinguishing water‐rock interaction phenomena, that is, mineral dissolution, clay swelling and chemical osmosis, and accounts for multi‐phase flow in a fractured shale reservoir. The coupling and solution method and a numerical simulator were developed. Numerical simulation indicates that the swelling volume of clay minerals occupies the pores and leads to a decline in matrix porosity, while mineral dissolution increases both the matrix porosity and the solute concentration in the aqueous phase in matrix pores. Clay swelling mainly affects the shape of the porosity ratio profiles. The effect of mineral dissolution becomes increasingly stronger as flowback progresses. Mineral dissolution mainly affects the relative positions of the porosity ratio curves with the progress of flowback. The water‐rock interaction coupled flowback modeling and the numerical simulation results in this study quantify the effects of chemical osmosis, clay swelling and mineral dissolution. Results from this study provide new insights into the mechanisms of fracturing‐fluid flowback and value to flowback transient analysis.

Correlation between Natural Seepage and Ionic Diffusion in Shales

“…Chen (2021) carried out mathematical modeling of shale water-rock interactions under fracturing fluid retentate conditions, and this model was able to quantitatively analyze the extent of the effects of chemical permeability, clay expansion, and mineral dissolution. A multiphase flow algorithm was incorporated . Liu (2022) initiated a particle tracking method study based on the proppant flowback effect during fracturing fluid flowback, which eventually led to an improved steady-state Navier–Stokes equation approach .…”

Experimental Water Activity Suppression and Numerical Simulation of Shale Pore Blocking

“…The results show that the swelling volume of clay minerals occupies the pore space, leading to a decrease in matrix porosity, while mineral dissolution increases matrix porosity and solute concentration of the aqueous phase in the matrix pore space. Clay expansion mainly affected the shape of the porosity curve [18]. In a study conducted by Wang Yepeng [19], an investigation of water-shale interactions was carried out.…”