一种油气检测方法、装置、电子设备及存储介质与流程

本技术涉及油汽地球物理，具体地涉及一种油气检测方法、装置、电子设备及存储介质。

背景技术：

1、地震波在储层中往往会发生较明显的衰减。这种衰减和储层的流体性质密切相关。在含油气的储层中地震波衰减往往相对较大，而在含水的储层中地震波衰减相对较小。因此，可以利用地震波的衰减来预测地下储层流体。

2、许多学者研究了地层吸收衰减的提取方法，并根据提取出的地震波衰减属性来预测地下油气藏。1983年，在大量实验数据的基础上，学者taner通过连续小波变换法估算出了地层的吸收参数值。1991年，r.tonn对于前人所提出的不同方法进行了整理，如上升时间法、解析信号法、振幅衰减法等，而且，他还利用合成的vsp资料数据对这10种方法进行了详细的比较。得出了地震资料品质严重影响吸收系数计算精度的结论。1997年，jerrym.harris等通过实验研究了品质因子q值与质心频率之间的某种关系，根据q值与质心频率之间的这种关系，他们提出了另一种计算q值的方法，即质心频率偏移法(cfs)。进入新世纪以后，出现了许多以时频分析以及反演为理论基础的新的地层吸收特征参数的评价方法。2002年，张长军等人在研究了地震波峰值频率和主频之间的关系之后，提出了剥层法原理，并研究了cmp道集中的峰值频率求取算法(pfvo法)，利用此算法预测出地震波反射位置的峰值频率，通过与零落散步的反射系数的结合，然后给出了任意一个反射点的峰值频率，最终利用地震信号峰值频率的改变求取出地层品质因子q值，此法称为rgqa(reflectivity guided qanalysis)方法。2002年，王西文等人在理论研究的基础上，利用小波域中的分频处理技术估算出了子波频率谱的瞬时振幅的高频与低频的比值，从而可以利用此比值计算地层品质因子或吸收衰减系数。2003年，james等人对s变换分析其原理进行研究，通过测量雷达数据的中心频率的偏移计算吸收参数。2004年，李宏兵等给出了地震波吸收衰减的数学公式并提取了不同尺度下的衰减剖面，并将其用来刻画储层，该关系式是在理论研究的前提下，通过推导地震波的传播公式而得到的。地面地震的全波形广义反演技术是得到q的一种有效方法，这是sheng jie sun在他的博士论文中提到的。2008年，高静怀通过研究给出了一种函数，该函数包括四个待定的参数，并用其对地震源的子波进行了逼近，并结合粘弹性介质中的单程波原理，推导出小波包络峰值处瞬时频率与品质因子之间的解析关系。

3、近年来，一种时频域内的频率衰减梯度(董宁和杨立强，2008；张固澜，2011；付勋勋等，2012)被广泛用来表针地震衰减属性，并用来预测地下含油气情况，且取得不错的应用效果。频率衰减梯度一般采用高频端两点斜率来表征地震波衰减。对于频谱波动较大的地震信号，该方法提取的吸收衰减属性效果较差，从而导致油气检测的可靠性较低。另一方面，当地震记录中的高频或低频成分受干扰较强时，频率衰减梯度则不能很好地用于检测油气。

技术实现思路

1、有鉴于此，本技术提供一种油气检测方法、装置、电子设备及存储介质，以利于解决现有技术中频率衰减梯度则不能很好地用于检测油气的问题。

2、第一方面，本技术实施例提供了一种基于地震波综合频率衰减梯度的油气检测方法，包括：

3、在输入的叠后地震记录上逐道进行时频变换，获得地震道的时频域结果；

4、对所述地震道的时频域结果上的不同时刻点采用数值拟合方法得到不同时刻点的拟合子波谱；

5、在所述拟合子波谱上提取综合频率衰减梯度，所述综合频率衰减梯度由高频频率衰减梯度和低频频率衰减梯度共同确定；

6、根据所述综合频率梯度预测地下油气储层。

7、在一种可能的实现方式中，所述在输入的叠后地震记录上逐道进行时频变换，获得地震道的时频域结果，包括：

8、在输入的叠后地震记录上逐道进行广义s变换，获得地震道的时频域结果。

9、在一种可能的实现方式中，所述广义s变换的表达式为：

10、

11、其中，h(t)为待分析的地震信号，f为频率，t为时间，τ为分析的时间点，s(τ,f)为地震信号的时频谱，λ和p为参数，λ和p用于共同调节时间窗长度随频率变化的速度。

12、在一种可能的实现方式中，在所述在输入的叠后地震记录上逐道进行时频变换，获得地震道的时频域结果之后，所述方法还包括：

13、设定常数k和n，对所述地震道的时频域结果，采用最小二乘谱模拟方法进行子波谱振幅拟合，获得a和an，其中，a为拟合结果的参数，an为拟合结果中关于频率f的多项式系数。

14、在一种可能的实现方式中，所述数值拟合方法的表达式为：

15、

16、其中，|w(f)|为拟合子波谱。

17、在一种可能的实现方式中，所述在所述拟合子波谱上提取综合频率衰减梯度，包括：

18、分别在所述拟合子波谱上提取高频频率衰减梯度gh和低频频率衰减梯度gl；

19、将所述高频频率衰减梯度gh和所述低频频率衰减梯度gl相乘，获得综合频率衰减梯度g，g＝gh×gl。

20、在一种可能的实现方式中，所述高频频率衰减梯度gh的计算公式为：

21、

22、其中，a80，a20分别为80％，20％峰值振幅值；fh80，fh20分别为高频端80％，20％峰值振幅值对应的频率值。

23、在一种可能的实现方式中，所述低频频率衰减梯度gl的计算公式为：

24、

25、其中，a80，a20分别为80％，20％峰值振幅值；fl80，fl20分别为低频端80％，20％峰值振幅值对应的频率值。

26、第二方面，本技术实施例提供了一种基于地震波综合频率衰减梯度的油气检测装置，包括：

27、时频变换模块，用于在输入的叠后地震记录上逐道进行时频变换，获得地震道的时频域结果；

28、子波谱拟合模块，用于对所述地震道的时频域结果上的不同时刻点采用数值拟合方法得到不同时刻点的拟合子波谱；

29、综合频率衰减梯度提取模块，用于在所述拟合子波谱上提取综合频率衰减梯度，所述综合频率衰减梯度由高频频率衰减梯度和低频频率衰减梯度共同确定；

30、油气储层预测模块，用于根据所述综合频率梯度预测地下油气储层。

31、第三方面，本技术实施例提供了一种电子设备，包括：

32、处理器；

33、存储器；

34、以及计算机程序，其中所述计算机程序被存储在所述存储器中，所述计算机程序包括指令，当所述指令被所述处理器执行时，使得所述电子设备执行第一方面中任意一项所述的方法。

35、第四方面，本技术实施例提供了一种计算机可读存储介质，所述计算机可读存储介质包括存储的程序，其中，在所述程序运行时控制所述计算机可读存储介质所在设备执行第一方面中任意一项所述的方法。

36、在本技术实施例中，利用地震信号高、低频端信息来综合表征地震衰减的大小，提高衰减属性的提取精度和地下油气储层预测的精度。