一种基于特征反射非线性褶积模型的多次波压制方法

1.本发明涉及地球物理勘探技术领域，尤其涉及一种基于特征反射非线性褶积模型的多次波压制方法。


背景技术：

2.海底节点采集技术因其长偏移距、全方位观测，以及丰富的低频信号等特点，有利于实现高精度的地震成像，在复杂构造成像方面具有重要作用，但海底节点观测系统独特的采集方式也给后续的数据处理带来了挑战。海底节点观测系统中震源在海面密集激发，观测节点位于海底，允许稀疏不规则布设，炮检点位置的明显深度差异，使得海底节点数据具有双基准面特点。由于海水-空气界面、海水-海底界面这两个强特征反射面的存在，海底节点数据中不仅包含来自地层的上行波信号，还包含来自海面的下行波信号，水体及自由表面相关的多次波也广泛发育。
3.由于目前的地震数据处理流程及成像方法主要针对一次反射，在进行速度分析以及成像等处理之前，通常需要先将上下行波信号分离，多次波的存在也会影响速度分析以及成像的质量，多次波压制是海底节点数据处理的重要议题。传统的自由表面相关多次波压制方法要求采集到的数据是规则完整的，但在海底节点观测系统中，海底观测节点通常是稀疏布设的，采集到的地震数据在共炮道集中空间采样不足，这就使得传统的自由表面相关多次波压制方法不能适应海底节点观测方式采集的地震数据的问题。


技术实现要素：

4.本部分的目的在于概述本发明的实施例的一些方面以及简要介绍一些较佳实施例。在本部分以及本技术的说明书摘要和发明名称中可能会做些简化或省略以避免使本部分、说明书摘要和发明名称的目的模糊，而这种简化或省略不能用于限制本发明的范围。
5.鉴于上述现有存在的问题，提出了本发明。
6.因此，本发明提供了一种基于特征反射非线性褶积模型的多次波压制方法，能够解决传统的自由表面相关多次波压制方法不能适应海底节点观测方式采集的地震数据的问题。
7.为解决上述技术问题，本发明提供如下技术方案，一种基于特征反射非线性褶积模型的多次波压制方法，包括：
8.收集多分量海底节点数据，并对所述数据进行预处理；
9.根据预处理后的参数构建平面波基函数，将频率空间域数据转换为时间平面波域数据；
10.引入自适应匹配因子，实现时间平面波域数据的自适应匹配，分离出上行波和下行波海底节点数据；
11.分别将分离出的上行波和下行波海底节点数据转化到频率域，进行上下行波反褶积；
12.将反褶积结果转换回时空域，得到一次反射格林函数，也即一次反射海底节点数据，实现压制。
13.作为本发明所述的基于特征反射非线性褶积模型的多次波压制方法的一种优选方案，其中：所述预处理包括，向道头关键字程序输入多分量海底节点数据，获取观测系统信息，分选出压力分量和垂直速度分量的共检波点道集数据。
14.作为本发明所述的基于特征反射非线性褶积模型的多次波压制方法的一种优选方案，其中：所述观测系统信息包括，时间采样点数、时间采样间隔、观测节点相对于观测系统的空间位置、海面震源相对于观测系统的激发位置和地震数据的倾角采样间隔。
15.作为本发明所述的基于特征反射非线性褶积模型的多次波压制方法的一种优选方案，其中：所述平面波基函数包括，
16.在平面波域中根据压力分量和垂直速度分量数据量级的物理关系，将二者的振幅水平调整为同一量级，量级的调整关系如下：
[0017][0018]
其中，p
′
表示量级调整后的压力分量，kz表示垂直方向的波数分量，ρ表示海水密度，ω表示傅里叶变换后的信号波数，p表示量级调整前的压力分量，p
x
表示局部平面波信号横向的射线参数分量，v
sea_water
表示海水的声波传播速度。
[0019]
作为本发明所述的基于特征反射非线性褶积模型的多次波压制方法的一种优选方案，其中：所述自适应匹配因子包括，选择直达波并由最小二乘法求得，方程构建如下：
[0020][0021]
其中，p表示压力分量数据，vz表示垂向速度分量数据，w表示压力分量和垂向速度分量的最佳匹配滤波器，表示垂向速度分量数据组成的褶积矩阵，表示矩阵的转置，d
p
表示压力分量数据组成的矩阵，和d
p
的矩阵排列方式参考一般的褶积矩阵方程的定义。
[0022]
作为本发明所述的基于特征反射非线性褶积模型的多次波压制方法的一种优选方案，其中：所述直达波包括，直达波在时空域的轨迹表示为：
[0023][0024][0025]
其中，t表示直达波信号在时空域的到达时曲线，即震源信号到达检波器前在海水中的传播时间，v
sea_water
表示海水的声波传播速度，z表示检波器放置的海水深度，x0表示检波器放置的横向空间坐标，x表示震源激发的横向空间坐标，y表示检波器放置的纵向空间坐标，y0表示震源激发的纵向空间坐标，表示震源距离检波器的横向空间距离，表示震源距离检波器的纵向空间距离。
[0026]
作为本发明所述的基于特征反射非线性褶积模型的多次波压制方法的一种优选
方案，其中：所述直达波还包括，时间平面域直达波轨迹表示如下：
[0027][0028]
其中，t0表示直达波在平面波域的到达时曲线，p
x
表示直达波的横向射线参数分量，py表示直达波的纵向射线参数分量，横、纵向射线参数的定义如下：
[0029][0030][0031]
其中，各个变量的物理含义和前述公式中的定义相同。
[0032]
作为本发明所述的基于特征反射非线性褶积模型的多次波压制方法的一种优选方案，其中：所述上行波和下行波海底节点数据包括，
[0033]
海底节点数据上下行波即其中所含自由表面相关多次波的形成过程，在平面波域中分别表示为：
[0034][0035][0036]
其中，p表示信号的射线参数向量(包含各个方向上射线参数的分量信息)，ω表示傅里叶变换后的信号波数，zr表示表示检波器放置的海水深度，z0表示震源激发的水深，exp()表示自然指数函数，u(p,ω,zr)表示频率平面波域的上行波信号，p(p,ω,z0)表示在海面处观测到的全波场信号，i表示复信号的虚数单位，kz表示波场的垂向波数分量，g(p,ω)表示一次反射格林函数(即需要提取的一次波信息)，a(ω)表示海面的波阻抗(通常取为常数-1)，s(p,ω)表示震源的激发信号，d(p,ω,zr)表示频率平面波域的下行波信号，pd(p,ω,zr)表示直达波信号。
[0037]
作为本发明所述的基于特征反射非线性褶积模型的多次波压制方法的一种优选方案，其中：所述反射格林函数包括，
[0038][0039]
其中，p
(u)
(p,ω,zr)表示频率平面波域的上行波信号(某个射线参数上的复数据点)，p
(d)
(p,ω,zr)表示表示频率平面波域的下行波信号(某个射线参数上的复数据点)，g
(p,ω,zr)表示最终求得的一次反射格林函数(即需要提取的一次波信息)。
[0040]
作为本发明所述的基于特征反射非线性褶积模型的多次波压制方法的一种优选方案，其中：所述第反褶积包括，
[0041]
为了保证反褶积过程的稳定性，引入稳定因子，将反褶积过程改进如下：
[0042][0043]
其中，表示p
(d)
(p,ω,zr)的复共轭，ε表示一个小的正数用于稳定求解。
[0044]
本发明的有益效果：本发明方法基于反馈环模型，得到了海底节点数据上下行波的非线性褶积模型，该模型表明，可以通过上下行波的反褶积来提取地下的一次反射格林函数，该方法假设海底节点数据中不同方向的平面波成分的多次波预测关系是相互独立的，因而可以在平面波域逐频率、逐平面波方向进行反褶积，计算效率高，并且只需要对海底节点数据的共检波点道集进行局部平面波分解，避免了海底节点数据共炮道集数据空间采样稀疏的问题。
附图说明
[0045]
为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其它的附图。其中：
[0046]
图1为本发明一个实施例提供的一种基于特征反射非线性褶积模型的多次波压制方法流程示意图；
[0047]
图2为本发明一个实施例提供的一种基于特征反射非线性褶积模型的多次波压制方法的反馈环模型示意图。
[0048]
图3为本发明一个实施例提供的一种基于特征反射非线性褶积模型的多次波压制方法的海底节点观测系统示意图；
[0049]
图4为本发明一个实施例提供的一种基于特征反射非线性褶积模型的多次波压制方法的应用数据示意图；
[0050]
图5为本发明一个实施例提供的一种基于特征反射非线性褶积模型的多次波压制方法的上下行波分离示意图；
[0051]
图6为本发明一个实施例提供的一种基于特征反射非线性褶积模型的多次波压制方法的上下行波反褶积提取一次反射格林函数示意图；
[0052]
图7为本发明一个实施例提供的一种基于特征反射非线性褶积模型的多次波压制方法的处理结果与原海底节点数据共检波点道集的对比分析图；
具体实施方式
[0053]
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合说明书附图对本发明的具体实施方式做详细的说明，显然所描述的实施例是本发明的一部分实施例，而不是全部实施例。基于本发明中的实施例，本领域普通人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明的保护的范围。
[0054]
在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是本发明还可以采用其他不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本发明内涵的情况下做类似推广，因此本发明不受下面公开的具体实施例的限制。
[0055]
其次，此处所称的“一个实施例”或“实施例”是指可包含于本发明至少一个实现方式中的特定特征、结构或特性。在本说明书中不同地方出现的“在一个实施例中”并非均指同一个实施例，也不是单独的或选择性的与其他实施例互相排斥的实施例。
[0056]
本发明结合示意图进行详细描述，在详述本发明实施例时，为便于说明，表示器件结构的剖面图会不依一般比例作局部放大，而且所述示意图只是示例，其在此不应限制本发明保护的范围。此外，在实际制作中应包含长度、宽度及深度的三维空间尺寸。
[0057]
同时在本发明的描述中，需要说明的是，术语中的“上、下、内和外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一、第二或第三”仅用于描述目的，而不能理解为指示或暗示相对重要性。
[0058]
本发明中除非另有明确的规定和限定，术语“安装、相连、连接”应做广义理解，例如：可以是固定连接、可拆卸连接或一体式连接；同样可以是机械连接、电连接或直接连接，也可以通过中间媒介间接相连，也可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
[0059]
实施例1
[0060]
参照图1-7，为本发明的第一个实施例，该实施例提供了一种基于特征反射非线性褶积模型的多次波压制方法，包括：
[0061]
s1：收集多分量海底节点数据，并对所述数据进行预处理；
[0062]
更进一步的，所述预处理包括向道头关键字程序输入标准seg-y格式的多分量海底节点数据，获取观测系统信息，根据seg-y格式的道头关键字信息分选出压力分量和垂直速度分量的共检波点道集数据。
[0063]
应说明的是，所述观测系统参数为：时间采样点数、时间采样间隔、观测节点相对于观测系统的空间位置、海面震源相对于观测系统的激发位置和地震数据的倾角采样间隔。
[0064]
还应说明的是，图2中可以看出震源在海面处密集激发，而检波器(用三角形表示)在海底稀疏不规则布设。
[0065]
s2：根据预处理后的参数构建平面波基函数，将频率空间域数据转换为时间平面波域数据；
[0066]
更进一步的，按数据道进行一维快速傅里叶变换，将时间(t)-空间(x)域的压力分量和垂直速度分量共检波点道集地震数据d(x,t)转换到频率-空间域d(x,ω)。
[0067]
应说明的是，根据观测系统参数构建平面波基函数，将频率-空间域数据d(x,ω)转换为平面波域数据s(p,t)，在平面波域中根据压力分量p(x,t)和垂直速度分量vz(x,t)数据量级的物理关系，将二者的振幅水平调整为同一量级，量级的调整关系如下：
[0068][0069]
其中，p
′
表示量级调整后的压力分量，kz表示垂直方向的波数分量，ρ表示海水密度，ω表示傅里叶变换后的信号波数，p表示量级调整前的压力分量，p
x
表示局部平面波信号横向的射线参数分量，v
sea_water
表示海水的声波传播速度。
[0070]
更进一步的，海底节点地震数据的倾角范围可以根据海水速度计算，结合地震数据的倾角采样间隔构建局部平面波分解所需要的平面波基函数，进行地震数据的局部平面波分解。
[0071]
应说明的是，对时间-空间域数据按道进行快速傅里叶变换时，使用多线程并行计算提高处理效率。
[0072]
s3：引入自适应匹配因子，实现时间平面波域数据的自适应匹配，分离出上行波和下行波海底节点数据；
[0073]
更进一步的，引入自适应匹配因子w，实现压力分量数据和垂直速度分量数据的自适应匹配，分离出上行波和下行波海底节点数据。
[0074]
应说明的是，其中上行波分量为u
′
＝(p w*vz)，下行波分量为d
′
＝(p-w*vz)。
[0075]
更进一步的，自适应匹配因子w，使用最小二乘方法求得，方程构建如下：
[0076][0077]
其中，p表示压力分量数据，vz表示垂向速度分量数据，w表示压力分量和垂向速度分量的最佳匹配滤波器，表示垂向速度分量数据组成的褶积矩阵，表示矩阵的转置，d
p
表示压力分量数据组成的矩阵，和d
p
的矩阵排列方式参考一般的褶积矩阵方程的定义。
[0078]
应说明的是，，应选择海底节点数据中的直达波信号构建方程求解自适应匹配因子，这是由于在海底节点数据中直达波属于下行波，求解得到的自适应匹配因子可以使得数据中的下行波能量极小化，从而分离海底节点数据上下行波。同时在海底节点数据中直达波信号的轨迹相对固定，可以用解析公式求得，之后根据地震信号子波的长度估计直达波信号数据窗的范围。
[0079]
更进一步的，海底节点数据中直达波在时空域的轨迹表示为：
[0080][0081][0082]
其中，t表示直达波信号在时空域的到达时曲线(也即震源信号到达检波器前在海水中的传播时间)，v
sea_water
表示海水的声波传播速度，z表示检波器放置的海水深度，x0表示检波器放置的横向空间坐标，x表示震源激发的横向空间坐标，y表示检波器放置的纵向
空间坐标，y0表示震源激发的纵向空间坐标，表示震源距离检波器的横向空间距离，表示震源距离检波器的纵向空间距离。
[0083]
应说明的是，由于本发明专利是在时间平面波域进行数据处理的，将上述直达波轨迹转换到时间平面波域表示为：
[0084][0085]
其中，t0表示直达波在平面波域的到达时曲线，p
x
表示直达波的横向射线参数分量，py表示直达波的纵向射线参数分量，横、纵向射线参数的定义如下：
[0086][0087][0088]
其中，各个变量的物理含义和前述公式中的定义相同。
[0089]
应说明的是，时间平面波域的直达波轨迹可以由时间域的直达波轨迹t计算得到。将海底节点数据分离为上下行波后，结合反馈环模型和海底节点观测系统的特点，可以得到海底节点数据中上下行波的非线性褶积模型。反馈环模型是描述地震数据自由表面相关多次波的经典模型，其将来自地下的上行反射视为“二次震源”再次输入地下；
[0090]
应说明的是，上述反馈环模型通常假设地震数据观测系统中震源和信号接收器都位于自由表面处，但在海底节点观测系统中信号接收器布设在海底，并不满足该模型的基本假设。
[0091]
更进一步的，为此本发明专利推导得出了在海底节点观测系统下自由表面相关多次波的形成模型，该模型分别描述了海底节点数据上下行波即其中所含自由表面相关多次波的形成过程，在平面波域中分别表示为：
[0092][0093][0094]
其中，p表示信号的射线参数向量(包含各个方向上射线参数的分量信息)，ω表示傅里叶变换后的信号波数，zr表示表示检波器放置的海水深度，z0表示震源激发的水深，exp()表示自然指数函数，u(p,ω,zr)表示频率平面波域的上行波信号，p(p,ω,z0)表示在海面处观测到的全波场信号，i表示复信号的虚数单位，kz表示波场的垂向波数分量，g(p,ω)
表示一次反射格林函数(即需要提取的一次波信息)，a(ω)表示海面的波阻抗(通常取为常数-1)，s(p,ω)表示震源的激发信号，d(p,ω,zr)表示频率平面波域的下行波信号，pd(p,ω,zr)表示直达波信号。
[0095]
s4：分别将分离出的上行波和下行波海底节点数据转化到频率域，进行上下行波反褶积。
[0096]
更进一步的，按数据道进行一维快速傅里叶变换，分别将分离出的平面波域上行波u
′
(p,t)和下行波数据d
′
(p,t)转化到频率域s(p,ω)，逐平面波方向、逐频率进行上下行波反褶积。
[0097]
更进一步的，根据所述分离得到的平面波域上下行波数据，以及推导得到的上下行波的数学模型，可以得到地下一次所述反射格林函数包括，
[0098][0099]
其中，p
(u)
(p,ω,zr)表示频率平面波域的上行波信号(某个射线参数上的复数据点)，p
(d)
(p,ω,zr)表示表示频率平面波域的下行波信号(某个射线参数上的复数据点)，g(p,ω,zr)表示最终求得的一次反射格林函数(即需要提取的一次波信息)。
[0100]
应说明的是，根据上式，在频率平面波域中将分离的上下行波进行反褶积，即可获得期望的地下一次反射格林函数。处理时需要对所述获得的上下行波数据逐道进行快速傅里叶变换，转换到频率平面波域后逐平面波方向、逐频率的进行上述反褶积过程。
[0101]
更进一步的，为了保证反褶积过程的稳定性，引入稳定因子，将反褶积过程近似为如下公式：
[0102][0103]
其中，表示p
(d)
(p,ω,zr)的复共轭，ε表示一个小的正数用于稳定求解，通常取0.1，其他变量的含义前已述即。
[0104]
应说明的是，根据上式即可得到稳定的反褶积结果，其中的exp(2ikzzr)项表示一个向下延拓算子，对反褶积结果进行基准面矫正，该算子可以包含在反褶积结果中。
[0105]
s5：将反褶积结果转换回时空域，得到一次反射格林函数，也即一次反射海底节点数据，实现压制。
[0106]
将所述反褶积得到的结果转换回时间-空间域，它是与s2局部平面波分解对应的反过程。
[0107]
应说明的是，所述反褶积过程消除了震源子波的影响，扩展了频带，因而转换回时间-空间域时需要计算更高的频带范围，而不是同s2一样仅考虑震源子波的频带范围。
[0108]
实施例2
[0109]
参照图3-7，为本发明的一个实施例，提供了一种基于特征反射非线性褶积模型的多次波压制方法，为了验证本发明的有益效果，通过实验进行科学论证。
[0110]
以图4所示的海底节点地震数据共检波点道集为例作为输入(同时输入压力和垂
向速度分量，二者的处理方法相似)，数据空间方向采样点数均为900道，采样间隔5米；时间采样点数为5000，采样间隔0.3毫秒；数据中包含的平面波倾角方向的范围根据海水速度确定，海水速度取1500m/s，则数据平面波分解的最大射线参数为1/1500＝0.000667s/m，设置平面波分解的射线参数采样间隔为0.000003s/m，采样数量为601个方向。
[0111]
图4的左侧为共检波点道集数据，右侧为对应的垂直井观测数据(vsp记录)，该数据为根据图3所示的模型模拟得到，因而只有一条同相轴是我们需要的来自地下界面的一次反射，根据vsp记录可以分析得到，自上往下第三条同相轴即为期望提取的一次反射，其余的同相轴大都是自由表面相关的多次波。自上往下的第一条同相轴为直达波，它虽然不是多次波，但通常并不用于偏移成像。
[0112]
使用快速傅里叶变换算法将时间-空间域的地震数据d(x,t)转换到频率-空间域d(x,ω)，根据快速傅里叶变换算法的加速特性，将原数据的时间采样点数充零至8192个(2的幂指数次方)，并按道使用多线程并行计算。根据所述确定的观测系统和射线参数，对两个分量的共检波点道集数据分别进行平面波分解，并调整数据量级。
[0113]
根据最小二乘方法构建法方程求解自适应滤波因子，求解时选择平面波域数据的前700个时间采样点数据构建方程，滤波器长度取约50个时间采样点，对于不同时间采样间隔的海底节点数据，可根据地震子波的时间长度调整相应的调整滤波器长度。
[0114]
图5表示本发明的上下行波分离示意图。上下行波分离是后续上下行波反褶积的必要前提，左上图为输入的海底节点时空域共检波点道集数据，数据中同时包含有上行波和下行波；右上图为将时空域数据转换到平面波域的结果；左下图为分离得到的上行波；右下图为分离得到的下行波。
[0115]
分离得到的上下行波将用于后续的上下行波反褶积，对分离得到的平面波域上下行波数据分别逐道进行一维快速傅里叶变换，将数据转换到频率平面波域，之后逐数据点地进行上下行波反褶积，反褶积的稳定因子取约0.001。
[0116]
反褶积得到的平面波域的一次反射格林函数转换回时间空间域。反变换是逐频率进行的，需要覆盖数据的最大频率采样范围(1到4096个频率片)。
[0117]
图6为本发明的上下行波反褶积提取一次反射格林函数示意图，左图即为在平面波域进行上下行波反褶积后得到的平面波域的一次反射格林函数；右图为将左图转换回时间-空间域后得到的多次波压制结果，可以看到与图5左上图的原数据对比，处理结果中的同相轴明显减少，数据中的多次波被有效压制，而一次波对应的同相轴被保留下来，分辨率也有所提高。
[0118]
图7为本发明处理结果与原海底节点数据共检波点道集的对比分析，左图为海底节点数据共检波点道集，并显示了数据接收器所在道的波形。右图为基准面矫正后的本发明的反褶积处理结果，并显示了接收器所在道的波形，根据图3所示的海底节点观测系统，震源和检波点位于两个不同深度的基准面上，而上下行波反褶积结果消除了这一基准面差异，故将处理结果基准面矫正回海底节点观测系统，以方便与原数据对比。通过与原数据对比可以看到，自上往下第三条同相轴(即一次反射)被保留下来，而其余的自由表面相关多次波被有效压制，此外通过抽道的波形图对比发现，处理结果中反射界面位置处为零相位子波，说明处理结果可以一定程度上消除震源子波的影响。
[0119]
本发明是一种基于特征反射非线性褶积模型的海底节点数据自由表面相关多次
波压制技术，主要用于勘探地震中海底节点观测数据的自由表面相关多次波压制。本方法首先将海底节点数据的压力分量和垂向速度分量分选为共检波点道集，通过局部平面波分解将共检波点道集数据转换到平面波域，在平面波域进行双检数据合并分离得到上下行波记录，再将上下行波记录转换到频率域，逐频率、逐平面波方向的进行上下行波反褶积，之后将反褶积结果转换回时间-空间域用于后续处理。
[0120]
应说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。
[0121]
应说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。
[0122]
本领域内的技术人员应明白，本技术的实施例可提供为方法、系统、或计算机程序产品。因此，本技术可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本技术可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。本技术实施例中的方案可以采用各种计算机语言实现，例如，面向对象的程序设计语言java和直译式脚本语言javascript等。
[0123]
本技术是参照根据本技术实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0124]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0125]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0126]
尽管已描述了本技术的优选实施例，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例作出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本技术范围的所有变更和修改。
[0127]
显然，本领域的技术人员可以对本技术进行各种改动和变型而不脱离本技术的精神和范围。这样，倘若本技术的这些修改和变型属于本技术权利要求及其等同技术的范围之内，则本技术也意图包含这些改动和变型在内。