航空重力测量数据向下延拓球内Dirichlet方法

1、航空重力测量数据向下延拓球内dirichlet方法摘要航空重力测量是一种新型的重力测量技 术，它以飞机为载体，综合应用重力传感器、gps、高度传 感器、姿态传感器等设备测定近地空间重力加速度。由于航 空重力测量得到的是在航线高度上的重力异常值，而地重力 学惯用的是地面或大地水准面上的重力异常，要将航空重力 测量得到的数据归算到地面或大地水准面上然后才能使用。 本文主要是研究航空重力测量数据向下延拓的精度，先将地 面数据向上延拓到空中，得到航空重力测量数据，再用向下 延拓的球内dirichlet方法将空中数据延拓到地面，然后进 行精度评定，衡量航空重力测量数据向下延拓的球内 dirichlet方

2、法的可行性。关键字地球半径 泊松积分 边值问题 向下延拓 向 上延拓中图分类号p631. 1+25 文献码b 文章编号 1000-405x (2013) -5-163-21概述近几十年来，航空重力测量一直是国内外重力测量研究 的热点。航空重力测量得到的是在航线高度上的重力异常 值，而地重力学惯用的是地面或大地水准面上的重力异常， 要将航空重力测量得到的数据归算到地面或大地水准面上 然后才能使用。本文主要是研究航空重力测量数据向下延拓 的精度，方法是先将地面数据向上延拓到空中，得到航空重 力测量数据，再用向下延拓的球内dirichlet方法将空中数 据延拓到地面，然后进行精度评定，衡量航空重力测

3、量数据 向下延拓的球内dirichlet方法的可行性。2空中重力异常向下解析延拓的基本方法空中重力异常向下解析延拓的基本方法是求解poisson 方程，称为逆poisson方法。具体解法有最小二乘配置法、 迭代法、直接代表法、梯度法、点质量法等。这里主要讨论 航空重力测量数据向下延拓的球内dirichlet方法。将航空 测量数据归算到相同的航线高度上，得到测线上的重力异常 agho根据位理论的球外dirichlet问题解法，gh与地面 （过地面点的水准面）重力异常之间的关系可由poisson 积分公式表为：其中，o是半经为r的球面，r=r+h ,为空中极算点 到球面流动的距离。将空中及地面重力

4、异常转化为格网平均 值，方程式（1）可以表示成线性方程组：其中，表示已知的第i个空中重力异常，agj是 地面上规则分布的第j个重力异常，系数，空中计算点到 球面流动的距离，。这样即可由和 的坐标计算出aijo 当nm时，解线性方程组（2）就可以由空中重力异常求出地面上的重力异常gj。3实例分析3. 1 高度 100m3. 1. 1向上延拓100m航空重力异常向上延拓的公式，即为了求取100m处的具体值，需将(1.01)式转化成离 散的使用模型，即其中，这里取d。为2，d入为2。0-80°入二10。0 = 0+ (i-0. 5) d 0 入二入 + (i-0. 5) dx ；该式使用于

5、任意高度的向上延拓，因此取h为100m代 入公式(3)即可求得100m处空中点的重力异常，用程序解 决，空中重力异常数据画图如下：3. 1.2向下延拓100m考虑到乡下延拓100m时距离比较近，模型误差的影响相对较小，可以直接用dirichlet法为了求取向下延拓100m处的地面值，需将(4)式转化成具体的求解公式，即用程序解决将获得的地面数据画图得3. 2 高度 10000m3. 2. 1向上延拓10000m公式和理论与向上延拓100m都一样，只需要将100m改 成10000m即可，这里直接贴结果图3. 2.2向上延拓10000m用d让ichlet法，画图得4总结分析4. 1100m图形对比将100m的原始图形，经过向上延拓和向下延拓后得到 的图形以及它们的差值图形列出，如下：比较后可知，100m时向下延拓误差控制在0. 8毫伽范围 之内，符合精度要求。4. 210000m图形对比将10000m的原始图形，经过向上延拓和向下延拓后得 到图形以及它们的差值图形列出，如下：比较后可知，10000m时向下延拓误差达到20毫伽，延 拓造成的误差较大。5结束语根据以上讨论，得出结论，即航空重力测量数据向下延 拓的球内dirichlet方法对距离为100