CN104898155B - 一种基于中子平均能量的次临界系统次临界度测量方法 - Google Patents
一种基于中子平均能量的次临界系统次临界度测量方法 Download PDFInfo
- Publication number
- CN104898155B CN104898155B CN201510274557.9A CN201510274557A CN104898155B CN 104898155 B CN104898155 B CN 104898155B CN 201510274557 A CN201510274557 A CN 201510274557A CN 104898155 B CN104898155 B CN 104898155B
- Authority
- CN
- China
- Prior art keywords
- neutron
- energy
- average
- mrow
- subcriticality
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明提供一种基于中子平均能量的次临界系统次临界度测量方法，其包括：在堆芯布置单个或多个探测器，测量探测器位置所有中子的平均能量；通过SuperMC等程序或替代燃料等实验方法，得到探测器位置处外源中子和裂变中子各自的平均能量；通过SuperMC等程序计算外源中子相对裂变中子的中子价值；利用得到的外源中子平均能量、裂变中子平均能量和外源中子价值，从探测得到的所有中子平均能量计算得到待测的次临界度；也可以在测量之前得到中子平均能量与次临界度的刻度曲线，对于待测状态，在测得所有中子平均能量之后，利用刻度曲线得到待测状态的次临界度；在上述过程中的测量或计算中子平均能量过程，只测量或计算高能部分中子，可以显著提高测量精度。
Description
技术领域
本发明涉及测量次临界系统次临界度测量的技术领域，尤其涉及一种基于中子平均能量的次临界系统次临界度测量方法。
背景技术
加速器驱动次临界系统通过外源中子驱动裂变堆芯裂变产生裂变中子，同时释放能量。由于外源中子和裂变中子能谱和分布不同，其在堆芯不同位置处的平均能量也不相同；次临界系统中，与次临界度对应的系统有效增殖因子keff可以表示为裂变中子在所有中子的份额。
由于堆芯所有中子的平均能量除与外源中子平均能量Es和裂变中子平均能量Ef相关外，还与两种中子的比例相关，因此可能通过所有中子的平均能量求得次临界系统的次临界度。
在已有的手段中，可以通过SuperMC等程序计算给出外源中子和裂变中子单独存在时各自的平均能量，也可以通过替代燃料的方法测量出各自的平均能量。
发明内容
本发明的目的在于：提供一种基于中子平均能量的次临界系统次临界度测量方法，该方法避免传统测量方法由于其基于点堆假设而在次临界系统中不适用的弊端。本发明涉及的方法直接基于次临界度的含义，通过测量堆芯中子的平均能量进行测量。
本发明采用的技术方案为：一种基于中子平均能量的次临界系统次临界度测量方法，所述该次临界系统为带外源中子的加速器驱动次临界系统，所述方法包括以下步骤：
步骤(101)、通过程序计算或替代燃料的实验，得到堆芯中单个或多个位置P处的外源中子平均能量Es,ave；
步骤(102)、通过程序计算或临界实验，得到上述步骤中位置P处的裂变中子平均能量Ef,ave；
步骤(103)、通过程序计算，得到外源中子相对与裂变中子的中子价值
步骤(104)、在次临界系统运行过程中，测量得到上述步骤中P处位置所有中子的平均能量，通过以下关系式得到系统源有效增殖因子ks：
其中，Eave为整个堆芯内所有中子的平均能量；
步骤(105)、在得到源有效增殖因子ks之后，由以下关系式得到待测次临界度keff：
其中，其所需参数Es,ave、Ef,ave、未知情况下，通过如下步骤完成次临界度测量：
步骤一、多次改变次临界系统的状态，测量得到这些状态下探测器位置处的中子平均能量；
步骤二、对于上述步骤一中的次临界状态，通过其它测量手段得到对应的次临界度；
步骤三、将上述步骤一、二中一系列次临界状态的中子平均能量与次临界度建立刻度曲线；
步骤四、对任一次待测临界状态，测量得到探测器位置处的中子平均能量，之后利用上述步骤三中的刻度曲线得到系统的次临界度值。
其中，在计算或测量中子平均能量时，统计的中子可选择所有能量的中子或是仅选择高能量段例如大于等于7MeV中子并获得不同测量精度。
本发明与现有技术相比的优点在于：
(1)、本发明方法利用外源中子与裂变中子平均能量不同进行测量，对于加速器驱动次临界系统(ADS)具有天然适用性；
(2)、本发明方法原理上未做点堆假设，具有更高的测量精度，可以测量较深的次临界度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1表示的是本发明一实例的测量方法流程图；
图2表示的是本发明实例提供的铅冷快堆的堆芯布置示意图；
图3利用所有中子的平均能量建立刻度曲线并进行测量得到的测量值与真实值的比较；
图4表示的是利用高能中子的平均能量建立刻度曲线并进行测量得到的测量值与真实值的比较。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明提供的一种基于中子平均能量的次临界系统次临界度测量方法，其流程图如图1所示。包括如下步骤：
步骤(101)、通过SuperMC等程序关闭燃料裂变截面进行计算或通过替代燃料的实验，得到堆芯中单个或多个位置的外源中子平均能量；
步骤(102)、通过SuperMC等程序进行计算或临界实验，得到步骤(101)中同样位置的裂变中子平均能量；
步骤(103)、通过SuperMC等程序进行计算，得到外源中子相对于裂变中子的中子价值；
步骤(104)、在堆芯中与步骤(101)相同位置处，布置单个或多个探测器测量中子的平均能量；
步骤(105)、利用上述步骤中得到的参数，求得系统的次临界度。
可选的，上述步骤(101)至(103)的参数在无法得到的情况下，可以通过其它测量方法(如脉冲源法)测量多个状态(通过调节控制棒等)的次临界度，再由步骤(104)中测得的平均能量，反推得到上述步骤(101)至(103)的参数。
可选的，在不实现上述步骤(101)至(103)的前提下，可以通过其它测量方法(如脉冲源法)测量多个状态(通过调节控制棒等)的次临界度，再与步骤(104)中测得的平均能量建立刻度曲线。之后对于待测的状态，在通过步骤(104)测得中子的平均能量后，可以根据该刻度曲线进行测量。
可选的，在上述步骤(101)、(102)和(104)中，可以只得到高能中子的平均能量，并利用高能中子的平均能量进行测量，可以显著提高测量精度。
为解决传统测量方法在次临界系统中不适用的问题，本申请提供了一种不基于点堆方程而是基于平均能量的次临界度测量方法。该方法原理简单，易于实现，其原理为：
从已有的资料中，可以知道系统有效增殖因子keff与源有效增殖因子ks具有如下的关系：
上式中的表示的是外源中子相对于裂变中子的价值；
假设某一稳态时刻，堆芯内中子中有ns个由外源直接释放的中子，有nf个裂变产生的中子，其平均能量分别为Es,ave和Ef,ave。则整个堆芯内所有中子的平均能量Eave为：
由源有效增殖因子ks的定义可以知道：
由(2)式和(3)式可以得到：
Eave＝(1-ks)Es,ave+ksEf,ave＝Es,ave-ks(Es,ave-Ef,ave) (4)
在测得中子平均能量后，可由(5)式得到源有效增殖因子ks，再由(1)式可得到系统有效增殖因子keff，进而可以得到系统的次临界度。
可选的，在上述原理介绍中，可以通过只对高能中子进行分析以提高测量精度。
可选的，如果外源中子平均能量Es,ave、裂变中子平均能量Ef,ave以及外源中子价值不可知的情况下，假设：
(1)外源中子和裂变中子在次临界度改变时能谱不发生或发生很小变化，即可以认为Es,ave和Ef,ave保持不变；
(2)外源中子价值(效率)在各种工况下不发生变化，或退而求其次，只随次临界度改变而改变，这实际上是忽略了外源价值在不同工况时的改变；
在如上假设前提下，则可以通过建立中子平均能量与次临界度的刻度曲线进行测量。
为详细说明本发明申请的实施方式，下面在一个参考堆芯上进行了通过刻度曲线进行测量的模拟实验：
图2是一个铅冷次临界堆芯示意图，其中外源中子在中间组件产生，裂变中子在燃料区域产生。
(1)首先通过同步调节控制棒得到不同的状态，得到不同状态的理论次临界度：这里使用SuperMC计算得到的理论keff代替其它实验手段(如脉冲源法)的测量结果；
(2)对于步骤(1)中得到的不同状态，通过“模拟”探测器—1-6#探测器进行模拟测量，可以得到堆芯中子平均能量；
(3)使用步骤(1)和(2)中得到的理论次临界度与中子平均能量建立刻度曲线；
(4)对于次临界系统的任意不同状态(如异步调节各个控制棒的棒位)，利用1-6#探测器模拟测量中子平均能量，并利用步骤(3)中得到的刻度曲线进行测量，得到待测状态的“测量”次临界度。
对图2所示的模拟堆芯进行了模拟测量，期间利用程序直接计算模拟得到堆芯中子平均能量，结果如图3所示：其中横坐标为堆芯所有中子的平均能量，纵坐标为各个不同状态的理论次临界度。在图3中还标识出了五个待测状态相对于刻度曲线的位置。
本申请发明还提供了进一步提高测量精度的方法：
(1)首先通过同步调节控制棒得到不同的状态，得到不同状态的理论次临界度：这里使用SuperMC计算得到的理论keff代替其它实验手段(如脉冲源法)的测量结果；
(2)对于步骤(1)中得到的不同状态，通过“模拟”探测器—1-6#探测器进行模拟测量，可以得到堆芯高能中子的平均能量；
(3)使用步骤(1)和(2)中得到的理论次临界度与高能中子平均能量建立刻度曲线；
(4)对于次临界系统的任意不同状态(如异步调节各个控制棒的棒位)，利用1-6#探测器模拟测量高能中子平均能量，并利用步骤(3)中得到的刻度曲线进行测量，得到待测状态的“测量”次临界度。
对图2所示的模拟堆芯进行了模拟测量，结果如图4所示：其中横坐标为高能中子(7MeV以上)的平均能量，纵坐标为各个不同状态的理论次临界度。在图4中还标识出了五个待测状态相对于刻度曲线的位置。
Claims (3)
1.一种基于中子平均能量的次临界系统次临界度测量方法，所述该次临界系统为带外源中子的加速器驱动次临界系统，其特征在于：所述方法包括以下步骤：
步骤(101)、通过程序计算或替代燃料的实验，得到堆芯中单个或多个位置P处的外源中子平均能量Es,ave；
步骤(102)、通过程序计算或临界实验，得到上述步骤中位置P处的裂变中子平均能量Ef,ave；
步骤(103)、通过程序计算，得到外源中子相对于裂变中子的中子价值
步骤(104)、在次临界系统运行过程中，测量得到上述步骤中P处位置所有中子的平均能量，通过以下关系式得到系统源有效增殖因子ks：
<mrow>
<msub>
<mi>k</mi>
<mi>s</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>f</mi>
<mo>,</mo>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
其中，Eave为整个堆芯内所有中子的平均能量；
步骤(105)、在得到源有效增殖因子ks之后，由以下关系式得到待测次临界度keff：
2.根据权利要求1所述的基于中子平均能量的次临界系统次临界度测量方法，其特征在于：其所需参数Es,ave、Ef,ave、未知情况下，通过如下步骤完成次临界度测量：
步骤一、多次改变次临界系统的状态，测量得到这些状态下探测器位置处的中子平均能量；
步骤二、对于上述步骤一中的次临界状态，通过其它测量手段得到对应的次临界度；
步骤三、将上述步骤一、二中一系列次临界状态的中子平均能量与次临界度建立刻度曲线；
步骤四、对任一次待测临界状态，测量得到探测器位置处的中子平均能量，之后利用上述步骤三中的刻度曲线得到系统的次临界度值。
3.根据权利要求1或2所述的基于中子平均能量的次临界系统次临界度测量方法，其特征在于：在计算或测量中子平均能量时，统计的中子可选择所有能量的中子或是仅选择高能量段即大于等于7MeV中子并获得不同测量精度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510274557.9A CN104898155B (zh) | 2015-05-26 | 2015-05-26 | 一种基于中子平均能量的次临界系统次临界度测量方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510274557.9A CN104898155B (zh) | 2015-05-26 | 2015-05-26 | 一种基于中子平均能量的次临界系统次临界度测量方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104898155A CN104898155A (zh) | 2015-09-09 |
| CN104898155B true CN104898155B (zh) | 2018-03-27 |
Family
ID=54030915
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510274557.9A Expired - Fee Related CN104898155B (zh) | 2015-05-26 | 2015-05-26 | 一种基于中子平均能量的次临界系统次临界度测量方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104898155B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107689256B (zh) * | 2017-09-19 | 2019-06-07 | 中国核动力研究设计院 | 一种核反应堆次临界度测量方法 |
| KR102034830B1 (ko) * | 2018-01-05 | 2019-10-21 | 한국수력원자력 주식회사 | 중성자 고출력 검출기 고장 모니터링 방법 |
| CN110739091B (zh) * | 2019-10-22 | 2022-05-10 | 上海核工程研究设计院有限公司 | 一种核电厂次临界度直接测量系统及硬件平台和测量方法 |
| CN110988970A (zh) * | 2019-12-25 | 2020-04-10 | 中国工程物理研究院激光聚变研究中心 | 可同步控制的大动态脉冲中子测量系统 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101105986A (zh) * | 2007-08-03 | 2008-01-16 | 蔡光明 | 反应堆反应性测量方法 |
| US20100150295A1 (en) * | 2008-12-11 | 2010-06-17 | Westinghouse Electric Company Llc | Subcritical reactivity measurement method |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011120555A1 (en) * | 2010-03-29 | 2011-10-06 | Jacobs E&C Limited | Accelerator-driven nuclear system with control of effective neutron multiplication coefficent |
| JP5762840B2 (ja) * | 2011-06-17 | 2015-08-12 | 株式会社東芝 | 未臨界度測定方法および装置 |
-
2015
- 2015-05-26 CN CN201510274557.9A patent/CN104898155B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101105986A (zh) * | 2007-08-03 | 2008-01-16 | 蔡光明 | 反应堆反应性测量方法 |
| US20100150295A1 (en) * | 2008-12-11 | 2010-06-17 | Westinghouse Electric Company Llc | Subcritical reactivity measurement method |
Non-Patent Citations (1)
| Title |
|---|
| 启明星1号装置中子通量能谱与中子平均能量计算;栗再新 等;《原子能科学技术》;20060331;第40卷(第2期);168-171 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104898155A (zh) | 2015-09-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104898155B (zh) | 一种基于中子平均能量的次临界系统次临界度测量方法 | |
| CN104808168B (zh) | 一种基于分段电流技术的电能表快速校验方法 | |
| CN106898394B (zh) | 一种wwer六角形栅格堆芯的控制棒价值测量方法 | |
| CN108172312A (zh) | 核电站堆外核仪表系统的轴向功率偏差的校准方法 | |
| Fujii et al. | Comparison of improved TMD and CGC frameworks in forward quark dijet production | |
| CN103823183B (zh) | 一种同步发电机定子开路瞬变时间常数的测量方法 | |
| CN103106314B (zh) | 太阳能光伏电源输出功率的时序概率建模方法 | |
| Gonzalez-Alonso et al. | Updated determination of chiral couplings and vacuum condensates from hadronic τ decay data | |
| Ave et al. | A generalized description of the signal size in extensive air shower detectors and its applications | |
| CN106782710B (zh) | 一种反应堆测瞬发中子衰减参数的数据处理方法 | |
| CN104934083A (zh) | 一种测量缓发中子有效份额的方法 | |
| CN107103189B (zh) | 一种进行反应堆临界曲率搜索的方法 | |
| CN109192341A (zh) | 基于三维时空动力学的大反应性测量方法 | |
| Ivanov et al. | Proton recoil energy and angular distribution of neutron radiative β-decay | |
| Li et al. | Resonance elastic scattering and interference effects treatments in subgroup method | |
| Dias et al. | Measurements of k0 and Q0 values for 64Zn (n, γ) 65Zn and 68Zn (n, γ) 69mZn reactions with covariance analysis | |
| CN106126926B (zh) | 一种能够处理温度分布效应的子群共振计算方法 | |
| Pavón-Rodríguez et al. | Characterisation of the n_TOF 20 m beam line at CERN with the new spallation target | |
| Jain et al. | Optical-model interpretation of average total neutron cross sections in the keV region | |
| Zheng et al. | Generation of few-group constants by Monte Carlo code cosRMC | |
| Liang et al. | Application of dynamic rod worth measurement in a VVER | |
| Gustavino | THE 2H (α, γ) 6LI REACTION AT LUNA AND BIG BANG NUCLEOSYNTHETIS | |
| CN105046010A (zh) | 一种卫星杂散光的分析方法 | |
| Cherepanov et al. | Conserved vector current and decays τ−→ η (η′) π− π 0 ν τ | |
| Campbell | Theory overview of electroweak physics at hadron colliders |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Yang Yingkun Inventor after: Chang Bo Inventor after: Liu Chao Inventor after: Sun Yanting Inventor after: Zhang Zhenyu Inventor before: Yang Yingkun Inventor before: Chang Bo Inventor before: Liu Chao Inventor before: Sun Yanting |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180327 Termination date: 20210526 |