CN101832545A - Method for measuring temperatures of out-flowing water and discharged water of heater of turbine steam cooler - Google Patents
Method for measuring temperatures of out-flowing water and discharged water of heater of turbine steam cooler Download PDFInfo
- Publication number
- CN101832545A CN101832545A CN201010150270.2A CN201010150270A CN101832545A CN 101832545 A CN101832545 A CN 101832545A CN 201010150270 A CN201010150270 A CN 201010150270A CN 101832545 A CN101832545 A CN 101832545A
- Authority
- CN
- China
- Prior art keywords
- under
- heater
- temperature
- actual condition
- shell
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Control Of Turbines (AREA)
Abstract
The invention discloses a method for measuring the temperatures of the out-flowing water and the discharged water of a heater of a turbine steam cooler, which comprises the following steps: selecting a design condition or performance test condition of an assembling unit as reference conditions, and selecting the thermodynamic parameters of the j(th) heater under the reference conditions, wherein the thermodynamic parameters comprise a shell pressure, an extraction temperature, a shell extraction enthalpy, an out-flowing water temperature, a water inlet temperature and a power of the assembling unit; calculating the heat-transfer characteristic coefficient of a stream cooling section of the heater under the reference conditions; reading or calculating the actual shell pressure, extraction temperature, water inlet temperature of the j(th) heater and the actual power of the assembling unit from a supervisor information system (SIS) or decentralized control system (DCS) of a heat-engine plant; calculating and obtaining the heat-transfer characteristic coefficient of the heater according to the heat-transfer characteristic coefficient of the stream cooling section of the heater and the actual power of the assembling unit; and finally obtaining the out-flowing water temperature under the actual conditions by carrying out iterative computations, wherein the out-flowing water temperature is equal to the saturation temperature under the shell pressure.
Description
Technical field
The present invention relates to a kind of heater parameter measuring method of steam turbine carrying vapour cooler, relate in particular to a kind of heater water outlet and drain temperature measuring method of steam turbine carrying vapour cooler.
Background technology
Along with the continuous lifting of fired power generating unit parameter and capacity, receive publicity day by day to improve the unit operation economy by the performance of improving heat regenerative system.The measurement of bleeder heater water outlet and drain temperature has important effect for the thermally equilibrated calculating of heat regenerative system, unit performance monitoring and optimization, therefore is necessary it is carried out on-line monitoring.So far do not see the heater water outlet of carrying vapour cooler in the heat regenerative system and the report of drain temperature measuring method.
At present, in plant level supervisory information system SIS of thermal power plant (Supervisory Information System) or the scattered control system DCS of system (Distribution Control System), for the bleeder heater that has steam condenser, though be provided with water outlet and drain temperature measuring point, but reason such as abominable and repair and maintenance weakness because of its service condition, ubiquity is measured the situation of poor reliability, in addition, the general measure method of heater water side temperature also has the following disadvantages: at first, in the fired power generating unit thermal measurement system, the normal leaving water temperature that adopts the thermal resistance type sensor to monitor bleeder heater, data collecting system correspondingly need adopt the resistance value of active balancing bridge measurement sensor, measures the cost height; Secondly, water temperature changes thermal inertia greatly, and when power condition changing was big, the water temperature Response Table revealed bigger thermal inertia, thereby influences certainty of measurement; The 3rd, because on-the-spot installation site complexity is not easy to maintenance and maintenance.In case sensor fault or inefficacy often cause the wrong of measurement data or disappearance.
And calculate heater water outlet and drain temperature according to traditional heat transfer equation, need to calculate the heat transfer coefficient of heat transfer process.Need to understand the numerous structural parameters of heater in the calculating of heat transfer coefficient, for example: the area of each heat transfer segment of heater, flow process number, pipe side and shell-side structure, pipeline inner and outer diameter, pipeline material or the like.The disappearance of any heater data all can cause heat transfer coefficient to calculate, so traditional heat transfer equation is applicable to design and check calculates that and water outlet is calculated and monitored with drain temperature when being not easy to be used for unit operation or test.
The heat exchange of heater pure condensate knot section belongs to condensation heat transfer in power plant's heat regenerative system, and draw gas heating tube side feedwater and condensing of shell-side is characterized in that the shell-side heat transfer coefficient is very big, and gas keeps the saturation temperature of shell pressure correspondence constant in the process of condensation heat.And the heat exchange of steam condenser section belongs to steam-water heat exchanging, but the shell-side feedwater flow is much bigger with pipe side extraction flow relatively.The present invention is based on above-mentioned heat-transfer mechanism, defined the heat compensator conducting property coefficient of steam cooling section and pure condensate knot section, found the variable working condition response pattern of steam cooling section heat compensator conducting property coefficient, the heater water outlet of carrying vapour cooler and the measuring method of drain temperature have been proposed, this method does not have need understand structural parameters, measures fast, the reliable measuring data advantages of higher of low, the tested parameter response of cost.
Summary of the invention
Heater water outlet and the drain temperature measuring method a kind ofly calculating that model is simple, computational accuracy is high, measure the steam turbine carrying vapour cooler of the low and rapid dynamic response speed of cost have been the object of the present invention is to provide.
The present invention realizes by following technical solution:
A kind of heater water outlet of steam turbine carrying vapour cooler and drain temperature measuring method is characterized in that,
Step 1: the mid-transition point temperature t of heater steam cooling section and pure condensate knot section under the calculating benchmark operating mode
Wsj o:
Choose unit rated power design conditions or performance certification test operating mode as the benchmark operating mode, symbol adds that marking-up mother's " o " parameter represents that it is the parameter under the benchmark operating mode, and chooses the thermal parameter of j level heater under the benchmark operating mode: shell pressure p
Nj o, extraction temperature t
j o, the shell-side enthalpy h that draws gas
Nj o, leaving water temperature t
Wj o, inflow temperature t
W (j+1) oWith power of the assembling unit P
e o, and according to shell pressure p under the industrial water and steam thermodynamic properties of the IAPWS-IF97 model calculating benchmark operating mode
Nj oCorresponding saturation temperature t
Sj o, saturated air enthalpy h
Ssj oAnd saturation water enthalpy h
Swj o, by the shell pressure p under the benchmark operating mode
Nj o, extraction temperature t
j oAnd calculate shell-side under the benchmark operating mode enthalpy h that draws gas according to the industrial water and steam thermodynamic properties of IAPWS-IF97 model
Nj o,
By heater steam cooling section equation of heat balance:
Pure condensate knot section equation of heat balance:
In the formula: j is the heater numbering, is numbered respectively from high to low 1~n number according to the heater extraction pressure, and n is the positive integer greater than 1;
D
j oBe the j level heater amount of drawing gas, unit is kg/h;
h
Nj oBe j level heater shell-side extraction pressure enthalpy, unit is kJ/kg;
h
Ssj oBe the saturated air enthalpy of j level heater case lateral pressure correspondence, unit is kJ/kg;
h
Swj oBe the saturation water enthalpy of j level heater case lateral pressure correspondence, unit is kJ/kg;
D
Wj oBe j level heater feedwater flow, unit is kg/h;
C
pSpecific heat at constant pressure for feedwater is taken as definite value: 4.1868kJ (kg ℃);
t
Wj oBe the leaving water temperature of j level heater, unit is ℃;
t
W (j+1) oBe the inflow temperature of j level heater, unit is ℃;
Obtain the middle transition temperature t of steam cooling section and pure condensate knot section after the arrangement
Wsj oRelational expression with heater water inlet, leaving water temperature:
Step 2: the heat compensator conducting property coefficient that calculates heater steam cooling section under the benchmark operating mode
Wherein: subscript " SC " expression is the steam cooling section, (KF)
SC oBe the product of steam cooling section Coefficient K under the benchmark operating mode and heat transfer area F, unit is kJ/ (m
2℃ h) m
2
(D
wC
p)
SC oBe steam cooling section feedwater flow D under the benchmark operating mode
wSpecific heat at constant pressure C with feedwater
pProduct, unit is kg/hkJ/ (kg a ℃);
Obtain the heat compensator conducting property coefficient of steam cooling section under the benchmark operating mode
Step 3: the leaving water temperature t that calculates heater under the actual condition
WjWith drain temperature t
Dj:
Step 3.1: in the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, read the shell pressure p under the actual condition
Nj, extraction temperature t
j, j level heater inflow temperature t
W (j+1)With power of the assembling unit P
e, as if the shell pressure p that in the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, does not read under the actual condition
Nj, then by calculating the shell pressure p under the actual condition
Nj, calculate shell pressure p under the actual condition according to the industrial water and steam thermodynamic properties of IAPWS-IF97 model again
NjSaturation temperature t under the corresponding actual condition
Sj, saturated air enthalpy h
SsjAnd saturation water enthalpy h
Swj, as if the extraction temperature t that in the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, does not read under the actual condition
j, then by calculating the extraction temperature t under the actual condition
j, as if the j level heater inflow temperature t that in the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, does not read under the actual condition
W (j+1), then by calculating the j level heater inflow temperature t under the actual condition
W (j+1),
Shell pressure p under the described calculating actual condition
NjMethod be:
In the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, read the extraction pressure p under the actual condition
j, calculate the shell pressure p under the actual condition
Nj=p
j(1-δ p
j), δ p
jBe pipeline crushing rate, δ p
j=3%~5%;
The t of the extraction temperature under the described calculating actual condition
jMethod be:
Make under the actual condition shell-side enthalpy h that draws gas
NjBe taken as shell-side under the benchmark operating mode enthalpy h that draws gas
Nj o, shell-side enthalpy and shell-side under the benchmark operating mode enthalpy that draws gas that draws gas equates substantially under the variable working condition, again according to the industrial water and steam thermodynamic properties of IAPWS-IF97 model, by the shell-side enthalpy h that draws gas
NjWith the shell pressure p under the actual condition
NjCan calculate the extraction temperature t under the actual condition
j
J level heater inflow temperature t under the described calculating actual condition
W (j+1)Method be:
In the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, read the j+1 level heater case lateral pressure p under the actual condition
N (j+1), as if the j+1 level heater case lateral pressure p that in the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, does not read under the actual condition
N (j+1), then in the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, read the j+1 level heater extraction pressure p under the actual condition
J+1, calculate the j+1 level heater case lateral pressure p under the actual condition
N (j+1)=p
J+1(1-δ p
J+1), δ p
J+1Be the pipeline crushing rate of the j+1 level heater under the actual condition, δ p
J+1=3%~5%; Then according to the j+1 level heater case lateral pressure p under the industrial water and steam thermodynamic properties of the IAPWS-IF97 model calculating actual condition
N (j+1)Corresponding saturation temperature t
S (j+1), and deduct the end difference θ of j+1 level heater under design conditions with saturation temperature
J+1, and be j level heater inflow temperature t under the actual condition with this difference
W (j+1), i.e. t
W (j+1)=t
S (j+1)-θ
J+1,
Step 3.2: the iterative computation of leaving water temperature and step thereof:
Leaving water temperature t under the actual condition is set
WjIterative initial value, get the heater inflow temperature t under the actual condition
W (j+1)+ 15 as iteration initial value (t
Wj)
K=0, wherein subscript k is an iterations;
By leaving water temperature t
Wj(hypothesis) calculation of steam cooling section and pure condensate knot section middle transition temperature:
Can calculate middle transition temperature (t
Wsj)
kThen according to steam cooling section HEAT TRANSFER LAW, numerical experimentation with based on the identification of Model Parameters algorithm of sample, calculate heater heat compensator conducting property coefficient under the actual condition by the power of the assembling unit of benchmark operating mode respective heater steam cooling section heat compensator conducting property coefficient and actual condition, finally according to the leaving water temperature under this heat compensator conducting property coefficient and the middle transition temperature computation actual condition:
Wherein: m is the power of the assembling unit P under the actual condition
eWith the power of the assembling unit P under the benchmark operating mode
e oThe index of ratio, for high-pressure heater m=0.6, for low-pressure heater m=0.3,
If current leaving water temperature (t
Wj)
kDo not meet the condition of convergence and then leaving water temperature newly is worth substitution formula (1) continuation iteration, described iteration convergence condition is: Δ t
Wj=| (t
Wj)
K+1-(t
Wj)
k|≤0.01,
Satisfy the current leaving water temperature (t of the iterative computation condition of convergence
Wj)
K+1Leaving water temperature t as heater
WjEnd value,
Step 3.3: the calculating of drain temperature:
By a heater pure condensate knot section phase-change heat transfer mechanism, draw gas that to keep temperature in the condensation heat transfer process be saturation temperature under the shell pressure, draw drain temperature t
DjEqual the saturation temperature t under the shell pressure
Sj
The invention has the advantages that:
The present invention is based on heat-transfer mechanism and can survey parameter with operation, the heater steam cooling section of carrying vapour cooler and the heat compensator conducting property coefficient of pure condensate knot section have been defined, and a rule of utilizing the steam cooling section of newfound carrying vapour cooler bleeder heater and pure condensate knot section heat compensator conducting property coefficient random groups power to change, a kind of indirect, simple and direct, method that high accuracy is calculated based on heater water outlet of carrying vapour cooler and drain temperature in the heat compensator conducting property coefficient measuring and calculating fired power generating unit heat regenerative system proposed.This method only needs design (perhaps test) benchmark floor data, and need not understand structural parameters, and model is simple and direct; Only need can survey parameter according to operation, measuring and calculating water outlet and drain temperature can reduce the measurement cost; Because model uses dynamic response fast (as pressure) and high accuracy water and steam character model, can significantly improve tested parameter response speed and measure reliability.
1, the measuring and calculating model is simple, the computational accuracy height
The measuring and calculating model that the present invention set up, only need parameter and the extraction pressure and the unit load etc. of a spot of benchmark operating mode (design conditions or thermal test operating mode) can survey parameter on a small quantity, need not the structural parameters of heater and flow parameter (regenerative steam flow and feedwater or coagulate discharge), model is simple, it is simple and direct to calculate.
Than traditional heater variable working condition model, its computational accuracy height shows two aspects, the one, and the heat compensator conducting property coefficient can reflect the influence of load variations, improves response accuracy; The 2nd, what the key of water outlet, drain temperature computation model precision was the heat compensator conducting property coefficient calculates the precision of model with steam, and the heat compensator conducting property coefficient to calculate the precision of model all higher with steam, so guaranteed computation model precision of the present invention.
2, make full use of the relevant measurement result that can survey parameter, it is low to measure cost
The present invention utilizes the measurement result of pressure of extracted steam from turbine (the important measurement parameter of turbine system), realize the measuring and calculating of heater water outlet and drain temperature by model, only need pressure-measuring-point relevant in DCS or the SIS system, and need not special temperature point, measure cost by the shared reduction of metrical information.
3, use the measuring and calculating model, can significantly improve the dynamic responding speed of tested parameter
Utilize rule and the high-precision water and steam character model of the heat compensator conducting property coefficient of bleeder heater steam cooling section and pure condensate knot section with the power variation, the dynamic response of the measurement result of heater water outlet and drain temperature is equivalent to the dynamic responding speed of extraction pressure and unit load, thereby has improved the dynamic responding speed of heater water outlet and drain temperature results of measuring.
4, improved the measurement reliability of tested parameter
Extraction pressure that uses in the model and unit load are the important monitoring parameters of steam turbine, often take the redundant arrangement of measuring point to improve its reliability with measures such as making things convenient for repair and maintenance, adopt the measuring and calculating model can be, thereby improved heater water outlet and drain temperature measuring reliability with the measurement certainty equivalence of heater water outlet and drain temperature in extraction pressure and unit load measuring reliability.
Description of drawings
Fig. 1 is the principled thermal system figure of the surface heater of carrying vapour cooler
Fig. 2 is heater heat transfer process T-F (temperature-structure) figure
Fig. 3 is a calculation flow chart of the present invention
The specific embodiment
A kind of heater water outlet of steam turbine carrying vapour cooler and drain temperature measuring method is characterized in that,
Step 1: the mid-transition point temperature t of heater steam cooling section and pure condensate knot section under the calculating benchmark operating mode
Wsj o:
Choose unit rated power design conditions or performance certification test operating mode as the benchmark operating mode, symbol adds that marking-up mother's " o " parameter represents that it is the parameter under the benchmark operating mode, and chooses the thermal parameter of j level heater under the benchmark operating mode: shell pressure p
Nj o, extraction temperature t
j o, the shell-side enthalpy h that draws gas
Nj o, leaving water temperature t
Wj o, inflow temperature t
W (j+1) oWith power of the assembling unit P
e o, and according to shell pressure p under the industrial water and steam thermodynamic properties of the IAPWS-IF97 model calculating benchmark operating mode
Nj oCorresponding saturation temperature t
Sj o, saturated air enthalpy h
Ssj oAnd saturation water enthalpy h
Swj o, by the shell pressure p under the benchmark operating mode
Nj o, extraction temperature t
j oAnd calculate shell-side under the benchmark operating mode enthalpy h that draws gas according to the industrial water and steam thermodynamic properties of IAPWS-IF97 model
Nj o,
By heater steam cooling section equation of heat balance:
Pure condensate knot section equation of heat balance:
In the formula: j is the heater numbering, is numbered respectively from high to low 1~n number according to the heater extraction pressure, and n is the positive integer greater than 1;
D
j oBe the j level heater amount of drawing gas, unit is kg/h;
h
Nj oBe j level heater shell-side extraction pressure enthalpy, unit is kJ/kg;
h
Ssj oBe the saturated air enthalpy of j level heater case lateral pressure correspondence, unit is kJ/kg;
h
Swj oBe the saturation water enthalpy of j level heater case lateral pressure correspondence, unit is kJ/kg;
D
Wj oBe j level heater feedwater flow, unit is kg/h;
C
pSpecific heat at constant pressure for feedwater is taken as definite value: 4.1868kJ/ (kg ℃);
t
Wj oBe the leaving water temperature of j level heater, unit is ℃;
t
W (j+1) oBe the inflow temperature of j level heater, unit is ℃;
Obtain the middle transition temperature t of steam cooling section and pure condensate knot section after the arrangement
Wsj oRelational expression with heater water inlet, leaving water temperature:
Step 2: the heat compensator conducting property coefficient that calculates heater steam cooling section under the benchmark operating mode
By steam cooling section heat transfer equation under the benchmark operating mode:
Wherein: subscript " SC " expression is the steam cooling section, (KF)
SC oBe the product of steam cooling section Coefficient K under the benchmark operating mode and heat transfer area F, unit is kJ/ (m
2℃ h) m
2
(D
wC
p)
SC oBe steam cooling section feedwater flow D under the benchmark operating mode
wSpecific heat at constant pressure C with feedwater
pProduct, unit is kg/hkJ (kg a ℃);
Steam cooling section heat transfer temperature difference under the benchmark operating mode
Obtain the heat compensator conducting property coefficient of steam cooling section under the benchmark operating mode
Step 3: the leaving water temperature t that calculates heater under the actual condition
WjWith drain temperature t
Dj:
Step 3.1: in the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, read the shell pressure p under the actual condition
Nj, extraction temperature t
j, j level heater inflow temperature t
W (j+1)With power of the assembling unit P
e, as if the shell pressure p that in the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, does not read under the actual condition
Nj, then by calculating the shell pressure p under the actual condition
Nj, calculate shell pressure p under the actual condition according to the industrial water and steam thermodynamic properties of IAPWS-IF97 model again
NjSaturation temperature t under the corresponding actual condition
Sj, saturated air enthalpy h
SsjAnd saturation water enthalpy h
Swj, as if the extraction temperature t that in the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, does not read under the actual condition
j, then by calculating the extraction temperature t under the actual condition
j, as if the j level heater inflow temperature t that in the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, does not read under the actual condition
W (j+1), then by calculating the j level heater inflow temperature t under the actual condition
W (j+1),
Shell pressure p under the described calculating actual condition
NjMethod be:
In the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, read the extraction pressure p under the actual condition
j, calculate the shell pressure p under the actual condition
Nj=p
j(1-δ p
j), δ p
jBe pipeline crushing rate, δ p
j=3%~5%;
The t of the extraction temperature under the described calculating actual condition
jMethod be:
Make under the actual condition shell-side enthalpy h that draws gas
NjBe taken as shell-side under the benchmark operating mode enthalpy h that draws gas
Nj o, shell-side enthalpy and shell-side under the benchmark operating mode enthalpy that draws gas that draws gas equates substantially under the variable working condition, again according to the industrial water and steam thermodynamic properties of LAPWS-IF97 model, by the shell-side enthalpy h that draws gas
NjWith the shell pressure p under the actual condition
NjCan calculate the extraction temperature t under the actual condition
j
J level heater inflow temperature t under the described calculating actual condition
W (j+1)Method be:
In the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, read the j+1 level heater case lateral pressure p under the actual condition
N (j+1), as if the j+1 level heater case lateral pressure p that in the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, does not read under the actual condition
N (j+1), then in the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, read the j+1 level heater extraction pressure p under the actual condition
J+1, calculate the j+1 level heater case lateral pressure p under the actual condition
N (j+1)=p
J+1(1-δ p
J+1), δ p
J+1Be the pipeline crushing rate of the j+1 level heater under the actual condition, δ p
J+1=3%~5%; Then according to the j+1 level heater case lateral pressure p under the industrial water and steam thermodynamic properties of the IAPWS-IF97 model calculating actual condition
N (j+1)Corresponding saturation temperature t
S (j+1), and deduct the end difference θ of j+1 level heater under design conditions with saturation temperature
J+1, and be j level heater inflow temperature t under the actual condition with this difference
W (j+1), i.e. t
W (j+1)=t
S (j+1)-θ
J+1,
Step 3.2: the iterative computation of leaving water temperature and step thereof:
Leaving water temperature t under the actual condition is set
WjIterative initial value, get the heater inflow temperature t under the actual condition
W (j+1)+ 15 as iteration initial value (t
Wj)
K=0, wherein subscript k is an iterations;
By leaving water temperature t
Wj(hypothesis) calculation of steam cooling section and pure condensate knot section middle transition temperature:
Can calculate middle transition temperature (t
Wsj)
lThen according to steam cooling section HEAT TRANSFER LAW, numerical experimentation with based on the identification of Model Parameters algorithm of sample, calculate heater heat compensator conducting property coefficient under the actual condition by the power of the assembling unit of benchmark operating mode respective heater steam cooling section heat compensator conducting property coefficient and actual condition, finally according to the leaving water temperature under this heat compensator conducting property coefficient and the middle transition temperature computation actual condition:
Wherein: m is the power of the assembling unit P under the actual condition
eWith the power of the assembling unit P under the benchmark operating mode
e oThe index of ratio, for high-pressure heater m=0.6, for low-pressure heater m=0.3,
If current leaving water temperature (t
Wj)
kDo not meet the condition of convergence and then leaving water temperature newly is worth substitution formula (1) continuation iteration, described iteration convergence condition is: Δ t
Wj=| (t
Wj)
K+1-(t
Wj)
k|≤0.01,
Satisfy the current leaving water temperature (t of the iterative computation condition of convergence
Wj)
K+1Leaving water temperature t as heater
WjEnd value,
Step 3.3: the calculating of drain temperature:
By a heater pure condensate knot section phase-change heat transfer mechanism, draw gas that to keep temperature in the condensation heat transfer process be saturation temperature under the shell pressure, draw drain temperature t
DjEqual the saturation temperature t under the shell pressure
Sj
With certain 300MW unit is example, realizes the water outlet of carrying vapour cooler heater in the Steam Turbine Regenerative System and the measuring and calculating of drain temperature.The #1 heater of this unit is the high-pressure heater of carrying vapour cooler.
The detailed calculated step is as follows:
(1), the mid-transition point temperature t of heater steam cooling section and pure condensate knot section under the calculating benchmark operating mode
Ws1 o:
The design conditions of choosing unit rated power are the benchmark operating mode, and the shell pressure p of #1 heater is arranged according to design parameter
N1 oBe 5.712MPa, the shell-side enthalpy h that draws gas
N1 oBe 3136.3kJ/kg, leaving water temperature t
W1 oBe 272.1 ℃, inflow temperature t
W2 oIt is 242.5 ℃.Calculate the saturation temperature t of shell pressure correspondence according to the industrial water and steam thermodynamic properties of IAPWS-IF97 model
S1 oBe 272.4 ℃, saturated air enthalpy h
Ss1 oBe 2787.6kJ/kg, saturation water enthalpy h
Sw1 oBe 1197.3kJ/kg.
(2), calculate the heat compensator conducting property coefficient of heater steam cooling section under the benchmark operating mode
Calculate heater pure condensate knot section heat compensator conducting property coefficient:
(3), calculate the leaving water temperature t of heater under the actual condition
W1With drain temperature t
D1:
When level pressure 75% load, from the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, get shell pressure p
N1Be 4.28MPa, the shell-side enthalpy h that draws gas
N1Be taken as the analog value 3136.3kJ/kg under the benchmark operating mode, heater inflow temperature t
W2It is 226.2 ℃.Calculate shell pressure p according to the industrial water and steam thermodynamic properties of IAPWS-IF97 model
N1Corresponding saturation temperature t
S1Be 249.55 ℃, saturated air enthalpy h
Ss5Be 2799.4kJ/kg, saturation water enthalpy h
Sw5Be 1107.2kJ/kg.
Heater inflow temperature t is set
W2+ 15=241.2 ℃ is leaving water temperature t
W1Iterative initial value.
By the calculation process iterative computation of following formula according to Fig. 3,
Through 9 iteration, finally calculate the leaving water temperature t of heater
W1Be 257.57 ℃, with measured value 254.2 relative errors be :-1.325%.
By pure condensate knot section phase-change heat transfer mechanism, drain temperature t
D1Equal saturation temperature t
S1, be 249.55 ℃, with the relative error of 249.55 ℃ of measured values be 0.000%.
Claims (1)
1. the heater water outlet of a steam turbine carrying vapour cooler and drain temperature measuring method is characterized in that,
Step 1: the mid-transition point temperature t of heater steam cooling section and pure condensate knot section under the calculating benchmark operating mode
Wsj o:
Choose unit rated power design conditions or performance certification test operating mode as the benchmark operating mode, symbol adds that marking-up mother's " o " parameter represents that it is the parameter under the benchmark operating mode, and chooses the thermal parameter of j level heater under the benchmark operating mode: shell pressure p
Nj o, extraction temperature t
j o, the shell-side enthalpy h that draws gas
Nj o, leaving water temperature t
Wj o, inflow temperature t
W (j+1) oWith power of the assembling unit P
e o, and according to shell pressure p under the industrial water and steam thermodynamic properties of the IAPWS-IF97 model calculating benchmark operating mode
Nj oCorresponding saturation temperature t
Sj o, saturated air enthalpy h
Ssj oAnd saturation water enthalpy h
Swj o, by the shell pressure p under the benchmark operating mode
Nj o, extraction temperature t
j oAnd calculate shell-side under the benchmark operating mode enthalpy h that draws gas according to the industrial water and steam thermodynamic properties of IAPWS-IF97 model
Nj o,
By heater steam cooling section equation of heat balance:
Pure condensate knot section equation of heat balance:
In the formula: j is the heater numbering, is numbered respectively from high to low 1~n number according to the heater extraction pressure, and n is the positive integer greater than 1;
D
j oBe the j level heater amount of drawing gas, unit is kg/h;
h
Nj oBe j level heater shell-side extraction pressure enthalpy, unit is kJ/kg;
h
Ssj oBe the saturated air enthalpy of j level heater case lateral pressure correspondence, unit is kJ/kg;
h
Swj oBe the saturation water enthalpy of j level heater case lateral pressure correspondence, unit is kJ/kg;
D
Wj oBe j level heater feedwater flow, unit is kg/h;
C
pSpecific heat at constant pressure for feedwater is taken as definite value: 4.1868kJ/ (kg ℃);
t
Wj oBe the leaving water temperature of j level heater, unit is ℃;
t
W (j+1) oBe the inflow temperature of j level heater, unit is ℃;
Obtain the middle transition temperature t of steam cooling section and pure condensate knot section after the arrangement
Wsj oRelational expression with heater water inlet, leaving water temperature:
Step 2: the heat compensator conducting property coefficient that calculates heater steam cooling section under the benchmark operating mode
Wherein: subscript " SC " expression is the steam cooling section, (KF)
SC oBe the product of steam cooling section Coefficient K under the benchmark operating mode and heat transfer area F, unit is kJ/ (m
2℃ h) m
2
(D
wC
p)
SC oBe steam cooling section feedwater flow D under the benchmark operating mode
wSpecific heat at constant pressure C with feedwater
pProduct, unit is kg/hkJ (kg a ℃);
Obtain the heat compensator conducting property coefficient of steam cooling section under the benchmark operating mode
Step 3: the leaving water temperature t that calculates heater under the actual condition
WjWith drain temperature t
Dj:
Step 3.1: in the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, read the shell pressure p under the actual condition
Nj, extraction temperature t
j, j level heater inflow temperature t
W (j+1)With power of the assembling unit P
e, as if the shell pressure p that in the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, does not read under the actual condition
Nj, then by calculating the shell pressure p under the actual condition
Nj, calculate shell pressure p under the actual condition according to the industrial water and steam thermodynamic properties of IAPWS-IF97 model again
NjSaturation temperature t under the corresponding actual condition
Sj, saturated air enthalpy h
SsjAnd saturation water enthalpy h
Swj, as if the extraction temperature t that in the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, does not read under the actual condition
j, then by calculating the extraction temperature t under the actual condition
j, as if the j level heater inflow temperature t that in the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, does not read under the actual condition
W (j+1), then by calculating the j level heater inflow temperature t under the actual condition
W (j+1),
Shell pressure p under the described calculating actual condition
NjMethod be:
In the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, read the extraction pressure p under the actual condition
j, calculate the shell pressure p under the actual condition
Nj=p
j(1-δ p
j), δ p
jBe pipeline crushing rate, δ p
j=3%~5%;
The t of the extraction temperature under the described calculating actual condition
jMethod be:
Make under the actual condition shell-side enthalpy h that draws gas
NjBe taken as shell-side under the benchmark operating mode enthalpy h that draws gas
Nj o, shell-side enthalpy and shell-side under the benchmark operating mode enthalpy that draws gas that draws gas equates substantially under the variable working condition, again according to the industrial water and steam thermodynamic properties of IAPWS-IF97 model, by the shell-side enthalpy h that draws gas
NjWith the shell pressure p under the actual condition
NjCan calculate the extraction temperature t under the actual condition
j
J level heater inflow temperature t under the described calculating actual condition
W (j+1)Method be:
In the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, read the j+1 level heater case lateral pressure p under the actual condition
N (j+1), as if the j+1 level heater case lateral pressure p that in the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, does not read under the actual condition
N (j+1), then in the database of plant level supervisory information system SIS of thermal power plant or scattered control system DCS, read the j+1 level heater extraction pressure p under the actual condition
J+1, calculate the j+1 level heater case lateral pressure p under the actual condition
N (j+1)=p
J+1(1-δ p
J+1), δ p
J+1Be the pipeline crushing rate of the j+1 level heater under the actual condition, δ p
J+1=3%~5%; Then according to the j+1 level heater case lateral pressure p under the industrial water and steam thermodynamic properties of the IAPWS-IF97 model calculating actual condition
N (j+1)Corresponding saturation temperature t
S (j+1), and deduct the end difference θ of j+1 level heater under design conditions with saturation temperature
J+1, and be j level heater inflow temperature t under the actual condition with this difference
W (j+1), i.e. t
W (j+1)=t
S (j+1)-θ
J+1,
Step 3.2: the iterative computation of leaving water temperature and step thereof:
Leaving water temperature t under the actual condition is set
WjIterative initial value, get the heater inflow temperature t under the actual condition
W (j+1)+ 15 as iteration initial value (t
Wj)
K=0, wherein subscript k is an iterations;
By leaving water temperature t
Wj(hypothesis) calculation of steam cooling section and pure condensate knot section middle transition temperature:
Can calculate middle transition temperature (t
Wsj)
lThen according to steam cooling section HEAT TRANSFER LAW, numerical experimentation with based on the identification of Model Parameters algorithm of sample, calculate heater heat compensator conducting property coefficient under the actual condition by the power of the assembling unit of benchmark operating mode respective heater steam cooling section heat compensator conducting property coefficient and actual condition, finally according to the leaving water temperature under this heat compensator conducting property coefficient and the middle transition temperature computation actual condition:
Wherein: m is the power of the assembling unit P under the actual condition
eWith the power of the assembling unit P under the benchmark operating mode
e oThe index of ratio, for high-pressure heater m=0.6, for low-pressure heater m=0.3,
If current leaving water temperature (t
Wj)
kDo not meet the condition of convergence and then leaving water temperature newly is worth substitution formula (1) continuation iteration, described iteration convergence condition is: Δ t
Wj=| (t
Wj)
K+1-(t
Wj)
k|≤0.01,
Satisfy the current leaving water temperature (t of the iterative computation condition of convergence
Wj)
K+1Leaving water temperature t as heater
WjEnd value,
Step 3.3: the calculating of drain temperature:
By a heater pure condensate knot section phase-change heat transfer mechanism, draw gas that to keep temperature in the condensation heat transfer process be saturation temperature under the shell pressure, draw drain temperature t
DjEqual the saturation temperature t under the shell pressure
Sj
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101502702A CN101832545B (en) | 2010-04-16 | 2010-04-16 | Method for measuring temperatures of out-flowing water and discharged water of heater of turbine steam cooler |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101502702A CN101832545B (en) | 2010-04-16 | 2010-04-16 | Method for measuring temperatures of out-flowing water and discharged water of heater of turbine steam cooler |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101832545A true CN101832545A (en) | 2010-09-15 |
CN101832545B CN101832545B (en) | 2011-08-03 |
Family
ID=42716706
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101502702A Expired - Fee Related CN101832545B (en) | 2010-04-16 | 2010-04-16 | Method for measuring temperatures of out-flowing water and discharged water of heater of turbine steam cooler |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101832545B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103267539A (en) * | 2013-04-23 | 2013-08-28 | 东南大学 | Method for measuring upper terminal difference and lower terminal difference of horizontal-type three-section feed water heater |
CN103335538A (en) * | 2013-07-04 | 2013-10-02 | 国家电网公司 | Method for calculating power station condenser pressure and heat transmission terminal temperature difference |
CN104462840A (en) * | 2014-12-19 | 2015-03-25 | 东南大学 | Heat regeneration work ratio and heat regeneration gain ratio measurement method for non-reheat unit with water draining cold source losses |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4827429A (en) * | 1987-06-16 | 1989-05-02 | Westinghouse Electric Corp. | Turbine impulse chamber temperature determination method and apparatus |
CN1952590A (en) * | 2006-11-10 | 2007-04-25 | 东北电力大学 | On-line monitoring method for variation of through-flow gap of steam turbine |
-
2010
- 2010-04-16 CN CN2010101502702A patent/CN101832545B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4827429A (en) * | 1987-06-16 | 1989-05-02 | Westinghouse Electric Corp. | Turbine impulse chamber temperature determination method and apparatus |
CN1952590A (en) * | 2006-11-10 | 2007-04-25 | 东北电力大学 | On-line monitoring method for variation of through-flow gap of steam turbine |
Non-Patent Citations (1)
Title |
---|
《电力设备》 20040630 钱瑾,王培红,李琳 水和水蒸气性质通用算法研究及其应用 40-42 1 第5卷, 第6期 2 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103267539A (en) * | 2013-04-23 | 2013-08-28 | 东南大学 | Method for measuring upper terminal difference and lower terminal difference of horizontal-type three-section feed water heater |
CN103267539B (en) * | 2013-04-23 | 2015-06-17 | 东南大学 | Method for measuring upper terminal difference and lower terminal difference of horizontal-type three-section feed water heater |
CN103335538A (en) * | 2013-07-04 | 2013-10-02 | 国家电网公司 | Method for calculating power station condenser pressure and heat transmission terminal temperature difference |
CN103335538B (en) * | 2013-07-04 | 2016-01-20 | 国家电网公司 | The computational methods of a kind of condenser of power station pressure and terminal temperature difference |
CN104462840A (en) * | 2014-12-19 | 2015-03-25 | 东南大学 | Heat regeneration work ratio and heat regeneration gain ratio measurement method for non-reheat unit with water draining cold source losses |
CN104462840B (en) * | 2014-12-19 | 2018-05-18 | 东南大学 | Backheat work done ratio and backheat ratio of profit increase assay method when non-reheat unit has hydrophobic cold source energy |
Also Published As
Publication number | Publication date |
---|---|
CN101832545B (en) | 2011-08-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101825502B (en) | Effluent and drain temperature measurement and calculation method of heater with drain cooler on steam turbine | |
CN106249728B (en) | A kind of thermal power generation unit on-line performance monitoring method based on characteristics of components | |
CN100437015C (en) | On-line monitoring method for variation of through-flow gap of steam turbine | |
CN107201921B (en) | Steam turbine heat consumption rate online monitoring system and measuring method | |
CN102004460A (en) | Online monitoring method for fouling degree of flow passage of steam turbine | |
CN107577859A (en) | A kind of Condenser Fouling degree on-line monitoring method | |
CN111079302A (en) | Low-pressure cylinder efficiency measuring and calculating system and method | |
CN102338568A (en) | Online monitoring system and method for performance of condenser in power plant based on cleanness coefficient index | |
CN109388844B (en) | Correction calculation method for energy-saving effect of low-pressure economizer | |
CN101825503B (en) | Effluent and drain temperature measurement and calculation method of steam turbine surface-type heater | |
CN107543733A (en) | A kind of condenser duty on-line monitoring system and measuring method | |
Li et al. | On-line fouling monitoring model of condenser in coal-fired power plants | |
Álvarez-Fernández et al. | Thermal analysis of closed feedwater heaters in nuclear power plants | |
CN108446465B (en) | Method for measuring and calculating steam quantity for thermal power plant on line through working medium decomposition | |
CN105956370A (en) | Method for coal-fired power plant to calculate coal consumption and other economic indicators in real time | |
JP2006010229A (en) | Boiler deterioration diagnosing method, device, system and recording medium for recording program | |
WO2024061196A1 (en) | Intelligent monitoring method and system for nuclear power station steam generator | |
CN101832545B (en) | Method for measuring temperatures of out-flowing water and discharged water of heater of turbine steam cooler | |
CN104615857A (en) | Method for determining heat loads of condenser of condensing steam turbine | |
CN103063436A (en) | Turbine heat consumption rate index calculation result judgment method | |
CN109613429A (en) | A kind of PWR steam generator model time constant test macro and method | |
CN101852658B (en) | Method for measuring and calculating temperature of effluent and drain from heater with steam cooler and drain cooler of steam engine | |
CN102589824B (en) | Air leakage detection method of condenser | |
CN103267539B (en) | Method for measuring upper terminal difference and lower terminal difference of horizontal-type three-section feed water heater | |
CN104459542A (en) | Heat rate measurement method for reheating regeneration combined cycle unit |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110803 |