东北地区重力均衡异常特征的初步研究

本文利用地面实测重力资料和地形高程资料,采用普拉特-海福特(Pratt-ttayford)重力均衡理论模型,取1°×1°方格网,通过使用现成改正表格查取改正值与个别计算点用理论公式计算作校核的方法,计算了我国东北地区75个计算点的均衡重力异常值;并对局部第四系覆盖较厚地区作了第四系密度改正;在此基础上,构制了我国东北N39°—49°,E121°—131°大部分地区的均衡重力异常图;结合区域布格重力异常和区域空间重力异常特征以及莫霍界面的起伏特点作了对比分析和讨论     

Regional gravity investigation of Honduras, Central America

A regional gravity study of Honduras was performed as part of a major study of the geothermal resource potential of Honduras. This study was conducted by Los Alamos National Laboratory, in cooperation with the Honduras government. Regional offshore free-air and onshore Bouguer gravity maps, and residual/isostatic gravity maps of Honduras and surrounding regions were produced. From these data several regional crustal and upper mantle models were produced. These models pass through two local geothermal sites, Platanares and San Ignacio. The regional geologic and tectonic implications of the models and their relevance to the geothermal potential of Honduras and to six well known geothermal sites in particular are examined. No obvious regional structures observed in the gravity data can explain the thermal enhancement in general or the specific geothermal sites. More local tectonic or structural conditions must control the distribution of the thermally enhanced areas.     

A REGIONAL GRAVITY STUDY OF THE NORTHERN NORTH SEA (56–62° N)*

A Bouguer gravity anomaly map is presented of the North Sea and adjacent land areas in Norway and Denmark, covering an area situated between 56° and 62°N, 1°W and 10°E. The gravity data from the UK sector of the North Sea, the land and offshore areas of Denmark, and the land areas of Norway have been published before. However, the gravity data from the Norwegian sector of the North Sea are new. A large number (about 60) of individual gravity features can be defined in the mapped area. Most of those situated in the UK sector of the North Sea and on land in Norway have been discussed earlier; however, most of the anomalies found elsewhere which are qualitatively interpreted here have not been discussed before. An interpreted Bouguer anomaly map is presented which identifies all these features. The majority of the gravity anomalies encountered in the mapped area can be shown to be associated with one of the following geological features: (i) basement highs, (ii) large bodies of heavy basic or ultrabasic rock in the crystalline basement, (iii) large igneous intrusions within the sedimentary column and thick accumulations of volcanic rocks or their associated eruption centers, (iv) major basement faults. Large-scale geological structures such as the Central, Viking and Sogn Grabens and the East Shetland, Stord, Forth Approaches and Norwegian-Danish Basins are essentially in isostatic equilibrium and are only locally marked by relatively weak gravity minima. A residual gravity anomaly map has been produced by subtracting from the observed Bouguer anomalies the estimated gravity effect of an assumed thinned crust. This residual gravity anomaly map shows a number of features of the Bouguer anomaly field with greater clarity.     

Crustal structure and tectonic setting of the south central Andes from gravimetric analysis

Based on terrestrial gravity data, in this paper we prepared a map of Bouguer anomalies, which was filtered to separate shallow and deep gravity sources. Based on a density model and gravimetric inversion techniques, the discontinuous crust-mantle boundary and the top of crystalline basement were modeled. Subsequently, the equivalent elastic thickness (Te) was evaluated, considering information from the crust-mantle discontinuity and topographic load, finding high Te values in the eastern Andean foothills and west of the Velasco range. These results are consistent with the positive isostatic and residual Bouguer anomaly values, which suggest the presence of high-density rocks in the mid-to upper crust. In addition, petrographic and geochemical analysis conducted in surface outcrops suggest a mantle origin.     

The geoid-to-quasigeoid difference using an arbitrary gravity reduction model

Recently it was proved that the classical formula for computing the geoid to quasigeoid separation (GQS) by the Bouguer gravity anomaly needs a topographic correction. Here we generalize the modelling of the GQS not only to Bouguer types of anomalies, but also to arbitrary reductions of topographic gravity. Of particular interest for practical applications should be isostatic and Helmert types of reductions, which provide smaller and smoother components, more suitable for interpolation and calculation, than the Bouguer reduction.     

巴颜喀拉地块东部及邻区重力均衡与岩石圈有效弹性厚度

基于SIO(Scripps Institute of Oceanography)最新全球重力和高程模型,计算了巴颜喀拉地块东部及邻区的布格重力异常、均衡重力异常、岩石圈有效弹性厚度及荷载比.结合大地热流、地震速度结构、地震活动和断裂构造分布等,分析了地壳均衡状态和岩石圈有效弹性厚度、地质构造单元间的差异及与地震活动的相关性特征.研究结果表明,该区域布格重力变化范围约为-500～0mGal(1mGal=10~(-5)m·s~(-2),下同),在巴颜喀拉块体东部区域形成弧形重力梯度带,近年来的中强地震活动频发于该梯度带不同部位,应与其应力依次释放有关;均衡重力异常结果表明,其变化范围约为-80～+100mGal,且大部分区域处于±20mGal以内的被认为处于重力均衡的状态,重力非均衡(正或负)多出现于块体边界带附近,地震多发生在靠近块体边界的均衡重力异常(正或负,主要为正)区域内;巴颜喀拉地块东部及邻区岩石圈有效弹性厚度(T_e)为10～65km,不同构造单元之间T_e空间分布差异明显,较低的T_e值出现在龙门山构造带附近,T_e值为20km左右,岩石圈荷载加载比为0.5～0.8,表明现今的岩石圈挠曲状态主要由莫霍面加载形成.进一步分析表明,巴颜喀拉地块东部挤压增生与横向流动同时发生,是造成该区域地震发生与重力均衡异常高值重合、岩石圈有效弹性厚度和大地热流值较低的主要原因.本文获得的地壳均衡特征及岩石圈有效弹性强度结果,加深了对巴颜喀拉东部及邻区岩石圈构造演化过程的认识.     

Crustal Structure of the Flood Basalt Province of Ethiopia from Constrained 3-D Gravity Inversion

The Oligocene Afar mantle plume resulted in the eruption of a large volume of basaltic magma, including major sequences of rhyolitic ignimbrites, in a short span of time across Ethiopia. In order to assess the impact of these magmatic processes on the crust and to investigate the general crustal configuration beneath the Ethiopian plateau, northern part of the Main Ethiopian Rift and the Afar depression, analysis and modeling of the gravity field have been conducted. The Bouguer gravity map is dominated by long-wavelength anomalies that primarily arise from the isostatic compensation of the topography. Consequently, anomalies within the crust/upper mantle are masked and quantitative interpretation becomes difficult. The long-wavelength anomalies are approximated using admittance technique and subsequently removed from the Bouguer anomalies to obtain the residual isostatic anomalies. The residual map contains both short- and intermediate-wavelength anomalies related to geologic and tectonic features. The long-wavelength regional isostatic field is used to map the crust-mantle interface and the results are in good agreement with those determined by other geophysical methods. Seismic constrained gravity inversion was performed on the isostatic residual field and series of three-dimensional models have been constructed for the structures of the crust and upper mantle beneath the uplifted and rifted flood basalt province of northern Ethiopia. The inversion results have shown that the NW plateau has thick crust that rests on normal lithospheric mantle. Afar, On the other hand, is marked by thin stretched crust resting on a low-density upper mantle indicating a hotter thermal regime and partial melt. No lithospheric mantle is observed beneath Afar. The models further indicate the presence of an extensive sub-crustal thick (~12 km on average) and high-density (~3.06 gm/cc) mafic accreted igneous layer of fractionated cumulate (magmatic underplating) beneath the NW plateau. The study suggests that the underplate was fundamental to the accretion process and may have played a role in compensating most of the plateau uplift and in localizing stresses.     

青藏高原东南缘多尺度重力场及构造含义

本文基于EGM2008重力场模型研究了青藏高原东南缘均衡重力异常和多尺度的布格重力异常特征,以鲁甸和景谷地震为例,认识其深部构造环境和动力学过程,为该区域的构造运动和地震孕育环境研究提供依据.结果表明,研究区布格重力异常和均衡重力异常与地质构造格局相关性较好,川滇地块剧烈的区域布格重力异常和非均衡状态与其强烈的地壳变形...     

新疆精河6.6级地震周边地区密度构造、均衡异常以及岩石圈挠曲机理

依据EIGEN-6C4重力模型和ETOPO1高程模型数据,围绕新疆精河6.6级地震展开岩石圈均衡与挠曲机理研究,得到如下结论：（1）震中附近的布格与自由空气重力异常分别为-221和-92mGal(10~(-5 )m·s~(-2)）,震中位于重力异常高梯度带上;（2）震中周边地区地壳厚度约为50km,密度结构总体变化平缓,东西方向地壳厚度变化较小,但自南向北地壳厚度逐渐变薄,精河6.6级地震初始破裂发生在上中地壳分界面附近;（3）震中附近岩石圈承载的垂向构造应力为20MPa左右,震中位于岩石圈垂向构造应力极大值附近的高梯度带上;（4）地震周边地区岩石圈有效弹性厚度最优解为26km,加载比最优解为F_1=1,F_2=F_3=0,表明该区域岩石圈相对坚硬,且导致岩石圈变形的初始加载全部来自地表．     

TOPOGRAPHIC AND ISOSTATIC CORRECTIONS TO GRAVITY SURVEYS IN MOUNTAINOUS AREAS*

In mountainous areas it is an inadequate procedure to reduce gravity observations by merely subtracting the effect of an infinite flat slab of material between the station and sea-level, and adding a terrain correction. A programme is described which directly computes the effect of masses above sea-level, and mass-deficiencies below it; as well as the effect of compensating masses under the Pratt and Airy isostatic systems. As an example, the method has been applied to a regional gravity survey of Papua and New Guinea where it is seen to remove the usually high correlation of the Bouguer anomaly with local topography.     

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

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

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

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

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.     

近地表密度估计的重力贝叶斯分析方法及在云南地区的应用

基于布格重力异常相对于地形起伏光滑分布的约束条件,从一维自由空气重力异常数据出发,采用贝叶斯方法估算近地表岩石密度,同时采用三次B样条函数拟合布格重力异常,获取光滑分布的布格重力异常.数据拟合和光滑约束之间的权重采用Akaike贝叶斯准则(ABIC准则)自动确定.均匀剖分模型和不均匀剖分模型数据试验都验证了该方法的有效性.相关参数评价表明,足够多的样条系数可以提高估计结果的准确性,样条系数的个数接近测点数时可获得较稳定的估计结果.增大异常的噪声水平时,ABIC准则可有效地自动增大先验光滑约束的权重.云南地区两条重力剖面应用结果表明,剖面沿线的近地表密度值起伏变化明显(达2.45~2.8g·cm^-3),前寒武纪和古生代地层密度相对较高(主要为2.53~2.75g·cm^-3),而中生代密度较低(2.45~2.73g·cm^-3);本文估计的近地表密度结果与区域物性资料及地表地质特征较吻合;估计的剖面布格重力异常具有光滑性;红河断裂两侧近地表密度差异较大,可达0.4g·cm^-3.本文获得的两条剖面近地表密度结构和布格重力异常为该区深部结构与构造研究提供更可靠的重力基础数据.     

格莱尼改正及其在地球物理数据处理中的意义

格莱尼改正是重力改正的一种,是对球面半径约166.7 km以外地形和均衡造成的重力综合效应的改正,改正后所得格莱尼异常可对应研究区范围内地壳结构特征,具有重要的研究意义.但近年来国内的重力研究中,对格莱尼改正的认识和重要性还不是十分透彻,特别是在大区域研究中,忽略这一远区重力效应可能会对研究结果的解释造成误导.针对这一问题,本文介绍了格莱尼改正的含义、由来与研究历史,指出了在大区域研究中需要考虑格莱尼效应的重要性,以“扇形球壳块”法为例给出了格莱尼改正的计算原理,列举了格莱尼异常在地球物理中的常见应用,并通过重力改正与莫霍面深度反演实例,通过对反演结果的控制点检测与精度评价,指出格莱尼异常反演所得莫霍面埋深更接近地震方法结果,从而表明了格莱尼改正在地球物理研究中的不可忽视性,并指出完全布格改正与格莱尼改正在计算方法上可实现统一.     

天山及邻区Vening Meinesz均衡重力异常特征及其动力学意义

本文基于Vening Meinesz区域均衡模型,通过试验不同参数计算Vening Meinesz均衡补偿深度,将其与CRUST1.0模型给出的莫霍面深度进行拟合,得到适应于天山及邻区的平均补偿深度、"地区性指标"以及区域补偿半径.结合地球重力场模型EIGEN-6C4与地形数据,利用球冠体积分方法进行地形效应、沉积层效应计算和均衡校正,得到了研究区的Vening Meinesz均衡重力异常.结果显示天山及邻区的均衡重力异常幅值在-110~120 mGal之间,表明了天山及周边盆地岩石圈所处于的均衡状态,同时揭示了研究区的壳幔密度分布特征.天山、塔里木盆地、准噶尔盆地等块体的地壳垂向形变可能部分地由均衡调整引起,且均衡调整趋势与地面形变测量结果相契合.通过对均衡重力异常成因的解释,从地壳均衡角度分析了该地区复杂的构造背景及其新生代以来的演化历程.     

ESTIMATION OF THE DEGREE OF POLYNOMIAL FITTED TO GRAVITY ANOMALIES AND ITS APPLICATION1

A method for estimating the degree of polynomial fitted to gravity anomalies to evaluate the regional anomaly is presented. The anomaly can be fitted by polynomials of different degrees with the least-squares method and the optimum degree of the polynomial evaluating this anomaly can be estimated from the point of discontinuity of the gradient on a graph of variance against the polynomial degree. The Bouguer gravity is initially separated by upward continuation to a proper height and then the degree of regional polynomial to fit the Bouguer anomaly can be estimated. Theoretical and field examples show the effectiveness of the method.     

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

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

华北北部均衡重力异常的初步研究

本文利用地面实测的重力资料和地形高程资料,采用爱里-海斯于宁(Airy-Heiskenen)均衡理论模式,构制了华北北部N37°—41°,E115°—120°范围内的均衡重力异常图和1°×1°均衡重力异常图。并对该图进行了分析、解释,同时结合本区地质构造、现代构造运动以及地震活动性特点分别进行了讨论     

3D upper crustal density structure of the Deccan Syneclise,Central India

A constrained 3D density model of the upper crust along a part of the Deccan Syneclise is carried out based on the complete Bouguer anomaly data. Spectral analysis of the complete Bouguer gravity anomaly map of the study region suggests two major sources: short wavelength anomalies (<100 km) caused primarily due to the density inhomogeneities at shallow crustal level and long wavelength anomalies (>100 km) produced due to the sources deeper than the upper crust. A residual map of the short wavelength anomalies is prepared from the complete Bouguer anomaly using Butterworth high‐pass filter (100 km cut‐off wavelength). Utilizing the constraints from deep resistivity sounding, magnetotellurics and deep seismic sounding studies, 2.5D density models have been generated along 39 profiles of this region. The mismatch between the calculated response of the a priori 2.5D model with the residual (short wavelength) gravity anomalies is minimized by introducing high‐density intrusive bodies (≥2.81 g/cm³) in the basement. With these 2.5D density models, the initial geometry of our 3D density model, which includes alluvium, Deccan trap, Mesozoic sediment and high‐density intrusive bodies in the basement up to a depth of 7 km (upper crust), is generated. In the final 3D model, Deccan trap extends from 200 m to nearly 1700 m below the 90–150 m thick Quaternary sediment. Further down, the sub‐trappean Mesozoic sediment is present at a depth range of 600–3000 m followed by the basement. The derived 3D density model also indicates six intrusive bodies of density 2.83 g/cm³ in the basement at an average depth of about 4–7 km that best fits the residual gravity anomaly of the study area.