VSP上下行反射波联合成像方法研究

VSP资料上下行波场发育丰富.本文在分析VSP直达波、上行反射波、下行反射波传播路径及其照明范围的基础上,指出了常规VSP波动方程偏移方法缺陷,进而通过修改波场延拓方式,提出了上下行反射波联合成像方法,并在高频近似下分析了该方法的成像原理.该方法不需要进行VSP上下行反射波场分离,能够同时对VSP资料中的一次反射波、自由表面多次波、层间多次波进行成像,比常规成像剖面具有更宽的成像范围和更好的成像效果.该方法能够对下行一次反射波进行成像,从而可以实现常规偏移方法难以处理的高陡倾角构造成像.模拟资料和实际资料处理证明了本文方法的正确性.     

A critical analysis of the concept of planetary wave breaking

The concept of planetary wave breaking (McIntyre andPalmer, 1983; 1984) is critically reviewed. It is concluded that the wave breaking signature is not unique to any particular dynamic event in the stratosphere. Therefore, the classification of stratospheric transport events, such as wave breaking, groups fundamentally different events together. Better qualification of the wave breaking signature and a more solid theoretical basis of planetary wave breaking must be presented if the concept is to be of significant utility in describing stratospheric tracer transport.     

岩石受力状态对弹性波传播速度影响的初步研究

岩石所受的三向应力状态对其弹性波传播性质有一定的影响，此认识为所做的以泥质砂岩、泥岩和硅质胶结的坚硬砂岩为试样的实验所证实。为保证岩石中存在微隙，对部分试样还作了高温预热处理。实验结果表明，泥质砂岩和泥岩的Ｖ＿ｐ／Ｖ＿ｓ值几乎保持不变，但坚硬砂岩的Ｖ＿ｐ／Ｖ＿ｓ值曲线却在应力达破坏荷载的９０％左右时出现下凹现象，经热处理的硬砂岩的这一下凹现象则更突出。     

我国北方500hPa春、夏季平均长波年际异常的EOF分析

利用1958—2008年全球逐日2.5°×2.5°NCEP/NCAR再分析资料,采用经验正交函数分解(EOF)和谐波分析方法对我国北方地区春、夏季平均长波进行了分析。结果表明:(1)我国北方春季平均长波主要表现为亚洲北部从咸海到日本存在一个完整的长波:巴尔喀什湖到我国罗布泊地区以长波脊活动为主,我国东北地区以长波槽活动为主,长波槽、脊的强度均呈增强趋势;夏季平均长波主要表现为乌拉尔山、贝加尔湖及鄂霍次克海地区的长波脊、槽、脊分布,蒙古国到我国北部地区多为长波槽活动,近51年来长波槽存在减弱趋势。(2)春季平均长波的距平变化主要表现为我国西北、华北地区与其上游和下游邻近地区存在反相关关系;夏季平均长波距平变化和春季相比有所东移,即我国北方大部分地区距平变化与其上游和下游地区存在反相关关系。     

A simplified method for simulation of strong ground motion using finite rupture model of the earthquake source

We present a simplified method to simulate strong ground motion for a realistic representation of a finite earthquake source burried in a layered earth. This method is based on the stochastic simulation method of Boore (Boore, D. M., 1983, Bull. Seism. Soc. Am. 73, 1865–1894) and the Empirical Greens Function (EFG) method of Irikura (Irikura, K., 1986, Proceedings of the 7th Japan Earthquake symposium, pp. 151–156). The rupture responsible for an earthquake is represented by several subfaults. The geometry of subfaults and their number is decided by the similarity relationships. For simulation of ground motion using the stochastic simulation technique we used the shapping window based on the kinetic source model of the rupture plane. The shaping window deepens on the geometry of the earthquake source and the propagation characteristics of the energy released by various subfaults. The division of large fault into small subfaults and the method for accounting their contribution at the surface is identical to the EGF. The shapping window has been modified to take into account the effect of the transmission of energy released form the finite fault at various boundaries of the layered earth model above the source. In the present method we have applied the correction factor to adjust slip time function of small and large earthquakes. The correction factor is used to simulate strong motion records having basic spectral shape of ² source model in broad frequency range. To test this method we have used the strong motion data of the Geiyo earthquake of 24th March 2001, Japan recorded by KiK network. The source of this earthquake is modelled by a simple rectangular rupture of size 24 × 15 km, burried at a depth of 31 km in a multilayered earth model. This rupture plane is divided into 16 rectangular subfaults of size 6.0 × 3.75 km each. Strong motion records at eight selected near-field stations were simulated and compared with the observed records in terms of the acceleration and velocity records and their response spectrum. The comparison confirms the suitability of proposed rupture model responsible for this earthquake and the efficacy of the approach in predicting the strong motion scenario of earthquakes in the subduction zone. Using the same rupture model of the Geiyo earthquake, we compared the simulated records from our and the EGF techniques at one near-field station. The comparison shows that this technique gives records which matches in a wide frequency range and that too from simple and easily accessible parameters of burried rupture.     

基于多重互相关函数分析剪切波速

介绍了一种提高单孔法剪切波速测量精度的新方法-多重互相关函数法。考虑到判断场地剪切波到时差的困难,提出了用多重互相关函数来求得两剪切波的到时差。互相关分析可求取两相似信号的到时差,但当两信号相似性比较差时,所求得到时差的误差较大;自相关函数本身可进行多重自相关分析,互相关函数结合多重自相关函数后,也可以进行多重互相关分析来求取两相似信号的到时差。试验表明：相比互相关函数,多重互相关函数可改善触发信号和井中三分向探头输出信号的相似性及更好地抑制噪声,因此能更精确地计算出触发信号和井中三分向探头输出信号的的到时差。基于此原理设计了剪切波波速的测试系统,此系统能自动分析出测点的剪切波速和形成完整的报告。     

模拟地震波传播的三维逐元并行谱元法

高效地震波场正演模拟对于复杂模型中地震波传播与成像研究至关重要.本文在谱元法原理框架内,对已有逐元谱元法改进,提出一种新的逐元并行谱元法求解三维地震波运动方程,并得到地震波场.逐元并行谱元法的核心思想在于在单元上进行质量矩阵与解向量的乘积运算,并将此运算平均分配至每一个CPU计算核心,此处理有利提升谱元法的并行计算效率.同时,根据Gauss-Lobatto-Legendre (GLL)数值积分点与插值点重合的特点,将稠密单元刚度矩阵的存储转化成单元雅克比矩阵行列式的值及其逆的存储,大幅减少谱元法计算内存开销.此外,在模型边界上利用逐元并行谱元法求解二阶位移形式完美匹配层(PML)吸收边界条件,消除边界截断而引入的虚假反射.通过逐元并行谱元法得到的数值解与解析解对比,以及实际地震波场模拟,数值结果证实了逐元并行谱元法用于地震波场模拟的高效性.

     

岩质高边坡岩体变形参数及松弛带厚度研究

预测岩质高边坡开挖后岩体变形模量的变化及松弛带厚度，是分析岩质高边坡在开挖后变形（位移）和作好防护设计的重要资料，运用波动力学关于平均应力与体积模量、岩体纵波速度与弹性模量、变形模量间的关系，通过部分实测资料及边坡应力场有限元分析的资料，分别建立了纵波速度与岩体变形模、岩体应力间的关系，研究了开挖边坡岩体变形模量的变化，预测了岩体松弛带的厚度。     

A COMBINED REFRACTION-DIFFRACTION-DISSIPATION MODEL OF WAVE PROPAGATION

A numerical model based on the mild-slope equation of water wave propagation over complicated bathymetry,taking into account the combined effects of refraction,diffraction and dissipation due to wavebreaking is presented.Wave breaking is simulated by modifying the wave height probability density func-tion and the wave energy dissipation mechanism is parameterized according to that of the hydraulic jumpformulation.Solutions of the wave height,phase function,and the wave direction at every grid point areobtained by finite difference approximation of the governing equations,using Gauss-Seidel Iterative Method(GSIM)row by row.Its computational convenience allows it to be applied to large coast regions tostudy the wave transformation problem.Several case studies have been made and the results compare verywell with the experiment data and other model solutions.The capability and utility of the model forreal coast areas are illustrated by application to a shallow bay of northeast Australia.     

适用于声波方程数值模拟的时间-空间域隐式有限差分算子优化方法

声波方程数值模拟已广泛应用于理论地震计算,同时构成了地震逆时偏移成像技术的基础.对于有限差分法而言,在满足一定的稳定性条件时,普遍存在着因网格化而形成的数值频散效应.如何有效地缓解或压制数值频散是有限差分方法研究的关键所在.有限差分格式分为显式有限差分和隐式有限差分.隐式有限差分能够进一步压制数值频散效应.因此本文提出了给定频率范围满足时间-空间域隐式有限差分频散关系的方法,并根据震源频率、波速和网格间距确定波数范围,在此基础上建立方程确定了相应的隐式有限差分系数,使得差分系数能在更大频率范围符合波场传播规律.通过频散分析和正演模拟,验证了本文方法的有效性.