沉积旋回划分的方法和装置与流程

本申请涉及油气勘探技术领域，特别涉及一种沉积旋回划分的方法和装置。

背景技术：

在进行油气勘探的过程中，在对目标地层进行沉积旋回划分时，往往需要先对沉积旋回进行分析。其中，上述的沉积旋回指的是自然界周期性海平面升降的结果。具体地：在时间上具有周期性，在成因上具有相似性，在地质上表现出旋回性，在地球物理上表现出一定响应特征。例如，地震波场中地震波频谱特征的规律性变化。因此，通过对沉积旋回的分析，可以进一步实现对于地震相及沉积相的分析、层序地层研究乃至对储层的预测。

目前，为了划分沉积旋回曲线，常用的方法是根据测井资料获得沉积旋回曲线。或者在缺少测井资料的地方，根据地震资料获取沉积旋回曲线。其中，所述根据地震资料获取沉积旋回曲线的方法，通常是根据计算出曲线的属性，根据曲线的趋势，人为主观地进行折线化，进而得到沉积旋回曲线。应用这种沉积旋回曲线划分方法，具体实施时，往往会存在误差大、效率低，只能定性或半定量分析的技术问题。

针对上述问题，目前尚未提出有效的解决方案。

技术实现要素：

本申请实施例提供了一种沉积旋回划分的方法和装置，以解决现有沉积旋回划分方法中划分沉积旋回曲线时存在的误差大、效率低的技术问题。

本申请实施例提供了一种沉积旋回划分的方法，包括：

对地震资料进行时频分析，得到时频数据体；

根据所述时频数据体，得到峰值频率属性曲线；

对所述峰值频率属性曲线进行分频滤波，得到低频属性的带通频率分量曲线和高频属性的带通频率分量曲线；

对所述低频属性的带通频率分量曲线进行折线化处理，得到第一沉积旋回曲线，对所述高频属性的带通频率分量曲线进行折线化处理，得到第二沉积旋回曲线；

根据所述第一沉积旋回曲线和所述第二沉积旋回曲线，进行沉积旋回划分。

在一个实施方式中，对所述峰值频率属性曲线进行分频滤波，得到低频属性的带通频率分量曲线和高频属性的带通频率分量曲线，包括：

根据待研究的沉积旋回级别，确定中心频率；

根据所述中心频率和通带带宽，设计带通滤波器；

通过所述带通滤波器对所述峰值频率属性曲线进行分频滤波，得到低频属性的带通频率分量曲线和高频属性的带通频率分量曲线。

在一个实施方式中，根据所述中心频率和通带带宽，设计带通滤波器，包括：

以所述中心频率与所述通带带宽的一半之间的差值作为低通频率，以所述中心频率与所述通带带宽的一半的相加值作为高通频率；

根据所述低通频率和所述高通频率，设计所述带通滤波器。

在一个实施方式中，对所述低频属性的带通频率分量曲线进行折线化处理，得到第一沉积旋回曲线，对所述高频属性的带通频率分量曲线进行折线化处理，得到第二沉积旋回曲线，包括：

根据所述低频属性的带通频率分量曲线，得到低频属性最值序列；

对所述低频属性最值序列进行数值处理，得到处理后的低频属性最值序列；

根据所述处理后的低频属性最值序列，通过线性插值得到所述第一沉积旋回曲线；

根据所述高频属性的带通频率分量曲线，得到高频属性最值序列；

对所述高频属性最值序列进行数值处理，得到处理后的高频属性最值序列；

根据所述处理后的高频属性最值序列，通过线性插值得到所述第二沉积旋回曲线。

在一个实施方式中，对所述低频属性最值序列进行数值处理，得到处理后的低频属性最值序列，包括：

根据所述待研究的沉积旋回级别，确定第一旋回阈值；

根据所述低频属性最值序列，确定低频属性最值序列的平均值；

按照低频属性最值序列的排列顺序，分别以所述低频属性最值序列中两个相邻的最值作为低频处理子单元，分别对所述低频属性最值序列中各个低频处理子单元的两个最值进行作差，得到对应的低频处理子单元的差值的绝对值；

根据所述低频处理子单元的差值的绝对值、所述低频属性最值序列的平均值和所述第一旋回阈值，对所述低频属性最值序列进行数值处理，得到处理后的低频属性最值序列。

在一个实施方式中，根据所述低频处理子单元的差值的绝对值、所述低频属性最值序列的平均值和所述第一旋回阈值，对所述低频属性最值序列进行数值处理，包括：

在所述低频处理子单元的差值的绝对值大于等于所述第一旋回阈值的情况下，则保留所述低频处理子单元中的两个最值和所述两个最值所分别对应的时间坐标；

在所述低频处理子单元的差值的绝对值小于所述第一旋回阈值，且所述低频处理子单元中的两个最值都小于等于所述低频属性最值序列的平均值的情况下，则保留所述低频处理子单元中的两个最值中较小的最值和所述较小的最值所对应的时间坐标；

在所述低频处理子单元的差值的绝对值小于所述第一旋回阈值，且所述低频处理子单元中的两个最值都大于所述低频属性最值序列的平均值的情况下，则保留所述低频处理子单元中的两个最值中较大的最值和所述较大的最值所对应的时间坐标；

在所述低频处理子单元的差值的绝对值小于所述第一旋回阈值，且所述低频处理子单元中的两个最值中一个较大最值大于所述低频属性最值序列的平均值，一个较小最值小于所述低频属性最值序列的平均值的情况下，则保留所述低频处理子单中的两个最值中与所述低频属性最值序列的平均值差值的绝对值较大的最值和对应的时间坐标。

在一个实施方式中，对所述高频属性最值序列进行数值处理，得到处理后的高频属性最值序列，包括：

根据所述待研究的沉积旋回级别，确定第二旋回阈值；

根据所述高频属性最值序列，确定高频属性最值序列的平均值；

按照高频属性最值序列的排列顺序，分别以所述高频属性最值序列中两个相邻的最值作为高频处理子单元，分别对所述高频属性最值序列中各个高频处理子单元的两个最值进行作差，得到对应的高频处理子单元的差值的绝对值；

根据所述高频处理子单元的差值的绝对值、所述高频属性最值序列的平均值和所述第二旋回阈值，对所述高频属性最值序列进行数值处理，得到处理后的高频属性最值序列。

在一个实施方式中，根据所述高频处理子单元的差值的绝对值、所述高频属性最值序列的平均值和所述第二旋回阈值，对所述高频属性最值序列进行数值处理，包括：

在所述高频处理子单元的差值的绝对值大于等于所述第二旋回阈值的情况下，则保留所述高频处理子单元中的两个最值和所述两个最值所对应的时间坐标；

在所述高频处理子单元的差值的绝对值小于所述第二旋回阈值，且所述高频处理子单元中的两个最值都小于等于所述高频属性最值序列的平均值的情况下，则保留所述高频处理子单元中的两个最值中较小的最值和所述较小的最值所对应的时间坐标；

在所述高频处理子单元的差值的绝对值小于所述第二旋回阈值，且所述高频处理子单元中的两个最值都大于所述高频属性最值序列的平均值的情况下，则保留所述高频处理子单元中的两个最值中较大的最值和所述较大的最值所对应的时间坐标；

在所述高频处理子单元的差值的绝对值小于所述第二旋回阈值，且所述高频处理子单元中的两个最值中一个较大最值大于所述高频属性最值序列的平均值，一个较小最值小于所述高频属性最值序列的平均值的情况下，则保留所述高频处理子单中的两个最值中与所述高频属性最值序列的平均值差值的绝对值较大的最值和对应的时间坐标。

在一个实施方式中，对所述低频属性的带通频率分量曲线进行折线化处理，得到第一沉积旋回曲线，对所述高频属性的带通频率分量曲线进行折线化处理，得到第二沉积旋回曲线之后，所述方法还包括：

对所述第一沉积旋回曲线进行多次折线化处理，直到处理后的第一旋回曲线符合预设要求，将符合预设要求的第一沉积旋回曲线作为所述第一沉积旋回曲线，

对所述第二沉积旋回曲线进行多次折线化处理，直到处理后的第二旋回曲线符合预设要求，将符合预设要求的第二沉积旋回曲线作为所述第二沉积旋回曲线。

基于相同的发明构思，本申请实施例还提供了一种沉积旋回划分的装置，包括：

第一获取模块，用于对地震资料进行时频分析，得到时频数据体；

第二获取模块，用于根据所述时频数据体，得到峰值频率属性曲线；

分频滤波模块，用于对所述峰值频率属性曲线进行分频滤波，得到低频属性的带通频率分量曲线和高频属性的带通频率分量曲线；

折线化处理模块，用于对所述低频属性的带通频率分量曲线进行折线化处理，得到第一沉积旋回曲线，对所述高频属性的带通频率分量曲线进行折线化处理，得到第二沉积旋回曲线；

沉积旋回划分模块，用于根据所述第一沉积旋回曲线和所述第二沉积旋回曲线，进行沉积旋回划分。

在一个实施方式中，所述分频滤波模块包括：

确定滤波器单元，用于根据待研究的沉积旋回级别，确定中心频率，并根据所述中心频率和通带带宽，确定带通滤波器；

分频滤波单元，用于通过所述带通滤波器对所述峰值频率属性曲线进行分频滤波，得到低频属性的带通频率分量曲线和高频属性的带通频率分量曲线。

在一个实施方式中，所述折线化处理模块包括：

低频最值序列获取单元，用于根据所述低频属性的带通频率分量曲线，得到低频属性最值序列；

低频最值序列处理单元，用于对所述低频属性最值序列进行数值处理，得到处理后的低频属性最值序列；

低频折线化处理单元，用于根据所述处理后的低频属性最值序列，通过线性插值得到所述第一沉积旋回曲线；

高频最值序列获取单元，用于根据所述高频属性的带通频率分量曲线，得到高频属性最值序列；

高频最值序列处理单元，用于对所述高频属性最值序列进行数值处理，得到处理后的高频属性最值序列；

高频折线化处理单元，用于根据所述处理后的高频属性最值序列，通过线性插值得到所述第二沉积旋回曲线。

在本申请实施例中，先通过滤波器进行分频处理，再对分频后的曲线分别进行折线化处理，得到相应的沉积旋回曲线。通过折线化处理的方式解决了现有沉积旋回划分方法中进行沉积旋回曲线划分时存在的误差大、效率低，只能定性或半定量分析的技术问题，达到了充分利用地震资料，快速、准确进行沉积旋回划分的技术效果。

附图说明

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

图1是根据本申请实施例的沉积旋回划分的方法的处理流程图；

图2是根据本申请实施例的沉积旋回划分的装置的组成结构图；

图3是应用本申请实施例提供沉积旋回划分的方法/装置的过程示意图；

图4是应用本申请实施例提供沉积旋回划分的方法/装置处理的试验区地震剖面图；

图5是应用本申请实施例提供沉积旋回划分的方法/装置获得的单道数据、单道时频分析和峰值频率属性示意图；

图6是应用本申请实施例提供沉积旋回划分的方法/装置获得的分频滤波后的时频属性图；

图7是应用本申请实施例提供沉积旋回划分的方法/装置获得的沉积旋回曲线。

具体实施方式

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

考虑到现有的沉积旋回划分方法，具体实施时，通常需要先划分相应的沉积旋回曲线。划分沉积旋回曲线一般都是根据计算出的频率属性曲线，根据曲线的趋势，人为主观地进行折线化。因此，导致现有沉积旋回曲线划分方法中划分沉积旋回曲线时往往存在误差大、效率低，只能定性或半定量分析的技术问题。针对产生上述技术问题的根本原因，本申请考虑可以先根据地震资料获得峰值频率属性曲线，再通过滤波器对该曲线进行分频滤波，对滤波后得到曲线分别进行相应的折线化处理，得到对应的沉积旋回曲线。从而解决了现有的沉积旋回曲线划分方法中划分沉积旋回曲线时存在的误差大、效率低，只能定性或半定量分析的技术问题，实现了充分利用地震资料，快速、准确进行沉积旋回划分的技术效果，达到了定量分析沉积旋回的目的。

基于上述思考思路，本申请提供了一种沉积旋回划分的方法。请参阅图1。本申请提供的沉积旋回划分的方法，具体可以包括以下步骤。

步骤101：对地震资料进行时频分析，得到时频数据体。

在一个实施方式中，为了根据地震资料获得时频数据，具体可以是：对叠后地震数据逐道进行广义S变换，得到时频数据体。需要说明的是，本实施方式中所提到的广义S变换只是为了更好的说明本申请实施例中的时频分析，所列举的一个示意说明。除通过广义S变换进行时频分析外，还可以通过短时傅里叶变化、小波变化和S变换等方法进行时频分析，获取时频数据体。具体实施时，可以根据具体情况，按照具体要求，灵活选择合适方法对地震数据进行时频分析，获取对应的时频数据体。对此，本申请不作限定。在本申请所涉及的实施方式中，主要通过广义S变换对地震数据进行时频分析。需要说明的是，在本实施方式中，从时间域地震数据体到时频域数据体可以是一个升维的过程。

步骤102：根据所述时频数据体，得到峰值频率属性曲线。

在本实施方式中，得到瞬时谱峰值频率属性数据体可以等同于得到所述峰值频率属性曲线。即当得到瞬时谱峰值频率属性数据体时，可以通过简单的方法，获得相应的峰值频率属性曲线。例如，可以根据瞬时谱峰值频率属性数据体绘制得到峰值频率属性曲线。也可以通过计算软件在计算得到瞬时谱峰值频率属性数据体的同时，自动生成相应的峰值频率属性曲线。对此，本申请不作限定。

在一个实施方式中，为了根据时频数据体获得瞬时谱峰值频率属性数据体，即峰值频率属性曲线。可以：先根据时频数据体获取瞬时谱峰值频率属性，其中，所述瞬时谱峰值频率属性可以是时间方向上单点振幅谱最大值所对应的频率。因此，可以通过对单点振幅谱序列进行比较求取瞬时谱峰值频率属性，进而根据求得的瞬时谱峰值频率属性得到瞬时谱峰值频率属性数据体，即峰值频率属性曲线。需要说明的是，在本实施方式中，相对的，从时频域数据体到瞬时谱峰值频率数据体是一个降维的过程。

步骤103：对所述峰值频率属性曲线进行分频滤波，得到低频属性的带通频率分量曲线和高频属性的带通频率分量曲线。

在一个实施方式中，为了得到低频属性的带通频率分量曲线和高频属性的带通频率分量曲线，具体实施时，可以按照以下步骤执行：

S1：根据待研究的沉积旋回级别，确定中心频率。

在本实施方式中，需要说明的是，所述旋回级别可以是根据使用者的需求而定的，有时候需要研究第一旋回，则设定相应的第一旋回级别；有时候需要研究第二旋回，则设定相应的第二旋回级别。对此，本申请不作限定。此外，还需要说明的是，中心频率可以依据旋回级别确定。例如，第一旋回选择的中心频率较低，第二旋回选择中心频率较高。在具体施工中，通常应用地震资料来研究旋回往往最多能达到三级。比如一级旋回选择的中心频率为10HZ，二级旋回选择的中心频率为15HZ，以此类推。具体如何根据旋回级别确定中心频率，可以根据具体实施情况和具体要求，灵活设定。对此，本申请不作限定。

S2：根据所述中心频率和通带带宽，设计带通滤波器。

在一个实施方式中，为了设计符合要求的带通滤波器进行分频滤波，具体可以：以所述中心频率与所述通带带宽的一半之间的差值作为低通频率，以所述中心频率与所述通带带宽的一半的相加值作为高通频率；根据所述低通频率和所述高通频率，设计所述带通滤波器。

在一个实施方式中，所述滤波器为巴特沃斯带通滤波器。当然还可以使用其他的滤波器。具体实施时，可以根据具体情况选择合适的滤波器。对此，本申请不作限定。

S3：通过所述带通滤波器对所述峰值频率属性曲线进行分频滤波，得到低频属性的带通频率分量曲线和高频属性的带通频率分量曲线。

在本实施方式中，需要说明的是，所述高频和所述低频在此都是相对的概念。并不特指某个具体的频率数值范围。结合本实施方式，可以举例说明：当通过滤波器分频滤波得到两条带通频率分量曲线时，可以以其中中心频率相对较高的一条带通频率分量曲线为高频属性的带通频率分量曲线；可以以其中中心频率相对较低的一条带通频率分量曲线为低频属性的带通频率分量曲线。对于具体的频率范围，本申请在此不作限定。

需要说明的是，在本申请实施例中，所述第一旋回具体指的是低级旋回，所述第二旋回具体指的是高级旋回。相应的，所述第一沉积旋回曲线具体指的是低级沉积旋回曲线，所述第二沉积旋回曲线具体指的是高级沉积旋回曲线。所述第一旋回级别具体指的是低级旋回级别，所述第二旋回级别具体指的是高级旋回级别。所述第一旋回阈值具体指的是低级旋回阈值，所述第二旋回阈值具体指的是高级旋回阈值。

步骤104：对所述低频属性的带通频率分量曲线进行折线化处理，得到第一沉积旋回曲线，对所述高频属性的带通频率分量曲线进行折线化处理，得到第二沉积旋回曲线。

在一个实施方式中，为了根据低频属性的带通频率分量曲线获得第一沉积旋回曲线，为了根据高频属性的带通频率分量曲线获得第二沉积旋回曲线，具体实施时可以按照以下步骤执行：

S1：根据所述低频属性的带通频率分量曲线，得到低频属性最值序列。

在一个实施方式中，为了获得低频属性最值序列，具体实施时，可以是：沿所述低频属性的带通频率分量曲线，每次从低频属性的带通频率分量曲线上第一预设个数的点中取出最大值和最小值；记录取出的最大值和最小值所对应的时间坐标；将所取出的多组最大值和最小值按照时间坐标表示的时间顺序排列，即可得到低频属性最值序列。在本实施方式中，需要说明的是第一预设个数的设置可以是用1除以中心频率得到的商乘以1000所得到的数值。例如，一条低频属性的带通频率分量曲线上按时间顺序排列有如下的点：1、3、4、6、7、9、1、4、7、10、5、6、3、9、8、5、3、6、2、1。假设计算得到的第一预设个数的数值为4，则每次取出4个点，从4个点中选出最大值和最小值。第一次，从1、3、4、6这4个点中选出最大值6和最小值1，并记录所选的6和1对应的时间坐标；第二次，从7、9、1、4这4个点中选出最大值9和最小值1，并记录这里的9和1对应的时间坐标；第三次，从7、10、5、6这4个点中选出最大值10和最小值5，并记录这里的10和5所对应的时间坐标；第四次，从3、9、8、5这4个点中选出最大值9和最小值3，并记录这里的9和3所对应的时间坐标；第五次，从3、6、2、1这4个点中选出最大值6和最小值1，并记录这里6和1所对应的时间坐标。根据选出的各个点的时间坐标，将选出的各个最大值点或最小值点按照各个点所对应的时间坐标表示的时间顺序进行排列，得到低频属性最值序列：1、6、9、1、10、5、3、9、6、1。

S2：对所述低频属性最值序列进行数值处理，得到处理后的低频属性最值序列。

在一个实施方式中，为了获得符合要求的处理后的低频属性最值序列，需要对低频属性最值序列进行数值处理，具体实施可以按照以下步骤执行：

S2-1：根据所述待研究的沉积旋回级别，确定第一旋回阈值。

在一个实施方式中，所述第一旋回阈值一般相对较大，例如，可以将第一旋回阈值设置为1。当然，这里设为1只是为了更好地说明本申请实施例。具体实施时，还可以设置为与1接近的其他数值，可以根据具体情况设置合适的数值作为第一旋回阈值。对此，本申请不作限定。

S2-2：根据所述低频属性最值序列，确定低频属性最值序列的平均值。

在一个实施方式中，例如低频属性最值序列为：1、6、9、1、10、5、3、9、6、1，对所述序列中各点求平均值为5.1，则5.1即为该低频属性最值序列的平均值。

S2-3：按照低频属性最值序列的排列顺序，分别以所述低频属性最值序列中两个相邻的最值作为低频处理子单元，分别对所述低频属性最值序列中各个低频处理子单元的两个最值进行作差，得到对应的低频处理子单元的差值的绝对值。

在一个实施方式中，为了对低频属性最值序列进行数值处理，需要先按照低频属性最值序列的排列顺序，依次取两个相邻的最值作为一个低频处理子单元，分别计算低频处理子单元的差值的绝对值。例如，对于低频属性最值序列1、6、9、1、10、5、3、9、6、1。可以按照该低频属性最值序列的排列顺序，依次取出1和6作为第一个低频处理子单元，取出9和1作为第二个低频处理子单元，取出10和5作为第三个低频处理子单元，取出3和9作为第四个低频处理子单元，取出6和1作为第五个低频处理子单元。分别计算各个低频处理子单元中两个数值的差值的绝对值，得到第一个低频处理子单元的差值的绝对值为5，第二个低频处理子单元的差值的绝对值为8，第三个低频处理子单元的差值的绝对值为5，第四个低频处理子单元的差值的绝对值为6，第五个低频处理子单元的差值的绝对值为5。

S2-4：根据所述低频处理子单元的差值的绝对值、所述低频属性最值序列的平均值和所述第一旋回阈值，对所述低频属性最值序列进行数值处理，得到处理后的低频属性最值序列。

在一个实施方式中，为了按要求根据所述低频处理子单元的差值的绝对值、所述低频属性最值序列的平均值和所述第一旋回阈值，对所述低频属性最值序列进行数值处理，具体可以按照以下方法执行：

在所述低频处理子单元的差值的绝对值大于等于所述第一旋回阈值的情况下，则保留所述低频处理子单元中的两个最值和所述两个最值所分别对应的时间坐标；

在所述低频处理子单元的差值的绝对值小于所述第一旋回阈值，且所述低频处理子单元中的两个最值都小于等于所述低频属性最值序列的平均值的情况下，则保留所述低频处理子单元中的两个最值中较小的最值和所述较小的最值所对应的时间坐标；

在所述低频处理子单元的差值的绝对值小于所述第一旋回阈值，且所述低频处理子单元中的两个最值都大于所述低频属性最值序列的平均值的情况下，则保留所述低频处理子单元中的两个最值中较大的最值和所述较大的最值所对应的时间坐标；

在所述低频处理子单元的差值的绝对值小于所述第一旋回阈值，且所述低频处理子单元中的两个最值中一个较大最值大于所述低频属性最值序列的平均值，一个较小最值小于所述低频属性最值序列的平均值的情况下，则保留所述低频处理子单中的两个最值中与所述高频属性最值序列的平均值差值的绝对值较大的最值和对应的时间坐标。

举例而言，对于低频属性最值序列1、6、9、1、10、5、3、9、6、1，经过计算得到第一个低频处理子单元的差值的绝对值为5，第二个低频处理子单元的差值的绝对值为8，第三个低频处理子单元的差值的绝对值为5，第四个低频处理子单元的差值的绝对值为6，第五个低频处理子单元的差值的绝对值为5。此外，该低频属性最值序列的平均值为5.1，为了方便进行示例说明，此处第一旋回阈值假设为6。由于第二、第四个低频处理子单元的差值的分别为8和6均大于等于第一旋回阈值，则保留第二、第四个低频处理子单元中的各个最值，即保留9、1和3、9，以及这4个最值对应的时间坐标。由于第一、第三、第五个低频处理子单元的差值的绝对值都为5，小于第一旋回阈值6，且上述各个低频处理子单元中的两个最值均符合其中一个较大最值大于所述低频属性最值序列的平均值，一个较小最值小于所述低频属性最值序列的平均值的情况。例如，第一个低频处理子单元中1小于5.1，而6大于5.1，同样第三个低频处理子单元中，10大于5.1，而5小于5.1，第五个低频处理子单元中，6大于5.1，而1小于5.1。因此，对上述各个低频处理子单元分别保留各个低频处理子单中的两个最值中与所述高频属性最值序列的平均值差值的绝对值较大的最值和对应的时间坐标。例如对于第一个低频处理子单元中最值1和5.1差值的绝对值为4.9，最值6和5.1的差值的绝对值为0.9，由于4.9大于0.9，因此第一个低频处理子单元保留最值1和对应的时间坐标。同理第三个低频处理子单元保留最值10和对应的时间坐标，第五个低频处理子单元保留最值1和对应的时间坐标。再将最后保留的最值按照所对应的时间坐标所表示的时间顺序排列，得到处理后的低频属性最值序列：1、9、1、10、3、9、1。

S3：根据所述处理后的低频属性最值序列，通过线性插值得到所述第一沉积旋回曲线。

在一个实施方式中，为了对处理后的低频属性最值序列进行线性插值，进而得到第一旋回曲线，具体可以按照以下步骤执行：

S3-1：按照处理后的低频属性最值序列的排列顺序，从所述处理后的低频属性最值序列中选取连续的三个最值作为一个低频折线处理组，得到多个低频折线处理组，其中，相邻的两个低频折线处理组共用一个最值；

S3-2：对所述多个低频折线处理组中的各个低频折线处理组分别按照以下步骤执行，得到各个处理后的低频折线组：

S3-2-1：确定所述低频折线处理组中三个最值中的最大值，并将所述最大值赋值为100，所述最大值的时间坐标保持不变；

S3-2-2：确定所述低频折线处理组中的三个最值中的最小值，并将所述最小值赋值为0，所述最小值的时间坐标保持不变；

S3-2-3：确定所述低频折线处理组中的三个最值中的中间值，并去掉所述中间值；

S3-3：根据所处理后的低频折线组，通过线性插值得到所述第一沉积旋回曲线。

举例而言，处理后的低频属性最值序列为：1、4、9、1、3、4、7、9、8。按照处理后的低频属性最值序列的排列顺序，依次从处理后的低频属性最值序列中取出3个最值得到一组低频折线处理组，取出1、4、9作为第一低频折线处理组，取出9、1、3作为第二低频折线处理组，取出3、4、7作为第三低频折线处理组，取出7、9、8作为第四低频折线处理组。需要说明的是，这里相邻的两个低频折线处理组共用一个最值。例如，第一低频折线处理组和第二低频折线处理组共用最值9，第二低频折线处理组与第三低频折线处理组共用最值3。参照规则，对上述四组低频折线处理组分别进行如下处理：将第一低频折线处理组中的数值最大值9赋值为100，但保留该9的时间坐标，即将原最值9改为100，而100的时间坐标还是原最值9的时间坐标。类似，将第一低频折线处理组中的数值最小值赋值为1，并保留原最值1的时间坐标作为该组数值为0的时间坐标。此外，在第一低频折线处理组中，最值4的数值大小既不是该组中的最大值，也不是该组中的最小值，是该组中的中间值，因此将最值4直接删除。得到处理后的第一低频折线处理组为0、100，其中0和100按照原来对应最值的时间坐标所表示的时间顺序排列。参照上述对第一低频折线处理组处理方法，可以得到处理后的第二低频折线处理组为100、0，处理后的第三低频折线处理组为0、100，处理后的第四低频折线处理组为0、100。进而可以根据处理后的第一低频折线组、处理后的第二低频折线组、处理后的第三低频折线处理组和处理后的第四低频折线处理组，对上述各个处理后的低频折线处理组中赋值后的值根据时间坐标位置进行线性插值，即可得到对应的第一沉积旋回曲线。

S4：根据所述高频属性的带通频率分量曲线，得到高频属性最值序列。

在一个实施方式中，为了获得高频属性最值序列，具体实施时，可以是：沿所述高频属性的带通频率分量曲线，每次从高频属性的带通频率分量曲线上第一预设个数的点中取出最大值和最小值；记录取出的最大值和最小值所对应的时间坐标；将所取出的多组最大值和最小值按照时间坐标表示的时间顺序排列，即可得到高频属性最值序列。在本实施方式中，需要说明的是第一预设个数的设置可以是用1除以中心频率得到的商乘以1000所得到的数值。具体实施可以参照低频属性最值序列的获取。在此，本申请不再赘述。

S5：对所述高频属性最值序列进行数值处理，得到处理后的高频属性最值序列。

在一个实施方式中，为了获得符合要求的处理后的高频属性最值序列，需要对高频属性最值序列进行数值处理，具体实施可以按照以下步骤执行：

S5-1：根据所述待研究的沉积旋回级别，确定第二旋回阈值。

在一个实施方式中，所述第二旋回阈值一般相对较小，例如，可以将第二旋回阈值设置为0.5。当然，这里设为0.5只是为了更好地说明本申请实施例。具体实施时，还可以设置为与0.5接近的其他数值，可以根据具体情况设置合适的数值作为第二旋回阈值。对此，本申请不作限定。

S5-2：根据所述高频属性最值序列，确定高频属性最值序列的平均值。

在本实施方式中，确定高频属性最值序列的平均值的方法和确定低频属性最值序列的平均值的方法类似，具体可以参照确定低频属性最值序列的平均值的方法。在此，本申请不再赘述。

S5-3：按照高频属性最值序列的排列顺序，分别以所述高频属性最值序列中两个相邻的最值作为高频处理子单元，分别对所述高频属性最值序列中各个高频处理子单元的两个最值进行作差，得到对应的高频处理子单元的差值的绝对值。

在一个实施方式中，为了对高频属性最值序列进行数值处理，需要先按照高频属性最值序列的排列顺序，依次取两个相邻的最值作为一个高频处理子单元，分别计算高频处理子单元的差值的绝对值。具体实施与计算低频处理子单元的差值的绝对值类似，可以参照计算低频处理子单元的差值的绝对值的过程。在此，本申请不再赘述。

S5-4：根据所述高频处理子单元的差值的绝对值、所述高频属性最值序列的平均值和所述第二旋回阈值，对所述高频属性最值序列进行数值处理，得到处理后的高频属性最值序列。

在一个实施方式中，为了按要求根据所述高频处理子单元的差值的绝对值、所述高频属性最值序列的平均值和所述第二旋回阈值，对所述高频属性最值序列进行数值处理，具体可以按照以下方法执行：

在所述高频处理子单元的差值的绝对值大于等于所述第二旋回阈值的情况下，则保留所述高频处理子单元中的两个最值和所述两个最值所对应的时间坐标；

在所述高频处理子单元的差值的绝对值小于所述第二旋回阈值，且所述高频处理子单元中的两个最值都小于等于所述高频属性最值序列的平均值的情况下，则保留所述高频处理子单元中的两个最值中较小的最值和所述较小的最值所对应的时间坐标；

在所述高频处理子单元的差值的绝对值小于所述第二旋回阈值，且所述高频处理子单元中的两个最值都大于所述高频属性最值序列的平均值的情况下，则保留所述高频处理子单元中的两个最值中较大的最值和所述较大的最值所对应的时间坐标；

在所述高频处理子单元的差值的绝对值小于所述第二旋回阈值，且所述高频处理子单元中的两个最值中一个较大最值大于所述高频属性最值序列的平均值，一个较小最值小于所述高频属性最值序列的平均值的情况下，则保留所述低频处理子单中的两个最值中与所述高频属性最值序列的平均值差值的绝对值较大的最值和对应的时间坐标。

具体实施时，对所述高频属性最值序列进行数值处理与对所述低频属性最值序列进行数值处理类似，可以参照对所述低频属性最值序列进行数值处理的过程。在此，本申请不再赘述。

S6：根据所述处理后的高频属性最值序列，通过线性插值得到所述第二沉积旋回曲线。

在一个实施方式中，为了对处理后的高频属性最值序列进行线性插值，进而得到第二旋回曲线，具体可以按照以下步骤执行：

S6-1：按照处理后的高频属性最值序列的排列顺序，从所述处理后的高频属性最值序列中选取连续的三个最值作为一个高频折线处理组，得到多个高频折线处理组，其中，相邻的两个高频折线处理组共用一个最值；

S6-2：对所述多个高频折线处理组中的各个高频折线处理组分别按照以下步骤执行，得到各个处理后的高频折线组：

S6-2-1：确定所述高频折线处理组中三个最值中的最大值，并将所述最大值赋值为100，所述最大值的时间坐标保持不变；

S6-2-2：确定所述高频折线处理组中的三个最值中的最小值，并将所述最小值赋值为0，所述最小值的时间坐标保持不变；

S6-2-3：确定所述高频折线处理组中的三个最值中的中间值，并去掉所述中间值；

S6-3：根据所处理后的高频折线组，通过线性插值得到所述第二沉积旋回曲线。

具体实施时，根据所述处理后的高频属性最值序列，通过线性插值得到所述第二沉积旋回曲线的具体处理过程与根据低频属性最值序列求得第一沉积旋回曲线类似，可以参照求解第一沉积旋回曲线的过程。在此，本申请不再赘述。

步骤105：根据所述第一沉积旋回曲线和所述第二沉积旋回曲线，进行沉积旋回划分。

在一个实施方式中，可以根据沉积旋回划分的结果，进行地质勘探。其中，所述根据沉积旋回划分的结果，进行地质勘探，可以包括：根据沉积旋回划分的结果，对地震相及沉积相进行分析，分析层序地层或者预测储层。具体地，可以根据沉积旋回划分的结果，进行井位部署。需要说明的是，上述所列举几种利用沉积旋回划分结果进行地质勘探的方式只是为了更好地说明本申请实施例，具体实施时，也可以根据具体情况或者施工要求，利用所述沉积旋回划分的结果进行其他相应的地质勘探工作。对此，本申请不作限定。

在本申请实施例中，由于先通过滤波器对峰值频率属性曲线进行分频处理，再对分频后的曲线分别进行折线化处理，得到相应的沉积旋回曲线，解决了现有沉积旋回划分方法中沉积旋回曲线划分时存在的误差大、效率低，只能定性或半定量分析的技术问题，实现了充分利用地震资料，快速、准确划分沉积旋回的技术效果，达到了定量分析沉积旋回的目的，进而可以为井位部署提供部署依据。

在一个实施方式中，为了获得更加准确的第一沉积旋回曲线和第二沉积旋回曲线，在根据所述处理后的低频属性最值序列，通过线性插值得到所述第一沉积旋回曲线，根据所述处理后的高频属性最值序列，通过线性插值得到所述第二沉积旋回曲线之后，所述方法还可以包括：重复线性差值处理，对所述第一沉积旋回曲线进行多次折线化处理，直到处理后的第一旋回曲线符合预设要求，将符合预设要求的第一沉积旋回曲线作为所述第一沉积旋回曲线，对所述第二沉积旋回曲线进行多次折线化处理，直到处理后的第二旋回曲线符合预设要求，将符合预设要求的第二沉积旋回曲线作为所述第二沉积旋回曲线。进而可以提高最后得到的第一沉积旋回曲线和第二沉积旋回曲线的准确度，进一步减少误差。

在一个实施方式，可以进一步扩展，可以根据具体要求或者施工情况，对3级旋回级别进行分析。由于是分析对象有3个旋回级别，对应的在设计滤波器进行分频时需要考虑到3个旋回级别，设计滤波器时，分别设定一级旋回中心频率和二级旋回中心频率，根据一级旋回中心频率和二级旋回中心频率设计所相应的分频滤波器，通过该滤波器得到3条带通滤波分量曲线，再对这3条带通滤波分量曲线分别进行相应的数值处理和折线化处理，具体可以参照上述根据低频带通滤波分量曲线获得第一沉积旋回曲线的过程进行处理，进而可以得到对应的3条沉积旋回曲线，实现对3级旋回级别的分析处理。当然对于1级旋回级别的分析处理，得到1条沉积旋回曲线，可以参照上述处理2级旋回级别，得到2条沉积旋回曲线和处理3级旋回级别，得到3条沉积旋回曲线的情况。在此，本申请不再赘述。

基于同一发明构思，本发明实施例中还提供了一种沉积旋回划分的装置，如下面的实施例所述。由于装置解决问题的原理与沉积旋回划分的方法相似，因此沉积旋回划分的装置的实施可以参见沉积旋回划分的方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。请参阅图2，是本发明实施例的沉积旋回划分的装置的一种组成结构图，该装置可以包括：第一获取模块201、第二获取模块202、分频滤波模块203、折线化处理模块204和沉积旋回划分模块205，下面对该结构进行具体说明。

第一获取模块201，用于对地震资料进行时频分析，得到时频数据体；

第二获取模块202，用于根据所述时频数据体，得到峰值频率属性曲线；

分频滤波模块203，用于对所述峰值频率属性曲线进行分频滤波，得到低频属性的带通频率分量曲线和高频属性的带通频率分量曲线；

折线化处理模块204，用于对所述低频属性的带通频率分量曲线进行折线化处理，得到第一沉积旋回曲线，对所述高频属性的带通频率分量曲线进行折线化处理，得到第二沉积旋回曲线；

沉积旋回划分模块205，用于根据所述第一沉积旋回曲线和所述第二沉积旋回曲线，进行沉积旋回划分。

在一个实施方式中，为了按照旋回级别，进行分频滤波，得到对应的多条带通滤波分量曲线，所述分频滤波模块204可以包括：

确定滤波器单元，用于根据待研究的沉积旋回级别，确定中心频率，并根据所述中心频率和通带带宽，确定带通滤波器；

分频滤波单元，用于通过所述带通滤波器对所述峰值频率属性曲线进行分频滤波，得到低频属性的带通频率分量曲线和高频属性的带通频率分量曲线。

在一个实施方式中，为了对带通频率分量曲线分别进行相应的数值处理和折线化处理，得到对应的沉积旋回曲线，所述折线化处理模块205可以包括：

低频最值序列获取单元，用于根据所述低频属性的带通频率分量曲线，得到低频属性最值序列；

低频最值序列处理单元，用于对所述低频属性最值序列进行数值处理，得到处理后的低频属性最值序列；

低频折线化处理单元，用于根据所述处理后的低频属性最值序列，通过线性插值得到所述第一沉积旋回曲线；

高频最值序列获取单元，用于根据所述高频属性的带通频率分量曲线，得到高频属性最值序列；

高频最值序列处理单元，用于对所述高频属性最值序列进行数值处理，得到处理后的高频属性最值序列；

高频折线化处理单元，用于根据所述处理后的高频属性最值序列，通过线性插值得到所述第二沉积旋回曲线。

在一个实施方式中，为了获得更加准确的第一沉积旋回曲线和第二沉积旋回曲线，在通过折线化处理模块205获得所述第一沉积旋回曲线和所述第二沉积旋回曲线之后，还可以通过折线化处理模块205对所述第一沉积旋回曲线和所述第二沉积旋回曲线分别进行多次的折线化处理，直到处理后的第二旋回曲线和第一沉积旋回曲线符合预设要求，将符合预设要求的第二沉积旋回曲线和第一沉积旋回曲线分别作为所述第二沉积旋回曲线和第一沉积旋回曲线。

在一个实施方式，可以进一步扩展，可以根据具体要求或者施工情况，对3级旋回级别进行分析。由于是分析对象有3个旋回级别，分频滤波模块204中的确定滤波器单元在执行具体操作时，考虑到3个旋回级别，确定相对应的滤波器：先分别设定一级旋回中心频率和二级旋回中心频率，再根据一级旋回中心频率和二级旋回中心频率设计所相应的分频滤波器，进而分频滤波模块204中的分频滤波单元可以通过该滤波器得到3条对应的带通滤波分量曲线，再通过折线化处理模块205这3条带通滤波分量曲线分别进行相应的数值处理和折线化处理，具体可以参照上述说明。在此，本申请不再赘述。当然对于1级旋回级别的分析处理，得到1条沉积旋回曲线，可以参照上述处理2级旋回级别，得到2条沉积旋回曲线和处理3级旋回级别，得到3条沉积旋回曲线的情况。在此，本申请同样也不再赘述。

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

需要说明的是，上述实施方式阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。为了描述的方便，在本说明书中，描述以上装置时以功能分为各种单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。

此外，在本说明书中，诸如第一和第二这样的形容词仅可以用于将一个元素或动作与另一元素或动作进行区分，而不必要求或暗示任何实际的这种关系或顺序。在环境允许的情况下，参照元素或部件或步骤(等)不应解释为局限于仅元素、部件、或步骤中的一个，而可以是元素、部件、或步骤中的一个或多个等。

从以上的描述中，可以看出，本申请实施例提供的沉积旋回划分的方法和装置由于先通过滤波器对峰值频率属性曲线进行分频处理，再对分频后的曲线分别进行折线化处理，得到相应的沉积旋回曲线，因此解决了现有方法中沉积旋回曲线划分时存在的误差大、效率低，只能定性或半定量分析的技术问题；由于本方法和装置是根据地震资料获取峰值频率属性曲线，再根据峰值频率属性曲线进行后续的处理，因此充分地利用了地震资料，解决了在测井资料缺少的地区无法根据测井资料沉积旋回曲线划分的技术问题，达到了扩展适用范围的技术效果，进而可以实现沉积旋回划分，为井位部署提供依据；又由于本方法和装置通过对分频过滤的曲线分别进行相对应的数值处理，再进行相应的折线化处理，提高了沉积旋回曲线划分的准确性，进而实现了快速、准确确定沉积相旋回曲线，进而进行沉积旋回划分的技术效果。

在一个具体的实施例中，应用本申请提供的沉积旋回划分的方法/装置对某地进行具体的沉积旋回的研究。具体可以按照以下方法执行。其中，具体的施工流程可以参阅图3。

请先参阅图4，待处理的是某地实际时间域纯波地震资料，其中，采样率为4ms。

步骤1：对时间域的地震数据体进行广义S变换时频分析，得到基于广义S变换时频分析的时频域数据体，即时频数据体。其中原始地震剖面可以参阅图4所示。选择过该地震剖面中某一点的单道地震数据参阅图5中的左图所示，对该单道数据进行广义S变换时频分析结果可以参阅图5中的中间图所示。

步骤2：对图5中的中间图所示的单道时频析结果，求取时间方向上单点振幅谱最大值所对应的频率，即瞬时谱峰值频率属性。具体结果可以参阅图5中的右图所示。

步骤3：根据所需要研究的沉积旋回级别(Ta～Tb),确定所要求取的沉积旋回曲线的中心频率(F＝1/T)。由此设计带通滤波器，使得低通频率和高通频率分别为：fa＝F-M和fb＝F+M(M一般取5～10)；对峰值频率属性曲线进行带通滤波，得到相应频率成分的带通频率分量曲线：低频属性的带通频率分量曲线CF1和高频属性的带通频率分量曲线CF2。具体可以参阅图6。

步骤4：分别对低频属性的带通频率分量曲线和高频属性的带通频率分量曲线进行折线化，具体可以根据所设计的带通滤波器中心频率F，对带通滤波后的带通频率分量曲线CF上每间隔T个样点上分别求取最大值Amax和最小值Amin并记录其对应的时间坐标C_Amax及C_Amin。按照时间坐标顺序排列所求取的最值序列。求取最值序列的平均值M_mean。根据待研究的沉积旋回级别(Ta～Tb)确定阈值，一般地，第一旋回阈值设置较大如1.0左右，第二旋回阈值设置较小如0.5左右。对最值序列中相邻两个最值进行做差，保留插值大于等于阈值的最值序列。对于差值小于阈值的最值做如下处理：如果相邻两最值小于平均值M_mean，则保留较小的最值及对应的时间坐标；反之，如果相邻两最值大于平均值M_mean，则保留较大的最值及对应的时间坐标；此外，如果相邻两最值一个大于平均值M_mean，一个小于平均值M_mean，则保留与平均值M_mean相差较大的最值及对应的时间坐标。对经过上述处理后的最值序列，每选择连续3个作为一组，其中下一组要包含上一组的最后一个值，然后确定这三个值中的最大值和最小值并且将最大值赋值为100，最小值赋值为0。对赋值后的值根据坐标位置进行线性插值，便得到了所要求取的沉积旋回曲线：第一沉积旋回曲线和第二沉积旋回曲线。具体结果可以参阅图7所示。

通过应用本申请提出的沉积旋回划分的方法/装置根据试验区内的一种无井地区叠后地震资料进行沉积旋回曲线的确定。由于是针对无井地区，所以通过工区内待钻井点的测井解释结果来验证其在该井点周围沉积旋回划分。通过比较，验证了本申请提供的沉积旋回划分的方法和装置解决了现有方法中存在的沉积旋回曲线划分误差大、效率低，只能定性或半定量分析的技术问题，实现了充分利用地震资料，快速、准确确定沉积相旋回曲线的技术效果，达到了定量分析沉积旋回的目的。

尽管本申请内容中提到不同的沉积旋回划分的方法或装置，但是，本申请并不局限于必须是行业标准或实施例所描述的情况等，某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取、处理、输出、判断方式等的实施例，仍然可以属于本申请的可选实施方案范围之内。

虽然本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。

上述实施例阐明的装置或模块等，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述模块的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个模块或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。

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

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

通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，移动终端，服务器，或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。

本说明书中的各个实施例采用递进的方式描述，各个实施例之间相同或相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。本申请可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。

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