一种地震子波的确定方法、装置、电子设备以及存储介质

技术领域

本发明涉及海洋地震技术领域，尤其涉及一种地震子波的确定方法、装置、电子设备以及存储介质。

背景技术

地震子波是地震信号处理中的一项重要参数，贯穿了反褶积、子波整形、多次波预测、鬼波压制、地震正演模拟、偏移成像、全波形反演、叠前叠后属性反演、宽频处理等处理过程。

目前获得地震子波的方式主要是自相关法提取地震子波，在反射系数为白噪声的情况下，地震数据自相关近似等于地震子波自相关，因而，可以直接根据地震数据的振幅谱得到地震子波的振幅谱，但是，对于地震子波的相位谱，通常只能假设为零相位或者最小相位，这样与地震子波实际为混合相位的事实相违背，导致合成的地震子波与实际地震子波之间存在差异。

发明内容

本发明提供了一种地震子波的确定方法、装置、电子设备以及存储介质，通过对非零偏移距直达波进行偏移距滤波，获得零偏移距直达波，其相位谱近似等于合成子波的相位谱，再结合子波振幅谱，实现快速、高效、准确的合成地震子波。

第一方面，本发明实施例提供了一种地震子波的确定方法，该方法包括：

根据地震数据的振幅谱，确定待确定地震子波的振幅谱；

对地震直达波数据进行偏移距滤波处理，得到零偏移距直达波；

根据所述零偏移距直达波的相位谱，确定待确定地震子波的相位谱；

根据待确定地震子波的振幅谱和相位谱，确定目标地震子波。

第二方面，本发明实施例还提供了一种地震子波的确定装置，该装置包括：

振幅谱确定模块，用于根据地震数据的振幅谱，确定待确定地震子波的振幅谱；

零偏移直达波确定模块，用于对地震直达波数据进行偏移距滤波处理，得到零偏移距直达波；

相位谱确定模块，用于根据所述零偏移距直达波的相位谱，确定待确定地震子波的相位谱；

地震子波确定模块，用于根据待确定地震子波的振幅谱和相位谱，确定目标地震子波。

第三方面，本发明实施例还提供了一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如本发明实施例中任一所述的地震子波的确定方法。

第四方面，本发明实施例还提供了一种存储计算机可执行指令的存储介质，所述计算机可执行指令在由计算机处理器执行时用于执行如本发明实施例中任一所述的地震子波的确定方法。

本发明实施例的技术方案，通过对非零偏移距直达波进行偏移距滤波，获得零偏移距直达波，其相位谱近似等于合成子波的相位谱，再结合子波振幅谱，实现快速、高效、准确的合成地震子波。

应当理解，本部分所描述的内容并非旨在标识本发明的实施例的关键或重要特征，也不用于限制本发明的范围。本发明的其它特征将通过以下的说明书而变得容易理解。

附图说明

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

图1是本发明实施例一提供的一种地震子波的确定方法的流程图；

图2是本发明实施例一提供的一种基于相位替换的海洋地震子波合成技术流程图；

图3是本发明实施例二提供的另一种地震子波的确定方法的流程图；

图4是本发明实施例二提供的一个实际地震数据二维测线的局部偏移剖面图；

图5(a)是本发明实施例二提供的实际地震数据二维测线的局部偏移剖面的振幅谱；

图5(b)是本发明实施例二提供的实际地震数据二维测线的局部偏移剖面的相位谱；

图5(c)是本发明实施例二提供的零偏移距滤波获得的地震子波的振幅谱；

图5(d)本发明实施例二提供的零偏移距滤波获得的地震子波的相位谱；

图5(e)是本发明实施例二提供的测线零偏移距滤波获得的地震子波时域信号谱；

图5(f)是本发明实施例二提供的采用图5(a)和图5(d)反傅里叶变换计算的合成子波；

图5(g)是本发明实施例二提供的对图5(f)预测反褶积气泡压制后的期望输出曲线图；

图5(h)本发明实施例二提供的对图5(f)预测反褶积气泡压制后的滤波算子曲线图；

图6是本发明实施例二提供的一种将滤波因子应用到地震数据处理的效果图；

图7是本发明实施例三提供的一种地震子波的确定装置的结构示意图；

图8是本发明实施例四提供的一种地震子波的确定方法电子设备的结构示意图。

具体实施方式

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

需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。

实施例一

图1为本发明实施例一提供了一种地震子波的确定方法的流程图，本实施例可适用于地震子波的确定情况，该方法可以由地震子波的确定装置来执行，该地震子波的确定装置可以采用硬件和/或软件的形式实现，该地震子波的确定装置可配置于地震波探测系统中。

如图1所示，该方法包括：

S110、根据地震数据的振幅谱，确定待确定地震子波的振幅谱。

本申请实施例中，以实测地震数据作为已知信息，需要说明的是，自相关法提取地震子波在工业界得到广泛应用，是一种从反射地震数据本身提取地震子波的方法，基本原理是在反射系数为白噪的前提下，地震数据自相关近似等于地震子波自相关，进而直接由实测记录的地震数据的振幅谱得到地震子波的振幅谱。

作为一种可选但非限定的实现方式，根据地震数据的振幅谱，确定待确定地震子波的振幅谱，包括以下步骤A1-A2：

步骤A1、在确定反射系数为白噪声时，对地震数据进行傅里叶变换，得到地震数据的振幅谱。

步骤A2、将所述地震数据的振幅谱，作为待确定地震子波的振幅谱。

本申请实施例中，白噪声是指功率谱密度在整个频域内是常数的噪声，所有频率具有相同能量密度的随机噪声称为白噪声，它在频率谱上具有均匀分布的能量。地震数据信号通常是复杂的波形，通过傅里叶变换可以将其分解成不同频率的成分，得到地震数据的频谱图像，其中包括振幅谱。使用反射地震数据本身自相关提取地震子波方法，在反射系数为白噪的前提下，地震数据自相关近似等于地震子波自相关，进而直接由实测记录的地震数据的振幅谱作为待测地震子波的振幅谱。

在反射系数为白噪情况下，地震子波自相关r

假如地震子波是我们的目标合成子波，那么地震子波振幅谱等于地震数据振幅谱。地震数据是实测已知的，其地震数据振幅谱可直接对地震数据傅里叶变换获得，至此，可以获得合成子波的振幅谱。

S120、对地震直达波数据进行偏移距滤波处理，得到零偏移距直达波。

本申请实施例中，需要说明的是，地震数据的直达波是地震震源激发后由激发源直接传播到接收电缆(拖缆)被记录的波形，仅包含震源鬼波，偏移距越小越接近“远场子波”。因未曾在地下界面发生反射、折射，故只携带最浅层的速度信息，不携带地下地层界面的信息。即直达波对应的震源信号子波波形没有因在界面处发生反射、折射而产生形态改变，因此相对于反射波等其它地震波型来说，利用直达波提取震源信号子波具有一定的优势。此外，直达波传播路径最短，在多数地震道上通常最早出现，不受反射波、折射波等其它波形的干扰，相对容易识别。

对地震直达波数据进行偏移距滤波处理得到零偏移距的直达波，此时的零偏移距直达波可视为“远场子波”，叠加检波点鬼波，可视为地震子波。

S130、根据所述零偏移距直达波的相位谱，确定待确定地震子波的相位谱。

本申请实施例中，海洋地震直达波数据包含了地震子波信息，对其进行偏移距滤波获取零偏移距直达波，此时的零偏移距直达波叠加震源鬼波已十分接近地震子波，可以将零偏移距直达波的相位谱近似作为合成地震子波的相位谱。

本申请实施例中，海洋地震勘探中，涉及子波的概念较为复杂，为明确本申请中各类子波名称的确切含义，需要说明的是，“震源信号”通常指的是震源激发的不含鬼波的初始信号，“远场子波”通常在震源信号的基础上叠加了一个震源鬼波，而采集的数据中的“地震子波”则在远场子波的基础上又包含了检波点鬼波。实际地震数据含震源鬼波和检波点鬼波，“地震子波”也含两个鬼波。在上述中已说明对地震直达波数据进行偏移距滤波处理得到零偏移距的直达波，此时的零偏移距直达波可视为远场子波，叠加检波点鬼波，可视为地震子波，通过公式说明，即：

z

其中，这里将地震直达波数据记为z(t)是实测已知的，取其中一个地震道z

作为一种可选但非限定的实现方式，根据所述零偏移距直达波的相位谱，确定待确定地震子波的相位谱，包括以下步骤B1-B2：

步骤B1、若确定所述地震数据已进行检波点鬼波消除处理，则对所述零偏移距直达波进行傅里叶变换，得到所述零偏移距直达波的相位谱。

步骤B2、将所述零偏移距直达波的相位谱，作为待确定地震子波的相位谱。

如果实际地震数据做了检波点鬼波消除，则上述过程直接采用远场子波(也即零偏移距直达波)的相位谱代替合成子波相位谱，总之，保证合成的地震子波包含的鬼波效应与地震数据保持一致。

作为一种可选但非限定的实现方式，根据所述零偏移距直达波的相位谱，确定待确定地震子波的相位谱，还包括以下步骤C1-C2：

步骤C1、若确定所述地震数据未进行检波点鬼波消除处理，则确定检波点鬼波，并将所述零偏移距直达波和所述检波点鬼波进行叠加，获得合成子波。

步骤C2、对所述合成子波进行傅里叶变换，得到所述合成子波的相位谱。

步骤C3、将所述合成子波的相位谱，作为待确定地震子波的相位谱。

本申请实施例中，是否叠加检波点鬼波视实际地震数据情况而定。若实际地震数据没有进行检波点鬼波消除，则合成子波由零偏移距直达波和检波点鬼波叠加，根据公式b(t)＝z

S140、根据待确定地震子波的振幅谱和相位谱，确定目标地震子波。

本申请实施例中，从地震子波两要素振幅谱和相位谱为立足点，根据傅里叶变换解析公式，通过地震子波的振幅谱和相位谱就能合成地震子波。

作为一种可选但非限定的实现方式，根据待确定地震子波的振幅谱和相位谱，确定目标地震子波，包括以下步骤D1-D2：

步骤D1、根据待确定地震子波的振幅谱和相位谱，确定目标地震子波的复数谱。

步骤D2、对所述复数谱进行反傅里叶变换，得到目标地震子波。

本申请实施例中，地震子波就是目标合成子波，时间域地震子波由傅里叶变换得到地震子波复数谱可表示为：

其中w＝2πf表示圆频率，FFT表示傅里叶变换函数，b(t)表示地震子波(也即合成子波)，

同理，地震子波b(t)相应的地震数据x(t)的复数谱

其中

对合成子波的复数谱

参见图2是本申请基于相位替换的海洋地震子波合成技术流程图，地震数据自相关近似等于地震子波自相关，进而直接由地震数据记录的振幅谱得到子波的振幅谱，对非零偏移距直达波进行偏移距滤波，获得零偏移距直达波，其相位谱近似等于合成子波的相位谱，再结合子波振幅谱，合成子波的复数谱，对该复数谱反傅里叶变换就可计算出合成子波。需要进行说明的是，本实施例中还提供了一种直接根据理论模拟的地震子波的相位谱合成地震子波的示例。具体的，直接由地震数据记录的振幅谱得到子波的振幅谱，通过理论模型，根据震源类型以及地震参数等，模拟得到模拟地震子波，提取模拟地震子波的相位谱，再结合子波振幅谱，合成子波的复数谱，对该复数谱反傅里叶变换就可计算出合成子波。以上是本申请合成子波计算全过程，实际叠前、叠后地震数据均可适应，该合成子波可直接用于后续子波整形、地震波正反演等处理环节。

本发明公开了一种地震子波的确定方法，该方法包括：根据地震数据的振幅谱，确定待确定地震子波的振幅谱；对地震直达波数据进行偏移距滤波处理，得到零偏移距直达波；根据零偏移距直达波的相位谱，确定待确定地震子波的相位谱；根据待确定地震子波的振幅谱和相位谱，确定目标地震子波。本发明实施例的技术方案，通过对非零偏移距直达波进行偏移距滤波，获得零偏移距直达波，其相位谱近似等于合成子波的相位谱，再结合子波振幅谱，实现快速、高效、准确的合成地震子波。

实施例二

图3为本发明实施例二提供的另一种地震子波的确定方法的流程图，本发明实施例在上述实施例的基础上，对地震直达波数据进行偏移距滤波处理，得到零偏移距直达波进行了进一步的具体化，未在本申请实施例中详尽描述的方案见上述实施例。如图3所示，本发明实施例二提供了一种地震子波的确定方法的流程图，该方法包括：

S210、根据地震数据的振幅谱，确定待确定地震子波的振幅谱。

S220、确定地震直达波数据对应的震源位置，并根据所述震源位置确定滤波因子，所述滤波因子用于表示当前地震道的地震直达波数据与零偏移距地震道中的地震直达波数据之间的差异。

本申请实施例中，滤波因子为体现当前地震道的地震直达波与已知地震道中子波数据差异的干扰因素的综合因子，其中干扰因素包括：偏移距、震源水中深度、检波点沉放深度、水中声波速以及水面反射系数等。

在实际应用中，不同震源由于位置信息(位于水中深度、水中声波速度以及水面反射系数等条件)的不一致，形成的滤波因子也不尽相同，因此，需要根据震源位置信息确定对应的滤波因子。

作为一种可选但非限定的实现方式，所述滤波因子的确定过程包括以下步骤E1-E3：

步骤E1、根据当前地震道的地震直达波数据与震源鬼波的到达时差，模拟得到第一脉冲响应。

具体的，先构建待测偏移距地震道直达波的时距曲线方程，然后构建偏移距地震道的震源鬼波的时距曲线方程，基于待测偏移距地震道的直达波和待测偏移距地震道的震源鬼波时距曲线方程计算得到待测偏移距地震道中的第一到达时差，根据第一到达时差进行数值模拟，得到所述第一脉冲响应。

步骤E2、根据零偏移距地震道中的地震直达波数据与震源鬼波的到达时差，模拟得到第二脉冲响应。

具体的，先构建已知偏移距地震道的直达波的时距曲线方程，然后构建已知偏移距地震道的震源鬼波的时距曲线方程，基于已知偏移距地震道的直达波和震源鬼波的时距曲线方程计算得到已知偏移距地震道中的第二到达时差，根据第二到达时差进行数值模拟，得到所述第二脉冲响应。

步骤E3、将第二脉冲响应与第一脉冲响应的比值，作为滤波因子。

具体的，根据直达波与震源鬼波的时距方程，计算待测偏移距和已知偏移距地震道的直达波与震源鬼波第一到达时差和第二到达时差，根据该时差模拟待测偏移距和已知偏移距带震源鬼波的地震道第一脉冲响应和第二脉冲响应，滤波因子为第一脉冲响应和第二脉冲响应的比值。

S230、根据所述滤波因子对地震直达波数据进行偏移距滤波处理，得到零偏移距直达波。

本申请实施例中，实际地震直达波数据是震源激发后由震源直接传播到接收电缆被记录的波形，仅包含震源鬼波，偏移距越小越接近远场子波。对地震直达波数据进行零偏移距滤波处理，通过对已构建偏移距地震道直达波时距曲线方程和偏移距地震道震源鬼波时距曲线方程，令偏移距为0时(也即零偏移距直达波)，计算对应的零偏移距到达时差，对零偏移距到达时差进行数值模拟，得到零偏移距地震道脉冲响应，将第一脉冲响应与零偏移距地震道脉冲响应的比值作为对应的滤波因子，从而计算得到远场子波数据。

S240、根据所述零偏移距直达波的相位谱，确定待确定地震子波的相位谱。

S250、根据待确定地震子波的振幅谱和相位谱，确定目标地震子波。

为进一步验证方法的有效性，以某海洋二维拖缆地震数据为例，图4是该实际地震数据二维测线的局部偏移剖面。按照本申请的方案步骤进行处理，图5给出了合成地震子波的全过程。图5(a)和图5(b)分别是图4的振幅谱和相位谱，相位谱显示实际的地震数据是混合相位的，图5(a)即作为合成子波的振幅谱；图5(e)是该测线直达波数据零偏移距滤波获得的“地震子波”(零偏移距直达波叠加检波点鬼波)；图5(c)和图5(d)分别是图5(e)的振幅谱和相位谱，对比图5(c)和图5(a)发现二者高度一致(地震数据的叠加效应振幅谱更加光滑)，证明直达波数据零偏移距滤波提取地震子波是可行的，图5(d)即作为合成子波的相位谱；图5(f)是采用图5(a)和图5(d)反傅里叶变换计算的合成子波，与图5(e)形状十分相似。

为了进一步验证本申请的合成地震子波正确性，应用图5(f)合成子波对地震数据进行气泡压制。图5(g)和图5(h)是对图5(f)预测反褶积气泡压制后的期望输出和滤波算子，图6展示了将滤波算子应用到地震数据处理效果，剪头标识的位置显示气泡效应得到有效压制，气泡压制之后微幅构造更加清晰，偏移剖面更加透亮。

本发明公开了一种地震子波的确定方法，该方法包括：根据地震数据的振幅谱，确定待确定地震子波的振幅谱；确定地震直达波数据对应的震源位置，并根据震源位置确定滤波因子，滤波因子用于表示当前地震道的地震直达波数据与零偏移距地震道中的地震直达波数据之间的差异；根据滤波因子对地震直达波数据进行偏移距滤波处理，得到零偏移距直达波；根据零偏移距直达波的相位谱，确定待确定地震子波的相位谱；根据待确定地震子波的振幅谱和相位谱，确定目标地震子波。本发明实施例的技术方案，通过对非零偏移距直达波进行偏移距滤波，获得零偏移距直达波，其相位谱近似等于合成子波的相位谱，再结合子波振幅谱，实现快速、高效、准确的合成地震子波。

实施例三

图7为本发明实施例三提供的一种地震子波的确定装置的结构示意图。如图7所示，该装置包括：

振幅谱确定模块710，用于根据地震数据的振幅谱，确定待确定地震子波的振幅谱；

零偏移直达波确定模块720，用于对地震直达波数据进行偏移距滤波处理，得到零偏移距直达波；

相位谱确定模块730，用于根据所述零偏移距直达波的相位谱，确定待确定地震子波的相位谱；

地震子波确定模块740，用于根据待确定地震子波的振幅谱和相位谱，确定目标地震子波。

可选的，振幅谱确定模块710包括：

地震数据振幅谱确定单元，用于在确定反射系数为白噪声时，对地震数据进行傅里叶变换，得到地震数据的振幅谱；

地震子波振幅谱确定单元，用于将地震数据的振幅谱，作为待确定地震子波的振幅谱。

可选的，零偏移直达波确定模块720包括：

滤波因子确定单元，用于确定地震直达波数据对应的震源位置，并根据所述震源位置确定滤波因子，其中，所述滤波因子用于表示当前地震道的地震直达波数据与零偏移距地震道中的地震直达波数据之间的差异；

偏移距滤波处理单元，用于根据所述滤波因子对地震直达波数据进行偏移距滤波处理，得到零偏移距直达波。

可选的，滤波因子确定单元包括：

根据当前地震道的地震直达波数据与震源鬼波的到达时差，模拟得到第一脉冲响应；

根据零偏移距地震道中的地震直达波数据与震源鬼波的到达时差，模拟得到第二脉冲响应；

将第二脉冲响应与第一脉冲响应的比值，作为滤波因子。

可选的，相位谱确定模块730包括：

鬼波消除的零偏移直达波相位谱确定单元，用于若确定所述地震数据已进行检波点鬼波消除处理，则对所述零偏移距直达波进行傅里叶变换，得到所述零偏移距直达波的相位谱；

地震子波相位谱确定单元，用于将零偏移距直达波的相位谱，作为待确定地震子波的相位谱。

可选的，相位谱确定模块730还包括：

鬼波未消除的合成子波确定单元，用于若确定所述地震数据未进行检波点鬼波消除处理，则确定检波点鬼波，并将零偏移距直达波和检波点鬼波进行叠加，获得合成子波；

合成子波相位谱确定单元，用于对合成子波进行傅里叶变换，得到合成子波的相位谱；

地震子波相位谱确定单元，用于将合成子波的相位谱，作为待确定地震子波的相位谱。

可选的，地震子波确定模块740包括：

复数谱确定单元，用于根据待确定地震子波的振幅谱和相位谱，确定目标地震子波的复数谱；

地震子波确定单元，用于对复数谱进行反傅里叶变换，得到目标地震子波。

本发明实施例所提供的地震子波的确定装置可执行本发明任意实施例所提供的地震子波的确定方法，具备执行方法相应的功能模块和有益效果，详细过程参见前述实施例中地震子波的确定方法的相关操作。

实施例四

图8示出了可以用来实施本发明的实施例的电子设备10的结构示意图。电子设备旨在表示各种形式的数字计算机，诸如，膝上型计算机、台式计算机、工作台、个人数字助理、服务器、刀片式服务器、大型计算机、和其它适合的计算机。电子设备还可以表示各种形式的移动装置，诸如，个人数字处理、蜂窝电话、智能电话、可穿戴设备(如头盔、眼镜、手表等)和其它类似的计算装置。本文所示的部件、它们的连接和关系、以及它们的功能仅仅作为示例，并且不意在限制本文中描述的和/或者要求的本发明的实现。

如图8所示，电子设备10包括至少一个处理器11，以及与至少一个处理器11通信连接的存储器，如只读存储器(ROM)12、随机访问存储器(RAM)13等，其中，存储器存储有可被至少一个处理器执行的计算机程序，处理器11可以根据存储在只读存储器(ROM)12中的计算机程序或者从存储单元18加载到随机访问存储器(RAM)13中的计算机程序，来执行各种适当的动作和处理。在RAM 13中，还可存储电子设备10操作所需的各种程序和数据。处理器11、ROM 12以及RAM 13通过总线14彼此相连。输入/输出(I/O)接口15也连接至总线14。

电子设备10中的多个部件连接至I/O接口15，包括：输入单元16，例如键盘、鼠标等；输出单元17，例如各种类型的显示器、扬声器等；存储单元18，例如磁盘、光盘等；以及通信单元19，例如网卡、调制解调器、无线通信收发机等。通信单元19允许电子设备10通过诸如因特网的计算机网络和/或各种电信网络与其他设备交换信息/数据。

处理器11可以是各种具有处理和计算能力的通用和/或专用处理组件。处理器11的一些示例包括但不限于中央处理单元(中央处理器)、图形处理单元(GPU)、各种专用的人工智能(AI)计算芯片、各种运行机器学习模型算法的处理器、数字信号处理器(DSP)、以及任何适当的处理器、控制器、微控制器等。处理器11执行上文所描述的各个方法和处理，例如地震子波的确定方法。

在一些实施例中，地震子波的确定方法可被实现为计算机程序，其被有形地包含于计算机可读存储介质，例如存储单元18。在一些实施例中，计算机程序的部分或者全部可以经由ROM 12和/或通信单元19而被载入和/或安装到电子设备10上。当计算机程序加载到RAM 13并由处理器11执行时，可以执行上文描述的地震子波的确定方法的一个或多个步骤。备选地，在其他实施例中，处理器11可以通过其他任何适当的方式(例如，借助于固件)而被配置为执行地震子波的确定方法。

本文中以上描述的系统和技术的各种实施方式可以在数字电子电路系统、集成电路系统、场可编程门阵列(FPGA)、专用集成电路(ASIC)、专用标准产品(ASSP)、芯片上系统的系统(SOC)、负载可编程逻辑设备(CPLD)、计算机硬件、固件、软件、和/或它们的组合中实现。这些各种实施方式可以包括：实施在一个或者多个计算机程序中，该一个或者多个计算机程序可在包括至少一个可编程处理器的可编程系统上执行和/或解释，该可编程处理器可以是专用或者通用可编程处理器，可以从存储系统、至少一个输入装置、和至少一个输出装置接收数据和指令，并且将数据和指令传输至该存储系统、该至少一个输入装置、和该至少一个输出装置。

用于实施本发明的方法的计算机程序可以采用一个或多个编程语言的任何组合来编写。这些计算机程序可以提供给通用计算机、专用计算机或其他可编程数据处理装置的处理器，使得计算机程序当由处理器执行时使流程图和/或框图中所规定的功能/操作被实施。计算机程序可以完全在机器上执行、部分地在机器上执行，作为独立软件包部分地在机器上执行且部分地在远程机器上执行或完全在远程机器或服务器上执行。

在本发明的上下文中，计算机可读存储介质可以是有形的介质，其可以包含或存储以供指令执行系统、装置或设备使用或与指令执行系统、装置或设备结合地使用的计算机程序。计算机可读存储介质可以包括但不限于电子的、磁性的、光学的、电磁的、红外的、或半导体系统、装置或设备，或者上述内容的任何合适组合。备选地，计算机可读存储介质可以是机器可读信号介质。机器可读存储介质的更具体示例会包括基于一个或多个线的电气连接、便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦除可编程只读存储器(EPROM或快闪存储器)、光纤、便捷式紧凑盘只读存储器(CD-ROM)、光学储存设备、磁储存设备、或上述内容的任何合适组合。

为了提供与用户的交互，可以在电子设备上实施此处描述的系统和技术，该电子设备具有：用于向用户显示信息的显示装置(例如，CRT(阴极射线管)或者LCD(液晶显示器)监视器)；以及键盘和指向装置(例如，鼠标或者轨迹球)，用户可以通过该键盘和该指向装置来将输入提供给电子设备。其它种类的装置还可以用于提供与用户的交互；例如，提供给用户的反馈可以是任何形式的传感反馈(例如，视觉反馈、听觉反馈、或者触觉反馈)；并且可以用任何形式(包括声输入、语音输入或者、触觉输入)来接收来自用户的输入。

可以将此处描述的系统和技术实施在包括后台部件的计算系统(例如，作为数据服务器)、或者包括中间件部件的计算系统(例如，应用服务器)、或者包括前端部件的计算系统(例如，具有图形用户界面或者网络浏览器的用户计算机，用户可以通过该图形用户界面或者该网络浏览器来与此处描述的系统和技术的实施方式交互)、或者包括这种后台部件、中间件部件、或者前端部件的任何组合的计算系统中。可以通过任何形式或者介质的数字数据通信(例如，通信网络)来将系统的部件相互连接。通信网络的示例包括：局域网(LAN)、广域网(WAN)、区块链网络和互联网。

计算系统可以包括客户端和服务器。客户端和服务器一般远离彼此并且通常通过通信网络进行交互。通过在相应的计算机上运行并且彼此具有客户端-服务器关系的计算机程序来产生客户端和服务器的关系。服务器可以是云服务器，又称为云计算服务器或云主机，是云计算服务体系中的一项主机产品，以解决了传统物理主机与VPS服务中，存在的管理难度大，业务扩展性弱的缺陷。

应该理解，可以使用上面所示的各种形式的流程，重新排序、增加或删除步骤。例如，本发明中记载的各步骤可以并行地执行也可以顺序地执行也可以不同的次序执行，只要能够实现本发明的技术方案所期望的结果，本文在此不进行限制。

上述具体实施方式，并不构成对本发明保护范围的限制。本领域技术人员应该明白的是，根据设计要求和其他因素，可以进行各种修改、组合、子组合和替代。任何在本发明的精神和原则之内所作的修改、等同替换和改进等，均应包含在本发明保护范围之内。