二维地震波折射层析成象

给出了一种用地震折射波初至走时确定二维速度的迭代层析成象反演法。该反演法适用于源－接收器间距比常规剖面法密的折射剖面，它是以线性化问题的迭代解为依据，不仅能用于近水平界面的几何形态，而且也能用于连续速度变化的情况。每次迭代中用射击法进行两点间的射线追踪，来建立线性系统。速度场用速度梯度为常数的三角形单元定义，这样射线路径可解析计算。根据正问题的不同线性化公式对两种不同的反演方法作了研究。用线性化走时－速度雅可比行列式反演得出的结果比用更近似于一般级数展开法的慢度公式的结果好。分辨率的例子不仅揭示了源－接收器组合方式及速度结构对分辨率的影响，而且也揭示了射线的几何形态造成的水平向模糊不清和分辨率随深度降低。反演例子表明，使用较好的初始模型时，整体范数产生的解比用最小扰动法计算的解更接近真实模型。对不均匀射线覆盖产生的条纹作用及其排除作了实例说明。     

多震相层析成像与在介质中的多波型射线跟踪

本文介绍了多震相的层析成像的思路和算法，通过穿透和反射走时可以同时作出２维和３维慢度（速度的倒数）重建。我们分析了在穿透和反射数据中确定速度和深度的不确定性，并认识到深度扰动对反射走时异常比慢度扰动更敏感。由不同波类型所提供的对速度和深度的约束，这个算法实际上减少了在一般反射层析成像在速度和深度之间的不确定性，并且也避免了在穿透层析成像中的不确定问题。线性化反演是通过从反射界面深度由分离速度参数迭代进行的。使用一个快速的２－Ｄ和３－Ｄ射线跟踪算法来计算穿透和反射走时和对幔度及反射界面深度的偏导数。深度和速度都用立方Ｂ样条函数来进行参数化。合成例子表明，当同时考虑穿透和反射时间，层析成像的结果得到改进。这个方法也应用到英国煤炭测量局（ＢｒｉｔｉｓｈＣｏａｌＭｅａｓｕｒｅｓ）沿跨线排列所记录的逆ＶＳＰ数据组。通过使用波形配合技术，用同时确定时间延迟和叠加权，可以自动拾取旅行时间。所观察到的逆ＶＳＰ层析成像可比周围介质具有较低速度的两个断层区域成像。断层的位置由附近的反射测线所确定。本文还讨论了在复杂２－Ｄ和３－Ｄ非均匀各向同性介质中地震射线跟踪方法。界面的几何形状和水平速度场都通过使用非均匀步长立方Ｂ－样条节点     

共炮检距二维弹性波地震剖面的偏移方法

共均匀各向同性介质的分界面上入射纵波时，除产生纵波反射Ｐ－Ｐ波外，还会产生转换横波的反射Ｐ－ＳＶ波。同样，当入射波为ＳＶ横波时，在分界面上会同时产生ＳＶ－ＳＶ反射波和ＳＶ－Ｐ反射转换波。对这种二维弹性波地震记录如何进行偏移是探讨的问题。为了能够处理实际地震数据，本文提出在共炮检距剖面上进行二维弹性波记录的偏移方法。     

滇西地区地壳浅部基底速度细结构的研究

本文论述滇西地区思茅－中甸剖面基底导速度细结构的研究结果。基底层项部Ｐｇ界面深度为０－３．５ｋｍ，而底部Ｐ１＾０界面深度为１１．０－１７．０ｋｍ。基底层速度在金河－洱海断裂以南为５．７０－６．３０ｋｍ／ｓ，断裂以北增至６．３０－６．５０ｋｍ／ｓ，其过渡带位于剑川附近，剖面上断裂附近，除界面深度变化外，Ｐｇ面速度横向变化也较明显。速度等值线较为稀疏的景云桥炮南侧、大仓炮附近与支梯炮南侧地区，估计发生     

利用三维弯曲界面的反演方法重建唐山震区莫霍界面三维构造形态

研究利用反射波的走时反演三维弯曲界面和介质层速度的计算方法．各层界面利用分段非完全三次多项式描述．正问题采用一种快速的三维射线追踪方法．反演过程采用变阻尼最小二乘法．数值模拟结果表明，解很快收敛到真模型．处理了通过唐山震区的实测资料，重建了该震区莫霍界面三维构造形态，并揭示了该区域构造与地震活动的关系．      

应县￣淄博剖面重力,地震测深联合解释

利用庆县￣淄博地震测深解释的二维速度结构，并将速度值转换为密度值，提供初始基本模型。运用任意多边形截面的块状二度体模型计算重力异常值，反演二维块状密度分布。计算结果表明：上、中地壳之间的滑脱面是一较明显的密度界面；下地壳底部的强速度梯度层－壳幔过渡带是一高密度层；五台山下部中、下地壳内的低速度块体是低密度块体。     

乌鲁木齐——库尔勒地震转换波测深

利用天然地震转换波研究了天山地区乌鲁木齐－－库尔勒剖面的深部构造特征。地壳中存在５个国煌速度界面，并被数条断裂切割，相应地分为５个构造单元；天山中部的界面均有没和蔼的抬升，在地壳中形成逆部断裂，而地的厚度小于南北两则地区，表明地幔的物质 隆起，有可能造成地壳平均密度增加；天山南缘的和静、巴伦台一带断裂错综复杂，界面严重变形，是未来构造活动比较强烈的地区。     

界面和速度反射联合成像──理论与方法

依据波逆行原理重新推导了由地震波反射走时资料反演界面深度的关系式──走时对界面偏导数关系，它不仅适用于反射波，也适用于透射波．还给出了任意多个复杂界面情况下走时对界面偏导数关系的离散形式．另外，正交算子投影法被推广用来解决速度与界面的同时成像问题．数值模拟的结果表明，本文的方法是有效的．     

共炮检距二维弹性波地震剖面的偏移方法

在均匀各向同性介质的分界面上入射纵波时，除产生纵波反射Ｐ－Ｐ波外，还会产生转换横波的反射Ｐ－ＳＶ波。同样，当入射波为ＳＶ横波时，在分界面上会同时产生ＳＶ－ＳＶ反射波和ＳＶ－Ｐ反射转换波，对这种二维弹性波地震记录如何进行偏移是正在探讨的问题。为了能够处理实际地震数据，本文提出在共炮检距剖面上进行二维弹性波记录的偏移方法。具体做法是，首先用弹性波方程将两分量的地震记录分解为Ｐ－Ｐ（或ＳＶ－ＳＶ）波和Ｐ－ＳＶ（或ＳＶ－Ｐ）波，然后用我们已导出的适合于处理共炮检距剖面的叠前偏移方法对它们进行处理。用模拟数据进行的试验说明方法正确，效果良好。     

利用地震回以资料反演界面位置与速度分布

论述了利用地震回折波资料以界面位置与速度分布的方法，推导了地震波走时对于界面位置偏导数的计算公式。数值模拟和实测资料的计算结果表明了该方法的有效性和编制的计算程序的实用性，该方法最的特点是充分地利用了透射波资料中所含的界面位置的信息。界面位置的分辩率与界面两边的速度反差有关，速度判别越大，则分辨率越高。     

Synthetic Seismograms and Wide-angle Seismic Attributes from the Gaussian Beam and Reflectivity Methods for Models with Interfaces and Velocity Gradients

— The effects of interfaces and velocity gradients on wide-angle seismic attributes are investigated using synthetic seismograms. The seismic attributes considered include envelope amplitude, pulse instantaneous frequency, and arrival time of selected phases. For models with interfaces and homogeneous layers, head waves can propagate which have lower amplitudes, as well as frequency content, compared to the direct arrivals. For media with interfaces and velocity gradients, higher amplitude diving waves and interference waves can also occur. The Gaussian beam and reflectivity methods are used to compute synthetic seismograms for simple models with interfaces and gradients. From the results of these methods, seismic attributes are obtained and compared. It was found that both methods were able to simulate wide-angle seismic attributes for the simple models considered. The advantage of using the Gaussian beam method for seismic modeling and inversion is that it is fast and also asymptotically valid for laterally varying media.     

用微动中的面波信息探测地壳浅部的速度结构

本文介绍了利用地脉动中的面波信息探测地壳浅部横波速度结构的观测方法、仪器设备、野外数据采集和用空间自相关方法测定相速度及相关的数据处理方法，并给出了在两个场地的探测实例. 通过与钻探结果对比表明，由该方法推测的地壳浅部1～2 km 深度范围的浅层横波速度结构，对于确定物性差异大的地层界面是有效和可靠的. 由于该方法使用的是地脉动噪声中的面波信息，不需要专门的震源， 因此具有经济快捷、应用范围广泛等优点，将可能成为一种新的具有很好应用前景的地球物理勘探方法.     

Using the pseudospectral technique on curved grids for 2D acoustic forward modelling1

The Fourier pseudospectral method has been widely accepted for seismic forward modelling because of its high accuracy compared to other numerical techniques. Conventionally, the modelling is performed on Cartesian grids. This means that curved interfaces are represented in a ‘staircase fashion‘causing spurious diffractions. It is the aim of this work to eliminate these non-physical diffractions by using curved grids that generally follow the interfaces. A further advantage of using curved grids is that the local grid density can be adjusted according to the velocity of the individual layers, i.e. the overall grid density is not restricted by the lowest velocity in the subsurface. This means that considerable savings in computer storage can be obtained and thus larger computational models can be handled. One of the major problems in using the curved grid approach has been the generation of a suitable grid that fits all the interfaces. However, as a new approach, we adopt techniques originally developed for computational fluid dynamics (CFD) applications. This allows us to put the curved grid technique into a general framework, enabling the grid to follow all interfaces. In principle, a separate grid is generated for each geological layer, patching the grid lines across the interfaces to obtain a globally continuous grid (the so-called multiblock strategy). The curved grid is taken to constitute a generalised curvilinear coordinate system, where each grid line corresponds to a constant value of one of the curvilinear coordinates. That means that the forward modelling equations have to be written in curvilinear coordinates, resulting in additional terms in the equations. However, the subsurface geometry is much simpler in the curvilinear space. The advantages of the curved grid technique are demonstrated for the 2D acoustic wave equation. This includes a verification of the method against an analytic reference solution for wedge diffraction and a comparison with the pseudospectral method on Cartesian grids. The results demonstrate that high accuracies are obtained with few grid points and without extra computational costs as compared with Cartesian methods.     

基于Hinge损失函数的垂向全变差约束全波形反演

全波形反演可以为叠前深度偏移成像提供更高精度的速度模型,但该方法具有较强的非线性,对初始速度模型的依赖性较强,尤其是在实际应用中,地质条件复杂多变,速度变化不连续,增加了反演非线性程度,常常使反演陷入局部极小值,影响反演的精度.全变差约束在图像去噪领域应用广泛,属于非光滑约束,在去噪过程中能有效的保留图像的不连续界面和边缘信息.本文提出基于Hinge损失函数的垂向全变差约束全波形反演方法,在全变差约束的基础上,利用Hinge损失函数控制模型的更新方向,并使用原-对偶混合梯度算法进行求解,给出这一优化问题的迭代格式,有效提高了对地下不连续界面的重构精度,同时也降低反演对初始速度模型的依赖程度.数值算例证明：与常规全波形反演方法相比,基于全变差约束的全波形反演方法可以有效的重构速度模型中的不连续界面,尤其对高速体边缘的重构效果更明显,但该方法对初始速度模型的依赖性仍然较强;基于Hinge损失函数的垂向全变差约束全波形反演方法降低了对初始速度模型的依赖程度,可以从一个较差的初始模型通过循环迭代的方式最终得到同样精确的速度模型,较好的重构了高速体边缘和不连续界面.     

接收函数方法及研究进展

远震P波波形数据中包含了大量在台站下方地壳上地幔速度间断面所产生的P-S转换波及其多次反射波的信息，是研究台站下方局部区域S速度分布理想的震相，由此产生的接收函数方法是反演台站下方S波速度结构的有效手段。接收函数方法可以通过波形反演拟合接收函数的径向分量，对观测台站下方地球介质的S波速度结构进行估计，也可以通过偏移叠加获得的接收函数道集（地震剖面图）追踪速度间断面。这种方法避免了对天然地震震源及其附近结构混响效应等复杂因素的影响，对S波速度的垂向分布敏感，垂向分辨率高。由于宽频带流动地震台阵的发展，用此方法还可获得研究区域速度结构的横向变化，横向分辨能力主要取决于台站的间距。本文回顾20年来接收函数研究的进展，探讨了方法研究的发展趋势，介绍了对地壳-上地幔结构的部分研究结果。     

Piecewise‐continuous velocity inversion

We suggest a new method to determine the piecewise‐continuous vertical distribution of instantaneous velocities within sediment layers, using different order time‐domain effective velocities on their top and bottom points. We demonstrate our method using a synthetic model that consists of different compacted sediment layers characterized by monotonously increasing velocity, combined with hard rock layers, such as salt or basalt, characterized by constant fast velocities, and low velocity layers, such as gas pockets. We first show that, by using only the root‐mean‐square velocities and the corresponding vertical travel times (computed from the original instantaneous velocity in depth) as input for a Dix‐type inversion, many different vertical distributions of the instantaneous velocities can be obtained (inverted). Some geological constraints, such as limiting the values of the inverted vertical velocity gradients, should be applied in order to obtain more geologically plausible velocity profiles. In order to limit the non‐uniqueness of the inverted velocities, additional information should be added. We have derived three different inversion solutions that yield the correct instantaneous velocity, avoiding any a priori geological constraints. The additional data at the interface points contain either the average velocities (or depths) or the fourth‐order average velocities, or both. Practically, average velocities can be obtained from nearby wells, whereas the fourth‐order average velocity can be estimated from the quartic moveout term during velocity analysis. Along with the three different types of input, we consider two types of vertical velocity models within each interval: distribution with a constant velocity gradient and an exponential asymptotically bounded velocity model, which is in particular important for modelling thick layers. It has been shown that, in the case of thin intervals, both models lead to similar results. The method allows us to establish the instantaneous velocities at the top and bottom interfaces, where the velocity profile inside the intervals is given by either the linear or the exponential asymptotically bounded velocity models. Since the velocity parameters of each interval are independently inverted, discontinuities of the instantaneous velocity at the interfaces occur naturally. The improved accuracy of the inverted instantaneous velocities is particularly important for accurate time‐to‐depth conversion.     

利用地震回折波资料反演界面位置与速度分布

论述了利用地震回折波资料反演界面位置与速度分布的方法，推导了地震波走时对于界面位置偏导数的计算公式。数值模拟和实测资料的计算结果表明了该方法的有效性和编制的计算程序的实用性。该方法最突出的特点是充分地利用了透射波资料中所含的界面位置的信息。界面位置的分辨率与界面两边的速度反差有关，速度差别越大，则分辨率越高。     

SOME SOURCES OF DISTORTION IN TOMOGRAPHIC VELOCITY IMAGES1

Two particular sources of distortion, which may be encountered when applying tomographic imaging techniques to crosshole seismic data, have been investigated. Errors in survey locations of the shots and receivers can produce significant distortions in the images obtained. A simple method for solving simultaneously for the velocity field and shot and receiver location errors is presented and applied to synthetic and real data. Reflection and refraction of rays at velocity interfaces may produce poor density and angular coverage of the rays within the region of interest. It is shown that the effect of the velocity field on the ray coverage can significantly affect the resolution in the velocity image, even if ray bending is taken into account. One consequence of this effect is that, in some cases, little improvement in image quality is achieved by using curvi-ray rather than straight-ray inversion techniques, despite the occurrence of pronounced ray bending.     

朝鲜半岛南部三维地壳速度结构成像

对多震相走时层析成像方法进行改进和优化, 使用朝鲜半岛陆地和近海的地震记录, 对朝鲜半岛南部(345°N～39°N，125°E～130°E)进行地震波走时反演， 获得了本研究区的三维速度结构成像结果.分析了沉积层、基底、上地壳、康腊面、下地壳和莫霍面的特征.研究了不同地质区的深部结构的主要差别、主要断裂特征、 海陆交界的地壳特征等.可以发现京畿山地与沃川褶皱带、太白山褶皱带、临津江褶皱带深部结构的差异尤为明显.     

Characteristic of crustal structure beneath the rifts in southern Tibetan plateau

The method of 2-D travel time inversion,which can be applied to determining 2-D velocity structure and in-terfaces simultaneously,is used in this paper to reprocess the data of Paiku Co-Pumoyingcuo seismic profile across the Ny-ima-Tingri rift and Shenzha-Dinggye rift. P-wave velocity structure and interfaces beneath the profile are obtained. The in-terfaces in the crust near Tingri and Dinggye which are located on rifts have a tendency to uplift,and velocities of middle and lower crusts are high. Low veloc...