地震数据倾角约束的多道谱反演算法

反射系数反演是联结地下储层和地震数据的桥梁,奇偶分解算法的出现使得子波间调谐效应减弱,这使得基于压缩感知的谱反演得到进一步应用．由于谱反演算法的不稳定性,所得到的反射系数横向连续性较差．因此,提出叠后地震数据倾角约束的多道谱反演算法,算法认为地震数据沿倾角方向具有一定连续性,在常规单道谱反演的基础上,推导了多道谱反演算法,基于局部倾角增加沿地层倾向的平滑约束,解决大角度地层反演横向连续性差的问题．算法继承了谱反演的高分辨率特性,并且有效增强了横向连续性,适用于地震数据的反射系数反演．模型和实际数据测试证明了倾角约束的多道谱反演算法得到的反演结果不仅能识别薄层,还能保持原始地层模型的横向连续性特征,并且具有一定的抗噪性,为地震地层学精细解释提供依据．     

基于波动方程的叠前地震反演

随着地震勘探和开发的不断深入,面向地质目标的精细储层预测技术变得越来越重要.由于透射损失、层间多次波、波模式转换以及随机噪声等的影响,观测地震数据和待反演的地下介质属性之间呈现出很强的非线性.考虑到这些非线性,本文基于积分波动方程开展叠前地震反演,从观测地震数据中恢复出介质属性和整体波场,其中反演参数是波动方程中的压缩系数、剪切柔度和密度的对比度,相比于常规线性AVO反演的波阻抗弹性参数,它们对流体指示有更强的敏感性.在反演过程中,从平滑的低频背景场出发,交替迭代求解数据方程和目标方程.采用乘性正则化方法于共轭梯度框架下求解反演参数,采用优化的散射级数Neumann序列获得整体波场,这种方法不易陷入局部极值,能收敛到正确解.测井资料和典型山前带模型测试表明,利用上述反演方法能获得高分辨率的深度域地下介质属性,可直接进行储层预测和解释.     

地震数据互相关驱动的多道反演方法

地震反演是储层定量描述和地震油气识别的关键技术,反演结果在复杂构造区域的横向连续性和保真性是影响地震资料定量解释精度的重要因素.基于此,本文发展了地震数据互相关驱动的多道反演方法.考虑地层反射系数与地震数据在结构上具有相似性的特点,基于地震数据互相关描述地层反射系数的结构特征,并将其作为多道地震反演的横向约束条件;此外,为改善地震数据本身横向连续性差对反演结果的影响,在目标泛函的惩罚项中引入局部优化算子,构建了一个易于求解的多道地震反演目标泛函.与常规多道地震反演方法相比,本文方法能够设计更合理、更符合实际情况的横向约束算子,提高反演结果的横向连续性,并且能有效降低地震资料质量对反演结果的影响.模型测试和实际应用验证了本方法的可靠性和稳定性.     

利用地震数据反演海水温盐结构

利用地震剖面获取海水层温度、盐度、密度等物理参数成为地震海洋学研究的一个重要问题.本文提出了以CTD(Conductivity-Temperature-Depth)温盐深剖面仪观测资料为约束的波阻抗、温-盐结构反演方法.该方法包括两个步骤:首先把少量的CTD作为"约束井"进行地震数据的波阻抗反演;然后利用从CTD资料获得的研究海区的温-盐关系式,结合波阻抗数据反演得到温度和盐度剖面.通过合成数据的试算表明,基于少量的CTD资料控制,利用地震数据可以反演得到高分辨率的二维温度、盐度结构.基于地震数据的温-盐结构反演方法有望弥补传统物理海洋学观测方法的不足,为海洋学研究提供大量的基础数据,有广泛的应用前景.     

地震随机噪声压缩感知迭代压制方法

复杂地表和复杂介质条件下,随机噪声往往严重影响着复杂地震信号的信噪比,同时深层地球物理目标探查中弱地震信号总是被随机噪声所掩盖,如何有效地压制随机噪声干扰、恢复有效地震信号仍然是高精度地震勘探中的关键问题.压缩感知理论突破了奈奎斯特采样定理的限制,利用有效地震信号的可压缩性和稀疏性,提供了从不可压缩随机噪声中进行有效信号分离的数据原理.本文系统分析压缩感知框架下地震随机噪声压制的稀疏优化反问题,提出了基于迭代软阈值算法的"采集-重建-修复"方案对该问题进行求解.在实现高度稀疏表征的基础上进行地震数据的压缩感知随机观测,通过迭代反演对有效地震信号进行重构,有效提高复杂地震数据的信噪比,同时,当求解稀疏优化问题时,如果出现正则化项引起重构信号衰减现象,可以匹配除偏对衰减的有效信号进行修复.通过与工业标准 f-x预测滤波方法进行比较,理论模型和实际数据处理的结果表明,压缩感知迭代噪声压制方法对复杂地震数据中的随机噪声有较好的压制效果,可以有效恢复出被较强非平稳随机噪声干扰的时空变同相轴信息.     

子波相位不准对反演结果的影响(英文)

本文重点讨论在振幅谱估计准确的情况下,采用不同相位谱子波作为实际估计子波进行线性最小二乘反演,并对结果进行分析。除子波相位外,所有其它影响反演结果的因素均忽略。稀疏反射系数模型(块状波阻抗模型)反演结果表明:(1)使用不同相位谱子波进行反演,其反演结果合成的记录与原始记录都非常匹配,但反演的反射系数和声波阻抗结果与真实模型有差异;(2)反演结果的可靠程度主要与不同相位子波z变换的根的分布有关,当估计子波与真实子波Z变换的根的分布仅在单位圆附近有差异时,反演的反射系数和声波阻抗与真实模型很接近;(3)尽管反演前后地震记录都匹配了,并且评价反演结果好坏的柯西准则或改进柯西准则(反演参数没有进行自适应处理)已经达到了最优(最小),但反演结果与真实模型仍存在较大差异。最后,针对子波相位估计不准可能导致反演效果较差这个问题,我们提出采用求L1范数、丰度、变分、柯西准则(反演参数进行了自适应处理)或/和改进柯西准则(反演参数进行了自适应处理)的最优值或次优值作为评价准则的一种解决办法,理论上得到了好的效果。     

用平面波延拓方程进行地震数据的叠前速度反演

本文讨论地震勘探数据的叠前速度反演方法及其在海洋地震勘探数据上的反演试验.反演主要的计算步骤是:1.采用Fourier-Hankel变换把球面波分解为平面谐波;2.用平面谐波的延拓方程将上行波与下行波同时向下延拓,并计算每一层底部的反射系数和下一层的波阻抗;3.用最小二乘法从波阻抗中确定该层的声波速度.重复第2步与第3步,直到某一预定深度时结束.通过反演试验,对地震振幅比例的改变,子波变形,以及第1层速度和密度的误差对反演方法的稳定性及其精度的影响进行了分析.还通过实际海洋地震勘探数据的反演试验,对这一方法在地震勘探中的应用前景作了论述.     

地震波反射系数偏导矩阵的精确解及对近似计算的精度分析

目前叠前反演方法大多是基于Zoeppritz方程近似式实现的,它仅适应于弱反射介质界面、中小角度(或小偏移距)的地震数据反演,不能满足勘探开发的地质需求.本文建立了基于zoeppritz方程精确求解反射系数的梯度矩阵,分析了矩阵特点和精度,为研究利用反射系数梯度精确解反演地震参数奠定了基础.     

用对称映射ARMA模型的零极点研究子波相位对反射系数序列反演的影响

为研究地震子波相位对反射系数序列反演的影响,在自回归滑动平均(ARMA)模型描述子波的基础上,提出采用z域对称映射ARMA模型零极点的方法构造了一系列相同振幅谱、不同相位谱的地震子波,并结合谱除法对人工合成地震记录进行反射系数序列反演.理论分析表明,子波相位估计不准时反射系数序列反演结果中残留一个纯相位滤波器,该纯相位滤波器的相位谱为真实子波和构造子波的相位谱之差.采用丰度和变分作为评价方法,在反演结果中确定出真实的或准确的反射系数序列.仿真实验和实际数据处理结果也验证了子波相位对反射系数序列反演的影响规律和评价方法的有效性,为进一步提高反射系数序列反演结果精度指明了研究方向.     

预条件共轭梯度法在地震数据重建方法中的应用

基于最小平方的Fourier地震数据重建方法最终转化为求解一个线性方程组, 其系数矩阵是Toeplitz矩阵,可以用共轭梯度法求解该线性方程组.共轭梯度法的迭代次数受系数矩阵病态程度的影响,地震数据的非规则采样程度越高,所形成的系数矩阵病态程度越高,就越难收敛和得到合理的计算结果.本文研究了基于Toeplitz矩阵的不同预条件的构造方法,以及对共轭梯度法收敛性的影响.通过预条件的使用,加快了共轭梯度法的迭代速度, 改进了共轭梯度算法的收敛性,提高了计算的效率.数值算例和实际地震数据重建试验证明了预条件共轭梯度法对计算效率有很大的提高.     

Accurate estimation of acoustic impedance based on spectral inversion

Acoustic impedance is one of the best attributes for seismic interpretation and reservoir characterisation. We present an approach for estimating acoustic impedance accurately from a band‐limited and noisy seismic data. The approach is composed of two stages: inverting for reflectivity from seismic data and then estimating impedance from the reflectivity inverted in the first stage. For the first stage, we achieve a two‐step spectral inversion that locates the positions of reflection coefficients in the first step and determines the amplitudes of the reflection coefficients in the second step under the constraints of the positions located in the first step. For the second stage, we construct an iterative impedance estimation algorithm based on reflectivity. In each iteration, the iterative impedance estimation algorithm estimates the absolute acoustic impedance based on an initial acoustic impedance model that is given by summing the high‐frequency component of acoustic impedance estimated at the last iteration and a low‐frequency component determined in advance using other data. The known low‐frequency component is used to restrict the acoustic impedance variation tendency in each iteration. Examples using one‐ and two‐dimensional synthetic and field seismic data show that the approach is flexible and superior to the conventional spectral inversion and recursive inversion methods for generating more accurate acoustic impedance models.     

An application of the autoregressive extrapolation technique to enhance deconvolution results: a 2D marine data example

Interpreting a post‐stack seismic section is difficult due to the band‐limited nature of the seismic data even post deconvolution. Deconvolution is a process that is universally applied to extend the bandwidth of seismic data. However, deconvolution falls short of this task as low and high frequencies of the deconvolved data are either still missing or contaminated by noise. In this paper we use the autoregressive extrapolation technique to recover these missing frequencies, using the high signal‐to‐noise ratio (S/N) portions of the spectrum of deconvolved data. I introduce here an algorithm to extend the bandwidth of deconvolved data. This is achieved via an autoregressive extrapolation technique, which has been widely used to replace missing or corrupted samples of data in signal processing. This method is performed in the spectral domain. The spectral band to be extrapolated using autoregressive prediction filters is first selected from the part of the spectrum that has a high signal‐to‐noise ratio (S/N) and is then extended. As there can be more than one zone of good S/N in the spectrum, the results of prediction filter design and extrapolation from three different bands are averaged. When the spectrum of deconvolved data is extended in this way, the results show higher vertical resolution to a degree that the final seismic data closely resemble what is considered to be a reflectivity sequence of the layered medium. This helps to obtain acoustic impedance with inversion by stable integration. The results show that autoregressive spectral extrapolation highly increases vertical resolution and improves horizon tracking to determine continuities and faults. This increase in coherence ultimately yields a more interpretable seismic section.     

Frequency‐domain waveform modelling and inversion for coupled media using a symmetric impedance matrix

To simulate the seismic signals that are obtained in a marine environment, a coupled system of both acoustic and elastic wave equations is solved. The acoustic wave equation for the fluid region simulates the pressure field while minimizing the number of degrees of freedom of the impedance matrix, and the elastic wave equation for the solid region simulates several elastic events, such as shear waves and surface waves. Moreover, by combining this coupled approach with the waveform inversion technique, the elastic properties of the earth can be inverted using the pressure data obtained from the acoustic region. However, in contrast to the pure acoustic and elastic cases, the complex impedance matrix for the coupled media does not have a symmetric form because of the boundary (continuity) condition at the interface between the acoustic and elastic elements. In this study, we propose a manipulation scheme that makes the complex impedance matrix for acoustic–elastic coupled media to take a symmetric form. Using the proposed symmetric matrix, forward and backward wavefields are identical to those generated by the conventional approach; thus, we do not lose any accuracy in the waveform inversion results. However, to solve the modified symmetric matrix, LDLT factorization is used instead of LU factorization for a matrix of the same size; this method can mitigate issues related to severe memory insufficiency and long computation times, particularly for large‐scale problems.     

河流相砂泥岩薄互层地震反射特征研究

A sedimentary geological model is established in order to study the seismic reflection characteristics of channel sand bodies. Synthetic seismic shot gathers are simulated using the acoustic wave equation and then are prestack time migrated. On the imaged data, the reflection characteristics and instantaneous attributes are analyzed and log-constrained impedance inversion is tested. Because of wave field interference, the experimental results show that seismic events do not definitely correspond to the channel sand bodies and that seismic modes of occurrence do not represent the actual ones. The seismic events formed by wave interference may lead to errors and pitfalls in sand body interpretation. The corresponding relations between instantaneous seismic attributes and sedimentary sands are not well established. Log-constrained impedance inversion improves the resolution of channel sands. However, if the inverted resolution is forced to be too high, artifacts related to the initial model may occur.     

宽频数据在致密砂岩储层预测中的应用

西湖凹陷深层致密砂岩储层具有良好的勘探开发前景,受埋深影响,目的层地震资料品质较差.A构造通过斜缆宽频采集和处理获取宽频地震数据,提升了资料品质,然而应用常规子波提取方法对宽频数据进行子波提取并反演计算纵波阻抗,结果与井上实测数值差异较大,影响储层的定量解释.针对这一问题,提出统计性子波和确定性子波相结合的长短子波合并宽频子波提取方法,提取的宽频子波比常规子波低频丰富、旁瓣小,能更真实地反映地震信息,约束稀疏脉冲反演的纵波阻抗结果与测井曲线吻合度更高.基于宽频数据和常规数据分别进行约束稀疏脉冲弹性波阻抗反演,预测A构造优质储层分布,经已钻井证实,宽频数据比常规数据储层预测精度高,预测的储层展布特征与研究区地质沉积认识一致.结果表明:这种基于宽频子波提取的宽频资料应用方法有效降低了致密砂岩储层预测的多解性.     

基于带状混合矩阵ICA实现地震盲反褶积

基于对地震反褶积本质上是一个盲过程的认识,引入高阶统计学盲源分离技术——独立分量分析（ICA）实现地震盲反褶积.在无噪声假设条件下,利用地震记录时间延迟矩阵和地震子波带状褶积矩阵,将地震褶积模型转化为一般线性混合ICA模型,采用FastICA算法,将带状性质作为先验信息,实现所谓带状ICA算法（B\|ICA）,得到个数与子波算子长度相等的多个估计反射系数序列和估计子波序列,最后利用褶积模型提供的附加信息从中优选出最佳的反射系数序列及相应的地震子波.模型数据和实际二维地震道数值算例表明：对于统计性反褶积,在不对反射系数作高斯白噪假设,不对子波作最小相位假设的所谓“全盲”条件下,基于ICA方法（反射系数非高斯分布,地震子波非最小相位）可以较好解决地震盲反褶积问题,是基于二阶统计特性的地震信号统计性反褶积方法的提升,具有可行性和应用前景.     

Exact elastic impedance tensor for isotropic media

The existing expressions of elastic impedance,as the generalized form of acoustic impedance,represent the resistance of subsurface media to seismic waves of non-normal incidence,and thus include information on the shear-wave velocity.In this sense,conventional elastic impedance is an attribute of the seismic reflection and not an intrinsic physical property of the subsurface media.The derivation of these expressions shares the approximations made for reflectivity,such as weak impedance contrast andisotropic or weakly anisotropic media,which limits the accuracy of reflectivity reconstruction and seismic inversion.In this paper,we derive exact elastic impedance tensors of seismic P-and S-waves for isotropic media based on the stress-velocity law.Each componentof the impedance tensor represents a unique mechanical property of the medium.Approximations of P-wave elastic impedance tensor components are discussed for seismic inversion and interpretation.Application to synthetic data and real data shows the accuracy and robust interpretation capability of the derived elastic impedance in lithology characterizations.     

Estimation of a non-stationary prior covariance from seismic data

Non-stationarity in statistical properties of the subsurface is often ignored. In a classical linear Bayesian inversion setting of seismic data, the prior distribution of physical parameters is often assumed to be stationary. Here we propose a new method of handling non-stationarity in the variance of physical parameters in seismic data. We propose to infer the model variance prior to inversion using maximum likelihood estimators in a sliding window approach. A traditional, and a localized shrinkage estimator is defined for inferring the prior model variance. The estimators are assessed in a synthetic base case with heterogeneous variance of the acoustic impedance in a zero-offset seismic cross section. Subsequently, this data is inverted for acoustic impedance using a non-stationary model set up with the inferred variances. Results indicate that prediction as well as posterior resolution is greatly improved using the non-stationary model compared with a common prior model with stationary variance. The localized shrinkage predictor is shown to be slightly more robust than the traditional estimator in terms of amplitude differences in the variance of acoustic impedance and size of local neighbourhood. Finally, we apply the methodology to a real data set from the North Sea basin. Inversion results show a more realistic posterior model than using a conventional approach with stationary variance.