JPH0219919B2 - - Google Patents
Info
- Publication number
- JPH0219919B2 JPH0219919B2 JP56073574A JP7357481A JPH0219919B2 JP H0219919 B2 JPH0219919 B2 JP H0219919B2 JP 56073574 A JP56073574 A JP 56073574A JP 7357481 A JP7357481 A JP 7357481A JP H0219919 B2 JPH0219919 B2 JP H0219919B2
- Authority
- JP
- Japan
- Prior art keywords
- axial
- core
- output
- power
- xenon
- 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 - Lifetime
Links
- 229910052724 xenon Inorganic materials 0.000 claims description 25
- FHNFHKCVQCLJFQ-UHFFFAOYSA-N xenon atom Chemical compound [Xe] FHNFHKCVQCLJFQ-UHFFFAOYSA-N 0.000 claims description 25
- 238000000034 method Methods 0.000 claims description 14
- 238000012544 monitoring process Methods 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 239000000284 extract Substances 0.000 claims 1
- 238000004364 calculation method Methods 0.000 description 34
- 238000009792 diffusion process Methods 0.000 description 13
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 8
- 229910052740 iodine Inorganic materials 0.000 description 8
- 239000011630 iodine Substances 0.000 description 8
- 238000010521 absorption reaction Methods 0.000 description 6
- 230000004992 fission Effects 0.000 description 5
- 230000009257 reactivity Effects 0.000 description 5
- 230000001052 transient effect Effects 0.000 description 5
- 238000001514 detection method Methods 0.000 description 3
- 238000012806 monitoring device Methods 0.000 description 3
- 230000004907 flux Effects 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 239000002826 coolant Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- PNDPGZBMCMUPRI-UHFFFAOYSA-N iodine Chemical compound II PNDPGZBMCMUPRI-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- 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
- Monitoring And Testing Of Nuclear Reactors (AREA)
Description
【発明の詳細な説明】
この発明は、原子炉炉心の軸方向出力を炉心外
で監視する方法に関するものである。
で監視する方法に関するものである。
原子炉を安全に運転するためには原子炉炉心の
出力分布を常時監視することが必要である。
出力分布を常時監視することが必要である。
従来、このための原子炉炉心の軸方向出力分布
の把握は、その一例として、１次元拡散計算コー
ドによる計算によつてなされていた。原子炉の各
種定常条件下における出力分布は、原子炉出力、
冷却材温度、燃焼度、制御棒位置を決定すれば、
１次元拡散計算によつて精度よく把握されること
が確認されている。また、出力変動状態のような
キセノン過渡変動時には、よう素およびキセノン
の時間的変化を計算によつて追跡することを１次
元拡散計算コードに組み込んでおけば、やはり精
度よく出力分布が把握できることが知られてい
る。しかし、この１次元拡散計算コードによる把
握のための計算はくり返し計算を多用した、かな
り複雑なもので、多くのインプツト・データが必
要であり、計算機の計算時間も多いという問題点
もあるので原子力発電所の運転に直接利用するに
は困難であつた。そこで、原子力発電所サイトで
簡単に原子炉の出力分布を１次元拡散計算コード
により計算するために、軸方向炉心の計算の分点
数を少なくするなど多少精度を落して処理時間を
短くすることも考えられるが、たとえそのように
してサイトの計算機で計算できるようにしても、
各種の出力変動時の予測計算に即時に使うにはな
お時間がかかりすぎて不適であつた。このような
ことから、簡単かつ手較で、計算に多くの時間を
かけてることなく、各種出力変動時などの出力分
布変化を把握できる軸方向出力分布監視方法の開
発が望まれている。
の把握は、その一例として、１次元拡散計算コー
ドによる計算によつてなされていた。原子炉の各
種定常条件下における出力分布は、原子炉出力、
冷却材温度、燃焼度、制御棒位置を決定すれば、
１次元拡散計算によつて精度よく把握されること
が確認されている。また、出力変動状態のような
キセノン過渡変動時には、よう素およびキセノン
の時間的変化を計算によつて追跡することを１次
元拡散計算コードに組み込んでおけば、やはり精
度よく出力分布が把握できることが知られてい
る。しかし、この１次元拡散計算コードによる把
握のための計算はくり返し計算を多用した、かな
り複雑なもので、多くのインプツト・データが必
要であり、計算機の計算時間も多いという問題点
もあるので原子力発電所の運転に直接利用するに
は困難であつた。そこで、原子力発電所サイトで
簡単に原子炉の出力分布を１次元拡散計算コード
により計算するために、軸方向炉心の計算の分点
数を少なくするなど多少精度を落して処理時間を
短くすることも考えられるが、たとえそのように
してサイトの計算機で計算できるようにしても、
各種の出力変動時の予測計算に即時に使うにはな
お時間がかかりすぎて不適であつた。このような
ことから、簡単かつ手較で、計算に多くの時間を
かけてることなく、各種出力変動時などの出力分
布変化を把握できる軸方向出力分布監視方法の開
発が望まれている。
この発明は上記の如き事情に鑑みてなされたも
のであつて、サイトの運転員が出力変動時等に原
子炉の出力分布変化を簡単かつ迅速で精度も良く
把握することの出来る軸方向出力分布監視方法を
提供することを目的とするものである。
のであつて、サイトの運転員が出力変動時等に原
子炉の出力分布変化を簡単かつ迅速で精度も良く
把握することの出来る軸方向出力分布監視方法を
提供することを目的とするものである。
この目的に対応して、この発明の軸方向出力分
布監視方法は、炉心の軸方向上半分の出力を検出
する出力領域中性子検出器と、炉心の軸方向下半
分の出力を検出する出力領域中性子検出器とを炉
心の外側に配設し、それぞれの前記出力領域中性
子検出器による出力検出に基づいて炉心上半分の
相対出力と下半分の相対出力の差を表わすアキシ
ヤルオフセツト信号を取り出し、かつ、炉心底部
の軸方向位置をZb、炉心の下半分における出力
ピークの軸方向位置をＺ１、炉心の上半分におけ
る出力ピークの軸方向位置をＺ２、炉心頂部の軸
方向位置をZtとするとき、前記炉心を軸方向に
Zb≦Ｚ≦Ｚ１、Ｚ１≦Ｚ≦Ｚ２及びＺ２≦Ｚ≦
Ztの三つの部分に分け、それぞれの部分について
前記アキシヤルオフセツトA.O.信号を使用して
関数形表示により軸方向出力分布を導いて監視す
る場合に、炉心の軸方向下半分及び炉心軸方向上
半分のキセノンによる反応差をΔKeff，ｇ１及び
ｇ２を定数とするとき、前記アキシヤルオフセツ
トの変化ΔA.O.を ΔA.O.＝g1・ΔKeff＋g2 として前記軸方向出力分布を導くことを特徴とし
ている。
布監視方法は、炉心の軸方向上半分の出力を検出
する出力領域中性子検出器と、炉心の軸方向下半
分の出力を検出する出力領域中性子検出器とを炉
心の外側に配設し、それぞれの前記出力領域中性
子検出器による出力検出に基づいて炉心上半分の
相対出力と下半分の相対出力の差を表わすアキシ
ヤルオフセツト信号を取り出し、かつ、炉心底部
の軸方向位置をZb、炉心の下半分における出力
ピークの軸方向位置をＺ１、炉心の上半分におけ
る出力ピークの軸方向位置をＺ２、炉心頂部の軸
方向位置をZtとするとき、前記炉心を軸方向に
Zb≦Ｚ≦Ｚ１、Ｚ１≦Ｚ≦Ｚ２及びＺ２≦Ｚ≦
Ztの三つの部分に分け、それぞれの部分について
前記アキシヤルオフセツトA.O.信号を使用して
関数形表示により軸方向出力分布を導いて監視す
る場合に、炉心の軸方向下半分及び炉心軸方向上
半分のキセノンによる反応差をΔKeff，ｇ１及び
ｇ２を定数とするとき、前記アキシヤルオフセツ
トの変化ΔA.O.を ΔA.O.＝g1・ΔKeff＋g2 として前記軸方向出力分布を導くことを特徴とし
ている。
以下この発明の詳細を一実施例について説明す
る。
る。
第１図において、１は炉心であつて、圧力容器
２内に配設されている。圧溶容器２の外側には炉
心１を囲んで、４個の出力領域中性子検出器３が
炉心１の軸に関して対称に配設されている。出力
領域中性子検出器３は第２図に示すように各々上
部出力領域検出器３ａならびに下部出力領域検出
器３ｂとによつて構成されており、上部出力領域
検出器３ａは炉心１の上半分に対応して配設され
て炉心１の上半分の出力を検出することができ、
また、下部出力領域検出器３ｂは炉心１の下半分
の出力を検出することができる。上部出力領域検
出器３ａの出力をqT、下部出力領域検出器３ｂの
出力をqBとすると、アキシヤルオフセツト（A.
O.）は、 A.O.＝qT−qB／qT＋qB×100（％） …(1) で定義される。この発明では、まず、このアキシ
ヤルオフセツトを取り出して軸方向出力分布を導
く。すなわち、第２図に示す如く、炉心の高さを
1.0とし、かつ、 Ｐ（ｚ）：軸方向出力 Ｚ：軸方向高さ（０〜1.0） Z1：下部ピーク位置 Z2：上部ピーク位置 Fz1：下部ピーク出力 Fz2：上部ピーク出力 Ai（ｉ＝１〜４）：軸方向各領域の出力分布の
表示式 とすると、炉心の軸方向に(1)０ＺZ1、(2)Z1
ＺZ2、及び(3)Z2Ｚ1.0の領域に３分割し、 (1) ０ＺZ1では、 Ｐ（ｚ）＝A1＝Fz1｛cos〔2π（Ｚ−Z1）×a1〕 ＋c1（Ｚ−Z1）｝ …(2) (2) Z1ＺZ2では Ｐ（ｚ）＝A2＋A3 A2＝Fz1e-(z-z 1 )b1×cos〔−2π（Ｚ −Z1）×a2〕 …(3) A3＝Fz2e(z-z 2 )b2×cos〔2π（Ｚ −Z2）×a3〕 …(4) (3) Z2Ｚ1.0では Ｐ（ｚ）＝A4＝Fz2｛cos〔−2π（Ｚ−Z2）×a4〕 ＋c4（Z2−Ｚ）｝ …(5) とする。ここで、a2，a3はa2＝a3＝１／４（Z2−Z1） である。
２内に配設されている。圧溶容器２の外側には炉
心１を囲んで、４個の出力領域中性子検出器３が
炉心１の軸に関して対称に配設されている。出力
領域中性子検出器３は第２図に示すように各々上
部出力領域検出器３ａならびに下部出力領域検出
器３ｂとによつて構成されており、上部出力領域
検出器３ａは炉心１の上半分に対応して配設され
て炉心１の上半分の出力を検出することができ、
また、下部出力領域検出器３ｂは炉心１の下半分
の出力を検出することができる。上部出力領域検
出器３ａの出力をqT、下部出力領域検出器３ｂの
出力をqBとすると、アキシヤルオフセツト（A.
O.）は、 A.O.＝qT−qB／qT＋qB×100（％） …(1) で定義される。この発明では、まず、このアキシ
ヤルオフセツトを取り出して軸方向出力分布を導
く。すなわち、第２図に示す如く、炉心の高さを
1.0とし、かつ、 Ｐ（ｚ）：軸方向出力 Ｚ：軸方向高さ（０〜1.0） Z1：下部ピーク位置 Z2：上部ピーク位置 Fz1：下部ピーク出力 Fz2：上部ピーク出力 Ai（ｉ＝１〜４）：軸方向各領域の出力分布の
表示式 とすると、炉心の軸方向に(1)０ＺZ1、(2)Z1
ＺZ2、及び(3)Z2Ｚ1.0の領域に３分割し、 (1) ０ＺZ1では、 Ｐ（ｚ）＝A1＝Fz1｛cos〔2π（Ｚ−Z1）×a1〕 ＋c1（Ｚ−Z1）｝ …(2) (2) Z1ＺZ2では Ｐ（ｚ）＝A2＋A3 A2＝Fz1e-(z-z 1 )b1×cos〔−2π（Ｚ −Z1）×a2〕 …(3) A3＝Fz2e(z-z 2 )b2×cos〔2π（Ｚ −Z2）×a3〕 …(4) (3) Z2Ｚ1.0では Ｐ（ｚ）＝A4＝Fz2｛cos〔−2π（Ｚ−Z2）×a4〕 ＋c4（Z2−Ｚ）｝ …(5) とする。ここで、a2，a3はa2＝a3＝１／４（Z2−Z1） である。
次に、
∫1.0 0Ｐ（ｚ）dz＝1.0 …(6)
と規格化する。(1)式のアキシヤルオフセツト
（A.O.）は次のように表される。
（A.O.）は次のように表される。
A.O.＝∫1.0 z2A4dz＋∫z2 0.5（A2＋A3）dz
−∫0.5 z1（A2＋A3）dz−∫z1 0A1dz …(7)
一方、軸方向各点における出力は第３図に例示
する如く、アキシヤルオフセツト変化に対して直
線となる傾向がある。したがつて、一般にA.O.
と上部ピークFz2は次のように表示できる。
する如く、アキシヤルオフセツト変化に対して直
線となる傾向がある。したがつて、一般にA.O.
と上部ピークFz2は次のように表示できる。
Fz1＝α1×（A.O.）＋α0
Fz2＝β1×（A.O.）＋β0 …(8)
先の(2)〜(5)までの式を(6)，(7)，(8)と組み合せる
と(2)〜(5)式中の定数a1，a4，b1，b2が決まつてく
る。
と(2)〜(5)式中の定数a1，a4，b1，b2が決まつてく
る。
一方、定数c1，c4は、炉心上下端での外挿距離
λ1，λ2が一定値であるとの仮定を使用して、次の
式で決められる。
λ1，λ2が一定値であるとの仮定を使用して、次の
式で決められる。
前記のように式(2)〜(5)の各定数は次のデータを
与えれば自動的に決定され、出力分布を合成計算
できることになる。
与えれば自動的に決定され、出力分布を合成計算
できることになる。
上部ピーク位置 ：Z2
下部ピーク位置 ：Z1
(7)式の係数 ：α1，α0，β1，β0
上端での外挿距離 ：λ2
下端での外挿距離 ：λ1
これらのデータを予め計算にて求めておけば、
任意のアキシヤルオフセツトに相当する軸方向出
力分布が(2)〜(4)式を使用して合成計算される。
任意のアキシヤルオフセツトに相当する軸方向出
力分布が(2)〜(4)式を使用して合成計算される。
以上の技術は昭和55年特許願第182107号の明細
書、図面に記載されているところである。
書、図面に記載されているところである。
ところで、こうした求められる軸方向出力分布
のうち、キセノン過渡変動時のものを一層精度を
良くするためには、先の軸方向出力分布にキセノ
ン、よう素の時間変化の影響を組込むことが必要
である。
のうち、キセノン過渡変動時のものを一層精度を
良くするためには、先の軸方向出力分布にキセノ
ン、よう素の時間変化の影響を組込むことが必要
である。
キセノン、よう素の時間的変化は下記の式の表
わされる。
わされる。
dI（ｔ，Ｚ）／dt＝λIＩ（ｔ，Ｚ）＋YI（
）…(10) dx（ｔ，Ｚ）／dt＝λIＩ（ｔ，Ｚ）＋Yx（
） −（λx＋（））・ｘ（ｔ，Ｚ）…(11
) ただし、 Ｉ（ｔ，Ｚ）：時間ｔにおける軸方向位置Ｚにお
けるよう素の個数密度 ｘ（ｔ，Ｚ）：時間ｔにおける軸方向位置Ｚにお
けるキセノンの個数密度 λI：よう素の崩壊定数 λx：キセノンの崩壊定数 YI：よう素の核分裂収率 Yx：キセノンの核分裂収率 （）：軸方向位置Ｚにおける核分裂数で次
のように計算される。
）…(10) dx（ｔ，Ｚ）／dt＝λIＩ（ｔ，Ｚ）＋Yx（
） −（λx＋（））・ｘ（ｔ，Ｚ）…(11
) ただし、 Ｉ（ｔ，Ｚ）：時間ｔにおける軸方向位置Ｚにお
けるよう素の個数密度 ｘ（ｔ，Ｚ）：時間ｔにおける軸方向位置Ｚにお
けるキセノンの個数密度 λI：よう素の崩壊定数 λx：キセノンの崩壊定数 YI：よう素の核分裂収率 Yx：キセノンの核分裂収率 （）：軸方向位置Ｚにおける核分裂数で次
のように計算される。
（）＝（Σf1φ1＋Σf2φ2）
×Ｐ（ｚ）×PREL …(12)
x（）：軸方向位置Ｚにおけるキセノンの
吸収断面積と中性子束の積で次のように計算
される。x （）＝（σa1 xφ1＋σa2 xφ2） ×Ｐ（ｚ）×PREL …(13) Σf1，Σf2：高速組、熱組の炉心平均核分裂断面
積 φ1，φ2：高速組、熱組の炉心平均中性子束 σa1 x，σa2 x：高速組、熱組のキセノンのミクロ
吸収断面積 PREL：原子炉の炉心相対出力 Ｐ（ｚ）：軸方向位置Ｚにおける出力で平均が
1.0に規格化されたもの。前述のＰ（ｚ）と同
じもの。
吸収断面積と中性子束の積で次のように計算
される。x （）＝（σa1 xφ1＋σa2 xφ2） ×Ｐ（ｚ）×PREL …(13) Σf1，Σf2：高速組、熱組の炉心平均核分裂断面
積 φ1，φ2：高速組、熱組の炉心平均中性子束 σa1 x，σa2 x：高速組、熱組のキセノンのミクロ
吸収断面積 PREL：原子炉の炉心相対出力 Ｐ（ｚ）：軸方向位置Ｚにおける出力で平均が
1.0に規格化されたもの。前述のＰ（ｚ）と同
じもの。
である。
(10)式と(11)式を解くと次のようになる。
ただし
I0（ｚ）：軸方向位置Ｚにおけるよう素の個数密
度の初期値 X0（ｚ）：軸方向位置Ｚにおけるキセノンの個
数密度の初期値 ΔT：経過時間 である。
度の初期値 X0（ｚ）：軸方向位置Ｚにおけるキセノンの個
数密度の初期値 ΔT：経過時間 である。
Ｐ（ｚ）に前述の(2)〜(5)式の関数形による表示
法を使用すればアキシヤルオフセツトから任意の
時点のキセノン軸方向分布が求まることになる
（なお、よう素については中性子吸収断面積が小
さいのでこれを無視する。）。キセノン分布が求ま
れば、炉心上半分および炉心下半分のキセノンに
よる反応度差を次のように算出できる。
法を使用すればアキシヤルオフセツトから任意の
時点のキセノン軸方向分布が求まることになる
（なお、よう素については中性子吸収断面積が小
さいのでこれを無視する。）。キセノン分布が求ま
れば、炉心上半分および炉心下半分のキセノンに
よる反応度差を次のように算出できる。
ΔKeffx＝lnKeffT／KeffB×100 …(17)
ただし、
D1，Σa1，ΣR，νΣf1：高速組の炉心平均核定数
Ｄ；拡散係数
Σa；吸収断面積
ΣR；除去断面積
νΣf；核分裂断面積と中性子発生数の積
D2，Σa2，νΣf2：熱組の炉心平均核定数
B2：バツクリング
Σax 1T，Σax 2T：炉心半分のキセノンの高速組、
熱組の吸収断面積 Σax 1B，Σax 2B：炉心下半分のキセノンの高速
組、熱組の吸収断面積 である。ここで反応度差ΔKeffxとアキシヤルオ
フセツト変化との関係を計算し、グラフにすると
第４図に示すようになり、両者の関係は良好な直
線状となる。すなわち、炉心上半分と下半分の反
応度差とアキシヤルオフセツト変化とには次のよ
うな関係があるとすることができる。
熱組の吸収断面積 Σax 1B，Σax 2B：炉心下半分のキセノンの高速
組、熱組の吸収断面積 である。ここで反応度差ΔKeffxとアキシヤルオ
フセツト変化との関係を計算し、グラフにすると
第４図に示すようになり、両者の関係は良好な直
線状となる。すなわち、炉心上半分と下半分の反
応度差とアキシヤルオフセツト変化とには次のよ
うな関係があるとすることができる。
ΔA.O.＝g1・ΔKeff＋g2 …(18)
ただしg1，g2は定数である。したがつて、炉心
上半分と下半分のキセノンによる反応度差にわか
れば、（18）式よりアキシヤルオフセツト変化幅
がわかるので新しい出力分布が計算できる。キセ
ノン過渡動時のアキシヤルオフセツト変化は、こ
のようにして時間経過にしたがい順次計算をして
行うことにより把握できる。出力レベル変化によ
つてアキシヤルオフセツト変化および制御棒挿入
によるアキシヤルオフセツト変化は予め計算した
ものを表形式でインプツトしておき利用すれば、
出力変動時などのキセノン過渡状態の出力分布が
上記のように簡単に把握できる。
上半分と下半分のキセノンによる反応度差にわか
れば、（18）式よりアキシヤルオフセツト変化幅
がわかるので新しい出力分布が計算できる。キセ
ノン過渡動時のアキシヤルオフセツト変化は、こ
のようにして時間経過にしたがい順次計算をして
行うことにより把握できる。出力レベル変化によ
つてアキシヤルオフセツト変化および制御棒挿入
によるアキシヤルオフセツト変化は予め計算した
ものを表形式でインプツトしておき利用すれば、
出力変動時などのキセノン過渡状態の出力分布が
上記のように簡単に把握できる。
以上の如き、操作及び演算による軸方向出力分
布の検出は、第５図に示す装置によつて行う。す
なわち、第５図において、１１は軸方向出力分布
監視装置であり、軸方向出力分布監視装置１１は
上部出力領域中性子検出器３ａ及び下部出力領域
中性子検出器３ｂを備え、それぞれの検出器３
ａ，３ｂはアキシヤルオフセツト演算装置１３に
接続し、アキシヤルオフセツト演算装置１３の出
力信号は、軸方向出力分布演算装置１４に接続す
る。検出器３ａ，３ｂはそれぞれ炉心の上半分出
力及び下半分出力を検出し、この信号をアキシヤ
ルオフセツト演算装置１３に送り、ここで(1)式に
示す初期アキシヤルオフセツト信号が作成され、
これが軸方向出力分布演算装置に送られる。軸方
向出力分布演算装置１４では、入力したアキシヤ
ルオフセツト信号と、及び別に入力された出力レ
ベル変化、制御棒位置の信号及び時間幅ΔTによ
つて、かつ、予め準備されて入力されたZ1，Z2，
α1，α0，β1，β0，λ1，λ2等のＰ（ｚ）計算用デー
タ、キセノン、よう素のデータ、２群核定数、ア
キシヤルオフセツト変化のデータ等によつて、
（18）式による演算を行い、軸方向出力分布変化
を演算する。この合成計算された軸方向出力分布
は、必要に応じて他の機器への入力信号とするこ
とができる。
布の検出は、第５図に示す装置によつて行う。す
なわち、第５図において、１１は軸方向出力分布
監視装置であり、軸方向出力分布監視装置１１は
上部出力領域中性子検出器３ａ及び下部出力領域
中性子検出器３ｂを備え、それぞれの検出器３
ａ，３ｂはアキシヤルオフセツト演算装置１３に
接続し、アキシヤルオフセツト演算装置１３の出
力信号は、軸方向出力分布演算装置１４に接続す
る。検出器３ａ，３ｂはそれぞれ炉心の上半分出
力及び下半分出力を検出し、この信号をアキシヤ
ルオフセツト演算装置１３に送り、ここで(1)式に
示す初期アキシヤルオフセツト信号が作成され、
これが軸方向出力分布演算装置に送られる。軸方
向出力分布演算装置１４では、入力したアキシヤ
ルオフセツト信号と、及び別に入力された出力レ
ベル変化、制御棒位置の信号及び時間幅ΔTによ
つて、かつ、予め準備されて入力されたZ1，Z2，
α1，α0，β1，β0，λ1，λ2等のＰ（ｚ）計算用デー
タ、キセノン、よう素のデータ、２群核定数、ア
キシヤルオフセツト変化のデータ等によつて、
（18）式による演算を行い、軸方向出力分布変化
を演算する。この合成計算された軸方向出力分布
は、必要に応じて他の機器への入力信号とするこ
とができる。
このよう構成された軸方向出力分布監視方法で
は、出力変動時などのキセノン過渡状態の軸方向
出力分布を精度よく、簡単に導くことができる。
すなわち、第６図に全出力状態から制御棒バンク
Ｄの挿入および引抜により外乱を作つた後のキセ
ノン振動時のアキシヤルオフセツト変化を１次元
拡散計算とこの発明の方法による予測値と比較し
て示している。この発明の方法によつて導いた出
力値は一次元拡散計算結果と良く一致しているこ
とがわかる。また、第７図には全出力状態から50
％出力状態に出力低下後のキセノン振動時のアキ
シヤルオフセツトをこの発明の方法による予測値
を１次元拡散計算法による計算結果と比較してい
る。第７図からこの場合にもこの発明による予測
値は一次元拡散計算結果と良く一致していること
がわかる。しかも、この発明の方法によれば、従
来の一次元拡散計算コードによる計算時間の約1/
１００程度で済んでおり、サイト運転員が利用する
方法として適切である。
は、出力変動時などのキセノン過渡状態の軸方向
出力分布を精度よく、簡単に導くことができる。
すなわち、第６図に全出力状態から制御棒バンク
Ｄの挿入および引抜により外乱を作つた後のキセ
ノン振動時のアキシヤルオフセツト変化を１次元
拡散計算とこの発明の方法による予測値と比較し
て示している。この発明の方法によつて導いた出
力値は一次元拡散計算結果と良く一致しているこ
とがわかる。また、第７図には全出力状態から50
％出力状態に出力低下後のキセノン振動時のアキ
シヤルオフセツトをこの発明の方法による予測値
を１次元拡散計算法による計算結果と比較してい
る。第７図からこの場合にもこの発明による予測
値は一次元拡散計算結果と良く一致していること
がわかる。しかも、この発明の方法によれば、従
来の一次元拡散計算コードによる計算時間の約1/
１００程度で済んでおり、サイト運転員が利用する
方法として適切である。
以上の説明から明らかな通り、この発明によれ
ば、サイトで運転員が出力変動時に原子炉の出力
分布変化を簡単かつ迅速で精度良く把握すること
ができる軸方向出力分布監視方法を得ることがで
きる。
ば、サイトで運転員が出力変動時に原子炉の出力
分布変化を簡単かつ迅速で精度良く把握すること
ができる軸方向出力分布監視方法を得ることがで
きる。
第１図は中性子検出器の配置を示す平面説明
図、第２図は中性子検出器の配置位置を示す側面
説明図、第３図は炉心軸方向各部出力とアキシヤ
ルオフセツト信号との関係を示すグラフ、第４図
は反応度差によるアキシヤルオフセツトの変化例
を示すグラフ、第５図は軸方向出力分布検出装置
の構成を示すブロツク図、第６図は一次元拡散計
算で求めたアキシヤルオフセツト変化とこの発明
の方法で求めたアキシヤルオフセツト変化との比
較を示すグラフ、及び第７図は一次元拡散計算で
求めたアキシヤルオフセツト変化とこの発明の方
法で求めたアキシヤルオフセツト変化との比較の
他の例を示すグラフである。 ３（３ａ，３ｂ）……出力領域中性子検出器、
４……中間領域中性子検出器、１１……軸方向出
力分布監視装置、１３……アキシヤルオフセツト
演算装置、１４……軸方向出力分布演算装置。
図、第２図は中性子検出器の配置位置を示す側面
説明図、第３図は炉心軸方向各部出力とアキシヤ
ルオフセツト信号との関係を示すグラフ、第４図
は反応度差によるアキシヤルオフセツトの変化例
を示すグラフ、第５図は軸方向出力分布検出装置
の構成を示すブロツク図、第６図は一次元拡散計
算で求めたアキシヤルオフセツト変化とこの発明
の方法で求めたアキシヤルオフセツト変化との比
較を示すグラフ、及び第７図は一次元拡散計算で
求めたアキシヤルオフセツト変化とこの発明の方
法で求めたアキシヤルオフセツト変化との比較の
他の例を示すグラフである。 ３（３ａ，３ｂ）……出力領域中性子検出器、
４……中間領域中性子検出器、１１……軸方向出
力分布監視装置、１３……アキシヤルオフセツト
演算装置、１４……軸方向出力分布演算装置。
Claims (1)
- 【特許請求の範囲】 １ 炉心の軸方向上半分の出力を検出する出力領
域中性子検出器と、炉心の軸方向下半分の出力を
検出する出力領域中性子検出器とを炉心の外側に
配設し、それぞれの前記出力領域中性子検出器に
よる出力検出に基づいて炉心上半分の相対出力と
下半分の相対出力の差を表わすアキシヤルオフセ
ツト信号を取り出し、かつ、炉心底部の軸方向位
置をZb、炉心の下半分における出力ピークの軸
方向位置をＺ１、炉心の上半分における出力ピー
クの軸方向位置をＺ２、炉心頂部の軸方向位置を
Ztとするとき、前記炉心を軸方向にZb≦Ｚ≦Ｚ
１、Ｚ１≦Ｚ≦Ｚ２及びＺ２≦Ｚ≦Ztの三つの部
分に分け、それぞれの部分について前記アキシヤ
ルオフセツトA.O.信号を使用して関数形表示に
より軸方向出力分布を導いて監視する場合に、炉
心の軸方向下半分及び炉心軸方向上半分のキセノ
ンによる反応差をΔKeff，ｇ１及びｇ２を定数と
するとき、前記アキシヤルオフセツトの変化ΔA.
O.を ΔA.O.＝g1・ΔKeff＋g2 として前記軸方向出力分布を導くことを特徴とす
る軸方向出力分布監視方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP56073574A JPS57189100A (en) | 1981-05-18 | 1981-05-18 | Method of monitoring axial power distribution |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP56073574A JPS57189100A (en) | 1981-05-18 | 1981-05-18 | Method of monitoring axial power distribution |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS57189100A JPS57189100A (en) | 1982-11-20 |
| JPH0219919B2 true JPH0219919B2 (ja) | 1990-05-07 |
Family
ID=13522183
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP56073574A Granted JPS57189100A (en) | 1981-05-18 | 1981-05-18 | Method of monitoring axial power distribution |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS57189100A (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0662218U (ja) * | 1993-02-09 | 1994-09-02 | 株式会社ベンリー | ねじ付釘 |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6283696A (ja) * | 1985-10-09 | 1987-04-17 | 三菱原子力工業株式会社 | 炉内軸方向出力分布の監視方法 |
| JP2008175692A (ja) * | 2007-01-18 | 2008-07-31 | Nuclear Fuel Ind Ltd | 炉心の軸方向出力分布の測定方法 |
-
1981
- 1981-05-18 JP JP56073574A patent/JPS57189100A/ja active Granted
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0662218U (ja) * | 1993-02-09 | 1994-09-02 | 株式会社ベンリー | ねじ付釘 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPS57189100A (en) | 1982-11-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4999222B2 (ja) | 原子炉の炉心の少なくとも１つの動作パラメタを監視する方法 | |
| Elias et al. | The reflooding phase of LOCA in PWRs, Part II, Rewetting and liquid entrainment | |
| KR970004354B1 (ko) | 노심반응도 검증장치 및 방법 | |
| SE513236C3 (sv) | System för effektövervakning i en kärnreaktor | |
| SE513236C2 (sv) | System för effektövervakning i en kärnreaktor | |
| JPH0219919B2 (ja) | ||
| Atta et al. | Determination of void fraction profile in a boiling water reactor channel using neutron noise analysis | |
| Mogensen | Determination of fission gas yields from isotope ratios | |
| JPH0219918B2 (ja) | ||
| Nijsing et al. | The hybrid method, a new accurate computation method for transient temperature distributions in fuel rods | |
| JPS6337916B2 (ja) | ||
| Katona | Investigation of the noise source characteristics related to the theory of subcooled boiling neutron noise | |
| JP2975654B2 (ja) | 炉心監視装置 | |
| Raj et al. | Analysis of conduction controlled rewetting of hot surfaces based on two-region model | |
| JPS58124987A (ja) | 原子力発電所の制御棒操作ガイド装置 | |
| JP2008175561A (ja) | 沸騰水型原子炉の炉心流量計測装置および方法 | |
| JPS59164987A (ja) | 高速炉の反応度監視装置 | |
| JPS59195192A (ja) | 原子炉炉心状態監視装置 | |
| JPS6132639B2 (ja) | ||
| JP2895082B2 (ja) | 原子炉におけるキセノン振動の制御支援方法 | |
| Hodde | FUEL CYCLE PROGRAM. A BOILING WATER REACTOR RESEARCH AND DEVELOPMENT PROGRAM. Quarterly Report No. 5, July 1961-September 1961 | |
| Tzanos | A semianalytic method for the solution of the steady-state steam generator equations | |
| JPS6243596A (ja) | 原子炉の臨界接近監視方法 | |
| JPS5951391A (ja) | 原子炉炉心状態監視装置 | |
| Estes et al. | Sandia Pulsed Reactor Facility (SPRF) calculator-assisted pulse analysis and display system |