叠前道集能量均衡处理方法、装置、设备及存储介质

技术领域

本申请涉及地球物理勘探领域，特别地涉及一种叠前道集能量均衡处理方法、装置、设备及存储介质。

背景技术

叠前CRP道集的质量好坏直接影响叠前弹性反演的信噪比和精度。如果叠前CRP道集并不能满足叠前弹性反演的要求，即使运用先进的反演算法，也会影响最终的反演效果。因此，在开展叠前弹性反演之前，必须对叠前CRP道集开展资料品质评价和保幅去噪处理，为后续叠前弹性反演提供高质量的地震资料基础。

通常的叠前道集优化处理步骤通常包括一下几个方面：①去噪，主要目的是为了提高叠前道集的信噪比。②转角度道集，主要是基于速度场将偏移距域叠前CRP道集转换到角度域用于叠前反演。③剩余时差校正，主要是为了将道集同相轴拉平，消除反演过程的虚假AVO现象。④超道集处理，目的同去噪类似，通过面元叠加的方式实现信噪比的提高。⑤角度部分叠加，通过将单个CDP角度道集约40道的数据，叠加成3-5角度的部分叠加数据，为叠前反演做输入，通过叠加可以提高叠前反演的稳定性。有的时候也会针对采集观测系统不规则的情况，针对性的采取一些补偿的方法。

但很少有人注意到，在高陡构造或者其他复杂构造区，叠前道集的能量存在系统性差异。这种能量差异通常是在地震资料处理阶段产生的。针对高陡构造的地震资料处理通常会在地震数据叠加之后进行射线补偿，从而达到对叠后地震数据能量均衡的处理，但叠前CRP道集上仍然存在由于构造形态导致的能量差异性。这种由于构造导致的能量差异性如果不进行均衡处理，带入到叠前反演中，会最终导致我们对储层的认识发生偏差。

发明内容

针对上述问题，本申请提供一种叠前道集能量均衡处理方法、装置、设备及存储介质。

本申请提供了一种叠前道集能量均衡处理方法，包括：

S1：数据准备：收集已知复杂构造区的叠后地震数据、叠前CRP道集和标志性地质层位；

S2：叠前CRP道集构造能量校正因子计算：应用所述叠后地震数据、叠前CRP道集和标志性地质层位计算叠前CRP道集构造能量校正因子；

S3：校正后的叠前CRP道集计算：应用叠前CRP道集构造能量校正因子和叠前CRP道集，得到校正后的叠前CRP道集。

在一些实施例中，所述应用所述叠后地震数据、叠前CRP道集和标志性地质层位计算叠前CRP道集构造能量校正因子的具体方法包括：

(1)提取叠后地震数据的大时窗均方根振幅属性；

(2)提取叠前CRP道集的大时窗均方根振幅属性；

(3)应用所述叠后地震数据的大时窗均方根振幅属性和叠前CRP道集的大时窗均方根振幅属性，得到叠前CRP道集构造能量校正因子。

在一些实施例中，所述提取叠后地震数据的大时窗均方根振幅属性的具体方法包括：

基于标志性地质层位，提取叠后地震数据的大时窗均方根振幅属性。

在一些实施例中，所述提取叠前CRP道集的大时窗均方根振幅属性的具体方法包括：

(1)对叠前CRP道集进行全叠加，得到基于叠前CRP道集直接叠加的叠后数据；

(2)基于标志性地质层位，对基于叠前CRP道集直接叠加的叠后数据，提取大时窗均方根振幅属性，得到叠前CRP道集的大时窗均方根振幅属性。

在一些实施例中，所述大时窗的选择原则为，时窗不小于200ms，且在全工区范围内，所述时窗内均为有效数据。

在一些实施例中，所述应用所述叠后地震数据的大时窗均方根振幅属性和叠前CRP道集的大时窗均方根振幅属性，得到叠前CRP道集构造能量校正因子的具体方法包括：

叠前CRP道集构造能量校正因子＝叠后地震数据的大时窗均方根振幅属性/叠前CRP道集的大时窗均方根振幅属性。

在一些实施例中，所述应用叠前CRP道集构造能量校正因子和叠前CRP道集，得到校正后的叠前CRP道集的具体方法包括：

校正后的叠前CRP道集＝叠前CRP道集构造能量校正因子×叠前CRP道集。

本申请实施例提供一种叠前道集能量均衡处理装置，包括：

数据准备模块、叠前CRP道集构造能量校正因子计算模块和校正后的叠前CRP道集计算模块；

所述数据准备模块：收集已知复杂构造区的叠后地震数据、叠前CRP道集和标志性地质层位；

所述叠前CRP道集构造能量校正因子计算模块：应用所述叠后地震数据、叠前CRP道集和标志性地质层位计算叠前CRP道集构造能量校正因子；

所述叠前CRP道集和标志性地质层位：应用叠前CRP道集构造能量校正因子和叠前CRP道集，得到校正后的叠前CRP道集。

本申请实施例提供一种叠前道集能量均衡处理设备，包括存储器和处理器，所述存储器上存储有计算机程序，该计算机程序被所述处理器执行时，执行上述任意一项所述叠前道集能量均衡处理方法。

本申请实施例提供一种存储介质，该存储介质存储的计算机程序，能够被一个或多个处理器执行，能够用来实现上述任一项所述叠前道集能量均衡处理方法。

本申请提供的一种叠前道集能量均衡处理方法、装置、设备及存储介质，

针对复杂构造、高陡构造地质条件下，叠前CRP道集存在与构造相关的系统性能量差异，在不需要重处理的基础上，提出的一种能够实现叠前道集能量校正的方法，从而提高叠前地震反演的精度，进一步更加准确的反映储层真实的弹性、物性及流体性质。

附图说明

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

图1为本申请实施例提供的一种叠前道集能量均衡处理方法的实现流程示意图；

图2为本申请实施例提供的标志层五峰-龙马溪组层图；

图3为本申请实施例提供的短轴方向基于叠前CRP道集直接叠加的剖面图；

图4为本申请实施例提供的叠前CRP道集叠加数据基于标志性地质层位的大时窗均方根振幅属性图；

图5为本申请实施例提供的叠后成果数据基于标志性地质层位的大时窗均方根振幅属性图；

图6为本申请实施例提供的叠前道集能量校正因子平面属性图；

图7为本申请实施例提供的能量校正后的叠前道集叠加数据基于标志性地质层位的大时窗均方根振幅属性图；

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

具体实施方式

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

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

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

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

在介绍本申请实施例提供的一种叠前道集能量均衡处理方法之前，对相关技术中存在的问题进行简单介绍：

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

实施例一

本申请实施例提供一种叠前道集能量均衡处理方法，图1为本申请实施例提供的一种叠前道集能量均衡处理方法的实现流程示意图，如图1所示，包括：

S1：数据准备：收集已知复杂构造区的叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)；

其中，x,y为地震观测系统的纵横坐标，t为地震数据纵向时间域采样，a记录了叠前CRP道集的偏移距信息或者入射角信息；

S2：叠前CRP道集构造能量校正因子λ(x,y)计算：应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)；

S3：校正后的叠前CRP道集SEIS_CRP_out(x,y,a,t)计算：应用叠前CRP道集构造能量校正因子λ(x,y)和叠前CRP道集SEIS_CRP(x,y,a,t)，得到校正后的叠前CRP道集SEIS_CRP_out(x,y,a,t)。

本申请提供的一种叠前道集能量均衡处理方法，针对复杂构造、高陡构造地质条件下，叠前CRP道集存在与构造相关的系统性能量差异，在不需要重处理的基础上，提出的一种能够实现叠前道集能量校正的方法，从而提高叠前地震反演的精度，进一步更加准确的反映储层真实的弹性、物性及流体性质。

实施例二

基于前述的实施例，本申请实施例再提供一种叠前道集能量均衡处理方法，包括：

S1：数据准备：收集已知复杂构造区的叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)；

其中，x,y为地震观测系统的纵横坐标，t为地震数据纵向时间域采样，a记录了叠前CRP道集的偏移距信息或者入射角信息；

S2：叠前CRP道集构造能量校正因子λ(x,y)计算：应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)；

在一些实施例中，所述应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)的具体方法包括：

(1)提取叠后地震数据SEIS_STK(x,y,t)的大时窗均方根振幅属性RMS_STK(x,y)；

(2)提取叠前CRP道集SEIS_CRP(x,y,a,t)的大时窗均方根振幅属性RMS_CDPstk(x,y)；

(3)应用所述叠后地震数据的大时窗均方根振幅属性RMS_STK(x,y)和叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)，得到叠前CRP道集构造能量校正因子λ(x,y)；

S3：校正后的叠前CRP道集SEIS_CRP_out(x,y,a,t)计算：应用叠前CRP道集构造能量校正因子λ(x,y)和叠前CRP道集SEIS_CRP(x,y,a,t)，得到校正后的叠前CRP道集SEIS_CRP_out(x,y,a,t)。

本申请提供的一种叠前道集能量均衡处理方法，针对复杂构造、高陡构造地质条件下，叠前CRP道集存在与构造相关的系统性能量差异，在不需要重处理的基础上，提出的一种能够实现叠前道集能量校正的方法，从而提高叠前地震反演的精度，进一步更加准确的反映储层真实的弹性、物性及流体。

实施例三

基于前述的实施例，本申请实施例再提供一种叠前道集能量均衡处理方法，包括：

S1：数据准备：收集已知复杂构造区的叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)；

其中，x,y为地震观测系统的纵横坐标，t为地震数据纵向时间域采样，a记录了叠前CRP道集的偏移距信息或者入射角信息；

S2：叠前CRP道集构造能量校正因子λ(x,y)计算：应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)；

在一些实施例中，所述应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)的具体方法包括：

(1)提取叠后地震数据SEIS_STK(x,y,t)的大时窗均方根振幅属性RMS_STK(x,y)；

(2)提取叠前CRP道集SEIS_CRP(x,y,a,t)的大时窗均方根振幅属性RMS_CDPstk(x,y)；

(3)应用所述叠后地震数据的大时窗均方根振幅属性RMS_STK(x,y)和叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)，得到叠前CRP道集构造能量校正因子λ(x,y)；

在一些实施例中，所述提取叠后地震数据的大时窗均方根振幅属性的具体方法包括：

基于标志性地质层位H(x,y)，提取叠后地震数据的大时窗均方根振幅属性RMS_STK(x,y)；

S3：校正后的叠前CRP道集SEIS_CRP_out(x,y,a,t)计算：应用叠前CRP道集构造能量校正因子λ(x,y)和叠前CRP道集SEIS_CRP(x,y,a,t)，得到校正后的叠前CRP道集SEIS_CRP_out(x,y,a,t)。

本申请提供的一种叠前道集能量均衡处理方法，针对复杂构造、高陡构造地质条件下，叠前CRP道集存在与构造相关的系统性能量差异，在不需要重处理的基础上，提出的一种能够实现叠前道集能量校正的方法，从而提高叠前地震反演的精度，进一步更加准确的反映储层真实的弹性、物性及流体。

实施例四

基于前述的实施例，本申请实施例再提供一种叠前道集能量均衡处理方法，包括：

S1：数据准备：收集已知复杂构造区的叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)；

其中，x,y为地震观测系统的纵横坐标，t为地震数据纵向时间域采样，a记录了叠前CRP道集的偏移距信息或者入射角信息；

S2：叠前CRP道集构造能量校正因子λ(x,y)计算：应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)；

在一些实施例中，所述应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)的具体方法包括：

(1)提取叠后地震数据SEIS_STK(x,y,t)的大时窗均方根振幅属性RMS_STK(x,y)；

(2)提取叠前CRP道集SEIS_CRP(x,y,a,t)的大时窗均方根振幅属性RMS_CDPstk(x,y)；

(3)应用所述叠后地震数据的大时窗均方根振幅属性RMS_STK(x,y)和叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)，得到叠前CRP道集构造能量校正因子λ(x,y)；

在一些实施例中，所述提取叠后地震数据的大时窗均方根振幅属性的具体方法包括：

基于标志性地质层位H(x,y)，提取叠后地震数据的大时窗均方根振幅属性RMS_STK(x,y)；

在一些实施例中，所述提取叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)的具体方法包括：

(1)对叠前CRP道集SEIS_CRP(x,y,a,t)进行全叠加，得到基于叠前CRP道集直接叠加的叠后数据SEIS_CRPstk(x,y,t)；

(2)基于标志性地质层位H(x,y)，对基于叠前CRP道集直接叠加的叠后数据SEIS_CRPstk(x,y,t)，提取大时窗均方根振幅属性，得到叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)；

S3：校正后的叠前CRP道集SEIS_CRP_out(x,y,a,t)计算：应用叠前CRP道集构造能量校正因子λ(x,y)和叠前CRP道集SEIS_CRP(x,y,a,t)，得到校正后的叠前CRP道集SEIS_CRP_out(x,y,a,t)。

本申请提供的一种叠前道集能量均衡处理方法，针对复杂构造、高陡构造地质条件下，叠前CRP道集存在与构造相关的系统性能量差异，在不需要重处理的基础上，提出的一种能够实现叠前道集能量校正的方法，从而提高叠前地震反演的精度，进一步更加准确的反映储层真实的弹性、物性及流体。

实施例五

基于前述的实施例，本申请实施例再提供一种叠前道集能量均衡处理方法，包括：

S1：数据准备：收集已知复杂构造区的叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)；

其中，x,y为地震观测系统的纵横坐标，t为地震数据纵向时间域采样，a记录了叠前CRP道集的偏移距信息或者入射角信息；

S2：叠前CRP道集构造能量校正因子λ(x,y)计算：应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)；

在一些实施例中，所述应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)的具体方法包括：

(1)提取叠后地震数据SEIS_STK(x,y,t)的大时窗均方根振幅属性RMS_STK(x,y)；

(2)提取叠前CRP道集SEIS_CRP(x,y,a,t)的大时窗均方根振幅属性RMS_CDPstk(x,y)；

(3)应用所述叠后地震数据的大时窗均方根振幅属性RMS_STK(x,y)和叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)，得到叠前CRP道集构造能量校正因子λ(x,y)；

在一些实施例中，在一些实施例中，所述提取叠后地震数据的大时窗均方根振幅属性的具体方法包括：

基于标志性地质层位H(x,y)，提取叠后地震数据的大时窗均方根振幅属性RMS_STK(x,y)；

在一些实施例中，所述提取叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)的具体方法包括：

(1)对叠前CRP道集SEIS_CRP(x,y,a,t)进行全叠加，得到基于叠前CRP道集直接叠加的叠后数据SEIS_CRPstk(x,y,t)；

(2)基于标志性地质层位H(x,y)，对基于叠前CRP道集直接叠加的叠后数据SEIS_CRPstk(x,y,t)，提取大时窗均方根振幅属性，得到叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)；

在一些实施例中，所述大时窗的选择原则为，时窗不小于200ms，且在全工区范围内，所述时窗内均为有效数据；

S3：校正后的叠前CRP道集SEIS_CRP_out(x,y,a,t)计算：应用叠前CRP道集构造能量校正因子λ(x,y)和叠前CRP道集SEIS_CRP(x,y,a,t)，得到校正后的叠前CRP道集SEIS_CRP_out(x,y,a,t)。

本申请提供的一种叠前道集能量均衡处理方法，针对复杂构造、高陡构造地质条件下，叠前CRP道集存在与构造相关的系统性能量差异，在不需要重处理的基础上，提出的一种能够实现叠前道集能量校正的方法，从而提高叠前地震反演的精度，进一步更加准确的反映储层真实的弹性、物性及流体。

实施例六

基于前述的实施例，本申请实施例再提供一种叠前道集能量均衡处理方法，包括：

S1：数据准备：收集已知复杂构造区的叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)；

其中，x,y为地震观测系统的纵横坐标，t为地震数据纵向时间域采样，a记录了叠前CRP道集的偏移距信息或者入射角信息；

S2：叠前CRP道集构造能量校正因子λ(x,y)计算：应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)；

在一些实施例中，所述应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)的具体方法包括：

(1)提取叠后地震数据SEIS_STK(x,y,t)的大时窗均方根振幅属性RMS_STK(x,y)；

(2)提取叠前CRP道集SEIS_CRP(x,y,a,t)的大时窗均方根振幅属性RMS_CDPstk(x,y)；

(3)应用所述叠后地震数据的大时窗均方根振幅属性RMS_STK(x,y)和叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)，得到叠前CRP道集构造能量校正因子λ(x,y)；

在一些实施例中，所述提取叠后地震数据的大时窗均方根振幅属性的具体方法包括：

基于标志性地质层位H(x,y)，提取叠后地震数据的大时窗均方根振幅属性RMS_STK(x,y)；

在一些实施例中，所述提取叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)的具体方法包括：

(1)对叠前CRP道集SEIS_CRP(x,y,a,t)进行全叠加，得到基于叠前CRP道集直接叠加的叠后数据SEIS_CRPstk(x,y,t)；

(2)基于标志性地质层位H(x,y)，对基于叠前CRP道集直接叠加的叠后数据SEIS_CRPstk(x,y,t)，提取大时窗均方根振幅属性，得到叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)；

在一些实施例中，所述大时窗的选择原则为，时窗不小于200ms，且在全工区范围内，所述时窗内均为有效数据；

在一些实施例中，所述应用所述叠后地震数据的大时窗均方根振幅属性和叠前CRP道集的大时窗均方根振幅属性，得到叠前CRP道集构造能量校正因子的具体方法包括：

叠前CRP道集构造能量校正因子＝叠后地震数据的大时窗均方根振幅属性/叠前CRP道集的大时窗均方根振幅属性，即，

SEIS_CRP_out(x,y,a,t)＝λ(x,y)·SEIS_CRP(x,y,a,t)；

S3：校正后的叠前CRP道集SEIS_CRP_out(x,y,a,t)计算：应用叠前CRP道集构造能量校正因子λ(x,y)和叠前CRP道集SEIS_CRP(x,y,a,t)，得到校正后的叠前CRP道集SEIS_CRP_out(x,y,a,t)。

本申请提供的一种叠前道集能量均衡处理方法，针对复杂构造、高陡构造地质条件下，叠前CRP道集存在与构造相关的系统性能量差异，在不需要重处理的基础上，提出的一种能够实现叠前道集能量校正的方法，从而提高叠前地震反演的精度，进一步更加准确的反映储层真实的弹性、物性及流体。

实施例七

基于前述的实施例，本申请实施例再提供一种叠前道集能量均衡处理方法，包括：

S1：数据准备：收集已知复杂构造区的叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)；

其中，x,y为地震观测系统的纵横坐标，t为地震数据纵向时间域采样，a记录了叠前CRP道集的偏移距信息或者入射角信息；

S2：叠前CRP道集构造能量校正因子λ(x,y)计算：应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)；

在一些实施例中，所述应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)的具体方法包括：

(1)提取叠后地震数据SEIS_STK(x,y,t)的大时窗均方根振幅属性RMS_STK(x,y)；

(2)提取叠前CRP道集SEIS_CRP(x,y,a,t)的大时窗均方根振幅属性RMS_CDPstk(x,y)；

(3)应用所述叠后地震数据的大时窗均方根振幅属性RMS_STK(x,y)和叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)，得到叠前CRP道集构造能量校正因子λ(x,y)；

在一些实施例中，所述提取叠后地震数据的大时窗均方根振幅属性的具体方法包括：

基于标志性地质层位H(x,y)，提取叠后地震数据的大时窗均方根振幅属性RMS_STK(x,y)；

在一些实施例中，所述提取叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)的具体方法包括：

(1)对叠前CRP道集SEIS_CRP(x,y,a,t)进行全叠加，得到基于叠前CRP道集直接叠加的叠后数据SEIS_CRPstk(x,y,t)；

(2)基于标志性地质层位H(x,y)，对基于叠前CRP道集直接叠加的叠后数据SEIS_CRPstk(x,y,t)，提取大时窗均方根振幅属性，得到叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)；

在一些实施例中，所述大时窗的选择原则为，时窗不小于200ms，且在全工区范围内，所述时窗内均为有效数据；

在一些实施例中，所述应用所述叠后地震数据的大时窗均方根振幅属性和叠前CRP道集的大时窗均方根振幅属性，得到叠前CRP道集构造能量校正因子的具体方法包括：

叠前CRP道集构造能量校正因子＝叠后地震数据的大时窗均方根振幅属性/叠前CRP道集的大时窗均方根振幅属性，即，

SEIS_CRP_out(x,y,a,t)＝λ(x,y)·SEIS_CRP(x,y,a,t)；

S3：校正后的叠前CRP道集SEIS_CRP_out(x,y,a,t)计算：应用叠前CRP道集构造能量校正因子λ(x,y)和叠前CRP道集SEIS_CRP(x,y,a,t)，得到校正后的叠前CRP道集SEIS_CRP_out(x,y,a,t)；

在一些实施例中，所述应用叠前CRP道集构造能量校正因子和叠前CRP道集，得到校正后的叠前CRP道集的具体方法包括：

校正后的叠前CRP道集＝叠前CRP道集构造能量校正因子×叠前CRP道集，即，SEIS_CRP_out(x,y,a,t)＝λ(x,y)·SEIS_CRP(x,y,a,t)。

对每个CDP点(x,y)的任意一道、任意时间(a,t)都加上增益λ(x,y)得到能量校正后的CRP道集SEIS_CRP_out(x,y,a,t)。

本申请提供的一种叠前道集能量均衡处理方法，针对复杂构造、高陡构造地质条件下，叠前CRP道集存在与构造相关的系统性能量差异，在不需要重处理的基础上，提出的一种能够实现叠前道集能量校正的方法，从而提高叠前地震反演的精度，进一步更加准确的反映储层真实的弹性、物性及流体。

实施例八

基于实施例七方法，本申请根据真实数据，给出中国南方某双复杂地质条件下非常规页岩气储层叠前道集能量一致性校正的实施例：

该工区目的层龙马溪组海相页岩构造起伏剧烈，埋深从500米至4000米变化.图2为龙马溪组时间域埋深图。图3为沿短轴方向基于CRP道集叠加得到的叠后地震剖面，可以清楚的看到再构造高陡部位能量变弱，剖面整体能量不均衡。图4为基于叠前CRP道集进行全叠加，并基于叠前CRP道集直接叠加的叠后数据基于标志层，提取相同时窗的大时窗均方根振幅属性。从图4可以清楚的看到平面振幅属性与标志层龙马溪组埋深存在很好的对应关系，说明能量的差异性是受复杂构造的影响。

图5为叠后成果数据基于标志性地质层位提取相同时窗的大时窗均方根振幅属性RMS_STK(x,y)，可以看到，叠后成果数据的能量更均衡、更合理，不受构造的影响。基于提取的RMS_CDPstk(x,y)属性和RMS_STK(x,y)属性，利用公式，计算得到能量校正因子λ(x,y)：

λ(x,y)＝RMS_STK(x,y)/RMS_CRPstk(x,y)；

图6为得到的能量校正因子平面属性。

基于公式SEIS_CRP_out(x,y,a,t)＝λ(x,y)·SEIS_CRP(x,y,a,t)，可以计算得到校正后的CRP道集为SEIS_CRP_out(x,y,a,t)。并对该道集进行叠加，基于标志层提取大时窗均方根振幅属性RMS_CRPOUTstk(x,y)，如图7所示。可以看到能量校正后的道集，能量更均衡、更合理，不受构造的影响。基于该道集进行预处理之后进行叠前弹性、物性反演，可以得到更加合理可靠的储层及流体信息。

实施例九

基于前述的实施例，本申请实施例提供一种叠前道集能量均衡处理装置，包括：

数据准备模块、叠前CRP道集构造能量校正因子计算模块和校正后的叠前CRP道集计算模块；

所述数据准备模块：收集已知复杂构造区的叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)；

其中，x,y为地震观测系统的纵横坐标，t为地震数据纵向时间域采样，a记录了叠前CRP道集的偏移距信息或者入射角信息；

所述叠前CRP道集构造能量校正因子计算模块：应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)；

所述校正后的叠前CRP道集计算模块：应用叠前CRP道集构造能量校正因子λ(x,y)和叠前CRP道集SEIS_CRP(x,y,a,t)，得到校正后的叠前CRP道集SEIS_CRP_out(x,y,a,t)。

本申请提供的一种叠前道集能量均衡处理装置，针对复杂构造、高陡构造地质条件下，叠前CRP道集存在与构造相关的系统性能量差异，在不需要重处理的基础上，提出的一种能够实现叠前道集能量校正的方法，从而提高叠前地震反演的精度，进一步更加准确的反映储层真实的弹性、物性及流体性质。

实施例十

基于前述的实施例，本申请实施例提供一种叠前道集能量均衡处理装置，包括：

数据准备模块、叠前CRP道集构造能量校正因子计算模块和校正后的叠前CRP道集计算模块；

所述数据准备模块：收集已知复杂构造区的叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)；

其中，x,y为地震观测系统的纵横坐标，t为地震数据纵向时间域采样，a记录了叠前CRP道集的偏移距信息或者入射角信息；

所述叠前CRP道集构造能量校正因子计算模块：应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)；

在一些实施例中，所述应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)的具体方法包括：

(1)提取叠后地震数据SEIS_STK(x,y,t)的大时窗均方根振幅属性RMS_STK(x,y)；

(2)提取叠前CRP道集SEIS_CRP(x,y,a,t)的大时窗均方根振幅属性RMS_CDPstk(x,y)；

(3)应用所述叠后地震数据的大时窗均方根振幅属性RMS_STK(x,y)和叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)，得到叠前CRP道集构造能量校正因子λ(x,y)；

所述校正后的叠前CRP道集计算模块：应用叠前CRP道集构造能量校正因子λ(x,y)和叠前CRP道集SEIS_CRP(x,y,a,t)，得到校正后的叠前CRP道集SEIS_CRP_out(x,y,a,t)。

本申请提供的一种叠前道集能量均衡处理装置，针对复杂构造、高陡构造地质条件下，叠前CRP道集存在与构造相关的系统性能量差异，在不需要重处理的基础上，提出的一种能够实现叠前道集能量校正的方法，从而提高叠前地震反演的精度，进一步更加准确的反映储层真实的弹性、物性及流体性质。

实施例十一

基于前述的实施例，本申请实施例提供一种叠前道集能量均衡处理装置，包括：

数据准备模块、叠前CRP道集构造能量校正因子计算模块和校正后的叠前CRP道集计算模块；

所述数据准备模块：收集已知复杂构造区的叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)；

其中，x,y为地震观测系统的纵横坐标，t为地震数据纵向时间域采样，a记录了叠前CRP道集的偏移距信息或者入射角信息；

所述叠前CRP道集构造能量校正因子计算模块：应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)；

在一些实施例中，所述应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)的具体方法包括：

(1)提取叠后地震数据SEIS_STK(x,y,t)的大时窗均方根振幅属性RMS_STK(x,y)；

(2)提取叠前CRP道集SEIS_CRP(x,y,a,t)的大时窗均方根振幅属性RMS_CDPstk(x,y)；

(3)应用所述叠后地震数据的大时窗均方根振幅属性RMS_STK(x,y)和叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)，得到叠前CRP道集构造能量校正因子λ(x,y)；

在一些实施例中，所述提取叠后地震数据的大时窗均方根振幅属性的具体方法包括：

基于标志性地质层位H(x,y)，提取叠后地震数据的大时窗均方根振幅属性RMS_STK(x,y)；

所述校正后的叠前CRP道集计算模块：应用叠前CRP道集构造能量校正因子λ(x,y)和叠前CRP道集SEIS_CRP(x,y,a,t)，得到校正后的叠前CRP道集SEIS_CRP_out(x,y,a,t)。

本申请提供的一种叠前道集能量均衡处理装置，针对复杂构造、高陡构造地质条件下，叠前CRP道集存在与构造相关的系统性能量差异，在不需要重处理的基础上，提出的一种能够实现叠前道集能量校正的方法，从而提高叠前地震反演的精度，进一步更加准确的反映储层真实的弹性、物性及流体性质。

实施例十二

基于前述的实施例，本申请实施例提供一种叠前道集能量均衡处理装置，包括：

数据准备模块、叠前CRP道集构造能量校正因子计算模块和校正后的叠前CRP道集计算模块；

所述数据准备模块：收集已知复杂构造区的叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)；

其中，x,y为地震观测系统的纵横坐标，t为地震数据纵向时间域采样，a记录了叠前CRP道集的偏移距信息或者入射角信息；

所述叠前CRP道集构造能量校正因子计算模块：应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)；

在一些实施例中，所述应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)的具体方法包括：

(1)提取叠后地震数据SEIS_STK(x,y,t)的大时窗均方根振幅属性RMS_STK(x,y)；

(2)提取叠前CRP道集SEIS_CRP(x,y,a,t)的大时窗均方根振幅属性RMS_CDPstk(x,y)；

(3)应用所述叠后地震数据的大时窗均方根振幅属性RMS_STK(x,y)和叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)，得到叠前CRP道集构造能量校正因子λ(x,y)；

在一些实施例中，所述提取叠后地震数据的大时窗均方根振幅属性的具体方法包括：

基于标志性地质层位H(x,y)，提取叠后地震数据的大时窗均方根振幅属性RMS_STK(x,y)；

在一些实施例中，所述提取叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)的具体方法包括：

(1)对叠前CRP道集SEIS_CRP(x,y,a,t)进行全叠加，得到基于叠前CRP道集直接叠加的叠后数据SEIS_CRPstk(x,y,t)；

(2)基于标志性地质层位H(x,y)，对基于叠前CRP道集直接叠加的叠后数据SEIS_CRPstk(x,y,t)，提取大时窗均方根振幅属性，得到叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)；

所述校正后的叠前CRP道集计算模块：应用叠前CRP道集构造能量校正因子λ(x,y)和叠前CRP道集SEIS_CRP(x,y,a,t)，得到校正后的叠前CRP道集SEIS_CRP_out(x,y,a,t)。

本申请提供的一种叠前道集能量均衡处理装置，针对复杂构造、高陡构造地质条件下，叠前CRP道集存在与构造相关的系统性能量差异，在不需要重处理的基础上，提出的一种能够实现叠前道集能量校正的方法，从而提高叠前地震反演的精度，进一步更加准确的反映储层真实的弹性、物性及流体性质。

实施例十三

基于前述的实施例，本申请实施例提供一种叠前道集能量均衡处理装置，包括：

数据准备模块、叠前CRP道集构造能量校正因子计算模块和校正后的叠前CRP道集计算模块；

所述数据准备模块：收集已知复杂构造区的叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)；

其中，x,y为地震观测系统的纵横坐标，t为地震数据纵向时间域采样，a记录了叠前CRP道集的偏移距信息或者入射角信息；

所述叠前CRP道集构造能量校正因子计算模块：应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)；

在一些实施例中，所述应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)的具体方法包括：

(1)提取叠后地震数据SEIS_STK(x,y,t)的大时窗均方根振幅属性RMS_STK(x,y)；

(2)提取叠前CRP道集SEIS_CRP(x,y,a,t)的大时窗均方根振幅属性RMS_CDPstk(x,y)；

(3)应用所述叠后地震数据的大时窗均方根振幅属性RMS_STK(x,y)和叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)，得到叠前CRP道集构造能量校正因子λ(x,y)；

在一些实施例中，所述提取叠后地震数据的大时窗均方根振幅属性的具体方法包括：

基于标志性地质层位H(x,y)，提取叠后地震数据的大时窗均方根振幅属性RMS_STK(x,y)；

在一些实施例中，所述提取叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)的具体方法包括：

(1)对叠前CRP道集SEIS_CRP(x,y,a,t)进行全叠加，得到基于叠前CRP道集直接叠加的叠后数据SEIS_CRPstk(x,y,t)；

(2)基于标志性地质层位H(x,y)，对基于叠前CRP道集直接叠加的叠后数据SEIS_CRPstk(x,y,t)，提取大时窗均方根振幅属性，得到叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)；

在一些实施例中，所述大时窗的选择原则为，时窗不小于200ms，且在全工区范围内，所述时窗内均为有效数据；

所述校正后的叠前CRP道集计算模块：应用叠前CRP道集构造能量校正因子λ(x,y)和叠前CRP道集SEIS_CRP(x,y,a,t)，得到校正后的叠前CRP道集SEIS_CRP_out(x,y,a,t)。

本申请提供的一种叠前道集能量均衡处理装置，针对复杂构造、高陡构造地质条件下，叠前CRP道集存在与构造相关的系统性能量差异，在不需要重处理的基础上，提出的一种能够实现叠前道集能量校正的方法，从而提高叠前地震反演的精度，进一步更加准确的反映储层真实的弹性、物性及流体性质。

实施例十四

基于前述的实施例，本申请实施例提供一种叠前道集能量均衡处理装置，包括：

数据准备模块、叠前CRP道集构造能量校正因子计算模块和校正后的叠前CRP道集计算模块；

所述数据准备模块：收集已知复杂构造区的叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)；

其中，x,y为地震观测系统的纵横坐标，t为地震数据纵向时间域采样，a记录了叠前CRP道集的偏移距信息或者入射角信息；

所述叠前CRP道集构造能量校正因子计算模块：应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)；

在一些实施例中，所述应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)的具体方法包括：

(1)提取叠后地震数据SEIS_STK(x,y,t)的大时窗均方根振幅属性RMS_STK(x,y)；

(2)提取叠前CRP道集SEIS_CRP(x,y,a,t)的大时窗均方根振幅属性RMS_CDPstk(x,y)；

(3)应用所述叠后地震数据的大时窗均方根振幅属性RMS_STK(x,y)和叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)，得到叠前CRP道集构造能量校正因子λ(x,y)；

在一些实施例中，所述提取叠后地震数据的大时窗均方根振幅属性的具体方法包括：

基于标志性地质层位H(x,y)，提取叠后地震数据的大时窗均方根振幅属性RMS_STK(x,y)；

在一些实施例中，所述提取叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)的具体方法包括：

(1)对叠前CRP道集SEIS_CRP(x,y,a,t)进行全叠加，得到基于叠前CRP道集直接叠加的叠后数据SEIS_CRPstk(x,y,t)；

(2)基于标志性地质层位H(x,y)，对基于叠前CRP道集直接叠加的叠后数据SEIS_CRPstk(x,y,t)，提取大时窗均方根振幅属性，得到叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)；

在一些实施例中，所述大时窗的选择原则为，时窗不小于200ms，且在全工区范围内，所述时窗内均为有效数据；

在一些实施例中，所述应用所述叠后地震数据的大时窗均方根振幅属性和叠前CRP道集的大时窗均方根振幅属性，得到叠前CRP道集构造能量校正因子的具体方法包括：

叠前CRP道集构造能量校正因子＝叠后地震数据的大时窗均方根振幅属性/叠前CRP道集的大时窗均方根振幅属性，即，

SEIS_CRP_out(x,y,a,t)＝λ(x,y)·SEIS_CRP(x,y,a,t)；

所述校正后的叠前CRP道集计算模块：应用叠前CRP道集构造能量校正因子λ(x,y)和叠前CRP道集SEIS_CRP(x,y,a,t)，得到校正后的叠前CRP道集SEIS_CRP_out(x,y,a,t)。

本申请提供的一种叠前道集能量均衡处理装置，针对复杂构造、高陡构造地质条件下，叠前CRP道集存在与构造相关的系统性能量差异，在不需要重处理的基础上，提出的一种能够实现叠前道集能量校正的方法，从而提高叠前地震反演的精度，进一步更加准确的反映储层真实的弹性、物性及流体性质。

实施例十五

基于前述的实施例，本申请实施例提供一种叠前道集能量均衡处理装置，包括：

数据准备模块、叠前CRP道集构造能量校正因子计算模块和校正后的叠前CRP道集计算模块；

所述数据准备模块：收集已知复杂构造区的叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)；

其中，x,y为地震观测系统的纵横坐标，t为地震数据纵向时间域采样，a记录了叠前CRP道集的偏移距信息或者入射角信息；

所述叠前CRP道集构造能量校正因子计算模块：应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)；

在一些实施例中，所述应用所述叠后地震数据SEIS_STK(x,y,t)、叠前CRP道集SEIS_CRP(x,y,a,t)和标志性地质层位H(x,y)计算叠前CRP道集构造能量校正因子λ(x,y)的具体方法包括：

(1)提取叠后地震数据SEIS_STK(x,y,t)的大时窗均方根振幅属性RMS_STK(x,y)；

(2)提取叠前CRP道集SEIS_CRP(x,y,a,t)的大时窗均方根振幅属性RMS_CDPstk(x,y)；

(3)应用所述叠后地震数据的大时窗均方根振幅属性RMS_STK(x,y)和叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)，得到叠前CRP道集构造能量校正因子λ(x,y)；

在一些实施例中，所述提取叠后地震数据的大时窗均方根振幅属性的具体方法包括：

基于标志性地质层位H(x,y)，提取叠后地震数据的大时窗均方根振幅属性RMS_STK(x,y)；

在一些实施例中，所述提取叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)的具体方法包括：

(1)对叠前CRP道集SEIS_CRP(x,y,a,t)进行全叠加，得到基于叠前CRP道集直接叠加的叠后数据SEIS_CRPstk(x,y,t)；

(2)基于标志性地质层位H(x,y)，对基于叠前CRP道集直接叠加的叠后数据SEIS_CRPstk(x,y,t)，提取大时窗均方根振幅属性，得到叠前CRP道集的大时窗均方根振幅属性RMS_CDPstk(x,y)；

在一些实施例中，所述大时窗的选择原则为，时窗不小于200ms，且在全工区范围内，所述时窗内均为有效数据；

在一些实施例中，所述应用所述叠后地震数据的大时窗均方根振幅属性和叠前CRP道集的大时窗均方根振幅属性，得到叠前CRP道集构造能量校正因子的具体方法包括：

叠前CRP道集构造能量校正因子＝叠后地震数据的大时窗均方根振幅属性/叠前CRP道集的大时窗均方根振幅属性，即，

SEIS_CRP_out(x,y,a,t)＝λ(x,y)·SEIS_CRP(x,y,a,t)；

所述校正后的叠前CRP道集计算模块：应用叠前CRP道集构造能量校正因子λ(x,y)和叠前CRP道集SEIS_CRP(x,y,a,t)，得到校正后的叠前CRP道集SEIS_CRP_out(x,y,a,t)；

在一些实施例中，所述应用叠前CRP道集构造能量校正因子和叠前CRP道集，得到校正后的叠前CRP道集的具体方法包括：

校正后的叠前CRP道集＝叠前CRP道集构造能量校正因子×叠前CRP道集，即，SEIS_CRP_out(x,y,a,t)＝λ(x,y)·SEIS_CRP(x,y,a,t)。

对每个CDP点(x,y)的任意一道、任意时间(a,t)都加上增益λ(x,y)得到能量校正后的CRP道集SEIS_CRP_out(x,y,a,t)。

本申请提供的一种叠前道集能量均衡处理装置，针对复杂构造、高陡构造地质条件下，叠前CRP道集存在与构造相关的系统性能量差异，在不需要重处理的基础上，提出的一种能够实现叠前道集能量校正的方法，从而提高叠前地震反演的精度，进一步更加准确的反映储层真实的弹性、物性及流体性质。

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

相应地，本申请实施例提供一种存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现上述实施例中提供的叠前道集能量均衡处理方法中的步骤。

实施例十六

本申请实施例提供一种叠前道集能量均衡处理设备包括存储器和处理器，所述存储器上存储有计算机程序，该计算机程序被所述处理器执行时，所述处理器配置为执行存储器中存储的叠前道集能量均衡处理方法的程序，以实现以上述实施例提供的叠前道集能量均衡处理方法中的步骤。

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

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

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

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

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

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

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

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

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

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