阅读说明：本技术 基于saa监测的地裂缝活跃性分析方法及装置 (Ground fracture activity analysis method and device based on SAA monitoring ) 是由 刘祥磊 苏珊 马静 于 2019-09-10 设计创作，主要内容包括：本发明公开了一种基于SAA监测的地裂缝活跃性分析方法及装置,所述方法将由SAA监测取得地裂缝活动的原始时序位移信号分解成一系列简单子信号、再利用相关系数对简单子信号降噪处理,得到对应的有效子信号；对获取得到的有效子信号进行叠加得到重构的有效时序位移信号,并对有效时序位移信号进行分解、得到一系列对应的有效简单子信号；计算并获取每一个有效简单子信号的瞬时能量；对所有简单子信号的瞬时能量进行叠加,得到代表地裂缝活动整体强度的瞬时总能量。根据本发明的基于SAA监测的地裂缝活跃性分析方法实施例提高了地裂缝活动有效监测的精度,还提高了对其活跃性分析的准确性,从而降低地裂缝对城市建设和人居安全的危害性。(The invention discloses a method and a device for analyzing the activity of a ground fracture based on SAA monitoring, wherein the method decomposes an original time sequence displacement signal of the ground fracture activity obtained by the SAA monitoring into a series of simple sub-signals, and then carries out noise reduction processing on the simple sub-signals by utilizing a correlation coefficient to obtain corresponding effective sub-signals; superposing the obtained effective sub-signals to obtain reconstructed effective time sequence displacement signals, and decomposing the effective time sequence displacement signals to obtain a series of corresponding effective simple sub-signals; calculating and obtaining the instantaneous energy of each effective simple sub-signal; and superposing the instantaneous energy of all the simple sub-signals to obtain the instantaneous total energy representing the integral intensity of the ground fracture activity. According to the method for analyzing the activity of the ground fissure based on the SAA monitoring, the precision of effective monitoring of the activity of the ground fissure is improved, and the accuracy of analyzing the activity of the ground fissure is also improved, so that the harm of the ground fissure to urban construction and human settlement safety is reduced.)

基于SAA监测的地裂缝活跃性分析方法及装置

技术领域

本发明涉及地裂缝活跃性分析技术，尤其涉及一种基于SAA监测的地裂缝活跃性分析方法及装置。

背景技术

地裂缝是一种地质灾害，即地面在自然或者人为因素的影响下，产生裂缝。地裂缝具有隐蔽性强、危害性大的特点，可以直接破坏各种工程建设，例如线性工程、水利设施和城市建筑，与此同时严重制约城市规划、土地有效利用、地下水开发利用和地下空间发展。

由于地裂缝的规模、成因等不尽相同，现有技术中，对地裂缝的监测手段主要有三种，即水准测量、全球定位系统(GPS)监测、干涉合成孔径雷达测量(InSAR)技术。其中，水准测量是一种传统的地面变形监测技术，它可以在选定的位置提供高精度的地面沉降测量，具有操作简单、成本低的优点。但是，人工野外勘察测量的工作量大、周期长，无论是从时间还是空间上很难获得完整的地裂缝形变信息。其次，GPS监测在三维定位方面具有高精度的优势，可以获得监测区监测点的三维绝对形变量，然而，由于各种因素限制，布置的GPS监测点总是有限的，难以有效获取整个区域的地裂缝形变量。

InSAR技术具有高空间分辨率、全天时和全天候等优点，已经成功应用于许多地质灾害监测工程中。InSAR技术可以以年平均变形率图的形式提供监测到的地裂变形序列结果，并进一步检测到地裂活动变化。然而，由于InSAR技术在监测时缺乏绝对参考位置，并且缺乏消除大气等干扰误差的有效手段，所以其获取的地裂缝形变量精度不高。

由上述可见，现有技术中以上三种地裂缝活动监测技术受地形和土地范围的影响较大，导致监测结果不够可靠。因此，本发明基于上述问题提供一种地裂缝活跃性分析方法，以提高地裂缝活动有效监测的精度，并提高对其活跃性分析的准确性，从而降低地裂缝对城市建设和人居安全的危害性。

发明内容

本发明旨在至少在一定程度上解决相关技术中的技术问题之一。为此，本发明的第一个目的在于提出一种基于SAA监测的地裂缝活跃性分析方法。

本发明的第二个目的在于提出一种基于SAA监测的地裂缝活跃性分析装置。

本发明的第三个目的在于提出一种计算机设备。

本发明的第四个目的在于提出一种计算机存储介质。

为实现上述目的，第一方面，根据本发明实施例的基于SAA监测的地裂缝活跃性分析方法，所述方法包括：

获取由SAA监测取得的地裂缝上、下盘原始时序位移信号；

将所述原始时序位移信号分解成一系列简单子信号；

利用相关系数对每一个所述简单子信号进行降噪处理，得到一系列对应的有效子信号；

对获取得到的所有所述有效子信号进行叠加，得到重构的有效时序位移信号，并将所述有效时序位移信号分解成一系列对应的有效简单子信号；

根据所述有效简单子信号、计算并获取每一个所述有效简单子信号的瞬时能量；

对获取得到的所有所述简单子信号的瞬时能量进行叠加，得到代表地裂缝活动整体强度的瞬时总能量。

第二方面，根据本发明实施例的基于SAA监测的地裂缝活跃性分析装置，包括：

原始信号获取模块，用于获取由SAA监测取得的地裂缝上、下盘原始时序位移信号；

分解模块，用于将所述原始时序位移信号分解成一系列简单子信号；

降噪处理模块，用于利用相关系数对每一个所述简单子信号进行降噪处理，得到一系列对应的有效子信号；

重构及再分解模块，用于对获取得到的所有所述有效子信号进行叠加，得到重构的有效时序位移信号，并将所述有效时序位移信号分解成一系列对应的有效简单子信号；

瞬时能量计算模块，用于根据所述有效简单子信号、计算并获取每一个所述有效简单子信号的瞬时能量；

瞬时总能量计算模块，用于对获取得到的所有所述简单子信号的瞬时能量进行叠加，得到代表地裂缝活动整体强度的瞬时总能量。

第三方面，根据本发明实施例的计算机设备，包括存储器、处理器以及存储在所述存储器上并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如上所述的基于SAA监测的地裂缝活跃性分析方法。

第四方面，根据本发明实施例的计算机存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现如上所述的基于SAA监测的地裂缝活跃性分析方法。

根据本发明实施例提供的基于SAA监测的地裂缝活跃性分析方法及装置，通过将由SAA监测取得地裂缝活动的原始时序位移信号分解成一系列简单子信号、再利用相关系数对简单子信号降噪处理，得到对应的有效子信号；对获取得到的有效子信号进行叠加得到重构的有效时序位移信号，并对有效时序位移信号进行分解、得到一系列对应的有效简单子信号；计算并获取每一个有效简单子信号的瞬时能量；对所有简单子信号的瞬时能量进行叠加，得到代表地裂缝活动整体强度的瞬时总能量。根据本发明的基于SAA监测的地裂缝活跃性分析方法实施例提高了地裂缝活动有效监测的精度，还提高了对其活跃性分析的准确性，从而降低地裂缝对城市建设和人居安全的危害性。

本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。

附图说明

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

图1是本发明实施例基于SAA监测的地裂缝活跃性分析方法的流程图；

图2是本发明实施例基于SAA监测的地裂缝活跃性分析方法中步骤S102的流程图；

图3是本发明实施例基于SAA监测的地裂缝活跃性分析方法中步骤S103的流程图；

图4是本发明实施例基于SAA监测的地裂缝活跃性分析装置的结构示意图；

图5是本发明计算机设备一个实施例的结构示意图。

本发明目的的实现、功能特点及优点将结合实施例，参照附图做进一步说明。

具体实施方式

下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。

受地形和土地范围的影响，现有技术中的三种主要地裂缝活动监测技术：水准测量、全球定位系统(GPS)监测、干涉合成孔径雷达测量(InSAR)技术，监测结果都不够可靠。

一般而言，地裂缝活动会受多种外力影响，当不同方向的外力同时作用到地裂缝时，地裂缝可能不会产生位移活动，但此时地裂缝是处于活跃状态的。因此，直接使用时序位移数据进行地裂缝的活跃性分析是不可靠的。为进一步剖析地裂缝的活动规律，本发明提供了一种基于时序位移数据瞬时总能量的地裂缝活跃性分析方法。该瞬时能量活跃性分析方法主要包括：实测信号分解、相关系数降噪、信号重构再次分解、本征模函数(IMF)瞬时能量求解、叠加瞬时总能量。因此，本发明的目的在于提高地裂缝活动有效监测的精度，并提高对其活跃性分析的准确性，从而降低地裂缝对城市建设和人居安全的危害性。

参照图1所示，图1示出了本发明实施例提供的基于SAA监测的地裂缝活跃性分析方法一个实施例的流程图，为了便于描述，仅示出了与本发明实施例相关的部分。

本发明具体实施时，该基于SAA监测的地裂缝活跃性分析方法具体包括：

S101、获取由SAA监测取得的地裂缝上、下盘原始时序位移信号。

S102、将所述原始时序位移信号分解成一系列简单子信号。

S103、利用相关系数对每一个所述简单子信号进行降噪处理，得到一系列对应的有效子信号。

S104、对获取得到的所有所述有效子信号进行叠加，得到重构的有效时序位移信号，并将所述有效时序位移信号分解成一系列对应的有效简单子信号。

S105、根据所述有效简单子信号、计算并获取每一个所述有效简单子信号的瞬时能量。

S106、对获取得到的所有所述简单子信号的瞬时能量进行叠加，得到代表地裂缝活动整体强度的瞬时总能量。

具体地，本发明使用与以往地裂缝监测不同的方法，将阵列式位移传感器SAA(Shape Acceleration Array，简称SAA)可以埋设于地裂缝内部，获取地裂缝内部实时的形变信息，能够更加有效的分析地裂缝的活动规律及其影响因素，以降低地裂缝灾害带来的经济损失，保障人身财产安全。阵列式位移传感器SAA是基于微电子机械系统的阵列式位移传感器，内部由三轴微电子机械系统(MEMS)加速度计组成，可以连续测量被测目标的三维坐标，数据采样频率可以达到一小时获取一次，其精度可以达到±1.5mm/32m，本发明使用阵列式位移传感器SAA作为地裂缝监测的工具，将其布设于地裂缝内部可以获得精确、连续、高时间分辨率的地裂缝上盘、下盘三维时序位移信号，提高了地裂缝活动有效监测的精度，具有重要应用价值和意义。

进一步地，如图2所示，所述步骤S102具体包括：

S201、根据原始时序位移信号、找出所有局部的极值点。

S202、获取相邻所述极值点的中点，并采用插值法获取左右边界的中点。

S203、分别用奇数中点、偶数中点构造对应的两条插值曲线，并按时序取得所述两条插值曲线中的每一个平均值。

S204、根据设定的允许误差或最大筛选次数作为终止判断条件，按时序依次判断所述平均值的绝对值是否大于所述允许误差或当前筛选次数是否小于所述最大筛选次数。

S205、当判定出所述平均值对应的绝对值大于所述允许误差或当前筛选次数小于所述最大筛选次数，则产生对应的本征模函数。

S206、当判定出所述平均值对应的绝对值小于所述允许误差或当前筛选次数大于所述最大筛选次数，则将其对应的原始时序位移信号减去所述平均值得到的差值作为第一时序位移信号进行迭代，直至满足所述设定的允许误差或最大筛选次数的终止判断条件，从而产生对应的本征模函数。

S207、使用所述原始时序位移信号与产生的第一个本征模函数的差值作为第二时序位移信号进行迭代，直到最后一个残差项小于任何一个极值点值。迭代完成后，产生对应的一系列本征模函数作为原始时序位移信号分解得到的一系列简单子信号。

通过上述步骤S201至S207的算法，将由SAA获取得到的原始时序位移信号分解成一系列简单子信号。具体地，将原始时序位移信号分解成代表各段频率的简单信号时，简单信号在一定程度上反映原始时序位移信号的本质特征。

本发明具体实施时，具体算法步骤如下：

S1021、开始；

S1022、获取原始时序位移信号曲线S(t)，并获取S(t)曲线的所有局部极值点；

S1023、获取相邻所述极值点的中点，并采用插值法获取S(t)曲线左右边界的(n+1)中点；

S1024、分别用奇数中点、偶数中点构造对应的两条插值曲线，并按时序取得所述两条插值曲线中的每一个平均值L*；

S1025、设定作为终止判断条件的允许误差ε和最大筛选次数κ；

S1026、判断所述平均值L的绝对值|L*|是否大于允许误差ε或当前筛选次数是否小于设定的最大筛选次数κ；若是则转入步骤S1027、若否则转入步骤S1028；

S1027、产生对应的第j个本征模函数IMF；

S1028、计算差值S(t)＝S(t)-L*，将计算得到的S(t)代入步骤S1022中，进行迭代；

S1029、使用所述原始时序位移信号与产生的第j个本征模函数M的差值S(t)＝S(t)-M，代入步骤S1022中，进行迭代，直到最后一个残差项小于任何一个极值点值；

S1030、结束。

经上述算法迭代完成后，产生对应的m个本征模函数IMF作为原始时序位移信号分解得到的一系列简单子信号、以及一个趋势项。本发明将原始时序位移信号分解成代表各段频率的简单子信号时，简单子信号在一定程度上反映原始时序位移信号的本质特征。

进一步地，如图3所示，所述步骤S103具体还包括：

S301、计算得到每一个所述简单子信号与其对应原始时序位移信号之间的相关系数。

S302、根据所得到的相关系数，找出所述相关系数由逐渐减小到逐渐增大的临界点。

S303、将所述临界点之前简单子信号的分量作为噪声去除，得到一系列对应的有效子信号。

具体地，在利用SAA传感器进行地裂缝监测的过程中，由于交通、地面运动等多种环境因素的影响，不可避免地会产生噪声，从而导致所获得的原始时间序列位移信号的精度有所下降。因此，由复杂原始时间序列位移信号分解成的多个简单子信号中，一定会包含对应的有效信号与噪声信号。本发明采用相关性分析对得到的所有简单子信号进行降噪处理。相关性分析是噪声处理中的一个重要应用，本发明利用相关系数确定噪声占有的本征模函数IMF的分量，并将其去除达到降噪效果。

具体实施时，原始信号分解所得的本征模函数IMF与原始信号之间的相关系数在[-1,1]之间，噪声占主导地位的本征模函数IMF与原始信号之间的相关系数较小。每一个所述简单子信号与其对应原始时序位移信号之间的相关系数可以由以下公式计算得到：

其中，ρ指的是本征模函数IMF与原始信号之间的相关系数。

S(t)指的是原始时序位移信号曲线。

n指的是时刻；

IMF_j(t)指的是t时刻的本征模函数(1≦t≦n)。

进一步地，具体实施时，根据所得到的相关系数ρ，找出所述相关系数由逐渐减小到逐渐增大的临界点。将所述临界点之前本征模函数IMF的分量作为噪声去除，得到一系列对应的有效子信号。

进一步地，对获取得到的所有所述有效子信号进行叠加，得到重构的有效时序位移信号；并将所述有效时序位移信号按上述步骤S1021至1030的迭代算法，将所述有效时序位移信号分解成一系列对应的有效简单子信号。通过重构再分解生成的有效简单子信号既多尺度地还原了原始时序位移信号的本质特征，同时又能提高各个有效简单子信号随时间变化的瞬时能量的准确性，继而确保基于瞬时总能量分析地裂缝的活跃性的可靠性。

进一步地，在步骤S105中，根据所述有效简单子信号、计算并获取每一个所述有效简单子信号的瞬时能量。

具体地，以往的能量模型，将总能量视为定值，但事实上，由于分解出的每个有效简单子信号IMF分量的频率和瞬时能量是随时间变化的，因此，这种能量的求解方法是不合理的。在上述数据分解与降噪的基础上，本发明提出了一种基于SAA时序位移信号分解的瞬时能量模型，随着SAA信号数据的分解，可以把每个有效简单子信号本征模函数IMF看成是振幅随时间变化的曲线，因而每个本征模函数的瞬时能量可以以动能的形式表示，本发明由以下公式计算：

其中，E_j(t)指的是第j(1≦t≦m)个本征模函数IMF随时间t变化的瞬时能量。

G指的是监测的地裂缝质量。

Aj(t)指的是在t时刻的幅值(1≦t≦m)。

Δt指的是数据采集的时间间隔。

为了简化各个本征模函数IMF的瞬时能量，将G和Δt均设为1，故上述每个有效简单子信号本征模函数IMF的瞬时能量可以由以下公式计算：

其中，E_j(t)指的是第j(1≦t≦m)个本征模函数IMF随时间t变化的瞬时能量；Aj(t)指的是在t时刻的幅值(1≦t≦m)。

进一步地，由于地裂缝的位移活动受外力因素影响，但由于力的相互作用，地裂缝可能不会产生位移活动。因此，将SAA获取的原始时序位移信号分解成一系列简单子信号时，这些简单子信号能够多尺度还原原始时序位移信号的特征，可以将其看作是不同的外力激励作用。因此，可以将上述经降噪处理、并经重构再分解生成的每个有效简单子信号的瞬时能量进行叠加、形成瞬时总能量，将该瞬时总能量看作成地裂缝活动的整体强度。

本发明实施例中，瞬时总能量的计算公式如下：

E_j(t)指的是第j(1≦t≦m)个本征模函数IMF随时间t变化的瞬时能量。

SAA可以监测地裂缝不同位置的位移情况，通过瞬时能量的量化，能够判断地裂缝不同位置的相对活跃性。由于地裂缝处于正常活动时，瞬时总能量应该处于平稳状态，而当地裂缝处于形变活跃期或受外力严重影响时，瞬时总能量会出现上升趋势甚至是突变，因此叠加形成的瞬时总能量能指示地裂缝活动的规律、用于分析地裂缝的活跃性。

本发明实施获取由SAA监测取得地裂缝活动的原始时序位移信号、并将其分解成一系列简单子信号，再利用相关系数对简单子信号降噪处理、得到对应的有效子信号；对获取得到的有效子信号进行叠加得到重构的有效时序位移信号，并对有效时序位移信号进行分解、得到一系列对应的有效简单子信号；计算并获取每一个有效简单子信号的瞬时能量；对所有简单子信号的瞬时能量进行叠加，得到代表地裂缝活动整体强度的瞬时总能量。本发明为剖析地裂缝的活动规律，提供了一种基于时序瞬时总能量指示地裂缝活跃性的分析方法，以提高对地裂缝活跃性分析的准确性。

参照图4所示，图4示出了本发明实施例提供的基于SAA监测的地裂缝活跃性分析装置一个实施例的结构示意图，为了便于描述，仅示出了与本发明实施例相关的部分。具体的，该基于SAA监测的地裂缝活跃性分析装置10包括：

原始信号获取模块11，用于获取由SAA监测取得的地裂缝上、下盘原始时序位移信号。

分解模块12，用于将所述原始时序位移信号分解成一系列简单子信号。

降噪处理模块13，用于利用相关系数对每一个所述简单子信号进行降噪处理，得到一系列对应的有效子信号。

重构及再分解模块14，用于对获取得到的所有所述有效子信号进行叠加，得到重构的有效时序位移信号，并将所述有效时序位移信号分解成一系列对应的有效简单子信号。

瞬时能量计算模块15，用于根据所述有效简单子信号、计算并获取每一个所述有效简单子信号的瞬时能量。

瞬时总能量计算模块16，用于对获取得到的所有所述简单子信号的瞬时能量进行叠加，得到代表地裂缝活动整体强度的瞬时总能量。

进一步地，所述分解模块12具体还包括：

极值点获取单元，用于根据原始时序位移信号、找出所有局部的极值点。

中点获取单元，用于获取相邻所述极值点的中点，并采用插值法获取左右边界的中点。

插值曲线构造与平均值获取单元，用于分别用奇数中点、偶数中点构造对应的两条插值曲线，并按时序取得所述两条插值曲线中的每一个平均值。

判断单元，用于根据设定的允许误差或最大筛选次数作为终止判断条件，按时序依次判断所述平均值的绝对值是否大于所述允许误差或当前筛选次数是否小于所述最大筛选次数。

本征模函数生成单元，用于当判定出所述平均值对应的绝对值大于所述允许误差或当前筛选次数小于所述最大筛选次数，则产生对应的本征模函数。

第一迭代单元，用于当判定出所述平均值对应的绝对值小于所述允许误差或当前筛选次数大于所述最大筛选次数，则将其对应的原始时序位移信号减去所述平均值得到的差值作为第一时序位移信号进行迭代，直至满足所述设定的允许误差或最大筛选次数的终止判断条件。从而产生对应的本征模函数。

第二迭代单元，用于使用所述时序位移信号与产生的第一个本征模函数的差值作为第二时序位移信号进行迭代，直到最后一个残差项小于任何一个极值点值。迭代完成后，产生对应的一系列本征模函数作为原始时序位移信号分解得到的一系列简单子信号。

进一步地，所述降噪处理模块13具体还包括：

相关系数计算单元，用于计算得到每一个所述简单子信号与其对应原始时序位移信号之间的相关系数。

临界点查找单元，用于根据所得到的相关系数，找出所述相关系数由逐渐减小到逐渐增大的临界点。

噪声去除单元，用于将所述临界点之前简单子信号的分量作为噪声去除，得到一系列对应的有效子信号。

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

参照图5所示，图5示出了本发明实施例提供的计算机设备实施例的结构示意图，为了便于描述，仅示出了与本发明实施例相关的部分。具体的，该计算机设备500包括存储器502、处理器501以及存储在所述存储器502中并可在所述处理器501上运行的计算机程序5021，所述处理器501执行所述计算机程序时实现如上述实施例所述方法的步骤，例如图1所示的S101至S106的步骤。或者，所述处理器501执行所述计算机程序时实现上述实施例所述装置中的各模块/单元的功能，例如图4所示模块11至16的功能。

示例性的，所述计算机程序可以被分割成一个或多个模块/单元，所述一个或者多个模块/单元被存储在所述存储器502中，并由所述处理器501执行，以完成本发明。所述一个或多个模块/单元可以是能够完成特定功能的一系列计算机程序指令段，该指令段用于描述所述计算机程序在所述计算机设备500中的执行过程。例如，所述计算机程序可以被分割成原始信号获取模块11、分解模块12、降噪处理模块13、重构及再分解模块14、瞬时能量计算模块15、瞬时总能量计算模块16。

其中，原始信号获取模块11，用于获取由SAA监测取得的地裂缝上、下盘原始时序位移信号。

分解模块12，用于将所述原始时序位移信号分解成一系列简单子信号。

降噪处理模块13，用于利用相关系数对每一个所述简单子信号进行降噪处理，得到一系列对应的有效子信号。

重构及再分解模块14，用于对获取得到的所有所述有效子信号进行叠加，得到重构的有效时序位移信号，并将所述有效时序位移信号分解成一系列对应的有效简单子信号。

瞬时能量计算模块15，用于根据所述有效简单子信号、计算并获取每一个所述有效简单子信号的瞬时能量。

瞬时总能量计算模块16，用于对获取得到的所有所述简单子信号的瞬时能量进行叠加，得到代表地裂缝活动整体强度的瞬时总能量。

所述计算机设备500可包括，但不仅限于处理器501、存储器502。本领域技术人员可以理解，图仅仅是计算机设备500的示例，并不构成对计算机设备500的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件，例如所述计算机设备500还可以包括输入输出设备、网络接入设备、总线等。

所称处理器501可以是中央处理单元(Central Processing Unit，CPU)，还可以是其他通用处理器501、数字信号处理器501(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现场可编程门阵列(FieldProgrammable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立预设硬件组件等。通用处理器501可以是微处理器501或者该处理器501也可以是任何常规的处理器501等。

所述存储器502可以是所述计算机设备500的内部存储单元，例如计算机设备500的硬盘或内存。所述存储器502也可以是所述计算机设备500的外部存储设备，例如所述计算机设备500上配备的插接式硬盘，智能存储卡(SmartMedia Card,SMC)，安全数字(SecureDigital,SD)卡，闪存卡(Flash Card)等。进一步地，所述存储器502还可以既包括所述计算机设备500的内部存储单元也包括外部存储设备。所述存储器502用于存储所述计算机程序5021以及所述计算机设备500所需的其他程序和数据。所述存储器502还可以用于暂时地存储已经输出或者将要输出的数据。

本发明实施例还提供了一种计算机可读存储介质，计算机可读存储介质存储有计算机程序，计算机程序被处理器501执行时实现如上述实施例中所述方法中的步骤，例如图1所示的步骤S101至S106。或者，所述计算机程序被处理器501执行时实现上述实施例中所述装置中的各模块/单元的功能，例如图4所示的模块11至16的功能。

所述的计算机程序可存储于一计算机可读存储介质中，该计算机程序在被处理器501执行时，可实现上述各个方法实施例的步骤。其中，所述计算机程序包括计算机程序代码，所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括：能够携带所述计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、电载波信号、电信信号以及软件分发介质等。

需要说明的是，所述计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减，例如在某些司法管辖区，根据立法和专利实践，计算机可读介质不包括是电载波信号和电信信号。

在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述或记载的部分，可以参见其它实施例的相关描述。

本发明实施例方法中的步骤可以根据实际需要进行顺序调整、合并和删减。

本发明实施例系统中的模块或单元可以根据实际需要进行合并、划分和删减。

本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子预设硬件、或者计算机软件和电子预设硬件的结合来实现。这些功能究竟以预设硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。

本发明所提供的实施例中，应该理解到，所揭露的装置/计算机设备500和方法，可以通过其它的方式实现。例如，以上所描述的装置/计算机设备500实施例仅仅是示意性的，例如，所述模块或单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通讯连接可以是通过一些接口，装置或单元的间接耦合或通讯连接，可以是电性，机械或其它的形式。

以上所述实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围，均应包含在本发明的保护范围之内。