共查询到20条相似文献，搜索用时 31 毫秒
1.
Geostatistical seismic inversion for non‐stationary patterns using direct sequential simulation and co‐simulation 
下载免费PDF全文
下载免费PDF全文 Hamid Sabeti Ali Moradzadeh Faramarz Doulati Ardejani Leonardo Azevedo Amilcar Soares Pedro Pereira Ruben Nunes 《Geophysical Prospecting》2017,65(Z1):25-48
Geostatistical seismic inversion methods are routinely used in reservoir characterisation studies because of their potential to infer the spatial distribution of the petro‐elastic properties of interest (e.g., density, elastic, and acoustic impedance) along with the associated spatial uncertainty. Within the geostatistical seismic inversion framework, the retrieved inverse elastic models are conditioned by a global probability distribution function and a global spatial continuity model as estimated from the available well‐log data for the entire inversion grid. However, the spatial distribution of the real subsurface elastic properties is complex, heterogeneous, and, in many cases, non‐stationary since they directly depend on the subsurface geology, i.e., the spatial distribution of the facies of interest. In these complex geological settings, the application of a single distribution function and a spatial continuity model is not enough to properly model the natural variability of the elastic properties of interest. In this study, we propose a three‐dimensional geostatistical inversion technique that is able to incorporate the reservoir's heterogeneities. This method uses a traditional geostatistical seismic inversion conditioned by local multi‐distribution functions and spatial continuity models under non‐stationary conditions. The procedure of the proposed methodology is based on a zonation criterion along the vertical direction of the reservoir grid. Each zone can be defined by conventional seismic interpretation, with the identification of the main seismic units and significant variations of seismic amplitudes. The proposed method was applied to a highly non‐stationary synthetic seismic dataset with different levels of noise. The results of this work clearly show the advantages of the proposed method against conventional geostatistical seismic inversion procedures. It is important to highlight the impact of this technique in terms of higher convergence between real and inverted reflection seismic data and the more realistic approximation towards the real subsurface geology comparing with traditional techniques. 相似文献
2.
针对利用地震道进行相对波阻抗反演中遇到的横向连续性难以保持、初始子波容错度差以及随机噪声干扰影响反演结果等问题,提出了一种基于矩阵Toeplitz稀疏分解的相对波阻抗反演方法.该方法将地震数据剖面的Toeplitz稀疏分解问题分解为两个子反演问题,其一以Toeplitz子波矩阵元素为待反演的参数,用Fused Lasso方法求解,可保证子波具有紧支集且是光滑的;其二以稀疏反射系数矩阵元素为待反演参数,用基于回溯的快速萎缩阈值迭代算法求解,大大降低了目标函数中参数选择的难度.通过交替迭代求解上述两个子反演问题可将地震数据剖面因式分解为一个Toeplitz子波矩阵和一个稀疏反射系数矩阵;然后由反射系数矩阵递推反演可以得到高分辨率的相对波阻抗剖面;利用测井资料加入低频分量后,也可得到高分辨率的绝对波阻抗剖面.Marmousi2模型生成的合成记录算例和实际地震资料算例均表明:本文方法可以从带限地震数据中有效地反演相对波阻抗,反演结果分辨率高并且能够很好地保持地震数据的横向连续性;即使在初始估计子波存在误差和地震数据被随机噪声污染的情况下也能取得较好的效果. 相似文献
3.
To reduce drilling uncertainties, zero-offset vertical seismic profiles can be inverted to quantify acoustic properties ahead of the bit. In this work, we propose an approach to invert vertical seismic profile corridor stacks in Bayesian framework for look-ahead prediction. The implemented approach helps to successfully predict density and compressional wave velocity using prior knowledge from drilled interval. Hence, this information can be used to monitor reservoir depth as well as quantifying high-pressure zones, which enables taking the correct decision during drilling. The inversion algorithm uses Gauss–Newton as an optimization tool, which requires the calculation of the sensitivity matrix of trace samples with respect to model parameters. Gauss–Newton has quadratic rate of convergence, which can speed up the inversion process. Moreover, geo-statistical analysis has been used to efficiently utilize prior information supplied to the inversion process. The algorithm has been tested on synthetic and field cases. For the field case, a zero-offset vertical seismic profile data taken from an offshore well were used as input to the inversion algorithm. Well logs acquired after drilling the prediction section was used to validate the inversion results. The results from the synthetic case applications were encouraging to accurately predict compressional wave velocity and density from just a constant prior model. The field case application shows the strength of our proposed approach in inverting vertical seismic profile data to obtain density and compressional wave velocity ahead of a bit with reasonable accuracy. Unlike the commonly used vertical seismic profile inversion approach for acoustic impedance using simple error to represent the prior covariance matrix, this work shows the importance of inverting for both density and compressional wave velocity using geo-statistical knowledge of density and compressional wave velocity from the drilled section to quantify the prior covariance matrix required during Bayesian inversion. 相似文献
4.
波阻抗反演中的不确定性分析 总被引：6，自引：2，他引：4
在波阻抗反演中引入一种新的计算机统计分析方法,通过实际地震资料作Ｂｏｏｔａｔｚａｐ试验,用重采样的方法获取大量的地震数据集,对数据集中所有地震剖面的反演结果作统计分析,获取反演结果的均值,方差,量信区间等统计结果,分析波阻抗反演的不确定性。 相似文献
5.
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. 相似文献
6.
The seismic reflection method provides high-resolution data that are especially useful for discovering mineral deposits under deep cover. A hindrance to the wider adoption of the seismic reflection method in mineral exploration is that the data are often interpreted differently and independently of other geophysical data unless common earth models are used to link the methods during geological interpretation. Model-based inversion of post-stack seismic data allows rock units with common petrophysical properties to be identified and permits increased bandwidth to enhance the spatial resolution of the acoustic-impedance model. However, as seismic reflection data are naturally bandlimited, any inversion scheme depends upon an initial model, and must deal with non-unique solutions for the inversion. Both issues can be largely overcome by using constraints and integrating prior information. We exploit the abilities of fuzzy c-means clustering to constrain and to include prior information in the inversion. The use of a clustering constraint for petrophysical values pushes the inversion process to select models that are primarily composed of several discrete rock units and the fuzzy c-means algorithm allows some properties to overlap by varying degrees. Imposing the fuzzy clustering techniques in the inversion process allows solutions that are similar to the natural geologic patterns that often have a few rock units represented by distinct combinations of petrophysical characteristics. Our tests on synthetic models, with clear and distinct boundaries, show that our methodology effectively recovers the true model. Accurate model recovery can be obtained even when the data are highly contaminated by random noise, where the initial model is homogeneous, or there is minimal prior petrophysical information available. We demonstrate the abilities of fuzzy c-means clustering to constrain and to include prior information in the acoustic-impedance inversion of a challenging magnetotelluric/seismic data set from the Carlin Gold District, USA. Using fuzzy c-means guided inversion of magnetotelluric data to create a starting model for acoustic-impedance proved important in obtaining the best result. Our inversion results correlate with borehole data and provided a better basis for geological interpretation than the seismic reflection images alone. Low values of the acoustic impedance in the basement rocks were shown to be prospective by geochemical analysis of rock cores, as would be predicted for later gold mineralization in weak, decalcified rocks. 相似文献
7.
8.
Guido Baeten Jan Willem de Maag René‐Edouard Plessix Rini Klaassen Tahira Qureshi Maren Kleemeyer Fons ten Kroode Zhang Rujie 《Geophysical Prospecting》2013,61(4):701-711
Velocity model building and impedance inversion generally suffer from a lack of intermediate wavenumber content in seismic data. Intermediate wavenumbers may be retrieved directly from seismic data sets if enough low frequencies are recorded. Over the past years, improvements in acquisition have allowed us to obtain seismic data with a broader frequency spectrum. To illustrate the benefits of broadband acquisition, notably the recording of low frequencies, we discuss the inversion of land seismic data acquired in Inner Mongolia, China. This data set contains frequencies from 1.5–80 Hz. We show that the velocity estimate based on an acoustic full‐waveform inversion approach is superior to one obtained from reflection traveltime inversion because after full‐waveform inversion the background velocity conforms to geology. We also illustrate the added value of low frequencies in an impedance estimate. 相似文献
9.
《中国科学:地球科学(英文版)》2015,(8)
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. 相似文献
10.
子波相位不准对反演结果的影响(英文) 总被引：5，自引：1，他引：4
本文重点讨论在振幅谱估计准确的情况下,采用不同相位谱子波作为实际估计子波进行线性最小二乘反演,并对结果进行分析。除子波相位外,所有其它影响反演结果的因素均忽略。稀疏反射系数模型(块状波阻抗模型)反演结果表明:(1)使用不同相位谱子波进行反演,其反演结果合成的记录与原始记录都非常匹配,但反演的反射系数和声波阻抗结果与真实模型有差异;(2)反演结果的可靠程度主要与不同相位子波z变换的根的分布有关,当估计子波与真实子波Z变换的根的分布仅在单位圆附近有差异时,反演的反射系数和声波阻抗与真实模型很接近;(3)尽管反演前后地震记录都匹配了,并且评价反演结果好坏的柯西准则或改进柯西准则(反演参数没有进行自适应处理)已经达到了最优(最小),但反演结果与真实模型仍存在较大差异。最后,针对子波相位估计不准可能导致反演效果较差这个问题,我们提出采用求L1范数、丰度、变分、柯西准则(反演参数进行了自适应处理)或/和改进柯西准则(反演参数进行了自适应处理)的最优值或次优值作为评价准则的一种解决办法,理论上得到了好的效果。 相似文献
11.
河流相砂泥岩薄互层地震反射特征研究 总被引：10，自引：3，他引：7
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. 相似文献
12.
澳大利亚西北大陆架卡拉汶盆地含有丰富的油气资源, 但其地质条件复杂, 普遍存在异常地层高压分布,特别是该区区域盖层-中深层巨厚的Muderong页岩层内存在局部异常高压带,给石油钻井带来巨大困难.本文根据卡拉汶盆地北部地区10余口井的实测地层压力资料、声波测井数据和联井地震剖面, 研究该页岩层内异常压力千赫兹尺度声波响应的衰减特征,及其与地层有效压力之间的关联计算模型.通过联井地震数据的井-震相关性分析进行井-震过渡,研究井旁地震数据在异常高压带的衰减特征,及其与地层有效压力之间的关联计算模型,从而实现井震资料联合进行地层压力地震外推反演.研究表明,Muderong页岩层地层有效压力与井中声波速度/阻抗品质因子和井旁地震阻抗品质因子具有良好的正向相关性,根据二者交汇数据拟合经验公式反算得到的各井有效压力结果基本反映出实测有效压力曲线的变化特征.卡拉汶盆地北部地区地震数据具有良好的井-震相关性,地震资料品质相对可靠,确保了声-压相关分析的井-震过渡.最后,通过地震阻抗品质因子与地层有效压力之间的关联计算模型,进行该区异常地层压力地震外推反演,预测结果与钻井证实的异常压力分布特征基本符合.本研究为卡拉汶盆地中深部区域盖层局部异常高压带的地层压力外推反演提供了一种有效途径. 相似文献
13.
The generalized Radon transform (GRT) inversion contains an explicit relationship between seismic amplitude variations, the reflection angle and the physical parameters which can be used to describe the earth efficiently for inversion purposes. Using this relationship, we have derived parametrizations for acoustic and P–P scattering so that the variations in seismic amplitude with reflection angle for each parameter are sufficiently independent. These parametrizations show that small offset and large offset amplitudes are related to different physical parameters. In the case of acoustic scattering, the small-offset amplitudes are related to impedance variations while large-offset amplitudes are related to velocity variations. A similar result has been established for P–P scattering. The Born approximation (which is used to derive the GRT inversion) does not correctly predict the amplitude due to velocity variations at large offsets, and thus the inversion of velocity is not as satisfactory as the inversion of impedance. 相似文献
14.
总结了牛顿类地震反演方法中Hessian算子的作用,对其在地震反演成像中的数学物理含义进行了分析.Hessian算子是误差泛函对模型参数的二阶导数,反映了误差泛函对模型变化的二次型特征.分析声波方程下的Hessian算子的格林函数表达形式,发现其表达了整个观测系统和子波频带等因素对地震数据空间到模型空间投影过程的影响.提出了两种分别适用于最小二乘偏移和全波形反演的Hessian算子简化格式.平面波Hessian算子应用于最小二乘偏移能够得到相对保真的成像结果,改善了地震偏移成像的精度.地下偏移距Hessian算子应用于全波形反演能够加快反演迭代的计算效率.最后,对Hessian算子在地震反演成像中的价值进行了讨论和评价. 相似文献
15.
随机地震反演关键参数优选和效果分析(英文) 总被引：2，自引：0，他引：2
随机地震反演技术是将地质统计理论和地震反演相结合的反演方法,它将地震资料、测井资料和地质统计学信息融合为地下模型的后验概率分布,利用马尔科夫链蒙特卡洛(MCMC)方法对该后验概率分布采样,通过综合分析多个采样结果来研究后验概率分布的性质,进而认识地下情况。本文首先介绍了随机地震反演的原理,然后对影响随机地震反演效果的四个关键参数,即地震资料信噪比、变差函数、后验概率分布的样本个数和井网密度进行分析并给出其优化原则。资料分析表明地震资料信噪比控制地震资料和地质统计规律对反演结果的约束程度,变差函数影响反演结果的平滑程度,后验概率分布的样本个数决定样本统计特征的可靠性,而参与反演的井网密度则影响反演的不确定性。最后通过对比试验工区随机地震反演和基于模型的确定性地震反演结果,指出随机地震反演可以给出更符合地下实际情况的模型。 相似文献
16.
Joint AVO inversion, wavelet estimation and noise-level estimation using a spatially coupled hierarchical Bayesian model 总被引：1，自引：0，他引：1
The main objective of the AVO inversion is to obtain posterior distributions for P-wave velocity, S-wave velocity and density from specified prior distributions, seismic data and well-log data. The inversion problem also involves estimation of a seismic wavelet and the seismic-noise level. The noise model is represented by a zero mean Gaussian distribution specified by a covariance matrix. A method for joint AVO inversion, wavelet estimation and estimation of the noise level is developed in a Bayesian framework. The stochastic model includes uncertainty of both the elastic parameters, the wavelet, and the seismic and well-log data. The posterior distribution is explored by Markov-chain Monte-Carlo simulation using the Gibbs' sampler algorithm. The inversion algorithm has been tested on a seismic line from the Heidrun Field with two wells located on the line. The use of a coloured seismic-noise model resulted in about 10% lower uncertainties for the P-wave velocity, S-wave velocity and density compared with a white-noise model. The uncertainty of the estimated wavelet is low. In the Heidrun example, the effect of including uncertainty of the wavelet and the noise level was marginal with respect to the AVO inversion results. 相似文献
17.
Yuanqiang Li Jingye Li Xiaohong Chen Jian Zhang Xin Bo 《Geophysical Prospecting》2020,68(7):2009-2026
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. 相似文献
18.
Traditional coherence algorithms are often based on the assumption that seismic traces are stationary and Gaussian. However, seismic traces are actually non-stationary and non-Gaussian. A constant time window and the canonical correlation analysis in traditional coherence algorithms are not optimal for non-stationary seismic traces and cannot describe the similarity between adjacent seismic traces in detail. To overcome this problem, a new coherence algorithm using the high-resolution time–time transform and the feature matrix is designed. The high-resolution time–time transform used to replace the constant time window can produce a frequency-dependent time local series to analyse non-stationary seismic traces. The feature matrix, constructed by the frequency-dependent time local series and the related local gradients, defines a new correlation metric that enhances more details of the geological discontinuities in seismic images than does the canonical correlation analysis. Additionally, the Riemannian metric is introduced for related calculations because the feature matrices are not defined in a Euclidean space but rather in a manifold space. Application to field data illustrates that the proposed method reveals more details of structural and stratigraphic features. 相似文献
19.
Exploring seismic inversion methodologies for non‐stationary geological environments: a benchmark study between deterministic and geostatistical seismic inversion 下载免费PDF全文
下载免费PDF全文 This paper presents a comparison between subsurface impedance models derived from different deterministic and geostatistical seismic inversion methodologies applied to a challenging synthetic dataset. Geostatistical seismic inversion methodologies nowadays are common place in both industry and academia, contrasting with traditional deterministic seismic inversion methodologies that are becoming less used as part of the geo‐modelling workflow. While the first set of techniques allows the simultaneous inference of the best‐fit inverse model along with the spatial uncertainty of the subsurface elastic property of interest, the second family of inverse methodology has proven results in correctly predicting the subsurface elastic properties of interest with comparatively less computational cost. We present herein the results of a benchmark study performed over a realistic three‐dimensional non‐stationary synthetic dataset in order to assess the performance and convergence of different deterministic and geostatistical seismic inverse methodologies. We also compare and discuss the impact of the inversion parameterisation over the exploration of the model parameter space. The results show that the chosen seismic inversion methodology should always be dependent on the type and quantity of the available data, both seismic and well‐log, and the complexity of the geological environment versus the assumptions behind each inversion technique. The assessment of the model parameter space shows that the initial guess of traditional deterministic seismic inversion methodologies is of high importance since it will determine the location of the best‐fit inverse solution. 相似文献
20.
在利用地震波数据进行地球物理反演时,地层对地震波的吸收衰减效应会对地层物性参数的准确反演产生较大的影响,因此利用黏弹性声波方程进行反演更符合实际情形.本文在考虑地层衰减效应进行频率空间域正演模拟的基础上,提出基于黏弹性声波方程的频率域逆散射反演算法并对地震波传播速度进行反演重建,在反演过程中分别用地震波传播复速度和实速度来表征是否考虑地层吸收衰减效应.基于反演参数总变差的正则化处理使反演更加稳定,在反演中将低频反演速度模型作为高频反演的背景模型进行逐频反演,由于单频反演过程中背景模型保持不变,故该方法不需要在每次迭代中重新构造正演算子,具有较高的反演效率;此外本文在反演过程中采用了基于MPI的并行计算策略,进一步提高了反演计算的效率.在二维算例中分别对是否考虑地层吸收衰减效应进行了地震波速度反演,反演结果表明考虑衰减效应可以得到与真实模型更加接近的速度分布结果,相反则无法得到正确的地震波速度重建结果.本文算法对复杂地质模型中浅层可以反演得到分辨率较高的速度模型,为其他地震数据处理提供比较准确的速度信息,在地层深部由于地震波能量衰减导致反演分辨率不太理想. 相似文献