CN117391240B - Thermal power generation and coal combustion blending scheme optimization method based on similarity calculation - Google Patents

Thermal power generation and coal combustion blending scheme optimization method based on similarity calculation Download PDF

Info

Publication number
CN117391240B
CN117391240B CN202311286367.XA CN202311286367A CN117391240B CN 117391240 B CN117391240 B CN 117391240B CN 202311286367 A CN202311286367 A CN 202311286367A CN 117391240 B CN117391240 B CN 117391240B
Authority
CN
China
Prior art keywords
electricity
cost
unit
coal
operation condition
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202311286367.XA
Other languages
Chinese (zh)
Other versions
CN117391240A (en
Inventor
王辉
张依依
甘玮
韩金涛
王翔
薛泽彬
王晴
张子涵
费依蕃
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
North China Electric Power University
Original Assignee
North China Electric Power University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by North China Electric Power University filed Critical North China Electric Power University
Priority to CN202311286367.XA priority Critical patent/CN117391240B/en
Publication of CN117391240A publication Critical patent/CN117391240A/en
Application granted granted Critical
Publication of CN117391240B publication Critical patent/CN117391240B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/22—Matching criteria, e.g. proximity measures
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Economics (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Business, Economics & Management (AREA)
  • Health & Medical Sciences (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • General Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Quality & Reliability (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • Operations Research (AREA)
  • Primary Health Care (AREA)
  • Development Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Artificial Intelligence (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention relates to a thermal power generation and coal burning blending scheme optimizing method based on similarity calculation, and belongs to the technical field of thermal power generation. Comprising the following steps: acquiring a plurality of groups of historical operation condition data of a thermal power enterprise and a plurality of historical coal blending burning schemes at corresponding time points; respectively calculating the similarity between the current operation condition data and each group of historical operation condition data of the thermal power enterprise to obtain a plurality of similarity calculation results; determining an optimal similarity result set based on all similarity calculation results; determining a historical coal blending combustion scheme corresponding to each similarity calculation result in the optimal similarity result set as an alternative blending combustion scheme; calculating the electricity-measuring production cost of each alternative doping burning scheme to obtain an electricity-measuring production cost set; and determining an optimal coal blending combustion scheme of the current operation condition based on the electricity production cost set. The method can calculate the accurate electricity generation cost in real time, and can match the fire coal blending combustion scheme which can reach the optimal cost under the current working condition based on the historical data, thereby effectively guiding production.

Description

Thermal power generation and coal combustion blending scheme optimization method based on similarity calculation
Technical Field
The invention belongs to the technical field of thermal power generation, and particularly relates to a thermal power generation and coal burning blending scheme optimizing method based on similarity calculation.
Background
In the power generation process of a thermal power enterprise, different fire coal blending schemes can influence each link of the production process, so that different blending schemes can be caused to correspond to different power generation costs. Therefore, on the basis of ensuring the safe production of stable and full firing, a reasonable coal burning blending scheme is formulated, and the power generation cost of a thermal power enterprise can be effectively reduced.
In the research of the prior art, the calculation time period of the power generation cost has coarse granularity (usually calculated in months or days), and the calculation method usually considers insufficient factors, so that the real-time performance and accuracy of the cost calculation are not ideal, the real-time effective tuning and optimization of the coal burning blending scheme are influenced, and the aim of timely reducing the cost of a thermal power enterprise can not be met.
Disclosure of Invention
In view of the above analysis, the invention aims to provide a thermal power generation coal blending combustion scheme optimizing method based on similarity calculation, which specifically comprises the following steps:
Acquiring historical operation condition data, historical coal blending burning schemes and related parameters of a plurality of time points of a thermal power enterprise;
Respectively calculating the similarity between the current operation condition data of the thermal power enterprise and each group of the historical operation condition data to obtain a plurality of similarity calculation results;
determining an optimal similarity result set based on all the similarity calculation results;
Taking a historical coal blending combustion scheme corresponding to each time point in the optimal similarity result set as an alternative blending combustion scheme;
Calculating corresponding electricity-measuring production cost based on each alternative doping burning scheme and the related parameters of the corresponding time point to obtain electricity-measuring production cost set;
And determining an optimal coal blending combustion scheme of the current operation condition based on the electricity production cost set.
Further, the operating condition operation data includes: the unit output, the maximum output of the unit in the first 8 hours, the minimum output of the unit in the first 8 hours, the output average value of the unit in the first 8 hours, the output variance of the unit in the first 8 hours, the air temperature and the air humidity; the first 8 hours refer to the current working condition time point tracing back to the history for 8 hours.
Further, calculating the similarity between the current operation condition data and each set of the historical operation condition data by using a cosine similarity calculation method comprises the following steps:
obtaining a current operation condition vector based on the current operation condition data;
obtaining corresponding historical operating condition operating vectors based on each set of historical operating condition data;
And calculating cosine similarity for the current operation condition vector and each set of history operation condition vector to obtain similarity between the current operation condition data and each set of history operation condition data.
Further, the electricity-to-electricity production cost comprises electricity-to-electricity coal cost and electricity-to-electricity operation cost; the electricity-measuring operation cost comprises electricity-measuring electricity generation fuel cost, electricity-measuring water cost, electricity-measuring station electricity cost, electricity-measuring maintenance cost and electricity-measuring discharge cost.
Further, the calculation method of the electricity-measuring coal cost comprises the following steps:
Ccoal＝P1·Xsum·X1+P2·Xsum·X2+…+Pi·Xsum·Xi…+PN·Xsum·XN/(Lf×1000×T);
Wherein N represents the type of coal in the coal blending combustion scheme; x sum represents the total coal amount; x i represents the ratio of the ith coal in the total coal amount, an P i represents the price of the ith coal; l f is a power generation load; t is the power generation time.
Further, the calculation formula of the electricity maintenance cost is as follows:
Cwh＝(Ccl+Cjx)/Pmf；
Wherein, C wh represents the degree electricity maintenance cost, and the unit is Yuan/(kW.h); c cl represents the daily maintenance material cost of the current month, and the unit is yuan; c jx represents the current month overhaul fee, and the unit is yuan; p mf represents the total power generation in the month.
Further, the electricity water cost comprises brine water production cost, circulating water production cost and reclaimed water cost, and the brine water production cost, the circulating water production cost and the reclaimed water cost are respectively calculated based on the flow of the demineralized water outlet pipe per minute, the flow of the circulating water replenishing water per minute and the flow of the reclaimed water main pipe per minute.
Further, the electricity discharge cost comprises limestone powder cost, environment-friendly discharge cost and electricity urea cost; the calculation method of the limestone powder cost comprises the following steps:
Wherein, The converted concentration of the raw flue gas SO 2 is expressed in kg/Nm3; c yq,in represents the standard dry flow of the original flue gas, the unit is m 3/min;ηds represents the desulfurization efficiency, and the unit is; /(I)Represents the molar mass (100 g/mol) of CaCO 3; /(I)Represents the molar amount of SO 2 (64 g/mol); r Ca/S represents the calcium-sulfur ratio in units of; p lime represents the purity of limestone powder, the unit is; p shs represents the unit price of limestone powder in yuan per ton.
Further, the emission environmental protection cost comprises a power NO x pollution discharge tax, a power SO 2 pollution discharge tax and a power smoke dust pollution discharge tax, and the calculation methods are as follows:
Cpws＝Cyc,out×Cyq,out/Vyc×Pyc/(Lf×1000×T)；
Wherein, The concentration conversion of the NO x in the clean flue gas is expressed, the unit is kg/Nm 3;Cyq,out, the standard dry flow of the clean flue gas is expressed, and the unit is m 3/min; /(I)Represents NO x pollution equivalent value; the unit is m 3/min; /(I)Represents NO x unit price in yuan per kilogram; l f is the power generation load, and the unit is MW;
represents the conversion of the SO 2 concentration of the clean flue gas, and the unit is kg/Nm 3; /(I) The value of SO 2 pollution equivalent is expressed in m 3/min; /(I)Represents the unit price of SO 2 in yuan per kilogram;
C yc,out represents the conversion of the concentration of the clean flue gas and the smoke dust, the unit is kg/Nm 3;Vyc represents the equivalent value of the smoke dust pollution, the unit is m 3/min;Pyc represents the unit price of the smoke dust, and the unit is yuan/kg.
Further, the determining the optimal coal blending combustion scheme for the current operating condition based on the electricity-to-electricity production cost set comprises: and determining the fire coal blending combustion scheme corresponding to the lowest cost in the electricity-to-electricity production cost set as the optimal fire coal blending combustion scheme of the current operation working condition.
The invention can realize at least one of the following beneficial effects:
According to the method, historical data are used as a basis for supporting, and the thermal power generation production cost can be scientifically guided through matching the operation condition data: the method comprises the steps of acquiring a large number of minute-level historical operation condition data, acquiring a corresponding coal burning scheme and related parameters, performing similarity matching on the current operation condition data and the historical operation condition data to obtain a plurality of most similar historical operation conditions with the current operation condition, calculating the corresponding electricity-measuring production cost of each historical operation condition time point according to the plurality of most similar historical operation conditions, selecting the coal burning scheme with the lowest cost, optimizing the current production, and accurately guiding the thermal power generation production in real time.
By using the cosine similarity calculation method and reasonably selecting the operation condition data, the historical operation condition and the current operation condition are matched, the historical operation condition with the highest similarity with the current operation condition can be matched, and reasonable basis and data support are provided for the tuning and optimization of the coal blending combustion scheme.
In the cost calculation of the fire coal winding scheme corresponding to the historical operation working condition time point, the cost of pollutant treatment is calculated according to the flue gas flow rate and the pollutant concentration value by using a chemical formula based on the treatment principle of sulfur dioxide and nitrogen oxide, so that the emission cost accounting of the 'minute' level can be realized; compared with the prior art, the method for calculating the discharge cost of one test period or the discharge cost taking the month as the period has better instantaneity, and can realize the discharge cost accounting of the minute level.
In the cost calculation of the fire coal winding scheme corresponding to the historical operation working condition time point, the water consumption of each minute of a thermal power enterprise is calculated by collecting the circulating water, desalted water and the flow value of the reclaimed water pipeline, so that the electricity consumption cost of the 'minute' level is calculated, and the problems that the electricity consumption can only be calculated based on the power plant water taking coefficient in the prior art and the electricity consumption cost cannot be calculated in real time are solved.
In the cost calculation of the fire coal blending combustion scheme corresponding to the historical operation working condition time point, the electricity consumption water cost calculation and the electricity discharge cost calculation of the minute level are realized, the electricity generation fuel cost, the electricity plant consumption cost and the electricity maintenance cost of the corresponding time point are calculated, the minute level electricity operation cost can be calculated, and real-time effective data support is provided for the adjustment and optimization of the fire coal blending combustion scheme.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
Drawings
The drawings are only for purposes of illustrating particular embodiments and are not to be construed as limiting the invention, like reference numerals being used to designate like parts throughout the drawings;
FIG. 1 is a flow chart of the method of the present invention.
Detailed Description
The following detailed description of preferred embodiments of the application is made in connection with the accompanying drawings, which form a part hereof, and together with the description of the embodiments of the application, are used to explain the principles of the application and are not intended to limit the scope of the application.
Example 1
The invention discloses a thermal power generation and coal burning blending scheme optimizing method based on similarity calculation, which comprises the following steps:
Step S01, acquiring historical operation condition data, a historical coal blending burning scheme and related parameters of a thermal power enterprise at a plurality of time points;
step S02, calculating the similarity between the current operation condition data of the thermal power enterprise and each group of the historical operation condition data respectively to obtain a plurality of similarity calculation results;
Step S03, determining an optimal similarity result set based on all the similarity calculation results;
Step S04, taking a historical coal blending combustion scheme corresponding to each similarity calculation result in the optimal similarity result set as an alternative blending combustion scheme;
Step S05, calculating corresponding electricity-measuring production cost based on each alternative doping burning scheme and the related parameters of the corresponding time point, and obtaining an electricity-measuring production cost set;
and step S06, determining an optimal coal blending combustion scheme of the current operation condition based on the electricity production cost set.
Specifically, in step S01, the plurality of time points may be time points on the order of minutes; preferably, historical operation condition data, a historical coal blending combustion scheme and related parameters corresponding to each time point of the minute level in the last two months are selected.
Further, the operating condition operation data of each time point includes: the unit output, the maximum output of the unit in the first 8 hours, the minimum output of the unit in the first 8 hours, the output average value of the unit in the first 8 hours, the output variance of the unit in the first 8 hours, the air temperature and the air humidity; wherein the first 8 hours refer to the current working condition time point tracing back to the history for 8 hours; the air temperature and the air humidity refer to the current air temperature and the current air humidity of the thermal power plant. Exemplary, the operating data for a time point is: 194.399MW, 144.610MW, 167.908MW, 12.607, 8.40℃and 2.128g/kg.
Further, the unit output refers to the unit power as the unit power generation load.
Further, the coal blending combustion scheme comprises total coal quantity, coal types and blending proportion of various coals.
Further, the related parameters comprise the coal burning parameters (including the parameters of the above-mentioned coal burning blending burning scheme), the fuel oil parameters and water consumption parameters, the equipment electricity consumption parameters of the thermal power enterprises, the equipment maintenance parameters and the emission parameters of the thermal power enterprises; the details of the relevant parameters will be set forth in the following detailed description of step S05.
Specifically, in step S02, a cosine similarity calculation method is used to calculate the similarity between the current operating condition data and each set of the historical operating condition data, including:
obtaining a current operation condition vector based on the current operation condition data;
obtaining corresponding historical operating condition operating vectors based on each set of historical operating condition data;
And calculating cosine similarity for the current operation condition vector and each set of history operation condition vector to obtain similarity between the current operation condition data and each set of history operation condition data.
Further, defining the current operating condition vector based on the current operating condition data is expressed as:
α= [ ω 11,ω12,...,ω1n ]; wherein ω 1k (k ε [1, n) represents current operating condition operating data.
Defining the set of historical operating condition vectors based on each set of operating condition data is expressed as:
Beta= [ omega 21,ω22,...,ω2n ], where omega 2k (k e [1, n ]) represents historical operating condition operating data.
The current/historical operation condition vector is formed by sequentially arranging the current/historical operation condition data to form a one-dimensional vector serving as the current/historical operation condition vector.
The formula for calculating cosine similarity for the current operating condition vector and each set of historical operating condition operating vectors is:
Specifically, in step S03, an optimal similarity result set is determined based on all the similarity calculation results. Specifically, the top N similarity calculation result values of the numerical rank and the corresponding time points of the corresponding historical operation conditions are selected as the optimal similarity result set, and preferably, n=3.
By way of example, table 1 gives time points 2022-07-23 12:00 and each group of history working condition operation data, wherein three optimal similarity calculation results and corresponding history working condition operation time points are selected from calculation results of minute-level history working condition operation data two months before the current operation working condition time point to form an optimal similarity result set.
TABLE 1 optimal similarity result set example
Specifically, in step S04, the historical coal blending combustion scheme corresponding to each similarity calculation result in the optimal similarity result set is used as an alternative blending combustion scheme. For example, the historical coal blending combustion schemes corresponding to each time point in table 1 are alternative blending combustion schemes.
Specifically, in step S05, the electricity-generating production costs include electricity-generating coal costs and electricity-generating operation costs; the electricity-measuring operation cost of each time point comprises electricity-measuring electricity generation fuel cost, electricity-measuring water cost, electricity-measuring station electricity cost, electricity-measuring maintenance cost and electricity-measuring discharge cost.
Further, the calculation method of the electricity-to-fire coal cost and related fire coal parameters are as follows:
Ccoal＝P1·Xsum·X1+P2·Xsum·X2+…+Pi·Xsum·Xi…+PN·Xsum·XN/(Lf×1000×T);
Wherein N represents the type of coal in the coal blending combustion scheme; x sum represents the total coal amount, in units of: t is; x i represents the ratio of the ith coal in the total coal amount, an P i represents the price of the ith coal in units of: meta/t; l f is the power generation load, and the unit is MW; t is the power generation time, and the unit is h, and it is to be noted that T in all the formulas in the invention is the corresponding power generation time when each formula is calculated.
Further, the calculation formula of the fuel cost of the power generation and the related fuel parameters are as follows:
Cry＝Poil×Moil/(Lf×1000×T)；
Wherein, C ry represents the fuel cost of the electric power generation, and the unit is Yuan/(kW.h); p oil represents the fuel price in yuan/t; m oil represents the current fuel consumption, and the unit is t (ton); l f is the power generation load in MW.
Further, the electricity water consumption cost comprises a brine water production cost, a circulating water production cost and a reclaimed water cost, wherein the brine water production cost, the circulating water production cost and the reclaimed water cost are respectively calculated based on the flow of a demineralized water outlet pipe per minute, the flow of circulating water replenishing water per minute and the flow of a reclaimed water main pipe per minute.
Further, the brine water production cost, the circulating water production cost and the reclaimed water cost are respectively calculated by the following methods and related water consumption parameters:
Ccys＝Gcys×T×Pcys/(Lf×1000×T)；
Wherein, C cys represents the cost of producing water by the electric demineralized water, and the unit is Yuan/(kW.h); g cys represents the flow rate of the desalted water outlet mother pipe, and the unit is t/h; p cys represents the unit price of desalted water for preparing water, and the unit is yuan/t; l f is the power generation load, and the unit is MW;
Cxhs＝Gxhs×T×Pxhs/(Lf×1000×T)；
Wherein, C xhs represents the cost of producing water by the electric circulating water, and the unit is Yuan/(kW.h); g xhs represents the water supplementing flow rate of the circulating water, and the unit is t/h; p xhs represents the unit price of water production of circulating water, and the unit is yuan/t;
Czs＝Gzs×Pzs/(Lf×1000×T)；
Wherein, C zs represents the cost of the electric reclaimed water, and the unit is Yuan/(kW.h); g zs represents the flow rate of the water main pipe, and the unit is t/h; p zs represents the price per unit of water in yuan/t.
Specifically, the electricity consumption cost of the power plant comprises the equipment electricity consumption cost of all production links of a thermal power enterprise, and the calculation formula and the relevant equipment electricity consumption parameters are as follows:
Ccyd＝Fe×0.01×Pe；
Wherein, C cyd represents the electricity cost of the power plant, and the unit is Yuan/(kW.h); f e represents the station service power consumption; p e represents the power price on line, and the unit is Yuan/(kW.h).
Specifically, the calculation formula of the electricity maintenance cost and the related equipment maintenance parameters are as follows:
Cwh＝(Ccl+Cjx)/Pmf；
Wherein, C wh represents the degree electricity maintenance cost, and the unit is Yuan/(kW.h); c cl represents a daily maintenance material fee (replacement) in units of yuan for the month; c jx represents the current month overhaul fee, and the unit is yuan; p mf represents the total power generation in the month.
Further, the electricity discharge cost comprises limestone powder cost, environment-friendly discharge cost and electricity urea cost.
Specifically, the calculation method of limestone powder cost and the emission parameters of the related thermal power enterprises are as follows:
Wherein, The converted concentration of the raw flue gas SO 2 is expressed in kg/Nm3; c yq,in represents the standard dry flow of the original flue gas, the unit is m 3/min;ηds represents the desulfurization efficiency, and the unit is; /(I)Represents the molar mass (100 g/mol) of CaCO 3; /(I)Represents the molar amount of SO 2 (64 g/mol); r Ca/S represents the calcium-sulfur ratio in units of; p lime represents the purity of limestone powder, the unit is; p shs represents the unit price of limestone powder in yuan per ton.
Further, the emission environmental protection cost comprises a power NO x pollution discharge tax, a power SO 2 pollution discharge tax and a power smoke dust pollution discharge tax, and the emission parameters of the calculation method and related thermal power enterprises are respectively as follows:
Cpws＝Cyc,out×Cyq,out/Vyc×Pyc/(Lf×1000×T)；
Wherein, The concentration conversion of the clean flue gas NO x is expressed in kg/Nm3; c yq,out represents the dry flow of the clean flue gas, and the unit is m 3/min; /(I)Represents NO x pollution equivalent value; the unit is m 3/min; /(I)Represents NO x unit price in yuan per kilogram; l f is the power generation load, and the unit is MW;
represents the conversion of the SO 2 concentration of the clean flue gas, and the unit is kg/Nm 3; /(I) The value of SO 2 pollution equivalent is expressed in m 3/min; /(I)Represents the unit price of SO 2 in yuan per kilogram;
C yc,out represents the conversion of the concentration of the clean flue gas and the smoke dust, the unit is kg/Nm 3;Vyc represents the equivalent value of the smoke dust pollution, the unit is m 3/min;Pyc represents the unit price of the smoke dust, and the unit is yuan/kg;
The flue gas standard dry flow is the flow obtained by converting the actually measured flue gas flow into the standard state (0 ℃ =273K, 101.325 Kpa), and the calculation method of the original flue gas standard dry flow and the net flue gas standard dry flow and the emission parameters of the thermal power enterprises are respectively as follows:
Wherein C yq,in represents the standard dry flow of the raw flue gas, the unit is m 3/min;Syyq represents the cross section area of a flue gas inlet, the unit is m 2;Vyyq represents the flow rate of the raw flue gas, and the unit is m/s; f a represents the atmospheric pressure in Pa; f yyq represents the static pressure of the raw flue gas at the measuring point, and the unit is Pa; The unit of the moisture content in the original flue gas is the volume percentage; t yyq represents the temperature of raw flue gas, and the unit is DEG C; c yq,out represents the dry flow of the clean flue gas, the unit is m 3/min;Sjyq represents the cross-sectional area of a flue gas outlet, the unit is m 2;Vjyq represents the flow rate of the clean flue gas, and the unit is m/s; f jyq represents the net smoke static pressure at the measuring point, and the unit is Pa; /(I) The unit of the water content volume percentage in the clean flue gas is shown as follows; t jyq represents the net flue gas temperature in degrees Celsius.
Specifically, the calculation method of the electricity-to-electricity urea cost and the emission parameters of the related thermal power enterprises are as follows:
Wherein, C ns represents the electricity urea cost, and the unit cell/(kW.h); The ammonia amount (ammonia supply flow rate) is expressed in kg/h; p ns represents the urea unit price in yuan/kg.
Further, the electricity-to-electricity running cost is the sum of the electricity-to-electricity power generation fuel cost, electricity-to-electricity water cost, electricity-to-electricity plant electricity cost, electricity-to-electricity maintenance cost and electricity-to-electricity discharge cost, and is expressed as:
Cop＝Cry+Cys+Ccyd+Cwh+Cpf；
Wherein, C ys represents the cost of electricity water, C ys＝Ccys+Cxhs+Czs;
C pf represents the electricity emission cost, C pf＝Cshs+Chbs+Cns;
c hbs represents the cost of the environmental protection of the emissions,
Further, the electricity-to-electricity production cost is the sum of electricity-to-electricity coal cost and electricity-to-electricity operation cost, and is expressed as:
C＝Ccoal+Cop。
for example, table 2 shows the calculated corresponding electricity production costs for each time point corresponding to each history in table 1.
TABLE 2 electric production costs for each respective degree calculated based on the optimal similarity result set
Specifically, in step S06, an optimal coal blending combustion scheme for the current operating condition is determined based on the electricity production cost set. For example, according to the table 2, in the optimal similarity result set, the lowest value in the electricity production cost is 0.3906 yuan, so that the historical fire coal blending combustion scheme corresponding to the time point is selected as the optimal fire coal blending combustion scheme of the current operation condition; illustratively, the historical coal blending firing schedule includes: the total coal amount is 92.22t, and the types of coal comprise high-heat high-volatile coal, high-heat low-volatile coal and economic coal, wherein the high-heat high-volatile coal accounts for 43.31%, the high-heat low-volatile coal accounts for 10.07% and the economic coal accounts for 46.62%.
The embodiment discloses a thermal power generation fire coal blending combustion scheme optimization method based on similarity calculation, which uses historical data in the operation process of a thermal power generation enterprise as a basis support, scientifically guides thermal power generation production to reduce cost through matching operation condition data, and comprises the following steps: the method comprises the steps of acquiring a large number of minute-level historical operation condition data, acquiring a corresponding coal burning scheme and related parameters, performing similarity matching on the current operation condition data and the historical operation condition data to obtain a plurality of most similar historical operation conditions with the current operation condition, calculating the corresponding electricity-measuring production cost of each historical operation condition time point according to the plurality of most similar historical operation conditions, selecting the coal burning scheme with the lowest cost, optimizing the current production, and accurately guiding the thermal power generation production in real time.
By using the cosine similarity calculation method and reasonably selecting the operation condition data, the historical operation condition and the current operation condition are matched, the historical operation condition with the highest similarity with the current operation condition can be matched, and reasonable basis and data support are provided for the tuning and optimization of the coal blending combustion scheme.
In the cost calculation of the fire coal winding scheme corresponding to the historical operation working condition time point, the cost of pollutant treatment is calculated according to the flue gas flow rate and the pollutant concentration value by using a chemical formula based on the treatment principle of sulfur dioxide and nitrogen oxide, so that the emission cost accounting of the 'minute' level can be realized; compared with the prior art for realizing the discharge cost calculation of one test period or the discharge cost calculation method taking 'month' as the period, the method has the real-time performance for realizing the discharge cost accounting of 'minute' level.
In the cost calculation of the fire coal winding scheme corresponding to the historical operation working condition time point, the water consumption of each minute of a thermal power enterprise is calculated by collecting the circulating water, desalted water and the flow value of the reclaimed water pipeline, so that the electricity consumption cost of the 'minute' level is calculated, and the problems that the electricity consumption can only be calculated based on the power plant water taking coefficient in the prior art and the electricity consumption cost cannot be calculated in real time are solved.
In the cost calculation of the fire coal winding scheme corresponding to the historical operation working condition time point, the electricity consumption water cost calculation and the electricity discharge cost calculation of the minute level are realized, the electricity generation fuel cost, the electricity plant electricity consumption cost and the electricity maintenance cost of the corresponding time point are calculated, the minute level electricity operation cost can be calculated, and real-time effective data support is provided for the adjustment and optimization of the fire coal blending scheme.
The present invention is not limited to the above-mentioned embodiments, and any changes or substitutions that can be easily understood by those skilled in the art within the technical scope of the present invention are intended to be included in the scope of the present invention.

Claims (5)

1. The thermal power generation and coal burning blending scheme optimizing method based on similarity calculation is characterized by comprising the following steps of:
Acquiring historical operation condition data, a historical coal blending burning scheme and related parameters of a plurality of time points of a thermal power enterprise; the operating condition data includes: the unit output, the maximum output of the unit in the first 8 hours, the minimum output of the unit in the first 8 hours, the output average value of the unit in the first 8 hours, the output variance of the unit in the first 8 hours, the air temperature and the air humidity; wherein the first 8 hours refer to the current working condition time point tracing back to the history for 8 hours; the plurality of time points are time points on the order of minutes;
Respectively calculating the similarity between the current operation condition data of the thermal power enterprise and each group of the historical operation condition data to obtain a plurality of similarity calculation results;
determining an optimal similarity result set based on all the similarity calculation results;
Taking a historical coal blending combustion scheme corresponding to each time point in the optimal similarity result set as an alternative blending combustion scheme;
Calculating corresponding electricity-measuring production cost based on each alternative doping burning scheme and the related parameters of the corresponding time point to obtain electricity-measuring production cost set; the electricity-to-electricity production cost comprises electricity-to-electricity coal cost and electricity-to-electricity operation cost; the electricity-measuring operation cost comprises electricity-measuring electricity generation fuel cost, electricity-measuring water cost, electricity-measuring station electricity cost, electricity-measuring maintenance cost and electricity-measuring discharge cost; the electricity water cost comprises brine water production cost, circulating water production cost and reclaimed water cost, and the brine water production cost, the circulating water production cost and the reclaimed water cost are respectively calculated based on the flow of a demineralized water outlet pipe per minute, the flow of circulating water replenishing water per minute and the flow of a reclaimed water main pipe per minute; the electricity-metering discharge cost comprises limestone powder cost, environment-friendly discharge cost and electricity-metering urea cost; the calculation method of the limestone powder cost comprises the following steps: Wherein/> The concentration of the converted raw flue gas SO 2 is expressed, the unit is kg/Nm 3;Cyq,in, the standard dry flow of the raw flue gas is expressed, the unit is m 3/min;ηds, the desulfurization efficiency is expressed; /(I)Represents the molar mass (100 g/mol) of CaCO 3; /(I)Represents the molar amount of SO 2 (64 g/mol); r Ca/S represents the calcium-sulfur ratio in units of; p lime represents the purity of limestone powder, the unit is; p shs represents the unit price of limestone powder in yuan per ton; l f is a power generation load; t is the power generation time; the emission environmental protection cost comprises a power NO x pollution discharge tax, a power SO 2 pollution discharge tax and a power smoke dust pollution discharge tax, and the calculation methods are as follows:
C pws＝Cyc,out×Cyq,out/Vyc×Pyc/(Lf ×1000×T); wherein, The concentration conversion of the NO x in the clean flue gas is expressed, the unit is kg/Nm 3;Cyq,out, the standard dry flow of the clean flue gas is expressed, and the unit is m 3/min; /(I)Represents NO x pollution equivalent value; the unit is m 3/min; /(I)Represents NO x unit price in yuan per kilogram; l f is the power generation load, and the unit is MW; /(I)Represents the conversion of the SO 2 concentration of the clean flue gas, and the unit is kg/Nm 3; /(I)The value of SO 2 pollution equivalent is expressed in m 3/min; /(I)Represents the unit price of SO 2 in yuan per kilogram; c yc,out represents the conversion of the concentration of the clean flue gas and the smoke dust, the unit is kg/Nm 3;Vyc represents the equivalent value of the smoke dust pollution, the unit is m 3/min;Pyc represents the unit price of the smoke dust, and the unit is yuan/kg;
And determining an optimal coal blending combustion scheme of the current operation condition based on the electricity production cost set.
2. The method of claim 1, wherein calculating the similarity between the current operating condition data and each set of the historical operating condition data using a cosine similarity calculation method comprises:
obtaining a current operation condition vector based on the current operation condition data;
obtaining corresponding historical operating condition operating vectors based on each set of historical operating condition data;
and calculating cosine similarity for the current operation condition vector and each group of the history operation condition vectors to obtain similarity of the current operation condition data and each group of the history operation condition data.
3. The method for optimizing the blending combustion scheme of the fire coal according to claim 2, wherein the calculation method of the electricity-measuring fire coal cost is as follows:
Ccoal＝P1·Xsum·X1+P2·Xsum·X2+…+Pi·Xsum·Xi…+PN·Xsum·XN/(Lf×1000×T);
Wherein N represents the type of coal in the coal blending combustion scheme; x sum represents the total coal amount; x i represents the ratio of the ith coal in the total coal amount, an P i represents the price of the ith coal; l f is a power generation load; t is the power generation time.
4. The method for optimizing a fire coal blending combustion scheme according to claim 3, wherein the calculation formula of the electricity-to-electricity maintenance cost is as follows:
Cwh＝(Ccl+Cjx)/Pmf；
Wherein, C wh represents the degree electricity maintenance cost, and the unit is Yuan/(kW.h); c cl represents the daily maintenance material cost of the current month, and the unit is yuan; c jx represents the current month overhaul fee, and the unit is yuan; p mf represents the total power generation in the month.
5. The method of claim 4, wherein determining an optimal coal blending combustion scheme for the current operating condition based on the set of electrical production costs comprises: and determining the fire coal blending combustion scheme corresponding to the lowest cost in the electricity-to-electricity production cost set as the optimal fire coal blending combustion scheme of the current operation working condition.
CN202311286367.XA 2023-10-07 2023-10-07 Thermal power generation and coal combustion blending scheme optimization method based on similarity calculation Active CN117391240B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311286367.XA CN117391240B (en) 2023-10-07 2023-10-07 Thermal power generation and coal combustion blending scheme optimization method based on similarity calculation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311286367.XA CN117391240B (en) 2023-10-07 2023-10-07 Thermal power generation and coal combustion blending scheme optimization method based on similarity calculation

Publications (2)

Publication Number Publication Date
CN117391240A CN117391240A (en) 2024-01-12
CN117391240B true CN117391240B (en) 2024-05-03

Family

ID=89435296

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311286367.XA Active CN117391240B (en) 2023-10-07 2023-10-07 Thermal power generation and coal combustion blending scheme optimization method based on similarity calculation

Country Status (1)

Country Link
CN (1) CN117391240B (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105069540A (en) * 2015-09-02 2015-11-18 华中科技大学 Method for designing fuel procurement decision support system for thermal power enterprises
CN105117808A (en) * 2015-09-22 2015-12-02 华润电力登封有限公司 Coal blending combustion optimizing method
CN106247395A (en) * 2016-06-30 2016-12-21 华润电力登封有限公司 A kind of feeder control method
CN107274027A (en) * 2017-06-22 2017-10-20 湖南华润电力鲤鱼江有限公司 A kind of many coal coal mixing combustion optimization methods of coal unit
CN107301477A (en) * 2017-06-22 2017-10-27 湖南华润电力鲤鱼江有限公司 A kind of coal-fired procurement decisions method based on many coal coal mixing combustion optimizing models
CN110111030A (en) * 2019-06-18 2019-08-09 华电国际电力股份有限公司技术服务分公司 A kind of coal mixing combustion evaluation method and system
CN115496280A (en) * 2022-09-20 2022-12-20 华能国际电力江苏能源开发有限公司 Coal yard blending combustion data analysis method and system
CN116739630A (en) * 2023-05-06 2023-09-12 华能临沂发电有限公司 Real-time cost accounting system for fire coal of power plant

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105069540A (en) * 2015-09-02 2015-11-18 华中科技大学 Method for designing fuel procurement decision support system for thermal power enterprises
CN105117808A (en) * 2015-09-22 2015-12-02 华润电力登封有限公司 Coal blending combustion optimizing method
CN106247395A (en) * 2016-06-30 2016-12-21 华润电力登封有限公司 A kind of feeder control method
CN107274027A (en) * 2017-06-22 2017-10-20 湖南华润电力鲤鱼江有限公司 A kind of many coal coal mixing combustion optimization methods of coal unit
CN107301477A (en) * 2017-06-22 2017-10-27 湖南华润电力鲤鱼江有限公司 A kind of coal-fired procurement decisions method based on many coal coal mixing combustion optimizing models
CN110111030A (en) * 2019-06-18 2019-08-09 华电国际电力股份有限公司技术服务分公司 A kind of coal mixing combustion evaluation method and system
CN115496280A (en) * 2022-09-20 2022-12-20 华能国际电力江苏能源开发有限公司 Coal yard blending combustion data analysis method and system
CN116739630A (en) * 2023-05-06 2023-09-12 华能临沂发电有限公司 Real-time cost accounting system for fire coal of power plant

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
火电厂发电污染成本核算及敏感度分析;闫顺林;谷兵;艾书剑;田东旭;;电力科学与工程;20160728(07);全文 *
面向电力市场竞价的燃煤电厂实时发电成本计算模型;李林;周大为;;湖北电力;20171228(12);全文 *

Also Published As

Publication number Publication date
CN117391240A (en) 2024-01-12

Similar Documents

Publication Publication Date Title
CN108173282B (en) A kind of consideration electricity turns gas operating cost integrated energy system Optimization Scheduling
Wolsky et al. Recovering CO2 from large− and medium-size stationary combustors
CN111799777A (en) Comprehensive energy planning method considering coupling of natural gas and electric power
CN107276080A (en) A kind of computational methods of energy-saving and emission-reduction developing goal contribution rate
CN114548584A (en) Optimization method of comprehensive energy system containing electric conversion gas and carbon capture equipment
CN114836776B (en) Economic evaluation method for full life cycle of new energy coupling coal chemical industry multi-energy system
CN116247719A (en) Micro-grid two-stage robust optimal configuration method based on ladder carbon transaction
CN114925328A (en) Coal blending and burning calculation method and system for coal-fired unit
CN117391240B (en) Thermal power generation and coal combustion blending scheme optimization method based on similarity calculation
CN114444785A (en) Deep scheduling method and system for deep peak shaving of coal-fired unit
CN113240205A (en) Regional energy utilization system substitution optimization method based on multi-energy comprehensive utilization
Fernando et al. Techno-economic analysis of a small-scale power-to-green urea plant
CN112418664B (en) Particle swarm optimization-based binning combination blending method and system
CN117010728B (en) Comprehensive power generation cost optimization method for thermal power enterprises
CN112862632B (en) Method and system for blending and burning coal in thermal power plant
CN105823071B (en) A kind of method for reducing boiler ultra-clean discharge operation cost
CN117391734B (en) Power generation cost prediction method based on support vector machine regression model
CN114818269A (en) Industrial park comprehensive energy system and factory production plan collaborative optimization method
CN117575658A (en) Coal-fired unit depth peak shaving cost prediction method based on life cycle and environmental protection compensation
CN116379464B (en) Automatic optimizing method for total cost of NOx under full load of coal-fired unit
CN110894980A (en) Economical evaluation method based on heat accumulating type electric heating load
CN115907362A (en) Industrial park multifunctional complementary system optimization scheduling model construction method
Kavuma et al. Reducing Carbon Dioxide Emission Generated by Thermal Power Plants with Solar Power on the Uganda's Electricity Grid
CN114754353B (en) Circulating fluidized bed boiler combustion optimization method integrating neighborhood rough set machine learning
Suo et al. Analysis of carbon verification calculation methods and emission reduction measures for thermal power plants

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant