CN111290017B

CN111290017B - 一种震电波场联合提取瑞雷波频散特征的面波勘探方法

Info

Publication number: CN111290017B
Application number: CN202010142633.1A
Authority: CN
Inventors: 任恒鑫; 杨振涛; 陈晓非; 袁士川; 郑旭桢; 何展翔
Original assignee: Southern University of Science and Technology
Current assignee: Southern University of Science and Technology
Priority date: 2020-03-04
Filing date: 2020-03-04
Publication date: 2020-12-15
Anticipated expiration: 2040-03-04
Also published as: CN111290017A

Abstract

本发明公开了一种震电波场联合提取瑞雷波频散特征的面波勘探方法，包括步骤：获取联合采集的数据；其中，所述联合采集的数据包括：地震波数据和电场数据；将所述联合采集的数据进行联合成像处理得到叠加的频散谱；对所述叠加的频散谱进行提取处理得到频散曲线，并对所述频散曲线进行反演处理得到地层结构剖面。由于采用地震波数据和电场数据进行联合成像处理得到叠加的频散谱，并提取得到多模式频散曲线，在进行反演时大大降低反演的多解性，从而大大提高面波勘探的精度和稳定性。

Description

一种震电波场联合提取瑞雷波频散特征的面波勘探方法

技术领域

本发明涉及地质、地球物理勘探技术领域，尤其涉及的是一种震电波场联合提取瑞雷波频散特征的面波勘探方法。

背景技术

瑞雷波(Rayleigh wave)由英国科学家Lord Rayleigh首先于1887年在理论上确定，故此命名。自二十世纪五十年代，随着科学家研究不断深入，发现瑞雷波(Rayleighwave)在层状介质中传播时携带了各层介质的P波速度、S波速度、密度等介质信息，呈现明显的频散特性(即速度随频率变化而变化)，且主要取决于层状介质的S波速度结构。所以瑞雷波在传播过程中能量和速度的变化特征携带了大量地层信息，因此在工程领域常常通过研究Rayleigh波频散特征来解决工程勘察、场地和地基处理评价、障碍物和空洞探测等浅层地质问题。现阶段面波勘探已成为工程物探领域应用最广泛的物探方法之一。

在实际应用中，大量学者发现面波勘探如果仅采用基阶模式面波频散信息进行反演得到的地层模型具有很大的不确定性，但将基阶、高阶模式频散联合反演将会大大降低这种不确定性；而且在工程勘探中，经常遇到某一频率范围内的瑞雷波高阶模式比基阶模式具有更强的能量，这就意味着在某些频率范围内我们只能得到高阶模式频散曲线。现阶段面波勘探方法均是从面波的波动信息中提取频散信息，虽然基阶模式能较好的成像，但高阶尤其较高模式的高阶模式频散成像精度有限。

因此，现有技术还有待于改进和发展。

发明内容

本发明要解决的技术问题在于，针对现有技术的上述缺陷，提供一种震电波场联合提取瑞雷波频散特征的面波勘探方法，旨在解决现有技术中高阶模式频散成像精度低的问题。

本发明解决技术问题所采用的技术方案如下：

一种震电波场联合提取瑞雷波频散特征的面波勘探方法，其中，包括步骤：

获取联合采集的数据；其中，所述联合采集的数据包括：地震波数据和电场数据；

将所述联合采集的数据进行联合成像处理得到叠加的频散谱；

对所述叠加的频散谱进行提取处理得到频散曲线，并对所述频散曲线进行反演处理得到地层结构剖面。

所述的震电波场联合提取瑞雷波频散特征的面波勘探方法，其中，所述地震波数据包括：地震波分量；所述电场数据包括：电场分量；

所述将所述联合采集的数据进行联合成像处理得到叠加的频散谱，包括：

根据所述地震波分量、所述电场分量分别得到地震波分量的频散谱、电场分量的频散谱；

将所述地震波分量的频散谱和所述电场分量的频散谱叠加得到叠加的频散谱。

所述的震电波场联合提取瑞雷波频散特征的面波勘探方法，其中，所述根据所述地震波分量、所述电场分量分别得到地震波分量的频散谱、电场分量的频散谱，包括：

分别提取所述地震波分量的虚部、所述电场分量的虚部后进行频率扫描得到地震波分量的频散谱、电场分量的频散谱。

所述的震电波场联合提取瑞雷波频散特征的面波勘探方法，其中，所述地震波分量包括：地震波径向分量和/或地震波垂向分量；所述电场分量包括：电场径向分量和/或电场垂向分量。

所述的震电波场联合提取瑞雷波频散特征的面波勘探方法，其中，所述将所述地震波分量的频散谱和所述电场分量的频散谱叠加得到叠加的频散谱，包括：

将地震波垂向分量的频散谱和电场径向分量的频散谱叠加得到叠加的频散谱；或者

将地震波径向分量的频散谱和电场径向分量的频散谱叠加得到叠加的频散谱；或者

将地震波垂向分量的频散谱、地震波径向分量的频散谱和电场径向分量的频散谱叠加得到叠加的频散谱。

所述的震电波场联合提取瑞雷波频散特征的面波勘探方法，其中，所述对所述叠加的频散谱进行提取处理得到频散曲线，包括：

通过能量极大值来提取所述叠加的频散谱得到频散曲线；所述频散曲线为含有基阶模式、高阶模式的多模式频散曲线。

所述的震电波场联合提取瑞雷波频散特征的面波勘探方法，其中，所述对所述频散曲线进行反演处理得到地层结构剖面，包括：

建立初始反演地层模型，利用多种反演方法对所述频散曲线进行拟合得到地层结构剖面；所述反演方法包括：遗传算法、群类算法、拟牛顿算法中的一种或多种。

一种震电波场联合提取瑞雷波频散特征的面波勘探装置，其中，包括：

采集装置，用于获取联合采集的数据；其中，所述联合采集的数据包括：地震波数据和电场数据；

联合成像装置，用于将所述联合采集的数据进行联合成像处理得到叠加的频散谱；

反演装置，用于对所述叠加的频散谱进行提取处理得到频散曲线，并对所述频散曲线进行反演处理得到地层结构剖面。

一种终端设备，其中，包括：处理器，以及与所述处理器连接的存储器，

所述存储器存储有震电波场联合提取瑞雷波频散特征的面波勘探程序，所述震电波场联合提取瑞雷波频散特征的面波勘探程序被所述处理器执行时实现以下步骤：

一种存储介质，其中，其上存储有震电波场联合提取瑞雷波频散特征的面波勘探程序，所述震电波场联合提取瑞雷波频散特征的面波勘探程序被处理器执行时实现以下步骤：

有益效果：由于采用地震波数据和电场数据进行联合成像处理得到叠加的频散谱，并提取得到多模式频散曲线，在进行反演时大大降低反演的多解性，从而大大提高面波勘探的精度和稳定性。

附图说明

图1是本发明中线性观测系统接收到地震波分量u_r、u_z以及电场分量E_r、E_z的炮集记录。

图2是本发明中地震波分量u_r、u_z以及电场分量E_r、E_z的炮集记录得到的频散谱。

图3a是本发明中由地震波分量u_r、u_z的频散谱得到的叠加的频散谱D_Seismic。

图3b是本发明中由径向电场E_r的频散谱与地震波分量u_r、u_z的频散谱得到的叠加的频散谱D_{Seismoelectirc}。

图4a是本发明中u_z、E_r叠加得到的叠加的频散谱。

图4b是本发明中u_z、u_r、E_r叠加得到的叠加的频散谱。

图5是本发明中频散曲线的示意图。

图6是本发明中多模式频散曲线进行反演得到的地层速度结构。

图7是本发明中将测线中多个面波反演的速度结构插值绘制成地层横波速度结构剖面。

图8是本发明中震、电联合采集示意图。

图9是本发明中震电波场联合提取瑞雷波频散特征的面波勘探方法的流程图。

具体实施方式

为使本发明的目的、技术方案及优点更加清楚、明确，以下参照附图并举实施例对本发明进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

请同时参阅图1-图9，本发明提供了一种震电波场联合提取瑞雷波频散特征的面波勘探方法的一些实施例。

震电效应最早由前苏联科学家Ivanov在1939年发现，在不加电压的情况下，测量到了由地震波导致的电场，并指出这可能与双电层动电效应有关。此后，Frenkel建立了固相运动与孔隙流体运动相耦合的孔隙介质波动理论，并得到了Biot的进一步完善，形成了著名的Biot孔隙介质弹性波理论。上世纪90年代，随着电子科学技术的发展以及微弱信号提取和信号处理手段的提高，震电效应的研究得到了长足发展。1994年Pride在Frenkel和Biot的工作基础上提出了震电耦合的宏观控制方程组，该方程组中将Biot空隙弹性方程组与Maxwell电磁方程组通过动电耦合系数耦合在一起，描述了在多孔饱和介质中波动场与电磁场的耦合关系，现今已成为震电理论研究的基础。

在层状孔隙介质模型的波场模拟研究中发现，当震波超过临界角入射孔隙介质界面时，震电转换产生的电磁波振幅在离开界面时快速衰减，称其为隐失电磁波。研究发现垂向单力点源所激发的瑞雷面波产生的隐失电磁波(EM)具有与瑞雷波一致的频散特性。这样为通过接收面波产生的隐失电磁波信号来提取面波的频散信息提供了理论基础。研究表明，通过隐失电磁波计算得到的频散谱，包含有丰富的、高质量的高阶模式信息，并且这些高阶模式频散是直接利用地震波(瑞雷波)计算得到的频散谱无法提供的，另外对于基阶模式，隐失电磁波和地震波在频散谱中具有显著能量的频率范围不尽相同，两者正好形成互补。所以，震电转换产生的隐失电磁波(EM)具有提高瑞雷波频散成像质量的巨大潜力，这对于进一步提高瑞雷面波勘探精度具有重大意义。

如图9所示，本发明的一种震电波场联合提取瑞雷波频散特征的面波勘探方法，包括以下步骤：

S100、获取联合采集的数据；其中，所述联合采集的数据包括：地震波数据和电场数据。

具体地，所述联合采集的数据还可以包括：磁场数据。所述地震波数据包括：地震波分量；所述电场数据包括：电场分量。更具体地，所述地震波分量包括：地震波径向分量u_r和/或地震波垂向分量u_z；所述电场分量包括：电场径向分量E_r和/或电场垂向分量E_z。

实际工程应用中瞬态瑞雷波勘探所采用的震源一般为锤击、落锤等震源，可近似为一垂直方向的单点力源，理论上仅能产生纵波(P波)和垂直分量横波(SV波)以及横磁模式电磁场。激发的P波和SV波相干产生瑞雷(Rayleigh)波并产生相应的隐失电磁波。所以在实际勘探中，我们能观测到瑞雷波的径向分量u_r和垂向分量u_z(柱坐标系)，以及电场的径向分量E_r和垂向分量E_z，以及磁场的横向分量B_θ，由于磁场分量微弱，故本实施例中仅考虑地震波分量和电场分量。

由于电场分量信号较地震波信号微弱，故采集装置采用不同的降噪和前段放大系统。当锤击激发震源后，震、电两种不同的传感器同时进行高频数字采样。

S200、将所述联合采集的数据进行联合成像处理得到叠加的频散谱。

具体地，步骤S200包括步骤：

S210、根据所述地震波分量、所述电场分量分别得到地震波分量的频散谱、电场分量的频散谱。

S210包括步骤：

S211、分别提取所述地震波分量的虚部、所述电场分量的虚部后进行频率扫描得到地震波分量的频散谱、电场分量的频散谱。

在水平层状地层模型中，脉冲震源时间函数的响应波场，例如G(r,θ,z,ω)，能够表示为：

其中，k′为所有可能的水平波数，

为一组完备正交的矢量基函数，具体表示如下：

其中，

i为虚数单位，m＝0,±1,±2,...，k′∈(0,+∞)；J_m(k′r)表示m阶第一类贝塞尔函数。式(1)中G_T,m，G_S,m和G_R,m是与基函数相对应的展开系数，具体表示如下：

其中，符号*表示复共轭。在柱坐标系中，G(r,θ,z,ω)也能表示为：

G(r,θ,z,ω)＝G_r(r,θ,z,ω)e_r+G_θ(r,θ,z,ω)e_θ+G_z(r,θ,z,ω)e_z. (4)

结合(1)、(2)、(3)、(4)，并考虑0阶第一类贝塞尔函数的导数性质：J′₀(k′r)＝-J₁(k′r)，我们可以得到：

根据贝塞尔函数的正交性质：

(其中k为观测的水平波数)，由式(5)可得：

定义如下的核函数：

将式(9)代入式(6)-(8)可得：

因为脉冲震源时间函数的响应波场G(r,θ,z,ω)与实际震源时间函数F(ω)以及实际波场U(r,θ,z,ω)之间存在如下关系：

U(r,θ,z,ω)＝F(ω)G(r,θ,z,ω)， (11)

所以，地震波的径向分量u_r和垂向分量u_z(采用柱坐标系)，以及电场的径向分量E_r和垂向分量E_z，根据式(10)、(11)可以表示为：

其中u_S(k,ω,z)，u_R(k,ω,z)，E_S(k,ω,z)，E_R(k,ω,z)为波动和电场的核函数；F(ω)表示实际震源时间函数的傅里叶变换；u_r(r,ω,z)，u_z(r,ω,z)，E_r(r,ω,z)，E_z(r,ω,z)为时间域波场分量u_r(r,t,z)，u_z(r,t,z)，E_r(r,t,z)，E_z(r,t,z)的傅里叶变换；r和z在实际观测中分别表示一个接收器与震源的径向距离以及接收器的深度。

这样，在面波勘探中，可将(12)中的积分近似为求和，并取其虚部进行频率扫描从而得到频散谱，具体公式如下：

其中，

和

表示近似得到的频散谱；Im[·]表示求取一个复变量的虚部；k为观测的水平波数，ω为角频率，N表示记录道的总数；r_j表示第j个检波器到震源的径向距离，由于面波勘探基于水平层状地层，所有传感器与震源无需要多个传感器沿线性等间距排列或其他规则形状排列，也无需与震源保持一定的偏移距，故检波器可根据实际情况任意布设。也就是说，本实施例中，通过地震波、电场数据的联合采集和成像(现有面波勘探仅用地震波分量进行探测法)，采集系统可根据探测区域任意布设，无需线性或其他规则形状布设，且无需偏移距。所以可根据探测场地环境按任意形状排列，这样能更加综合地反映探测区域的地下速度结构，同时相较于传统的直线等间距排列且需要一定最小偏移距的面波观测系统具有更高的适用性。

S220、将所述地震波分量的频散谱和所述电场分量的频散谱叠加得到叠加的频散谱。

地球浅地表介质中时常存在地下水位面，在水位面以上介质，尤其是对于接近地表的介质，可以近似为空气填充的孔隙介质。因此，我们根据实际地层情况建立一个四层孔隙半空间模型，参数详见表1。

表1.一个四层含高速夹层半空间孔隙介质模型的参数

该模型第二层被设置为饱含水的孔隙介质，其他层都被设置为饱含空气的孔隙介质，潜水面位于地下10米深度的位置。震源为地表激发的垂向单力点源(锤击、落锤震源)，震源时间函数采用主频为10赫兹、延迟时间为0.5秒的雷克子波(Ricker wavelet)；为了便于数据采集，如图8所示，观测系统采用线性排列，道间距2米，共91道，接收到的地震波分量炮集记录和电场分量炮集记录见图1；根据地震波各分量炮集记录、电场各分量炮集记录得到的频散谱见图2。图2中点线为由地层参数计算得到的瑞雷波理论频散曲线。

S220包括步骤：

S221、将地震波垂向分量的频散谱和电场径向分量的频散谱叠加得到叠加的频散谱；或者

S222、将地震波径向分量的频散谱和电场径向分量的频散谱叠加得到叠加的频散谱；或者

S223、将地震波垂向分量的频散谱、地震波径向分量的频散谱和电场径向分量的频散谱叠加得到叠加的频散谱。

从图1，我们可以看到炮集记录上瑞雷波占据主导地位，并且在面波区域电场分量E_r、E_z比地震波分量u_r、u_z高阶模式更加发育。图2频散谱中的黑色点线为理论频散曲线，可以看到频散谱中能量最大值与理论频散曲线位置相吻合，所以我们可通过能量极大值来提取瑞雷波频散曲线。从图2中我们可以看到在电场分量E_r、E_z的频散谱

中，瑞雷波高阶模式的成像质量明显高于地震分量u_r、u_z的频散谱

u_r、u_z的频散谱中基阶模式的成像频段明显宽于E_r、E_z的频散谱。所以地震波分量和电场分量联合提取多模式频散曲线进行反演将大大降低反演的多解性，从而大大提高反演的精度和稳定性。

从图3a，我们可以看到仅用地震波u_r、u_z分量叠加得到的频散谱D_Seismic中高阶模式成像依然没有电场分量得到的高阶模式多。考虑到电场垂直分量E_z信号较弱并且在实际探测中数据采集较难，所以实际应用中我们将电场径向分量E_r的频散谱与D_Seismic叠加得到新的频散谱D_{Seismoelectirc}(如图3b所示)，可以看到叠加后的频散谱中，基阶模式和高阶模式的成像质量均得到了提高。

地震波径向分量u_r的采集需要水平两个分量进行转换，实际探测中需要三分量检波器进行采集，而地震波垂向分量u_z的采集仅需要垂直检波器，不但经济而且高效，所以我们将u_z、E_r频散谱叠加与u_z、u_r、E_r三个分量叠加频散谱进行比较。图4a为u_z、E_r叠加得到的频散谱，图4b为u_z、u_r、E_r叠加得到的频散谱，我们可以看到两者没有显著差别，所以也可以只用u_z、E_r进行叠加，从而避免多个地震波、电场分量的采集和叠加。这样可以在保证频散谱成像质量的前提下，有效地降低采集成本。

S300、对所述叠加的频散谱进行提取处理得到频散曲线，并对所述频散曲线进行反演处理得到地层结构剖面。

具体地，S300包括步骤：

S310、通过能量极大值来提取所述叠加的频散谱得到频散曲线。

具体地，所述频散曲线为含有基阶模式、高阶模式的多模式频散曲线。通过设置区域，然后通过计算机程序自动搜索区域内的极值点的方法从频散谱上提取含有基阶、高阶模式的多模式频散曲线(见图5)，以便用于后期反演。

S320、建立初始反演地层模型，利用多种反演方法对所述频散曲线进行拟合得到地层结构剖面。

具体地，所述反演方法包括：遗传算法、群类算法、拟牛顿算法中的一种或多种。根据基础资料建立初始反演地层模型，然后可利用遗传算法、群类算法(粒子群、蜂群等)、拟牛顿等多种反演方法对图5中提取得到的多模式频散曲线进行拟合，从而反演地层结构，见图6。

将观测系统移动就可以得到多个地层位置的速度结构，将不同位置的速度结构进行插值，根据空间位置即可绘制地层剖面，如图7。

基于上述任意一实施例所述的震电波场联合提取瑞雷波频散特征的面波勘探方法，本发明还提供了一种震电波场联合提取瑞雷波频散特征的面波勘探装置的较佳实施例：

本发明实施例所述一种震电波场联合提取瑞雷波频散特征的面波勘探装置，包括：采集装置，用于获取联合采集的数据；其中，所述联合采集的数据包括：地震波数据和电场数据；

所述地震波数据包括：地震波分量；所述电场数据包括：电场分量。

所述联合成像装置具体用于根据所述地震波分量、所述电场分量分别得到地震波分量的频散谱、电场分量的频散谱；并将所述地震波分量的频散谱和所述电场分量的频散谱叠加得到叠加的频散谱。

所述联合成像装置还用于分别提取所述地震波分量的虚部、所述电场分量的虚部后进行频率扫描得到地震波分量的频散谱、电场分量的频散谱。

所述地震波分量包括：地震波径向分量和/或地震波垂向分量；所述电场分量包括：电场径向分量和/或电场垂向分量。

所述联合成像装置还用于将地震波垂向分量的频散谱和电场径向分量的频散谱叠加得到叠加的频散谱；或者

所述反演装置具体用于通过能量极大值来提取所述叠加的频散谱得到频散曲线；所述频散曲线为含有基阶模式、高阶模式的多模式频散曲线。

所述反演装置还用于建立初始反演地层模型，利用多种反演方法对所述频散曲线进行拟合得到地层结构剖面；所述反演方法包括：遗传算法、群类算法、拟牛顿算法中的一种或多种。

基于上述任意一实施例所述的震电波场联合提取瑞雷波频散特征的面波勘探方法，本发明还提供了一种终端设备的较佳实施例：

本发明实施例所述一种终端设备，包括：处理器，以及与所述处理器连接的存储器，

所述震电波场联合提取瑞雷波频散特征的面波勘探程序被所述处理器执行时，还实现以下步骤：

将地震波径向分量的频散谱、地震波垂向分量的频散谱和电场径向分量的频散谱叠加得到叠加的频散谱。

通过能量极大值来提取所述叠加的频散谱得到频散曲线。

所述频散曲线为含有基阶模式、高阶模式的多模式频散曲线。

建立初始反演地层模型，利用多种反演方法对所述频散曲线进行拟合得到地层结构剖面。

所述反演方法包括：遗传算法、群类算法、拟牛顿算法中的一种或多种。

基于上述任意一实施例所述的震电波场联合提取瑞雷波频散特征的面波勘探方法，本发明还提供了一种存储介质的较佳实施例：

本发明实施例所述一种存储介质，其上存储有震电波场联合提取瑞雷波频散特征的面波勘探程序，所述震电波场联合提取瑞雷波频散特征的面波勘探程序被处理器执行时实现以下步骤：

应当理解的是，本发明的应用不限于上述的举例，对本领域普通技术人员来说，可以根据上述说明加以改进或变换，所有这些改进和变换都应属于本发明所附权利要求的保护范围。

Claims

1.一种震电波场联合提取瑞雷波频散特征的面波勘探方法，其特征在于，包括步骤：

对所述叠加的频散谱进行提取处理得到频散曲线，并对所述频散曲线进行反演处理得到地层结构剖面；

所述地震波数据包括：地震波分量；所述电场数据包括：电场分量；

2.根据权利要求1所述的震电波场联合提取瑞雷波频散特征的面波勘探方法，其特征在于，所述根据所述地震波分量、所述电场分量分别得到地震波分量的频散谱、电场分量的频散谱，包括：

3.根据权利要求1-2任意一项所述的震电波场联合提取瑞雷波频散特征的面波勘探方法，其特征在于，所述地震波分量包括：地震波径向分量和/或地震波垂向分量；所述电场分量包括：电场径向分量和/或电场垂向分量。

4.根据权利要求3所述的震电波场联合提取瑞雷波频散特征的面波勘探方法，其特征在于，所述将所述地震波分量的频散谱和所述电场分量的频散谱叠加得到叠加的频散谱，包括：

5.根据权利要求1所述的震电波场联合提取瑞雷波频散特征的面波勘探方法，其特征在于，所述对所述叠加的频散谱进行提取处理得到频散曲线，包括：

6.根据权利要求1所述的震电波场联合提取瑞雷波频散特征的面波勘探方法，其特征在于，所述对所述频散曲线进行反演处理得到地层结构剖面，包括：

7.一种震电波场联合提取瑞雷波频散特征的面波勘探装置，其特征在于，包括：

反演装置，用于对所述叠加的频散谱进行提取处理得到频散曲线，并对所述频散曲线进行反演处理得到地层结构剖面；

所述地震波数据包括：地震波分量；所述电场数据包括：电场分量；所述联合成像装置具体用于根据所述地震波分量、所述电场分量分别得到地震波分量的频散谱、电场分量的频散谱；并将所述地震波分量的频散谱和所述电场分量的频散谱叠加得到叠加的频散谱。

8.一种终端设备，其特征在于，包括：处理器，以及与所述处理器连接的存储器，

所述存储器存储有震电波场联合提取瑞雷波频散特征的面波勘探程序，所述震电波场联合提取瑞雷波频散特征的面波勘探程序被所述处理器执行时实现权利要求1-6任意一项所述的震电波场联合提取瑞雷波频散特征的面波勘探方法的步骤。

9.一种存储介质，其特征在于，其上存储有震电波场联合提取瑞雷波频散特征的面波勘探程序，所述震电波场联合提取瑞雷波频散特征的面波勘探程序被处理器执行时实现权利要求1-6中任意一项所述的震电波场联合提取瑞雷波频散特征的面波勘探方法的步骤。