本发明涉及军控核查技术领域，特别是涉及一种密闭容器核材料质量的计算方法。
背景技术：
在关系世界安全的核裁军、核安全、核保障三个重要领域里，都面临着一个共同的科学技术难题，即对贮存在密闭容器内核材料的识别、认证、性能检测和属性测量，而质量属性是其最重要的属性之一。基于3he探测器的热中子多重性探测方法(neutronmultiplicitycounting，nmc)是较为经典的方法，但是由于3he材料价格昂贵等原因，逐渐出现了许多新的探测方法。由于近些年3he气体的产量降低，需求量上升，因此其价格上升较快，目前价格大体为4万元/升，生产一套nmc探测设备价格大体在百万元到千万元级别。但是液体闪烁体材料相对价格便宜，因此此替代方法的经济适应性较好。而且nmc方法在探测中子时，需要先将中子慢化而后测量，这样就不能保留中子探测的时间信息。
上述基于3he探测器的热中子多重性探测方法在实践中已经得到了一定的应用，但是由于其所用材料3he方面的高成本的缺点以及方法本身慢化中子导致无法保存中子信号的时间信息，因此该方法实用性有待提高。
技术实现要素：
本发明的目的是提供一种密闭容器核材料质量的计算方法，以克服现有基于3he探测器的热中子多重性探测方法存在的成本高、无法保存中子信号的时间信息的问题，进而提高密闭容器核材料质量的测量精度。
为实现上述目的，本发明提供了一种密闭容器核材料质量的计算方法，所述计算方法包括：
获取液体闪烁体探测器探测到中子数的中子脉冲时间序列；
利用前景计数器对所述中子脉冲时间序列进行前景重数计数，得到前景重数分布，所述前景重数分布包含与触发中子信号相关的中子信号的计数和与所述触发中子信号无关的中子信号的计数；
利用背景计数器对所述中子脉冲时间序列进行背景重数计数，得到背景重数分布，所述背景重数分布仅包含与所述触发中子信号无关的中子信号的计数；
根据所述前景重数分布和所述背景重数分布确定真符合重数分布，所述真符合重数分布仅包含与触发中子信号相关的中子信号的计数；
根据样品参数和所述真符合重数分布确定真符合重数的总体矩，所述总体矩为所述样品参数与所述真符合重数分布之间的函数关系，所述样品参数包括自发裂变率、样品增殖系数和(α，n)反应中子数与自发裂变中子数的比例系数，所述(α，n)反应为重原子核α衰变产生的α粒子碰撞轻原子核发生的反应；
根据所述前景重数分布、所述背景重数分布和所述真符合重数分布确定真符合重数的样品矩，所述样品矩为对所述中子脉冲时间序列的计数与所述所述真符合重数分布之间的函数关系；
根据所述总体矩和所述样品矩确定总中子计数率、二重符合中子计数率和三重符合计数率；
根据所述总中子计数率、所述二重符合中子计数率和所述三重符合计数率求解所述自发裂变率、所述样品增殖系数和所述比例系数；
根据所述自发裂变率和核材料每秒的自发裂变次数计算核材料的质量。
可选的，所述利用前景计数器对所述中子脉冲时间序列进行前景重数计数，得到前景重数分布，具体包括：
获取前景时间内所述中子脉冲时间序列中中子信号个数，记为前景重数，所述前景重数对应一前景计数器，所述前景时间为自触发信号后时刻t+短延时pd至所述时刻t+所述短延时pd+符合计数门宽g的时间段；
每次计数结束所述前景重数对应的前景计数器的计数加1，得到前景计数器组；
根据所述前景计数器组中前景重数的分布规律，确定前景重数分布。
可选的，所述利用背景计数器对所述中子脉冲时间序列进行背景重数计数，得到背景重数分布，具体包括：
获取背景时间内所述中子脉冲时间序列中中子信号个数，记为背景重数，所述背景重数对应一背景计数器，所述背景时间为自触发信号后时刻t+长延时ld至所述时刻t+所述长延时ld+所述符合计数门宽g的时间段，所述长延时ld大于所述短延时pd+所述符合计数门宽g；
每次计数结束所述背景重数对应的背景计数器的计数加1，得到背景计数器组；
根据所述背景计数器组中背景重数的分布规律，确定背景重数分布。
可选的，所述根据样品参数和所述真符合重数分布确定真符合重数的总体矩，具体包括：
根据所述前景重数分布和所述背景重数分布确定样品发射中子重数的前一阶阶乘矩h′(1)、二阶阶乘矩h″(1)和三阶阶乘矩h″′(1)：
其中，m是增殖系数，α是(α，n)反应中子数与自发裂变中子数的比例系数，vsf,1自发裂变重数的一阶阶乘矩，vsf,2自发裂变重数的二阶阶乘矩，vsf,3自发裂变重数的三阶阶乘矩，vi1诱发裂变重数的一阶阶乘矩，vi2诱发裂变重数的二阶阶乘矩，vi3诱发裂变重数的三阶阶乘矩；
根据所述样品发射中子重数与所述样品发射中子重数的前一阶阶乘矩h′(1)、二阶阶乘矩h″(1)和三阶阶乘矩h″′(1)，确定一个初级事件探测中子重数的前三阶阶乘矩，所述一个初级事件探测中子重数的前三阶阶乘矩为单次源事件探测重数的前一阶阶乘矩d′(1)、二阶阶乘矩d″(1)和三阶阶乘矩d″′(1)：
d′(1)＝h′(1)e′(1)＝ε(1+κ)h′(1)
d″(1)＝h″(1)(e′(1))2+h′(1)e″(1)＝ε2(1+κ)2h″(1)+2εκh′(1)
d″′(1)＝h″′(1)(e′(1))3+3h″(1)e″(1)e′(1)+h′(1)e″′(1)＝ε3(1+κ)3h″′(1)+6ε2κ(1+κ)h″(1)
其中，e′(1)表示每一个样品发射中子产生信号数的一阶阶乘矩，e″(1)表示每一个样品发射中子产生信号数的二阶阶乘矩，e″′(1)表示每一个样品发射中子产生信号数的三阶阶乘矩，ε表示探测效率，κ表示散射串扰率，所述源事件为自发裂变和(α，n)反应；
根据所述样品发射中子重数的前一阶阶乘矩h′(1)、二阶阶乘矩h″(1)、三阶阶乘矩h″′(1)和所述初级事件探测中子重数的前三阶阶乘矩确定真符合重数的总体矩，所述真符合重数的总体矩为：
其中，f(t)表示探测的中子信号时间分布，f(s)也表示探测的中子信号时间分布，s表示时间，fd表示双重计数率的门宽因子，
可选的，所述根据所述前景重数分布、所述背景重数分布和所述真符合重数分布确定真符合重数的样品矩，具体包括：
获取总触发率s＝(f+fα)d′(1)+sbkg＝(f+fα)ε(1+κ)h′(1)+sbkg；
根据所述前景重数分布获取前景重数的阶乘矩，所述前景重数的阶乘矩fk为：
根据所述总触发率和所述所述前景重数的阶乘矩确定真符合重数的样品矩，所述真符合重数的样品矩为：
可选的，所述根据所述总体矩和所述样品矩确定总中子计数率、二重符合中子计数率和三重符合计数率，具体包括：
令所述总体矩等于所述样品矩，得到总中子计数率s、二重符合中子计数率d和三重符合计数率t分别为：
s＝fε(1+κ)νsf,1(1+α)m
可选的，所述根据所述总中子计数率、所述二重符合中子计数率和所述三重符合计数率求解所述自发裂变率、所述样品增殖系数和所述比例系数，具体包括：
利用所述所述总中子计数率、所述二重符合中子计数率和所述三重符合计数率计算所述样品增殖系数m；
利用公式
利用公式
根据本发明提供的具体实施例，本发明公开了以下技术效果：本发明中采用了成本相对更低的有机液闪探测器对触发中子进行探测，降低了整个技术的成本。并且本发明所依据的探测中子重数数据是未经慢化的中子信号，因此极大地保留了中子信号的时间信息，克服了现有技术中采用基于3he探测器的热中子多重性探测方法进行密闭容器内核材料质量测量的缺点。本发明提供的基于液闪的快中子多重性分析方法的密闭容器核材料质量的计算方法的准确性、经济性更好，具有更好的推广和应用价值。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的密闭容器核材料质量的计算方法的流程图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的目的是提供一种密闭容器核材料质量的计算方法，以克服现有基于3he探测器的热中子多重性探测方法存在的成本高、无法保存中子信号的时间信息的问题，进而提高密闭容器核材料质量的测量精度。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
本发明的基本思想是：首先，以nmc的六条基本假设为依据，即(1)假设所有诱发裂变中子与其对应的源自发裂变或(α，n)反应是同时发生的，这条假设被称为超裂变假设；(2)假设对于整个样品空间，探测效率和发生裂变的可能性是均匀的；(3)假设(α，n)反应与自发裂变以及诱发裂变的中子能谱是一样的，即它们产生的中子有着相同的探测效率和诱发裂变率；(4)假设所有被俘获的中子都会引起诱发裂变；(5)假设中子多重性分析与裂变事件发出的中子能量无关联；(6)假设在样品探测器结构中的中子衰退时间近似服从单一的e指数分布。逐条进行分析，得出前四条假设可以直接套用，后两条假设进行修改，即(5)假设中子在测量装置中的衰退时间全部分布于100ns以内。这是因为快中子的探测没有聚乙烯慢化中子的过程，1mev的中子飞行1m大约需要70ns，而快中子多重性测量(fastneutronmultiplicitycounting，fnmc)测量装置样品与探测器距离一般小于30cm，所以可以假设中子衰退时间全部分布于100ns以内，这样100ns的符合门宽足够大。(6)假设每个从样品中出射的中子在第一个探测器中产生信号后都有相同的机率再在第二个探测器中产生信号，但不会引起第三个探测器的响应。然后，通过方程推导得出在fnmc中s、d、t与f、m、α的关系。使用液体闪烁体探测器对未知质量核素进行探测，得到其中子的探测时间序列，由于中子脉冲序列包含的中子数目较多，因此需要对其进行处理。经过处理后得到前景计数器和背景计数器(说明参见资料)的计数重数分布。而后经过处理，得到中子俘获时间的s、d、t方程，求解出f、m、α，从而求解出有效质量m。由于fnmc方法直接探测未经慢化的快中子，因此保留了中子信号的时间信息。
如图1所示，本发明提供的密闭容器核材料质量的计算方法，包括：
步骤101：获取液体闪烁体探测器探测到中子数的中子脉冲时间序列；
步骤102：利用前景计数器对所述中子脉冲时间序列进行前景重数计数，得到前景重数分布，所述前景重数分布包含与触发中子信号相关的中子信号的计数和与所述触发中子信号无关的中子信号的计数；
步骤103：利用背景计数器对所述中子脉冲时间序列进行背景重数计数，得到背景重数分布，所述背景重数分布仅包含与所述触发中子信号无关的中子信号的计数；
步骤104：根据所述前景重数分布和所述背景重数分布确定真符合重数分布，所述真符合重数分布仅包含与触发中子信号相关的中子信号的计数；
步骤105：根据样品参数和所述真符合重数分布确定真符合重数的总体矩，所述总体矩为所述样品参数与所述真符合重数分布之间的函数关系，所述样品参数包括自发裂变率、样品增殖系数和(α，n)反应中子数与自发裂变中子数的比例系数，所述(α，n)反应为重原子核α衰变产生的α粒子碰撞轻原子核发生的反应；
步骤106：根据所述前景重数分布、所述背景重数分布和所述真符合重数分布确定真符合重数的样品矩，所述样品矩为对所述中子脉冲时间序列的计数与所述所述真符合重数分布之间的函数关系；
步骤107：根据所述总体矩和所述样品矩确定总中子计数率、二重符合中子计数率和三重符合计数率；
步骤108：根据所述总中子计数率、所述二重符合中子计数率和所述三重符合计数率求解所述自发裂变率、所述样品增殖系数和所述比例系数；
步骤109：根据所述自发裂变率和核材料每秒的自发裂变次数计算核材料的质量。
下面详细描述核材料的质量的计算过程：
1.1fnmc方程推导
在fnmc的基本假设下，影响测量结果的未知参数可以简化为三个：自发裂变率f、样品增殖系数m和(α，n)反应中子数与自发裂变中子数的比例系数α。本部分将推导这三个样品参数与重数计数的结果的数学关系，建立fnmc的分析方程。运用fnmc分析方程，可以由测量结果反推三个样品参数。
由测量结果分析裂变材料的问题可看成参数估计问题。参数估计采用矩估计方法，本实施方式使用阶乘矩，分别推出真符合重数的总体矩(样品参数与真符合重数的关系)和样本矩(探测重数与真符合重数的关系)的表达式，令总体矩与样本矩相等，即得到参数估计方程。具体来说就是：由自发裂变中子重数分布和中子在样品内输运过程中的中子数变化性质，求得一次自发裂变的样品发射中子重数分布，再根据探测器探测和计数方式推导出真符合重数的总体矩rk(f,α,m)；依据前景计数分布、背景计数分布与真符合重数分布之间的关系，推导出由探测重数表示的真符合重数样本矩rk(f,b)。
推导中通常假定所有探测到的中子信号，除本底信号和噪声外，都可以追溯到一次自发裂变或者是一次由重原子核(heavynuclei)α衰变产生的α粒子碰撞轻原子核发生的(α，n)反应，所以把样品中发生的自发裂变和(α，n)反应称作源事件。并假定中子经过一个可能发生增殖的输运过程后从样品中发射出去，然后以一定的概率被探测系统探测到。
串扰会导致探测系统给出的中子计数大于实际的中子计数，从而影响对放射源参数的分析。根据钚裂变发射中子重数分布的统计规律，若探测系统的探测效率在10％左右，二重计数大约为一重的十分之一，三重也大约为二重的十分之一。那么若串扰率为1％，就意味着二重计数有10％的偏差，三重计数有20％的偏差(二重计数的中的每个中子都可能产生串扰)。这对于方程求解分析的准确性可能带来很大影响。所以散射串扰就是一个原来的表达式中没有考虑，而在fnmc中又无法避免的因素，要把它考虑进分析使用的表达式中。
本文采用概率母函数(probabilitygeneratingfunction，pgf)作为数学工具，使用h(z)表示发射中子重数的概率母函数，对fnmc方程的进行构建。根据概率母函数的性质，x的k阶阶乘矩e(x(x-1)…(x-k+1))就是x的pgf在1处的k阶导数。可用h′(1),h″(1),h″′(1)代替单次裂变中子的“一阶矩”、“二阶矩”、“三阶矩”，即h′(1),h″(1),h″′(1)是发射中子重数的前三阶阶乘矩。
1.1.1总体矩
即样品参数与真符合重数的关系。样品增殖系数m采用阶乘矩表达式表示为：
m＝h′(1)(1)
其中，h′(1)表示核材料中每个初级中子所致样品发射中子数的一阶矩。
本部分从物理过程出发，推导f，α，m三个样品参数与真符合重数之间的函数关系——真符合重数的总体矩rk(f,α,m)。因为由概率母函数(pgf，probabilitygeneratingfunction)求阶乘矩只需做求导运算，且母函数处理复合函数的高阶导数也特别方便。所以这里将使用概率母函数完成对rk(f,α,m)的推导。
1)源事件的概率母函数
每次发生自发裂变时发射的中子重数是独立同分布的非负整值随机变量，用概率母函数pgf表示，表达示为
式中，pi－表示发射i个中子时的概率。
核材料中自发裂变是泊松过程，泊松过程的强度f称为自发裂变率。这样，自发裂变产生中子的过程就是一个复合泊松过程。
(α，n)反应与自发裂变无关，也是泊松过程，设其强度为fα，每次(α，n)反应发射中子数为1，用概率母函数pgf表示为
sα(z)＝z(3)
将自发裂变与(α，n)反应统称为源事件，则源事件可以看作是强度为f+fα的泊松过程，每次发射的中子数的概率母函数pgf为
根据pgf的性质和参数α的定义式，一个源事件发射中子数的一阶矩、二阶矩、三阶矩分别为
由于自发裂变中子阶乘矩是可以根据核数据计算的常量，所有公式中将使用νsf,1,νsf,2,νsf,3代替ssf′(1),ssf″(1),ssf″′(1)。
2)样品发射中子数的概率母函数
若源事件产生1个中子，从样品中发射出来的中子数的pgf记为h1(u)。那么，如果源事件中产生n个中子，则根据pgf的性质，发射中子数的pgf为
hn(z)＝[h1(z)]n(8)
将样品中一个中子发生诱发裂变的概率记为pif，诱发裂变产生i个中子的可能性记为pifν(i)，则这i个中子所致样品发射中子数的pgf也为hi(z)。一次诱发裂变产生中子数的pgf为
假设中子在样品中只要被吸收就会诱发裂变，忽略中子在样品中被吸收而没有诱发裂变的可能性，根据pgf的性质
为书写方便，后文统一用h(z)代替h1(z)，即
h(z)≡h1(z)＝(1-pif)z+pifu[h(z)](11)
式(11)给出了当源事件发射一个中子时，样品所发射中子数h(z)、诱发裂变概率pif和诱发裂变发射中子数u(z)之间的关系。考虑源事件发射的中子数分布，根据pgf的性质，一次源事件产生的发射中子数pgf为
h(z)＝s[h(z)](12)
3)源事件产生的发射中子数的阶乘矩
根据pgf的性质，由式(12)可得源事件对应的样品发射中子数的一阶阶乘矩、二阶阶乘矩、三阶阶乘矩为
h′(1)＝s′(1)h′(1)(13)
h″(1)＝s″(1)[h′(1)]2+s′(1)h″(1)(14)
h″′(1)＝s″′(1)[h′(1)]3+3s″(1)h″(1)h′(1)+s′(1)h″′(1)(15)
对式(11)，求z＝1处的导数，可得
考虑到诱发裂变中子阶乘矩是可以由核数据计算得到的常量，式中用νi1代替u′(1)，下文中同样使用νi2,νi3代替u″(1),u″′(1)。
h′(1)是每个源中子在输运过程中的增殖期望，所以设样品发射中子数与源中子数之比增殖参数m为
对式(11)，求z＝1处的二阶导数，可获得中子增殖的二阶矩：
对式(11)，求z＝1处的三阶导数，可得中子增殖的三阶矩：
将式(17)、式(18)、式(19)代入式(13)、式(14)、式(15)，得到样品发射中子数前一阶矩、二阶矩、三阶矩为
4)真符合重数的总体矩
先假定同一个源事件的中子中有n个中子被探测到，在其中随机选取一个中子作为触发信号，推导这个中子触发信号的真符合计数重数的概率分布p(i|n)。
假设源事件发生于t＝0时刻，探测的中子信号时间分布为f(t)。符合计数门宽为g，短延时为pd。
在(t,t+dt)内能够探测到1个中子的可能性为
p(t|n)＝nf(t)dt(23)
以它作为触发信号符合计数门为从t+pd时刻到t+pd+g时刻，其余n-1个中子信号中的任意一个处于此符合计数门中(记为事件a)，对此次触发的计数重数产生贡献)的可能性为
因而在(t,t+dt)内探测到1个中子，并且余下的n-1个中子在符合计数门中有i个的可能性为
在整个时间范围对t做积分，即可得到随机选取n个信号中的一个作为触发信号，所记重数为i(符合计数门中有i个中子信号)的可能性p(i|n)，
快中子测量的探测重数与热中子的不同。对于热中子探测，一个发射中子被探测并产生一个信号的可能为ε，不被探测的可能性是(1-ε)，所以经典nmc测量中每个发射中子产生的中子信号重数的pgf为e(z)＝(1-ε)+εz，而快中子多重性探测还需要考虑散射串扰的影响，根据fnmc第(6)条基本假设，fnmc测量中每一个样品发射中子产生信号数的pgf为
e(z)＝(1-ε)+ε(1-κ)z+εcz2(27)
根据pgf的性质，一次源事件的探测重数为：
d(z)＝h(e(z))(28)
记一次源事件探测到的中子数为n的可能性为d(n)，则对于每个探测到的中子都会触发一次计数，则某一个触发中子属于探测到n个中子的源事件的可能性为
式中，d(n)―探测重数分布
由公式(28)、(29)，推导出单次源事件探测重数的前三阶矩测量方程。即一个初级事件探测中子重数的前三阶矩为：
某次计数重数为i的概率为
由式(31)、式(32)、式(33)、式(34)，可求得出真符合重数总体矩r(i)(信号触发关联重数)的一、二阶矩分别为：
式中，fd―双重计数率(doubles)门宽因子；ft―三重计数率(triples)的门宽因子。可以由式(36)、(37)计算。
1.1.2样品矩
即测量计数与真符合重数的关系。探测到的中子既包括源事件直接或间接诱发裂变产生的中子，又有与源事件无关的本底中子。前景计数器以探测到的中子为触发，总触发率为：
s＝(f+fα)d′(1)+sbkg＝(f+fα)ε(1+κ)h′(1)+sbkg(38)
前景重数的阶乘矩为：
式中bk―本底中子触发重数，即偶然符合重数，与背景重数分布相同；gk―源事件中子触发重数，包含真符合重数rj与偶然符合重数bk。
所以源事件中子触发重数的阶乘矩公式为，
对式(39)两边同时减bk，再将式(40)代入得
依次k＝1,2将代入式，可以获得真符合重数的一、二阶矩r'1，r'2使用fk,bk的形式表达。对应的真符合重数样本矩为：
1.1.3参数估计方程组
令真符合重数的总体矩(式(34)、(35))，与样本矩(式(42)、(43))相等，即可得到两个矩估计方程，而式反应了探测到总中子数(触发数)与样本参数之间的关系。联立三个方程，并把测量量与物理参数量分别放在等式两边，即可得到考虑散射串扰推导出的fnmc测量方程为：
s-sbkg＝(f+fα)ε(1+κ)h′(1)(44)
式(44)、式(45)、式(46)左侧就是fnmc中常用到的测量量singles，doubles和triples，简写为s，d，t。
将式(20)、式(21)、式(22)代入式(44)、式(45)、式(46)得：
s＝fε(1+κ)νsf,1(1+α)m(47)
1.2fnmc方程求解
通过对公式(47)、(48)、(49)的相互代入消除f和α，得到方程，m的解可以通过求解该一元三次方程获得。
a+bm+cm2+m3＝0(50)
其中，
求解出m后，自发裂变率f由下式求出
(α，n)反应系数α为
一旦求出裂变率f后，根据每克240pu平均每秒发生473次裂变的物理性质，可计算等效240pu质量。钚样品的有效质量可以由式给出。然后用有效质量除以同位素组成得到总质量。
下面给出本发明的计算方法在具体实例中的应用。
现以bc501-a液体闪烁体探测器并结合g4模拟仿真工具箱为例对本发明方法进一步说明。
使用g4仿真工具箱模拟bc501-a液体闪烁体探测器对240pu中子源进行探测，模拟测量82.9g金属钚(90％240pu，10％239pu，19.8g/cm3)。模拟测量中模拟了4,000,000次初级行为(全部为240pu自发裂变)，输出的数据包含922,178个探测到中子的时间。按时间先后将数据排序后，为保证数据的准确性，取用中间90％的数据(从第46,108个到第875,522个)，对应的测量时间长度为101.9s。按照上述步骤对其探测结果结果进行处理。
根据中子的能谱模拟发射中子，使其发射概率符合中子的能谱。
根据概率母函数，计算单次源事件探测重数的前三阶矩。并得到金属钚模拟测量重数计数分布。(为简化计算且在误差允许范围内，fd和ft简单取为1)
金属钚模拟测量重数计数分布表
由fnmc测量方程及相关参数求解s、d、t的值。
d＝s(f1-b1)＝1.03×103
计算得到有效质量
因样品是82.9g金属钚(90％240pu，10％239pu)的，考虑240pu的同位素组成为90％，最后可求出钚的质量相对偏差为8.1％，解算效果较好。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。