一种地震和直流电法巷道超前探测的数据融合方法与流程

本发明涉及工程地球物理勘探技术领域，更具体的说是涉及一种地震反射波法与直流电阻率法超前探测的数据融合方法。

背景技术：

在巷道掘进过程中，探明前方的地质条件是十分必要的，这也是保证安全掘进的前提。长期以来，巷道掌子面前方的构造异常、采空区、岩溶空洞、陷落柱等不良地质体，是困扰巷道安全掘进的几大难题，如何更准确的探测前方不良地质体，是现在技术研究的热点。

目前巷道探测多使用地震反射波法、直流电阻率法、瞬变电磁法等物探方法，分别对单一物性数据进行处理，然后进行结果的相互验证。考虑到巷道有限的空间，单独使用上述某一种方法探测时，探测精度相对较低。

因此，如何提供一种高精度巷道超前探测的方法是本领域技术人员亟需解决的问题。

技术实现要素：

有鉴于此，本发明提供了一种地震和直流电法巷道超前探测的数据融合方法，基于交叉梯度约束对地震反射波法与直流电阻率法超前探测的数据进行融合计算，提高了超前探测的探测精度。

为了实现上述目的，本发明采用如下技术方案：

一种地震和直流电法巷道超前探测的数据融合方法，具体包括以下步骤：

(1)首先进行地震反射波法速度的单一反演，得到速度数据集dms；同时进行直流电阻率法电阻率的单一反演，得到电阻率数据集dmr；

(2)选取地震反射波法和直流电阻率法探测的公共区域为数据融合区域，得到数据融合区域内的速度数据集d'ms和电阻率数据集d'mr；

(3)对数据融合区域内的速度数据集d'ms和电阻率数据集d'mr进行数据的归一化处理；

(4)对归一化处理后的数据集计算交叉梯度值t(x,z)，对交叉梯度值t(x,z)进行离散化，得到离散化后的数据ti,j；

(5)基于离散化后的数据ti,j和步骤(2)中数据融合区域内的速度数据集d'ms计算速度更新量v″i,j；基于离散化后的数据ti,j和步骤(2)中数据融合区域内的电阻率数据集d'mr计算电阻率更新量r″i,j；

(6)基于速度更新量v″i,j构成新的速度数据集d″ms，并基于新的速度数据集d″ms进行速度成图；基于电阻率更新量r″i,j构成新的电阻率数据集d″mr，并基于新的电阻率数据集d″mr进行电阻率成图。

优选的，所述步骤(3)中，利用min-max标准法对速度数据集d'ms和电阻率数据集d'mr进行数据归一化处理，具体包括：

设a为速度数据集或者电阻率数据集，mina和maxa分别为速度数据集或电阻率数据集的最大值和最小值，将a中每个原始数据x通过min-max标准化映射成在区间[0,1]中的新数据x’，其公式为：新数据x’＝(原始数据x-最小值mina)/(最大值maxa-最小值mina)，所得新数据构成归一化后的速度数据集和电阻率数据集。

优选的，所述步骤(4)具体包括如下步骤：

a、基于二维速度模型与电阻率模型，定义两者之间的交叉梯度值为式中为梯度运算符，ms(x,z)表示地震模型速度数据，mr(x,z)表示电法模型电阻率数据，t(x,z)为定义的交叉梯度值；

b、采用中心差分法，对a中交叉梯度值t(x,z)进行数据离散化，得到：

式中i,j＝1,2,...,n

其中，ms(i,j)表示当前网格节点速度数值，ms(i,j+1)表示当前网格节点沿垂直方向的网格节点数加1位置的速度数值，ms(i,j-1)表示当前网格节点沿垂直方向的网格节点数减1位置的速度数值，mr(i,j)表示当前网格节点电阻率数值，mr(i+1,j)表示当前网格节点沿水平方向的网格节点数加1位置的电阻率数值，mr(i-1,j)表示当前网格节点沿水平方向的网格节点数减1位置的电阻率数值；δxj表示沿水平方向每次计算时网格节点序号相差的个数，即当前网格节点左右相邻的网格节点数值的差值，δzi表示沿垂直方向每次计算时网格节点序号相差的个数，即当前网格节点上下相邻的网格节点数值的差值。

优选的，所述步骤(4)中，速度更新量其中i,j＝1,2,3……n；式中，αs为地震参数的权重因子，表示数据融合区域的速度集数据；电阻率更新量其中i,j＝1,2,3……n；式中，βr为电阻率参数的权重因子，表示数据融合区域的电阻率集数据。

优选的，在所述步骤(1)之前还包括：初始化速度模型和初始化电阻率模型。

经由上述的技术方案可知，与现有技术相比，本发明公开提供了一种地震和直流电法巷道超前探测的数据融合方法，首先进行地震反射波法速度的单一反演和直流电阻率法电阻率的单一反演，得到速度数据集和电阻率数据集；选择地震反射波法和直流电阻率法探测的公共区域为数据融合区域，对数据融合区域内的速度数据集和电阻率数据集的数据进行归一化处理；对归一化处理后的数据集计算交叉梯度值并进行离散化；基于离散化的交叉梯度值和数据融合区域内的速度数据集和电阻率数据集计算速度更新量和电阻率更新量；基于速度更新量和电阻率更新量得到新的速度数据集和电阻率数据集，并基于新的速度数据集进行速度成像，基于新的电阻率数据集进行电阻率成像。本发明的技术方案基于交叉梯度约束对地震反射波法与直流电阻率法超前探测的数据进行融合计算，提高了超前探测的探测精度。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。

图1为本发明提供的地震和直流电法巷道超前探测的数据融合方法的流程图。

图2为本发明提供的交叉梯度离散化网格示意图。

图3为本发明提供的数据融合区域示意图。

图4a为本发明提供的模拟数据单一反演速度剖面图。

图4b为本发明提供的模拟数据数据融合后速度剖面图。

图5a为本发明提供的模拟数据单一反演电阻率剖面图。

图5b为本发明模拟数据数据融合后电阻率剖面图。

其中，图3中的1表示巷道掌子面；2表示巷道；3表示电法探测区域；4表示地震探测区域；5表示数据融合区域。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

参见图1，本发明实施例公开了一种地震和直流电法巷道超前探测的数据融合方法，具体步骤如下：

(1)参见图4a和图5a，首先进行地震反射波法速度的单一反演，得到速度数据集dms；同时进行直流电阻率法电阻率的单一反演，得到电阻率数据集dmr；

(2)参见图3，选取地震反射波法和直流电阻率法探测的公共区域为数据融合区域，得到数据融合区域内的速度数据集d'ms和电阻率数据集d'mr；

(3)对数据融合区域内的速度数据集d'ms和电阻率数据集d'mr进行数据的归一化处理；

(4)参见图2，对归一化处理后的数据集计算交叉梯度值t(x,z)，对交叉梯度值t(x,z)进行离散化，得到离散化后的数据ti,j；

(5)基于离散化后的数据ti,j和步骤(2)中数据融合区域内的速度数据集d'ms计算速度更新量v″i,j；基于离散化后的数据ti,j和步骤(2)中数据融合区域内的电阻率数据集d'mr计算电阻率更新量r″i,j；

(6)参见图4b和图5b，基于速度更新量v″i,j构成新的速度数据集d″ms，并基于新的速度数据集d″ms进行速度成图；基于电阻率更新量r″i,j构成新的电阻率数据集d″mr，并基于新的电阻率数据集d″mr进行电阻率成图。

优选的，步骤(3)中，利用min-max标准法对速度数据集d'ms和电阻率数据集d'mr进行数据归一化处理，具体包括：

设a为速度数据集或者电阻率数据集，mina和maxa分别为速度数据集或电阻率数据集的最大值和最小值，将a中每个原始数据x通过min-max标准化映射成在区间[0,1]中的新数据x’，其公式为：新数据x’＝(原始数据x-最小值mina)/(最大值maxa-最小值mina)，所得新数据构成归一化后的速度数据集和电阻率数据集。

优选的，参见图2，步骤(4)具体包括如下步骤：

a、基于二维速度模型与电阻率模型，定义两者之间的交叉梯度值为式中为梯度运算符，ms(x,z)表示地震模型速度数据，mr(x,z)表示电法模型电阻率数据，t(x,z)为定义的交叉梯度值；

b、采用中心差分法，对a中交叉梯度值t(x,z)进行数据离散化，得到：

式中i,j＝1,2,...,n

其中，ms(i,j)表示当前网格节点速度数值，ms(i,j+1)表示当前网格节点沿垂直方向的网格节点数加1位置的速度数值，ms(i,j-1)表示当前网格节点沿垂直方向的网格节点数减1位置的速度数值，mr(i,j)表示当前网格节点电阻率数值，mr(i+1,j)表示当前网格节点沿水平方向的网格节点数加1位置的电阻率数值，mr(i-1,j)表示当前网格节点沿水平方向的网格节点数减1位置的电阻率数值；δxj表示沿水平方向每次计算时网格节点序号相差的个数，即当前网格节点左右相邻的网格节点数值的差值，δzi表示沿垂直方向每次计算时网格节点序号相差的个数，即当前网格节点上下相邻的网格节点数值的差值。

优选的，步骤(4)中，速度更新量其中i,j＝1,2,3……n；式中，αs为地震参数的权重因子，表示数据融合区域的速度集数据；电阻率更新量其中i,j＝1,2,3……n；式中，βr为电阻率参数的权重因子，表示数据融合区域的电阻率集数据。

优选的，在步骤(1)之前还包括：初始化速度模型和初始化电阻率模型。

经由上述的实施例可知，与现有技术相比，本发明公开提供了一种地震和直流电法巷道超前探测的数据融合方法，首先进行地震反射波法速度的单一反演和直流电阻率法电阻率的单一反演，得到速度数据集和电阻率数据集；选择地震反射波法和直流电阻率法探测的公共区域为数据融合区域，对数据融合区域内的速度数据集和电阻率数据集的数据进行归一化处理；对归一化处理后的数据集计算交叉梯度值并进行离散化；基于离散化的交叉梯度值和数据融合区域内的速度数据集和电阻率数据集计算速度更新量和电阻率更新量；基于速度更新量和电阻率更新量得到新的速度数据集和电阻率数据集，并基于新的速度数据集进行速度成像，基于新的电阻率数据集进行电阻率成像。本发明的实施例基于交叉梯度约束对地震反射波法与直流电阻率法超前探测的数据进行融合计算，提高了超前探测的探测精度。

本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。

对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。