大地测量Stokes边值问题的自然边界元解法武汉理工大学学报.PDFVIP

第 34 卷 　第 2 期 ( ) Vol . 34 　No . 2 武汉理工大学学报 交通科学与工程版 20 10 年 4 月 J our nal of Wuhan U niver sit y of Technolo gy Ap r . 20 10 ( Tran sport ation Science Engineering) 大地测量 Sto ke s 边值问题的自然边界元解法 1 ,2) 2) 1) 3) 唐元义 　许厚泽 　胡清峰 　骆有隆 ( 1) ) 湖北大学数学与计算机科学学院 　武汉 　430062 ( 2) ) 中国科学院测量与地球物理研究动力大地测量学重点实验室 　武汉 　430077 ( 3) ) 武汉理工大学管理学院 　武汉 　430070 摘要 :对球近似的大地测量边值问题 ,从扰动位的 Poi sson 积分公式出发 ,推导了大地测量 Stokes 边值问题的重力异常球谐函数级数表达式 ,再基于虚功原理 ,将重力异常球谐函数级数表达式写 为变分方程 ,即大地测量边值问题的自然积分方程 ,然后再用离散化技术 ,将所得到的自然积分方 程转化为代数线性方程组 ,通过对代数线性方程组求解 , 即可得大地测量 Stokes 边值问题的数值 解. 关键词 :边值问题 ; 自然边界元解法 ;数值解 中图法分类号 : TP242 DOI :10 . 3963/ j . i ssn .. 20 10 . 02 . 039 0 　引　　言 假定通过重力归算, 已经把大地水准面外的 地球质量去掉, 并假定正常椭球中心与地球质心 大地测量边值问题 ,一般利用 Sto ke s 公式 , 重合, 旋转角速度等于地球自转角速度, 则大地水 ( ) 准面外空间的扰动位满足 L ap lace 微分方程 , 即 采用有限差分法或快速富氏变换 F F T 等来计 Δ 算. 唐元义等研究了用 自然边界元法对大地测量 t = 0 , 而且在无穷远处是正则的. 的第二 、第三边值问题进行求解[ 12 ] ,所采用的方 已知观测量是重力异常 Δg , 求扰动位的问 题, 就是大地测量 Sto ke s 边值问题 法是由N eumann 外问题得大地测量边值问题的 ΔT = 0 在边界面外部 自然积分方程 ,然后再离散求解.