测井约束地震反演综述

阐述了测井约束地震反演研究的必要性、技术内涵及研究思路，回顾了测井约束地震反演技术的发晨历程，比较了三种非线性测井约束地震反演方法的优缺点。小波变换以其独特的变时窗特性，在测井约束地震反演技术发展中潜力较大，值得深入研究。同时对测井约束地震反演技术的研究现状做了小结，指出当前的研究对象已转向了复杂地质体。特别是海相碳酸盐岩复杂地质体，由于海相和陆相油气田存在许多不同之处，所以测井约束地震反演技术将面临巨大的机遇和挑战。最后，指出了测井约束地震反演技术的发展方向。     

基于FCM聚类算法的二维RMT电阻率与介电常数联合反演

射频大地电磁(RMT)是新兴的浅地表勘探法, 其工作频段(10 kHz~1 MHz)处于音频大地电磁(AMT)与探地雷达(GPR)之间, 随频率的增加, 位移电流对RMT观测的影响逐渐增大, 单一电阻率参数反演会影响数据解释的准确性.本文在同时考虑传导电流与位移电流影响的条件下, 开展二维RMT数据电阻率和介电常数联合反演研究.为提高电阻率与介电常数之间的联合约束效果, 在经典最小结构模型正则化约束的基础上, 引入模糊C均值(FCM)聚类算法进行联合约束.模型试算表明, 基于FCM聚类的二维RMT电阻率与介电常数联合反演可以在丰富解释参数的同时提高反演效果.

     

三维电阻率测深数据Zohdy近似反演方法

Zohdy方法近似反演三维电阻率测深数据。正演计算采用有限单元法。反演初始模型由测量视电阻率数据给出。通过比较实测视电阻率值和预测模型计算的视电阻率值对数差来修改模型网格电阻率．为了解决任意电极距测深数据的反演,采用大、小双网格剖分。大网格反映地下电性分布情况。小网格用于实际有限元正演计算．在电阻率调整公式中加入一个迭代系数,能够加快收敛速度．并对加5％随机噪声的模型理论视电阻率测深断面数据进行反演,得到的电阻率分布与模型电阻率基本一致．     

Estimating thermal conductivity from core and well log data

The aim of the presented work was to introduce a method of estimating thermal conductivity using well log data. Many petrophysical properties of rocks can be determined both by laboratory measurements and well-logs. It is thus possible to apply geophysical data to empirical models based on relationships between laboratory measured parameters and derive continuous thermal conductivity values in well profiles. Laboratory measurements were conducted on 62 core samples of Meso-Paleozoic rocks from the Carpathian Foredeep. Mathematical models were derived using multiple regression and neural network methods. Geophysical data from a set of seven well logs: density, sonic, neutron, gamma ray, spectral gamma ray, caliper and resistivity were applied to the obtained models. Continuous thermal conductivity values were derived in three well profiles. Analysis of the obtained results shows good consistence between laboratory data and values predicted from well log data.     

测井约束地震波形反演的同伦摄动法

测井数据和地震数据是地震勘探中两种最重要的资料. 测井约束地震波形反演是在非线性波形反演的基础上,利用已知测井资料详细的垂直分辨能力和地震资料均匀密集的水平采样特点, 通过迭代反演来求取一个具有较高分辨率的速度参数.本文建立了测井约束反演模型,研究了测井约束下地震波形反演的同伦摄动求解方法.同伦摄动法作为一种新的、求解数学物理中各种非线性问题的有效方法,具有计算速度快、计算精度高的优点.这对于提高反演的精度和效率是十分有益的. 为了表征该方法的有效性和稳定性,分别对水平层状介质模型和逆冲断层带模型进行了数值模拟,并与Landweber迭代法相对比,结果表明该算法具有更好的收敛性,能够取得更为满意的反演效果.     

测井约束地震波形反演的同伦摄动法

测井数据和地震数据是地震勘探中两种最重要的资料. 测井约束地震波形反演是在非线性波形反演的基础上,利用已知测井资料详细的垂直分辨能力和地震资料均匀密集的水平采样特点, 通过迭代反演来求取一个具有较高分辨率的速度参数.本文建立了测井约束反演模型,研究了测井约束下地震波形反演的同伦摄动求解方法.同伦摄动法作为一种新的、求解数学物理中各种非线性问题的有效方法,具有计算速度快、计算精度高的优点.这对于提高反演的精度和效率是十分有益的. 为了表征该方法的有效性和稳定性,分别对水平层状介质模型和逆冲断层带模型进行了数值模拟,并与Landweber迭代法相对比,结果表明该算法具有更好的收敛性,能够取得更为满意的反演效果.     

Two-dimensional joint inversion of radiomagnetotelluric and direct current resistivity data

An algorithm for the two-dimensional (2D) joint inversion of radiomagnetotelluric and direct current resistivity data was developed. This algorithm can be used for the 2D inversion of apparent resistivity data sets collected by multi-electrode direct current resistivity systems for various classical electrode arrays (Wenner, Schlumberger, dipole-diplole, pole-dipole) and radiomagnetotelluric measurements jointly. We use a finite difference technique to solve the Helmoltz and Poisson equations for radiomagnetotelluric and direct current resistivity methods respectively. A regularized inversion with a smoothness constrained stabilizer was employed to invert both data sets. The radiomagnetotelluric method is not particularly sensitive when attempting to resolve near-surface resistivity blocks because it uses a limited range of frequencies. On the other hand, the direct current resistivity method can resolve these near-surface blocks with relatively greater accuracy. Initially, individual and joint inversions of synthetic radiomagnetotelluric and direct current resistivity data were compared and we demonstrated that the joint inversion result based on this synthetic data simulates the real model more accurately than the inversion results of each individual method. The developed 2D joint inversion algorithm was also applied on a field data set observed across an active fault located close to the city of Kerpen in Germany. The location and depth of this fault were successfully determined by the 2D joint inversion of the radiomagnetotelluric and direct current resistivity data. This inversion result from the field data further validated the synthetic data inversion results.     

测井数据岩石物理参数敏感性的应用研究

利用测井资料所获得的基本参数进行岩石物理参数敏感性分析,可以较大提高含烃岩石的识别能力。本文讨论了对地震属性敏感性高的参数的选择及其使用方法,以便提高流体和岩性的识别能力。为了更好的与地震数据解释相结合,我们在敏感性参数空间构造了研究区的岩石物理模板,用以指导实际地震数据定量解释。     

Piecewise 1D laterally constrained inversion of resistivity data

In a sedimentary environment, layered models are often capable of representing the actual geology more accurately than smooth minimum structure models. Furthermore, interval thicknesses and resistivities are often the parameters to which non‐geophysicist experts can relate and base decisions on when using them in waste site remediation, groundwater modelling and physical planning. We present a laterally constrained inversion scheme for continuous resistivity data based on a layered earth model (1D). All 1D data sets and models are inverted as one system, producing layered sections with lateral smooth transitions. The models are regularized through laterally equal constraints that tie interface depths and resistivities of adjacent layers. Prior information, e.g. originating from electric logs, migrates through the lateral constraints to the adjacent models, making resolution of equivalences possible to some extent. Information from areas with well‐resolved parameters will migrate through the constraints in a similar way to help resolve the poorly constrained parameters. The estimated model is complemented by a full sensitivity analysis of the model parameters, supporting quantitative evaluation of the inversion result. Examples from synthetic 2D models show that the model recognition of a sublayered 2D wedge model is improved using the laterally constrained inversion approach when compared with a section of combined 1D models and when compared with a 2D minimum structure inversion. Case histories with data from two different continuous DC systems support the conclusions drawn from the synthetic example.     

测井约束反演过程中测井资料统计分析研究

本文通过对测井、地质等资料的分析,选择分布稳定、岩性单一的地层作为标准化处理的标准地层.通过对声波曲线进行误差校正和重构,使储层的速度和波阻抗特征更加明显,以此来区分砂岩与泥岩,提高识别储层的精确度.同时,通过对测井资料的统计分析使反演剖面的横向变化比较均衡,而且与地下地层岩性的实际变化更加吻合,解决了井间能量分布不均匀这一问题.     

Shear wave velocity prediction using seismic attributes and well log data

Formation’s properties can be estimated indirectly using joint analysis of compressional and shear wave velocities. Shear wave data is not usually acquired during well logging, which is most likely for cost saving purposes. Even if shear data is available, the logging programs provide only sparsely sampled one-dimensional measurements: this information is inadequate to estimate reservoir rock properties. Thus, if the shear wave data can be obtained using seismic methods, the results can be used across the field to estimate reservoir properties. The aim of this paper is to use seismic attributes for prediction of shear wave velocity in a field located in southern part of Iran. Independent component analysis (ICA) was used to select the most relevant attributes to shear velocity data. Considering the nonlinear relationship between seismic attributes and shear wave velocity, multi-layer feed forward neural network was used for prediction of shear wave velocity and promising results were presented.     

基于印模法的地面磁电阻率数据三维非线性共轭梯度反演

本文提出了能提高异常体分辨能力,同时得到绝对电导率的地面磁电阻率数据三维反演方法.磁电阻率响应用准直流的低频磁场代替;数值模拟由频率域电场满足的Helmholtz方程出发,采用三维交错网格有限差分法;长直导线源作为发射源,其中源的计算包含在背景场中;结合地面磁电阻率数据各分量的特点,选择y分量进行反演研究;反演采用三维非线性共轭梯度反演技术,为了提高异常体的深度分辨能力,进行迭代重构反演;用印模法对初始模型进行重构,采用的是辅模型在浅部,元模型在深部的组合方式.从合成数据和实际数据的反演结果可以得到以下的认识：（1）由频率域麦克斯韦方程组出发,低频磁场数据反演可以直接得到电导率,而不是相对电导率之比;（2）采用印模法组合初始模型,进行迭代重构反演,可以提高地面磁电阻率数据反演对异常体的分辨能力,确定埋深位置,同时不会丧失对于浅部异常体的分辨能力;（3）在结合印模法的地面磁电阻率数据三维反演中,深部异常体的分辨能力受地表不均匀导电体影响较小;（4）确定印模深度可以采用上一次重构反演结束时的模型变化量,通过相邻两次重构反演结束时的模型变化量之差来确定迭代重构是否终止.因为静磁场与重力场在数学上的相似性,本文的反演方法可以被运用到重力场等位场的地面数据的反演中.

     

Practical techniques for 3D resistivity surveys and data inversion1

Techniques to reduce the time needed to carry out 3D resistivity surveys with a moderate number (25 to 100) of electrodes and the computing time required to interpret the data have been developed. The electrodes in a 3D survey are normally arranged in a square grid and the pole-pole array is used to make the potential measurements. The number of measurements required can be reduced to about one-third of the maximum possible number without seriously degrading the resolution of the resulting inversion model by making measurements along the horizontal, vertical and 45° diagonal rows of electrodes passing through the current electrode. The smoothness-constrained least-squares inversion method is used for the data interpretation. The computing time required by this technique can be greatly reduced by using a homogeneous half-space as the starting model so that the Jacobian matrix of partial derivatives can be calculated analytically. A quasi-Newton updating method is then used to estimate the partial derivatives for subsequent iterations. This inversion technique has been tested on synthetic and field data where a satisfactory model is obtained using a modest amount of computer time. On an 80486DX2/66 microcomputer, it takes about 20 minutes to invert the data from a 7 by 7 electrode survey grid. using the techniques described below, 3D resistivity surveys and data inversion can be carried out using commercially available field equipment and an inexpensive microcomputer.     

二维直流电阻率与射频大地电磁双参数聚类联合反演

考虑到射频大地电磁（RMT）法反演介电常数具有可行性,但目前RMT与直流电阻率（DC）联合反演研究中均只进行了电阻率参数的反演,本文基于模糊C均值（FCM）聚类算法,开展了RMT与DC的电阻率与介电常数联合反演.首先,检验了介电常数对RMT观测的影响,并通过DC与RMT法的电阻率参数反演,分析了联合反演相较于单一方法反演的优势;基于双数据集反演电阻率的优势,辅以FCM聚类的双参数联合约束,进一步实现最佳介电常数模型反演.理论试算表明,DC与RMT法的电阻率联合反演,可补足RMT浅部高阻分辨率不理想的缺点,同时改善DC法深部反演能力的不足;此外,通过FCM聚类的联合约束以及真实物性的引导,相较于单一RMT方法,DC与RMT的联合反演可恢复更准确的地下介电常数分布.

     

Two‐dimensional joint inversions of cross‐hole resistivity data and resolution analysis of combined arrays

In this study, a new two‐dimensional inversion algorithm was developed for the inversion of cross‐hole direct current resistivity measurements. In the last decades, various array optimisation methods were suggested for resistivity tomography. However, researchers have still collected data by using classical electrode arrays in most cross‐hole applications. Therefore, we investigated the accuracy of both the individual and the joint inversion of the classical cross‐hole arrays by using both synthetic and field data with the developed algorithm. We showed that the joint inversion of bipole–bipole, pole–bipole, bipole–pole, and pole–tripole electrode arrays gives inverse solutions that are closer to the real model than the individual inversions of the electrode array datasets for the synthetic data inversion. The model resolution matrix of the suggested arrays was used to analyse the inversion results. This model resolution analysis also showed the advantage of the joint inversion of bipole–bipole, pole–bipole, bipole–pole, and pole–tripole arrays. We also used sensitivity sections from each of the arrays and their superpositions to explain why joint inversion gives better resolution than the any individual inversion result.     

Deep resistivity structure of the Dikili-Bergama region, west Anatolia, revealed by two- dimensional inversion of vertical electrical sounding data

Within the framework of the National Marine Geological and Geophysical Program, we re‐examined deep vertical electrical sounding (VES) data. The data, measured in 1968 by the General Directorate of Mineral Research and Exploration (MTA) of Turkey with the aim of exploring the deep resistivity structure of the Dikili–Bergama region, focus on the geothermal potential. The geoelectrical resistivity survey was conducted using a Schlumberger array with a maximum electrode half‐spacing of 4.5 km. The two‐dimensional (2D) inversion was utilized to interpret the VES data that were collected along 15‐ to 30‐km profiles. The 2D resistivity–depth cross‐sections obtained show very low resistivity values near the Dikili and Kaynarca hot springs. The 2D inversion results also indicate the presence of fault zones striking nearly N–S and E–W, and fault‐bounded graben‐horst structures that show promising potential for geothermal field resources. The 2D gravity model, which is in good agreement with the density variation of the region, supports the resistivity structure revealed by 2D inversion. The lithology information obtained from the borehole near Kaynarca also confirms the results of the resistivity interpretation and the density model.     

基于测井和地震资料的神狐海域天然气水合物资源量估算

饱和度是估算资源量的一个重要参数,尤其是深水油气和天然气水合物共存的白云凹陷.本文利用神狐地区测井的速度、电阻率、氯离子浓度和地震反演的纵波速度异常,对该区水合物饱和度进行估算.利用二维高分辨率地震资料,分析了气烟囱、断裂、砂体、滑塌体等特殊构造对气体运移的控制作用和对水合物形成的作用.基于双相介质中纵波速度,利用迭代正演模拟反演法估算出神狐地区的天然气水合物饱和度占孔隙空间的5～30%,局部地区饱和度为44%.水合物中甲烷气体资源量达到千亿立方米,气体储量达到中等气藏的储量丰度.     

Three‐dimensional inversion of towed streamer electromagnetic data

A towed streamer electromagnetic system capable of simultaneous seismic and electromagnetic data acquisition has recently been developed and tested in the North Sea. We introduce a 3D inversion methodology for towed streamer electromagnetic data that includes a moving sensitivity domain. Our implementation is based on the 3D integral equation method for computing responses and Fréchet derivatives and uses the re‐weighted regularized conjugate gradient method for minimizing the objective functional with focusing regularization. We present two model studies relevant to hydrocarbon exploration in the North Sea. First, we demonstrate the ability of a towed electromagnetic system to detect and characterize the Harding field, a medium‐sized North Sea hydrocarbon target. We compare our 3D inversion of towed streamer electromagnetic data with 3D inversion of conventional marine controlled‐source electromagnetic data and observe few differences between the recovered models. Second, we demonstrate the ability of a towed streamer electromagnetic system to detect and characterize the Peon discovery, which is representative of an infrastructure‐led shallow gas play in the North Sea. We also present an actual case study for the 3D inversion of towed streamer electromagnetic data from the Troll field in the North Sea and demonstrate our ability to image all the Troll West Oil and Gas Provinces and the Troll East Gas Province. We conclude that 3D inversion of data from the current generation of towed streamer electromagnetic systems can adequately recover hydrocarbon‐bearing formations to depths of approximately 2 km. We note that by obviating the need for ocean‐bottom receivers, the towed streamer electromagnetic system enables electromagnetic data to be acquired over very large areas in frontier and mature basins for higher acquisition rates and relatively lower cost than conventional marine controlled‐source electromagnetic methods.     

Assessment of the reliability of 2D inversion of apparent resistivity data

The reliability of inversion of apparent resistivity pseudosection data to determine accurately the true resistivity distribution over 2D structures has been investigated, using a common inversion scheme based on a smoothness‐constrained non‐linear least‐squares optimization, for the Wenner array. This involved calculation of synthetic apparent resistivity pseudosection data, which were then inverted and the model estimated from the inversion was compared with the original 2D model. The models examined include (i) horizontal layering, (ii) a vertical fault, (iii) a low‐resistivity fill within a high‐resistivity basement, and (iv) an upfaulted basement block beneath a conductive overburden. Over vertical structures, the resistivity models obtained from inversion are usually much sharper than the measured data. However, the inverted resistivities can be smaller than the lowest, or greater than the highest, true model resistivity. The substantial reduction generally recorded in the data misfit during the least‐squares inversion of 2D apparent resistivity data is not always accompanied by any noticeable reduction in the model misfit. Conversely, the model misfit may, for all practical purposes, remain invariant for successive iterations. It can also increase with the iteration number, especially where the resistivity contrast at the bedrock interface exceeds a factor of about 10; in such instances, the optimum model estimated from inversion is attained at a very low iteration number. The largest model misfit is encountered in the zone adjacent to a contact where there is a large change in the resistivity contrast. It is concluded that smooth inversion can provide only an approximate guide to the true geometry and true formation resistivity.     

Combined straightforward inversion of resistivity and induced polarization (time-domain) sounding data

It is proposed that the Straightforward Inversion Scheme (SIS) developed by the authors for 1D inversion of resistivity sounding and magneto-telluric sounding data can also be used in similar fashion for time-domain induced polarization sounding data. The necessary formulations based on dynamic dipole theory are presented. It is shown that by using induced polarization potential, measured at the instant when steady state current is switched off, an equation can be developed for apparent ‘chargeability–resistivity’ which is similar to the one for apparent resistivity. The two data sets of apparent resistivity and apparent chargeability–resistivity can be inverted in a combined manner, using SIS for a common uniform thickness layer earth model to estimate the respective subsurface distributions of resistivity and chargeability–resistivity. The quotient of the two profiles will give the sought after chargeability profile. A brief outline of SIS is provided for completeness. Three theoretical models are included to confirm the efficacy of SIS software by inverting only the synthetic resistivity sounding data. Then one synthetic data set based on a geological model and three field data sets (combination of resistivity and IP soundings) from diverse geological and geographical regions are included as validation of the proposal. It is hoped that the proposed scheme would complement the resistivity interpretation with special reference to shaly sand formations.