利用高阶交错网格有限差分法模拟地震波在非均匀孔隙介质中的传播

烃类储集层是一种复合多相介质，在固体颗粒的空隙中含有气体或液体. 研究弹性波在该类地层中的传播规律对于油气勘探开发，特别对于全波列声波测井有重要意义. 为了提高孔隙弹性介质数值模拟的计算效率，本文采用改进显式交错网格有限差分算法取代常用的空间域四阶和时间域二阶的速度 - 应力有限差分算法，算法的空间域为八阶、时间域为二阶. 虽然计算的时间步长略小于空间域四阶的情形，但高阶有限差分算法可以选择较粗糙的网格，因此补偿了计算的低效；同时高阶交错网格有限差分算法的空间频散性比低阶算法小. 利用该算法计算了一个两层模型的波场，同时还模拟了等效弹性和孔隙弹性模型中波的传播. 结果表明慢波及其影响明显，尽管慢波衰减很快，但被某一界面反射后，转换形成的P波和S波仍以正常的方式传播，且比慢波衰减小.     

地震台阵滤波及子台分布对聚束效果的影响(Ⅰ)——滤波与相关性的规律

以上海余山实验地震台阵的测震记录为基础，研究了近震事件、远震事件及噪声不同的幅频特征，分析了影响台阵记录之间相关性质的一些主要因素：当子台间距增大时，记录之间相关性降低；当间距固定时，所记录事件的相关性质受信号的频率、震中距(波阵面的入射角)及噪声振动背景的影响，提出了一维时间序列滤波函数的优化结果并不是二维台阵的最佳频率窗，欲取得较好的聚束信噪比，滤波函数不仅考虑压制噪音振动，还要考虑压制信号的非相关频段。     

求解波动方程的任意差分精细积分法

在地震偏移成像技术中，常常要反向外推波场，因此涉及到波动方程的数值求解问题。本文提出了一种求解偏微分方程的新的半解析方法——任意差分精细积分(ADPI)法。其大体思路是：空间域上作坐标离散，但不采用古典差分法的等分离散方式，而是一种相对自由的、任意的离散法；时间域上则采用子域精细积分的方法，既保留了精细积分法的高精度，又克服了工作量大、占用存储大等缺点。该方法具有精度高、带宽小、稳定性好等多项突出优点，并且可以灵括处理各类边界条件。本文从简单一维、较普遍一维、以及二维3种情形讨论该方法对波动方程的具体应用，结合实际的算例，详细分析了各种算法的可行性和精度特点。     

Microscopic observation and contact stress analysis of granite under compression

To study the damage process of microscale to macroscale in coarse-grained granite specimen under uniaxial compressive stress, we have observed micro-damage localization and propagation by using a newly developed experimental system that allows us to observe the damaging process continuously.

The results showed that pre-existing microcracks lead to macroscopic shear fracture through the damage development process. The mechanism of micro-damage initiation in a granite specimen under uniaxial compressive stress may be considered for two cases. One is that two grains such as quartz and feldspar contact each other in the same direction as the axial stress, and the other is that a biotite grain inclined to the axial stress direction is surrounded by feldspar grains. The homogenization theory was applied to verify numerically the micromechanics of stress-induced damage in the mineral contacts. Local stress distribution in the periodic-micro structure was also calculated by the homogenization theory. It is shown that this analysis, which takes into account the initial state of the specimen, is well adapted to the behavior of two grains for which microcracking is the fundamental mechanism of damage.     

Nonlinear surface ground analysis via statistical approach

A statistical approach is proposed for nonlinear surface ground analysis. In contrast to the conventional method which deals with only a single ground motion for equivalent linearization of soil properties, a design response spectrum defined at the upper level (bottom of the surface ground) of an engineering bedrock can be handled as the target design earthquake in the present paper. The effective shear strain in each soil layer is evaluated by means of a statistical procedure in which the mean peak shear strain is computed in terms of its standard deviation and the corresponding peak factor. The stiffness and damping ratio of each soil layer are obtained iteratively from the nonlinear relation of stiffness reduction factors and damping ratios with respect to the strain level. After the evaluation of the equivalent stiffness and damping ratio of every soil layer, the ground surface response spectrum is transformed from the design response spectrum defined at the upper level of the engineering bedrock via the one-dimensional wave propagation theory. The reliability and accuracy of the proposed analysis method is examined through the comparison with the results by the conventional method (represented by the program) for many simulated spectrum-compatible ground motions.     

Magnetic clouds and interplanetary particle transport: a numerical model

Magnetic clouds modify the structure of the interplanetary magnetic field on spatial scales of tenth of AU. Their influence on the transport of energetic charged particles is studied with a numerical model that treats the magnetic cloud as an outward propagating modification of the focusing length. As a rule of thumb, the influence of the magnetic cloud on particle intensity and anisotropy profiles increases with decreasing particle mean free path and decreasing particle speed. Three cases are considered: (1) when the magnetic cloud is the driver of a shock that accelerates particles as it propagates outward, (2) when the magnetic cloud interacts with a prior solar energetic particle event, and (3) when a magnetic cloud already is present in interplanetary space at the time of a solar energetic particle event. In the latter case the cloud acts as a barrier, storing the bulk of the particles in its downstream medium.