KR20220119939A - Apparatus and method for estimating main feed water flow rate - Google Patents

Apparatus and method for estimating main feed water flow rate Download PDF

Info

Publication number
KR20220119939A
KR20220119939A KR1020210023434A KR20210023434A KR20220119939A KR 20220119939 A KR20220119939 A KR 20220119939A KR 1020210023434 A KR1020210023434 A KR 1020210023434A KR 20210023434 A KR20210023434 A KR 20210023434A KR 20220119939 A KR20220119939 A KR 20220119939A
Authority
KR
South Korea
Prior art keywords
estimation
state quantity
model
quantity data
flow rate
Prior art date
Application number
KR1020210023434A
Other languages
Korean (ko)
Inventor
김성근
양학진
Original Assignee
호서대학교 산학협력단
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 호서대학교 산학협력단 filed Critical 호서대학교 산학협력단
Priority to KR1020210023434A priority Critical patent/KR20220119939A/en
Publication of KR20220119939A publication Critical patent/KR20220119939A/en

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/68Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using thermal effects
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/05Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects
    • G01F1/20Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects by detection of dynamic effects of the flow
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/68Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using thermal effects
    • G01F1/684Structural arrangements; Mounting of elements, e.g. in relation to fluid flow
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N5/00Computing arrangements using knowledge-based models
    • G06N5/04Inference or reasoning models
    • 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

Abstract

An estimation apparatus is provided. The estimation apparatus may include: an acquisition unit that acquires state quantity data of a turbine cycle; and an estimation unit that estimates a main feed-water flow for each time domain using the state quantity data. The present invention is to provide the estimation apparatus and estimation method capable of accurately estimating the main feed-water flow of a turbine cycle.

Description

주급수 유량 추정 장치 및 추정 방법{Apparatus and method for estimating main feed water flow rate}Apparatus and method for estimating main feed water flow rate

본 발명은 터빈 사이클의 주급수 유량을 추정하는 장치 및 방법에 관한 것이다.The present invention relates to an apparatus and method for estimating the main feedwater flow rate of a turbine cycle.

터빈 사이클의 성능 상태량을 결정하기 위한 열 성능 분석은 발전소의 향상된 경제성 운전을 위해 요구된다.Thermal performance analysis to determine the amount of performance state of a turbine cycle is required for improved economic operation of a power plant.

열 성능 분석에는 터빈 사이클을 형성하는 단위 기기의 노후화 등의 각종 인자가 반영될 수 있다. 일부 인자의 반영은 특정 단위 기기의 성능 평가의 정확도를 저하시키는 요인으로 작용할 수 있다.In the thermal performance analysis, various factors such as aging of unit devices forming a turbine cycle may be reflected. The reflection of some factors may act as a factor that lowers the accuracy of performance evaluation of a specific unit device.

한국등록특허공보 제1185581호에는 원전의 주급수 제어 밸브 전환시 증기발생기 수위를 제어하는 주급수 제어 시스템이 나타나 있다.Korean Patent Publication No. 1185581 discloses a main water supply control system that controls the water level of a steam generator when switching the main water control valve of a nuclear power plant.

한국등록특허공보 제1185581호Korean Patent Publication No. 1185581

본 발명은 터빈 사이클의 주급수 유량을 정확하게 추정할 수 있는 추정 장치 및 추정 방법을 제공하기 위한 것이다.An object of the present invention is to provide an estimating apparatus and an estimating method capable of accurately estimating the main feedwater flow rate of a turbine cycle.

본 발명의 추정 장치는 터빈 사이클의 상태량 데이터를 획득하는 획득부; 상기 상태량 데이터를 이용하여 주급수 유량(main feed-water flow)을 시간 영역별로 추정하는 추정부;를 포함할 수 있다.An estimating apparatus of the present invention includes: an acquisition unit for acquiring state quantity data of a turbine cycle; and an estimator for estimating a main feed-water flow for each time domain using the state quantity data.

본 발명의 추정 방법은 터빈 사이클의 상태량 데이터 중에서 상관 관계가 설정값을 만족하는 인자를 입력받는 커널 회귀 모델을 이용하여 주급수 유량을 시간 영역별로 추정할 수 있다.The estimation method of the present invention may estimate the main feedwater flow rate for each time domain using a kernel regression model that receives a factor satisfying a set value for correlation among state quantity data of a turbine cycle.

본 발명의 추정 방법은 터빈 사이클의 상태량 데이터를 획득하는 단계; 상태량 데이터의 상관 계수를 계산하여 강인성을 부여하는 단계; 상기 강인성을 비교하여 상기 상태량 데이터의 형상(feature)을 분류하는 단계; 분류된 상기 형상과 분류 테이블을 이용하여 상기 터빈 사이클의 주급수 유량을 추정하는 추정 모델의 인자를 선정하는 단계; 선정된 상기 인자와 상기 추정 모델을 이용하여 상기 주급수 유량을 시간 영역별로 추정하는 단계;를 포함할 수 있다.The estimation method of the present invention includes: acquiring state quantity data of a turbine cycle; calculating a correlation coefficient of the state quantity data to give robustness; classifying a feature of the state quantity data by comparing the toughness; selecting a factor of an estimation model for estimating the main feedwater flow rate of the turbine cycle using the classified shape and classification table; It may include; estimating the main water flow rate for each time domain using the selected factor and the estimation model.

본 발명에서는 유용하고 정확한 성능 분석을 위해서 산업 표준인 ASME PTC를 기본으로 하여 성능 데이터를 사용하여 주급수 유량의 영역별 판정 알고리즘을 개발하고 각 영역별 추정 알고리즘을 개발하였다. 추정 알고리즘은 측정 상태량의 상관관계를 기반으로 형상 분류를 제시하고, 이를 기반으로 커널 회귀모델을 이용하여 추정 모델을 구성하였으며, 서포트 벡터 머신 모델링의 우수성을 검증하기 위하여 신경 회로망 모델과 비교하였다. 주급수 유량의 형상 분류 및 추정 모델은 터빈 사이클에서 정확한 보정 열 성능 분석을 제공함으로써 향상된 성능 분석에 기여할 것이다.In the present invention, for useful and accurate performance analysis, a determination algorithm for each area of the main water flow rate was developed using performance data based on the industry standard ASME PTC, and an estimation algorithm for each area was developed. The estimation algorithm suggested shape classification based on the correlation of measured state quantities, and based on this, an estimation model was constructed using a kernel regression model, and compared with a neural network model to verify the superiority of support vector machine modeling. The shape classification and estimation model of main feedwater flow will contribute to improved performance analysis by providing accurate calibrated thermal performance analysis in turbine cycle.

도 1은 본 발명의 추정 장치를 나타낸 개략도이다.
도 2는 본 발명의 추정 방법을 나타낸 흐름도이다.
도 3은 터빈 사이클에 해당하는 발전소 계통도를 나타낸 개략도이다.
도 4는 열평형 성능(heat balance performance)의 계산 절차를 나타낸 개략도이다.
도 5는 터빈 사이클의 동작 열평형을 계산하는 절차를 나타낸 개략도이다.
도 6은 주급수 유량의 경향을 나타낸 개략도이다.
도 7은 모델부에서 학습 데이터를 이용해 추정 모델을 생성하는 과정을 나타낸 개략도이다.
도 8은 빅데이터를 이용해 추출된 상관 인자를 나타낸 표이다.
도 9는 서포트 벡터 머신의 알고리즘을 나타낸 개략도이다.
도 10은 서포트 벡터 머신의 계산 결과를 나타낸 그래프이다.
도 11은 서포트 벡터 머신에 대한 오류 비교 결과를 나타낸 표이다.
도 12는 신경 회로망(ANN)과 서포트 벡터 머신에 대한 오류 비교 결과를 나타낸 표이다.
도 13은 본 발명의 실시예에 따른, 컴퓨팅 장치를 나타내는 도면이다.
1 is a schematic diagram showing an estimation apparatus of the present invention.
2 is a flowchart illustrating an estimation method of the present invention.
3 is a schematic diagram showing a power plant system diagram corresponding to a turbine cycle.
4 is a schematic diagram showing a calculation procedure of heat balance performance.
5 is a schematic diagram illustrating a procedure for calculating the operating thermal equilibrium of a turbine cycle.
6 is a schematic diagram showing the trend of the main feed water flow rate.
7 is a schematic diagram illustrating a process of generating an estimation model using training data in the model unit.
8 is a table showing correlation factors extracted using big data.
9 is a schematic diagram illustrating an algorithm of a support vector machine.
10 is a graph showing the calculation results of the support vector machine.
11 is a table showing error comparison results for support vector machines.
12 is a table showing error comparison results for a neural network (ANN) and a support vector machine.
13 is a diagram illustrating a computing device according to an embodiment of the present invention.

아래에서는 첨부한 도면을 참고로 하여 본 발명의 실시예에 대하여 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 자가 용이하게 실시할 수 있도록 상세히 설명한다. 그러나 본 발명은 여러 가지 상이한 형태로 구현될 수 있으며 여기에서 설명하는 실시예에 한정되지 않는다. 그리고 도면에서 본 발명을 명확하게 설명하기 위해서 설명과 관계없는 부분은 생략하였으며, 명세서 전체를 통하여 유사한 부분에 대해서는 유사한 도면 부호를 붙였다. Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings so that those of ordinary skill in the art to which the present invention pertains can easily implement them. However, the present invention may be embodied in several different forms and is not limited to the embodiments described herein. And in order to clearly explain the present invention in the drawings, parts irrelevant to the description are omitted, and similar reference numerals are attached to similar parts throughout the specification.

본 명세서에서, 동일한 구성요소에 대해서 중복된 설명은 생략한다.In this specification, duplicate descriptions of the same components will be omitted.

또한 본 명세서에서, 어떤 구성요소가 다른 구성요소에 '연결되어' 있다거나 '접속되어' 있다고 언급된 때에는, 그 다른 구성요소에 직접적으로 연결되어 있거나 또는 접속되어 있을 수도 있지만, 중간에 다른 구성요소가 존재할 수도 있다고 이해되어야 할 것이다. 반면에 본 명세서에서, 어떤 구성요소가 다른 구성요소에 '직접 연결되어' 있다거나 '직접 접속되어' 있다고 언급된 때에는, 중간에 다른 구성요소가 존재하지 않는 것으로 이해되어야 할 것이다.Also, in this specification, when a component is referred to as 'connected' or 'connected' to another component, it may be directly connected or connected to the other component, but other components in the middle It should be understood that there may be On the other hand, in the present specification, when it is mentioned that a certain element is 'directly connected' or 'directly connected' to another element, it should be understood that the other element does not exist in the middle.

또한, 본 명세서에서 사용되는 용어는 단지 특정한 실시예를 설명하기 위해 사용되는 것으로써, 본 발명을 한정하려는 의도로 사용되는 것이 아니다.In addition, the terms used herein are used only to describe specific embodiments, and are not intended to limit the present invention.

또한 본 명세서에서, 단수의 표현은 문맥상 명백하게 다르게 뜻하지 않는 한, 복수의 표현을 포함할 수 있다. Also in this specification, the singular expression may include the plural expression unless the context clearly dictates otherwise.

또한 본 명세서에서, '포함하다' 또는 '가지다' 등의 용어는 명세서에 기재된 특징, 숫자, 단계, 동작, 구성요소, 부품, 또는 이들을 조합한 것이 존재함을 지정하려는 것일 뿐, 하나 또는 그 이상의 다른 특징, 숫자, 단계, 동작, 구성요소, 부품 또는 이들을 조합한 것의 존재 또는 부가 가능성을 미리 배제하지 않는 것으로 이해되어야 할 것이다.Also, in this specification, terms such as 'include' or 'have' are only intended to designate that the features, numbers, steps, operations, components, parts, or combinations thereof described in the specification exist, and one or more It should be understood that the existence or addition of other features, numbers, steps, acts, components, parts, or combinations thereof, is not precluded in advance.

또한 본 명세서에서, '및/또는' 이라는 용어는 복수의 기재된 항목들의 조합 또는 복수의 기재된 항목들 중의 어느 항목을 포함한다. 본 명세서에서, 'A 또는 B'는, 'A', 'B', 또는 'A와 B 모두'를 포함할 수 있다.Also in this specification, the term 'and/or' includes a combination of a plurality of listed items or any of a plurality of listed items. In this specification, 'A or B' may include 'A', 'B', or 'both A and B'.

또한 본 명세서에서, 본 발명의 요지를 흐리게 할 수 있는 공지 기능 및 구성에 대한 상세한 설명은 생략될 것이다.Also, in this specification, detailed descriptions of well-known functions and configurations that may obscure the gist of the present invention will be omitted.

원자력 발전 열성능 분석 시스템은 사이클 열평형을 계산할 수 있으며, 신속하고 정확한 단위기기의 성능 계산이 이루어져야 한다. 즉, 성능분석 프로그램은 필요한 계산 모듈과 단위기기들의 특성에 대한 정보를 내장하고 개별기기들의 모델을 구축하며, 모델끼리 연계·통합함으로서 개별기기의 성능변화가 전체 시스템에 미치는 영향을 쉽고 정확하게 분석할 수 있어야 한다. 이러한 상용화된 주요 발전설비용 프로그램은 보일러와 증기터빈 사이클의 열평형 계산을 주목적으로 개발된 PEPSE(Performance Evaluation of Power System Efficiency), 가스터빈, 복합 및 열병합발전설비의 다양한 운전조건에서 성능계산이 가능한 GR-PRO, GT-MASTER, PEPSE의 기능과 유사한 STEAM-PRO, STEAM-MASTER, 유연성이 높은 Thermoflex, Thermoflex와 유사한 패키지인 GateCycle 등이 있다.Nuclear power generation thermal performance analysis system can calculate cycle thermal equilibrium, and performance calculation of unit equipment must be performed quickly and accurately. In other words, the performance analysis program embeds necessary calculation modules and information on the characteristics of unit devices, builds models of individual devices, and connects and integrates models to easily and accurately analyze the impact of performance changes of individual devices on the overall system. should be able These commercialized programs for major power generation facilities are PEPSE (Performance Evaluation of Power System Efficiency), which was developed for the main purpose of calculating the thermal equilibrium of boiler and steam turbine cycle, and it is possible to calculate performance under various operating conditions of gas turbine, combined cycle and cogeneration facilities. There are GR-PRO, GT-MASTER, STEAM-PRO, STEAM-MASTER similar to PEPSE, Thermoflex with high flexibility, and GateCycle, a package similar to Thermoflex.

PEPSE 모델은 영광 3, 4호기 운전 열성능 모듈의 기본 자료로 활용하고 있으며, 미측정 항목 및 열성능에 대한 예측 모듈 개발에 활용하여 열평형 분석을 위한 성능 계산 절차 및 분석 시스템을 개발하였다. 이러한 터빈 사이클 성능 분석은 ASME(American Society of Mechanical Engineers) PTC(Performance Test Code) 6, PTC 6A, PTC 12를 토대로 터빈 사이클 열평형 계산 및 보정 계산을 하게 된다. 본 발명에서는 이러한 절차를 기반으로 보정 성능을 정확히 결정할 수 있도록, 측정 데이터의 입력 검증 모델을 개발하여, 열 성능 분석시스템의 성능을 업그레이드하고자 한다.The PEPSE model is used as basic data for the thermal performance module of Yeonggwang Units 3 and 4, and the performance calculation procedure and analysis system for thermal equilibrium analysis were developed by using it to develop predictive modules for unmeasured items and thermal performance. The turbine cycle performance analysis is based on ASME (American Society of Mechanical Engineers) PTC (Performance Test Code) 6, PTC 6A, and PTC 12 to calculate turbine cycle thermal balance and correction. The present invention intends to upgrade the performance of the thermal performance analysis system by developing an input verification model of the measurement data so that the correction performance can be accurately determined based on this procedure.

도 1은 본 발명의 추정 장치를 나타낸 개략도이다.1 is a schematic diagram showing an estimation apparatus of the present invention.

열 성능 분석시스템의 성능을 업그레이드하기 위한 일환으로, 도 1에 도시된 추정 장치는 획득부(110), 분류부(130), 선정부(150), 추정부(170), 모델부(190)를 포함할 수 있다.As part of upgrading the performance of the thermal performance analysis system, the estimation apparatus shown in FIG. 1 includes an acquisition unit 110 , a classification unit 130 , a selection unit 150 , an estimation unit 170 , and a model unit 190 . may include.

획득부(110)는 터빈 사이클의 상태량 데이터를 획득할 수 있다.The acquisition unit 110 may acquire state quantity data of the turbine cycle.

터빈 사이클을 발전소의 계통도 또는 발전 시스템의 계통도에 해당하는 것으로, 발전에 필요한 각종 단위 기기, 예를 들어 보일러, 터빈, 복수기, 저압 급수 가열기, 고압 급수 가열기, 각 요소를 연결하는 관, 밸브 등을 포함할 수 있다.A turbine cycle corresponds to a schematic diagram of a power plant or a schematic diagram of a power generation system, and includes various unit devices required for power generation, such as boilers, turbines, condensers, low-pressure feedwater heaters, high-pressure feedwater heaters, and pipes and valves connecting each element. may include

상태량 데이터는 센서를 이용해 각 단위 기기를 대상으로 측정된 측정값, 측정값을 이용해서 산출된 산출값을 포함할 수 있다. 일 예로, 측정값은 유량, 온도, 엔탈피 등을 포함할 수 있다. 측정값은 발전소별로 차이가 있지만 대량 1000개 이상 획득될 수 있다. 산출값은 터빈 효율, 습분 분리기의 습분 분리율, 복수기의 진공도, 급수 가열기의 온도 성능 등의 성능 정보를 포함할 수 있다.The state quantity data may include a measured value measured for each unit device by using a sensor, and a calculated value calculated using the measured value. For example, the measured value may include flow rate, temperature, enthalpy, and the like. Measurements vary from plant to plant, but can be obtained in bulk of 1000 or more. The calculated value may include performance information such as turbine efficiency, a moisture separation rate of a moisture separator, a vacuum degree of a condenser, and a temperature performance of a feedwater heater.

획득부(110)는 센서로부터 상태량 데이터를 입수하거나, 키보드 등의 입력 수단을 이용한 사용자의 입력에 의해 상태량 데이터를 입수할 수 있다. 획득부(110)는 센서 또는 입력 수단과 유무선 통신하는 각종 통신 모듈 또는 인터페이스를 포함할 수 있다.The acquisition unit 110 may obtain the state quantity data from the sensor or may obtain the state quantity data by a user input using an input means such as a keyboard. The acquisition unit 110 may include various communication modules or interfaces for wired/wireless communication with sensors or input means.

분류부(130)는 상태량 데이터의 상관 관계를 기반으로 형상(feature)을 분류할 수있다.The classification unit 130 may classify a feature based on the correlation of the state quantity data.

모델부(190)는 분류된 형상을 기반으로 커널 회귀 모델을 이용하여 주급수 유량(main feed-water flow)을 추정하는 추정 모델을 생성할 수 있다. The model unit 190 may generate an estimation model for estimating a main feed-water flow using a kernel regression model based on the classified shape.

분류부(130) 또는 모델부(190)는 서포트 벡터 머신 모델링을 이용할 수 있다.The classification unit 130 or the model unit 190 may use support vector machine modeling.

구체적으로, 분류부(130)는 서포트 벡터 머신(SVM, Support Vector Machine) 모델링을 이용하여 주급수 유량의 형상을 분류할 수 있다. 모델부(190)는 서포트 벡터 머신 모델링을 이용하여 추정 모델을 생성할 수 있다. 이때, 서포트 벡터 머신 모델링은 선형 함수를 사용할 수 있다.Specifically, the classification unit 130 may classify the shape of the main water flow rate by using a support vector machine (SVM) modeling. The model unit 190 may generate an estimated model using support vector machine modeling. In this case, the support vector machine modeling may use a linear function.

추정부(170)는 상태량 데이터를 이용하여 주급수 유량(main feed-water flow)을 시간 영역별로 추정할 수 있다. 일 예로, 추정부(170)는 모델부(190)에서 생성된 추정 모델을 이용해서 주급수 유량을 시간 영역별로 추정할 수 있다.The estimator 170 may estimate a main feed-water flow for each time domain by using the state quantity data. As an example, the estimator 170 may estimate the main water flow rate for each time domain using the estimation model generated by the model unit 190 .

터빈 사이클과 관련된 상태량 데이터의 종류가 매우 많은 관계로, 추정 모델에서 사용하는 인자는 일부의 상태량 데이터로 제한될 수 있다. 주급수 유량을 정확하게 추정하기 위해 추정 모델에서 사용되는 인자를 선별할 필요가 있다. 인자의 선별 또는 선정에 선정부(150)가 사용될 수 있다.Since there are many types of state quantity data related to the turbine cycle, factors used in the estimation model may be limited to some state quantity data. In order to accurately estimate the main water flow rate, it is necessary to select the factors used in the estimation model. The selection unit 150 may be used for selection or selection of factors.

일 예로, 선정부(150)는 복수의 상기 상태량 데이터 중에서 상기 추정 모델의 인자를 설정 개수만큼 선정할 수 있다. 선정부(150)는 분류부(130)에서 수행된 형상 분류를 기초로 상태량 데이터 중에서 추정 모델의 인자를 선정할 수 있다.For example, the selector 150 may select a set number of factors of the estimation model from among the plurality of state quantity data. The selection unit 150 may select a factor of the estimation model from among the state quantity data based on the shape classification performed by the classification unit 130 .

분류부(130)는 모든 상태량 데이터의 상관 계수를 계산하여 각 상태량 데이터에 강인성을 부여할 수 있다. 분류부(130)는 상태량 데이터의 강인성을 비교하여 형상을 분류할 수 있다.The classification unit 130 may calculate correlation coefficients of all state quantity data to impart robustness to each state quantity data. The classification unit 130 may classify the shape by comparing the toughness of the state quantity data.

선정부(150)는 분류부(130)에서 분류된 형상과 빅데이터를 이용해 작성된 도 8의 분류 테이블을 이용하여 추정 모델의 인자를 그룹 단위로 판별할 수 있다. 선정부(150)는 판별된 특정 그룹을 포함하는 상태량 데이터를 추정 모델의 인자로 선정할 수 있다.The selection unit 150 may determine the factors of the estimation model in groups by using the classification table of FIG. 8 created using the shapes classified by the classification unit 130 and big data. The selection unit 150 may select the state quantity data including the determined specific group as a factor of the estimation model.

본 발명의 추정 방법은 도 1에 도시된 추정 장치에 의해 수행될 수 있다.The estimation method of the present invention may be performed by the estimation apparatus shown in FIG. 1 .

본 발명의 추정 방법은 터빈 사이클의 상태량 데이터 중에서 상관 관계가 설정값을 만족하는 인자를 입력받는 커널 회귀 모델을 이용하여 주급수 유량을 시간 영역별로 추정할 수 있다.The estimation method of the present invention may estimate the main feedwater flow rate for each time domain using a kernel regression model that receives a factor satisfying a set value for correlation among state quantity data of a turbine cycle.

도 2는 본 발명의 추정 방법을 나타낸 흐름도이다.2 is a flowchart illustrating an estimation method of the present invention.

도 2에 도시된 추정 방법은 추정 장치에 의해 수행될 수 있다.The estimation method shown in FIG. 2 may be performed by an estimation apparatus.

터빈 사이틀의 상태량 데이터를 획득하는 단계(S 510)가 수행될 수 있다. 해당 단계는 획득부(110)에 의해 수행될 수 있다.A step ( S510 ) of obtaining the state quantity data of the turbine cycle may be performed. The corresponding step may be performed by the acquisition unit 110 .

상태량 데이터의 상관 계수를 계산하여 강인성을 부여하는 단계(S 520)가 수행될 수 있다. 해당 단계는 분류부(130)에 의해 수행될 수 있다.A step ( S520 ) of providing robustness by calculating a correlation coefficient of the state quantity data may be performed. The corresponding step may be performed by the classification unit 130 .

강인성을 비교하여 상태량 데이터의 형상(feature)을 분류하는 단계(S 530)가 수행될 수 있다. 해당 단계는 분류부(130)에 의해 수행될 수 있다.A step ( S530 ) of classifying features of the state quantity data by comparing the toughness may be performed. The corresponding step may be performed by the classification unit 130 .

분류된 형상과 분류 테이블을 이용하여 터빈 사이클의 주급수 유량을 추정하는 추정 모델의 인자를 선정하는 단계(S 540)가 수행될 수 있다. 해당 단계는 선정부(150)에 의해 수행될 수 있다.Using the classified shape and classification table to select a factor of the estimation model for estimating the main feedwater flow rate of the turbine cycle (S540) may be performed. The corresponding step may be performed by the selector 150 .

선정된 인자와 추정 모델을 이용하여 주급수 유량을 시간 영역별로 추정하는 단계(S 550)가 수행될 수 있다. 해당 단계는 추정부(170)에 의해 수행될 수 있다.The step of estimating the main water flow rate for each time domain using the selected factor and the estimation model (S550) may be performed. This step may be performed by the estimator 170 .

발전소의 계통도를 살펴보면(도 3), 연료의 연소열을 이용하여 터빈이 필요로 하는 증기를 발생하는 보일러(증기 발생기), 보일러에서 나온 고온고압의 증기가 갖고 있는 열에너지를 일로 전환시키는 터빈, 증기를 냉각응축시켜 물로 회수하는 복수기, 복수 펌프에서 공급된 응축수를 보일러에서 다시 가열하여 증기로 변환하는 저압급수가열기, 저압급수가열기와 비슷한 역할을 하는 고압급수가열기을 순환하면서 발전기 출력을 생성하는 사이클로 구성되어 있다. 이러한 터빈 사이클의 열평형(Heat Balance) 분석은 유량, 온도, 엔탈피 등의 상태량을 결정하여 발전 효율, 열 소비율(Heat Rate) 및 기기 성능을 분석하기 위한 것으로 발전소별 차이가 있으나 유동 트레인별로 약 250~300개의 압력, 터빈 사이클 내의 상태량 데이터가 얻어지는데, 발전소별로 차이가 있지만 약 1,500개 정도의 상태량 데이터가 얻어지며, 이 데이터를 토대로 터빈 사이클의 대표 성능 지수인 열 소비율과 터빈 사이클 효율이 얻어진다. 그리고, 터빈 사이클을 형성하는 각 단위 기기에 대해서는 터빈의 효율, 습분 분리기의 습분 분리율, 복수기의 진공도, 급수 가열기의 온도 성능 등의 성능 데이터가 얻어진다(도 4). 성능 데이터는 측정값에 해당하는 상태량 데이터를 통해 산출될 수 있다.Looking at the schematic diagram of the power plant (FIG. 3), a boiler (steam generator) that generates the steam required by the turbine using the combustion heat of fuel, a turbine that converts the thermal energy of high-temperature and high-pressure steam from the boiler into work, and steam A cycle that generates generator output while circulating a condenser that recovers water by cooling and condensing, a low-pressure feedwater heater that converts the condensed water supplied from the condensate pump into steam by heating it again in a boiler, and a high-pressure feedwater heater that plays a similar role to the low-pressure feedwater heater Consists of. Heat balance analysis of such a turbine cycle is to analyze power generation efficiency, heat rate, and device performance by determining state quantities such as flow rate, temperature, and enthalpy. ~300 pressure and state quantity data in the turbine cycle are obtained. Although there is a difference for each power plant, about 1,500 state quantity data are obtained. Based on this data, the heat consumption rate and turbine cycle efficiency, which are representative performance indices of the turbine cycle, are obtained. . Then, for each unit device forming the turbine cycle, performance data such as the efficiency of the turbine, the moisture separation rate of the moisture separator, the vacuum degree of the condensers, and the temperature performance of the feedwater heater are obtained ( FIG. 4 ). The performance data may be calculated through state quantity data corresponding to the measured value.

본 명세서에서 성능 데이터는 상태량 데이터에 포함될 수 있다.In this specification, performance data may be included in state quantity data.

열평형 계산은 터빈 사이클에 질량 및 에너지 보존 법칙(Law of Conservation of Mass and Energy)을 기본개념으로 하여 사이클의 성능 변수(열소비율 및 발전기 출력)를 예측함과 동시에 사이클 내에 모든 지점에 대하여 유동 유체의 상태량들을 결정한다. 에너지 평형의 오차를 10 ~ 15개의 독립변수로 구성된 목적함수로 놓고 이를 최소화하는 최적화 문제로 구성되며 터빈 사이클은 모든 발전소를 포함하는 최대구조의 네트워크 모델을 사용하여 구성하게 된다. 계산을 반복하는 과정 중에 독립변수의 공간에서 다음 위치를 설정하는 알고리즘은 ASME PTC에서 제공되는 휴리스틱과 터빈 추기단 압력에 대한 선형 및 비선형 회귀 모델을 적용하여 구성된다.Thermal equilibrium calculation is based on the Law of Conservation of Mass and Energy in the turbine cycle to predict the performance variables (heat consumption rate and generator output) of the cycle, and at the same time, the flow fluid for all points in the cycle. determine the state quantities of It is composed of an optimization problem that minimizes the energy balance error by setting it as an objective function composed of 10 to 15 independent variables, and the turbine cycle is constructed using a network model of the maximum structure including all power plants. The algorithm for setting the next position in the space of the independent variable during the process of repeating the calculation is constructed by applying the heuristic provided by ASME PTC and linear and non-linear regression models for the turbine extractor pressure.

운전 열평형 계산을 통하여 산출된 터빈 사이클의 성능변수는 터빈 사이클의 현재 성능을 나타내며, 이는 터빈 자체뿐만 아니라 터빈 사이클 내 부속기기의 영향이 복합적으로 나타난 것이므로 터빈 자체의 현재 성능을 파악하는 데는 어려움이 있다. 이에 따라 터빈의 현재 성능을 파악할 수 있는 별도의 보정 계산이 필요하게 된다. 보정 계산은 터빈을 제외한 사이클내 주요기기들의 성능을 설계조건 기준으로 변환시켜 주면서 사이클 열평형 계산을 행하는 그룹 Ⅰ 보정과 터빈 제작사가 제공하는 보정 곡선을 이용하여 외부로부터 터빈에 직접적으로 영향을 미치는 항목에 대하여 각각의 영향을 개별적으로 보정해주는 그룹 Ⅱ 보정으로 이루어진다. 성능 분석에서는 이러한 보정을 설계 성능 모듈을 사용하는 대신 커브 피팅(Curve Fitting)을 사용해서 처리한다.The performance variable of the turbine cycle calculated through the calculation of the operating thermal equilibrium represents the current performance of the turbine cycle, and it is difficult to understand the current performance of the turbine itself because the effects of the turbine itself as well as the accessories in the turbine cycle are compounded. have. Accordingly, a separate calibration calculation is required to determine the current performance of the turbine. Calibration calculation is an item that directly affects the turbine from the outside using the group I calibration, which calculates cycle thermal equilibrium while converting the performance of major devices in the cycle to the design condition standard, except for the turbine, and the calibration curve provided by the turbine manufacturer. It consists of a group II correction that compensates for each effect individually. In performance analysis, these corrections are handled using Curve Fitting instead of using the Design Performance module.

열평형 수렴문제는 터빈 팽창선 끝점을 독립 변수로 하고 열평형 오차 또는 습분 분리율 오차를 목적 함수로 하는 최적화 문제가 되므로 몰리에르 차트 상에서 팽창선의 곡률을 이용한 근사해법을 응용하여 도 5와 같은 계산 절차를 적용하였다. 저압 터빈 ELEP(Expansion Line End Point)의 반복 계산이 이루어지며, 새로운 팽창선에서 계산 ELEP와 추정 ELEP의 차가 ASME PTC의 허용 편차 0.1 Btu/lb를 초과하면 추정 ELEP 대신 계산 ELEP를 이용하여 터빈 팽창선을 재구성하고 반복 열평형 계산 과정에 들어가게 된다. 이러한 열평형 시스템의 계산 절차는 이미 많은 검증이 이루어졌다.The thermal equilibrium convergence problem is an optimization problem using the end point of the turbine expansion line as an independent variable and the thermal equilibrium error or moisture separation rate error as an objective function. did. Iterative calculation of the low pressure turbine Expansion Line End Point (ELEP) is made, and if the difference between the calculated ELEP and the estimated ELEP at the new expansion line exceeds the ASME PTC allowable deviation of 0.1 Btu/lb, the calculated ELEP is used instead of the estimated ELEP to reconstruct the turbine expansion line and enters the iterative thermal equilibrium calculation process. The calculation procedure of such a thermal equilibrium system has already been verified a lot.

도 6의 (a)는 터빈 사이클의 성능 분석을 위하여 현장에서 측정하는 데이터 변수 중에서 주요 항목인 발전기 출력, 주급수 유량의 장기간(약 160주기)에 걸친 실제 변화이다. 이러한 발전기 출력은 터빈 사이클 효율과 함께 터빈 사이클 성능을 분석할 수 있는 중요한 지표로서 터빈이 공급하는 일량에서 베어링 마찰 손실 등과 같은 터빈까지 발생하는 각종 손실과 발전기 절연 손실 등과 같은 발전기에서 발생하는 각종 손실을 감하여 구해지지만, 각 항들은 복잡한 식으로 이루어져서 실제 계산식에 의한 결과 도출은 쉽지 않다. 현장 설비의 노후화 같은 환경 요인에 의하여 이상 측정 상태가 발생할 수가 있다. 예를 들어, 운전 기간이 경과함에 따라 오리피스에 생기는 불순물의 축척 등으로 실제 값보다 큰 값이 측정될 수 있다. 이에 따라 열 출력이 과다 계산되므로 안전 기준에 따라 운전 제어 변수를 조정해야 하기 때문에, 열 출력이 실제 가능한 최대값보다 축소된 상태에서 운전되는 비효율적인 결과를 만들어 낸다. 이와 같은 상황이 지속된다고 판단되면 계측기 보정으로 유량 측정을 강제 보정하게 되며, 그 후에는 보정된 계측에 따라 성능을 측정하고 분석하는 것이 일반적이다. 도 6의 (a)는 주급수 유량의 동향에 따른 발전기 출력량을 보여주고 있다. 특별한 흐름을 보이지 않다가 앞에서 언급한 것처럼 측정 오류로 인한 강제 보정에 따른 변동 구간이 발생한다. 하지만 발전기 출력의 동향은 뚜렷한 변화가 포착되지 않으나 강제 보정 이후 영역에서 급증하는 것을 알 수 있다.Figure 6 (a) is the actual change over a long period of time (about 160 cycles) of the main items of the generator output and main feed water flow rate among the data variables measured in the field for the performance analysis of the turbine cycle. This generator output is an important indicator to analyze turbine cycle performance along with turbine cycle efficiency. It is an important indicator to analyze the various losses that occur from the work supplied by the turbine to the turbine, such as bearing friction loss, and various losses that occur in the generator, such as the generator insulation loss. However, it is not easy to derive the result by the actual calculation formula because each term is composed of a complex formula. An abnormal measurement condition may occur due to environmental factors such as deterioration of field equipment. For example, a value larger than the actual value may be measured due to the accumulation of impurities generated in the orifice as the operation period elapses. As a result, heat output is over-calculated, and operation control variables must be adjusted according to safety standards, resulting in inefficient operation with heat output reduced from the maximum possible value. If it is determined that such a situation continues, the flow rate measurement is forcibly corrected by calibrating the instrument. After that, it is common to measure and analyze the performance according to the calibrated instrument. Figure 6 (a) shows the generator output according to the trend of the main feed water flow rate. Although it does not show any special flow, as mentioned above, there is a variation section due to forced correction due to measurement error. However, it can be seen that there is no clear change in the generator output trend, but it sharply increases in the area after the forced correction.

도 6의 (b)와 같이, 주급수 유량과의 상관관계가 높은 3개의 인자를 추출하여, 신경 회로망 모델의 입력 인자로 설정하고 학습시킨 후 추정값을 살펴보았다. 정상 상태에서의 평균 오차는 1.52%이며, 최대 오차가 발생하는 추정값도 3.87% 로서 ASME PTC 6의 데이터 유효기준이내에 만족하고 있지만, 강제 보정 구간 및 보정 이후 구간을 포함한 전체 평균 오차는 3.16%로서, 몇몇 데이터를 제외한 강제 보정 구간의 데이터에서는 유효기준을 만족하지 못하고 있으며, 보정 구간 이후에서도 이러한 오차 데이터가 발생함을 알 수 있었다. 이러한 이유는 강제 보정 구간에서의 높은 오차는 정상 상태의 학습 데이터에 따른 신경 회로망의 일반화 능력이 미약하기 때문에 발생하는 것으로 판단되며, 보정 이후 구간에서는 강제 보정 구간의 데이터를 학습 데이터로 활용하면서 학습 데이터와 추정 값의 상관성이 약화되어서 발생하는 것으로 판단된다.As shown in (b) of FIG. 6 , three factors having a high correlation with the main water flow rate were extracted, set as input factors of the neural network model, and after training, the estimated values were examined. The average error in the steady state is 1.52%, and the estimated value at which the maximum error occurs is 3.87%, which is within the data validity criteria of ASME PTC 6, but the overall average error including the forced correction section and the post-correction section is 3.16%, It was found that the data in the forced correction section except for some data did not satisfy the validity criteria, and this error data occurred even after the correction section. The reason for this is that the high error in the forced correction section is judged to occur because the generalization ability of the neural network according to the learning data in the steady state is weak. It is judged that this occurs because the correlation between and the estimated value is weakened.

상관 인자 선정시 고려해야 할 사항은 (1) 인자의 수가 많아지면 계산이 복잡해져서 학습 및 인식 속도가 느려지며, (2) 불필요한 인자는 학습을 어렵게 만든다. (3) 또한 부적합한 인자의 추가로 인하여 모델의 정밀도가 악화될 수 있다.Things to consider when selecting a correlation factor are: (1) if the number of factors increases, the calculation becomes complicated, slowing down learning and recognition, and (2) unnecessary factors make learning difficult. (3) Also, the precision of the model may deteriorate due to the addition of inappropriate factors.

일반적으로 복잡한 데이터를 단순한 데이터로 정제하는 단계로서 추정 모델의 학습 단계 이전의 전처리 과정으로, 원용 데이터(Raw Data)에 대한 적절한 가지치기 작업을 해야 한다. In general, it is a step of refining complex data into simple data. It is a pre-processing process before the training step of the estimation model. Appropriate pruning of raw data should be performed.

데이터 축소란 대량의 변수 군을 소량의 변수 군으로 병합하는 과정을 의미하는 데, 일반적으로 많이 사용하는 인자 축소 방식인 주성분 분석은 인자들의 상관 분석을 통하여 고유 벡터를 도출하고 큰 고유 벡터를 사상시킴으로서 축소된 인자를 추출하는 방식이다. 하지만 원용 데이터에서 모든 데이터를 사용할 경우 부적합한 인자도 적용되기 때문에 모델의 정밀도는 악화될 수 있다. Data reduction refers to the process of merging a large number of variable groups into a small number of variable groups. Principal component analysis, which is a commonly used factor reduction method, derives eigenvectors through correlation analysis of factors and maps large eigenvectors. This is a method of extracting reduced factors. However, when all data is used in the raw data, the precision of the model may deteriorate because inappropriate factors are also applied.

상관 인자의 선정은 크게 (1) 모델을 기반으로 하는 방식(Model-based Method)과 (2) 모델을 기반으로 하지 않는 방식(Model-free Method)으로 나눌 수 있다. 모델 기반 방식은 선택된 인자를 이용하여 학습을 통한 최적의 설정 변수를 찾는 방식으로서, 인자를 변경시키면서 반복 과정에 의하여 최적의 인자를 선택하는 과정이다. 그리고 모델을 기반으로 하지 않는 방식은 학습 데이터의 인자와 이에 대한 출력의 통계적 관계를 바탕으로 수행하는 방식이다. The selection of correlation factors can be largely divided into (1) a model-based method and (2) a model-free method. The model-based method is a method of finding an optimal set variable through learning using a selected factor, and is a process of selecting an optimal factor through an iterative process while changing the factor. And the method that is not model-based is a method that is performed based on the statistical relationship between the factors of the training data and the output.

본 발명에서는 상관 인자와 출력의 강도를 알아보기 위하여 피어슨 상관 계수를 사용하여 변수간의 영향력을 파악하였다.In the present invention, the influence between variables was identified using the Pearson correlation coefficient in order to examine the strength of the correlation factor and the output.

데이터 기반 모델에서의 학습 단계란 상관 인자로 구성된 학습 데이터로부터 객체의 특색(Traits) 혹은 특징(Feature)을 정의하는 과정으로서 학습을 통하여 속이 빈(Blank) 벡터를 완전한 특징 벡터로 만드는 단계이며, 이러한 벡터는 학습 데이터의 반복 계산에 의하여 이루어지기 때문에 학습 데이터의 영향력이 크다. 이렇게 완성된 특징 벡터를 이용하여 예측값을 추정하게 되며, 추정값이 실제 값과 흡사할 때 일반화 능력이 우수하다고 한다. 데이터 기반 모델에 대한 일반화 능력은 특징 벡터에 의하여 핸들링되며, 학습 데이터에 의존하게 된다. 즉 추정 값에 대한 상관 인자의 데이터들이 학습 데이터의 특징과 비슷한 경향이 있어야 실제 값과의 오차를 줄일 수 있다고 판단된다. The learning stage in the data-based model is the process of defining the traits or features of an object from the learning data composed of correlation factors. Since the vector is made by iterative calculation of the training data, the influence of the training data is large. The predicted value is estimated using the completed feature vector, and it is said that the generalization ability is excellent when the estimated value is similar to the actual value. The generalization ability for data-driven models is handled by feature vectors and depends on the training data. That is, it is judged that the error with the actual value can be reduced when the data of the correlation factor for the estimated value have a tendency similar to the characteristics of the training data.

본 발명에서는 추정값에 대한 추세를 알아보고, 이에 대한 강인한 상관 인자로 조합된 학습 데이터를 구성하고자 한다. 추정 값에 대한 경향은 기존에 수집된 데이터를 이용하여 추정 항목에 대한 특징 영역을 분류하고, 추정값의 우수한 상관 인자와 특징 영역에서의 강인한 상관 인자를 비교하여, 추정값에 대한 특징 영역, 즉 형상을 파악하는 것이다(도 7). 강인한 상관 인자란 도 8처럼 기존의 빅데이터를 이용하여 미리 추출한 특징 영역의 상관 인자들이다. 이러한 상관 인자들은 형상별 인자 및 강인성을 갖고 있는 인자수가 다르다. 상관 계수가 0.7 이상인 경우, 강인성을 갖고 있는 인자수가 40개가 넘는 경우도 발생하며, 적을 때는 10개 미만인 경우도 나타난다. 상관 인자 수가 많은 경우는 가동 시간이 길어지는 단점이 있어서, 신경회로망 모델에서는 입력 인자수를 정해서 진행하는 연구가 많았다. 본 발명에서는 학습 데이터를 이용하여 모든 인자의 상관 계수를 계산하고, 인자의 강인성(0, 1)을 비교하여 형상을 구별하며, 인자의 순위별 그룹(Low, Middle, High)을 정하여, 해당 그룹까지 포함되는 인자를 입력 인자로 설정한다. 예를 들어 학습 데이터의 모든 인자별 상관 계수를 계산하여, 강인성을 부여하고 도 8과 비교하였을 때, #5 형상으로 판별되었고, 임의의 인자수를 3으로 하였을 때 순위에 포함되는 그룹은 Low 그룹이다. 그런 경우에는 Low 그룹까지 포함하는 인자(factor)는 모두 입력 인자로 설정될 수 있다.In the present invention, the trend of the estimated value is investigated, and training data combined with a strong correlation factor therefor is to be constructed. The trend for the estimated value is determined by classifying the feature region for the estimated item using previously collected data, and comparing the excellent correlation factor of the estimated value with the strong correlation factor in the feature region to determine the feature region for the estimated value, that is, the shape. to understand (FIG. 7). The strong correlation factors are correlation factors of a feature region previously extracted using existing big data, as shown in FIG. 8 . These correlation factors have different factors for each shape and the number of factors having toughness. When the correlation coefficient is 0.7 or more, there are cases where the number of factors with toughness exceeds 40, and when the correlation coefficient is small, there are cases where there are fewer than 10. In the case of a large number of correlation factors, the operation time is long, so there have been many studies in which the number of input factors is determined in the neural network model. In the present invention, correlation coefficients of all factors are calculated using learning data, shapes are distinguished by comparing the toughness (0, 1) of factors, and groups (Low, Middle, High) are determined by rank of factors, and the corresponding group Sets the argument including up to as the input argument. For example, when the correlation coefficient for each factor of the learning data is calculated, toughness is given, and when compared with FIG. 8, it is determined as #5 shape, and when the arbitrary factor number is 3, the group included in the ranking is the Low group to be. In such a case, all factors including the Low group may be set as input factors.

서포트 벡터 머신은 두 부류를 나누는 기준의 여백(margin)을 최대화하여 일반화 능력을 최대로 하는 것을 목적으로 한다. 이는 추가 발생된 임의의 패턴을 분류할 때에 오분류될 가능성이 작아진다는 의미이다. 기본적으로 SVM은 이진 분류를 위하여 초평면(Hyper-plane)을 기준으로 분류하며 다음과 같은 수학식1로 나타낼 수 있다.The support vector machine aims to maximize the generalization ability by maximizing the margin of the criterion that divides the two classes. This means that the probability of misclassification is reduced when classifying any additionally generated patterns. Basically, the SVM classifies based on a hyper-plane for binary classification, and it can be expressed by Equation 1 below.

Figure pat00001
Figure pat00001

X는 특징 벡터로서 X=(x1,...,xd)T이며, W는 법선 벡터, b는 중심에서 초평면까지의 거리이다.X is the feature vector, where X=(x 1 ,...,x d ) T , W is the normal vector, and b is the distance from the center to the hyperplane.

여기서 초평면 d(X)를 기준으로 0보다 크고 작음에 따라 특징 벡터의 분할이 가능하다. 임의의 점 x에서 초평면까지의 거리는

Figure pat00002
로 나타낼 수 있다.Here, it is possible to divide the feature vector as it is greater than or less than 0 based on the hyperplane d(X). The distance from any point x to the hyperplane is
Figure pat00002
can be expressed as

초평면으로부터 여백의 크기를 좌우하는 벡터를 서프트 벡터(Support Vector)라 한다. 여백(Margin)은 초평면으로부터 가장 가까운 벡터까지의 거리에 2배 값을 가지므로 다음의 수학식 2와 같이 나타낼 수 있다.A vector that determines the size of the blank space from the hyperplane is called a support vector. Since the margin has a value twice the distance from the hyperplane to the nearest vector, it can be expressed as in Equation 2 below.

Figure pat00003
Figure pat00003

조건부 최적화 문제로 최대 여백을 갖는 결정 초평면을 구할 수 있다.As a conditional optimization problem, a crystal hyperplane with maximum margin can be obtained.

Figure pat00004
의 조건일 때,
Figure pat00005
를 최대화하여야 한다.
Figure pat00004
When the condition of
Figure pat00005
should be maximized.

비용함수 J(.)를 사용하여 바꾸어 나타내면,If expressed using the cost function J(.),

Figure pat00006
의 조건일 때,
Figure pat00007
을 최소화하는 것과 같다.
Figure pat00006
When the condition of
Figure pat00007
is equivalent to minimizing

따라서, J(W)를 최소로 하는 W를 구해야 하며, 이는 라그랑제 승수(Lagrange multiplier)를 이용하여 구한다. 조건식마다 라그랑제 승수 ai를 넣어 그 벡터를Therefore, it is necessary to find W that minimizes J(W), which is obtained using a Lagrange multiplier. Put the Lagrange multiplier a i in each conditional expression to obtain the vector.

Figure pat00008
라 표기하고, 라그랑제 함수를 L(.)로 두면, 다음의 수학식 3과 같이 정의될 수 있다.
Figure pat00008
, and the Lagrange function is set to L(.), it can be defined as in Equation 3 below.

Figure pat00009
Figure pat00009

라그랑제 함수로 표현한 조건부 최적화 문제는 KKT(Karush-Kuhn-Tucker) 조건을 이용하여 풀 수 있다. KKT 조건에 의하여, 수학식 4와 같이 나타낼 수 있다.The conditional optimization problem expressed by the Lagrange function can be solved using the Karush-Kuhn-Tucker (KKT) condition. According to the KKT condition, it can be expressed as Equation (4).

Figure pat00010
Figure pat00010

여기서,

Figure pat00011
일 때,
Figure pat00012
을 만족하는 벡터가 서포터 벡터이다. KKT조건의 첫 번째와 두 번째 식을 라그랑제 함수에 대입하여 풀면 W와 b가 사라져 라그랑제 승수만 가지게 되는데 이를
Figure pat00013
로 두고 다시 정리하면,here,
Figure pat00011
when,
Figure pat00012
A vector that satisfies is a supporter vector. By substituting the first and second expressions of the KKT condition into the Lagrange function, W and b disappear and only the Lagrange multiplier is obtained.
Figure pat00013
If you put it back and rearrange it,

Figure pat00014
를 최대화시키는 것으로 정리될 수 있다.
Figure pat00014
It can be arranged by maximizing .

최대 마진을 갖는 초평면을 찾는 순서는 최대값을 얻는 값을 찾고, W,b순으로 구하여 d(X)를 얻는다. 듀얼 최적화 모델과 최적 서포트 벡터 머신 회귀 모델이 적용되었고, 대표적인 커널 함수인 방사 기저 함수 커널(Radial Basis Function Kernel)을 적용하였다.The order of finding the hyperplane with the maximum margin is to find the value that obtains the maximum value, and then obtain d(X) in the order of W and b. A dual optimization model and an optimal support vector machine regression model were applied, and a radial basis function kernel, a representative kernel function, was applied.

서포트 벡터 머신 모델은 의 범위 내에서 적절한 값을 설정하는 것으로서, 가장 간단한 방법은 학습시 사용한 학습 데이터의 예측오차를 최소로 하는 값을 찾는 것이다. 하지만, 이러한 설정은 과대적합(Over-fitting) 문제가 발생하여 학습되지 않은 데이터에 대한 예측성능이 오히려 저하되는 현상이 발생한다는 것이 알려져 있다. 이와 같은 현상을 방지하고 일반화 성능을 개선시키기 위하여 여러 가지 기법이 제시되고 있는데, 대표적인 방법으로 VC 차원 기법, 베이시안 학습기법, 교차검증기법 등이 있다. 이들 중에서 상대적으로 간단하면서도 우수한 성능을 나타내는 방법이 교차검증이다. 교차검증으로, k-folds 방식, Leave-One-Out 방식, Bootstrap 방식 등이 있다. 본 발명에서는 교차검증기법 중에서 Leave-One-Out 기법을 사용하여 최적의 λ를 선정하고자 하며 도 9에 도시된 다음과 같은 절차에 의하여 이루어진다.The support vector machine model sets an appropriate value within the range of , and the simplest method is to find a value that minimizes the prediction error of the training data used for training. However, it is known that such a setting causes an over-fitting problem, resulting in a decrease in predictive performance on untrained data. In order to prevent such a phenomenon and improve generalization performance, various techniques have been proposed. Representative methods include the VC dimension technique, Bayesian learning technique, and cross-validation technique. Among them, cross-validation is a relatively simple and excellent method. As for cross-validation, there are k-folds method, Leave-One-Out method, Bootstrap method, and the like. In the present invention, the optimal λ is selected by using the Leave-One-Out method among the cross-validation methods, and is performed by the following procedure shown in FIG. 9 .

(1) 모든 학습 데이터만큼 임의의 λ를 초기화시킨다.(1) Initialize arbitrary λ as much as all training data.

(2) 학습 데이터 중에서 하나의 데이터를 하위 검증 데이터(Sub-validation data)로 귀속시키고, 나머지 데이터를 하위 학습 데이터(Sub-learning data)로 이용하여 데이터를 랜덤으로 재배치시킨다. 또한, 하위 학습 데이터에 상속된 λ는 상수화하고, 하위 검증 데이터에 상속된 λ는 변수화하여 계산한다.(2) One data among the training data is attributed to sub-validation data, and the data is randomly rearranged using the remaining data as sub-learning data. In addition, λ inherited from the sub-learning data is constant, and λ inherited from the sub-verification data is calculated as a variable.

(3) KKT 조건에 맞는 y(x) 및 w0를 계산하고, 설정 오차와 비교한다.(3) Calculate y(x) and w 0 according to the KKT condition and compare it with the setting error.

(4) 설정 오차에 의하여 만족하지 못할 때는 심플렉스 방법(Simplex Method)을 이용하여 하위 검증 데이터에 상속된 λ를 재설정한다.(4) If it is not satisfied due to the setting error, the λ inherited from the lower verification data is reset using the Simplex Method.

여기서, 하위 검증 데이터의 상속된 λ가 만족하면 (2) 항목의 하위 학습 데이터와 하위 검증 데이터를 재설정한다. 모든 데이터와 상관된 λ를 갱신할 때까지 (2) ~ (4)을 반복하고 검증 데이터를 통한 학습 평가를 한다.Here, if the inherited λ of the sub-verification data is satisfied, the sub-learning data and sub-verification data of item (2) are reset. Repeat (2) to (4) until λ correlated with all data is updated, and learning evaluation is performed through validation data.

서포트 벡터 머신 모델에 의한 추정 값은 도 10과 같으며, 평균 오차는 도 11처럼 나타났다. 모든 커널 함수들이 대체적으로 오차가 1% 미만으로 나타났으며, 측정값을 정확하게 추정하는 값(0.1% 미만)들도 상당수가 있었다. 하지만 몇 개의 경우는 다소 큰 오차가 발생하였다. 먼저, 선형 커널 함수로 구현하였을 때, #1 형상에서 1% 이상인 오차율을 갖는 추정 값은 2 주기이고, #2 형상 및 #5 형상에서는 모두 1 주기가 나타났으며, #3 형상 및 #4 형상에서는 모든 추정 값이 %1 미만으로 수행하였다. 특히 #2 형상에서의 오차율이 상대적으로 너무 큰 값을 갖고 있어서, 평균 오차율이 높아진 것을 알 수 있었으며, #5 형상에서는 1% 이상인 오차율이 1 주기가 나타났지만, 다른 추정 값에서 신뢰성이 높았기 때문에, 형상별 평균 오차율이 가장 작게 나타날 수 있었다.The estimated value by the support vector machine model is shown in FIG. 10, and the average error is shown as shown in FIG. 11. All the kernel functions showed an error of less than 1% in general, and there were many values that accurately estimated the measured value (less than 0.1%). However, in some cases, a rather large error occurred. First, when implemented as a linear kernel function, the estimated value with an error rate of 1% or more in the #1 shape is 2 cycles, and 1 cycle appears in both shapes #2 and #5, and shapes #3 and #4 In , all estimated values were less than %1. In particular, it can be seen that the error rate in shape #2 has a relatively large value, so the average error rate is increased. , the average error rate for each shape was the smallest.

다항식 함수로 구현했을 때, #1 형상 및 #5 형상에서 1% 이상인 오차율을 갖는 추정 값은 3 주기이고, #2 형상에서는 2 주기, #3 형상에서는 1주기가 나타났다. 이렇듯, 1% 이상인 오차율을 갖는 추정 값이 다른 함수보다 많아서, 전체 평균 오차율이 가장 높았지만, #2 형상의 시작 구간에서는 가장 낮은 오차율을 보였다.When implemented as a polynomial function, the estimated values with an error rate of 1% or more for shape #1 and #5 were 3 cycles, 2 cycles for shape #2, and 1 cycle for shape #3. As such, there were more estimated values with an error rate of 1% or more than other functions, so the overall average error rate was the highest, but the lowest error rate was shown in the starting section of shape #2.

RBF 함수로 구현하였을 때는, #1 형상에서 1% 이상인 오차율을 갖는 추정 값은 1 주기이고, #2 형상 및 #5 형상에서는 모두 2 주기가 나타났으며, #3 형상 및 #4 형상에서는 모든 추정 값이 %1 미만으로 수행하였다. 특히 #2 형상에서의 오차율이 상대적으로 너무 큰 값을 갖고 있어서, 평균 오차율이 높아짐을 알 수 있었다.When implemented as an RBF function, the estimated value with an error rate of 1% or more for shape #1 is 1 period, for shape #2 and shape #5, 2 cycles appeared, and for shape #3 and #4, all estimates Values were less than %1. In particular, it was found that the error rate in shape #2 had a relatively large value, so that the average error rate was increased.

서포트 벡터 머신 모델에서는 선형 함수가 높은 신뢰성을 보임으로서 일반화 능력이 가장 우수하였다.In the support vector machine model, the linear function showed high reliability, so the generalization ability was the best.

도 12는 서포트 벡터 머신 모델을 이용한 추정 모델 결과의 한 예로서, 오차가 크게 발생한 몇몇 데이터도 포함하였다. 이러한 값들은 모든 함수에서 대부분이 #2 형상과 #5 형상에서 나타났는데, 오차가 가장 큰 값은 측정값이 5.9786 Mkg/hr 이었으나 추정 값이 5.91251 Mkg/hr (선형 함수), 5.9244 Mkg/hr (다항식 함수), 5.91382 Mkg/hr (RBF 함수)으로 나타난 곳으로서 강제 보정이 시작되는 영역이었다. 이렇게 크게 발생하는 오차에 의하여 #2 형상에서 평균 오차가 커진 것을 알 수 있었으며, #5 형상에서는 전반적으로 오차가 크게 나타나서 평균 오차율이 커진 것을 알 수 있었다.12 is an example of an estimation model result using a support vector machine model, and some data with large errors are also included. Most of these values appeared in #2 shape and #5 shape in all functions. The measured value was 5.9786 Mkg/hr, but the estimated value was 5.91251 Mkg/hr (linear function), 5.9244 Mkg/hr ( polynomial function), 5.91382 Mkg/hr (RBF function), which is the region where the forced correction starts. Due to this large error, it was found that the average error was increased in the #2 shape, and the overall error was large in the #5 shape, indicating that the average error rate was increased.

본 발명에서는 이러한 측정 데이터의 정확성을 얻기 위하여 유효화 검증을 하였으며, 신뢰성 있는 측정 데이터를 얻을 수 있었다.In the present invention, validation was performed to obtain the accuracy of such measurement data, and reliable measurement data could be obtained.

1) 열평형 분석을 하기 위한 운전 성능 데이터 중에서 측정 데이터와 상관성이 있다고 판단되는 성능 인자를 측정 데이터 및 성능 지침서의 엔지니어링 정보를 활용하여 인자 추출 및 데이터 쌍을 구성하였으며, 상관관계를 통하여 우수한 상관 인자 및 데이터 쌍을 추출하였다.1) Among the driving performance data for thermal equilibrium analysis, the performance factors judged to be correlated with the measured data were extracted and data pairs were constructed using the measured data and engineering information in the performance guide, and excellent correlation factors through correlation and data pairs were extracted.

2) 기존의 빅데이터에서 데이터 집단의 특징을 활용하여 형상을 분리할 수 있었으며 각각의 형상별 강인한 상관 인자를 도출시켜 형상에 맞는 강인한 상관 인자를 활용하여 윈도우를 설정하였다. 또한 적절한 윈도우 수는 각각의 인자별 데이터 쌍에 대한 상관계수를 통하여 파악할 수 있었으며, 이렇게 추출된 윈도우를 활용하여, 추정 모델에 맞는 학습 데이터로 재구성할 수 있었다. 2) In the existing big data, shapes could be separated by utilizing the characteristics of the data group, and a strong correlation factor for each shape was derived, and a window was set using the strong correlation factor suitable for the shape. In addition, the appropriate number of windows could be identified through the correlation coefficient for each data pair for each factor, and the extracted windows could be used to reconstruct the training data suitable for the estimation model.

3) 서포트 벡터 머신 모델에서는 선형 함수가 높은 신뢰성을 보임으로서 일반화 능력이 가장 우수하다는 것을 알 수 있었다.3) In the support vector machine model, the linear function showed high reliability, indicating that the generalization ability was the best.

4) 도 12에서처럼 신경 회로망보다 서포트 벡터 머신 함수의 추정 값이 신뢰할 수 있었다. 형상별 집입 구간에서는 이전 형상에 대한 학습 데이터로 추정값에 반영하지만 윈도우가 지나가면서 빠르게 오차를 줄여서 추정값을 도출하는 것을 알 수 있었다.4) As shown in FIG. 12, the estimated value of the support vector machine function was more reliable than that of the neural network. In the collection section for each shape, the learning data for the previous shape is reflected in the estimated value, but it was found that the estimated value was derived by rapidly reducing the error as the window passed.

도 13은 본 발명의 실시예에 따른, 컴퓨팅 장치를 나타내는 도면이다. 도 13의 컴퓨팅 장치(TN100)는 본 명세서에서 기술된 장치(예, 추정 장치 등) 일 수 있다. 13 is a diagram illustrating a computing device according to an embodiment of the present invention. The computing device TN100 of FIG. 13 may be a device (eg, an estimation device, etc.) described herein.

도 13의 실시예에서, 컴퓨팅 장치(TN100)는 적어도 하나의 프로세서(TN110), 송수신 장치(TN120), 및 메모리(TN130)를 포함할 수 있다. 또한, 컴퓨팅 장치(TN100)는 저장 장치(TN140), 입력 인터페이스 장치(TN150), 출력 인터페이스 장치(TN160) 등을 더 포함할 수 있다. 컴퓨팅 장치(TN100)에 포함된 구성 요소들은 버스(bus)(TN170)에 의해 연결되어 서로 통신을 수행할 수 있다.In the embodiment of FIG. 13 , the computing device TN100 may include at least one processor TN110 , a transceiver device TN120 , and a memory TN130 . Also, the computing device TN100 may further include a storage device TN140 , an input interface device TN150 , an output interface device TN160 , and the like. Components included in the computing device TN100 may be connected by a bus TN170 to communicate with each other.

프로세서(TN110)는 메모리(TN130) 및 저장 장치(TN140) 중에서 적어도 하나에 저장된 프로그램 명령(program command)을 실행할 수 있다. 프로세서(TN110)는 중앙 처리 장치(CPU: central processing unit), 그래픽 처리 장치(GPU: graphics processing unit), 또는 본 발명의 실시예에 따른 방법들이 수행되는 전용의 프로세서를 의미할 수 있다. 프로세서(TN110)는 본 발명의 실시예와 관련하여 기술된 절차, 기능, 및 방법 등을 구현하도록 구성될 수 있다. 프로세서(TN110)는 컴퓨팅 장치(TN100)의 각 구성 요소를 제어할 수 있다.The processor TN110 may execute a program command stored in at least one of the memory TN130 and the storage device TN140. The processor TN110 may mean a central processing unit (CPU), a graphics processing unit (GPU), or a dedicated processor on which methods according to an embodiment of the present invention are performed. The processor TN110 may be configured to implement procedures, functions, and methods described in connection with an embodiment of the present invention. The processor TN110 may control each component of the computing device TN100 .

메모리(TN130) 및 저장 장치(TN140) 각각은 프로세서(TN110)의 동작과 관련된 다양한 정보를 저장할 수 있다. 메모리(TN130) 및 저장 장치(TN140) 각각은 휘발성 저장 매체 및 비휘발성 저장 매체 중에서 적어도 하나로 구성될 수 있다. 예를 들어, 메모리(TN130)는 읽기 전용 메모리(ROM: read only memory) 및 랜덤 액세스 메모리(RAM: random access memory) 중에서 적어도 하나로 구성될 수 있다. Each of the memory TN130 and the storage device TN140 may store various information related to the operation of the processor TN110 . Each of the memory TN130 and the storage device TN140 may be configured as at least one of a volatile storage medium and a nonvolatile storage medium. For example, the memory TN130 may include at least one of a read only memory (ROM) and a random access memory (RAM).

송수신 장치(TN120)는 유선 신호 또는 무선 신호를 송신 또는 수신할 수 있다. 송수신 장치(TN120)는 네트워크에 연결되어 통신을 수행할 수 있다.The transceiver TN120 may transmit or receive a wired signal or a wireless signal. The transceiver TN120 may be connected to a network to perform communication.

한편, 본 발명의 실시예는 지금까지 설명한 장치 및/또는 방법을 통해서만 구현되는 것은 아니며, 본 발명의 실시예의 구성에 대응하는 기능을 실현하는 프로그램 또는 그 프로그램이 기록된 기록 매체를 통해 구현될 수도 있으며, 이러한 구현은 상술한 실시예의 기재로부터 본 발명이 속하는 기술 분야의 통상의 기술자라면 쉽게 구현할 수 있는 것이다. On the other hand, the embodiment of the present invention is not implemented only through the apparatus and/or method described so far, and a program for realizing a function corresponding to the configuration of the embodiment of the present invention or a recording medium in which the program is recorded may be implemented. And, such an implementation can be easily implemented by those skilled in the art from the description of the above-described embodiments.

이상에서 본 발명의 실시예에 대하여 상세하게 설명하였지만 본 발명의 권리범위는 이에 한정되는 것은 아니고 다음의 청구범위에서 정의하고 있는 본 발명의 기본 개념을 이용한 통상의 기술자의 여러 변형 및 개량 형태 또한 본 발명의 권리범위에 속하는 것이다.Although the embodiments of the present invention have been described in detail above, the scope of the present invention is not limited thereto, and various modifications and improvements by those skilled in the art using the basic concept of the present invention as defined in the following claims are also present. It belongs to the scope of the invention.

110...획득부 130...분류부
150...선정부 170...추정부
190...모델부
110...Acquisition Department 130...Classification Department
150...selection unit 170...estimation unit
190...Model Department

Claims (10)

터빈 사이클의 상태량 데이터를 획득하는 획득부;
상기 상태량 데이터를 이용하여 주급수 유량(main feed-water flow)을 시간 영역별로 추정하는 추정부;
를 포함하는 추정 장치.
an acquisition unit for acquiring state quantity data of a turbine cycle;
an estimator for estimating a main feed-water flow for each time domain using the state quantity data;
Estimation device comprising a.
제1항에 있어서,
상기 상태량 데이터의 상관 관계를 기반으로 형상(feature)을 분류하는 분류부, 분류된 상기 형상을 기반으로 커널 회귀 모델을 이용하여 상기 주급수 유량을 추정하는 추정 모델을 생성하는 모델부가 마련되고,
상기 추정부는 상기 추정 모델을 이용해서 상기 주급수 유량을 시간 영역별로 추정하는 추정 장치.
According to claim 1,
A classification unit for classifying a feature based on the correlation of the state quantity data, a model unit for generating an estimation model for estimating the main feedwater flow rate using a kernel regression model based on the classified shape is provided,
The estimation unit is an estimation device for estimating the main water flow rate for each time domain using the estimation model.
제2항에 있어서,
상기 분류부는 서포트 벡터 머신 모델링을 이용하여 상기 주급수 유량의 형상을 분류하는 추정 장치.
3. The method of claim 2,
The classification unit is an estimation device for classifying the shape of the main feed water flow rate using support vector machine modeling.
제2항에 있어서,
상기 모델부는 서포트 벡터 머신 모델링을 이용하여 상기 추정 모델을 생성하는 추정 장치.
3. The method of claim 2,
and the model unit generates the estimated model by using support vector machine modeling.
제2항에 있어서,
상기 분류부 또는 상기 모델부는 서포트 벡터 머신 모델링을 이용하고,
상기 서포트 벡터 머신 모델링은 선형 함수를 사용하는 추정 장치.
3. The method of claim 2,
The classification unit or the model unit uses support vector machine modeling,
The support vector machine modeling is an estimation device using a linear function.
제1항에 있어서,
상기 추정부는 추정 모델을 이용해서 상기 주급수의 유량을 추정하고,
복수의 상기 상태량 데이터 중에서 상기 추정 모델의 인자를 설정 개수만큼 선정하는 선정부가 마련된 추정 장치.
According to claim 1,
The estimator estimates the flow rate of the main water supply using an estimation model,
An estimator having a selector configured to select a set number of factors of the estimation model from among the plurality of state quantity data.
제1항에 있어서,
상기 추정부는 추정 모델을 이용해서 상기 주급수의 유량을 추정하고,
상기 상태량 데이터의 상관 관계를 기반으로 형상(feature)을 분류하는 분류부, 상기 형상 분류를 기초로 상기 상태량 데이터 중에서 상기 추정 모델의 인자를 선정하는 선정부가 마련된 추정 장치.
According to claim 1,
The estimator estimates the flow rate of the main water supply using an estimation model,
An estimation apparatus provided with a classification unit for classifying a feature based on the correlation of the state quantity data, and a selection unit for selecting a factor of the estimation model from among the state quantity data based on the shape classification.
제7항에 있어서,
상기 분류부는 모든 상태량 데이터의 상관 계수를 계산하여 강인성을 부여하며,
상기 분류부는 상기 상태량 데이터의 강인성을 비교하여 상기 형상을 분류하고,
상기 선정부는 분류된 상기 형상과 분류 테이블을 이용하여 상기 추정 모델의 인자를 그룹 단위로 판별하며,
상기 선정부는 판별된 특정 그룹을 포함하는 상태량 데이터를 상기 추정 모델의 인자로 선정하는 추정 장치.
8. The method of claim 7,
The classification unit gives robustness by calculating the correlation coefficient of all state quantity data,
The classification unit classifies the shape by comparing the toughness of the state quantity data,
The selector determines the factors of the estimation model in groups using the classified shape and classification table,
The selection unit is an estimation device for selecting the state quantity data including the determined specific group as a factor of the estimation model.
추정 장치에 의해 수행되는 추정 방법에 있어서,
터빈 사이클의 상태량 데이터 중에서 상관 관계가 설정값을 만족하는 인자를 입력받는 커널 회귀 모델을 이용하여 주급수 유량을 시간 영역별로 추정하는 추정 방법.
An estimation method performed by an estimation device, comprising:
An estimation method for estimating the main feedwater flow rate for each time domain using a kernel regression model that receives factors that satisfy a set value for correlation among the state quantity data of the turbine cycle.
추정 장치에 의해 수행되는 추정 방법에 있어서,
터빈 사이클의 상태량 데이터를 획득하는 단계;
상태량 데이터의 상관 계수를 계산하여 강인성을 부여하는 단계;
상기 강인성을 비교하여 상기 상태량 데이터의 형상(feature)을 분류하는 단계;
분류된 상기 형상과 분류 테이블을 이용하여 상기 터빈 사이클의 주급수 유량을 추정하는 추정 모델의 인자를 선정하는 단계;
선정된 상기 인자와 상기 추정 모델을 이용하여 상기 주급수 유량을 시간 영역별로 추정하는 단계;
를 포함하는 추정 방법.
An estimation method performed by an estimation device, comprising:
acquiring state quantity data of the turbine cycle;
calculating a correlation coefficient of the state quantity data to give robustness;
classifying a feature of the state quantity data by comparing the toughness;
selecting a factor of an estimation model for estimating the main feedwater flow rate of the turbine cycle using the classified shape and classification table;
estimating the main water flow rate for each time domain using the selected factor and the estimation model;
Estimation method comprising
KR1020210023434A 2021-02-22 2021-02-22 Apparatus and method for estimating main feed water flow rate KR20220119939A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020210023434A KR20220119939A (en) 2021-02-22 2021-02-22 Apparatus and method for estimating main feed water flow rate

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020210023434A KR20220119939A (en) 2021-02-22 2021-02-22 Apparatus and method for estimating main feed water flow rate

Publications (1)

Publication Number Publication Date
KR20220119939A true KR20220119939A (en) 2022-08-30

Family

ID=83114313

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020210023434A KR20220119939A (en) 2021-02-22 2021-02-22 Apparatus and method for estimating main feed water flow rate

Country Status (1)

Country Link
KR (1) KR20220119939A (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101185581B1 (en) 2010-12-30 2012-09-24 한국전력기술 주식회사 System of controlling steam generator level during main feed-water control valve transfer for the nuclear power plant

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101185581B1 (en) 2010-12-30 2012-09-24 한국전력기술 주식회사 System of controlling steam generator level during main feed-water control valve transfer for the nuclear power plant

Similar Documents

Publication Publication Date Title
JP6686151B2 (en) Model parameter value estimating device and method, program, recording medium storing program, model parameter value estimating system
CN108492000B (en) Fault diagnosis method for nonstationary characteristic of million-kilowatt ultra-supercritical unit
US7931041B2 (en) System and method for controlling liquid level in a vessel
US9200591B2 (en) Automatic identification of operating parameters for power plants
US7933754B2 (en) System and method for damage propagation estimation
US7890296B2 (en) Method of analyzing the performance of gas turbine engines
CN103759769A (en) Soft measurement method and device for reheat steam flow of coal-fired power plant
US11512848B2 (en) Systems and methods for real-time steam quality estimation
KR20220119939A (en) Apparatus and method for estimating main feed water flow rate
CN112329191A (en) Modal identification method and transition mode fault detection method in multi-mode process
CN102967393B (en) Supercritical unit boiler separator outlet temperature on-line calibration method
JP4008431B2 (en) Turbine rotor inlet temperature synthesis system and synthesis method
JPWO2005041209A1 (en) Nuclear power plant thermal efficiency diagnostic system, nuclear power plant thermal efficiency diagnostic program, and nuclear power plant thermal efficiency diagnostic method
CN115290218A (en) Soft measurement method and system for wall temperature of boiler water wall of thermal generator set
CN116029433A (en) Energy efficiency reference value judging method, system, equipment and medium based on grey prediction
CN115014454A (en) Soft measurement method, system, equipment and medium for main steam flow of thermal power generating unit
CN112734158B (en) Thermoelectric load distribution method and device of generator set and electronic equipment
KR102603020B1 (en) Plant monitoring devices, plant monitoring methods, and programs
Gu et al. A modified clustering procedure for energy consumption monitoring in the steam turbine considering volume effect
JP7142545B2 (en) Boiler tube leak diagnostic system and boiler tube leak diagnostic method
Kim et al. Development of performance analysis system (NOPAS) for turbine cycle of nuclear power plant
Tadić et al. Fault diagnosis for steam separators based on parameter identification and CUSUM classification
WO2022181574A1 (en) Abnormality factor estimation method for power plant
US11886157B2 (en) Operational optimization of industrial steam and power utility systems
JP2004171425A (en) Plant reliability evaluation method and performance evaluation system

Legal Events

Date Code Title Description
E902 Notification of reason for refusal