Temperature Transient Analysis of Naturally Fractured Geothermal Reservoirs

A potential approach to develop geothermal energy is by producing low - to-medium temperature fluids from naturally fractured geother-mal reservoirs (NFGRs). Pressure transient analysis (PTA) is the most common approach to characterize such reservoirs for improv-ing development efficiency. However, pressure inversion leads to nonuniqueness and cannot be used to estimate thermal properties. Moreover, reliable methods to evaluate the development potential of fractured geothermal reservoirs are lacking. To address the gap, this work aims to study the temperature behavior and explore a suitable analysis method for characterizing geothermal reservoirs and evaluating development potential. We developed numerical and analytical models to analyze the temperature behavior in NFGRs. The developed models account for the Joule-Thomson [J -T effect (mu(JT))], adiabatic heat expansion/compression effect (zeta), reservoir formation damage, heat conduction, and convection effects. The developed numerical solution is verified and found to agree with the proposed analytical solutions. The results show that temperature transient analysis (TTA) with constant or temperature-dependent mu(JT) and zeta assumption leads to a minor difference when reservoir temperature changes significantly. Moreover, three heat radial flow regimes (HRFR) and a thermal interporosity regime with a V -shape characteristic have been identified. The results also show that temperature data provide information not accessible by PTA. The results reveal that temperature derivative curves signify a "hump " when formation around a wellbore is dam-aged, and the temperature data can be used to characterize the skin -zone radius and permeability. It is demonstrated that the properties such as J -T coefficient, effective adiabatic heat expansion coefficient, and fracture intrinsic porosity can be estimated using TTA. The results indicate that fracture thermal storativity (omega(T)) and matrix thermal interporosity coefficient (alpha(T)) can be estimated from the thermal interporosity regime exhibited on the temperature derivative curve. The results also suggest that commercial geothermal energy harness is more difficult when the omega T is high or the alpha T is very small. Finally, we introduced an integrated workflow of combining PTA and TTA to characterize NFGRs. Simulated test examples are interpreted to demonstrate the applicability of the developed workflow. This work aids in better understanding the potentials of temperature data on geothermal reservoir characterization.