Experimental and numerical study on the fracture of rocks during injection of CO<Subscript>2</Subscript>-saturated water

Geological sequestration of CO₂ into depleted hydrocarbon reserviors or saline aquifers presents the enormous potential to reduce greenhouse gas emission from fossil fuels. However, it may give rise to a complicated coupling physical and chemical process. One of the processes is the hydro-mechanical impact of CO₂ injection. During the injection project, the increase of pore pressures of storing formations can induce the instability, which finally results in a catastrophic failure of disposal sites. This paper focuses mainly on the role of CO₂-saturated water in the fracturing behavior of rocks. To investigate how much the dissolved CO₂ can influence the pore pressure change of rocks, acoustic emission (AE) experiments were performed on sandstone and granite samples under triaxial conditions. The main innovation of this paper is to propose a time dependent porosity method to simulate the abrupt failure process, which is observed in the laboratory and induced by the pore pressure change due to the volume dilatancy of rocks, using a finite element scheme associated with two-phase characteristics. The results successfully explained the phenomena obtained in the physical experiments.     

Paleoseismological evidence for non-characteristic behavior of surface rupture associated with the 2004 Mid-Niigata Prefecture earthquake, central Japan

The 2004 Mid-Niigata Prefecture earthquake sequence (mainshock magnitude, M_JMA 6.8), which occurred in an active fold-and-thrust belt in northern central Japan, generated a small thrust surface rupture (< 20 cm of vertical displacement) along a previously unmapped northern extension of the active Muikamachi–Bonchi–Seien fault zone, on the eastern margin of the epicentral region. To better understand past seismic behavior of the rupture, we conducted a paleoseismic trenching study across the 10-cm-high west-side-up surface rupture at the foot of a pre-existing 1.8-m-high east-facing scarp, which probably resulted from past earthquake(s). A well-defined west-dipping thrust fault zone accompanied by drag folding and displacing the upper Pliocene to lower Pleistocene strata and the unconformably overlying upper Pleistocene (?) to Holocene strata was exposed. The principal fault zone is connected directly to the 2004 surface rupture. From the deformational characteristics of the strata and radiocarbon dating, we inferred that two large paleoseismic events occurred during the past 9000 years prior to the 2004 event. These two pre-2004 events have a nearly identical fault slip (at minimum, 1.5 m), which is ≥ 15 times that of the 2004 event (∼ 10 cm). These paleoseismic data, coupled with the geological and geomorphological features, suggest that the 2004 event represented non-characteristic behavior of the fault, which can potentially generate a more destructive earthquake accompanied by meter-scale surface displacement. This study provides insight into the interpretation of past faulting events and increases our understanding of rupture behavior.     

Evidence for ocean water invasion into the Chicxulub crater at the Cretaceous/Tertiary boundary

Abstract The possibility of ocean water invasion into the Chicxulub crater following the impact at the Cretaceous/Tertiary boundary was investigated based on examination of an impactite between approximately 794.63 and 894.94 m in the Yaxcopoil‐1 (Yax‐1) core. The presence of cross lamination in the uppermost part of the impactite suggests the influence of an ocean current at least during the sedimentation of this interval. Abundant occurrence of nannofossils of late Campanian to early Maastrichtian age in the matrices of samples from the upper part of the impactite suggests that the carbonate sediments deposited on the inner rim margin and outside the crater were eroded and transported into the crater most likely by ocean water that invaded the crater after its formation. The maximum grain size of limestone lithics and vesicular melt fragments, and grain and bulk chemical compositions show a cyclic variation in the upper part of the impactite. The upward fining grain size and the absence of erosional contact at the base of each cycle suggest that the sediments were derived from resuspension of units elsewhere in the crater, most likely by high energy currents association with ocean water invasion.     

Seasonal variation in the energy and water exchanges above and below a larch forest in eastern Siberia

The water and energy exchanges in forests form one of the most important hydro‐meteorological systems. There have been far fewer investigations of the water and heat exchange in high latitude forests than of those in warm, humid regions. There have been few observations of this system in Siberia for an entire growing season, including the snowmelt and leaf‐fall seasons. In this study, the characteristics of the energy and water budgets in an eastern Siberian larch forest were investigated from the snowmelt season to the leaf‐fall season. The latent heat flux was strongly affected by the transpiration activity of the larch trees and increased quickly as the larch stand began to foliate. The sensible heat dropped at that time, although the net all‐wave radiation increased. Consequently, the seasonal variation in the Bowen ratio was clearly ‘U’‐shaped, and the minimum value (1·0) occurred in June and July. The Bowen ratio was very high (10–25) in early spring, just before leaf opening. The canopy resistance for a big leaf model far exceeded the aerodynamic resistance and fluctuated over a much wider range. The canopy resistance was strongly restricted by the saturation deficit, and its minimum value was 100 s m^?1 (10 mm s^?1 in conductance). This minimum canopy resistance is higher than values obtained for forests in warm, humid regions, but is similar to those measured in other boreal conifer forests. It has been suggested that the senescence of leaves also affects the canopy resistance, which was higher in the leaf‐fall season than in the foliated season. The mean evapotranspiration rate from 21 April 1998 to 7 September 1998 was 1·16 mm day^?1, and the maximum rate, 2·9 mm day^?1, occurred at the beginning of July. For the growing season from 1 June to 31 August, this rate was 1·5 mm day^?1. The total evapotranspiration from the forest (151 mm) exceeded the amount of precipitation (106 mm) and was equal to 73% of the total water input (211 mm), including the snow water equivalent. The understory evapotranspiration reached 35% of the total evapotranspiration, and the interception evaporation was 15% of the gross precipitation. The understory evapotranspiration was high and the interception evaporation was low because the canopy was sparse and the leaf area index was low. Copyright © 2001 John Wiley & Sons, Ltd.     

Drainage density, slope angle, and relative basin position in Japanese bare lands from high-resolution DEMs

Relationships between drainage density and slope angle for three bare lands in Japan were analyzed with special attention to channels at early erosion stages and channels in a badland-type terrain. Two of the bare lands were caused by volcanic eruptions 1 or 30–40 years ago, and the other one is a landslide scar formed more than 100 years ago. Raster digital elevation models (DEMs) with a 1-m resolution and ortho aerial photos were generated using digital photogrammetry to enable detailed stream-net extraction and topographic analyses. Data for drainage density, slope angle, and relative height for 88 subwatersheds were obtained from the DEMs and derived stream-nets. The relationship between drainage density and slope angle for each subwatershed can be divided into two types: downward sloping and convex upward. Although previous studies suggested that drainage density positively correlates with slope angle if overland flow is dominant, this correlation seldom occurs in the study areas. The two types of drainage density–slope angle relationships correspond to differing channelization stages that reflect the extension and integration of existing channels, as well as the formation of new low-order streams in response to base-level lowering. The location of subwatersheds within each study area seems to play a major role in determining the stages of channel development and, in turn, the types of drainage density–slope angle relationships.     

Seismicity changes related to a water injection experiment in the Nojima Fault Zone

Abstract A water injection experiment was carried out by the scientific drilling program named the 'Nojima Fault Zone Probe' during the two periods 9–13 February and 16–25 March 1997. The pumping pressure at the surface was approximately 4 MPa. The total amount of injected water was 258 m³. The injection was made between depths of 1480 m and 1670 m in the Disaster Prevention Research Institute, Kyoto University (DPRI) 1800 m borehole drilled into the Nojima Fault zone. A seismic observation network was deployed to monitor seismic activity related to the water injections. Seismicity suddenly increased in the region not far from the injection hole 4 or 5 days after the beginning of each water injection. These earthquakes were likely to be induced by the water injections. Most of the earthquakes had magnitudes ranging from −2 to +1. Numerous earthquakes occurred during the first injection, but only one could be reliably located and it was approximately 2 km north of the injection site. Between the two injection periods, earthquakes concentrated in the region approximately 1 km northwest of the injection site. During and after the second injection experiment, earthquakes were located approximately 1.5 km west of the injection site. Those earthquakes were located approximately 3 km or 4 km from the injection point and between 2 km and 4 km in depth. Values of intrinsic permeability of 10⁻¹⁴–10⁻¹⁵ m² were estimated from the time lapse of the induced seismic activity. The coefficient of friction in the area where the induced earthquakes occurred was estimated to be less than 0.3.     

Occurrence of Free-Gas Associated with Methane-Hydrate-Bearing Formations in the MITI Nankai Trough Wells, Offshore Tokai, Japan

Abstract. The MITI Nankai Trough wells were drilled for exploration of methane-hydrate-bearing sediments in association with seismic inferred bottom simulating reflectors (BSRs). In this project, log data showed low velocity compressional-wave (P-wave) layers below methane-hydrate-bearing formations. Dipole shear sonic acoustic tools (DSI) could not acquire accurate compres-sional velocity in this zone, thus it was not possible to accurately correlate between logging, VSP and surface seismic profiles.
Small amount of gas was presumed to cause the problem in obtaining the low velocity P-wave data. VSP interval velocity data was used to assess the DSI inferred low-velocity layer, which showed lower values than the velocity of the drilling muds. Synthetic seismogram was created by VSP-compensated velocity to compare against corridor stack of VSP. As a result, the depths above and below the methane-hydrate-bearing interval were correlated with synthetic seismograms and reflectivity events on the VSP profiles. By using this correlation technique, distribution of methane-hydrate-bearing formations and free-gas-bearing formations can be determined.     

Deformation history and exhumation process of the Horoman Peridotite Complex, Hokkaido, Japan

The Horoman Peridotite Complex is an Alpine-type orogenic peridotite massif in the Hidaka metamorphic belt, Hokkaido, Northern Japan. Because of wide exposure and extremely limited serpentinization, the complex provides important information on uplift and emplacement processes of an Alpine-type peridotite massif into the crust. Based on microstructures, the massif can be divided into five structural units parallel to the lithological layering as follows; (1) Equigranular Zone, (2) Internal Shear Zone (ISZ), (3) Transition Zone, (4) Porphyroclastic Zone and (5) Basal Shear Zone (BSZ). A top-to-the-north sense of shear deformation in the Porphyroclastic Zone and the Basal Shear Zone implies that the Horoman Peridotite Complex had uplifted from the upper mantle to the lower crust along a northward dipping extensional shear-zone systems. After incorporation of the mantle peridotite with lower crustal rocks, the upper part of the massif (i.e. the Equigranular Zone and the Internal Shear Zone) was overprinted by a top-to-the-south sense of shear deformation that was comparable with the sub-horizontal displacement of the crustal granulite sequences in the Hidaka metamorphic belt under transpressive tectonic environment.