台站地磁观测数据精度提高方法.pdfVIP

地 震 地 磁 观 测 与 研 究 第 35卷 第 5／6期 SEISM 0L0GICAL AND GE0MAGNETIC Vo1．35 NO．5／6 2014年 12月 0BSERVAT10N AND RESEARCH Dec． 2O14 doi：10．3969／j．issn．1003—3246．2014．os／o6．042 台站地磁观测数据精度提高方法 徐 涛 吴宝元 宁百齐 (中国北京 100029中国科学 院地质与地球物理研究所) 摘要 在当前多种地磁 日常数据处理方法基础上 ，探讨提高地磁观测数据精度 的问题 。通过使用 一 种改进 的地磁观测基线值计算方法 ，能有效减小外磁场变化对观测结果的影响，考虑到基线值 日变化等因素，使用高阶多项式拟合基线值 日变化 曲线 。利用该方法处理的地磁数据与实测数据 进行结果 比对后 ，发现此方法能有效提高地磁观测数据 的精度和质量 。 关键词 地磁观测 ；地磁基线值 ；曲线拟合 ；相关性分析 ；正态分布检验 0 引言 目前地磁观测实现数字化和 自动化 ，能组建更高效 的观测系统 ，提供更多观测手段 ，对数 据质量随之提出更高要求 ，这就需要数据处理方法与时俱进 。在数据预处理阶段 ，各要素基线 值的计算及其Et变化 曲线绘制有许多技巧和注意事项值得讨论。以往通过 DI仪得到一组倾 角 值后 ，使用 G856手动采集总场 F值 ，从而得到水平分量 H 和垂直分量 z的观测值 。只 是 I和F(H、Z)的观测时间不在 同一时段 ，如果此时外界磁场变化剧烈 ，H 分量和Z分量的 观测值必然会存在误差 (陆镜辉 ，2008)。同样的，相对记录仪器得到的D的变化值是观测时 间段的均值 ，与读数时刻的真实变化值有一定误差 。这些误差可以通过改进基线值算法来规 避 。本文还利用数理统计方法，比对中国科学院地质与地球物理研究所空间探测实验室的北 京十三陵地磁台多套地磁仪器观测结果 ，分析误差来源 。 1 北京十三陵地磁 台的地磁观测系统 北京十三陵地磁 台坐落在北京十三陵风景 区，四周为 国有 山林地 ，环境开 阔，磁场分布 均匀 ，无磁污染源。由3套 D、H、Z相对变化量记录仪器 、两套 F标量记录仪器、两套 D、绝 对观测仪器以及一套记录D、H、Z增量变化 的仪器等组成若干套数字化地磁观测系统 ，见表 1。主要观测系统由MAGO1H和 Overhauser对观察室标准墩 的D、、F三要素进行绝对观 测 ，以此对相对记录仪 FGE、LEM1018和 GM4基线值进行标度 ，构建联合观测体系。绝对观 测系统均采用两套仪器 同时进行观测和记录，通过观测数据的直接 比对 ，保证数据的完整性和 可靠性 。相对记录仪器和Overhauser的探头安置在半地下室 ，室内 日温差小于0．5℃，年温差 约 12℃，相对湿度小于 70 。 作者简介 ：徐涛(1982一)，男，湖北武汉人，工程师，主要从事地磁观测与资料分析工作 本文收到 日期 ：2014—07—15 238 地 震 地 磁 观 测 与 研 究 35卷 2 基线值处理方法 地磁相对记录仪器连续记录地磁三要素的变化值 ，需要绝对观测仪器定期对其基线值进 行定标 。基线值计算公式为 PBa5 一Pob 一 △P (P表示D、、F、H、Z等要素) (1) 式中，P。 为绝对观测值 ，每周 2次 ，每次 2组观测值 ；△P为相对记录仪器采集的连续且高 分辨率的地磁三要素相对变化值 ；P 为地磁观测基线值 ，采集频率为一周 2次 ，共 4组 数据 。 观测基线值是一组离散 的大小不一的点的集合 。通常的处理方法是 ，将这些离散 的点，按 照 日期排序后 ，描绘在绘 图纸上 ，其 中横坐标是 以 “天”为单位 ，纵坐标是各要素的基线值 。 用曲尺描绘出一条平滑的曲线 ，使得这些采样点尽可能分布在曲线周 围。这条 曲线即为基线 值的 日变化曲线 。即可将每周两天的观测基线值，拟合全年每天