CN114758736A - Method for calculating activation energy of catalyst of scr device - Google Patents
Method for calculating activation energy of catalyst of scr device Download PDFInfo
- Publication number
- CN114758736A CN114758736A CN202210492416.4A CN202210492416A CN114758736A CN 114758736 A CN114758736 A CN 114758736A CN 202210492416 A CN202210492416 A CN 202210492416A CN 114758736 A CN114758736 A CN 114758736A
- Authority
- CN
- China
- Prior art keywords
- aging
- catalyst
- aging process
- activation energy
- calculating
- 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
Links
- 239000003054 catalyst Substances 0.000 title claims abstract description 50
- 230000004913 activation Effects 0.000 title claims abstract description 24
- 238000000034 method Methods 0.000 title claims abstract description 22
- 230000032683 aging Effects 0.000 claims abstract description 73
- 230000000694 effects Effects 0.000 claims abstract description 11
- 238000006243 chemical reaction Methods 0.000 claims abstract description 7
- 230000006866 deterioration Effects 0.000 claims abstract description 5
- 238000009825 accumulation Methods 0.000 claims description 6
- SYHGEUNFJIGTRX-UHFFFAOYSA-N methylenedioxypyrovalerone Chemical compound C=1C=C2OCOC2=CC=1C(=O)C(CCC)N1CCCC1 SYHGEUNFJIGTRX-UHFFFAOYSA-N 0.000 claims description 3
- 230000015556 catabolic process Effects 0.000 abstract description 4
- 238000006731 degradation reaction Methods 0.000 abstract description 4
- 238000000746 purification Methods 0.000 abstract description 4
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 abstract description 2
- 239000003546 flue gas Substances 0.000 abstract description 2
- 238000012512 characterization method Methods 0.000 abstract 1
- 239000007789 gas Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 229910044991 metal oxide Inorganic materials 0.000 description 1
- 150000004706 metal oxides Chemical class 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002808 molecular sieve Substances 0.000 description 1
- 229910000510 noble metal Inorganic materials 0.000 description 1
- URGAHOPLAPQHLN-UHFFFAOYSA-N sodium aluminosilicate Chemical compound [Na+].[Al+3].[O-][Si]([O-])=O.[O-][Si]([O-])=O URGAHOPLAPQHLN-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/30—Prediction of properties of chemical compounds, compositions or mixtures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Crystallography & Structural Chemistry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Operations Research (AREA)
- Chemical & Material Sciences (AREA)
- Algebra (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computing Systems (AREA)
- Investigating Or Analyzing Non-Biological Materials By The Use Of Chemical Means (AREA)
- Catalysts (AREA)
Abstract
The invention discloses a method for calculating activation energy of a scr device catalyst, which relates to the technical field of flue gas purification, and the technical scheme is as follows: s1: carrying out single-point temperature steady-state aging on the same SCR catalyst or two catalysts with the same batch and the same model at different temperatures, and recording the single-point temperature steady-state aging as an aging process A and an aging process B; the temperatures of the aging process A and the aging process B are respectively TAAnd TBAging times are respectively tAAnd tCThe deterioration rates of the conversion efficiency before and after the aging of the catalyst are X% and Y% respectively; s2: calculating T based on the theoretical basis of linear agingBAging time t at temperatureB(ii) a S3: respectively calculating equivalent thermal effect damage of the aging process A and the aging process B: s4: and (5) performing characterization by an Arrhenius equation to finally obtain an expression of the activation energy. The catalyst degradation rate before and after the catalyst aging can be calculated by an equivalent back-stepping method by acquiring the temperature and the time of two aging processesThe activation energy has important significance for the calibration, aging judgment and life estimation of the SCR.
Description
Technical Field
The invention relates to the technical field of flue gas purification, in particular to a method for calculating the activation energy of a catalyst of an scr device.
Background
With the upgrading of emission regulations, the requirements of the regulations cannot be met only by the in-cylinder purification technology. At present, NH3SCR catalyst is currently the most mature technology for reducing NOx emissions from diesel vehicles, NH3The catalyst of the SCR catalyst undergoes mainly 3 material changes: noble metal catalysts, metal oxide catalysts and molecular sieve catalysts, the performance of the catalyst directly determines the purification effect of NOx emission and the durability of the catalyst, and the performance of the catalyst is the activation energy. The activation energy of each catalyst is different due to the complex and diverse catalyst type formulas, and the conventional determination of the activation energy of the catalyst requires special machine equipment and complicated steps, which increases the difficulty for the calibration, aging determination and service life estimation of the SCR. There is a need for a convenient means to obtain the catalyst activation energy during the aging of the catalyst.
Disclosure of Invention
The invention aims to provide a method for calculating the activation energy of a catalyst of an SCR device, which can calculate the activation energy of the catalyst by an equivalent backstepping method by acquiring the temperature and time of two aging processes and the degradation rate of the catalyst before and after aging, and has important significance for the calibration, aging judgment and service life estimation of SCR.
The technical purpose of the invention is realized by the following technical scheme: a method for calculating the activation energy of a catalyst of an scr device specifically comprises the following steps:
s1: carrying out single-point temperature steady-state aging on the same SCR catalyst or two catalysts with the same batch and the same model at different temperatures, and recording the single-point temperature steady-state aging as an aging process A and an aging process B; the temperatures of the aging process A and the aging process B are respectively TAAnd TBAging times are respectively tAAnd tCThe deterioration rates of the conversion efficiency before and after the catalyst aging are X% and Y%, respectively;
s2: calculating T based on the theoretical basis of linear agingBDeterioration of catalyst conversion efficiency at temperatureAging time t required for the rate X%B：
S3: respectively calculating equivalent thermal effect damage of the aging process A and the aging process B:
wherein k isAThe rate of accumulation of thermal effects, k, of Process ABThe rate of accumulation of thermal effects for process B;
s4: characterized by the Arrhenius equation, there are:
wherein A isAIs a pre-exponential factor of aging process A, ABIs a pre-exponential factor of the aging process B; r is ideal gas constant, R is 8.314J mol-1·K-1，EaJ.mol for activation energy-1；
S5: according to the aging process A and the aging process B, T is respectivelyAAnd TBSingle point temperature steady state aging at temperature, the formula in step S4 is transformed as:
s6: from the formula of step S5:
s7: and f, taking logarithm with e as a base on two sides of the equal sign of the formula of the step S6 to obtain:
s8: from the formula of step S7:
s9: the temperature difference between the catalyst according to aging A and aging B is less than 50 deg.C, so that A is smallA＝ABTo obtain:
s10: will be known as AA、AB、TA、TBSubstituting R into the formula in step S9 to obtain activation energy Ea。
In conclusion, the invention has the following beneficial effects: by combining the research theory that the aging process of the SCR catalyst is basically linear, the activation energy of the catalyst can be accurately calculated by an equivalent backstepping method by acquiring the temperature and the time of the two aging processes and the degradation rate of the catalyst before and after aging, and the method has important significance for the calibration, aging judgment and service life estimation of the SCR.
Drawings
Fig. 1 is a flow chart of a method for calculating the activation energy of a catalyst of an scr device in an embodiment of the invention.
Detailed Description
The present invention is described in further detail below with reference to FIG. 1.
Example (b): a method for calculating the activation energy of a catalyst of a scr device is shown in figure 1, and specifically comprises the following steps:
s1: carrying out single-point temperature steady-state aging on the same SCR catalyst or two catalysts of the same batch and the same type at different temperatures, and recording the single-point temperature steady-state aging as an aging process A and an aging process B; the temperatures of the aging process A and the aging process B are respectively TAAnd TBAging times are respectively tAAnd tCThe deterioration rates of the conversion efficiency before and after the aging of the catalyst are X% and Y% respectively;
s2: calculating T based on the theoretical basis of linear agingBAging time t at temperatureB：
S3: respectively calculating equivalent thermal effect damage of the aging process A and the aging process B:
wherein k isAThe rate of accumulation of thermal effects, k, of Process ABThe rate of accumulation of thermal effects for process B;
s4: characterized by the Arrhenius equation, there are:
wherein A isAIs a pre-exponential factor of aging process A, ABIs a pre-exponential factor of the aging process B; r is ideal gas constant, R is 8.314J mol-1·K-1，EaJ.mol for activation energy-1；
S5: according to the aging process A and the aging process B, T is respectivelyAAnd TBSingle point temperature steady state aging at temperature, the formula in step S4 is transformed as:
s6: from the formula of step S5:
s7: taking e as a base for logarithm on two sides of the equation equal sign of the step S6 to obtain:
s8: from the formula of step S7:
s9: the temperature difference between the catalyst according to aging A and aging B is less than 50 deg.C, so that A is smallA＝ABTo obtain:
s10: will be known as AA、AB、TA、TBSubstituting R into the formula in step S9 to obtain activation energy Ea。
In this embodiment, the approximate ratios of X% and Y% are not equal, and in order to perform the equivalent back-stepping, the time conversion of step S2 needs to be performed; arrhenius equation in S4 is
Combining with the equation in step S3, S4 can be obtained.
The working principle is as follows: by combining the research theory that the aging process of the SCR catalyst is basically linear, the activation energy of the catalyst can be accurately calculated by an equivalent backstepping method by acquiring the temperature and the time of the two aging processes and the degradation rate of the catalyst before and after aging, and the method has important significance for the calibration, aging judgment and service life estimation of the SCR.
The present embodiment is only for explaining the present invention, and it is not limited to the present invention, and those skilled in the art can make modifications without inventive contribution to the present embodiment as required after reading the present specification, but all of them are protected by patent law within the scope of the present invention.
Claims (1)
1. A method for calculating the activation energy of a catalyst of a scr device is characterized by comprising the following steps: the method specifically comprises the following steps:
s1: carrying out single-point temperature steady-state aging on the same SCR catalyst or two catalysts with the same batch and the same model at different temperatures, and recording the single-point temperature steady-state aging as an aging process A and an aging process B; the temperatures of the aging process A and the aging process B are respectively TAAnd TBAging times are respectively tAAnd tCThe deterioration rates of the conversion efficiency before and after the catalyst aging are X% and Y%, respectively;
s2: calculating T based on the theoretical basis of linear agingBAging time t required for deterioration rate of conversion efficiency of catalyst at temperature to be X%B：
S3: respectively calculating equivalent thermal effect damage of the aging process A and the aging process B:
wherein k isAThe rate of accumulation of thermal effects, k, of Process ABThe rate of accumulation of thermal effects for process B;
s4: characterized by the Arrhenius equation, there are:
wherein A isAIs a pre-exponential factor of aging process A, ABIs a pre-exponential factor of the aging process B; r is ideal gas constant, R is 8.314J mol-1·K-1，EaJ.mol for activation energy-1；
S5: according to the aging process A and the aging process B, T is respectivelyAAnd TBSingle point temperature steady state aging at temperature, the formula in step S4 is transformed as:
s6: from the formula of step S5:
s7: and f, taking logarithm with e as a base on two sides of the equal sign of the formula of the step S6 to obtain:
s.: from the formula of step S7:
s9: the temperature difference between the catalyst according to aging A and aging B is less than 50 deg.C, so that A is smallA＝ABTo obtain:
s10: will be known as AA、AB、TA、TBSubstituting R into the formula in step S9 to obtain activation energy Ea。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210492416.4A CN114758736A (en) | 2022-05-07 | 2022-05-07 | Method for calculating activation energy of catalyst of scr device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210492416.4A CN114758736A (en) | 2022-05-07 | 2022-05-07 | Method for calculating activation energy of catalyst of scr device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114758736A true CN114758736A (en) | 2022-07-15 |
Family
ID=82335063
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210492416.4A Pending CN114758736A (en) | 2022-05-07 | 2022-05-07 | Method for calculating activation energy of catalyst of scr device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114758736A (en) |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130008149A1 (en) * | 2011-07-06 | 2013-01-10 | Ford Global Technologies, Llc | Estimation of efficiency for aged catalysts |
| CN104234802A (en) * | 2014-07-14 | 2014-12-24 | 浙江大学 | SCR (Selective Catalytic Reduction) catalyst aging judgment method based on NOx feedback and ammonia storage prediction |
| EP2966275A1 (en) * | 2014-07-08 | 2016-01-13 | Peugeot Citroën Automobiles Sa | Exhaust gas purification device of a motor vehicle engine, comprising a selective catalytic reduction catalyst |
| CN106055840A (en) * | 2016-06-30 | 2016-10-26 | 重庆长安汽车股份有限公司 | Three-way catalytic converter lectotype design method |
| CN106709255A (en) * | 2016-12-29 | 2017-05-24 | 无锡威孚力达催化净化器有限责任公司 | SCR catalyst aging computing method |
| US20180232499A1 (en) * | 2017-02-10 | 2018-08-16 | Phillips 66 Company | Configuring multiple catalytic beds |
| CN111157049A (en) * | 2020-01-13 | 2020-05-15 | 北京理工大学 | System and method for rapidly evaluating aging performance of SCR (selective catalytic reduction) postprocessor of diesel vehicle |
| CN111999610A (en) * | 2020-08-11 | 2020-11-27 | 国网天津市电力公司电力科学研究院 | Dry-type insulation equipment aging evaluation and service life prediction method based on activation energy |
| WO2020259790A1 (en) * | 2019-06-24 | 2020-12-30 | Toyota Motor Europe | Computer-implemented method for simulating chemical reactions in an aged catalytic converter |
| US20210095590A1 (en) * | 2019-05-09 | 2021-04-01 | Robert Bosch Gmbh | Method of monitoring an scr catalyst |
| CN113670389A (en) * | 2021-09-28 | 2021-11-19 | 潍柴动力股份有限公司 | Three-way catalyst aging test method, device, equipment, storage medium and program |
| CN113969819A (en) * | 2021-10-21 | 2022-01-25 | 奇瑞商用车（安徽）有限公司 | Preparation method of critical part of catalytic converter based on rack aging |
| CN114216921A (en) * | 2021-05-08 | 2022-03-22 | 中国科学院上海微系统与信息技术研究所 | Method for testing activation energy of catalyst |
-
2022
- 2022-05-07 CN CN202210492416.4A patent/CN114758736A/en active Pending
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130008149A1 (en) * | 2011-07-06 | 2013-01-10 | Ford Global Technologies, Llc | Estimation of efficiency for aged catalysts |
| EP2966275A1 (en) * | 2014-07-08 | 2016-01-13 | Peugeot Citroën Automobiles Sa | Exhaust gas purification device of a motor vehicle engine, comprising a selective catalytic reduction catalyst |
| CN104234802A (en) * | 2014-07-14 | 2014-12-24 | 浙江大学 | SCR (Selective Catalytic Reduction) catalyst aging judgment method based on NOx feedback and ammonia storage prediction |
| CN106055840A (en) * | 2016-06-30 | 2016-10-26 | 重庆长安汽车股份有限公司 | Three-way catalytic converter lectotype design method |
| CN106709255A (en) * | 2016-12-29 | 2017-05-24 | 无锡威孚力达催化净化器有限责任公司 | SCR catalyst aging computing method |
| US20180232499A1 (en) * | 2017-02-10 | 2018-08-16 | Phillips 66 Company | Configuring multiple catalytic beds |
| US20210095590A1 (en) * | 2019-05-09 | 2021-04-01 | Robert Bosch Gmbh | Method of monitoring an scr catalyst |
| WO2020259790A1 (en) * | 2019-06-24 | 2020-12-30 | Toyota Motor Europe | Computer-implemented method for simulating chemical reactions in an aged catalytic converter |
| CN111157049A (en) * | 2020-01-13 | 2020-05-15 | 北京理工大学 | System and method for rapidly evaluating aging performance of SCR (selective catalytic reduction) postprocessor of diesel vehicle |
| CN111999610A (en) * | 2020-08-11 | 2020-11-27 | 国网天津市电力公司电力科学研究院 | Dry-type insulation equipment aging evaluation and service life prediction method based on activation energy |
| CN114216921A (en) * | 2021-05-08 | 2022-03-22 | 中国科学院上海微系统与信息技术研究所 | Method for testing activation energy of catalyst |
| CN113670389A (en) * | 2021-09-28 | 2021-11-19 | 潍柴动力股份有限公司 | Three-way catalyst aging test method, device, equipment, storage medium and program |
| CN113969819A (en) * | 2021-10-21 | 2022-01-25 | 奇瑞商用车（安徽）有限公司 | Preparation method of critical part of catalytic converter based on rack aging |
Non-Patent Citations (2)
| Title |
|---|
| 杜永强;郑坚;彭威;张晓;顾志旭;: "基于修正Arrhenius活化能方法的HTPB推进剂贮存寿命预估", 固体火箭技术, no. 01, 28 February 2017 (2017-02-28) * |
| 王志鹏;李鹏飞;陈媛;禚玉群;: "SCR脱硝反应活化能和指前因子分析", 环境工程, no. 10, 22 October 2018 (2018-10-22) * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103670633B (en) | Particulate filter based on pressure reduction is diagnosed | |
| EP3775513B1 (en) | Method for controlling urea mixers to reduce nox emission from combustion engines | |
| CN112761757B (en) | DPF initialization self-learning method and device | |
| CN104948280B (en) | Reducing agent quality system including rationality diagnostic | |
| CN110486130B (en) | Diesel engine DPF carbon loading capacity initialization self-learning method | |
| EP2246294B1 (en) | Control method for fuel reformer | |
| JP5884920B2 (en) | Exhaust gas purification system for internal combustion engine | |
| RU2013142080A (en) | METHOD BASED ON THE MODEL OF DETERMINING THE DISTRIBUTION OF THE TEMPERATURE OF THE EMISSION GAS UNIT | |
| JP2012154216A (en) | Detection apparatus | |
| CN104234799A (en) | Method and system for adapting the clean filter correction map for a selective catalyst reduction filter | |
| CN110145387A (en) | Stationary gas generator tail gas processing method and its system | |
| CN114758736A (en) | Method for calculating activation energy of catalyst of scr device | |
| CN110872979B (en) | Method for predicting performance of emission control system | |
| CN107597104A (en) | Oxide carried Pt catalyst of core shell structure and its preparation method and application | |
| CN107983345A (en) | Exhaust emissions of diesel engine DOC oxidation catalyst preparation methods | |
| CN113969819A (en) | Preparation method of critical part of catalytic converter based on rack aging | |
| CN112324547B (en) | Selective catalytic reduction control method and system suitable for mixed fuel | |
| JP2007502388A5 (en) | ||
| CN102575556B (en) | Exhaust catalyst temperature estimating apparatus and exhaust catalyst temperature estimation method for internal combustion engine | |
| JP7114887B2 (en) | Exhaust purification device for internal combustion engine | |
| CN113914982A (en) | System and method for detecting passive regeneration efficiency of particle trap | |
| JP2011069296A (en) | Deterioration predicting method of denitration catalyst, deterioration countermeasures method, deterioration countermeasures system, and design method of exhaust gas treatment system | |
| CN113202605A (en) | Method for calculating ammonia leakage amount of SCR (Selective catalytic reduction) aftertreatment system | |
| JP2012168062A (en) | Simulation device | |
| CN111749769A (en) | DPF active regeneration efficiency detection system and method |
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 |