JPH0484799A - Water quality control method for plant - Google Patents

Water quality control method for plant

Info

Publication number
JPH0484799A
JPH0484799A JP2197964A JP19796490A JPH0484799A JP H0484799 A JPH0484799 A JP H0484799A JP 2197964 A JP2197964 A JP 2197964A JP 19796490 A JP19796490 A JP 19796490A JP H0484799 A JPH0484799 A JP H0484799A
Authority
JP
Japan
Prior art keywords
water quality
input
layer
output
pattern
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.)
Pending
Application number
JP2197964A
Other languages
Japanese (ja)
Inventor
Masanori Sakai
政則 酒井
Takuya Takahashi
卓也 高橋
Katsumi Mabuchi
勝美 馬渕
Hiroshi Yamauchi
博史 山内
Noriyuki Onaka
大中 紀之
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 Ltd
Original Assignee
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 Ltd filed Critical Hitachi Ltd
Priority to JP2197964A priority Critical patent/JPH0484799A/en
Publication of JPH0484799A publication Critical patent/JPH0484799A/en
Pending 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 contrive water quality control by allowing the neural circuit model of a hierarchic structure to learn the typical pattern of an input parameter value of various time points as a teaching signal to input a water quality pattern which is not learned by the neural model so as to operate. CONSTITUTION:An input value is input into an input layer and this signal is output into the neurone element model of an intermediate layer 2. The multiplication and addition of these values and weighting factor are calculated by the use of specified equations on the neurone element model of the intermediate layer 2, and an output value for an output layer 3 is determined according to this compensation. Also the output value of the intermediate layer successively calculates the multiplication and addition of a next intermediate layer and a next output layer to obtain the calculated value of the output layer 3. In order to execute the learning of a neural network, after the output layer 3, a comparator layer 4 and a teacher signal layer 5 are provided further, the signal of the output layer 3 and a teaching signal are input into the comparator layer 4 and the output signal is compared with the teaching signal. The size of the weighting factor is modified so as to lessen the errors. With the use of the modified value, recalculation and the teaching signal are compared with each other to successively repeat modification so as to lessen the errors.

Description

【発明の詳細な説明】 〔産業上の利用分野〕 本発明はプラントの各部に設けられたセンサシステムに
より得られた情報をニューニルネットワ−’7を用いて
処理することにより最適な水質制御をする方法に関する
[Detailed Description of the Invention] [Industrial Application Field] The present invention provides optimal water quality control by processing information obtained by sensor systems installed in each part of a plant using New Nil Network '7. Regarding how to.

〔従来の技術〕[Conventional technology]

従来の分析技術へのニューニルネットワークの適用例が
ぶんせき4.(1988年)第280頁から第285頁
の第284頁において論じられている。
An example of the application of New Nil Network to conventional analysis techniques is shown in 4. (1988), pp. 280-285, p. 284.

〔発明が解決しようとする課題〕[Problem to be solved by the invention]

上記従来技術はニューラルネットワークそのものを用い
て定量分析に必要な重み関数そのものを求めてゆくとい
うものであった。本発明はこのニューラルネットワーク
をプラントの水質診断制御に適用した。即ち、センサシ
ステムにより得られた水質情報をニューラルネットに入
力することによりそれらの入力情報に対する最適水質制
御手段を提供し、プラントの水質制御をすることを目的
とする。
The above-mentioned conventional technology uses a neural network itself to obtain the weighting function itself necessary for quantitative analysis. The present invention applies this neural network to plant water quality diagnostic control. That is, the purpose is to provide optimal water quality control means for the input information by inputting water quality information obtained by a sensor system into a neural network, thereby controlling the water quality of a plant.

〔課題を解決するための手段〕[Means to solve the problem]

上記目的を達成するため、本発明はセンサからの水質情
報の入力層から出力層までの間に少なくとも一層の中間
層を設ける階層構造の神経回路モデルを用い、過去の水
質履歴情報のうち、異なる時点の複数の入力変数値の代
表的なパターンを入力信号とすると共に代表的なパター
ンに対応する制御変数値を教師信号として、神経回路モ
デルに学習させる。この学習した神経回路モデルに入力
変数値として未学習の水質パターンを入力することによ
り目的とする水質制御のための変数値を求め、制御シス
テムを運転するようにした。
In order to achieve the above object, the present invention uses a neural circuit model with a hierarchical structure in which at least one intermediate layer is provided between the input layer of water quality information from the sensor and the output layer. A neural circuit model is caused to learn by using a representative pattern of a plurality of input variable values at a time as an input signal and a control variable value corresponding to the representative pattern as a teacher signal. By inputting an unlearned water quality pattern as an input variable value to this learned neural circuit model, variable values for the target water quality control were obtained and the control system was operated.

〔作用〕[Effect]

ニューラルネットはセンサシステムにより得られた過去
の水質情報、あるいは実験室データによりどのような制
御信号を発生させたら良いか、数数の入力信号パターン
に対応させ制御システム運転用の出力信号を教師信号と
して予め学習させることができる。従って未学習の水質
情報パターンを入力することにより、その入カバターン
に対応する最適の出力信号を学習レベルに応じて出力し
てくれるため、充分に学習されたニューラルネットを用
いれば人が対応不能と思えるレベルの水質パターン入力
信号に対しても最適な水質制御運転可能な出力制御信号
を得ることができるので信頼性の高いプラントの水質診
断制御ができる。
The neural network uses past water quality information obtained by the sensor system or laboratory data to determine what kind of control signal should be generated, and generates output signals for control system operation in response to several input signal patterns. It can be learned in advance as Therefore, by inputting an unlearned water quality information pattern, the optimal output signal corresponding to that input pattern will be output according to the learning level. Since it is possible to obtain an output control signal that enables optimal water quality control operation even for a water quality pattern input signal of a seemingly low level, highly reliable water quality diagnostic control of the plant can be performed.

〔実施例〕〔Example〕

〈実施例1〉 本実施例では原子カプラント水質に関し、溶存酸素(D
o)、溶存水質(DH)、溶存過酸化水素(DHzOz
)構造材腐食電位(ECP)のモニタパターンをニュー
ニルネットワークを用いて学習させ、そのデータパター
ンが得られる原因を推定し、水質制御用のシステムを可
動させる例について示す6第1図の原子炉周辺のブロッ
ク図中、4の計装管内にDo、DH,DHzOz、EC
P各種の電気化学センサー1がセットされている。電気
化学センサの詳細を第3図に示す。この電気化学センサ
からの情報は5のポテンショスタット、6のデータ処理
、電気化学計測コントローラを通し、8のニューニルネ
ットに送り込まれ、ニューラルネット処理された情報パ
ターンに応じ9のガス、及び、薬品注入系が指令を受け
、水質をコントロール管理する仕組みになっている。得
られた水質情報のパターンを与えている水質の状態、即
ち、平常状態か、パターンの変動原因が海水リークか等
の情報は7の遠隔制御指令装置及び分析結果表示、出力
装置に、随時1表示される。10は原子炉給水配管、1
1は原子炉再循環系、12は原子炉圧力容器、13はタ
ービン、2はドライヤ、3は炉心である。
<Example 1> In this example, dissolved oxygen (D
o), dissolved water quality (DH), dissolved hydrogen peroxide (DHzOz
) The reactor shown in Figure 1 shows an example of learning the monitoring pattern of structural material corrosion potential (ECP) using a new network, estimating the cause of the data pattern, and activating the water quality control system. In the surrounding block diagram, Do, DH, DHzOz, EC are in the instrumentation pipe 4.
P Various electrochemical sensors 1 are set. Details of the electrochemical sensor are shown in Figure 3. Information from this electrochemical sensor is sent through the potentiostat (5), data processing (6), and electrochemical measurement controller (8) to the neural net (8), and depending on the information pattern processed by the neural network, gas and chemicals (9) are sent. The injection system receives commands to control and manage water quality. Information such as the state of the water quality that gives the pattern of the obtained water quality information, i.e., whether it is in a normal state or whether the cause of the fluctuation in the pattern is a seawater leak, is sent to the remote control command device and the analysis result display and output device in 7 at any time. Is displayed. 10 is the reactor water supply pipe, 1
1 is a reactor recirculation system, 12 is a reactor pressure vessel, 13 is a turbine, 2 is a dryer, and 3 is a reactor core.

第2図にニューラルネットの階層構造を示しである。第
4図は第1図の電気化学センサ1により得られた水質情
報の内、D O、D H、D H2O2の三種の化学種
濃度が表示システム7に示された。
Figure 2 shows the hierarchical structure of the neural network. In FIG. 4, among the water quality information obtained by the electrochemical sensor 1 of FIG. 1, the concentrations of three chemical species, D O, D H, and D H2O2, are shown on the display system 7.

予め学習された入圧カバターンから第2図に示すニュー
ラルネットを用いて処理した結果DOの上昇、DHの上
昇原因(測定回か月日)として原子炉の出力上昇が原因
と考えられ対応処置として9のガス注入系から水素ガス
を注入したところDOは、再び、測定表示開始時レベル
まで戻すことに成功し、ニューラルネットを通したセン
サ情報パターン出力表示は平常水質の表示になった。
As a result of processing using the neural network shown in Figure 2 from the pre-learned input pressure cover pattern, it was thought that the cause of the increase in DO and DH (measurement times) was an increase in reactor output, and countermeasures were taken. When hydrogen gas was injected from the gas injection system No. 9, the DO was successfully returned to the level at the start of the measurement display, and the sensor information pattern output display through the neural network became a display of normal water quality.

第3図の電気化学センサ拡大図で14が直径50μmの
白金ミクロワイヤ作用極、15が参照電極、又は、EC
P計、16がpHセンサ、17か白金対極、18は炉水
通過可能レベルの液絡用孔、19はタングステンメタラ
イジングアルミナ、20はセンサケースで5US316
L鋼製、21はMIケーブル(Mineral In5
ulated Cable) 22はMIケーブル集積
管、23は溶接部である。
In the enlarged view of the electrochemical sensor in Figure 3, 14 is a platinum microwire working electrode with a diameter of 50 μm, and 15 is a reference electrode, or EC.
P meter, 16 is pH sensor, 17 is platinum counter electrode, 18 is hole for liquid junction at a level that allows reactor water to pass through, 19 is tungsten metallized alumina, 20 is sensor case, 5US316
Made of L steel, 21 is MI cable (Mineral In5
22 is an MI cable integration pipe, and 23 is a welded part.

次に第4図に示される電気化学センサ出力結果を得る分
析手法について詳細に示す。この分析手法にもニューラ
ルネットを適用しである。
Next, the analysis method for obtaining the electrochemical sensor output results shown in FIG. 4 will be described in detail. Neural networks are also applied to this analysis method.

第3図に示した電気化学センサを用い、このセンサを第
1図に示す計装管4に設置し、ランダムパルスポルタン
メトリ及び第2図に示したニューラルネットワークを適
用することにより溶存酸素。
Using the electrochemical sensor shown in FIG. 3, installing this sensor in the instrumentation tube 4 shown in FIG. 1, and applying random pulse portammetry and the neural network shown in FIG. 2, dissolved oxygen can be detected.

過酸化水素、水素濃度(C0(021v C0(H2O
21C01H2))を同時定量分析したものである。
Hydrogen peroxide, hydrogen concentration (C0 (021v C0 (H2O
21C01H2)) was simultaneously quantitatively analyzed.

ランダムパルスポルタンメトリはバタン認識的手法を導
入したポルタンメトリである。複雑な外乱信号(ランダ
ムなパルス電位信号)を用いて被定量分析系を多様に揺
さぶることにより、それぞれの被測定系の電気化学反応
の情報を含む多量の電解電流データ(バタン情報)を得
ることを試みた。即ち、上記化学種の様に反応の電位が
近接し、分離定量が困難な場合にも、それぞれの被測定
物質の電気化学反応のプロセスの内、電子移動反応の速
度、パラメータ等の違いに基づき、下記に示すデータ圧
縮処理により得られたデータマトリックス中の特徴ある
バタン(特徴ベクトル)に着目し、特徴ベクトルの大き
さと検量線から定量分析、データマトリックス中に現れ
る特徴ベクトルの出現バタンの違いにより定性分析を行
う。
Random pulse portammetry is a portammetry that introduces a batan recognition method. By shaking the quantitative analysis system in various ways using a complex disturbance signal (random pulse potential signal), a large amount of electrolytic current data (bang information) including information on the electrochemical reaction of each system to be measured can be obtained. I tried. In other words, even when the reaction potentials of the above chemical species are close to each other and separation and quantification are difficult, it is possible to analyze the electrochemical reaction process of each analyte based on the differences in the rate of electron transfer reaction, parameters, etc. , we focused on the characteristic patterns (feature vectors) in the data matrix obtained by the data compression process shown below, and conducted quantitative analysis based on the size of the feature vectors and the calibration curve. Perform qualitative analysis.

ランダムパルスポルタンメトリは、■バタン情報(電解
電流データ)の取得、■バタン情報の圧縮、■圧縮バタ
ン情報からの目的情報(C’ +oz+ 。
Random pulse portammetry consists of: 1. Obtaining baton information (electrolytic current data), 2. Compressing baton information, and 2. Purpose information (C' + oz+ ) from the compressed baton information.

COu+zoz+ + C’ u+zt値)の抽出、と
いう手順で実行される。■、■、■の各ステップの概要
は次の通りである。すなわち、■では、まず、測定系に
入力されるランダムパルス電位に対する応答信号として
酸素、過酸化水素、水素の電解電流を測定する。本実施
例の場合、ランダムパルスの電位信号は二つの電位波形
パラメータ(m、n)を用いて記述する。即ち、初めの
ランダムパルスのパルス電位に対応する電位をパラメー
タmで、二番目のランダムパルスのパルス電位に対応す
る電位をパラメータnで表わす。従って5mをX軸にと
り、nをy軸にとり、各(xt y)座標上の入力電位
に対応する応答信号である電解電流をそのパラメータ値
に応じてマトリックスの形に測定データを整理する。■
では、■で得られた64 (8x8)ポイントの入力電
位信号の応答信号である電流成分のデータマトリックス
を、数個の成分を持つデータベクトルに変換、圧縮する
。このデータ圧縮方法には、フーリエ変換や、画像処理
の分野で用いられているアダマール変換、カルーネンー
レーブ変換、バール変換等がある0本実施例ではこの中
で一部元アダマール変換を用いた。■では、■で得られ
たデータベクトルの中の特徴ベクトルを抽出し、C0+
oz+* C0+uzoz+、 C0utz+  を決
定する。すなわち、得られたデータペルトルの各成分値
を入力とし、C’ (021、C’(8202) + 
C’+o2+  の各値が出力となる関数を探してセッ
トする。本実施例では、この関数をニューラルネットワ
ークを用いて決定した。以下、本実施例における■、■
This is executed in the following steps: extraction of COu+zoz++C'u+zt value). The outline of each step of ■, ■, and ■ is as follows. That is, in (2), first, the electrolytic currents of oxygen, hydrogen peroxide, and hydrogen are measured as response signals to the random pulse potential input to the measurement system. In the case of this embodiment, the random pulse potential signal is described using two potential waveform parameters (m, n). That is, the potential corresponding to the pulse potential of the first random pulse is represented by a parameter m, and the potential corresponding to the pulse potential of the second random pulse is represented by a parameter n. Therefore, 5m is set on the X axis, n is set on the y axis, and the electrolytic current, which is a response signal corresponding to the input potential on each (xt y) coordinate, is arranged in the form of a matrix according to its parameter value. ■
Now, the data matrix of the current component, which is the response signal of the input potential signal of 64 (8x8) points obtained in step (2), is converted and compressed into a data vector having several components. This data compression method includes Fourier transform, Hadamard transform, Karhunen-Löve transform, Barr transform, etc. which are used in the field of image processing.In this example, among these, the partial Hadamard transform was used. . In ■, extract the feature vector from the data vector obtained in ■, and
oz+*C0+uzoz+, C0utz+ are determined. That is, input each component value of the obtained data pertle, and calculate C' (021, C' (8202) +
Find and set a function whose output is each value of C'+o2+. In this example, this function was determined using a neural network. Below, ■, ■ in this example
.

■の各ステップの詳細を述べる。③ Describe the details of each step.

まず、■の「バタン情報の取得」の詳細を述べる。First, details of ``obtaining slam information'' will be described.

第15図はランダムパルス電位信号の印加方法の説明図
である。ランダムパルス電位信号といっても、全く任意
に印加したのではC0+oz++C0tHxox+ 、
 C0(821に関する情報を得ることができない場合
もありえる0本実施例では、ノーマルパルスとリバース
パルスを併用したポルタンメトリがC0tOZl+ C
0u+zoz+、 C’+Hz+  を分離決定できる
ことに着目し、このパルス印加方法を基本型としたラン
ダムパルス電位信号を変えた。第5図の上部は酸素、過
酸化水素、水素の三種混合系における電流(i)−電位
(E)曲線の概略図である。電位がポジティブな領域で
は過酸化水素と水素の酸化電流が、ネガティブな領域で
は酸素と過酸化水素の還元電流がそれぞれIi!測され
る。酸化電流及び還元電流が観測される電位領域を、そ
れぞれE^ ECとし、2’  (Q=1.2.3・・
・)個のレベルに分割する。21という数を採用するの
は、コンピュータによる処理を実行しやすくするためで
ある。本実施例ではそれぞれ八個CQ=3)のレベルに
分割し、各レベルをE^(n)、EC(n)(n=1.
2.3・・8)と表した。また、酸素、過酸化水素、水
素の電極反応が起こらない電位をEc 、過酸化水素と
水素の酸化反応の拡散限界電流が得られる電位をElと
する。
FIG. 15 is an explanatory diagram of a method of applying a random pulse potential signal. Although it is called a random pulse potential signal, if it is applied completely arbitrarily, C0+oz++C0tHxox+,
In this example, portammetry using both a normal pulse and a reverse pulse detects C0tOZl+C.
Noting that Ou+zoz+ and C'+Hz+ can be determined separately, we changed the random pulse potential signal using this pulse application method as the basic type. The upper part of FIG. 5 is a schematic diagram of a current (i)-potential (E) curve in a ternary mixed system of oxygen, hydrogen peroxide, and hydrogen. In the positive potential region, the oxidation currents of hydrogen peroxide and hydrogen are Ii!, and in the negative potential region, the reduction currents of oxygen and hydrogen peroxide are Ii! be measured. The potential regions where oxidation current and reduction current are observed are respectively E^EC, and 2' (Q=1.2.3...
・Divide into ) levels. The reason why the number 21 is adopted is to facilitate processing by a computer. In this embodiment, each level is divided into eight levels (CQ=3), and each level is divided into eight levels: E^(n), EC(n) (n=1.
2.3...8). Furthermore, the potential at which the electrode reaction of oxygen, hydrogen peroxide, and hydrogen does not occur is Ec, and the potential at which the diffusion limit current for the oxidation reaction of hydrogen peroxide and hydrogen is obtained is El.

第5図の下部のB1〜B3は電気化学センサの作用電極
に印加する三種類のランダムパルス電圧信号を表わす。
B1-B3 at the bottom of FIG. 5 represent three types of random pulse voltage signals applied to the working electrode of the electrochemical sensor.

第一種のランダムパルスは、Ec−+EC(m)→EC
(n)の順に印加される。本実施例では、Ecを1秒、
Eo(m)、Eo(n)をそれぞれ0.05秒とし、E
C(n)の終了時に電流値をサンプリングした。m、n
はそれぞれECの分割されたレベル(1〜8)を表わす
。Eo(m)とEC(n)の電位の組み合わせパターン
は8X8=64通りあり、8×8のデータマトリックス
XCNPにmとnの組に対応した電流値を格納する。こ
のようにして得られたXCNPは被検水中の酸素、過酸
化水素の還元反応の電流データから或るので、CO+o
z+。
The first type of random pulse is Ec-+EC(m)→EC
(n) are applied in this order. In this example, Ec is 1 second,
Eo(m) and Eo(n) are each 0.05 seconds, and E
The current value was sampled at the end of C(n). m, n
each represents a divided level (1 to 8) of EC. There are 8×8=64 combination patterns of the potentials of Eo(m) and EC(n), and current values corresponding to the pairs of m and n are stored in the 8×8 data matrix XCNP. Since the XCNP obtained in this way is based on the current data of the reduction reaction of oxygen and hydrogen peroxide in the test water, CO+o
z+.

COn+zo2+ に関する情報とほぼ同質の情報を持
つ。
It has almost the same information as the information regarding CON+zo2+.

第二種のランダムパルスEc−)E^(m)→E^(n
)の順に印加される。本実施例では、ECを1秒E^(
m)、E^(n)をそれぞれ0.05 秒とし、E^(
n)の終了時に電流値をサンプリングした。
Second kind of random pulse Ec-)E^(m)→E^(n
) are applied in this order. In this example, EC is set to 1 second E^(
m) and E^(n) are each 0.05 seconds, and E^(
Current values were sampled at the end of step n).

この場合も第一種の場合と同様、mとnの組に対応した
8×8のデータマトリックスXANPが得ら九るe X
ANPは被検水中の過酸化水素、水素の酸化反応の電流
データから或るので、C0(ozr。
In this case, as in the case of the first type, an 8×8 data matrix XANP corresponding to the pair m and n is obtained.
Since ANP is determined from the current data of the oxidation reaction of hydrogen peroxide and hydrogen in the test water, C0 (ozr.

c0utz+に関する情報を持つ。Contains information regarding c0utz+.

第三種のランダムパルスは、Ec→E1→EC(m)→
E’(n)の順に印加される0本実施例では、Ec、 
EC(m) + EC(n )  を第一種のランダム
パルスの場合とそれぞれ同時間、Elを0.1秒とし、
EC(n)の終了時に電流値をサンプリングした。この
場合も第一種、第二種の場合と同様、mとnの組に対応
した8×8のデータマトリックスX CRPが得られる
。第三種のランダムパルスが第一種と異なる点は、EC
(m)、EC(n)の前段に順にElが印加されること
である。すなわち、X CRPには、被検水中の酸素、
過酸化水素の還元反応の情報に加えて、Elで過酸化水
素の酸化により生成した酸素の還元反応の情報が含まれ
ることになる。
The third type of random pulse is Ec→E1→EC(m)→
In this example, Ec,
EC(m) + EC(n) is the same time as in the case of the first type random pulse, El is 0.1 second,
Current values were sampled at the end of EC(n). In this case, as in the case of the first type and the second type, an 8×8 data matrix X CRP corresponding to the set of m and n is obtained. The third type of random pulse differs from the first type in that EC
(m) and EC(n) are sequentially applied with El. In other words, X CRP includes oxygen in the test water,
In addition to information on the reduction reaction of hydrogen peroxide, information on the reduction reaction of oxygen generated by oxidation of hydrogen peroxide with El is included.

次に、■の「バタン情報の圧縮」の詳細を述べる。Next, the details of ``compression of button information'' will be described.

データマトリックスベクトルに変形、圧縮する前に、マ
トリックスXCRP NP を(1)式で定義する。
Before transforming and compressing into a data matrix vector, the matrix XCRP NP is defined by equation (1).

Xo(RP−NPI =X0Rp−X”pip    
    −(1)この操作により、″被検水中の酸素、
過酸化水素の還元反応の情報″がほとんど相殺されるの
で、XC(RP−11PIには、゛過酸化水素の酸化に
より生成した酸素の還元反応の情報”が濃縮される。
Xo(RP-NPI =X0Rp-X"pip
-(1) Through this operation, ``oxygen in the test water,''
Since the information about the reduction reaction of hydrogen peroxide is almost canceled out, the information about the reduction reaction of oxygen produced by the oxidation of hydrogen peroxide is concentrated in XC(RP-11PI).

データマトリックスXCRP NP t XCNP、 
XANPは、(2)式によってそれぞれマトリックスY
’(RP−NPI、 Y’NPI YANPに変換され
る。
Data Matrix XCRP NP t XCNP,
XANP is determined by matrix Y according to equation (2), respectively.
'(RP-NPI, Y'NPI Converted to YANP.

Y=H−X−H・・・(2) ただし、 ・・・(3) である。マトリックスYの各成分を とおくと、Xの持つ情報はYの64個の成分の内ある特
定の成分に現れる。すなわち、Xの64個の成分に分散
していた情報がYの特定の数個の成分に圧縮されたこと
になる。本実施例の場合はy11+ y11a+ 3’
fit+ y17+ y71の各成分に情報が圧縮され
た。ここでYo(RP−NPI 、 YONP、 YA
NPのy11+ y151 ’/ 61+ y17+ 
y7を成分からそれぞれZC(RP−NPI 、 ZC
NPI ZANPなるベクトルを作る。
Y=H−X−H...(2) However,...(3). For each component of the matrix Y, information possessed by X appears in a specific component among the 64 components of Y. In other words, the information that was dispersed over 64 components of X has been compressed into several specific components of Y. In this example, y11+ y11a+ 3'
Information was compressed into each component of fit+y17+y71. Here Yo(RP-NPI, YONP, YA
NP's y11+ y151'/61+ y17+
ZC (RP-NPI, ZC
Create a vector NPI ZANP.

Z”  (y11+  yi13+  y51t  y
17t  3’71)  ”’(5)以上の説明により
、ベクトルZ C+pp−NPI 。
Z” (y11+ yi13+ y51t y
17t 3'71) ''(5) According to the above explanation, vector Z C+pp-NPI.

2 ONP I Z ANPには、それぞれC’(Hz
oz+ + C0+oz+十C0(H2O21、C’o
+zo2+ +C01H21に関する情報が高密度に圧
縮されていることになる。
2 ONP I Z ANP has C' (Hz
oz+ + C0+oz+ten C0 (H2O21, C'o
This means that the information regarding +zo2+ +C01H21 is compressed with high density.

次に、■の「圧縮バタン情報からの目的情報の抽出」の
詳細を述べる。
Next, details of ``extraction of target information from compressed button information'' will be described.

任意のC01H21すC0tHzoz+、 C0(82
1に被検水に対して得られたベクトルZC(RP−NP
Iに対し、C’ n+zoz+ = f (Z 0ut
p−NPI )      −(6)を満たすスカラー
場fを準備する。スカラー場fは、予め、C0(H2O
2) が既知でC0toz+ + C0(H2)が任意
の何種類かの被検水に対して求めておく。
Any C01H21C0tHzz+, C0(82
1, the vector ZC (RP-NP
For I, C' n+zoz+ = f (Z 0ut
p-NPI ) - Prepare a scalar field f that satisfies (6). The scalar field f is defined in advance as C0(H2O
2) is known and C0toz+ + C0(H2) is determined for several arbitrary types of water to be tested.

本実施例では、ベクトルZC+RP−NPIの各成分を
入力とじCO+ozoz ) を出力とするニューラル
ネットワークを用いて、このスカラー場fを決定した。
In this example, this scalar field f was determined using a neural network that inputs each component of the vector ZC+RP-NPI and outputs CO+ozoz ).

C0(Hzoz+ 未知の被検水に対しても(6)式の
演算によってC0(Hzoz+ が求められる・同様に
任意のC’ +oz+ 、 C0tHzoz+ 、 C
01H21の被検水中で得られたZ ANPに対して、
C0+oz+ = g (Z’NP) −C0(Hzo
z+   −(7)C0(Hz+ = h (ZANP
) −C’ +82021  − (8)を満たすスカ
ラー場g、hを準備する。スカラー場g+hは予め、C
0ntzoz+とC0(021が既知、またはC’ +
Hzoz+とCO+oz+が既知の被検水中に対してそ
れぞれ求めておく。本実施例では、ベクトルZ ONp
の各成分を入力としC0toz+ + C0<ozoz
+を出力するニューラルネットワーク、および、ベクト
ルZANPの各成分を入力とじC0Toz++C0(H
202+  を出力とするニューラルネットワークをそ
れぞれ用いて、このスカラー場gvhを決定した。C0
(ov 、 COH4zoz+ * C’+nz+未知
の被検水に対しても、(6)式と(7)式、(8)式に
よりC0+oz+ e C0+oz+z+ 、 C0+
oz+ がそれぞれ求められる。
C0(Hz+) Even for unknown test water, C0(Hz+) can be found by calculating equation (6).Similarly, any C'+oz+, C0tHz+, C
For Z ANP obtained in the test water of 01H21,
C0+oz+ = g (Z'NP) -C0(Hzo
z+ −(7)C0(Hz+ = h (ZANP
) -C' +82021 - Prepare scalar fields g and h that satisfy (8). The scalar field g+h is previously written as C
0ntzoz+ and C0 (021 is known, or C' +
Hz+ and CO+oz+ are determined in advance for the known test water. In this example, vector Z ONp
Input each component of C0toz+ + C0<ozoz
A neural network that outputs
This scalar field gvh was determined using each neural network with 202+ outputs. C0
(ov, COH4zoz+ * C'+nz+ Also for unknown test water, C0+oz+ e C0+oz+z+ , C0+ by equations (6), (7), and (8).
oz+ is determined respectively.

次にニューラルネットワークでの基本的な計算方法を示
す。例として前述のランダムパルスに適用した例で示す
Next, we will show the basic calculation method using neural networks. As an example, an example in which the above-mentioned random pulse is applied will be shown.

まず、設定された信号値Zl−Znの各々に重み係数W
IJを乗じ、さらにこれらを加算する演算(積和演算)
で計算される。
First, a weighting coefficient W is assigned to each of the set signal values Zl−Zn.
An operation that multiplies IJ and then adds these (product-sum operation)
is calculated.

CJ(2)= ΣW J + (2←1)・zt(i)
   ・・(1)n=1 ここで、Zi(1):入力層(第−層)の値、Wo(2
←l)二人力層(第−層)の1番目の変数から中間層(
第二層)のj番目のニューロン素子モデルへの重み係数
、CJ(2):中間層(第二層)のj番目のニューロン
素子モデルへの入力総和値である。
CJ (2) = ΣW J + (2←1)・zt(i)
...(1) n=1 Here, Zi(1): value of input layer (-th layer), Wo(2
←l) From the first variable of the two-manpower layer (layer -) to the middle layer (
Weighting coefficient for the j-th neuron element model of the second layer (second layer), CJ(2): is the total input value to the j-th neuron element model of the intermediate layer (second layer).

入力層のニューロン素子モデルでは、CJ(2)の代償
に応じてここでの出力値が次式で計算される。
In the neuron element model of the input layer, the output value here is calculated according to the compensation of CJ(2) using the following equation.

ZJ(2)= 1/(1e−”j(2))    ・・
12)(2)式の計算内容は第20図のような関係であ
る。すなわち、同図のグラフに示すように、CJ(2)
の値に応じて0”からパ1”までの間の値がZJ(2)
として得られる。計算値ZJ(2)は、さらに出力層へ
送られ、出力層でも同様の計算が実行される。
ZJ(2)=1/(1e-"j(2))...
12) The calculation contents of equation (2) are as shown in FIG. 20. That is, as shown in the graph of the same figure, CJ(2)
The value between 0'' and Pa1'' is ZJ(2) depending on the value of
obtained as. The calculated value ZJ(2) is further sent to the output layer, where similar calculations are performed.

次に、ニューラルネットワークでの計算方法の概要につ
いて説明する。
Next, an overview of the calculation method using a neural network will be explained.

前述した入力値Zl(1)は第20図の入力層に入力さ
れ、この信号(値)は中間層のニューロン素子モデルに
出力される。中間層のニューロン素子モデルでは、これ
ら出力値ztD)と重み係数w+、(2←1)との積和
Cj(2)を(1)式で計算し、この代償に応じて出力
層への出力値Zl(2)を(2)式で決定する。同様に
して、中間層の出力値z、(2)は、更に、中間層(第
二層)と出力層(第三層)との重み係数w lJ (s
←2)との積和CJ(3)を次式で計算する。
The input value Zl(1) described above is input to the input layer of FIG. 20, and this signal (value) is output to the neuron element model of the intermediate layer. In the neuron element model of the intermediate layer, the sum of products Cj (2) of these output values ztD) and the weighting coefficient w+, (2←1) is calculated using equation (1), and the output to the output layer is adjusted according to this compensation. The value Zl(2) is determined using equation (2). Similarly, the output value z, (2) of the intermediate layer is further calculated by the weighting coefficient w lJ (s
The sum of products CJ(3) with ←2) is calculated using the following formula.

ここで、Zl(2):中間層(第二層)の値、WJI(
3←2):中間層(第二層)のi番目の変数から出力層
(第三層)のj番目のニューロン素子モデルへの重み係
数、C,(3):出力層(第三層)のj番目のニューロ
ン素子モデルへの入力総和値である。
Here, Zl(2): value of the middle layer (second layer), WJI(
3←2): Weighting coefficient, C, from the i-th variable of the intermediate layer (second layer) to the j-th neuron element model of the output layer (third layer), (3): Output layer (third layer) is the total input value to the j-th neuron element model.

さらに、CJ(3)の大小に応じて出力層へ出力値ZJ
(3)を次式で計算する。
Furthermore, the output value ZJ is sent to the output layer depending on the size of CJ(3).
(3) is calculated using the following formula.

ZJ(3)= 1 / (l  e−cj (3)) 
   −(4)このようにして、出力層の計算値ZJ(
3)が得られる。
ZJ(3)=1/(le-cj(3))
-(4) In this way, the calculated value ZJ(
3) is obtained.

ニューラルネットでの学習を実行するには、出力層の後
にさらに第20図に示すように、比較層と教師信号層と
を設け、出力層の信号と教師信号とが比較層に入力され
、ここで出力信号と教師信号とが比較される。この誤差
が小さくなるように、重み係数Wハ(3←2)およびW
 J I (2←1)の大きさを修正する。
In order to perform learning in a neural network, a comparison layer and a teacher signal layer are further provided after the output layer as shown in FIG. 20, and the output layer signal and teacher signal are input to the comparison layer, and the The output signal and the teacher signal are compared. In order to reduce this error, the weighting coefficient W(3←2) and W
Correct the size of J I (2←1).

この修正値を用いて、再度(1)ないしく4)式の計算
並びに教師信号との比較を行なうと、同様に誤差がでて
くる。この誤差が小さくなるように、再度重み係数W 
J i(3←2)及び”tLi(2←1)の大きさを修
正する。
If the calculation of equations (1) to 4) and comparison with the teacher signal are performed again using this corrected value, errors will similarly occur. In order to reduce this error, the weighting coefficient W
The magnitudes of J i (3←2) and “tLi (2←1) are corrected.

このようにして、重み係数WJ iを繰り返して修正し
ていき、誤差が十分ちいさくなるまで続ける。
In this way, the weighting coefficient WJ i is repeatedly corrected until the error becomes sufficiently small.

最初は重み係数はランダム(乱数で発生)与えるので、
誤差は大きいが、出力信号は次第に教師信号値に近づい
ていく。ここで記号c、zはそれぞれ濃度、特徴ベクト
ルに対応する。
Initially, the weighting coefficients are given randomly (generated by random numbers), so
Although the error is large, the output signal gradually approaches the teacher signal value. Here, symbols c and z correspond to density and feature vector, respectively.

このようにして誤差を修正していく方法は、誤差逆伝搬
法と呼ばれ、ルーメルハート(Ru+++elhart
)らによって考案された技術を利用する。
The method of correcting errors in this way is called the error backpropagation method.
) and others.

〔発明の効果〕〔Effect of the invention〕

本発明によれば、得られた水質情報パターンに対応した
水質制御方法において、過去の水質情報パターンに対応
する最適な出力信号パターンをニューラルネットに学習
させることができ、複雑な水質情報パターンにも既学習
レベルに応じ最適な水質制御に必要な出力信号パターン
を得ることができるので、過去の水質情報を基に、現在
の水質情報のパターンに対応した最適な水質制御方法を
得るという効果がある。また、ニューラルネットの学習
レベルを、随時、上げることにより、水質制御をさらに
高精度化できる。
According to the present invention, in a water quality control method that corresponds to the obtained water quality information pattern, it is possible to make a neural network learn the optimal output signal pattern corresponding to the past water quality information pattern, and even with complex water quality information patterns. Since it is possible to obtain the output signal pattern necessary for optimal water quality control according to the learning level, it has the effect of obtaining the optimal water quality control method corresponding to the current water quality information pattern based on past water quality information. . Furthermore, by increasing the learning level of the neural network as needed, water quality control can be made even more precise.

【図面の簡単な説明】[Brief explanation of the drawing]

第1図は本発明の一実施例の原子カプラントに採用した
プラントの水質診断制御システムを示す系統図、第2図
はニューラルネットの階層構造を示すブロック図、第3
図は電気化学センサの断面図、第4図は溶存化学種濃度
とニューラルネットを用いた水質制御結線を示す特性図
、第5図は酸素、過酸化水素及び水素の三種混合系にお
ける電流(i)−電位(E)曲線図および電気化学セン
サの作用電極に印加するランダムパルス信号を表わすグ
ラフである。 1・・・電気化学センサ、2・・・ドライヤ、3・・・
炉心、4・・・計装管、5・・・ポテンショスタット、
6・・・電気化学計測コントローラ、7・・・遠隔制御
指令装置及び分析結果表示・出力装置、8・・・ニュー
ラルネット、9・・・ガス及び薬品注入系、1o・・・
原子炉給水配管、11・・・原子炉再li環系、12・
・・原子炉圧力容器、13・・・タービン、14・・・
作用電極、15・・・参照電極又はECP計、16・・
・pH計、17・・・対極、18・・・液絡用孔、19
・・・タングステンメタライズアルミナ、 20・・・センサケース。 21・・・MI ケーブル、 22・・・M■ケーブル集積管、 23・・・溶 接部。 宅 瞥 図 8−−−二1−う1しキ、7ト 范3図 高50 時間(助
Fig. 1 is a system diagram showing a plant water quality diagnosis control system adopted in an atomic couplant according to an embodiment of the present invention, Fig. 2 is a block diagram showing the hierarchical structure of a neural network, and Fig. 3
Figure 4 is a cross-sectional view of the electrochemical sensor, Figure 4 is a characteristic diagram showing the concentration of dissolved chemical species and water quality control connections using a neural network, and Figure 5 is the current (i )-potential (E) curve diagram and a graph representing a random pulse signal applied to a working electrode of an electrochemical sensor. 1... Electrochemical sensor, 2... Dryer, 3...
Reactor core, 4... Instrumentation tube, 5... Potentiostat,
6... Electrochemical measurement controller, 7... Remote control command device and analysis result display/output device, 8... Neural network, 9... Gas and chemical injection system, 1o...
Reactor water supply piping, 11...Reactor reli ring system, 12.
...Reactor pressure vessel, 13...Turbine, 14...
Working electrode, 15... Reference electrode or ECP meter, 16...
・pH meter, 17...Counter electrode, 18...Liquid junction hole, 19
...Tungsten metallized alumina, 20...Sensor case. 21... MI cable, 22... M■ cable accumulation pipe, 23... welded part. Takubetsu Map 8 --- 21-U1 Shiki, 7 Toto 3 Diagrams High 50 Hours (Support)

Claims (1)

【特許請求の範囲】 1、プラントの監視対象部に設置された水質センサによ
つて得られる各種の水質情報を基にプラントの水質制御
方法において、 各種の前記水質センサの情報を入力する入力層をもち、
前記入力層から出力層までの間に少なくとも一層の中間
層をもつ階層構造の神経回路モデルを用い、過去のプラ
ント水質情報のうち、異なる時点の複数の入力変数値の
代表的なパターンを入力信号とし、前記代表的なパター
ンに対応する水質制御用の変数値を教師信号として前記
神経回路モデル学習させ、学習した前記神経回路モデル
に前記入力変数値として未学習のパターンを入力するこ
とにより目的の前記水質制御用の変数値を求めることを
特徴とするプラント水質診断制御システム。 2、請求項1において、時間とともに変化する複数の前
記プラントの前記水質情報に応じて、少なくとも一個の
水質制御対象が目標状態となるような制御対象の水質制
御変数値を求め、プラントの水質を制御するもので、各
種水質センサの情報を入力する入力層をもち、前記入力
層から出力層までの間に少なくとも一層の中間層をもつ
た階層構造の前記神経回路モデルを用い、過去の前記プ
ラントの前記水質情報のうち、異なる時点の複数の入力
変数値の代表的なパターンを入力信号とすると共に前記
代表的なパターンに対応する水質制御用の変数値を前記
教師信号として前記神経回路モデルに学習させ、学習し
た前記神経回路モデルに前記入力変数値として未学習の
パターンを入力して目的の前記水質制御変数値を求める
プラント水質診断制御方法。 3、請求項1または2において、前記神経回路モデルの
学習時に、或る時点の前記入力変数値パターンと或る時
点の一定時間前の時点の前記入力変数値パターンとを同
時に入力信号として用いるとともに前記或る時点の制御
変数値を前記教師信号として用い、複数の入力信号につ
いて学習した後、前記神経回路モデルに現時点の前記入
力変数値パターンおよび現時点より一定時間前の時点の
前記入力変数値パターンを同時に入力することにより現
時点の制御変数値を求めるプラント水質診断制御方法。 4、請求項1、2または3において、時間とともに変化
する複数の前記水質情報の前記入力変数値に応じて制御
対象が目標状態となるような制御対象の制御変数値を求
めるプラント水質診断制御方法。 5、請求項1、2、3または4において、前記神経回路
モデルの学習した結果に基づいて、前記各入力変数の各
制御変数との間の因果関係を抽出し、前記因果関係に基
づいて前記プロセスの運転を支援するプラント水質診断
制御方法。
[Scope of Claims] 1. In a water quality control method for a plant based on various types of water quality information obtained by water quality sensors installed in monitoring target parts of a plant, an input layer for inputting information from various water quality sensors. have,
Using a neural circuit model with a hierarchical structure having at least one intermediate layer between the input layer and the output layer, representative patterns of multiple input variable values at different points in time among past plant water quality information are input signals. Then, the neural circuit model is trained using variable values for water quality control corresponding to the representative pattern as a teacher signal, and the unlearned pattern is input as the input variable value to the learned neural circuit model to achieve the objective. A plant water quality diagnosis control system, characterized in that the variable value for water quality control is determined. 2. In claim 1, water quality control variable values of the controlled objects are determined so that at least one water quality controlled object reaches a target state according to the water quality information of the plurality of plants that changes over time, and the water quality of the plant is controlled. The neural circuit model has an input layer that inputs information from various water quality sensors, and has a hierarchical neural circuit model with at least one intermediate layer between the input layer and the output layer. Among the water quality information, a representative pattern of a plurality of input variable values at different times is used as an input signal, and a variable value for water quality control corresponding to the representative pattern is used as the teacher signal to the neural circuit model. A plant water quality diagnosis and control method that involves learning and inputting an unlearned pattern as the input variable value to the learned neural circuit model to obtain the target water quality control variable value. 3. According to claim 1 or 2, when learning the neural circuit model, the input variable value pattern at a certain point in time and the input variable value pattern at a certain point in time before the certain point in time are simultaneously used as input signals, and After learning about a plurality of input signals using the control variable value at a certain time as the teacher signal, the neural circuit model is given the input variable value pattern at the current time and the input variable value pattern at a certain time before the current time. A plant water quality diagnosis control method that calculates the current control variable values by inputting the following simultaneously. 4. A plant water quality diagnosis control method according to claim 1, 2 or 3, for determining a control variable value of a controlled object such that the controlled object reaches a target state in accordance with the input variable values of the plurality of water quality information that change with time. . 5. In claim 1, 2, 3 or 4, a causal relationship between each of the input variables and each control variable is extracted based on the learned result of the neural circuit model, and based on the causal relationship, the A plant water quality diagnosis control method that supports process operation.
JP2197964A 1990-07-27 1990-07-27 Water quality control method for plant Pending JPH0484799A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2197964A JPH0484799A (en) 1990-07-27 1990-07-27 Water quality control method for plant

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2197964A JPH0484799A (en) 1990-07-27 1990-07-27 Water quality control method for plant

Publications (1)

Publication Number Publication Date
JPH0484799A true JPH0484799A (en) 1992-03-18

Family

ID=16383252

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2197964A Pending JPH0484799A (en) 1990-07-27 1990-07-27 Water quality control method for plant

Country Status (1)

Country Link
JP (1) JPH0484799A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110618245A (en) * 2019-09-30 2019-12-27 河南沃海水务有限公司 Water quality data preprocessing circuit
JP2021081385A (en) * 2019-11-22 2021-05-27 日立Geニュークリア・エナジー株式会社 Device for predicting furnace water quality of nuclear reactor

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110618245A (en) * 2019-09-30 2019-12-27 河南沃海水务有限公司 Water quality data preprocessing circuit
JP2021081385A (en) * 2019-11-22 2021-05-27 日立Geニュークリア・エナジー株式会社 Device for predicting furnace water quality of nuclear reactor
US11942228B2 (en) 2019-11-22 2024-03-26 Hitachi-Ge Nuclear Energy, Ltd. Device for prediction of reactor water quality of nuclear reactor

Similar Documents

Publication Publication Date Title
De Carvalho et al. Simultaneous determination of phenol isomers in binary mixtures by differential pulse voltammetry using carbon fibre electrode and neural network with pruning as a multivariate calibration tool
Lindner et al. A systematic workflow for oscillation diagnosis using transfer entropy
CN108562811A (en) Complicated electrical energy power quality disturbance analysis method based on the memory of two-way shot and long term
CN111433689A (en) Generation of control system for target system
JP4190517B2 (en) Process model parameter adjustment support apparatus and method
JPH0484799A (en) Water quality control method for plant
CN117710802A (en) Gravity field direction suitability analysis method based on image texture features
Cotrim et al. Neural network meta-models for FPSO motion prediction from environmental data with different platform loads
CN115936079B (en) Calibration concentration selection method suitable for gas sensor array
JP3740382B2 (en) Process model parameter adjusting apparatus and method
Walter Attribution and detection of anthropogenic climate change using a backpropagation neural network
Pérez-Pérez et al. ANFIS and Takagi–Sugeno interval observers for fault diagnosis in bioprocess system
Sumalatha et al. Novel method of system identification
CN115167129B (en) Robust fault detection method, computing device, and storage medium
CN116008379B (en) Automatic titration system, method and device based on model fitting and machine learning
Lastres et al. Use of neural networks in solving interferences caused by formation of intermetallic compounds in anodic stripping voltammetry
CN117992860B (en) Method and device for early warning and identifying defects of oil immersed power transformer
CN118378037B (en) Bridge acceleration data denoising method based on long-short-term memory depth neural network
Trifunov et al. Time Series Causal Link Estimation under Hidden Confounding using Knockoff Interventions
EP2350628B1 (en) Method and apparatus for on-line estimation and forecasting of species concentration during a reaction
JP2821580B2 (en) Neural network applied driving support system
Baklouti et al. Monitoring of wastewater treatment plants using improved univariate statistical strategy
Selisteanu et al. Some on adaptive control of a wastewater biodegradation process
González-Garcı́a et al. Characterization of a two-parameter mixed-mode electrochemical behavior regime using neural networks
JPH0463758B2 (en)