基于Metropolis抽样的非线性反演方法
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
基于Metropolis抽样的非线性反演应用贝叶斯理论框架,是一种基于蒙特卡洛的非线性反演方法,能够有效地融合测井资料中的高频信息,提高反演结果的分辨率。首先通过快速傅里叶滑动平均模拟算法(FFTMA)和逐渐变形算法(GDM)得到基于地质统计学的先验信息;进而构建似然函数;最后利用Metropolis算法对后验概率密度进行抽样,得到反演问题的解。其中FFT-MA模拟作为一种高效的频率域模拟方法,融入GDM更新算法之后,可以在保持模拟空间结构不变的前提下,连续修改储层模型,保证反演结果有效地收敛,直至满足实际观测地震记录。模型试算和实际数据处理结果表明:基于Metropolis抽样的非线性反演可以提供合理的弹性参数信息,尤其是提高纵波速度的分辨率,即使信噪比较小时,仍然可以反演出合理的弹性参数信息,从而证明了该方法的有效性;当不考虑噪声时,纵、横波阻抗的反演分辨率较弹性参数本身的反演分辨率更高。
Nonlinear inversion based on Metropolis sampling algorithm is formulated in the Bayesian framework.As one kind of Monte Carl non-linear inversions,it can effectively integrate high frequency information of well logging data,and obtain inversion results with a higher resolution.Firstly,we get the priori information through fast Fourier transform moving average(FFT-MA)and gradual deformation method(GDM).Second,we structure likelihood function.Then we apply Metropolis algorithm in order to obtain an exhaustive characterization of the posteriori probability density.FFT-MA is a kind of efficient simulation method.Combined with GDM,it can constantly modify reservoir model and keep the spatial structure unchanged until it matches the observed seismic data.According to the model trial and real data processing,we can conclude that nonlinear inversion based on Metropolis sampling algorithm provide reasonable elastic parameter information,especially it improves the resolution of P-wave velocity.Even when the signal noise ratio(SNR)is relatively low,it can still show reasonable elastic parameter information,which proves the effectiveness of the proposed method.The inversion resolution of Pwave and S-wave impedances is higher than elastic parameters inversion if we do not consider the noise.
引文
[1]Sancevero S S,Remacre A Z,de Souza Portugal R et al.Comparing deterministic and stochastic seismic inversion for thin-bed reservoir characterization in a turbidite synthetic reference model of Campos Basin,Brazil.The Leading Edge,2005,24(11):1168-1172.
    [2]Moyen R,Doyen P M.Reservoir connectivity uncertainty from stochastic seismic inversion.SEG Technical Program Expanded Abstracts,2009,28:2378-2382.
    [3]Sams M S,Saussus D.Comparison of uncertainty estimates from deterministic and geostatistical inversion.The 70th EAGE Conference&Exhibition,2008.
    [4]Francis A.Limitations of deterministic and advantages of stochastic seismic inversion.Canadian Society of Exploration Geophysicists Recorder,2005,30(2):5-11.
    [5]Francis A.Understanding stochastic inversion:Part 1.First Break,2006,24(11):69-77.
    [6]Francis A.Understanding stochastic inversion:Part 2.First Break,2006,24(12):79-84.
    [7]Dubrule O.Workshop report:“Uncertainty in reserve estimates”EAGE Conference,Amsterdam,2 June1996.Petroleum Geoscience,1996,2(4):351-352.
    [8]Kjnsberg H,Hauge R,Kolbjrnsen O et al.Bayesian Monte Carlo method for seismic predrill prospect assessment.Geophysics,2010,75(2):O9-O19.
    [9]张繁昌,肖张波,印兴耀.地震数据约束下的贝叶斯随机反演.石油地球物理勘探,2014,49(1):176-182.Zhang Fanchang,Xiao Zhangbo,Yin Xingyao.Bayesian stochastic inversion constrained by seismic data.OGP,2014,49(1):176-182.
    [10]Le Ravalec M,Noetinger B,Hu L Y.The FFT moving average(FFT-MA)generator:An efficient numerical method for generating and conditioning Gaussian simulations.Mathematical Geology,2000,32(6):701-723.
    [11]Hu L Y.Gradual deformation and iterative calibration of Gaussian-related stochastic models.Mathematical Geology,2000,32(1):87-108.
    [12]John A S,Martin L S,Sven T.Introductory Geophysical Inverse Theory.Colorado and New England:Samizdat Press,2001.
    [13]Oliver D S.Moving averages for Gaussian simulation in two and three dimensions.Mathematical Geology,1995,27(8):939-960.
    [14]Duijndam A J W.Bayesian estimation in seismic inversion,Part II:Uncertainty analysis.Geophysical Prospecting,1988,36(8):899-918.
    [15]Journel A G,Huijbregts C J.Mining Geostatistics.London,New York and San Francisco:Academic Press,1978.
    [16]Hu L Y,Blanc G.Constraining a reservoir facies model to dynamic data using agradual deformation method.The 6th European Conference on the Mathematics of Oil Recovery,1998.
    [17]Hu L Y,Le Ravalec M,Blanc G et al.Reducing uncertainties in production forecasts by constraining geological modeling to dynamic data.SPE Annual Technical Conference and Exhibition,Society of Petroleum Engineers,1999.
    [18]Mosegaard K,Tarantola A.Monte Carlo sampling of solutions to inverse problems.Journal of Geophysical Research:Solid Earth(1978-2012),1995,100(B7):12431-12447.
    [19]Metropolis N,Rosenbluth A W,Rosenbluth M N et al.Equation of state calculations by fast computing machines.The Journal of Chemical Physics,1953,21(6):1087-1092.
    [20]Tarantola A.Inverse Problem Theory and Methods for Model Parameter Estimation.Paris:SIAM,2005.

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心