CN113378353B - Method for evaluating similarity of NOx flux distribution based on gridding test - Google Patents

Method for evaluating similarity of NOx flux distribution based on gridding test Download PDF

Info

Publication number
CN113378353B
CN113378353B CN202110533814.1A CN202110533814A CN113378353B CN 113378353 B CN113378353 B CN 113378353B CN 202110533814 A CN202110533814 A CN 202110533814A CN 113378353 B CN113378353 B CN 113378353B
Authority
CN
China
Prior art keywords
nox
matrix
flue
flux distribution
flux
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
CN202110533814.1A
Other languages
Chinese (zh)
Other versions
CN113378353A (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.)
Huadian Electric Power Research Institute Co Ltd
Original Assignee
Huadian Electric Power Research Institute 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 Huadian Electric Power Research Institute Co Ltd filed Critical Huadian Electric Power Research Institute Co Ltd
Priority to CN202110533814.1A priority Critical patent/CN113378353B/en
Publication of CN113378353A publication Critical patent/CN113378353A/en
Application granted granted Critical
Publication of CN113378353B publication Critical patent/CN113378353B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]

Abstract

The invention discloses a method for evaluating NOx flux distribution similarity based on a gridding test, which comprises the following steps: step 1: dividing grids on the cross section of the flue, respectively testing the speed and NOx concentration of each grid under different working conditions by using a pitot tube and a flue gas analyzer, and correspondingly calculating NOx flux distribution under different working conditions; step 2: calculating the total flux of NOx of the cross section of the flue and the standard deviation of a NOx flux distribution matrix under different working conditions, calculating a conversion coefficient according to the principle that the total flux is consistent, and solving a point division matrix; and step 3: respectively defining 3 indexes of a relation degree, a concentration degree and a contrast degree according to a mathematical statistics theory; and 4, step 4: and establishing a similarity comprehensive evaluation function by combining the indexes, and giving a similar criterion. According to the invention, by constructing the NOx flux distribution similarity evaluation method, the similarity of the internal flue gas conditions of the flue gas denitration device under different operation conditions is evaluated, and a guiding direction is provided for the multi-condition collaborative optimization and reconstruction of the flue gas conditions.

Description

Method for evaluating similarity of NOx flux distribution based on gridding test
Technical Field
The invention relates to a grid test-based NOx flux distribution similarity evaluation method, and belongs to the technical field of coal-fired unit flue gas treatment.
Background
With the normalization of the ultra-low emission operation management mode of the flue gas of the coal-fired unit, the strict real-time monitoring and control of the emission of the atmospheric pollutants and the improvement of the working enthusiasm of energy conservation and consumption reduction promoted by each enterprise, the seeking of the high-efficiency, economic and stable operation of the flue gas denitration facility becomes the topic of wide attention of each coal-fired power generation enterprise. The realization of the uniform distribution of the ammonia nitrogen molar ratio at the inlet of the denitration catalyst is the key for solving the problem of cooperative optimization of reaching the emission standard and reducing the consumption of reducing agent. Usually, a mixing device is additionally arranged behind an ammonia injection device and in front of a denitration catalyst to generate vortex to enhance the mixing of the denitration reducing agent and NOx in flue gas, and the matching relation between the denitration reducing agent amount injected into a flue by each partition and the NOx content in the flue gas is adjusted by adjusting the opening degree of each partition valve of the ammonia injection device. Tests show that the upstream flue gas flow rate of each partition of the ammonia injection device and the concentration of NOx in the flue gas can be adjusted to change along with the change of coal quality and combustion. At present, the flow field distribution characteristics of all working conditions are analyzed by adopting the flue gas flow velocity, the NOx concentration or the relative standard deviation of the ammonia nitrogen molar ratio in an isolated manner, the time-varying relation of the NOx flux distribution in the flue gas of all subareas of the ammonia spraying device is ignored, namely the relation degree of the composite index NOx flux distribution characteristics of the flue gas flow velocity, the NOx concentration and the subarea area among different working conditions, and the possibility of realizing the uniform distribution of the ammonia nitrogen molar ratio under all working conditions is reflected more visually and more accurately. It is therefore a meaningful task to study and evaluate the similarity of NOx flux distribution upstream of the ammonia injection grid under various conditions.
Disclosure of Invention
The invention aims to overcome the defects in the prior art, and provides a method for evaluating NOx flux distribution similarity based on a gridding test, which is used for comprehensively analyzing the similar conditions of NOx flux distribution under multiple working conditions.
The technical scheme adopted by the invention for solving the problems is as follows: a NOx flux distribution similarity evaluation method based on a gridding test is characterized by comprising the following steps:
step 1: dividing grids on the cross section of the flue, respectively testing the speed and NOx concentration of each grid under different working conditions by using a pitot tube and a flue gas analyzer, and correspondingly calculating NOx flux distribution under different working conditions;
step 2: calculating the total flux of NOx of the cross section of the flue and the standard deviation of a NOx flux distribution matrix under different working conditions, calculating a conversion coefficient according to the principle that the total flux is consistent, and solving a point division matrix;
and step 3: respectively defining 3 indexes of a relation degree, a concentration degree and a contrast degree according to a mathematical statistics theory;
and 4, step 4: and establishing a similarity comprehensive evaluation function by combining the indexes, and giving a similar criterion.
In the step 1, the cross section of the flue is divided into grids with m multiplied by n (length multiplied by width) areas a, the geometric center of each grid is taken as a measurement point, and the speed V of each grid is testedijAnd NOx concentration Cij(i-1 … m, j-1 … n); the distribution expression matrix of the NOx flux of the flue section of the working condition 1 is as follows:
Figure BDA0003067462330000021
the distribution expression matrix of NOx flux of the flue section of the working condition 2 is as follows:
Figure BDA0003067462330000022
in step 2, the matrix Φ is expressed according to the NOx flux distribution1Calculating the total flux phi' and standard deviation sigma of NOx in the section of the flue under the working condition 11
Figure BDA0003067462330000023
Figure BDA0003067462330000024
Representation matrix phi from NOx flux distribution2Calculating the total NOx flux phi' and standard deviation sigma of the section of the flue under the working condition 22
Figure BDA0003067462330000025
Figure BDA0003067462330000026
Calculating a conversion coefficient alpha and a point division matrix psi according to the total NOx flux under different working conditions:
α=Φ′Φ″
Figure BDA0003067462330000031
in the step 3: satisfaction in the elements of the point division matrix Ψ
Figure BDA0003067462330000032
Counting the number of the first and second groups, namely the same number K;
satisfaction in the elements of the point division matrix Ψ
Figure BDA0003067462330000033
Counting the number of the adjacent numbers P;
satisfaction in the elements of the point division matrix Ψ
Figure BDA0003067462330000034
Counting the number of the pairs, namely counting the pairs Q;
obviously, the sum of the same number, the adjacent number and the opposite number is the number of the dot division matrix elements, i.e., K + P + Q ═ mn;
introducing degree of relation which is a linear function of identity number, adjacent number and opposite number, i.e.
Figure BDA0003067462330000035
Introducing a concentration degree which is an exponential function on the point divided matrix elements, i.e.
Figure BDA0003067462330000036
Contrast is introduced as an exponential function representing the standard deviation of the matrix with respect to NOx flux distribution
Figure BDA0003067462330000037
In step 4, establishing a similar comprehensive evaluation function f ═ r + s + t)/3;
when the value range of the comprehensive evaluation function is within [0.7,1], the NOx flux distribution under two working conditions is determined to be similar, and the design requirement that the flue gas condition of the flue section is basically unchanged is met.
Compared with the prior art, the invention has the following advantages and effects: according to the invention, by constructing the NOx flux distribution similarity evaluation method, the similarity of the internal flue gas conditions of the flue gas denitration device under different operation conditions is evaluated, and a guiding direction is provided for the multi-condition collaborative optimization and reconstruction of the flue gas conditions.
Drawings
FIG. 1 is a schematic diagram of a cross-sectional meshing of a flue in an embodiment of the present invention.
Detailed Description
The present invention will be described in further detail below by way of examples with reference to the accompanying drawings, which are illustrative of the present invention and are not to be construed as limiting the present invention.
Examples are given.
Referring to fig. 1, in this embodiment, a method for evaluating NOx flux distribution similarity based on a gridding test includes the following specific steps:
step 1: carrying out grid division on the cross section of the flue, respectively testing the speed and the NOx concentration of the central point of each grid under different working conditions by using a pitot tube and a flue gas analyzer, and correspondingly calculating the NOx flux distribution under different working conditions;
and 2, step: calculating the total flux of NOx of the cross section of the flue and the standard deviation of a NOx flux distribution matrix under different working conditions, calculating a conversion coefficient according to the principle that the total flux is consistent, and solving a point division matrix;
and step 3: respectively defining 3 indexes of a relation degree, a concentration degree and a contrast degree according to a mathematical statistics theory;
and 4, step 4: and establishing a similar comprehensive evaluation function by combining the indexes, and giving a similar criterion.
The invention is further illustrated with reference to fig. 1 as follows:
under the working condition 1, the actual measurement results of the flue section speed distribution, the NOx concentration distribution, the NOx flux distribution, the flux distribution matrix standard deviation and the total flux are respectively as follows:
Figure BDA0003067462330000041
Figure BDA0003067462330000042
Figure BDA0003067462330000043
σ1=0.58
Φ′=179
under the working condition 2, the actual measurement results of the flue section speed distribution, the NOx concentration distribution, the NOx flux distribution, the flux distribution matrix standard deviation and the total flux are respectively as follows:
Figure BDA0003067462330000051
Figure BDA0003067462330000052
Figure BDA0003067462330000053
σ2=0.51
Φ″=120
next, the conversion coefficients, the dot division matrix, and the respective analysis indexes are calculated.
Conversion coefficient:
α=179/120=1.49
dot division matrix:
Figure BDA0003067462330000054
the number of elements in the dot division matrix is 30;
satisfies | Φ in the elements of the dot division matrixij-14 identical numbers K | < 0.1;
the element of the difference matrix satisfies | phi more than or equal to 0.1ij-12, the number of neighbors P ═ 1| < 0.3;
satisfies | Φ in elements of the difference matrixij-4 as many as 1 ≧ 0.3 or the number of opponents Q;
and (3) calculating the contact degree:
Figure BDA0003067462330000055
and (3) calculating the concentration ratio:
Figure BDA0003067462330000061
calculating the contrast ratio:
Figure BDA0003067462330000062
similar comprehensive evaluation function:
f=(r+s+t)/3=0.75
when the value range of the comprehensive evaluation function is within [0.7,1], the NOx flux distributions under the two working conditions are considered to have high similarity. From the similar comprehensive evaluation function calculated value of this embodiment, the flue gas distribution condition under the different operating modes is similar, and is comparatively stable to operating mode change flue gas distribution promptly, and the fluctuation is less.
Those not described in detail in this specification are well within the skill of the art.
Although the present invention has been described with reference to the above embodiments, it should be understood that the scope of the present invention is not limited thereto, and that various changes and modifications can be made by those skilled in the art without departing from the spirit and scope of the present invention.

Claims (1)

1. A NOx flux distribution similarity evaluation method based on a gridding test is characterized by comprising the following steps:
step 1: dividing grids on the cross section of the flue, respectively testing the speed and NOx concentration of each grid under different working conditions by using a pitot tube and a flue gas analyzer, and correspondingly calculating NOx flux distribution under different working conditions;
step 2: calculating the total flux of NOx of the cross section of the flue and the standard deviation of a NOx flux distribution matrix under different working conditions, calculating a conversion coefficient according to the principle that the total flux is consistent, and solving the matrix;
and step 3: respectively defining 3 indexes of a link ratio, a concentration ratio and a contrast ratio according to a mathematical statistics theory;
and 4, step 4: establishing a similarity comprehensive evaluation function by combining the indexes, and giving a similar criterion;
in the step 1, the cross section of the flue is divided into grids with m multiplied by n areas of a, the geometric center of each grid is taken as a measuring point, and the speed V of each grid is testedijAnd NOx concentration Cij(i-1 … m, j-1 … n); the distribution expression matrix of NOx flux of the flue section of the working condition 1 is as follows:
Figure FDA0003561133650000011
under the working condition 2, the distribution expression matrix of the NOx flux of the flue section is as follows:
Figure FDA0003561133650000012
in step 2, the matrix Φ is expressed according to the NOx flux distribution1Calculating the total flux phi' and standard deviation sigma of NOx in the section of the flue under the working condition 11
Figure FDA0003561133650000013
Figure FDA0003561133650000014
Representation matrix phi from NOx flux distribution2Calculating the total NOx flux phi' and standard deviation sigma of the section of the flue under the working condition 22
Figure FDA0003561133650000015
Figure FDA0003561133650000021
Calculating a conversion coefficient alpha and a matrix psi according to the total NOx flux under different working conditions:
α=Φ′/Φ″
Figure FDA0003561133650000022
in the step 3:
for satisfaction of the elements of the matrix Ψ
Figure FDA0003561133650000023
Counting the number of the first and second groups, namely the same number K;
for satisfaction of the elements of the matrix Ψ
Figure FDA0003561133650000024
Counting the number of the adjacent numbers P;
for satisfaction of the elements of the matrix Ψ
Figure FDA0003561133650000025
Counting the number of the pairs, namely counting the pairs Q;
the sum of the same number, the adjacent number and the opposite number is the number of matrix elements, namely K + P + Q ═ mn;
the degree of relation is introduced, which is a linear function of the same number, the adjacent number and the opposite number,namely, it is
Figure FDA0003561133650000026
Introducing a concentration which is an exponential function with respect to the matrix elements, i.e.
Figure FDA0003561133650000027
Contrast is introduced as an exponential function representing the standard deviation of the matrix with respect to NOx flux distribution
Figure FDA0003561133650000028
In step 4, establishing a similar comprehensive evaluation function f ═ r + s + t)/3;
when the value range of the comprehensive evaluation function is within [0.7,1], the NOx flux distribution under two working conditions is determined to be similar, and the design requirement that the flue gas condition of the flue section is basically unchanged is met.
CN202110533814.1A 2021-05-14 2021-05-14 Method for evaluating similarity of NOx flux distribution based on gridding test Active CN113378353B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110533814.1A CN113378353B (en) 2021-05-14 2021-05-14 Method for evaluating similarity of NOx flux distribution based on gridding test

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110533814.1A CN113378353B (en) 2021-05-14 2021-05-14 Method for evaluating similarity of NOx flux distribution based on gridding test

Publications (2)

Publication Number Publication Date
CN113378353A CN113378353A (en) 2021-09-10
CN113378353B true CN113378353B (en) 2022-05-27

Family

ID=77571051

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110533814.1A Active CN113378353B (en) 2021-05-14 2021-05-14 Method for evaluating similarity of NOx flux distribution based on gridding test

Country Status (1)

Country Link
CN (1) CN113378353B (en)

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5413477A (en) * 1992-10-16 1995-05-09 Gas Research Institute Staged air, low NOX burner with internal recuperative flue gas recirculation
GB2481433A (en) * 2010-06-24 2011-12-28 Gm Global Tech Operations Inc Determining NOx concentration upstream of an SCR catalyst
CN109493250B (en) * 2018-11-06 2021-09-07 大唐南京环保科技有限责任公司 Method for evaluating denitration capability of SCR reactor
CN111359438A (en) * 2020-04-09 2020-07-03 国电环境保护研究院有限公司 Ammonia spraying optimization control method and control system for SCR flue gas denitration device
CN112221323A (en) * 2020-08-11 2021-01-15 华电电力科学研究院有限公司 NOx partition measuring method based on flow velocity weighting
CN112446000B (en) * 2020-10-26 2023-11-28 华电电力科学研究院有限公司 Flow field similarity evaluation method based on speed distribution

Also Published As

Publication number Publication date
CN113378353A (en) 2021-09-10

Similar Documents

Publication Publication Date Title
CN103499671B (en) Method for testing activity of denitration catalyst of selective catalytic reduction (SCR) system of thermal power plant on site
CN111523246B (en) Composite flue gas flow field measurement method and measurement system
CN106248864A (en) A kind of SCR denitration life-span prediction method based on magnanimity service data
CN110094251B (en) SCR catalyst performance degradation analysis method based on time-interval multi-model modeling
CN103473408A (en) Method for restoring missing air temperature records on basis of spatial-temporal information fusion
CN203216921U (en) Smoke distribution measuring device in SCR (selective catalytic reduction) smoke desulfurizing system
CN114036736B (en) Causal network learning method based on local Granges causal analysis
CN113689917A (en) Visual ammonia injection optimization method and device based on SCR outlet NOx concentration
CN109740301A (en) A kind of accounting method of the Gas Generator Set carbon emission amount based on BP neural network
CN104504238A (en) Method for predicting nitrogen oxide emission concentration of coal-fired unit based on neural network
CN110580936A (en) Method and system for predicting service life of medium-low temperature SCR denitration catalyst
CN113378353B (en) Method for evaluating similarity of NOx flux distribution based on gridding test
CN106971076A (en) A kind of water quality of river Monitoring Data sequential encryption method
CN105975789A (en) Ammonia-escaping-rate online obtaining method for desulfurization and denitrification control
CN116664360A (en) Carbon emission accounting and carbon flow tracking method for key industrial enterprises
CN112446000B (en) Flow field similarity evaluation method based on speed distribution
CN108593850A (en) The representative measurement point characteristic parameter that flue gas composition concentration total cross-section measures determines method
CN109260948B (en) Denitration uniformity measurement optimization system and method under fluctuating load
CN112379035B (en) System and method for evaluating performance of SCR (selective catalytic reduction) denitration device catalyst on line
CN112191101B (en) Based on NO x Ammonia spraying grid leveling method for quickly measuring ammonia escape concentration
CN109766577B (en) SCR (selective catalytic reduction) reactor information rapid acquisition method based on PLS (partial least squares) local fitting
Elmi Effect of Compensation, Job Satisfaction & Work Stress on Turnover Intention
CN212134259U (en) Online multipoint monitoring sampling device for denitration flue gas parameters of thermal power plant
CN212060151U (en) Multipoint itinerant switching sampling analysis system based on kriging interpolation method
CN216955914U (en) A gas turbine outlet flue NOx concentration rapid survey system for deNOx systems

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