METHOD FOR CALCULATING THICKNESS OF OXIDE FILM OF MARTENSITE HEAT-RESISTANT STEEL UNDER SUPERCRITICAL HIGH-TEMPERATURE STEAM

Abstract

A method for calculating a thickness of an oxide film of a martensite heat-resistant steel under supercritical high-temperature steam is disclosed, which includes following steps: the martensite heat-resistant steel is a 9% Cr martensite heat-resistant steel; and a formula for calculating the thickness of the oxide film is

[00001] $X = A .Math. .Math. \exp (\frac{- Q}{R .Math. T}) .Math. t^{n},$

which X is the thickness of the oxide film (μm), A is a constant coefficient, Q is an activation energy (J.Math.mol.sup.−1), R is a gas constant, T is temperature (° C.), and t is time (h).

Claims

1. A method for calculating a thickness of an oxide film of a martensite heat-resistant steel under supercritical high-temperature steam, comprising: the martensite heat-resistant steel is a 9% Cr martensite heat-resistant steel; and $X = A .Math. \exp (\frac{- Q}{R .Math. T}) .Math. t^{n}$ a formula for calculating the thickness of the oxide film is wherein X is the thickness of the oxide film (μm), A is a constant coefficient, Q is an activation energy (J.Math.mol.sup.−1), R is a gas constant, T is temperature (° C.), and t is time (h).

2. The method according to claim 1, wherein the method comprises: the temperature of the high-temperature steam T is in a range of 550-700° C., the steam pressure is in a range of 23.0-35.0 MPa, and the time t is in a range of 200-20000 h.

3. The method according to claim 1, wherein n is 0.5.

4. The method according to claim 1, wherein a mathematical relationship between the activation energy Q and the time t is: Q=39659.32t.sup.0.0092.

5. The method according to claim 1, wherein the constant coefficient is 11616.83.

6. The method according to claim 1, wherein Q is 41628.07 when oxidation time is 200 h.

7. The method according to claim 1, wherein Q is 42057.09 when oxidation time is 600 h.

8. The method according to claim 1, wherein Q is 42256.60 when oxidation time is 1000 h.

9. The method according to claim 1, wherein Q is 42414.76 when oxidation time is 1500 h.

10. The method according to claim 1, wherein Q is 43426.66 when oxidation time is 2000 h.

Description

DETAILED DESCRIPTION

[0020] The metal oxidation kinetic model applied in the present disclosure includes following steps: linear velocity law, parabolic velocity law, logarithmic velocity law, and cubic velocity law. And the parabolic law may be applied when a ratio of a volume of the oxide to that of the metal is close to 1 but not more than 15%, the oxide layer may cover the metallic surface compactly meanwhile not crack due to excessive internal stress. Further, oxidation reaction is performed by mass diffusion within the oxide layer, namely, metallic cations diffuse outwards and oxygen anions diffuse inward, or cations and anions bi-directional diffuse, then new oxides are formed within the oxide layer. When an interface-reaction speed is faster than a diffusion speed, an oxide generation speed depends on the speed of the mass diffusion. The oxide film of 9% Cr heat-resistant steel is mainly divided into an inner layer and an outer layer, and growth of the oxide film disposed on the inner layer depends on O.sup.2− ions diffusing inward, and the growth of the oxide film disposed on the outer layer depends on Fe.sup.2+ ions diffusing outward, thus a growth rate of the oxide film may be described by the parabolic rate law. Researches show that the 9% Cr martensite heat-resistant steel oxidation kinetic model is basically in accordance with a parabolic model.

[0021] Assuming that the growth of the oxide film/layer is controlled only by the inward diffusion of oxygen ions, the oxygen concentration at the metal/oxide layer interface is el. The oxygen concentration at the oxide layer/water vapor interface is c.sub.0. The thickness of the oxide layer is X, accordingly a concentration gradient in the oxide layer is

[00003] $\frac{d .Math. c}{d .Math. X}$

and a diffusion coefficient of oxygen is D, the oxygen flux per unit time at the interface S is

[00004] $d .Math. n = D .Math. S .Math. \frac{d .Math. c}{d .Math. X} .Math. d .Math. t$

determined according to the Fick's first law.

[0022] Under a condition of steady-state diffusion,

[00005] $\frac{d .Math. c}{d .Math. X} = \frac{c_{0} - c_{1}}{X}$

is constant, the diffusion speed at an unit interface is

[00006] $v = \frac{d .Math. n}{s .Math. d .Math. t} = D .Math. \frac{c_{0} - c_{1}}{X} .$

Because the oxygen concentration at the interface of the metal/oxide layer is very low, the interface reaction is fast, and oxygen is not accumulated, therefore c.sub.1 approaching 0 which is negligible; when the ambient oxygen concentration (oxygen partial pressure) is constant, c.sub.0 may be regarded as a constant, the oxidation rate is inversely proportional to the oxide layer thickness. Similarly, the process of metallic cations diffusing outward may also be analyzed according to the above process, as long as the concentration difference of metal cations at the two interfaces is a constant value, the growth speed of the oxide layer is only inversely proportional to the thickness of the oxide layer:

[00007] $\begin{matrix} v = \frac{d .Math. X}{d .Math. t} = D .Math. \frac{c_{0} + c_{0}^{'}}{X} = \frac{k_{p}}{X} & (1) \end{matrix}$

[0023] wherein c.sub.0′ is the concentration difference of the metallic cations at the two interfaces.

[0024] Integrating the formula (1) to obtain X.sup.2=2k.sub.pt, or is written as

[00008] $X = {(2 .Math. k_{p})}^{\frac{1}{2}} .Math. t^{\frac{1}{2}} .$

In the formula, k.sub.p is a constant related to the diffusion coefficient, generally abided by the Arrhenius equation:

[00009] $k_{p} = k_{0} .Math. \exp (\frac{- Q}{R .Math. T}),$

wherein Q is an activation energy, R is a gas constant, k.sub.0 is a constant, k.sub.p is brought into the foregoing formula to obtain:

[00010] $X = {(2 .Math. k_{0})}^{\frac{1}{2}} .Math. \exp .Math. .Math. (\frac{- Q}{2 .Math. R .Math. T}) .Math. t^{\frac{1}{2}}$

[0025] the constant coefficient is expressed by A, further it is simplified as follows:

[00011] $\begin{matrix} X = A .Math. .Math. \exp .Math. .Math. (\frac{- Q}{R .Math. T}) .Math. t^{\frac{1}{2}} & (2) \end{matrix}$

[0026] The experimental result shows that, an expression formula of the metal high-temperature oxidation kinetics generally is X=At.sup.n, which the index n in the formula is no need to be equal with ½ which is different from formula (2), but fluctuates within a certain range, such as a cubic law when n is ⅓. The research shows that the high-temperature steam oxidation model of the 9% Cr martensite heat-resistant steel also meets same regularity that X=At.sup.n. For different steel types and working conditions, values of the coefficient A and the index n are different. Therefore, the formula of the high-temperature steam oxidation film thickness of the 9% Cr martensite heat-resistant steel can be initially defined as follows:

[00012] $\begin{matrix} X = A .Math. \exp (\frac{- Q}{R .Math. T}) .Math. t^{n} & (3) \end{matrix}$

[0027] The present disclosure collects a large amount of practical experimental data that include the temperature of 550-700° C., the steam pressure of 23.0-35.0 MPa, the thickness data of the oxide film of the 9% Cr martensitic heat-resistant steel at the oxidation time of 200-20000 h, the parameters n, Q and A in equation (3) are calculated using the above data:

[0028] The calculation method is as follows:

[0029] step 1, solving n, and taking logarithm of two sides of the formula (3) to obtain:

[00013] $\ln .Math. X = \ln .Math. A + (\frac{- Q}{R .Math. T}) + n .Math. l .Math. n .Math. t$

[0030] letting the temperature T be a constant value, then

[00014] $\ln .Math. A + (\frac{- Q}{R .Math. T})$

is a constant, denoted C, the above formula can be simplified to InX=C+nInt; Substituting the experimental data at various temperatures into the formula for linear fitting to obtain a fitting formula as follows:

[0031] when t=550° C.: ln X=0.017+0.501 ln t;

[0032] when t=600° C.: ln K=0.37+0.55 ln t;

[0033] when t=650° C.: ln X=3.14+0.255 ln t;

[0034] when t=700° C.: ln X=1.04+0.64 ln t.

[0035] It can be found that the value n is basically close to 0.5, which means the oxidation kinetics of the 9% Cr martensitic heat-resistant steel is basically in accordance with the parabolic law, so that if n is 0.5, the formula (3) is modified as follows:

[00015] $\begin{matrix} X = A .Math. .Math. \exp .Math. .Math. (\frac{- Q}{R .Math. T}) .Math. t^{0.5} & (4) \end{matrix}$

[0036] step 2, solving the activation energy Q, and taking logarithm of two sides of the formula (4) to obtain:

[00016] $\ln .Math. X = \ln (A .Math. t^{0.5}) - \frac{Q}{R .Math. T}$

[0037] letting t be a constant, then In (At.sup.0.5) is recorded as a constant C, and the above formula is simplified into:

[00017] $\begin{matrix} \ln .Math. X = C - \frac{Q}{R .Math. T} & (5) \end{matrix}$

[0038] When t is 200h, the obtained experimental data is substituted into the fitting formula obtained from the step 1 to calculate ln X values at 550° C., 600° C., 650° C. and 700° C. respectively. Then substitute ln X to formula (5), Q is calculated to be 41628.07 when t is 200h. Similarly, Q can be calculated at other time as shown in Table 1, Q may be different at different time which indicates that the activation energy of the oxidation reaction is different at different time periods. The research shows that the oxidation reaction of 9% Cr is a complex and dynamically changing process; the reaction mechanism, generated products, and compositions and structures of the oxide materials are different at different stages of oxidation. Therefore, a mathematical model is used to fit the change of the activation energy along with the time to obtain an exponential model therein the activation energy Q and the time t are highly consistent, the obtained fitting formula is as follows:

Q=(39659.32+10.29)t.sup.(0.0092±3.52E4) (6)

[0039] and substituting the formula (6) back to the formula (4) to obtain a corrected thickness formula as follows:

[00018] $\begin{matrix} X = A .Math. .Math. \exp .Math. .Math. (\frac{- 3 .Math. 9 .Math. 6 .Math. 5 .Math. 9.3 .Math. 2 .Math. t^{0.0 .Math. 0 .Math. 9 .Math. 2}}{R .Math. T}) .Math. t^{0.5} & (7) \end{matrix}$

TABLE 1 shows different activation energy at different times

TABLE-US-00001 t/h Q/J .Math. mol.sup.−1 200 41628.07 600 42057.09 1000 42256.60 1500 42414.76 20000 43426.66

[0040] Step 3, solving a constant A, substituting the experimental data into the formula (7) to perform nonlinear surface fitting, to calculate A to be 11616.83, so that the obtained fitting formula is finally as follows:

[00019] $\begin{matrix} X = 1 .Math. 1 .Math. 6 .Math. 1 .Math. 6.8 .Math. 3 .Math. .Math. \exp .Math. .Math. (\frac{- 3 .Math. 9 .Math. 6 .Math. 5 .Math. 9.3 .Math. 2 .Math. t^{0.0 .Math. 0 .Math. 9 .Math. 2}}{R .Math. T}) .Math. t^{0.5} & (8) \end{matrix}$

[0041] In the above formula, unit of the temperature T is ° C., unit of time t is h, and unit of the thickness X of the oxide film is μm, the applicable time range of the formula is 200-20000 h, the temperature range is 550-700° C., and the steam pressure range is 23.0-35.0 MPa.

Embodiment One

[0042] Comparison of the calculation methods involved in the present application with the experimental results of oxidation of T91.

[0043] In 2013, the oxidation situation of T91 steel under the conditions of 26 MPa, 600° C./650° C./700° C. is reported by Yunhai M A et al, the experimental conditions are respectively substituted into the formula provided by the instant application. The thickness of the oxide film calculated by the instant application is compared with the thickness obtained by experimental measurement, and the comparison results are shown in Table 2. It can be seen that the calculated thickness is very close to the thickness as experimentally measured.

TABLE-US-00002 TABLE 2 comparison of calculated thickness in the instant application to the thickness in actual measurement. Tem- Measured Calculated Error perature/ Time/ thickness/ thickness/ Absolute percentage/ ° C. h μm μm error/μm % 600 1100 75 80.2 5.2 6.9 650 500 106 109.8 3.8 3.6 700 1000 238 258.2 20.2 8.5

Embodiment Two

[0044] The calculation method provided by the instant application is compared with the experimental result of T/P92.

[0045] The oxidation experiment of P92 steel at 550° C. and 25 MPa, for 600h was reported by Zhongliang Zhu et al in 2013, and the thickness of the oxide film measured from the sectional SEM image was about 28 μm. By substituting the experimental conditions into the fitting formula of the instant application, the thickness of the oxide film obtained is calculated to be 28.8 μm, which is very close to the measurement result, and the error percentage is only 2.8%.

Embodiment Three

[0046] The calculation method in the instant application is applied in a real power plant environment.

[0047] Data of the steam oxide-scale thickness in the tubeline operated in the power plant which is recorded in a design guidance for preventing steam oxidation, gas corrosion and erosion on the heated surface of the pulverized coal boiler of the large power station shows that, under a condition of 600° C. and 25 MPa, the thickness of the oxidation scale of the T92 tube is 376 μm after 22981 h of operation. The above parameters are substituted into the formula of the instant application to obtain the fitting formula, the thickness of the oxide film is calculated to be 367.14 μm, and the deviation is only 2.4 percent compared with a measurement result, which means that the fitting formula performs well in practical application.

Embodiment Four

[0048] The calculation method of the instant application is applied in an actual power plant environment.

[0049] The steam pressure of an ultra-supercritical unit used in a certain power plant at abroad is about 28.4 MPa. Under a condition that a superheater tube is made of T92 material and operates for about 15000 hours at the temperature of 600° C., the thickness of the oxide scale in the tube is measured to be about 215 μm. And the thickness of the oxide scale is calculated to be about 241 μm by substituting operating parameters into the formula as disclosed in the instant application, with a deviation of 26 μm and an error percentage of 12.1%.

Embodiment Five

[0050] The calculation method as disclosed in the instant application is applied in the actual power plant environment.

[0051] The tube near outlet area of a high-temperature superheater of a 600 MW supercritical once-through boiler is made of a T91 material in a certain domestic power generation company, the steam temperature near outlet area is about 580° C., the steam pressure is about 26 MPa, and the thickness of oxide scale in the tube is about 214 μm after running for 20000 hours. Substituting above operating parameters into the formula as disclosed in the instant application, the thickness of the oxide scale is about 201.5 μm, and the error percentage is 5.8%.

[0052] The above embodiments all show that the thickness of the oxide scale/film/layer of the 9% Cr martensitic steel calculated by the aforementioned method is in good accordance with the actual measurement result, and the error percentage is controlled within 15%.

[0053] The technical solutions of the instant application are not limited to the above embodiments, and all technical solutions obtained by using equivalent substitution methods fall within the scope of the instant application.

METHOD FOR CALCULATING THICKNESS OF OXIDE FILM OF MARTENSITE HEAT-RESISTANT STEEL UNDER SUPERCRITICAL HIGH-TEMPERATURE STEAM

Inventors

Cpc classification

Classification Explorer

G01B21/085

PHYSICS

Classification Explorer

G16C20/10

PHYSICS

Classification Explorer

G01B21/08

PHYSICS

Classification Explorer

G16C60/00

PHYSICS

International classification

Classification Explorer

G16C60/00

PHYSICS

Classification Explorer

G01B21/08

PHYSICS

Classification Explorer

G16C20/10

PHYSICS

Abstract

Claims

Description