储层弹性与物性参数地震叠前同步反演的确定性优化方法
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
储层弹性与物性参数可直接应用于储层岩性预测和流体识别,是储层综合评价和油气藏精细描述的基本要素之一.现有的储层弹性与物性参数地震同步反演方法大都基于Gassmann方程,使用地震叠前数据,通过随机优化方法反演储层弹性与物性参数;或基于Wyllie方程,使用地震叠后数据,通过确定性优化方法反演储层弹性与物性参数.本文提出一种基于Gassmann方程、通过确定性优化方法开展储层弹性和物性参数地震叠前反演的方法,该方法利用Gassmann方程建立储层物性参数与叠前地震观测数据之间的联系,在贝叶斯反演框架下以储层弹性与物性参数的联合后验概率为目标函数,通过将目标函数的梯度用泰勒公式展开得到储层弹性与物性参数联合的方程组,其中储层弹性参数对物性参数的梯度用差分形式表示,最后通过共轭梯度算法迭代求解得到储层弹性与物性参数的最优解.理论试算与实际资料反演结果证明了方法的可行性.
Elastic and petrophysical parameters of reservoir can be directly applied to lithology prediction and fluid identification,they are the basic elements in both comprehensive evaluation and fine characterization of hydrocarbon reservoir.Existing seismic joint inversion methods for estimating elastic and petrophysical parameters of reservoir are mainly based on either Gassmann equation,with which these parameters are inversed from prestack seismic data by stochastic optimization methods;or Wyllie equation,with which these parameters are inverted from poststack seismic data by deterministic optimization methods.The purpose of this work is todevelop a strategy for estimating elastic and petrophysical parameters of reservoir based on Gassmann equation with deterministic seismic prestack inversion.We use Gassmann equation to build up the relationship between prestack seismic data and petrophysical parameters.We treat the joint posterior probability of elastic and petrophysical parameters as the objective function under the Bayesian architecture,by expanding the the objective function with Taylor formula,the joint equations composed of physical and petrophysical parameters can thus be obtained,and the derivatives of the elastic parameters with respect to petrophysical parameters are obtained by differentiating.The conjugate gradient method is then used to find the optimal solutions of P-wave velocity,shear velocity,density,porosity,water saturation and clay content.We apply this inversion method to prestack seismic data set in Sulige gas field,which is one of the largest gas fields in central Ordos Basin in China.The main gas-producing interval is the H8 sandstone reservoir at the bottom of Lower Shihezi Formation,Lower Permian,Upper Paleozoic.Drilling results indicate that this region is a composite sandstone gas reservoir and its heterogeneous reservoirs are characterized by low permeability,low pressure and low abundance,which bring great challenge to the exploitation of this field.From the comparison of the inversion result at the well location using our method and two-step method,we can see that there are slight but noticeable differences between these two methods,especially for the elastic parameters.Our joint inversion is significantly closer to the "true"elastic and reservoir model than the two-step inversion.Later we use 2-D field prestack seismic data to assess the performance of the proposed method,and selected a 160 ms time interval of interest for the reservoir description.Three wells that intercept the section are previously characterized in two-way travel time,and used to calibrate the rock physics model and generate low-frequency starting model.The major seismic events in the inversion result show more continuity along the reference structural horizon direction,and the estimated property have a good match with the corresponding well-log curves,and commensurate with the frequency content of the seismic data,as expected.We have developed a method for inverting prestack seismic data under well-log constraints derived from Gassmann fluid substitution relations and prior petrophysical models.In our model,we use Gassmann fluid substitution relations calibrated to the well-log data for more robust modeling,which is consistent with common petrophysical knowledge.In this specific setting,the extended application to invert prestack seismic data can estimate water saturation,porosity,clay volume,P-wave velocity,shear wave velocity and density jointly.Compared with the stochastic optimization method,our method utilize gradient information of objective function,the optimal solution can be obtained more quickly in favor of large-scale production mission.Compared with Bosch′s poststack deterministic optimization,prestack seismic amplitude have more information varying with the angle of incidence.And compared with Wyllie equation,Gassmann equation is better for establishing the relationship between the reservoir parameters and prestack seismic observation data,which helps a lot to understand the impact of reservoir parameters on the elastic properties.
引文
Aki K,Richards P G.1980.Quantitative Seismology:Theory andMethods.San Francisco:Company.
    Avseth P,Mukerji T,Jrstad A,et al.2001.Seismic reservoirmapping from 3-D AVO in a North Sea turbidite system.Geophysics,66(4):1157-1176.
    Bachrach R.2006.Joint estimation of porosity and saturation usingstochastic rock-physics modeling.Geophysics,71(5):O53-O63.
    Bosch M.2004.The optimization approach to lithological tomography:Combining seismic data and petrophysics for porosity prediction.Geophysics,69(5):1272-1282.
    Bosch M,Cara L,Rodrigues J,et al.2007.A Monte Carloapproach to the joint estimation of reservoir and elasticparameters from seismic amplitudes.Geophysics,72(6):O29-O39.
    Bosch M,Carvajal C,Rodrigues J,et al.2009.Petrophysicalseismic inversion conditioned to well-log data:Methods andapplication to a gas reservoir.Geophysics,74(2):O1-O15.
    Bosch M,Mukerji T,Gonzalez E F.2010.Seismic inversion forreservoir properties combining statistical rock physics andgeostatistics:A review.Geophysics,75(5):75A165-75A176.
    Buland A,Kolbjrnsen O,Hauge R,et al.2008.Bayesian lithologyand fluid prediction from seismic prestack data.Geophysics,73(3):C13-C21.
    Eidsvik J,Avseth P,Omre H,et al.2004.Stochastic reservoircharacterization using prestack seismic data.Geophysics,69(4):978-993.
    Elliott S E,Wiley B F.1975.Compressional velocities of partiallysaturated,unconsolidated sands.Geophysics,40(6):949-954.
    Gidlow P M,Smith G C,Vail P.1993.Hydrocarbon detectionusing fluid factor traces:A case study.//3rd SAGA BiennialConference and Exhibition.Expanded Abstracts.
    González E F,Mukerji T,Mavko G.2007.Seismic inversioncombining rock physics and multiple-point geostatistics.Geophysics,73(1):R11-R21.
    Grana D,Rossa E D.2010.Probabilistic petrophysical-propertiesestimation integrating statistical rock physics with seismicinversion.Geophysics,75(3):O21-O37.
    Grana D,Mukerji T,Dvorkin J,et al.2012.Stochastic inversion offacies from seismic data based on sequential simulations andprobability perturbation method.Geophysics,77(4):M53-M72.
    Gunning J,Glinsky M E.2007.Detection of reservoir quality usingBayesian seismic inversion.Geophysics,72(3):R37-R49.
    Hampson D P,Russell B H,Bankhead B.2005.Simultaneousinversion of pre-stack seismic data.//75th Annual InternationalMetting,Society of Exploration Geophysicists,ExpandedAbstracts.
    Kjnsberg H,Hauge R,Kolbjrnsen O,et al.2010.BayesianMonte Carlo method for seismic predrill prospect assessment.Geophysics,75(2):O9-O19.
    Larsen A L,Ulvmoen M,Omre H,et al.2006.Bayesian lithology/fluid prediction and simulation on the basis of a Markov-chainprior model.Geophysics,71(5):R69-R78.
    Loures L G L,Moraes F S.2006.Porosity inference andclassification of siliciclastic rocks from multiple data sets.Geophysics,71(5):O65-O76.
    Mavko G,Mukerji T,Dvorkin J.2009.The Rock PhysicsHandbook:Tools for Seismic Analysis of Porous Media.Cambridge:Cambridge University Press.
    Mukerji T,Jrstad A,Avseth P,et al.2001.Mapping lithofaciesand pore-fluid probabilities in a North Sea reservoir:Seismicinversions and statistical rock physics.Geophysics,66(4):988-1001.
    Rimstad K,Omre H.2010.Impact of rock-physics depth trends andMarkov random fields on hierarchical Bayesian lithology/fluidprediction.Geophysics,75(4):R93-R108.
    Saleh A A,Castagna J P.2004.Revisiting the Wyllie time averageequation in the case of near-spherical pores.Geophysics,69(1):45-55.
    Saltzer R,Finn C,Burtz O.2005.Predicting V shale and porosityusing cascaded seismic and rock physics inversion.The LeadingEdge,24(7):732-736.
    Sengupta M,Bachrach R.2007.Uncertainty in seismic-based payvolume estimation:Analysis using rock physics and Bayesianstatistics.The Leading Edge,26(2):184-189.
    Spikes K,Mukerji T,Dvorkin J,et al.2007.Probabilistic seismicinversion based on rock-physics models.Geophysics,72(5):R87-R97.
    Ulvmoen M,Omre H.2010.Improved resolution in Bayesianlithology/fluid inversion from prestack seismic data and wellobservations:Part 1—Methodology.Geophysics,75(2):R21-R35.
    Ulvmoen M,Omre H,Buland A.2010.Improved resolution inBayesian lithology/fluid inversion from prestack seismic data andwell observations:Part 2—Real case study.Geophysics,75(2):B73-B82.

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