中国布格重力异常水平梯度图的判读和构造解释

分析和对比中国大陆布格重力异常水平梯度的彩色图像,得到一些新的构造特征信息. 大致以105°E为界,图像的东部以蓝色低梯度为背景,出现中、高梯度带状构造和大梯度差值边界的线状构造;西部以高、中梯度的带状构造为格架,镶嵌着低梯度区,区内有阴影线状构造. 总体上,重力异常水平梯度构造可概括为三大组. 第一组是NNE向的大兴安岭、太行山、武陵山高梯度带与斜切东部全区的郯庐大断裂带和NNW向的斜切西部全区的三条大断裂带;第二组为西部的NWW向与NEE向高、中梯度带的交叉构造格局,以及东部近于等间距的三条NWW向断续的阴影构造带和一条单一的构造线,而华南则仅有NE向的低梯度构造线;第三组是横贯全区的的几条近EW向和中部SN向宽而复杂的构造带,在后者的东、西两侧散臣着大间距、小长度、断续的SN向阴影构造线. 这三组重力梯度构造不仅反映了地貌起伏的信息,而且显示了大区域构造,特别是深部构造的密度差. 各组线状构造带的内部结构和各组的形成时间及演变历史可能不同,但它们相互复合构成现今总体的构造格局,对大尺度现代构造应力场作用都会有不同的响应,因此导致主应力方向的区域性变化.     

Hidden equivalence and efficiency of gravity data interpretation

Analytical properties of the solutions of the inverse problem of gravimetry are studied in the context of the approximative approach and method of linear integral representations. A new effect of the inheritance of specific analytical properties in the solutions, which is named the effect of hidden equivalence, is revealed and investigated. This effect significantly influences the interpretation informativity in the complex models of a medium, and it should be taken into account in the geological interpretation of gravity data. Hidden equivalence is studied for both linear and nonlinear inverse problems. As an example of a nonlinear problem, the inverse problem of structural gravimetry is analyzed. The correlation between the boundary equations and boundary values of the harmonic functions is demonstrated. Methods are suggested to allow for the effects that occur by expanding the approximative approach for complex conditions (the models of the media with spatially distributed parameters) by the dynamical and criterion principles.     

Analysis and interpretation of gravity data from the Aluto-Langano geothermal field of Ethiopia

The Aluto-Langano geothermal field is located in the central southern portion of Ethiopia within the Ethiopian Rift Valley. The gravity of the area was surveyed in an attempt to delineate the subsurface structure and to better understand the relationship between the geothermal systems and the subsurface structure. The gravity data were analyzed using integrated gradient interpretation techniques, such as the Horizontal Gradient (HG), Source Edge Detection (SED), and Euler Deconvolution (ED) methods. These techniques detected many faults that were compared with the mapped faults in the surface geology. The results of the present study will lead to an improved understanding of the geothermal system in the study area and aid the future geothermal exploration of the area.     

离散时变重力数据的可视化、指标量定义与解释

针对流动重力观测获得的数据成果特点,提出一种基于流动重力段差变化的时变重力数据可视化方法,并定义了两个指标量G值和C值,用以评价区域性重力场变化的显著性程度.在此基础上,应用该方法处理和分析了首都圈地区的流动重力数据.结果表明:与传统采用等值线方式来刻画时变重力场特征不同,该方法更能突出发生变化的重力测点位置、测量误差及其变化的显著性程度,可为研究与地震孕育、发生过程有关的重力场变化提供更多的定量依据.      

Wavelet analysis and interpretation of gravity data in Sichuan-Yunnan region, China

Introduction Sichuan-Yunnan region, located in the east margin of Qinghai-Xizang (Tibetean) Plateau, is a transitional zone between the rapidly upheaving Tibetean Plateau and relatively steady Yangtze Platform. Under the pressure exerted by the northward movement of Indian Plate, Sichuan-Yunnan region has been undergone strong deformation and regmagenesis, becoming one of the regions with the most intensive seismicity in the world. The research on the tectonics and seismicity there is alw…     

Quantitative interpretation of three dimensional finite mass distribution from gravity data

Summary It is well known that the problem of determination of subsurface masses from gravity observations does not offer any unique solution. Yet, as has been shown in this paper, under some simplified but realistic assumptions some parameters of the subterrainean body are uniquely determinable in terms of some weighted double summations of gravity values (residuals). Physical significances of these parameters are also discussed. They themselves suggest about the possible mass distribution, more over they can be utilised in «cut and try method» of the determining masses.Published with the kind permission of the Director, Geological Survey of India, Calcutta.     

Morphometric analysis and tectonic interpretation of digital terrain data: a case study

Tectonic movement along faults is often re?ected by characteristic geomorphological features such as linear valleys, ridgelines and slope‐breaks, steep slopes of uniform aspect, regional anisotropy and tilt of terrain. Analysis of digital elevation models, by means of numerical geomorphology, provides a means of recognizing fractures and characterizing the tectonics of an area in a quantitative way. The objective of this study is to investigate the use of numerical geomorphometric methods for tectonic geomorphology through a case study. The methodology is based on general geomorphometry. In this study, the basic geometric attributes (elevation, slope, aspect and curvatures) are complemented with the automatic extraction of ridge and valley lines and surface speci?c points. Evans' univariate and bivariate methodology of general geomorphometry is extended with texture (spatial) analysis methods, such as trend, autocorrelation, spectral, and network analysis. Terrain modelling is implemented with the integrated use of: (1) numerical differential geometry; (2) digital drainage network analysis; (3) digital image processing; and (4) statistical and geostatistical analysis. Application of digital drainage network analysis is emphasized. A simple shear model with principal displacement zone with an NE–SW orientation can account for most of the the morphotectonic features found in the basin by geological and digital tectonic geomorphology analyses. Copyright © 2003 John Wiley & Sons, Ltd.     

应用多种来源重力异常编制中国海陆及邻区空间重力异常图及重力场解读

本文通过分析陆地实测空间重力异常数据、海洋船载测量空间重力异常数据、卫星测高重力异常,布格重力异常数据、EGM2008地球重力模型数据等多种来源数据的性质和精度,并对相关数据进行对比,研究了编制1:500万中国海陆空间重力异常图的数据使用方案和技术方法.在地形较为平坦、实测数据分布均匀的陆区,使用实测数据,在地形复杂,实测数据稀少以及没有实测数据的陆区或岛屿,利用布格重力异常反推空间异常的方法合成平均空间重力数据,西藏地区的数据对比实验证明合成平均空间重力异常数据是一种有效的数据补充.利用三观测列方差分解法在南海地区对船载测量空间重力数据和美国SS系列及丹麦DNSC08GRA卫星重力数据进行了方差分解计算,结果表明不同来源的卫星测高重力数据具有很大的一致性,数据精度较以往有了很大的提高.海区空间重力数据使用原则是在船载重力测量数据校准下,全面使用卫星测高重力数据进行编图.海陆过渡区的异常处理应以EGM2008地球重力模型重力场为基准参考场,实现海陆异常平缓过渡,无缝连接.对中国海陆空间重力异常场进行了小波变换处理,对空间重力异常场进行了解读,勾画出三横四竖的一级重力梯级带及其所围限的8个一级重力异常区,并划分了二级重力异常区和梯级带,为块体构造学体系中大地构造格架的建立提供了地球物理证据.     

The R-approximation method and its application to the interpretation of detailed gravity and magnetic data

The application of the R-approximation method based on the use of the Radon integral transform within the framework of the method of linear integral representations to the interpretation of detailed gravity and magnetic data is considered in this paper. The results are presented for some research areas located in Russia and in the Crimea (Ukraine). The method determines the distribution of anomalous field elements, construction of the analytical continuation of fields, and their linear transformations with a sufficiently high degree of accuracy.     

A generalized method for various data processing techniques in gravity interpretation

Summary The various data processing techniques (downward continuation, first and second derivatives and their downward continuation) used in gravity interpretation, are analogous to different types of linear filtering operations whose theoretical filter (amplitude) responses can be derived from by suitably choosingN andd, whereu andv are angular frequencies in two perpendicular directions,d the height or depth of continuation in unit of grid interval; andN denotes the order of the vertical derivative. By incorporating a mathematical smoothing function, ( being the smoothing parameter) in the theoretical filter response function, it has been possible, by selecting a suitable value of smoothing parameter, to establish an approximate equivalence of the effect of the mathematical smoothing with the inherent smoothing introduced, because of the numerical approximation (approximation error) for practically all data-processing techniques. This approximate equivalence leads to a generalized method of computing sets of weight coefficients for various data-processing techniques from filter response matching method. Several sets of weight coefficients thus have been computed with different smoothing parameter. The amplitude response curves of the various existing sets of weight coefficients have also been calculated for assessing the quality of the approximation in achieving the desired filtering operation.     

重磁异常解释断裂构造的处理方法及图示技术

利用重磁异常解释断裂是地质构造研究的主要手段之一.重磁异常解释断裂通常是在等值线图或剖面平面图上进行的,而重磁异常在进行等值线网格化成图时会造成微弱信息丢失,这些原因造成重磁异常解释断裂的多解性.本文对重磁异常数据在化极、曲面延拓处理的基础上,采用水平梯度法提取延拓曲面上的重磁异常梯度带,之后对断裂带进行窄化处理,通过图示技术将重磁异常数据转换成灰度值,图像的灰度值以变密度显示,形成彩色的变密度图像,这样就提高了数据图像识别断裂的视觉效果.该方法应用在鄂尔多斯盆地不同层次(时间序列)的断裂研究中,提取和识别重磁异常特征所反映的断裂信息,效果较好.     

基于流动重力资料的平原地区深部构造解释

以河北省流动重力测点覆盖区域内邢台地区4条断裂为例,将卫星布格重力异常与实测布格重力异常进行对比,发现流动重力资料可更为详细地反映测点周围构造的空间分布特征,为精细刻画华北平原地区深部构造背景提供可靠资料。采用向上延拓方法,对实测布格重力异常进行处理,获得不同上延高度的布格重力异常图像,结合当地地震地质与构造资料进行分析,并对流动重力资料在构造解释上的应用前景进行展望,认为该方法可提高现有资料的利用价值,适合在流动重力测量地区推广。     

利用重磁资料研究莺-琼盆地构造分界及其两侧断裂特征

莺-琼盆地(即莺歌海盆地和琼东南盆地)位于印支地块与华南大陆交汇处,该区域地质构造复杂,为南海西北走滑型和伸展型陆缘的交汇区,也是印澳-欧亚板块碰撞“挤出-逃逸构造区”和“古南海俯冲-拖曳构造区”.莺-琼盆地的北部以红河断裂带中1号断裂为界,前人对其认识比较统一,但莺-琼盆地南部由于地质和地球物理资料较少、研究程度低,以致其构造格局及分界位置不明确.通过对重、磁资料的处理,认为在莺-琼盆地分界的北段存在明显的重力异常和磁力异常梯级带,该梯级带与1号断裂位置相对应;在其南段存在明显的重力异常和磁力异常高值带,推断该高值带为一“中建凸起”,结合重、磁对应分析技术确定了该“中建凸起”的分布范围.采用重、磁异常归一化总水平导数垂向导数技术识别出控制莺-琼盆地分界两侧的断裂构造,在分界以西断裂走向主要为北北西向,分界以东主要为北东向.通过地质和地球物理综合研究认为莺-琼盆地的北段以1号断裂为界,南段以“中建凸起”为界.     

Application of the F-approximation to gravity data interpretation. II. The results of testing with the use of detailed gravity and magnetic surveys

The basic theoretical foundations of the spectral analysis of gravity and magnetic data based on the approximation approach have been considered and corresponding expressions implementing various statements of the F-approximation method have been derived [Strakhov and Kerimov, 1999; 2000; 2001]. The author considered the algorithms and computer implementation of the F-approximation of anomalous gravity and magnetic fields presented in the results of testing with the use of a number of model examples and data on the accuracy of the model field approximation earlier on in [Kerimov, 2003]. In this paper, I consider methodological features of constructing the F-approximation from the detailed gravity and magnetic surveys, anomalous gravity and magnetic fields reconstructed with the use of the F -approximation, and examples of the field transforms obtained with the use of this approach.     

重力与地形数据揭示的巴颜喀拉块体东缘垂向构造应力场

利用基于消去-恢复原理的最小二乘配置方法,对2009—2013年相对重力/GPS联合观测数据与EGM2008模型数据进行融合,更新了巴颜喀拉块体东缘地区的自由空气与布格重力异常场.基于该布格重力异常数据,以CRUST1.0地壳密度模型为初始条件,使用二维多边形棱柱体正演与非线性最小二乘反演方法,获取了巴颜喀拉块体东缘地壳分层密度结构.基于地壳不可压缩和均衡调整原理提出了计算垂向构造应力新方法,并结合上述地壳分层密度结构和地形数据计算了巴颜喀拉块体东缘垂向构造应力分布.结果表明,龙门山断裂带中南段蓄积了较高的正向构造应力（约40 MPa）,马尔康周边地区蓄积了较高的负向构造应力（约-30 MPa）.对研究区域1970年以来5级（M_S）以上地震进行统计发现,地震多发生在垂向构造应力梯度带上,垂向构造应力为正的地区易触发浅源地震,为负的地区易触发深源地震.在地壳横向变形强烈的区域,垂向构造应力与地震深度的对应关系减弱.

     

Future of the interpretation theory of gravity and magnetic anomalies

The general methodology of the interpolation theory of potential fields (gravity and magnetic anomalies), and, primarily, the three fundamental ideas underlying this methodology—(I) approximation, (II) linearization, and (III) optimization-are the matter of the utmost importance. The rationale of these ideas and the ways of their practical implementation are discussed in this paper.     

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

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

Application of satellite imagery in tectonic interpretation

The advabtages of remote sensing are widely used in geologic research. In Hungary, however, in the deep sedimentary basin and densely covered mountain areas, the methods which have been developed and which have led to significant result in arid and semiarid regions cannot be used effectively.In this case, our interest should be focused on detecting structural features instead of on lithologic discrimination. Consequently, this approach can be applied for only morphologically performed areas. A better understanding of the tectonics can be achieved by the evaluation of its morphological expressions.A block-shift method is described for the structural evaluation of linear features, and it is applied for one of the most controversial and least known tectonic units: the Mecsek Mts. (S-Hungary). This method is based on the calculation of several conditional frequencies of length, azimuth and shape for overlapping blocks.Our experience leads to the conclusion that the widely accepted East-West segmentation of the Mecsek Mts. should be changed. We have found an Eastern-Central-Western Mecsek structure which correlates well with some up-to-date geologic mapping results. The lithologic differences are enhanced by morphological and structural features.Besides its scientific achievements, the economic importance of this research is supported by the perspective of the utilization of non-renewable natural resources in the Mecsek Mts.     

Application of the Hilbert transform for gravity and magnetic interpretation

Either the magnetic anomaly or the horizontal gradient of gravity, if plotted against its Hilbert transform and then joined by a smooth curve, the relation figure is obtained. Properties of relation figures for some two-dimensional sources of simple geometry are shown to be useful in distinguishing the source geometry and also in estimating the source parameters, even in the presence of remanence.