一种瞬变电磁探测方法、装置和系统
技术领域
本发明涉及地质与勘探地球物理领域,尤其涉及一种瞬变电磁探测方法、装置和系统。
背景技术
瞬变电磁法(TEM,Transient Electromagnetic Method)是一种常见的基于电磁感应的地球物理方法,能够提供地下地电信息。大回线中心方式瞬变电磁法装置广泛应用于矿产、工程和环境调查,以及其它地质研究,在传统瞬变电磁研究中,为了使信噪比达到最大多使用大发射回线。一般情况下,根据地质目标的深度,这种装置的发送回线边长为100米到300米,发送电流10A到20A。在我国的山区工作时,很难布置足够大的发射回路进行传统的测量。为解决上述问题,研发出了小尺度、大电流发射回路系统中的设备,并得到了一些应用。比如:利用小回线装置(回线边长3米,电流10-20A)进行地下隧道掌子面探测;利用5米线圈进行煤田水文地质探测。申请号为03262603.7的中国专利将小回线的发射电流增大到1000-2000A,进一步提升了小回线大电流的探测深度.但由于这种装置的大电流问题的影响,在进行数据处理时,目前还没有成熟的处理方法。目前,采用较多的处理方式是利用已知地质条件分段进行响应模板校正的方法,视电阻率计算多采用重叠回线晚期视电阻率定义方式,深度计算采用理论扩散深度进行分析。对于利用地质条件进行分段校正的方式,由于对测区的地质已知信息的要求较高,并不适用于所有的探测区域,尤其是山区复杂地形地区,而且,分段模拟校正的速度较慢,不利于瞬变电磁资料的快速高效的解释。采用重叠回线晚期视电阻率晚期公式对早期数据的处理会引起较大的计算误差,影响瞬变电磁早期探测精度。在不考虑实际地电结构的情况下,探测深度对于无干扰的均匀半空间介质是正确的,对于野外实测资料的处理采用理论探测深度会引起深度认知的误区。
发明内容
本发明要解决的技术问题是克服目前小回线大电流的探测存在视电阻率计算误差大的问题,本发明提出一种瞬变电磁探测方法、装置和系统。
为了解决上述技术问题,本发明提供了一种瞬变电磁探测方法,包括:采集待评估区域第一回线的瞬变电磁数据;所述第一回线的边长小于长度阈值;
将所述的第一回线的瞬变电磁数据等效转换成第二回线中心方式瞬变电磁数据;所述第二回线的边长大于或者等于长度阈值;
按照所述第二回线中心方式瞬变电磁数据的解释方法,确定待评估区域的视电阻率和探测深度。
进一步地,将所述的第一回线的瞬变电磁数据等效转换成第二回线中心方式瞬变电磁数据包括:
将所述的第一回线的瞬变电磁数据等效转换成相同磁矩大小的第二回线中心方式瞬变电磁数据。
进一步地,确定待评估区域的视电阻率包括:
根据中心回线的视电阻率计算公式确定所述待评估区域的视电阻率:
所述中心回线的视电阻率计算公式为
其中,M是源发射磁矩,大小为发射电流与发射回线面积的乘积;q代表观测位置探头的有效接收面积,V(t)是测量的二次感应电压值,t是观测时间,μ0为自由空间的磁导率;
或者,采用二分法确定全时域视电阻率。
进一步地,所述采用二分法确定全时域视电阻率包括:
将实际观测的二次感应电压转换成的垂直磁场,预设单调函数Hz=Y(ρ),并确定视电阻率ρ的预设值域(ρ1,ρ2),确定对应观测值Hz0值的反函数ρ0,包括:设定函数η(ρ)=Hz0-Y(ρ),选择预设值域(ρ1,ρ2)中点ρ3,确定η(ρ3),利用函数的单调性,当η(ρ3)*η(ρ1)<0时,则判断ρ0在(ρ1,ρ3)内,当η(ρ3)*η(ρ2)<0时,则判断ρ0在(ρ3,ρ2)内,如此循环,直到η(ρ3)小于设定的精度ε时,则此时的ρ3即为对应观测值Hz0值所对应的视电阻率值ρ0。
进一步地,确定待评估区域的探测深度包括:
利用公式或者确定所述探测深度。
为了解决上述技术问题,本发明还提供了一种瞬变电磁探测装置,包括:
采集模块,用于采集待评估区域第一回线的瞬变电磁数据;所述第一回线的边长小于长度阈值;
等效模块,用于将所述的第一回线的瞬变电磁数据等效转换成第二回线中心方式瞬变电磁数据;所述第二回线的边长大于或者等于长度阈值;
计算模块,按照所述第二回线中心方式瞬变电磁数据的解释方法,确定待评估区域的视电阻率和探测深度。
优选地,所述等效模块用于:
将所述的第一回线的瞬变电磁数据等效转换成相同磁矩大小的第二回线中心方式瞬变电磁数据。
优选地,所述计算模块包括:
视电阻率计算单元,用于根据中心回线的视电阻率计算公式确定所述待评估区域的视电阻率:
所述中心回线的视电阻率计算公式为
其中,M是源发射磁矩,大小为发射电流与发射回线面积的乘积;q代表观测位置探头的有效接收面积,V(t)是测量的二次感应电压值,t是观测时间,μ0为自由空间的磁导率;
或者,采用二分法确定全时域视电阻率。
优选地,视电阻率计算单元采用二分法确定全时域视电阻率包括:
将实际观测的二次感应电压转换成的垂直磁场,预设单调函数Hz=Y(ρ),并确定视电阻率ρ的预设值域(ρ1,ρ2),确定对应观测值Hz0值的反函数ρ0,包括:设定函数η(ρ)=Hz0-Y(ρ),选择预设值域(ρ1,ρ2)中点ρ3,确定η(ρ3),利用函数的单调性,当η(ρ3)*η(ρ1)<0时,则判断ρ0在(ρ1,ρ3)内,当η(ρ3)*η(ρ2)<0时,则判断ρ0在(ρ3,ρ2)内,如此循环,直到η(ρ3)小于设定的精度ε时,则此时的ρ3即为对应观测值Hz0值所对应的视电阻率值ρ0。
优选地,所述计算模块包括:
探测深度计算单元,用于利用公式或者确定所述探测深度。
为了解决上述技术问题,本发明还提供了一种瞬变电磁探测系统,包括上述的瞬变电磁探测装置和第一回线,所述第一回线的边长小于长度阈值。
与现有技术相比,本发明的方法和装置,从探测深度和磁矩两个方面对小回线和中心回线的等效性进行对比分析,采用中心回线资料的处理和解释方法代替以往的依据地质条件制作校正模型的方法,提高了资料的解释和精度和成果的可靠性。
附图说明
图1为本发明实施例的一种瞬变电磁探测方法的流程图;
图2为本发明实施例的一种瞬变电磁探测装置的结构示意图;
图3为本发明实施例的一种瞬变电磁探测的等效转换示意图;
图4为本发明实施例一的发射回线1的等效电场分布示意图;
图5为本发明实施例一的发射回线2的等效电场分布示意图;
图6为本发明实施例二的瞬变电磁探测的等效转换感应电压对比示意图;
图7为本发明实施例三的瞬变电磁探测数据的曲线图;
图8为本发明实施例三的瞬变电磁探测感应电压的衰减曲线图;
图9为本发明实施例三的瞬变电磁探测视电阻率等值曲线图;
图10为本发明实施例三的瞬变电磁探测视电阻率等值示意图;
图11为本发明实施例三的瞬变电磁探测视电阻率异常区域的示意图;
图12为本发明实施例三的瞬变电磁探测视电阻率等值线切片的示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下文中将结合附图对本发明的实施例进行详细说明。需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互任意组合。
如图1所示,本发明实施例提供一种瞬变电磁探测方法,包括:
A、采集待评估区域第一回线的瞬变电磁数据;所述第一回线的边长小于长度阈值;
B、将所述的第一回线的瞬变电磁数据等效转换成第二回线中心方式瞬变电磁数据;所述第二回线的边长大于或者等于长度阈值;
C、按照所述第二回线中心方式瞬变电磁数据的解释方法,确定待评估区域的视电阻率和探测深度。
当复杂山区工作时,很难布置足够大的发射回路进行传统的测量。常规大回线(第二回线)方法难以开展。本发明实施例把边长为3米的小线圈(第一回线)大电流瞬变电磁装置引入到山区进行探测。即向边长为3米的回线中通入1000A至2000A的大电流的,这种探测可以减少地形影响,并减小探测时的体积效应。但由于大电流问题的影响,到目前为至,小回线大电流瞬变电磁资料的处理还没有形成成熟的处理体系,本发明实施例把第一回线的瞬变电磁数据等效转换成第二回线中心方式瞬变电磁数据,按照中心方式瞬变电磁法的数据处理和解释方法对数据进行处理解释。本发明实施例把边长3米的小线圈1000A-2000A大电流瞬变电磁装置引入到山区进行探测,可以减少地形影响,并减小探测的体积效应。
其中,如图2所示,步骤B包括:
将所述的第一回线的瞬变电磁数据转换成相同磁矩大小的第二回线中心方式瞬变电磁数据。
小回线(第一回线)和中心回线(第二回线)具有磁矩和探测深度的等效性,小回线(第一回线)对应探测方法的能够得出和大回线(第二回线)探测方法相似的探测深度;在小线圈瞬变电磁数据处理时,首先把小线圈大电流瞬变电磁数据转换成相同磁矩大小的大回线小电流等效中心回线瞬变电磁数据,然后按照常规中心回线瞬变电磁资料解释方法进行解释。本发明实施例探测方法的结果和解释还可用于探测隐伏的煤层缝隙和寻找充水地区。
步骤C中确定待评估区域的视电阻率包括:
根据中心回线的视电阻率计算公式确定所述待评估区域的视电阻率:
所述中心回线的视电阻率计算公式为
其中,M是源发射磁矩,大小为发射电流与发射回线面积的乘积;q代表观测位置探头的有效接收面积,V(t)是测量的二次感应电压值,t是观测时间,μ0为自由空间的磁导率;
或者,采用二分法确定全时域视电阻率。
小回线(第一回线)对应探测方法为收发一体,实际观测位置为第一回线的中心点,因此,小回线装置的视电阻率计算可以采用第二回线的视电阻率计算公式,将第二回线公式中的发射磁矩用等效的小回线装置代替。
中心回线/重叠回线/第二回线的视电阻率晚期公式
其中,M是源发射磁矩,大小为发射电流与发射回线面积的乘积;q代表观测位置探头的有效接收面积,V(t)是测量的二次感应电压值,t是观测时间,μ0为自由空间的磁导率。
采用上述计算瞬变电磁的视电阻率对于晚期数据的处理是正确的,而对于瞬变电磁早期资料,采用晚期视电阻率会带来较大的计算误差,往往需要引入瞬变电磁早期信号的校正公式进行处理。本发明实施例还提出适用于中心回线的全程视电阻率计算方式:用迭代的方法:采用二分法确定全时域视电阻率包括:
将实际观测的二次感应电压转换成的垂直磁场,预设单调函数Hz=Y(ρ),并确定视电阻率ρ的预设值域(ρ1,ρ2),确定对应观测值Hz0值的反函数ρ0,包括:设定函数η(ρ)=Hz0-Y(ρ),选择预设值域(ρ1,ρ2)中点ρ3,确定η(ρ3),利用函数的单调性,当η(ρ3)*η(ρ1)<0时,则判断ρ0在(ρ1,ρ3)内,当η(ρ3)*η(ρ2)<0时,则判断ρ0在(ρ3,ρ2)内,如此循环,直到η(ρ3)小于设定的精度ε时,则此时的ρ3即为对应观测值Hz0值所对应的视电阻率值ρ0。用这种方法可以算出任意有限区间内的单调函数的反函数值,精度人为控制,程序易于在微机上实现,可用于实际资料处理。
步骤C中确定待评估区域的探测深度包括:
计算出视电阻率后,按照公式计算探测深度。
探测深度是任意时间瞬变电场的最大值作为阶跃脉冲到达的时间对应的深度,即场到达的深度。对于理想的观测情形(均匀介质无干扰),计算探测深度是有效的,而对于实际观测的大部分情况,计算探测深度都是不准确的。为了提高深度计算的有效性,本发明实施例提出一种计算探测深度(有效深度)的方法,场的传播和反射所需要时间是一致的.所以对地面T E M探测方法而言,要探测到地下一定深度的异常体,场需要用时间t到达该异常体然后反射回地面,即需要至少2t的时间将地下的地质信息带到地面被接收仪器所观察到,因此,探测深度(有效深度)表示为:
具体地,采用中心回线资料的处理和解释方法提高了资料的解释效率和精度;采用中心回线全程视电阻率,解决了传统定义引起的早期视电阻率计算精度较差的问题,提高了浅层目标体的探测精度;提出有效深度的深度计算方式,进一步提高了小回线大电流瞬变电磁的探测精度。
如图2所示,本发明实施例提供一种瞬变电磁探测装置,包括:
采集模块,用于采集待评估区域第一回线的瞬变电磁数据;所述第一回线的边长小于长度阈值;
等效模块,用于将所述的第一回线的瞬变电磁数据等效转换成第二回线中心方式瞬变电磁数据;所述第二回线的边长大于或者等于长度阈值;
计算模块,按照所述第二回线中心方式瞬变电磁数据的解释方法,确定待评估区域的视电阻率和探测深度。
其中,所述等效模块用于:
将所述的第一回线的瞬变电磁数据等效转换成相同磁矩大小的第二回线中心方式瞬变电磁数据。
所述计算模块包括:
探测深度计算单元,用于利用公式或者确定所述探测深度;
视电阻率计算单元,用于根据中心回线的视电阻率计算公式确定所述待评估区域的视电阻率:
所述中心回线的视电阻率计算公式为
其中,M是源发射磁矩,大小为发射电流与发射回线面积的乘积;q代表观测位置探头的有效接收面积,V(t)是测量的二次感应电压值,t是观测时间,μ0为自由空间的磁导率;
或者,采用二分法确定全时域视电阻率。
视电阻率计算单元采用二分法确定全时域视电阻率包括:
将实际观测的二次感应电压转换成的垂直磁场,预设单调函数Hz=Y(ρ),并确定视电阻率ρ的预设值域(ρ1,ρ2),确定对应观测值Hz0值的反函数ρ0,包括:设定函数η(ρ)=Hz0-Y(ρ),选择预设值域(ρ1,ρ2)中点ρ3,确定η(ρ3),利用函数的单调性,当η(ρ3)*η(ρ1)<0时,则判断ρ0在(ρ1,ρ3)内,当η(ρ3)*η(ρ2)<0时,则判断ρ0在(ρ3,ρ2)内,如此循环,直到η(ρ3)小于设定的精度ε时,则此时的ρ3即为对应观测值Hz0值所对应的视电阻率值ρ0。
本发明实施例提供一种瞬变电磁探测系统,包括上述瞬变电磁探测装置和第一回线,所述第一回线的边长小于长度阈值。
实施例一
小线圈和中心回线的探测深度等效性验证
时域电磁探测方法证明时间域TEM方法的探测深度主要由观测时长决定,同时与信噪比和接收机灵敏度有关。回线的边长与发射电流决定信号的强度,只要信号有足够的强度,那么对于小回线TEM方法的探测深度,在保证信噪比和接收机灵敏度的前提下,有和大回线同样的探测深度和有效性。在地面障碍无密集、场地狭小、和山地施工中,小回线表现出了比大回线更好的适应性。
图4和图5分别是发射回线边长为600米和10米,相同时刻的地下电场的数值模拟结果。可见大回线和小回线的场在同样的时间内扩散到了同样的深度。在深度1000米左右处,大、小回线激发的电场强度分别约为3V/m和0.05V/m。如果可以得到较大发射电流,小回线探测能够达到与大回线相当的探测深度。
实施例二
小线圈和中心回线的磁矩等效性说明
为了对两种回线的等效性进行定量分析,对于具有等效磁矩的中心回线和小回线进行对比工作,在同一点进行感应电压的观测工作,观测结果如图6所示。图中菱形曲线代表传统中心回线(回线边长100米,发送电流1A)观测的感应电压曲线,方形曲线为具有等效磁矩的小回线大电流(回线边长3米,发送电流1111A)观测的感应电压。在早期,小回线观测的响应大于中心回线的响应值,这是由两种观测装置的关断时间的不同引起的,小回线装置关断时间更早,响应也更强。曲线对比结果说明,可以把小回线大电流探测数据等效转换成相同磁矩的大回线小电流数据。并进一步按大回线中心方式瞬变电磁进行资料处理和解释。
实施例三
在本发明实施例三的待评估区域存在许多问题,因此,有必要探测地下富水区,以保证采煤过程的安全性。对于煤田富水区的探测,瞬变电磁方法是首选的地球物理方法。
本发明实施例三主要目的是探测出富水区在含煤地层中的位置。调查区地层主要包括奥陶系马家沟组、石炭系本溪太原组、二叠系山西下石河子组。岩性主要为砂岩,泥岩。煤层的电阻率范围为40-360Ω·m。石灰岩的电阻率较大(>500Ω·m)。由于在含煤地层的底板与顶板之间富水区分布不均匀,裂缝和断层的存在引起电性差异的变化,导致低电阻率异常的出现。
本发明实施例三的待评估区域为山区,平均海拔1387m,相对高度342m,地形条件复杂。在这一区域很难布置200m边长的线圈开展TEM勘探。因此,选择本发明实施例提供的小回线TEM系统(3m×3m)。
运用小回线TEM系统进行了煤矿充水区数据采集工作。采集的数据以V(t)/I表示,单位是μV/A。
数据采集的主要参数:边长为3m的发射线圈,发射频率为25Hz,时间窗的范围介于0.087-10ms之间,以及20个时间道的时间延时(表1)。
测网密度为40m×10m。共布置39条测线,线距40m,总长度13.7km,测点间距10m。
为了保证本发明实施例提供的小回线TEM系统的可靠性,在该区域勘探前,运用小回线TEM系统在同一位置反复采集瞬变测深数据。图7显示了在同一位置采集的五次TEM数据。衰减曲线几乎一致,这意味着小回线TEM系统TEM操作的稳定性。
图8是第一待评估区域衰减曲线。其中虚线显示衰减较慢,可以解释为地下含水丰富,实线衰减较快,可以解释为正常的层状介质。
第一待评估区域的视电阻率等值线以及相应的地质解释如图9所示。视电阻率计算采用二分法,探测深度计算采用公式
由图9可见,第一待评估区域存在三个相对低阻异常区(120-200Ω·m),充水区位于0m,300m-330m,和400m-460m处,在图中用虚线连接起来。异常带位于海拔1140m到1240m之间,倾角为10°。
根据瞬变电磁测深结果(图9),在第一待评估区域180m处进行打钻验证。钻孔深度为187.5m。打钻结果证实了在130.25m深度处存在一个宽度为5.1m的充水区。打钻结果和TEM测量结果一致。
为了对整个测区的富水采空区进行分析,图10绘出了海拔为1170m处的视电阻率等值线平面图。将视电阻率相对较低的区域定义为富水采空区。地质解释如11所示。在研究区域有四个相对较大的采空区,分别位于测区的中部,东北部和南部(阴影区域)。
图12描绘出了深度分别为-150m,-250m,和-350m的三个电阻率等值线切片。白色圆圈(低电阻率)代表富水区,浅灰色阴影区域(高电阻率)代表无采空的区域。在三处富水区进行打钻,所有的钻探结果和TEM解释结果相一致。
以上实施例仅用以说明本发明的技术方案而非限制,仅仅参照较佳实施例对本发明进行了详细说明。本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,均应涵盖在本发明的权利要求范围当中。