一种针对岩石力学参数的叠前反演方法及装置制造方法及图纸

本发明专利技术公开一种针对岩石力学参数的叠前反演方法及装置。该方法包括：建立岩石力学参数与合成数据的关系模型；基于岩石力学参数和合成数据都满足高斯分布，利用最小二乘法分别得到岩石力学参数的概率密度函数和合成数据的概率密度函数；根据贝叶斯理论，得到岩石力学参数的后验概率密度函数；基于所述岩石力学参数的后验概率密度函数满足高斯分布求取岩石力学参数的平均值和协方差，由此得到地下介质中岩石力学参数的分布。本发明专利技术利用数学反演理论为指导，实现了从叠前地震数据来对地下岩石力学参数的高精度反演。

【技术实现步骤摘要】
一种针对岩石力学参数的叠前反演方法及装置
本专利技术涉及油气地球物理勘探领域，更具体地，涉及一种针对岩石力学参数的叠前反演方法和一种针对岩石力学参数的叠前反演装置。
技术介绍
叠前反演技术在油气勘探开发中起到了非常重要的作用，它是储层弹性参数与地表地震响应之间的桥梁。基于平面波假设，在纵波入射情况下纵横波反射系数可以通过Zoeppritz方程表示，这就是所有叠前反演的基础。这种方法最早由地球物理学家提出，因此几乎所有的叠前反演方法都是针对地下介质中的速度场和密度场信息。专利技术人发现，对于非常规油气藏，例如致密型储层或页岩气储层，其岩石力学性质对开发方案的制定有着重要影响。杨氏模量和泊松比作为表征岩石力学性质的两个重要参数，需着重分析。目前，尚未有成熟的、完整的针对岩石力学参数进行反演的叠前反演方法。因此，有必要开发一种针对岩石力学参数的叠前反演方法和装置。公开于本专利技术
技术介绍
部分的信息仅仅旨在加深对本专利技术的一般
技术介绍
的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
技术实现思路
本专利技术的目的是通过叠前地震数据反演得到岩石力学参数，包本文档来自技高网...

【技术保护点】
一种针对岩石力学参数的叠前反演方法，其特征在于，所述方法包括以下步骤：建立岩石力学参数m与合成数据d的关系模型；基于岩石力学参数m和合成数据d都满足高斯分布，利用最小二乘法分别得到岩石力学参数m的概率密度函数ρ(m)和合成数据d的概率密度函数ρ(d)；根据贝叶斯理论，基于概率密度函数ρ(m)和ρ(d)得到岩石力学参数m的后验概率密度函数p(m|d)；基于所述岩石力学参数m的后验概率密度函数p(m|d)满足高斯分布求取岩石力学参数m的平均值

【技术特征摘要】
1.一种针对岩石力学参数的叠前反演方法，其特征在于，所述方法包括以下步骤：建立岩石力学参数m与合成数据d的关系模型；基于岩石力学参数m和合成数据d都满足高斯分布，利用最小二乘法分别得到岩石力学参数m的概率密度函数ρ(m)和合成数据d的概率密度函数ρ(d)；根据贝叶斯理论，基于概率密度函数ρ(m)和ρ(d)得到岩石力学参数m的后验概率密度函数p(m|d)；基于所述岩石力学参数m的后验概率密度函数p(m|d)满足高斯分布求取岩石力学参数m的平均值和协方差Cm，由此得到地下介质中岩石力学参数的分布。2.根据权利要求1所述的针对岩石力学参数的叠前反演方法，其特征在于，岩石力学参数m与合成数据d的关系模型为：d＝Gm+e其中，G为褶积模型算子，e代表噪声。3.根据权利要求1所述的针对岩石力学参数的叠前反演方法，其特征在于，岩石力学参数m和合成数据d的概率密度函数分别为：其中，mprior是岩石力学参数m的先验模型的期望值，Cm是岩石力学参数m的协方差，dobs是所述合成数据d的实际数据，Cd是实际数据的协方差，L为高斯分布函数的系数。4.根据权利要求3所述的针对岩石力学参数的叠前反演方法，其特征在于，岩石力学参数m的后验概率密度函数为：p(m|d)＝L·exp(-S(m))其中，S(m)的表达式为：其中，G(m)为褶积算子G作用在模型参数m上的结果。5.根据权利要求4所述的针对岩石力学参数的叠前反演方法，其特征在于，基于岩石力学参数m的后验概率密度函数p(m|d)满足高斯分布，将p(m|d)表达为：其中，为岩石力学参数m的平均值，Cm为岩石力学参数m的协方差。6.一种针对岩石...

【专利技术属性】
技术研发人员：滕龙，胡华锋，周单，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京,11