一种局部海域扰动重力数据反...底地形的非线性序列相关方法_翟振和.pdf
《一种局部海域扰动重力数据反...底地形的非线性序列相关方法_翟振和.pdf》由会员分享,可在线阅读,更多相关《一种局部海域扰动重力数据反...底地形的非线性序列相关方法_翟振和.pdf(10页珍藏版)》请在咨信网上搜索。
1、2023 年第 38 卷 第2期2023,38(2):0631-0640地球物理学进展Progress in Geophysicshttp:/wwwprogeophyscnISSN 1004-2903CN 11-2982/P翟振和,孙中苗,马健,等2023 一种局部海域扰动重力数据反演海底地形的非线性序列相关方法 地球物理学进展,38(2):0631-0640,doi:10 6038/pg2023GG0240ZHAI ZhenHe,SUN ZhongMiao,MA Jian,et al 2023 Nonlinear sequence correlation method for inversi
2、on of seabed topography from disturbed gravitydata in local sea area Progress in Geophysics(in Chinese),38(2):0631-0640,doi:106038/pg2023GG0240一种局部海域扰动重力数据反演海底地形的非线性序列相关方法Nonlinear sequence correlation method for inversion of seabed topography fromdisturbed gravity data in local sea area翟振和1,孙中苗1,马健
3、1,李楠2,管斌1,欧阳明达1ZHAI ZhenHe1,SUN ZhongMiao1,MA Jian1,LI Nan2,GUAN Bin1,OUYANG MingDa1收稿日期2022-08-09;修回日期2022-11-10投稿网址http:/www progeophys cn基金项目国家自然科学基金项目(41674082)资助第一作者简介翟振和,男,1980 年生,副研究员,主要从事卫星测高、高程基准方向研究 E-mail:wuji41163 com1 地理信息工程国家重点实验室,西安7100542 宇航动力学国家重点实验室,西安7100431 State Key Laboratory o
4、f Geo-Information Engineering,Xian 710054,China2 State Key Laboratory of Astronautics Dynamics,Xian 710043,China摘要在不考虑地壳均衡条件下,从重力场基本理论出发推导了重力异常/扰动与海深信息之间的近似解析关系,其中扰动重力反演海底地形的模型较之重力异常更加简洁 通过对反演模型级数展开式的深入分析表明,在重力数据分辨率数值小于海深分辨率时并不明显收敛,因此在用重力数据反演高分辨率海底地形时应顾及展开级数的高阶项并采取相应处理方法 在大量统计分析基础上,构建了局部海域扰动重力数据反演高分
5、辨率海底地形的非线性序列相关方法,该方法利用少量船测重力/水深测线数据获取非线性相关参数,而后利用卫星测高重力数据按非线性相关函数模型反演海底地形 试验分析表明,2 次非线性相关函数反演效果表现较优,对于海底地形平缓海域,1 1区域内单条船载测线数据获得的相关参数即可达到 2%相对精度 对于海底地形复杂海域,论文方法适用范围应该尽量缩小,30 30区域范围基于单条船载测线数据反演的相对精度可优于 7%从论文试验结果分析,海底地形与重力数据之间存在 2 次非线性相关函数特征,且这种特征在高分辨率情况下仍然适用,利用非线性函数模型并结合稀疏船载测量数据、密集卫星测高重力数据可为局部海域高分辨率海底
6、地形反演提供一种可行的解算方法AbstractBased on the basic theory of gravity field,theapproximate analytic relationship between gravity anomaly/disturbance and bathymetric information is deduced withoutconsidering the isostasy of crust By analyzing the seriesexpansion of the inversion model,it is shown that theresol
7、ution of gravity data does not converge when theresolution of gravity data is less than that of bathymetricdata,therefore,the high-order terms of the expansionseriesshouldbeconsideredandthecorrespondingprocessing methods should be adopted in the inversion ofhigh-resolution seabed topography from gra
8、vity data Onthe basis of a large number of statistical analysis,anonlinear sequence correlation method for retrieving high-resolution seabed topography from disturbed gravity data inlocal sea area is proposed,which can be used to inversethe seabed topography by using satellite altimetry gravitydataa
9、ndsmallship-bornesurveyingdataTheexperimental resultsshowthatthe 2-ordernonlinearcorrelation functionhasbetterperformance,andtherelative accuracy can reach 2%for the sea area with flatbottom topography based on the parameters obtained fromthe single ship-borne line survey data in the area of 1 1 For
10、 the sea area with complex seabed topography,thescope of application of the method in this paper should bereduced as far as possible,the relative accuracy of 30 30 area inversion based on single shipborne survey linedata can be better than 7%From the analysis of the testresults in this paper,there i
11、s a quadratic nonlinearcorrelation function feature between seabed terrain andgravity data,and this feature is still applicable in the case地球物理学进展www progeophys cn2023,38(2)关键词海底地形;卫星测高;船载重力/水深测量;扰动重力;2 次非线性相关函数中图分类号P223文献标识码Adoi:10 6038/pg2023GG0240of high resolution Using the nonlinear function mo
12、del,combined with sparse shipborne measurement data anddense satellite altimetry gravity data,a feasible solutionmethod for high-resolution seabed terrain inversion in localsea areas can be providedKeywordsSeabed topography;Satellite altimetry;Ship-borne gravimetry/bathymetry;Disturbed gravity;2 ord
13、ernonlinear function model0引言全球宽阔海域海底地形的确定目前主要采用船载多波束测量、卫星雷达测高反演等技术(许厚泽等,1999;李 建 成 等,2003;Becker et al,2009;Sandwell et al,2014;文汉江等,2017),狭窄水道、近海及滩涂区域海深一般采用船载多波束、机载激光测深、多光谱反演等技术确定(赵建虎和刘经南,2008;张熠星等,2018;曹斌等,2016)船载多波束测量和重力测量同时进行则可获得海洋区域某一测线的高精度重力/水深数据,但受限于测量平台和海洋客观环境条件,船载测量技术很难达到全球均匀覆盖 海洋卫星测高通过双星串
14、飞(鲍李峰和许厚泽,2014;翟振和等,2018)、合成孔径高度计等手段可在2.5 年左右时间获取全球高达 1分辨率的重力数据(翟振和和史灵卫,2017),进而可进一步提升海底地形反演精细度 利用最新的 CryoSat-2 和 Jason-1卫星测量数据获得的全球海洋重力场模型在墨西哥湾发现了一条已灭绝的扩张海脊,在南大西洋发现了一条扩张裂谷,在缓慢扩张的海脊上发现了深海丘陵结构以及数千座以前未知的海山(Sandwell etal,2014)目前,各国研究机构广泛应用包括船载测深、卫星测高等数据构建海底地形模型,美国国家环境信息中心(National Centers for Environme
15、ntalInformation,NCEI)、丹 麦 科 技 大 学(TechnicalUniversity of Denmark,DTU)、斯克里普斯海洋研究所(Scripps Institution of Oceanography,SIO)等多个研究机构相继发布了 ETOPO-1、DTU18、SIO 32.1 等全球海底地形格网模型(Amante and Eakins,2009;Smith and Sandwell,1997),模型格网大小可达到1 国内学者结合测高重力异常模型和全球船测水深数据,构建了新版的全球海底地形模型 BAT_WHU2020(胡敏章等,2020)在海底地形反演方法上,
16、利用卫星测高重力数据可通过解析法、频域法、最小二乘配置法、梯度法等反演宽阔海域海底地形(Parker,1973;王勇等,2001;黄谟涛等,2002;罗佳等,2002;聂琳娟等,2012;李丽丽和马国庆,2014;Fan et al,2020),但反演精度较之船载测量结果偏低,同时上述方法受限于扰动重力场与海深的复杂函数关系,在实际解算中只能在区域假设先验信息条件下进行,因此理论上只适合于局部海域计算 重力地质法其原理是采用船测海深数据进行控制,对区域内其他重力数据进行长波改正而后根据线性模型进行反算,由于引入了船测数据控制是目前反演精度较高的方法(Hsiao et al,2011;Kim e
17、t al,2011;胡敏章等,2012;欧阳明达等,2014)在目前构建全球海底地形模型的方法中,主要测高卫星的海面高真实分辨率达不到 1,再加上反演过程通常采用 10 20 km 波长进行低通滤波,因此海底地形模型的真实分辨率往往达不到 1(Smith and Sandwell,1997;胡敏章等,2020),在缺乏船测数据的区域表现尤为明显 通过以上分析可以看出,在全球宽阔海域水深确定中我们既不能完全依靠卫星数据,也不能仅使用高质量的船测数据,较好的思路是将船载重力/水深数据和卫星测高反演获得的重力数据进行高效融合,以提高海底地形反演的准确性和效率,达到数据资源的优化使用 基于以上考虑本文
18、深入分析了重力扰动与海深之间的解析函数关系,从重力/地形相关统计特性出发提出了一种适合局部海域海底地形反演的非线性相关方法,并在南海海域开展了计算试验1重力扰动/异常与海深的解析函数关系分析由地球重力场理论可知,地球质体对空间一点P 的扰动引力位 Tp(xp,yp,zp)在空间直角坐标系下可表示为(Hofmann-Wellenhof and Moritz,2005):Tp(xp,yp,zp)=VP VE=G(D(x,y,z)DE(x,y,z)(x xP)2+(y yP)2+(z zP)2dxdydz,(1)2362023,38(2)翟振和,等:一种局部海域扰动重力数据反演海底地形的非线性序列相
19、关方法(www progeophys cn)式中,Vp表示实际地球对 P 点引力位,VE表示水准椭球对 P 点的引力位,G 表示引力常数,D(x,y,z)表示实际地球单位质体的密度,DE(x,y,z)表示正常椭球单位质体的密度忽略径向与垂直方向差异并假设海面一点重力异常/扰动由实际地球海水部分与正常椭球对应的海洋空间区域质量差异引起,则可得海面一点的重力异常 gp:gP=GDh(0z hP(x xP)2+(y yP)2+(z hP)2)3/)2dxdydz 2GDh(01(x xP)2+(y yP)2+(z hP)2dxdydz,(2)式中 D 表示海水密度与 DE(x,y,z)之差,h 表示
20、单位质体的海深,hp表示 P 点的海深,进一步积分得:gP=GD1(x xp)2+(y yP)2+(h hP)2)1/21(x xP)2+(y yP)2+h2P)1/2dxdy 2GD ln(x xP)2+(y yP)2+(h hP)2)1/2+h hP)ln(hP+(x xP)2+(y yP)2+(hP)2)1/2)dxdy(3)按照上述推导思路,也可以得到扰动重力与海深之间的近似解析关系:gp=GD1(x xp)2+(y yp)2+(h hp)2)1/21(x xp)2+(y yp)2+h2p)1/2dxdy(4)对比式(3)与式(4)可以看出,对于重力场反演海底地形的工程实际而言,扰动重
21、力较重力异常而言更加简洁 由式(4)可以看出,海洋表面一点的重力扰动由整个海洋区域的海深、计算点位置及密度差决定 实际地球内部由地壳、地幔及地核组成,而水准椭球的质量分布较为复杂,仅有少量文献(陆仲连,1996)对水准椭球密度进行了假设分析 然而,即便密度差 D 是常数,式(4)的求解依然很复杂,其本质上可归为一个大型多元非线性方程组的求解 由以上推导可以看出,在多个近似假设下,利用重力扰动/异常求解海深信息依然是一个非常复杂的问题为了便于进一步简化求解,可将式(4)关于海深 h 进行麦克劳林级数展开,设 L2=(x xp)2+(y yp)2,则:gp=GD1(L2+(h hp)2)1/21(
22、L2+h2p)dxdy(5)式(5)关于海深 h 的零阶项为 0,1 阶项为:g1p=GDhph(L2+(hp)2)3/2dxdy,(6)2 阶项为:g2p=12GDh2(L2+(hp)2)3/2dxdy+32GD(hp)2h2(L2+(hp)2)5/2dxdy,(7)3 阶项为:g3p=32GD hph3(L2+(hp)2)5/2dxdy 54GD h3ph3(L2+(hp)2)7/2dxdy,(8)4 阶项展开并忽略h4L5以上项得到:g4p=38GDh4(L2+(hp)2)5/2dxdy 154GD(hp)2h4(L2+(hp)2)7/2dxdy+(9)最终得到 gp展开至 4 阶项且忽
23、略h4L5以上项的表达式为:gp=GDhph(L2+(hp)2)3/2dxdy 12GDh2(L2+(hp)2)3/2dxdy+32GD(hp)2h2(L2+(hp)2)5/2dxdy 32GDhph3(L2+(hp)2)5/2dxdy+38GDh4(L2+(hp)2)5/2dxdy,(10)式(10)虽然包含了 4 阶项但并不能明显判断其收敛性,文献(Zhu and Jekeli,2009)对其类似表达式收336地球物理学进展www progeophys cn2023,38(2)敛性的处理进行了详细讨论,这里不再赘述 本文重点考虑 1 4 阶项对反演的影响,分别以 ETOPO-1模型 1、5
24、分辨率数据计算扰动重力为例来进行讨论,试验区位于 15N16N,115E117E,将 1 4阶分别单独计算获得扰动重力,统计所获得数值绝对值的平均值,以此衡量各阶对扰动重力的贡献,总的量级利用 1 4 阶求和得到,具体结果见表 1、表 2表 11 4 阶海底地形(1分辨率)在计算扰动重力中所占量级(单位:105m/s2)Table 1The magnitude of 1 4 order seafloor topography in calculatingdisturbing gravity(1 resolution)(unit:105m/s2)计算区域1 阶项量级2 阶项量级3 阶项量级4 阶
25、项量级总的量级15N16N,115E116E43.413.2714.836.2538.1115N16N,116E117E107.549.0636.9815.8193.94表 21 4 阶海底地形(5分辨率)在计算扰动重力中所占量级(单位:105m/s2)Table 2The magnitude of 1 4 order seafloor topography in calculatingdisturbing gravity(5 resolution)(unit:105m/s2)计算区域1 阶项量级2 阶项量级3 阶项量级4 阶项量级总的量级15N16N,115E116E60.2831.2222
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 一种 局部 海域 扰动 重力 数据 地形 非线性 序列 相关 方法
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。