一种重构反应堆时空分布模型的方法和系统技术方案

本发明专利技术提供一种重构反应堆时空分布模型的方法和系统，该方法包括如下步骤:S1、根据各探测器样本的权重，对给定的探测器样本集合进行多次抽样，获得抽样探测器样本集合；S2、对抽样探测器样本集合使用重构函数，获得抽样探测器样本集合的时空分布函数；S3、利用时空分布函数计算获得该分布下的给定的探测器样本集合中的各探测器样本的计算值；S4、根据各探测器样本的计算值计算所述重构函数的可靠性；S5、在可靠性小于预设的参考值时，更新所述给定探测器样本集合中的各探测器样本的权重，直至重构函数的可靠性大于参考值；S6、根据所述历次迭代的重构函数获得最终的重构函数。本发明专利技术能够显著提高重构函数的精度。

【技术实现步骤摘要】
一种重构反应堆时空分布模型的方法和系统
本专利技术涉及核反应堆的监测测量领域，尤其涉及一种重构反应堆时空分布模型的方法和系统。
技术介绍
目前堆芯监测或者堆芯运行诊断(或者事故诊断等)都依赖于各种测量装置，或称为探测器，典型探测器包括：核电厂反应堆一、二回路系统的各种温度、压力、液位或中子通量或中子能谱等测量，还包括其他的物理量测量仪器，例如电厂辐射剂量监测所需的相关测量仪器等。首先，各类探测器本身都有其局限性。例如探测器的测量精度问题，以及探测器的空间分布等。即探测器不可能百分百准确无误，而且不可能在所知影响区域的所有位置都布置探测器。其次，探测器表征的物理量，并非核电厂真正关注的物理量。例如探测器通常得到测量电流或者电压信号等，而电厂关注的可能是其信号的本质来源如中子通量水平等。因此物理量之间的转换因子需要定义。显然，该转换因子，只能通过理论计算或者大量的经验数据拟合公式计算得到。再次，所需探测的物理量本身的具体属性差异较大而且往往未能事先可知。例如有些物理量本身在空间或时间上变化比较缓慢，就有光滑或者其特征函数的偏导数是连续等等特征。然而有些物理量本身在空间或时间上呈现跳跃性或剧烈变化。这种属性的未知性对于监测要求大幅提高，典型的大部分的重构函数等能有效重构在变化缓慢的时空分布，但对剧烈变化的时空分布却需要做更多的处理(典型地，需要进行时空变换等，使其更平缓)。最后，在综合考虑探测器局限性、物理量的内在属性和物理量转换关系等特点后，需要定义全局重构函数。该算法将探测器的信号转化为堆芯监测或者诊断的物理量的时空分布。例如事故诊断，在获取局部的某些特征物理量的时空测量值后，根据算法推算潜在来源事故特征的时空分布或概率分布等。又如，在堆芯运行监测时，在获取堆芯局部测量点的探测器测量后，经过重构函数，计算获得关注特征量的三维时空分布等等。为提高堆芯运行监测或事故监测相关物理量的时空分布的精度，往往从两个角度出发：(1)从探测器角度：提高探测器本身的精度，获取精度更高测量值。探测器布置更合理。(2)从重构函数的角度，采用更高效率或更高精度的优化算法等等，获得更精确或者经过权衡后更符合实际的算法。在工程应用上，基于各种原理和假设，已经开发了多种类型的重构函数，这些算法各有利弊：有些算法能够有效的减少探测器的不确定度和模拟全时空区域的整体变化趋势，但却忽略了局部高阶精细的变化；有些算法能够模拟高度扭曲变化的局部趋势，却对探测器的布置和探测器的失效特别敏感；有些算法能够实现高精度的时空区域的模拟，却需要耗费大量的计算资源或者计算收敛性不够好；有些算法适用于反应堆的趋势变化不明显的场景；有些算法，却能适用于反应堆的异常运行工况的变化。
技术实现思路
本专利技术所要解决的技术问题在于，提供一种重构反应堆时空分布模型方法和系统，该方法和系统通过逐渐加大误差较大的探测器的权重，降低误差较小的探测器的权重，从而提高重构结果时空分布的精度。为了解决上述技术问题，本专利技术提供一种重构反应堆时空分布模型的方法，包括如下步骤：步骤S1、根据各探测器样本的权重，对给定的探测器样本集合进行多次抽样，获得抽样探测器样本集合；步骤S2、对所述抽样探测器样本集合使用重构函数，获得所述抽样探测器样本集合的时空分布函数；步骤S3、利用所述时空分布函数计算获得所述分布下的所述给定的探测器样本集合中的各探测器样本的计算值；步骤S4、根据所述各探测器样本的计算值计算所述重构函数的可靠性；步骤S5、在所述可靠性小于预设的参考值时，更新所述给定样本集合中的各探测器的权重，并重复所述步骤S1-S4，直至所述重构函数的可靠性大于所述参考值；步骤S6、根据所述历次迭代的重构函数获得最终的重构函数。其中，所述方法还包括：获取所述给定的探测器样本集合，并假设所述给定的探测器样本集合中的每个探测器样本的初始权重相等。其中，所述步骤S3具体包括：采用下式获得所述给定的探测器样本集合中的各探测器样本的计算值：其中，为给定的探测器样本集合中的第k个探测器样本的计算值，Z(s,t)为所述时空分布函数，f1为所述第k个探测器的裂变反应截面分布。其中，所述步骤S4具体包括，步骤S41、根据所述各探测器样本的计算值计算所述给定探测器样本集合中的各探测器样本的计算值的可靠性；步骤S42、根据所述各探测器样本的计算值的可靠性计算所述重构函数的可靠性。其中，所述步骤S41具体包括：其中，sigk是给定探测器样本集合中的第k个探测器样本计算值的可靠性，Uk是给定的探测器样本集合中的第k个探测器样本的值，为给定的探测器样本集合中的第k个探测器样本的计算值，δ为一给定的拟合误差值。其中，所述步骤S42具体包括，采用下式计算所述重构函数的可靠性其中，εn是第n个重构函数的可靠性，pk为给定探测器样本中的第k个探测器样本的权重。其中，所述优化所述重构函数具体包括：根据所述给定的探测器样本集合中的各探测器样本的可靠性更新所述给定样本集合中的各探测器样本的权重。其中，所述根据所述各探测器样本的可靠性更新所述各探测器样本的权重具体包括：其中，αn是权重函数，是给定的探测器样本集合中的第k个探测器样本更新后的权重，是给定的探测器样本集合中第k个探测器样本的当前权重。其中，采用下式计算所述最终的重构函数：其中，f为最终的重构函数，N为历次迭代次数，fn为第n次迭代的重构函数。本专利技术还提供一种重构反应堆时空分布模型的系统，包括：采集模块，用于根据各探测器样本的权重，对给定的探测器样本集合进行多次抽样，获得抽样探测器样本集合；时空分布函数计算模块，用于对所述抽样探测器样本集合使用重构函数，获得所述抽样探测器样本集合的时空分布函数；第一计算模块，用于利用所述时空分布函数计算获得所述分布下的所述给定的探测器样本集合中的各探测器样本的计算值；第二计算模块，用于根据所述各探测器样本的计算值计算所述重构函数的可靠性；重构函数优化模块，用于在所述可靠性小于预设的参考值时，优化所述重构函数，将优化后的重构函数代入所述第一计算模块中进行计算，直至所述重构函数的可靠性大于所述参考值；重构函数计算模块，用于根据所述历次迭代的重构函数获得最终的重构函数。本专利技术实施例的有益效果在于：本专利技术实施例在获得重构函数的过程中，逐渐加强了误差较大的探测器的权重，降低了误差较小的探测器的权重，从而从全局上降低了整体的误差，其次，从探测器的本质出发，该算法增加了优质探测器的权重，降低了较大确定度探测器的权重，从探测器的空间分布来讲，该算法增加了合理探测器布置区域的权重，最后，通过多次取样的方式，能够通过平均的手段有效降低探测器的测量不确定性，减轻重构函数空间拟合时的过拟合现象的影响，从而使得最终重构结果时空分布具有更高精度。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1是本专利技术实施例一种基于多种重构函数重构反应堆时空分布模型的方法的主流程示意图。具体实施方式以下各实施例的说明是参考附图，用以示例本专利技术可以用以实施的特定实施例。本文档来自技高网...

【技术保护点】
1.一种重构反应堆时空分布模型的方法，包括如下步骤:步骤S1、根据给定的探测器样本集合中的各探测器样本的权重，对给定的探测器样本集合进行多次抽样，获得抽样探测器样本集合；步骤S2、对所述抽样探测器样本集合使用重构函数，获得所述抽样探测器样本集合的时空分布函数；步骤S3、利用所述时空分布函数计算获得所述时空分布函数下的所述给定的探测器样本集合中的各探测器样本的计算值；步骤S4、根据所述各探测器样本的计算值计算所述重构函数的可靠性；步骤S5、在所述可靠性小于预设的参考值时，更新所述给定的探测器样本集合中的各探测器样本的权重，并重复所述步骤S1‑S4，直至所述重构函数的可靠性大于所述参考值；步骤S6、根据所述历次迭代的重构函数获得最终的重构函数。

【技术特征摘要】
1.一种重构反应堆时空分布模型的方法，包括如下步骤:步骤S1、根据给定的探测器样本集合中的各探测器样本的权重，对给定的探测器样本集合进行多次抽样，获得抽样探测器样本集合；步骤S2、对所述抽样探测器样本集合使用重构函数，获得所述抽样探测器样本集合的时空分布函数；步骤S3、利用所述时空分布函数计算获得所述时空分布函数下的所述给定的探测器样本集合中的各探测器样本的计算值；步骤S4、根据所述各探测器样本的计算值计算所述重构函数的可靠性；步骤S5、在所述可靠性小于预设的参考值时，更新所述给定的探测器样本集合中的各探测器样本的权重，并重复所述步骤S1-S4，直至所述重构函数的可靠性大于所述参考值；步骤S6、根据所述历次迭代的重构函数获得最终的重构函数。2.根据权利要求1所述的方法，其特征在于，所述方法还包括：获取所述给定的探测器样本集合，并假设所述给定的探测器样本集合中的每个探测器样本的初始权重相等。3.根据权利要求2所述的方法，其特征在于，所述步骤S3具体包括：采用下式获得所述给定的探测器样本集合中的各探测器样本的计算值：其中，为给定的探测器样本集合中的第k个探测器样本的计算值，Z(s,t)为所述时空分布函数，f1为所述第k个探测器的裂变反应截面分布。4.根据权利要求3所述的方法，其特征在于，所述步骤S4具体包括，步骤S41、根据所述各探测器样本的计算值计算所述给定探测器样本集合中的各探测器样本的计算值的可靠性；步骤S42、根据所述各探测器样本的计算值的可靠性计算所述重构函数的可靠性。5.根据权利要求4所述的方法，其特征在于，所述步骤S41具体包括：其中，sigk是给定探测器样本集合中的第k个探测器样本计算值的可靠性，Uk是给定的探测器样本集合中的第k个探测器样本的值，为给定的...

【专利技术属性】
技术研发人员：李文淮，张香菊，仇若翔，李晓，丁鹏，王军令，段承杰，郑颖，
申请(专利权)人：中广核研究院有限公司，中国广核集团有限公司，中国广核电力股份有限公司，
类型：发明
国别省市：广东,44