g856和gsm-19质子磁力仪数据采集程序的设计与实现.docxVIP

g856和gsm-19质子磁力仪数据采集程序的设计与实现 0 配套日变道德 在同一测区，使用g856或gm-19t两种类型的质子磁力仪来采集磁力仪数据。两种磁力仪的测量数据的输出格式不同，并且两种磁力仪数据的日变校正软件不同。G856质子磁力仪数据输出文件格式为stn, 日变改正软件为Mag Map2000;GSM-19T质子磁力仪数据输出文件格式为txt, 日变改正软件为GEMLink W。因此, 这两种型号的磁测观测数据用配套软件相互日变改正特别困难和繁琐。此外, 根据《地面高精度磁测技术规程》 (DZ/T0071-93) , 为了提高日变改正精度, 先要对日变观测原始数据作五点或七点滑动平均, 压制噪声后, 再对测点观测数据进行日变改正, 但G856和GSM-19T质子磁力仪配套日变改正软件日变改正过程中没有五点或七点滑动平均。针对这些问题, 查阅了相关的文献和软件, 管志宁提出日变改正值从日变曲线上查得, 在日变曲线上量得某时刻相对早基点时间的日变值并取反号, 即为该时刻的日变改正值, 这种日变改正方法速度慢、精度不高。于宝显等提出不同型号质子磁力仪之间数据共享, 经转换后的数据可以利用其中一种型号质子磁力仪所配套软件做日变改正, 这种方法有些繁琐, 并且没有做五点或七点滑动平均。中国地质大学 (武汉) 刘天佑教授开发的磁法勘探软件系统MAGS 2.0没有实现这两种型号质子磁力仪观测数据之间的日变改正处理。对此, 通过Microsoft Office Excel平台VBA编程解决。 1 dat格式文件数据输出 为了统一G856和GSM-19T质子磁力仪观测数据格式, 用Mag Map2000导入G856质子磁力仪stn格式文件数据, 然后输出dat格式文件数据, 最后把dat文件格式数据复制到Microsoft Office Excel文件中;GSM-19T质子磁力仪txt文件格式数据可以直接导入Microsoft Office Excel文件中。 2 日变难以消除的误差 根据《地面高精度磁测技术规程》 (DZ/T0071-93) , 地磁场短周期变化的振幅与质子磁力仪噪声均方根值是近于同一数量级, 噪声是随机的, 地磁场脉动变化是有规律的, 先要对日变观测原始数据先作五点或七点滑动平均, 压制噪声水平后, 再对测点观测数据进行日变改正, 即可提高日变改正的精度。下面介绍日变观测数据五点和七点滑动平均原理。设在一个范围内某点经平滑后的日变数值可表示: (1) 式中:a0和a1为待定系数。设日变观测的数值为T (xi) , 它的平滑值为, 则可得出: (2) 式中δ为偏差的平方和。利用微分求极值的方法将 (2) 式对a0和a1求导数, 然后令其为零得: 若xi以剖面上的点距为单位, 步长为1, 则式 (3) 中的xi=0, ±1, ±2…±m。代入式 (3) 中, 可解得a0和a1分别为 由此可见, (5) 式就是平滑公式, 当m=±1, 三点平滑公式为 同理可得五点平滑公式为 七点平滑公式为 根据五点和七点滑动平均公式 (7) 和 (8) 编写Microsoft Office Excel VBA程序, 代码如下, 利用此程序代码把日变原始观测数据 (见图1) 五点和七点滑动平均, 见图2和图3。 3 日变观测分段线性插值 根据《地面高精度磁测技术规程》 (DZ/T0071-93) 质子磁力仪日变观测时间间隔为5秒至20秒, 为了求出测点观测与日变观测同一时刻的日变数据, 先要求出日变观测数据每一秒的数值, 因此对日变数据进行分段线性插值, 公式为 (9) 式, 公式中T (n) 为前一时刻的日变观测数据, T (n+m) 为后一时刻的日变观测数据, m为日变观测时间间隔, j为日变观测第j秒, T (j) 为分段线性插值后第j秒数值, Microsoft Office Excel VBA程序代码如下: 4 观测数据的处理 求出日变观测每一秒数值后, 要提取测点观测与日变观测同一时刻的日变数据, Microsoft Office Excel VBA程序如下: 5 日的变化是固定 提取测点与日变观测同一时刻的日变数据后, 用测点观测的数据减去同一时刻的日变数值就是日变改正后的值。 6 日变难以消除的情况 为了验证VBA程序的正确性, 在野外磁测过程中做了实验验证。用一台GSM-19T质子磁力仪做日变观测, 两台GSM-19T质子磁力仪做测点观测, 点距约为2米, 点数为64个。测点观测数据用GEMLink W软件和VBA程序进行日变改正, 用VBA程序日变改正时没有做五点和七点滑动平均。日变改正后的数据见表1。从表1可以看出:GEMLink W软件和VBA程序日变改正后的值相等 (用GEMLink W软件日变改正后的值软