大气颗粒切割器性能的模拟方法和装置与流程

本发明涉及大气颗粒切割器技术领域，尤其是涉及一种大气颗粒切割器性能的模拟方法和装置。

背景技术：

大气颗粒物是形成大气污染的重要成分，并且已经成为诸多城市的首要大气污染物；大气颗粒物对空气质量的能见度、酸沉降、云和降水、大气的辐射平衡等因素均有重要影响。PM10为空气动力学直径小于10um的颗粒物，在TSP(total suspended particulate，总悬浮颗粒物)中占有比例高达60-80％。为了能精确测量大气中PM10、PM2.5的浓度，可以采用切割器将这些细颗粒物从空气中分离出来，送至检测装置。

为了评价切割器的性能，通常采用实验方法可以对PM10切割器的性能进行研究、评价；但需要较为繁复的准备工作和昂贵的设备资金，且数据结果的精确度较差。

技术实现要素：

有鉴于此，本发明的目的在于提供一种大气颗粒切割器性能的模拟方法和装置，以数值模拟的方式获得切割器的多种场分布，得到切割器颗粒分离性能的模拟结果，使操作便捷且数据精确度较高。

第一方面，本发明实施例提供了一种大气颗粒切割器性能的模拟方法，包括：建立大气颗粒切割器的三维模型；根据预设的气相湍流模型和边界条件，对三维模型进行六面体网格划分；根据切割器的额定工作流量，确定压力场分布、速度场分布和速度矢量场分布；根据预设的离散相模型，确定颗粒粒径和捕集效率的对应关系；根据对应关系，确定切割器的颗粒分离性能。

结合第一方面，本发明实施例提供了第一方面的第一种可能的实施方式，其中，上述根据预设的气相湍流模型和边界条件，对三维模型进行六面体网格划分的步骤，包括：采用RNGκ-ε双方程湍流模型和边界条件，对三维模型进行六面体网格划分；其中，RNGκ-ε双方程湍流模型中，控制方程的压力项采用PRESTO差分格式求解，除压力项之外的剩余项采用QUICK差分格式求解；压力与速度耦合采用SIMPLEC算法求解，颗粒项选用DPM。

结合第一方面，本发明实施例提供了第一方面的第二种可能的实施方式，其中，上述根据预设的离散相模型，确定颗粒粒径和捕集效率的对应关系的步骤，包括：对每种粒径的颗粒进行多次采样，对采样数据进行分析，获得对应的初始平均捕集效率和平均捕集效率的相对偏差；将相对偏差小于设定阈值的初始平均捕集效率确定为当前粒径对应的捕集效率。

结合第一方面的第二种可能的实施方式，本发明实施例提供了第一方面的第三种可能的实施方式，其中，上述根据对应关系，确定切割器的颗粒分离性能的步骤，包括：将多种粒径对应的粒径和捕集效率的对应关系输入至回归方程η＝A+BD_a¹+CD_a²+DD_a³，计算获得系数A、B、C和D；计算η＝16％、50％和84％对应的颗粒空气动力学当量直径D_a16、D_a50和D_a84；计算空气动力学当量直径的几何标准差；根据D_a50和几何标准差，计算切割器的分离效率。

结合第一方面，本发明实施例提供了第一方面的第四种可能的实施方式，其中，上述方法还包括：计算获得多种气体流速下的，颗粒运动轨迹和逃逸率曲线；根据颗粒运动轨迹和逃逸率曲线，对切割器内的压降进行定量分析；根据分析结果，优化切割器的分离效率。

结合第一方面的第四种可能的实施方式，本发明实施例提供了第一方面的第五种可能的实施方式，其中，上述根据分析结果，优化切割器的分离效率的步骤，包括：根据环境空气中颗粒物的粒径分布数据，模拟粒径分布函数；根据季节、气相条件、颗粒物主要成因调整切割器相关系数；根据逃逸率曲线，调整切割器的摆布、采气流速，以优化切割器的分离效率和采集样品质量。

第二方面，本发明实施例提供了一种大气颗粒切割器性能的模拟装置，包括：模型建立模块，用于建立大气颗粒切割器的三维模型；网格划分模块，用于根据预设的气相湍流模型和边界条件，对三维模型进行六面体网格划分；分布场确定模块，用于根据切割器的额定工作流量，确定压力场分布、速度场分布和速度矢量场分布；对应关系确定模块，用于根据预设的离散相模型，确定颗粒粒径和捕集效率的对应关系；性能确定模块，用于根据对应关系，确定切割器的颗粒分离性能。

结合第二方面，本发明实施例提供了第二方面的第一种可能的实施方式，其中，上述网格划分模块，还用于：采用RNGκ-ε双方程湍流模型和边界条件，对三维模型进行六面体网格划分；其中，RNGκ-ε双方程湍流模型中，控制方程的压力项采用PRESTO差分格式求解，除压力项之外的剩余项采用QUICK差分格式求解；压力与速度耦合采用SIMPLEC算法求解，颗粒项选用DPM。

结合第二方面，本发明实施例提供了第二方面的第二种可能的实施方式，其中，上述对应关系确定模块，还用于：对每种粒径的颗粒进行多次采样，对采样数据进行分析，获得对应的初始平均捕集效率和平均捕集效率的相对偏差；将相对偏差小于设定阈值的初始平均捕集效率确定为当前粒径对应的捕集效率。

结合第二方面的第二种可能的实施方式，本发明实施例提供了第二方面的第三种可能的实施方式，其中，上述性能确定模块，还用于：将多种粒径对应的粒径和捕集效率的对应关系输入至回归方程η＝A+BD_a¹+CD_a²+DD_a³，计算获得系数A、B、C和D；计算η＝16％、50％和84％对应的颗粒空气动力学当量直径D_a16、D_a50和D_a84；计算空气动力学当量直径的几何标准差；根据D_a50和几何标准差，计算切割器的分离效率。

本发明实施例带来了以下有益效果：

本发明实施例提供的一种大气颗粒切割器性能的模拟方法和装置，建立大气颗粒切割器的三维模型，根据预设的气相湍流模型和边界条件，对该三维模型进行六面体网格划分；根据切割器的额定工作流量，确定压力场分布、速度场分布和速度矢量场分布；再根据预设的离散相模型，确定颗粒粒径和捕集效率的对应关系，进而确定切割器的颗粒分离性能。该方式基于流体动力学，以数值模拟的方式获得切割器的多种场分布，在结合粒子运动轨迹规律，得到切割器颗粒分离性能的模拟结果，操作便捷且数据精确度较高，可以帮助进一步优化切割器的分离效果。

本发明的其他特征和优点将在随后的说明书中阐述，或者，部分特征和优点可以从说明书推知或毫无疑义地确定，或者通过实施本发明的上述技术即可得知。

为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施方式，并配合所附附图，作详细说明如下。

附图说明

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

图1为本发明实施例提供的一种PM10切割器的三维结构示意图；

图2为本发明实施例提供的一种PM10切割器的二维结构示意图；

图3为本发明实施例提供的一种大气颗粒切割器性能的模拟方法的流程图；

图4为本发明实施例提供的流量与捕集效率的对应关系曲线图；

图5为本发明实施例提供的参数X＝0mm情况下的压力场分布示意图；

图6为本发明实施例提供的参数Z＝0mm情况下的压力场分布示意图；

图7为本发明实施例提供的参数X＝0mm情况下的速度场分布示意图；

图8为本发明实施例提供的参数Y＝141mm情况下的速度场分布示意图；

图9为本发明实施例提供的参数X＝0mm情况下的速度矢量场分布示意图；

图10为本发明实施例提供的另一种大气颗粒切割器性能的模拟方法的流程图；

图11为本发明实施例提供的颗粒粒径和捕集效率的对应关系曲线图；

图12为本发明实施例提供的一种大气颗粒切割器性能的模拟装置的结构示意图。

具体实施方式

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

考虑到现有的大气颗粒切割器的性能评价方式操作繁琐，且数据精确度较低的问题，本发明实施例提供了一种大气颗粒切割器性能的模拟方法和装置；该技术可以应用于PM10切割器、冲击式切割器等多种切割器的数值模拟、仿真以及性能评价中；该技术可以采用相关的软件或硬件实现，下面通过实施例进行描述。

为便于对本实施例进行理解，首先对本发明实施例所公开的一种PM10切割器进行详细介绍，参见图1所示的一种PM10切割器的三维结构示意图和图2所示的一种PM10切割器的二维结构示意图。该PM10切割器包括入口、冲击板、集尘板、滤膜、滤膜架和连接杆等部分；其中，滤膜采用玻璃纤维新型材料；冲击板和集尘板均采用PL102材料。

通常，切割器的切割效率和能耗受到了入口面积、冲击孔半径、集尘板面积、滤膜架高度和连接杆高度的影响。基于流体动力学技术，模拟出冲击式切割器内部气流的速度矢量场、压力场、速度场等，结合粒子运动轨迹规律，得到准确的数据模拟结果。因此，CFD(computational fluid dynamics，计算流体动力学)技术被广泛应用于数据模拟中。

该PM10切割器是根据冲击原理设计的中流量多级颗粒物切割器；由于PM10(一种可吸入颗粒物)的冲击惯性，使PM10在切割器中能够精确的分离出来。具体的分离原理描述如下：如图2所示，含尘气体由入口进入；含尘气体途径冲击板时，由于冲击板孔径大小的限制对颗粒进行阻碍分离，较大粒径的颗粒由于粒径大沉积在集尘板上，细小粒径的颗粒由于惯性力与压力力作用继续向下运动到达滤膜，细小粒径的颗粒由于滤膜的分离作用沉积于滤膜之上，达到气固分离，分离后的气体经过连接杆流出PM10切割器。

通常，采用实验方法可以对PM10切割器的性能进行研究、评价；但需要较为繁复的准备工作和昂贵的设备资金，且数据结果的精确度较差；基于此，参见图3所示的一种大气颗粒切割器性能的模拟方法的流程图；该方法包括如下步骤：

步骤S302，建立大气颗粒切割器的三维模型；

步骤S304，根据预设的气相湍流模型和边界条件，对三维模型进行六面体网格划分；

该大气颗粒切割器通常为冲击式切割器；在实际实现时，可以采用CFD软件，例如，GAMBIT，模拟内流流体力学，从而建立该切割器的三维模型、对该三维模型进行六面体网格划分并设置边界条件；后续的数值计算步骤可以采用FLUENT软件实现。

在进行网格划分时，通常需要考虑下述三个因素：

(1)设置时间，实际工程中，当涉及到比较复杂的几何形状时，建立非结构或多块网格时十分耗费时间的；对于复杂的几何形状问题，设置网格的时间是使用四面体或六面体单元的结构网格的主要考虑因素。

(2)计算成本，当几何形状复杂或流程的长度尺度的范围比较大的时候，选用六面体网格可以提供更多的优越条件，六面体网格能够允许一个比四面体单元大的多的纵横比，而四面体单元中一个大的纵横比会影响单元的斜率，会妨碍到计算的精确与收敛。

(3)数值耗散，在多维情形中，一个主要的错误来源在于数值耗散，即对一个流动计算的影响近似于增加真实耗散系数的影响。数值耗散起于切断错误，而切断错误是一个表达离散形式的流体流动方程的结果。数值耗散的总数反过来与网格的分解有关。因此，处理数值耗散的一个方法就是改进网格。当流动与网格相吻合一致时，数值耗散减到最小。

参见图4所示的流量与捕集效率的对应关系曲线图；在四面体和六面体网格各占不同比例的情况下，模拟所得的捕集效率基本相同。在0-40L/min阶段，捕集效率会上升到最高水平；在40-110L/min阶段，捕集效率水平稳定在90％左右；在110-200L/min阶段，捕集效率会随着流量升高而降低。将四面体网格的捕集效率与六面体网格的捕集效率进行对比可知在流率相同的情况下，六面体网格对于颗粒物的捕集效率高于四面体网格。同时，对于PM10切割器六面体网格划分更加精密，对区域的覆盖面积更大，使非结构化网格的计算更加精确，提高了运算效率。

综合上述考虑因素，本实施例中对切割器的三维模型进行六面体网格划分。

步骤S306，根据切割器的额定工作流量，确定压力场分布、速度场分布和速度矢量场分布；

参见图5所示的参数X＝0mm情况下的压力场分布示意图和图6所示的参数Z＝0mm情况下的压力场分布示意图；压降是表征冲击式切割器性能的一项重要参数，不仅直接影响到能耗的大小，而且还会影响到操作条件的选择。切割器整体的总压降中80％都来自冲击孔处与滤膜处，是切割器内部气体遇到孔径突然缩小的情况和致密的滤膜来过滤含尘气体，从而引起能量损失，剩下的20％分别是气体在出口处的收缩、进口处的膨胀和与器壁表面的摩擦引起的损失。由图5和图6可知，滤膜造成了巨大的压降，大约达到了500pa的压力差，巨大的压降差能够使得含尘气体中的PM10更好的脱附下来，累积到滤膜上，形成气固分离。

参见图7所示的参数X＝0mm情况下的速度场分布示意图和图8所示的参数Y＝141mm情况下的速度场分布示意图；在冲击式切割器的额定工作流量(100L/min)下，每个入口的含尘气体流入速度均为0.621315m/s，在该标准流入速度的流入情况下，可以得到切割器内各处压力与速度之间的关系，由图5和图6可知，滤膜以上区域压降大约为10pa，并且主要压降区域为冲击孔到集尘板之间，这段区域压力差达到了5pa。因此，含尘气体在冲击孔中能够得到足够大的冲击力，较大的压力差加大了气体流速，使得含尘气体在冲击板上能够达到更加有效的分离效果。

由图7可知，通过X＝0mm轴截面图可得切割器内部整体速度分布，根据图中切割器上的标注，在切割器的四个入口速度均为0.6213m/s,向下到达冲击孔时，速度达到2.164m/s，速度突然增大，形成巨大的冲击力，使得气体能够成功冲击到冲击板。气体从冲击板的缝隙向下流动，速度逐渐减小到0.96m/s。通过滤膜后在流动到连接杆的过程中速度均匀,达到了1.983m/s。在流入连接杆中，速度逐渐增大，从大约2m/s增加到17m/s。

由图8可知，在切割器内，速度分布为中间以及四周的速度最大，达到了1.3m/s以上，冲击孔处速度大约在0.5-1m/s。由于切割器内部速度变化较大，在积尘板与滤膜之间由于四周空隙处速度过大，因此形成了湍流。其余部分速度梯度变化小，从而保证了颗粒物的竖直流动，确保了分离效率。

参见图9所示的参数X＝0mm情况下的速度矢量场分布示意图；从图9的流场中可以看出，对于X＝0mm的平面，图中显示的是内孔，而流入孔中的气流主要为孔上内孔之间的气流。在气流流入孔的三分之一后，气流方向就已经与孔平行，并且垂直于集尘板，这保证了切割时的最佳来流方向。出孔后，两侧的气流向外流出，若此时颗粒随气流流出，则造成无法撞击；而两侧以内的大部分气流则会抵达集尘板，发生撞击现象。发生撞击后的气流，由中央和四周的通道向下流向滤膜。

步骤S308，根据预设的离散相模型，确定颗粒粒径和捕集效率的对应关系；

步骤S310，根据对应关系，确定切割器的颗粒分离性能。

本发明实施例提供的一种大气颗粒切割器性能的模拟方法，建立大气颗粒切割器的三维模型，根据预设的气相湍流模型和边界条件，对该三维模型进行六面体网格划分；根据切割器的额定工作流量，确定压力场分布、速度场分布和速度矢量场分布；再根据预设的离散相模型，确定颗粒粒径和捕集效率的对应关系，进而确定切割器的颗粒分离性能。该方式基于流体动力学，以数值模拟的方式获得切割器的多种场分布，在结合粒子运动轨迹规律，得到切割器颗粒分离性能的模拟结果，操作便捷且数据精确度较高，可以帮助进一步优化切割器的分离效果。

参见图10所示的另一种大气颗粒切割器性能的模拟方法的流程图；该方法在图3中所示方法基础上实现；该方法包括如下步骤：

步骤S1002，建立大气颗粒切割器的三维模型；

步骤S1004，采用RNGκ-ε双方程湍流模型和边界条件，对三维模型进行六面体网格划分；其中，RNGκ-ε双方程湍流模型中，控制方程的压力项采用PRESTO差分格式求解，除压力项之外的剩余项采用QUICK差分格式求解；压力与速度耦合采用SIMPLEC算法求解，颗粒项选用DPM。

上述κ代表湍动能，ε代表耗散率；上述PRESTO差分格式、QUICK差分格式、SIMPLEC算法和DPM均可以在FLUENT软件中进行选择和参数设置。根据一般情况，假定颗粒物密度为1.225g/cm，模型精度为1×10³，冲击式切割器的流场相当复杂，为了验证模型的有效性，可以对流场中三位分布的速度值进行验证。

步骤S1006，根据切割器的额定工作流量，确定压力场分布、速度场分布和速度矢量场分布；

步骤S1008，对每种粒径的颗粒进行多次采样，对采样数据进行分析，获得对应的初始平均捕集效率和平均捕集效率的相对偏差；

步骤S1010，将相对偏差小于设定阈值的初始平均捕集效率确定为当前粒径对应的捕集效率。

步骤S1012，将多种粒径对应的粒径和捕集效率的对应关系输入至回归方程η＝A+BD_a¹+CD_a²+DD_a³，计算获得系数A、B、C和D；

步骤S1014，计算η＝16％、50％和84％对应的颗粒空气动力学当量直径D_a16、D_a50和D_a84；

步骤S1016，计算空气动力学当量直径的几何标准差；根据D_a50和几何标准差，计算切割器的分离效率。

在进行数值模拟时，通常采用七种不同粒径单分散相(σg≤1.1)荧光素铵(C20H15O5N)球形粒子作为实验粒子，对切割器的颗粒捕捉性能进行检测；在额定电压下，待切割器工作稳定后，对每一种粒径的实验粒子，采样三次，经分析后得出切割器对于该粒径的三个捕集效率η_Daik，再分别计算出η_Daik的平均捕集效率(η_Dai)和相对偏差(C_v)。当相对偏差C_v≤15％时，η_Dai即为切割器对于该粒径的平均捕集效率。

将切割器的七种粒子的粒径及所对应的切割器的七个捕集效率采用多项式回归方法给出；该回归方程即η＝A+BD_a¹+CD_a²+DD_a³，求出A、B、C、D四个系数；再由回归方程求出切割器的捕集效率η为16％、50％、84％时，所对应的粒子空气动力学当量直径D_a16、D_a50、D_a84；计算切割粒径(D_a50)和几何标准差(σg)，通过D_a50和σg评价切割器的性能，并判断该性能是否满足预设的检测标准条件。

参见图11所示的颗粒粒径和捕集效率的对应关系曲线图；其中应用本实施例方法获得的模拟数值表明，PM10的捕集效率达到了93％；而切割器的实验数据表示切割器对PM10的捕集效率为50％。由于大气中颗粒污染物中PM10所占份数达到了99％，应用本实施例方法获得的模拟数值，采集的粉尘质量更加接近于实际值。在一些污染较为严重的工业区，重污染天气将抵消上述误差，污染越严重，误差越小，模拟数值数据更加接近真实值。

步骤S1018，计算获得多种气体流速下的，颗粒运动轨迹和逃逸率曲线；根据颗粒运动轨迹和逃逸率曲线，对切割器内的压降进行定量分析；

步骤S1020，根据分析结果，优化切割器的分离效率。

上述步骤S1020，具体可以通过下述方式实现：

(1)根据环境空气中颗粒物的粒径分布数据，模拟粒径分布函数；

(2)根据季节、气相条件、颗粒物主要成因调整切割器相关系数；

(3)根据逃逸率曲线，调整切割器的摆布、采气流速，以优化切割器的分离效率和采集样品质量。

本发明实施例提供的一种大气颗粒切割器性能的模拟方法，基于流体动力学，以数值模拟的方式获得切割器的多种场分布，在结合粒子运动轨迹规律，得到切割器颗粒分离性能的模拟结果，操作便捷且数据精确度较高，可以帮助进一步优化切割器的分离效果。

对应于上述方法实施例，参见图12所示的一种大气颗粒切割器性能的模拟装置的结构示意图；该装置包括：

模型建立模块120，用于建立大气颗粒切割器的三维模型；

网格划分模块121，用于根据预设的气相湍流模型和边界条件，对三维模型进行六面体网格划分；

分布场确定模块122，用于根据切割器的额定工作流量，确定压力场分布、速度场分布和速度矢量场分布；

对应关系确定模块123，用于根据预设的离散相模型，确定颗粒粒径和捕集效率的对应关系；

性能确定模块124，用于根据对应关系，确定切割器的颗粒分离性能。

上述网格划分模块，还用于：采用RNGκ-ε双方程湍流模型和边界条件，对三维模型进行六面体网格划分；其中，RNGκ-ε双方程湍流模型中，控制方程的压力项采用PRESTO差分格式求解，除压力项之外的剩余项采用QUICK差分格式求解；压力与速度耦合采用SIMPLEC算法求解，颗粒项选用DPM。

上述对应关系确定模块，还用于：对每种粒径的颗粒进行多次采样，对采样数据进行分析，获得对应的初始平均捕集效率和平均捕集效率的相对偏差；将相对偏差小于设定阈值的初始平均捕集效率确定为当前粒径对应的捕集效率。

上述性能确定模块，还用于：将多种粒径对应的粒径和捕集效率的对应关系输入至回归方程η＝A+BD_a¹+CD_a²+DD_a³，计算获得系数A、B、C和D；计算η＝16％、50％和84％对应的颗粒空气动力学当量直径D_a16、D_a50和D_a84；计算空气动力学当量直径的几何标准差；根据D_a50和几何标准差，计算切割器的分离效率。

上述装置还包括：计算模块，用于计算获得多种气体流速下的，颗粒运动轨迹和逃逸率曲线；定量分析模块，用于根据颗粒运动轨迹和逃逸率曲线，对切割器内的压降进行定量分析；优化模块，用于根据分析结果，优化切割器的分离效率。

本发明实施例提供的大气颗粒切割器性能的模拟装置，与上述实施例提供的大气颗粒切割器性能的模拟方法具有相同的技术特征，所以也能解决相同的技术问题，达到相同的技术效果。

本发明实施例提供的大气颗粒切割器性能的模拟方法和装置，通过计算机数值模拟的方法进行各种参数的设置，可以得到切割器内部压力场分布、速度场分布、速度矢量分布、不同粒径与捕集效率的关系，并对切割器的性能进行评价。具体地，首先运用CFD技术模拟出湍流场，分析速度矢量场中湍流分布对对颗粒运动的影响；再结合离散相模型对粒子进行追踪，得到粒子的捕集和离散数量，通过逃逸率来衡量切割器的分离效率；最后，在不同流速情况下，分析粒子的运动轨迹和逃逸率曲线，对压降进行定量分析，考察流速对能耗的影响，在避免压降提高的情况下降低临界粒径。从而优化切割器的分离效果。

本发明实施例提供的大气颗粒切割器性能的模拟方法和装置，通过Fluent软件对PM10切割器进行数值模拟，从而得出PM10切割器的压力场，速度场，速度矢量场，颗粒粒径与收集效率关系等一系列参数结果。对比实际实验结果可以发现，本次模拟的PM10切割器对于PM10的实际收集效率可以达到93％，而根据现有标准，对于PM10采样器的收集效率要求为采样器对颗粒物的捕集效率为50％时，所对应的粒子空气动力学当量直径D_a50＝10±0.5μm，捕集效率的几何偏差数σg＝1.5±0.1。根据本实施例的数据模拟可知，该PM10切割器对于PM10的捕集效率可达到93％远高于标准规定的50％。

本发明实施例所提供的大气颗粒切割器性能的模拟方法和装置的计算机程序产品，包括存储了程序代码的计算机可读存储介质，所述程序代码包括的指令可用于执行前面方法实施例中所述的方法，具体实现可参见方法实施例，在此不再赘述。

所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。

最后应说明的是：以上所述实施例，仅为本发明的具体实施方式，用以说明本发明的技术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。