Carbonization reaction solving method, solver, computer equipment and medium

Abstract

The embodiment of the invention relates to the technical field of carbonization tower process data processing, in particular to a carbonization reaction solving method, a solver, computer equipment and a computer storage medium. The method mainly comprises the following steps: circularly calling a reaction solution equation, an absorption mass transfer equation, a crystallization reaction equation and a concentration solution equation in a carbonization reaction process until the carbon dioxide absorption rate and the substance mass transfer rate converge, and outputting the newest reactant component concentration value, the newest crystallization reaction rate and the newest crystal particle size distribution data, and outputting the reactant component concentration value, the crystallization reaction rate and the crystal particle size distribution data after the reaction is balanced, so that control equipment corresponding to the carbonization reaction is provided for parameter regulation and control of the carbonization reaction equipment. According to the method, automatic solving of the carbonization reaction process data is realized, the carbonization reaction process data can be accurately obtained, and data support is provided for parameter regulation and control of control equipment.

Description

Technical Field

[0001] This disclosure relates to the field of carbonization tower process data processing technology, and in particular to a carbonization reaction solution method, solver, computer equipment and computer storage medium. Background Technology

[0002] The carbonization reaction process includes multiple processes such as mass transfer, heat transfer, chemical reaction, and crystallization, and includes gas, liquid, and solid multiphase systems. Solving the carbonization reaction process data is of great significance for carbonization reaction research and production guidance.

[0003] However, currently there is a lack of multi-parameter coupling solution capability for the carbonization production process, a lack of automated solution for carbonization reaction process data, an inability to accurately obtain carbonization reaction process data, and the parameter adjustment of the control equipment corresponding to the carbonization reaction relies on manual experience, making it difficult to accurately match the dynamic requirements of the reaction process. Summary of the Invention

[0004] Therefore, it is necessary to provide an automated solution method, solver, computer device, and computer storage medium for at least one of the above-mentioned technical problems.

[0005] In a first aspect, embodiments of this disclosure provide a method for solving carbonization reactions, the method comprising the following steps:

[0006] Based on the preset reaction equations and the initial concentrations, temperatures, reaction rate constants, and ionic strengths of the reactant components, the absorption reaction rate and the hydrolysis reaction rate are determined.

[0007] Based on the preset absorption mass transfer equation and absorption reaction rate, determine the carbon dioxide absorption rate and mass transfer rate corresponding to the concentration values ​​of reactant components obtained during the carbonization reaction.

[0008] Based on the preset crystallization reaction equation and crystallization reaction parameters, determine the crystallization reaction rate and crystal grain size distribution data;

[0009] Based on the preset concentration solution equation, crystallization reaction rate, hydrolysis reaction rate, and reactant component parameters, determine the concentration values ​​of reactant components after carbonization reaction;

[0010] During the carbonization reaction, the reaction solution equation, absorption mass transfer equation, crystallization reaction equation, and concentration solution equation are called in a loop until the carbon dioxide absorption rate and mass transfer rate converge. The latest reactant component concentration values, the latest crystallization reaction rate, and the latest crystal particle size distribution data are output to provide the corresponding control equipment for parameter adjustment of the carbonization reaction.

[0011] In some embodiments, the carbonation reaction is a carbonation reaction that occurs inside a soda ash carbonation tower, which includes multiple trays. The solution method further includes the following steps:

[0012] Based on the preset heat balance equation, the heat data of a single tray is calculated, and the plate-by-plate calculation is performed according to the preset overall tower balance equation to obtain the temperature distribution data, carbon dioxide absorption rate distribution data, crystal particle size distribution data, and supersaturation distribution data of each tray in the soda ash carbonization tower.

[0013] In some embodiments, the total balance equations include a pressure balance equation for material components, a moment balance equation for crystal size distribution, and a mass balance equation for material components.

[0014] In some embodiments, determining the crystallization reaction rate and crystal grain size distribution data based on a preset crystallization reaction equation and crystallization reaction parameters includes the following steps:

[0015] Based on the preset supersaturation calculation equation, the crystal growth rate and nucleation rate are determined;

[0016] Based on the crystallization reaction equation, crystal growth rate, nucleation rate, and crystallization reaction parameters, the crystallization reaction rate and crystal size distribution data are determined.

[0017] In some embodiments, the solution method further includes: determining the kinetic parameters and equilibrium constant of the carbonization reaction based on the reaction solution equation.

[0018] In a second aspect, embodiments of this disclosure provide a solver, the solver code of which is written in MATLAB, and the solver is used to implement the steps of the carbonization reaction solving method provided in any embodiment of the first aspect of this disclosure.

[0019] In a third aspect, embodiments of this disclosure provide a computer device including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the steps of the carbonization reaction solving method provided in any embodiment of the first aspect of this disclosure.

[0020] In a fourth aspect, embodiments of the present disclosure provide a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the carbonization reaction solving method provided in any embodiment of the first aspect of the present disclosure.

[0021] The aforementioned carbonization reaction solution method, solver, computer equipment, and storage medium, through cyclically calling the reaction solution equation, absorption mass transfer equation, crystallization reaction equation, and concentration solution equation during the carbonization reaction process until the carbon dioxide absorption rate and mass transfer rate converge, output the latest reactant component concentration values, the latest crystallization reaction rate, and the latest crystal particle size distribution data. It also outputs the reactant component concentration values, crystallization reaction rate, and crystal particle size distribution data after reaction equilibrium to provide the corresponding control equipment for parameter adjustment of the carbonization reaction equipment. This achieves automated solution of carbonization reaction process data, accurately obtains carbonization reaction process data, and provides data support for parameter adjustment of control equipment. Attached Figure Description

[0022] Figure 1 This is a flowchart illustrating the carbonization reaction solution method in some embodiments;

[0023] Figure 2 This is a diagram showing the internal structure of a computer device in some embodiments. Detailed Implementation

[0024] To make the technical solutions and advantages of this disclosure clearer, the embodiments and related technical content of this disclosure will be further described in detail below with reference to the accompanying drawings and text description. It should be understood that the embodiments described below are only used to explain the technical solutions of the embodiments of this disclosure and are not intended to limit more possible implementations of this disclosure.

[0025] It should be noted that relational terms such as "first" and "second" appearing in this document are used only to distinguish things, states, or actions, and do not necessarily indicate or imply relative importance or order. The terms "including," "comprising," or any other variations thereof are used to indicate non-exclusive inclusion, and the included objects may not be limited to those listed in this document. The terms "multiple" or other variations are used to indicate that the number of objects is two or more.

[0026] In a first aspect, embodiments of this disclosure provide a method for solving a carbonization reaction. The solving device applied to the carbonization reaction can be, but is not limited to, various personal computers, laptops, smartphones, tablets, etc. The processor included in the solving device can be implemented using at least one hardware form selected from programmable logic arrays (PLAs), field-programmable gate arrays (FPGAs), digital signal processors (DSPs), application-specific integrated circuits (ASICs), general-purpose processors, or other programmable logic devices. The solving device is used to solve reaction data of the carbonization reaction process within a carbonization tower. Specifically, the solving device includes a solver for solving the reaction data.

[0027] In related technologies, the carbonation tower is the main equipment used in the carbonation process of soda ash (sodium carbonate) production. Its function is to absorb carbon dioxide from the Solvay process's ammonia brine or the Hou's process's ammonia mother liquor to produce sodium bicarbonate, which is then calcined to decompose into sodium carbonate. The large amount of heat released during the carbonation process needs to be removed from the tower, and the tower's structure must meet the requirements of absorption, crystallization, and cooling processes. This can be achieved by increasing the contact area between carbon dioxide and the ammonia brine through the tower's packing or sieve plate structure, thereby increasing the reaction rate. The temperature gradient and fluid dynamics design within the tower (such as stirring and countercurrent operation) ensure uniform growth of sodium bicarbonate crystals and prevent agglomeration.

[0028] Carbonization towers involve multiple processes, including mass transfer, heat transfer, chemical reactions, and crystallization, forming a multiphase system of gas, liquid, and solid. The presence of supersaturated solutions and a large amount of suspended solid crystals within the tower inevitably leads to scaling on internal components and heat exchange surfaces. Furthermore, crystal accumulation can easily clog fluid channels. Solvay carbonization towers are traditional equipment suitable for carbonization operations.

[0029] Typically, carbonization towers are vertical cylindrical structures, with packing or sieve plates inside to enhance gas-liquid contact. The reaction process in the carbonization tower is divided into three stages:

[0030] Phase 1: Absorption Section (Upper part of the carbonization tower), ammonia brine flows from top to bottom, absorbing CO2 to generate HCO3-. The temperature is controlled at 30–40℃ to prevent NH3 volatilization.

[0031] Stage Two: Reaction Crystallization Section (Middle of the Carbonization Tower), NaHCO3 becomes supersaturated and precipitates, forming fine crystals. Crystal growth is controlled by stirring or fluid design.

[0032] Phase 3: Separation section (lower part), NaHCO3 slurry is discharged from the bottom of the tower, and unreacted CO2 and NH3 are recovered from the top of the tower.

[0033] The main substances involved in the carbonization reaction in the carbonization tower include:

[0034] (1) CO2 (carbon dioxide): It participates in the reaction as a carbon source, and reacts with ammonia (NH3) under alkaline conditions to produce carbamate (NH2COO). - );

[0035] (2) NH3 (ammonia): provides an alkaline environment, promoting CO2 dissolution and subsequent reactions;

[0036] (3)NH2COO - (Carbamate): An intermediate product that further hydrolyzes to produce bicarbonate (HCO3-). - ) and ammonia (NH3);

[0037] (4)NH4 +(Ammonium ion): Formed by the protonation of NH3, reacting with HCO3- - They can combine to form ammonium bicarbonate (NH4HCO3);

[0038] (5) HCO3 - (Bicarbonate): One of the final products, which can crystallize into sodium bicarbonate (NaHCO3) or ammonium bicarbonate (NH4HCO3).

[0039] The application of the carbonization reaction solution method to the solution equipment is illustrated as an example. In some embodiments, such as... Figure 1 As shown, the carbonization reaction solution method includes steps S101 to S105 that can be performed by the solution device. Each step is described in detail below.

[0040] Step S101: Determine the absorption reaction rate and hydrolysis reaction rate based on the preset reaction solution equation and the initial values ​​of the concentration of the reactant components, temperature, reaction rate constant and ionic strength.

[0041] The reaction equations are used to solve for the absorption and hydrolysis rates of the carbonization reaction. The absorption rate refers to the rate of absorption within the liquid film near the gas-liquid interface, while the hydrolysis rate refers to the rate of hydrolysis occurring in the bulk liquid phase. The absorption rate is considered a fast reaction, while the hydrolysis rate is considered a slow reaction.

[0042] The chemical expressions related to the absorption reaction are as follows:

[0043]

[0044] The chemical formulas related to the hydrolysis reaction are:

[0045] Carbonization reactions also include crystallization reactions that occur in the liquid bulk, and their relevant chemical formulas are as follows:

[0046] The reactant components include CO2, NH3, and NH2COO. - NH4 + and HCO3 - By inputting the initial concentrations of the reactant components, temperature, reaction rate constants, and ionic strengths into the reaction equation, the absorption reaction rate and hydrolysis reaction rate can be obtained.

[0047] The reaction rate constant includes the equilibrium constants K1, K3, and Kh. In some specific cases, K1 = 1.18e-19, %K3 = 6.697e-13, and %Kh = 4.634e-7, where e is the natural constant.

[0048] In some specific embodiments, the initial concentration values ​​of reactant components can be obtained in different ways depending on the application scenario. In experimental simulation scenarios, during the simulated carbonization reaction process in a carbonization tower, the initial concentration values ​​can be set and input by the experimenter or automatically generated by a computer program. In actual carbonization tower production scenarios, the initial concentration values ​​can be obtained through sensor devices, such as pH meters, conductivity meters, infrared spectroscopy (IR), or Raman spectroscopy, to monitor CO2, NH3, and HCO3 in the liquid phase in real time. - The concentration can also be determined by sampling the reactants of the carbonization reaction, or by calculating based on the feed ratio.

[0049] Temperature can be obtained using temperature acquisition equipment, or it can be manually input or automatically input by a computer.

[0050] Ionic strength can be calculated from ion concentration, or it can be obtained by measuring the conductivity of the solution with a conductivity meter and combining it with empirical formulas. In the carbonization tower production process summary, ionic strength can also be determined based on the feed composition (the concentration value of the components).

[0051] Step S102: Based on the preset absorption mass transfer equation and absorption reaction rate, determine the carbon dioxide absorption rate and mass transfer rate corresponding to the concentration values ​​of reactant components obtained during the carbonization reaction.

[0052] The reactant components obtained during the carbonization reaction include CO2 and NH2COO. - NH3 and NH4 + The determined mass transfer rate includes NH2COO. - The mass transfer rate of NH3 and NH4 + The mass transfer rate.

[0053] In some specific embodiments, determining the carbon dioxide absorption rate and the mass transfer rate further includes determining the liquid film thickness based on the molecular diffusion coefficient, the liquid phase mass transfer coefficient, and / or the gas phase mass transfer coefficient.

[0054] In some specific embodiments, determining the carbon dioxide absorption rate and the mass transfer rate may further include:

[0055] The carbon dioxide absorption rate is determined based on the liquid film thickness, carbon dioxide partial pressure, and absorption reaction rate.

[0056] The mass transfer rate is determined based on the liquid film thickness and the absorption reaction rate.

[0057] In some specific embodiments, the absorption mass transfer equation may include an absorption equation and a mass transfer equation, wherein the absorption equation is used to solve for the absorption rate of CO2, and the mass transfer equation is used to solve for the mass transfer rate of ions in the reactant components from the membrane to the bulk.

[0058] The carbon dioxide absorption process during the carbonization reaction takes place within the liquid membrane, and the mass transfer and reaction within the membrane can be described by the following differential equation:

[0059]

[0060] rA represents the absorption reaction rate; correspondingly, the mass transfer rate of NH2COO during the carbonization reaction can be determined. - NH3, NH4 + The reaction rate equations for equal substances are determined. Based on the preset boundary conditions, the absorption equation and mass transfer equation are determined.

[0061] The preset boundary conditions are:

[0062]

[0063] In some specific embodiments, the parameters input to the absorption mass transfer equation may include: molecular diffusion coefficient, liquid phase mass transfer coefficient, gas phase mass transfer coefficient, CO2 partial pressure at the CO2 bulk gas phase, CO2 partial pressure at the CO2 gas-liquid interface, CO2 concentration at the liquid bulk phase, and NH2COO at the liquid bulk phase. - Concentration, NH3 concentration in the bulk liquid phase, NH4 concentration in the bulk liquid phase + Concentration and absorption reaction rate.

[0064] By inputting the molecular diffusion coefficient, liquid-phase mass transfer coefficient, gas-phase mass transfer coefficient, carbon dioxide partial pressure, and reactant component concentrations into the absorption mass transfer equation, the absorption rate of carbon dioxide and the mass transfer rate can be obtained. The mass transfer rate includes the concentrations of NH₂COO₂. - Mass transfer rate, NH3 mass transfer rate and NH4 + Mass transfer rate.

[0065] Specifically, the parameters input into the equation can be obtained through experimental calculations or actual data collection. There are no restrictions on this, and those skilled in the art can obtain them according to the actual needs of the scenario.

[0066] Step S103: Determine the crystallization reaction rate and crystal size distribution data based on the preset crystallization reaction equation and crystallization reaction parameters.

[0067] Crystallization reaction parameters may include volume, crystal growth rate, liquid flow rate, particle size distribution function handle at the inlet, nucleation rate, gas holdup, crystal surface shape factor, crystal density, NaHCO3 molecular weight, and liquid phase residence time. Crystallization particle size distribution data may include particle size / length values ​​and particle size density distribution data at the corresponding particle size.

[0068] By inputting the crystallization reaction parameters into the preset crystallization reaction equation, data on the crystallization reaction rate and crystal size distribution can be obtained.

[0069] In some specific embodiments, the parameters input to the crystallization reaction equation include: volume, crystal growth rate, liquid flow rate, particle size distribution function handle at the inlet, nucleation rate, gas holdup, crystal surface shape factor, crystal density, sodium bicarbonate molecular weight, and liquid phase residence time.

[0070] In some specific embodiments, the preset crystallization reaction equations may include crystallization rate equations, particle number balance equations, and crystallization kinetic equations. Solving these equations using the relevant solution modules of the solving equipment yields the crystallization reaction rate and crystal particle size distribution data. The supersaturation used in the solution process is obtained from the input temperature, sodium ion concentration, and bicarbonate ion concentration. The ion concentration can be calculated from the feed components or from subsequent component concentration modules.

[0071] In some specific embodiments, the expression for the crystallization reaction rate can be determined based on the hydrolysis and crystallization reaction equations within the liquid phase, as well as the hydrolysis reaction rate: After performing the Laplace transform, the moment method can be used to express the moments of each order of particle size distribution as functions of the moments of each order of nin at the inlet, from which the expression for the crystallization rate can be derived:

[0072] Step S104: Determine the concentration values ​​of reactant components after carbonization reaction based on the preset concentration solution equation, crystallization reaction rate, hydrolysis reaction rate, and reactant component parameters.

[0073] Reactant component parameters may include liquid film volume, specific surface area, gas content, flow rate, and mass transfer rate of each reactant component. In some specific embodiments, the concentration of components with higher concentrations, such as NH₂COO, is calculated based on the moment of moment conservation equation, mass conservation equation, and charge conservation equation. - The concentrations of components such as NH3, as well as the concentrations of components with smaller concentrations calculated based on the chemical equilibrium equation.

[0074] In some specific embodiments, the input parameters of the preset concentration solution equation include specific surface area, gas holdup, flow rate, hydrolysis reaction rate, crystallization reaction rate, mass transfer rate of each reactant component, and inlet concentration value of each reactant component.

[0075] By inputting the crystallization reaction rate, hydrolysis reaction rate, and reactant component parameters into a preset concentration calculation equation, the concentration values ​​of the reactant components after the carbonization reaction can be determined.

[0076] Specifically, the parameters of the input equation can be obtained through experimental calculations, actual data collection, or by solving the equation through the aforementioned steps. There are no restrictions here, and those skilled in the art can obtain them according to the actual needs of the scenario.

[0077] Step S105: During the carbonization reaction, the reaction solution equation, absorption mass transfer equation, crystallization reaction equation, and concentration solution equation are called in a loop until the carbon dioxide absorption rate and mass transfer rate converge. The latest reactant component concentration values, the latest crystallization reaction rate, and the latest crystal particle size distribution data are output to provide the corresponding control equipment for parameter adjustment of the carbonization reaction.

[0078] In some specific embodiments, the solution method further includes:

[0079] When the carbon dioxide absorption rate and mass transfer rate converge, the entire carbonization process is simulated to solve for the concentration values ​​of each component, the crystallization reaction rate, and the crystal grain size distribution data after the carbonization reaction reaches equilibrium.

[0080] In the above-mentioned carbonization reaction solution method, the reaction solution equation, absorption mass transfer equation, crystallization reaction equation, and concentration solution equation are cyclically called during the carbonization reaction process until the carbon dioxide absorption rate and mass transfer rate converge. The latest reactant component concentration values, the latest crystallization reaction rate, and the latest crystal particle size distribution data are output. The reactant component concentration values, crystallization reaction rate, and crystal particle size distribution data after the reaction equilibrium is output to provide the corresponding control equipment for the carbonization reaction to adjust the parameters of the carbonization reaction equipment. This realizes the automated solution of carbonization reaction process data, can accurately obtain carbonization reaction process data, and provides data support for the parameter adjustment of control equipment.

[0081] In some embodiments, the carbonation reaction is a carbonation reaction that occurs inside a soda ash carbonation tower, which includes multiple trays. The method for solving the carbonation reaction also includes:

[0082] Based on the preset heat balance equation, the heat data of a single tray is calculated, and the plate-by-plate calculation is performed according to the preset overall tower balance equation to obtain the temperature distribution data, carbon dioxide absorption rate distribution data, crystal particle size distribution data, and supersaturation distribution data of each tray in the soda ash carbonization tower.

[0083] Specifically, the preset heat balance equation can be a heat balance equation determined after considering the heat released by the carbonization tower reaction, the heat of hydrolysis, the heat of crystallization, the heat transfer due to cooling, the heat dissipation from the environment, the sensible heat of the liquid phase, and the sensible heat of the gas phase. The relevant expressions may include:

[0084]

[0085] By using the heat balance equation to perform heat balance on each plate of the carbonization tower, the temperature increment after obtaining the exothermic reaction and heat exchange can be obtained, which can be used for plate-by-plate heat calculation.

[0086] In some specific embodiments, the total tower balance equation can be determined based on the component pressure balance equation, the moment balance of crystal grain size distribution, and the component mass balance equation.

[0087] The expressions related to the total tower balance equation can include:

[0088] ΔP=ΔP ad +10 -5 ρhg(1-ε)

[0089] Δμ (0) =Bτ,Δμ (k) =kGτμ (k-1)

[0090] qΔc i =N i aV+r i V(1-ε)

[0091]

[0092] By calculating the temperature distribution, carbon dioxide absorption rate distribution, crystal size distribution, and supersaturation distribution of each tray in the soda ash carbonization tower using the overall tower balance equation, the relationship between the carbonization reaction operating conditions and the final reaction result can be established, providing a reference for the parameter adjustment of the control equipment.

[0093] In some embodiments, determining the crystallization reaction rate and crystal grain size distribution data based on a preset crystallization reaction equation and crystallization reaction parameters includes:

[0094] Based on the preset supersaturation calculation equation, the crystal growth rate and nucleation rate are determined;

[0095] Based on the crystallization reaction equation, crystal growth rate, nucleation rate, and crystallization reaction parameters, the crystallization reaction rate and crystal size distribution data are determined.

[0096] In some embodiments, the method further includes: determining the kinetic parameters and equilibrium constant of the carbonization reaction based on the reaction solving equation.

[0097] In some embodiments, the carbonization reaction solution method can be applied to the carbonization tower control system to solve for key reaction data during the carbonization process, thereby controlling relevant reaction parameters. Key reaction data include, for example, absorption rate, reaction rate, crystallite size distribution data, and concentration values ​​of reactant components.

[0098] In some embodiments, the carbonization reaction solution method can be applied to a monitoring system for the carbonization tower reaction process, which includes a data acquisition unit, a server, and a display unit.

[0099] Specifically, the data acquisition unit acquires the input parameter data involved in solving the equations in the solution method. These input parameters can come from real-time data directly collected by sensors, experimentally calculated data, or data manually entered by staff based on experience. The input parameter data acquired by the data acquisition unit is transmitted to the server for solution calculation, generating monitoring results which are then displayed on the display unit. The key reaction data obtained by the server in the carbonization process can serve as the basis for adjusting the parameters of the carbonization tower reaction process.

[0100] In some cases, the monitoring system can also communicate with the carbonization tower control system. The server in the monitoring system can issue corresponding operation commands based on the solution results, allowing the carbonization tower control equipment to adjust process equipment or components according to these commands, thereby optimizing the carbonization production process. In other cases, the server can also display the key reaction data obtained from the solution in real time on the monitoring system's display interface. The display method is not limited to text, charts, or animations, to prompt staff to adjust the corresponding process equipment parameters.

[0101] In some specific embodiments, the code corresponding to the equations involved in the carbonization reaction solution method can be written in MATLAB (Matrix Laboratory, a high-level programming language and numerical computing environment).

[0102] It should be understood that, although Figure 1 The steps in the flowchart are shown sequentially according to the arrows, but these steps are not necessarily executed in the order indicated by the arrows. Figure 1Unless otherwise expressly stated herein, the steps illustrated and other steps involved in the embodiments are not subject to strict order restrictions and may be performed in other orders. Furthermore, at least some steps in the foregoing embodiments may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily completed at the same time, but may be performed at different times. The execution order of these sub-steps or stages is not necessarily sequential, but may be performed alternately or in turn with other steps or at least a portion of the sub-steps or stages of other steps.

[0103] In a second aspect, embodiments of this disclosure provide a solver whose code is written in MATLAB, the solver being used to implement the steps of the carbonization reaction solving method provided in any embodiment of the first aspect of this disclosure.

[0104] In some specific embodiments, the solver code is written using MATLAB, and the code file includes...

[0105] In some embodiments, the solver includes a carbonation reaction kinetics module, an absorption and mass transfer module, a crystallization process simulation module, a component concentration calculation module, a carbonation reaction module, a heat balance module, and a plate-by-plate balance module.

[0106] The carbonation reaction kinetics module is used to call reaction equations to solve for the absorption reaction rate and the hydrolysis reaction rate. It is also used to solve for kinetic parameters and equilibrium constants for each reaction.

[0107] The absorption and mass transfer module is used to call the absorption mass transfer equation to calculate the carbon dioxide absorption rate and mass transfer rate, as well as to determine the mass transfer and reaction processes of various substances within the membrane.

[0108] The crystallization process simulation module is used to call the crystallization reaction equation to calculate the crystal grain size distribution and crystallization reaction rate, as well as to determine the expression for the crystallization reaction rate.

[0109] The component concentration calculation module is used to call the concentration solving equation to calculate the concentration values ​​of reactant components after carbonization.

[0110] The carbonation reaction module is used to cyclically call the carbonation reaction kinetics module, absorption and mass transfer module, crystallization process simulation module, and component concentration calculation module to solve the problem until the carbon dioxide absorption rate and mass transfer rate converge and the entire carbonation process is simulated.

[0111] The heat balance module is used to call the heat balance equation to perform heat balance calculations on the carbonization tower plates.

[0112] The plate-by-plate balance module is used to call the overall tower balance equation to calculate the distribution of temperature, component concentration, and CO2 absorption rate for each plate of the carbonization tower.

[0113] For further specific limitations regarding the solver, please refer to the limitations of the carbonization reaction solution method above. The solver can also be used to perform more steps of the carbonization reaction solution method in the embodiments of this disclosure, which will not be repeated here. Each module in the solver described above can be implemented entirely or partially through software, hardware, or a combination thereof. These modules can be embedded in or independent of the processor in a computer device in hardware form, or stored in the memory of a computer device in software form, so that the processor can call and execute the operations corresponding to each module.

[0114] In a third aspect, embodiments of this disclosure provide a computer device including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the steps of the carbonization reaction solving method provided in any embodiment of the first aspect of this disclosure.

[0115] In some embodiments, the computer device may be a server, and its internal structure diagram may be as follows: Figure 2 As shown. The computer device includes a processor, memory, network interface, and database connected via a system bus. The processor provides computing and control capabilities. The memory includes a non-volatile storage medium and internal memory. The non-volatile storage medium stores an operating system, computer programs, and a database. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage medium. The database stores code files. The network interface communicates with external terminals via a network connection. When the computer program is executed by the processor, it implements the carbonization reaction solution method in any embodiment of this document.

[0116] Those skilled in the art will understand that Figure 2 The structures shown are merely block diagrams of some structures related to the embodiments of this disclosure and do not constitute a limitation on the computer devices to which the embodiments of this disclosure are applied. Specific computer devices may include more or fewer components than those shown in the figures, or combine certain components, or have different component arrangements.

[0117] In a third aspect, embodiments of this disclosure provide a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the carbonization reaction solution method provided in any embodiment of the first aspect of this disclosure.

[0118] The computer-readable storage medium may be Figure 2The computer-readable storage medium in the computer device shown.

[0119] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The aforementioned computer program can be stored in a non-volatile computer-readable storage medium. When executed, the computer program can include the processes of the embodiments of the above methods. Any references to memory, storage, databases, or other media used in the embodiments of this disclosure can include non-volatile and / or volatile memory. Non-volatile memory may include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory may include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in a variety of forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), dual data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.

[0120] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this disclosure.

[0121] The above embodiments merely illustrate several implementation methods of this disclosure, and their descriptions are relatively specific and detailed, but they should not be construed as limiting the scope of protection of this disclosure. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of this disclosure, and these all fall within the scope of protection of this disclosure. Therefore, the scope of protection of this disclosure should be determined by the appended claims.

Claims

1. A method for solving carbonization reactions, characterized in that, The method includes: Based on the preset reaction equations and the initial concentrations, temperatures, reaction rate constants, and ionic strengths of the reactant components, the absorption reaction rate and the hydrolysis reaction rate are determined. Based on the preset absorption mass transfer equation and the absorption reaction rate, determine the carbon dioxide absorption rate and mass transfer rate corresponding to the concentration values ​​of reactant components obtained during the carbonization reaction. Based on the preset crystallization reaction equation and crystallization reaction parameters, determine the crystallization reaction rate and crystal grain size distribution data; Based on the preset concentration calculation equation, the crystallization reaction rate, the hydrolysis reaction rate, and the reactant component parameters, the concentration values ​​of the reactant components after the carbonization reaction are determined. During the carbonization reaction, the reaction solution equation, the absorption mass transfer equation, the crystallization reaction equation, and the concentration solution equation are called in a loop until the carbon dioxide absorption rate and the mass transfer rate converge. The latest reactant component concentration values, the latest crystallization reaction rate, and the latest crystal particle size distribution data are output to provide the corresponding control equipment for parameter adjustment of the carbonization reaction.

2. The method according to claim 1, characterized in that, The carbonation reaction is a carbonation reaction that occurs inside a soda ash carbonation tower, which includes multiple trays. The method further includes: According to the preset heat balance equation, the heat data of a single tray is calculated, and according to the preset whole-tower balance equation, the tray-by-tray calculation is performed to obtain the temperature distribution data, carbon dioxide absorption rate distribution data, crystal particle size distribution data and supersaturation distribution data of each tray in the soda ash carbonization tower.

3. The method according to claim 2, characterized in that, The total balance equations include the pressure balance equations for material components, the moment balance equations for crystal grain size distribution, and the mass balance equations for material components.

4. The method according to claim 1, characterized in that, The step of determining the crystallization reaction rate and crystal grain size distribution data based on the preset crystallization reaction equation and crystallization reaction parameters includes: Based on the preset supersaturation calculation equation, the crystal growth rate and nucleation rate are determined; Based on the crystallization reaction equation, the crystal growth rate, the nucleation rate, and the crystallization reaction parameters, the crystallization reaction rate and crystal size distribution data are determined.

5. The method according to claim 1, characterized in that, The method further includes: determining the kinetic parameters and equilibrium constant of the carbonization reaction by solving the reaction equation.

6. A solver, characterized in that, The solver code is written in MATLAB, and the solver is used to implement the steps of the method according to any one of claims 1 to 5.

7. A computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the steps of the method according to any one of claims 1 to 5.

8. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 5.

Images