Post-stack impedance blocky inversion based on analytic solution of viscous acoustic wave equation

The conventional impedance inversion method ignores the attenuation effect, transmission loss and inter-layer multiple waves; the smooth-like regularization approach makes the corresponding impedance solution excessively smooth. Both fundamentally limit the resolution of impedance result and lead to the inadequate ability of boundary characterization. Therefore, a post-stack impedance blocky inversion method based on the analytic solution of viscous acoustic equation is proposed. Based on the derived recursive formula of reflections, the 1D viscous acoustic wave equation is solved analytically to obtain zero-offset full-wave field response. Applying chain rule, the analytical expression of the Fréchet derivative is derived for gradient-descent non-linear inversion. Combined with smooth constraints, the blocky constraints can be introduced into the Bayesian inference framework to obtain stable and well-defined inversion results. According to the above theory, we firstly use model data to analyse the influence of incompleteness of forward method on seismic response, and further verify the effectiveness of the proposed method. Then the Q-value sensitivity analysis of seismic trace is carried out to reduce the difficulty of Q-value estimation. Finally, the real data from Lower Congo Basin in West Africa indicate that the proposed approach provide the high-resolution and well-defined impedance result. As a supplement and development of linear impedance inversion method, the non-linear viscous inversion could recover more realistic and reliable impedance profiles.     

基于三维正演的音频大地电磁阻抗相位不变量校正技术

复杂地形、地质条件的大地电磁数据解释容易出现假象,采用三维正演技术模拟地形和地表不均匀体的背景响应,对实测数据阻抗相位不变量进行校正,实现更准确的定性分析;对三维异常体模型的合成数据进行一维、二维多参数反演试算,以确定地形剧变区选择反演技术的最佳方案.合成数据的试反演结果显示一维反演水平切片假异常较多,二维反演能压制测向假异常,但不能压制走向的假异常,水平切片多出现测向条带.使用本文提出的阻抗相位不变量校正法扣除地形、地表背景响应,结合一维、二维反演,能使实际资料解释成果更加可靠.     

叠前反演和地震吸收技术在复杂天然气藏地震预测中的应用

针对某复杂断块天然气目标储层,在岩石物理分析的指导下,综合利用地质、地震、测井等资料,提出了一套面向复杂天然气藏的叠前地震预测技术.首先基于地震岩石物理分析得到的初始横波信息,采用叠前贝叶斯非线性三参数反演得到了井旁控制点处精确纵横波速度和密度信息,然后通过叠前/叠后联合反演技术实现了面向目标的弹性阻抗体反演及含气储层敏感参数直接提取,最后结合小波变换时频谱分析的方法从叠前地震资料中估算地层吸收参数值,提高天然气藏识别精度.实际应用表明,综合各种叠前地震预测技术,可以大大提高对复杂天然气藏的识别精度,降低勘探风险.     

火山地层结构特征及其对波阻抗反演的约束

火山地层是由似层状、层状和块状结构组成的复合体,其地层结构的拟合是火山岩储层地震预测的难点所在,火山地层结构模型的建立是火山岩勘探开发中首先需要解决的关键问题.通过对九台和大屯典型火山岩地层结构的解剖,分析似层状结构和层状结构特征.似层状结构地层内部的岩层倾向和倾角变化大,与地层顶底面表现为斜交.从火山机构喷发中心到远...     

Impedance matrices for circular foundation embedded in layered medium

A numerical scheme is developed in the paper for calculating torsional, vertical, horizontal, coupling and rocking impedances in frequency domain for axial-symmetric foundations embedded in layered media. In the scheme, the whole soil domain is divided into interior and exterior domains. For the exterior domain, the analytic solutions with unknown coefficients are obtained by solving three-dimensional (3D) wave equations in cylindrical coordinates satisfying homogeneous boundary conditions. For the interior domain, the analytical solutions are also obtained by solving the same 3D wave equations satisfying the homogeneous boundary conditions and the prescribed boundary conditions. The prescribed conditions are the interaction tractions at the interfaces between embedded foundation and surrounding soil. The interaction tractions are assumed to be piecewise linear. The piecewise linear tractions at the bottom surface of foundation will be decomposed into a series of Bessel functions which can be easily fitted into the general solutions of wave equations in cylindrical coordinates. After all the analytic solutions with unknown coefficients for both interior and exterior domains are found, the variational principle is employed using the continuity conditions (both displacements and stresses) at the interfaces between interior and exterior domains, interior domain and foundation, and exterior domain and foundation to find impedance functions.     

利用岩性阻抗识别砂岩、泥岩以及砂体叠置——以北部湾盆地涠洲A油田为例

针对涠洲A油田储层平面非均质性强、存在砂体叠置,以及单纯的波阻抗反演不能十分有效地识别砂泥岩和叠置砂体,提出了将纵波阻抗与横波阻抗通过一定的数学变换获得岩性阻抗的方法,证实了该岩性阻抗反演结果与实钻井吻合程度高,有效区分了砂岩和泥岩,并成功地预测了涠洲A油田砂体叠置的关系。     

长白山天池火山区大地电磁测深初步观测

在长白山火山区沿着两条剖面进行了１５个测点的大地电磁测深观测，对地磁感应矢量和地下电性分布维数特征进行了分析，并使用阻抗张量分解和二维自动反演等新技术对资料进行了解释。结果表明，火山区及其周围地区的地下电性结构具有二维特征，但在火山口及其附近地下存在局部三维异常体，异常体埋深约１２ｋｍ，初步推测可能是火山岩浆囊     

表面局部三维大地电磁曲线畸变校正——MT畸变校正阻抗张量分解方法

针对三维／二维地质结构开发了大地电磁阻抗张量分解技术,利用该技术仅需对各个独立测点的大地电磁资料进行分析即可分离和消除由三维局部异常体引起的其附近测点的ＭＴ响应函数的形态扭曲畸变。     

大地电磁(MT)阻抗张量的正则分解及其初步应用

鉴于大地电磁阻抗张量传统分析方法的局限性,提出了普遍适用的正则分解方法。在用正则分解方法分析时,主坐标系提供的８个标量参数,描述了地球系统传输特性和场源极化状态及其联系,因此是物理意义明确、易于解释的参数体系。讨论了正则分解方法在各种构造条件下的特点和场源的极化问题,并给出了简单的应用实例     

Assessing the Validity of Seismo-Acoustic Predictor Equations for Obtaining Seabed and Sub-Surface Sediment Physical Properties

The interdependence between the seismo-acoustic properties of a marine sediment and its geotechnical/physical parameters has been known for many years, and it has been postulated that this should allow the extraction of geotechnical information from seismic data. Though in the literature many correlations have been published for the surficial layer, there is a lack of information for greater sediment depths. In this article, a desktop study on a synthetic seafloor model illustrates how the application of published near-surface prediction equations to subsurface sediments (up to several tens of meters burial depth) can lead to spurious predictions. To test this further, acoustic and geotechnical properties were measured on a number of sediment core samples, some of which were subjected to loading in acoustically-equipped consolidation cells (oedometers) to simulate greater burial depth conditions. For low effective pressures (representing small burial depths extending to around 10 meters subsurface), the general applicability of established relationships was confirmed: the prediction of porosity, bulk density, and mean grain size from acoustic velocity and impedance appears generally possible for the investigated sedimentary environments. As effective pressure increases through, the observed relationships deviate more and more from the established ones for the near-surface area. For the samples tested in this study, in some instances increasing pressure even resulted in decreasing velocities. There are several possible explanations for this abnormal behavior, including the presence of gas, overconsolidation, or bimodal grain size distribution. The results indicate that an appropriate depth correction must be introduced into the published prediction equations in order to obtain reliable estimates of physical sediment properties for greater subsurface depths.