CN108572389A - 频变粘弹性流体因子叠前地震反演方法 - Google Patents

Abstract

本发明提供一种频变粘弹性流体因子叠前地震反演方法，该频变粘弹性流体因子叠前地震反演方法包括：步骤1，构建粘弹性介质反射特征方程并分析该方程的合理性；步骤2，进行地震资料多尺度分解；步骤3，针对每个尺度或频率内的地震反射记录，开展频变弹性阻抗反演方法；步骤4，确定频变加权系数，并在此基础上提取频变粘弹性流体因子；步骤5，基于频变粘弹性流体因子，对流体的分布状态进行预测。该频变粘弹性流体因子叠前地震反演方法，结合研究目标岩石物理分析，实现基于频变粘弹性流体因子叠前地震反演的储层流体识别，降低了流体检测的多解性，消除了仅仅考虑振幅进行流体检测时出现的识别假象。

Description

频变粘弹性流体因子叠前地震反演方法

技术领域

本发明涉及勘探地球物理地震流体识别领域，特别是涉及到一种频变粘弹性流体因子叠前地震反演方法。

背景技术

储层流体识别是油气藏勘探与储层评价的重要环节。在人类油气勘探活动已经经历了两个多世纪的今天，我们面临的是日益复杂的地下条件，隐蔽性强的剩余油气资源分布，识别和描述难度较大的岩性地层油气藏以及高额的勘探开发成本。与此同时，随着地震采集、处理以及解释技术的不断进步与发展，以地质知识和油气富集规律为指导，基于地震资料蕴含的振幅、频率等信息实现储层含流体识别已成为可能。以地震资料为主体进行储层流体识别可以在一定程度上增加勘探成功率，提高生产效率以及减少勘探开发成本，是现阶段油气勘探领域的研究热点之一。

为了降低勘探开发风险，提高勘探准确率成为地球物理工作者近年来着重研究的问题，因此，基于地震资料的储层流体识别技术成为研究的热点。Smith和Gidlow(1987)提出可以利用叠前数据通过不同加权函数进行叠加得到流体因子和伪泊松比剖面来预测岩性和流体，并首次提出了流体因子这个概念，该技术促进了用于流体异常识别的AVO截距-梯度交会技术的发展(Verm和Hilterman，1995；Castagna，1998)。Goodway等人(1997)提出了lambda-mu-rho技术，利用拉梅弹性参数进行储层流体识别。Hilterman(2001)在Goodway和Hedlin等人研究成果的基础上，对基于AVO的流体识别技术进行了总结。Batzle(2001)对流体指示因子进行了比较，提出对碎屑岩来说拉梅参数组合属性对流体类型最为敏感，且特别强调在实际应用中要根据区域特性对流体因子进行敏感性选择。George(2003)根据叠前AVO分析，提出了流体因子角度和交会图角度的概念，通过模型试算和实际应用发现这两种属性对储层流体类型有较强的识别能力。宁忠华等人(2006)等人在总结分析前人方法的基础上，提出了高灵敏度流体因子的概念。Mark等人(2006)提出了泊松阻抗的概念。在考虑多孔饱和弹性介质的前提下，Russell等人(2003，2006)总结了前人的观点，利用Biot-Gassmann方程对饱和流体条件下的的纵波速度方程进行了改写，提出用ρf作为流体指示因子，并且指出流体项f可以直接作为一项流体因子参与流体检测和储层预测。Hilterman(2009)结合实例应用重点研究了非固结砂岩的流体因子敏感性，并指出决定非固结砂岩储层流体识别敏感性的关键是纵横波速度的选取。

现阶段常用流体因子的计算方式是基于弹性参数的间接组合运算，此类流体识别的质量主要取决于两方面，一是基础弹性参数反演是否可靠；二是构建的流体因子对孔隙流体类型是否敏感。考虑到叠前地震反演是提取基础弹性参数的主要手段，通过改善反演方法可以较好的提高弹性参数可靠性；另一方面，基于间接组合的流体因子不可避免的会造成累计误差，而流体因子直接反演可以较好的解决这个问题，基于双相介质岩石物理理论构建敏感流体因子，通过研究其与地震动反射特征的内在联系，利用叠前地震资料即可实现流体因子直接反演，从而提高了流体因子的指示敏感性与其估算的可靠性。为此我们发明了一种新的频变粘弹性流体因子叠前地震反演方法，解决了以上技术问题。

发明内容

本发明的目的是提供一种频变粘弹性流体因子叠前地震反演方法，降低了流体检测的多解性，消除了仅仅考虑振幅进行流体检测时出现的识别假象。

本发明的目的可通过如下技术措施来实现：频变粘弹性流体因子叠前地震反演方法，该频变粘弹性流体因子叠前地震反演方法包括：步骤1，构建粘弹性介质反射特征方程并分析该方程的合理性；步骤2，进行地震资料多尺度分解；步骤3，针对每个尺度或频率内的地震反射记录，开展频变弹性阻抗反演方法；步骤4，确定频变加权系数，并在此基础上提取频变粘弹性流体因子；步骤5，基于频变粘弹性流体因子，对流体的分布状态进行预测。

本发明的目的还可通过如下技术措施来实现：

在步骤1中，通过粘弹性介质参数来表征与储层流体相关的特性，结合岩石弹性模量估算方法获取频变地震岩石物理模型，在此基础上研究含流体孔隙介质地震波粘弹性衰减，从而明确敏感反映储层及孔隙流体信息的有效频变流体敏感参数，构建频变粘弹性流体因子；在岩石物理理论的指导下，获取与孔隙流体特征相关并具有表征流体性质的频变弹性参数组合，构建粘弹性介质反射特征方程并分析该方程的合理性。

在步骤1中，所构建的粘弹性介质反射特征方程为：

其中，为粘弹性情况下角度为θ、频率为ω时的反射系数，α和β分别为频率ω对应的纵波相速度和横波相速度，ρ代表密度，ω_r为参考频率，表示饱和流体岩石纵横波速度比的平方，表示干岩石纵横波速度比的平方，f_anela,μ和ρ分别表示粘弹性介质两侧的频变粘弹流体因子，剪切模量和密度的平均值；Δf_anela(ω),ΔμΔρ则分别表示界面两侧的频变粘弹流体因子，剪切模量和密度的差值，Q_P为纵波的品质因子。

在步骤2中，为了利用多尺度地震资料进行频散属性的反演方法研究，通过井旁道地震记录时频分析，结合目的层的频谱关系选择合适的尺度因子，选择优势频率并通过分频处理将地震记录分解为不同频带的地震信息；借助小波变换分频方法具有一定自适应分辨分析特性，将地震记录分解为一系列具有质心频率的窄带剖面，实现不同频带范围的地震振幅信息分离。

在步骤3中，采用稀疏约束弹性阻抗反演策略，考虑叠前地震资料包含更加丰富的振幅和频率信息，利用频变粘弹性流体因子弹性阻抗方程，在贝叶斯反演框架下，建立了频变弹性阻抗反演目标函数，发展弹性阻抗叠前地震反演优化算法，实现不同角度频变弹性阻抗反演估计。

在步骤4中，通过实际测井曲线与其对应位置的频变弹性阻抗曲线，确定频变加权系数；在此基础上结合频变弹性阻抗方程与多频段角度弹性阻抗数据，实现整个工区的频变属性提取，最终实现频变粘弹性流体因子的提取。

通过地震岩石物理研究表明，地下介质孔隙含流体类型的差异会引起地震波的衰减和频散特征，尤其在含烃储层中传播时，这种衰减现象更明显。本发明联合地震振幅、频率及偏移距等信息，直接利用地震波粘弹性参数的频散程度进行储层预测和流体表征。本发明中的频变粘弹性流体因子叠前地震反演方法，首先考虑复杂地下介质的粘弹性质构建频变粘弹性流体因子，并推导利用该流体敏感性强的频变粘弹性流体识别因子表征的地震AVOF反射特征方程；此外联合现代谱分解方法挖掘地震数据中包含的不同频带地震响应；最终提出了基于频变弹性阻抗反演的频变粘弹性流体因子直接提取方法，对提高储层流体性质的判识精度具有较强的指导意义。

本发明中的基于频变粘弹性流体因子的叠前储层流体识别方法，构建频变粘弹性流体因子敏感参数，在考虑介质粘弹性的情况下，更为精确地模拟了地震波在复杂地下介质中传播的过程。该方法利用叠前地震资料中蕴含的振幅和频率变化关系，将谱分解方法和叠前弹性阻抗反演理论相结合，能够实现频变粘弹性特征参数的提取，且表征不同程度衰减和频散规律的频变粘弹性流体因子，一定程度上降低了流体检测的多解性，消除了仅仅考虑振幅进行流体检测时出现的识别假象。

附图说明

图1为本发明的频变粘弹性流体因子叠前地震反演方法的一具体实施例的流程图；

图2为本发明的频变弹性阻抗的粘弹性频变流体因子与不同流体因子的敏感性分析图；

图3为本发明的基于频变弹性阻抗的粘弹性频变流体因子反射系数方程与不同反射系数方程的敏感性分析图；

图4为本发明的过埕北252埕北255连井线馆陶组小角度地震数据的示意图；

图5为本发明的过埕北252埕北255连井线馆陶组中角度地震数据的示意图；

图6为本发明的过埕北252埕北255连井线馆陶组大角度地震数据的示意图；

图7为本发明的过埕北252埕北255连井线馆陶组叠前小角度单频地震数据20Hz单频资料的示意图；

图8为本发明的过埕北252埕北255连井线馆陶组叠前小角度单频地震数据30Hz单频资料的示意图；

图9为本发明的过埕北252埕北255连井线馆陶组叠前小角度单频地震数据40Hz单频资料的示意图；

图10为本发明的过埕北252埕北255连井线馆陶组叠前小角度单频地震数据50Hz单频资料的示意图；

图11为本发明的过埕北252埕北255连井线馆陶组叠前小角度20Hz单频弹性阻抗的示意图；

图12为本发明的过埕北252埕北255连井线馆陶组叠前小角度30Hz单频弹性阻抗的示意图；

图13为本发明的过埕北252埕北255连井线馆陶组叠前小角度40Hz单频弹性阻抗的示意图；

图14为本发明的过埕北252埕北255连井线馆陶组叠前小角度50Hz单频弹性阻抗的示意图；

图15为本发明的过埕北252埕北255连井线馆陶组叠频变流体因子反演结果的示意图。

具体实施方式

为使本发明的上述和其他目的、特征和优点能更明显易懂，下文特举出较佳实施例，并配合附图所示，作详细说明如下。

如图1所示，图1为本发明基于频变粘弹性流体因子的储层流体识别的流程图，主要分为粘弹性介质特征方程的构建、地震资料的多尺度分解、频变弹性阻抗反演、频变粘弹性流体因子提取以及应用这五个步骤。

步骤101，通过粘弹性介质参数来表征与储层流体相关的特性，结合岩石弹性模量估算方法获取频变地震岩石物理模型，在此基础上研究含流体孔隙介质地震波粘弹性衰减，从而明确敏感反映储层及孔隙流体信息的有效频变流体敏感参数，构建频变粘弹性流体因子。在岩石物理理论的指导下，获取与孔隙流体特征相关并具有表征流体性质的频变弹性参数组合，构建粘弹性介质反射特征方程并分析该方程的合理性。

近似常数Q模型描述的非弹性介质中纵波复速度和横波复速度分别表示为：

其中，V_P和V_S分别为参考频率ω_r对应的纵波相速度和横波相速度，α和β分别为频率ω对应的纵波相速度和横波相速度，Q_P和Q_S分别为纵波和横波的品质因子，i为虚数单位。

弹性介质流体因子表示为：或(I_p为纵波阻抗，I_s为横波阻抗，ρ代表密度，γ_dry代表干岩石纵横波速度比)。类似地，将非弹性介质的流体因子称为衰减流体因子，其形式为：

其中，f_anela即所求粘弹性流体因子，ρ代表密度，γ_dry代表干岩石纵横波速度比。与弹性部分对应的流体因子为非弹性对流体因子的影响相当于在弹性背景上施加的微小扰动，因此用Δf_Q表示衰减对流体因子的扰动。非弹性介质弹性阻抗和衰减流体因子为从地震数据提取非弹性介质参数提供了新的途径。频变流体因子是指在岩石物理理论的指导下，构建的与孔隙流体特征相关并具有表征流体性质的频变弹性参数组合。图2为8种弹性参数随着频率的变化特征，为了明确弹性参数随频率的变化梯度大小，采用归一化弹性参数的频变程度对比分析，图中可见，Gassman流体项在粘弹性介质中的频散程度在众多弹性参数中最大，因此考虑利用粘弹性Gassman流体因子作为储层流体识别的敏感参数，进而增强流体识别的可靠性和稳定性。

基于频变粘弹性流体因子，推导粘弹性介质特征方程，研究地震反射特征随频率、入射角、介质参数等之间的变化特征，分析不同因素对地震响应的贡献，进一步明确根据叠前地震反射数据获取频变粘弹性流体因子的合理性。

基于Biot-Gassmann理论，Russell等人对饱含流体多孔介质的AVO理论进行了研究，具体推导过程如下所示：

已知粘弹性介质中的Aki-Richard近似表示为：

其中，为粘弹性情况下角度为θ、频率为ω时的反射系数，α和β分别为频率ω对应的纵波相速度和横波相速度，ρ代表密度，Δα、Δβ和Δρ分别为纵波相速度、横波相速度和密度的变化值，表示饱和流体岩石纵横波速度比的平方。公式(4)两边同乘以α²ρ，可以其变形为：

考虑到粘弹性介质中Gassmann流体项的表达公式为：

因此：

即：

将(8)式代入(5)式，可以得到：

化简之后可得：

又因为：

所以：

其中μ为剪切模量，Δμ为界面两侧剪切模型的变化值。

将(12)式代入(10)式，方程两边同除以α²ρ，可以得到：

由于：

所以：

又因为：

所以：

将(15)和(17)代入(13)，最终得到如下所示反射系数近似公式：

其中，f_anela(ω),μ和ρ分别表示粘弹性介质两侧的Gassmann流体项，剪切模量和密度的平均值；Δf_anela(ω),ΔμΔρ则分别表示界面两侧的Gassmann流体项，剪切模量和密度的差值。

将常Q衰减模型中的复速度用相速度表示，代入上式，得：

为粘弹介质下饱含流体岩石的纵横波速度比平方，假设干岩石纵横波速度比不受衰减影响。将饱和岩石的纵横波速度比代入反射系数表达式，得

粘弹性纵波反射系数中含有虚部项，由于假设弱非弹性性质，因此，虚部项为相对小量，将其忽略可得到反射系数为：

式(23)即为用频变粘弹性流体因子来替换弹性流体因子推导所得的粘弹性介质特征方程。为研究非弹性性质对反射系数的影响，通过实验模拟验证。图3为粘弹性近似反射系数与精确反射系数方程随着频率和入射角度变化的平面图，其中蓝线为粘弹性介质中精确Rpp-Zoeppritz方程频变反射信息，虚线为粘弹性介质中Rpp-Aki近似方程频变反射信息，红线为本文粘弹性介质Rpp-Approxi近似方程频变反射信息。图中可见，在粘弹性介质中(近似常Q模型)，精确Zoeppritz方程、Aki-Rhichard近似方程、Aki-approxi近似方程随着角度和频率的变化趋势基本保持一致，验证了粘弹性近似方程Aki-approxi的有效性，因此利用本专利推导的近似方程发展粘弹性介质频变流体因子叠前反演具有较强的合理性。

在推导出频变粘弹性流体因子的基础上，利用叠前地震频变响应特征方程，获得频变粘弹性流体因子弹性阻抗方程，研究弹性阻抗随频率及入射角的变化特征。从而在贝叶斯反演框架下，建立频变弹性阻抗反演目标函数，发展弹性阻抗叠前地震反演优化算法,实现基于频变粘弹性流体因子叠前地震反演的储层流体识别。

借鉴Connolly推导弹性阻抗的思想，用频变弹性阻抗表示反射系数，得到：

其中，QEI(θ,ω)为单频弹性阻抗，ΔQEI为界面两侧单频弹性阻抗的相对变化，Δln[QEI(θ,ω)]为粘弹性参数相对变化的对数形式，将粘弹性反射特征方程(23)代入上式，得到：

将式(25)简化，可以得到：

其中：

根据对频变流体因子的频散特征分析可知，频变流体因子随着频率的变化呈现出近似线性变化，因此假设频变流体因子在参考频率附近可以进行一阶泰勒展开:

因此参考频率ω₀附近的频变弹性阻抗为：

相应任意频率ω附近的频变弹性阻抗为

将式(31)-式(30)可近似得到：

将粘弹性弹性参数的相对变化量用对数形式表示，

由于：

根据式(34)和式(35)，式(33)可以简化得到，

进一步写成指数函数的形式，可得

两边取积分并将其指数化，消掉等式两边的微分项和对数项，进一步取积分常数为0，得到：

QEI(θ,ω)≈QEI(θ,ω₀)I_f ^a(θ,ω)ΔωI_μ ^b(θ,ω)Δω (38)

即：

其中，a(θ,ω)，b(θ,ω)分别表示为与频率相关的系数，因此在考虑频变粘弹性流体因子反演的过程中，需要输入不同频带的弹性阻抗反演信息，只有如此才能体现频率差异性响应，由公式(39)知，粘弹性频变流体因子与阻抗之间是非线性指数关系，为了简化求解方法，对公式进行对数变换，将其转化为线性形式：

该公式即所求频变弹性阻抗反演目标函数。

步骤102，地震资料的多尺度分解，为了能够较好地利用多尺度地震资料进行频散属性的反演方法研究，本次研究通过井旁道地震记录时频分析，结合目的层的频谱关系选择合适的尺度因子，选择优势频率并通过分频处理将地震记录分解为不同频带的地震信息。借助小波变换分频方法具有一定自适应分辨分析特性，可将地震记录分解为一系列具有质心频率的窄带剖面，实现不同频带范围的地震振幅信息分离。

如图4-6所示，为过埕北252埕北255连井线馆陶组的远中近三种角度的叠前地震资料，而图7-10则为过埕北252埕北255连井线馆陶组叠前小角度单频地震数据(12)(a)20Hz单频资料(b)30Hz单频资料(c)40Hz单频资料(d)50Hz单频资料。通过井旁道地震记录，结合目的层频谱确定优势频率，通过小波变换分频方法将地震记录分解为一系列具有质心频率的窄频带剖面，进而实现不同频带范围的地震振幅信息分离，为频变粘弹性流体因子叠前地震反演识别地下流体打下坚实的数据基础。

步骤103，针对每个尺度(频率)内的地震反射记录，开展频变弹性阻抗反演方法，与常规全频带弹性阻抗反演方法相同，本专利采用稀疏约束弹性阻抗反演策略。从而在多频段响应的基础上获取多频段弹性阻抗。

考虑叠前地震资料包含更加丰富的振幅和频率信息，利用频变粘弹性流体因子弹性阻抗方程，在贝叶斯反演框架下，建立了频变弹性阻抗反演目标函数，发展弹性阻抗叠前地震反演优化算法，实现不同角度频变弹性阻抗反演估计。如图11-14即过埕北252埕北255连井线馆陶组叠前小角度不同单频弹性阻抗反演结果。

步骤104，即频变粘弹性流体因子的提取步骤。通过实际测井曲线与其对应位置的频变弹性阻抗曲线，确定频变加权系数。在此基础上结合频变弹性阻抗方程与多频段角度弹性阻抗数据，实现整个工区的频变属性提取，最终实现频变粘弹性流体因子的提取。

由公式(40)知道，要得到两个粘弹性频变参数的对数表达式，至少需要两个不同角度的弹性阻抗体和两个不同频率的弹性阻抗体，因此项目组假设N个入射角度，M个频率信息，可构建如下方程组。

受工区纵横波速度比以及频率选择差异的影响，若直接求解该方程组，得到的lnI_fane和lnI_μ数值极有可能与弹性参数的意义相悖，为了提高参数计算的可靠性，可以采用如下方法首先计算方程中弹性参数的系数。

考虑到入射角度和频率选择相同时，各采样点所对应的系数是相同的，其不随时间发生变化。故可以考虑同一道地震数据的不同采样点(假设有n个采样点)，借助井旁道频变弹性阻抗反演结果与测井数据计算得到的频变粘弹性参数构建如下矩阵：

求解以上矩阵，即可获得适用于研究工区的系数a(θ,ω)Δω和b(θ,ω)Δω的数值，针对不同的入射角度和不同的频率选择分别计算，便可得到N个角度(分别为θ₁、θ₂和θ_N)和M个频率(分别为ω₁、ω₂和ω_M)对应的M*N个系数。将M*N个系数代入上述方程组，即可得到如下方程组：

其中，QEI(θ,ω)为单频弹性阻抗，lnI_fane、lnI_μ是将粘弹性流体因子与剪切模量的相对变化量用对数形式表示，而a(θ,ω)，b(θ,ω)分别表示为与频率相关的系数。通过求解以上方程组即可得到任意采样点处的频变粘弹性流体因子和频变粘弹性剪切模量。如图15，即过埕北252埕北255连井线馆陶组叠频变流体因子反演结果。

步骤105，基于频变粘弹性流体因子，对流体的分布状态进行预测，从而提高储层和油藏描述的准确性和可靠性，刻画储层含流体空间变化、降低勘探风险。

综上所述，粘弹性介质中频变粘弹性流体因子叠前地震直接反演方法以及流体识别应用流程，如图1所示。

常规频变流体因子的研究多是基于完全弹性假设展开的，但是由于地下介质往往并非完全弹性的，因此考虑介质粘弹性质可以更好的模拟地震波在复杂介质中的传播过程，对提高储层流体识别精度具有实际意义。该方法首先，结合Futterman近似恒Q模型构建频变特定的粘弹性频变流体因子，并对众多频变弹性参数的流体敏感性进行了分析；然后，推导了利用频变粘弹性流体因子表征的地震波反射系数方程，进而有助于研究地震反射特征随频率、入射角、介质参数等的变化特征；其次，综合谱分解方法将地震记录分解为不同频带的地震信息，为开展叠前频变粘弹性流体因子反演奠定数据基础；最后，推导出基于频变粘弹性流体因子的弹性阻抗方程，并在贝叶斯框架下建立起地震-弹性阻抗-频变粘弹性流体因子的直接反演流程，提出了频变粘弹性流体因子叠前地震反演方法。

本发明中的基于频变粘弹性流体因子的叠前储层流体识别方法，考虑含流体孔隙介质地震波衰减，构建频变粘弹性流体因子，结合研究区内关键井的岩性、物性和流体特征，研究频变粘弹性流体因子表征储层及孔隙流体变化的能力，明确频变粘弹性流体因子的流体敏感性。推导用频变粘弹性流体因子表示的地震波反射系数方程，研究地震反射特征随频率、入射角、介质参数等之间的变化特征，分析不同因素对地震响应的贡献，明确从叠前地震资料中获取频变粘弹性流体因子的合理性。通过井旁道地震记录时频分析，结合目的层的频谱关系选择合适的尺度因子，通过分频处理将地震记录分解为不同频带的地震信息，为开展叠前地震频变反演奠定理论及资料基础。利用叠前地震频变响应特征方程，推导出频变粘弹性流体因子弹性阻抗方程，研究弹性阻抗随频率及入射角的变化特征，在贝叶斯反演框架下，建立频变弹性阻抗反演目标函数，发展弹性阻抗叠前地震反演优化算法，研究频变粘弹性流体因子叠前地震反演方法，结合研究目标岩石物理分析，实现基于频变粘弹性流体因子叠前地震反演的储层流体识别。

Claims (6)

1.频变粘弹性流体因子叠前地震反演方法，其特征在于，该频变粘弹性流体因子叠前地震反演方法包括：

步骤1，构建粘弹性介质反射特征方程并分析该方程的合理性；

步骤2，进行地震资料多尺度分解；

步骤3，针对每个尺度或频率内的地震反射记录，开展频变弹性阻抗反演方法；

步骤4，确定频变加权系数，并在此基础上提取频变粘弹性流体因子；

步骤5，基于频变粘弹性流体因子，对流体的分布状态进行预测。

2.根据权利要求1所述的频变粘弹性流体因子叠前地震反演方法，其特征在于，在步骤1中，通过粘弹性介质参数来表征与储层流体相关的特性，结合岩石弹性模量估算方法获取频变地震岩石物理模型，在此基础上研究含流体孔隙介质地震波粘弹性衰减，从而明确敏感反映储层及孔隙流体信息的有效频变流体敏感参数，构建频变粘弹性流体因子；在岩石物理理论的指导下，获取与孔隙流体特征相关并具有表征流体性质的频变弹性参数组合，构建粘弹性介质反射特征方程并分析该方程的合理性。

3.根据权利要求2所述的频变粘弹性流体因子叠前地震反演方法，其特征在于，在步骤1中，所构建的粘弹性介质反射特征方程为：

其中，为粘弹性情况下角度为θ、频率为ω时的反射系数，α和β分别为频率ω对应的纵波相速度和横波相速度，ρ代表密度，ω_r为参考频率，表示饱和流体岩石纵横波速度比的平方，表示干岩石纵横波速度比的平方，f_anela,μ和ρ分别表示粘弹性介质两侧的频变粘弹流体因子，剪切模量和密度的平均值；Δf_anela(ω),ΔμΔρ则分别表示界面两侧的频变粘弹流体因子，剪切模量和密度的差值，Q_P为纵波的品质因子。

4.根据权利要求1所述的频变粘弹性流体因子叠前地震反演方法，其特征在于，在步骤2中，为了利用多尺度地震资料进行频散属性的反演方法研究，通过井旁道地震记录时频分析，结合目的层的频谱关系选择合适的尺度因子，选择优势频率并通过分频处理将地震记录分解为不同频带的地震信息；借助小波变换分频方法具有一定自适应分辨分析特性，将地震记录分解为一系列具有质心频率的窄带剖面，实现不同频带范围的地震振幅信息分离。

5.根据权利要求1所述的频变粘弹性流体因子叠前地震反演方法，其特征在于，在步骤3中，采用稀疏约束弹性阻抗反演策略，考虑叠前地震资料包含更加丰富的振幅和频率信息，利用频变粘弹性流体因子弹性阻抗方程，在贝叶斯反演框架下，建立了频变弹性阻抗反演目标函数，发展弹性阻抗叠前地震反演优化算法，实现不同角度频变弹性阻抗反演估计。

6.根据权利要求1所述的频变粘弹性流体因子叠前地震反演方法，其特征在于，在步骤4中，通过实际测井曲线与其对应位置的频变弹性阻抗曲线，确定频变加权系数；在此基础上结合频变弹性阻抗方程与多频段角度弹性阻抗数据，实现整个工区的频变属性提取，最终实现频变粘弹性流体因子的提取。