Geological structure of the yellow sea area from regional gravity and magnetic interpretation

Regional gravity and aeromagnetic data covering the area of 32°- 38° N, 118°-127° E at the scale of 1：1,000,000 are coordinated and integrated in a synthetic study of the South China Yellow Sea and adjacent areas. Depth to magnetic crystalline basement and its structure are determined by magnetic anomaly inversion. Depth to and thickness of the Paleozoic rock are also revealed by gravity anomaly inversion with constrains of the basement and known seismic information from several profiles. Structure units, main faults, basin boundaries, and sub-suppressions are outlined on the basis of gravity data interpretation.     

Estimating Field Source Parameters of Gravity Change in North China Based on the Euler Deconvolution Method

Based on the absolute and relative gravity observations in North China from 2009 to 2014, spatial dynamic variations of the regional gravity field are obtained. We employed the Euler deconvolution method and the theoretical model to get the best estimates of parameters. Gravity field change caused by the depth and distribution in North China is calculated by back analysis. The results show the structural index that equals 1 is suitable for inversion of the gravity variation data. The inversion results indicate that the depths of anomaly field sources are spread over the Hetao fault. The research method of this paper can be used in the quantitative study on the field source and may shed new light on the interpretations of gravity change, and also provide quantitative basis for earthquake prediction index criterions based on the gravity change.     

The inversion of density structure by graphic processing unit （GPU） and identification of igneous rocks in Xisha area

Organic reefs, the targets of deep-water petroleum exploration, developed widely in Xisha area. However, there are concealed igneous rocks undersea, to which organic rocks have nearly equal wave impedance. So the igneous rocks have become interference for future exploration by having similar seismic reflection characteristics. Yet, the density and magnetism of organic reefs are very different from igneous rocks. It has obvious advantages to identify organic reefs and igneous rocks by gravity and magnetic data. At first, frequency decomposition was applied to the free-air gravity anomaly in Xisha area to obtain the 2D subdivision of the gravity anomaly and magnetic anomaly in the vertical direction. Thus, the distribution of igneous rocks in the horizontal direction can be acquired according to high-frequency field, low-frequency field, and its physical properties. Then, 3D forward modeling of gravitational field was carried out to establish the density model of this area by reference to physical properties of rocks based on former researches. Furthermore, 3D inversion of gravity anomaly by genetic algorithm method of the graphic processing unit(GPU) parallel processing in Xisha target area was applied, and 3D density structure of this area was obtained. By this way, we can confine the igneous rocks to the certain depth according to the density of the igneous rocks. The frequency decomposition and 3D inversion of gravity anomaly by genetic algorithm method of the GPU parallel processing proved to be a useful method for recognizing igneous rocks to its 3D geological position. So organic reefs and igneous rocks can be identified, which provide a prescient information for further exploration.     

THE DISPERSION RELATION, OF EARTH''S ELECTROMAGNETIC RESPONSE FUNCTIONS AND THE JOINT INVERSION OF DATA

On the basis of the dispersion relations of MT field, the necessity and applied prospects of the joint inversions using a pair of MT response functions which are correlative with the dispersion relations, are infered. A filter coefficient algorithm is made, with which the corresponding impedance phase data can be estimated using a set of apparent resistivities. The tests for the observed MT data show that when comparing the impedance phase estimated using the dispersion relation with the ob served phase, it can be checked whether the dispersion relation between observed apparent resistivity and phase data is satisfied or not, and that the use of the phase data corrected using the dispersion relation in the joint inversion is advantageous to obtain more confident results. It is shown that joint inversions are more advantageous than single parameter inversions, and that in the most case the joint inversion using the apparent resistivities of impedance real and imaginary parts is more advantageous than the jointinversion using the normal apparent resistivity and impedance phase. The existence of the dipersion relations between the ratio apparent resistivity and corresponding impedance phase of the orthogonal electric and magnetic field horizontal Components in the frequency EM sounding with horizontal electric dipole(FEMS) are discussed, the better effect of the joint inversion using the pair of EM response functions is obtained. The problems on the one-dimensional joint inversion for the MT and FEMS apparent resistivities, for which the observed frequency bands partly overlape each other, are studied. It is shown that this joint inversion is applicable and effective:the joint inversions of the practical data for two kinds of EM methods at two sites give the results well corresponding to the drilling data. The simulated MT inversions for the data of two kinds of EM methods are made, and more confident results also are obtained.     

依据密度分布频率特性的重力反演方法研究

Three-dimensional gravity inversion based on the mass property model is very popular in recent years. The time and efficiency of inversion algorithms is relative to the magnitude of the target mesh. One approach is to search over the entire solution space for a more refined result. However, the inversion will be difficult with the increased parameters in the large search space and the number of computations increases exponentially. ｜n this paper, we propose a novel approach based on the frequency characteristics of the density distribution over the mesh. The purposes of our study are to reduce the parameters of three- dimensional gravity inversion and to lighten the image quality of the inversion result. The results show that the new method can expedite the inversion processing and get a better geological interpretation than tradition methods.     

玉树MS 7. 1 地震同震效应模拟与绝对重力检验

Using the results of aftershock relocation, inversion on seismic waves and InSAR results, and surface rupture displacements obtained by geological survey after the earthquake, this paper constructs a fault model of the Yushu Ms7. 1 earthquake. Based on rectangular dislocation theory in an elastic-viscoelastic layered medium, we have simulated the co- seismic deformation and gravity change with gravitational effect considered. The pictures show that the absolute gravity measuring point is beside the extremum of coseismic gravity change, and the numerical value reaches 25.02 x 10-Sm. s-2. After a discussion about the gravity changes before the earthquake and the coherence consistency between two FG-5 absolute gravimeters, we think that the measured value 27.2 × 10^-8 m· s^-2 at Yushu station is coseismic gravity change. It＇s coincident with the simulation results based on dislocation theory. Therefore it is a good tool to test the near-field changes found by dislocation theory.     

Wavefield continuation datuming using a near surface model

When topography and low velocity zone differences vary greatly, conventional vertical static time shifts will cause wavefield distortion and influence wave equation seismic imaging for seismic data acquired on a complex near surface. In this paper, we propose an approach to datum correction that combines a joint tomography inversion with wavefield continuation to solve the static problem for seismic data on rugged acquisition topography. First, the near surface model is obtained by refracted wave tomography inversion. Second, the wavefield of sources and receivers are continued downward and upward to accomplish datum correction starting from a flat surface and locating the datum above topography. Based on the reciprocal theorem, Huygens＇ and Fresnel principles, the location of sources and receivers, and regarding the recorded data on the surface as a secondary emission, the sources and receivers are upward-continued to the datum above topography respectively. Thus, the datum correction using joint tomography inversion and wavefield continuation with the condition of a complex near surface is accomplished.     

Review of Mercury's dynamic magnetosphere:Post-MESSENGER era and comparative magnetospheres

This review paper summarizes the research of Mercury’s magnetosphere in the Post-MESSENGER era and compares its dynamics to those in other planetary magnetospheres,especially to those in Earth’s magnetosphere.This review starts by introducing the planet Mercury,including its interplanetary environment,magnetosphere,exosphere,and conducting core.The frequent and intense magnetic reconnection on the dayside magnetopause,which is represented by the flux transfer event"shower",is reviewed on how they depend on magnetosheath plasma β and magnetic shear angle across the magnetopause,followed by how it contributes to the flux circulation and magnetosphere-surface-exosphere coupling.In the next,Mercury’s magnetosphere under extreme solar events,including the core induction and the reconnection erosion on the dayside magnetosphere,as well as the responses of the nightside magnetosphere,are reviewed.Then,the dawn-dusk properties of the plasma sheet,including the features of the ions,the structure of the current sheet,and the dynamics of magnetic reconnection,are summarized.The last topic is devoted to the particle energization in Mercury’s magnetosphere,which includes the energization of the Kelvin-Helmholtz waves on the magnetopause boundaries,reconnection-generated magnetic structures,and the cross-tail electric field.In each chapter,the last section discusses the open questions related to each topic,which can be considered by the simulations and the future spacecraft mission.We end this paper by summarizing the future Bepi Colombo opportunities,which is a joint mission of ESA and JAXA and is en route to Mercury.     

Precollisional tectonics and terrain amalgamation offshore southern Taiwan:Characterizations from reflection seismic and potential field data

Sponsored by the Chinese National Fundamental Research and Development Program in 2001,Guang-zhou Marine Geological Survey launched out a long geophysical survey from the northeastern part of the South China Sea (SCS),through the Luzon Arc,to the Huatung Basin and the Gagua Ridge. Based on high-resolution seismic data from this survey,combined with gravimetric and magnetic modeling,a systematic effort is made to the study of the regional geodynamics offshore southern Taiwan. By focusing particularly on precollisional tectonic interactions between adjacent geological units and their tectonic affiliations,this study can help reveal early arc-continent collisional processes that formed the Taiwan orogen. The construction of the Manila accretionary prism and its eastward progressive deformation indicate that the subduction of SCS have experienced multiple phases of increased activity. Active precollisional crustal shortening within the Northern Luzon Trough resulted in tilting of sedimentary layers at angles between 6° and 13°. But the shortening induced by tilting accounts for only a tiny part of regional total crustal compression. The eastern flank of the Luzon Arc appears to be more active than the rest,evidenced by active faulting and folding in the intra-arc basins on the eastern flank. Magnetic modeling/inversion shows that the Luzon Arc may have experienced multiple phases of magmatic activities,causing lateral magnetic inhomogeneity. Bouguer gravity anomalies and gravity modeling indicate that the Huatung Basin has anomalously higher crustal and upper mantle densities than those of SCS and the Luzon Arc. In addition,there is a large bathymetric difference between the Huatung Basin and the northeastern part of SCS basin. These observations argue against early hypothesis that the Huatung Basin and the northeastern part of SCS basin may once have belonged to one single oceanic crust,in part or in whole. The Gagua Ridge,as a sliver of uplifted oceanic crust,may be related to a transient northwestward subduction of the western Philippine plate. All evidences point to the argument that the region offshore southern Taiwan is experiencing multiple terrain amalgamation,which is a classical model for continental growth.     

Tectonic and climatic controls on late Quaternary bajada sedimentation along Narmada-Son Fault(NSF),Gujarat,Western India

An alluvial bajada has been recognized along a part of the active Narmada-Son Fault(NSF) and confined by the Karjan River on the eastern side and by the Madhumati River on western side.The bajada sequence exposed along the incised cliffs of various north flowing parallel streams has been studied in their proximal,medial and distal part. The sediments are characterized on the basis of grain size,fabric,sedimentary structures,bed geometry and sorting, and lithologs prepared from the mapped cliff sections.Detailed study of sedimentary characteristics,lithofacies analysis and facies associations indicate that sediments are mainly of three types-sediments that were deposited by debris flows and sediment gravity flows,and as extensive bar deposits.Three major aggredational sequences are recognized.Each sequence is composed of coarsening-upward sequence of proximal facies overlain by fining-upward sequence of distal facies.Coarsening upward sequence record periods of tectonic activity related to uplift along the NSF and fan progradation,whereas fining-upward sequence results from tectonic quiescence periods.The presence of rhizocretions,calcium carbonate nodules and calcite sheets within the basal debris flow and sediment gravity flow indicate semi-arid climate,whereas the formation of pedoginized paleosol indicate relatively wetter climate prevailing in the study area.Tectonic activity along the NSF has played dominant role by controlling the geometry and volume of bajada sediments.Climate is found to be responsible for compositional and temporal distribution of bajada sediments.OSL dating suggests that the bajada sediments were deposited during the later part of late Pleistocene.The incision of the sediments is attributed to uplift due to inversion of the lower Narmada basin during the early Holocene.     

基于交叉梯度约束的重力、磁法和大地电磁三维联合反演

为了降低单一地球物理方法反演的多解性及受噪声的影响程度,本文围绕重力、磁法和大地电磁法开展了三维联合反演的研究.重、磁采用基于对数障碍法的正则化反演算法,大地电磁使用limited-memory BroydenFletcher-Goldfarb-Shanno(L-BFGS)反演算法,引入交叉梯度函数实现了三种物性结构的相互耦合,最终开发出一套重磁电三维联合反演算法,并实现MPI并行加速计算.通过理论模型算例验证了算法的准确性,结果表明:不论是单棱柱体模型还是组合棱柱体模型,联合反演结果相较单独反演对于异常体的空间形态刻画以及物性数值恢复具有较好的提升;单棱柱体模型算例使得异常体的物性参数(密度、磁化率和电阻率)更加接近于真实的物性参数;组合棱柱体模型的联合反演结果不仅仅消除了围岩物性参数的假异常,而且还增强了异常体边界结构的恢复程度.     

南北地震带南段地壳厚度重震联合最优化反演

重力反演方法是研究地壳结构和物性界面起伏的有效地球物理手段之一.本文收集了南北地震带南段67个已有的固定台站接收函数反演的Moho面深度结果,并使用基于EGM2008重力异常模型计算的布格重力异常,验证了本文提出的重震联合密度界面反演方法的有效性.利用接收函数对台站下方Moho面深度估计作为先验约束,定义了一类评价函数,通过对重力反演算法中尺度因子,平移因子和稳定性因子的最优选择,最小化重力反演结果与接收函数模型之间的差异.结果表明,本文提出的方法,可以有效地同化不同地球物理方法获得的反演模型,且通过重震联合反演可以改进由于对空间分布不均匀的接收函数结果插值可能而引起的误差.本文还通过引入Crust1.0的Moho面深度为初值,同时考虑地壳密度的横向不均匀分布,通过模型之间的联合反演有效改善了地球物理反演模型间的不一致性问题.本文反演得到的最优化Moho面深度模型与已知67个台站位置接收函数模型之间的标准差约1.9km,小于Crust1.0与接收函数结果模型之间标准差为3.73km的统计结果.本文研究结果对于同化重震反演结果、精化地壳密度界面模型,都具有十分重要的参考意义.     

空-地-井重力异常正则化协同密度反演方法

围绕综合利用空-地-井多源重力异常联合反演提升精度的目标,提出正则化协同密度反演方法,从而有效利用多维数据的横向和纵向变化特征提高反演精度,且无需先验信息的约束.此外,还提出利用奇异值谱和深度分辨率工具评价航空数据观测高度和层数对反演分辨率的影响.通过理论模型试验本文方法在空-地、地-井等不同情况下的应用效果,结果表明空-地-井重力异常正则化协同密度反演方法能获得更高分辨率、高精度的反演结果,且证明不同的钻孔位置对反演结果有较大的影响,从而可指导实际的空-地-井联合勘探.最后,利用文顿盐丘实际勘探数据进行空地和空地井重力数据协同反演,反演结果垂向分辨率明显优于地面观测数据反演结果,盐盖的顶面及中心埋深与前人研究和地质资料解释相吻合,验证了方法正确性及实用性,为推进我国空-地-井立体重磁勘探提供了重要的技术手段.     

基于多台站的接收函数和重力联合反演确定莫霍面起伏和地壳平均波速比

地壳厚度和波速比是研究地壳结构和组分的两个重要参数,可为区域构造研究提供重要约束。接收函数被广泛地用于确定地壳厚度和波速比,例如H-κ方法或H-κ-c方法,但是该方法只能确定台站下方的地壳厚度和速度比,当地震台站分布稀疏时,很难约束台站间的横向不均匀性。另一方面,重力数据也可用于莫霍面的起伏变化研究,它在横向上覆盖很好,有较高的分辨率,但在纵向上分辨率相对较低。为此,本研究提出了一种联合反演算法求解莫霍面深度和地壳波速比参数。联合反演算法综合考虑了接收函数在纵向上的较高分辨率和重力数据在横向上的较高分辨率,同时拟合区域内所有台站上的接收函数和区域重力数据。模型测试表明联合反演算法较单一的接收函数反演更精确,特别是对于地壳厚度的确定。     

基于深度学习的重力异常与重力梯度异常联合反演

高效高精度的反演算法在重力大数据时代背景下显得尤为重要,受深度学习卓越的非线性映射能力的启发,本文提出了一种基于深度学习的重力异常及重力梯度异常的联合反演方法.文中首先提出了一种基于网格点几何格架的重力异常及重力梯度异常的空间域快速正演算法,这为本文深度学习反演算法的实现奠定了基础;其次对大量的不同密度模型进行正演计算获得样本数据集;然后设计了一种端到端的深度学习网络结构(GraInvNet),再利用样本数据对该网络结构进行训练;最后进行反演预测.组合模型试验表明,多维度数据联合反演相比单一分量反演其结果更“聚焦”,且与模型边界高度吻合,并且对于复杂模型的姿态与物性预测具有极为显著的优势,以及对于含噪声数据的反演,其质量也不会降低;Vinton岩丘实测重力数据也验证了文中方法的有效性;从而证明了深度学习在重力数据的高效高精度反演方面具有的巨大潜力.     

密度和速度随机分布共网格模型的重力与地震联合反演

针对重力与地震联合反演存在的问题,结合已有的研究成果,本文研究实现了速度和密度随机分布共网格单元模型的建模技术,以适应密度和速度剧烈变化的复杂模型及联合反演的计算要求.重力正演利用了该网格的二度半体模型,并进一步改进了地震走时的二维射线追踪计算方法,以适用于速度随机分布的网格介质.结合改进的模拟退火算法,实现了这种共网格条件下的重力与地震资料的同步联合反演.模型试验证明了重力与地震联合反演可以准确确定复杂物性界面的密度和速度结构,适用于物性界面不完全一致和物性变化剧烈的复杂模型,并且联合反演结果要优于单独的重力反演.带先验信息约束下的实际资料的联合反演,进一步证明了该方法的适用性和效果,可提高反演精度并减少多解性.     

部分区域约束下的交叉梯度多重地球物理数据联合反演

交叉梯度联合反演方法通过对多种地球物理模型实现结构耦合,在岩石物性关系不确定的情况下,既能提高反演结果的可靠性,又能减少反演的多解性,还能减少不同方法解释结果之间的矛盾.当不同的模型观测数据覆盖范围不一致时,交叉梯度联合反演通常需要取出重叠区域数据进行联合反演,并且建模时还要扩展一些模型范围.本文首先提出并实现了部分区域约束下的交叉梯度多重地球物理数据联合反演算法;接着进行了算法的模型试算;最后,我们将该反演算法用于本溪—集安深部地质调查重磁电综合地质地球物理解释中.结果表明:该算法不但能在重叠区域内很好地恢复结构相似的模型,而且在非重叠区域与重叠区域的边界处仍然可以得到平滑变化的模型;在本溪—集安10号剖面所获得的结构上相似的电阻率、密度及磁化率模型较好地反映了该区的深部地质结构,对于确定深部地质体的性质提供了有力的证据.     

三维多层介质重力-地震同步联合反演

联合反演是地球物理勘探的重要解释手段,能够提高模型参数的反演精度.本文在归纳和分析重力与地震资料联合反演的研究和应用现状的基础上,利用三维多层介质模型的地震走时和重力正演公式,推导了地震走时和重力异常对界面深度的雅可比矩阵,实现了三维重力-地震同步联合反演界面成像.最后进行了数值理论模型模拟和实例计算,结果表明地震走时和重力同步联合反演很好的重建了三维多层介质界面.     

大地电磁、重力、磁法和地震初至波走时的交叉梯度二维联合反演研究

为了准确的探测和描绘地下复杂的地质结构,同时克服地球物理单一方法反演的多解性和单一参数反演模型的不一致性等问题,近年来基于交叉梯度联合反演的综合地球物理解释已经得到了广泛的关注和应用.本文首先研究了两种地球物理方法的交叉梯度联合反演算法,在此基础上,推导并实现了多种地球物理方法(大地电磁,重力,磁法,地震初至波走时)的多交叉梯度约束的二维联合反演算法;其次,我们设计了结构不一致模型和复杂模型,针对多物性联合反演算法的准确性和有效性进行了模拟试算,并对复杂模型的单独反演结果和联合反演结果进行了交叉梯度值和物性交会图的对比;最后,本文将成熟的卫星资料多光谱综合分析技术应用到联合反演中,将多物性参数反演模型结果图通过RGB(红-绿-蓝)模式进行合成,得到融合的RGB合成图.结果表明:通过对结构不一致模型和复杂模型的联合反演结果和单独反演结果的对比分析,可以得出联合反演得到的结果更接近真实模型,并从得到的交叉梯度值进一步证明了联合反演模型相似度高,也从物性交会图中得到联合反演的物性相关性更好的结论,反向证明了算法的正确性.最终从得到的RGB合成图像,我们可以更直观的分析反演结果,更有利于准确划分地下模型结构.     

The distribution of deep source rocks in the GS Sag: joint MT–gravity modeling and constrained inversion

The coal-bearing strata of the deep Upper Paleozoic in the GS Sag have high hydrocarbon potential. Because of the absence of seismic data, we use electromagnetic (MT) and gravity data jointly to delineate the distribution of deep targets based on well logging and geological data. First, a preliminary geological model is established by using three-dimensional (3D) MT inversion results. Second, using the formation density and gravity anomalies, the preliminary geological model is modified by interactive inversion of the gravity data. Then, we conduct MT-constrained inversion based on the modified model to obtain an optimal geological model until the deviations at all stations are minimized. Finally, the geological model and a seismic profile in the middle of the sag is analysed. We determine that the deep reflections of the seismic profile correspond to the Upper Paleozoic that reaches thickness up to 800 m. The processing of field data suggests that the joint MT–gravity modeling and constrained inversion can reduce the multiple solutions for single geophysical data and thus improve the recognition of deep formations. The MT-constrained inversion is consistent with the geological features in the seismic section. This suggests that the joint MT and gravity modeling and constrained inversion can be used to delineate deep targets in similar basins.