面向页岩气储层的叠前多参数地震反演方法研究
时间：2023-03-10 08:15:05 来源：柠檬阅读网 本文已影响 人 
陈 勇,孙振涛,许 凯
(中石化石油物探技术研究院有限公司,江苏南京211103)
页岩气是蕴藏于页岩地层中可供开采的天然气资源,中国的页岩气可采储量较大,具有良好的勘探、开发潜力。随着油气勘探开发逐渐由常规油气向非常规油气尤其是页岩气发展,以及逐步走向精细化,储层精细描述至关重要。地震反演技术可揭示地下油气藏的精细分布特征,因此受到广泛关注[1]。页岩气储层的水平薄互层构造和水平层理具有明显的VTI各向异性特征[2],大量研究表明,页岩气储层的水平层理缝引起的各向异性影响地震反射特征[3-7],裂缝弱度参数作为非常规页岩气储层水平层理缝表征的关键参数[8],有助于页岩气储层各向异性特征的识别与描述,因此十分有必要开展基于VTI等效介质的叠前多参数反演方法研究。
建立VTI等效介质弹性参数及各向异性参数随地震振幅的变化关系是地震反演的基础。针对VTI等效介质,众多学者推导了反射系数和透射系数的精确表达式[9-11],由于方程形式复杂,限制了VTI等效介质反射系数精确方程的实际应用。为此,众多学者针对不同类型的储层特征参数推导出相应的VTI等效介质反射系数近方程[12-15]。THOMSEN[12]基于弱各向异性假设,推导了VTI等效介质的P波反射系数线性化近似方程。URSIN等[13]基于弱阻抗差和弱各向异性假设,推导了VTI等效介质的反射系数和透射系数线性化近似方程。RÜGER[14]进一步完善了THOMSEN[12]提出的P波反射系数近似方程,以适应大入射角度的情况,该公式即为经典的VTI等效介质反射系数方程。对于具有任意对称性的弱各向异性介质,SHAW等[16]基于稳相法推导出了PP波和PS波反射系数方程。本文基于波恩积分法,推导了VTI等效介质纵波反射系数近似方程,构建了储层弹性模量、裂缝弱度参数与地震反射系数之间的量化关系,为后续叠前反演奠定基础。
稳定、可靠的地震反演是页岩气储层精细表征的关键。贝叶斯反演理论常被用于求解多元目标函数[17-27]。从目前的研究趋势看,地震反演方法研究主要包括叠前AVO反演和全波形反演。基于精确Zoeppritz方程的地震反演方法能够利用丰富的叠前地震信息预测模型参数,但其非线性强、稳定性差,全波形反演方法则计算量大、耗时长,反演尺度和计算效率不能满足储层精细表征要求。针对VTI等效介质理论模型,传统的近似反射系数方程为包含弹性参数和各向异性参数的五项模型参数表达式,由于反演参数较多,造成的反演病态问题限制了方法的应用[28-33],且未考虑裂缝诱导的各向异性效应。
针对传统VTI等效介质的叠前多参数反演存在的问题,本文开展了面向页岩气储层的叠前多参数反演方法研究,提出了一种稳健的基于VTI等效介质AVA反演的储层弹性模量与裂缝弱度参数评价方法。首先,基于散射理论以及弱各向异性近似假设条件,推导出含纵、横波模量、密度和裂缝弱度参数的VTI等效介质纵波反射系数近似方程,经过整合与简化,获得含4个组合属性参数(A、B、C、D属性)的线性化纵波反射系数近似方程。进一步,基于贝叶斯反演理论提出了VTI等效介质四项组合属性参数反演方法,实现了组合属性参数的稳定、可靠预测。然后,通过建立组合属性参数与储层参数之间的内在关系,获得弹性模量和裂缝弱度等储层参数,实现页岩气储层弹性特征与裂缝发育特征的地震准确描述。最后,利用VTI等效介质反射系数方程合成AVA叠前道集和实际数据,验证所提出方法的正确性和有效性。
针对页岩气储层精细表征与描述问题,开展基于VTI等效介质AVA反演的储层弹性模量与裂缝弱度参数反演方法研究。首先基于推导的五项式VTI介质纵波反射系数近似方程,经过四项组合属性参数表征的新VTI介质纵波反射系数方程,以减少反演的模型参数,提高反演的稳定性;进一步结合贝叶斯反演理论以及属性之间的数学关系,实现储层弹性模量和裂缝参数的稳定、可靠反演,有效刻画页岩气储层的弹性特征和裂缝发育特征。
1.1 含裂缝弱度的VTI等效介质P波反射系数线性化推导
VTI对称介质存在5个独立的弹性刚度分量,其6×6的弹性刚度矩阵CVTI可以表示为:
CVTI=
(1)
式中:C11,C13,C33,C44,C66分别为弹性刚度分量。
在VTI等效介质的弱各向异性假设下,Thomsen简化了上述复杂关系式,定义了5个弹性常数:
(2)
式中:vP,vS和ρ分别为各向同性背景部分的纵、横波速度和密度;ε,δ,γ为弱各向异性参数,也称为Thomsen各向异性参数。根据弱各向异性近似假设条件,RÜGER[14]推导了VTI等效介质纵波反射系数近似方程,具体形式如下:
(3)
式中:θ为入射角;Δε,Δδ分别为上下两层介质的Thomsen弱各向异性参数差值。
根据线性滑移理论,裂缝诱导的VTI等效介质刚度系数可以表示为[8]:
(4)
式中:M和μ分别表示纵波、横波弹性模量,M=λ+2μ,ξ=λ/M;δN和δT分别表示法向和切向裂缝弱度参数。
基于散射理论和玻恩近似,在弱各向异性假设条件下,推导出基于VTI等效介质的纵波反射系数五项式近似方程(具体推导过程见附录A):
(5)
式中:ΔM/M和Δμ/μ是纵、横波弹性模量的反射率;Δρ/ρ是密度反射率;g是纵横波速度比。
为了减少待反演模型参数的数量,将公式(5)进一步改写为如下形式(具体简化过程见附录B):
(6)
在界面两侧弹性特征差异较小的情况下,对公式(6)进行BORTFELD[34]和OLDENBURG等[35]提出的近似,也即:
(7)
式中i+1和i分别表示第i+1层和第i层。将公式(7)代入公式(6),得到如下表达式:
(8)
设:
(9)
则有以下表达式:
(10)
公式(10)被称为“含四项组合属性参数的新方程”,新属性A,B,C和D是用于VTI等效介质的弹性模量和裂缝弱度预测的四项模型组合参数,A=Mρe-δN,B=μe-1/2gΔδN-1/2ΔδT,C=(Me-δN)/ρ,D=eδN,A为纵波模量、密度与裂缝弱度参数相组合的属性,B为剪切模量与裂缝弱度参数相组合的属性,C为剪切模量与裂缝弱度参数、密度相组合的属性,D为指数的裂缝弱度参数。4个属性的物理意义可以解释为:A和C属性代表弹性性质和裂缝弱度的项,两者相除与密度相关,两者的乘积与纵波模量和裂缝弱度有关;在裂缝弱度较小的情况下,属性B近似等于剪切模量;D属性与裂缝弱度直接相关。
设计两层反射系数模型(表1),对比近似公式与精确公式的纵波反射系数。反射系数误差计算公式为:error=|R精确-R近似|/R精确,由图1可算出,在入射角达到30°时,反射系数计算误差为3.26%,可见近似公式的精度完全满足反演的要求。
表1 双层反射系数模型参数
图1 近似公式与精确公式反射系数
1.2 基于贝叶斯反演理论的VTI等效介质新属性参数估计
根据褶积理论,地震振幅响应可表示为如下矩阵形式:
d=Gm
(11)
式中:d是观测数据;模型向量m=[A,B,C,D]T,T表示转置;正演算子G=WP,其中,W是子波矩阵;P是敏感度矩阵,具体表现形式如下:
P=
(12)
式中:θ1,θ2,…,θm表示m个入射角。基于贝叶斯定理,密度p(m|d)与先验概率密度p(m)和p(d|m)的乘积成正比[17]:
p(m|d)∝p(m)p(d|m)
(13)
当观测地震资料包含高斯随机噪声时,似然函数p(d|m)可表示为:
(14)
(15)
式中:Cm为模型参数的方差。根据贝叶斯理论,求解反演问题可以转化为求解观测数据与拟合数据之间的残差最小值,即:
(16)
(17)
(18)
因此,m的解可表示为[23,36-37]:
(19)
式中:μh是模型参数扰动的权重。
2.1 合成数据算例
对实际测井曲线进行预处理后,首先根据构建的页岩储层岩石物理模型,以矿物组分、孔隙度、含水饱和度等储层参数为输入,通过岩石物理正演,获得弹性参数和裂缝弱度参数(图2),通过计算进一步获得待估新属性参数(图3)。
图2 测井曲线a 纵波速度; b 横波速度; c 密度; d 法向弱度参数δN; e 切向弱度参数δT
图3 四项新属性组合参数测井曲线a A属性; b B属性; c C属性; d D属性
利用雷克子波(主频为30Hz)与由公式(10)计算的反射系数进行褶积得到合成叠前角度道集(图4)。为了验证反演方法在含噪情况下的可靠性和稳定性,在无噪合成地震数据上添加信噪比分别为5和2的高斯随机噪声,采用频谱估算法计算随机噪声,如图4a和图4b所示,叠前道集的入射角范围为3°～30°,每隔3°为一道,共有10道。反演结果如图5所示,从图5中可看出,在信噪比分别为5和2的情况下,反演值(黑色)与真实值(红色)之间均具有良好的一致性。因此,验证了本文所提出的反演方法的正确性和有效性。
图4 合成角道集a 信噪比为5; b 信噪比为2
图5 2种信噪比下A属性、B属性、C属性和D属性的真实曲线(红色)和反演曲线(黑色)(蓝色曲线表示通过平滑真实曲线生成的初始模型)a 信噪比为5; b 信噪比为2
2.2 实际数据反演
过井地震剖面如图6所示,X井位于CDP为2122处,在采样时刻2200～2400ms位置钻遇含气储层。基于实际工区构建的页岩岩石物理模型,通过岩石物理正演获得表征储层弹性特征和裂缝特征的参数,测井曲线如图7所示,在采样时刻2235ms处,速度和密度的值相对较低,裂缝弱度的值相对较高,结合测井解释成果,明确该处为含气储层。利用上述测井曲线进一步计算待反演的组合属性参数,为后续反演奠定模型基础。
图6 工区的叠后地震剖面(红线表示CDP2122处X井穿过目标储层)
图7 位于CDP 2122处X井的测井曲线a 纵波速度; b 横波速度; c 密度; d 法向裂缝弱度参数; e 切向裂缝弱度参数
图8为基于贝叶斯反演理论获得的四项组合属性参数的反演剖面。叠合测井曲线显示,可看出反演结果与测井曲线吻合度高、一致性好。为了进一步对比反演方法的效果,抽取过井处的单道反演结果,开展井震对比与分析(图9),可看出在目的层(采样时刻为2235ms处)反演的组合属性参数(红线)与测井真实值(黑线)吻合度高。为了进一步验证反演的合理性和稳定性,对比CDP为2122处的原始道集、经预处理后的叠前道集和利用估计的新属性参数合成道集(图10),可看出三者具有很好的一致性。
图8 新参数属性反演剖面(剖面上曲线表示相应的测井曲线)a A属性; b B属性; c C属性; d D属性
图9 各属性反演结果(红线)与测井曲线(黑线)比较(蓝线表示通过平滑测井曲线生成的初始模型)a A属性; b B属性; c C属性; d D属性
图10 原始道集与反演合成道集对比a CDP 2122处的原始角道集; b 经K-L变换处理后的拉平原始道集; c 采用反演的新属性A、B、C和D合成的角道集
在获得四项组合属性参数反演剖面的基础上,首先预测纵波模量和剪切模量,通过对比反演属性B(红线)、实际剪切模量(黑线)和实际测井资料计算获得的属性B(蓝线)(图11a),可看出属性B(蓝线)和真实剪切模量(黑线)之间差异较小,因此,反演属性B可近似表征各向同性背景介质的剪切模量。进一步利用反演属性A,C和D预测纵波模量,计算A与C的乘积为A·C=M2e-2δN,代入由属性D计算获得的δN,则可计算得到纵波模量,对比反演(红线)和实际(黑线)的纵波模量(图11b)可见,两者吻合度高。进一步根据反演属性A和属性C,计算过渡属性(A/C=ρ2),开方即可获得密度预测结果,对比真实值(黑线)与反演值(红线)(图11c)可见,两者趋势一致。利用反演获得的D属性预测裂缝弱度,对比反演(红线)和实际(黑线)裂缝弱度(图11d)也发现,两者一致性较好。
图11 叠前反演结果与实际值的对比a 属性B的反演(红线)、计算(蓝线)和真实(黑线)剪切模量; b 反演(红线)和真实(黑线)纵波模量; c 反演(红线)与真实(黑线)密度; d 反演(红线)与真实(黑线)裂缝弱度参数
将上述计算单井的方法应用于地震数据,利用反演获得的属性A和属性C(图8a和图8c)计算纵波模量和密度(图12a和图12b);利用反演获得的属性D(图8d)计算裂缝弱度(图12c)。从反演剖面可看出,目的层高产储层的纵波模量和密度表现为低值异常,裂缝弱度表现为高值异常,与测井解释结论一致,同时X井处的预测结果与测井曲线具有良好的一致性,充分验证了该方法的合理性和有效性。
图12 过井剖面的反演结果(曲线表示相应的测井数据)a 纵波模量; b 密度; c 裂缝弱度δN
针对页岩气储层精细表征与描述问题,提出了一种稳定的基于VTI等效介质AVA反演的储层弹性模量与裂缝弱度参数评价方法,并利用合成数据和实际资料进行了测试与应用,得出如下结论和认识。
1) 基于五项式VTI等效介质纵波反射系数近似方程,构建了四项组合属性参数表征的新VTI等效介质纵波反射系数近似方程,有效减少了反演的模型参数,为VTI等效介质叠前AVA反演奠定基础。
2) 基于贝叶斯反演理论及属性之间的数学关系,实现储层弹性模量和裂缝弱度参数的稳定、可靠反演,有效刻画页岩气储层的弹性特征和裂缝发育特征。
3) 合成数据和实际资料的测试与应用显示了良好的效果,因此本文提出的方法在页岩气的勘探开发中具有广泛的潜在应用前景。
需要指出的是,本文的方法基于VTI等效介质理论假设,实际储层中常发育高角度裂缝,因此仍需进一步研究基于正交各向异性介质的叠前多参数反演方法。
致谢:感谢中南大学潘新朋教授在论文撰写过程中的帮助。
猜你喜欢 反射系数纵波测井 本期广告索引测井技术(2022年3期)2022-11-25花岗岩物理参数与纵波波速的关系分析矿产与地质(2022年2期)2022-09-06自由界面上SV波入射的反射系数变化特征*石油管材与仪器(2022年3期)2022-06-10垂直发育裂隙介质中PP波扰动法近似反射系数研究科学技术与工程(2021年28期)2021-10-20多道随机稀疏反射系数反演石油物探(2020年6期)2020-11-25氮化硅陶瓷的空气耦合超声纵波传播特性研究陶瓷学报(2015年4期)2015-12-17基于测井响应评价煤岩结构特征中国煤层气(2015年4期)2015-08-22随钻电阻率测井的固定探测深度合成方法西南石油大学学报(自然科学版)(2015年4期)2015-08-20变截面阶梯杆中的纵波传播特性实验浙江大学学报（工学版）(2015年6期)2015-03-01中石油首个全国测井行业标准发布中国质量与标准导报(2015年2期)2015-02-28
