Gravity interpretation of nonoutcropping sedimentary basins in which the density contrast decreases parabolically with depth

The decrease in density contrast of sedimentary rocks with depth in many sedimentary basins can be approximated by a parabolic density function. Analytical gravity expression of an outcropping two-dimensional vertical step along which the density contrast decreases parabolically with depth is derived in the space domain. A modification ofBott's (1960) method of gravity interpretation is proposed by considering two outcropping vertical steps on either side of the first and last observation points in addition toN outcropping vertical prisms in order to interpret the gravity anomalies of nonoutcropping basins. The thicknesses of the two outcropping vertical steps are made equal to the thicknesses of the two outcropping vertical prisms placed below the first and last observation points. The initial depth estimates of the sedimentary basin are calculated by the infinite slab formula ofVisweswara Rao et al. (1993). The gravity effects of theN outcropping prisms and the two outcropping vertical steps are calculated at each anomaly point and the depth to the floor of the basin are adjusted based on the differences between the observed and calculated anomalies. A gravity anomaly profile of Los Angeles basin, California is interpreted.     

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.     

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.     

Gravity interpretation of sedimentary basins with hyperbolic density contrast1

The variation in the density of sediments with depth in a sedimentary basin can be represented by a hyperbolic function. Gravity anomaly expressions for a 2D vertical prism and an asymmetric trapezium with a hyperbolic density distribution are derived in a closed form. These are used in inverting the gravity anomaly of a sedimentary basin with variable density. Firstly, the basin is viewed as a series of prisms juxtaposed with each other. The initial thickness of each prism is obtained from the gravity anomaly at its centre, based on the gravity anomaly of an infinite slab with a hyperbolic density contrast. These thicknesses are improved, based on the differences between the observed and the calculated anomalies. For an improved rate of convergence of the solution, these thicknesses may alternatively be refined using the well-known ridge regression technique. Secondly, the basin is approximated by an asymmetric trapezium and its anomalies are inverted for the parameters of the trapezium using the ridge regression. Since this approximation serves to oversimplify the floor of the basin, it must be used only when the sediment-basement interface has minor undulations. The results of a hypothetical case and two field cases (the San Jacinto Graben, California and the Godavari Graben, southern India) are presented. In both field cases, the interpreted depths are comparable with the real ones, proving the validity of the assumption of a hyperbolic density distribution of the sediments in the two basins considered.     

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.     

3D AND 2½ D MODELLING OF GRAVITY ANOMALIES WITH VARIABLE DENSITY CONTRAST1

The decrease of density contrast in sedimentary basins may be approximated by a quadratic function. A sedimentary basin may be viewed as a number of prisms placed in juxtaposition. Equations in closed form for the gravity anomalies of 3D and 2½ D prismatic models are derived. Approximate equations for these models are also derived for rapid calculations. Efficient methods are developed for anomaly calculation by an appropriate use of the exact and approximate equations, and hence, for 3D and 2½ D modelling. The depths to the basement are adjusted iteratively by comparing the calculated anomalies with the observed anomalies. These methods are applied for analysis of the residual anomaly map of the Los Angeles Basin, California.     

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

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

THREE-DIMENSIONAL INTERPRETATION OF GRAVITY DATA FROM SEDIMENTARY BASINS USING AN EXPONENTIAL DENSITY-DEPTH FUNCTION*

In mapping the topography of the basement of deep sedimentary basins by gravity modelling, the accuracy can be improved by incorporating an exponential increase in density with depth. For calculating the gravity effect of a three-dimensional (3D) structure with such an exponential density-depth relation a frequency-domain forward algorithm based on series expansion is presented, the numerical evaluation of which can be performed efficiently by fast Fourier transform. The algorithm can be applied in a recursive procedure to give the inverse solution in terms of basement relief. The inversion procedure is satisfactorily tested on a 2D synthetic example and a 3D field example of gravity data from the western margin of the Pannonian Basin in eastern Austria, where up to 2.2 km of Tertiary sediments overlie an igneous or metamorphic basement. The results are confirmed by basement intersections in several wells.     

A new technique of interpreting gravity and vertical magnetic anomalies due to infinite horizontal cylindrical ore bodies

Summary A novel method of interpreting gravity and magnetic anomalies is presented here. Two diagrams of master curves, one for gravity and the other for vertical magnetic interpretation are presented. They are useful in calculating the depth of burial and the radius of infinite horizontal cylindrical bodies irrespective of density contrast or strength and direction of magnetisation. This method also enables us to infer the direction of magnetisation, resultant intensity of magnetisation and from them the susceptibility contrast; and density contrast. ThoughHenderson [8]²) reported the applicability of continuations and derivatives in the interpretation of magnetic anomalies for the first time, the authors have treated them in a more exhaustive manner in the present paper to obtain valuable relationships.     

History of faulting and magmatism in the Galilee (Israel) and across the Levant continental margin inferred from potential field data

The Galilee study area, northern Israel, is at present an uplifted, steep continental margin that formed mainly during the Jurassic and has a large positive isostatic anomaly. Since the Jurassic, it was modified by several tectonomagmatic events, which this study attempts to define and classify by updating, reprocessing and reinterpreting gravity, aeromagnetic and geological data. The prominent Rehovot-Carmel N–S positive reduced-to-pole (RTP) magnetic anomaly caused by the Gevim Volcanics, as well as the coexisting Helez-Gaash high Bouguer gravity and the Pleshet low Bouguer gravity, represent the deep (>5 km) Permo-Triassic dominant horst and graben structure of Israel. The Jonah Ridge and Beirut high SW–NE RTP magnetic anomalies in the Levant basin delineate the Levant continental edge that is marked by a deeply buried horst covered by a Late Cretaceous volcanic complex. The Asher and Devora Jurassic volcanics appear to be responcible for the Atlit and Galilee negative magnetic anomalies and for significant negative gravity anomalies which became clear after removing gravity effect of the upper (post-Turonian) light density sediments from the observed gravity. The volcanics extend along a SW–NE belt parallel to the strike of the Moho. It is suggested here that the Carmel-Gilboa fault propagated during the Late Cretaceous from the Levant basin across the Galilee area southeastward to form the Azraq-Sirhan graben in Jordan. As such, it forms a right-step, en echelon, dextral strike-slip fault with associated tectonic basins of various shapes. During the Oligocene and before formation of the Dead Sea transform (DST), the reactivation of the Azraq-Sirhan graben was accompanied by tectonic driven rift propagation in the opposite direction, from Azraq-Sirhan to northwest. It dispersed into many faults and terminated ∼10 km west of the present DST. During the Miocene it propagated in the same direction and includes internal volcanic activity. The numerous Miocene-Pliocene volcanic centers on the margins of the DST indicate that the preferred pathway for magmas at that time was not within the deep basins of the DST.     

东海陆架盆地及其周边海域地质、地球物理场特征

分东海陆架盆地及其周边海域的重、磁场特征、对研究区域的地球物理场资料进行了小波变换与多尺度分布与计算，计算了研究了海域重，磁资料的1－4阶小波变换逼近、细节、根据分析与计算可知，在东海陆架盆地及其周边域，重、磁异常的主要特征没有内陆盆地那样明显，但盆地内的绝大部分地区均位于异常相对缓变区，在盆地内部有时有局部异常的圈闭及剧变区，但总的特征为缓变的，而周边区域多为相对剧变区，但也存在局部区域的缓变, 、磁场特征表明，盆地基底相对盆地周边区域主要由密度较大，磁性中的岩石组成，因此盆地中部以高正磁异常为主，布格重力异常 对周边地区高，东海陆架盆地及其周边海域的重力场的主要走向为NNE向，局部重力异常较为平缓，磁异常走向没有重力异常明显，但主要也为NNE向，重、磁场的分布特征表明，东海陆架盆地及其邻域的地质构造为大陆架构造的延伸。     

Marquardt optimization of gravity anomalies of anticlinal and synclinal structures with prescribed depth-dependent density

An inversion technique using the Marquardt optimization is developed to interpret the gravity anomalies due to anticlinal and synclinal structures with density contrast varying continuously with depth. The algorithm simultaneously estimates the parameters of the respective models, in addition to the regional gravity background that is invariably associated with the residual gravity anomaly. Forward modelling is realized through analytically derived gravity expressions for the respective models in the space domain. The efficacy of the inversion is demonstrated with the gravity anomaly due to a theoretical model, in each case with and without the regional background. In addition, the applicability is illustrated using the gravity anomalies of the Pays De Bray anticline, situated north‐west of Paris, France. The interpreted depth of the Pays De Bray anticline using the present inversion compares well with the drilling depth.     

A geophysical approach to the granite batholith under the eastern Southern Uplands,Scotland

From our interpretation of the Bouguer gravity and aeromagnetic anomalies in south-east Scotland, we conclude that a massive granite batholith underlies the greater part of the eastern Southern Uplands. The granite model which we computed earlier from gravity anomalies in the Tweeddale area fits the observed magnetic anomalies closely, if a normal magnetization of 0.095 A m^–1 is assigned, similar to values found for exposed local granites. Further gravity modelling shows that, apart from the Tweeddale boss, the granite shallows to less than 1 km near Lammer Law in East Lothian and extends north of the Lammermuir Fault. A model for the East Lothian volcanics was computed from their aeromagnetic anomalies, then their gravitational effect was combined with that estimated for the Devonian and Carboniferous sediments and the result stripped off the observed gravity field. The residual gravity anomalies were used to generate a two-dimensional model for the granite north of the Lammermuir Fault. The expected tectonic consequences of a massive granite batholith in the eastern Southern Uplands are compared with the known development of faults and sedimentary basins around its margins.     

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

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

Consequences of foreland basin development on thinned continental lithosphere: Application to the Aquitaine basin (SW France)

We present an analysis of the consequences of foreland basin development on thinned continental lithosphere, inherited from pre-orogenic phases of extension. Bathymetry at the transition from pre-orogenic extensional basin to foreland basin and compaction of pre-orogenic sediments contribute to the accommodation space for foreland basin sediments and thrust loads. In addition, the extension-induced transient thermal state of the lithosphere, results in ongoing thermal subsidence, and a flexural rigidity which changes through time. Quantitative modelling of the phase of extension and the foreland basin stage of the Aquitaine basin (southern France) shows that the inherited transient thermal state of the lithosphere contributes significantly to (1) the total foreland basin depth and width, (2) the post-compressional subsidence history, and (3) the cratonward onlap pattern. Accounting for the thermo-mechanical effects of pre-orogenic extension significantly reduces the estimates of both the flexural rigidity (30–43% for the Aquitaine basin) and the required topographic or thrust load (40% for the Aquitaine basin) at foreland basins. Emplacement of thrust loads below sea level, as expected in a pre-orogenic extensional basin setting, further reduces the required topographic load. This sheds light on the wide range of flexural rigidity values reported for continental lithosphere from foreland basin modelling studies, and explains, in many instances, the inferred ‘hidden load’ or subsurface load in flexural modelling studies at foreland basins. The present study has shown that pre-orogenic extension phases significantly affect the record of vertical motion and the stratigraphy of the Aquitaine basin and is probably important for foreland basin evolution in general.     

南海北部陆架西区盆地地质、地球物理场特征及其深部结构

分析了南海北部陆架西区盆地的地质、地球物理场特征，计算了研究海域重、磁资料的1阶小波细节、4阶小波逼近变换。根据分析与计算可知，研究区的布格重力异常以北西低的负值，东南高的正值为特征。在东部及东南部异常等值线走向为北东；西部异常等值线以北西走向为特征；西北地区异常以北东东、北东走向的局部等值线圈闭为特征。磁场的展布十分复杂，按磁异常的变化程度可分为三个变化区，即磁异常平静区、剧变区及缓变区。磁异常的平静区位于研究区的西部，即莺歌海盆地所在位置，这一带磁异常等值线极为稀疏，异常值为负背景异常。剧变区位于海南岛，该地区的磁异常变化极为剧烈，异常特征以局部小圈闭为特征，等值线分布密集。磁异常的平缓区位于平静区及剧变区之外的其它地区。琼东南盆地、北部湾盆地的磁异常具有此特征。根据重、磁场资料以及南海北部盆地钻井取样的测试结果、同时参考穿越南海地学断面的结果，对研究区的地壳结构进行了反演计算，计算表明南海陆架盆地区域地壳结构较为复杂，研究区的地壳厚度在22－33km之间，总的趋势由陆向洋地壳厚度逐渐减薄，反映出该区域地壳具有陆壳、拉伸陆壳、过渡壳的性质，同时存在有上地幔隆起区及凹陷区。磁性底界面厚度在17－24km之间变化，其中在莺歌海盆地较深，在海南岛地区磁性界面较浅。     

盆地基底岩性的综合地球物理预测方法——以松辽盆地滨北地区基底岩性预测为例

预测盆地基岩岩性不仅对于研究盆地的深部地质结构及盆地的形成演化具有重要的意义,而且也对基岩风化壳油气藏的勘探具有一定的指导作用.本文通过对盆地重、磁异常成因的综合分析,提出了一系列盆地基底岩性综合预测研究的综合地球物理资料处理解释方法技术.指出在地震构造界面的约束下采用重力剥皮技术可以较为可靠地获取基底岩性重力异常并分...     

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.     

大庆外围盆地地球物理场与盆地基底特征

根据收集到的大庆外围盆地的布格重力异常和航磁异常数据,对外围盆地中的8个盆地进行了地球物理场特征的分析.对各盆地的重力异常进行了由下地壳与上地幔的密度差引起的重力效应的剥离,对磁力异常进行了由居里面起伏引起的磁力效应的剥离.在此基础上应用调和级数法和遗传算法分别反演计算了这8个盆地的重力基底和磁性顶界面,并分析了其基底特征. 各盆地重力基底在0.2~9.0 km之间变化,磁性顶界面在1.8~9.8 km之间变化,其特征反映了各盆地的基本现状.     

基于重磁震资料的南海新生代盆地分布综合研究

作为西太平洋最大的边缘海,南海分布有30多个新生代沉积盆地,其蕴含着丰富的油气资源.但由于资料的限制,南海存在不同区域盆地研究程度不同,不同区域盆地面积差别较大,部分盆地只是坳陷而没有达到盆地的级别以及盆地外围可能存在凹陷等问题.南海新生代盆地分布问题制约了其油气分布规律、储量等基础地质问题的研究.本文以地震剖面数据为约束,以重力资料为主、辅以磁力资料,研究了南海新生代盆地分布及构造区划.通过提取新生代盆地及其构造单元引起的重力异常,结合地震剖面等资料反演了新生界底界面深度及新生界厚度.在充分调研已有盆地和构造单元划分方案的基础上,根据南海的地质及地球物理特征,确定了盆地及构造单元划分标准.以新生界厚度为基础并结合重、磁、震、地质等资料,进行地质-重磁震联合解释,将南海原有的36个盆地重新划分为24个盆地,盆地总面积扩大了约15万km².研究表明,南海新生代盆地沉积层厚度在1.5~16 km之间,有6个北东东/北东向沉积坳陷带、2个近南北向沉积坳陷带以及1个三角沉积坳陷区;盆地展布方向主要为北东和北东东向,其次为北西和近南北向,呈现"南三北三"的分布特征.