CN107024722B - 一种基于异常环的低温超导磁源瞬变电磁地形校正方法 - Google Patents
一种基于异常环的低温超导磁源瞬变电磁地形校正方法 Download PDFInfo
- Publication number
- CN107024722B CN107024722B CN201710256983.9A CN201710256983A CN107024722B CN 107024722 B CN107024722 B CN 107024722B CN 201710256983 A CN201710256983 A CN 201710256983A CN 107024722 B CN107024722 B CN 107024722B
- Authority
- CN
- China
- Prior art keywords
- abnormal ring
- ring
- low
- abnormal
- temperature superconducting
- 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/38—Processing data, e.g. for analysis, for interpretation, for correction
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
本发明公开了一种基于异常环的低温超导磁源瞬变电磁地形校正方法,其方法为:步骤一、在单倾斜山体进行地面磁源瞬变电磁实验;步骤二、发射机激发梯形波;步骤三、利用低温超导传感器在发射线圈中心点测量两次;步骤四、理论计算得到异常环的理论响应Btheory;步骤五、得到校正系数。本发明的有益效果:根据本发明提出的方法,实现地形校正,消除地形影响导致的大地电阻率解释误差,最终实现了起伏地形下的大地真实电阻率解释成像。
Description
技术领域
本发明涉及一种基于异常环的地形校正方法,特别涉及一种基于异常环的低温超导磁源瞬变电磁地形校正方法。
背景技术
目前,瞬变电磁法(Transient Electromagnetic Methods)是一种时间域的物探方法。它的测量原理是,介质在一次电流脉冲场激励下会产生涡流,在脉冲间断期间涡流不会立刻消失,其周围空间会形成随时间衰减的二次磁场。二次磁场随时间衰减规律主要取决于导电体的导电性、体积规模、埋深,以及发射电流的形态和频率,因此,可以通过接收系统测量的二次场空间分布形态来了解异常体的空间分布。
超导量子干涉器(superconducting quantum interference device,SQUID)直接测量磁场,噪声低(fT,10-15T,量级)、带宽大、低频响应特性好,可提升瞬变电磁法(TEM)晚期接收信号质量,实现大深度探测。2011年,中国科学院上海微系统与信息技术研究所、德国于利希超导研究所与吉林大学三方合作,开始研究将Low-Tc SQUID应用于地面瞬变电磁探测中。低温SQUID需要工作在液氦温区,无磁杜瓦可以储存液氦,为SQUID提供低温工作环境。
中国专利CN201510124047.3公开了一种由井出发的瞬变电磁勘探地形校正处理方法,首先对瞬变电磁采集的原始资料进行处理得到瞬变电磁数据集,再对瞬变电磁数据集进行高频、低频分量的分离,得到低频分量数据集;将井资料低频分量与实际观测资料低频分量进行匹配滤波,最后将匹配滤波后的低频分量与实际观测到的高频分量进行叠加完成地形校正。解决了现有的比值法和带地形的二维或三维反演方法将具有复杂地质结构的地层转换为均匀介质容易出现转换过程中的计算误差的问题。但是在瞬变电磁实际勘探矿山过程中,很多测区的地质资料不全面或者没有地质资料,所以没办法从地质资料中提取低频分量完成地形校正。
中国专利CN201610139371.7公布了一种消除瞬变电磁观测数据地形影响的方法和装置,利用瞬变电磁二次感应电压实测数据和观测参数确定观测点对应的视电阻率数据,从而确定视电阻率随时间变化的斜率曲线;根据观测点对应视电阻率随时间变化的斜率曲线值确定获得消除地形影响后的斜率数据等值线剖面图,消除了地形给瞬变电磁数据带来的影响。
以上所述方法公布了电磁勘探中地形校正的方法,但均未涉及低温超导磁传感器的测量数据的地形校正。在低温超导传感器进行测量时,杜瓦必须竖直放置于发射线圈中心点,则在测区地形不平缓的情况下,接收机接收到的是实际磁场响应的垂直分量,而实际磁场响应的方向是垂直于倾斜地形的,所以直接利用实测数据进行视电阻率成像,是不准确的,因此,有必要对实测数据进行地形校正。
发明内容
本发明的目的是为了解决在现有的电磁勘探的地形校正方法中,均未涉及低温超导磁传感器的测量数据的地形校正的问题,而提供的一种基于异常环的低温超导磁源瞬变电磁地形校正方法。
本发明提供的基于异常环的低温超导磁源瞬变电磁地形校正方法,其方法如下所述:
步骤一、在单倾斜山体进行地面磁源瞬变电磁实验,采用中心方式,发射线圈依山体铺设成边长为a的正方形,低温超导传感器保存在真空杜瓦盛装的液氦里,为确保实验安全,杜瓦必须竖直放置于发射线圈中心点,不随山体倾斜角度变化;
步骤二、发射机激发梯形波,利用采集器记录发射波形得到关断时间toff,发射电流为I;
步骤三、利用低温超导传感器在发射线圈中心点测量两次,分别是中心点有异常环和无异常环时的情况,测量的响应分别是Btotal和Bground,Btotal同时包含了大地的响应Bground和异常环的响应Banomaly,即Banomaly=Btotal-Bground;
步骤四、理论计算得到异常环的理论响应Btheory;异常环等效为一个电阻和一个电感串联,当发射电流为梯形波时,得到等效回路中的感应电流表达式如下:
t<0 i(t)=0;
上式中,τ是异常环的时间常数,是等效回路切割一次场的磁通量,toff是发射电流的关断时间,L是异常环的等效电感;Btheory结果如下式:
式中,μ0是真空磁导率,n是异常环的匝数,r是异常环的半径;
步骤五、得到校正系数cosθ=Banomaly/Btheory,该系数用于大地的电磁数据校正,即把异常环的实测数据除以cosθ,得到与实际地形垂直方向上的中心磁场响应。
本发明的工作原理:
在起伏地形下采用低温超导传感器,首先测量大地的磁场响应垂直分量,然后放置已知电性参数的异常环,再测量含有异常环和大地的磁场响应垂直分量,将两次测量数据进行作差,提取异常环的磁场响应垂直分量;然后通过理论计算获得异常环的总磁场响应,最后计算磁场总响应和垂直分量的比值,进而间接获得了大地倾角θ;结合磁场响应之间的几何关系,利用大地的倾角参数对实测数据进行地形校正,从而消除地形影响导致的大地电阻率解释误差,最终实现了起伏地形下的大地真实电阻率解释成像。
本发明的有益效果:
本发明针对高阻地区、低温超导磁场传感器测量系统的数据,在单斜地形或者发射线圈铺设的地形条件是单斜的测区引起的数据误差,提出了只需在野外测区测量大地磁场响应的同时,在发射线圈的中心点放置已知电性参数的异常环,测量异常环和大地的总磁场响应,得到两组测量数据,再根据本发明提出的方法,实现地形校正,消除地形影响导致的大地电阻率解释误差,最终实现了起伏地形下的大地真实电阻率解释成像。
附图说明
图1为本发明所述方法流程示意图。
图2为本发明所述TEM野外实验实施配置示意图。
图3为本发明所述异常环等效回路示意图。
图4为本发明所述单斜地形条件下TEM测量数据的几何校正原理示意图。
图5为本发明实施例测区磁场传感器测量系统的数据地形校正效果示意图。
具体实施方式
请参阅图1、图2、图3、图4和图5所示:
本发明提供的基于异常环的低温超导磁源瞬变电磁地形校正方法,其方法如下所述:
步骤一、在单倾斜山体进行地面磁源瞬变电磁实验,采用中心方式,发射线圈依山体铺设成边长为200的正方形,杜瓦必须竖直放置于发射线圈中心点;
步骤二、发射机激发梯形波,采用低温超导磁传感器在测区内进行磁场测量,测量结束后,在发射线圈中心放置闭合的异常环,再次进行磁场测量。每一次测量时,利用采集器记录发射波形得到关断时间toff,发射电流为I;异常环的半径是0.3m,匝数是100匝;
步骤三、首先通过多次叠加消除噪声对测量数据进行初步的处理,然后用有异常环时的大地响应与没有异常环时的大地响应作差,得到异常环引起的磁场响应;
步骤四、根据异常环的几何参数,计算出异常环的理论上的中心磁场;
异常环可以等效为一个电阻和一个电感串联,如图3所示。当发射电流为梯形波时,得到等效回路中的感应电流表达式为:
t<0 i(t)=0;
式中,τ是异常环的时间常数,是等效回路切割一次场的磁通量,toff是发射电流的关断时间,L是异常环的等效电感;
式中,μ0是真空磁导率,n是异常环的匝数,r是异常环的半径。
给定异常环的参数,包括时间常数,电感,电阻,带入上述公式,即可计算出异常环理论上的磁场响应。
步骤五、得到校正系数cosθ=Banomaly/Btheory,该系数用于该测点的电磁数据校正,即把异常环的实测数据除以cosθ,得到与实际地形垂直方向上的中心磁场响应。利用校正之后的数据进行数据解释,可以提高解释精度。
图5是本发明一个实施例测区磁场传感器测量系统的数据地形校正效果图;从图中可以看出,校正后的曲线与理论曲线更加贴近,充分验证了基于异常环的磁场传感器测量系统的地形校正方法的有效性,为地面瞬变电磁法野外高精度测量提供了新的思路和方法。
Claims (1)
1.一种基于异常环的低温超导磁源瞬变电磁地形校正方法,其特征在于:其方法如下所述:
步骤一、在单倾斜山体进行地面磁源瞬变电磁实验,采用中心方式,发射线圈依山体铺设成边长为a的正方形,低温超导传感器保存在真空杜瓦盛装的液氦里,为确保实验安全,杜瓦必须竖直放置于发射线圈中心点,不随山体倾斜角度变化;
步骤二、发射机激发梯形波,利用采集器记录发射波形得到关断时间toff,发射电流为I;
步骤三、利用低温超导传感器在发射线圈中心点测量两次,分别是中心点有异常环和无异常环时的情况,测量的响应分别是Btotal和Bground,Btotal同时包含了大地的响应Bground和异常环的响应Banomaly,即Banomaly=Btotal-Bground;
步骤四、理论计算得到异常环的理论响应Btheory,异常环等效为一个电阻和一个电感串联,当发射电流为梯形波时,得到等效回路中的感应电流表达式如下:
t<0 i(t)=0;
0≤t<toff
t≥toff
上式中,τ是异常环的时间常数,是等效回路切割一次场的磁通量,toff是发射电流的关断时间,L是异常环的等效电感;Btheory结果如下式:
式中,μ0是真空磁导率,n是异常环的匝数,r是异常环的半径;
步骤五、得到校正系数cosθ=Banomaly/Btheory,该系数用于大地的电磁数据校正,即把异常环的实测数据除以cosθ,得到与实际地形垂直方向上的中心磁场响应。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710256983.9A CN107024722B (zh) | 2017-04-19 | 2017-04-19 | 一种基于异常环的低温超导磁源瞬变电磁地形校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710256983.9A CN107024722B (zh) | 2017-04-19 | 2017-04-19 | 一种基于异常环的低温超导磁源瞬变电磁地形校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107024722A CN107024722A (zh) | 2017-08-08 |
CN107024722B true CN107024722B (zh) | 2018-09-28 |
Family
ID=59528058
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710256983.9A Active CN107024722B (zh) | 2017-04-19 | 2017-04-19 | 一种基于异常环的低温超导磁源瞬变电磁地形校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107024722B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110133733B (zh) * | 2019-04-28 | 2020-07-24 | 吉林大学 | 一种基于粒子群优化算法的电导-极化率多参数成像方法 |
CN110068870B (zh) * | 2019-04-28 | 2020-08-28 | 中国科学院上海微系统与信息技术研究所 | 一种超导瞬变电磁信号的测量装置及方法 |
CN110501748B (zh) * | 2019-09-09 | 2021-06-01 | 中煤科工集团西安研究院有限公司 | 一种半航空瞬变电磁接收探头姿态矫正系统与方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102096113A (zh) * | 2010-12-03 | 2011-06-15 | 吉林大学 | 时间域地空电磁探测系统及标定方法 |
US7977940B2 (en) * | 2005-03-09 | 2011-07-12 | Andreas Chwala | Low temperature squid transient electromagnetic receiver system |
CN102495431A (zh) * | 2011-11-16 | 2012-06-13 | 中煤科工集团西安研究院 | 瞬变电磁数据静校正的数据处理方法 |
CN103064124A (zh) * | 2013-01-31 | 2013-04-24 | 江苏大学 | 一种校正电磁勘探中地形影响的比值方法 |
CN104698502A (zh) * | 2015-03-22 | 2015-06-10 | 山西煤炭进出口集团科技发展有限公司 | 由井出发的瞬变电磁勘探地形校正处理方法 |
CN104714254A (zh) * | 2015-02-02 | 2015-06-17 | 中国科学院地质与地球物理研究所 | 一种基于瞬变电磁法探测煤矿积水采空区的方法和装置 |
CN105676295A (zh) * | 2016-01-30 | 2016-06-15 | 吉林大学 | 基于squid的磁源激发极化-感应的联合探测系统与方法 |
CN105759316A (zh) * | 2016-02-04 | 2016-07-13 | 中国科学院地质与地球物理研究所 | 一种矩形回线源瞬变电磁探测的方法和装置 |
-
2017
- 2017-04-19 CN CN201710256983.9A patent/CN107024722B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7977940B2 (en) * | 2005-03-09 | 2011-07-12 | Andreas Chwala | Low temperature squid transient electromagnetic receiver system |
CN102096113A (zh) * | 2010-12-03 | 2011-06-15 | 吉林大学 | 时间域地空电磁探测系统及标定方法 |
CN102495431A (zh) * | 2011-11-16 | 2012-06-13 | 中煤科工集团西安研究院 | 瞬变电磁数据静校正的数据处理方法 |
CN103064124A (zh) * | 2013-01-31 | 2013-04-24 | 江苏大学 | 一种校正电磁勘探中地形影响的比值方法 |
CN104714254A (zh) * | 2015-02-02 | 2015-06-17 | 中国科学院地质与地球物理研究所 | 一种基于瞬变电磁法探测煤矿积水采空区的方法和装置 |
CN104698502A (zh) * | 2015-03-22 | 2015-06-10 | 山西煤炭进出口集团科技发展有限公司 | 由井出发的瞬变电磁勘探地形校正处理方法 |
CN105676295A (zh) * | 2016-01-30 | 2016-06-15 | 吉林大学 | 基于squid的磁源激发极化-感应的联合探测系统与方法 |
CN105759316A (zh) * | 2016-02-04 | 2016-07-13 | 中国科学院地质与地球物理研究所 | 一种矩形回线源瞬变电磁探测的方法和装置 |
Non-Patent Citations (3)
Title |
---|
Quantifying AEM system characteristics using a ground loop;Dacis A 等;《geophysics》;20081231;第73卷(第4期);第179-188页 * |
基于瞬变电磁法的铁峰煤矿采空区探测研究;牟铁超;《中国优秀硕士学位论文全文数据库 工程科技I辑》;20170315(第03(2017)期);第B021-20页 * |
接地电性源空地瞬变电磁探测传感器的研制;林君 等;《吉林大学学报(工学版)》;20140531;第44卷(第3期);第888-893页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107024722A (zh) | 2017-08-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102736114B (zh) | 电磁及其综合勘探的设备及方法 | |
CN105044792B (zh) | 地-井时频电磁勘探数据采集装置及方法 | |
CN106772644A (zh) | 矿井瞬变电磁三分量探测方法 | |
CN103064120B (zh) | 煤矿井下磁电综合探测方法 | |
CN107024722B (zh) | 一种基于异常环的低温超导磁源瞬变电磁地形校正方法 | |
Zhen‐Zhu et al. | Opposing coils transient electromagnetic method for shallow subsurface detection | |
CN108008451A (zh) | 一种瞬变电磁探测装置及消除电感干扰的方法 | |
CN105353428B (zh) | 一种地面参考区磁场延拓的地空协同电磁数据校正方法 | |
CN104035137B (zh) | 地下全空间瞬变电磁探测仪及探测方法 | |
CN103837900A (zh) | 一种基于矢量磁场探测的地下电缆定位方法及装置 | |
CN104407392A (zh) | 一发三收式对充水采空区的探测装置及探测方法 | |
CN104360401A (zh) | 一种瞬变电磁b场确定地下目标体地质信息方法 | |
CN104612671A (zh) | 一种测量垂直方向地层电导率的阵列感应线圈系 | |
CN107191181A (zh) | 一种基于电磁散射的井周界面探测方法 | |
CN105549100B (zh) | 基于u形螺线源的矿井瞬变电磁探测装置与方法 | |
CN109490968A (zh) | 一种可控源高频大地电磁仪器系统及测量方法 | |
CN103499844B (zh) | 全程全空间瞬变电磁装置及测量方法 | |
CN108957560A (zh) | 一种适用于城镇地下空间探测的电法勘探方法 | |
CN113608267A (zh) | 一种双回线零磁通瞬变电磁探测装置 | |
CN108287367A (zh) | 一种基于时域电磁法的地面-孔中联合探测系统及应用方法 | |
CN112068211A (zh) | 一种无人机半航空时间域电磁勘探系统 | |
US3321700A (en) | Method and apparatus for geophysical exploration by measuring time variations in theearth's magnetic field at a plurality of locations | |
CN204855829U (zh) | 地-井时频电磁勘探数据采集装置 | |
CN114217354A (zh) | 基于光纤电磁传感器的电磁数据采集系统及采集方法 | |
CN103499838A (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 |