利用扰动重力数据反演海底地形的高斯曲面函数估计方法

doi:10.11947/j.AGCS.2024.20230204

摘要/Abstract

摘要： 利用2017年1月—2020年12月的SARAL卫星测高数据反演得到南海海域1'×1'格网扰动重力数据，通过与船载重力数据比较，精度达到5.5 mGal。提出了利用重力数据及高斯曲面函数求解区域特征参数进而估计海底地形模型的方法，利用SARAL卫星反演获得的重力数据在南海海域开展了计算试验。结果表明，利用ETOPO-1先验模型，在10×10格网一组划分条件下估计得到的1'×1'海底地形精度相对于先验模型提高了约10 m；利用DTU18先验模型，在9×9格网一组划分条件下估计得到的1'×1'海底地形精度相对于先验模型提高了约9 m。计算结果表明，利用卫星测高反演得到的重力数据，通过高斯曲面函数解算获得的5个特征参数，可以在一定程度上代表相应区域海底地形的曲面特征，进而可在不依赖外部实测数据条件下对先验海底模型进行精化。对于曲面估计而言，格网划分越小，曲面函数就越能反映区域变化特征。因此，对于未来卫星测高技术发展而言，更高分辨率的重力场探测技术有望继续提升海底地形细节的反演能力。

关键词: 扰动重力, 海底地形, SARAL卫星, 高斯曲面估计, ETOPO-1, DTU18

Abstract: In this paper, 1'×1' grid disturbing gravity data in the South China Sea are derived from the satellite altimetry data of SARAL from January 2017 to December 2020. By comparing with the ship-borne gravity data, the accuracy is 5.5 mGal. A method is proposed to estimate the terrain model of the seabed by using gravity data and the regional characteristic parameters solved by Gauss surface function. For ETOPO-1 prior model, the accuracy of 1'×1' grid seafloor terrain estimated under a group of 10×10 grids is improved by about 10 meters. For DTU18 prior model, compared with the prior model, the 1'×1' grid seafloor terrain accuracy estimated under a set of 9×9 grids is improved by about 9 meters. To a certain extent, the results show that the five characteristic parameters obtained from gravity data retrieved from satellite altimetry can represent the surface features of the seabed topography in the corresponding area by solving the Gauss surface function, furthermore, the prior seafloor model can be refined by iterative cycle without relying on ship-borne data. Theoretically, for surface estimation, the smaller the mesh, the better the surface function can reflect the regional variation. Therefore, for the future development of satellite altimetry technology, it is expected that the gravity field detection technology with higher resolution will continue to improve the retrieval capability of seabed topographic details.

Key words: disturbing gravity, seabed topography, SARAL satellite, Gaussian surface estimate, ETOPO-1, DTU18

中图分类号:

P258

翟振和, 孙中苗, 管斌, 马健, 李端. 利用扰动重力数据反演海底地形的高斯曲面函数估计方法[J]. 测绘学报, 2024, 53(2): 231-238.

ZHAI Zhenhe, SUN Zhongmiao, GUAN Bin, MA Jian, LI Duan. An estimation method of seabed topography based on Gauss surface function using ocean gravity data[J]. Acta Geodaetica et Cartographica Sinica, 2024, 53(2): 231-238.

参考文献

[1] 许厚泽, 王海瑛, 陆洋, 等. 利用卫星测高数据推求中国近海及邻域大地水准面起伏和重力异常研究[J]. 地球物理学报, 1999, 42(4):465-471. XU Houze, WANG Haiying, LU Yang, et al. Geoid undulations and gravity anomalies from t/p and ERS-1 altimeter data in the China Sea and vicinity[J]. Chinese Journal of Geophysics, 1999, 42(4):465-471.
[2] 李建成,陈俊勇,宁津生,等.地球重力场逼近理论与中国2000似大地水准面的确定[M].武汉:武汉大学出版社, 2003. LI Jiancheng, CHEN Junyong, NING Jinshen, et al. Approximation theory of Earth gravity field and determination of China 2000 quasi geoid [M]. Wuhan:Wuhan University Press, 2003.
[3] SANDWELL D T, MVLLER R D, SMITH W H F, et al. New global marine gravity model from CryoSat-2 and Jason-1 reveals buried tectonic structure[J]. Science, 2014, 346(6205):65-67.
[4] 赵建虎,刘经南.多波束测深及图像数据处理[M].武汉:武汉大学出版社, 2008. ZHAO Jianhu, LIU Jingnan. Multi beam sounding and image data processing [M]. Wuhan:Wuhan University Press, 2008.
[5] 张熠星,尚建华,贺岩.机载激光测深技术的研究进展[J].激光技术,2018,42(5):588-592. ZHANG Yixing, SHANG Jianhua,HE Yan. Research progress of airborne laser sounding technology[J].Laser Technology, 2018,42(5):588-592.
[6] 曹斌, 邱振戈, 曹彬才. 四种遥感浅海水深反演算法的比较[J]. 测绘科学技术学报, 2016, 33(4):388-393. CAO Bin, QIU Zhenge, CAO Bincai. Comparison among four inverse algorithms of water depth[J]. Journal of Geomatics Science and Technology, 2016, 33(4):388-393.
[7] 文汉江, 金涛勇, 朱广彬.卫星测高原理及应用[M].北京:测绘出版社,2017. WEN Hanjiang, JIN Taoyong, ZHU Guangbin. Principle and application of satellite altimetry(in Chinese) [M]. Beijing:Surveying and Mapping Press, 2017.
[8] 鲍李峰, 许厚泽. 双星伴飞卫星测高模式及其轨道设计[J]. 测绘学报, 2014, 43(7):661-667. BAO Lifeng, XU Houze. Twin-satellites altimetry mode and its orbit design[J]. Acta Geodaetica et Cartographica Sinica, 2014, 43(7):661-667.
[9] 翟振和, 孙中苗, 肖云, 等. 自主海洋测高卫星串飞模式的设计与重力场反演精度分析[J]. 武汉大学学报(信息科学版), 2018, 43(7):1030-1035, 1128. ZHAI Zhenhe, SUN Zhongmiao, XIAO Yun, et al. Two-satellites tandem mode design and accuracy analysis of gravity field inversion for independent marine altimetry satellite[J]. Geomatics and Information Science of Wuhan University, 2018, 43(7):1030-1035, 1128.
[10] AMANTE C, EAKINS B W. ETOPO1 global relief model converted to PanMap layer format[R]. [S.l.]:LSLF NOAA-National Geophysical Data Center, 2009.
[11] BECKER J J, SANDWELL D T, SMITH W H F, et al. Global bathymetry and elevation data at 30 arc seconds resolution:SRTM30_PLUS[J]. Marine Geodesy, 2009, 32(4):355-371.
[12] 胡敏章, 张胜军, 金涛勇, 等. 新一代全球海底地形模型BATWHU2020[J]. 测绘学报, 2020, 49(8):939-954. DOI:10.11947/j.AGCS.2020.20190526. HU Minzhang, ZHANG Shengjun, JIN Taoyong, et al. A new generation of global bathymetry model BATWHU2020[J]. Acta Geodaetica et Cartographica Sinica, 2020, 49(8):939-954. DOI:10.11947/j.AGCS.2020.20190526.
[13] PARKER R L. The rapid calculation of potential anomalies[J]. Geophysical Journal International, 1973, 31(4):447-455.
[14] 王勇, 许厚泽, 詹金刚. 中国海及其邻近海域高分辨率海底地形[J]. 科学通报, 2001, 46(11):956-960. WANG Yong, XU Houze, ZHAN Jingang. High resolution seabed topography in the Sea of China and its adjacent waters[J]. Chinese Science Bulletin, 2001, 46(11):956-960.
[15] 黄谟涛, 翟国君, 欧阳永忠, 等. 卫星测高资料在反演海底地形中的应用[J]. 海洋测绘, 2002, 22(1):3-7. HUANG Motao, ZHAI Guojun, OUYANG Yongzhong, et al. Application of satellite altimetry data in inversion of seabed topography[J]. Hydrographic Surveying and Charting, 2002, 22(1):3-7.
[16] 罗佳, 李建成, 姜卫平. 利用卫星资料研究中国南海海底地形[J]. 武汉大学学报(信息科学版), 2002, 27(3):256-260. LUO Jia, LI Jiancheng, JIANG Weiping. Bathymetry prediction of South China Sea from satellite data[J]. Editoral Board of Geomatics and Information Science of Wuhan University, 2002, 27(3):256-260.
[17] HSIAO Y S, KIM J W, KIM K B, et al. Bathymetry estimation using the gravity-geologic method:an investigation of density contrast predicted by the downward continuation method[J]. Terrestrial, Atmospheric and Oceanic Sciences, 2011, 22(3):347.
[18] KIM J W, VON FRESE R R B, LEE B Y, et al. Altimetry-derived gravity predictions of bathymetry by the gravity-geologic method[J]. Pure and Applied Geophysics, 2011, 168(5):815-826.
[19] SMITH W H F, SANDWELL D T. Global sea floor topography from satellite altimetry and ship depth soundings[J]. Science, 1997, 277(5334):1956-1962.
[20] 聂琳娟, 吴云孙, 金涛勇, 等. 基于海水质量亏损引起的重力异常反演南海海底地形[J]. 大地测量与地球动力学, 2012, 32(1):43-46. NIE Linjuan, WU Yunsun, JIN Taoyong, et al. Inversion of submarine topography of South China Sea by using gravity anomaly caused by mass deficiency[J]. Journal of Geodesy and Geodynamics, 2012, 32(1):43-46.
[21] 胡敏章, 李建成, 金涛勇. 应用重力地质方法反演皇帝海山的海底地形[J]. 武汉大学学报(信息科学版), 2012, 37(5):610-612, 629. HU Minzhang, LI Jiancheng, JIN Taoyong. Bathymetry inversion with gravity-geologic method in emperor seamount[J]. Geomatics and Information Science of Wuhan University, 2012, 37(5):610-612, 629.
[22] 欧阳明达, 孙中苗, 翟振和. 基于重力地质法的南中国海海底地形反演[J]. 地球物理学报, 2014, 57(9):2756-2765. OUYANG Mingda, SUN Zhongmiao, ZHAI Zhenhe. Predicting bathymetry in South China Sea using the gravity-geologic method[J]. Chinese Journal of Geophysics, 2014, 57(9):2756-2765.
[23] 范雕, 李姗姗, 孟书宇, 等. 联合多源重力数据反演菲律宾海域海底地形[J]. 测绘学报, 2018, 47(10):1307-1315. DOI:10.11947/j.AGCS.2018.20170423. FAN Diao, LI Shanshan, MENG Shuyu, et al. Recovery of bathymetry over philippine sea by combination of multi-source gravity data[J]. Acta Geodaetica et Cartographica Sinica, 2018, 47(10):1307-1315. DOI:10.11947/j.AGCS.2018.20170423.
[24] FAN Diao, LI Shanshan, MENG Shuyu, et al. Applying iterative method to solving high-order terms of seafloor topography[J]. Marine Geodesy, 2020, 43(1):63-85.
[25] XU Huan, YU Jinhai, WAN Xiaoyun, et al. An expression for gravity generated by an anomalous geological body and its application in bathymetry inversion[J]. Journal of Geodesy and Geoinformation Science, 2021, 4(4):63-73.
[26] 翟振和, 孙中苗, 马健, 等. 一种局部海域扰动重力数据反演海底地形的非线性序列相关方法[J]. 地球物理学进展, 2023, 38(2):631-640. ZHAI Zhenhe, SUN Zhongmiao, MA Jian, et al. Nonlinear sequence correlation method for inversion of seabed topography from disturbed gravity data in local sea area[J]. Progress in Geophysics, 2023, 38(2):631-640.
[27] 翟振和, 孙中苗, 王兴涛. 全球及局部海洋扰动重力反演的快速解析方法[J]. 测绘学报, 2015, 44(8):827-832. DOI:10.11947/j.AGCS.2015.20140482. ZHAI Zhenhe, SUN Zhongmiao, WANG Xingtao. The analytical and quick computation method of disturbing gravity in global and local ocean area[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(8):827-832. DOI:10.11947/j.AGCS.2015.20140482.