Characterization of fractured porous media using multivariate statistics and hydraulic travel time tomography

Fractured hard rocks are important in engineering, geotechnical, and hydrogeological practice because they act as hydraulic conductors providing pathways for fluid flow or barriers that prevent flow across them. In this thesis fractured porous media such as sandstones, which show a significant permeability of the matrix, are treated. In comparison to low permeable hard rocks such as granite and gneiss, fractured porous media show significant storage, flow, and transport e.g. of contaminants within the matrix. This has to be taken into account by the evaluation of flow and transport experiments. A fundamental step in understanding and predicting the behavior of fractured porous aquifers involves the identification of the spatial distribution of hydraulically significant features. Therefore, high resolution flow and transport data sets are required. Two different approaches are applied in order to evaluate high resolution flow and transport experiments arranged in a tomographical array. First, a multivariate statistical approach for the analysis of meso-scale flow and transport experiments is proposed. Within the statistical approach flow and transport measurements are classified into quasihomogeneous groups in order to determine properties and parameter zonation, which control flow and transport. The methodology was applied to a data set of gas flow and tracer experiments conducted to a gas saturated, fissured sandstone block of 60 × 60 × 60 cm3 in size. The agreement of the classification results with a geological surface mapping and numerical simulations show that it is possible to characterize and simplify complex systems without losing essential information about the investigated system. In the second part a travel time based hydraulic tomographic approach is proposed. The approach provides the inversion of travel times of hydraulic or pneumatic tests conducted in a tomographic array. The inversion is based on the relation between the peak time of a recorded transient pressure curve and the diffusivity of the investigated system. The development of a transformation factor enables the inversion of further travel times besides the peak time of a transient curve. As the early travel times of the curve are mainly related to preferential flow features while the inversion based on late travel times reflect an integral behavior, it can be assumed that the different inversion results reflect the properties of the overall system. The method has been applied to data from a set of pneumatic short term tests conducted to a gas saturated fractured sandstone cylinder with a height of 34 cm and a diameter of 31 cm. The three-dimensional reconstructions of the diffusivity distribution are found to be highly reliable and robust. In particular, the mapped fracture of the sandstone cylinder coincides with the reconstructed diffusivity distribution. In order to get a better insight of the hydraulic travel time approach, synthetic data sets are used. The synthetic data sets comprise simulated slug tests conducted in a tomographical array. The main focus of the synthetic study is to appraise the influence of the various travel times on the inversion results in dependence of the permeability distribution of the investigated geological medium. The inversion results are showing a strong dependence on the used travel times. The comparison of the different reconstructions based on different travel times with the forward model allows the verification and appraisal of the quality of the inversion. Hence, it is possible to suggest which travel time is most suited to reconstruct a certain subsurface structure.Geklüftete Festgesteine sind wichtig im Rahmen der ingenieurtechnischen, geotechnischen und hydrogeologischen Praxis, da sie entweder als hydraulische Leiter agieren und somit bevorzugte Fließpfade für Fluide darstellen oder sie die Eigenschaften von Geringleiter aufweisen und somit hydraulische Barrieren darstellen. In dieser Arbeit werden geklüftet poröse Festgesteine untersucht, die eine signifikante Matrixpermeabilität aufweisen. Im Vergleich zu gering durchlässigen Festgesteinen wie Granit oder Gneis weisen geklüftet poröse Medien bedeutende Speicher-, Strömung- und Transportkapazitäten auf, was sich z.B. auf den Transport von Schadstoffen innerhalb der Matrix auswirkt. Dies muss bei der Auswertung von Strömungs- und Transportexperimenten berücksichtigt werden. Ein grundlegender Schritt zum Verständnis und zur Vorhersage des Verhaltens von geklüftet porösen Aquiferen ist die Identifizierung der räumlichen Verteilung hydraulisch signifikanter Merkmale. Dafür werden hochauflösende Strömungs- und Transportdatensätze benötigt. Zur Auswertung der hochauflösenden Strömungs- und Transportexperimente, die unter Verwendung einer tomographischen Messanordnung aufgezeichnet worden sind, werden zwei unterschiedliche Ansätze verfolgt. Zuerst wird ein multivariater statistischer Ansatz für die Analyse von Strömungs- und Transportexperimenten vorgestellt. In diesem statistischen Ansatz werden Strömungs- und Transportmessungen in quasi homogene Gruppen eingeteilt, um die Eigenschaften und Parameterzonierung, die Strömung und Transport kontrollieren, zu bestimmen. Mit dieser Methode sind Strömungs- und Transportexperimente ausgewertet worden, die an einem gasgesättigten, geklüfteten Sandsteinblock mit einer Kantenlänge von 60 cm durchgeführt wurden. Die Klassifikationsergebnisse stimmen mit der geologischen Oberflächenkartierung und den numerischen Simulationen überein. Dies zeigt, dass es mit Hilfe des verwendeten Ansatzes möglich ist, komplexe Systeme zu charakterisieren und zu vereinfachen, ohne dass wichtige Systeminformationen verloren gehen. Im zweiten Teil wird ein auf Laufzeiten basierender Ansatz zur hydraulischen Tomographie vorgestellt. Dieser Ansatz umfasst die Inversion von Laufzeiten hydraulischer oder pneumatischer Tests, die unter Verwendung einer tomographischen Messanordnung aufgezeichnet worden sind. Die Inversion beruht dabei auf dem Zusammenhang zwischen der Laufzeit der maximalen Signaländerung eines aufgezeichneten instationären Signals und der Diffusivität des untersuchten Systems. Die Entwicklung eines Transformationsfaktors ermöglicht die Inversion von weiteren Laufzeiten, neben der Laufzeit der maximalen Änderung eines instationären Signals. Da frühe Laufzeiten einer aufgezeichneten Druckkurve überwiegend bevorzugte Flieswege und spätere Laufzeiten das integrale Verhalten des untersuchten Systems widerspiegeln, kann davon ausgegangen werden, dass die unterschiedlichen Inversionsergebnis das Verhalten des Gesamtsystems repräsentieren. Mit diese Methode wurden Daten von pneumatischen Kurzzeittests ausgewertet. Diese Daten sind an einem gasgesättigten, geklüfteten Sandsteinzylinder mit einer Höhe von 34 cm und einem Durchmesser von 31 cm aufgezeichnet worden. Die Ergebnisse der dreidimensionalen Diffusivitätsrekonstruktionen haben sich als zuverlässig und stabil erwiesen. Insbesondere zeigt die dreidimensionale Rekonstruktion der Diffusivitätsverteilung eine hohe Übereinstimmung mit der Lage der Kluft innerhalb des Zylinders. Um neue Erkenntnisse auf dem Gebiet der hydraulische Laufzeittomographie zu gewinnen, sind synthetische Datensätze verwendet worden. Die synthetischen Datensätze umfassen Slug-Tests, die unter Verwendung einer tomographischen Messanordnung modelliert worden sind. Ziel der synthetischen Fallstudie ist, den Einfluss unterschiedlicher Laufzeiten auf das Inversionsergebnis, in Abhängigkeit von der Permeabilitätsverteilung des untersuchten geologischen Mediums, zu bewerten. Die Inversionsergebnisse zeigen eine deutliche Abhängigkeit von den verwendeten Laufzeiten. Durch einen Vergleich zwischen den Diffusivitätsrekonstruktionen, die auf unterschiedlichen Laufzeiten basieren, und dem verwendeten Vorwärtsmodell ist es möglich, die Qualität der Inversionsergebnisse zu bewerten und zu überprüfen. Somit kann eine Inversionsstrategie vorgeschlagen werden, die für den jeweiligen Untergrund am besten geeignet ist

Bestimmung von hydraulischen Parametern in Lockergesteinen: Ein Vergleich unterschiedlicher Feldmethoden

Zusammenfassung: In dieser Feldstudie werden die laufzeitbasierte tomographische Inversion von Daten aus Kurzzeitpumpversuchen mit der analytischen Auswertung verglichen und die ermittelten hydraulischen Parameter hinsichtlich ihrer räumlichen Auflösung diskutiert und bewertet. Als Datenbasis dienen Messergebnisse aus Kurzzeitpumpversuchen, die in einer tomographischen Messanordnung in einem zwei Meter mächtigen, gut charakterisierten Sand- und Kiesgrundwasserleiter unter Verwendung eines 2"-Brunnens und eines Multikammerbrunnens, beide mit Direct-Push-Technik installiert, durchgeführt wurden. Die analytische Auswertung der Kurzzeitpumpversuche hat gezeigt, dass es nicht möglich ist, Bereiche mit unterschiedlichen hydraulischen Eigenschaften voneinander abzugrenzen. Entsprechend einem Vergleich mit den Ergebnissen von Multilevel-Slug-Tests werden die ermittelten hydraulischen Parameter, trotz einer geringen Pumpdauer von 200 Sekunden und hydraulisch isolierten Pump- und Beobachtungsintervallen, von einem hydraulisch höher durchlässigen Bereich am unteren Rand des Grundwasserleiters dominiert. Die laufzeitbasierte tomographische Inversion ermöglicht hingegen, vertikale und laterale Änderungen der Diffusivitätsverteilung zwischen Pump- und Beobachtungsbrunnen hochaufgelöst zu rekonstruiere

An extended trajectory-mechanics approach for calculating the path of a pressure transient: Travel-time tomography

The application of a technique from quantum dynamics to the governing equation for hydraulic head leads to a trajectory-based solution that is valid for a general porous medium. The semi-analytic expressions for propagation path and velocity of a change in hydraulic head form the basis of a travel-time tomographic imaging algorithm. An application of the imaging algorithm to synthetic arrival times reveals that a cross-well inversion based upon the extended trajectories correctly reproduces the magnitude of a reference model, improving upon an existing asymptotic approach. An inversion of hydraulic head arrival times from cross-well slug tests at the Widen field site in northern Switzerland captures a general decrease in permeability with depth, which is in agreement with previous studies, but also indicates the presence of a high-permeability feature in the upper portion of the cross-well plane.ISSN:1027-5606ISSN:1607-793

Hydrogeological characterization of the Heletz Sands reservoir, Heletz (Israel) as a preliminary step towards CO2 injection experiments

Bensabat, Jacob et al.Ponencia presentada en la European Geosciences Union General Assembly, celebrada en Viena del 7 al 12 de abril de 2013.One the major components of the EU-FP7 funded MUSTANG project is to conduct a highly controlled se- ries of CO 2 injection experiments, aimed at determining field values of key CO 2 trapping mechanisms such as dissolution and residual trapping and to establish a comprehensive and consistent dataset for model validation. Prior to injecting CO 2 there is a need to achieve a sufficient degree of hydrogeological characterization of the reservoir. In what follows we present a sequence of hydrologic tests to be conducted at Heletz and their expected contribution to the understanding relevant hydrogeology. These include: 1) Chemical characterization of the formation fluid; 2) Flowing Fluid Electrical Conductivity log, aimed at determining the vertical variability of the reservoir permeability in the near well vicinity; 3) Water pulse and pumping tests, aimed at determining the reservoir scale hydraulic properties; 4) Thermal test, aimed at determining the value of the heat transfer coefficient from the reservoir to the borehole fluid, which is responsible for the heating of injected fluid in the borehole; 5) two-well injection and pumping of water and tracers test, in order to determine the impact of heterogeneity on the hydraulic parameters and to identify preferential flow paths in the reservoir. This paper presents the design and planning of the experiments, the results obtained in field and a preliminary interpretation.Peer reviewe

Bestimmung von hydraulischen Parametern in Lockergesteinen: Ein Vergleich unterschiedlicher Feldmethoden

ISSN:1432-1165ISSN:1430-483

Widen crosswell slug test data

<p>Cross well slug test data gathers by Ralf Brauchler at the Widen test site in northern Switzerland.</p> <p>The program WidenReduction.f may be used to read the data.</p

Time-lapse pressure tomography for characterizing CO2 plume evolution in a deep saline aquifer

A time-lapse pressure tomography inversion approach is applied to characterize the CO2 plume development in a virtual deep saline aquifer. Deep CO2 injection leads to flow properties of the mixed-phase, which vary depending on the CO2 saturation. Analogous to the crossed ray paths of a seismic tomographic experiment, pressure tomography creates streamline patterns by injecting brine prior to CO2 injection or by injecting small amounts of CO2 into the two-phase (brine and CO2) system at different depths. In a first step, the introduced pressure responses at observation locations are utilized for a computationally rapid and efficient eikonal equation based inversion to reconstruct the heterogeneity of the subsurface with diffusivity (D) tomograms. Information about the plume shape can be derived by comparing D-tomograms of the aquifer at different times. In a second step, the aquifer is subdivided into two zones of constant values of hydraulic conductivity (K) and specific storage (Ss) through a clustering approach. For the CO2 plume, mixed-phase K and Ss values are estimated by minimizing the difference between calculated and “true” pressure responses using a single-phase flow simulator to reduce the computing complexity. Finally, the estimated flow property is converted to gas saturation by a single-phase proxy, which represents an integrated value of the plume. This novel approach is tested first with a doublet well configuration, and it reveals a great potential of pressure tomography based concepts for characterizing and monitoring deep aquifers, as well as the evolution of a CO2 plume. Still, field-testing will be required for better assessing the applicability of this approach