Groundwater Modeling in Integrated Water Resources Management—Visions for 2020

Groundwater modeling is undergoing a change from traditional stand-alone studies toward being an integrated part of holistic water resources management procedures. This is illustrated by the development in Denmark, where comprehensive national databases for geologic borehole data, groundwater-related geophysical data, geologic models, as well as a national groundwater-surface water model have been established and integrated to support water management. This has enhanced the benefits of using groundwater models. Based on insight gained from this Danish experience, a scientifically realistic scenario for the use of groundwater modeling in 2020 has been developed, in which groundwater models will be a part of sophisticated databases and modeling systems. The databases and numerical models will be seamlessly integrated, and the tasks of monitoring and modeling will be merged. Numerical models for atmospheric, surface water, and groundwater processes will be coupled in one integrated modeling system that can operate at a wide range of spatial scales. Furthermore, the management systems will be constructed with a focus on building credibility of model and data use among all stakeholders and on facilitating a learning process whereby data and models, as well as stakeholders' understanding of the system, are updated to currently available information. The key scientific challenges for achieving this are (1) developing new methodologies for integration of statistical and qualitative uncertainty; (2) mapping geological heterogeneity and developing scaling methodologies; (3) developing coupled model codes; and (4) developing integrated information systems, including quality assurance and uncertainty information that facilitate active stakeholder involvement and learning.     

三维非均匀地质模型中的逐段迭代射线追踪

地震射线追踪是地震定位、层析成像、偏移等领域的重要正演环节.随着这些领域研究的深入,针对传统的网格结构和层状结构在描述复杂地质模型遇到的很大困难,我们采用大小不等、形状各异的地质块组成的集合体来描述三维复杂地质模型,并用三角形面片来描述地质块之间的物性间断面,理论上可以描述任意复杂的地质模型.为适应任意非均匀速度分布的地质模型,基于费马原理,本文发展了与之相适应的逐段迭代射线追踪方法.该方法属于弯曲法范畴,对路径点采用一阶显式增量修正,相对于传统的迭代法,高效省时.数值试验表明,联合逐段迭代法和伪弯曲法的射线追踪扰动修正方案在三维复杂非均匀块状模型中有适用性和高效性.     

基于计算岩石物理方法的页岩储层弹性参数提取

发展了应用数值计算方法获取页岩储层的速度、各向异性参数的计算岩石物理系列方法.该系列方法包括了大尺度精细地质模型数值建模、计算网格尺度的地球物理建模和地震波数值模拟提取岩石物理弹性参数.本文方法利用储层的统计数据而不是具体岩心的测量数据,可获得储层岩石物理弹性参数的变化规律.相比于基于岩心测试的岩石物理方法,本文方法可精细考虑实际储层的非均匀特征,可得到岩心测试难以求取的与尺寸效应高度相关的弹性参数,也避免了求取弹性参数变化规律时获取不同地质特征岩心的困难.本文发展了计算岩石物理方法,为计算岩石物理面临的大尺度地质建模和计算能力限制问题提供了有效的解决方案.文中以胜利罗家的页岩储层为例,求得了储层TOC含量从3%到21%变化情况下储层的P波、S波速度以及各向异性参数变化规律.     

Integration of ground-penetrating radar and microgravimetric methods to map shallow caves

Ground-penetrating radar (GPR) and microgravimetric surveys have been conducted in the southern Jura mountains of western Switzerland in order to map subsurface karstic features. The study site, La Grande Rolaz cave, is an extensive system in which many portions have been mapped. By using small station spacing and careful processing for the geophysical data, and by modeling these data with topographic information from within the cave, accurate interpretations have been achieved.The constraints on the interpreted geologic models are better when combining the geophysical methods than when using only one of the methods, despite the general limitations of two-dimensional (2D) profiling. For example, microgravimetry can complement GPR methods for accurately delineating a shallow cave section approximately 10×10 m² in size. Conversely, GPR methods can be complementary in determining cavity depths and in verifying the presence of off-line features and numerous areas of small cavities and fractures, which may be difficult to resolve in microgravimetric data.     

航空电磁勘查技术发展现状及展望

航空电磁作为一种高效的地球物理勘查技术手段,其发展在国外(加拿大、澳大利亚等国家)已十分成熟.然而,在我国该项技术仍处于发展当中,在国内目前尚未形成具有实际探测能力的航空电磁系统和解释手段.这一现状严重影响了我国对地形地质条件复杂区域(比如广大西部地区)矿产资源勘查的需求.本文旨在通过系统介绍航空电磁勘查技术中的基础理论、关键技术、仪器系统、数据处理、解释及应用,并对未来我国航空电磁勘查技术的发展提出建议,使读者了解该技术未来发展方向和研究热点,以期该项技术在我国得到快速发展并获得广泛应用.     

复杂构造建模技术综述

目前复杂构造建模技术是地震资料处理中的技术瓶颈.本文主要介绍了目前国内外常用的复杂构造建模技术,包括层位标定技术、地层追踪技术、面块切片技术,以及相干体、方差体和三维可视化技术等,并分析了各自的适应性.     

A methodology for 3D modeling and visualization of geological objects

Geological body structure is the product of the geological evolution in the time dimension, which is presented in 3D configuration in the natural world. However, many geologists still record and process their geological data using the 2D or 1D pattern, which results in the loss of a large quantity of spatial data. One of the reasons is that the current methods have limitations on how to express underground geological objects. To analyze and interpret geological models, we present a layer data model to organize different kinds of geological datasets. The data model implemented the unification expression and storage of geological data and geometric models. In addition, it is a method for visualizing large-scaled geological datasets through building multi-resolution geological models rapidly, which can meet the demand of the operation, analysis, and interpretation of 3D geological objects. It proves that our methodology is competent for 3D modeling and self-adaptive visualization of large geological objects and it is a good way to solve the problem of integration and share of geological spatial data. Supported by National High Technology Research and Development Program of China (Grant Nos. 2006AA12Z220, 2006AA12Z114, 2007AA12Z226), and Open Fund of State Key Laboratory of Information Engineering in Surveying, Mapping and Remote Sensing (Grant No. WKL(06)0304)     

APPLICATION OF 3D HIGH PRECISION GRAVITY METHOD TO THE EXPLORATION OF DEEP BURIED HILL

Due to the influence of volcanic rock or weathering crust coverage, the quality of deep seismic data is often not reliable and it is difficult to identify small deep buried hill structure beneath the cover. However, the gravity exploration technology can use the remarkable density differences between the object and the overlying strata to identify this special geological reservoir. Although recently several rounds of seismic exploration have been carried out in Raoyang depression, Huabei Oilfield, North China to determine the existence of Hu 8 north deep buried hill, whether the buried hill really exists or not as well as its scale is still in controversial. In this paper, based on the information of seismic data and formation density, deep processing of 3D high precision ground gravity data has been carried out for 3D forward and inversion computation. The dimensional gravity forward calculation results show that the density of the deep anomaly body forms a relatively low gravity anomaly in the earth's surface. By using the potential processing method of vertical second derivative, and sliding filtering, the residual anomaly is separated from regional or background field, which illustrates the existence of a local high gravity anomaly at Hu 8 north area. According to the amplitude of residual gravity anomaly and formation density difference modeling, through a number of 3D forward calculations and 3D inversion of gravity data, the existence of Hu 8 north buried hill and its possible scale are basically determined. The results prove that the 3D high precision gravity method is effective in determining the deep buried hill structure in case that seismic data is not reliable.     

地下工程灾害水源的磁共振探测研究

地面磁共振探测(Surface Magnetic Resonance Sounding,SMRS)是近年来发展起来的一种地球物理新方法,这种在地面直接探测地下介质中氢核丰度的技术,不仅可以用于缺水地区的地下水资源勘查与评价,还可以在地下水引起的堤坝渗漏、滑坡、海水入侵等地质灾害水源的探测预警中发挥独特的作用.本文首次提出了地下磁共振探测方法(Underground Magnetic Resonance Sounding,UMRS),将SMRS方法引入到地下工程领域,实现隧道工程和煤矿开采等地下狭窄空间极端环境的探测.为应用UMRS方法,需要深入研究地下水超前探测理论、准全空间处理与反演方法、旋转多匝小线圈探测模式,强电磁干扰环境自适应噪声压制策略、以及复杂地质环境磁共振与瞬变电磁联合探测关键技术等难题.论文还简要介绍了超导磁探测技术和工程盾构及掘进实时探测等新技术在地下工程生产安全探测预警中的应用前景.     

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

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

测井技术在地震目标追踪应用中的方法探讨

目前常见测井-地震综合应用技术主要是测井约束地震反演技术,该技术存在两方面的研究精度问题.一是初始模型等因素影响反演的精度和分辨率,二是波阻抗的多解性难免造成反演的多解性.由于存在以上原因,有必要探索其他方式的"井震结合"研究思路,根据测井、地震信息具有地质背景成因一致的特点,提出寻找二者具有的、可相互辨识的共同地质属性的分析方法,研究认为,测井、地震共同地质属性的辨识可表现在两方面,一是测井、地震信息在地质界面上下地质成因相同.深入研究两者可相互辨识的地质属性信息,可提供地质界面的识别与追踪依据;二是对于一个具体的地质事件或地质体,也理应找到具有同一成因基础的、可相互辨识的信息响应依据.利用该方法对羊二庄某区的实践表明,在测井、地震信息成因一致性分析手段的基础上,开展测井-地质共同地质属性的辨识、追踪研究,对复杂油气区有利目标的追踪具有有效性,该研究成果被钻井结果证实.     

大型复杂构造地震物理模型设计制作及实验精度分析(英文)

大型、构造复杂的三维物理模型可用于模拟油气勘探。构造逼近实际地质状况的模拟具有制作技术难度大、质量控制严格等特点,可用于采集宽方位、多方位和全方位的地震数据,从而进行多种三维处理、解释方法验证。本文针对中国西部前陆盆地地表条件复杂地下构造复杂,导致成像不理想等问题,基于复杂的地下构造,设计制作了目前世界上模拟施工面积最大、构造最复杂的KS(塔里木盆地克深勘探工区)物理模型。本文的模型技术的进步主要涉及3个方面:模型的设计方法、模型的浇铸流程和数据采集,首次给出了物理模型的三维真实速度模型,定量分析了物理模型的制作精度,绝对误差小于3mm,可以满足方法试验的需要。该模型基于三维形态测量技术建立了三维真实速度模型,可作为方法试验的基础数据。因此,该模型可作为地震物理模拟技术的标准。     

地下电磁波法多重信息平面的场值分布模式

地下电磁波法的解释工作，原用的直射线代数重建层析方法，从理论分析即能看出。存在着对小目标反映不明显的局限性．另一方面，实际观测到的场值是叠加场，包含目标的散射信息，在多重情况下有明显的规律，可以表示成一定的模式．借助正演，反映出这些场值的规律，这是建立模型模式的过程，它为一种新的解释方法──模式识别，打下了基础．     

Multiple-point geostatistical simulation of dykes: application at Sungun porphyry copper system,Iran

Variogram-based methods are not capable of capturing high (>2) order statistics since the variogram measures the relationship between two points at a time only. Multiple-point geostatistics (MPS) has brought new insights into many geological modeling problems. The application of MPS methods has been well documented in realizing complex geological patterns. These methods have often been used in reservoir characterization since their advent in recent decades. The frequent non-linear behaviors of geologic continuity are not limited to reservoirs, but mineral deposits bear complicated formations in many cases. Relying on the power of MPS methods and considering the complexity of geological scenarios in mineral deposits, we have applied MPS in the modeling of mineral deposits. A training image (TI) is produced using geological data from upper horizons of a porphyry copper ore deposit which have been mined out during the previous mining operations. In this study, the SNESIM algorithm has been used. A number of realizations are produced using this multiple-point geostatistical method. Extensive validation steps are performed considering the TI as the reference model. These validations first show that the TI is representative for the domain under study and also illustrates some degrees of similarity between the TI and the realizations. Despite simplifications made to the problem, the application of MPS in mineral deposit modeling still faces many challenges.     

Numerical modeling of zero-offset laboratory data in a strong topographic environment: results for a spectral-element method and a discretized Kirchhoff integral method

Accurate simulation of seismic wave propagation in complex geological structures is of particular interest nowadays. However conventional methods may fail to simulate realistic wavefields in environments with great and rapid structural changes, due for instance to the presence of shadow zones, diffractions and/or edge effects. Different methods, developed to improve seismic modeling, are typically tested on synthetic configurations against analytical solutions for simple canonical problems or reference methods, or via direct comparison with real data acquired in situ. Such approaches have limitations, especially if the propagation occurs in a complex environment with strong-contrast reflectors and surface irregularities, as it can be difficult to determine the method which gives the best approximation of the “real” solution, or to interpret the results obtained without an a priori knowledge of the geologic environment. An alternative approach for seismics consists in comparing the synthetic data with high-quality data collected in laboratory experiments under controlled conditions for a known configuration. In contrast with numerical experiments, laboratory data possess many of the characteristics of field data, as real waves propagate through models with no numerical approximations. We thus present a comparison of laboratory-scaled measurements of 3D zero-offset wave reflection of broadband pulses from a strong topographic environment immersed in a water tank with numerical data simulated by means of a spectral-element method and a discretized Kirchhoff integral method. The results indicate a good quantitative fit in terms of time arrivals and acceptable fit in amplitudes for all datasets.     

基于Kriging方法的三维地质可视化系统构建

针对传统三维地质建模方法存在的诸多问题,对基于Kriging方法的三维地质体构建进行研究。选取某地震工程项目钻孔数据为主要研究对象,在Arc GIS平台上通过地统计分析,实现了复杂三维地质体的快速、准确建立。在此基础上,采用Arc Engine作为二次开发组件,完成该地质体三维展示系统的开发,实现三维展示、专题图出图、三维分析等功能。     

中海油多分量地震处理技术的发展与应用实例

20世纪末中海油在南海油气勘探中成功试验了二维多分量地震技术.其后十多年,中海油进一步在不同近海油气盆地采集了三维多分量OBC地震资料.但由于针对三维转换波的处理未能形成有效的关键技术,除双检叠加利用了水压分量和陆检垂直分量外,大量水平分量数据未被使用.鉴于此,中海油在十二五期间开展三维多分量地震数据处理的关键技术攻关：针对弹性波场的矢量特征、转换波射线路径不对称以及双程旅行时多时间尺度等问题,成功研发了矢量化的信号处理技术、弹性波速度建模以及叠前时间偏移成像等技术;集成了相对完善的三维弹性波成像EWI软件系统;完成了多个试验区的三维四分量OBC地震转换波的处理,取得了好于以往技术的效果;建立了海上多分量OBC地震资料处理流程.但相对于成熟的纵波处理技术,海洋多分量地震处理技术仍需要不断完善与发展.     

地震数据叠前偏移成像及其在浅层地震勘探中的应用研究

文中讨论了在浅层地震勘探中地质结构和构造的复杂性以及对浅层地震勘探成果要求精度较高的特殊性,指出了在浅层地震勘探中发展叠前深度偏移的重要性,并对相移法叠前深度偏移的基本原理和实现方法进行了详细论述,编制了相应的软件。通过使用该方法对 2种理论试验模型进行试算,获得的深度偏移剖面层位清晰,构造明显,较准确地体现了理论模型地质结构和构造的特征,效果良好。然后又使用实际浅层地震反射资料进行了叠前偏移计算,所得结果与已知地质结构和构造相吻合,也取得了较好的结果,证实了该方法及软件的可靠性和实用性,可用于实际浅层地震反射勘探资料的叠前深度偏移处理     

Three-dimensional geological modeling and its application in Digital Mine

The 3D geological modeling is the prerequisite and core foundation for Digital Mine. Although this new technology brings new opportunities and motivation for the mineral exploration industry, it still has many difficulties to be solved in this area. Based on the characteristics of mine data and the aim of Digital Mine construction, this paper introduces a theory including multi-source data coupling, multi-modeling methods integration, multi-resolution visualization and detection, and multidimensional data analysis and application. By analyzing problems such as the uncertainty in each step of the modeling process, we designed a novel modeling method that can be applied to the complex geological body modeling, mineral resource/reserve estimation, and the mining exploration engineering. Along with the process of mine exploration, development, and reclamation, 3D modeling undergoes the process of “construction-simulation-revision” during which the 3D model is able to be dynamically updated and gradually improved. Based on the result of practical utilization, it is proven that the methodology introduced by this paper can be used to build an effective 3D model by fully using the mining data under the control of spatial information quality evaluation. Our experiments show that such a 3D model can be used to evaluate the mine resource and provide the scientific evidence to improve mining efficiency during the various stages of evolvement process in mine.     

Simulation of the hydrocarbon structures using the P-SV wave solution

Numerical solution of the scalar and elastic wave equations has greatly aided geophysicists in modeling seismic wave fields in the complicated geologic structures containing hydrocarbons and hence increases the geologic interpretation. Finite-difference method offers a versatile approach to compute synthetic seismograms numerically for given subsurface complex geological structures. To avoid the spatial derivative of the elastic parameters and density, elastodynamic equation (first-order hyperbolic equation) has been solved using the Lax-Wendroff scheme. A numerical finite-difference modeling program has been developed for the P-SV wave using the above solution. A line source with a time delay of 0.015s and dominant frequency of 120 Hz has been utilized in the simulation. In order to avoid the large values of the displacement vector in the source region,Alterman andKaral's method (1968) has been utilized. Horizontal and vertical component synthetic seismograms have been computed for two different geological models with and without oil and gas bearing zones. It has been concluded from the response that a finite-difference technique not only yields the relative arrival times but also accounts for the variation in amplitude and phase according to the elastic impedance contrast across the interfaces. It should come as no surprise to learn that in spite of the limitation of this numerical method, the scheme has provided a valid response for the thin layer, high acoustic impedance contrast and the pinch out.