Alfredo Battistelli

3D Voronoi grid dedicated software for modeling gas migration in deep layered sedimentary formations with TOUGH2-TMGAS

As is known, a full three-dimensional (3D) unstructured grid permits a great degree of flexibility when performing accurate numerical reservoir simulations. However, when the Integral Finite Difference Method (IFDM) is used for spatial discretization, constraints (arising from the required orthogonality between the segment connecting the blocks nodes and the interface area between blocks) pose difficulties in the creation of grids with irregular shaped blocks. The full 3D Voronoi approach guarantees the respect of IFDM constraints and allows generation of grids conforming to geological formations and structural objects and at the same time higher grid resolution in volumes of interest. In this work, we present dedicated pre- and post-processing gridding software tools for the TOUGH family of numerical reservoir simulators, developed by the Geothermal Research Group of the DICAM Department, University of Bologna. VORO2MESH is a new software coded in C++, based on the voro++ library, allowing computation of the 3D Voronoi tessellation for a given domain and the creation of a ready to use TOUGH2 MESH file. If a set of geological surfaces is available, the software can directly generate the set of Voronoi seed points used for tessellation. In order to reduce the number of connections and so to decrease computation time, VORO2MESH can produce a mixed grid with regular blocks (orthogonal prisms) and irregular blocks (polyhedron Voronoi blocks) at the point of contact between different geological formations. In order to visualize 3D Voronoi grids together with the results of numerical simulations, the functionality of the TOUGH2Viewer post-processor has been extended. We describe an application of VORO2MESH and TOUGH2Viewer to validate the two tools. The case study deals with the simulation of the migration of gases in deep layered sedimentary formations at basin scale using TOUGH2-TMGAS. A comparison between the simulation performances of unstructured and structured grids is presented. VORO2MESH is a 3D Voronoi gridding tool for TOUGH2.VORO2MESH allows generation of 3D Voronoi grids in compliance with the geological shape.TOUGH2Viewer is improved for 3D Voronoi grid visualization.

Coupled Processes Modeling in Rock Salt and Crushed Salt Including Halite Solubility Constraints: Application to Disposal of Heat-Generating Nuclear Waste

This paper presents numerical modeling of coupled thermal, hydraulic and mechanical processes in rock salt and crushed salt considering halite solubility constraints. The TOUGH-FLAC simulator is used, with a recently enhanced Equation-Of-State module that includes the thermodynamic properties of aqueous fluids of variable salinity. Laboratory and field scale tests performed on rock salt and crushed salt under temperature gradients are modeled first to evaluate the capabilities of the simulator to reproduce important features, such as porosity changes induced by halite dissolution/precipitation, and brine and heat migration. Since the results are quite satisfactory, the simulator is used to predict the long-term response of a generic salt repository for heat-generating nuclear waste. To evaluate the impacts of halite solubility on the predictions, two simulations that respectively consider or neglect solubility constraints are performed. In the scenario studied, the results are not significantly affected by dissolution/precipitation, and only some differences are observed due to changes in porosity, but the dominating processes remain the same. With the new provisions, TOUGH-FLAC is more complete in terms of processes occurring around a heat-releasing nuclear waste package and can therefore provide more accurate predictions of the long-term performance of a nuclear waste repository in salt formations.