Fault zones exhibit 3D variable thickness, a feature that remains inadequately explored, particularly with regard to the impact on fluid flow. Upon analyzing an analytic solution, we examine 3D thermal-hydraulic (TH) dynamical models through a benchmark experiment, which incorporates a fault zone with thickness variations corresponding to realistic orders of magnitude. The findings emphasize an area of interest where vigorous convection drives fluid flow, resulting in a temperature increase to 150°C at a shallow depth of 2.7 km in the thickest sections of the fault zone. Moreover, by considering various tectonic regimes (compressional, extensional, and strike-slip) within 3D thermal-hydraulic-mechanical (THM) models and comparing them to the benchmark experiment, we observe variations in fluid pressure induced by poroelastic forces acting on fluid flow within the area of interest. These tectonic-induced pressure changes influence the thermal distribution of the region and the intensity of temperature anomalies. Outcomes of this study emphasize the impact of poroelasticity-driven forces on transfer processes and highlight the importance of addressing fault geometry as a crucial parameter in future investigations of fluid flow in fractured systems. Such research has relevant applications in geothermal energy, CO2 storage, and mineral deposits.
