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

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

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

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

Estimation of the mass density contrasts and the 3D geometrical shape of the source bodies in the Yilgarn area, Eastern Goldfields, Western Australia

We invert 2D surface gravity data constrained both by geological and seismic information. We use a number of pre-processing tools in order to reduce the general multi-body inversion into several single-body inversions, whereby we can reduce the overall complexity of the inversion task. This is done with as few assumptions as possible. Furthermore, for a single-body inversion we uncouple the determination of the shape of the causative sources from the determination of their mass density contrast to the surroundings. The inversion for the geometrical shape of the source body is done in steps. Firstly, a rough 3D shape of the source is modelled—a model consisting of the vertical mass columns of equal height. The horizontal extension is implied by the surface gravity signal. Subsequently, the shape of each source body is modified to obtain a better fit to the surface gravity data. In each modification step, the overall change of the shape of the source body is followed by an update of the mass density contrast to the surroundings. The technique was applied to a set of gravity data from the Eastern Goldfield area in Western Australia. The area has been widely studied in the past. In 1999, two seismic profiles that cross-sect the area were measured. Furthermore, an extensive geological modelling for the area has been conducted. The practical goal of this work was to verify the geological interpretation using the potential field data (mainly the gravity data although magnetic data were also available) and only weakly constrained by the seismic information. The result was the reconstruction of the ‘rough’ 3D geometry of the source bodies and the estimation of a constant mass density contrast to the surroundings. A possible extension of this technique for detailed studies of the geological model is briefly discussed.     

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.     

A magnetotelluric investigation of the structural geology beneath Charlevoix Crater, Quebec

The results of ten tensor magnetotelluric soundings recorded within the Charlevoix crater are presented. The geological complexity of the region (the crater is located at the junction of two major geological provinces, i.e. the crystalline basement of the Canadian shield and the more conductive sediments of the Appalachian system) and the amount of cultural noise present in our data prevent a highly definitive interpretation of the shallow structure. A subhorizontal conductive fault related to Rondot's (1971) structural sketch of the original impact structure may exist within the upper 5 km, but it is not completely constrained by the data. Pseudo 3D modelling using stacked thin sheets indicates that superficial patches of high conductivity may explain the large anisotropy observed at the long periods and suggests that our deeper structures fit well with the 2D model previously proposed by Kurtz (1982) for the region.     

先验地质信息约束下的三维重磁反演建模研究 ——以安徽泥河铁矿为例

安徽泥河铁(硫)矿床是近年发现的大深度隐伏矿床,它的发现再次引发了庐枞地区深部找矿的热潮.建立矿区的三维精细地质模型,对寻找深、边部隐伏矿体,深入认识深部成矿、控矿规律意义重大.本文以泥河矿区为实例,开展先验地质信息约束的三维重磁建模研究,获得了矿区面积5.6 km²(2.8 km×2.0 km)、深度1.2 km内的三维地质模型.在三维可视化平台上对该模型进行了地质解释,全面分析了矿体、地层与次火山岩之间的空间分布及对应关系,发现铁矿主要赋存于闪长玢岩与砖桥组火山岩之间,且在玢岩穹窿地段矿体厚大,这对认识"玢岩型"铁矿的成矿模式具有重要的实际意义;三维模型的重磁正演响应基本拟合了实际重磁异常,说明在建模区域内已没有新的矿体;本文提出的三维地质建模流程可为其它地区开展类似工作提供借鉴,同时研究表明地质信息约束下的三维重磁建模研究在深、边部找矿和重磁异常的精细解剖等方面具有潜在的价值和广阔的应用前景.     

Detailed reservoir inversion addressing geological problems in reservoir development

垦71开发区块的主要产油层系为新近系的馆陶组河流相和古近系的东营组三角洲相地层,具有被断层复杂化的多油水系统。该区井网密集,储层厚度薄,横向关系复杂,储层与非储层速度差异小,后续开发方案调整所要求的精细储层描述难度较大。本文展示了以精细储层描述为目标的地震约束反演实例,利用多井地震、地质标定,结合三维构造解释成果,建立三维多断层精细约束地质模型。在此基础上,采用随机地震反演方法,进行了波阻抗、岩性等多种地层参数反演,较好地解决了该区薄互层储层精细描述问题,证实了精细储层反演技术在油田开发中的重要作用。     

面向油藏开发地质问题的精细储层反演研究

垦71开发区块的主要产油层系为新近系的馆陶组河流相和古近系的东营组三角洲相地层,具有被断层复杂化的多油水系统。该区井网密集,储层厚度薄,横向关系复杂,储层与非储层速度差异小,后续开发方案调整所要求的精细储层描述难度较大。本文展示了以精细储层描述为目标的地震约束反演实例,利用多井地震、地质标定,结合三维构造解释成果,建立三维多断层精细约束地质模型。在此基础上,采用随机地震反演方法,进行了波阻抗、岩性等多种地层参数反演,较好地解决了该区薄互层储层精细描述问题,证实了精细储层反演技术在油田开发中的重要作用。     

二度半长方体组合模型的重力模拟退火反演

针对重力异常反演的现状,采用矩形单元体网格剖分的二度半体组合模型,模拟任意的多层密度界面或形体,实现了重力正演.该建模方法适应性强,适合于实际工作中常见的密度横向变化的复杂模型.使用改进的全局寻优的快速模拟退火算法,对重力异常进行反演,结合这种灵活的密度建模方法,反演过程中只需要反演各矩形单元体的密度参数,即可同时得到地质体的界面或形态以及密度值分布.通过复杂界面和形体的模型试验,证明了方法的效果.在江苏下扬子地质构造复杂区的实际应用说明,在地震等先验信息约束下,该重力反演方法提高了反演精度并减少了多解性,可有效解决古生界泥盆系-志留系目的层分布等地质问题.     

基于先验信息约束的重磁三维交互反演建模技术——以铜陵矿集区为例

三维地质建模是实现深部矿产勘查突破的重要途径,其通过控矿地质体的三维建模,直观刻画控矿要素之间的空间、成因和演化关系,帮助理解成矿系统,开展深部找矿预测.当前三维建模主要采用地质资料构建,在缺少钻孔等已知资料的情况下,难于构建出可靠的三维地质模型.为了克服常规三维地质建模方法可信度低、精度差的缺点,本文将重磁交互反演技术引入到三维地质建模中,提出了基于先验信息约束,通过二度半剖面交互反演、三维物性反演联合修正的三维地质模型方法.采用该方法建立了铜陵矿集区的三维地质模型,并开展了深部找矿预测,取得以下主要认识:1)基于先验信息约束的重磁交互反演建模技术能大幅提高三维地质模型的可信度,是实现地下地质体"透明化"的重要途径;2)铜陵地区重要控矿地层(C-P-T2)主要分布于向斜区和火山岩覆盖区,深部岩浆岩条件优越,深部找矿工作应从传统的背斜隆起区转向向斜区和火山岩覆盖区以及部分凹陷区,这将极大拓展铜陵地区找矿空间;3)铜陵地区传统的五大矿田的侵入岩均具有复合岩体特征,特别是与辉石闪长岩共生的矿产以铁铜矿为主,而与金矿密切相关的侵入岩则多为中酸性岩体,这一推断如果成立,将对铜陵地区找矿突破产生重大影响;4)铜陵地区的推覆构造十分发育,木镇凹陷区的早古生代基底和南陵盆地北西缘的戴公山背斜就是一套区域性推覆构造的产物,三维地质模型还揭示铜陵地区侵入岩具有深、中、浅三重结构,这为铜陵隆起的推覆成因提供了证据,也指示了在铜陵之外再寻找类似铜陵矿集区的良好前景.     

A parallel,scalable and memory efficient inversion code for very large‐scale airborne electromagnetics surveys

Over the past decade the typical size of airborne electromagnetic data sets has been growing rapidly, along with an emerging need for highly accurate modelling. One‐dimensional approximate inversions or data transform techniques have previously been employed for very large‐scale studies of quasi‐layered settings but these techniques fail to provide the consistent accuracy needed by many modern applications such as aquifer and geological mapping, uranium exploration, oil sands and integrated modelling. In these cases the use of more time‐consuming 1D forward and inverse modelling provide the only acceptable solution that is also computationally feasible. When target structures are known to be quasi layered and spatially coherent it is beneficial to incorporate this assumption directly into the inversion. This implies inverting multiple soundings at a time in larger constrained problems, which allows for resolving geological layers that are undetectable using simple independent inversions. Ideally, entire surveys should be inverted at a time in huge constrained problems but poor scaling properties of the underlying algorithms typically make this challenging. Here, we document how we optimized an inversion code for very large‐scale constrained airborne electromagnetic problems. Most importantly, we describe how we solve linear systems using an iterative method that scales linearly with the size of the data set in terms of both solution time and memory consumption. We also describe how we parallelized the core region of the code, in order to obtain almost ideal strong parallel scaling on current 4‐socket shared memory computers. We further show how model parameter uncertainty estimates can be efficiently obtained in linear time and we demonstrate the capabilities of the full implementation by inverting a 3327 line km SkyTEM survey overnight. Performance and scaling properties are discussed based on the timings of the field example and we describe the criteria that must be fulfilled in order to adapt our methodology for similar type problems.     

郯庐断裂带江苏段新生界三维地质构造建模

本文基于SKUA-GOCAD三维建模软件平台的构造-地层实体建模流程,通过整合研究区大量地质、地球物理资料和数据利用地质时间域坐标系转换和离散光滑插值（DSI）技术,并结合非定量化的地质基础利用地质网格功能建立了郯庐断裂带江苏段新生界的三维地质构造模型,该模型包括地表、第四系底部、新近系底部和基岩顶部等四个层面的地层模型和断层模型。该模型定量化地给出了研究区第四系、新近系底界的三维结构和特征,并刻画了主要活动断裂的三维空间分布,通过该模型可以更直观地显示研究区新生代地层单元和断裂的变化规律,对深入分析活动断裂带的沉积与构造具有重要的意义。      

重磁数据三维物性反演方法进展

综述了重磁数据三维物性反演方法中的几个关键问题.主要包括正演快速算法、反演框架、约束因子讨论、反演算法实现等方面.正演快速算法主要讨论了等效存储几何格架技术、基于GPU加速的并行计算技术以及小波压缩技术.三维物性反演则是在最小二乘意义下使目标函数达到极小的线性或非线性反演.指出,对于特定地质问题需要谨慎选择不同且合适的约束方法乃至反演算法,才能达到好的效果.最后讨论了重磁数据三维物性反演较好的应用前景及发展方向.     

Using constrained inversion of gravity and magnetic field to produce a 3D litho‐prediction model

Geologically constrained inversion of gravity and magnetic field data of the Victoria property (located in Sudbury, Canada) was undertaken in order to update the present three‐dimensional geological model. The initial and reference model was constructed based on geological information from over 950 drillholes to constrain the inversion. In addition, downhole density and magnetic susceptibility measured in six holes were statistically analysed to derive lower and upper bounds on the physical properties attributed to the lithological units in the reference model. Constrained inversion of the ground gravity and the airborne magnetic data collected at the Victoria property were performed using GRAV3D and MAG3D, respectively. A neural network was trained to predict lithological units from the physical properties measured in six holes. Then, the trained network was applied on the three‐dimensional distribution of physical properties derived from the inversion models to produce a three‐dimensional litho‐prediction model. Some of the features evident in the lithological model are remnants of the constraints, where the data did not demand a significant change in the model from the initial constraining model (e.g., the thin pair of diabase dykes). However, some important changes away from the initial model are evident; for example, a larger body was predicted for quartz diorite, which may be related to the prospective offset dykes; a new zone was predicted as sulfide, which may represent potential mineralisation; and a geophysical subcategory of metabasalt was identified with high magnetic susceptibility and high density. The litho‐prediction model agrees with the geological expectation for the three‐dimensional structure at Victoria and is consistent with the geophysical data, which results in a more holistic understanding of the subsurface lithology.     

Pseudo-3D constrained inversion of transient electromagnetic data for a polarizable SMS hydrothermal system in the Deep Sea

Seafloor massive sulfide (SMS) deposits are generated by high-temperature hydrothermal systems. Their precious resources have attracted global interest. A number of investigations with controlled-source electromagnetic (CSEM) methods have been implemented in recent years. There are three major problems with SMS surveying using EM methods. First, SMS imaging techniques for hydrothermal systems have a limited range. Simulations and applications have validated only simple layered models. Second, their inversion efficiencies must be improved further. Laterally constrained inversions and spatially constrained inversions are usually used to map geological structures. However, choosing their suitable weighting parameters is inefficient. Third, the effects of induced polarization (IP) on ore deposits are not considered in such inversions. A non-polarizable model is unable to accurately depict a polarizable model. To resolve these problems, an advanced strategy is used to improve the efficiency of the pseudo-3D inversion process. The proposed imaging method has the ability to map complex 3D geoelectrical structures, and therefore, it can both obtain information regarding surface ore deposits and distinguish between active and inactive hydrothermal systems. However, this method can also be used to depict the distributions of alteration zones and buried deposits. Furthermore, the influences of IP on the inversion are discussed with respect to the Cole- Cole model, and it is shown that the effects of IP on polarizable deposits cannot be ignored during the inversion.     

A fracture network model and open fracture analysis of a tight sandstone gas reservoir in Dongpu Depression, Bohaiwan Basin, eastern China

Fracture systems can significantly influence productivity in tight sandstone reservoirs. Understanding the in situ fracture network is of importance in exploration and development of such reservoirs. This paper illustrates that stress-strain analysis, based on geological processes, is a credible prediction method for fracture networks and enables an estimation of which fractures are likely to contribute most to flow based on fracture orientation. This study shows the workflow for fracture modelling and its application in a sandstone reservoir using structural validation tools and interactive fracture simulations constrained with lithology, well log and core data. The fracture network model can help us understand not only what a fractured reservoir looks like but also how the existing fracture network behaves. An important consideration is that we need to understand and constrain the geological history to assess the fracture network properly.     

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.     

Three-dimensional geophysical characterization of the La Rambla and Zafra de Záncara anticlines (Loranca Basin,Central Spain)

The Zafra de Záncara anticline (also known as the El Hito anticline), located in the Loranca Cenozoic Basin (part of the Tagus Basin, Central Spain), had been studied by several oil companies during the late 1960s and early 1970s. In 2009, within the ‘Plan for selection and characterization of suitable structures of CO₂ geological storage’, this anticline was selected as a potential CO₂ storage site. A preliminary three-dimensional geological model, based on five geological cross sections that were constrained with the interpretation of the available seismic profiles (that are rather old and do not have very good quality), was created. With the aim of improving the geological knowledge of the Zafra de Záncara anticline and helping to investigate the suitability of a nearby anticline, namely La Rambla, as another structural closure that might make it a possible CO₂ storage site, a local gravity survey (1 station every km²) was carried out in the area, seven new geological cross sections, based on these existing seismic profiles and field geology, were build, and a new three-dimensional geological model that included both anticlines, improved through three-dimensional stochastic gravity inversion, was constructed. The densities needed for the geological formations of the model come from the analysis of rock samples, logging data from El Hito-1 drillhole and petrophysical information from Instituto Geológico y Minero de España database. The inversion has improved the knowledge about the geometry of the anticlines’ traps and seals as well as the geometry of the basement relief and the structural relationship between basement and cover.     

Fast and dynamic generation of linear octrees for geological bodies under hardware acceleration

In the application of 3D Geoscience Modeling,we often need to generate the volumetric representations of geological bodies from their surface representations.Linear octree,as an efficient and easily operated volumetric model,is widely used in 3D Geoscience Modeling.This paper proposes an algorithm for fast and dynamic generation of linear octrees of geological bodies from their surface models under hardware acceleration.The Z-buffers are used to determine the attributes of octants and voxels in a fast way,a...     

Mapping deep targets based on integrated 3D MT-gravity interpretation: a case study

Mapping deep geological hydrocarbon targets is of significant importance in basin exploration. In areas lacking reliable seismic data, magnetotelluric (MT) and gravity explorations are helpful to delineate the distribution of potential deep geological hydrocarbon targets. Here we investigate the effectiveness of the integrated 3D MT and gravity explorations for mapping the potential deep hydrocarbon source rocks. The result based on the data from the W Basin (part of the Ordes Basin) of China demonstrates that the method is efficient and economical for basin exploration. The method is particularly useful in target areas which are of great interest for oil and gas exploration but lack high quality seismic data. In our method, we first use the high-precision 3D small-bin MT data acquisition to improve the data accuracy. Then we perform datum static correction method and apply 3D inversion to obtain the3D resistivity distribution. We also develop a layered resistivity model based on resistivity logging to assist the interpretation of the inverted 3D resistivity data so as to derive an initial 3D geological model. Starting from the initial model, we use 2D gravity data to update the model via 2D inversion line by line, and then pass the updated model for the next round of the 3D MT inversion. The integrated inversion is implemented iteratively so the model converges to satisfy the need of final geological analysis. The application to the W Basin shows that we could successfully delineate the geological distribution of the potential deep hydrocarbon source rocks within the basin and map the thickness of the upper Paleozoic.