JP2585703B2 - Reactor performance prediction device - Google Patents
Reactor performance prediction deviceInfo
- Publication number
- JP2585703B2 JP2585703B2 JP63088101A JP8810188A JP2585703B2 JP 2585703 B2 JP2585703 B2 JP 2585703B2 JP 63088101 A JP63088101 A JP 63088101A JP 8810188 A JP8810188 A JP 8810188A JP 2585703 B2 JP2585703 B2 JP 2585703B2
- Authority
- JP
- Japan
- Prior art keywords
- output
- data
- core
- thermal efficiency
- thermal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
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
【発明の詳細な説明】 [発明の目的] (産業上の利用分野) 本発明は、原子炉発電プラントにおける運転管理およ
び運転計画に使用する原子炉性能予測装置に関する。DETAILED DESCRIPTION OF THE INVENTION [Object of the Invention] (Field of Industrial Application) The present invention relates to a reactor performance prediction device used for operation management and operation planning in a nuclear power plant.
(従来の技術) 一般に、原子炉発電プラントでは、原子炉の炉心シミ
ュレータを備えた原子炉性能予測装置を使用して、原子
炉の運転計画が行なわれている。(Related Art) Generally, in a nuclear power plant, a reactor operation plan is performed using a reactor performance prediction device provided with a reactor core simulator.
このような原子炉性能予測装置では、定められた運転
制限値を守って要求された炉出力を得るために、炉心に
装荷する燃料の性能に基づいて、制御棒の位置,炉心冷
却材料,および燃料の燃焼度を主な運転パラメータとし
て、原子炉の運転状態が模擬される。In such a reactor performance prediction device, in order to obtain a required reactor power while observing a specified operation limit value, the position of the control rod, the core cooling material, and the The operating state of the reactor is simulated using the burnup of fuel as a main operating parameter.
(発明が解決しようとする課題) 従来の原子炉性能予測装置は、運転日,要求された炉
心熱出力,制御棒の位置,炉心流量などを入力すること
により、炉心内の熱出力分布等が予測計算されて、その
結果が出力されていた。(Problems to be Solved by the Invention) The conventional reactor performance prediction device is designed to calculate the thermal power distribution in the reactor core by inputting the operation date, the required core thermal output, the position of the control rod, the core flow rate, etc. The prediction was calculated and the result was output.
ところで、一般に、原子力発電所を運転する際には、
電力系統における電気出力の値が運転指令として規定さ
れる。By the way, in general, when operating a nuclear power plant,
The value of the electric output in the power system is defined as an operation command.
このため、原子力発電所の運転計画を行なうために、
上記原子炉性能予測装置に、炉心熱出力を入力する際に
は、運転員は、タービン等のプラントの熱効率に基づい
て、規定された電気出力を炉出力に換算する必要があっ
た。また、上記原子炉性能予測装置からは、熱出力の値
が出力されるので、その熱出力を発電機からの電気出力
に換算して読み取る必要があった。Therefore, in order to plan the operation of the nuclear power plant,
When the core thermal output is input to the reactor performance prediction device, the operator needs to convert a specified electrical output into a reactor output based on the thermal efficiency of a plant such as a turbine. Further, since the value of the thermal output is output from the reactor performance prediction device, it is necessary to convert the thermal output into the electrical output from the generator and read it.
ところで、プラントの熱効率は、例えば、夏期には、
冬期に比べ復水器の冷却効率が低下するため、1〜2%
も熱効率が低下するというように、季節により変動する
ことが知られている。By the way, the thermal efficiency of the plant, for example, in summer,
1% to 2% because the cooling efficiency of the condenser is lower than in winter
It is known that thermal efficiency varies depending on the season, such as a decrease in thermal efficiency.
従って、原子炉性能予測装置を使用して精度の高い運
転計画を行なうためには、運転期日における上記熱効率
を正確に求めなければならない。Therefore, in order to perform a highly accurate operation plan using the reactor performance prediction device, the above-mentioned thermal efficiency on the operation date must be accurately obtained.
しかしながら、上記正確な熱効率を求めるためには、
発電プラントの長期的な運転データを採取して、それら
に基づいて判定しなければならない上、原子炉性能予測
装置を使用する際に、その判定した熱効率に基づいて、
前記のような電気出力と熱出力との換算を行なわなけれ
ばならないため、従来は、これらの作業が煩わしいとい
う問題があった。However, in order to obtain the above accurate thermal efficiency,
In addition to collecting long-term operation data of the power plant and making decisions based on them, when using the reactor performance prediction device, based on the determined thermal efficiency,
Conventionally, there is a problem that these operations are troublesome because the conversion between the electric output and the heat output has to be performed as described above.
本発明は、煩わしい作業なしに精度の高い運転計画が
行なえる原子炉性能予測装置を提供することを目的とす
る。An object of the present invention is to provide a reactor performance prediction device capable of performing an accurate operation plan without troublesome work.
[発明の構成] (課題を解決するための手段) このために本発明は、原子力発電プラントの運転パラ
メータを設定する運転パラメータ設定手段と、設定され
た前記運転パラメータに応じて原子炉の運転状態の予測
を行なう炉心シミュレータと、予測された前記運転状態
を表示する表示手段とを備えた原子炉性能予測装置にお
いて、前記原子力発電プラントの実運転時における炉心
熱出力、電気出力を示す各データを日付と共に採取する
プラントデータ収集手段と、採取した前記炉心熱出力と
電気出力とに基づいて前記原子力発電プラントの熱効率
を予測時期に応じて算出する熱効率学習手段と、前記運
転パラメータとして設定される電気出力を前記算出した
熱効率に基づいて炉心熱出力に変換する熱・電気変換手
段とを備え、その変換された炉心熱出力に応じて原子炉
の運転状態の予測を行なうことを特徴とするものであ
る。[Constitution of the Invention] (Means for Solving the Problems) For this purpose, the present invention provides an operation parameter setting means for setting an operation parameter of a nuclear power plant, and an operation state of a nuclear reactor according to the set operation parameter. In a reactor performance prediction apparatus comprising a core simulator for performing prediction of the above, and a display means for displaying the predicted operation state, the core thermal output during actual operation of the nuclear power plant, each data indicating the electrical output A plant data collection unit that collects together with the date, a thermal efficiency learning unit that calculates the thermal efficiency of the nuclear power plant according to the predicted time based on the collected core thermal output and electrical output, and electricity set as the operation parameter. Thermo-electric conversion means for converting the output to a core thermal output based on the calculated thermal efficiency. It is characterized in that to predict the operating conditions of the reactor in accordance with the core thermal power.
(作用) 原子力発電プラントの実運転時における炉心熱出力お
よび電気出力に基づいて自動的に熱効率が算出されるの
で、運転員は運転データを採取する手間をかけずに、長
期的な運転実績に基づいて正確な熱効率を求められると
共に、運転パラメータとして、電気出力を設定した場
合、自動的に電気出力から熱出力への換算が行なわれて
処理されるので、運転員は、運転データの採取や電気出
力と熱出力との換算という煩わしい作業なしに、精度の
高い運転計画が行なえるようになる。(Operation) The thermal efficiency is automatically calculated based on the core thermal output and electrical output during actual operation of the nuclear power plant, so operators can save long-term operation results without having to collect operation data. In addition to obtaining an accurate thermal efficiency based on the above, when an electric output is set as an operation parameter, the electric output is automatically converted into a heat output and the processing is performed. A highly accurate operation plan can be performed without the troublesome work of converting between the electric output and the heat output.
(実施例) 以下、本発明の実施例を、添付図面を参照しながら詳
細に説明する。Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings.
第1図は、本発明の原子炉性能予測装置の一実施例に
係る炉心管理システムのブロック構成図を示したもので
ある。この炉心管理システムは、原子力発電所サイトに
設置され、オンラインで炉心管理を行なう機能を備えて
いる。図において、プロセスデータ入力装置1は、主と
してアナログ信号入力回路よりなり、プラント2に配設
された各種センサからの検出値を入力し、各種データを
熱効率観測データとして収集するものである。電子計算
機3は、プラントデータ入力装置1により収集された上
記熱効率観測データを入力して処理する一方、炉心内の
熱出力分布等の予測計算を行なうものである。FIG. 1 shows a block diagram of a core management system according to an embodiment of a reactor performance prediction device of the present invention. This core management system is installed at a nuclear power plant site and has a function of performing core management online. In the figure, a process data input device 1 mainly comprises an analog signal input circuit, inputs detection values from various sensors provided in a plant 2, and collects various data as thermal efficiency observation data. The electronic computer 3 inputs and processes the above-mentioned thermal efficiency observation data collected by the plant data input device 1, and performs a prediction calculation of a heat output distribution and the like in the reactor core.
電子計算機3には、不揮発性メモリ上に形成された上
記熱効率観測データを記憶する熱効率データファイル
4、および上記予測計算に必要な炉心に関する各種デー
タが記憶される炉心データファイル5を備えている。こ
の炉心データファイル5内には、燃料の核定数データ,
燃焼度・核分裂生成物組成などの運転履歴データ,およ
びその時の運転パラメータ・出力分布等の瞬時データな
ど、炉心状態の初期設定と計算の初期条件設定ための各
種データが格納される。また、炉心データファイル5に
は、電子計算機3内のシミュレーション処理により更新
される運転履歴依存データと共に、計算結果一式が保存
される。The electronic computer 3 includes a thermal efficiency data file 4 formed on a nonvolatile memory for storing the thermal efficiency observation data, and a core data file 5 for storing various data relating to the core required for the prediction calculation. The core data file 5 contains nuclear constant data of fuel,
Various data for initial setting of the core state and initial conditions for calculation, such as operating history data such as burn-up and fission product composition, and instantaneous data such as operating parameters and power distribution at that time, are stored. Further, in the core data file 5, a set of calculation results is stored together with operation history dependent data updated by the simulation processing in the electronic computer 3.
また、電子計算機3を操作したり運転パラータ等を入
力する入力装置61,各種情報を表示する表示装置62や記
録紙に出力する印字装置63など、各種対話装置6が配設
されている。In addition, various interactive devices 6 such as an input device 61 for operating the computer 3 and inputting operation parameters, a display device 62 for displaying various information, and a printing device 63 for outputting on recording paper are provided.
電子計算機3内の各部は、ソフトウェアで形成される
もので、プラントデータ入力部31は、プロセスデータ入
力装置1からの熱効率観測データを入力するもので、熱
効率学習部32は、その熱効率観測データを熱効率データ
ファイル4に格納する一方、それらのデータに基づいて
熱効率計算を行なうものである。Each unit in the electronic computer 3 is formed by software, the plant data input unit 31 inputs thermal efficiency observation data from the process data input device 1, and the thermal efficiency learning unit 32 outputs the thermal efficiency observation data. While storing the data in the thermal efficiency data file 4, the thermal efficiency is calculated based on the data.
熱・電気変換部33は、算出された熱効率に基づいて、
発電機出力と炉出力との相互変換を行なうものである。The thermal-electrical conversion unit 33, based on the calculated thermal efficiency,
It performs mutual conversion between generator output and furnace output.
運転パラメータ設定部34は、入力装置61より入力され
る運転パラメータを入力するものである。炉心シミュレ
ータ35は、その入力された運転パラメータに基づいて炉
心内の熱出力分布等の予測計算を行なうものである。The operation parameter setting unit 34 inputs operation parameters input from the input device 61. The core simulator 35 performs a prediction calculation of a heat output distribution and the like in the core based on the input operation parameters.
運転データ出力部36は、得られた熱出力分布を入力さ
れた運転パラメータと共に、表示装置62あるいは印字装
置63に出力するものである。The operation data output unit 36 outputs the obtained heat output distribution to the display device 62 or the printing device 63 together with the input operation parameters.
以上の構成で、プロセスデータ入力装置1は、例え
ば、1時間という一定周期で、プラント2から、熱効率
観測データとして、電気出力,炉心熱出力およびタービ
ン復水入口海水温度の各データを収集する。なお、炉心
熱出力は、直接検出できないので、実際には、炉心に供
給される水の温度・流量,炉心から出ていく上記の温度
・流量などを検出し、炉心の熱の収支計算により求めら
れる。With the above configuration, the process data input device 1 collects, as thermal efficiency observation data, each data of the electrical output, the core thermal output, and the turbine condensing inlet seawater temperature from the plant 2 at a constant period of, for example, one hour. Since the core thermal output cannot be directly detected, the temperature and flow rate of the water supplied to the core and the above-mentioned temperature and flow rate exiting the core are actually detected and calculated by calculating the heat balance of the core. Can be
上記3種のデータは、プラントデータ入力部31を介し
て熱効率学習部32に転送される。熱効率学習部32は、上
記3種のデータに日付・時刻のデータを付加し、第2図
(a)に示すように、1組の熱効率観測データDsを作成
して、熱効率データファイル4に順次格納する。The above three types of data are transferred to the thermal efficiency learning unit 32 via the plant data input unit 31. The thermal efficiency learning unit 32 adds date / time data to the three types of data, creates a set of thermal efficiency observation data Ds as shown in FIG. Store.
同図(b)は、熱効率データファイル4内の上記各熱
効率観測データDsを記憶するための熱効率観測データテ
ーブルT1〜T12を示したもので、例えば、1月中に観測
された熱効率観測データDsは熱効率観測データテーブル
T1に、2月中に観測された熱効率観測データDsは熱効率
観測データテーブルT2というように、熱効率観測データ
テーブルT1〜T12には、それぞれ各月ごとの熱効率観測
データDsが格納される。FIG. 4B shows thermal efficiency observation data tables T 1 to T 12 for storing the respective thermal efficiency observation data Ds in the thermal efficiency data file 4. For example, the thermal efficiency observation data observed during January is shown. Data Ds is thermal efficiency observation data table
To T 1, the thermal efficiency observed data Ds observed in February and so thermal efficiency observed data table T 2, the thermal efficiency observed data table T 1 through T 12 are respectively stored thermal efficiency observation data Ds for each month You.
また、各熱効率観測データテーブルT1〜T12は、熱出
力レベル10±5%〜110±5%まで10段階のエリアに分
割されると共に、その各エリアは、さらにNo.1〜No.10
までのスロットに分割されている。Each of the thermal efficiency observation data tables T 1 to T 12 is divided into 10 areas with heat output levels of 10 ± 5% to 110 ± 5%, and each area is further divided into No. 1 to No. 10.
It is divided into slots.
熱効率データファイル4に転送された上記熱効率観測
データDsは、その炉心熱出力に対応する上記エリア内の
No.1スロットから順次格納されると共に、そのエリア内
に登録されたデータ数が同時に記憶される。また、No.1
0スロットまで格納されてそのエリアが一杯になると、
観測日時の古いNo.1スロットの熱効率観測データDsから
順に消去され、そこに新しい熱効率観測データDsが順次
格納される。The thermal efficiency observation data Ds transferred to the thermal efficiency data file 4 is stored in the area corresponding to the core thermal output.
The data is sequentially stored from the No. 1 slot, and the number of data registered in the area is simultaneously stored. No.1
When the area is filled up to 0 slots,
The thermal efficiency observation data Ds of No. 1 slot with the oldest observation date and time are deleted in order, and new thermal efficiency observation data Ds are sequentially stored therein.
ところで、海水温度は季節により変化するが、本実施
例では、海水温度が、一定温度Ts以上のとき該当する熱
効率観測データDsを夏のデータとし、海水温度が一定温
度Tw以下のとき冬のデータとして区別するために、一定
温度Ts,Twが設定されている。By the way, the seawater temperature changes depending on the season.In this embodiment, when the seawater temperature is equal to or higher than the fixed temperature Ts, the corresponding thermal efficiency observation data Ds is set to summer data, and when the seawater temperature is equal to or lower than the fixed temperature Tw, winter data is set. The constant temperatures Ts and Tw are set in order to distinguish them from each other.
熱効率学習部32は、観測されて電子計算機3より各デ
ータが入力されるごとに、その海水温度が、上記一定温
度Ts以上であった回数Ns、および上記一定温度Tw以下で
あった回数Nwを計数している。そして、各月において、
全入力回数をNiとすると、夏期係数ksおよび冬期係数kw
を、 Ks=Ns÷Ni Kw=Nw÷Ni として、算出すると共に、季節係数Knを、 Kn=Ks−Kw として算出して記憶する。The thermal efficiency learning unit 32 calculates the number Ns of times that the seawater temperature was equal to or higher than the certain temperature Ts and the number Nw of times that the seawater temperature was equal to or less than the certain temperature Tw each time the observed data was input from the computer 3. Counting. And in each month,
Assuming that the total number of inputs is Ni, the summer coefficient ks and the winter coefficient kw
Is calculated as Ks = Ns ÷ Ni Kw = Nw ÷ Ni, and the seasonal coefficient Kn is calculated and stored as Kn = Ks−Kw.
つまり、上記夏期係数ks,冬期係数kwは、その月の全
観測回数の内、海水温度が“夏温度”あるいは“冬温
度”となった割合を示すことになる。従って、季節係数
Knは、その月が、全て“夏温度”のとき、kn=1、“夏
温度”と“冬温度”とが同数のとき、kn=0、全て“冬
温度”のとき、kn=−1となる。In other words, the summer coefficient ks and the winter coefficient kw indicate the ratio of the seawater temperature at “summer temperature” or “winter temperature” in the total number of observations in the month. Therefore, the season coefficient
Kn is kn = 1 when the month is all “summer temperature”, kn = 0 when the “summer temperature” and “winter temperature” are the same number, and kn = −1 when the month is all “winter temperature”. Becomes
このようにして、最大1年間に対応する熱効率観測デ
ータDsと、各月ごとの季節係数Knとが記憶される。In this way, the thermal efficiency observation data Ds corresponding to a maximum of one year and the seasonal coefficient Kn for each month are stored.
運転員は、プラント2の運転計画のために炉心内の熱
出力分布等をシミュレーションする場合、入力装置61で
所定の起動操作を行なうと、運転パラメータ設定部34に
より表示装置62に、第3図に示すような運転パラメータ
入力画面Gが表示される。すなわち、画面の上部には、
予測したい年月日,時刻を入力するための表示部Ga、電
気出力または炉出力を入力するための表示部Gb、炉心流
量を入力するための表示部Gcが表示され、さらに、その
下部に制御棒位置を指定するための表示部Gdが表示され
る。When the operator performs a predetermined start-up operation with the input device 61 to simulate the heat output distribution or the like in the core for the operation plan of the plant 2, the operation parameter setting unit 34 causes the display device 62 to display on FIG. Is displayed as shown in FIG. That is, at the top of the screen,
A display section Ga for inputting a date and time to be predicted, a display section Gb for inputting an electrical output or a furnace output, and a display section Gc for inputting a core flow rate are displayed. A display section Gd for specifying the bar position is displayed.
運転員は、入力装置61の操作により、上記各表示部Ga
〜Gcに対して所望のデータを入力する一方、表示部Gdで
制御棒位置を指定する。この場合、表示部Gbに対応する
出力値は、第4図(a)に示すように、熱出力の絶対値
に文字“T"を付けて入力しても、同図(b)に示すよう
に、電気出力の絶対値に文字“E"を付けて入力してもよ
い。また、同図(c),(d)に示すように、その熱出
力または電気出力を、定格終力に対する比で、同様に文
字“T"または“E"付けて入力してもよい。The operator operates the input device 61 to operate the display units Ga.
While inputting desired data to .about.Gc, the control rod position is designated on the display unit Gd. In this case, as shown in FIG. 4A, the output value corresponding to the display unit Gb can be obtained by adding the letter “T” to the absolute value of the heat output as shown in FIG. Alternatively, the absolute value of the electrical output may be input with the letter “E” added. Further, as shown in FIGS. 3 (c) and 3 (d), the heat output or the electric output may be similarly input with the letter “T” or “E” in a ratio to the rated end force.
また、運転パラメータは、表示部Ga〜Gdに対応して4
種あるが、その内の1種を未知数として、“?"を入力
し、そのパラメータに対して探索計算を指示することが
できる。In addition, the operation parameters are set to 4 corresponding to the display units Ga to Gd.
Although there are species, one of them can be set as an unknown number and "?" Can be input to instruct a search calculation for the parameter.
ここでは、いま、運転員により、第4図(b)に示す
ように、電気出力の絶対値が入力されたものとする。Here, it is assumed that the operator has input the absolute value of the electric output as shown in FIG. 4 (b).
ところで、第4図(b)または(d)に示したよう
に、電気出力の指定が行なわれた場合、電気出力を熱出
力に換算する必要がある。By the way, as shown in FIG. 4 (b) or (d), when an electric output is designated, it is necessary to convert the electric output into a heat output.
前記したように、冬期は夏期に比較して、熱効率が高
い。ここで、電気出力をPe,炉心熱出力をPtとすると、
電気出力Peは、冬期には、第5図ラインW、また夏期に
は、ラインSで示すように、それぞれ炉心熱出力Ptの上
昇に対してほぼ直線的に上昇する。このため、炉心管理
上、両者の相関は一時式で近似できる。As described above, the thermal efficiency is higher in winter than in summer. Here, assuming that the electric output is Pe and the core heat output is Pt,
The electric power Pe rises almost linearly with the increase in the core thermal output Pt, as shown by the line W in FIG. 5 in winter and by the line S in summer, respectively. Therefore, in core management, the correlation between the two can be approximated by a temporary equation.
従って、係数をa〜dとおくと、電気出力Peと炉心熱
出力Ptとは、次のような相関式で表される。Accordingly, if the coefficients are a to d, the electric power Pe and the core heat output Pt are expressed by the following correlation equation.
Pe=a・Pt+b …(1) Pt=c・Pe+d …(2) また、係数a〜dは、次式のように表される。Pe = a · Pt + b (1) Pt = c · Pe + d (2) The coefficients a to d are represented by the following equations.
c=1/a…(3),d=−b/a…(4) 従って、上記係数a〜bが分かれば、電気出力Peと熱
出力Ptとの換算が行なえる。c = 1 / a ... (3), d = -b / a ... (4) Therefore, if the above coefficients a and b are known, the conversion between the electric output Pe and the thermal output Pt can be performed.
本実施例では、熱効率データファイル4に記憶してい
る熱効率観測データDs内のプラント2の実運転時データ
である炉心熱出力と電気出力とに基づいて、上記係数a
〜bを求めるようにしている。In the present embodiment, the coefficient a is calculated based on the core thermal output and the electrical output, which are the actual operation data of the plant 2 in the thermal efficiency observation data Ds stored in the thermal efficiency data file 4.
To b.
この場合、前記したようにプラントの熱効率は、海水
温度により変化するため、係数a〜bも海水温度により
変化する。このため、本実施例では、熱効率データファ
イル4内の全ての熱観測データDsに基づいて、通年用係
数a1〜b1を算出する一方、熱観測データDs内の海水温度
が、前記一定温度Ts以上である“夏温度”の熱観測デー
タDsに基づいて、夏期用係数a2〜b2を算出すると共に、
その海水温度が、前記一定温度Tw以下である“冬温度”
の熱観測データDsに基づいて、冬期用係数a3〜b3を算出
するようにしている。In this case, as described above, since the thermal efficiency of the plant changes according to the seawater temperature, the coefficients a and b also change according to the seawater temperature. Therefore, in this embodiment, on the basis of all the heat observation data Ds of the thermal efficiency data file 4, while calculating the year coefficient a 1 ~b 1, sea water temperature in the heat observation data Ds is, the certain temperature Based on the heat observation data Ds of “summer temperature” which is equal to or more than Ts, the summer coefficients a 2 to b 2 are calculated,
"Winter temperature" where the seawater temperature is equal to or lower than the constant temperature Tw.
Based on the thermal observation data Ds, and to calculate the winter coefficient a 3 ~b 3.
熱効率学習部32は、各係数a1〜b1、a2〜b2、a3〜b3を
次に示すように算出する。すなわち、第6図(a)は、
通年用係数a1〜b1の算出動作を示したもので、この場
合、熱効率データファイル4の熱効率観測データテーブ
ルT1〜T12の全熱観測データDsを読み出す(処理701)。
ついで、読み出した熱観測データDsの個数をチェックす
る(処理702)。ここで、統計処理するに十分な個数が
あれば(処理702のY)、さらに読み出した熱観測デー
タDs内の炉心熱出力の観測域の幅を判定する(処理70
3)。The thermal efficiency learning unit 32 calculates the coefficients a 1 to b 1 , a 2 to b 2 , and a 3 to b 3 as follows. That is, FIG. 6 (a)
Shows the operation of calculating the year coefficient a 1 ~b 1, in this case, reads out the total heat observed data Ds of the thermal efficiency data file 4 in the thermal efficiency observed data table T 1 through T 12 (process 701).
Next, the number of the read thermal observation data Ds is checked (process 702). Here, if there is a sufficient number for performing statistical processing (Y in step 702), the width of the observation area of the core heat output in the read heat observation data Ds is further determined (step 70).
3).
例えば、熱観測データDs内に最低80%、最高100%の
炉心熱出力データがある場合のように、その炉心熱出力
が20%以上の幅で観測されている熱観測データDsである
場合(処理703のY)、各熱観測データDs内の熱出力と
電気出力とのデータを使用して、最少自乗法により、前
記(1)式における係数a〜bを算出する(処理70
4)。For example, when the core heat output is observed at a width of 20% or more, such as when the core heat output data has a minimum of 80% and a maximum of 100% in the heat observation data Ds ( Y) of the process 703, the coefficients a and b in the equation (1) are calculated by the least squares method using the data of the thermal output and the electrical output in each thermal observation data Ds (process 70).
Four).
いま、例えば、上記データをグラフ上に×マークでプ
ロットしたとすると、第7図に示すように、各データ
は、ある程度分散したものとなり、上記算出された係数
a〜bにより、同図に示すような相関Fが求められる。Now, for example, if the above data is plotted on the graph with an X mark, as shown in FIG. 7, each data is scattered to some extent, and is shown in FIG. 7 by the calculated coefficients a and b. Such a correlation F is obtained.
ここで、求めた相関Fに対して、遠く離れた観測デー
タが生じるが、そのようなデータは、海水温度以外の要
因、例えばタービンバイパス弁開状態などが原因してい
る。そこで、このようなデータを排除するために、統計
的手法により、標準偏差σを用いて、相関Fに対して±
3σの範囲外のデータ、同図ではデータd1,d2を排除す
る(処理705)。Here, observation data that is far away from the obtained correlation F is generated. Such data is caused by factors other than the seawater temperature, for example, a turbine bypass valve open state. Therefore, in order to exclude such data, the correlation F is ±
Data out of the range of 3σ, that is, data d 1 and d 2 in the figure are excluded (process 705).
次いで、上記±3σの範囲内のデータのみを用いて、
再度、係数a〜bを算出する(処理706)。Then, using only the data within the above ± 3σ range,
The coefficients a and b are calculated again (process 706).
一方、熱効率データファイル4から読み出した熱観測
データDsの個数が、少なく統計処理するに不十分な場合
(処理702のN)、このプラント設計時において想定さ
れたいる設計値を係数a〜bとする(処理707)。On the other hand, when the number of pieces of heat observation data Ds read from the thermal efficiency data file 4 is small and insufficient for statistical processing (N in process 702), the design values assumed at the time of this plant design are represented by coefficients a and b. (Step 707).
また、炉心熱出力の観測幅がが20%未満である場合
(処理703のN)、係数aについては、上記同様に設計
値を採用する一方(処理708)、係数bは、前記最少自
乗法により算出する(処理709)。When the observation width of the core heat output is less than 20% (N in Step 703), the design value is adopted for the coefficient a in the same manner as described above (Step 708), while the coefficient b is the least squares method. (Step 709).
このように算出した係数a,bを通年用係数a1〜b1とす
る(処理701)。The coefficients a and b calculated in this manner are set as year-round coefficients a 1 to b 1 (process 701).
次に、第6図(b)に示すように、夏期用係数a2〜b2
の算出動作を実行する。すなわち、この場合、熱効率デ
ータファイル4の熱効率観測データテーブルT1〜T12内
の熱観測データDsの内、海水温度が前記一定温度Ts以上
である“夏温度”の熱観測データDsのみを、夏期データ
組として読み出す(処理801)。Next, as shown in FIG. 6 (b), in summer coefficient a 2 ~b 2
Is performed. That is, in this case, among the heat observation data Ds thermal efficiency data file thermal efficiency observed data table T 1 through T 12 4, sea water temperature is the predetermined temperature Ts or only thermal observation data Ds for "summer temperature", The data is read out as a summer data set (process 801).
次に、前記と同様に、熱観測データDsの個数をチェッ
クして(処理802)、統計処理するに十分な個数があれ
ば(処理802のY)、さらに熱観測データDs内の炉心熱
出力の観測域の幅を判定する(処理803)。Next, as described above, the number of heat observation data Ds is checked (process 802), and if there is a sufficient number for performing statistical processing (Y in process 802), the core heat output in the heat observation data Ds is further determined. The width of the observation area is determined (step 803).
ここで、20%以上の幅で観測されている場合(処理80
3のY)、各熱観測データDs内の熱出力と電気出力との
データを使用して、最少自乗法により、係数a〜bを算
出する(処理804)。Here, when the observation is made with a width of 20% or more (processing 80
3), the coefficients a and b are calculated by the least square method using the data of the heat output and the electric output in each heat observation data Ds (process 804).
ついで、算出された係数a〜bにより、熱出力と電気
出力との相関Fを求め、これに対して±3σの範囲外の
データを排除する(処理805)。そして、その±3σの
範囲内のデータのみを用いて、再度、係数a〜bを算出
する(処理806)。Next, a correlation F between the thermal output and the electrical output is obtained from the calculated coefficients a and b, and data outside the range of ± 3σ is excluded (step 805). Then, the coefficients a and b are calculated again using only the data within the range of ± 3σ (step 806).
一方、熱効率データファイル4から読み出した熱観測
データDsの個数が、少なく統計処理するに不十分な場合
(処理802のN)、前記で算出した通年用係数a1〜b
1を、ここで求める係数a,bとする(処理807)。On the other hand, when the number of pieces of heat observation data Ds read from the thermal efficiency data file 4 is small and insufficient for statistical processing (N in processing 802), the year-round coefficients a 1 to b calculated above are used.
Let 1 be the coefficients a and b found here (process 807).
また、炉心熱出力の観測幅がが20%未満である場合
(処理803のN)、通年用係数a1をここで求める係数a
とする一方(処理808)、前記最少自乗法により係数b
を算出する(処理809)。Moreover, (N of step 803) when the observation width of core thermal power is less than 20%, the coefficient determining the year for coefficients a 1 here a
(Process 808), the coefficient b is calculated by the least squares method.
Is calculated (step 809).
このように算出した係数a,bを夏期用係数a2〜b2とす
る(処理810)。Thus calculated coefficient a, b and summer coefficient a 2 ~b 2 to (operation 810).
次に、第6図(c)に示すように、冬期用係数a3〜b3
の算出動作を実行する。すなわち、この場合、熱効率デ
ータファイル4の熱効率観測データテーブルT1〜T12内
の熱観測データDsの内、海水温度が前記一定温度Tw以下
である“冬温度”の熱観測データDsのみを、冬期データ
組として読み出す(処理901)。Next, as shown in FIG. 6 (c), winter coefficient a 3 ~b 3
Is performed. That is, in this case, of the thermal observation data Ds in the thermal efficiency observation data tables T 1 to T 12 of the thermal efficiency data file 4, only the thermal observation data Ds of “winter temperature” in which the seawater temperature is equal to or lower than the constant temperature Tw, The data is read out as a winter data set (process 901).
次いで、上記処理802〜処理809と同様に、読み出した
熱観測データDs内の熱出力と電気出力とのデータに基づ
いて、係数a〜bを算出する(処理902)。そして、求
めた係数a〜bを冬期用係数a3〜b3とする(処理90
3)。Next, coefficients a and b are calculated based on the data of the thermal output and the electrical output in the read thermal observation data Ds, as in the processes 802 to 809 (process 902). Then, the obtained coefficients a and b are set as winter coefficients a 3 to b 3 (processing 90
3).
さて、前記運転員により、入力された各運転パラメー
タは、運転パラメータ設定部34より炉心シミュレータ35
に転送され、その内の月日を示すデータが炉心シミュレ
ータ35から熱・電気変換部33に転送される。The operating parameters inputted by the operator are transmitted from the operating parameter setting unit 34 to the core simulator 35.
The data indicating the date is transferred from the core simulator 35 to the thermoelectric converter 33.
熱・電気変換部33は、そのデータにより月を判別し、
その月に対応する季節係数Knを、熱効率学習部32より読
み出す。The heat / electricity conversion unit 33 determines the month based on the data,
The seasonal coefficient Kn corresponding to the month is read from the thermal efficiency learning unit 32.
次いで、その季節係数Knが、Kn≧0の場合、前記通年
用係数a1,夏期用係数a2に基づいて、係数a,bを、次式の
ように算出する。Next, when the season coefficient Kn is Kn ≧ 0, the coefficients a and b are calculated as in the following equation based on the coefficient a 1 for the whole year and the coefficient a 2 for the summer season.
a=a2×Kn+a1×(1−Kn) …(5) b=b2×Kn+b1×(1−Kn) …(6) また、上記季節係数Knが、Kn<0の場合、前記通年用
係数a1,冬期用係数a3に基づいて、係数a,bを、次式のよ
うに算出する。a = a 2 × Kn + a 1 × (1-Kn) (5) b = b 2 × Kn + b 1 × (1-Kn) (6) When the seasonal coefficient Kn is Kn <0, the whole year Based on the coefficient for use a 1 and the coefficient for winter a 3 , coefficients a and b are calculated as in the following equation.
a=a3×Kn+a1×(1−Kn) …(7) b=b3×Kn+b1×(1−Kn) …(8) 従って、例えば、その月内に観測された海水温度が、
全て「夏」温度の場合、前記したように、knは、kn=1
となるので、(5),(6)上式により、係数a,bは、
a=a2,b=b2となる。また、全て“冬温度”の場合、kn
=−1となるので、(7),(8)上式により、係数a,
bは、a=a3,b=b3となる。さらに、“夏温度”と“冬
温度”とが同数のとき、kn=0となるので、a=a1,b=
b1となる。すならち、このように、海水温度の“夏温
度”と“冬温度”との出現率より、運転月の季節が判定
され、その季節に応じた係数a,bが設定されるようにな
る。a = a 3 × Kn + a 1 × (1-Kn) (7) b = b 3 × Kn + b 1 × (1-Kn) (8) Therefore, for example, the seawater temperature observed during the month is
For all "summer" temperatures, kn is kn = 1, as described above.
(5), (6) From the above equations, the coefficients a and b are
a = a 2 and b = b 2 . In the case of all “winter temperatures”, kn
= −1, the coefficients a,
b is a = a 3 and b = b 3 . Further, when “summer temperature” and “winter temperature” are the same number, kn = 0, so that a = a 1 , b =
b becomes 1 . In other words, the season of the operating month is determined from the appearance rates of the “summer temperature” and “winter temperature” of the seawater temperature, and the coefficients a and b according to the season are set. .
いま、運転員により、運転パラメータの出力値は、第
4図(b)に示した形式で、電気出力が指定されてい
る。Now, the operator designates the output value of the operation parameter as an electrical output in the format shown in FIG. 4 (b).
ここで、熱効率学習部32は、上記(5),(6)式ま
たは(7),(8)式により算出した係数a,bから、前
記(3),(4)式により、係数c,dを求める。そし
て、その係数c,dを前記(2)式に適用して、その式に
より上記電気出力を炉心熱出力に換算する。次いで、求
めたその炉心熱出力の値を炉心シミュレータ35に転送す
る。Here, the thermal efficiency learning unit 32 calculates the coefficients c, c from the coefficients a, b calculated by the above equations (5) and (6) or the equations (7) and (8) by the above equations (3) and (4). Find d. Then, the coefficients c and d are applied to the above equation (2), and the above electric output is converted into the core thermal output according to the equation. Next, the calculated core heat output value is transferred to the core simulator 35.
炉心シミュレータ35は、その炉心熱出力と、先に運転
パラメータ設定部34から転送された各運転パラメータに
基づいた所定のシミュレーション処理により、炉心内の
熱出力分布等を予測計算する。The core simulator 35 predicts and calculates the heat output distribution and the like in the core by a predetermined simulation process based on the core heat output and each operation parameter previously transferred from the operation parameter setting unit 34.
このシミュレーション処理は、既知のもので、物理モ
デルに比較的計算時間の短い1群粗格子拡散モデルを使
用した3次元核熱水力結合計算法に基づいて行ない、オ
ンラインの運転管理に十分な応答性と精度を確保するよ
うにしている。この計算法については、例えば、文献M.
Tuiki etal;“Convergence and Acceleration of Void
Iteration in Boiling Water Reacter Core Calculatio
ns",Nucl.Sci.Eng.,64,724−732(1977)に記載されて
いる。This simulation process is based on a known three-dimensional nuclear thermal-hydraulic coupling calculation method using a first-group coarse lattice diffusion model, which has a relatively short calculation time as a physical model, and has a sufficient response to online operation management. We try to ensure the quality and accuracy. About this calculation method, for example, reference M.
Tuiki etal; “Convergence and Acceleration of Void
Iteration in Boiling Water Reacter Core Calculatio
ns ", Nucl. Sci. Eng., 64, 724-732 (1977).
このように予測計算された熱出力分布を示すデータ
は、炉心シミュレータ35から熱・電気変換部33に転送さ
れる。The data indicating the heat output distribution thus predicted and calculated is transferred from the core simulator 35 to the heat-electricity converter 33.
熱・電気変換部33は、前記上記(5),(6)式また
は(7),(8)式により算出した係数a,bを、前記
(1)式に適用して、その式により熱出力を上記電気出
力に換算する。次いで、換算した求めた電気出力データ
を炉心シミュレータ35に転送する。The thermoelectric converter 33 applies the coefficients a and b calculated by the above formulas (5) and (6) or the formulas (7) and (8) to the formula (1), and calculates the heat by the formula. The output is converted to the above electrical output. Next, the converted electrical output data is transferred to the core simulator 35.
炉心シミュレータ35は、予測計算した上記熱出力分布
等のデータ、および上記電気出力データを運転データ出
力部36に転送する。運転データ出力部36は、それらのデ
ータに基づいた所定の運転状態の予測表示、および入力
装置61で入力された運転パラメータなどを表示装置62で
表示する。The core simulator 35 transfers the predicted and calculated data such as the heat output distribution and the electric output data to the operation data output unit 36. The operation data output unit 36 displays on the display device 62 a prediction display of a predetermined operation state based on the data, an operation parameter input by the input device 61, and the like.
これにより、表示装置62に、熱出力分布等が表示され
るが、このとき、熱出力の値と共に、それらが電気出力
に換算された値も同時に表示される。As a result, the heat output distribution and the like are displayed on the display device 62. At this time, the value of the heat output and the value converted into the electric output are also displayed at the same time.
一方、以上の例では、運転員により、第4図(b)に
示したように、電気出力が指定された場合を説明した
が、同図(d)に示したように、電気出力が定格出力に
対する比で指定された場合には、運転パラメータ設定部
34において、電気出力が絶対値に変換されて同様に処理
される。On the other hand, in the above example, the case where the electric output is designated by the operator as shown in FIG. 4B has been described. However, as shown in FIG. If specified as a ratio to output, the operation parameter setting section
At 34, the electrical output is converted to an absolute value and similarly processed.
また、第4図(a)に示したように、熱出力で指定さ
れた場合には、熱・電気変換部33による上記換算処理は
実行されないで、そのまま炉心シミュレータ35により、
所定のシミュレーション処理が実行される。また、この
場合、予測計算された上記熱出力分布等のデータは、熱
・電気変換部33で電気出力データに換算されて、上記と
同様に表示される。さらに、第4図(c)に示したよう
に、熱出力が、定格出力に対する比で指定された場合、
運転パラメータ設定部34で熱出力が絶対値に変換されて
同様に処理される。Further, as shown in FIG. 4 (a), when the heat output is designated, the conversion process by the thermoelectric converter 33 is not executed, and the core simulator 35
A predetermined simulation process is executed. In addition, in this case, the data such as the heat output distribution that is predicted and calculated is converted into electric output data by the thermoelectric converter 33 and displayed in the same manner as described above. Further, as shown in FIG. 4 (c), when the heat output is specified by a ratio to the rated output,
The operation parameter setting unit 34 converts the heat output into an absolute value and performs the same processing.
また、一方、第4図(e)に示したように、未知数と
して“?"と入力された場合、あるいは、第3図の他の運
転パラメータが未知数として“?"で入力された場合、原
子炉運転時における臨界条件より既知の方法で算出し、
その算出値に基づいて、運転状態の予測処理処理を実行
する。On the other hand, as shown in FIG. 4 (e), when "?" Is input as an unknown number, or when other operation parameters in FIG. 3 are input as "?" Calculated by a known method from the critical conditions during furnace operation,
Based on the calculated value, an operation state prediction process is executed.
例えば、運転パラメータとして、電気出力が入力さ
れ、日付・時刻が未知数として指定された場合、炉心シ
ミュレータ35は、指定された電気出力になる日時を探索
して、所定の運転状態の予測処理を実行する。この場
合、月の経過に応じて熱効率が変化するので、各月ごと
に、前記相関係数a,bを算出して、電気出力を炉心熱出
力に換算しながら、上記日時を探索する。また、炉心熱
出力が未知数として指定された場合、炉心シミュレータ
35により算出された炉心熱出力が電気出力に換算され
て、表示装置62で表示される。For example, when an electric output is input as an operation parameter and the date / time is specified as an unknown number, the core simulator 35 searches for a date and time at which the specified electric output is obtained, and executes a predetermined operation state prediction process. I do. In this case, since the thermal efficiency changes with the passage of the month, the above date and time are searched for while calculating the correlation coefficients a and b for each month and converting the electric output to the core heat output. If the core thermal output is specified as unknown, the core simulator
The core thermal output calculated by 35 is converted into an electrical output and displayed on the display device 62.
以上のように、この実施例によれば、プラントデータ
入力部31により、プラント2から炉心熱出力,電気出
力,海水温度のデータを採取して日付と共に熱効率観測
データDsとして、過去1年間に相当するデータを熱効率
ファイル4で記憶している。また、このとき、採取した
海水温度の“夏温度”と“冬温度”の出現の度合いを、
各月毎に季節係数Knとして記憶している。そして、熱効
率学習部32により、採取した炉心熱出力と電気出力とに
基づいて、プラント2の実運転時の熱効率を算出すると
共に、原子炉の運転状態予測時期の熱効率を、過去の同
月における季節係数Knに基づいて、上記算出した熱効率
を補正することにより求めるようにしている。As described above, according to this embodiment, the plant data input unit 31 collects core thermal output, electrical output, and seawater temperature data from the plant 2 and, together with the date, obtains thermal efficiency observation data Ds corresponding to the past year. Is stored in the thermal efficiency file 4. At this time, the degree of appearance of “summer temperature” and “winter temperature” of the collected seawater temperature was
It is stored as a seasonal coefficient Kn for each month. Then, the thermal efficiency learning unit 32 calculates the thermal efficiency during the actual operation of the plant 2 based on the collected core thermal output and the electrical output, and calculates the thermal efficiency at the predicted operating state of the reactor in the past month in the same month. The calculated thermal efficiency is obtained by correcting the calculated thermal efficiency based on the coefficient Kn.
これにより、運転員は運転データを採取する手間をか
けずに、長期的な運転実績に基づいた正確な熱効率が求
められるようになる。As a result, the operator is required to have accurate thermal efficiency based on long-term operation results without taking the trouble of collecting operation data.
また、運転パラメータとして電気出力が入力された場
合に、その電気出力を予測した上記熱効率に基づいて炉
心熱出力に変換して、その炉心熱出力に応じて原子炉の
運転状態の予測を行なうと共に、その予測結果である炉
心熱出力を、上記判定した熱効率に基づいて電気出力に
変換して、表示させるようにしたので、運転員は、運転
データの採取や電気出力と熱出力との換算という煩わし
い作業なしに、精度の高い運転計画が行なえるようにな
るようになる。Further, when an electric output is input as an operation parameter, the electric output is converted into a core heat output based on the predicted thermal efficiency, and an operation state of the reactor is predicted according to the core heat output. The core thermal output, which is the result of the prediction, is converted into an electrical output based on the determined thermal efficiency, and is displayed.Therefore, the operator is required to collect operating data and to convert the electrical output and the thermal output. An operation plan with high accuracy can be performed without troublesome work.
この場合、プラントの運転期間が長くなり、熱効率観
測データテーブルT1〜T12内の熱効率観測データDsの個
数が多くなるほど、正確な係数a,bが算出され、精度の
良い運転計画が行なえる。In this case, the operation period of the plant is increased, as the number of thermal efficiency observed data Ds of the thermal efficiency observed data table T 1 through T 12 increases, accurate coefficients a, b are calculated, enables good operation plan accuracy is .
次に、第8図は、炉心管理システムの他の実施例を示
したものである。図において、第1図と同一符号は同一
ブロックを示し、第1図と異なる点は、熱・電気変換部
33は、炉心シミュレータ35と接続されず、運転パラメー
タ設定部34および運転データ出力部36と接続されている
点である。Next, FIG. 8 shows another embodiment of the core management system. In the figure, the same reference numerals as those in FIG. 1 denote the same blocks, and the difference from FIG.
Reference numeral 33 denotes a point not connected to the core simulator 35 but connected to the operation parameter setting unit 34 and the operation data output unit 36.
この場合、入力装置61で、運転パラメータとして電気
出力が入力されると、その電気出力のデータが運転パラ
メータ設定部34から熱・電気変換部33に転送されて熱出
力に変換され、運転パラメータ設定部34に戻される。そ
して、その熱出力と他の運転パラメータが炉心シミュレ
ータ35に転送され、所定の運転状態の予測処理が実行さ
れる。また、その予測結果である炉心内の熱出力分布の
データは、運転データ出力部36に転送された後、さらに
熱・電気変換部33に転送され、ここで電気出力に変換さ
れる。この電気出力のデータが、運転データ出力部36を
介して、表示装置62で表示される。In this case, when an electric output is input as an operation parameter by the input device 61, the data of the electric output is transferred from the operation parameter setting unit 34 to the thermoelectric conversion unit 33 and converted into a heat output, and the operation parameter setting is performed. Returned to section 34. Then, the heat output and other operation parameters are transferred to the core simulator 35, and a predetermined operation state prediction process is executed. Further, the data of the heat output distribution in the core, which is the result of the prediction, is transferred to the operation data output unit 36, and further transferred to the heat / electricity conversion unit 33, where it is converted to electric power. The electrical output data is displayed on the display device 62 via the operation data output unit 36.
これによれば、炉心シミュレータ35は、第1図の場合
のように、熱・電気変換部33とのデータの入出力処理が
ないため、従来機能で済むようになる。なお、この場
合、運転パラメータの電気出力を入力して、日付を未知
数として指定した場合、炉心シミュレータ35は、各月ご
とに電気出力を熱出力に変換する処理が必要なので、日
付の探索が実行できなくなる。According to this, the core simulator 35 does not need to perform data input / output processing with the thermoelectric converter 33 as in the case of FIG. In this case, when the electric output of the operation parameter is input and the date is specified as an unknown number, the core simulator 35 needs to perform a process of converting the electric output to the heat output every month, so that the date search is executed. become unable.
なお、以上の実施例では、プラント2から実運転時の
プラントデータを観測するようにしたが、予めプラント
のプロセスデータ,設計値データをデータファイルに格
納しておき、それらのデータに基づいて熱効率を判定す
るというように、オフライン処理を行なうようにしても
よい。In the above embodiment, the plant data during actual operation is observed from the plant 2. However, the process data and design value data of the plant are stored in advance in a data file, and the thermal efficiency is determined based on the data. In such a case, the offline processing may be performed.
また、季節係数Knは、各月毎に算出しておき、運転予
測する日付をm月d日とした場合、m月に対応する季節
係数Knに基づいて、前記(5)〜(8)式により前記
(1)式で使用する係数a,bを算出するようにしたが、
季節係数Knは、各月をさらに上旬,中旬,下旬、週毎、
あるいは日毎に分け、それぞれの季節係数Knを算出して
おき、運転予測日であるd日に対応する季節係数Knに基
づいて、係数a,bをを算出するようにしてもよい。Further, the seasonal coefficient Kn is calculated for each month, and when the date on which the operation is predicted is m month d day, the above-mentioned equations (5) to (8) are calculated based on the seasonal coefficient Kn corresponding to m month. Calculates the coefficients a and b used in the above equation (1).
The seasonal coefficient Kn is calculated for each month in the early, middle, late, weekly,
Alternatively, the seasonal coefficients Kn may be calculated for each day, and the coefficients a and b may be calculated based on the seasonal coefficient Kn corresponding to the d day which is the predicted operation day.
また、季節に応じた係数a,bを算出するために、熱効
率観測データテーブルT1〜T12内の熱効率観測データDs
を海水温度により分類したが、単に観測した月日に対応
する熱効率観測データDsより、係数a,bを算出するよう
にしてもよい。The coefficient corresponding to the season a, in order to calculate the b, thermal efficiency observed data Ds of the thermal efficiency observed data table T 1 through T 12
Are classified according to the seawater temperature, but the coefficients a and b may be calculated simply from the thermal efficiency observation data Ds corresponding to the observed date.
さらに、電気出力と炉心熱出力との換算を行なう場
合、前記(1),(2)式のように1次式にしたが、2
次式にして精度を向上させたり、炉圧や主蒸気温度もパ
ラメータとする他の相関式を適用したりするなど、各種
の換算方法が考えられる。Further, when converting the electric power and the core heat output, a linear equation was used as in the above equations (1) and (2).
Various conversion methods are conceivable, such as using the following formula to improve the accuracy or applying another correlation formula that also uses the furnace pressure and the main steam temperature as parameters.
[発明の効果] 以上のように、本発明によれば、運転員は運転データ
を採取する手間をかけずに、長期的な運転実績に基づい
て正確な熱効率を求められると共に、運転パラメータと
して設定される電気出力を、算出して得た上記熱効率に
基づいて炉心熱出力に変換するようにしたので、運転員
は、運転データの採取や電気出力と熱出力との換算とい
う煩わしい作業なしに、精度の高い運転計画が行なえる
ようになる。[Effects of the Invention] As described above, according to the present invention, an operator can obtain accurate thermal efficiency based on long-term operation results without taking the trouble of collecting operation data, and can set the operation parameters as operation parameters. Since the electrical output to be converted to the core thermal output based on the thermal efficiency obtained by calculation, the operator, without the cumbersome work of collecting operating data and converting the electrical output and thermal output, A highly accurate operation plan can be performed.
第1図は本発明の本発明の原子炉性能予測装置の一実施
例に係る炉心管理システムのブロック構成図、第2図
(a)は熱効率観測データの説明図、同図(b)は熱効
率観測データテーブルの説明図、第3図は運転パラメー
タ入力時の画面図、第4図(a)〜(e)は出力値の入
力方法の説明図、第5図は炉心熱出力と発電機の電気出
力との相関を示すグラフ図、第6図(a)は通年用係数
の算出動作を示す動作フローチャート、同図(b)は夏
期用係数の算出動作を示す動作フローチャート、同図
(c)は冬期用係数の算出動作を示す動作フローチャー
ト、第7図は最小自乗法による相関係数の算出法を示す
グラフ図、第8図は炉心管理システムの他の実施例を示
すブロック構成図である。 1……プロセスデータ入力装置、2……プラント、3…
…電子計算機、4……熱効率データファイル、5……炉
心データファイル、31……プラントデータ入力部、32…
…熱効率学習部、33……熱・電気変換部、34……運転パ
ラメータ設定部、35……炉心シミュレータ、36……運転
データ出力部、61……入力装置、62……表示装置、63…
…印字装置。FIG. 1 is a block diagram of a core management system according to one embodiment of the reactor performance prediction apparatus of the present invention, FIG. 2 (a) is an explanatory diagram of thermal efficiency observation data, and FIG. FIG. 3 is an explanatory view of an observation data table, FIG. 3 is a screen view when operating parameters are input, FIGS. 4 (a) to (e) are explanatory views of an input method of output values, and FIG. FIG. 6 (a) is an operation flowchart showing a whole-year coefficient calculation operation, FIG. 6 (b) is an operation flowchart showing a summer-time coefficient calculation operation, FIG. 6 (c) FIG. 7 is an operation flowchart showing a coefficient calculation operation for winter, FIG. 7 is a graph showing a method of calculating a correlation coefficient by the least squares method, and FIG. 8 is a block diagram showing another embodiment of the core management system. . 1 ... process data input device, 2 ... plant, 3 ...
... Electronic computer, 4 ... Thermal efficiency data file, 5 ... Core data file, 31 ... Plant data input section, 32 ...
… Thermal efficiency learning unit, 33… thermal-electric conversion unit, 34… operating parameter setting unit, 35… core simulator, 36… operating data output unit, 61… input device, 62… display device, 63…
... Printing device.
Claims (2)
定する運転パラメータ設定手段と、設定された前記運転
パラメータに応じて原子炉の運転状態の予測を行なう炉
心シミュレータと、予測された前記運転状態を表示する
表示手段とを備えた原子炉性能予測装置において、前記
原子力発電プラントの実運転時における炉心熱出力、電
気出力を示す各データを日付と共に採取するプラントデ
ータ収集手段と、採取した前記炉心熱出力と電気出力と
に基づいて前記原子力発電プラントの熱効率を予測時期
に応じて算出する熱効率学習手段と、前記運転パラメー
タとして設定される電気出力を前記算出した熱効率に基
づいて炉心熱出力に変換する熱・電気変換手段とを備
え、その変換された炉心熱出力に応じて原子炉の運転状
態の予測を行なうことを特徴とする原子炉性能予測装
置。An operating parameter setting means for setting an operating parameter of a nuclear power plant, a core simulator for predicting an operating state of a nuclear reactor in accordance with the set operating parameter, and displaying the predicted operating state. A reactor performance prediction device comprising: a reactor core thermal output during actual operation of the nuclear power plant; and plant data collecting means for collecting data indicating electric power along with a date; and the collected core thermal output. Thermal efficiency learning means for calculating the thermal efficiency of the nuclear power plant according to the predicted time based on the power output and the electrical output; and heat converting the electrical output set as the operating parameter into a core thermal output based on the calculated thermal efficiency.・ Electric conversion means are provided, and the operation state of the reactor is predicted according to the converted core heat output. Reactor performance prediction apparatus according to claim.
熱・電気変換手段により、前記炉心シミュレータの予測
結果である炉心熱出力を、前記算出された熱効率に基づ
いて電気出力に変換すると共に、その電気出力を表示す
ることを特徴とする原子炉性能予測装置。2. The method according to claim 1, wherein the thermal / electrical conversion means converts a core thermal output, which is a prediction result of the core simulator, into an electrical output based on the calculated thermal efficiency. A reactor performance prediction device for displaying the electric output of the reactor.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP63088101A JP2585703B2 (en) | 1988-04-12 | 1988-04-12 | Reactor performance prediction device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP63088101A JP2585703B2 (en) | 1988-04-12 | 1988-04-12 | Reactor performance prediction device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH01260399A JPH01260399A (en) | 1989-10-17 |
| JP2585703B2 true JP2585703B2 (en) | 1997-02-26 |
Family
ID=13933478
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP63088101A Expired - Lifetime JP2585703B2 (en) | 1988-04-12 | 1988-04-12 | Reactor performance prediction device |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2585703B2 (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6748348B1 (en) * | 1999-12-30 | 2004-06-08 | General Electric Company | Design method for nuclear reactor fuel management |
| KR100415487B1 (en) * | 2001-04-30 | 2004-01-24 | 한국수력원자력 주식회사 | Nuclear heating method for integral reactor |
| JP5122071B2 (en) * | 2005-09-09 | 2013-01-16 | 中国電力株式会社 | Power supply plan creation system and power supply plan creation method |
| JP4986938B2 (en) * | 2008-06-17 | 2012-07-25 | 中国電力株式会社 | Apparatus, method and program for calculating total weight of fissile elements in nuclear reactor |
| JP4937234B2 (en) * | 2008-11-20 | 2012-05-23 | 中国電力株式会社 | Seawater temperature rise prediction system, seawater temperature rise prediction method |
| JP5710282B2 (en) | 2011-01-12 | 2015-04-30 | 株式会社東芝 | Core monitoring system |
-
1988
- 1988-04-12 JP JP63088101A patent/JP2585703B2/en not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JPH01260399A (en) | 1989-10-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Pellegrinetti et al. | Nonlinear control oriented boiler modeling-a benchmark problem for controller design | |
| DE69520097T2 (en) | SYSTEM FOR REAL-TIME OPTIMIZATION AND PRESENTATION OF THE PROFIT | |
| CN107871181A (en) | A kind of cogeneration units generation schedule formulating method, apparatus and system | |
| CN116579149B (en) | Reliability analysis method of electric-thermal comprehensive energy system based on chaos polynomial | |
| CN105957568A (en) | System and method for simulating nuclear heat releasing by electric heating | |
| JP2585703B2 (en) | Reactor performance prediction device | |
| CN120450612A (en) | A digital intelligent management method and system for water conservancy project construction | |
| JP3215279B2 (en) | Method and apparatus for operating energy equipment | |
| US7366273B2 (en) | Method of determining margins to operating limits for nuclear reactor operation | |
| Sumner et al. | BOP-301 and BOP-302R: Test Definitions and Analyses | |
| JP3999467B2 (en) | Combustion vibration prediction device | |
| Lin et al. | The sub-metered HVAC implemented for demand response dataset | |
| KR101877627B1 (en) | Apparatus and method for testing power performance of generator | |
| CN101609480A (en) | Inter-node phase relation identification method of electric system based on wide area measurement noise signal | |
| Rinker et al. | Using the reactive current profile of a feeder to determine optimal capacitor placement | |
| CN111913444B (en) | A chemical process monitoring method based on time series multi-block modeling strategy | |
| Aguirre | Computer-aided analysis and design of control systems using model approximation techniques | |
| CN104361020A (en) | Curve overlapping method based on configuration tool | |
| Sharp | PWR steady-state capability of WRAP-A water reactor analysis package | |
| JPS6252604A (en) | Simulation device | |
| JPS6151208A (en) | Prediction simulator | |
| JP2004126832A (en) | Dam inflow prediction device and dam inflow prediction program | |
| Tseng et al. | Validation of the dynamic simulation of a natural circulation loop | |
| CN119106524A (en) | A method, system and device for real-time combustion fault identification of coal-fired boilers | |
| Adorján et al. | VERONA-plus: extended core-monitoring system for WWER-440 type nuclear power plants |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20071205 Year of fee payment: 11 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081205 Year of fee payment: 12 |
|
| EXPY | Cancellation because of completion of term | ||
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081205 Year of fee payment: 12 |