流体流度预测方法、装置、终端设备和存储介质与流程

1.本发明涉及油气地震勘探数据处理与解释领域，特别涉及一种流体流度预测方法、装置、终端设备和存储介质。


背景技术：

2.如今的油气勘探目标已面向具有更加复杂的地质条件的地层，地下复杂的构造条件给储层预测带来了巨大的挑战，提高储层预测精度是如今油气勘探领域的关键。常规的弹性参数反演技术通常很难直接得到与储层直接相关的识别因子，并且受弹性异常的干扰较大、因此，直接从地震资料中提取与储层流体相关的参数来进行储层预测是提高储层预测精度的重要手段。
3.流体流度(渗透率与流体粘度之比)是直接反应储层渗透性与流体性质的参数，陈学华等提出了利用低频端地震资料的瞬时谱计算储层流体流度的方法。近年来的研究发现，依赖频率的avo(amplitude variation with offset,振幅随偏移距的变化)反演技术为利用叠前地震资料反演弹性参数的频散属性提供了有效手段，该技术充分考虑了含流体孔隙介质引起的频散衰减效应，可提取不同流体在不同频率下的频散异常，进而实现储层预测。
4.利用低频端地震资料的瞬时谱计算储层流体流度的方法，该方法仅利用了单个频率下的地震信号的信息，未充分利用依赖频率的地震反射信息，在反演时会遗漏其余频率成分的反射信息，进而导致反演结果的精确度降低；通常的叠前频散反演方法需要道集优化处理后高品质的叠前地震道集，叠前地震资料的品质会直接影响反演结果，并且计算量大，效率低。


技术实现要素：

5.为了至少部分地解决现有技术存在的技术问题，发明人做出本发明，通过具体实施方式，提供一种流体流度预测方法、装置、终端设备和存储介质。
6.第一方面，本发明实施例提供一种流体流度预测方法，包括：
7.根据获取的参考层位数据和目的层段叠后地震数据中的单个地震道数据，使用时频变换得到所述单个地震道的时频振幅谱；
8.从所述单个地震道的时频振幅谱中，选定预设数量个角频率，并根据地震信号主频，确定参考角频率；
9.对所述单个地震道的时频振幅谱进行谱均衡处理，得到单个地震道谱均衡后的数据；
10.根据所述单个地震道谱均衡后的数据和所述参考角频率，得到所述参考角频率时的纵波速度变化率；
11.根据所述纵波速度变化率、所述预设数量个角频率和所述参考角频率，得到所述参考角频率时的纵波频散数据；
12.根据所述纵波频散数据，得到所述目的层段的单个地震道的流体流度。
13.可选的，所述根据所述单个地震道谱均衡后的数据和所述参考角频率，得到所述参考角频率时的纵波速度变化率，包括：
14.利用最小二乘法计算参考频率下纵波变化率d
p
，计算结果为：
15.d
p
＝8/5*xb(t,f0)
16.式中，xb表示单个地震道谱均衡后的数据,t表示时间，f0表示和参考角频率对应的频率，ω0表示参考角频率。
17.可选的，所述根据所述纵波速度变化率、所述预设数量个角频率和所述参考角频率，得到所述参考角频率时的纵波频散数据，包括：
18.使用下式计算所述参考角频率时的纵波频散数据，
[0019][0020]
式中，xb表示单个地震道谱均衡后的数据,t表示时间，ω1，...，ωm表示所选定的m个角频率,ω0表示参考角频率,d
p
表示参考角频率ω0时的纵波速度变化率，i
p
表示参考角频率ω0时的纵波频散数据。
[0021]
可选的，所述根据所述纵波频散数据，得到所述目的层段的单个地震道的流体流度，包括：
[0022]
使用下式得到所述目的层段的单个地震道的流体流度，
[0023][0024]
式中，f为流体流度结果，ω0表示参考角频率，c为与流体密度相关的常数，i
p
表示参考角频率ω0时的纵波频散数据。
[0025]
可选的，所述使用下式得到所述目的层段的单个地震道的流体流度，
[0026][0027]
式中，f为流体流度结果，ω0表示参考角频率，c为与流体密度相关的常数，i
p
表示参考角频率ω0时的纵波频散数据，包括：
[0028]
基于流体流度与反射系数之间的关系建立依赖频率的流体流度的表达式：
[0029][0030]
式中，ω为角频率，v
p
为纵波速度，c为与流体密度相关的常数，f为流体流度；
[0031]ip
的表达式为：
[0032]
根据所述和所述得到
[0033]
式中，f为流体流度结果，ω0表示参考角频率，c为与流体密度相关的常数，i
p
表示参考角频率ω0时的纵波频散数据。
[0034]
可选的，所述基于流体流度与反射系数之间的关系建立依赖频率的流体流度的表达式：
[0035][0036]
式中，ω为角频率，v
p
为纵波速度，c为与流体密度相关的常数，f为流体流度，包括：
[0037]
当储层岩石包含流体时，依赖频率的反射系数公式可表示为：
[0038][0039]
式中，θ为入射角，v
p
为纵波速度，vs为横波速度，ω为角频率；
[0040]
当θ＝0时，b(θ)＝0，则有：
[0041]
对其求导可得：
[0042][0043]
流体流度的计算公式为：
[0044]
式中，κ为储层的渗透率，η为流体粘度，c为与流体密度相关的常数，ω为角频率，r为ω的函数；
[0045]
根据所述和所述得到依赖频率的流体流度的表达式：
[0046][0047]
式中，ω为角频率，v
p
为纵波速度，c为与流体密度相关的常数，f为流体流度。
[0048]
可选的，还包括：
[0049]
分别根据所述获取的参考层位数据和目的层段叠后地震数据中的每个地震道数据，分别得到所述目的层段的每个地震道的流体流度，从而获得目的层段的流体流度剖面。
[0050]
可选的，还包括：
[0051]
根据整个工区的所有叠后地震数据和参考层位数据，得到整个工区不同目的层段的储层流体流度剖面；
[0052]
根据所述整个工区不同目的层段的储层流体流度剖面，进行储层预测，并抽取沿层切片进行分析，得到整个工区不同目的层段的储层预测空间展布。
[0053]
第二方面，本发明实施例提供一种流体流度预测装置，包括：
[0054]
频率选定模块，用于对获取的参考层位数据和目的层段叠后地震数据中的单个地震道数据，使用时频变换得到所述单个地震道的时频振幅谱；从所述单个地震道的时频振幅谱中，选定预设数量个角频率，并根据地震信号主频，确定参考角频率；
[0055]
流体流度预测模块，用于对所述单个地震道的时频振幅谱进行谱均衡处理，得到单个地震道谱均衡后的数据；根据所述单个地震道谱均衡后的数据和所述参考角频率，得到所述参考角频率时的纵波速度变化率；根据所述纵波速度变化率、所述预设数量个角频率和所述参考角频率，得到所述参考角频率时的纵波频散数据；根据所述纵波频散数据，得到所述目的层段的单个地震道的流体流度。
[0056]
可选的，所述频率选定模块，包括：
[0057]
时频变换模块，用于对获取的参考层位数据和目的层段叠后地震数据中的单个地震道数据，使用时频变换得到所述单个地震道的时频振幅谱；
[0058]
角频率选定模块，用于从所述单个地震道的时频振幅谱中，选定预设数量个角频率；
[0059]
参考角频率确定模块，用于根据地震信号主频，确定参考角频率。
[0060]
可选的，所述流体流度预测模块，包括：
[0061]
谱均衡处理模块，用于对所述单个地震道的时频振幅谱进行谱均衡处理，得到单个地震道谱均衡后的数据；
[0062]
纵波速度变化率确定模块，用于利用最小二乘法计算参考频率下纵波变化率d
p
，计算结果为：
[0063]dp
＝8/5*xb(t,f0)
[0064]
式中，xb表示单个地震道谱均衡后的数据,t表示时间，f0表示和参考角频率对应的频率，ω0表示参考角频率；
[0065]
纵波频散数据确定模块，用于使用下式计算所述参考角频率时的纵波频散数据，
[0066][0067]
式中，xb表示单个地震道谱均衡后的数据,t表示时间，ω1，...，ωm表示所选定的m个角频率,ω0表示参考角频率,d
p
表示参考角频率ω0时的纵波速度变化率，i
p
表示参考角
频率ω0时的纵波频散数据；
[0068]
流体流度确定模块，用于使用下式得到所述目的层段的单个地震道的流体流度，
[0069][0070]
式中，f为流体流度结果，ω0表示参考角频率，c为与流体密度相关的常数，i
p
表示参考角频率ω0时的纵波频散数据。
[0071]
可选的，还包括：
[0072]
流体流度剖面建立模块，用于分别根据所述获取的参考层位数据和目的层段叠后地震数据中的每个地震道数据，分别得到所述目的层段的每个地震道的流体流度，从而获得目的层段的流体流度剖面；
[0073]
储层预测空间展布建立模块，用于根据整个工区的所有叠后地震数据和参考层位数据，得到整个工区不同目的层段的储层流体流度剖面；根据所述整个工区不同目的层段的储层流体流度剖面，进行储层预测，并抽取沿层切片进行分析，得到整个工区不同目的层段的储层预测空间展布。
[0074]
基于同一发明构思，本发明实施例还提供一种终端设备，包括：存储器、处理器及存储于存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现前述的流体流度预测方法。
[0075]
基于同一发明构思，本发明实施例还提供一种计算机存储介质，所述计算机存储介质中存储有计算机可执行指令，所述计算机可执行指令执行时实现前述的流体流度预测方法。
[0076]
本发明实施例提供的上述技术方案的有益效果至少包括：
[0077]
本发明实施例提供的流体流度预测方法、装置、存储介质、设备，充分利用依赖频率的反演技术从叠后地震数据中提取储层流体流度信息，本发明直接利用叠后数据进行反演可避免叠前道集处理的步骤，并且直接利用叠后数据，使得计算量大大减小，提升了流体流度预测效率，降低了流体流度预测的时间成本和叠前道集处理成本，适用于依据大规模的三维地震数据进行预测；流体流度预测方法充分利用了地震数据的多种频率信息，该信息可反映地震波穿过含流体储层时引起的频散效应，反演时得到的结果也包含有纵波频散数据的贡献，可进一步提高储层预测的精度；充分挖掘了地震信号的频率信息，得到的流体流度结果受与储层无关的异常干扰很弱，对储层具有很强的敏感性，可进一步提高储层预测的精度；利用流体流度预测数据，可以得到流体流度剖面，并进一步得到储层的空间展布。
[0078]
本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
[0079]
下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
附图说明
[0080]
附图用来提供对本发明的进一步理解，并且构成说明书的一部分，与本发明的实施例一起用于解释本发明，并不构成对本发明的限制。在附图中：
[0081]
图1为本发明实施例中一种流体流度预测方法流程图；
[0082]
图2为本发明实施例中x井的叠后地震剖面图；
[0083]
图3为本发明实施例中x井的流体流度属性剖面图；
[0084]
图4为本发明实施例中一种流体流度预测装置框图；
[0085]
图5为本发明实施例中一种终端设备的结构示意图。
具体实施方式
[0086]
下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例，然而应当理解，可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更透彻地理解本公开，并且能够将本公开的范围完整的传达给本领域的技术人员。
[0087]
为了解决现有技术中存在的的问题，本发明实施例提供一种流体流度预测方法、装置、终端设备和存储介质。
[0088]
实施例一
[0089]
本发明实施例一提供一种流体流度预测方法，其流程如图1所示，包括如下步骤：
[0090]
步骤s101：对获取的参考层位数据和目的层段的叠后地震数据中的单个地震道数据，使用时频变换得到所述单个地震道的时频振幅谱；
[0091]
步骤s102：从所述单个地震道的时频振幅谱中，选定预设数量个角频率，并根据地震信号主频，确定参考角频率；
[0092]
步骤s103：对所述单个地震道的时频振幅谱进行谱均衡处理，得到单个地震道谱均衡后的数据；
[0093]
步骤s104：根据所述单个地震道谱均衡后的数据和所述参考角频率，得到所述参考角频率时的纵波速度变化率；
[0094]
可选的，利用最小二乘法计算参考频率下纵波变化率d
p
，计算结果为：
[0095]dp
＝8/5*xb(t,f0)
[0096]
式中，xb表示单个地震道谱均衡后的数据,t表示时间，f0表示和参考角频率对应的频率，ω0表示参考角频率。
[0097]
步骤s105：根据所述纵波速度变化率、所述预设数量个角频率和所述参考角频率，得到所述参考角频率时的纵波频散数据；
[0098]
可选的，使用下式计算所述参考角频率时的纵波频散数据，
[0099][0100]
式中，xb表示单个地震道谱均衡后的数据,t表示时间，ω1，...，ωm表示所选定的m个角频率,ω0表示参考角频率,d
p
表示参考角频率ω0时的纵波速度变化率，i
p
表示参考角频率ω0时的纵波频散数据。
[0101]
步骤s106：根据所述纵波频散数据，得到所述目的层段的单个地震道的流体流度。
[0102]
可选的，使用下式得到所述目的层段的单个地震道的流体流度，
[0103][0104]
式中，f为流体流度结果，ω0表示参考角频率，c为与流体密度相关的常数，i
p
表示参考角频率ω0时的纵波频散数据。
[0105]
可选的，基于流体流度与反射系数之间的关系建立依赖频率的流体流度的表达式：
[0106][0107]
式中，ω为角频率，v
p
为纵波速度，c为与流体密度相关的常数，f为流体流度；
[0108]ip
的表达式为：
[0109]
根据所述和所述得到
[0110]
式中，f为流体流度结果，ω0表示参考角频率，c为与流体密度相关的常数，i
p
表示参考角频率ω0时的纵波频散数据。
[0111]
可选的，当储层岩石包含流体时，依赖频率的反射系数公式可表示为：
[0112][0113]
式中，θ为入射角，v
p
为纵波速度，vs为横波速度，ω为角频率；
[0114]
当θ＝0时，b(θ)＝0，则有：
[0115]
对其求导可得：
[0116][0117]
流体流度的计算公式为：
[0118]
式中，κ为储层的渗透率，η为流体粘度，c为与流体密度相关的常数，ω为角频率，r为ω的函数；
[0119]
根据所述和所述得到依赖频率的流体流度的表达式：
[0120][0121]
式中，ω为角频率，v
p
为纵波速度，c为与流体密度相关的常数，f为流体流度。
[0122]
步骤s107：分别根据所述获取的参考层位数据和目的层段叠后地震数据中的每个地震道数据，分别得到所述目的层段的每个地震道的流体流度，从而获得目的层段的流体流度剖面。
[0123]
步骤s108：根据整个工区的所有叠后地震数据和参考层位数据，得到整个工区不同目的层段的储层流体流度剖面；根据所述整个工区不同目的层段的储层流体流度剖面，进行储层预测，并抽取沿层切片进行分析，得到整个工区不同目的层段的储层预测空间展布。
[0124]
例如，获取x井的叠后地震数据x(t,n)和参考层位数据h，其中t表示时间，n代表地震道序号；
[0125]
获取得到的叠后地震数据x(t,n)中序号为i的单个地震道x(t)，使用时频变换求取时频振幅谱x(t,f)，其中f表示频率，频率的范围为1～100hz；
[0126]
从时频振幅谱x(t,f)中，选定m个角频率ω1,...,ωm，并根据地震信号主频，确定参考角频率ω0，对时频振幅谱x(t,f)进行谱均衡处理得到单个地震道谱均衡后的数据xb(t,f)；
[0127]
利用谱均衡后的数据xb(t,f)，利用最小二乘法计算参考角频率ω0时的纵波速度变化率d
p
，
[0128]dp
＝8/5*xb(t,f0)
[0129]
式中，f0表示和参考角频率对应的频率，
[0130]
根据纵波速度变化率d
p
，进行依赖频率的反演，计算参考角频率ω0时纵波频散结果i
p
；
[0131]
根据参考角频率ω0时纵波频散结果i
p
，计算参考角频率ω0时的流体流度结果f；
[0132]
使用上述方法对叠后剖面的每个地震道数据进行处理，分别得到所述目的层段的每个地震道的流体流度，从而获得目的层段的流体流度剖面。图2为x井的叠后地震剖面图，横轴表示地震道的道号，纵轴表示时间，单位秒，图中的位置颜色越深代表振幅能量越强，该数据共有131道，时间范围为0.3～0.8秒，时间采样间隔为0.002秒，图中黑色虚线圈标注了含气储层的位置(约0.6～0.7s之间)，插入的曲线为纵波速度曲线，纵波速度曲线在图中越往左表示纵波速度越大，越往右表示纵波速度越小，可见，在含气储层位置处，纵波速度表现出明显的低速异常。低速异常指的是，地下介质中含油气时往往会造成地震波传播速度的下降,形成低速异常。地震资料中由于浅层气云、气烟囱、泥浆带、低速异常体等高吸收异常带的存在,会经常引起地震波传播速度异常现象。
[0133]
图3为利用本发明技术得到的流体流度属性剖面图，横轴表示地震道的道号，纵轴
表示时间，单位秒，在被黑线圈标注的含气储层位置处，流体流度属性表现出明显的高流度特征，很好地刻画了储层的位置和展布，与纵波速度曲线测井曲线的异常吻合度很高，纵波速度曲线是测井曲线的一种，可用来进行储层油气分布分析。图3排除了无关弹性层的图像，说明受无关弹性层的异常干扰很弱并且对储层预测的精度很高。因此，可利用该技术对含流体储层进行检测，进而为后续的地震解释和钻井部署提供有力支撑。
[0134]
更进一步的，根据整个工区的所有叠后地震数据，得到整个工区不同目的层段的储层流体流度剖面；根据所述整个工区不同目的层段的储层流体流度剖面，进行储层预测，并抽取沿层切片进行分析，得到整个工区不同目的层段的储层预测平面展布。处理研究工区所有的三维地震数据，对其进行储层预测并抽取沿层切片进行分析，得到整个工区不同目的层段的储层预测空间展布，也可以根据需要，基于二维地震数据，得到整个工区不同目的层段的储层预测平面展布。
[0135]
更进一步的，所述参考角频率ω0选取地震信号的主频。
[0136]
通过下式得到参考角频率ω0时纵波频散结果i
p
：
[0137]
通过下式得到参考角频率ω0时的流体流度结果f：
[0138]
式中，f为流体流度结果，c为与流体密度相关的常数。
[0139]
更进一步地，所述计算参考角频率ω0时的流体流度结果f的计算式，是通过下述过程得到的：
[0140]
基于流体流度与反射系数之间的关系建立依赖频率的流体流度的表达式：
[0141][0142]
式中，ω为角频率，v
p
为纵波速度，c为与流体密度相关的常数，f为流体流度；
[0143]ip
的表达式为：
[0144]
根据所述和所述得到
[0145]
上述的推导过程如下所示：
[0146]
当储层岩石包含流体时，地震波的速度随频率发生变化，则依赖频率的反射系数公式可表示为：
[0147][0148]
式中，θ为入射角，v
p
为纵波速度，vs为横波速度；
[0149]
当θ＝0时，b(θ)＝0，则有：
[0150]
对其求导可得：
[0151][0152]
流体流度的计算公式为：
[0153]
式中，κ为储层的渗透率，η为流体粘度，c为与流体密度相关的常数，ω为角频率，r为ω的函数；
[0154]
根据所述和所述得到依赖频率的流体流度的表达式：
[0155][0156]
式中，ω为角频率，v
p
为纵波速度，c为与流体密度相关的常数，f为流体流度。
[0157]
本实施例的上述方法中，充分利用依赖频率的反演技术从叠后地震数据中提取储层流体流度信息，本发明直接利用叠后数据进行反演可避免叠前道集处理的步骤，并且直接利用叠后数据，使得计算量大大减小，提升了流体流度预测效率，降低了流体流度预测的时间成本和叠前道集处理成本，适用于依据大规模的三维地震数据进行预测；流体流度预测方法充分利用了地震数据的多种频率信息，该信息可反映地震波穿过含流体储层时引起的频散效应，反演时得到的结果也包含有纵波频散数据的贡献，可进一步提高储层预测的精度；充分挖掘了地震信号的频率信息，得到的流体流度结果受与储层无关的异常干扰很弱，对储层具有很强的敏感性，可进一步提高储层预测的精度；利用流体流度预测数据，可以得到流体流度剖面，并进一步得到储层的空间展布。
[0158]
实施例二
[0159]
本发明实施例二提供一种流体流度预测装置，其结构如图4所示，包括如下：
[0160]
频率选定模块101，用于对获取的目的层段的叠后地震数据中的单个地震道数据，使用时频变换得到所述单个地震道的时频振幅谱；从所述单个地震道的时频振幅谱中，选定预设数量个角频率，并根据地震信号主频，确定参考角频率；
[0161]
流体流度预测模块102，用于对所述单个地震道的时频振幅谱进行谱均衡处理，得到单个地震道谱均衡后的数据；根据所述单个地震道谱均衡后的数据和所述参考角频率，得到所述参考角频率时的纵波速度变化率；根据所述纵波速度变化率、所述预设数量个角频率和所述参考角频率，得到所述参考角频率时的纵波频散数据；根据所述纵波频散数据，得到所述目的层段的单个地震道的流体流度；
[0162]
可选的，所述频率选定模块，包括：
[0163]
时频变换模块，用于对获取的参考层位数据和目的层段叠后地震数据中的单个地震道数据，使用时频变换得到所述单个地震道的时频振幅谱；
[0164]
角频率选定模块，用于从所述单个地震道的时频振幅谱中，选定预设数量个角频率；
[0165]
参考角频率确定模块，用于根据地震信号主频，确定参考角频率。
[0166]
可选的，所述流体流度预测模块，包括：
[0167]
谱均衡处理模块，用于对所述单个地震道的时频振幅谱进行谱均衡处理，得到单个地震道谱均衡后的数据；
[0168]
纵波速度变化率确定模块，用于利用最小二乘法计算参考频率下纵波变化率d
p
，计算结果为：
[0169]dp
＝8/5*xb(t,f0)
[0170]
式中，xb表示单个地震道谱均衡后的数据,t表示时间，f0表示和参考角频率对应的频率，ω0表示参考角频率；
[0171]
纵波频散数据确定模块，用于使用下式计算所述参考角频率时的纵波频散数据，
[0172][0173]
式中，xb表示单个地震道谱均衡后的数据,t表示时间，ω1，...，ωm表示所选定的m个角频率,ω0表示参考角频率,d
p
表示参考角频率ω0时的纵波速度变化率，i
p
表示参考角频率ω0时的纵波频散数据；
[0174]
流体流度确定模块，用于使用下式得到所述目的层段的单个地震道的流体流度，
[0175][0176]
式中，f为流体流度结果，ω0表示参考角频率，c为与流体密度相关的常数，i
p
表示参考角频率ω0时的纵波频散数据。
[0177]
可选的，所述流体流度预测装置，还包括：
[0178]
流体流度剖面建立模块103，用于分别根据所述获取的目的层段的叠后地震数据中的每个地震道数据，分别得到所述目的层段的每个地震道的流体流度，从而获得目的层段的流体流度剖面；
[0179]
储层预测空间展布建立模块104，用于根据整个工区的所有叠后地震数据，得到整
个工区不同目的层段的储层流体流度剖面；根据所述整个工区不同目的层段的储层流体流度剖面，进行储层预测，并抽取沿层切片进行分析，得到整个工区不同目的层段的储层预测空间展布。也可以根据需要，基于二维地震数据，得到整个工区不同目的层段的储层预测平面展布。
[0180]
基于同一发明构思，本发明实施例还提供一种终端设备，其结构如图5所示，包括：存储器、处理器及存储于存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现前述的流体流度预测方法。
[0181]
基于同一发明构思，本发明实施例还提供一种计算机存储介质，所述计算机存储介质中存储有计算机可执行指令，所述计算机可执行指令被处理器执行时实现前述的流体流度预测方法。
[0182]
关于上述实施例中的装置，其中各个模块执行操作的具体方式已经在有关该方法的实施例中进行了详细描述，此处将不做详细阐述说明。
[0183]
本发明实施例的上述方法和装置，充分利用依赖频率的反演技术从叠后地震数据中提取储层流体流度信息，本发明直接利用叠后数据进行反演可避免叠前道集处理的步骤，并且直接利用叠后数据，使得计算量大大减小，提升了流体流度预测效率，降低了流体流度预测的时间成本和叠前道集处理成本，适用于依据大规模的三维地震数据进行预测；流体流度预测方法充分利用了地震数据的多种频率信息，该信息可反映地震波穿过含流体储层时引起的频散效应，反演时得到的结果也包含有纵波频散数据的贡献，可进一步提高储层预测的精度；充分挖掘了地震信号的频率信息，得到的流体流度结果受与储层无关的异常干扰很弱，对储层具有很强的敏感性，可进一步提高储层预测的精度；利用流体流度预测数据，可以得到流体流度剖面，并进一步得到储层的空间展布。

技术特征：
1.一种流体流度预测方法，其特征在于，包括：根据获取的参考层位数据和目的层段叠后地震数据中的单个地震道数据，使用时频变换得到所述单个地震道的时频振幅谱；从所述单个地震道的时频振幅谱中，选定预设数量个角频率，并根据地震信号主频，确定参考角频率；对所述单个地震道的时频振幅谱进行谱均衡处理，得到单个地震道谱均衡后的数据；根据所述单个地震道谱均衡后的数据和所述参考角频率，得到所述参考角频率时的纵波速度变化率；根据所述纵波速度变化率、所述预设数量个角频率和所述参考角频率，得到所述参考角频率时的纵波频散数据；根据所述纵波频散数据，得到所述目的层段的单个地震道的流体流度。2.如权利要求1所述的方法，其特征在于，所述根据所述单个地震道谱均衡后的数据和所述参考角频率，得到所述参考角频率时的纵波速度变化率，包括：利用最小二乘法计算参考频率下纵波变化率d
p
，计算结果为：d
p
＝8/5*x
b
(t,f0)式中，x
b
表示单个地震道谱均衡后的数据,t表示时间，f0表示和参考角频率对应的频率，ω0表示参考角频率。3.如权利要求1所述的方法，其特征在于，所述根据所述纵波速度变化率、所述预设数量个角频率和所述参考角频率，得到所述参考角频率时的纵波频散数据，包括：使用下式计算所述参考角频率时的纵波频散数据，式中，x
b
表示单个地震道谱均衡后的数据,t表示时间，ω1，...，ω
m
表示所选定的m个角频率,ω0表示参考角频率,d
p
表示参考角频率ω0时的纵波速度变化率，i
p
表示参考角频率ω0时的纵波频散数据。4.如权利要求1所述的方法，其特征在于，所述根据所述纵波频散数据，得到所述目的层段的单个地震道的流体流度，包括：使用下式得到所述目的层段的单个地震道的流体流度，式中，f为流体流度结果，ω0表示参考角频率，c为与流体密度相关的常数，i
p
表示参考角频率ω0时的纵波频散数据。5.如权利要求4所述的方法，其特征在于，所述使用下式得到所述目的层段的单个地震道的流体流度，
式中，f为流体流度结果，ω0表示参考角频率，c为与流体密度相关的常数，i
p
表示参考角频率ω0时的纵波频散数据，包括：基于流体流度与反射系数之间的关系建立依赖频率的流体流度的表达式：式中，ω为角频率，v
p
为纵波速度，c为与流体密度相关的常数，f为流体流度；i
p
的表达式为：根据所述和所述得到式中，f为流体流度结果，ω0表示参考角频率，c为与流体密度相关的常数，i
p
表示参考角频率ω0时的纵波频散数据。6.如权利要求5所述的方法，其特征在于，所述基于流体流度与反射系数之间的关系建立依赖频率的流体流度的表达式：式中，ω为角频率，v
p
为纵波速度，c为与流体密度相关的常数，f为流体流度，包括：当储层岩石包含流体时，依赖频率的反射系数公式可表示为：式中，θ为入射角，v
p
为纵波速度，v
s
为横波速度，ω为角频率；当θ＝0时，b(θ)＝0，则有：对其求导可得：流体流度的计算公式为：
式中，κ为储层的渗透率，η为流体粘度，c为与流体密度相关的常数，ω为角频率，r为ω的函数；根据所述和所述得到依赖频率的流体流度的表达式：式中，ω为角频率，v
p
为纵波速度，c为与流体密度相关的常数，f为流体流度。7.如权利要求1所述的方法，其特征在于，还包括：分别根据所述获取的参考层位数据和目的层段叠后地震数据中的每个地震道数据，分别得到所述目的层段的每个地震道的流体流度，从而获得目的层段的流体流度剖面。8.如权利要求7所述的方法，其特征在于，还包括：根据整个工区的所有叠后地震数据和参考层位数据，得到整个工区不同目的层段的储层流体流度剖面；根据所述整个工区不同目的层段的储层流体流度剖面，进行储层预测，并抽取沿层切片进行分析，得到整个工区不同目的层段的储层预测空间展布。9.一种流体流度预测装置，其特征在于，包括：频率选定模块，用于对获取的参考层位数据和目的层段叠后地震数据中的单个地震道数据，使用时频变换得到所述单个地震道的时频振幅谱；从所述单个地震道的时频振幅谱中，选定预设数量个角频率，并确定参考角频率；流体流度预测模块，用于对所述单个地震道的时频振幅谱进行谱均衡处理，得到单个地震道谱均衡后的数据；根据所述单个地震道谱均衡后的数据和所述参考角频率，得到所述参考角频率时的纵波速度变化率；根据所述纵波速度变化率、所述预设数量个角频率和所述参考角频率，得到所述参考角频率时的纵波频散数据；根据所述纵波频散数据，得到所述目的层段的单个地震道的流体流度。10.如权利要求9所述的装置，其特征在于，所述频率选定模块，包括：时频变换模块，用于对获取的参考层位数据和目的层段叠后地震数据中的单个地震道数据，使用时频变换得到所述单个地震道的时频振幅谱；角频率选定模块，用于从所述单个地震道的时频振幅谱中，选定预设数量个角频率；参考角频率确定模块，用于根据地震信号主频，确定参考角频率。11.如权利要求9所述的装置，其特征在于，所述流体流度预测模块，包括：谱均衡处理模块，用于对所述单个地震道的时频振幅谱进行谱均衡处理，得到单个地震道谱均衡后的数据；纵波速度变化率确定模块，用于利用最小二乘法计算参考频率下纵波变化率d
p
，计算结果为：d
p
＝8/5*x
b
(t,f0)
式中，x
b
表示单个地震道谱均衡后的数据,t表示时间，f0表示和参考角频率对应的频率，ω0表示参考角频率；纵波频散数据确定模块，用于使用下式计算所述参考角频率时的纵波频散数据，式中，x
b
表示单个地震道谱均衡后的数据,t表示时间，ω1，...，ω
m
表示所选定的m个角频率,ω0表示参考角频率,d
p
表示参考角频率ω0时的纵波速度变化率，i
p
表示参考角频率ω0时的纵波频散数据；流体流度确定模块，用于使用下式得到所述目的层段的单个地震道的流体流度，式中，f为流体流度结果，ω0表示参考角频率，c为与流体密度相关的常数，i
p
表示参考角频率ω0时的纵波频散数据。12.如权利要求9所述的装置，其特征在于，还包括：流体流度剖面建立模块，用于分别根据所述获取的参考层位数据和目的层段叠后地震数据中的每个地震道数据，分别得到所述目的层段的每个地震道的流体流度，从而获得目的层段的流体流度剖面；储层预测空间展布建立模块，用于根据整个工区的所有叠后地震数据和参考层位数据，得到整个工区不同目的层段的储层流体流度剖面；根据所述整个工区不同目的层段的储层流体流度剖面，进行储层预测，并抽取沿层切片进行分析，得到整个工区不同目的层段的储层预测空间展布。13.一种终端设备，其特征在于，包括：存储器、处理器及存储于存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现权利要求1-8任一所述的流体流度预测方法。14.一种计算机存储介质，其特征在于，所述计算机存储介质中存储有计算机可执行指令，所述计算机可执行指令执行时实现权利要求1-8任一所述的流体流度预测方法。

技术总结
本发明公开了一种流体流度预测方法、装置、终端设备和存储介质。包括：根据获取的参考层位数据和目的层段叠后地震数据中的单个地震道数据，使用时频变换得到所述单个地震道的时频振幅谱；选定预设数量个角频率，并确定参考角频率；进行谱均衡处理，得到单个地震道谱均衡后的数据；根据所述单个地震道谱均衡后的数据和所述参考角频率，得到所述参考角频率时的纵波速度变化率；根据所述纵波速度变化率、所述预设数量个角频率和所述参考角频率，得到所述参考角频率时的纵波频散数据；根据所述纵波频散数据，得到流体流度。能够降低计算量，提高预测效率和精度，适于对大量地震数据的处理，受异常干扰很弱。受异常干扰很弱。受异常干扰很弱。


技术研发人员：张本健 严威 陈骁 周刚 钟原 田兴旺 李堃宇 杨岱林
受保护的技术使用者：中国石油天然气股份有限公司
技术研发日：2022.03.10
技术公布日：2023/9/20