Influence of source region properties on scaling relations forM<0 events

Excavation induced seismic events with moment magnitudesM<0 are examined in an attempt to determine the role geology, excavation geometry, and stress have on scaling relations. Correlations are established based on accurate measurements of excavation geometry and methodology, stress regime, rock mass structure, local tectonics, and seismic locations. Scaling relations incorporated seismic moments and source radii obtained by spectral analysis, accounting for source, propagation, and site effects, and using Madariaga's dynamic circular fault model. Observations suggest that the interaction of stresses with pre-existing fractures, fracture complexity and depth of events are the main factors influencing source characteristics and scaling behaviour. Self-similar relationships were found for events at similar depths or for weakly structured rock masses with reduced clamping stresses, whereas a non-similar behaviour was found for events with increasing depth or for heavily fractured zones under stress confinement. Additionally, the scaling behaviour for combined data sets tended to mask the non-similar trends. Overall, depth and fracture complexity, initially thought as second order effects, appear to significantly influence source characteristics of seismic events withM<0 and consequently favour a non-similar earthquake generation process.     

The experimental research on b value of AE for the rock specimens with pre-existing crack or notch under uniaxial compression

A lot of researches onb value have been made in seismology. Since the 1960’s Mogi, Scholz and others have studied AE of rock specimens in laboratory and discovered that it is related to natural earthquakes. All former researchers used integral specimens to studyb value in the laboratory. However a major earthquake is usually related to a existing seismic-fault in that area. For this reason, a series of fracture experiments with rock and glass specimens having pre-existing crack or notch is performed in order to examine the effect of preexisting crack tob value. The experimental results show that theb value begins to decrease as soon as the initiation of the crack and finally drop to a very low value when the specimen breaks unstably. Based on these, a brief discussion on the possible mechanism ofb value change for natural earthquakes is given. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,9, 393–400, 1987. Projects sponsored by the Chinese Joint Seismological Science Foundation.     

断裂带热释光的初步研究

本文对湖南水口山矿田断裂带岩石进行了热释光分析,断裂作用可对岩石的热释光特征产生显著变化。断裂带岩石天然热释光不存在250℃以下的低温峰;从围岩到断裂带中心,热释光峰值大幅度增加;同一断裂带从地表向深部热释光峰值增加。碎裂断裂带天然热释光峰值为10~3—10~4级,高峰温度位于290—310℃范围内;硅化断裂带的天然热释光峰值大多为10~2—10~3级,高峰温度多在330—350℃范围内。     

四川省青川断裂的特征及形成的物理条件

本文详细地讨论了四川省青川断裂的断裂特征及形成的物理条件。对断裂构造岩的微构造和宏观构造进行了细致的描述分析,将断裂岩分为四个带。结合区域资料,讨论了断裂的形成时代及与邻区的关系。     

非等间距序列的灰色预测模型及其在岩石蠕变断裂中的应用

给出了非等间距序列的灰色预测模型，讨论了该模型的灰色预测平面和最优预测函数。最一将该模型应用于符合Ｋｅｌｖｉｎ模型的岩石材料的蠕变断裂研究，给出了应用实例，并验证了该方法的可行性。     

从钠溪气田看阳新统裂缝型气藏分布复杂性

本文以纳溪气田的实际资料，特别是纳６井裂缝系统的典型资料，裂缝系统在纵横向上的分布具明显的方向性和系统中流体分布异常，剖析川南地区二叠系阳新统裂缝型气藏分布的复杂性，为裂缝型气藏的勘探提供借鉴。     

Hydrofracture mechanisms in rock during pressure grouting

Summary Hydrofracture Mechanisms in Rock During Pressure Grouting. The paper examines the basic meachnisms controlling the initiation of fractures in rocks and layered soils during pressure grouting, and their subsequent propagation into the ground mass. Previous analyses of fracture initiation have tended to concentrate on simplified models in which the ground is treated as an impervious elastic or Mohr-Coulomb continuum. The present method allows for the porous or fissured nature of the ground by considering the effect of seepage forces induced by the pore pressure gradient. The effect is quantified by use of a parameterN such that the ratio of fluid force used in expanding the injection hole, to that used in forcing fluid through void spaces, isN to (1—N).Analysis of hydrofracture propagation is based on stress analysis of a borehole in an elastic continuum, the propagating fracture zone around the borehole being represented as a non-elastic material governed by the Mohr-Coulomb failure criterion. This is supplemented by an energy approach which equates the energy supplied to the ground from the injection pump, with the energy stored in the ground and the energy necessary to fracture it.

---

Zusammenfassung Brucherscheinungen im Fels bei Verprearbeiten. In der Arbeit wird der grundlegende Mechanismus untersucht, welcher bei Verpreßarbeiten in Fels und geschichtetem Boden zur Einleitung und Ausbreitung von Brüchen führt. Frühere Untersuchungen des Bruchbeginns stützten sich im wesentlichen auf vereinfachte Modelle, in denen der Untergrund als undurchlässiges Kontinuum angesehen wird, das entweder elastisch ist oder der Mohr-Coulombschen Bruchbedingung genügt. Die neue Methode berücksichtigt dagegen eine Porosität oder Klüftung des Untergrundes durch Ansatz der vom strömenden Medium auf das Gebirge ausgeübten Belastung. Diese Belastung wird aufgeteilt in einen Druckverlust an der Bohrlochwand (gleich ParameterN mal Verpreßdruck) und die entsprechende, über den gesamten durchströmten Bereich verteilte Belastung.Die Untersuchung der Bruchausbreitung geht von der Spannungsermittlung um ein Bohrloch in einem elastischen Kontinuum aus, wobei in der sich ausbreitenden Bruchzone um das Bohrloch herum nichtelastisches Material angenommen wird, das dem Mohr-Coulombschen Bruchkriterium genügt. Zur Ergänzung dient eine Energie-Betrachtung, bei der die von der Injektionspumpe abgegebene Energie gleichgesetzt wird der im Untergrund gespeicherten Energie und der aufgewendeten Brucharbeit.

---

Résumé Les mécanismes de la fracturation hydraulique dans les roches pendant les injections sous pression. Le mémoire examine les mécanismes fondamentaux qui gouvernent l'initiation des ruptures dans les roches et les sols stratifiés, au cours des injections et leur propagation dans les massifs. Les analyses antérieures de l'initiation de la rupture, se sont concentrées sur des modèles simplifiés où l'on considérait la roche comme un milieu élastique et imperméable, ou comme un milieu de Mohr-Coulomb. La présente méthode admet que le massif est poreux ou fissuré, en considérant l'action des forces de percolation engendrées par le gradient de pression interstitielle. Cette action est quantifiée par un paramètreN, tel que le rapport de la force du liquide employée à dilater le forage d'injection, à celle employée pour forcer le coulis à travers les vides soitN/(1—N). L'analyse de la propagation des ruptures se base sur l'analyse des contraintes autour d'un forage dans un milieu élastique, alors que la zone de la rupture qui se propage autour du forage est représentée par un milieu non-élastique admettant le critère de rupture de Mohr-Coulomb. Cette analyse est complétée par une approche énergétique, où l'énergie qui est fournie au massif par la pompe d'injection est égalée à l'énergie emmagasinée dans la roche et à l'énergie de rupture.

---

Nomenclature A area of new cracks created per unit volume of time t - E total work done by injection fluid - E _i irrecoverable component of energy - E _{i c} work done in fracturing rock or soil - E _{i l} work done to overcome various frictional forces in grouting system - E _{i p} work done to cause plastic deformation of fractured zone - E _{i s} work done to overcome shear strength of fluid during flow - E _{i v} work done to overcome frictional drag between fluid and rock in soil surfaces during flow - E _r recoverable component of energy - E _{r f} elastic strain energy stored in fluid - E _{r s} elastic strain energy stored in rock or soil - h height of overburden - i j 1, 2, 3 - K ₀ coefficient of horizontal earth pressure - k permeability of ground to grout - L length of cylindrical grout source - n rock or soil porosity - p average fluid pressure between timet and (t + t) - p ₀ injection pressure - R radius of grout front - r radial distance from borehole axis - r ₀ radius of borehole - r ₁ radius of fractured zone - S specific surface area of rock or soil - S _T tensile strength of rock or soil - t time - u grout seepage velocity - V volume of grout injected - v volumetric strain - specific surface energy of rock - bulk density of rock or soil - _{i j} ^e elastic strain increment tensor - _{i j} plastic strain increment tensor - v Poisson's ratio - _{i j} average stress tensor in the ground during timet and (t + t) - _R , _T, _Z radial, tangential and vertical stresses induced by grouting - _r , _t , _z radial, tangential and vertical stress around borehole before grouting - grout shear strength - angle of internal friction of rock or soil With 7 Figures     

Monitoring of fault healing after the 1999 Kocaeli, Turkey, earthquake

The North Anatolian fault zone that ruptured during the mainshock of theM 7.4 Kocaeli (Izmit) earthquake of 17 August 1999 has beenmonitored using S wave splitting, in order to test a hypothesisproposed by Tadokoro et al. (1999). This idea is based on the observationof the M 7.2 1995 Hyogo-ken Nanbu (Kobe) earthquake, Japan.After the Hyogo-ken Nanbu earthquake, a temporal change was detectedin the direction of faster shear wave polarization in 2–3 years after the mainshock (Tadokoro, 1999). Four seismic stations were installed within andnear the fault zone at Kizanlik where the fault offset was 1.5 m, about80 km to the east of the epicenter of the Kocaeli earthquake. Theobservation period was from August 30 to October 27, 1999. Preliminaryresult shows that the average directions of faster shear wave polarization attwo stations were roughly parallel to the fault strike. We expect that thedirection of faster shear wave polarization will change to the same directionas the regional tectonic stress reflecting fault healing process. We havealready carried out a repeated aftershock observation at the same site in2000 for monitoring the fault healing process.     

A new scheme for numerical modelling of flow and transport processes in 3D fractured porous media

The objective of this work is to develop a new numerical approach for the three-dimensional modelling of flow and transient solute transport in fractured porous media which would provide an accurate and efficient treatment of 3D complex geometries and inhomogeneities. For this reason, and in order to eliminate as much as possible the number of degrees of freedom, the fracture network, fractures and their intersections, are solved with a coupled 2D–1D model while the porous matrix is solved independently with a 3D model. The interaction between both models is accounted for by a coupling iterative technique. In this way it is possible to improve efficiency and reduce CPU usage by avoiding 3D mesh refinements of the fractures. The approach is based on the discrete-fracture model in which the exact geometry and location of each fracture in the network must be provided as an input. The formulation is based on a multidimensional coupling of the boundary element method-multidomain (BEM-MD) scheme for the flow and boundary element dual reciprocity method-multidomain (BE-DRM-MD) scheme for the transport. Accurate results and high efficiency have been obtained and are reported in this paper.     

Cohesive Crack Analysis of Toughness Increase Due to Confining Pressure

Apparent fracture toughness in Mode I of microcracking materials such as rocks under confining pressure is analyzed based on a cohesive crack model. In rocks, the apparent fracture toughness for crack propagation varies with the confining pressure. This study provides analytical solutions for the apparent fracture toughness using a cohesive crack model, which is a model for the fracture process zone. The problem analyzed in this study is a fluid-driven fracture of a two-dimensional crack with a cohesive zone under confining pressure. The size of the cohesive zone is assumed to be negligibly small in comparison to the crack length. The analyses are performed for two types of cohesive stress distribution, namely the constant cohesive stress (Dugdale model) and the linearly decreasing cohesive stress. Furthermore, the problem for a more general cohesive stress distribution is analyzed based on the fracture energy concept. The analytical solutions are confirmed by comparing them with the results of numerical computations performed using the body force method. The analytical solution suggests a substantial increase in the apparent fracture toughness due to increased confining pressures, even if the size of the fracture process zone is small.