地面瞬变电磁法数据处理和成像系统TEM_SW

0.1
0.05
0
-0.05 0 20 40 Sampling time (ms) 60
划分出工频干扰
0.2 1 Norm alized induc ctive EM (uV/A) F 2 0.15 3 0.1
0.05
0
-0.05 0 20 40 Sampling time (ms) 60
相继压制随机干扰和工频干扰结果
理论数据成像结果
理 理 理 理 TEM数 数 数 数 数 数
全 全 全 全 全 全 全 全 全 全 全
Profile / m
全 拟 拟 拟 拟 全 全 全 全 全 全
Profile / m
100 200 300 400 500 600 700 800 900 1000 100 200 300 400 500 600 700 800 900 1000 10 欧 欧 欧 10 欧 欧 欧 100 100 0 -100 -200 Elevation / m -300 -400 -500 -600 -700 -800 -900 -1000
输入文件的例子DEMO1.DAT 输入文件的例子
中心回线装置理论数据 20 1 6000 10.00 1 测点数，供电波形，采样时间道数，回线边长，装置类型） （测点数，供电波形，采样时间道数，回线边长，装置类型） 100.0000 10000.00 100.0000 （发送框面积，接收框面积，收发距） 发送框面积，接收框面积，收发距） 10.00000 0.0000000E+00 供电脉宽，断电后沿宽度） （供电脉宽，断电后沿宽度） DH X Y Z 1 50.00 0.00 78.00 2 100.00 0.00 92.30 3 150.00 0.00 102.00 …… 19 950.00 0.00 75.00 20 1000.00 0.00 70.00 0.1000E-01 0.2000E-01 0.3000E-01 0.4000E-01 0.5000E-01 0.6000E-01 0.7000E-01 0.9000E-01 0.1000E+00 0.1100E+00 0.1200E+00 0.1300E+00 0.1400E+00 0.1500E+00 0.1700E+00 0.1800E+00 0.1900E+00 0.2000E+00 0.2100E+00 0.2200E+00 0.2300E+00 …… …… （以上为采样时间 ms） ） …… 0.1485E+07 0.2764E+06 0.1016E+06 0.4911E+05 0.2762E+05 0.1714E+05 0.1139E+05 0.5806E+04 0.4364E+04 0.3368E+04 0.2657E+04 0.2134E+04 0.1741E+04 0.1441E+04 0.1021E+04 0.8726E+03 0.7518E+03 0.6530E+03 0.5711E+03 0.5025E+03 0.4447E+03 …… …… 以上为归一化感应电动势µ （以上为归一化感应电动势µV/A） ）
具体功能
• 剔除单点突变的过失误差校正 剔除单点突变的过失误差校正。 • 压制随机起伏的偶然误差校正和工频干扰 校正。 校正。 • 筛选合理采样时间的实测数据：供绘制多 筛选合理采样时间的实测数据： 测道剖面图。 • 计算全时域视电阻率：实现视电阻率拟断 计算全时域视电阻率： 面图成像。 • 瞬变电磁反演： 瞬变电磁反演： 1. 拟二维反演 拟二维反演：实现本征电阻率断面图成像 2. 逐点一维反演： 逐点一维反演：
10 欧 欧 欧
紫 紫 紫 紫 紫 紫 紫 理 理 理 理
100 200 300 400 500 600 700 800 900 1000 ρ ( m) 55 50 45 40 35 30 25 20 15 10 5
100 200 300 400 500 600 700 800 900 1000
很愿意和同行们交流！ 很愿意和同行们交流！ 谢谢！ 谢谢！
六、成像结果举例 压制随机干扰和工频干扰 的理论算例
加有干扰的理论瞬变响应曲线
0.2 1 Normalized induct ctive EMF (uV/A) 0.15
0.1
0.05
0
-0.05 0 20 40 Sampling time (ms) 60
多项式拟和压制随机干扰后
0.2 1 N orm alized induc ctive EM (uV/A) F 2 0.15
0
-0.05 0 20 40 Sampling time (ms) 60
压制随机干扰的实例
（山西某煤矿CUGTEM） 山西某煤矿 ）
N rm liz d in u tiv E F(u /A o a e dc e M V ) -1 0 2 3 1 0 4 8 12 Sampling time (ms) 1 16 4
三、输入文件
数据处理和反演所需的实际观测条件参数和实 测归一化感应电动势数据，采用数据文件输入。 测归一化感应电动势数据，采用数据文件输入。 输入数据文件的内容和格式： 输入数据文件的内容和格式： LINE （说明测区、测线的字符串） 说明测区、测线的字符串） ND, N, NTsc, AA, NARRAY TXAREA, RXAREA, TRDIST TD, TRMP 'DH X Y Z ' （测深点的测点号和坐标） 测深点的测点号和坐标） DH(i),X(i),Y(i),Z(i), i=1, ND (Tsc(j),j=1,NTsc) (Vsc(i,j),j=1,NTsc), i=1,ND
适用范围
• 适用于中心回线、重叠回线和偶极－ 适用于中心回线、重叠回线和偶极－ 中心回线 偶极三种装置类型 三种装置类型。 偶极三种装置类型。 多种波形， • 适用于多种波形，特别是适用于有一 适用于多种波形 断电后沿的供电波形 的供电波形。 定断电后沿的供电波形。 • 一次能处理的测点数上限为 一次能处理的测点数上限为500个； 个 • 每个测点的最大采样时间数为 每个测点的最大采样时间数为6000； ； • 筛选出的合理采样时间数不大于 ； 筛选出的合理采样时间数不大于60； • 逐点一维反演的最多层数为 层。 逐点一维反演的最多层数为10层
0.8000E 0.1600E+00 0.2400E+00
0.7972E+04 0.1207E+04 0.3956E+03
……
四、程序的运行 演示） （演示）
五、输出文件（共五个） 输出文件（共五个）
• 第一个输出文件“Vsx*** +输入文件名” ，是供 第一个输出文件“ 输入文件名” 输入文件名 绘制多测道剖面图的数据文件。 绘制多测道剖面图的数据文件。 • 第二个输出文件“Rs ***+输入文件名” ，是全 第二个输出文件“ 输入文件名” 输入文件名 时域视电阻率的计算结果。 时域视电阻率的计算结果。由它可借助专业绘图 软件Surfer，绘制视电阻率拟断面图。 软件 ，绘制视电阻率拟断面图。 • 第三个输出文件“R***+输入文件名” ，是拟二 第三个输出文件“ 输入文件名” 输入文件名 维反演结果。由它可借助专业绘图软件Surfer， 维反演结果。由它可借助专业绘图软件 ， 绘制反演所得本征电阻率断面图。 绘制反演所得本征电阻率断面图。 • 第四个输出文件“测深点号 输入文件名”，是 第四个输出文件“测深点号+输入文件名 输入文件名” 逐点一维反演的结果。 逐点一维反演的结果。 • 第五个输出文件“bln+输入文件名” ，输出观 第五个输出文件“ 输入文件名” 输入文件名 测剖面上各测深点的编号、横向坐标和地形高程。 测剖面上各测深点的编号、横向坐标和地形高程。
紫 紫 紫 紫 紫 紫 紫 理 理 理 理
20 欧 欧 欧
20 欧 欧 欧
0 20 欧 -100 -200
欧 欧
20 欧 欧 欧
50 欧 欧 欧
50 欧 欧 欧
-300 -400 -500
50 欧 欧 欧
50 欧 欧 欧
10 欧 欧 欧
Байду номын сангаас
10 欧 欧 欧 -600 -700 -800 -900
10 欧 欧 欧
实测瞬变响应曲线
N rm lize in u tiv E F (u /A o a d dc e M V ) 0 1 -1 2 3 0 4 8 12 Sampling time (ms) 1 2 16 4
压制随机干扰后
N rm liz d in u tiv E F (u /A o a e dc e M V ) 2 3 4 0 1 -1 0 4 8 12 Sampling time (ms) 1 2 3 16
二、程序运行的硬件和软件环境 • 可以在任何个人计算机（PC机） 可以在任何个人计算机（ 机 上运行。 上运行。 • 运行环境要求有 硬盘富余空间 运行环境要求有1G硬盘富余空间 和内存不低于1G。 和内存不低于 。 • 系统以 系统以WINDOWS为支撑平台， 为支撑平台， 为支撑平台 要求WINDOWS 98以上的版本。 以上的版本。 要求 以上的版本
压制随机干扰后
1E-007 N orm alized indu uctive E F (uV/A M ) 1 8E-008 2 6E-008
4E-008
2E-008
0
-2E-008 0 10 20 Sampling time (ms) 30
划分出工频干扰
1E-007 N orm alized indu uctive EMF (uV/A ) 1 8E-008 2 6E-008 3
0.2 1 Normalized induc ctive EMF (uV/A) 2 0.15 3 0.1 4
0.05
0
-0.05 0 20 40 Sampling time (ms) 60
无干扰信号和压制干扰结果几乎重合 无干扰信号和压制干扰结果几乎重合
0.2 1 Normalized induct ctive EMF (uV/A) 2 0.15 3 0.1 4 5 0.05
筛选采样时间后
压制随机干扰和 工频干扰的实例
广西南丹某锡矿GDP-32 广西南丹某锡矿
实测瞬变响应曲线
1E-007 N orm alized indu uctive E F (uV ) M /A 1 8E-008
6E-008
4E-008
2E-008
0
-2E-008 0 10 20 Sampling time (ms) 30
4E-008
2E-008
0
-2E-008 0 10 20 Sampling time (ms) 30
工频干扰校正后的结果
1E-007 N orm alized indu uctive E F (uV ) M /A 1 8E-008 2 6E-008 3 4 4E-008
2E-008
0
-2E-008 0 10 20 Sampling time (ms) 30
地面瞬变电磁法数据处理和成像 系统TEM_SW 系统
中国地质大学 电法科研组 RGEP,CUG 罗延钟
一、系统的功能 二、程序运行的硬件和软件环境 三、输入文件 四、程序的运行 五、输出文件 六、成像结果举例
一、系统的功能
用于对一条观测剖面的 用于对一条观测剖面的 一条观测剖面 实测感应电动势数据作 数据处理和二维成像。 数据处理和二维成像。