海洋磁力测量的垂直空间改正

  • 格式:pdf
  • 大小:1012.70 KB
  • 文档页数:4

)
A b s t r a c t T h en e wc o n c e p t o f s p a t i a l c o r r e c t i o no f s e n s o r h a s b e e nb r o u g h t o u t . T h e s p h e r i c a l m a g n e t i c f i e l di s
, 并使观测面是一个起伏的曲面, 有时还可能
是间断的曲面( 相同海区不同时期的测量资料的衔
2 0 0 7 0 5 3 1 收稿日期: 基金项目: 国家自然科学基金( 4 0 6 7 1 1 6 1 )
作者简介: 于波, 男, 1 9 8 1年生, 博士生, 主要从事海洋磁力测量理论与数据处理研究. E- m a i l : y u b o z x p @1 6 3 . c o m
M/ 2 - 1 N / 2 - 1
T( x , y )= ∑ ∑ ( A c o s ( 2 i x / L ) c o s ( 2 j y / L )+ Δ π π i j x y
i = 0 j = 0
B c o s ( 2 i x / L ) s i n ( 2 j y / L )+C s i n ( 2 i x / L ) π π π i j x y i j x c o s ( 2 j y / L )+ D s i n ( 2 i x / L ) s i n ( 2 j y / L ) ) π π π y i j x y 3 ) (
1 ~ 4 ] 时所用的磁异常计算公式为 [ :
1 引言
在海洋磁力测量中, 由于风、 流、 浪、 涌、 潮汐、 船 速等因素的影响, 拖鱼在测量过程中是不规则上下 起伏的; 此外, 大面积的海洋磁力测量都是两艘或两 艘以上的测量船同时测量, 由于测量船的规格不同, 拖 缆 长 度、 船速的不同导致拖鱼的入水深度不 同
i = 0 j = 0 2 2 2 x
c o s ( 2 j y / L )+ B c o s ( 2 i x / L ) s i n ( 2 j y / L )+ π π π y i j x y C s i n ( 2 i x / L ) c o s ( 2 j y / L )+D s i n ( 2 i x / L ) π π π i j x y i j x s i n ( 2 j y / L ) ) π y ( 2 ) 式( 2 ) 中, x , y 为测点的平面坐标; M、 N分别为测点 A , B 、 C 、 D 为傅立叶级数; L M d x , L 和测线数; i j i j i j i j x= y = N d y , d x 、 d y 分别为点距和线距。 第二步, 把傅立叶系数代入式( 3 ) , 求出基准平 面上的磁异常。
图2 基准平面的磁异常等值线图 F i g . 2 M a g n e t i ca n o m a l yi s o l i n e s o f r e f e r e n c ep l a n e
第 1期
于波等: 海洋磁力测量的垂直空间改正
1 2 5
表2 改正前原始交叉点不符值统计( 单位: n T ) T a b . 2 S t a t i s t i c so f c r o s s i n g p o i n t d i s c r e p a n c yb r e c t i o n 不符值 Z 1 Z 2 Z 3 Z 4 Z 5 Z 6 图3 观测曲面变化立体图 F i g . 3 S t e r e o s c o p i cd i a g r a mo f c h a n g eo f u n d u l a t i n go b s e r v a t i o ns u r f a c e … Z 9 8 Z 9 9 J 1 J 2 J 3 J 4 J 5 J 6
V E R T I C A LF R E E A I RC O R R E C T I O NI NMA R I N EMA G N E T I CS U R V E Y
1 ) 2 ) 1 ) 1 ) 1 ) Y uB o ,Z h a i G u o j u n ,L i uY a n c h u n ,B i a nG a n g a n dB i a nG u a n g l a n g
图1 垂直空间改正示意图 F i g . 1 S k e t c ho f t h ev e r t i c a l f r e e a i r c o r r e c t i o n
值线间隔为 1 0n T , 最大磁异常为 1 5 0n T , 最小磁异 1 0 0n T 。仿真区域中, 主测线共 1 0 0条, 检查 常为 - 线 6条。根据测线轨迹仿真拖鱼在水中运动起伏的 x 、 y 为 高度, 观测曲面立体图如图 3所示。图 3中, 测线的平面坐标, z 为拖鱼相对于基准平面的起伏 高度变化量, 最大变化量约为 4 . 1 1m 。 采用调和级数法对观测曲面上的磁场值进行垂 直空间改正, 将曲面上的测点归算到基准平面上, 改 正前、 后的磁异常等值线比较如图 4所示( 实线表 示改正前 的 等 值 线, 虚线表示改正后的平面等值 线) 。从图 4中可以看出, 改正前与改正后的等值 线在磁异常变化剧烈的地方差别较大, 并且改正后 的等值线较改正前光滑。
[ 1 , 2 ]
T= T- T T T ( 1 ) Δ d- s- o T为海洋地磁异常值, T为地磁场总场测量 式中 Δ 值, T T T o 为地磁正常场值, d 地磁日变改正值, s为船 磁改正值。 因此测量得到的磁场值经过式( 1 ) 中的各项改 正后只是平面位置的变换, 垂直位置并没有改变, 即 处理后的测量数据面仍然是曲面。磁场数据不在同
第一步, 根据观测曲面上的测量值代入式( 2 ) , 建立线性方程组, 求出不同的傅立叶系数:
M/ 2 - 1 N / 2 - 1
T( x , y , h )= ∑ ∑ Δ
2 1 2 e x p2 h ( i/ L+ j / L )/ A c o s ( 2 i x / L ) [ π ]( Δ π y i j x
(
1 ) D e p t .o f h y d r o g r a p h ya n dM a p p i n g ,D a l i a nN a v a l A c a d e m y ,D a l i a n 1 1 6 0 1 8 ) I n s t i t u t e o f H y d r o g r a p h i c S u r v e y i n ga n dC h a r t i n g,T i a n j i n 3 0 0 0 6 1 2
第2 8 卷 第1 期 2 0 0 8年 2月
大地测量与地球动力学 J O U R N A LO FG E O D E S YA N DG E O D Y N A M I C S
V o l . 2 8N o . 1 F e b . , 2 0 0 8
文章编号: 1 6 7 1 5 9 4 2 ( 2 0 0 8 ) 0 1 0 1 2 3 0 4
1 2 4
大地测量与地球动力学
2 8卷
一平面上会出现交叉点不符值过大的问题, 即测量 精度较差( 测量精度的评价指标是交叉点不符值中 误差 的
[ 2 ]
y R s i n 2 I s i n A ′ ] 3
( 4 )
) 。目前, 仅仅通过“ 调平” 的方法
[ 5 , 6 ]
降低交
式中 μ m 为磁矩; R为球体中心埋 0 为真空磁导率; 深; x 、 y 为空间任意一点的平面坐标; I 为磁化倾角; A ′ 为磁化偏角。 本文所用的球体参数如表 1所示。
[ 7 ]
。本文仿真了一个起伏曲面的球体磁场, 采用
[ 8 , 9 ]
调和级数法
对仿真数据进行计算, 研究了进行
垂直空间改正对提高测量精度的影响量级。
2 方法原理
图 1中, 为基准平面, h为 Γ Γ Δ s为观测曲面, s ′ 观测曲面相对于基准平面的高度变化量, 观测曲面 h> 0 ; 在基准平面下 上的点位于基准平面上面时, Δ 面时, h < 0 。 Δ
表1 仿真球体参数 T a b . 1 P a r a me t e r s o f e mu l a t i v es p h e r e 球体半径 ( m ) 5 0 磁化强度 ( A/ m ) 1 0 0 磁化倾角 ( ° ) 4 5 磁化偏角 ( ° ) 4 5
叉点不符 值 中 误 差 以 满 足 测 量 精 度, 这是不合理
,v e r t i c a l f r e e a i r c o r r e c t i o ni s c o m p u t e df r o mh a r m o n i c s e r i e s .T h e s i m u l a t e do nt h eu n d u l a t i n go b s e r v a t i o ns u r f a c e m e a ns q u a r ee r r o r o f c r o s s i n g p o i n t d i s c r e p a n c yi sa l s oc o m p u t e d .I t i sc o n c l u d e dt h a t t h em e a ns q u a r ee r r o r o f c r o s s i n g p o i n t i s i m p r o v e da b o u t 1 0n Ta f t e r t h ev e r t i c a l f r e e a i r c o r r e c t i o n . K e yw o r d s : m a r i n em a g n e t i cs u r v e y ;v e r t i c a l f r e e a i r c o r r e c t i o n ;c r o s s i n g p o i n t d i s c r e p a n c y ;m e a ns q u a r ee r r o r ; s u r v e y i n ga c c u r a c y 接处) 。目前国内外在进行海洋磁力测量数据处理