一种方位频变流体因子识别储层的方法与流程

1.本发明属于油气地震勘探数据处理与解释领域，是一种通过新的与频率相关的方位avo公式实现新频散流体因子的反演，利用与储层敏感的新频散参数和宽方位地震数据实现高精度储层识别的解释技术。


背景技术：

2.频散是含流体储层的固有特性，可充分挖掘该属性进行储层识别，依赖频率的avo反演方法为从实际地震数据中提取与含流体储层有关的频散异常提供了一种有效技术。而如今依赖频率的avo反演技术是强行将频率引入常规avo公式进行频散反演，具有不合理性，需要利用直接与频率相关的avo公式进行反演，并且需要构建与储层敏感的流体因子来提高储层预测的精度。
3.不同频散参数对流体的敏感性存在差异，考虑到频散属性对流体的敏感性问题，基于敏感性分析的频散参数反演可获得与储层流体敏感的频散因子，进而提高储层预测精度。随着宽方位地震采集技术的发展，可充分利用宽方位叠前地震数据进行频散参数的反演进而进行储层的高精度预测。为此，本发明技术提供了一种通过构建对储层敏感的频散因子并结合宽方位地震数据进行储层识别的方法。


技术实现要素：

4.本发明是要提供一种方位频变流体因子识别储层的方法，通过推导的一种与频率直接相关的方位avo公式，实现五类频散因子的同步反演，并利用构建的与储层渗透率、流体粘度等参数有关的敏感频散流体因子直接识别储层。
5.本发明的目的在于提供一种方位频变流体因子识别储层的方法，所述方法主要包括以下步骤：
6.⑴
对叠前角道集数据进行时频变换，并构建与新频散流体因子相关的目标函数；
7.r＝ed
8.式中，r为m个频率下反射系数的列向量，共m
×
n行，e为m个频率下反射系数公式中呈m
×
n行、5列排列的系数矩阵，d
m
、d
μ
、和分别为纵波模量、横波模量、法向弱度、新流体因子和切向弱度的频散程度。m和μ分别表示纵波模量和横波模量，δ
n
和δ
t
分别为岩石的法向和切向弱度，ψ
n
为与裂缝弱度与松弛时间(受流体粘度、岩石渗透率、裂隙纵横比等参数的影响)的平方相关的参数，该参数即为构建的一种新的流体识别因子，表达式为：ψ
n
＝δ
n
γ2，
9.式中：
[0010][0011]
式中，η
f
为流体的粘度，k为岩石的渗透率，k
f
为流体的体积模量，a为裂隙纵横比，
γ为松弛时间。
[0012]
⑵
选取参考频率ω0，利用最小二乘求解该参考频率下五类模量对应的频散结果：p、s、n、k、l，并将该结果带入向量r中，利用下式进行五类频散因子的同步反演：
[0013]
d＝(e
t
e)
‑1e
t
r
[0014]
⑶
对研究工区的不同方位的叠前地震数据利用步骤
⑴
至
⑵
的计算，获取不同方位新频散流体因子的结果，进而利用其结果识别进行储层识别；
[0015]
⑷
重复上述步骤
⑴
至
⑶
，直到处理完整个工区的所有叠前地震角道集，即可获得整个工区的三维储层预测结果，该预测结果具有很高的精度，能够压制与储层无关的背景异常，对储层具有很强的敏感性，可对后续的储层精细评价、钻井部署等高精度勘探提供支撑。
[0016]
本发明的具体实现原理如下：
[0017]
基于hudson et al(1996)的等效介质模型，通过简化刚度矩阵参数，构建与裂缝弱度和新流体指示因子有关的新刚度矩阵参数，再利用散射函数与反射系数间的关系推导与频率相关的avoa反射系数公式。对于hti介质，其刚度矩阵c可表示为：
[0018][0019]
式中，m＝λ 2μ，e为裂隙密度，λ和μ为岩石的拉梅系数，和u
33
为与流体参数(流体粘度，体积模量)和裂隙性质(渗透率，裂隙纵横比)有关的重要参数。此处，我们考虑裂隙中充填有流体，则和u
33
的表达式为：
[0020][0021]
式中，ω＝2πf，g＝μ/m，η
f
为流体的粘度，k为岩石的渗透率，k
f
为流体的体积模量，a为裂隙纵横比，γ为松弛时间。
[0022]
经过分析表明，u
33
的虚部响度实部很小，可以忽略。进一步求取u
33
实部的近似表达式：
[0023][0024]
此外，研究表明，在实际地震应用中，的衰减和频散是可以忽略的，可以将填充
流体裂隙的简写为：
[0025][0026]
则可进一步得到刚度矩阵元素的表达式：
[0027]
c
11
＝m
‑
mδ
n
‑
ω2mψ
n
,
[0028]
c
12
＝λ
‑
λδ
n
‑
ω2λψ
n
,
[0029]
c
23
＝λ
‑
λ(1
‑
2g)δ
n
‑
ω2λ(1
‑
2g)ψ
n
,
[0030]
c
33
＝m
‑
λ(1
‑
2g)δ
n
‑
ω2λ(1
‑
2g)ψ
n
,
[0031]
c
44
＝μ,c
55
＝μ
‑
μδ
t
,
[0032]
c
13
＝c
21
＝c
31
＝c
12
,c
32
＝c
23
,c
66
＝c
55
,c
22
＝c
33
[0033]
式中，ψ
n
＝δ
n
γ2，ψ
n
为与裂缝弱度与松弛时间(受流体粘度、岩石渗透率、裂隙纵横比等参数的影响)的平方相关的参数，该参数可以作为一种新的流体识别因子，δ
n
＝4e/[3g(1
‑
g)]和δ
t
＝16e/[3(3
‑
2g)]分别为岩石的法向和切向弱度。
[0034]
可进一步得到刚度矩阵元素中扰动：
[0035]
δc
11
＝(m δm)
‑
(m δm)(δ
n
δδ
n
)
‑
ω2(m δm)(ψ
n
δψ
n
)
‑
(m
‑
mδ
n
‑
ω2mψ
n
)
[0036]
ꢀꢀꢀꢀ
≈δm
‑
δ
n
δm
‑
mδδ
n
‑
ω2mψ
n
‑
ω2ψ
n
δm
[0037]
δc
12
≈δλ
‑
δλδ
n
‑
λδδ
n
‑
ω2λδψ
n
‑
ω2ψ
n
δλ,
[0038]
δc
23
≈δλ
‑
(1
‑
2g)δ
n
δλ
‑
λ(1
‑
2g)δδ
n
‑
ω2λ(1
‑
2g)δψ
n
‑
ω2(1
‑
2g)ψ
n
δλ,
[0039]
δc
33
≈δm
‑
(1
‑
2g)δ
n
δλ
‑
λ(1
‑
2g)δδ
n
‑
ω2λ(1
‑
2g)δψ
n
‑
ω2(1
‑
2g)ψ
n
δλ,
[0040]
δc
44
＝δμ,
[0041]
δc
55
≈δμ
‑
δ
t
δμ
‑
μδδ
t
,
[0042]
δc
13
＝δc
21
＝δc
31
＝δc
12
,
[0043]
δc
32
＝δc
23
,δc
66
＝δc
55
,δc
22
＝δc
33
[0044]
式中，δm和δμ分别为反射界面两侧纵波模量和剪切模量的变化量，δδ
n
和δψ
n
为界面两侧δ
n
和新流体因子ψ
n
的变化量。此处假设界面两侧的扰动量较小，则可忽略δmδδ
n
、δmδψ
n
项。
[0045]
散射函数的一般表达式为：
[0046][0047]
式中，p
i
和p
s
分别为入射波和散射波的极化矢量，可表示为：
[0048][0049]
式中，θ为p波的入射角，为方位角，在弱散射建设条件下，pp波反射系数与散射函数之间有如下关系：
[0050][0051]
则有上述公式可得：
[0052][0053]
上式即为推导的新的纵波反射系数公式，它同时与入射角、方位角和频率有关，包含了与流体特性直接相关的流体因子ψ
n
。该公式可用于分析依赖频率的方位avo反演。
[0054]
再将上述推导的公式进行拓展，忽略密度项并将在参考频率ω0处进行泰勒展开，则有：
[0055][0056]
式中，
[0057]
a
′
(θ)＝a(θ),b
′
(θ)＝b(θ),
[0058][0059][0060][0061]
令：
[0062][0063][0064]
式中，d
m
、d
μ
、和分别为纵波模量、横波模量、法向弱度、新流体因子和切向弱度的频散程度。
[0065]
当方位角时，考虑到m个频率ω1,ω2,
…
,ω
m
的情况，定义m
×
n行的列向量r为：
[0066][0067]
定义e为m
×
n行、5列的矩阵：
[0068][0069]
则可得到目标函数：
[0070]
r＝ed
[0071]
式中，
[0072]
求解目标函数则需求取向量r中的δδ
n
(t,ω0)，δψ
n
(t,ω0)，δδ
t
(t,ω0)，分别令其为p、s、n、k、l。
[0073]
其中，p、s、n、k、l的求解过程为：
[0074]
a)将a(θ
i
)，b(θ
i
)，c(θ
i
)，d(θ
i
)和e(θ
i
)看成采样时间t和接收道n(每个采样点对应一个系数a，b、c、d和e)的函数，则：
[0075]
a(θ
i
)＝a(t,n),b(θ
i
)＝b(t,n),c(θ
i
)＝c(t,n),d(θ
i
)＝d(t,n),e(θ
i
)＝e(t,n)
[0076]
并对叠前avo道集进行谱均衡处理：
[0077][0078]
式中，w(f
i
,n)为谱均衡算子，为原始叠前角道集的时频谱，为谱均衡之后的结果。
[0079]
b)对于某一参考频率ω0，有如下表达形式：
[0080]
[0081]
令a
i
＝a(t,i)，b
i
＝b(t,i)，c
i
＝c(t,i)，d
i
＝d(t,i)，e
i
＝e(t,i)，n＝δδ
n
(t,f0)，k＝δψ
n
(t,f0)，l＝δδ
t
(t,f0)，则定义道集的拟合误差为：
[0082][0083]
为了寻找使上式的误差达到最小值的变量p，s，n，k和l，利用上式分别对变量p，s，n，k和l，求偏导数如下：
[0084][0085]
令上式等于0，则有：
[0086][0087]
求解上式即可得p、s、n、k、l。
[0088]
同样，给定不同的方位角，即可获得不同方位角时的p、s、n、k、l，进而可求得不同方位的频散因子。
[0089]
本发明的一种方位频变流体因子识别储层的方法，具有如下优越性：
[0090]
(1)推导的新的方位avo反射系数公式直接与频率相关，可直接用于依赖频率的
avo反演，避免了常规依赖频率avo反演中强行将频率引入avo公式的不合理性；
[0091]
(2)构建了直接与储层岩石渗透率、流体粘度有关的新流体因子，该因子对含流体储层具有更强的敏感性，能够用于高精度储层预测；
[0092]
(3)可同时实现五类频散因子的同步反演，得到的新频散因子对储层具有很好的识别能力，能够压制与储层无关的异常干扰，直接突出储层异常，并且分辨率高，横向连续性好，能够解决薄储层难以识别的问题；
[0093]
(4)可实现宽方位的频散因子反演，利用宽方位地震数据能够更好的分析裂缝型储层的特征，裂缝弱度等参数对识别致密裂缝型储层具有很好的效果。
附图说明
[0094]
图1是本发明技术利用的与频率和方位相关的反射系数对比图。
[0095]
图2是本发明技术反演得到的五类频散因子的单道结果对比图。
[0096]
图3是本发明实施例工区反演得到的新频散因子剖面。
具体实施方式
[0097]
本发明的具体实施方式如下：
[0098]
⑴
输入叠前角道集，并对角道集进行时频谱分析，生成不同频率下的振幅谱；
[0099]
⑵
对不同频率下的振幅谱进行振幅谱均衡处理，并获取地震数据的层速度；
[0100]
⑶
选定参考频率利用最小二乘反演该频率下模量的变化率，即求取p、s、n、k、l；
[0101]
⑷
给定方位角将求取的p、s、n、k、l带入目标函数，并进行反演求解，生成该方位角下的最终频散属性结果；
[0102]
⑸
重复步骤
⑴
至
⑷
即可获得所有方位时的频散因子结果；
[0103]
⑹
利用步骤
⑸
得到的新频散流体因子结果进行储层的识别；
[0104]
⑺
重复步骤
⑴
至
⑹
，处理研究工区所有的三维数据并对其进行储层预测并抽取沿层切片进行分析，即可得到整个工区不同目的层段的储层预测平面展布。
[0105]
为了更清楚的表达本发明的技术优势，结合附图度对本发明实施例做进一步的详细说明。本发明的实施例说明：
[0106]
图1为不同频率时的反射系数随入射角和方位角的变化对比。分别对频率为10hz、20hz、30hz、40hz、50hz情况下的反射系数进行对比，由图可知，反射系数随方位角的变化在180
°
位置处对称，方位角在0～180
°
的反射系数也在方位角为90
°
的位置对称，此外，可以明显看出不同频率成分的反射系数随方位角和入射角变化的差异较大，说明反射系数与频率有较强的敏感性，可利用此特性进行依赖频率的方位avo反演。
[0107]
图2为利用本发明技术反演得到的五类频散因子单道结果对比图。由图可知，图a
‑
d中的四个频散因子均受明显的无关异常干扰(黑色箭头所示)，难以区分有利储层，而图e中新构建的频散因子最大限度地压制了其余层的异常干扰，仅在两个储层位置处有很强的频散异常。
[0108]
图3为实施例工区1反演得到的新频散因子剖面。由测井曲线可知，在此处具有两套储层，而由预测的剖面结果可知，新的频散因子能够很好地识别出这两套储层，具有较高的分辨率，并且横向连续性很好，压制了无关的背景异常干扰。
[0109]
上述实施例仅用于说明本发明，其中方法的各实施步骤等都是可以有所变化的，凡是在本发明技术方案的基础上进行的等同变换和改进，均不应排除在本发明的保护范围之外。