KR101896115B1 - Method and apparatus for estimating corrosion rate of oil plant facility - Google Patents
Method and apparatus for estimating corrosion rate of oil plant facility Download PDFInfo
- Publication number
- KR101896115B1 KR101896115B1 KR1020160142752A KR20160142752A KR101896115B1 KR 101896115 B1 KR101896115 B1 KR 101896115B1 KR 1020160142752 A KR1020160142752 A KR 1020160142752A KR 20160142752 A KR20160142752 A KR 20160142752A KR 101896115 B1 KR101896115 B1 KR 101896115B1
- Authority
- KR
- South Korea
- Prior art keywords
- corrosion
- data
- rate
- corrosion rate
- predicting
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N17/00—Investigating resistance of materials to the weather, to corrosion, or to light
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/20—Metals
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/20—Metals
- G01N33/204—Structure thereof, e.g. crystal structure
- G01N33/2045—Defects
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/66—Subsurface modeling
- G01V2210/665—Subsurface modeling using geostatistical modeling
- G01V2210/6652—Kriging
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Ecology (AREA)
- Environmental & Geological Engineering (AREA)
- Biodiversity & Conservation Biology (AREA)
- Environmental Sciences (AREA)
- Engineering & Computer Science (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Crystallography & Structural Chemistry (AREA)
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
Abstract
유류 설비의 부식 속도 예측 방법으로서, 상기 유류 설비의 부식 데이터를 획득하는 단계; 획득된 부식 데이터로부터 이상점을 제거하는 단계; 및 이상점이 제거된 부식 데이터를 크리깅 모델(Kriging model)에 적용하여 부식 속도를 예측하는 단계를 포함한다. 상기 이상점을 제거하는 단계는, 부식 속도와 부식 변수들의 측정 시간을 고려하여 상기 부식 데이터에 구분적 선형화(piecewise linearization) 및 시간 가중 평균화 기법을 적용하는 단계; 및 적용된 결과에 Q-test를 수행하여 소정 임계를 벗어나는 부식 데이터를 제거하는 단계를 포함하는 방법이 제공된다. 본 발명에 의하여, 정유 공장에서 부식 속도를 예측할 수 있으므로, 부식 예측 결과를 토대로 부식 평가법 구축, 신규 화학 공정 설비 재질 선정, 부식 관리, 설비 진단 평가 일정 계획 수립 등에 활용할 수 있어서 설비의 안전성 향상을 도모할 수 있다.A method of predicting erosion rate of an oil plant, comprising: obtaining corrosion data of the oil plant; Removing an abnormal point from the obtained corrosion data; And applying the corrosion data from which the anomalous points have been removed to a kriging model to predict the corrosion rate. The step of removing the anomaly point includes applying piecewise linearization and time weighted averaging to the corrosion data in consideration of the corrosion rate and the measurement time of the corrosion parameters; And performing a Q-test on the applied results to remove erosion data that deviates from a predetermined threshold. According to the present invention, since the corrosion rate can be predicted in the refinery, it is possible to utilize the corrosion evaluation method, the material selection of the new chemical processing facility, the corrosion management, and the schedule of the facility diagnosis evaluation based on the corrosion prediction result. can do.
Description
본 발명은 정유 플랜트에서 사용되는 유류 설비의 부식 속도를 예측하기 위한 것으로서, 특히 크리깅 모델(Kriging model)을 사용하여 용이하면서도 정확하게 유류 설비의 부식 속도를 예측하고, 예측된 결과를 다양한 기법들과 비교하여 검증하기 위한 기술에 관한 것이다.The present invention is intended to predict the corrosion rate of oil installations used in refinery plants. In particular, the kriging model can be used to easily and accurately predict the corrosion rate of oil installations and to compare the predicted results with various techniques And to a technique for verifying the same.
정유 공장에 사용되는 유류 설비에는 다양한 종류가 사용되는데, 이러한 설비가 부식되면 정유 공장의 가동에 큰 제한이 된다. 특히, 증류탑의 상부는 탄소강 소재로 제작되어 부식에 취약하며, 장치 산업의 특성상 수년 간 가동 정지(shut-down) 없이 운전하기 때문에 증류탑 상부에 부식이 발생되면 정유 공장의 안전에도 치명적인 영향을 미치게 된다. 따라서, 이러한 설비의 부식 정도를 사전에 예측하여, 사고가 발생하기 이전에 조치하는 것이 필요하다.Various types of oil facilities used in oil refineries are used. When these facilities are corroded, the operation of oil refineries is very limited. Particularly, since the upper part of the distillation tower is made of carbon steel and is vulnerable to corrosion and corrosion of the upper part of the distillation tower is caused by the operation of the apparatus without shut-down for several years due to the characteristics of the equipment industry, . Therefore, it is necessary to anticipate the extent of corrosion of these facilities and take action before the accident occurs.
그러나, 부식은 매우 느린 반응인 것은 물론, 그 메커니즘이 복잡하여 이론적 모델을 실제 정유 공장에 적용하기에는 많은 어려움이 있다.However, corrosion is a very slow reaction, and the mechanism is complicated, which makes it difficult to apply the theoretical model to a real refinery.
부식 반응에 대한 기존 모델들에는 전기 화학 모델, 물질 전달 모델, 경험 모델들이 있다. 이 중 전기 화학 모델은 철 표면에서의 부식 반응에 대한 반응 속도에 중점을 둔 모델이다. 또한, 물질 전달 모델은 철 표면에 생성된 산화철로 된 얇은 층이 부식 유발 물질이 철 표면으로 확산되는 것을 막는 것에 중점을 둔다. 즉, 물질 전달 모델은 물질 전달이 부식 속도를 결정하는 것이라고 가정한 모델이다. 또한, 경험적 모델은 실험 결과를 바탕으로 부식 속도를 예측하는 모델이다. 특히, 경험적 모델에서는 시편이나 전류 측정 등을 통하여 간접적으로 부식 속도 또는 부식 정도를 예측하고, 실험 데이터 등을 통계적으로 처리하여 선형이나 다항식, 이론식 등으로 회귀 분석을 하여 부식을 예측하기도 한다. Existing models for corrosion reactions include electrochemical models, mass transfer models, and empirical models. Among these, the electrochemical model focuses on the reaction rate to the corrosion reaction on the iron surface. In addition, the mass transfer model focuses on preventing the thin layer of iron oxide produced on the iron surface from spreading corrosion-causing substances to the iron surface. That is, the mass transfer model assumes that mass transfer determines the rate of corrosion. The empirical model is a model for predicting corrosion rate based on experimental results. Especially, in the empirical model, the corrosion rate or the degree of corrosion is indirectly predicted through the specimen or current measurement, and the corrosion is predicted by regression analysis of the linearity, polynomial, and theoretical expression by statistically processing the experimental data.
하지만, 정유 공장의 경우, 정확한 물질의 조성을 알기가 힘들고 여러 가지 불순물들이 포함되어 있으므로 실험이 어렵고, 이론식에 필요한 모든 값들을 측정하는데 제약이 있다. 또한, 단순한 선형이나 다항식 형태의 회귀분석은 복잡한 부식 메커니즘을 나타내는데 한계가 있다. 즉, 기존 모델들에는 실험적으로 결정해야 하는 많은 파라미터들이 포함되기 때문에 그 결과를 예측하기가 어렵다. 또한, 부식 반응은 매우 느린 반응이므로, 실험을 하더라도 며칠에서 최대 몇 주 정도 밖에 수행하지 못하는 경우가 많아, 몇 년 동안 가동되는 정유 공장에 적용하기가 쉽지 않다. 더 나아가, 정유 공정에 사용되는 원유 등은 정확한 조성을 알 수 없는 복합한 혼합물인 데다가, 똑같은 조건으로 실험 조건을 설정하는 것이 곤란하다. 그리고, H2S가 극소량이 첨가되어도 CO2만 존재할 때에 비해 90%나 부식 속도가 감소할 정도로 부식 반응의 메커니즘은 매우 복잡하기 때문에, 이를 이론적으로 예측하기란 더욱 어렵다.However, refining plants have difficulties in knowing the exact composition of the substance and contain various impurities, so it is difficult to experiment and there is a limitation in measuring all the necessary values in the formula. In addition, simple linear or polynomial regression analysis has limitations in representing complex corrosion mechanisms. In other words, it is difficult to predict the results because existing models contain many parameters that must be experimentally determined. In addition, the corrosion reaction is very slow, so it is not easy to apply to a refinery operating over a few years, since experiments often only take a few days to several weeks. Furthermore, the crude oil used in the refining process is a complex mixture in which the exact composition is unknown, and it is difficult to set the experimental conditions under the same conditions. It is more difficult to theoretically predict the mechanism of the corrosion reaction to such an extent that the corrosion rate is reduced by 90% compared with the case where only H 2 S is added even when only CO 2 is present.
대한민국 등록 특허 번호 제 10-1199462 호 "탄소강 파이프 부식 속도 예측 장치, 그 예측 방법 및 그 기록 매체"(이하 선행기술1)에 따르면, 다양한 부식 시간에 대하여 수집된 특정 금속 구조물의 부식량 데이터를 통계적으로 처리하여 특정 금속 구조물이 시간의 경과에 따라 부식되는 양을 예측할 수 있는 일반화된 실험식이 제공되어 있다. 선행 기술1은, 예측된 부식량에 근거하여 수리, 교체, 폐기 등의 각종 조치를 계획 및 실행할 수 있도록 하는 부식 환경에서의 시간 경과에 따른 부식량 예측 방법에 관한 기술이다. 이를 위하여, 이를 위해 선행 기술에서는 다양한 부식 시간에 다양한 부식 노출 시간에 대하여 수집된 특정 금속 구조물의 부식량 데이터를 입력하고, 입력된 데이터에 대해 평균값은 최대로 하고 COV는 최소값이 되도록 하는 부식량 간격을 결정하며, 결정된 부식량 간격에 대한 각 부식 노출 시간별 부식량 데이터에 대해 형상 파라미터(α) 및 치수 파라미터(β)를 계산하고, 계산된 복수의 부식 노출 시간별 형상 파라미터(α) 및 치수 파라미터(β)로부터 각 파라미터(α,β)를 근사적으로 표현하는 부식 노출 시간에 대한 선형식 또는 다항식을 계산하며, 이를 이용하여 시간종속 부식량 예측 실험식을 완성한다. According to Korean Patent No. 10-1199462 entitled " Apparatus for Predicting Carbon Steel Pipe Corrosion Rate, its Prediction Method and its Recording Medium "(hereinafter referred to as Prior Art 1), corrosion data of a specific metal structure collected for various corrosion times are statistically To provide a generalized empirical formula that can predict the amount of corrosion of a particular metallic structure over time.
그러나, 선행 기술1에 따르면 금속 구조물별 부식량 데이터가 필요하고, 다양한 부식시간 및 부식 노출 시간 별로 상이한 데이터가 필요하므로, 부식량을 측정하는 과정이 복잡하고 실제 적용하기 어려운 단점이 있다.However, according to Prior
따라서, 유류 설비의 부식 속도를 상대적으로 용이하면서도 정확하게 예측하기 위한 기술이 절실히 요구된다.Therefore, there is a desperate need for a technique for relatively easily and accurately predicting the erosion rate of oil installations.
본 발명의 목적은 실제 공장 데이터를 바탕으로 지구 통계학에서 사용되는 크리깅 모델을 사용하여 부식 속도를 예측함으로써, 부식 속도를 용이하고 정확하게 예측할 수 있는 방법을 제공하는 것이다. It is an object of the present invention to provide a method for easily and accurately predicting the corrosion rate by predicting the corrosion rate using a kriging model used in geostatistics based on actual plant data.
또한, 본 발명의 목적은 예측된 결과를 다중 선형 회귀 분석 및 인공 신경망의 결과와 비교, 분석함으로써 예측 결과를 검증할 수 있는 유류 설비의 부식 속도 예측 장치를 제공하는 것이다.
It is also an object of the present invention to provide an apparatus for predicting a corrosion rate of a petroleum plant capable of verifying a predicted result by comparing and analyzing a predicted result with the results of multiple linear regression analysis and artificial neural network.
상기와 같은 목적들을 달성하기 위한 본 발명의 일면은, 유류 설비의 부식 속도 예측 방법으로서, 상기 유류 설비의 부식 데이터를 획득하는 단계; 획득된 부식 데이터로부터 이상점을 제거하는 단계; 및 이상점이 제거된 부식 데이터를 크리깅 모델(Kriging model)에 적용하여 부식 속도를 예측하는 단계를 포함한다. 상기 이상점을 제거하는 단계는, 부식 속도와 부식 변수들의 측정 시간을 고려하여 상기 부식 데이터에 구분적 선형화(piecewise linearization) 및 시간 가중 평균화 기법을 적용하는 단계; 및 적용된 결과에 Q-test를 수행하여 소정 임계를 벗어나는 부식 데이터를 제거하는 단계를 포함한다. 상기 부식 속도를 예측하는 단계는, 상기 크리깅 모델에 적용될 부식 변수들을 획득하는 단계; 획득된 부식 변수와 부식 속도의 상관 관계를 유의 확률로 분석하는 단계; 및 분석 결과에 따라 유의도가 낮은 부식 변수를 제거하여 크리깅 모델의 초평면의 차원을 낮추는 단계를 포함한다. 상기 부식 속도를 예측하는 단계는, 를 사용하여 부식 속도를 예측하는 단계를 더 포함하며, 여기에서, z*는 예측하려는 부식 속도, z i 는 예측하려는 환경과 유사한 환경에서 측정된 부식 데이터, n은 상기 유사한 환경에서 측정된 부식 데이터의 수, λi는 z i 에 적용되는 가중치를 나타낸다. 또한, 본 발명에 의한 유류 설비의 부식 속도 예측 방법은 예측된 결과를 다중 선형 회귀 분석 결과 및 인공 신경망 분석 결과 중 적어도 하나와 함께 제공하는 단계를 더 포함한다. 상기 유류 설비는 정류 공장의 증류탑 상부일 수 있다.According to an aspect of the present invention, there is provided a method for predicting a corrosion rate of an oil facility, the method comprising: acquiring corrosion data of the oil facility; Removing an abnormal point from the obtained corrosion data; And applying the corrosion data from which the anomalous points have been removed to a kriging model to predict the corrosion rate. The step of removing the anomaly point includes applying piecewise linearization and time weighted averaging to the corrosion data in consideration of the corrosion rate and the measurement time of the corrosion parameters; And performing a Q-test on the applied result to remove corrosion data that deviates from a predetermined threshold. The step of predicting the erosion rate comprises: obtaining corrosion parameters to be applied to the kriging model; Analyzing the correlation between the obtained corrosion parameter and the corrosion rate with a significant probability; And lowering the dimension of the hyperplane of the kriging model by removing the low-significance corrosion parameter according to the analysis result. The step of predicting the corrosion rate comprises: Using further comprising the step of predicting the erosion rate, where, z * is corrosive to predict speed, z i is the corrosion data obtained in similar to predict environmental environment, n is the corrosion data obtained in the similar environment And? I represents a weight applied to z i . The method of predicting erosion rate of the oil plant according to the present invention further includes the step of providing predicted results together with at least one of the results of the polynomial regression analysis and the ANN analysis. The oil facility may be above the distillation tower of the rectification plant.
상기와 같은 목적들을 달성하기 위한 본 발명의 다른 면은 유류 설비의 부식 속도 예측 장치로서, 상기 유류 설비의 부식 데이터를 획득하는 부식 데이터 획득부; 획득된 부식 데이터로부터 이상점을 제거하는 이상점 제거부; 및 이상점이 제거된 부식 데이터를 크리깅 모델(Kriging model)에 적용하여 부식 속도를 예측하는 부식 속도 예측부를 포함한다. 상기 이상점 제거부는, 부식 속도와 부식 변수들의 측정 시간을 고려하여 상기 부식 데이터에 구분적 선형화(piecewise linearization) 및 시간 가중 평균화 기법을 적용하고, 적용된 결과에 Q-test를 수행하여 소정 임계를 벗어나는 부식 데이터를 제거하도록 구성된다. 상기 부식 속도 예측부는, 상기 크리깅 모델에 적용될 부식 변수들을 획득하고, 획득된 부식 변수와 부식 속도의 상관 관계를 유의 확률로 분석하며, 분석 결과에 따라 유의도가 낮은 부식 변수를 제거하여 크리깅 모델의 초평면의 차원을 낮추도록 구성된다. 상기 부식 속도 예측부는, 를 사용하여 부식 속도를 예측하고, 여기에서, z*는 예측하려는 부식 속도, z i 는 예측하려는 환경과 유사한 환경에서 측정된 부식 데이터, n은 상기 유사한 환경에서 측정된 부식 데이터의 수, λi는 z i 에 적용되는 가중치를 나타낸다. 본 발명의 유류 설비의 부식 속도 예측 장치는 예측된 결과를 다중 선형 회귀 분석 결과 및 인공 신경망 분석 결과 중 적어도 하나와 함께 제공하는 결과 제공부를 더 포함하며, 상기 유류 설비는 정류 공장의 증류탑 상부일 수 있다.
According to another aspect of the present invention, there is provided an apparatus for predicting a corrosion rate of a oil facility, comprising: a corrosion data obtaining unit obtaining corrosion data of the oil facility; An abnormal point removing unit for removing an abnormal point from the obtained corrosion data; And a corrosion rate predicting unit for predicting a corrosion rate by applying corrosion data from which the anomalous points are removed to a kriging model. The anomaly removing unit applies piecewise linearization and time-weighted averaging to the corrosion data in consideration of the corrosion rate and the measurement time of the corrosion parameters, performs a Q-test on the result of the application, And is configured to remove corrosion data. The corrosion rate predicting unit obtains the corrosion parameters to be applied to the kriging model, analyzes the correlation between the obtained corrosion parameter and the corrosion rate with a significant probability, removes the corrosion parameter having a low degree of significance according to the analysis result, It is configured to lower the dimension of the hyperplane. The corrosion rate predicting unit may include: Where z * is the corrosion rate to be predicted, z i is the corrosion data measured in an environment similar to the environment to be predicted, n is the number of corrosion data measured in the similar environment, and λ i Represents a weight applied to z i . The apparatus for predicting erosion rate of a plant according to the present invention further comprises a result providing unit for providing a predicted result together with at least one of a result of a polynomial regression analysis and an artificial neural network analysis, have.
본 발명에 의하여, 정유 공장에서 부식 속도를 예측할 수 있으므로, 부식 예측 결과를 토대로 부식 평가법 구축, 신규 화학 공정 설비 재질 선정, 부식 관리, 설비 진단 평가 일정 계획 수립 등에 활용할 수 있어서 설비의 안전성 향상을 도모할 수 있다.According to the present invention, since the corrosion rate can be predicted in the refinery, it is possible to utilize the corrosion evaluation method, the material selection of the new chemical processing facility, the corrosion management, and the schedule of the facility diagnosis evaluation based on the corrosion prediction result. can do.
또한, 본 발명에 의하여, 정유 공장에서의 부식 속도를 예측하기 위하여 고등 통계 기법인 크리깅 모델을 이용하기 때문에 개별 정유 공장에 적합한 부식 예측이 가능해진다.In addition, according to the present invention, since the kriging model, which is a statistical technique, is used to predict the corrosion rate in an oil refinery, corrosion prediction suitable for an individual refinery can be made.
더 나아가, 본 발명에 의하여, 다중 선형 회귀 분석과 같은 일반적인 통계 기법과 본 발명의 예측 결과를 비교함으로써 보다 정확한 결과를 얻을 수 있고, 궁극적으로는 정유 공장의 안전한 운전과 더욱 긴 시간동안 연속 운전을 가능케 하여 수익을 극대화할 수 있다.
Further, according to the present invention, it is possible to obtain more accurate results by comparing general statistical methods such as multiple linear regression analysis with the prediction results of the present invention, and ultimately, to operate safely and continuously for a longer time Thereby maximizing profits.
도 1 은 본 발명의 일 실시예에 의한 유류 설비의 부식 속도 예측 방법을 개략적으로 나타내는 흐름도이다.
도 2 는 본 발명에서 사용되는 구분적 선형화를 나타내는 그래프이다.
도 3 은 본 발명에서 사용되는 크리깅 모델을 설명하는 그래프이다.
도 4 는 부식 변수가 6 개인 경우의 실험적 베리오그램과 이론적 베리오그램을 나타내는 그래프이다.
도 5 는 부식 변수가 6 개인 경우의 예측 결과들을 나타내는 그래프이다.
도 6 은 부식 변수가 4 개인 경우의 실험적 베리오그램과 이론적 베리오그램을 나타내는 그래프이다.
도 7 은 부식 변수가 4 개인 경우의 예측 결과들을 나타내는 그래프이다.
도 8 은 본 발명의 다른 실시예에 의한 유류 설비의 부식 속도 예측 장치를 개략적으로 나타내는 블록도이다.FIG. 1 is a flowchart schematically illustrating a corrosion rate predicting method of an oil plant according to an embodiment of the present invention.
FIG. 2 is a graph showing the segmental linearization used in the present invention. FIG.
3 is a graph illustrating a kriging model used in the present invention.
4 is a graph showing an experimental variogram and a theoretical variogram in the case of six corrosion variables.
FIG. 5 is a graph showing the prediction results in the case where the corrosion parameter is six.
6 is a graph showing an experimental variogram and a theoretical variogram in the case of four corrosion variables.
FIG. 7 is a graph showing prediction results in the case where four corrosion variables are used.
8 is a block diagram schematically showing an apparatus for predicting corrosion rate of a oil plant according to another embodiment of the present invention.
본 발명과 본 발명의 동작상의 이점 및 본 발명의 실시에 의하여 달성되는 목적을 충분히 이해하기 위해서는 본 발명의 바람직한 실시예를 예시하는 첨부 도면 및 첨부 도면에 기재된 내용을 참조하여야만 한다. In order to fully understand the present invention, operational advantages of the present invention, and objects achieved by the practice of the present invention, reference should be made to the accompanying drawings and the accompanying drawings which illustrate preferred embodiments of the present invention.
이하, 첨부한 도면을 참조하여 본 발명의 바람직한 실시예를 설명함으로서, 본 발명을 상세히 설명한다. 그러나, 본 발명은 여러 가지 상이한 형태로 구현될 수 있으며, 설명하는 실시예에 한정되는 것이 아니다. 그리고, 본 발명을 명확하게 설명하기 위하여 설명과 관계없는 부분은 생략되며, 도면의 동일한 참조부호는 동일한 부재임을 나타낸다. Hereinafter, the present invention will be described in detail with reference to the preferred embodiments of the present invention with reference to the accompanying drawings. However, the present invention can be implemented in various different forms, and is not limited to the embodiments described. In order to clearly describe the present invention, parts that are not related to the description are omitted, and the same reference numerals in the drawings denote the same members.
도 1 은 본 발명의 일 실시예에 의한 유류 설비의 부식 속도 예측 방법을 개략적으로 나타내는 흐름도이다.FIG. 1 is a flowchart schematically illustrating a corrosion rate predicting method of an oil plant according to an embodiment of the present invention.
우선, 정유 공장에서의 부식 예측을 위해서 정유 공장에서 측정 가능한 부식 인자들과 부식 속도에 대한 과거의 운전 데이터들을 입력받는다(S110).First, in order to predict the corrosion at the refinery, the past operation data on corrosion factors measurable in the refinery and the corrosion rate are inputted (S110).
정유 공장에서 측정가능한 부식 인자들의 예로는 pH, Cl-, Fe2+, H2S 농도, NH3 농도, 유량 등이 있다. CO2와 H2S에 의한 부식은 정유 공장에서의 주요한 부식 반응으로 주요 부식인자이다. 그 외에도 다양한 부식에 영향을 미치는 인자들이 있으나 정확한 부식 메커니즘이 밝혀지지 않아 이론적 모델의 적용에 한계가 있다. 또한, 실제 정유 공장에서 측정할 수 있는 변수들의 수와 종류에 제약이 따르며 정유 공장마다 운전 환경이 달라 각 정유 공장에 맞는 부식 예측 모델이 필요하다. 그러므로, 본 발명에서는 실제 측정 데이터를 사용하여 부식 속도를 예측한다.Examples of corrosive factors measurable in refineries are pH, Cl - , Fe 2+ , H 2 S concentration, NH 3 concentration, and flow rate. Corrosion by CO2 and H 2 S is a major corrosion factor in the refinery due to major corrosion reactions. There are other factors affecting various corrosion, but the precise erosion mechanism is not clarified, so the application of the theoretical model is limited. In addition, there are restrictions on the number and types of variables that can be measured in a real refinery. Different refinery plants have different operating environments and require a corrosion prediction model for each refinery. Therefore, in the present invention, corrosion rate is predicted using actual measurement data.
데이터가 획득되면, 획득된 부식 데이터로부터 이상점을 제거한다(S130). Once the data is obtained, the abnormal point is removed from the obtained corrosion data (S130).
이상점을 제거하기 위해서는 Q-test를 수행해야 하며, 본 발명에서는 Q-test를 수행하기 전에 구분적 선형화(piecewise linearization) 및 시간 가중 평균화 기법을 적용할 수 있다.In order to remove the anomalous points, a Q-test must be performed. In the present invention, piecewise linearization and time-weighted averaging techniques can be applied before performing the Q-test.
이와 같이 이상점이 제거되면, 나머지 데이터를 적용할 부식 변수를 획득한다(S140). 본 발명에서는 부식 변수들로서, H2S, Cl-, Fe2+, NH3, pH, 유량의 6 개의 변수를 사용하였으나, 본 발명은 이에 한정되는 것이 아니고 부식 속도에 영향을 미칠 수 있는 다양한 부식 변수들을 도입하는 것이 가능함은 물론이다.When the abnormal point is removed as described above, a corrosion parameter to which the remaining data is to be applied is obtained (S140). In the present invention, six parameters of H 2 S, Cl - , Fe 2+ , NH 3 , pH, and flow rate were used as the corrosion parameters. However, the present invention is not limited thereto and various corrosions It is of course possible to introduce variables.
부식 변수가 결정되면, 획득된 부식 변수와 부식 속도의 상관 관계를 유의 확률로 분석하고 유의도가 낮은 부식 변수를 제거한다(S160). 이와 같이 부식 변수를 제거하면 크리깅 모델의 초평면의 차원을 낮추는 효과를 얻게 된다.Once the corrosion parameters are determined, the correlation between the obtained corrosion parameters and the corrosion rate is analyzed with a significant probability, and the corrosion parameters with low significance are removed (S160). Removal of the corrosion parameters in this way can reduce the dimension of the hyperplane of the kriging model.
초평면의 차원이 결정되면, 이상점이 제거된 부식 데이터를 크리깅 모델(Kriging model)에 적용하여 부식 속도를 예측한다(S170). 이와 같이, 본 발명에서는 고등 통계 기법인 크리깅 모델을 사용하여 정유 공장의 운전 데이터를 입력받아 해당 정유 공장의 운전 환경에 맞는 부식 속도를 예측한다. 크리깅 모델은 일반적으로 공간 좌표계에서 사용되지만, 본 발명에서는 이러한 모델을 부식 인자들로 이루어진 비공간 좌표계의 초평면에 적용하여 부식 예측을 수행한다.If the dimension of the hyperplane is determined, the corrosion rate is estimated by applying the corrosion data from which the abnormal point is removed to the kriging model (S170). As described above, in the present invention, the operation data of the refinery is input using the Kriging model, which is an advanced statistical method, and the corrosion rate suitable for the operation environment of the refinery is predicted. The kriging model is generally used in a spatial coordinate system, but in the present invention, this model is applied to a hyperplane of a non-spatial coordinate system made up of corrosion factors to perform corrosion prediction.
부식 속도가 예측되면, 다중 선형 회귀 분석 및 인공 신경망 분석을 수행하고, 수행된 결과를 크리깅 모델을 사용하여 예측된 결과와 함께 사용자에게 제공한다(S190).When the corrosion rate is predicted, a polynomial linear regression analysis and an artificial neural network analysis are performed, and the result is provided to the user along with the predicted result using the kriging model (S190).
즉, 본 발명에서는 비공간 좌표인 부식 변수들(H2S, Cl-, Fe2+, NH3, pH, 유량)로 이루어진 초평면을 구성하고, 유의도에 따라 부식 변수의 개수를 감축함으로써 연산 부담을 줄이는 한편, 종래 기술보다 용이하게 부식 속도를 예측한다.That is, in the present invention, a hyperplane composed of the non-spatial coordinates of the corrosion parameters (H 2 S, Cl - , Fe 2+ , NH 3 , pH and flow rate) is constituted and the number of corrosion parameters is reduced While reducing the burden, predict corrosion rate more easily than in the prior art.
이하, 관련 도면을 참조하여 본원 발명을 더욱 상세히 설명한다.Hereinafter, the present invention will be described in more detail with reference to the related drawings.
우선, 운전 데이터는 수 년간 측정된다. 전술된 바와 같이 변수의 세트는 pH, Cl-, Fe2+, H2S, NH3, 및 유량을 포함하며, 그 외의 변수들도 포함할 수 있다. 부식 속도는 비정규적으로 매주 측정되지만, 그 외의 모든 변수들은 비정규적으로 매일 측정될 수 있다. First, the operation data is measured for several years. A set of variables as described above, is pH, Cl - includes, Fe 2+, H 2 S, NH 3, and flow rate, and may include also other parameters. The rate of corrosion is measured irregularly weekly, but all other variables can be measured daily irregularly.
6 개의 부식 변수들과 부식 속도를 조절하기 위하여, 구분적 선형화(piecewise linearization)가 이용되고, 각 변수들의 구간 별 평균값을 구하기 위하여 시간-가중 평균 기법이 사용된다. To control the six corrosion parameters and corrosion rate, piecewise linearization is used, and a time-weighted average technique is used to obtain the mean value of each interval.
도 2 는 본 발명에서 사용되는 구분적 선형화를 나타내는 그래프이다.FIG. 2 is a graph showing the segmental linearization used in the present invention. FIG.
도 2와 같이 측정 시간이 다른 경우에도 구분적 선형화를 통해 변수들을 선형 처리할 수 있다.As shown in FIG. 2, even when the measurement time is different, the variables can be linearly processed through the segmented linearization.
예측의 정확도를 높이고 처리를 용이하게 하기 위하여, 특정 이상점(outliner)이 배제될 수 있다. 예를 들어, 95% 신뢰도를 벗어나는 값들이 운전 데이터 중에서 배제될 수 있다. 본 명세서에서는 42개의 측정된 데이터 중 35개의 데이터가 학습용 샘플로 사용되었고, 나머지 7 개의 데이터는 검증용 샘플로 사용되었다. 하지만, 이는 본 발명의 설명을 위한 것이고, 더 적거나 더 많은 수의 데이터들이 사용될 수 있음은 물론이다. In order to increase the accuracy of the prediction and to facilitate processing, certain outliners may be excluded. For example, values outside the 95% confidence may be excluded from the driving data. In the present specification, 35 data out of 42 measured data were used as learning samples, and the remaining 7 data were used as verification samples. However, this is for the purpose of explanation of the present invention, and it goes without saying that fewer or more data can be used.
도 3 은 본 발명에서 사용되는 크리깅 모델을 설명하는 그래프이다.3 is a graph illustrating a kriging model used in the present invention.
도 3 에 나타나는 바와 같이, 크리깅 모델은 인접 지점 주위의 값들의 선형 조합을 이용하여 관심점의 값을 예측하는 기법이다. 크리깅 모델은 회귀 분석법에 대한 대안이라고 간주될 수 있다. 다른 통계적 모델과는 달리, 크리깅 모델은, 측정된 값들로부터 공간 상관 함수인 베리오그램(variogram)을 이용하여 값을 예측하는데, 베리오그램은 지구통계학에서 흔히 사용되는 관측값들의 공간적 연관성을 표현하는 함수로서, 베리오그램을 이용하여 공분산 값이 계산된다.As shown in FIG. 3, the kriging model is a technique for predicting the value of a point of interest using a linear combination of values around the adjacent points. The kriging model can be considered as an alternative to regression analysis. Unlike other statistical models, the kriging model predicts values from measured values using a spatial correlation function, variogram, which is a function that expresses the spatial association of observations commonly used in geostatistics , A covariance value is calculated using a variogram.
일반적으로 크리깅 모델은 지구 통계학에서 사용되는 기법으로, 공간 좌표 상에서 예측하려는 지점과 가까운 지점들의 값들로 선형 결합을 하여 예측하는 기법으로서, 지질학, 지리학, 기상학 외에도 전산 유체 역학(Computational Fluid Dynamics)과 같은 공간 정보와 관련된 분야로도 응용이 되고 있다. 예를 들어, 크리깅 모델에 대한 추가 설명을 본 명세서에 인용된 [논문 1] 내지 [논문 4]에서 찾을 수 있다. 이러한 논문들은 그 내용이 본원 발명에 원용되어 포함된다.In general, the Kriging model is a technique used in geostatistics. It is a technique for predicting linear combinations of values at points close to the point to be predicted on the spatial coordinates. In addition to geology, geography, and meteorology, It is also applied to fields related to spatial information. For example, further discussion of the kriging model can be found in [Thesis 1] through [Thesis 4] cited herein. These papers are incorporated herein by reference in their entirety.
크리깅 모델의 수학식 1과 같이 표현될 수 있다.Can be expressed as Equation (1) of the kriging model.
수학식 1에서, z*는 추정하려는 부식 속도이고, zi 는 i지점에서의 값이며, λi는 가중 파라미터이고, n은 인접 지점들의 개수이다. 가중 함수를 결정하기 위해서 여러 방정식을 풀어내야 한다. 가중 함수를 결정하기 위한 방법들은, 심플 크리깅(simple Kriging), 오디너리 크리깅(ordinary Kriging), 블록 크기링(block Kriging), 코-크리깅(co-Kriging), 또는 전역 크리깅(universal Kriging) 등 다양한 크리깅 모델 중 어느 것을 사용하는지에 따라 달라진다. 이용될 수 있는 크리깅 모델은 본 발명을 한정하지 않으며, 본 발명의 출원 전 공지된 모든 기법들이 사용될 수 있는 것이 이해되어야 한다. In Equation (1), z * is a corrosion rate to be estimated, z i is a value at an i- th point, λ i is a weighting parameter, and n is the number of adjacent points. Several equations must be solved to determine the weighting function. The methods for determining the weighting function may be implemented using various kriging techniques such as simple kriging, ordinary kriging, block kriging, co-kriging, or universal kriging. It depends on which model you use. It should be understood that the kriging model that may be used does not limit the invention and that all techniques known before the invention can be used.
이와 같이, 크리깅 모델을 이용한 부식 속도 예측은 입력받은 과거의 데이터들 중 유사한 환경에서의 부식 속도 값들의 선형 결합으로 구할 수 있다.Thus, the corrosion rate prediction using the Kriging model can be obtained by linear combination of the corrosion rate values in the similar environment among the inputted historical data.
설명의 편의를 위하여 심플 크리깅의 다양한 방법 중 심플 크리깅 모델을 예시하면, 수학식 2와 같다.For convenience of explanation, a simple kriging model among various methods of simple kriging is expressed by Equation (2).
수학식 2에서, C(i, j)는 인접 지점 i와 j 사이의 거리에서의 공분산값이고, c(i)는 목표지점과 인접 지점 i 사이의 거리에서의 공분산 값이다. 이러한 공분산 값들은 수학식 3 에서와 같이 베리오그램으로부터 유도된 공분산 함수를 사용해서 얻어질 수 있다.In Equation (2), C (i, j) is the covariance value at the distance between the adjacent points i and j, and c (i) is the covariance value at the distance between the target point and the adjacent point i. These covariance values can be obtained using a covariance function derived from the variogram as in Equation (3).
수학식 3에서, h는 두 지점 사이의 이격거리이고, γ는 베리오그램이다. Sill은 제로의 거리에서의 공분산 값 또는 점들의 쌍의 최대 변이이다. 즉, 이것은 거리가 증가함에 따라 제로로 감소하는 공분산 값(C(h)) 때문에 무한대의 거리로 커지는 베리오그램의 상한이다. 가중 파라미터 is가 베리오그램으로부터 얻어지기 때문에, 베리오그램의 정확성이 예측 성능 전체에 큰 영향을 미친다. In Equation (3), h is the separation distance between two points, and? Is the variogram. Sill is the covariance value at a distance of zero or the maximum variation of a pair of points. That is, this is the upper limit of the variogram that grows to an infinite distance due to the covariance value C (h), which decreases to zero as the distance increases. Since the weighted parameter is is obtained from the variogram, the accuracy of the variogram greatly affects the predicted performance.
베리오그램은 데이터의 상관을 거리로 기술한다. 실험적인 베리오그램은 다음 수학식 4와 같이 거리가 h인 점들의 쌍들 사이의 평균 제곱차(mean square difference)의 절반이다.The variogram describes the correlation of data as distance. The experimental variogram is half of the mean square difference between pairs of points with distance h, as shown in equation (4).
수학식 4에서, |N(h)|는 쌍들의 개수이고, z는 어느 지점에서의 종속 변수값이며, (i, j)는 거리 h만큼 떨어진 점들의 쌍이다. 이론적 베리오그램은 실험적인 베리오그램으로부터 얻을 수 있다. Is the number of pairs, z is the dependent variable value at some point, and (i, j) is a pair of points spaced by distance h. Theoretical variogram can be obtained from an experimental variogram.
일 예로서, 본 명세서에서는 수학식 5 와 같은 오디너리 크리깅을 사용하여 본 발명을 설명하지만, 본 발명은 이러한 기법으로 한정되지 않고, 다양한 크리깅 기법들이 채용될 수 있음이 명백하다.As an example, although the present invention is described herein using an orininal kriging as in
또한, 본 명세서에서는 전술한 바와 같이, 부식 변수로 이루어진 초평면(hyperplane)을 사용하여 크리깅 모델을 적용한다.Also, in this specification, as described above, a kriging model is applied using a hyperplane made of a corrosion parameter.
크리깅 기법으로 예측한 부식 속도의 정확도를 다중 선형 회귀 분석(Multi Linear Regression; MLR) 및 인공 신경망(Artificial Neural Network; ANN)과 같은 기법들과 비교한다.The accuracy of the predicted corrosion rate with the Kriging technique is compared with techniques such as MLR and ANN.
다중 선형 회귀 분석을 사용하면, 종속 변수에 대한 각각의 독립 변수의 영향을 분석하고, 독립 변수와 종속 변수 사이의 관련성을 인식하는 것이 수월하게 이루어질 수 있다. 하지만, 다중 선형 회귀 분석은 비선형성과 복잡성을 가지는 거동을 가지는 자연 현상을 예측하는 데에는 충분하지 않다. Using multiple linear regression analysis, it is easy to analyze the effect of each independent variable on the dependent variable and to recognize the relationship between the independent variable and the dependent variable. However, multiple linear regression analysis is not sufficient to predict natural phenomena with nonlinear and complex behavior.
인공 신경망 기법은 인간 두뇌의 신경망을 모사하는 인공 지능 방법이다. 이러한 방법은 복잡한 비선형계를 풀어내는 데에 강력하기 때문에 에너지 소비, 유가, 열전도율, 공정 최적화 등의 분야에서 널리 사용된다. 그러나, 인공 신경망은 몇 가지 단점을 가지는데, 첫 번째는 과맞춤(over-fitting)의 문제점이다. 데이터가 샘플 구간에서는 잘 맞지만 다른 지역에 대해서는 양호한 예측 성능을 보이지 않는 것이 과맞춤 문제점이다. 다른 문제점은, 자유도가 너무 높아서 선택 변수가 너무 많다는 것이다.Artificial neural network technique is an artificial intelligence method that simulates the neural network of the human brain. Because these methods are robust in solving complex nonlinear systems, they are widely used in energy consumption, oil prices, thermal conductivity, and process optimization. However, artificial neural networks have some disadvantages, the first being the problem of over-fitting. The problem is that the data fit well in the sample interval but do not show good prediction performance in other regions. Another problem is that the degree of freedom is too high and there are too many optional variables.
따라서, 변수의 개수를 줄이는 것이 필요하다. 그 첫 번째 이유는 부식 속도 예측을 위해서 6 개의 부식 변수를 모두 사용하면 만족할만한 예측 결과가 나오지 않는다는 것이다. 어떤 변수는 다른 것들에 비하여 부식 속도에 미치는 영향이 적을 수 있는데도, 무리하게 해당 변수를 포함하여 모델링을 수행할 필요가 없다. Therefore, it is necessary to reduce the number of variables. The first reason is that using all six corrosion variables to predict corrosion rate does not yield satisfactory predictions. Some variables may have less impact on corrosion rate than others, but it is not necessary to model them involuntarily.
크리깅 모델의 경우에, 6차원 초평면에 42 개의 데이터 세트를 적용하였지만, 차원이 높아질수록 적절한 베리오그램을 찾기 위해서 더 많은 샘플 데이터가 필요하기 때문에 6차원은 적합하지 않을 수 있다. 그러므로, 본 발명에서는 데이터 세트를 두 개 감소시켜서 초평면을 구성하였다. 사용되는 변수의 개수는 신뢰도가 높아질수록 많아지게 되는데, 일반적으로 4차원 초평면을 구성하여 수행한 결과를 가지고 신뢰 구간을 측정해 본 결과 충분히 양호한 결과를 얻을 수 있다는 것을 알 수 있다. In the case of the kriging model, the 42 data sets were applied to the 6-dimensional hyperplane, but the higher the dimension, the more dimensional data may be not suitable because more sample data is needed to find the proper variogram. Therefore, in the present invention, the hyperplane is constructed by reducing the number of data sets by two. The number of variables used increases as the reliability increases. In general, it can be seen that satisfactory results can be obtained by measuring the confidence interval with the results obtained by constructing the 4-dimensional hyperplane.
35 개의 데이터 세트를 사용하여 6 개의 부식 변수가 개발되고 검증되었다. 모든 변수들은 정규화되어 모델에 적용된다.Six corrosion parameters were developed and validated using 35 data sets. All variables are normalized and applied to the model.
도 4 는 실험적이고 이론적인 베리오그램을 나타낸다.Figure 4 shows an experimental and theoretical variogram.
그래프에서 이격 거리는 학습용 샘플지점들 사이에서의 6 개의 부식 변수들의 유사도를 나타낸다. 35 개의 데이터 세트만을 사용하였기 때문에, 0부터 0.2 까지의 이격 거리 범위에서는 베리오그램 실험값이 존재하지 않는다. 그러므로, 이러한 구간에 베리오그램 실험값이 존재하지 않는다는 것은 이론적인 베리오그램이 정확하지 않다는 것을 나타낸다. 도 4 는 0.6 의 이격 거리에 도달할 때까지 초평면 상의 이격 거리에는 공간적 관련성이 존재한다는 것을 나타낸다(이것은 range라고 불린다). 이러한 range 밖에서의 베리오그램 값은 이격 거리와는 아무런 관련성을 가지지 않는다. 최적화 툴을 사용하여 파라미터를 추정하면 이론적인 베리오그램이 얻어진다. 이론적 베리오그램으로서 지수 모델이 사용되는데, 이것은 각각 0.574의 range와 0.087의 sill값이 된다. The spacing in the graph represents the similarity of the six corrosion variables between the training sample points. Since only 35 data sets are used, there is no variogram experiment in the range of 0 to 0.2. Therefore, the absence of a variogram experiment in this interval indicates that the theoretical variogram is not correct. Figure 4 shows that there is a spatial relationship (this is called a range) in the separation distance on the hyperplane until a separation distance of 0.6 is reached. The variogram values outside this range have no relation to the separation distance. Estimating parameters using an optimization tool yields a theoretical variogram. An exponential model is used as the theoretical variogram, which has a range of 0.574 and sill of 0.087, respectively.
가장 가까운 3개, 5개, 그리고 7 개의 지점과 검증 구간 내의 모든 점들을 사용한 크리깅 모델의 평균 자승 오차(mean square error; MSE)는 각각 0.541, 0.625, 0.572, 0.351 이다. 가장 가까운 3개, 5개, 그리고 7 개의 지점에서 얻은 결과가 더 안 좋은 이유는 0.2 의 거리 내에 근접 지점이 없기 때문이다. 도 5 의 (a)에 표시되는 바와 같이, 학습 구간에서의 MSE는, 다른 모델들과 달리 해당 구간 내에서의 측정치와 예측치가 크리깅 모델에서 정확하게 일치하기 때문에 정확히 0이 된다. The mean square error (MSE) of the Kriging model using the closest 3, 5, and 7 points and all points within the validation interval are 0.541, 0.625, 0.572, and 0.351, respectively. The results from the nearest 3, 5, and 7 points are worse because there is no near point within a distance of 0.2. As shown in (a) of FIG. 5, unlike the other models, the MSE in the learning section is exactly zero because the measured and predicted values in the corresponding section exactly match in the kriging model.
도 5 의 (b)는 MLR 모델 결과를 나타낸다. 검증 샘플들의 MSE는 0.420이고, 학습 샘플들의 MSE는 1.126이다. 비록 파라미터의 세트가 학습 샘플들에 의해서 결정되었지만, 검증 구간의 결과는 MSE보다 더 양호한 결과를 보여 준다. 그 이유는, 부식 속도의 범위 때문이다: 학습 구간에서는 0.456 MPY 내지 4.867 MPY 이지만, 검증 구간에서는 1.217 MPY 내지 2.920 MPY 이다. 더 나아가, 예측 결과에는 측정값들 사이에 트렌드도 잘 보이지 않는다.Figure 5 (b) shows the MLR model results. The MSE of the verification samples is 0.420 and the MSE of the training samples is 1.126. Although the set of parameters is determined by the training samples, the results of the verification interval show better results than the MSE. The reason is due to the range of corrosion rates: 0.456 MPY to 4.867 MPY in the learning period, but 1.217 MPY to 2.920 MPY in the verification period. Furthermore, trends in measured values are not well visible in the prediction results.
40 개의 상이한 타입의 인공 신경망 모델들은 은닉 뉴런, 바이어스 뉴런의 존재/비존재, 및 에포크의 조합으로 이루어진다. 표 1은 MSE에 대한 인공 신경망 모델의 결과를 나타낸다. Forty different types of artificial neural network models consist of concealed neurons, presence / absence of bias neurons, and a combination of epochs. Table 1 shows the results of the artificial neural network model for MSE.
은닉 뉴런의 개수가 많아질수록, 학습 구역에서의 MSE는 모든 경우에 감소하는 경향성을 가진다. 반면에, 4 개 내지 6 개의 은닉 뉴런까지는, 검증 구역에서의 MSE는 증가되는 경향을 가지며, 은닉 뉴런이 더 많아지면 경향성에 일관성이 없어진다. 그러나, 은닉 뉴런의 개수가 많아지면 더 많은 파라미터가 결정되어야 한다. 바이어스 뉴런이 포함된 모델은 일반적으로 바이어스 뉴런이 없는 모델보다 열약한 예측 성능을 보여준다. 또한, 에포크의 개수도 과맞춤 문제점에 영향을 준다. As the number of hidden neurons increases, the MSE in the learning zone has a tendency to decrease in all cases. On the other hand, up to four to six concealed neurons, the MSE in the probe area tends to increase, and the more concealed neurons become inconsistent in the tendency. However, as the number of concealed neurons increases, more parameters must be determined. Models with bias neurons typically exhibit poor predictive performance than models without bias neurons. It also affects the number and alignment of epochs.
가장 바람직한 인공 신경망 모델은 바이어스 뉴런이 없이, 3 개의 은닉 뉴런과 10,000 개의 에포크를 가지는 모델이라고 볼 수 있으며, 이는 도 5의 (c)에 표시된다. 다중 선형 회귀 분석모델과는 달리 측정된 부식 속도에 대한 양호한 예측값을 제공하고, 다중 선형 회귀 분석과 비교할 때 주어진 학습 구간과 검증 구간 모두에서 MSE 값이 절반 밖에 되지 않는다. The most desirable artificial neural network model is a model with 3 hidden neurons and 10,000 epochs without bias neurons, which is shown in Figure 5 (c). Unlike the multiple linear regression model, it provides a good predictor of the measured corrosion rate, and compared to polynomial regression analysis, the MSE value is only half that of the given learning and verification periods.
도 5 의 다중 선형 회귀 분석 모델과 인공 신경망 모델과 비교하여, 크리깅 모델은 다중 선형 회귀 분석 모델과 유사하다고 볼 수 있다. 또한, 크리깅 모델을 살펴보면, 검증 구간 내의 약 1.5 MPY의 작은 부식 속도점에서, 측정값과 예측값 사이에 큰 차이가 있다. 그러므로, 6 개의 부식 변수가 있는 경우에는 인공 신경망 모델이 최적의 결과를 나타낸다. 하지만, 크리깅 모델은 더 정확한 베리오그램이 얻어진다면 개선될 수 있다. 그러려면 데이터 세트를 더 많이 사용하거나 변수들의 개수를 줄여야 한다.Compared with the multiple linear regression analysis model and the artificial neural network model of FIG. 5, the kriging model is similar to the multiple linear regression model. Also, looking at the Kriging model, there is a large difference between the measured and predicted values at a small corrosion rate of about 1.5 MPY in the validation period. Therefore, an artificial neural network model shows optimal results when there are six corrosion variables. However, the kriging model can be improved if more accurate variograms are obtained. To do this, you need to use more data sets or reduce the number of variables.
이와 같이, 6 개의 부식 변수들이 4 개로 감소되면 5% 미만의 유의 확률(significance probability)이 얻어진다. 6 개의 변수들 중에서 NH3와 Fe2+가 제거될 수 있다. 그리고, 나머지 변수인 H, Cl-, H2S, 유량들이 사용된다.Thus, if the six corrosion parameters are reduced to four, a significance probability of less than 5% is obtained. Of the six variables, NH 3 and Fe 2+ can be removed. The remaining variables, H, Cl - , H 2 S, and flow rates are used.
4 개의 변수가 사용될 경우의 크리깅 모델에 대한 실험 베리오그램과 이론적 베리오그램은 도 6 에 표시된다. 도 6 에는 이격 거리 0.2 내에 인접 점들이 존재하고, range 내의 점들의 개수가 변수가 6 개일 경우에 비하여 증가되었다. 결과적으로, 실험 베리오그램은 이격 거리가 증가함에 따라 range 구간 내에 명백한 증가 경향을 나타낸다. 이론적 베리오그램인 가우스 모델이 실험 베리오그램과 잘 맞는다. 모델의 파라미터 세트는 파라미터 추정을 통해 얻어진다. 6 개의 부식 변수가 사용되는 것에 비하여, range는 0.574에서 0.585로 증가했다. 또한, sill도 0.087에서 0.116으로 증가했다.The experimental variogram and the theoretical variogram for the Kriging model when four variables are used are shown in FIG. In Fig. 6, adjacent points exist within a distance of 0.2, and the number of points in the range is increased as compared with the case where the number of points is 6. As a result, the experimental variogram exhibits an apparent increasing tendency within the range interval as the separation distance increases. The Gaussian model, the theoretical variogram, fits well with the experimental variogram. The parameter set of the model is obtained through parameter estimation. The range increased from 0.574 to 0.585, compared to the use of six corrosion variables. The sill also increased from 0.087 to 0.116.
도 7 의 (a)는 변수가 네 개의 경우에 크리깅 모델들 중 최적인 크리깅 모델을 도시한다. 이 때, 가장 가까운 3개, 5개, 7 개 점들과 관련된 MSE는 각각 0.377, 0.183, 0.375이다. 가장 가까운 5 개의 점들을 사용한 크리깅 모델이 가장 낮은 MSE 값을 나타낸다. 도 5의 (a)와 도 7 의 (a)에서 볼 수 있는 것처럼, 검증 구간에서의 예측된 값들의 범위는 변수가 6 개 인 경우에 비하여 더 커진다. 그러므로, 이러한 크리깅 모델은 예측 성능을 개선하고, 부식 속도가 낮은 경우에도 신뢰성이 높다.Figure 7 (a) shows the optimal kriging model among the kriging models with four variables. The MSEs associated with the nearest 3, 5, and 7 points are 0.377, 0.183, and 0.375, respectively. The Kriging model using the closest 5 points shows the lowest MSE value. As can be seen in Figures 5 (a) and 7 (a), the range of predicted values in the verification interval is larger than when the number of variables is six. Therefore, this kriging model improves predictive performance and is highly reliable even at low rates of corrosion.
변수가 4 개의 경우의 다중 선형 회귀 분석 모델은 검증 구역과 학습 구역에서 각각 0.448과 1.126의 MSE를 나타낸다. 4 개의 변수를 가진 모델의 예측 성능은 변수가 6 개인 모델보다 차이가 적은 이유는, 부식 속도와 상관성이 적은 데이터 세트를 얻기 위하여, 제거된 변수들의 세트가 유의도 확률을 적용하여 결정되었기 때문이다.The multiple linear regression model with four variables shows MSEs of 0.448 and 1.126 in the verification and learning zones, respectively. The reason that the predictive performance of a model with four variables is smaller than the model with six variables is that the set of removed variables was determined by applying the significance probability in order to obtain a data set with little correlation with the corrosion rate .
4 개의 부식 변수를 가지는 인공 신경망 모델의 결과가 표 2 에 표시된다. 변수가 6 개인 경우에 비하여, 변수가 4 개일 경우 매우 개선된 결과를 보여준다. 그러나, MSE에 대해서 살펴보면, 최적의 결과는 변수가 6 개인 경우와 유사하며 그 차이는 0.209:0.232와 같다.The results of the artificial neural network model with four corrosion variables are shown in Table 2. Compared to the case of six variables, four variables show very improved results. However, looking at the MSE, the optimal result is similar to the case of six variables, with a difference of 0.209: 0.232.
각각의 경우에 모델별로 최소 MSE 값을 구한 것이 표 3 에 표시된다. 부식 변수가 4 개이고 가장 가까운 5 개의 점들을 사용하는 크리깅 모델이 이들 중에서 가장 양호한 결과를 나타낸다.Table 3 shows the minimum MSE value for each model in each case. The Kriging model, which has four corrosion variables and uses the closest five points, has the best results.
도 8 은 본 발명의 다른 실시예에 의한 유류 설비의 부식 속도 예측 장치를 개략적으로 나타내는 블록도이다.8 is a block diagram schematically showing an apparatus for predicting corrosion rate of a oil plant according to another embodiment of the present invention.
본 발명의 다른 면에 의한 유류 설비의 부식 속도 예측 장치(800)는, 부식 데이터 획득부(810), 이상점 제거부(820), Q-test 부(840), 중앙 제어부(850), 부식 속도 예측부(860), 상관 관계 분석부(880), 및 결과 제공부(890)를 포함한다.The apparatus for predicting erosion rate of a plant according to another aspect of the present invention includes a corrosion
부식 데이터 획득부(810)는 유류 설비의 부식에 대해 측정된 부식 데이터를 획득한다. 전술된 바와 같이, 부식 데이터는 다양한 시간 간격으로 측정된 데이터일 수 있다.The corrosion
그러면, 이상점 제거부(820)가 획득된 부식 데이터로부터 이상점을 제거한다. 이를 위하여, 이상점 제거부(820)는 부식 속도와 부식 변수들의 측정 시간을 고려하여 부식 데이터에 구분적 선형화 및 시간 가중 평균화 기법을 적용할 수 있다. 이와 같이 처리된 부식 데이터는 Q-test 부(840)로 제공되고, Q-test 부(840)는 이러한 데이터에 Q-test를 수행하여 소정 임계를 벗어나는 부식 데이터를 제거한다.Then, the
중앙 제어부(850)는 유류 설비의 부식 속도 예측 장치(800)의 구성 요소들의 동작을 제어한다.The central control unit 850 controls the operation of the components of the corrosion
이상점이 제거되면, 부식 속도 예측부(860)는 이상점이 제거된 부식 데이터를 크리깅 모델(Kriging model)에 적용하여 부식 속도를 예측한다. 부식 속도의 예측을 용이화하기 위하여, 상관 관계 분석부(880)는 크리깅 모델에 적용될 부식 변수할 부식 변수와 부식 속도의 상관 관계를 유의 확률로 분석하고, 분석 결과에 따라 유의도가 낮은 부식 변수를 제거하여 크리깅 모델의 초평면의 차원을 낮출 수 있다. 이와 같이 초평면의 차원이 낮아지면, 부식 속도를 더욱 신속하고 용이하게 예측할 수 있는 장점이 있다.If the abnormal point is removed, the corrosion
부식 속도가 예측되면, 결과 제공부(890)가 예측된 결과를 다중 선형 회귀 분석 결과 및 인공 신경망 분석 결과와 함께 사용자에게 제공한다.When the erosion rate is predicted, the
이와 같이, 본 발명에서는 이상점이 제거된 42개 데이터를 중 35개 데이터를 학습에 사용하였고, 7개 데이터를 검증하는데 사용하였다. 그 결과 크리깅 모델 모델은 단순한 다중선형회귀분석보다는 정확하지만 인공신경망보다는 평균제곱오차가 좋지 않은 것으로 나타났으며, 정확성 향상을 위하여 부식 변수와 부식 속도의 선형 관계를 유의 확률로 분석한 결과, 상관 관계가 적은 부식 변수인 Fe2+, NH3를 제거한 4차원의 초평면에서 크리깅 모델을 적용하였다.As described above, in the present invention, 35 data out of 42 data whose anomalous points have been removed are used for learning and 7 data are used for verification. As a result, Kriging model model is more correct than simple polynomial regression analysis but mean square error is worse than artificial neural network. In order to improve accuracy, linear correlation of corrosion rate and corrosion rate is analyzed with significance, The Kriging model was applied to the 4 - dimensional hyperplane where Fe 2+ and NH 3 were removed.
이와 같이, 차원의 축소로 35개 데이터가 초평면을 더 잘 나타낼 수 있게 됨에 따라 크리깅 모델의 예측 결과는 상당히 개선이 되어 다중 선형 회귀 분석 및 인공 신경망보다 정확한 결과를 나타냈다. 따라서, 크리깅 모델을 비공간 좌표계에 적용을 할 때에는 데이터 수에 비해 변수의 수가 적어 정확한 예측이 힘든 경우에 효과적이라는 것을 확인할 수 있다.As a result of the reduction of dimensions, 35 data can better represent the hyperplane, so the predictions of the Kriging model have been improved significantly and more accurate than multilinear regression analysis and artificial neural network. Therefore, it can be confirmed that when the kriging model is applied to the non-spatial coordinate system, the number of variables is smaller than the number of data, which is effective when the accurate prediction is difficult.
본 발명은 도면에 도시된 실시예를 참고로 설명되었으나 이는 예시적인 것에 불과하며, 본 기술 분야의 통상의 지식을 가진 자라면 이로부터 다양한 변형 및 균등한 타 실시예가 가능하다는 점을 이해할 것이다.While the present invention has been particularly shown and described with reference to exemplary embodiments thereof, it is evident that many alternatives, modifications and variations will be apparent to those skilled in the art.
또한, 본 발명에 따르는 방법은 컴퓨터로 읽을 수 있는 기록 매체에 컴퓨터가 읽을 수 있는 코드로서 구현하는 것이 가능하다. 컴퓨터가 읽을 수 있는 기록 매체는 컴퓨터 시스템에 의하여 읽혀질 수 있는 데이터가 저장되는 모든 종류의 기록 장치를 포함할 수 있다. 컴퓨터가 읽을 수 있는 기록 매체의 예로는 ROM, RAM, CD-ROM, 자기 테이프, 플로피 디스크, 광 데이터 저장 장치 등이 있으며, 또한 캐리어 웨이브(예를 들어 인터넷을 통한 전송)의 형태로 구현되는 것도 포함한다. 또한 컴퓨터가 읽을 수 있는 기록 매체는 네트워크로 연결된 분산 컴퓨터 시스템에 의하여 분산 방식으로 실행될 수 있는 컴퓨터가 읽을 수 있는 코드를 저장할 수 있다. In addition, the method according to the present invention can be embodied as computer-readable code on a computer-readable recording medium. A computer-readable recording medium may include any type of recording device that stores data that can be read by a computer system. Examples of the computer-readable recording medium include a ROM, a RAM, a CD-ROM, a magnetic tape, a floppy disk, an optical data storage device, and the like, and may be implemented in the form of a carrier wave (for example, transmission via the Internet) . The computer readable recording medium may also store computer readable code that may be executed in a distributed manner by a distributed computer system connected to the network.
본 명세서에서 사용되는 용어에서 단수의 표현은 문맥상 명백하게 다르게 해석되지 않는 한 복수의 표현을 포함하는 것으로 이해되어야 하고, "포함한다" 등의 용어는 설시된 특징, 수, 단계, 동작, 구성요소, 부분품 또는 이들을 조합한 것이 존재함을 의미하는 것이지, 하나 또는 그 이상의 다른 특징들이나 개수, 단계 동작 구성요소, 부분품 또는 이들을 조합한 것들의 존재 또는 부가 가능성을 배제하지 않는 것으로 이해되어야 한다. 그리고, 명세서에 기재된 "...부", "...기", "모듈", "블록" 등의 용어는 적어도 하나의 기능이나 동작을 처리하는 단위를 의미하며, 이는 하드웨어나 소프트웨어 또는 하드웨어 및 소프트웨어의 결합으로 구현될 수 있다. As used herein, the singular " include " should be understood to include a plurality of representations unless the context clearly dictates otherwise, and the terms "comprises & , Parts or combinations thereof, and does not preclude the presence or addition of one or more other features, integers, steps, components, components, or combinations thereof. The terms "part", "unit", "module", "block", and the like described in the specification mean units for processing at least one function or operation, And a combination of software.
따라서, 본 실시예 및 본 명세서에 첨부된 도면은 본 발명에 포함되는 기술적 사상의 일부를 명확하게 나타내고 있는 것에 불과하며, 본 발명의 명세서 및 도면에 포함된 기술적 사상의 범위 내에서 당업자가 용이하게 유추할 수 있는 변형 예와 구체적인 실시예는 모두 본 발명의 권리범위에 포함되는 것이 자명하다고 할 것이다.Therefore, it should be understood that the present invention and the drawings attached hereto are only a part of the technical idea included in the present invention, and that those skilled in the art will readily understand the technical ideas included in the specification and drawings of the present invention It is obvious that all the variations and concrete examples that can be deduced are included in the scope of the present invention.
본 발명은 정유 공장의 유류 시설들의 부식 속도를 예측하기 위하여 사용될 수 있다.The present invention can be used to predict the erosion rate of oil installations in refineries.
800 유류 설비의 부식 속도 예측 장치
810 부식 데이터 획득부 820 이상점 제거부
840 Q-test 부 850 중앙 제어부
860 부식 속도 예측부 880 상관 관계 분석부
890 결과 제공부800 Corrosion rate prediction device of oil facility
810 Corrosion
840 Q-test unit 850 Central control unit
860 Corrosion
890 Results provided
Claims (12)
상기 유류 설비의 부식 데이터를 획득하는 부식 데이터 획득부;
획득된 부식 데이터로부터 이상점을 제거하는 이상점 제거부; 및
이상점이 제거된 부식 데이터를 크리깅 모델(Kriging model)에 적용하여 부식 속도를 예측하는 부식 속도 예측부를 포함하여 구성되고,
상기 이상점 제거부는,
부식 속도와 부식 변수들의 측정 시간을 고려하여 상기 부식 데이터에 구분적 선형화(piecewise linearization) 및 시간 가중 평균화 기법을 적용하고,
적용된 결과에 Q-test를 수행하여 소정 임계를 벗어나는 부식 데이터를 제거하도록 구성되며, 상기 부식 속도 예측부는,
상기 크리깅 모델에 적용될 부식 변수들을 획득하고,
획득된 부식 변수와 부식 속도의 상관 관계를 유의 확률로 분석하며,
분석 결과에 따라 유의도가 낮은 부식 변수를 제거하여 크리깅 모델의 초평면의 차원을 낮추도록 구성되고,
상기 부식 속도 예측부는,
를 사용하여 부식 속도를 예측하고,
여기에서, z*는 예측하려는 부식 속도, zi 는 예측하려는 환경과 유사한 환경에서 측정된 부식 데이터, n은 상기 유사한 환경에서 측정된 부식 데이터의 수, λi는 zi 에 적용되는 가중치를 나타내며,
예측된 결과를 다중 선형 회귀 분석 결과 및 인공 신경망 분석 결과 중 적어도 하나와 함께 제공하는 결과 제공부를 더 포함하고, 상기 유류 설비는 정류 공장의 증류탑 상부에 해당하는 것을 특징으로 하는 유류 설비의 부식 속도 예측 장치.
An apparatus for predicting erosion rate of an oil plant,
A corrosion data obtaining unit for obtaining corrosion data of the oil facility;
An abnormal point removing unit for removing an abnormal point from the obtained corrosion data; And
And a corrosion rate predicting unit for predicting a corrosion rate by applying the corrosion data from which the abnormal point is removed to a kriging model,
Wherein the abnormal point removing unit comprises:
Applying piecewise linearization and time weighted averaging techniques to the corrosion data taking into account the corrosion rate and the measurement time of the corrosion parameters,
Wherein the corrosion rate estimator is configured to perform a Q-test on the applied result to remove corrosion data that deviates from a predetermined threshold,
Obtaining the corrosion parameters to be applied to the kriging model,
The correlation between the obtained corrosion parameters and the corrosion rate is analyzed with significance probability,
According to the results of the analysis, it is constructed to lower the dimension of the hyperplane of the kriging model by removing the low -
The corrosion rate predicting unit may include:
The corrosion rate is predicted by using < RTI ID = 0.0 >
Where z * is the corrosion rate to be predicted, z i is the corrosion data measured in an environment similar to the environment to be predicted, n is the number of corrosion data measured in the similar environment, and λ i is the weight applied to z i ,
Further comprising a result providing unit for providing the predicted result together with at least one of the result of the multi-linear regression analysis and the result of the artificial neural network analysis, and the oil facility corresponds to the upper part of the distillation tower of the rectification plant. Device.
상기 유류 설비의 부식 데이터를 획득하는 단계;
획득된 부식 데이터로부터 이상점을 제거하는 단계; 및
이상점이 제거된 부식 데이터를 크리깅 모델(Kriging model)에 적용하여 부식 속도를 예측하는 단계를 포함하여 구성되고,
상기 이상점을 제거하는 단계는,
부식 속도와 부식 변수들의 측정 시간을 고려하여 상기 부식 데이터에 구분적 선형화(piecewise linearization) 및 시간 가중 평균화 기법을 적용하는 단계; 및
적용된 결과에 Q-test를 수행하여 소정 임계를 벗어나는 부식 데이터를 제거하는 단계를 포함하며,
상기 부식 속도를 예측하는 단계는,
상기 크리깅 모델에 적용될 부식 변수들을 획득하는 단계;
획득된 부식 변수와 부식 속도의 상관 관계를 유의 확률로 분석하는 단계; 및
분석 결과에 따라 유의도가 낮은 부식 변수를 제거하여 크리깅 모델의 초평면의 차원을 낮추는 단계를 포함하고,
상기 부식 속도를 예측하는 단계는,
을 사용하여 부식 속도를 예측하는 단계를 더 포함하며,
여기에서, z*는 예측하려는 부식 속도, zi 는 예측하려는 환경과 유사한 환경에서 측정된 부식 데이터, n은 상기 유사한 환경에서 측정된 부식 데이터의 수, λi는 zi 에 적용되는 가중치를 나타내며, 예측된 결과를 다중 선형 회귀 분석 결과 및 인공 신경망 분석 결과 중 적어도 하나와 함께 제공하는 단계를 더 포함하며,
상기 유류 설비는 정류 공장의 증류탑 상부인 것을 특징으로 하는 유류 설비의 부식 속도 예측 방법.A method for predicting corrosion rate of an oil plant using the corrosion rate predicting device of the oil plant according to claim 1,
Obtaining corrosion data of the oil facility;
Removing an abnormal point from the obtained corrosion data; And
And a step of applying the corrosion data from which the abnormal point is removed to the kriging model to predict the corrosion rate,
The step of removing the abnormal point may include:
Applying piecewise linearization and time weighted averaging techniques to the corrosion data in consideration of the corrosion rate and the measurement time of the corrosion parameters; And
Performing a Q-test on the applied result to remove corrosion data that deviates from a predetermined threshold,
The step of predicting the corrosion rate comprises:
Obtaining corrosion parameters to be applied to the kriging model;
Analyzing the correlation between the obtained corrosion parameter and the corrosion rate with a significant probability; And
And a step of lowering the dimension of the hyperplane of the kriging model by removing the corrosion parameter having a low degree of significance according to the analysis result,
The step of predicting the corrosion rate comprises:
Further comprising the step of predicting the rate of corrosion,
Where z * is the corrosion rate to be predicted, z i is the corrosion data measured in an environment similar to the environment to be predicted, n is the number of corrosion data measured in the similar environment, and λ i is the weight applied to z i , And providing the predicted result together with at least one of a polynomial regression analysis result and an artificial neural network analysis result,
Wherein the oil facility is above the distillation tower of the rectification plant.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020160142752A KR101896115B1 (en) | 2016-10-31 | 2016-10-31 | Method and apparatus for estimating corrosion rate of oil plant facility |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020160142752A KR101896115B1 (en) | 2016-10-31 | 2016-10-31 | Method and apparatus for estimating corrosion rate of oil plant facility |
Publications (2)
Publication Number | Publication Date |
---|---|
KR20180047073A KR20180047073A (en) | 2018-05-10 |
KR101896115B1 true KR101896115B1 (en) | 2018-09-06 |
Family
ID=62184907
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR1020160142752A KR101896115B1 (en) | 2016-10-31 | 2016-10-31 | Method and apparatus for estimating corrosion rate of oil plant facility |
Country Status (1)
Country | Link |
---|---|
KR (1) | KR101896115B1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110175740A (en) * | 2019-04-15 | 2019-08-27 | 大唐环境产业集团股份有限公司 | A kind of scene analysis implementation method based on Kriging agent model |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6566172B1 (en) * | 2018-01-05 | 2019-08-28 | Jfeスチール株式会社 | Corrosion amount prediction method of metal material, selection method of metal material and corrosion amount prediction device of metal material |
WO2021100341A1 (en) * | 2019-11-18 | 2021-05-27 | Jfeスチール株式会社 | Metal material corrosion amount mapping method, metal material selection method, and metal material corrosion amount mapping device |
CN118096135B (en) * | 2024-04-28 | 2024-07-30 | 山东港源管道物流有限公司 | Corrosion prevention management system for port tank field equipment |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004252781A (en) | 2003-02-20 | 2004-09-09 | Japan Energy Corp | Modeling data forming method of corrosion velocity estimation system and device |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101084637B1 (en) * | 2009-04-01 | 2011-11-17 | 연세대학교 산학협력단 | Method for preventing corrosion in chemical plants |
KR101199462B1 (en) * | 2010-09-09 | 2012-11-09 | 중앙대학교 산학협력단 | Forecasting apparatus for carbon steel pipe, forecasting method therefor and recording medium thereof |
-
2016
- 2016-10-31 KR KR1020160142752A patent/KR101896115B1/en active IP Right Grant
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004252781A (en) | 2003-02-20 | 2004-09-09 | Japan Energy Corp | Modeling data forming method of corrosion velocity estimation system and device |
Non-Patent Citations (2)
Title |
---|
Kyungjae Tak, Kriging models for forecasting crude unit overhead corrosion, Korean J. Chem. Eng., Vol.33, No.7, 2016년06월, pp.1999-2006 |
이선규 외, 가스전의 추출속도를 고려한 C3MR 공정의 동력기 선택 전략, 한국가스학회지제20권제1호, 한국가스학회, 2016.02.pp7~12 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110175740A (en) * | 2019-04-15 | 2019-08-27 | 大唐环境产业集团股份有限公司 | A kind of scene analysis implementation method based on Kriging agent model |
Also Published As
Publication number | Publication date |
---|---|
KR20180047073A (en) | 2018-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR101896115B1 (en) | Method and apparatus for estimating corrosion rate of oil plant facility | |
Kratzert et al. | NeuralHydrology–interpreting LSTMs in hydrology | |
Fortuna et al. | Soft sensors for product quality monitoring in debutanizer distillation columns | |
US20190318288A1 (en) | Computer Systems And Methods For Performing Root Cause Analysis And Building A Predictive Model For Rare Event Occurrences In Plant-Wide Operations | |
Liu et al. | Development of a novel adaptive soft-sensor using variational Bayesian PLS with accounting for online identification of key variables | |
Zhang et al. | Integration of machine learning algorithms with Gompertz Curves and Kriging to estimate resources in gold deposits | |
Li et al. | Adaptive selective ensemble-independent component analysis models for process monitoring | |
Maldonado-Cruz et al. | Fast evaluation of pressure and saturation predictions with a deep learning surrogate flow model | |
Geiges | Efficient concepts for optimal experimental design in nonlinear environmental systems | |
EP3058338B1 (en) | Detection of corrosion rates in process piping and vessels | |
Diaz et al. | Time series modeling of methane gas in underground mines | |
Najafabadipour et al. | The innovative combination of time series analysis methods for the forecasting of groundwater fluctuations | |
Shao et al. | Soft sensor development for multimode processes based on semisupervised Gaussian mixture models | |
Borghesan et al. | Forecasting of process disturbances using k-nearest neighbours, with an application in process control | |
Tan | Statistical methods for the analysis of corrosion data for integrity assessments | |
US20200400010A1 (en) | Downhole and near wellbore reservoir state inference through automated inverse wellbore flow modeling | |
Ruiz-Díaz et al. | Two-phase oil and water flow pattern identification in vertical pipes applying long short-term memory networks | |
Reineking et al. | Computer-intensive methods in the analysis of species-habitat relationships | |
Tak et al. | Kriging models for forecasting crude unit overhead corrosion | |
Abfertiawan et al. | The application of artificial neural network model to predicting the acid mine drainage from long-term lab scale kinetic test | |
Kokatnoor et al. | Deducing water quality index (WQI) by comparative supervised machine learning regression techniques for India region | |
Alzraiee et al. | A probabilistic approach to training machine learning models using noisy data | |
Curreri et al. | Application of data distribution metrics for soft sensors in industrial scenarios | |
Pereira et al. | Using deep learning models to predict the electrical conductivity of the influent in a wastewater treatment plant | |
Ni et al. | Research on the Application of Integrated Learning Models in Oilfield Production Forecasting |
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) |