岩相和弹性参数同步反演方法、装置、电子设备及介质制造方法及图纸

本发明专利技术提出了一种岩相和弹性参数同步反演方法、装置、电子设备及介质。反演方法包括：初始化待反演工区的岩相；利用岩相约束的贝叶斯线性反演计算弹性参数，并计算弹性参数关于岩相的似然函数；利用修正的Viterbi算法计算岩相关于弹性参数的后验概率最大值，并更新岩相，反复迭代直至收敛。本发明专利技术基于测井数据建立不同岩相的弹性低频趋势，可有效避免后续岩相预测的弹性压实“陷阱”。通过利用岩相反演结果约束弹性参数反演结果的低频和弹性参数之间的相关性，提高弹性参数反演结果的精度，进而更新弹性参数关于岩相的似然函数，并在贝叶斯网络和马尔科夫随机场先验约束下，更新岩相反演结果，提高岩相反演结果的稳定性和空间连续性。续性。续性。

【技术实现步骤摘要】
岩相和弹性参数同步反演方法、装置、电子设备及介质


[0001]本专利技术属于油气探勘领域，涉及地震资料解释与反演技术，更具体地涉及一种基于马尔科夫模型的岩相和弹性参数同步反演方法、反演装置、电子设备及存储介质。

技术介绍

[0002]基于地震数据的岩相预测一直是油气地球物理在勘探开发应用中的重点和难点，通过结合测井数据的岩相分类和油气解释结论，可以实现基于地震数据的岩性、物性以及含油气性的定量分类，因此也可以更加直观的给出勘探开发的有利“甜点区域”。
[0003]现阶段基于地震数据的岩相预测大致分为两类：一类是基于监督或者无监督学习的属性融合方法，该方法主要是依靠地震数据或地震属性的空间结构的相似性，并借助机器学习对其挖掘，实现岩相分类。这种分类方法并不需要借助明确的地球物理意义，与其说是岩相预测，不如说是大套的地震相预测，因此其垂向分辨率可能较低。例如，波形聚类、PCA属性融合等方法属于此类的岩相预测。第二类是基于叠前弹性参数反演结果的基础上，并结合测井上获取的岩石物理关系，利用统计的方法给出岩相关于弹性参数的似然关系，进而对岩相进行本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于马尔科夫模型的岩相和弹性参数同步反演方法，其特征在于，包括：初始化待反演工区的岩相；利用岩相约束的贝叶斯线性反演计算弹性参数，并计算弹性参数关于岩相的似然函数；利用修正的Viterbi算法计算岩相关于弹性参数的后验概率最大值，并更新岩相，反复迭代直至收敛。2.根据权利要求1所述的基于马尔科夫模型的岩相和弹性参数同步反演方法，其特征在于，所述初始化待反演工区的岩相包括：结合测井数据中的岩相曲线，统计岩相的先验分布和岩相概率转移矩阵P，将岩相先验分布中概率最大的岩相赋值到工区中每一道的每一个采样点中，作为初始化岩相数据体。3.根据权利要求1所述的基于马尔科夫模型的岩相和弹性参数同步反演方法，其特征在于，利用岩相约束的贝叶斯线性反演计算弹性参数包括：结合测井数据中的弹性参数，利用指数函数从测井数据中拟合所有岩相的纵波阻抗低频趋势
k
I
T
、纵横波速度比低频趋势
k
γ
T
、密度低频趋势
k
ρ
T
，并计算所有岩相的弹性参数协方差矩阵
k
Σ
m
；输入M个角度叠加地震数据和M个角度子波，并构建角度叠加地震数据向量d和角度子波褶积矩阵；利用岩相约束的贝叶斯线性反演计算本次迭代的弹性参数反演结果。4.根据权利要求3所述的基于马尔科夫模型的岩相和弹性参数同步反演方法，其特征在于，岩相约束的贝叶斯线性反演获取的弹性参数反演结果为：
π
μ＝
π
μ
T
+(G
·
π
Σ
m
)
T
(G
·
π
Σ
m
·
G
T
+σ
n2
I)-1
(d-G
·
π
μ
T
)其中，
π
I为纵波阻抗自然对数反演结果，
π
γ为纵横波速度比自然对数反演结果，
π
ρ为密度自然对数反演结果；
π
I
T
、
π
γ
T
和
π
ρ
T
分别表示利用前一次迭代更新的岩相π确定的纵波阻抗自然对数低频模型、纵横波速度比自然对数低频模型和密度自然对数低频模型；
π
Σ
m
表示利用前一次迭代更新的岩相π确定的参数协方差矩阵；G表示由子波褶积矩阵、AVO近似公式系数以及差分矩阵确定的AVO正演矩阵；表示噪声的方...

【专利技术属性】
技术研发人员：张丰麒，张宏，莫延钢，
申请(专利权)人：中国石油化工股份有限公司石油勘探开发研究院，
类型：发明
国别省市：