Geomechanical model for fracture deformation under hydraulic, mechanical and thermal loads     

Chris McDermott  Olaf Kolditz 《Hydrogeology Journal》2006,14(4):485-498

Hydraulic flow and transport (heat and solute) within crystalline rocks is dominated by the fracture systems found within them. In situ stress conditions have a significant impact on the hydraulic, mechanical and thermal coupled processes, and quantification of these processes provides a key to understanding the often transient time-dependent behaviour of crystalline rocks. In this paper, a geomechanical model is presented which describes fracture closure as a function of effective stress and the changes in parameters such as storage, permeability, porosity and aperture. Allowing the fracture closure to be defined by the change in normal effective stress provides a link to the numerical consideration of parametrical changes due to rock stress alterations caused for example by changes in fracture fluid pressure, stress release, tectonic stress, thermal stress, orientation of the natural fracture in the pervasive stress system and local changes in a rock mass due to stress alteration. The model uses geometrical considerations based on a fractal distribution of apertures on the fracture surface, and applies well-established analytical elastic deformation solutions to calculate the deformation response to changes in effective stress. Analysis of the fractal generation method allows a standard normal distribution of fracture apertures to be predicted for all common fractal dimensions relating to a 2D surface. Changes in the fracture aperture are related to hydraulic functions such as permeability, storage and porosity of the fracture. The geomechanical model is experimentally validated against laboratory scale experimental data gained from the closure of a fractured sample recovered at a depth of 3,800 m from the KTB pilot borehole. Parameters for matching the experimental data were established externally, the only fitting parameters applied were the minimum and maximum contact area between the surfaces and the number of allowable contacts. The model provides an insight into the key processes determining the closure of a fracture, and can act as a material input function for numerical models linking the effects of changes in the stress field, hydraulic or thermal conditions, to the flow and transport parameters of a fractured system.

---

Résumé L’écoulement et le transport (chaleur et soluté) dans les roches cristallines sont dominés par les systèmes de fracture. Les conditions de stress in-situ ont un impact significatif sur l’hydraulique, les processus couplés de mécanique et thermique et la quantification de ces processus apportent une clé pour comprendre le comportement transitoire des roches cristallines. Dans cet article un modèle géomécanique est présenté, modèle qui décrit la fermeture des fractures comme une fonction de la contrainte effective et des changements de paramètres tels le coefficient d’emmagasinement, la perméabilité, la porosité et l’ouverture. En s’accordant que la fermeture des fractures est définit par les changements de la contrainte effective normale, on apporte le lien avec la considération numérique des changements paramétriques dus aux altérations de la contrainte des roches, causés par exemple par des variations de la pression des fluides dans les fractures, du dégagement de la contrainte, des contraintes tectoniques et thermiques, des orientations des fractures naturelles dans le système de contraintes pénétrantes, et des changements locaux dans un massif de roches dus à l’altération des contraintes. Le modèle utilise des considérations géométriques basées sur une distribution fractale des ouvertures à la surface des fractures, et permet d’établir des solutions analytiques de la déformation élastique pour calculer la réponse de la déformation à la contrainte effective. L’analyse de la méthode par génération fractale permet de prédire une distribution normale standard de l’ouverture des fractures, pour toutes les dimensions fractales en relation avec les surfaces 2D. Les changements dans l’ouverture des fractures sont mis en relation avec les fonctions hydrauliques tels la perméabilité, l’emmagasinement et la porosité de la fracture. Le modèle géoméchanique est expérimentalement validé à l’échelle du laboratoire sur un échantillon fracturé récupéré à une profondeur de 3,800 mètres sur le puits du site pilote KTB. Les paramètres du calibrage des données expérimentales ont été établies extérieurement, les seuls paramètres utilisés étant les surfaces de contact minimum et maximum, et le nombre de contacts permis. Le modèle apporte une connaissance perspicace sur le processus clé déterminant la fermeture des fractures, et peut servir de fonction input dans les modèles numériques reliant les effets des variations de la contrainte du terrain, les conditions hydrauliques ou thermales, les paramètres de l’écoulement et du transport et les systèmes de fracture.

---

Resumen El flujo hidráulico y transporte (de calor y solutos) dentro de rocas cristalinas está dominado por los sistemas de fracturas que se encuentran en ellas. Las condiciones de esfuerzos in-situ tienen un impacto significativo en los procesos aparejados termales, mecánicos e hidráulicos y la cuantificación de estos procesos aporta una clave para entender el frecuente comportamiento transitorio dependiente de las rocas cristalinas. En este artículo se presenta un modelo geomecánico que describe el cierre de fracturas en función del esfuerzo efectivo y los cambios en parámetros tal como almacenamiento, permeabilidad, porosidad y apertura. El definir el cierre de fractura mediante el cambio en esfuerzo normal efectivo aporta un vínculo con la consideración numérica de cambios paramétricos ocasionados por alteraciones de esfuerzos en la roca causadas, por ejemplo, por cambios en presión de fluidos en fractura, liberación de esfuerzo, esfuerzo tectónico, esfuerzo termal, orientación de fracturas naturales en el sistema de esfuerzos penetrante, y cambios locales en una masa rocosa ocasionados por alteración de esfuerzos. El modelo utiliza consideraciones geométricas basadas en la distribución fractal de aperturas en la superficie de fractura y aplica soluciones analíticas bien establecidas de deformación elástica para calcular la respuesta de deformación a cambios en el esfuerzo efectivo. Los análisis del método de generación fractal permiten predecir una distribución normal standard para la distribución de aperturas de fracturas para todas las dimensiones fractales comunes que se relacionan con una superficie 2D. Los cambios en la apertura de fractura se relacionan con funciones hidráulicas tal como permeabilidad, almacenamiento y porosidad de la fractura. El modelo geomecánico se ha validado experimentalmente en contra de datos experimentales a escala de laboratorio obtenidos a partir del cierre de una muestra fracturada recuperada a una profundidad de 3,800 m en el pozo piloto KTB. Se establecieron externamente parámetros que se ajustan a los datos experimentales, con los parámetros de ajuste aplicados que fueron el área máxima y mínima de contacto entre las superficies y el número de contactos permisibles. El modelo arroja luz sobre los procesos clave que determinan el cierre de una fractura y puede actuar como un material de función de entrada para modelos numéricos que vinculan los efectos de cambios en el campo de esfuerzos, condiciones termales o hidráulicas, con los parámetros de flujo y transporte de un sistema fracturado.

  相似文献   

Saprolite Remnants as Indicators of Pre-Glacial Landform Genesis in Southeast Sweden   总被引：3，自引：0，他引：3  

Mats Olvmo  Karna Lidmar-Bergström  Kerstin Ericson  Johan M. Bonow 《Geografiska Annaler: Series A, Physical Geography》2005,87(3):447-460

Twenty‐six sites with remnants of gravelly saprolites (grus) have been located in southeast Sweden. Joint block hills (castle kopjes) and steep rock walls with weathered joints as well as rounded boulders are documented to have an origin in deep weathering and subsequent stripping of saprolites. The saprolite remnants and landforms result from the fragmentation of the re‐exposed sub‐Cambrian peneplain along fracture systems. Only shallow saprolites occur on the elevated intact parts of the sub‐Cambrian peneplain, while saprolites up to 20 m thick are encountered in areas where the sub‐Cambrian peneplain is fractured and dissected. Neogene uplift with reactivation of the weathering system is thought to be the main cause of saprolite formation. Deep weathering is thus judged to have been the major agent of landform formation in the study area, while glacial and glaciofluvial erosion has contributed mainly by stripping saprolites, detaching corestones, and plucking joint blocks along weathered joints.  相似文献   

Metal binding to dissolved organic matter and adsorption to ferrihydrite in shallow peat groundwaters: Application to diamond exploration in the James Bay Lowlands,Canada     

Jamil A. Sader  Keiko Hattori  Stewart Hamilton  Kerstin Brauneder 《Applied Geochemistry》2011,26(9-10):1649-1664

The speciation and solubility of kimberlite pathfinder metals (Ni, Nd, Ba and K) in shallow peat groundwaters is investigated over the Yankee, Zulu and Golf kimberlites in the Attawapiskat region, James Bay Lowlands, Canada. The purpose of this study is to examine the relationship between dissolved organic matter (DOM) complexation with kimberlite pathfinder metals and determine the spatial distribution of those metals in shallow peat groundwaters along sampling transects over subcropping kimberlites. Nickel, Nd, Ba and K complexation with DOM and the adsorption of these metals onto ferrihydrite were calculated using Visual MINTEQ 3.0 and the NICA-Donnan database. Calculations predict almost 100% of soluble Nd, Ni and Ba form complexes with DOM at sampling sites with little to no contribution from upwelling groundwater (i.e., dissolved organic C (DOC) concentrations = 40–132 mg/L, pH = 3.9–5.5, and log ionic strength ??3). In only the most ombrotrophic peat groundwater conditions does a majority fraction of K bind to DOM. By contrast, under conditions with large contributions from upwelling groundwaters (i.e., DOC concentrations ?40 mg/L, pH = 5.5–6.5, and log ionic strength = ?3 to ?2), as little as 10% of Nd and Ni, and 0% K and Ba are predicted to complex with DOM. The modeling calculations suggest the dominant control on metal–DOM complexation, particularly with respect to Ni and Nd, is competitive effects for DOM binding sites due to elevated ionic strength where there is evidence of strong groundwater upwelling. Visual MINTEQ modeling of metal adsorption on ferrihydrite surfaces predicts that under strong upwelling conditions, Ni and Nd are scavenged from solution due to increased ferrihydrite precipitation and decreased fractions of metals complexed with DOM. Analytical geochemical data are consistent with model predictions of metal adsorption on ferrihydrite. Total dissolved Ni and Nd concentrations at sites of strong upwelling are up to five times lower than waters with little to no upwelling and log ferrihydrite saturation indices (logSI_ferr) indicate precipitation (values up to 5) at sites of strong groundwater upwelling. Where the majority of Ni and Nd complex with DOM and ferrihydrite is highly under saturated (logSI_ferr = ?18 to ?5), the concentrations of total Ni and Nd are elevated compared to other sites along sampling transects. Metal complexation with DOM effectively inhibits metal scavenging from solution via adsorption and/or from forming secondary mineral precipitates. Also, because alkaline earth metals do not compete strongly with Ni and Nd for adsorption sites on ferrihydrite surfaces, but do compete strongly for insoluble organic sites, Ni and Nd are more likely to adsorb onto ferrihydrite.  相似文献   

Is there a superior conceptual groundwater model structure for baseflow simulation?          下载免费PDF全文

Michael Stoelzle  Markus Weiler  Kerstin Stahl  Andreas Morhard  Tobias Schuetz 《水文研究》2015,29(6):1301-1313

Previous work has shown that streamflow response during baseflow conditions is a function of storage, but also that this functional relationship varies among seasons and catchments. Traditionally, hydrological models incorporate conceptual groundwater models consisting of linear or non‐linear storage–outflow functions. Identification of the right model structure and model parameterization however is challenging. The aim of this paper is to systematically test different model structures in a set of catchments where different aquifer types govern baseflow generation processes. Nine different two‐parameter conceptual groundwater models are applied with multi‐objective calibration to transform two different groundwater recharge series derived from a soil‐atmosphere‐vegetation transfer model into baseflow separated from streamflow data. The relative performance differences of the model structures allow to systematically improve the understanding of baseflow generation processes and to identify most appropriate model structures for different aquifer types. We found more versatile and more aquifer‐specific optimal model structures and elucidate the role of interflow, flow paths, recharge regimes and partially contributing storages. Aquifer‐specific recommendations of storage models were found for fractured and karstic aquifers, whereas large storage capacities blur the identification of superior model structures for complex and porous aquifers. A model performance matrix is presented, which highlights the joint effects of different recharge inputs, calibration criteria, model structures and aquifer types. The matrix is a guidance to improve groundwater model structures towards their representation of the dominant baseflow generation processes of specific aquifer types. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

Analysis of the THM behaviour in a clay-based engineered barrier system (EBS): modelling of the HE-E experiment (Mont Terri URL)     

Xuerui Wang  Hua Shao  Jürgen Hesser  Olaf Kolditz 《Environmental Earth Sciences》2016,75(20):1350

  相似文献   

Comparative modelling of laboratory experiments for the hydro-mechanical behaviour of a compacted bentonite–sand mixture     

A. Millard  N. Mokni  J. D. Barnichon  K. E. Thatcher  A. E. Bond  A. Fraser-Harris  C. Mc Dermott  R. Blaheta  Z. Michalec  M. Hasal  T. S. Nguyen  O. Nasir  R. Fedors  H. Yi  O. Kolditz 《Environmental Earth Sciences》2016,75(19):1311

  相似文献   

Similarity-based approaches in hydrogeology: proposal of a new concept for data-scarce groundwater resource characterization and prediction     

Barthel  Roland  Haaf  Ezra  Giese  Markus  Nygren  Michelle  Heudorfer  Benedikt  Stahl  Kerstin 《Hydrogeology Journal》2021,29(5):1693-1709

Hydrogeology Journal - A new concept is proposed for describing, analysing and predicting the dynamic behaviour of groundwater resources based on classification and similarity. The concept makes...  相似文献   

Long-term effects of a toxic algal bloom on subtidal soft-sediment macroinvertebrate communities in Wellington Harbour, New Zealand   总被引：3，自引：1，他引：2  

Kerstin Krger  Jonathan P.A. Gardner  Ashley A. Rowden  Robert G. Wear 《Estuarine, Coastal and Shelf Science》2006,67(4):589-604

The long-term effects (>1 year) of a naturally occurring toxic plankton bloom (Karenia brevisulcata) on subtidal benthic macroinvertebrate communities were investigated in Wellington Harbour, a semi-enclosed temperate embayment in New Zealand. For 3 years communities were sampled at three different sites in the harbour. Analyses revealed that community recovery following the bloom was site-specific. Multivariate analyses indicated that at one site community composition was approaching recovery 3 years post-bloom. At the second site, a sequential recovery process was indicated, whereas at the third site the community composition oscillated from year to year, but did not show any signs of a sequential recovery process. The nature of the hydrodynamic regime was identified as a major factor influencing the observed recovery processes. Communities exposed to an active hydrodynamic regime were less affected by the bloom and differed little in their composition pre- and post-bloom, as they were naturally in a perpetual state of recovery as indicated by a dominance of r-selected species. The community at the hydrodynamically less active site was more affected by the bloom and exhibited temporal differences in composition consistent with successional models. Complete recovery to a pre-disturbance climax community dominated by K-selected species is likely to take 4–5 years, if not interrupted by other disturbances. Given the increased occurrence of harmful algal blooms worldwide, more monitoring and manipulative studies are needed to further evaluate the effects of such disturbances on macrobenthic communities.  相似文献   

Multi-angle backscatter classification and sub-bottom profiling for improved seafloor characterization     

Evangelos Alevizos  Mirjam Snellen  Dick Simons  Kerstin Siemes  Jens Greinert 《Marine Geophysical Researches》2018,39(1-2):289-306

This study applies three classification methods exploiting the angular dependence of acoustic seafloor backscatter along with high resolution sub-bottom profiling for seafloor sediment characterization in the Eckernförde Bay, Baltic Sea Germany. This area is well suited for acoustic backscatter studies due to its shallowness, its smooth bathymetry and the presence of a wide range of sediment types. Backscatter data were acquired using a Seabeam1180 (180 kHz) multibeam echosounder and sub-bottom profiler data were recorded using a SES-2000 parametric sonar transmitting 6 and 12 kHz. The high density of seafloor soundings allowed extracting backscatter layers for five beam angles over a large part of the surveyed area. A Bayesian probability method was employed for sediment classification based on the backscatter variability at a single incidence angle, whereas Maximum Likelihood Classification (MLC) and Principal Components Analysis (PCA) were applied to the multi-angle layers. The Bayesian approach was used for identifying the optimum number of acoustic classes because cluster validation is carried out prior to class assignment and class outputs are ordinal categorical values. The method is based on the principle that backscatter values from a single incidence angle express a normal distribution for a particular sediment type. The resulting Bayesian classes were well correlated to median grain sizes and the percentage of coarse material. The MLC method uses angular response information from five layers of training areas extracted from the Bayesian classification map. The subsequent PCA analysis is based on the transformation of these five layers into two principal components that comprise most of the data variability. These principal components were clustered in five classes after running an external cluster validation test. In general both methods MLC and PCA, separated the various sediment types effectively, showing good agreement (kappa >0.7) with the Bayesian approach which also correlates well with ground truth data (r²?>?0.7). In addition, sub-bottom data were used in conjunction with the Bayesian classification results to characterize acoustic classes with respect to their geological and stratigraphic interpretation. The joined interpretation of seafloor and sub-seafloor data sets proved to be an efficient approach for a better understanding of seafloor backscatter patchiness and to discriminate acoustically similar classes in different geological/bathymetric settings.  相似文献   

North Atlantic climate impact on early late‐glacial climate oscillations in the south‐eastern Alps inferred from a multi‐proxy lake sediment record     

Roland Schmidt  Kaarina Weckström  Stefan Lauterbach  Richard Tessadri  Kerstin Huber 《第四纪科学杂志》2012,27(1):40-50

High‐resolution multi‐proxy analyses of a sediment core section from Lake Jeserzersee (Saissersee) in the piedmont lobe of the Würmian Drau glacier (Carinthia, Austria) reveal pronounced climatic oscillations during the early late glacial (ca. 18.5–16.0k cal a BP). Diatom‐inferred epilimnetic summer water temperatures show a close correspondence with temperature reconstructions from the adjacent Lake Längsee record and, on a hemispheric scale, with fluctuations of ice‐rafted debris in the North Atlantic. This suggests that North Atlantic climate triggered summer climate variability in the Alps during the early late glacial. The expansion of pine (mainly dwarf pine) between ca. 18.5 and 18.1k cal a BP indicates warming during the so‐called ‘Längsee oscillation’. The subsequent stepwise climate deterioration between ca. 18.1 and 17.6k cal a BP culminated in a tripartite cold period between ca. 17.6 and 16.9k cal a BP with diatom‐inferred summer water temperatures 8.5–10 °C below modern values and a shift from wet to dry conditions. This period probably coincides with a major Alpine glacier advance termed the Gschnitz stadial. A warmer interval between ca. 16.9 and 16.4k cal a BP separates this cold phase from a second, shorter and less pronounced cold phase between ca. 16.4 and 16.0k cal a BP, which is thought to correlate with the Clavadel/Senders glacier advance in the Alps. The following temperature increase, coupled with wet (probably snow‐rich) conditions, caused the expansion of birch during the transition period to the late glacial interstadial. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献