Solid material pyrolysis kinetic parameter acquisition method and system and computer device
By collecting and analyzing the experimental curves of pyrolyte samples, building and traversing the kinetic mechanism function library, calculating and sorting the root mean square difference RMSE, the problem of inaccurate acquisition of kinetic parameters in the existing technology is solved, and more efficient parameter acquisition is achieved.
Patent Information
- Application Number
- CN202510122873.8
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2024-02-23
- Publication Date
- 2025-08-15
AI Technical Summary
In the prior art, the most general mechanism function obtained by subjective or approximate inference methods is inaccurate, resulting in low efficiency in obtaining pyrolysis kinetic parameters of solid materials and requires repeated verification.
The experimental curve of the percentage of mass loss of pyrolyte samples changes with temperature or time was collected, a library of kinetic mechanism functions was constructed, and each kinetic mechanism function was traversed, the root mean square difference RMSE was calculated, and the kinetic mechanism function corresponding to the smallest root mean square difference RMSE was selected as the kinetic parameters of the pyrolysis of solid materials.
The accuracy and globality of the kinetic mechanism function are improved, and the accuracy of the most general mechanism function is ensured, thereby improving the acquisition efficiency of the pyrolysis kinetic parameters of solid materials.
Smart Images
Figure CN120496693A_ABST
Abstract
Description
[0001] This case is a divisional application for the method for obtaining the pyrolysis kinetic parameters of solid materials and the acquisition system and electronic equipment in patent application 2024102038863. Technical Field
[0002] The present invention belongs to the technical field of pyrolysis kinetic parameter analysis, and specifically relates to a method and system for obtaining pyrolysis kinetic parameters of solid materials, and electronic equipment. Background Art
[0003] Pyrolysis kinetics is a branch of science that uses thermal analysis and calorimetry to study the kinetics of physical changes or chemical reactions. In the study of pyrolysis kinetics, the most probable mechanism function is one of the important pyrolysis kinetic parameters of solid materials.
[0004] In related technologies, the most probable mechanism function is mostly obtained through subjective inference or approximate inference. This results in the obtained most probable mechanism function being often inaccurate and requiring repeated verification and inference, which seriously affects the efficiency of obtaining the pyrolysis kinetic parameters of solid materials. Summary of the Invention
[0005] The purpose of the present invention is to provide a method and system for obtaining kinetic parameters of pyrolysis of solid materials, and electronic equipment.
[0006] In order to solve the above technical problems, the present invention provides a method for obtaining pyrolysis kinetic parameters of solid materials, comprising:
[0007] Collect experimental curves showing the percentage of mass loss of pyrolysis samples changing with temperature or time;
[0008] Construct a kinetic mechanism function library and traverse each kinetic mechanism function to obtain a simulation curve of the mass loss percentage corresponding to each kinetic mechanism function as a function of temperature or time;
[0009] Calculate the root mean square difference (RMSE) between the simulation curve of each mass loss percentage varying with temperature or time and the experimental curve of each mass loss percentage varying with temperature or time; and
[0010] The root mean square differences (RMSEs) were ranked, and the kinetic mechanism function corresponding to the minimum root mean square difference (RMSE) was used as the kinetic parameter of the pyrolysis kinetics of the solid material.
[0011] In another aspect, the present invention further provides a system for obtaining kinetic parameters of pyrolysis of solid materials, comprising: a computer device, wherein the computer device is configured to include:
[0012] a storage unit configured to store an experimental curve of the mass loss percentage of the collected pyrolysis product sample changing with temperature or time;
[0013] A traversal simulation unit is configured to build a kinetic mechanism function library and traverse each kinetic mechanism function to obtain a simulation curve of the mass loss percentage corresponding to each kinetic mechanism function changing with temperature or time;
[0014] a calculation unit configured to calculate the root mean square difference (RMSE) between the simulation curve of each mass loss percentage varying with temperature or time and the experimental curve of each mass loss percentage varying with temperature or time; and
[0015] The sorting unit is configured to sort the root mean square differences (RMSEs) and use the kinetic mechanism function corresponding to the minimum root mean square difference (RMSE) as the pyrolysis kinetic parameter of the solid material.
[0016] The beneficial effect of the present invention is that the method for obtaining the kinetic parameters of pyrolysis of solid materials of the present invention traverses each kinetic mechanism function in sequence according to the function sequence number, and in the process of traversal, the kinetic parameter activation energy E and pre-exponential factor A of each kinetic mechanism function are optimized based on the quasi-Newton method, thereby reducing the error between the simulation curve of each mass loss percentage changing with temperature or time and the experimental curve of the mass loss percentage changing with temperature or time, that is, the root mean square error (RMSE) value is optimized to the minimum, approaching the theoretical minimum value, thereby ensuring the accuracy of the kinetic mechanism function sorting, and at the same time, since all the kinetic mechanism functions in the function library are traversed in sequence, the globality is guaranteed, and the most probable mechanism function obtained thereby has high accuracy, thereby improving the efficiency of obtaining the kinetic parameters of pyrolysis of solid materials.
[0017] Other features and advantages of the present invention will be described in the following description, and in part will become apparent from the description, or understood by practicing the present invention. The purpose and other advantages of the present invention are realized and obtained by the structures particularly pointed out in the description and the drawings.
[0018] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, preferred embodiments are given below and described in detail with reference to the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS
[0019] In order to more clearly illustrate the specific embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the specific embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.
[0020] Figure 1 A diagram showing the steps of a method for obtaining kinetic parameters of pyrolysis of solid materials according to some embodiments;
[0021] Figure 2 Schematic diagram showing the comparison between the most probable mechanism function No. 3 and experimental data in cases involved in some embodiments;
[0022] Figure 3 Schematic diagram showing the comparison between 10 non-most probable mechanism functions and experimental data in cases involved in some embodiments;
[0023] Figure 4 A principle block diagram of a solid material pyrolysis kinetic parameter acquisition system according to some embodiments is shown;
[0024] Figure 5 The figure shows the principle block diagram of the electronic device involved in some embodiments. DETAILED DESCRIPTION
[0025] To make the objectives, technical solutions, and advantages of the embodiments of the present invention more clear, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the embodiments described are only part of the embodiments of the present invention, not all of them. All other embodiments obtained by ordinary technicians in this field based on the embodiments of the present invention without making any creative efforts shall fall within the scope of protection of the present invention.
[0026] In related technologies, the most probable mechanism function in the pyrolysis kinetic parameters of solid materials is often obtained through subjective inference or approximate inference. This results in the obtained most probable mechanism function being often inaccurate and requiring repeated verification and inference, which seriously affects the efficiency of obtaining the pyrolysis kinetic parameters of solid materials.
[0027] Therefore, at least one embodiment provides a method for obtaining kinetic parameters of pyrolysis of solid materials, comprising:
[0028] Collect experimental curves showing the percentage of mass loss of pyrolysis samples changing with temperature or time;
[0029] Construct a kinetic mechanism function library and traverse each kinetic mechanism function to obtain a simulation curve of the mass loss percentage corresponding to each kinetic mechanism function as a function of temperature or time;
[0030] Calculating the root mean square error (RMSE) between the simulation curve of the mass loss percentage changing with temperature or time and the experimental curve of the mass loss percentage changing with temperature or time; and
[0031] The root mean square differences (RMSEs) were ranked, and the kinetic mechanism function corresponding to the minimum root mean square difference (RMSE) was used as the kinetic parameter of the pyrolysis kinetics of the solid material.
[0032] Various non-limiting implementations of the embodiments of the present disclosure are described in detail below with reference to the accompanying drawings.
[0033] like Figure 1 As shown, some embodiments provide a method for obtaining kinetic parameters of pyrolysis of solid materials, including:
[0034] Step S101, collecting an experimental curve showing the percentage of mass loss of a pyrolysis sample changing with temperature or time;
[0035] Step S102: constructing a kinetic mechanism function library and traversing each kinetic mechanism function to obtain a simulation curve of the mass loss percentage corresponding to each kinetic mechanism function changing with temperature or time;
[0036] Step S103, respectively calculating the root mean square difference (RMSE) between the simulation curve of each mass loss percentage varying with temperature or time and the experimental curve of each mass loss percentage varying with temperature or time; and
[0037] Step S104 , sorting the root mean square differences (RMSEs) and taking the kinetic mechanism function corresponding to the minimum root mean square difference (RMSE) as the kinetic parameter of the solid material pyrolysis kinetics.
[0038] Specifically, at least one embodiment constructs a function library including but not limited to the following 41 commonly used kinetic mechanisms, as shown in Table 1 below, where α is the conversion rate, which refers to the percentage or fraction of conversion of a certain reactant; G(α) and f(α) are the mechanism functions in integral and differential forms, respectively.
[0039] Table 1 Commonly used dynamic mechanism function libraries
[0040]
[0041]
[0042] In some embodiments, a method of traversing each kinetic mechanism function to obtain a simulation curve of the mass loss percentage corresponding to each kinetic mechanism function as a function of temperature or time includes:
[0043] Set the initial activation energy E and pre-exponential factor A;
[0044] Each kinetic mechanism function is traversed in turn, and the initial activation energy E and pre-exponential factor A are optimized based on the quasi-Newton method during the traversal to obtain the optimized activation energy E and pre-exponential factor A corresponding to each kinetic mechanism function; the optimization formula is:
[0045]
[0046] Each kinetic mechanism function is used as a thermogravimetric simulation model, and the optimized activation energy E and pre-exponential factor A are used as parameters of each thermogravimetric simulation model for simulation. The simulation curve of the mass loss percentage corresponding to each kinetic mechanism function changing with temperature or time is output; where a is the number of iterations in the optimization process, a≥1.
[0047] Specifically, the initial activation energy E and the pre-exponential factor A are both artificially given an initial value, for example but not limited to an initial value of 70000 for E and an initial value of e^10 for A, which may also be ±10%.
[0048] The method for obtaining the kinetic parameters of pyrolysis of solid materials in this embodiment traverses each kinetic mechanism function in sequence according to the function number, and in the process of traversal, optimizes the kinetic parameters activation energy E and pre-exponential factor A based on the quasi-Newton method, and then uses each kinetic mechanism function as a simulation model, and the optimized activation energy E and pre-exponential factor A as parameters of the simulation model to obtain a simulation curve of the mass loss percentage changing with temperature or time, and then uses the root mean square difference (RMSE) of the experimental curve of the mass loss percentage changing with temperature or time and the simulation curve of the mass loss percentage changing with temperature or time as the objective function, and obtains 41 root mean square difference (RMSE) values. The root mean square difference (RMSE) values are sorted, and the kinetic mechanism function corresponding to the minimum root mean square difference (RMSE) value is the most probable mechanism function. Since the parameters activation energy E and pre-exponential factor A of the simulation model of each kinetic mechanism function have been optimized, the error between the simulation curve of each mass loss percentage changing with temperature or time and the experimental curve of the mass loss percentage changing with temperature or time is reduced, that is, the root mean square error RMSE value has been optimized to the minimum, approaching the theoretical minimum value, thereby ensuring the accuracy of the sorting of the kinetic mechanism function. At the same time, since all the kinetic mechanism functions in the function library are traversed in sequence, the globality of the most probable mechanism function is guaranteed, and the most probable mechanism function obtained in this way has high accuracy, thereby improving the efficiency of obtaining the kinetic parameters of the pyrolysis of solid materials.
[0049] The final output is the most probable mechanism function, the kinetic parameter activation energy E of the most probable mechanism function, the kinetic parameter pre-exponential factor A of the most probable mechanism function, and the root mean square error (RMSE) of the experimental curve of the mass loss percentage changing with temperature or time as the kinetic parameters of the solid material pyrolysis. More accurate kinetic mechanism functions and kinetic parameters can help build more accurate application models for industrial processes such as pyrolysis, gasification, and combustion, and then use CFD (Computational Fluid Dynamics) to more accurately grasp the internal working conditions of the reactor and promote the development of related reactors.
[0050] In some embodiments, the method of collecting an experimental curve of the mass loss percentage of a pyrolysis sample versus temperature or time includes:
[0051] In the non-isothermal pyrolysis mode, the mass loss data of the pyrolysis sample corresponding to different pyrolysis heating rates are collected as samples; and
[0052] At each pyrolysis heating rate, the pyrolysis time and temperature values are used as input, and the mass loss percentage is used as output to obtain the experimental curve of the mass loss percentage changing with temperature or time.
[0053] In some embodiments, a method for calculating the root mean square difference (RMSE) between a simulation curve of each mass loss percentage varying with temperature or time and an experimental curve of each mass loss percentage varying with temperature or time includes:
[0054]
[0055] in:
[0056] i is the serial number of the kinetic mechanism function in the function library, 1≤i≤I, I is the total number of kinetic mechanism functions;
[0057] t n is the nth moment of the thermogravimetric experiment;
[0058] n is the time sequence of the thermogravimetric experiment process, and 1≤n≤N;
[0059] N is the total number of moments in the thermogravimetric experiment process;
[0060] m exp (t n ) is the t during the thermogravimetric experiment n The weight of the pyrolyzate sample at the time;
[0061] m E,A (t n ) is the optimized activation energy E, pre-exponential factor A parameters, during the thermogravimetric simulation process, t n The weight of the pyrolyzate sample at the time.
[0062] In some embodiments, the method of sorting the root mean square differences (RMSEs) and using the kinetic mechanism function corresponding to the minimum root mean square difference (RMSE) as the kinetic parameter of the pyrolysis kinetics of the solid material includes:
[0063] All mechanism functions are sorted in descending order according to the root mean square error (RMSE), and the mechanism function ranked first is the most probable mechanism function; and the most probable mechanism function is used as the pyrolysis kinetic parameter of the solid material.
[0064] In some embodiments, the method for obtaining kinetic parameters of pyrolysis of solid materials further comprises:
[0065] The activation energy E of the most probable mechanism function, the pre-exponential factor A of the most probable mechanism function and the minimum root mean square error RMSE are used as the kinetic parameters of the pyrolysis of solid materials.
[0066] Assuming that the mechanism function No. 3 obtained by the method for obtaining the kinetic parameters of the solid material pyrolysis in this embodiment is the most probable mechanism function, the comparison between it and the experimental data is as follows: Figure 2 As shown, it is obvious that its root mean square error RMSE is small; the comparison between the No. 10 mechanism function, which is not the most probable mechanism function, and the experimental data is as follows Figure 3 As shown in the figure, it is obvious that its root mean square error (RMSE) is large.
[0067] like Figure 4 As shown, some embodiments further provide a system for obtaining kinetic parameters of pyrolysis of solid materials, comprising: a computer device, wherein the computer device is configured to include:
[0068] a storage unit configured to store an experimental curve of the mass loss percentage of the collected pyrolysis product sample changing with temperature or time;
[0069] A traversal simulation unit is configured to build a kinetic mechanism function library and traverse each kinetic mechanism function to obtain a simulation curve of the mass loss percentage corresponding to each kinetic mechanism function changing with temperature or time;
[0070] a calculation unit configured to calculate the root mean square difference (RMSE) between the simulation curve of each mass loss percentage varying with temperature or time and the experimental curve of each mass loss percentage varying with temperature or time; and
[0071] The sorting unit is configured to sort the root mean square differences (RMSEs) and use the kinetic mechanism function corresponding to the minimum root mean square difference (RMSE) as the pyrolysis kinetic parameter of the solid material.
[0072] In some embodiments, the storage unit is further configured to store a kinetic mechanism function library, a simulation curve of mass loss percentage changing with temperature or time, a root mean square error (RMSE), and a ranking result of the root mean square error (RMSE).
[0073] Among them, the specific implementation functions of the storage unit, traversal simulation unit, calculation unit and sorting unit are realized in the computer device. The details can be referred to the content of the aforementioned method for obtaining the pyrolysis kinetic parameters of solid materials, which will not be repeated here.
[0074] The electronic device in the embodiment of the present disclosure is described below from the perspective of hardware processing:
[0075] The embodiments of the present disclosure do not limit the specific implementation of the electronic device.
[0076] like Figure 5As shown, some embodiments further provide a computer device / equipment / system, which may include but is not limited to: a processor, and a readable storage medium that communicates with the processor; wherein the readable storage medium is used to store a program for executing the method for obtaining the kinetic parameters of pyrolysis of solid materials as described above, and the program enables the processor to execute operations corresponding to the method for obtaining the kinetic parameters of pyrolysis of solid materials.
[0077] In some embodiments, a computer device or an industrial computer can also be used as a type of electronic device.
[0078] Figure 5 The structure shown does not limit the electronic device and may include fewer or more components than shown in the figure, or combine some components, or arrange the components differently.
[0079] In some embodiments, the readable storage medium or computer-readable storage medium includes at least one type of memory, including flash memory, hard disk, multimedia card, card-type memory (such as SD memory, etc.), magnetic memory, magnetic disk, optical disk, etc. In some embodiments, it can be an internal storage unit of a computer device, such as the hard disk of the computer device. In other embodiments, the memory can also be an external storage device of the computer device, such as a plug-in hard disk equipped on the computer device, a smart memory card (Smart Media Card, SMC), a secure digital (Secure Digital, SD) card, a flash card (Flash Card), etc. Furthermore, the memory can also include both an internal storage unit of the computer device and an external storage device. The memory can be used not only to store application software and various types of data installed in the computer device, such as computer program code, but also to temporarily store data that has been output or is to be output.
[0080] In some embodiments, the processor may be a central processing unit (CPU), a controller, a microcontroller, a microprocessor, or other data processing chip, used to run program codes stored in a memory or process data, such as executing a computer program.
[0081] Optionally, the computer device may further include a user interface, which may include a display and an input unit such as a keyboard. Optionally, the user interface may also include a standard wired interface or a wireless interface. Optionally, in some embodiments, the display may be an LED display, a liquid crystal display, a touch-sensitive liquid crystal display, or an OLED (Organic Light-Emitting Diode) touchscreen. The display may also be appropriately referred to as a display screen or display unit, and is used to display information processed by the computer device and to display a visual user interface.
[0082] When the processor executes the program, the above Figure 1 The steps in the embodiment of the method for obtaining the pyrolysis kinetic parameters of solid materials shown are, for example, Figure 1 Alternatively, when the processor executes the computer program, the functions of the modules or units in the above-mentioned device embodiments are realized.
[0083] Some embodiments further provide a computer-readable storage medium configured to store any of the above-described possible methods for obtaining kinetic parameters of solid material pyrolysis.
[0084] Some embodiments further provide a computer-readable storage medium storing computer-readable instructions that, when executed by at least one processor, cause the aforementioned method for obtaining kinetic parameters of pyrolysis of solid materials to be executed, specifically including: storing experimental curves showing the percentage of mass loss of collected pyrolysis samples as a function of temperature or time; constructing a library of kinetic mechanism functions and traversing each kinetic mechanism function to obtain a simulation curve showing the percentage of mass loss as a function of temperature or time corresponding to each kinetic mechanism function; calculating the root mean square difference (RMSE) between each simulation curve showing the percentage of mass loss as a function of temperature or time and the experimental curve showing the percentage of mass loss as a function of temperature or time; and sorting the RMSEs, selecting the kinetic mechanism function corresponding to the smallest RMSE as the kinetic parameter of pyrolysis of solid materials. For a detailed description of the method for obtaining kinetic parameters of pyrolysis of solid materials, please refer to the detailed description of the method; this description will not be repeated here.
[0085] Some embodiments further provide a computer program product, comprising a computer program or instructions, wherein when the computer program or instructions are executed on a computer, the computer is enabled to execute any of the above-mentioned possible methods for obtaining kinetic parameters of pyrolysis of solid materials.
[0086] Some embodiments also provide a computer program product comprising a computer-readable storage medium having computer-readable program code stored thereon, the computer-readable program code comprising instructions that cause at least one processor or at least one computer device to perform the following operations: store an experimental curve of the mass loss percentage of a collected pyrolysis sample as it changes with temperature or time; construct a kinetic mechanism function library and traverse each kinetic mechanism function to obtain a simulation curve of the mass loss percentage corresponding to each kinetic mechanism function as it changes with temperature or time; calculate the root mean square difference (RMSE) between the simulation curve of the mass loss percentage as it changes with temperature or time and the experimental curve of the mass loss percentage as it changes with temperature or time; and sort the root mean square differences (RMSE) and use the kinetic mechanism function corresponding to the minimum root mean square difference (RMSE) as the pyrolysis kinetic parameter of the solid material.
[0087] In the several embodiments provided in this application, it should be understood that the disclosed devices and methods can also be implemented in other ways. The device embodiments described above are merely illustrative. For example, the flowcharts and block diagrams in the accompanying drawings show the possible architectures, functions and operations of the devices, methods and computer program products according to multiple embodiments of the present invention. In this regard, each box in the flowchart or block diagram can represent a module, a program segment or a portion of code, and the module, program segment or a portion of code contains one or more executable instructions for implementing the specified logical function. It should also be noted that in some alternative implementations, the functions marked in the box can also occur in an order different from that marked in the accompanying drawings. For example, two consecutive boxes can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each box in the block diagram and / or flowchart, and the combination of boxes in the block diagram and / or flowchart, can be implemented using a dedicated hardware-based system that performs the specified function or action, or can be implemented using a combination of dedicated hardware and computer instructions.
[0088] In addition, the functional modules in the various embodiments of the present invention may be integrated together to form an independent part, or each module may exist independently, or two or more modules may be integrated to form an independent part.
[0089] If the functions are implemented in the form of software modules and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the portion that contributes to the prior art, or the portion of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes a number of instructions for causing a computer device (which can be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in various embodiments of the present invention.
[0090] With the above-described preferred embodiments of the present invention as a guide, and with reference to the above description, relevant personnel are fully capable of making various changes and modifications without departing from the technical scope of this invention. The technical scope of this invention is not limited to the contents of the specification and must be determined according to the scope of the claims.
Claims
1. A method for obtaining kinetic parameters of pyrolysis of solid materials, characterized in that: include: Construct a kinetic mechanism function library and traverse each kinetic mechanism function to obtain a simulation curve of the mass loss percentage corresponding to each kinetic mechanism function as a function of temperature or time; Calculate the root mean square difference (RMSE) between the simulation curve of each mass loss percentage changing with temperature or time and the experimental curve of each mass loss percentage changing with temperature or time; as well as Sort the root mean square differences (RMSEs) and use the kinetic mechanism function corresponding to the minimum RMSE as the kinetic parameter of the solid material pyrolysis; or When the kinetic mechanism function is the most probable mechanism function, the activation energy E corresponding to the most probable mechanism function, the pre-exponential factor A corresponding to the most probable mechanism function and the minimum root mean square error RMSE are also used as the kinetic parameters of the solid material pyrolysis.
2. The method for obtaining kinetic parameters of solid material pyrolysis according to claim 1, characterized in that: Methods for traversing each kinetic mechanism function to obtain a simulation curve of the mass loss percentage corresponding to each kinetic mechanism function changing with temperature or time include: Set the initial activation energy E and pre-exponential factor A; Each kinetic mechanism function is traversed in turn, and the initial activation energy E and pre-exponential factor A are optimized based on the quasi-Newton method during the traversal to obtain the optimized activation energy E and pre-exponential factor A corresponding to each kinetic mechanism function; the optimization goal is to minimize the function RMSE (E, A), and the optimization formula is as follows: in: E a+1 Represents the activation energy value of the iterative step a+1 in the optimization process, which depends on the activation energy value of step a and the gradient information; E a Represents the activation energy value of the iterative step a during the optimization process. When a=1, the initial activation energy value is taken. A a+1 Represents the pre-exponential factor value of the iterative step a+1 in the optimization process, which depends on the pre-exponential factor value and gradient information of the step a; A a Represents the pre-exponential factor value of the iterative step a in the optimization process. When a=1, the initial pre-exponential factor value is taken; RMSE(E a ,A a ) represents the root mean square difference between the simulated curve and the experimental curve of the mass loss percentage changing with temperature or time when the activation energy is Ea and the pre-exponential factor is Aa; and Each kinetic mechanism function is used as a thermogravimetric simulation model, and the optimized activation energy E and pre-exponential factor A are used as parameters of each thermogravimetric simulation model for simulation, and the simulation curve of the mass loss percentage corresponding to each kinetic mechanism function changing with temperature or time is output; a is the number of iterations in the optimization process, and a≥1.
3. The method for obtaining kinetic parameters of solid material pyrolysis according to claim 2, characterized in that: The calculation method of the root mean square difference RMSE between the simulation curve of each mass loss percentage varying with temperature or time and the experimental curve of each mass loss percentage varying with temperature or time includes: in: i is the serial number of the kinetic mechanism function in the function library, 1≤i≤I, I is the total number of kinetic mechanism functions; t n is the nth moment of the thermogravimetric experiment; n is the time sequence of the thermogravimetric experiment process, and 1≤n≤N; N is the total number of moments in the thermogravimetric experiment process; m exp (t n ) is the t during the thermogravimetric experiment n The weight of the pyrolyzate sample at the time; m E,A (t n ) is the optimized activation energy E, pre-exponential factor A parameters, during the thermogravimetric simulation process, t n The weight of the pyrolyzate sample at the time.
4. The method for obtaining kinetic parameters of solid material pyrolysis according to any one of claims 1 to 3, characterized in that: Also includes: An experimental curve of the percentage of mass loss of a pyrolysis sample versus temperature or time is collected, and the method includes: In the non-isothermal pyrolysis mode, the mass loss data of the pyrolysis sample corresponding to different pyrolysis heating rates are collected as samples; as well as At each pyrolysis heating rate, the pyrolysis time and temperature values are used as input, and the mass loss percentage is used as output to obtain the experimental curve of mass loss percentage changing with temperature or time.
5. A system for acquiring kinetic parameters of pyrolysis of solid materials, characterized in that: include: A computer device, the computer device being configured to include: A traversal simulation unit is configured to build a kinetic mechanism function library and traverse each kinetic mechanism function to obtain a simulation curve of the mass loss percentage corresponding to each kinetic mechanism function changing with temperature or time; a calculation unit configured to calculate the root mean square difference (RMSE) between the simulation curve of each mass loss percentage varying with temperature or time and the experimental curve of each mass loss percentage varying with temperature or time; and The sorting unit is configured to sort the root mean square differences (RMSEs) and use the kinetic mechanism function corresponding to the minimum root mean square difference (RMSE) as the pyrolysis kinetic parameter of the solid material.
6. The solid material pyrolysis kinetic parameter acquisition system according to claim 5, characterized in that: The traversal simulation unit constructs a kinetic mechanism function library and traverses each kinetic mechanism function to obtain a simulation curve of the mass loss percentage corresponding to each kinetic mechanism function changing with temperature or time, that is, Methods for traversing each kinetic mechanism function to obtain a simulation curve of the mass loss percentage corresponding to each kinetic mechanism function changing with temperature or time include: Set the initial activation energy E and pre-exponential factor A; Each kinetic mechanism function is traversed in turn, and the initial activation energy E and pre-exponential factor A are optimized based on the quasi-Newton method during the traversal to obtain the optimized activation energy E and pre-exponential factor A corresponding to each kinetic mechanism function; the optimization goal is to minimize the function RMSE (E, A), and the optimization formula is as follows: in: E a+1 Represents the activation energy value of the iterative step a+1 in the optimization process, which depends on the activation energy value of step a and the gradient information; E a Represents the activation energy value of the iterative step a during the optimization process. When a=1, the initial activation energy value is taken. A a+1 Represents the pre-exponential factor value of the iterative step a+1 in the optimization process, which depends on the pre-exponential factor value and gradient information of the step a; A a Represents the pre-exponential factor value of the iterative step a in the optimization process. When a=1, the initial pre-exponential factor value is taken; RMSE(E a ,A a ) represents the root mean square difference between the simulated curve and the experimental curve of the mass loss percentage changing with temperature or time when the activation energy is Ea and the pre-exponential factor is Aa; and Each kinetic mechanism function is used as a thermogravimetric simulation model, and the optimized activation energy E and pre-exponential factor A are used as parameters of each thermogravimetric simulation model for simulation, and the simulation curve of the mass loss percentage corresponding to each kinetic mechanism function changing with temperature or time is output; a is the number of iterations in the optimization process, and a≥1.
7. The solid material pyrolysis kinetic parameter acquisition system according to claim 6, characterized in that: The calculation unit calculates the root mean square difference RMSE of the simulation curve of each mass loss percentage changing with temperature or time and the experimental curve of each mass loss percentage changing with temperature or time, that is, The calculation method of the root mean square difference RMSE between the simulation curve of each mass loss percentage varying with temperature or time and the experimental curve of each mass loss percentage varying with temperature or time includes: in: i is the serial number of the kinetic mechanism function in the function library, 1≤i≤I, I is the total number of kinetic mechanism functions; t n is the nth moment of the thermogravimetric experiment; n is the time sequence of the thermogravimetric experiment process, and 1≤n≤N; N is the total number of moments in the thermogravimetric experiment process; m exp (t n ) is the t during the thermogravimetric experiment n The weight of the pyrolyzate sample at the time; m E,A (t n ) is the optimized activation energy E, pre-exponential factor A parameters, during the thermogravimetric simulation process, t n The weight of the pyrolyzate sample at the time.
8. The solid material pyrolysis kinetic parameter acquisition system according to claim 7, characterized in that: The ranking unit ranks the root mean square differences RMSE and uses the kinetic mechanism function corresponding to the minimum root mean square difference RMSE as the kinetic parameter of the solid material pyrolysis, that is, The kinetic mechanism function ranked first among all kinetic mechanism functions according to the root mean square error (RMSE) is defined as the most probable mechanism function; and Taking the most probable mechanism function as the kinetic parameter of the pyrolysis of solid materials; or, When the kinetic mechanism function is the most probable mechanism function, the activation energy E corresponding to the most probable mechanism function, the pre-exponential factor A corresponding to the most probable mechanism function and the minimum root mean square error RMSE are also used as the kinetic parameters of the solid material pyrolysis.
9. The solid material pyrolysis kinetic parameter acquisition system according to any one of claims 5 to 8, characterized in that: Also includes: a storage unit configured to store an experimental curve of the mass loss percentage of the collected pyrolysis product sample changing with temperature or time; Wherein, the storage unit stores the experimental curve, that is, In the non-isothermal pyrolysis mode, the mass loss data of the pyrolysis sample corresponding to different pyrolysis heating rates are collected as samples; as well as At each pyrolysis heating rate, the pyrolysis time and temperature values are used as input, and the mass loss percentage is used as output. The experimental curve of the mass loss percentage changing with temperature or time is obtained and stored.
10. The solid material pyrolysis kinetic parameter acquisition system according to claim 9, characterized in that: The storage unit is further configured to store a kinetic mechanism function library, a simulation curve of mass loss percentage changing with temperature or time, a root mean square error (RMSE), and a ranking result of the root mean square error (RMSE).
11. A method for traversing various kinetic mechanism functions to obtain a simulation curve of the mass loss percentage corresponding to each kinetic mechanism function as a function of temperature or time, comprising: Set the initial activation energy E and pre-exponential factor A; Traverse each kinetic mechanism function in turn, and optimize the initial activation energy E and pre-exponential factor A based on the quasi-Newton method while traversing, and obtain the optimized activation energy E and pre-exponential factor A corresponding to each kinetic mechanism function; The optimization goal is to minimize the function RMSE(E,A). The optimization formula is as follows: in: E a+1 Represents the activation energy value of the iterative step a+1 in the optimization process, which depends on the activation energy value of step a and the gradient information; E a Represents the activation energy value of the iterative step a during the optimization process. When a=1, the initial activation energy value is taken. A a+1 Represents the pre-exponential factor value of the iterative step a+1 in the optimization process, which depends on the pre-exponential factor value and gradient information of the step a; A a Represents the pre-exponential factor value of the iterative step a in the optimization process. When a=1, the initial pre-exponential factor value is taken; RMSE(E a ,A a ) represents the root mean square difference between the simulated curve and the experimental curve of the mass loss percentage changing with temperature or time when the activation energy is Ea and the pre-exponential factor is Aa; and Each kinetic mechanism function is used as a thermogravimetric simulation model, and the optimized activation energy E and pre-exponential factor A are used as parameters of each thermogravimetric simulation model for simulation, and the simulation curve of the mass loss percentage corresponding to each kinetic mechanism function changing with temperature or time is output; a is the number of iterations in the optimization process, and a≥1.
12. A computer-readable storage medium having a computer program / instruction stored thereon, characterized in that: When the computer program / instruction is executed by a processor, the method for obtaining pyrolysis kinetic parameters of solid materials according to any one of claims 1 to 4 or the method according to claim 11 is implemented.
13. A computer device / equipment / system, characterized in that: A memory, a processor, and a computer program stored in the memory, characterized in that the processor executes the computer program to implement the method for obtaining the pyrolysis kinetic parameters of solid materials according to any one of claims 1 to 4 or the method according to claim 11.
14. A computer program product comprising a computer program / instructions, characterized in that When the computer program / instruction is executed by a processor, the method for obtaining pyrolysis kinetic parameters of solid materials according to any one of claims 1 to 4 or the method according to claim 11 is implemented.