WO2023065719A1 - Numerical simulation and hydrodynamic collaborative integration calculation method for boiler furnace - Google Patents

Numerical simulation and hydrodynamic collaborative integration calculation method for boiler furnace Download PDF

Info

Publication number
WO2023065719A1
WO2023065719A1 PCT/CN2022/103073 CN2022103073W WO2023065719A1 WO 2023065719 A1 WO2023065719 A1 WO 2023065719A1 CN 2022103073 W CN2022103073 W CN 2022103073W WO 2023065719 A1 WO2023065719 A1 WO 2023065719A1
Authority
WO
WIPO (PCT)
Prior art keywords
calculation
furnace
hydrodynamic
numerical simulation
boiler
Prior art date
Application number
PCT/CN2022/103073
Other languages
French (fr)
Chinese (zh)
Inventor
聂鑫
谢贝贝
刘辉
柳宏刚
Original Assignee
西安热工研究院有限公司
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 西安热工研究院有限公司 filed Critical 西安热工研究院有限公司
Publication of WO2023065719A1 publication Critical patent/WO2023065719A1/en

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Definitions

  • the application belongs to the boiler field of thermal power units, and in particular relates to a boiler furnace numerical simulation collaborative hydrodynamic integration calculation method.
  • this application proposes a calculation method using furnace numerical simulation in conjunction with hydrodynamic integration, which can obtain the combustion status in the furnace, the flow distribution of the water wall and the metal temperature distribution of the water wall tubes under various working conditions, which can be used for the operation and adjustment of power plant boilers. Provide guidance.
  • the purpose of this application is to provide a numerical simulation and hydrodynamic collaborative integrated calculation method for the boiler furnace.
  • This method aims at the problems of overheating of the water wall and bursting tubes in the current boiler, and obtains the combustion status in the furnace through numerical simulation calculation.
  • the simulated combustion heat load is hydrodynamically calculated, so as to obtain the flow distribution of the water-cooled wall tube, the metal wall temperature distribution, etc., analyze and propose design and operation countermeasures to solve problems such as overheating.
  • a numerical simulation and hydrodynamic collaborative integrated calculation method for a boiler furnace comprising the following steps:
  • hydrodynamic calculation modeling is carried out, including the determination of the working condition parameters of the water wall system, the division of the calculation circuit of the water wall system, the division of the circuit pipe section, etc., and the test is carried out simultaneously according to the burner structure data, fuel characteristics and test air volume.
  • the basic data is used for furnace numerical simulation modeling, and the entire furnace is meshed and calculated for boundary condition settings.
  • the numerical simulation of the whole furnace is carried out. Combustion numerical simulation is carried out according to the test conditions or optimized conditions, the combustion conditions of the whole furnace are observed and the heat flux distribution along the height direction of the furnace wall and the heat flux distribution of the horizontal section of the furnace are output.
  • the outlet temperature distribution of the water-cooled wall tube obtained from the actual boiler test is compared with the hydrodynamic calculation results. If the error between the calculated and measured temperature data is greater than 10%, the measured temperature data is used to correct the heat-absorbing deviation coefficient to obtain the corrected heat-absorbing deviation coefficient Then re-calculate the hydrodynamic force until the calculation error is less than 10%, then proceed to the next calculation.
  • the single-pipe flow distribution, pressure drop distribution, outlet temperature distribution of each flow circuit, enthalpy value distribution and furnace heating metal wall temperature distribution in the furnace are output (denoted as t metal wall surface temperature k+1 ). Compare “t metal wall temperature k+1 " and “t metal wall temperature k “, if the difference between "t metal wall temperature k+1 " and “t metal wall temperature k " is less than the set value ⁇ , then all calculations end, The furnace numerical simulation and hydrodynamic calculation results are obtained.
  • Step 1 of this application it specifically includes:
  • the purpose of the test is to determine the current operating conditions and characteristics of the boiler, and use this as a relative benchmark for subsequent adjustments and optimization. Observe the main performance parameters of the steam-water system, denitrification system, heating surface, blower fan, primary fan, induced draft fan, air preheater, feed water pump, condensate pump and control system; record the temperature measurement point data of the upper and lower furnace walls of the boiler to ensure The temperature is within the safe range.
  • the test records the main operating parameters of the boiler, measured the wind speed of the powder pipe and the distribution of pulverized coal, measured the oxygen content at the inlet of the SCR and the outlet of the air preheater, the flue temperature, the concentration of CO, the concentration of NO and atmospheric parameters, etc., and collected raw coal, fly ash and slag samples.
  • the specific test rules shall be carried out in accordance with the relevant test standards for power plant boiler tests.
  • step 2 of this application the hydrodynamic calculation modeling is performed according to relevant calculation standards.
  • Furnace numerical simulation modeling is carried out according to the calculation method of common combustion simulation commercial software.
  • Step 3 of the present application the numerical simulation of the whole furnace is carried out, and the grid division is performed to verify the grid quality and grid independence.
  • the calculation model selection, boundary conditions and other settings are determined in combination with the actual test parameters.
  • Step 4 of the present application the hydrodynamic calculation is performed according to relevant calculation standards.
  • Step 5 of the present application after the basic simulation test is completed, the furnace width and the heat absorption deviation coefficient in the depth direction are fitted according to the temperature distribution of the recorded water wall temperature measurement points.
  • the specific fitting method is as follows:
  • each measuring point represents the enthalpy value of the working fluid at the corresponding position in the water-cooled wall tube, but the pressure of the working medium in the water-cooled wall tube at each measuring point is unknown, and the pressure value at the position of the low-load steam-water separator is used as a reference.
  • the distribution of the enthalpy value of the working medium of the water wall tube at each measuring point is obtained, and then according to the enthalpy value at the entrance of the water wall, the difference in the increase of the enthalpy value at each measuring point is obtained , the ratio of the enthalpy increment at each measuring point to the average enthalpy increment is taken as the endothermic deviation coefficient;
  • ⁇ h i is the enthalpy increment of the water-cooled wall tube corresponding to the i-th measuring point
  • t i is the measured temperature of the water-cooled wall tube corresponding to the i-th measuring point
  • i is the number of the measuring point
  • p is the measured pressure at the outlet of the water-cooled wall
  • h The inlet is the enthalpy value of the working fluid at the inlet of the water wall, which is determined by the parameters of the working fluid at the outlet of the economizer;
  • i represents the number of the measuring point
  • N is the total number of measuring points
  • ⁇ h i is the enthalpy value increment of the i-th measuring point corresponding to the water-cooled wall tube
  • ⁇ i is the heat absorption deviation coefficient of the water-cooled wall tube at the position of the i-th measuring point, and i represents the number of the measuring point;
  • Step 5 of the present application the hydrodynamic calculation result is obtained, and the outlet temperature of each water-cooled wall tube is obtained according to the calculation, and the error between the measured temperature and the calculated value of the low-load water-cooled wall tube is compared. If the error is greater than 10%, then the calculated The flow deviation correction and heat absorption deviation correction heat load deviation are then substituted into and recalculated until the temperature error at the outlet of the water wall tube is less than 10%.
  • step six of the present application after the temperature error is less than 10%, according to the output calculation result, compare the difference between the metal wall temperature distribution used in step three and the calculated metal wall temperature distribution, if it is greater than the set error , then reassign the output metal wall temperature distribution of the hydrodynamic calculation to the numerical simulation calculation boundary conditions, and return to step 3 to start the iterative calculation. Until the error is less than the set value, output all numerical simulation calculation results and hydrodynamic calculation results.
  • the thermal power units in active service respond to the call of national policies, and the boilers of the units participate in deep peak shaving and deviate from the design operating conditions, and the mixed combustion of complex fuels causes problems such as overheating and tube explosion in boilers to varying degrees.
  • the key factors affecting boiler overheating and other problems can be obtained, and practical improvement measures can be put forward to ensure the safe operation of the boiler under deep peak load regulation and mixed combustion of complex fuels.
  • the application has the following advantages:
  • boundary conditions can be provided for the numerical simulation of the furnace wall through hydrodynamic calculation.
  • the wall surface of the water-cooled wall was often assumed to be a boundary condition of a certain temperature.
  • the temperature distribution of the water-cooled wall output by the hydrodynamic calculation results can be used as the wall boundary condition for numerical simulation calculations, which further improves the numerical simulation of furnace combustion. accuracy.
  • Fig. 1 is a schematic flowchart of a numerical simulation and hydrodynamic collaborative integrated calculation method for a boiler furnace of the present application.
  • Figure 2 is a schematic diagram of the numerical simulation modeling grid.
  • Fig. 3 is a schematic diagram of the water wall flow system for hydrodynamic calculation.
  • Figure 4 is a schematic diagram of the furnace temperature distribution.
  • Figure 5 is a schematic diagram of the furnace heat load distribution.
  • Figure 6 is a schematic diagram of calculation error analysis.
  • Fig. 7 is a schematic diagram of the distribution of metal wall temperature by hydrodynamic calculation.
  • the purpose of the test is to determine the current operating conditions and characteristics of the boiler, and use this as a relative benchmark for subsequent adjustments and optimization. Observe the main performance parameters of the steam-water system, denitrification system, heating surface, blower fan, primary fan, induced draft fan, air preheater, feed water pump, condensate pump and control system; record the temperature measurement point data of the upper and lower furnace walls of the boiler to ensure The temperature is within the safe range.
  • the test records the main operating parameters of the boiler, measured the wind speed of the powder pipe and the distribution of pulverized coal, measured the oxygen content at the inlet of the SCR and the outlet of the air preheater, the flue temperature, the concentration of CO, the concentration of NO and atmospheric parameters, etc., and collected raw coal, fly ash and slag samples.
  • the specific test rules shall be carried out in accordance with the relevant test standards for power plant boiler tests.
  • the application provides a boiler furnace numerical simulation and hydrodynamic collaborative integrated calculation method, as shown in Figure 1, which specifically includes the following steps:
  • hydrodynamic calculation modeling is carried out, including the determination of the working condition parameters of the water wall system, the division of the calculation circuit of the water wall system, the division of the circuit pipe section, etc., and the test is carried out simultaneously according to the burner structure data, fuel characteristics and test air volume.
  • the basic data is used for furnace numerical simulation modeling, and the entire furnace is meshed and calculated for boundary condition settings.
  • the numerical simulation of the whole furnace is carried out. Combustion numerical simulation is carried out according to the test conditions or optimized conditions, the combustion conditions of the whole furnace are observed and the heat flux distribution along the height direction of the furnace wall and the heat flux distribution of the horizontal section of the furnace are output.
  • the outlet temperature distribution of the water-cooled wall tube obtained from the actual boiler test is compared with the hydrodynamic calculation results. If the error between the calculated and measured temperature data is greater than 10%, the measured temperature data is used to correct the heat-absorbing deviation coefficient to obtain the corrected heat-absorbing deviation coefficient Then re-calculate the hydrodynamic force until the calculation error is less than 10%, then proceed to the next calculation.
  • the single-pipe flow distribution, pressure drop distribution, outlet temperature distribution of each flow circuit, enthalpy value distribution and furnace heating metal wall temperature distribution in the furnace are output (denoted as t metal wall surface temperature k+1 ). Compare “t metal wall temperature k+1 " and “t metal wall temperature k “, if the difference between "t metal wall temperature k+1 " and “t metal wall temperature k " is less than the set value ⁇ , then all calculations end, The furnace numerical simulation and hydrodynamic calculation results are obtained.
  • Step 1 of this application it specifically includes:
  • the purpose of the test is to determine the current operating conditions and characteristics of the boiler, and use this as a relative benchmark for subsequent adjustments and optimization. Observe the main performance parameters of the steam-water system, denitrification system, heating surface, blower fan, primary fan, induced draft fan, air preheater, feed water pump, condensate pump and control system; record the temperature measurement point data of the upper and lower furnace walls of the boiler to ensure The temperature is within the safe range.
  • the test records the main operating parameters of the boiler, measured the wind speed of the powder pipe and the distribution of pulverized coal, measured the oxygen content at the inlet of the SCR and the outlet of the air preheater, the flue temperature, the concentration of CO, the concentration of NO and atmospheric parameters, etc., and collected raw coal, fly ash and slag samples.
  • the specific test rules shall be carried out in accordance with the relevant test standards for power plant boiler tests.
  • step 2 of this application the hydrodynamic calculation modeling is performed according to relevant calculation standards.
  • Furnace numerical simulation modeling is carried out according to the calculation method of common combustion simulation commercial software.
  • Step 3 of the present application the numerical simulation of the whole furnace is carried out, and the grid division is performed to verify the grid quality and grid independence.
  • the calculation model selection, boundary conditions and other settings are determined in combination with the actual test parameters.
  • Step 4 of the present application the hydrodynamic calculation is performed according to relevant calculation standards.
  • Step 5 of the present application after the basic simulation test is completed, the furnace width and the heat absorption deviation coefficient in the depth direction are fitted according to the temperature distribution of the recorded water wall temperature measurement points.
  • the specific fitting method is as follows:
  • each measuring point represents the enthalpy value of the working fluid at the corresponding position in the water-cooled wall tube, but the pressure of the working medium in the water-cooled wall tube at each measuring point is unknown, and the pressure value at the position of the low-load steam-water separator is used as a reference.
  • the distribution of the enthalpy value of the working medium of the water wall tube at each measuring point is obtained, and then according to the enthalpy value at the entrance of the water wall, the difference in the increase of the enthalpy value at each measuring point is obtained , the ratio of the enthalpy increment at each measuring point to the average enthalpy increment is taken as the endothermic deviation coefficient;
  • ⁇ h i is the enthalpy increment of the i-th measuring point corresponding to the water-cooled wall tube
  • t i is the actual measured temperature of the water-cooled wall tube corresponding to the i-th measuring point
  • i is the number of the measuring point
  • p is the actual measured pressure at the outlet of the water-cooled wall
  • h The inlet is the enthalpy value of the working fluid at the inlet of the water wall, which is determined by the parameters of the working fluid at the outlet of the economizer;
  • i represents the number of the measuring point
  • N is the total number of measuring points
  • ⁇ h i is the enthalpy value increment of the i-th measuring point corresponding to the water-cooled wall tube
  • ⁇ i is the heat absorption deviation coefficient of the water-cooled wall tube at the position of the i-th measuring point, and i represents the number of the measuring point;
  • Step 5 of the present application the hydrodynamic calculation result is obtained, and the outlet temperature of each water-cooled wall tube is obtained according to the calculation, and the error between the measured temperature and the calculated value of the low-load water-cooled wall tube is compared. If the error is greater than 10%, then the calculated The flow deviation correction and heat absorption deviation correction heat load deviation are then substituted into and recalculated until the temperature error at the outlet of the water wall tube is less than 10%.
  • step six of the present application after the temperature error is less than 10%, according to the output calculation result, compare the difference between the metal wall temperature distribution used in step three and the calculated metal wall temperature distribution, if it is greater than the set error , then reassign the output metal wall temperature distribution of the hydrodynamic calculation to the numerical simulation calculation boundary conditions, and return to step 3 to start the iterative calculation. Until the error is less than the set value, output all numerical simulation calculation results and hydrodynamic calculation results.
  • Step 1 Boiler combustion basic test
  • test is carried out in accordance with the "Boiler Performance Test Regulations” (ASME PTC 4-1998) or “Power Plant Boiler Performance Test Regulations” (GB/T 10184-2015), and refer to the "Powdered Coal Boiler Combustion Adjustment Test Method”; the pulverizing system test is in accordance with “Performance Test Procedures for Power Plant Coal Mill and Pulverizing System” (DL/T467-2004).
  • a subcritical natural circulation boiler is tested for working conditions.
  • the working condition test is the test condition under the customary operation mode of the power plant. The purpose of this test is to understand the actual operation status of the boiler.
  • the test records the main operating parameters of the boiler, measures the oxygen content, smoke temperature, CO concentration, NO concentration and atmospheric parameters at the SCR inlet and the air preheater outlet, and collects raw coal, fly ash and slag samples, masters the basic fuel data, and calculates the thermal efficiency of the boiler and NOx production concentration.
  • the cold primary air volume at the inlet of the coal mill was calibrated, and the test was carried out under the condition that the coal mill did not feed coal.
  • the hot air door was closed, the cold air door was opened, the air volume of the coal mill was adjusted, the primary air velocity of the 8 pulverized coal pipes at the outlet of each mill was actually measured, and the actual air volume at the entrance of the coal mill was calculated.
  • the wind speed deviation of the 8 primary air pipes corresponding to each coal mill is adjusted to within ⁇ 5% by using the adjustment shrinkage cavity installed on each powder pipe.
  • Numerical simulation uses 3D modeling software to conduct full-scale geometric modeling of the prototype furnace, and then uses numerical simulation commercial software to mesh the computational domain.
  • the grid is encrypted to avoid the pseudo-diffusion phenomenon caused by the drastic change of the flow field, and the grid division in the upper part of the combustion area is gradually sparse to reduce the number of grids and improve the calculation speed.
  • the number of grids is 1.89 million, and the number of grids can meet the requirements of calculation accuracy through the verification of grid independence.
  • the furnace modeling and grid division are shown in Figure 2.
  • the furnace wall is regarded as no turbulence motion and no slip boundary conditions.
  • the inlets of the burner and the overfired air are set as the mass flow inlet boundary, the turbulence intensity is 10%, and the hydraulic diameter of each inlet is set according to the specific calculation results of the nozzle size.
  • the wind speed of each wind is given respectively, and the slip coefficient of pulverized coal particles is 0.8.
  • the outlet of the furnace is set as the pressure outlet boundary condition, the outlet pressure is given according to the actual operating parameters of the boiler, and the influence of gravity on the pulverized coal particles is considered.
  • the wind boundary conditions of the pulverized coal burner are as follows:
  • Step 5 the wall temperature distribution of the hydrodynamic calculation is shown, and it is judged whether the wall temperature of the water-cooled wall and the wall temperature boundary condition of the numerical simulation are within the error range, and if it is within the error range, output all The calculation results guide the further operation adjustment and design modification of the boiler. If it is not within the error range, repeat step 3 with the temperature distribution of the water wall calculated by hydrodynamics as the boundary condition.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Regulation And Control Of Combustion (AREA)

Abstract

The present application discloses a numerical simulation and hydrodynamic collaborative integration calculation method for a boiler furnace. According to the numerical simulation and hydrodynamic collaborative integration calculation method for the furnace, the combustion condition and the thermal load distribution condition in the furnace can be accurately obtained by means of furnace combustion numerical simulation calculation. Numerical simulation calculation may provide a reliable calculation basis for hydrodynamic calculation, and the accuracy of data such as hydrodynamic input heat flux density is ensured. Calculation boundary conditions can be provided by means of hydrodynamic calculation for the furnace wall surface subjected to numerical simulation. The water cooled wall temperature distribution output by a hydrodynamic calculation result can serve as a wall surface boundary condition of numerical simulation calculation, and the accuracy of furnace combustion numerical simulation is further improved. The calculation result can guide further operation adjustment and design transformation of the boiler, practical improvement measures are put forward, and safe operation of the boiler during depth peak load regulation and complex fuel blending combustion is ensured.

Description

一种锅炉炉膛数值模拟与水动力协同集成计算方法An Integrated Calculation Method of Boiler Furnace Numerical Simulation and Hydrodynamic Synergy
本申请要求在2021年10月21日提交中国专利局、申请号为202111227957.6、发明名称为“一种锅炉炉膛数值模拟与水动力协同集成计算方法”的中国专利申请的优先权，其全部内容通过引用结合在本申请中。This application claims the priority of the Chinese patent application submitted to the China Patent Office on October 21, 2021, with the application number 202111227957.6, and the title of the invention is "A Calculation Method for Boiler Furnace Numerical Simulation and Hydrodynamic Collaborative Integration", the entire content of which is passed References are incorporated in this application.
技术领域technical field
本申请属于火电机组的锅炉领域，具体涉及一种锅炉炉膛数值模拟协同水动力集成计算方法。The application belongs to the boiler field of thermal power units, and in particular relates to a boiler furnace numerical simulation collaborative hydrodynamic integration calculation method.
背景技术Background technique
在“双碳目标”下，国内燃煤机组频繁参与调峰，以提高新能源消纳能力，频繁调峰要求机组负荷响应速率快，并具备快速启停调峰能力，使得机组运行工况复杂，深度调峰期间锅炉运行偏离设计工况。加之燃料成本引起的煤种复杂多变，同时为了实现生物质减量化、资源化和无害化，更多电厂参与燃煤耦合生物质发电技改，这使电站锅炉炉内燃烧和热负荷情况更为复杂。Under the "dual carbon target", domestic coal-fired units frequently participate in peak regulation to improve the capacity of new energy consumption. Frequent peak regulation requires the units to respond quickly to load and have the ability to quickly start and stop peak regulation, which makes the operating conditions of the units complex. , the boiler operation deviates from the design condition during the deep peak shaving period. In addition, the coal types caused by fuel costs are complex and changeable. At the same time, in order to achieve biomass reduction, resource utilization and harmlessness, more power plants are involved in coal-fired coupled biomass power generation technology transformation, which makes the combustion and heat load of power plant boilers The situation is more complicated.
上述现象造成锅炉出现不同程度的燃烧恶化、结焦结渣、水冷壁工质流量偏差等问题，易导致受热面超温甚至爆管。针对上述问题，本申请提出了采用炉膛数值模拟协同水动力集成计算方法，可以得到多种工况条件下炉内燃烧状况和水冷壁流量分布及水冷壁管金属温度分布，为电站锅炉运行和调整提供指导意见。The above phenomena lead to different degrees of combustion deterioration of the boiler, coking and slagging, flow deviation of working fluid in the water wall, etc., which can easily lead to overheating of the heating surface and even tube explosion. In view of the above problems, this application proposes a calculation method using furnace numerical simulation in conjunction with hydrodynamic integration, which can obtain the combustion status in the furnace, the flow distribution of the water wall and the metal temperature distribution of the water wall tubes under various working conditions, which can be used for the operation and adjustment of power plant boilers. Provide guidance.
发明内容Contents of the invention
本申请的目的在于提供了一种锅炉炉膛数值模拟与水动力协同集成计算方法，该方法针对目前锅炉出现的水冷壁超温、爆管等问题，通过数值模拟计算得到炉内燃烧状况，利用数值模拟的燃烧热负荷进行水动力计算，从而得到水冷壁管的流量分布、金属壁温分布等，分析并提出解决超温等问题的设计及运行对策。The purpose of this application is to provide a numerical simulation and hydrodynamic collaborative integrated calculation method for the boiler furnace. This method aims at the problems of overheating of the water wall and bursting tubes in the current boiler, and obtains the combustion status in the furnace through numerical simulation calculation. The simulated combustion heat load is hydrodynamically calculated, so as to obtain the flow distribution of the water-cooled wall tube, the metal wall temperature distribution, etc., analyze and propose design and operation countermeasures to solve problems such as overheating.
为达到上述目的，本申请采用如下技术方案来实现的：In order to achieve the above object, the application adopts the following technical solutions to achieve:
一种锅炉炉膛数值模拟与水动力协同集成计算方法，包括以下步骤：A numerical simulation and hydrodynamic collaborative integrated calculation method for a boiler furnace, comprising the following steps:
步骤一、锅炉燃烧摸底基础试验 Step 1. Boiler combustion basic test
在燃煤电站锅炉进行燃烧摸底基础试验，在试验阶段主要测试锅炉运行总体情况，包括制粉***摸底、锅炉性能等试验内容，记录水冷壁出口温度数据，了解锅炉存在的主要问题，并检查各***运行状态是否正常，得到数值计算的边界条件。Carry out the basic test of combustion in coal-fired power plant boilers. In the test stage, the overall operation of the boiler is mainly tested, including the test content of the pulverization system and boiler performance, and the temperature data at the outlet of the water wall is recorded to understand the main problems of the boiler. Whether the system is running normally or not, the boundary conditions for numerical calculation are obtained.
步骤二、数值模拟建模和水动力计算建模Step 2. Numerical simulation modeling and hydrodynamic calculation modeling
根据锅炉燃烧摸底基础试验，进行水动力计算建模，包含水冷壁***工作条件参数确定、水冷壁***计算回路的划分、回路管段划分等，同步根据燃烧器结构数据、燃料特性和测试风量等试验基础数据进行炉膛数值模拟建模，对整个炉膛进行网格划分和计算边界条件设置。According to the basic test of boiler combustion, hydrodynamic calculation modeling is carried out, including the determination of the working condition parameters of the water wall system, the division of the calculation circuit of the water wall system, the division of the circuit pipe section, etc., and the test is carried out simultaneously according to the burner structure data, fuel characteristics and test air volume. The basic data is used for furnace numerical simulation modeling, and the entire furnace is meshed and calculated for boundary condition settings.
步骤三、全炉膛数值模拟Step 3. Numerical simulation of the whole furnace
假定锅炉炉膛水冷壁面温度分布(记为t 金属壁面温度k)，进行全炉膛数值模拟。按照试验工 况或优化工况进行燃烧数值模拟，观察全炉膛燃烧状况并输出沿炉膛墙面高度方向的热流密度分布和炉膛水平截面热流密度分布。 Assuming the temperature distribution of the water-cooled wall surface of the boiler furnace (denoted as t metal wall temperature k ), the numerical simulation of the whole furnace is carried out. Combustion numerical simulation is carried out according to the test conditions or optimized conditions, the combustion conditions of the whole furnace are observed and the heat flux distribution along the height direction of the furnace wall and the heat flux distribution of the horizontal section of the furnace are output.
步骤四、锅炉水动力计算 Step 4. Boiler hydrodynamic calculation
根据全炉膛数值模拟结果及输出的沿炉膛墙面高度方向的热流密度分布和炉膛水平截面热流密度分布，拟定锅炉水冷壁沿炉膛高度方向热负荷和吸热偏差系数，进行水动力计算。According to the numerical simulation results of the whole furnace and the output heat flux density distribution along the height direction of the furnace wall and the heat flux density distribution of the horizontal section of the furnace, the heat load and heat absorption deviation coefficient of the boiler water wall along the height direction of the furnace are drawn up, and the hydrodynamic calculation is carried out.
步骤五、吸热偏差系数修正 Step 5. Correction of endothermic deviation coefficient
由锅炉实测试验所得的水冷壁管出口温度分布和水动力计算结果进行比较，如果计算与实测温度数据误差大于10％，则用实测温度数据对吸热偏差系数进行修正，得到修正吸热偏差系数后重新进行水动力计算，直至计算误差小于10％后，进行下一步计算。The outlet temperature distribution of the water-cooled wall tube obtained from the actual boiler test is compared with the hydrodynamic calculation results. If the error between the calculated and measured temperature data is greater than 10%, the measured temperature data is used to correct the heat-absorbing deviation coefficient to obtain the corrected heat-absorbing deviation coefficient Then re-calculate the hydrodynamic force until the calculation error is less than 10%, then proceed to the next calculation.
步骤六、输出数值计算结果 Step 6. Output numerical calculation results
完成水动力计算后，则输出炉膛内受热管单管流量分布、压降分布、各个流动回路出口温度分布、焓值分布和炉膛受热金属壁温分布(记为t 金属壁面温度k+1)。比较“t 金属壁面温度k+1”和“t 金属 壁面温度k”，若“t 金属壁面温度k+1”和“t 金属壁面温度k”之差小于设定值ε，则所有计算结束，得到炉膛数值模拟和水动力计算结果。若“t 金属壁面温度k+1”和“t 金属壁面温度k”之差大于设定值ε，则将t 金属壁面温度k+1赋值给数值计算的壁温边界条件，重复步骤三、四、五、六，直至“t 金属壁面温度k+1”和“t 金属壁面温度k”之差小于设定值ε。 After the hydrodynamic calculation is completed, the single-pipe flow distribution, pressure drop distribution, outlet temperature distribution of each flow circuit, enthalpy value distribution and furnace heating metal wall temperature distribution in the furnace are output (denoted as t metal wall surface temperature k+1 ). Compare "t metal wall temperature k+1 " and "t metal wall temperature k ", if the difference between "t metal wall temperature k+1 " and "t metal wall temperature k " is less than the set value ε, then all calculations end, The furnace numerical simulation and hydrodynamic calculation results are obtained. If the difference between “t metal wall temperature k+1 ” and “t metal wall temperature k ” is greater than the set value ε, then assign the t metal wall temperature k+1 to the numerically calculated wall temperature boundary condition, and repeat steps 3 and 4 , 5, 6, until the difference between "t metal wall temperature k+1 " and "t metal wall temperature k " is less than the set value ε.
可选地，本申请步骤一中，具体包括：Optionally, in Step 1 of this application, it specifically includes:
试验目的在于测定锅炉目前运行状况及特性，并以此作为后续调整和优化改造的相对比较基准。观测汽水***、脱硝***、受热面、送风机、一次风机、引风机、空预器、给水泵、凝结水泵和控制***主要性能参数；记录锅炉上、下炉膛壁温测点数据，保证各受热面温度在安全范围内。试验记录锅炉主要运行参数，实测粉管风速和煤粉分配，实测SCR进口及空预器出口氧量、烟温、CO浓度、NO浓度及大气参数等，并采集原煤、飞灰、炉渣样品。具体试验细则根据电站锅炉试验相关试验标准执行。The purpose of the test is to determine the current operating conditions and characteristics of the boiler, and use this as a relative benchmark for subsequent adjustments and optimization. Observe the main performance parameters of the steam-water system, denitrification system, heating surface, blower fan, primary fan, induced draft fan, air preheater, feed water pump, condensate pump and control system; record the temperature measurement point data of the upper and lower furnace walls of the boiler to ensure The temperature is within the safe range. The test records the main operating parameters of the boiler, measured the wind speed of the powder pipe and the distribution of pulverized coal, measured the oxygen content at the inlet of the SCR and the outlet of the air preheater, the flue temperature, the concentration of CO, the concentration of NO and atmospheric parameters, etc., and collected raw coal, fly ash and slag samples. The specific test rules shall be carried out in accordance with the relevant test standards for power plant boiler tests.
可选地，本申请步骤二中，所述水动力计算建模，按照相关计算标准执行。炉膛数值模拟建模按照常用燃烧模拟商业软件计算方法进行。Optionally, in step 2 of this application, the hydrodynamic calculation modeling is performed according to relevant calculation standards. Furnace numerical simulation modeling is carried out according to the calculation method of common combustion simulation commercial software.
可选地，本申请步骤三中，所述全炉膛数值模拟，网格划分进行网格质量、网格无关性等验证。计算模型选取、边界条件等设定结合实际试验参数确定。Optionally, in Step 3 of the present application, the numerical simulation of the whole furnace is carried out, and the grid division is performed to verify the grid quality and grid independence. The calculation model selection, boundary conditions and other settings are determined in combination with the actual test parameters.
可选地，本申请步骤四中，所述水动力计算按照相关计算标准执行。Optionally, in Step 4 of the present application, the hydrodynamic calculation is performed according to relevant calculation standards.
可选地，本申请步骤五中，完成基础模拟试验后，根据记录水冷壁壁温测点温度大小分布拟合炉膛宽度、深度方向吸热偏差系数，具体拟合方法如下所示：Optionally, in Step 5 of the present application, after the basic simulation test is completed, the furnace width and the heat absorption deviation coefficient in the depth direction are fitted according to the temperature distribution of the recorded water wall temperature measurement points. The specific fitting method is as follows:
根据壁温测点的温度数值推算各测点代表水冷壁管内对应位置工质焓值大小，但各测点位置水冷壁管内工质压力未知，根据低负荷汽水分离器位置压力值作为参考，统一假设水冷壁测点位置处工质压力，得到各测点处水冷壁管工质焓值大小分布，再根据水冷壁入口位置焓值大小，得到每个测点位置焓值增大的差值大小，将每个测点位置的焓值增量与平均焓值增量之比作为吸热偏差系数；According to the temperature value of the wall temperature measuring point, each measuring point represents the enthalpy value of the working fluid at the corresponding position in the water-cooled wall tube, but the pressure of the working medium in the water-cooled wall tube at each measuring point is unknown, and the pressure value at the position of the low-load steam-water separator is used as a reference. Assuming the pressure of the working medium at the measuring point of the water wall, the distribution of the enthalpy value of the working medium of the water wall tube at each measuring point is obtained, and then according to the enthalpy value at the entrance of the water wall, the difference in the increase of the enthalpy value at each measuring point is obtained , the ratio of the enthalpy increment at each measuring point to the average enthalpy increment is taken as the endothermic deviation coefficient;
Δh i＝f(t i,p)-h 入口 Δh i =f(t i ,p)-h entry
式中Δh i为第i测点对应水冷壁管的焓值增量，t i为第i测点对应的水冷壁管实测温度，i表示测点编号，p为水冷壁出口位置实测压力，h 入口为水冷壁入口位置工质焓值，由省煤 器出口工质参数决定； In the formula, Δh i is the enthalpy increment of the water-cooled wall tube corresponding to the i-th measuring point, t i is the measured temperature of the water-cooled wall tube corresponding to the i-th measuring point, i is the number of the measuring point, p is the measured pressure at the outlet of the water-cooled wall, h The inlet is the enthalpy value of the working fluid at the inlet of the water wall, which is determined by the parameters of the working fluid at the outlet of the economizer;
Figure PCTCN2022103073-appb-000001
Figure PCTCN2022103073-appb-000001
式中
Figure PCTCN2022103073-appb-000002
为水冷壁管焓值增量的平均值，i表示测点编号，N为测点总数，Δh i为第i测点对应水冷壁管的焓值增量； In the formula
Figure PCTCN2022103073-appb-000002
is the average value of the enthalpy value increment of the water-cooled wall tube, i represents the number of the measuring point, N is the total number of measuring points, and Δh i is the enthalpy value increment of the i-th measuring point corresponding to the water-cooled wall tube;
Figure PCTCN2022103073-appb-000003
Figure PCTCN2022103073-appb-000003
式中η i为第i测点位置处水冷壁管吸热偏差系数，i表示测点编号； In the formula, ηi is the heat absorption deviation coefficient of the water-cooled wall tube at the position of the i-th measuring point, and i represents the number of the measuring point;
而吸热偏差η i又是与热负荷偏差η t、各管流量偏差η q以及结构偏差η s有关： And the heat absorption deviation η i is related to the heat load deviation η t , each pipe flow deviation η q and the structural deviation η s :
η i＝η tη q/η s η i =η t η q /η s
在初步计算时，假设各管流量偏差和结构偏差均相同，则吸热偏差与热负荷偏差近似相等。In the preliminary calculation, assuming that the flow deviation and structure deviation of each tube are the same, the heat absorption deviation and heat load deviation are approximately equal.
可选地，本申请步骤五中，得到水动力计算结果，根据计算得到各水冷壁管出口温度，比较低负荷水冷壁管实测温度和计算值的误差，如果误差大于10％，则用计算得到的流量偏差修正和吸热偏差修正热负荷偏差，再代入重新计算，直至水冷壁管出口温度误差小于10％。Optionally, in Step 5 of the present application, the hydrodynamic calculation result is obtained, and the outlet temperature of each water-cooled wall tube is obtained according to the calculation, and the error between the measured temperature and the calculated value of the low-load water-cooled wall tube is compared. If the error is greater than 10%, then the calculated The flow deviation correction and heat absorption deviation correction heat load deviation are then substituted into and recalculated until the temperature error at the outlet of the water wall tube is less than 10%.
可选地，本申请步骤六中，所述温度误差小于10％后，根据输出计算结果，比较步骤三中所用金属壁温分布和计算所得金属壁温分布的差值大小，若大于设定误差，则再重新将水动力计算输出金属壁温分布赋值给数值模拟计算边界条件，重新返回步骤三开始迭代计算。直至误差小于设定值，输出所有数值模拟计算结果和水动力计算结果。Optionally, in step six of the present application, after the temperature error is less than 10%, according to the output calculation result, compare the difference between the metal wall temperature distribution used in step three and the calculated metal wall temperature distribution, if it is greater than the set error , then reassign the output metal wall temperature distribution of the hydrodynamic calculation to the numerical simulation calculation boundary conditions, and return to step 3 to start the iterative calculation. Until the error is less than the set value, output all numerical simulation calculation results and hydrodynamic calculation results.
本申请至少具有如下有益的技术效果：The application has at least the following beneficial technical effects:
现役火电机组响应国家政策号召，机组锅炉参与深度调峰偏离设计工况运行，掺烧复杂燃料等造成锅炉出现不同程度的超温、爆管等问题。通过炉膛数值模拟与水动力协同集成计算方法，能够得到影响锅炉超温等问题的关键因素，并提出切合实际的改进措施，确保锅炉在深度调峰负荷和掺烧复杂燃料时安全运行。The thermal power units in active service respond to the call of national policies, and the boilers of the units participate in deep peak shaving and deviate from the design operating conditions, and the mixed combustion of complex fuels causes problems such as overheating and tube explosion in boilers to varying degrees. Through the numerical simulation of the furnace and the hydrodynamic collaborative integrated calculation method, the key factors affecting boiler overheating and other problems can be obtained, and practical improvement measures can be put forward to ensure the safe operation of the boiler under deep peak load regulation and mixed combustion of complex fuels.
具体而言，本申请具有以下优点：Specifically, the application has the following advantages:
(1)通过燃烧摸底试验，可以准确得到锅炉存在的问题所在，同时掌握锅炉实际运行基础数据，为数值计算提供可靠的依据。(1) Through the combustion test, the problems of the boiler can be accurately obtained, and at the same time, the basic data of the actual operation of the boiler can be mastered to provide a reliable basis for numerical calculation.
(2)通过炉膛燃烧数值模拟计算，可以准确得到炉内燃烧状况和热负荷分布情况。因炉内热负荷数据实测比较困难，而且准确性难以保证，数值模拟计算可以为水动力计算提供可靠的计算依据，确保水动力输入热流密度等数据的准确性。(2) Through the numerical simulation calculation of furnace combustion, the combustion status and heat load distribution in the furnace can be accurately obtained. Because the actual measurement of heat load data in the furnace is difficult and the accuracy is difficult to guarantee, the numerical simulation calculation can provide a reliable calculation basis for the hydrodynamic calculation and ensure the accuracy of the hydrodynamic input heat flux and other data.
(3)通过水动力计算可以为数值模拟的炉膛壁面提供计算边界条件。在以往的炉膛燃烧数值计算中，水冷壁墙面往往假定为某一定温度的边界条件，水动力计算结果输出的水冷壁壁温分布可以作为数值模拟计算的壁面边界条件，进一步提高炉膛燃烧数值模拟的准确性。(3) The calculation of boundary conditions can be provided for the numerical simulation of the furnace wall through hydrodynamic calculation. In previous numerical calculations of furnace combustion, the wall surface of the water-cooled wall was often assumed to be a boundary condition of a certain temperature. The temperature distribution of the water-cooled wall output by the hydrodynamic calculation results can be used as the wall boundary condition for numerical simulation calculations, which further improves the numerical simulation of furnace combustion. accuracy.
(4)通过协同计算研究，可针对燃煤电站深度调峰锅炉，能够校核低负荷炉内燃烧和水 动力安全特性，并提出锅炉水冷壁优化设计方案及运行对策；针对水冷壁壁温偏差较大的燃煤电站锅炉，能够校核计算不同工况燃烧和水动力特性，并提出锅炉水冷壁优化设计方案及运行对策。针对燃煤机组生物质耦合发电项目，可以研究分析掺烧生物质对锅炉燃烧和水动力影响，并提出设计运行对策。总体实现锅炉安全稳定运行的目标。(4) Through collaborative calculation research, it is possible to check the combustion and hydrodynamic safety characteristics of low-load furnaces for deep peak-shaving boilers in coal-fired power plants, and propose optimal design schemes and operational countermeasures for water-cooled walls of boilers; for temperature deviation of water-cooled walls Larger coal-fired power plant boilers can check and calculate combustion and hydrodynamic characteristics under different working conditions, and propose optimal design schemes and operation countermeasures for boiler water walls. For the biomass coupled power generation project of coal-fired units, it is possible to study and analyze the influence of biomass blending on boiler combustion and hydrodynamics, and propose design and operation countermeasures. Overall realization of the goal of safe and stable operation of the boiler.
附图说明Description of drawings
图1为本申请的一种锅炉炉膛数值模拟与水动力协同集成计算方法的流程示意图。Fig. 1 is a schematic flowchart of a numerical simulation and hydrodynamic collaborative integrated calculation method for a boiler furnace of the present application.
图2为数值模拟建模网格示意图。Figure 2 is a schematic diagram of the numerical simulation modeling grid.
图3为水动力计算水冷壁流动***示意图。Fig. 3 is a schematic diagram of the water wall flow system for hydrodynamic calculation.
图4为炉膛温度分布示意图。Figure 4 is a schematic diagram of the furnace temperature distribution.
图5为炉膛热负荷分布意图。Figure 5 is a schematic diagram of the furnace heat load distribution.
图6为计算误差分析示意图。Figure 6 is a schematic diagram of calculation error analysis.
图7为水动力计算金属壁温分布示意图。Fig. 7 is a schematic diagram of the distribution of metal wall temperature by hydrodynamic calculation.
具体实施方式Detailed ways
下面将参照附图更详细地描述本申请的示例性实施例。虽然附图中显示了本申请的示例性实施例，然而应当理解，可以以各种形式实现本申请而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更透彻地理解本申请，并且能够将本申请的范围完整的传达给本领域的技术人员。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。Exemplary embodiments of the present application will be described in more detail below with reference to the accompanying drawings. Although exemplary embodiments of the present application are shown in the drawings, it should be understood that the present application may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that the present application can be more thoroughly understood, and the scope of the present application can be fully conveyed to those skilled in the art. It should be noted that, in the case of no conflict, the embodiments in the present application and the features in the embodiments can be combined with each other. The present application will be described in detail below with reference to the accompanying drawings and embodiments.
如图1所示，本申请所述的一种锅炉炉膛数值模拟与水动力协同集成计算方法，需要在计算之前进行燃煤机组锅炉燃烧摸底基础试验，在试验阶段主要测试锅炉运行总体情况，得到数值计算的边界条件。As shown in Figure 1, the numerical simulation and hydrodynamic collaborative integrated calculation method of the boiler furnace described in this application requires a basic test of the combustion of the coal-fired unit boiler before the calculation. Boundary conditions for numerical calculations.
试验目的在于测定锅炉目前运行状况及特性，并以此作为后续调整和优化改造的相对比较基准。观测汽水***、脱硝***、受热面、送风机、一次风机、引风机、空预器、给水泵、凝结水泵和控制***主要性能参数；记录锅炉上、下炉膛壁温测点数据，保证各受热面温度在安全范围内。试验记录锅炉主要运行参数，实测粉管风速和煤粉分配，实测SCR进口及空预器出口氧量、烟温、CO浓度、NO浓度及大气参数等，并采集原煤、飞灰、炉渣样品。具体试验细则根据电站锅炉试验相关试验标准执行。The purpose of the test is to determine the current operating conditions and characteristics of the boiler, and use this as a relative benchmark for subsequent adjustments and optimization. Observe the main performance parameters of the steam-water system, denitrification system, heating surface, blower fan, primary fan, induced draft fan, air preheater, feed water pump, condensate pump and control system; record the temperature measurement point data of the upper and lower furnace walls of the boiler to ensure The temperature is within the safe range. The test records the main operating parameters of the boiler, measured the wind speed of the powder pipe and the distribution of pulverized coal, measured the oxygen content at the inlet of the SCR and the outlet of the air preheater, the flue temperature, the concentration of CO, the concentration of NO and atmospheric parameters, etc., and collected raw coal, fly ash and slag samples. The specific test rules shall be carried out in accordance with the relevant test standards for power plant boiler tests.
本申请提供的一种锅炉炉膛数值模拟与水动力协同集成计算方法，如图1具体包括以下步骤：The application provides a boiler furnace numerical simulation and hydrodynamic collaborative integrated calculation method, as shown in Figure 1, which specifically includes the following steps:
步骤一、锅炉燃烧摸底基础试验 Step 1. Boiler combustion basic test
在燃煤电站锅炉进行燃烧摸底基础试验，在试验阶段主要测试锅炉运行总体情况，包括制粉***摸底、锅炉性能等试验内容，记录水冷壁出口温度数据，了解锅炉存在的主要问题，并检查各***运行状态是否正常，得到数值计算的边界条件。Carry out the basic test of combustion in coal-fired power plant boilers. In the test stage, the overall operation of the boiler is mainly tested, including the test content of the pulverization system and boiler performance, and the temperature data at the outlet of the water wall is recorded to understand the main problems of the boiler. Whether the system is running normally or not, the boundary conditions for numerical calculation are obtained.
步骤二、数值模拟建模和水动力计算建模Step 2. Numerical simulation modeling and hydrodynamic calculation modeling
根据锅炉燃烧摸底基础试验，进行水动力计算建模，包含水冷壁***工作条件参数确定、水冷壁***计算回路的划分、回路管段划分等，同步根据燃烧器结构数据、燃料特性和测试风量等试验基础数据进行炉膛数值模拟建模，对整个炉膛进行网格划分和计算边界条件设置。According to the basic test of boiler combustion, hydrodynamic calculation modeling is carried out, including the determination of the working condition parameters of the water wall system, the division of the calculation circuit of the water wall system, the division of the circuit pipe section, etc., and the test is carried out simultaneously according to the burner structure data, fuel characteristics and test air volume. The basic data is used for furnace numerical simulation modeling, and the entire furnace is meshed and calculated for boundary condition settings.
步骤三、全炉膛数值模拟Step 3. Numerical simulation of the whole furnace
假定锅炉炉膛水冷壁面温度分布(记为t 金属壁面温度k)，进行全炉膛数值模拟。按照试验工况或优化工况进行燃烧数值模拟，观察全炉膛燃烧状况并输出沿炉膛墙面高度方向的热流密度分布和炉膛水平截面热流密度分布。 Assuming the temperature distribution of the water-cooled wall surface of the boiler furnace (denoted as t metal wall temperature k ), the numerical simulation of the whole furnace is carried out. Combustion numerical simulation is carried out according to the test conditions or optimized conditions, the combustion conditions of the whole furnace are observed and the heat flux distribution along the height direction of the furnace wall and the heat flux distribution of the horizontal section of the furnace are output.
步骤四、锅炉水动力计算 Step 4. Boiler hydrodynamic calculation
根据全炉膛数值模拟结果及输出的沿炉膛墙面高度方向的热流密度分布和炉膛水平截面热流密度分布，拟定锅炉水冷壁沿炉膛高度方向热负荷和吸热偏差系数，进行水动力计算。According to the numerical simulation results of the whole furnace and the output heat flux density distribution along the height direction of the furnace wall and the heat flux density distribution of the horizontal section of the furnace, the heat load and heat absorption deviation coefficient of the boiler water wall along the height direction of the furnace are drawn up, and the hydrodynamic calculation is carried out.
步骤五、吸热偏差系数修正 Step 5. Correction of endothermic deviation coefficient
由锅炉实测试验所得的水冷壁管出口温度分布和水动力计算结果进行比较，如果计算与实测温度数据误差大于10％，则用实测温度数据对吸热偏差系数进行修正，得到修正吸热偏差系数后重新进行水动力计算，直至计算误差小于10％后，进行下一步计算。The outlet temperature distribution of the water-cooled wall tube obtained from the actual boiler test is compared with the hydrodynamic calculation results. If the error between the calculated and measured temperature data is greater than 10%, the measured temperature data is used to correct the heat-absorbing deviation coefficient to obtain the corrected heat-absorbing deviation coefficient Then re-calculate the hydrodynamic force until the calculation error is less than 10%, then proceed to the next calculation.
步骤六、输出数值计算结果 Step 6. Output numerical calculation results
完成水动力计算后，则输出炉膛内受热管单管流量分布、压降分布、各个流动回路出口温度分布、焓值分布和炉膛受热金属壁温分布(记为t 金属壁面温度k+1)。比较“t 金属壁面温度k+1”和“t 金属 壁面温度k”，若“t 金属壁面温度k+1”和“t 金属壁面温度k”之差小于设定值ε，则所有计算结束，得到炉膛数值模拟和水动力计算结果。若“t 金属壁面温度k+1”和“t 金属壁面温度k”之差大于设定值ε，则将t 金属壁面温度k+1赋值给数值计算的壁温边界条件，重复步骤三、四、五、六，直至“t 金属壁面温度k+1”和“t 金属壁面温度k”之差小于设定值ε。 After the hydrodynamic calculation is completed, the single-pipe flow distribution, pressure drop distribution, outlet temperature distribution of each flow circuit, enthalpy value distribution and furnace heating metal wall temperature distribution in the furnace are output (denoted as t metal wall surface temperature k+1 ). Compare "t metal wall temperature k+1 " and "t metal wall temperature k ", if the difference between "t metal wall temperature k+1 " and "t metal wall temperature k " is less than the set value ε, then all calculations end, The furnace numerical simulation and hydrodynamic calculation results are obtained. If the difference between “t metal wall temperature k+1 ” and “t metal wall temperature k ” is greater than the set value ε, then assign the t metal wall temperature k+1 to the numerically calculated wall temperature boundary condition, and repeat steps 3 and 4 , 5, 6, until the difference between "t metal wall temperature k+1 " and "t metal wall temperature k " is less than the set value ε.
可选地，本申请步骤一中，具体包括：Optionally, in Step 1 of this application, it specifically includes:
试验目的在于测定锅炉目前运行状况及特性，并以此作为后续调整和优化改造的相对比较基准。观测汽水***、脱硝***、受热面、送风机、一次风机、引风机、空预器、给水泵、凝结水泵和控制***主要性能参数；记录锅炉上、下炉膛壁温测点数据，保证各受热面温度在安全范围内。试验记录锅炉主要运行参数，实测粉管风速和煤粉分配，实测SCR进口及空预器出口氧量、烟温、CO浓度、NO浓度及大气参数等，并采集原煤、飞灰、炉渣样品。具体试验细则根据电站锅炉试验相关试验标准执行。The purpose of the test is to determine the current operating conditions and characteristics of the boiler, and use this as a relative benchmark for subsequent adjustments and optimization. Observe the main performance parameters of the steam-water system, denitrification system, heating surface, blower fan, primary fan, induced draft fan, air preheater, feed water pump, condensate pump and control system; record the temperature measurement point data of the upper and lower furnace walls of the boiler to ensure The temperature is within the safe range. The test records the main operating parameters of the boiler, measured the wind speed of the powder pipe and the distribution of pulverized coal, measured the oxygen content at the inlet of the SCR and the outlet of the air preheater, the flue temperature, the concentration of CO, the concentration of NO and atmospheric parameters, etc., and collected raw coal, fly ash and slag samples. The specific test rules shall be carried out in accordance with the relevant test standards for power plant boiler tests.
可选地，本申请步骤二中，所述水动力计算建模，按照相关计算标准执行。炉膛数值模拟建模按照常用燃烧模拟商业软件计算方法进行。Optionally, in step 2 of this application, the hydrodynamic calculation modeling is performed according to relevant calculation standards. Furnace numerical simulation modeling is carried out according to the calculation method of common combustion simulation commercial software.
可选地，本申请步骤三中，所述全炉膛数值模拟，网格划分进行网格质量、网格无关性等验证。计算模型选取、边界条件等设定结合实际试验参数确定。Optionally, in Step 3 of the present application, the numerical simulation of the whole furnace is carried out, and the grid division is performed to verify the grid quality and grid independence. The calculation model selection, boundary conditions and other settings are determined in combination with the actual test parameters.
可选地，本申请步骤四中，所述水动力计算按照相关计算标准执行。Optionally, in Step 4 of the present application, the hydrodynamic calculation is performed according to relevant calculation standards.
可选地，本申请步骤五中，完成基础模拟试验后，根据记录水冷壁壁温测点温度大小分布拟合炉膛宽度、深度方向吸热偏差系数，具体拟合方法如下所示：Optionally, in Step 5 of the present application, after the basic simulation test is completed, the furnace width and the heat absorption deviation coefficient in the depth direction are fitted according to the temperature distribution of the recorded water wall temperature measurement points. The specific fitting method is as follows:
根据壁温测点的温度数值推算各测点代表水冷壁管内对应位置工质焓值大小，但各测点位置水冷壁管内工质压力未知，根据低负荷汽水分离器位置压力值作为参考，统一假设水冷壁测点位置处工质压力，得到各测点处水冷壁管工质焓值大小分布，再根据水冷壁入口位置焓值大小，得到每个测点位置焓值增大的差值大小，将每个测点位置的焓值增量与平均焓值增量之比作为吸热偏差系数；According to the temperature value of the wall temperature measuring point, each measuring point represents the enthalpy value of the working fluid at the corresponding position in the water-cooled wall tube, but the pressure of the working medium in the water-cooled wall tube at each measuring point is unknown, and the pressure value at the position of the low-load steam-water separator is used as a reference. Assuming the pressure of the working medium at the measuring point of the water wall, the distribution of the enthalpy value of the working medium of the water wall tube at each measuring point is obtained, and then according to the enthalpy value at the entrance of the water wall, the difference in the increase of the enthalpy value at each measuring point is obtained , the ratio of the enthalpy increment at each measuring point to the average enthalpy increment is taken as the endothermic deviation coefficient;
Δh i＝f(t i,p)-h 入口 Δh i =f(t i ,p)-h entry
式中Δh i为第i测点对应水冷壁管的焓值增量，t i为第i测点对应的水冷壁管实测温度， i表示测点编号，p为水冷壁出口位置实测压力，h 入口为水冷壁入口位置工质焓值，由省煤器出口工质参数决定； In the formula, Δh i is the enthalpy increment of the i-th measuring point corresponding to the water-cooled wall tube, t i is the actual measured temperature of the water-cooled wall tube corresponding to the i-th measuring point, i is the number of the measuring point, p is the actual measured pressure at the outlet of the water-cooled wall, h The inlet is the enthalpy value of the working fluid at the inlet of the water wall, which is determined by the parameters of the working fluid at the outlet of the economizer;
Figure PCTCN2022103073-appb-000004
Figure PCTCN2022103073-appb-000004
式中
Figure PCTCN2022103073-appb-000005
为水冷壁管焓值增量的平均值，i表示测点编号，N为测点总数，Δh i为第i测点对应水冷壁管的焓值增量； In the formula
Figure PCTCN2022103073-appb-000005
is the average value of the enthalpy value increment of the water-cooled wall tube, i represents the number of the measuring point, N is the total number of measuring points, and Δh i is the enthalpy value increment of the i-th measuring point corresponding to the water-cooled wall tube;
Figure PCTCN2022103073-appb-000006
Figure PCTCN2022103073-appb-000006
式中η i为第i测点位置处水冷壁管吸热偏差系数，i表示测点编号； In the formula, ηi is the heat absorption deviation coefficient of the water-cooled wall tube at the position of the i-th measuring point, and i represents the number of the measuring point;
而吸热偏差η i又是与热负荷偏差η t、各管流量偏差η q以及结构偏差η s有关： And the heat absorption deviation η i is related to the heat load deviation η t , each pipe flow deviation η q and the structural deviation η s :
η i＝η tη q/η s η i =η t η q /η s
在初步计算时，假设各管流量偏差和结构偏差均相同，则吸热偏差与热负荷偏差近似相等。In the preliminary calculation, assuming that the flow deviation and structure deviation of each tube are the same, the heat absorption deviation and heat load deviation are approximately equal.
可选地，本申请步骤五中，得到水动力计算结果，根据计算得到各水冷壁管出口温度，比较低负荷水冷壁管实测温度和计算值的误差，如果误差大于10％，则用计算得到的流量偏差修正和吸热偏差修正热负荷偏差，再代入重新计算，直至水冷壁管出口温度误差小于10％。Optionally, in Step 5 of the present application, the hydrodynamic calculation result is obtained, and the outlet temperature of each water-cooled wall tube is obtained according to the calculation, and the error between the measured temperature and the calculated value of the low-load water-cooled wall tube is compared. If the error is greater than 10%, then the calculated The flow deviation correction and heat absorption deviation correction heat load deviation are then substituted into and recalculated until the temperature error at the outlet of the water wall tube is less than 10%.
可选地，本申请步骤六中，所述温度误差小于10％后，根据输出计算结果，比较步骤三中所用金属壁温分布和计算所得金属壁温分布的差值大小，若大于设定误差，则再重新将水动力计算输出金属壁温分布赋值给数值模拟计算边界条件，重新返回步骤三开始迭代计算。直至误差小于设定值，输出所有数值模拟计算结果和水动力计算结果。Optionally, in step six of the present application, after the temperature error is less than 10%, according to the output calculation result, compare the difference between the metal wall temperature distribution used in step three and the calculated metal wall temperature distribution, if it is greater than the set error , then reassign the output metal wall temperature distribution of the hydrodynamic calculation to the numerical simulation calculation boundary conditions, and return to step 3 to start the iterative calculation. Until the error is less than the set value, output all numerical simulation calculation results and hydrodynamic calculation results.
实施例Example
本申请具体实施例如下：The specific embodiment of this application is as follows:
步骤一：锅炉燃烧摸底基础试验Step 1: Boiler combustion basic test
试验遵循《锅炉性能试验规程》(ASME PTC 4-1998)或《电站锅炉性能试验规程》(GB/T 10184-2015)进行，并参照《煤粉锅炉燃烧调整试验方法》；制粉***试验按照《电站磨煤机及制粉***性能试验规程》(DL/T467-2004)进行。The test is carried out in accordance with the "Boiler Performance Test Regulations" (ASME PTC 4-1998) or "Power Plant Boiler Performance Test Regulations" (GB/T 10184-2015), and refer to the "Powdered Coal Boiler Combustion Adjustment Test Method"; the pulverizing system test is in accordance with "Performance Test Procedures for Power Plant Coal Mill and Pulverizing System" (DL/T467-2004).
在某亚临界自然循环锅炉进行摸底工况试验，摸底工况试验为电厂习惯操作运行方式下的试验工况，该项试验的目的在于了解锅炉目前的实际运行状况。试验记录锅炉主要运行参数，实测SCR进口及空预器出口氧量、烟温、CO浓度、NO浓度及大气参数等，并采集原煤、飞灰、炉渣样品，掌握基础的燃料数据，计算锅炉热效率及NOx生成浓度。A subcritical natural circulation boiler is tested for working conditions. The working condition test is the test condition under the customary operation mode of the power plant. The purpose of this test is to understand the actual operation status of the boiler. The test records the main operating parameters of the boiler, measures the oxygen content, smoke temperature, CO concentration, NO concentration and atmospheric parameters at the SCR inlet and the air preheater outlet, and collects raw coal, fly ash and slag samples, masters the basic fuel data, and calculates the thermal efficiency of the boiler and NOx production concentration.
此外，为了掌握燃烧制粉***基本情况，对磨煤机进口冷态一次风量进行标定，该项试验在磨煤机不投煤的状况下进行。试验时关闭热风门，开启冷风门，调节磨煤机风量，对各台磨出口8根煤粉管的一次风速进行实际测量，并计算磨煤机进口实际通风量。在上述试验的同时，根据试验测得的一次风速偏差情况，利用安装在各粉管上的调节缩孔将各台磨煤机 对应的8根一次风管风速偏差调整至±5％以内。对磨煤机热态一次风量和风粉偏差进行测量，使用标准靠背管对热态带粉时的一次风速进行实际测量，并计算磨入口通风量，校核磨煤机热态运行时的表盘风量指示是否准确，并综合冷态试验结果进行最终修改。同时对各试验工况各根粉管的煤粉样进行称重，了解各根煤粉管道的粉量分配情况并计算热态一次风速偏差。In addition, in order to grasp the basic situation of the combustion pulverization system, the cold primary air volume at the inlet of the coal mill was calibrated, and the test was carried out under the condition that the coal mill did not feed coal. During the test, the hot air door was closed, the cold air door was opened, the air volume of the coal mill was adjusted, the primary air velocity of the 8 pulverized coal pipes at the outlet of each mill was actually measured, and the actual air volume at the entrance of the coal mill was calculated. At the same time as the above test, according to the deviation of the primary wind speed measured in the test, the wind speed deviation of the 8 primary air pipes corresponding to each coal mill is adjusted to within ±5% by using the adjustment shrinkage cavity installed on each powder pipe. Measure the primary air volume and air powder deviation of the coal mill in the hot state, use the standard back tube to actually measure the primary air speed in the hot state with powder, and calculate the ventilation volume at the mill inlet, and check the air volume of the dial when the coal mill is running in hot state Indicate whether it is accurate, and make final revisions based on the results of the cold test. At the same time, the pulverized coal samples of each pulverized coal pipe in each test working condition were weighed to understand the powder volume distribution of each pulverized coal pipe and calculate the deviation of the primary wind speed in the thermal state.
330MW负荷试验基础数据如下：The basic data of the 330MW load test are as follows:
Figure PCTCN2022103073-appb-000007
Figure PCTCN2022103073-appb-000007
Figure PCTCN2022103073-appb-000008
Figure PCTCN2022103073-appb-000008
步骤二、数值模拟建模和水动力计算建模Step 2. Numerical simulation modeling and hydrodynamic calculation modeling
数值模拟采用三维建模软件对原型炉膛进行全尺寸几何建模，然后通过数值模拟商业软件对计算域进行网格划分。在燃烧器喷口与炉膛交界的区域内，对网格进行了加密，以避免由于流场变化剧烈引起的伪扩散现象，在燃烧区域上部网格划分逐渐稀疏以减少网格数量，提高计算速度，网格数量为189万，网格数量通过网格无关性验证能够达到计算精度的要求，炉膛建模及网格划分如图2所示。Numerical simulation uses 3D modeling software to conduct full-scale geometric modeling of the prototype furnace, and then uses numerical simulation commercial software to mesh the computational domain. In the area at the junction of the burner nozzle and the furnace, the grid is encrypted to avoid the pseudo-diffusion phenomenon caused by the drastic change of the flow field, and the grid division in the upper part of the combustion area is gradually sparse to reduce the number of grids and improve the calculation speed. The number of grids is 1.89 million, and the number of grids can meet the requirements of calculation accuracy through the verification of grid independence. The furnace modeling and grid division are shown in Figure 2.
在数值计算边界条件设定中，将炉膛壁面视为无湍流运动且无滑移边界条件。燃烧器、燃尽风的入口均设置为质量流量入口边界，湍流强度为10％，各个入口的水力直径根据喷口尺寸具体计算结果设定。各次风风速分别给定，煤粉颗粒的滑移系数为0.8。炉膛出口设为压力出口边界条件，出口压力按照锅炉实际运行参数给定，考虑重力对煤粉颗粒的影响。In the setting of numerical calculation boundary conditions, the furnace wall is regarded as no turbulence motion and no slip boundary conditions. The inlets of the burner and the overfired air are set as the mass flow inlet boundary, the turbulence intensity is 10%, and the hydraulic diameter of each inlet is set according to the specific calculation results of the nozzle size. The wind speed of each wind is given respectively, and the slip coefficient of pulverized coal particles is 0.8. The outlet of the furnace is set as the pressure outlet boundary condition, the outlet pressure is given according to the actual operating parameters of the boiler, and the influence of gravity on the pulverized coal particles is considered.
煤粉燃烧器风边界条件如下：The wind boundary conditions of the pulverized coal burner are as follows:
Figure PCTCN2022103073-appb-000009
Figure PCTCN2022103073-appb-000009
在数值模拟计算得到的热流密度和热负荷不均匀系数分布的基础上进行水动力计算建模。如图3，采用将自然循环锅炉水冷壁等效为流动网络***的方法，将水冷壁划分为流量回路、压力节点、连接管等三类元件如图。根据质量守恒方程、动量守恒方程、能量守恒方程，建立了自然循环锅炉水冷壁流量分配计算模型。根据炉膛热负荷均流系数分布模型和等截面直肋导热控制方程，建立了内壁温度、中间点壁温、外壁温度、鳍根温度和鳍端温度沿炉高方向分布的计算模型。具体计算参考水动力计算步骤按照《JB/Z 201—83电站锅炉水动力计算方法》和专利《一种超超临界锅炉通用水动力计算方法》(CN106897547B,2019-04-12)。Based on the distribution of heat flux density and heat load non-uniformity coefficient calculated by numerical simulation, hydrodynamic calculation modeling is carried out. As shown in Figure 3, the method of equating the water wall of the natural circulation boiler to a flow network system is adopted, and the water wall is divided into three types of components, such as flow loop, pressure node, and connecting pipe, as shown in the figure. According to the mass conservation equation, momentum conservation equation, and energy conservation equation, a calculation model for the flow distribution of the water wall of the natural circulation boiler is established. According to the distribution model of furnace heat load flow sharing coefficient and the heat conduction control equation of equal cross-section straight ribs, the calculation model for the distribution of inner wall temperature, middle point wall temperature, outer wall temperature, fin root temperature and fin end temperature along the furnace height direction is established. The specific calculation refers to the hydrodynamic calculation steps in accordance with "JB/Z 201-83 Power Station Boiler Hydrodynamic Calculation Method" and the patent "A Universal Hydrodynamic Calculation Method for Ultra-Supercritical Boilers" (CN106897547B, 2019-04-12).
步骤三、全炉膛数值模拟Step 3. Numerical simulation of the whole furnace
根据输入数值模拟边界条件和计算模型，利用商业软件进行全炉膛数值，得到全炉膛温度分布等计算结果如图，导出水冷壁壁面热流密度分布和热负荷不均匀系数分布。According to the input numerical simulation boundary conditions and calculation model, commercial software is used to carry out the numerical value of the whole furnace, and the calculation results such as the temperature distribution of the whole furnace are obtained as shown in the figure, and the heat flux density distribution of the water-cooled wall surface and the distribution of the non-uniform coefficient of heat load are derived.
步骤四、锅炉水动力计算 Step 4. Boiler hydrodynamic calculation
根据试验阶段锅炉水冷壁工作条件数据，按照上文所述数值模拟结果拟定图5所示的热负荷和不均匀系数分布，作为水动力计算依据。水动力计算步骤按照《JB/Z 201—83电站锅炉水动力计算方法》和专利《一种超超临界锅炉通用水动力计算方法》(CN106897547B,2019-04-12)。According to the working condition data of the boiler water wall in the test stage, and according to the numerical simulation results mentioned above, the distribution of heat load and non-uniformity coefficient shown in Figure 5 is drawn up as the basis for hydrodynamic calculation. The hydrodynamic calculation steps follow the "JB/Z 201-83 Power Station Boiler Hydrodynamic Calculation Method" and the patent "A Universal Hydrodynamic Calculation Method for Ultra-Supercritical Boilers" (CN106897547B, 2019-04-12).
步骤五、吸热偏差系数修正 Step 5. Correction of endothermic deviation coefficient
初步计算试验阶段，比较水动力计算得到的回路出口温度与实测的温度误差是否小于10％，如图6误差分析所示。In the preliminary calculation test stage, compare whether the error between the circuit outlet temperature obtained by hydrodynamic calculation and the measured temperature is less than 10%, as shown in the error analysis in Figure 6.
步骤六、输出数值计算结果 Step 6. Output numerical calculation results
根据步骤五计算得到的水动力结果，如图7为水动力计算壁温分布，判断水冷壁壁温与数值模拟的壁温边界条件是否在误差范围内，，如果在误差范围内，则输出所有计算结果指导锅炉进一步运行调整和设计改造，如不在误差范围内，则以水动力计算水冷壁温度温度分布为边界条件重复步骤三。According to the hydrodynamic results calculated in Step 5, as shown in Figure 7, the wall temperature distribution of the hydrodynamic calculation is shown, and it is judged whether the wall temperature of the water-cooled wall and the wall temperature boundary condition of the numerical simulation are within the error range, and if it is within the error range, output all The calculation results guide the further operation adjustment and design modification of the boiler. If it is not within the error range, repeat step 3 with the temperature distribution of the water wall calculated by hydrodynamics as the boundary condition.
虽然，上文中已经用一般性说明及具体实施方案对本申请作了详尽的描述，但在本申请基础上，可以对之作一些修改或改进，这对本领域技术人员而言是显而易见的。因此，在不偏离本发明精神的基础上所做的这些修改或改进，均属于本申请要求保护的范围。Although the present application has been described in detail with general descriptions and specific implementations above, it is obvious to those skilled in the art that some modifications or improvements can be made on the basis of the present application. Therefore, the modifications or improvements made on the basis of not departing from the spirit of the present invention all belong to the protection scope of the present application.

Claims (6)

  1. 一种锅炉炉膛数值模拟与水动力协同集成计算方法，其特征在于，包括以下步骤：A boiler furnace numerical simulation and hydrodynamic collaborative integrated calculation method, characterized in that it includes the following steps:
    步骤一、锅炉燃烧摸底基础试验Step 1. Boiler combustion basic test
    在燃煤电站锅炉进行燃烧摸底基础试验，得到数值计算的边界条件；Carry out basic combustion experiments in coal-fired power plant boilers to obtain boundary conditions for numerical calculations;
    步骤二、数值模拟建模和水动力计算建模Step 2. Numerical simulation modeling and hydrodynamic calculation modeling
    根据锅炉燃烧摸底基础试验，进行水动力计算建模和炉膛数值模拟建模；According to the basic test of boiler combustion, carry out hydrodynamic calculation modeling and furnace numerical simulation modeling;
    步骤三、全炉膛数值模拟Step 3. Numerical simulation of the whole furnace
    假定锅炉炉膛水冷壁面温度分布(记为t 金属壁面温度k)，进行全炉膛数值模拟；输出沿炉膛墙面高度方向的热流密度分布和炉膛水平截面热流密度分布； Assuming the temperature distribution of the water-cooled wall surface of the boiler furnace (denoted as t metal wall surface temperature k ), the numerical simulation of the whole furnace is carried out; the heat flux distribution along the height direction of the furnace wall and the heat flux distribution of the horizontal section of the furnace are output;
    步骤四、锅炉水动力计算Step 4. Boiler hydrodynamic calculation
    根据全炉膛数值模拟结果及输出的沿炉膛墙面高度方向的热流密度分布和炉膛水平截面热流密度分布，拟定锅炉水冷壁沿炉膛高度方向热负荷和吸热偏差系数，进行水动力计算；According to the numerical simulation results of the whole furnace and the output heat flux density distribution along the height direction of the furnace wall and the heat flux density distribution of the horizontal section of the furnace, the heat load and heat absorption deviation coefficient of the boiler water wall along the height direction of the furnace are drawn up, and the hydrodynamic calculation is carried out;
    步骤五、吸热偏差系数修正Step 5. Correction of endothermic deviation coefficient
    由锅炉实测试验所得的水冷壁管出口温度分布和水动力计算结果进行比较，如果计算与实测温度数据误差大于10％，则用实测温度数据对吸热偏差系数进行修正，得到修正吸热偏差系数后重新进行水动力计算，直至计算误差小于10％后，进行下一步计算；The outlet temperature distribution of the water-cooled wall tube obtained from the actual boiler test is compared with the hydrodynamic calculation results. If the error between the calculated and measured temperature data is greater than 10%, the measured temperature data is used to correct the heat-absorbing deviation coefficient to obtain the corrected heat-absorbing deviation coefficient Then re-calculate the hydrodynamic force until the calculation error is less than 10%, then proceed to the next calculation;
    步骤六、输出数值计算结果Step 6. Output numerical calculation results
    完成水动力计算后，则输出炉膛内受热管单管流量分布、压降分布、各个流动回路出口温度分布、焓值分布和炉膛受热金属壁温分布(记为t 金属壁面温度k+1)；比较“t 金属壁面温度k+1”和“t 金属 壁面温度k”，若“t 金属壁面温度k+1”和“t 金属壁面温度k”之差小于设定值ε，则所有计算结束，得到炉膛数值模拟和水动力计算结果；若“t 金属壁面温度k+1”和“t 金属壁面温度k”之差大于设定值ε，则将t 金属壁面温度k+1赋值给数值计算的壁温边界条件，重复步骤三、四、五、六，直至“t 金属壁面温度k+1”和“t 金属壁面温度k”之差小于设定值ε。 After the hydrodynamic calculation is completed, the single-pipe flow distribution, pressure drop distribution, outlet temperature distribution of each flow circuit, enthalpy value distribution and furnace heating metal wall temperature distribution in the furnace are output (denoted as t metal wall surface temperature k+1 ); Compare "t metal wall temperature k+1 " and "t metal wall temperature k ", if the difference between "t metal wall temperature k+1 " and "t metal wall temperature k " is less than the set value ε, then all calculations end, Obtain the furnace numerical simulation and hydrodynamic calculation results; if the difference between "t metal wall temperature k+1 " and "t metal wall temperature k " is greater than the set value ε, then assign the t metal wall temperature k+1 to the numerical calculation Wall temperature boundary condition, repeat steps 3, 4, 5, and 6 until the difference between "t metal wall temperature k+1 " and "t metal wall temperature k " is less than the set value ε.
  2. 根据权利要求1所述的一种锅炉炉膛数值模拟与水动力协同集成计算方法，其特征在于，步骤一中，具体包括：According to claim 1, a boiler furnace numerical simulation and hydrodynamic collaborative integrated calculation method is characterized in that, in step 1, it specifically includes:
    试验目的在于测定锅炉目前运行状况及特性，并以此作为后续调整和优化改造的相对比较基准；观测汽水***、脱硝***、受热面、送风机、一次风机、引风机、空预器、给水泵、凝结水泵和控制***主要性能参数；记录锅炉上、下炉膛壁温测点数据，保证各受热面温度在安全范围内；试验记录锅炉主要运行参数，实测粉管风速和煤粉分配，实测SCR进口及空预器出口氧量、烟温、CO浓度、NO浓度及大气参数等，并采集原煤、飞灰、炉渣样品；具体试验细则根据电站锅炉试验相关试验标准执行。The purpose of the test is to determine the current operating conditions and characteristics of the boiler, and use it as a relative benchmark for subsequent adjustments and optimizations; observe the steam-water system, denitrification system, heating surface, blower, primary fan, induced draft fan, air preheater, feed water pump, The main performance parameters of the condensate pump and control system; record the temperature measurement points of the upper and lower furnace walls of the boiler to ensure that the temperature of each heating surface is within a safe range; test and record the main operating parameters of the boiler, the actual measurement of the wind speed of the powder pipe and the distribution of pulverized coal, and the actual measurement of the SCR inlet And air preheater outlet oxygen, flue temperature, CO concentration, NO concentration and atmospheric parameters, etc., and collect raw coal, fly ash, slag samples; specific test rules are implemented in accordance with the relevant test standards for power plant boiler tests.
  3. 根据权利要求1所述的一种锅炉炉膛数值模拟与水动力协同集成计算方法，其特征在于，步骤二中，所述水动力计算建模，按照相关计算标准执行；炉膛数值模拟建模按照常用燃烧模拟商业软件计算方法进行。According to claim 1, a boiler furnace numerical simulation and hydrodynamic collaborative integrated calculation method is characterized in that, in step 2, the hydrodynamic calculation modeling is performed according to relevant calculation standards; the furnace numerical simulation modeling is based on commonly used Combustion simulation is carried out by commercial software calculation methods.
  4. 根据权利要求1所述的一种锅炉炉膛数值模拟与水动力协同集成计算方法，其特征在于，步骤三中，所述全炉膛数值模拟，网格划分进行网格质量、网格无关性等验证；计算模型选取、边界条件等设定结合实际试验参数确定。According to claim 1, a boiler furnace numerical simulation and hydrodynamic collaborative integrated calculation method is characterized in that, in step 3, the numerical simulation of the whole furnace, the grid division is carried out to verify the grid quality, grid independence, etc. ; Calculation model selection, boundary conditions and other settings are determined in combination with actual test parameters.
  5. 根据权利要求4所述的一种锅炉炉膛数值模拟与水动力协同集成计算方法，其特征在于，步骤四、五中，所述水动力计算方法按照相关计算标准执行，得到初步计算结果，根据计算得到各水冷壁管出口温度，比较低负荷水冷壁管实测温度和计算值的误差，如果误差 大于10％，则用计算得到的流量偏差修正和吸热偏差修正热负荷偏差，再代入重新计算，直至水冷壁管出口温度误差小于10％。According to claim 4, a boiler furnace numerical simulation and hydrodynamic collaborative integrated calculation method is characterized in that, in steps 4 and 5, the hydrodynamic calculation method is executed according to relevant calculation standards to obtain preliminary calculation results, and according to the calculation Get the outlet temperature of each water-cooled wall tube, and compare the error between the measured temperature and the calculated value of the low-load water-cooled wall tube. If the error is greater than 10%, use the calculated flow deviation correction and heat absorption deviation to correct the heat load deviation, and then substitute it into the calculation again. Until the outlet temperature error of the water wall tube is less than 10%.
  6. 根据权利要求5所述的一种锅炉炉膛数值模拟与水动力协同集成计算方法，其特征在于，步骤六中，所述温度误差小于10％后，根据输出计算结果，比较步骤三中所用金属壁温分布和计算所得金属壁温分布的差值大小，若大于设定误差，则再重新将水动力计算输出金属壁温分布赋值给数值模拟计算边界条件，重新返回步骤三开始迭代计算；直至误差小于设定值，输出所有数值模拟计算结果和水动力计算结果。According to claim 5, a boiler furnace numerical simulation and hydrodynamic collaborative integrated calculation method is characterized in that in step six, after the temperature error is less than 10%, the metal wall used in step three is compared according to the output calculation result If the difference between the temperature distribution and the calculated metal wall temperature distribution is greater than the set error, then re-assign the metal wall temperature distribution output from the hydrodynamic calculation to the numerical simulation calculation boundary conditions, and return to step 3 to start iterative calculation; until the error If it is less than the set value, output all numerical simulation calculation results and hydrodynamic calculation results.
PCT/CN2022/103073 2021-10-21 2022-06-30 Numerical simulation and hydrodynamic collaborative integration calculation method for boiler furnace WO2023065719A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202111227957.6A CN113887155A (en) 2021-10-21 2021-10-21 Boiler furnace numerical simulation and hydrodynamic force collaborative integrated calculation method
CN202111227957.6 2021-10-21

Publications (1)

Publication Number Publication Date
WO2023065719A1 true WO2023065719A1 (en) 2023-04-27

Family

ID=79004202

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2022/103073 WO2023065719A1 (en) 2021-10-21 2022-06-30 Numerical simulation and hydrodynamic collaborative integration calculation method for boiler furnace

Country Status (2)

Country Link
CN (1) CN113887155A (en)
WO (1) WO2023065719A1 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116882173A (en) * 2023-07-11 2023-10-13 南方电网电力科技股份有限公司 Comprehensive evaluation method and device for boiler performance
CN116976162A (en) * 2023-06-16 2023-10-31 武汉光谷蓝焰新能源股份有限公司 Denitration design simulation method based on biomass boiler finite element analysis
CN116992189A (en) * 2023-07-24 2023-11-03 南方电网电力科技股份有限公司 Wall temperature analysis method and related device for boiler screen type superheater
CN117828797A (en) * 2024-02-21 2024-04-05 昆明理工大学 Design method and cooling system for water-cooling movable hood of side-blown converter
CN117852447A (en) * 2024-02-04 2024-04-09 南方电网电力科技股份有限公司 Universal boiler water wall rapid and accurate hydrodynamic force and wall temperature calculation method

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113887155A (en) * 2021-10-21 2022-01-04 西安热工研究院有限公司 Boiler furnace numerical simulation and hydrodynamic force collaborative integrated calculation method
CN115374636B (en) * 2022-08-23 2024-01-12 南方电网电力科技股份有限公司 Wall temperature prediction method for coupling numerical simulation and performance calculation of boiler
CN116796666B (en) * 2023-08-21 2023-11-07 中国航发上海商用航空发动机制造有限责任公司 Axial-flow compressor measuring point arrangement method
CN117823882B (en) * 2024-01-02 2024-05-17 中国电力工程顾问集团有限公司 Low-load hydrodynamic early warning method and device for coal-fired power plant boiler
CN117704349B (en) * 2024-01-02 2024-06-04 中国电力工程顾问集团有限公司 Low-load hydrodynamic self-adaptive regulation and control method and device for coal-fired power plant boiler

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102799775A (en) * 2012-07-06 2012-11-28 广东电网公司电力科学研究院 Method for calculating temperature of tube wall of platen superheater of ultra supercritical boiler by assistance of numerical simulation
CN106897547A (en) * 2017-01-12 2017-06-27 西安交通大学 A kind of general hydrodynamic calculation method of ultra-supercritical boiler
CN113887155A (en) * 2021-10-21 2022-01-04 西安热工研究院有限公司 Boiler furnace numerical simulation and hydrodynamic force collaborative integrated calculation method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102799775A (en) * 2012-07-06 2012-11-28 广东电网公司电力科学研究院 Method for calculating temperature of tube wall of platen superheater of ultra supercritical boiler by assistance of numerical simulation
CN106897547A (en) * 2017-01-12 2017-06-27 西安交通大学 A kind of general hydrodynamic calculation method of ultra-supercritical boiler
CN113887155A (en) * 2021-10-21 2022-01-04 西安热工研究院有限公司 Boiler furnace numerical simulation and hydrodynamic force collaborative integrated calculation method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
ZHANG, SHIHONG ET AL.: "Water Wall Temperature Deviation of An Ultra-supercritical 1,000 MW Unit Boiler: Analysis and Countermeasures", THERMAL POWER GENERATION, vol. 46, no. 11, 30 November 2017 (2017-11-30), XP009545735 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116976162A (en) * 2023-06-16 2023-10-31 武汉光谷蓝焰新能源股份有限公司 Denitration design simulation method based on biomass boiler finite element analysis
CN116976162B (en) * 2023-06-16 2024-01-30 武汉光谷蓝焰新能源股份有限公司 Denitration design simulation method based on biomass boiler finite element analysis
CN116882173A (en) * 2023-07-11 2023-10-13 南方电网电力科技股份有限公司 Comprehensive evaluation method and device for boiler performance
CN116992189A (en) * 2023-07-24 2023-11-03 南方电网电力科技股份有限公司 Wall temperature analysis method and related device for boiler screen type superheater
CN117852447A (en) * 2024-02-04 2024-04-09 南方电网电力科技股份有限公司 Universal boiler water wall rapid and accurate hydrodynamic force and wall temperature calculation method
CN117828797A (en) * 2024-02-21 2024-04-05 昆明理工大学 Design method and cooling system for water-cooling movable hood of side-blown converter

Also Published As

Publication number Publication date
CN113887155A (en) 2022-01-04

Similar Documents

Publication Publication Date Title
WO2023065719A1 (en) Numerical simulation and hydrodynamic collaborative integration calculation method for boiler furnace
CN105276563B (en) A kind of furnace outlet gas temperature flexible measurement method based on the real-time slagging scorification situation of burner hearth
CN108716664B (en) A kind of method and apparatus of on-line measurement burner hearth ash fouling coefficient
CN111339716B (en) Boiler high-temperature flue gas flow field online agent model construction method
CN103955599B (en) Acquisition method of combustion characteristic under different circumference air quantities after boiler improvement
CN115374636B (en) Wall temperature prediction method for coupling numerical simulation and performance calculation of boiler
Qiao et al. Combustion and NOx formation characteristics from a 330 MWe retrofitted anthracite-fired utility boiler with swirl burner under deeply-staged-combustion
CN116861737A (en) Method for relieving wall overtemperature of ultralow-load heating surface
CN111289278A (en) Method for evaluating effect of hedging coal-fired boiler after secondary air box flow field transformation
WO2022142321A1 (en) Air volume control method suitable for pulverized coal boiler
CN105509035B (en) A kind of method, apparatus and automatic control system of determining opposed firing intake
CN116776770B (en) CFD numerical simulation coupling BP neural network based method for relieving high-temperature corrosion of boiler
CN110864855A (en) Device and method for measuring air leakage rate of air preheater of power station boiler under medium and low loads
CN110298502A (en) Based on the boiler optimum oxygen calculation method that efficiency is optimal
Hernik et al. Numerical research of combustion with a minimum boiler load
Djayanti Energy Efficiency Improvement Strategies for Boilers: A Case Study in Pharmacy Industry
CN104806995B (en) Method for optimizing boiler operation under condition of varied coal quality
WO2022188375A1 (en) Pulverized coal boiler integrated comprehensive adjustment method based on water-cooled wall high-temperature corrosion prevention and control
CN104102842A (en) Nitrogen oxide emission predicting method and device based on economic index parameters
CN212319677U (en) Wide-load denitration and stable combustion system for auxiliary coal-fired unit
CN114659099A (en) Flue gas circulation system for improving combustion characteristics and thermal characteristics of pulverized coal fired boiler
CN112255978A (en) Intelligent operation and maintenance system and method for coal-fired boiler
Du et al. Numerical Study of Supercritical Opposed Wall-Fired Boiler Furnace Temperature and High-Temperature Heating Surface Stress under Variable Load Operation
FUa et al. Intelligent Prediction and Experimental Study on Deep Peak-Regulation Capability of Thermal Power Units
CN211502748U (en) Waste heat absorption flue gas recirculation system for improving efficiency of power plant unit

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 22882341

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE