Method and device for adjusting kinetic parameters, storage medium and electronic device

Abstract

The application discloses a method and device for adjusting a kinetic parameter, a storage medium and an electronic device. The method comprises: determining the polymerization heat of a reaction in a polymerization reactor according to a heat balance formula, wherein the heat balance formula comprises the mass flow of inlet material of the polymerization reactor, the mass enthalpy of the inlet material of the polymerization reactor, the heat taken away by the cooling medium in the polymerization reactor, the mass flow of outlet material of the polymerization reactor, the mass enthalpy of the outlet material in the polymerization reactor and the polymerization heat; determining the activity of a target catalyst according to the polymerization heat and the feeding amount of the target catalyst in the polymerization reactor, and determining the activity of a mixed catalyst in the polymerization reactor according to the activity of the target catalyst; and adjusting the pre-exponential factor in the kinetic parameter according to the activity of the mixed catalyst. The application solves the technical problem that the kinetic parameter in the model cannot be adjusted when the catalyst changes in the related art using the steady-state mechanism model.

Description

Technical Field

[0001] This application relates to the field of industrial automation, and more specifically, to a method, apparatus, storage medium, and electronic device for adjusting dynamic parameters. Background Technology

[0002] Polyolefins are important polymer materials with wide applications in the automotive industry, home appliances, electronics, packaging, building materials, and furniture. The most crucial aspect of polyolefin production is the control of the polymerization reaction stage, and key control parameters for this process include catalyst feed, which is directly related to catalyst activity.

[0003] In related technologies, catalyst activity is characterized by calculating polymer yield per unit catalyst feed using material balance. However, this method requires chromatographic data to calculate the polymer output from the reactor. Chromatographic data has a certain time lag, and online chromatography may produce inaccurate data due to factors such as long sampling pipelines, uneven heating, and component condensation. Therefore, using this method to characterize catalyst activity online may not be accurate.

[0004] Furthermore, the modeling of polymerization reactions in related technologies mainly focuses on steady-state mechanisms, neglecting changes in catalyst activity, especially when a new batch of catalyst is added to the reactor. This method is applicable when the catalyst remains unchanged, but during continuous online switching of catalyst grades, the addition of new catalyst and the removal of old catalyst occur simultaneously in the polymerization reactor. As the reaction progresses, the amount of new catalyst gradually becomes dominant, and the activity of the catalyst in the reactor gradually approaches that of a reactor entirely composed of new catalyst. If there is a significant difference in activity between the new and old catalysts, it will lead to large fluctuations in the component concentration within the reactor. When the activity difference between the new and old catalysts is significant, even if a dynamic model is constructed based on the original steady-state kinetic parameters, it cannot describe the changes in reactor component concentration.

[0005] There is currently no effective solution to the above problems. Summary of the Invention

[0006] This application provides a method, apparatus, storage medium, and electronic device for adjusting kinetic parameters, to at least solve the technical problem in related technologies where steady-state mechanism models are not applicable when the catalyst changes, and the kinetic parameters in the model cannot be adjusted.

[0007] According to one aspect of the embodiments of this application, a method for adjusting kinetic parameters is provided, comprising: determining the heat of polymerization in a polymerization reactor based on a heat balance formula, wherein the heat balance formula includes the mass flow rate of the inlet material of the polymerization reactor, the mass enthalpy of the inlet material of the polymerization reactor, the heat carried away by the cooling medium in the polymerization reactor, the mass flow rate of the outlet material of the polymerization reactor, the mass enthalpy of the outlet material of the polymerization reactor, and the heat of polymerization; determining the activity of a target catalyst based on the heat of polymerization and the feed rate of a target catalyst in the polymerization reactor, wherein the target catalyst includes a new catalyst and a used catalyst; determining the activity of a mixed catalyst in the polymerization reactor based on the activity of the target catalyst; and adjusting the pre-exponential factor in the kinetic parameters based on the activity of the mixed catalyst.

[0008] Optionally, the heat carried away by the cooling medium in the polymerization reactor is determined by: obtaining the mass of the cooling medium, the output temperature of the cooling medium in the polymerization reactor, the input temperature of the cooling medium in the polymerization reactor, and the average heat capacity of the cooling medium between the output temperature and the input temperature; and determining the heat carried away by the cooling medium in the polymerization reactor based on the mass of the cooling medium, the output temperature, the input temperature, and the average heat capacity.

[0009] Optionally, the activity of the target catalyst is determined based on the heat of polymerization and the feed rate of the target catalyst in the polymerization reactor, including: obtaining the feed rate of the monomer in the polymerization reactor; if the target catalyst is a new catalyst, obtaining the first feed rate of the new catalyst in the polymerization reactor; if the target catalyst is an old catalyst, obtaining the second feed rate of the old catalyst in the polymerization reactor; determining the first activity corresponding to the new catalyst based on the feed rate of the monomer, the heat of polymerization, and the first feed rate; and determining the second activity corresponding to the old catalyst based on the feed rate of the monomer, the heat of polymerization, and the second feed rate.

[0010] Optionally, the activity of the mixed catalyst in the polymerization reactor is determined based on the activity of the target catalyst, including: obtaining the average residence time of the polymerization reactor; determining the target proportion of the target catalyst in the polymerization reactor within a preset time period based on the average residence time, wherein the preset time period is the time from the start of the addition of new catalyst to the polymerization reactor to the preset time; and determining the activity of the mixed catalyst in the polymerization reactor based on the target proportion and the activity of the target catalyst.

[0011] Optionally, after determining the activity of the mixed catalyst in the polymerization reactor, the method further includes: determining the maximum activity of the old catalyst; and determining the multiplier factor of the kinetic parameters based on the activity of the mixed catalyst and the maximum activity of the old catalyst.

[0012] Optionally, based on the activity of the mixed catalyst, the pre-exponential factor in the kinetic parameters is adjusted, including: obtaining the multiplier factor and the pre-exponential factor of chain growth at the target active site; and determining the adjusted pre-exponential factor of chain growth at the target active site based on the multiplier factor and the pre-exponential factor of chain growth at the target active site.

[0013] Optionally, the chain growth pre-exponential factor at the target active site is determined by: obtaining the chain growth pre-exponential factor of a single active site, the mass ratio of the polymer generated at the target active site to the total polymer generated at all active sites, and the number of active sites of the new catalyst; and determining the chain growth pre-exponential factor at the target active site based on the chain growth pre-exponential factor of the single active site, the mass ratio, and the number of active sites of the new catalyst.

[0014] According to another aspect of the embodiments of this application, a kinetic parameter adjustment device is also provided, comprising: a first determining module, configured to determine the heat of polymerization in the polymerization reactor based on a heat balance formula, wherein the heat balance formula includes the mass flow rate of the inlet material of the polymerization reactor, the mass enthalpy of the inlet material of the polymerization reactor, the heat carried away by the cooling medium in the polymerization reactor, the mass flow rate of the outlet material of the polymerization reactor, the mass enthalpy of the outlet material of the polymerization reactor, and the heat of polymerization; a second determining module, configured to determine the activity of the target catalyst based on the heat of polymerization and the feed rate of the target catalyst in the polymerization reactor, wherein the target catalyst includes a new catalyst and a used catalyst; a third determining module, configured to determine the activity of the mixed catalyst in the polymerization reactor based on the activity of the target catalyst; and an adjusting module, configured to adjust the pre-exponential factor in the kinetic parameters based on the activity of the mixed catalyst.

[0015] According to another aspect of the embodiments of this application, an electronic device is also provided, including a memory for storing program instructions; and a processor connected to the memory for executing program instructions to perform the following functions: determining the heat of polymerization in the polymerization reactor based on a heat balance formula, wherein the heat balance formula includes the mass flow rate of the inlet material of the polymerization reactor, the mass enthalpy of the inlet material of the polymerization reactor, the heat carried away by the cooling medium in the polymerization reactor, the mass flow rate of the outlet material of the polymerization reactor, the mass enthalpy of the outlet material of the polymerization reactor, and the heat of polymerization; determining the activity of a target catalyst based on the heat of polymerization and the feed rate of a target catalyst in the polymerization reactor, wherein the target catalyst includes a new catalyst and a used catalyst; determining the activity of a mixed catalyst in the polymerization reactor based on the activity of the target catalyst; and adjusting the pre-exponential factor in the kinetic parameters based on the activity of the mixed catalyst.

[0016] According to another aspect of the embodiments of this application, a non-volatile storage medium is also provided, the non-volatile storage medium including a stored computer program, wherein the device where the non-volatile storage medium is located executes the above-mentioned method for adjusting dynamic parameters by running the computer program.

[0017] In this embodiment, the heat of polymerization in the polymerization reactor is determined based on a heat balance formula, which includes the mass flow rate of the inlet material, the mass enthalpy of the inlet material, the heat carried away by the cooling medium in the polymerization reactor, the mass flow rate of the outlet material, the mass enthalpy of the outlet material, and the heat of polymerization. The activity of the target catalyst is determined based on the heat of polymerization and the feed rate of the target catalyst in the polymerization reactor, where the target catalyst includes both new and old catalysts. The activity of the mixed catalyst in the polymerization reactor is determined based on the activity of the target catalyst. Based on the activity of the mixed catalyst, the pre-exponential factor in the kinetic parameters is adjusted, achieving the goal of rapidly and accurately calculating the catalyst activity. This realizes the technical effect of timely and dynamic adjustment of kinetic parameters, thereby solving the technical problem in related technologies where steady-state mechanism models are not applicable when the catalyst changes, and the kinetic parameters in the model cannot be adjusted. Attached Figure Description

[0018] The accompanying drawings, which are included to provide a further understanding of this application and form part of this application, illustrate exemplary embodiments and are used to explain this application, but do not constitute an undue limitation of this application. In the drawings:

[0019] Figure 1 This is a hardware structure block diagram of a computer terminal for implementing a method for adjusting dynamic parameters according to an embodiment of this application;

[0020] Figure 2 This is a flowchart of a method for adjusting dynamic parameters according to an embodiment of this application;

[0021] Figure 3a This is a flowchart illustrating the construction of a dynamic mechanism model according to an embodiment of this application;

[0022] Figure 3b This is a flowchart of an HDPE slurry process according to an embodiment of this application;

[0023] Figure 3c This is a graph showing the relationship between the relative activity of a novel catalyst and time according to an embodiment of this application;

[0024] Figure 3d This is a schematic diagram showing the particle ratio of new and old catalysts in a reactor according to an embodiment of this application;

[0025] Figure 3e This is a schematic diagram showing the relative activity of a mixed catalyst in a reactor according to an embodiment of this application;

[0026] Figure 3f This is a schematic diagram of the actual measured value and simulated calculated value of the change in ethylene concentration during a grade switching process, according to an embodiment of this application.

[0027] Figure 4 This is a structural diagram of a dynamic parameter adjustment device according to an embodiment of this application. Detailed Implementation

[0028] 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 with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present application, and not all embodiments. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative effort should fall within the scope of protection of the present application.

[0029] It should be noted that the terms "first," "second," etc., in the specification, claims, and accompanying drawings of this application are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that the embodiments of this application described herein can be implemented in orders other than those illustrated or described herein. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover non-exclusive inclusion; for example, a process, method, system, product, or apparatus that comprises a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or apparatus.

[0030] First, some nouns or terms that appear in the explanation of the embodiments of this application shall be interpreted as follows:

[0031] Steady-state mechanism model: The component concentration in the reactor does not change with time. At the same time, from a kinetic perspective, a reaction rate equation model covering the main elementary reactions of polymerization (including catalyst activation, chain initiation, chain growth, chain termination, etc.) is constructed to establish the relationship between component concentration and kinetic parameters, temperature, pressure, etc.

[0032] Dynamic mechanism model: The component concentration in the reactor changes over time, and the constructed reaction rate equation model includes time variables, while the others are consistent with the steady-state mechanism model.

[0033] Phillips catalysts: The main active component of the catalyst is chromium oxide, also known as chromium-based catalysts, and they have a certain induction period.

[0034] Induction period: The catalyst induction period refers to the initial stage of the reaction when the catalyst has low activity and the reaction rate approaches 0.

[0035] Monomer: A monomer is a small molecule that can be covalently linked with the same or other molecules to form a polymer. For example, the monomer of polyethylene is ethylene, and the main monomer of polypropylene is propylene.

[0036] A fully mixed-flow reactor: The reaction fluids are completely mixed inside the reactor, and have a uniform temperature and composition within the reactor, consistent with the temperature and composition of the material flowing out of the reactor.

[0037] Chain-growing elementary reactions: One of the main steps in polymerization, where active molecules and monomer molecules interact to generate long-chain active molecules. Almost all monomers are consumed in this process during polymerization, and this process determines the molecular structure of the resulting polymer.

[0038] The method for adjusting dynamic parameters provided in this application can be executed on a mobile terminal, computer terminal, or similar computing device. Figure 1 A hardware block diagram of a computer terminal for implementing a method for adjusting dynamic parameters is shown. Figure 1 As shown, the computer terminal 10 may include one or more processors (shown as 102a, 102b, ..., 102n in the figure) (the processor may include, but is not limited to, a microprocessor MCU or a programmable logic device FPGA, etc.), a memory 104 for storing data, and a transmission module 106 for communication functions. In addition, it may also include: a display, an input / output interface (I / O interface), a universal serial bus (USB) port (which may be included as one of the ports of the I / O interface), a network interface, a power supply, and / or a camera. Those skilled in the art will understand that... Figure 1 The structure shown is for illustrative purposes only and does not limit the structure of the aforementioned electronic device. For example, computer terminal 10 may also include... Figure 1 The more or fewer components shown, or having the same Figure 1 The different configurations shown.

[0039] It should be noted that the aforementioned one or more processors and / or other data processing circuits are generally referred to herein as "data processing circuits". These data processing circuits may be implemented wholly or partially as software, hardware, firmware, or any other combination thereof. Furthermore, the data processing circuits may be a single, independent processing module, or may be wholly or partially integrated into any other element in the computer terminal 10. As involved in the embodiments of this application, the data processing circuits serve as processor control (e.g., selection of a variable resistor termination path connected to an interface).

[0040] The memory 104 can be used to store software programs and modules of application software, such as the program instructions / data storage device corresponding to the dynamic parameter adjustment method in the embodiments of this application. The processor executes various functional applications and data processing by running the software programs and modules stored in the memory 104, thereby realizing the aforementioned dynamic parameter adjustment method. The memory 104 may include high-speed random access memory, and may also include non-volatile memory, such as one or more magnetic storage devices, flash memory, or other non-volatile solid-state memory. In some instances, the memory 104 may further include memory remotely located relative to the processor, and these remote memories can be connected to the computer terminal 10 via a network. Examples of such networks include, but are not limited to, the Internet, corporate intranets, local area networks, mobile communication networks, and combinations thereof.

[0041] The transmission module 106 is used to receive or send data via a network. Specific examples of the network described above may include a wireless network provided by the communication provider of the computer terminal 10. In one example, the transmission module 106 includes a Network Interface Controller (NIC), which can connect to other network devices via a base station to communicate with the Internet. In another example, the transmission module 106 may be a Radio Frequency (RF) module, used for wireless communication with the Internet.

[0042] The display may be, for example, a touchscreen liquid crystal display (LCD) that allows the user to interact with the user interface of the computer terminal 10.

[0043] It should be noted here that, in some optional embodiments, the above... Figure 1 The computer terminal shown may include hardware elements (including circuitry), software elements (including computer code stored on a computer-readable medium), or a combination of both hardware and software elements. It should be noted that... Figure 1 This is only one instance of a specific particular instance, and is intended to illustrate the types of components that may exist in the aforementioned computer terminal.

[0044] In the above operating environment, this application provides an embodiment of a method for adjusting dynamic parameters. It should be noted that the steps shown in the flowchart in the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions. Although a logical order is shown in the flowchart, in some cases, the steps shown or described may be executed in a different order than that shown here.

[0045] Figure 2 This is a flowchart of a method for adjusting dynamic parameters according to an embodiment of this application, such as... Figure 2 As shown, the method includes the following steps:

[0046] Step S202: Determine the heat of polymerization in the polymerization reactor based on the heat balance formula. The heat balance formula includes the mass flow rate of the inlet material of the polymerization reactor, the mass enthalpy of the inlet material of the polymerization reactor, the heat carried away by the cooling medium in the polymerization reactor, the mass flow rate of the outlet material of the polymerization reactor, the mass enthalpy of the outlet material of the polymerization reactor, and the heat of polymerization.

[0047] In the method for adjusting kinetic parameters provided in the embodiments of this application, the heat carried away by the cooling medium in the polymerization reactor is determined by: obtaining the mass of the cooling medium, the output temperature of the cooling medium in the polymerization reactor, the input temperature of the cooling medium in the polymerization reactor, and the average heat capacity of the cooling medium between the output temperature and the input temperature; and determining the heat carried away by the cooling medium in the polymerization reactor based on the mass of the cooling medium, the output temperature, the input temperature, and the average heat capacity.

[0048] In step S202 above, the material at the reactor inlet may include, for example, monomer feed, comonomer feed, diluent feed, etc. The heat carried away by the cooling medium in the polymerization reactor can be determined by the temperature difference between the cooling medium entering and exiting the reactor jacket, the mass of the cooling medium, and the average heat capacity of the cooling medium within the temperature range of the reactor jacket inlet and outlet.

[0049] In some embodiments of this application, the heat balance formula is shown in the following formula (1):

[0050] ∑F ini H ini +ΔH=Q+∑F oi H oi (1)

[0051] Where Q represents the heat carried away by the cooling medium in the polymerization reactor, and Q is determined by formula (2):

[0052]

[0053] In the above formula (1), F iniH represents the mass flow rate of component i in the inlet material of the polymerization reactor, in kg / h. ini ΔH represents the mass enthalpy of component i in the inlet material of the polymerization reactor, in kJ / (kg / h); ΔH represents the heat of polymerization, in kJ / h; Q represents the heat removed by the cooling medium in the polymerization reactor, in kJ / h; F oi H represents the mass flow rate of component i in the effluent from the polymerization reactor, expressed in kg / h. oi The mass enthalpy of component i in the effluent from the polymerization reactor is expressed in kJ / (kg / h).

[0054] In the above formula (2), t1 is the average heat capacity of the cooling medium within the inlet and outlet temperature range of the polymerization reactor jacket, expressed in kJ / (kg / h) / ℃; m is the mass of the cooling medium, expressed in kg / h; t2 is the temperature of the cooling medium exiting the polymerization reactor jacket, i.e., the output temperature of the cooling medium in the polymerization reactor, expressed in ℃; and t1 is the temperature of the cooling medium entering the polymerization reactor jacket, i.e., the input temperature of the cooling medium in the polymerization reactor, expressed in ℃.

[0055] Step S204: Determine the activity of the target catalyst based on the heat of polymerization and the feed rate of the target catalyst in the polymerization reactor, wherein the target catalyst includes both new and old catalysts.

[0056] In step S204 above, determining the activity of the new catalyst requires obtaining the dynamic characteristics of the catalyst just added to the polymerization reactor. Therefore, the activity of the new catalyst is determined using data from the start-up period, which is the time from when the new catalyst is added to the polymerization reactor until the activity of the catalyst is completely stable. The activity of the old catalyst is basically very stable in the reactor, and the activity is fully released. Therefore, determining the activity of the old catalyst can basically be done by selecting data from stable production conditions.

[0057] Step S206: Determine the activity of the mixed catalyst in the polymerization reactor based on the activity of the target catalyst.

[0058] In step S206 above, for a continuously producing reactor, new catalyst is constantly being added while old catalyst is constantly being discharged. Their proportions in the fully mixed flow reactor conform to the residence time distribution law of the fully mixed flow reactor. Therefore, the activity of the mixed catalyst in the polymerization reactor can be determined by using the proportion data based on the activity of the target catalyst.

[0059] Step S208: Adjust the pre-exponential factor in the kinetic parameters according to the activity of the mixed catalyst.

[0060] In step S208 above, the adjustment coefficient of the pre-exponential factor is determined by the mixed catalyst activity, and then the pre-exponential factor in the kinetic parameters is adjusted to construct a new dynamic mechanism model. Figure 3a This is a flowchart illustrating the construction of a dynamic mechanism model according to an embodiment of this application, such as... Figure 3a As shown, the physical properties (including density, heat capacity, enthalpy, etc.) of each component in the polymerization reactor are first calculated using the PC-SAFT equation of state. Then, based on the molecular weight distribution data obtained from gel permeation chromatography (GPC) analysis of the product sample, the number of catalyst active sites, the molecular weight of polymer produced at each catalyst active site, and the mass ratio of polymer produced at each catalyst active site to the total polymer are obtained by deconvolution of the molecular weight distribution. Combined with the conversion data under typical reaction conditions, the kinetic parameters (including pre-exponential factor and activation energy) of each elementary reaction at multiple active sites are comprehensively adjusted and determined, thereby constructing the original dynamic model. For dynamic processes such as grade switching or load changes, the adjustment coefficient of the pre-exponential factor in the kinetic parameters is determined based on the activity of the mixed catalyst. The kinetic parameters are then adjusted, and a set of time differential equations for the concentration of reaction components is constructed based on material balance, thus completing the construction of a new dynamic mechanism model under conditions of large changes in catalyst activity.

[0061] Based on steps S202 to S208 above, the goal of rapidly and accurately calculating the catalyst activity is achieved, thereby realizing the technical effect of timely and dynamic adjustment of kinetic parameters. This solves the technical problem in related technologies where steady-state mechanism models are inapplicable when catalysts change, and the kinetic parameters in the model cannot be adjusted. A detailed explanation follows.

[0062] Optionally, in the kinetic parameter adjustment method provided in the embodiments of this application, the activity of the target catalyst is determined based on the heat of polymerization and the feed amount of the target catalyst in the polymerization reactor. Specifically, this includes the following steps: obtaining the feed amount of the monomer in the polymerization reactor; if the target catalyst is a new catalyst, obtaining the first feed amount of the new catalyst in the polymerization reactor; if the target catalyst is an old catalyst, obtaining the second feed amount of the old catalyst in the polymerization reactor; determining the first activity corresponding to the new catalyst based on the feed amount of the monomer, the heat of polymerization, and the first feed amount; and determining the second activity corresponding to the old catalyst based on the feed amount of the monomer, the heat of polymerization, and the second feed amount.

[0063] In some embodiments of this application, the activities of both the new and old catalysts are determined according to formula (3):

[0064]

[0065] Where Act represents the activity of the catalyst. When calculating the activity of a new catalyst, Act represents the first activity mentioned above, and when calculating the activity of an old catalyst, Act represents the second activity mentioned above; ΔH represents the heat of polymerization of the polymerization reaction, in kJ / h, calculated from the above formulas (1) to (2), F in This refers to the feed rate of the monomer, expressed in kg / h; F cat F represents the catalyst feed rate, expressed in kg / h. This value is used when calculating the activity of a new catalyst. cat F represents the feed rate of the new catalyst, i.e., the first feed rate mentioned above. When calculating the activity of the old catalyst, F... cat This indicates the feed rate of the old catalyst, i.e., the second feed rate mentioned above; monomer refers to the monomer in the polymerization reaction (e.g., the monomer of polyethylene is ethylene). The materials at the inlet of the polymerization reactor include monomer feed, comonomer feed, diluent feed, etc. For example, the Innovene S process refers to the activity of catalysts such as ethylene feed, hexene feed, and isobutane feed, respectively. Here, it represents the heat of polymerization under a unit catalyst feed and a unit monomer feed (ethylene feed).

[0066] It should be noted that when calculating the activity of the new and old catalysts, the data for the new catalyst is selected based on the start-up time of the unit, while the data for the old catalyst is selected based on the stable production conditions. The calculation formula is formula (3).

[0067] The reason for choosing data from the plant's start-up period is that the exothermic data at this time is not affected by the old catalyst. For example, catalyst A is considered the old catalyst (producing grade C), and catalyst B is considered the new catalyst (producing grade D). During the online switch from grade C to grade D production, the new catalyst B is added to the reactor, and the old catalyst A is no longer added. Simultaneously, the catalyst at the reactor outlet flows out with the polymer. However, the catalyst has a certain residence time and distribution during the reaction. Therefore, the proportion of the new catalyst in the reactor will gradually increase, and the proportion of the old catalyst will gradually decrease. The activity of the mixed catalyst is related to the ratio of the new and old catalysts and the activity of each individual catalyst. When calculating the activity of each individual catalyst, the situation of other catalysts is not considered. Therefore, when calculating the activity of the new catalyst, the actual selected time period is the plant start-up data during the production of grade D. The start-up time operating data is chosen because there is no catalyst in the reactor before start-up. After start-up, the catalyst is slowly added to the reactor, and there is an induction period, which can fully reflect the dynamic characteristics of the polymerization reaction when the catalyst is inactive.

[0068] Optionally, in the method for adjusting kinetic parameters provided in the embodiments of this application, the activity of the mixed catalyst in the polymerization reactor is determined based on the activity of the target catalyst, specifically including the following steps: obtaining the average residence time of the polymerization reactor; determining the target proportion of the target catalyst in the catalyst in the polymerization reactor within a preset time period based on the average residence time, wherein the preset time period is the time corresponding to the start of the addition of new catalyst to the polymerization reactor to a preset time; and determining the activity of the mixed catalyst in the polymerization reactor based on the target proportion and the activity of the target catalyst.

[0069] In some embodiments of this application, the ratio of new to old catalyst can be obtained based on the residence time distribution of materials in a fully mixed reactor. Figure 3b This is a flowchart of an HDPE slurry process according to an embodiment of this application. In the specific production process of high-density polyethylene (HDPE) using the slurry method, the monomer component is ethylene, the comonomer is hexene, the diluent is isobutane, and the catalyst is a Phillips catalyst (chromium-based catalyst). The molecular weight and melt index of HDPE are controlled by temperature adjustment, and the catalyst has a certain induction period. The polymerization reactor is a loop reactor, which can be regarded as a fully mixed flow reactor. The heat released by polyethylene is carried away by the jacket water in the loop reactor. In this process flow of this case, data from the start-up stage of the unit (from ethylene feed of 0 to normal production load) are selected to calculate the heat released by the reactor Q and the catalyst activity Act. i The relative activity of the new catalyst can be calculated using formula (4), which normalizes the activity of the new catalyst and characterizes the activity of the new catalyst at the maximum relative activity of the new catalyst. Figure 3c The graph showing the relationship between the relative activity of the new catalyst and time is presented.

[0070]

[0071] In formula (4), Act1 represents the activity of the new catalyst, Act... max This represents the maximum activity of the new catalyst. The above formula (4) is a normalization of the new catalyst activity; it should originally be:

[0072]

[0073] However, for dynamic processes, Act min The lowest activity of the new catalyst is the activity when the new catalyst is first added and is in the induction period, at which point the activity can be approximated as 0. Therefore, the normalization of the activity of the new catalyst can be achieved by using the above formula (4).

[0074] Based on the residence time distribution pattern in a fully mixed reactor (formulas (5) to (6)), the proportion of new and old catalyst particles within the time interval 0 to τ can be calculated. Figure 3d The diagram illustrates the particle ratio of new and old catalysts in the reactor. Based on the particle ratio of each catalyst in the fully mixed-flow reactor and the activities of the new and old catalysts, the weighted activity of the mixed catalyst in the reactor is calculated. Figure 3e A schematic diagram showing the relative activities of the mixed catalysts in the reactor is shown.

[0075] Specifically, the proportion of new catalyst in the polymerization reactor during the time period 0 to t is shown in formula (5), and the proportion of old catalyst in the polymerization reactor during the time period 0 to t is shown in formula (6):

[0076] θ1=1-e -tτ (5)

[0077] θ2=e -tτ (6)

[0078] In formulas (5) to (6), t is the time elapsed from the addition of the new catalyst to time t, which is the preset time mentioned above. τ is the average residence time of the polymerization reactor. θ1 represents the proportion corresponding to the new catalyst, and θ2 represents the proportion corresponding to the old catalyst. θ1 and θ2 are collectively referred to as the target proportion mentioned above.

[0079] After determining the activity and proportion of the new and old catalysts according to formulas (1) to (6), the activity of the mixed catalyst in the reactor can be determined by formula (7):

[0080] Act t =θ1γ1+θ2γ2 (7)

[0081] In formula (7), γ1 is the relative activity of the new catalyst in the dynamic process reactor, and γ2 is the relative activity of the old catalyst in the dynamic process reactor.

[0082] Optionally, in the method for adjusting kinetic parameters provided in the embodiments of this application, after determining the activity of the mixed catalyst in the polymerization reactor, the method further includes the following steps: determining the maximum activity of the old catalyst; and determining the multiplier factor of the kinetic parameters based on the activity of the mixed catalyst and the maximum activity of the old catalyst.

[0083] In some embodiments of this application, the maximum activity of the old catalyst is determined based on data from three times the residence time under past stable production conditions, and the multiplier factor can be determined according to formula (8):

[0084]

[0085] In formula (8), Actt For the activity of mixed catalysts, Act 2max λ represents the maximum activity of the old catalyst. λ is the multiplier factor that adjusts the pre-exponential factor of the kinetic parameters due to changes in catalyst activity during the dynamic process.

[0086] Optionally, in the method for adjusting kinetic parameters provided in the embodiments of this application, the pre-exponential factor in the kinetic parameters is adjusted according to the activity of the mixed catalyst, specifically including the following steps: obtaining the multiplier factor and the pre-exponential factor of chain growth on the target active site; and determining the adjusted pre-exponential factor of chain growth on the target active site based on the multiplier factor and the pre-exponential factor of chain growth on the target active site.

[0087] In some embodiments of this application, high-temperature gel permeation chromatography (GPC) was used to analyze the molecular weight distribution (MWD) of the reactor powder to determine that the number of active sites of the catalyst was 5, and the kinetic parameters at each active site were shown in Table 1. Based on the calculated catalyst activity, the pre-exponential factor of ethylene on the ethylene chain was adjusted, and the adjusted pre-exponential factor is shown in Table 2.

[0088] Table 1. Kinetic parameters of multiple active sites (pre-exponential factors)

[0089] Catalyst activation l / mol / s 200 200 200 200 200 Ethylene chain trigger l / mol / s 374 374 374 374 374 Hexene chain initiation l / mol / s 22 22 22 22 22 Ethylene in the ethylene chain growth l / mol / s 694.4 235.2 1064 1070.3 235.2 Ethylene chain growth l / mol / s 13.7 17.2 2.5 5 8.1 Hexene in ethylene chain growth l / mol / s 118.4 41.8 135 161.2 38.6 ethylene chain transfer to ethylene l / mol / s 1.2 1.6 0.4 0.3 3e-2 ethylene chain transfer to hexene l / mol / s 0.1 0.1 0.1 0.1 0.1 Hexene chain transfer to ethylene l / mol / s 0.01 0.01 0.01 0.01 0.01 Ethylene chain self-transfer l / mol / s 0.1 0.1 0.01 0.1 0.01 Hexene chain self-transfer l / mol / s 1.2 1.2 0.01 1.2 0.01 Chain self-deactivation 1 / s 1.2e-6 1.2e-6 1.2e-6 1.2e-6 1.2e-6

[0090] Table 2. Adjusted pre-exponential factors of ethylene on the ethylene chain for multi-active-site catalysts.

[0091]

[0092]

[0093] Based on Table 3, the generation and consumption rates of each species in the polymerization reactor were calculated, and material balance was performed accordingly.

[0094] Table 3. Generation and consumption rates of various species in the reactor.

[0095]

[0096] In Table 3, the physical meaning of each letter is as follows: These represent the 0th, first, and second moments of the active chain, respectively, where:

[0097]

[0098]

[0099]

[0100] n represents the chain length, [P] n ] represents an active chain of length n.

[0101] , respectively, represent the rate constants for catalyst activation reaction, monomer i-initiated reaction, self-transfer reaction, chain transfer to monomer i, self-deactivation reaction, and chain growth reaction.

[0102] [Mi] represents the concentration of monomer component i, [S] p [] indicates the concentration of the catalyst; [P0] indicates the concentration of vacant active sites.

[0103] For catalysts, based on material balance:

[0104] For ethylene, according to the material balance:

[0105] For hexene, according to the material balance:

[0106] For isobutane, according to the material balance:

[0107] For the 0th moment of the active chain, we have:

[0108] With respect to the first moment of the active chain:

[0109] With respect to the second moment of the active chain:

[0110] For the zeroth moment of the dead chain, we have:

[0111] For the first moment of the dead chain, we have:

[0112] For the second moment of the dead chain, we have:

[0113] Regarding the feed:

[0114] For the discharge, based on the material balance, we have:

[0115]

[0116] (F cat ) out =(F cat ) in

[0117] In the above formula,

[0118] V represents the volume of the polymerization reactor, in cubic meters (m³). 3 ;

[0119] m cat This is the mass flow rate of the catalyst, expressed in kg / h.

[0120] m M1 This is the mass flow rate of monomer (ethylene), expressed in kg / h.

[0121] m M2 It is the mass flow rate of the comonomer (hexene), expressed in kg / h;

[0122] MW cat It is the molar mass of the catalyst, expressed in kg / kmol;

[0123] MW M1 It is the molar mass of the monomer (ethylene), expressed in kg / kmol;

[0124] MW M2 It is the molar mass of the comonomer (hexene), expressed in kg / kmol;

[0125] (F M1 ) in 、(F M1 ) out These are the molar flow rates of monomer (ethylene) at the inlet and outlet of the polymerization reactor, respectively.

[0126] (F M2 ) in 、(F M2 ) out These are the molar flow rates of the comonomer (hexene) at the reactor inlet and outlet, respectively.

[0127] (F x ) in 、(F x ) out This represents the molar flow rate of isobutane at the reactor inlet and outlet.

[0128] R M1 It is the rate at which monomer (ethylene) is consumed in the polymerization reaction, expressed in mol / m³ / h.

[0129] R M2 It is the rate at which the comonomer (hexene) is consumed in the polymerization reaction, expressed in mol / m³ / h.

[0130] In summary, a dynamic mechanism model for HDPE in a fully mixed-flow reactor can be constructed. Furthermore, by adjusting the kinetic parameters of ethylene chain growth on the ethylene chain based on changes in catalyst activity, a dynamic mechanism model can be constructed. This model can then calculate the trend of ethylene concentration changes during grade switching, thereby optimizing the catalyst feed. Figure 3fThis is a schematic diagram showing the actual measured values ​​and simulated calculated values ​​of ethylene concentration changes during the grade switching process. Figure 3f It can be seen that the model simulation effect is good.

[0131] After determining the multiplier factor, the pre-exponential chain growth factor on the target active site can be determined using formula (9):

[0132] k′ p (j)=λ*k p (j) (9)

[0133] Wherein, λ represents the activity of the mixed catalyst, which can be calculated using formula (8). p (j) is the chain growth pre-exponential factor at the active site j, that is, the chain growth pre-exponential factor at the target active site mentioned above. As shown in Table 1-2, there are 5 active sites, that is, the value of j ranges from 1 to 5. k′ p (j) is the pre-exponential chain growth factor on the adjusted active site i.

[0134] Optionally, in the method for adjusting kinetic parameters provided in the embodiments of this application, the chain growth pre-exponential factor on the target active site is determined in the following manner: obtaining the chain growth pre-exponential factor of a single active site, the mass ratio of the polymer generated at the target active site to the total polymer generated at all active sites, and the number of active sites of the new catalyst; and determining the chain growth pre-exponential factor on the target active site based on the chain growth pre-exponential factor of the single active site, the mass ratio, and the number of active sites of the new catalyst.

[0135] In the embodiments of this application, k p (j) can be determined by the following formula (10):

[0136] k p (j)=k p *mass i *N (10)

[0137] In formula (10), k p It is a pre-exponential factor for chain growth at a single active site, mass i It is the mass ratio of polymer generated at active site i to the total polymer generated at all active sites, where N is the number of active sites in the new catalyst, and mass is... i N was obtained by Flory peak fitting of the molecular weight distribution (MWD) measured by GPC (gel permeation chromatography). k was adjusted. p To ensure that the polymer yield and molecular weight are consistent with the actual polymer yield and the molecular weight obtained from GPC testing, k was finally determined. p .

[0138] The dynamic parameter adjustment method according to the embodiments of this application has the following advantages:

[0139] (1) The instantaneous activity of the catalyst can be calculated quickly and accurately, especially after the addition of a new catalyst. The calculation of the overall activity of the catalyst is more accurate and faster than the previous method of using polymer yield per unit catalyst.

[0140] (2) Based on the instantaneous catalyst activity calculation method proposed in the embodiments of this application and the proportion of new and old catalysts in the fully mixed reactor, the overall activity of the catalyst in the polymerization reactor and its relative activity relative to when the polymerization reactor is entirely composed of old catalyst are calculated. Based on the relative activity coefficient of the catalyst, the pre-exponential factor in the chain growth elementary reaction rate constant is adjusted to complete the construction of the dynamic mechanism model, which can accurately describe the change in the concentration of components in the reaction after the addition of new catalyst.

[0141] Figure 4 This is a structural diagram of a dynamic parameter adjustment device according to an embodiment of this application. Figure 4 As shown, the device includes:

[0142] The first determining module 402 is used to determine the heat of polymerization of the reaction in the polymerization reactor based on the heat balance formula, wherein the heat balance formula includes the mass flow rate of the inlet material of the polymerization reactor, the mass enthalpy of the inlet material of the polymerization reactor, the heat carried away by the cooling medium in the polymerization reactor, the mass flow rate of the outlet material of the polymerization reactor, the mass enthalpy of the outlet material of the polymerization reactor, and the heat of polymerization.

[0143] The second determining module 404 is used to determine the activity of the target catalyst based on the heat of polymerization and the feed rate of the target catalyst in the polymerization reactor, wherein the target catalyst includes new catalyst and old catalyst.

[0144] The third determining module 406 is used to determine the activity of the mixed catalyst in the polymerization reactor based on the activity of the target catalyst.

[0145] Adjustment module 408 is used to adjust the pre-exponential factor in the kinetic parameters based on the activity of the mixed catalyst.

[0146] Through the first determining module 402, the second determining module 404, the third determining module 406, and the adjusting module 408 in the aforementioned kinetic parameter adjustment device, the activity of the catalyst can be calculated quickly and accurately, thereby achieving the technical effect of timely and dynamic adjustment of kinetic parameters. This solves the technical problem in related technologies where steady-state mechanism models are not applicable when the catalyst changes, and the kinetic parameters in the model cannot be adjusted.

[0147] Optionally, in the kinetic parameter adjustment device provided in the embodiments of this application, the first acquisition module is further used to acquire the mass of the cooling medium, the output temperature of the cooling medium in the polymerization reactor, the input temperature of the cooling medium in the polymerization reactor, and the average heat capacity of the cooling medium between the output temperature and the input temperature; and to determine the heat carried away by the cooling medium in the polymerization reactor based on the mass of the cooling medium, the output temperature, the input temperature, and the average heat capacity.

[0148] Optionally, in the kinetic parameter adjustment device provided in the embodiments of this application, the second determining module is further configured to obtain the feed amount of monomer in the polymerization reactor; if the target catalyst is a new catalyst, obtain the first feed amount of the new catalyst in the polymerization reactor; if the target catalyst is an old catalyst, obtain the second feed amount of the old catalyst in the polymerization reactor; determine the first activity corresponding to the new catalyst based on the feed amount of monomer, the heat of polymerization and the first feed amount; and determine the second activity corresponding to the old catalyst based on the feed amount of monomer, the heat of polymerization and the second feed amount.

[0149] Optionally, in the kinetic parameter adjustment device provided in the embodiments of this application, the third determining module is further used to obtain the average residence time of the polymerization reactor; determine the target proportion of the target catalyst in the polymerization reactor within a preset time period based on the average residence time, wherein the preset time period is the time corresponding to the start of the addition of new catalyst to the polymerization reactor to a preset time; and determine the activity of the mixed catalyst in the polymerization reactor based on the target proportion and the activity of the target catalyst.

[0150] Optionally, in the kinetic parameter adjustment device provided in the embodiments of this application, the third determining module is further used to determine the maximum activity of the old catalyst; and to determine the multiplier factor of the kinetic parameters based on the activity of the mixed catalyst and the maximum activity of the old catalyst.

[0151] Optionally, in the kinetic parameter adjustment device provided in this application embodiment, the processing module is further configured to obtain the multiplier factor and the chain growth pre-exponential factor on the target active site; and determine the adjusted chain growth pre-exponential factor on the target active site based on the multiplier factor and the chain growth pre-exponential factor on the target active site.

[0152] Optionally, in the kinetic parameter adjustment device provided in the embodiments of this application, the adjustment module is further used to obtain the chain growth pre-exponential factor of a single active site, the mass ratio of polymer generated at the target active site to the total polymer generated at all active sites, and the number of active sites of the new catalyst; and to determine the chain growth pre-exponential factor at the target active site based on the chain growth pre-exponential factor of the single active site, the mass ratio, and the number of active sites of the new catalyst.

[0153] It should be noted that, Figure 4 The adjustment device for the dynamic parameters shown is used to perform... Figure 2 The method for adjusting the dynamic parameters shown above also applies to the device for adjusting these dynamic parameters, and will not be repeated here.

[0154] This application also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements the above-mentioned method for adjusting dynamic parameters.

[0155] This application also provides a non-volatile storage medium, which includes a stored computer program. The device containing the non-volatile storage medium executes the above-mentioned method for adjusting dynamic parameters by running the computer program, which will not be described in detail here.

[0156] The sequence numbers of the embodiments in this application are for descriptive purposes only and do not represent the superiority or inferiority of the embodiments.

[0157] In the above embodiments of this application, the descriptions of each embodiment have different focuses. For parts not described in detail in a certain embodiment, please refer to the relevant descriptions of other embodiments.

[0158] In the several embodiments provided in this application, it should be understood that the disclosed technical content can be implemented in other ways. The device embodiments described above are merely illustrative; for example, the division of units can be a logical functional division, and in actual implementation, there may be other division methods. For instance, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the displayed or discussed mutual coupling, direct coupling, or communication connection may be through some interfaces; the indirect coupling or communication connection between units or modules may be electrical or other forms.

[0159] The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple units. Some or all of the units can be selected to achieve the purpose of this embodiment according to actual needs.

[0160] Furthermore, the functional units in the various embodiments of this application can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. The integrated unit can be implemented in hardware or as a software functional unit.

[0161] If the integrated unit is implemented as a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, or all or part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as a USB flash drive, read-only memory (ROM), random access memory (RAM), portable hard drive, magnetic disk, or optical disk.

[0162] The above description is only a preferred embodiment of this application. It should be noted that for those skilled in the art, several improvements and modifications can be made without departing from the principle of this application, and these improvements and modifications should also be considered within the scope of protection of this application.

Claims

1. A method for adjusting dynamic parameters, characterized in that, include: The heat of polymerization in the polymerization reactor is determined based on the heat balance formula, wherein the heat balance formula includes the mass flow rate of the inlet material of the polymerization reactor, the mass enthalpy of the inlet material of the polymerization reactor, the heat carried away by the cooling medium in the polymerization reactor, the mass flow rate of the outlet material of the polymerization reactor, the mass enthalpy of the outlet material of the polymerization reactor, and the heat of polymerization. The activity of the target catalyst is determined based on the heat of polymerization and the feed rate of the target catalyst in the polymerization reactor, wherein the target catalyst includes both new and old catalysts; Based on the activity of the target catalyst, the activity of the mixed catalyst in the polymerization reactor is determined; Based on the activity of the mixed catalyst, the pre-exponential factor in the kinetic parameters is adjusted.

2. The method according to claim 1, characterized in that, The heat carried away by the cooling medium in the polymerization reactor is determined in the following way: The mass of the cooling medium, the output temperature of the cooling medium in the polymerization reactor, the input temperature of the cooling medium in the polymerization reactor, and the average heat capacity of the cooling medium between the output temperature and the input temperature are obtained. The amount of heat carried away by the cooling medium in the polymerization reactor is determined based on the mass of the cooling medium, the output temperature, the input temperature, and the average heat capacity.

3. The method according to claim 1, characterized in that, The activity of the target catalyst is determined based on the heat of polymerization and the feed rate of the target catalyst in the polymerization reactor, including: Obtain the feed rate of the monomer in the polymerization reactor; When the target catalyst is the new catalyst, the first feed rate of the new catalyst in the polymerization reactor is obtained; In the case that the target catalyst is the old catalyst, a second feed rate of the old catalyst in the polymerization reactor is obtained; The first activity of the new catalyst is determined based on the feed rate of the monomer, the heat of polymerization, and the first feed rate. The second activity corresponding to the old catalyst is determined based on the feed rate of the monomer, the heat of polymerization, and the second feed rate.

4. The method according to claim 1, characterized in that, Determining the activity of the mixed catalyst in the polymerization reactor based on the activity of the target catalyst includes: Obtain the average residence time of the polymerization reactor; Based on the average residence time, a target proportion of the target catalyst in the polymerization reactor is determined within a preset time period, wherein the preset time period is the time from the start of the addition of the new catalyst to the polymerization reactor to the preset time. The activity of the mixed catalyst in the polymerization reactor is determined based on the target ratio and the activity of the target catalyst.

5. The method according to claim 1, characterized in that, After determining the activity of the mixed catalyst in the polymerization reactor, the method further includes: Determine the maximum activity of the old catalyst; The multiplier factor of the kinetic parameter is determined based on the activity of the mixed catalyst and the maximum activity of the old catalyst.

6. The method according to claim 5, characterized in that, Based on the activity of the mixed catalyst, the pre-exponential factor in the kinetic parameters is adjusted, including: Obtain the multiplier factor and the chain growth pre-exponential factor on the target active site; Based on the multiplier factor and the chain growth pre-exponential factor on the target active site, the adjusted chain growth pre-exponential factor on the target active site is determined.

7. The method according to claim 6, characterized in that, The chain growth pre-exponential factor at the target active site is determined in the following manner: Obtain the chain growth pre-exponential factor of a single active site, the mass ratio of polymer generated at the target active site to the total polymer generated at all active sites, and the number of active sites of the new catalyst; The chain growth pre-exponential factor at the target active site is determined based on the chain growth pre-exponential factor of the single active site, the mass ratio, and the number of active sites of the new catalyst.

8. A device for adjusting dynamic parameters, characterized in that, include: The first determining module is used to determine the heat of polymerization in the polymerization reactor based on the heat balance formula, wherein the heat balance formula includes the mass flow rate of the inlet material of the polymerization reactor, the mass enthalpy of the inlet material of the polymerization reactor, the heat carried away by the cooling medium in the polymerization reactor, the mass flow rate of the outlet material of the polymerization reactor, the mass enthalpy of the outlet material of the polymerization reactor, and the heat of polymerization. The second determining module is used to determine the activity of the target catalyst based on the polymerization heat and the feed amount of the target catalyst in the polymerization reactor, wherein the target catalyst includes a new catalyst and an old catalyst; The third determining module is used to determine the activity of the mixed catalyst in the polymerization reactor based on the activity of the target catalyst. The adjustment module is used to adjust the pre-exponential factor in the kinetic parameters based on the activity of the mixed catalyst.

9. An electronic device, characterized in that, include: Memory, used to store program instructions; A processor, connected to the memory, is configured to execute program instructions to perform the following functions: determining the heat of polymerization in the polymerization reactor based on a heat balance formula, wherein the heat balance formula includes the mass flow rate of the inlet material of the polymerization reactor, the mass enthalpy of the inlet material of the polymerization reactor, the heat carried away by the cooling medium in the polymerization reactor, the mass flow rate of the outlet material of the polymerization reactor, the mass enthalpy of the outlet material of the polymerization reactor, and the heat of polymerization; determining the activity of the target catalyst based on the heat of polymerization and the feed rate of the target catalyst in the polymerization reactor, wherein the target catalyst includes a new catalyst and a used catalyst; determining the activity of the mixed catalyst in the polymerization reactor based on the activity of the target catalyst; and adjusting the pre-exponential factor in the kinetic parameters based on the activity of the mixed catalyst.

10. A non-volatile storage medium, characterized in that, The non-volatile storage medium includes a stored computer program, wherein the device containing the non-volatile storage medium executes the method for adjusting the dynamic parameters according to any one of claims 1 to 7 by running the computer program.

Images