Corresponding author: Olga M. Gulina ( omgulina18@mail.ru ) Academic editor: Yury Korovin
© 2020 Valery I. Baranenko, Olga M. Gulina, Sergey A. Mironov, Nikolay L. Salnikov.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Baranenko VI, Gulina OM, Mironov SA, Salnikov NL (2020) Repeated measurements and quality of estimates in the analysis of NPP pipeline erosion-corrosion wear. Nuclear Energy and Technology 6(4): 281-287. https://doi.org/10.3897/nucet.6.60459
|
The article presents a study carried out on carbon steel pipe components subjected to erosion-corrosion wear (ECW). Based on the repeated control data, the authors present the calculated ECW characteristics, i.e., the wall-thinning and rate. It is shown that such estimates contain great uncertainty due to corrosion products deposited on the pipeline inner surface and their migration during operation. In addition, with an increase in the operating time, for example, when the lifetime is extended, the difference between the forecast and the results of control becomes larger. This means that the error in the estimates of the residual lifetime also increases. The study is based on the data of wall thickness measurements of the feedwater pipeline (273×16 mm) and steam pipeline (465×16mm) of nuclear power plants with the VVER-440 reactors, for which a sufficient number of repeated measurements were performed over a large time interval. An analysis is made of the error in estimating the pipeline wall-thinning and ECW rate using Chexal-Horowitz Flow-Accelerated Corrosion (FAC) Model (EKI-02 and EKI-03 software tools). The estimate of the ECW rate according to the above forecast model differs from the estimate according to the current control data by no more than 12.5%, since the corrosion products deposited on the pipeline inner surface wall are leveled at a large time base. When calculating the wall-thinning, due to the obvious filtering of the control data, it is possible to achieve an acceptable accuracy of estimates, i.e., about 16% without upgrading the model.
Erosion-corrosion wear (ECW), lifetime forecasting, Keller coefficient, data filtering, repeated measurements, ECW (Erosion-Corrosion Wear) rate, Chexal-Horowitz Flow-Accelerated Corrosion (FAC) Model
The results of the metal condition monitoring show a significant unevenness of erosion-corrosive wear (ECW) of the equipment and pipelines due to a large number of factors. Knowing the ECW mechanism, its dependence on the metal properties, geometric characteristics or operating conditions, it is possible to make the most effective design decisions at the stage of designing technological systems of the NPP secondary circuit, optimize the scope of operational control and reduce the number of incidents due to ECW and unscheduled shutdowns of power units as well as control costs.
The ECW rate of a component of equipment or a pipeline, calculated from the results of one-time control of wall thicknesses, contains many uncertainties (
These factors contribute to the fact that it is very difficult to estimate the error in calculating the ECW rate based on the operational control data.
There is a wide variety of methods for estimating the ECW characteristics, taking into account the above factors (
The simplest model for estimating the ECW rate WECW (
W ECW = (Snom – Smin)/∆τ, (1)
where Snom is the nominal thickness, mm; Smin is the minimum thickness according to the control data, mm; Δτ is the time interval from the date of putting the component into operation until the date of control, years; (Snom – Smin) is the thinning.
The ECW rate is used to determine the pipeline remaining lifetime τ (or its operating period until the allowable wall thickness is reached), years,
τ = (Smin – Sallow)/WECW, (2)
where Sallow is the the minimum allowable thickness (determined by (
To refine the estimates of the ECW characteristics, various methods of processing repeated measurements (
One of the ways to adequately estimate the condition of a pipeline component, as well as to solve the problem of optimizing the scope and frequency of control, is to use predictive models. One of the most common models, embedded in various software tools, is Chexal-Horowitz Flow-Accelerated Corrosion (FAC) Model (
The ECW rate in this model is determined by the following empirical dependence
W ECW = C0×F1(T)×F2(ХС)×F3(W)× F4(O2)×F5(pH)×F6(Ке)×F7(α)×F8(A), (3)
where WECW is the ECW rate, mm/year; С0 is the coefficient equal to 1 mm/year; F1(T) is the coefficient taking into account temperature; F2(ХС) is the coefficient taking into account the composition of the metal (content of chromium, copper and molybdenum); F3(W) is the coefficient determined by the fluid vlocity; F4(O2) is the coefficient taking into account the oxygen concentration; F5(pH) is the coefficient depending on рН; F6(Ке) is the coefficient taking into account the geometry of the pipeline (its standard size and Keller coefficient); F7(α) is the coefficient taking into account the steam moisture (for a single-phase environment F7(α) = 1); F8(А) is the coefficient taking into account the amine used (ammonia, morpholine, ethanolamine).
Thus, Chexal-Horowitz model takes into account the main parameters affecting the ECW rate. Using the software tools, it is possible to carry out forecast calculations and rank pipelines according to the degree of ECW exposure.
However, it was noticed that, with long operating periods, for example, when the lifetime is extended, the difference between the forecast and the results of control becomes larger. This means that the error in the estimates of the residual lifetime also increases. This situation requires additional research and possible modernization of the forecast model, if necessary.
The work is focused on the analysis of the wall thickness control data, their interpretation and estimation of the discrepancy with the forecast results depending on the operating period.
The aim of the study is to reduce the conservatism of the ECW forecast based on estimates of the informativeness of repeated measurements. As part of this study, calculations were made for components of feedwater pipelines and steam pipelines of NPPs with the RBMK-1000 reactors and NPPs with the VVER-440 reactors.
When analyzing the repeated control data, it is possible to reveal the influence of corrosion products deposited on the pipeline inner surface on the estimates of the ECW characteristics. The deposits of corrosion products lead to the fact that repeated measurements can show an increased or decreased wall thinning at a certain point, which contradicts the physical meaning of wear. However, the control data look exactly like this. Figures
Cr, % | Тstart, years | Тend, years | ECWstart, mm/yr | ECWfin, mm/yr | ECWav, mm/yr | Thinning, mm |
0.0674 | 0 | 30 | 0.4700 | 0.0535 | 0.0916 | 2.7489 |
0.0685 | 0 | 30 | 0.4641 | 0.0528 | 0.0905 | 2.7145 |
0.0696 | 0 | 30 | 0.4583 | 0.0521 | 0.0893 | 2.6809 |
0.0708 | 0 | 30 | 0.4528 | 0.0515 | 0.0883 | 2.6482 |
0.0719 | 0 | 30 | 0.4473 | 0.0509 | 0.0872 | 2.6162 |
0.073 | 0 | 30 | 0.4419 | 0.0503 | 0.0861 | 2.585 |
0.0741 | 0 | 30 | 0.4367 | 0.0497 | 0.0852 | 2.5545 |
0.0752 | 0 | 30 | 0.4317 | 0.0491 | 0.0842 | 2.5248 |
0.0764 | 0 | 30 | 0.4267 | 0.0485 | 0.0832 | 2.4957 |
0.0775 | 0 | 30 | 0.4218 | 0.0480 | 0.0822 | 2.4673 |
0.0786 | 0 | 30 | 0.4171 | 0.0474 | 0.0813 | 2.4395 |
It can be seen that, at practically the same average thickness, there was a significant migration of the deposited corrosion products. Estimating the ECW rate according to the usual formula (1) will lead to the fact that over time the estimated residual life will increase. To obtain adequate estimates of the ECW characteristics, it is necessary to use forecast models or refine the estimates using the methods for processing repeated measurements (
To make calculations using the forecast model (3), it is required to know the chromium content in the metal. Since this indicator is not known a priori, it was taken as a variable to adapt the model to the control data.
To ensure the conservatism and additional reliability, when calculating according to the EKI-02 software tool, we take three thinning values that are slightly higher than the actual value obtained in 2011, namely, 2.5 mm. Thus, the chromium values 0.0752, 0.0741 and 0.073% correspond to the thinning values 2.5248, 2.5545 and 2.585 mm, which are close to the control values (Tab.
The obtained values of the chromium content are used as parameters to forecast the ECW process for 2013 and 2014. The result of comparing the forecast and control data is shown in Fig.
According to the operational control data (
A decreased thinning at a certain point in time is associated with the deposited corrosion products, due to which the monotony of the wear process is visually disturbed. To adequately predict the ECW process, it is necessary to either modify (upgrade) the model or take to data filtering.
To adapt the EKI-02 and EKI-03 software tools, the chromium concentration in the pipeline metal and the Keller coefficient reflecting the component geometry, including, possibly, the condition of the inner surface, were selected as the most significant but unknown parameters. Since the chromium content in the pipeline metal is not known a priori, this indicator was taken as a variable. The adaptation of these tools is carried out by forecasting the thinning on the components under consideration and comparing the results obtained with the actual values. The chromium concentration is selected with the smallest discrepancy between the calculated and measured thinning for the feedwater pipeline section for 1984 and the steam pipeline section for 1985, respectively.
The ECW rate and thinning take the maximum values at a chromium content of 0.037%, and decrease at its higher values. Further, the chromium content equal to 0.037% is used, since the thinning according to the control data is significant.
The results of the thinning forecast at a chromium content of 0.037% for the feedwater pipeline section in the range from 3 to 28 years and the steam pipeline section in the range from 4 to 25 years are shown in Fig.
It can be seen from the graphs that the forecast results do not agree well with the control data. One of the reasons for this may be the Keller coefficient, which, during long-term operation, may turn out to be time-dependent.
Tables
Despite the variety of combinations of the Keller coefficient and Cr content, general trends in the development of the ECW process can be observed. The graph in Figure
Values of thinning and ECW rate of the feedwater pipeline and steam pipeline at NPPs with VVER-440 in different operating periods.
Feedwater, 273×16 mm | Steam pipelines, 465×16 mm | ||||||||
---|---|---|---|---|---|---|---|---|---|
Year | Тoper, years | Smin, mm | ∆Sthin, mm | WECW, mm/yr | Year | Тoper, years | Smin, mm | ∆Sthin, mm | WECW, mm/yr |
1984 | 3 | 13.3 | 2.7 | 0.9 | 1985 | 4 | 13.2 | 2.8 | 0.7 |
1988 | 7 | 13.2 | 2.8 | 0.4 | 1985 | 4 | 13.7 | 2.3 | 0.575 |
1989 | 8 | 14.1 | 1.9 | 0.238 | 1995 | 14 | 13.2 | 2.8 | 0.2 |
1995 | 14 | 12.0 | 4.0 | 0.286 | 1995 | 14 | 13.2 | 2.8 | 0.2 |
2001 | 20 | 12.3 | 3.7 | 0.185 | 2000 | 19 | 13.0 | 3.0 | 0.158 |
2005 | 24 | 12.2 | 3.8 | 0.158 | 2005 | 24 | 13.7 | 2.3 | 0.096 |
2006 | 25 | 12.6 | 3.4 | 0.136 | 2005 | 24 | 13.7 | 2.3 | 0.096 |
2007 | 26 | 13.1 | 2.9 | 0.112 | 2006 | 25 | 13.7 | 2.3 | 0.092 |
2008 | 27 | 12.0 | 4.0 | 0.148 | 2006 | 25 | 13.7 | 2.3 | 0.092 |
2009 | 28 | 11.7 | 4.3 | 0.154 |
Cr content,% | Ke | Thinning, mm |
---|---|---|
0.03 | 0.4699 | 2.801 |
0.04 | 0.5617 | 2.8021 |
0.05 | 0.6535 | 2.8029 |
0.06 | 0.7444 | 2.8001 |
0.07 | 0.8362 | 2.801 |
0.08 | 0.928 | 2.8017 |
Combinations of Cr content and Keller coefficient for feedwater pipelines.
Cr content,% | Ke | Thinning, mm |
---|---|---|
0.03 | 0.271 | 2.701 |
0.04 | 0.325 | 2.7109 |
0.05 | 0.379 | 2.718 |
0.06 | 0.433 | 2.7234 |
0.07 | 0.487 | 2.7276 |
0.08 | 0.541 | 2.7309 |
Based on the analysis performed, it is possible to state the discrepancy with the time of the predicted and control values. At the same time, the forecast model gives a more conservative estimation. To reduce the conservatism, it is necessary to test two possible hypotheses for improving the forecast model:
It should be expected that the intensity of the ECW process is influenced by a change in the roughness of the inner surface of the pipeline walls, caused by the deposition and washout of corrosion products with an increase in the operating life due to a change in the flow regime near the walls. However, the height roughness is an order of magnitude less than the height of the deposited corrosion products; therefore, it has a significantly smaller effect on the change in the geometry of the pipeline inner surface. Let us assume that the influence of deposits is represented by the change over time of the Keller coefficient as the only metric responsible for the pipeline geometry in the Chexal-Horowitz model (3).
The use of a statistical approach to assessing the thinning trend did not lead to the desired result, since the trend, for example, in the case of the steam pipeline generally tended to decrease with time; for the feedwater pipeline, the regression shape also changed to decrease with time. This suggests that the processes of deposition of corrosion products significantly noise the real picture. In this case, the forecast results exceeded the control data, and the regulation by the Keller coefficient was reduced to a decrease in this coefficient with time.
For the regression of the latter, the confidence interval was calculated by the formula
,
where ∆s,t is the confidence interval; β is the confidence probability (β = 0.95); n is the number of measurements (sample size, in our case: 10 and 8); tn,β is the quantile of Student’s t-distribution (for n = 10 and β = 0.95 tn,β = 2.262 and for n = 8 and β = 0.95 tn,β = 2.365); S2 is the sample variance:
,
where f (ti) are the values of the Keller coefficient by regression at time ti; yi are the values of the Keller coefficient obtained from the forecast at time ti; ti is the operating life of the component at the time of measurement.
Figure
In this case, the constructed approximations describe satisfactorily the ECW process of the component under consideration, especially the lower limit of the confidence interval (less than 35%). However, a large number of repeated measurements are required to ensure the accuracy of the estimates.
For the steam pipeline, the situation is similar, only the thinning decrease rate over time is higher for the steam pipeline than for the feedwater pipeline, which corresponds to the control trends, i.e., the deposits of corrosion products are more intensive.
The graph of changes in the wall thickness (thinning), for example, for the feedwater pipeline is described above (see Fig.
Thus, the thinning value for this area should be taken equal to the largest previous value. Any subsequent thinning values, greater than the previous ones, should be left unchanged.
For the feedwater pipeline, screened values and the corresponding regression were found. According to this regression equation, graphs were plotted of the correspondence of the thinning values obtained during the control, the forecast data for the values of the Keller coefficient in the range from 0.16 (the value according to the Keller table for a given component) to 0.3 (assuming that over time the Keller coefficient will only increase) and the thinning by regression of the screened values (Fig.
To select the best Keller coefficient, the sum of the squared forecast deviations (for different values of the Keller coefficient) were calculated depending on the regression of the screened control data:
,
where f (ti) is the screened thinning at time ti; yi is the thinning at a fixed value of the Keller coefficient at time ti; n is the number of measurements (sample size n = 10); ti is the operating life of the component at the time of measurement.
Table
For the steam line, similar results were obtained, only the Keller coefficient changes from 0.23 (the value of the coefficient according to the table for a given component) to 0.5 (under the assumption that the Keller coefficient will only increase over time), and the regression equation has the form y = 0,011t + 2,708, where t is the operating time.
The residual sum of the squared deviations is also minimal for the original geometry of the component. Thus, the control data filtering makes it possible, without upgrading the model, to achieve an acceptable accuracy of the estimated thinning and residual life (16%).
A comparison of the ECW rates calculated using the software tools and those determined from the operational control data on feed water pipelines of NPPs with the VVER-440 reactors showed (Fig.
On a long time base, the values of deposits are leveled in both directions.
Sums of the squared deviations for the shielded thinning and thinning at fixed values of the Keller coefficient.
Operating life, years | Screened thinning mm | Thinning, mm (Ке = 0.16) | Thinning, mm (Ке = 0.18) | Thinning, mm (Ке = 0.2) | Thinning, mm (Ке = 0.22) | Thinning, mm (Ке = 0.24) |
3 | 2.7357 | 1.5497 | 1.794 | 1.9933 | 2.1927 | 2.392 |
7 | 2.9789 | 2.3844 | 2.6825 | 2.9805 | 3.2786 | 3.5767 |
8 | 3.0397 | 2.529 | 2.8452 | 3.1613 | 3.4774 | 3.7936 |
14 | 3.4045 | 3.2798 | 3.6898 | 4.0998 | 4.5098 | 4.9198 |
20 | 3.7693 | 3.928 | 4.419 | 4.9101 | 5.4011 | 5.8921 |
24 | 4.0125 | 4.3309 | 4.8723 | 5.4136 | 5.955 | 6.4964 |
25 | 4.0733 | 4.4313 | 4.9852 | 5.5391 | 6.093 | 6.6469 |
26 | 4.1341 | 4.5298 | 5.096 | 5.6622 | 6.2284 | 6.7946 |
27 | 4.1949 | 4.6263 | 5.2046 | 5.7829 | 6.3612 | 6.9395 |
28 | 4.2557 | 4.7235 | 5.3139 | 5.9043 | 6.4948 | 7.0852 |
Sum of squared deviations | 2.85 | 6.15 | 14.04 | 26.41 | 43.26 |
The control data were analyzed for components of feed water pipelines and steam pipelines of NPPs with the RBMK-1000 reactors and NPPs with the VVER-440 reactors. Based on the data of operational control, it was shown that the estimates of the ECW characteristics contained great uncertainty due to the corrosion products deposited on the pipeline inner surface and their migration during operation. At the same time, the forecast thinning data differed significantly from the control data for long operating lives, which was caused by the influence of deposited corrosion products recorded during control together with the base metal.
To adapt the forecast models, EKI-02 and EKI-03, to control data under uncertainty, two approaches were considered: (1) upgrading the forecast model by changing the Keller coefficient and (2) filtering control data with the model remaining unchanged.
It was shown that data filtering made it possible to improve the quality of the thinning forecast according to the model under consideration without its modernization with a discrepancy between the forecast and control data of no more than 16%.
It was also shown that calculating the ECW rate based on the control data adequately reflected reality and, with an increase in the operating time, the forecast accuracy increased (the error was no more than 13%).
The results obtained make it possible to confidently recommend the use of the forecast models, EKI-02.1 and EKI-03.1, both at the design stage of NPPs and in the process of operation to make decisions on the timing and scope of control.