The limitation of curvature gravity gradient tensor for edge detection and a method for overcoming it

Eigenvalues of the curvature gravity gradient tensor (CGGT) have been utilized to detect and outline edges of geologic bodies. In this paper, we analyzed and discussed the application scope and the limitation of CGGT using a complex synthetic model. We found that the large eigenvalue only can be used to outline edges of positive density bodies, and the small eigenvalue only can be used to delineate edges of negative density bodies. However, in the actual geological situation, the positive and negative density bodies usually coexist simultaneously. Therefore, it has certain limitations to detect edges using this method directly. In view of the limitations of the method, we combined the gravity anomaly with large eigenvalues, and proposed a new improved method for the eigenvalues of CGGT. The improved method can be utilized to outline edges of causative sources in any case. It was tested on synthetic model data and real field data. All of the results have shown that the new improved method is effective for edge detection.     

喜玛拉雅“东构造结”地区特异重力场的探讨

跨越中、印、缅三国交界的青藏高原东南的喜玛拉雅“东构造结”地区（92°E~97°E，26°N~30°N）一半以上的面积尚没有重力测点，是重力数据空白区,故无法直接研究其重力场特征与深部地壳结构（构造）.本文分析了卫星重力异常的特性，提出应用卫星重力异常作为近似空间重力异常，并作布格改正后，得到的布格重力异常具有与该地区地形高程呈镜像相关的特征，可用以研究深部地壳结构.据三条重力剖面计算得到该地区三个地壳深部结构剖面的结果，给出青藏高原地壳厚度>70 km；喜马拉雅造山带为55 km左右；布拉马普特拉河谷盆地为33~35 km；那加山山脉地区为40~45 km，显示出三者为三个不同的构造单元.同时给出布拉马普特拉构造单元为相对高密度的刚性物质构成，随着印度洋板块向北运移，在碰撞、挤压下，插入青藏高原东南缘一带.导致该地带的强烈构造运动，和频发大、小地震.最后提出了几点认识和建议.     

南海重力异常的沉积层密度改正及其对区域构造特征分析的意义

南海位于太平洋板块、印澳板块和欧亚板块交汇处,自晚中生代以来历经张裂作用、海底扩张以及印藏碰撞、菲律宾海板块西向运动等构造事件的叠加改造,不仅形成了复杂多样的构造格局,而且堆积了厚薄不均的沉积层.为了考察沉积层密度改正对利用重力资料分析南海不同尺度构造特征的影响,本文利用南海各区域不同深度沉积层的地震波速度及钻孔密度等数据,建立了沉积层与沉积基底密度差随深度变化的二次函数关系式,并基于该关系式,计算了南海沉积层相对基底密度低而产生的重力异常值.结果显示,南海沉积层的重力异常值在海盆区介于-40~-60 mGal,而在堆积巨厚沉积物的莺歌海盆地可达到-135 mGal;相对于空间重力异常、布格重力异常,经沉积层重力异常改正后的地壳布格重力异常更能突出深部不同尺度的密度结构和莫霍面的起伏特征,其总水平导数模更突显了南海西北部红河断裂带的海上延伸;利用谱分析技术估算岩石圈强度时,经沉积层重力异常改正的地壳布格重力异常数据获得的岩石圈有效弹性厚度值更为符合地质实际,特别是在长条形的巨厚沉积区如莺歌海盆地和马来盆地.分析表明,重力异常的沉积层密度改正对揭示南海构造特征具有重要的意义.     

鄂尔多斯—中秦岭—四川东部的重力异常场与深部地壳结构

本文基于在跨越鄂尔多斯盆地、中秦岭造山带与渭河盆地、四川盆地东部长达1010km的综合地球物理探测剖面,实地采集的最新高精度重力位场数据,分析了不同构造单元的重力位场特征,构建了该剖面区域的地壳密度结构模型.进而通过分析地壳内部不同密度界面的横向差异特征及分区,确定了该剖面区域内的断裂分布,提出了鄂尔多斯盆地—中秦岭造山带—四川盆地的"盆-山-盆"型构造体系的新认识.     

Lineaments in the Shamakhy–Gobustan and Absheron hydrocarbon containing areas using gravity data

In this study, we purposed to investigate the edge of geostructures and position of existing faults of the Shamakhy–Gobustan and Absheron hydrocarbon containing regions in Azerbaijan. For this purpose, the horizontal gradient, analytic signal, tilt angle, and hyperbolic of tilt angle methods were applied to the first vertical derivative of gravity data instead of Bouguer gravity data. We obtained the maps that show the previous lineaments which were designated by considering the maximum contours of horizontal gradient, analytic signal maps, and zero values of tilt angle, hyperbolic of tilt angle maps. The geometry of basement interface was also modeled utilizing the Parker–Oldenburg algorithm to understand the sediment thickness and coherency or incoherency between the gravity values and basement topography. The lineaments were held a candle to most current tectonic structure map of the study area. It was seen that the techniques used in this study are very effective to determine the old and new lineaments in the Shamakhy–Gobustan and Absheron regions. The epicenter distribution of earthquakes within the study area supports the new lineaments which are extracted by our interpretation. We concluded that better comprehension of Azerbaijan geostructures and its effect on the large scale works will be provided by means of this study.     

黄骅坳陷横向构造转换带与基底三分结构的重磁证据

黄骅坳陷盆地结构及构造单元展布受深断裂及盆地基底的控制与影响,利用重磁异常可以较大范围地揭示出控制盆地构造格局的深断裂及基底结构特征.本文综合考虑研究区盆地基底的宏观磁性差异以及存在的剩磁影响,采用对磁化方向依赖性小的磁异常模量数据,研究了黄骅盆地区域横向构造转换带的展布与磁性基底结构.识别出五条NW向与优势构造走向高角度相交的区域横向构造转换带和两条呈“T” 型展布的隐伏深断裂;研究区基底具有三分性结构特征,隐伏深断裂可能为三分基底的拼合线;基底结构差异及区域横向构造转换带共同控制了黄骅盆地“南北分区,东西分带”的宏观构造格局及沉积构造单元的展布.     

基于卫星重磁资料的南海深部构造研究

南海深部构造对研究南海构造演化和油气勘探具有重要意义.本文对南海地区的自由空气重力异常进行布格校正、海水层校正和沉积层校正,得到布格重力异常,再对布格重力异常进行区域异常和局部异常分离,利用位场界面反演方法对区域布格异常进行反演计算得到研究区域的莫霍面深度分布;采用全变倾角化极方法对研究区域的卫星磁异常数据进行化极处理,并进一步对化极磁异常作向上延拓,得到延拓后化极磁异常结果.分析布格重力异常、莫霍面深度及化极磁异常特征,结合天然地震层析成像的证据,得到以下结论：推测南海北部陆缘的古俯冲带位置是从118.5°E,24°N沿北东向延伸至109°E,15°N;红河断裂入海后经过莺歌海盆地在海南岛南部转为南北向与越东断裂相接并延伸至万安盆地;推测中特提斯洋的部分闭合位置是从110°E,2°N到101°E,21°N.     

渭河盆地及邻区重力异常小波多尺度分解与解释

基于EGM2008重力场模型计算获得了渭河盆地及邻区布格重力异常。采用小波多尺度分解方法对布格重力异常进行了4阶小波逼近和小波细节分解,同时基于平均径向对数功率谱方法定量化地计算出1~4阶小波细节和小波逼近所对应的场源平均埋深。结合区域地质和地震资料,对获得的重力场结果进行分析,得到如下结论:①鄂尔多斯地块、渭河盆地、秦岭造山带3个一级构造单元的布格重力异常之间存在明显差异;构造区内部重力异常也存在横向的显著差异。布格重力异常的走向、规模、分布特征与二级构造区及主要的断裂具有一定的对应关系。②渭河盆地及邻区布格重力异常1~4阶细节对应4~23 km不同深度的场源信息,鄂尔多斯地块南缘东、西部的地壳结构存在明显的差异;渭河盆地凹陷、凸起构造区边界清晰,断裂边界与重力异常边界具有较好的一致性;秦岭造山带重力异常连贯性不好,东、西部重力异常变化特征表现出明显的差异。③渭河盆地及邻区布格重力异常分布与莫霍面埋深具有非常明显的镜像关系。渭河盆地及邻区地震主要分布在六盘山—陇县—宝鸡断裂带、渭河断裂与渭南塬前断裂交汇处、韩城断裂与双泉—临猗断裂交汇处。渭河盆地及邻区重力异常主要由中上地壳剩余密度体所影响,这可能是该区地震以浅源地震为主的主要原因。     

喜马拉雅“东构造结”地区的特异重力场与深部地壳结构

跨越中、印、缅三国交界的喜马拉雅“东构造结”地区（92°E～97°E，26°N～30°N）有一半以上的面积尚没有重力测点，是重力数据空白区,故无法直接研究其重力场特征与深部地壳结构（构造）.本文应用卫星重力异常资料作为近似空间重力异常，经计算给出的布格重力异常，其特征与该地区的地形高程呈很好的镜像相关.据此得到该区不同方位的3个地壳深部结构剖面.重力异常反演求得青藏高原地壳厚度>70 km；喜马拉雅造山带为55 km左右；布拉马普特拉河谷盆地为33～35 km；那加山山脉地区为40～45 km，即呈现出3个不同构造单元的展布.同时求得“东构造结"区由高密度的刚性物质构成，在印度洋板块的碰撞、挤压作用下呈向北运移，并插入青藏高原东缘.基于这样的构造格局和深层动力过程，导致了青藏高原东南缘和东北缘的强烈构造运动，大、小地震的频频发生和矿产资源的聚集.     

基于布格重力异常研究运城盆地的构造特征

利用运城及周边地区的布格重力资料,对布格重力数据进行小波分解得到不同阶次的重力异常信息。依据各个阶次的小波变换结果并结合该区域的地质构造环境,对小波变换细节图中存在的特征较为明显的布格重力异常带进行深入分析,对异常所揭示的构造现象在地壳不同深度的分布特性等进行分析研究。结果表明:小波变换细节与区内地壳内部的隆起、凹陷及断裂构造具有很好的对应关系,可以对运城盆地的构造特征及地壳结构做更加深入的研究。     

裂陷盆地基底双界面模式二维重力反演

裂陷盆地基底的起伏表现为非光滑的几何形态,传统的重力反演结果并不能很好地反映这种特点.此外,大多数情况下,重力观测面并不位于盆地上界面,应为单独的起伏观测面,盆地应为上界面和基底组成的双界面模式.基于此,本文研究了起伏观测面上裂陷盆地基底双界面模式二维重力反演方法.研究中假设沉积盆地的沉积层与基底的密度差随深度按双曲线规律变化.将沉积盆地的沉积层剖分成相邻的垂直柱体,其水平尺寸是已知的,顶面与沉积层上界面重合,底面深度代表基底的深度,即为要反演的参数.反演中引入全变差函数作为盆地模型的约束,使得反演结果呈现非光滑形态,符合裂陷盆地基底特征.为减小反演多解性,引入已知深度点作为约束.建立由重力数据拟合、已知深度约束及全变差函数组成的目标函数,采用非线性共轭梯度算法使目标函数最小化.模型试算结果表明该方法可反演裂陷盆地基底起伏,并通过调整正则化参数的值可反演坳陷盆地基底起伏.将该反演方法用于珠江口盆地惠州凹陷和运城-临汾裂陷盆地实际资料处理,其结果较好地反映了裂陷盆地基底起伏特征,为研究盆地构造、油气勘探等提供重要参考.     

The gravity and isostatic Moho in North China Craton and their implications to seismicity

The Bouguer gravity is the combination of field sources in different depths. Based on the multi-scale analysis of the Bouguer gravity, we can get the gravity anomaly caused by the Moho undulation. This study presents the various orders of approximation of gravity anomaly in North China Craton (NCC), the possible source depths with radial logarithmic power spectrum, and the relationship between the deep structure and gravity anomaly. Furthermore, we discuss the isostatic compensation about the Moho depth from gravity and deep seismic sounding profiles (DSS). The results show that: (1) the fourth approximation could have resulted from the Moho undulation, (2) in contrast to the isostatic Moho, the inverted gravity Moho and the DSS Moho show that most of NCC has been isostatically compensated, and (3) the isostatic compensation rate has some close relation to the seismicity.     

Optimal Design of a Gravity Survey Network and its Application to Delineate the Jabera-Damoh Structure in the Vindhyan Basin,Central India

Fractal dimension analysis was carried out for optimal designing of 2-D gravity survey network and to determine an optimum range of gridding interval to generate least aliased Bouguer anomaly maps. As a test case, this method has been successfully applied to the Jabera-Damoh region of the Vindhyan Basin, which is considered as a potential hydrocarbon bearing area. In particular, we aim to delineate accurately the lateral extent of a possible hydrocarbon bearing structure. To achieve this aim, fractal dimension of survey network was computed using 2-D distributions of observation points in the planning phase of the survey so that the optimum station spacing for gravity survey can be obtained. A range of optimum gridding interval for the gravity data set was suggested using the box-counting method of fractal dimension determination. Bouguer anomaly maps of the region are prepared utilizing the optimum gridding interval. For the first time, these anomaly maps clearly outline the gravity evidence of an anomalous rifted structure, which is bounded by parallel faults on either side. This structure is interpreted as a favorable basin for the occurrence of hydrocarbons. Another finding of this study has been the delineation of an apparently small ridge-like structure running east-west, dividing this basin in two parts. A subsurface geological model along a profile across the Jabera structure has also been presented.     

Interpretation of gravimetric and magnetic anomalies in the Cameros Basin (North Spain): combination of deep and shallow sources

Outstanding potential field anomalies (gravimetric and magnetic) in the Cameros Basin (N Spain) follow a WNW-ESE trend, parallel to the geological structures resulting from Mesozoic extension and Tertiary basin inversion. The positive Bouguer gravity anomaly (15 mGal) is interpreted as the result of a strong contrast between the density of Tertiary rocks of the foreland basin and the Paleozoic and Mesozoic rocks, combined with crustal thickening in the Iberian Chain with respect to the Ebro Basin. The dipolar magnetic anomaly, slightly shifted to the south with respect to the relative maximum of the Bouguer anomaly, can be interpreted as related to volcanic rocks within the basement, which are linked to Triassic rifting as witnessed by outcrops of basalts along the basin margins. An exhaustive analysis of rock properties (density, magnetic susceptibility and remanence) and basin geometry from other sources (seismic reflection profiles) allow to constrain variations in crustal thickness and the location of large-scale basement faults.     

南海深部构造特征及其地质意义:来自重磁位场反演的认识

莫霍面和居里面是认识深部过程重要的地质与地球物理界面.为了进一步理解南海深部构造活动与洋盆扩张的关系,本文以OBS剖面和深反射地震剖面作为约束,对卫星测高重力异常进行海水、沉积层影响校正,采用最小曲率位场分离方法消除局部密度体的重力影响,获取了反映莫霍面起伏的重力异常,并利用双界面模型重力场快速反演方法计算得到了南海地区莫霍面深度值.通过与居里面起伏的对比研究,发现南海莫霍面和居里面整体均表现为"洋盆浅、周缘深"的菱形特征,两者在洋陆转换区呈现明显的窄梯级带特征,反映了南海扩张期岩石圈的强烈伸展减薄、南北向构造拉张作用等深部构造过程.洋盆莫霍面和居里面的西南向楔形形态是对南海由东向西渐进式扩张的深部构造响应.洋盆南部莫霍面浅于北部,这与扩张中心逐渐向南迁移的特征一致,而洋盆居里面南深北浅的特征则可能与洋盆的简单剪切扩张方式以及洋盆北部的岩浆活动更活跃有关.南海地区莫霍面和居里面呈现交错叠置关系,南、北陆缘表现为明显的深部构造差异,说明南海为非对称式扩张.北部陆缘区居里面深度浅于莫霍面,而洋盆区和南部陆缘区居里面深于莫霍面,这与南、北陆缘性质的差异和南部陆缘复杂的中-新生代俯冲碰撞等构造演化相关,而洋盆区居里面深于莫霍面的现象推测与大洋上地幔橄榄岩蛇纹石化导致的岩石磁性增强有关.     

Magnetic basement: gravity‐guided magnetic source depth analysis and interpretation

We present a new integrated approach to the interpretation of magnetic basement that is based on recognition of characteristic patterns in distributions and alignments of magnetic source depth solutions above and below the surface of magnetic basement. This approach integrates a quantitative analysis of depth solutions, obtained by 2D Werner deconvolution of the magnetic data, with a qualitative evaluation of the Bouguer gravity anomalies. The crystalline/metamorphic basement and sedimentary cover have different origins, tectonic histories, lithologies and magnetic properties. These differences result in different geometries of magnetic sources associated with faults, fracture zones, igneous intrusions, erosional truncations, subcrop edges and other structural discontinuities. Properly tuned, 2D Werner deconvolution is able to resolve the intra‐sedimentary and intra‐basement magnetic source geometries into distinctly different distributions and alignments of calculated depth solutions. An empirical set of criteria, basement indicators, was developed for identification and correlation of the basement surface. The ambiguity of basement correlation with limited or non‐existent well control, which is common for onshore frontier and offshore explorations, can be reduced by incorporating the Bouguer gravity data into the process of correlation.     

利用重力/GPS联合观测数据计算地壳垂向构造应力的新方法

本文提出一种基于重力/GPS联合观测数据计算垂向构造应力的新方法.计算步骤如下:(1)通过重力/GPS联合观测数据计算布格重力异常;(2)依据布格重力异常数据推算莫霍面深度;(3)依据GPS观测数据,通过均衡理论计算均衡面深度;(4)依据莫霍面与均衡面之间剩余物质(壳幔物质密度差)所承受的附加浮力,计算地壳承载的垂向构造应力.本文利用上述构造应力新算法,计算了巴颜喀拉块体东边界及周边地区垂向构造应力分布,发现龙泉山断裂带以东地区垂向构造应力基本为零,龙泉山断裂带与龙门山断裂带之间地区垂向构造应力为正值,巴颜喀拉地块东部垂向构造应力为负值.鲜水河断裂带东南段周边蓄积了-40~-50 MPa的垂向构造应力,且梯度变化剧烈;松潘高原蓄积的垂向构造应力大约为-10~-20 MPa,相对较小.     

中国东北地区重力场研究——利用改进的三方向小子域滤波划分主构造线及大地构造单元

利用重力异常进行大地构造单元的划分是区域重力资料处理与解释的重要任务.为了精细地划分中国东北地区大地构造单元及其分界线位置,对三方向小子域滤波进行了改进,探讨了改进的三方向小子域滤波的合理性.利用小子域及改进的三方向小子域滤波对研究区布格重力异常进行处理,并与欧拉反褶积反演结果进行对比分析.结果表明,改进的三方向小子域滤波在识别平面断裂的应用中具有更好的效果.综合分析东北地区布格重力异常特征、小子域和改进的三方向小子域滤波结果,划分出延续性较好、规模较大的断裂18条.对研究区重力异常进行了分区,描述了一、二级异常的分区特点,给出了东北地区大地构造单元的基本格架.黑河-齐齐哈尔-白城等重要构造线平面位置的确定突破了传统的重力异常分区模式,指出东北地区NE向构造具有走滑特征,且走滑距离具有由西向东逐渐增大的特点.     

Integrated Analysis on Gravity and Magnetic Fields of the Hailar Basin, NE China: Implications for Basement Structure and Deep Tectonics

The Hailar Basin is one of the typical basins among the NE China Basin Groups, which is situated in the east of East Asia Orogene between the Siberia Plate and the North China Plate. Based on the detailed analysis of magnetic, gravity, petrophysical, geothermal and seismological data, we separate the Gravity and Magnetic Anomalies (GMA) into four orders using Wavelet Multi-scale Decomposition (WMD). The apparent depths of causative sources were then assessed by Power Spectrum Analysis (PSA) of each order. Low-order wavelet detail anomalies were used to study the basin’s basement structure such as major faults, the basement lithology, uplifts and depressions. High-order ones were used for the inversion of Moho and Curie discontinuities using the Parker method. The results show that the Moho uplifting area of the Hailar Basin is located at the NE part of the basin, the Curie uplifting area is at the NW part, and neither of them is consistent with the basin’s sedimentary center. This indicates that the Hailar Basin may differ in basin building pattern from other middle and eastern basins of the basin groups, and the Hailar Basin might be of a passive type. When the Pacific Plate was subducting to NE China, the frontier of the plate lying on the mantle transition zone didn’t pass through the Great Khingan Mountains region, so there is not an obvious magma upwelling or lithospheric extension in the Hailar Basin area. Finally, based on the seismological data and results of WMD, a probable 2D crust model is derived from an across-basin profile using the 2D forward modeling of the Bouguer gravity anomaly. The results agree with those from seismic inversion, suggesting WMD is suitable for identifying major crustal density interfaces.     

Fast mapping of magnetic basement depth,structure and nature using aeromagnetic and gravity data: combined methods and their application in the Paris Basin

Assessment of deep buried basin/basement relationships using geophysical data is a challenge for the energy and mining industries as well as for geothermal or CO₂ storage purposes. In deep environments, few methods can provide geological information; magnetic and gravity data remain among the most informative and cost‐effective methods. Here, in order to derive fast first‐order information on the basement/basin interface, we propose a combination of existing and original approaches devoted to potential field data analysis. Namely, we investigate the geometry (i.e., depth and structure) and the nature of a deep buried basement through a case study SW of the Paris Basin. Joint processing of new high‐resolution magnetic data and up‐to‐date gravity data provides an updated overview of the deep basin. First, the main structures of the magnetic basement are highlighted using Euler deconvolution and are interpreted in a structural sketch map. The new high‐resolution aeromagnetic map actually offers a continuous view of regional basement structures and reveals poorly known and complex deformation at the junction between major domains of the Variscan collision belt. Second, Werner deconvolution and an ad hoc post‐processing analysis allow the extraction of a set of magnetic sources at (or close to) the basin/basement interface. Interpolation of these sources together with the magnetic structural sketch provides a Werner magnetic basement map displaying realistic 3D patterns and basement depths consistent with data available in deep petroleum boreholes. The last step of processing was designed as a way to quickly combine gravity and magnetic information and to simply visualize first‐order petrophysical patterns of the basement lithology. This is achieved through unsupervised classification of suitably selected gravity and magnetic maps and, as compared to previous work, provides a realistic and updated overview of the cartographic distribution of density/magnetization of basement rocks. Altogether, the three steps of processing proposed in this paper quickly provide relevant information on a deep buried basement in terms of structure, geometry and nature (through petrophysics). Notwithstanding, limitations of the proposed procedure are raised: in the case of the Paris Basin for instance, this study does not provide proper information on Pre‐Mesozoic basins, some of which have been sampled in deep boreholes.