发明名称 获取用于核电站辐射仿真系统的多核素等效参数的方法
摘要 本发明涉及核反应堆工程技术领域,特别涉及一种获取用于核电站辐射仿真系统的多核素等效参数的方法。本发明从基本的核素衰变理论出发,经过推理和计算,获得适用于多核素等效为单一核素的精确计算方法,即给等效核素赋予两个等效常数,一个等效常数用于描述原子核数的衰减概率,另一个等效常数用于描述放射性活度衰减的动态特性,同时考虑了具有任意运行历史的核反应堆源项数据和影响核素等效参数的时间因素,并对仿真问题的时间跨度可采用分段等效或整段等效,从而获得了用于核电站辐射仿真系统的多核素等效参数,使辐射仿真更加精确。
申请公布号 CN105373665A 申请公布日期 2016.03.02
申请号 CN201510801761.1 申请日期 2015.11.19
申请人 厦门大学 发明人 倪星河;郑剑香;谢珊;李宁
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 厦门市精诚新创知识产权代理有限公司 35218 代理人 方惠春
主权项 获取用于核电站辐射仿真系统的多核素等效参数的方法,其特征在于:包括S10,收集可靠的放射性源项数据,进入步骤S20;S20,处理收集到的放射性源项数据,明确每一大类核素所含的具体核素及其放射性活度或者比活度,进入步骤S30;S30,由核数据库获取各相应核素的半衰期及相对原子质量,并计算出相应的衰变常数λ<sub>i</sub>,进入步骤S40;S40,针对所要仿真的问题及其最大时间TP,判定是否需要采用分段等效,如果否,则进入步骤S50;如果是,则进入步骤S60;S50,计算问题时间t=TP所对应的等效活度衰减率λ<sub>ae</sub>,问题初始时刻t=0所对应的等效核衰变常数λ<sub>ne</sub>和问题时间t=TP所对应的等效核相对原子质量M<sub>e</sub>,进入步骤S80;S60,确定问题时间的分段节点t<sub>j</sub>,进入步骤S70;S70,计算各分段问题时间节点t=t<sub>j</sub>所对应的等效活度衰减率(λ<sub>ae</sub>)<sub>j</sub>,等效核衰变常数(λ<sub>ne</sub>)<sub>j</sub>和等效核相对原子质量(M<sub>e</sub>)<sub>j</sub>,进入步骤S80;S80,将等效活度衰减率λ<sub>ae</sub>和等效核衰变常数λ<sub>ne</sub>代入等效核活度函数<img file="FDA0000851676580000011.GIF" wi="493" he="79" />中,或将等效活度衰减率(λ<sub>ae</sub>)<sub>j</sub>和等效核衰变常数(λ<sub>ne</sub>)<sub>j</sub>代入等效核活度函数<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>A</mi><mi>e</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mrow><mo>&lsqb;</mo><mrow><msub><mrow><mo>(</mo><msub><mi>&lambda;</mi><mrow><mi>n</mi><mi>e</mi></mrow></msub><mo>)</mo></mrow><mn>1</mn></msub><mo>&CenterDot;</mo><msub><mrow><mo>(</mo><msub><mi>N</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>1</mn></msub></mrow><mo>&rsqb;</mo></mrow><mo>&CenterDot;</mo><mrow><mo>&lsqb;</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mrow><mo>(</mo><msub><mi>&lambda;</mi><mrow><mi>a</mi><mi>e</mi></mrow></msub><mo>)</mo></mrow><mn>1</mn></msub><mo>&CenterDot;</mo><mi>t</mi></mrow></msup><mo>&rsqb;</mo></mrow></mrow></mtd><mtd><mrow><mo>(</mo><mn>0</mn><mo>&le;</mo><mi>t</mi><mo>&lt;</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mtable><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr></mtable></mtd><mtd><mtable><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr></mtable></mtd></mtr><mtr><mtd><mrow><mrow><mo>&lsqb;</mo><mrow><msub><mrow><mo>(</mo><msub><mi>&lambda;</mi><mrow><mi>n</mi><mi>e</mi></mrow></msub><mo>)</mo></mrow><mi>j</mi></msub><mo>&CenterDot;</mo><msub><mrow><mo>(</mo><msub><mi>N</mi><mn>0</mn></msub><mo>)</mo></mrow><mi>j</mi></msub></mrow><mo>&rsqb;</mo></mrow><mo>&CenterDot;</mo><mrow><mo>&lsqb;</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mrow><mo>(</mo><msub><mi>&lambda;</mi><mrow><mi>a</mi><mi>e</mi></mrow></msub><mo>)</mo></mrow><mi>j</mi></msub><mo>&CenterDot;</mo><mi>t</mi></mrow></msup><mo>&rsqb;</mo></mrow></mrow></mtd><mtd><mrow><mo>(</mo><msub><mi>t</mi><mrow><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>&le;</mo><mi>t</mi><mo>&lt;</mo><msub><mi>t</mi><mi>j</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mtable><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr></mtable></mtd><mtd><mtable><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr></mtable></mtd></mtr><mtr><mtd><mrow><mrow><mo>&lsqb;</mo><mrow><msub><mrow><mo>(</mo><msub><mi>&lambda;</mi><mrow><mi>n</mi><mi>e</mi></mrow></msub><mo>)</mo></mrow><mi>Y</mi></msub><mo>&CenterDot;</mo><msub><mrow><mo>(</mo><msub><mi>N</mi><mn>0</mn></msub><mo>)</mo></mrow><mi>Y</mi></msub></mrow><mo>&rsqb;</mo></mrow><mo>&CenterDot;</mo><mrow><mo>&lsqb;</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mrow><mo>(</mo><msub><mi>&lambda;</mi><mrow><mi>a</mi><mi>e</mi></mrow></msub><mo>)</mo></mrow><mi>Y</mi></msub><mo>&CenterDot;</mo><mi>t</mi></mrow></msup><mo>&rsqb;</mo></mrow></mrow></mtd><mtd><mrow><mo>(</mo><msub><mi>t</mi><mrow><mi>Y</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>&le;</mo><mi>t</mi><mo>&le;</mo><mi>T</mi><mi>P</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000851676580000012.GIF" wi="982" he="382" /></maths>中,判断等效核活度函数与各核素的叠加活度函数的误差是否满足设定的误差要求,如果是,则进入步骤S90;如果否,则判断前述计算过程是否是采用整段等效,如果否,则返回步骤S60重新进行时段划分,如果是,则判断是否仍采用整段等效,如果否,则返回步骤S60,如果是,则进入步骤S90;其中,N<sub>0</sub>为问题初始时刻t=0时的放射性原子核数目,(N<sub>0</sub>)<sub>j</sub>为第j时段内的初始时刻的放射性原子核数目,(N<sub>0</sub>)<sub>1</sub>=N<sub>0</sub>;S90,输出最终的等效活度衰减率λ<sub>ae</sub>或(λ<sub>ae</sub>)<sub>j</sub>,等效核衰变常数λ<sub>ne</sub>或(λ<sub>ne</sub>)<sub>j</sub>和等效核相对原子质量M<sub>e</sub>或(M<sub>e</sub>)<sub>j</sub>,即获得多核素等效参数。
地址 361000 福建省厦门市思明南路422号