JPS58205894A - 原子炉における中性子検出器測定データの処理方法 - Google Patents
原子炉における中性子検出器測定データの処理方法Info
- Publication number
- JPS58205894A JPS58205894A JP57088785A JP8878582A JPS58205894A JP S58205894 A JPS58205894 A JP S58205894A JP 57088785 A JP57088785 A JP 57088785A JP 8878582 A JP8878582 A JP 8878582A JP S58205894 A JPS58205894 A JP S58205894A
- Authority
- JP
- Japan
- Prior art keywords
- tip
- neutron
- neutron absorber
- measurement value
- 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.)
- Granted
Links
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)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
本発明は、特に沸騰水型動力炉において、走行型中性子
検出器（以下本明細書においてはＴＩＰと略称する。）
の測定データを処理する方法に関する。
検出器（以下本明細書においてはＴＩＰと略称する。）
の測定データを処理する方法に関する。
ＴＩＰは、第１図に示す機構で炉心内の導管を移動し、
炉心高さ方向の中性子束分布を測定する。
炉心高さ方向の中性子束分布を測定する。
図中、１は原子炉容器、２は炉心、３は炉心の頂上、４
は炉心の底、５はＴＩＰ、６は導管、７はワイヤ、８は
駆動装置、９はＸ−Ｙレコーダ、１０はプロセス計算機
で、ＡはＴＩＰ測定値、ＢはＴＩＰ５の炉心高さ方向位
置を表わす信号の通路を意味する。測定結果は、第２図
に示すように、Ｘ−Ｙレコーダ９上に、縦軸がＴＩＰに
よる中性子束分布の測定値（以下本明細書においては１
’−ＴＩＰ測定値」と略称する。）、横軸が炉心高さ方
向の位置（以下本明細書においては「軸方向位置」と略
称する。）として出力されるとともに、プロセス計算機
１０に取込まれて処理され、炉心内の出力分布を計算す
るための基準データとして使用される。
は炉心の底、５はＴＩＰ、６は導管、７はワイヤ、８は
駆動装置、９はＸ−Ｙレコーダ、１０はプロセス計算機
で、ＡはＴＩＰ測定値、ＢはＴＩＰ５の炉心高さ方向位
置を表わす信号の通路を意味する。測定結果は、第２図
に示すように、Ｘ−Ｙレコーダ９上に、縦軸がＴＩＰに
よる中性子束分布の測定値（以下本明細書においては１
’−ＴＩＰ測定値」と略称する。）、横軸が炉心高さ方
向の位置（以下本明細書においては「軸方向位置」と略
称する。）として出力されるとともに、プロセス計算機
１０に取込まれて処理され、炉心内の出力分布を計算す
るための基準データとして使用される。
従来性なわれてきた典型的な処理の手順はつぎの通りで
ある。
ある。
（１）ＴＩＰ５を導管６内に挿入して行き、炉心の頂上
３に到達したことを、導管６内のリミット・スイッチで
検出する。
３に到達したことを、導管６内のリミット・スイッチで
検出する。
（１＋）ＴＩＰ５を引抜きながら、炉心の頂上３から底
４までの間、軸方向位置１インチ毎にＴＩＰ測定値を収
集する。炉心の萬さが１４４インチの場合、ＴＩＰ測定
値の個数は１４５となる。
４までの間、軸方向位置１インチ毎にＴＩＰ測定値を収
集する。炉心の萬さが１４４インチの場合、ＴＩＰ測定
値の個数は１４５となる。
０ｉｉ）ｉ４５個のＴＩＰ測定値から基準データを作成
する。例えば、炉心高さ方向を４８個のノードに分割し
、各ノードの基準データを作成する場合、次式のような
方法がある。
する。例えば、炉心高さ方向を４８個のノードに分割し
、各ノードの基準データを作成する場合、次式のような
方法がある。
ただし、
Ｎ＝３・Ｋ−２・・・・・・・・・（２）ここで、Ｘ（
Ｋ） ：炉心の底からに番目の基準データ ｘ（Ｎ）：炉心の底から（Ｎ−１）イ ンチの軸方向位置でのＴＩＰ 測定値 に＝１．４８ 上記処理の主眼は、式（１）の処理によｊｏ、ＴＩＰ測
定値のノイズを除去するこ≧である。第２図から判るよ
うに、ＴＩＰ測定値は細かく振動している。この振動成
分（すなわちノイズ）を、＠接するＴＩＰ測定値間での
平均操作により、除去し、出力分布計算用の基準データ
を作成する。
Ｋ） ：炉心の底からに番目の基準データ ｘ（Ｎ）：炉心の底から（Ｎ−１）イ ンチの軸方向位置でのＴＩＰ 測定値 に＝１．４８ 上記処理の主眼は、式（１）の処理によｊｏ、ＴＩＰ測
定値のノイズを除去するこ≧である。第２図から判るよ
うに、ＴＩＰ測定値は細かく振動している。この振動成
分（すなわちノイズ）を、＠接するＴＩＰ測定値間での
平均操作により、除去し、出力分布計算用の基準データ
を作成する。
しかし、ＴＩＰ測定値には、ノイズの他に、軸方向位置
のずれの問題がある。第１図において、走行中のＴＩＰ
５の軸方向位置は、駆動装置８においてワイヤ７の送出
し・引込み量を計数することによって判定される。した
がって、ワイヤ７にたるみが生じた場合、また導管６に
たわみが生じた場合には、実際にＴＩＰ測定値を得た軸
方向位置と、駆動装置８から出力される軸方向位置との
間にずれが生ずる。例えば、第２図の軸方向位置１８〜
３６インチの部分では、ずれが１インチあると、ＴＩＰ
測定値は約３Ｗ／ｚ”変化する。これＦｉ’ｌ’ＩＰ測
定値の平均の８％に相尚する。
のずれの問題がある。第１図において、走行中のＴＩＰ
５の軸方向位置は、駆動装置８においてワイヤ７の送出
し・引込み量を計数することによって判定される。した
がって、ワイヤ７にたるみが生じた場合、また導管６に
たわみが生じた場合には、実際にＴＩＰ測定値を得た軸
方向位置と、駆動装置８から出力される軸方向位置との
間にずれが生ずる。例えば、第２図の軸方向位置１８〜
３６インチの部分では、ずれが１インチあると、ＴＩＰ
測定値は約３Ｗ／ｚ”変化する。これＦｉ’ｌ’ＩＰ測
定値の平均の８％に相尚する。
このＴＩＰのずれに基因する、系統的に現れる中性子測
定値の誤差を補正する方法はこれまで知られていなかっ
た。
定値の誤差を補正する方法はこれまで知られていなかっ
た。
１１１１１１
本発明の目的は、ＴＩＰ測定値の軸方向位置のずれを自
動的に検出し、これを補正するＴＩＰ測定値処理方法を
提供することによｊ９、ＴＩＰ測定値から作成した基準
データを使用した出力分布計算の精度を向上することに
ある。
動的に検出し、これを補正するＴＩＰ測定値処理方法を
提供することによｊ９、ＴＩＰ測定値から作成した基準
データを使用した出力分布計算の精度を向上することに
ある。
上記目的を達成す名ために、本発明による原子炉におけ
る中性子検出測定データの処理方法は、走行型中性子検
出器の測定値から中性子吸収体の位置を検出する処理と
、実際に中性子吸収体の存在する位置と測定値から検出
された位置との差から測定値に対応する上記走行型中性
子検出器の駆動装置から算出される炉心高さ方向の位置
の誤差を検出し、測定値を補正する処理とを含むことを
要旨とする。より具体的には、本発明による原子炉にお
ける中性子検出測定データの処理方法においては、走行
型中性子検出器の正しい位置を検出するために、燃料集
合体スペーサおよび局所中性子検出器（以下本明Ｍ８Ｎ
においてａＬＰＲＭと略称する。）の軸方向位置が利用
される。
る中性子検出測定データの処理方法は、走行型中性子検
出器の測定値から中性子吸収体の位置を検出する処理と
、実際に中性子吸収体の存在する位置と測定値から検出
された位置との差から測定値に対応する上記走行型中性
子検出器の駆動装置から算出される炉心高さ方向の位置
の誤差を検出し、測定値を補正する処理とを含むことを
要旨とする。より具体的には、本発明による原子炉にお
ける中性子検出測定データの処理方法においては、走行
型中性子検出器の正しい位置を検出するために、燃料集
合体スペーサおよび局所中性子検出器（以下本明Ｍ８Ｎ
においてａＬＰＲＭと略称する。）の軸方向位置が利用
される。
ＴＩＰの導管に隣接し、炉心に固定されている、燃料ビ
ン支持用のスペーサおよびＬＰＲＭは、中性子吸収材の
役割をする。したがって、スペーサおよびＬＰＲＭのあ
る軸方向位置では、ＴＩＰ測定値は減少する。たとえば
、第２図において、５ｔ−ｓｔｏはスペーサおよびＬＰ
ＲＭの配置例である（破線がスペーサ、一点鎖線がＬＰ
ＲＭの位置を示し、ｓｌでは、スペーサとＬＰＲＭとが
近接しているので、まとめて−個所と見なした）。
ン支持用のスペーサおよびＬＰＲＭは、中性子吸収材の
役割をする。したがって、スペーサおよびＬＰＲＭのあ
る軸方向位置では、ＴＩＰ測定値は減少する。たとえば
、第２図において、５ｔ−ｓｔｏはスペーサおよびＬＰ
ＲＭの配置例である（破線がスペーサ、一点鎖線がＬＰ
ＲＭの位置を示し、ｓｌでは、スペーサとＬＰＲＭとが
近接しているので、まとめて−個所と見なした）。
ＴＩＰ測定値のグラフは、８１〜８１Ｇの位置付近で、
局所的に減少してくぼみを持っている。もし吸収材がな
ければ、各吸収材の近傍で点線で平滑化されている箇所
における実線表示のようなくぼみの発生はない。このと
き、くぼみの深さはＴＩＰの測定値が大きければ大きい
程大きくなる。したがって、図ではＳ！近傍のくぼみが
最も大きい。
局所的に減少してくぼみを持っている。もし吸収材がな
ければ、各吸収材の近傍で点線で平滑化されている箇所
における実線表示のようなくぼみの発生はない。このと
き、くぼみの深さはＴＩＰの測定値が大きければ大きい
程大きくなる。したがって、図ではＳ！近傍のくぼみが
最も大きい。
特に、ＳｌにスペーサとＬＰＲＭの２種の吸収材がある
ことがくぼみを一層大きくしている。そこで、実際の８
１〜８１Ｇの位置と、ＴＩＰ測定値のくぼみの位置との
差を求めれば、ずれの検出・補正が可能である。
ことがくぼみを一層大きくしている。そこで、実際の８
１〜８１Ｇの位置と、ＴＩＰ測定値のくぼみの位置との
差を求めれば、ずれの検出・補正が可能である。
これが本発明の処理方法の基本的考え方であるが、ずれ
の検出・修正を精度よく、かつ効率よ〈実施するために
、以下の５ステツプが必要である（第３図）。
の検出・修正を精度よく、かつ効率よ〈実施するために
、以下の５ステツプが必要である（第３図）。
（１）平滑化処理：中性子吸収体によるくぼみが消えな
いようにしながら、ノイズを除去する。
いようにしながら、ノイズを除去する。
（２）準備処理二平滑化したＴＩＰ測定値の軸方向の差
分を計算するなど、Ｔ工Ｐ測定値のくぼみを使った中性
子吸収体の位置決めに必要な準備をする。
分を計算するなど、Ｔ工Ｐ測定値のくぼみを使った中性
子吸収体の位置決めに必要な準備をする。
（３）中性子吸収体の概略位置推定処理：中性子吸収体
の大まかな位置を選び出す。
の大まかな位置を選び出す。
（４） 中性子吸収体の位置決定処理：概略位置の近
傍から、中性子吸収体の位置を決定する。
傍から、中性子吸収体の位置を決定する。
（５）ＴＩＰ測定値の補正処理：（１）〜（４）のステ
ップで、ＴＩＰ測定値から得た中性子吸収体の位置と、
実際に中性子吸収体のある位置との差から、ＴＩＰ測定
値を補正する。
ップで、ＴＩＰ測定値から得た中性子吸収体の位置と、
実際に中性子吸収体のある位置との差から、ＴＩＰ測定
値を補正する。
以下、本発明の処理方法を実施例により詳細に説明する
。説明に際しては、第２図のＴＩＰ測定データが得られ
た場合を想定し、第３図の５ステツプの処理ごとに、実
施例を示す。
。説明に際しては、第２図のＴＩＰ測定データが得られ
た場合を想定し、第３図の５ステツプの処理ごとに、実
施例を示す。
（１）平滑化処理１００
ノイズを除去する代表的な方法は、フィルタを使用する
ことである。第２図から判るように、ＴＩＰ測定値のノ
イズは、中性子吸収体８１〜Ｓ１０によるくぼみに比べ
て、周波数が大きい。したがって、高調波を遮断するフ
ィルタを付加すれば、くぼみを残したまま、ノイズを除
去できる。
ことである。第２図から判るように、ＴＩＰ測定値のノ
イズは、中性子吸収体８１〜Ｓ１０によるくぼみに比べ
て、周波数が大きい。したがって、高調波を遮断するフ
ィルタを付加すれば、くぼみを残したまま、ノイズを除
去できる。
しかし、より簡便な方法として、従来から使用してきた
隣接するＴＩＰ測定値間での平均操作により、ノイズが
除去できる。この時、中性子吸収体によるくぼみを保存
するためには、平均をとるＴＩＰ測定値の個数を、中性
子吸収体による影響の及ぶ範囲を考慮して、決定すれば
よい。例えば、第２図の場合には、従来から使用してき
た式（１）を用いて、計４８個の平滑化したＴＩＰ測定
値Ｘ（ＩＯを求めれば、くぼみを保存したままノイズが
除去できる。その実施例を第４図に、結果を第５図に示
す。
隣接するＴＩＰ測定値間での平均操作により、ノイズが
除去できる。この時、中性子吸収体によるくぼみを保存
するためには、平均をとるＴＩＰ測定値の個数を、中性
子吸収体による影響の及ぶ範囲を考慮して、決定すれば
よい。例えば、第２図の場合には、従来から使用してき
た式（１）を用いて、計４８個の平滑化したＴＩＰ測定
値Ｘ（ＩＯを求めれば、くぼみを保存したままノイズが
除去できる。その実施例を第４図に、結果を第５図に示
す。
全長が１４４インチあるので、４８個の平滑化したＴＩ
Ｐ測定値Ｘ（Ｋ）は、３インチ毎の平均となる。３イン
チ毎の平均で中性子吸収体にくぼみが保存されるのは、
中性子吸収体の大きさを設計変更しない限シ、同様であ
る。将来、中性子吸収体の大きさが変更された場合には
、平均をとるＴＩＰ測定値の個数を増減することにより
、本実施例の方法で対応可能である。
Ｐ測定値Ｘ（Ｋ）は、３インチ毎の平均となる。３イン
チ毎の平均で中性子吸収体にくぼみが保存されるのは、
中性子吸収体の大きさを設計変更しない限シ、同様であ
る。将来、中性子吸収体の大きさが変更された場合には
、平均をとるＴＩＰ測定値の個数を増減することにより
、本実施例の方法で対応可能である。
従来の方法では、式（１）で求めたＸ （Ｋ）を、その
まま出力分布計算用の基準データとして使用していたが
、本発明の方法では、さらに以下の処理を加える。
まま出力分布計算用の基準データとして使用していたが
、本発明の方法では、さらに以下の処理を加える。
（２）準備処理２００
このステップでは、平滑化したＴＩＰ測定値Ｘ （Ｋ）
を処理して、後に続く２ステツプ（中性子吸収体の概略
位置推定３００および位置決定４００）で必要なデータ
を提供をする。したがって、処理の内容は後続の２ステ
ツプに依存する。
を処理して、後に続く２ステツプ（中性子吸収体の概略
位置推定３００および位置決定４００）で必要なデータ
を提供をする。したがって、処理の内容は後続の２ステ
ツプに依存する。
ここでは、後続の２ステツプで、ＴＩＰ測定値Ｘ（Ｋ）
の差分を必要とする場合を例にとシ、第６図に処理内容
を示す。まず、次式によ！１７、ＴＩＰ測定値Ｘ（Ｋ）
の差分を計算する。
の差分を必要とする場合を例にとシ、第６図に処理内容
を示す。まず、次式によ！１７、ＴＩＰ測定値Ｘ（Ｋ）
の差分を計算する。
ΔＸ（Ｌ）＝Ｘ（Ｋ＋１３−Ｘ（Ｋ） ・・・・
・・・・・・・・（３）（９）゛ ただし、 Ｋ＝Ｌ ・・・・・・・・・・・・（
４）ここで、ΔＸ（Ｌ）：ノードにと（Ｋ＋１）との間
の差分 Ｌ＝１．４７ 第７図は、ΔＸ（Ｌ）を図示したものである。各点はノ
ードにと（Ｋ＋１）との境界線上にプロットしである。
・・・・・・・・（３）（９）゛ ただし、 Ｋ＝Ｌ ・・・・・・・・・・・・（
４）ここで、ΔＸ（Ｌ）：ノードにと（Ｋ＋１）との間
の差分 Ｌ＝１．４７ 第７図は、ΔＸ（Ｌ）を図示したものである。各点はノ
ードにと（Ｋ＋１）との境界線上にプロットしである。
つぎに、谷差分ΔＸ（Ｌ）の差分を次式により計算する
。
。
Δ”Ｘ（Ｍ）＝ΔＸ（Ｌ＋１）−ΔＸ（Ｌ） ・・・
・・・（５）ただし、 Ｌ＝Ｍ ・・・・・・・・・・・・・・
・（６）ここで、Δ意Ｘ（Ｍ）：ΔＸ（Ｌ）とΔＸ（Ｌ
＋１） との間の差分 Ｍ＝１．４６ 第８図はΔ”Ｘ（Ｍ） を図示したものである。
・・・（５）ただし、 Ｌ＝Ｍ ・・・・・・・・・・・・・・
・（６）ここで、Δ意Ｘ（Ｍ）：ΔＸ（Ｌ）とΔＸ（Ｌ
＋１） との間の差分 Ｍ＝１．４６ 第８図はΔ”Ｘ（Ｍ） を図示したものである。
（３）中性子吸収体の概略位置推定処理３００このステ
ップでは、準備処理２００の結果を使用して、中性子吸
収体のある位置をノード単位で推定する。すなわち、第
５図に示すＴＩＰ測定値（１０） Ｘ （Ｋ）のグラフにおいて、くぼみを探す処理が必要
である。
ップでは、準備処理２００の結果を使用して、中性子吸
収体のある位置をノード単位で推定する。すなわち、第
５図に示すＴＩＰ測定値（１０） Ｘ （Ｋ）のグラフにおいて、くぼみを探す処理が必要
である。
第９図に示す実施例は、準備処理２００で求めたΔ”Ｘ
（Ｍ）を使って、くぼみを探す処理の手順である。第５
図から判るように、くぼみがあると、その付近でＴＩＰ
測定値Ｘ（ＩＯの勾配（すなわちΔＸ（Ｌ））は一時減
少し、再び上昇する。したがって、Δ”Ｘ（Ｍ） （す
なわちΔＸ（Ｌ）の勾配）は負から正に大きく変化する
。第９図は、このような場所を炉心の底から頂上に向っ
て順次探し出し、次式の処理をするための手順である。
（Ｍ）を使って、くぼみを探す処理の手順である。第５
図から判るように、くぼみがあると、その付近でＴＩＰ
測定値Ｘ（ＩＯの勾配（すなわちΔＸ（Ｌ））は一時減
少し、再び上昇する。したがって、Δ”Ｘ（Ｍ） （す
なわちΔＸ（Ｌ）の勾配）は負から正に大きく変化する
。第９図は、このような場所を炉心の底から頂上に向っ
て順次探し出し、次式の処理をするための手順である。
ＬＯＣ（Ｎ）＝Ｍ ・・・・・・・・・・
・・・・・（７）ここで、ＬＯＣ（Ｎ）：Ｎ番目の中性
子吸収体の位置の候補 Ｍ：Δ”Ｘ（Ｍ）が負から正に 大きく変化した場所Ｍ ただし、現実にはΔ”Ｘ（Ｎ）は第８図のように変化す
るので、つぎのような処理を付加しである。
・・・・・（７）ここで、ＬＯＣ（Ｎ）：Ｎ番目の中性
子吸収体の位置の候補 Ｍ：Δ”Ｘ（Ｍ）が負から正に 大きく変化した場所Ｍ ただし、現実にはΔ”Ｘ（Ｎ）は第８図のように変化す
るので、つぎのような処理を付加しである。
（１） Δ”Ｘ（Ｍ）が正に変化した場所を拾い出す
と、第８図の点ａ、ｃ、ｅのように、中性子吸収体（１
１） のない位置までＬＯＣ（Ｎ）に選択される。これを除外
するために、拾い出す基準を次式のようにする。
と、第８図の点ａ、ｃ、ｅのように、中性子吸収体（１
１） のない位置までＬＯＣ（Ｎ）に選択される。これを除外
するために、拾い出す基準を次式のようにする。
Δ”Ｘ（Ｍ）＞α ・・・・・・・・・・・
・・・・（８）ここで、αは正の定数である。この値は
、原子炉に応じて（すなわち、炉心の出力密度や中性子
吸収体の仕様に応じて）、適切に設定できる。第８図の
例では、 α＝ １ （Ｗ／ｃｍ”／ （ノード）″）・・・・
・・・・・（９）となる。
・・・・（８）ここで、αは正の定数である。この値は
、原子炉に応じて（すなわち、炉心の出力密度や中性子
吸収体の仕様に応じて）、適切に設定できる。第８図の
例では、 α＝ １ （Ｗ／ｃｍ”／ （ノード）″）・・・・
・・・・・（９）となる。
（１１）第８図の点２と５１点４と６１点５とｆのよう
に、式（８）を満足する点が２個以上連続する場合があ
る。このような場合には、Δ”Ｘ （Ｍ）が大きい方を
Ｌ ＯＣ（Ｎ）として選択する。
に、式（８）を満足する点が２個以上連続する場合があ
る。このような場合には、Δ”Ｘ （Ｍ）が大きい方を
Ｌ ＯＣ（Ｎ）として選択する。
以上の処理によシ、第８図に示すＮ＝１〜１００１０点
が選択できる。第２図と比較すると、中性子吸収体Ｓ１
〜ＳｔＯのすべてに相当する位置が選択されている。
が選択できる。第２図と比較すると、中性子吸収体Ｓ１
〜ＳｔＯのすべてに相当する位置が選択されている。
なお、くぼみを探す処理としては、ΔＲＸ （Ｍ）を使
用する方法の他に、第５図においてＴＩＰ測定（１２） 値をなめらかに結ぶ包絡線を求め、包絡線とＸ（Ｋ）と
の差の大きい個所を、中性子吸収体の位置として選択す
る方法も考えられる。
用する方法の他に、第５図においてＴＩＰ測定（１２） 値をなめらかに結ぶ包絡線を求め、包絡線とＸ（Ｋ）と
の差の大きい個所を、中性子吸収体の位置として選択す
る方法も考えられる。
（４）中性子吸収体の位置決定処理４００このステップ
では、概略位置推定処理３００の結果得られた中性子吸
収体の位置ＬＯＣ（Ｎ）（第８図の例ではＮ＝１〜１０
）の各々について、ＴＩＰ測定値上での正確な中性子吸
収体の位置を決定する。すなわち、Ｌ ＯＣ（Ｎ）の近
傍で、くぼみの底となる位置を探し出す。
では、概略位置推定処理３００の結果得られた中性子吸
収体の位置ＬＯＣ（Ｎ）（第８図の例ではＮ＝１〜１０
）の各々について、ＴＩＰ測定値上での正確な中性子吸
収体の位置を決定する。すなわち、Ｌ ＯＣ（Ｎ）の近
傍で、くぼみの底となる位置を探し出す。
第１０図に示す実施例では、このくぼみの底を探す方法
として、Ｔ工Ｐ測定値の勾配（すなわちΔＸ（Ｌ） ）
を利用している。くぼみの部分が極小値になっていれば
、その点でΔＸ（Ｌ）は０になる。
として、Ｔ工Ｐ測定値の勾配（すなわちΔＸ（Ｌ） ）
を利用している。くぼみの部分が極小値になっていれば
、その点でΔＸ（Ｌ）は０になる。
しかし、第５図から判るように、ＴＩＰ測定値Ｘ （Ｋ
）自身が増加あるいは減少している部分にくぼみが重っ
ている場合には、極小値があられれない。そこで、第１
０図では、以下のような処理をする。
）自身が増加あるいは減少している部分にくぼみが重っ
ている場合には、極小値があられれない。そこで、第１
０図では、以下のような処理をする。
（１） ＬＯＣ（Ｎ）に対応するノードには、次式で
与（１３） えられる。
与（１３） えられる。
Ｋ＝ＬＯＣ（Ｎ）＋１ ・・・・・・・・・・
・・αＯこのノードにの前後で、ＴＩＰ測定値Ｘ（ＩＯ
の平均的な勾配を求める。第８図のＮ＝２の点について
、平均的な勾配を求めた例を第１１図（第５図を部分的
に拡大したもの）に示す。この場合には、ノード）（（
＝ １３ ）の前後に２ノードずつ、計５ノードをとっ
て、最小自乗法によシ直線回帰式（図中の破線）を計算
し、平均的な勾配を求めている。その値は−Ｚ９２であ
る。
・・αＯこのノードにの前後で、ＴＩＰ測定値Ｘ（ＩＯ
の平均的な勾配を求める。第８図のＮ＝２の点について
、平均的な勾配を求めた例を第１１図（第５図を部分的
に拡大したもの）に示す。この場合には、ノード）（（
＝ １３ ）の前後に２ノードずつ、計５ノードをとっ
て、最小自乗法によシ直線回帰式（図中の破線）を計算
し、平均的な勾配を求めている。その値は−Ｚ９２であ
る。
（１１）つぎに、ノードにの近傍で、ΔＸ（Ｌ）が平均
的勾配と一致する場所を求める。第１２図は、Ｋ＝１３
の近傍に関して、第７図を拡大したものである。この場
合、Ｌとしては（Ｋ−２）。
的勾配と一致する場所を求める。第１２図は、Ｋ＝１３
の近傍に関して、第７図を拡大したものである。この場
合、Ｌとしては（Ｋ−２）。
（Ｋ−１） 、 Ｋ、 （Ｋ＋１ ） の４点とる。
この間でΔＸ（Ｌ）を内挿しく図中の太線部分）、β＝
−２，９２との交点を求めると、これがくぼみの底すな
わち中性子吸収体の位置である。一般的に、中性子吸収
体ＳＮについて得た交点をＬＨ（１４） とおくと、炉心の底からの距離（インチ）Ｚ（ＳＮ）は
次式で計算される。
−２，９２との交点を求めると、これがくぼみの底すな
わち中性子吸収体の位置である。一般的に、中性子吸収
体ＳＮについて得た交点をＬＨ（１４） とおくと、炉心の底からの距離（インチ）Ｚ（ＳＮ）は
次式で計算される。
Ｚ（８Ｎ）＝ＬＮ・３ ・・・・・・・・・・・
・・・・αη第１２図の場合（Ｎ＝２）は、Ｌ２＝１２
．６５であシ、これを式αυに代入すると、Ｚ （Ｓ
ｓ ）＝３７．９５となる。
・・・・αη第１２図の場合（Ｎ＝２）は、Ｌ２＝１２
．６５であシ、これを式αυに代入すると、Ｚ （Ｓ
ｓ ）＝３７．９５となる。
（５）ＴＩＰ測定値の補正処理５００
このステップでは、位置決定処理４００でＴＩＰ測定値
Ｘ（Ｋ、）から求めた中性子吸収体８Ｎの位置Ｚ （Ｓ
Ｎ） と、実際の中性子吸収体ＳＮの位置Ｚ （Ｓｗ
） との差に基づいて、ＴＩＰ測定値を補正する。
Ｘ（Ｋ、）から求めた中性子吸収体８Ｎの位置Ｚ （Ｓ
Ｎ） と、実際の中性子吸収体ＳＮの位置Ｚ （Ｓｗ
） との差に基づいて、ＴＩＰ測定値を補正する。
補正したＴＩＰ測定値は、出力分布計算の基準データと
して使用される。
して使用される。
補正処理の実施例として第１３図に示したのは、中性子
吸収体ＳＮの位置のずれ分（Ｚ（ＳＮ）−Ｚ （８Ｎ）
）だけずらして、’ｒＩＰ測定値Ｘ （Ｋ）を読取る
方法である。
吸収体ＳＮの位置のずれ分（Ｚ（ＳＮ）−Ｚ （８Ｎ）
）だけずらして、’ｒＩＰ測定値Ｘ （Ｋ）を読取る
方法である。
（１） ノードにの補正したＴＩＰ測定値Ｘ（Ｋ）を
読取るために、まずノードにの炉心底からの距離Ｚ （
Ｋ）を次式で計算する。
読取るために、まずノードにの炉心底からの距離Ｚ （
Ｋ）を次式で計算する。
（１５）
Ｚ（Ｋ）＝Ｋ・３−１．５ ・・・川・・・Ｑ
路ここで、３はノード１個の長さくインチ）である。ノ
ード中央のＴＩＰ測定値を読取るので、１．５減じであ
る。式（６）で与えられるＺ （Ｋ）は実際の位置であ
り、ＴＩＰ測定値を読取る場合には、中性子吸収体位置
のずれに応じて、Ｚ（Ｋ）ドにの下および上にある最も
近い中性子吸収体ＳｏｌおよびＳＮ＋ｔの炉心底からの
距離を使って、つぎのように求められる。
路ここで、３はノード１個の長さくインチ）である。ノ
ード中央のＴＩＰ測定値を読取るので、１．５減じであ
る。式（６）で与えられるＺ （Ｋ）は実際の位置であ
り、ＴＩＰ測定値を読取る場合には、中性子吸収体位置
のずれに応じて、Ｚ（Ｋ）ドにの下および上にある最も
近い中性子吸収体ＳｏｌおよびＳＮ＋ｔの炉心底からの
距離を使って、つぎのように求められる。
ただし、
ΔＺ ＝ Ｚ （Ｋ） Ｚ （ＳＮ） ・・・・
・・・・・・・・α尋ここで、Ｚ（ＳＮ）−、Ｚ’（Ｓ
Ｎ＋１） ：中性子吸収体の実際の位置 Ｚ（ＳＮ）、 Ｚ（８Ｎ−１） ：中性子吸収体のＴ
ＩＰ測定値上 の位置 （１１）弐〇３で求めたＺ （Ｋ）を用いて、次式にょ
シ補（１６） 正したＴＩＰ測定値マ（Ｋ）を求める。
・・・・・・・・α尋ここで、Ｚ（ＳＮ）−、Ｚ’（Ｓ
Ｎ＋１） ：中性子吸収体の実際の位置 Ｚ（ＳＮ）、 Ｚ（８Ｎ−１） ：中性子吸収体のＴ
ＩＰ測定値上 の位置 （１１）弐〇３で求めたＺ （Ｋ）を用いて、次式にょ
シ補（１６） 正したＴＩＰ測定値マ（Ｋ）を求める。
ただし、
ΔＺ’ ＝Ｚ（Ｋ）−Ｚ（Ｋ’） ・旧旧・・α
・ここで、Ｚ（Ｋ勺 はＺ （Ｋ）より炉心の下側にあ
って最も近いノードに′の実際の位置であり、式（６）
のＫの代わシにに′を代入することによシ得られる。
・ここで、Ｚ（Ｋ勺 はＺ （Ｋ）より炉心の下側にあ
って最も近いノードに′の実際の位置であり、式（６）
のＫの代わシにに′を代入することによシ得られる。
この処理手順を、第２図の例において、ノード１０にお
ける補正したＴＩＰ測定値を求める場合に適用してみる
。まず、式（イ）よシ Ｚ（１０）＝ １０．３−１．５ ＝２８．５（インチ） である。第２図よｐｌこの両側にある中性子吸収体の実
際の位置がつぎのように得られる。
ける補正したＴＩＰ測定値を求める場合に適用してみる
。まず、式（イ）よシ Ｚ（１０）＝ １０．３−１．５ ＝２８．５（インチ） である。第２図よｐｌこの両側にある中性子吸収体の実
際の位置がつぎのように得られる。
Ｚ（Ｓｌ）＝１８．３
Ｚ（Ｓｓ）＝３９．０
これに対して、ＴＩＰ測定値のくぼみから決定した中性
子吸収体の位置は、ｓｇに関しては、前（１７） 項で導出したように、 Ｚ（Ｓｇ） ＝ａ ｓ、。
子吸収体の位置は、ｓｇに関しては、前（１７） 項で導出したように、 Ｚ（Ｓｇ） ＝ａ ｓ、。
である。Ｓｌに関しても、中性子吸収体の概略位置推定
処理３００、位置決定処理４０００手順によって同様に
求めると、つぎのようになる。
処理３００、位置決定処理４０００手順によって同様に
求めると、つぎのようになる。
Ｚ（Ｓｔ） ＝ １８．８
これらを式（至）１０４に代入すると、っぎのようにな
る。
る。
＝２７．８（インチ）
Ｚ（１０）と比べると、炉心の底の方向に０． フイン
チずれている。したがって、式（ト）、α呻において、
Ｋ′＝９となシ、 Ｚ（９）＝９．３−１．５ ＝ ２５．５ さらに、第５図によシ Ｘ（９） ＝ ６１．９２ Ｘ（１０） ＝ ｓｇ、２５ （Ｗ／ｃｍ２）となる。
チずれている。したがって、式（ト）、α呻において、
Ｋ′＝９となシ、 Ｚ（９）＝９．３−１．５ ＝ ２５．５ さらに、第５図によシ Ｘ（９） ＝ ６１．９２ Ｘ（１０） ＝ ｓｇ、２５ （Ｗ／ｃｍ２）となる。
最終的に求めるノード１ｏの補正した（１８）
ＴＩＰ測定値は、つぎのようになる。
＝５９．１１（Ｗ／ｃｒｎ” ）
なお、第１３図の実施例においては、式に）に示すよう
に、平滑化したノード毎のＴＩＰ測定値Ｘ （Ｋ）を内
挿して、補正した測定値Ｘ （Ｋ）を求めている。この
部分は、１インチ毎のＴＩＰ測定値Ｘ （Ｎ）を用いて
、式（１）と同様に、計算し直す方式にしてもよい。
に、平滑化したノード毎のＴＩＰ測定値Ｘ （Ｋ）を内
挿して、補正した測定値Ｘ （Ｋ）を求めている。この
部分は、１インチ毎のＴＩＰ測定値Ｘ （Ｎ）を用いて
、式（１）と同様に、計算し直す方式にしてもよい。
以上説明したように、本発明の処理方法を使用すれば、
より正確なＴＩＰ測定値が、出力分布計算の基準データ
として、利用可能となる。その結果、出力分布の計算精
度が向上し、よシ効率的に沸騰水型動力炉を運転するこ
とができる。
より正確なＴＩＰ測定値が、出力分布計算の基準データ
として、利用可能となる。その結果、出力分布の計算精
度が向上し、よシ効率的に沸騰水型動力炉を運転するこ
とができる。
第１図はＴＩＰ駆動機構を示す図、第２図はＴＩＰ測定
値の一例を示す図、第３図は本発明の処理方法の手順を
示す図、第４．６，９．１０および１３図はそれぞれ、
平滑化処理、準備処理、（１９） 概算位置推定処理、位置決定処理、および補正処理の手
順を示すフロー・シート、第５．７，８゜１１および１
２図は本発明による中性子検出測定データの処理方法の
いろいろな段階における処理されたＴＩＰ測定値と軸方
向位置の間の関係を示すダイヤグラムである。 １・・・原子炉容器、２・・・炉心、３・・・炉心の頂
上、４・・・炉心の底、５・・・ＴＩＰ、６・・・導管
、７・・・ワイヤ、８・・・駆動装置、９・・・Ｘ−Ｙ
レコーダ、１０・・・プロセス計算機、１００・・・平
滑化処理、２００・・・準備処理、３００・・・概略位
置推定処理、４００・・・位置決定処理、５００・・・
補正処理。 代理人 弁理士 秋本正実 （２０） 子４１１ａ ３θθ〜 第７ 目 θ ６ １２ ＩＩ／ ２４
Ｊθ ３１４２ （ｉＢ軸方向位量く）−ド
ラ 早９図 第 １０ 図 茅Ｉｌ 目 ２ 茅Ｉ２図 た 軸方ｔ１体１（ツートノ 早 ７３ 目
値の一例を示す図、第３図は本発明の処理方法の手順を
示す図、第４．６，９．１０および１３図はそれぞれ、
平滑化処理、準備処理、（１９） 概算位置推定処理、位置決定処理、および補正処理の手
順を示すフロー・シート、第５．７，８゜１１および１
２図は本発明による中性子検出測定データの処理方法の
いろいろな段階における処理されたＴＩＰ測定値と軸方
向位置の間の関係を示すダイヤグラムである。 １・・・原子炉容器、２・・・炉心、３・・・炉心の頂
上、４・・・炉心の底、５・・・ＴＩＰ、６・・・導管
、７・・・ワイヤ、８・・・駆動装置、９・・・Ｘ−Ｙ
レコーダ、１０・・・プロセス計算機、１００・・・平
滑化処理、２００・・・準備処理、３００・・・概略位
置推定処理、４００・・・位置決定処理、５００・・・
補正処理。 代理人 弁理士 秋本正実 （２０） 子４１１ａ ３θθ〜 第７ 目 θ ６ １２ ＩＩ／ ２４
Ｊθ ３１４２ （ｉＢ軸方向位量く）−ド
ラ 早９図 第 １０ 図 茅Ｉｌ 目 ２ 茅Ｉ２図 た 軸方ｔ１体１（ツートノ 早 ７３ 目
Claims (1)
- １、走行型中性子検出器の測定値から中性子吸収体の位
置を検出する処理と、実際に中性子吸収体の存在する位
置と測定値から検出された位置との差から測定値に対応
する上記走行型中性子検出器の駆動装置から算出される
炉心高さ方向の位置の誤差を検出し、測定値を補正する
処理とを含むことを特徴とする、原子炉における中性子
検出測定データの処理方法。′
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP57088785A JPS58205894A (ja) | 1982-05-27 | 1982-05-27 | 原子炉における中性子検出器測定データの処理方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP57088785A JPS58205894A (ja) | 1982-05-27 | 1982-05-27 | 原子炉における中性子検出器測定データの処理方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS58205894A true JPS58205894A (ja) | 1983-11-30 |
| JPH0317117B2 JPH0317117B2 (ja) | 1991-03-07 |
Family
ID=13952497
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP57088785A Granted JPS58205894A (ja) | 1982-05-27 | 1982-05-27 | 原子炉における中性子検出器測定データの処理方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS58205894A (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6275294A (ja) * | 1985-09-30 | 1987-04-07 | 株式会社東芝 | 中性子束測定装置 |
| JP2009145332A (ja) * | 2007-12-13 | 2009-07-02 | Global Nuclear Fuel Americas Llc | 原子力プラントの横方向インコアプローブ監視および較正デバイスならびにその方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS53107588A (en) * | 1977-03-01 | 1978-09-19 | Nippon Atom Ind Group Co Ltd | Traveling incoreprobe position correction method and device |
-
1982
- 1982-05-27 JP JP57088785A patent/JPS58205894A/ja active Granted
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS53107588A (en) * | 1977-03-01 | 1978-09-19 | Nippon Atom Ind Group Co Ltd | Traveling incoreprobe position correction method and device |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6275294A (ja) * | 1985-09-30 | 1987-04-07 | 株式会社東芝 | 中性子束測定装置 |
| JP2009145332A (ja) * | 2007-12-13 | 2009-07-02 | Global Nuclear Fuel Americas Llc | 原子力プラントの横方向インコアプローブ監視および較正デバイスならびにその方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0317117B2 (ja) | 1991-03-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6219196B1 (en) | Method and apparatus for suppressing mechanical resonance in a disk drive storage device using a notch filter | |
| US7271977B1 (en) | Method to compensate for microjog error induced by localized track squeeze | |
| EP1408496A2 (en) | Off-track data read retry method for hard disk drive | |
| JPS58205894A (ja) | 原子炉における中性子検出器測定データの処理方法 | |
| JP4521181B2 (ja) | 磁気テープ・ヘッド用の書き込みトラック幅測定装置 | |
| EP0997901A1 (en) | Wave-shaping apparatus and reproduction signal processing apparatus including the same | |
| JP2005243160A (ja) | ディスク装置 | |
| JPS5985930A (ja) | 形状検出装置 | |
| CN116597691A (zh) | 一种航线数据异常值检测与消除的方法及航行监视系统 | |
| CN111580456B (zh) | 一种对小线段连续轨迹进行刀具半径补偿的实现方法 | |
| JPH05296754A (ja) | エッジ検出方法 | |
| JP2005271682A (ja) | トロリ線摩耗量の測定方法及び測定装置 | |
| JP4465857B2 (ja) | コイルエンドの形状判定方法 | |
| JP4905549B2 (ja) | コイルのテレスコープ量算出方法 | |
| US20230307002A1 (en) | Magnetic disk apparatus and method | |
| US8976473B1 (en) | Inter-track interference cancellation based on predetermined data patterns in adjacent tracks | |
| US20250102329A1 (en) | Magnetic sensor system, and distance measuring method for the magnetic sensor system | |
| US20040125496A1 (en) | Methods for setting thresholds for rotatable storage media | |
| JP3100098B2 (ja) | 形状処理方法 | |
| US20140362462A1 (en) | Systems and Methods for Media Defect Detection With Pattern Qualification | |
| JPS5819042B2 (ja) | 位置検出法 | |
| CN121594982A (zh) | 一种超声流量计的增益值调整方法、计算机存储介质及超声流量计 | |
| JPH0116383B2 (ja) | ||
| JPS6316900B2 (ja) | ||
| JPS5897609A (ja) | マ−ク位置検出装置 |