球冠谐分析方法进行球面卫星重力场的向上向下延拓.pdfVIP

球冠谐分析方法进行球面卫星重力场的向上向下延拓 尹航王喜臣徐宝慈贾建秀 吉林大学地探学院 吉林省长春市 130026 G。V。Haines于1985年首次提出的适用于在球坐标中求解局部区域拉普拉斯方程边值问题的解法，即球冠 谐分析方法，该方法除具有与球谐分析方法一样的优点，即具有坚实的物理基础；满足位理论的物理限制，进行球 冠谐分析外，亦能进行位场的向上向下延拓，可方便的表示位场的三维结构：可提取不同波长位场异常，达到分 离场的目的。它与传统的向上向下延拓方法的区别在于：传统的向上向下延拓方法通常是建立在适应于平面延拓 的矩坐标上，而冠谐分析的方法则在较大的区域内，更适用于真实反映地球表面的一个有限区域的球坐标上进行 曲面的向上向下延拓。因为球谐分析中的M和N(边值问题本征值)是正整数，而冠谐分析中，则用整阶M但 非整次Nk的缔合勒让德函数来表达球函数的模型。 在无源区域，重力场可表示为满足拉普拉斯方程的重力位的梯度，在球坐标中，拉普拉斯方程的解为： O) (1) cos聊p+办?sinm伊)p：’(cos V=∑∑W7=口∑(形)”1∑Lg。111 式中：r为研究点至地心的距离，a为地球半径，口为余纬度，c,o为经度，g：’与办：’为球谐系数。 实际分析中，本征值m和，?是由边界条件决定的，重力位关于精度的连续性要求m是整数，而位于o：0“ 和o：180·，处的条件决定胛也是实整数。在关于p的常微分方程中，具有整数的m和门的解是缔合勒让德多项 0)。而对地球表面一部分的球冠的边界条件，在经度及纬度p=00处与全球分析中的相同，但在 式p?(cos 0=Oo处边界上，重力位及其导数必须满足 V(r，Oo，妒)=f(r，伊) (2) 型：如缈)…7 (3) 00 (2)式与(3)式右端项的f与g均为与目无关的函数。Haines担¨31已经说明，方程(2)可以通过选择基函数使得来满 足。这可以通过决定 —clp：(c—osO)：0 (4) dO ttj根n来完成，其中m=0，1，2，3，…。同样，方程(3)可以通过使基函数r?’在球冠的边界上为零来满足。也 就是通过决定rl值来完成，即 。 p：’(cos0)=0 (5) ·680· 能使方程(4)和(5)成立的n值一般是非整阶的，用‰表示。这样，对于一个给定的聊值，通过求(4)式雨 (5)式所决定的两组解，形成了两组基函数。每一组基函数在球冠上是正交的，他与传统的球谐分析中勒让德多项 式秒=00到目=1800区间上正交相同． 但两组基函数之间不正交，这样，重力位可以表示为 七 矿 = (6) ∑口(∥加卅1(簖cosm6p+向?’sinm(p)p,il．㈣)(cose) ‰∑㈨ m=O (6)式与(1)式类似，但(1)式中n的求和已由新的角标k的求和所代替，k值被用于高斯系数的下角标，聆值由 ‰(聊)所取代，甩+1由／-／t(m)+1所取代。 给定‰和m值，可以计算pZ(肿)(．cosO)(见参考文献)来完成上式。 由于重力异常是重力位的梯度，根据g=一dy／咖便可得到球冠谐和分析中重力异常的表达式： k II一 (7) g ‰∑脚 in=O 其中r为研究点至地心的距离，a为地球半径，口为余纬度，缈为经度，g≯与办：’为球谐系数，此式即为计算 不同高度重力异常的基本公式，当r不同时，亦即实现了不同高度延拓的目的。由此可见，利用冠谐分析的方法 进行球面卫星重力场的向上向下延拓的过程也就是