Method and system for optimizing coal gasification reaction rate correction coefficient

Abstract

The embodiment of the invention provides an optimization method and system for a coal gasification reaction rate correction coefficient. The optimization method comprises the following steps: acquiring the initial concentration of reaction components in a coal gasification reaction of a reactor to be treated; based on the initial concentration, through a kinetic parameter equation of the coal gasification reaction, the final concentration of reaction components in the coal gasification reaction of the to-be-treated reactor after a preset time t is obtained; and constructing an objective function based on the final concentration of the to-be-corrected component in the final concentrations and the actual concentration of the to-be-corrected component, and optimizing the reaction rate correction coefficient of the coal gasification reaction to obtain an optimal reaction rate correction coefficient.

Description

Technical Field

[0001] The embodiments of the present application belong to the technical field of reactor dynamics simulation, and in particular, to a method and system for optimizing a coal gasification reaction rate correction coefficient. Background Art

[0002] The simulation of coal gasification process can help engineers optimize the design and operation parameters, thus improving the economy and environmental friendliness of coal gasification. However, in situations where coal gasification reaction needs to be simulated, due to factors such as furnace size and coal properties, there will be a certain gap between the simulated data and the actual production data. How to quickly adjust the kinetic parameters to make the simulation data fit the actual production, so as to accurately predict the reaction product situation is an important task.

[0003] In the prior art, the optimization of the coal gasification reaction rate correction coefficient is mainly performed manually. This method of manually correcting the coal gasification reaction rate correction system is time-consuming and labor-intensive, and cannot guarantee accuracy and cannot be automated.

[0004] Based on this, a new optimization method and system for the coal gasification reaction rate correction coefficient is needed. Summary of the invention

[0005] The embodiments of this specification provide a method and system for optimizing a coal gasification reaction rate correction coefficient to solve some or all of the problems: In the prior art, the optimization of the coal gasification reaction rate correction coefficient is mainly performed manually. This method of manually calibrating the coal gasification reaction rate correction system is time-consuming and labor-intensive, and cannot guarantee accuracy and cannot be automated.

[0006] To solve the above technical problems, the embodiments of this specification are implemented as follows:

[0007] The embodiment of this specification provides a method for optimizing a coal gasification reaction rate correction coefficient, and the optimization method includes:

[0008] Obtaining the initial concentration of reaction components in the coal gasification reaction of the reactor to be treated;

[0009] Based on the initial concentration, the final concentration of the reaction component in the coal gasification reaction of the reactor to be treated after a preset time t is obtained by using a kinetic parameter equation of the coal gasification reaction;

[0010] An objective function is constructed based on the final concentration of the component to be corrected in the final concentration and the actual concentration of the component to be corrected, and an optimal reaction rate correction coefficient is obtained by optimizing the reaction rate correction coefficient of the coal gasification reaction.

[0011] The present specification also provides an optimization system for a coal gasification reaction rate correction coefficient, the optimization system comprising:

[0012] An initial concentration acquisition module is used to acquire the initial concentration of the reaction components in the coal gasification reaction of the reactor to be processed;

[0013] A final concentration calculation module, based on the initial concentration, obtains the final concentration of the reaction component in the coal gasification reaction of the reactor to be treated after a preset time t through a kinetic parameter equation of the coal gasification reaction;

[0014] The optimization module constructs an objective function based on the final concentration of the component to be corrected in the final concentration and the actual concentration of the component to be corrected, and obtains the optimal reaction rate correction coefficient by optimizing the reaction rate correction coefficient of the coal gasification reaction.

[0015] At least one of the above-mentioned technical solutions adopted in the embodiments of this specification can achieve the following beneficial effects: by obtaining the initial concentration of the reaction components in the coal gasification reaction of the reactor to be treated; based on the initial concentration, the final concentration of the reaction components in the coal gasification reaction of the reactor to be treated after a preset time t is obtained through the kinetic parameter equation of the coal gasification reaction; constructing an objective function based on the final concentration of the component to be corrected in the final concentration and the actual concentration of the component to be corrected, and optimizing the reaction rate correction coefficient of the coal gasification reaction to obtain the optimal reaction rate correction coefficient, which can realize automatic optimization of the coal gasification reaction rate correction coefficient, and has a sensitive response, and can improve the accuracy of the coal gasification reaction simulation results. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] The drawings described herein are used to provide a further understanding of the present application and constitute a part of the present application. The illustrative embodiments of the present application and their descriptions are used to explain the present application and do not constitute an improper limitation on the present application. Some specific embodiments of the present application will be described in detail in an illustrative and non-restrictive manner with reference to the drawings. The same reference numerals in the drawings indicate the same or similar components or parts. It should be understood by those skilled in the art that these drawings are not necessarily drawn to scale. In the drawings:

[0017] Figure 1 A schematic flow chart of a method for optimizing a coal gasification reaction rate correction coefficient provided in an embodiment of this specification;

[0018] Figure 2 A flow chart for obtaining the final concentration provided in the embodiments of this specification;

[0019] Figure 3 A schematic diagram of a system for optimizing a coal gasification reaction rate correction coefficient provided in an embodiment of this specification. DETAILED DESCRIPTION

[0020] In order to enable those skilled in the art to better understand the present application, the technical solutions in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are only embodiments of a part of the present application, not all embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in the field without creative work should fall within the scope of protection of the present application.

[0021] Figure 1 The following is a flow chart of a method for optimizing the correction coefficient of the coal gasification reaction rate provided in the embodiment of this specification. From the program perspective, the execution subject of the method for optimizing the correction coefficient of the coal gasification reaction rate is a program installed on an application server or an application terminal. It can be understood that the method can be executed by any device, equipment, platform, or equipment cluster with computing and processing capabilities. Figure 1 As shown, the optimization method includes:

[0022] Step S101: obtaining the initial concentration of the reaction components in the coal gasification reaction of the reactor to be treated;

[0023] Pyrolysis of coal gasification is an important step in the coal gasification process, which involves decomposing solid fuel (coal) into gas, liquid and solid products under high temperature and oxygen-deficient conditions. Pyrolysis usually occurs in the drying and pyrolysis zone of the coal gasifier, with temperatures ranging from approximately 300°C to 600°C.

[0024] Coal gasification pyrolysis considers coal pyrolysis into volatile matter and coke, wherein the coke contains fixed carbon, ash and residual volatile matter, the volatile matter will be partially burned, and the partial combustion products and coke will be further gasified.

[0025] In the embodiments of this specification, the reactor to be treated is a reactor for performing coal gasification reaction, and the specific structure and type of the reactor to be treated does not constitute a limitation to the present application.

[0026] During the coal gasification process, the reaction of combustion of part of the volatile matter is:

[0027] C6H6+7.5O2→6CO2+3H2O

[0028] H2+0.5O2→H2O

[0029] CO+0.5O2→CO2

[0030] CH4+2O2→CO2+2H2O

[0031] In the coal gasification reaction, the gasification reaction is divided into homogeneous reaction and heterogeneous reaction. The reactions involved in the gasification reaction are shown in Table 1.

[0032] Table 1 Gasification reaction equation of coal gasification reaction

[0033] reaction Gas phase reaction (gas) reaction Heterogeneous reaction (gas-solid) Reaction 1 <![CDATA[CO+0.5O2→CO2]]> Reaction 7 <![CDATA[C+H2O→CO+H2]]> Reaction 2 <![CDATA[H2+0.5O2→H2O]]> Reaction 8 <![CDATA[C+2H2→CH4]]> Reaction 3 <![CDATA[CO+H2O→CO2+H2]]> Reaction 9 <![CDATA[C+1 / ΦO2→(2-2 / Φ)CO+(2 / Φ-1)CO2]]> Reaction 4 <![CDATA[CO2+H2→CO+H2O]]> Reaction 10 <![CDATA[C+CO2→2CO]]> React 5 <![CDATA[CH4+2O2→CO2+2H2O]]> Reaction 11 <![CDATA[CO+H2O→CO2+H2]]> Reaction 6 <![CDATA[CH4+H2O→CO+3H2]]>

[0034] It should be particularly noted that, in the gasification reaction equation of the coal gasification reaction, the reaction equations of reaction 3 and reaction 11 are similar, but reaction 11 contains the case where Fe2O3 is present in the coal ash.

[0035] According to the gasification reaction equation of the coal gasification reaction, in the embodiments of this specification, the reaction components in the coal gasification reaction in the reactor to be treated include: oxygen, methane, carbon monoxide, carbon dioxide, hydrogen sulfide, hydrogen, nitrogen, and water.

[0036] Step S103: Based on the initial concentration, the final concentration of the reaction components in the coal gasification reaction of the reactor to be treated after a preset time t is obtained by using a kinetic parameter equation of the coal gasification reaction.

[0037] In the examples of this specification, the kinetic equation of the coal gasification reaction is:

[0038]

[0039] R H2S =R VH2S,i

[0040]

[0041] R N2 =R vN2,i

[0042]

[0043] in,

[0044] R O2 Indicates the rate at which oxygen is produced or consumed;

[0045] R CH4 Indicates the rate of methane production or consumption;

[0046] R CO Indicates the rate of production or consumption of carbon monoxide;

[0047] R CO2 Indicates the rate of production or consumption of carbon dioxide;

[0048] R H2S Indicates the rate of generation or consumption of hydrogen sulfide;

[0049] R H2 Indicates the rate of generation or consumption of hydrogen;

[0050] R N2 Indicates the rate of generation or consumption of nitrogen;

[0051] R H2 Indicates the rate at which water is produced or consumed;

[0052] r1 represents the reaction rate of reaction 1;

[0053] r2 represents the reaction rate of reaction 2;

[0054] r3 represents the reaction rate of reaction 3;

[0055] r4 represents the reaction rate of reaction 4;

[0056] r5 represents the reaction rate of reaction 5;

[0057] r6 represents the reaction rate of reaction 6;

[0058] r7 represents the reaction rate of reaction 7;

[0059] r8 represents the reaction rate of reaction 8;

[0060] r9 represents the reaction rate of reaction 9;

[0061] r 10 represents the reaction rate of reaction 10;

[0062] r 11 represents the reaction rate of reaction 11;

[0063] V represents the volume of the reactor to be treated;

[0064] Φ represents the mechanical factor;

[0065] N P Indicates the total number of coke particles in the reactor to be treated;

[0066] d p represents the diameter of coke particles;

[0067] R VO2,i It indicates the rate at which coke releases oxygen;

[0068] R VCH4,i It indicates the rate of methane release from coke;

[0069] R VCO,i Indicates the rate at which coke releases carbon monoxide;

[0070] R VCO2,i Indicates the rate at which coke releases carbon dioxide;

[0071] R VH2S,i Indicates the rate at which coke releases hydrogen sulfide;

[0072] R VH2,i Indicates the rate at which coke releases hydrogen;

[0073] R VH2,i It indicates the rate of nitrogen release from coke;

[0074] R VH2O,i It indicates the rate at which coke releases water;

[0075] π stands for pi.

[0076] It should be noted that regarding the mechanical factor Φ, when Φ=1, CO2 is the only product, and when Φ=2, CO is the only product. Therefore, Φ is usually between 1 and 2.

[0077] In the examples of this specification, r1, r2, r3, r4 and r5 are all homogeneous reactions;

[0078] The expressions of r1, r2, r3, r4 and r5 are:

[0079]

[0080] in,

[0081] k0 represents the pre-exponential factor;

[0082] Cor k0 represents the correction coefficient for k0;

[0083] E represents the activation energy of chemical reaction kinetics;

[0084] Cor E represents the correction factor for E;

[0085] R represents a constant;

[0086] T g Indicates the reaction gas temperature;

[0087] C A , C B Represents r 1-5 The concentrations of the corresponding reactants;

[0088] r 1-5 Represents reaction 1, reaction 2, reaction 3, reaction 4, reaction 5; the expression of r6 is:

[0089]

[0090] in,

[0091] represents the concentration of reactant methane;

[0092] R7, R8, R9, R10 and r 11 All are heterogeneous reactions;

[0093] R7, R8, R9, R 10 The expression is:

[0094]

[0095] in

[0096] k diff It indicates the rate at which gas diffuses to the coke surface;

[0097] Cor kd Indicates k diff The correction factor of

[0098] r 7-10 It indicates reaction 7, reaction 8, reaction 9, and reaction 10;

[0099] Y represents the ratio of the unreacted carbon core diameter to the coke diameter;

[0100] k dash Indicates the rate at which gas diffuses through the outer ash layer to the coke surface; k s represents the reaction rate on the surface of the carbon core;

[0101]

[0102] T p Indicates the surface temperature of coke particles;

[0103] d p represents the diameter of coke particles;

[0104] Indicates the effective partial pressure corresponding to reactions 7-10

[0105] The r 11 The expression is:

[0106]

[0107] in,

[0108] r 11 Indicates reaction 11;

[0109] F w represents the catalyst activity correction factor;

[0110] P f represents the pressure correction factor;

[0111] T f represents the temperature correction factor;

[0112] represents the effective partial pressure corresponding to reaction 11.

[0113] In the examples of this specification, the catalyst activity correction factor F w The value of is 0.2.

[0114] Although the reaction equations of reaction 3 and reaction 11 are similar, reaction 11 involves the situation where Fe2O3 is present in the coal ash, and therefore, the kinetic parameter equations of reaction 3 and reaction 11 are different.

[0115] In the examples of this specification, the expression for the effective partial pressure corresponding to the reaction 7 is:

[0116]

[0117] in,

[0118] represents the effective partial pressure corresponding to reaction 7;

[0119] Indicates the pressure of water in the coal gasifier to be treated;

[0120] Indicates the pressure of hydrogen in the coal gasifier to be treated;

[0121] P CO Indicates the pressure of carbon monoxide in the coal gasifier to be treated;

[0122] Cor a represents the correction coefficient for a;

[0123] Cor b represents the correction coefficient for b;

[0124] a represents the equilibrium coefficient of the first reaction;

[0125] b represents the equilibrium coefficient of the second reaction;

[0126] c represents the equilibrium coefficient of the third reaction;

[0127] T p Indicates the surface temperature of coke particles;

[0128] The expression of the effective partial pressure corresponding to the reaction 8 is:

[0129]

[0130] in,

[0131] represents the effective partial pressure corresponding to reaction 8;

[0132] Indicates the pressure of methane in the coal gasifier to be treated;

[0133] The expression of the effective partial pressure corresponding to the reaction 9 is:

[0134]

[0135] in,

[0136] represents the effective partial pressure corresponding to reaction 9;

[0137] Indicates the pressure of oxygen in the coal gasifier to be treated;

[0138] The expression of the effective partial pressure corresponding to the reaction 10 is:

[0139]

[0140] in,

[0141] represents the effective partial pressure corresponding to reaction 10;

[0142] Indicates the pressure of carbon dioxide in the coal gasifier to be treated;

[0143] The expression of the effective partial pressure corresponding to the reaction 11 is:

[0144]

[0145] in,

[0146] represents the effective partial pressure corresponding to reaction 11.

[0147] It should be noted that the first reaction equilibrium coefficient a, the second reaction equilibrium coefficient b, and the third reaction equilibrium coefficient c are all constants, and the first reaction equilibrium coefficient a, the second reaction equilibrium coefficient b, and the third reaction equilibrium coefficient c are related to the effective partial pressure in the heterogeneous reaction, and in different reactions, the first reaction equilibrium coefficient a, the second reaction equilibrium coefficient b, and the third reaction equilibrium coefficient c have different values. For example, in reaction 7, reaction 8, and reaction 10, the first reaction equilibrium coefficient a, the second reaction equilibrium coefficient b, and the third reaction equilibrium coefficient c have different values. For the effective partial pressures of different components, the first reaction equilibrium coefficient a, the second reaction equilibrium coefficient b, and the third reaction equilibrium coefficient c also have different values.

[0148] The kinetic equation of the coal gasification reaction is used to represent the reaction generation or consumption rate of each component of the coal gasification reaction. Therefore, based on the kinetic equation, the final reactant concentration can be calculated.

[0149] In the embodiment of this specification, based on the initial concentration, the final concentration of the reaction components in the coal gasification reaction of the reactor to be treated after the preset time t is obtained by the kinetic parameter equation of the coal gasification reaction, specifically including:

[0150] Based on the initial concentration and the kinetic parameter equation of the coal gasification reaction, the preset time t is iterated according to dt to calculate the concentration of each reaction component in the coal gasification reaction of the reactor to be treated, and the calculated value of each reaction component in the coal gasification reaction of the reactor to be treated when the final dt≥the preset time t is used as the final concentration of the reaction component in the coal gasification reaction of the reactor to be treated after the preset time t.

[0151] In the embodiment of this specification, the kinetic parameter equation based on the initial concentration and the coal gasification reaction, iterating the preset time t according to dt, and calculating the concentration of each reaction component in the coal gasification reaction of the reactor to be treated, specifically includes:

[0152] Based on the initial concentration and the kinetic parameter equation of the coal gasification reaction, a calculation formula for obtaining the concentration of each reaction component in the coal gasification reaction of the reactor to be treated;

[0153] Based on the calculation formula of the concentration of each reaction component in the coal gasification reaction of the reactor to be treated, the concentration of each reaction component in the coal gasification reaction of the reactor to be treated is iterated according to dt for the preset time t;

[0154] The calculation formula for the concentration of each reaction component in the coal gasification reaction of the reactor to be treated is:

[0155]

[0156] in,

[0157] R i An expression representing the rate at which component i is generated or consumed by a reaction;

[0158] is the initial concentration of component i;

[0159] c i is the concentration of component i after the preset time t.

[0160] In order to further understand the generation of the final concentration, it will be explained with reference to a specific example below.

[0161] Figure 2 The final concentration acquisition flow chart provided in the embodiments of this specification. Figure 2 As shown, within the preset time t of the reactor to be treated, the final concentration can be calculated by the coal gasification reaction kinetic parameter equation. The iterative calculation process is to assume a very small time dt, calculate the real-time concentration of each component according to the reaction kinetic parameter equation, then increase dt, repeat the above calculation, and finally make the calculated reaction time greater than or equal to t, at which time the concentration of each component is the final reaction concentration.

[0162] Step S105: constructing an objective function based on the final concentration of the component to be corrected in the final concentration and the actual concentration of the component to be corrected, and obtaining an optimal reaction rate correction coefficient by optimizing the reaction rate correction coefficient of the coal gasification reaction.

[0163] According to the kinetic equation of coal gasification reaction, there is an adjustable reaction rate correction factor in the coal gasification reaction. Specifically, the reaction rate correction factor includes: the correction factor Cor of the rate at which gas diffuses to the coke surface kd , Correction coefficient Cor of pre-exponential factor of chemical reaction kinetic rate k0 , Correction coefficient Cor of chemical reaction kinetics activation energy E , the correction coefficient Cor for a a , the correction factor Cor for b b , and the correction factor Cor for the rate at which gas diffuses to the coke surface kd , Correction coefficient Cor of pre-exponential factor of chemical reaction kinetic rate k0 , Correction coefficient Cor of chemical reaction kinetics activation energy E , the correction coefficient Cor for a a , the correction factor Cor for b b The default value of both is 1.

[0164] The relationship between the reaction equation and the reaction rate correction coefficient is shown in Table 2.

[0165] Table 2: Adjustable kinetic parameters in coal gasification reaction

[0166]

[0167] It can be seen that since the reaction rate correction coefficient can be adjusted, the concentration of each reaction component will be adjusted accordingly as the reaction rate correction coefficient is adjusted, so that the concentration of the reaction component obtained by simulation calculation is closer and closer to the actual concentration. In other words, by adjusting the reaction rate correction coefficient, the optimal reaction rate correction coefficient can be obtained, which can achieve a more accurate concentration of the simulated reaction component.

[0168] Since certain components are often of interest in actual coal gasification reactions, in the embodiments of this specification, components to be corrected may be selected according to business needs, and then reaction rate correction coefficients to be optimized may be determined based on the reactions involved in the components to be corrected.

[0169] The component to be corrected is a component that needs to be corrected, and is generally a component that the user is interested in. In a specific embodiment, the component to be corrected is one or more components in the coal gasification reaction components, that is, the component to be corrected is one or more of oxygen, methane, carbon monoxide, carbon dioxide, hydrogen sulfide, hydrogen, nitrogen, and water.

[0170] In the embodiments of this specification, optimization is performed using the BFGSBMinimizer algorithm.

[0171] In the embodiment of this specification, the objective function is constructed based on the final concentration of the component to be corrected in the final concentration and the actual concentration of the component to be corrected, and the optimal reaction rate correction coefficient is obtained by optimizing the reaction rate correction coefficient of the coal gasification reaction, specifically including:

[0172] Based on the final concentrations, respectively determining the final concentrations of the components to be calibrated;

[0173] Based on the actual concentrations of oxygen, methane, carbon monoxide, carbon dioxide, hydrogen sulfide, hydrogen, nitrogen, and water in the reactor to be treated, respectively determining the actual concentrations of the components to be corrected;

[0174] constructing an objective function based on the actual concentration of the component to be corrected and the final concentration of the component to be corrected;

[0175] The reaction rate correction coefficient of coal gasification reaction is optimized by BFGSBMinimizer algorithm to obtain the optimal reaction rate correction coefficient.

[0176] In the embodiment of this specification, the expression of the objective function is:

[0177]

[0178] in,

[0179] c i is the actual concentration of component i to be corrected.

[0180] c c,i is the final concentration of component i to be calibrated;

[0181] N represents the number of groups to be corrected;

[0182] f(x) represents the objective function.

[0183] BFGSBMinimizer is an algorithm for numerical optimization, belonging to the Quasi-Newton method. It is a variant of the Broyden–Fletcher–Goldfarb–Shanno (BFGS) algorithm, used to solve constrained optimization problems. The BFGS algorithm is an iterative method for finding the local minimum of a multivariable function. It does not require the calculation of the second-order derivative of the objective function (i.e., the Hessian matrix), but approximates it through information from the iterative process.

[0184] In the embodiments of this specification, the reaction rate correction coefficient of the coal gasification reaction is optimized by the BFGSBMinimizer algorithm to obtain the optimal reaction rate correction coefficient, which specifically includes:

[0185] The reaction rate correction coefficient of the coal gasification reaction is adjusted by the BFGSBMinimizer algorithm so that the objective function is minimized and the optimal reaction rate correction coefficient is obtained.

[0186] In a specific embodiment, when the objective function is less than or equal to a preset tolerance or the number of iterations reaches a preset number of iterations, the objective function reaches a minimum value and the iteration stops. The reaction rate correction coefficient at this time is the optimal reaction rate correction coefficient.

[0187] In one embodiment, the preset tolerance is 0.001, and the preset number of iterations is 10. It should be noted that the preset tolerance and / or the preset number of iterations may be determined according to specific scenarios.

[0188] In the examples of this specification, it is illustrated that multiple reaction rate correction coefficients are optimized. In specific embodiments, one or more reaction rate correction coefficients can be optimized. At the same time, in the process of optimizing the reaction rate correction coefficients, reactions that need to be calibrated and / or components that need to be calibrated can be selected.

[0189] The examples of this specification further verify the optimization effect of the method for optimizing the coal gasification reaction rate correction coefficient provided by this application.

[0190] In one embodiment of the present specification, the coal gasification process of the treatment reactor (coal gasification reactor) is simulated, and the feed and product components required by the device are input in the coal gasification reaction. Since the concentrations of carbon monoxide CO and hydrogen H2 in the coal gasification simulation software used are quite different from the actual concentrations, and the water gas shift reaction has a relatively large impact on the concentrations of CO and H2 (CO+H2O→CO2+H2), in this embodiment, a correction water gas shift is used, and a carbon dioxide target is set to determine the reaction rate correction coefficient that needs to be corrected, that is, the correction coefficient Cor for the pre-exponential factor of the chemical reaction kinetic ratek0 , Correction coefficient Cor of chemical reaction kinetics activation energy E , the correction coefficient Cor for a a , the correction factor Cor for b b Before the correction, each reaction rate correction coefficient was 1, but after the correction by the method provided in the examples of this specification, the reaction rate correction coefficient changed, as shown in Table 3 for details.

[0191] Table 3: Optimization results of reaction rate correction coefficient

[0192] Reaction rate correction factor Before calibration After correction <![CDATA[Cor k0 ]]> 1 1.005 <![CDATA[Cor E ]]> 1 0.96424 <![CDATA[Cor a ]]> 1 0.77750 <![CDATA[Cor b ]]> 1 1.134

[0193] Furthermore, the concentrations of the reaction components before and after the correction of the reaction rate correction coefficient were compared, as shown in Table 4 for details.

[0194] Table 4: Comparison of reactant concentrations before and after reaction rate correction factor correction

[0195] Concentration before correction Corrected concentration Actual concentration CO 0.54364 0.43769 0.44257 H2 0.30832 0.35574 0.36349

[0196] According to Table 4, after the reaction rate correction coefficient is corrected, the concentration of the reactant is closer to the actual concentration, which improves the accuracy of the simulation results.

[0197] The method for optimizing the coal gasification reaction rate correction coefficient provided in the embodiments of this specification can realize automatic optimization of the coal gasification reaction rate correction coefficient, and has a sensitive response, and can improve the accuracy of the coal gasification reaction simulation results.

[0198] The aforementioned embodiment of this specification provides a method for optimizing a coal gasification reaction rate correction coefficient. Based on the same idea, the embodiment of this specification also provides a system for optimizing a coal gasification reaction rate correction coefficient. Figure 3 A schematic diagram of a system for optimizing a coal gasification reaction rate correction coefficient provided in an embodiment of this specification. Figure 3 As shown, the optimization system includes:

[0199] An initial concentration acquisition module 301 is used to acquire the initial concentration of the reaction components in the coal gasification reaction of the reactor to be processed;

[0200] The final concentration calculation module 303 obtains the final concentration of the reaction component in the coal gasification reaction of the reactor to be treated after a preset time t based on the initial concentration and through the kinetic parameter equation of the coal gasification reaction;

[0201] The optimization module 305 constructs an objective function based on the final concentration of the component to be corrected in the final concentration and the actual concentration of the component to be corrected, and obtains the optimal reaction rate correction coefficient by optimizing the reaction rate correction coefficient of the coal gasification reaction.

[0202] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present application, rather than to limit it. Although the present application has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or replace some or all of the technical features therein with equivalents. However, these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present application.

Claims

1. A method for optimizing a coal gasification reaction rate correction coefficient, characterized in that: The optimization method comprises: Obtaining the initial concentration of reaction components in the coal gasification reaction of the reactor to be treated; Based on the initial concentration, the final concentration of the component in the coal gasification reaction of the reactor to be treated after a preset time t is obtained by using a kinetic parameter equation of the coal gasification reaction; An objective function is constructed based on the final concentration of the component to be corrected in the final concentration and the actual concentration of the component to be corrected, and an optimal reaction rate correction coefficient is obtained by optimizing the reaction rate correction coefficient of the coal gasification reaction.

2. The optimization method according to claim 1, characterized in that: The reaction components in the coal gasification reaction in the reactor to be treated include: oxygen, methane, carbon monoxide, carbon dioxide, hydrogen sulfide, hydrogen, nitrogen and water.

3. The optimization method according to claim 1, characterized in that: The kinetic equation of the coal gasification reaction is: R H2S =R VH2S,i R N2 =R VN2,i in, R O2 Indicates the rate at which oxygen is produced or consumed; R CH4 Indicates the rate of methane production or consumption; R CO Indicates the rate of production or consumption of carbon monoxide; R CO2 Indicates the rate of production or consumption of carbon dioxide; R H2S Indicates the rate of generation or consumption of hydrogen sulfide; R H2 Indicates the rate of generation or consumption of hydrogen; R N2 Indicates the rate of generation or consumption of nitrogen; R H2 Indicates the rate at which water is produced or consumed; r1 represents the reaction rate of reaction 1; r2 represents the reaction rate of reaction 2; r3 represents the reaction rate of reaction 3; r4 represents the reaction rate of reaction 4; r5 represents the reaction rate of reaction 5; r6 represents the reaction rate of reaction 6; r7 represents the reaction rate of reaction 7; r8 represents the reaction rate of reaction 8; r9 represents the reaction rate of reaction 9; r 10 represents the reaction rate of reaction 10; r 11 represents the reaction rate of reaction 11; V represents the volume of the reactor to be treated; Φ represents the mechanical factor; N P Indicates the total number of coke particles in the reactor to be treated; d p represents the diameter of coke particles; R VO2,i It indicates the rate at which coke releases oxygen; R VCH4,i It indicates the rate of methane release from coke; R VCO,i Indicates the rate at which coke releases carbon monoxide; R VCO2,i Indicates the rate at which coke releases carbon dioxide; R VH2S,i Indicates the rate at which coke releases hydrogen sulfide; R VH2,I Indicates the rate at which coke releases hydrogen; R VH2,i It indicates the rate of nitrogen release from coke; R VH2O,i It indicates the rate at which coke releases water; π stands for pi.

4. The optimization method according to claim 3, characterized in that: The r1, r2, R3, R4 and r5 are all homogeneous reactions; The expressions of R1, R2, R3, R4 and R5 are: in, k0 represents the pre-exponential factor; Cor k0 represents the correction coefficient for k0; E represents the activation energy of chemical reaction kinetics; Cor E represents the correction factor for E; R represents a constant; T g Indicates the reaction gas temperature; C A , C B Represents R 1-5 The concentrations of the corresponding reactants; r 1-5 Represents reaction 1, reaction 2, reaction 3, reaction 4, reaction 5; the expression of r6 is: in, represents the concentration of reactant methane; R7, R8, R9, R 10 and r 11 All are heterogeneous reactions; R7, R8, R9, R 10 The expression is: in k diff It indicates the rate at which gas diffuses to the coke surface; Cor kd Indicates k diff The correction factor of r 7-10 It indicates reaction 7, reaction 8, reaction 9, and reaction 10; Y represents the ratio of the unreacted carbon core diameter to the coke diameter; k dash k represents the rate at which gas diffuses through the outer ash layer to the coke surface; s represents the reaction rate on the surface of the carbon core; T p Indicates the surface temperature of coke particles; d p represents the diameter of coke particles; Indicates the effective partial pressure corresponding to reactions 7-10 The r 11 The expression is: in, r 11 Indicates reaction 11; F w represents the catalyst activity correction factor; P f represents the pressure correction factor; T f represents the temperature correction factor; represents the effective partial pressure corresponding to reaction 11.

5. The optimization method according to claim 4, characterized in that: The expression of the effective partial pressure corresponding to the reaction 7 is: in, represents the effective partial pressure corresponding to reaction 7; Indicates the pressure of water in the coal gasifier to be treated; Indicates the pressure of hydrogen in the coal gasifier to be treated; P CO Indicates the pressure of carbon monoxide in the coal gasifier to be treated; Cor a represents the correction coefficient for a; Cor b represents the correction coefficient for b; a represents the equilibrium coefficient of the first reaction; b represents the equilibrium coefficient of the second reaction; c represents the equilibrium coefficient of the third reaction; T p Indicates the surface temperature of coke particles; The expression of the effective partial pressure corresponding to the reaction 8 is: in, represents the effective partial pressure corresponding to reaction 8; Indicates the pressure of methane in the coal gasifier to be treated; The expression of the effective partial pressure corresponding to the reaction 9 is: in, represents the effective partial pressure corresponding to reaction 9; Indicates the pressure of oxygen in the coal gasifier to be treated; The expression of the effective partial pressure corresponding to the reaction 10 is: in, represents the effective partial pressure corresponding to reaction 10; Indicates the pressure of carbon dioxide in the coal gasifier to be treated; The expression of the effective partial pressure corresponding to the reaction 11 is: in, represents the effective partial pressure corresponding to reaction 11.

6. The optimization method according to claim 1, characterized in that: The method of obtaining the final concentration of the reaction components in the coal gasification reaction of the reactor to be treated after a preset time t based on the initial concentration through the kinetic parameter equation of the coal gasification reaction specifically includes: Based on the initial concentration and the kinetic parameter equation of the coal gasification reaction, the preset time t is iterated according to dt to calculate the concentration of each reaction component in the coal gasification reaction of the reactor to be treated, and the calculated value of each reaction component in the coal gasification reaction of the reactor to be treated when the final dt≥the preset time t is used as the final concentration of the reaction component in the coal gasification reaction of the reactor to be treated after the preset time t.

7. The optimization method according to claim 1, characterized in that: The method of iterating the preset time t according to dt based on the initial concentration and the kinetic parameter equation of the coal gasification reaction to calculate the concentration of each reaction component in the coal gasification reaction of the reactor to be treated specifically includes: Based on the initial concentration and the kinetic parameter equation of the coal gasification reaction, a calculation formula for obtaining the concentration of each reaction component in the coal gasification reaction of the reactor to be treated; Based on the calculation formula of the concentration of each reaction component in the coal gasification reaction of the reactor to be treated, the concentration of each reaction component in the coal gasification reaction of the reactor to be treated is iterated according to dt for the preset time t; The calculation formula for the concentration of each reaction component in the coal gasification reaction of the reactor to be treated is: in, r i An expression representing the rate at which component i is generated or consumed by a reaction; is the initial concentration of component i; c i is the concentration of component i after the preset time t.

8. The optimization method according to claim 1, characterized in that: The objective function is constructed based on the final concentration of the component to be corrected in the final concentration and the actual concentration of the component to be corrected, and the optimal reaction rate correction coefficient is obtained by optimizing the reaction rate correction coefficient of the coal gasification reaction, specifically including: Based on the final concentration, determining the final concentration of the component to be calibrated; Based on the actual concentrations of oxygen, methane, carbon monoxide, carbon dioxide, hydrogen sulfide, hydrogen, nitrogen, and water in the reactor to be treated, respectively determining the actual concentrations of the components to be corrected; constructing an objective function based on the actual concentration of the component to be corrected and the final concentration of the component to be corrected; The reaction rate correction coefficient of coal gasification reaction is optimized by BFGSBMinimizer algorithm to obtain the optimal reaction rate correction coefficient.

9. The optimization method according to claim 1, characterized in that: The expression of the objective function is: in, c i is the actual concentration of component i to be corrected. c c,i is the final concentration of component i to be calibrated; N represents the number of groups to be corrected; f(x) represents the objective function.

10. A system for optimizing coal gasification reaction rate correction coefficient, characterized in that: The optimization system comprises: An initial concentration acquisition module is used to acquire the initial concentration of the reaction components in the coal gasification reaction of the reactor to be processed; A final concentration calculation module, based on the initial concentration, obtains the final concentration of the reaction component in the coal gasification reaction of the reactor to be treated after a preset time t by using a kinetic parameter equation of the coal gasification reaction; The optimization module constructs an objective function based on the final concentration of the component to be corrected in the final concentration and the actual concentration of the component to be corrected, and obtains the optimal reaction rate correction coefficient by optimizing the reaction rate correction coefficient of the coal gasification reaction.

Images