Inﬂuence of Fuel Type and Water Content Variation on Pollutant Emission Characteristics of a Biomass Circulating Fluidized Bed Boiler

: In general, the biomass raw materials burned by biomass power plants generally have the characteristics of variable fuel types, high moisture content, and high volatile content. In this paper, a 130 t/h biomass circulating ﬂuidized bed (BCFB) model was established on the MWorks platform with Modelica language. The inﬂuence of biomass type changes on operation parameters, the corresponding steady-state characteristics, and the dynamic characteristics of the BCFB were carried out. The temperature corresponding to the combustion of pine was overall higher than that of the other fuels, and the ﬂue gas from the combustion of pine had the highest concentration of SO 2 , up to 520.49 mg/Nm 3 . The ﬂue gas from the combustion of pure cotton sticks had the highest concentration of NO, up to 254.34 mg/Nm 3 . The changes of fuel type and moisture content all have a great inﬂuence on the operation of BCFBs. The emission of pollutants was not only related to the element content of fuel, but also closely related to the furnace temperature. The fuel moisture content also indirectly affects the pollutant emission concentration and the steam-water system.


Introduction
Energy is an important material basis for the survival and development of human society. The extensive use of fossil energy has greatly promoted the development of society, but has also caused serious environmental pollution and greenhouse gas emissions [1]. Therefore, it has become the consensus of many countries to regard renewable energy and clean energy as an important part of the national energy structure. Among them, the utilization of biomass resources is considered to be one of the important ways to achieve carbon neutrality or negative CO 2 emissions, which is of great significance to ecological environmental protection, climate change mitigation and energy security [2][3][4]. At present, the direct combustion is the most economical and effective technology for the large-scale utilization of biomass energy [5].
The global energy demand is expected to grow by 53% from current levels by 2030, with 70% coming from China, India and other developing countries, while carbon dioxide emissions will increase by 55% [6]. In the long run, replacing fossil fuels with more biomass fuels can significantly reduce the greenhouse effect and reduce the impact of energy shortages [7]. Replacing coal power generation with biomass power generation can significantly reduce environmental pollution. Many countries have also formulated relevant goals and roadmaps to gradually realize the transformation from coal-fired power plants to biomass power generation [8]. China is also trying to reduce the proportion of coal-fired power generation, encouraging the use of biomass instead of coal-fired power generation. By the end of 2020, the total installed capacity of biomass power generation in China has reached 15 million kilowatts [9].
With its unique advantages of fuel flexibility, low NO x emissions and high sulfur capture efficiency, the circulating fluidized bed (CFB) technology is expected to continue to be one of the main furnace types for biomass direct-fired power generation units in the future [10]. However, the biomass burned in the furnace generally has the characteristics of variable fuel types, high moisture content, high volatile fraction, high alkali metal and chlorine content, and low calorific value [11]. In particular, the variation of fuel type and quality will cause drastic changes in combustion process variables (such as the furnace temperature, main steam pressure and temperature, flue gas composition, emission pollutants, etc.). This easily leads to combustion instability, heating surface corrosion and slagging, and even bed material agglomeration causing fluidization state failure. It brings serious challenges to the safe and stable operation of the generation units and low pollutant emission. A thorough understanding of the dynamic characteristics of CFB boilers is the basis of operation optimization and parameter adjustment. However, due to the restriction of the safe operation and test costs of power plant units, modeling simulation has become an important tool to study the characteristics of power plant systems.
In the modeling of biomass circulating fluidized beds, the current research mainly focuses on the system modeling of biomass in coal mixed fuel, and the research on pure biomass direct combustion circulating fluidized bed systems is less. Krzywansk et al. [12][13][14] mainly studied the mixed combustion of coal-fired and biomass, and studied the effect of different proportions of the material and pulverized coal on the combustion characteristics and the generation of pollutants. It was concluded that the SO2 concentration increases with increasing the proportion of coal combustion, and it was pointed out that the precipitation of volatile fraction from biomass combustion is not instantaneous and varies with the combustion temperature, as well as the height of the furnace. Gungor [15] established a 2-D model of a circulating fluidized bed biomass combustion chamber with three aspects of flow, heat transfer, and combustion. The axial and radial distributions of porosity, particle size distribution, pressure drop, gas emission, and temperature of gas phase and solid phase in the bottom and upper regions at different time intervals were simulated. The experimental results of small and industrial biomass combustors with different biomass fuels were compared.
Modelica is an object-oriented declarative physical modeling language based on C language. In addition, Modelica inherits the advantages of many previous domain modeling languages, and integrates the object-oriented mechanism of Java language and the array expression mechanism of Matlab. Based on the idea of non-causal modeling, the Modelica language uses mathematical equations (groups) and object-oriented structures to promote the reuse of model knowledge [16]. Therefore, the mathematical model of a biomass circulating fluidized bed (BCFB) boiler system was constructed by using Modelica language in this paper.
Compared with coal, the content of chlorine in biomass is high, and the chloride gas generated during the combustion process causes the high-temperature corrosion of various alloys. However, many previous modeling work did not consider the precipitation and release of chlorine during biomass combustion. Therefore, the overall model of a 130 t/h BCFB was established on the Mworks with Modelica language in this paper. The release of the chlorine element was considered. The effects of different types and quality of biomass fuel on the furnace temperature, steam and water system, and pollutant emissions in the flue gas were studied.

Combustion System Model
The furnace of a BCFB is a very complex system, including combustion, flow, heat transfer, and other physical and chemical processes in the furnace. Through reasonable simplification and assumptions, a mathematical model reflecting its dynamic characteristics was established, including a flow process model, combustion model, heat transfer model, desulfurization model, and an equal sub-model of mass and energy conservation. The fluid dynamic characteristics determine the fluidization velocity in the circulating fluidized bed and the limit of variable working conditions, the energy consumption of the auxiliary machine, the heat transfer in the bed, the temperature distribution, the bed memory material, and the wear of the heating surface. Therefore, the gas-solid, two-phase fluid dynamics is the basis of the heat and mass transfer in the CFB.
The average diameter of bubbles at the interface of the dense phase zone is treated by bubble bed [17], and the calculation formula is as follows: where, h den is the height of the dense phase region. N nos is the number of nozzles for distributors. A B is the cross-sectional area of the furnace in zone. U m f Critical fluidization velocity. U 0 is the apparent velocity of dense-phase bed gas. µ g is the dynamic viscosity of gas. The critical fluidization velocity (U mf ) is the minimum velocity that the fixed bed begins to transform to the fluidized bed. The calculation formula is [18]: where, d p , ϑ g , ρ p , ρ g are the bed material particle size, kinematic viscosity of fluidized gas, bed material particle density, and gas phase density, respectively. The axial void fraction distribution in the dilute phase region proposed by Kunii et al is [19].
The product of the superficial gas velocity and the decay index is a constant.
Hiller gives an empirical correlation for calculating the decay index, which is considered to be related only to the apparent gas velocity and the diameter of the particles.
The volume fraction of solid particles above transport disengaging height (TDH) is:

Combustion Model
Biomass combustion is divided into volatile combustion and coke combustion. The volatiles mainly include CO, CO 2 , H 2 , CH 4 , and C 2 H 6 . The combustion reaction equation and its kinetic equation are shown in Table 1 [20][21][22][23].

Ks
Where, N C is the number of particles, Φ is a chemical reaction mechanical factor, which determines the concentration balance between CO and CO 2 , and is related to coke particle diameter dc and coke surface temperature Tc. K d is diffusion reaction rate. K s is the surface chemical reaction rate of coke.
The conversion rate of fuel NO can be calculated by the following empirical formula [24]: where, Z = 100X VP − 30.637. X NO is the conversion of NO, %. X VP is the volatile content.
Assuming that all the sulfur elements contained in biomass are oxidized to SO 2 and released in the dense phase region. The reaction processes are as follows [25]: The SO 2 removal rate is as follows: where, r SO 2 is removal rate of SO 2 , mol/(m 3 •s). M Ca is the limestone accumulation amount, kg. ρ Ca is the density of limestone, kg/m 3 >. γ 1 is the activity coefficient of limestone, γ 1 = 0.035. C SO 2 is the molar of SO 2 in flue gas, mol/m 3 .
At present, a large number of studies [26][27][28][29] show that the chlorine mainly enters the gas phase in the form of HCl and KCl during the biomass combustion. At 700-900 • C, almost 100% chlorine enters the gas phase, while the potassium part enters the gas phase [30][31][32][33].
It is assumed that the chlorine is released in the dense phase region of the BCFB. The kinetic equation of solid chlorine conversion to gaseous chlorine is: In this formula, X is the precipitation amount of chlorine in fuel of unit mass. τ represents time. R is the ideal gas constant, 8.314 J/(mol•K). T B is the temperature of dense phase region, K.
The ratios of solid potassium to gaseous potassium in biomass fuels were 5% (700 • C), 10% (800 • C) and 30% (900 • C), respectively. The proportion of gaseous potassium at Energies 2021, 14, 5962 5 of 17 different temperatures was determined by the interpolation method. The yields of HCl and KCl in gas phase were calculated by the mass conservation method of each element.

Heat Transfer Model
The heat transfer coefficient between the fluidized bed and the heated surface is usually divided into three components.
where h gc is the convective heat transfer coefficient of gas. h pc is the convective heat transfer coefficient of particles, and h rad is the radiation heat transfer coefficient. Gas convective heat transfer coefficient is calculated by the following formula [34]: Prandt number Pr and Archimedes number Ar are calculated as follows: where d p is the equivalent average particle size of gas-solid two-phase flow material. C pg is the constant pressure specific heat of flue gas. Solid particle convective heat transfer below the formula calculation [35].
where, ε is the average porosity. C A is a constant, and the value is 2.85. C p is the specific heat of the material. λ g is the flue gas thermal conductivity. σ is the average free travel of gas molecules. C K is constant, 2.6 in dense phase and 2.0 in dilute phase. R is a gas constant. T is the gas temperature. M m is the average molecular weight. P is pressure. The radiative heat transfer between the bed and the heated surface occurs mainly between the high temperature particles and the heated surface. The radiative heat transfer in the CFB is not only related to the temperature of the heated surface, but also to the bed particle density, bed temperature, and other related and other physical quantities [36].
where, T B is bed temperature. T W is the wall temperature of heating surface. σ is the Boltzmann constant, and the value is 5.672 × 10 −8 W/m 2 K 4 . ε b is the equivalent effective emissivities of the bed, 0.8. ε w is the effective emissivities of the wall, 0.78.

Mass and Energy Conservation Models
The mass conservation of solid in each interval is as follows: where, M i is the accumulation of i interval solid. W in is the solid flow entering the i interval and W out,i is the solid flow leaving the i interval. W r,i is the solid amount consumed by the reaction in interval i, and W g,i is the amount of reaction production of solid in the i interval.
The following are the energy conservation equations for each region.
where, T k is the temperature in each zone, K. M p,k is the inventory of solids in each zone, kg. c p,k is specific heat at constant pressure in each zone, J/(kg•k). Q f uel,k is the energy that the fuel brings into each zone, J. Q air,k is the energy that air brings into each zone, J. Q in,k is the energy entering each zone, J. Q is the heat released by fuel combustion, J. Q out,k is the energy leaving each zone, J. Q exc,k is the amount of heat exchanged between each zone and the heat exchange equipment, J.

Economizer Model
The economizer, superheater, and reheater of the BCFB boiler are important parts of the steam-water system. The working fluids in the heating surface are a unidirectional medium. The heat transfer process and characteristics are similar, that is, no phase transition occurs in the heat transfer process [37]. In the modeling process, a simplified modeling method was adopted. The working fluid is regarded as a homogeneous medium, and the system modeling is carried out by using the piecewise lumped parameter method. The modeling process of the economizer is as follows: The mass conservation equation is as follows: where V sm represents the total volume of economizer, m 3 . dρ is the differential of the water supply density. τ represents time, s. D sm1 and D sm2 represent the import flow and export flow of economizer respectively, kg/s. The following is the energy conservation equation: where, H sm1 and H sm2 represent the enthalpy of economizer inlet and outlet, respectively, kJ/kg. Q sm is the heat absorption of water supply in economizer from the wall metal, kJ/s. The heat transfer formula is shown in (28) where, H sd is the set value of economizer outlet under 100% BMCR, kJ/kg.

Second-Order Model of Evaporation Zone
The pressure change model is a typical second-order model, that is, the model contains two state variables. When the second-order model is used for the evaporation zone, the pressure dynamic characteristics under water supply disturbance, fuel disturbance, and load disturbance can be simulated in a real way. When the evaporation zone is described by a more complex model, it can be refined on the basis of the above second-order model [38]. The mass and energy conservation equations for two working fluids in the evaporation zone are listed below.
The mass conservation equation is as follows: where, D sm represents the working fluid flow rate of the economizer into the evaporation zone, kg/s. D qb is the working fluid flow out of the evaporation zone, kg/s. V l and V v are the liquid volume in the evaporation zone and the gas volume in the evaporation zone, respectively, m 3 .
The following is the energy conservation equation: where, Q qb is the total heat transfer in the furnace to the evaporation zone, kJ/s. D l and D v are flow rates of water and steam, kg/s, respectively. H l and H v are the enthalpy of water and steam, kJ/kg, respectively. V qb is the total volume of evaporation zone, m 3 . P is drum pressure, MPa. M tm is effective metal mass, kg. C p is the metal specific heat, kJ/(kg K), the value is 0.45. T is the temperature of evaporation zone, K.
There are also the following volume relationships:

Calculation of Thermodynamic Properties of Water and Steam
The solution of the model also requires the calculation of the thermodynamic properties of water and steam. Modelica provides a set of standard working fluid library (Modelica. Media. Water) with a high calculation accuracy and fast calculation speed. Based on the IF97 formula (The IAPWS Industry Formulation 1997 for the thermodynamic properties of water and steam), the library provides a function library for solving the thermodynamic properties of water and steam [39,40]. The method of least squares fitting and other methods to solve the physical parameters may exceed the fitting range and lead to larger errors in the solution process. It is more efficient and accurate to call the IF97 formula in the standard library.

Simulation Model
The 130 t/h BCFB boiler (manufactured by Jinan Boiler Group Co., Ltd.No.8 Huanggang Road, Tianqiao District, Jinan City, Shandong Province, China) in a power plant was modeled. The combustion system is composed of a combustion chamber, a cyclone separator, and a U-shaped feed-back device. The structure and model structure of the circulating fluidized bed boiler are shown in Figure 1. The calorific value of the combustion system is the coupling point with the steamwater system. The steam-water system matching with the 30 MW high-temperature and high-pressure extraction condensing turbine generator is composed of an economizer, evaporation zone, and a superheater. The boiler structure and parameters are introduced in reference [41]. In this paper, the composite biomass fuel was used as a raw material, and the steady-state output of the model was compared with the actual measured value, as shown in Table 2.  The calorific value of the combustion system is the coupling point with the steamwater system. The steam-water system matching with the 30 MW high-temperature and high-pressure extraction condensing turbine generator is composed of an economizer, evaporation zone, and a superheater. The boiler structure and parameters are introduced in reference [41]. In this paper, the composite biomass fuel was used as a raw material, and the steady-state output of the model was compared with the actual measured value, as shown in Table 2. From Table 2, the relative error between the model output values and the actual measured values were within 7%. The relative errors of the economizer outlet temperature, drum pressure, superheater export temperature, outlet pressure of the superheater and main steam flow were within 0.4%. This means that the model established in this paper has a high reliability. The reason for the relatively large errors in the combustion-side parameters, such as bed temperature and oxygen content, is that the chemical reactions in the furnace are complex. In order to ensure the real-time nature of the model, a reasonable Energies 2021, 14, 5962 9 of 17 simplification of the model is made, and a control system will be added for error compensation in subsequent studies. Moreover, the reaction in the steam-water system is relatively simple, and the IF97 formula also makes the calculation more accurate.

Results and Discussion
The following five kinds of fuels were selected for the test. There were 40% bark + 30% sawdust + 20% split wood + 5% wheat straw + 5% corn straw (#1), pure broken wood (pine) (#2), pure cotton stalk (#3), pure corn stalk (#4), undried composite biomass (#5). Proximate and elemental analysis results of the five biomass fuels are shown in Table 3. It can be seen from Table 3 that the nitrogen and chlorine contents in the cotton stalk fuel are the highest. Corn straw contains the highest sulfur potassium and ash content. The difference in the proximate analysis between fuel (#1) and fuel (#5) is mainly due to the change of water. The effect of moisture will be reported in Section 3.2.3.

Steady-State Test
After the system is stable, the furnace temperature distribution is shown in Figure 2a. It can be seen that the temperature corresponding to the combustion of the biomass fuel pine was generally higher than that of other fuels, and the temperature difference in the suspension zone was more obvious. The reason may be that the combustion of fuel is mainly concentrated in the suspension zone. Figure 2b shows the output results of CO 2 and O 2 volume fraction in the flue gas when the four fuels are burned separately. Combined with Figure 2a, it can be seen that the CO 2 and O 2 volume fractions are closely related to the furnace chamber temperature. The (#2) fuel is burned with a high chamber temperature, high CO 2 concentration in the flue gas, and a low oxygen content. It can be seen from Figure 2c that the concentrations of NO, HCl, and SO 2 in the flue gas generated by the combustion of pure corn stalk (#4) were the highest. There were 199.47 mg/Nm 3 , 329.30 mg/Nm 3 and 412.42 mg/Nm 3 , respectively. The concentration of KCl in the flue gas generated by the combustion of the cotton stalk (#3) was the highest, up to 232.56 mg/Nm 3 . The SO 2 generated by the pine (#2) combustion was as high as 350.56 mg/Nm 3 . Combined with Table 3, the SO 2 in the flue gas generated by burning the fuel with more sulfur elements was high, and the transformation of the NO from the cotton stalk to corn stalk has the opposite change. This is because the generation of NO is not only related to the nitrogen content of fuel, but also related to temperature, fuel combustion state, and combustion atmosphere in the furnace. When burning pure corn stalk (#4), the pressure dropped corresponding to the whole furnace pressure drop was the largest, which was 5966 Pa (Figure 2d). Figure 2e shows the comparison of the main steam pressure and the main steam temperature for the four fuels. The main steam pressure and the main steam temperature were corresponding to the furnace temperature. When the pine (#2) was used as fuel, the main steam pressure and temperature were higher than those of the other three fuels. The change of fuel has a more obvious effect on the main steam temperature.
2d). Figure 2e shows the comparison of the main steam pressure and the main steam temperature for the four fuels. The main steam pressure and the main steam temperature were corresponding to the furnace temperature. When the pine (#2) was used as fuel, the main steam pressure and temperature were higher than those of the other three fuels. The change of fuel has a more obvious effect on the main steam temperature.

Dynamic Test
The BCFB boiler combustion system, which is a complex system with significant nonlinearity, has a strong parameter coupling and a large combustion response hysteresis

Dynamic Test
The BCFB boiler combustion system, which is a complex system with significant nonlinearity, has a strong parameter coupling and a large combustion response hysteresis [42,43]. There are many factors that affect the BCFB boiler, such as primary air, secondary air, feed rate, slag discharge, and circulating ratio. The changes of these operating parameters have complex and staggered effects on the BCFB. It is difficult to obtain the coupling relationship between the above factors and the bed temperature, bed pressure, main steam pressure, temperature, and pollutant concentration in the flue gas through simple experi- ments. Therefore, it is necessary to carry out the dynamic characteristics experiment based on a numerical model. Figure 3 shows the dynamic output of the model after reaching a steady state for the experiment using fuel (#1), where the increase in fuel volume is 5% and the increase in both primary and secondary air is 2.5%. [42,43]. There are many factors that affect the BCFB boiler, such as primary air, secondary air, feed rate, slag discharge, and circulating ratio. The changes of these operating parameters have complex and staggered effects on the BCFB. It is difficult to obtain the coupling relationship between the above factors and the bed temperature, bed pressure, main steam pressure, temperature, and pollutant concentration in the flue gas through simple experiments. Therefore, it is necessary to carry out the dynamic characteristics experiment based on a numerical model. Figure 3 shows the dynamic output of the model after reaching a steady state for the experiment using fuel (#1), where the increase in fuel volume is 5% and the increase in both primary and secondary air is 2.5%. As shown in Figure 3a, the furnace temperature increased with the increase of the fuel quantity step. In Figure 3b,c, with the increase of the furnace temperature, the flue gas temperature and metal wall temperature were also increased. Therefore, the main steam pressure and the main steam temperature would increase. When the valve opening was a constant, the evaporation in the evaporation zone would also increase with the increase of temperature. Therefore, the flow rate entering the superheater increased instantaneously. The outlet temperature of the superheater would decrease instantaneously, As shown in Figure 3a, the furnace temperature increased with the increase of the fuel quantity step. In Figure 3b,c, with the increase of the furnace temperature, the flue gas temperature and metal wall temperature were also increased. Therefore, the main steam pressure and the main steam temperature would increase. When the valve opening was a constant, the evaporation in the evaporation zone would also increase with the increase of temperature. Therefore, the flow rate entering the superheater increased instantaneously. The outlet temperature of the superheater would decrease instantaneously, and then the temperature gradually increased. The output results were basically consistent with the previous research results [44].

Dynamic Response of Fuel Variety Change
In order to explore the dynamic response of the boiler when the fuel type changes, it is assumed that the composite biomass fuel (#1) was used in the initial operation stage. After the system stabilizes for 1000 s, pure broken wood (pine) fuel (#2) was used due to the change of the biomass raw material source. Pure cotton rod fuel (#3) was used as the BCFB fuel after 6000 s, and pure corn stalk fuel (#4) was used after 11,000 s. The results are shown in Figure 4. and then the temperature gradually increased. The output results were basically consistent with the previous research results [44].

Dynamic Response of Fuel Variety Change
In order to explore the dynamic response of the boiler when the fuel type changes, it is assumed that the composite biomass fuel (#1) was used in the initial operation stage. After the system stabilizes for 1000 s, pure broken wood (pine) fuel (#2) was used due to the change of the biomass raw material source. Pure cotton rod fuel (#3) was used as the BCFB fuel after 6000 s, and pure corn stalk fuel (#4) was used after 11,000 s. The results are shown in Figure 4.  Figure 4a shows that the composite biomass (#1) was used in the initial operation stage, and pure broken wood (pine) (#2) was used after 1000 s, the furnace temperature decreased first and then increased due to the change of the biomass raw material source. After 6000 s, the whole temperature first increased and then decreased. After 11,000 s, the whole furnace temperature slowly decreased. Combined with Table 3, it is inferred that this phenomenon is related to the calorific value and moisture content of the fuel. When replacing biomass fuels, power plants should choose fuels with a similar calorific value and moisture content to avoid excessive fluctuations in operating conditions. It can be seen from Figure 4b-d that the change of fuel type has a complex impact on the composition of the flue gas. The change trend of the CO 2 concentration in the flue gas was consistent with that of the furnace temperature. The change of the oxygen content in the flue gas was opposite to the change trend of the CO 2 concentration, except in the instantaneous positions. Figure 4c shows that there is a certain correlation between the variation of the NO concentration and oxygen content in the flue gas, indicating that the generation of NO is closely related to oxygen. When pure pine fuel (#2) was used as the fuel, the NO concentration in the flue gas decreased rapidly first, then slightly increased, and then slowly decreased to equilibrium. The SO 2 concentration increased rapidly and then slowly decreased. When the fuel changes from pure broken pine (#2) to pure cotton stalk (#3), the NO concentration in the flue gas increased rapidly first, then decreased slowly and then increased slowly to reach equilibrium. The SO 2 concentration decreased rapidly and then decreased slowly to reach equilibrium. When pure corn stalk fuel (#4) was used to replace pure cotton stalk fuel (#3), the NO concentration in the flue gas decreased rapidly first and then increased slowly to reach equilibrium, and the SO 2 concentration increased rapidly and then increased slowly to reach equilibrium. It can be seen from Figure 4d that the variation trend of the HCl concentration in the flue gas is opposite to that of the KCl concentration, and the variation trend of the KCl concentration is generally consistent with that of the furnace temperature. In addition to the influence of the fuel element content, the combustion temperature of the furnace also has a great influence on the KCl concentration. Figure 4e shows that the saturated steam flow rate changes with the furnace temperature. With the increase of the furnace temperature, the heat absorption in the evaporation zone increased, and the evaporation capacity increased. Under the condition of a constant main steam valve and water supply, the water volume in the drum decreased. When the pure corn stalk (#4) was replaced, the furnace temperature began to decrease. The evaporation capacity decreased, and the saturated water volume began to rise. The change of the furnace temperature has a great influence on the drum water level, so the real-time monitoring and control of the water level is of great significance for the safe operation of the boiler.
The changes of the main steam pressure and temperature were consistent with the changes of the furnace temperature (Figure 4f). In the dynamic process, the maximum main steam pressure was 10.42 MPa, and the maximum main steam temperature was 862 K. With the increase of temperature, the outlet flow of the evaporation zone would increase. The fluctuation of the main steam temperature is the result of the combined effect of the flue gas temperature and superheater inlet flow.

Dynamic Response of Fuel Moisture Content Change
In order to explore the dynamic response of the boiler when the moisture content changes, the composite biomass fuel (#1) and the undried composite biomass fuel (#5) were used to carry out the experiment. It is assumed that the fuel (#1) was used in the initial stage, the fuel (#5) was replaced after 1000 s, and the fuel (#1) and fuel (#5) 1:1 mixed combustion was used after 5000 s (denoted as #6). The results are shown in Figure 5. It can be seen from Figure 5a that after the conversion from fuel (#1) to undried composite biomass (#5), the furnace temperature began to decrease. This effect on the temperature of the dense phase region was the most obvious, and the temperature decreased by 50 K. The effect on the temperature of the dilute phase region was small, and the temperature decreased by 33 K. After 6000 s, the mixed fuel (#6) was used instead, and the furnace It can be seen from Figure 5a that after the conversion from fuel (#1) to undried composite biomass (#5), the furnace temperature began to decrease. This effect on the temperature of the dense phase region was the most obvious, and the temperature decreased by 50 K. The effect on the temperature of the dilute phase region was small, and the temperature decreased by 33 K. After 6000 s, the mixed fuel (#6) was used instead, and the furnace temperature began to rise. The change of the fuel moisture content has a significant impact on the furnace temperature.
When the fuel moisture content increases, due to the decrease of the furnace temperature, the combustion becomes mild, the concentration of CO 2 decreased, and the concentration of O 2 increased (Figure 5b). With the decrease of the furnace temperature, the concentrations of the KCl and SO 2 decreased significantly, and the concentrations of the HCl and NO increased significantly (Figure 5c,d). The concentration of the pollution gas is closely related to the furnace temperature. The change of the fuel moisture content affects the emission of pollution gas by affecting the furnace temperature.
It can be seen from Figure 5e that with the decrease of the furnace temperature, the heat absorption in the evaporation zone decreased, and the evaporation decreased. Under the condition of a constant water supply and valve opening, the saturated water volume in the drum increased, and the water level raised. When the fuel (#6) was used as a combustion fuel after 6000 s, the furnace temperature raised and the saturated steam flow increased, but it was still less than the flow of fuel (#1). The volume growth of the saturated water in the drum slowed down.
As shown in Figure 5f, the main steam pressure varies with temperature, while the response of the main steam temperature is somewhat complex. When the fuel (#1) suddenly changed into the fuel (#5), the saturated steam flow into the superheater suddenly decreased and the main steam temperature increased. But then the heat absorption of the superheater decreased and the main steam temperature decreased. When the fuel was replaced (#6), the saturated steam flow at the inlet of the superheater increased and the main steam temperature decreased. Then with the increase of the heat absorption of the superheater, the main steam temperature increased and finally stabilized. Due to the saturated steam flow being less than the initial value, the final main steam temperature was higher than the temperature of fuel (#1), and the change of the fuel moisture content had an obvious influence on the steam-water system.

Conclusions
Based on Modelica language, a 130 t/h biomass circulating fluidized bed model was established on an MWorks platform. The relative error between the model output values with the actual measured values was within 7%. The effect on the boiler when there is a sudden change in fuel type and moisture content was studied. The concentration of SO 2 in the flue gas generated by pine combustion is the highest, up to 520.49 mg/Nm 3 , and the concentration of NO in the flue gas generated by pure cotton rod combustion is the highest, up to 254.34 mg/Nm 3 . The change trend of CO 2 in the flue gas is consistent with that of the temperature, while the change trend of the oxygen content in the flue gas is opposite to that of CO 2 in the flue gas. The change of NO in the flue gas is consistent with the change trend of oxygen content in the flue gas, indicating that the generation of NO is closely related to oxygen and furnace temperature. The main steam pressure is consistent with the main steam temperature and the furnace temperature, and the main steam temperature is greatly affected by the flow change. The change of fuel types and fuel moisture content have a great influence on the operation of the biomass circulating fluidized bed. The variation of the fuel moisture content has the greatest effect on the furnace chamber temperature in the dense phase zone.