CN108303753A - 地磁矢量场观测数据一致性校正方法 - Google Patents
地磁矢量场观测数据一致性校正方法 Download PDFInfo
- Publication number
- CN108303753A CN108303753A CN201810101718.8A CN201810101718A CN108303753A CN 108303753 A CN108303753 A CN 108303753A CN 201810101718 A CN201810101718 A CN 201810101718A CN 108303753 A CN108303753 A CN 108303753A
- Authority
- CN
- China
- Prior art keywords
- data
- components
- curves
- observation
- component
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V13/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
本发明公开了一种地磁矢量场观测数据一致性校正方法,通过高斯滤波法将高采样率的观测设备数据转换为低采样率数据,接着,通过简单归一化处理使各个分量数据的起始点均为0值,通过被校正物理量的差值曲线与地磁矢量场各物理量的日变化曲线的相关性分析,确定出主要误差因素,再通过差值曲线的峰峰值与相关系数最高的物理量日变化曲线的峰峰值的比例系数,确定出主要误差修正系数并进行第一次数据误差修正,再接着,对修正后的数据再进行第二次数据误差修正,最后,完成各分量地磁矢量场观测数据一致性校正。该方法可大大提高观测数据的准确性,并且校正过程易懂,计算公式简单,适合在科研观测数据校正过程中使用。
Description
技术领域
本发明涉及一种观测数据修正方法,特别涉及地磁矢量场观测数据一致性的校正方法,属于地球物理研究领域。
背景技术
地磁矢量场测量一般包括地磁南北水平强度(水平分量)H,向北为正,地磁东西水平强度(水平分量)D,向东为正,垂直强度(垂直分量)Z,向下为正。目前,中国地震局地磁台网经过“九五”、“十五”、“十一五”的升级改造,地磁相对记录已经实现数字化、网络化,同时,实现了秒数据或分数据的采集、存储。
地磁台站的相对记录一般为两套同类型仪器并行观测,以便及时鉴别观测设备故障和环境干扰,并确保观测数据的连续性和完整性。经过近几年的观测,发现有些台站同类型仪器的H、D、Z日变观测曲线不是完全重合,即存在不一致性问题,研究人员认为导致数据不一致主要有三个误差源:1、观测仪器本身性能导致误差产生,由于受加工和装配工艺水平的限制,难以保证仪器的正交度,在仪器调试过程中,格值、温度系数、零点漂移等情况引入产生误差;2、仪器安装过程引入误差,其包括仪器D分量定向、底座水平调节难以保证测量坐标系XOY坐标平面的绝对水平;3、观测环境变化产生误差,由于季节变化等因素引起的观测墩漂移以及台站供电***不稳定等因素产生的误差。其中,观测仪器的姿态和格值误差为主要的误差因素,也是本发明想要校正的主要因素。
发明内容
鉴于上述情况,本发明旨在提供一种地磁矢量场观测数据一致性校正方法,以避免数据误差,满足提供高质量地磁矢量场观测数据的需要,为科学研究提供丰富、准确的观测数据。
本发明通过差值法寻找地磁日变化观测数据曲线不一致的误差因素,并给出具体的校正比例系数分别进行H、D、Z分量的日变观测曲线误差修正,以实现对应观测数据的一致性和准确性。其中,差值法是利用两套观测设备的同一观测物理量之间的数据差值,并根据差值与标准观测设备中各观测分量的比例关系进行计算的一种简捷而快速的修正方法,利用差值校正的关键在于寻求差量与标准地磁日变化各观测分量之间的比例关系,以差量作为修正的突破口。
具体方法是:首先,因为各观测设备的采样率可能存在不同,即存在分采用率和秒采样率,通过高斯滤波法将高采样率的观测设备数据转换为低采样率;其次,为了直观表示观测误差形态,对转换后的数据做简单归一化处理后使得各分量的起始点均为0值;然后,通过被校正物理量的差值曲线与标准设备的地磁矢量场各物理分量的日变化观测曲线进行相关性分析,确定主要的误差因素,若被校正物理量的差值曲线与标准设备本身物理量的日变化观测曲线相关系数较,则判断存在格值误差,否则,判断存在姿态误差;接着,通过差值曲线的峰峰值与相关系数最高的物理量日变化观测曲线的峰峰值的比例系数确定出第一次修正系数,进行第一次数据修正,然后对校正后的数据重复上述过程进行第二次数据修正,完成该物理量的最终数据修正;最后,依次对其它物理量进行上述过程修正,完成所有物理量的地磁矢量场观测数据一致性校正。
本发明是通过以下技术来实现的:
一种地磁矢量场观测数据一致性校正方法,具体步骤包括:
步骤1、如两个以上同类型地磁矢量场观测设备采样率不一致,应先对高采样率观测设备的H分量、D分量和Z分量的日变观测数据曲线进行对应分量的低采样率数据转换,即采用高斯滤波法将秒采样率观测数据转换为分采样率观测数据,计算第i分钟(im)的分钟值时,取im00s及其前后各M秒的秒采样数据进行高斯滤波计算,公式如下:
其中,是第i分钟的分钟值,是第i分钟前后取秒采样数据个数,是第i前一分钟第60+n秒的秒采样数据,是第i分钟00秒的秒采样数据,是第i分钟第n秒的秒采样数据,和为高斯系数。
如两个以上同类型地磁矢量场观测设备的采样率一致,则直接进行步骤2。
步骤2、对转换后数据或采样率一致的H分量、D分量和Z分量的日变观测数据曲线分别做简单归一化处理,使H分量、D分量和Z分量的日变观测数据曲线的起始点均为0值,同时,设定同类型地磁矢量场标准观测设备的各日变观测分量数据曲线为、、。
步骤3、对归一化后的两个以上D分量的日变观测数据曲线进行差值化处理,得到被校正物理量的差值数据曲线。
步骤4、将差值数据曲线与标准观测设备的数据曲线、、分别调用Spearman函数进行相关性分析;若差值数据曲线与或数据曲线相似,则判断存在姿态误差;若差值数据曲线与数据曲线相似,则判断存在格值误差。
步骤5、若存在姿态误差,且与数据曲线相似,第一次数据修正系数为:
。
D分量第一次修正后的数据为:
。
若存在姿态误差,且与数据曲线相似,第一次数据修正系数为:
。
D分量第一次修正后的数据为:
。
若存在格值误差,则第一次数据修正系数为:
。
D分量第一次修正后的数据为:
。
其中,、、分别是标准观测设备H、Z、D分量日变观测数据曲线的峰峰值,是D分量的差值数据曲线的峰峰值。
步骤6、将修正后的数据曲线与标准观测设备的数据曲线做差值化处理,得到被校正物理量的修正后的差值数据曲线,并重复上述步骤4、5,进行第二次D分量数据一致性修正,得到最终修正后的D分量数据。
步骤7、重复上述步骤3、4、5、6,分别对H分量和Z分量进行数据一致性修正,完成最终修正后的H分量和Z分量数据。
本发明所述的地磁矢量场观测数据一致性校正方法的有益效果包括:
1、利用差值法简捷、快速的判断地磁日变化观测数据曲线不一致的误差因素,并给出具体的校正比例系数,通过两次误差修正得到最终校正后的数据,大大提高了数据的准确性。
2、数据校正方法易懂、计算公式简单,适合在科研观测数据修正中使用。
附图说明
图1为本发明所述地磁矢量场观测数据一致性校正方法的数据流程图。
具体实施方式
本发明的中心思想是:将地磁台站进行矢量场观测的多套设备的一个物理量通过差值法计算校正比例系数,修正到一套标准观测设备的同一物理量上,完成该物理量的观测数据修正,并经上述步骤的重复使用,最终完成多个物理量的数据修正。
下面以同一观测站上两个同类型地磁矢量场观测设备的观测数据为例,其采样率分别为秒采样和分采用,结合附图1对本发明所述方法做进一步的详细描述:
本发明所述的地磁矢量场观测数据一致性校正方法,具体步骤包括:
步骤1、对两套同类型地磁矢量场观测设备的H分量、D分量和Z分量的日变观测数据曲线进行对应分量数据的转换,即采用高斯滤波法将秒采样率观测数据转换为分采样率观测数据,计算第i分钟(im)的分钟值时,取im00s及其前后各45秒共91秒的秒采样数据进行高斯滤波计算,公式如下:
其中,是第i分钟的分钟值,为采样数据个数45,是第i前一分钟第60+n秒的秒采样数据,是第i分钟00秒的秒采样数据,是第i分钟第n秒的秒采样数据,和为高斯系数。
转换后数据单位一致,从而可以形成相互对应的分量数据,为后续差值化处理提供方便。当然,如果两个同类型地磁矢量场观测设备的采样率一致,则可直接进行步骤2。
步骤2、对转换后数据的H分量、D分量和Z分量的日变观测数据曲线分别做简单归一化处理,使H分量、D分量和Z分量的日变观测数据曲线的起始点均为0值,同时,设定同类型地磁矢量场标准观测设备的各日变观测分量数据曲线为、、。
通常情况下,通过地磁台站的基线值来确定标准观测设备,进而得到标准的H分量、D分量和Z分量的数据曲线、、。
步骤3、对归一化后的两个D分量的日变观测数据曲线进行差值化处理,得到被校正物理量的差值数据曲线。
本步骤首先选择D分量进行数据一致性校正,当然也可先选择H分量或Z分量进行校正。
步骤4、将差值数据曲线与标准观测设备的数据曲线、、分别调用Spearman函数进行相关性分析;如果差值数据曲线与或数据曲线相似,则判断存在姿态误差;如果差值数据曲线与数据曲线相似,则判断存在格值误差。
因为斯皮尔曼Spearman等级相关对数据条件的要求没有积差相关系数严格,只要两个变量的观测值是成对的等级评定资料,其特别适用于两个变量的总体分布形态比较,进而完成数据相似性判断。
步骤5、若存在姿态误差,且与数据曲线相似,第一次数据修正系数为:
D分量第一次修正后的数据为:
;
若存在姿态误差,且与数据曲线相似,第一次数据修正系数为:
D分量第一次修正后的数据为:
;
若存在格值误差,则第一次数据修正系数为:
,
D分量第一次修正后的数据为:
;
其中,、、分别是标准观测设备H、Z、D分量日变观测数据曲线的峰峰值,是D分量的差值数据曲线的峰峰值。
上述过程得到的修正系数为导致观测数据不一致的主要误差因素,也是修正数据的关键。
步骤6、将修正后的数据曲线与标准观测设备的数据曲线做差值化处理,得到被校正物理量的修正后的差值数据曲线,并将作为差值数据曲线重复上述步骤4、5,进行第二次D分量数据一致性修正,得到最终修正后的D分量数据。
为进一步提高修正效果和修正的准确性,本方法采用两次修正的方式,对第一次修正完成的数据再进行第二次修正,即通过得出第二次修正的次要误差因素,实现对观测数据的最终修正。
步骤7、重复上述步骤3、4、5、6,依次对H分量和Z分量进行数据一致性修正,完成最终修正后的H分量和Z分量数据。
通过对D分量、H分量和Z分量的依次修正,可准确实现同一观测站上两个同类型地磁矢量场观测设备各分量日变观测数据曲线的一致性修正,为科学研究提供丰富、准确的观测数据。
Claims (1)
1.一种地磁矢量场观测数据一致性校正方法,其特征在于,具体步骤包括:
步骤1、如两个以上同类型地磁矢量场观测设备采样率不一致,应先对高采样率观测设备的H分量、D分量和Z分量的日变观测数据曲线进行对应分量的低采样率数据转换,即采用高斯滤波法将秒采样率观测数据转换为分采样率观测数据,计算第i分钟(im)的分钟值时,取im00s及其前后各M秒的秒采样数据进行高斯滤波计算,公式如下:
其中,是第i分钟的分钟值,是第i分钟前后取秒采样数据个数,是第i前一分钟第60+n秒的秒采样数据,是第i分钟00秒的秒采样数据,是第i分钟第n秒的秒采样数据,和为高斯系数;
如两个以上同类型地磁矢量场观测设备的采样率一致,则直接进行步骤2;
步骤2、对转换后数据或采样率一致的H分量、D分量和Z分量的日变观测数据曲线分别做简单归一化处理,使H分量、D分量和Z分量的日变观测数据曲线的起始点均为0值,同时,设定同类型地磁矢量场标准观测设备的各日变观测分量数据曲线为、、;
步骤3、对归一化后的两个以上D分量的日变观测数据曲线进行差值化处理,得到被校正物理量的差值数据曲线;
步骤4、将差值数据曲线与标准观测设备的数据曲线、、分别调用Spearman函数进行相关性分析;
若差值数据曲线与或数据曲线相似,则判断存在姿态误差;
若差值数据曲线与数据曲线相似,则判断存在格值误差;
步骤5、若存在姿态误差,且与数据曲线相似,第一次数据修正系数为:
,
D分量第一次修正后的数据为:
;
若存在姿态误差,且与数据曲线相似,第一次数据修正系数为:
,
D分量第一次修正后的数据为:
;
若存在格值误差,则第一次数据修正系数为:
,
D分量第一次修正后的数据为:
;
其中,、、分别是标准观测设备H、Z、D分量日变观测数据曲线的峰峰值,是D分量的差值数据曲线的峰峰值;
步骤6、将修正后的数据曲线与标准观测设备的数据曲线做差值化处理,得到被校正物理量的修正后的差值数据曲线,并重复上述步骤4、5,进行第二次D分量数据一致性修正,得到最终修正后的D分量数据;
步骤7、重复上述步骤3、4、5、6,分别对H分量和Z分量进行数据一致性修正,完成最终修正后的H分量和Z分量数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810101718.8A CN108303753B (zh) | 2018-02-01 | 2018-02-01 | 地磁矢量场观测数据一致性校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810101718.8A CN108303753B (zh) | 2018-02-01 | 2018-02-01 | 地磁矢量场观测数据一致性校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108303753A true CN108303753A (zh) | 2018-07-20 |
CN108303753B CN108303753B (zh) | 2019-07-26 |
Family
ID=62851099
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810101718.8A Expired - Fee Related CN108303753B (zh) | 2018-02-01 | 2018-02-01 | 地磁矢量场观测数据一致性校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108303753B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109061753A (zh) * | 2018-10-26 | 2018-12-21 | 中国人民解放军61540部队 | 一种纬度和经度方向双因子定权的地磁数据通化方法 |
CN109188543A (zh) * | 2018-08-15 | 2019-01-11 | 中国人民解放军61540部队 | 一种地磁台站测量数据通化处理的双因子定权方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003166825A (ja) * | 2001-12-04 | 2003-06-13 | Nec Tokin Corp | 磁気センサ装置用地球磁極方位計算方法 |
CN103389517A (zh) * | 2013-07-30 | 2013-11-13 | 高建东 | 一种高精度地磁矢量测量方法及其装置 |
CN105722211A (zh) * | 2016-03-22 | 2016-06-29 | 中国科学院电子学研究所 | 利用地磁相关性的磁场传感器数据时间同步方法 |
CN106033131A (zh) * | 2015-03-20 | 2016-10-19 | 阿里巴巴集团控股有限公司 | 一种地磁传感器校准方法、装置及智能设备 |
CN106772683A (zh) * | 2017-01-12 | 2017-05-31 | 中国地震局地球物理研究所 | 一种简易测量矢量磁力仪中分量正交线圈正交度的方法 |
-
2018
- 2018-02-01 CN CN201810101718.8A patent/CN108303753B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003166825A (ja) * | 2001-12-04 | 2003-06-13 | Nec Tokin Corp | 磁気センサ装置用地球磁極方位計算方法 |
CN103389517A (zh) * | 2013-07-30 | 2013-11-13 | 高建东 | 一种高精度地磁矢量测量方法及其装置 |
CN106033131A (zh) * | 2015-03-20 | 2016-10-19 | 阿里巴巴集团控股有限公司 | 一种地磁传感器校准方法、装置及智能设备 |
CN105722211A (zh) * | 2016-03-22 | 2016-06-29 | 中国科学院电子学研究所 | 利用地磁相关性的磁场传感器数据时间同步方法 |
CN106772683A (zh) * | 2017-01-12 | 2017-05-31 | 中国地震局地球物理研究所 | 一种简易测量矢量磁力仪中分量正交线圈正交度的方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109188543A (zh) * | 2018-08-15 | 2019-01-11 | 中国人民解放军61540部队 | 一种地磁台站测量数据通化处理的双因子定权方法 |
CN109061753A (zh) * | 2018-10-26 | 2018-12-21 | 中国人民解放军61540部队 | 一种纬度和经度方向双因子定权的地磁数据通化方法 |
CN109061753B (zh) * | 2018-10-26 | 2020-04-28 | 中国人民解放军61540部队 | 一种纬度和经度方向双因子定权的地磁数据通化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108303753B (zh) | 2019-07-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110058236B (zh) | 一种面向三维地表形变估计的InSAR和GNSS定权方法 | |
CN109580003B (zh) | 一种静止气象卫星热红外数据估算近地面大气温度方法 | |
Shapiro et al. | Measurement of the Solar Gravitational Deflection of Radio Waves using Geodetic<? format?> Very-Long-Baseline Interferometry Data, 1979–1999 | |
Li et al. | Assessment of the uncertainties in temperature change in China during the last century | |
CN104280049B (zh) | 一种高精度星敏感器外场精度测试方法 | |
Wang et al. | Comparison of geographically weighted regression and regression kriging for estimating the spatial distribution of soil organic matter | |
CN103115624B (zh) | 一种基于地磁匹配的地磁日变修正方法 | |
CN110146839A (zh) | 一种移动平台磁梯度张量***校正方法 | |
Dunn et al. | Pairwise homogeneity assessment of HadISD | |
CN104406610B (zh) | 一种磁力计实时校正装置及方法 | |
Huterer et al. | No evidence for bulk velocity from type Ia supernovae | |
CN106197405B (zh) | 地磁日变影响下的惯性地磁匹配定位方法 | |
CN110146924B (zh) | 基于水波初至偏振方位的海底地震仪位置及方位反演方法 | |
CN104992054B (zh) | 基于时间序列二维化的电离层垂直总电子含量预报方法 | |
CN114253962B (zh) | 一种考虑非线性因素的区域格网速度场构建方法及*** | |
CN106959105A (zh) | 罗盘校准方法和装置 | |
Henderson et al. | A probabilistic method of paleobiogeographic analysis | |
CN108303753A (zh) | 地磁矢量场观测数据一致性校正方法 | |
CN110568458B (zh) | 一种基于gnss的电离层vtec闭环测试***及方法 | |
CN109521444A (zh) | 一种地壳运动gps水平速度场自适应最小二乘拟合推估算法 | |
Rigoti et al. | Characteristics of the equatorial electrojet determined from an array of magnetometers in N-NE Brazil | |
CN105547237A (zh) | 一种基于最小二乘法的影子定位技术 | |
CN105588569A (zh) | 一种利用影子定位的方法 | |
CN107255837B (zh) | 一种二维地形正演和改正方法 | |
Feng et al. | Data fitting and modeling of regional geomagnetic field |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190726 Termination date: 20200201 |