CN114880919A - Method for calculating optimal furnace internal and external desulfurization proportion of circulating fluidized bed unit - Google Patents

Method for calculating optimal furnace internal and external desulfurization proportion of circulating fluidized bed unit Download PDF

Info

Publication number
CN114880919A
CN114880919A CN202210298473.9A CN202210298473A CN114880919A CN 114880919 A CN114880919 A CN 114880919A CN 202210298473 A CN202210298473 A CN 202210298473A CN 114880919 A CN114880919 A CN 114880919A
Authority
CN
China
Prior art keywords
furnace
desulfurization
formula
cost
limestone
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.)
Pending
Application number
CN202210298473.9A
Other languages
Chinese (zh)
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.)
Tsinghua University
North China Electric Power University
Huaneng Group Technology Innovation Center Co Ltd
Original Assignee
Tsinghua University
North China Electric Power University
Huaneng Group Technology Innovation Center Co Ltd
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 Tsinghua University, North China Electric Power University, Huaneng Group Technology Innovation Center Co Ltd filed Critical Tsinghua University
Priority to CN202210298473.9A priority Critical patent/CN114880919A/en
Publication of CN114880919A publication Critical patent/CN114880919A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/27Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/12Computing arrangements based on biological models using genetic models
    • G06N3/126Evolutionary algorithms, e.g. genetic algorithms or genetic programming
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/10Analysis or design of chemical reactions, syntheses or processes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/08Probabilistic or stochastic CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Evolutionary Computation (AREA)
  • Biophysics (AREA)
  • Health & Medical Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Artificial Intelligence (AREA)
  • Chemical & Material Sciences (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Biology (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Analytical Chemistry (AREA)
  • Physiology (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Genetics & Genomics (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Biomedical Technology (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Mathematical Physics (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Treating Waste Gases (AREA)

Abstract

The invention discloses a method for calculating the optimal furnace internal and external desulfurization proportion of a circulating fluidized bed unit, which comprises the steps of firstly establishing a CFB unit furnace internal and external desulfurization comprehensive cost model, and selecting load, coal quality, coal feeding quantity, total air quantity, bed temperature, rated power of a limestone conveying fan, rated power of a slurry circulating pump, feeding flow rate of limestone in a furnace, ammonia injection quantity and raw flue gasSO 2 Concentration, clean flue gas SO 2 Concentration as an input variable; determining relational expressions of the desulfurization efficiency in the furnace and load, the molar ratio of calcium to sulfur, the bed temperature and the air-coal ratio when establishing the comprehensive cost model; secondly, determining the generation concentration of SO2 by using the input of the comprehensive cost model selected in the step 1, and establishing an in-furnace and out-furnace desulphurization comprehensive cost model under the typical load working condition of the CFB unit after assuming the in-furnace desulphurization proportion; and finally, solving the optimal in-furnace desulfurization proportion under the typical load working condition by using an intelligent optimization algorithm, and fitting to obtain the optimal in-furnace desulfurization proportion under each load working condition.

Description

Method for calculating optimal furnace internal and external desulfurization proportion of circulating fluidized bed unit
Technical Field
The invention belongs to the field of pollutant control optimization of thermal power generating units, and relates to a method for calculating an optimal furnace internal and external desulfurization proportion of a circulating fluidized bed unit.
Background
The Circulating Fluidized Bed (CFB) unit has the advantages of strong combustion stability and low pollutant treatment cost, can reach the standard of low emission through in-furnace desulfurization, and has simple process flow. With the ultralow emission standard of pollutant emission of coal-fired power generating units in China, in order to realize ultralow emission, the circulating fluidized bed unit is additionally provided with external desulfurization (flue gas desulfurization) equipment, and ultralow emission, SO, is realized by a combined desulfurization mode inside and outside the furnace 2 The ultra-low emission standard is less than or equal to 35mg/m 3 . At present, the research on the optimal furnace internal and external desulfurization proportion of the CFB unit is less, and the field operation is lack of guidance. On-site operators often ensure the original SO of the flue gas at the outlet of the hearth by adjusting the calcium-sulfur ratio of the desulfurization in the furnace 2 The concentration is in a certain fixed range, SO is ensured after the external desulfurization of the furnace 2 The emission concentration is less than 35mg/m 3 . Although the operation mode is simple, the change rules of the in-furnace/out-furnace desulfurization efficiency and the operation cost under different load working conditions of the unit are not considered, so that the desulfurization material consumption of the CFB unit is increased, the desulfurization cost is increased, and the operation economy of the unit is reduced.
Object of the Invention
The invention aims to solve the problems of poor distribution of the internal and external desulfurization proportion of the CFB unit furnace and poor desulfurization operation economy in the prior art, and provides a method for calculating the optimal internal and external desulfurization proportion of the circulating fluidized bed unit furnace.
Disclosure of Invention
The invention provides a method for calculating the optimal furnace internal and external desulfurization proportion of a circulating fluidized bed unit, which comprises the following steps:
step 1, establishing a CFB unit furnace internal and external desulfurization comprehensive cost model, and selecting load, coal quality, coal feeding quantity, total air quantity, bed temperature, limestone conveying fan rated power, slurry circulating pump rated power, limestone feeding flow rate in a furnace, ammonia injection quantity and raw flue gas SO 2 Concentration, clean flue gas SO 2 Concentration is used as an input variable of the comprehensive cost model;
step 2, determining relational expressions of desulfurization efficiency in the furnace and load, calcium-sulfur molar ratio, bed temperature and air-coal ratio when establishing the comprehensive cost model;
step 3, determining the SO2 generation concentration by using the input of the comprehensive cost model selected in the step 1, and establishing an in-furnace and out-furnace desulfurization comprehensive cost model under the typical load working condition of the CFB unit after assuming the in-furnace desulfurization proportion;
and 4, solving the optimal in-furnace desulfurization proportion under the typical load working condition by using an intelligent optimization algorithm, and fitting to obtain the optimal in-furnace desulfurization proportion under each load working condition.
Preferably, the process of establishing the comprehensive cost model of the internal and external desulfurization of the CFB unit in step 1 is as follows: defining the total combined desulfurization cost inside and outside the furnace as comprising: the comprehensive cost is the sum of the total desulfurization cost minus the benefit generated by gypsum; wherein the limestone consumption cost comprises the consumption cost of the in-furnace desulfurized limestone and the consumption cost of the out-furnace desulfurized limestone; the power consumption cost of the equipment comprises the power consumption of a limestone conveying fan in the furnace and the power consumption of a slurry circulating pump outside the furnace; the denitration cost in the furnace is increased along with the increase of the molar ratio of the desulfurized calcium to the desulfurized calcium in the furnace and the denitration cost in the furnace is increased; when the desulfurization efficiency in the furnace is unchanged, the lower the molar ratio of calcium to sulfur is, the lower the heat loss cost is, and the lower the converted coal consumption is; the dosage of the in-furnace desulfurization limestone is determined by the in-furnace desulfurization calcium-sulfur molar ratio, the in-furnace desulfurization efficiency and the coal quality, and the in-furnace desulfurization efficiency is related to load, the calcium-sulfur molar ratio, bed temperature and air-coal ratio; the amount of limestone desulfurized outside the furnace is related to the desulfurization efficiency outside the furnace;
efficiency of desulfurization in the furnace
Figure BDA0003562524650000031
Expressed as shown in formula (1):
Figure BDA0003562524650000032
in the formula (1), W c For coal supply, S ar Is a sulfur content, A ir Total air volume, k f For the dimensionless conversion coefficient of the smoke,
Figure BDA0003562524650000033
is the raw flue gas SO under standard conditions 2 Concentration;
efficiency of external desulfurization of said furnace
Figure BDA0003562524650000034
Expressed as shown in formula (2):
Figure BDA0003562524650000035
in the formula (2), under the standard condition
Figure BDA0003562524650000036
Is clean flue gas SO 2 And (4) concentration.
Preferably, in the step 2, a least square method is adopted to fit a relational expression of the in-furnace desulfurization efficiency, the load, the calcium-sulfur molar ratio and the bed temperature under a typical load working condition, and the coefficient of the relational expression is corrected by using the air-coal ratio, and the method specifically comprises the following substeps:
substep S21, selecting bed temperature, coal supply quantity, total air quantity, coal quality and raw flue gas SO under typical load working condition 2 Concentration, in-furnace limestone feed flow rate;
and a substep S22, calculating the in-furnace desulfurization efficiency, the air-coal ratio and the calcium-sulfur molar ratio under the typical load working condition, wherein the calculation of the in-furnace desulfurization calcium-sulfur molar ratio is shown by an equation (3):
Figure BDA0003562524650000037
in the formula (3), the reaction mixture is,
Figure BDA0003562524650000038
in order to obtain the purity of the limestone,
Figure BDA0003562524650000039
a feed flow rate for limestone in the furnace;
and a substep S23 of fitting a relational expression of the desulfurization efficiency in the furnace, the molar ratio of calcium to sulfur and the bed temperature under the typical load working condition by adopting a least square method, wherein the relational expression is expressed as a formula (4):
Figure BDA00035625246500000310
in the formula (4), A is a dimensionless coefficient, B is a function of bed temperature, m CFB The molar ratio of calcium to sulfur in the furnace is shown;
and a substep S24, correcting the A by using the air-coal ratio under the typical load working condition according to the operation data, and finally obtaining a relational expression of the in-furnace desulfurization efficiency, the calcium-sulfur molar ratio and the bed temperature under the typical load working condition.
Preferably, in step 3, after assuming the ratio of desulfurization in the furnace, determining the power consumption of a limestone conveying fan, the power consumption of a slurry circulating pump, and the limestone consumption of desulfurization inside and outside the furnace, wherein the power consumption of the limestone conveying fan is determined by the feeding flow rate of limestone in the furnace and the rated power of the limestone conveying fan, and the power consumption of the slurry circulating pump is determined by the number of circulating pumps in operation and the rated power, specifically comprising the following substeps:
substep S31, fitting the external desulfurization efficiency and the raw flue gas SO by a least square method 2 The relationship between the concentration and the total air volume is expressed by the following formula (6):
Figure BDA0003562524650000041
in the formula (6), a1, a2 and a3 are model coefficients;
and a substep S32, assuming that the desulfurization ratio inside the furnace is x under a certain typical load working condition, the desulfurization ratio outside the furnace is 1-x, and combining the formula (1) and the formula (4), determining the molar ratio m of calcium and sulfur desulfurized inside the furnace CFB
Substep S33, determining the limestone consumption cost W in the furnace 1 Expressed as shown in formula (7):
Figure BDA0003562524650000042
in the formula (7), MCaCO3 and MCaO are the relative molecular masses of CaCO3 and CaO respectively, and the unit is g/mol; u1 is the unit price of limestone in the furnace, and the unit is Yuan/kg;
the heat loss cost W 2 Expressed as shown in formula (8):
W 2 =W c [η/(η-Δη)-1]u 2 (8),
in the formula (8), η is the boiler design efficiency; Δ η is boiler heat loss; u2 is the unit price of the fire coal, and the unit is Yuan/kg;
the denitration cost W 3 Expressed as shown in formula (9):
W 3 =k 1 m CFB W c u 3 (9),
in the formula (9), k1 is a cost coefficient; u3 is the urea unit price, unit is yuan/kg;
electrically connecting the conveying fanCost W 4 Expressed as shown in formula (10):
Figure BDA0003562524650000051
in the formula (10), α is a compressed air coefficient; u. of 4 The unit is yuan/kWh for the price of the power on the internet;
the dosage cost W of the limestone outside the furnace 5 Expressed as shown in formula (11):
Figure BDA0003562524650000052
in the formula (11), mCFB and wet are in the molar ratio of calcium to sulfur outside the furnace; u5 is the unit price of limestone outside the furnace, and the unit is Yuan/t.
The power consumption W6 of the circulating pump is expressed as shown in the formula (12):
W 6 =nP w u 4 (12),
in the formula (12), n is the number of the slurry circulating pumps started and is the raw flue gas SO 2 Determining concentration and load; p w The power of a single slurry circulating pump is in kW;
the gypsum income V 7 Expressed as shown in formula (13):
Figure BDA0003562524650000053
wherein eta (H) 2 O) is the water content of the gypsum; u. of 8 Is the unit price of gypsum, and the unit is yuan/kg;
and (3) expressing the comprehensive cost f (x) of the combined desulfurization inside and outside the furnace as shown in a formula (14):
f(x)=W 1 +W 2 +W 3 +W 4 +W 5 +W 6 +W 7 -V 8 0≤x≤x max (14),
in the formula (14), W 1 、W 2 、W 3 、W 4 、W 5 、W 6 、V 7 The amount cost of the limestone in the furnace, the cost of heat loss, the cost of denitration, the cost of power consumption of a conveying fan, the cost of the amount of the limestone outside the furnace, the power consumption of a circulating pump and the benefit of gypsum are respectively determined by the desulfurization proportion x in the furnace; x is the number of max Determined by the in-furnace desulfurization capacity of the CFB unit.
Preferably, in step 4, the optimal in-furnace desulfurization proportion under the typical load condition is a solution obtained by the genetic algorithm optimizing the minimum integrated cost of inside and outside desulfurization of the furnace, wherein the genetic algorithm optimizing process comprises the following substeps:
substep S41, encoding: selecting unsigned binary integers to represent individuals x i
Substep S42, generating an initial population: randomly generating N individuals as an initial population, and setting the iteration number as N;
substep S43, fitness calculation: using the value of the integrated cost function g (x) i) As an individual x i Selecting a fitness function as shown in formula (15):
g(x i )=minf(x i ) (15);
substep S44, selection, crossover, mutation operation: the individuals with higher fitness in the current group are inherited to the next generation; the method adopts a single-point crossing method to carry out crossing operation, adopts a basic bit variation method to carry out variation operation, saves the next generation of group, and increases the iteration times by 1;
substep S45, termination condition judgment: if the iteration times are more than or equal to n, stopping the calculation, and outputting the obtained individual with the minimum fitness as an optimal solution, namely the optimal furnace internal and external desulfurization proportion; otherwise, return to substep S43.
Further preferably, the number of population N is 20, the number of termination iterations N is 80, the crossover probability is 0.4, and the mutation probability is 0.001.
Drawings
FIG. 1 is a schematic flow chart of the method for calculating the optimal furnace internal and external desulfurization proportion of the circulating fluidized bed unit.
Detailed Description
The following detailed description of the embodiments of the invention refers to the accompanying drawings. It will be understood by those skilled in the art that the present description is illustrative of the preferred embodiments of the invention and is not to be construed as limiting the scope of the invention in any way, and that variations or modifications can be made without departing from the spirit and scope of the invention.
The invention discloses a method for calculating the optimal furnace internal and external desulfurization proportion of a circulating fluidized bed unit, and fig. 1 is a flow schematic diagram of the method for calculating the optimal furnace internal and external desulfurization proportion of the circulating fluidized bed unit, as shown in fig. 1, the method comprises the following steps:
step 1, establishing a CFB unit furnace internal and external desulfurization comprehensive cost model, and selecting load, coal quality, coal feeding quantity, total air quantity, bed temperature, limestone conveying fan rated power, slurry circulating pump rated power, limestone feeding flow rate in a furnace, ammonia injection quantity and raw flue gas SO 2 Concentration, clean flue gas SO 2 Concentration is used as an input of the comprehensive cost model;
step 2, determining relational expressions of desulfurization efficiency in the furnace and load, calcium-sulfur molar ratio, bed temperature and air-coal ratio when establishing the comprehensive cost model;
step 3, determining the SO2 generation concentration by using the input of the comprehensive cost model selected in the step 1, and establishing an in-furnace and out-furnace desulfurization comprehensive cost model under the typical load working condition of the CFB unit after assuming the in-furnace desulfurization proportion;
and 4, solving the optimal in-furnace desulfurization proportion under the typical load working condition by using an intelligent optimization algorithm, and fitting to obtain the optimal in-furnace desulfurization proportion under each load working condition.
In this embodiment, in step 2, a least square method is used to fit a relational expression of the in-furnace desulfurization efficiency, the load, the calcium-sulfur molar ratio, and the bed temperature under typical load conditions, and the coefficients of the relational expression are corrected by using the air-coal ratio. And 3, after the desulfurization proportion in the furnace is assumed, determining the power consumption of a limestone conveying fan, the power consumption of a slurry circulating pump and the limestone consumption of desulfurization inside and outside the furnace, wherein the power consumption of the limestone conveying fan is determined by the feeding flow rate of limestone in the furnace and the rated power of the limestone conveying fan, and the power consumption of the slurry circulating pump is determined by the number of circulating pumps in operation and the rated power. In step 4, the optimal in-furnace desulfurization proportion under the typical load working condition is the solution of the minimum comprehensive cost of in-furnace and out-furnace desulfurization optimized by the genetic algorithm.
The method comprises the following specific steps:
1. selection of model input variables
The circulating fluidized bed unit realizes ultralow emission by a combined desulfurization operation mode inside and outside the furnace. The total cost of the combined desulfurization inside and outside the furnace comprises: the power consumption cost of equipment, the denitration cost in the furnace, the heat loss cost and the consumption cost of limestone, and the comprehensive cost is the total desulfurization cost minus the benefit generated by gypsum. Wherein: the limestone consumption cost comprises the consumption cost of the desulfurized limestone inside the furnace and the consumption cost of the desulfurized limestone outside the furnace. The equipment power consumption cost comprises the following steps: the power consumption of a limestone conveying fan in the furnace and the power consumption of a slurry circulating pump outside the furnace; the denitration cost in the furnace is influenced by the molar ratio of calcium to sulfur of desulfurization in the furnace, and the denitration cost in the furnace is increased along with the increase of the molar ratio of calcium to sulfur; the heat loss cost is also influenced by the molar ratio of calcium to sulfur desulfurized in the furnace, and when the desulfurization efficiency in the furnace is unchanged, the lower the molar ratio of calcium to sulfur is, the lower the heat loss cost is, and the lower the converted coal consumption is; the dosage of the in-furnace desulfurization limestone is determined by the in-furnace desulfurization calcium-sulfur molar ratio, the in-furnace desulfurization efficiency and the coal quality, and the in-furnace desulfurization efficiency is related to load, the calcium-sulfur molar ratio, bed temperature and air-coal ratio; the amount of limestone desulfurized outside the furnace is related to the desulfurization efficiency outside the furnace.
According to the operation data of the CFB unit, the desulfurization efficiency in the furnace
Figure BDA0003562524650000081
Can be obtained by the following formula:
Figure BDA0003562524650000082
wherein, W c Is the coal feeding amount; s ar The sulfur content is determined by the coal quality; a. the ir The total air volume; k is a radical of f Converting the dimensionless smoke into coefficients;
Figure BDA0003562524650000083
is raw flue gas SO 2 Concentration (under standard conditions).
External desulfurization efficiency
Figure BDA0003562524650000084
Can be obtained by the following formula:
Figure BDA0003562524650000085
wherein the content of the first and second substances,
Figure BDA0003562524650000086
is clean flue gas SO 2 Concentration (under standard conditions).
Therefore, the load, the coal quality, the coal feeding amount, the total air volume, the bed temperature, the rated power of the limestone conveying fan, the rated power of the slurry circulating pump, the feeding flow rate of the limestone in the furnace, the ammonia injection amount, and the raw flue gas SO are selected 2 Concentration, clean flue gas SO 2 Concentration is used as an input variable for the model.
2. Determining the relation between the desulfurization efficiency in the furnace and the load, bed temperature and air-coal ratio
The method specifically comprises the following substeps:
substep S21, selecting bed temperature, coal supply quantity, total air quantity, coal quality and raw flue gas SO under typical load working condition 2 Concentration, in-furnace limestone feed flow rate;
and a substep S22, calculating the in-furnace desulfurization efficiency, the air-coal ratio and the calcium-sulfur molar ratio under the typical load working condition, wherein the calculation of the in-furnace desulfurization calcium-sulfur molar ratio is shown by an equation (3):
Figure BDA0003562524650000091
in the formula (3), the reaction mixture is,
Figure BDA0003562524650000092
in order to achieve the purity of the limestone,
Figure BDA0003562524650000093
a feed flow rate for limestone in the furnace;
and a substep S23, fitting a relational expression of the desulfurization efficiency in the furnace, the molar ratio of calcium to sulfur and the bed temperature under a typical load working condition by adopting a least square method, wherein the relational expression is shown as a formula (4):
Figure BDA0003562524650000094
in the formula (4), A is a dimensionless coefficient, B is a function of bed temperature, m CFB The molar ratio of calcium to sulfur in the furnace is shown;
and a substep S24, correcting the A by using the air-coal ratio under the typical load working condition according to the operation data, and finally obtaining a relational expression of the in-furnace desulfurization efficiency, the calcium-sulfur molar ratio and the bed temperature under the typical load working condition.
3. Synthetic cost model and solution
External desulfurization efficiency
Figure BDA0003562524650000095
With raw flue gas SO 2 The concentration and the total air volume are related, and the external desulfurization efficiency and the raw flue gas SO can be fitted by a least square method 2 The relation between the concentration and the total air volume specifically comprises the following substeps:
substep S31, fitting the external desulfurization efficiency and the raw flue gas SO by a least square method 2 The relationship between the concentration and the total air volume is expressed by the following formula (6):
Figure BDA0003562524650000096
in the formula (6), a1, a2 and a3 are model coefficients;
and a substep S32, assuming that the desulfurization ratio inside the furnace is x under a certain typical load working condition, the desulfurization ratio outside the furnace is 1-x, and combining the formula (1) and the formula (4), determining the molar ratio m of calcium and sulfur desulfurized inside the furnace CFB
Substep S33, determining the limestone consumption cost W in the furnace 1 Expressed as shown in formula (7):
Figure BDA0003562524650000097
in the formula (7), MCaCO3 and MCaO are the relative molecular masses of CaCO3 and CaO respectively, and the unit is g/mol; u1 is the unit price of limestone in the furnace, and the unit is Yuan/kg;
the heat loss cost W 2 Expressed as shown in formula (8):
W 2 =W c [η/(η-Δη)-1]u 2 (8),
in the formula (8), η is the boiler design efficiency; Δ η is boiler heat loss; u2 is the unit price of the fire coal, and the unit is Yuan/kg;
the denitration cost W 3 Expressed as shown in formula (9):
W 3 =k 1 m CFB W c u 3 (9),
in the formula (9), k1 is a cost coefficient; u3 is the urea unit price, unit is yuan/kg;
the power consumption cost W of the conveying fan 4 Expressed as shown in formula (10):
Figure BDA0003562524650000101
in the formula (10), α is a compressed air coefficient; u. of 4 The unit is yuan/kWh for the price of the power on the internet;
the dosage cost W of the limestone outside the furnace 5 Expressed as shown in formula (11):
Figure BDA0003562524650000102
in the formula (11), mCFB and wet are in the molar ratio of calcium to sulfur outside the furnace; u5 is the unit price of limestone outside the furnace, and the unit is Yuan/t.
The power consumption W6 of the circulating pump is expressed as shown in the formula (12):
W 6 =nP w u 4 (12),
in the formula (12), n is the number of the slurry circulating pumps started and is the raw flue gas SO 2 Determining concentration and load; p w The power of a single slurry circulating pump is in kW;
the gypsum income V 7 Expressed as shown in formula (13):
Figure BDA0003562524650000103
wherein eta (H) 2 O) is the water content of the gypsum; u. of 8 The gypsum is monovalent, yuan/kg.
And (3) expressing the comprehensive cost f (x) of the combined desulfurization inside and outside the furnace as shown in a formula (14):
f(x)=W 1 +W 2 +W 3 +W 4 +W 5 +W 6 +W 7 -V 8 0≤x≤x max (14),
in the formula (14), W 1 、W 2 、W 3 、W 4 、W 5 、W 6 、V 7 The desulfurization ratio x in the furnace is determined, and the cost is the consumption cost of limestone in the furnace, the cost of heat loss, the cost of denitration, the cost of power consumption of a conveying fan, the consumption cost of limestone outside the furnace, the power consumption of a circulating pump and the benefit of gypsum; x is the number of max Determined by the in-furnace desulfurization capacity of the CFB unit.
The process of optimizing the optimal in-furnace desulfurization proportion by a genetic algorithm specifically comprises the following substeps:
substep S41, encoding: selecting unsigned binary integers to represent individuals x i
Substep S42, generating an initial population: randomly generating N individuals as an initial population, and setting the iteration number as N;
substep S43, fitness calculation: using the value of the integrated cost function g (x) i) As an individual x i Selecting a fitness function as shown in formula (15):
g(x i )=minf(x i ) (15);
substep S44, selection, crossover, mutation operation: the individuals with higher fitness in the current group are inherited to the next generation; the method adopts a single-point crossing method to carry out crossing operation, adopts a basic bit variation method to carry out variation operation, saves the next generation of group, and increases the iteration times by 1;
substep S45, termination condition judgment: if the iteration times are more than or equal to n, stopping the calculation, and outputting the obtained individual with the minimum fitness as an optimal solution, namely the optimal furnace internal and external desulfurization proportion; otherwise, return to substep S43.
In a specific embodiment, the number of population N is 20, the number of termination iterations N is 80, the crossover probability is 0.4, and the mutation probability is 0.001.
The invention has the following beneficial effects:
(1) and selecting corresponding variables according to the operating characteristics of the CFB unit, and determining the relationship between the in-furnace desulfurization efficiency and the bed temperature and the air-coal ratio by adopting a least square method.
(2) On the basis, after the desulfurization proportion in the furnace is assumed, a desulfurization comprehensive cost model under a typical load working condition is respectively established.
(3) And (3) optimizing the optimal furnace internal and external desulfurization ratio under the typical load working condition by adopting a genetic algorithm and taking the minimum desulfurization comprehensive cost as an objective function.

Claims (6)

1. A method for calculating the optimal furnace inside and outside desulfurization proportion of a circulating fluidized bed unit is characterized by comprising the following steps of:
step 1, establishing a CFB unit furnace internal and external desulfurization comprehensive cost model, and selecting load, coal quality, coal feeding quantity, total air quantity, bed temperature, limestone conveying fan rated power, slurry circulating pump rated power, limestone feeding flow rate in a furnace, ammonia injection quantity and raw flue gas SO 2 Concentration, clean flue gas SO 2 Concentration is used as an input variable of the comprehensive cost model;
step 2, determining relational expressions of desulfurization efficiency in the furnace and load, calcium-sulfur molar ratio, bed temperature and air-coal ratio when establishing the comprehensive cost model;
step 3, determining the SO2 generation concentration by using the input of the comprehensive cost model selected in the step 1, and establishing an in-furnace and out-furnace desulfurization comprehensive cost model under the typical load working condition of the CFB unit after assuming the in-furnace desulfurization proportion;
and 4, solving the optimal in-furnace desulfurization proportion under the typical load working condition by using an intelligent optimization algorithm, and fitting to obtain the optimal in-furnace desulfurization proportion under each load working condition.
2. The method for calculating the optimal furnace internal and external desulfurization proportion of the circulating fluidized bed unit according to claim 1, wherein the process for establishing the CFB unit furnace internal and external desulfurization comprehensive cost model in the step 1 is as follows:
defining the total combined desulfurization cost inside and outside the furnace as comprising: the comprehensive cost is the sum of the total desulfurization cost minus the benefit generated by gypsum; wherein the limestone consumption cost comprises the consumption cost of the in-furnace desulfurized limestone and the consumption cost of the out-furnace desulfurized limestone; the power consumption cost of the equipment comprises the power consumption of a limestone conveying fan in the furnace and the power consumption of a slurry circulating pump outside the furnace; the denitration cost in the furnace is increased along with the increase of the molar ratio of the desulfurized calcium to the desulfurized calcium in the furnace and the denitration cost in the furnace is increased; when the desulfurization efficiency in the furnace is unchanged, the lower the molar ratio of calcium to sulfur is, the lower the heat loss cost is, and the lower the converted coal consumption is; the dosage of the in-furnace desulfurization limestone is determined by the in-furnace desulfurization calcium-sulfur molar ratio, the in-furnace desulfurization efficiency and the coal quality, and the in-furnace desulfurization efficiency is related to load, the calcium-sulfur molar ratio, bed temperature and air-coal ratio; the amount of limestone desulfurized outside the furnace is related to the desulfurization efficiency outside the furnace;
the desulfurization efficiency in the furnace
Figure FDA0003562524640000021
Expressed as shown in formula (1):
Figure FDA0003562524640000022
in the formula (1), the reaction mixture is,W c for coal supply, S ar Is a sulfur content, A ir Total air volume, k f For the dimensionless conversion coefficient of the smoke,
Figure FDA0003562524640000023
is the raw flue gas SO under standard conditions 2 Concentration;
efficiency of external desulfurization of said furnace
Figure FDA0003562524640000024
Expressed as shown in formula (2):
Figure FDA0003562524640000025
in the formula (2), under the standard condition
Figure FDA0003562524640000026
Is clean flue gas SO 2 And (4) concentration.
3. The method for calculating the optimal furnace internal and external desulfurization proportion of the circulating fluidized bed unit according to claim 2, wherein in the step 2, a least square method is adopted to fit a relational expression of the furnace desulfurization efficiency, the load, the calcium-sulfur molar ratio and the bed temperature under the typical load working condition, and the coefficient of the relational expression is corrected by using a wind-coal ratio, and the method specifically comprises the following substeps:
substep S21, selecting bed temperature, coal supply quantity, total air quantity, coal quality and raw flue gas SO under typical load working condition 2 Concentration, in-furnace limestone feed flow rate;
and a substep S22, calculating the in-furnace desulfurization efficiency, the air-coal ratio and the calcium-sulfur molar ratio under the typical load working condition, wherein the calculation of the in-furnace desulfurization calcium-sulfur molar ratio is shown by an equation (3):
Figure FDA0003562524640000031
formula (3)) In the step (1), the first step,
Figure FDA0003562524640000032
in order to obtain the purity of the limestone,
Figure FDA0003562524640000033
a feed flow rate for limestone in the furnace;
and a substep S23 of fitting a relational expression of the desulfurization efficiency in the furnace, the molar ratio of calcium to sulfur and the bed temperature under the typical load working condition by adopting a least square method, wherein the relational expression is expressed as a formula (4):
Figure FDA0003562524640000034
in the formula (4), A is a dimensionless coefficient, B is a function of bed temperature, m CFB The molar ratio of calcium to sulfur in the furnace is shown;
and a substep S24, correcting the A by using the air-coal ratio under the typical load working condition according to the operation data, and finally obtaining a relational expression of the in-furnace desulfurization efficiency, the calcium-sulfur molar ratio and the bed temperature under the typical load working condition.
4. The method according to claim 3, wherein in step 3, assuming the ratio of desulfurization in the furnace, the power consumption of the limestone conveying fan, the power consumption of the slurry circulating pump, and the limestone consumption for desulfurization inside and outside the furnace are determined, wherein the power consumption of the limestone conveying fan is determined by the limestone feeding flow rate in the furnace and the rated power of the limestone conveying fan, and the power consumption of the slurry circulating pump is determined by the number of circulating pumps in operation and the rated power, and the method specifically comprises the following substeps:
substep S31, fitting the external desulfurization efficiency and the raw flue gas SO by a least square method 2 The relationship between the concentration and the total air volume is expressed by the following formula (6):
Figure FDA0003562524640000035
in the formula (6), a1, a2 and a3 are model coefficients;
and a substep S32, assuming that the desulfurization ratio inside the furnace is x under a certain typical load working condition, the desulfurization ratio outside the furnace is 1-x, and combining the formula (1) and the formula (4), determining the molar ratio m of calcium and sulfur desulfurized inside the furnace CFB
Substep S33, determining the limestone consumption cost W in the furnace 1 Expressed as shown in formula (7):
Figure FDA0003562524640000041
in the formula (7), MCaCO3 and MCaO are the relative molecular masses of CaCO3 and CaO respectively, and the unit is g/mol; u1 is the unit price of limestone in the furnace, and the unit is Yuan/kg;
the heat loss cost W 2 Expressed as shown in formula (8):
W 2 =W c [η/(η-Δη)-1]u 2 (8),
in the formula (8), η is the boiler design efficiency; Δ η is boiler heat loss; u2 is the unit price of the fire coal, and the unit is Yuan/kg;
the denitration cost W 3 Expressed as shown in formula (9):
W 3 =k 1 m CFB W c u 3 (9),
in the formula (9), k1 is a cost coefficient; u3 is the urea unit price, unit is yuan/kg;
the power consumption cost W of the conveying fan 4 Expressed as shown in formula (10):
Figure FDA0003562524640000042
in the formula (10), α is a compressed air coefficient; u. of 4 The unit is yuan/kWh for the price of the power on the internet;
the dosage cost W of the limestone outside the furnace 5 Expressed as shown in formula (11):
Figure FDA0003562524640000043
in the formula (11), mCFB and wet are in the molar ratio of calcium to sulfur outside the furnace; u5 is the unit price of limestone outside the furnace, and the unit is Yuan/t;
the power consumption W6 of the circulating pump is expressed as shown in the formula (12):
W 6 =nP w u 4 (12),
in the formula (12), n is the number of the slurry circulating pumps started and is the raw flue gas SO 2 Determining concentration and load; p w The power of a single slurry circulating pump is in kW;
the gypsum income V 7 Expressed as shown in formula (13):
Figure FDA0003562524640000051
wherein eta (H) 2 O) is the water content of the gypsum; u. of 8 Is the unit price of gypsum, and the unit is yuan/kg;
and (3) expressing the comprehensive cost f (x) of the combined desulfurization inside and outside the furnace as shown in a formula (14):
f(x)=W 1 +W 2 +W 3 +W 4 +W 5 +W 6 +W 7 -V 8 0≤x≤x max (14),
in the formula (14), W 1 、W 2 、W 3 、W 4 、W 5 、W 6 、V 7 The desulfurization ratio x in the furnace is determined, and the cost is the consumption cost of limestone in the furnace, the cost of heat loss, the cost of denitration, the cost of power consumption of a conveying fan, the consumption cost of limestone outside the furnace, the power consumption of a circulating pump and the benefit of gypsum; x is the number of max Determined by the in-furnace desulfurization capacity of the CFB unit.
5. The method for calculating the optimal furnace internal and external desulfurization proportion of the circulating fluidized bed unit according to claim 4, wherein in the step 4, the optimal furnace internal desulfurization proportion under the typical load working condition is a solution obtained by genetic algorithm optimization of the minimum integrated cost of the furnace internal and external desulfurization, and the process of genetic algorithm optimization of the minimum integrated cost of the furnace internal and external desulfurization comprises the following sub-steps:
substep S41, encoding: selecting unsigned binary integers to represent an individual x i
Substep S42, generating an initial population: randomly generating N individuals as an initial population, and setting the iteration number as N;
substep S43, fitness calculation: using the value of the integrated cost function g (x) i ) As an individual x i Selecting a fitness function as shown in formula (15):
g(x i )=min f(x i ) (15);
substep S44, selection, crossover, mutation operation: the individuals with higher fitness in the current group are inherited to the next generation; the method adopts a single-point crossing method to carry out crossing operation, adopts a basic bit variation method to carry out variation operation, saves the next generation of group, and increases the iteration times by 1;
substep S45, termination condition determination: if the iteration times are more than or equal to n, stopping the calculation, and outputting the obtained individual with the minimum fitness as an optimal solution, namely the optimal furnace internal and external desulfurization proportion; otherwise, return to substep S43.
6. The method of claim 5, wherein the number of groups N is 20, the number of termination iterations N is 80, the crossover probability is 0.4, and the variation probability is 0.001.
CN202210298473.9A 2022-03-24 2022-03-24 Method for calculating optimal furnace internal and external desulfurization proportion of circulating fluidized bed unit Pending CN114880919A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210298473.9A CN114880919A (en) 2022-03-24 2022-03-24 Method for calculating optimal furnace internal and external desulfurization proportion of circulating fluidized bed unit

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210298473.9A CN114880919A (en) 2022-03-24 2022-03-24 Method for calculating optimal furnace internal and external desulfurization proportion of circulating fluidized bed unit

Publications (1)

Publication Number Publication Date
CN114880919A true CN114880919A (en) 2022-08-09

Family

ID=82667299

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210298473.9A Pending CN114880919A (en) 2022-03-24 2022-03-24 Method for calculating optimal furnace internal and external desulfurization proportion of circulating fluidized bed unit

Country Status (1)

Country Link
CN (1) CN114880919A (en)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100049369A1 (en) * 2008-08-22 2010-02-25 Alstom Technology Ltd. Modeling and control optimization system for integrated fluidized bed combustion process and air pollution control system
CN105066121A (en) * 2015-07-29 2015-11-18 华北电力大学 Dynamic bed temperature prediction system and method of circulating fluidized bed boiler
CN111085091A (en) * 2020-01-02 2020-05-01 神华国能集团有限公司 Method for determining desulfurization proportion inside and outside circulating fluidized bed boiler and desulfurization system
CN111461391A (en) * 2020-01-22 2020-07-28 杭州电子科技大学 Desulfurization collaborative optimization operation method for coal slime fluidized bed boiler
CN113393059A (en) * 2021-07-27 2021-09-14 新奥数能科技有限公司 Method and device for adjusting output distribution of desulfurization system in circulating fluidized bed group
CN114186708A (en) * 2021-10-13 2022-03-15 华北电力大学 Circulating fluidized bed unit SO based on PSO-ELM2Concentration prediction method
CN115169632A (en) * 2022-05-10 2022-10-11 华北电力大学 SO after desulfurization in circulating fluidized bed boiler 2 Emission concentration prediction method

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100049369A1 (en) * 2008-08-22 2010-02-25 Alstom Technology Ltd. Modeling and control optimization system for integrated fluidized bed combustion process and air pollution control system
CN105066121A (en) * 2015-07-29 2015-11-18 华北电力大学 Dynamic bed temperature prediction system and method of circulating fluidized bed boiler
CN111085091A (en) * 2020-01-02 2020-05-01 神华国能集团有限公司 Method for determining desulfurization proportion inside and outside circulating fluidized bed boiler and desulfurization system
CN111461391A (en) * 2020-01-22 2020-07-28 杭州电子科技大学 Desulfurization collaborative optimization operation method for coal slime fluidized bed boiler
CN113393059A (en) * 2021-07-27 2021-09-14 新奥数能科技有限公司 Method and device for adjusting output distribution of desulfurization system in circulating fluidized bed group
CN114186708A (en) * 2021-10-13 2022-03-15 华北电力大学 Circulating fluidized bed unit SO based on PSO-ELM2Concentration prediction method
CN115169632A (en) * 2022-05-10 2022-10-11 华北电力大学 SO after desulfurization in circulating fluidized bed boiler 2 Emission concentration prediction method

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
AMIT KUMAR 等: "Design methodology of a highly loaded tandem rotor and its performance analysis under clean and distorted inflows", MECHANICAL ENGINEERING SCIENCE, 3 April 2021 (2021-04-03), pages 1 - 24 *
WEI ZHANG 等: "A novel operation cost optimization system for mix-burning coal slime circulating fluidized bed boiler unit", APPLIED THERMAL ENGINEERING 148 (2019), 22 November 2018 (2018-11-22), pages 620 *
侯祥松, 李金平, 张海, 赵石铁, 吕俊复, 岳光溪: "石灰石脱硫对循环流化床中NO_x生成和排放的影响", 电站系统工程, no. 01, 25 January 2005 (2005-01-25), pages 7 - 9 *
张洪福: "循环流化机组负荷响应特性与灵活运行控制研究", 知网研学, 1 May 2023 (2023-05-01) *
张磊;王少臣;苑广存;李书元;: "CFB锅炉内外脱硫系统超低排放容量分配优化", 洁净煤技术, no. 04, 15 July 2018 (2018-07-15), pages 124 - 129 *
李锦萍;靳智平;王艳玲;: "一种耦合支持向量机遗传算法的燃烧优化方法研究", 山西电力, no. 01, 28 February 2018 (2018-02-28), pages 46 - 50 *
范常浩;侯鹏飞;: "基于SOPSO算法的CFB机组联合脱硫系统经济性优化研究", 计算机与应用化学, no. 05, 28 October 2019 (2019-10-28), pages 130 - 135 *
谭波;王传志;司硕;刘忠攀;蓝天;: "循环流化床锅炉高效炉内脱硫理论和关键技术", 洁净煤技术, no. 04, 31 December 2020 (2020-12-31), pages 171 - 177 *

Similar Documents

Publication Publication Date Title
CN105116855B (en) A kind of optimal control method of flue gas circulating fluidized bed desulfurization
CN107014217B (en) Coke oven gas utilization and flue gas treatment system of coking plant and treatment method thereof
CN106753489A (en) A kind of pyrolysis of coal steam based on coal-powder boiler, burnt oil and gas co-generation system and technique
CN104696948B (en) Sugar refinery bagasse boiler front end flue gas three-step approach denitrating technique
CN217763522U (en) Natural gas ammonia-doped combustion system for coupling green hydrogen to produce ammonia
CN201714432U (en) Oxygen-enriched combustion coal-fired generating system of solar energy-wind energy integrated power generating set
CN113713599A (en) Circulating fluidized bed semi-dry method low-load desulfurization system and control method
CN114880919A (en) Method for calculating optimal furnace internal and external desulfurization proportion of circulating fluidized bed unit
CN107906553B (en) A kind of Boiler Air Heater's optimization progress control method and device
CN210186850U (en) Synergistic semi-dry desulfurization system
CN112239797A (en) Large blast furnace top combustion type hot blast furnace vault temperature control technology
CN115164610B (en) CO/O-based2Cement decomposing furnace combustion optimization method and system
CN110865623A (en) NO in SCR denitration controlxMeasurement signal substitution system and control method thereof
CN110452717A (en) A kind of method and system for distributing comprehensive treatment coking tail gas rationally
CN113019086B (en) Precise feeding method of desulfurization synergist
CN212283563U (en) Denitration system suitable for biomass boiler
CN209865734U (en) Sintering flue gas circulation is united ozone preoxidation's SOx/NOx control system
CN202849334U (en) System for regulating smoke temperature of water-coal-slurry-fired cracking furnace by recirculation of excess air and smoke gas
CN107957079B (en) The control method of corner tangential firing pulverized-coal fired boiler
CN116379464B (en) Automatic optimizing method for total cost of NOx under full load of coal-fired unit
CN113280322A (en) Process and system for waste electricity hydrogen production combined circulating fluidized bed boiler combustion
CN104373933A (en) Circulating fluidized bed boiler system and combustion air-supply method for low nitrogen oxide discharge
CN117010728B (en) Comprehensive power generation cost optimization method for thermal power enterprises
CN209978023U (en) Flue gas circulating system for circulating fluidized bed boiler
CN220852108U (en) Oil field steam injection system with fossil energy and clean energy coupled

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