基于BISQ模型的储层物性参数贝叶斯反演方法.pdf

2023 年 8 月第 58 卷 第 4 期基于BISQ模型的储层物性参数贝叶斯反演方法曲志鹏1，温瑨*2，韩宏伟1，步帆2，王兴谋1，朱剑兵1（1.中国石化胜利油田分公司物探研究院，山东东营 257022；2.清华大学数学科学系，北京 100084）摘要：对于物性参数反演这类存在多解性的问题，常见的方法往往难以高效、高精度地获得最优解。因此，根据贝叶斯统计学，采用高斯先验模型引入测井数据约束，并发展储层物性参数贝叶斯反演方法，以克服反演的多解性。其次，在 BISQ 理论和等效流体理论的基础上建立了岩石物理正演模型，该模型反映了孔隙介质中固体和流体的多种物理作用，可精细地刻画储层中地震波的传播过程，进而提升反演结果的可靠性。此外，由于确定性优化算法通常只有局部收敛能力，无法收敛到反演问题的全局最优解，所以引入具有全局收敛能力的杂交遗传算法，以提高反演方法的精度和效率。最终，利用所提反演方法预测中国东部 A 油田目标区域的孔隙度和含水饱和度，获得了较好效果。关键词：BISQ模型，贝叶斯反演，储层弹性参数，储层物性参数中图分类号：P631 文献标志码：A doi：10.13810/ki.issn.1000-7210.2023.04.019Bayesian inversion based on BISQ model for reservoir physical propertiesQU Zhipeng1，WEN Jin2，HAN Hongwei1，BU Fan2，WANG Xingmou1，ZHU Jianbing1（1.Shengli Geophysical Research Institute of SINOPEC，Dongying，Shandong 257022，China；2.Department of Mathematical Sciences，Tsinghua University，Beijing 100084，China）Abstract：For problems with multiple solutions such as the physical property inversion，common methods are usually difficult to obtain optimal solutions with high efficiency and accuracy.Therefore，according to the Bayesian statistics，we use the Gaussian prior model to introduce the constraints from welllog data and develop a Bayesian inversion method for reservoir physical properties to overcome the multisolution problem.Next，a petrophysical forward model is developed based on the BISQ and effective fluid theories，which can reflect the multiple physical effects of solids and fluid in the porous medium and finely characterize the propagation process of seismic waves in the reservoir，thus enhancing the reliability of inversion results.In addition，because deterministic optimization algorithms usually have only local convergence capability and fail to converge to the global optimal solution of the inversion problem，hybrid genetic algorithms with global convergence capability are introduced to improve the accuracy and efficiency of the inversion method.Finally，the proposed inversion method is used to predict the porosity and water saturation in oil field A in eastern China，and positive results are obtained.Keywords：BISQ model，Bayesian inversion，reservoir elastic parameters，reservoir physical properties曲志鹏，温瑨，韩宏伟，等.基于 BISQ 模型的储层物性参数贝叶斯反演方法 J.石油地球物理勘探，2023，58（4）：942948.QU Zhipeng，WEN Jin，HAN Hongwei，et al.Bayesian inversion based on BISQ model for reservoir physical properties J.Oil Geophysical Prospecting，2023，58（4）：942948.0 引 言储层物性参数是描述储层特征的关键信息，是预测储层的含油气性、估计油气储量等的基础1。在实际工作中，通常从两个方面预测储层物性参数。一方面，仅根据有限的钻井数据和（或）测井资料估算储层物性参数，无法描述整个研究区的地质特征。另一方 综合研究 文章编号：1000-7210（2023）04-0942-07*北京市海淀区清华大学数学科学系，100084。Email： 本文于 2022年 7月 18日收到，最终修改稿于 2023年 5月 10日收到。本项研究受中国石化重点科技攻关课题“储层物性定量预测与精细评价方法研究”（P21018）资助。第 58 卷 第 4 期曲志鹏，等：基于 BISQ模型的储层物性参数贝叶斯反演方法面，虽然采集的大量地震资料蕴含了丰富的地下信息，但由地震数据直接反演物性参数往往是困难的。为此，人们引入弹性参数作为中间变量，发展了一系列物性参数的间接地震反演方法2。这种方法需要建立一个正演模型描述弹性和物性参数之间的关系，进而求解对应的反问题预测物性参数。针对物性参数预测问题，统计学方法是一种容易理解且方便的方法1，3，但需要大量高质量的测井或岩心数据，具有明显的数据依赖性，而且缺乏物理过程模拟，难以反映普遍规律。因此，引入考虑实际物理过程的岩石物理模型作为正演模型，可以有效弥补纯统计学方法的不足45。自 20世纪中叶以来，人们针对实际油气储层，发展了一系列岩石物理模型描述弹性参数和物性参数之间的关系，并模拟了储层介质中的地震波传播现象。Biot/Squirt（BISQ）模型6 是典型的岩石物理模型，起源于20世纪中叶对孔隙介质中的固液耦合机制的研究。Biot7 基于数理方程和拉格朗日力学理论，率先指出孔隙双相介质存在一种动态的固流耦合机制Biot机制，并建立了相应的弹性波方程。然而，众所周知，Biot模型并不能充分描述在某些储层中地震波的高频散和强衰减现象8。为此，Dvorkin等6 基于 Biot模型引入喷射流机制，建立了最初的各向同性介质BISQ模型。Yang等910 系统研究了BISQ模型，建立了包含两种力学机制（固流耦合各向异性效应和孔隙流体流动）的各向异性BISQ模型。近年来，通过引入等效介质、黏弹性等理论描述复杂油气储层特征1115。此外，大量的岩心实验数据和油田实际数据进一步证明了BISQ模型的可靠性和有效性1617。但是，在实际应用中基于岩石物理模型的反演易受噪声干扰，甚至得到完全错误的预测结果4。在地球物理反演领域，贝叶斯理论具有重要的理论和应用价值18，包括量化地质结构预测的不确定性、提高地震反演精度、添加空间约束等。针对物性参数预测问题，贝叶斯反演方法也备受关注。一方面，由于贝叶斯反演是基于统计学框架的反演方法，在反演目标参数的同时给出模型参数的后验概率，进而分析预测结果的不确定性4。但是，由于地下模型的特殊性，无法通过反复试验验证物性参数的不确定性估计，而且也缺乏可靠的统计检验。另一方面，根据最大后验估计理论，贝叶斯反演通过引入先验信息实现物性参数反演的正则化，进而提高了物性参数反演精度和抗噪能力19。物性参数反演经常涉及多参数优化问题，其中的最优解就是物性参数预测结果。因此，稳定、高效的优化算法也是物性参数反演的重要研究内容。如单纯形法、最速下降法、共轭梯度法等确定性算法是求解线性优化问题的常用算法。然而，对于物性参数反演这类存在多解性的优化问题，上述确定性算法往往难以找到全局最优解。因此，如遗传算法、模拟退火算法等具有全局搜索能力的随机算法通常是更好的选择，其中遗传算法的应用最广泛。Eiben 等20 基于马尔科夫链理论，在数学上证明了遗传算法具全局收敛性，并成功应用于地球物理反问题2122。Fang等16 基于BISQ模型利用混合遗传算法可有效反演储层物性参数。本文旨在发展一种间接地震反演方法，进而可靠地预测储层物性参数，其核心是根据储层弹性参数预测多种物性参数的联合非线性反演。虽然现有方法取得了一定效果，但仍存在几个问题亟待解决。一是建立可靠的岩石物理正演模型刻画弹性参数和物性参数的非线性关系。二是在岩石物理模型的基础上引入测井数据约束，以提高反演精度。为此，本文首先在BISQ理论和等效流体理论的基础上建立了岩石物理正演模型，该模型反映了孔隙介质中固体和流体的多种物理作用，可精细地刻画储层中地震波的传播过程。其次，在贝叶斯统计学的框架下，采用高斯先验模型，引入测井数据约束，进而发展了基于BISQ模型的储层物性参数贝叶斯反演方法。再次，由于反演问题的多解性，引入具有全局搜索能力的杂交遗传算法提高反演精度和效率。利用所提反演方法预测中国东部A油田目标区域的孔隙度和含水饱和度，获得了较好效果。1基于 BISQ 模型的岩石物理建模将孔隙流体等效模型和固体骨架模量孔隙度模型引入部分饱和的 BISQ模型11，建立相应的岩石物理模型刻画储层弹性参数与物性参数之间的非线性关系。储层介质中的孔隙流体通常是石油、天然气和地层水的混合物。由于混合成分是气体或液体，其剪切模量为零。在液体充分混合的情况下，通过 Reuss平均公式8 943石 油 地 球 物 理 勘 探2023 年1Kf=swKw+1-swKo（1）计算油水混合的等效流体体积模量Kf。式中：sw为含水饱和度；Kw和Ko分别为水和油的体积模量。当孔隙连通性较差时，多孔岩石被分为一系列的斑块，此时孔隙流体分布不均匀。在这种复杂情况下，虽然简单的平均公式不能精确地描述液体混合物的等效作用，但是仍可用式（1）计算Kf下界Kdf。此外，Voigt平均公式8 Kuf=swKw+(1-sw)Ko（2）给出了Kf上界Kuf。在实际应用中，因为Kuf和Kdf是接近的，所以可将两者视为Kf的有效近似值。固体骨架是一种理想的等效模型，而不是真实的饱和气体岩石。固体骨架模量只与固体组分及其孔隙特性相关，不受孔隙流体弹性模量影响8。一般情况下，人们根据现有测井资料或岩心资料建立适当的数学模型近似地计算固体骨架模量4，16。Dvorkin等6 结合含流体孔隙岩石在宏观尺度的Biot流作用和在微观尺度的喷射流作用建立了BISQ模型。Biot流特征频率为fb=2f（3）式中：为流体黏度；为孔隙度；为渗透率；f为孔隙流体密度。由于实际地震波频率f远小于fb（fb106 Hz），Dvorkin等23 基于 BISQ 模型提出了更简洁的 P波速度公式VP=1Re(Y)（4）其中 Y=s(1-)+fM+Fsq2Fsq=F 1-2J1()J0()=Ri2fF=1-KKs（5）式中：J0和J1分别为 0阶和 1阶 Bessel函数；s为固体基质密度；R为特征喷射流长度6；F为Biot流动系数；K 为固体骨架的体积模量；Ks为固体基质的体积模量；M为固体骨架的单轴应变模量，即 M=K+43G（6）式中G为固体骨架的剪切模量。由下式7 1F=1Kf+-Ks（7）计算F。在BISQ模型中，S波速度为11 VS=GRe(x)x=(1-)s+(1-)f 1=a/f+ifbf（8）式中a为固流耦合密度7。根据质量守恒定律得到含流体岩石密度和f=f+(1-)s f=sww+(1-sw)o（9）式中w和o分别为地层水和油的密度。2储层物性参数的贝叶斯估计为方便书写，基于前述BISQ模型，可以将密度和纵、横波阻抗IP、IS视为关于模型参数l（待反演物性参数）和辅助参数m的非线性函数 =f1(l，m)IP=VP=f2(l，m)IS=VS=f3(l，m)（10）其中 l=(，sw)Tm=(s，w，o，Ks，Kw，Ko，)（11）在实际应用中，密度和弹性阻抗作为观测数据，可以视为模型计算的理论值和观测误差之和 =f1(l，m)+e1IP=f2(l，m)+e2IS=f3(l，m)+e3（12）式中 ei（i=1，2，3）为观测误差。设观测数据d=(，IP，IS)，将式（12）改为更简洁的向量形式d=f(l，m)+e（13）在贝叶斯统计框架下，通过地震属性的条件分布944第 58 卷 第 4 期曲志鹏，等：基于 BISQ模型的储层物性参数贝叶斯反演方法的概率密度p(d|l，m)和模型参数的先验分布概率密度p(l)计算后验分布概率密度p(l|d，m)=p(d|l，m)p(l)p(d|m)（14）式中p(d|m)可视为p(l|d，m)的概率归一化常数p(d|m)=p(d|l，m)p(l)dl（15）最大后验（MAP）估计是贝叶斯错误率最小的预测值。根据式（14）可以计算目标参数的MAP估计 l=arg maxl Rnp(l|d，m)=arg maxl Rnln p(l|d，m)=arg maxl Rnlnp(d|l，m)p(l)p(d|m)=arg maxl Rn ln p(d|l，m)+ln p(l)-ln p(d|m)式中：Rn为n维实数空间，n为待反演的物性参数个数；ln p(d|m)和l无关，在反演中可以忽略。假设式（13）的随机误差 e 服从高斯分布 N（i，2i），因此p(d|l，m)=C exp-12i=13wi di-fi(l，m)-i2（17）式中：i为平均偏差；wi=1/2i为权重，2i为方差；C为p(d|l，m)的概率归一化常数。通常情况下，p（l）作为已知信息约束目标参数的反演，从而有效地避免物性参数反演的多解性并抑制噪声。高斯分布是一种描述先验信息的简单且有效的统计模型19，其概率密度为p(l)=C0exp -12D2(l)（18）式中：C0为p(l)的概率归一化常数；D(l)为模型参数的马氏距离D(l)=(l-0)TW0(l-0)（19）式中：W0=-10为模型参数先验分布协方差矩阵的逆；0为先验分布的期望。实际问题中，模型参数空间通常是有界的，即lL（L表示模型参数所属的有界空间）。因此，p(l)应该是截断的高斯分布，而不是定义在无界空间的高斯分布。根据式（18），通过修正C0可获得正确的先验分布。将式（17）、式（18）代入式（16），得到关于模型参数的正则化反演问题l=arg maxl Rni=13wi di-fi(l，m)-i2+D2(l)（20）C和C0与模型参数无关，可以在反演中忽略。BISQ 模型是一种非线性岩石物理模型，所以式（20）呈多局部最优解。然而，传统的基于梯度下降的反演算法通常强烈依赖初值，因此全局收敛能力弱。同时，由于BISQ模型的复杂性，梯度往往难以计算。为此，文中引入一种具有全局收敛能力的杂交遗传算法16，该算法模仿了生物进化机制，将优化问题的候选解视为生物种群中的个体，并根据目标函数量化种群个体的生存能力（适应值）。该方法主要包括选择、交叉和变异运算。在遗传算法运行过程中，种群整体向适应值高的方向进化，而不同个体在不同的局部最优解附近聚集，进而找到全局最优解。3实际应用实际数据来自中国东部A油田，图1为 A井、B井及 C井测井数据。由图可见，A1、A2和 C1主要为含油层，B1主要为含水层。采用叠前AVO反演24 得到、VP和 VS（图 2）。首先采用本文提出的贝叶斯反演方法初步反演，再联合反演储层区域（0.1）的 和sw。在初步反演 的过程中，的先验分布为正态分布N（0，2103）。图3为C井测井数据直方图和先验分布。在对储层区域的反演过程中，根据物性参数测井数据得到 0.1 时的 和 sw的期望和协方差矩阵 0=(15 90)T 0=()2.7 10-35.0 10-35.0 10-31.5 10-2图 4为 和 sw的先验分布。由图可见，先验分布呈截断的高斯分布，其中 的变化范围为 0，0.3，sw的变化范围为 0，1。在两步反演过程中，各组分密度和体积模量数据见表 1。根据文献 8，19，25 选取其他辅助参数，并根据C井物性参数和井旁道弹性参数（图5）进行了校准：/=1 Pasm-2，f=30 Hz，a=400 kg/m3，R=2 mm。选择C井为校准井，给出固体骨架模量的计算公式(16)945石 油 地 球 物 理 勘 探2023 年K=233.32-53.06+16.54G=84.452-9.558+5.315（21）将 C井的物性参数数据代入式（13）获得观测误差样本。根据观测误差样本关于正态分布的分位数分位数图（图 6）分析样本的正态性。可见，若样本数据点（蓝色）整体越接近参考值（红色），则表示样本分布越接近正态分布。由于样本整体不满足正态分图 1A井（a）、B井（b）、C井（c）测井数据DT为声波曲线，SP为自然电位曲线，GR为伽马曲线。在综合解释中，黄色为含油层，绿色为含水层，蓝色为非储层。图 2弹性参数剖面图 3C井 测井数据直方图和先验分布 红线表示 的先验概率密度。图 4和 sw的先验分布（储层区域）946第 58 卷 第 4 期曲志鹏，等：基于 BISQ模型的储层物性参数贝叶斯反演方法布假设，所以在反演中仅使用储层区域（0.1）的密度数据。根据观测误差样本，得到平均偏差为0，IP平均偏差为-0.6 106 kgs-1m-2，IS平均偏差为-0.2106 kgs-1m-2，并给定权重为310-2(kg-1m-3)-1、IP权 重 为2 10-6(kgs-1m-2)-1，IS权 重 为310-6(kgs-1m-2)-1。在遗传算法方面，选取迭代次数为20，种群个体数为20，交叉概率为0.9，变异概率为0.1。图 7为 和 sw反演结果。由图可见：含油区域（A1、A2和 C1）呈高、低 sw特征，含水区域（B1）呈高、高sw特征。与图2相比，图7的储层和非储层的差别更明显，储层位置和形态更清晰；与图2相比，图7的含油区域和含水区域的sw差异更明显，能有效识别流体类别。4结束语本文基于BISQ模型提出了一种新的物性参数贝叶斯反演方法，其中 BISQ模型描述了储层物性参数和弹性参数之间的确定性物理关系。在贝叶斯统计理论框架下，采用观测误差和先验概率模型刻画实际问题存在的随机性，并为基于 BISQ模型的反演问题引入正则化约束，从而提高了反演精度。此外，引入的杂交遗传算法可以有效改善基于复杂岩石物理模型的物性参数反演的多解性。实际应用结果表明，所提反演方法具有广阔的应用前景。参 考 文 献1 DOYEN P M.Seismic Reservoir Characterization:An Earth Modelling PerspectiveM.EAGE Publications,2007.2 韩宏伟,程远锋,张云银,等.储层物性的地震预测技术综述J.地球物理学进展,2021,36(2):595610.HAN Hongwei,CHENG Yuanfeng,ZHANG Yunyin,et al.Review of seismic prediction of reservoir geophysical propertiesJ.Progress in Geophysiscs,2021,36(2):595610.3 魏国华,韩宏伟,刘浩杰,等.基于半监督高斯混合模型与梯度提升树的砂岩储层相控孔隙度预测J.石油地球物理勘探,2023,58(1):4655.表 1各组分密度和体积模量组分固体基质地层水石油体积模量GPa39.002.601.65密度 gcm-32.9501.1450.945图 5C井物性参数曲线及其井旁道弹性参数Ip、Is的量纲均为106 kg s-1m-2，的量纲为g cm-3图 6观测误差样本关于正态分布的分位数分位数图左侧为全体样本，右侧为储层区域样本。图 7和 sw反演结果947石 油 地 球 物 理 勘 探2023 年WEI Guohua,HAN Hongwei,LIU Haojie,et al.Faciescontrolled porosity prediction of sandstone reservoirs based on semisupervised Gaussian mixture model and gradient boosting treeJ.Oil Geophysical Prospecting,2023,58(1):4655.4 BACHRACH R.Joint estimation of porosity and saturation using stochastic rockphysics modelingJ.Geophysics,2006,71(5):O53O63.5 田军,刘永雷,徐博,等.深埋储层孔隙度迭代反演方法J.石油地球物理勘探,2022,57(3):666675.TIAN Jun,LIU Yonglei,XU Bo,et al.A method for porosity prediction of deeply buried reservoirs based on iterative inversionJ.Oil Geophysical Prospecting,2022,57(3):666675.6 DVORKIN J,NUR A.Dynamic poroelasticity:A unified model with the squirt and the Biot mechanismsJ.Geophysics,1993,58(4):524533.7 BIOT M A.Theory of propagation of elastic waves in a fluidsaturated porous solid.I.Lowfrequency rangeJ.Journal of the Acoustical Society of America,1956,28(2):168178.8 MAVKO G,MUKERJI T,DVORKIN J.The Rock Physics Handbook(2nd Edition)M.Cambridge University Press,2009.9 YANG D,ZHANG Z.Effects of the Biot and the squirtflow coupling interaction on anisotropic elastic wavesJ.Chinese Science Bulletin,2000,45(23):21302138.10 YANG D,ZHANG Z.Poroelastic wave equation including the Biot/squirt mechanism and the solid/fluid coupling anisotropyJ.Wave Motion,2002,35(3):223245.11 NIE J,YANG D,YANG H.Wave dispersion and attenuation in partially saturated sandstonesJ.Chinese Physics Letters,2004,21(3):572575.12 YANG L,YANG D,NIE J.Wave dispersion and attenuation in viscoelastic isotropic media containing multiphase flow and its applicationJ.Science China Physics,Mechanics&Astronomy,2014,57(6):10681077.13 ZHANG B,YANG D,CHENG Y,et al.A unified poroviscoelastic model with mesoscopic and microscopic heterogeneitiesJ.Science Bulletin,2019,64(17):12461254.14 YANG J,YANG D,HAN H,et al.A wave propagation model with the Biot and the fractional viscoelastic mechanismsJ.Science China Earth Sciences,2021,64(3):364376.15 孙卫涛,熊繁升,曹宏,等.非牛顿流体孔隙介质弹性波动方程J.石油地球物理勘探,2021,56(3):519527.SUN Weitao,XIONG Fansheng,CAO Hong,et al.Elastic wave equation for porous media saturated with nonNewtonian fluidJ.Oil Geophysical Prospecting,2021,56(3):519527.16 FANG Z,YANG D.Inversion of reservoir porosity,saturation,and permeability based on a robust hybrid genetic algorithmJ.Geophysics,2015,80(5):R265R280.17 NIE J,YANG D,YANG H.Inversion of reservoir parameters based on the BISQ model in partially saturated porous mediumJ.Chinese Journal of Geophysics,2004,47(6):12411246.18 BOSCH M,MUKERJI T,GONZALEZ E F.Seismic inversion for reservoir properties combining statistical rock physics and geostatistics:A reviewJ.Geophysics,2010,75(5):75A16575A176.19 WEN J,CHENG Y,YANG D,et al.Joint twostep inversion of porosity and saturation based on the BISQ model with mutliphase fluidsJ.Geophysics,2022,87(4):M99M110.20 EIBEN A E,AARTS E,VAN HEE K M.Global convergence of genetic algorithms:a Markov chain analysisC.International Conference on Parallel Problem Solving from Nature,1990,312.21 周琦,印兴耀,李坤.煤系地层地震岩石物理建模及横波预测方法J.石油地球物理勘探,2022,57(2):357366.ZHOU Qi,YIN Xingyao,LI Kun.Seismic rock physics modeling and shear wave velocity prediction method of coal measure strataJ.Oil Geophysical Prospecting,2022,57(2):357366.22 顾雯,印兴耀,巫芙蓉,等.波形驱动下多参数约束高分辨率反演方法以四川盆地渝西地区龙马溪组页 岩 气 为 例 J.石 油 地 球 物 理 勘 探,2021,56(6):13111321.GU Wen,YIN Xingyao,WU Furong,et al.Multiparameter constrained highresolution inversion method driven by waveform:A case study of Longmaxi Formation shale gas in western Chongqing,Sichuan BasinJ.Oil Geophysical Prospecting,2021,56(6):13111321.23 DVORKIN J,NOLENHOEKSEMA R,NUR A.The squirtflow mechanism:Macroscopic descriptionJ.Geophysics,1994,59(3):428438.24 RUSSELL B,HAMPSON D P,BANKHEAD B.An inversion primerJ.CSEG Recorder,2006,31(Special):96103.25 张永刚.复杂介质地震波场模拟分析与应用M.北京:石油工业出版社,2007.（本文编辑：刘勇）作 者 简 介曲志鹏研究员，1981 年生；2004 年获中国石油大学（华东）地球科学与技术学院矿产普查与勘探专业学士学位；现就职于中国石化胜利油田分公司物探研究院，长期从事油气地球物理勘探工作。948