水平界面上P-SV转换波转换点的精确解

转换波共转换点的叠加和道集选取都需要准确地计算转换点的位置.Tessmer和Behle、Taylor分别给出了水平单层介质中P SV转换波转换点坐标的解析解,由于其表达式的复杂性,在应用中几乎不被采用.在本文中,运用Snell定律重新建立了在水平反射界面上反射的P SV转换波的转换点坐标的四次方程,并严格地推导出与纵波速度、横波速度、炮检距和反射深度有关的转换点坐标的解析解,确定了惟一的解析表达式.将这一结果应用于P SV转换波的速度分析和叠加处理中.简化的公式有较好的应用价值.     

垂向非均匀介质中转换波转换点计算

转换点位置的计算是转换波资料处理中的一个关键问题. 本文提出了分别基于速度随深度线性变化、速度随垂直走时线性变化、慢度随深度线性变化和慢度随垂直走时线性变化四种等效垂向非均匀介质情况下转换点位置的计算方法. 研究了通过速度拟合、走时近似和相似系数谱三种方式选择合适的等效速度方法. 结合理论模型对非均匀介质转换点计算方法、渐进转换点计算方法、Thomsen近似公式和均匀介质解析计算方法的误差进行了分析,结果表明非均匀介质转换点计算方法能更准确地计算转换点位置.     

多层倾斜层状介质P-SV转换波转换点计算的最佳角度搜索法

转换波转换点的位置对转换波道集的抽取和叠加都非常重要.目前对于单一倾斜反射界面或水平层状介质模型,已经能较好地计算转换波转换点的位置.本文针对多层倾斜层状介质,提出最佳角度搜索法求取P-SV转换波转换点的位置,此方法是通过搜索最符合Snell定律的入射角和反射角来确定最佳的P-SV转换波转换点的位置.在搜索最佳转换点位...     

3D3C陆地反射PS转换波CMP成像影响因素分析

三维三分量(3D3C)陆地反射PS转换波共中心点(CMP)叠加成像方法,虽然抽道集简单,但是对实际资料处理结果往往不理想.尤其当反射界面为三维倾斜界面时,其成像质量较差.本文提出有三个主要因素影响其成像质量:第一,转换点离散.运用实例计算得出,转换点离散度随着纵横波速度比、偏移距和界面倾角的增大而增大.相同界面倾角,不同测线方位的转换点离散度不同,视倾角的绝对值越大离散度也越大;第二,道集内静校正量差异增大.CMP道集中,由于转换点离散使得转换点横向跨度较大,经倾斜界面反射转换的S波出射到近地表地层时的角度差异也较大,导致静校突出;第三,加大动校叠加复杂性.三维倾斜界面PS波CMP道集近炮检距时距方程可表示为双曲形式,但是曲线的顶点位置和动校速度同时随测线方位变化,使得CMP道集同相轴很难校平,动校叠加过程很复杂.     

快速有效的转换波共转换点叠加技术

转换波共转换点(Common Conversion Point简称CCP)叠加的关键在于CCP抽道集和非双曲线正常时差校正NMO(Normal MoveOut). 目前方法的精度限制了其在中-浅层或大炮检距情况下的应用. 我们对CCP叠加技术进行了系统研究,导出了新的CCP位置计算公式和非双曲线时距关系式,并给出了具体的CCP抽道集方法. 理论模型试验和实验资料处理表明本文方法精度高、简便易行,特别是对于中-浅层和大炮检距情况也能得到良好的叠加效果.     

VTI介质中P-SV波转换点与各向异性参数关系

对于转换波地震勘探中的转换点位置这个重要问题,提出转换点位置不仅与纵横波速度比,偏移距深度比以及源检距有关,还与地下介质的各向异性的性质有关,计算了忽略地下介质的各向异性影响对转换点的确定带来的严重误差,从而影响地下地质体的精确成像.通过对层状VTI介质中的转换点近似方程的推导过程,引入该方程不同于传统方程的导出是对层状各向同性介质而言,该方程通过引入各向异性参数,使我们对转换波可以有进一步的认识,拓展了转换波处理中各向异性的应用.该方程对于偏移距深度比小于3.0的情况是比较准确的,这对于大偏移距转换波勘探具有实际意义.     

浅析接收函数共转换点叠加成像方法

回顾了接收函数叠加成像方法的发展过程,发现现在仍被广泛用于接收函数成像的共转换点(CCP)叠加成像方法存在着一定的缺点.通过计算转换波射线路径,得到不同深度、不同震中距转换波射线参数,发现共转换点叠加方法适合于P波接收函数叠加成像,对S波接收函数叠加成像却可能产生错误.结果显示在震中距为60°左右时,真正的转换点位置与共转换点叠加方法追踪到的转换点位置在深度为140km处水平方向偏离约200km,在震中距70°-80°之间时,真正的转换点位置与传统共转换点叠加方法追踪到的转换点位置之间的偏差小于S波接收函数的水平分辨率.利用不同深度产生的转换波射线参数进行反投影,能将转换点归位到更接近于真实的转换点位置上,从而使得通过S波接收函数获取的地球内部结构更加可靠.在岩石圈厚度变化剧烈的地区,建议优先采用震中距在70°-80°之间的远震S波事件进行岩石圈结构的研究.     

火山岩地震屏蔽层的转换波叠前时间偏移成像

在反射地震转换波资料处理中,准确求取共转换点一直是一个难题,采用叠前时间偏移技术能避免共转换点道集的抽取,而且能够使转换波归位到真正的反射点上,实现准确成像.本文针对火山岩地震屏蔽层的转换波成像问题,通过对转换波共近似转换点道集进行速度分析,建立了转换波叠前时间偏移的初始速度场,通过速度扫描和纵、横波速度比值扫描确定最佳的偏移速度场和纵、横波速度比值,实现了在火山岩高速层覆盖区域的转换波偏移成像.实际资料的成像结果表明,本文采用的近似转换点计算以及转换波叠前时间偏移方法是有效的.     

自适应转换波炮域偏移方法的研究与应用

抛开常规的转换点位置估计及极性校正,把检波点作为虚拟震源,利用单程波方程实现转换波记录的炮域偏移处理.由于该方法具有无需寻找转换点位置和进行极性校正的自适应特性,尤其适合于复杂构造转换波记录成像,并且在整个偏移处理流程中减少了估算转换点位置这一环节,有利于偏移精度的进一步提高,同时提高了运算效率.理论模型及实际资料的处理效果论证了该方法的正确性和有效性.     

高分辨率P-SV波速度分析

本文在常规转换波速度分析算法基础上,引入两个速度检测因子,提出了两种高分辨率的转换波速度分析方法,称之为转换波协方差速度分析法和转换波ECM（Enhanced Complex Matched Filter Method）速度分析法,同时对转换点计算公式、共转换点道集选排进行了研究,选取了最佳的转换点精确解计算公式,应用到基于时窗的CCP道集分选方法中,并与高分辨率速度分析同时进行,实现了转换波的高分辨率速度分析的算法,理论模型的计算结果表明,二者在时间及速度上都具有较高的分辨率.     

CONVERSION POINTS AND TRAVELTIMES OF CONVERTED WAVES IN PARALLEL DIPPING LAYERS1

For converted waves, stacking as well as AVO analysis requires a true common reflection point gather which, in this case, is also a common conversion point (CCP) gather. The coordinates of the conversion points for PS or SP waves, in a single homogeneous layer can be calculated exactly as a function of the offset, the reflector depth and the ratio v_p/v_s. An approximation of the conversion point on a dipping interface as well as for a stack of parallel dipping layers is given. Numerical tests show that the approximation can be used for offsets smaller than the depth of the reflector under consideration. The traveltime of converted waves in horizontal layers can be expanded into a power series. For small offsets a two-term truncation of the series yields a good approximation. This approximation can also be used in the case of dipping reflectors if a correction is applied to the traveltimes. This correction can be calculated from the approximated conversion point coordinates.     

COMMON REFLECTION POINT DATA-STACKING TECHNIQUE FOR CONVERTED WAVES1

For converted waves stacking requires a true common reflection point gather which, in this case, is also a common conversion point (CCP) gather. We consider converted waves of the PS- and SP-type in a stack of horizontal layers. The coordinates of the conversion points for waves of PS- or SP-type, respectively, in a single homogeneous layer are calculated as a function of the offset, the reflector depth and the velocity ratio v_p/v_s. Knowledge of the conversion points enables us to gather the seismic traces in a common conversion point (CCP) record. Numerical tests show that the CCP coordinates in a multilayered medium can be approximated by the equations given for a single layer. In practical applications, an a priori estimate of v_p/v_s is required to obtain the CCP for a given reflector depth. A series expansion for the traveltime of converted waves as a function of the offset is presented. Numerical examples have been calculated for several truncations. For small offsets, a hyperbolic approximation can be used. For this, the rms velocity of converted waves is defined. A Dix-type formula, relating the product of the interval velocities of compressional and shear waves to the rms velocity of the converted waves, is presented.     

Dipping-interface Mapping Using Mode-separated Rayleigh Waves

Multichannel analysis of surface waves (MASW) method is a non-invasive geophysical technique that uses the dispersive characteristic of Rayleigh waves to estimate a vertical shear (S)-wave velocity profile. A pseudo-2D S-wave velocity section is constructed by aligning 1D S-wave velocity profiles at the midpoint of each receiver spread that are contoured using a spatial interpolation scheme. The horizontal resolution of the section is therefore most influenced by the receiver spread length and the source interval. Based on the assumption that a dipping-layer model can be regarded as stepped flat layers, high-resolution linear Radon transform (LRT) has been proposed to image Rayleigh-wave dispersive energy and separate modes of Rayleigh waves from a multichannel record. With the mode-separation technique, therefore, a dispersion curve that possesses satisfactory accuracy can be calculated using a pair of consecutive traces within a mode-separated shot gather. In this study, using synthetic models containing a dipping layer with a slope of 5, 10, 15, 20, or 30 degrees and a real-world example, we assess the ability of using high-resolution LRT to image and separate fundamental-mode Rayleigh waves from raw surface-wave data and accuracy of dispersion curves generated by a pair of consecutive traces within a mode-separated shot gather. Results of synthetic and real-world examples demonstrate that a dipping interface with a slope smaller than 15 degrees can be successfully mapped by separated fundamental waves using high-resolution LRT.     

选择相关法提高转换波速度分析精度

利用转换波速度分析拾取的横波速度是转换波共转换点（Common Conversion Point简称CCP）叠加技术中很重要的参数.它的拾取难度在于转换波时距关系在中—浅层或大偏移距情况下误差明显增大,以及在相邻薄反射层上叠加速度谱分辨率很低.本文对选择相关速度谱分析技术进行了研究,即在速度分析中不是利用所有道集,而是有选择的利用部分（可调）道集进行速度分析,研制了相应的实用软件.理论模型和实际资料均表明,该方法在中—浅层或大偏移距情况下,比传统方法有较大改善,分辨率明显提高,能大大提高速度谱估算的精度.     

利用共转换点叠加方法研究华北地区地壳结构

利用华北地震台阵L测线的35个台站记录的895个远震数据进行了接收函数的计算,并利用H-κ叠加方法得到华北克拉通西部陆块东侧和中部陆块内基岩台站下方的地壳结构.利用得到的基岩台站下方的地壳结构和通过波形模拟方法得到的渤海湾盆地的沉积层结构作为背景模型对测线进行共转换点(CCP)叠加成像.在渤海湾盆地,通过增大CCP叠加...     

图像数字去噪处理方法研究及其应用

本文基于地震剖面的数字化图像特征,提出用光栅算子将地震剖面按产状分离成产状相对单一的多幅数字图像剖面。对每一幅产状相对单一的数字图像剖面,用向量分解技术的向量最大特征值对应基向量投影的原理获取有效相干信息的办法实现信噪分离;在结合地质、物探分析的基础上进行最终输出剖面的滤波加权重建,获得信噪比高、波组特征保持好、分辨率相对保证的剖面。实际资料处理结果表明本文提出方法的去噪应用效果显著。     

P‐wave stacking‐velocity tomography for VTI media

A major complication caused by anisotropy in velocity analysis and imaging is the uncertainty in estimating the vertical velocity and depth scale of the model from surface data. For laterally homogeneous VTI (transversely isotropic with a vertical symmetry axis) media above the target reflector, P‐wave moveout has to be combined with other information (e.g. borehole data or converted waves) to build velocity models for depth imaging. The presence of lateral heterogeneity in the overburden creates the dependence of P‐wave reflection data on all three relevant parameters (the vertical velocity V_P0 and the Thomsen coefficients ε and δ) and, therefore, may help to determine the depth scale of the velocity field. Here, we propose a tomographic algorithm designed to invert NMO ellipses (obtained from azimuthally varying stacking velocities) and zero‐offset traveltimes of P‐waves for the parameters of homogeneous VTI layers separated by either plane dipping or curved interfaces. For plane non‐intersecting layer boundaries, the interval parameters cannot be recovered from P‐wave moveout in a unique way. Nonetheless, if the reflectors have sufficiently different azimuths, a priori knowledge of any single interval parameter makes it possible to reconstruct the whole model in depth. For example, the parameter estimation becomes unique if the subsurface layer is known to be isotropic. In the case of 2D inversion on the dip line of co‐orientated reflectors, it is necessary to specify one parameter (e.g. the vertical velocity) per layer. Despite the higher complexity of models with curved interfaces, the increased angle coverage of reflected rays helps to resolve the trade‐offs between the medium parameters. Singular value decomposition (SVD) shows that in the presence of sufficient interface curvature all parameters needed for anisotropic depth processing can be obtained solely from conventional‐spread P‐wave moveout. By performing tests on noise‐contaminated data we demonstrate that the tomographic inversion procedure reconstructs both the interfaces and the VTI parameters with high accuracy. Both SVD analysis and moveout inversion are implemented using an efficient modelling technique based on the theory of NMO‐velocity surfaces generalized for wave propagation through curved interfaces.     

Double plane‐wave reverse‐time migration

We develop a new time‐domain reverse‐time migration method called double plane‐wave reverse‐time migration that uses plane‐wave transformed gathers. Original shot gathers with appropriate data acquisition geometry are double slant stacked into the double plane‐wave domain with minimal slant stacking artefacts. The range of plane‐wave components needed for migration can be determined by estimating the maximum time dips present in shot gathers. This reduces the total number of input traces for migration and increases migration efficiency. Unlike the pre‐stack shot‐profile reverse‐time migration where the number of forward propagations is proportional to the number of shots, the number of forward propagations needed for the proposed method remains constant and is relatively small even for large seismic datasets. Therefore, the proposed method can improve the efficiency of the migration and be suitable for migrating large datasets. Double plane‐wave reverse‐time migration can be performed for selected plane‐wave components to obtain subsurface interfaces with different dips, which makes the migration method target oriented. This feature also makes the method a useful tool for migration velocity analysis. For example, we are able to promptly obtain trial images with nearly horizontal interfaces and adjust velocity models according to common image gathers. Seismic signal coming from steeply dipping interfaces can be included into the migration to build images with more detailed structures and higher spatial resolution as better velocity models become available. Illumination compensation imaging conditions for the proposed method are also introduced to obtain images with balanced amplitudes.