多道瑞利波勘探中频散曲线的提取方法.pdfVIP

多道瑞利波勘探中频散曲线的提取方法 ‘ 鲁来玉，张碧星 中国科学院声学研究所，北京100080 actthedispersioncurvesinMulti-channel RayleighWaveexploration LULai-yu,MANGBi-xing InstituteofAcoustics,TheChineseAcademyofSciences,Beijing,100080,China 1 引 言 由于瑞利波在非均匀介质中的高度频散特征，在表面接收到的不同波长的瑞利波就携带了不 同深度处的信息，野外瑞利波勘探正是基于这一原理，通过发射不同波长的瑞利波信息，由接收 到的瑞利波信号，提取频散曲线，从而根据频散曲线的特征或进行反演运算来推断地表下的介质 情况。因此，获取频散曲线是瑞利波勘探中的一个关键步骤 2 互相关法求取频散曲线的缺点 传统的瑞利波勘探同常是利用两个接收道的信号，采用相关计算来获取瑞利波的频散曲线。 这种方法存在两个明显的缺陷，首先，在层状介质中、瑞利波一般有多个模式传播、并且存在严 重的频散现象，而由相关法的原理可知，它是求取单个非频散模式相速度的一种方法，如果仍用 相关法处理接收到的由各个频散模式叠加形成的信号，得到的结果只能是各个模式传播速度的综 合反映，不能将每一个模式的频散曲线分离开来，严格说，这条频散曲线既不是群速度也不是相 速度;另外，利用相关法求取频散曲线，首先要计算各个频率的瑞利波信号在两道之间的相位差， 其结果是随频率的增加在一1800到++1800度之间成周期性变化的，必须对这样的周期性相位按频 率的增加，每过一个周期加上一个 3600，进行展开，才能进行频散曲线的计算，而对于那些相 位出现周期性变化频点的判断，需要人为的干预选择，这对于更精确的自动反演程序是无法做到 的。为此，人们采用多通道瑞利波勘探，利用所有通道的信息来求取瑞利波的频散曲线。 3 利用波场 f-k域 变换法求取瑞利波频散曲线 接收到的时域波场可由下式给出 N kul er c- dkdou 1 U x,t 仃 D k,al 这里N k,网是与源有关的函数，D k,叻 0即是瑞利波的频散方程，k为波数，。 2可为角 频率。由 1 式可以看出，对多个接收道的时域信号进行双重傅立叶逆变换，就可以得到f一k 域中的波场，在f一k域中，能量的最大值点即对应着瑞利波的解，对其进行处理便可以得到瑞 利波的频散曲线。按照上述方法，我们对一个含有低速夹层的三层介质模型进行了数值模拟运算， 时域中的信号是应用传递矩阵法正演计算的结果。图1给出了利用256道信号计算得到的结果， 图1 a 为厂一k域中的波场，图1 b 是与其对应的位移等高线分布，图中给出了斜率分别为 本项目由国家自然科学重点基金 和声学所所长基金资助 作者简介:鲁来玉 1976- ，男，2001年3月获中南大学地球探测与信息技术专业硕士学位，现为中国科学院声 学研究所博士研究生，主要从事瑞利波勘探理论与应用方面的研究工作 4 0 0 第一层纵波速度vP，横波速度 2400 v、，瑞利波速度V;的直线，可 0 0 0 2000 以看出f一k域中位移的最大 6 0 0 一1600 0 0 2 N 墓 丢1200 、、 0 o e 值点连线，在斜率为v。的直线 o o ， 附近，能量的次极大值点连线 ::: 与斜率为V，的直线重合，由于 。 0 2 s 已 6 计算的量为垂向位移，所以不 k Wavenumber k WavenumbeO ︵ a b 存在SH波型的能量。对f一k ︶ 一 熟 . 、 ﹄ ‘ E 靶 妙 域中的波场进行处理可以得到 0 。 宾 奎至夏垂至至鑫鑫璧重三鑫羞扣 1 . ; A v一f域中的波场图，如图1 … : : : : ! 。 d 所示，可以看出在 心 600Hz以前，位移极大值点与 几 400 日oo izoo 1600 2000 2400 蘸鬓 基阶模式一致，在此频率之后 f Hz 公o. I200 1.00 2000 2400 I HZ 与速度为 3213m/s的直线重 d 合，这正是只有第一层介质时 图1.由256道时域记录得到的.f一k波场和V一f波场 的瑞利波速度，这个结果与正 演计算得到的各个模式的能量最大的频段是一致的，在图1 d 中用虚线给出了各个模式的频散 曲线，在能量比较小的那些频段，可以证明它们都是能陷波，能量集中在低速层内，在表面不能 被接收到。根据图1 d 中的能量极大点可以绘出瑞利波的频散曲线，以此频散曲线的特征或者 进行反演计算，即可对地下介质信息进行推断。 4 实际应用中的问题 扣 三