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

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

A data fusion method of seismic and DC electrical method in advance detection of roadway

【技术实现步骤摘要】
一种地震和直流电法巷道超前探测的数据融合方法
本专利技术涉及工程地球物理勘探
，更具体的说是涉及一种地震反射波法与直流电阻率法超前探测的数据融合方法。
技术介绍
在巷道掘进过程中，探明前方的地质条件是十分必要的，这也是保证安全掘进的前提。长期以来，巷道掌子面前方的构造异常、采空区、岩溶空洞、陷落柱等不良地质体，是困扰巷道安全掘进的几大难题，如何更准确的探测前方不良地质体，是现在技术研究的热点。目前巷道探测多使用地震反射波法、直流电阻率法、瞬变电磁法等物探方法，分别对单一物性数据进行处理，然后进行结果的相互验证。考虑到巷道有限的空间，单独使用上述某一种方法探测时，探测精度相对较低。因此，如何提供一种高精度巷道超前探测的方法是本领域技术人员亟需解决的问题。
技术实现思路
有鉴于此，本专利技术提供了一种地震和直流电法巷道超前探测的数据融合方法，基于交叉梯度约束对地震反射波法与直流电阻率法超前探测的数据进行融合计算，提高了超前探测的探测精度。为了实现上述目的，本专利技术采用如下技术方案：一种地震和直流电法巷道超前探测的数据融合方法，具体包括以下步骤：(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进行数据的归一化处本文档来自技高网...

【技术保护点】
1.一种地震和直流电法巷道超前探测的数据融合方法，其特征在于，包括以下步骤：/n(1)首先进行地震反射波法速度的单一反演，得到速度数据集Dm

【技术特征摘要】
1.一种地震和直流电法巷道超前探测的数据融合方法，其特征在于，包括以下步骤：
(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进行电阻率成图。


2.根据权利要求1所述的一种地震和直流电法巷道超前探测的数据融合方法，其特征在于，所述步骤(3)中，利用Min-max标准法对速度数据集D'ms和电阻率数据集D'mr进行数据归一化处理，具体包括：
设A为速度数据集或者电阻率数据集，minA和maxA分别为速度数据集或电阻率数据集的最大值和最小值，将A中每个原始数据x通过Min-max标准化映射成在区间[0,1]中的新数据x’，其公式为：新数据x’＝(原始数据x-最小值minA)/(最大值maxA-最小值minA)，所得新数据构成归一化后的速度数据集和电阻率数...

【专利技术属性】
技术研发人员：张平松，邱实，李圣林，郭立全，胡雄武，
申请(专利权)人：安徽理工大学，
类型：发明
国别省市：安徽;34