An Approach for Modeling Rock Discontinuous Mechanical Behavior Under Multiphase Fluid Flow Conditions

In this paper, the two computer codes TOUGH2 and RDCA (for “rock discontinuous cellular automaton”) are integrated for coupled hydromechanical analysis of multiphase fluid flow and discontinuous mechanical behavior in heterogeneous rock. TOUGH2 is a well-established code for geohydrological analysis involving multiphase, multicomponent fluid flow and heat transport; RDCA is a numerical model developed for simulating the nonlinear and discontinuous geomechanical behavior of rock. The RDCA incorporates the discontinuity of a fracture independently of the mesh, such that the fracture can be arbitrarily located within an element, while the fluid pressure calculated by TOUGH2 can be conveniently applied to fracture surfaces. We verify and demonstrate the coupled TOUGH–RDCA simulator by modeling a number of simulation examples related to coupled multiphase flow and geomechanical processes associated with the deep geological storage of carbon dioxide—including modeling of ground surface uplift, stress-dependent permeability, and the coupled multiphase flow and geomechanical behavior of fractures intersecting the caprock.     

PAH in fish bile detected by fixed wavelength fluorescence

Fixed wavelength fluorescence (FF) for detection of polyaromatic hydrocarbon (PAH) metabolites in fish bile has been studied, and is presented through a laboratory and a field experiment. The method, which can be regarded as a biomarker of exposure to PAH, has proved to be a simple, rapid and cost effective method for monitoring PAH contamination.     

Multi-scale permeability estimation for a tropical catchment

Physically based and spatially distributed modelling of catchment hydrology involves the estimation of block or whole-hillslope permeabilities. Invariably these estimates are derived by calibration against rainfall–runoff response. Rarely are these estimates rigorously compared with parameter measurements made at the small scale. This study uses a parametrically simple model, TOPMODEL, and an uncertainty framework to derive permeability at the catchment scale. The utility of expert knowledge of the internal catchment dynamics (i.e. extent of saturated area) in constraining parameter uncertainty is demonstrated. Model-derived estimates are then compared with core-based measurements of permeability appropriately up-scaled. The observed differences between the permeability estimates derived by the two methods might be attributed to the role of intermediate scale features (natural soil pipes). An alternative method of determining block permeabilities at the intermediate or hillslope scale is described. This method uses pulse-wave tests and explicitly incorporates the resultant effects of phenomena such as soil piping and kinematic wave migration. The study aims to highlight issues associated with parameterizing or validating distributed models, rather than to provide a definitive solution. The fact that the permeability distribution within the Borneo study catchment is comparatively simple, assists the comparisons. The field data were collected in terrain covered by equatorial rainforest. Combined field measurement and modelling programmes are rare within such environments. © 1998 John Wiley & Sons, Ltd.     

Thermo-hydrogeomechanical modeling of vertical ground deformation during the operation of the Mutnovskii Geothermal Field

Vertical ground displacements at the Mutnovskii Geothermal Field were measured annually during the 2004–2013 period using a network of borehole markers. Three areas have been identified where the vertical displacements were different in character. Positive vertical deformation (2–5 mm/yr) was recorded in the middle of the Dachnyi Area in 2005–2006, giving way to a later stabilization. The North Test Area where the exhaust heat carrier was reinjected did not show any substantial deformation during the 2003–2006 period. This was followed by an uplift in 2006–2008 (6–7 mm/yr) and then by a subsidence (5–8 mm/yr) in 2009–2013. No vertical deformation worth mentioning was recorded in the Verkhne-Mutnovskii Area prior to 2008, but a rapid subsidence began at a rate of 6–18 mm/yr after 2008. We used the TOUGH-FLAC software to analyze vertical ground deformation. This program is used to perform thermo-hydrogeomechanical (THM) modeling based on the previous TOUGH2 model. THM modeling can explain relative vertical deformation during the exploitation of the Mutnovskii Geothermal Field by separating the geothermal field into two compartments by the Osnovnoi fault, as the two have different tectonic settings and petrophysical properties.     

Constitutive Relationships for Elastic Deformation of Clay Rock: Data Analysis

Geological repositories have been considered a feasible option worldwide for storing high-level nuclear waste. Clay rock is one of the rock types under consideration for such purposes, because of its favorable features to prevent radionuclide transport from the repository. Coupled hydromechanical processes have an important impact on the performance of a clay repository, and establishing constitutive relationships for modeling such processes are essential. In this study, we propose several constitutive relationships for elastic deformation in indurated clay rocks based on three recently developed concepts. First, when applying Hooke’s law in clay rocks, true strain (rock volume change divided by the current rock volume), rather than engineering strain (rock volume change divided by unstressed rock volume), should be used, except when the degree of deformation is very small. In the latter case, the two strains will be practically identical. Second, because of its inherent heterogeneity, clay rock can be divided into two parts, a hard part and a soft part, with the hard part subject to a relatively small degree of deformation compared with the soft part. Third, for swelling rock like clay, effective stress needs to be generalized to include an additional term resulting from the swelling process. To evaluate our theoretical development, we analyze uniaxial test data for core samples of Opalinus clay and laboratory measurements of single fractures within macro-cracked Callovo-Oxfordian argillite samples subject to both confinement and water reduced swelling. The results from this evaluation indicate that our constitutive relationships can adequately represent the data and explain the related observations.     

High matrix interference affecting detection of PAH metabolites in bile of Atlantic hagfish (Myxine glutinosa) used for biomonitoring of deep-water oil production

The characteristic biology and wide distribution of hagfish species makes them relevant for use in pollution biomonitoring at great water depths, particularly in regions where deep-water oil production may take place. The exposure of fish to petrogenic contaminants can normally be detected from the level of polycyclic aromatic hydrocarbon (PAH) metabolites in bile fluid. Some of these metabolites are strong fluorophores, allowing analytical detection by means of simple fluorometric techniques such as fixed wavelength fluorescence (FF) and synchronous fluorescence scanning (SFS). In the present study bile from Atlantic hagfish (Myxine glutinosa) collected in pristine areas (Barents Sea and southwestern Norway) displayed strong bile fluorescence levels, suggesting the presence of PAH contaminants. However, gas-chromatography-mass spectrometry (GC-MS) analyses ruled out PAHs as the origin for this fluorescence signal. Rather, the bile of Myxine contains components resulting in unusually strong background fluorescence interfering at the wavelength pairs used for detection of PAH metabolites. Possible background for the observed matrix interference and implications for detection of PAH metabolites in hagfish is discussed.     

Probabilistic Analysis of Fracture Reactivation Associated with Deep Underground CO2 Injection

In the context of carbon capture and storage, deep underground injection of CO₂ induces the geomechanical changes within and around the injection zone and their impact on CO₂ storage security should be evaluated. In this study, we conduct coupled multiphase fluid flow and geomechanical modeling to investigate such geomechanical changes, focusing on probabilistic analysis of injection-induced fracture reactivation (such as shear slip) that could lead to enhanced permeability and CO₂ migration across otherwise low-permeability caprock formations. Fracture reactivation in terms of shear slip was analyzed by implicitly considering the fracture orientations generated using the Latin hypercube sampling method, in one case using published fracture statistics from a CO₂ storage site. The analysis was conducted by a coupled multiphase fluid flow and geomechanical simulation to first calculate the three-dimensional stress evolution during a hypothetical CO₂ injection operation and then evaluate the probability of shear slip considering the statistical fracture distribution and a Coulomb failure analysis. We evaluate the probability of shear slip at different points within the injection zone and in the caprock just above the injection zone and relate this to the potential for opening of new flow paths through the caprock. Our analysis showed that a reverse faulting stress field would be most favorable for avoiding fracture shear reactivation, but site-specific analyses will be required because of strong dependency of the local stress field and fracture orientations.     

Uptake and tissue distribution of C4-C7 alkylphenols in Atlantic cod (Gadus morhua): relevance for biomonitoring of produced water discharges from oil production

The sensitivity of different tissues for assessment of chronic low-dose environmental exposure of fish to alkylphenols (APs) was investigated. We exposed Atlantic cod (Gadus morhua) in the laboratory to tritium labelled 4-tert-butylphenol, 4n-pentylphenol, 4n-hexylphenol, and 4n-heptylphenol via seawater (8 ng/l) and via contaminated feed (5 μg/kg fish per day). Measurements of different fish tissues during eight days of exposure and eight subsequent days of recovery revealed that APs administered via spiked seawater were readily taken up whereas the uptake was far less efficient when APs were administered in spiked feed. AP residues were mainly located in the bile fluid whereas the concentrations in liver were very low, indicating a rapid excretion and the liver-bile axis to be the major route of elimination. The biological half-life of APs in the exposed cod was short, between 10 and 20 h. Our study shows that in connection with biomonitoring of AP exposure in fish, assessment of AP metabolites in bile fluid is a more sensitive tool than detection of parent AP levels in liver or other internal tissues.