一种去除炮点虚反射干扰的方法及装置与流程

本申请涉及地质勘探技术领域，特别涉及一种去除炮点虚反射干扰的方法及装置。

背景技术

随着地震勘探技术的发展，海上油气勘探的难度和深度也越来越大，对地震资料的信噪比和分辨率要求也越来越高。海上拖缆采集技术，把检波器固定在拖缆上，可以获得高分辨率的三维地震数据。具体的，海上拖缆采集技术可以是把固定排列的检波器固定在拖缆上接收地震波，同一艘航船拖拽气枪震源进行地震波激发。可以通过全球定位系统(globalpositioningsystem，简称为gps)或者基于海岸的无线电定位来提供炮船的精确位置，同时也要对拖缆上检波器的位置进行精确定位以及记录。

在进行数据采集时，气枪震源发出的一个反射地震子波从震源位置向下传播到达海底，通过海底这个强反射界面来反射这个地震子波。反射后的地震子波继续向上前进，到达拖缆位置的检波器处，检波器感应并记录下这个反射地震子波，这是第一次反射波，是地震数据中的有效反射信号。同时，气枪震源发出的这个反射地震子波还会从震源位置向上前进到达海面，由于海面是一个较强的反射界面，该反射地震子波受到海面的反射，改变方向向下传播到达海底后，由于海底是一个强反射界面，它会发射这个地震子波，反射子波继续向上前进再次到达拖缆中的检波器，此时检波器再次感应并记录下这个反射地震子波，而这个通过海平面反射后再经过海底反射并由检波器记录下的反射地震子波是无效反射信号，即海平面炮点虚反射多次波。海平面炮点虚反射多次波不但干涉了地震数据中的有效反射信号，而且会造成地震反射子波产生非常深的频率凹陷。

目前，可以通过常规的反褶积等技术来消除海平面虚反射多次波产生的噪声干扰。然而，采用反褶积很难完全消除海水地震数据中海平面虚反射多次波的噪声干扰，也难以完全恢复虚反射多次波噪声干扰产生的频率凹陷。因此，亟需一种有效消除地震数据中海平面产生的炮点虚反射多次波干扰影响的方法。

技术实现要素：

本申请实施例的目的是提供一种去除炮点虚反射干扰的方法及装置，以有效消除地震数据中海平面产生的炮点虚反射多次波干扰影响，从而达到提高海上工区内地震数据信噪比和分辨率的目的。

为解决上述技术问题，本申请实施例提供一种去除炮点虚反射干扰的方法及装置是这样实现的：

一种去除炮点虚反射干扰的方法，提供有海上工区内目标检波点在时间域中的地震数据、所述目标检波点到所述海上工区内的炮点的距离、所述目标检波点和所述炮点分别到海平面的距离，以及所述海平面至海底反射界面的距离；所述方法包括：

根据所述目标检波点到所述海上工区内的炮点的距离、所述目标检波点和所述炮点分别到海平面的距离，以及所述海平面至海底反射界面的距离，确定所述目标检波点对应的炮点虚反射延迟的延迟时间；

基于所述延迟时间，确定所述目标检波点对应的去除炮点虚反射的算子；

根据所述地震数据和所述算子，确定所述目标检波点对应的去除炮点虚反射干扰后地震数据的时域值。

优选方案中，所述方法还提供有海水速度；确定所述目标检波点对应的炮点虚反射延迟的延迟时间，包括：

根据所述目标检波点到所述海上工区内的炮点的距离、所述目标检波点和所述炮点分别到海平面的距离、所述海平面至海底反射界面的距离，以及所述海水速度，分别确定所述炮点到所述目标检波点的一次反射时间和虚反射时间；

将所述虚反射时间减去所述一次反射时间，得到所述延迟时间。

优选方案中，采用下述公式分别确定所述一次反射时间和所述虚反射时间：

其中，和分别表示所述炮点到所述海上工区内的第n个检波点的一次反射时间和虚反射时间，所述第n个检波点表示所述目标检波点；xn表示所述炮点到所述目标检波点的距离，d表示所述海平面至海底反射界面的距离，zr和zs分别表示所述目标检波点和所述炮点分别到海平面的距离，v表示所述海水速度。

优选方案中，所述方法还提供有所述海上工区对应的多个海平面反射系数；确定所述目标检波点对应的去除炮点虚反射的算子，包括：

根据所述多个海平面反射系数，分别确定目标反射系数和目标平方和反射系数；

基于所述延迟时间和所述目标反射系数，确定所述目标检波点对应的平均复共轭虚反射因子；基于所述延迟时间、所述目标反射系数和所述目标平方和反射系数，确定所述目标检波点对应的平均平方和虚反射因子；

根据所述平均复共轭虚反射因子和所述平均平方和虚反射因子，确定所述目标检波点对应的去除炮点虚反射的算子。

优选方案中，采用下述公式确定所述目标反射系数：

其中，r表示所述目标反射系数，ri表示所述多个海平面反射系数中第i个海平面反射系数，nk表示海平面反射系数的数量；

采用下述公式确定所述目标平方和反射系数：

其中，rr表示所述目标平方和反射系数，ri表示所述多个海平面反射系数中第i个海平面反射系数，nk表示海平面反射系数的数量。

优选方案中，所述方法还提供有对所述目标检波点在时间域中的地震数据进行傅里叶变换后得到的所述目标检波点在频率域中的地震数据；采用下述公式确定所述平均复共轭虚反射因子：

其中，表示所述平均复共轭虚反射因子，k表示所述目标检波点在频率域中的地震数据中的样点顺序号，k＝1,2,3,…,nl，nl表示所述目标检波点在频率域中的地震数据中的样点总个数，r表示所述目标反射系数，δf表示频率采样间隔，δt表示所述目标检波点在时间域中的地震数据的采样间隔，j是虚数单位，j²＝-1，表示所述海上工区内的第n个检波点对应的炮点虚反射延迟的延迟时间；其中，所述第n个检波点表示所述目标检波点。

优选方案中，所述方法还提供有对所述目标检波点在时间域中的地震数据进行傅里叶变换后得到的所述目标检波点在频率域中的地震数据；采用下述公式确定所述平均平方和虚反射因子：

其中，表示所述平均平方和虚反射因子，k表示所述目标检波点在频率域中的地震数据中的样点顺序号，k＝1,2,3,…,nl，nl表示所述目标检波点在频率域中的地震数据中的样点总个数，r表示所述目标反射系数，rr表示所述目标平方和反射系数，δf表示频率采样间隔，δt表示所述目标检波点在时间域中的地震数据的采样间隔，表示所述海上工区内的第n个检波点对应的炮点虚反射延迟的延迟时间；其中，所述第n个检波点表示所述目标检波点。

优选方案中，采用下述公式确定所述目标检波点对应的去除炮点虚反射的算子：

wnl＝e^-j2π/nl

其中，gn[m]表示所述海上工区内的第n个检波点对应的去除炮点虚反射的算子，所述第n个检波点表示所述目标检波点；m表示在时间域中的去除炮点虚反射的算子中的样点顺序号，m＝1,2,3,…,ng，ng表示在时间域中的去除炮点虚反射的算子中的样点总个数，g[k]表示所述目标检波点对应的去除炮点虚反射因子，k表示所述目标检波点在频率域中的地震数据中的样点顺序号，k＝1,2,3,…,nl，nl表示所述目标检波点在频率域中的地震数据中的样点总个数，表示所述平均复共轭虚反射因子，表示所述平均平方和虚反射因子，j是虚数单位，j²＝-1。

优选方案中，采用下述公式确定所述目标检波点对应的去除炮点虚反射干扰后地震数据的时域值：

其中，yn[l]表示所述海上工区内的第n个检波点对应的去除炮点虚反射干扰后地震数据的时域值，所述第n个检波点表示所述目标检波点，l表示所述目标检波点对应的去除炮点虚反射干扰后地震数据中的样点顺序号，l＝1,2,3,…,nn，nn表示所述目标检波点在时间域中的地震数据中的样点总个数，gn[m]表示所述目标检波点对应的去除炮点虚反射的算子，m＝1,2,3,…,ng，ng表示在时间域中的去除炮点虚反射的算子中的样点总个数，sn[l-m]表示所述目标检波点在时间域中的地震数据。

一种去除炮点虚反射干扰的装置，所述装置提供海上工区内目标检波点在时间域中的地震数据、所述目标检波点到所述海上工区内的炮点的距离、所述目标检波点和所述炮点分别到海平面的距离，以及所述海平面至海底反射界面的距离；所述装置包括：延迟时间确定模块、算子确定模块和目标地震数据确定模块；其中，

所述延迟时间确定模块，用于根据所述目标检波点到所述海上工区内的炮点的距离、所述目标检波点和所述炮点分别到海平面的距离，以及所述海平面至海底反射界面的距离，确定所述目标检波点对应的炮点虚反射延迟的延迟时间；

所述算子确定模块，用于基于所述延迟时间，确定所述目标检波点对应的去除炮点虚反射的算子；

所述目标地震数据确定模块，用于根据所述地震数据和所述算子，确定所述目标检波点对应的去除炮点虚反射干扰后地震数据的时域值。

由以上本申请实施例提供的技术方案可见，本申请实施例提供的去除炮点虚反射干扰的方法及装置，可以根据所述目标检波点到所述海上工区内的炮点的距离、所述目标检波点和所述炮点分别到海平面的距离，以及所述海平面至海底反射界面的距离，确定所述目标检波点对应的炮点虚反射延迟的延迟时间；然后，可以基于所述延迟时间，确定所述目标检波点对应的去除炮点虚反射的算子；最后，可以根据所述地震数据和所述算子，确定所述目标检波点对应的去除炮点虚反射干扰后地震数据的时域值。如此，本申请方法可以实现海上拖缆数据中去除炮点虚反射处理，以消除地震数据中海平面产生的炮点虚反射多次波干扰影响，从而达到提高海上工区内地震数据信噪比和分辨率的目的。

附图说明

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

图1是本申请一种去除炮点虚反射干扰的方法实施例的流程图；

图2是本申请实施例中炮点虚反射多次波的传播路径示意图；

图3(a)是本申请实施例中去除炮点虚反射处理前的共炮点道集数据的示意图；

图3(b)是本申请实施例中采用现有技术的方法进行去除炮点虚反射处理后的共炮点道集数据的示意图；

图3(c)是本申请实施例中采用本申请的方法进行去除炮点虚反射处理后的共炮点道集数据的示意图；

图4(a)是图3(a)中的共炮点道集数据的振幅谱示意图；

图4(b)是图3(b)中的共炮点道集数据的振幅谱示意图；

图4(c)是图3(c)中的共炮点道集数据的振幅谱示意图；

图5是本申请去除炮点虚反射干扰的装置的一种实施例的组成结构示意图；

图6是本申请去除炮点虚反射干扰的装置的另一种实施例的组成结构示意图。

具体实施方式

本申请实施例提供一种去除炮点虚反射干扰的方法及装置。

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

本申请实施例提供了一种去除炮点虚反射干扰的方法。所述去除炮点虚反射干扰的方法可以提供有海上工区内目标检波点在时间域中的地震数据、所述目标检波点到所述海上工区内的炮点的距离、所述目标检波点和所述炮点分别到海平面的距离，以及所述海平面至海底反射界面的距离。不仅如此，所述去除炮点虚反射干扰的方法还可以提供有所述海上工区对应的多个海平面反射系数

在本实施方式中，可以获取所述海上工区内目标检波点在时间域中的地震数据。还可以采用傅里叶变换的方法将在时间域中的地震数据变换至频率域，得到在频率域中的地震数据。其中，在本实施方式中，所述目标检波点可以是指所述海上工区内的任意一个检波点。例如，所述海上工区内的第n个检波点。所述检波点也可以称为接收点。即，上述获取海上工区内目标检波点在频率域中的地震数据可以等同于获取海上工区内目标接收点在频率域中的地震数据。

在本实施方式中，可以利用海上拖缆采集技术，将检波器固定在拖缆上来采集地震数据。在实际应用中也可以采用其他方法采集地震数据，本申请对此不作限定。

在本实施方式中，可以在采集地震数据的过程中，直接测量所述目标检波点到所述海上工区内的炮点的距离、所述目标检波点和所述炮点分别到海平面的距离，以及所述海平面至海底反射界面的距离。

在本实施方式中，可以对海上工区内目标检波点所采集到的地震数据进行预处理。其中，所述预处理可以是为地震数据设置标签、定义地震数据的观测系统、对采集到的地震数据进行速度分析等处理。在本实施方式中，在进行地震数据预处理的过程中，可以将采集到的地震道数据按照segy、cgg和grisys等任意一种文件格式进行保存。

在本实施方式中，可以通过资料收集的方式，获取所述海上工区对应的多个海平面反射系数。在实际应用中，每一位地震勘探人员均会根据各自经验测得该海上工区中不同区域的海平面反射系数，如此，可以通过收集这些数据，得到所述海上工区对应的多个海平面反射系数。

图1是本申请一种去除炮点虚反射干扰的方法实施例的流程图。如图1所示，所述去除炮点虚反射干扰的方法，包括以下步骤。

步骤s101：根据所述目标检波点到所述海上工区内的炮点的距离、所述目标检波点和所述炮点分别到海平面的距离，以及所述海平面至海底反射界面的距离，确定所述目标检波点对应的炮点虚反射延迟的延迟时间。

虚反射是多次反射中的一种，可以指由爆炸点向上传播，在遇到低速带底面或者地面之后，又向下反射传播，最后又从下面的反射界面再反射至地面的现象。因此，虚反射常和一次反射相伴随(故又可以称为伴随波)，或者和一次反射的尾部混合在一起，或者形成单独的波。

进一步的，对于海平面虚反射而言，反射地震子波从震源位置向下传播到达海底，由于海底是一个强反射界面，它会反射这个地震子波。反射子波继续向上前进到达拖缆位置的检波器处，检波器感应并记录这个反射地震子波，这是一次反射波，是地震数据中的有效反射信号。同时，反射地震子波还可以从震源位置向上前进到达海面，由于海面是一个较强的反射界面，该反射地震子波受到海面的反射，改变方向向下传播到达海底后，由于海底是一个强反射界面，它会发射这个地震子波，反射子波继续向上前进再次到达拖缆中的检波器，此时检波器再次感应并记录下这个反射地震子波，而这个通过海平面反射后再经过海底反射并由检波器记录下的反射地震子波是无效反射信号，即海平面炮点虚反射多次波。海平面虚反射多次波是海上工区内地震勘探数据中最大的噪声干扰。在本申请实施方式中所描述的去除炮点虚反射干扰可以指的是海平面炮点虚反射多次波所造成的噪声干扰。

在本实施方式中，可以将同一检波点所接收到的来自同一炮点(即震源)的一次反射波和虚反射多次波之间的时间间隔称为该炮点虚反射延迟的时间。

在本实施方式中，在获取所述目标检波点到所述海上工区内的炮点的距离、所述目标检波点和所述炮点分别到海平面的距离，以及所述海平面至海底反射界面的距离之后，可以根据所述目标检波点到所述海上工区内的炮点的距离、所述目标检波点和所述炮点分别到海平面的距离，以及所述海平面至海底反射界面的距离，确定所述目标检波点对应的炮点虚反射延迟的延迟时间。具体地，可以根据所述目标检波点到所述海上工区内的炮点的距离、所述目标检波点和所述炮点分别到海平面的距离、所述海平面至海底反射界面的距离，以及所述海水速度，分别确定所述炮点到所述目标检波点的一次反射时间和虚反射时间，并可以将所述虚反射时间减去所述一次反射时间，得到所述延迟时间。例如，根据图2所示的炮点虚反射多次波的传播路径示意图，可以采用下述公式确定所述海上工区内的第n个检波点对应的炮点虚反射延迟的延迟时间

其中，和分别表示所述炮点到所述海上工区内的第n个检波点的一次反射时间和虚反射时间，所述第n个检波点表示所述目标检波点。

对于如图2所示的炮点虚反射多次波的传播路径示意图而言，d表示海平面到地下反射界面的距离，即所述海平面至海底反射界面的距离；zs和zr分别表示海平面到炮点水平面和检波点水平面的距离，即所述炮点和所述目标检波点分别到海平面的距离；xn表示第n个检波点或接收点到炮点的距离，v表示海水速度。

可以采用下述公式分别确定所述炮点到所述海上工区内的第n个检波点的一次反射时间和虚反射时间：

其中，和分别表示所述炮点到所述海上工区内的第n个检波点的一次反射时间和虚反射时间。

步骤s102：基于所述延迟时间，确定所述目标检波点对应的去除炮点虚反射的算子。

在本实施方式中，基于所述延迟时间，确定所述目标检波点对应的去除炮点虚反射的算子，具体可以包括，可以根据所述多个海平面反射系数，分别确定目标反射系数和目标平方和反射系数；可以基于所述延迟时间和所述目标反射系数，确定所述目标检波点对应的平均复共轭虚反射因子；可以基于所述延迟时间、所述目标反射系数和所述目标平方和反射系数，确定所述目标检波点对应的平均平方和虚反射因子；可以根据所述平均复共轭虚反射因子和所述平均平方和虚反射因子，确定所述目标检波点对应的去除炮点虚反射的算子。

在本实施方式中，可以采用下述公式确定所述目标反射系数：

其中，r表示所述目标反射系数，ri表示所述多个海平面反射系数中第i个海平面反射系数，i＝1,2,3,…,nk，nk表示海平面反射系数的数量。

可以采用下述公式确定所述目标平方和反射系数：

其中，rr表示所述目标平方和反射系数，ri表示所述多个海平面反射系数中第i个海平面反射系数，nk表示海平面反射系数的数量。

在本实施方式中，作为一次反射延迟的伴随反射，虚反射会出现在一个地震道中。虚反射可以包括：炮点虚反射和接收点虚反射。炮点虚反射是由向下移动的一次反射能量的二次震源产生的，而这个一次反射是由向上移动的能量产生的；接收点虚反射是由向上移动的一次反射能量的二次震源产生的，而这个一次反射是由向下移动的能量产生的。使用最小相位脉冲来模拟虚反射，在时间域中可以表示为：

其中，sn(t)表示第n个检波点到震源的包含一次反射和虚反射的地震道，rn(t)表示炮点到第n个检波点的一次反射地震道，t表示地震数据的传播时间，ri表示所述多个海平面反射系数中第i个海平面反射系数，其取值范围为0.9≤ri≤1；为第n个检波点对应的炮点虚反射延迟时间，为第n个检波点对应的接收点虚反射延迟时间，为第n个检波点对应的炮点接收点虚反射延迟时间。

所述第n个检波点到震源的包含一次反射和虚反射的地震道在频率域中可以表示为：

其中，s(f)表示第n个检波点到炮点的包含一次反射和虚反射的地震道频谱，r(f)表示第n个检波点到炮点的一次反射地震道频谱，f表示频率。对于零偏移距地震数据，有

将公式(8)代入公式(7)中，有

这样，令

显然，有

将公式(10)、(11)和(12)代入公式(9)，有

s(f)＝u(f)+d(f)(13)

这样，所述第i个海平面反射系数和所述第n个检波点对应的炮点虚反射算子可以表示为：

则有

即

其中，符号h表示复共轭运算。所述第i个海平面反射系数对应的去除炮点虚反射后的地震道数据可以表示为：

对于不同海平面反射系数ri，其中，i＝1,2,3,…,nk，nk表示海平面反射系数的数量，由公式(17)，有

由公式(19)，有

这样，去除炮点虚反射因子可以表示为：

将公式(15)离散化处理后，有所述第i个海平面反射系数对应的去除炮点虚反射因子

对应的共轭去除炮点虚反射因子

其中，k表示所述目标检波点在频率域中的地震数据中的样点顺序号，k＝1,2,3,…,nl，nl表示所述目标检波点在频率域中的地震数据中的样点总个数，r表示所述目标反射系数，j是虚数单位，j²＝-1，表示所述海上工区内的第n个检波点对应的炮点虚反射延迟的延迟时间。其中，可以采用下述公式计算频率采样间隔δf：

其中，δt表示所述目标检波点在时间域中的地震数据的采样间隔。

将公式(21)进行离散化处理后，有

由公式(23)，有

如此，可以采用下述公式确定平均复共轭虚反射因子

其中，k表示所述目标检波点在频率域中的地震数据中的样点顺序号，k＝1,2,3,…,nl，nl表示所述目标检波点在频率域中的地震数据中的样点总个数，r表示所述目标反射系数，δf表示频率采样间隔，j是虚数单位，j²＝-1，表示所述海上工区内的第n个检波点对应的炮点虚反射延迟的延迟时间；其中，所述第n个检波点表示所述目标检波点。

由公式(22)，有

如此，可以采用下述公式确定所述平均平方和虚反射因子

其中，r表示所述目标反射系数，rr表示所述目标平方和反射系数。

由公式(25)，可以采用下述公式确定所述目标检波点对应的去除炮点虚反射因子g[k]：

其中，表示所述平均复共轭虚反射因子，表示所述平均平方和虚反射因子。

在本实施方式中，可以采用下述公式确定所述目标检波点对应的去除炮点虚反射的算子：

其中，wnl＝e^-j2π/nl(32)

在公式(31)和(32)中，gn[m]表示所述海上工区内的第n个检波点对应的去除炮点虚反射的算子，所述第n个检波点表示所述目标检波点；m表示在时间域中的去除炮点虚反射的算子中的样点顺序号，m＝1,2,3,…,ng，ng表示在时间域中的去除炮点虚反射的算子中的样点总个数，g[k]表示所述目标检波点对应的去除炮点虚反射因子，k表示所述目标检波点在频率域中的地震数据中的样点顺序号，k＝1,2,3,…,nl，nl表示所述目标检波点在频率域中的地震数据中的样点总个数，j是虚数单位，j²＝-1。

步骤s103：根据所述地震数据和所述算子，确定所述目标检波点对应的去除炮点虚反射干扰后地震数据的时域值。

在本实施方式中，可以采用下述公式确定所述目标检波点对应的去除炮点虚反射干扰后地震数据的时域值：

其中，yn[l]表示所述海上工区内的第n个检波点对应的去除炮点虚反射干扰后地震数据的时域值，所述第n个检波点表示所述目标检波点，l表示所述目标检波点对应的去除炮点虚反射干扰后地震数据中的样点顺序号，l＝1,2,3,…,nn，nn表示所述目标检波点在时间域中的地震数据中的样点总个数，gn[m]表示所述目标检波点对应的去除炮点虚反射的算子，m＝1,2,3,…,ng，ng表示在时间域中的去除炮点虚反射的算子中的样点总个数，sn[l-m]表示所述目标检波点在时间域中的地震数据。

在步骤s103后，本申请实施例还可以包括绘制去除炮点虚反射后的地震数据剖面，以及存储去除炮点虚反射后的地震数据。

所述去除炮点虚反射干扰的方法实施例，可以根据所述目标检波点到所述海上工区内的炮点的距离、所述目标检波点和所述炮点分别到海平面的距离，以及所述海平面至海底反射界面的距离，确定所述目标检波点对应的炮点虚反射延迟的延迟时间；然后，可以基于所述延迟时间，确定所述目标检波点对应的去除炮点虚反射的算子；最后，可以根据所述地震数据和所述算子，确定所述目标检波点对应的去除炮点虚反射干扰后地震数据的时域值。如此，本申请方法可以实现海上拖缆数据中去除炮点虚反射处理，以消除地震数据中海平面产生的炮点虚反射多次波干扰影响，从而达到提高海上工区内地震数据信噪比和分辨率的目的。不仅如此，采用本申请方法不仅可以恢复地震反射子波低频频率成分，还恢复高频频率成分，因此可以有效拓宽了有效频带，提高了拖缆地震数据分辨率。

同时，本申请方法采用直接计算炮点虚反射延迟时间，可以直接设计去除炮点虚反射干扰因子和去除炮点虚反射的算子，具有计算量小、计算速度快、稳定性好和计算精度高的特点。本申请方法引入平均平方和虚反射因子，使得虚反射算子计算更加平稳，具有很强的抗噪能力。

为了清除的说明本申请实施例的有益效果，下面结合附图进行说明：

试验数据是一实际深海拖缆采集数据，拖缆沉放深度9.5米(m)，拖缆总长6000m，工作道数480，道间距12.5m，枪沉放深度7.5m，炮点间隔是50m。实际处理时炮点深度使用7.5m，接收点深度从数据道头字中读取，海水速度取值1545米/秒(m/s)。数据时间采样率是2毫秒(ms)，最大频率250赫兹(hz)。图3(a)是本申请实施例中去除炮点虚反射处理前的共炮点道集数据的示意图；图3(b)是本申请实施例中采用现有技术的方法进行去除炮点虚反射处理后的共炮点道集数据的示意图；图3(c)是本申请实施例中采用本申请的方法进行去除炮点虚反射处理后的共炮点道集数据的示意图。图3(a)、图3(b)和图3(c)中的横坐标和纵坐标分别为道号和采样时间，采样时间的单位为ms。如图3(a)、图3(b)和图3(c)所示，采用本申请方法去除炮点虚反射处理后的地震数据中，海底反射同相轴和子波旁瓣连续性和一致性，明显好于现有常规方法，而且反射子波长度明显小于现有常规方法，因此，本申请方法可以有效地去除炮点虚反射，压缩地震子波，从而提高了地震数据分辨率。图4(a)是图3(a)中的共炮点道集数据的振幅谱示意图；图4(b)是图3(b)中的共炮点道集数据的振幅谱示意图；图4(c)是图3(c)中的共炮点道集数据的振幅谱示意图。图4(a)、图4(b)和图4(c)中的横坐标和纵坐标分别为频率和振幅，频率的单位为hz。如图4(a)、图4(b)和图4(c)所示，本申请方法去除炮点虚反射处理后数据，低频比现有常规方法更低，高频比常规方法更高，拓宽了地震数据的有效带宽。利用本申请方法可以更加有效地去除了炮点虚反射，压缩地震子波，不但恢复了低频频率成分，而且也恢复了高频频率成分，因此有效拓宽了有效频带，提高了拖缆地震数据分辨率。

图5是本申请去除炮点虚反射干扰的装置的一种实施例的组成结构示意图。所述去除炮点虚反射干扰的装置提供海上工区内目标检波点在时间域中的地震数据、所述目标检波点到所述海上工区内的炮点的距离、所述目标检波点和所述炮点分别到海平面的距离，以及所述海平面至海底反射界面的距离。如图5所示，所述去除炮点虚反射干扰的装置可以包括：延迟时间确定模块100、算子确定模块200和目标地震数据确定模块300。

所述延迟时间确定模块100，用于根据所述目标检波点到所述海上工区内的炮点的距离、所述目标检波点和所述炮点分别到海平面的距离，以及所述海平面至海底反射界面的距离，确定所述目标检波点对应的炮点虚反射延迟的延迟时间。

所述算子确定模块200，用于基于所述延迟时间，确定所述目标检波点对应的去除炮点虚反射的算子。

所述目标地震数据确定模块300，用于根据所述地震数据和所述算子，确定所述目标检波点对应的去除炮点虚反射干扰后地震数据的时域值。

图6是本申请去除炮点虚反射干扰的装置的另一种实施例的组成结构示意图。如图6所示，所述去除炮点虚反射干扰的装置可以包括存储器、处理器、以及存储在所述存储器上的计算机程序，所述存储器中存储有海上工区内目标检波点在时间域中的地震数据、所述目标检波点到所述海上工区内的炮点的距离、所述目标检波点和所述炮点分别到海平面的距离，以及所述海平面至海底反射界面的距离，所述计算机程序被所述处理器运行时执行如下步骤：

步骤s101：根据所述目标检波点到所述海上工区内的炮点的距离、所述目标检波点和所述炮点分别到海平面的距离，以及所述海平面至海底反射界面的距离，确定所述目标检波点对应的炮点虚反射延迟的延迟时间；

步骤s102：基于所述延迟时间，确定所述目标检波点对应的去除炮点虚反射的算子；

步骤s103：根据所述地震数据和所述算子，确定所述目标检波点对应的去除炮点虚反射干扰后地震数据的时域值。

所述去除炮点虚反射干扰的装置实施例与所述去除炮点虚反射干扰的方法实施例相对应，可以实现去除炮点虚反射干扰的方法实施例，并取得方法实施例的技术效果。

在20世纪90年代，对于一个技术的改进可以很明显地区分是硬件上的改进(例如，对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而，随着技术的发展，当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此，不能说一个方法流程的改进就不能用硬件实体模块来实现。例如，可编程逻辑器件(programmablelogicdevice,pld)(例如现场可编程门阵列(fieldprogrammablegatearray，fpga))就是这样一种集成电路，其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片pld上，而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且，如今，取代手工地制作集成电路芯片，这种编程也多半改用“逻辑编译器(logiccompiler)”软件来实现，它与程序开发撰写时所用的软件编译器相类似，而要编译之前的原始代码也得用特定的编程语言来撰写，此称之为硬件描述语言(hardwaredescriptionlanguage，hdl)，而hdl也并非仅有一种，而是有许多种，如abel(advancedbooleanexpressionlanguage)、ahdl(alterahardwaredescriptionlanguage)、confluence、cupl(cornelluniversityprogramminglanguage)、hdcal、jhdl(javahardwaredescriptionlanguage)、lava、lola、myhdl、palasm、rhdl(rubyhardwaredescriptionlanguage)等，目前最普遍使用的是vhdl(very-high-speedintegratedcircuithardwaredescriptionlanguage)与verilog2。本领域技术人员也应该清楚，只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中，就可以很容易得到实现该逻辑方法流程的硬件电路。

本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。

上述实施例阐明的装置、模块，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。

为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现。

通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，在一个典型的配置中，计算设备包括一个或多个处理器(cpu)、输入/输出接口、网络接口和内存。该计算机软件产品可以包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。该计算机软件产品可以存储在内存中，内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(ram)和/或非易失性内存等形式，如只读存储器(rom)或闪存(flashram)。内存是计算机可读介质的示例。计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(pram)、静态随机存取存储器(sram)、动态随机存取存储器(dram)、其他类型的随机存取存储器(ram)、只读存储器(rom)、电可擦除可编程只读存储器(eeprom)、快闪记忆体或其他内存技术、只读光盘只读存储器(cd-rom)、数字多功能光盘(dvd)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括短暂电脑可读媒体(transitorymedia)，如调制的数据信号和载波。

本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于装置实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。

本申请可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络pc、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。

本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。

虽然通过实施例描绘了本申请，本领域普通技术人员知道，本申请有许多变形和变化而不脱离本申请的精神，希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。