Time filter

Source Type

Mottaghy D.,Geophysica Beratungsgesellschaft GmbH | Pechnig R.,Geophysica Beratungsgesellschaft GmbH | Vogt C.,RWTH Aachen
Geothermics | Year: 2011

The proposed Den Haag Zuidwest district heating system of the city of The Hague consists of a deep doublet in a Jurassic sandstone layer that is designed for a production temperature of 75°C and a reinjection temperature of 40°C at a flow rate of 150m 3h -1. The prediction of reservoir temperature and production behavior is crucial for success of the proposed geothermal doublet. This work presents the results of a study of the important geothermal and geohydrological issues for the doublet design. In the first phase of the study, the influences of the three-dimensional (3D) structures of anticlines and synclines on the temperature field were examined. A comprehensive petrophysical investigation was performed to build a large scale 3D-model of the reservoir. Several bottomhole temperatures (BHTs), as well as petrophysical logs were used to calibrate the model using thermal conductivity measurements on 50 samples from boreholes in different lithological units in the study area. Profiles and cross sections extracted from the calculated temperature field were used to study the temperature in the surrounding areas of the planned doublet. In the second phase of the project, a detailed 3D numerical reservoir model was set up, with the aim of predicting the evolution of the producer and injector temperatures, and the extent of the cooled area around the injector. The temperature model from the first phase provided the boundary conditions for the reservoir model. Hydraulic parameters for the target horizons, such as porosity and permeability, were taken from data available from the nearby exploration wells. The simulation results are encouraging as no significant thermal breakthrough is predicted. For the originally planned location of the producer, the extracted water temperature is predicted to be around 79°C, with an almost negligible cooling in the first 50 years of production. When the producer is located shallower parts of the reservoir, the yield water temperatures is lower, starting at ≈76°C and decreasing to ≈74°C after 50 years of operation. This comparatively larger decrease in temperature with time is caused by the structural feature of the reservoir, namely a higher dip causes the cooler water to easily move downward. In view of the poor reservoir data, the reservoir simulation model is constructed to allow iterative updates using data assimilation during planned drilling, testing, and production phases. Measurements during an 8h pumping test carried out in late 2010 suggest that a flow rate of 150m 3h -1 is achievable. Fluid temperatures of 76.5°C were measured, which is very close to the predicted value. © 2011 Elsevier Ltd.

Mottaghy D.,Geophysica Beratungsgesellschaft GmbH | Dijkshoorn L.,RWTH Aachen
Renewable Energy | Year: 2012

We present an effective finite difference formulation for implementing and modeling multiple borehole heat exchangers (BHE) in the general 3-D coupled heat and flow transport code SHEMAT. The BHE with arbitrary length can be either coaxial or double U-shaped. It is particularly suitable for modeling deep BHEs which contain varying pipe diameters and materials. Usually, in numerical simulations, a fine discretization of the BHE assemblage is required, due to the large geometric aspect ratios involved. This yields large models and long simulation times. The approach avoids this problem by considering heat transport between fluid and the soil through pipes and grout via thermal resistances. Therefore, the simulation time can be significantly reduced. The coupling with SHEMAT is realized by introducing an effective heat generation. Due to this connection, it is possible to consider heterogeneous geological models, as well as the influence of groundwater flow. This is particularly interesting when studying the long term behavior of a single BHE or a BHE field. Heating and cooling loads can enter the model with an arbitrary interval, e.g. from hourly to monthly values. When dealing with large BHE fields, computing times can be further significantly reduced by focusing on the temperature field around the BHEs, without explicitly modeling inlet and outlet temperatures. This allows to determine the possible migration of cold and warm plumes due to groundwater flow, which is of particular importance in urban areas with a high BHE installation density.The model is validated against the existing BHE modeling codes EWS and EED. A comparison with monitoring data from a deep BHE in Switzerland shows a good agreement. Synthetic examples demonstrate the field of application of this model. © 2012 Elsevier Ltd.

Mottaghy D.,Geophysica Beratungsgesellschaft GmbH | Schwamborn G.,Alfred Wegener Institute for Polar and Marine Research | Rath V.,Complutense University of Madrid
Climate of the Past | Year: 2013

This study focuses on the temperature field observed in boreholes drilled as part of interdisciplinary scientific campaign targeting the El'gygytgyn Crater Lake in NE Russia. Temperature data are available from two sites: the lake borehole 5011-1 located near the center of the lake reaching 400 m depth, and the land borehole 5011-3 at the rim of the lake, with a depth of 140 m. Constraints on permafrost depth and past climate changes are derived from numerical simulation of the thermal regime associated with the lake-related talik structure. The thermal properties of the subsurface needed for these simulations are based on laboratory measurements of representative cores from the quaternary sediments and the underlying impact-affected rock, complemented by further information from geophysical logs and data from published literature.

The temperature observations in the lake borehole 5011-1 are dominated by thermal perturbations related to the drilling process, and thus only give reliable values for the lowermost value in the borehole. Undisturbed temperature data recorded over more than two years are available in the 140 m deep land-based borehole 5011-3. The analysis of these observations allows determination of not only the recent mean annual ground surface temperature, but also the ground surface temperature history, though with large uncertainties. Although the depth of this borehole is by far too insufficient for a complete reconstruction of past temperatures back to the Last Glacial Maximum, it still affects the thermal regime, and thus permafrost depth. This effect is constrained by numerical modeling: assuming that the lake borehole observations are hardly influenced by the past changes in surface air temperature, an estimate of steady-state conditions is possible, leading to a meaningful value of 14 ± 5 K for the post-glacial warming. The strong curvature of the temperature data in shallower depths around 60 m can be explained by a comparatively large amplitude of the Little Ice Age (up to 4 K), with low temperatures prevailing far into the 20th century. Other mechanisms, like varying porosity, may also have an influence on the temperature profile, however, our modeling studies imply a major contribution from recent climate changes. © Author(s) 2013.

Kurten S.,RWTH Aachen | Mottaghy D.,Geophysica Beratungsgesellschaft GmbH | Ziegler M.,RWTH Aachen
Acta Geotechnica | Year: 2015

In this paper, a thermal resistance model for an energy wall using the example of thermo-active seal panels is presented. In the developed model, the resistances of the pipes as well as the resistance of the structure itself are considered. The resistance model is transferred to a 2D finite difference model, which itself is implemented into the general 3D subsurface heat and flow transport code SHEMAT-Suite. This coupling of a semi-analytical model with a numerical code avoids a complete discretisation of the model domain and thus enables fast computing times. This new approach has been verified by pure finite element simulations and by laboratory tests. © 2014, Springer-Verlag Berlin Heidelberg.

Chen T.,RWTH Aachen | Clauser C.,RWTH Aachen | Marquart G.,RWTH Aachen | Willbrand K.,RWTH Aachen | And 2 more authors.
Advances in Water Resources | Year: 2015

We present a method to determine equivalent permeability of fractured porous media. Inspired by the previous flow-based upscaling methods, we use a multi-boundary integration approach to compute flow rates within fractures. We apply a recently developed multi-point flux approximation Finite Volume method for discrete fracture model simulation. The method is verified by upscaling an arbitrarily oriented fracture which is crossing a Cartesian grid. We demonstrate the method by applying it to a long fracture, a fracture network and the fracture network with different matrix permeabilities. The equivalent permeability tensors of a long fracture crossing Cartesian grids are symmetric, and have identical values. The application to the fracture network case with increasing matrix permeabilities shows that the matrix permeability influences more the diagonal terms of the equivalent permeability tensor than the off-diagonal terms, but the off-diagonal terms remain important to correctly assess the flow field. © 2015 Elsevier Ltd.

Discover hidden collaborations