CN102998703B - 基于地表一致性反褶积进行储层预测的方法及设备 - Google Patents

Abstract

本发明提供一种基于地表一致性反褶积进行储层预测的方法及设备，所述的方法包括：采集来自地下的反射地震信号形成的地震记录；根据地震记录生成目标函数；对地震记录进行地表一致性分解，得到炮点分量和检波点分量；对炮点分量和检波点分量进行最小相位化处理，得到最小相位炮检响应；根据目标函数以及最小相位炮检响应得到消除非地质因素影响后的地震记录；根据消除非地质因素影响后的地震记录绘制地震剖面图像；对地震剖面图像进行综合解释，得到储层预测的结果。消除了采集和近地表因素对地震记录的影响，提高了地震资料反映真实地下结构的能力，进一步提高了储层预测的准确性。

Description

基于地表一致性反褶积进行储层预测的方法及设备

技术领域

本发明关于地球物理勘探技术领域，特别是关于油气勘探中地震资料的处理技术，具体的讲是一种基于地表一致性反褶积进行储层预测的方法及设备。

背景技术

地震勘探是一种利用人工地震技术探测地下结构的勘探方法，其按照一定的方式人工激发地震波，检波器接收来自地下的反射信号，通过对反射信号的处理和分析探测地下结构，进而根据探测到的地下结构进行油气的储层预测。

根据地震数据的反射特征对地下介质的构造形态和储层物性进行分析和预测是地震资料解释工作的主要研究内容。在实际地震勘探中，地震反射特征的变化除了受地下构造和岩性等地质因素的影响之外，还受到野外采集因素和地表因素的影响，如何在地震资料处理阶段，有效地消除地震采集和地表因素对地震反射特征的影响，凸显地下介质的真实反射特征，对提高地震勘探、储层预测精度具有十分重要的作用。

Taner于1981年提出利用地表一致性方法对近地表影响进行校正的基本思想，所述的地表一致性是指非地质因素（采集和地表因素）对地震记录的影响只与炮点和检波点的位置有关；共炮点地震道具有相同的炮点校正量，共检波点地震道具有相同的检波点校正量；校正量与地震波的路径和反射时间无关。Levin于1989年将地表一致性的思想应用于校正激发和接收因素对地震子波的影响，提出了地表一致性反褶积方法，该方法提出不久就在工业界取得了迅速的推广和应用，在地震资料处理和解释中发挥了重要作用。Cambois于1992年对数谱域实现了地震记录的分解，将地表一致性反褶积和其它地表一致性处理纳入统一的理论框架。Cary于1993年将地震记录分解为炮点、检波点、中心点和炮检距四个分量，提高了由地震记录估算炮点和检波点分量的精度。李生杰于2002年采用共轭梯度法在时间域对地震记录进行谱分解处理，改善了地表一致性的处理质量。

上述的研究工作大部分集中在地震记录的分解上，即如何由地震记录中更加准确地估算出非地质因素对地震子波振幅谱的影响。但是上述的各项研究涉及如何消除非地质因素对地震记录的影响时，基本上都是沿用预测反褶积的做法，即：将炮点分量和检波点分量褶积得到炮点和检波点的综合响应，然后对每个地震道的综合响应进行相同步长的预测反褶积。当采集和地表差异变化较大时，这种做法不能取得理想的处理效果。针对上述问题，李国发等人于2011年提出了基于目标函数的地表一致性反褶积方法，并对理论模型进行了实验处理，取得了较为理想的处理效果。但由于该方法没有解决实际地震记录的目标函数估算和非地质因素影响的最小相位化两个关键问题，未能对实际地震记录进行处理。

发明内容

为了克服现有技术存在的上述问题，本发明提供了一种基于地表一致性反褶积进行储层预测的方法及设备，解决了利用实际地震记录估计目标函数和非地质因素影响最小相位化两个关键问题，消除了采集和近地表因素对地震记录的影响，提高了地震资料反映真实地下结构的能力，进一步提高了储层预测的准确性。

本发明的目的之一是，提供一种基于地表一致性反褶积进行储层预测的方法，包括：采集来自地下的反射地震信号形成的地震记录；根据所述的地震记录生成目标函数；对所述的地震记录进行地表一致性分解，得到炮点分量和检波点分量；对所述的炮点分量和检波点分量进行最小相位化处理，得到最小相位炮检响应；根据所述的目标函数以及最小相位炮检响应得到消除非地质因素影响后的地震记录；根据所述的消除非地质因素影响后的地震记录绘制地震剖面图像；对所述的地震剖面图像进行综合解释，得到储层预测的结果。

本发明的目的之一是，提供了一种基于地表一致性反褶积进行储层预测的设备，包括：采集装置，用于采集来自地下的反射地震信号形成的地震记录；目标函数生成装置，用于根据所述的地震记录生成目标函数；分解装置，用于对所述的地震记录进行地表一致性分解，得到炮点分量和检波点分量；最小相位化处理装置，用于对所述的炮点分量和检波点分量进行最小相位化处理，得到最小相位炮检响应；影响消除装置，用于根据所述的目标函数以及最小相位炮检响应得到消除非地质因素影响后的地震记录；剖面图像生成装置，用于根据所述的消除非地质因素影响后的地震记录绘制地震剖面图像；储层预测装置，用于对所述的地震剖面图像进行综合解释，得到储层预测的结果。

本发明的有益效果在于，解决了利用实际地震记录估计目标函数和非地质因素影响最小相位化两个关键问题，消除了采集和近地表因素对地震记录的影响，恢复地下结构的真实反射特征，在子波一致性的基础上，提高地震数据处理质量，增强地震信号反映地下结构的能力，提高地震勘探精度。

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

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例提供的一种基于地表一致性反褶积进行储层预测的方法的实施方式一的流程图；

图2为图1中的步骤S102的具体流程图；

图3为图1中的步骤S103的具体流程图；

图4为图1中的步骤S104的具体流程图；

图5为本发明实施例提供的一种基于地表一致性反褶积进行储层预测的方法的实施方式二的流程图；

图6为本发明实施例提供的一种基于地表一致性反褶积进行储层预测的设备的实施方式一的结构框图；

图7为本发明实施例提供的一种基于地表一致性反褶积进行储层预测的设备中目标函数生成装置的结构框图；

图8为本发明实施例提供的一种基于地表一致性反褶积进行储层预测的设备中的分解装置的结构框图；

图9为本发明实施例提供的一种基于地表一致性反褶积进行储层预测的设备中的最小相位化处理装置的机构框图；

图10为本发明实施例提供的一种基于地表一致性反褶积进行储层预测的设备的实施方式二的结构框图；

图11是处理之前的叠前时间偏移地震记录图像；

图12是利用本发明的方法及设备处理之后的叠前时间偏移地震记录图像；

图13是某工区采集的共炮点道集地震记录图；

图14是本发明处理之前的某工区的叠前时间偏移地震记录图像；

图15是利用本发明提供的方法及设备处理之后的某工区的叠前时间偏移地震记录图像。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

图1为本发明实施例提供的一种基于地表一致性反褶积进行储层预测的方法的实施方式一的流程图，由图1可知，所述方法包括：

S101：采集来自地下的反射地震信号形成的地震记录。在具体的实施方式中，首先可由人工激发地震波，利用检波器探测来自地下的反射地震信号，形成地震记录x_ij(t),i＝1,2,…n_s;j=1,2,…n_c，其中，t是时间，单位是毫秒，n_s为炮点的个数，n_c是记录道数。图11是叠前时间偏移的地震记录图像，图中箭头指向的部位出现了由于采集和地表因素引起的振幅异常和构造异常。

S102：根据所述的地震记录生成目标函数；

S103：对所述的地震记录进行地表一致性分解，得到炮点分量和检波点分量；

S104：对所述的炮点分量和检波点分量进行最小相位化处理，得到最小相位炮检响应；

S105：根据所述的目标函数以及最小相位炮检响应得到消除非地质因素影响后的地震记录；

S106：根据所述的消除非地质因素影响后的地震记录绘制地震剖面图像；

S107：对所述的地震剖面图像进行综合解释，得到储层预测的结果。

图2为图1中的步骤S 102的具体流程图，由图2可知，步骤S102具体包括：

S201：确定所述的地震记录对应的平均振幅谱，平均振幅谱为f为频率，单位是赫兹；

S202：对所述的平均振幅谱进行平滑处理，得到平滑后的平均振幅谱，平滑后的平均振幅谱为通过如下公式得到：

$\stackrel{\‾}{B}\left(f\right)=\frac{1}{T}\underset{\mathrm{\τ}=-T}{\overset{T}{\mathrm{\Σ}}}\cos \left(\frac{\mathrm{\π\τ}}{2T}\right)\stackrel{\‾}{A}(f-\mathrm{\τ})$

其中，T为平滑算子长度，单位是赫兹，3.0≤T≤7.0。

S203：判断所述的地震记录中是否存在面波干扰，当判断为是时，执行步骤S204，否则执行步骤S205。面波干扰具有低频、低速、频散、强能量、扫帚状等特征，因此很容易判断地震记录中是否存在面波干扰。

S204：将所述的平滑后的平均振幅谱的第二峰值所对应的频率作为主频；

S205：将所述的平滑后的平均振幅谱的第一峰值所对应的频率作为主频；

S206：根据所述的主频确定目标函数。在具体的实施方式中，主频为f_p，目标函数为v(t)，目标函数通过如下公式得到：

$v\left(t\right)=[1-2{\left(\mathrm{\π}{f}_{p}t\right)}^2]e^{-{\left(\mathrm{\π}{f}_{p}t\right)}^2}$

图3为图1中的步骤S103的具体流程图，由图3可知，所述的步骤S103具体包括：

S301：对所述的地震记录进行傅里叶变换，得到对应的振幅谱，在具体的实施方式中，地震记录的振幅谱为A_ij(f)，通过如下公式得到：

A_ij(f)=|∫x_ij(t)e^-i2πftdt|

S302：对所述的振幅谱取对数，得到与所述的地震记录对应的对数谱，在具体的实施方式中，对数谱为B_ij(f)；

S303：将所述的对数谱分解为炮点项、检波点项、炮检距项以及中心点项。在具体的实施方式中，将地震记录的对数谱B_ij(f)分解为炮点项S_i(f)、检波点项G_j(f)、炮检距项H_k(f)和中心点项C_l(f)，即

B_ij(f)=S_i(f)+G_j(f)+H_k(f)+E_l(f)

其中，i,j,k,l分别表示炮点、检波点、炮检距和中心点的序号。

S304：将所有的地震记录进行分解，得到矩阵方程。即把所有的地震记录都进行步骤S303的分解，构成一个线性代数方程组，其矩阵形成表示为DP=B，其中，P为由炮点项S_i(f)、检波点项G_j(f)、炮检距项H_k(f)和中心点项C_l(f)构成的向量；B为对数谱B_ij(f)构成的向量；D为系数矩阵，它是一个每一行只有4个元素为1，其它元素均为0的大型稀疏矩阵，其维数为(n_s×n_c)×(n_s+n_g+n_h+n_ε)，其中，n_s为炮点个数，n_c为记录道数，n_g为检波点个数，n_h为炮检距个数，n_ε为中心点个数。

S305：求解所述的矩阵方程得到每个地震道对应的炮点分量的对数谱以及检波点分量对应的对数谱。在具体的实施方式中，可利用高斯－赛得尔（Gauss-Seidel）方法求解矩阵方程，得到每个地震道所对应的炮点分量的对数谱S_i(f)以及检波点分量的对数谱G_j(f)；

S306：对所述的炮点分量对应的对数谱以及检波点分量对应的对数谱分别取指数；

S307：对取指数后的所述的炮点分量对应的对数谱以及检波点分量对应的对数谱进行傅里叶反变换，得到时间域的所述的地震记录对应的炮点分量以及检波点分量。在具体的实施方式中，时间域的炮点分量可为s_i(t)，检波点分量可为g_j(t)。

图4为图1中的步骤S104的具体流程图，由图4可知，所述的步骤S104具体包括：

S401：将所述的炮点分量与所述的检波点分量进行褶积，并计算褶积结果的自相关函数，在具体的实施方式中，炮点分量为s_i(t)，检波点分量为g_j(t)，褶积结果的自相关函数为r_ij(t)；

S402：根据所述的自相关函数确定脉冲反褶积算子，在具体的实施方式中，脉冲反褶积算子为a_ij(t)；

S403：确定脉冲反褶积算子对应的自相关函数，在具体的实施方式中，脉冲反褶积算子a_ij(t)对应的自相关函数为

S404：根据所述的脉冲反褶积算子对应的自相关函数确定所述的地震记录对应的最小相位炮检响应。在具体的实施方式中，地震记录x_ij(t)所对应的最小相位炮检响应为w_ij(t)。

图5为本发明实施例提供的一种基于地表一致性反褶积进行储层预测的方法的实施方式二的流程图，由图5可知，图1中的步骤S105具体包括：

S505：将所述的目标函数作为希望输出；

S506：对所述的最小相位炮检响应进行整型滤波处理，得到反褶积滤波器，在具体的实施方式中，反褶积滤波器为d_ij(t)；

S507：将所述的反褶积滤波器与所述的地震记录进行褶积，得到消除非地质因素影响后的地震记录。在具体的实施方式中，消除非地质因素影响后的地震记录为y_ij(t)。图12是处理之后的叠前时间偏移的地震记录图像，由图12可知，本发明很好地消除了采集因素和地表因素引起的振幅异常和构造异常，恢复了地下结构的真实反射特征。

图6为本发明实施例提供的一种基于地表一致性反褶积进行储层预测的设备的实施方式一的结构框图，由图1可知，所述的设备包括：

采集装置100，用于采集来自地下的反射地震信号形成的地震记录。在具体的实施方式中，首先可由人工激发地震波，利用检波器探测来自地下的反射地震信号，形成地震记录x_ij(t),i＝1,2,…n_s;j=1,2,…n_c，其中，t是时间，单位是毫秒，n_s为炮点的个数，n_c是记录道数。图11是叠前时间偏移的地震记录图像，图中箭头指向的部位出现了由于采集和地表因素引起的振幅异常和构造异常。

目标函数生成装置200，用于根据所述的地震记录生成目标函数；

分解装置300，用于对所述的地震记录进行地表一致性分解，得到炮点分量和检波点分量；

最小相位化处理装置400，用于对所述的炮点分量和检波点分量进行最小相位化处理，得到最小相位炮检响应；

影响消除装置500，用于根据所述的目标函数以及最小相位炮检响应得到消除非地质因素影响后的地震记录；

剖面图像生成装置600，用于根据所述的消除非地质因素影响后的地震记录绘制地震剖面图像；

储层预测装置700，用于对所述的地震剖面图像进行综合解释，得到储层预测的结果。

图7为本发明实施例提供的一种基于地表一致性反褶积进行储层预测的设备中目标函数生成装置的结构框图，由图7可知，所述的目标函数生成装置200具体包括：

平均振幅谱确定单元201，用于确定所述的地震记录对应的平均振幅谱，平均振幅谱为f为频率，单位是赫兹；

平滑处理单元202，用于对所述的平均振幅谱进行平滑处理，得到平滑后的平均振幅谱平滑后的平均振幅谱为通过如下公式得到：

$\stackrel{\‾}{B}\left(f\right)=\frac{1}{T}\underset{\mathrm{\τ}=-T}{\overset{T}{\mathrm{\Σ}}}\cos \left(\frac{\mathrm{\π\τ}}{2T}\right)\stackrel{\‾}{A}(f-\mathrm{\τ})$

其中，T为平滑算子长度，单位是赫兹，3.0≤T≤7.0。

判断单元203，用于判断所述的地震记录中是否存在较强的面波干扰，面波干扰具有低频、低速、频散、强能量、扫帚状等特征，据此特征很容易判断地震记录中是否存在面波干扰。；

第一主频确定单元204，用于当判断单元的判断结果为是时，将所述的平滑后的平均振幅谱的第二峰值所对应的频率作为主频；

第二主频确定单元205，用于当判断单元的判断结果为否时，将所述的平滑后的平均振幅谱的第一峰值所对应的频率作为主频；

目标函数确定单元206，用于根据所述的主频确定目标函数。在具体的实施方式中，主频为f_p，目标函数为v(t)，目标函数通过如下公式得到：

$v\left(t\right)=[1-2{\left(\mathrm{\π}{f}_{p}t\right)}^2]e^{-{\left(\mathrm{\π}{f}_{p}t\right)}^2}$

图8为本发明实施例提供的一种基于地表一致性反褶积进行储层预测的设备中的分解装置的结构框图，由图8可知，所述的分解装置300具体包括：

振幅谱确定单元301，用于对所述的地震记录进行傅里叶变换，得到对应的振幅谱，在具体的实施方式中，地震记录的振幅谱为A_ij(f)，通过如下公式得到：

A_ij(f)=|∫x_ij(t)e^-i2πftdt|

对数谱确定单元302，用于对所述的振幅谱取对数，得到与所述的地震记录对应的对数谱，在具体的实施方式中，对数谱为B_ij(f)；

分解单元303，用于将所述的对数谱分解为炮点项、检波点项、炮检距项以及中心点项，在具体的实施方式中，将地震记录的对数谱B_ij(f)分解为炮点项S_i(f)、检波点项G_j(f)、炮检距项H_k(f)和中心点项C_l(f)，即

B_ij(f)=S_i(f)+G_j(f)+H_k(f)+E_l(f)

其中，i,j,k,l分别表示炮点、检波点、炮检距和中心点的序号。

矩阵方程生成单元304，用于将所有的地震道进行分解，得到矩阵方程。即把所有地震记录都进行步骤S303的分解，构成一个线性代数方程组，其矩阵形成表示为DP=B，其中，P为由炮点项S_i(f)、检波点项G_j(f)、炮检距项H_k(f)和中心点项C_l(f)构成的向量；B为对数谱B_ij(f)构成的向量；D为系数矩阵，它是一个每一行只有4个元素为1，其它元素均为0的大型稀疏矩阵，其维数为(n_s×n_c)×(n_s+n_g+n_h+n_ε)，其中，n_s为炮点个数，n_c为记录道数，ng为检波点个数，n_h为炮检距个数，n_εn_e为中心点个数。

矩阵方程求解单元305，用于求解所述的矩阵方程得到每个地震道对应的炮点分量的对数谱以及检波点分量对应的对数谱；在具体的实施方式中，可利用高斯－赛得尔（Gauss-Seidel）方法求解矩阵方程，得到每个地震道所对应的炮点分量的对数谱S_i(f)以及检波点分量的对数谱G_j(f)；

取指数单元306，用于对所述的炮点分量对应的对数谱以及检波点分量对应的对数谱分别取指数；

分量确定单元307，用于对取指数后的所述的炮点分量对应的对数谱以及检波点分量对应的对数谱进行傅里叶反变换，得到时间域的所述的地震记录对应的炮点分量以及检波点分量。在具体的实施方式中，时间域的炮点分量可为s_i(t)，检波点分量可为g_j(t)。

图9为本发明实施例提供的一种基于地表一致性反褶积进行储层预测的设备中的最小相位化处理装置的机构框图，由图9可知，所述的最小相位化处理装置400具体包括：

褶积单元401，用于将所述的炮点分量与所述的检波点分量进行褶积，并计算褶积结果自相关函数，在具体的实施方式中，炮点分量为s_i(t)，检波点分量为g_j(t)，褶积结果的自相关函数为r_ij(t)；

脉冲反褶积算子确定单元402，用于根据所述的自相关函数确定脉冲反褶积算子，在具体的实施方式中，脉冲反褶积算子为a_ij(t)；

自相关函数确定单元403，用于确定脉冲反褶积算子对应的自相关函数，在具体的实施方式中，脉冲反褶积算子a_ij(t)对应的自相关函数为

最小相位炮检响应确定单元404，用于根据所述的脉冲反褶积算子对应的自相关函数确定所述的地震记录对应的最小相位炮检响应。在具体的实施方式中，地震记录x_ij(t)所对应的最小相位炮检响应为w_ij(t)。

图10为本发明实施例提供的一种基于地表一致性反褶积进行储层预测的设备的实施方式二的结构框图，由图10可知，所述的影响消除装置500具体包括：

希望输出确定单元501，用于将所述的目标函数作为希望输出；

整形滤波处理单元502，用于对所述的最小相位炮检响应进行整型滤波处理，得到反褶积滤波器，在具体的实施方式中，反褶积滤波器为d_ij(t)；

影响消除单元503，用于将所述的反褶积滤波器与所述的地震记录进行褶积，得到消除非地质因素影响后的地震记录。在具体的实施方式中，消除非地质因素影响后的地震记录为y_ij(t)。图12是处理之后的叠前时间偏移的地震记录图像，由图12可知，本发明很好地消除了采集因素和地表因素引起的振幅异常和构造异常，恢复了地下结构的真实反射特征。

常规的地表一致性反褶积方法在进行炮点分量和检波点分量分解之后，对两个分量进行褶积得到采集和地表因素对每个地震道的整体响应，然后对每个地震道的响应进行相同步长的预测反褶积，由此消除非地质因素对地震子波的影响。当野外采集和地表因素变化较大时，由于预测反褶积没有明确的目标函数，很难彻底消除非地质因素的影响。本发明利用目标函数驱动反褶积滤波器对地震数据进行地表一致性处理，根据地震记录上地表一致性特征的结构差异，自适应地产生匹配滤波器，消除非地质因素对地震反射特征的影响，较常规方法具有更强的地表一致性处理能力。本发明的关键点包括：

（1）由地震记录估算反褶积目标函数；

（2）非地质因素影响最小相位化处理；

（3）目标函数驱动的地表一致性反褶积。

下面结合具体的实施例，详细介绍本发明的技术方案。将本发明应用在大港油田和华北油田多个勘探区块，实验结果表明本发明很好地消除了由于激发、接收和近地表因素对地震反射特征的影响，凸显了地震资料的真实反射特征，改善了地震资料处理质量，提高了地震资料反映真实地下结构的能力。以下仅以将本发明应用在其中一个勘探区块为例，进行详细说明。

本实施例为某油田A区块的应用实例，该区块属于海陆交互带地震勘探，施工区域横跨滩海和陆地两种地表条件，两种地表条件采用了不同类型的激发和接收系统，野外采集因素变化较大。共采集了2000炮地震记录，每炮有960个检波器进行接收，图13是一张典型的原始地震数据共炮点道集，左侧的地震道为海上检波器接收，右侧的地震道为陆地检波器接收，两者的差异十分明显。

在本实施例中，人工激发地震波，时间t=6000ms，炮点的个数n_s＝2000，记录道数n_c＝960，平滑算子长度T＝5Hz，主频f_p＝30Hz。

图14是本发明处理之前的叠前时间偏移地震记录剖面图像，箭头左侧是滩海区域地震记录，箭头右侧是陆地区域的地震记录，由于野外采集和地表因素的影响，箭头两侧的反射特征存在明显差异。图15是本发明处理之后叠前时间偏移地震记录的剖面图像，消除了由于采集和地表因素对地震记录的影响，箭头两侧滩海和陆地区域的地震反射特征趋于一致，地震剖面图像真实地反映了地下构造的结构特征。

综上所述，本发明的有益成果是：提供了一种基于地表一致性反褶积进行储层预测的方法及设备，解决了利用实际地震记录估计目标函数和非地质因素影响最小相位化两个关键问题，消除了采集和近地表因素对地震记录的影响，恢复地下结构的真实反射特征，在子波一致性的基础上，提高地震数据处理质量，增强地震信号反映地下结构的能力，提高地震勘探精度。

本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。

Claims (8)

1.一种基于地表一致性反褶积进行储层预测的方法，其特征是，所述的方法包括：

采集来自地下的反射地震信号形成的地震记录；

根据所述的地震记录生成目标函数；

对所述的地震记录进行地表一致性分解，得到炮点分量和检波点分量；

对所述的炮点分量和检波点分量进行最小相位化处理，得到最小相位炮检响应；

根据所述的目标函数以及最小相位炮检响应得到消除非地质因素影响后的地震记录；

根据所述的消除非地质因素影响后的地震记录绘制地震剖面图像；

对所述的地震剖面图像进行综合解释，得到储层预测的结果；

其中，对所述的炮点分量和检波点分量进行最小相位化处理，得到最小相位炮检响应包括：将所述的炮点分量与所述的检波点分量进行褶积，并计算褶积结果的自相关函数；根据所述的自相关函数确定脉冲反褶积算子；确定脉冲反褶积算子对应的自相关函数；根据所述的脉冲反褶积算子对应的自相关函数确定所述的地震记录对应的最小相位炮检响应。

2.根据权利要求1所述的方法，其特征是，根据所述的地震记录生成目标函数具体包括：

确定所述的地震记录对应的平均振幅谱；

对所述的平均振幅谱进行平滑处理，得到平滑后的平均振幅谱；

判断所述的地震记录中是否存在面波干扰；

当判断为是时，将所述的平滑后的平均振幅谱的第二峰值所对应的频率作为主频；

否则，将所述的平滑后的平均振幅谱的第一峰值所对应的频率作为主频；

根据所述的主频确定目标函数；

其中，对所述的平均振幅谱进行平滑处理通过如下公式进行：

$\stackrel{\‾}{B}\left(f\right)=\frac{1}{T}\underset{\mathrm{\τ}=-T}{\overset{T}{\mathrm{\Σ}}}\cos \left(\frac{\mathrm{\π\τ}}{2T}\right)\stackrel{\‾}{A}(f-\mathrm{\τ})$

T为平滑算子长度，其单位为赫兹，为平均振幅谱，f为频率，τ为平滑算子长度变量；

所述的目标函数为：

$v\left(t\right)=[1-2{\left(\mathrm{\π}{f}_{p}t\right)}^2]e^{-{\left(\mathrm{\π}{f}_{p}t\right)}^2}$

v(t)为目标函数，f_p为主频，t为时间变量。

3.根据权利要求1所述的方法，其特征是，对所述的地震记录进行地表一致性分解，得到炮点分量和检波点分量具体包括：

对所述的地震记录进行傅里叶变换，得到对应的振幅谱；

对所述的振幅谱取对数，得到与所述的地震记录对应的对数谱；

将所述的对数谱分解为炮点项、检波点项、炮检距项以及中心点项；

将所有的地震道进行分解，得到矩阵方程；

求解所述的矩阵方程得到每个地震道对应的炮点分量的对数谱以及检波点分量对应的对数谱；

对所述的炮点分量对应的对数谱以及检波点分量对应的对数谱分别取指数；

对取指数后的所述的炮点分量对应的对数谱以及检波点分量对应的对数谱进行傅里叶反变换，得到时间域的所述的地震记录对应的炮点分量以及检波点分量。

4.根据权利要求1所述的方法，其特征是，根据所述的目标函数以及最小相位炮检响应，得到消除非地质因素影响后的地震记录具体包括：

将所述的目标函数作为希望输出；

对所述的最小相位炮检响应进行整型滤波处理，得到反褶积滤波器；

将所述的反褶积滤波器与所述的地震记录进行褶积，得到消除非地质因素影响后的地震记录。

5.一种基于地表一致性反褶积进行储层预测的设备，其特征是，所述的设备包括：

采集装置，用于采集来自地下的反射地震信号形成的地震记录；

目标函数生成装置，用于根据所述的地震记录生成目标函数；

分解装置，用于对所述的地震记录进行地表一致性分解，得到炮点分量和检波点分量；

最小相位化处理装置，用于对所述的炮点分量和检波点分量进行最小相位化处理，得到最小相位炮检响应，包括褶积单元，用于将所述的炮点分量与所述的检波点分量进行褶积，并计算褶积结果的自相关函数；脉冲反褶积算子确定单元，用于根据所述的自相关函数确定脉冲反褶积算子；自相关函数确定单元，用于确定脉冲反褶积算子对应的自相关函数；最小相位炮检响应确定单元，用于根据所述的脉冲反褶积算子对应的自相关函数确定所述的地震记录对应的最小相位炮检响应；

影响消除装置，用于根据所述的目标函数以及最小相位炮检响应得到消除非地质因素影响后的地震记录；

剖面图像生成装置，用于根据所述的消除非地质因素影响后的地震记录绘制地震剖面图像；

储层预测装置，用于对所述的地震剖面图像进行综合解释，得到储层预测的结果。

6.根据权利要求5所述的设备，其特征是，所述的目标函数生成装置具体包括：

平均振幅谱确定单元，用于确定所述的地震记录对应的平均振幅谱；

平滑处理单元，用于对所述的平均振幅谱进行平滑处理，得到平滑后的平均振幅谱；

判断单元，用于判断所述的地震记录中是否存在面波干扰；

第一主频确定单元，用于当判断单元的判断结果为是时，将所述的平滑后的平均振幅谱的第二峰值所对应的频率作为主频；

第二主频确定单元，用于当判断单元的判断结果为否时，将所述的平滑后的平均振幅谱的第一峰值所对应的频率作为主频；

目标函数确定单元，用于根据所述的主频确定目标函数；

其中，所述的平滑处理单元通过如下公式进行：

$\stackrel{\‾}{B}\left(f\right)=\frac{1}{T}\underset{\mathrm{\τ}=-T}{\overset{T}{\mathrm{\Σ}}}\cos \left(\frac{\mathrm{\π\τ}}{2T}\right)\stackrel{\‾}{A}(f-\mathrm{\τ})$

T为平滑算子长度，其单位为赫兹，为平均振幅谱，f为频率，τ为平滑算子长度变量；

所述的目标函数为：

$v\left(t\right)=[1-2{\left(\mathrm{\π}{f}_{p}t\right)}^2]e^{-{\left(\mathrm{\π}{f}_{p}t\right)}^2}$

v(t)为目标函数，f_p为主频，t为时间变量。

7.根据权利要求5所述的设备，其特征是，所述的分解装置具体包括：

振幅谱确定单元，用于对所述的地震记录进行傅里叶变换，得到对应的振幅谱；

对数谱确定单元，用于对所述的振幅谱取对数，得到与所述的地震记录对应的对数谱；

分解单元，用于将所述的对数谱分解为炮点项、检波点项、炮检距项以及中心点项；

矩阵方程生成单元，用于将所有的地震道进行分解，得到矩阵方程；

矩阵方程求解单元，用于求解所述的矩阵方程得到每个地震道对应的炮点分量的对数谱以及检波点分量对应的对数谱；

取指数单元，用于对所述的炮点分量对应的对数谱以及检波点分量对应的对数谱分别取指数；

分量确定单元，用于对取指数后的所述的炮点分量对应的对数谱以及检波点分量对应的对数谱进行傅里叶反变换，得到时间域的所述的地震记录对应的炮点分量以及检波点分量。

8.根据权利要求5所述的设备，其特征是，所述的影响消除装置具体包括：

希望输出确定单元，用于将所述的目标函数作为希望输出；

整形滤波处理单元，用于对所述的最小相位炮检响应进行整型滤波处理，得到反褶积滤波器；

影响消除单元，用于将所述的反褶积滤波器与所述的地震记录进行褶积，得到消除非地质因素影响后的地震记录。