Method, device and equipment for evaluating service life of denitration catalyst

Abstract

The embodiment of the invention provides a denitration catalyst service life evaluation method, device and equipment. The method comprises the following steps: acquiring current operation parameters of the denitration system; according to the current operation parameters of the denitration system, the flue gas flow entering the denitration system is corrected, and the corrected flue gas flow is obtained; calculating the surface velocity of a denitration catalyst of the denitration system according to the corrected flue gas flow; calculating the dynamic decay time of the denitration catalyst in unit time according to the surface velocity; and evaluating the service life of the denitration catalyst according to the dynamic decay time in the unit time. In this way, the service life of the denitration catalyst can be efficiently, rapidly and accurately evaluated through the current operation parameters of the denitration system, and it is avoided that the service life of the denitration catalyst needs to be evaluated according to human experience.

Description

Technical Field

[0001] This disclosure relates to the field of denitrification catalysts, and more particularly to the field of denitrification catalyst lifetime assessment technology. Background Technology

[0002] Currently, the main method for assessing the remaining effective life of denitrification catalysts is through sampling testing. This involves waiting until the denitrification system is shut down, then using hoisting equipment to extract the denitrification catalyst from the denitrification reactor and sending it to a professional testing institution for activity evaluation. After the laboratory activity evaluation is completed, the results are sent back to the user.

[0003] The current method for assessing the remaining effective life of denitrification catalysts, which combines physical sampling, precise laboratory testing, and data modeling, has the following drawbacks: First, it lacks timeliness, meaning users cannot know the remaining effective life of the catalyst in real time, especially when problems are exposed in the reactor, making it impossible to immediately rule out whether the catalyst has deactivated. Second, it is costly, requiring significant time and direct costs, with the total cost of a single sampling test reaching tens of thousands of yuan, and the entire monitoring cycle taking approximately one month, resulting in a substantial cumulative cost burden for users.

[0004] Therefore, how to efficiently and quickly assess the lifespan of denitrification catalysts has become an urgent problem to be solved. Summary of the Invention

[0005] This disclosure provides a method, apparatus, equipment, and storage medium for assessing the lifespan of a denitrification catalyst.

[0006] According to a first aspect of this disclosure, a method for assessing the lifetime of a denitrification catalyst is provided. The method includes: Obtain the current operating parameters of the denitrification system; The flue gas flow rate entering the denitrification system is corrected based on the current operating parameters of the denitrification system to obtain the corrected flue gas flow rate; Based on the corrected flue gas flow rate, calculate the surface velocity of the denitrification catalyst in the denitrification system; The dynamic decay time of the denitration catalyst per unit time is calculated based on the surface velocity. The lifetime of the denitrification catalyst is evaluated based on the dynamic decay time per unit time.

[0007] In addition to the aspects and any possible implementations described above, a further implementation is provided that evaluates the lifetime of the denitration catalyst based on the dynamic decay time per unit time, including: Calculate the denitrification efficiency of the denitrification system under different operating conditions and for different operating times; Based on the denitrification efficiency of the denitrification system under different operating conditions and different operating times, calculate the theoretical activity value under the different operating conditions; Calculate the theoretical decay time under the different operating conditions based on the theoretical activity values ​​under the different operating conditions. The lifetime of the denitrification catalyst is evaluated based on the dynamic decay time per unit time and the theoretical decay time under different operating conditions.

[0008] As described above and in any possible implementation, a further implementation is provided in which the calculation of the denitrification efficiency of the denitrification system under different operating conditions and for different operating times includes: The operating load and operating temperature of the denitrification system under different operating conditions and at different operating times were obtained; Obtain the preset theoretical denitrification efficiency formula; Based on the preset theoretical denitrification efficiency formula and the operating load and operating temperature of the denitrification system under different operating conditions and different operating times, the denitrification efficiency of the denitrification system under different operating conditions and different operating times is calculated. The step of calculating the theoretical activity value under different operating conditions based on the denitrification efficiency of the denitrification system under different operating conditions and different operating times includes: The theoretical activity values ​​under different operating conditions are calculated based on the denitrification efficiency of the denitrification system under different operating conditions and operating times, the preset value of ammonia slip, the concentration of nitrogen oxides at the inlet of the denitrification system, and the surface velocity of the denitrification catalyst.

[0009] In addition to the aspects and any possible implementations described above, a further implementation is provided, wherein the preset theoretical denitrification efficiency formula includes:

[0010] This represents the theoretical denitrification efficiency after operating time t under different operating conditions. and These represent the operating load and operating temperature under different working conditions, respectively. T represents the theoretical denitrification efficiency. B This indicates the flue gas temperature under BMCR operating conditions. This represents the operating load under BMCR conditions, where α, β, and γ represent correction factors for operating load, operating temperature, and operating time, respectively, and t is time.

[0011] In addition to the aspects and any possible implementations described above, a further implementation is provided, which calculates the theoretical decay time under the different operating conditions based on the theoretical activity values ​​under the different operating conditions, including: The design life and stable operating time of the denitrification catalyst were obtained. The theoretical activity values ​​of the denitrification catalyst at its designed life and at its stable operating time under different operating conditions were obtained. Based on the theoretical activity values ​​of the denitrification catalyst during its designed lifespan and the theoretical activity values ​​during its stable operating time under different operating conditions, the theoretical decay rate of the denitrification catalyst under different operating conditions is calculated. The attenuation rate weighting coefficient is determined based on the theoretical attenuation rate of the denitrification catalyst under full load conditions. The theoretical decay time of the denitrification catalyst under different operating conditions is calculated based on the theoretical decay rate of the denitrification catalyst under different operating conditions and the decay rate weighting coefficient.

[0012] In addition to the aspects and any possible implementations described above, a further implementation is provided, wherein correcting the flue gas flow rate entering the denitrification system based on the current operating parameters of the denitrification system to obtain a corrected flue gas flow rate includes: Call the preset flue gas flow rate conversion formula; The inlet flue gas flow rate of the denitrification system is calculated based on the current operating parameters of the denitrification system and the preset flue gas flow rate conversion formula. Call the preset flue gas flow correction formula; Based on the converted flue gas flow rate at the inlet of the denitrification system, the flue gas flow rate reading at the inlet of the denitrification system, and the preset flue gas flow rate correction formula, the corrected flue gas flow rate is obtained, wherein the preset flue gas flow rate conversion formula is:

[0013] This represents the inlet flue gas flow rate calculated based on operational conditions at the i-th data acquisition time. , , These represent the inlet flue gas temperature reading, inlet flue gas pressure reading, and inlet operating load reading of the denitrification system at the i-th data acquisition time, respectively. , , , These represent the inlet flue gas temperature, flue gas pressure, operating load, and flue gas flow rate readings of the denitrification system at the i-th data acquisition time under BMCR operating conditions.

[0014] In addition to the aspects and any possible implementations described above, a further implementation is provided in which the preset flue gas flow correction formula includes:

[0015] in, This represents the corrected flue gas flow rate at the i-th data acquisition time. This represents the inlet flue gas flow rate reading of the denitrification system at the i-th data acquisition time, and ε represents the correction coefficient, which is calculated as follows:

[0016] Where ε1 represents the deviation factor, ε2 represents the time series factor, and ε3 represents the operating condition factor:

[0017]

[0018]

[0019] Where Qread(i) represents the inlet flue gas flow rate reading of the denitrification system at the i-th data acquisition time, and Qread(max) and Qread(min) represent the maximum and minimum values ​​of the inlet flue gas flow rate reading of the denitrification system within n data acquisition times, respectively.

[0020] In addition to the aspects and any possible implementations described above, a further implementation is provided, which calculates the dynamic decay time of the denitrification catalyst per unit time based on the surface velocity, including: Obtain the nitrogen oxide concentration readings at the inlet and outlet of the denitrification system; The real-time denitrification efficiency of the denitrification system is calculated based on the nitrogen oxide concentration readings at the inlet and outlet of the denitrification system. Obtain the preset value of ammonia slip in the denitrification system; The ammonia-nitrogen molar ratio of the denitrification system is calculated based on the preset ammonia slip value and the real-time denitrification efficiency of the denitrification system. The current real-time activity of the denitrification catalyst is calculated based on the ammonia-nitrogen molar ratio of the denitrification system, the real-time denitrification efficiency of the denitrification system, and the surface velocity. Obtain the current theoretical activity of the denitrification system; The dynamic decay time of the denitrification catalyst per unit time is calculated based on the current real-time activity of the denitrification catalyst and the current theoretical activity of the denitrification system.

[0021] According to a second aspect of this disclosure, a denitrification catalyst lifetime assessment device is provided. The device includes: The acquisition module is used to acquire the current operating parameters of the denitrification system; The correction module is used to correct the flue gas flow rate entering the denitrification system according to the current operating parameters of the denitrification system, so as to obtain the corrected flue gas flow rate. The first calculation module is used to calculate the surface velocity of the denitrification catalyst in the denitrification system based on the corrected flue gas flow rate. The second calculation module is used to calculate the dynamic decay time of the denitrification catalyst per unit time based on the surface velocity; An evaluation module is used to evaluate the lifetime of the denitrification catalyst based on the dynamic decay time per unit time.

[0022] According to a third aspect of this disclosure, an electronic device is provided. The electronic device includes a memory and a processor, wherein the memory stores a computer program, and the processor executes the program to implement the method described above.

[0023] According to a fourth aspect of this disclosure, a computer-readable storage medium is provided having a computer program stored thereon that, when executed by a processor, implements the method according to a first aspect of this disclosure.

[0024] In this disclosure, after obtaining the current operating parameters of the denitrification system, the flue gas flow rate entering the denitrification system can be corrected based on the current operating parameters to obtain the corrected flue gas flow rate. Then, based on the corrected flue gas flow rate, the surface velocity of the denitrification catalyst in the denitrification system is calculated. Subsequently, the dynamic decay time of the denitrification catalyst per unit time is calculated based on the surface velocity. Finally, the lifetime of the denitrification catalyst is evaluated based on the dynamic decay time per unit time. In this way, the lifetime of the denitrification catalyst can be evaluated efficiently, quickly, and accurately using the current operating parameters of the denitrification system, avoiding the need to evaluate the lifetime of the denitrification catalyst based on human experience. Moreover, this evaluation method does not require sampling and is fast, thus reducing the cost of catalyst lifetime evaluation.

[0025] It should be understood that the description in the Summary of the Invention section is not intended to limit the key or essential features of the embodiments of this disclosure, nor is it intended to restrict the scope of this disclosure. Other features of this disclosure will become readily apparent from the following description. Attached Figure Description

[0026] The above and other features, advantages, and aspects of the embodiments of this disclosure will become more apparent from the accompanying drawings and the following detailed description. The drawings are provided for a better understanding of the invention and are not intended to limit the scope of this disclosure. In the drawings, the same or similar reference numerals denote the same or similar elements, wherein: Figure 1 A flowchart of a method for evaluating the lifetime of a denitrification catalyst according to an embodiment of the present disclosure is shown; Figure 2 A block diagram of a denitrification catalyst lifetime assessment apparatus according to an embodiment of the present disclosure is shown; Figure 3 A block diagram of an exemplary electronic device capable of implementing embodiments of the present disclosure is shown. Detailed Implementation

[0027] To make the objectives, technical solutions, and advantages of the embodiments of this disclosure clearer, the technical solutions of the embodiments of this disclosure will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this disclosure, and not all embodiments. Based on the embodiments of this disclosure, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this disclosure.

[0028] Furthermore, the term "and / or" in this article is merely a description of the relationship between related objects, indicating that three relationships can exist. For example, A and / or B can represent: A existing alone, A and B existing simultaneously, or B existing alone. Additionally, the character " / " in this article generally indicates that the preceding and following related objects have an "or" relationship.

[0029] Figure 1 A flowchart of a denitrification catalyst lifetime assessment method 100 according to an embodiment of the present disclosure is shown. Method 100 may include: Step 110: Obtain the current operating parameters of the denitrification system; To ensure the accuracy and real-time nature of the denitrification catalyst life assessment, it is necessary to read some operating data from the denitrification side of the system at 3-second intervals. The current operating parameters to be read include, but are not limited to: unit load, NOx concentration at the inlet of the denitrification reactor (i.e., the denitrification system), NOx concentration at the outlet of the denitrification reactor, flue gas temperature at the inlet of the denitrification reactor, flue gas flow rate at the inlet of the denitrification reactor, and flue gas pressure at the inlet of the denitrification reactor. Data preprocessing: To avoid the impact of outliers on the evaluation results, thresholds need to be set for each current operating parameter during the data preprocessing stage to directly remove data that obviously exceeds the reasonable range. That is, when a certain current operating parameter exceeds the set range, the corresponding entire data will be deleted. Data verification: Since the denitrification data is continuous, the time correlation of the data is used to identify outliers with short-term sudden changes. For each parameter, the current value is compared with the average of the previous 5 time points. If the deviation exceeds the set ratio of ±30%, the value is determined to be an outlier, and the corresponding entire data is deleted.

[0030] Step 120: Based on the current operating parameters of the denitrification system, correct the flue gas flow rate entering the denitrification system to obtain the corrected flue gas flow rate; Step 130: Calculate the surface velocity of the denitrification catalyst in the denitrification system based on the corrected flue gas flow rate;

[0031] Among them, A V V represents the surface velocity (m / h), and V represents the total volume of the catalyst (m³ / h). 3 A p Geometric specific surface area (m²) 2 / m 3 The geometric specific surface area corresponds to the number of pores in the catalyst. For the corrected flue gas flow rate, A V The corresponding table 1 is as follows: Table 1

[0032] Step 140: Calculate the dynamic decay time of the denitrification catalyst per unit time based on the surface velocity; Step 150: Evaluate the lifetime of the denitrification catalyst based on the dynamic decay time per unit time.

[0033] After obtaining the current operating parameters of the denitrification system, the flue gas flow rate entering the system can be corrected based on these parameters to obtain the corrected flue gas flow rate. Then, based on the corrected flue gas flow rate, the surface velocity of the denitrification catalyst in the denitrification system can be calculated. Subsequently, the dynamic decay time of the denitrification catalyst per unit time can be calculated based on the surface velocity. Finally, the lifespan of the denitrification catalyst can be evaluated based on the dynamic decay time per unit time. In this way, the lifespan of the denitrification catalyst can be evaluated efficiently, quickly, and accurately using the current operating parameters of the denitrification system, avoiding the need to evaluate the lifespan of the denitrification catalyst based on human experience. Moreover, this evaluation method does not require sampling and is fast, thus reducing the cost of catalyst lifespan evaluation.

[0034] In some embodiments, the lifetime of the denitration catalyst is evaluated based on the dynamic decay time per unit time, including: Calculate the denitrification efficiency of the denitrification system under different operating conditions and for different operating times; Based on the denitrification efficiency of the denitrification system under different operating conditions and different operating times, calculate the theoretical activity value under the different operating conditions; Calculate the theoretical decay time under the different operating conditions based on the theoretical activity values ​​under the different operating conditions. The lifetime of the denitrification catalyst is evaluated based on the dynamic decay time per unit time and the theoretical decay time under different operating conditions.

[0035] After calculating the denitrification efficiency (the theoretical denitrification efficiency) of the denitrification system under different operating conditions and for different operating times, the theoretical activity value under different operating conditions can be calculated based on the denitrification efficiency of the denitrification system under different operating conditions and for different operating times. Then, based on the theoretical activity value under different operating conditions, the theoretical decay time of the catalyst under different operating conditions can be calculated. Furthermore, based on the dynamic decay time per unit time and the theoretical decay time under different operating conditions, the lifetime of the denitrification catalyst can be quickly and efficiently evaluated.

[0036] In some embodiments, calculating the denitrification efficiency of the denitrification system under different operating conditions and for different operating times includes: The operating load and operating temperature of the denitrification system under different operating conditions and at different operating times were obtained; Obtain the preset theoretical denitrification efficiency formula; Based on the preset theoretical denitrification efficiency formula and the operating load and operating temperature of the denitrification system under different operating conditions and different operating times, the denitrification efficiency of the denitrification system under different operating conditions and different operating times is calculated. The step of calculating the theoretical activity value under different operating conditions based on the denitrification efficiency of the denitrification system under different operating conditions and different operating times includes: The theoretical activity values ​​under different operating conditions are calculated based on the denitrification efficiency of the denitrification system under different operating conditions and operating times, the preset value of ammonia slip, the concentration of nitrogen oxides at the inlet of the denitrification system, and the surface velocity of the denitrification catalyst.

[0037] By substituting the operating load and operating temperature of the denitrification system under different operating conditions and operating times into the aforementioned preset theoretical denitrification efficiency formula, the denitrification efficiency of the denitrification system under different operating conditions and operating times can be accurately calculated. Furthermore, based on the denitrification efficiency of the denitrification system under different operating conditions and operating times, the preset ammonia slip value, the inlet nitrogen oxide concentration of the denitrification system, and the surface velocity of the denitrification catalyst, the theoretical activity value under different operating conditions can be accurately and quickly calculated.

[0038] In some embodiments, the preset theoretical denitrification efficiency formula includes:

[0039] This represents the theoretical denitrification efficiency after operating time t under different operating conditions. and These represent the operating load and operating temperature under different working conditions, respectively. T represents the theoretical denitrification efficiency. B This indicates the flue gas temperature under BMCR operating conditions. This represents the operating load under BMCR conditions, where α, β, and γ represent correction factors for operating load, operating temperature, and operating time, respectively, and t is time.

[0040] In some embodiments, calculating the theoretical decay time under different operating conditions based on the theoretical activity values ​​under different operating conditions includes: The design life and stable operating time of the denitrification catalyst were obtained. The theoretical activity values ​​of the denitrification catalyst at its designed life and at its stable operating time under different operating conditions were obtained. For example, if a certain denitrification catalyst has a designed lifespan of 24,000 hours and its stable operating time is 4,400 hours, then we need to calculate the theoretical activity value (i.e., K) of the denitrification catalyst under each operating condition at 24,000 hours. 24000 * ) and the theoretical activity value at 4400 hours (i.e., K) 4400 * (as shown in Table 2).

[0041] Table 2

[0042] Based on the theoretical activity values ​​of the denitrification catalyst during its designed lifespan and the theoretical activity values ​​during its stable operating time under different operating conditions, the theoretical decay rate of the denitrification catalyst under different operating conditions is calculated. Using the above theoretical activity value as a reference, the theoretical attenuation rate θ under different operating loads is calculated (i.e., 6 operating loads correspond to 6 working conditions, and m ranges from 1 to 6, resulting in 6 values). The calculation formula is as follows:

[0043] The attenuation rate weighting coefficient is determined based on the theoretical attenuation rate of the denitrification catalyst under full load conditions. Full load operation means operating at 100% load. The operating load is the heat supply from the boiler connected upstream of the denitrification system.

[0044] Referring to the table above, the attenuation rate weighting coefficient = 1 / catalyst design life / The attenuation rate weighting coefficient k = 1 / 24000 / .

[0045] The theoretical decay time of the denitrification catalyst under different operating conditions is calculated based on the theoretical decay rate of the denitrification catalyst under different operating conditions and the decay rate weighting coefficient.

[0046] Theoretical decay time of denitrification catalyst under the different operating conditions , = N represents the total number of operating conditions, and 24000 represents the design life of the catalyst.

[0047] In some embodiments, correcting the flue gas flow rate entering the denitrification system based on the current operating parameters of the denitrification system to obtain the corrected flue gas flow rate includes: Call the preset flue gas flow rate conversion formula; The inlet flue gas flow rate of the denitrification system is calculated based on the current operating parameters of the denitrification system and the preset flue gas flow rate conversion formula. Call the preset flue gas flow correction formula; Based on the converted flue gas flow rate at the inlet of the denitrification system, the flue gas flow rate reading at the inlet of the denitrification system, and the preset flue gas flow rate correction formula, the corrected flue gas flow rate is obtained, wherein the preset flue gas flow rate conversion formula is:

[0048] This represents the inlet flue gas flow rate calculated based on operational conditions at the i-th data acquisition time. , , These represent the inlet flue gas temperature reading, inlet flue gas pressure reading, and inlet operating load reading of the denitrification system at the i-th data acquisition time, respectively. , , , These represent the inlet flue gas temperature, flue gas pressure, operating load, and flue gas flow rate readings of the denitrification system at the i-th data acquisition time under BMCR operating conditions.

[0049] In some embodiments, the preset flue gas flow correction formula includes:

[0050] in, This represents the corrected flue gas flow rate at the i-th data acquisition time. This represents the inlet flue gas flow rate reading of the denitrification system at the i-th data acquisition time, and ε represents the correction coefficient, which is calculated as follows:

[0051] Where ε1 represents the deviation factor, ε2 represents the time series factor, and ε3 represents the operating condition factor:

[0052]

[0053]

[0054] Where Qread(i) represents the inlet flue gas flow rate reading of the denitrification system at the i-th data acquisition time, and Qread(max) and Qread(min) represent the maximum and minimum values ​​of the inlet flue gas flow rate reading of the denitrification system within n data acquisition times, respectively.

[0055] In some embodiments, calculating the dynamic decay time of the denitrification catalyst per unit time based on the surface velocity includes: Obtain the nitrogen oxide concentration readings at the inlet and outlet of the denitrification system; The real-time denitrification efficiency of the denitrification system is calculated based on the nitrogen oxide concentration readings at the inlet and outlet of the denitrification system. Real-time denitrification efficiency :

[0056] Obtain the preset value of ammonia slip in the denitrification system; The ammonia-nitrogen molar ratio of the denitrification system is calculated based on the preset ammonia slip value and the real-time denitrification efficiency of the denitrification system. Ammonia-nitrogen molar ratio MR:

[0057] This is the preset value for ammonia escape.

[0058] The current real-time activity of the denitrification catalyst is calculated based on the ammonia-nitrogen molar ratio of the denitrification system, the real-time denitrification efficiency of the denitrification system, and the surface velocity. Real-time activity calculation: The current real-time activity of the catalyst is calculated according to the following formula:

[0059] Where Kt represents the current real-time activity, and MR represents the ammonia-nitrogen molar ratio. Indicates real-time denitrification efficiency:

[0060]

[0061] Among them, NH3 des This indicates the guaranteed ammonia slip value (ppm) (the guaranteed value is a preset value, usually 3), NOx in and NOX outThese represent the NOx readings (mg / m³) at the inlet and outlet of the denitrification reactor, respectively. 3 ).

[0062] Obtain the current theoretical activity of the denitrification system; Theoretical activity calculation: Under current operating conditions, the theoretical activity value of the catalyst needs to be calculated. By comparing the theoretical value with the actual value, the relative activity curve can be obtained.

[0063] Where K0 represents the current theoretical activity, and η0 represents the theoretical denitrification efficiency. When calculating the theoretical activity, the ammonia-nitrogen molar ratio MR is constant at 1. Since the guaranteed ammonia slip value NH3des and the inlet NOX reading NOXin are known, the theoretical denitrification efficiency can be deduced from the formula:

[0064] The dynamic decay time of the denitrification catalyst per unit time is calculated based on the current real-time activity of the denitrification catalyst and the current theoretical activity of the denitrification system.

[0065] Dynamic decay time of denitrification catalyst per unit time :

[0066] Let t represent the current real-time activity of the denitration catalyst. Let t be the current theoretical activity of the denitration catalyst. This represents the current real-time activity of the denitration catalyst at time t-1. This represents the current theoretical activity of the denitration catalyst at time t-1. This indicates the catalyst design life, and γ represents the dynamic decay time adjustment parameter, which defaults to 240.

[0067] This invention is a method for evaluating the lifetime of denitrification catalysts based on a hyperbolic decay curve, comprising: 1. Denitrification data acquisition and preprocessing (1) Data acquisition: In order to ensure the accuracy and real-time nature of the denitrification catalyst life assessment, it is necessary to read some operating data of the denitrification side from the user's operating system. The data is read at 3-second intervals. The data includes: unit load, NOx concentration at the inlet of the denitrification reactor, NOx concentration at the outlet of the denitrification reactor, flue gas temperature at the inlet of the denitrification reactor, flue gas flow rate at the inlet of the denitrification reactor, and flue gas pressure at the inlet of the denitrification reactor. (2) Data preprocessing: In order to avoid the impact of outliers on the evaluation results, hard thresholds need to be set for each reading during the data preprocessing stage to directly remove data that obviously exceeds the reasonable range. That is, when a parameter exceeds the set range, the corresponding entire data will be deleted. (3) Data verification: Since the denitrification data is continuous, the time correlation of the data is used to identify outliers with short-term sudden changes. For each parameter, the current value is compared with the average value of the previous 5 times. If the deviation exceeds the set ratio ±30%, the value is determined to be an outlier and the corresponding whole data is deleted.

[0068] 2. Calculation of dynamic decay time (1) Flue gas flow rate correction: Due to uncertainties in the field environment, the readings of the flue gas flow sensor may not be accurate enough. Therefore, it is necessary to correct the flue gas flow rate by taking into account the flue gas temperature and operating load. The formula is as follows:

[0069] in, This represents the inlet flue gas flow rate calculated based on operational conditions at the i-th data acquisition time. , , These represent the inlet flue gas temperature reading, inlet flue gas pressure reading, and inlet operating load reading of the denitrification system at the i-th data acquisition time, respectively. , , , These represent the inlet flue gas temperature, flue gas pressure, operating load, and flue gas flow rate readings of the denitrification system at the i-th data acquisition time under BMCR operating conditions.

[0070] The flue gas flow rate is corrected by combining the converted flue gas flow rate and the flue gas flow rate reading. The correction formula is as follows:

[0071] in, This represents the corrected flue gas flow rate at the i-th data acquisition time. This represents the inlet flue gas flow rate reading of the denitrification system at the i-th data acquisition time, and ε represents the correction coefficient, which is calculated as follows:

[0072] Where ε1 represents the deviation factor, ε2 represents the time series factor, and ε3 represents the operating condition factor:

[0073]

[0074]

[0075] Where Qread(i) represents the flue gas flow rate reading at time i, and Qread(max) and Qread(min) represent the maximum and minimum values ​​of the flue gas flow rate reading at time i, respectively.

[0076] (2) Surface velocity calculation: After the flue gas flow rate is corrected, the surface velocity can be calculated in conjunction with the catalyst design parameters. The calculation formula is as follows:

[0077] Among them, A V V represents the surface velocity (m / h), and V represents the total volume of the catalyst (m³ / h). 3 A p Geometric specific surface area (m²) 2 / m 3 The geometric specific surface area corresponds to the number of catalyst pores, as shown in Table 1 below: Table 1

[0078] Real-time activity calculation: The current real-time activity of the catalyst is calculated according to the following formula:

[0079] Where Kt represents the current real-time activity, MR represents the ammonia-nitrogen molar ratio, and ηt represents the real-time denitrification efficiency:

[0080]

[0081] Wherein, NH3des represents the guaranteed ammonia slip value (ppm), and NOXin and NOXout represent the NOX readings (mg / m3) at the inlet and outlet of the denitrification reactor, respectively.

[0082] Theoretical activity calculation: Under current operating conditions, the theoretical activity value of the catalyst needs to be calculated. By comparing the theoretical value with the actual value, the relative activity curve can be obtained.

[0083] Where K0 represents the current theoretical activity, and η0 represents the theoretical denitrification efficiency. When calculating the theoretical activity, the ammonia-nitrogen molar ratio MR is constant at 1, due to the guaranteed ammonia slip value NH3. des Given the NOx reading (NOXin) at the outlet, the theoretical denitrification efficiency can be deduced using the formula:

[0084] Dynamic decay time calculation: Due to the long change cycle of the catalyst, the dynamic decay time is calculated every hour. The calculation formula is as follows:

[0085] Where tdy represents the dynamic decay time per hour, tdes represents the catalyst design life, and γ represents the dynamic decay time adjustment parameter, which defaults to 240. 3. Calculation of theoretical decay time (1) Based on the design parameters of the denitrification reactor and the denitrification catalyst, the denitrification efficiency under different operating conditions and for different operating times can be calculated. The calculation method is as follows:

[0086] Where ηt* represents the theoretical denitrification efficiency after running for time t under the current operating conditions, L* and T* represent the operating load and temperature under the current operating conditions, respectively, and α, β, and γ represent the correction coefficients for operating load, temperature, and operating time, respectively, with default values ​​of 0.2, 40, and 0.05. (2) The activity after different operating times can be calculated based on the denitrification efficiency at different operating times:

[0087]

[0088] (3) Under normal circumstances, it is necessary to run the activity values ​​for 4400 hours and 24000 hours respectively under operating loads of 30%, 50%, 75%, 90%, 100%, and 110%, and calculate the corresponding K / K0 values, as shown in the table below.

[0089] After obtaining the corresponding tables, the theoretical attenuation rate θ under different operating loads can be calculated (i.e., 6 operating loads correspond to 6 working conditions, and the value of m ranges from 1 to 6, resulting in 6 values). The calculation formula is as follows:

[0090] Since the denitrification reactor operates for 1 hour at 100% load, the ideal lifespan of the denitrification catalyst should be reduced by 1 hour. Therefore, the weight k is calculated based on 100% load (i.e., m=5 when referring to the table above):

[0091] The weighted attenuation rate θ° is obtained by weighting each theoretical attenuation rate. = (where m ranges from 1 to 6), the theoretical decay time tth under different operating loads can be calculated:

[0092] (4) Curve fitting: By iteratively calculating the operating load minus the theoretical decay time using the least squares method, an exponential fitting function can be obtained. That is, a fitting function with the operating load as the independent variable and the theoretical decay time as the dependent variable. By substituting the current operating load reading into the function, the theoretical decay time under the current operating load can be obtained. .

[0093] 4. Decay Time Calculation: Finally, by weighting and combining the dynamic decay time with the theoretical decay time, the actual lifespan decay time (s) corresponding to 1 hour of operation under different working conditions can be obtained:

[0094] For example, if it runs for 1 hour, the theoretical lifespan decay time is 1 hour, but the actual lifespan decay time can be 40 minutes or 1 hour and 10 minutes, etc.

[0095] a and b represent the weights of the theoretical decay time and the dynamic decay time, respectively, with default values ​​of 0.85 and 0.15.

[0096] Thus, under different operating conditions, after operating time t, the remaining life of the catalyst = The value of p ranges from 0 to t. (p) represents the actual lifetime decay time in the p-th unit of time.

[0097] The denitration catalyst lifetime assessment method based on the hyperbola of decay developed in this invention can quantify the catalyst performance decay mechanism and provide users with sufficient time to deal with problems. When calculating the dynamic decay time, this invention uses design parameters to correct the operating conditions, which can effectively eliminate the influence of the dispersion of operating data on the results. The method of combining dynamic decay time and theoretical decay time developed in this invention can effectively characterize the decay law of the catalyst and effectively avoid the problems of being affected by short-term extreme operating conditions or being out of touch with actual operating scenarios.

[0098] It should be noted that, for the sake of simplicity, the foregoing method embodiments are all described as a series of actions. However, those skilled in the art should understand that this disclosure is not limited to the described order of actions, because according to this disclosure, some steps can be performed in other orders or simultaneously. Furthermore, those skilled in the art should also understand that the embodiments described in the specification are all optional embodiments, and the actions and modules involved are not necessarily essential to this disclosure.

[0099] The above is an introduction to the method embodiments. The following describes the solution described in this disclosure further through device embodiments.

[0100] Figure 2 A block diagram of a denitrification catalyst lifetime assessment device 200 according to an embodiment of the present disclosure is shown. Figure 2 As shown, the device 200 includes: Module 210 is used to acquire the current operating parameters of the denitrification system; The correction module 220 is used to correct the flue gas flow rate entering the denitrification system according to the current operating parameters of the denitrification system, so as to obtain the corrected flue gas flow rate. The first calculation module 230 is used to calculate the surface velocity of the denitrification catalyst in the denitrification system based on the corrected flue gas flow rate. The second calculation module 240 is used to calculate the dynamic decay time of the denitrification catalyst per unit time based on the surface velocity; Evaluation module 250 is used to evaluate the lifetime of the denitrification catalyst based on the dynamic decay time per unit time.

[0101] Those skilled in the art will understand that, for the sake of convenience and brevity, the specific working process of the described module can be referred to the corresponding process in the foregoing method embodiments, and will not be repeated here.

[0102] According to embodiments of the present disclosure, the present disclosure also provides an electronic device and a non-transitory computer-readable storage medium storing computer instructions.

[0103] Figure 3 A schematic block diagram of an electronic device 800 that can be used to implement embodiments of the present disclosure is shown. The electronic device is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device may also represent various forms of mobile devices, such as personal digital processors, cellular phones, smartphones, wearable devices, and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely illustrative and are not intended to limit the implementation of the present disclosure described and / or claimed herein.

[0104] Device 800 includes a computing unit 801, which can perform various appropriate actions and processes based on a computer program stored in read-only memory (ROM) 802 or a computer program loaded from storage unit 808 into random access memory (RAM) 803. RAM 803 may also store various programs and data required for the operation of device 800. The computing unit 801, ROM 802, and RAM 803 are interconnected via bus 804. Input / output (I / O) interface 805 is also connected to bus 804.

[0105] Multiple components in device 800 are connected to I / O interface 805, including: input unit 806, such as keyboard, mouse, etc.; output unit 807, such as various types of monitors, speakers, etc.; storage unit 808, such as disk, optical disk, etc.; and communication unit 809, such as network card, modem, wireless transceiver, etc. Communication unit 809 allows device 800 to exchange information / data with other devices through computer networks such as the Internet and / or various telecommunications networks.

[0106] The computing unit 801 can be a variety of general-purpose and / or special-purpose processing components with processing and computing capabilities. Some examples of the computing unit 801 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various special-purpose artificial intelligence (AI) computing chips, various computing units running machine learning model algorithms, a digital signal processor (DSP), and any suitable processor, controller, microcontroller, etc. The computing unit 801 performs the various methods and processes described above, such as method 100. For example, in some embodiments, method 100 may be implemented as a computer software program tangibly contained in a machine-readable medium, such as storage unit 808. In some embodiments, part or all of the computer program may be loaded and / or installed on device 800 via ROM 802 and / or communication unit 809. When the computer program is loaded into RAM 803 and executed by the computing unit 801, one or more steps of method 100 described above may be performed. Alternatively, in other embodiments, the computing unit 801 may be configured to perform method 100 by any other suitable means (e.g., by means of firmware).

[0107] Various embodiments of the systems and techniques described above herein can be implemented in digital electronic circuit systems, integrated circuit systems, field-programmable gate arrays (FPGAs), application-specific integrated circuits (ASICs), application-specific standard products (ASSPs), systems-on-a-chip (SoCs), payload-programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various embodiments may include implementations in one or more computer programs that can be executed and / or interpreted on a programmable system including at least one programmable processor, which may be a dedicated or general-purpose programmable processor, capable of receiving data and instructions from a storage system, at least one input device, and at least one output device, and transmitting data and instructions to the storage system, the at least one input device, and the at least one output device.

[0108] The program code used to implement the methods of this disclosure may be written in any combination of one or more programming languages. This program code may be provided to a processor or controller of a general-purpose computer, special-purpose computer, or other programmable data processing apparatus, such that when executed by the processor or controller, the program code causes the functions / operations specified in the flowcharts and / or block diagrams to be implemented. The program code may be executed entirely on a machine, partially on a machine, as a standalone software package partially on a machine and partially on a remote machine, or entirely on a remote machine or server.

[0109] In the context of this disclosure, a machine-readable medium can be a tangible medium that may contain or store a program for use by or in conjunction with an instruction execution system, apparatus, or device. A machine-readable medium can be a machine-readable signal medium or a machine-readable storage medium. A machine-readable medium can be, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, apparatus, or devices, or any suitable combination of the foregoing. More specific examples of machine-readable storage media include electrical connections based on one or more wires, portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fiber, portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination of the foregoing.

[0110] To provide interaction with a user, the systems and techniques described herein can be implemented on a computer having: a display device for displaying information to the user (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor); and a keyboard and pointing device (e.g., a mouse or trackball) through which the user provides input to the computer. Other types of devices can also be used to provide interaction with the user; for example, feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form (including sound input, voice input, or tactile input).

[0111] The systems and technologies described herein can be implemented in computing systems that include backend components (e.g., as a data server), or computing systems that include middleware components (e.g., an application server), or computing systems that include frontend components (e.g., a user computer with a graphical user interface or web browser through which a user can interact with implementations of the systems and technologies described herein), or any combination of such backend, middleware, or frontend components. The components of the system can be interconnected via digital data communication of any form or medium (e.g., a communication network). Examples of communication networks include local area networks (LANs), wide area networks (WANs), and the Internet.

[0112] Computing systems can include clients and servers. Clients and servers are generally located far apart and typically interact via communication networks. Client-server relationships are created by computer programs running on the respective computers and having a client-server relationship with each other. Servers can be cloud servers, servers in distributed systems, or servers incorporating blockchain technology.

[0113] It should be understood that the various forms of processes shown above can be used to rearrange, add, or delete steps. For example, the steps described in this disclosure can be executed in parallel, sequentially, or in different orders, as long as the desired result of the technical solution disclosed in this disclosure can be achieved, and this is not limited herein.

[0114] The specific embodiments described above do not constitute a limitation on the scope of protection of this disclosure. Those skilled in the art should understand that various modifications, combinations, sub-combinations, and substitutions can be made according to design requirements and other factors. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of this disclosure should be included within the scope of protection of this disclosure.

Claims

1. A method for evaluating the lifespan of a denitrification catalyst, characterized in that, include: Obtain the current operating parameters of the denitrification system; The flue gas flow rate entering the denitrification system is corrected based on the current operating parameters of the denitrification system to obtain the corrected flue gas flow rate; Based on the corrected flue gas flow rate, calculate the surface velocity of the denitrification catalyst in the denitrification system; The dynamic decay time of the denitration catalyst per unit time is calculated based on the surface velocity. The lifetime of the denitrification catalyst is evaluated based on the dynamic decay time per unit time.

2. The method as described in claim 1, characterized in that, The lifetime of the denitration catalyst is evaluated based on the dynamic decay time per unit time, including: Calculate the denitrification efficiency of the denitrification system under different operating conditions and for different operating times; Based on the denitrification efficiency of the denitrification system under different operating conditions and different operating times, calculate the theoretical activity value under the different operating conditions; Calculate the theoretical decay time under the different operating conditions based on the theoretical activity values ​​under the different operating conditions. The lifetime of the denitrification catalyst is evaluated based on the dynamic decay time per unit time and the theoretical decay time under different operating conditions.

3. The method as described in claim 2, characterized in that, The calculation of the denitrification efficiency of the denitrification system under different operating conditions and for different operating times includes: The operating load and operating temperature of the denitrification system under different operating conditions and at different operating times were obtained; Obtain the preset theoretical denitrification efficiency formula; Based on the preset theoretical denitrification efficiency formula and the operating load and operating temperature of the denitrification system under different operating conditions and different operating times, the denitrification efficiency of the denitrification system under different operating conditions and different operating times is calculated. The step of calculating the theoretical activity value under different operating conditions based on the denitrification efficiency of the denitrification system under different operating conditions and different operating times includes: The theoretical activity values ​​under different operating conditions are calculated based on the denitrification efficiency of the denitrification system under different operating conditions and operating times, the preset value of ammonia slip, the concentration of nitrogen oxides at the inlet of the denitrification system, and the surface velocity of the denitrification catalyst.

4. The method as described in claim 3, characterized in that, The preset theoretical denitrification efficiency formula includes: This represents the theoretical denitrification efficiency after operating time t under different operating conditions. and These represent the operating load and operating temperature under different working conditions. T represents the theoretical denitrification efficiency. B This indicates the flue gas temperature under BMCR operating conditions. This represents the operating load under BMCR conditions, where α, β, and γ represent correction factors for operating load, operating temperature, and operating time, respectively, and t is time.

5. The method as described in claim 2, characterized in that, Based on the theoretical activity values ​​under the different operating conditions, the theoretical decay time under the different operating conditions is calculated, including: The design life and stable operating time of the denitrification catalyst were obtained. The theoretical activity values ​​of the denitrification catalyst at its designed life and at its stable operating time under different operating conditions were obtained. Based on the theoretical activity values ​​of the denitrification catalyst during its designed lifespan and the theoretical activity values ​​during its stable operating time under different operating conditions, the theoretical decay rate of the denitrification catalyst under different operating conditions is calculated. The attenuation rate weighting coefficient is determined based on the theoretical attenuation rate of the denitrification catalyst under full load conditions. The theoretical decay time of the denitrification catalyst under different operating conditions is calculated based on the theoretical decay rate of the denitrification catalyst under different operating conditions and the decay rate weighting coefficient.

6. The method as described in claim 1, characterized in that, The step of correcting the flue gas flow rate entering the denitrification system based on the current operating parameters of the denitrification system to obtain the corrected flue gas flow rate includes: Call the preset flue gas flow rate conversion formula; The inlet flue gas flow rate of the denitrification system is calculated based on the current operating parameters of the denitrification system and the preset flue gas flow rate conversion formula. Call the preset flue gas flow correction formula; Based on the converted flue gas flow rate at the inlet of the denitrification system, the flue gas flow rate reading at the inlet of the denitrification system, and the preset flue gas flow rate correction formula, the corrected flue gas flow rate is obtained, wherein the preset flue gas flow rate conversion formula is: This represents the inlet flue gas flow rate calculated based on operational conditions at the i-th data acquisition time. , , These represent the inlet flue gas temperature reading, inlet flue gas pressure reading, and inlet operating load reading of the denitrification system at the i-th data acquisition time, respectively. , , , These represent the inlet flue gas temperature, flue gas pressure, operating load, and flue gas flow rate readings of the denitrification system at the i-th data acquisition time under BMCR operating conditions.

7. The method as described in claim 6, characterized in that, The preset flue gas flow correction formula includes: in, This represents the corrected flue gas flow rate (i.e., the corrected flue gas flow rate) at the i-th data acquisition time. This represents the inlet flue gas flow rate reading of the denitrification system at the i-th data acquisition time, and ε represents the correction coefficient, which is calculated as follows: Where ε1 represents the deviation factor, ε2 represents the time series factor, and ε3 represents the operating condition factor: Where Qread(i) represents the inlet flue gas flow rate reading of the denitrification system at the i-th data acquisition time, and Qread(max) and Qread(min) represent the maximum and minimum values ​​of the inlet flue gas flow rate reading of the denitrification system within n data acquisition times, respectively.

8. The method according to any one of claims 1 to 7, characterized in that, The dynamic decay time of the denitrification catalyst per unit time is calculated based on the surface velocity, including: Obtain the nitrogen oxide concentration readings at the inlet and outlet of the denitrification system; The real-time denitrification efficiency of the denitrification system is calculated based on the nitrogen oxide concentration readings at the inlet and outlet of the denitrification system. Obtain the preset value of ammonia slip in the denitrification system; The ammonia-nitrogen molar ratio of the denitrification system is calculated based on the preset ammonia slip value and the real-time denitrification efficiency of the denitrification system. The current real-time activity of the denitrification catalyst is calculated based on the ammonia-nitrogen molar ratio of the denitrification system, the real-time denitrification efficiency of the denitrification system, and the surface velocity. Obtain the current theoretical activity of the denitrification system; The dynamic decay time of the denitrification catalyst per unit time is calculated based on the current real-time activity of the denitrification catalyst and the current theoretical activity of the denitrification system.

9. A device for evaluating the lifespan of a denitrification catalyst, characterized in that, include: The acquisition module is used to acquire the current operating parameters of the denitrification system; The correction module is used to correct the flue gas flow rate entering the denitrification system according to the current operating parameters of the denitrification system, so as to obtain the corrected flue gas flow rate. The first calculation module is used to calculate the surface velocity of the denitrification catalyst in the denitrification system based on the corrected flue gas flow rate. The second calculation module is used to calculate the dynamic decay time of the denitrification catalyst per unit time based on the surface velocity; An evaluation module is used to evaluate the lifetime of the denitrification catalyst based on the dynamic decay time per unit time.

10. An electronic device, characterized in that, include: Memory and processor The memory stores a computer program, and when the processor executes the program, it implements the method as described in any one of claims 1-8.

Images