The distribution of deep source rocks in the GS Sag: joint MT–gravity modeling and constrained inversion

The coal-bearing strata of the deep Upper Paleozoic in the GS Sag have high hydrocarbon potential. Because of the absence of seismic data, we use electromagnetic (MT) and gravity data jointly to delineate the distribution of deep targets based on well logging and geological data. First, a preliminary geological model is established by using three-dimensional (3D) MT inversion results. Second, using the formation density and gravity anomalies, the preliminary geological model is modified by interactive inversion of the gravity data. Then, we conduct MT-constrained inversion based on the modified model to obtain an optimal geological model until the deviations at all stations are minimized. Finally, the geological model and a seismic profile in the middle of the sag is analysed. We determine that the deep reflections of the seismic profile correspond to the Upper Paleozoic that reaches thickness up to 800 m. The processing of field data suggests that the joint MT–gravity modeling and constrained inversion can reduce the multiple solutions for single geophysical data and thus improve the recognition of deep formations. The MT-constrained inversion is consistent with the geological features in the seismic section. This suggests that the joint MT and gravity modeling and constrained inversion can be used to delineate deep targets in similar basins.     

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.     

Crustal structure of eastern Cuba,derived by constrained 3D gravity inversion

A three-dimensional crustal model for Eastern Cuba, obtained through a process of gravity data inversion is presented. The study area cover a rectangular area of 64 600 km². The initial model for the inversion was constrained by surface geology, seismic and drilling data. The inversion algorithm uses gravity data to estimate 3-D topographies from the main geological units. The model provides quantitative information on the depths and thicknesses of the geological formations. The resulting model provides new information about the regional composition of the crust. Alien sequences are observed with different compositions and origin over the basement of Bahamas carbonate platform. Most of the maximum gravity anomalies are associated with presence of dense shallow ophiolite sheets. The most remarkable detail is the gravity “southwest” maximum, related to the presence of denser oceanic crust generated in the Cayman spreading center.     

九江-瑞昌矿集区的3D结构及对区域找矿的启示

本文分析了九瑞地区重、磁场的分布特征.运用最新的方法技术对已有的重磁资料重新进行处理:采用重磁多尺度边缘检测方法,对九瑞矿集区区域重力和航磁数据进行了边缘检测,并根据检测结果重新厘定了断裂系统的展布位置.在整理、分析九瑞地区地质、地层物性资料的基础上,对实测的1∶5万重磁数据进行较细致的准三维反演.同时将重磁三维物性反演应用到岩浆岩空间结构研究中,获得了矿集区地层结构及岩浆岩三维空间形态特征.根据反演所得磁化率强弱,分析了岩体的基性程度,为寻找与火山岩、侵入岩体有关的金属矿产提供了指示信息.最终建立的模型给出了地下地层的分布特征、控矿构造的展布规律、与成矿相关岩体的三维形态以及已知矿点的空间分布特征,为在九瑞矿集区的深部寻找隐伏矿体提供了新的信息.     

3-D gravity and magnetic interpretation for the Haifa Bay area (Israel)

Recently observed features in the subsurface geology of the Haifa Bay area (northern Israel) have been evaluated using 3-D forward gravity and magnetic modeling and inversion schemes. The interpretation is based on updated petrophysical data of the Jurassic, Cretaceous and Tertiary sedimentary layers and volcanics. It has been shown that the Bouguer gravity anomalies correspond mainly to thickness variations in the Senonian to Tertiary sediments. The gravity effect of these sediments was calculated using their actual densities and structural setting as interpreted from seismic reflection data. This effect was removed from the Bouguer gravity in order to study the pre-Senonian geological structures. The pattern of residual gravity anomalies (named “stripped gravity”) is essentially different from the pattern of the Bouguer gravity. The prominent Carmel gravity high, clearly seen on the Bouguer gravity map, completely vanishes on the “stripped” gravity map. That suggests that this relatively positive anomaly is caused by the considerable thickness of the low-density young sediments in the surrounding areas and does not correspond to high-density magmatic rocks or crystalline basement uplift as previously suggested. The average densities of the Jurassic and Cretaceous volcanics are generally lower then those of the background sedimentary rocks. Volcanics are the main cause for magnetic anomalies onshore and offshore northern Israel. The magmatic root of the Asher volcanics is, most probably, located close to the Yagur fault. A large, deep-seated gabbroic intrusion is assumed to be located under the Mediterranean abyssal plain in the NW part of the study area. The Atlit marine gravity low appears to be caused by a thick Mesozoic and Tertiary sedimentary accumulation. The results presented should be of considerable assistance in delineating some aspects of hydrocarbon exploration in the area.     

Magnetic polarity during the Middle Jurassic as recorded in the Summerville and Curtis Formations

Three-quarters of the Summerville Formation and a small portion of the Curtis Formation of central Utah have been sampled at close stratigraphic spacing. The two formations constitute an intertonguing marine-marginal marine, tidal flat sequence representing middle Callovian time. The magnetic polarity sequence observed is dominantly reversed with a few relatively short normal intervals. The age of these formations corresponds to some portion of time within the Jurassic oceanic quiet zone. The fact that normal polarity was not prominent in these formations casts some doubt upon the hypothesis that a long period of normal polarity was the cause of the quiet zone. Furthermore, a literature compilation indicates that normal polarity did not dominate the Jurassic as had been thought. The pole position determined from the Middle Jurassic Summerville Formation slightly revises earlier results and suggests a tiny amount of apparent polar motion relative to North American between Late Triassic and Middle Jurassic time.     

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.     

Feasibility study on application of satellite formations for eliminating the influence from aliasing error of ocean tide model

Currently,aliasing error of temporal signal model becomes the main factor constraining the accuracy of temporal gravity field.In provision of three types of satellite formations,i.e.,GRACE-type,Pendulum-type and n-s-Cartwheel-type,which are suitable for gravity mission and composed of observation in different directions,here we design two cases and conduct a simulation experiment on the feasibility to apply satellite formations for eliminating the influence from the aliasing error of ocean tide models.The result of our experiment shows that,when the aliasing error is disregarded,n-s-Cartwheel formation can provide the best conditions for gravity field determination,which,compared with GRACE-type,can improve the accuracy by 43%.When aliasing error of the ocean tide model acts as the main source of error,the satellite formation applied in dynamic method for gravity field inversion cannot eliminate aliasing or improve the accuracy of gravity field.And due to its higher sensitivity to the high-degree variation of gravity field,the Cartwheel-type formation,which includes the radial observation,can result in the gravity field containing more high-frequency signals for the ocean tide model error,and lead to a dramatically larger error.     

Accumulation characteristics and the main controlling factors of Lunnan multilayer oil province,Tarim Basin,China

Lunnan region is a large-scale paleohigh with many coexisting oil and gas bearing series.At present, about 2 billions tons of proved,probable and possible oil and gas reverses have been proved there.Eight oil and gas bearing series have been found in the Ordovician,Carboniferous,Triassic and Jurassic of Lunnan region,they all bear the characteristics of large-scale multilayer oil-gas province.Ordovician is the main reservoir series where over 0.8 billion tons of oil geologic reserves were discovered,and a super large-scale marine carbonate oil and gas field has formed.Reservoir space of the carbonate reservoirs is mainly composed of dissolved hole,dissolved pore and fracture in Lunnan paleo-burial hill.Generally, dissolved holes are widely distributed among them.Reservoir developments are mainly controlled by karstification and tectonic disruption.Due to the similar geochemical characters,the Ordovician,Carboniferous,Triassic and Jurassic oil and gas reservoirs present the same oil source rock of Mid-Upper Ordovician,the latter except Ordovician are mostly of secondary oil and gas reservoirs migrated vertically by faults during the process of multiple phase tectonic movement,adjustment and reconstruction. Lunnan composite oil and gas accumulation region is situated in the vicinity of large-scale hydrocarbon generation depressions in three directions,ample oil and gas from hydrocarbon generation depressions supplied the adjacent oil and gas reservoirs once.Hereby,the succeed paleohigh is the long-term hydrocarbon accumulation region,which is favor for the formations of high quality reservoirs,fault systems and huge-scale composite oil and gas accumulation.     

Accumulation characteristics and the main controlling factors of Lunnan multilayer oil province,Tarim Basin,China

Lunnan region is a large-scale paleohigh with many coexisting oil and gas bearing series. At present, about 2 billions tons of proved, probable and possible oil and gas reverses have been proved there. Eight oil and gas bearing series have been found in the Ordovician, Carboniferous, Triassic and Jurassic of Lunnan region, they all bear the characteristics of large-scale multilayer oil-gas province. Ordovician is the main reservoir series where over 0.8 billion tons of oil geologic reserves were discovered, and a super large-scale marine carbonate oil and gas field has formed. Reservoir space of the carbonate reservoirs is mainly composed of dissolved hole, dissolved pore and fracture in Lunnan paleo-burial hill. Generally, dissolved holes are widely distributed among them. Reservoir developments are mainly controlled by karstification and tectonic disruption. Due to the similar geochemical characters, the Ordovician, Carboniferous, Triassic and Jurassic oil and gas reservoirs present the same oil source rock of Mid-Upper Ordovician, the latter except Ordovician are mostly of secondary oil and gas reservoirs migrated vertically by faults during the process of multiple phase tectonic movement, adjustment and reconstruction. Lunnan composite oil and gas accumulation region is situated in the vicinity of large-scale hydrocarbon generation depressions in three directions, ample oil and gas from hydrocarbon generation depressions supplied the adjacent oil and gas reservoirs once. Hereby, the succeed paleohigh is the long-term hydrocarbon accumulation region, which is favor for the formations of high quality reservors, fault systems and huge-scale composite oil and gas accumulation.     

重力和重力梯度数据联合聚焦反演方法

重力数据包含较多的低频信息,重力梯度数据包含较多的高频信息,将重力数据和重力梯度数据进行联合反演得到的结果更加可信.本文基于聚焦反演方法,实现了这一过程.因为联合反演中分量种类增加,所以计算灵敏度矩阵所需要的时间增加,为此,本文提出了一种快速计算灵敏度矩阵的方法.因为联合反演对内存的要求增大,本文选择有限内存BFGS拟牛顿法求解反演问题.本文通过再加权的方法实现深度加权.文中利用单一分量的反演结果来预测异常体的埋深信息,随后将埋深信息结合到深度加权函数中,将其用于多分量组合反演计算.给出了模型试验,发现预测得到的异常体的埋深信息与其实际埋深存在偏差,但是将这一信息应用到反演计算,能够得到与真实模型一致的结果.之后,本文通过模型试验来探究重力和重力梯度联合反演的优势,发现将重力和重力梯度数据联合,能够识别出额外的噪声,反演得到的模型更加合理.但是,对于不同分量组合得到的反演结果是相近的,反演模型的提高很小.最后,将联合反演方法应用到美国路易斯安那州Vinton岩丘的实际数据中,结果显示,将重力和重力梯度数据联合反演,反演模型得到了提高,反演得到的结果与地质资料吻合.     

Accumulation characteristics and the main controlling factors of Lunnan multilayer oil province,Tarim Basin,China

Lunnan region is a large-scale paleohigh with many coexisting oil and gas bearing series. At present, about 2 billions tons of proved, probable and possible oil and gas reverses have been proved there. Eight oil and gas bearing series have been found in the Ordovician, Carboniferous, Triassic and Jurassic of Lunnan region, they all bear the characteristics of large-scale multilayer oil-gas province. Ordovician is the main reservoir series where over 0.8 billion tons of oil geologic reserves were discovered, and a super large-scale marine carbonate oil and gas field has formed. Reservoir space of the carbonate reservoirs is mainly composed of dissolved hole, dissolved pore and fracture in Lunnan paleo-burial hill. Generally, dissolved holes are widely distributed among them. Reservoir developments are mainly controlled by karstification and tectonic disruption. Due to the similar geochemical characters, the Ordovician, Carboniferous, Triassic and Jurassic oil and gas reservoirs present the same oil source rock of Mid-Upper Ordovician, the latter except Ordovician are mostly of secondary oil and gas reservoirs migrated vertically by faults during the process of multiple phase tectonic movement, adjustment and reconstruction. Lunnan composite oil and gas accumulation region is situated in the vicinity of large-scale hydrocarbon generation depressions in three directions, ample oil and gas from hydrocarbon generation depressions supplied the adjacent oil and gas reservoirs once. Hereby, the succeed paleohigh is the long-term hydrocarbon accumulation region, which is favor for the formations of high quality reservors, fault systems and huge-scale composite oil and gas accumulation.

     

三维多层介质重力-地震同步联合反演

联合反演是地球物理勘探的重要解释手段,能够提高模型参数的反演精度.本文在归纳和分析重力与地震资料联合反演的研究和应用现状的基础上,利用三维多层介质模型的地震走时和重力正演公式,推导了地震走时和重力异常对界面深度的雅可比矩阵,实现了三维重力-地震同步联合反演界面成像.最后进行了数值理论模型模拟和实例计算,结果表明地震走时和重力同步联合反演很好的重建了三维多层介质界面.     

基于模型降阶的贝叶斯方法在三维重力反演中的实践

三维重力反演是地质工作者了解地球深部构造,认知地下结构的重要手段.按照反演单元划分,三维重力反演有离散多面体(Discrete)反演和网格节点(Voxels)反演两种方式.离散多面体反演由于易于吸收先验地质信息得到的理论场能够很好地拟合观测场,因此,在实际重力反演中更受欢迎.目前离散多面体重力反演中初始模型的建立方法繁杂不一,实际应用受到很大的限制.本文本着充分挖掘利用先验信息和重力观测数据得到丰富可靠的反演结果这一原则,以离散多面体反演技术为基础,改进建模过程.在初始模型的建立中,吸收贝叶斯算法优势,采用隐马尔科夫链改善朴素贝叶斯方法的分类效果,通过最大似然函数算法求解,再采取模型降阶技术,固定所建模型中几何体的形态或密度,达到在几何体形态(x,y,z)、密度(σ)和重力值(g)五个参数中降低维数目的,从而减小高维不确定性和正演的计算量,由此反演计算的地质体密度和分布范围相对更准确,更利于重现重力模型结构.通过单位球体和任意形态几何体模拟实验,以及安徽省泥河矿区三维重力反演实践,得到非常接近实际的密度或重力值,大幅提高了三维重力反演的精度和效率,说明该方法是有效、实用的.     

基于深度学习的重力异常与重力梯度异常联合反演

高效高精度的反演算法在重力大数据时代背景下显得尤为重要,受深度学习卓越的非线性映射能力的启发,本文提出了一种基于深度学习的重力异常及重力梯度异常的联合反演方法.文中首先提出了一种基于网格点几何格架的重力异常及重力梯度异常的空间域快速正演算法,这为本文深度学习反演算法的实现奠定了基础;其次对大量的不同密度模型进行正演计算获得样本数据集;然后设计了一种端到端的深度学习网络结构(GraInvNet),再利用样本数据对该网络结构进行训练;最后进行反演预测.组合模型试验表明,多维度数据联合反演相比单一分量反演其结果更“聚焦”,且与模型边界高度吻合,并且对于复杂模型的姿态与物性预测具有极为显著的优势,以及对于含噪声数据的反演,其质量也不会降低;Vinton岩丘实测重力数据也验证了文中方法的有效性;从而证明了深度学习在重力数据的高效高精度反演方面具有的巨大潜力.     

Gravity compression forward modeling and multiscale inversion based on wavelet transform

The main problems in three-dimensional gravity inversion are the non-uniqueness of the solutions and the high computational cost of large data sets. To minimize the high computational cost, we propose a new sorting method to reduce fluctuations and the high frequency of the sensitivity matrix prior to applying the wavelet transform. Consequently, the sparsity and compression ratio of the sensitivity matrix are improved as well as the accuracy of the forward modeling. Furthermore, memory storage requirements are reduced and the forward modeling is accelerated compared with uncompressed forward modeling. The forward modeling results suggest that the compression ratio of the sensitivity matrix can be more than 300. Furthermore, multiscale inversion based on the wavelet transform is applied to gravity inversion. By decomposing the gravity inversion into subproblems of different scales, the non-uniqueness and stability of the gravity inversion are improved as multiscale data are considered. Finally, we applied conventional focusing inversion and multiscale inversion on simulated and measured data to demonstrate the effectiveness of the proposed gravity inversion method.     

空-地-井重力异常正则化协同密度反演方法

围绕综合利用空-地-井多源重力异常联合反演提升精度的目标,提出正则化协同密度反演方法,从而有效利用多维数据的横向和纵向变化特征提高反演精度,且无需先验信息的约束.此外,还提出利用奇异值谱和深度分辨率工具评价航空数据观测高度和层数对反演分辨率的影响.通过理论模型试验本文方法在空-地、地-井等不同情况下的应用效果,结果表明空-地-井重力异常正则化协同密度反演方法能获得更高分辨率、高精度的反演结果,且证明不同的钻孔位置对反演结果有较大的影响,从而可指导实际的空-地-井联合勘探.最后,利用文顿盐丘实际勘探数据进行空地和空地井重力数据协同反演,反演结果垂向分辨率明显优于地面观测数据反演结果,盐盖的顶面及中心埋深与前人研究和地质资料解释相吻合,验证了方法正确性及实用性,为推进我国空-地-井立体重磁勘探提供了重要的技术手段.     

背景噪声面波与布格重力异常联合反演:山西断陷带三维地壳结构

讨论了利用面波与布格重力异常联合反演三维地壳速度结构的新方法,并利用该方法联合反演获得山西断陷带地壳S波速度结构.通过建立速度与密度之间的经验关系,利用非线性迭代反演方法获得最终速度模型.结果显示,联合反演获得的速度模型可以同时提高对面波及重力数据的观测拟合程度,而面波单独反演得到的速度模型则无法很好的拟合重力观测数据.相比较,联合反演速度模型中的大同火山区中下地壳的低速异常幅值小于面波单独反演模型中低速异常体的幅值.联合反演速度模型结果揭示,吕梁山地区在中下地壳存在低速异常,并且和北部的大同火山区低速异常相连接,说明可能导致新生代以来大同火山区岩浆活动的上地幔构造活动(上地幔局部上涌,地幔柱)可能对山西断陷带的形成和构造活动起到了一定的控制作用,并且导致了吕梁山地区中下地壳的低速异常.     

基于3DDelaunay剖分算法的重力建模与分析

地质体的重力建模是正确解释和应用重力资料的关键问题之一.针对非规则形状变密度的三度体,本文提出了基于3D Delaunay剖分算法的重力建模方法.采用3D Delaunay剖分算法将三维目标地质体分解为若干变密度四面体体元,推导了基于四面体体元的重力正演公式,建立了剩余密度值与重力异常值的线性方程组;以变密度的长方体和倾斜台阶组合体为例,比较分析了常规块体算法和3D Delaunay剖分算法应用于重力正演的有效性,并采用共轭梯度法加密度约束条件对非规则形状变密度的倾斜台阶组合体进行了密度反演.计算结果验证了本文方法的正确性和有效性.基于3D Delaunay剖分算法的重力建模可应用于存在褶皱、断层、裂缝等复杂地质体的重力正反演计算.     

基于地质体空间位置优化约束的航空重力梯度数据三维物性反演

重力数据的物性反演面临着严重的多解性问题,降低多解性的有效手段是加入约束条件.而边界识别、深度估计及成像方法可获取地质体的水平位置、深度范围等几何参数信息,本文将基于数据本身挖掘的地质体几何参数信息约束到物性反演中,以降低反演的多解性.通过引入基于深度信息的深度加权函数及基于水平位置的水平梯度加权函数建立优化约束条件,有效地提高了反演结果的横向及纵向分辨率.重力梯度数据包含更多的地质体空间特征信息,将优化约束反演方法应用到全张量数据的反演中,模型试验表明本文方法反演结果与理论模型更加吻合.最后对美国路易斯安那州文顿盐丘实测航空重力梯度数据的应用表明,本文方法在其他地球物理、地质资料不足的情况下获得更可靠的反演结果.