Rock Stress,Rock Stress Measurements,and the Integrated Stress Determination Method (ISDM)

The primary objectives of this work are to (1) improve the understanding of the prevailing stress distribution at the Äspö Hard Rock Laboratory (HRL) in SE Sweden by employing an integrated stress determination approach, and in order to accomplish this (2) extend the existing stress integration methodology denominated integrated stress determination method (ISDM; Cornet in Comprehensive Rock Engineering vol 3, Pergamon Press, Oxford, pp 413–432, 1993a). The new developments of the ISDM involve a 12-parameter representation of the regional stress field in the rock mass (i.e., the full stress tensor and its variation with depth) that is applicable to hydraulic stress data (sleeve fracturing, hydraulic fracturing, and hydraulic tests on pre-existing fractures), overcoring data (CSIR- and CSIRO-type of devices), and to combinations of hydraulic and overcoring stress data. For the latter case, the elastic parameters of the overcoring technique may be solved in situ by allowing the hydraulic stress data to constrain them. As a result, the problem then involves 14 model parameters. Results from the study show that the ISDM effectively improves the precision of the prevailing stress field determination and that it is especially powerful for identification of consistencies/inconsistencies in an existing data set. Indeed, this is the very basic premise and goal of stress integration; combine all available data to achieve as complete a characterization of the mechanical stress model as possible, and not to identify a solution that fits only loosely the maximum amount of stress data.     

Stress Measurement by Overcoring at Shallow Depths in a Rock Slope: the Scattering of Input Data and Results

This paper describes a field experiment of stress measurement using the overcoring method performed in a rock slope, called Rochers de Valabres (located in France’s Southern Alps Region), a field laboratory site prone to rockfalls. Six measurements were conducted at shallow depths from the surface, moving deeper along a sub-horizontal borehole. The experiment was conducted in heterogeneous and anisotropic gneiss, with the overcored rock elastic properties, as evaluated by biaxial and uniaxial tests, being widely scattered. Since stress calculations are sensitive to all input data uncertainties, strain inversion was, thus, performed using an experimental device and Monte Carlo simulations. The experimental device allows the assessment of rather broad confidence intervals for both stress magnitude and orientation. The results indicate that the stress state in the surface area is quite heterogeneous and may be correlated with topography. The measurements show a nonlinear stress distribution with distance to the free surface, along with high values of principal stresses, despite the vicinity of the surface. Although influenced by local topography, orientations of the principal computed stresses are characterized by a high turnover due to local heterogeneities. The results are roughly in accordance to a 2D finite element model of the site.     

Effects of petrophysical parameters on attenuation and dispersion of seismic waves in the simplified poroelastic theory: Comparison between the Biot's theory and viscosity‐extended theory

The simplified macro‐equations of porous elastic media are presented based on Hickey's theory upon ignoring effects of thermomechanical coupling and fluctuations of porosity and density induced by passing waves. The macro‐equations with definite physical parameters predict two types of compressional waves (P wave) and two types of shear waves (S wave). The first types of P and S waves, similar to the fast P wave and S wave in Biot's theory, propagate with fast velocity and have relatively weak dispersion and attenuation, while the second types of waves behave as diffusive modes due to their distinct dispersion and strong attenuation. The second S wave resulting from the bulk and shear viscous loss within pore fluid is slower than the second P wave but with strong attenuation at lower frequencies. Based on the simplified porous elastic equations, the effects of petrophysical parameters (permeability, porosity, coupling density and fluid viscosity) on the velocity dispersion and attenuation of P and S waves are studied in brine‐saturated sandstone compared with the results of Biot's theory. The results show that the dispersion and attenuation of P waves in simplified theory are stronger than those of Biot's theory and appear at slightly lower frequencies because of the existence of bulk and shear viscous loss within pore fluid. The properties of the first S wave are almost consistent with the S wave in Biot's theory, while the second S wave not included in Biot's theory even dies off around its source due to its extremely strong attenuation. The permeability and porosity have an obvious impact on the velocity dispersion and attenuation of both P and S waves. Higher permeabilities make the peaks of attenuation shift towards lower frequencies. Higher porosities correspond to higher dispersion and attenuation. Moreover, the inertial coupling between fluid and solid induces weak velocity dispersion and attenuation of both P and S waves at higher frequencies, whereas the fluid viscosity dominates the dispersion and attenuation in a macroscopic porous medium. Besides, the heavy oil sand is used to investigate the influence of high viscous fluid on the dispersion and attenuation of both P and S waves. The dispersion and attenuation in heavy oil sand are stronger than those in brine‐saturated sandstone due to the considerable shear viscosity of heavy oil. Seismic properties are strongly influenced by the fluid viscosity; thus, viscosity should be included in fluid properties to explain solid–fluid combination behaviour properly.     

页岩气储层脆性评价方法研究现状

随着常规油气资源的不断开发,其存量资源越来越少。美国非常规油气资源的开发成功,为页岩气勘探开发奠定了技术基础。储层改造技术是决定页岩气开发成败的两大关键技术之一。而页岩气储层的脆性指数是评价页岩气储层可改造（可压裂）的关键指标。本文对国内外的页岩气储层脆性评价方法进行了广泛深入的研究,总结了岩石脆性的各种力学表征方法以及矿物成分对岩石脆性的影响,为我国页岩气储层脆性评价方法的研究提供了基础性参考资料。     

利用掩膜和多因子逐步正交变换区分遥感数据中的岩性信息

主成分变换对于区分多光谱遥感数据中岩性信息有一定效果。利用主成分变换的２ 个因子散点图分类提取多光谱遥感数据中较弱岩性信息时,发现这些弱信息并不是出现在散点图重合部分的重心位置,说明多光谱原始遥感数据之间除了存在着起支配作用的主特征值外,还存在反应较弱信息的特征因子。而主成分变换只表示了特征空间系列变换的一种结果,为了提取和显示多因子变化空间中出现的较弱特征信息,开发了多因子逐步正交变换算法软件,在峪耳崖金矿带应用时取得了较好的效果。研究中还采取了掩膜技术排除水体和河道的干扰信息。该方法已用在华北地台北缘其它地区的找矿。     

Geoscientific R&D for high level radioactive waste disposal in Sweden — current status and future plans

SKB (Svensk Kärnbränslehantering AB) is responsible for all handling, transport and storage of the nuclear wastes outside the Swedish nuclear power stations. According to Swedish law, SKB is responsible for an R&D-programme needed to take care of the radwastes. The programme comprises, among others, a general supportive geo-scientific R&D and the Äspö Hard Rock Laboratory (HRL) for more in-situ specific tasks.

Sweden is geologically located in the Fennoscandian shield which is dominated by gneisses and granitoids of Precambrian age. The Swedish reference repository concept thus considers an excavated vault at ca. 500 m depth in crystalline rocks. In this concept (KBS-3), copper canisters with high level waste will be emplaced in deposition holes from a system of tunnels. Blocks of highly compacted swelling bentonite clay are placed in the holes leaving ample space for the canisters. At the final closure of the repository, the galleries are backfilled with a mixture of sand and bentonite. This repository design aims to make the disposal system as redundant as possible. Although the KBS-3 concept is the reference concept, alternative concepts and/or repository lay-outs are also studied. The main alternative, currently under development at SKB, is disposal in boreholes with depths of 4–5 km. The geoscientific research will to a great extent be guided by the demands posed by the performance and safety assessments, as well as the constuctability issues. Some main functions of the geological barrier are fundamental for the long-term safety of a repository. These are: bedrock mechanical stability, a chemically stable environment as well as a slow and stable groundwater flux. The main time-table for the final disposal of long-lived radioactive waste in Sweden foresees the final selection of the disposal system and site during the beginning of next decade.     

THE STABILITY OF ROCK-VENEERED HILLSLOPES

1　ＩＮＴＲＯＤＵＣＴＩＯＮ　Ｕｎｄｅｒｓｔａｎｄｉｎｇｈｉｌｌｓｌｏｐｅｄｅｖｅｌｏｐｍｅｎｔｒｅｍａｉｎｓａｎｅｌｕｓｉｖｅｇｏａｌｆｏｒｅａｒｔｈｓｃｉｅｎｔｉｓｔｓ.Ｔｈｅｏｒｉｅｓｂａｓｅｄｕｐｏｎｏｂｓｅｒｖａｔｉｏｎ,ｌｏｇｉｃ,ａｎｄｓｐ...     

Timing and seismic origin of the historic Luanshibao rock avalanche in the Maoyaba basin,SE Tibetan Plateau: New evidence from 10Be exposure-ages

Rock avalanche is one of the most notable geological disasters in the mountain areas, such as the southeastern Tibetan Plateau. A typical one therein is the Luanshibao (LSB) rock avalanche that occurred in the Maoyaba basin. This rock avalanche has attracted a great deal of attentions, as it has a potential threat to the construction of Sichuan-Tibet Railway. It has been widely accepted that the LSB rock avalanche was caused by a seismic event. However, it is still an open question as to the timing of the earthquake-triggered rock avalanche. Here, we report twenty new ¹⁰Be exposure-ages obtained from the deposition zone. These tightly clustered exposure-ages, combined with geomorphic evidence, indicate that the LSB rock avalanche occurred during the mid-Holocene, possibly at 5.2 ± 0.2 ka. A comparison between the timing of rock avalanche and seismic events suggests a close correlation of the LSB rock avalanche with recurrent earthquakes around ∼5 ka BP. Such a correlation is well supported by the view from previous studies.     

Statistical facies classification from multiple seismic attributes: comparison between Bayesian classification and expectation–maximization method and application in petrophysical inversion

We present here a comparison between two statistical methods for facies classifications: Bayesian classification and expectation–maximization method. The classification can be performed using multiple seismic attributes and can be extended from well logs to three‐dimensional volumes. In this work, we propose, for both methods, a sensitivity study to investigate the impact of the choice of seismic attributes used to condition the classification. In the second part, we integrate the facies classification in a Bayesian inversion setting for the estimation of continuous rock properties, such as porosity and lithological fractions, from the same set of seismic attributes. The advantage of the expectation–maximization method is that this algorithm does not require a training dataset, which is instead required in a traditional Bayesian classifier and still provides similar results. We show the application, comparison, and analysis of these methods in a real case study in the North Sea, where eight sedimentological facies have been defined. The facies classification is computed at the well location and compared with the sedimentological profile and then extended to the 3D reservoir model using up to 14 seismic attributes.     

Effects of fracture intersections on seismic dispersion: theoretical predictions versus numerical simulations

The detection and characterisation of domains of intersecting fractures are important goals in several disciplines of current interest, including exploration and production of unconventional reservoirs, nuclear waste storage, CO₂ sequestration, and groundwater hydrology, among others. The objective of this study is to propose a theoretical framework for quantifying the effects of fracture intersections on the frequency‐dependent elastic properties of fluid‐saturated porous and fractured rocks. Three characteristic frequency regimes for fluid pressure communication are identified. In the low‐frequency limit, fractures are in full pressure communication with the embedding porous matrix and with other fractures. Conversely, in the high‐frequency limit, fractures are hydraulically isolated from the matrix and from other fractures. At intermediate frequencies, fractures are hydraulically isolated from the matrix porosity but can be in hydraulic communication with each other, depending on whether fracture sets are intersecting. For each frequency regime, the effective stiffness coefficients are derived using the linear‐slip theory and anisotropic Gassmann equations. Explicit mathematical expressions for the two characteristic frequencies that separate the three frequency regimes are also determined. Theoretical predictions are then applied to two synthetic 2D samples, each containing two orthogonal fracture sets: one with and another without intersections. The resulting stiffness coefficients, Thomsen‐style anisotropy parameters, and the transition frequencies show good agreement with corresponding numerical simulations. The theoretical results are applicable not only to 2D but also to 3D fracture systems and are amenable to being employed in inversion schemes designed to characterise fracture systems.