柴达木盆地低頻地震探測的“迴折”波和多次波

在柴达木盆地,利用低頻折射地震的方法,曾經得到許多不同类型的波,其中包括有: （i）沉积岩首波; （ii）基岩首波; （iii）大角度反射波; （iv）“迴折”波; （v）不同类型的多次波; （vi）地壳深界面的反射波. 本文专門討論“迴折”波和不同类型的多次波的性貭,以及辨认他們的主要标志;其他波的性貭,将在另文中討論.     

Reflection and refraction at an imperfectly bonded interface between poroelastic solid and cracked elastic solid

Porous solid is in contact with a cracked elastic solid at a plane interface between them. For the presence of vertically aligned microcracks, the elastic solid behaves transversely isotropic to wave propagation. The coefficients of elastic anisotropy depend on the crack density and crack porosity in the medium. A loose bonding is considered between the two solids so that a limiting case could be the welded contact. At the plane interface, the imperfection in welded bonding is represented by tangential slipping and, hence, results in the dissipation of a part of strain energy. Three types of waves propagate in an isotropic fluid-saturated porous medium, which are considered for incidence at the interface. Incidence of a wave results in three reflected waves and two refracted waves. Partition of incident energy among the reflected and refracted waves is studied for each incidence, varying from normal to grazing directions. Numerical example calculates the energy shares of reflected and refracted waves at the plane interface between water-saturated sandstone and basalt. These energy shares are computed and analyzed for different values of crack parameters as well as loose bonding parameter.     

A METHOD TO DETERMINE THE VELOCITIES OF THE SEAFLOOR AND NEAR-SURFACE SEDIMENTS*

Shotpoint gathers from conventional reflection seismic surveys contain both reflected and refracted waves. In this study shot records were processed and analyzed, and the data were modeled with reflected, refracted, and reflected-refracted waves to fit the recorded data. The result is a detailed velocity model. The inverse problem for refracted waves was solved by using the Wiechert-Herglotz inversion. A 500-km-long 26-fold reflection seismic line from the Barents Sea, north of Norway, has been investigated. The data show high velocities, multiple reflections, and various types of noise. To test the method a total of 34 shot gathers were analyzed along this line. The aim of the interpretation was to determine the velocity in the seafloor and the near-surface sediments. It is possible to map the vertical as well as the lateral velocity distribution in detail. Depending on the length of the streamer and the velocity gradient in the sediments, the calculated depth varies between 300 and 500 m below the seafloor. These velocities were also compared to the stacking velocities obtained from the reflection seismic data to see how the velocities determined by different methods were related. The velocity distribution in the sediments is one of the key factors in seismic interpretation. The technique discussed in this paper can contribute to velocity information both in the processing and interpretation of seismic data.     

REFLECTIONS AND REFRACTIONS FROM CURVED INTERFACES: MODEL STUDY *

Reflections and refractions from curved interfaces were studied on two dimensional scale models. Time of arrival, amplitude and character of reflected, converted, and refracted waves were mainly used for this study. Some reflected refractions, refracted reflections and diffractions were also considered. It was possible to separate PS and SP waves and to study their amplitude and character separately. From the amplitude study of refracted arrivals it was concluded that the refracted ray path penetrates into the high velocity layer rather than propagating along the interface. Although most of the results are interesting from the theoretical point of view, a few applications to exploration problems are suggested.     

Sensitivity analysis and application of time‐lapse full‐waveform inversion: synthetic testing and field data example from the North Sea,Norway

Time‐lapse refraction can provide complementary seismic solutions for monitoring subtle subsurface changes that are challenging for conventional P‐wave reflection methods. The utilization of refraction time lapse has lagged behind in the past partly due to the lack of robust techniques that allow extracting easy‐to‐interpret reservoir information. However, with the recent emergence of the full‐waveform inversion technique as a more standard tool, we find it to be a promising platform for incorporating head waves and diving waves into the time‐lapse framework. Here we investigate the sensitivity of 2D acoustic, time‐domain, full‐waveform inversion for monitoring a shallow, weak velocity change (?30 m/s, or ?1.6%). The sensitivity tests are designed to address questions related to the feasibility and accuracy of full‐waveform inversion results for monitoring the field case of an underground gas blowout that occurred in the North Sea. The blowout caused the gas to migrate both vertically and horizontally into several shallow sand layers. Some of the shallow gas anomalies were not clearly detected by conventional 4D reflection methods (i.e., time shifts and amplitude difference) due to low 4D signal‐to‐noise ratio and weak velocity change. On the other hand, full‐waveform inversion sensitivity analysis showed that it is possible to detect the weak velocity change with the non‐optimal seismic input. Detectability was qualitative with variable degrees of accuracy depending on different inversion parameters. We inverted, the real 2D seismic data from the North Sea with a greater emphasis on refracted and diving waves’ energy (i.e., most of the reflected energy was removed for the shallow zone of interest after removing traces with offset less than 300 m). The full‐waveform inversion results provided more superior detectability compared with the conventional 4D stacked reflection difference method for a weak shallow gas anomaly (320 m deep).     

深地震测深震相识别的一些问题

本文论述了深地震测深工作中一些震相的运动学特点,描述了在水平均匀速度模型的各种不同参数条件下,迴折波、反射波及首波几种震相的走时曲线之间的相互位置及关系,并讨论了界面倾斜时的情况。由此说明对这些震相识别的可能性及依据。最后,列举了随县爆破及永平爆破的应用例子。      

井间地震物理模拟研究

井间地震资料具有极高的分辨率,但是,其波场十分复杂,特别是有效一次反射波场能量较弱,且往往被管波等强相干干扰淹没,波场的识别与分离难度较大.地震物理模拟技术是认识复杂地震波场的有效手段.我们通过单一地质体模型、复杂地质体模型和真实井间模拟模型等多种逼真地质模型的物理模拟,揭示了井间地震观测下直达波、透射波、反射波、折射波,多次波和导波等各类波型的特征.物理模拟试验表明,在炮检域平面上分析直达波的信噪比与能量分布,可以使处理人员对旅行时的误差有更直观的认识.所以识别初至与拾取旅行时间时,不仅要在共震源点道集(CSG)上分析拾取,而且还要在共炮检距道集(COG)上分析拾取,这样就提高了直达波识别的可靠性与拾取时间的精度.     

柴达木东盆地的基岩首波和反射波

一、引言 目前,研究结晶基底,多采用地震首波的方法,利用反射波只能探测地震地质条件较好的浅层基岩面。由于柴达木盆地区内基岩面埋藏很深,界面性质很复杂,过去曾利用重力、磁力、电测方法和中频地震探测结晶基底,均未能取得满意的结果。也曾对西盆     

TheoreticalPKP seismograms near the caustic calculated by sum of reflected waves

Summary This paper gives some examples of theoretical seismograms of PKP waves near the caustic. Seismograms of refracted waves for the original medium are compared with seismograms composed as a sum of the reflected waves, generated at boundaries of a substitute medium. All seismograms are calculated by zero approximation of the ray theory. The influence of some parameters of the source function and of the substitute medium on the results is shown.     

DETERMINATION OF LATERAL INHOMOGENEITIES IN REFLECTION SEISMICS BY INVERSION OF TRAVELTIME RESIDUALS*

Lateral inhomogeneities generate fluctuations in the traveltime of seismic waves. By evaluation of these traveltime fluctuations from different source and receiver positions, lateral inhomogeneities can be located using a pseudo inverse matrix method (Aki, Christoffersson and Husebye 1977). The formulation of the problem is possible for transmitted waves as well as for reflected and refracted waves. In reflection seismics this method is of importance, if no reflections from the inhomogeneities themselves, but only reflections from lower boundaries can be observed. The basic assumptions for the mathematical formulation are (1) the average velocities and depths of the reflecting horizons are known already from standard processing methods, and (2) the traveltime residuals are due to lateral velocity changes between different reflectors or between reflectors and the surface. The area of the earth to be considered is divided into layers and the layers into rectangular blocks. The parallel displacement of a ray after passing a disturbed block is neglected, only the traveltime residual is taken into account. In this paper the method and its application to data obtained with two-dimensional models are described.     

折射波法探测断层的模拟试验——记录图的特征

用折射波法对断层进行勘探的实际工作中,不仅观察到上、下盘的正常折射波及棱上的绕射波,而且还可以看到因绕射、透过所产生的一些次生异常波.这些波郁同时在记录上出现,并构成一幅复杂的波形图.如果不能一一地识别它们,并掌握它们的运动学及动力学特征,就无法正确地识别及解释断层.以往一些文章多半讨论板状介质边缘上的绕射现象,而未涉及上、下盘同时存在时,折射法记录上可能出现的次生绕射、折射波. 本文为了解决这一问题,利用了超声波脉冲地震模拟仪进行了实验,对不同大小断距情况下的折射异常波及次生波进行了研究,回答了有关断层勘探中必须阐明的一些问题.     

Ray path of head waves with irregular interfaces

Head waves are usually considered to be the refracted waves propagating along flat interfaces with an underlying higher velocity. However, the path that the rays travel along in media with irregular interfaces is not clear. Here we study the problem by simulation using a new approach of the spectral-element method with some overlapped elements (SEMO) that can accurately evaluate waves traveling along an irregular interface. Consequently, the head waves are separated from interface waves by a time window. Thus, their energy and arrival time changes can be analyzed independently. These analyses demonstrate that, contrary to the case for head waves propagating along a flat interface, there are two mechanisms for head waves traveling along an irregular interface: a refraction mechanism and transmission mechanism. That is, the head waves may be refracted waves propagating along the interface or transmitted waves induced by the waves propagating in the higher-velocity media. Such knowledge will be helpful in constructing a more accurate inversion method, such as head wave travel-time tomography, and in obtaining a more accurate model of subsurface structure which is very important for understanding the formation mechanism of some special areas, such as the Tibetan Plateau.     

Reflected and refracted waves from a linear transition layer

Summary The problem of a transition layer lying between two homogeneous liquid media is discussed. After obtaining the formal solution for a periodic point source lying in the upper layer, the integrand in the expression for the displacement potential of the upper layer is expanded into series of negative powers of exponentials. Some of the terms in the Bromwich expansion are then evaluated along the Sommerfeld loops which give various reflected and refracted waves. The results for the refracted waves are discussed for the two extreme cases, when frequency is extremely low and extremely high. In both the cases it is found that the frequency dependence for refraction arrivals is the same as expected from a sharp boundary, viz., ^–1. And that for high frequency the travel-time of refraction arrivals is the same as expected from geometric ray theory. Both first- order and second-order discontinuities in density and bulk modulus are considered at the boundaries of the transition layer.     

Using refracted shear waves for velocity estimation

The most difficult part of multicomponent processing is the estimation of the shear-wave velocity map for migration. We used refracted shear waves and a simple iterative method called wavefield continuation (WFC) to evaluate the shallow shear-wave velocity profile on a real data example. The WFC was developed in 1981 by Clayton and McMechan to determine compressional-wave velocity profiles from refracted compressional waves. The application to refracted shear waves is straightforward. The real data example shows that shear structure can be easily determined independently of the compressional structure.     

含有倾斜薄裂缝孔隙地层中的井孔声场

应用三维交错网格应力-速度有限差分方法,数值模拟了含有倾斜裂缝孔隙介质地层中点声源所激发的井孔声场问题.为满足薄裂缝计算需求,开发了不均匀网格有限差分算法,提高了计算精度及计算速度.利用将孔隙介质方程参数取为流体极限的办法来处理裂缝中的流体,实现了流体-孔隙介质界面处的差分方程统一,使界面处的计算更加灵活方便.在验证了方法正确性的基础上,分别考察了单裂缝宽度、裂缝带宽度、裂缝倾斜角度以及孔隙介质渗透率等参数的变化对井轴上阵列波形的影响并进行了分析.结果表明,声波经过裂缝时可能产生反射横波及斯通利波,后者随裂缝宽度的减小而减小,而前者随裂缝宽度的改变,变化不大,在裂缝很小(20μm)时依然存在;裂缝带的宽度、密度越大,反射斯通利波越强;当裂缝(裂缝带)倾斜时,反射横波消失,但反射斯通利波受裂缝倾斜角度的影响较小;渗透率的改变对斯通利波的衰减影响较为明显.     

奇异值分解(SVD)实现地震波场分离与去噪新思路

依据不同性质的地震信号（反射波、折射波、直达波、面波、VSP上\下行波、多次波、随机干扰等）之间在运动学、视速度和相干性上的差异,借助某种数学变换、SVD分解与重构联合、时域与频域结合的方式,通过这些间接的处理手段,把要提取的目标信号或要剔除的干扰信号转换到一种相干性更好的空间域中,再进行SVD分解与重构,最充分利用SVD滤波技术特点,实现地震波场分离与去噪,而不是直接对信号进行SVD分解与重构来实现地震波场分离与去噪,这样做可有效地避免以往对SVD波场分离与去噪技术应用空间狭窄、有效信号损失严重等缺陷性.     

SEISMIC AND GRAVITY TWO-DIMENSIONAL INHOMOGENEOUS MODEL OF THE EARTH'S CRUST (ON THE EXAMPLE OF CENTRAL KAZAKHSTAN)*

The interpretation of deep seismic sounding (DSS) data has been made on the basis of a two-dimensional inhomogeneous model. The refracted first arrivals as well as reflected and diffracted waves on the seismic records have been utilized. The seismic section was modeled in the iso-veolcity lines v(x, y) = const, taking into account the zones of diffraction associated with deep faults. Gravity observations have been used to construct a block model of the Earth's crust with vertical boundaries. It is suggested to define the base of the crust as the zone with velocities between 7.8 and 8.2 km/s. The reflecting boundaries of different length occurring in this zone can be conformal or unconformal with the iso-velocity lines near the base of the crust. As an example of our approach to the interpretation of DSS data the folded-blocky structure of the crust with horizontal inhomogeneities of velocity and density is shown in the Kzyl-Orda-Dzheskazgan profile in Central Kazakhstan.     

柱状多层准弹性介质中声压波形的数值计算与分析——声法测井理论研究（Ⅱ）

为逼近油井实际情况,本文处理圆柱状（任意）多层介质中弹性波的传播问题。采用与文[1]相同的物理模型和二维FFT算法,给出了任意层数柱状准弹性介质中全波声压波形的计算方法和程序。 对多种参数进行了声压波形的数值计算,考察了串槽对全波波形的影响,着重分析了两个胶结面串槽引起初至波的变化。结果表明,串槽导致套管滑行纵波幅度明显降低,固、流界面的反射波是影响串槽时初至波特征的一个重要机制。     

Finite-difference migration of the field of refracted waves in studies of the deep structure of the Earth’s crust and the upper mantle based on the DSS (on the example of the <Emphasis Type="Italic">DOBRE</Emphasis> profile)

The main results of deep seismic sounding (DSS) are usually presented in the form of high-velocity models of the medium. Some model examples and the international DOBRE profile have shown that the informativeness of the data obtained can be significantly enhanced by the construction of wave images of the Earth’s crust, based on the migration of refracted and wide-angle reflected waves. The Donets Basin Refraction/Reflection Experiment (DOBRE) profile crosses the Dnieper-Donets paleorift zone in the Donbas region. Along the profile, refracted waves from the basement and the upper mantle and the reflections from the crust basement (the M boundary) are reliably traced. This wave migration has been used to construct a wave image of the structure of the Earth’s crust. As a result, a clear seismic image of the basement surface, whose depth changes along the profile from 0 to 20 km, was obtained. In near-slope parts of the basin, several major faults were identified that had not been identified previously during standard kinematic data processing. It is shown that the crust-upper mantle transition zone is a clearly reflective horizon only within the crystalline massifs; under a depression, it is represented by a lens-shaped highly-heterogeneous area. As shown in the model examples, the images obtained using such a migration accurately reflect the structural features of the medium, in spite of its complicated structure.     

Combination of common-midpoint-refraction seismics with the generalized reciprocal method

A useful method for increasing the signal/noise ratio of refracted waves is Common-Midpoint (CMP)-refraction seismics. With this technique the shallow underground can be described in detail using all information (amplitude, frequency, phase characteristics) of the wavetrain following the first break (first-break phase). Thus, the layering can be determined and faults, weak zones, and clefts can be identified. This paper deals with the optimization of CMP-refraction seismics used in combination with the Generalized Reciprocal Method (GRM). Theoretical studies show a close relationship of both methods to the kinematics of wave propagation. Velocities and optimum offsets determined by the GRM can be used directly in the partial Radon transformation in CMP-refraction seismics. The integration of refracted waves leads to an increase in the signal/noise ratio but simultaneously the integration boundaries must be restricted to deal only with selective parts of the investigated refractor. The result of this process is an intercept-time section which can be converted directly to a depth section using standard refraction seismic techniques. Another possibility of depth conversion is the transformation of this intercept-time section to a `pseudo-zero-offset section', known from reflection seismics. Thus, zero-offset sections can be migrated using wave-equation techniques such as Kirchhoff migration.