CN106782709A - 一种零功率物理试验等温温度系数测量值修正方法 - Google Patents

一种零功率物理试验等温温度系数测量值修正方法 Download PDF

Info

Publication number
CN106782709A
CN106782709A CN201611052690.0A CN201611052690A CN106782709A CN 106782709 A CN106782709 A CN 106782709A CN 201611052690 A CN201611052690 A CN 201611052690A CN 106782709 A CN106782709 A CN 106782709A
Authority
CN
China
Prior art keywords
boron concentration
temperature coefficient
correction
isothermal temperature
critical
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.)
Granted
Application number
CN201611052690.0A
Other languages
English (en)
Other versions
CN106782709B (zh
Inventor
蔡光明
刘国明
高鑫
何子帅
郑东佳
程宏亚
胡娟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CNNC Fujian Nuclear Power Co Ltd
Original Assignee
CNNC Fujian Nuclear Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by CNNC Fujian Nuclear Power Co Ltd filed Critical CNNC Fujian Nuclear Power Co Ltd
Priority to CN201611052690.0A priority Critical patent/CN106782709B/zh
Publication of CN106782709A publication Critical patent/CN106782709A/zh
Application granted granted Critical
Publication of CN106782709B publication Critical patent/CN106782709B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21C—NUCLEAR REACTORS
    • G21C17/00—Monitoring; Testing ; Maintaining
    • G21C17/10—Structural combination of fuel element, control rod, reactor core, or moderator structure with sensitive instruments, e.g. for measuring radioactivity, strain
    • G—PHYSICS
    • G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21C—NUCLEAR REACTORS
    • G21C17/00—Monitoring; Testing ; Maintaining
    • G21C17/10—Structural combination of fuel element, control rod, reactor core, or moderator structure with sensitive instruments, e.g. for measuring radioactivity, strain
    • G21C17/112—Measuring temperature
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00—Energy generation of nuclear origin
    • Y02E30/30—Nuclear fission reactors

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Plasma & Fusion (AREA)
  • General Engineering & Computer Science (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

本发明属于核反应堆物理技术领域，具体涉及一种零功率物理试验等温温度系数测量值控制棒棒位‑临界硼浓度修正方法，利用该方法可以有效地减小等温温度系数测量中的修正误差。本发明的方法采用三维计算模型更加接近堆芯运行状态，计算出各堆芯状态参数同时考虑了棒位‑硼浓度的影响，大大减小等温温度系数测量中的修正误差。

Description

一种零功率物理试验等温温度系数测量值修正方法
技术领域
本发明属于核反应堆物理技术领域，具体涉及一种零功率物理试验等温温度系数测量值控制棒棒位-临界硼浓度修正方法，利用该方法可以有效地减小等温温度系数测量中的修正误差。
背景技术
压水堆核电站要求，反应堆在各种功率水平下运行时，慢化剂温度系数必须为负值或零，使反应堆具有负反馈特性。零功率物理试验期间，要求进行慢化剂温度系数的测量和验证，但是由于测量上的困难，慢化剂温度系数不是直接测量得到，而是通过测量等温温度系数间接得到。等温温度系数定义为慢化剂、燃料包壳以及燃料芯块的温度单位变化所引起的反应性变化，定义式如下：
αiso＝Δρ/ΔT (1)
其中，△ρ是由于温度变化引起的反应性变化量，单位pcm；△T是温度变化量，单位℃。
反应堆处于热态零功率时，一二回路处于热平衡状态，为了测量等温温度系数，通过调节大气释放阀开度来控制蒸汽排放量，使得慢化剂温度呈线性变化。反应性的变化通过反应性仪测量、记录仪记录，绘制反应性和慢化剂温度变化的曲线，曲线斜率(Δρ/ΔT)就是等温温度系数。
为最大限度地减小试验误差，将试验分成对称的几部分，包括四个过程，分别为冷却降温、加热升温、冷却降温、加热升温四个过程，使冷却剂平均温度在参考温度(Reference Temperature)Tref±1℃的范围变化，如图1所示。每个过程可得到一个等温温度系数测量值，由于每个测量过程的控制棒棒位、冷却剂平均温度、临界硼浓度与理论状态有偏差，需要进行修正。经过修正后得到测量修正值，最后对4个测量修正值取平均得到平均值。
试验时的堆芯状态与设计给出的标准状态有一定的偏离，为了与设计状态下的理论值进行比较，需要对每个测量值用理论修正公式进行修正。传统的修正方法为两维核设计计算的结果，对控制棒棒位和临界硼浓度进行简单线性插值计算，与堆芯实际运行不符合，尤其是控制棒的插入对等温温度系数的影响实际上是非线性的，因此此修正带来的误差极大。修正引入的误差，不仅包括线性插值计算方法的误差，还包括程序计算模型与堆芯实际运行的误差。
发明内容
本发明要解决的技术问题是提供一种更加精确的修正方法进一步减小等温温度系数测量中的修正误差，通过使用这种方法得到更加准确的等温温度系数。
为了实现这一目的，本发明采取的技术方案是：
一种零功率物理试验等温温度系数测量值修正方法，修正采用线性插值，在临界状态下，硼浓度和控制棒棒位相互影响，对控制棒棒位的修正包含相应棒位下临界硼浓度的修正，用三维程序直接计算控制棒在不同位置及临界硼浓度下的等温温度系数计算值，从而包含控制棒位和硼浓度的偏差修正；包括以下步骤：
步骤一：使用三维程序计算得到控制棒棒位-临界硼浓度对等温温度系数的综合修正量
采用三维计算以避免堆芯曲率变化对修正量的影响；
根据堆芯中不同控制棒棒位的临界硼浓度与计算的等温温度系数ISO以及棒组不同位置时，对应的临界硼浓度下的等温温度系数ITC计算值及修正量；
实际修正时，假设控制棒位在A时，实测临界硼浓度为Xt，计算值为Xc，控制棒位在B时，实测硼浓度Yt，计算值为Yc，设定不同控制棒位时，相同的临界硼浓度差，其修正值相等，从而得出下关系：(Xt-Yt)＝(Xc-Yc)；
步骤二：采用插值方法得到温度对等温温度系数测量值的修正量
在具体的零功率物理实验中，对温度的修正不用与棒位及硼浓度的修正耦合，进行单独修正；
步骤三：结合步骤一和步骤二，综合堆芯温度与控制棒棒位-临界硼浓度对等温温度系数测量值的修正量，得到修正后的等温温度系数测量值；
再通过多次的测量与修正取平均值减小误差；
本修正方法总结为关系式：
ITC测量修正值＝试验状态下的棒位下ITC测量值+(控制棒棒位-临界硼浓度修正量)+温度修正量。
进一步的，如上所述的一种零功率物理试验等温温度系数测量值修正方法，A＝225，B＝184。
本发明技术方案的有益效果在于：①修正方法简化，传统的修正方法需要对控制棒棒位、硼浓度分开修正，本发明的方法对控制棒棒位-硼浓度进行综合修正；②本发明的方法采用三维计算模型更加接近堆芯运行状态，计算出各堆芯状态参数同时考虑了棒位-硼浓度的影响，大大减小等温温度系数测量中的修正误差。
附图说明
图1是两维计算ITC随各控制棒位的变化；
图2是三维计算ITC随各控制棒位的变化。
具体实施方式
下面通过附图和具体实施例对本发明技术方案进行进一步详细说明。
以某电厂某机组某循环零功率ARO等温温度系数测量试验为例，本发明技术方案包括以下步骤：
步骤一：使用三维程序计算得到控制棒棒位-临界硼浓度对等温温度系数的综合修正量。
由于控制棒棒位和硼浓度的改变对堆芯曲率影响较大，采用三维计算可以避免堆芯曲率变化对综合修正量的影响。
如表1所示，给出了堆芯中不同控制棒棒位的临界硼浓度与计算的等温温度系数(Isothermal temperature coefficient,ISO)。表1给出了R棒组不同位置，对应的临界硼浓度下的ITC计算值及修正量。实际修正时，假设R＝225时，实测临界硼浓度为Xt，计算值为Xc，R棒插入在184时，实测硼浓度Yt，计算值为Yc，有以下关系：(Xt-Yt)＝(Xc-Yc)。R棒插入在184且临界硼时，实测值计算值R棒在225且临界硼时，实测值计算值根据表2，不同控制棒位时，相同的临界硼浓度差，其修正值接近相等。在计算时，可以通过计算得到，即关于硼浓度的修正刚好可以抵消。表1给出的修正量即包含了控制棒位的修正，也包含了硼浓度的修正。
表1各R棒组及临界硼浓度下的ITC计算结果
*相对于R插入限184步
表2各R棒组位置不同硼浓度的ITC
步骤二：采用插值方法得到温度对等温温度系数测量值的修正量。
由于在具体的零功率物理实验中，由于温度导致的堆芯曲率变化较小，因此，对温度的修正可以不用与棒位及硼浓度的修正耦合，可以进行单独修正。在表3中，试验中R棒棒位为213步，根据步骤一的表格修正值为0.214pcm/℃。
表3各插棒状态下等温温度系数随温度变化的修正量
步骤三：综合堆芯温度与控制棒棒位-临界硼浓度对等温温度系数测量值的修正量，即得到修正后的等温温度系数测量值，再通过多次的测量与修正取平均值减小误差，如表4所示。本修正方法可以总结为关系式：
ITC测量修正值＝试验状态下的棒位下ITC测量值+(控制棒棒位-临界硼浓度修正量)+温度修正量。
试验结果处理如下：
表4采用修正方法对等温温度系数的修正
在以下的具体实施例中，展现传统方法与本发明技术方案的对比：
情况1：使用传统方法计算某电厂某机组某循环零功率物理试验等温温度系数测量值。
表5使用两维修正方法计算等温温度系数的修正
图1是这种情况下两维计算ITC随各控制棒位的变化示意图。
情况2：使用三维修正方法计算某电厂某机组某循环零功率物理试验等温温度系数测量值。
表6使用三维修正方法计算等温温度系数的修正
图2是三维计算ITC随各控制棒位的变化情况示意图。

Claims (2)

1.一种零功率物理试验等温温度系数测量值修正方法，修正采用线性插值，在临界状态下，硼浓度和控制棒棒位相互影响，对控制棒棒位的修正包含相应棒位下临界硼浓度的修正，用三维程序直接计算控制棒在不同位置及临界硼浓度下的等温温度系数计算值，从而包含控制棒位和硼浓度的偏差修正；其特征在于，包括以下步骤：
步骤一：使用三维程序计算得到控制棒棒位-临界硼浓度对等温温度系数的综合修正量
采用三维计算以避免堆芯曲率变化对修正量的影响；
根据堆芯中不同控制棒棒位的临界硼浓度与计算的等温温度系数ISO以及棒组不同位置时，对应的临界硼浓度下的等温温度系数ITC计算值及修正量；
实际修正时，假设控制棒位在A时，实测临界硼浓度为Xt，计算值为Xc，控制棒位在B时，实测硼浓度Yt，计算值为Yc，设定不同控制棒位时，相同的临界硼浓度差，其修正值相等，从而得出下关系：(Xt-Yt)＝(Xc-Yc)；
步骤二：采用插值方法得到温度对等温温度系数测量值的修正量
在具体的零功率物理实验中，对温度的修正不用与棒位及硼浓度的修正耦合，进行单独修正；
步骤三：结合步骤一和步骤二，综合堆芯温度与控制棒棒位-临界硼浓度对等温温度系数测量值的修正量，得到修正后的等温温度系数测量值；
再通过多次的测量与修正取平均值减小误差；
本修正方法总结为关系式：
ITC测量修正值＝试验状态下的棒位下ITC测量值+(控制棒棒位-临界硼浓度修正量)+温度修正量。
2.如权利要求1所述的一种零功率物理试验等温温度系数测量值修正方法，其特征在于：A＝225，B＝184。
CN201611052690.0A 2016-11-25 2016-11-25 一种零功率物理试验等温温度系数测量值修正方法 Active CN106782709B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611052690.0A CN106782709B (zh) 2016-11-25 2016-11-25 一种零功率物理试验等温温度系数测量值修正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611052690.0A CN106782709B (zh) 2016-11-25 2016-11-25 一种零功率物理试验等温温度系数测量值修正方法

Publications (2)

Publication Number Publication Date
CN106782709A true CN106782709A (zh) 2017-05-31
CN106782709B CN106782709B (zh) 2019-04-16

Family

ID=58912391

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611052690.0A Active CN106782709B (zh) 2016-11-25 2016-11-25 一种零功率物理试验等温温度系数测量值修正方法

Country Status (1)

Country Link
CN (1) CN106782709B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107221365A (zh) * 2017-06-26 2017-09-29 中广核研究院有限公司 一种商用压水堆物理试验系统
CN107887043A (zh) * 2017-10-10 2018-04-06 福建福清核电有限公司 一种反应堆功率系数测量方法
CN109727690A (zh) * 2019-03-18 2019-05-07 中国原子能科学研究院 双堆芯零功率实验装置仪控系统
WO2023226394A1 (zh) * 2022-05-25 2023-11-30 西安热工研究院有限公司 一种高温气冷堆堆芯等温温度系数测量方法
CN117198573A (zh) * 2023-11-07 2023-12-08 福建福清核电有限公司 堆芯次临界度的修正方法、堆芯次临界刻棒方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103594131A (zh) * 2012-08-15 2014-02-19 广东核电合营有限公司 百万千瓦级压水堆核电站慢化剂温度系数测量方法及装置
CN105336382A (zh) * 2014-08-15 2016-02-17 中国广核集团有限公司 核电站动态刻棒现场实施方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103594131A (zh) * 2012-08-15 2014-02-19 广东核电合营有限公司 百万千瓦级压水堆核电站慢化剂温度系数测量方法及装置
CN105336382A (zh) * 2014-08-15 2016-02-17 中国广核集团有限公司 核电站动态刻棒现场实施方法

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107221365A (zh) * 2017-06-26 2017-09-29 中广核研究院有限公司 一种商用压水堆物理试验系统
CN107887043A (zh) * 2017-10-10 2018-04-06 福建福清核电有限公司 一种反应堆功率系数测量方法
CN107887043B (zh) * 2017-10-10 2019-09-17 福建福清核电有限公司 一种反应堆功率系数测量方法
CN109727690A (zh) * 2019-03-18 2019-05-07 中国原子能科学研究院 双堆芯零功率实验装置仪控系统
WO2023226394A1 (zh) * 2022-05-25 2023-11-30 西安热工研究院有限公司 一种高温气冷堆堆芯等温温度系数测量方法
CN117198573A (zh) * 2023-11-07 2023-12-08 福建福清核电有限公司 堆芯次临界度的修正方法、堆芯次临界刻棒方法及装置
CN117198573B (zh) * 2023-11-07 2024-03-01 福建福清核电有限公司 堆芯次临界度的修正方法、堆芯次临界刻棒方法及装置

Also Published As

Publication number Publication date
CN106782709B (zh) 2019-04-16

Similar Documents

Publication Publication Date Title
CN106782709B (zh) 一种零功率物理试验等温温度系数测量值修正方法
CN107066745B (zh) 获取快中子堆堆芯瞬态过程三维中子通量密度分布的方法
Chang et al. Generalized measurable ignition criterion for inertial confinement fusion
CN111446014B (zh) 利用源量程探测器信号作为信号源的次临界刻棒方法
CN109859867B (zh) 一种基于本征正交分解的堆芯三维中子通量实时重构方法
CN105895174B (zh) 一种压水堆控制棒价值计算方法
CN104122071A (zh) 一种聚变堆增殖包层球床吹氚氦气热工水力特性实验装置及方法
CN107910083A (zh) 一种核电厂控制棒价值测量方法
CN104091036A (zh) 一种自然循环蒸汽发生器的热传导建模与计算方法
CN105895175A (zh) 一种标定核反应堆堆芯出口热电偶的方法
CN103335538B (zh) 一种电站凝汽器压力和传热端差的计算方法
CN107766641A (zh) 一种计算uo2燃料裂变气体热释放率的方法
CN114783637B (zh) 一种压水堆多燃料循环动态刻棒空间修正因子计算方法
CN107103189B (zh) 一种进行反应堆临界曲率搜索的方法
CN109119177A (zh) 一种压水堆稳压器压力和水位的耦合强度的确定方法、装置
CN109741840A (zh) 一种等温温度系数测量的优化方法
CN107122564A (zh) 一种压水堆pin‑by‑pin计算少群常数参数化方法
CN107423542A (zh) 一种适用于逐棒计算的非均匀泄漏修正方法
CN114242280B (zh) 反应堆反应性测量方法、装置、设备及存储介质
CN106683723B (zh) 一种反应堆钐毒在线测量方法
CN112632873B (zh) 一种核反应堆热腿测温热电偶测点位置布置方式确定方法
JP2014163803A (ja) 反応度温度係数推定装置及び方法
CN106610917B (zh) 模拟启堆时堆芯流道内冷却剂温度变化的新型计算方法
RU2786517C1 (ru) Способ определения коэффициентов реактивности по температуре топлива и плотности теплоносителя для области малых уровней мощности для ядерных реакторов большой мощности типа водо-водяных энергетических реакторов
JPS6211317B2 (zh)

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant