JP4901631B2 - ドップラー反応度係数の測定方法 - Google Patents
ドップラー反応度係数の測定方法 Download PDFInfo
- Publication number
- JP4901631B2 JP4901631B2 JP2007197130A JP2007197130A JP4901631B2 JP 4901631 B2 JP4901631 B2 JP 4901631B2 JP 2007197130 A JP2007197130 A JP 2007197130A JP 2007197130 A JP2007197130 A JP 2007197130A JP 4901631 B2 JP4901631 B2 JP 4901631B2
- Authority
- JP
- Japan
- Prior art keywords
- time
- series data
- reactivity
- temperature
- reactor
- 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.)
- Active
Links
Images
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
- G21C17/104—Measuring reactivity
-
- G—PHYSICS
- G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
- G21D—NUCLEAR POWER PLANT
- G21D3/00—Control of nuclear power plant
- G21D3/001—Computer implemented control
-
- 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)
Description
さらに、各運転サイクルの間に行なわれる定期検査において、次のサイクルを運転する炉心の炉物理的な特性を測定し、評価するために原子炉特性調査工事(起動時炉物理試験)が実施されている。この工事により、例えば設計された炉心が所定の操作で臨界に達するか否か、炉心に対する制御棒動作時の反応度の変化や減速材温度が変化したときの反応度変化等が測定され、炉心設計の妥当性の確認がなされている。
この方法は、基本的には等温温度反応度係数測定法と動特性同定法を組合せて燃料温度の反応度係数の測定(炉出力が低い状態であり、ドップラー反応度係数の測定となる)を行うものであり、凡そ以下の手順(ステップ、工程)による。
なお、ここに「等温温度反応度係数」とは、燃料温度のみの(燃料温度で偏微分した)反応度係数と減速材温度のみの(減速材温度で偏微分した)反応度係数の和をいう。
次いで、制御棒を引抜くことにより外部反応度を添加し、炉出力を所定量、例えば1%程度増加させる。
この際、外部添加反応度、中性子炉外検出器応答、冷却材(減速材)入口温度、冷却材(減速材)平均温度の時間の経過による変化を測定し、これらについての時系列データを採取する。
また、中性子炉外検出器の応答をディジタル反応度計に入力して原子炉の反応度の時間応答を求める。
具体的には、数値フーリエ変換し、周波数応答伝達関数に代入し、この応答関係を充たすドップラー反応度係数を最小自乗フィッティング法(基本的には最小自乗法)により求める。
なお、減速材温度反応度係数は、等温温度反応度係数からドップラー反応度係数を差引く事で求められる。
また、冷却材(減速材)入口温度を測定する必要があるが、現在の炉物理試験時の測定項目には無いためふさわしくない。
なおここに、炉周期とは、原子炉の出力がe(約2.718)倍となる時間のことである。
未臨界、従って炉の出力が極めて低い状態では、原子炉の自己制御が作用しない。
一定炉周期反応度は、未臨界の状態における反応度に無関係にほぼ一定である。なおこれについては、実測とシミュレーション解析が最も一致する組合せを調べることにより確認した。
さらに、反応度フィードバックの寄与を読取る手段として、ρinとρpを、フィードバックが効きにくい低出力領域でシミュレーション解析し、実測したNIS信号を再現させる値を探索して決定する。
測定の前提として、その対象となる原子炉は、出力領域の炉外中性子束検出器、冷却水平均温度等が時系列データとして測定可能であり、炉心特性解析を行うのに必要なデータは揃っているものとする。従って、例えば一点炉動特性パラメータβi、λi(ここに、iは6群の遅発中性子。i=1、・・・、6)、燃料棒温度のインポータンス出力加重平均補正係数{一点炉動特性モデルから算出される燃料棒の体積加重平均値から中性子束分布と随伴中性子束(中性子インポータンス)分布で加重された平均温度を換算する係数}、制御棒バンク操作履歴等は、別途の理論解析等からかなり正確な値または凡その値が判明していたり、運転記録から取得可能であったりする。さらに、本測定に先立って実施される等温温度係数測定試験から、等温温度反応度係数(=ドップラー反応度係数+減速材温度反応度係数)も既知である。
以上の他、測定対象である燃料のドップラー反応度係数を含む各種の係数、初期の未臨界度等も、凡その値は設計時の理論解析、従来からの経験等で判明している。またこのため、誤差関数を基本的には試行錯誤法で評価する際にも、これらの凡その値を初期値として入力したり、その近辺の値を入力したりすることが多くなる。
未臨界の状態の原子炉(PWR)から制御棒を所定量引抜いて臨界かつ僅かな出力状態とし、この際の中性子束と減速材平均温度の一連の変化を時系列データとして連続採取する。なお、この際の未臨界度ρsubは、制御棒を所定量引抜くことにより与えた印加反応度ρinと一定炉周期反応度ρpから逆算可能である。
また、商用炉においては、未臨界時においても多少の出力がある。
さらに、臨界状態とは、核分裂を原因として炉内で発生する中性子数と炉内での吸収と炉外への漏れを原因として消滅する中性子数が等しい(実効増倍率が1)平衡状態を指し、炉内の(熱)出力はこの平衡状態がどの中性子数レベルで実現されるかによって決まる。
前記の如く、炉の出力が定格時に比べて極めて小さく、このため中性子束が小さい状態の原子炉における中性子束の変化を時系列データとして採取する場合には、そのままでは中性子測定器からの制約が生じる。即ち、炉内には使用済み核燃料が在るが、一般的に使用されているNIS(中性子検出器)は、この使用済み核燃料に由来するγ線に対しても反応する。そして、この使用済み燃料に由来するγ線の線量は、零出力試験時においてもほぼ一定のレベルで存在する。一方、零出力試験時に炉内で発生する中性子束は、炉出力が小さいため低レベルである。その結果、零出力試験時においては、採取したデータのバックグラウンド、あるいはγ線による雑音(ノイズ)を無視できないからである。なお、炉出力が多少大きくなれば、炉内で発生する中性子が増加するため、γ線の影響は小さくなり、無視してよいこととなる。
そこで、零出力試験時においては、電離箱を使用して得られた中性子のデータからγ線による影響、即ち誤って中性子として採取されたγ線による成分(雑音)を、γ線による成分は炉出力に対して一定という性質を用いて除洗する。但し、将来建設されるPWR等であって、中性子をγ線の影響なく直接測定することが可能となる措置がなされていれば、この操作は不必要となる。
なおここに、誤差関数が対数(ln)となっているのは、炉出力は時間に対して指数的に増大することを考慮したものである。
また、(1)式におけるパラメータおよび誤差関数の表記は、この限りではない。
なお、上記解析に必要な炉出力応答が反応度フィードバックの寄与が小さい低出力領域であること、また初期出力に対する炉出力応答の相対変化のみが必要であることから、解析ではこの条件を満足させる限り初期炉出力Ps 0は任意に設定してよい。炉出力の絶対値応答は後述する処理法で決定する。
γ線が除洗された出力データを実測のNIS信号から再構成すると、以下の(2)式のように求められる。
なお、前記のごとく商用炉では、未臨界の状態でも炉出力が多少あるため、(3)式の分母は0ではない。
γ線が除洗された炉出力の時系列データから、一点炉動特性方程式に対する逆動特性法より反応度ρの時系列データを求める。反応度フィードバック寄与成分Δρfdは、(5)式のように反応度変化ρ(t)からρinを差し引くことで求める。
しかし、NIS信号をγ線除洗した信号からは、初期出力と最大到達出力の比Rzmは決定できるものの、その出力の絶対値は決定できない。この場合、減速材炉心入口温度が一定であるか、あるいは測定されているならば、実測された減速材平均温度Tc,avから炉心入口温度と出口温度の差を評価することで、炉出力を求めることができる。しかし、炉出力変更時には、炉心側からの熱供給と冷却ループ蒸気発生器2次側からの除熱の一時的なアンバランス状態が生じるため、この影響が炉心入口温度の変化となって現れ、炉心入口温度一定の仮定を適用できない。
また、既に建造され、運転履歴がある原子炉の炉心入口に新しく温度センサーを取付け、減速材の炉心入口の温度を時系列データとして精度良く測定するためには、新たな計測装置の導入などが必要となる。
図1において、10は炉心であり、20は蒸気発生器であり、30は冷却水循環ポンプであり、41は原子炉出口側配管であり、42は原子炉入口側配管であり、矢印は減速材(冷却水)の流れを示し、太い白矢印は熱の流れを示す。
一次冷却閉ループでの冷却特性を決定づける最も重要なパラメータは、蒸気発生器の1次側から2次側への伝熱に関する時定数τsg,12であり、これを決定しなければならない。
材料試験炉、臨界実験装置等の小型の原子炉と異なり、発電用の大型の原子炉、例えばPWRは蒸気発生器の様な熱交換器を有しているため、中性子束データのピークの到達時間と冷却水等の減速材の温度のピークの到達時間とには、時間差が生じるのが普通である。この時間差と時定数τsg,12とに直接的な関係が存在すること、さらに到達上限出力と減速材平均温度Tc,avの最大到達温度に強い相関が見られることに着目して、到達上限出力とτsg,12を決定する。
この誤差関数E(τsg,12,P0)が最小となる条件を、即ち減速材の最大温度到達時刻と減速材の平均温度の最大上昇幅が共に実測値と等しくなる時定数τsg,12および初期出力P0の値を、前記(1)式の場合と同様に求める。到達上限炉出力Pmaxは、誤差関数Eが最小となる(τsg,12,P0)のときのシミュレーション解析から求められた最大到達炉出力とする。P0は(3)式から再評価する。
なお、(9)式のパラメータおよび誤差関数の表記はこの限りではない。
上記の方法で決定した最大到達炉出力Pmaxと出力変更幅Rzmを用いて、γ線雑音除洗処理をしたNIS信号から初期出力から最大到達炉出力までの炉出力応答を決定する。燃料棒平均温度に関する熱伝導方程式に、炉出力応答と実測した減速材平均温度Tc,avを代入することで、燃料棒平均温度変化ΔTf,av(t)を決定する。
燃料の平均温度は炉出力に応じて昇降し、さらにその温度変化は減速材に比較して大きくかつ速やかである。燃料の温度変化に対する反応度の応答を調べる一次摂動理論を導入した場合、平均燃料温度変化ΔTf,av(t)は、以下の(10)式に示す様なインポータンス出力加重平均ΔTl,ip f,av(t)として与えられる。
なおここで一次摂動理論を用いたのは、摂動理論は、微小な変化を与え、その変化の影響を考慮する理論であるため、摂動の無い基本式に補正量を適用する手法として優れており、さらに補正は補正項として一次の項を近似した一次摂動理論が一般的であり、炉物理試験時においては、燃料の温度変化が微小であるので、摂動理論がふさわしいと判断したことによる。
ドップラー反応度係数から定義される誤差関数を、(13)式に示す。この誤差関数が最小となるドップラー反応度係数αfが、測定されたドップラー反応度係数である。
未臨界または臨界になった状態の炉心に反応度を印加して炉出力を一定炉周期で所定量上昇させ、その際の中性子束を時系列データとして測定する中性子束時系列データ測定ステップと、
未臨界または臨界になった状態の炉心に反応度を印加して炉出力を一定炉周期で所定量上昇させ、その際の炉内の減速材の平均温度を所定の手順で時系列データとして取得する炉内減速材平均温度時系列データ取得ステップと、
測定された中性子束の時系列データから、一点炉動特性方程式に対する逆動特性法を用いて反応度の時系列データを取得する反応度時系列データ取得ステップと、
前記取得された炉内減速材平均温度の時系列データと中性子束の時系列データを基に、所定の手順で前記2つの時系列データに整合する炉出力の時系列データを取得する炉出力時系列データ取得ステップと、
取得した炉出力の時系列データと所定の動特性モデルを用いて取得された所定平均の燃料温度の時系列データを得る燃料温度時系列データ取得ステップと、
取得した反応度の時系列データと印加した反応度を用いて反応度フィードバック寄与成分の時系列データを求める反応度フィードバック寄与成分の時系列データ取得ステップと、
前記炉内の減速材平均温度の時系列データと、前記所定平均の燃料温度の時系列データと、等温温度反応度係数と、前記反応度フィードバック寄与成分の時系列データを用いて、所定の手順でドップラー反応度係数を求めるドップラー反応度係数取得ステップを有していることを特徴とするドップラー反応度係数の測定方法である。
また、「炉出力を一定炉周期で所定量上昇させ」の所定量とは、実際に原子炉を定格で運転している状態におけるドップラー反応度係数を測定するという面からは、出来るだけ定格出力に近い、即ち大きければ大きいほど好ましいが、原子炉特性調査工事中等の試験であること、計器や測定そのものから生じる制約等のため実際には1%以内である。
また、「時系列データ」とは、起動時から所定量の出力までの時間の経過に沿って測定されたデータであるが、必ずしも全時間のデータでなくてもよく、中性子束測定器のレンジ切替え後30秒程度、制御棒バンクの移動後最大100秒間等における好ましくないデータが排除されていてもよい。さらに、サンプリング間隔は、解析の精度と計算するデータ量の兼ね合い等から0.001秒間隔でサンプリングするのが望ましいが、必ずしもこれに限定されず、アナログデータを排除するものでもない。
また、炉出力時系列データ取得ステップにおける「所定の動特性モデル」とは、通常用いられている一点炉動特性モデルや炉心解析コード等を指す。
また、「所定平均の燃料温度」は、一次摂動理論やその他の解析で得られた値、経験値等を使用する。
また、「所定平均」とは、例えば「インポータンス加重平均」等である。
また、「等温温度反応度係数」とは、零出力試験時において、燃料温度のみの(で偏微分した)反応度係数と減速材温度のみの(で偏微分した)反応度係数の和を言う。
前記中性子束測定ステップにおける中性子束の時系列データの測定は、中性子束と共にγ線を測定するものであり、
前記反応度時系列データ取得ステップは、測定された中性子束の時系列データからγ線の影響を除洗する除洗手順を有し、さらにγ線の影響が除洗された中性子束の時系列データから一点炉動特性方程式に対する逆動特性法を用いて反応度の時系列データを求めるものであることを特徴とするドップラー反応度係数の測定方法である。
前記除洗手順は、反応度フィードバックの寄与が小さい低出力領域での炉出力応答について、一定炉周期反応度とγ線の混入割合をパラメータとして所定の原子炉動特性方程式を用いて求めた時間的変化解析値と、実測された中性子束の時系列データ中のこの炉出力応答に相当する時間的変化部分とを用いて定義されるものであって、両者の差を対数値で表す誤差関数を、最小自乗法を用いて評価し、誤差関数の値が最も小さくなる一定炉周期反応度とγ線の混入割合の組合せを求め、その組合せを構成するγ線の混入割合を真のγ線の混入割合とするものであることを特徴とするドップラー反応度係数の測定方法である。
なお、「反応度フィードバックの寄与が小さい低出力領域の時系列データ」とは、定格出力の1%以内の出力における時系列データであり、その様な領域のデータを採用するのは、反応度フィードバックの影響を補正せずにγ線の混入割合を正確に求めることが出来るからである。
また、炉内減速材の入口(低い)と出口(高い)の温度等、他の温度を求めてもよい。
前記炉出力時系列データ取得ステップにおける所定の手順は、蒸気発生器の1次側から2次側への伝熱に関する時定数と初期の炉出力をパラメータとし、{1−(減速材平均温度の最大温度への到達時間の解析値/減速材平均温度の最大温度への到達時間の測定値)}2+{1−(減速材平均温度の最大到達温度の解析値/減速材平均温度の最大到達温度の測定値)}2で表される誤差関数の値が最小となる前記時定数と初期の炉出力の組合せを求めることを特徴とするドップラー反応度係数の測定方法である。
また、到達上限出力値の絶対値と初期出力から到達上限出力までの出力変更幅が決定できるため、γ線を除洗したNIS信号から絶対値で表される炉出力応答を求めることができ、さらに求められた炉出力応答と実測された減速材温度から正確な燃料棒平均温度の時系列データを求めることができる。
前記燃料温度時系列データ取得ステップは、炉出力の時系列データと燃料棒平均温度に関する熱伝導方程式を用いて算出された体積加重平均の燃料温度に、零出力状態の減速材流路方向の中性子束および随伴中性子束(中性子インポータンス)の分布を考慮して作成した補正係数を用いて修正を行ない、一次摂動理論に適合した所定平均の燃料温度の時系列データを得るものであることを特徴とするドップラー反応度係数の測定方法である。
前記所定平均とはインポータンス出力加重平均であり、
前記ドップラー反応度係数取得ステップにおける所定の手順とは、「反応度寄与成分=ドップラー反応度係数×(インポータンス出力加重平均の燃料温度の時系列データを用いて得られた燃料温度の変化量−炉内減速材の平均温度の変化量)+等温温度反応度係数×減速材温度変化」という数式を用いるものであることを特徴とするドップラー反応度係数の測定方法である。
なお、「インポータンス出力加重平均」とは、中性子束分布と出力分布とがほぼ比例するとして、中性子インポータンス分布と中性子束分布とによる加重平均を指す。
前記ドップラー反応度係数取得ステップにおける所定の手順はさらに、ドップラー反応度係数をパラメータとし、採取したデータに対して、{1.0−ドップラー反応度係数×(インポータンス出力加重平均の燃料温度の時系列データを用いて得られた燃料温度の変化量−炉内減速材の平均温度の変化量)/ドップラー反応度係数に関する反応度寄与成分}2で定義される誤差関数の値が最小となるドップラー反応度係数を実際のドップラー反応度係数と推定するものであることを特徴とするドップラー反応度係数の測定方法である。
炉心特性の解析は、一点炉動特性シミュレーションモデル、一次元(冷却水の流路方向)動特性シミュレーションモデルで行い、この際炉出力が小さいため炉心の半径方向出力分布は平坦と仮定した。
(測定系)
図2に、本実施の形態における機器構成を含めた測定系を概念的に示す。図2において、11と12は炉外中性子束検出用の電離箱であり、21と22は熱電対等の温度センサーあり、31は電流計であり、50は微小電流計であり、51は直流増幅器であり、52は端子台であり、53はA/D変換ボード(ノートPC)である。
また、破線は計測用の信号線を示す。
データのサンプリングの時間間隔は、0.001秒であり、測定時間は2600秒であった。
以上の他、ローパスフィルター、アンプ等も使用した。
制御棒バンクの移動後最大100秒の間は、制御棒の移動による出力分布の空間変化が生じ、NIS信号にも影響が生じるため、この影響がある時間領域の応答は、フィッティングの(解析の)対象としていない。
また、未臨界から定格の1%程度まで炉出力を上昇させるため、NISの測定レンジ(測定対象の桁数)を切替える必要があるが、切替え後約30秒はその影響が出るため、フィッティングの対象としていない。
さらに、蒸気発生器2次系の動作環境が比較的安定していると考えられる上限値確認のピーク前までのデータを使用した。
制御棒を操作して、初期未臨界のPWRを臨界にし、さらに僅かに出力を上げ、その際の中性子束と減速材温度を測定した。
図3に、実測で得られた中性子束の時系列データを示す。図3において、縦軸が電流(A)に変換された中性子束であり、横軸はデータ取り込み開始からの経過時間である。なお、以降の物理量の時系列データを示す図においても、縦軸に物理量を、横軸にデータ取り込み開始からの経過時間を示す。
図4に、同じく減速材平均温度の時系列データを示す。図3と比較すれば、最大値のピークに約50秒の遅れが発生しているのが判る。理論解析から、この時間遅れは蒸気発生器の徐熱特性に拠るものであり、蒸気発生器の1次側から2次側への熱伝達(伝熱)に関する定数τsg、12が支配的影響を持ち、τsg、12が大きければ遅れも増大することが判った。また、τsg、12は34秒程度と判明した。
得られた中性子束データから、(1)式より、γ線の混入割合gc=0.78、一定炉周期反応度ρp=46.6pcmと言う数値が求められた。求められたγ線の混入割合gcを用いて、(2)式から得られた出力の時系列データを、図5に示す。図5において、縦軸は定格出力に対する炉出力Pの比を示し、実線はγ線を除洗後の中性子束を基に得られた炉出力Pの時系列データであり、破線は除洗前の中性子束を基に得られた炉出力Pの時系列データである。なお、破線は、基本的には電流値で示した図3と同じである。
図5において、炉出力Pが小さい領域では、測定した中性子束には、γ線の影響が大きく出ていることが判る。
また、炉出力の変化幅は、約220倍であった。
図1に示すモデルを対象にして、前記(一次冷却ループ除熱モデル)で説明した手法を用いて、τsg,12=34s、P0=3.77×10−4%/定格出力と求められ、到達上限炉出力Pmax=8.35×10−2%/定格出力となった。また、それらの値を用いて減速材の平均温度Tc,avの時系列データを求めた。結果を、図7に示す。図7において、実線が計算値であり、破線が測定値である。フィッティング対象となるピーク到達時間以前の応答において、ピーク到達時刻およびピーク値を含めて、測定値と計算値に差異は認められない。
測定した中性子束分布から、(12)式を用いてインポータンス加重平均の燃料の(燃料棒の)平均温度変化の補正係数値cipを求めた結果、1.296となった。なお、中性子束は、ドップラー効果であるため高速群を用いたが、熱中性子を用いても結果に差は生じなかった。さらに、この値を用いて、インポータンス加重平均の燃料の温度を求めた。
(13)式により、誤差関数が最小となるドップラー反応度係数αfを推定した結果、αf=−3.2[pcm/K]となり、設計値と有効数字2桁で一致する。
11 電離箱
12 電離箱
20 蒸気発生器
21 温度センサー
22 温度センサー
30 冷却水循環ポンプ
31 電流計
41 原子炉出口側配管
42 原子炉入口側配管
50 微小電流計
51 直流増幅器
52 端子台
53 A/D変換ボード(ノートPC)
Claims (8)
- 未臨界または臨界になった状態の
炉心に反応度を印加して炉出力を一定炉周期で所定量上昇させ、その際の中性子束を時系列データとして測定する中性子束時系列データ測定ステップと、
未臨界または臨界になった状態の炉心に反応度を印加して炉出力を一定炉周期で所定量上昇させ、その際の炉内の減速材の平均温度を所定の手順で時系列データとして取得する炉内減速材平均温度時系列データ取得ステップと、
測定された中性子束の時系列データから、一点炉動特性方程式に対する逆動特性法を用いて反応度の時系列データを取得する反応度時系列データ取得ステップと、
前記取得された炉内減速材平均温度の時系列データと中性子束の時系列データを基に、所定の手順で前記2つの時系列データに整合する炉出力の時系列データを取得する炉出力時系列データ取得ステップと、
取得した炉出力の時系列データと所定の動特性モデルを用いて取得された所定平均の燃料温度の時系列データを得る燃料温度時系列データ取得ステップと、
取得した反応度の時系列データと印加した反応度を用いて反応度フィードバック寄与成分の時系列データを求める反応度フィードバック寄与成分の時系列データ取得ステップと、
前記炉内の減速材平均温度の時系列データと、前記所定平均の燃料温度の時系列データと、等温温度反応度係数と、前記反応度フィードバック寄与成分の時系列データを用いて、所定の手順でドップラー反応度係数を求めるドップラー反応度係数取得ステップを有していることを特徴とするドップラー反応度係数の測定方法。 - 前記中性子束測定ステップにおける中性子束の時系列データの測定は、中性子束と共にγ線を測定するものであり、
前記反応度時系列データ取得ステップは、測定された中性子束の時系列データからγ線の影響を除洗する除洗手順を有し、さらにγ線の影響が除洗された中性子束の時系列データから一点炉動特性方程式に対する逆動特性法を用いて反応度の時系列データを求めるものであることを特徴とする請求項1に記載のドップラー反応度係数の測定方法。 - 前記除洗手順は、反応度フィードバックの寄与が小さい低出力領域での炉出力応答について、一定炉周期反応度とγ線の混入割合をパラメータとして所定の原子炉動特性方程式を用いて求めた時間的変化解析値と、実測された中性子束の時系列データ中のこの炉出力応答に相当する時間的変化部分とを用いて定義されるものであって、両者の差を対数値で表す誤差関数を、最小自乗法を用いて評価し、誤差関数の値が最も小さくなる一定炉周期反応度とγ線の混入割合の組合せを求め、その組合せを構成するγ線の混入割合を真のγ線の混入割合とするものであることを特徴とする請求項2に記載のドップラー反応度係数の測定方法。
- 前記炉内減速材平均温度時系列データ取得ステップは、未臨界または臨界になった状態の炉心に反応度を印加して炉出力を一定炉周期で所定量上昇させる際の減速材平均温度を時系列データの形で取得するものであることを特徴とする請求項1ないし請求項3のいずれかに記載のドップラー反応度係数の測定方法。
- 前記炉出力時系列データ取得ステップにおける所定の手順は、蒸気発生器の1次側から2次側への伝熱に関する時定数と初期の炉出力をパラメータとし、{1−(減速材平均温度の最大温度への到達時間の解析値/減速材平均温度の最大温度への到達時間の測定値)}2+{1−(減速材平均温度の最大到達温度の解析値/減速材平均温度の最大到達温度の測定値)}2で表される誤差関数の値が最小となる前記時定数と初期の炉出力の組合せを求めることを特徴とする請求項4に記載のドップラー反応度係数の測定方法。
- 前記燃料温度時系列データ取得ステップは、炉出力の時系列データと燃料棒平均温度に関する熱伝導方程式を用いて算出された体積加重平均の燃料温度に、零出力状態の減速材流路方向の中性子束および随伴中性子束(中性子インポータンス)の分布を考慮して作成した補正係数を用いて修正を行ない、一次摂動理論に適合した所定平均の燃料温度の時系列データを得るものであることを特徴とする請求項1ないし請求項5のいずれかに記載のドップラー反応度係数の測定方法。
- 前記所定平均とはインポータンス出力加重平均であり、
前記ドップラー反応度係数取得ステップにおける所定の手順とは、「反応度寄与成分=ドップラー反応度係数×(インポータンス出力加重平均の燃料温度の時系列データを用いて得られた燃料温度の変化量−炉内減速材の平均温度の変化量)+等温温度反応度係数×減速材温度変化」という数式を用いるものであることを特徴とする請求項1ないし請求項6のいずれかに記載のドップラー反応度係数の測定方法。 - 前記ドップラー反応度係数取得ステップにおける所定の手順はさらに、ドップラー反応度係数をパラメータとし、採取したデータに対して、{1.0−ドップラー反応度係数×(インポータンス出力加重平均の燃料温度の時系列データを用いて得られた燃料温度の変化量−炉内減速材の平均温度の変化量)/ドップラー反応度係数に関する反応度寄与成分}2で定義される誤差関数の値が最小となるドップラー反応度係数を実際のドップラー反応度係数と推定するものであることを特徴とする請求項7に記載のドップラー反応度係数の測定方法。
Priority Applications (7)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007197130A JP4901631B2 (ja) | 2007-07-30 | 2007-07-30 | ドップラー反応度係数の測定方法 |
RU2010107215/07A RU2491664C2 (ru) | 2007-07-30 | 2008-07-29 | Способ измерения доплеровского коэффициента реактивности |
CN2008801008323A CN101809674B (zh) | 2007-07-30 | 2008-07-29 | 多普勒反应性系数测量方法 |
EP08791792.8A EP2172943B1 (en) | 2007-07-30 | 2008-07-29 | Doppler reactivity coefficient measuring method |
PCT/JP2008/063557 WO2009017103A1 (ja) | 2007-07-30 | 2008-07-29 | ドップラー反応度係数の測定方法 |
KR1020107003539A KR101488549B1 (ko) | 2007-07-30 | 2008-07-29 | 도플러 반응도계수의 측정방법 |
US12/694,729 US8462909B2 (en) | 2007-07-30 | 2010-01-27 | Doppler reactivity coefficient measuring method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007197130A JP4901631B2 (ja) | 2007-07-30 | 2007-07-30 | ドップラー反応度係数の測定方法 |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2009031186A JP2009031186A (ja) | 2009-02-12 |
JP2009031186A5 JP2009031186A5 (ja) | 2010-06-24 |
JP4901631B2 true JP4901631B2 (ja) | 2012-03-21 |
Family
ID=40304333
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2007197130A Active JP4901631B2 (ja) | 2007-07-30 | 2007-07-30 | ドップラー反応度係数の測定方法 |
Country Status (7)
Country | Link |
---|---|
US (1) | US8462909B2 (ja) |
EP (1) | EP2172943B1 (ja) |
JP (1) | JP4901631B2 (ja) |
KR (1) | KR101488549B1 (ja) |
CN (1) | CN101809674B (ja) |
RU (1) | RU2491664C2 (ja) |
WO (1) | WO2009017103A1 (ja) |
Families Citing this family (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7894565B2 (en) * | 2008-12-11 | 2011-02-22 | Westinghouse Electric Company Llc | Subcritical reactivity measurement method |
JP5491879B2 (ja) * | 2009-02-16 | 2014-05-14 | 株式会社東芝 | 中性子増倍体系の未臨界度判定装置、及び未臨界度判定プログラム |
JP5172777B2 (ja) * | 2009-05-15 | 2013-03-27 | 株式会社東芝 | 温度反応度変化測定装置および方法 |
US8596837B1 (en) | 2009-07-21 | 2013-12-03 | Cooper Technologies Company | Systems, methods, and devices providing a quick-release mechanism for a modular LED light engine |
WO2011011323A1 (en) | 2009-07-21 | 2011-01-27 | Cooper Technologies Company | Interfacing a light emitting diode (led) module to a heat sink assembly, a light reflector and electrical circuits |
JP5725792B2 (ja) * | 2010-10-25 | 2015-05-27 | 三菱重工業株式会社 | 減速材の温度係数測定装置および減速材の温度係数測定方法 |
CN102279901B (zh) * | 2011-05-17 | 2013-01-09 | 湖北省电力公司电力科学研究院 | 一种针对第三代压水堆核电机组的建模方法 |
JP6091765B2 (ja) * | 2012-03-26 | 2017-03-08 | 三菱重工業株式会社 | シミュレーション装置、シミュレーション方法および連係プログラム |
JP6091764B2 (ja) * | 2012-03-26 | 2017-03-08 | 三菱重工業株式会社 | シミュレーション装置、シミュレーション方法および連係プログラム |
CN102800372B (zh) * | 2012-08-17 | 2015-03-11 | 中国原子能科学研究院 | 一种反应堆反应性测量值的外源修正方法 |
JP6139175B2 (ja) * | 2013-02-25 | 2017-05-31 | 三菱重工業株式会社 | 反応度温度係数推定装置及び方法 |
CN103268728B (zh) * | 2013-04-27 | 2016-01-20 | 国家电网公司 | 一种电力系统动态仿真压水反应堆控制系统模型构建方法 |
CA2955978C (en) | 2014-08-28 | 2021-10-19 | Terrapower, Llc | Doppler reactivity augmentation device |
CN106156460B (zh) * | 2015-03-27 | 2020-03-17 | 国核(北京)科学技术研究院有限公司 | 获得核燃料棒内部燃料芯块温度分布的方法 |
CN105448362B (zh) * | 2015-11-23 | 2017-06-30 | 中核核电运行管理有限公司 | 利用功率量程探测器和信号作为信号源的动态刻棒方法 |
CN106297920B (zh) * | 2016-08-23 | 2018-01-05 | 中国科学院合肥物质科学研究院 | 一种无空间效应的次临界度逆动态测量方法 |
CN111149175B (zh) * | 2017-08-18 | 2023-04-07 | 西屋电气有限责任公司 | 核仪表隔离输出信号标度方法和采用相同方法的系统 |
CN107763610B (zh) * | 2017-09-28 | 2019-05-24 | 哈尔滨工程大学 | 一种基于结构的核电厂高压给水加热器精细化仿真方法 |
CN107887043B (zh) * | 2017-10-10 | 2019-09-17 | 福建福清核电有限公司 | 一种反应堆功率系数测量方法 |
EP3747027B1 (en) | 2018-02-02 | 2022-11-23 | Westinghouse Electric Company Llc | Nuclear fuel failure protection method |
CN110739091B (zh) * | 2019-10-22 | 2022-05-10 | 上海核工程研究设计院有限公司 | 一种核电厂次临界度直接测量系统及硬件平台和测量方法 |
CN113466922B (zh) * | 2020-03-31 | 2024-03-08 | 刘畅源 | 一种核截面多普勒展宽方法和装置 |
CN111723506B (zh) * | 2020-06-22 | 2022-03-25 | 中国核动力研究设计院 | 一种系统级分析模型各部件动力贡献度分析方法及系统 |
CN113806941B (zh) * | 2021-09-22 | 2024-01-05 | 上海核星核电科技有限公司 | 一种具有氙瞬态模拟能力的压水堆燃耗跟踪计算方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3916444A (en) * | 1973-02-23 | 1975-10-28 | Westinghouse Electric Corp | Training simulator for nuclear power plant reactor monitoring |
US4075059A (en) * | 1976-04-28 | 1978-02-21 | Combustion Engineering, Inc. | Reactor power reduction system and method |
FR2546329B1 (fr) * | 1983-05-19 | 1985-08-23 | Framatome Sa | Procede de detection des variations de la reactivite du coeur d'un reacteur nucleaire a eau sous pression et dispositif de mise en oeuvre de ce procede |
JP3137569B2 (ja) * | 1995-10-11 | 2001-02-26 | 三菱重工業株式会社 | 原子炉の中性子源強度及びガンマ線強度を評価する方法 |
JP2006084181A (ja) * | 2004-09-14 | 2006-03-30 | Hokkaido Univ | 加圧水型軽水炉の温度反応度係数分離測定方法 |
JP2006105814A (ja) * | 2004-10-06 | 2006-04-20 | Mitsubishi Heavy Ind Ltd | 減速材温度係数の測定方法及び減速材温度係数測定装置 |
RU2280908C1 (ru) * | 2004-12-31 | 2006-07-27 | Открытое акционерное общество "Всероссийский научно-исследовательский институт по эксплуатации атомных электростанций" (ОАО "ВНИИАЭС") | Способ определения парового коэффициента реактивности ядерного энергетического реактора типа рбмк |
-
2007
- 2007-07-30 JP JP2007197130A patent/JP4901631B2/ja active Active
-
2008
- 2008-07-29 RU RU2010107215/07A patent/RU2491664C2/ru active
- 2008-07-29 WO PCT/JP2008/063557 patent/WO2009017103A1/ja active Application Filing
- 2008-07-29 CN CN2008801008323A patent/CN101809674B/zh not_active Expired - Fee Related
- 2008-07-29 KR KR1020107003539A patent/KR101488549B1/ko active IP Right Grant
- 2008-07-29 EP EP08791792.8A patent/EP2172943B1/en not_active Not-in-force
-
2010
- 2010-01-27 US US12/694,729 patent/US8462909B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
RU2010107215A (ru) | 2011-09-10 |
EP2172943B1 (en) | 2013-04-24 |
US20100128832A1 (en) | 2010-05-27 |
RU2491664C2 (ru) | 2013-08-27 |
CN101809674B (zh) | 2013-09-25 |
KR101488549B1 (ko) | 2015-02-02 |
KR20100042281A (ko) | 2010-04-23 |
JP2009031186A (ja) | 2009-02-12 |
EP2172943A1 (en) | 2010-04-07 |
CN101809674A (zh) | 2010-08-18 |
EP2172943A4 (en) | 2012-01-18 |
US8462909B2 (en) | 2013-06-11 |
WO2009017103A1 (ja) | 2009-02-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4901631B2 (ja) | ドップラー反応度係数の測定方法 | |
JP5542150B2 (ja) | 臨界未満反応度測定方法 | |
JP6771553B2 (ja) | 迅速応答自己給電型炉内検出器を用いる臨界未満反応度監視装置 | |
JP6400685B2 (ja) | 原子炉運転停止時におけるホウ素希釈監視方法 | |
EP2296153A2 (en) | Neutron monitoring systems including gamma thermometers and methods of calibrating nuclear instruments using gamma thermometers | |
JP2015148524A (ja) | 原子炉の制御棒価値を求める方法、プログラム、記録媒体、およびシステム | |
CN111587460A (zh) | 亚临界堆芯反应性偏差预测技术 | |
Saito et al. | Measurement and evaluation on pulsing characteristics and experimental capability of NSRR | |
Kim et al. | A high-fidelity Monte Carlo evaluation of CANDU-6 safety parameters | |
Demazière et al. | Development of a method for measuring the moderator temperature coefficient by noise analysis and its experimental verification in Ringhals-2 | |
JP2012163379A (ja) | 燃料集合体ガンマ線測定装置 | |
JP5491879B2 (ja) | 中性子増倍体系の未臨界度判定装置、及び未臨界度判定プログラム | |
Buell et al. | A neutron scatterometer for void-fraction measurement in heated rod-bundle channels under CANDU LOCA conditions | |
Tonoike et al. | Real time a value measurement with Feynman-α method utilizing time series data acquisition on low enriched uranium system | |
JP4299995B2 (ja) | 沸騰水型原子炉の炉心流量計測演算システム | |
JP3785847B2 (ja) | 原子炉出力測定装置 | |
Lepore et al. | On Capabilities and Limitations of Current Fast Neutron Flux Monitoring Instrumentation for the Demo LFR ALFRED | |
Mutneja et al. | Design of discrete-time Kalman filter based delay compensator for Rhodium-103 SPND | |
Fehrenbach et al. | Description of the blowdown test facility COG program on in-reactor fission product release, transport, and deposition under severe accident conditions | |
Demaziere et al. | Analysis of an MTC noise measurement performed in ringhals-2 using gamma-thermometers and in-core neutron detectors | |
Tsuji et al. | Estimating temperature reactivity coefficients by experimental procedures combined with isothermal temperature coefficient measurements and dynamic identification | |
Z Mesquita et al. | An innovative method for online power monitoring in nuclear reactors | |
Tsuji et al. | Estimation of the Doppler Coefficient from a Lower Power Transient Observed in a Zero-Power Reactor Physics Test of a PWR (I)—Methodology— | |
Ammon et al. | Monitoring the Integrity of Control Rods On-Line with a Helium Leak Detector | |
TSUJI et al. | Identification of the Doppler Coefficient from a Low Power Transient Observed in a Zero-Power Reactor Physics Test of PWRs |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20100507 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20100507 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20111220 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20111227 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4901631 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20150113 Year of fee payment: 3 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |