陕西地区地壳速度结构的层析成象研究

本文利用CT成象技术，根据陕西及邻省的地震波到时资料，反演出本区域地壳的速度结构图，并从中找出速度异常分布区。高速区为陕北地区，处于鄂尔多斯地台内；低速区为关中南部大部分区域及陕南的汉中、石泉地区。最后，对反演结果进行了讨论。     

三维复杂地形近地表速度估算及地震层析静校正

在地表一致性模型的基础上提出一种可适用于宽线剖面、弯曲测线、传统的二维和目前广泛使用的三维地震观测.在地形及近地表低降速带地质结构复杂的探区，低降速带厚度及速度估算的精度是静校正处理的关键.本研究根据三维地震观测的初至走时数据，利用最小平方与QR分解相结合的算法，在三维空间重建近地表低降速带速度模型，根据重建速度模型实现了静校正长波长分量与短波长分量的同步计算.分析了复杂的近地表低降速带模型初至波的性质，在观测值的自动拾取以及理论值的计算中充分考虑了可能成为初至波的直达波、折射波和反射波的利用，提高了低降速带速度模型反演的精度.在初至走时观测数据的拾取中，本研究采用分形算法克服了初至波波形差异以及折射波相位反转导致的拾取误差，实现了三维初至拾取的大规模全自动化运算.在射线路径与初至波理论走时的计算中，本研究采用一种计算量与模型复杂程度无关的三维射线追踪方法，该方法以最小走时射线路径保证了与观测数据有同等意义的初至波的射线追踪及理论走时的计算.野外实际资料的处理结果表明了方法的有效性.     

三维复杂地形近地表速度估算及地震层析静校正

在地表一致性模型的基础上提出一种可适用于宽线剖面、弯曲测线、传统的二维和目前广泛使用的三维地震观测.在地形及近地表低降速带地质结构复杂的探区,低降速带厚度及速度估算的精度是静校正处理的关键.本研究根据三维地震观测的初至走时数据,利用最小平方与QR分解相结合的算法,在三维空间重建近地表低降速带速度模型,根据重建速度模型实现了静校正长波长分量与短波长分量的同步计算.分析了复杂的近地表低降速带模型初至波的性质,在观测值的自动拾取以及理论值的计算中充分考虑了可能成为初至波的直达波、折射波和反射波的利用,提高了低降速带速度模型反演的精度.在初至走时观测数据的拾取中,本研究采用分形算法克服了初至波波形差异以及折射波相位反转导致的拾取误差,实现了三维初至拾取的大规模全自动化运算.在射线路径与初至波理论走时的计算中,本研究采用一种计算量与模型复杂程度无关的三维射线追踪方法,该方法以最小走时射线路径保证了与观测数据有同等意义的初至波的射线追踪及理论走时的计算.野外实际资料的处理结果表明了方法的有效性.     

用于速度图象重建的层析成象法

本文提出了一种用于速度图象重建的层析成象法。特别注意到地震学上的ST与医学上的CT的不尽相同之处,并处理了由此引起的困难。与速度图象重建有关的正问题的分析表明,不同频率范围的资料其分辨能力是不同的,因此,区域地震与远震资料的解释应予注意。本文有关速度图象重建的讨论不仅涉及到与成象方法有关的反演理论和方法,还给出了进行数值计算的内存量和运算量的分析。关于重建图象的质量评价,还给出了反演解的可靠性分析,为分辨特性的描述提供了一种既简便又形象直观的方案。 本文的方法与现行主要方法的比较表明:1.本文的方法允许存在速度间断面,有利于揭示不同构造区域的地壳厚度的差异;2.在常用的方法中,作为正问题假定的常速度块同反演结果解释时非常速度块之间不协调。由于我们把给定网格点的速度值的插值函数作为速度的空间函数,这一困难已被解决;3.相对于Aki等的一般方法,本文的内存量节省约一个量级、运算量减少约一半;4.相对于ART类方法本身不能给出成象结果的可靠性分析,本文的方法则具有显著的优点。 作者在分析了天然地震资料的状况之后相信,对充分利用我国地震台网的现有观测资料,本文的方法是有效的。     

用于速度图象重建的层析成象法

本文提出了一种用于速度图象重建的层析成象法。特别注意到地震学上的ST与医学上的CT的不尽相同之处,并处理了由此引起的困难。与速度图象重建有关的正问题的分析表明,不同频率范围的资料其分辨能力是不同的,因此,区域地震与远震资料的解释应予注意。本文有关速度图象重建的讨论不仅涉及到与成象方法有关的反演理论和方法,还给出了进行数值计算的内存量和运算量的分析。关于重建图象的质量评价,还给出了反演解的可靠性分析,为分辨特性的描述提供了一种既简便又形象直观的方案。 本文的方法与现行主要方法的比较表明:1.本文的方法允许存在速度间断面,有利于揭示不同构造区域的地壳厚度的差异;2.在常用的方法中,作为正问题假定的常速度块同反演结果解释时非常速度块之间不协调。由于我们把给定网格点的速度值的插值函数作为速度的空间函数,这一困难已被解决;3.相对于Aki等的一般方法,本文的内存量节省约一个量级、运算量减少约一半;4.相对于ART类方法本身不能给出成象结果的可靠性分析,本文的方法则具有显著的优点。 作者在分析了天然地震资料的状况之后相信,对充分利用我国地震台网的现有观测资料,本文的方法是有效的。     

中国南北带地壳和上地幔的三维速度图象

本文采用作者提出的地震层析成象法得到了中国南北带地壳和上地幔的三维速度图象。通过误差和分辨分析,以及同爆炸地震测深剖面的比较,证明了成象结果的可靠性。 成象结果表明:1.南北带的地壳和上地幔存在显著的横向不均匀性,深达450km这种不均匀性还依然存在;2.地壳上部的速度图象与地表的已知地质特征明显相关:四川盆地显著低速,康滇地轴显著高速;3.中地壳在很大范围内存在低速层,其最低速度值达5.60km/s;4.在25°N-38°N和100.0°E-103.2°E的长条带内,上地幔顶部出现低速异常,异常速度值约为7.49km/s。 成象结果还清晰地勾划了各块体间的焊接边界。120km深度的速度图象表明,扬子准地台自秦岭以南以龙门-大巴和盐源-丽江台缘褶带为其西部边界;西南以哀劳山褶皱带为界;东南则以右江褶皱带内的南盘江为界。 统计表明,地震活动与南北带的速度结构相关:从20km以上的速度图象发现,大地震大都发生在高速和低速间的过渡条带上。     

中国南北带地壳和上地幔的三维速度图象

本文采用作者提出的地震层析成象法得到了中国南北带地壳和上地幔的三维速度图象。通过误差和分辨分析,以及同爆炸地震测深剖面的比较,证明了成象结果的可靠性。 成象结果表明:1.南北带的地壳和上地幔存在显著的横向不均匀性,深达450km这种不均匀性还依然存在;2.地壳上部的速度图象与地表的已知地质特征明显相关:四川盆地显著低速,康滇地轴显著高速;3.中地壳在很大范围内存在低速层,其最低速度值达5.60km/s;4.在25°N—38°N和100.0°E—103.2°E的长条带内,上地幔顶部出现低速异常,异常速度值约为7.49km/s。 成象结果还清晰地勾划了各块体间的焊接边界。120km深度的速度图象表明,扬子准地台自秦岭以南以龙门—大巴和盐源—丽江台缘褶带为其西部边界;西南以哀劳山褶皱带为界;东南则以右江褶皱带内的南盘江为界。 统计表明,地震活动与南北带的速度结构相关:从20km以上的速度图象发现,大地震大都发生在高速和低速间的过渡条带上。     

用波速层析成像方法提高储层横向预测精度

以数值模拟方法为基础,考察了反射地震层析成像技术对油气储层波速横向检测的效果,然后用渤海油田的地震资料进行了实用试验．对数值模拟结果和实际数据的试验结果进行了分析与评价,认为油气勘探中的反射地震资料具有地震射线近似垂直分布的特点,适合用反射地震层析成像技术进行目标层横向非均匀性检测．     

用波速层析成像方法提高储层横向预测精度

以数值模拟方法为基础，考察了反射地震层析成像技术对油气储层波速横向检测的效果，然后用渤海油田的地震资料进行了实用试验．对数值模拟结果和实际数据的试验结果进行了分析与评价，认为油气勘探中的反射地震资料具有地震射线近似垂直分布的特点，适合用反射地震层析成像技术进行目标层横向非均匀性检测．     

利用地震层析成象技术确定横向非均匀介质中的速度分布和界面位置

本文提出一种横向非均匀介质中的地震层析成象方法，由地震走时资料确定介质的速度分布和界面位置。计算中采用了参数分离、变阻尼、奇异值分解等技术，方程的解较稳定，且可得到解的分办。除反射波外，该方法还可同时利用分析波的资料。因此，适宜于折射及宽角度反射地震测深剖面的资料处理。数字试验结果表明，在模型内部，该方法能得到可靠的反演结果。在边界附近效果较差，受初始值影响较大。     

PULSE DISTORTION IN CAUSTIC REGIONS1

It is well known that there is a phase shift of π associated with the propagation of a seismic pulse through a focus. But this simple behaviour is true only if the two points at which the pulse is observed lie sufficiently far from the focus, one on either side. Near the focus matters are more complicated and the shape of the pulse may be severely distorted because of diffraction. To examine this problem we consider a pulse generated by a point source and focused by a spherical reflector. Using a general combined method of ray tracing and diffraction, we compute seismograms for receivers lying both inside and outside the focal region, and determine how the shape of the pulse changes on propagation through the focus. To assist the physical interpretation of the numerical results, we also present relevant exact and asymptotic results for comparison.     

A STUDY OF MAGNETOTELLURIC STATIC DISTORTION IN THE CONTEXT OF INTRUSIVE VOLCANISM1

MT data from a detailed 84-site grid array (28 × 12 km) in the Paraná Basin (Brazil) reveal a high degree of frequency-independent and parallel amplitude distortion. The Palaeozoic sediments across the survey area are covered by 1.2 km of flood basalts. A deep well-log provides some control regarding the emplacement of thin diabase sills but information regarding vertical feeder dikes is non-existent. The degree of parallel behaviour is identified using anisotropy ratios which quantify the extent and bandwidth of the distortion characteristics. A 2D modelling study is carried out using the concept of a horizontally layered (1D) basin with superimposed, small-scale inhomogeneities. Two likely geological distortion structures are considered. The first represents at- or near-surface (0 to 100 m) inhomogeneities. The second represents thin but vertically elongate dikes. The two distortion effects are found to produce different and characteristic behaviour in the H-polarization (TM) mode model data. The modelled E-polarization (TE) mode data are undistorted and provide control. The model developed to represent vertical dikes gives rise to a characteristic frequency dependence which is distinct from that of the undistorted response. This model also generates subparallel distortion effects which give rise to a form of variance in the response estimates as a function of location. The behaviour occurs, importantly, in both the amplitude and phase of the distorted data. A comparison of the normalized amplitude data, observed and modelled, indicates a high degree of correspondence down to a low frequency limit (< 0.1 Hz) where additional and deep contributions to the observed response become apparent. The modelling results and data characteristics indicate that the survey area is underlain by a series of closely parallel, thin intrusive dikes (a dike swarm).     

砂岩在不同条件下的强度变形特点和P波速度特征

本文介绍了砂岩样品在干燥和饱和状态下,沿平行层理和垂直层理两个不同方向进行单轴加载时的抗压强度、变形特点和垂直于应力方向的P波走时的变化。结果表明,水对砂岩强度、P波速度的影响较大;具有原生层理构造的岩石样品的强度、变形和P波速度显示一定的方向性。所得结果对水库地震、工程地震、矿山开采及地震预报有一定的参考价值。     

TWO METHODS FOR CONTINUOUS MONITORING OF HARMONIC DISTORTION IN VIBROSEIS SIGNALS1

The quality of Vibroseis survey data can be improved by continuously monitoring the vibrator's baseplate and reaction mass accelerations. Equipment failures can be detected as they occur, rather than relying on similarity trials at the beginning and end of the day's production. Equipment faults can then be corrected as they happen and thus would not have a detrimental effect on the quality of the survey data. Source efficiency can be optimized by monitoring the amount of harmonic distortion generated by the vibrator at different drive levels on the different surfaces which may be encountered during a survey. Phase problems introduced by poor coupling of the baseplate to the ground can also be identified and addressed in the field. Rapid analysis of vibrator signals is required if continuous monitoring is to be useful. Frequency-time (f-t) analyses of vibrator signals are often used in processing centres, but are slow and require a large storage capacity which makes the technique unsuitable for a field analysis system. The two methods proposed to analyse vibrator signals entail the use of hodograms and time-varying notch filters. Hodograms provide a qualitative analysis of harmonic distortion and vibrator performance. A fast, time-varying notch filter gives quantitative and qualitative information about the harmonic distortion present in the signal and can be used to identify problems with vibrator behaviour. Both the hodogram and fast, time-varying notch filter methods can analyse the vibrator's reaction mass and baseplate accelerations as it progresses through its sweep and can present automatically interpreted results to the operator before moving to the next vibrator point.     

SOME ASPECTS OF HANDLING VELOCITY INVERSION AND HIDDEN LAYER PROBLEMS IN SEISMIC REFRACTION WORK*

There are two types of masked layers in seismic refraction work: the velocity reversal (low-velocity layer) and the hidden layer (insufficient velocity contrast or layer thickness). On the basis of an analytical formulation of the general case of a masked layer under an overburden of plane and parallel multiple refractors the two limiting cases are discussed: the solution resulting from an uncritical interpretation of the measured time-distance curve and the blind zone solution. Between these two limiting cases there is a variety of possible masked layer solutions. These no-blind zone solutions—as well as the blind zone solution itself—are formulated separately for the velocity inversion and the hidden layer case. For the evaluation of some no-blind zone solution a diagram is presented which can be used for any case of multiple refractors in the overburden of the masked layer. However, it is only for the three- and the four-layer case that a blind zone interpretation by use of diagrams is advisable. Such diagrams are presented together with the basic sets of formulae which contain as parameters only ratios of velocities and layer thicknesses. As the velocity of the masked layer is usually unknown the diagrams are principally constructed to show the dependence on the masked layer velocity. This is useful for estimation of the largest possible error.     

SOME EFFECTS OF MULTIPLE LATERAL IN HOMOGENEITIES IN MAGNETOTELLURICS1

The analytical solution, for the H-polarization magnetotelluric impedance, of a series of multiple, vertical, conducting slabs (dikes) embedded in a host medium is extended to an infinite array in order to model an anisotropic layer. The solution is used to study the effects of such strongly anisotropic media on the surface impedance. At low frequencies such vertically laminated structures behave as a bulk material. It is shown that the effective bulk parameters are those expected from d.c. theory. However, practical electrode separations may not be long enough or adequately positioned to correctly obtain these bulk parameters from the measured impedance. Thus, such structures can masquerade as quite different one-dimensional structures. A haphazard use of long electrode spacings will not necessarily produce correct results.     

REDUCTION OF HARMONIC DISTORTION IN VIBRATORY SOURCE RECORDS*

Due to non-linear effects, the swept frequency signals (sweeps) transmitted into the subsurface by vibrators are contaminated by harmonics. Upon correlation of the recorded seismograms, these harmonics lead to noise trains which are particularly disturbing in the case of down-sweeps. The method described in this paper—which can be regarded as a generalization of Sorkin's approach to the suppression of even order harmonics—allows elimination, from the final vibratory source seismogram, of harmonics of the sweep up to any desired order. It requires that not one single signal but rather a series of M signals is employed where each signal has an initial phase differing from that of the previous one of the series by the phase angle 2πM. Prior to stacking, the seismograms generated with the different signals have to be brought into the form they would have if they had been generated with the same signal. The method seems also to be capable of reducing the correlation noise if sign-bit recording techniques are used.     

SEISMIC TOMOGRAPHIC IMAGING OF A BURIED CONCRETE TARGET1

This paper describes a field evaluation of an algebraic reconstruction technique for the tomographic imaging of sub-surface velocity anomalies. We describe the construction of a three-dimensional concrete model and the acquisition and processing of seismic traveltime data through the model. Image reconstructions of the data sets, using an algebraic reconstruction technique and incorporating prior knowledge are presented and these are compared with the actual model. Reconstructions show that it is essential that accurate data are obtained as we demonstrate that relatively small errors in the traveltime data can seriously degrade the reconstruction. We also show that raypath effects are very important limiting factors to the analysis.     

我国西部地区一些强震及中强震前后波速异常的初步研究（一）——波速比异常

本文对近年来我国西部地区发生的某些强震及中强震（M_S>5.0）的波速比异常作了一些初步探讨.文中给出并讨论了这些地震的波速比异常的一些共同特征,即随时间变化的曲线形态、异常特征量（异常持续时间与异常幅度）和异常的统计分布特征.此外还讨论了波速比异常的特性与地震序列类型、震源深度及震源机制的关系,探索了波速比异常与波速异常间的关系.最后还指出了地震活动性与波速比异常的某种对应关系,并作了简要的讨论.