Method for calculating non-uniform heat absorption coefficient of high-temperature heating surface of boiler based on monitoring data

Abstract

This invention discloses a method for calculating the heat absorption non-uniformity coefficient of a boiler's high-temperature heating surface based on monitoring data. It relates to the field of heat absorption non-uniformity coefficient calculation and solves the problem of inaccurate calculation of the heat absorption non-uniformity coefficient of a boiler's high-temperature heating surface. The method includes acquiring real-time monitoring data of the boiler, calculating steam parameters at the inlet and outlet of the same tube coil in different tube panels based on the real-time monitoring data, calculating the mass flow rate of steam at the inlet of the distribution header and the outlet of the collection header, calculating the enthalpy increase, enthalpy increase deviation, and average specific volume of steam in the same tube coil in different tube panels based on the steam enthalpy and specific volume at the inlet and outlet, calculating the static pressure difference at the inlet and outlet of the same tube coil in different tube panels, calculating the steam flow rate deviation in the same tube coil in different tube panels, and calculating the heat absorption non-uniformity coefficient of the same tube coil in different tube panels based on the steam flow rate deviation and the steam enthalpy increase deviation. This invention achieves accurate calculation of the heat absorption non-uniformity coefficient of a boiler's high-temperature heating surface.

Description

Technical Field

[0001] This invention belongs to the field of heat absorption non-uniformity coefficient calculation, specifically a method for calculating the heat absorption non-uniformity coefficient of a boiler high-temperature heating surface based on monitoring data. Background Technology

[0002] In a П-type boiler, screen-type high-temperature heating surfaces such as rear superheaters, high-temperature (final stage) superheaters, and high-temperature (final stage) reheaters are generally arranged above the furnace and in the width direction (i.e., left-right direction) of the horizontal flue. Due to factors such as combustion method, fuel, burner input method, and furnace type, there is a significant non-uniformity in the distribution of flue gas temperature and velocity along the width direction of the furnace, and this non-uniformity will continue to varying degrees in the convection flue.

[0003] In existing literature, there are few calculations on the heat absorption non-uniformity coefficient of high-temperature heating surfaces in boilers. It is usually obtained from numerical simulation or commissioning test data. Both of these methods are relatively complex. The accuracy of the former depends on the accuracy of the modeling, while the latter is time-consuming and labor-intensive and cannot achieve the effect of online monitoring. At the same time, some boiler manufacturers' technical documents often directly give the maximum heat absorption non-uniformity coefficient. For example, a boiler manufacturer's technical document states that the furnace outlet heat load non-uniformity coefficient of a tangentially circular boiler is generally controlled below 1.4, but in actual operation, it may reach 1.5 to 1.6 due to residual airflow rotation. This direct reference method does not take into account the actual combustion and operation of the boiler and is only suitable for rough estimation.

[0004] To address this, the present invention proposes a method for calculating the heat absorption non-uniformity coefficient of the high-temperature heating surface of a boiler based on monitoring data. Summary of the Invention

[0005] The purpose of this invention is to propose a method for calculating the heat absorption non-uniformity coefficient of a boiler's high-temperature heating surface based on monitoring data, so as to solve the problem mentioned in the background art of being unable to accurately calculate the heat absorption non-uniformity coefficient of a boiler's high-temperature heating surface.

[0006] The objective of this invention can be achieved through the following technical solutions:

[0007] The calculation method for the uneven heat absorption coefficient of the high-temperature heating surface of a boiler based on monitoring data includes:

[0008] Step S1: Obtain real-time monitoring data of the boiler;

[0009] Step S2: Calculate the steam parameters at the inlet and outlet of the same tube coil in different tube panels based on the real-time monitoring data of the boiler;

[0010] Step S3: Calculate the mass flow rate of steam at the inlet of the distribution header and the outlet of the collection header.

[0011] Step S4: Calculate the enthalpy increase, enthalpy increase deviation, and average specific volume of steam in the same tube coil in different tube screens based on the steam enthalpy and steam specific volume at the inlet and outlet.

[0012] Step S5: Calculate the static pressure difference at the inlet and outlet of the same tube coil in different tube panels;

[0013] Step S6: Calculate the steam flow deviation within the same pipe coil in different pipe panels;

[0014] Step S7: Calculate the heat absorption non-uniformity coefficient of the same tube coil in different tube screens based on the steam flow deviation and the steam enthalpy increase deviation.

[0015] Furthermore, the real-time monitoring data in step S1 includes the total steam flow rate F, the steam pressure Pj at the inlet of the distribution header, the steam pressure Pc at the outlet of the collection header, the steam temperature tzj at the left inlet of the distribution header, the steam temperature tyj at the right inlet of the distribution header, and the wall temperature tb above the ceiling at the outlet of the same pipe coil in different pipe panels. i .

[0016] Furthermore, the steam parameters in step S2 are the pressure and steam temperature at the inlet and outlet of the same tube coil in different tube panels.

[0017] Furthermore, the calculation process in step S2 includes the following sub-steps:

[0018] Step S21: Obtain the total number N of tube panels in the width direction of the boiler heating surface and the number of tube rings in a single tube panel;

[0019] Obtain the steam temperature tzj at the left inlet of the distribution header and the steam temperature tyj at the right inlet of the distribution header, and then use the formula tj i =tzj-(tzj-tyj)×(i-1) / (N-1)The steam temperature tj at the inlet of the same tube coil in different tube screens is calculated. i i is the leftmost tube screen number;

[0020] Step S22: Obtain the wall temperature tb above the ceiling at the outlet of the same tube coil in different tube panels. i Through the formula tc i ≈tb i +1 is used to calculate the steam temperature tc at the outlet of the same tube coil in different tube screens. i ;

[0021] Step S23: Obtain the steam pressure Pj at the inlet of the distribution header, and use the formula Pj i The pressure Pj at the inlet of the same tube coil in different tube panels is calculated by ≈Pj. i ;

[0022] Step S24: Obtain the steam pressure Pc at the outlet of the collection tank, and use the formula Pc i The pressure Pc at the outlet of the same tube coil in different tube panels is calculated by ≈Pc. i .

[0023] Furthermore, the calculation process in step S3 is as follows:

[0024] Step S31: Collect the inner diameter Dj of the distribution header and the inner diameter Dc of the collection header to obtain the total steam flow rate F;

[0025] Step S32, using the formula fj=2F / (Π×Dj) 2 The mass flow rate fj of the steam at the inlet of the distribution header is calculated, where π is the coefficient of pi.

[0026] Step S33, using the formula fc=2F / (Π×Dc) 2 The mass flow rate fc of the steam at the outlet of the collection tank is calculated.

[0027] Furthermore, the calculation process in step S4 is as follows:

[0028] Step S41, using formula Hj i =(a1×Pj 2 +a2×Pj+a3)ln(tj i )-(b1×Pj 2 The steam enthalpy Hj at the inlet of the same tube coil in different tube panels is calculated by adding b2×Pj+b3). i Similarly, through the formula Hc i =(a1×Pc 2 +a2×Pc+a3)ln(tc i )-(b1×Pc 2 The steam enthalpy Hc at the outlet of the same tube coil in different tube panels is calculated by adding b2×Pc+b3). i Where a1, a2, a3, b1, b2, and b3 are all coefficients;

[0029] Step S42: Subtract the steam enthalpy at the inlet from the steam enthalpy at the outlet to obtain the enthalpy increase ΔH of the steam in the same tube coil in different tube panels. i ;

[0030] Step S43, using the formula The enthalpy increase deviation (HZP) of steam in the same tube coil in different tube panels was calculated. i ;

[0031] Step S44, using formula vj i =[e1×(Pj / 100) 2+e2×(Pj / 100)+e3]ln(tj i )-[d1×(Pj / 100) 2 The specific steam volume vj at the inlet of the same tube coil in different tube screens is calculated by adding d2×(Pj / 100)+d3]. i Similarly, through the formula vc i =[e1×(Pc / 100) 2 +e2×(Pc / 100)+e3]ln(tc i )-[d1×(Pc / 100) 2 The specific steam volume vc at the outlet of the same tube coil in different tube panels is calculated by [+d2×(Pc / 100)+d3]. i Where d1, d2, d3, e1, e2, and e3 are all coefficients;

[0032] Step S45: Add the steam specific volume at the inlet and the steam specific volume at the outlet, and take the average value to obtain the average steam specific volume v in the same pipe ring in different pipe panels. i .

[0033] Furthermore, the calculation process of the pressure difference at the inlet and outlet of the same tube coil in different tube panels in step S5 includes the following sub-steps:

[0034] Step S51, firstly through the formula The average specific volume of steam at the inlet of the distribution header, vj, is calculated.

[0035] Similarly, through the formula The average specific volume vc of the steam at the outlet of the collection tank was calculated.

[0036] In step S52, the steam velocity wj at the inlet of the distribution header is calculated using the formula wj=vj×fj. Similarly, the steam velocity wc at the outlet of the collection header is calculated using the formula wc=vc×fc.

[0037] Furthermore, the calculation process of the pressure difference at the inlet and outlet of the same tube coil in different tube panels in step S5 also includes the following sub-steps:

[0038] Step S53, using the formula △Pjmax=kj×(wj) 2 The maximum static pressure difference ΔPjmax at the inlet of the distribution header is calculated using / 2vj).

[0039] Similarly, using the formula △Pcmax=kc×(wc) 2 The maximum static pressure difference ΔPcmax at the outlet of the collection tank is calculated using / 2vc; where kj is the pressure variation coefficient of the distribution tank and kc is the pressure variation coefficient of the collection tank.

[0040] Step S54, using the formula △P i =Pj-Pc+(△Pcmax-△Pjmax)×4(i-1)×(iN) / (N-1) 2 The static pressure difference ΔP at the inlet and outlet of the same tube coil in different tube panels was calculated. i .

[0041] Furthermore, the calculation process in step S6 is as follows:

[0042] Through formula Calculate the steam flow deviation (LLP) within the same tube coil in different tube panels. i .

[0043] Furthermore, the calculation process in step S7 is as follows:

[0044] Through formula Q i =LLP i ×HZP i The heat absorption non-uniformity coefficient Q of the same tube coil in different tube screens was calculated. i .

[0045] In summary, due to the adoption of the above technical solution, the beneficial effects of the present invention are:

[0046] 1. This invention first obtains real-time monitoring data of the boiler, then uses the real-time monitoring data of the boiler to calculate the steam parameters at the inlet and outlet of the same tube ring in different tube panels, and then calculates the mass flow rate of steam at the inlet of the distribution header and the outlet of the collection header. This invention realizes the calculation of the mass flow rate of steam at the inlet of the distribution header and the outlet of the collection header.

[0047] 2. This invention also calculates the enthalpy increase, enthalpy increase deviation, and average specific volume of steam in the same tube coil in different tube screens based on the steam enthalpy value and steam specific volume at the inlet and outlet. Then, it calculates the static pressure difference at the inlet and outlet of the same tube coil in different tube screens, calculates the steam flow deviation in the same tube coil in different tube screens, and finally uses the steam flow deviation and steam enthalpy increase deviation to calculate the heat absorption non-uniformity coefficient of the same tube coil in different tube screens. This invention achieves accurate calculation of the heat absorption non-uniformity coefficient of the high-temperature heating surface of the boiler. Attached Figure Description

[0048] To facilitate understanding by those skilled in the art, the present invention will be further described below with reference to the accompanying drawings.

[0049] Figure 1 This is a flowchart of the method of the present invention;

[0050] Figure 2 This is a schematic diagram of the tube screen structure in this invention;

[0051] Figure 3 This is a schematic diagram of the steam process in this invention;

[0052] Figure 4 This is a schematic diagram of the pressure difference at different positions in the left and right directions of the header in this invention;

[0053] Figure 5 This is a distribution diagram of the heat absorption non-uniformity coefficient in the width direction of the boiler in this invention;

[0054] Figure 6 This is a schematic diagram of the electronic device in this invention. Detailed Implementation

[0055] The technical solution of the present invention will be clearly and completely described below with reference to the embodiments. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0056] Example 1: Please refer to Figures 1-5 As shown, the technical solution provided by this invention is: a method for calculating the heat absorption non-uniformity coefficient of a boiler's high-temperature heating surface based on monitoring data. In this embodiment, taking the calculation of the heat absorption non-uniformity coefficient of the rear screen heating surface of a boiler of model HG-2038 / 26.15-YM3 in a power plant as an example, the specific method is as follows:

[0057] Step S1: Obtain real-time monitoring data of the boiler;

[0058] For details, please refer to Figure 2 and Figure 3 As shown, the boiler consists of a distribution header, a collection header, tube panels, tube coils, and a roof. Real-time monitoring data includes the total steam flow rate F, the steam pressure Pj at the inlet of the distribution header, the steam pressure Pc at the outlet of the collection header, the steam temperature tzj at the left inlet of the distribution header, the steam temperature tyj at the right inlet of the distribution header, and the wall temperature tb above the roof at the outlet of the same tube coil in different tube panels. i ;

[0059] In this embodiment, real-time monitoring data of the boiler can be obtained from the SIS system (Safety Instrumented System);

[0060] For example, the real-time monitoring data of the boiler collected at 19:51 on August 22, 2024 is shown in Table 1, as detailed below:

[0061]

[0062] It should be explained that, in this embodiment, a wall temperature measuring point is installed at the outlet section above the 6th ring of the ceiling from the outside in for each tube panel, where i is the sequence number of the corresponding tube panel from left to right.

[0063] Step S2: Calculate the steam parameters at the inlet and outlet of the same tube coil in different tube panels based on the real-time monitoring data of the boiler;

[0064] Among them, the steam parameters are the pressure and steam temperature at the inlet and the pressure and steam temperature at the outlet of the same tube coil in different tube panels;

[0065] In this embodiment, the calculation process in step S2 includes the following sub-steps:

[0066] Step S21: Obtain the total number N of the boiler heating surface width direction tube panels and the number of tube coils in a single tube panel; in this embodiment, the total number of boiler heating surface width direction tube panels is 35 panels, and a single tube panel has 18 tube coils.

[0067] Obtain the steam temperature tzj at the left inlet of the distribution header and the steam temperature tyj at the right inlet of the distribution header, and then use the formula tj i =tzj-(tzj-tyj)×(i-1) / (N-1)The steam temperature tj at the inlet of the same tube coil in different tube screens is calculated. i i is the leftmost tube screen number;

[0068] It should be explained that the steam in the distribution manifold enters from the desuperheater pipes on the left and right sides of the distribution manifold. Due to the adjustment of the desuperheating water, the steam temperature on both sides may be significantly different, resulting in a temperature gradient at different locations in the left and right directions.

[0069] Step S22: Obtain the wall temperature tb above the ceiling at the outlet of the same tube coil in different tube panels. i Through the formula tc i ≈tb i +1 is used to calculate the steam temperature tc at the outlet of the same tube coil in different tube screens. i ;

[0070] It should be explained that the wall temperature measuring point at the outlet is located inside the large enclosure above the ceiling. This pipe section is wrapped with thick insulation. According to heat transfer calculations, the steam temperature exceeds the wall temperature by less than 1°C.

[0071] Step S23: Obtain the steam pressure Pj at the inlet of the distribution header, and use the formula Pj i The pressure Pj at the inlet of the same tube coil in different tube panels is calculated by ≈Pj. i ;

[0072] Step S24: Obtain the steam pressure Pc at the outlet of the collection tank, and use the formula Pc iThe pressure Pc at the outlet of the same tube coil in different tube panels is calculated by ≈Pc. i For example, the calculation results of steam parameters at the inlet and outlet of the same tube coil in different tube panels are shown in Table 2, and Table 2 is as follows:

[0073]

[0074] It should be explained that, theoretically, the static pressure at each location within the header is different, but this difference is only on the order of one-thousandth of the average pressure, which is insufficient to cause a significant change in enthalpy and specific volume. Therefore, an approximation is made here.

[0075] Step S3: Calculate the mass flow rate of steam at the inlet of the distribution header and the outlet of the collection header.

[0076] The calculation process in step S3 is as follows:

[0077] Step S31: Collect the inner diameter Dj of the distribution header and the inner diameter Dc of the collection header to obtain the total steam flow rate F;

[0078] Step S32, using the formula fj=2F / (Π×Dj) 2 The mass flow rate fj of the steam at the inlet of the distribution header is calculated, where π is the coefficient of pi, π≈3.1415926;

[0079] Step S33, using the formula fc=2F / (Π×Dc) 2 The mass flow rate fc of the steam at the outlet of the collection header is calculated. For example, the calculation results of the mass flow rates of the steam at the inlet of the distribution header and the outlet of the collection header are shown in Table 3, which is as follows:

[0080]

[0081] Step S4: Calculate the enthalpy increase, enthalpy increase deviation, and average specific volume of steam in the same tube coil in different tube screens based on the steam enthalpy and steam specific volume at the inlet and outlet.

[0082] Specifically, the calculation process in step S4 is as follows:

[0083] Step S41, using formula Hj i =(a1×Pj 2 +a2×Pj+a3)ln(tj i )-(b1×Pj 2 The steam enthalpy Hj at the inlet of the same tube coil in different tube panels is calculated by adding b2×Pj+b3). i Similarly, through the formula Hc i =(a1×Pc 2 +a2×Pc+a3)ln(tci )-(b1×Pc 2 The steam enthalpy Hc at the outlet of the same tube coil in different tube panels is calculated by adding b2×Pc+b3). i Where a1, a2, a3, b1, b2, and b3 are all coefficients;

[0084] Step S42: Subtract the steam enthalpy at the inlet from the steam enthalpy at the outlet to obtain the enthalpy increase ΔH of the steam in the same tube coil in different tube panels. i ;

[0085] Step S43, using the formula The enthalpy increase deviation (HZP) of steam in the same tube coil in different tube panels was calculated. i ;

[0086] Step S44, using formula vj i =[e1×(Pj / 100) 2 +e2×(Pj / 100)+e3]ln(tj i )-[d1×(Pj / 100) 2 The specific steam volume vj at the inlet of the same tube coil in different tube screens is calculated by adding d2×(Pj / 100)+d3]. i Similarly, through the formula vc i =[e1×(Pc / 100) 2 +e2×(Pc / 100)+e3]ln(tc i )-[d1×(Pc / 100) 2 The specific steam volume vc at the outlet of the same tube coil in different tube panels is calculated by [+d2×(Pc / 100)+d3]. i Where d1, d2, d3, e1, e2, and e3 are all coefficients;

[0087] Step S45: Add the steam specific volume at the inlet and the steam specific volume at the outlet, and take the average value to obtain the average steam specific volume v in the same pipe ring in different pipe panels. i For example, the calculation results of the enthalpy increase, enthalpy increase deviation, and average specific volume of steam in the same tube coil in different tube panels are shown in Table 4, as detailed below:

[0088]

[0089] In this embodiment, the value of a1 is 0.5581, the value of a2 is 8.6019, the value of a3 is 1283.9, the value of b1 is 3.6182, the value of b2 is 62.84, the value of b3 is 4494, the value of e1 is 0.4927, the value of e2 is 0.2606, the value of e3 is 0.0512, the value of d1 is 2.4001, the value of d2 is 1.2537, and the value of d3 is 0.2573.

[0090] Step S5: Calculate the static pressure difference at the inlet and outlet of the same tube coil in different tube panels;

[0091] The calculation process for the pressure difference at the inlet and outlet of the same tube coil in different tube panels in step S5 includes the following sub-steps:

[0092] Step S51, firstly through the formula The average specific volume of steam at the inlet of the distribution header, vj, is calculated.

[0093] Similarly, through the formula The average specific volume vc of the steam at the outlet of the collection tank was calculated.

[0094] For step S52, please refer to... Figure 4 As shown, the steam velocity wj at the inlet of the distribution header is calculated using the formula wj=vj×fj. Similarly, the steam velocity wc at the outlet of the collection header is calculated using the formula wc=vc×fc.

[0095] Step S53, using the formula △Pjmax=kj×(wj) 2 The maximum static pressure difference ΔPjmax at the inlet of the distribution header is calculated using / 2vj).

[0096] Similarly, using the formula △Pcmax=kc×(wc) 2 The maximum static pressure difference ΔPcmax at the outlet of the collection tank is calculated using / 2vc.

[0097] Where kj is the pressure variation coefficient of the distribution header and kc is the pressure variation coefficient of the collection header; for example, the value of kj is 0.75 and the value of kc is 2.

[0098] Step S54, using the formula △P i =Pj-Pc+(△Pcmax-△Pjmax)×4(i-1)×(iN) / (N-1) 2 The static pressure difference ΔP at the inlet and outlet of the same tube coil in different tube panels was calculated. i ;

[0099] For example, the calculation results of the static pressure difference at the inlet and outlet of the same pipe coil in different pipe panels are shown in Table 5, and the details of Table 5 are as follows:

[0100]

[0101] It should be explained that the static pressure in the left and right directions of the header is distributed parabolically. Therefore, it is necessary to first calculate the maximum static pressure difference ΔPjmax at the inlet of the distribution header and the maximum static pressure difference ΔPcmax at the outlet of the collection header.

[0102] Step S6, calculate the steam flow deviation within the same pipe coil in different pipe panels, specifically:

[0103] Through formula Calculate the steam flow deviation (LLP) within the same tube coil in different tube panels. i ;

[0104] For example, the calculation results of the steam flow deviation within the same pipe ring in different pipe panels are shown in the 'Flow Deviation' column.

[0105] Step S7: Calculate the heat absorption non-uniformity coefficient of the same tube coil in different tube screens based on the steam flow deviation and the steam enthalpy increase deviation.

[0106] Among them, through formula Q i =LLP i ×HZP i The heat absorption non-uniformity coefficient Q of the same tube coil in different tube screens was calculated. i ;

[0107] For example, the calculation results of the heat absorption non-uniformity coefficient for the same tube coil in different tube panels are shown in the 'Heat Absorption Non-uniformity Coefficient' column, and the distribution diagram of the heat absorption non-uniformity coefficient in the width direction of the boiler is shown in the figure. Figure 5 As shown.

[0108] In this application, if a corresponding calculation formula appears, the above calculation formula is a dimensionless calculation. The weighting coefficient, proportional coefficient and other coefficients in the formula are set to quantify each parameter to obtain a result value. The size of the weighting coefficient and proportional coefficient is only required to not affect the proportional relationship between the parameter and the result value.

[0109] Example 2: Figure 6This is a schematic diagram of an electronic device, which may include a processor, a communication interface, a memory, and a communication bus. The processor, communication interface, and memory communicate with each other via the communication bus. The processor can call logical instructions in the memory to execute a method for calculating the heat absorption non-uniformity coefficient of a boiler's high-temperature heating surface based on monitoring data. This method includes: acquiring real-time monitoring data of the boiler; calculating steam parameters at the inlet and outlet of the same tube coil in different tube panels based on the real-time monitoring data; calculating the mass flow rate of steam at the inlet of the distribution header and the outlet of the collection header; calculating the enthalpy increase, enthalpy increase deviation, and average specific volume of steam in the same tube coil in different tube panels based on the steam enthalpy and specific volume at the inlet and outlet; calculating the static pressure difference at the inlet and outlet of the same tube coil in different tube panels; calculating the steam flow rate deviation in the same tube coil in different tube panels; and calculating the heat absorption non-uniformity coefficient of the same tube coil in different tube panels based on the steam flow rate deviation and the steam enthalpy increase deviation.

[0110] Furthermore, when the logical instructions in the aforementioned memory can be implemented as software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, or a portion of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0111] On the other hand, this application also provides a computer program product, which includes a computer program stored on a computer-readable storage medium. The computer program includes program instructions. When the program instructions are executed by a computer, the computer can execute the boiler high-temperature heating surface heat absorption non-uniformity coefficient calculation method based on monitoring data provided by the above methods. The method includes: acquiring real-time monitoring data of the boiler; calculating steam parameters at the inlet and outlet of the same tube coil in different tube panels based on the real-time monitoring data of the boiler; calculating the mass flow rate of steam at the inlet of the heating surface distribution header and the outlet of the collection header; calculating the enthalpy increase, enthalpy increase deviation, and average specific volume of steam in the same tube coil in different tube panels based on the steam enthalpy value and steam specific volume at the inlet and outlet; calculating the static pressure difference at the inlet and outlet of the same tube coil in different tube panels; calculating the steam flow deviation in the same tube coil in different tube panels; and calculating the heat absorption non-uniformity coefficient of the same tube coil in different tube panels based on the steam flow deviation and the steam enthalpy increase deviation.

[0112] In another aspect, this application also provides a computer-readable storage medium storing a computer program thereon, which, when executed by a processor, implements the above-mentioned methods for calculating the heat absorption non-uniformity coefficient of a boiler high-temperature heating surface based on monitoring data. The method includes: acquiring real-time monitoring data of the boiler; calculating steam parameters at the inlet and outlet of the same tube coil in different tube panels based on the real-time monitoring data of the boiler; calculating the mass flow rate of steam at the inlet of the distribution header and the outlet of the collection header; calculating the enthalpy increase, enthalpy increase deviation, and average specific volume of steam in the same tube coil in different tube panels based on the steam enthalpy value and steam specific volume at the inlet and outlet; calculating the static pressure difference at the inlet and outlet of the same tube coil in different tube panels; calculating the steam flow deviation in the same tube coil in different tube panels; and calculating the heat absorption non-uniformity coefficient of the same tube coil in different tube panels based on the steam flow deviation and the steam enthalpy increase deviation.

[0113] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs. Those skilled in the art can understand and implement this without any creative effort.

[0114] Through the above description of the embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus necessary general-purpose hardware platforms, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solutions, in essence or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute the methods described in the various embodiments or some parts of the embodiments.

[0115] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of this application, and are not intended to limit them. Although this application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of this application.

Claims

1. A method for calculating the heat absorption non-uniformity coefficient of a boiler's high-temperature heating surface based on monitoring data, characterized in that, include: Step S1: Obtain real-time monitoring data of the boiler; Step S2: Calculate the steam parameters at the inlet and outlet of the same tube coil in different tube panels based on the real-time monitoring data of the boiler; The calculation process in step S2 includes the following sub-steps: Step S21: Obtain the total number N of tube panels in the width direction of the boiler heating surface and the number of tube rings in a single tube panel; The steam temperature tzj at the left inlet of the distribution header and the steam temperature tyj at the right inlet of the distribution header are obtained, and the steam temperature tj at the inlet of the same tube ring in different tube panels is calculated by the formula tj = tzj - (tzj - tyj) × (i - 1) / (N - 1) i , wherein i is the left number of tube panel sequence i . Step S22: Obtain the wall temperature tb above the ceiling at the outlet of the same tube coil in different tube panels. i Through the formula tc i ≈tb i +1 is used to calculate the steam temperature tc at the outlet of the same tube coil in different tube screens. i ; Step S23: Obtain the steam pressure Pj at the inlet of the distribution header, and use the formula Pj i The pressure Pj at the inlet of the same tube coil in different tube panels is calculated by ≈Pj. i ; Step S24: Obtain the steam pressure Pc at the outlet of the collection tank, and use the formula Pc i The pressure Pc at the outlet of the same tube coil in different tube panels is calculated by ≈Pc. i ; Step S3: Calculate the mass flow rate of steam at the inlet of the distribution header and the outlet of the collection header. The calculation process in step S3 is as follows: Step S31: Collect the inner diameter Dj of the distribution header and the inner diameter Dc of the collection header to obtain the total steam flow rate F; Step S32, using the formula fj=2F / (Π×Dj) 2 The mass flow rate fj of the steam at the inlet of the distribution header is calculated, where π is the coefficient of pi. Step S33, using the formula fc=2F / (Π×Dc) 2 The mass flow rate fc of the steam at the outlet of the collection tank is calculated. Step S4: Calculate the enthalpy increase, enthalpy increase deviation, and average specific volume of steam in the same tube coil in different tube screens based on the steam enthalpy and steam specific volume at the inlet and outlet. The calculation process in step S4 is as follows: Step S41, using formula Hj i =(a1×Pj 2 +a2×Pj+a3)ln(tj i )-(b1×Pj 2 The steam enthalpy Hj at the inlet of the same tube coil in different tube panels is calculated by adding b2×Pj+b3). i ; Similarly, using the formula Hc i =(a1×Pc 2 +a2×Pc+a3)ln(tc i )-(b1×Pc 2 The steam enthalpy Hc at the outlet of the same tube coil in different tube panels is calculated by adding b2×Pc+b3). i Where a1, a2, a3, b1, b2, and b3 are all coefficients; Step S42: Subtract the steam enthalpy at the inlet from the steam enthalpy at the outlet to obtain the enthalpy increase ΔH of the steam in the same tube coil in different tube panels. i ; Step S43, using the formula The enthalpy increase deviation (HZP) of steam in the same tube coil in different tube panels was calculated. i ; Step S44, using formula vj i =[e1×(Pj / 100) 2 +e2×(Pj / 100)+e3]ln(tj i )-[d1×(Pj / 100) 2 The specific steam volume vj at the inlet of the same tube coil in different tube screens is calculated by adding d2×(Pj / 100)+d3]. i ; Similarly, using the formula vc i =[e1×(Pc / 100) 2 +e2×(Pc / 100)+e3]ln(tc i )-[d1×(Pc / 100) 2 The specific steam volume vc at the outlet of the same tube coil in different tube panels is calculated by [+d2×(Pc / 100)+d3]. i Where d1, d2, d3, e1, e2, and e3 are all coefficients; Step S45: Add the steam specific volume at the inlet and the steam specific volume at the outlet, and take the average value to obtain the average steam specific volume v in the same pipe ring in different pipe panels. i ; Step S5: Calculate the static pressure difference at the inlet and outlet of the same tube coil in different tube panels; The calculation process for calculating the pressure difference at the inlet and outlet of the same tube coil in different tube panels in step S5 includes the following sub-steps: Step S51, firstly through the formula The average specific volume of steam at the inlet of the distribution header, vj, is calculated. Similarly, through the formula The average specific volume vc of the steam at the outlet of the collection tank was calculated. Step S52: The steam velocity wj at the inlet of the distribution header is calculated using the formula wj=vj×fj. Similarly, the steam velocity wc at the outlet of the collection header is calculated using the formula wc=vc×fc. Step S53, using the formula △Pjmax=kj×(wj) 2 The maximum static pressure difference ΔPjmax at the inlet of the distribution header is calculated using / 2vj). Similarly, using the formula △Pcmax=kc×(wc) 2 The maximum static pressure difference ΔPcmax at the outlet of the collection tank is calculated using / 2vc. Where kj is the pressure variation coefficient of the distribution header and kc is the pressure variation coefficient of the collection header; Step S54, using the formula △P i =Pj-Pc+(△Pcmax-△Pjmax)×4(i-1)×(iN) / (N-1) 2 The static pressure difference ΔP at the inlet and outlet of the same tube coil in different tube panels was calculated. i ; Step S6: Calculate the steam flow deviation within the same pipe coil in different pipe panels; The calculation process in step S6 is as follows: Through formula Calculate the steam flow deviation (LLP) within the same tube coil in different tube panels. i ; Step S7: Calculate the heat absorption non-uniformity coefficient of the same tube coil in different tube screens based on the steam flow deviation and the steam enthalpy increase deviation. The calculation process in step S7 is as follows: Through formula Q i =LLP i ×HZP i The heat absorption non-uniformity coefficient Q of the same tube coil in different tube screens was calculated. i .

2. The method for calculating the heat absorption non-uniformity coefficient of a boiler high-temperature heating surface based on monitoring data according to claim 1, characterized in that, The real-time monitoring data in step S1 includes the total steam flow rate F, the steam pressure Pj at the inlet of the distribution header, the steam pressure Pc at the outlet of the collection header, the steam temperature tzj at the left inlet of the distribution header, the steam temperature tyj at the right inlet of the distribution header, and the wall temperature tb above the ceiling at the outlet of the same pipe coil in different pipe panels. i .

3. The method for calculating the heat absorption non-uniformity coefficient of a boiler high-temperature heating surface based on monitoring data according to claim 2, characterized in that, The steam parameters in step S2 are the pressure and steam temperature at the inlet and outlet of the same tube coil in different tube panels.

Images