含气性检测方法、装置、设备及存储介质

技术领域

本申请涉及石油勘探领域，特别地涉及一种含气性检测方法、装置、设备及存储介质。

背景技术

目前含气性检测方法大多为叠后衰减属性做定性预测，叠前反演做定量预测。对于叠后衰减属性进含气性定性预测，仅仅使用叠后地震数据缺少了叠前地震数据丰富的信息。叠前反演的含气性定量预测技术，是基于建立弹性参数与含气饱和度之间的关系进而完成含气饱和度的推算，或者通过神经网络算法进行学习，由于弹性参数中含气饱和度能引起的变化很小，因此叠前含气性定量预测的效果也并不理想。在进行含气性定量预测时，测井解释的含气饱和度解释结果的稳定性较差，因此在含气性预测时如果参考了测井解释的含气饱和度数据，预测结果的准确率也会降低。

发明内容

针对上述问题，本申请提供一种含气性检测方法、装置、设备及存储介质。

本申请提供了一种含气性检测方法，包括：

S1：收集叠前角度域道集，确定有效入射角范围以及确定叠前角度域道集有效的频带范围；

S2：按照入射角范围进行角度叠加，获得n个角度叠加后的地震数据，对应的入射角叠加范围的中心角度分别为angle1，angle2，angle3…anglen；

S3：对n个角度叠加后的地震数据的每一道计算振幅频率谱，获得地震数据有效的频带范围内的振幅频率谱面积area和振幅频率谱中能量最大值max_amp；

S4：应用振幅频率谱面积area和入射角叠加范围的中心角度，获得地震数据频谱面积在入射角方向的变化率数据体area_rate；

S5：应用振幅频率谱中能量最大值max_amp和入射角叠加范围的中心角度，获得振幅频率谱中能量最大值在入射角方向的变化率数据体max_amp_rate；

S6：将数据体area_rate和数据体max_amp_rate与实际工区已知产气井进行标定，完成含气性预测。

在一些实施例中，所述确定叠前角度域道集有效的频带范围之前，所述方法还包括：

进行剩余时间校正。

在一些实施例中，所述n不大于5。

在一些实施例中，所述应用振幅频率谱面积area和入射角叠加范围的中心角度，获得地震数据频谱面积在入射角方向的变化率数据体area_rate的具体方法包括：

令Y轴为n个振幅频率谱面积area，X轴为n个入射角叠加范围的中心角度angle1，angle2，angle3…anglen，通过曲线拟合，获得地震数据频谱面积在入射角方向的变化率数据体area_rate。

在一些实施例中，所述应用振幅频率谱中能量最大值max_amp和入射角叠加范围的中心角度，获得振幅频率谱中能量最大值在入射角方向的变化率数据体max_amp_rate的具体方法包括：

令Y轴为n个振幅频率谱中能量最大值max_amp，X轴为n个入射角叠加范围的中心角度angle1，angle2，angle3…anglen，通过曲线拟合，获得振幅频率谱中能量最大值在入射角方向的变化率数据体max_amp_rate。

在一些实施例中，所述有效入射角范围为0-36。

在一些实施例中，所述叠前角度域道集有效的频带范围为7-75Hz。

本申请实施例提供一种含气性检测装置，包括：

叠前角度域道集收集模块、中心角度计算模块、振幅频率谱数据计算模块、数据体area_rate计算模块、数据体max_amp_rate计算模块和含气性预测模块；

叠前角度域道集收集模块：收集叠前角度域道集，确定有效入射角范围以及确定叠前角度域道集有效的频带范围；

中心角度计算模块：按照入射角范围进行角度叠加，获得n个角度叠加后的地震数据，对应的入射角叠加范围的中心角度分别为angle1，angle2，angle3…anglen；

振幅频率谱数据计算模块：对n个角度叠加后的地震数据的每一道计算振幅频率谱，获得地震数据有效的频带范围内的振幅频率谱面积area和振幅频率谱中能量最大值max_amp；

数据体area_rate计算模块：应用振幅频率谱面积area和入射角叠加范围的中心角度，获得地震数据频谱面积在入射角方向的变化率数据体area_rate；

数据体max_amp_rate计算模块：应用振幅频率谱中能量最大值max_amp和入射角叠加范围的中心角度，获得振幅频率谱中能量最大值在入射角方向的变化率数据体max_amp_rate；

含气性预测模块：将数据体area_rate和数据体max_amp_rate与实际工区已知产气井进行标定，完成含气性预测。

本申请实施例提供一种含气性检测设备，包括存储器和处理器，所述存储器上存储有计算机程序，该计算机程序被所述处理器执行时，执行上述任意一项所述含气性检测方法。

本申请实施例提供一种存储介质，该存储介质存储的计算机程序，能够被一个或多个处理器执行，能够用来实现上述任一项所述含气性检测方法。

本申请提供的一种含气性检测方法、装置、设备及存储介质，具有如下有益效果：

(1)本发明完全基于地震数据进行计算，获得的结果确定性强；

(2)本发明在计算过程中不参考测井解释结果，不受人为因素影响。

附图说明

在下文中将基于实施例并参考附图来对本申请进行更详细的描述。

图1为本申请实施例提供的一种含气性检测方法的实现流程图；

图2为本申请实施例提供的剩余时差校正后的叠前角度域道集；

图3为本申请实施例提供的叠前角度域道集振幅频率谱；

图4为本申请实施例提供的实际工区数据体area_rate连井剖面图；

图5为本申请实施例提供的实际工区数据体max_amp_rate连井剖面图；

图6为本申请实施例提供的一种含气性检测方法的实现流程图。

在附图中，相同的部件使用相同的附图标记，附图并未按照实际的比例绘制。

具体实施方式

为了使本申请的目的、技术方案和优点更加清楚，下面将结合附图对本申请作进一步地详细描述，所描述的实施例不应视为对本申请的限制，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本申请保护的范围。

在以下的描述中，涉及到“一些实施例”，其描述了所有可能实施例的子集，但是可以理解，“一些实施例”可以是所有可能实施例的相同子集或不同子集，并且可以在不冲突的情况下相互结合。

如果申请文件中出现“第一第二第三”的类似描述则增加以下的说明，在以下的描述中，所涉及的术语“第一第二第三”仅仅是是区别类似的对象，不代表针对对象的特定排序，可以理解地，“第一第二第三”在允许的情况下可以互换特定的顺序或先后次序，以使这里描述的本申请实施例能够以除了在这里图示或描述的以外的顺序实施。

除非另有定义，本文所使用的所有的技术和科学术语与属于本申请的技术领域的技术人员通常理解的含义相同。本文中所使用的术语只是为了描述本申请实施例的目的，不是旨在限制本申请。

在介绍本申请实施例提供的一种含气性检测方法之前，对相关技术中存在的问题进行简单介绍：

目前含气性检测方法大多为叠后衰减属性做定性预测，叠前反演做定量预测。对于叠后衰减属性进含气性定性预测，仅仅使用叠后地震数据缺少了叠前地震数据丰富的信息。叠前反演的含气性定量预测技术，是基于建立弹性参数与含气饱和度之间的关系进而完成含气饱和度的推算，或者通过神经网络算法进行学习，由于弹性参数中含气饱和度能引起的变化很小，因此叠前含气性定量预测的效果也并不理想。在进行含气性定量预测时，测井解释的含气饱和度解释结果的稳定性较差，因此在含气性预测时如果参考了测井解释的含气饱和度数据，预测结果的准确率也会降低。

基于相关技术中存在的问题，本申请实施例提供一种含气性检测方法，所述方法应用于含气性检测设备，所述含气性检测设备可以为电子设备，例如计算机、移动终端等。本申请实施例提供的含气性检测方法所实现的功能可以通过电子设备的处理器调用程序代码来实现，其中，程序代码可以保存在计算机存储介质中。

实施例一

本申请实施例提供一种含气性检测方法，图6为本申请实施例提供的一种含气性检测方法的实现流程示意图，如图6所示，包括：

S1：收集叠前角度域道集，确定有效入射角范围以及确定叠前角度域道集有效的频带范围；

S2：按照入射角范围进行角度叠加，获得n个角度叠加后的地震数据，对应的入射角叠加范围的中心角度分别为angle1，angle2，angle3…anglen；

S3：对n个角度叠加后的地震数据的每一道计算振幅频率谱，获得地震数据有效的频带范围内的振幅频率谱面积area和振幅频率谱中能量最大值max_amp；

S4：应用振幅频率谱面积area和入射角叠加范围的中心角度，获得地震数据频谱面积在入射角方向的变化率数据体area_rate；

S5：应用振幅频率谱中能量最大值max_amp和入射角叠加范围的中心角度，获得振幅频率谱中能量最大值在入射角方向的变化率数据体max_amp_rate；

数据体area_rate代表的意义：由于受到储层含气的影响，地震数据的能量产生衰减，这种能量衰减的变化在振幅频率谱中体现为，随着入射角的增加能量衰减越强，即随着入射角的增加叠加数据频谱面积逐渐变小，受到含气性的影响地震数据的总能量随入射角增加减小，这里只考虑地震数据能量衰减与含气有关并未考虑储层厚度；数据体max_amp_rate代表的意义：由于受到储层含气的影响，地震数据的能量产生衰减，这种能量衰减的变化在振幅频率谱中体现为，随着入射角的增加能量衰减越强，即随着入射角的增加叠加数据振幅频率谱中能量最大值逐渐变小，受到含气的影响地震数据的峰值能量随入射角增加减小；

S6：将数据体area_rate和数据体max_amp_rate与实际工区已知产气井进行标定，完成含气性预测。

本申请提供的一种含气性检测方法，

(1)本发明完全基于地震数据进行计算，获得的结果确定性强；

(2)本发明在计算过程中不参考测井解释结果，不受人为因素影响。

实施例二

基于前述的实施例，本申请实施例再提供一种含气性检测方法，包括：

S21：收集叠前角度域道集，确定有效入射角范围以及确定叠前角度域道集有效的频带范围；

在一些实施例中，所述确定叠前角度域道集有效的频带范围之前，所述方法还包括：

进行剩余时间校正；

S22：按照入射角范围进行角度叠加，获得n个角度叠加后的地震数据，对应的入射角叠加范围的中心角度分别为angle1，angle2，angle3…anglen；

S23：对n个角度叠加后的地震数据的每一道计算振幅频率谱，获得地震数据有效的频带范围内的振幅频率谱面积area和振幅频率谱中能量最大值max_amp；

S24：应用振幅频率谱面积area和入射角叠加范围的中心角度，获得地震数据频谱面积在入射角方向的变化率数据体area_rate；

S25：应用振幅频率谱中能量最大值max_amp和入射角叠加范围的中心角度，获得振幅频率谱中能量最大值在入射角方向的变化率数据体max_amp_rate；

数据体area_rate代表的意义：由于受到储层含气的影响，地震数据的能量产生衰减，这种能量衰减的变化在振幅频率谱中体现为，随着入射角的增加能量衰减越强，即随着入射角的增加叠加数据频谱面积逐渐变小，受到含气性的影响地震数据的总能量随入射角增加减小，这里只考虑地震数据能量衰减与含气有关并未考虑储层厚度；数据体max_amp_rate代表的意义：由于受到储层含气的影响，地震数据的能量产生衰减，这种能量衰减的变化在振幅频率谱中体现为，随着入射角的增加能量衰减越强，即随着入射角的增加叠加数据振幅频率谱中能量最大值逐渐变小，受到含气的影响地震数据的峰值能量随入射角增加减小；

S26：将数据体area_rate和数据体max_amp_rate与实际工区已知产气井进行标定，完成含气性预测。

本申请提供的一种含气性检测方法，

(1)本发明完全基于地震数据进行计算，获得的结果确定性强；

(2)本发明在计算过程中不参考测井解释结果，不受人为因素影响。

实施例三

基于前述的实施例，本申请实施例再提供一种含气性检测方法，包括：

S31：收集叠前角度域道集，确定有效入射角范围以及确定叠前角度域道集有效的频带范围；

在一些实施例中，所述确定叠前角度域道集有效的频带范围之前，所述方法还包括：

进行剩余时间校正；

S32：按照入射角范围进行角度叠加，获得n个角度叠加后的地震数据，对应的入射角叠加范围的中心角度分别为angle1，angle2，angle3…anglen；

在一些实施例中，由于之后需要对角度叠加的地震数据提取属性，所以需要保证地震数据的信噪比不能太低，因此在进行角度叠加的时候，叠加的数据体个数一般不大于5个，所述n不大于5；

S33：对n个角度叠加后的地震数据的每一道计算振幅频率谱，获得地震数据有效的频带范围内的振幅频率谱面积area和振幅频率谱中能量最大值max_amp；

S34：应用振幅频率谱面积area和入射角叠加范围的中心角度，获得地震数据频谱面积在入射角方向的变化率数据体area_rate；

S35：应用振幅频率谱中能量最大值max_amp和入射角叠加范围的中心角度，获得振幅频率谱中能量最大值在入射角方向的变化率数据体max_amp_rate；

数据体area_rate代表的意义：由于受到储层含气的影响，地震数据的能量产生衰减，这种能量衰减的变化在振幅频率谱中体现为，随着入射角的增加能量衰减越强，即随着入射角的增加叠加数据频谱面积逐渐变小，受到含气性的影响地震数据的总能量随入射角增加减小，这里只考虑地震数据能量衰减与含气有关并未考虑储层厚度；数据体max_amp_rate代表的意义：由于受到储层含气的影响，地震数据的能量产生衰减，这种能量衰减的变化在振幅频率谱中体现为，随着入射角的增加能量衰减越强，即随着入射角的增加叠加数据振幅频率谱中能量最大值逐渐变小，受到含气的影响地震数据的峰值能量随入射角增加减小；

S36：将数据体area_rate和数据体max_amp_rate与实际工区已知产气井进行标定，完成含气性预测。

本申请提供的一种含气性检测方法，

(1)本发明完全基于地震数据进行计算，获得的结果确定性强；

(2)本发明在计算过程中不参考测井解释结果，不受人为因素影响。

实施例四

基于前述的实施例，本申请实施例再提供一种含气性检测方法，包括：

S41：收集叠前角度域道集，确定有效入射角范围以及确定叠前角度域道集有效的频带范围；

在一些实施例中，所述确定叠前角度域道集有效的频带范围之前，所述方法还包括：

进行剩余时间校正；

S42：按照入射角范围进行角度叠加，获得n个角度叠加后的地震数据，对应的入射角叠加范围的中心角度分别为angle1，angle2，angle3…anglen；

在一些实施例中，由于之后需要对角度叠加的地震数据提取属性，所以需要保证地震数据的信噪比不能太低，因此在进行角度叠加的时候，叠加的数据体个数一般不大于5个，所述n不大于5；

S43：对n个角度叠加后的地震数据的每一道计算振幅频率谱，获得地震数据有效的频带范围内的振幅频率谱面积area和振幅频率谱中能量最大值max_amp；

S44：应用振幅频率谱面积area和入射角叠加范围的中心角度，获得地震数据频谱面积在入射角方向的变化率数据体area_rate；

在一些实施例中，所述应用振幅频率谱面积area和入射角叠加范围的中心角度，获得地震数据频谱面积在入射角方向的变化率数据体area_rate的具体方法包括：

令Y轴为n个振幅频率谱面积area，X轴为n个入射角叠加范围的中心角度angle1，angle2，angle3…anglen，通过曲线拟合，获得地震数据频谱面积在入射角方向的变化率数据体area_rate；

S45：应用振幅频率谱中能量最大值max_amp和入射角叠加范围的中心角度，获得振幅频率谱中能量最大值在入射角方向的变化率数据体max_amp_rate；

数据体area_rate代表的意义：由于受到储层含气的影响，地震数据的能量产生衰减，这种能量衰减的变化在振幅频率谱中体现为，随着入射角的增加能量衰减越强，即随着入射角的增加叠加数据频谱面积逐渐变小，受到含气性的影响地震数据的总能量随入射角增加减小，这里只考虑地震数据能量衰减与含气有关并未考虑储层厚度；数据体max_amp_rate代表的意义：由于受到储层含气的影响，地震数据的能量产生衰减，这种能量衰减的变化在振幅频率谱中体现为，随着入射角的增加能量衰减越强，即随着入射角的增加叠加数据振幅频率谱中能量最大值逐渐变小，受到含气的影响地震数据的峰值能量随入射角增加减小；

S46：将数据体area_rate和数据体max_amp_rate与实际工区已知产气井进行标定，完成含气性预测。

本申请提供的一种含气性检测方法，

(1)本发明完全基于地震数据进行计算，获得的结果确定性强；

(2)本发明在计算过程中不参考测井解释结果，不受人为因素影响。

实施例五

基于前述的实施例，本申请实施例再提供一种含气性检测方法，包括：

S51：收集叠前角度域道集，确定有效入射角范围以及确定叠前角度域道集有效的频带范围；

在一些实施例中，所述确定叠前角度域道集有效的频带范围之前，所述方法还包括：

进行剩余时间校正；

S52：按照入射角范围进行角度叠加，获得n个角度叠加后的地震数据，对应的入射角叠加范围的中心角度分别为angle1，angle2，angle3…anglen；

在一些实施例中，由于之后需要对角度叠加的地震数据提取属性，所以需要保证地震数据的信噪比不能太低，因此在进行角度叠加的时候，叠加的数据体个数一般不大于5个，所述n不大于5；

S53：对n个角度叠加后的地震数据的每一道计算振幅频率谱，获得地震数据有效的频带范围内的振幅频率谱面积area和振幅频率谱中能量最大值max_amp；

S54：应用振幅频率谱面积area和入射角叠加范围的中心角度，获得地震数据频谱面积在入射角方向的变化率数据体area_rate；

在一些实施例中，所述应用振幅频率谱面积area和入射角叠加范围的中心角度，获得地震数据频谱面积在入射角方向的变化率数据体area_rate的具体方法包括：

令Y轴为n个振幅频率谱面积area，X轴为n个入射角叠加范围的中心角度angle1，angle2，angle3…anglen，通过曲线拟合，获得地震数据频谱面积在入射角方向的变化率数据体area_rate；

S55：应用振幅频率谱中能量最大值max_amp和入射角叠加范围的中心角度，获得振幅频率谱中能量最大值在入射角方向的变化率数据体max_amp_rate；

在一些实施例中，所述应用振幅频率谱中能量最大值max_amp和入射角叠加范围的中心角度，获得振幅频率谱中能量最大值在入射角方向的变化率数据体max_amp_rate的具体方法包括：

令Y轴为n个振幅频率谱中能量最大值max_amp，X轴为n个入射角叠加范围的中心角度angle1，angle2，angle3…anglen，通过曲线拟合，获得振幅频率谱中能量最大值在入射角方向的变化率数据体max_amp_rate；

数据体area_rate代表的意义：由于受到储层含气的影响，地震数据的能量产生衰减，这种能量衰减的变化在振幅频率谱中体现为，随着入射角的增加能量衰减越强，即随着入射角的增加叠加数据频谱面积逐渐变小，受到含气性的影响地震数据的总能量随入射角增加减小，这里只考虑地震数据能量衰减与含气有关并未考虑储层厚度；数据体max_amp_rate代表的意义：由于受到储层含气的影响，地震数据的能量产生衰减，这种能量衰减的变化在振幅频率谱中体现为，随着入射角的增加能量衰减越强，即随着入射角的增加叠加数据振幅频率谱中能量最大值逐渐变小，受到含气的影响地震数据的峰值能量随入射角增加减小；

S56：将数据体area_rate和数据体max_amp_rate与实际工区已知产气井进行标定，完成含气性预测。

本申请提供的一种含气性检测方法，

(1)本发明完全基于地震数据进行计算，获得的结果确定性强；

(2)本发明在计算过程中不参考测井解释结果，不受人为因素影响。

实施例六

基于前述的实施例，本申请实施例再提供一种含气性检测方法，包括：

S61：收集叠前角度域道集，确定有效入射角范围以及确定叠前角度域道集有效的频带范围；

在一些实施例中，所述有效入射角范围为0-36；

在一些实施例中，所述确定叠前角度域道集有效的频带范围之前，所述方法还包括：

进行剩余时间校正；

S62：按照入射角范围进行角度叠加，获得n个角度叠加后的地震数据，对应的入射角叠加范围的中心角度分别为angle1，angle2，angle3…anglen；

在一些实施例中，由于之后需要对角度叠加的地震数据提取属性，所以需要保证地震数据的信噪比不能太低，因此在进行角度叠加的时候，叠加的数据体个数一般不大于5个，所述n不大于5；

S63：对n个角度叠加后的地震数据的每一道计算振幅频率谱，获得地震数据有效的频带范围内的振幅频率谱面积area和振幅频率谱中能量最大值max_amp；

S64：应用振幅频率谱面积area和入射角叠加范围的中心角度，获得地震数据频谱面积在入射角方向的变化率数据体area_rate；

在一些实施例中，所述应用振幅频率谱面积area和入射角叠加范围的中心角度，获得地震数据频谱面积在入射角方向的变化率数据体area_rate的具体方法包括：

令Y轴为n个振幅频率谱面积area，X轴为n个入射角叠加范围的中心角度angle1，angle2，angle3…anglen，通过曲线拟合，获得地震数据频谱面积在入射角方向的变化率数据体area_rate；

S65：应用振幅频率谱中能量最大值max_amp和入射角叠加范围的中心角度，获得振幅频率谱中能量最大值在入射角方向的变化率数据体max_amp_rate；

在一些实施例中，所述应用振幅频率谱中能量最大值max_amp和入射角叠加范围的中心角度，获得振幅频率谱中能量最大值在入射角方向的变化率数据体max_amp_rate的具体方法包括：

令Y轴为n个振幅频率谱中能量最大值max_amp，X轴为n个入射角叠加范围的中心角度angle1，angle2，angle3…anglen，通过曲线拟合，获得振幅频率谱中能量最大值在入射角方向的变化率数据体max_amp_rate；

数据体area_rate代表的意义：由于受到储层含气的影响，地震数据的能量产生衰减，这种能量衰减的变化在振幅频率谱中体现为，随着入射角的增加能量衰减越强，即随着入射角的增加叠加数据频谱面积逐渐变小，受到含气性的影响地震数据的总能量随入射角增加减小，这里只考虑地震数据能量衰减与含气有关并未考虑储层厚度；数据体max_amp_rate代表的意义：由于受到储层含气的影响，地震数据的能量产生衰减，这种能量衰减的变化在振幅频率谱中体现为，随着入射角的增加能量衰减越强，即随着入射角的增加叠加数据振幅频率谱中能量最大值逐渐变小，受到含气的影响地震数据的峰值能量随入射角增加减小；

S66：将数据体area_rate和数据体max_amp_rate与实际工区已知产气井进行标定，完成含气性预测。

本申请提供的一种含气性检测方法，

(1)本发明完全基于地震数据进行计算，获得的结果确定性强；

(2)本发明在计算过程中不参考测井解释结果，不受人为因素影响。

实施例七

基于前述的实施例，本申请实施例再提供一种含气性检测方法，包括：

S71：收集叠前角度域道集，确定有效入射角范围以及确定叠前角度域道集有效的频带范围；

在一些实施例中，所述有效入射角范围为0-36；

在一些实施例中，所述叠前角度域道集有效的频带范围为7-75Hz；

在一些实施例中，所述确定叠前角度域道集有效的频带范围之前，所述方法还包括：

进行剩余时间校正；

S72：按照入射角范围进行角度叠加，获得n个角度叠加后的地震数据，对应的入射角叠加范围的中心角度分别为angle1，angle2，angle3…anglen；

在一些实施例中，由于之后需要对角度叠加的地震数据提取属性，所以需要保证地震数据的信噪比不能太低，因此在进行角度叠加的时候，叠加的数据体个数一般不大于5个，所述n不大于5；

S73：对n个角度叠加后的地震数据的每一道计算振幅频率谱，获得地震数据有效的频带范围内的振幅频率谱面积area和振幅频率谱中能量最大值max_amp；

S74：应用振幅频率谱面积area和入射角叠加范围的中心角度，获得地震数据频谱面积在入射角方向的变化率数据体area_rate；

在一些实施例中，所述应用振幅频率谱面积area和入射角叠加范围的中心角度，获得地震数据频谱面积在入射角方向的变化率数据体area_rate的具体方法包括：

令Y轴为n个振幅频率谱面积area，X轴为n个入射角叠加范围的中心角度angle1，angle2，angle3…anglen，通过曲线拟合，获得地震数据频谱面积在入射角方向的变化率数据体area_rate；

S75：应用振幅频率谱中能量最大值max_amp和入射角叠加范围的中心角度，获得振幅频率谱中能量最大值在入射角方向的变化率数据体max_amp_rate；

在一些实施例中，所述应用振幅频率谱中能量最大值max_amp和入射角叠加范围的中心角度，获得振幅频率谱中能量最大值在入射角方向的变化率数据体max_amp_rate的具体方法包括：

令Y轴为n个振幅频率谱中能量最大值max_amp，X轴为n个入射角叠加范围的中心角度angle1，angle2，angle3…anglen，通过曲线拟合，获得振幅频率谱中能量最大值在入射角方向的变化率数据体max_amp_rate；

数据体area_rate代表的意义：由于受到储层含气的影响，地震数据的能量产生衰减，这种能量衰减的变化在振幅频率谱中体现为，随着入射角的增加能量衰减越强，即随着入射角的增加叠加数据频谱面积逐渐变小，受到含气性的影响地震数据的总能量随入射角增加减小，这里只考虑地震数据能量衰减与含气有关并未考虑储层厚度；数据体max_amp_rate代表的意义：由于受到储层含气的影响，地震数据的能量产生衰减，这种能量衰减的变化在振幅频率谱中体现为，随着入射角的增加能量衰减越强，即随着入射角的增加叠加数据振幅频率谱中能量最大值逐渐变小，受到含气的影响地震数据的峰值能量随入射角增加减小；

S76：将数据体area_rate和数据体max_amp_rate与实际工区已知产气井进行标定，完成含气性预测。

本申请提供的一种含气性检测方法，

(1)本发明完全基于地震数据进行计算，获得的结果确定性强；

(2)本发明在计算过程中不参考测井解释结果，不受人为因素影响。

实施例八

基于实施例七方法，如图1和图6所示本申请根据真实数据给出的实施例：

S81：收集叠前角度域道集，确定有效入射角范围以及确定叠前角度域道集有效的频带范围；

在一些实施例中，如图2所示，所述有效入射角范围为0-36；

在一些实施例中，如图3所示，所述叠前角度域道集有效的频带范围为7-75Hz；

在一些实施例中，所述确定叠前角度域道集有效的频带范围之前，所述方法还包括：

进行剩余时间校正；

S82：按照入射角范围进行角度叠加，获得n个角度叠加后的地震数据，对应的入射角叠加范围的中心角度分别为angle1，angle2，angle3…anglen；

在一些实施例中，由于之后需要对角度叠加的地震数据提取属性，所以需要保证地震数据的信噪比不能太低，因此在进行角度叠加的时候，叠加的数据体个数一般不大于5个，所述n不大于5；

在一些实施例中，按照入射角范围进行角度叠加，获得3个角度叠加后的地震数据，每个角度叠加数据的入射角范围分别为0-12,12-24,24-36，3个角度叠加后的地震数据按照中心角可以表示为，seis_6,seis_18,seis_30；

S83：对n个角度叠加后的地震数据的每一道计算振幅频率谱，获得地震数据有效的频带范围内的振幅频率谱面积area和振幅频率谱中能量最大值max_amp；

在一些实施例中，对3个角度叠加后的地震数据每一道计算振幅频率谱，获得获得地震数据有效频带内的振幅频率谱面积area和振幅频率谱中能量最大值max_amp，area_6,area_18,area_30这些数据体代表的是不同入射角范围内的地震数据振幅频率谱面积信息，max_amp_6,max_amp_18,max_amp_30这些数据体代表的是不同入射角范围内的地震数据振幅频率谱中能量最大值的信息；

S84：应用振幅频率谱面积area和入射角叠加范围的中心角度，获得地震数据频谱面积在入射角方向的变化率数据体area_rate；

在一些实施例中，所述应用振幅频率谱面积area和入射角叠加范围的中心角度，获得地震数据频谱面积在入射角方向的变化率数据体area_rate的具体方法包括：

令Y轴为n个振幅频率谱面积area，X轴为n个入射角叠加范围的中心角度angle1，angle2，angle3…anglen，通过曲线拟合，获得地震数据频谱面积在入射角方向的变化率数据体area_rate；

在一些实施例中，令Y轴为中心角6、18、30的振幅频率谱面积area_fre，x轴为6、18、30，通过线性回归的方法计算出对应的斜率值，形成随着入射角增加的振幅频率谱面积的变化斜率体area_rate

S85：应用振幅频率谱中能量最大值max_amp和入射角叠加范围的中心角度，获得振幅频率谱中能量最大值在入射角方向的变化率数据体max_amp_rate；

在一些实施例中，所述应用振幅频率谱中能量最大值max_amp和入射角叠加范围的中心角度，获得振幅频率谱中能量最大值在入射角方向的变化率数据体max_amp_rate的具体方法包括：

令Y轴为n个振幅频率谱中能量最大值max_amp，X轴为n个入射角叠加范围的中心角度angle1，angle2，angle3…anglen，通过曲线拟合，获得振幅频率谱中能量最大值在入射角方向的变化率数据体max_amp_rate；

在一些实施例中，根据步骤5令Y轴为中心角6、18、30的振幅频率谱中能量最大值max_fre，x轴为6、18、30，通过线性回归的方法计算出对应的斜率值，形成随着入射角增加的振幅频率谱中能量最大值的变化斜率体max_amp_rate；

数据体area_rate代表的意义：由于受到储层含气的影响，地震数据的能量产生衰减，这种能量衰减的变化在振幅频率谱中体现为，随着入射角的增加能量衰减越强，即随着入射角的增加叠加数据频谱面积逐渐变小，受到含气性的影响地震数据的总能量随入射角增加减小，这里只考虑地震数据能量衰减与含气有关并未考虑储层厚度；数据体max_amp_rate代表的意义：由于受到储层含气的影响，地震数据的能量产生衰减，这种能量衰减的变化在振幅频率谱中体现为，随着入射角的增加能量衰减越强，即随着入射角的增加叠加数据振幅频率谱中能量最大值逐渐变小，受到含气的影响地震数据的峰值能量随入射角增加减小；

S86：将数据体area_rate和数据体max_amp_rate与实际工区已知产气井进行标定，完成含气性预测；

在一些实施例中，数据体area_rate、数据体max_amp_rate与实际工区已知产气井进行标定。工区内已知3口井，A井为产气井(含气饱和度70％)、B井为产水井(含气饱和度40％)，C井为特高产气井，D井为差产气井。A、B井作为标定井，C、D井作为验证井；

图4为数据体area_rate，图中A、B、C、D四口井的值分别为，-0.14、-0.04、-0.2和-0.1，将A井为产气井(含气饱和度70％)、B井为产水井(含气饱和度40％)作为标定井，因此可以确定C井含气饱和度明显高于70％，D井含气饱和度介于40-70％之间，这与实际井的产量完全吻合。

图5为数据体max_amp_rate，图中A、B、C、D四口井的值分别为，-0.1、-0.04、-0.14、-0.08，将A井为产气井(含气饱和度70％)、B井为产水井(含气饱和度40％)作为标定井，因此可以确定C井含气饱和度明显高于70％，D井含气饱和度介于40-70％之间，这与实际井的产量完全吻合。该方法可以有效得进行含气性预测研究。

本申请提供的一种含气性检测方法，

(1)本发明完全基于地震数据进行计算，获得的结果确定性强；

(2)本发明在计算过程中不参考测井解释结果，不受人为因素影响。

实施例九

基于前述的实施例，本申请实施例提供一种含气性检测装置，包括：

叠前角度域道集收集模块、中心角度计算模块、振幅频率谱数据计算模块、数据体area_rate计算模块、数据体max_amp_rate计算模块和含气性预测模块；

叠前角度域道集收集模块：收集叠前角度域道集，确定有效入射角范围以及确定叠前角度域道集有效的频带范围；

中心角度计算模块：按照入射角范围进行角度叠加，获得n个角度叠加后的地震数据，对应的入射角叠加范围的中心角度分别为angle1，angle2，angle3…anglen；

振幅频率谱数据计算模块：对n个角度叠加后的地震数据的每一道计算振幅频率谱，获得地震数据有效的频带范围内的振幅频率谱面积area和振幅频率谱中能量最大值max_amp；其中，所述确定叠前角度域道集有效的频带范围之前，还包括：进行剩余时间校正，所述n不大于5。

数据体area_rate计算模块：应用振幅频率谱面积area和入射角叠加范围的中心角度，获得地震数据频谱面积在入射角方向的变化率数据体area_rate；

所述应用振幅频率谱面积area和入射角叠加范围的中心角度，获得地震数据频谱面积在入射角方向的变化率数据体area_rate的具体方法包括：

令Y轴为n个振幅频率谱面积area，X轴为n个入射角叠加范围的中心角度angle1，angle2，angle3…anglen，通过曲线拟合，获得地震数据频谱面积在入射角方向的变化率数据体area_rate。

数据体max_amp_rate计算模块：应用振幅频率谱中能量最大值max_amp和入射角叠加范围的中心角度，获得振幅频率谱中能量最大值在入射角方向的变化率数据体max_amp_rate；

所述应用振幅频率谱中能量最大值max_amp和入射角叠加范围的中心角度，获得振幅频率谱中能量最大值在入射角方向的变化率数据体max_amp_rate的具体方法包括：

令Y轴为n个振幅频率谱中能量最大值max_amp，X轴为n个入射角叠加范围的中心角度angle1，angle2，angle3…anglen，通过曲线拟合，获得振幅频率谱中能量最大值在入射角方向的变化率数据体max_amp_rate。

含气性预测模块：将数据体area_rate和数据体max_amp_rate与实际工区已知产气井进行标定，完成含气性预测。

所述有效入射角范围为0-36。

所述叠前角度域道集有效的频带范围为7-75Hz。

需要说明的是，本申请实施例中，如果以软件功能模块的形式实现上述的含气性检测方法，并作为独立的产品销售或使用时，也可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请实施例的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机、服务器、或者网络设备等)执行本申请各个实施例所述方法的全部或部分。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read Only Memory)、磁碟或者光盘等各种可以存储程序代码的介质。这样，本申请实施例不限制于任何特定的硬件和软件结合。

相应地，本申请实施例提供一种存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现上述实施例中提供的含气性检测方法中的步骤。

实施例十

本申请实施例提供一种含气性检测设备存储器和处理器，所述存储器上存储有计算机程序，该计算机程序被所述处理器执行时，所述处理器配置为执行存储器中存储的多尺度电磁场分量去噪方法的程序，以实现以上述实施例提供的多尺度电磁场分量去噪方法中的步骤。

以上显示设备和存储介质实施例的描述，与上述方法实施例的描述是类似的，具有同方法实施例相似的有益效果。对于本申请计算机设备和存储介质实施例中未披露的技术细节，请参照本申请方法实施例的描述而理解。

这里需要指出的是：以上存储介质和设备实施例的描述，与上述方法实施例的描述是类似的，具有同方法实施例相似的有益效果。对于本申请存储介质和设备实施例中未披露的技术细节，请参照本申请方法实施例的描述而理解。

应理解，说明书通篇中提到的“一个实施例”或“一实施例”意味着与实施例有关的特定特征、结构或特性包括在本申请的至少一个实施例中。因此，在整个说明书各处出现的“在一个实施例中”或“在一实施例中”未必一定指相同的实施例。此外，这些特定的特征、结构或特性可以任意适合的方式结合在一个或多个实施例中。应理解，在本申请的各种实施例中，上述各过程的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本申请实施例的实施过程构成任何限定。上述本申请实施例序号仅仅为了描述，不代表实施例的优劣。

需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者装置不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者装置所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括该要素的过程、方法、物品或者装置中还存在另外的相同要素。

在本申请所提供的几个实施例中，应该理解到，所揭露的设备和方法，可以通过其它的方式实现。以上所描述的设备实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，如：多个单元或组件可以结合，或可以集成到另一个系统，或一些特征可以忽略，或不执行。另外，所显示或讨论的各组成部分相互之间的耦合、或直接耦合、或通信连接可以是通过一些接口，设备或单元的间接耦合或通信连接，可以是电性的、机械的或其它形式的。

上述作为分离部件说明的单元可以是、或也可以不是物理上分开的，作为单元显示的部件可以是、或也可以不是物理单元；既可以位于一个地方，也可以分布到多个网络单元上；可以根据实际的需要选择其中的部分或全部单元来实现本实施例方案的目的。

另外，在本申请各实施例中的各功能单元可以全部集成在一个处理单元中，也可以是各单元分别单独作为一个单元，也可以两个或两个以上单元集成在一个单元中；上述集成的单元既可以采用硬件的形式实现，也可以采用硬件加软件功能单元的形式实现。

本领域普通技术人员可以理解：实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成，前述的程序可以存储于计算机可读取存储介质中，该程序在执行时，执行包括上述方法实施例的步骤；而前述的存储介质包括：移动存储设备、只读存储器(ROM，Read Only Memory)、磁碟或者光盘等各种可以存储程序代码的介质。

或者，本申请上述集成的单元如果以软件功能模块的形式实现并作为独立的产品销售或使用时，也可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请实施例的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台控制器执行本申请各个实施例所述方法的全部或部分。而前述的存储介质包括：移动存储设备、ROM、磁碟或者光盘等各种可以存储程序代码的介质。

以上所述，仅为本申请的实施方式，但本申请的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本申请揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本申请的保护范围之内。因此，本申请的保护范围应以所述权利要求的保护范围为准。