基于杨氏模量、泊松比和各向异性梯度的裂缝介质AVAZ反演方法
杜炳毅;杨午阳;王恩利;张广智;高建虎
【摘 要】为了全面地描述裂缝介质的储层特征、流体性质和各向异性特征,提出了由宽方位叠前道集数据直接获取裂缝介质弹性参数和各向异性梯度参数的振幅随方位角变化(Amplitude Versus Azimuth,AVAZ)反演方法.首先推导了基于杨氏模量、泊松比和各向异性梯度的各向异性AVO方程,通过与Ruger近似进行对比,分析了该方程的近似精度.然后利用实际测井资料和二维逆掩断层模型建立了不同方位角的合成叠前角度道集,对未加噪声和信噪比为3的宽方位叠前角度道集进行了AVAZ反演方法测试,结果表明,未加噪声和信噪比为3的合成宽方位叠前道集均能反演得到符合地震反演精度的杨氏模量(E),泊松比(σ)和各向异性梯度(T,噪声对杨氏模量的反演影响较小,不含噪声时泊松比和各向异性梯度的反演精度更高,与真实值吻合程度较高.
【期刊名称】《石油物探》
【年(卷),期】2015(054)002
【总页数】8页(P218-225)
【关键词】杨氏模量;泊松比;各向异性梯度;裂缝介质;AVAZ反演
【作 者】杜炳毅;杨午阳;王恩利;张广智;高建虎
【作者单位】中国石油勘探开发研究院西北分院,甘肃兰州730020;中国石油勘探开发研究院西北分院,甘肃兰州730020;中国石油勘探开发研究院西北分院,甘肃兰州730020;中国石油大学(华东)地球科学与技术学院,山东青岛266580;中国石油勘探开发研究院西北分院,甘肃兰州730020
【正文语种】中 文
【中图分类】P631
Thomsen[1]提出了弱各向异性介质理论,引入了弱各向异性参数,为各向异性的研究奠定了基础;Ruger[2]对HTI介质的反射与透射特征进行了详细的研究,推导了纵波反射系数近似公式,并分析了方位各向异性与入射角和方位角的关系;Bakulin等[3]和Sayers等[4]利用带有线性滑动边界条件的有限平行裂缝和硬币型裂缝理论来估计各向异性参数与裂缝结构之间的关系,进而估测裂缝的发育方向和裂缝密度,利用各向异性差值判断裂缝的含油气性。AVAZ反演中弹性
参数和各向异性参数的反演具有很多的不确定性,不同的学者对此进行了相关的研究。Subhashis等[5]利用纵波和转换波随入射角和方位角的变化特征进行AVAZ反演，得到稳健的各向异性参数估计值;Downtown等[6-7]运用贝叶斯理论估计HTI介质AVAZ反演的弹性参数预测值的不确定性,提出了同步弹性参数反演,加入岩石物理模型的约束,克服了AVAZ反演的固有缺陷,从而定量地解释AVAZ反演的稳定性。
为了更好地识别裂缝储层特征,孔丽云等[8]深入分析了裂缝孔隙、裂缝弱度、裂缝渗透率等参数对双孔隙裂缝介质地震各向异性的影响。P波反射振幅与弹性参数之间的关系是裂缝性储层反演研究的主要问题之一,刘军迎等[9]认为裂缝储层的反射振幅随着观测方位的变化而变化,表现出明显的方位各向异性特征,可以应用该特征来预测碳酸盐岩储层中的裂缝密度及发育方向。张广智等[10]研究了裂缝弹性参数的各向异性AVO反演方法,准确地得到了反映储层特征的纵波阻抗、横波阻抗以及各向异性梯度。陈怀震等[11-12]研究了纵波动力学特征与裂缝流体因子、缝隙流体因子、裂缝密度、入射角和方位角等因素之间的关系特征,通过合成的模型数据对该方法的适用性进行了详细的分析。宗兆云等[13]分析了纵波模量、横波模量、杨氏模量和泊松比对储层反射特征的影响,提出了利用叠前AVO反演方法直接估算纵波模量、横波模量、杨氏模量和泊松比的算法,避免了间接求取带来的误差累积,为裂缝型储层的识别
提供了更加可靠的理论依据;宗兆云等[14]通过建立纵横波模量、密度以及各向异性差值与纵波反射系数之间的关系,进一步估测地层的地应力,这对页岩气储层的评价是非常有利的。裂缝储层的弹性阻抗也表现为强方位各向异性特征,通过方位各向异性弹性阻抗反演可以获取表征裂缝性质的弹性参数和各向异性参数[15-16]。
杨氏模量是表示岩石抗压缩能力的物理量,表征了储层的岩性特性;泊松比是横向应变与纵向应变的比值绝对值,与储层流体性质有关。为了全面地描述裂缝介质的储层特征、流体性质和各向异性特征,本文在上述研究的基础上,建立了裂缝介质的纵波反射系数与杨氏模量、泊松比和各向异性梯度之间的关系,利用模型数据对新的AVAZ反演算法进行测试,取得了较好的应用效果。
裂缝介质的各向异性特征可以由Thomsen提出的弱各向异性参数ε(v),δ(v),γ来表征,它们与弹性矩阵之间的关系为[1]：
Ruger研究了各向异性半空间界面的地震波反射特征,推导了纵波反射系数随入射角、方位角变化的关系式[2],如(2)式所示：
式中：θ和φ分别是纵波入射角和方位角ρ;vP和vS分别是纵、横波速度;vP0和vS0分别是上、下介质的平均纵波速度和横波速度。
按照扰动理论,分为背景场和扰动场,可以将纵波的AVO近似方程表示为：
式中(θ,φ)代表各向异性介质的反射系数；(θ)代表各向同性介质的反射系数；(θ,φ)代表各向异性介质反射系数的扰动量。且有：
式中为横、纵波速度比值和分别被称为纵、横波速度反射系数及密度反射系数。
宗兆云推导了纵波近似反射系数与杨氏模量(E),泊松比(σ)和密度反射系数之间的变化关系式[14]：
式中和分别被称为杨氏模量反射系数和泊松比反射系数。由于密度和纵波速度之间存在一定的关系,即Gardner公式：
式中:C是一个常数;g为纵波速度的幂指数。
对其求一阶偏导,求出密度反射系数与纵波速度反射系数的关系：
最后,可以求得密度反射系数与杨氏模量反射系数、泊松比反射系数的函数表示式：
将其代入到(6)式中,进一步得到：
Chen等[11]指出,在入射角小于30°的情况下,可以令sin2θtan2θ≈0,(5)式可以进一步简化为：
Ruger定义各向异性梯度Γ为：
最终,可以得到基于杨氏模量、泊松比和各向异性梯度的P波AVAZ近似方程：
在上述推导的基础上,选取两层的裂缝介质模型,相关参数如表1所示,对新推导的纵波AVAZ近似方程进行精度分析。由Ruger近似得到曲线和新推导公式的近似曲线对比结果(图1),对比发现,在入射角较大的情况下,新推导方程的近似精度与Ruger近似吻合程度较高,因此,利用该方程进行裂缝介质AVAZ反演可以提高裂缝预测的精度。
在不同方位角和入射角的情况下,可以将(13)式表示为矩阵的形式：
令：
cracks什么意思中文
(14)式可以表示为：
d=Gm
利用最小二乘算法求取出模型参数的最佳估算值：
2.1 实际工区测井资料合成叠前角度道集的AVAZ反演方法测试
为了验证反演算法的有效性,本文选取某实际工区的纵波速度、横波速度和各向异性参数等测井曲线(图2),利用新推导的公式求取不同入射角(0～30°)和方位角(0,30°,60°,90°)的P波AVO反射系数,并且与40Hz的雷克子波褶积得到不同方位角的叠前合成角度道集,对其添加信噪比为3的随机噪声。分别对不含噪声和信噪比为3的合成数据进行算法测试。图3是密度和纵波速度的幂指数拟合,拟合点基本分布在公式(7)所表示的曲线附近,说明密度和纵波速度之间满足关系式(7)。
图4是不含噪声时合成的不同方位角(0,30°,60°,90°)的叠前角度道集;图5是未加噪声合成数据经过基于杨氏模量、泊松比和各向异性梯度反演得到的反演结果与真实结果的对比;图6是信噪比为3时合成的不同方位角(0,30°,60°,90°)的叠前角度道集;图7是信噪比为3的合成数据经过基于杨氏模量、泊松比和各向异性梯度反演得到的反演结果与真实结果的对比。对比图7
与图5可以看出，噪声对杨氏模量的反演影响较小；不含噪声时泊松比和各向异性梯度的反演精度更高，与真实值吻合程度较高。
2.2 二维逆掩断层模型测试
再应用上述AVAZ反演方法对二维逆掩断层模型进行算法测试,进一步验证该算法的合理性。图8给出了二维逆掩断层的原始模型剖面,分别是杨氏模量剖面(图8a)、泊松比剖面(图8b)和各向异性梯度剖面(图8c)。通过正演模拟得到二维逆掩断层模型不同方位角(0,30°,60°,90°)的合成地震角度道集,并且加入信噪比为3的随机噪声。利用含噪声合成地震记录直接通过AVAZ反演获取二维逆掩断层模型的杨氏模量、泊松比和各向异性梯度,反演结果如图9所示。由图9可以看到，杨氏模量、泊松比和各向异性梯度均能准确反映该模型的构造特征,清晰地展示了地层界面,具有很强的横向连续性。测试结果进一步验证了新的P波AVAZ近似公式的正确性和反演算法的有效性。
本文推导了基于杨氏模量、泊松比和各向异性梯度的方位各向异性AVO方程,提出了一种AVAZ反演算法；选用某实际工区的测井数据和二维逆掩断层模型合成不同方位角的叠前角度道集测试了该AVAZ反演算法,为裂缝介质参数反演的发展提供了新的思路。通过以上研究
工作可以得出以下结论与认识。
1) 对新推导的方位各向异性AVO方程进行精度分析,通过与Ruger近似方程对比,可以看出新推导的方程能够满足裂缝介质AVAZ反演的精度要求。
2) 选择实际井资料模型和二维逆掩断层模型,合成不同方位角的叠前角度道集对AVAZ反演算法进行测试,结果表明,杨氏模量、泊松比和各向异性梯度的反演结果与真实值吻合程度较高,能够反演裂缝储层的构造特征和岩性特征等信息。
3) 值得指出的是，裂缝介质中流体的识别是裂缝介质研究的主要目的之一,因此，下一步需要对基于上述AVAZ反演结果的裂缝流体识别进行深入的研究。
【相关文献】
[1] Thomsen L.Weak elastic anisotropy[J].Geophysics,1986,51(10):1954-1966
[2] Ruger A.Reflection coefficients and azimuthal AVO analysis in anisotropic media[D].Colorado:Colorado School of Mines Golden,1996
[3] Bakulin A,Grechka V,Tsvankin I.Estimation of fracture parameters from reflection seismic data—part I:HTI model due to a single fracture set[J].Geophysics,2000,65(6):1788-1802
[4] Sayers C M,Simon D.Azimuth-dependent AVO in reservoirs containing non-orthogonal fracture sets[J].Geophysical Prospecting,2001,49(1):100-106

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。