JPH08327776A - Core control system of boiling water reactor, and core control method - Google Patents

Core control system of boiling water reactor, and core control method

Info

Publication number
JPH08327776A
JPH08327776A JP7133640A JP13364095A JPH08327776A JP H08327776 A JPH08327776 A JP H08327776A JP 7133640 A JP7133640 A JP 7133640A JP 13364095 A JP13364095 A JP 13364095A JP H08327776 A JPH08327776 A JP H08327776A
Authority
JP
Japan
Prior art keywords
core
boiling water
neutron flux
water reactor
signal
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
JP7133640A
Other languages
Japanese (ja)
Other versions
JP3226003B2 (en
Inventor
Kenichi Takeuchi
健一 竹内
Yasuhiro Masuhara
康博 増原
Hideo Soneda
秀夫 曽根田
Hidetaka Sakamoto
英隆 坂本
Takanori Fukahori
貴憲 深堀
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.)
Hitachi Engineering Co Ltd
Hitachi Ltd
Original Assignee
Hitachi Engineering Co Ltd
Hitachi 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 Hitachi Engineering Co Ltd, Hitachi Ltd filed Critical Hitachi Engineering Co Ltd
Priority to JP13364095A priority Critical patent/JP3226003B2/en
Publication of JPH08327776A publication Critical patent/JPH08327776A/en
Application granted granted Critical
Publication of JP3226003B2 publication Critical patent/JP3226003B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00Energy generation of nuclear origin
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00Energy generation of nuclear origin
    • Y02E30/30Nuclear fission reactors

Landscapes

  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

PURPOSE: To provide a core control system and method for boiling water reactor which can accurately predict the drift component varying with time, and estimate the stability margin during a transient event, and immediately after the end of the event. CONSTITUTION: Provided are a delay filter 7 which takes in neutron flux signal data of an APRM or an LPRM which are monitoring signals of a boiling water reactor core and predicts drift component, a self-correlater 8 for calculating a self-correlating value sequentially based on the prediction, an stability evaluator 22 for calculating a width reduction ratio based on the self-correlation value and evaluating stability margin of the core, and an indication device 11 for displaying the stability margin in a core monitoring device 24. Further a controlling device 23 is provided with an SPI operator 14 receiving a control rod selection and insertion signal from a judging device 13, an SPI device 15 directing insertion based on the signal, a control rod drive mechanism 16 inserting control rod in the core by the direction, and a flow controller 31 controlling the core flow rate.

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【産業上の利用分野】本発明は、沸騰水型原子炉炉心の
制御システムに係り、特に、不安定振動の発生を未然に
防止するのに好適な監視装置を有する沸騰水型原子炉炉
心制御システム及び炉心制御方法に関するものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a boiling water nuclear reactor core control system, and more particularly to a boiling water nuclear reactor core control having a monitoring device suitable for preventing the occurrence of unstable vibrations. The present invention relates to a system and a core control method.

【0002】[0002]

【従来の技術】従来の、沸騰水型原子炉を対象とした安
定性に対する監視装置は、特開昭59-79894号公報に記載
のように、原子炉の中性子束検出器LPRM(Local Pow
er Range Monitor)からの中性子束信号の必要周波数成
分の自己相関関数のパターンを決定する自己相関値演算
装置を用いている。ここで、自己相関関数g(t)は次式
(数3)で定義される。
2. Description of the Related Art A conventional stability monitoring device for a boiling water reactor is disclosed in Japanese Patent Laid-Open No. 59-79894, which is a neutron flux detector LPRM (Local Power) for a nuclear reactor.
er Range Monitor) is used to determine the pattern of the autocorrelation function of the required frequency component of the neutron flux signal. Here, the autocorrelation function g (t) is defined by the following equation (Equation 3).

【0003】[0003]

【数3】 (Equation 3)

【0004】ただし、φは中性子束信号からドリフト成
分を差し引いたもの、tは遅れ時間、Tは時間を表す。
However, φ is the neutron flux signal minus the drift component, t is the delay time, and T is the time.

【0005】また、自己相関関数の関数形は、入力信号
のノイズが白色性である場合には、fを周波数として、
次式(数4)のように表せる。
The function form of the autocorrelation function is such that, when the noise of the input signal is white, f is the frequency,
It can be expressed as the following equation (Equation 4).

【0006】[0006]

【数4】 [Equation 4]

【0007】ここで、aは状態量の減衰定数を表し、自
己相関関数は減衰振動の関数となることがわかる。
Here, it is understood that a represents the damping constant of the state quantity, and the autocorrelation function is a function of damping vibration.

【0008】従来、安定余裕度は、減幅比と呼ばれる指
標を用いて、図3及び次式(数5)に示すように、自己
相関の第一ピークg(0)と一周期後の第二ピークg(1/f)
の比で表される。
Conventionally, the stability margin is calculated by using an index called a width reduction ratio, as shown in FIG. 3 and the following equation (Equation 5), and the first peak g (0) of the autocorrelation and the first peak after one cycle. Two peaks g (1 / f)
It is expressed by the ratio of.

【0009】[0009]

【数5】 (Equation 5)

【0010】減幅比γ>1 ならば 不安定 減幅比γ<1 ならば 安定 図4に、従来の安定性監視装置による実機のプラントデ
ータによる結果を示す。 安定性監視装置の入力信号
は、炉内監視信号である中性子束信号APRM(Average
Power Range Monitor)を使用している。(a)はAP
RM信号の生信号である。(b)はAPRMの変動成分
を示す。これは、(a)のAPRM信号よりドリフト成
分を差し引いたものである。また、(c)は(b)のA
PRM信号の変動分の自己相関値を示したものであり、
中性子束信号の自己相関値は、減衰振動波形になり、
(数3)式の予想とよく一致することがわかる。したがっ
て、上記に説明したように、自己相関の第一ピークと第
二ピークより、安定性の指標である減幅比が求められ
る。
If the width reduction ratio γ> 1, it is unstable. If the width reduction ratio γ <1, it is stable. FIG. 4 shows the result of the plant data of the actual machine by the conventional stability monitoring device. The input signal of the stability monitoring device is the neutron flux signal APRM (Average
Power Range Monitor) is being used. (A) is AP
It is a raw signal of the RM signal. (B) shows the fluctuation component of APRM. This is obtained by subtracting the drift component from the APRM signal in (a). Also, (c) is A in (b)
It shows the autocorrelation value of the fluctuation of the PRM signal,
The autocorrelation value of the neutron flux signal becomes a damped oscillation waveform,
It turns out that it agrees well with the expectation of the expression (3). Therefore, as described above, the width reduction ratio, which is an index of stability, is obtained from the first peak and the second peak of the autocorrelation.

【0011】また、別の方法としては、自己相関より、
減衰係数a、周波数fを求めると、次式(数6)により
減幅比を求めることもできる。
As another method, autocorrelation
When the damping coefficient a and the frequency f are obtained, the reduction ratio can be obtained by the following equation (Equation 6).

【0012】[0012]

【数6】 (Equation 6)

【0013】図5に、通常の運転領域と安定性との関係
を示す。沸騰水型原子炉では、すべての運転領域で安定
となるように設計されているため、運転領域内には、不
安定領域は存在しなく、また、余裕をもって設計されて
いる。
FIG. 5 shows the relationship between the normal operating range and stability. Since the boiling water reactor is designed to be stable in all operating regions, there is no unstable region in the operating region, and it is designed with a margin.

【0014】安定性監視上は、運転領域の内、低流量、
高出力側で安定性の余裕が最小となることがわかる。そ
のため、この領域に対しては、安定性用の監視装置によ
る監視が特に重要となることがわかる。
For stability monitoring, low flow rate within the operating range,
It can be seen that the stability margin is minimized on the high output side. Therefore, it can be seen that the monitoring by the stability monitoring device is particularly important in this area.

【0015】安定性用の監視装置による監視が重要とな
る低流量、高出力での運転は、通常の運転では、起動運
転時、制御棒パターン交換時、ポンプトリップ等の過渡
事象時に行われる可能性がある。そのため、従来は、設
計計算により安定性評価を行い、余裕度を考慮して、起
動運転ライン及び制御棒交換ラインを決定している。ま
た、ポンプトリップ等の過渡事象に対しては、設計計算
を基に余裕を考慮し、選択制御棒の作動領域を決定し、
安定性の余裕が少ない運転領域に入らないようにブロッ
クしている。
The operation at low flow rate and high output, which is important to be monitored by the stability monitoring device, can be performed during normal operation, during start-up operation, control rod pattern replacement, or transient event such as pump trip. There is a nature. Therefore, conventionally, the stability is evaluated by design calculation, and the starting operation line and the control rod exchange line are determined in consideration of the margin. For transient events such as pump trips, consider the margin based on the design calculation and determine the operating area of the selected control rod.
It is blocked so as not to enter the operating area where there is little margin of stability.

【0016】また、図5に示すように、A点(減幅比
a)から不安定度を減少させる(減幅比を小さくする)
方向としては、制御棒(SRI)によって出力低下させ
る方法(B点へ)と、流量を増加させる方法(C点へ)
と、両方を組合わせて行う方法(D点へ)がある。
As shown in FIG. 5, the degree of instability is reduced from point A (width reduction ratio a) (the width reduction ratio is reduced).
As for the direction, a method of decreasing the output by the control rod (SRI) (to the point B) and a method of increasing the flow rate (to the point C)
There is a method (to point D) in which both are combined.

【0017】[0017]

【発明が解決しようとする課題】しかしながら、従来技
術は、設計計算に基づき、保守的に運転ラインを決定す
るため、安定性の余裕は大きく、起動、制御棒パターン
交換等に多くの時間が要するという問題があった。
However, in the prior art, since the operation line is conservatively determined based on the design calculation, there is a large margin of stability, and much time is required for startup, control rod pattern replacement, etc. There was a problem.

【0018】また、ポンプトリップ等の過渡事象におい
ては、中性子束信号が急激に変化するため、従来の安定
性監視装置では、信号の変動量の処理により安定性を評
価するが、急激な信号の変化には追従できないため、誤
評価をする懸念があった。
In a transient event such as a pump trip, the neutron flux signal changes abruptly. Therefore, in the conventional stability monitoring device, the stability is evaluated by processing the fluctuation amount of the signal. There was a concern of making an erroneous evaluation because it cannot follow changes.

【0019】図6に、ポンプトリップ直前、直後の安定
性監視装置による結果を示す。(a)はAPRM信号の
生信号値、(b)はドリフト成分を除去したもの、
(c)は減幅比の結果を示す。ポンプトリップ前は、定
格運転であり、減幅比も小さく安定であることわかる。
ポンプトリップが生じると、流量が低減し、出力も急速
に低下する。その結果、中性子束信号APRMも急激に
低減する。ところが、この急激な過渡事象には、現状の
フィルター(ドリフト成分除去)ではドリフト成分の予測
が追従できないため、(b)に示すように時間遅れによ
る大きなドリフト成分が残り、その結果、状態が落ちつ
いてフィルターによるドリフト成分予測が追いつくま
で、減幅比は不信頼な値を表示する。すなわち、過渡事
象発生後、50秒程度は安定性指標、すなわち減幅比を
表示できないことがわかる。
FIG. 6 shows the results obtained by the stability monitoring device immediately before and after the pump trip. (A) is the raw signal value of the APRM signal, (b) is the one with the drift component removed,
(C) shows the result of the reduction ratio. It can be seen that the pump was rated before the trip, and the reduction ratio was small and stable.
When a pump trip occurs, the flow rate is reduced and the output is reduced rapidly. As a result, the neutron flux signal APRM also sharply decreases. However, since the current filter (drift component removal) cannot follow the drift component prediction with this sudden transient event, a large drift component due to a time delay remains as shown in (b), and as a result, the state stabilizes. The reduction ratio displays an unreliable value until the prediction of the drift component by the filter catches up. That is, it is understood that the stability index, that is, the width reduction ratio cannot be displayed for about 50 seconds after the occurrence of the transient event.

【0020】言い替えれば、従来の安定性監視装置で
は、過渡事象直後は、ドリフト成分が急激に変化するた
め、安定性余裕度を精度よく推定できないという課題が
ある。
In other words, the conventional stability monitoring device has a problem that the stability margin cannot be accurately estimated because the drift component changes rapidly immediately after the transient event.

【0021】本発明の目的は、ポンプトリップ等の過渡
事象が生じても、時事刻々変化するドリフト成分を精度
良く予測でき、過渡事象中及び事象終了直後にも、安定
性余裕度が推定できる沸騰水型原子炉炉心制御システム
及び炉心制御方法を提供することにある。
An object of the present invention is to make it possible to accurately predict a drift component that changes from time to time even when a transient event such as a pump trip occurs, and to estimate the stability margin during and immediately after the transient event. An object of the present invention is to provide a water reactor core control system and core control method.

【0022】[0022]

【課題を解決するための手段】上記目的を達成するため
に、本発明は、沸騰水型原子炉炉心の監視信号であるA
PRM、あるいはLPRMの中性子束信号データ等を用
いて前記炉心の安定余裕度を評価監視する炉心監視装置
と、前記安定余裕度を判定する判定装置と、前記判定結
果に基づき前記炉心の運転を制御する制御装置を有する
沸騰水型原子炉炉心制御システムにおいて、前記炉心監
視装置は、前記中性子束信号データを取り込んでドリフ
ト成分を予測する遅延フィルターと、前記予測に基づき
遂次的に自己相関値を求める自己相関器と、前記自己相
関値に基づき減幅比を算出し、前記炉心の安定余裕度を
評価する安定性評価器を備えることを特徴とする。
In order to achieve the above object, the present invention provides a monitoring signal A for a boiling water reactor core.
A core monitoring device that evaluates and monitors the stability margin of the core using PRM or neutron flux signal data of LPRM, a determination device that determines the stability margin, and controls the operation of the core based on the determination result. In a boiling water reactor core control system having a control device, the core monitoring device is a delay filter that predicts a drift component by capturing the neutron flux signal data, and an autocorrelation value is sequentially obtained based on the prediction. It is characterized by comprising an autocorrelator to be obtained and a stability evaluator which calculates a reduction ratio based on the autocorrelation value and evaluates the stability margin of the core.

【0023】また、本発明の他の特徴は、沸騰水型原子
炉炉心の監視信号であるAPRM、あるいはLPRMの
中性子束信号データ等を用いて、前記炉心の安定余裕度
を評価監視しながら前記安定余裕度を判定し、前記判定
結果に基づき前記炉心の運転を制御する沸騰水型原子炉
炉心制御方法において、遅延フィルターに前記中性子束
信号データを取り込んでドリフト成分を予測し、前記予
測に基づき遂次的に自己相関値を求め、前記自己相関値
に基づき減幅比を算出して前記炉心の安定余裕度を評価
することにより、前記炉心を監視しながら前記安定余裕
度を判定し、判定結果により、制御棒選択挿入信号を受
信し、前記信号に基づき前記制御棒挿入を指示し、前記
炉心に挿入することにより、前記炉心の運転を制御する
ことにある。
Another feature of the present invention is that while monitoring and monitoring the stability margin of the core by using neutron flux signal data of APRM or LPRM which is a monitoring signal of the boiling water reactor core, Stability margin is determined, in the boiling water reactor core control method for controlling the operation of the core based on the determination result, the drift component is predicted by incorporating the neutron flux signal data into the delay filter, based on the prediction Sequentially obtain the autocorrelation value, by calculating the reduction ratio based on the autocorrelation value to evaluate the stability margin of the core, to determine the stability margin while monitoring the core, the determination As a result, the control rod selection insertion signal is received, the control rod insertion is instructed based on the signal, and the control rod is inserted into the core to control the operation of the core.

【0024】また、本発明の他の特徴は、沸騰水型原子
炉炉心の監視信号であるAPRM、あるいはLPRMの
中性子束信号データ等を用いて、前記炉心の安定余裕度
を評価監視しながら前記安定余裕度を判定し、前記判定
結果に基づき前記炉心の運転を制御する沸騰水型原子炉
炉心制御方法において、遅延フィルターに前記中性子束
信号データを取り込んでドリフト成分を予測し、前記予
測に基づき遂次的に自己相関値を求め、前記自己相関値
に基づき減幅比を算出して前記炉心の安定余裕度を評価
することにより、前記炉心を監視しながら前記安定余裕
度を判定し、判定結果により、炉心流量調整信号を受信
し、前記信号に基づきポンプシステムへポンプ回転数変
化の信号を発信し、前記ポンプ回転数を変化させて前記
炉心の流量を調整することにより、前記炉心の運転を制
御することにある。
Further, another feature of the present invention is that the stability margin of the reactor core is evaluated and monitored by using the neutron flux signal data of APRM or LPRM which is the monitoring signal of the boiling water reactor core. Stability margin is determined, in the boiling water reactor core control method for controlling the operation of the core based on the determination result, the drift component is predicted by incorporating the neutron flux signal data into the delay filter, based on the prediction Sequentially obtain the autocorrelation value, by calculating the reduction ratio based on the autocorrelation value to evaluate the stability margin of the core, to determine the stability margin while monitoring the core, the determination As a result, a core flow rate adjustment signal is received, a pump rotation speed change signal is transmitted to the pump system based on the signal, and the pump rotation speed is changed to adjust the core flow rate. The Rukoto is to control the operation of the reactor core.

【0025】[0025]

【作用】本発明によれば、遅延フィルターは、ドリフト
成分の急激な変化に対応するため、時刻T1の前後の時
系列データを用いて、データを取込んでいる時刻T0よ
りτだけ過去の時刻T1(=T0-τ)で、ドリフト成分の
予測を行う。
According to the present invention, since the delay filter responds to the abrupt change of the drift component, the time series data before and after the time T1 is used, and the time τ past the time T0 at which the data is captured is τ. The drift component is predicted at T1 (= T0-τ).

【0026】従来のオンライン安定性監視装置で使用さ
れているドリフト成分除去フィルターでは、図7に示す
ように、データが取り込まれている時刻T0でのドリフ
ト成分yiを予想し、元の信号siよりこのドリフト成
分yiを差し引くことにより、APRM信号の変動成分
xiを評価する。このドリフト成分yiの予想には、次
式(数7)及び(数8)のような、算術平均や、一次遅
れのフィルターが使用される。
In the drift component removal filter used in the conventional online stability monitoring device, as shown in FIG. 7, the drift component yi at the time T0 when the data is captured is predicted, and the original signal si is used. By subtracting this drift component yi, the fluctuation component xi of the APRM signal is evaluated. For the prediction of the drift component yi, an arithmetic mean or a filter of first-order lag, which is expressed by the following equations (7) and (8), is used.

【0027】例1(xiの算術平均値)Example 1 (arithmetic mean value of xi)

【0028】[0028]

【数7】 (Equation 7)

【0029】例2(一次遅れのフィルター)Example 2 (first-order lag filter)

【0030】[0030]

【数8】 (Equation 8)

【0031】(ここで、Cはカット周波数に依存した係
数であるが、サンプルタイム0.1秒で、 カット周波数
を数Hzの時はC=0.95程度) 従来の方法を一般化し、まとめると、次式(数9)のよ
うな式に縮約される。
(Here, C is a coefficient depending on the cut frequency, but C = 0.95 when the cut time is several seconds and the cut frequency is several Hz.) The conventional method is generalized and summarized as follows. It is reduced to an equation like (Equation 9).

【0032】[0032]

【数9】 [Equation 9]

【0033】言い替えると、時刻T0より過去のデータ
により、時刻T0でのドリフト成分を外挿的に評価・予
測している。そのため、xiが大きく変化した場合に
は、T0より過去のデータの荷重が大きいため、時刻T
0より過去のデータの重みが付き、xiの変化に追従で
きないという問題がある。そこで、時刻T0以前のデー
タとT0以降のデータより、その間にある時刻T0での
ドリフト成分を内挿的に予想すると、データの荷重が分
散されため、予測信頼性の向上が期待できる。この方法
が本発明の基本原理である。そのため、ドリフト成分の
予想は、データを取り込んでいる時刻T0より、多少遅
れて評価するため、安定度も同じ時間だけ遅れる。しか
し、遅れ時間は数秒程度であり、安定性の監視上は問題
はない。さらに、急激な過渡時に対しても、ドリフト成
分予想が追従するため、安定余裕度を表示できる。
In other words, the drift component at time T0 is extrapolated and estimated from the data past time T0. Therefore, when xi changes greatly, the weight of the past data is larger than T0, so that time T
There is a problem that the weight of the data past 0 is added and the change of xi cannot be followed. Therefore, when the drift component at the time T0 between them is interpolated from the data before the time T0 and the data after the time T0, the load of the data is dispersed, so that the prediction reliability can be improved. This method is the basic principle of the present invention. Therefore, the prediction of the drift component is evaluated with a slight delay from the time T0 when the data is captured, and the stability is also delayed by the same time. However, the delay time is about several seconds, and there is no problem in monitoring the stability. Further, the drift margin prediction follows even during a sudden transition, so that the stability margin can be displayed.

【0034】また、ドリフト成分予想の一般形は、次式
(数10)のようになるが、一例として、最小自乗法の
フィット式によるドリフト成分予測で説明する。
Further, the general form of the drift component prediction is as in the following formula (Equation 10), but as an example, the drift component prediction by the fit formula of the least square method will be explained.

【0035】[0035]

【数10】 [Equation 10]

【0036】最小自乗法によるフィット式も上記一般式
に含まれ、ドリフト成分を精度よく予測できる。図8
は、過渡事象直後のデータ処理の一例を示したものであ
る。データの左端がデータを取り込んでいる時刻T0と
なるが、最小自乗法のフィット式の左端(時刻T0)での
値は、従来技術と同様な手法となるので、大きな差を生
じる。しかし、データを取り込んでいる時刻T0よりτ
秒(2波長程度、4秒程度)遅れた時刻T1(T1=T0ー
τ)では、ドリフト成分は精度良く予想できることがわ
かる。
The fitting formula by the least squares method is also included in the above general formula, and the drift component can be accurately predicted. FIG.
Shows an example of data processing immediately after the transient event. The left end of the data is the time T0 when the data is taken in, but the value at the left end (time T0) of the least-squares fit expression is the same as that of the conventional technique, and thus a large difference occurs. However, τ from the time T0 when the data is acquired
It can be seen that the drift component can be accurately predicted at time T1 (T1 = T0-τ) delayed by a second (about 2 wavelengths, about 4 seconds).

【0037】したがって、上記(数10)で定義した遅
延フィルターを用いると、急激な過渡事象に対しても、
常時、安定性余裕度を評価できる。
Therefore, by using the delay filter defined in (Equation 10), even for a sudden transient event,
The stability margin can always be evaluated.

【0038】つぎに、遅延フィルターの遅れ時間τにつ
いての検討結果を説明する。図9に遅延フィルターの遅
れ時間τをパラメータとした時に過渡事象直後の計算結
果を示す。図の横軸は遅れ時間τを中性子束信号の基本
周波数fを掛け無次元化したものであり、縦軸は安定性
指標の減幅比を示す。図9より明らかなように、τ=0
の時は精度が悪くなるが、τが大きくなるに従い、減幅
比の値は一定値に近づくことがわかる。すなわち、減幅
比の予測精度向上の観点より、τは基本周波数の逆数1
/f以上に確保する必要があることがわかる。
Next, the examination result of the delay time τ of the delay filter will be described. FIG. 9 shows the calculation result immediately after the transient event when the delay time τ of the delay filter is used as a parameter. The horizontal axis of the figure is a dimensionless product obtained by multiplying the delay time τ by the fundamental frequency f of the neutron flux signal, and the vertical axis shows the reduction ratio of the stability index. As is clear from FIG. 9, τ = 0
It can be seen that the accuracy becomes worse at the time of, but the value of the reduction ratio approaches a constant value as τ becomes larger. That is, from the viewpoint of improving the prediction accuracy of the width reduction ratio, τ is the reciprocal 1 of the fundamental frequency.
It is understood that it is necessary to secure at least / f.

【0039】また、過渡事象として、ポンプトリップ事
象に対する本発明を適用した結果を図10に示す。
(a)はAPRM信号の生信号値、(b)はドリフト成
分を除去したAPRM変動値、(c)は安定性の指標の
減幅比の値を示す。図10から明らかなように、ポンプ
トリップの急激な変化に対して、本発明の炉心監視装置
では、ドリフト成分を精度良く予測するため、(b)の
中性子束の変動成分は精度よく評価でき、その結果、減
幅比は、過渡事象発生中も、過渡事象終了後も、精度よ
い値を評価できる。しかし、データを取り込んでいる時
刻T0に対して数秒遅れた(過去の)時刻T1での結果を
評価することにはなるが、遅れ時間も数秒であり、従来
例での50秒間も表示できないことと比較しても、安定
性評価上問題はない。
FIG. 10 shows the result of applying the present invention to a pump trip event as a transient event.
(A) shows the raw signal value of the APRM signal, (b) shows the APRM fluctuation value with the drift component removed, and (c) shows the value of the reduction ratio of the stability index. As is clear from FIG. 10, in the core monitoring device of the present invention with respect to a sudden change in pump trip, since the drift component is accurately predicted, the fluctuation component of the neutron flux in (b) can be accurately evaluated, As a result, the decrement ratio can be evaluated with high accuracy both during the occurrence of the transient event and after the end of the transient event. However, although the result at the time T1 (past time) delayed by several seconds with respect to the time T0 at which the data is fetched is evaluated, the delay time is also several seconds, and it cannot be displayed for 50 seconds in the conventional example. There is no problem in stability evaluation even when compared with.

【0040】このように、急激な過渡事象に対しても、
時事刻々変化するドリフト成分を精度良く予測でき、過
渡事象中、事象終了直後も、安定性余裕度が推定でき、
常時、安定性余裕度を評価することができる。
Thus, even for a sudden transient event,
The drift component that changes from time to time can be accurately predicted, and the stability margin can be estimated during the transient event and immediately after the event.
The stability margin can always be evaluated.

【0041】[0041]

【実施例】以下、本発明の一実施例に係る沸騰水型原子
炉炉心制御システム及び炉心制御方法を、図を用いて説
明する。
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS A boiling water reactor core control system and core control method according to an embodiment of the present invention will be described below with reference to the drawings.

【0042】図1及び図2に、本発明による第1の実施
例を示す。図1は沸騰水型原子炉炉心制御システムの全
体構成を示し、図2は沸騰水型原子炉炉心制御システム
に備えられている炉心監視装置の構成の詳細を示す。
1 and 2 show a first embodiment according to the present invention. FIG. 1 shows the entire configuration of a boiling water reactor core control system, and FIG. 2 shows the details of the configuration of a core monitoring device provided in the boiling water reactor core control system.

【0043】炉心監視装置24は、沸騰水型原子炉炉心
の監視信号であるAPRM21,LPRM20の中性子
束信号データを取込んでドリフト成分を予測する遅延フ
ィルタ7と、BWRの安定度を評価するための自己相関
器8と、周波数評価器9及び減幅比評価器10を備える
安定性評価器22と、安定性余裕度を表示する表示装置
11で構成されている。
The core monitoring device 24 evaluates the stability of the BWR and the delay filter 7 for predicting the drift component by taking in the neutron flux signal data of the APRM 21 and LPRM 20 which are the monitoring signals of the boiling water reactor core. The auto-correlator 8, the stability evaluator 22 including the frequency evaluator 9 and the width reduction ratio evaluator 10, and the display device 11 for displaying the stability margin.

【0044】遅延フィルタ7は、遅延レジスタ1、乗算
器2、加算器3で構成され、入力信号Xnに対して、
(数10)式の関係を満たす。(数10)式を満たす場
合、評価を行う信号は、入力信号Xnより遅延時間Dだ
け過去のXn−Dについて行われる。(数10)式をX
とYの一次式で表したとすると、次式(数11)のよう
The delay filter 7 is composed of a delay register 1, a multiplier 2 and an adder 3, and with respect to the input signal Xn,
The relationship of Expression (10) is satisfied. When the expression (10) is satisfied, the signal to be evaluated is Xn-D past the input signal Xn by the delay time D. X in the equation (10)
And a linear expression of Y, the following expression (Equation 11) is obtained.

【0045】[0045]

【数11】 [Equation 11]

【0046】となる。It becomes

【0047】入力信号は、遅延レジスタ1によりあるX
i,Xi−1,Xi−2,…,Xkとして記憶される。
この記憶された信号は、乗算器2によりai・Xi,b
i・Yi、加算器3によりΣai・Xi,Σbi・Yi
の信号となる。このΣai・Xi,Σbi・Yiにより
ドリフト成分Yn−Dを算出する。このドリフト成分Y
n−Dの符号を逆転させ、加算器によってXn−Dに加
算することにより入力信号の変動成分のみの信号とな
る。
The input signal is the X signal from the delay register 1.
It is stored as i, Xi-1, Xi-2, ..., Xk.
This stored signal is output by the multiplier 2 to ai · Xi, b.
i · Yi, Σai · Xi, Σbi · Yi by the adder 3
Signal. The drift component Yn-D is calculated from Σai · Xi and Σbi · Yi. This drift component Y
By reversing the sign of n-D and adding it to Xn-D by the adder, only the fluctuation component of the input signal becomes a signal.

【0048】自己相関器8は、乗算器5と積分器6及び
遅延器4で構成されている。乗算器5と積分器6により
構成される演算手段では、自己相関器8への入力信号に
対して、遅れ時間0秒の自己相関値g(0)を出力す
る。また、乗算器5と積分器6及び遅延器4により構成
される演算手段は、それぞれ、遅れ時間△t、2△t、
…、n△t秒の自己相関値g(△t)、g(2△t)、
…、g(n△t)を出力する。
The autocorrelator 8 is composed of a multiplier 5, an integrator 6 and a delay device 4. The computing means composed of the multiplier 5 and the integrator 6 outputs the autocorrelation value g (0) with a delay time of 0 seconds in response to the input signal to the autocorrelator 8. Further, the calculation means composed of the multiplier 5, the integrator 6 and the delay device 4 respectively have delay times Δt, 2Δt,
..., the autocorrelation value of nΔt seconds g (Δt), g (2Δt),
,, g (nΔt) is output.

【0049】この出力は、安定評価器22の周波数評価
器9に入力され、周波数f及び減衰係数aの評価の後、
減幅比評価器10によって、減幅比を算出する。算出さ
れた減幅比は、表示装置11により可視化される。
This output is input to the frequency evaluator 9 of the stability evaluator 22, and after evaluation of the frequency f and the attenuation coefficient a,
The reduction ratio evaluator 10 calculates the reduction ratio. The calculated reduction ratio is visualized by the display device 11.

【0050】このように、炉心監視装置24は、過渡事
象においても、時々刻々変化するドリフト成分を精度よ
く予測でき、過渡事象中、事象整定後に精度よく安定性
余裕度を監視することが可能となる。
As described above, the core monitoring device 24 can accurately predict a drift component that changes from moment to moment even in a transient event, and can accurately monitor the stability margin during the transient event after the event settling. Become.

【0051】また、定常時においても、判定装置13に
より安定性余裕度が少なくなったと判断した場合は、制
御装置23に含まれ、制御棒を挿入する一連のSRIや
ジェットポンプやインターナルポンプなど炉心流量を確
保するポンプシステム30を制御する流量調整装置31
によって、安定性余裕度を確保することができる。
Further, even in the steady state, when the determination device 13 determines that the stability margin has decreased, the control device 23 includes a series of SRI for inserting the control rod, a jet pump, an internal pump, etc. Flow rate adjusting device 31 for controlling the pump system 30 for ensuring the core flow rate
Therefore, the stability margin can be secured.

【0052】制御装置23は、選択制御棒挿入信号を判
定装置13から受信するSRI(Selected Rod Inserti
on)作動装置14と、信号に基づき挿入指示をするSR
I装置15と、指示により炉心へ制御棒を挿入する制御
棒駆動機構16及び判定装置13から炉心流量調整信号
を受信し、ポンプシステム30へポンプ回転数変化の信
号を発信し、ポンプの回転数を変化させて炉心の流量を
調整する流量調整装置31で構成されている。
The control device 23 receives the selected control rod insertion signal from the determination device 13 and selects an SRI (Selected Rod Inserti).
on) Actuating device 14 and SR for instructing insertion based on a signal
I device 15, a control rod drive mechanism 16 for inserting a control rod into the core according to an instruction, and a determination device 13 receive a core flow rate adjustment signal, and send a pump rotation speed change signal to the pump system 30 to change the pump rotation speed. Is configured by a flow rate adjusting device 31 for adjusting the flow rate of the core by changing

【0053】図11に、本発明による第2の実施例を示
す。本実施例は、原子炉の起動を、炉心監視装置24を
用いて行うものである。本発明の炉心監視装置24を用
いていない原子炉においては、図11の(a)、(b)
に示すように、予め炉心熱出力と炉心流量で保守的に設
定された安定性制限曲線を避けて起動を行う必要があ
る。本実施例では、リアルタイムに直接安定性余裕度を
監視できるため、保守的に設定された安定性制限曲線に
おける過剰な余裕度をとることなく、安定性余裕度が大
きい場合には、図11の(c)に示すように、若干、安
定性制限曲線を超える起動を行うことができる。これに
より、起動操作の簡略化、起動時間の短縮が可能とな
る。
FIG. 11 shows a second embodiment according to the present invention. In this embodiment, the reactor is started by using the core monitoring device 24. In a nuclear reactor which does not use the core monitoring device 24 of the present invention, (a) and (b) of FIG.
As shown in, it is necessary to avoid the stability limit curve that is conservatively set in advance by the core heat output and the core flow rate and perform the startup. In this embodiment, since the stability margin can be directly monitored in real time, if the stability margin is large without taking an excessive margin in the stability limit curve set conservatively, the stability margin shown in FIG. As shown in (c), it is possible to start the vehicle slightly beyond the stability limit curve. This makes it possible to simplify the startup operation and shorten the startup time.

【0054】図12に、本発明による第3の実施例を示
す。本実施例は、原子炉の制御棒パターン交換を、炉心
監視装置24を用いて行うものである。本発明の炉心監
視装置24を用いていない原子炉においては、図12の
(a)、(b)に示すように、予め炉心熱出力と炉心流
量で保守的に設定された安定性制限曲線を避けて制御棒
パターン交換を行う必要がある。本実施例では、リアル
タイムに直接安定性余裕度を監視できるため、保守的に
設定された安定性制限曲線における過剰な余裕度をとる
ことなく、安定性余裕度が大きい場合には図12の
(c)に示すように、若干、安定性制限曲線を超える制
御棒パターン交換を行うことができる。これにより、制
御棒パターン交換操作の簡略化、制御棒パターン交換時
間の短縮化が可能となる。
FIG. 12 shows a third embodiment according to the present invention. In this embodiment, the control rod pattern of the nuclear reactor is exchanged by using the core monitoring device 24. In a nuclear reactor that does not use the core monitoring device 24 of the present invention, as shown in FIGS. 12 (a) and 12 (b), a stability limiting curve that is conservatively set in advance by the core heat output and the core flow rate is used. It is necessary to avoid it and replace the control rod pattern. In the present embodiment, since the stability margin can be directly monitored in real time, if the stability margin is large without taking an excessive margin in the stability limit curve that is set conservatively, (( As shown in c), a control rod pattern exchange that slightly exceeds the stability limit curve can be performed. As a result, the control rod pattern exchange operation can be simplified and the control rod pattern exchange time can be shortened.

【0055】図13に、本発明による第4の実施例を示
す。本実施例の原子炉炉心制御システムは、炉心監視装
置24と、安定性モード判定器17と、不安定発生位置
判別器18と、不安定発生領域判別器19と、制御装置
23で構成されている。また、炉心監視装置24は、並
列に配置された遅延フィルタ7、自己相関器8、安定性
評価器22と1個の表示装置で構成されている。
FIG. 13 shows a fourth embodiment according to the present invention. The reactor core control system according to the present embodiment includes a core monitoring device 24, a stability mode determiner 17, an instability occurrence position determiner 18, an unstable occurrence region determiner 19, and a controller 23. There is. The core monitoring device 24 is composed of a delay filter 7, an autocorrelator 8, a stability evaluator 22 and one display device arranged in parallel.

【0056】本実施例においては、安定性余裕度が少な
くなった場合、安定性モード判定器17で、安定性のモ
ード(全炉心不安定モード、領域不安定モード、チャン
ネル不安定モード)を判断し適切な安定性対策をガイド
するものである。LPRM20、APRM21の信号か
ら安定性モード判定器17において判別された安定性の
モードに対して、全炉心不安定モードに対しては炉心全
体の出力を低下させるように炉心へ均等に制御棒の挿
入、領域不安定モードに対しては安定性余裕度が少なく
なっている領域に制御棒の挿入、チャンネル不安定モー
ドに対しては安定性余裕度が少なくなっているチャンネ
ルに隣接する制御棒の挿入を行うものである。
In the present embodiment, when the stability margin decreases, the stability mode determiner 17 determines the stability mode (all core unstable mode, region unstable mode, channel unstable mode). It guides appropriate stability measures. For the stability mode determined by the stability mode determiner 17 from the signals of the LPRM20 and the APRM21, the control rods are evenly inserted into the core so as to reduce the output of the entire core for the all core unstable mode. , Insert the control rod into the region where the stability margin is small for the region unstable mode, and insert the control rod adjacent to the channel where the stability margin is small for the channel unstable mode. Is to do.

【0057】また、本実施例によれば、安定性余裕度が
少なくなった場合に、適切な制御棒を挿入することがで
き、安定性余裕度をとるために、過剰な制御棒挿入を行
うことがなくなる。これにより、制御棒挿入後安定性余
裕度が確保された時点からの復旧操作の簡略化、復旧時
間の短縮化が可能となる。
Further, according to the present embodiment, when the stability margin decreases, an appropriate control rod can be inserted, and excessive control rod insertion is performed in order to secure the stability margin. Will disappear. As a result, it becomes possible to simplify the recovery operation from the time when the stability margin after the control rod is inserted and to shorten the recovery time.

【0058】図14に、本発明による第5の実施例を示
す。本実施例の炉心監視装置24は、遅延フィルタ7、
自己相関器8、周波数評価器9、積分範囲評価器12、
減幅比評価器10、表示装置11で構成される。また、
本実施例における動作原理は、基本的には図1の実施例
と同様であるが、遅延レジスタ1内において、ドリフト
成分を規定するために用いる信号データ数を、周波数評
価器9において算出された周波数fを用いて決定する。
これは、(数10)式におけるΣを行うデータ数(積分
範囲)を規定するものであり、周波数fから数波長分の
時間に相当するデータ数で規定する。この様な範囲のデ
ータで規定されるドリフト成分は精度良く求まり、ま
た、評価するデータの範囲を制限することから、処理時
間の短縮も可能となる。
FIG. 14 shows a fifth embodiment according to the present invention. The core monitoring device 24 of the present embodiment includes the delay filter 7,
Auto-correlator 8, frequency evaluator 9, integration range evaluator 12,
The reduction ratio evaluator 10 and the display device 11 are included. Also,
The operating principle of this embodiment is basically the same as that of the embodiment of FIG. 1, but the number of signal data used for defining the drift component in the delay register 1 is calculated by the frequency evaluator 9. Determined using frequency f.
This defines the number of data (integration range) for performing Σ in the equation (10), and is defined by the number of data corresponding to the time of several wavelengths from the frequency f. The drift component defined by the data in such a range can be accurately obtained, and the range of the data to be evaluated is limited, so that the processing time can be shortened.

【0059】図15に、本発明による第6の実施例を示
す。本実施例の原子炉炉心制御システムは、炉心監視装
置24と、判定装置13と、制御装置23で構成されて
いる。また、炉心監視装置24は、遅延フィルタ7、自
己相関器8、周波数評価器9、積分範囲評価器12、減
幅比評価器10、表示装置11で構成されている。
FIG. 15 shows a sixth embodiment according to the present invention. The nuclear reactor core control system according to the present embodiment includes a reactor core monitoring device 24, a determination device 13, and a control device 23. The core monitoring device 24 is composed of a delay filter 7, an autocorrelator 8, a frequency evaluator 9, an integration range evaluator 12, a reduction ratio evaluator 10, and a display device 11.

【0060】本実施例における減幅比算出までの動作原
理は、図1の実施例と同様であるが、本実施例では、算
出された減幅比を判定装置13で基準値と比較し、安定
性余裕が少ないと判断した場合は、制御装置23のSR
I作動装置14に信号を送る。SRI作動装置14に信
号が送られると、SRI作動のための信号がSRI装置
15に送られ、制御棒駆動機構16によりある決められ
た制御棒パターンに従って、制御棒が挿入される。制御
棒が挿入されると出力が低下し、これにより安定性余裕
度が回復する。また、本実施例の信号処理手法を用いる
ことで、SRI作動信号自体の信頼性も高くなる。
The principle of operation up to the calculation of the width reduction ratio in this embodiment is the same as that of the embodiment of FIG. 1, but in the present embodiment, the calculated width reduction ratio is compared with the reference value by the judging device 13, When it is determined that the stability margin is small, the SR of the control device 23
Send a signal to the I-actuator 14. When a signal is sent to the SRI actuating device 14, a signal for SRI actuation is sent to the SRI device 15, and the control rod driving mechanism 16 inserts the control rod according to a predetermined control rod pattern. When the control rod is inserted, the output decreases, which restores the stability margin. Further, by using the signal processing method of the present embodiment, the reliability of the SRI actuation signal itself is increased.

【0061】図16に、本発明による第7の実施例を示
す。従来の処理手法と本発明の処理手法を同時に表示し
て監視する例である。本発明の処理手法を用いる場合は
精度良く処理することが可能となるが、処理に遅れ時間
が伴う。(a)に示すように、従来の処理手法は、応答
は早いため、算出される減幅比を予測値として表示し、
また、(b)に示すように、本発明の処理手法は、算出
される減幅比を確定値として表示するものである。この
ように表示することにより、安定性余裕度を精度よく監
視でき、遅れ時間に伴う安定性余裕不監視時間を予測値
として概略的に把握することによって、安定性余裕度監
視の信頼性を向上することができる。
FIG. 16 shows a seventh embodiment according to the present invention. It is an example of simultaneously displaying and monitoring the conventional processing method and the processing method of the present invention. When the processing method of the present invention is used, the processing can be performed with high accuracy, but the processing involves a delay time. As shown in (a), since the conventional processing method has a quick response, the calculated reduction ratio is displayed as a predicted value,
Further, as shown in (b), the processing method of the present invention displays the calculated reduction ratio as a fixed value. By displaying in this way, the stability margin can be accurately monitored, and the stability margin non-monitoring time associated with the delay time can be roughly understood as a predicted value to improve the reliability of stability margin monitoring. can do.

【0062】図17に、本発明による第8の実施例を示
す。本実施例は、従来のフィルタ処理手法と本発明のフ
ィルタ処理手法を組合せた例である。この方法は、従来
の手法を用いる処理は予測値、本発明の手法を用いる処
理は確定値として表示を行うものである。本発明の手法
を用いる処理は、精度良く処理することが可能となる
が、遅れ時間の間は減幅比の算出ができなくなるため、
従来の手法の値を予測値として表示し、遅れ時間以降
は、本発明の手法を用いて算出した値を確定値として表
示するものである。本実施例により、安定性余裕度を精
度よく監視でき、遅れ時間に伴う安定性余裕不監視時間
を予測値として概略的に把握することによって、安定性
余裕度監視の信頼性を向上することができる。
FIG. 17 shows an eighth embodiment according to the present invention. The present embodiment is an example in which the conventional filter processing method and the filter processing method of the present invention are combined. In this method, a process using the conventional method is displayed as a predicted value, and a process using the method of the present invention is displayed as a fixed value. Although the processing using the method of the present invention can be processed with high accuracy, the width reduction ratio cannot be calculated during the delay time.
The value of the conventional method is displayed as the predicted value, and after the delay time, the value calculated using the method of the present invention is displayed as the fixed value. According to the present embodiment, the stability margin can be accurately monitored, and the stability margin non-monitoring time associated with the delay time can be roughly grasped as a predicted value to improve the reliability of the stability margin monitoring. it can.

【0063】図18に、本発明による第9の実施例を示
す。従来、過渡事象の発生の判断は、中性子束の絶対値
の変化量に基づくものであり、そのため、中性子束の変
化レベルがある設定値に達するまでは時間を要する懸念
がある。本実施例は、遅延フィルター7、信号変換器2
8、自己相関器8、安定性評価器22、表示装置11を
備えた炉心監視装置24において、遅延フィルター7の
下流側に、τ秒過去の時刻T1で急激なドリフト変化が
発生したか否かを判定するドリフト判定器27を設け、
ドリフト判定器27の下流側を選択制御棒駆動機構29
あるいは警報器に連結させる構成としたものである。
FIG. 18 shows a ninth embodiment according to the present invention. Conventionally, the determination of the occurrence of a transient event is based on the amount of change in the absolute value of the neutron flux, so there is a concern that it will take time before the level of change in the neutron flux reaches a certain set value. In this embodiment, the delay filter 7 and the signal converter 2 are used.
8. In the core monitoring device 24 including the autocorrelator 8, the stability evaluator 22, and the display device 11, whether or not a drastic drift change occurs at the time T1 τ seconds past on the downstream side of the delay filter 7. A drift determiner 27 for determining
The downstream side of the drift determiner 27 is a selection control rod drive mechanism 29.
Alternatively, it is configured to be connected to an alarm device.

【0064】かかる構成にしたことにより、急激なドリ
フト成分変化時に、遅延フィルター7の遅れ時間のみの
遅れで、選択制御棒駆動機構29あるいは警報器を作動
させることができる。遅延フィルター7の遅れ時間は、
中性子束の固有振動モードの1周期程度である。中性子
束の固有振動モードの周期は2秒程度であるので、2秒
程度で、過渡事象発生の判定ができる。そのため、十分
に余裕を保持しながら、選択制御棒の作動、警報を早期
に発するので、安全性余裕度がさらに増加する。
With this configuration, when the drift component changes abruptly, the selective control rod drive mechanism 29 or the alarm device can be activated with a delay of only the delay time of the delay filter 7. The delay time of the delay filter 7 is
It is about one cycle of the natural vibration mode of the neutron flux. Since the cycle of the natural vibration mode of the neutron flux is about 2 seconds, the occurrence of a transient event can be determined in about 2 seconds. Therefore, the operation of the selection control rod and the alarm are issued at an early stage while maintaining a sufficient margin, which further increases the safety margin.

【0065】図19に、本発明による第10の実施例を
示す。本実施例は、遅延フィルター7、信号変換器2
8、自己相関器8、安定性評価器22、表示装置11を
備えた炉心監視装置24において、遅延フィルター7と
自己相関器8の間に、τ秒過去の時刻T1で急激なドリ
フト変化の否かを判定するドリフト判定器27と、急激
なドリフト成分変化時にはドリフト成分を現信号成分と
置き換え、信号処理成分を零とみなす信号変換器28を
設けたものである。かかる構成にすると、いかなる高速
の過渡事象に対しても、ドリフト成分は追従できるの
で、ドリフト成分予測の遅れによる精度低下はなくな
る。
FIG. 19 shows a tenth embodiment of the present invention. In this embodiment, the delay filter 7 and the signal converter 2 are used.
8. In the core monitoring device 24 including the autocorrelator 8, the stability evaluator 22, and the display device 11, between the delay filter 7 and the autocorrelator 8, there is a rapid drift change at time T1 τ seconds past. It is provided with a drift determiner 27 for determining whether or not, and a signal converter 28 for replacing the drift component with the current signal component when the drift component changes rapidly and regarding the signal processing component as zero. With this configuration, the drift component can follow any high-speed transient event, so that the accuracy is not reduced due to the delay in drift component prediction.

【0066】[0066]

【発明の効果】本発明によれば、安定性上重要な事象に
おいても、安定性余裕度の監視、運転制御ができ、その
信頼性を向上することができる。
According to the present invention, the stability margin can be monitored and the operation can be controlled even in an event that is important for stability, and its reliability can be improved.

【図面の簡単な説明】[Brief description of drawings]

【図1】本発明の一実施例に係る沸騰水型原子炉炉心制
御システムの全体構成図である。
FIG. 1 is an overall configuration diagram of a boiling water reactor core control system according to an embodiment of the present invention.

【図2】図1の沸騰水型原子炉炉心制御システムに備え
られている炉心監視装置の構成の詳細図である。
FIG. 2 is a detailed view of the configuration of a core monitoring device provided in the boiling water reactor core control system of FIG.

【図3】自己相関値を示す図である。FIG. 3 is a diagram showing an autocorrelation value.

【図4】従来の安定性監視装置による実機評価結果を示
す図である。
FIG. 4 is a diagram showing an evaluation result of an actual machine by a conventional stability monitoring device.

【図5】運転領域と安定性との関係を示す図である。FIG. 5 is a diagram showing a relationship between an operating region and stability.

【図6】従来の安定性監視装置によるポンプトリップ時
の評価結果を示す図である。
FIG. 6 is a diagram showing an evaluation result at the time of pump trip by a conventional stability monitoring device.

【図7】従来のドリフト成分予測方法での結果を示す図
である。
FIG. 7 is a diagram showing a result of a conventional drift component prediction method.

【図8】本発明の方法によるポンプトリップ時のドリフ
ト成分予測結果を示す図である。
FIG. 8 is a diagram showing a drift component prediction result at the time of pump trip by the method of the present invention.

【図9】遅延フィルターの遅れ時間の感度解析を示す図
である。
FIG. 9 is a diagram showing a sensitivity analysis of delay time of a delay filter.

【図10】本発明の方法によるポンプトリップ時の評価
結果を示す図である。
FIG. 10 is a diagram showing an evaluation result at the time of pump trip according to the method of the present invention.

【図11】原子炉の起動を、炉心監視装置を用いて行う
第2の実施例を示す図である。
FIG. 11 is a diagram showing a second embodiment in which the reactor is started up using a core monitoring device.

【図12】原子炉の制御棒パターン交換を、炉心監視装
置を用いて行う第3の実施例を示す図である。
FIG. 12 is a diagram showing a third embodiment in which the control rod pattern of the nuclear reactor is exchanged using a core monitoring device.

【図13】本発明による第4の実施例を示す図である。FIG. 13 is a diagram showing a fourth embodiment according to the present invention.

【図14】本発明による第5の実施例を示す図である。FIG. 14 is a diagram showing a fifth embodiment according to the present invention.

【図15】本発明とSRIを組み合わせた第6の実施例
を示す図である。
FIG. 15 is a diagram showing a sixth embodiment in which the present invention and SRI are combined.

【図16】従来の処理手法と本発明の処理手法を同時に
表示して監視する第7の実施例を示す図である。
FIG. 16 is a diagram showing a seventh embodiment for simultaneously displaying and monitoring the conventional processing method and the processing method of the present invention.

【図17】従来のフィルタ処理手法と本発明のフィルタ
処理手法を組合せた第8の実施例を示す図である。
FIG. 17 is a diagram showing an eighth embodiment in which the conventional filter processing method and the filter processing method of the present invention are combined.

【図18】本発明による第9の実施例を示す図である。FIG. 18 is a diagram showing a ninth embodiment according to the present invention.

【図19】本発明による第10の実施例を示す図であ
る。
FIG. 19 is a diagram showing a tenth embodiment according to the present invention.

【符号の説明】[Explanation of symbols]

1…遅延レジスタ、2,5…乗算器、3…加算器、4…
遅延器、6…積分器、7…遅延フィルター、8…自己相
関器、9…周波数評価器、10…減幅比評価器、11…
表示装置、12…積分範囲評価器、13…判定装置、1
4…SRI作動装置、15…SRI装置、16…制御棒
駆動機構、17…安定性モード判定器、18…不安定発
生位置判別器、19…不安定発生領域判別器、20…L
PRM、21…APRM、22…安定性評価器、23…
制御装置、24…炉心監視装置、25…符号逆転器、2
7…ドリフト判定器、28…信号変換器、29…選択制
御棒駆動装置、30…ポンプシステム、31…流量調整
装置
1 ... Delay register, 2, 5 ... Multiplier, 3 ... Adder, 4 ...
Delay device, 6 ... Integrator, 7 ... Delay filter, 8 ... Autocorrelator, 9 ... Frequency evaluator, 10 ... Reduction ratio evaluator, 11 ...
Display device, 12 ... Integration range evaluator, 13 ... Judgment device, 1
4 ... SRI actuating device, 15 ... SRI device, 16 ... control rod drive mechanism, 17 ... stability mode determiner, 18 ... unstable occurrence position discriminator, 19 ... unstable occurrence region discriminator, 20 ... L
PRM, 21 ... APRM, 22 ... Stability evaluator, 23 ...
Control device, 24 ... Core monitoring device, 25 ... Sign reversal device, 2
7 ... Drift judgment device, 28 ... Signal converter, 29 ... Selection control rod drive device, 30 ... Pump system, 31 ... Flow control device

───────────────────────────────────────────────────── フロントページの続き (72)発明者 曽根田 秀夫 茨城県日立市幸町三丁目1番1号 株式会 社日立製作所日立工場内 (72)発明者 坂本 英隆 茨城県日立市幸町三丁目1番1号 株式会 社日立製作所日立工場内 (72)発明者 深堀 貴憲 茨城県日立市幸町三丁目2番1号 日立エ ンジニアリング株式会社内 ─────────────────────────────────────────────────── ─── Continuation of the front page (72) Hideo Soneda 3-1-1, Saiwai-cho, Hitachi-shi, Ibaraki Hitachi Ltd. Hitachi factory (72) Inventor Hidetaka Sakamoto 3-chome, Saiwai-cho, Hitachi, Ibaraki No. 1 Stock company Hitachi Ltd. Hitachi factory (72) Inventor Takanori Fukahori 32-1 No. 2 Sachimachi, Hitachi City, Ibaraki Hitachi Engineering Co., Ltd.

Claims (16)

【特許請求の範囲】[Claims] 【請求項1】沸騰水型原子炉炉心の監視信号であるAP
RM、あるいはLPRMの中性子束信号データ等を用い
て前記炉心の安定余裕度を評価監視する炉心監視装置
と、前記安定余裕度を判定する判定装置と、前記判定結
果に基づき前記炉心の運転を制御する制御装置を有する
沸騰水型原子炉炉心制御システムにおいて、 前記炉心監視装置は、前記中性子束信号データを取り込
んでドリフト成分を予測する遅延フィルターと、前記予
測に基づき遂次的に自己相関値を求める自己相関器と、
前記自己相関値に基づき減幅比を算出し、前記炉心の安
定余裕度を評価する安定性評価器を備えることを特徴と
する沸騰水型原子炉炉心制御システム。
1. An AP which is a supervisory signal for a boiling water reactor core.
A core monitoring device that evaluates and monitors the stability margin of the core by using neutron flux signal data of RM or LPRM, a determination device that determines the stability margin, and controls the operation of the core based on the determination result. In a boiling water reactor core control system having a control device, the core monitoring device is a delay filter that predicts a drift component by capturing the neutron flux signal data, and an autocorrelation value is sequentially obtained based on the prediction. The desired autocorrelator,
A boiling water reactor core control system, comprising a stability evaluator that calculates a reduction ratio based on the autocorrelation value and evaluates a stability margin of the core.
【請求項2】請求項1において、前記遅延フィルター
は、現在の時刻T0よりτ秒過去の時刻T1(T1=T0
ーτ)で前記ドリフト成分を予測し、前記τ秒過去の時刻
T1での前記中性子束信号より前記ドリフト成分を差し
引くことにより前記中性子束信号の変動量を評価するこ
とを特徴とする沸騰水型原子炉炉心制御システム。
2. The delay filter according to claim 1, wherein the delay filter has a time T1 (T1 = T0) τ seconds past the present time T0.
-Τ) is used to predict the drift component, and the fluctuation amount of the neutron flux signal is evaluated by subtracting the drift component from the neutron flux signal at the time T1 in the past τ seconds. Reactor core control system.
【請求項3】請求項1または請求項2において、前記遅
延フィルターでのτ(過去へのさかのぼり時間)を前記
中性子束信号変動分の基本周波数の逆数1/f以上にし
たことを特徴とする沸騰水型原子炉炉心制御システム。
3. The τ according to claim 1 or 2, wherein the delay filter has a reciprocal time 1 / f or more of the fundamental frequency of the fluctuation of the neutron flux signal. Boiling water reactor core control system.
【請求項4】請求項1ないし請求項3において、前記遅
延フィルター内での前記中性子束信号のドリフト成分予
測を次式(数1)で行うことを特徴とする沸騰水型原子
炉炉心制御システム。 【数1】 ここで、Yk; 時刻Tkでの中性子束信号のドリフト
成分値 Xi; 時刻Tiでの中性子束信号データ値 ai; 係数、 bi; 係数 である。ただし、1<k<mとなり、中性子束信号デー
タの間隔をΔtとすると、(k−m)=τ/Δtとな
る。
4. The boiling water reactor core control system according to claim 1, wherein drift component prediction of the neutron flux signal in the delay filter is performed by the following equation (Equation 1). . [Equation 1] Here, Yk; drift component value of neutron flux signal at time Tk, Xi; neutron flux signal data value at time Ti, ai; coefficient, bi; coefficient. However, if 1 <k <m, and the interval of the neutron flux signal data is Δt, then (k−m) = τ / Δt.
【請求項5】請求項1において、前記制御装置は、選択
制御棒挿入信号を前記判定装置から受信するSRI作動
装置と、前記信号に基づき挿入指示をするSRI装置
と、前記指示により前記炉心へ前記制御棒を挿入する制
御棒駆動機構を備えることを特徴とする沸騰水型原子炉
炉心制御システム。
5. The control device according to claim 1, wherein the control device receives an SRI actuation device for receiving a selection control rod insertion signal from the determination device, an SRI device for instructing insertion based on the signal, and an SRI device for instructing insertion to the core according to the instruction. A boiling water reactor core control system comprising a control rod drive mechanism for inserting the control rod.
【請求項6】請求項1において、前記制御装置は、炉心
流量調整信号を前記判定装置から受信し、前記信号に基
づきポンプシステムへポンプ回転数変化の信号を発信し
前記ポンプ回転数を変化させて前記炉心の流量を調整す
る流量調整装置を備えることを特徴とする沸騰水型原子
炉炉心制御システム。
6. The control device according to claim 1, wherein the control device receives a core flow rate adjustment signal from the determination device, and transmits a pump rotation speed change signal to a pump system based on the signal to change the pump rotation speed. And a flow rate adjusting device for adjusting the flow rate of the reactor core.
【請求項7】請求項1において、前記制御装置は、前記
判定装置で前記安定余裕度が少ないと判定された時、前
記炉心の運転を制御することを特徴とする沸騰水型原子
炉炉心制御システム。
7. The boiling water reactor core control according to claim 1, wherein the control device controls the operation of the reactor core when the judging device judges that the stability margin is small. system.
【請求項8】請求項1において、前記制御装置は、前記
原子炉の起動時に、前記炉心の運転を制御することを特
徴とする沸騰水型原子炉炉心制御システム。
8. The boiling water reactor core control system according to claim 1, wherein the control device controls the operation of the reactor core when the reactor is started.
【請求項9】請求項1において、前記制御装置は、前記
制御棒パータン交換時に、前記炉心の運転を制御するこ
とを特徴とする沸騰水型原子炉炉心制御システム。
9. The boiling water reactor core control system according to claim 1, wherein the control device controls the operation of the core when the control rod pattern is replaced.
【請求項10】請求項1ないし請求項9において、前記
安定余裕度を少なくしている複数の不安定モードを識別
する安定性モード判定器を備え、前記識別された各々の
不安定モードに対して、前記制御装置は、前記炉心の運
転をそれぞれ制御することを特徴とする沸騰水型原子炉
炉心制御システム。
10. The stability mode determiner according to claim 1, wherein the stability mode determiner discriminates a plurality of unstable modes in which the stability margin is reduced, and a stability mode determiner is provided for each of the identified unstable modes. The controller controls the operation of the reactor core, respectively, and the boiling water reactor core control system is characterized.
【請求項11】請求項1ないし請求項10において、前
記炉心監視装置は、前記遅延フィルターに基づく、τ秒
過去の時刻T1での安定余裕度の表示の他に、前記中性
子束信号を前記遅延フィルターを介さず直接前記自己相
関器、前記安定性評価器に導く系統を設け、前記系統に
基づく時刻T0での予想値も同時に表示することを特徴
とする沸騰水型原子炉炉心制御システム。
11. The core monitoring device according to claim 1, wherein the core monitoring device delays the neutron flux signal in addition to displaying a stability margin at time T1 τ seconds past based on the delay filter. A boiling water reactor core control system characterized in that a system is provided directly to the autocorrelator and the stability evaluator without passing through a filter, and an expected value at time T0 based on the system is also displayed at the same time.
【請求項12】沸騰水型原子炉炉心の監視信号であるA
PRM、あるいはLPRMの中性子束信号データ等を用
いて、前記炉心の安定余裕度を評価監視しながら前記安
定余裕度を判定し、前記判定結果に基づき前記炉心の運
転を制御する沸騰水型原子炉炉心制御方法において、 遅延フィルターに前記中性子束信号データを取り込んで
ドリフト成分を予測し、前記予測に基づき遂次的に自己
相関値を求め、前記自己相関値に基づき減幅比を算出し
て前記炉心の安定余裕度を評価することにより、前記炉
心を監視しながら前記安定余裕度を判定し、判定結果に
より、選択制御棒挿入信号を受信し、前記信号に基づき
前記制御棒挿入を指示し、前記炉心に挿入することによ
り、前記炉心の運転を制御することを特徴とする沸騰水
型原子炉炉心制御方法。
12. A supervisory signal for a boiling water reactor core
A boiling water reactor for determining the stability margin while evaluating and monitoring the stability margin of the core by using neutron flux signal data of PRM or LPRM, and controlling the operation of the core based on the determination result. In the core control method, the neutron flux signal data is taken into the delay filter to predict the drift component, the autocorrelation value is sequentially obtained based on the prediction, and the decrement ratio is calculated based on the autocorrelation value to By evaluating the stability margin of the core, to determine the stability margin while monitoring the core, depending on the determination result, receiving a selection control rod insertion signal, instructing the control rod insertion based on the signal, A boiling water reactor core control method comprising controlling the operation of the core by inserting the core into the core.
【請求項13】沸騰水型原子炉炉心の監視信号であるA
PRM、あるいはLPRMの中性子束信号データ等を用
いて、前記炉心の安定余裕度を評価監視しながら前記安
定余裕度を判定し、前記判定結果に基づき前記炉心の運
転を制御する沸騰水型原子炉炉心制御方法において、 遅延フィルターに前記中性子束信号データを取り込んで
ドリフト成分を予測し、前記予測に基づき遂次的に自己
相関値を求め、前記自己相関値に基づき減幅比を算出し
て前記炉心の安定余裕度を評価することにより、前記炉
心を監視しながら前記安定余裕度を判定し、判定結果に
より、炉心流量調整信号を受信し、前記信号に基づきポ
ンプシステムへポンプ回転数変化の信号を発信し前記ポ
ンプ回転数を変化させて前記炉心の流量を調整すること
により、前記炉心の運転を制御することを特徴とする沸
騰水型原子炉炉心制御方法。
13. A supervisory signal for a boiling water reactor core
A boiling water reactor for determining the stability margin while evaluating and monitoring the stability margin of the core using neutron flux signal data of PRM or LPRM, and controlling the operation of the core based on the determination result. In the core control method, the drift component is predicted by incorporating the neutron flux signal data into the delay filter, the autocorrelation value is obtained sequentially based on the prediction, and the width reduction ratio is calculated based on the autocorrelation value. By evaluating the stability margin of the core, the stability margin is determined while monitoring the core, the core flow rate adjustment signal is received according to the determination result, and a signal indicating a change in pump speed to the pump system based on the signal. To control the operation of the core by changing the pump rotation speed to adjust the flow rate of the core. Way.
【請求項14】請求項12または請求項13において、
前記遅延フィルターは、現在の時刻T0よりτ秒過去の
時刻T1(T1=T0ーτ)で前記ドリフト成分を予測し、
前記τ秒過去の時刻T1での前記中性子束信号より前記
ドリフト成分を差し引くことにより前記中性子束信号の
変動量を評価することを特徴とする沸騰水型原子炉炉心
制御方法。
14. The method according to claim 12 or 13,
The delay filter predicts the drift component at time T1 (T1 = T0−τ) τ seconds past the current time T0,
A boiling water reactor core control method, characterized in that the fluctuation amount of the neutron flux signal is evaluated by subtracting the drift component from the neutron flux signal at time T1 in the past τ seconds.
【請求項15】請求項12ないし請求項14において、
前記遅延フィルターでのτ(過去へのさかのぼり時間)
を前記中性子束信号変動分の基本周波数の逆数1/f以
上にしたことを特徴とする沸騰水型原子炉炉心制御方
法。
15. The method according to any one of claims 12 to 14,
Τ in the delay filter (time to go back to the past)
Is set to the reciprocal number 1 / f or more of the fundamental frequency of the fluctuation of the neutron flux signal, and a boiling water reactor core control method.
【請求項16】請求項12ないし請求項15において、
前記遅延フィルター内での前記中性子束信号のドリフト
成分予測を次式(数2)で行うことを特徴とする沸騰水
型原子炉炉心制御方法。 【数2】 ここで、Yk; 時刻Tkでの中性子束信号のドリフト
成分値 Xi; 時刻Tiでの中性子束信号データ値 ai; 係数、 bi; 係数 である。ただし、1<k<mとなり、中性子束信号デー
タの間隔をΔtとすると、(k−m)=τ/Δtとな
る。
16. The method according to claim 12 to claim 15,
A boiling water reactor core control method, wherein the drift component prediction of the neutron flux signal in the delay filter is performed by the following equation (Equation 2). [Equation 2] Here, Yk; drift component value of neutron flux signal at time Tk, Xi; neutron flux signal data value at time Ti, ai; coefficient, bi; coefficient. However, if 1 <k <m, and the interval of the neutron flux signal data is Δt, then (k−m) = τ / Δt.
JP13364095A 1995-05-31 1995-05-31 Boiling water reactor core control system and core control method Expired - Fee Related JP3226003B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP13364095A JP3226003B2 (en) 1995-05-31 1995-05-31 Boiling water reactor core control system and core control method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP13364095A JP3226003B2 (en) 1995-05-31 1995-05-31 Boiling water reactor core control system and core control method

Publications (2)

Publication Number Publication Date
JPH08327776A true JPH08327776A (en) 1996-12-13
JP3226003B2 JP3226003B2 (en) 2001-11-05

Family

ID=15109551

Family Applications (1)

Application Number Title Priority Date Filing Date
JP13364095A Expired - Fee Related JP3226003B2 (en) 1995-05-31 1995-05-31 Boiling water reactor core control system and core control method

Country Status (1)

Country Link
JP (1) JP3226003B2 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002181984A (en) * 2000-12-11 2002-06-26 Global Nuclear Fuel-Japan Co Ltd Monitoring controller for boiling water reactor
CN104778983A (en) * 2015-04-09 2015-07-15 中国核动力研究设计院 H2/H<infinity> mixed filtration-based signal delay elimination method for rhodium self-powered detector
CN104778985A (en) * 2015-04-09 2015-07-15 中国核动力研究设计院 H<infinity> filtration-based signal delay elimination method for rhodium self-powered detector

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002181984A (en) * 2000-12-11 2002-06-26 Global Nuclear Fuel-Japan Co Ltd Monitoring controller for boiling water reactor
CN104778983A (en) * 2015-04-09 2015-07-15 中国核动力研究设计院 H2/H<infinity> mixed filtration-based signal delay elimination method for rhodium self-powered detector
CN104778985A (en) * 2015-04-09 2015-07-15 中国核动力研究设计院 H<infinity> filtration-based signal delay elimination method for rhodium self-powered detector

Also Published As

Publication number Publication date
JP3226003B2 (en) 2001-11-05

Similar Documents

Publication Publication Date Title
US5406598A (en) System for monitoring power of nuclear reactor
JP4279763B2 (en) Power system stability diagnosis device, power system stabilization device, and power system contraction support device
JPH08327776A (en) Core control system of boiling water reactor, and core control method
SE513236C3 (en) System for power monitoring in a nuclear reactor
JP3256079B2 (en) Neutron flux monitoring method and device for reactor, and control rod operation control system using the same
JP4791993B2 (en) Core coolant flow rate measuring device and core coolant flow rate measuring method
JPH10221476A (en) Reactor power control device
JP3583399B2 (en) Method and apparatus for measuring void reactivity coefficient
JPH08101710A (en) Operation controller for plant
JP2005207944A (en) Reactor output control method and its system
JPH05134080A (en) Monitoring apparatus of reactor power generation plant
JPH0844415A (en) Plant stability monitor
JP4363788B2 (en) Reactor core coolant flow rate measurement method
JPH11142588A (en) Reactor automatic start device
JP2574576B2 (en) Dam gate control device
JP2003315484A (en) System for measuring and computing quantity of flow of core of boiling water reactor
JPH05134079A (en) Monitoring apparatus of reactor power distribution
RU2328780C1 (en) Method for automatic control of nuclear reaction of undercritical nuclear reactor
JPH05249271A (en) Reactor core monitor
JP2001147286A (en) Measurement method for void reactivity coefficient and device
Dudzic et al. Multivariate statistical monitoring of a continuous steel slab caster
JP2676714B2 (en) Demand control device
CN116759124A (en) Nuclear power plant reactor power control system input signal correction method
JPS6316718B2 (en)
JPS62228199A (en) Reactivity monitor device for nuclear reactor

Legal Events

Date Code Title Description
S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313115

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20070831

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080831

Year of fee payment: 7

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080831

Year of fee payment: 7

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313115

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080831

Year of fee payment: 7

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080831

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090831

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100831

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100831

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110831

Year of fee payment: 10

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120831

Year of fee payment: 11

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130831

Year of fee payment: 12

LAPS Cancellation because of no payment of annual fees