地震数据处理方法、装置、电子设备及计算机存储介质
阅读：435发布：2021-12-11
专利汇可以提供地震数据处理方法、装置、电子设备及计算机存储介质专利检索，专利查询，专利分析的服务。并且本 发明 的目的是提供一种 地震 数据处理 方法、装置、 电子 设备及计算机存储介质，所述方法包括：获取第一地震数据，其中，所述第一地震数据包括一次波与多次波；根据 指定 速度对所述地震数据动校正处理，获得第二地震数据；其中，所述指定速度包括所述一次 波速 度的数值、所述多次波速度的数值或一次波速度的数值与多次波速度的数值之间的数值；对所述第二地震数据进行预设处理，得到第三地震数据；所述预设处理包括，根据指定规则缩小所述第二地震数据中的所述炮检距，和/或，根据指定规则增大所述第二地震数据中的所述时间。通过本发明的实施方式使得一次波与多次波更好区分，从而获得更好的多次波滤波效果。，下面是地震数据处理方法、装置、电子设备及计算机存储介质专利的具体信息内容。
1.一种地震数据处理方法，其特征在于，所述方法包括：
获取第一地震数据，其中，所述第一地震数据包括一次波与多次波；
根据指定速度对所述地震数据动校正处理，获得第二地震数据；其中，所述指定速度包括所述一次波速度的数值、所述多次波速度的数值或一次波速度的数值与多次波速度的数值之间的数值；所述第二地震数据包括多个维度，所述多个维度包括炮检距和时间；
对所述第二地震数据进行预设处理，得到第三地震数据；所述预设处理包括，根据指定规则缩小所述第二地震数据中的所述炮检距，和/或，根据指定规则增大所述第二地震数据中的所述时间。
2.如权利要求1所述的方法，其特征在于，根据所述指定速度对所述地震数据动校正处理，获得所述第二地震数据，包括：
根据所述指定速度，获得不同炮检距处检波点对应的动校正时差数据；
将不同炮检距处检波点接收的旅行时间与所述动校正时差数据进行校正，获得所述第二地震数据；其中，所述旅行时间是地震波产生后至所述检波点所需时间。
3.如权利要求1所述的方法，其特征在于，根据所述指定规则缩小所述第二地震数据中的所述炮检距，和/或，根据所述指定规则增大所述第二地震数据中的所述时间包括：按比例缩小所述第二地震数据中的所述炮检距，和/或，按比例增大所述第二地震数据中的所述时间。
4.如权利要求1所述的方法，其特征在于，所述第三地震数据至少符合第一特征，所述第一特征为所述一次波的曲率图形特征参数与所述多次波的曲率图形特征参数不处于同一正数区间或负数区间；在获得所述第三地震数据后，所述方法还包括：根据所述第一特征，对所述第三地震数据进行压制多次波处理，获得压制多次波后的第三地震数据。
5.如权利要求4所述的方法，其特征在于，对所述第三地震数据进行压制多次波处理，获得所述压制多次波后的第三地震数据包括：
将所述第三地震数据映射至指定坐标系，获得第四地震数据；所述指定坐标系包括曲率图形特征维度和时间维度；
对所述第四地震数据进行滤波处理，获得压制多次波后的第四地震数据；
将所述压制多次波后的第四地震数据反映射处理，获得所述压制多次波后的第三地震数据。
6.如权利要求5所述的方法，其特征在于，在获得压制多次波后的第三地震数据后，还包括：
根据第二指定规则，放大所述压制多次波后的第三地震数据后中的所述炮检距，和/或，根据第二指定规则，缩小所述压制多次波后的第三地震数据后中的所述时间，获得压制多次波后的第二地震数据；其中，所述第二指定规则为对应于所述指定规则的反变换。
7.如权利要求6所述的方法，其特征在于，在获得所述压制多次波后的第二地震数据后还包括：
对所述压制多次波后的第二地震数据进行反动校正处理，获得压制多次波后的第一地震数据。
8.一种地震数据处理装置，其特征在于，所述装置包括：
获取单元，所述获取单元用于获取第一地震数据，其中，所述第一地震数据包括一次波与多次波；
第一处理单元，第一处理单元用于根据指定速度对所述地震数据动校正处理，获得第二地震数据；其中，所述指定速度包括所述一次波速度的数值、所述多次波速度的数值或一次波速度的数值与多次波速度的数值之间的数值；所述第二地震数据包括多个维度，所述多个维度包括炮检距和时间；
第二处理单元，第二处理单元用于对所述第二地震数据进行预设处理，得到第三地震数据；所述预设处理包括，根据指定规则缩小所述第二地震数据中的所述炮检距，和/或，根据指定规则增大所述第二地震数据中的所述时间。
9.一种电子设备，其特征在于，包括：输入设备，存储器，处理器；
所述输入设备用于获取第一地震数据，其中，所述第一地震数据包括一次波与多次波；
所述处理器用于根据指定速度对所述地震数据动校正处理，获得第二地震数据；其中，所述指定速度包括所述一次波速度的数值、所述多次波速度的数值或一次波速度的数值与多次波速度的数值之间的数值；所述第二地震数据包括多个维度，所述多个维度包括炮检距和时间；对所述第二地震数据进行预设处理，得到第三地震数据；所述预设处理包括，根据指定规则缩小所述第二地震数据中的所述炮检距，和/或，根据指定规则增大所述第二地震数据中的所述时间；
所述存储器用于存取数据，所述数据包括所述输入或被获取的地震数据，所述处理器处理的中间数据以及结果数据。
10.一种计算机存储介质，其特征在于，所述计算机存储介质存储有计算机程序指令，在所述计算机程序指令被执行时实现：
获取第一地震数据，其中，所述第一地震数据包括一次波与多次波；
根据指定速度对所述地震数据动校正处理，获得第二地震数据；其中，所述指定速度包括所述一次波速度的数值、所述多次波速度的数值或一次波速度的数值与多次波速度的数值之间的数值；所述第二地震数据包括多个维度，所述多个维度包括炮检距和时间；
对所述第二地震数据进行预设处理，得到第三地震数据；所述预设处理包括，根据指定规则缩小所述第二地震数据中的所述炮检距，和/或，根据指定规则增大所述第二地震数据中的所述时间。
地震数据处理方法、装置、电子设备及计算机存储介质
技术领域
背景技术
[0002] 一般情况下，在地震数据处理中多次波被视为一种严重干扰一次波的相干噪声。多次波与一次波叠加在一起，难以分离，这不仅降低了地震资料的信噪比，而且严重影响了地震成像的真实性和可靠性，对后续的解释工作带来不利影响。
[0003] 目前压制多次波的方法一般根据多次波与一次波的差异，采用滤波的方法进行压制，如Radon变换，频率-波数(F-K frequency-wavenumber)变换和预测反褶积等。由于多次波和一次波的时差可能很小，尤其是层间多次波与一次波的差异非常小，影响滤波效果。
发明内容
[0005] 为实现上述目的，本说明书实施方式提供一种地震数据处理方法，所述方法包括：获取第一地震数据，其中，所述第一地震数据包括一次波与多次波；根据指定速度对所述地震数据动校正处理，获得第二地震数据；其中，所述指定速度包括所述一次波速度的数值、所述多次波速度的数值或一次波速度的数值与多次波速度的数值之间的数值；所述第二地震数据包括多个维度，所述多个维度包括炮检距和时间；对所述第二地震数据进行预设处理，得到第三地震数据；所述预设处理包括，根据指定规则缩小所述第二地震数据中的所述炮检距，和/或，根据指定规则增大所述第二地震数据中的所述时间。在所述第二地震数据中，所述一次波与所述多次波清楚分离。
[0006] 为实现上述目的，本说明书实施方式还提供一种地震数据处理装置，所述装置包括：获取单元，所述获取单元用于获取第一地震数据，其中，所述第一地震数据包括一次波与多次波；第一处理单元，第一处理单元用于根据指定速度对所述地震数据动校正处理，获得第二地震数据；其中，所述指定速度包括所述一次波速度的数值、所述多次波速度的数值或一次波速度的数值与多次波速度的数值之间的数值；所述第二地震数据包括多个维度，所述多个维度包括炮检距和时间；第二处理单元，第二处理单元用于对所述第二地震数据进行预设处理，得到第三地震数据；所述预设处理包括，根据指定规则缩小所述第二地震数据中的所述炮检距，和/或，根据指定规则增大所述第二地震数据中的所述时间。
[0007] 为实现上述目的，本说明书实施方式还提供一种电子设备，包括：输入设备，存储器，处理器；所述输入设备用于获取第一地震数据，其中，所述第一地震数据包括一次波与多次波；所述处理器用于根据指定速度对所述地震数据动校正处理，获得第二地震数据；其中，所述指定速度包括所述一次波速度的数值、所述多次波速度的数值或一次波速度的数值与多次波速度的数值之间的数值；所述第二地震数据包括多个维度，所述多个维度包括炮检距和时间；对所述第二地震数据进行预设处理，得到第三地震数据；所述预设处理包括，根据指定规则缩小所述第二地震数据中的所述炮检距，和/或，根据指定规则增大所述第二地震数据中的所述时间；所述存储器用于存取数据，所述数据包括所述输入或被获取的地震数据，所述处理器处理的中间数据以及结果数据。
[0008] 为实现上述目的，本说明书实施方式还提供一种计算机存储介质，所述计算机存储介质存储有计算机程序指令，在所述计算机程序指令被执行时实现：获取第一地震数据，其中，所述第一地震数据包括一次波与多次波；根据指定速度对所述地震数据动校正处理，获得第二地震数据；其中，所述指定速度包括所述一次波速度的数值、所述多次波速度的数值或一次波速度的数值与多次波速度的数值之间的数值；所述第二地震数据包括多个维度，所述多个维度包括炮检距和时间；对所述第二地震数据进行预设处理，得到第三地震数据；所述预设处理包括，根据指定规则缩小所述第二地震数据中的所述炮检距，和/或，根据指定规则增大所述第二地震数据中的所述时间。
[0009] 为实现上述目的，本说明书实施方式还提供一种地震数据处理方法，所述方法包括：获取第一地震数据，其中，所述第一地震数据包括一次波与多次波；根据所述一次波的波速与所述多次波的波速，确定指定速度；所述指定速度介于所述一次波的波速与所述多次波的波速之间；根据所述指定速度对所述地震数据动校正处理，获得第二地震数据。在所述第二地震数据中，一次波与多次波大致分离。
[0010] 为实现上述目的，本说明书实施方式还提供一种地震数据处理装置，所述装置包括：获取单元，所述获取单元用于获取第一地震数据，其中，所述第一地震数据包括一次波与多次波；第一确定单元，所述第一确定单元用于根据所述一次波的波速与所述多次波的波速，确定指定速度，所述指定速度介于所述一次波的波速与所述多次波的波速之间；第二确定单元，所述第一确定单元用于根据所述指定速度对所述地震数据动校正处理，获得第二地震数据。
[0011] 由以上本申请实施方式提供的技术方案可见，本申请实施方式主要通过确定指定速度，所述指定速度介于所述一次波波速与所述多次波波速之间，来进行动校正，大致分离所述一次波与所述多次波，然后通过缩小所述第二地震数据中的炮检距，和/或，增大所述第二地震数据中的时间参数，进一步分离所述一次波与所述多次波，为后续压制多次波的方法提供条件。以上做法利用了在常规动校正基础上，利用指定速度对所述一次波过校正，对所述多次波欠校正，分离所述一次波与多次波，并通过炮检距或时间的改变，进一步分离多次波，可以实现对多次波更好的压制效果。附图说明
[0012] 为了更清楚地说明本说明书实施方式或现有技术中的技术方案，下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本说明书中记载的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
[0013] 图1是本说明书提供的一种地震数据处理的流程图；
[0014] 图2为本说明书提供的一个含两组一次波和多次波的地震数据示意图；
[0015] 图3为本说明书提供的一个对地震数据进行动校正后的地震数据示意图；
[0016] 图4为本说明书提供的一个在指定坐标系的一次波与多次波示意图；
[0017] 图5为本说明书提供的另一个在指定坐标系的一次波与多次波示意图；
[0018] 图6为本说明书提供的一个反动校正后的地震数据示意图；
[0019] 图7为本说明书提供的另一个反动校正后的地震数据示意图；
[0020] 图8为本说明书提供的一个实际情况中共反射点道集(CRP common reflection point)表征图；
[0021] 图9为本说明书提供的一个实际情况中对地震数据进行动校正后的地震数据表征图；
[0022] 图10为本说明书提供的一个实际情况中压制多次波后的地震数据表征图；
[0023] 图11为本说明书提供的一个实际情况中反动校正后地震数据的表征图；
[0024] 图12为本说明书提供的一个对多次波压制前数据计算的速度谱；
[0025] 图13为本说明书提供的一个对多次波压制后数据计算的速度谱；
[0026] 图14为本说明书提供的一个直接进行叠前偏移结果的表征图；
[0027] 图15是本说明书提供的一个压制多次波后叠前偏移结果的表征图；
[0028] 图16是本说明书提供的一种电子设备示意图。
具体实施方式
[0029] 下面将结合本说明书实施方式中的附图，对本说明书实施方式中的技术方案进行清楚、完整地描述，显然，所描述的实施方式仅仅是本说明书一部分实施方式，而不是全部的实施方式。基于本说明书中的实施方式，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施方式，都应当属于本申请保护的范围。
[0030] 请参阅图1。本说明书提供的一种地震数据处理方法的流程图。所述方法可以包括以下步骤。
[0031] 在本实施方式中，执行所述地震数据处理的客体可以是具有逻辑运算功能的电子设备。所述电子设备可以是服务器和客户端。所述客户端可以为台式电脑、平板电脑、笔记本电脑、工作站等。当然，客户端并不限于上述具有一定实体的电子设备，其还可以为运行于上述电子设备中的软体。还可以是一种通过程序开发形成的程序软件，该程序软件可以运行于上述电子设备中。
[0032] 步骤S10：获取第一地震数据，其中，所述第一地震数据包括一次波与多次波。
[0033] 在本实施方式中，所述地震数据是指地震勘探中通过各种传感器采集的原始地震数据，也可以是中间处理过的地震数据，例如直接从地震勘探获得的资料中获取。可以通过对地震数据进行处理和解释，可以推断地下岩层的性质和形态。所述地震数据包括一次波与多次波。所述一次波表示从震源出发，到达接收点(检波点)时，在低下界面之间发生了一次反射的波，所述多次波表示从震源出发，到达接收点(检波点)时，在地下界面之间发生了多次反射或折射的波。其中，所述一次波的速度大于所述多次波的速度。
[0034] 在本实施方式中，获取所述第一地震数据，可以通过输入设备的输入，所述输入设备可以是键盘、鼠标、语音输入设备等；也可以通过外接存储设备的输入，所述外接存储设备可以是U盘、机械硬盘等；也可以通过网络接收，例如因特网、局域网；也可以通过读取本地数据等。
[0035] 步骤S12：根据指定速度对所述地震数据动校正处理，获得第二地震数据；其中，所述指定速度包括所述一次波速度的数值、所述多次波速度的数值或一次波速度的数值与多次波速度的数值之间的数值；所述第二地震数据包括多个维度，所述多个维度包括炮检距和时间。
[0036] 在本实施方式中，所述指定速度可以是所述一次波的波速或所述多次波的波速或两者之间的速度。优选的，所述指定速度可以是所述一次波波速与所述多次波波速的平均速度。当指定速度为所述一次波波速与所述多次波波速的平均速度时，动校正处理后，更有利于所述一次波与所述多次波的分离。
[0037] 在本实施方式中，请参阅图2，所述第一地震数据处于所述第一坐标系，时间—炮检距坐标系。所述第一地震数据包括，所述一次波的接收时间、所述一次波能量团、所述一次波的跑检距，所述多次波的接收时间、所述多次波能量团、所述多次波的跑检距等。所述动校正是对第一地震数据中的时间进行时差校正，所以所述第二地震数据依旧处于时间—炮检距坐标系，所述第二地震数据包括多个维度，所述多个维度包括炮检距和时间。
[0038] 在本实施方式中，所述动校正是指把炮检距不同的各道上来自同一界面、同一点的反射波到达时间(接收时间)，校正为共中心点处的回声时间(接收时间)。即常规动校正将不同炮检距的反射时间校正到零炮检距的反射时间的过程，所述根据所述指定速度对所述地震数据动校正处理是将所述指定速度替换所述常规动校正中的速度进行动校正。
[0039] 步骤S14：对所述第二地震数据进行预设处理，得到第三地震数据；所述预设处理包括，根据指定规则缩小所述第二地震数据中的所述炮检距，和/或，根据指定规则增大所述第二地震数据中的所述时间。
[0040] 在本实施方式中，所述预设处理是指放大所述一次波与所述多次波差异的处理。可以增大所述一次波与所述多次波的分离度，为后期压制多次波提供良好基础。
[0041] 在本实施方式中，所述预设处理包括根据指定规则缩小所述第二地震数据中的所述炮检距，和/或，根据指定规则增大所述第二地震数据中的所述时间。所述指定规则可以使按一定比例放大缩小。例如同时将第二地震数据中的炮检距乘以0.3或0.5，当然也可以通过指数方式缩小炮检距或时间，但是，放的太大，曲率范围就非常大，一者会增加计算量，二者容易产生假频所以只要放大到能够分开的程度就可以，并不是越大越好，优选的是按照一定比例进行放大和缩小。在一个具体的实施方式中，所述缩小所述第二地震数据中的所述炮检距为按一定比例缩小所述第二地震数据中的所述炮检距。在另一个具体的实施方式中，增大所述第二地震数据中的所述时间为按一定比例增大所述第二地震数据中的所述时间。在另一个具体的实施方式中，根据指定规则缩小所述第二地震数据中的所述炮检距，和，根据指定规则增大所述第二地震数据中的所述时间同时进行。通过所述预设处理，增大所述一次波与所述多次波的曲率的绝对值，进一步分离所述一次波与所述多次波。下面结合图3具体说明，图3中，所述一次波曲线上扬，曲率为负，所述多次波曲线向下，曲率为正。通过缩小图3坐标系中对应点的跑检距，则曲线向坐标原点压缩，所述一次波的曲率图形特征参数与所述多次波的曲率图形特征参数绝对值增大，进一步分离所述一次波与所述多次波。同理，通过增大图3坐标系中的时间，所述一次波的曲率图形特征参数与所述多次波的曲率图形特征参数绝对值也可以增大，或者将两种预设处理的方式结合使用，都能达到增大所述一次波的曲率图形特征参数与所述多次波的曲率图形特征参数的绝对值的效果。
[0042] 在一个具体的场景示例中，本说明书实施方式提供的地震数据处理方法，可以通过运行于电子设备中的软件实现。
[0043] 在本场景示例中，获取第一地震数据，其中，所述第一地震数据包括一次波与多次波；可以通过输入设备的输入，也可以通过外接存储设备的输入，也可以通过网络接收，也可以通过读取本地数据等获取第一地震数据。根据指定速度对所述地震数据动校正处理，获得第二地震数据；其中，所述指定速度包括所述一次波速度的数值、所述多次波速度的数值或一次波速度的数值与多次波速度的数值之间的数值。动校正是在常规动校正的基础上用指定速度代替常规动校正速度来实施的动校正。对所述第二地震数据进行预设处理，得到第三地震数据；所述预设处理包括，根据指定规则缩小所述第二地震数据中的所述炮检距，和/或，根据指定规则增大所述第二地震数据中的所述时间。例如，将第二的证数据中的炮检距都乘以0.5，或第二地震数据中的时间都乘以2。所述预设处理可以进一步分离所述一次波与所述多次波的曲率，便于所述多次波的滤除。以上做法利用了所述一次波与所述多次波具有不同的速度，并利用了动校正原来进行初步分离，又通过预设处理，即利用所述一次波与所述多次波在坐标系的两个维度的性质，处理放大所述一次波与所述多次波的分离程度。
[0044] 在一个实施方式中，根据所述指定速度对所述地震数据动校正处理，获得所述第二地震数据，包括：根据所述指定速度，获得不同炮检距处检波点对应的动校正时差数据；将不同炮检距处检波点接收的旅行时间与所述动校正时差数据进行校正，获得所述第二地震数据；其中，所述旅行时间是地震波产生后至所述检波点所需时间。
[0045] 在本实施方式中，所述旅行时间是指反射波到达时间，也是检波器的接收时间。
[0046] 在本实施方式中，常规水平动校正获得的动校正时差根据如下公式获得：
[0047]
[0048] 其中，所述Δt为动校正时差，x为炮检距，V为震源在介质中的传播速度，即表述对应的一次波波速或对应的多次波波速，t0为自激自收点的接收时间即零炮检距处的反射时间，t为对应炮检距处的接收时间。所述根据所述指定速度对所述地震数据动校正处理及将所述指定速度替换所述V。在其他情况下的动校正采用同样的方式或原理，此处不再复述。当所述指定速度小于一次波波速，大于多次波波速时，所述一次波Δt动校正时差会增大，所述多次波Δt动校正时差会减小，一次波是过校正的，而多次波是欠校正的。下面结合图2和图3具体说明。所述图2表示含有两组一次波与多次波的模型数据，每一组中，所述一次波与所述多次波接近，不利于压制多次波。根据所述指定速度对图2中的数据进行动校正处理，输出结果如图3所示，所述一次波是过校正的，而多次波是欠校正的。当指定速度为一次波波速时，所述一次波校正正常，所述多次波欠校正；当指定速度为多次波波速时，所述一次波过校正，方法原理同采用中间速度的动校正处理，此处不再复述。
[0049] 在一个实施方式中，根据所述指定规则缩小所述第二地震数据中的所述炮检距，和/或，根据所述指定规则增大所述第二地震数据中的所述时间包括：按比例缩小所述第二地震数据中的所述炮检距，和/或，按比例增大所述第二地震数据中的所述时间。
[0050] 在本实施方式中，所述按比例缩小所述第二地震数据中的所述炮检距，和/或，按比例增大所述第二地震数据中的所述时间，例如将第二地震数据中的炮检距乘以0.3或0.5等或将时间乘以2或3等的方式，进行放大曲率，也可以按比例缩小所述第二地震数据中的所述炮检距和按比例增大所述第二地震数据中的所述时间同时进行的方式，进行放大曲率。本方法不至于将曲率图形特征参数放太大，避免计算量过大和产生假频。
[0051] 在一个实施方式中，所述第三地震数据至少符合第一特征，所述第一特征为所述一次波的曲率图形特征参数与所述多次波的曲率图形特征参数不处于同一正数区间或负数区间；在获得所述第三地震数据后，所述方法还包括：根据所述第一特征，对所述第三地震数据进行压制多次波处理，获得压制多次波后的第三地震数据。
[0052] 参阅图3，在本实施方式中，当所述指定速度采用所述一次波与所述多次波之间的中间速度进行动校正时，所述一次波过校正，所述多次波欠校正，所述一次波曲线上扬，所述多次波曲线向下。所述一次波曲率为负数，所述多次波曲率为正数。同理，当所述指定速度采用所述一次波的速度时，所述多次波欠校正，所述多次波的曲率图形特征参数处于正数区间，所述一次波的曲率图形特征参数不处于正数区间；当所述指定速度采用所述多次波的速度时，所述一次波过校正，所述一次波的曲率图形特征参数处于负数区间，所述多次波的曲率图形特征参数不处于负数区间。因此所述一次波的曲率图形特征参数与所述多次波的曲率图形特征参数不处于同一正数区间或负数区间。
[0053] 在本实施方式中，在获得所述第三地震数据后，所述方法还包括：根据所述第一特征，对所述第三地震数据进行压制多次波处理，获得压制多次波后的第三地震数据。
[0054] 在本实施方式中，通过所述一次波与所述多次波的曲率图形特征参数不处于同一区间，对其中分离出来的多次波进行滤波处理。在一个具体实施例中，可以通过Radon变换，在Radon域中对所述多次波进行滤波。
[0055] 在一个场景示例中，获得压制多次波后的第三地震数据可以包括以下步骤。
[0056] 步骤S20：将所述第三地震数据映射至指定坐标系，获得第四地震数据；所述指定坐标系包括曲率图形特征维度和时间维度。在一个具体的实施例中，对所述第三地震数据根据如下公式进行抛物Radon变换。
[0057] m(τ，q)＝Σd(x，t＝τ+qx2)
[0058] 其中，m(τ,q)为Radon域(τ,q)中的数据，d(x,t)是炮检距—时间域的数据，此处即为所述第一坐标系中的所述第三地震数据，τ为截距时间，q为曲率参数。炮检距—时间域的同相轴在Radon域的曲率为 其中Δt＝tmax-τ，tmax为同相轴在最大炮检距处的时间。下面结合图5具体说明，图5是在对所述第二地震数据进行预设处理获得的所述第三地震数据进行Radon变换，如图所示，所述指定坐标系的横坐标代表曲率，纵坐标代表时间。所示一次波曲率图形特征参数为负，经过Radon变换后，所述一次波的能量团在负数区间，所示多次波曲率图形特征参数为正，经过Radon变换后，所述多次波的能量团在正数区间，便于对多次波的滤波处理。
[0059] 步骤S22：对所述第四地震数据进行滤波处理，获得压制多次波后的第四地震数据。请参阅图5，在一个具体的实施例中，所述一次波能量团在负数区间，所述多次波能量团在正数区间，对正数区间的能量团进行滤波，获得较为纯净的一次波能量团。
[0060] 步骤S24：将所述压制多次波后的第四地震数据反映射处理，获得所述压制多次波后的第三地震数据。具体的，对所述第四地震数据根据如下公式进行反Radon变换。
[0061] d(x,t)＝∑m(q,τ＝t-qx2)
[0062] 其中，m(τ,q)为Radon域(τ,q)中的数据，d(x,t)是炮检距—时间域的数据，即所述第一坐标系中的所述第二数据体，τ为截距时间，q为曲率参数。炮检距—时间域的同相轴在Radon域的曲率为 其中Δt＝tmax-τ，tmax为同相轴在最大炮检距处的时间。
[0063] 在一个实施方式中，在获得压制多次波后的第三地震数据后，还包括：根据第二指定规则，放大所述压制多次波后的第三地震数据后中的所述炮检距，和/或，根据第二指定规则，缩小所述压制多次波后的第三地震数据后中的所述时间，获得压制多次波后的第二地震数据；其中，所述第二指定规则为对应于所述指定规则的反变换。例如，当所述预设处理的方式为将所述第一数据体的炮检距缩小为原来的三分之一，则此时，将所述压制多次波后的第二数据体的炮检距放大三倍，恢复至炮检距。同理，如果所述预设处理的方式为放大时间，则所述指定规则为缩小对应的时间。通过指定规则的处理，可以使所述第二数据体中的数据恢复至实际情况下的炮检距和时间。
[0064] 在一个实施方式中，在获得所述压制多次波后的第二地震数据后还包括：对所述压制多次波后的第二地震数据进行反动校正处理，获得压制多次波后的第一地震数据。其中，所述反动校正处理对应于所述动校正处理，即按照所述指定速度将不同炮检距的反射时间从零炮检距的反射时间恢复至不同炮检距各自对应的反射时间的过程，以方便后期对地震数据的处理。
[0065] 下面结合图2至图7具体说明本说明书的一个具体的场景示例。本说明书实施方式提供的地震数据处理方法，可以通过运行于电子设备中的软件实现。该地震数据处理方法使得一次波与多次波更好区分，从而获得更好的多次波滤波效果，方便后续作业。
[0066] 在本实施场景中，获取第一地震数据，其中，所述第一地震数据包括一次波与多次波。参阅图2，获取两组一次波与多次波的地震数据示意图。获取所述第一地震数据，可以通过输入设备的输入，所述输入设备可以是键盘、鼠标、语音输入设备等；也可以通过外接存储设备的输入，所述外接存储设备可以是U盘、机械硬盘等；也可以通过网络接收，例如因特网、局域网；也可以通过读取本地数据等。
[0067] 在本实施场景中，根据指定速度对所述地震数据动校正处理，获得第二地震数据；其中，所述指定速度包括所述一次波速度的数值、所述多次波速度的数值或一次波速度的数值与多次波速度的数值之间的数值；所述第二地震数据包括多个维度，所述多个维度包括炮检距和时间。所述动校正是指把炮检距不同的各道上来自同一界面、同一点的反射波到达时间(接收时间)，校正为共中心点处的回声时间(接收时间)。即常规动校正将不同炮检距的反射时间校正到零炮检距的反射时间的过程，所述根据所述指定速度对所述地震数据动校正处理是将所述指定速度替换所述常规动校正中的速度进行动校正。参阅图3，通过确定的所述指定速度对图2中所述地震数据进行动校正，在这里，所述指定速度为所述一次波速度和所述多次波速度的平均速度，动校正后所述一次波过校正，曲率为负，所述多次波欠校正，曲率为正。
[0068] 在本实施场景中，对所述第二地震数据进行预设处理，得到第三地震数据；所述预设处理包括，根据指定规则缩小所述第二地震数据中的所述炮检距，和/或，根据指定规则增大所述第二地震数据中的所述时间。参阅图5，对图3中数据进行预设处理，所述预设处理为缩小炮检距和/或增大时间，并进行高分辨Radon变换。具体的，同时将第二地震数据中的炮检距乘以0.3或0.5。通过所述预设处理，增大所述一次波与所述多次波的曲率的绝对值，进一步分离所述一次波与所述多次波。通过缩小图3坐标系中对应点的跑检距，则曲线向坐标原点压缩，所述一次波的曲率图形特征参数与所述多次波的曲率图形特征参数绝对值增大，进一步分离所述一次波与所述多次波。同理，通过增大图3坐标系中的时间，所述一次波的曲率图形特征参数与所述多次波的曲率图形特征参数绝对值也可以增大，或者将两种预设处理的方式结合使用，都能达到增大所述一次波的曲率图形特征参数与所述多次波的曲率图形特征参数的绝对值的效果。对比图5与图4，所述一次波能量团与所述多次波能量团分离效果更好。参阅图7，对图5Radon域中的多次波进行滤波处理，然后进行Radon反变换，根据对应于所述预设处理的指定规则，放大炮检距和/或减小时间进行恢复炮检距和/或时间，然后进行了反动校正处理，获得压制多次波后的地震数据。对比图7与图6，多次波压制效果优于图6的多次波压制效果。
[0069] 下面结合图8至图14，具体本说明书实施方式对实际数据的处理效果。参阅图8，图8表示具体的共反射点道集(CRP common reflection point)，当然，本实施方式也可以运用于共中心点道集(CMP common middle point)等道集。对图8所示的CRP道集进行动校正，得到图9所示的道集，可以看出存在欠校正的多次波信息。对所述第二地震数据进行预设处理，得到第三地震数据。对所述第三地震数据进行压制多次波处理，获得压制多次波后的第三地震数据。对所述压制多次波后的第三地震数据进行处理，获得压制多次波后的第二地震数据，获得图10所示的压制多次波后的表征图。参阅图11，是对图10中的数据进行反动校正，获得了压制了多次波后地震数据表征图。对比图11与图8，可以看出，所述多次波被很好压制。参阅图12、图13，图12是对应于图8所示的速度谱，图13是对应于图11的速度谱，图中虚线表示所述一次波波速的拾取线，拾取线左下方表示多次波能量，图13对比于图12左下方的多次波能量被显著滤除。参阅图14、图15，图14是对多次波压制前的数据即图8中的数据，直接进行叠前时间偏移的结果，图15是对多次波压制后的数据即图11中的数据，进行叠前时间偏移的结果，从椭圆形和矩形框内可以看出，多次波压制后，偏移剖面中的同相轴连续性变好，分辨率也得到提高。例如，矩形框内的同相轴，多次波压制前，两个同相轴存在一定程度的干涉；多次波压制后，这两个同相轴得到了较好地分离。
[0070] 本说明书实施方式中还提供了一种地震数据处理装置，如下面的实施方式所述。由于一种地震数据处理装置解决问题的原理与一种地震数据处理方法相似，因此一种地震数据处理装置的实施可以参见一种地震数据处理方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。该装置具体可以包括：获取单元、第一确定单元、第一处理单元、第二处理单元。下面对该结构进行具体说明。
[0071] 获取单元，所述获取单元用于获取第一地震数据，其中，所述第一地震数据包括一次波与多次波。
[0072] 在本实施方式中，所述地震数据包括一次波速度与多次波速度。其中，所述地震数据可以从地震勘探获得的资料中获取；所述一次波表示从震源出发，到达接收点时，在低下界面之间发生了一次反射的波，所述多次波表示从震源出发，到达接收点时，在地下界面之间发生了多次反射或折射的波。
[0073] 第一处理单元，第一处理单元用于根据指定速度对所述地震数据动校正处理，获得第二地震数据；其中，所述指定速度包括所述一次波速度的数值、所述多次波速度的数值或一次波速度的数值与多次波速度的数值之间的数值；所述第二地震数据包括多个维度，所述多个维度包括炮检距和时间。
[0074] 在本实施方式中，常规动校正是将不同炮检距的反射时间校正到零炮检距的反射时间的过程，所述根据所述指定速度对所述地震数据动校正处理是将所述指定速度替换所述常规动校正中的速度。
[0075] 第二处理单元，第二处理单元用于对所述第二地震数据进行预设处理，得到第三地震数据；所述预设处理包括，根据指定规则缩小所述第二地震数据中的所述炮检距，和/或，根据指定规则增大所述第二地震数据中的所述时间。
[0076] 在一个本实施方式中，在第二处理单元后，所述装置还包括：第三处理单元，用于获得压制多次波后的第三地震数据。所述第三地震数据至少符合第一特征，所述第一特征为所述一次波的曲率图形特征参数与所述多次波的曲率图形特征参数不处于同一正数区间或负数区间；在获得所述第三地震数据后，所述方法还包括：根据所述第一特征，对所述第三地震数据进行压制多次波处理，获得压制多次波后的第三地震数据。
[0077] 在一个本实施方式中，在第三处理单元后，还包括：第四处理单元。所述第四处理单元根据第二指定规则，放大所述压制多次波后的第三地震数据后中的所述炮检距，和/或，根据第二指定规则，缩小所述压制多次波后的第三地震数据后中的所述时间，获得压制多次波后的第二地震数据；其中，所述第二指定规则为对应于所述指定规则的反变换。
[0078] 在一个本实施方式中，在第四处理单元后，还包括：第五处理单元；所述第五处理单元对所述压制多次波后的第二地震数据进行反动校正处理，获得压制多次波后的第一地震数据。
[0079] 请参阅图16。本说明书实施方式中还提供了一种电子设备，包括：输入设备，存储器，处理器；所述输入设备用于获取第一地震数据，其中，所述第一地震数据包括一次波与多次波；所述处理器用于根据指定速度对所述地震数据动校正处理，获得第二地震数据；其中，所述指定速度包括所述一次波速度的数值、所述多次波速度的数值或一次波速度的数值与多次波速度的数值之间的数值；所述第二地震数据包括多个维度，所述多个维度包括炮检距和时间；对所述第二地震数据进行预设处理，得到第三地震数据；所述预设处理包括，根据指定规则缩小所述第二地震数据中的所述炮检距，和/或，根据指定规则增大所述第二地震数据中的所述时间；所述存储器用于存取数据，所述数据包括所述输入或被获取的地震数据，所述处理器处理的中间数据以及结果数据。
[0080] 所述输入设备是用户和计算机系统之间进行信息交换的主要装置之一。所述输入设备包括键盘、鼠标、摄像头、扫描仪、光笔、手写输入板、语音输入装置等；输入设备用于把原始数据和处理这些数的程序输入到计算机中。所述输入设备还可以获取接收其他模块、单元、设备传输过来的数据。
[0081] 所述处理器可以按任何适当的方式实现。例如，处理器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific Integrated Circuit，ASIC)、可编程逻辑控制器和嵌入微控制器的形式等等。
[0082] 所述存储器现代信息技术中用于保存信息的记忆设备。所述存储器包括很多层次，在数字系统中，只要能保存二进制数据的都可以是存储器；在集成电路中，一个没有实物形式的具有存储功能的电路也叫存储器，如RAM、FIFO等；在系统中，具有实物形式的存储设备也叫存储器，如内存条、TF卡等。
[0083] 在本实施方式中，该电子设备具体实现的功能和效果，可以与其它实施方式对照解释，在此不再赘述。
[0084] 本说明书实施方式中还提供了一种计算机存储介质，所述计算机存储介质存储有计算机程序指令，在所述计算机程序指令被执行时实现：获取第一地震数据，其中，所述第一地震数据包括一次波与多次波；根据指定速度对所述地震数据动校正处理，获得第二地震数据；其中，所述指定速度包括所述一次波速度的数值、所述多次波速度的数值或一次波速度的数值与多次波速度的数值之间的数值；所述第二地震数据包括多个维度，所述多个维度包括炮检距和时间；对所述第二地震数据进行预设处理，得到第三地震数据；所述预设处理包括，根据指定规则缩小所述第二地震数据中的所述炮检距，和/或，根据指定规则增大所述第二地震数据中的所述时间。
[0085] 在本实施方式中，所述存储器包括但不限于随机存取存储器(Random Access Memory,RAM)、只读存储器(Read-Only Memory,ROM)、缓存(Cache)、硬盘(Hard Disk Drive,HDD)或者存储卡(Memory Card)。所述存储器可以用于存储计算机程序指令。网络通信单元可以是依照通信协议规定的标准设置的，用于进行网络连接通信的接口。
[0086] 在本实施方式中，该计算机存储介质存储的程序指令具体实现的功能和效果，可以与其它实施方式对照解释，在此不再赘述。
[0087] 本说明书还提供的另一种地震数据处理方法，所述方法包括：获取第一地震数据，其中，所述第一地震数据包括一次波与多次波；根据所述一次波的波速与所述多次波的波速，确定指定速度；所述指定速度介于所述一次波的波速与所述多次波的波速之间；根据所述指定速度对所述地震数据动校正处理，获得第二地震数据。
[0088] 在一个实施方式中，所述第二地震数据至少符合第一特征，所述第一特征为所述一次波的曲率图形特征参数与所述多次波的曲率图形特征参数不处于同一正数区间或负数区间；在获得所述第二地震数据后，所述方法还包括：根据所述第一特征，对所述第二地震数据进行压制多次波处理，获得压制多次波后的第二地震数据。
[0089] 在一个实施方式中，在获得所述压制多次波后的第二地震数据后还包括：对所述压制多次波后的第二地震数据进行反动校正处理，获得压制多次波后的第一地震数据。
[0090] 下面结合图2、图3、图4、图6具体说明本说明书的一个实施例。参阅图2，获取两组一次波与多次波的地震数据示意图。参阅图3，通过确定的所述指定速度对图2中所述地震数据进行动校正，所述一次波过校正，曲率为负，所述多次波欠校正，曲率为正。参阅图4，图4所示为对图3中数据进行高分辨Radon变换的结果，如图中所述，多次波与一次波相分离。
参阅图6，对图4Radon域中的多次波进行滤波处理，然后进行Radon反变换，并进行了反动校正处理，获得了压制多次波后的地震数据，从中可以看出，所述多次波被显著压制。
参阅图6，对图4Radon域中的多次波进行滤波处理，然后进行Radon反变换，并进行了反动校正处理，获得了压制多次波后的地震数据，从中可以看出，所述多次波被显著压制。
[0091] 参照本说明书提供的第一种地震数据处理方法，本说明书提供的另一种地震数据处理方法无需经过预设处理和对应的指定规则，方法更简洁，但是所述多次波的压制效果不如经过预设处理获得的压制多次波后的地震数据，近偏移距的多次波仍有残留。
[0092] 本说明书实施方式提出地震数据处理方法最本质的特点是，在一种实施方式中，采用所述一次波的波速与所述多次波的波速中间速度进行动校正，分离所述一次波与多次波。在另一种实施方式中，通过缩小所述炮检距参数，和/或，增大所述时间参数，增大一次波与多次波的分离程度。在一种优选的实施例中，同时采用所述一次波的波速与所述多次波的波速中间速度进行动校正与通过缩小所述炮检距参数，和/或，增大所述时间参数，对所述一次波与多次波进行分离，为地震数据中多次波的压制提供改进方法。
[0093] 本说明书还提供一种地震数据处理装置，所述装置包括：获取单元，所述获取单元用于获取第一地震数据，其中，所述第一地震数据包括一次波与多次波；第一确定单元，所述第一确定单元用于根据所述一次波的波速与所述多次波的波速，确定指定速度，所述指定速度介于所述一次波的波速与所述多次波的波速之间；第二确定单元，所述第一确定单元用于根据所述指定速度对所述地震数据动校正处理，获得第二地震数据。
[0094] 尽管本申请内容中提到一种地震数据处理方法、装置、电子设备或计算机存储介质，但是，本申请并不局限于必须是行业标准或实施例所描述的情况等，某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取、处理、输出、判断方式等的实施例，仍然可以属于本申请的可选实施方案范围之内。
[0095] 虽然本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。
[0096] 上述实施例阐明的装置或模块等，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述模块的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个模块或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。
[0097] 本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
[0098] 本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构、类等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
[0099] 通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，移动终端，服务器，或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。
[0100] 本说明书中的各个实施例采用递进的方式描述，各个实施例之间相同或相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。本申请可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 一种地震波形特征分解的方法及装置 | 2020-05-12 | 40 |
| 一种地震波检波器测试系统及方法 | 2020-05-12 | 745 |
| 一种高分辨率地震波成像方法和装置 | 2020-05-13 | 962 |
| 一种高分辨率地震波成像方法和装置 | 2020-05-13 | 644 |
| 基于地震波的地下入侵监测系统 | 2020-05-13 | 467 |
| 一种利用斜井直达地震波走时计算地层地震波速度的方法 | 2020-05-12 | 481 |
| 地震波传播展示方法 | 2020-05-11 | 575 |
| 地震波传播速度获取方法、测深方法及地震波观测系统 | 2020-05-12 | 883 |
| 勘测三维地震波的方法及设备 | 2020-05-12 | 120 |
| VTI介质中地震波的走时计算方法 | 2020-05-13 | 76 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
分析报告
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
QQ群二维码
意见反馈
意见反馈