高精度航空重力测量系统分项指标设计分析
1.
2.
3.
Design and analysis of sub-index of high precision Airborne Gravity System
1.
2.
3.
责任编辑: 王萌
收稿日期: 2018-09-11 修回日期: 2019-01-14 网络出版日期: 2019-04-20
基金资助: |
|
Received: 2018-09-11 Revised: 2019-01-14 Online: 2019-04-20
作者简介 About authors
张虹(1963-),女,辽宁锦县人,副教授,北京信息科技大学,现主要从事地球物理测量方法技术、数理统计应用技术等领域的研究工作。Email:
为研制高精度航空重力仪,开展了高精度航空重力仪分项指标设计与分析,为高精度航空重力测量系统设计提供依据。笔者基于航空重力测量的数学模型,归纳出影响航空重力测量的主要因素;结合航空重力测量的相关理论公式,使用理论模型,推导出各影响因素的误差模型;开展了各项误差的分析研究,并通过设定合理的分项指标精度来有效地控制分项误差,确保航空重力测量精度优于0.6×10 -5 m/s 2。
关键词:
In order to meet the needs of development of high-precision airborne gravimeter , the sub-index design and analysis of high-precision airborne gravimeter were carried out, which provided the basis for the design of high-precision airborne gravimeter system. Based on the mathematical model of airborne gravimetry, this paper summarizes the main errors affecting airborne gravimetry, uses the theoretical model formula, deduces the error models of various influencing factors, and carries out the analysis and study of various errors, and sets reasonable sub-index by combining the relevant theoretical formula of airborne gravimetry. Accuracy can effectively control the sub item errors and ensure the accuracy of airborne gravity measurement is better than 0.6×10 -5 m/s 2.
Keywords:
本文引用格式
张虹, 屈进红, 姜作喜, 王萌, 李行素.
ZHANG Hong, QU Jin-Hong, JIANG Zuo-Xi, WANG Meng, LI Xing-Su.
0 引言
我国从2006年开始进行航空重力测量技术研究和航空重力仪研制的工作。到目前为止,已经研制出三轴惯性稳定平台型和捷联式航空重力仪样机,测量精度达到1.0×10-5 m/s2,取得了突破性进展[6]。为了进一步提高我国自主研制航空重力勘查系统的测量精度,使其达到国际先进水平,“十三五”期间开始了新一轮航空重力测量系统的研制工作。
航空重力测量精度取决于传感器误差、传感器姿态误差、卫星导航误差等各方面[7],为了研制测量精度优于0.6×10-5 m/s2的航空重力测量系统,须控制好各分项指标的误差,从而保证航空重力测量精度满足要求。笔者通过建立航空重力(标量重力)测量的误差模型,分析各分项误差对航空重力测量精度的影响贡献,由此通过设定合理的分项指标有效控制分项误差,实现航空重力测量精度优于0.6×10-5 m/s2的研究目标。
1 航空重力标量测量的数学模型
假设在重力传感器处于垂向时,航空重力标量测量数学模型[10]可表示为:
其中: ga是航空自由空间重力值,m/s2; gb是停机坪处的基点重力值,由陆地重力仪从重力基准点联测得到,m/s2;fZ是重力传感器测点观测值,m/s2;
为了获得航空自由空间重力值,须对垂向重力传感器观测值进行载体垂向加速度改正、厄特渥斯改正、水平加速度改正、零漂改正、正常重力场改正和高度改正。
航空重力测量误差可表示为各项改正误差传递[10],由式(1)可得:
其中:dδgb为停机坪处重力基点误差;dδfz为重力传器观测误差;d$\delta_{\ddot{Z}}$为载体垂向加速度测量误差;dδaE为厄特渥斯改正误差;dδaH为水平加速度改正误差;dδaF为高度改正误差;dδaK为零漂改正误差;dδaI为重力仪和差分GPS的随机误差;dδγ0为椭球面上的正常重力场改正误差。
2 航空重力测量各项误差分析
2.1 载体垂向加速度测量误差
载体飞行高度、测量精度直接影响着航空重力的测量精度,可以说决定了航空重力的测量精度[9]。
当沿着某一测线进行航空重力测量飞行时,假设:飞行高度变化幅值为A,飞行高度变化周期(T)为60 s,且飞行高度按正弦规律变化,因此载体飞行高度(Z)变化可表示为[10]:
其中:ω=2π/T=2×3.14/60≈0.1 s-1。
从飞行高度的位移可得到载体垂向加速度为:
$\ddot{Z}=\omega^{2}Asin\omega t$(4)
当测得的飞行高度幅度A=1 mm时,由此可得载体最大垂向加速度为:
$\ddot{Z}_{max}=0.1^{2} · 0.001=1 \times 10^{-5}m/s^{2}$
由此可得,高程1 mm的误差可引起载体垂向加速度最大误差1×10-5 m/s2,同时载体垂向加速度误差与飞行起伏周期T2反比。因此,飞行起伏周期60 s时GPS高程测量误差须小于1 mm,而当飞行起伏周期100 s以上时GPS高程测量误差可放宽到2~3 mm,才能保证载体垂向加速度最大误差控制在1×10-5 m/s2之内。
2.2 同步精度引起的测量误差
仍采用2.1的飞行方法,假设GPS记录的载体加速度滞后或超前Δt s,因此由GPS解算得到的载体垂向加速度可表示为[7]:
由此可得到由于采样不同步造成载体垂向加速度与GPS解算垂向加速度之间的误差:
由于Δt≪t且Δt值很小,所以式(6)可简化为:
将A=1、ω=0.1代入式(7),可得:
当Δt=±1 ms时,Δ$\ddot{Z}$max=±0.1(10-5 m/s2);
当Δt=±10 ms 时,Δ$\ddot{Z}$max=±1.0(10-5 m/s2)。
由此可知:当航空重力仪与GPS时间同步误差在1 ms时,将产生0.1×10-5 m/s2的重力测量误差,因此设计的航空重力仪与GPS的同步精度应好于 1 ms。实际上,同步精度可达微秒级,其造成的误差可以忽略不计。
2.3 厄特渥斯改正误差
因载体的运动而使航空重力仪产生了附加离心力,这种影响称为厄特渥斯效应,消除该影响的方法即为厄特渥斯改正。载体飞行速度的测量精度直接影响了厄特渥斯的改正精度[16]。
以RN+h=RM+h=R代入式(9),即近似情况下的厄特渥斯改正公式:
其中:vN、vE、v分别为载体速度的N向、E向分量和水平速度,单位m/s;φ、h分别为测点的地理纬度和载体的飞行高度,单位分别为rad、m; R、RM、RN分别为地球半径、参考椭球上各点的子午圈半径和卯酉圈半径,单位m; ω是地球自转角速度,单位rad/s; A是飞行航向角(或方位角),单位rad。
从式(10)可以看出,厄特渥斯改正的计算精度主要取决于载体速度、大地纬度以及方位角的确定精度。按误差传播定理可得厄特渥斯改正误差:
取地球半径R为6 378 137 m,平均自转角速度ω为7.292 1×10-5 rad/s。
1) 位置误差引起的厄特渥斯改正误差[7]
由式(11)右端第三项先分析测点纬度误差对δaE的影响,令|sinA|=1,|sinφ|=1,计算测点纬度误差对厄特渥斯改正影响(如表1)。目前利用GPS测量纬度的精度远优于0.1’,不同航行速度下该改正项精度将优于0.07×10-5 m/s2,通常可以忽略不计。
表1 纬度误差对厄特渥斯改正的影响 (单位:10-5 m/s2)
Table 1
速度/(km/h) 纬度误差/(‘) | 100 | 200 | 300 | 400 | 500 | 600 |
---|---|---|---|---|---|---|
0.05 | 0.0059 | 0.0118 | 0.0177 | 0.0236 | 0.0295 | 0.0354 |
0.1 | 0.0118 | 0.0236 | 0.0354 | 0.0471 | 0.0589 | 0.0707 |
0.2 | 0.0236 | 0.0471 | 0.0707 | 0.0943 | 0.1178 | 0.1414 |
2) 速度误差引起的厄特渥斯改正误差
表2 速度误差对厄特渥斯改正的影响 (单位:10-5 m/s2)
Table 2
速度误差/(m/s) 飞行速度/(m/s) | 0.005 | 0.01 | 0.015 | 0.02 | 0.025 | 0.03 |
---|---|---|---|---|---|---|
50 | 0.0808 | 0.1615 | 0.2423 | 0.3230 | 0.4038 | 0.4846 |
60 | 0.0823 | 0.1647 | 0.2470 | 0.3293 | 0.4116 | 0.4940 |
70 | 0.0839 | 0.1678 | 0.2517 | 0.3356 | 0.4195 | 0.5034 |
80 | 0.0855 | 0.1709 | 0.2564 | 0.3419 | 0.4273 | 0.5128 |
90 | 0.0870 | 0.1741 | 0.2611 | 0.3481 | 0.4352 | 0.5222 |
3) 航向角误差引起的厄特渥斯改正误差
由式(11)右端第二项分析航向误差对δaE的影响。令|cosA|=1,|cosφ|=1,计算出航向误差对厄特渥斯改正的影响(如表3)。
表3 航向误差对厄特渥斯改正的影响 (单位:10-5 m/s2)
Table 3
飞行速度/(m/s) 航向误差/(‘) | 200 | 220 | 300 | 400 | 500 | 600 |
---|---|---|---|---|---|---|
0.1 | 0.0236 | 0.0259 | 0.0354 | 0.0471 | 0.0589 | 0.0707 |
0.5 | 0.1178 | 0.1296 | 0.1768 | 0.2357 | 0.2946 | 0.3535 |
1 | 0.2357 | 0.2593 | 0.3535 | — | — | — |
目前航向精度优于1’,航向误差将引起最大0.36×10-5 m/s2厄特渥斯的改正误差。同时由式(11)可以看出在EW向上航向误差对厄特渥斯改正影响最小,在SN航向上航向误差对其影响最大。
4) 航空重力厄特渥斯改正综合误差
测点纬度误差对厄特渥斯改正影响小,可以忽略不计。载体在赤道附近EW向飞行时,速度测量误差对厄特渥斯改正影响最大,但航向误差影响接近0,厄特渥斯改正的精度优于0.35×10-5 m/s2。载体在赤道附近SN向飞行时,航向误差对厄特渥斯改正影响最大,速度测量误差影响很小,厄特渥斯改正的精度优于0.36×10-5 m/s2。
因此,综合考虑到测点纬度误差、速度测量误差和航向误差对厄特渥斯改正的影响,厄特渥斯改正最大误差参考值取为(0.20~0.30)×10-5 m/s2,此时位置纬度精度应优于0.1’、速度精度应优于0.02 m/s和航向精度应优于1’。
2.4 水平改正误差
图1
图1
平台非水平时重力传感器感知的重力值示意
Fig.1
The gravity value perceived by the gravity sensor when the platform is not horizontal
其中:g为当地重力值,单位m/s2;aH为水平加速度,单位m/s2;
通常情况下,稳定平台(或姿态)偏离角θ很小,于是水平加速度改正值可简化为:
由式(13)可得,平台偏离角(或姿态角)测量误差dθ引起水平加速度改正误差dδaH可表示为:
取g=9.8 m/s2,aH=0.1 m/s2,利用式(14)计算因平台姿态角测量误差引起的水平加速度改正误差(如表4)。
表4 平台测量误差对水平加速度改正的影响 (单位:10-5 m/s2)
Table 4
θ dθ | -30″ | -1’ | -2’ | -5’ |
---|---|---|---|---|
1″ | -0.0235 | -0.0229 | -0.0215 | -0.0173 |
5″ | -0.1177 | -0.1143 | -0.1074 | -0.0867 |
10″ | -0.2355 | -0.2286 | -0.2148 | -0.1733 |
15″ | -0.3532 | -0.3429 | -0.3221 | -0.2600 |
20″ | -0.4710 | -0.4572 | -0.4295 | -0.3466 |
由表4可见,在外加水平加速度0.1 m/s2(测量时通常小于该值)时,平台测量角(或姿态角)误差为10″(角秒)所造成的水平加速度改正误差小于0.24×10-5 m/s2。
2.5 正常重力场改正改正误差
国际大地测量协会(IAG)1980年[9]大地测量参考系统中的正常重力场改正γ0(m/s2):
其中:φ是观测点地理纬度,单位rad。
由式(15)可得到正常重力场值改正误差dγ0(m/s2):
式中:
目前差分GPS定位的水平定位精度好于1.0 m,由此引起的正常重力场测量误差小于0.001×10-5 m/s2,可以忽略不计。
2.6 高度改正误差
航空重力的空间(高度)改正δaF(m/s2)的计算公式为[14]:
其中:φ是观测点的地理纬度,单位rad;h是观测点的飞行高度,单位m;N是大地水准面高,单位m;Δh是大地高的偏心改正值,单位m。
由式(17)可得航空重力高度改正误差dδaF(m/s2)为:
其中:
2.7 零漂改正误差
零漂改正误差取决于重力仪零漂线性度[17],一般重力仪漂移非线性引起的改正误差约0.10×10-5 m/s2,因此,重力仪零漂改正误差约小于0.1×10-5 m/s2。
2.8 重力基点误差
停机坪处的基点重力值gb,由陆地重力仪从重力基准点联测得到[4],其测量精度一般优于0.05×10-5 m/s2,因此,航空重力测量误差可以不考虑基点引点误差的影响。
2.9 重力传感器观测误差
航空重力传感器观测值fZ来自于重力传感器的测量值。目前重力传感器有两种类型:通用的加速度计和弹簧类的重力传感器,其测量精度一般都优于0.3×10-5 m/s2,并可通过地面静态测量精度来衡量[5]。因此,航空重力传感器静态测量精度可优于0.3×10-5 m/s2。
3 航空重力仪各分项精度要求
忽略重力基点误差、正常重力场改正误差、重力仪和差分GPS的随机误差的影响[19],利用误差传递公式计算式(2)各项误差的总误差,并将前文分析的各项误差值代入计算,即获得航空重力测量精度为:
因此,要实现航空重力测量精度优于0.6×10-5 m/s2,分项技术设计指标应优于表5所示。在保证总精度指标的前提下,可以根据需要适当调配各分项精度指标。
表5 航空重力测量系统分项指标精度要求
Table 5
名称 | 重力传感器 观测精度(静态) | 载体高程 测量精度 | 载体位置 测量精度 | 载体速度 测量精度 | 稳定平台角(或姿 态角)测量精度 | 同步精度 |
---|---|---|---|---|---|---|
指标 | 优于0.3(10-5 m/s2) | 优于3 mm | 优于1 m | 优于0.02 m/s | 优于10″ | 优于1 ms |
4 总结
为满足新研制航空重力测量系统总体设计指标的要求,针对影响航空重力测量精度的分项指标,建立航空重力各分项误差模型,并开展了分析研究,总结归纳出航空重力测量系统所涉及的各分项指标的精度要求,为航空重力测量系统设计奠定了基础。主要成果为:
1) 建立了航空重力测量的数学模型,给出了影响航空重力测量精度的主要因素;
2) 逐一建立了航空重力各影响因素的数学模型,或使用理论模型公式,推导出各影响因素的误差计算公式。
3) 通过设定一定的分项指标,计算出各分项所产生的误差大小,分析各因素对测量结果的影响程度;
4) 总结出要实现航空重力测量精度优于0.6×10-5 m/s2各分项误差的控制范围,从而得出引起误差的各分项指标的精度要求。
总之,要想设计出高精度的航空重力测量系统,各分项指标的控制十分重要,论文研究成果一定程度上为读者开展航空重力测量系统设计提供了各项指标的分析方法和技术支撑数据,同时为十三五国家重点研发计划项目航空重力测量系统设计和研制打下坚实的理论基础。
参考文献
我国航空重磁技术现状与发展趋势
[J].<FONT face=Verdana>本文概要介绍了我国航空磁力测量方法技术现状,近年来测量技术与解释方法方面的最新进展.首次介绍了我国引进国外航空重力测量系统以来的最新测量技术成果.分析了我国航空物探技术与国际上先进国家的差距,并对我国航空物探的发展趋势进行了分析与预测.</FONT>
The present situation and development of airborne gravity and magnetic survey techniques in China
[J].
航空重力测量系统研究进展
[J].<p>简要回顾了航空重力测量的发展过程,在介绍了航空重力测量基本原理的基础上,系统地叙述了航空标量重力测量系统的研发历程、研究现状及发展动态。最后,阐述了我国航空重力测量的研究现状和未来的发展方向。</p>
The advances in the study of the airborne gravimetry system
[J].
A comparison of stable platform and stapdown airborne gravity
[J].
DOI:10.1007/s001900000082
URL
[本文引用: 1]
To date, operational airborne gravity results have been obtained using either a damped two-axis stable platform gravimeter system such as the LaCoste and Romberg (LCR) S-model marine gravimeter or a strapdown inertial navigation system (INS), showing comparable accuracies. In June 1998 three flight tests were undertaken which tested an LCR gravimeter and a strapdown INS gravity system side by side. To the authors' knowledge, this was the first time such a comparison flight was undertaken. The flights occurred in Disko Bay, off the west coast of Greenland. Several of the flight lines were partly flown along existing shipborne gravity profiles to allow for an independent source of comparison of the results. The results and analysis of these flight tests are presented. The measurement method and error models for both the stable platform and strapdown INS gravity systems are presented and contrasted. An intercomparison of gravity estimates from both systems is given, along with a comparison of the individual estimates with existing shipborne gravity profiles. The results of the flight tests show that the gravity estimates from the two systems agree at the 2 3 Gal level, after the removal of a linear bias. This is near the combined noise level of the two systems. It appears that a combination of both systems would provide an ideal airborne gravity survey system, combining the excellent bias stability of the LCR gravimeter with the higher dynamic range and increased spatial resolution of the strapdown INS.
State-of-the-art Airborne Gravimetry [A]. In Gravity and Geoid: Joint Symposium of the International Gravity Commission and the International Geoid Commission, Graz, Austria: Springer,
The GT-1A mobile gravimeter
[A].
SGA-WZ: A New Strapdown Airborne Gravimeter
[J].
航空重力测量理论方法及应用研究
[D].
Theory, methods and applications of airborne gravimetry
[D].
航空重力测量数学模型及其测量精度分析
[J].
Mathematical models and accuracy evaluation for the airborne gravimetry
[J].
基于捷联惯导系统/差分全球定位系统的航空重力测量技术
[J].
Technology of airborne gravimetry based on SINS/DGPS
[J].
Theory, method and application of the airborne gravity prospecting
[M].
基于经验模态分解的航空重力测量动态误差分离
[J].
Dynamic errors separation of airborne gravimetry based on empirical mode decomposition
[J].
Accuracy analysis of airborne gravity when gravimeter parameters are identified in flight
[J].
航空重力数据的等波纹FIR低通滤波试验
[J].
Experimental researches on FIR low pass filter based on equiripple
[J].
航空重力数据窗函数法FIR低通数字滤波试验
[J].
利用Kalman平滑技术估算航空重力测量中的载体垂直加速度
[J].
Estimation of the vertical acceleration for the airborne gravimetry using Kalman smoothing
[J].
航空矢量重力测量中科里奥利加速度对测速精度要求分析
[J].
Analysing the accuracy of measuring velocity for the airborne vectror gravity survey according to coriolis acceleration
[J].
相关性分析在GT航空重力仪零漂改正中的应用
[J].
The correlation analytical method and its application to GT aerogravimeter drift correction
[J].
Control theoretic approach to inertial navigation systems
[J].
GT-1A inertial gravimeter system design consideration and results of flight tests
[C]//
航空重力测量中载体运动加速度的确定
[J].
Determination of acceleration of carrier motion in aeronautical gravity measurement
[J].
航空矢量重力测量误差模型的建立及初步分析
[J].
Establishment and preliminary analysis of error model of aeronautical vector gravity measurement
[J].
航空重力测量精度估计
[J].
Accuracy estimation of airborne gravimetryJ]
基于GPS双差相位率确定运动载体的速度和加速度
[J].
Determination of velocity and acceleration of moving carrier based on GPS double difference phase rate
[J].
/
〈 |
|
〉 |
