Kinematics of throughgoing fractures in jointed rocks

Throughgoing fractures play a major role in subsurface fluid flow yet the kinematics of their formation, which directly impact rock flow properties, are often difficult to establish. We investigate throughgoing fractures in the Monterey Formation of California that developed by the coalescence of pre-existing joints. At Lompoc Landing, throughgoing fractures fall into three main groups: linked, linked with aperture, and breccia zones. The segmented nature of their walls provides numerous piercing points to firmly establish the sense of displacement. Analysis of displacement vectors derived from piercing points demonstrates that the NW–SE trending throughgoing fractures, often interpreted as strike–slip faults, are in fact extensional structures in origin. We suggest that this method may be applied to throughgoing fractures that form by the same mechanism in other geologic settings. Establishing kinematics of throughgoing fractures will lead to a better understanding of their contribution to subsurface fluid flow.     

基于统一强度理论的柱形孔扩张问题研究

考虑了岩土材料的软化和剪胀特性,利用空间轴对称问题的统一强度理论分析了柱形孔扩张问题,求出了岩土材料在柱形孔扩张时的应力场、应变场、位移场和最终扩张压力的统一解表达式。对蒋明镜-沈珠江扩孔公式从不同的方面作了论述。通过算例分析了岩土材料的软化和剪胀特性以及统一强度理论参数对应力场、位移场和扩张压力的影响。已有的关于该问题的Tresca解和Mohr-Coulomb解均为本文的一个特例。     

Dynamic nucleation process of shallow earthquake faulting in a fault zone

Two distinct phases are commonly observed at the initial part of seismograms of large shallow earthquakes: low-frequency and low-amplitude waves following the onset of a P wave ( P ₁) are interrupted by the arrival of the second impulsive phase P₂ enriched with high-frequency components. This observation suggests that a large shallow earthquake involves two qualitatively different stages of rupture at its nucleation.
We propose a theoretical model that can naturally explain the above nucleation behaviour. The model is 2-D and the deformation is assumed to be anti-plane. A key clement in our model is the assumption of a zone in which numbers of pre-existing cracks are densely distributed; this cracked zone is a model for the fault zone. Dynamic crack growth nucleated in such a zone is intensely affected by the crack interactions, which exert two conflicting effects: one tends to accelerate the crack growth, and the other tends to decelerate it. The accelerating and decelerating effects are generally ascribable to coplanar and non-coplanar crack interactions, respectively. We rigorously treat the multiple interactions among the cracks, using the boundary integral equation method (BIEM), and assume the critical stress fracture criterion for the analysis of spontaneous crack propagation.
Our analysis shows that a dynamic rupture nucleated in the cracked zone begins to grow slowly due to the relative predominance of non-coplanar interactions. This process radiates the P₁ phase. If the crack continues to grow, coalescence with adjacent coplanar cracks occurs after a short time. Then, coplanar interactions suddenly begin to prevail and crack growth is accelerated; the P₂ phase is emitted in this process. It is interpreted that the two distinct phases appear in the process of the transition from non-coplanar to coplanar interaction predominance.     

Modelling the compliance of crustal rock—II. Response to temporal changes before earthquakes

There have been several claims that seismic shear waves respond to changes in stress before earthquakes. The companion paper develops a stress-sensitive model (APE) for the behaviour of low-porosity low-permeability crystalline rocks containing pervasive distributions of fluid-filled intergranular microcracks, and this paper uses APE to model the behaviour before earthquakes. Modelling with APE shows that the microgeometry and statistics of distributions of such fluid-filled microcracks respond almost immediately to changes in stress, and that the behaviour can be monitored by analysing seismic shear-wave splitting. The physical reasons for the coupling between shear-wave splitting and differential stress are discussed.
In this paper, we extend the model by using percolation theory to show that large crack densities are limited at the grain-scale level by the percolation threshold at which interacting crack clusters lead to pronounced increases in rock-matrix permeability. In the simplest formulation, the modelling is dimensionless and almost entirely constrained without free parameters. Nevertheless, APE modelling of the evolution of fluid-saturated rocks under stress reproduces the observed fracture criticality and the narrow range of shear-wave azimuthal anisotropy in crustal rocks. It also reproduces the behaviour of temporal variations in shear-wave splitting observed before and after the 1986, M = 6, North Palm Springs earthquake, Southern California, and several other smaller earthquakes.
The agreement of APE modelling with a wide range of observations confirms that fluid-saturated crystalline rocks are stress-sensitive and respond to changes in stress by critical fluid-rock interactions at the microscale level. This means that the effects of changes in stress and other parameters can be numerically modelled and monitored by appropriate observations of seismic shear waves.     

利用遥感监测和分析珠江三角洲的城市扩张过程——以东莞市为例

伴随着经济的快速增长,珠江三角洲近年来出现了迅速的城市化以及乡村城市化的过程。沿着城镇边缘和公路两旁,大片的农田被侵占来单方面地满足城市开发的需要,忽略了粮食产量减少以及环境质量恶化等问题。利用遥感可以有效地监测珠江三角洲的城市扩张过程。在遥感多时相监测的基础上通过熵的计算,可以定量地描述城市扩张的空间规律及扩散过程。     

瞬变电磁场的奇点展开法

在分析电磁系统自然谐振现象的基础上,简述了瞬变电磁场奇点展开法的物理基础、基本理论、极点提取方法,并考查了奇点展开法在地球物理中的应用可能性。     

Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

变差函数参数的计算

为了求出地质统计学中变差函数理论模型的参数，首先把相应模型化成线性函数的形式，多面手应用线性方程组非负解的理论计算出它的参数。