CN106772625B - 一字形三维电法勘探方法 - Google Patents
一字形三维电法勘探方法 Download PDFInfo
- Publication number
- CN106772625B CN106772625B CN201611213018.5A CN201611213018A CN106772625B CN 106772625 B CN106772625 B CN 106772625B CN 201611213018 A CN201611213018 A CN 201611213018A CN 106772625 B CN106772625 B CN 106772625B
- Authority
- CN
- China
- Prior art keywords
- electrical prospecting
- exploration
- linear type
- dimensional electrical
- dimensional
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/02—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with propagation of electric current
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
一字形三维电法勘探方法,本发明基于电法勘探区的勘探目标体特征、地形地貌起伏特征,通过布置一字形的三维电法勘探供电点,利用电法勘探发射机通过供电点的电极向地下供电,采用电法勘探接收机同时接收所有测量点的电场信号,求取一字形的每个供电点供电时所有测量点的视电阻率和激电参数,从而实现一字形三维电法勘探。本发明可克服二维电法勘探效果较差、常规三维电法勘探效率低的缺点,又继承二维电法勘探效率高、常规三维电法勘探效果优的优点,在勘探效率和勘探效果两者之间找到一个更优的平衡点。
Description
技术领域
本发明涉及一种勘查地球物理领域的三维电法勘探方法。
背景技术
在勘查地球物理领域,电法勘探分为一维、二维和三维电法勘探。一维电法勘探由于勘探精度低,效果差,目前基本不使用;二维电法勘探精度和效果均优于一维,但勘探效率低于后者;常规三维电法勘探分为完整形、十字交叉和Γ形布置测量电极的三维电法勘探,其精度和效果均优于一维和二维,但由于供电点和测量点数量上远远多于一维和二维,从而导致常规三维电法勘探效率远远低于一维和二维,因此,该方法虽然有很好的勘探效果,但仅在部分非常关键勘探区且经费非常充足的情况下少量采用。
二维电法勘探需要在每条勘探线上布置电法勘探供电系统和测量系统,且后期数据处理过程中均把每条勘探线的数据作为二维数据进行单独处理,由于实际勘探区很少存在真实的二维地质异常体,更多的是三维地质异常体,故二维电法勘探成果均为近似成果,从而导致其勘探效果较差。但由于二维电法勘探在一条勘探线的供电点上供电,且只在同一条勘探线上开展电法数据采集,其勘探效率较高,勘探效果优于一维,故是目前常用的电法勘探方法。但由于该方法每条勘探线均需要布置供电系统,导致布置供电系统的工作量较大;每布置一条供电线仅测量一条测量线,若勘探区地形复杂、植被发育,则会导致难以甚至无法布置供电线,从而导致二维电法勘探无法开展的结果,限制了该方法在复杂地形的应用。
假设某个二维电法勘探区每条测线有m(m为大于1的正整数,以下相同)个供电点和m个测量点,有n(n为大于1的正整数,以下相同)条测线,采用单极-单极工作装置进行勘探,则二维电法勘探中需要布置m×n个供电点,以便采集n×m×(m-1)个二维电法勘探数据(供电点本身不予测量)。若基于电法勘探互换原理,可以对上述勘探工作进行简化,供电点数为n×(m-1)个,采集二维电法勘探数据为n×(m-1)×(m-1)/2个。若为获取更深部的电性特征,则供电点需要外延到测线端点以外,供电点数量会相应增加,进一步降低了该方法的勘探效率。
完整形三维电法勘探方法是对整个勘探区的供电点均分别供电一次,且每个供电点供电时,需要对整个勘探区所有测量点采集数据。该方法把实际勘探区内的地质异常体均作为三维地质异常体考虑,更符合实际地质情况。但由于完整形三维电法勘探的供电点需要覆盖整个三维电法勘探区内的所有测量点,且供电点的数量等于测量点数量,需要对整个勘探区内的所有供电点进行分别供电,每个供电点供电时需要对整个勘探区内的所有测量点采集数据,由于供电点数量过多,从而导致工作量成倍增加。
假设某个三维电法勘探矩形区网格数为m×n(其中m为X方向一条测量线的网格数,n为Y方向一条测量线的网格数),采用单极-单极工作装置进行勘探,则完整形三维电法勘探中需要布置m×n个供电点,采集m×n×(m×n-1)个测量数据(供电点本身不予测量)。基于电法勘探互换原理,可以对上述野外采集数据量减少一半,即为m×n×(m×n-1)/2个数据。
完整形三维电法勘探数据量和工作量远大于二维电法勘探,从而导致在实际勘探过程中需要投入大量的人力、物力和时间,效率低,成本高。另外,在野外实际勘探过程中,由于三维电法勘探工作量大,很难在一天之内完成整个勘探区的工作量,一旦需要跨天完成,则每天均需要布置和回收整个勘探区的供电系统和测量系统,存在很大的人力、物力和时间的浪费,也会导致其勘探效率更低,适用性不强。而且,由于完整形三维电法勘探对于野外勘探区的地形要求高,若地形复杂、植被发育,则给布置电法勘探供电系统和测量系统带来非常大的困难,有可能导致无法开展完整的三维电法勘探工作。因此,完整形三维电法勘探虽然在效果和精度上具有明显的优势,但由于其效率过低,成本过高,对地形要求过高,从而导致该方法在野外实际勘探工作中很少采用。
基于完整形三维电法勘探效率非常低、适用性不强的缺点,目前有人研究了十字交叉和Γ形布置测量电极的三维电法勘探方法。其着眼点主要是通过减少某个供电点供电时所测量的测量点数量,从而实现减少野外工作量,达到提高勘探效率的目的。但十字交叉形测量电极和Γ形测量电极的三维电法勘探方法没有减少供电点,其三维电法勘探的供电点仍需要覆盖整个三维电法勘探区内的所有测量点,供电点的物理点数量等于测量点的物理点数量,由于供电点数量过多,且目前供电技术采用的是分时分点供电,每次仅能对一个供电点供电,从而导致供电系统的工作量无法减少。该方法仅采用某个供电点供电时选择部分测量点采集电法数据,但为完成所有供电点供电时的数据采集工作,整个勘探区每个测量点的物理点均需要布置,只是某些供电点供电时,有部分测量点未采集电法数据,并没有减少测量点的物理点布置。故上述十字交叉形测量电极和Γ形测量电极三维电法勘探方法相对于完整形三维电法勘探方法的勘探效率有所提高,但由于是通过减少某个供电点供电时所测量的测量点数量,因此仍需要对整个勘探区的所有供电点进行分时供电,且需要布置整个勘探区测量点的物理点,故其勘探效率提高不明显。
假设某个三维电法勘探矩形区网格数为m×n(其中m为X方向一条测量线的网格数,n为Y方向一条测量线的网格数),采用单极-单极工作装置进行电法勘探,则十字交叉形测量电极三维电法勘探中需要布置(m×n-1)个供电点,采集个测量数据;Γ形测量电极三维电法勘探中同样需要布置(m×n-1)个供电点,采集个电法数据。以上两种方法的数据采集量依然很大,效率提高不明显。
发明内容
本发明的目的是基于电法勘探区的勘探目标体特征、地形地貌起伏特征,通过布置一字形的三维电法勘探供电点,减少三维电法勘探供电点数量,实现一字形三维电法勘探方法,提高三维电法勘探效率和效果。
由于电法勘探供电系统每次仅能对一个供电点供电,目前无法通过增加多套供电系统同时供电,实现减少布置电法勘探供电系统的时间和成本。对于常规三维电法勘探而言,布置供电系统的劳动强度大,效率低,故若能实现减少供电点的数量,则能明显提高三维电法勘探效率。另由于供电系统为满足电法勘探采集信号的质量要求,其供电线和供电极均比电法勘探测量系统复杂和笨重,故减少供电点的数量,相对于减少测量点的数据采集量而言,其勘探效率改善更加明显。
由于电法勘探测量系统可以实现同时采集,通过多通道电法勘探测量系统或多个电法勘探测量系统,可以实现减少整个三维电法勘探数据采集时间和成本,从而提高野外勘探效率,降低勘探成本;而且由于不减少三维电法勘探所有测量点的数据采集量,又能保证三维电法勘探效果。
本发明利用二维电法勘探效率高,常规三维电法勘探精度和效果优的优点,基于电法勘探区的勘探目标体特征、地形地貌起伏特征,通过一字形三维电法勘探供电点布置规则,减少三维电法勘探供电点数量,实施对勘探区内所有测量点的三维电法勘探,实现提高三维电法勘探效率和改善勘探效果的目的。本发明可克服二维电法勘探效果较差、常规三维电法勘探效率低的缺点,又继承二维电法勘探效率高、常规三维电法勘探效果优的优点,在勘探效率和勘探效果两者之间找到一个更优的平衡点。
本发明提出的一字形三维电法勘探方法具体步骤为:
(1)先行分析勘探区的前期地质成果,推断区内主要地质异常体的走向或主轴方向,近似垂直于走向或主轴方向设定一字形三维电法勘探的X方向,近似平行于走向或主轴方向设定一字形三维电法勘探的Y方向,其中X方向与Y方向相互垂直,且均为水平方向;
(2)根据勘探要求,确定勘探区的X方向与Y方向的勘探长度,并确定勘探区内X方向测量点间距c和Y方向测量点间距d的数值,从而确定勘探区网格参数a和b的数值;
(3)基于勘探区的勘探目标体特征、地形地貌起伏特征,采用一字形的电法勘探供电点网格坐标(x,y)和实际相对坐标(XC,YC)计算规则,确定一字形三维电法勘探所有供电点网格坐标(x,y)和实际相对坐标(XC,YC),并布置供电点;
(4)基于勘探区的勘探目标体特征、地形地貌起伏特征,整个勘探区内按照勘探要求及X方向测量点间距c、Y方向测量点间距d、网格参数a和b的数值,布置完整的三维电法勘探测量点;
(5)在上述选定的一字形三维电法勘探供电线上的每个供电点分别供电,采用多通道电法勘探测量系统或多个电法勘探测量系统对整个勘探区内的一字形三维电法勘探测量点进行电法数据采集,获取到每个供电点供电时勘探区所有测量点的一字形三维电法勘探数据;
(6)把所有电法勘探测量数据整合为一字形三维电法勘探数据,进行一字形三维电法数据处理,获取整个勘探区的一字形三维电法勘探成果,从而实现整个勘探区的一字形三维电法勘探。
其中步骤(1)中的近似垂直走向或主轴方向设定一字形三维电法勘探的X方向,近似平行走向或主轴方向设定一字形三维电法勘探的Y方向的具体含义为X方向与主要地质异常体的走向或主轴方向夹角在60~120度之间,并尽量接近或等于90度;Y方向与主要地质异常体的走向或主轴方向夹角在-30~30度之间,并尽量接近或等于0度;其中X方向与Y方向相互垂直,且均为水平方向。
其中步骤(2)中勘探区网格参数a和b的数值的确定规则如下:设一字形三维电法勘探矩形区X方向长L米,Y方向长M米,勘探区内X方向测量点间距为c米,Y方向测量点间距为d米,以整个勘探区的中心点为原点,则整个勘探区X方向和Y方向的网格参数a和b计算过程如下:设
其中步骤(3)中一字形三维电法勘探供电点布置规则分为两种情况:第一种情况为,三维电法勘探矩形区X方向测量线长度达到勘探目标体预测宽度的2倍及以上;第二种情况为,三维电法勘探矩形区X方向测量线长度为勘探目标体预测宽度的2倍以下。
1、第一种情况,当一字形三维电法勘探矩形区X方向测量线长度达到勘探目标体预测宽度的2倍及以上时,本发明的一字形三维电法勘探供电点布置规则如下:
设该勘探矩形区的网格左下角坐标为(-a,-b),右上角坐标为(a,b),则一字形三维电法勘探供电点位置的网格坐标(x,y)公式设置为:
其中公式(2)中的x、y、a、i均为整数。
根据以上规则,得出当一字形三维电法勘探矩形区X方向测量线长度为勘探目标体的预测宽度2倍及以上时,以整个勘探区的中心点为原点,一字形三维电法勘探供电点位置的实际相对坐标(XC,YC)公式如下:
(XC,YC)=(x×c,y×d)=(x×c,0) (3)
其中公式(3)中的x、y为公式(2)的计算结果。
2、第二种情况,当三维电法勘探矩形区X方向测量线长度为勘探目标体的预测宽度2倍以下时,本发明的一字形三维电法勘探供电点布置规则如下:
若勘探区地形限制较多,一字形三维电法勘探矩形区X方向测量线长度为勘探目标体的预测宽度2倍以下,则为进一步减少勘探区边界测量点的边界效应,提高整个勘探区的勘探精度,所有供电点采取与勘探区内测量点相同的间隔距离布置。
设该勘探区的网格左下角坐标为(-a,-b),右上角坐标为(a,b),则一字形三维电法勘探供电点位置的网格坐标(x,y)公式设置为:
(x,y)=(x,0);|x|≤2a (4)
其中公式(4)中的x、a均为整数。
根据以上规则,得出当X方向测量线长度为勘探目标体预测宽度的2倍以下时,以整个勘探区的中心点为原点,一字形三维电法勘探供电点位置的实际相对坐标(XC,YC)公式如下:
(XC,YC)=(x×c,y×d)=(x×c,0) (5)
其中公式(5)中的x、y为公式(4)的计算结果。
本发明的一字形三维电法勘探供电装置采用单极装置,即采用电法勘探意义上的一个无穷远供电电极作为三维电法勘探的一个固定供电电极,另一个供电电极则为按照本发明提出的一字形三维电法勘探供电点布置规则布置的移动供电电极。
本发明的一字形三维电法勘探测量装置可以采用单极和偶极测量装置两种。单极测量装置采用电法勘探意义上的一个无穷远测量电极作为一个固定的三维电法勘探测量电极,另一个测量电极为勘探区内移动的三维电法勘探测量电极。偶极测量装置则是两个电法勘探测量电极均为勘探区中移动的三维电法勘探测量电极。
本发明提出的一字形三维电法勘探方法为减少勘探区边界测量点数据量少所导致的边界效应问题,提出将一字形三维电法勘探的供电点布置到测量线以外,以便增加勘探区边界测量点的数据量,从而减弱边界效应引起的勘探区边界测量点的勘探精度降低的问题。
本发明提出的一字形三维电法勘探方法由于把所有供电点位置布置在直线上,有利于在实际勘探工作中进行供电点位置的改变,也能节省供电导线,提高改变供电点位置的工作效率和节省勘探成本。
本发明提出的一字形三维电法勘探方法,通过选取一字形三维电法勘探供电点,达到减少供电点的数量,规避电法勘探供电系统无法同时供电且笨重的劣势;提出不减少三维电法勘探测量点数量的方式,充分挖掘电法勘探测量系统能同时采集的优势,并保留三维电法勘探效果优的优点,从而实现一字形三维电法勘探的高效率和更优勘探效果。该方法同时具有二维的高勘探效率和三维的高精度勘探效果的优点。
假设某个三维电法勘探矩形区网格数为m×n(其中m为X方向一条测量线的网格数,n为Y方向一条测量线的网格数),采用单极-单极工作装置进行勘探,当三维电法勘探矩形区X方向测量线长度达到勘探目标体的预测宽度2倍及以上时,则一字形三维电法勘探只需布置个供电点,共采集个三维电法勘探数据。
本发明提出的一字形三维电法勘探方法由于大量减少了供电点数,其勘探效率远远高于常规的三维电法勘探,但又采用三维电法勘探模式,其勘探效果优于二维电法勘探。
若勘探区过大,则对勘探区进行分区处理,每个小分区按照上述原则进行勘探区网格设计,并开展每个小分区的一字形三维电法勘探,获取所有小分区的一字形三维电法勘探数据,最后进行整个勘探区的一字形三维电法勘探数据处理,得到整个勘探区的一字形三维电法勘探成果。
若勘探区已知地质情况少,难以判断主要地质异常体的走向或主轴方向,则先对勘探区进行初步分析,适当加大X方向测量点间距c和Y方向测量点间距d,选择合适位置开展一字形三维电法勘探试探性工作,获取初步成果后,再对主要异常区进行减小X方向测量点间距c和Y方向测量点间距d的加密型一字形三维电法勘探,以便实现勘探效率和勘探效果的最优化。
附图说明
图1为本发明的一字形三维电法勘探野外布置示意图;
图2为1号地电模型单极-单极数值模拟结果的二维电法勘探电阻率反演成果;
图3为1号地电模型单极-偶极数值模拟结果的二维电法勘探电阻率反演成果;
图4为1号地电模型单极-单极数值模拟结果的一字形三维电法勘探电阻率反演成果;
图5为1号地电模型单极-偶极数值模拟结果的一字形三维电法勘探电阻率反演成果;
图6为1号地电模型单极-单极数值模拟结果的完整形三维电法勘探电阻率反演成果;
图7为1号地电模型单极-偶极数值模拟结果的完整形三维电法勘探电阻率反演成果;
图8为1号勘探线单极-偶极二维电法勘探实验的电阻率反演成果;
图9为1号勘探线单极-偶极一字形三维电法勘探实验的电阻率反演成果。
具体实施方式
以下结合图1和实施例对本发明做进一步说明,以下说明采用一字形三维电法勘探单极供电-单极测量装置。
本方法的供电系统和测量系统能完全兼容目前已有的电法勘探设备,且测量系统能采用多通道电法勘探测量系统或多套电法勘探测量系统。
实施例1
图1中的1为无穷远供电电极,2为无穷远供电线,3为电法勘探发射机,4为连接勘探区供电点的供电线,5所指的十字形符号为电法勘探移动供电极的供电点,6所指的圆形符号为电法勘探测量电极的测量点,部分供电点与测量点位置重叠。
假设三维电法勘探矩形区X方向长L=400米,Y方向长M=400米,勘探区内X方向测量点间距c=40米,Y方向测量点间距d=40米,勘探区内勘探目标体的预测宽度为100米,由于三维电法勘探矩形区X方向测量线长度达到勘探目标体的预测宽度4倍,则根据公式(1)
求得a=5,b=5,则该勘探区的网格左下角坐标为(-5,-5),右上角坐标为(5,5)。
根据公式(2)
代入a、b数值,求得一字形三维电法勘探供电点位置的网格坐标(x,y)为
其中公式(6)中的x、i均为整数。
根据以上所得,求出以整个勘探区的中心点为原点的一字形三维电法勘探供电点位置的实际相对坐标(XC,YC)为:
(XC,YC)=(x×40,y×40)=(x×40,0) (7)
其中公式(7)中的x为公式(6)的计算结果。
根据以上的计算结果,一字形三维电法勘探单极供电-单极测量方式的具体操作步骤如下:
(1)分析勘探区的前期地质成果,推断区内主要地质异常体的走向或主轴方向,垂直于走向或主轴方向设定为一字形三维电法勘探的X方向,平行于走向或主轴方向设定为一字形三维电法勘探的Y方向,其中X方向与Y方向相互垂直,且均为水平方向;
(2)根据以上求取的一字形三维电法勘探供电点网格坐标(x,y)和实际相对坐标(XC,YC),并结合勘探区实际情况,在勘探区及周边合理布置供电线上的移动供电电极的供电点5;
(3)在勘探区外按照电法勘探原理和要求选择无穷远供电点,并布设无穷远供电电极1和供电线2;
(4)选择某一个一字形三维电法勘探移动供电电极的供电点5(如图1所示的A-13)布设供电导线4和移动供电电极;
(5)在勘探区外按照电法勘探原理和要求选择电法勘探无穷远测量点,并布设无穷远测量电极和导线;
(6)在勘探区内按网格方式布置所有一字形三维电法勘探测量电极的测量点6,并布置测量电极和导线;
(7)通过电法勘探发射机3向无穷远供电电极1和移动供电电极的供电点5(如图1所示的A-13)上的电极供电,记录供电电流值,采用多通道或多套电法勘探测量系统记录下该供电点供电时勘探区内每一个测量电极的测量点6和无穷远测量电极之间的电法数据,待采集完毕,则通知供电系统改变供电电极的供电点5位置(如图1所示的A-9、A-7、A-5、A-4、……A13),再重复上述电法数据采集过程,直至采集完每个供电点供电时所有测量点的电法数据;
(8)获取到上述一字形三维电法勘探的所有点的电法勘探数据后,求取一字形三维电法勘探的视电阻率或视激电等电法参数,通过后期数据处理,获取到勘探区的一字形三维电法勘探成果。
本发明经勘探效率对比实验、数值模拟勘探效果对比实验、勘探效率和勘探效果野外对比实验,可以发现该发明提出的方法在勘探效率和勘探效果上都有明显的提高。
勘探效率对比实验
假设对11×11的网格划分勘探区开展二维、完整形三维、十字交叉形测量电极三维、Γ形测量电极三维和一字形三维电法勘探,则每种方法所开展的工作量对比如表1所示。
表1多种电法勘探工作量对比表
从上表可以明显看出,一字形三维电法勘探的电法数据采集量虽然接近二维电法勘探的两倍,但由于该方法仅需要在勘探区布置一条供电线,就可实现勘探区所有测量点的数据采集,且仅布设17个供电点,而二维电法勘探需要在11条测线上布置11次供电系统,共需布设121个供电点,故一字形三维电法勘探勘探效率反而比二维电法勘探效率要高。更重要的是,一字形三维电法勘探采用三维电法勘探模式,其勘探效果比二维勘探更真实可靠。
数值模拟勘探效果对比实验
为对比一字形三维电法勘探与二维电法勘探,以及一字形三维电法勘探与完整形三维电法勘探的勘探效果,特设定1号地电模型:300×300米的勘探区,X方向测量点间距c和Y方向测量点间距d均为20米,采用单极供电-单极测量的电法勘探装置进行数值模拟和反演。围岩电阻率为1000欧姆·米,低阻异常体电阻率为10欧姆·米。低阻异常体为一边长为100米的立方体,低阻异常体水平放置,其平面坐标范围为(-50,-50)~(50,50),顶部埋深20米,厚度100米。
本次仅对X方向进行室内数值模拟和反演对比。
对1号地电模型通过数值模拟获取单极供电-单极测量和单极供电-偶极测量电法勘探数据,再对获取的数据整理为相应的二维电法勘探数据、一字形三维电法勘探数据以及完整形三维电法勘探数据,并分别采用二维电法勘探反演程序和三维电法勘探反演程序进行反演,获取到三种方法两种装置的反演成果。
从三种方法两种装置的反演成果图(图2、3、4、5、6、7,图中的等值线为电阻率等值线)可以得出,二维电法勘探反演成果对于低阻异常体形态反应基本准确,但对于围岩的电阻率的反应出现明显的畸变现象。而一字形三维电法勘探反演成果对于低阻异常体形态的反应优于二维电法勘探成果,与完整形三维电法勘探成果类似,获得的低阻异常形态基本上为一矩形低阻异常体,与1号地电模型实际情况类似;对于围岩电阻率的反应比二维电法勘探成果更真实,未出现明显畸变现象,与1号地电模型吻合度明显高于二维电法勘探成果。
通过本次对1号地电模型的二维电法勘探、一字形三维电法勘探和完整形三维电法勘探进行的数值模拟及反演,可以发现,一字形三维电法勘探成果与1号地电模型实际情况的吻合度比二维电法勘探成果要高,说明前者的勘探效果更优。一字形三维电法勘探成果与完整形三维电法勘探成果类似,与1号地电模型吻合度高,但前者勘探效率高于后者。
勘探效率和勘探效果野外对比实验
为验证一字形三维电法勘探的效果,在某地质情况比较明了的实验区开展了单极-偶极二维电法勘探和单极-偶极一字形三维电法勘探对比实验。本次实验共布设了1条供电线,3条勘探线,测量点距40米,勘探线距100米,每条测量线长度320米。由于该实验区地形复杂、植被发育,仅有实验区中间1号实验线适宜布置供电线,2号和3号实验线的大部分测量点位于植被发育、陡峭的山头,根据前期地质成果,分析该区勘探目标体宽度约为200米,由于勘探区测量点较少,勘探线布置测量点的长度为勘探目标体的预测宽度2倍以下,为减少勘探区边界测量点的边界效应影响,故本次供电点全部采用等间距布置。本次二维电法勘探实验共布置27个供电点,1号勘探线布置7个测量点。本次一字形三维电法勘探共布置27个供电点,每条勘探线7个测量点,3条勘探线共21个测量点,其中1号勘探线测量点和供电点位置与二维电法勘探完全重合。由于1号勘探线开展了二维电法勘探和一字形三维电法勘探,故选择1号勘探线为对比实验线。
由于本次实验中一字形三维电法勘探数据完全包括二维电法勘探数据,故本次实验仅在后期针对不同需要对实验数据整理出单极-偶极二维电法勘探数据和单极-偶极一字形三维电法勘探数据。一字形三维电法勘探增加了2名电法勘探接收机技术员和2台电法勘探接收机作为电法勘探测量系统。二维电法勘探获取了1号勘探线的成果,一字形三维电法勘探则同时获取了1号、2号、3号共3条勘探线的成果。
根据对1号勘探线上两种方法的电法勘探成果图(图8、9,图中的等值线为电阻率等值线)对比可以得出,两种方法对于该区的电阻率变化情况基本类似,比较符合该区存在一个往大号测量点方向倾斜的断层的地质认识。但一字形三维电法勘探成果对于断层的反应和倾向更准确,通过对1号勘探线一字形三维电法勘探成果的小号和大号测量点存在的高阻体位置和埋深,可以初步判断该断层为往大号测量点方向倾斜的正断层,与已知地质情况相符。而二维电法勘探成果对于断层的倾向和断层性质的判断则相对模糊,且在500~700号测量点间出现电阻率超过10000欧姆·米的高阻异常,而一字形三维电法勘探成果中该位置电阻率仅为1000欧姆·米左右,结合实验区的实际地质情况和前述二维电法勘探数值模拟结果中围岩电阻率出现明显畸变的现象,推断二维电法勘探所发现的高阻异常存在明显假象。通过本次野外实验对比,可以发现一字形三维电法勘探在勘探效率与二维电法勘探类似的情况下,能获取到精度更高、更符合地质成果的勘探成果。
通过以上对多种方法的勘探效率对比,数值模拟勘探效果对比,以及勘探效率和勘探效果野外实验对比,结果都证明了本发明提出的一字形三维电法勘探方法在勘探效率和勘探效果上都有明显改进。
Claims (4)
1.一字形三维电法勘探方法,其特征在于:基于电法勘探区的勘探目标体特征、地形地貌起伏特征,布置一字形的电法勘探供电点,利用电法勘探发射机通过供电点的电极向地下供电,采用电法勘探接收机同时接收多条勘探线上每个测量点的电场信号,求取每个供电点供电时所有测量点的视电阻率和激电参数,实现一字形三维电法勘探,具体步骤为:
(1)先行分析勘探区的前期地质成果,推断区内主要地质异常体的走向或主轴方向,近似垂直于走向或主轴方向设定一字形三维电法勘探的X方向,近似平行于走向或主轴方向设定一字形三维电法勘探的Y方向,其中X方向与Y方向相互垂直,且均为水平方向;
(2)根据勘探要求,确定勘探区的X方向与Y方向的勘探长度,并确定勘探区内X方向测量点间距c和Y方向测量点间距d的数值,从而确定勘探区网格参数a和b的数值;
(3)基于勘探区的勘探目标体特征、地形地貌起伏特征,采用一字形的电法勘探供电点网格坐标(x,y)和实际相对坐标(XC,YC)计算规则,确定一字形三维电法勘探所有供电点网格坐标(x,y)和实际相对坐标(XC,YC),并布置供电点;
(4)基于勘探区的勘探目标体特征、地形地貌起伏特征,整个勘探区内按照勘探要求及X方向测量点间距c、Y方向测量点间距d、网格参数a和b的数值,布置完整的三维电法勘探测量点;
(5)给选定的一字形三维电法勘探供电线上的每个供电点分别供电,采用多通道电法勘探测量系统或多个电法勘探测量系统对整个勘探区内的一字形三维电法勘探测量点进行电法数据采集,获取到每个供电点供电时勘探区所有测量点的一字形三维电法勘探数据;
(6)把所有电法勘探测量数据整合为一字形三维电法勘探数据,进行一字形三维电法数据处理,获取整个勘探区的一字形三维电法勘探成果,从而实现整个勘探区的一字形三维电法勘探。
2.如权利要求1所述的一字形三维电法勘探方法,其特征在于:步骤(1)中的近似垂直于走向或主轴方向设定一字形三维电法勘探的X方向,近似平行于走向或主轴方向设定一字形三维电法勘探的Y方向的具体含义为X方向与主要地质异常体的走向或主轴方向夹角在60~120度之间,并尽量接近或等于90度;Y方向与主要地质异常体的走向或主轴方向夹角在-30~30度之间,并尽量接近或等于0度;其中X方向与Y方向相互垂直,且均为水平方向。
3.如权利要求1或2所述的一字形三维电法勘探方法,其特征在于:步骤(2)中勘探区网格参数a、b的数值的确定规则如下:
设三维电法勘探矩形区X方向长L米,Y方向长M米,勘探区内X方向测量点间距为c米,Y方向测量点间距为d米,以整个勘探区的中心点为原点,则整个勘探区X方向和Y方向的网格参数a和b计算过程如下:
4.如权利要求3所述的一字形三维电法勘探方法,其特征在于:步骤(3)中的一字形三维电法勘探供电点网格坐标(x,y)和实际相对坐标(XC,YC)计算规则如下:
A)若电法勘探区的地形地貌起伏特征不影响勘探线的布置,一字形三维电法勘探矩形区X方向测量线的长度达到勘探目标体预测宽度的2倍及以上,则一字形的三维电法勘探供电点网格坐标(x,y)和实际相对坐标(XC,YC)计算规则如下:设该勘探矩形区的网格左下角坐标为(-a,-b),右上角坐标为(a,b),则一字形三维电法勘探供电点位置的网格坐标(x,y)设置为:
其中公式(2)中的x、y、a、i均为整数;
根据以上规则,得出当一字形三维电法勘探矩形区X方向测量线长度为勘探目标体预测宽度的2倍及以上时,以整个勘探区的中心点为原点,一字形三维电法勘探供电点位置的实际相对坐标(XC,YC)公式如下:
(XC,YC)=(x×c,y×d)=(x×c,0) (3)
其中公式(3)中的x、y为公式(2)的计算结果;
B)若电法勘探区的地形地貌起伏特征影响勘探线的布置,导致一字形三维电法勘探矩形区X方向测量线长度小于勘探目标体预测宽度的2倍,则一字形三维电法勘探供电点网格坐标(x,y)和实际相对坐标(XC,YC)计算规则如下:设该勘探区的网格左下角坐标为(-a,-b),右上角坐标为(a,b),则一字形三维电法勘探供电点位置的网格坐标(x,y)设置为:
(x,y)=(x,0);|x|≤2a (4)
其中公式(4)中的x、a均为整数;
根据以上规则,得出当一字形三维电法勘探矩形区X方向测量线长度小于勘探目标体预测宽度的2倍时,以整个勘探区的中心点为原点,一字形三维电法勘探供电点位置的实际相对坐标(XC,YC)公式如下:
(XC,YC)=(x×c,y×d)=(x×c,0) (5)
其中公式(5)式中的x、y为公式(4)的计算结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611213018.5A CN106772625B (zh) | 2016-12-25 | 2016-12-25 | 一字形三维电法勘探方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611213018.5A CN106772625B (zh) | 2016-12-25 | 2016-12-25 | 一字形三维电法勘探方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106772625A CN106772625A (zh) | 2017-05-31 |
CN106772625B true CN106772625B (zh) | 2018-06-22 |
Family
ID=58920861
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611213018.5A Active CN106772625B (zh) | 2016-12-25 | 2016-12-25 | 一字形三维电法勘探方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106772625B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110850484B (zh) * | 2019-08-22 | 2021-09-07 | 广东电网有限责任公司 | 一种线圈海缆探测装置及其探测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102010050143A1 (de) * | 2010-11-03 | 2012-05-03 | Qumon Gmbh | Elektrisches Verfahren zur zerstörungsfreien Erkundung und Überwachung unterirdischen Raums, Vorrichtung und Verwendung dazu |
CN105116452A (zh) * | 2015-08-24 | 2015-12-02 | 中国石油天然气集团公司 | 一种确定地质异常体电阻率和极化率的方法和装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101495836B1 (ko) * | 2014-08-23 | 2015-02-25 | 주식회사 혜안지반 | 3차원 전기비저항 탐사시스템 및 이를 이용한 터널의 붕락구간 탐사방법 |
-
2016
- 2016-12-25 CN CN201611213018.5A patent/CN106772625B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102010050143A1 (de) * | 2010-11-03 | 2012-05-03 | Qumon Gmbh | Elektrisches Verfahren zur zerstörungsfreien Erkundung und Überwachung unterirdischen Raums, Vorrichtung und Verwendung dazu |
CN105116452A (zh) * | 2015-08-24 | 2015-12-02 | 中国石油天然气集团公司 | 一种确定地质异常体电阻率和极化率的方法和装置 |
Non-Patent Citations (2)
Title |
---|
三维全空间坑道直流聚焦超前探测数值模拟;张力等;《物探与化探》;20110630;第35卷(第3期);第419-422页 * |
三维电阻率测量方法试验研究;汤洪志等;《地质与勘探》;20031031;第39卷;第126-129页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106772625A (zh) | 2017-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103454686B (zh) | 基于地层切片的小尺度沉积相进行储层预测的方法及系统 | |
CN106772620B (zh) | 十字形三维电法勘探方法 | |
CN105911603B (zh) | 基于天然电场的四维物探方法 | |
CN104597519B (zh) | 一种深部金矿阶梯式找矿方法 | |
CN104636980B (zh) | 针对河道砂油藏类型油气汇集条件的地球物理表征方法 | |
CN104330832A (zh) | 一种瞬变电磁快速三维人机交互反演方法 | |
CN206960673U (zh) | 一种探查河道溶洞的多维探测装置 | |
Calixto et al. | 3-D soil stratification methodology for geoelectrical prospection | |
CN106802432B (zh) | 土层结构的探查方法和装置 | |
CN206378448U (zh) | 基于综合物探方法的污染土检测测线布置结构 | |
CN106646624B (zh) | 非对称电法测深方法 | |
CN110333543A (zh) | 基于反射系数分析的低阻体解释及成像方法与系统 | |
CN107191181A (zh) | 一种基于电磁散射的井周界面探测方法 | |
CN107065019A (zh) | 应用于道路灾害和塌陷检测的三维电磁成像装置及使用方法 | |
CN107037492A (zh) | 一种地质数据分析建模方法 | |
CN106019386A (zh) | 一种回线源内比值视电阻率测定方法 | |
CN112415602A (zh) | 基于深度分辨率的隧道电阻率超前探测优化方法及系统 | |
CN203849425U (zh) | 地裂缝三维电法探测装置 | |
CN106772625B (zh) | 一字形三维电法勘探方法 | |
CN108873073B (zh) | 一种基于网络并行电法的三维跨孔电阻率层析成像方法 | |
CN106842332B (zh) | 一种综合电法勘探方法 | |
CN106707347B (zh) | 一种综合电法勘探方法 | |
CN106646622B (zh) | 井字形三维电法勘探方法 | |
CN106597552B (zh) | 米字形三维电法勘探方法 | |
CN113341467A (zh) | 基于多插值方法的矿井瞬变电磁三维显示方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |