Abstract. Petrophysical and mechanical rock properties are key parameters for the characterization of the deep subsurface in different disciplines such as geothermal heat extraction, petroleum reservoir engineering or mining. They are commonly used for the interpretation of geophysical data and the parameterization of numerical models and thus are the basis for economic reservoir assessment. However, detailed information regarding petrophysical and mechanical rock properties for each relevant target horizon is often scarce, inconsistent or distributed over multiple publications. Therefore, subsurface models are often populated with generalized or assumed values resulting in high uncertainties. Furthermore, diagenetic, metamorphic and hydrothermal processes significantly affect the physiochemical and mechanical properties often leading to high geological variability. A sound understanding of the controlling factors is needed to identify statistical and causal relationships between the properties as a basis for a profound reservoir assessment and modeling. Within the scope of the GEMex project (EU H2020, grant agreement no. 727550), which aims to develop new transferable exploration and exploitation approaches for enhanced and super-hot unconventional geothermal systems, a new workflow was applied to overcome the gap of knowledge of the reservoir properties. Two caldera complexes located in the northeastern Trans-Mexican Volcanic Belt – the Acoculco and Los Humeros caldera – were selected as demonstration sites. The workflow starts with outcrop analog and reservoir core sample studies in order to define and characterize the properties of all key units from the basement to the cap rock as well as their mineralogy and geochemistry. This allows the identification of geological heterogeneities on different scales (outcrop analysis, representative rock samples, thin sections and chemical analysis) enabling a profound reservoir property prediction. More than 300 rock samples were taken from representative outcrops inside the Los Humeros and Acoculco calderas and the surrounding areas and from exhumed “fossil systems” in Las Minas and Zacatlán. Additionally, 66 core samples from 16 wells of the Los Humeros geothermal field and 8 core samples from well EAC1 of the Acoculco geothermal field were collected. Samples were analyzed for particle and bulk density, porosity, permeability, thermal conductivity, thermal diffusivity, and heat capacity, as well as ultrasonic wave velocities, magnetic susceptibility and electric resistivity. Afterwards, destructive rock mechanical tests (point load tests, uniaxial and triaxial tests) were conducted to determine tensile strength, uniaxial compressive strength, Young's modulus, Poisson's ratio, the bulk modulus, the shear modulus, fracture toughness, cohesion and the friction angle. In addition, X-ray diffraction (XRD) and X-ray fluorescence (XRF) analyses were performed on 137 samples to provide information about the mineral assemblage, bulk geochemistry and the intensity of hydrothermal alteration. An extensive rock property database was created (Weydt et al., 2020; https://doi.org/10.25534/tudatalib-201.10), comprising 34 parameters determined on more than 2160 plugs. More than 31 000 data entries were compiled covering volcanic, sedimentary, metamorphic and igneous rocks from different ages (Jurassic to Holocene), thus facilitating a wide field of applications regarding resource assessment, modeling and statistical analyses.

Abstract. The Los Humeros geothermal system is steam dominated and currently under exploration with 65 wells (23 producing). Having temperatures above 380 ∘C, the system is characterized as a super hot geothermal system (SHGS). The development of such systems is still challenging due to the high temperatures and aggressive reservoir fluids which lead to corrosion and scaling problems. The geothermal system in Acoculco (Puebla, Mexico; so far only explored via two exploration wells) is characterized by temperatures of approximately 300 ∘C at a depth of about 2 km. In both wells no geothermal fluids were found, even though a well-developed fracture network exists. Therefore, it is planned to develop an enhanced geothermal system (EGS). For better reservoir understanding and prospective modeling, extensive geological, geochemical, geophysical and technical investigations are performed within the scope of the GEMex project. Outcrop analogue studies have been carried out in order to identify the main fracture pattern, geometry and distribution of geological units in the area and to characterize all key units from the basement to the cap rock regarding petro- and thermo-physical rock properties and mineralogy. Ongoing investigations aim to identify geological and structural heterogeneities on different scales to enable a more reliable prediction of reservoir properties. Beside geological investigations, physical properties of the reservoir fluids are determined to improve the understanding of the hydrochemical processes in the reservoir and the fluid-rock interactions, which affect the reservoir rock properties.

Against this background, a large-sample triaxial setup was developed at RWTH Aachen University under a project funded by BMWi (Federal Ministry of Economics and Technology). The planning and construction of the testing facility, including the analytical and numerical validation of the final experimental procedure, is thoroughly discussed in 2,3. The setup allows for experiments on rocks with size 30 cm × 30 cm × 45 cm under defined and controlled boundary conditions, yielding reproducible and repeatable datasets, well-suited for code benchmarking. In this paper we present hydraulic stimulation datasets from experiments where the borehole axis was parallel to the minimum horizontal stress direction and therefore the plane of crack propagation was parallel to the maximum horizontal stress direction. Investigated rock samples represent the reservoir rock types of a potential EGS site in Mexico. The datasets represent hydraulic stimulation responses in quasi-homogenous and extremely heterogeneous crystalline rock types, such as, a very fine-grained granite and a coarse-grained marble, respectively. Additionally, a dense network of 32 acoustic sensors, comprising of 28 GmuG standard ultrasonic sensors (http://www.gmugmbh.de/prod.html) and 4 Glaser amplitude-calibrated sensors 4,5 , was used to track the fracture propagation in real time. The validation dataset includes fluid injection rate, pressure in the injection interval, confining stresses, mechanical and petrophysical properties of the rock specimens, properties of the injection fluid, mechanical details of the experimental setup , and acoustic emission data. Additionally, we provide a Python-based code which can be used for processing the seismic data and visualizing of the experimental results. Methods The main components of the experimental setup , sample requirements and experimental protocol are described separately in the following subsections .

In this work, we aim to verify the predictions of the numerical simulators, which are used for designing field-scale hydraulic stimulation experiments. Although a strong theoretical understanding of this process has been gained over the past few decades, numerical predictions of fracture propagation in low-permeability rocks still remains a challenge. Against this background, we performed controlled laboratory-scale hydraulic fracturing experiments in granite samples, which not only provides high-quality experimental data but also a well-characterized experimental set-up. Using the experimental pressure responses and the final fracture sizes as benchmark, we compared the numerical predictions of two coupled hydraulic fracturing simulators—CSMP and GEOS. Both the simulators reproduced the experimental pressure behavior by implementing the physics of Linear Elastic Fracture Mechanics (LEFM) and lubrication theory within a reasonable degree of accuracy. The simulation results indicate that even in the very low-porosity (1–2 %) and low-permeability ($${10}^{-18}\ {\mathrm{m}}^{2}- {10}^{-19}\ {\mathrm{m}}^{2}$$ 10 - 18 m 2 - 10 - 19 m 2 ) crystalline rocks, which are usually the target of EGS, fluid-loss into the matrix and unsaturated flow impacts the formation breakdown pressure and the post-breakdown pressure trends. Therefore, underestimation of such parameters in numerical modeling can lead to significant underestimation of breakdown pressure. The simulation results also indicate the importance of implementing wellbore solvers for considering the effect of system compressibility and pressure drop due to friction in the injection line. The varying injection rate as a result of decompression at the instant of fracture initiation affects the fracture size, while the entry friction at the connection between the well and the initial notch may cause an increase in the measured breakdown pressure.

No abstract