CN109063890B - Thermal load distribution method based on thermal power plant whole-plant peak shaving capacity maximization - Google Patents
Thermal load distribution method based on thermal power plant whole-plant peak shaving capacity maximization Download PDFInfo
- Publication number
- CN109063890B CN109063890B CN201810644679.6A CN201810644679A CN109063890B CN 109063890 B CN109063890 B CN 109063890B CN 201810644679 A CN201810644679 A CN 201810644679A CN 109063890 B CN109063890 B CN 109063890B
- Authority
- CN
- China
- Prior art keywords
- steam extraction
- unit
- power
- heat supply
- plant
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 20
- 238000000605 extraction Methods 0.000 claims abstract description 144
- 238000005457 optimization Methods 0.000 claims abstract description 15
- 238000004364 calculation method Methods 0.000 claims abstract description 13
- 238000010586 diagram Methods 0.000 claims description 12
- 230000014509 gene expression Effects 0.000 claims description 12
- 239000008235 industrial water Substances 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 238000013178 mathematical model Methods 0.000 abstract description 2
- 238000010438 heat treatment Methods 0.000 description 6
- 238000010248 power generation Methods 0.000 description 4
- 238000005485 electric heating Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000003303 reheating Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009833 condensation Methods 0.000 description 1
- 230000005494 condensation Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000002209 hydrophobic effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01K—STEAM ENGINE PLANTS; STEAM ACCUMULATORS; ENGINE PLANTS NOT OTHERWISE PROVIDED FOR; ENGINES USING SPECIAL WORKING FLUIDS OR CYCLES
- F01K17/00—Using steam or condensate extracted or exhausted from steam engine plant
- F01K17/02—Using steam or condensate extracted or exhausted from steam engine plant for heating purposes, e.g. industrial, domestic
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P80/00—Climate change mitigation technologies for sector-wide applications
- Y02P80/10—Efficient use of energy, e.g. using compressed air or pressurized fluid as energy carrier
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- General Engineering & Computer Science (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Geometry (AREA)
- Chemical & Material Sciences (AREA)
- Primary Health Care (AREA)
- Public Health (AREA)
- Mechanical Engineering (AREA)
- Combustion & Propulsion (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Engine Equipment That Uses Special Cycles (AREA)
Abstract
The invention discloses a thermal load distribution method based on the maximization of the peak shaving capacity of the whole thermal power plant. The invention starts from the research of thermoelectric relationship differences among different thermoelectric units of a thermal power plant, establishes a mathematical model of the steam extraction amount of the units and the upper limit and the lower limit of electric power, further provides a target function of the upper limit and the lower limit of total power of the whole plant, performs optimization calculation by using an fmincon function of an MATLAB tool on the premise of ensuring that the heat supply amount is not changed, and provides an optimal heat load distribution scheme, so that the total peak regulation capacity of the whole plant is maximized. The research result of the invention can optimize the operation condition of the thermal power plant on line, improve the peak regulation upper limit of the whole plant and reduce the peak regulation lower limit of the whole plant by reasonably distributing the steam extraction amount of each unit.
Description
Technical Field
The invention relates to the technical field of thermal power generation, in particular to a thermal load distribution method based on the maximization of the whole plant peak shaving capacity of a thermal power plant.
Background
In recent years, with the increase of environmental protection pressure, the power generation of new energy resources is developed by the nation, wherein the new energy resources comprise conventional hydropower, wind power, solar power generation and the like. However, the new energy unit in a part of areas develops too fast, and due to the particularity of the new energy unit in power generation, especially the inverse peak regulation characteristic of the wind turbine generator, severe challenges are brought to peak regulation of a power grid, and the new energy consumption problem is increasingly prominent. The characteristic of peak-reversal regulation of the wind turbine generator is represented by the characteristics of output of the wind turbine generator which fluctuates continuously in a certain range in the daytime and gradually rises at night, and output of the wind turbine generator which is highest in the winter and lowest in the summer, so that the peak-reversal regulation difficulty of the wind turbine generator is increased to a considerable extent when the wind turbine generator is on line, and therefore the wind curtailment of the wind turbine generator is serious, the national curtailment rate in 2015 is as high as 15%, the wind curtailment rate in 2012016 is more 21%, and the wind curtailment rate in 2017 is improved but is 12%. Meanwhile, the number of utilization hours of the thermal power generating unit is gradually reduced, the number of utilization hours of thermal power in 2015 is 4329h, the number of utilization hours of thermal power in 2016 is 4165h, and about 4000h is estimated in 2017, so that the thermal power generating unit basically operates in a peak regulation state; on the other hand, the proportion of heat supply units in northern areas is higher and higher, the thermoelectric units as main heat sources need to meet heat load preferentially, operating constraints of working conditions of 'fixing electricity by heat' exist, peak regulation capacity in heating seasons is severely restricted, even peak regulation cannot be achieved, and peak regulation difficulty is aggravated by 'wind-heat conflict'.
The research on the peak regulation capacity of the thermoelectric units in the past is mostly concentrated on the research on the electric heating characteristics of the units, the research on the peak regulation capacity of the whole plant is less, the running characteristics and the peak regulation capacities of different units are different for the thermoelectric plant covering various types of thermoelectric units, the maximum peak regulation capacity of the whole plant can be obtained by reasonably distributing heat loads under the condition that the heat supply load of the whole plant is certain, the greater the individual difference of the thermoelectric relationship among the units is, and the more obvious the effect of the heat load distribution on improving the whole peak regulation performance of the whole plant is. The heat load distribution of the whole plant must fully consider the difference of the peak shaving capacities of the different types of units, and determine the heat load distribution scheme among the different types of units, so that the total peak shaving capacity of the whole plant units is the maximum.
Disclosure of Invention
In order to solve the problems, the invention provides a plant heat load distribution method to maximize the peak regulation capacity of the whole plant by considering the peak regulation capacity difference among the heat supply units from the electric heating characteristics of the heat supply units.
The invention adopts the following technical scheme: a heat load distribution method based on the maximization of the peak regulation capacity of the whole thermal power plant is characterized by comprising the following steps of:
step 1: recording the steam extraction amount, the steam extraction pressure, the steam extraction temperature and the steam extraction drainage temperature of each unit under the current operation condition;
step 2: according to the steam extraction pressure, the steam extraction temperature and the steam extraction drainage temperature recorded in the step 1, calculating by using an industrial water and steam thermal property model to obtain steam extraction enthalpy and steam extraction drainage enthalpy of each unit;
and step 3: calculating the heat supply of each unit according to the steam extraction amount of each unit in the step 1 and the steam extraction enthalpy and the steam extraction drainage enthalpy in the step 2, and further obtaining the current total heat supply of the whole plant;
and 4, step 4: collecting heat supply working condition graphs of all units by taking a whole plant as a unit, and drawing a thermoelectric relation curve by utilizing the heat supply working condition graphs of all the units; namely, a curve of the upper power limit of each unit along with the change of the steam extraction amount and a curve of the lower power limit of each unit along with the change of the steam extraction amount;
wherein, the heat supply working condition diagram is a coupling relation diagram of the steam inlet quantity, the steam extraction quantity and the active power; each steam extraction amount has a corresponding maximum power and a corresponding minimum power, and different steam extraction amounts and corresponding maximum and minimum power values are respectively input into EXCEL to generate a thermoelectric relation curve graph;
and 5: respectively adding trend lines to the curve of the upper power limit of each unit along with the steam extraction amount and the curve of the lower power limit of each unit along with the steam extraction amount, and fitting a formula to obtain a mathematical relation formula of the upper power limit of each unit, the lower power limit and the steam extraction amount;
step 6: further accumulating to obtain a mathematical relation expression of the upper limit, the lower limit and the steam extraction quantity of the total power of the whole plant;
and 7: and (4) according to the mathematical relation among the upper limit and the lower limit of the total power of the whole plant, the steam extraction amount of each unit and the set constraint conditions, taking the mathematical relation among the upper limit and the lower limit of the total power of the whole plant and the steam extraction amount of each unit as an objective function, and solving by using an fmincon function in an MATLAB optimization toolbox to obtain the steam extraction amount of each unit based on the upper limit of the peak regulation of the whole plant, the steam extraction amount based on the lower limit of the peak regulation of the whole plant, and the corresponding maximum value and the minimum value of the upper limit and the lower limit of the total power of the whole plant, so that the heat load distribution with the maximum.
The constraint conditions are as follows: and 3, the total heat supply of each unit is equal to the current total heat supply of the whole plant calculated in the step 3, and the steam extraction amount of each unit does not exceed the maximum steam extraction amount allowed by the unit.
The step 2 is calculated using the IFC-97 industrial water and water vapor thermal property model loaded in EXCEL.
In step 3, the current heat supply amount calculation formula of each unit is as follows:
Qn′＝Gn′(hn′-hsn′)
in the formula: qn' is the current heat supply of the nth unit, Gn' is the current heat supply and steam extraction quantity of the nth unit, hn' is the current heat supply extraction enthalpy of the nth unit, hsn' is the current heat supply steam extraction and drainage enthalpy of the nth unit.
Qa＝Q1′+Q2′+…+Qn′
In the formula: qaIs the current total heat supply of the whole plant, Q1′～Qn' are the current heat supply amounts of the 1 st to nth machine sets respectively.
In step 5, the calculation formulas of the upper power limit and the lower power limit of each unit are respectively as follows:
PMAXn＝fmax(Gn)
in the formula: PMAXnIs the upper power limit of the nth unit, fmax(Gn) Is a unitary function of the heat supply steam extraction quantity and is obtained by EXCEL fitting;
PMINn＝fmin(Gn)
in the formula: PMINnIs the lower power limit, f, of the nth unitmax(Gn) Is a unitary function of the heat supply steam extraction quantity and is obtained by EXCEL fitting;
Gnthe extraction steam quantity of the nth unit.
In the step 6, the calculation formulas of the total power upper limit and the total power lower limit of the whole plant are respectively as follows:
PMAXa＝PMAX1+PMAX2+…+PMAXn
＝fmax(G1)+fmax(G2)+…+fmax(Gn)
in the formula: PMAXaIs the total power upper limit of the whole plant;
PMINa＝PMIN1+PMIN2+…+PMINn
＝fmin(G1)+fmin(G2)+…+fmin(Gn)
in the formula: PMINaIs the total power lower limit of the whole plant;
in step 7, the expression for solving the extreme values of the upper limit and the lower limit of the total power of the whole plant is as follows: because the fmincon function can only solve the minimum value, the problem of solving the maximum value of the upper limit of the power is converted into the problem of solving the minimum value of the negative function of the upper limit of the power, and the objective function is as follows:
min(-PMAXa)＝min(-(fmax(G1)+fmax(G2)+…+fmax(Gn)))
the total plant power lower limit objective function is as follows:
minPMINa＝min(fmin(G1)+fmin(G2)+…+fmin(Gn))
the constraint condition expression for solving the extreme values of the upper limit and the lower limit of the total power of the whole plant is as follows:
Qa＝G1(h1-hs1)+G2(h2-hs2)+…+Gn(hn-hsn)
0≤G1≤G1max
0≤G2≤G2max
0≤Gn≤Gn max
in the formula: qaThe current total heat supply of the whole plant calculated in the step 3 is a fixed value; g1,G2,…,GnRespectively corresponding steam extraction amounts of the first unit to the nth unit under the upper power limit and the lower power limit of the whole plant; g1max，G2max，Gn maxThe maximum steam extraction amounts of the first unit to the nth unit are fixed values and are obtained by checking heat supply working condition diagrams of the units; h is1,h2……hnThe heat supply steam extraction enthalpies of the first unit to the nth unit are respectively; h iss1、hs2……hsnThe heat supply, steam extraction and drainage enthalpies of the first unit to the nth unit are respectively.
The invention has the beneficial effects that: the invention starts from the research of thermoelectric relationship differences among different thermoelectric units of a thermal power plant, establishes a mathematical model of the steam extraction amount of the units and the upper limit and the lower limit of electric power, further provides a target function of the upper limit and the lower limit of total power of the whole plant, performs optimization calculation by using an fmincon function of an MATLAB tool on the premise of ensuring that the heat supply amount is not changed, and provides an optimal heat load distribution scheme, so that the total peak regulation capacity of the whole plant is maximized. The research result of the invention can optimize the operation condition of the thermal power plant on line, improve the peak regulation upper limit of the whole plant and reduce the peak regulation lower limit of the whole plant by reasonably distributing the steam extraction amount of each unit.
Through the online optimization of each thermal power plant, the peak regulation capacity of a power grid is improved, the consumption of new energy is promoted, remarkable economic benefits and social benefits are brought, and the application prospect is very wide.
Drawings
FIG. 1 is a diagram of the heating condition of a certain unit;
FIG. 2 is a graph showing the thermoelectric relationship of a unit;
FIG. 3 is a graph of a thermoelectric relationship curve fitting equation for a unit.
Detailed Description
The basic conditions of a certain power plant are as follows, the first-stage project #1 and #2 units are two subcritical, once intermediate reheating, steam extraction and steam condensation type units, the model is NC330/242-16.7/0.8/535/535, the steam extraction position is the steam exhaust of a medium pressure cylinder, the designed rated heat supply pressure is 0.8MPa, the rated steam extraction flow is 340t/h, and the maximum steam extraction flow is 400 t/h; the second stage project #3 and #4 units are two supercritical, one-time intermediate reheating, steam extraction and condensing units with the model number of C350/294-24.2/0.43/566/566, the steam extraction position is the steam exhaust of a medium pressure cylinder, the designed rated steam extraction pressure is 0.43MPa, the rated steam extraction flow is 400t/h, and the maximum steam extraction flow is 500 t/h.
The current operation working conditions of the whole power plant are that the #1, #3 and #4 units operate, the #2 unit stops, and the current steam extraction amount and peak regulation upper and lower limits of each unit and the whole plant are shown in the following table:
TABLE 1 current extraction and peak regulation limits for each unit and whole plant
The thermal load distribution method provided by the invention is used for carrying out on-line optimization on the peak regulation capacity of the whole plant on the thermal power plant, and the total peak regulation upper limit of the whole plant is improved and the total peak regulation lower limit is reduced by redistributing the steam extraction amount of each unit. The heat load distribution method comprises the following steps:
step 1: recording the steam extraction amount, the steam extraction pressure, the steam extraction temperature and the steam extraction drainage temperature of each unit under the current operation condition from a DCS (distributed control system) of the thermal power plant, wherein the specific table is shown in Table 2;
TABLE 2 steam extraction parameters of each unit
| Type (B) | Extraction pressure (MPa) | Extraction temperature (. degree.C.) | Extraction flow (t/h) | Steam extraction drainage temperature (. degree.C.) |
| #1 | 0.64 | 340.81 | 152.84 | 90.28 |
| #3 | 0.2 | 227.11 | 252.45 | 84.14 |
| #4 | 0.16 | 272.14 | 379 | 80.32 |
The steam extraction drainage temperature is the drainage temperature of the heat supply network heater, and the steam inlet pressure of the heat supply network heater is the default saturated pressure.
Step 2: according to the steam extraction pressure, the steam extraction temperature and the steam extraction drainage temperature recorded in the step 1, calculating by using IFC-97 industrial water and water vapor thermodynamic property calculation models loaded in EXCEL to obtain the steam extraction enthalpy and the steam extraction drainage enthalpy of each unit, as listed in Table 3;
TABLE 3 enthalpy of extraction for each unit
And step 3: and (3) calculating the heat supply of each unit according to the steam extraction amount of each unit in the step (1) and the steam extraction enthalpy and the steam extraction drainage enthalpy in the step (2), and further obtaining the current total heat supply of the whole plant.
The current heat supply calculation formula of each unit is as follows:
Qn′＝Gn′(hn′-hsn′)
in the formula: qn' is the current heat supply (GJ/h) of the nth unit, Gn' is the current heat supply and steam extraction amount (t/h), h of the nth unitn' is the current heat supply extraction enthalpy (kJ/kg) of the nth unit, hsn' is the current heat supply steam extraction hydrophobic enthalpy (kJ/kg) of the nth unit.
Qa＝Q1′+Q2′+…+Qn′
In the formula: qaIs the current total heat supply (GJ/h), Q of the whole plant1′～Qn' are the current heat supply amounts of the 1 st to nth machine sets respectively.
The heat supply amount of each unit and the total heat supply amount of the whole plant can be obtained according to the above calculation formula and the data listed in tables 2 and 3, as shown in table 4.
TABLE 3 Heat supply value of each unit
| Type (B) | Heat supply (GJ/h) |
| #1 | 643064.26 |
| #3 | 1009915.57 |
| #4 | 429456.91 |
| Whole plant | 2082436.73 |
And 4, step 4: and collecting heat supply working condition graphs of all the thermoelectric units by taking the whole plant as a unit, and drawing a thermoelectric relation curve by using the heat supply working condition graphs of all the thermoelectric units. The heating condition diagram is usually a coupling relationship diagram of the steam inlet amount, the steam extraction amount and the active power, and the corresponding relationship between the heating amount and the active power cannot be visually displayed, as shown in fig. 1. It can be seen from the figure that each extraction amount has a corresponding maximum and minimum power, and we can input different extraction amounts and corresponding maximum and minimum power values into EXCEL respectively to generate a thermoelectric relationship graph as shown in fig. 2, where the upper curve in the graph is a curve of the upper limit of the unit power with the extraction amount, and the lower curve is a curve of the lower limit of the unit power with the extraction amount.
And 5: respectively adding trend lines to the generated power upper limit curve and the power lower limit curve in the EXCEL and fitting a formula, wherein the power upper limit curve is generally an approximate straight line and can adopt linear fitting; the power lower limit curve is generally a broken line, polynomial fitting can be adopted, in most cases, the fourth-order polynomial fitting can meet the precision requirement, the trend line and the fitting formula are shown in figure 3, and therefore the mathematical relation among the power upper limit, the power lower limit and the steam extraction amount of each unit is obtained.
The calculation formulas of the upper power limit and the lower power limit of each unit are respectively as follows:
PMAXn＝fmax(Gn)
in the formula: PMAXnIs the upper power limit (MW), G of the nth unitnThe current heat supply steam extraction quantity (t/h) of the nth unit is obtained; f. ofmax(Gn) Is a unitary function of the heat supply steam extraction quantity and is obtained by EXCEL fitting;
PMINn＝fmin(Gn)
in the formula: PMINnIs the lower power limit, f, of the nth unitmax(Gn) Is a unitary function of the heat supply steam extraction quantity and is obtained by EXCEL fitting;
fitting a relation between the upper power limit and the heat supply steam extraction quantity of the #1 unit according to the heat supply working condition diagram of the #1 unit as follows:
PMAX1＝-0.2261G1+341.96
the relation between the lower power limit and the heat supply steam extraction amount of the #1 unit is as follows:
PMIN1＝1.47*10-8G1 4-1.1069*10-5G1 3+2.7366*10-3G1 2-0.47985G1+220.03
the #3 and #4 units are the same type, so that the relational expressions of the upper power limit, the lower power limit and the heat supply steam extraction amount, which are fitted according to the heat supply working condition diagram of the units, are the same, and the relational expressions of the upper power limit and the heat supply steam extraction amount are as follows:
PMAX3＝-0.2003G3+391.57
the lower power limit and the heating steam extraction quantity are in the following relation:
PMIN3＝1.2083*10-8G3 4+1.1009*10-5G3 3-1.9486*10-3G3 2-0.17934G3+235.44
step 6: and (5) accumulating to obtain a mathematical relation expression of the total power upper limit, the total power lower limit and the steam extraction amount of the whole plant according to the mathematical relation expressions of the power upper limit, the power lower limit and the steam extraction amount of each unit obtained in the step (5).
Wherein, the calculation formulas of the total power upper limit and the total power lower limit of the whole plant are respectively as follows:
PMAXa＝PMAX1+PMAX2+…+PMAXn
＝fmax(G1)+fmax(G2)+…+fmax(Gn)
in the formula: PMAXaIs the total power upper limit of the whole plant;
PMINa＝PMIN1+PMIN2+…+PMINn
＝fmin(G1)+fmin(G2)+…+fmin(Gn)
in the formula: PMINaIs the lower limit of the total power of the whole plant.
And 5, obtaining a fitting relation between the upper and lower limits of the total power of the whole plant and the heat supply steam extraction amount of each unit according to the fitting relation between the upper and lower limits of the power of each unit and the heat supply steam extraction amount in the step 5, wherein the relation between the upper limit of the total power of the whole plant and the heat supply steam extraction amount of each unit is as follows:
PMAXa＝PMAX1+PMAX2+PMAX3＝-0.2261G1+341.96-0.2003G3+391.57-0.2003G4+391.57
the relationship between the lower limit of total power of the whole plant and the heat supply steam extraction quantity of each unit is as follows:
PMINa＝PMIN1+PMIN3+PMIN4＝1.47*10-8G1 4-1.1069*10-5G1 3+2.7366*10-3G1 2-0.47985G1+220.03+1.2083*10-8G3 4+1.1009*10-5G3 3-1.9486*10-3G3 2-0.17934G3+
235.44+1.2083*10-8G4 4+1.1009*10-5G4 3-1.9486*10-3G4 2-0.17934G4+235.44
and 7: according to the mathematical relation of the upper limit and the lower limit of the total power of the whole plant, the steam extraction amount of each unit obtained in the step 6 and a set constraint condition: and (3) the total heat supply of each unit is equal to the current total heat supply of the whole plant calculated in the step (3), the steam extraction amount of each unit does not exceed the maximum steam extraction amount allowed by the unit, the mathematical relation between the upper limit and the lower limit of the total power of the whole plant and the steam extraction amount of each unit is taken as an objective function, and the fmincon function in the MATLAB optimization tool box is used for solving to obtain the steam extraction amount of each unit based on the peak regulation upper limit of the whole plant, the steam extraction amount based on the peak regulation lower limit of the whole plant, and the corresponding maximum value and the minimum value of the upper limit and the lower.
The expression for solving the extreme values of the upper limit and the lower limit of the total power of the whole plant is as follows:
because the fmincon function can only solve the minimum value, the problem of solving the maximum value of the upper limit of the power is converted into the problem of solving the minimum value of the negative function of the upper limit of the power, and the objective function is as follows:
min(-PMAXa)＝min(-(fmax(G1)+fmax(G2)+…+fmax(Gn)))
the total plant power lower limit objective function is as follows:
minPMINa＝min(fmin(G1)+fmin(G2)+…+fmin(Gn))
the constraint condition expression for solving the extreme values of the upper limit and the lower limit of the total power of the whole plant is as follows:
Qa＝G1(h1-hs1)+G2(h2-hs2)+…+Gn(hn-hsn)
0≤G1≤G1max
0≤G2≤G2max
0≤Gn≤Gn max
in the formula: qaThe current total heating load (GJ/h) of the whole plant calculated in the step 3 is a fixed value; g1,G2,…,GnRespectively corresponding steam extraction amount (t/h) of the first unit to the nth unit under the upper power limit and the lower power limit of the whole plant; g1max，G2max，Gn maxThe maximum extraction steam volume (t/h) from the first unit to the nth unit is a fixed value and is searched from a heat supply working condition diagram of each unit; h is1,h2……hnThe heat supply steam extraction enthalpies of the first unit to the nth unit are respectively; h iss1、hs2……hsnThe heat supply, steam extraction and drainage enthalpies of the first unit to the nth unit are respectively.
The method is used for carrying out online optimization on the thermal power plant, the steam extraction amount of each unit is redistributed, and the steam extraction amount and the peak regulation upper and lower limits of each unit and the whole plant after optimization are shown in a table 4:
table 4 optimized steam extraction and peak regulation upper and lower limits of each unit and whole plant
The ratio of the upper limit and the lower limit of the total peak regulation of the whole plant before and after optimization is shown in Table 5:
TABLE 5 Total Peak-shaving upper and lower limits comparison table before and after optimization
| Whole plant | Peak regulation lower limit (MW) | Peak regulation upper limit (MW) |
| Before optimization | 606.72 | 873.65 |
| After optimization | 512.61 | 964.28 |
As can be seen from Table 5, the method provided by the invention is used for carrying out online optimization on the thermal power plant, and after the thermal load of each unit is redistributed, the total peak regulation lower limit of the whole plant is obviously reduced, the total peak regulation upper limit is also obviously increased, and the total peak regulation capacity of the whole plant is obviously improved.
The above description is only a preferred embodiment of the present application and is not intended to limit the present application, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, improvement and the like made within the spirit and principle of the present application shall be included in the protection scope of the present application.
Claims (10)
1. A heat load distribution method based on the maximization of the peak regulation capacity of the whole thermal power plant is characterized by comprising the following steps of:
step 1: recording the steam extraction amount, the steam extraction pressure, the steam extraction temperature and the steam extraction drainage temperature of each unit under the current operation condition;
step 2: according to the steam extraction pressure, the steam extraction temperature and the steam extraction drainage temperature recorded in the step 1, calculating by using an industrial water and steam thermal property model to obtain steam extraction enthalpy and steam extraction drainage enthalpy of each unit;
and step 3: calculating the heat supply of each unit according to the steam extraction amount of each unit in the step 1 and the steam extraction enthalpy and the steam extraction drainage enthalpy in the step 2, and further obtaining the current total heat supply of the whole plant;
and 4, step 4: collecting heat supply working condition graphs of all units by taking a whole plant as a unit, and drawing a thermoelectric relation curve by utilizing the heat supply working condition graphs of all the units; the thermoelectric relation curve is a curve of the upper power limit of each unit along with the change of the steam extraction amount and a curve of the lower power limit of each unit along with the change of the steam extraction amount;
and 5: respectively adding trend lines to the curve of the upper power limit of each unit along with the steam extraction amount and the curve of the lower power limit of each unit along with the steam extraction amount, and fitting a formula to obtain a mathematical relation formula of the upper power limit of each unit, the lower power limit and the steam extraction amount;
step 6: further accumulating to obtain a mathematical relation expression of the upper limit, the lower limit and the steam extraction quantity of the total power of the whole plant;
and 7: according to the mathematical relation formulas of the upper limit and the lower limit of the total power of the whole plant, the steam extraction amount of each unit and the set constraint conditions, which are obtained in the step 6, the mathematical relation formulas of the upper limit and the lower limit of the total power of the whole plant and the steam extraction amount of each unit are used as objective functions, and the optimization tool box is used for solving, so that the steam extraction amount of each unit based on the upper limit of the peak regulation of the whole plant, the steam extraction amount based on the lower limit of the peak regulation of the whole plant, the corresponding maximum value of the upper limit and the minimum value of the lower limit of the total power of the whole plant are obtained, and;
the constraint condition is as follows: and 3, the total heat supply of each unit is equal to the current total heat supply of the whole plant calculated in the step 3, and the steam extraction amount of each unit does not exceed the maximum steam extraction amount allowed by the unit.
2. The method for heat load distribution based on peak shaving capacity maximization of the whole thermal power plant as claimed in claim 1, wherein the heat supply working condition diagram of the step 4 is a coupling relation diagram of steam intake, steam extraction and active power; each extraction steam quantity has a corresponding maximum and minimum power, and different extraction steam quantities and corresponding maximum and minimum power values are respectively input into EXCEL to generate a thermoelectric relation curve graph.
3. The method for heat load distribution based on peak shaving capacity maximization of a whole thermal power plant as claimed in claim 1, wherein the step 2 is calculated by using IFC-97 industrial water and water vapor thermal property models loaded in EXCEL.
4. The method for distributing heat load based on peak shaving capacity maximization of the whole thermal power plant as claimed in claim 1, wherein in the step 3, the current heat supply calculation formula of each unit is as follows:
Qn′＝Gn′(hn′-hsn′)
in the formula: qn' is the current heat supply of the nth unit, Gn' is the current heat supply and steam extraction quantity of the nth unit, hn' is the current heat supply extraction enthalpy of the nth unit, hsn' is the current heat supply steam extraction and drainage enthalpy of the nth unit.
5. The method as claimed in claim 4, wherein in the step 3, the calculation formula of the total current heat supply capacity of the whole thermal power plant is as follows:
Qa＝Q1′+Q2′+…+Qn′
in the formula: qaIs the current total heat supply of the whole plant, Q1′～Qn' are the current heat supply amounts of the 1 st to nth machine sets respectively.
6. The method according to claim 1, wherein in the step 5, the calculation formulas of the upper power limit, the lower power limit and the extraction steam volume of each unit are as follows:
PMAXn＝fmax(Gn)
in the formula: PMAXnIs the upper power limit of the nth unit, fmax(Gn) Is a unitary function of the heat supply steam extraction quantity and is obtained by EXCEL fitting;
PMINn＝fmin(Gn)
in the formula: PMINnIs the lower power limit, f, of the nth unitmax(Gn) Is a unitary function of the heat supply steam extraction quantity and is obtained by EXCEL fitting.
7. The method according to claim 1, wherein in step 6, the total power upper limit and the total power lower limit of the plant are calculated by the following formulas:
PMAXa＝PMAX1+PMAX2+…+PMAXn
＝fmax(G1)+fmax(G2)+…+fmax(Gn)
in the formula: PMAXaIs the total power upper limit of the whole plant;
PMINa＝PMIN1+PMIN2+…+PMINn
＝fmin(G1)+fmin(G2)+…+fmin(Gn)
in the formula: PMINaIs the lower limit of the total power of the whole plant.
8. The method for heat load distribution based on peak shaving capacity maximization of the whole thermal power plant as claimed in claim 1, wherein in the step 7, fmincon function in MATLAB optimization toolbox is used for solving.
9. The method according to claim 8, wherein in the step 7, the expressions for solving the extreme values of the upper limit and the lower limit of the total power of the plant are as follows: because the fmincon function can only solve the minimum value, the problem of solving the maximum value of the upper limit of the power is converted into the problem of solving the minimum value of the negative function of the upper limit of the power, and the objective function is as follows:
min(-PMAXa)＝min(-(fmax(G1)+fmax(G2)+…+fmax(Gn)))
the total plant power lower limit objective function is as follows:
minPMINa＝min(fmin(G1)+fmin(G2)+…+fmin(Gn))。
10. the method for distributing the heat load based on the maximization of the peak shaving capacity of the whole thermal power plant as claimed in claims 1 to 9, wherein the constraint condition expression is as follows:
Qa＝G1(h1-hs1)+G2(h2-hs2)+…+Gn(hn-hsn)
0≤G1≤G1max
0≤G2≤G2max
0≤Gn≤Gnmax
in the formula: qaThe current total heat supply of the whole plant calculated in the step 3 is a fixed value; g1，G2，...，GnRespectively corresponding steam extraction amounts of the first unit to the nth unit under the upper power limit and the lower power limit of the whole plant; g1max，G2max，GnmaxThe maximum steam extraction amounts of the first unit to the nth unit are fixed values and are obtained by checking heat supply working condition diagrams of the units; h is1，h2......hnThe heat supply and steam extraction enthalpies of the first unit to the nth unit are respectively set; h iss1、hs2......hsnThe first to the nth machine sets respectively supply heat, extract steam and drain enthalpy.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810644679.6A CN109063890B (en) | 2018-06-21 | 2018-06-21 | Thermal load distribution method based on thermal power plant whole-plant peak shaving capacity maximization |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810644679.6A CN109063890B (en) | 2018-06-21 | 2018-06-21 | Thermal load distribution method based on thermal power plant whole-plant peak shaving capacity maximization |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109063890A CN109063890A (en) | 2018-12-21 |
| CN109063890B true CN109063890B (en) | 2020-12-15 |
Family
ID=64821369
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810644679.6A Active CN109063890B (en) | 2018-06-21 | 2018-06-21 | Thermal load distribution method based on thermal power plant whole-plant peak shaving capacity maximization |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109063890B (en) |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110007595B (en) * | 2019-03-29 | 2022-04-22 | 常州英集动力科技有限公司 | Real-time load optimization scheduling method, unit model, unit and system for heat supply system |
| CN110566295B (en) * | 2019-07-27 | 2021-11-26 | 华电电力科学研究院有限公司 | Double-unit coupling peak shaving method and device based on power plant electricity and heat load cooperative scheduling |
| CN110717622A (en) * | 2019-09-19 | 2020-01-21 | 威立雅（哈尔滨）热电有限公司 | Optimization control method for energy management and control comprehensive indexes |
| CN110580663B (en) * | 2019-09-24 | 2023-05-12 | 国网辽宁省电力有限公司沈阳供电公司 | Calculation method for upper limit of power generation output of cogeneration unit |
| CN111047463B (en) * | 2019-12-02 | 2023-06-23 | 国网河北省电力有限公司电力科学研究院 | Peak regulating capability evaluation method for heat supply unit after heat supply modification by adopting low-temperature waste heat pump |
| CN111475951B (en) * | 2020-04-09 | 2023-06-02 | 国网山东省电力公司电力科学研究院 | Working condition analysis method of thermoelectric unit |
| CN111522365B (en) * | 2020-04-30 | 2023-04-18 | 国网山东省电力公司电力科学研究院 | Method and system for regulating and controlling steam extraction flow in multi-unit combined heat supply mode |
| CN112053027B (en) * | 2020-07-17 | 2024-05-24 | 国网山东省电力公司电力科学研究院 | Online verification method and system for minimum operation mode of thermal power plant |
| CN113033103B (en) * | 2021-03-30 | 2023-04-21 | 吉林松花江热电有限公司 | Method for determining heat consumption curve of turbine unit containing two sections of extraction steam |
| CN113537795B (en) * | 2021-07-22 | 2024-08-02 | 国网山东省电力公司电力科学研究院 | Analysis method and system for flexibility adjustment space of thermal power plant |
| CN113882917B (en) * | 2021-08-16 | 2024-08-13 | 华能国际电力股份有限公司大连电厂 | Deep peak regulation method for high back pressure unit |
| CN114046183B (en) * | 2021-09-23 | 2024-08-13 | 华能国际电力股份有限公司大连电厂 | Determination method of high back pressure operation marginal condition |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103530504A (en) * | 2013-09-27 | 2014-01-22 | 广东电网公司电力科学研究院 | System and method for calculating feasible operation range of BHKW under condition of ordering power by heat |
| CN105046064A (en) * | 2015-07-01 | 2015-11-11 | 国网天津市电力公司 | Calculation method for electric load adjustable range of heat and power cogeneration unit in heating period |
| CN105512800A (en) * | 2015-11-27 | 2016-04-20 | 国家电网公司 | Method for determining peak adjustment scheduling of heat supply unit according to mode of ordering power by heat |
| CN108182206A (en) * | 2017-12-14 | 2018-06-19 | 国电南瑞科技股份有限公司 | A kind of thermodynamic conditions thermodynamic property method for monitoring and analyzing |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103296679B (en) * | 2013-05-20 | 2016-08-17 | 国家电网公司 | The medium-term and long-term long-term wind power run that optimizes of power system is exerted oneself model modelling approach |
-
2018
- 2018-06-21 CN CN201810644679.6A patent/CN109063890B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103530504A (en) * | 2013-09-27 | 2014-01-22 | 广东电网公司电力科学研究院 | System and method for calculating feasible operation range of BHKW under condition of ordering power by heat |
| CN105046064A (en) * | 2015-07-01 | 2015-11-11 | 国网天津市电力公司 | Calculation method for electric load adjustable range of heat and power cogeneration unit in heating period |
| CN105512800A (en) * | 2015-11-27 | 2016-04-20 | 国家电网公司 | Method for determining peak adjustment scheduling of heat supply unit according to mode of ordering power by heat |
| CN108182206A (en) * | 2017-12-14 | 2018-06-19 | 国电南瑞科技股份有限公司 | A kind of thermodynamic conditions thermodynamic property method for monitoring and analyzing |
Non-Patent Citations (1)
| Title |
|---|
| 热电机组电功率限值计算方法及其在节能调度中的应用;许政强;《中国优秀硕士学位论文全文数据库(工程科技Ⅱ辑)》;20170715(第07期);第2章，第3章第3.3节 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109063890A (en) | 2018-12-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109063890B (en) | Thermal load distribution method based on thermal power plant whole-plant peak shaving capacity maximization | |
| Yang et al. | Design and simulation of gas turbine-based CCHP combined with solar and compressed air energy storage in a hotel building | |
| CN112633560B (en) | Power station optimal scheduling method containing coal-fired cogeneration unit | |
| Wang et al. | Peak shaving and heat supply flexibility of thermal power plants | |
| CN111159624B (en) | Method for calculating heat supply coal consumption rate of new steam and steam extraction combined heat supply unit | |
| Zhang et al. | Optimal clean heating mode of the integrated electricity and heat energy system considering the comprehensive energy-carbon price | |
| Wang et al. | Combined heat and power plants integrated with steam turbine renovations: Optimal dispatch for maximizing the consumption of renewable energy | |
| Ding et al. | Simulation study on a novel solar aided combined heat and power system for heat-power decoupling | |
| CN111799772A (en) | Electric heating system optimal scheduling method considering unit deep peak shaving | |
| CN106712033A (en) | Wind curtailment absorption method in thermal power plant | |
| CN109038589B (en) | Multi-provincial power grid coordinated operation production simulation method | |
| CN113283121A (en) | Flow and capacity design method and system for molten salt heat storage industrial steam supply system | |
| Guédez et al. | Techno-economic comparative analysis of innovative combined cycle power plant layouts integrated with heat pumps and thermal energy storage | |
| Wu et al. | Day-ahead optimal dispatch with CHP and wind turbines based on room temperature control | |
| CN204476487U (en) | Secondary turbine system is used in intermediate extraction heating heating | |
| Huan et al. | Retracted: Research on Optimized Operation of Electrothermal Combined System for Enhancing Wind Power Consumption | |
| CN109931581B (en) | System combining oxygen-enriched combustion with auxiliary peak regulation and frequency modulation equipment of boiler | |
| CN116105222A (en) | Heat load overall distribution method for improving flexibility of cogeneration power station | |
| CN206094020U (en) | Improve device of boiler feed water temperature | |
| CN105391085B (en) | Consider the nuclear power wind-powered electricity generation cooperation analysis method of nuclear power generating sets optimized maintenance | |
| Fang et al. | Optimal operation strategy considering wind power accommodation in heating district | |
| CN112465228A (en) | User-side comprehensive energy system optimal configuration method | |
| CN105244879A (en) | Nuclear power water pumping energy-storage power station combined operation optimization method considering nuclear power unit maintenance | |
| CN204552948U (en) | The gas handling system of Combined cycle gas-steam turbine unit | |
| CN221629852U (en) | Waste heat and electricity continuous heating system of photo-thermal power station |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |