一种基于双相介质固液解耦的叠前地震反演方法及系统技术方案

本发明专利技术提供一种基于双相介质固液解耦的叠前地震反演方法及系统，所述方法包括：采集目的储层的叠前地震资料、测井资料以及地质成果资料；构建固液解耦近似方程；根据所述的叠前地震资料、测井资料、地质成果资料以及所述的固液解耦近似方程提取子波；根据所述的叠前地震资料、测井资料、地质成果资料以及子波确定绝对弹性阻抗数据体；根据所述的绝对弹性阻抗数据体确定流体体积模量数据体；根据所述的流体体积模量数据体识别当前储层孔隙中流体的类型。得到的流体体积模量作为流体因子不仅具有更高的流体指示灵敏性，而且有效的提高了储层流体识别和储层特征描述的可靠性。

【技术实现步骤摘要】

本专利技术关于石油地球物理勘探领域，特别是关于叠前地震的反演技术，具体的讲是一种基于双相介质固液解耦的叠前地震反演方法及系统。
技术介绍
在石油工业中，利用流体因子进行地下储层流体识别是业界常用的烃类描述与预测方式。随着全球地球物理勘探的重点从构造油气藏向岩性油气藏转变，早先在流体识别领域应用广泛的“亮点”及AVO属性分析技术已无法满足现阶段复杂岩性油气藏的储层流体识别需求。针对上述问题，现有技术中许多学者提出了不同类型的流体因子，主要包括两类：(1)通过岩石物理参数的权差运算来突出体现流体异常效应，如交会图角、流体因子角、泊松阻抗等，但是此类流体因子的岩石物理意义模糊且普适性较差；(2)将对岩石孔隙流体敏感的弹性参数直接作为流体因子，如拉梅参数，泊松比，Russell流体因子以及Gassmann流体项f等，该类流体因子岩石物理意义明确且普适性较好，在现阶段储层流体识别中占据主要地位，但是在岩石固结程度较为复杂的地区应用该类流体因子时，受岩石固液两相组分的耦合效应影响，特别是受孔隙度影响，在实际储层的孔隙流体判识中会出现流体识别假象。另外，常规流体因子的计算多依赖弹性参数之间的间接代数运算，受地震反问题固有的“病态”问题影响，最终计算的流体因子无法避免会引入累计误差，降低了储层流体指示作用的可靠性。因此，随着以复杂岩性因素为主导作用的岩性地层油气藏勘探程度的加深，急需研究一种能够准确反应孔隙流体类型的高灵敏流体因子类型以及更为可靠的流体因子直接反演方法。
技术实现思路
为了克服现有技术存在的上述问题，本专利技术提供了一种基于双相介质固液解耦的叠前地震反演方法及系统，以岩石固液双相效应解耦并反演完全体现了流体效应的流体体积模量参数Kf为出发点，利用叠前弹性阻抗反演实现了流体因子的直接提取，得到的流体体积模量作为流体因子不仅具有更高的流体指示灵敏性，而且有效的提高了储层流体识别和储层特征描述的可靠性。本专利技术的目的之一是，提供一种基于双相介质固液解耦的叠前地震反演方法，包括：采集目的储层的叠前地震资料、测井资料以及地质成果资料；构建固液解耦近似方程；根据所述的叠前地震资料、测井资料、地质成果资料以及所述的固液解耦近似方程提取子波；根据所述的叠前地震资料、测井资料、地质成果资料以及子波确定绝对弹性阻抗数据体；根据所述的绝对弹性阻抗数据体确定流体体积模量数据体；根据所述的流体体积模量数据体识别当前储层孔隙中流体的类型。本专利技术的目的之一是，提供了一种基于双相介质固液解耦的叠前地震反演系统，包括：资料采集装置，用于采集目的储层的叠前地震资料、测井资料以及地质成果资料；固液解耦近似方程构建装置，用于构建固液解耦近似方程；子波提取装置，用于根据所述的叠前地震资料、测井资料、地质成果资料以及所述的固液解耦近似方程提取子波；弹性阻抗数据体确定装置，用于根据所述的叠前地震资料、测井资料、地质成果资料以及子波确定绝对弹性阻抗数据体；流体体积模量确定装置，用于根据所述的绝对弹性阻抗数据体确定流体体积模量数据体；流体类型识别装置，用于根据所述的流体体积模量数据体识别当前储层孔隙中流体的类型。本专利技术的有益效果在于，提供了一种基于双相介质固液解耦的叠前地震反演方法及系统，以岩石固液双相效应解耦并反演完全体现了流体效应的流体体积模量参数Kf为出发点，利用叠前弹性阻抗反演实现了流体因子的直接提取，得到的流体体积模量作为流体因子不仅具有更高的流体指示灵敏性，而且有效的提高了储层流体识别和储层特征描述的可靠性，常规的流体因子岩石物理意义模糊且普适性较差，而且受岩石双相特性的耦合效应影响，在表征岩石流体效应的时候易受岩石固体效应(特别是孔隙度参数)影响，造成流体识别假象，而流体体积模量Kf只与岩石孔隙流体有关，不受其他岩石固体效应影响，作为一项对储层孔隙流体类型敏感的弹性参数在流体识别领域的效果更佳，利用叠前弹性阻抗反演实现了流体因子的直接提取，可以避免引入累计误差，提高了储层流体指示作用的可靠性。为让本专利技术的上述和其他目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附图式，作详细说明如下。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1是利用某工区实际数据计算得到的常用流体因子的流体指示系数比较图；图2是本专利技术提供的基于固液解耦近似方程的流体因子叠前地震反演方法在具体实施例中的流程图；图3是根据某工区实际地震数据进行基于固液解耦近似方程的流体因子叠前地震反演后得到的流体因子Kf的剖面示意图；图4为本专利技术实施例提供的一种基于双相介质固液解耦的叠前地震反演方法的实施方式一的流程图；图5为本专利技术实施例提供的一种基于双相介质固液解耦的叠前地震反演方法的实施方式二的流程图；图6为本专利技术实施例提供的一种基于双相介质固液解耦的叠前地震反演方法的实施方式三的流程图；图7为本专利技术实施例提供的一种基于双相介质固液解耦的叠前地震反演方法的实施方式四的流程图；图8为图4中的步骤S103的具体流程图；图9为图4中的步骤S104的具体流程图；图10为本专利技术实施例提供的一种基于双相介质固液解耦的叠前地震反演系统的实施方式一的结构框图；图11为本专利技术实施例提供的一种基于双相介质固液解耦的叠前地震反演系统的实施方式二的结构框图；图12为本专利技术实施例提供的一种基于双相介质固液解耦的叠前地震反演系统的实施方式三的结构框图；图13为本专利技术实施例提供的一种基于双相介质固液解耦的叠前地震反演系统的实施方式四的结构框图；图14为本专利技术实施例提供的一种基于双相介质固液解耦的叠前地震反演系统中的子波提取装置300的具体结构框图；图15为本专利技术实施例提供的一种基于双相介质固液解耦的叠前地震反演系统中的弹性阻抗数据体确定装置400的具体结构框图。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。在现阶段，储层流体识别中占据主要地位的流本文档来自技高网...

【技术保护点】
一种基于双相介质固液解耦的叠前地震反演方法，其特征是，所述的方法具体包括：采集目的储层的叠前地震资料、测井资料以及地质成果资料；构建固液解耦近似方程；根据所述的叠前地震资料、测井资料、地质成果资料以及所述的固液解耦近似方程提取子波；根据所述的叠前地震资料、测井资料、地质成果资料以及子波确定绝对弹性阻抗数据体；根据所述的绝对弹性阻抗数据体确定流体体积模量数据体；根据所述的流体体积模量数据体识别当前储层孔隙中流体的类型。

【技术特征摘要】
1.一种基于双相介质固液解耦的叠前地震反演方法，其特征是，所述的方法具体包括：
采集目的储层的叠前地震资料、测井资料以及地质成果资料；
构建固液解耦近似方程；
根据所述的叠前地震资料、测井资料、地质成果资料以及所述的固液解耦近似方程提
取子波；
根据所述的叠前地震资料、测井资料、地质成果资料以及子波确定绝对弹性阻抗数据
体；
根据所述的绝对弹性阻抗数据体确定流体体积模量数据体；
根据所述的流体体积模量数据体识别当前储层孔隙中流体的类型。
2.根据权利要求1所述的方法，其特征是，所述的方法在构建固液解耦近似方程之
前还包括：
对所述的叠前地震资料进行保幅处理；
对经过保幅处理的叠前地震资料进行去噪处理；
根据目的储层对经过去噪处理的叠前地震资料进行角度叠加，生成四个入射角度的角
度叠加道集。
3.根据权利要求2所述的方法，其特征是，所述的方法在构建固液解耦近似方程之
前还包括：
对所述的测井资料进行环境校正；
对经过环境校正后的测井资料进行奇异值消除；
根据经过奇异值消除后的测井资料确定所述目的储层的干岩石纵横波速度比、饱和岩
石纵横波速度比；
根据所述的目的储层的干岩样纵横波速度比以及所述的测井资料确定与所述测井资
料对应的流体体积模量曲线以及剪切模量曲线。
4.根据权利要求3所述的方法，其特征是，所述的方法在构建固液解耦近似方程之
前还包括：
根据所述的叠前地震资料、测井资料、地质成果资料对目的储层进行精细层位追踪，
得到地震层位数据。
5.根据权利要求1所述的方法，其特征是，所述的固液解耦近似方程为：
RPP(θ)=[(1-γdry2γsat2)sec2θ4]ΔKfKf+[γdry24γsat2sec2θ-2γsat2sin2θ]Δ(fm)fm[12-sec2θ4]Δρρ+(sec2θ4-γdry22γsat2sec2θ+2γsat2sin2θ)Δφφ]]>其中，fm＝φμ为地下岩石的固体项，θ为入射角度，Kf为地下岩石的流体体积模量，
fm为地下岩石的固体项，φ为地下岩石的孔隙度，ρ为地下岩石的密度，μ为地下岩石的
剪切模量，△Kf为界面两侧的流体体积模量的差值，△fm为界面两侧固体项的差值，△ρ为
界面两侧的密度的差值，△φ为界面两侧的孔隙度的差值，为
干岩石纵横波速度比的平方，为饱和岩石纵横波速度比的平方。
6.根据权利要求4所述的方法，其特征是，根根据所述的叠前地震资料、测井资料、
地质成果资料以及固液解耦近似方程提取子波具体包括：
从所述的测井资料中提取出地下岩石的密度以及孔隙度；
根据所述的剪切模量曲线以及所述的孔隙度确定地下岩石的固体项；
基于所述的固液解耦近似方程推导出弹性阻抗方程；
根据所述的入射角度、地下岩石的密度、固体项、孔隙度、流体体积模量曲线以及弹
性阻抗方程确定弹性阻抗伪测井曲线；
针对不同的入射角度，根据入射角度对应的角度叠加道集、弹性阻抗伪测井曲线、地
震层位数据提取不同的子波。
7.根据权利要求6所述的方法，其特征是，所述的弹性阻抗方程为：
EI(θ)=EI0(KfKf0)a(fmfm0)b(φφ0)c(ρρ0)d]]>其中，a=(1-γdry2γsat2)sec2θ2,b=γdry22γsat2sec2θ-4γsat2sin2θ;]]>c=sec2θ2-γdry2γsat2sec2θ+4γsat2sin2θ,d=1-sec2θ2;]]>Kf0、fm0、ρ0和φ0分别为Kf、fm、ρ和φ的平
均值。
8.根据权利要求6所述的方法，其特征是，根据所述的叠前地震资料、测井资料、

\t地质成果资料以及子波确定绝对弹性阻抗数据体具体包括：
根据所述的叠前地震资料、测井资料、地震层位数据以及子波确定出相对弹性阻抗数
据体；
根据所述的弹性阻抗伪测井曲线以及地震层位数据构建低频阻抗数据体；
根据所述的低频阻抗数据体对所述的弹性阻抗数据体进行补偿，得到绝对弹性阻抗数
据体。
9.一种基于双相介质固液解耦的叠前地震反演系统，其特征是，所述的系统具体包括：
资料采集装置，用于采集目的储层的叠前地震资料、测井资料以及地质成果资料；<...

【专利技术属性】
技术研发人员：印兴耀，张世鑫，李超，曹丹平，吴国忱，宗兆云，
申请(专利权)人：中国石油天然气股份有限公司，中国石油大学华东，
类型：发明
国别省市：北京;11