物理大地测量边值问题的自然边界元解法.pdfVIP

维普资讯 第 32卷 第 1期 武汉理工大学学报(鸯 差) Vo1．32 No．1 2008年 2月 JournalofWuhanUniversityofTechnology Feb．2008 (TransportationScienceEngineering) 物理大地测量边值问题的自然边界元解法* 唐元义 许厚泽 (湖北大学数学与计算机科学学院” 武汉 430062) (中国科学院测量与地球物理研究所” 武汉 430077) (中国科学院研究生院” 北京 100049) 摘要 ：物理大地测量边值问题是物理大地测量学的核心问题和理论基础．在数学上表示为Laplace 方程的以大地水准面为边界面的第三外部边值 问题 ，即Stokes边值 问题 ．文 中在物理大地测量 Stokes边值 问题 的球近似下，将 Stokes边值问题转化为 Neumann外问题 ，用 自然边界元法 ，对 Neumann外问题进行 自然边界归化 ，得到 自然边界积分方程 ，再求 自然边界积分方程 的数值解 ，并 给 出了一个算例． 关键词 ：Stokes边值问题 ；自然边界元法 ；数值解 中图法分类号 ：TP242 在大地测量中，地球重力场的研究在理论上 函数．大地测量边值问题 ，就是利用观测数据给定 归结为解算大地测量边值 问题．在数学上表示为 边值条件，来推求扰动位及外部重力场的问题．当 Laplace方程的 以大地水准面为边界面的第三外 已知观测量是重力异常Ag时，求扰动位的问题归 部边值问题，即Stokes边值问题 ，它要求将地面 结为解下列微分方程 的边值问题． 上的实测重力值归算到大地水准面上．由于大地 f △丁一 0 在边界面外部 水准面是未知的，因此实际上是以同大地水准面 Ja 一专ya一一Ag 在边界。面-上L‘ 十分接近 的旋转椭球面作为边界面．边界面为旋 转椭球面的边值问题 ，解有解析的表达式 ，但作 了 【T一0(专) y—oo 某种情况 的近似．在实际求解物理大地测量边值 第 l式为扰动位丁的Laplace方程；第2式为物理 问题时 ，通常是用球面近似代替旋转椭球面 ，即球 大地测量学的基本微分关系 ，式 中的y为正常重 近似，然后再作椭球改正． 3 11 力 ，左端是丁和 的线性组合 (为边界面的外 若需要求解得到重力数据值，事实上都是作 法线方向)，右端的重力异常是已知的． 数值计算．常用的方法有有限差分法、有限元法[1 物理大地测量 Stokes理论是重力场逼近中 等．本文拟对于球近似的物理大地测量Stokes边 的经典理论．根据这种理论，大地水准面外不能有 值问题用 自然边界元法进行数值求解． 物质存在 ，所以在归算时，必须把大地水准面外的 物质去掉或移到大地水准面内部去．重力归算过 1 Stokes边值 问题 程不可避免地牵涉到地壳的内部构造——密度分