KR20110043277A - A method of failure prediction at power plant using bayesian statics - Google Patents

A method of failure prediction at power plant using bayesian statics Download PDF

Info

Publication number
KR20110043277A
KR20110043277A KR1020090100320A KR20090100320A KR20110043277A KR 20110043277 A KR20110043277 A KR 20110043277A KR 1020090100320 A KR1020090100320 A KR 1020090100320A KR 20090100320 A KR20090100320 A KR 20090100320A KR 20110043277 A KR20110043277 A KR 20110043277A
Authority
KR
South Korea
Prior art keywords
failure
rank
year
value
power generation
Prior art date
Application number
KR1020090100320A
Other languages
Korean (ko)
Other versions
KR101086027B1 (en
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 KR1020090100320A priority Critical patent/KR101086027B1/en
Publication of KR20110043277A publication Critical patent/KR20110043277A/en
Application granted granted Critical
Publication of KR101086027B1 publication Critical patent/KR101086027B1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • 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/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Economics (AREA)
  • Data Mining & Analysis (AREA)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • General Business, Economics & Management (AREA)
  • Mathematical Analysis (AREA)
  • Health & Medical Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Marketing (AREA)
  • Operations Research (AREA)
  • Pure & Applied Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Databases & Information Systems (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Game Theory and Decision Science (AREA)
  • Software Systems (AREA)
  • Probability & Statistics with Applications (AREA)
  • Quality & Reliability (AREA)
  • Development Economics (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Algebra (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Testing And Monitoring For Control Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

PURPOSE: A method of failure prediction at a power plant using Bayesian statistics is provided to efficiently maintain and manage a power facility by determining the efficient repair and maintenance period based on Bayesian statistics in consideration of a TVFR(Time-Varying Failure Rate) and the maintenance frequency. CONSTITUTION: The information on the failure history of each component is inputted to a database(201,202). Based on the information on the failure history, the years of failure and the number of failure times depending on each failure reason are analyzed(203). Based on the analyzed result, the rank and reverse rank are calculated(204). Through the calculated reverse rank, an adjusted rank and a media rank are calculated(205). In consideration of the calculated media rank and the failure years, regression analysis is performed.

Description

베이지안 통계기법을 이용한 발전설비 고장예측방법{A Method of Failure Prediction at Power Plant Using Bayesian Statics} A Method of Failure Prediction at Power Plant Using Bayesian Statics

본 발명은 발전설비 각 부품들에 대한 고장이력을 데이터베이스화하고, 데이터베이스에 저장된 데이터를 바탕으로 베이지안 통계방법을 이용하여 발전설비 각 부품에 대한 고장을 예측하는 알고리즘 및 개발 프로그램을 설계 제작하되, 상기 발전설비 각 부품의 고장을 예측하는 알고리즘은 베이지안 통계기법 중 하나인 기기 고장 예측에 사용되는 웨이블(Weibull)함수를 본 발명에 적용하고, 이 방법을 이용하여 일정주기로 교체되는 발전설비의 구성품들의 수명을 예측ㆍ진단하여 불필요하게 미리 주기적으로 부품을 교체함에 따른 경제적인 손실을 미연에 방지하여 부품구입에 따른 비용증가 및 발전 중지에 의한 경제적 손실을 방지하고, 발전설비의 신뢰성을 향상시킬 수 있는 베이지안 통계기법을 이용한 발전설비 고장예측 방법에 관한 것이다.The present invention design and manufacture an algorithm and a development program for predicting the failure of each component of the power generation plant by using the Bayesian statistical method based on the data stored in the database of the failure history for each component of the power plant; The algorithm for predicting the failure of each component of the power plant is applied to the present invention the Weibull function, which is one of the Bayesian statistical techniques, for the equipment failure prediction. By predicting and diagnosing the service life, it is possible to prevent the economic loss caused by replacing the parts periodically and unnecessary in advance, thereby preventing the cost increase due to the purchase of parts and the economic loss caused by stopping the power generation, and improving the reliability of the power generation equipment. The present invention relates to a power plant failure prediction method using Bayesian statistical technique.

통상적인 기계제품에서는 기능설계, 강도설계, 최적화 설계방법 등을 통하여 신뢰성과 안정성을 확보할 수 있으나, 다품종 소량으로 생산되는 제품에서는 통계학 처리에 필요한 데이터를 충분히 얻지 못하는 경우가 많아 품질의 확률적 특성을 평가하는 것이 어려운 문제점이 있다.In general mechanical products, reliability and stability can be secured through functional design, strength design, and optimization design methods. However, in the case of products produced in small quantities, many products do not have enough data for statistical processing. There is a difficult problem to evaluate.

통상적인 전자제품의 경우는 단품종, 대량생산으로 통계기법에 의한 예측이 가능하지만 기계 부품에 대해서는 신뢰성과 안정성 두 가지를 동시에 고려하여야 하기 때문에 안정성 평가의 잣대인 리스크를 기준으로 리스크 관리기준(Ri나 Based Design), 리스크 기준 보수관리(Risk Based Maintenance)등의 기준으로 평가하여야 한다. In the case of general electronic products, it is possible to predict by statistical technique with single product type and mass production. However, since both reliability and stability should be considered at the same time for mechanical parts, risk management criteria (Ri Or based on criteria such as Based Design) or Risk Based Maintenance.

현재까지 발전설비의 유지보수는 TBM(Time Based Maintenance)에 따라 일정 기간이 경과하면 해당 부품들을 분해 교체 및 정비하는 방법을 시행하여 왔으나, 신속하게 변화하는 발전산업의 주위 환경은 TBM의 유지보수 방법으로는 최대효율과 비용대비 이익을 보장 할 수 없는 문제점이 있다.Until now, maintenance of power generation facilities has been implemented by disassembling, replacing and maintaining the parts after a certain period of time according to TBM (Time Based Maintenance) .However, the environment around the rapidly changing power generation industry is how to maintain TBM. There is a problem that cannot guarantee the maximum efficiency and cost-benefit.

본 발명과 관련된 종래기술로, 대한민국 공개특허공보 제특10-2004-0015011호에는 네트워크혈관 건강 분석평가 시스템 및 방법을 이루기 위하여 사용된 베이지안 통계추론에 대한 기술적 구성이 개시되어 있고, 대한민국 공개특허공보 제10-2006-0044772호에는 베이지안 네트워크를 트리 학습을 위한 테이블 사용 방법에 대한 기술적 구성이 개시되어 있으나, 각각의 발명에서 사용하고 있는 베이지안 추론 및 베이지안 네트워크의 기술적 사상 및 구성은 발전설비 각 부품에 대한 고장을 예측하는 통계기법과는 현저한 차이가 있다.In the prior art related to the present invention, Korean Patent Publication No. 10-2004-0015011 discloses a technical configuration of Bayesian statistical inference used to achieve the network vascular health analysis evaluation system and method, and the Republic of Korea 10-2006-0044772 discloses a technical configuration of a Bayesian network using a table for tree learning, but the technical concept and configuration of Bayesian inference and Bayesian network used in each invention are described for each component of the power plant. There is a significant difference from statistical techniques for predicting failure.

본 발명이 해결하고자 하는 과제는 TBM(Time Based Maintenance) 방식의 발전설비의 정비를 리스크 기준 보수관리(Risk Based Maintenance) 등의 규정에 따라 발전설비를 유지관리하기 위한 것이며, 발전설비 부품들의 고장예측을 위하여 각 부품들의 고장이력을 데이터베이스화하고 이를 이용하여 통계기법 중의 하나인 베이지안 통계방법을 이용하여 발전설비의 구성 부품들의 고장확률을 예측 제공하여 부품의 잦은 교체에 따른 경제적 손실을 줄이고, 교체 시에 발전 중단에 따른 경제적인 손실을 줄이는데 있다.The problem to be solved by the present invention is to maintain the power generation equipment according to the risk-based maintenance, such as maintenance of the power generation equipment of the TBM (Time Based Maintenance) method, the failure prediction of the power generation equipment components For this purpose, the failure history of each part is made into a database and the Bayesian statistical method, which is one of the statistical techniques, is used to predict the failure probability of the components of the power generation facilities to reduce the economic loss due to frequent replacement of parts. This is to reduce the economic loss of stopping the development.

본 발명이 해결하고자 하는 또 다른 과제는 베이지안 통계기법을 기초로 하는 유지보수의 시변고장율(TVFR : Time-Varying Failure Rate) 및 유지보수의 횟수 등을 고려하여 효과적인 유지보수 주기를 결정하고, 리스크 기준 보수관리(Risk Based Maintenance)등의 기준을 설정하여 발전설비를 효율적으로 유지 관리하여 경제적 손실을 최소화하고 신뢰성을 향상시키는데 있다.Another problem to be solved by the present invention is to determine the effective maintenance interval in consideration of the Time-Varying Failure Rate (TVFR) of maintenance based on Bayesian statistical techniques and the number of maintenance, risk criteria By setting standards such as Risk Based Maintenance, it is possible to efficiently maintain power generation facilities to minimize economic losses and improve reliability.

본 발명의 과제해결 수단은 발전설비의 용량, 효율, 운전정보와 운전시간, 고장 및 정비이력을 입력하여 데이터베이스에 저장하고, 부품이나 설비의 고장확률 분포를 작성하며, 이를 바탕으로 웨이블 분포함수에 의하여 고장을 예측하되, 발전설비 각 부품의 고장 및 정비이력을 조사 분석하여 부품별 고장이력에 대한 정보(데이터)를 데이터베이스에 저장하는 단계를 거쳐서, 데이터베이스에 저장된 데이 터를 이용하여 발전설비 부품의 원인 별 각 고장주기, 고장년도 및 고장건수 등을 분석하는 단계를 거쳐서, 분석한 결과를 이용하여 발전설비 각 부품에 대하여 랭크(RANK) 및 역 랭크(REVERSE RANK)를 결정하는 단계를 거쳐서, 역 랭크(REVERSE RANK)를 이용하여 조절랭크(ADJUSTED RANK)와 중앙랭크(MEDIAN RANK)를 결정하는 단계를 거치며, 고장년도와 중앙랭크(MEDIAN RANK)를 이용하여 회귀분석 수행 후 웨이블(WEIBULL) 분포함수의 두 인자인 α,β 값을 결정하는 단계를 거쳐서, 웨이블(WEIBULL) 분포함수를 연산하여 고장예측을 수행하는 단계로 구성된 베이지안 통계기법을 이용한 발전설비 고장예측 방법을 구현하는데 있다.The problem solving means of the present invention inputs the capacity, efficiency, operation information and operation time, failure and maintenance history of the power generation facility in the database, and prepare the probability distribution of parts or equipment, based on the wave distribution function Predict failure by the step, and analyze and analyze the failure and maintenance history of each part of the power generation equipment, and store the information (data) of the failure history for each part in the database, and use the data stored in the database to generate the power generation equipment parts. After analyzing each cycle, failure year, and number of failures of each cause, and determining the rank and reverse rank for each part of the power generation equipment by using the analysis result, REVERSE RANK is used to determine the ADJUSTED RANK and the MEDIAN RANK.The failure year and the MEDIAN Bayesian consisting of a step of calculating the wavelet distribution function and performing the failure prediction by performing the regression analysis using RANK) and determining the α, β values of the two factors of the wavelet distribution function. To implement a method for predicting failure of a power plant using statistical techniques.

본 발명의 또 다른 과제해결 수단은 베이지안 통계기법을 기초로 하는 유지보수의 시변 고장율(TVFR : Time-Varying Failure Rate) 및 유지보수의 횟수 등을 고려하여 효과적인 유지보수 주기를 결정하고, 리스크 기준 보수관리(Risk Based Maintenance)등의 기준을 설정하여 발전설비를 효율적으로 유지 관리하여 경제적 손실을 최소화하고 신뢰성을 향상시킬 수 있는 발전설비 고장예측 방법을 구현하는데 있다.Another problem solving means of the present invention is to determine the effective maintenance interval in consideration of the time-varying failure rate (TVFR) of maintenance and the number of maintenance based on Bayesian statistical techniques, risk-based maintenance It is to implement power generation equipment failure prediction method that can minimize the economic loss and improve the reliability by setting up the criteria such as Risk Based Maintenance.

본 발명은 TBM 방식의 발전설비의 정비를 리스크 기준 보수관리 등의 기준으로 발전설비를 유지 관리하기 위한 것이며, 발전설비의 부품들의 고장예측을 위하여 각 부품들의 고장이력을 데이터베이스화하고 이를 이용하여 통계기법 중의 하나인 베이지안 통계방법으로 발전설비의 구성 부품들의 고장확률을 예측 제공하여 부품의 잦은 교체에 따른 경제적 손실을 줄이고, 부품교체를 위하여 발전을 중단함에 따른 경제적인 손실을 줄이는 작용효과가 있다.The present invention is to maintain the maintenance of the power generation facilities on the basis of the risk-based maintenance management, such as maintenance of the TBM-type power generation facilities, database the failure history of each component to predict the failure of the components of the power generation facilities and statistics Bayesian statistical method, which is one of the techniques, predicts the failure probability of the components of power generation facilities, thereby reducing the economic loss due to frequent replacement of components and reducing the economic loss due to stopping the power generation for parts replacement.

본 발명의 또 다른 효과는 베이지안 통계기법을 기초로 하는 유지보수의 시변고장율 및 유지보수의 횟수 등을 고려하여 효과적인 유지보수 주기를 결정하고, 리스크 기준 보수관리 등의 기준을 설정하여 발전설비를 효율적으로 유지 관리하여 경제적 손실을 최소화하고 신뢰성을 높이는데 있다.Another effect of the present invention is to determine the effective maintenance cycle in consideration of the time-varying failure rate of maintenance and the number of maintenance based on Bayesian statistical techniques, and to set the criteria such as risk-based maintenance management, etc. To minimize economic losses and increase reliability.

본 발명의 실시를 위한 구체적인 내용에 대하여 살펴본다. 본 발명은 발전설비의 용량, 효율, 운전정보와 운전시간, 고장 및 정비이력에 대한 정보를 데이터베이스에 입력 저장하는 단계를 거쳐서, 상기 데이터베이스에 저장된 데이터를 이용하여 발전설비 각 부품의 원인 별 고장주기, 고장년도 및 고장건수 등을 분석하는 단계를 거친다. Hereinafter, the present invention will be described in detail. The present invention through the step of inputting and storing the capacity, efficiency, operation information and operation time, failure and maintenance history information of the power plant in the database, using the data stored in the database failure cycle for each component of the power plant In this case, the system analyzes the failure year and the number of failures.

상기 데이터베이스에 저장된 데이터를 이용하여 발전설비 부품의 원인 별 각 고장주기, 고장년도 및 고장건수 등을 분석하는 단계를 거쳐서, 분석한 결과를 이용하여 발전설비 각 부품에 대하여 랭크(RANK) 및 역 랭크(REVERSE RANK)를 계산하는 단계를 거친다.Analyze each failure cycle, failure year and number of failures of each power plant component by using the data stored in the database. Rank and reverse rank for each power plant component using the analyzed result Calculate the (REVERSE RANK).

상기 역 랭크(REVERSE RANK)를 이용하여 조절랭크(ADJUSTED RANK)와 중앙랭크(MEDIAN RANK)를 계산하는 단계를 거치며, 고장년도와 중앙랭크(MEDIAN RANK)를 이용하여 회귀분석 수행 후 웨이블(WEIBULL) 분포함수의 두 인자인 α,β 값을 연산하는 단계를 거쳐서, 웨이블(WEIBULL) 분포함수를 연산하여 고장예측을 수행하는 단계를 거쳐서 베이지안 통계기법을 이용한 발전설비 고장예측 방법을 구현한다.Calculate the ADJUSTED RANK and the MEDIAN RANK using the REVERSE RANK, and perform the regression analysis using the failure year and the MEDIAN RANK. ) The fault prediction method of power generation facilities using Bayesian statistical technique is implemented through the step of calculating the value of α, β, which is two factors of the distribution function, and the fault prediction by calculating the WEIBULL distribution function.

본 발명은 국내 발전 설비 뿐만 아니라, 다양한 산업설비 등에도 범용으로 적용될 수 있으며, 국내 산업 전반의 경제적인 설비 운영에도 적용할 수 있다. 현재 약 50기의 국내화력발전소와 국내 전반에 걸쳐있는 산업설비 및 완제품의 고장 및 불량확률에 대한 예측을 수행할 수 있으므로 이 기술의 수요성 및 적용분야가 매우 넓다. 국내외 제작사 및 연구소를 중심으로 발전설비 고장예측에 대한 요구와 연구가 활발히 이루어지고 있다.The present invention can be applied not only to domestic power generation facilities, but also to various industrial facilities. Currently, about 50 domestic thermal power plants and industrial facilities and finished products are predicted for the probability of failure and defects throughout the country, so the demand and application of this technology are very wide. There are active demands and researches on power plant failure prediction, mainly from domestic and overseas manufacturers and research institutes.

이하, 첨부된 도면을 참조하여 본 발명의 실시 예의 구성과 작용을 설명하며, 도면에 도시되고 설명되는 본 발명의 구성과 작용은 적어도 하나 이상의 실시 예로서 설명되는 것이며, 이것에 의해 상기한 본 발명의 기술적 사상과 그 핵심 구성 및 작용이 제한되지는 않는다.Hereinafter, the configuration and operation of the embodiments of the present invention with reference to the accompanying drawings, the configuration and operation of the present invention shown and described in the drawings will be described by at least one embodiment, whereby the present invention described above The technical idea and its core composition and operation are not limited.

본 발명의 이해를 용이하게 하는 도면을 살펴본다. 도1은 고장 DB을 구축하는 방법과 이를 이용하여 고장확률 분포를 계산하고, 웨이블(WEIBULL) 분포함수를 이용하여 고장예측을 계산하는 방법을 나타낸 것이며, 도2는 발전설비 각 부품에 대한 고장예측을 위한 알고리즘을 수행하는데 필요한 각 단계별 흐름도를 나타낸 것이다. 도3은 고장예측을 위한 기본 알고리즘을 수행하기 위한 흐름도를 나타낸 것이다. 도4는 본 발명에 따라 설계 제작된 데이터베이스의 구성 중에서 하나의 화면을 캡쳐하여 나타낸 것이며, 도5는 회귀분석을 통한 웨이블 인자를 도출하기 위한 일차함수를 구하는 방법을 도시한 것이다. 본 발명에 따른 구체적인 실시 예를 살펴본다. Look at the drawings to facilitate understanding of the present invention. FIG. 1 shows a method of constructing a failure DB, a method of calculating a probability of failure using this method, and a method of calculating a failure prediction using a wavelet distribution function, and FIG. 2 shows a failure of each component of a power plant. Each step required to perform the algorithm for prediction is shown. 3 shows a flowchart for performing a basic algorithm for fault prediction. 4 is a screen capture of one of the configurations of a database designed and manufactured according to the present invention, and FIG. 5 illustrates a method for obtaining a first-order function for deriving a wave factor through regression analysis. A specific embodiment according to the present invention will be described.

<실시 예><Example>

본 발명에 따른 구체적인 실시 예를 도면에 기초하여 살펴본다. 도4는 본 발명에 따라 설계 제작된 데이터베이스의 한 화면을 캡쳐한 도면이고, 발전설비 부품 중에서 터빈을 선택한 경우의 화면을 나타낸 것이다. 터빈을 선택한 경우에 좌측 화면에는 고압(HP), 중압(IP) 및 저압(LP)로 나누어지고, 각각의 압력에 따른 부품이 나열되어 있으며, 우측에는 좌측에 나열된 각각의 부품에 대한 부품코드, 부품명, 준공년도, 해당 부품의 고장년도, 고장의 종류, 고장원인 및 고장에 대한 조치사항을 데이터베이스에 입력할 수 있도록 구성되어 있다. 상기 데이터베이스에는 발전설비 부품 중에서 터빈 외에 보일러, 발전기 등 발전설비에 채용되는 부품이 모두 해당될 수 있다. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS A preferred embodiment of the present invention will now be described with reference to the drawings. FIG. 4 is a view of capturing a screen of a database designed and manufactured according to the present invention, and shows a screen when a turbine is selected among power plant components. When turbine is selected, the left screen is divided into high pressure (HP), medium pressure (IP) and low pressure (LP), and the parts according to each pressure are listed, and on the right, the part code for each part listed on the left, It is configured to input part name, year of completion, year of failure, type of failure, cause of failure and corrective action in the database. The database may correspond to all components employed in a power generation facility, such as a boiler and a generator, in addition to a turbine.

도2는 발전설비 각 부품에 대한 고장예측을 위한 알고리즘을 수행하는데 필요한 각 단계별 흐름도를 나타낸 것이다. 각 단계별로 살펴보면, 발전설비 각 부품의 고장 및 정비이력을 조사하고 분석하여 부품별 고장이력에 대한 정보를 찾아내어 데이터베이스에 입력하는 단계를 거치고, 데이터베이스에 입력된 각 부품에 대한 정보를 바탕으로 발전설비 부품의 고장원인별 각 고장년도와 고장건수를 분석하는 단계를 거친다.Figure 2 shows a flow chart for each step required to perform the algorithm for failure prediction for each component of the power plant. Looking at each stage, it analyzes and analyzes the failure and maintenance history of each part of the power generation equipment, finds out the information on the failure history of each part, inputs it into the database, and develops based on the information about each part entered in the database. It analyzes each failure year and number of failures by cause of equipment parts.

상기 분석된 결과를 이용하여 랭크(RANK) 및 역 랭크(REVERSE RANK)를 계산하는 단계를 거치며, 결정된 역 랭크(REVERSE RANK)를 이용하여 조절랭크(ADJUSTED RANK)와 중앙랭크(MEDIAN RANK)를 계산하는 단계를 거친다.Calculating the rank (RANK) and the reverse rank (REVERSE RANK) using the analyzed result, and calculates the adjusted rank (ADJUSTED RANK) and the median rank (MEDIAN RANK) using the determined reverse rank (REVERSE RANK) Go through the steps.

상기 결정된 중앙랭크(MEDIAN RANK)와 고장년도를 이용하여 회귀분석 수행하여 웨이블 분포함수의 두 인자인 척도모수(α)와 형상모수(β)를 결정하는 단 계(206)를 거치며, 결정된 척도모수(α)와 형상모수(β)를 웨이블 분포함수에 대입하여 각 부품에 대한 고장예측을 수행하는 단계(207)로 이루어져 있다.The regression analysis is performed using the determined median rank and failure year, and the step 206 of determining the scale parameter (α) and the shape parameter (β), which are two factors of the wavelet distribution function, is determined. Step 207 is performed by substituting the parameter α and the shape parameter β into the wavelet distribution function to perform failure prediction for each component.

상기 회귀분석은 통계학에서 관찰된 연속형 변수들에 대해 독립변수와 종속변수 사이의 인과관계에 따른 수학적 모델인 선형적 관계식을 구하여 어떤 독립변수가 주어졌을 때 이에 따른 종속변수를 예측하는 방법이며, 또한 이 수학적 모델이 얼마나 잘 설명하고 있는지를 판별하기 위한 적합도를 측정하는 분석방법이다.The regression analysis is a method of predicting a dependent variable when a given independent variable is obtained by obtaining a linear relation, which is a mathematical model of a causal relationship between the independent variable and the dependent variable, for the continuous variables observed in statistics. It is also an analytical method that measures fitness to determine how well this mathematical model explains.

본 발명은 상기 데이터베이스에 발전설비 각 부품의 고장이력 및 정보를 입력할 때, 도4에서와 같이 발전설비를 터빈, 발전기, 보일러 등의 큰 부품으로 나누어 구성하고(401), 대상설비(402) 선택 후 부품코드(404) 및 부품명(405)을 입력하며, 해당 부품이 채용된 발전설비의 준공년도(406), 부품의 고장발생년도(407), 고장종류(408), 고장원인(409) 및 조치사항(410)등을 차례로 입력하되, 정확하게 입력하여야 고장예측에 대한 신뢰도가 높아진다. In the present invention, when the failure history and information of each component of the power generation equipment is inputted to the database, the power generation equipment is divided into large parts such as a turbine, a generator, and a boiler as shown in FIG. 4 (401), and the target equipment (402). After the selection, the part code 404 and the part name 405 are inputted, and the year of completion (406), year of occurrence of failure (407), type of failure (408), and cause of failure (409) of the power generation equipment employing the corresponding parts are entered. ) And measures (410) in sequence, but the correct input increases the reliability of failure prediction.

본 발명에 따른 고장예측 방법은 베이지안 통계기법을 기초로 하는 유지보수의 시변 고장율(TVFR : Time-Varying Failure Rate)과 유지보수의 횟수 등을 고려하여 효과적인 유지보수 주기를 결정하고, 리스크 기준 보수관리(Risk Based Maintenance)등의 기준을 설정하여 발전설비를 효율적으로 유지 관리하여 경제적 손실을 최소화하고 발전설비의 신뢰성을 향상시킬 수 있도록 구성되어 있다.The failure prediction method according to the present invention determines an effective maintenance interval in consideration of the time-varying failure rate (TVFR) of maintenance based on Bayesian statistical technique and the number of maintenance, and risk-based maintenance management. It is designed to minimize the economic loss and improve the reliability of power generation facilities by setting the standards such as (Risk Based Maintenance).

본 발명에 사용된 고장을 예측하는 알고리즘은 베이지안 통계기법 중 하나인 기기 고장 예측에 사용되는 웨이블(Weibull) 분포함수{f(t), 식(1)}를 이용하였다.The algorithm for predicting failure used in the present invention uses a Weibull distribution function {f (t), equation (1)}, which is one of the Bayesian statistical techniques, for predicting device failure.

Figure 112009064522811-PAT00001
Figure 112009064522811-PAT00001

여기에서 α 는 척도모수이고, β 는 형상모수이다.Where α is the scale parameter and β is the shape parameter.

각 부품의 고장 데이터베이스를 이용하여 상기 척도모수 α와 형상모수 β를 다음과 같은 단계를 통하여 구할 수 있다.Using the failure database of each component, the scale parameter α and the shape parameter β can be obtained through the following steps.

a. (고장년도 - 준공년도)를 계산하고, 같은 년도에 발생한 고장건수를 계산하는 단계를 거친다.a. (Failure Year-Completion Year) is calculated, and the number of failures occurring in the same year is calculated.

b. 년도 별 계산된 고장건수(n)를 이용하여 랭크(Rank)를 부여하는 단계를 거친다. 상기 랭크(Rank)는 한해에 발생한 고장건수(n)가 n=1이면 1이고, 고장건수(n)가 1보다 크면 n/2가 된다.b. Rank is given by using the number of failures (n) calculated for each year. The rank is 1 when the number of failures n occurred per year is n = 1, and n / 2 when the number of failures n is greater than one.

c. 상기 랭크(Rank)를 이용하여 역 랭크(Reverse Rank)를 계산하는 단계를 거친다.c. The reverse rank is calculated using the rank.

상기 역 랭크(Reverse Rank)는 랭크(Rank)의 역수이다.  The reverse rank is an inverse of rank.

d. 역 랭크(Reverse Rank)를 이용하여 조절랭크(Adjusted Rank)를 계산하는 단계를 거치며, 보다 구체적으로 조절랭크는 역 랭크를 이전조절랭크로 나누고, 이 값에 고장건수에 1을 더한 후, 이 값을 역 랭크에 1을 더한 값으로 나누어 계산한다. d. The adjusted rank is calculated by using reverse rank. More specifically, the adjusted rank divides the reverse rank by the previous adjustment rank, adds 1 to the number of failures, and then Calculate by dividing by the reverse rank plus one.

수식으로 표현하면, 상기 조절랭크(Adjusted Rank)=[(Reverse Rank)(Previous Adjusted Rank)+(n+1)]/(Reverse Rank+1)Expressed as a formula, the adjusted rank (Adjusted Rank) = [(Reverse Rank) (Previous Adjusted Rank) + (n + 1)] / (Reverse Rank + 1)

이다.to be.

상기 이전 조절랭크(Previous Adjusted Rank)의 초기치는 0이다.The initial value of the previous adjusted rank is zero.

e. 조절랭크(Adjusted Rank)를 이용하여 중앙랭크(Median Rank)를 계산하는 단계를 거치며, 보다 구체적으로 중앙랭크는 상기 조절랭크에서 0.3을 뺀 값을 역 랭크에 0.4를 더한 값으로 나누어 계산한다. 수식으로 표현하면, e. The median rank is calculated using the adjusted rank, and more specifically, the median rank is calculated by dividing the value obtained by subtracting 0.3 from the control rank by the value of 0.4 added to the reverse rank. In formulas,

중앙랭크(Median Rank)=[(Adjusted Rank-0.3)]/[Reverse Rank+0.4]      Median Rank = [(Adjusted Rank-0.3)] / [Reverse Rank + 0.4]

이다.to be.

f. 고장년도에서 준공년도를 뺀 값에 로그(ln)를 취하여 계산하는 단계를 거친다.f. Logarithm (ln) is calculated by subtracting the year of completion from the year of failure.

수식으로 표현하면, In formulas,

x = ln(고장년도-준공년도)x = ln (failure year-completion year)

이다.to be.

g. 1에서 중앙랭크를 뺀 후 역수를 취하고, 역수를 취하여 얻은 값에 로그를 취하고, 로그를 취하여 얻은 값에 다시 로그를 취하여 계산하는 단계를 거친다.g. After subtracting the median rank from 1, take the reciprocal, take the logarithm of the value obtained by taking the reciprocal, and take the log to the value obtained by taking the log and then calculate it.

수식으로 표현하면,In formulas,

y = lnㆍln(1/(1-Median Rank))y = ln · ln (1 / (1-Median Rank))

이다.to be.

상기 f 단계에서 계산한 다수의 값을 x축으로 하고, g 단계에서 각각의 x 값에 대응시켜 계산하여 얻은 다수의 값을 y축으로 하여 도5에서와 같이 그래프상에 도시한다. A plurality of values calculated in step f are represented by the x-axis, and a plurality of values obtained by calculating corresponding to each x value in step g are represented by the y-axis as shown in FIG.

상기 그래프의 회귀분석을 통하여 y = αx + β 의 1차 방정식(도5의 직선)을 구할 수 있으며, 이때 일차방정식의 기울기가 척도모수(α)가 된다.Through the regression analysis of the graph, it is possible to obtain a linear equation of y = αx + β (a straight line in FIG. 5), where the slope of the linear equation becomes the scale parameter α.

또한 웨이블(Weibull) 분포함수의 두 번째 특성인자인 형상모수(β)는 수식(2)로부터 계산하여 얻을 수 있다.In addition, the shape parameter β, which is the second characteristic factor of the Weibull distribution function, can be obtained by calculating from Equation (2).

Figure 112009064522811-PAT00002
Figure 112009064522811-PAT00002

상기 수식(2)에서 형상모수(β)는 다수의 x 값의 평균값에서 다수의 y 값의 평균값을 α 값으로 나눈 값을 뺀 값으로 지수함수를 취하여 계산하여 얻어진다. In Equation (2), the shape parameter β is obtained by calculating the exponential function by subtracting the average value of the plurality of y values from the average value of the plurality of y values divided by the α value.

본 발명은 발전설비 각 구성품을 통계적 방법인 웨이블(Weibull) 분포함수를 이용하여 사용년도 별로 고장예측을 수행할 수 있도록 설계된 알고리즘을 이용하므로 현재까지는 발전설비 각 구성품의 고장을 정확하게 예측할 수 없기 때문에 일정주기가 되면 구성품들을 일방적으로 교체함에 따른 경제적 손실이 크게 있었으나, 이 방법을 통하여 고장 예측을 수행할 수 있어 불필요한 구성품의 교체를 예방할 수 있다. Since the present invention uses an algorithm designed to perform failure prediction for each use year by using the Weibull distribution function, which is a statistical method, it is impossible to accurately predict the failure of each power generation component until now. Although there was a great economic loss from unilateral replacement of components at a certain period, failure prediction can be performed through this method, which prevents unnecessary replacement of components.

본 발명의 현장 적용을 통해 지속적으로 발전설비 각 부품의 고장 예측에 적용 사용할 수 있고, 통계를 이용한 발전설비 구성품의 고장예측 방법은 지속적인 고장 데이터베이스의 업데이트를 통하여 신뢰성 높은 데이터를 확보할 수 있으며 지속적으로 발전설비의 경제적인 운영에 크게 기여할 수 있다.Through on-site application of the present invention can be used continuously to predict failure of each component of the power plant, the failure prediction method of the power plant components using statistics can secure reliable data through the continuous failure database update and continuously It can greatly contribute to the economic operation of power generation facilities.

<산업상 이용가능성>Industrial Applicability

본 발명은 TBM 방식의 발전설비의 정비를 리스크 기준 보수 관리 등의 기준으로 발전설비를 유지 관리하기 위한 것이며, 발전설비의 부품들의 고장예측을 위하여 각 부품들의 고장이력을 데이터베이스화하고 이를 이용하여 통계기법 중의 하나인 베이지안 통계방법을 이용하여 발전설비의 구성 부품들의 고장확률을 예측하는 베이지안 통계기법을 이용한 발전설비 고장예측 방법을 제공하므로 부품의 잦은 교체에 따른 경제적 손실을 줄이고, 교체 시에 발전 중단에 따른 경제적인 손실을 줄일 수 있으므로 산업상 이용가능성이 매우 높다. The present invention is to maintain the power generation facilities on the basis of the risk-based maintenance management, such as maintenance of the TBM-type power generation facilities, the database of failure history of each part to predict the failure of the parts of the power generation facility statistics Bayesian statistical method, which is one of the techniques, provides a method for predicting the failure of power plant equipment using Bayesian statistical technique that predicts the probability of failure of the components of power plant. Industrial applicability is very high as it can reduce the economic loss.

도1 : 고장 데이터베이스를 구축하는 방법과 웨이블 분포함수를 이용하여 고장예측을 계산하는 흐름도1 is a flowchart for calculating a failure prediction using a method of building a failure database and a wavelet distribution function.

도2 : 고장예측 계산을 수행하기 위한 일련의 단계별 흐름도Figure 2: A series of step-by-step flowcharts for performing fault prediction calculations

도3 : 고장예측을 위한 기본 알고리즘Figure 3: Basic algorithm for fault prediction

도4 : 데이터베이스의 구성도의 하나의 예를 나타낸 화면을 캡쳐한 도면4: Screen capture showing an example of the configuration diagram of a database

도5 : 회귀분석을 통한 웨이블 인자 도출Figure 5: Derivation of wavelet factor through regression analysis

Claims (6)

베이지안 통계기법을 이용한 발전설비 고장예측 방법에 있어서,In the failure prediction method of power generation equipment using Bayesian statistical technique, 발전설비 각 부품의 고장 및 정비이력을 조사 분석하여 부품별 고장이력에 대한 정보를 찾아내어 데이터베이스에 입력하는 단계;Investigating and analyzing failure and maintenance history of each component of power generation equipment, finding information on failure history for each component, and inputting the same into a database; 상기 데이터베이스에 입력된 각 부품별 고장이력에 대한 정보를 바탕으로 발전설비 부품의 고장원인 별 각 고장년도와 고장건수를 분석하는 단계;Analyzing each failure year and the number of failures for each failure cause of power generation equipment parts based on the information on the failure history of each component inputted in the database; 상기 고장년도와 고장건수의 분석된 결과를 이용하여 랭크 및 역 랭크를 계산하는 단계;Calculating a rank and an inverse rank using the analyzed result of the failure year and the number of failures; 계산된 역 랭크를 통하여 조절 랭크와 중앙 랭크를 계산하는 단계;Calculating an adjustment rank and a median rank based on the calculated inverse rank; 계산된 중앙랭크와 고정년도를 이용하여 회귀분석 수행하여 웨이블 분포함수의 두 인자인 척도모수(α)와 형상모수(β) 값을 연산하는 단계; 및Performing regression analysis using the calculated median rank and the fixed year to calculate the values of the scale parameter α and the shape parameter β, which are two factors of the wavelet distribution function; And 연산된 척도모수와 형상모수를 웨이블 분포함수에 대입하여 각 부품에 대한 고장예측을 수행하는 단계로 구성된 발전설비 고장예측 방법.A method for predicting a failure of a power generation facility, comprising: performing a failure prediction for each part by substituting the calculated scale parameter and the shape parameter into a wavelet distribution function. 청구항1에 있어서, The method according to claim 1, 상기 데이터베이스에 입력된 각 부품에 대한 고장 및 정비이력은 부품이 채용된 발전설비의 준공년도, 부품의 고장발생년도, 고장종류, 고장원인 및 조치사항으로 구성됨을 특징으로 하는 발전설비 고장예측 방법.Breakdown and maintenance history for each part entered into the database is a generation equipment failure prediction method, characterized in that consisting of the year of completion of the power generation equipment employing the parts, the year of occurrence of the failure, the type of failure, the cause of the failure and measures. 청구항1에 있어서, The method according to claim 1, 상기 척도모수(α)는 고장년도에서 준공년도를 뺀 값에 로그(ln)를 취한 값으로 x 축 값을 구하고, The scale parameter (α) is obtained by subtracting the log year (ln) from the failure year minus the completion year and obtaining the x-axis value. 1에서 중앙 랭크를 뺀 후 역수를 취하고, 역수를 취하여 얻은 값에 로그를 취하고, 로그를 취하여 얻은 값에 다시 로그를 취하여 y 값을 구하며,Subtract the central rank from 1, take the inverse, take the logarithm of the value obtained by taking the inverse, take the logarithm of the value obtained by taking the logarithm, and find the y value, 상기 다수의 x 값 및 x 에 대응하는 다수의 y 값으로부터 구한 일차함수의 기울기를 연산하여 획득함을 특징으로 하는 발전설비 고장예측 방법. And calculating the slope of the first function obtained from the plurality of x values and the plurality of y values corresponding to x. 청구항3에 있어서, The method according to claim 3, 상기 형상모수(β)는 y 의 평균값을 α 값으로 나누고, 이 값에서 x 의 평균값에서 뺀 값에 지수함수 취하여 획득함을 특징으로 하는 발전설비 고장예측 방법. The shape parameter β is obtained by dividing an average value of y by an α value and taking an exponential function on a value obtained by subtracting an average value of x from this value. 청구항1내지 청구항4 중 어느 한 항에 있어서, The method according to any one of claims 1 to 4, 상기 랭크는 한해에 발생한 고장건수(n)가 n=1이면 1이며, 고장건수(n)가 1보다 크면 n/2로 하여 계산함을 특징으로 하는 발전설비 고장예측 방법.The rank is 1 when the number of failures (n) occurred in a year is n = 1, and when the number of failures (n) is greater than 1, it is calculated as n / 2. 청구항1내지 청구항4 중 어느 한 항에 있어서, The method according to any one of claims 1 to 4, 상기 조절랭크는 역 랭크를 이전조절랭크로 나누고, 나눈 값에 고장건수에 1을 더한 후, 이 값을 역 랭크에 1을 더한 값으로 나누어 계산함을 특징으로 하는 발전설비 고장예측 방법.The control rank is calculated by dividing the reverse rank by the previous control rank, add 1 to the number of failures to the divided value, and calculates this value by dividing the value by adding 1 to the reverse rank.
KR1020090100320A 2009-10-21 2009-10-21 A Method of Failure Prediction at Power Plant Using Bayesian Statics KR101086027B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020090100320A KR101086027B1 (en) 2009-10-21 2009-10-21 A Method of Failure Prediction at Power Plant Using Bayesian Statics

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020090100320A KR101086027B1 (en) 2009-10-21 2009-10-21 A Method of Failure Prediction at Power Plant Using Bayesian Statics

Publications (2)

Publication Number Publication Date
KR20110043277A true KR20110043277A (en) 2011-04-27
KR101086027B1 KR101086027B1 (en) 2011-11-22

Family

ID=44048576

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020090100320A KR101086027B1 (en) 2009-10-21 2009-10-21 A Method of Failure Prediction at Power Plant Using Bayesian Statics

Country Status (1)

Country Link
KR (1) KR101086027B1 (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103198219A (en) * 2013-04-03 2013-07-10 中国南方电网有限责任公司超高压输电公司广州局 Time varying failure rate estimation method of relay protection device and based on failure mode distinction and norm distribution
KR101369341B1 (en) * 2012-05-03 2014-03-06 한국수력원자력 주식회사 System for constructing integrated management database of electronic circuit board
KR20160026811A (en) 2014-09-01 2016-03-09 경희대학교 산학협력단 Apparatus and method for goal attainment diagnostic
CN110110967A (en) * 2019-04-09 2019-08-09 华能山东石岛湾核电有限公司 A kind of risk evaluating method suitable for the change of specific power plant's licensing benchmark
KR20220150696A (en) * 2021-05-04 2022-11-11 한국남동발전 주식회사 Failure prediction system based on failure prediction model
KR20230032799A (en) * 2021-08-31 2023-03-07 주식회사 나루 Leisure boat sharing safety prediction management system

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5819251B2 (en) 2012-05-14 2015-11-18 住友重機械工業株式会社 Operation control system for circulating fluidized bed boiler.
JP5925587B2 (en) 2012-05-14 2016-05-25 住友重機械工業株式会社 Operation diagnosis method and operation diagnosis device for circulating fluidized bed boiler
KR101869884B1 (en) 2016-10-06 2018-06-21 더블유에이테크놀러지 주식회사 Server for managing generating equipment
US11703846B2 (en) 2019-05-30 2023-07-18 Cameron International Corporation Equipment failure diagnostics using Bayesian inference
KR102110601B1 (en) 2019-11-27 2020-05-13 (주)대한기전 Facility failure prediction system

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101369341B1 (en) * 2012-05-03 2014-03-06 한국수력원자력 주식회사 System for constructing integrated management database of electronic circuit board
CN103198219A (en) * 2013-04-03 2013-07-10 中国南方电网有限责任公司超高压输电公司广州局 Time varying failure rate estimation method of relay protection device and based on failure mode distinction and norm distribution
KR20160026811A (en) 2014-09-01 2016-03-09 경희대학교 산학협력단 Apparatus and method for goal attainment diagnostic
US10417573B2 (en) 2014-09-01 2019-09-17 University-Industry Cooperation Foundation Of Kyung Hee University Goal-attainment assessment apparatus and method
CN110110967A (en) * 2019-04-09 2019-08-09 华能山东石岛湾核电有限公司 A kind of risk evaluating method suitable for the change of specific power plant's licensing benchmark
CN110110967B (en) * 2019-04-09 2022-12-20 华能山东石岛湾核电有限公司 Risk evaluation method suitable for nuclear power station license benchmark change
KR20220150696A (en) * 2021-05-04 2022-11-11 한국남동발전 주식회사 Failure prediction system based on failure prediction model
KR20230032799A (en) * 2021-08-31 2023-03-07 주식회사 나루 Leisure boat sharing safety prediction management system

Also Published As

Publication number Publication date
KR101086027B1 (en) 2011-11-22

Similar Documents

Publication Publication Date Title
KR101086027B1 (en) A Method of Failure Prediction at Power Plant Using Bayesian Statics
Hu et al. A prediction method for the real-time remaining useful life of wind turbine bearings based on the Wiener process
US10521490B2 (en) Equipment maintenance management system and equipment maintenance management method
AU2016201724B2 (en) Method and system for predicting equipment failure
CN107423414B (en) Information transfer model-based process industry complex electromechanical system fault tracing method
CN105844050A (en) Time correlation based replacement time calculation method for numerical control machine tool components
Bangalore et al. An approach for self evolving neural network based algorithm for fault prognosis in wind turbine
Blancke et al. A holistic multi-failure mode prognosis approach for complex equipment
JP7053152B2 (en) Systems and methods for optimizing recommended inspection intervals
CN111461450A (en) Hydroelectric generating set degradation degree prediction method based on EEMD and L STM
Zhang et al. Age-dependent maintenance strategies of medium-voltage circuit-breakers and transformers
JP2017151980A5 (en)
KR102366922B1 (en) heavy electric equipment deterioration prediction system based on time-series isolation diagnostic information and method therefor
Nichenametla et al. Optimizing life cycle cost of wind turbine blades using predictive analytics in effective maintenance planning
Ahmed et al. Throughput-based importance measures of multistate production systems
KR20140068521A (en) System and method for equipment stability and engineering process effectiveness assessment
Asjad et al. A life cycle cost based approach of O&M support for mechanical systems
Pariaman et al. AVAILABILITY ANALYSIS OF THE INTEGRATED MAINTENANCE TECHNIQUE BASED ON RELIABILITY, RISK, AND CONDITION IN POWER PLANTS.
RU2703874C1 (en) Method of monitoring and predicting operation of a gas turbine plant using a matrix of defects
Nielsen et al. Risk-based decision making for deterioration processes using POMDP
Friederich et al. A Framework for Validating Data-Driven Discrete-Event Simulation Models of Cyber-Physical Production Systems
CN117897712A (en) Dimension reduction on a wafer
Sadiki et al. Impact of intelligent wireless sensor network on predictive maintenance cost
Beduschi et al. Optimizing rotating equipment maintenance through machine learning algorithm
CN113591266A (en) Method and system for analyzing fault probability of electric energy meter

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
AMND Amendment
E601 Decision to refuse application
AMND Amendment
X701 Decision to grant (after re-examination)
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20141114

Year of fee payment: 4

FPAY Annual fee payment

Payment date: 20151117

Year of fee payment: 5

FPAY Annual fee payment

Payment date: 20161115

Year of fee payment: 6

FPAY Annual fee payment

Payment date: 20181101

Year of fee payment: 8