叠前地震的粘弹性介质振幅随频率变化的反演方法和系统技术方案

公开了一种叠前地震的粘弹性介质振幅随频率变化的反演方法和系统。该方法可以包括以下步骤：基于先验信息建立反演模型；对地震资料进行预处理以获得振幅随频率变化的道集数据；以及基于反演模型，对道集数据进行振幅随频率变化反演以获得粘弹性介质的参数。

Method and system for inversion of viscoelastic medium amplitude with frequency variation in pre stack seismic

A method and system for inversion of viscoelastic media amplitude versus frequency variation in pre stack seismic. The method comprises the following steps: the inversion model is constructed based on the prior information; preprocessing of seismic data to obtain amplitude changes with the frequency of the gather data; and based on the inversion model, parameters of amplitude to gather data varies with frequency inversion to obtain viscoelastic medium.

【技术实现步骤摘要】

本专利技术涉及石油地球物理勘探领域，更具体地，涉及一种叠前地震的粘弹性介质振幅随频率变化(AVF)的反演方法和系统。
技术介绍
在石油地球物理勘探领域，叠前地震反演是油藏精细描述的重要手段。其中，地震波衰减通常被认为是直接与流体相关的参数，可直接用于储层预测。目前，常用的叠前地震反演主要有两种，一种是Connolly提出的弹性阻抗反演，另一种是基于Zoeppritz方程纵波反射系数近似式的叠前AVA三参数同时反演。弹性阻抗反演采用的是部分角度叠加后的数据，不可避免地会损失部分叠前信息，影响了反演的精度和有效性。AVA三参数同时反演有多种形式，包括利用Aki-Richard近似或Gidlow近似等，可反演纵波速度、横波速度和密度或者反演纵波阻抗，横波阻抗和密度。专利技术人发现，这些常规方法反演结果已很难满足隐蔽油气藏勘探的要求。尤其是作为叠前反演的重要参数的泊松比，其对地层岩性与流体的变化十分敏感，但上述方法都只能间接地得到泊松比，即先通过反演得到纵波速度(阻抗)与横波速度(阻抗)，再计算泊松比。由于反演的纵、横波速度或纵、横波阻抗本身存在误差，因此计算出的泊松比误差较大，有时甚至会得到没有地质意义的结果。因此，开发一种有效地反演地下介质的纵波衰减和横波衰减的方法和系统是非常有必要的。公开于本专利技术
技术介绍
部分的信息仅仅旨在加深对本专利技术的一般
技术介绍
的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
技术实现思路
本专利技术提供了一种叠前地震的粘弹性介质振幅随频率变化(AVF)的反演方法和系统，其能够有效地反演地下介质的纵波衰减和横波衰减，从而提高对储层含油气性预测的可靠性和准确性。根据本专利技术的一方面，提出了一种叠前地震的粘弹性介质振幅随频率变化的反演方法，该方法可以包括以下步骤：基于先验信息建立反演模型；对地震资料进行预处理以获得振幅随频率变化的道集数据；以及基于反演模型，对道集数据进行振幅随频率变化反演以获得粘弹性介质的参数。根据本专利技术的一方面，提出了一种叠前地震的粘弹性介质振幅随频率变化的反演系统，该系统可以包括以下单元：建模单元，其基于先验信息建立反演模型；预处理单元，其对地震资料进行预处理以获得振幅随频率变化的道集数据；以及反演单元，其基于反演模型，对道集数据进行振幅随频率变化反演以获得粘弹性介质的参数。本专利技术的方法和系统具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方案中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方案中进行详细陈述，这些附图和具体实施方案共同用于解释本专利技术的特定原理。附图说明通过结合附图对本专利技术的示例性实施方案进行更详细的描述，本专利技术的上述以及其它目的、特征和优势将变得更加明显，其中，在本专利技术的示例性实施方案中，相同的参考标号通常代表相同部件。图1是示出了根据本专利技术的示例性实施方案的粘弹性介质振幅随频率变化的反演方法的步骤的框图。图2a是示出了某工区实际测井曲线及解释中的纵波速度Vp的示意图。图2b是示出了某工区实际测井曲线及解释中的横波速度Vs的示意图。图2c是示出了某工区实际测井曲线及解释中的密度的示意图。图3是示出了根据本专利技术的示例性实施方案的频率域AVF道集的示意图。图4是示出了根据本专利技术的示例性实施方案的粘弹性介质振幅随频率变化的反演结果与某工区实际测井曲线的比较的示意图。具体实施方式下面将参照附图更详细地描述本专利技术的优选实施方案。虽然附图中显示了本专利技术的优选实施方案，然而应该理解，可以以各种形式实现本专利技术而不应被这里阐述的实施方案所限制。相反，提供这些实施方案是为了使本专利技术更加透彻和完整，并且能够将本专利技术的范围完整地传达给本领域的技术人员。实施方案1图1是示出了根据本专利技术的示例性实施方案的粘弹性介质振幅随频率变化的反演方法的步骤的框图。在该实施方案中，粘弹性介质振幅随频率变化的反演方法可以包括以下步骤：步骤101，基于先验信息建立反演模型；步骤102，对地震资料进行预处理以获得叠前地震振幅随频率变化的道集数据；以及步骤103，基于反演模型，对道集数据进行振幅随频率变化反演以获得粘弹性介质的参数。该实施方案通过建立反演模型，对地震资料进行预处理并进行振幅随频率变化反演等步骤，实现了有效地反演地下介质的纵波衰减和横波衰减的技术效果。建立反演模型在本实施方案的一个示例中，反演模型可以为弹性参数先验信息的统计学模型，可以通过概率密度函数的形式表示。首先可以获得反演模型所需的先验信息及似然函数。在一个示例中，先验信息可以是独立于地震数据的先验信息，并可以通过已知的渠道(比如，测井信息和地质知识等)获得，可以表示为P(m)。其中，P可以表示概率，m可以为地下介质的粘弹性参数，m可以表示为m＝[QP,QS]，QP,QS可以为粘弹性参数m的两个分量。在一个示例中，似然函数可以为在给定粘弹性参数m的情况下得到地震观测数据dobs的概率，可以表示为P(dobs|m)。其中，dobs可以为振幅随频率变 化的数据，也即频率域的AVF道集数据，dobs可以通过时频域分析方法获得。在一个示例中，基于贝叶斯理论的振幅随频率变化反演模型可表达为： P ( m | d o b s ) = P ( m ) P ( d o b s | m ) P ( d ) - - - ( 1 ) ]]>也即，在给定dobs的情况下的m的后验概率分布。P(m)为先验信息，d为地震数据，P(d)为地震数据概率表达。通常，可以假设粘弹性参数的分布满足高斯分布： m | d o b s ~ N n m ( μ m | d o b s 本文档来自技高网...

【技术保护点】
一种叠前地震的粘弹性介质振幅随频率变化的反演方法，所述方法包括以下步骤：基于先验信息建立反演模型；对地震资料进行预处理以获得振幅随频率变化的道集数据；以及基于所述反演模型，对所述道集数据进行振幅随频率变化反演以获得粘弹性介质的参数。

【技术特征摘要】
1.一种叠前地震的粘弹性介质振幅随频率变化的反演方法，所述方法包括以下步骤：基于先验信息建立反演模型；对地震资料进行预处理以获得振幅随频率变化的道集数据；以及基于所述反演模型，对所述道集数据进行振幅随频率变化反演以获得粘弹性介质的参数。2.根据权利要求1所述的叠前地震的粘弹性介质振幅随频率变化的反演方法，其中，所述反演模型为基于贝叶斯理论的振幅随频率变化反演模型。3.根据权利要求2所述的叠前地震的粘弹性介质振幅随频率变化的反演方法，其中，基于贝叶斯理论的振幅随频率变化反演模型表示为： P ( m | d obs ) = P ( m ) P ( d obs | m ) P ( d ) ]]>其中，m为粘弹性参数，dobs为振幅随频率变化的道集数据，P(dobs|m)为在给定的粘弹性参数m的情况下得到振幅随频率变化的道集数据dobs的概率，P(m)为先验信息，d为地震数据，P(d)为地震数据概率表达，P(m|dobs)为在给定的振幅随频率变化的道集数据dobs的情况下粘弹性参数m的后验概率分布。4.根据权利要求1所述的叠前地震的粘弹性介质振幅随频率变化的反演方法，其中，对地震资料进行预处理包括，基于粘弹性介质的Zeoppritz方程，根据粘弹性参数计算随入射角和频率变化的反射系数；以及将所述反射系数与子波褶积，以获得振幅随频率变化的道集数据。5.根据权利要求4所述的叠前地震的粘弹性介质振幅随频率变化的反演方法，其中，对地震资料进行预处理包括，将所述反射系数与30Hz的雷克子波褶积，从而获得振幅随频率变化道集数据。6.一种叠前地震的粘弹性...

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