CN104793268A - 一种瞬变电磁探测的盲深度测量方法和装置 - Google Patents

一种瞬变电磁探测的盲深度测量方法和装置 Download PDF

Info

Publication number
CN104793268A
CN104793268A CN201510210281.8A CN201510210281A CN104793268A CN 104793268 A CN104793268 A CN 104793268A CN 201510210281 A CN201510210281 A CN 201510210281A CN 104793268 A CN104793268 A CN 104793268A
Authority
CN
China
Prior art keywords
mrow
self
signal voltage
msub
transient electromagnetic
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.)
Granted
Application number
CN201510210281.8A
Other languages
English (en)
Other versions
CN104793268B (zh
Inventor
薛国强
钟华森
崔江伟
李海
底青云
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Geology and Geophysics of CAS
Original Assignee
Institute of Geology and Geophysics of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Geology and Geophysics of CAS filed Critical Institute of Geology and Geophysics of CAS
Priority to CN201510210281.8A priority Critical patent/CN104793268B/zh
Publication of CN104793268A publication Critical patent/CN104793268A/zh
Application granted granted Critical
Publication of CN104793268B publication Critical patent/CN104793268B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本申请涉及一种瞬变电磁探测的盲深度测量方法和装置,盲深度测量方法包括:获取或接收通过在待测区域进行瞬变电磁探测所得到的多个时刻对应的实测信号电压幅值;计算瞬变电磁探测所采用的发射线圈在多个时刻对应的自感信号电压幅值;根据多个时刻对应的自感信号电压幅值和实测信号电压幅值,在多个时刻之中确定在待测区域进行瞬变电磁探测的盲深度的最小可分辨时间;以及根据最小可分辨时间及其对应的实测信号电压幅值,计算在待测区域进行瞬变电磁探测的盲深度。本发明的盲深度测量方法和装置,能够提高瞬变电磁探测盲深度的识别准确度。

Description

一种瞬变电磁探测的盲深度测量方法和装置
技术领域
本发明涉及电磁技术领域,具体涉及一种瞬变电磁探测的盲深度测量方法和装置。
背景技术
近年来,瞬变电磁法在金属矿勘探、工程地质、煤矿采空区勘察等领域的应用越来越广泛,取得了良好的经济效益。
瞬变电磁法(Transient electromagnetic method,TEM)是一种时间域电磁勘探方法,一般采用不接地的矩形或圆形回线发射阶跃电流,电流关断后,在地下介质中会感应出二次场,通过测量二次场随时间的变化而获知地质信息。由于瞬变电磁法在一次场关断后采集纯二次场,采集到的信号不受一次场干扰,因此可以在近区进行观测。此外,瞬变电磁法具有受一次场干扰小、施工方便、受地形影响小等特点。
像其它电法勘探方法一样,瞬变电磁探测深度问题一直是地球物理勘探的研究内容之一。研究探测深度的重要性是不言而喻的,因为无论是施工设计还是资料解释,都需要估算所用的方法、装置和仪器等能否达到目的层,以完成地质任务;对于野外作业中观测到的不同极距、不同频率、不同采样时间的数据,需要确定它们包含了哪一深度范围内的地质信息,以便给出正确的解释结果。
通常所说的瞬变电磁探测深度,一般是指瞬变电磁法的最大探测深度,即在仪器可分辨的前提下,电磁场的最大扩散深度。当一定延迟时间的信号强度太低,导致仪器无法分辨时,该时刻之后的数据就失去了意义,探测深度就停留在该时刻所对应的探测深度。
然而,由于发射线圈本身存在自感,使得接收器所接收的信号,一种是来自地下的电磁感应信号,一种是线圈本身的自感信号。在瞬变电磁的超早期,电磁感应信号幅值很小,而自感信号幅值很大,在通常情况下,自感信号要比电磁感应信号大三个数量级,所以自感信号会淹没电磁感应信号。并且电磁感应信号常常衰减很慢,而线圈的自感信号由高值向低值却衰减很快。所以,瞬变电磁有一个探测盲深度,该深度之内的地质信息是无法获知的。目前用于识别瞬变电磁探测盲深度的方法的识别准确度较低。
发明内容
本发明提供一种瞬变电磁探测的盲深度测量方法和装置,其目的是解决现有的用于识别瞬变电磁探测盲深度的技术所识别的盲深度的准确度较低的问题。
为了实现上述发明目的,本发明采取的技术方案如下:
一种瞬变电磁探测的盲深度测量方法,该方法包括:获取或接收通过在待测区域进行瞬变电磁探测所得到的多个时刻对应的实测信号电压幅值;计算瞬变电磁探测所采用的发射线圈在多个时刻对应的自感信号电压幅值;根据多个时刻对应的自感信号电压幅值和实测信号电压幅值,在多个时刻之中确定在待测区域进行瞬变电磁探测的盲深度的最小可分辨时间;以及根据最小可分辨时间及其对应的实测信号电压幅值,计算在待测区域进行瞬变电磁探测的盲深度。
优选地,最小可分辨时间通过如下方式确定:针对多个时刻中的每一个时刻,用该时刻对应的实测信号电压幅值减去该时刻对应的自感信号电压幅值后所得的值,作为该时刻对应的有用电压幅值;计算多个时刻中每个时刻对应的有用电压幅值与该时刻对应的自感信号电压幅值之比值;以及在多个时刻中,将比值与预定值最接近的时刻确定为在待测区域进行瞬变电磁探测的盲深度的最小可分辨时间。
优选地,预定值为1.5。
优选地,计算瞬变电磁探测所采用的发射线圈在多个时刻对应的自感信号电压幅值的步骤包括:针对多个时刻中的每一个时刻,根据公式UL(t)=2IRmte-mt,计算发射线圈中的每个单匝线圈在该时刻的自感信号电压幅值,以通过对该时刻的各单匝线圈的自感信号电压幅值进行叠加来获得发射线圈在该时刻对应的自感信号电压幅值;其中,t为时间,UL(t)表示t时刻的自感信号电压,I为发射电流,R为发射线圈的电阻;以及 L = L ex + L in , L in = μ 8 π × 4 a = μa 2 π , L ex = 2 μa π , 其中,Lin为内自感强度,Lex为外自感强度,单位为mL/H,a为线圈边长,C为导线电容,C=4LC1,C1为每米长度的线圈电容,μ=μ0,μr≈μ0,μ0为真空磁导率。
优选地,在待测区域进行瞬变电磁探测的盲深度根据如下公式获得: δ LTD = 12.64 × 10 - 3 a 4 / 3 q 2 / 3 [ V ( t d ) / I ] - 2 / 3 t d - 2 / 3 μ 0 , 其中,δTD表示盲深度,td为最小可分辨时间,q为进行瞬变电磁探测的接收器探头的有效接收面积,V(td)为td时刻对应的实测信号电压幅值,a为线圈边长,I为发射电流,μ0为真空磁导率。
一种瞬变电磁探测的盲深度测量装置,该装置包括:感应电压获得单元,其用于获取或接收通过在待测区域进行瞬变电磁探测所得到的多个时刻对应的实测信号电压幅值;自感电压计算单元,其用于计算瞬变电磁探测所采用的发射线圈在多个时刻对应的自感信号电压幅值;最小可分辨时间确定单元,其用于根据多个时刻对应的自感信号电压幅值和实测信号电压幅值,在多个时刻之中确定在待测区域进行瞬变电磁探测的盲深度的最小可分辨时间;以及盲深度计算单元,其用于根据最小可分辨时间及其对应的实测信号电压幅值,计算在待测区域进行瞬变电磁探测的盲深度。
优选地,最小可分辨时间确定单元用于:针对多个时刻中的每一个时刻,用该时刻对应的实测信号电压幅值减去该时刻对应的自感信号电压幅值后所得的值,作为该时刻对应的有用电压幅值;计算多个时刻中每个时刻对应的有用电压幅值与该时刻对应的自感信号电压幅值之比值;以及在多个时刻中,将比值与预定值最接近的时刻确定为在待测区域进行瞬变电磁探测的盲深度的最小可分辨时间。
优选地,预定值为1.5。
优选地,自感电压计算单元用于:针对多个时刻中的每一个时刻,根据公式UL(t)=2IRmte-mt,计算发射线圈中的每个单匝线圈在该时刻的自感信号电压幅值,以通过对该时刻的各单匝线圈的自感信号电压幅值进行叠加来获得发射线圈在该时刻对应的自感信号电压幅值;其中,t为时间,UL(t)表示t时刻的自感信号电压,I为发射电流,R为发射线圈的电阻;以及 L = L ex + L in , L in = μ 8 π × 4 a = μa 2 π , L ex = 2 μa π , 其中,Lin为内自感强度,Lex为外自感强度,单位为mL/H,a为线圈边长,C为导线电容,C=4LC1,C1为每米长度的线圈电容,μ=μ0,μr≈μ0,μ0为真空磁导率。
优选地,盲深度计算单元用于根据如下公式获得在待测区域进行瞬变电磁探测的盲深度: δ LTD = 12.64 × 10 - 3 a 4 / 3 q 2 / 3 [ V ( t d ) / I ] - 2 / 3 t d - 2 / 3 μ 0 , 其中,δTD表示盲深度,td为最小可分辨时间,q为进行瞬变电磁探测的接收器探头的有效接收面积,V(td)为td时刻对应的实测信号电压幅值,a为线圈边长,I为发射电流,μ0为真空磁导率。
本发明和现有技术相比,具有如下有益效果:
本发明的一种瞬变电磁探测的盲深度测量方法和装置,能够通过对比线圈自感信号的计算值与野外实际测试信号值,获得瞬变电磁识别盲深度的最小可分辨时间,进而根据该最小可分辨时间来获得瞬变电磁探测盲深度大小。相比于现有技术,利用该方法或装置所识别的瞬变电磁探测盲深度的准确度较高,对于实际应用意义重大。利用本发明的上述技术所获得的瞬变电磁探测盲深度,其准确度能够达到80%以上,相比于现有技术能够大大提高识别准确度。
附图说明
图1是本发明实施例的一种瞬变电磁探测的盲深度测量方法的一个示例的流程图;
图2是野外回线源测量装备示意图;
图3是本发明的实施例的一种瞬变电磁探测的盲深度测量装置的框图。
具体实施方式
为使本发明的发明目的、技术方案和有益效果更加清楚明了,下面结合附图对本发明的实施例进行说明,需要说明的是,在不冲突的情况下,本申请中的实施例和实施例中的特征可以相互任意组合。
本发明的实施例提供了一种瞬变电磁探测的盲深度测量方法,该方法包括:获取或接收通过在待测区域进行瞬变电磁探测所得到的多个时刻对应的实测信号电压幅值;计算瞬变电磁探测所采用的发射线圈在多个时刻对应的自感信号电压幅值;根据多个时刻对应的自感信号电压幅值和实测信号电压幅值,在多个时刻之中确定在待测区域进行瞬变电磁探测的盲深度的最小可分辨时间;以及根据最小可分辨时间及其对应的实测信号电压幅值,计算在待测区域进行瞬变电磁探测的盲深度。
图1示出了根据本发明实施例的一种瞬变电磁探测的盲深度测量方法的一个示例处理的流程图。如图1所示,该处理流程开始之后,首先执行步骤S110。
在步骤S110中,获取或接收通过在待测区域进行瞬变电磁探测所得到的多个时刻对应的实测信号电压幅值。然后,执行步骤S120。
其中,步骤S110中获取或接收的实测信号电压幅值可以通过如图2所示的实际测量设备来测得。
如图2所示,在地面上布设一个矩形的发射回线线圈(即步骤S120中所述的发射线圈),在发射回线正中心放置一个接收磁探头(即下文所述的接收器探头),然后向发射回线内输入电流。在电流关断的瞬间,利用接收磁探头采集地下介质感应出的二次感应电压信号,作为实测信号电压幅值,完成感应电压野外实际测量工作。其中,实测信号电压幅值是在断电(即电流关断)后预定的多个时刻测量的,也即,实测信号电压幅值是包括多个测量值的一组值。
在步骤S120中,计算瞬变电磁探测所采用的发射线圈在多个时刻对应的自感信号电压幅值。然后,执行步骤S130。
优选地,在步骤S120中,可以根据实际工作中发射线材的特性来计算瞬变电磁探测所采用的发射线圈在断电后的多个时刻对应的自感信号电压幅值。例如,针对多个时刻中的每一个时刻,可以根据如下描述的公式一来计算发射线圈中的每个单匝线圈在该时刻的自感信号电压幅值,然后,通过对该时刻的各单匝线圈的自感信号电压幅值进行叠加,将叠加后得到的值作为发射线圈在该时刻对应的自感信号电压幅值。
公式一:UL(t)=2IRmte-mt
其中,t为时间,UL(t)表示t时刻的自感信号电压,I为发射电流,R为发射线圈的电阻。
其中,上述公式一的各参数可以根据公式二至公式五来获得。
公式二: m = 1 / LC .
公式三:L=Lex+Lin
公式四: L in = μ 8 π × 4 a = μa 2 π .
公式五: L ex = 2 μa π .
其中,Lin为内自感强度,Lex为外自感强度,单位为mL/H,a为线圈边长,C为导线电容,C=4LC1,C1为每米长度的线圈电容,μ=μ0,μr≈μ0,μ0为真空磁导率。
由于线圈与大地是绝缘的,因此在公式四和公式五中,磁导率近似等于真空中的磁导率。
在步骤S130中,根据多个时刻对应的自感信号电压幅值和实测信号电压幅值,在多个时刻之中确定在待测区域进行瞬变电磁探测的盲深度的最小可分辨时间。然后,执行步骤S140。
优选地,在步骤S130中,最小可分辨时间通过如下方式确定:针对多个时刻中的每一个时刻,用该时刻对应的实测信号电压幅值减去该时刻对应的自感信号电压幅值后所得的值,作为该时刻对应的有用电压幅值;然后,针对多个时刻中的每一个时刻,计算该时刻对应的有用电压幅值与该时刻对应的自感信号电压幅值之比值;在多个时刻对应的比值之中,选择与预定值最接近的比值,然后将该比值对应的时刻确定为在待测区域进行瞬变电磁探测的盲深度的最小可分辨时间。
例如,用UL(t)表示t时刻对应的自感信号电压幅值,V(t)表示实测信号电压幅值,UV(t)表示t时刻对应的有用电压,则UV(t)=V(t)-UL(t)。其中,t时刻可以是上述多个时刻中的任一个时刻。这样,通过计算获得t时刻对应的UV(t)/UL(t)比值,进而可以在多个时刻对应的UV(t)/UL(t)比值之中选择与预定值最接近的那个比值,然后将该比值对应的时刻作为上述最小可分辨时间。
在一个优选例子中,上述预定值可以为1.5。经过实验验证,比值大于等于时,能够明显识别出有用信息。此外,最小可分辨时间也会随着该比值的增大而增大,进而使得探测的盲深度随着该比值的增大而增大,当选择较大的比值时,盲深度过大将对瞬变电磁探测不利。这样,当预定值选取1.5时,通过比对每个时刻(即每个探测时间点)的有用电压与实测信号电压的幅值之比值,将比值最接近1.5的那个时刻确定为上述最小可分辨时间,由此识别的最小可分辨时间的准确度非常高。经过实验验证,当预定值选取1.5时,后续获得的探测盲深度的准确度较高。
实际测量的电压信号由地下介质感应出的有用电压UV(t)和线圈的自感电压UL(t)构成。当有用电压幅值约为自感电压幅值的1.5倍时,对应时间为最小可分辨时间,其对应的深度为最小探测深度。这样,在考虑了仪器本身灵敏度的情况下,当有用电压UV(t)大于(特别是远远大于)自感电压UL(t)时,瞬变电磁探测设备(例如图2所示的设备)能够较好地识别出有用电压,从而较准确地确定出地下目标体的位置。
在步骤S140中,根据最小可分辨时间及其对应的实测信号电压幅值,计算在待测区域进行瞬变电磁探测的盲深度。然后,结束处理。
优选地,在待测区域进行瞬变电磁探测的盲深度根据如下公式六获得。
公式六: δ LTD = 12.64 × 10 - 3 a 4 / 3 q 2 / 3 [ V ( t d ) / I ] - 2 / 3 t d - 2 / 3 μ 0 .
其中,δTD表示盲深度,td为最小可分辨时间,q为进行瞬变电磁探测的接收器探头的有效接收面积,V(td)为td时刻对应的实测信号电压幅值。
经过实验验证,利用上述公式六所获得的瞬变电磁探测盲深度,其准确度能够达到80%以上,相比于现有技术能够大大提高瞬变电磁探测盲深度的识别准确度。
本发明的实施例还提供了一种瞬变电磁探测的盲深度测量装置,如图3所示,该装置包括:感应电压获得单元1,其用于获取或接收通过在待测区域进行瞬变电磁探测所得到的多个时刻对应的实测信号电压幅值;自感电压计算单元2,其用于计算瞬变电磁探测所采用的发射线圈在多个时刻对应的自感信号电压幅值;最小可分辨时间确定单元3,其用于根据多个时刻对应的自感信号电压幅值和实测信号电压幅值,在多个时刻之中确定在待测区域进行瞬变电磁探测的盲深度的最小可分辨时间;以及盲深度计算单元4,其用于根据最小可分辨时间及其对应的实测信号电压幅值,计算在待测区域进行瞬变电磁探测的盲深度。
优选地,最小可分辨时间确定单元3可以用于:针对多个时刻中的每一个时刻,用该时刻对应的实测信号电压幅值减去该时刻对应的自感信号电压幅值后所得的值,作为该时刻对应的有用电压幅值;计算多个时刻中每个时刻对应的有用电压幅值与该时刻对应的自感信号电压幅值之比值;以及在多个时刻中,将比值与预定值最接近的时刻确定为在待测区域进行瞬变电磁探测的盲深度的最小可分辨时间。
优选地,自感电压计算单元2可以用于:针对多个时刻中的每一个时刻,根据上文所述的公式一来计算发射线圈中的每个单匝线圈在该时刻的自感信号电压幅值,以通过对该时刻的各单匝线圈的自感信号电压幅值进行叠加来获得发射线圈在该时刻对应的自感信号电压幅值;其中,公式一中各参数意义和表达与上文相同,这里不再重述。
优选地,盲深度计算单元4可以用于根据上文所述的公式六来获得在待测区域进行瞬变电磁探测的盲深度;其中,公式六中各参数意义和表达与上文相同,这里不再重述。
需要说明的是,感应电压获得单元1、自感电压计算单元2、最小可分辨时间确定单元3和盲深度计算单元4可以分别对应地执行上文所描述的一种瞬变电磁探测的盲深度测量方法的步骤S110-S140中的处理,并可以达到相类似的功能和效果,这里不再详述。
通过以上描述可知,上述根据本发明实施例的一种瞬变电磁探测的盲深度测量方法或装置,其通过对比线圈自感信号的计算值与野外实际测试信号值,获得瞬变电磁识别盲深度的最小可分辨时间,进而根据该最小可分辨时间来获得瞬变电磁探测盲深度大小。相比于现有技术,利用该方法或装置所识别的瞬变电磁探测盲深度的准确度较高,对于实际应用意义重大。
下面描述本发明实施例的一种瞬变电磁探测的盲深度测量方法的应用示例。
采用图2所示测量方式,按照本发明实施例所描述的一种瞬变电磁探测的盲深度测量方法在某矿区(作为待测区域的示例)进行瞬变电磁探测。
其中,选用中心回线方式进行探测,采用600m×600m回线框进行发射,发射电流为6安培,采用较常用的线材:导线横截面半径r0=0.73mm(作为发射线圈的示例),其电阻为:0.3Ω/100m,等效电容C为:C=4LC1,C1为每米长度的线圈电容,计算中取C1=160×10-9F。
首先根据公式一至公式五计算断电后的自感信号电压幅值UL(t),计算结果见表1。这里,选择1μs、1.5μs、2μs、......、12.5μs、13μs作为上文所述的多个时刻。
表1不同时刻的自感信号电压幅值
表2是该矿区某测点实测瞬变电磁数据(即上文所述的实测信号电压幅值),采用与表1中数字计算所用的材料相同的发射线材进行信号发射。探测时,采用SB70K磁探头进行数据采集,有效接收面积2000m2,发射电流为6安培。
表2瞬变电磁实测数据
对比表1和表2可知,在瞬变早期,实测信号中自感信号占主导,随着时间的推移,自感信号在实测信号中所占的比例越来越小,在td=9.5μs时,有用信号明显占优势,约为自感信号的1.5倍,之后自感信号在实测信号中所占的比例越来越小,因此可以将td=9.5μs确定为最小可分辨时间。
这样,根据公式六,计算瞬变电磁探测的盲深度,将td=9.5μs和V(t)=0.056mV代入公式六,可得δLTD=47.95m,可知本次瞬变电磁探测盲深度为47.95m,小于该深度的地质信息无法获得。
虽然本发明所揭示的实施方式如上,但其内容只是为了便于理解本发明的技术方案而采用的实施方式,并非用于限定本发明。任何本发明所属技术领域内的技术人员,在不脱离本发明所揭示的核心技术方案的前提下,可以在实施的形式和细节上做任何修改与变化,但本发明所限定的保护范围,仍须以所附的权利要求书限定的范围为准。

Claims (10)

1.一种瞬变电磁探测的盲深度测量方法,其特征在于,所述盲深度测量方法包括:
获取或接收通过在待测区域进行瞬变电磁探测所得到的多个时刻对应的实测信号电压幅值;
计算所述瞬变电磁探测所采用的发射线圈在所述多个时刻对应的自感信号电压幅值;
根据所述多个时刻对应的所述自感信号电压幅值和所述实测信号电压幅值,在所述多个时刻之中确定在所述待测区域进行所述瞬变电磁探测的盲深度的最小可分辨时间;以及
根据所述最小可分辨时间及其对应的实测信号电压幅值,计算在所述待测区域进行瞬变电磁探测的盲深度。
2.根据权利要求1所述的盲深度测量方法,其特征在于,所述最小可分辨时间通过如下方式确定:
针对所述多个时刻中的每一个时刻,用该时刻对应的实测信号电压幅值减去该时刻对应的自感信号电压幅值后所得的值,作为该时刻对应的有用电压幅值;
计算所述多个时刻中每个时刻对应的有用电压幅值与该时刻对应的自感信号电压幅值之比值;以及
在所述多个时刻中,将比值与预定值最接近的时刻确定为在所述待测区域进行所述瞬变电磁探测的盲深度的最小可分辨时间。
3.根据权利要求2所述的盲深度测量方法,其特征在于,所述预定值为1.5。
4.根据权利要求1-3中任一项所述的盲深度测量方法,其特征在于,所述计算所述瞬变电磁探测所采用的发射线圈在所述多个时刻对应的自感信号电压幅值的步骤包括:
针对所述多个时刻中的每一个时刻,根据公式UL(t)=2IRmte-mt,计算所述发射线圈中的每个单匝线圈在该时刻的自感信号电压幅值,以通过对该时刻的各单匝线圈的自感信号电压幅值进行叠加来获得所述发射线圈在该时刻对应的自感信号电压幅值;
其中,t为时间,UL(t)表示t时刻的自感信号电压,I为发射电流,R为所述发射线圈的电阻;以及
m = 1 / LC ,
L=Lex+Lin
L in = μ 8 π × 4 a = μa 2 π ,
L ex = 2 μa π ,
其中,Lin为内自感强度,Lex为外自感强度,单位为mL/H,a为线圈边长,C为导线电容,C=4LC1,C1为每米长度的线圈电容,μ=μ0,μr≈μ0,μ0为真空磁导率。
5.根据权利要求1-3中任一项所述的盲深度测量方法,其特征在于,所述在所述待测区域进行瞬变电磁探测的盲深度根据如下公式获得:
δ LTD = 12.64 × 10 - 3 a 4 / 3 q 2 / 3 [ V ( t d ) / I ] - 2 / 3 t d - 2 / 3 μ 0 ,
其中,δTD表示所述盲深度,td为最小可分辨时间,q为进行瞬变电磁探测的接收器探头的有效接收面积,V(td)为td时刻对应的实测信号电压幅值,a为线圈边长,I为发射电流,μ0为真空磁导率。
6.一种瞬变电磁探测的盲深度测量装置,其特征在于,所述盲深度测量装置包括:
感应电压获得单元,其用于获取或接收通过在待测区域进行瞬变电磁探测所得到的多个时刻对应的实测信号电压幅值;
自感电压计算单元,其用于计算所述瞬变电磁探测所采用的发射线圈在所述多个时刻对应的自感信号电压幅值;
最小可分辨时间确定单元,其用于根据所述多个时刻对应的所述自感信号电压幅值和所述实测信号电压幅值,在所述多个时刻之中确定在所述待测区域进行所述瞬变电磁探测的盲深度的最小可分辨时间;以及
盲深度计算单元,其用于根据所述最小可分辨时间及其对应的实测信号电压幅值,计算在所述待测区域进行瞬变电磁探测的盲深度。
7.根据权利要求6所述的盲深度测量装置,其特征在于,所述最小可分辨时间确定单元用于:
针对所述多个时刻中的每一个时刻,用该时刻对应的实测信号电压幅值减去该时刻对应的自感信号电压幅值后所得的值,作为该时刻对应的有用电压幅值;
计算所述多个时刻中每个时刻对应的有用电压幅值与该时刻对应的自感信号电压幅值之比值;以及
在所述多个时刻中,将比值与预定值最接近的时刻确定为在所述待测区域进行所述瞬变电磁探测的盲深度的最小可分辨时间。
8.根据权利要求7所述的盲深度测量方法,其特征在于,所述预定值为1.5。
9.根据权利要求6-8中任一项所述的盲深度测量装置,其特征在于,所述自感电压计算单元用于:
针对所述多个时刻中的每一个时刻,根据公式UL(t)=2IRmte-mt,计算所述发射线圈中的每个单匝线圈在该时刻的自感信号电压幅值,以通过对该时刻的各单匝线圈的自感信号电压幅值进行叠加来获得所述发射线圈在该时刻对应的自感信号电压幅值;
其中,t为时间,UL(t)表示t时刻的自感信号电压,I为发射电流,R为所述发射线圈的电阻;以及
m = 1 / LC ,
L=Lex+Lin
L in = μ 8 π × 4 a = μa 2 π ,
L ex = 2 μa π ,
其中,Lin为内自感强度,Lex为外自感强度,单位为mL/H,a为线圈边长,C为导线电容,C=4LC1,C1为每米长度的线圈电容,μ=μ0,μr≈μ0,μ0为真空磁导率。
10.根据权利要求6-8中任一项所述的盲深度测量装置,其特征在于,所述盲深度计算单元用于根据如下公式获得在所述待测区域进行瞬变电磁探测的盲深度:
δ LTD = 12.64 × 10 - 3 a 4 / 3 q 2 / 3 [ V ( t d ) / I ] - 2 / 3 t d - 2 / 3 μ 0 ,
其中,δTD表示所述盲深度,td为最小可分辨时间,q为进行瞬变电磁探测的接收器探头的有效接收面积,V(td)为td时刻对应的实测信号电压幅值,a为线圈边长,I为发射电流,μ0为真空磁导率。
CN201510210281.8A 2015-04-28 2015-04-28 一种瞬变电磁探测的盲深度测量方法和装置 Active CN104793268B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510210281.8A CN104793268B (zh) 2015-04-28 2015-04-28 一种瞬变电磁探测的盲深度测量方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510210281.8A CN104793268B (zh) 2015-04-28 2015-04-28 一种瞬变电磁探测的盲深度测量方法和装置

Publications (2)

Publication Number Publication Date
CN104793268A true CN104793268A (zh) 2015-07-22
CN104793268B CN104793268B (zh) 2017-07-25

Family

ID=53558229

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510210281.8A Active CN104793268B (zh) 2015-04-28 2015-04-28 一种瞬变电磁探测的盲深度测量方法和装置

Country Status (1)

Country Link
CN (1) CN104793268B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105629314A (zh) * 2016-03-24 2016-06-01 中国科学院地质与地球物理研究所 一种瞬变电磁发射线圈的自感信号的剔除方法和装置
CN110133737A (zh) * 2019-06-26 2019-08-16 中国科学院地质与地球物理研究所 一种隐伏矿的电磁学预测方法
CN110220944A (zh) * 2019-05-16 2019-09-10 太原理工大学 应用于混凝土搅拌罐内料浆监测及称重装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5955884A (en) * 1994-08-15 1999-09-21 Western Atlas International, Inc. Method and apparatus for measuring transient electromagnetic and electrical energy components propagated in an earth formation
CN103675926A (zh) * 2012-09-24 2014-03-26 成都理工大学 一种航空瞬变电磁资料的电导率-深度转换方法
CN104502987A (zh) * 2014-12-09 2015-04-08 天津华勘地质勘查有限公司 物探瞬变电磁系统测量数据精细化处理方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5955884A (en) * 1994-08-15 1999-09-21 Western Atlas International, Inc. Method and apparatus for measuring transient electromagnetic and electrical energy components propagated in an earth formation
CN103675926A (zh) * 2012-09-24 2014-03-26 成都理工大学 一种航空瞬变电磁资料的电导率-深度转换方法
CN104502987A (zh) * 2014-12-09 2015-04-08 天津华勘地质勘查有限公司 物探瞬变电磁系统测量数据精细化处理方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
F. KAMENETSKY,ET AL.: "Distortions of EM transients in coincident loops at short time-delays", 《GEOPHYSICAL PROSPECTING》 *
于景邨等: "用瞬变电磁法探查综放工作面顶板水体的研究", 《中国矿业大学学报》 *
应征: "土石坝隐患瞬变电磁检测方法研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *
薛国强: "论瞬变电磁测深法的探测深度", 《石油地球物理勘探》 *
薛国强等: "大回线源瞬变电磁法最小探测深度的分析与估算", 《工程勘察》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105629314A (zh) * 2016-03-24 2016-06-01 中国科学院地质与地球物理研究所 一种瞬变电磁发射线圈的自感信号的剔除方法和装置
CN110220944A (zh) * 2019-05-16 2019-09-10 太原理工大学 应用于混凝土搅拌罐内料浆监测及称重装置
CN110133737A (zh) * 2019-06-26 2019-08-16 中国科学院地质与地球物理研究所 一种隐伏矿的电磁学预测方法

Also Published As

Publication number Publication date
CN104793268B (zh) 2017-07-25

Similar Documents

Publication Publication Date Title
CN105044792B (zh) 地-井时频电磁勘探数据采集装置及方法
CN101520517B (zh) 一种能准确评价碎屑岩盆地含油气目标的方法
CN101382599B (zh) 一种确定储层孔隙各向异性的瞬变电磁方法
CN104360401B (zh) 一种瞬变电磁b场确定地下目标体地质信息方法
CN112415615B (zh) 基于分布式光纤传感的时频电磁压裂监测系统及监测方法
CN109541695B (zh) 人工场源频率域电场梯度远区视电阻率快速成像方法
CN103852794A (zh) 烃类污染浅层地下水磁共振检测装置及检测方法
CN111538093A (zh) 一种用于浅层地表探测方法及瞬变电磁仪器
CN104020496A (zh) 一种轴向同线方式的地面可控源电磁勘探方法
CN104375197B (zh) 一种电磁探测方法及装置
CN104793268B (zh) 一种瞬变电磁探测的盲深度测量方法和装置
JP2939334B2 (ja) 可聴周波以下の低周波磁気測定装置
CN106199708B (zh) 一种城区含煤构造地球物理预测方法与装置
CN104749643A (zh) 一种用于电磁干扰地区的csamt测量方法和装置
Thiesson et al. Characterization of buried cables and pipes using electromagnetic induction loop-loop frequency-domain devices
CN104391335A (zh) 瞬变电磁纯异常信号提取方法
CN111290037B (zh) 磁共振无损原位探测地下有机污染物的装置及探测方法
Kang et al. Forward modeling and analyzing for frequency domain semi-airborne EM method
RU2733095C2 (ru) Способ поиска трехмерных объектов методами геоэлектрики тм-поляризации
CN115522914A (zh) 一种套后储层径向远距离高精度探测方法和系统
CN104749645B (zh) 一种煤层基底奥灰顶界面起伏形态探测的方法和装置
CN106610509A (zh) 一种瞬变电磁资料的时域处理方法
RU2231089C1 (ru) Способ геоэлектроразведки (варианты)
CA2809328A1 (en) Magnetic variation method
EA036449B1 (ru) Устройство для геоэлектрического профилирования почвенно-мерзлотного комплекса

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant