JPS5810717B2 - 原子炉炉心異常診断方法 - Google Patents
原子炉炉心異常診断方法Info
- Publication number
- JPS5810717B2 JPS5810717B2 JP50099530A JP9953075A JPS5810717B2 JP S5810717 B2 JPS5810717 B2 JP S5810717B2 JP 50099530 A JP50099530 A JP 50099530A JP 9953075 A JP9953075 A JP 9953075A JP S5810717 B2 JPS5810717 B2 JP S5810717B2
- Authority
- JP
- Japan
- Prior art keywords
- nuclear reactor
- abnormality
- dispersion range
- spectral density
- predetermined dispersion
- 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
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)
Description
【発明の詳細な説明】
本発明は、運転中の原子炉における炉心内の異常を診断
する方法に関するものである。
する方法に関するものである。
本発明の目的は、異常の原因を精度良く判別できる原子
炉炉心異常診断方法を提供することにある。
炉炉心異常診断方法を提供することにある。
本発明の特徴は、炉心状態を示す複数であるm個の測定
値からなる時系列信号を多変数自己回帰モデルに基づい
てパワースペクトル密度に変換し、パワースペクトル密
度のm個の要素を順次基準スペクトル密度と比較してm
個の要素が所定の分散範囲内に存在するか否かを判定し
、所定の分散範囲から外れている要素がある場合にはそ
の要素の所定の分散範囲から外れている周波数領域を求
め、この周波数領域で所定の分散範囲を外れる他の要素
の存在の有無を判別することにある。
値からなる時系列信号を多変数自己回帰モデルに基づい
てパワースペクトル密度に変換し、パワースペクトル密
度のm個の要素を順次基準スペクトル密度と比較してm
個の要素が所定の分散範囲内に存在するか否かを判定し
、所定の分散範囲から外れている要素がある場合にはそ
の要素の所定の分散範囲から外れている周波数領域を求
め、この周波数領域で所定の分散範囲を外れる他の要素
の存在の有無を判別することにある。
上記の目的を達成するための総括的な説明を述べる。
運転中の原子炉から計測される実測データは、原子炉内
の動的な物理現象の不規則現象の見本過程の結果である
と考えられる。
の動的な物理現象の不規則現象の見本過程の結果である
と考えられる。
計測されたデータ記録が定常正規確率過程にしたがうと
、その確率過程は平均値と相関関数(またはスペクトル
密度)によって完全に記述できることになる。
、その確率過程は平均値と相関関数(またはスペクトル
密度)によって完全に記述できることになる。
このとき任意の確率過程(時系列)は正規性白色雑音を
入力とする定係数線形フィルタの出力として記述できる
。
入力とする定係数線形フィルタの出力として記述できる
。
このため時系列解析を行なって、線形フィルタの特性を
決定することによって、原子炉の動的な特性の推定が可
能になる。
決定することによって、原子炉の動的な特性の推定が可
能になる。
いま原子炉から測定される不規側に変動する各パラメー
タから平均値成分を除去したm次元(原子炉の動的特性
を記述しているm個の計測パラメータに対応)の確率過
程(X(t))の数学的モデルは一般的に次のように表
わされる。
タから平均値成分を除去したm次元(原子炉の動的特性
を記述しているm個の計測パラメータに対応)の確率過
程(X(t))の数学的モデルは一般的に次のように表
わされる。
ここにξ(1)は正規性白色雑音系列であり、m次元ベ
クトルを意味していて、各モーメントはなる性質をもつ
。
クトルを意味していて、各モーメントはなる性質をもつ
。
ここで遅延作用素z−1を用いると式(1)から
で示される多変数の有理型モデルが得られる。
ここでp=qとして
とすると式(2)はまた
と伝達関数でもって書くことができる。
いま式(3)で示されるm次元の伝達関数A(Z)−1
B(Z)および雑音入力の分散Σ2が得られたとすると
X(t)のスペクトル密度は、で与えられる。
B(Z)および雑音入力の分散Σ2が得られたとすると
X(t)のスペクトル密度は、で与えられる。
いまm個の白色雑音系列ξ(1)が互いに独立であると
Σ2は対角行列となり、任意のパラメータ(i)のパワ
ースペクトル度はで表わされる。
Σ2は対角行列となり、任意のパラメータ(i)のパワ
ースペクトル度はで表わされる。
ここでσii2はΣ2の対角項であり、の線形結合で表
現される関係である。
現される関係である。
式(5)に従がって、原子炉の任意の計測パラメータの
パワースペクトルはそれぞれのパラメータの時系列を生
成している白色雑音の分散の線形結合でもって表わされ
る。
パワースペクトルはそれぞれのパラメータの時系列を生
成している白色雑音の分散の線形結合でもって表わされ
る。
別な見方をすると式(5)に従って、原子炉の任意のパ
ラメータのPSDは、m個のそれぞれのパラメータに依
存する割合に分離されることになる。
ラメータのPSDは、m個のそれぞれのパラメータに依
存する割合に分離されることになる。
したがって、任意のパラメータのPSDを求め、ある特
性の周波数帯のパターンが通常のパターンと異なる場合
に、そのパターンこの相異がどのパラメータに依存する
かを識別することが可能となる。
性の周波数帯のパターンが通常のパターンと異なる場合
に、そのパターンこの相異がどのパラメータに依存する
かを識別することが可能となる。
以下本発明を実施例によって詳細に説明する。
第1図は本発明の原子炉異常診断装置の一例である。
原子炉炉心の計測量は信号取り込み装置1でもって、本
異常診断装置に入力される。
異常診断装置に入力される。
ここで計測パラメータをm個選択し、同時に計測上の特
定の雑音をろ液し、アナログ計測信号を時系列に変換す
る。
定の雑音をろ液し、アナログ計測信号を時系列に変換す
る。
時系列信号9は定常性検定装置2によって、適当な時間
間隔で移動平均を求め、1次のモーメントの定常性を検
定し、さらにあらかじめ設定されたサンプル数Nで平均
化を行ない、入力時系列からこの平均値を除いた時系列
を求める。
間隔で移動平均を求め、1次のモーメントの定常性を検
定し、さらにあらかじめ設定されたサンプル数Nで平均
化を行ない、入力時系列からこの平均値を除いた時系列
を求める。
この時系列信号10は統計処理装置3に導かれ、共分散
行列を計算される。
行列を計算される。
得られた共分散行列信号11は定常性検定装置2にもど
されて、2次のモーメントの定常性が検定される。
されて、2次のモーメントの定常性が検定される。
定常な場合には検定された共分散行列信号12は有理型
モデル推定装置4に入力される。
モデル推定装置4に入力される。
いま原子炉の計測信号による時系列がとくに自己回帰モ
デルでもって表わされる場合を例にとると、式(1)は で表わされる。
デルでもって表わされる場合を例にとると、式(1)は で表わされる。
このとき、パラメータA(i)は一般化されたYule
−Walker方程式を解くことになって得られること
が知られている。
−Walker方程式を解くことになって得られること
が知られている。
C(j−k)は遅れ数がj−にの共分散行列であり、信
号12で入力される。
号12で入力される。
第2図は自己回帰モデルの場合の有理型モデル推定装置
4の構成例を示すものであり、式(7)で示される正準
方程式を解く演算装置18と白色雑音時系列分散演算装
置19と次数決定装置20から構成されている。
4の構成例を示すものであり、式(7)で示される正準
方程式を解く演算装置18と白色雑音時系列分散演算装
置19と次数決定装置20から構成されている。
正準方程式演算装置18では次数pに関する繰り返し解
法(P、Wittle。
法(P、Wittle。
Biometrica、50,1963参照〕が利用で
きる。
きる。
またある次数pのもとて演算装置18により得られたパ
ラメータAp1、……、Ap(p)より白色雑音時系列
分散Σ2は となりΣ2は演算装置19で演算される。
ラメータAp1、……、Ap(p)より白色雑音時系列
分散Σ2は となりΣ2は演算装置19で演算される。
ここでC’(i)は分散行列C(i)の転置行列である
。
。
自己回帰モデルのパラメータ決定には例えば次式で示さ
れるFinal Prediction Error(
FPE)〔H。
れるFinal Prediction Error(
FPE)〔H。
Akaike、Ann、of the In5t、 o
f 5tatist。
f 5tatist。
Math、vol、22.No、2.1970参照〕が
最小になるようにすればよい。
最小になるようにすればよい。
ここで1Σ21は演算装置19から与えられる白色雑音
行列の行列式であり、次数決定装置20では、式(9)
によりFPE(p)が遂次増加するpに対して演算され
て、FPE(p)を最小にするpの値が自動的に選ばれ
る。
行列の行列式であり、次数決定装置20では、式(9)
によりFPE(p)が遂次増加するpに対して演算され
て、FPE(p)を最小にするpの値が自動的に選ばれ
る。
このpに対するパラメータA1……、A<p)、および
Σは演算装置18より出力信号13として出力され、ス
ペクトル演算装置5の入力となる。
Σは演算装置18より出力信号13として出力され、ス
ペクトル演算装置5の入力となる。
自己回帰モデルの場合のスペクトル密度は(4)式に応
じて次式から求まる。
じて次式から求まる。
いまm個の白色雑音時系列が独立とすると、任意の計測
時系列iのパワースペクトル密度はで与えられる。
時系列iのパワースペクトル密度はで与えられる。
ここで(A(e−J2πf)−1)iKは行列A(eJ
2πf)−1のfj型要素示し、σ2KKは行列Σ2の
は要素を示す。
2πf)−1のfj型要素示し、σ2KKは行列Σ2の
は要素を示す。
ついで、スペクトル演算装置5では式0υに従って各パ
ワースペクトルを計算する。
ワースペクトルを計算する。
式(11)の左辺の値および右辺の各項は原子炉スペク
トル信号14として異常判定装置6に入力される。
トル信号14として異常判定装置6に入力される。
第3図は異常判定装置6の構成例を示す。
スペクトル信号は異常判定制御装置21に入力される。
該制御装置21では、定期的にあるいは適当な時点で、
基準となるスペクトル密度を基準スペクトル記憶装置2
3に出力する。
基準となるスペクトル密度を基準スペクトル記憶装置2
3に出力する。
さらに周期的、あるいは診断に必要な時点でスペクトル
を密度を診断スペクトル記憶装置22に出力する。
を密度を診断スペクトル記憶装置22に出力する。
スペクトル比較装置24は面記憶装置22.23に記憶
されたスペクトルの比較のための同期信号15を出し、
まず弐(11)で定義されたパワースペクトル密度をそ
れぞれ(m個)順次比較する機能をもつ。
されたスペクトルの比較のための同期信号15を出し、
まず弐(11)で定義されたパワースペクトル密度をそ
れぞれ(m個)順次比較する機能をもつ。
このm個のパワースペクトル密度が、あらかじめ設定し
た分散の範囲内にあるときは、診断した時点での炉心の
時系列は基準となる定常状態と一致して、正常であると
判定する。
た分散の範囲内にあるときは、診断した時点での炉心の
時系列は基準となる定常状態と一致して、正常であると
判定する。
m個のpsdのうち、あるi番目(計測変数1)のps
dが、あらかじめ設定した分散値以上に偏りがでると、
まずアラーム信号16を出し計測変数iのpsdに異常
があることをアラーム表示装置7で指示する。
dが、あらかじめ設定した分散値以上に偏りがでると、
まずアラーム信号16を出し計測変数iのpsdに異常
があることをアラーム表示装置7で指示する。
同時にi番目のpsdのどの周波数帯で、偏りがあるか
を示すために、スペクトルパターン出力信号17により
CRT等の表示装置8に出力する。
を示すために、スペクトルパターン出力信号17により
CRT等の表示装置8に出力する。
また面記憶装置に再度同期信号を出し、i番目のpsd
をさらに詳細に比較するため弐(11)の右辺のm個の
スペクトルを順次面記憶装置から引き出し比較する。
をさらに詳細に比較するため弐(11)の右辺のm個の
スペクトルを順次面記憶装置から引き出し比較する。
この結果例えばm個のうちのl計測変数のパターンが同
じ周波数帯で異常であると判定されるとその結果を、両
スペクトルパターンとともに表示する。
じ周波数帯で異常であると判定されるとその結果を、両
スペクトルパターンとともに表示する。
つぎにi番目のpsdのうち圓式の右辺の1(A(e−
J2πf)−t)1112σ112を比較するための周
期信号を出す。
J2πf)−t)1112σ112を比較するための周
期信号を出す。
このスペクトルが対象となっている周波数帯で偏りがあ
るかどうかを判定する。
るかどうかを判定する。
偏りがあると判定されると、計測チャンネルiの、ある
周波数帯における異常は計測チャンネルlとの関係で生
じるものであり、さらに計測チャンネルl自体(他のチ
ャンネルと相関がない)の異常な現象にもとづいいてい
ることが示される。
周波数帯における異常は計測チャンネルlとの関係で生
じるものであり、さらに計測チャンネルl自体(他のチ
ャンネルと相関がない)の異常な現象にもとづいいてい
ることが示される。
また偏りがないと判定される計測チャンネルiの異常は
、計測チャンネルlから1への情報の伝達径路において
異常が生じていると判断できる。
、計測チャンネルlから1への情報の伝達径路において
異常が生じていると判断できる。
本発明によれば、所定の分散範囲から外れる周波数領域
に基づく判定を行っているので、炉心の異常原因が1つ
または関連している複数である場合を問わず、その異常
原因を明確に判別できる。
に基づく判定を行っているので、炉心の異常原因が1つ
または関連している複数である場合を問わず、その異常
原因を明確に判別できる。
さらに、前述したように周波数領域を考慮して判定して
いるので、原子炉自体に起因する異常原因および計測チ
ャンネル等の計装系自体に起因する異常原因をも区別で
きる。
いるので、原子炉自体に起因する異常原因および計測チ
ャンネル等の計装系自体に起因する異常原因をも区別で
きる。
これは、原子炉自体に基づく異常原因と計装系自体に基
づく異常原因では異常を示す周波数領域が異っているか
らである。
づく異常原因では異常を示す周波数領域が異っているか
らである。
第1図は本発明による原子炉異常診断装置の一実施例の
構成を示すブロック図、第2図は自己回帰の有理型モデ
ル推定装置の構成を示すブロック図、第3図は異常判定
装置の構成を示すブロック図である。
構成を示すブロック図、第2図は自己回帰の有理型モデ
ル推定装置の構成を示すブロック図、第3図は異常判定
装置の構成を示すブロック図である。
Claims (1)
- 1 運転中の原子炉より同時に測定された炉心状態を示
す複数であるm個の測定値からなる時系列信号を多変数
自己回帰モデルに基づいてパワースペクトル密度に変換
し、前記パワースペクトル密度のm個の要素を順次基準
スペクトル密度と比較してm個の前記要素が所定の分散
範囲内に存在するか否かを判定し、前記所定の分散範囲
から外れている前記要素がある場合には、その要素の前
記所定の分散範囲から外れている周波数領域を求め、こ
の周波数領域で前記所定の分散範囲を外れる他の前記要
素の存在の有無を判別する原子炉炉心異常診断方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP50099530A JPS5810717B2 (ja) | 1975-08-18 | 1975-08-18 | 原子炉炉心異常診断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP50099530A JPS5810717B2 (ja) | 1975-08-18 | 1975-08-18 | 原子炉炉心異常診断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPS5224696A JPS5224696A (en) | 1977-02-24 |
JPS5810717B2 true JPS5810717B2 (ja) | 1983-02-26 |
Family
ID=14249766
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP50099530A Expired JPS5810717B2 (ja) | 1975-08-18 | 1975-08-18 | 原子炉炉心異常診断方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JPS5810717B2 (ja) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6163213U (ja) * | 1984-09-29 | 1986-04-28 | ||
JPS61502168A (ja) * | 1984-05-22 | 1986-10-02 | サージカル レーザー テクノロジーズ インコーポレイテッド | 内科および外科用レ−ザ探針 |
JPS625341A (ja) * | 1985-07-02 | 1987-01-12 | オリンパス光学工業株式会社 | レ−ザプロ−ブ |
JPH0279923U (ja) * | 1988-12-09 | 1990-06-20 | ||
US8947760B2 (en) | 2009-04-23 | 2015-02-03 | Ravenbrick Llc | Thermotropic optical shutter incorporating coatable polarizers |
US9116302B2 (en) | 2008-06-19 | 2015-08-25 | Ravenbrick Llc | Optical metapolarizer device |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6010279B2 (ja) * | 1979-01-24 | 1985-03-15 | 日本原子力事業株式会社 | 原子力プラントの診断装置 |
-
1975
- 1975-08-18 JP JP50099530A patent/JPS5810717B2/ja not_active Expired
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS61502168A (ja) * | 1984-05-22 | 1986-10-02 | サージカル レーザー テクノロジーズ インコーポレイテッド | 内科および外科用レ−ザ探針 |
JPH0586225B2 (ja) * | 1984-05-22 | 1993-12-10 | Saajikaru Reezaa Tekunorojiizu | |
JPS6163213U (ja) * | 1984-09-29 | 1986-04-28 | ||
JPS625341A (ja) * | 1985-07-02 | 1987-01-12 | オリンパス光学工業株式会社 | レ−ザプロ−ブ |
JPH0279923U (ja) * | 1988-12-09 | 1990-06-20 | ||
US9116302B2 (en) | 2008-06-19 | 2015-08-25 | Ravenbrick Llc | Optical metapolarizer device |
US8947760B2 (en) | 2009-04-23 | 2015-02-03 | Ravenbrick Llc | Thermotropic optical shutter incorporating coatable polarizers |
Also Published As
Publication number | Publication date |
---|---|
JPS5224696A (en) | 1977-02-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Worden et al. | Novelty detection in a changing environment: regression and interpolation approaches | |
Wise et al. | A theoretical basis for the use of principal component models for monitoring multivariate processes | |
US5251151A (en) | Method and apparatus for diagnosing the state of a machine | |
EP1416348B1 (en) | Methodology for temporal fault event isolation and identification | |
US5586066A (en) | Surveillance of industrial processes with correlated parameters | |
US20140365179A1 (en) | Method and Apparatus for Detecting and Identifying Faults in a Process | |
JP4046309B2 (ja) | プラント監視装置 | |
US6477485B1 (en) | Monitoring system behavior using empirical distributions and cumulative distribution norms | |
Isliker et al. | A test for stationarity: finding parts in time series apt for correlation dimension estimates | |
Chatterjee et al. | Estimation of close sinusoids in colored noise and model discrimination | |
JPS5810717B2 (ja) | 原子炉炉心異常診断方法 | |
KR19990082532A (ko) | 이상 검출 방법 및 이상 검출 시스템 | |
CN111145895A (zh) | 异常数据的检测方法及终端设备 | |
Sugimoto et al. | Stationarity and normality test for biomedical data | |
Kawabata | Test of statistical stability of the electroencephalogram | |
CN113557414A (zh) | 用于设备监测的装置 | |
CN112237433B (zh) | 一种脑电信号异常监测***及方法 | |
CN114090949A (zh) | 一种鲁棒的振动信号特征值计算方法 | |
CN114002574A (zh) | 半导体结构的测试方法 | |
Cao et al. | Deterministic structure in multichannel physiological data | |
JP4049331B2 (ja) | 診断対象物の評価方法および評価装置 | |
CN113254874B (zh) | 一种面向不确定性非平稳工业过程的异常监测方法 | |
JPH04316198A (ja) | プラント異常検出装置 | |
Tatge | Acoustic techniques for machinery diagnostics | |
SU813472A1 (ru) | Анализатор случайных процессов |