-
20世纪70年代，卫星测高技术的出现以及多个卫星测高任务的顺利实施为人们获取全球海域高精度、高分辨率的海面高数据提供了极大的便利。利用卫星测高数据推求海域重力异常是卫星测高在大地测量研究中的主要应用，国内外许多学者对此进行了大量的研究，取得了丰富的研究成果[1-9]。其中，逆Vening-Meinesz公式法的输入量为垂线偏差，而垂线偏差是由测高观测值的一次差分求得，可以消除与地理位置相关的径向轨道误差和长波海面地形等类似系统误差，同时含有丰富的重力场高频成分，非常有利于高分辨率海洋重力场的恢复。因此，逆Vening-Meinesz公式在测高重力反演中的应用较为广泛[10]。
在计算点本身及其附近区域计算重力异常时，积分区域包含计算点，导致逆Vening-Meinesz公式中的积分奇异，本文称该积分区域为中央区。Hwang[5]将中央区视为圆域，推导出了中央区重力异常的计算公式；常晓涛等[9]将中央区视为方形域，推导出了中央区效应的计算公式，指出中央区效应与奇异区内垂线偏差分量的梯度及奇异区面积的大小有关。由于实际数据通常为网格化分布，受子午线收敛的影响，中央区更接近于矩形域，因此，将中央区视为圆域和方形域的处理方法与数据真实分布情况并不相符，由此产生的误差在高精度重力反演中能否忽略值得深入研究。鉴于此，本文推导出了视中央区为矩形域时的重力异常计算公式，并对导出公式和圆域与方形域下传统公式的误差进行了分析比较。
HTML
-
略去详细的推导过程，根据文献[5]，直接写出逆Vening-Meinesz公式如下：
式中，γ0=979.8Gal为地球平均重力；αQP代表流动点Q至计算点P的方位角；ψPQ为P、Q之间的球面距离；ξQ和ηQ分别代表流动点Q处的垂线偏差在子午圈方向和卯酉圈方向上的分量；H′(ψPQ)为积分核函数，其表达式为：
为对中央区进行计算，首先定义以计算点P为原点的局部切平面直角坐标系，如图 1所示，x轴指向北极，y轴指向东，且xy平面在P点与地球表面相切，Q(x, y)为中央区内流动的积分点，l为计算点与流动点的平面距离。由图 1可知：
Figure 1. The Sketch Map of the Local Coordinate System
式中，$l = \sqrt {{x^2} + {y^2}} $。
当积分区域比较小时，球面角距ψPQ可以简化为：
则式(1)中的积分核函数可以简化为：
将式(3)和式(5)代入式(1)，顾及平面近似下积分面元为R2dσ=dxdy，可得中央区重力异常的计算公式为：
可以看出，式(6)中的积分在计算点P处奇异。为解决这一问题，文献[5]将中央区视为圆域，将中央区内的垂线偏差分量ξQ和ηQ展开为泰勒级数形式：
式中，${\xi _x} = \frac{{\partial \xi }}{{\partial x}}, {\xi _y} = \frac{{\partial \xi }}{{\partial y}}$，${\xi _{xx}} = \frac{{{\partial ^2}\xi }}{{\partial {x^2}}}, {\xi _{yy}} = \frac{{{\partial ^2}\xi }}{{\partial {y^2}}}$，${\xi _{xy}} = \frac{{{\partial ^2}\xi }}{{\partial x\partial y}}, {\eta _x} = \frac{{\partial \eta }}{{\partial x}}$，${\eta _y} = \frac{{\partial \eta }}{{\partial y}}, {\eta _{xx}} = \frac{{{\partial ^2}\eta }}{{\partial {x^2}}}$，${\eta _{yy}} = \frac{{{\partial ^2}\eta }}{{\partial {y^2}}}, {\eta _{xy}} = \frac{{{\partial ^2}\eta }}{{\partial x\partial y}}$。
导出的中央区重力异常计算公式为：
式中，Δx、Δy分别表示x方向和y方向的网格间距。
文献[9]将中央区视为方形域，将垂线偏差分量同样表示为式(7)和式(8)，导出的重力异常计算公式为：
式中，s为方形域长度的一半。
-
首先将中央区垂线偏差分量表示为双二次多项式插值形式，之后利用文献[11-13]提出的非奇异变换对逆Vening-Meinesz公式中的奇异积分进行了处理，推导出了中央区重力异常的精密计算公式。
-
如图 2所示，设中央区为σ[－a < x < a, －b < y < b]，其中a=1，b=cosφ。中央区共包含4个网格单元，9个网格节点。将垂线偏差子午分量ξQ、卯酉分量ηQ表示成双二次多项式插值形式：
Figure 2. The Innermost Area when the Components of Deflections of the Vertical are Expressed as Bi-quad- ratic Polynomials
为确定出待定系数αij和βij，将式(11)和式(12)改写为：
则有：
将图 2中网格节点处的垂线偏差子午分量ξij=ξ(i, jcosφ)和卯酉分量ηij=η(i, jcosφ)(i, j=－1, 0, 1)作为插值条件代入式(13)和式(14)，可得：
式(17)、(18)可以简化为：
其中，
系数矩阵A非奇异、可逆，并且：
因此有：
将式(23)和式(24)代入式(15)和式(16)即可确定出待定系数αij和βij。
-
为推导方便，记式(6)中的奇异积分为：
其中，
将式(11)和式(12)分别代入式(26)和式(27)，并考虑到奇偶函数的积分性质，则有：
由计算点P向右上顶点作连线分右上象限为σ1[0 < x < 1, 0 < y < bx]和σ2[0 < x < b－1y, 0 < y < b](见图 2)。
则三角形σ1映射为矩形σ′1[0 < x < 1, 0 < k < b]，三角形σ2映射为矩形σ′2[0 < λ < b－1, 0 < y < b]。将式(30)和式(31)分别代入式(28)和式(29)，并注意到两式中的被积函数均为偶函数，因此有：
可以看出，在非奇异变换下，原来含有x、y两变量的二维积分转换为只含k(λ)变量的一维积分，利用Mathematica计算机代数系统[14-15]可求得：
因此，式(6)中的奇异积分可以表示为：
将式(38)代入式(6)，可得中央区包含4个网格时，垂线偏差双二次多项式插值表示下重力异常的计算公式为：
若将积分区域取为$\sigma '\left[{-\frac{1}{2} < x < \frac{1}{2}, -\frac{1}{2}b < y < \frac{1}{2}b} \right]$，则计算点P所在的1个网格对重力异常的贡献可采用类似方法得到：
2.1. 中央区垂线偏差分量双二次多项式插值表示
2.2. 中央区重力异常的解析计算公式
-
选定中国南海12°N~18°N，112°E~118°E海域作为试算区，以EGM2008地球重力场模型计算得到的2′×2′分辨率的垂线偏差数据作为背景场，计算了该区域的中央区效应。
试算区共179×179个网格，假定中央区为4′×4′，即包含4个网格，利用式(9)、式(10)、式(39)计算了该区域的中央区效应，结果分别为Δg1、Δg2、Δg3, 计算结果之间的比较情况如表 1所示。Δg1－Δg3、Δg2－Δg3计算结果分别如图 3和图 4所示。
差值 最小值 最大值 标准差 Δg1－Δg3 -1.229 1.276 0.142 Δg2－Δg3 -1.192 1.156 0.129 Δg2－Δg1 -0.199 0.128 0.024 Table 1. Comparisons of the Innermost Effects Computed by Three Methods/mGal
Figure 3. The Sketch Map of Δg1－Δg3
Figure 4. The Sketch Map of Δg2－Δg3
由表 1、图 3和图 4可以看出，当中央区包含4个网格时，文献[5]给出的式(9)与本文导出的式(39)计算得到的中央区重力异常差值的标准差均为0.142 mGal，最大值达1.276 mGal；文献[9]给出的式(10)与本文导出的式(39)计算得到的中央区重力异常差值的标准差均为0.129 mGal，最大值达1.156 mGal。
-
为提高利用逆Vening-Meinesz公式反演中央区重力异常的精度，本文将中央区垂线偏差表示成双二次多项式插值形式，推导出了与实际数据分布更为相符的矩形域中央区效应精密计算公式，并利用模型数据计算分析了导出公式和传统公式的误差。研究表明：
1) 非奇异变换可化逆Vening-Meinesz公式中的奇异积分为非奇异积分，可在其他地球重力场奇异积分问题中推广使用。
2) 低纬度试算区2′×2′分辨率垂线偏差数据下的中央区效应分析表明，当中央区包含4个网格时，传统公式与本文导出公式计算结果差值的均方差和标准差大于0.1 mGal，最大值大于1 mGal，需要在高精度测高重力计算中加以考虑。
DownLoad:
