Introduction to ground surface self-potential tomography

A new approach to self-potential (SP) data interpretation for the recognition of a buried causative SP source system is presented. The general model considered is characterized by the presence of primary electric sources or sinks, located within any complex resistivity structure with a flat air-earth boundary. First, using physical considerations of the nature of the electric potential generated by any arbitrary distribution of primary source charges and the related secondary induced charges over the buried resistivity discontinuity planes, a general formula is derived for the potential and the electric field component along any fixed direction on the ground surface. The total effect is written as a sum of elementary contributions, all of the same simple mathematical form. It is then demonstrated that the total electric power associated with the standing natural electric field component can be written in the space domain as a sum of cross-correlation integrals between the observed component of the total electric field and the component of the field due to each single constitutive elementary charge. By means of the cross-correlation bounding inequality, the concept of a scanning function is introduced as the key to the new interpretation procedure. In the space domain, the scanning function is the unit strength electric field component generated by an elementary positive charge. Next, the concept of charge occurrence probability is introduced as a suitable function for the tomographic imaging of the charge distribution geometry underground. This function is defined as the cross-correlation product of the total observed electric field component and the scanning function, divided by the square root of the product of the respective variances. Using this physical scheme, the tomographic procedure is described. It consists of scanning the section, through any SP survey profile, by the unit strength elementary charge, which is given a regular grid of space coordinates within the section, at each point of which the charge occurrence probability function is calculated. The complete set of calculated grid values can be used to draw contour lines in order to single out the zones of highest probability of concentrations of polarized, primary and secondary electric charges. An extension to the wavenumber domain and to three-dimensional tomography is also presented and discussed. A few simple synthetic examples are given to demonstrate the resolution power of the new SP inversion procedure.     

据断层附近布格重力异常分布反推地质时期断层运动的探讨

本文应用变形介质位移场产生的附加重力场的积分表达式,对麦克斯威尔介质半空间中断层长期活动引起的永久变形与现今布格重力异常分布的关系作了定量分析,并提出了由断层附近布格异常分布特征反推地质时期断层运动趋势的一般方法。其主要应用于平原地区断层运动(特别是扭动方式)的判别。文中以滦南断裂为实例,通过数值积分,计算了几种断层运动可能产生的地表重力异常。经理论计算与实测重力资料相对比,该断层曾发生过下盘上升的正倾滑运动,并具有顺扭走滑运动分量     

磁总场异常垂直梯度三维相关成像

本文提出了磁总场异常垂直梯度三维相关成像方法,用于成像地下等效磁源分布.它首先将地下待成像空间划分为三维规则网格,然后直接计算每个网格节点磁偶板子在观测面理论磁总场异常垂直梯度与实测磁总场异常垂直梯度的互相关,其相关系数值表征等效磁偶极子分布(即磁偶极子发生的概率).理论长方体组合模型数据和实际某矿区磁测资料试验结果表明本文方法计算得到的相关系数值能基本反映地下的磁源分布,且分辨率明显高于磁总场异常三维相关成像的分辨率,也高于基于熵滤波分离异常的磁总场异常三维相关成像的分辨率.     

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

基于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.     

Role of Geological and Geophysical Data in Modeling a Southwestern Alluvial Basin

A parameter-estimation technique based on existing hydrological, geophysical, and geological data was developed to approximate transmissivity values for use in a ground-water flow model of the Animas Valley, southwest New Mexico. Complete Bouguer gravity anomaly maps together with seismic-refraction profiles, geologic maps, geologic, geophysical, and drillers' logs, water levels, and pumping-test data provide insight into the transmissivity of bolson deposits throughout the basin. The transmissivity distribution was primarily based on reported pumping and specific-capacity tests in conjunction with complete Bouguer gravity anomaly maps and well log data. Reported transmissivity values were characterized by gravity values and well log data. In grid blocks lacking pumping and specific-capacity tests, transmissivity values were assigned based on the relationship of gravity values and well log data within the grid block to gravity values and well log data within other grid blocks for which transmissivity values are available. A two-dimensional, finite-difference, ground-water flow computer code was used to evaluate the effectiveness of the parameter-estimation technique. Although the trial-and-error method of calibration was employed, the actual computer time necessary for model calibration was minimal. The conceptually straightforward approach for parameter estimation utilizing existing hydrological, geophysical, and geological data provides realistic parameter estimates.     

ANALYSIS ON STRUCTURAL FEATURES:INSIGHTS FROM SATELLITE-DERIVED GRAVITY MODEL AND SEISMIC PATTERN

In this paper, based on a large number of cumulative observational data from the seismic monitoring network in China, we grid the research area to calculate the density values at each grid node and convert the qualitative earthquake epicenter distribution to quantitative seismic pattern. Minimum magnitude of completeness(M_C)is determined by magnitude-rank analysis, which provides lower limit earthquake and original time. New satellite-derived gravity model v23.1, which is based on satellites CryoSat-2 and Jason-1 data, is used to determine the Bouguer gravity anomaly derived from free-air gravity anomaly and elevation database sets SRTM30, and ultimately, the complete Bouguer correction is obtained. In this paper, the Xingtai earthquake zone and Tanlu fault zone (Anhui segment) are selected for case study. Bouguer gravity anomaly presents a NE-trending U-shaped narrow strip in the Xingtai earthquake zone, and its location is consistent with Shulu Fault Basin. Grid density value contours are restricted by the U-shaped strip, and the extreme value of seismic activity density lies in the bottom of the U-shaped strip as shown in the cross section. The results of Bouguer gravity anomaly and upward continuations to the different heights show good linearity and gradient in the Tanlu fault zone (Anhui segment); and both long-axis direction of seismic pattern and nodal plane strike of seismogenic fault from focal mechanism solutions trend NNE. In short, the Tanlu fault zone(Anhui segment)is a large deep-seated fault that still has the ability to control seismic activity along it. Based on the measured gravity and magmatic data, using the edge detection TDX method to interpret the concealed boundary of the Anqing M4.8 earthquake near the Tanlu fault, and combining with the results from deep seismic reflection profiles of the study area, we discussed the causative fault of the Anqing earthquake.     

宜川—泰安重力剖面布格重力异常及归一化总梯度的构造意义

利用宜川—泰安剖面高精度重力观测及GPS测量结果,计算得到该剖面较为精细的垂向一阶导数梯度和布格重力异常归一化总梯度图像.结合前人相关研究结果,本文认为垂向一阶导数梯度主要反映浅部异常场源的位置和低密度地质体的特性.利用该梯度图像,对宜川—泰安剖面经过的7条主要断裂及其空间分布特征进行分析解释,并推测各断裂的可能埋深和性质.      

基于抛物线密度模型的频率域三维界面反演及其在川滇地区的应用

密度界面反演作为了解地球内部结构的一种重要方法,长期以来都是重力学研究的主要内容.本文结合抛物线密度模型及频率域算法的优点,将抛物线密度函数应用于Parker-Oldenburg算法,经过理论推导得到了抛物线密度模型的频率域公式,从而建立了基于抛物线密度模型的三维密度界面重力异常正反演的算法和流程.理论模型数据试验表明本方法快速、有效,适用于大多数浅部比深部增加更快的实际地壳密度.研究中还利用该方法对川滇地区重力异常进行了反演,获得了该区的莫霍面深度分布,并与接收函数研究结果进行对比分析,进一步验证了本文方法的正确性和有效性.     

临汾强震区地壳结构及发震背景研究

以诸城—宜川深地震测深速度剖面为约束, 对沿该剖面得到的高精度重力数据进行拟合, 并对临汾强震区平面布格重力异常进行处理, 得到了该地区地壳密度结构及平面重力异常分布． 利用上述结果分析了临汾强震区的地壳结构及构造环境, 结合前人相关研究成果, 认为临汾强震区地壳中存在塑性相对较强的介质, 洪洞和临汾两次历史地震皆发生在其与周边弹性介质的转换边界上． 另外, 临汾凹陷南北两侧局部构造环境存在差异, 在区域应力场作用下, 导致了洪洞地震和临汾地震的发生. 两次地震在发震时间、 地点和震级等地震要素上有所不同．      

Gravity contributions to the understanding of salt tectonics from the Jebel Cheid area (dome zone, Northern Tunisia)

Detailed gravity measurements integrated with geological data were computed to constrain the mechanisms that were active during the emplacement of the Triassic evaporite-bearing folds of Jebel Cheid from the salt-dome zone in the Atlassic region. The gravity analysis consists in mapping the contrasting gravity responses: complete Bouguer anomaly, residual anomaly and derivative maps; the main results display a positive amplitude gravity anomaly as the response of Triassic evaporite bodies and important NE–SW-trending features at the boundaries between the Triassic outcrops and their enveloping strata. In contrast with gravity calculations of a salt dome structure usually resulting in negative gravity anomaly models, the Jebel Cheid clearly expresses a positive gravity anomaly; furthermore, this result is supported by synthetic gravity interpretation.     

大别造山带东段重力异常多尺度分解及其构造意义

文中通过多源数据融合、模型构建、数据试验、二维离散小波变换和功率谱分析等方法获取了大别造山带东段深、浅部场源布格异常及其场源似深度,并结合地壳结构、地质构造、岩石圈有效弹性厚度和地震活动等资料,讨论了地壳深、浅部的结构特征及地震活动构造背景。结果表明,低频布格异常显示大别造山带东段与华北地块间深部构造缝合带在东部应位于青山-晓天断裂前缘,在落儿岭-土地岭断裂和商城-麻城断裂之间向N偏移至梅山-龙河口断裂之下,造山带南侧与扬子地块间深部构造缝合带位于襄樊-广济断裂以北约20km,造山带东侧与扬子地块间的深部构造转换带位于郯庐断裂带之下,造山带东段腹地显著的低频布格异常低值表明对应部位的莫霍面存在明显下凹,造山带内部的布格异常高梯度带表明其深部结构不完整;高频布格异常揭示肥中断裂、六安-合肥断裂、肥西-韩摆渡断裂和郯庐断裂带等主要断裂对地壳中上部密度结构的影响明显,落儿岭-土地岭断裂对地壳中上部密度结构的影响范围向N延伸至肥西-韩摆渡断裂前缘。结合地震活动资料进一步分析认为,大别造山带东段与华北地块在青山-晓天断裂前缘附近接触和相互作用,且大别造山带东段地壳深、浅部结构均不完整,不利于应力积累,趋向于在断裂交错的脆弱部位频繁释放应力,是霍山地区小地震活动频繁的主要原因。     

Focusing iterative migration of gravity gradiometry data acquired in the Nordkapp Basin,Barents Sea

Geological interpretation based on gravity gradiometry data constitutes a very challenging problem. Rigorous 3D inversion is the main technique used in quantitative interpretation of the gravity gradiometry data. An alternative approach to the quantitative analysis of the gravity gradiometry data is based on 3D smooth potential field migration. This rapid imaging approach, however, has the shortcomings of providing smooth images since it is based on direct integral transformation of the observed gravity tensor data. Another limitation of migration transformation is related to the fact that, in a general case, the gravity data generated by the migration image do not fit the observed data well. In this paper, we describe a new approach to rapid imaging that allows us to produce the density distribution which adequately describes the observed data and, at the same time, images the structures with anomalous densities having sharp boundaries. This approach is based on the basic theory of potential field migration with a focusing stabilizer in the framework of regularized scheme, which iteratively transfers the observed gravity tensor field into an image of the density distribution in the subsurface formations. The results of gravity migration can also be considered as an a priori model for conventional inversion subsequently. We demonstrate the practical application of migration imaging using both synthetic and real gravity gradiometry data sets acquired for the Nordkapp Basin in the Barents Sea.     

Two examples of spectral approach to source depth estimation in gravity and magnetics

Fourier transformation of gravity and magnetic anomalies from space to the frequency domain provides aready method for source depth estimation since progressively deeper sources are indicated by the lower frequency components of an anomaly. Two examples are presented here to demonstrate the efficacy of this approach. One example pertains to the interpretation of vertical intensity magnetic anomalies over a layered ultramafic body for estimating the thickness of relatively non-magnetic layer in it and map the dispositions of the layers. The other example is the estimation of crustal thickness in the Indian region from the Bouguer anomaly map. The spectral method is shown, by these two examples, to provide a rapid and elegant tool for the source depth estimation for magnetic and gravity data.     

东北地区构造分区与深断裂研究——基于重力场小波多尺度分解

基于东北地区的布格重力异常资料,通过小波多尺度分解,将重力异常分解为四阶.提取主要反映深断裂的重力异常2、3阶小波细节,利用不同方向的水平梯度,计算并强化深断裂的线性特征,解译出12条深断裂,同时,结合功率谱法、东北地区地震分布、震源机制解和地质信息推断断裂的下切深度和性质,结果表明:12条深断裂分别对应拼贴地块的缝合带和对东北地区盆地群的分布与演化有着明显控制作用的岩石圈断裂两类;深断裂体系主要方向为北东、北北东向且大多为走滑断层;提取主要反映深层块体异常的3、4阶小波细节并结合深断裂体系将研究区划分为五个重力场一级分区,并首次将额尔古纳块体和兴安块体分开.五个大区分别为额尔古纳块体、兴安块体、松嫩块体、佳木斯块体和华北板块;其中兴安块体又分为二连盆地、大兴安岭两个二级分区;松嫩块体又分为小兴安岭、松辽盆地、张广才岭三个二级分区;佳木斯块体又分为那丹哈达地体、三江盆地、桦南—穆棱三个二级分区,华北板块又分为辽东半岛、华北板块、北部断裂带三个二级分区,结果表明了研究区的基本构造格局.     

THE DEPTH OF MOHO INTERFACE AND CRUSTAL THICKNESS IN SICHUAN-YUNNAN REGION,CHINA

By using moving average method to separate Bouguer gravity anomaly field in Sichuan-Yunnan region, we got the low-frequency Bouguer gravity anomaly field which reflects the undulating of Moho interface. The initial model is obtained after seismic model transformation and elevation correction. Then, we used Parker method to invert the low-frequency Bouguer gravity anomaly field to obtain the depth of Moho interface and crustal thickness in the area. The results show that the Qinghai-Tibet block in the northwest of the study area deepens and thickens from the edge to the interior, with the depth of Moho interface and the crust thickness of about 52~62km and 54~66km, respectively. The depth of Moho interface in Sichuan Basin is about 38~42km. In Sichuan-Yunnan block, the depth of Moho interface is about 42~62km from southeast to northwest. Beneath the West Yunnan block, west of the Red River fault zone, the Moho depth is about 34~52km from south to north. The Longmen Mountains and Red River fault zone are the gradient zone of the Moho depth change. Along the Red River fault zone, the depth difference of Moho interface is increasing gradually from north to south. No obvious uplift is found on the Moho interface of Panzhihua rift valley. The depth of Moho interface distribution in Sichuan and Yunnan is obviously restricted by the collision between the Indian plate and the Eurasian plate and the lateral subduction of the Indo-China peninsula. The mean square error of the depth of Moho interface is less than 1.7km between the result of divisional density interface inversion and artificial seismic exploration. At the same time, we compared the integral with divisional inversion result. It shows that:in areas where there is obvious difference between the crust velocity and density structure in different tectonic blocks, the use of high resolution seismic exploration data as the constraints to the divisional density interface inversion can effectively improve the reliability of inversion results.     

TWO STRUCTURAL MODELS FOR THE WESTERN FLANK OF THE BRABANT MASSIF*

This paper is concerned with an interpretation of the Bouguer anomaly on the Western flank of the Brabant Massif (Belgium). The position, shape, and density contrast of elementary bodies in the upper part of the earth's crust are determined in a purely numerical manner. A batholith-like body and a basin-type structure both adequately account for the observed surface field, at least from a one-sided geophysical and mathematical point of view. Since seismological control is lacking the two models are representative of the ambiguity and indeterminacy of the definitions of the geologic cause of the anomalous features, if the data are restricted to gravity information alone.     

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

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

应用卫星重力信息对横断山系地区布格重力异常特异分布的纠正

在云南省西部,跨越中、缅两国交界的横断山系地区（97°E~102°E,24°N~30°N）有近一半的面积尚没有重力测点、即重力数据空白区和重力测点稀少的普查级测区.以前的有关文献、图集中所给出对此地区的重力场都是十分模糊的结果与图件.因此应用这些资料无法详细地研究该地区重力场特征与深部地壳结构（构造）.本文应用卫星重力异常资料作为“近似空间重力异常”,经中间层改正后给出的“计算布格重力异常”,其分布特征与该地区的地形高程呈很好的镜像相关.对相应山脉、河谷以及断裂构造都有所反映.特别是在横断山系地区该布格重力异常呈现为近南北的走向.为此,据该“计算布格重力异常”,并选定对该区有代表性的一条重力异常剖面作正反演计算,以得到其地壳深部结构剖面.结果表明,在横断山脉地区的地壳厚度在51~56 km间起伏变化;滇西北云岭山系以及玉龙山区的地壳厚度约在60 km以上. 最后,对所得结果与图件进行了讨论,并提出了几点认识和纠正的建议.     

Comparison of different gravity field implied density models of the topography

Density within the Earth crust varies between 1.0 and 3.0 g/cm³. The Bouguer gravity field measured in south Iran is analyzed using four different regional-residual separation techniques to obtain a residual map of the gravity field suitable for density modeling of topography. A density model of topography with radial and lateral distribution of density is required for an accurate determination of the geoid, e.g., in the Stokes-Helmert approach. The apparent density mapping technique is used to convert the four residual Bouguer anomaly fields into the corresponding four gravity im-plied subsurface density (GRADEN) models. Although all four density models showed good correlation with the geological density (GEODEN) model of the region, the GRADEN models obtained by high-pass filter-ing and GGM high-pass filtering show better numerical correlation with GEODEN model than the other models.