借助EGM2008的海域重力异常反距离加权插值算法

构建高精度高分辨率的海域重力异常模型,是将海洋重力测量数据应用于全球高程基准统一和水下重力匹配导航等技术领域的关键步骤之一。针对反距离加权插值算法仅利用了海域重力观测值的空间信息,没有考虑重力观测值物理特性的缺陷,提出一种借助EGM2008重力场模型构建海域重力异常模型的反距离插值算法。以DTU10重力异常数值模型为基础数据,设计了海域重力测量数据格网化实验,结果表明本文算法能有效改进海域重力异常格网模型的计算精度。     

利用Kmeans-RBF神经网络提高重力格网化数据精度

针对重力测量数据在格网化过程中精度会被降低的问题,顾及空间重力异常和地形的强相关性,提出了三维Kmeans-RBF神经网络方法,该方法利用神经网络的复杂非线性映射学习能力进行推估建模,并在模型训练和推估时加入地形数据作为物理控制。最后基于美国爱达荷州地区的实测重力数据进行验证,实验结果表明:该方法相对于二维Kmeans-RBF神经网络方法和直接进行拟合推估的Kriging方法,实验区内精度分别提高了24.85%和44.84%。     

测高重力场反演海底地形方法比较

根据重力异常与海底地形的相关性,利用日本西之岛附近海域内的重力异常和船测水深,分别采用SandwellSmith方法(SS,1994)和重力地质法(GGM)反演了该海域内1'×1'海底地形模型,并将反演的模型与格网化船测水深模型和ETOPO1模型进行了对比分析。在GGM中,使用了迭代法向下延拓确定密度差常数为5.7g/cm~3。与SS模型相比,GGM模型与格网化船测水深模型的差值的标准差达到了37.10m。     

基于少量地面控制的航空重力测量

提出6种少量地面重力测量布设方案和4种基于少量地面控制的航空重力测量数据处理方法,建立基于最小二乘配置的格网平均重力异常的融合模型,有效地估计出了航空重力数据可能存在的系统偏差,提高了格网数据的质量。     

几种插值法在海洋磁力测量数据格网化中的应用

为获得反映海区磁异常分布的磁异常图,海洋磁力测量数据处理中需要将离散数据通过插值方法进行格网化。给出了几种常用的数据格网化方法,并采用磁异常均匀分布和复杂分布情况对方法有效性进行了评价。结论表明:磁异常分布越复杂,对插值方法的要求越严格,无论是磁异常均匀分布还是复杂分布情况,Kriging法的效果最好,建议在海洋磁力测量数据处理中使用。     

利用卫星测高垂线偏差计算其他重力场元的特性分析

利用卫星测高技术确定海洋重力场,垂线偏差数据作为导出观测量在实际工作中被普遍采用。利用物理大地测量边值问题的定义以及扰动位在球面边界条件下的解,给出了由垂线偏差计算大地水准面高、重力异常和扰动重力的公式。分析了不同积分计算公式在重力场阶谱表达形式下对垂线偏差误差的抑制作用,也分析了不同积分核函数的变化特性,得出基本结论:在利用卫星测高数据求解海洋重力场时,当以格网化海面垂线偏差数据计算重力场参数时,求解的大地水准面高的有效性和稳定性优于重力异常和扰动重力。     

Surfer软件在海洋倾倒区地形监测与管理中的应用——以连云港倾倒区为例

针对海洋倾倒区地形测量频次高、数据多的特点,利用Surfer软件对连云港倾倒区多次地形测量数据进行处理、计算、分析,以了解不同时段内水动力影响下连云港倾倒区地形的冲淤变化过程和特点,根据Surfer软件的特点和倾倒区管理的现状,探讨Surfer软件在海洋倾倒区监测与管理中的应用.     

基于GIS组件和Surfer软件的海洋调查数据等值线绘制方法

论文讨论了利用海洋调查数据自动绘制等值线图的方法。该方法采用组件化GIS实现研究站位的交互选取,根据交互选取的站位信息访问后台数据库取得站位数据,使用Surfer软件作为后台服务器接收站位数据并绘制等值线图。该方法集成了GIS组件的交互选取功能和Surfer软件的空间插值和绘图功能,实现了海洋空间要素分布图的自动绘制,而且操作简单、直观,等值线图效果好。文章给出了该方法的主要步骤和调用Surfer实现插值和绘图的关键VB语句,并给出了该方法的一个具体应用实例。     

重力数据融合与重力垂直梯度异常反演

重力垂直梯度异常反应了重力异常的空间变化率,在地球物理勘探等多学科中得到越来越多的应用。利用南海局部区域实测重力异常数据和Sandwell测高重力异常数据,将搜索范围、距离和精度多种因素融合考虑并对Shepard算法进行改进,给出了南海局部区域(19°N~20.5°N,114°E~115.5°E)分辨率1'×1'的重力异常并反演了对应分辨率的重力垂直梯度异常。结果表明,基于Shepard改进算法的高精度船测重力和测高重力的有机融合,增强了单一测高重力数据反演重力垂直梯度异常的细节纹理,提高了反演重力垂直梯度异常的分辨率和精度。     

卫星测高数据格网化方法研究

用连续曲率张力样条方法和Shepard方法分别对多颗卫星和单颗卫星的海面高数据进行格网化,得到平均海面高。将平均海面高与DNSC08、KMSS04、CLS01平均海面高模型进行了比较,分析了两种格网化方法的优劣。结果表明,连续曲率张力样条方法比Shepard方法更适合于卫星测高数据的格网化。     

不同前处理方法对东沙底质沉积物土粒相对密度的影响

通过采用不同方法对东沙底质沉积物进行了土粒相对密度分析前处理,以探讨东沙底质沉积物中可溶盐、碳酸盐和有机质对其土粒相对密度的影响。结果表明,去除盐、碳酸盐和有机质后,东沙底质沉积物土粒相对密度明显降低,其中可溶盐对其土粒相对密度的影响最大。有机质含量的变化与东沙底质沉积物土粒相对密度变化的关系较大,呈负相关;而可溶盐和碳酸盐含量的变化与东沙底质沉积物土粒相对密度变化的关系则很小。     

Comparison of Satellite-Only Gravity Field Models Constructed with All and Parts of the GOCE Gravity Gradient Dataset

The impact of GOCE Satellite Gravity Gradiometer data on gravity field models was tested. All models were constructed with the same Laser Geodynamics Satellite (LAGEOS) and Gravity Recovery and Climate Experiment (GRACE) data, which were combined with one or two of the diagonal gravity gradient components for the entire GOCE mission (November 2009 to October 2013). The Stokes coefficients were estimated by solving large normal equation (NE) systems (i.e., the direct numerical approach). The models were evaluated through comparisons with the European Space Agency's (ESA) gravity field model DIR-R5, by GPS/Leveling, GOCE orbit determination, and geostrophic current evaluations. Among the single gradient models, only the model constructed with the vertical ZZ gradients gave good results that were in agreement with the formal errors. The model based only on XX gradients is the least accurate. The orbit results for all models are very close and confirm this finding. All models constructed with two diagonal gradient components are more accurate than the ZZ-only model due to doubling the amount of data and having two complementary observation directions. This translates also to a slower increase of model errors with spatial resolution. The different evaluation methods cannot unambiguously identify the most accurate two-component model. They do not always agree, emphasizing the importance of evaluating models using many different methods. The XZ gravity gradient gives a small positive contribution to model accuracy.     

增加近岸船载重力测量数据覆盖的推估方法研究

针对海洋区域离岸距离5～30km的范围内船载重力测量数据覆盖空白的现状,基于已有测线数据,对其进行不同空间距离采样形成对应的采样序列。利用动态时间规整算法计算其与初始测线数据的相关系数,依据相关系数与采样距离之间的关系,确定了最优重采样空间距离新方法。以最优重采样空间距离对测线数据进行重采样,利用拉格朗日插值算法,沿测线方向将测线数据向陆地推估。经过不同测线的内外部检核,结果表明船载重力测量向陆地方向扩展的保守距离约为5～10km,减少了船载重力测量数据在近岸海域覆盖空白的面积。本研究成果可为建立陆海一致垂直基准工作提供更全面的基础数据,技术方法可为航空重力、地磁等测线数据的精细处理及应用提供参考。     

基于频域法的留尼汪海域海底地形反演

海底地形是全球地形的重要组成部分,对地球物理科学研究、经济活动等具有重要作用。基于Parker公式,利用卫星测高重力异常和船测水深数据,采用频域的方法反演了疑似马航MH370失事区域的留尼汪海域的10°×10°的海底地形。最后将反演的水深和船测水深、国际通用的海深模型ETOPO1作比较进行精度评估,结果表明:本文反演结果与船测水深相比误差平均值为-26.038 m,标准差为176.588 m;与ETOPO1相比,差异平均值为-33.541 m,标准差为160.769 m。这表明采用重力异常数据,结合船测数据能较高精度地反演海底地形。     

不同来源重力数据在海洋重力编图中的拼接——以南通幅海洋重力编图为例

以南通幅海洋重力编图为例,论述了如何应用海洋重力与陆地重力以及卫星重力数据等不同来源的重力数据进行海洋重力编图。讨论了编图中关键问题,如各种数据的统一改算问题,包括统一的重力基准网、统一的正常重力公式、统一的投影方式;不同来源数据间的拼接调平问题;卫星重力数据使用范围问题等,并对最新的南通幅编图结果与以往相同范围内的编图进行了比较。     

Comparison of Satellite Altimetric Gravity and Global Geopotential Models with Shipborne Gravity in the Red Sea

The determination of high-resolution geoid for marine regions requires the integration of gravity data provided by different sources, e.g. global geopotential models, satellite altimetry, and shipborne gravimetric observations. Shipborne gravity data, acquired over a long time, comprises the short-wavelengths gravitation signal. This paper aims to produce a consistent gravity field over the Red Sea region to be used for geoid modelling. Both, the leave-one-out cross-validation and Kriging prediction techniques were chosen to ensure that the observed shipborne gravity data are consistent as well as free of gross-errors. A confidence level equivalent to 95.4% was decided to filter the observed shipborne data, while the cross-validation algorithm was repeatedly applied until the standard deviation of the residuals between the observed and estimated values are less than 1.5 mGal, which led to the elimination of about 17.7% of the shipborne gravity dataset. A comparison between the shipborne gravity data with DTU13 and SSv23.1 satellite altimetry-derived gravity models is done and reported. The corresponding results revealed that altimetry models almost have identical data content when compared one another, where the DTU13 gave better results with a mean and standard deviation of ?2.40 and 8.71 mGal, respectively. A statistical comparison has been made between different global geopotential models (GGMs) and shipborne gravity data. The Spectral Enhancement Method was applied to overcome the existing spectral gap between the GGMs and shipborne gravity data. EGM2008 manifested the best results with differences characterised with a mean of 1.35 mGal and a standard deviation of 11.11 mGal. Finally, the least-squares collocation (LSC) was implemented to combine the shipborne gravity data with DTU13 in order to create a unique and consistent gravity field over the Red Sea with no data voids. The combined data were independently tested using a total number of 95 randomly chosen shipborne gravity stations. The comparison between the extracted shipborne gravity data and DTU13 altimetry anomalies before and after applying the LSC revealed that a significant improvement is procurable from the combined dataset, in which the mean and standard deviation of the differences dropped from ?3.60 and 9.31 mGal to ?0.39 and 2.04 mGal, respectively.     

多类资料计算南海重力垂直梯度的方法比较

重力垂直梯度在解决和解释地球表层地质和地球物理问题中的作用日益明显,因而获得其模型和分布是非常必要的。利用测高卫星可以得到空间大范围高精度、高分辨率的垂线偏差、重力异常以及大地水准面数据,利用测高重力资料和地球重力场模型,采用不同方法分别计算了南海海域重力垂直梯度,并对它们进行了比较。     

KSS32-M型海洋重力仪动态性能分析

采用测网交点差、重复测线和与KSS31-M型海洋重力仪重合测线对比的方法,利用近年来KSS32-M海洋重力仪的实测数据对KSS32-M海洋重力仪测量稳定性和数据可靠性进行分析。利用机动转向法验证重力仪阻尼延迟时间为70 s,基于70 s阻尼延迟时间计算的重力测网的测量准确度为0.65 m Gal,与KSS31-M型海洋重力仪采集的重力剖面对比结果看,重合测线相关性为高度相关,4条重合测线网的交点差绝对值最大为1.66 m Gal,准确度为0.59 m Gal,均达到国家标准要求的近海重力测网交点差均方根小于2 m Gal的技术指标。重复测线的幅值接近,相位吻合,匹配测点异常差的平均值小于0.9,均方根均小于0.8,相关性均在0.98以上。本研究表明KSS32-M型海洋重力仪动态测量性能稳定、测量数据可靠。     

Combined satellite altimetry and shipborne gravimetry data processing

The recovery of quantities related to the gravity field (i.e., geoid heights and gravity anomalies) is carried out in a test area of the central Mediterranean Sea using 5' × 5' marine gravity data and satellite altimeter data from the Geodetic Mission (GM) of ERS‐J. The optimal combination of the two heterogeneous data sources is performed using (1) the space‐domain least‐squares collocation (LSC) method, and (2) the frequency‐domain input‐output system theory (IOST). The results derived by these methods agree at the level of 2 cm in terms of standard deviation in the case of the geoid height prediction. The gravity anomaly prediction results by the same methods vary between 2.18 and 2.54 mGal in terms of standard deviation. In all cases, the spectral techniques have a much higher computational efficiency than the collocation procedure. In order to investigate the importance of satellite altimetry for gravity field modeling, a pure gravimetric geoid solution, carried out in a previous study for our lest area by the fast collocation approach (FCOL), is used in comparison with the combined geoid models. The combined solutions give more accurate results, at the level of about 15 cm in terms of standard deviation, than the gravimetric geoid solution, when the geoid heights derived by each method are compared with TOPEX altimeter sea surface heights (SSHs). Moreover, nonisotropic power spectral density functions (PSDs) can be easily used by IOST, while LSC requires isotropic covariance functions. The results show that higher prediction accuracies are always obtained when using a priori nonisotropic information instead of isotropic information.     

Applying Iterative Method to Solving High-Order Terms of Seafloor Topography

Abstract

We introduce an iterative inversion method to address the problems in high-order seafloor topography inversion using gravity data (gravity anomaly and vertical gravity gradient anomaly), such as the difficulty in computing the equation and the uniqueness of the calculation results. A part of the South China Sea is selected as the experimental area. Considering the coherence and admittance function of gravity topography and vertical gravity gradient topography, the inversion band of the gravity anomaly and vertical gravity gradient anomaly in the study area is 30?km–120?km. Seafloor topography models of different orders are constructed using an iterative method, and the performance of each seafloor topography model is analyzed against ETOPO1 and other seafloor topography models. The experimental results show that as the inversion order increases, the clarity and richness of seafloor topographic expression continuously improve. However, the accuracy of seafloor topography inversion does not improve significantly when the inversion order exceeds a certain value, which is related to the contribution of high-order seafloor topography to gravity information. The results show that the accuracy of BGT4 (inversion model constructed by the gravity anomaly) is slightly poorer than that of BVGGT4 (inversion model constructed by the vertical gravity gradient anomaly) in areas with complex topography, such as multi-seamounts and trenches, and the results are generally better in areas with flat seafloor topography.