各向异性介质中反射率法波场模拟及体波辐射图案

本文从一般各向异性的弹性动力学方程出发，采用反射率法来计算任意各向异性的黏弹性层状介质中的全波场合成地震记录。系统矩阵的特征解是在三维空间计算的，考虑了偏离径向平面的波传播。时间假频由复频率技术来压制，空间假频则通过使用Ｆｉｌｏｎ法估计两个波数的积分来避免。文中以裂隙介质和叠层－裂隙复合介质为例，给出了各向异性介质中合成地震记录计算的实用方法，并分析了横波分裂现象及振幅特征。     

垂直地震剖面的KIRCHHOFF积分偏移

把垂直地震剖面(VSP)资料转换到偏移距—深度域,可以通过 VSP CDP叠加来实现,也可用偏移技术来达到。但是,由于 VSP 为共炮点道集,所含数据量较少,所以用常规偏移方法很难获得良好的偏移效果。本文把传统的 Kirchhoff 积分算法推广到变速的叠前 VSP 资料偏移。要把Kirchhoff 积分算法用于 VSP 资料,其关键在于确定 Kirchhoff 轨迹,这相当于确定以源点和接收点为焦点的一等时轨迹。此轨迹与常规地面地震资料所用的常速 Kirchhoff 积分偏移形成的椭圆轨迹不同。本文采用射线方程来确定这种变速 Kirchhoff 轨迹。这种偏移方法的具体实现步骤和一般叠前偏移方法是一样的,主要包括成像时间的确定、近似的 Kirchhoff 积分记录波场外推和求和成像处理等内容。     

由裂隙引起的三分量资料中的泛张各向异性

文中利用有限差分正演模拟方法分别对合裂隙的双层模型（EDA—各向同性）、横向各向同性双层模型（PTL—各自同性）、不含裂隙的双层模型（各向同性—各向同性）等三种情况下的三分量记录进行分析，可以看到三分量记录上存在着由裂隙引起的泛张各向异性。第一种情况，可以看到明显的横渡分裂现象，且在Y分量上出现很强的P-SH波；第二种情况，也有横渡分裂现象，在X分量上可观测到P－SV波，在Y分量上只有单一的横渡信息，却观测不到P－SH波；第三种情况，三分量记录无横渡分裂现象，只在Y分量上观测到独立的横波信息。文中提供了一个含煤层地区的X和Y分量记录，证实在Y分量上的P－SH波是由煤层的裂缝引起的。     

振幅谱比法求取吸收参数的误差

本文分析了当吸收系数与频率不成线性关系时,品质因素Q随频率f变化的情况,并考察了用振幅频比法求取吸收参数的误差。分析结果表明,当吸收系数α=Kf而r≠1(即吸收系数与频率不成线性关系)时,用振幅谱比法求取的Q_m随f的变化比较大,因此而引入的误差也比较大。在这方面,作者对振幅谱比法提出了改进,并在合成地震记录上进行了试算,所得Q_m的质量得到明显改善。     

EDA介质中的弹性波

EDA介质是地壳中广泛存在的一种方位各向异性介质,它是由定向的裂隙所引起。本文通过坐标变换推导出EDA介质弹性张量在一般坐标系下的表达式,并以Chriosteff方程和弹性波方程出发,讨论了平面波在EDA介质传播时的幔度面特征。文中利用有限差分数值计算的结果展示了弹性波在EDA介质中的传播特征,揭示了横波分裂现象及几种弹性参数的地质意义。     

一种新的连续自适应格点反褶积方法

在“格点滤波--子波处理的一种新技术”(《石油地球物理勘探》,21(2)1986:111～124)一文中,已经建立了比较完整的格点滤波理论。由此出发,提出一种连续自适应估算格点滤波参数的新方法。实现此方法的关键在于连续自适应地估计偏相关系数K_m。Makhoul和Griffth分别提出的自适应方法都存在以下三个问题:①由于时间方向上采用逐点往前递推的方法,无法使用高速阵列机;②在估计n时刻的Km时,只利用n时刻以前的信息,所以自适应跟踪速度总是慢于信号的时变速度;③不能对信噪比与分辨率作一折衷的选择。本方法由于采用与输入信息有关的和式比来表示K_m,而且宜于以数据流为单位进行运算,所以这种新的自适应估算方法能够采用阵列机运算,实现快速跟踪的连续自适应格点反褶积。在此反褶积方法中引入加白噪的方法,还能有效地控制信噪比和分辨率。理论研究与实际资料处理表明,这种新方法具有简单、易行、快速、效果稳定的特点。     

关于波动方程偏移——与马在田同志商榷

贵刊从1977年以来,陆续介绍了许多有关波动方程偏移的文章,虽源皆出自克莱鲍特(claerbout)等的文章及有关文献,但理解不一,造成了混乱。有几个问题提出,供讨论。 1.关于用于水平迭加剖面向下延拓的基本方程问题。马氏等在贵刊1977年第5期、第6期和1978年第1期发表的文章中,都用了下述变换来导出上行波方程或在三维偏移中用下述变换来导出上行波方程并用以作为将地震水平迭加剖面向下延拓的基本方程,笔者认为这是不合适的,因为我们看到,在所有文章中的边界条件都是以地面(水平迭加)记录作为τ=0时的边界波场     

对从各向同性到各向异性带来的几个问题的初探

文中根据应力与位移关系,导出了各向异性介质的波动方程(i=1,2,3;k=1,2,3)在二维平面(x-z平面)与铅直方向异性、横向同性(TI)介质的构造走向垂直,且TI介质的对称轴与z轴平行的情况下,可得出二维椭圆各向异性介质中SH波的波动方程,根据平面波理论,由此方程导出了速度函数。文中参考声波方程相移法正、反演原理,提出了各向异性波动方程正、反演的方法,并在各向同性和各向异性情况下,讨论了速度对偏移结果的影响。文中通过实例分析证实,如果被处理数据的假设前提和使用的处理软件的假设前提一致,则用各向异性波动方程偏移,其偏移成像结果是正确的;反之,将会导致明显的误差。即使在弱各向异性情况下,此种误差也是不容忽视的。     

用VSP资料求取吸收参数方法的研究

VSP 资料的直达下行波能量强,信噪比高,为求取吸收参数提供了有利条件。然而,现有的一些求取吸收参数的方法均不能有效地压制上行波及短程多次下行波的影响,因此利用下行波求得的吸收参数精度不高。相比而言,“振幅谱比法”应用较广,而且不受震源能量变化的影响。本文对振幅谱比法作了修正,即利用声波、密度测井曲线构成的反射系数序列对直达下行波的非弹性吸收的振幅进行校正,以消除上行波的干涉和短程多次波的影响,然后再用振幅谱比法求品质因数。这种改进方法称“校后谱比法”。在有高质量的声波、密度测井资料时,使用 VSP 资料可比较准确地求取介质内的吸收参数,并可改善吸收参数的空间分辨率。通过在 SVSP 资料上试算表明,改进的方法是行之有效的。     

横向各向同性介质中的AVO

介质各向异性是影响振幅随炮检距变化的重要因素之一。本文将Ａｋｉ和Ｒｒｃｈａｒｄａ的方法进行推广，推导出两个弹性横向各向同性介质密接条件下平面波反射和透射系数，并给出弱弹性ＴＩ假设下的近似式。近似式的物理意义明确，而且精度明显于前人的结果。