CN111856589A - 一种地表综合物探的三维立体成像方法及系统 - Google Patents

Abstract

本发明公开了一种地表综合物探的三维立体成像方法及系统，包括以下步骤：采集地表探测现场的多个二维剖面的探测数据；对探测数据通过地球物理反演形成二维剖面电阻率数据；将二维剖面电阻率数据进行三维坐标转换，获取三维坐标系的电阻率数据；利用克里金插值法将三维坐标系的电阻率数据转化成三维立体模型。

Description

一种地表综合物探的三维立体成像方法及系统

技术领域

本发明属于物探技术领域，具体涉及一种地表综合物探的三维立体成像方法及系统。

背景技术

这里的陈述仅提供与本发明相关的背景技术，而不必然地构成现有技术。

在当今的物探领域，在进行地表探测时，诸如跨孔法、井地法、高密度电法或者是综合物探方法等，都是在工程现场较为常见的探测手段。发明人发现，上述物探方法的探测成果多为二维剖面，现场布置的测线数量众多，形成的二维剖面十分冗杂，得到的单个二维剖面电阻率数据仅仅只能反映一个地质剖面的信息，难以反映地下复杂结构体的实际情况，也同样难以建立能直观展示成果的模型，可视性较差，而在进行实际问题分析解释时往往需要运用到三维模型的建立。

发明内容

针对现有技术存在的不足，本发明的目的是提供一种地表综合物探的三维立体成像方法及系统，该方法能够在完成现场检测之后将现场探测获取到的二维平面的数据经过处理后转换到三维空间坐标系中建立三维电阻率剖面模型，方便后期进行数据分析时三维模型的建立。

为了实现上述目的，本发明是通过如下的技术方案来实现：

第一方面，本发明的实施例提供了一种地表综合物探的三维立体成像方法，包括以下步骤：

采集地表探测现场的多个二维剖面的探测数据；

对探测数据通过地球物理反演形成二维剖面电阻率数据；

将二维剖面电阻率数据进行三维坐标转换，获取三维坐标系的电阻率数据；

利用克里金插值法将三维坐标系的电阻率数据转化成三维立体模型。

作为进一步的技术方案，采集二维剖面的探测数据时，采用跨孔法、井地法、高密度电法、瞬变电磁法、地表综合物探方法中的一种进行地表探测。

作为进一步的技术方案，三维坐标转换的过程为：

确定原点，建立三维坐标轴；

获取测线位置坐标及二维剖面电阻率数据点在测线平面的位置；

通过三维坐标转换公式获得二维剖面电阻率数据点的三维坐标。

作为进一步的技术方案，所述三维坐标轴X方向为水平距离方向，以水平向右为正方向；Y方向为探测深度方向，以平行地面朝探测方向为正方向；Z方向为垂直距离方向，以垂直地面往下为正方向。

作为进一步的技术方案，所述三维坐标转换公式为：

其中X、Y、Z为以O为原点的最终三维坐标，X1为测线起点距离坐标原点的水平距离，Y1为测线起点距离坐标原点的纵向埋深，Z1为测线起点距离坐标原点的垂直高度，X＇为原始数据点的水平长度，Y＇为原始数据点的探测深度，Z＇初值为0，R1为位置矩阵，R2为数据点矩阵。

作为进一步的技术方案，所述R1为

作为进一步的技术方案，若测线平面平行于三维坐标系XOZ平面，则R2为

作为进一步的技术方案，若测线平面平行于三维坐标系YOZ平面，则R2为

作为进一步的技术方案，若测线平面垂直于XOY平面，且与YOZ平面Y轴正方向夹角夹角为α，则R2为

若测线平面垂直于XOY平面，且与YOZ平面Y轴负方向夹角夹角为α，则R2为

第二方面，本发明实施例还提供了一种地表综合物探的三维立体成像系统，包括：

采集模块，用于采集地表探测现场的多个二维剖面的探测数据；

反演模块，用于对探测数据通过地球物理反演形成二维剖面电阻率数据；

坐标转换模块，用于将二维剖面电阻率数据进行三维坐标转换，获取三维坐标系的电阻率数据；

转化模块，用于利用克里金插值法将三维坐标系的电阻率数据转化成三维立体模型。

上述本发明的实施例的有益效果如下：

本发明的三维立体成像方法具有良好的兼容性，无论是跨孔法、井地法、高密度电法、瞬变电磁法或者是综合物探方法,都可以通过该方法将获取的二维电阻率剖面数据转化为三维空间直角坐标系数据，再形成三维立体模型。

本发明的三维立体成像方法具有良好的可视性，可以将施工现场数百条杂乱无序的测线形成的多个二维平面的剖面数据集合成一个三维立体模型，能够非常直观的反映出探测区域异常体的真实情况，也便于后期的解释、分析与指导。

本发明的三维立体成像方法中，三维坐标转换方法，能将地表探测获取到的二维剖面数据转化为三维坐标系下的电阻率数据，方便后期进行数据分析以及三维建模，为后期建模分析解释工作的开展提供了极大的便利。

附图说明

构成本发明的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。

图1是本公开实施例的地表探测三维成像方法流程示意图；

图2(a)是本公开实施例的三维坐标转化方法中测线平面平行XOZ平面时的原理图；

图2(b)是本公开实施例的三维坐标转化方法中测线平面平行YOZ平面时的原理图；

图2(c)是本公开实施例的三维坐标转化方法中测线平面倾斜时的原理图。

具体实施方式

应该指出，以下详细说明都是例示性的，旨在对本发明提供进一步的说明。除非另有指明，本发明使用的所有技术和科学术语具有与本发明所属技术领域的普通技术人员通常理解的相同含义。

需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本发明的示例性实施方式。如在这里所使用的，除非本发明另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合；

为了方便叙述，本发明中如果出现“上”、“下”、“左”“右”字样，仅表示与附图本身的上、下、左、右方向一致，并不对结构起限定作用,仅仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的设备或元件必须具有特定的方位，以特定的方位构造和操作，因此不能理解为对本发明的限制。

术语解释部分:本发明中如出现术语“安装”、“相连”、“连接”、“固定”等，应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或为一体；可以是机械连接，也可以是电连接，可以是直接连接，也可以是通过中间媒介间接相连，可以是两个元件内部连接，或者两个元件的相互作用关系，对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明的具体含义。

正如背景技术所介绍的，现有技术中存在不足，为了解决如上的技术问题，本发明提出了一种地表综合物探的三维立体成像方法及系统，可以应用于岩溶、孤石、桩基、矿产等勘查工作中。该方法将现场多条测线形成的冗杂的二维剖面的电阻率数据通过三维坐标转换以及克里金插值法整合成像到了一个三维的立体成像系统中，使勘探数据具有良好的直观性与可视性。

本发明的一种典型的实施方式中，如图1所示，提供一种地表综合物探的三维立体成像方法，其包括以下步骤：

先通过跨孔法、井地法、高密度电法、瞬变电磁法或者地表综合物探方法，获取得到探测现场的多个二维剖面的探测数据，通过现场采集到的数据作为正演基础，通过反演对数据进行解释，得到多个二维剖面电阻率数据；

通过探测以及反演获得了多个二维剖面电阻率数据之后，再通过特定的算法，将二维剖面的数据点的平面直角坐标转换为特定三维直角坐标系下的三维直角坐标，获取三维空间直角坐标系下的数据点的电阻率数据；在获得三维坐标系下的电阻率数据之后通过克里金插值法，将获取的多个电阻率数据点转化为一个三维立体模型；

在获取三维立体电阻率模型之后，可以对三维立体电阻率模型的异常区进行解释，从而指导现场施工。

具体的，如图1所示，本实施例的地表探测三维成像方法流程，包括以下步骤：

1:根据地表探测现场具体的施工情况，选择探测方案以及测线布设方案。可以采用跨孔法、井地法、高密度电法、瞬变电磁法、地表综合物探方法中的一种进行地表探测，探测后可得到地表探测现场的多个二维剖面的电阻率探测数据。

如果采用跨孔法，则发射孔和接收孔按一定跨孔间距平行分布，且发射孔和接收孔的深度均大于探测目标深度，在发射孔中布置一定数量的发射电极，在接收孔中布置一定数量的接收电极，通过给发射电极通电，接收电极接收来获取探测数据；

如果采用井地法，则发射电极位于发射孔中，接收电极位于地表，通过给发射电极通电，接收电极接收来获取探测数据；

如果采用高密度电法，则需要将全部电极置于观测剖面的各测点上，然后利用程控电极转化装置和微机工程电测仪便可实现数据的快速采集和自动采集；

如果采用瞬变电磁法，则需要利用不接地回线或接地电极向地下发送脉冲式一次电磁场，用线圈或接地电极观测该脉冲电磁场感应的地下涡流产生的二次电磁场的空间和时间分布；

如果采用地表综合物探方法，则需要结合上述多种测线布置方案以及探测方式。

2:对采集到的二维剖面的探测数据进行地球物理反演。通过步骤1中现场采集到的数据作为正演基础，不断调整模型参数，使其模型响应向观测数据响应，即通过反演对数据进行解释，从而得到探测平面的二维剖面的电阻率数据。

3:将获取的二维剖面的电阻率数据点坐标通过特定的公式进行三维坐标转化，获得空间三维直角坐标系下每个数据点的三维直角坐标，即得到三维坐标系下的电阻率数据。

4:将获取的多个剖面下的三维直角坐标系下的电阻率数据点通过克里金插值法形成模型。该插值法考虑的是空间属性在空间位置上的变异分布，确定对一个待插点值有影响的距离范围，然后用此范围内的采样点来估计待插点的属性值。根据样品空间位置不同、样品间相关程度的不同，对每个样品品位赋予不同的权，进行滑动加权平均，以估计中心块段平均品位，最后将所有三维坐标系下的数据点集合成像成为一个三维立体模型。

5:通过对建立的三维模型进行分析，对异常区域进行解释，从而指导工程的进行。

对于上述步骤中的三维坐标转化的过程，如图2(a)-图2(c)所示，包括以下步骤：

本实施例构建了三种常规测线方式下的坐标转换方式。首先进行三维坐标轴的建立，根据施工现场具体情况选择原点后，X方向为水平距离方向，以水平向右为正方向；Y方向为探测深度方向，以平行地面朝探测方向为正方向；Z方向为垂直距离方向，以垂直地面往下为正方向。

测线是在探测现场人工按照特定方向在地表铺设的，将人工铺设测线的三种可能的铺设方向分为三种测线布置情况，测线平面在三维坐标中的位置与倾角就是地表测线的位置与倾角。

测线平面即指地表测线正下方的二维剖面，通过人工铺设的地表测线激发出的电信号可以获取位于测线正下方的二位剖面上的电阻率数据。

第一种测线情况，即测线情况(a)，如图2(a)所示，为测线平面平行于XOZ平面的情况，测线情况(a)又分为测线平面在坐标原点附近以及测线平面在原点远处两种类型；

第二种测线情况，即测线情况(b)，如图2(b)所示，为测线平面平行于YOZ平面的情况，测线情况(b)又分为测线平面在坐标原点附近以及测线平面在原点远处两种类型；

第三种测线情况，即测线情况(c)，如图2(c)所示，为测线平面倾斜于YOZ平面(垂直于XOY平面)，且与YOZ平面夹角为α，测线情况(c)又分为测线平面与Y轴正方向夹角为α以及测线平面与Y轴负方向夹角为α两种情况。

本实施例构建了三种常规测线方式下的坐标转换方式，在获得所布置的测线位置坐标以及数据点在测线平面上的位置之后，通过数学公式运算，即可获得所测区域内任意点的三维坐标，为后期三维成图奠定了方法基础，使得实际工程的探测成果三维立体化，解译更加便捷可行。

三维坐标转换通过以下公式进行转换：

其中X、Y、Z为以O为原点的最终三维坐标，X1为测线起点距离坐标原点的水平距离，Y1为测线起点距离坐标原点的纵向埋深，Z1为测线起点距离坐标原点的垂直高度，X＇为原始数据点的水平长度，Y＇为原始数据点的探测深度，Z＇初值为0，R1为位置矩阵，R2为数据点矩阵。

R1为位置矩阵，

R2为数据点矩阵，根据测线平面情况的不同，R2取值不同。

针对前述三种不同测线情况下的R1、R2按下表选取：

若测线平面平行于三维坐标系XOZ平面，则R2为

若测线平面平行于三维坐标系YOZ平面，则R2为

若测线平面垂直于XOY平面，测线平面向右倾斜，且与YOZ平面Y轴负方向夹角夹角为α，则R2为

若测线平面垂直于XOY平面，测线平面向左倾斜，且与YOZ平面Y轴正方向夹角夹角为α，则R2为

以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种地表综合物探的三维立体成像方法，其特征是，包括以下步骤：

采集地表探测现场的多个二维剖面的探测数据；

对探测数据通过地球物理反演形成二维剖面电阻率数据；

将二维剖面电阻率数据进行三维坐标转换，获取三维坐标系的电阻率数据；

利用克里金插值法将三维坐标系的电阻率数据转化成三维立体模型。

2.如权利要求1所述的地表综合物探的三维立体成像方法，其特征是，采集二维剖面的探测数据时，采用跨孔法、井地法、高密度电法、瞬变电磁法、地表综合物探方法中的一种进行地表探测。

3.如权利要求1所述的地表综合物探的三维立体成像方法，其特征是，三维坐标转换的过程为：

确定原点，建立三维坐标轴；

获取测线位置坐标及二维剖面电阻率数据点在测线平面的位置；

通过三维坐标转换公式获得二维剖面电阻率数据点的三维坐标。

4.如权利要求3所述的地表综合物探的三维立体成像方法，其特征是，所述三维坐标轴X方向为水平距离方向，以水平向右为正方向；Y方向为探测深度方向，以平行地面朝探测方向为正方向；Z方向为垂直距离方向，以垂直地面往下为正方向。

5.如权利要求3所述的地表综合物探的三维立体成像方法，其特征是，所述三维坐标转换公式为：

其中X、Y、Z为以O为原点的最终三维坐标，X1为测线起点距离坐标原点的水平距离，Y1为测线起点距离坐标原点的纵向埋深，Z1为测线起点距离坐标原点的垂直高度，X＇为原始数据点的水平长度，Y＇为原始数据点的探测深度，Z＇初值为0，R1为位置矩阵，R2为数据点矩阵。

6.如权利要求5所述的地表综合物探的三维立体成像方法，其特征是，所述R1为

7.如权利要求5所述的地表综合物探的三维立体成像方法，其特征是，若测线平面平行于三维坐标系XOZ平面，则R2为

8.如权利要求5所述的地表综合物探的三维立体成像方法，其特征是，若测线平面平行于三维坐标系YOZ平面，则R2为

9.如权利要求5所述的地表综合物探的三维立体成像方法，其特征是，若测线平面垂直于XOY平面，且与YOZ平面Y轴正方向夹角夹角为α，则R2为

若测线平面垂直于XOY平面，且与YOZ平面Y轴负方向夹角夹角为α，则R2为

10.一种地表综合物探的三维立体成像系统，其特征是，包括：

采集模块，用于采集地表探测现场的多个二维剖面的探测数据；

反演模块，用于对探测数据通过地球物理反演形成二维剖面电阻率数据；

坐标转换模块，用于将二维剖面电阻率数据进行三维坐标转换，获取三维坐标系的电阻率数据；

转化模块，用于利用克里金插值法将三维坐标系的电阻率数据转化成三维立体模型。

Images