WO2023000927A1 - Heat supply load distribution optimization method for multi-unit and multi-heat-supply-mode thermal power plant - Google Patents

Abstract

Provided in the present disclosure is a heat supply load distribution optimization method for a multi-unit and multi-heat-supply-mode thermal power plant. The method comprises: determining an optimizing objective function for the heat supply load distribution optimization method for a multi-unit and multi-heat-supply-mode thermal power plant, wherein the optimizing objective function comprises the total standard coal consumption of the whole plant; determining association characteristics of electrical load-heat supply load and electrical load-heat supply load-standard coal consumption in different heat supply modes of each thermal cogeneration unit; and according to the total standard coal consumption of the whole plant, performing plant-level operation optimization under given conditions for a heat supply load and electrical loads of different sub-units.

Description

Optimization method of heating load distribution in thermal power plants with multiple units and multiple heating modes

Cross References to Related Applications

This application is based on a Chinese patent application with application number 202110831724.0 and a filing date of July 22, 2021, and claims the priority of this Chinese patent application. The entire content of this Chinese patent application is hereby incorporated by reference into this application.

technical field

The disclosure belongs to the field of cogeneration and heat supply of coal-fired generating sets, and in particular relates to an optimization method for heat supply load distribution of thermal power plants with multiple units and multiple heat supply modes.

Background technique

For cogeneration enterprises, the rational distribution of power generation and heat supply is of great significance to reduce the operating cost of cogeneration, while cogeneration enterprises with multiple units and multiple heating modes are faced with plant-level operation optimization problems, among which heat supply The distribution of load among multiple units and multiple heating modes has a great influence on the energy saving and consumption reduction indicators of thermal power plants.

At present, relevant professionals have conducted some research on the operation optimization of cogeneration units, but there is no in-depth study on the optimal distribution of heat load among heating units with 2 or more units and 2 or more heating modes. The literature "Operation Optimization of Cogeneration Units with Heat Storage Tanks" introduces the problem of operation optimization of cogeneration units with heat storage tanks. The research object is two units in thermal power plants. The low-pressure connecting pipe is drilled to extract steam, and the research object is too single. The master's thesis "Study on Optimizing the Distribution of Heat Loads during Multi-unit Combined Steam Extraction and Heat Supply" studies the heat load distribution principles of the case plant under different power generation loads, but the heat supply load of each calculation condition is fixed, and the actual heat supply time of the thermal power plant is not considered. Changes in heating load during the severe cold period, the average period, and the beginning and end of the period. The document "Analysis and Optimization of Steam Extraction-High Back Pressure Combined Heating Supply with Two-machine Joint Commissioning" studies the technology of extraction heat supply, steam extraction-high back pressure combined heat supply, and dual-machine joint commissioning extraction-high back pressure combined heat supply technology The heating capacity of the three, but the focus of the research is on the comparison of the heating capacity of the three heating technologies, and no in-depth research has been conducted on the optimal distribution of the heating load. Master thesis "Energy-saving Analysis and Operation Optimization of High Backpressure Heating Units" studied the influence of power generation load, heating load, supply and return water temperature and back pressure on the performance of high backpressure heating units. The research object was a 300MW cogeneration For the unit, the heating method is only high back pressure heating, and the research scope is relatively narrow.

Therefore, in order to guide the optimization of the heat supply load distribution of multi-unit (2 and above) and multi-heat supply modes (2 and above) cogeneration enterprises, it is urgent to propose a multi-unit and multi-heat supply model. Optimization method for heating load distribution in thermal power plants.

Contents of the invention

The purpose of this disclosure is to provide a thermal power plant heating load distribution optimization method with multiple units and multiple heating modes, so as to reduce the total consumption of standard coal at the plant level of the cogeneration unit under the given conditions of external heating load and sub-unit electrical load The working condition corresponding to the lowest value is regarded as the optimal operating condition, and the operating back pressure of the high back pressure heating unit, the heat load borne by other units and the corresponding heating mode are listed.

According to the first aspect of an embodiment of the present disclosure, a multi-unit, multi-heating mode thermal power plant heating load distribution optimization method is provided, including: determining the search for the multi-unit, multi-heating mode thermal power plant heating load distribution optimization method The optimal objective function, the optimal objective function includes the total consumption of standard coal in the whole plant; determine the correlation characteristics of electricity load-heat supply load, electricity load-heat supply load-standard coal consumption of each cogeneration unit under different heating modes; According to the total standard coal consumption of the whole plant, the plant-level operation optimization is carried out under the given conditions of heating load and electrical load of different sub-units.

In some embodiments, when determining the optimization objective function of the thermal power plant heating load distribution optimization method with multiple units and multiple heating modes, the optimization objective function includes the total amount of standard coal consumption of the entire plant, the total amount of the entire plant The calculation method of profit value M is as follows:

M=H+EC=Q _t ×h+(N _ge1 +N _ge2 +N _ge3 )×e-(B ₁ +B ₂ +B ₃ )×c

Among them, H is heat sales revenue; Q _t is total external heat supply load, h is ex-factory heat price, E is electricity sales revenue, N _ge1 , N _ge2 , N _ge3 are real-time power generation loads of three units respectively, and e is grid connection Electricity price, C is the cost of standard coal consumption, B ₁ , B ₂ , B ₃ are the real-time standard coal consumption of the three units respectively, and c is the unit price of standard coal.

In some embodiments, when determining the correlation characteristics of electric load-heat supply load, electric load-heat supply load-standard coal consumption under different heating modes of each cogeneration unit, when the unit adopts high back pressure cascade heating , the calculation method of heating load Q is as follows:

t _s =f(P _c )=-0.029×P _c ² +2.28×P _c +26.13

Among them, m _cw is the flow rate of circulating water in and out of the high back pressure heating network condenser, C _p is the constant pressure specific heat capacity of circulating water in the heating network, t ₀ and t ₁ are the heat flow in and out of the high back pressure heating network condenser respectively Network circulating water temperature, δt is the heat exchange end difference of the high back pressure heating network condenser, t _s is the condensate temperature of the exhaust steam after the heat release of the high back pressure heating network condenser, P _c is the high back pressure machine The group's running backpressure.

The heating load Q is a multivariate function of the electric load N _ge1 , the circulating water temperature t ₀ and t ₁ of the heating network circulating water entering and leaving the high back pressure heating network condenser, the circulating water flow m _cw of the heating network and the operating back pressure P _c , through The following formula characterizes:

Q＝f ₁ (N _ge1 ,P _c ,m _cw ,t ₀ )

The standard coal consumption B under the condition of dual supply of electricity and heat of the high back pressure heating unit is a binary function of the electric load N _ge1 and the operating back pressure P _c , which is characterized by the following formula:

B=F ₁ (N _ge1 ,P _c ).

In some embodiments, in determining the correlation characteristics of electricity load-heat supply load, electricity load-heat supply load-standard coal consumption under different heating modes of each cogeneration unit, when the unit adopts the middle exhaust steam extraction mode , the electrical load N _ge is determined, and the heating load is between Q between 0 and the maximum value Q _max :

0≤Q≤Q _max ＝f ₂ (N _ge )

The standard coal consumption B is a binary function of the electrical load N _ge and the heating load Q, which is characterized by the following formula:

B=F ₂ (N _ge ,Q).

In some embodiments, in determining the correlation characteristics of electricity load-heat supply load, electricity load-heat supply load-standard coal consumption under different heat supply modes of each cogeneration unit, when the unit uses low-pressure cylinders to supply heat with zero output In mode, the unit operates with heat constant power, and the heating load Q is a one-variable linear function of the electric load N _ge , which is represented by the following formula:

Q＝f ₃ (N _ge )

The standard coal consumption B is a binary function of the electrical load N _ge and the heating load Q, which is characterized by the following formula:

B=F ₃ (N _ge )

Under the condition of boiler design output D _ms0 , the unit adopts low-pressure cylinder zero output heating mode, and the operating range of electric load is:

0≦N _ge ≦N _ge,max =f ₄ (D _ms0 ).

According to the total standard coal consumption of the whole plant, the plant-level operation optimization under the given conditions of heating load and different sub-generator loads includes:

S301, input boundary parameters, the boundary parameters include heating return water temperature t ₀ , heating network circulating water flow m _cw , total heating load Q _t , generating load N _ge1 of Unit 1, generating load N _ge2 of Unit 2, and generating load of Unit 3 Unit generating load N _ge3 ;

S302, formulating an initial value of back pressure P _c0 for the operation of the high back pressure heating unit, and using the preset threshold value of the temperature rise of the high back pressure heating network condenser corresponding to the back pressure P _c0 as the back pressure reference for iterative optimization;

S303, calculate Q ₁ =f ₁ (N _ge ,P _c0 ,t ₀ ,m _cw ), where Q ₁ is the heat supply load of unit 1, N _ge is the real-time power generation load, P _c0 is the running back pressure, t ₀ is the circulating water temperature of the heating network entering the high back pressure heating network condenser, m _cw is the circulating water flow rate of the heating network entering and leaving the high back pressure heating network condenser;

S304, for the No. 2 unit, judge whether the electric load conditions of the middle exhaust steam extraction mode and the zero output mode of the low-pressure cylinder are met, if yes, go to S305, if not, go to S315;

S305, for unit No. 3, judge whether the electric load conditions of the two modes of mid-exhaust steam extraction or low-pressure cylinder zero output are met, if yes, go to S306, if not, go to S311 and S315;

S306, determine the heat load of the No. 2 unit, and judge whether the No. 2 unit has the heat load conditions of the two modes of mid-exhaust steam extraction and low-pressure cylinder zero output, if yes, enter S307, and if not, execute S315;

S307, determine the heat load of the No. 3 unit, and judge whether the No. 3 unit has the heat load conditions of the two modes of mid-exhaust steam extraction and low-pressure cylinder zero output, if yes, enter S308, and if not, execute S315;

S308. Carry out thermal load judgment, judge whether f _2-2 (N _ge2 )+f _3-3 (N _ge3 ) is greater than or equal to Q _t -Q ₁ , if yes, enter S309, if not, terminate the optimization iterative process;

S309, determine the heat load, and judge whether f _2-3 (N _ge3 )+f _3-2 (N _ge2 ) is greater than or equal to Q _t -Q ₁ , if yes, enter S310, if not, terminate the optimization iterative process;

S310, No. 2 unit and No. 3 unit have two heat supply modes: mid-exhaust steam extraction or low-pressure cylinder zero output, enter the first iterative optimization mode, and output the optimal result;

S311, judging whether the No. 2 unit has the electric load and heat load conditions of the mid-exhaust steam extraction mode or the low-pressure cylinder zero output mode for operation, and at the same time, whether the No. 3 unit has the electric load and thermal load conditions of the mid-exhaust steam extraction mode for operation condition, then enter S312, if not, then terminate the optimization iterative process;

S312, enter the second iterative optimization mode, and output the optimal result;

S313, judging whether the No. 2 unit has the electric load and thermal load conditions of the mid-exhaust steam extraction mode in operation, and at the same time, whether the No. 3 unit has the electric load and thermal load of the mid-exhaust steam extraction mode and low-pressure cylinder zero output mode. condition, if so, then enter S314, if not, then terminate the optimization iterative process;

S314, performing a third iterative optimization mode, and outputting an optimal result;

S315, judging whether the No. 2 unit and the No. 3 unit have the electric load and heat load conditions of the exhaust steam extraction mode in operation, and if so, proceed to S316, and if not, terminate the optimization iterative process;

S316. Perform a fourth iterative optimization mode, and output an optimal result.

The first iteration optimization mode includes:

S10-1-1: Establish the initial value of the back pressure P _c0 of the high back pressure heating unit operation, and take the back pressure P _c0 corresponding to the temperature rise of the high back pressure heating network condenser at 5°C as the back pressure benchmark for iterative optimization.

S10-1-2: Calculate Q ₁₀ =f ₁ (N _ge ,P _c0 ,t ₀ ,m _cw );

S10-1-3: let Q ₂₀ =f _3-2 (N _ge2 ), then Q ₃₀ =Q _t -Q ₁₀ -f _3-2 (N _ge2 );

S10-1-4: judge whether Q ₃₀ is greater than 0, if so, then enter S10-1-5; if not, then transfer to S10-1-11;

S10-1-5: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _3-2 (N _ge2 ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ), B ₀ =B ₁₀ +B ₂₀ +B ₃₀ , enter S10-1-6;

S10-1-6: Let Q ₂₁ =Q ₂₀ -(f _3-2 (N _ge2 )-f _2-2 (N _ge2 )), then Q ₃₁ =Q ₃₀ +(f _3-2 (N _ge2 )- f _2-2 (N _ge2 )), enter S10-1-7;

S10-1-7: Judging whether Q ₂₁ is less than 0, or whether Q ₃₁ is greater than or equal to f _3-3 (N _ge3 ), if yes, then the optimization iteration process is terminated; if not, then transfer to S10-1-8;

S10-1-8: Judging whether Q ₃₁ is less than or equal to f _2-3 (N _ge3 ), if yes, proceed to S10-1-9; if not, proceed to S10-1-15;

S10-1-9: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₁ =F _2-2 (N _ge2 ,Q ₂₁ ), B ₃₁ =F _2-3 (N _ge3 ,Q ₃₁ ) , B ₁ =B ₁₀ +B ₂₁ +B ₃₁ , enter S10-1-10;

S10-1-10: Determine whether B ₁ is less than or equal to B ₀ , if not, determine that the original reference working condition is still the reference working condition; if so, determine that the new working condition is the optimal working condition, set B ₁ =B _b0 , Transfer to S10-1-6;

S10-1-11: Determine whether f _2-2 (N _ge2 ) is greater than or equal to Q _t -Q ₁₀ , if yes, go to S10-1-12; if not, go to S10-1-14;

S10-1-12: set Q ₂₀ =Q _t -Q ₁₀ , then Q ₃₀ =0, enter S10-1-13;

S10-1-13: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ) , B ₀ =B ₁₀ +B ₂₀ +B ₃₀ , enter S10-1-6;

S10-1-14: Let Q ₂₀ =f _2-2 (N _ge2 ), then Q ₃₀ =Q _t -Q ₁₀ -f _2-2 (N _ge2 ), go to S10-1-13;

S10-1-15: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₁ =F _2-2 (N _ge2 ,Q ₂₁ ), B ₃₁ =F _3-3 (N _ge3 ), B ₁ =B ₁₀ +B ₂₁ +B ₃₁ , enter S10-1-10;

S10-1-16: Based on the lowest value B _b0 of the total coal consumption of the whole plant under the back pressure of P _c0 , the heat load distribution Q ₂ and Q ₃ of the No. 2 unit and No. 3 unit are obtained as the next back pressure iterative optimization The benchmark of comparison, enter S10-1-17;

S10-1-17: let P _c1 =P _c0 +1kPa, enter S10-1-18;

S10-1-18: Determine whether P _c1 is less than or equal to min(P _c,s, P _c,max ), if yes, increase the back pressure of Unit 1 by 1kPa, and enter S10-1-2; if not, Then the optimization iterative process is terminated, where P _c,s and P _c,max are the maximum back pressure value of the safe operation of the high back pressure unit and the actual maximum value that can be achieved by the back pressure adjustment means;

S10-1-19: Taking the lowest value B _b of the total coal consumption of the whole plant as the optimization target, obtain the operating back pressure P _cb of No. 1 unit, the heat load distribution Q _2b and Q _3b of No. 2 and No. 3 units;

S10-2-1: Formulate the initial value of the back pressure P _c0 of the high back pressure heating unit operation, and use the preset threshold of the temperature rise of the high back pressure heating network condenser corresponding to the back pressure P _c0 as the back pressure benchmark for iterative optimization;

S10-2-2: Calculate Q ₁₀ =f ₁ (N _ge ,P _c0 ,t ₀ ,m _cw );

S10-2-3: let Q ₃₀ =f _3-3 (N _ge3 ), then Q ₂₀ =Q _t -Q ₁₀ -f _3-3 (N _ge3 );

S10-2-4: Determine whether Q ₂₀ is greater than or equal to 0, if yes, then enter S10-2-5; if not, then proceed to S10-2-11;

S10-2-5: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _3-3 (N _ge3 ), B ₀ =B ₁₀ +B ₂₀ +B ₃₀ , enter S10-2-6;

S10-2-6: Let Q ₃₁ =Q ₃₀ -(f _3-3 (N _ge3 )-f _2-3 (N _ge3 )), then Q ₂₁ =Q ₂₀ +(f _3-3 (N _ge3 )- f _2-3 (N _ge3 )), enter S10-2-7;

S10-2-7: Determine whether Q ₃₁ is less than or equal to 0, or whether Q ₂₁ is greater than or equal to f _3-2 (N _ge2 ), if yes, terminate the optimization iterative process; if not, transfer to S10-2-8;

S10-2-8: Judging whether Q ₂₁ is less than or equal to f _2-2 (N _ge2 ), if yes, proceed to S10-2-9; if not, proceed to S10-2-15;

S10-2-9: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₁ =F _2-2 (N _ge2 ,Q ₂₁ ), B ₃₁ =F _2-3 (N _ge3 ,Q ₃₁ ) , B ₁ =B ₁₀ +B ₂₁ +B ₃₁ , enter S10-2-10;

S10-2-10: Determine whether B ₁ is less than or equal to B ₀ , if not, determine that the original reference working condition is still the reference working condition; if so, determine that the new working condition is the optimal working condition, set B ₁ =B _b0 , Transfer to S10-2-6;

S10-2-11: Determine whether f _2-3 (N _ge3 ) is greater than or equal to Q _t -Q ₁₀ , if yes, go to S10-2-12; if not, go to S10-2-14;

S10-2-12: set Q ₃₀ =Q _t -Q ₁₀ , then Q ₂₀ =0, enter S10-2-13;

S10-2-13: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ) , B ₀ =B ₁₀ +B ₂₀ +B ₃₀ , enter S10-2-6;

S10-2-14: Let Q ₃₀ =f _2-3 (N _ge3 ), then Q ₂₀ =Q _t -Q ₁₀ -f _2-3 (N _ge3 ), go to S10-2-13;

S10-2-15: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₁ =F _3-2 (N _ge2 , B ₃₁ =F _2-3 (N _ge3 ,Q ₃₁ ), B ₁ = B ₁₀ +B ₂₁ +B ₃₁ , enter S10-2-10;

S10-2-16: Based on the lowest value B _b0 of the total coal consumption of the whole plant under the back pressure of P _c0 , the heat load distribution Q ₂ and Q ₃ of the No. Compare benchmarks, enter S10-2-17;

S10-2-17: let P _c1 =P _c0 +1kPa, enter S10-2-18;

S10-2-18: Determine whether P _c1 is less than or equal to min(P _c,s, P _c,max ), if yes, increase the back pressure of Unit 1 by 1kPa, and enter S10-2-2; if not, Then the optimization iterative process is terminated, where P _c,s and P _c,max are the maximum back pressure value of the safe operation of the high back pressure unit and the actual maximum value that can be achieved by the back pressure adjustment means;

S10-2-19: Taking the lowest value B _b of the total coal consumption of the whole plant as the optimization target, obtain the operating back pressure P _cb of No. 1 unit, the heat load distribution Q _2b and Q _3b of No. 2 and No. 3 units;

The specific methods of the second iterative optimization mode include:

S12-1: Formulate the initial value of the back pressure P _c0 of the high back pressure heating unit operation, and use the preset threshold of the temperature rise of the high back pressure heating network condenser corresponding to the back pressure P _c0 as the back pressure benchmark for iterative optimization;

S12-2: Calculate Q ₁₀ =f ₁ (N _ge ,P _c0 ,t ₀ ,m _cw );

S12-3: Let Q ₂₀ =f _3-2 (N _ge2 ), then Q ₃₀ =Q _t -Q ₁₀ -f _3-2 (N _ge2 );

S12-4: judge whether Q ₃₀ is greater than or equal to 0, if yes, then enter S12-5; if not, then transfer to S12-10;

S12-5: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _3-2 (N _ge2 ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ), B ₀ =B ₁₀ +B ₂₀ +B ₃₀ ;

S12-6: Let Q ₂₁ =Q ₂₀ -(f _3-2 (N _ge2 )-f _2-2 (N _ge2 )), then Q ₃₁ =Q ₃₀ +(f _3-2 (N _ge2 )-f _{2 -2} (N _ge2 ));

S12-7: Determine whether Q ₂₁ is less than or equal to 0, or whether Q ₃₁ is greater than or equal to f _2-3 (N _ge3 ), if so, terminate the optimization iterative process; if not, enter S12-8;

S12-8: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₁ =F _2-2 (N _ge2 ,Q ₂₁ ), B ₃₁ =F _2-3 (N _ge3 ,Q ₃₁ ), B ₁ = B ₁₀ +B ₂₁ +B ₃₁ , enter S12-9;

S12-9: Determine whether B ₁ is less than or equal to B ₀ , if not, determine that the original reference working condition is still the reference working condition; if so, determine that the new working condition is the optimal working condition, set B ₁ = B _b0 , and transfer to S12-6;

S12-10: Determine whether f _2-2 (N _ge2 ) is greater than or equal to Q _t -Q ₁₀ , if yes, go to S12-11; if not, go to S12-13;

S12-11: set Q ₂₀ =Q _t -Q ₁₀ , then Q ₃₀ =0, enter S12-12;

S12-12: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ), B ₀ = B ₁₀ +B ₂₀ +B ₃₀ , enter S12-6;

S12-13: Let Q ₂₀ =f _2-2 (N _ge2 ), then Q ₃₀ =Q _t -Q ₁₀ -f _2-2 (N _ge2 ), go to S12-12;

S12-14: Based on the lowest value B _b0 of the total coal consumption of the whole plant under the back pressure of P _c0 , the heat load distribution Q ₂ and Q ₃ of the No. 2 unit and No. 3 unit are obtained, which will be used as the comparison benchmark for the next back pressure iterative optimization , enter S12-15;

S12-15: Let P _c1 =P _c0 +1kPa, go to S12-16;

S12-16: Determine whether P _c1 is less than or equal to min(P _c,s, P _c,max ), if yes, increase the back pressure of Unit 1 by 1kPa, and enter S12-2; if not, terminate the optimization Iterative process, where P _c,s and P _c,max are the highest back pressure value for safe operation of high back pressure units and the actual highest value that can be achieved by back pressure adjustment means;

S12-17: Taking the lowest value B _b of the total coal consumption of the whole plant as the optimization target, obtain the operating back pressure P _cb of No. 1 unit, and the heat load distribution Q _2b and Q _3b of No. 2 and No. 3 units.

The specific methods of the third iterative optimization mode include:

S14-1: Formulate the initial value of the back pressure P _c0 of the high back pressure heating unit operation, and use the preset threshold of the temperature rise of the high back pressure heating network condenser corresponding to the back pressure P _c0 as the back pressure benchmark for iterative optimization;

S14-2: Calculate Q ₁₀ =f ₁ (N _ge ,P _c0 ,t ₀ ,m _cw );

S14-3: Let Q ₃₀ =f _3-3 (N _ge3 ), then Q ₂₀ =Q _t -Q ₁₀ -f _3-3 (N _ge3 );

S14-4: Determine whether Q ₂₀ is greater than or equal to 0, if yes, then enter S14-5; if not, then proceed to S14-10;

S14-5: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _3-3 (N _ge3 ), B ₀ =B ₁₀ +B ₂₀ +B ₃₀ ;

S14-6: Let Q ₃₁ =Q ₃₀ -(f _3-3 (N _ge3 )-f _2-3 (N _ge3 )), then Q ₂₁ =Q ₂₀ +(f _3-3 (N _ge3 )-f _{2 -3} (N _ge3 ));

S14-7: Determine whether Q ₃₁ is less than or equal to 0, or whether Q ₂₁ is greater than or equal to f _2-2 (N _ge2 ), if so, terminate the optimization iterative process; if not, enter S14-8;

S14-8: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₁ =F _2-2 (N _ge2 ,Q ₂₁ ), B ₃₁ =F _2-3 (N _ge3 ,Q ₃₁ ), B ₁ = B ₁₀ +B ₂₁ +B ₃₁ , enter S14-9;

S14-9: Determine whether B ₁ is less than or equal to B ₀ , if not, determine that the original reference working condition is still the reference working condition; if so, determine that the new working condition is the optimal working condition, set B ₁ = B _b0 , and transfer to S14-6;

S14-10: Determine whether f _2-3 (N _ge3 ) is greater than or equal to Q _t -Q ₁₀ , if yes, go to S14-11; if not, go to S14-13;

S14-11: set Q ₃₀ =Q _t -Q ₁₀ , then Q ₂₀ =0, enter S14-12;

S14-12: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ), B ₀ = B ₁₀ +B ₂₀ +B ₃₀ , enter S14-6;

S14-13: Let Q ₃₀ =f _2-3 (N _ge3 ), then Q ₂₀ =Q _t -Q ₁₀ -f _3-3 (N _ge3 ), enter S14-12;

S14-14: Based on the lowest value B _b0 of the total coal consumption of the whole plant under the back pressure of P _c0 , the heat load distribution Q ₂ and Q ₃ of the second unit and the third unit are obtained, which will be used as the comparison benchmark for the next back pressure iterative optimization, Enter S14-15;

S14-15: let P _c1 =P _c0 +1kPa, go to S14-16;

S14-16: Determine whether P _c1 is less than or equal to min(P _c,s, P _c,max ), if yes, increase the back pressure of Unit 1 by 1kPa, and enter S14-2; if not, terminate the optimization Iterative process, where P _c,s and P _c,max are the highest back pressure value for safe operation of high back pressure units and the actual highest value that can be achieved by back pressure adjustment means;

S14-17: Taking the lowest value B _b of the total coal consumption of the whole plant as the optimization goal, obtain the operating back pressure P _cb of the No. 1 unit, and the heat load distribution Q _2b and Q _3b of the No. 2 unit and No. 3 unit.

The specific methods of the fourth iterative optimization mode include:

S16-1: Formulate the initial value of the back pressure P _c0 of the high back pressure heating unit operation, and use the preset threshold of the temperature rise of the high back pressure heating network condenser corresponding to the back pressure P _c0 as the back pressure benchmark for iterative optimization;

S16-2: Calculate Q ₁₀ =f ₁ (N _ge ,P _c0 ,t ₀ ,m _cw );

S16-3: let Q ₂₀ =f _2-2 (N _ge2 ), then Q ₃₀ =Q _t -Q ₁₀ -f _2-2 (N _ge2 );

S16-4: judge whether Q ₃₀ is greater than or equal to 0, if yes, then enter S16-5; if not, then transfer to S16-10;

S16-5: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ), B ₀ = B ₁₀ +B ₂₀ +B ₃₀ ;

S16-6: Let Q ₂₁ =Q ₂₀ -0.1×Q ₂₀ , then Q ₃₁ =Q ₃₀ +0.1×Q ₂₀ ;

S16-7: Determine whether Q ₂₁ is equal to 0, or whether Q ₃₁ is equal to f _2-3 (N _ge3 ), if so, terminate the optimization iterative process; if not, enter S16-8;

S16-8: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₁ =F _2-2 (N _ge2 ,Q ₂₁ ), B ₃₁ =F _2-3 (N _ge3 ,Q ₃₁ ), B ₁ = B ₁₀ +B ₂₁ +B ₃₁ , enter S16-9;

S16-9: Determine whether B ₁ is less than or equal to B ₀ , if not, determine that the original reference working condition is still the reference working condition; if so, determine that the new working condition is the optimal working condition, set B ₁ = B _b0 , and transfer to S16-6;

S16-10: Set Q ₂₀ =Q _t -Q ₁₀ , Q ₃₀ =0, enter S16-11;

S16-11: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ), B ₀ = B ₁₀ +B ₂₀ +B ₃₀ , enter S16-6;

S12-12: Based on the lowest value B _b0 of the total coal consumption of the whole plant under the back pressure of P _c0 , the thermal load distribution Q ₂ and Q ₃ of the No. 2 unit and No. 3 unit are obtained, which will be used as the comparison benchmark for the next back pressure iterative optimization , enter S16-13;

S16-13: let P _c1 =P _c0 +1kPa, go to S16-14;

S16-14: Determine whether P _c1 is less than or equal to min(P _c,s, P _c,max ), if yes, increase the back pressure of Unit 1 by 1kPa, and enter S16-2; otherwise, terminate the optimization iteration process, where P _c,s and P _c,max are the maximum back pressure value for the safe operation of the high back pressure unit and the actual maximum value that can be achieved by the back pressure adjustment means;

S16-15: Taking the lowest value B _b of the total coal consumption of the whole plant as the optimization goal, obtain the operating back pressure P _cb of No. 1 unit, and the heat load distribution Q _2b and Q _3b of No. 2 and No. 3 units.

According to the second aspect of the embodiments of the present disclosure, there is provided a multi-unit, multi-heating mode thermal power plant heating load distribution optimization device, including: a first determination module configured to determine the multi-unit, multi-heating mode The optimization objective function of the thermal power plant heating load distribution optimization method, the optimization objective function includes the total standard coal consumption of the whole plant; the second determination module is configured to determine the electric load of each cogeneration unit under different heating modes -Heating load, electricity load-heating load-standard coal consumption correlation characteristics; and an optimization module configured to perform given conditions for heating load and electrical load of different sub-units according to the total standard coal consumption of the whole plant The next plant level operation optimization.

According to a third aspect of an embodiment of the present disclosure, there is provided an electronic device, including: a processor; and a memory for storing instructions executable by the processor; wherein the processor is configured to execute the instructions to The method for optimizing the heat supply load distribution of thermal power plants with multiple units and multiple heat supply modes in the first aspect of the present disclosure is realized.

According to the fourth aspect of the embodiments of the present disclosure, there is provided a computer-readable storage medium, on which a computer program is stored, which is characterized in that, when the program is executed by a processor, the multi-unit, multi-heating system of the first aspect of the present disclosure can be realized Model-based thermal power plant heating load distribution optimization method.

According to the fifth aspect of the embodiments of the present disclosure, there is provided a computer program product, including a computer program, when the computer program is executed by a processor, it realizes the multi-unit, multi-heating mode thermal power plant heating load distribution of the first aspect of the present disclosure Optimization.

The present disclosure proposes a real-time dispatching of sub-unit electric load and total heating load, as well as heating return water temperature, supply temperature, etc. The plant-level operation optimization method under multi-variable constraint conditions such as thermal circulation water flow, takes the minimum value of the total standard coal consumption at the plant level under the conditions of heating and power supply loads as the optimization objective function, and obtains the high back pressure heating unit operation back pressure, the heat load borne by other units and the corresponding heating mode. Existing traditional optimization methods that use heat rate and coal consumption as targets are only suitable for scenarios where the number of units does not exceed 2 and the heating mode is single. The calculation process is cumbersome and difficult for technical personnel in related fields to understand and accept. Compared with the traditional method, the optimization target of the present disclosure is intuitive and conforms to the production reality, greatly simplifies the optimization process, expands the number of cogeneration units in operation and the heating mode, and has broad application prospects.

Description of drawings

Figure 1 is a schematic diagram of the heating system of the whole plant;

Fig. 2 is a flow chart of a method for optimizing heating load distribution of thermal power plants with multiple units and multiple heating modes according to an embodiment of the present disclosure;

Fig. 3 is a flow chart of a method for optimizing heating load distribution of a thermal power plant with multiple units and multiple heating modes according to an embodiment of the present disclosure;

FIG. 4 is a flowchart of a first iterative optimization mode according to an embodiment of the present disclosure;

FIG. 5 is a flowchart of a first iterative optimization mode according to an embodiment of the present disclosure;

FIG. 6 is a flowchart of a second iterative optimization mode according to an embodiment of the present disclosure;

FIG. 7 is a flowchart of a third iterative optimization mode according to an embodiment of the present disclosure;

FIG. 8 is a flowchart of a fourth iterative optimization mode according to an embodiment of the present disclosure;

FIG. 9 is a block diagram of an electronic device according to an embodiment of the disclosure.

Among them, 1. High and medium pressure cylinder, 2. Low pressure cylinder, 3. High back pressure heat network condenser, 4. Heat network circulating water pump, 5. Heat network heater, 6. Heating butterfly valve.

detailed description

The present disclosure will be further described below in conjunction with the accompanying drawings.

The heat supply power plant consists of four units, numbered No. 1, No. 2, No. 3 and No. 4. Among them, No. 1 and No. 2 units have the same parameters and capacity, which is called the first-stage unit; No. 3 unit The parameters and capacity parameters of the No. 4 unit are the same, and it is called the second-stage unit; the unit parameters and capacity of the first-stage unit and the second-stage unit can be the same or different. In order to make the method proposed in the present disclosure more practical, it is agreed that the unit parameters and capacities of the first-stage unit and the second-stage unit in the implementation of the present disclosure are different.

Unit 1 has high back pressure heat supply, unit 2 has two modes of mid-exhaust steam extraction and low-pressure cylinder zero output heating, unit 3 has two modes of mid-exhaust steam extraction and low-pressure cylinder zero output heating, unit 4 Pure condensing operation, not responsible for external heating. The flow chart of the heating system of the whole plant is shown in Figure 1. After the heat supply return water is boosted by the heating network circulating water pump 4, it first enters the high back pressure heating network condenser 3, absorbs the exhaust heat of the steam turbine low-pressure cylinder 2 of the No. 1 unit, and then enters the heating network heater 5. After the peak temperature of the exhaust steam of the high and medium pressure cylinder 1 of unit No. 1 is raised, it is supplied to the outside.

The multi-unit, multi-heating mode thermal power plant heating load distribution optimization method proposed in the embodiment of the present disclosure does not involve No. 4 unit.

As shown in FIG. 2 , the multi-unit, multi-heating mode thermal power plant heating load distribution optimization method provided by an embodiment of the present disclosure includes the following steps 210 to 230 .

Step 210: Determine the optimization objective function of the thermal power plant heating load distribution optimization method with multiple units and multiple heating modes.

The production profit value M of the coal-fired thermal power plant is calculated according to formula (1).

M=H+EC=Q _t ×h+(N _ge1 +N _ge2 +N _ge3 )×e-(B ₁ +B ₂ +B ₃ )×c (1)

In the formula, H is heat sales revenue; Q _t is the total external heat supply load, MW; h is ex-factory heat price, yuan/MW; E is electricity sales revenue; N _ge1 , N _ge2 , N _ge3 are respectively The real-time power generation load of No. 2 and No. 3 units, MW; e is the on-grid electricity price, yuan/MW; C is the cost of standard coal consumption; B ₁ , B ₂ , and B ₃ are respectively The real-time standard coal consumption of the unit, t; c is the unit price of standard coal, yuan/t; B=B ₁ +B ₂ +B ₃ .

_Coal - _fired thermal power plants supply heat and power to the outside world. The total heat supply load Q _t is _dispatched by heat users in real time according to demand. The heat and electricity load scheduling authority.

Economically related data such as ex-factory heat price h, on-grid electricity price e, and standard coal unit price c are affected by market supply and demand, and generally do not change in the short term.

To sum up, the optimization objective function of the thermal power plant heating load distribution _{optimization method} with multiple units and multiple heating modes is the total standard coal consumption B of the whole _plant . , N _ge3 constraint conditions, by changing the distribution of the total heating load Q _t in each unit, the total standard coal consumption B of the whole plant is the lowest, and the highest profitability can be achieved.

Step 220: Obtain the correlation characteristics of electricity load-heat supply load and electricity load-heat supply load-standard coal consumption of each cogeneration unit under different heat supply modes.

Through the method of combining performance test and theoretical calculation, the high back pressure cascade heating of No. 1 unit, mid-exhaust steam extraction and low-pressure cylinder zero output heating of No. Electric load-heat supply load and electric load-heat supply load-standard coal consumption correlation characteristics of output heat supply.

Unit 1 provides cascade heating with high back pressure, and the heating load Q is calculated according to formula (2).

In the formula, m _cw is the circulating water flow of the heating network entering and leaving the high back pressure heating network condenser, t/h; C _p is the specific heat capacity of the circulating water of the heating network at constant pressure, kJ/kg·K; t ₀ and t ₁ are respectively The temperature of circulating water in and out of the high back pressure heat network condenser, °C; δt is the heat exchange end difference of the high back pressure heat network condenser, °C; t _s is the exhaust steam in the high back pressure heat network condenser Condensate temperature after heat release, °C; P _c is the operating back pressure of the high back pressure unit, kPa.

The heating load Q is a multivariate function of the electric load N _ge1 , the circulating water temperature t ₀ and t ₁ of the heating network circulating water entering and leaving the high back pressure heating network condenser, the circulating water flow m _cw of the heating network and the operating back pressure P _c Formula (3) representation.

Q＝f ₁ (N _ge1 ,P _c ,m _cw ,t ₀ ) (3)

The standard coal consumption B of the high back pressure heating unit under the condition of dual supply of electricity and heat is a binary function of the electric load N _ge1 and the operating back pressure P _c , which is represented by formula (4).

B＝F ₁ (N _ge1 ,P _c ) (4)

Units No. 2 and No. 3 also have two modes of mid-exhaust steam extraction and low-pressure cylinder zero output heating.

In the mid-exhaust steam extraction mode, the electric load N _ge is given, and the heating load Q is flexibly adjustable between 0 and the maximum value Q _max : 0≤Q≤Q _max =f ₂ (N _ge ), at this time the standard coal consumption The quantity B is a binary function of the electrical load N _ge and the heating load Q: B=F ₂ (N _ge ,Q).

When switching to the low-pressure cylinder zero-output heating mode, the unit operates with constant heat and electricity, and the heating load Q is a linear function of the electric load N _ge : Q=f ₃ (N _ge ), at this time the standard coal consumption B is the electric load Binary function of load N _ge and heating load Q: B=F ₃ (N _ge ). Under the condition of boiler design output D _ms0 , the unit adopts the low-pressure cylinder zero output heating mode, and the electric load operating range: 0≤N _ge ≤N _ge,max = f ₄ (D _ms0 ).

In order to distinguish Unit 2 and Unit 3, subscripts are added to the functions f ₂ , F ₂ , f ₃ , and F ₃ . The second unit is: f _2-2 , F _2-2 , f _3-2 , F _3-2 . Unit 3 is: f _2-3 , F _2-3 , f _3-3 , F _3-3 .

Step 230: With the lowest total standard coal consumption B of the whole plant, perform plant-level operation optimization under the given conditions of total heating load Q _t and sub-generator loads N _ge1 , N _ge2 , and N _ge3 .

Fig. 3 to Fig. 8 show the implementation process of the method for optimizing the heat supply load distribution of thermal power plants with multiple units and multiple heating modes in the embodiment of the present disclosure.

Fig. 3 shows a general diagram of the implementation of the optimization method for heat supply load distribution of a thermal power plant with multiple units and multiple heat supply modes in an embodiment of the present disclosure.

S1: Input boundary parameters, including heating return water temperature t ₀ , heating network circulating water flow m _cw , total heating load Q _t , generating load N _ge1 of unit 1, generating load N _ge2 of unit 2, and generating load of unit 3 Unit generating load N _ge3 ;

S2: Formulate the initial value of the back pressure P _c0 of the high back pressure heating unit operation, and take the back pressure P _c0 corresponding to the temperature rise of the high back pressure heating network condenser at 5°C as the back pressure benchmark for iterative optimization.

S3: Calculate Q ₁ =f ₁ (N _ge ,P _c0 ,t ₀ ,m _cw );

S4: For Unit No. 2, judge whether N _ge2 is less than or equal to f _4-2 (D _ms0-2 ), if so, determine that Unit No. 2 has electric loads in two modes of mid-exhaust steam extraction and low-pressure cylinder zero output. condition, and go to S5; if not, go to S15;

S5: For Unit No. 3, judge whether N _ge3 is less than or equal to f _4-3 (D _ms0-3 ), if so, determine that Unit No. 3 has electric loads in two modes of mid-exhaust steam extraction or low-pressure cylinder zero output. condition, and go to S6; if not, go to S11 and S15;

S6: Heat load determination: whether f _3-2 (N _ge2 ) is less than or equal to Q _t -Q ₁ , if yes, it is determined that the No. And go to S7; No go to S15;

S7: Heat load determination: whether f _3-3 (N _ge3 ) is less than or equal to Q _t -Q ₁ , if yes, it is determined that unit No. And go to S8; No go to S11 and S15;

S8: Carry out heat load judgment, judge whether f _2-2 (N _ge2 )+f _3-3 (N _ge3 ) is greater than or equal to Q _t -Q ₁ , if yes, go to S309, if not, stop the optimization iterative process;

S9: Carry out heat load judgment, judge whether f _2-3 (N _ge3 )+f _3-2 (N _ge2 ) is greater than or equal to Q _t -Q ₁ , if yes, enter S310, if not, terminate the optimization iterative process;

S10: Units No. 2 and No. 3 have two heat supply modes of mid-exhaust steam extraction or low-pressure cylinder zero output, enter iterative optimization mode 1-1 and 1-2, and output the optimal result;

S10-1-1: Formulate the initial value of the back pressure P _c0 of the high back pressure heating unit operation, and take the back pressure P _c0 corresponding to the temperature rise of the high back pressure heating network condenser at 5°C as the back pressure benchmark for iterative optimization;

S10-1-2: Calculate Q ₁₀ =f ₁ (N _ge ,P _c0 ,t ₀ ,m _cw );

S10-1-3: let Q ₂₀ =f _3-2 (N _ge2 ), then Q ₃₀ =Q _t -Q ₁₀ -f _3-2 (N _ge2 );

S10-1-4: judge whether Q ₃₀ is greater than 0, if so, then enter S10-1-5; if not, then transfer to S10-1-11;

S10-1-5: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _3-2 (N _ge2 ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ), B ₀ =B ₁₀ +B ₂₀ +B ₃₀ , enter S10-1-6;

S10-1-6: Let Q ₂₁ =Q ₂₀ -(f _3-2 (N _ge2 )-f _2-2 (N _ge2 )), then Q ₃₁ =Q ₃₀ +(f _3-2 (N _ge2 )- f _2-2 (N _ge2 )), enter S10-1-7;

S10-1-7: Judging whether Q ₂₁ is less than 0, or whether Q ₃₁ is greater than or equal to f _3-3 (N _ge3 ), if yes, then the optimization iteration process is terminated; if not, then transfer to S10-1-8;

S10-1-8: Judging whether Q ₃₁ is less than or equal to f _2-3 (N _ge3 ), if yes, proceed to S10-1-9; if not, proceed to S10-1-15;

S10-1-9: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₁ =F _2-2 (N _ge2 ,Q ₂₁ ), B ₃₁ =F _2-3 (N _ge3 ,Q ₃₁ ) , B ₁ =B ₁₀ +B ₂₁ +B ₃₁ , enter S10-1-10;

S10-1-10: Determine whether B ₁ is less than or equal to B ₀ , if not, determine that the original reference working condition is still the reference working condition; if so, determine that the new working condition is the optimal working condition, set B ₁ =B _b0 , Transfer to S10-1-6;

S10-1-11: Determine whether f _2-2 (N _ge2 ) is greater than or equal to Q _t -Q ₁₀ , if yes, go to S10-1-12; if not, go to S10-1-14;

S10-1-12: set Q ₂₀ =Q _t -Q ₁₀ , then Q ₃₀ =0, enter S10-1-13;

S10-1-13: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ) , B ₀ =B ₁₀ +B ₂₀ +B ₃₀ , enter S10-1-6;

S10-1-14: Let Q ₂₀ =f _2-2 (N _ge2 ), then Q ₃₀ =Q _t -Q ₁₀ -f _2-2 (N _ge2 ), go to S10-1-13;

S10-1-15: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₁ =F _2-2 (N _ge2 ,Q ₂₁ ), B ₃₁ =F _3-3 (N _ge3 ), B ₁ =B ₁₀ +B ₂₁ +B ₃₁ , enter S10-1-10;

S10-1-16: Based on the lowest value B _b0 of the total coal consumption of the whole plant under the back pressure of P _c0 , the heat load distribution Q ₂ and Q ₃ of the No. 2 unit and No. 3 unit are obtained as the next back pressure iterative optimization The benchmark of comparison, enter S10-1-17;

S10-1-17: let P _c1 =P _c0 +1kPa, enter S10-1-18;

S10-1-18: Determine whether P _c1 is less than or equal to min(P _c,s, P _c,max ), if yes, increase the back pressure of Unit 1 by 1kPa, and enter S10-1-2; if not, Then the optimization iterative process is terminated, where P _c,s and P _c,max are the maximum back pressure value of the safe operation of the high back pressure unit and the actual maximum value that can be achieved by the back pressure adjustment means;

S10-1-19: Taking the lowest value of the total coal consumption B _b of the whole plant as the optimization target, obtain the operating back pressure P _cb of No. 1 unit, and the heat load distribution Q _2b and Q _3b of No. 2 and No. 3 units.

S10-2: See Figure 5 for the implementation process of optimization mode 1-2.

S10-2-1: Formulate the initial value of the back pressure P _c0 of the high back pressure heating unit operation, and use the preset threshold of the temperature rise of the high back pressure heating network condenser corresponding to the back pressure P _c0 as the back pressure benchmark for iterative optimization;

S10-2-2: Calculate Q ₁₀ =f ₁ (N _ge ,P _c0 ,t ₀ ,m _cw );

S10-2-3: let Q ₃₀ =f _3-3 (N _ge3 ), then Q ₂₀ =Q _t -Q ₁₀ -f _3-3 (N _ge3 );

S10-2-4: Determine whether Q ₂₀ is greater than or equal to 0, if yes, then enter S10-2-5; if not, then proceed to S10-2-11;

S10-2-5: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _3-3 (N _ge3 ), B ₀ =B ₁₀ +B ₂₀ +B ₃₀ , enter S10-2-6;

S10-2-6: Let Q ₃₁ =Q ₃₀ -(f _3-3 (N _ge3 )-f _2-3 (N _ge3 )), then Q ₂₁ =Q ₂₀ +(f _3-3 (N _ge3 )- f _2-3 (N _ge3 )), enter S10-2-7;

S10-2-7: Determine whether Q ₃₁ is less than or equal to 0, or whether Q ₂₁ is greater than or equal to f _3-2 (N _ge2 ), if yes, terminate the optimization iterative process; if not, transfer to S10-2-8;

S10-2-8: Judging whether Q ₂₁ is less than or equal to f _2-2 (N _ge2 ), if yes, proceed to S10-2-9; if not, proceed to S10-2-15;

S10-2-9: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₁ =F _2-2 (N _ge2 ,Q ₂₁ ), B ₃₁ =F _2-3 (N _ge3 ,Q ₃₁ ) , B ₁ =B ₁₀ +B ₂₁ +B ₃₁ , enter S10-2-10;

S10-2-10: Determine whether B ₁ is less than or equal to B ₀ , if not, determine that the original reference working condition is still the reference working condition; if so, determine that the new working condition is the optimal working condition, set B ₁ =B _b0 , Transfer to S10-2-6;

S10-2-11: Determine whether f _2-3 (N _ge3 ) is greater than or equal to Q _t -Q ₁₀ , if yes, go to S10-2-12; if not, go to S10-2-14;

S10-2-12: set Q ₃₀ =Q _t -Q ₁₀ , then Q ₂₀ =0, enter S10-2-13;

S10-2-13: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ) , B ₀ =B ₁₀ +B ₂₀ +B ₃₀ , enter S10-2-6;

S10-2-14: Let Q ₃₀ =f _2-3 (N _ge3 ), then Q ₂₀ =Q _t -Q ₁₀ -f _2-3 (N _ge3 ), go to S10-2-13;

S10-2-15: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₁ =F _3-2 (N _ge2 , B ₃₁ =F _2-3 (N _ge3 ,Q ₃₁ ), B ₁ = B ₁₀ +B ₂₁ +B ₃₁ , enter S10-2-10;

S10-2-16: Based on the lowest value B _b0 of the total coal consumption of the whole plant under the back pressure of P _c0 , the heat load distribution Q ₂ and Q ₃ of the No. Compare benchmarks, enter S10-2-17;

S10-2-17: let P _c1 =P _c0 +1kPa, enter S10-2-18;

S10-2-18: Determine whether P _c1 is less than or equal to min(P _c,s, P _c,max ), if yes, increase the back pressure of Unit 1 by 1kPa, and enter S10-2-2; if not, Then the optimization iterative process is terminated, where P _c,s and P _c,max are the maximum back pressure value of the safe operation of the high back pressure unit and the actual maximum value that can be achieved by the back pressure adjustment means;

S10-2-19: Taking the lowest value B _b of the total coal consumption of the whole plant as the optimization target, obtain the operating back pressure P _cb of No. 1 unit, and the heat load distribution Q _2b and Q _3b of No. 2 and No. 3 units.

S10-3: Compare the lowest value of the total coal consumption of the whole plant in the iterative optimization process of S10-1 and S10-2, and take the working condition corresponding to the small value as the optimal operating condition: the operating back pressure of unit 1 P _cb , 2 and Unit 3 heat load distribution Q _2b , Q _3b .

S11: The judgment of S4 and the judgment of S6 are transferred to S11 together to judge the heat load: whether f _2-2 (N _ge2 )+f _2-3 (N _ge3 ) is greater than or equal to Q _t -Q ₁ , if so, Then it is determined that the No. 2 unit has the electric load and thermal load conditions of the middle exhaust steam extraction mode and the low-pressure cylinder zero output mode, and the No. 3 unit has the electrical load and thermal load conditions of the mid-exhaust steam extraction mode, and enters S12; As shown in the dotted line in Fig. 3; if not, the optimization iterative process is terminated;

S12: Enter the iterative optimization mode 2, and output the optimal result.

The implementation process of optimization mode 2 is shown in Figure 6.

S12-1: Formulate the initial value of the back pressure P _c0 of the high back pressure heating unit operation, and use the preset threshold of the temperature rise of the high back pressure heating network condenser corresponding to the back pressure P _c0 as the back pressure benchmark for iterative optimization;

S12-2: Calculate Q ₁₀ =f ₁ (N _ge ,P _c0 ,t ₀ ,m _cw );

S12-3: Let Q ₂₀ =f _3-2 (N _ge2 ), then Q ₃₀ =Q _t -Q ₁₀ -f _3-2 (N _ge2 );

S12-4: judge whether Q ₃₀ is greater than or equal to 0, if yes, then enter S12-5; if not, then transfer to S12-10;

S12-5: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _3-2 (N _ge2 ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ), B ₀ =B ₁₀ +B ₂₀ +B ₃₀ ;

S12-6: Let Q ₂₁ =Q ₂₀ -(f _3-2 (N _ge2 )-f _2-2 (N _ge2 )), then Q ₃₁ =Q ₃₀ +(f _3-2 (N _ge2 )-f _{2 -2} (N _ge2 ));

S12-7: Determine whether Q ₂₁ is less than or equal to 0, or whether Q ₃₁ is greater than or equal to f _2-3 (N _ge3 ), if so, terminate the optimization iterative process; if not, enter S12-8;

S12-8: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₁ =F _2-2 (N _ge2 ,Q ₂₁ ), B ₃₁ =F _2-3 (N _ge3 ,Q ₃₁ ), B ₁ = B ₁₀ +B ₂₁ +B ₃₁ , enter S12-9;

S12-9: Determine whether B ₁ is less than or equal to B ₀ , if not, determine that the original reference working condition is still the reference working condition; if so, determine that the new working condition is the optimal working condition, set B ₁ = B _b0 , and transfer to S12-6;

S12-10: Determine whether f _2-2 (N _ge2 ) is greater than or equal to Q _t -Q ₁₀ , if yes, go to S12-11; if not, go to S12-13;

S12-11: set Q ₂₀ =Q _t -Q ₁₀ , then Q ₃₀ =0, enter S12-12;

S12-12: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ), B ₀ = B ₁₀ +B ₂₀ +B ₃₀ , enter S12-6;

S12-13: Let Q ₂₀ =f _2-2 (N _ge2 ), then Q ₃₀ =Q _t -Q ₁₀ -f _2-2 (N _ge2 ), go to S12-12;

S12-14: Based on the lowest value B _b0 of the total coal consumption of the whole plant under the back pressure of P _c0 , the heat load distribution Q ₂ and Q ₃ of the No. 2 unit and No. 3 unit are obtained, which will be used as the comparison benchmark for the next back pressure iterative optimization , enter S12-15;

S12-15: Let P _c1 =P _c0 +1kPa, go to S12-16;

S12-16: Determine whether P _c1 is less than or equal to min(P _c,s, P _c,max ), if yes, increase the back pressure of Unit 1 by 1kPa, and enter S12-2; if not, terminate the optimization Iterative process, where P _c,s and P _c,max are the highest back pressure value for safe operation of high back pressure units and the actual highest value that can be achieved by back pressure adjustment means;

S12-17: Taking the lowest value B _b of the total coal consumption of the whole plant as the optimization target, obtain the operating back pressure P _cb of No. 1 unit, and the heat load distribution Q _2b and Q _3b of No. 2 and No. 3 units.

S13: The judgment of S5 and the judgment of S7 are transferred to S13 together to judge the heat load: whether f _2-2 (N _ge2 )+f _2-3 (N _ge3 ) is greater than or equal to Q _t -Q ₁ , if so, Then it is determined that Unit 2 has the electrical load and thermal load conditions for the mid-exhaust steam extraction mode in operation, and Unit 3 has the electrical load and thermal load conditions for the mid-exhaust steam extraction mode or low-pressure cylinder zero output mode, and enters S14; See the dotted line in Figure 3; if not, the optimization iterative process is terminated

S14: Enter the iterative optimization mode 3, and output the optimal result.

The implementation process of optimization mode 3 is shown in Figure 7.

S14-1 Formulate the initial value of the back pressure P _c0 of the high back pressure heating unit operation, and use the preset threshold of the temperature rise of the high back pressure heating network condenser corresponding to the back pressure P _c0 as the back pressure benchmark for iterative optimization;

S14-2 Calculate Q ₁₀ =f ₁ (N _ge ,P _c0 ,t ₀ ,m _cw );

S14-3 Let Q ₃₀ =f _3-3 (N _ge3 ), then Q ₂₀ =Q _t -Q ₁₀ -f _3-3 (N _ge3 );

_S14-4 judges whether Q20 is greater than or equal to 0, if so, then enters S14-5; if not, then proceeds to S14-10;

S14-5: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _3-3 (N _ge3 ), B ₀ =B ₁₀ +B ₂₀ +B ₃₀ ;

S14-6: Let Q ₃₁ =Q ₃₀ -(f _3-3 (N _ge3 )-f _2-3 (N _ge3 )), then Q ₂₁ =Q ₂₀ +(f _3-3 (N _ge3 )-f _{2 -3} (N _ge3 ));

S14-7: Determine whether Q ₃₁ is less than or equal to 0, or whether Q ₂₁ is greater than or equal to f _2-2 (N _ge2 ), if so, terminate the optimization iterative process; if not, enter S14-8;

S14-8: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₁ =F _2-2 (N _ge2 ,Q ₂₁ ), B ₃₁ =F _2-3 (N _ge3 ,Q ₃₁ ), B ₁ = B ₁₀ +B ₂₁ +B ₃₁ , enter S14-9;

S14-9: Determine whether B ₁ is less than or equal to B ₀ , if not, determine that the original reference working condition is still the reference working condition; if so, determine that the new working condition is the optimal working condition, set B ₁ = B _b0 , and transfer to S14-6;

S14-10: Determine whether f _2-3 (N _ge3 ) is greater than or equal to Q _t -Q ₁₀ , if yes, go to S14-11; if not, go to S14-13;

S14-11: set Q ₃₀ =Q _t -Q ₁₀ , then Q ₂₀ =0, enter S14-12;

S14-12: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ), B ₀ = B ₁₀ +B ₂₀ +B ₃₀ , enter S14-6;

S14-13: Let Q ₃₀ =f _2-3 (N _ge3 ), then Q ₂₀ =Q _t -Q ₁₀ -f _3-3 (N _ge3 ), enter S14-12;

S14-14: Based on the lowest value B _b0 of the total coal consumption of the whole plant under the back pressure of P _c0 , the heat load distribution Q ₂ and Q ₃ of the second unit and the third unit are obtained, which will be used as the comparison benchmark for the next back pressure iterative optimization, Enter S14-15;

S14-15: let P _c1 =P _c0 +1kPa, go to S14-16;

S14-16: Determine whether P _c1 is less than or equal to min(P _c,s, P _c,max ), if yes, increase the back pressure of Unit 1 by 1kPa, and enter S14-2; if not, terminate the optimization Iterative process, where P _c,s and P _c,max are the highest back pressure value for safe operation of high back pressure units and the actual highest value that can be achieved by back pressure adjustment means;

S14-17: Taking the lowest value B _b of the total coal consumption of the whole plant as the optimization goal, obtain the operating back pressure P _cb of the No. 1 unit, and the heat load distribution Q _2b and Q _3b of the No. 2 unit and No. 3 unit.

S15: Carry out thermal load judgment: whether f _2-2 (N _ge2 )+f _2-3 (N _ge3 ) is greater than or equal to Q _t -Q ₁ , if yes, it is determined that Unit 2 and Unit 3 are equipped with mid-exhaust pumping equipment for operation Electric load and thermal load conditions of the steam mode, enter S16; if not, the optimization iteration process is terminated

S16: Enter iterative optimization mode 4, and output the optimal result

The implementation process of optimization mode 4 is shown in Figure 8.

S16-1 Formulate the initial value of the back pressure P _c0 of the high back pressure heating unit operation, and use the preset threshold of the temperature rise of the high back pressure heating network condenser corresponding to the back pressure P _c0 as the back pressure benchmark for iterative optimization;

S16-2 Calculate Q ₁₀ =f ₁ (N _ge ,P _c0 ,t ₀ ,m _cw );

S16-3 Let Q ₂₀ =f _2-2 (N _ge2 ), then Q ₃₀ =Q _t -Q ₁₀ -f _2-2 (N _ge2 );

S16-4 judges whether Q ₃₀ is greater than or equal to 0, if so, then enters S16-5; if not, then proceeds to S16-10;

S16-5: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ), B ₀ = B ₁₀ +B ₂₀ +B ₃₀ ;

S16-6: Let Q ₂₁ =Q ₂₀ -0.1×Q ₂₀ , then Q ₃₁ =Q ₃₀ +0.1×Q ₂₀ ;

S16-7: Determine whether Q ₂₁ is equal to 0, or whether Q ₃₁ is equal to f _2-3 (N _ge3 ), if so, terminate the optimization iterative process; if not, enter S16-8;

S16-8: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₁ =F _2-2 (N _ge2 ,Q ₂₁ ), B ₃₁ =F _2-3 (N _ge3 ,Q ₃₁ ), B ₁ = B ₁₀ +B ₂₁ +B ₃₁ , enter S16-9;

S16-9: Determine whether B ₁ is less than or equal to B ₀ , if not, determine that the original reference working condition is still the reference working condition; if so, determine that the new working condition is the optimal working condition, set B ₁ = B _b0 , and transfer to S16-6;

S16-10: let Q ₂₀ =Q _t -Q ₁₀ , Q ₃₀ =0, enter into S16-11;

S16-11: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ), B ₀ = B ₁₀ +B ₂₀ +B ₃₀ , enter S16-6;

S12-12: Based on the lowest value B _b0 of the total coal consumption of the whole plant under the back pressure of P _c0 , the thermal load distribution Q ₂ and Q ₃ of the No. 2 unit and No. 3 unit are obtained, which will be used as the comparison benchmark for the next back pressure iterative optimization , enter S16-13;

S16-13: let P _c1 =P _c0 +1kPa, go to S16-14;

S16-14: Determine whether P _c1 is less than or equal to min(P _c,s, P _c,max ), if yes, increase the back pressure of Unit 1 by 1kPa, and enter S16-2; otherwise, terminate the optimization iteration process, where P _c,s and P _c,max are the maximum back pressure value for the safe operation of the high back pressure unit and the actual maximum value that can be achieved by the back pressure adjustment means;

S16-15: Taking the lowest value B _b of the total coal consumption of the whole plant as the optimization goal, obtain the operating back pressure P _cb of No. 1 unit, and the heat load distribution Q _2b and Q _3b of No. 2 and No. 3 units.

The output results of S10, S12, S14, and S16 are the optimal operation mode: the operating back pressure of the high back pressure heating unit, the heat load borne by other units and the corresponding heating mode.

In order to realize the above-mentioned embodiments, the embodiments of the present disclosure propose a thermal power plant heating load distribution optimization device with multiple units and multiple heating modes. The device includes a first determination module, a second determination module and an optimization module. The first determination module is configured to determine an optimization objective function of the multi-unit, multi-heat supply mode thermal power plant heating load distribution optimization method, and the optimization objective function includes the total standard coal consumption of the whole plant. The second determination module is configured to determine the correlation characteristics of electric load-heat supply load and electric load-heat supply load-standard coal consumption of each cogeneration unit under different heat supply modes. The optimization module is configured to perform plant-level operation optimization under given conditions of heat supply load and different sub-generator loads according to the total standard coal consumption of the whole plant.

Regarding the apparatus in the foregoing embodiments, the specific manner in which each module executes operations has been described in detail in the embodiments related to the method, and will not be described in detail here.

In order to implement the above embodiments, the embodiments of the present disclosure further provide an electronic device. An electronic device includes a processor and a memory for storing instructions executable by the processor. Wherein, the processor is configured to execute the instructions, so as to realize the method for optimizing heat supply load distribution of thermal power plants with multiple units and multiple heating modes in the first aspect of the present disclosure.

As an example, FIG. 9 is a block diagram of an electronic device 900 according to an embodiment of the present disclosure. As shown in FIG. 9, the above-mentioned electronic device 900 may further include:

The memory 910 and the processor 920, the bus 930 connecting different components (including the memory 910 and the processor 920), the memory 910 stores a computer program, and when the processor 920 executes the program, the multi-unit, Optimization method for heat load distribution of thermal power plant with multiple heat supply modes.

Bus 930 represents one or more of several types of bus structures, including a memory bus or memory controller, a peripheral bus, an accelerated graphics port, a processor, or a local bus using any of a variety of bus structures. These architectures include, by way of example, but are not limited to Industry Standard Architecture (ISA) bus, Micro Channel Architecture (MAC) bus, Enhanced ISA bus, Video Electronics Standards Association (VESA) local bus, and Peripheral Component Interconnect ( PCI) bus.

Electronic device 900 typically includes a variety of electronic device readable media. These media can be any available media that can be accessed by electronic device 900 and include both volatile and nonvolatile media, removable and non-removable media.

Memory 910 may also include computer system readable media in the form of volatile memory, such as random access memory (RAM) 940 and/or cache memory 950 . The server 900 may further include other removable/non-removable, volatile/nonvolatile computer system storage media. By way of example only, storage system 960 may be used to read from and write to non-removable, non-volatile magnetic media (not shown in FIG. 9, commonly referred to as a "hard drive"). Although not shown in FIG. 9, a disk drive for reading and writing to removable nonvolatile disks (e.g., "floppy disks") may be provided, as well as for removable nonvolatile optical disks (e.g., CD-ROM, DVD-ROM). or other optical media) CD-ROM drive. In these cases, each drive may be connected to bus 930 through one or more data media interfaces. The memory 910 may include at least one program product having a set (eg, at least one) of program modules configured to perform the functions of various embodiments of the present disclosure.

Program/utility 980 having a set (at least one) of program modules 970, such as may be stored in memory 910, such program modules 970 including - but not limited to - an operating system, one or more application programs, other program Modules and program data, each or some combination of these examples may include the implementation of the network environment. The program modules 970 generally perform the functions and/or methods of the embodiments described in the present disclosure.

The electronic device 900 can also communicate with one or more external devices 990 (such as keyboards, pointing devices, displays 991, etc.), and can also communicate with one or more devices that enable the user to interact with the electronic device 900, and/or communicate with Any device (eg, network card, modem, etc.) that enables the electronic device 900 to communicate with one or more other computing devices. Such communication may occur through input/output (I/O) interface 992 . Moreover, the electronic device 900 can also communicate with one or more networks (such as a local area network (LAN), a wide area network (WAN) and/or a public network such as the Internet) through the network adapter 993 . As shown, the network adapter 993 communicates with other modules of the electronic device 900 through the bus 930 . It should be appreciated that although not shown in FIG. 9 , other hardware and/or software modules may be used in conjunction with electronic device 900, including but not limited to: microcode, device drivers, redundant processing units, external disk drive arrays, RAID systems, tape Drives and data backup storage systems, etc.

The processor 920 executes various functional applications and data processing by executing programs stored in the memory 910 .

It should be noted that, for the implementation process and technical principle of the electronic equipment in this embodiment, please refer to the above-mentioned explanation of the multi-unit, multi-heating mode thermal power plant heating load distribution optimization method in the embodiment of the present disclosure, and will not be repeated here.

In order to realize the above-mentioned embodiments, the embodiments of the present disclosure also propose a computer-readable storage medium on which a computer program is stored, and the feature is that, when the program is executed by a processor, the multi-unit, multi-unit Optimization method for heating load distribution of thermal power plants in heating mode.

In order to realize the above-mentioned embodiments, the embodiment of the present disclosure also proposes a computer program product, including a computer program, when the computer program is executed by a processor, it realizes the multi-unit, multi-heating mode thermal power plant heat supply in the first aspect of the present disclosure Load distribution optimization method.

It should be understood that steps may be reordered, added or deleted using the various forms of flow shown above. For example, each step described in the present disclosure may be executed in parallel, sequentially, or in a different order, as long as the desired result of the technical solution disclosed in the present disclosure can be achieved, no limitation is imposed herein.

The specific implementation manners described above do not limit the protection scope of the present disclosure. It should be apparent to those skilled in the art that various modifications, combinations, sub-combinations and substitutions may be made depending on design requirements and other factors. Any modifications, equivalent replacements and improvements made within the spirit and principles of the present disclosure shall be included within the protection scope of the present disclosure.

Claims (14)

1. A multi-unit, multi-heating mode thermal power plant heating load distribution optimization method, including:

   Determine the optimization objective function of the heat supply load distribution optimization method of the thermal power plant with multiple units and multiple heating modes, the optimization objective function includes the total standard coal consumption of the whole plant;

   Determine the correlation characteristics of electricity load-heat supply load, electricity load-heat supply load-standard coal consumption of each cogeneration unit under different heating modes;

   According to the total standard coal consumption of the whole plant, the plant-level operation optimization is carried out under the given conditions of heating load and electrical load of different sub-units.
2. The method according to claim 1, wherein, when determining the optimization objective function of the heat supply load distribution optimization method of the thermal power plant with multiple units and multiple heating modes, the optimization objective function includes the total consumption of standard coal of the whole plant , the calculation method of the total profit value M of the whole plant is as follows:

   M=H+EC=Q _t ×h+(N _ge1 +N _ge2 +N _ge3 )×e-(B ₁ +B ₂ +B ₃ )×c

   Among them, H is heat sales revenue; Q _t is total external heat supply load, h is ex-factory heat price, E is electricity sales revenue, N _ge1 , N _ge2 , N _ge3 are real-time power generation loads of three units respectively, and e is grid connection Electricity price, C is the cost of standard coal consumption, B ₁ , B ₂ , B ₃ are the real-time standard coal consumption of the three units respectively, and c is the unit price of standard coal.
3. The method according to claim 1, wherein, in determining the correlation characteristics of electric load-heat supply load, electric load-heat supply load-standard coal consumption under different heating modes of each cogeneration unit, when the unit adopts high When the back pressure cascade heating is used, the calculation method of the heating load Q is as follows:

   

   t _s =f(P _c )=-0.029×P _c ² +2.28×P _c +26.13

   Among them, m _cw is the flow rate of circulating water in and out of the high back pressure heating network condenser, C _p is the constant pressure specific heat capacity of circulating water in the heating network, t ₀ and t ₁ are the heat flow in and out of the high back pressure heating network condenser respectively Network circulating water temperature, δt is the heat exchange end difference of the high back pressure heating network condenser, t _s is the condensate temperature of the exhaust steam after the heat release of the high back pressure heating network condenser, P _c is the high back pressure machine the running backpressure of the group;

   The heat supply load Q is a multivariate function of the electric load N _ge1 , the temperature t ₀ and t ₁ of the heating network circulating water entering and leaving the high back pressure heating network condenser, the flow rate m _cw of the heating network circulating water and the operating back pressure P _c , represented by the following formula:

   Q＝f ₁ (N _ge1 ,P _c ,m _cw ,t ₀ )

   The standard coal consumption B under the condition of dual supply of electricity and heat of the high back pressure heating unit is a binary function of the electric load N _ge1 and the operating back pressure P _c , which is characterized by the following formula:

   B=F ₁ (N _ge1 ,P _c ).
4. The method according to claim 1, wherein, in determining the correlation characteristics of electric load-heat supply load, electric load-heat supply load-standard coal consumption under different heating modes of each cogeneration unit, when the unit adopts In the exhaust extraction mode, the electric load N _ge is determined, and the heating load is between Q between 0 and the maximum value Q _max :

   0≤Q≤Q _max ＝f ₂ (N _ge )

   The standard coal consumption B is a binary function of the electrical load N _ge and the heating load Q, which is characterized by the following formula:

   B=F ₂ (N _ge ,Q).
5. The method according to claim 1, wherein, in determining the correlation characteristics of electric load-heat supply load, electric load-heat supply load-standard coal consumption under different heating modes of each cogeneration unit, when the unit adopts low-pressure In the heating mode with zero cylinder output, the unit operates with constant heat and electricity, and the heating load Q is a one-element linear function of the electric load N _ge , which is represented by the following formula:

   Q＝f ₃ (N _ge )

   The standard coal consumption B is a binary function of the electrical load N _ge and the heating load Q, which is characterized by the following formula:

   B=F ₃ (N _ge )

   Under the condition of boiler design output D _ms0 , the unit adopts low-pressure cylinder zero output heating mode, and the operating range of electric load is:

   0≦N _ge ≦N _ge,max =f ₄ (D _ms0 ).
6. The method according to any one of claims 1 to 5, wherein, according to the total standard coal consumption of the whole plant, the plant-level operation optimization under the given conditions of heating load and different sub-generator loads includes:

   S301, input boundary parameters, the boundary parameters include heating return water temperature t ₀ , heating network circulating water flow m _cw , total heating load Q _t , generating load N _ge1 of Unit 1, generating load N _ge2 of Unit 2, and generating load of Unit 3 Unit generating load N _ge3 ;

   S302, formulating an initial value of back pressure P _c0 for the operation of the high back pressure heating unit, and using the preset threshold value of the temperature rise of the high back pressure heating network condenser corresponding to the back pressure P _c0 as the back pressure reference for iterative optimization;

   S303, calculate Q ₁ =f ₁ (N _ge ,P _c0 ,t ₀ ,m _cw ), where Q ₁ is the heat supply load of unit 1, N _ge is the real-time power generation load, P _c0 is the running back pressure, t ₀ is the circulating water temperature of the heating network entering the high back pressure heating network condenser, m _cw is the flow rate of the circulating water of the heating network entering and leaving the high back pressure heating network condenser;

   S304, for the No. 2 unit, judge whether the electric load conditions of the middle exhaust steam extraction mode and the zero output mode of the low-pressure cylinder are met, if yes, go to S305, if not, go to S315;

   S305, for unit No. 3, judge whether the electric load conditions of the two modes of mid-exhaust steam extraction or low-pressure cylinder zero output are met, if yes, go to S306, if not, go to S311 and S315;

   S306, determine the heat load of the No. 2 unit, and judge whether the No. 2 unit has the heat load conditions of the two modes of mid-exhaust steam extraction and low-pressure cylinder zero output, if yes, enter S307, and if not, execute S315;

   S307, determine the heat load of the No. 3 unit, and judge whether the No. 3 unit has the heat load conditions of the two modes of mid-exhaust steam extraction and low-pressure cylinder zero output, if yes, enter S308, and if not, execute S315;

   S308, determine the heat load, and judge whether f _2-2 (N _ge2 )+f _3-3 (N _ge3 ) is greater than or equal to Q _t -Q ₁ , if yes, enter S309, if not, terminate the optimization iterative process;

   S309, determine the heat load, and judge whether f _2-3 (N _ge3 )+f _3-2 (N _ge2 ) is greater than or equal to Q _t -Q ₁ , if yes, enter S310, if not, terminate the optimization iterative process;

   S310, No. 2 unit and No. 3 unit have two heat supply modes: mid-exhaust steam extraction or low-pressure cylinder zero output, enter the first iterative optimization mode, and output the optimal result;

   S311, judging whether the No. 2 unit has the electric load and heat load conditions of the mid-exhaust steam extraction mode or the low-pressure cylinder zero output mode for operation, and at the same time, whether the No. 3 unit has the electric load and thermal load conditions of the mid-exhaust steam extraction mode for operation condition, then enter S312, if not, then terminate the optimization iterative process;

   S312, enter the second iterative optimization mode, and output the optimal result;

   S313, judging whether the No. 2 unit has the electric load and thermal load conditions of the mid-exhaust steam extraction mode in operation, and at the same time, whether the No. 3 unit has the electric load and thermal load of the mid-exhaust steam extraction mode and low-pressure cylinder zero output mode. condition, if so, then enter S314, if not, then terminate the optimization iterative process;

   S314, performing a third iterative optimization mode, and outputting an optimal result;

   S315, judging whether the No. 2 unit and the No. 3 unit have the electric load and heat load conditions of the exhaust steam extraction mode in operation, and if so, proceed to S316, and if not, terminate the optimization iterative process;

   S316. Perform a fourth iterative optimization mode, and output an optimal result.
7. The method according to claim 6, wherein the first iterative optimization mode comprises:

   S10-1-1: Establish the initial value of the back pressure P _c0 of the high back pressure heating unit operation, and take the back pressure P _c0 corresponding to the temperature rise of the high back pressure heating network condenser at 5°C as the back pressure benchmark for iterative optimization.

   S10-1-2: Calculate Q ₁₀ =f ₁ (N _ge ,P _c0 ,t ₀ ,m _cw );

   S10-1-3: let Q ₂₀ =f _3-2 (N _ge2 ), then Q ₃₀ =Q _t -Q ₁₀ -f _3-2 (N _ge2 );

   S10-1-4: judge whether Q ₃₀ is greater than 0, if so, then enter S10-1-5; if not, then transfer to S10-1-11;

   S10-1-5: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _3-2 (N _ge2 ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ), B ₀ =B ₁₀ +B ₂₀ +B ₃₀ , enter S10-1-6;

   S10-1-6: Let Q ₂₁ =Q ₂₀ -(f _3-2 (N _ge2 )-f _2-2 (N _ge2 )), then Q ₃₁ =Q ₃₀ +(f _3-2 (N _ge2 )- f _2-2 (N _ge2 )), enter S10-1-7;

   S10-1-7: Judging whether Q ₂₁ is less than 0, or whether Q ₃₁ is greater than or equal to f _3-3 (N _ge3 ), if yes, then the optimization iteration process is terminated; if not, then transfer to S10-1-8;

   S10-1-8: Judging whether Q ₃₁ is less than or equal to f _2-3 (N _ge3 ), if yes, proceed to S10-1-9; if not, proceed to S10-1-15;

   S10-1-9: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₁ =F _2-2 (N _ge2 ,Q ₂₁ ), B ₃₁ =F _2-3 (N _ge3 ,Q ₃₁ ) , B ₁ =B ₁₀ +B ₂₁ +B ₃₁ , enter S10-1-10;

   S10-1-10: Determine whether B ₁ is less than or equal to B ₀ , if not, determine that the original reference working condition is still the reference working condition; if so, determine that the new working condition is the optimal working condition, set B ₁ =B _b0 , Transfer to S10-1-6;

   S10-1-11: Determine whether f _2-2 (N _ge2 ) is greater than or equal to Q _t -Q ₁₀ , if yes, go to S10-1-12; if not, go to S10-1-14;

   S10-1-12: set Q ₂₀ =Q _t -Q ₁₀ , then Q ₃₀ =0, enter S10-1-13;

   S10-1-13: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ) , B ₀ =B ₁₀ +B ₂₀ +B ₃₀ , enter S10-1-6;

   S10-1-14: Let Q ₂₀ =f _2-2 (N _ge2 ), then Q ₃₀ =Q _t -Q ₁₀ -f _2-2 (N _ge2 ), go to S10-1-13;

   S10-1-15: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₁ =F _2-2 (N _ge2 ,Q ₂₁ ), B ₃₁ =F _3-3 (N _ge3 ), B ₁ =B ₁₀ +B ₂₁ +B ₃₁ , enter S10-1-10;

   S10-1-16: Based on the lowest value B _b0 of the total coal consumption of the whole plant under the back pressure of P _c0 , the thermal load distribution Q ₂ and Q ₃ of the No. 2 unit and No. 3 unit are obtained as the next back pressure iterative optimization The benchmark of comparison, enter S10-1-17;

   S10-1-17: let P _c1 =P _c0 +1kPa, enter S10-1-18;

   S10-1-18: Determine whether P _c1 is less than or equal to min(P _c,s, P _c,max ), if yes, increase the back pressure of Unit 1 by 1kPa, and enter S10-1-2; if not, Then the optimization iterative process is terminated, where P _c,s and P _c,max are the maximum back pressure value of the safe operation of the high back pressure unit and the actual maximum value that can be achieved by the back pressure adjustment means;

   S10-1-19: Taking the lowest value B _b of the total coal consumption of the whole plant as the optimization target, obtain the operating back pressure P _cb of No. 1 unit, the heat load distribution Q _2b and Q _3b of No. 2 and No. 3 units;

   S10-2-1: Formulate the initial value of the back pressure P _c0 of the high back pressure heating unit operation, and use the preset threshold of the temperature rise of the high back pressure heating network condenser corresponding to the back pressure P _c0 as the back pressure benchmark for iterative optimization;

   S10-2-2: Calculate Q ₁₀ =f ₁ (N _ge ,P _c0 ,t ₀ ,m _cw );

   S10-2-3: let Q ₃₀ =f _3-3 (N _ge3 ), then Q ₂₀ =Q _t -Q ₁₀ -f _3-3 (N _ge3 );

   S10-2-4: Determine whether Q ₂₀ is greater than or equal to 0, if yes, then enter S10-2-5; if not, then proceed to S10-2-11;

   S10-2-5: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _3-3 (N _ge3 ), B ₀ =B ₁₀ +B ₂₀ +B ₃₀ , enter S10-2-6;

   S10-2-6: Let Q ₃₁ =Q ₃₀ -(f _3-3 (N _ge3 )-f _2-3 (N _ge3 )), then Q ₂₁ =Q ₂₀ +(f _3-3 (N _ge3 )- f _2-3 (N _ge3 )), enter S10-2-7;

   S10-2-7: Determine whether Q ₃₁ is less than or equal to 0, or whether Q ₂₁ is greater than or equal to f _3-2 (N _ge2 ), if yes, terminate the optimization iterative process; if not, transfer to S10-2-8;

   S10-2-8: Judging whether Q ₂₁ is less than or equal to f _2-2 (N _ge2 ), if yes, proceed to S10-2-9; if not, proceed to S10-2-15;

   S10-2-9: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₁ =F _2-2 (N _ge2 ,Q ₂₁ ), B ₃₁ =F _2-3 (N _ge3 ,Q ₃₁ ) , B ₁ =B ₁₀ +B ₂₁ +B ₃₁ , enter S10-2-10;

   S10-2-10: Determine whether B ₁ is less than or equal to B ₀ , if not, determine that the original reference working condition is still the reference working condition; if so, determine that the new working condition is the optimal working condition, set B ₁ =B _b0 , Transfer to S10-2-6;

   S10-2-11: Determine whether f _2-3 (N _ge3 ) is greater than or equal to Q _t -Q ₁₀ , if yes, go to S10-2-12; if not, go to S10-2-14;

   S10-2-12: set Q ₃₀ =Q _t -Q ₁₀ , then Q ₂₀ =0, enter S10-2-13;

   S10-2-13: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ) , B ₀ =B ₁₀ +B ₂₀ +B ₃₀ , enter S10-2-6;

   S10-2-14: Let Q ₃₀ =f _2-3 (N _ge3 ), then Q ₂₀ =Q _t -Q ₁₀ -f _2-3 (N _ge3 ), go to S10-2-13;

   S10-2-15: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₁ =F _3-2 (N _ge2 , B ₃₁ =F _2-3 (N _ge3 ,Q ₃₁ ), B ₁ = B ₁₀ +B ₂₁ +B ₃₁ , enter S10-2-10;

   S10-2-16: Based on the lowest value B _b0 of the total coal consumption of the whole plant under the back pressure of P _c0 , the heat load distribution Q ₂ and Q ₃ of the No. Compare benchmarks, enter S10-2-17;

   S10-2-17: let P _c1 =P _c0 +1kPa, enter S10-2-18;

   S10-2-18: Determine whether P _c1 is less than or equal to min(P _c,s, P _c,max ), if yes, increase the back pressure of Unit 1 by 1kPa, and enter S10-2-2; if not, Then the optimization iterative process is terminated, where P _c,s and P _c,max are the maximum back pressure value of the safe operation of the high back pressure unit and the actual maximum value that can be achieved by the back pressure adjustment means;

   S10-2-19: Taking the lowest value B _b of the total coal consumption of the whole plant as the optimization target, obtain the operating back pressure P _cb of No. 1 unit, the heat load distribution Q _2b and Q _3b of No. 2 and No. 3 units;
8. The method according to claim 6 or 7, wherein the second iterative optimization mode comprises:

   S12-1: Formulate the initial value of the back pressure P _c0 of the high back pressure heating unit operation, and use the preset threshold of the temperature rise of the high back pressure heating network condenser corresponding to the back pressure P _c0 as the back pressure benchmark for iterative optimization;

   S12-2: Calculate Q ₁₀ =f ₁ (N _ge ,P _c0 ,t ₀ ,m _cw );

   S12-3: Let Q ₂₀ =f _3-2 (N _ge2 ), then Q ₃₀ =Q _t -Q ₁₀ -f _3-2 (N _ge2 );

   S12-4: judge whether Q ₃₀ is greater than or equal to 0, if yes, then enter S12-5; if not, then transfer to S12-10;

   S12-5: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _3-2 (N _ge2 ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ), B ₀ =B ₁₀ +B ₂₀ +B ₃₀ ;

   S12-6: Let Q ₂₁ =Q ₂₀ -(f _3-2 (N _ge2 )-f _2-2 (N _ge2 )), then Q ₃₁ =Q ₃₀ +(f _3-2 (N _ge2 )-f _{2 -2} (N _ge2 ));

   S12-7: Determine whether Q ₂₁ is less than or equal to 0, or whether Q ₃₁ is greater than or equal to f _2-3 (N _ge3 ), if so, terminate the optimization iterative process; if not, enter S12-8;

   S12-8: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₁ =F _2-2 (N _ge2 ,Q ₂₁ ), B ₃₁ =F _2-3 (N _ge3 ,Q ₃₁ ), B ₁ = B ₁₀ +B ₂₁ +B ₃₁ , enter S12-9;

   S12-9: Determine whether B ₁ is less than or equal to B ₀ , if not, determine that the original reference working condition is still the reference working condition; if so, determine that the new working condition is the optimal working condition, set B ₁ = B _b0 , and transfer to S12-6;

   S12-10: Determine whether f _2-2 (N _ge2 ) is greater than or equal to Q _t -Q ₁₀ , if yes, go to S12-11; if not, go to S12-13;

   S12-11: set Q ₂₀ =Q _t -Q ₁₀ , then Q ₃₀ =0, enter S12-12;

   S12-12: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ), B ₀ = B ₁₀ +B ₂₀ +B ₃₀ , enter S12-6;

   S12-13: Let Q ₂₀ =f _2-2 (N _ge2 ), then Q ₃₀ =Q _t -Q ₁₀ -f _2-2 (N _ge2 ), go to S12-12;

   S12-14: Based on the lowest value B _b0 of the total coal consumption of the whole plant under the back pressure of P _c0 , the heat load distribution Q ₂ and Q ₃ of the No. 2 unit and No. 3 unit are obtained, which will be used as the comparison benchmark for the next back pressure iterative optimization , enter S12-15;

   S12-15: Let P _c1 =P _c0 +1kPa, go to S12-16;

   S12-16: Determine whether P _c1 is less than or equal to min(P _c,s , P _c,max ), if yes, increase the back pressure of Unit 1 by 1kPa, and enter S12-2; if not, terminate the optimization Iterative process, where P _c,s and P _c,max are the highest back pressure value for safe operation of high back pressure units and the actual highest value that can be achieved by back pressure adjustment means;

   S12-17: Taking the lowest value B _b of the total coal consumption of the whole plant as the optimization target, obtain the operating back pressure P _cb of No. 1 unit, and the heat load distribution Q _2b and Q _3b of No. 2 and No. 3 units.
9. The method according to any one of claims 6 to 8, wherein the third iterative optimization mode comprises:

   S14-1: Formulate the initial value of the back pressure P _c0 of the high back pressure heating unit operation, and use the preset threshold of the temperature rise of the high back pressure heating network condenser corresponding to the back pressure P _c0 as the back pressure benchmark for iterative optimization;

   S14-2: Calculate Q ₁₀ =f ₁ (N _ge ,P _c0 ,t ₀ ,m _cw );

   S14-3: Let Q ₃₀ =f _3-3 (N _ge3 ), then Q ₂₀ =Q _t -Q ₁₀ -f _3-3 (N _ge3 );

   S14-4: Determine whether Q ₂₀ is greater than or equal to 0, if yes, then enter S14-5; if not, then proceed to S14-10;

   S14-5: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _3-3 (N _ge3 ), B ₀ =B ₁₀ +B ₂₀ +B ₃₀ ;

   S14-6: Let Q ₃₁ =Q ₃₀ -(f _3-3 (N _ge3 )-f _2-3 (N _ge3 )), then Q ₂₁ =Q ₂₀ +(f _3-3 (N _ge3 )-f _{2 -3} (N _ge3 ));

   S14-7: Determine whether Q ₃₁ is less than or equal to 0, or whether Q ₂₁ is greater than or equal to f _2-2 (N _ge2 ), if so, terminate the optimization iterative process; if not, enter S14-8;

   S14-8: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₁ =F _2-2 (N _ge2 ,Q ₂₁ ), B ₃₁ =F _2-3 (N _ge3 ,Q ₃₁ ), B ₁ = B ₁₀ +B ₂₁ +B ₃₁ , enter S14-9;

   S14-9: Determine whether B ₁ is less than or equal to B ₀ , if not, determine that the original reference working condition is still the reference working condition; if so, determine that the new working condition is the optimal working condition, set B ₁ = B _b0 , and transfer to S14-6;

   S14-10: Determine whether f _2-3 (N _ge3 ) is greater than or equal to Q _t -Q ₁₀ , if yes, go to S14-11; if not, go to S14-13;

   S14-11: set Q ₃₀ =Q _t -Q ₁₀ , then Q ₂₀ =0, enter S14-12;

   S14-12: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ), B ₀ = B ₁₀ +B ₂₀ +B ₃₀ , enter S14-6;

   S14-13: Let Q ₃₀ =f _2-3 (N _ge3 ), then Q ₂₀ =Q _t -Q ₁₀ -f _3-3 (N _ge3 ), enter S14-12;

   S14-14: Based on the lowest value B _b0 of the total coal consumption of the whole plant under the back pressure of P _c0 , the heat load distribution Q ₂ and Q ₃ of the second unit and the third unit are obtained, which will be used as the comparison benchmark for the next back pressure iterative optimization, Enter S14-15;

   S14-15: let P _c1 =P _c0 +1kPa, enter S14-16;

   S14-16: Determine whether P _c1 is less than or equal to min(P _c,s , P _c,max ), if yes, increase the back pressure of Unit 1 by 1kPa, and enter S14-2; if not, terminate the optimization Iterative process, where P _c,s and P _c,max are the highest back pressure value for safe operation of high back pressure units and the actual highest value that can be achieved by back pressure adjustment means;

   S14-17: Taking the lowest value B _b of the total coal consumption of the whole plant as the optimization goal, obtain the operating back pressure P _cb of the No. 1 unit, and the heat load distribution Q _2b and Q _3b of the No. 2 unit and No. 3 unit.
10. The method according to any one of claims 6 to 9, wherein the fourth iterative optimization mode comprises:

    S16-1: Formulate the initial value of the back pressure P _c0 of the high back pressure heating unit operation, and use the preset threshold of the temperature rise of the high back pressure heating network condenser corresponding to the back pressure P _c0 as the back pressure benchmark for iterative optimization;

    S16-2: Calculate Q ₁₀ =f ₁ (N _ge ,P _c0 ,t ₀ ,m _cw );

    S16-3: let Q ₂₀ =f _2-2 (N _ge2 ), then Q ₃₀ =Q _t -Q ₁₀ -f _2-2 (N _ge2 );

    S16-4: judge whether Q ₃₀ is greater than or equal to 0, if yes, then enter S16-5; if not, then transfer to S16-10;

    S16-5: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ), B ₀ = B ₁₀ +B ₂₀ +B ₃₀ ;

    S16-6: Let Q ₂₁ =Q ₂₀ -0.1×Q ₂₀ , then Q ₃₁ =Q ₃₀ +0.1×Q ₂₀ ;

    S16-7: Determine whether Q ₂₁ is equal to 0, or whether Q ₃₁ is equal to f _2-3 (N _ge3 ), if so, terminate the optimization iterative process; if not, enter S16-8;

    S16-8: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₁ =F _2-2 (N _ge2 ,Q ₂₁ ), B ₃₁ =F _2-3 (N _ge3 ,Q ₃₁ ), B ₁ = B ₁₀ +B ₂₁ +B ₃₁ , enter S16-9;

    S16-9: Determine whether B ₁ is less than or equal to B ₀ , if not, determine that the original reference working condition is still the reference working condition; if so, determine that the new working condition is the optimal working condition, set B ₁ = B _b0 , and transfer to S16-6;

    S16-10: let Q ₂₀ =Q _t -Q ₁₀ , Q ₃₀ =0, enter into S16-11;

    S16-11: Calculation: B ₁₀ =F ₁ (N _ge ,P _c0 ), B ₂₀ =F _2-2 (N _ge2 ,Q ₂₀ ), B ₃₀ =F _2-3 (N _ge3 ,Q ₃₀ ), B ₀ = B ₁₀ +B ₂₀ +B ₃₀ , enter S16-6;

    S12-12: Based on the lowest value B _b0 of the total coal consumption of the whole plant under the back pressure of P _c0 , the thermal load distribution Q ₂ and Q ₃ of the No. 2 unit and No. 3 unit are obtained, which will be used as the comparison benchmark for the next back pressure iterative optimization , enter S16-13;

    S16-13: let P _c1 =P _c0 +1kPa, go to S16-14;

    S16-14: Determine whether P _c1 is less than or equal to min(P _c,s , P _c,max ), if yes, increase the back pressure of Unit 1 by 1kPa, and enter S16-2; otherwise, terminate the optimization iteration process, where P _c,s and P _c,max are the maximum back pressure value for the safe operation of the high back pressure unit and the actual maximum value that can be achieved by the back pressure adjustment means;

    S16-15: Taking the lowest value B _b of the total coal consumption of the whole plant as the optimization goal, obtain the operating back pressure P _cb of No. 1 unit, and the heat load distribution Q _2b and Q _3b of No. 2 and No. 3 units.
11. A thermal power plant heating load distribution optimization device with multiple units and multiple heating modes, including:

    The first determination module is configured to determine the optimization objective function of the thermal power plant heating load distribution optimization method with multiple units and multiple heating modes, and the optimization objective function includes the total standard coal consumption of the whole plant;

    The second determination module is configured to determine the correlation characteristics of electric load-heat supply load, electric load-heat supply load-standard coal consumption of each cogeneration unit under different heating modes; and

    The optimization module is configured to perform plant-level operation optimization under the given conditions of heating load and different sub-generator loads according to the total standard coal consumption of the whole plant.
12. An electronic device, characterized in that it comprises:

    processor; and

    memory for storing said processor-executable instructions;

    Wherein, the processor is configured to execute the instructions, so as to realize the multi-unit, multi-heating mode thermal power plant heating load distribution optimization method according to any one of claims 1 to 10.
13. A computer-readable storage medium, on which a computer program is stored, characterized in that, when the program is executed by a processor, the heat supply load of a thermal power plant with multiple units and multiple heating modes according to any one of claims 1 to 10 is realized Allocation optimization method.
14. A computer program product, characterized in that it includes a computer program, and when the computer program is executed by a processor, it realizes the heat supply load distribution of a thermal power plant with multiple units and multiple heating modes according to any one of claims 1 to 10 Optimization.

Images