A VGGNet-based correction for satellite altimetry-derived gravity anomalies to improve the accuracy of bathymetry to depths of 6500 m

Understanding the topographic patterns of the seafloor is a very important part of understanding our planet.Although the science involved in bathymetric surveying has advanced much over the decades,less than 20%of the seafloor has been precisely modeled to date,and there is an urgent need to improve the accuracy and reduce the uncertainty of underwater survey data.In this study,we introduce a pretrained visual geometry group network(VGGNet)method based on deep learning.To apply this method,we input gravity anomaly data derived from ship measurements and satellite altimetry into the model and correct the latter,which has a larger spatial coverage,based on the former,which is considered the true value and is more accurate.After obtaining the corrected high-precision gravity model,it is inverted to the corresponding bathymetric model by applying the gravity-depth correlation.We choose four data pairs collected from different environments,i.e.,the Southern Ocean,Pacific Ocean,Atlantic Ocean and Caribbean Sea,to evaluate the topographic correction results of the model.The experiments show that the coefficient of determination(R~2)reaches 0.834 among the results of the four experimental groups,signifying a high correlation.The standard deviation and normalized root mean square error are also evaluated,and the accuracy of their performance improved by up to 24.2%compared with similar research done in recent years.The evaluation of the R^(2) values at different water depths shows that our model can achieve performance results above 0.90 at certain water depths and can also significantly improve results from mid-water depths when compared to previous research.Finally,the bathymetry corrected by our model is able to show an accuracy improvement level of more than 21%within 1%of the total water depths,which is sufficient to prove that the VGGNet-based method has the ability to perform a gravity-bathymetry correction and achieve outstanding results.

Robust inversion analysis of local gravity anomalies caused by geological dislocation model of faults

Theoretical analysis and practical observations show that fault dislocations can change the gravity field around the fault. Gravity changes which were caused by the repeated dislocations over a long period of time were superimposed on the Bougeur gravity anomalies. These anomalies became the evidence of historical movement of fault as well as provide a way for the study of paleo earthquakes. This paper investigates inversion methods for the geological dislocation modeling of faults using the local Bouguer's gravity anomalies. To remove the effects of the irrelevant part of gravity anomalies to fault movements, we propose the robust nonlinear inversion method and set up the corresponding algorithm. Modeling examples indicate that the Marquardt's and Baye's least squares solutions depart from the true solution due to the attraction of gross errors in the data. The more seriously the data is contaminated, the more seriously the solutions are biased. In contrast, the proposed robust Marquardt's and Baye's inversion solutions can still maintain consistency with the solution without gross errors, even though 50 percent of the data is contaminated. This indicates that the proposed robust methods are effective. Using the proposed methods, we invert the geological dislocation models of the faults around the Erhai Lake in West Yunnan. The results show that the Northern Cangdong fault and the Erhai fault are normal dip slip faults with about 4 to 5 km dislocations; and that the Southern Cangdong fault has a less dip slip compared with the former two. A satisfactory fitting between the theoretical values of the inversion solution and the actual local gravity field is achievable.

Gravity anomalies determined from mean sea surface model data over the Gulf of Mexico

With the improvements in the density and quality of satellite altimetry data,a high-precision and high-resolution mean sea surface model containing abundant information regarding a marine gravity field can be calculated from long-time series multi-satellite altimeter data.Therefore,in this study,a method was proposed for determining marine gravity anomalies from a mean sea surface model.Taking the Gulf of Mexico(15°–32°N,80°–100°W)as the study area and using a removal-recovery method,the residual gridded deflections of the vertical(DOVs)are calculated by combining the mean sea surface,mean dynamic topography,and XGM2019e_2159 geoid,and then using the inverse Vening-Meinesz method to determine the residual marine gravity anomalies from the residual gridded DOVs.Finally,residual gravity anomalies are added to the XGM2019e_2159 gravity anomalies to derive marine gravity anomaly models.In this study,the marine gravity anomalies were estimated with mean sea surface models CNES_CLS15MSS,DTU21MSS,and SDUST2020MSS and the mean dynamic topography models CNES_CLS18MDT and DTU22MDT.The accuracy of the marine gravity anomalies derived by the mean sea surface model was assessed based on ship-borne gravity data.The results show that the difference between the gravity anomalies derived by DTU21MSS and CNES_CLS18MDT and those of the ship-borne gravity data is optimal.With an increase in the distance from the coast,the difference between the gravity anomalies derived by mean sea surface models and ship-borne gravity data gradually decreases.The accuracy of the difference between the gravity anomalies derived by mean sea surface models and those from ship-borne gravity data are optimal at a depth of 3–4 km.The accuracy of the gravity anomalies derived by the mean sea surface model is high.

Three-dimensional gravity inversion based on 3D U-Net++

The gravity inversion is to restore genetic density distribution of the underground target to be explored for explaining the internal structure and distribution of the Earth.In this paper,we propose a new 3D gravity inversion method based on 3D U-Net++.Compared with two-dimensional gravity inversion,three-dimensional(3D)gravity inversion can more precisely describe the density distribution of underground space.However,conventional 3D gravity inversion method input is two-dimensional,the input and output of the network proposed in our method are three-dimensional.In the training stage,we design a large number of diversifi ed simulation model-data pairs by using the random walk method to improve the generalization ability of the network.In the test phase,we verify the network performance by using the model-data pairs generated by the simulation.To further illustrate the eff ectiveness of the algorithm,we apply the method to the inversion of the San Nicolas mining area,and the inversion results are basically consistent with the borehole measurement results.Moreover,the results of the 3D U-Net++inversion and the 3D U-Net inversion are compared.The density models of the 3D U-Net++inversion have higher resolution,more concentrated inversion results,and a clearer boundary of the density model.

Joint inversion of gravity and seismic data along a profile across the seismogenic fault of 2010 Yushu Ms7. 1 earthquake

Yushu Ms7.1 earthquake occurred on the Ganzi-Yushu fault zone, across which we carried out a joint relative-gravity and seismic-reflection survey, and then performed a gravity inversion constrained by the seismic-reflection result. Based on the data of complete Bouguer gravity anomaly and seismic reflection, we obtained a layered interface structure in deep crust down to Moho. Our study showed that the inversion could reveal the interfaces of strata along the survey profile and the directions of regional faults in two-dimension. From the characteristics of the observed topography of the Moho basement, we tentatively confirmed that the uplift of eastern edge of Qinghai-Tibet plateau was caused by the subduetion of the Indian plate.

Gravity inversion of deep-crust and mantle interfaces in the Three Gorges area

To better understand the heterogeneity of deep-crust and mantle interfaces in the region of the Three Gorges, China, we used the Parker-Oldenburg iterative inversion method to invert existing Bouguer gravity data from the Three Gorges area （ 1 ： 500000） , a new gravity map of the Three Gorges Dam （ 1 ： 200000） , and the results of deep seismic soundings. The inversion results show a Moho depth of 42 km be- tween Badong and Zigui and the depth of the B2 lower-crustal interface beneath the Jianghan Plain and sur- rounding areas at 21 -25 km. The morphology of crustal interfaces and the surface geology present an over- pass structure. The mid-crust beneath the Three Gorges Dam is approximately 9 km thick, which is the thin- nest in the Three Gorges area and may be related to the shallow low-density body near the Huangling anti- cline. The upper crust is seismogenic, and there is a close relationship between seismicity and the deep- crust and mantle interfaces. For example, the MS. 1 Zigui earthquake occurred where the gradients of the Moho and the B2 interface are the steepest, showing that deep structure has a very important effect on re- gional seismicity.

3D joint inversion of surface and borehole gravity data using zeroth-order minimum entropy regularization

Surface and borehole gravity data contain complementary information.Thus,the joint inversion of these two data types can help retrieve the real spatial distributions of density bodies.When a sharp boundary exists between an anomalous density body and its surrounding rock,the interface recovered by smooth inversion with Tikhonov regularization is not clear,leading to difficulties in the subsequent geological interpretation.In this work,we develop a joint inversion of surface and borehole gravity data using zeroth-order minimum entropy regularization.The method takes advantage of the complementary information from surface and borehole gravity data to enhance the imaging resolution of density bodies.It also produces a focused imaging of bodies through the zeroth-order minimum entropy regularization without requiring a preselection of a proper focusing parameter.We apply the developed joint inversion approach to three diff erent synthetic data sets.Inversion results show that the focusing inversion with the zeroth-order minimum entropy regularization provides a good description of the true spatial extent of anomalous density bodies.Meanwhile,the joint focusing inversion reconstructs a more reliable density model with a relatively high resolution when a density body is passed through by one or more boreholes.

A NEW METHOD TO DETERMINE THE START DEPTH IN GRAVITY ANOMALY INVERSION

The start depth in gravity anomaly inversion is discussed in this paper, and a newmethod called "test comparing method" is presented to determine the start depth when both boring dataand seismic data are not available. Anticline and synline density interface models are designed to test themethed; the calculated resall agreed well with the designed value, which Proved tha the new method canbe used to determine the start depth in the guity anomal inversion when boring data and seisndc dataare not available.

Effective elastic thickness of the lithosphere from joint inversion in western China and its implications

The western China lies in the convergence zone between Eurasian and Indian plates.It is an ideal place to study the lithosphere dynamics and tectonic evolutions on the continental Earth.The lithospheric strength is a key factor in controlling the lithosphere dynamics and deformations.The effective elastic thickness(T_(e))of the lithosphere can be used to address the lithospheric strength.Previous researchers only used one of the admittance or coherence methods to investigate the T_(e) in the western China.Moreover,most of them ignored the internal loads of the lithosphere during the T_(e) calculation,which can produce large biases in the T_(e) estimations.To provide more reliable T_(e) estimations,we used a new joint inversion method that integrated both admittance and coherence techniques to compute the T_(e) in this study,with the WGM2012 gravity data,the ETOPO1 topographic data,and the Moho depths from the CRUST1.0 model.The internal loads are considered and investigated using the load ratio(F).Our results show that the joint inversion method can yield reliable T_(e) and F values.Based on the analysis of T_(e) and F distributions,we suggest(1)the northern Tibetan Plateau could be the front edge of the plate collision of Eurasian and Indian plates;(2)the southern and part of central Tibetan Plateau have a strong lithospheric mantle related to the rigid underthrusting Indian plate;(3)the southeastern Tibetan Plateau may be experiencing the delamination of lithosphere and upwelling of asthenosphere.

Using WGM2012 to Compute Gravity Anomaly Corrections of Leveling Observations in China

Gravity Anomaly Correction(GAC)is a very important term in leveling data processing.In most cases,it is troublesome for field surveyors to measure gravity when leveling.In this paper,based on the complete Bouguer Gravity Anomaly(BGA)map of WGM2012,the feasibility of replacing in-situ gravity surveying in China is investigated.For leveling application,that is to evaluate the accuracy of WGM2012 in China.Because WGM2012 is organized with a standard rectangle grid,two interpolation methods,bilinear interpolating and Inverse Distance Weighted(IDW)interpolating,are proposed.Four sample areas in China,i.e.,Hanzhong,Chengdu,Linzhi and Shantou,are selected to evaluate the systems bias and precision of WGM2012.Numerical results show the average system bias of WGM2012 BGA in west China is about-100.1 mGal(1 mGal=10^(-5) m/s^(2))and the standard deviation is about 30.7 mGal.Tests in Shantou indicate the system bias in plain areas is about-130.4 mGal and standard deviation is about 6.8 mGal.All these experiments means the accuracy of WGM2012 is limited in high mountain areas of western China,but in plain areas,such as Shantou,WGM2012 BGA map is quite good for most leveling applications after calibrating the system bias.

高精度重磁方法寻找隐伏矽卡岩型铁矿

位于华北克拉通东部的鲁中莱芜地区是矽卡岩型富铁矿重要产区,矿床主要产于矿山岩体与中奥陶统碳酸盐岩地层的接触带中。本文利用最新重磁测量成果,对矿山岩体外围西部石家泉—刘家庙一带重磁异常特征进行了细致研究,结合区内已知铁矿重磁场特征,圈定了深部找矿靶区。在成矿有利地段布设了大比例尺重磁剖面,以已知钻孔为约束条件,利用2.5D重磁联合反演技术,对重磁异常进行了定性及定量解释。解释结果为后续钻孔位置布设及深度预测提供了依据,该钻孔揭露了15.8 m厚的富铁矿,找矿效果显著,为今后本区进一步寻找矽卡岩型铁矿提供了重要指示依据。

华南陆块西南缘上地幔顶部三维密度结构及其地质意义

华南陆块自新远古代拼合以来,经历了多期构造-岩浆活动,岩石圈及其深部结构较为复杂,前人对西南缘扬子与华夏地块缝合带位置的厘定仍有不同认识.研究该区域上地幔顶部密度结构,有助于进一步认识扬子—华夏地块拼合对岩石圈的改造,对确定扬子—华夏地块在广西地区的边界位置有重要意义.基于卫星布格重力异常,利用地震波速度-密度转换关系构建的上地幔密度参考模型作为约束,反演剩余的岩石圈和上地幔重力异常,获得了岩石圈至300 km深度的上地幔顶部三维密度结构.密度结构揭示了华南陆块西南部的岩石圈—软流圈界面起伏和上地幔物质密度的不均匀分布.分析结果表明广西中部存在高密度岩石圈物质下沉,南部低密度的软流圈物质上涌,取代已拆沉的岩石圈根并向上侵入地壳,是地表出露花岗岩的来源.结合区内基底特征和碎屑锆石陆缘来源,推断沿广西富川—武宣—隆安—凭祥一线为扬子—华夏地块在华南陆块西南缘(广西段)缝合位置.

Multi-scale inversion of density structure from gravity anomalies in Tarim Basin

Based on the power spectra of gravity anomalies in Tarim Basin, the anomalies can be decomposed to the following three components: a sub-anomaly formed mainly by the basin crystallized basement, a sub-anomaly formed mainly by deep sedimentary layers, and that by shallow sedimentary layers. A special wavelet transform analysis scheme and a density inversion method are designed and applied to the decomposition and inversion of gravity sub-anomalies, which are correlated with regional geology and drilling data. The results indicate that the deep and the shallow sub-anomalies show some relations with ancient fluid active zones. The negative density disturbances inversed from the shallow sub-anomaly are mainly caused by Mesozoic fluid active zones, whereas the negative density disturbances form the deep sub-anomalies are mainly correlated with Paleozoic fluid active zones. As the ancient fluid movement was good for formation of oil/gas fields, the multi-scale inversion methods for locating the ancient fluid movement zones seem to be helpful and valuable to oil/gas exploration.

四川长宁及周缘重力数据融合及深浅构造特征研究

高精度地球物理数据是获取地壳内部物质变化的基础,2019年6月17日四川长宁发生MS6.0级地震,近期该区域多次发生5级左右地震,且区域内地震的发震诱因存在争议。本研究以实测重力剖面数据为基准,首先寻找重力场模型数据与实测重力剖面数据之间的映射关系式;然后对重力场模型数据进行修正,采用双线性插值法获取研究区内高精度网格数据;以长宁MS6.0地震震中为中心点,提取4条重力剖面并采用密度反演获取跨长宁MS6.0地震震中的地壳密度差变化特征。研究结果表明:通过重力场模型数据与实测重力数据间的映射关系式计算所获取的融合数据精度小于5%;长宁及周边震群发震区域对应三叠纪地层且震群多沿背斜、向斜发生;长宁MS6.0地震位于剩余布格重力异常0~30 mGal的正低值区域和密度反演高低转换带上0.048~0.058 g/cm^(3)的正低值区域,该区域位于长宁背斜上的大地湾断层与双河背斜轴部隐伏断层的交汇处。本研究可为多源重力数据融合提供技术基础,服务于以动力学为基础的地震物理预测。

基于UNet++卷积神经网络的重力异常三维密度反演

三维密度反演是地球物理领域的研究热点,而在大数据及人工智能发展的时代背景下如何快速高效地实现重力数据反演显得更为重要.传统反演方法通常需要存储大型系数矩阵,内存占用大,耗费时间长,同时为约束反演结果而加入的正则化约束项参数难以确定;深度学习可以不依赖先验信息,也不需要计算及存储系数矩阵,使得计算效率大大提高.基于此,本文提出了一种基于UNet++网络的重力异常反演方法.首先将UNet++网络中部分参数进行更改:选择在输入数据绝对值较大时梯度更稳定的LeakyReLU作为激活函数;加入了Batch Normalization层,增强了网络的收敛速度及稳定性.然后为了提高网络的全局最优化能力,引入了基于余弦退火的学习率更新策略,使用梯度的一阶以及二阶矩估计的Adam最优化算法,利用数据集与标签集进行网络训练,实现了重力异常的三维密度反演.通过实验验证了UNet++、LeakyReLU更快速稳定的收敛能力,而余弦退火学习率更新策略具有更强的全局寻优能力.含噪模型实验及实际数据反演结果进一步证明该方法的正确性和有效性,及其良好的泛化能力与抗噪能力.

近震体波走时与重力联合反演华北地区地壳上地幔顶部P波速度结构

华北是克拉通破坏的典型案例,也是陆内强震活动最为强烈、频繁的地区之一,对该区壳幔速度结构进行研究,可为深化认识克拉通的形成演化、板内强震与火山发生机理等研究提供重要的深部约束.利用202个固定台站、360个流动台站所记录到的天然地震事件P波资料和实测布格重力异常,采用近震体波走时与重力联合反演方法获得了地壳、上地幔顶部P波速度结构.重力资料的引入不仅对壳幔结构进行了更合理的约束,而且有效提高了浅层和边缘射线分布稀疏区的成像能力.联合成像结果显示,鄂尔多斯地块和华北平原浅层表现为明显的低速异常,可能与其存在厚的松软沉积层有关;鄂尔多斯东部上地壳存在明显的P波高速异常,推测该高速异常与古元古代孔兹岩带的形成有关.大同火山区15~25 km和25~40 km表现为明显的低速异常,且这两个低速层是连通的,推测与地幔热物质上涌导致的多期次岩浆侵位有关.三河—平谷8.0级地震和邢台7.2级地震震源都位于脆性的上地壳,均下伏低速异常,推测在区域动力环境加载作用下,韧性的壳内蠕变导致脆性上地壳岩石中弹性应变能的局部积累,从而引起地震的发生.

Cascade Kalman filter for gravity anomaly distortion correction based on second order potential model

Using a gravity anomaly covariance function based on the second-order Ganssian Markov gravity anomaly potential model, the state equation of a gravity anomaly signal is obtained in marine gravimetry. Combined with the system state equation and the measurement equation, a new method of the cascade Kalman filter is proposed and applied to the correction of gravity anomaly distortion. In the signal processing procedure, an inverse Kalman filter is used to restore the gravity anomaly signal and high frequency noises first. Then an adaptive Kalman filter, which uses the gravity anomaly state equation as the system equation, is set to estimate the actual gravity anomaly data. Emulations and experiments indicate that both the cascade Kalman filter method and the single inverse Kalman filter method are effective in alleviating the distortion of the gravity anomaly signal, but the performance of the cascade Kalman filter method is better than that of the single inverse Kalman filter method.

基于地形单元分区的BP神经网络反演海底地形

海底地形和海洋重力场在一定波段频率域上存在相关性,传统海底地形反演方法主要考虑海底地形与重力数据之间的线性关系,而忽略非线性项作用对反演结果的影响。基于此,提出了基于海底地形单元分区的BP神经网络反演方法,改善海底地形反演精度。选取东印度洋北部东经90°海岭附近海域(83°E~92°E,10°S~10°N)为实验区,利用船测水深、卫星测高重力等数据,根据地形特征将实验区分为海盆区和海岭区,使用BP神经网络方法分别进行海底地形反演,构建了实验区1′×1′局部海底地形模型。结果表明,基于地形单元分区的反演结果模型平均相对误差精度可达1.45%,比未分区的反演结果均方根误差降低了42 m,验证了本文方法的有效性和可行性。

胶东招远-平度断裂地球物理异常特征及其控矿条件研究

招远-平度断裂是胶东金矿集区三条超大型控矿断裂(三山岛、焦家及招远-平度断裂)之一。受地表第四系覆盖影响,金矿体空间赋存位置与主断裂及其次级断裂地表出露位置的耦合关系尚存问题。通过开展精细化处理区域重力和航空磁测数据,如空间域不同高度延拓(100 m、200 m和500 m)及垂向一阶导数等方法,发现招远-平度断裂地表出露位置与深部地球物理梯度带存在位移差,且位移差向断裂倾向方向南东侧移动趋势愈发明显,进一步证明招远-平度断裂深部倾向南东,倾角由浅至深逐渐变缓,局部变化剧烈的地质推论。同时,开展地表高精度1/10000比例尺的重力结果显示,重力异常曲线自西向东整体呈快速上升趋势,局部存在凹陷区,表明花岗质碎裂岩及绢英岩化花岗质碎裂岩密度(2.63 g/cm^(3))介于玲珑二长花岗岩(2.57 g/cm^(3))和古元古界(2.82 g/cm^(3))之间,含矿层位与围岩存在较大密度差,金矿体主要赋存于花岗质碎裂岩及绢英岩化花岗质碎裂岩层位的倾角由陡变缓或由缓变陡的转折部位,可分别发现石英脉型金矿和蚀变岩型金矿。

基于残差神经网络的马里亚纳海沟地形反演

为提高利用重力异常数据反演马里亚纳海沟地形的精度,基于残差深度神经网络(residual deep neural network,RDNN)方法和重力异常等数据反演了马里亚纳海沟的1′×1′海底地形,通过实测检核点水深对RDNN模型的精度进行评估,并和重力地质法(gravity-geology method,GGM)模型结果进行对比。结果表明RDNN较GGM对马里亚纳海沟地形反演更为精细,实测水深检核表明RDNN模型均方根误差为128.98 m,优于GGM模型的150.14 m,且RDNN与船测检核水深也有更好一致性,RDNN深度学习模型为利用重力异常数据反演高精度海底地形提供了参考和依据。