技术领域
[0001] 本
发明涉及
地震勘探技术领域，尤其是涉及高精度反射系数求取方法和系统。
背景技术
[0002] 在
地震波传播过程中，
地层的粘弹特性会对
震源激发的地震波造成
能量衰减。对
煤炭勘探而言，比较普遍的一种现象是，上覆强反射地层阻碍地震波向下传播，影响目标
煤层的勘探和开发。通常的地震波处理方法是基于弹性
框架研究的，因此不能很好解决此类地质问题。实际地下介质并不是理想的弹性体，具有粘滞性，对地震波有很大的吸收衰减作用。
[0003] 地震成像的常用方法一般是通过求解地下介质的反射系数序列来进行成像，在计算地下介质反射系数序列的过程中必然要考虑来自地层对地震波的吸收衰减作用。目前，国内外通过计算衰减函数与地震子波获得反射系数的方法分为两大类，一种是反Q滤波，主要是沿时间方向做地震波的能量补偿，但不适用于沿波场偏移方向进行外推；另一种是反褶积方法，通过反褶积方法求解反射系数序列的过程主要是转换在
频率域中进行的，往往在计算的过程中缺少了某些时间域中的信息，不利于反射系数序列的精确求解。因此，需要一种方法结合反射系数在时间域与频率域中的变化特点来计算反射系数序列，伽伯反褶积便应运而生。传统的伽伯反褶积虽然同时考虑了反射系数在时间域与频率域中的变化特征，但是求解反射系数的精度不高，导致
地震勘探的
分辨率较低，影响地震剖面的成像精度。
发明内容
[0004] 有鉴于此，本发明的目的在于提供高精度反射系数求取方法和系统，可以提高地震勘探的分辨率和地震剖面成像精度，从而更好的为地震勘探服务。
[0005] 第一方面，本发明
实施例提供了高精度反射系数求取方法，所述方法包括：
[0006] 获取目标对象的地震记录和衰减子波；
[0007] 将所述地震记录和所述衰减子波转换到频域中，得到频域中的地震记录和频域中的衰减子波；
[0008] 根据所述频域中的地震记录和所述频域中的衰减子波，得到频域中的反射系数；
[0009] 对所述频域中的地震记录、所述频域中的衰减子波和所述频域中的反射系数做伽伯变换，得到伽伯域中的地震记录、伽伯域中的衰减子波和伽伯域中的反射系数；
[0010] 对所述衰减子波进行奇偶分解，得到衰减子波奇分量和衰减子波偶分量；
[0011] 结合高阶走时计算，对所述衰减子波奇分量和所述衰减子波偶分量中地震走时进行高阶走时计算，得到伽伯域中高阶走时衰减子波奇分量和伽伯域中高阶走时衰减子波偶分量；
[0012] 根据所述伽伯域中的反射系数、所述伽伯域中的地震记录、所述伽伯域中高阶走时衰减子波奇分量和所述伽伯域中高阶走时衰减子波偶分量，得到高阶走时反射系数奇分量和高阶走时反射系数偶分量；
[0013] 根据所述高阶走时反射系数奇分量和所述高阶走时反射系数偶分量，得到高精度反射系数序列；
[0014] 根据所述高精度反射系数序列进行地震成像，得到所述目标对象的地震剖面。
[0015] 结合第一方面，本发明实施例提供了第一方面的第一种可能的实施方式，其中，所述获取目标对象的地震记录和衰减子波，包括：
[0016] 根据地震子波、衰减函数和反射系数的褶积结果得到所述地震记录，根据所述地震子波和所述衰减函数的合并结果得到所述衰减子波；
[0017] 根据下式计算所述地震记录：
[0018] x(t)＝w(t)*a(t)*r(t)
[0019] 其中，x(t)为所述地震记录，w(t)为所述地震子波，a(t)为所述衰减函数，r(t)为所述反射系数；
[0020] 根据下式计算所述衰减子波：
[0021] wa(t)＝w(t)*a(t)
[0022] 其中，wa(t)为所述衰减子波，w(t)为所述地震子波，a(t)为所述衰减函数。
[0023] 结合第一方面的第一种可能的实施方式，本发明实施例提供了第一方面的第二种可能的实施方式，其中，所述将所述地震记录和所述衰减子波转换到频域中，得到频域中的地震记录、频域中的衰减子波，包括：
[0024] 根据下式计算所述频域中的地震记录：
[0025] x(f)＝w(f)a(f)r(f)
[0026] 其中，x(f)为所述频域中的地震记录，w(f)为频域中的地震子波，a(f)为频域中的衰减函数，r(f)为频域中的反射系数；
[0027] 根据下式计算所述频域中的衰减子波：
[0028] wa(f)＝w(f)a(f)
[0029] 其中，wa(f)为所述频域中的衰减子波，w(f)为频域中的地震子波，a(f)为频域中的衰减函数。
[0030] 结合第一方面的第二种可能的实施方式，本发明实施例提供了第一方面的第三种可能的实施方式，其中，所述根据所述频域中的地震记录和所述频域中的衰减子波，得到频域中的反射系数，包括：
[0031] 根据下式计算所述频域中的反射系数：
[0032]
[0033] 其中，r(f)为所述频域中的反射系数，x(f)为所述频域中的地震记录，wa(f)为所述频域中的衰减子波。
[0034] 结合第一方面，本发明实施例提供了第一方面的第四种可能的实施方式，其中，所述对所述频域中的地震记录、所述频域中的衰减子波和所述频域中的反射系数做伽伯变换，得到伽伯域中的地震记录、伽伯域中的衰减子波和伽伯域中的反射系数，包括：
[0035] 根据下式计算所述伽伯域中的地震记录：
[0036]
[0037] 其中，G[x](t,f)为所述伽伯域中的地震记录， 为对所述地震记录做小时窗傅里叶变换；
[0038] 根据下式计算所述伽伯域中的衰减子波：
[0039] G[wa](t,f)＝Gn[w(f)a(f)]
[0040] 其中，G[wa](t,f)为所述伽伯域中的衰减子波，w(f)为频域中的地震子波，a(f)为频域中的衰减函数；
[0041] 根据下式计算所述伽伯域中的反射系数：
[0042]
[0043] 其中，G[r](t,f)为所述伽伯域中的反射系数，G[x](t,f)为所述伽伯域中的地震记录，G[wa](t,f)为所述伽伯域中的衰减子波。
[0044] 结合第一方面，本发明实施例提供了第一方面的第五种可能的实施方式，其中，所述对所述衰减子波进行奇偶分解，得到衰减子波奇分量和衰减子波偶分量，包括：
[0045] 根据下式对所述衰减子波进行奇偶分解：
[0046] wa(t)＝2wo(t)+2we(t)
[0047] 其中，wo(t)为所述衰减子波奇分量，we(t)为所述衰减子波偶分量，wa(t)为所述衰减子波。
[0048] 结合第一方面的第五种可能的实施方式，本发明实施例提供了第一方面的第六种可能的实施方式，其中，所述结合高阶走时计算，对所述衰减子波奇分量和所述衰减子波偶分量中地震走时进行高阶走时计算，得到伽伯域中高阶走时衰减子波奇分量和伽伯域中高阶走时衰减子波偶分量，包括：
[0049] 根据下式得到所述伽伯域中高阶走时衰减子波奇分量和所述伽伯域中高阶走时衰减子波偶分量：
[0050] G[wa](t,f)＝G[wo](t,f)+G[we](t,f)
[0051] 其中，G[wo](t,f)为所述伽伯域中高阶走时衰减子波奇分量，G[we](t,f)为所述伽伯域中高阶走时衰减子波偶分量，G[wa](t,f)为所述伽伯域中的衰减子波。
[0052] 结合第一方面的第六种可能的实施方式，本发明实施例提供了第一方面的第七种可能的实施方式，其中，所述高阶走时计算是将地震波的走时做高阶展开得到，包括：
[0053] 根据下式得到所述高阶走时计算：
[0054]
[0055] 其中， 为射线参数，v为地下介质速度，t0为初始时间。
[0056] 结合第一方面，本发明实施例提供了第一方面的第八种可能的实施方式，其中，所述根据所述伽伯域中的反射系数、所述伽伯域中的地震记录、所述伽伯域中高阶走时衰减子波奇分量和所述伽伯域中高阶走时衰减子波偶分量，得到高阶走时反射系数奇分量和高阶走时反射系数偶分量，包括：
[0057] 根据下式计算所述高阶走时反射系数奇分量和所述高阶走时反射系数偶分量：
[0058]
[0059] 其中，G[ro](t,f)为所述高阶走时反射系数奇分量，G[re](t,f)为所述高阶走时反射系数偶分量，G[x](t,f)为所述伽伯域中的地震记录，G[wa](t,f)为所述伽伯域中的衰减子波，G[wo](t,f)为所述伽伯域中高阶走时衰减子波奇分量，G[we](t,f)为所述伽伯域中高阶走时衰减子波偶分量，G[r](t,f)为所述伽伯域中的反射系数。
[0060] 第二方面，本发明实施例还提供高精度反射系数求取系统，所述系统包括：
[0061] 获取单元，用于获取目标对象的地震记录和衰减子波；
[0062] 第一转换单元，用于将所述地震记录和所述衰减子波转换到频域中，得到频域中的地震记录和频域中的衰减子波；
[0063] 反射系数获取单元，用于根据所述频域中的地震记录和所述频域中的衰减子波，得到频域中的反射系数；
[0064] 第二转换单元，用于对所述频域中的地震记录、所述频域中的衰减子波和所述频域中的反射系数做伽伯变换，得到伽伯域中的地震记录、伽伯域中的衰减子波和伽伯域中的反射系数；
[0065] 奇偶分解单元，用于对所述衰减子波进行奇偶分解，得到衰减子波奇分量和衰减子波偶分量；
[0066] 第三转换单元，用于结合高阶走时计算，对所述衰减子波奇分量和所述衰减子波偶分量中地震走时进行高阶走时计算，得到伽伯域中高阶走时衰减子波奇分量和伽伯域中高阶走时衰减子波偶分量；
[0067] 第一计算单元，用于根据所述伽伯域中的反射系数、所述伽伯域中的地震记录、所述伽伯域中高阶走时衰减子波奇分量和所述伽伯域中高阶走时衰减子波偶分量，得到高阶走时反射系数奇分量和高阶走时反射系数偶分量；
[0068] 第二计算单元，用于根据所述高阶走时反射系数奇分量和所述高阶走时反射系数偶分量，得到高精度反射系数序列；
[0069] 成像单元，用于根据所述高精度反射系数序列进行地震成像，得到所述目标对象的地震剖面。
[0070] 本发明实施例提供了高精度反射系数求取方法和系统，包括：获取目标对象的地震记录和衰减子波；将地震记录和衰减子波转换到频域中，得到频域中的地震记录和频域中的衰减子波；根据频域中的地震记录和频域中的衰减子波，得到频域中的反射系数；对频域中的地震记录、频域中的衰减子波和频域中的反射系数做伽伯变换，得到伽伯域中的地震记录、伽伯域中的衰减子波和伽伯域中的反射系数；对衰减子波进行奇偶分解，得到衰减子波奇分量和衰减子波偶分量；结合高阶走时计算，对衰减子波奇分量和衰减子波偶分量中地震走时进行高阶走时计算，得到伽伯域中高阶走时衰减子波奇分量和伽伯域中高阶走时衰减子波偶分量；根据伽伯域中的反射系数、伽伯域中的地震记录、伽伯域中高阶走时衰减子波奇分量和伽伯域中高阶走时衰减子波偶分量，得到高阶走时反射系数奇分量和高阶走时反射系数偶分量；根据高阶走时反射系数奇分量和高阶走时反射系数偶分量，得到高精度反射系数序列；根据高精度反射系数序列进行地震成像，得到目标对象的地震剖面，可以提高地震勘探的分辨率和地震剖面成像精度，从而更好的为地震勘探服务。
[0071] 本发明的其他特征和优点将在随后的
说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点在说明书、
权利要求书以及
附图中所特别指出的结构来实现和获得。
[0072] 为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
[0073] 为了更清楚地说明本发明具体实施方式或
现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0074] 图1为本发明实施例一提供的反射系数求取方法
流程图；
[0075] 图2为本发明实施例二提供的常规伽伯反褶积偏移剖面示意图；
[0076] 图3为本发明实施例二提供的基于高精度奇偶分解伽伯反褶积偏移剖面示意图；
[0077] 图4为本发明实施例三提供的反射系数求取系统示意图。
[0078] 图标：
[0079] 10-获取单元；20-第一转换单元；30-反射系数获取单元；40-第二转换单元；50-奇偶分解单元；60-第三转换单元；70-第一计算单元；80-第二计算单元；90-成像单元。
具体实施方式
[0080] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0081] 为便于对本实施例进行理解，下面对本发明实施例进行详细介绍。
[0082] 实施例一：
[0083] 图1为本发明实施例一提供的反射系数求取方法流程图。
[0084] 参照图1，该方法包括以下步骤：
[0085] 步骤S101，获取目标对象的地震记录和衰减子波；
[0086] 步骤S102，将地震记录和衰减子波转换到频域中，得到频域中的地震记录和频域中的衰减子波；
[0087] 步骤S103，根据频域中的地震记录和频域中的衰减子波，得到频域中的反射系数；
[0088] 步骤S104，对频域中的地震记录、频域中的衰减子波和频域中的反射系数做伽伯变换，得到伽伯域中的地震记录、伽伯域中的衰减子波和伽伯域中的反射系数；
[0089] 步骤S105，对衰减子波进行奇偶分解，得到衰减子波奇分量和衰减子波偶分量；
[0090] 步骤S106，结合高阶走时计算，对衰减子波奇分量和衰减子波偶分量中地震走时进行高阶走时计算，得到伽伯域中高阶走时衰减子波奇分量和伽伯域中高阶走时衰减子波偶分量；
[0091] 步骤S107，根据伽伯域中的反射系数、伽伯域中的地震记录、伽伯域中高阶走时衰减子波奇分量和伽伯域中高阶走时衰减子波偶分量，得到高阶走时反射系数奇分量和高阶走时反射系数偶分量；
[0092] 步骤S108，根据高阶走时反射系数奇分量和高阶走时反射系数偶分量，得到高精度反射系数序列；
[0093] 步骤S109，根据高精度反射系数序列进行地震成像，得到目标对象的地震剖面。
[0094] 进一步的，步骤S101包括：
[0095] 根据地震子波、衰减函数和反射系数的褶积结果得到地震记录，根据地震子波和衰减函数的合并结果得到衰减子波。
[0096] 进一步的，步骤S101还包括：
[0097] 根据公式(1)计算地震记录：
[0098] x(t)＝w(t)*a(t)*r(t) (1)
[0099] 其中，x(t)为地震记录，w(t)为地震子波，a(t)为衰减函数，r(t)为反射系数；
[0100] 根据公式(2)计算衰减子波：
[0101] wa(t)＝w(t)*a(t) (2)
[0102] 其中，wa(t)为衰减子波，w(t)为地震子波，a(t)为衰减函数。
[0103] 进一步的，步骤S102包括：
[0104] 根据公式(3)计算频域中的地震记录：
[0105] x(f)＝w(f)a(f)r(f) (3)
[0106] 其中，x(f)为频域中的地震记录，w(f)为频域中的地震子波，a(f)为频域中的衰减函数，r(f)为频域中的反射系数；
[0107] 根据公式(4)计算频域中的衰减子波：
[0108] wa(f)＝w(f)a(f) (4)
[0109] 其中，wa(f)为频域中的衰减子波，w(f)为频域中的地震子波，a(f)为频域中的衰减函数。
[0110] 进一步的，步骤S103包括：
[0111] 根据公式(5)计算频域中的反射系数：
[0112]
[0113] 其中，r(f)为频域中的反射系数，x(f)为频域中的地震记录，wa(f)为频域中的衰减子波。
[0114] 进一步的，步骤S104包括：
[0115] 根据公式(6)计算伽伯域中的地震记录：
[0116]
[0117] 其中，G[x](t,f)为伽伯域中的地震记录， 为对地震记录做小时窗傅里叶变换；
[0118] 根据公式(7)计算伽伯域中的衰减子波：
[0119] G[wa](t,f)＝Gn[w(f)a(f)] (7)
[0120] 其中，G[wa](t,f)为伽伯域中的衰减子波，w(f)为频域中的地震子波，a(f)为频域中的衰减函数；
[0121] 根据公式(8)计算伽伯域中的反射系数：
[0122]
[0123] 其中，G[r](t,f)为伽伯域中的反射系数，G[x](t,f)为伽伯域中的地震记录，G[wa](t,f)为伽伯域中的衰减子波。
[0124] 进一步的，步骤S105包括：
[0125] 根据公式(9)对衰减子波进行奇偶分解：
[0126] wa(t)＝2wo(t)+2we(t) (9)
[0127] 其中， 为衰减子波奇分量， 为衰减子波偶分量，wa(t)为衰减子波。
[0128] 进一步的，步骤S106包括：
[0129] 根据公式(10)得到伽伯域中高阶走时衰减子波奇分量和伽伯域中高阶走时衰减子波偶分量：
[0130] G[wa](t,f)＝G[wo](t,f)+G[we](t,f) (10)
[0131] 其中，G[wo](t,f)为伽伯域中高阶走时衰减子波奇分量，G[we](t,f)为伽伯域中高阶走时衰减子波偶分量，G[wa](t,f)为伽伯域中的衰减子波。
[0132] 进一步的，步骤S106还包括：
[0133] 在
水平层状介质中，根据公式(11)计算地震波的走时t：
[0134]
[0135] 其中，t(p)为地震波的走时， 为射线参数，v为地下介质速度；
[0136] 根据公式(12)得到地震波的走时高阶展开：
[0137]
[0138] 其中，t(p)为地震波的走时， 为射线参数，v为地下介质速度，t0为初始时间。
[0139] 进一步的，步骤S107包括：
[0140] 根据公式(13)计算高阶走时反射系数奇分量和高阶走时反射系数偶分量：
[0141]
[0142] 其中，G[ro](t,f)为高阶走时反射系数奇分量，G[re](t,f)为高阶走时反射系数偶分量，G[x](t,f)为伽伯域中的地震记录，G[wa](t,f)为伽伯域中的衰减子波，G[wo](t,f)为伽伯域中高阶走时衰减子波奇分量，G[we](t,f)为伽伯域中高阶走时衰减子波偶分量，G[r](t,f)为伽伯域中的反射系数。
[0143] 实施例二：
[0144] 将奇偶分解伽伯反褶积方法所得到的地震剖面与常规伽伯反褶积方法得到的地震剖面进行对比，如图2和图3所示，走时计算经过高阶展开后，较常规计算方法可以很好地提高反射系数计算的精度；在高精度走时计算的
基础上，地震波同相轴的连续性得到了增强，地震剖面成像
质量和地震勘探分辨率可得到明显提高；经奇偶分解计算后，对消除强反射对下覆地层的屏蔽作用起到了明显改善作用。
[0145] 实施例三：
[0146] 图4为本发明实施例三提供的反射系数求取系统示意图。
[0147] 参照图4，该系统包括获取单元10、第一转换单元20、反射系数获取单元30、第二转换单元40、奇偶分解单元50、第三转换单元60、第一计算单元70、第二计算单元80和成像单元90。
[0148] 获取单元10，用于获取目标对象的地震记录和衰减子波；
[0149] 第一转换单元20，用于将地震记录和衰减子波转换到频域中，得到频域中的地震记录和频域中的衰减子波；
[0150] 反射系数获取单元30，用于根据频域中的地震记录和频域中的衰减子波，得到频域中的反射系数；
[0151] 第二转换单元40，用于对频域中的地震记录、频域中的衰减子波和频域中的反射系数做伽伯变换，得到伽伯域中的地震记录、伽伯域中的衰减子波和伽伯域中的反射系数；
[0152] 奇偶分解单元50，用于对衰减子波进行奇偶分解，得到衰减子波奇分量和衰减子波偶分量；
[0153] 第三转换单元60，用于结合高阶走时计算，对衰减子波奇分量和衰减子波偶分量中地震走时进行高阶走时计算，得到伽伯域中高阶走时衰减子波奇分量和伽伯域中高阶走时衰减子波偶分量；
[0154] 第一计算单元70，用于根据伽伯域中的反射系数、伽伯域中的地震记录、伽伯域中高阶走时衰减子波奇分量和伽伯域中高阶走时衰减子波偶分量，得到高阶走时反射系数奇分量和高阶走时反射系数偶分量；
[0155] 第二计算单元80，用于根据高阶走时反射系数奇分量和高阶走时反射系数偶分量，得到高精度反射系数序列；
[0156] 成像单元90，用于根据高精度反射系数序列进行地震成像，得到目标对象的地震剖面。
[0157] 本发明实施例提供了高精度反射系数求取方法和系统，包括：获取目标对象的地震记录和衰减子波；将地震记录和衰减子波转换到频域中，得到频域中的地震记录和频域中的衰减子波；根据频域中的地震记录和频域中的衰减子波，得到频域中的反射系数；对频域中的地震记录、频域中的衰减子波和频域中的反射系数做伽伯变换，得到伽伯域中的地震记录、伽伯域中的衰减子波和伽伯域中的反射系数；对衰减子波进行奇偶分解，得到衰减子波奇分量和衰减子波偶分量；结合高阶走时计算，对衰减子波奇分量和衰减子波偶分量中地震走时进行高阶走时计算，得到伽伯域中高阶走时衰减子波奇分量和伽伯域中高阶走时衰减子波偶分量；根据伽伯域中的反射系数、伽伯域中的地震记录、伽伯域中高阶走时衰减子波奇分量和伽伯域中高阶走时衰减子波偶分量，得到高阶走时反射系数奇分量和高阶走时反射系数偶分量；根据高阶走时反射系数奇分量和高阶走时反射系数偶分量，得到高精度反射系数序列；根据高精度反射系数序列进行地震成像，得到目标对象的地震剖面，可以提高地震勘探的分辨率和地震剖面成像精度，从而更好的为地震勘探服务。
[0158] 本发明实施例还提供一种
电子设备，包括
存储器、处理器及存储在存储器上并可在处理器上运行的
计算机程序，处理器执行计算机程序时实现上述实施例提供的高精度反射系数求取方法的步骤。
[0159] 本发明实施例还提供一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，计算机程序被处理器运行时执行上述实施例的高精度反射系数求取方法的步骤。
[0160] 本发明实施例所提供的计算机程序产品，包括存储了程序代码的计算机可读存储介质，所述程序代码包括的指令可用于执行前面方法实施例中所述的方法，具体实现可参见方法实施例，在此不再赘述。
[0161] 所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统和装置的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
[0162] 另外，在本发明实施例的描述中，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
[0163] 所述功能如果以
软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个处理器可执行的非易失的计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，
服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动
硬盘、
只读存储器(ROM，Read-Only Memory)、
随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
[0164] 在本发明的描述中，需要说明的是，术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或
位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”、“第三”仅用于描述目的，而不能理解为指示或暗示相对重要性。
[0165] 最后应说明的是：以上所述实施例，仅为本发明的具体实施方式，用以说明本发明的技术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉
本技术领域的技术人员在本发明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行
修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。
