Coupled THM and Matrix Stability Modeling of Hydroshearing Stimulation in a Coupled Fracture-Matrix Hot Volcanic System

A coupled thermal-hydraulic-mechanical (THM) model is developed to simulate the combined effect of fracture fluid flow, heat transfer from the matrix to injected fluid, and shearing dilation behaviors in a coupled fracture-matrix hot volcanic reservoir system. Fluid flows in the fracture are calculated based on the cubic law. Heat transfer within the fracture involved is thermal conduction, thermal advection, and thermal dispersion; within the reservoir matrix, thermal conduction is the only mode of heat transfer. In view of the expansion of the fracture network, deformation and thermal-induced stress model are added to the matrix node’s in situ stress environment in each time step to analyze the stability of the matrix. A series of results from the coupled THM model, induced stress, and matrix stability indicate that thermal-induced aperture plays a dominant role near the injection well to enhance the conductivity of the fracture. Away from the injection well, the conductivity of the fracture is contributed by shear dilation. The induced stress has the maximum value at the injection point; the deformation-induced stress has large value with smaller affected range; on the contrary, thermal-induced stress has small value with larger affected range. Matrix stability simulation results indicate that the stability of the matrix nodes may be destroyed; this mechanism is helpful to create complex fracture networks.