Total: 1
Simulating high-enthalpy fractured geothermal reservoirs is challenging due to the complex coupled processes of non-isothermal, multiphase, multicomponent flow, strongly nonlinear thermodynamics, and the dominant role of fractures. These complexities are amplified by mineral scaling, such as halite precipitation, which can impair reservoir permeability and well productivity. To address this, we present a new compositional flow model based on a persistent set of primary variables (pressure, enthalpy, and overall salt mass fraction). The formulation naturally handles phase transitions without manual switching, enhancing numerical stability. The model integrates a discrete fracture-matrix approach and employs an efficient, robust correlation-based phase-behaviour linearisation of saltwater thermodynamics, replacing expensive on-the-fly phase separation calculations. It incorporates the Kozeny-Carman relation to dynamically model porosity and permeability reduction from halite precipitation. Implemented in the open-source PorePy framework, the model is verified through a 1D salt dissolution benchmark against the established closed-source simulator CSMP++, showing strong agreement across geothermal conditions involving transitions between single- and multi-phase regions. Application to a 2D halite-saturated fractured reservoir with injection and production demonstrates the model's capability to predict halite precipitation patterns and their impact on permeability damage and energy recovery. Numerical results further show the model's value in predicting operational challenges such as wellbore blockage and the role of fracture connectivity. The model thus provides an open-source numerical tool for analysing complex heat and mass transport with mineral scaling in high-enthalpy fractured geothermal systems.