Search Results
Predictability of properties of a fractured geothermal reservoir: the opportunities and limitations of an outcrop analogue study
Minimizing exploration risk in deep geothermics is of great economic importance. Especially, knowledge about temperature and permeability of the reservoir is essential. We test the potential of an outcrop analogue study to minimize uncertainties in prediction of the rock properties of a fractured reservoir in the Upper Rhine Graben. Our results show that although mineralogical composition, clay content, grain size, and fabric type are basically comparable, porosity and quartz cementation are not. Young’s modulus, as observed in the outcrop closest to the reservoir is about twice as high (~ 64 GPa) as observed in the reservoir (~ 34 GPa). Most importantly, however, the parameters that describe the fracture system, which are essential to predict reservoir permeability, differ significantly. While the outcrops are dominated by perpendicular fracture sets (striking NE–SW and NW–SE), two different conjugate fracture sets (striking NW–SE and N–S) occur in the reservoir. Fracture apertures, as reported from the FMI, are one order of magnitude wider than in the outcrop. We conclude that our outcrop analogue study fails to predict important properties of the reservoir (such as permeability and porosity). This must be in part because of the tectonically complex setting of the reservoir. We propose that analogue studies are important, but they must be treated with care when attempting to predict the controlling parameters of a fractured reservoir.
Post-IR IRSL dating of K-feldspar from last interglacial marine terrace deposits on the Kamikita coastal plain, northeastern Japan
To establish a suitable luminescence dating protocol for marine terrace deposits in Japan, we tested the applicability of K-feldspar post-infrared (IR) infrared stimulated luminescence (IRSL) (pIRIR) dating using a marine isotope stage (MIS) 5e terrace deposit from the Kamikita coastal plain (NE Japan), where independent age control from a tephra is available. One of the most commonly used pIRIR signals, measured at 290°C with the first IR stimulation temperature at 50°C (pIRIR50/290), faded with a mean g2days value of 1.94 ± 0.19%/decade. In contrast, the pIRIR signal with a higher first IR stimulation temperature of 200°C (pIRIR200/290) had a much lower fading rate (g2days = 0.16 ± 0.49%/decade). The average fading-uncorrected and -corrected pIRIR200/290 ages of MIS 5e subtidal sediments obtained from two sampling sites were 126 ± 3 ka and 132 ± 2 ka, which is in good agreement with the independent age control. We conclude that is it is now possible to use pIRIR protocol to estimate the ages of not only marine terraces formed during MIS 5 substages (5a, 5c) but also of older marine terraces, for which age evidence is limited.
pyGIMLi: An open-source library for modelling and inversion in geophysics
Many tasks in applied geosciences cannot be solved by single measurements, but require the integration of geophysical, geotechnical and hydrological methods. Numerical simulation techniques are essential both for planning and interpretation, as well as for the process understanding of modern geophysical methods. These trends encourage open, simple, and modern software architectures aiming at a uniform interface for interdisciplinary and flexible modelling and inversion approaches. We present pyGIMLi (Python Library for Inversion and Modelling in Geophysics), an open-source framework that provides tools for modelling and inversion of various geophysical but also hydrological methods. The modelling component supplies discretization management and the numerical basis for finite-element and finite-volume solvers in 1D, 2D and 3D on arbitrarily structured meshes. The generalized inversion framework solves the minimization problem with a Gauss-Newton algorithm for any physical forward operator and provides opportunities for uncertainty and resolution analyses. More general requirements, such as flexible regularization strategies, time-lapse processing and different sorts of coupling individual methods are provided independently of the actual methods used. The usage of pyGIMLi is first demonstrated by solving the steady-state heat equation, followed by a demonstration of more complex capabilities for the combination of different geophysical data sets. A fully coupled hydrogeophysical inversion of electrical resistivity tomography (ERT) data of a simulated tracer experiment is presented that allows to directly reconstruct the underlying hydraulic conductivity distribution of the aquifer. Another example demonstrates the improvement of jointly inverting ERT and ultrasonic data with respect to saturation by a new approach that incorporates petrophysical relations in the inversion. Potential applications of the presented framework are manifold and include time-lapse, constrained, joint, and coupled inversions of various geophysical and hydrological data sets.
3D architecture of cyclic-step and antidune deposits in glacigenic subaqueous fan and delta settings: Integrating outcrop and ground-penetrating radar data
Bedforms related to supercritical flows are increasingly recognised as important constituents of many depositional environments, but outcrop studies are commonly hampered by long bedform wavelengths and complex three-dimensional geometries. We combined outcrop-based facies analysis with ground-penetrating radar (GPR) surveys to analyse the 3D facies architecture of subaqueous ice-contact fan and glacifluvial delta deposits. The studied sedimentary systems were deposited at the margins of the Middle Pleistocene Scandinavian ice sheets in Northern Germany. Glacifluvial Gilbert-type deltas are characterised by steeply dipping foreset beds, comprising cyclic-step deposits, which alternate with antidune deposits. Deposits of cyclic steps consist of lenticular scours infilled by backset cross-stratified pebbly sand and gravel. The GPR sections show that the scour fills form trains along the delta foresets, which can locally be traced for up to 15 m. Perpendicular and oblique to palaeoflow direction, these deposits appear as troughs with concentric or low-angle cross-stratified infills. Downflow transitions from scour fills into sheet-like low-angle cross-stratified or sinusoidally stratified pebbly sand, deposited by antidunes, are common. Cyclic steps and antidunes were deposited by sustained and surge-type supercritical density flows, which were related to hyperpycnal flows, triggered by major meltwater discharge or slope-failure events. Subaqueous ice-contact fan deposits include deposits of progradational scour fills, isolated hydraulic jumps, antidunes and (humpback) dunes. The gravel-rich fan succession consists of vertical stacks of laterally amalgamated pseudo-sheets, indicating deposition by pulses of waning supercritical flows under high aggradation rates. The GPR sections reveal the large-scale architecture of the sand-rich fan succession, which is characterised by lobe elements with basal erosional surfaces associated with scours filled with backsets related to hydraulic jumps, passing upwards and downflow into deposits of antidunes and (humpback) dunes. The recurrent facies architecture of the lobe elements and their prograding and retrograding stacking pattern are interpreted as related to autogenic flow morphodynamics.
Geophysical investigation of a freshwater lens on the island of Langeoog, Germany – Insights from combined HEM, TEM and MRS data
A multi-method geophysical survey, including helicopter-borne electromagnetics (HEM), transient electromagnetics (TEM), and magnetic resonance sounding (MRS), was conducted to investigate a freshwater lens on the North Sea island of Langeoog, Germany. The HEM survey covers the entire island and gives an overview of the extent of three freshwater lenses that reach depths of up to 45 m. Ground-based TEM and MRS were conducted particularly on the managed western lens to verify the HEM results and to complement the lithological information from existing boreholes. The results of HEM and TEM are in good agreement. Salt- and freshwater-bearing sediments can, as expected, clearly be distinguished due to their individual resistivity ranges. In the resistivity data, a large transition zone between fresh- and saltwater with a thickness of up to 20 m is identified, the existence of which is verified by borehole logging and sampling. Regarding lithological characterisation of the subsurface, the MRS method provides more accurate and reliable results than HEM and TEM. Using a lithological index derived from MRS water content and relaxation time, thin aquitard structures as well as fine and coarse sand aquifers can be distinguished. Complementing the existing borehole data with the lithology information estimated from MRS, we generate a map showing the occurrence of aquitard structures, which significantly improves the hydrogeological model of the island. Moreover, we demonstrate that the estimates of groundwater conductivity in the sand aquifers from geophysical data are in agreement with the fluid conductivity measured in the boreholes.
Quartz OSL and K-feldspar post-IR IRSL dating of sand accumulation in the Lower Liao Plain (Liaoning, NE China)
The timing of the formation of extensively distributed sand dunes in the Bohai coastal area and its forcing factors are poorly understood. In this study, the chronology of a well-preserved sand dune located in Panjin Forest Park (PJ) in the Lower Liao Plain (LLP) is investigated using quartz optically stimulated luminescence (OSL) and K-feldspar postinfrared (IR) infrared stimulated luminescence (IRSL) (pIRIR) dating. For the pIRIR measurements, the combination of preheating at 180°C and pIRIR stimulation at 150°C (pIRIR150) is exploited. The quartz results show that the sand dune accumulated from c. 120 a (1890 AD) to c. 70 a (1940 AD) before present, and the underlying sandy soil sediments deposited from c. 5.0 ka to c. 0.13 ka as marsh sediment after the sea level highstand since the mid-Holocene. From the evidence in historical coastline records, the PJ sand dune is an inland sand dune and not a coastal sand dune. Based on further information of climate and temperature change after the Little Ice Age (LIA) and human activity in northeastern China, we conclude that the PJ sand dune accumulation was very likely impacted by the immigrants and land reclamation at the end of Qing dynasty. The fading corrected IR50 ages, the apparent and fading corrected pIRIR150 ages are consistent with quartz ages for two sandy soil samples but overestimate those for six sand samples. The overestimation of the feldspar ages is derived from the residual signal which has not been bleached before burial. The offset obtained from the difference between the quartz OSL and the feldspar pIRIR150 ages are ~20-160 a (predicted residual dose: ~0.08-0.60 Gy), whereas the measured residual dose after bleaching 4 h in a solar simulator yielded age overestimation of ~10-40 a (~0.05-0.16 Gy). The age discrepancy calculated from the predicted residual was larger than those obtained from the laboratory measured residuals. We conclude that the pIRIR150 of aeolian sediment is applicable for samples older than ~1000 years where the effect of the residual dose become negligible. © 2016 Y. Li et al.
The environmental and evolutionary history of Lake Ohrid (FYROM/Albania): interim results from the SCOPSCO deep drilling project
This study reviews and synthesises existing information generated within the SCOPSCO (Scientific Collaboration on Past Speciation Conditions in Lake Ohrid) deep drilling project. The four main aims of the project are to infer (i) the age and origin of Lake Ohrid (Former Yugoslav Republic of Macedonia/Republic of Albania), (ii) its regional seismotectonic history, (iii) volcanic activity and climate change in the central northern Mediterranean region, and (iv) the influence of major geological events on the evolution of its endemic species. The Ohrid basin formed by transtension during the Miocene, opened during the Pliocene and Pleistocene, and the lake established de novo in the still relatively narrow valley between 1.9 and 1.3 Ma. The lake history is recorded in a 584 m long sediment sequence, which was recovered within the framework of the International Continental Scientific Drilling Program (ICDP) from the central part (DEEP site) of the lake in spring 2013. To date, 54 tephra and cryptotephra horizons have been found in the upper 460 m of this sequence. Tephrochronology and tuning biogeochemical proxy data to orbital parameters revealed that the upper 247.8 m represent the last 637 kyr. The multi-proxy data set covering these 637 kyr indicates long-term variability. Some proxies show a change from generally cooler and wetter to drier and warmer glacial and interglacial periods around 300 ka. Short-term environmental change caused, for example, by tephra deposition or the climatic impact of millennial-scale Dansgaard-Oeschger and Heinrich events are superimposed on the long-term trends. Evolutionary studies on the extant fauna indicate that Lake Ohrid was not a refugial area for regional freshwater animals. This differs from the surrounding catchment, where the mountainous setting with relatively high water availability provided a refuge for temperate and montane trees during the relatively cold and dry glacial periods. Although Lake Ohrid experienced significant environmental change over the last 637 kyr, preliminary molecular data from extant microgastropod species do not indicate significant changes in diversification rate during this period. The reasons for this constant rate remain largely unknown, but a possible lack of environmentally induced extinction events in Lake Ohrid and/or the high resilience of the ecosystems may have played a role. © Author(s) 2017.
Structural analysis of S-wave seismics around an urban sinkhole: evidence of enhanced dissolution in a strike-slip fault zone
In November 2010, a large sinkhole opened up in the urban area of Schmalkalden, Germany. To determine the key factors which benefited the development of this collapse structure and therefore the dissolution, we carried out several shear-wave reflection-seismic profiles around the sinkhole. In the seismic sections we see evidence of the Mesozoic tectonic movement in the form of a NW–SE striking, dextral strike-slip fault, known as the Heßleser Fault, which faulted and fractured the subsurface below the town. The strike-slip faulting created a zone of small blocks ( < 100 m in size), around which steep-dipping normal faults, reverse faults and a dense fracture network serve as fluid pathways for the artesian-confined groundwater. The faults also acted as barriers for horizontal groundwater flow perpendicular to the fault planes. Instead groundwater flows along the faults which serve as conduits and forms cavities in the Permian deposits below ca. 60 m depth. Mass movements and the resulting cavities lead to the formation of sinkholes and dissolution-induced depressions. Since the processes are still ongoing, the occurrence of a new sinkhole cannot be ruled out. This case study demonstrates how S-wave seismics can characterize a sinkhole and, together with geological information, can be used to study the processes that result in sinkhole formation, such as a near-surface fault zone located in soluble rocks. The more complex the fault geometry and interaction between faults, the more prone an area is to sinkhole occurrence.
Impacts of a capillary barrier on infiltration and subsurface stormflow in layered slope deposits monitored with 3-D ERT and hydrometric measurements
Identifying principles of water movement in the shallow subsurface is crucial for adequate process-based hydrological models. Hillslopes are the essential interface for water movement in catchments. The shallow subsurface on slopes typically consists of different layers with varying characteristics. The aim of this study was to draw conclusions about the infiltration behaviour, to identify water flow pathways and derive some general interpretations for the validity of the water movement on a hillslope with periglacial slope deposits (cover beds), where the layers differ in their sedimentological and hydrological properties. Especially the described varying influence of the basal layer (LB) as an impeding layer on the one hand and as a remarkable pathway for rapid subsurface stormflow on the other. We used a time lapse 3-D electrical resistivity tomography (ERT) approach combined with punctual hydrometric data to trace the spreading and the progression of an irrigation plume in layered slope deposits during two irrigation experiments. This multi-technical approach enables us to connect the high spatial resolution of the 3-D ERT with the high temporal resolution of the hydrometric devices. Infiltration through the uppermost layer was dominated by preferential flow, whereas the water flow in the deeper layers was mainly matrix flow. Subsurface stormflow due to impeding characteristic of the underlying layer occurs in form of "organic layer interflow" and at the interface to the first basal layer (LB1). However, the main driving factor for subsurface stormflow is the formation of a capillary barrier at the interface to the second basal layer (LB2). The capillary barrier prevents water from entering the deeper layer under unsaturated conditions and diverts the seepage water according to the slope inclination. With higher saturation, the capillary barrier breaks down and water reaches the highly conductive deeper layer. This highlights the importance of the capillary barrier effect for the prevention or activation of different flow pathways under variable hydrological conditions.
Drilling into an active mofette: pilot-hole study of the impact of CO2-rich mantle-derived fluids on the geo–bio interaction in the western Eger Rift (Czech Republic)
Microbial life in the continental "deep biosphere" is closely linked to geodynamic processes, yet this interaction is poorly studied. The Cheb Basin in the western Eger Rift (Czech Republic) is an ideal place for such a study because it displays almost permanent seismic activity along active faults with earthquake swarms up to ML 4.5 and intense degassing of mantle-derived CO2 in conduits that show up at the surface in form of mofettes. We hypothesize that microbial life is significantly accelerated in active fault zones and in CO2 conduits, due to increased fluid and substrate flow. To test this hypothesis, pilot hole HJB-1 was drilled in spring 2016 at the major mofette of the Hartoušov mofette field, after extensive pre-drill surveys to optimize the well location. After drilling through a thin caprock-like structure at 78.5 m, a CO2 blowout occurred indicating a CO2 reservoir in the underlying sandy clay. A pumping test revealed the presence of mineral water dominated by Na+, Ca2+, HCO3−, SO42− (Na-Ca-HCO3-SO4 type) having a temperature of 18.6 °C and a conductivity of 6760 µS cm−1. The high content of sulfate (1470 mg L−1) is typical of Carlsbad Spa mineral waters. The hole penetrated about 90 m of Cenozoic sediments and reached a final depth of 108.50 m in Palaeozoic schists. Core recovery was about 85 %. The cored sediments are mudstones with minor carbonates, sandstones and lignite coals that were deposited in a lacustrine environment. Deformation structures and alteration features are abundant in the core. Ongoing studies will show if they result from the flow of CO2-rich fluids or not.