基于复杂地形的重力场数值模拟方法、装置和计算机设备

技术领域

本申请涉及计算机数值模拟技术领域，特别是涉及一种基于复杂地形的重力场数值模拟方法、装置、计算机设备和存储介质。

背景技术

现阶段，由于重力观测数据的增加，重力场正演所需的内存和计算的急剧增长，对于普通计算机难以实现，已经成为三维密度快速成像的难点问题。

目前，空间域方法以其准确性仍然占据主导地位，但是随着标准FFT扩边法和NUFFT在重力异常正演数值模拟中的应用，频率域方法以其简洁性、准确性和高校性成为处理大规模正演的首选方法。文献(Chai Y , Hinze W J . Gravity inversion of aninterface above which the density contrast varies exponentially with depth.Geophysics, 1988, 53(6):837-845.)采用偏移抽样方法有效提高了傅里叶变换的精度，使频率域方法逐渐应用到复杂模型正演模拟中。文献(Tontini, F.C., L. Cocchi, C.Carmisciano. Rapid 3-D forward model of potential fields with application tothe Palinuro Seamount magnetic anomaly (southern Tyrrhenian Sea, Italy).Journal of Geophysical Research, 2009. 114.)给出了任意密度分布情况下的重力异常频率域表达式。

现有技术中，大尺度大规模的航空重力测量已经兴起，但实际航空测量中地形并非是一个平面，而是复杂的起伏地形，传统空间域方法和频率域方法难以满足海量高分辨率数据精细快速反演要求，存在适应性不佳的问题，因此研究一种基于复杂地形的三维重力异常高效高精度数值模拟方法具有现实意义。

发明内容

基于此，有必要针对上述技术问题，提供一种能够适用于复杂地形重力场数值模拟的基于复杂地形的重力场数值模拟方法、装置、计算机设备和存储介质。

一种基于复杂地形的重力场数值模拟方法，所述方法包括：

根据复杂地形的地形信息和勘探目标的目标信息，构建包含所述复杂地形和所述勘探目标的三维长方体模型，将所述三维长方体模型在x轴方向、y轴方向和z轴方向进行网格剖分，得到多个小长方体，每个所述小长方体在x轴方向和y轴方向的尺寸相同，并根据所述地形信息和所述目标信息设定所述小长方体的空间域密度值；

根据所述三维长方体模型网格剖分的网格参数和预设的高斯参数，得到偏移波数；

根据所述空间域密度值和所述网格参数得到所述三维长方体模型的空间域重力位积分表达式，对所述空间域重力位积分表达式沿水平方向进行二维高斯傅里叶变换得到所述三维长方体模型的频率域重力位积分表达式；

根据所述地形信息，将所述三维长方体模型分为下部分模型和上部分模型，根据所述偏移波数分别计算所述下部分模型的第一垂向一维积分和所述上部分模型的第二垂向一维积分；所述下部分模型是复杂地形最低点以下部分的长方体模型，所述上部分是所述最低点以上部分的长方体模型；

根据所述频率域重力位积分表达式、预设的频率域位场之间的关系信息、所述第一垂向一维积分和所述第二垂向一维积分，得到频率域重力异常表达式，根据所述频率域重力异常表达式计算多个观测面的频率域重力异常场；

对所述频率域重力异常场沿水平方向进行二维高斯傅里叶反变换，得到多个观测面的空间域重力异常场；所述观测面是z轴方向取定值时对应的水平面；

根据所述多个观测面的空间域重力异常场，通过样条插值算法得到所述复杂地形的起伏表面的空间域重力异常场的数值模拟结果。

在其中一个实施例中，还包括：根据所述地形信息和所述目标信息设定所述小长方体的空间域密度值；所述小长方体的密度为常值，不同小长方体的密度值不同，将空气部分小长方体的密度赋为零。

在其中一个实施例中，还包括：根据所述三维长方体模型网格剖分的网格参数和预设的高斯参数，得到偏移波数为：

其中，

当

在其中一个实施例中，还包括：根据所述空间域密度值和所述网格参数得到所述三维长方体模型的空间域重力位积分表达式为：

其中，

对所述空间域重力位积分表达式沿水平方向进行二维傅里叶变换得到所述三维长方体模型的频率域重力位积分表达式为：

其中，

在其中一个实施例中，还包括：根据所述地形信息，将所述三维长方体模型分为下部分模型和上部分模型，根据所述偏移波数分别计算所述下部分模型的第一垂向一维积分和所述上部分模型的第二垂向一维积分为：

其中，

在其中一个实施例中，还包括：根据所述频率域重力位积分表达式、预设的频率域位场之间的关系信息、所述第一垂向一维积分和所述第二垂向一维积分，得到频率域重力异常表达式为：

其中，

在其中一个实施例中，还包括：将不同的偏移波数对应的所述积分系数

计算多个观测面的频率域重力异常场时直接调用存储的值。

一种基于复杂地形的重力场数值模拟装置，所述装置包括：

模型建立模块，用于根据复杂地形的地形信息和勘探目标的目标信息，构建包含所述复杂地形和所述勘探目标的三维长方体模型，将所述三维长方体模型在x轴方向、y轴方向和z轴方向进行网格剖分，得到多个小长方体，每个所述小长方体在x轴方向和y轴方向的尺寸相同，并根据所述地形信息和所述目标信息设定所述小长方体的空间域密度值；

偏移波数计算模块，用于根据所述三维长方体模型网格剖分的网格参数和预设的高斯参数，得到偏移波数；

频率域重力位积分表达式获取模块，用于根据所述空间域密度值和所述网格参数得到所述三维长方体模型的空间域重力位积分表达式，对所述空间域重力位积分表达式沿水平方向进行二维傅里叶变换得到所述三维长方体模型的频率域重力位积分表达式；

垂向一维积分获取模块，用于根据所述地形信息，将所述三维长方体模型分为下部分模型和上部分模型，根据所述偏移波数分别计算所述下部分模型的第一垂向一维积分和所述上部分模型的第二垂向一维积分；所述下部分模型是复杂地形最低点以下部分的长方体模型，所述上部分是所述最低点以上部分的长方体模型；

频率域重力异常场获取模块，用于根据所述频率域重力位积分表达式、预设的频率域位场之间的关系信息、所述第一垂向一维积分和所述第二垂向一维积分，得到频率域重力异常表达式，根据所述频率域重力异常表达式计算多个观测面的频率域重力异常场；

空间域重力异常场获取模块，用于对所述频率域重力异常场沿水平方向进行二维傅里叶反变换，得到多个观测面的空间域重力异常场；所述观测面是z轴方向取特定值时对应的水平面；根据所述多个观测面的空间域重力异常场，通过样条插值算法得到所述复杂地形的起伏表面的空间域重力异常场。

一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现以下步骤：

根据复杂地形的地形信息和勘探目标的目标信息，构建包含所述复杂地形和所述勘探目标的三维长方体模型，将所述三维长方体模型在x轴方向、y轴方向和z轴方向进行网格剖分，得到多个小长方体，每个所述小长方体在x轴方向和y轴方向的尺寸相同，并根据所述地形信息和所述目标信息设定所述小长方体的空间域密度值；

根据所述三维长方体模型网格剖分的网格参数和预设的高斯参数，得到偏移波数；

根据所述空间域密度值和所述网格参数得到所述三维长方体模型的空间域重力位积分表达式，对所述空间域重力位积分表达式沿水平方向进行二维高斯傅里叶变换得到所述三维长方体模型的频率域重力位积分表达式；

根据所述地形信息，将所述三维长方体模型分为下部分模型和上部分模型，根据所述偏移波数分别计算所述下部分模型的第一垂向一维积分和所述上部分模型的第二垂向一维积分；所述下部分模型是复杂地形最低点以下部分的长方体模型，所述上部分是所述最低点以上部分的长方体模型；

根据所述频率域重力位积分表达式、预设的频率域位场之间的关系信息、所述第一垂向一维积分和所述第二垂向一维积分，得到频率域重力异常表达式，根据所述频率域重力异常表达式计算多个观测面的频率域重力异常场；

对所述频率域重力异常场沿水平方向进行二维高斯傅里叶反变换，得到多个观测面的空间域重力异常场；所述观测面是z轴方向取定值时对应的水平面；

根据所述多个观测面的空间域重力异常场，通过样条插值算法得到所述复杂地形的起伏表面的空间域重力异常场的数值模拟结果。

一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：

根据复杂地形的地形信息和勘探目标的目标信息，构建包含所述复杂地形和所述勘探目标的三维长方体模型，将所述三维长方体模型在x轴方向、y轴方向和z轴方向进行网格剖分，得到多个小长方体，每个所述小长方体在x轴方向和y轴方向的尺寸相同，并根据所述地形信息和所述目标信息设定所述小长方体的空间域密度值；

根据所述三维长方体模型网格剖分的网格参数和预设的高斯参数，得到偏移波数；

根据所述空间域密度值和所述网格参数得到所述三维长方体模型的空间域重力位积分表达式，对所述空间域重力位积分表达式沿水平方向进行二维高斯傅里叶变换得到所述三维长方体模型的频率域重力位积分表达式；

根据所述地形信息，将所述三维长方体模型分为下部分模型和上部分模型，根据所述偏移波数分别计算所述下部分模型的第一垂向一维积分和所述上部分模型的第二垂向一维积分；所述下部分模型是复杂地形最低点以下部分的长方体模型，所述上部分是所述最低点以上部分的长方体模型；

根据所述频率域重力位积分表达式、预设的频率域位场之间的关系信息、所述第一垂向一维积分和所述第二垂向一维积分，得到频率域重力异常表达式，根据所述频率域重力异常表达式计算多个观测面的频率域重力异常场；

对所述频率域重力异常场沿水平方向进行二维高斯傅里叶反变换，得到多个观测面的空间域重力异常场；所述观测面是z轴方向取定值时对应的水平面；

根据所述多个观测面的空间域重力异常场，通过样条插值算法得到所述复杂地形的起伏表面的空间域重力异常场的数值模拟结果。

上述基于复杂地形的重力场数值模拟方法、装置、计算机设备和存储介质，通过构建包含复杂地形和勘探目标的三维长方体模型，将三维长方体模型在x、y、z轴方向进行网格剖分，得到多个小长方体，根据地形信息和目标信息设定小长方体的空间域密度值；根据三维长方体模型网格剖分的网格参数和预设的高斯参数，得到偏移波数；根据空间域密度值和网格参数得到三维长方体模型的空间域重力位积分表达式，对空间域重力位积分表达式沿水平方向进行二维傅里叶变换得到三维长方体模型的频率域重力位积分表达式；根据地形信息，将三维长方体模型分为下部分模型和上部分模型，根据偏移波数分别计算下部分模型的第一垂向一维积分和上部分模型的第二垂向一维积分，得到频率域重力异常表达式，根据频率域重力异常表达式计算多个观测面的频率域重力异常场；对频率域重力异常场沿水平方向进行二维高斯傅里叶反变换，得到多个观测面的空间域重力异常场，通过样条插值算法得到复杂地形的起伏表面的空间域重力异常场。本发明将研究区域剖分为若干个小长方体，根据每个小长方体密度变化刻画复杂地形；该方法垂直方向保留在空间域，网格剖分灵活，适合模拟复杂地形及复杂地质体；并且多个偏移波数相互独立，便于并行计算，从而大大提高了重力异常计算效率、减少了计算所需的内存；根据多个观测平面上的重力异常采用样条插值的方法计算任意复杂地形上的重力异常，实现基于复杂地形的快速重力异常正演。

附图说明

图1为一个实施例中基于复杂地形的重力场数值模拟方法的流程示意图；

图2为另一个实施例中基于复杂地形的重力场数值模拟方法的流程示意图；

图3为一个具体实施例中复杂地形模型示意图；

图4为一个实施例中空间域重力异常场

图5为一个实施例中采用本发明所述方法和采用传统累加方法的算法计算时间随复杂地形观测面个数变化的趋势；

图6为一个实施例中基于复杂地形的重力场数值模拟装置的结构框图；

图7为一个实施例中计算机设备的内部结构图。

具体实施方式

为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本申请，并不用于限定本申请。

本申请提供的基于复杂地形的重力场数值模拟方法，可以应用于如下应用环境中。其中，通过终端执行一种基于复杂地形的重力场数值模拟方法。给定研究区域、三维地质体密度分布及复杂地形模型数据；根据研究区域、三维地质体密度分布和三维频率域重力异常正演方法，计算三维地质体在复杂地形表面产生的重力异常；当所述方法计算的重力异常和采用重力仪测量得到的三维地质体产生的重力异常相同，将所述方法用于勘探所述密度分布异常体。其中，终端可以但不限于是各种个人计算机、笔记本电脑、平板电脑和便携式设备。

在一个实施例中，如图1所示，提供了一种基于复杂地形的重力场数值模拟方法，包括以下步骤：

步骤102，根据复杂地形的地形信息和勘探目标的目标信息，构建包含复杂地形和勘探目标的三维长方体模型，将三维长方体模型在x轴方向、y轴方向和z轴方向进行网格剖分，得到多个小长方体，每个小长方体在x轴方向和y轴方向的尺寸相同，并根据地形信息和目标信息设定小长方体的空间域密度值。

根据勘探目标的形状、大小及密度分布及复杂地形构建包含所述勘探目标及复杂地形的三维长方体模型；将三维长方体模型进行网格剖分，

步骤104，根据三维长方体模型网格剖分的网格参数和预设的高斯参数，得到偏移波数。

给定高斯参数，高斯参数包括高斯点和对应的系数，利用高斯参数计算偏移波数时可以采用任意个高斯点。实际应用中，当算法精度不够时，可以通过增加高斯点个数来提高精度，适用性强，尤其是任意密度分布任意截面形状分布的复杂介质模型。

步骤106，根据空间域密度值和网格参数得到三维长方体模型的空间域重力位积分表达式，对空间域重力位积分表达式沿水平方向进行二维高斯傅里叶变换得到三维长方体模型的频率域重力位积分表达式。

高斯傅里叶变换方法减小了边界截断效应问题，可以提高计算精度。

步骤108，根据地形信息，将三维长方体模型分为下部分模型和上部分模型，根据偏移波数分别计算下部分模型的第一垂向一维积分和上部分模型的第二垂向一维积分；下部分模型是复杂地形最低点以下部分的长方体模型，上部分是最低点以上部分的长方体模型。

下部分模型是复杂地形最低点以下部分的长方体模型，计算起伏地形最低点与最高点之间不同高度上的重力异常时，下部分模型是复杂地形最低点以下部分的长方体模型部分积分需要不断重复计算，造成计算时间的叠加，而将起伏地形以下部分先计算存储，需要时直接调用，节约计算成本和内存空间；上部分是最低点以上部分的长方体模型，分为该观测点以上和以下两部分，依次计算存储，计算最低点与最高点异常时直接调用，提高计算效率。

步骤110，根据频率域重力位积分表达式、预设的频率域位场之间的关系信息、第一垂向一维积分和第二垂向一维积分，得到频率域重力异常表达式，根据频率域重力异常表达式计算多个观测面的频率域重力异常场。

观测面是z轴方向取定值时对应的水平面，设定多个z值计算多个观测面的频率域重力异常场，相当于观测重力异常场在z轴方向的变化，以用于之后计算复杂地形起伏表面的重力异常。

步骤112，对频率域重力异常场沿水平方向进行二维高斯傅里叶反变换，得到多个观测面的空间域重力异常场。

步骤114，根据多个观测面的空间域重力异常场，通过样条插值算法得到复杂地形的起伏表面的空间域重力异常场的数值模拟结果。

三次样条插值函数为：

上式中z表示观测面的z轴坐标值，系数

上述基于复杂地形的重力场数值模拟方法中，通过构建包含复杂地形和勘探目标的三维长方体模型，将三维长方体模型在x、y、z轴方向进行网格剖分，得到多个小长方体，根据地形信息和目标信息设定小长方体的空间域密度值；根据三维长方体模型网格剖分的网格参数和预设的高斯参数，得到偏移波数；根据空间域密度值和网格参数得到三维长方体模型的空间域重力位积分表达式，对空间域重力位积分表达式沿水平方向进行二维傅里叶变换得到三维长方体模型的频率域重力位积分表达式；根据地形信息，将三维长方体模型分为下部分模型和上部分模型，根据偏移波数分别计算下部分模型的第一垂向一维积分和上部分模型的第二垂向一维积分，得到频率域重力异常表达式，根据频率域重力异常表达式计算多个观测面的频率域重力异常场；对频率域重力异常场沿水平方向进行二维傅里叶反变换，得到多个观测面的空间域重力异常场，通过样条插值算法得到复杂地形的起伏表面的空间域重力异常场。本发明将研究区域剖分为若干个小长方体，根据每个小长方体密度变化刻画复杂地形；该方法垂直方向保留在空间域，网格剖分灵活，适合模拟复杂地形及复杂地质体；并且多个偏移波数相互独立，便于并行计算，从而大大提高了重力异常计算效率、减少了计算所需的内存；根据多个观测平面上的重力异常采用样条插值的方法计算任意复杂地形上的重力异常，实现基于复杂地形的快速重力异常正演。

在其中一个实施例中，还包括：根据地形信息和目标信息设定小长方体的空间域密度值；小长方体的密度为常值，不同小长方体的密度值不同，将空气部分小长方体的密度赋为零。

在其中一个实施例中，还包括：根据三维长方体模型网格剖分的网格参数和预设的高斯参数，得到偏移波数为：

其中，

当

在其中一个实施例中，还包括：根据空间域密度值和网格参数得到三维长方体模型的空间域重力位积分表达式为：

其中，

对空间域重力位积分表达式沿水平方向进行二维傅里叶变换得到三维长方体模型的频率域重力位积分表达式为：

其中，

在其中一个实施例中，还包括：根据地形信息，将三维长方体模型分为下部分模型和上部分模型，根据偏移波数分别计算下部分模型的第一垂向一维积分和上部分模型的第二垂向一维积分为：

其中，

在其中一个实施例中，还包括：根据频率域重力位积分表达式、预设的频率域位场之间的关系信息、第一垂向一维积分和第二垂向一维积分，得到频率域重力异常表达式为：

其中，

在其中一个实施例中，还包括：将不同的偏移波数对应的积分系数

根据下部分模型的第一垂向一维积分和上部分模型的第二垂向一维积分的计算公式，将公共部分

应该理解的是，虽然图1的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图1中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。

在一个具体实施例中，如图2所示，提供了一种基于复杂地形的重力场数值模拟方法，根据勘探目标体的大小和密度分布构建包含勘探目标的长方体模型，将大的长方体模型剖分成若干个小长方体模型、给定模型参数，并对每个小长方体进行密度赋值，给定高斯点个数并及其相应的高斯点和高斯系数，根据模型参数、高斯点数计算二维高斯傅里叶变换偏移波数，对空间域重力位三维卷积公式沿水平方向做二维高斯傅里叶变换，采用二维高斯傅里叶变换计算频率域密度，根据频率域求导的特性，直接得到频率域重力异常场公式，根据分块叠加计算复杂地形的多个水平观测面上重力异常场，对频率域重力异常沿水平方向进行二维高斯傅里叶反变换，对多个观测面上重力异常采用样条插值算法计算起伏地形上重力异常场。

在一个具体实施例中，设计了如图3所示的复杂地形模型如下：

模拟区域范围为：

本实施例中的三维重力异常计算方法利用Fortran语言编程实现，运行程序所用的个人笔记本电脑配置为：CPU-Inter Core i7-8700，主频为3.2GHz，运行内存为8.00GB。图4为本发明所述方法计算的复杂地形上的重力异常场三分量解析解、数值解以及相对误差等值线图，从图4相对误差图中可以看出最大相对误差为

在一个实施例中，如图6所示，提供了一种基于复杂地形的重力场数值模拟装置，包括：模型建立模块602、偏移波数计算模块604、频率域重力位积分表达式获取模块606、垂向一维积分获取模块608、频率域重力异常场获取模块610和空间域重力异常场获取模块612，其中：

模型建立模块602，用于根据复杂地形的地形信息和勘探目标的目标信息，构建包含复杂地形和勘探目标的三维长方体模型，将三维长方体模型在x轴方向、y轴方向和z轴方向进行网格剖分，得到多个小长方体，每个小长方体在x轴方向和y轴方向的尺寸相同，并根据地形信息和目标信息设定小长方体的空间域密度值；

偏移波数计算模块604，用于根据三维长方体模型网格剖分的网格参数和预设的高斯参数，得到偏移波数；

频率域重力位积分表达式获取模块606，用于根据空间域密度值和网格参数得到三维长方体模型的空间域重力位积分表达式，对空间域重力位积分表达式沿水平方向进行二维傅里叶变换得到三维长方体模型的频率域重力位积分表达式；

垂向一维积分获取模块608，用于根据地形信息，将三维长方体模型分为下部分模型和上部分模型，根据偏移波数分别计算下部分模型的第一垂向一维积分和上部分模型的第二垂向一维积分；下部分模型是复杂地形最低点以下部分的长方体模型，上部分是最低点以上部分的长方体模型；

频率域重力异常场获取模块610，用于根据频率域重力位积分表达式、预设的频率域位场之间的关系信息、第一垂向一维积分和第二垂向一维积分，得到频率域重力异常表达式，根据频率域重力异常表达式计算多个观测面的频率域重力异常场；

空间域重力异常场获取模块612，用于对频率域重力异常场沿水平方向进行二维傅里叶反变换，得到多个观测面的空间域重力异常场；观测面是z轴方向取特定值时对应的水平面；根据多个观测面的空间域重力异常场，通过样条插值算法得到复杂地形的起伏表面的空间域重力异常场。

偏移波数计算模块604还用于根据三维长方体模型网格剖分的网格参数和预设的高斯参数，得到偏移波数为：

其中，

当

频率域重力位积分表达式获取模块606还用于根据空间域密度值和网格参数得到三维长方体模型的空间域重力位积分表达式为：

其中，

对空间域重力位积分表达式沿水平方向进行二维傅里叶变换得到三维长方体模型的频率域重力位积分表达式为：

其中，

垂向一维积分获取模块608还用于根据地形信息，将三维长方体模型分为下部分模型和上部分模型，根据偏移波数分别计算下部分模型的第一垂向一维积分和上部分模型的第二垂向一维积分为：

其中，

频率域重力异常场获取模块610还用于根据频率域重力位积分表达式、预设的频率域位场之间的关系信息、第一垂向一维积分和第二垂向一维积分，得到频率域重力异常表达式为：

其中，

空间域重力异常场获取模块612还用于在根据频率域重力异常表达式计算多个观测面的频率域重力异常场之前，将不同的偏移波数对应的积分系数

关于基于复杂地形的重力场数值模拟装置的具体限定可以参见上文中对于基于复杂地形的重力场数值模拟方法的限定，在此不再赘述。上述基于复杂地形的重力场数值模拟装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。

在一个实施例中，提供了一种计算机设备，该计算机设备可以是终端，其内部结构图可以如图7所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口、显示屏和输入装置。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种基于复杂地形的重力场数值模拟方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置可以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。

本领域技术人员可以理解，图7中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。

在一个实施例中，提供了一种计算机设备，包括存储器和处理器，该存储器存储有计算机程序，该处理器执行计算机程序时实现上述方法实施例中的步骤。

在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现上述方法实施例中的步骤。

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器（ROM）、可编程ROM（PROM）、电可编程ROM（EPROM）、电可擦除可编程ROM（EEPROM）或闪存。易失性存储器可包括随机存取存储器（RAM）或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM（SRAM）、动态RAM（DRAM）、同步DRAM（SDRAM）、双数据率SDRAM（DDRSDRAM）、增强型SDRAM（ESDRAM）、同步链路（Synchlink） DRAM（SLDRAM）、存储器总线（Rambus）直接RAM（RDRAM）、直接存储器总线动态RAM（DRDRAM）、以及存储器总线动态RAM（RDRAM）等。

以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。

以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所附权利要求为准。