【技术实现步骤摘要】
本专利技术属于油气地球物理勘探，尤其涉及一种关于致密砂岩含气量识别因子与孔隙度的avo反演方法。
技术介绍
1、致密砂岩气储层具有低孔低渗、孔隙结构复杂和非均质性强等特点，这些特点导致了甜点识别困难。流体识别是评价致密砂岩气储层潜在生产能力和天然气资源的关键。smith和gidlow首先定义了由纵波和横波速度的相对变化加权组成的流体因子。fatti等定义了新的流体因子。之后，流体因子的定义扩展到阻抗域。goodway等分析了弹性模量比常规的速度和密度具有更明显的物理意义，并提出了拉梅阻抗来区分储层流体。虽然基于弹性参数提出了许多新的流体因子，但对于致密储层，由于其复杂的孔隙结构和强非均质性，应根据流体特征建立新的流体因子，以避免流体识别误差。russell将流体项和骨架项从各向同性地震速度公式中分离出来，并使用流体项(ρf)作为敏感的孔隙流体指示因子。russell等提出了流体因子、剪切模量和密度的三项avo(amplitude variation withoffset，振幅随偏移距的变化)表达式。然而，gassmann流体项f可能导致孔隙流体识别模糊，因为其敏感性受到孔隙流体和孔隙度的叠加影响。为了解决这一问题，提出了一个新的包含流体体积模量的线性化avo近似式(yin和zhang)。根据地层流体对弹性参数的敏感性，xu等根据不同的地质条件提出了致密气敏感参数λ/vs和扩展属性f/vs。
2、气体含量是影响致密砂岩气储层有效开发的因素之一。根据地震数据估算与流体和孔隙度相关的储层性质对于预测高质量储层至关重要，这在
技术实现思路
1、本专利技术实施例的目的在于提供一种关于致密砂岩含气量识别因子与孔隙度的avo反演方法，旨在解决上述
技术介绍
中提出的问题。
2、本专利技术实施例是这样实现的，一种关于致密砂岩含气量识别因子与孔隙度的avo反演方法，包括以下步骤：
3、步骤1、通过测井数据分析，提出了一个新的参数f＝1/(kf×μ)，因其与sg×φ具有明显的正相关性，因此将其作为含气量识别因子；
4、步骤2、通过岩石物理理论与avo理论，构建包含含气量指示因子f、横波阻抗和孔隙度乘积is×φ、密度ρ和孔隙度φ的avo表达式；
5、步骤3、与zoeppritz反射系数公式和fatti反射系数公式对比，对构建的avo表达式进行精度分析；
6、步骤4：获取目标储层的叠前地震资料，生成叠前角道集；
7、步骤5：基于步骤2的avo表达式，利用弹性阻抗对储层参数f进行直接反演，用于致密砂岩储层甜点区域的直接检测。
8、步骤6：应用测井数据标定反演结果，验证构建的f对致密砂岩储层高含气量区域预测识别的精度与可靠性。
9、进一步的技术方案，在所述步骤2中，murphy等(1991)提出了gassmann公式的速度形式：
10、
11、
12、式中，kdry为干岩的体模量，μ为岩石骨架的剪切模量，f为流体项，vp、vs和ρ分别为饱和岩石的纵波速度、横波速度和密度。
13、此外，russell(2011)给出关系：
14、
15、
16、公式(3)两侧同时乘μ后，将公式(2)代入，得到：
17、
18、将公式(5)带入公式(1)，得到：
19、
20、公式(6)两侧同时除ρvp2后，将公式(4)代入，得到：
21、
22、利用关系ip＝ρvp，公式(7)可以进一步写为：
23、
24、公式(8)两侧同时乘μ，并使用han和batzle(2003)给出的流体项f的简化形式：
25、f＝kfg(φ) (9)
26、则公式(8)可以表示为：
27、
28、公式(10)两侧同时乘并定义参数f：
29、
30、得到：
31、
32、利用表达式：
33、
34、以及关系公式(12)可以进一步写为：
35、
36、如果定义：
37、
38、
39、通过将多变量微积分规则应用于公式(15)和(16)中，可以得到跨地下界面对应值的差值：
40、
41、
42、将公式(15)代入公式(17)，得到：
43、
44、进一步地，结合公式(15)和关系ip＝isγsat，可以得到：
45、
46、此外，将公式(16)代入公式(18)，得到：
47、
48、进一步地，δb/b可以表示为：
49、
50、根据yin和zhang(2014)的研究结果：
51、
52、公式(22)可以表示为：
53、
54、当a＝b≠0,可以假设δa＝δb，并进一步得到δa/a＝δb/b。因此，结合公式(20)和(24)得到：
55、
56、公式(25)可以进一步写为：
57、
58、将公式(26)代入fatti的反射系数方程(fatti等，1994)：
59、
60、得到：
61、
62、公式(28)可以进一步写为：
63、
64、进一步的技术方案，在所述步骤5中，根据connolly(1999)提出的弹性阻抗反转理论，p波反射系数rpp(θ)与弹性阻抗ei(θ)之间的关系可以用与声阻抗相似的形式表示：
65、
66、其中，ei1(θ)和ei2(θ)为上下介质在反射界面上的弹性阻抗，而δei(θ)和为两者的差值和平均值。
67、将公式(29)代入为公式(30)，两边同时积分和指数：
68、ei(θ)＝fa(θ)(isφ)b(θ)ρc(θ)φd(θ) (31)
69、其中：
70、
71、
72、
73、
74、为了使弹性阻抗数据体的尺寸与声阻抗相同，根据whitcombe(2002)提出的弹性阻抗归一化方法，对(31)进行归一化：
75、
76、其中，f0、is0、φ0和ρ0为平均值；ei0为弹性阻抗归一化因子，可表示为：
77、
78、公式(36)两边除以ei0后取自然对数
79、
80、得到方程组：
81、
82、本专利技术实施例提供的一种关于致密砂岩含气量识别因子与孔隙度的avo反演方法，该方法提出了一个因子f＝1/(kf×本文档来自技高网...
【技术保护点】
1.一种关于致密砂岩含气量识别因子与孔隙度的AVO反演方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的关于致密砂岩含气量识别因子与孔隙度的AVO反演方法，其特征在于，所述步骤2中构建的AVO表达式为：
3.根据权利要求2所述的关于致密砂岩含气量识别因子与孔隙度的AVO反演方法，其特征在于，在所述步骤5中，利用弹性阻抗对储层参数F进行直接反演的公式为：
【技术特征摘要】
1.一种关于致密砂岩含气量识别因子与孔隙度的avo反演方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的关于致密砂岩含气量识别因子与孔隙度的avo反演方法，其特征在于，所述步...
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。