Solid material pyrolysis kinetic parameter acquisition method and system, electronic equipment

By constructing a dynamic mechanism function library and calculating the root mean square difference RMSE, and optimizing the kinetic parameters, the problem of low efficiency in obtaining pyrolysis kinetic parameters of solid materials in the prior art is solved, and higher accuracy and efficiency are achieved.

CN118072881BActive Publication Date: 2025-02-25SHANGHAI TOSUN TECH LTD

Patent Information

Application Number
CN202410203886.3
Authority / Receiving Office
CN · China
Patent Type
Patents(China)
Current Assignee / Owner
Filing Date
2024-02-23
Publication Date
2025-02-25
Estimated Expiration
2044-02-23

AI Technical Summary

Technical Problem

In the prior art, the acquisition efficiency of pyrolysis kinetic parameters of solid materials is low, mainly because the acquisition of the most general mechanism function depends on subjective inference or approximate inference, resulting in low accuracy and repeated verification.

Method used

By collecting the experimental curves of the percentage of mass loss of pyrolyte samples with temperature or time, a library of kinetic mechanism functions is constructed, and each kinetic mechanism function is traversed, the root mean square difference RMSE of each simulation curve and experimental curve is calculated, and after sorting, The kinetic mechanism function corresponding to the minimum RMSE is selected as the pyrolysis kinetic parameter of solid materials.

Benefits of technology

The accuracy and efficiency of obtaining kinetic parameters of pyrolysis of solid materials are improved. By optimizing the activation energy E and pre-finger factor A of the kinetic parameters, the error between the simulation curve and the experimental curve is reduced, and the theoretical minimum value is approximated.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN118072881B_ABST
    Figure CN118072881B_ABST
Patent Text Reader

Abstract

The present invention belongs to the technical field of pyrolysis kinetics, and specifically relates to a method for acquiring pyrolysis kinetic parameters of solid materials, an acquisition system, and electronic equipment, wherein the method for acquiring pyrolysis kinetic parameters of solid materials comprises: collecting an experimental curve of mass loss percentage of a pyrolysis sample varying with temperature or time; constructing a kinetic mechanism function library, and traversing each kinetic mechanism function to acquire a simulation curve of mass loss percentage corresponding to each kinetic mechanism function varying with temperature or time; 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 the mass loss percentage varying with temperature or time; and sorting each root mean square difference (RMSE), and taking the kinetic mechanism function corresponding to the minimum root mean square difference (RMSE) as the pyrolysis kinetic parameter of the solid material.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] The present invention belongs to the technical field of pyrolysis kinetic parameter analysis, and in particular relates to a method and system for acquiring pyrolysis kinetic parameters of solid materials, and electronic equipment. Background Art

[0002] 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.

[0003] In the related technologies, the most probable mechanism function is mostly obtained through subjective inference or approximate inference. This results in that the most probable mechanism function obtained is often inaccurate and requires repeated verification and inference, which seriously affects the efficiency of obtaining the pyrolysis kinetic parameters of solid materials. Summary of the invention

[0004] 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.

[0005] In order to solve the above technical problems, the present invention provides a method for obtaining pyrolysis kinetic parameters of solid materials, comprising:

[0006] Collect the experimental curve of the mass loss percentage of the pyrolysis sample as a function of temperature or time;

[0007] 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 changing with temperature or time;

[0008] 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

[0009] The root mean square differences (RMSE) are sorted, and the kinetic mechanism function corresponding to the minimum root mean square difference (RMSE) is used as the kinetic parameter of the pyrolysis kinetics of the solid material.

[0010] 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:

[0011] A storage unit configured to store an experimental curve of the mass loss percentage of the collected pyrolysis product sample versus temperature or time;

[0012] 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;

[0013] 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

[0014] The sorting unit is configured to sort the root mean square differences RMSE and use 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.

[0015] 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 the 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, each root mean square error RMSE value is optimized to the minimum, close to the theoretical minimum value, and the accuracy of the kinetic mechanism function sorting is guaranteed. 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.

[0016] Other features and advantages of the present invention will be described in the following description, and partly 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.

[0017] 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

[0018] In order to more clearly illustrate the specific implementation methods of the present invention or the technical solutions in the prior art, the drawings required for use in the specific implementation methods or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are some implementation methods of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0019] Figure 1 A step diagram showing a method for obtaining kinetic parameters of pyrolysis of solid materials according to some embodiments;

[0020] Figure 2 A schematic diagram showing the comparison between the most probable mechanism function No. 3 and experimental data in cases involved in some embodiments is shown;

[0021] Figure 3A schematic diagram showing the comparison between 10 non-most probable mechanism functions and experimental data in cases involved in some embodiments is shown;

[0022] Figure 4 A principle block diagram of a solid material pyrolysis kinetic parameter acquisition system according to some embodiments is shown;

[0023] Figure 5 Some electronic device principle block diagrams involved in the embodiments are shown. DETAILED DESCRIPTION

[0024] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution of the present invention will be clearly and completely described below in conjunction with the accompanying drawings. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0025] In the related art, the most probable mechanism function of the pyrolysis kinetic parameters of solid materials is often obtained through subjective inference or approximate inference. This results in the most probable mechanism function being inaccurate and requiring repeated verification and inference, which seriously affects the efficiency of obtaining the pyrolysis kinetic parameters of solid materials.

[0026] Therefore, at least one embodiment provides a method for obtaining kinetic parameters of pyrolysis of solid materials, comprising:

[0027] Collect the experimental curve of the mass loss percentage of the pyrolysis sample as a function of temperature or time;

[0028] 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 changing with temperature or time;

[0029] Calculate the root mean square error (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; and

[0030] The root mean square differences (RMSE) are sorted, and the kinetic mechanism function corresponding to the minimum root mean square difference (RMSE) is used as the kinetic parameter of the pyrolysis kinetics of the solid material.

[0031] Various non-limiting implementations of the embodiments of the present disclosure are described in detail below in conjunction with the accompanying drawings.

[0032] like Figure 1 As shown, some embodiments provide a method for obtaining kinetic parameters of pyrolysis of solid materials, including:

[0033] Step S101, collecting an experimental curve of the mass loss percentage of the pyrolysis sample changing with temperature or time;

[0034] 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;

[0035] Step S103, respectively calculating 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; and

[0036] Step S104, sorting the root mean square errors RMSE, and taking the kinetic mechanism function corresponding to the minimum root mean square error RMSE as the kinetic parameter of pyrolysis of the solid material.

[0037] Specifically, at least one embodiment constructs a function library including but not limited to the following 41 commonly used kinetic mechanism functions, as shown in Table 1 below, where α is the conversion rate, which refers to the percentage or fraction of a certain reactant converted; G(α) and f(α) are mechanism functions in integral and differential forms, respectively.

[0038] Table 1 Commonly used dynamic mechanism function libraries

[0039]

[0040]

[0041] In some embodiments, the method of 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 includes:

[0042] Set the initial activation energy E and pre-exponential factor A;

[0043] 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 while traversing, so as to obtain the optimized activation energy E and pre-exponential factor A corresponding to each kinetic mechanism function; the optimization formula is:

[0044]

[0045] Each kinetic mechanism function is used as the thermogravimetric simulation model, and the optimized activation energy E and pre-exponential factor A are used as the 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; where a is the number of iterations in the optimization process, a≥1.

[0046] Specifically, the initial activation energy E and the pre-exponential factor A are both given an initial value artificially, 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%.

[0047] 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 sequence 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 are used 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, and sorts the root mean square difference (RMSE) values, and the kinetic mechanism function corresponding to the minimum root mean square difference (RMSE) value is the most probable mechanism function. Since the parameters of the simulation model of each kinetic mechanism function, such as activation energy E and pre-exponential factor A, are optimized, the error 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 is reduced, that is, each root mean square error (RMSE) value is optimized to the minimum and approaches the theoretical minimum value, thereby ensuring the accuracy of the sorting of the kinetic mechanism functions. 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 thereby has high accuracy, thereby improving the efficiency of obtaining the kinetic parameters of the pyrolysis of solid materials.

[0048] The final outputs are 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 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.

[0049] In some embodiments, the method of collecting an experimental curve of the mass loss percentage of a pyrolysis sample as a function of temperature or time comprises:

[0050] In the non-isothermal pyrolysis mode, the mass loss data of the pyrolysis product sample corresponding to different pyrolysis heating rates are collected as samples; and

[0051] At each pyrolysis heating rate, the pyrolysis time and temperature 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.

[0052] 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:

[0053]

[0054] in:

[0055] i is the serial number of the kinetic mechanism function in the function library, 1≤i≤I, and I is the total number of kinetic mechanism functions;

[0056] t n is the nth moment of the thermogravimetric experiment;

[0057] n is the time sequence of the thermogravimetric experiment process, and 1≤n≤N.

[0058] N is the total number of moments in the thermogravimetric experiment process;

[0059] m exp (t n ) is the t during the thermogravimetric experiment n The weight of the pyrolysis sample at the time;

[0060] 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 pyrolysate sample at the time.

[0061] In some embodiments, the method of sorting the root mean square differences RMSE 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:

[0062] 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.

[0063] In some embodiments, the method for obtaining kinetic parameters of pyrolysis of solid materials further comprises:

[0064] 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 taken as the pyrolysis kinetic parameters of solid materials.

[0065] 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 smaller; while the comparison between the mechanism function No. 10, which is not the most probable mechanism function, and the experimental data is Figure 3 As shown, it is obvious that its root mean square error RMSE is large.

[0066] 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:

[0067] A storage unit configured to store an experimental curve of the mass loss percentage of the collected pyrolysis product sample versus temperature or time;

[0068] 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;

[0069] 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

[0070] The sorting unit is configured to sort the root mean square differences RMSE and use 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.

[0071] 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.

[0072] 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 specific contents of the aforementioned method for obtaining kinetic parameters of pyrolysis of solid materials can be referred to, and will not be elaborated here.

[0073] The electronic device in the embodiment of the present disclosure is described below from the perspective of hardware processing:

[0074] The embodiments of the present disclosure do not limit the specific implementation of the electronic device.

[0075] 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.

[0076] In some embodiments, a computer device or an industrial computer can also be used as a type of electronic device.

[0077] Figure 5 The structure shown does not constitute a limitation on the electronic device, and may include fewer or more components than shown in the figure, or combine certain components, or arrange the components differently.

[0078] In some embodiments, the readable storage medium or computer-readable storage medium includes at least one type of memory, and the memory includes a flash memory, a hard disk, a multimedia card, a card-type memory (such as an SD memory, etc.), a magnetic memory, a disk, an optical disk, etc. In some embodiments, it can be an internal storage unit of a computer device, such as a hard disk of the computer device. In other embodiments, the memory can also be an external storage device of a 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. Further, the memory can also include both an internal storage unit of a 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 the code of a computer program, etc., but also to temporarily store data that has been output or is to be output.

[0079] In some embodiments, the processor may be a central processing unit (CPU), a controller, a microcontroller, a microprocessor or other data processing chip for running program codes stored in a memory or processing data, such as executing a computer program.

[0080] Optionally, the computer device may further include a user interface, which may include a display (Display), an input unit such as a keyboard (Keyboard), and optionally, the user interface may also include a standard wired interface and a wireless interface. Optionally, in some embodiments, the display may be an LED display, a liquid crystal display, a touch-sensitive liquid crystal display, and an OLED (Organic Light-Emitting Diode) touch device, etc. Among them, the display may also be appropriately referred to as a display screen or a display unit, which is used to display information processed in the computer device and to display a visual user interface.

[0081] When the processor executes the program, the above Figure 1 The steps in the embodiment of the method for obtaining kinetic parameters of pyrolysis of solid materials shown are, for example, Figure 1 Alternatively, when the processor executes the computer program, the functions of each module or unit in the above-mentioned device embodiments are realized.

[0082] Some embodiments further provide a computer-readable storage medium configured to store any possible method for obtaining kinetic parameters of solid material pyrolysis.

[0083] Some embodiments further provide a computer-readable storage medium storing computer-readable instructions, which, when executed by at least one processor, enables the method for obtaining the kinetic parameters of pyrolysis of solid materials as described above to be executed, specifically, storing the experimental curve of the mass loss percentage of the collected pyrolysis sample as a function of temperature or time; constructing a kinetic mechanism function library, and traversing each kinetic mechanism function to obtain the simulation curve of the mass loss percentage corresponding to each kinetic mechanism function as a function of temperature or time; respectively calculating the root mean square difference RMSE between the simulation curve of each mass loss percentage as a function of temperature or time and the experimental curve of the mass loss percentage as a function of temperature or time; and sorting each root mean square difference RMSE, and taking the kinetic mechanism function corresponding to the minimum root mean square difference RMSE as the kinetic parameter of pyrolysis of solid materials. Please refer to the specific description of the method for obtaining the kinetic parameters of pyrolysis of solid materials, which will not be repeated here.

[0084] Some embodiments further provide a computer program product, including 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.

[0085] Some embodiments also provide a computer program product, including a computer-readable storage medium having a computer-readable program code stored thereon, the computer-readable program code including instructions, which 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 the collected pyrolysis sample as a function of 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 a function of temperature or time; calculate the root mean square difference (RMSE) between the simulation curve of each mass loss percentage as a function of temperature or time and the experimental curve of the mass loss percentage as a function of temperature or time; and 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.

[0086] In several embodiments provided in the present 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 schematic. For example, the flowcharts and block diagrams in the accompanying drawings show the possible architecture, 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 part of a code, and the module, a program segment or a part of a 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 a different order from the order 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 with a dedicated hardware-based system that performs a specified function or action, or can be implemented with a combination of dedicated hardware and computer instructions.

[0087] 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.

[0088] If the functions are implemented in the form of software function 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 part that contributes to the prior art, or the part of the technical solution, can be embodied in the form of a software product, which is stored in a storage medium and includes several instructions for enabling a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the methods described in various embodiments of the present invention.

[0089] Based on the above ideal embodiments of the present invention, the relevant staff can make various changes and modifications without departing from the technical concept of the present invention through the above description. The technical scope of the present invention is not limited to the contents of the specification, and its technical scope 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: Collect the experimental curve of the mass loss percentage of the pyrolysis sample as a function of 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 changing with 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 The root mean square errors (RMSE) are ranked, and the kinetic mechanism function corresponding to the minimum root mean square error (RMSE) is used as the kinetic parameter of the pyrolysis of solid materials. The method of 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 includes: 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 while traversing, and the optimized activation energy E and pre-exponential factor A corresponding to each kinetic mechanism function are obtained; 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 iteration step a+1 in the optimization process, which depends on the activation energy value of the 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 iteration 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 iteration 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 simulation 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 a 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.

2. The method for obtaining pyrolysis kinetic parameters of solid materials according to claim 1, characterized in that: Methods for collecting experimental curves of the percentage of mass loss of pyrolysis samples versus temperature or time include: In the non-isothermal pyrolysis mode, the mass loss data of the pyrolysis product sample corresponding to different pyrolysis heating rates are collected as samples; and 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.

3. The method for obtaining pyrolysis kinetic parameters of solid materials 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 changing with temperature or time and the experimental curve of the mass loss percentage changing with temperature or time includes: in: i is the serial number of the kinetic mechanism function in the function library, 1≤i≤I, and 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 pyrolysis 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 pyrolysate sample at the time.

4. The method for obtaining kinetic parameters of pyrolysis of solid materials according to claim 3, characterized in that: Methods for sorting the root mean square differences RMSE 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 include: All kinetic mechanism functions are sorted in descending order according to the root mean square error (RMSE), and the kinetic mechanism function ranked first is the most probable mechanism function; and The most probable mechanism function is taken as the kinetic parameter of the pyrolysis of solid materials.

5. The method for obtaining kinetic parameters of pyrolysis of solid materials according to claim 4, characterized in that: The method for obtaining kinetic parameters of pyrolysis of solid materials also includes: 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 pyrolysis kinetic parameters of solid materials.

6. 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 storage unit configured to store an experimental curve of the mass loss percentage of the collected pyrolysis product sample versus temperature or time; 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 A sorting unit is configured to sort the root mean square differences RMSE and use the kinetic mechanism function corresponding to the minimum root mean square difference RMSE as the kinetic parameter of pyrolysis of the solid material; Traverse each kinetic mechanism function to obtain the simulation curve of the mass loss percentage corresponding to each kinetic mechanism function changing with temperature or time, including: 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 while traversing, and the optimized activation energy E and pre-exponential factor A corresponding to each kinetic mechanism function are obtained; 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 iteration step a+1 in the optimization process, which depends on the activation energy value of the 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 iteration 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 iteration 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 simulation 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 a 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 storage unit is also 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.

8. 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 the pyrolysis kinetic parameters of a solid material according to any one of claims 1 to 5 is implemented.

9. 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 kinetic parameters of pyrolysis of solid materials as described in any one of claims 1-5.

10. 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 the pyrolysis kinetic parameters of a solid material according to any one of claims 1 to 5 is implemented.

Citation Information

Patent Citations

  • Method and system for measuring thermal resistance of dielectric film and metal interface

    CN113419120A

Cited By

  • Method and system for obtaining pyrolysis kinetics parameter of solid material and storage medium

    EP4607521A1