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

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

Three-dimensional analysis of gravity anomalies of sedimentary basins by polygonal prismatic model with a quadratic density function

The decrease of density contrast in sedimentary strata may be approximated by a quadratic function. A sedimentary basin may be considered as a number of horizontal polygonal prisms of finite thickness placed one over the other. Equations for the gravity anomaly of ann-sided polygon prism have been derived using a quadratic density function. A method has been developed for inversion of the gravity anomalies using a polygonal prismatic model. An analytical method has been employed for evaluating the necessary derivatives as the computing time by this method is much slower than by a numerical method. Approximate equations have been derived for rapid calculation of the anomalies and derivatives. Efficient computer programs have been developed for calculation of the anomalies and derivatives by an appropriate use of the exact and approximate equations. The depths to the polygon prisms constituting the basin are iteratively adjusted by minimizing the sum of the squares of the differences between the observed and calculated anomalies. These methods have been applied to analyze the Bouguer anomaly map of the Los Angeles basin, California.     

基于等效源的总强度磁异常非线性处理方法

总强度磁异常(ΔT1)常规处理方法通常将其近似当作磁异常矢量在地磁正常场方向上的投影(ΔT2).然而对于高磁环境如磁铁矿处的磁异常场幅值可达104 nT甚至更大的情况,上述近似假设则不再成立,如采用常规处理方法可能会带来明显的误差.本文针对该问题,提出一种基于等效源的总强度磁异常非线性处理方法;该方法根据总强度磁异常获取流程,直接反演实测地磁场总强度幅值与正常场强度幅值之差来求取等效场源.本文首先通过模型试验,分析ΔT1与ΔT2的差异;然后将ΔT1当作ΔT2采用常规处理方法以及ΔT1采用新方法得到的结果作对比分析,结果表明新方法处理结果与理论值的差异为常规方法处理结果的1/5甚至更小,充分说明了新方法的有效性;最后将该方法应用于铁矿实测总强度磁异常处理实例中来转换计算磁异常总模量,其实际应用效果进一步体现了高幅值总强度磁异常数据处理过程中采用新方法的必要性.     

Estimating free parameters in 2D inversion: example of gravity inversion in a rifted basement

We advance a principle directed to assigning numerical values to free parameters usually present in inversion methods. It may be formulated as: ‘Optimum estimates of free parameters in an inversion procedure must lead, in tests using synthetic data, to solutions whose geometrical expression reflects the main qualitative or semiquantitative geological characteristic of the study area.’ To this end, the interpreter should (i) specify a typical anomalous source geometry which incorporates the most relevant geological information for the study area, (ii) compute the corresponding gravity anomaly and (iii) invert the anomaly for the source geometry finding the numerical values of the free parameters that lead to a solution closest to the typical source. Application of the above methodology to synthetic and real data from the basement relief of a rift basin has asserted its efficacy.     

Three-dimensional structure of the Laguna Salada Basin and its thermal regime

A comprehensive reinterpretation of the available gravity, magnetic, geothermal, geological and borehole information has been made of the Laguna Salada Basin to establish a 3D model of the basement and sedimentary infill. According to statistical spectral analysis, the residual gravity anomaly is due to sources with a mean regional depth of 2.8 km. The topography of the basement was obtained from a three‐dimensional inversion carried out in the wavenumber domain using an iterative scheme. The maximum density contrast of ?300 kg/m³ estimated from previous studies and the mean depth of 2.5 km finally constrained this inversion. The resulting model indicated that the sedimentary infill is up to 4.2 km thick at its deepest point. According to the gravity‐derived basement topography, the basin presents an asymmetry (i.e. it is of the half‐graben type). It is deeper to the east, where it is delimited from the Sierra Cucapah by a step fault. By contrast, the limit with the Sierra de Juarez is a gently sloping fault (i.e. a listric fault). The basement is not even, but it comprises a series of structural highs and lows. N–S to NW–SE and E–W to NE–SW faults delimit these structural units. The magnetic modelling was constrained by (i) the gravity‐derived basement topography; (ii) a Curie isotherm assumed to be between 7 km and 10 km; (iii) assuming induced magnetization only; (iv) the available geological and borehole information. The magnetic anomalies were interpreted successfully using the gravity‐derived basement/sedimentary interface as the top of the magnetic bodies (i.e. the magnetic modelling supports the gravity basement topography). An elongated N–S to NW–SE trending highly magnetized body running from south to north along the basin is observed to the west of the basin. This magnetic anomaly has no gravity signature. Such a feature can be interpreted as an intrusive body emplaced along a fault running through the Laguna Salada Basin. Treatment of the gravity and magnetic information (and of their horizontal gradients) with satellite image processing techniques highlighted lineaments on the basement gravity topography correlating with mapped faults. Based on all this information, we derived detailed geological models along four selected profiles to simulate numerically the heat and fluid flow in the basin. We used a finite‐difference scheme to solve the coupled Darcy and Fourier differential equations. According to our results, we have fluid flow in the sedimentary layers and a redistribution of heat flow from the basin axis toward its rims (Sierra de Juárez and Sierra Cucapah). Our model temperatures agree within an error of 4% with the observed temperature profiles measured at boreholes. Our heat‐flow determinations agree within an error of ±15% with extrapolated observations. The numerical and chemical analyses support the hypothesis of fluid circulation between the clay–lutite layer and the fractured granitic basement. Thermal modelling shows low heat‐flow values along the Laguna Salada Basin. Deep fluid circulation patterns were observed that redistribute such flow at depth. Two patterns were distinguished. One displays the heat flow increasing from the basin axis towards its borders (temperature increase of 20°C). The second pattern shows an increasing heat flow from south to north of the basin. Such behaviour is confirmed by the temperature measurements in the thermometric boreholes.     

南海及邻区莫霍面深度分布特征及其与各构造单元的关系

为了研究南海及邻区莫霍面分布特征及其与边缘海盆、海沟、岛弧、新生代沉积盆地的关系等构造单元的关系,本文通过对研究区的空间重力异常数据进行全布格改正,得到研究区内的布格重力异常,并以近年来的声纳浮标探测与海底地震仪探测剖面所得到的莫霍面深度资料为控制点采用三维带控制点界面反演方法得到了研究区的莫霍面深度图和地壳厚度图.本...     

Interpretation of magnetic anomalies by horizontal and vertical derivatives of the analytic signal

Magnetic anomalies are often disturbed by the magnetization direction, so we can’t directly use the original magnetic anomaly to estimate the exact location and geometry of the source. The 2D analytic signal is insensitive to magnetization direction. In this paper, we present an automatic method based on the analytic signal horizontal and vertical derivatives to interpret the magnetic anomaly. We derive a linear equation using the analytic signal properties and we obtain the 2D magnetic body location parameters without giving a priori information. Then we compute the source structural index (expressing the geometry) by the estimated location parameters. The proposed method is demonstrated on synthetic magnetic anomalies with noise. For different models, the proposed technique can both successfully estimate the location parameters and the structural index of the sources and is insensitive to noise. Lastly, we apply it to real magnetic anomalies from China and obtain the distribution of unexploited iron ore. The inversion results are consistent with the parameters of known ore bodies.     

Determination of Sedimentary Basin Basement Depth: A Space Domain Based Gravity Inversion using Exponential Density Function

An automatic inversion using ridge regression algorithm is developed in the space domain to analyze the gravity anomalies of sedimentary basins, among which the density contrast decreases with depth following a prescribed exponential function. A stack of vertical prisms having equal widths, whose depths become the unknown parameters to be estimated, describes the geometry of a sedimentary basin above the basement complex. Because no closed form analytical equation can be derivable in the space domain using the exponential density-depth function, a combination of analytical and numerical approaches is used to realize forward gravity modeling. The depth estimates of sediment-basement interface are initiated and subsequently improved iteratively by minimizing the objective function between the observed and modeled gravity anomalies within the specified convergence criteria. Two gravity anomaly profiles, one synthetic and a real, are interpreted using the proposed technique to demonstrate its applicability.     

Gravity inversion of 2D bedrock topography for heterogeneous sedimentary basins based on line integral and maximum difference reduction methods

Based on the line integral (LI) and maximum difference reduction (MDR) methods, an automated iterative forward modelling scheme (LI‐MDR algorithm) is developed for the inversion of 2D bedrock topography from a gravity anomaly profile for heterogeneous sedimentary basins. The unknown basin topography can be smooth as for intracratonic basins or discontinuous as for rift and strike‐slip basins. In case studies using synthetic data, the new algorithm can invert the sedimentary basins bedrock depth within a mean accuracy better than 5% when the gravity anomaly data have an accuracy of better than 0.5 mGal. The main characteristics of the inversion algorithm include: (1) the density contrast of sedimentary basins can be constant or vary horizontally and/or vertically in a very broad but a priori known manner; (2) three inputs are required: the measured gravity anomaly, accuracy level and the density contrast function, (3) the simplification that each gravity station has only one bedrock depth leads to an approach to perform rapid inversions using the forward modelling calculated by LI. The inversion process stops when the residual anomalies (the observed minus the calculated) falls within an ‘error envelope’ whose amplitude is the input accuracy level. The inversion algorithm offers in many cases the possibility of performing an agile 2D gravity inversion on basins with heterogeneous sediments. Both smooth and discontinuous bedrock topography with steep spatial gradients can be well recovered. Limitations include: (1) for each station position, there is only one corresponding point vertically down at the basement; and (2) the largest error in inverting bedrock topography occurs at the deepest points.     

THE MAGNETIC ANOMALY OF 3D SOURCES HAVING ARBITRARY GEOMETRY AND MAGNETIZATION

Magnetic anomalies of complicated 3D sources can be calculated by using a combination of analytical and numerical integration. Two surfaces and the magnetization parameters (the amplitudes of the induced and remanent components and the direction cosines) of the source can be defined by arbitrary functions or by discrete data points in a plane. When combined with a polynomial magnetization function in the direction of the third axis, 3D magnetization distribution can also be modelled. The method gives very general equations for anomaly calculation. It can be used for direct modelling of sources interpreted by seismic or other methods and also for interactive interpretation with fast computers. It is possible to calculate anomalies of, for example, intrusives or folded sedimentary beds whose surfaces are functions of horizontal coordinates and which have polynomial magnetization variations in the vertical direction due to gravitational differentiation and arbitrarily varying magnetization in the horizontal direction due to regional metamorphosis. If the distribution of magnetization parameters in the vertical direction cannot be described satisfactorily by polynomials, models can be used whose surfaces are functions of the vertical coordinate and which can then have any arbitrary magnetization distribution in the vertical direction.     

Inversion of magnetic measurements of the CHAMP satellite over the Pannonian Basin

The Pannonian Basin is a deep intra-continental basin that formed as part of the Alpine orogeny. In order to study the nature of the crustal basement we used the long-wavelength magnetic anomalies acquired by the CHAMP satellite. The anomalies were distributed in a spherical shell, some 107,927 data recorded between January 1 and December 31 of 2008. They covered the Pannonian Basin and its vicinity. These anomaly data were interpolated into a spherical grid of 0.5° × 0.5°, at the elevation of 324 km by the Gaussian weight function. The vertical gradient of these total magnetic anomalies was also computed and mapped to the surface of a sphere at 324 km elevation. The former spherical anomaly data at 425 km altitude continued downward to 324 km. To interpret these data at the elevation of 324 km we used an inversion method. A polygonal prism forward model was used for the inversion. The minimum problem was solved numerically by the Simplex and Simulated annealing methods; a L₂ norm in the case of Gaussian distribution parameters and a L₁ norm was used in the case of Laplace distribution parameters. We interpret that the magnetic anomaly was produced by several sources and the effect of the sable magnetization of the exsolution of hemo-ilmenite minerals in the upper crustal metamorphic rocks.     

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.     

磁化强度矢量反演新方法及其在卡拉塔格铜多金属矿床找矿中的应用

磁异常的反演是地球物理勘探的重要手段,三维磁化率反演是磁异常定量解释中的一种重要方法.由于剩磁的存在使得磁化方向与地磁场方向产生偏差,从而影响了磁异常反演与解释的精度.本文基于磁异常模量反演和磁化强度矢量反演方法得到了一种新的磁化强度矢量反演方法.与以往的磁化强度矢量反演方法相比,该方法以磁异常模量反演得到的磁化率模型为约束,采用Lp范数正则化方法求解,提高了磁化强度矢量反演的精度和效率.本文通过模拟试验的反演计算,验证了这种磁化强度矢量反演方法的有效性.最后,将本文方法应用于新疆东天山卡拉塔格地区航磁数据的解释,获得了地下空间不同磁性差异的磁性体的空间分布特征,为进一步分析研究区隐伏矿床提供了重要信息.     

欧拉反褶积与解析信号相结合的位场反演方法

由于解析信号具有不受(二维)或少受磁化方向影响,能够较好反映磁性体边界的特性,因此受到人们的重视.欧拉反褶积法可以确定场源的位置和深度以及形状因子,具有较强的适应性.因此前人提出将二者相结合的方法.针对前人提出的方法中存在受高频干扰严重的问题,本文提出低阶的欧拉反褶积与解析信号相结合的位场反演方法.本方法在反演中只需计...     

磁异常揭示的峨眉山大火成岩省的深部结构

峨眉山大火成岩省位于中国西南部,在晚二叠纪约260 Ma喷发出巨量的大陆溢流型玄武岩.对于大火成岩省的岩浆喷发,在地下必定有一个相应的大规模岩浆聚集和运移系统.地球物理方法是探测岩石圈内部的有效方式.峨眉山大火成岩省为镁铁质岩浆喷发,由于镁铁质-超镁铁质岩石一般具有强磁性,因此,在喷发结束之后,地下岩浆系统如果被镁铁质岩浆填充,冷却固化成为岩石圈的一部分,很有可能会引起磁异常.本文使用区域磁异常数据来对峨眉山大火成岩省的深部构造进行研究.该区域的磁异常由一系列离散的异常组成,通过3D磁化率反演可以得到磁性体的空间分布.由于磁异常中具有明显的剩磁,直接使用经典的反演方法会有较大误差,我们首先将磁异常转换为弱敏感于磁化方向的磁异常模量,再使用模量数据进行3D反演,得到地下空间内磁异常源的分布.经过分析认为这些离散分布的磁异常源反映了岩石圈内部的镁铁质-超镁铁质侵入体.侵入体的位置可能反映了底侵和内侵的镁铁质岩浆固化形成的侵入体,代表镁铁质岩浆房位置或者岩浆运移的主要通道.     

Gravity Anomalies of 2.5-D Multiple Prismatic Structures with Variable Density: A Marquardt Inversion

We present an inversion technique based on the Marquardt algorithm to estimate the depth of a 2.5-D sedimentary basin in addition to the regional gravity anomaly that is associated with the residual gravity anomaly, wherein the density contrast varies parabolically with depth. Forward modeling is carried out through a derived analytical gravity expression of a 2.5-D vertical prism. Inversion of a theoretical gravity anomaly with and without a regional gravity anomaly illustrates the procedure that it is found to be insensitive to the regional gravity effect. Furthermore, the algorithm is exemplified with the gravity anomalies of the derived density-depth model of the Godavari subbasin, India with a parabolic density profile resulting in a more consistent geological model rather than a constant density profile. The main advantage of this method is that it works well even when the profile of interpretation does not bisect the strike length of the sedimentary basin.     

渤海残留盆地分布综合地球物理研究

以活动论构造历史观为指导进行渤海海域残留盆地宏观分布特征研究. 以岩石物性为基础建立地质界面与物性界面的联系,通过重、磁、震等方法组合圈画不同物性界面的分布特征. 浅层沉积盆地结构通过反射地震资料控制,深层地质结构信息利用重磁异常场分离、物性界面反演等方法得到. 在地层分布格架基础上分析前新生代残留盆地宏观分布范围. 并结合渤海海域实例研究介绍了残留盆地宏观分布研究的技术方法和流程,在计算上古生界-中生界及元古界-下古生界残余地层厚度的基础上给出了渤海海域残留盆地的宏观分布范围,并指出在歧口凹陷东南部、石臼坨隆起南部及东北部、渤南凸起南部、辽东湾等地区是较好的前新生代油气潜力区.     

鄂尔多斯盆地上地壳速度分布与沉积建造和结晶基底起伏的构造研究

基于在鄂尔多斯盆地与其北部造山带地域,即沿延川—包头—满都拉地带进行的地震宽角反射和折射波场探测,取得了高分辨率的Pg波震相.通过走时差分层析成像方法进行Pg波波场走时反演,给出了沿剖面辖区的上地壳速度分布,求得了沉积建造和结晶基底折射界面的起伏变化,并给予了解释.研究结果表明,鄂尔多斯盆地上地壳为双层结构,上下层之间存在明显的折射界面,上层速度低,纵向变化梯度大;下层速度高,变化较为均匀.基于沿剖面辖区上地壳的速度分布特征提出自南向北应分为：榆林南凹陷、榆林—刀兔隆起、刀兔北至鄂尔多斯北缘断裂为箕形凹陷及其内部的次级构造、呼包凹陷以及伴随的断裂等沉积建造和结晶基底的起伏变化.阴山造山带上地壳速度明显比两侧地区高,速度呈纵向条带状展布,故呈现出结晶基底结构的分布特征和乌拉山、色尔腾山、和教岩体及白云鄂博群陆壳拼合及增生现象.内蒙造山带上地壳亦为双层结构,但基底折射界面不如鄂尔多斯块体明显.本文基于对该区沉积建造和结晶基底的起伏,讨论了沿剖面各有关凹陷的沉积特征与油气前景.     

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.     

基于磁力资料的南海及邻区磁性基底和居里面深度特征研究

磁性基底和居里面是研究地壳和岩石圈的地质构造和热演化过程的两个重要磁性界面.为了研究南海及邻区磁性基底和居里面所反映的深部构造及其热活动的地质效应,本文在对磁异常进行化极处理的基础上,采用最小曲率位场分离方法,获得了磁性基底和居里面引起的化极磁异常,利用双界面模型快速反演方法,反演了南海及邻区的磁性基底和居里面深度,研究了磁性基底、居里面深度及其分布特征,讨论了磁性基底、居里面与新生界深度之间相关性特征及其地质意义.研究表明,磁性基底深度5~20 km,洋盆南北两侧磁性基底走向分别以NE、NEE向为主,中南半岛周缘磁性基底呈NW、NNW走向.居里面深度15~32 km,宏观表现为"洋壳浅、周缘深"及周缘"北浅南深"的特征,洋盆地区居里面深度呈现"西南浅、东部深",洋壳与陆壳接触带在居里面深度上表现为梯级带特征.新生界深度与磁性基底深度相关性（Correlation between the depth of magnetic basement and Cenozoic,CDMBC）多以不规则形状分布,在盆地的沉积中心呈现正相关;新生界深度与居里面深度相关性（Correlation between the depth of Curie surface and Cenozoic,CDCSC）多呈NE、NEE向带状正相关分布,走向与盆地走向一致;莺歌海盆地、琼东南盆地、万安盆地南部和曾母盆地CDMBC呈正相关、CDCSC呈负相关,莺歌海相关性特征推测为：居里面随岩石圈变形隆起而抬升,磁性基底张裂下沉,发生大规模沉降引起;琼东南盆地相关性特征推测为：居里面随岩石圈变形下坳而下降,沉积中心与磁性基底下沉方向一致;万安盆地和曾母盆地相关性特征推测为：深部流体沿南海西缘断裂直接进入地壳,引起该处居里面深度变浅.