您的当前位置:首页正文

大地水准面与似大地水准面之差及其模型检核

2024-05-30 来源:易榕旅网
维普资讯 http://www.cqvip.com

第29卷第l期 2002年2月 J0URNAL OF CHENGDU LR'qIVERSITY 0F TECHNOLOGY 成都理工学院学报 .29 No l Feb.2002 [文章编号]1005—9539(2002)0[0105 05 大地水准面与似大地水准面之差及其模型检核 张赤军边少锋 (中国科学院测地所动力大地测量学开放研究实验室,武汝430077) [摘要]大地测量学曲发展,要求能岱精确推得大地水准面高(Ⅳ),为此,在阐述顾及扰动重力 垂直梯度的大地水准面-q似大地水准面(f)之差曲公式基础上,对它们在三种模型上进行了数 值甘算却槛校。结果表 ,谤套式曲精度是相当高的.由此可以将正常高化算为正高(海拔高) 关建词]大地水准面i位大地水准面;正高;正常高i模型 分类号]P22 [文献标识码:A 中国的郭守敬早在1276年的水利测量中就 注意到以海平面作为高程基准面的问题_1j 到了 l8世纪,在法国的大地测量工作中,这一基准面 或辅助面也得到了应用;并于1819年L.Puissant 大地水准面(Ⅳ)与似大地水准面(})之问又 存在着一定的内在联系,例如它们都是位理论第 三边值问题的解,而Stokes解就是Molodensky 解的零次近fc上。为了进一步揭示它们之间的差异 与联系.近年来L E S5oberg(1995)¨ 和边少锋等 在自己写的大地测量教程中,将海潮最高水位和 最低水位的平均值称为绝对海平面,他还指出有 (1 999)口一先后采用不同方法,推导了顾及高程二 必要将大地网投影到向大陆延伸的静止的海平面 上 在德国,继K.F.Gauss和Bessel的大地测量 工作之后,L.Lichter ̄stelP,于1872年首次将上述 基准面称之为大地水准面口]。实际上在此之前的 1849年、英国的G G.Stokes发表了研究地球形 次项的(Ⅳ一})公式,这比文献Ea] ̄Eao]提供的 近似公式要精确得多;在文献E6]和ET]中.前者的 二次项用的地面点的混合重力异常的垂直梯度. 后者用的是纯重力异常(扰动重力)的垂直梯度, 两者本质上是一致的.数值上也相差甚微 他们的 状的论文,该形状就是大地水准面的形状,但真正 使Stokes理论得以实现的还在20世纪。在1945 ~贡献在于:在已知H (正常高)和}的情况下为精 确求取正高(Ho)和Ⅳ提供了理论依据。张赤军 等(1997,1 999)还利用上述公式推求了珠峰及野 猫山的大地水准面和正高一”J-L”j。可以认为,随着 地形测量和地质地球物理勘探的广泛开展,以及 1 960年期间前苏联的M.C.Molodensky避开 了大地水准面外不应存在质量的缺点,提出了研 究似大地水准面与真实地球形状的理论嘲 ̄【 ,并 在此后得到了应用和发展。到目前止,由于卫星测 高、海空重力测量、垒球定位系统(GPs)等新手段 的成功应用,以及卫星重力梯度仪等空间手段即 将实现,使大地测量边值问题的研究在不断向前 发展,也使大地水准面和似大地水准面的研究继 续引向深入。由于大地水准面具有明显的物理意 新技术合成孔径雷达(ISAR)等的应用,顾及高程 二次项的(Ⅳ一})的公式会进一步得到重视和运 用。实际上,R-H、Rapp于】997年已强调深入研 究这一问题的重要性u 为了进一步检核Sjoberg和Bian推导公式的 正确性.本文采用三个模型作了模拟计算,结果也 较理想 义(等位面)和甚为实用的价值(高程基准面),因 此研究和精化大地水准面仍很重要。 :收穑日期 2000-12-26 [基金项目]国家自然科学基金(49874016,49674208)和国家攀登计划资助项目 [作者简介]张赤军(1033一), ,研究黄,地球重力学专业(E mail:Fang@aseh.whigg.ac.cn) 维普资讯 http://www.cqvip.com

・lO6・ 成都理工学院学报 第29卷 1 顾及扰动重力垂直梯度的(Ⅳ一 一 f)公式及其解释 H( 308 ̄+4 G P一赘 ㈤ 又由文献[9]可知, 出文献[6]可知,正高(H )及正常高(H )有 7一 7d7—7。 3806 ㈤ 如下关系: —H0一N—f—Ho( 一7)/ (1) 式中 为参考椭球上的正常重力值.如此有 正常重力是一解析函数,故其平均值 可由 g一7一gp+0.3806 一2 H一 公式很易计算出来 而重力平均值的计算比较复 杂.根据一些岩石密度的实际样品及地形资料有 旦豢一 +。2 拼 .a806旦 (2 。)可能以相当高精度求得,这时大地水准面与似大 又据Bouguer异常的定义: 地水准面的偏差就能很好地求取。 —gp一2ⅡG PH+0.R086H一 由-J: (10j g 一7 +gp一7 —y 十 (2) 则 一7一 一 警 (11) 如此在地面以上空间(忽略大气影响或视之 将(】])式代人(1)式,最后得: 为无质量)有 N一;一 H一 警 (】2j ah ah— l0 hl +盘。 ah l c3) 式中H已略去其上下角标,且 公里为单位。因 在地面以下时,由于地形质量的影响,经过界 为即使正高和正常高相差几米,于计算精度无妨。 面时有一与质量有关的突变 顺便指出,上一公式与Sjoberg导出的公式 一 襄1 孔+警1 + G ㈨ Ⅳ一f一 H一 (13) 式中.(;为牛顿引力常数;p为地壳密度,一般可 稍有不同,两者互差甚微,其误差完垒可以不计i 取为 ̄=2670 kg/m ;+p表示相应导数从空间向 但(1 2)式中右边第二项为同一点上重力异常(或 地面求取;一声表示相应导数从地形内部向地面 扰动重力)的垂直梯度,而(1 3)式为空闻异常 求取 (ags)的垂直梯度。文献[73的推导过程也较为简 设实际重力在地形内部可按Taylor级数展 便。 开,仅取其一阶项,则在距地面为。的地内q点 式(12)和(1 3)中右端的第二项为二阶近似 有: 项,它与异常的梯度直接相关,坷此局部效应更为 一 一 l (H一) 显著。如果说重力垂直梯度刻划了它与大地水准 面的平均曲率的联系,那幺异常垂直梯度则描述 一 +ah +4 ̄G,o (H 了大地水准面平均曲率与正常水准面平均曲率的 (5) 偏离。这样就把大地测量中的几何概念和力学概 如略去J/ah随纬度及高度的变化(这些变 念具体地联系起来。 化在一定范围和有限高程(H)内均比较小), 0有 式中的第一项为高程一次项的系数(△ ),它 岛一 一(一o.3。86+鲁+4 G p)(H一。) 反映了区域场的影响;而式中的第二项即高程二 次项的系数则反映局部场的作用,只有一并考虑 (6) 区域场和局部场的贡献,这时由似大地水准面转 式中(H— )以公里为单位。 化成大地水准面才能在物理上和数学意义上比较 对gq在垂线方向求其积分平均值后,即得 垒面地反映出来。 g,这时, g一南1 ( )dz 2Ⅳ一f公式的模拟计算与检核 一-Lf? 一( + 呻+ (日 ]以 2.1模型I(半球体) 设地球表面为一平坦的表面,并取之参考面 维普资讯 http://www.cqvip.com

第1期 张赤军等:大地水准面与似大地水准面之差及其模型检棱 ・107・ 为AOA (见图1);其上覆盖一半球体的高山(异 常体).球之半径为R,顶点P距球心(坐标原点) 的高度为H,且令H=R=3 km;球体的密度。与 岩石的平均密度相同,即为2670 kg/m。;设参考 面上的重力y为9.8Ills,球顶点P处的引力位 可根据文献[14]推得,且 v(P)一2 G P: 一 R (14) 图2 Fig.2 Sphcro[d 一 a 一常密度( )为3000 kg/m。.半径为R一6 km的球 状异常体,它出露、埋藏于地球表面的上、下方而 没有挤出地球的质量,且视地球表面为一参考面. 而该参考面将异常球体平分。又设地球的平均密 度(02)为常密度且等于6000 kg/m ,半径R为 6000 km,由此算得参考面上的重力y一 10.058 123 ms一,在异常球顶点P处的扰动重力 或重力异常: {丌GP Rl一 一502.9 X 10 ms。” (20) 布格异常: Agb一△ 一2rtGo ̄R1一一号ⅡG R】 一一251.4 X 10一 ms一 (21) 似大地水准面高: }一 一 PI皤一3.00m(22) 大地水准面高(该点在体内): N一鲁丌GP1(碍一hr ̄)/7—4.50m(23) 这时Ⅳ一 一4+50—3.00一1.50m 异常球顶点P处的扰动重力垂直梯度 盘P 一一1676EOh 一 a 一一 3 … 垫OH 一(24)…  号√ =1319E,(1E一10- ̄s-2) 将Agb和斯/劬代人公式(12),得 Ⅳ一}一一1.50+3.00=1.50 m 这一结果与由公式(22)和(23)算得的结果完 全一致,也与文献[2]用水准重力算得的正常高与 正高之差的结果相同 2 3模型Ⅲ(均衡补偿的两个球体) 这里取的模型也与文献[2]相同(见图3) 设 异常体为半径( )分别等于2,0 km互为补偿的两 维普资讯 http://www.cqvip.com

・108・ 成都理工学院学报 第29卷 个均质球体.它们的密度分别为 和 .且 一 此值与由公式(25)和(26)算得之结果0.316 3000 kg/m。, 一一3000 kg/m ,正密度球体在 相比.互差为0.056 m,这也是相当小的。 上,负密度球体在下i它与参考面相切于 处;两 虽然以上三个模型不能完垒代表实际的地 球心(0, )的距离h。一3000 m,其联线与参考面 球.怛从一个侧面说明了顾及高程二次项的Ⅳ一f 垂直。取正常重力7=9.797 20 Ills.显然,从形 的级数展开式是具有较高精度的。 式上看.它们的质量得到完全的均衡补偿,这时由 扰动位【异常位)引起的似大地水准面高: 3正常高(H )与正高( ) 一 7 \ L R1一 十0h )』 一 R ̄---ho 996 m正常高与正高是高程系统中的两种高程,如 前所述,正常高是前苏联学者莫洛金斯基为了研 (25) 究地表面的真实形状而设置的高程,因此它在研 究地球形状中很有意义,同时还可以从已知的位 基数Igdh与正常重力y平均值之比而精确地求 出(当然由于用到水准和重力测量的结果也存在 一定的误差).前苏联、东欧和中国都采用这一高 程 但正常高是由似大地水准面起算的高程,因此 它仅有几何意义,没有具体的物理意义。 正高又称海拔高,它是从大地水准面(与垒球 乎均海面非常按近)到地面点的高程,它不仅具有 几何意义而且具有明显的物理意义.井在经济和 图3 国防上和科学上起重要作用。精确推求正高和确 Fig.3 Two spheroids compen ̄ted 定大地水准面与似大地水准面之间的偏差有直接 关系,因为上述之差就是正常高与正高之差。正是 式中Ag。一Gm/R 一2514×10 n1s~, :为小球 由于正高具有明显的几何和物理意义以及直观实 之质量。 用的原因.欧洲和美洲等许多国家都采用正高系 同理,由文献(2)推得大地水准面高: 统。又由于正高是位基数与地表面到大地水准面 Jv—A__ goi 募+ ro]一7,312m 重力乎均值之比,后者牵涉地球表层内部密度的 分布,而过去人们知之不多,以致许多学者认为无 (26) 这时Ⅳ一f一0.316m。 法精确求得.例如文献[8]指出“这种几何意义很 由于二小球引起上一小球顶点P处的重力 清楚的正高要领严格说是不能实现的”。然而事物 异常(扰动重力) 的发展并非像该文断言的那样,这是因为,一方面 &gp一4现令对许多地区的岩石密度已有较详细的了解. r ̄pro(1一 ) 并已具备了较好的地形图和重力资料,这样就可 一436×10 ms (27) 以相当精确地计算出Ⅳ一f值,并可使正常高化 在 减去层间改正377×10 n1S_”后便得△ , 为正高;另一方面,根据上述资料和相应的梯度公 且△ 一59×10一 rns。 式,可以推算出地面与大地水准面之间的重力乎 而重力异常垂直梯度: 均值 ,比Helmert方法更为精确,从而可由位基 亟“密度的误差对 Oh一一号 ( 一 j 数精确地求取正高。文献[9]指出,正高的影响很小”,这是对的;根据我们的估计,如 一(一0.1676+0.1259)×10 5s 密度的误差为士100 kg/m。,对于4000 m的高山 一一417E(28) 而言,由此影响布格异常将为17×10 rns,垂 将上述结果代入公式(12),可得 直梯度的影响将为10E,它们对Ⅳ一f的总影响 N—f一0.372m 将仅为0.07 m,由于 维普资讯 http://www.cqvip.com

第1期 月 —H 张赤军等:大地水准面与似大地水准面之差爰其模型橙棱 (Ⅳ一}) (29) ・109・ [3:Mo[odenskey M S,Eremeev V F.Yourkinna M J Methods for Study of the External Gravitational Fie Ld 对于中国西部4000 m左右的高山地区,当 H 的误差为厘米级时,Ho的精度稍逊于H.,但 它的误差也有可能达到厘米级或稍大一点。 and Figure of the Earth__M]Jeruslaem:Israel Pn gren for Scienlific Transations,i962. [4]管泽霖,宁津生.地球形状及外部重力场[M].北京: 删绘出版社,1981. 4结语与讨论 本文的工作在于:进一步指出了大地水准面 :5]许厚泽,陆仲连.中国地球重力场与大地水准面[M]. 北京:解放军出版社.1997. 与似大地水准面在物理和数学意义上的内在联 系.并且从理论上解释了它们之间的差值完垒可 6]Sjoberg I .On the 9ums[geoid to geoid separationEJ]. Manuscripta Geodetica,1995,20:i82—1 92. 以由近似值向精确化方向转变。为实现这一目标. 一[7]边少锋,张赤军.论大地水准面与似大地水准的偏差 [A].庆祝陈永龄院士9o寿辰文集[c:北京:测绘出 版社,I999. 方面可以利用实测的重力及其梯度资料,另一 方面也可以从地形及地球表层的岩石密度资料中 求取,这在科技日益发展的今天是容易做到的 需 说明的是,这里的工作性质和Stokes问题中的重 力归算不一样,后者是为了求取大地水准面上的 边值,因此这里不存在什么肯定Stokes问题和否 定IVlolodensky的问题。 [83 Pddin R P Advanced Geode-sy[M].Mosoow:Min- era]Press.1978 :9]Heiskanen W A,Moritz H.Physical Geode.,y_-M- San Frand ̄o:Freeman and Company,1967. El0:,Torge W Coeodesy[M].Ber[in:Waker de Grayter &Cu,1980. 文中列举的三个模型,仅对所推导的公式作 些初步计算和检核,但这些模型不一定符合地球 的客观实际,因此,对于更详细的计算和更多的模 型选取等问题仍需进一步研究。 奉项研究得到国家基金委(NO:49874016, Ei]]张赤军.珠穆郎玛峰大地水准面与高程的l删定[J]. 科学通报.1S97,42(23):2543—2545. [12:张赤军.野猫山的大地水准面与正高D:地壳形变 与地震,2000.20(3):lO2—103. E13]Rapp R H use of potentila coefficient models for geold undulation determinations using a spherieal harmonic representation of the height anoma y/geoid .No:49674208)和国家攀登计划中的“中国地区重 力场与大地水准面研究”课题的资助,特在此致 谢 undu[ation diference[J].Joura[of Geodesy.1997, 7l(5):274—282. [14] ̄rang z T Geold and Crustla Structure in Fen 参198l_ 考文献 noo ̄andia口].Kirkkonunm ̄i,Pub[ication of the Finnish Geodetic Ir ̄stitute.1998.126. [1 删绘辞典编委会.测绘词典[M].上海:辞书出版社 ElS]壬谦身,张赤军,蒋福珍,等.微重力测量——理论、 方法与应用[M]・北京:科学出版社,1995 [2 2 Yeremeeev B,Yourkina M N.Height T} eory Df Earth Gravity[M].Mcemow:Minerla Press,1972. DIFFERENCE BErWEEN GEOlD AND QUASI-GE0ID AND rI Ⅳ【oDEL VERIFICATION ZHANG Chi— un.BIAN Shao-feng (I ̄,orotory of Geodynarm'cal C-eodexy,Institute o/Geodesy ate4 Geophysics,Chi ̄ese Academy of Sciences,Wukan) Abstract:The development of geodesy calls for determining precisely height of geoid(Ⅳ).If the difference between geoid and quasl-geoid(})call be expressd by a precNce formula and is known,Ⅳcan be determined.Three models ale used and the(Ⅳorthometric height. Key wnrcls:geoid quasi—geoid;orthometric height;normal height 。})calculated in this paper.The numerical re.suits indicate that the accuracy of formula is very high,and the normal height can be reduced precisely to 

因篇幅问题不能全部显示,请点此查看更多更全内容