CN107908110B - Tubular reactor dynamic optimization system based on control grid refinement - Google Patents
Tubular reactor dynamic optimization system based on control grid refinement Download PDFInfo
- Publication number
- CN107908110B CN107908110B CN201711116200.3A CN201711116200A CN107908110B CN 107908110 B CN107908110 B CN 107908110B CN 201711116200 A CN201711116200 A CN 201711116200A CN 107908110 B CN107908110 B CN 107908110B
- Authority
- CN
- China
- Prior art keywords
- flow rate
- tubular reactor
- module
- rate control
- concentration
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 claims abstract description 52
- 238000011217 control strategy Methods 0.000 claims abstract description 40
- 230000000875 corresponding effect Effects 0.000 claims abstract description 27
- 239000000463 material Substances 0.000 claims abstract description 22
- 238000004519 manufacturing process Methods 0.000 claims abstract description 15
- 238000006243 chemical reaction Methods 0.000 claims abstract description 9
- 239000013598 vector Substances 0.000 claims description 35
- 238000004364 calculation method Methods 0.000 claims description 31
- 230000010354 integration Effects 0.000 claims description 27
- 239000000376 reactant Substances 0.000 claims description 5
- 238000007670 refining Methods 0.000 claims description 5
- 230000009191 jumping Effects 0.000 claims description 3
- 238000009412 basement excavation Methods 0.000 abstract description 3
- 239000000047 product Substances 0.000 description 34
- 238000007796 conventional method Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000006227 byproduct Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000012847 fine chemical Substances 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000006116 polymerization reaction Methods 0.000 description 1
- QQONPFPTGQHPMA-UHFFFAOYSA-N propylene Natural products CC=C QQONPFPTGQHPMA-UHFFFAOYSA-N 0.000 description 1
- 125000004805 propylene group Chemical group [H]C([H])([H])C([H])([*:1])C([H])([H])[*:2] 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/418—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM]
- G05B19/4185—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM] characterised by the network communication
- G05B19/41855—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM] characterised by the network communication by local area network [LAN], network structure
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/418—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM]
- G05B19/41865—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM] characterised by job scheduling, process planning, material flow
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/20—Pc systems
- G05B2219/25—Pc structure of the system
- G05B2219/25232—DCS, distributed control system, decentralised control unit
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Quality & Reliability (AREA)
- Manufacturing & Machinery (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
- Feedback Control In General (AREA)
Abstract
The invention discloses a tubular reactor dynamic optimization system based on control grid refinement, which consists of a tubular reactor, a flow rate sensor, an analog-to-digital converter, a field bus network, a DCS, a main control room flow rate and product concentration display, a digital-to-analog converter at the flow rate control valve end and a flow rate control valve. A control room engineer specifies the concentration of each reaction material in the feeding of the tubular reactor, the DCS obtains a flow rate control strategy for maximizing the concentration of a target product at the tail end of the tubular reactor through a control grid refinement optimization method, the flow rate control strategy is converted into an opening instruction of a flow rate control valve and is sent to a digital-to-analog converter at the end of the flow rate control valve through a field bus network, the flow rate control valve executes corresponding actions, a flow rate sensor collects the flow rate of the tubular reactor in real time and sends the flow rate to the DCS, and the control room engineer can master the production process at any time. The invention can maximize the concentration of the target product at the tail end of the tubular reactor and realize the potential excavation and efficiency improvement.
Description
Technical Field
The invention relates to the field of reactor control, in particular to a tubular reactor dynamic optimization system based on control grid refinement. The flow rate of the tubular reactor can be automatically controlled to maximize the concentration of the target product, thereby improving the production efficiency of the reactor.
Background
The tubular reactor belongs to a plug flow reactor, and is widely applied to petrochemical industry, fine chemical industry and other industrial production, such as propylene polymerization production.
For a tubular reactor with a fixed tube length, when material parameters such as a raw material quota ratio, concentration and the like are determined, a key factor influencing the product yield is the flow velocity of reactants, namely the flow rate. Because the production process requirements of different products are different, the automatic flow rate control of the tubular reactor according to the production process requirements has important significance.
At present, the control method of the domestic tubular reactor rarely adopts a dynamic optimization theory and a corresponding method, parameters in the controller are often set according to the existing experience, and the production efficiency needs to be further improved. The product concentration of the tubular reactor after the dynamic optimization method is adopted can be further improved, and the potential excavation and efficiency improvement are realized.
Disclosure of Invention
In order to improve the concentration of the target product of the tubular reactor, the invention provides a tubular reactor dynamic optimization system based on control grid refinement.
The purpose of the invention is realized by the following technical scheme: the tubular reactor dynamic optimization system based on control grid refinement can automatically control the flow rate of the tubular reactor to maximize the concentration of a target product at the tail end of the tubular reactor. The method is characterized in that: the system consists of a tubular reactor 11, a flow rate sensor 12, an analog-to-digital converter 13, a field bus network 14, DCS15, a main control room flow rate and product concentration display 16, a digital-to-analog converter 17 at the end of a flow rate control valve and a flow rate control valve 18. A control room engineer specifies the concentration of each reaction material in the feeding of the tubular reactor, the DCS obtains a flow rate control strategy for maximizing the concentration of a target product at the tail end of the tubular reactor through a control grid refinement optimization method, the flow rate control strategy is converted into an opening instruction of a flow rate control valve and is sent to a digital-to-analog converter at the end of the flow rate control valve through a field bus network, the flow rate control valve executes corresponding actions, a flow rate sensor collects the flow rate of the tubular reactor in real time and sends the flow rate to the DCS, and the control room engineer can master the production process at any time. The operation process of the system comprises the following steps:
step A1: a control room engineer specifies the concentration of each reactant in the feed to the tubular reactor, the target product for which the concentration needs to be maximized, and the flow rate control requirements;
step A2: the DCS executes an internal control grid refinement optimization method to obtain a flow rate control strategy for maximizing the concentration of the target product at the tail end of the tubular reactor;
step A3: the DCS converts the flow rate control strategy obtained through calculation into an opening instruction of the flow rate control valve, and sends the opening instruction to a digital-to-analog converter at the end of the flow rate control valve through a field bus network, so that the flow rate control valve executes corresponding action according to the received control instruction;
step A4: the flow rate sensor of the tubular reactor collects the flow rate in real time, the flow rate is returned to the DCS through the analog-to-digital converter by using a field bus network and is displayed in the main control room, so that an engineer in the control room can master the production process at any time.
The DCS comprises an information acquisition module, an initialization module, a control grid refinement module, an ODE solving module, a gradient calculation module, a Non-linear Programming (NLP) problem solving module, a refinement convergence judgment module and a control instruction output module. The information acquisition module comprises three submodules of feeding concentration state acquisition, target product acquisition and flow rate control requirement acquisition, and the NLP problem solving module comprises three submodules of optimizing direction calculation, optimizing step length calculation and NLP convergence judgment.
The production process of a tubular reactor can be described as:
wherein t represents a change in the tube length direction; u (t) represents the flow rate; x (t) represents the concentration of material in the tubular reactor varying along the length of the tube; f (-) is a system of differential equations established based on material balance, energy balance, etc. From this description, it can be seen that the process of producing the target product in the tubular reactor can be represented mathematically by a set of differential equations.
To maximize the concentration of the target product in the tubular reactor, x1(t) represents the concentration of the target product as it varies along the length of the tube, the final expression for this problem is:
wherein, t0Denotes the feed inlet of the tubular reactor, tfRepresenting the end of the tubular reactor and J represents the objective function to be maximized. The problem is essentially a dynamic optimization problem. However, conventional methods for solving such problems haveLow efficiency and poor precision, and is difficult to meet the requirement of high efficiency in actual production.
The technical scheme adopted by the invention for solving the technical problems is as follows: a control grid refinement optimization method is integrated in the DCS, and a set of dynamic optimization control system is constructed on the basis of the control grid refinement optimization method. The complete structure of the system comprises a tubular reactor 21, a flow rate sensor 22, an analog-to-digital converter 23, a field bus network 24, a DCS25, a main control room flow rate and product concentration display 26, a digital-to-analog converter 27 at the flow rate control valve end, and a flow rate control valve 28.
The operation process of the system comprises the following steps:
step A1: a control room engineer specifies the concentration of each reactant in the feed to the tubular reactor, the target product for which the concentration needs to be maximized, and the flow rate control requirements;
step A2: the DCS executes an internal control grid refinement optimization method to obtain a flow rate control strategy for maximizing the concentration of the target product at the tail end of the tubular reactor;
step A3: the DCS converts the flow rate control strategy obtained through calculation into an opening instruction of the flow rate control valve, and sends the opening instruction to a digital-to-analog converter at the end of the flow rate control valve through a field bus network, so that the flow rate control valve executes corresponding action according to the received control instruction;
step A4: the flow rate sensor of the tubular reactor collects the flow rate in real time, the flow rate is returned to the DCS through the analog-to-digital converter by using a field bus network and is displayed in the main control room, so that an engineer in the control room can master the production process at any time.
The DCS comprises an information acquisition module, an initialization module, a control grid refinement module, an ODE solving module, a gradient calculation module, a Non-linear Programming (NLP) problem solving module, a refinement convergence judgment module and a control instruction output module. The information acquisition module comprises three submodules of feeding concentration state acquisition, target product acquisition and flow rate control requirement acquisition, and the NLP problem solving module comprises three submodules of optimizing direction calculation, optimizing step length calculation and NLP convergence judgment.
In order to obtain a flow rate control strategy for maximizing the concentration of a target product at the tail end of the tubular reactor, the control grid fine optimization method executed by the DCS comprises the following operation steps:
step B1: the information acquisition module 31 acquires the initial concentration of the reaction materials, the target product of which the concentration needs to be maximized and the flow rate control requirement specified by an engineer;
step B2: the initialization module 32 starts to run, parameterizes with a piecewise constant, sets the number of segments N of the pipe length and the corresponding control grid toInitial guess of parameterized vector of flow rate control strategySetting the computational accuracy tol of NLP problem1Convergence accuracy tol of sum mesh refinement2The number of iterations k1And number of refinements k2Setting zero;
step B3: when k is2When 0, execute step B4; otherwise, the control mesh is refined by the control mesh refinement module 33Refining to obtain new control gridAnd corresponding parameterized vector
Step B4: obtaining the material concentration of the iteration through the ODE solving module 34And an objective function value
Step B5: the book is acquired by the gradient calculation module 35Sub-iterative gradient informationWhen k is1When 0, directly executing the step B7 by skipping the step B6;
step B6: the NLP problem solving module 36 operates to judge convergence through the NLP convergence judging module, if soObjective function value from last iterationThe absolute value of the difference is less than the accuracy tol1If yes, determining that the convergence is satisfied, and executing step B9; if convergence is not satisfied, continuing to step B7;
step B8: the NLP problem solving module 36 obtains a ratio by calculating the optimizing direction and the optimizing step size using the objective function values and gradient information obtained in steps B4 and B5More optimal new flow rate control strategyAfter the step is executed, jumping to the step B4 again;
step B9: the refinement convergence determination module 37 operates to recordWhen k is2When the value is 0, the step B10 is executed, otherwise, the judgment is madeAnd the last refined objective function valueWhether the absolute value of the difference is less than the accuracy tol2If yes, the convergence is judged to be satisfied, the flow rate control strategy of the iteration is converted into an opening instruction of the flow rate control valve to be output, otherwise, the convergence is not satisfied, and the refinement times k are set2:=k2+1, continue to execute step B3 until the refinement convergence determination module is satisfied.
The control grid refinement module is realized by adopting the following steps:
step C1: calculating the left slope at a grid node by the following equationAnd right slope(k=1,…,N-1):
Wherein u iskThe k-th component, t, of a parameterized vector u representing a flow rate control strategykRepresents ukAnd uk+1A mesh node in between.
Step C2: if the mesh node tkIf the left and right slopes meet the following requirements, the node is removed from the grid:
wherein epsiloneIs a small positive real number. Grid node tkAfter removal of ukAnd uk+1The corresponding grids are merged into a new grid, and the parameters on the new grid are updated to (u)k+uk+1)/2。
Step C3: if the mesh node tkThe left slope at (b) satisfies:
wherein epsiloniIs one greater than epsilonePositive real number of (1), then at [ tk-1,tk]Inserting grid nodes upwards; if the mesh node tkThe right slope at (b) satisfies:
then at [ tk-1,tk]And inserting grid nodes. In practical application, the number of the added nodes can be freely set according to the absolute value of the left slope and the right slope.
Step C4: and generating a new control grid and a corresponding parameterization vector according to the nodes removed and inserted in the steps C2 and C3.
The ODE solving module adopts a four-step Runge-Kutta method, and the calculation formula is as follows:
wherein t represents a change in the longitudinal direction of the tube, and tiDenotes the integration time, t, selected by the Runge-Kutta methodi+1Indicating that it is at time tiThe integration step h is the difference between any two adjacent integration moments, x (t)i) Denotes the distance t from the inlet of the tubular reactoriThe material concentration, F (-) is a function describing a state differential equation, and K1, K2, K3 and K4 respectively represent function values of 4 nodes in the Runge-Kutta method integration process.
The gradient calculation module adopts an accompanying method:
step D1: let λ (t) be the co-modal vector, whose value is determined by the adjoint equation:
wherein, tfDenotes the end of a tubular reactor, H denotes the Hamiltonian, and H ═ L + λ (t)TF, L is the integral term of the objective function, Φ x (t)f)]Is the steady state term of the objective function.
Step D2: for the adjoint equation, a four-step Runge-Kutta method is adopted to obtain the value of the co-modal vector lambda (t) at each integration moment, and the calculation formula is as follows:
wherein t represents a change in the longitudinal direction of the tube, and tiSolving the selected integration time, t, in the module for ODEi+1Indicating that it is at time tiThe latter integration time, and ti+1=ti+ h, h is the integration step, and Q1, Q2, Q3 and Q4 respectively represent the function values of 4 nodes in the integration process of the Runge-Kutta method.
Step D3: based on the obtained value of the co-modal vector λ (t), gradient information is obtained by the following formula
The NLP problem solving module is realized by adopting the following steps:
step E1: if it is notObjective function value from last iterationIs less than the accuracy tol1If yes, judging that the convergence is satisfied, and returning to the flow rate control strategy obtained by the iteration; if convergence is not satisfied, continuing to step E2;
step E3: control strategy for flow rateAs a point in vector space, denoted as P1,P1The corresponding objective function value is
Step E4: from point P1Starting from the selected NLP algorithm and point P1Information of the gradient ofConstructing a direction of optimization in vector spaceAnd step size
Step E5, by formulaConstructing correspondences in vector spaceAnother point P of2So that P is2Corresponding objective function valueRatio ofPreferably, wherein I isA vector of the same dimension.
The invention has the following beneficial effects: the tubular reactor dynamic optimization system based on control grid refinement can calculate the optimal flow rate control strategy of the tubular reactor, can adapt to the optimal control curve of the problem, particularly find the discontinuous point of the problem, and can obtain higher precision; after the adaptive strategy is adopted, the initial estimation value of the next optimal control curve is the optimal curve of the current iteration, so that the faster convergence speed can be obtained, and the calculation time required for obtaining the optimal control strategy is reduced. The invention can maximize the concentration of the target product in the tubular reactor and realize excavation and synergy.
Drawings
FIG. 1 is a functional schematic of the present invention;
FIG. 2 is a schematic structural view of the present invention;
FIG. 3 is a block diagram of the DCS internal module of the present invention;
Detailed Description
As shown in fig. 1, the production process of a tubular reactor can be described as:
wherein t represents a change in the tube length direction; u (t) represents the flow rate; x (t) represents the concentration of material in the tubular reactor varying along the length of the tube; f (-) is a system of differential equations established based on material balance, energy balance, etc. From this description, it can be seen that the process of producing the target product in the tubular reactor can be represented mathematically by a set of differential equations.
To maximize the concentration of the target product in the tubular reactor, x1(t) represents the concentration of the target product as it varies along the length of the tube, the final expression for this problem is:
wherein, t0Denotes the feed inlet of the tubular reactor, tfRepresenting the end of the tubular reactor and J represents the objective function to be maximized. The problem is essentially a dynamic optimization problem. However, the conventional method has the defects of low efficiency and poor precision in solving the problems, and is difficult to meet the requirement of high efficiency in actual production.
The technical scheme adopted by the invention for solving the technical problems is as follows: a control grid refinement optimization method is integrated in the DCS, and a set of dynamic optimization control system is constructed on the basis of the control grid refinement optimization method. The complete structure of the system is shown in fig. 2, and comprises a tubular reactor 21, a flow rate sensor 22, an analog-to-digital converter 23, a fieldbus network 24, a DCS25, a main control room flow rate and product concentration display 26, a digital-to-analog converter 27 at the flow rate control valve end, and a flow rate control valve 28.
The operation process of the system comprises the following steps:
step A1: a control room engineer specifies the concentration of each reactant in the feed to the tubular reactor, the target product for which the concentration needs to be maximized, and the flow rate control requirements;
step A2: the DCS executes an internal control grid refinement optimization method to obtain a flow rate control strategy for maximizing the concentration of the target product at the tail end of the tubular reactor;
step A3: the DCS converts the flow rate control strategy obtained through calculation into an opening instruction of the flow rate control valve, and sends the opening instruction to a digital-to-analog converter at the end of the flow rate control valve through a field bus network, so that the flow rate control valve executes corresponding action according to the received control instruction;
step A4: the flow rate sensor of the tubular reactor collects the flow rate in real time, the flow rate is returned to the DCS through the analog-to-digital converter by using a field bus network and is displayed in the main control room, so that an engineer in the control room can master the production process at any time.
The DCS comprises an information acquisition module, an initialization module, a control grid refinement module, an ODE solving module, a gradient calculation module, a Non-linear Programming (NLP) problem solving module, a refinement convergence judgment module and a control instruction output module. The information acquisition module comprises three submodules of feeding concentration state acquisition, target product acquisition and flow rate control requirement acquisition, and the NLP problem solving module comprises three submodules of optimizing direction calculation, optimizing step length calculation and NLP convergence judgment.
In order to obtain a flow rate control strategy for maximizing the concentration of a target product at the tail end of the tubular reactor, the control grid fine optimization method executed by the DCS comprises the following operation steps:
step B1: the information acquisition module 31 acquires the initial concentration of the reaction materials, the target product of which the concentration needs to be maximized and the flow rate control requirement specified by an engineer;
step B2: the initialization module 32 starts to run, parameterizes with a piecewise constant, sets the number of segments N of the pipe length and the corresponding control grid toInitial guess of parameterized vector of flow rate control strategySetting the computational accuracy tol of NLP problem1Convergence accuracy tol of sum mesh refinement2The number of iterations k1And number of refinements k2Setting zero;
step B3: when k is2When 0, execute step B4; otherwise, by controllingGrid refinement module 33 pairs control gridsRefining to obtain new control gridAnd corresponding parameterized vector
Step B4: obtaining the material concentration of the iteration through the ODE solving module 34And an objective function value
Step B5: obtaining the gradient information of the iteration through the gradient calculation module 35When k is1When 0, directly executing the step B7 by skipping the step B6;
step B6: the NLP problem solving module 36 operates to judge convergence through the NLP convergence judging module, if soObjective function value from last iterationThe absolute value of the difference is less than the accuracy tol1If yes, determining that the convergence is satisfied, and executing step B9; if convergence is not satisfied, continuing to step B7;
step B8: the NLP problem solving module 36 obtains a ratio by calculating the optimizing direction and the optimizing step size using the objective function values and gradient information obtained in steps B4 and B5More optimal new flow rate control strategyAfter the step is executed, jumping to the step B4 again;
step B9: the refinement convergence determination module 37 operates to recordWhen k is2When the value is 0, the step B10 is executed, otherwise, the judgment is madeAnd the last refined objective function valueWhether the absolute value of the difference is less than the accuracy tol2If yes, the convergence is judged to be satisfied, the flow rate control strategy of the iteration is converted into an opening instruction of the flow rate control valve to be output, otherwise, the convergence is not satisfied, and the refinement times k are set2:=k2+1, continue to execute step B3 until the refinement convergence determination module is satisfied.
The control grid refinement module is realized by adopting the following steps:
step C1: calculating the left slope at a grid node by the following equationAnd right slope(k=1,…,N-1):
Wherein u iskThe k-th component, t, of a parameterized vector u representing a flow rate control strategykRepresents ukAnd uk+1A mesh node in between.
Step C2: if the mesh node tkIf the left and right slopes meet the following requirements, the node is removed from the grid:
wherein epsiloneIs a small positive real number. Grid node tkAfter removal of ukAnd uk+1The corresponding grids are merged into a new grid, and the parameters on the new grid are updated to (u)k+uk+1)/2。
Step C3: if the mesh node tkThe left slope at (b) satisfies:
wherein epsiloniIs one greater than epsilonePositive real number of (1), then at [ tk-1,tk]Inserting grid nodes upwards; if the mesh node tkThe right slope at (b) satisfies:
then at [ tk-1,tk]And inserting grid nodes. In practical application, the number of the added nodes can be freely set according to the absolute value of the left slope and the right slope.
Step C4: and generating a new control grid and a corresponding parameterization vector according to the nodes removed and inserted in the steps C2 and C3.
The ODE solving module adopts a four-step Runge-Kutta method, and the calculation formula is as follows:
wherein t represents a change in the longitudinal direction of the tube, and tiDenotes the integration time, t, selected by the Runge-Kutta methodi+1Indicating that it is at time tiThe integration step h is the difference between any two adjacent integration moments, x (t)i) Denotes the distance t from the inlet of the tubular reactoriThe material concentration, F (-) is a function describing a state differential equation, and K1, K2, K3 and K4 respectively represent function values of 4 nodes in the Runge-Kutta method integration process.
The gradient calculation module adopts an accompanying method:
step D1: let λ (t) be the co-modal vector, whose value is determined by the adjoint equation:
wherein, tfDenotes the end of a tubular reactor, H denotes the Hamiltonian, and H ═ L + λ (t)TF, L is the integral term of the objective function, Φ x (t)f)]Is the steady state term of the objective function.
Step D2: for the adjoint equation, a four-step Runge-Kutta method is adopted to obtain the value of the co-modal vector lambda (t) at each integration moment, and the calculation formula is as follows:
wherein t represents a change in the longitudinal direction of the tube, and tiSolving the selected integration time, t, in the module for ODEi+1Indicating that it is at time tiThe latter integration time, and ti+1=ti+ h, h is the integration step, Q1, Q2, Q3 and Q4 respectively represent the integration process of Runge-Kutta methodFunction values of 4 nodes.
Step D3: based on the obtained value of the co-modal vector λ (t), gradient information is obtained by the following formula
The NLP problem solving module is realized by adopting the following steps:
step E1: if it is notObjective function value from last iterationIs less than the accuracy tol1If yes, judging that the convergence is satisfied, and returning to the flow rate control strategy obtained by the iteration; if convergence is not satisfied, continuing to step E2;
step E3: control strategy for flow rateAs a point in vector space, denoted as P1,P1The corresponding objective function value is
Step E4: from point P1Starting from the selected NLP algorithm and point P1Information of the gradient ofConstructing a direction of optimization in vector spaceAnd step size
Example 1
Parallel reactions occur in a tubular reactor: a → B and A → C, reaction rate constants are k1And k2B is the target product, C is a by-product, and the aim is to control the choked flow rate so that product B is in the reactorThe concentration at the end was maximal. By xA(t) and xB(t) represents the concentration of materials A and B along the length of the tube, respectively, and u (t) k1L/v (t) is a physical quantity directly related to the choked flow rate, where L is the reactor length and v (t) is the choked flow rate. Finally, the optimal control problem can be simplified as:
wherein J represents the concentration of the material B to be maximized, tfRepresenting the end of the tubular reactor, t represents the change in length along the length of the tube, x (0) the initial concentration of materials a and B at the inlet of the tubular reactor. In order to obtain a flow rate control strategy for maximizing the concentration of a target product B at the tail end of the tubular reactor, a DCS operation control grid fine optimization method is adopted, the operation process is shown in FIG. 3, and the execution steps are as follows:
step F1: the information acquisition module 31 obtains the initial concentration x (0) ═ 10 of the reaction material specified by the engineer]TThe target product B with the concentration needing to be maximized and the flow rate control requirement of 0 to u (t) to 5;
step F2: the initialization module 32 starts to run, parameterizes with a piecewise constant, sets the number of segments of the pipe length to 8, and sets the corresponding control gridInitial guess u of parameterized vector for uniform partitioning, flow rate control strategy(k)To 0.5, the calculation accuracy tol of the NLP problem is set1Convergence accuracy tol of sum mesh refinement2Are respectively 10-7And 10-6The number of iterations k1And number of refinements k2Setting zero;
step F3: when k is2When 0, step F4 is executed; otherwise, the control mesh is refined by the control mesh refinement module 33Refining to obtain new control gridAnd corresponding parameterized vector
Step F4: obtaining the material concentration of the iteration through the ODE solving module 34And an objective function value
Step F5: obtaining the gradient information of the iteration through the gradient calculation module 35When k is1Step F7 is directly performed skipping step F6 when 0;
step F6: the NLP problem solving module 36 operates to judge convergence through the NLP convergence judging module, if soObjective function value from last iterationThe absolute value of the difference is less than the accuracy tol1If yes, determining that the convergence is satisfied, and executing step F9; if the convergence is not satisfied, continuing to perform step F7;
step F8: the NLP problem solving module 36 uses the objective function values and gradient information obtained in steps F4 and F5 by calculating an optimizing methodStep size of direction and optimization, ratio of gainMore optimal new flow rate control strategyAfter the step is completed, the step goes to step F4 again;
step F9: the refinement convergence determination module 37 operates to recordWhen k is2When 0, the process proceeds to step F10, otherwise, it is judgedAnd the last refined objective function valueWhether the absolute value of the difference is less than the accuracy tol2If yes, the convergence is judged to be satisfied, the flow rate control of the iteration is converted into the opening instruction of the mixed flow rate valve and output, otherwise, the convergence is not satisfied, and the refining times k are set2:=k2+1, continue to execute step F3 until the refinement convergence determination module is satisfied.
And finally, the DCS converts the flow rate control strategy obtained by the control grid refinement optimization method into an opening instruction of the flow rate control valve, and sends the opening instruction to a digital-to-analog converter at the end of the flow rate control valve through the field bus network, so that the flow rate control valve executes corresponding actions according to the received control instruction, meanwhile, the flow rate sensor is used for collecting the flow rate distributed along the length of the pipe of the tubular reactor in real time, and the flow rate is returned to the DCS through the field bus network after passing through the analog-to-digital converter and is displayed in the main control chamber.
The foregoing is a more detailed description of the invention in connection with specific preferred embodiments and is not intended to limit the practice of the invention to these embodiments. For those skilled in the art to which the invention pertains, several simple deductions or substitutions may be made without departing from the inventive concept, which should be construed as falling within the scope of the present invention.
Claims (1)
1. The tubular reactor dynamic optimization system based on control grid refinement can automatically control the flow rate of the tubular reactor so as to maximize the concentration of a target product at the tail end of the tubular reactor; the method is characterized in that: the system is composed of a tubular reactor (11), a flow rate sensor (12), an analog-to-digital converter (13), a field bus network (14), a distributed control system (15), a main control room flow rate and product concentration display (16), a digital-to-analog converter (17) at the end of a flow rate control valve, and a flow rate control valve (18); a control room engineer specifies the concentration of each reaction material in the feeding of the tubular reactor, the distributed control system obtains a flow rate control strategy which enables the concentration of a target product to be maximum at the tail end of the tubular reactor through a control grid refinement optimization method, converts the flow rate control strategy into an opening instruction of a flow rate control valve, sends the opening instruction to a digital-to-analog converter at the end of the flow rate control valve through a field bus network, enables the flow rate control valve to execute corresponding actions, and enables a flow rate sensor to acquire the flow rate of the tubular reactor in real time and send the flow rate back to the distributed control system, so that the control room engineer can master the; the operation process of the system comprises the following steps:
step A1: a control room engineer specifies the concentration of each reactant in the feed to the tubular reactor, the target product for which the concentration needs to be maximized, and the flow rate control requirements;
step A2: executing an internal control grid refinement optimization method by the distributed control system to obtain a flow rate control strategy for maximizing the concentration of a target product at the tail end of the tubular reactor;
step A3: the distributed control system converts the flow rate control strategy obtained by calculation into an opening instruction of the flow rate control valve, and sends the opening instruction to a digital-to-analog converter at the end of the flow rate control valve through a field bus network, so that the flow rate control valve executes corresponding action according to the received control instruction;
step A4: the flow rate sensor of the tubular reactor collects the flow rate in real time, the flow rate passes through the analog-to-digital converter and is returned to the distributed control system by a field bus network, and the flow rate is displayed in the main control room, so that an engineer in the control room can master the production process at any time;
the distributed control system comprises an information acquisition module, an initialization module, a control grid refinement module, an ODE solving module, a gradient calculation module, a Non-linear Programming (NLP) problem solving module, a refinement convergence judgment module and a control instruction output module; the information acquisition module comprises three submodules of feeding concentration state acquisition, target product acquisition and flow rate control requirement acquisition, and the NLP problem solving module comprises three submodules of optimizing direction calculation, optimizing step length calculation and NLP convergence judgment;
in order to obtain a flow rate control strategy for maximizing the concentration of the target product at the tail end of the tubular reactor, the control grid fine optimization method executed by the distributed control system comprises the following operation steps:
step B1: the information acquisition module (31) acquires the initial concentration of the reaction materials, target products needing to maximize the concentration and flow rate control requirements specified by a control room engineer;
step B2: the initialization module (32) starts to operate, parameterizes by adopting a piecewise constant, and sets the number N of segments of the pipe length and the corresponding control grid asInitial guess of parameterized vector of flow rate control strategySetting the computational accuracy tol of NLP problem1Convergence accuracy tol of sum mesh refinement2The number of iterations k1And number of refinements k2Setting zero;
step B3: when k is2When 0, execute step B4; otherwise, the control grid is refined by a control grid refinement module (33)Refining to obtain new control gridAnd corresponding parameterized vector
Step B4: obtaining the material concentration of the iteration through an ODE solving module (34)And an objective function value
Step B5: obtaining gradient information of the iteration through a gradient calculation module (35)When k is1When 0, directly executing the step B7 by skipping the step B6;
step B6: the NLP problem solving module (36) operates, the convergence judgment is carried out through the NLP convergence judgment module, and if the convergence judgment module is used, the convergence judgment is carried outObjective function value from last iterationThe absolute value of the difference is less than the accuracy tol1If yes, determining that the convergence is satisfied, and executing step B9; if convergence is not satisfied, continuing to step B7;
step B8: NLP problem solving modelA block (36) obtains a ratio by calculating a seek direction and a seek step length using the objective function values and gradient information obtained in steps B4 and B5More optimal new flow rate control strategyAfter the step is executed, jumping to the step B4 again;
step B9: the refined convergence judgment module (37) runs and recordsWhen k is2When the value is 0, the step B10 is executed, otherwise, the judgment is madeAnd the last refined objective function valueWhether the absolute value of the difference is less than the accuracy tol2If yes, the convergence is judged to be satisfied, the flow rate control strategy of the iteration is converted into an opening instruction of the flow rate control valve to be output, otherwise, the convergence is not satisfied, and the refinement times k are set2:=k2+1, continuing to execute the step B3 until the refinement convergence judgment module is satisfied;
the control grid refinement module is realized by adopting the following steps:
Wherein u iskParameterized vector representing flow rate control strategyThe k component of (a), tkRepresents ukAnd uk+1A mesh node in between;
step C2: if the mesh node tkIf the left and right slopes meet the following requirements, the node is removed from the grid:
wherein epsiloneIs a small positive real number; grid node tkAfter removal of ukAnd uk+1The corresponding grids are merged into a new grid, and the parameters on the new grid are updated to (u)k+uk+1)/2;
Step C3: if the mesh node tkThe left slope at (b) satisfies:
wherein epsiloniIs one greater than epsilonePositive real number of (1), then at [ tk-1,tk]Inserting grid nodes upwards; if the mesh node tkThe right slope at (b) satisfies:
then at [ tk-1,tk]Inserting grid nodes upwards; in actual application, the number of the added nodes can be freely set according to the absolute value of the left slope and the right slope;
step C4: generating a new control grid and corresponding parameterized vectors according to the nodes removed and inserted in the steps C2 and C3;
the ODE solving module adopts a four-step Runge-Kutta method, and the calculation formula is as follows:
wherein t represents a change in the longitudinal direction of the tube, and tiDenotes the integration time, t, selected by the Runge-Kutta methodi+1Indicating that it is at time tiThe integration step h is the difference between any two adjacent integration moments, x (t)i) Denotes the distance t from the inlet of the tubular reactoriThe material concentration, F (-) is a function describing a state differential equation, and K1, K2, K3 and K4 respectively represent function values of 4 nodes in the Runge-Kutta method integration process;
the gradient calculation module adopts an accompanying method:
step D1: let λ (t) be the co-modal vector, whose value is determined by the adjoint equation:
wherein, tfDenotes the end of a tubular reactor, H denotes the Hamiltonian, and H ═ L + λ (t)TF, L is the integral term of the objective function, Φ x (t)f)]A steady state term that is an objective function;
step D2: for the adjoint equation, a four-step Runge-Kutta method is adopted to obtain the value of the co-modal vector lambda (t) at each integration moment, and the calculation formula is as follows:
wherein t represents a change in the longitudinal direction of the tube, and tiSolving the selected integration time, t, in the module for ODEi+1Indicating that it is at time tiThe latter integration time, and ti+1=ti+ h, h is the integration step length, and Q1, Q2, Q3 and Q4 respectively represent function values of 4 nodes in the Runge-Kutta method integration process;
step D3: based on the obtained value of the co-modal vector λ (t), gradient information is obtained by the following formula
the NLP problem solving module is realized by adopting the following steps:
step E1: performing steps B6-B7;
step E2: control strategy for flow rateAs a point in vector space, denoted as P1,P1The corresponding objective function value is
Step E3: from point P1Starting from the selected NLP algorithm and point P1Information of the gradient ofConstructing a direction of optimization in vector spaceAnd step size
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711116200.3A CN107908110B (en) | 2017-11-13 | 2017-11-13 | Tubular reactor dynamic optimization system based on control grid refinement |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711116200.3A CN107908110B (en) | 2017-11-13 | 2017-11-13 | Tubular reactor dynamic optimization system based on control grid refinement |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107908110A CN107908110A (en) | 2018-04-13 |
CN107908110B true CN107908110B (en) | 2020-03-03 |
Family
ID=61845175
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711116200.3A Expired - Fee Related CN107908110B (en) | 2017-11-13 | 2017-11-13 | Tubular reactor dynamic optimization system based on control grid refinement |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107908110B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110119085B (en) * | 2019-04-03 | 2022-03-15 | 杭州电子科技大学 | Manutec R3 type industrial robot dynamic optimization system |
CN111308893B (en) * | 2020-03-05 | 2021-06-18 | 江南大学 | Non-uniform grid optimization method for bioreactor fed-batch regulation |
CN116382209B (en) * | 2023-04-19 | 2023-10-03 | 扬州市管件厂有限公司 | Process optimization method and system for seamless elbow machining |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7536240B2 (en) * | 2005-07-22 | 2009-05-19 | Ut-Battelle, Llc | Real power regulation for the utility power grid via responsive loads |
CN101464664B (en) * | 2009-01-09 | 2010-12-01 | 浙江工业大学 | Batch reactor optimal control method based on single population and pre-crossed differential evolution algorithm |
CN101763086A (en) * | 2009-12-29 | 2010-06-30 | 浙江大学 | High-efficiency system and method for dynamic optimization of industrial process |
CN101763083A (en) * | 2009-12-29 | 2010-06-30 | 浙江大学 | An effective industrial process dynamic optimization system and method capable of controlling variable parameterization |
CN103324086B (en) * | 2013-06-08 | 2016-01-20 | 浙江大学 | A kind of batch reactor control system based on accurately punishment optimization |
CN103309234B (en) * | 2013-06-08 | 2015-12-09 | 浙江大学 | A kind of batch reactor control system optimized based on orthogonal configuration |
-
2017
- 2017-11-13 CN CN201711116200.3A patent/CN107908110B/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN107908110A (en) | 2018-04-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107908110B (en) | Tubular reactor dynamic optimization system based on control grid refinement | |
CN103439887B (en) | PI attitude conirol method and system optimal lower order system ITAE | |
CN108520325B (en) | Integrated life prediction method based on accelerated degradation data in variable environment | |
Peroni et al. | Optimal control of a fed-batch bioreactor using simulation-based approximate dynamic programming | |
CN108710299B (en) | Optimal control system and method for catalyst feed supplement of continuous stirred tank reactor | |
CN106406083A (en) | Fuzzy control coal dressing method | |
CN110309481B (en) | Model modeling method for farmland soil reaction dynamic process | |
CN107943119B (en) | Catalyst mixing dynamic optimization system based on control grid refinement | |
CN108919639A (en) | A kind of PID controller parameter best proportion method for establishing model | |
Patel et al. | Modified genetic algorithm using Box Complex method: Application to optimal control problems | |
CN102566426A (en) | Fractional order parameter adjustment controller algorithm of PI<alpha>D<beta> controller | |
CN110221540B (en) | Continuous stirred reactor system control method based on Hammerstein model | |
CN103177291A (en) | Variable-search-space ribonucleic acid (RNA) genetic algorithm modeling method for continuous stirred tank reactor | |
CN103324086A (en) | Batch reactor control system based on accurate punishment optimization | |
CN107807528A (en) | A kind of plug flow tubular reactor optimal control system based on adaptive congestion control algorithm node | |
CN109635330B (en) | Method for accurately and rapidly solving complex optimization control problem based on direct method | |
Schneider et al. | An iterative partition-based moving horizon estimator for large-scale linear systems | |
Würth et al. | On the numerical solution of discounted economic NMPC on infinite horizons | |
CN106647242A (en) | Multivariable PID controller parameter setting method | |
CN110794676A (en) | CSTR process nonlinear control method based on Hammerstein-Wiener model | |
CN107844056A (en) | A kind of plug flow tubular reactor optimal control system based on adaptive congestion control algorithm grid | |
CN107918277A (en) | A kind of plug flow tubular reactor optimal control system based on self adaptive control grid | |
CN113189871B (en) | Thermal power generating unit flexible operation-oriented coordination control strategy | |
CN107902555B (en) | A kind of shore container crane optimal control system of grid fining | |
CN107942967B (en) | Batch reactor optimal control system based on self-adaptive optimization control grid |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200303 Termination date: 20211113 |
|
CF01 | Termination of patent right due to non-payment of annual fee |