一种两道多模式瑞利波地下探测系统及方法技术方案

本发明专利技术提供了一种两道多模式瑞利波地下探测系统及方法；系统包括：人工锤击震源、多分量检波器、采集站主机、电脑以及若干数据传输线；方法为基于该系统采集分析和处理多模式瑞利波信号进行地下结构探测：首先利用垂向和水平向信号分析瑞利波不同模式的质点振动轨迹特征及其相位角曲线斜率差异，并据此对基阶模式和第一高阶模式瑞利波信号进行分离，进而在两个模式互不干扰的情况下提取基阶和第一高阶频散曲线，最后联合反演多模式频散曲线获得地层结构，实现对地下介质的探测。本发明专利技术提供的系统及方法利用两道多分量检波器采集瑞利波信号并对地下结构进行探测，减少了采集设备，降低了施工难度，同时提高了瑞利波探测的精度和准确性。

【技术实现步骤摘要】
一种两道多模式瑞利波地下探测系统及方法
本专利技术属于地震面波探测方法领域；尤其涉及一种两道多模式瑞利波地下探测系统及方法。
技术介绍
瑞利波是一种沿着地球自由表面扩散的地震波场，在传播过程中其不同频率成分的波动传播速度不同，具有频散的特性。理论与实践均已表明瑞利波的频散特征受其传播路径下方介质速度结构的控制，利用瑞利波的频散特性能够获得介质的内部结构，因此，瑞利波探测手段广泛应用于工程勘察、资源勘探、地球内部结构成像等不同领域中。瑞利波具有多模式发育的特点，即某一频率会对应两个甚至多个传播速度(相速度)，通常将不同频率对应的以最低速度传播的部分称为瑞利波的基阶模式，速度较高一些的称为瑞利波的高阶模式，若有多个较高的传播速度则依次称为第一高阶模式瑞利波、第二高阶模式瑞利波等等，实际中通常能够观测到基阶和第一高阶模式。近年来的研究表明，瑞利波基阶模式的质点振动轨迹为“逆进”的椭圆，而第一高阶模式的质点振动轨迹为“顺进”的椭圆，二者在时间域记录中混叠在一起，难以区分，这就为后续的瑞利波数据处理及地下结构探测带来了影响，造成相位计算误差、频散能量分辨率降低以及探测精度下降等问题。因此，在实际应用中，瑞利波探测方法往往采用多道记录的方式，利用f-k变换、高分辨率线性拉东变换、τ-p变换等方法将瑞利波各模式的“频散核”从时间-距离域变换至频率-速度域，使得基阶和第一高阶模式的瑞利波能量相互分离，从而提高后续数据处理的精度和探测的效果。然而上述方法的实施依赖于多道记录，在采集道较少时，无法很好的区分不同模式能量；并且实际中的瑞利波能量在不同模式之间是不均衡，当某个模式能量占主导时，较弱的模式在频散能量谱中往往不易识别(或在部分频段不易识别)。同时，多道采集方式也具有仪器设备较多、野外施工不便等问题。因此，如何基于更少的采集道数实现更高精度的瑞利波数据处理和地下探测是目前尚待解决的问题。
技术实现思路
本专利技术的目的是提供了一种两道多模式瑞利波地下探测系统及方法。本专利技术基于两道多分量记录实现基阶和第一高阶瑞利波频散曲线的提取，并利用其反演横波速度进行地下结构探测。本专利技术提出的方法减少了地震信号采集设备，降低了野外施工难度，同时能够分离瑞利波基阶和高阶模式，获取多模式频散曲线，提高了探测的精度。本专利技术是通过以下技术方案实现的：本专利技术涉及一种两道多模式瑞利波地下探测系统，包括：人工锤击震源、多分量检波器、采集站主机、电脑以及若干数据传输线；其中，所述多分量检波器的信号输出端通过数据传输线与所述采集站主机的输入端相连；所述采集站主机的信号输出端与所述电脑的信号输入端相连，信号从所述多分量检波器传入所述电脑；所述电脑对垂向和水平径向分量地震数据进行处理提取多模式瑞利波频散曲线并反演得到地下结构。优选地，所述多分量检波器用于数据采集，记录垂向和水平径向分量的地震记录；所述多分量检波器的数量为2个。优选地，所述人工锤击震源与两个多分量检波器布设于一条直测线上，多分量检波器的垂向分量沿铅垂方向竖直向下，水平分量与直测线远离震源的延展方向一致。优选地，人工锤击震源与较近的多分量检波器的间距为5-15米，2个多分量检波器的间距为5-20米，人工锤击震源激发时重锤垂向下落高度不低于2米。优选地，所述电脑安装有fortran或matlab数据处理软件。本专利技术还涉及基于前述两道多模式瑞利波地下探测系统的方法，包括如下步骤：步骤1：以常规一维直测线方式布设2道多分量检波器，并令检波器垂向分量沿重力方向竖直向下，检波器的一个水平向分量沿测线展布方向，利用人工锤击震源在测线端点外侧延长线一定距离处激发地震信号，则两道多分量检波器记录的信号为Uv(i,t)、Uh(i,t)，其中i＝1,2表示道号，t表示时间序列，Uv和Uh分别表示垂向记录和沿着测线方向的水平径向记录；步骤2：对检波器记录的数据进行初至切除、带通滤波预处理，将预处理后的垂向和水平径向记录组合记为复数，其中水平径向地震记录做为复数的实部，垂向地震记录做为复数的虚部，即Ucomplex(i,t)＝Uh(i,t)+j·Uv(i,t)，j表示单位虚数，并规定径向分量远离震源方向的振动位移记为正实部，反之为负实部；垂向分量向下的振动位移记为正虚部，反之为负虚部；步骤3：对步骤2得到的Ucomplex(i,t)计算相位角：式中Imag表示虚部，Real表示实部。将计算得到的相位角解缠，使其沿时间序列连续展开，则质点振动轨迹“逆进”的基阶模式瑞利波信号将呈现出相位角随时间减小的特征，而质点轨迹“顺进”的第一高阶模式瑞利波将呈现相位角随时间增大的特征；步骤4：利用Savitzky-Golay滤波器对步骤3得到的相位角数据进行平滑滤波，消除噪声引起的相位角数据高频扰动，然后绘制相位角随时间变化的曲线图，则基阶模式信号表现为负斜率、第一高阶模式表现为正斜率；步骤5：瑞利波模式分离，根据步骤4得到的两道记录的相位角数据，将斜率为正的部分所对应时刻的时间域波形记录进行充零处理，则得到基阶模式瑞利波记录，记为Uv基阶(i,t)、Uh基阶(i,t)；将斜率为负的部分所对应时刻的时间域波形记录进行充零处理，则得到第一高阶模式瑞利波记录记为Uv高阶(i,t)、Uh高阶(i,t)；步骤6：利用SASW或双台法对基阶模式瑞利波Uv基阶(1,t)和Uv基阶(2,t)进行处理计算得到基阶频散曲线；同样方法对高阶模式瑞利波Uv高阶(1,t)和Uv高阶(2,t)数据进行处理计算得到第一高阶模式频散曲线；步骤7：利用半波长法估计地下地层的初始横波速度模型，然后利用基于奇异值分解的阻尼最小二乘方法对基阶和第一高阶频散曲线进行联合反演，得到地下横波速度结构，实现对地下介质结构的探测。优选地，所述步骤2对检波器记录的数据进行预处理，预处理包括切除初至及带通滤波，滤波的通带频率设为5-60Hz，该频带内的瑞利波基阶和第一高阶模式质点振动轨迹分别表现为逆进和顺进的特征；再对预处理后的多分量检波器的垂向和水平径向记录数据进行组合，构建成实部为径向位移、虚部为垂向位移的复数形式的地震记录数据。优选地，所述步骤3为对步骤2所得的复数地震记录求相位角，并对求得的角度进行相位解缠(unwrap)，使其在时间序列上连续分布，避免反正切函数计算结果的2π周期间断跳跃现象。优选地，所述步骤4为利用Savitzky-Golay滤波器对相位角数据进行平滑滤波，设置滤波器阶数为3阶，若相位角数据点个数记为n，则设置滤波器框长度为大于阶数且处于范围内的奇数。本专利技术具有以下优点：(1)本专利技术提供的两道多模式瑞利波探测方法相比于常规两道探测方法，提高了信号处理的准确性和探测的精度，通过多分量检波器记录垂向和水平径向的地震信号，并利用多分量信号分析瑞利波不同模式的质点振动轨迹特征及其相位角曲线斜率差异，从而实现对基阶模式和第一高阶模式瑞利波信号的分离，最终在两个模式互不干扰的情况下提取了本文档来自技高网...

【技术保护点】
1.一种两道多模式瑞利波地下探测系统，其特征在于，包括：人工锤击震源、多分量检波器、采集站主机、电脑以及若干数据传输线；/n其中，/n所述多分量检波器的信号输出端通过数据传输线与所述采集站主机的输入端相连；/n所述采集站主机的信号输出端与所述电脑的信号输入端相连，信号从所述多分量检波器传入所述电脑；/n所述电脑对垂向和水平径向分量地震数据进行处理提取多模式瑞利波频散曲线并反演得到地下结构。/n

【技术特征摘要】
1.一种两道多模式瑞利波地下探测系统，其特征在于，包括：人工锤击震源、多分量检波器、采集站主机、电脑以及若干数据传输线；
其中，
所述多分量检波器的信号输出端通过数据传输线与所述采集站主机的输入端相连；
所述采集站主机的信号输出端与所述电脑的信号输入端相连，信号从所述多分量检波器传入所述电脑；
所述电脑对垂向和水平径向分量地震数据进行处理提取多模式瑞利波频散曲线并反演得到地下结构。


2.如权利要求1所述的两道多模式瑞利波地下探测系统，其特征在于，所述多分量检波器用于数据采集，记录垂向和水平径向分量的地震数据；所述多分量检波器的数量为2个。


3.如权利要求1所述的两道多模式瑞利波地下探测系统，其特征在于，所述人工锤击震源与两个多分量检波器布设于一条直测线上，多分量检波器的垂向分量沿铅垂方向竖直向下，水平分量与直测线远离震源的延展方向一致。


4.如权利要求3所述的两道多模式瑞利波地下探测系统，其特征在于，人工锤击震源与较近的多分量检波器的间距为5-15米，2个多分量检波器的间距为5-20米，人工锤击震源激发时重锤垂向下落高度不低于2米。


5.如权利要求1所述的两道多模式瑞利波地下探测系统，其特征在于，所述电脑安装有fortran或matlab数据处理软件。


6.一种基于权利要求1所述的两道多模式瑞利波地下探测的方法，其特征在于，包括如下步骤：
步骤1：以常规一维直测线方式布设2道多分量检波器，并令检波器垂向分量沿重力方向竖直向下，检波器的一个水平向分量沿测线展布方向，利用人工锤击震源在测线端点外侧延长线一定距离处激发地震信号，则两道多分量检波器记录的信号为Uv(i,t)、Uh(i,t)，其中i＝1,2表示道号，t表示时间序列，Uv和Uh分别表示垂向记录和沿着测线方向的水平径向记录；
步骤2：对检波器记录的数据进行初至切除、带通滤波预处理，将预处理后的垂向和水平径向记录组合记为复数，其中水平径向记录数据做为复数的实部，垂向记录数据做为复数的虚部，即Ucomplex(i,t)＝Uh(i,t)+j·Uv(i,t)，j表示单位虚数，并规定径向分量远离震源方向的振动位移记为正实部，反之为负实部；垂向分量向下的振动位移记为正虚部，反之为负虚部；
步骤3：对步骤2得到的Ucomplex(i,t)计算相位角：



式中Imag表示虚...

【专利技术属性】
技术研发人员：李欣欣，沈鸿雁，李萌，
申请(专利权)人：西安石油大学，
类型：发明
国别省市：陕西;61