阅读说明：本技术 基于单峰热解曲线的炭化可燃物热解动力学参数计算方法 (carbonized combustible pyrolysis kinetic parameter calculation method based on unimodal pyrolysis curve ) 是由 丁彦铭 张文龙 毛少华 黄必清 于 2019-07-17 设计创作，主要内容包括：本发明公开了一种基于单峰热解曲线的炭化可燃物热解动力学参数计算方法,包括：等转化率法FWO法求取可燃物单峰热解的活化能值E；假设多种反应机理,通过CR法求取活化能值E与FWO法得到的值对比,确定可燃物单峰热解的反应机理函数f(α),求取指前因子A；确定求解参数的取值范围；粒子群人工智能优化算法(PSO)对指前因子A、活化能值E以及炭生成率v进行寻优,并在最优解时输出计算结果。本发明降低了求可燃物单峰热解高阶偏微分方程的难度；计算精度好；优化效率高；使用范围广,可以同时应用到多组升温速率和多种可燃物的单峰热解求解中,有效预测实验结果。(the invention discloses a carbonized combustible pyrolysis kinetic parameter calculation method based on a unimodal pyrolysis curve, which comprises the following steps: obtaining the activation energy value E of unimodal pyrolysis of the combustible by an equal conversion method FWO method; assuming a plurality of reaction mechanisms, comparing the activation energy value E obtained by a CR method with a value obtained by a FWO method, determining a reaction mechanism function f (alpha) of unimodal pyrolysis of the combustible, and obtaining a pre-exponential factor A; determining the value range of the solving parameter; and optimizing the pre-pointing factor A, the activation energy value E and the carbon generation rate v by a particle swarm artificial intelligence optimization algorithm (PSO), and outputting a calculation result in the optimal solution. The method reduces the difficulty of solving the high-order partial differential equation of the unimodal pyrolysis of the combustible; the calculation precision is good; the optimization efficiency is high; the application range is wide, and the method can be simultaneously applied to solving unimodal pyrolysis of multiple groups of heating rates and multiple combustibles, so that the experimental result can be effectively predicted.)

1. A carbonized combustible pyrolysis kinetic parameter calculation method based on a unimodal pyrolysis curve is characterized by comprising the following steps:

s1, charred combustible pyrolysis single-step reaction formula based on unimodal pyrolysis curve: solid → v carbon + (1-v) volatile matter, v is carbon generation rate, the activation energy value E of pyrolysis of the carbonized combustible is obtained by adopting an equal conversion rate method FWO under different conversion rates alpha, and finally the average value of the activation energy values E under different conversion rates is obtained;

S2, under various reaction mechanisms, calculating an activation energy average value E of pyrolysis of the carbonized combustible materials under different mechanisms by using a CR method, comparing the activation energy average value E with a value obtained by using a FWO method, finding a reaction mechanism function f (alpha) corresponding to the CR method when the activation energy value calculated by using the FWO method is closest to the activation energy value calculated by using the CR method, and then calculating a pre-exponential factor A;

S3, floating the pre-finger factor A and the activation energy value E up and down by a certain percentage respectively to obtain the optimizing ranges of the pre-finger factor A and the activation energy value E respectively, and setting the optimizing range of the carbon generation rate v as (0, 1);

s4, optimizing the pre-indication factor A, the activation energy value E and the char generation rate v by adopting a particle swarm artificial intelligence optimization algorithm PSO program on the basis of the optimization range determined in the step S3, and outputting a conversion rate alpha, a conversion rate d alpha/dt, a mass loss m/m0 and a mass loss rate d (m/m0)/dt in the optimal solution, thereby completing the calculation of the pyrolysis kinetic parameters of the carbonized combustible; in the optimization process, the iterative update formula of the parameters is as follows:

α＝α+(dα/dt)×(t-t)，

(m/m)＝1-α×(1-v)，

(d(m/m)/dt)＝((m/m)-(m/m))/(t-t)，

dα/dt＝Af(α)exp(-E/RT)；

wherein T represents time, m0 represent instantaneous mass and initial mass of charred combustible pyrolysis, R is universal gas constant, T is absolute temperature, and i represents update iteration number.

2. the method for calculating pyrolysis kinetic parameters of carbonized combustible based on a unimodal pyrolysis curve according to claim 1, wherein the step S1 is specifically: single-step reaction formula of carbonized combustible unimodal pyrolysis: solid → v char + (1-v) volatile matter, v is char formation rate, under the condition that the mechanism of the charring combustible unimodal pyrolysis reaction is unknown, the temperature of fixed conversion rate corresponding to the experimental result is measured by the equal conversion rate method FWO method under different temperature rising rates, and the temperature is determined by the formula: the method comprises the following steps of (a) obtaining a linear slope-1.052 (E/R) of ln beta and 1/T under various conversion rates according to a formula, obtaining an activation energy value E of unimodal pyrolysis of the carbonized combustible, and finally obtaining an average value of the activation energy E under various conversion rates, wherein the ln beta is equal to the ln (AE/Rg (alpha)) -5.331-1.052(E/RT), R represents a universal gas constant, beta represents a heating rate, g (alpha) represents a differential reaction mechanism function, and T represents an absolute temperature.

3. the method for calculating pyrolysis kinetic parameters of charred combustible material based on a unimodal pyrolysis curve according to claim 1, wherein the step S2 is specifically as follows: respectively under a plurality of different reaction mechanisms assumed by the unimodal pyrolysis of the carbonized combustible, adopting a CR method to obtain an activation energy value E of the unimodal pyrolysis of the carbonized combustible, namely one reaction mechanism corresponds to one activation energy value, and obtaining the activation energy value E through a formula: ln (g (alpha)/T2) ═ ln (AR/beta E) - (E/RT), calculating slope of ln (g (alpha)/T2) and 1/T of a plurality of heating rates under each pyrolysis reaction mechanism to obtain activation energy value E of pyrolysis single peak of the carbonized combustible, calculating average value of a plurality of groups of heating rate activation energy values, comparing the average activation energy value calculated by CR method with the average activation energy value calculated by FWO method, finding out the mechanism corresponding to CR method when the activation energy value calculated by FWO method is closest to the activation energy value calculated by CR method, namely the reaction mechanism of pyrolysis of single peak of the carbonized combustible, determining reaction mechanism function f (alpha), after determining the mechanism, obtaining pre-index factor A through intercept of formula ln (AR/beta E).

4. The method for calculating pyrolysis kinetic parameters of carbonized combustible based on a unimodal pyrolysis curve as claimed in claim 1, wherein in step S3, the upper and lower floating by a certain percentage means that the upper and lower floating by 50% respectively.

5. the method for calculating pyrolysis kinetic parameters of carbonized combustible based on a unimodal pyrolysis curve according to claim 1, wherein the step S4 is specifically: when PSO is optimized, the fitness function value phi is the deviation degree of a predicted value and an experimental value, and the specific calculation formula is as follows:

wherein φ m, φ mlr, φ α and φ d α/dt represent the objective functions of mass loss, mass loss rate, conversion rate and conversion rate, respectively; n represents the number of experiments; n represents the number of experimental data points per experiment; CMLmod and CMLexp represent simulated and experimental values of accumulated mass loss; MLRmod and MLRexp represent a simulation value and an experimental value of the mass loss rate; α mod, α exp represent the simulated and experimental values of conversion; d α/dtmod, d α/dtexp represent the simulated and experimental values of the conversion rate; wCML, wMLR, w α, wd α/dt represent weighting factors for mass loss, rate of conversion.

6. the method for calculating pyrolysis kinetic parameters of charred combustible material based on unimodal pyrolysis curve according to claim 5, wherein each experimental value in step S4 means the calculation formula of the experimental value calculated by the FWO method in step S1 as follows:

CML＝m/m

MLR＝d(m/m)/dt＝

((m/m)-(m/m))/(t-t)

dα/dt＝Af(α)exp(-E/RT)；

Where m and mt represent the masses at the moment of reaction, m0 represents the initial mass of the sample, and m ∞ represents the final mass.

Technical Field

The invention relates to a carbonized combustible material pyrolysis kinetic parameter calculation method based on a single-peak pyrolysis curve, in particular to solving and optimization of carbonized combustible material single-peak pyrolysis kinetic parameters.

background

under the big background of the world energy crisis, the consumption of non-renewable energy needs to develop new energy technology urgently, and carbonized combustible, liquid and gas products can be obtained through the pyrolysis of the carbonized combustible, so that the renewable utilization of energy is realized. For example, pyrolysis can produce lignocellulosic biomass into fuel gas, bio-oil, bio-char, and the like. For example, the paper "Thermal degradation of char wood with Thermal analysis/Fourier transform induced analysis" with DOI of 10.1016/j.enconman.2016.05.007 indicates that charred combustible pyrolysis is an important part of energy conversion and has important research significance. Secondly, the life of people can not be separated from carbonized combustible substances such as synthetic materials, biomass, buildings and the like, the carbonized combustible substances bring convenience to people and certain dangerous hidden dangers, and fire disasters are the most common accidents. Pyrolysis, as the first step in the development of a fire, has led more and more researchers to study it. Finally, with the development of society, a large amount of carbonized waste combustible materials cannot be effectively and reasonably treated, and the current waste treatment mode comprises the following steps: landfill, incineration, pyrolysis is considered to be the most promising waste treatment method in view of land and environmental protection. At present, the solution of the pyrolysis dynamics of the carbonized combustible materials mostly exists in theoretical mathematical calculation, but the pyrolysis of most carbonized combustible materials needs a complex dynamics mechanism to simulate the pyrolysis process, and the difficulty in constructing a model, simulating and calculating is large, and the accuracy is low. Moreover, these mechanisms require a number of parameters that cannot be calculated by simple conventional pyrolysis methods.

The global optimization algorithm is an effective method for calculating and optimizing pyrolysis parameters, and PSO is an improved simulation method, and comprises the steps of firstly writing a PSO script main program, setting algorithm parameters and setting simulation output values. And then, the algorithm is operated, and the pre-exponential factor A, the activation energy E and the carbon generation rate v which are obtained on the basis of experimental data can be automatically searched in a set range, and the optimal solution and the deviation degree of each parameter are searched. The deviation degree provided by the algorithm represents the deviation degree of the prediction result and the experimental value, and the smaller the deviation degree is, the better the deviation degree is, and the best deviation degree is close to 0. For example, The paper "The accuracycacy and effectiveness of GA and PSO optimization schemes on estimation of reaction kinetics of biological pyrolysis" with DOI of 10.1016/j.energy.2019.04.030 constructs a calculation model on The basis of experimental values, calculates and optimizes pyrolysis kinetic parameters of carbonized combustible materials, and has short calculation time, high precision and highly consistent predicted values and experimental values. Therefore, the invention provides a chemical reaction kinetic parameter solving and optimizing method for the unimodal pyrolysis of the carbonized combustible, which effectively reduces the difficulty and complexity of solving a high-order partial differential equation of the unimodal pyrolysis of the carbonized combustible; the calculation precision is good; the optimization efficiency is high; the application range is wide, and the method can be simultaneously applied to unimodal pyrolysis solution of multiple groups of heating rates and multiple carbonized combustibles, so that the experimental result can be effectively predicted.

disclosure of Invention

The invention provides a carbonized combustible pyrolysis kinetic parameter calculation method based on a unimodal pyrolysis curve aiming at solving the problems existing in the existing carbonized combustible unimodal pyrolysis kinetic parameter, and the method reduces the difficulty and complexity of solving a carbonized combustible unimodal pyrolysis high-order partial differential equation; the calculation precision is good; the optimization efficiency is high; the application range is wide, and the method can be simultaneously applied to solving unimodal pyrolysis of multiple groups of heating rates and multiple carbonized combustibles, and the result is effectively calculated.

the invention solves the technical problem, and the adopted carbonized combustible material pyrolysis kinetic parameter calculation method based on the unimodal pyrolysis curve comprises the following steps:

s1, charred combustible pyrolysis single-step reaction formula based on unimodal pyrolysis curve: solid → v carbon + (1-v) volatile matter, v is carbon generation rate, the activation energy value E of pyrolysis of the carbonized combustible is obtained by adopting an equal conversion rate method FWO under different conversion rates alpha, and finally the average value of the activation energy values E under different conversion rates is obtained;

s2, under various reaction mechanisms, calculating an activation energy average value E of pyrolysis of the carbonized combustible materials under different mechanisms by adopting a CR method, comparing the activation energy average value E with a value obtained by an FWO method, finding out a reaction mechanism function f (alpha) of unimodal pyrolysis of the carbonized combustible materials when the activation energy value calculated by the FWO method is closest to the activation energy value calculated by the CR method, and then calculating an index factor A;

S3, floating the pre-pointing factor A and the closest activation energy value E up and down by a certain percentage respectively to obtain the optimizing ranges of the pre-pointing factor A and the closest activation energy value E respectively, and setting the optimizing range of the carbon generation rate v as (0, 1);

s4, optimizing the pre-pointing factor A, the closest activation energy value E and the char generation rate v by adopting a particle swarm artificial intelligence optimization algorithm PSO program on the basis of the optimization range determined in the step S3, and outputting a conversion rate alpha, a conversion rate d alpha/dt, a mass loss m/m0 and a mass loss rate d (m/m0)/dt in the optimal solution, thereby completing the calculation of the pyrolysis kinetic parameters of the carbonized combustible; in the optimization process, the iterative update formula of the parameters is as follows:

α＝α+(dα/dt)×(t-t)，

(m/m)＝1-α×(1-v)，

(d(m/m)/dt)＝((m/m)-(m/m))/(t-t)，

dα/dt＝Af(α)exp(-E/RT)；

Wherein T represents time, m0 represent instantaneous mass and initial mass of charred combustible pyrolysis, R is universal gas constant, T is absolute temperature, and i represents update iteration number.

Further, in the method for calculating pyrolysis kinetic parameters of carbonized combustible based on a unimodal pyrolysis curve of the invention, the step S1 is specifically: single-step reaction formula of carbonized combustible unimodal pyrolysis: solid → v char + (1-v) volatile matter, v is char formation rate, under the condition that the mechanism of the charring combustible unimodal pyrolysis reaction is unknown, the temperature of fixed conversion rate corresponding to the experimental result is measured by the equal conversion rate method FWO method under different temperature rising rates, and the temperature is determined by the formula: the method comprises the following steps of (a) obtaining a linear slope-1.052 (E/R) of ln beta and 1/T under various conversion rates according to a formula, obtaining an activation energy value E of unimodal pyrolysis of the carbonized combustible, and finally obtaining an average value of the activation energy E under various conversion rates, wherein the ln beta is equal to the ln (AE/Rg (alpha)) -5.331-1.052(E/RT), R represents a universal gas constant, beta represents a heating rate, g (alpha) represents a differential reaction mechanism function, and T represents an absolute temperature.

further, in the method for calculating pyrolysis kinetic parameters of carbonized combustible based on a unimodal pyrolysis curve according to the present invention, the step S2 specifically includes: respectively under a plurality of different reaction mechanisms of the unimodal pyrolysis of the carbonized combustible, adopting a CR method to obtain an activation energy value E of the unimodal pyrolysis of the carbonized combustible, namely one reaction mechanism corresponds to one activation energy value, and obtaining the activation energy value E through a formula: ln (g (alpha)/T2) ═ ln (AR/beta E) - (E/RT), calculating slopes of ln (g (alpha)/T2) and 1/T of a plurality of heating rates under each pyrolysis reaction mechanism to obtain an activation energy value E of a pyrolysis single peak of the carbonized combustible, calculating an average value of a plurality of groups of heating rate activation energies, comparing the average activation energy value calculated by the CR method with the average activation energy value calculated by the FWO method, finding a mechanism corresponding to the CR method when the activation energy value calculated by the FWO method is closest to the activation energy value calculated by the CR method, namely a reaction mechanism of pyrolysis of the single peak of the carbonized combustible, determining a reaction mechanism function f (alpha), and under the determined mechanism, obtaining an index factor A through the intercept of a formula and the ln (AR/beta E).

further, in the method for calculating pyrolysis kinetic parameters of carbonized combustible based on a unimodal pyrolysis curve according to the present invention, in step S3, the upper and lower floating by a certain percentage means that the upper and lower floating by 50% respectively.

Further, in the method for calculating pyrolysis kinetic parameters of carbonized combustible based on a unimodal pyrolysis curve of the invention, the step S4 is specifically: when PSO is optimized, the fitness function value phi is the deviation degree of a predicted value and an experimental value, and the specific calculation formula is as follows:

Wherein φ m, φ mlr, φ α and φ d α/dt represent the objective functions of mass loss, mass loss rate, conversion rate and conversion rate, respectively; n represents the number of experiments; n represents the number of experimental data points per experiment; CMLmod and CMLexp represent simulated and experimental values of accumulated mass loss; MLRmod and MLRexp represent a simulation value and an experimental value of the mass loss rate; α mod, α exp represent the simulated and experimental values of conversion; d α/dtmod, d α/dtexp represent the simulated and experimental values of the conversion rate; wCML, wMLR, w α, wd α/dt represent weighting factors for mass loss, rate of conversion.

further, in the method for calculating pyrolysis kinetic parameters of carbonized combustible based on a unimodal pyrolysis curve of the present invention, the respective experimental values in step S4 refer to the calculation formula of the experimental values measured by performing the FWO method in step S1 as follows:

CML＝m/m

MLR＝d(m/m)/dt＝

((m/m)-(m/m))/(t-t)

dα/dt＝Af(α)exp(-E/RT)；

Where m and mt represent the masses at the moment of reaction, m0 represents the initial mass of the sample, and m ∞ represents the final mass.

Compared with the existing pyrolysis solving method, the method has the following beneficial effects: 1) according to the method for solving the pyrolysis kinetic parameters of the carbonized combustible based on the single-peak pyrolysis curve, the reaction mechanism of the single-peak pyrolysis of the carbonized combustible is obtained on the basis of the single-peak pyrolysis experimental data of the carbonized combustible, a global optimization algorithm-Particle Swarm Optimization (PSO) is improved, and the difficulty in solving a high-order partial differential equation in the formula derivation process of theoretical analysis is reduced; the algorithm can automatically calculate and draw a chart, and has high calculation speed and small error in high-frequency iterative calculation; the reasonability of the experimental result is verified; the accuracy of the theoretical construction model is determined and applied to working condition research of larger size; manpower, material resources and financial resources are saved; 2) according to the method for solving the pyrolysis kinetic parameters of the carbonized combustible based on the unimodal pyrolysis curve, the accuracy of the analog value obtained by the algorithm is high; the application range is wide, and the method can be simultaneously applied to unimodal pyrolysis solution of multiple groups of heating rates and multiple carbonized combustibles, so that the experimental result can be effectively predicted; 3) according to the method for solving the pyrolysis kinetic parameters of the carbonized combustible based on the unimodal pyrolysis curve, the complex process of pyrolysis can be simulated on the basis of experimental data, actual condition optimization data is attached on the basis of experiments, and the accuracy of experimental results is effectively improved.

Drawings

The invention will be further described with reference to the accompanying drawings and examples, in which:

FIG. 1 is a flow chart of a method for solving pyrolysis kinetic parameters of carbonized combustible based on a unimodal pyrolysis curve according to the invention;

FIG. 2 is a graph comparing a simulation value and an experimental value of a conversion rate alpha of pyrolysis in a nitrogen environment, which is obtained by an algorithm in the first embodiment of the invention and takes a single-peak pyrolysis of an EPS thermal insulation material for an external wall as an example;

FIG. 3 is a graph comparing the simulated value and the experimental value of the conversion rate d α/dt of the single-peak pyrolysis of the EPS material in nitrogen environment according to the first embodiment of the present invention;

FIG. 4 is a graph comparing the simulated and experimental mass loss m/m0 values obtained by the algorithm in the first embodiment of the present invention, wherein the mass loss m/m0 is obtained by using a single-peak pyrolysis of an EPS insulation material for an external wall as an example;

FIG. 5 is a graph comparing the simulated and experimental mass loss rate d (m/m0)/dt obtained by the algorithm according to the first embodiment of the present invention, wherein the mass loss rate d (m/m0)/dt is obtained by pyrolysis in a nitrogen environment, for example, when the single-peak pyrolysis of the EPS as the exterior wall insulation material is performed;

FIG. 6 is a graph showing deviation degree obtained by pyrolysis in a nitrogen atmosphere for example of unimodal pyrolysis of the outer wall insulation material EPS in the first embodiment of the present invention.

Detailed Description

For a more clear understanding of the technical features, objects and effects of the present invention, embodiments of the present invention will now be described in detail with reference to the accompanying drawings.

The invention provides a method for solving pyrolysis kinetic parameters of carbonized combustible materials based on a unimodal pyrolysis curve, the flow of the method is shown in figure 1, and as can be seen from figure 1, the method for calculating the pyrolysis kinetic parameters of the carbonized combustible materials based on the unimodal pyrolysis curve comprises the following steps:

S1, charred combustible pyrolysis single-step reaction formula based on unimodal pyrolysis curve: solid → v carbon + (1-v) volatile matter, v is carbon generation rate, an equal conversion rate method FWO (Flynn-Wall-Ozawa method) is adopted to obtain the activation energy value E of pyrolysis of the carbonized combustible material under different conversion rates alpha, and finally the average value of the activation energy values E under different conversion rates is obtained.

In this embodiment, step S1 specifically includes: single-step reaction formula of carbonized combustible unimodal pyrolysis: solid → v char + (1-v) volatile matter, v is char formation rate, under the condition that the mechanism of the charring combustible unimodal pyrolysis reaction is unknown, the temperature of fixed conversion rate corresponding to the experimental result is measured by the equal conversion rate method FWO method under different temperature rising rates, and the temperature is determined by the formula: the method comprises the following steps of (a) obtaining a linear slope-1.052 (E/R) of ln beta and 1/T under various conversion rates according to a formula, obtaining an activation energy value E of unimodal pyrolysis of the carbonized combustible, and finally obtaining an average value of the activation energy values E under various conversion rates, wherein R represents a universal gas constant, beta represents a heating rate, g (alpha) represents a differential reaction mechanism function, and T represents an absolute temperature.

S2, under various reaction mechanisms, calculating an average value E of activation energy of pyrolysis of the carbonized combustible materials under different mechanisms by using a CR (coatings-Redfern method), comparing the average value E with a value obtained by a FWO method, finding a reaction mechanism function f (alpha) corresponding to the CR method when the activation energy value calculated by the FWO method is closest to the activation energy value calculated by the CR method, and then calculating a pre-index factor A.

In this embodiment, step S2 specifically includes: respectively under a plurality of different reaction mechanisms of the unimodal pyrolysis of the carbonized combustible, adopting a CR method to obtain an activation energy value E of the unimodal pyrolysis of the carbonized combustible, namely one reaction mechanism corresponding to one activation energy value is calculated by a formula: ln (g (alpha)/T2) ═ ln (AR/beta E) - (E/RT), under a plurality of assumed pyrolysis reaction mechanisms, calculating the slope of ln (g (alpha)/T2) and 1/T of a plurality of heating rates to obtain the activation energy value E of the pyrolysis single peak of the carbonized combustible, calculating the average value of a plurality of groups of heating rate activation energy values, comparing the average activation energy value calculated by the CR method with the average activation energy value calculated by the FWO method, finding the mechanism corresponding to the CR method when the activation energy value calculated by the FWO method is closest to the activation energy value calculated by the CR method, namely the reaction mechanism of the pyrolysis single peak of the carbonized combustible, determining a reaction mechanism function f (alpha), and after determining the mechanism, obtaining the pre-index factor A through the intercept of the formula ln (AR/beta E).

S3, floating the pre-index factor A and the activation energy value E (the activation energy value E calculated by FWO) up and down by a certain percentage respectively to obtain the optimization ranges of the pre-index factor A and the activation energy value E respectively, and according to the single-peak pyrolysis single-step reaction formula of the carbonized combustible, wherein the single-peak pyrolysis single-step reaction formula mainly generates carbon, liquid and gas products by pyrolysis of the carbonized combustible: solid → v char + (1-v) volatiles, so the char formation rate v is between 0 and 1, i.e. the optimum range of char formation rate v is set to (0,1), for example, a certain percentage of upper and lower fluctuation means 50% of each of upper and lower fluctuation, taking A as an example: the lower limit is 50% a and the upper limit is 150% a.

S4, optimizing the pre-indication factor A, the activation energy value E and the char generation rate v by adopting a particle swarm artificial intelligence optimization algorithm PSO program on the basis of the optimization range determined in the step S3, and outputting a conversion rate alpha, a conversion rate d alpha/dt, a mass loss m/m0 and a mass loss rate d (m/m0)/dt in the optimal solution, thereby completing the calculation of the pyrolysis kinetic parameters of the carbonized combustible; in the optimization process, the iterative update formula of the parameters is as follows:

α＝α+(dα/dt)×(t-t)，

(m/m)＝1-α×(1-v)，

(d(m/m)/dt)＝((m/m)-(m/m))/(t-t)，

dα/dt＝Af(α)exp(-E/RT)；

Wherein T represents time, m0 represent instantaneous mass and initial mass of charred combustible pyrolysis, R is universal gas constant, T is absolute temperature, and i represents update iteration number.

when PSO is optimized, the fitness function value phi is the deviation degree of a predicted value and an experimental value, and the specific calculation formula is as follows:

Wherein, the smaller the fitness function value phi is, the better, phi m, phi mlr, phi alpha and phi d alpha/dt respectively represent the objective functions of mass loss, mass loss rate, conversion rate and conversion rate; n represents the number of experiments; n represents the number of experimental data points per experiment; CMLmod and CMLexp represent simulated and experimental values of accumulated mass loss; MLRmod and MLRexp represent a simulation value and an experimental value of the mass loss rate; α mod, α exp represent the simulated and experimental values of conversion; d α/dtmod, d α/dtexp represent the simulated and experimental values of the conversion rate; wCML, wMLR, w α, wd α/dt represent weighting factors for mass loss, rate of conversion. These experimental values are the experimental values calculated by the FWO method in step S1 or the experimental values calculated by the CR method in step S2, and the formula of the experimental values is as follows:

CML＝m/m

MLR＝d(m/m)/dt＝

((m/m)-(m/m))/(t-t)

dα/dt＝Af(α)exp(-E/RT)；

Where m and mt represent the masses at the moment of reaction, m0 represents the initial mass of the sample, and m ∞ represents the final mass.