反射纵波与反射转换横波的振幅匹配方法及装置与流程

本发明涉及石油勘探领域，尤其是地震数据处理及解释技术，具体涉及一种反射纵波与反射转换横波的振幅匹配方法。

背景技术：

多波多分量地震资料联合解释技术是探测石油、天然气资源的一项新技术。可通过布置在地面检波器接收到的反射纵波信号和反射转换横波信号联合描述地下地层中石油、天然气储集体。当纵波震源激发，三分量检波器接收多波信号时，理论上，不考虑地层吸收衰减，转换波反射能量是纵波反射能量的50-60％，但实际地震资料采集中，由于受地表吸收衰减及噪声的影响，转换横波反射能量仅为纵波能量的10％-20％，同时后续地震资料处理也会造成纵波和转换横波叠后地震剖面不能保持两者振幅强弱关系。而地震振幅信息包含地震波动力学特征，是地震波最重要的属性之一。地层孔隙发育程度、烃类饱和度和储层厚度都会引起不同的纵波、转换横波地震振幅响应特征，解释人员通常利用纵波、转换横波叠后地震剖面储层段振幅信息建立气层反射纵波和反射转换横波振幅强弱关系，因此反射纵波、反射转换横波振幅信息直接影响了气层预测精度。

现有技术中对于多波资料的振幅关系校正处于探索阶段，一般是利用一个标准层的纵波和转换横波振幅匹配关系校正整个转换波数据体，该方法有着两个明显的缺陷：第一，利用一个标准层的纵波和转换横波振幅匹配关系校正整个转换波数据体，对整个数据体采用一个振幅匹配因子，能够满足标准层附近转换波地震振幅的校正，对远离标准层的同相轴恢复的地震振幅与理论振幅存在偏差，利用该数据进行后续储层预测，造成储层预测精度偏低。第二，该方法操作繁琐，选择一次标准层便得到一个转换波数据体。例如：，若要恢复五层振幅关系必须选择五次标准层，而且得到五个转换波数据体，故难以满足实际生产需要。

技术实现要素：

针对现有技术中的问题，本发明所提供的反射纵波与反射转换横波的振幅匹配方法，可以保证匹配后转换波振幅的可靠性，减少了单层振幅匹配时远离标准层振幅匹配结果与理论振幅值误差，提高实际生产工作效率，可以一次性恢复浅、中、深层转换波数据体地震振幅，避免了多次标准层的选择，为后续利用振幅信息进行储层预测提供可靠数据体。

为解决上述技术问题，本发明提供以下技术方案：

第一方面，本发明提供一种反射纵波与反射转换横波的振幅匹配方法，包括：

分别读取反射纵波的合成地震记录及反射转换横波的合成地震记录中井中多个标准层位的第一同相轴振幅；

分别读取反射纵波的实际地震剖面及反射转换横波的实际地震剖面中井中所述多个标准层的第二同相轴振幅；

根据多个第一同相轴振幅及多个第二同相轴振幅生成时变振幅匹配因子；

根据所述时变振幅匹配因子及所述反射转换横波的实际地震剖面，将所述反射转换横波的振幅与所述反射纵波的振幅进行匹配。

一实施例中，反射纵波与反射转换横波的振幅匹配方法还包括：对地震数据进行预处理，生成所述反射纵波和反射转换横波。

一实施例中，反射纵波与反射转换横波的振幅匹配方法还包括：分别对所述反射纵波目的时窗内的反射纵波及反射转换横波目的时窗内的反射转换横波进行频谱分析，生成所述反射纵波的主频数据及所述反射转换横波的主频数据。

一实施例中，反射纵波与反射转换横波的振幅匹配方法还包括：根据所述反射纵波的主频数据、所述反射转换横波的主频数据分别生成所述反射纵波的地震子波及所述反射转换横波的地震子波。

一实施例中，反射纵波与反射转换横波的振幅匹配方法还包括：根据所述井的测井数据计算所述反射纵波的反射系数及所述反射转换横波的反射系数。

一实施例中，所述测井数据包括纵波速度、横波速度及测井密度。

一实施例中，反射纵波与反射转换横波的振幅匹配方法还包括：将所述反射纵波的反射系数与所述反射纵波的地震子波进行褶积，生成所述反射纵波的合成地震记录；

将所述反射转换横波的反射系数与反射转换横波的地震子波进行褶积，生成所述反射转换横波的合成地震记录。

一实施例中，反射纵波与反射转换横波的振幅匹配方法还包括：所述根据多个第一同相轴振幅及多个第二同相轴振幅生成时变振幅匹配因子，包括：

利用最小二乘法，根据多个第一同相轴振幅及多个第二同相轴振幅生成所述井的多个标准层位的标准振幅匹配因子；

利用三次样条差值算法，根据所述井的多个标准层位的振幅匹配因子生成所述时变振幅匹配因子。

一实施例中，反射纵波与反射转换横波的振幅匹配方法还包括：所述根据所述时变振幅匹配因子及所述反射转换横波的实际地震剖面，将所述反射转换横波的振幅与所述反射纵波的振幅进行匹配，包括：

将所述反射转换横波的实际地震剖面的每个共转换点道集乘以所述时变振幅匹配因子，以将所述反射转换横波的振幅与所述反射纵波的振幅进行匹配。

第二方面，本发明提供一种反射纵波与反射转换横波的振幅匹配装置，该装置包括：

第一读取单元，用于分别读取反射纵波的合成地震记录及反射转换横波的合成地震记录中井中多个标准层位的第一同相轴振幅；

第二读取单元，用于分别读取反射纵波的实际地震剖面及反射转换横波的实际地震剖面中井中所述多个标准层的第二同相轴振幅；

因子生成单元，用于根据多个第一同相轴振幅及多个第二同相轴振幅生成时变振幅匹配因子；

匹配单元，用于根据所述时变振幅匹配因子及所述反射转换横波的实际地震剖面，将所述反射转换横波的振幅与所述反射纵波的振幅进行匹配。

一实施例中，反射纵波与反射转换横波的振幅匹配装置还包括：预处理单元，用于对地震数据进行预处理，生成所述反射纵波和反射转换横波。

一实施例中，反射纵波与反射转换横波的振幅匹配装置还包括：主频数据生成单元，用于分别对所述反射纵波目的时窗内的反射纵波及反射转换横波目的时窗内的反射转换横波进行频谱分析，生成所述反射纵波的主频数据及所述反射转换横波的主频数据。

一实施例中，反射纵波与反射转换横波的振幅匹配装置还包括：子波生成单元，用于根据所述反射纵波的主频数据、所述反射转换横波的主频数据分别生成所述反射纵波的地震子波及所述反射转换横波的地震子波。

一实施例中，反射纵波与反射转换横波的振幅匹配装置还包括反射系数生成单元，用于根据所述井的测井数据计算所述反射纵波的反射系数及所述反射转换横波的反射系数。

一实施例中，所述测井数据包括纵波速度、横波速度及测井密度。

一实施例中，反射纵波与反射转换横波的振幅匹配装置还包括：

纵波合成生成单元，用于将所述反射纵波的反射系数与所述反射纵波的地震子波进行褶积，生成所述反射纵波的合成地震记录；

横波合成生成单元，用于将所述反射转换横波的反射系数与反射转换横波的地震子波进行褶积，生成所述反射转换横波的合成地震记录。

一实施例中，因子生成单元包括：

标准因子生成模块，用于利用最小二乘法，根据多个第一同相轴振幅及多个第二同相轴振幅生成所述井的多个标准层位的标准振幅匹配因子；

时变因子生成模块，用于利用三次样条差值算法，根据所述井的多个标准层位的振幅匹配因子生成所述时变振幅匹配因子。

一实施例中，匹配单元具体用于将所述反射转换横波的实际地震剖面的每个共转换点道集乘以所述时变振幅匹配因子，以将所述反射转换横波的振幅与所述反射纵波的振幅进行匹配。

第三方面，本发明提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行程序时实现反射纵波与反射转换横波的振幅匹配方法的步骤。

第四方面，本发明提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现反射纵波与反射转换横波的振幅匹配方法的步骤。

从上述描述可知，本发明实施例提供的反射纵波与反射转换横波的振幅匹配方法及装置，首先读取反射纵波及反射转换横波地震记录上的多个标准层位的同相轴振幅，并读取该多个标准层位的实际地震剖面的同相轴振幅，以此来生成时变振幅匹配因子，最终将反射转换横波的振幅与反射纵波的振幅进行匹配。本方法保证了匹配后反射转换横波振幅的可靠性，减少了单层振幅匹配时远离标准层振幅匹配结果与理论振幅值误差，提高实际生产工作效率，一次性恢复浅中深层转换波数据体地震振幅，避免了多次标准层的选择，为后续利用振幅信息进行储层预测提供可靠数据体。

综上，本发明实施例提供的方法利用多个标准层反射纵波、反射转换横波合成地震记录振幅比例关系校正实际地震记录纵波和转换横波振幅关系，计算出单层振幅比例因子后，通过三次样条差值计算出整个数据体的振幅匹配因子(整个数据体纵向上采用不同的振幅匹配因子)，考虑了时变和空变。能够避免现有技术中振幅匹配后远离标准层振幅匹配结果与理论振幅值误差，能够一次性恢复浅中深层转换波数据体地震振幅，提高实际生产工作效率。具有的优点如下：采用多个标准层、三次样条差值算法进行时变和空变差值，匹配结果考虑了时变和空变。

附图说明

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

图1为本发明的实施例中的反射纵波与反射转换横波的振幅匹配方法的流程示意图一；

图2为本发明的实施例中的反射纵波与反射转换横波的振幅匹配方法的流程示意图二；

图3为本发明的实施例中的反射纵波与反射转换横波的振幅匹配方法的流程示意图三；

图4为本发明的实施例中的反射纵波与反射转换横波的振幅匹配方法的流程示意图四；

图5为本发明的实施例中反射纵波理论地震子波组成示意图；

图6为本发明的实施例中的反射纵波与反射转换横波的振幅匹配方法的流程示意图五；

图7为本发明的实施例中测井数据组成示意图；

图8为本发明的实施例中的反射纵波与反射转换横波的振幅匹配方法的流程示意图六；

图9为本发明的实施例中步骤300的流程示意图；

图10为本发明的实施例中步骤400的流程示意图；

图11为本发明的具体应用实例中反射纵波与反射转换横波的振幅匹配方法的流程示意图；

图12为本发明的具体应用实例中反射纵波与反射转换横波的振幅匹配方法的思路流程图；

图13为本发明的具体应用实例中蓬莱107井的pp波正演合成记录及实际pp波地震剖面图；

图14为本发明的具体应用实例中蓬莱111井的pp波正演合成记录及实际pp波地震剖面图；

图15为本发明的具体应用实例中蓬莱9井的pp波正演合成记录及实际pp波地震剖面图；

图16为本发明的具体应用实例中蓬莱107井的ps波正演合成记录及实际ps波地震剖面图；

图17为本发明的具体应用实例中蓬莱111井的ps波正演合成记录及实际ps波地震剖面图；

图18为本发明的具体应用实例中蓬莱9井的ps波正演合成记录及实际ps波地震剖面图；

图19为本发明的具体应用实例中目的层段的反射纵波频谱分析图；

图20为本发明的具体应用实例中目的层段的反射转换横波频谱分析图；

图21为本发明的具体应用实例中振幅匹配前反射纵波地震剖面图；

图22为本发明的具体应用实例中振幅匹配前反射转换横波地震剖面图；

图23为本发明的具体应用实例中振幅匹配后反射纵波地震剖面图；

图24为本发明的具体应用实例中振幅匹配后反射转换横波地震剖面图；

图25为本发明的实施例中的反射纵波与反射转换横波的振幅匹配装置的结构示意图；

图26为本发明的实施例中的电子设备的结构示意图。

具体实施方式

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

现有技术中，对多波资料的振幅关系校正一般是利用已知井目的层附近标准层纵波和转换横波理论振幅关系来校正实际地震数据振幅比例关系。该方法存在以下缺点：第一：该方法仅仅利用一个标准层的纵波和转换横波振幅匹配关系校正整个转换波数据体，对整个数据体采用一个振幅匹配因子，能够满足标准层附近转换波地震振幅的校正，对远离标准层的同相轴恢复的地震振幅与理论振幅存在偏差，利用该数据进行后续储层预测，造成储层预测精度偏低。第二：该方法恢复几层振幅关系必须选择几次标准层，(选择一次标准层得到一个转换波数据体，选择几次就得到几个转换波数据体)，不满足实际生产需要。

针对现有技术中存在的不足，本发明的实施例提供一种反射纵波与反射转换横波的振幅匹配方法的具体实施方式，参见图1，该方法具体包括如下内容：

步骤100：分别读取反射纵波的合成地震记录及反射转换横波的合成地震记录中井中多个标准层位的第一同相轴振幅。

可以理解的是，步骤100中的合成地震记录是用声波测井或垂直地震剖面资料经过人工合成转换成的地震记录(地震道)，是地震模型技术中的一种，也是层位标定、油藏描述等工作的基础，是把地质模型转化为地震信息的中间媒介。合成地震记录是联合高分辨率的测井信息与区域性的地震信息的桥梁，其精度直接影响到地质层位的准确标定。

步骤200：分别读取反射纵波的实际地震剖面及反射转换横波的实际地震剖面中井中所述多个标准层的第二同相轴振幅。

可以理解的是，步骤200中的反射转换横波是指常规爆炸震源激发的纵波以非零入射角入射到反射界面时产生反射转换横波，又称为转换横波或横波。

在步骤100及步骤200中，通过读取多个标准层的同相轴振幅，可以减少单层振幅匹配时，远离标准层振幅匹配结果与理论振幅值误差，可以一次性恢复浅、中、深层转换波数据体地震振幅，避免了多次标准层的选择。

步骤300：根据多个第一同相轴振幅及多个第二同相轴振幅生成时变振幅匹配因子。

步骤400：根据所述时变振幅匹配因子及所述反射转换横波的实际地震剖面，将所述反射转换横波的振幅与所述反射纵波的振幅进行匹配。

可以理解的是，将时变振幅匹配因子乘以反射转换横波的实际地震剖面的每个共转换点道集即可完成反射转换横波的振幅与反射纵波的振幅之间匹配。

从上述描述可知，本发明提供反射纵波与反射转换横波的振幅匹配方法，首先读取反射纵波及反射转换横波地震记录上的多个标准层位的同相轴振幅，并读取该多个标准层位的实际地震剖面的同相轴振幅，以此来生成时变振幅匹配因子，最终将反射转换横波的振幅与反射纵波的振幅进行匹配。本方法保证了匹配后反射转换横波振幅的可靠性，减少了单层振幅匹配时远离标准层振幅匹配结果与理论振幅值误差，提高实际生产工作效率，一次性恢复浅中深层转换波数据体地震振幅，避免了多次标准层的选择，为后续利用振幅信息进行储层预测提供可靠数据体。

一实施例中，参见图2，反射纵波与反射转换横波的振幅匹配方法还包括：

步骤500：对地震数据进行预处理，生成所述反射纵波和反射转换横波。

步骤500中的地震数据(或称为地震记录、地震资料、地震剖面)是指地面激发、地面接收的观测方式，即炮点设置在地面上且检波点也设在地面。也就是说，地震记录是通过地面激发，接收来自地下地层的反射信息，经过处理就得出一张张地震剖面图，随后通过地质解释，来推断地下岩层的性质和含油气情况。

一实施例中，参见图3，反射纵波与反射转换横波的振幅匹配方法还包括：

步骤600：分别对所述反射纵波目的时窗内的反射纵波及反射转换横波目的时窗内的反射转换横波进行频谱分析，生成所述反射纵波的主频数据及所述反射转换横波的主频数据。

可以理解的是，步骤600中的时窗是指由于地震资料通常以时间作为纵向的标尺展示(地震波接收地下岩层反射信号以时间标度，乘以速度，可以转换为地下的深度)，对目的层的纵向范围可以用时窗来表述。

一实施例中，参见图4，反射纵波与反射转换横波的振幅匹配方法还包括：

步骤700：根据所述反射纵波的主频数据、所述反射转换横波的主频数据分别生成所述反射纵波的地震子波及所述反射转换横波的地震子波。

步骤700中的地震子波是指一段具有确定的起始时间、能量有限且有一定延续长度的信号，它是地震记录中的基本单元。地震震源激发时所产生的地震波仅是一个延续时间极短的尖脉冲，随着尖脉冲在粘弹性介质中传播，尖脉冲的高频成分很快衰减，波形随之增长，便形成了地震子波，一个地震子波一般有2至3个相位的延续长度，大约有90ms左右，然后以地震子波的形式在地下传播。

优选的，步骤700为利用零相位子波、最小相位子波和混合相位子波参数中的一个以及主频和预设的子波长度制作反射纵波的地震子波(反射纵波理论地震子波)。其中，预设的子波长度不大于所述反射纵波目的层时窗长度的1/3。

利用零相位子波、最小相位子波和混合相位子波参数中的一个以及主频和预设的子波长度制作反射转换横波的地震子波(反射转换横波理论地震子波)。其中，预设的子波长度不大于反射转换横波目的层时窗长度的1/3。

需要说明的是，目的层可以是指为完成某项勘探任务所选定的探明层位。

参见图5，可以理解的是，反射纵波理论地震子波可以为反射纵波雷克子波、井子波、带通子波，统计子波或者克劳德子波。

一实施例中，参见图6，反射纵波与反射转换横波的振幅匹配方法还包括：

步骤800：根据所述井的测井数据计算所述反射纵波的反射系数及所述反射转换横波的反射系数。

测井数据又称测井资料，是指利用电、磁、声、热、核等物理原理制造的各种测井仪器，由测井电缆下入井内，使地面电测仪沿着井筒连续记录随深度变化的各种参数。

一实施例中，参见图7，步骤800中的测井数据包括纵波速度、横波速度及测井密度。

一实施例中，参见图8，反射纵波与反射转换横波的振幅匹配方法还包括：

步骤900a：将所述反射纵波的反射系数与所述反射纵波的地震子波进行褶积，生成所述反射纵波的合成地震记录。

步骤900b：将所述反射转换横波的反射系数与反射转换横波的地震子波进行褶积，生成所述反射转换横波的合成地震记录。

一实施例中，参见图9，步骤300包括：

步骤301：利用最小二乘法，根据多个第一同相轴振幅及多个第二同相轴振幅生成所述井的多个标准层位的标准振幅匹配因子。

步骤302：利用三次样条差值算法，根据所述井的多个标准层位的振幅匹配因子生成所述时变振幅匹配因子。

可以理解的是，步骤301及步骤302可以减少单层振幅匹配时，远离标准层振幅匹配结果与理论振幅值的误差，以提高实际生产工作效率，一次性恢复浅中深层转换波数据体地震振幅，避免了多次标准层的选择，为后续利用振幅信息进行储层预测提供可靠数据体。

一实施例中，参见图10，步骤400包括：

步骤401：将所述反射转换横波的实际地震剖面的每个共转换点道集乘以所述时变振幅匹配因子，以将所述反射转换横波的振幅与所述反射纵波的振幅进行匹配。

从上述描述可知，本发明实施例提供的反射纵波与反射转换横波的振幅匹配方法，首先读取反射纵波及反射转换横波地震记录上的多个标准层位的同相轴振幅，并读取该多个标准层位的实际地震剖面的同相轴振幅，以此来生成时变振幅匹配因子，最终将反射转换横波的振幅与反射纵波的振幅进行匹配。本方法保证了匹配后反射转换横波振幅的可靠性，减少了单层振幅匹配时远离标准层振幅匹配结果与理论振幅值误差，提高实际生产工作效率，一次性恢复浅中深层转换波数据体地震振幅，避免了多次标准层的选择，为后续利用振幅信息进行储层预测提供可靠数据体。

综上，本发明实施例提供的方法利用多个标准层反射纵波、反射转换横波合成地震记录振幅比例关系校正实际地震记录纵波和转换横波振幅关系，计算出单层振幅比例因子后，通过三次样条差值计算出整个数据体的振幅匹配因子(整个数据体纵向上采用不同的振幅匹配因子)，考虑了时变和空变。能够避免现有技术中振幅匹配后远离标准层振幅匹配结果与理论振幅值误差，能够一次性恢复浅中深层转换波数据体地震振幅，提高实际生产工作效率。具有的优点如下：采用多个标准层、三次样条差值算法进行时变和空变差值，匹配结果考虑了时变和空变。

为进一步地说明本方案，本发明以四川盆地某地区三维三分量工区为具体应用实例，运用本发明示例性实施例的反射纵波与反射转换横波的振幅匹配方法对纵波、转换横波实际地震数据振幅关系进行校正，参见图11及图12。其中，该工区有3口井，即n＝3。

s0：生成反射纵波和反射转换横波数据体。

以本领域公知的方法对该工区进行野外勘探采集，以得到原始采集的地震数据，经过地震资料预处理，得到纵波和转换横波地震偏移数据体(即实际地震数据)。

参见图13至图18，其中图13中的a、图14中的c、图15中的e、图16中的a、图17中的c及图18中的e为合成地震记录的重复8道显示图，图13中的a为蓬莱107井pp波正演合成记录，图13中的b为过蓬莱107井实际pp波地震剖面；图14中的c为蓬莱111井pp波正演合成记录，图14中的d为过蓬莱111井实际pp波地震剖面；图15中的e为蓬莱9井pp波正演合成记录，图15中的f为过蓬莱9井实际pp波地震剖面；图16中的a为蓬莱107井ps波正演合成记录，图16中的b为过蓬莱107井实际ps波地震剖面；图17中的c为蓬莱111井ps波正演合成记录，图17中的d为过蓬莱111井实际ps波地震剖面；图18中的e为蓬莱9井ps波正演合成记录，图18中的f为过蓬莱9井实际ps波地震剖面。可以理解的是，上述pp波为反射纵波，ps波为反射转换横波。过井地震剖面是指根据钻井实际大地坐标位置投影到地震测网上对应的地震剖面。该工区3口井的过井地震剖面如图13中的b、图14中的d、图15中的f、图16中的b、图17中的d及图18中的f所示。

需要说明的是，三维地震勘探分为inline、xline和time三个方向，地震剖面实际为inline或xline方向与time方向的切片。地震剖面上横轴数字是在地震资料处理中定义三维观测系统而进行的编号，即道号，在本具体应用实例中，图13-图18中的横坐标上的数字表示道号。

首先利用研究区所有钻井资料纵波速度、横波速度和密度信息制作纵波、转换横波合成地震记录；选取多个标准层，并读取每个标准层在合成地震记录和井位置处的实际地震记录振幅值，利用振幅匹配公式计算出h个标准层振幅匹配因子k(h)，并将k(h)作为插值节点，采用三次样条插值算法生成时变的振幅匹配因子kj，并将该时变匹配因子kj应用于转换横波，实现纵波和转换横波全数据的振幅匹配工作。

对工区原始采集的地震数据进行预处理，得到反射纵波地震数据和反射转换横波地震数据；对反射纵波目的层时窗内的反射纵波地震数据进行频谱分析，得到反射纵波目的层时窗内反射纵波地震数据的振幅谱，并获取第一主频；利用第一主频创建反射纵波理论地震子波；根据所述工区内测井数据中n口井的纵波速度和密度计算得到n口井的反射纵波反射系数，并将n口井的反射纵波反射系数分别与所述反射纵波理论地震子波进行褶积得到n口井的反射纵波合成地震记录，其中，n≥1且为正整数；对所述反射转换横波目的层时窗内的反射转换横波地震数据进行频谱分析，得到反射转换横波目的层时窗内反射转换横波地震数据的振幅谱，并获取第二主频；利用第二主频创建反射转换横波理论地震子波；根据所述工区内测井数据中n口井的纵波速度、横波速度和密度计算得到n口井的反射转换横波反射系数，并将n口井的反射转换横波反射系数分别与所述反射转换波理论地震子波进行褶积得到n口井的反射转换横波合成地震记录，其中，n≥1且为正整数；选取h个标准层(h≥2且为正整数)；读取h个标准层n口井合成地震记录上同相轴振幅apphwell(1)、apphwell(2)、…、apphwell(n)；apshwell(1)、apshwell(2)、…、apshwell(n)；在所述反射纵波和转换横波纯波地震数据上分别读取第h个标准层n口井井位处实际地震剖面同相轴振幅apphseis(1)、apphseis(2)、…、apphseis(n)；apshseis(1)、apshseis(2)、…、apshseis(n)。根据下面的等式1，使用最小二乘法求取第h个标准层反射纵波和反射转换横波的振幅匹配因子k(h)，采用通用的三次样条插值算法生成时变的振幅匹配因子kj；将实际转换横波地震数据每个ccp道集乘以时变匹配因子kj，完成纵波、转换横波振幅匹配工作。

s1：获取反射纵波及反射转换横波的实际地震记录主频。

根据反射纵波实际地震数据，选取目的层时窗范围内的纵波地震数据进行全数据的频谱分析，制作整个数据体振幅谱，得到纵波实际地震记录主频。具体地，根据纵波实际纯波地震数据，选取目的层以上1～10000ms至以下1～10000ms的时窗制作整个数据体振幅谱，得到纵波实际地震记录主频。图19示出了研究区纵波实际地震数据在纵波目的层段内的振幅谱，由图19可以看出，纵波地震记录主频为45hz。

同样地，对反射转换横波地震数据目的层时窗内全部转换横波地震数据进行频谱分析，制作整个数据体振幅谱，并从振幅谱上获取转换横波地震主频。图20是研究区转换横波目的层时窗内反射转换横波地震数据的振幅谱，由图20可以看出，反射转换横波地震记录主频为25hz。

s2：生成理论反射纵波的理论雷克子波及反射转换横波的理论雷克子波。

利用获得的反射纵波地震记录主频参数，选用零相位或最小相位或混合相位参数；选用1～10000ms的时间长度制作理论纵波地震子波,如雷克子波。具体地，利用步骤s1中获取的反射纵波地震记录主频参数(即45hz、零相位子波以及100ms的时间长度制作反射纵波的理论雷克子波。

利用获得的反射转换横波地震记录主频参数，选用零相位或最小相位或混合相位相位参数；选用1～10000ms的时间长度制作理论转换横波地震子波,如雷克子波。采用主频25hz、零相位子波以及100ms的子波时间长度制作得到的反射转换横波的理论雷克子波。

s3：分别生成反射纵波及反射转换横波的合成地震记录。

根据该工区3口井的测井数据中的纵波速度和密度分别计算3口井各自对应的反射纵波反射系数，并基于褶积模型，将3口井的反射纵波的反射系数分别纵波理论雷克子波进行褶积得到3口井的纵波合成地震记录，其中，图13中的a为蓬莱107井反射纵波合成地震记录；图14中的c为蓬莱111反射纵波合成地震记录；图15中的e为蓬莱9井反射纵波合成地震记录。需要说明的是，在图13中的a、图14中的c和图15中的e中以及图16中的a、图17中的c和图18中的e中，是将合成地震记录重复8道显示，例如，褶积得到的蓬莱107井反射纵波合成地震记录是1道，图13中的a是将褶积得到的蓬莱107井反射纵波合成地震记录重复显示8道而成。

根据该工区内3口井的测井数据中的纵波速度、横波速度和密度计算得到3口井各自的反射转换横波反射系数，并基于褶积模型，将3口井的反射转换横波的反射系数分别反射转换横波理论雷克子波进行褶积得到3口井的转换横波合成地震记录。图16中的a是蓬莱107的ps合成地震记录；图17中的c是蓬莱111的ps合成地震记录；图18中的e是蓬莱9井的ps合成地震记录。

s4：选取标准层。

按照石油天然气行业标准sy/t5481-2009选取标准层，在本具体应用实例中，该工区选取侏罗系下统东岳庙组(j1d)和三叠系上统须家河组三段底界(t3x³)为标准层，即图13至图18中t3x³所示例的黑色横线为标准层。

在本具体应用实例中，地震标准层是指波形特征明显、稳定，并在区域内大多数地段可连续追踪的与勘探目的层相联系的地震界面。可以理解的是，在本具体应用实例中，选取的标准层数量不小于2。

s5：读取合成地震记录上的标准层同相轴振幅。

根据3口井的反射纵波和反射转换横波合成地震记录，获取标准层同相轴振幅。其中，从地震合成记录上根据标准层所对应的时间刻度(深度)，获得相应的振幅。从图13中的a所示的蓬莱107井纵波合成地震记录得到的第一个标准层j1d的振幅值app1well(1)为0.13和第二个标准层t3x3的振幅值app2well(1)为0.07；从图14中的c所示蓬莱111井纵波合成地震记录得到的第一个标准层j1d的振幅值app1well(2)为0.16和第二个标准层t3x3的振幅值app2well(2)为0.04；从图15中的e所示蓬莱9井纵波合成地震记录得到的第一个标准层j1d的振幅值app1well(3)为0.14和第二个标准层t3x3的振幅值app2well(3)为0.04；从图16中的a所示的蓬莱107井转换横波合成地震记录得到的第一个标准层j1d的振幅值aps1well(1)为0.07和第二个标准层t3x3的振幅值aps2well(1)为0.08；从图17中的c所示的蓬莱111井转换横波合成地震记录得到的第一个标准层j1d的振幅值aps1well(2)为0.1和第二个标准层t3x3的振幅值aps2well(2)为0.06；从图18中的e所示的蓬莱9井转换横波合成地震记录得到的第一个标准层j1d的振幅值aps1well(3)为0.21和第二个标准层t3x3的振幅值aps2well(3)为0.03。

s6：读取实际地震数据剖面上的标准层同相轴振幅。

在反射纵波和反射转换横波实际地震数据上分别读取3口井井位处实际地震标准层同相轴振幅。从图13中的b所示过蓬莱107井的实际纵波地震剖面上获得第一个标准层j1d的振幅值app1seis(1)为21780和第二个标准层t3x3的振幅值app2seis(1)为9102；从图14中的d所示过蓬莱111井的实际纵波地震剖面上获得第一个标准层j1d的振幅值app1seis(2)为14154和第二个标准层t3x3的振幅值app2seis(2)7485；从图15中的f所示过蓬莱9井的实际纵波地震剖面获得第一个标准层j1d的振幅值app1seis(3)为15222和第二个标准层t3x3的振幅值app2seis(3)为9102；从图16中的b所示过蓬莱107井的实际转换横波地震剖面获得的第一个标准层j1d的振幅值aps1seis(1)为0.005和第二个标准层t3x3的振幅值aps2seis(1)为0.0033；从图17中的d所示过蓬莱111井的实际转换波地震剖面获得的第一个标准层j1d的振幅值aps1seis(2)为0.008和第二个标准层t3x3的振幅值aps2seis(2)为0.005；从图18中的f所示过蓬莱9井的实际转换波地震剖面获得的第一个标准层j1d的振幅值aps1seis(3)为0.005和第二个标准层t3x3的振幅值aps2seis(3)为0.002。

s7：求取反射纵及反射转换横波振幅匹配因子。

依据振幅匹配误差函数关系，使用最小二乘法分别求取反射纵波、反射转换横波振幅匹配因子k(1)＝2105781；k(2)为2902656。

具体地：

s(xj)＝kj,j＝0,1,...,m

j[k(h)]代表所选取的第h个标准层上所有井计算的振幅匹配因子总误差；

k(h)为第h个标准层上计算出的pp波、ps波振幅匹配因子；

n为工区总的井位个数；

i为井号；

apphwell(i)代表第h个标准层在第i口井pp波合成地震记录上标准层振幅值；

apshwell(i)代表第h个标准层在第i口井ps波合成地震记录上标准层振幅值；

apphseis(i)代表第h个标准层在第i口井井位处pp波实际地震记录标准层振幅值；

apshseis(i)代表第h个标准层在第i口井井位处ps波实际地震记录标准层振幅值；

s8：生成时变振幅匹配因子。

将2105781和2902656作为插值区间采用通用的三次样条插值算法生成时变的振幅匹配因子kj；具体地：

，sgn为符号函数；

s(xj)为三次样条插值函数；

j为插值个数，个数范围为[0，m]；

x取值范围为k(1)＝x0<x1<…<xm＝k(h)；k(1)为第一个标准层计算的振幅匹配因子；k(h)为第h个标准层计算的振幅匹配因子。

s9：将反射转换横波的振幅与反射纵波的振幅进行匹配。

将实际转换横波地震数据每个ccp道(共转换点道)乘以时变振幅因子kj，完成纵波、转换横波振幅匹配工作。

图21是振幅匹配前的反射纵波地震记录；图22是振幅匹配前的反射转换横波地震记录；图23是振幅匹配后的反射纵波地震记录；图24是振幅匹配后的反射转换横波地震记录。可见进行振幅能量匹配后，恢复了反射纵波和反射转换横波正常比例关系。

经实践证明，在四川盆地某地区三维三分量工区运用反射纵波与反射转换横波的振幅匹配方法，对反射纵波、反射转换横波实际地震数据振幅关系进行校正，恢复了正确的纵波和转换横波振幅比例关系，为后续多波储层预测工作中利用振幅信息建立储层段纵波、转换横波振幅强弱关系和定量属性预测奠定了正确的基础。

从上述描述可知，本发明实施例提供的反射纵波与反射转换横波的振幅匹配方法，首先读取反射纵波及反射转换横波地震记录上的多个标准层位的同相轴振幅，并读取该多个标准层位的实际地震剖面的同相轴振幅，以此来生成时变振幅匹配因子，最终将反射转换横波的振幅与反射纵波的振幅进行匹配。本方法保证了匹配后反射转换横波振幅的可靠性，减少了单层振幅匹配时远离标准层振幅匹配结果与理论振幅值误差，提高实际生产工作效率，一次性恢复浅中深层转换波数据体地震振幅，避免了多次标准层的选择，为后续利用振幅信息进行储层预测提供可靠数据体。

综上，本发明实施例提供的方法利用多个标准层反射纵波、反射转换横波合成地震记录振幅比例关系校正实际地震记录纵波和转换横波振幅关系，计算出单层振幅比例因子后，通过三次样条差值计算出整个数据体的振幅匹配因子(整个数据体纵向上采用不同的振幅匹配因子)，考虑了时变和空变。能够避免现有技术中振幅匹配后远离标准层振幅匹配结果与理论振幅值误差，能够一次性恢复浅中深层转换波数据体地震振幅，提高实际生产工作效率。具有的优点如下：采用多个标准层、三次样条差值算法进行时变和空变差值，匹配结果考虑了时变和空变。

基于同一发明构思，本申请实施例还提供了反射纵波与反射转换横波的振幅匹配装置，可以用于实现上述实施例所描述的方法，如下面的实施例所述。由于反射纵波与反射转换横波的振幅匹配装置解决问题的原理与反射纵波与反射转换横波的振幅匹配方法相似，因此反射纵波与反射转换横波的振幅匹配装置的实施可以参见反射纵波与反射转换横波的振幅匹配方法实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的系统较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。

本发明的实施例提供一种能够实现反射纵波与反射转换横波的振幅匹配方法的反射纵波与反射转换横波的振幅匹配装置的具体实施方式，参见图25，反射纵波与反射转换横波的振幅匹配装置具体包括如下内容：

第一读取单元10，用于分别读取反射纵波的合成地震记录及反射转换横波的合成地震记录中井中多个标准层位的第一同相轴振幅；

第二读取单元20，用于分别读取反射纵波的实际地震剖面及反射转换横波的实际地震剖面中井中所述多个标准层的第二同相轴振幅；

因子生成单元30，用于根据多个第一同相轴振幅及多个第二同相轴振幅生成时变振幅匹配因子；

匹配单元40，用于根据所述时变振幅匹配因子及所述反射转换横波的实际地震剖面，将所述反射转换横波的振幅与所述反射纵波的振幅进行匹配。

一实施例中，反射纵波与反射转换横波的振幅匹配装置还包括：预处理单元，用于对地震数据进行预处理，生成所述反射纵波和反射转换横波。

一实施例中，反射纵波与反射转换横波的振幅匹配装置还包括：主频数据生成单元，用于分别对所述反射纵波目的时窗内的反射纵波及反射转换横波目的时窗内的反射转换横波进行频谱分析，生成所述反射纵波的主频数据及所述反射转换横波的主频数据。

一实施例中，反射纵波与反射转换横波的振幅匹配装置还包括：子波生成单元，用于根据所述反射纵波的主频数据、所述反射转换横波的主频数据分别生成所述反射纵波的地震子波及所述反射转换横波的地震子波。

一实施例中，反射纵波与反射转换横波的振幅匹配装置还包括反射系数生成单元，用于根据所述井的测井数据计算所述反射纵波的反射系数及所述反射转换横波的反射系数。

一实施例中，所述测井数据包括纵波速度、横波速度及测井密度。

一实施例中，反射纵波与反射转换横波的振幅匹配装置还包括：

纵波合成生成单元，用于将所述反射纵波的反射系数与所述反射纵波的地震子波进行褶积，生成所述反射纵波的合成地震记录；

横波合成生成单元，用于将所述反射转换横波的反射系数与反射转换横波的地震子波进行褶积，生成所述反射转换横波的合成地震记录。

一实施例中，因子生成单元包括：

标准因子生成模块，用于利用最小二乘法，根据多个第一同相轴振幅及多个第二同相轴振幅生成所述井的多个标准层位的标准振幅匹配因子；

时变因子生成模块，用于利用三次样条差值算法，根据所述井的多个标准层位的振幅匹配因子生成所述时变振幅匹配因子。

一实施例中，匹配单元具体用于将所述反射转换横波的实际地震剖面的每个共转换点道集乘以所述时变振幅匹配因子，以将所述反射转换横波的振幅与所述反射纵波的振幅进行匹配。

从上述描述可知，本发明实施例提供的反射纵波与反射转换横波的振幅匹配装置，首先读取反射纵波及反射转换横波地震记录上的多个标准层位的同相轴振幅，并读取该多个标准层位的实际地震剖面的同相轴振幅，以此来生成时变振幅匹配因子，最终将反射转换横波的振幅与反射纵波的振幅进行匹配。本方法保证了匹配后反射转换横波振幅的可靠性，减少了单层振幅匹配时远离标准层振幅匹配结果与理论振幅值误差，提高实际生产工作效率，一次性恢复浅中深层转换波数据体地震振幅，避免了多次标准层的选择，为后续利用振幅信息进行储层预测提供可靠数据体。

综上，本发明实施例提供的装置利用多个标准层反射纵波、反射转换横波合成地震记录振幅比例关系校正实际地震记录纵波和转换横波振幅关系，计算出单层振幅比例因子后，通过三次样条差值计算出整个数据体的振幅匹配因子(整个数据体纵向上采用不同的振幅匹配因子)，考虑了时变和空变。能够避免现有技术中振幅匹配后远离标准层振幅匹配结果与理论振幅值误差，能够一次性恢复浅中深层转换波数据体地震振幅，提高实际生产工作效率。具有的优点如下：采用多个标准层、三次样条差值算法进行时变和空变差值，匹配结果考虑了时变和空变。

本申请的实施例还提供能够实现上述实施例中的反射纵波与反射转换横波的振幅匹配方法中全部步骤的一种电子设备的具体实施方式，参见图26，电子设备具体包括如下内容：

处理器(processor)1201、存储器(memory)1202、通信接口(communicationsinterface)1203和总线1204；

其中，处理器1201、存储器1202、通信接口1203通过总线1204完成相互间的通信；通信接口1203用于实现服务器端设备、采集设备以及用户端设备等相关设备之间的信息传输。

处理器1201用于调用存储器1202中的计算机程序，处理器执行计算机程序时实现上述实施例中的反射纵波与反射转换横波的振幅匹配方法中的全部步骤，例如，处理器执行计算机程序时实现下述步骤：

步骤100：分别读取反射纵波的合成地震记录及反射转换横波的合成地震记录中井中多个标准层位的第一同相轴振幅。

步骤200：分别读取反射纵波的实际地震剖面及反射转换横波的实际地震剖面中井中所述多个标准层的第二同相轴振幅。

步骤300：根据多个第一同相轴振幅及多个第二同相轴振幅生成时变振幅匹配因子。

步骤400：根据所述时变振幅匹配因子及所述反射转换横波的实际地震剖面，将所述反射转换横波的振幅与所述反射纵波的振幅进行匹配。

从上述描述可知，本申请实施例中的电子设备，首先读取反射纵波及反射转换横波地震记录上的多个标准层位的同相轴振幅，并读取该多个标准层位的实际地震剖面的同相轴振幅，以此来生成时变振幅匹配因子，最终将反射转换横波的振幅与反射纵波的振幅进行匹配。本方法保证了匹配后反射转换横波振幅的可靠性，减少了单层振幅匹配时远离标准层振幅匹配结果与理论振幅值误差，提高实际生产工作效率，一次性恢复浅中深层转换波数据体地震振幅，避免了多次标准层的选择，为后续利用振幅信息进行储层预测提供可靠数据体。

综上，本发明实施例提供的电子设备利用多个标准层反射纵波、反射转换横波合成地震记录振幅比例关系校正实际地震记录纵波和转换横波振幅关系，计算出单层振幅比例因子后，通过三次样条差值计算出整个数据体的振幅匹配因子(整个数据体纵向上采用不同的振幅匹配因子)，考虑了时变和空变。能够避免现有技术中振幅匹配后远离标准层振幅匹配结果与理论振幅值误差，能够一次性恢复浅中深层转换波数据体地震振幅，提高实际生产工作效率。具有的优点如下：采用多个标准层、三次样条差值算法进行时变和空变差值，匹配结果考虑了时变和空变。

本申请的实施例还提供能够实现上述实施例中的反射纵波与反射转换横波的振幅匹配方法中全部步骤的一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的反射纵波与反射转换横波的振幅匹配方法的全部步骤，例如，处理器执行计算机程序时实现下述步骤：

步骤100：分别读取反射纵波的合成地震记录及反射转换横波的合成地震记录中井中多个标准层位的第一同相轴振幅。

步骤200：分别读取反射纵波的实际地震剖面及反射转换横波的实际地震剖面中井中所述多个标准层的第二同相轴振幅。

步骤300：根据多个第一同相轴振幅及多个第二同相轴振幅生成时变振幅匹配因子。

步骤400：根据所述时变振幅匹配因子及所述反射转换横波的实际地震剖面，将所述反射转换横波的振幅与所述反射纵波的振幅进行匹配。

从上述描述可知，本申请实施例中的计算机可读存储介质，首先读取反射纵波及反射转换横波地震记录上的多个标准层位的同相轴振幅，并读取该多个标准层位的实际地震剖面的同相轴振幅，以此来生成时变振幅匹配因子，最终将反射转换横波的振幅与反射纵波的振幅进行匹配。本方法保证了匹配后反射转换横波振幅的可靠性，减少了单层振幅匹配时远离标准层振幅匹配结果与理论振幅值误差，提高实际生产工作效率，一次性恢复浅中深层转换波数据体地震振幅，避免了多次标准层的选择，为后续利用振幅信息进行储层预测提供可靠数据体。

综上，本发明实施例提供的计算机可读存储介质利用多个标准层反射纵波、反射转换横波合成地震记录振幅比例关系校正实际地震记录纵波和转换横波振幅关系，计算出单层振幅比例因子后，通过三次样条差值计算出整个数据体的振幅匹配因子(整个数据体纵向上采用不同的振幅匹配因子)，考虑了时变和空变。能够避免现有技术中振幅匹配后远离标准层振幅匹配结果与理论振幅值误差，能够一次性恢复浅中深层转换波数据体地震振幅，提高实际生产工作效率。具有的优点如下：采用多个标准层、三次样条差值算法进行时变和空变差值，匹配结果考虑了时变和空变。

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

上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。

虽然本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。

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

为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本说明书实施例时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。

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

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

在一个典型的配置中，计算设备包括一个或多个处理器(cpu)、输入/输出接口、网络接口和内存。

内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(ram)和/或非易失性内存等形式，如只读存储器(rom)或闪存(flashram)。内存是计算机可读介质的示例。

计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(pram)、静态随机存取存储器(sram)、动态随机存取存储器(dram)、其他类型的随机存取存储器(ram)、只读存储器(rom)、电可擦除可编程只读存储器(eeprom)、快闪记忆体或其他内存技术、只读光盘只读存储器(cd-rom)、数字多功能光盘(dvd)或其他光学存储、磁盒式磁带，磁带磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitorymedia)，如调制的数据信号和载波。

本领域技术人员应明白，本说明书的实施例可提供为方法、系统或计算机程序产品。因此，本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本说明书实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。

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

本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书实施例的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。

以上所述仅为本说明书实施例的实施例而已，并不用于限制本说明书实施例。对于本领域技术人员来说，本说明书实施例可以有各种更改和变化。凡在本说明书实施例的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本说明书实施例的权利要求范围之内。