共查询到17条相似文献，搜索用时 109 毫秒
1.
《地球物理学进展》2020,(2)
对弹性参数进行转换而获取物性参数的方法,适用范围有限且预测精度不高.本文基于岩石物理理论,提出一种带先验约束的碎屑岩储层物性参数预测方法,为储层预测、流体识别、储层描述、储量估计、烃类开发方案设计提供重要参数.应用岩石物理理论和物性参数先验信息建立了带先验约束的目标函数;并利用最优化方法对反演目标函数进行优化求解,对储层的物性及孔隙度和泥质含量进行同步迭代反演.经过正演数据及实际地震资料的测试应用,其结果表明本文方法反演得到的孔隙度和泥质含量为储层综合预测提供了愈加丰富的评价指标,同时孔隙度和泥质含量数据的联合应用对减小储层预测的多解性和降低勘探开发的风险有重要意义. 相似文献
2.
《中国科学:地球科学》2021,(7)
碳酸盐岩储层具有复杂的孔隙结构,其显著影响了岩石的弹性和地震响应特征,也影响物性参数预测精度.现有的岩石物理反演方法主要是针对碎屑岩储层,反演的对象主要是孔隙度和饱和度,采用的数据主要是弹性参数,使用的算法主要是线性近似法,缺少针对碳酸盐岩储层的孔隙结构与物性参数反演方法.针对这些问题,文章首次提出了一种基于弹性阻抗的适用于碳酸盐岩复杂孔隙储层的孔隙结构和物性参数贝叶斯非线性同步反演新方法,该方法把多孔介质岩石物理模型、Gassmann方程、AVO理论、贝叶斯理论和非线性反演算法结合起来,实现复杂孔隙储层的孔隙结构和物性参数的定量同步预测.模型正演揭示,孔隙结构参数即孔隙扁度对AVO反射系数和弹性阻抗的贡献度仅次于孔隙度,远远大于含水饱和度的影响.实际资料应用表明,文章提出的从叠前道集中直接反演物性参数和等效孔隙扁度新方法可以准确预测储层孔隙度和饱和度及其空间展布,并可评价有效储层的孔隙结构. 相似文献
3.
岩石物理是研究储层参数和岩石弹性参数之间关系的基础,在储层特性和地震特性之间起到了桥梁作用.为了更加精细的识别岩性,更好的进行砂泥岩储层油气的识别和预测,需要基于岩石物理寻找有效的识别方法.本文基于微分等效介质理论构建了砂泥岩储层岩石物理模型.为了考虑砂泥岩储层流体分布的非均匀性,模型中采用了斑块状饱和理论.文中分析了砂泥岩弹性参数和储层物性之间的关系,给出了各物性参数对弹性参数的影响结果,并根据储层物性的影响绘制了弹性参数交会模板.不同于常规交会分析,该弹性参数交会模板综合考虑了孔隙度、含水饱和度、泥质含量和孔隙纵横比等参数的影响,更加真实地反应了储层的地下特征.该交会模板的优势在于可以合理的界定出储层的砂、泥岩范围,而且根据泥质含量等参数的不同,可以进一步界定储层岩性的范围,为储层岩性的精细识别提供了新的方法.通过实际井资料的应用,该交会模板做到了砂泥岩储层岩性的有效区分,另外根据交会模板中泥质含量的不同,由0.5到0.25,我们可以在地震剖面中做出砂岩储层范围的精细界定,较好地识别了岩性,为寻找储层提供可靠依据. 相似文献
4.
叠前地震反演和岩石物理反演分别是获取弹性参数和物性参数的重要手段,两者结合有助于实现储层参数预测并精细刻画储层特征.储层物性参数的反演依赖于岩石物理模型,在进行物性参数反演时可以将复杂的岩石物理模型做泰勒展开,进而得到其一阶或高阶的近似表达式,然而这会降低模型的精确性并增加反演的误差.为了提高储层物性参数反演的稳定性和准确性,本文以碎屑岩储层为例,提出了岩石物理驱动的储层物性参数非线性地震反演方法.首先,基于贝叶斯框架和高斯分布约束条件,从叠前地震数据中实现纵、横波速度及密度等弹性参数的反演.其次,通过碎屑岩岩石物理模型建立起弹性参数与物性参数之间的联系.最后,利用粒子群算法进行全局寻优获得较为准确的孔隙度、泥质含量和含水饱和度等物性参数.合成数据和实际资料测试结果验证了所提方法的可行性和准确性,反演结果与测井数据吻合较好,可有效指示含气储层区域,本文方法在储层预测和评价方面具有广泛的应用前景. 相似文献
5.
数字岩心微观孔隙结构十分复杂,有限元模拟物性参数与弹性参数之间关系是非线性的,直接反演其物性参数准确度低、稳定性差.本文发展了一种数字岩石物理逆建模方法,实现了基于数字岩心的储层参数有效预测.从数字岩心基函数的构建出发,基于有限元方法,计算了一系列具有等间距物性参数值(孔隙度、泥质含量和含水饱和度)的数字岩心弹性参数(体积模量、剪切模量和密度),通过插值算法建立了数字岩心弹性参数三维数据集,从而实现了弹性模量的有限元数值解的快速构建;然后搜索弹性参数的单值等值面,通过等值面的空间交会得到交点,完成储层参数预测.测试结果表明:基于数字岩心逆建模理论的储层参数预测结果与实际模型一致,具有可行性,并且可以通过增加插值点数目提高预测的准确性;孔隙度和泥质含量预测结果稳定性很好,而含水饱和度对噪声的加入较为敏感. 相似文献
6.
复杂孔隙储层往往同时发育孔缝洞等多种孔隙类型，这种孔隙结构的复杂性使得岩石的速度与孔隙度之间的相关性很差.经典的二维岩石物理模版只研究弹性参数与孔隙度和饱和度之间的定量关系，而不考虑孔隙结构的影响，用这样的模版来预测复杂孔隙储层的物性参数时带来很大偏差.本文首先证明多重孔隙岩石的干骨架弹性参数可以用一个等效孔隙纵横比的单重孔隙岩石物理模型来模拟；进而基于等效介质岩石物理理论和Gassmann方程，建立一个全新的三维岩石物理模版，用它来建立复杂孔隙岩石的弹性性质与孔隙扁度及孔隙度和饱和度之间的定量关系；在此基础上，预测复杂储层的孔隙扁度、孔隙度以及孔隙中所包含的流体饱和度.实际测井和地震反演数据试验表明，三维岩石物理模版可有效提高复杂孔隙储层参数的预测精度. 相似文献
7.
《中国科学:地球科学》2020,(6)
地震岩石物理是连接岩石弹性参数与储层物性参数的桥梁,叠前地震反演是实现地下岩石弹性、物性、岩性及含流体性质定量表征的重要方法.文章构建了碎屑岩地震岩石物理高阶近似模型,推导了利用岩石模量高阶近似(Jacobian、Hessian矩阵)表征的叠前地震AVO反射特征方程,并分析了岩石孔隙度、泥质含量及流体饱和度对AVO反射率的贡献度,探讨了此方程在岩石物性参数直接预测方面的可行性.以此为基础,在待反演模型参数服从混合概率先验模型的前提下,文章提出了基于差分进化-马尔可夫链蒙特卡罗随机模型的相约束叠前地震概率化反演方法,兼具差分进化算法的全局寻优特性和马尔可夫链蒙特卡罗模型的不确定性分析能力;通过多条马尔可夫链的交叉并行,可以同步获得待反演模型参数的多个随机解,进而模拟待反演模型的后验概率密度分布,后验均值作为待反演模型的最优解,方差与置信区间用来评价反演结果的不确定性,实现储层弹性、物性、离散岩相及干岩石骨架等参数的同步预测.通过理论试验和实际资料处理验证了该理论方法的有效性. 相似文献
8.
储层弹性与物性参数可直接应用于储层岩性预测和流体识别,是储层综合评价和油气藏精细描述的基本要素之一.现有的储层弹性与物性参数地震同步反演方法大都基于Gassmann方程,使用地震叠前数据,通过随机优化方法反演储层弹性与物性参数;或基于Wyllie方程,使用地震叠后数据,通过确定性优化方法反演储层弹性与物性参数.本文提出一种基于Gassmann方程、通过确定性优化方法开展储层弹性和物性参数地震叠前反演的方法,该方法利用Gassmann方程建立储层物性参数与叠前地震观测数据之间的联系,在贝叶斯反演框架下以储层弹性与物性参数的联合后验概率为目标函数,通过将目标函数的梯度用泰勒公式展开得到储层弹性与物性参数联合的方程组,其中储层弹性参数对物性参数的梯度用差分形式表示,最后通过共轭梯度算法迭代求解得到储层弹性与物性参数的最优解.理论试算与实际资料反演结果证明了方法的可行性. 相似文献
9.
提出了各向异性页岩储层统计岩石物理反演方法.通过统计岩石物理模型建立储层物性参数与弹性参数的定量关系，使用测井数据及井中岩石物理反演结果作为先验信息，将地震阻抗数据定量解释为储层物性参数、各向异性参数的空间分布.反演过程在贝叶斯框架下求得储层参数的后验概率密度函数，并从中得到参数的最优估计值及其不确定性的定量描述.在此过程中综合考虑了岩石物理模型对复杂地下介质的描述偏差和地震数据中噪声对反演不确定性的影响.在求取最大后验概率过程中使用模拟退火优化粒子群算法以提高收敛速度和计算准确性.将统计岩石物理技术应用于龙马溪组页岩气储层，得到储层泥质含量、压实指数、孔隙度、裂缝密度等物性，以及各向异性参数的空间分布及相应的不确定性估计，为页岩气储层的定量描述提供依据. 相似文献
10.
传统基于单相介质理论的储层参数反演方法将孔隙流体与固体骨架等效为单一固体,弱化了孔隙流体的影响,反演结果精度不高.本文提出根据双相介质理论反演储层参数的方法.首先,在前人研究的基础上,利用岩石物理模型建立弹性参数与孔隙度、饱和度、泥质含量等储层参数间的关系,进而将双相介质反射系数推导为储层参数的函数;其次,根据贝叶斯反演理论,在高斯噪声假设的基础上,采用更加符合实际情况的修正柯西分布函数描述反射系数的稀疏性,推导出储层物性参数目标反演函数;最后,应用差分进化非线性全局寻优算法来求解目标反演函数,使得反演结果与实际资料间误差最小.新方法旨在突出流体对介质反射系数的影响,以期得到较高的储层参数反演精度.模型与实际资料测试均表明该方法可行、有效且反演精度较高. 相似文献
11.
含气饱和度预测是天然气储层地震解释工作的重要目标.本文将岩石物理分析与地震物理模拟技术相结合,构建了部分;饱和砂岩储层物理模型并进行含气饱和度预测分析.物理模型中设置了高孔渗常规砂岩和低孑孔渗致密砂岩两种模拟储层,每种储层都是由具有不同含水饱和度的气-水双相饱和砂体组成.岩石物理分析结果显示在低孔渗致密砂岩中气-水混合流体更加倾向于非均匀的斑块分布,而结合了Brie等效流体公式的Gassmann流体替换理论可以更准确地描述纵波速度随含水饱和度的变化趋势.对物理模型进行地震资料采集处理后,对比了AVO特征和叠前同步反演结果对两种砂岩储层含气饱和度预测能力的差异.AVO特征结果显示,对于混合流体均匀分布的高孔渗砂岩储层,AVO响应曲线和属性变化很难对含气饱和度进行估算;对于混合流体斑块分布的致密砂岩储层,AVO特征可以定性地分辨出储层是否为高、中、低含气情况.反演结果显示,密度及纵横波速度比分别对高孔渗及致密砂岩储层的含气饱和度有着较好的指示能力. 相似文献
12.
二氧化碳地质封存是减少温室气体排放和减缓温室效应的重要手段.二氧化碳封存的一个重要组成部分是地震监测,即用地震的方法监测封存后的二氧化碳的分布变化.为了实现这个目标,需要建立储层参数与地震性质之间的关系(岩石物理模型)和从地震监测数据中反演获得储层流体的饱和度等参数.首先,本文以Biot理论为基础,结合多相流模型研究了多个物理参数(孔隙度、二氧化碳饱和度、温度和压力等)对同时含有二氧化碳和水的孔隙介质的波速和衰减等属性的影响.结果表明:孔隙度和二氧化碳饱和度对岩石的频散和衰减属性影响强烈,而温度和压力通过孔隙流体性质对岩石的波速产生影响.然后,本文基于含多相流的Biot理论,应用抗干扰能力强、且具有更好的局部搜索能力和抗早熟能力的自适应杂交遗传算法对实际数据进行了反演研究.对岩心实验数据的反演研究表明了算法的有效性,而且表明含多相流的Biot理论能够很好地解释水和二氧化碳饱和岩石的波速特征.最后,我们将自适应杂交遗传算法应用于实际封存项目的地震监测数据,获得了封存后不同时期的二氧化碳饱和度,达到了用地震方法监测二氧化碳分布的目的. 相似文献
13.
《中国科学:地球科学(英文版)》2021,(6)
Carbonate reservoirs have complex pore structures, which not only significantly affect the elastic properties and seismic responses of the reservoirs but also affect the accuracy of the prediction of the physical parameters. The existing rockphysics inversion methods are mainly designed for clastic rocks, and the inversion objects are generally porosity and water saturation. The data used are primarily based on the elastic parameters, and the inversion methods are mainly linear approximations. To date, there has been a lack of a simultaneous pore structure and physical parameter inversion method for carbonate reservoirs. To solve these problems, a new Bayesian nonlinear simultaneous inversion method based on elastic impedance is proposed. This method integrates the differential effective medium model of multiple-porosity rocks, Gassmann equation,Amplitude Versus Offset(AVO) theory, Bayesian theory, and a nonlinear inversion algorithm to achieve the simultaneous quantitative prediction of the pore structure and physical parameters of complex porous reservoirs. The forward modeling indicates that the contribution of the pore structure, i.e., the pore aspect ratio, to the AVO response and elastic impedance is second only to that of porosity and is far greater than that of water saturation. The application to real data shows that the new inversion method for determining the pore structure and physical parameters directly from pre-stack data can accurately predict a reservoir's porosity and water saturation and can evaluate the pore structure of the effective reservoir. 相似文献
14.
CO2 saturations are estimated at Sleipner using a two-step imaging workflow. The workflow combines seismic tomography (full-waveform inversion) and rock physics inversion and is applied to a two-dimensional seismic line located near the injection point at Sleipner. We use baseline data (1994 vintage, before CO2 injection) and monitor data that was acquired after 12 years of CO2 injection (2008 vintage). P-wave velocity models are generated using the Full waveform inversion technology and then, we invert selected rock physics parameters using an rock physics inversion methodology. Full waveform inversion provides high-resolution P-wave velocity models both for baseline and monitor data. The physical relations between rock physics properties and acoustic wave velocities in the Utsira unconsolidated sandstone (reservoir formation) are defined using a dynamic rock physics model based on well-known Biot–Gassmann theories. For data prior to injection, rock frame properties (porosity, bulk and shear dry moduli) are estimated using rock physics inversion that allows deriving physically consistent properties with related uncertainty. We show that the uncertainty related to limited input data (only P-wave velocity) is not an issue because the mean values of parameters are correct. These rock frame properties are then used as a priori constraint in the monitor case. For monitor data, the Full waveform inversion results show nicely resolved thin layers of CO2–brine saturated sandstones under intra-reservoir shale layers. The CO2 saturation estimation is carried out by plugging an effective fluid phase in the rock physics model. Calculating the effective fluid bulk modulus of the brine–CO2 mixture (using Brie equation in our study) is shown to be the key factor to link P-wave velocity to CO2 saturation. The inversion tests are done with several values of Brie/patchiness exponent and show that the CO2 saturation estimates are varying between 0.30 and 0.90 depending on the rock physics model and the location in the reservoir. The uncertainty in CO2 saturation estimation is usually lower than 0.20. When the patchiness exponent is considered as unknown, the inversion is less constrained and we end up with values of exponent varying between 5 and 20 and up to 33 in specific reservoir areas. These estimations tend to show that the CO2–brine mixing is between uniform and patchy mixing and variable throughout the reservoir. 相似文献
15.
Seismic Rock physics plays a bridge role between the rock moduli and physical properties of the hydrocarbon reservoirs. Prestack seismic inversion is an important method for the quantitative characterization of elasticity, physical properties, lithology and fluid properties of subsurface reservoirs. In this paper, a high order approximation of rock physics model for clastic rocks is established and one seismic AVO reflection equation characterized by the high order approximation(Jacobian and Hessian matrix) of rock moduli is derived. Besides, the contribution of porosity, shale content and fluid saturation to AVO reflectivity is analyzed. The feasibility of the proposed AVO equation is discussed in the direct estimation of rock physical properties. On the basis of this, one probabilistic AVO inversion based on differential evolution-Markov chain Monte Carlo stochastic model is proposed on the premise that the model parameters obey Gaussian mixture probability prior model. The stochastic model has both the global optimization characteristics of the differential evolution algorithm and the uncertainty analysis ability of Markov chain Monte Carlo model. Through the cross parallel of multiple Markov chains, multiple stochastic solutions of the model parameters can be obtained simultaneously, and the posterior probability density distribution of the model parameters can be simulated effectively. The posterior mean is treated as the optimal solution of the model to be inverted.Besides, the variance and confidence interval are utilized to evaluate the uncertainties of the estimated results, so as to realize the simultaneous estimation of reservoir elasticity, physical properties, discrete lithofacies and dry rock skeleton. The validity of the proposed approach is verified by theoretical tests and one real application case in eastern China. 相似文献
16.
Seismic amplitudes contain important information that can be related to fluid saturation. The amplitude‐versus‐offset analysis of seismic data based on Gassmann's theory and the approximation of the Zoeppritz equations has played a central role in reservoir characterization. However, this standard technique faces a long‐standing problem: its inability to distinguish between partial gas and “fizz‐water” with little gas saturation. In this paper, we studied seismic dispersion and attenuation in partially saturated poroelastic media by using frequency‐dependent rock physics model, through which the frequency‐dependent amplitude‐versus‐offset response is calculated as a function of porosity and water saturation. We propose a cross‐plotting of two attributes derived from the frequency‐dependent amplitude‐versus‐offset response to differentiate partial gas saturation and “fizz‐water” saturation. One of the attributes is a measure of “low frequency”, or Gassmann, of reflectivity, whereas the other is a measure of the “frequency dependence” of reflectivity. This is in contrast to standard amplitude‐versus‐offset attributes, where there is typically no such separation. A pragmatic frequency‐dependent amplitude‐versus‐offset inversion for rock and fluid properties is also established based on Bayesian theorem. A synthetic study is performed to explore the potential of the method to estimate gas saturation and porosity variations. An advantage of our work is that the method is in principle predictive, opening the way to further testing and calibration with field data. We believe that such work should guide and augment more theoretical studies of frequency‐dependent amplitude‐versus‐offset analysis. 相似文献
17.
Bayesian inversion of time‐lapse seismic data for the estimation of static reservoir properties and dynamic property changes 
下载免费PDF全文
下载免费PDF全文 Seismic conditioning of static reservoir model properties such as porosity and lithology has traditionally been faced as a solution of an inverse problem. Dynamic reservoir model properties have been constrained by time‐lapse seismic data. Here, we propose a methodology to jointly estimate rock properties (such as porosity) and dynamic property changes (such as pressure and saturation changes) from time‐lapse seismic data. The methodology is based on a full Bayesian approach to seismic inversion and can be divided into two steps. First we estimate the conditional probability of elastic properties and their relative changes; then we estimate the posterior probability of rock properties and dynamic property changes. We apply the proposed methodology to a synthetic reservoir study where we have created a synthetic seismic survey for a real dynamic reservoir model including pre‐production and production scenarios. The final result is a set of point‐wise probability distributions that allow us to predict the most probable reservoir models at each time step and to evaluate the associated uncertainty. Finally we also show an application to real field data from the Norwegian Sea, where we estimate changes in gas saturation and pressure from time‐lapse seismic amplitude differences. The inverted results show the hydrocarbon displacement at the times of two repeated seismic surveys. 相似文献