CN114185369A - Chlorine-hydrogen ratio dynamic balance control method in hydrogen chloride preparation process - Google Patents

Chlorine-hydrogen ratio dynamic balance control method in hydrogen chloride preparation process Download PDF

Info

Publication number
CN114185369A
CN114185369A CN202111324006.0A CN202111324006A CN114185369A CN 114185369 A CN114185369 A CN 114185369A CN 202111324006 A CN202111324006 A CN 202111324006A CN 114185369 A CN114185369 A CN 114185369A
Authority
CN
China
Prior art keywords
flow
chlorine
hydrogen
value
loop
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.)
Granted
Application number
CN202111324006.0A
Other languages
Chinese (zh)
Other versions
CN114185369B (en
Inventor
苏宝玉
孔庆儒
臧春华
臧博丰
马传山
康朔
陈香归
张帅杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenyang Huakong Technology Development Co ltd
Original Assignee
Shenyang Huakong Technology Development Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shenyang Huakong Technology Development Co ltd filed Critical Shenyang Huakong Technology Development Co ltd
Priority to CN202111324006.0A priority Critical patent/CN114185369B/en
Publication of CN114185369A publication Critical patent/CN114185369A/en
Application granted granted Critical
Publication of CN114185369B publication Critical patent/CN114185369B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D11/00—Control of flow ratio
    • G05D11/02—Controlling ratio of two or more flows of fluid or fluent material
    • G05D11/13—Controlling ratio of two or more flows of fluid or fluent material characterised by the use of electric means
    • G05D11/139—Controlling ratio of two or more flows of fluid or fluent material characterised by the use of electric means by measuring a value related to the quantity of the individual components and sensing at least one property of the mixture
    • C—CHEMISTRY; METALLURGY
    • C01—INORGANIC CHEMISTRY
    • C01B—NON-METALLIC ELEMENTS; COMPOUNDS THEREOF; METALLOIDS OR COMPOUNDS THEREOF NOT COVERED BY SUBCLASS C01C
    • C01B7/00—Halogens; Halogen acids
    • C01B7/01—Chlorine; Hydrogen chloride
    • C01B7/012—Preparation of hydrogen chloride from the elements

Landscapes

  • Chemical & Material Sciences (AREA)
  • Organic Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Inorganic Chemistry (AREA)
  • Flow Control (AREA)

Abstract

A chlorine-hydrogen ratio dynamic balance control method in the process of preparing hydrogen chloride belongs to the technical field of automatic control, and comprises the following steps: step S01: calculating the loop characteristics of the chlorine and hydrogen flows, and determining the reference values of the changes of the chlorine and hydrogen flows through a loop flow model; step S02: calculating an optimum load change rate based on the reference value of the flow rate in step S01 in combination with the output characteristics of the controller; step S03: calculating a load increment target value according to the load change rate in step S02, and performing variable load control according to the calculated change rate; step S04: calculating an optimal control output value, and adjusting the set values of the chlorine flow and the hydrogen flow according to the optimal control output value when the actual ratio of the chlorine flow to the hydrogen flow measurement value deviates from the target ratio; step S05: and calculating the locking characteristic of the regulating valve, and performing valve locking control according to the locking characteristic of the regulating valve when the chlorine flow and hydrogen flow loops are disturbed.

Description

Chlorine-hydrogen ratio dynamic balance control method in hydrogen chloride preparation process
Technical Field
The invention relates to the technical field of automatic control, in particular to a chlorine-hydrogen ratio dynamic balance control method in a hydrogen chloride preparation process.
Background
In the process for preparing hydrogen chloride, a hydrogen chloride synthesis furnace is an important process in the production process, and chlorine Cl is used in the process2And hydrogen H2Feeding the mixture into a synthesis furnace according to a certain proportion for combustion, and then producing hydrogen chloride HCl. Hydrogen chloride is generally used as a raw material for producing vinyl chloride VCM, so the quality of hydrogen chloride product has a great influence on the safety and stability of VCM production process, for example, chlorine cannot be overchlorided in hydrogen chloride in vinyl chloride production, if the overchloride reacts with acetylene in VCM conversion section, chloroethyne which is extremely explosive is generated, and explosion danger exists.
In addition, in the production process, the influence of variable load and some disturbance always exists, so that the manual control and the traditional ratio control method are difficult to meet the requirement that the actual ratio of the chlorine flow to the hydrogen flow tracks the expected ratio, and the chlorine-hydrogen ratio is easily mismatched, thereby bringing potential safety hazard to the production.
Disclosure of Invention
In order to solve the technical problem that the existing manual control and traditional ratio control methods cannot meet the requirement that the actual ratio of chlorine flow and hydrogen flow tracks the expected ratio, the invention provides a dynamic balance control method for the chlorine-hydrogen ratio in the hydrogen chloride preparation process, which can ensure that the real-time dynamic balance chlorine-hydrogen ratio in the production process always meets the expected chlorine-hydrogen ratio and effectively solves the problem of mismatch of the chlorine-hydrogen ratio in the hydrogen chloride production process. In order to achieve the technical purpose, the technical scheme provided by the invention is as follows:
a chlorine-hydrogen ratio dynamic balance control method in the process of preparing hydrogen chloride comprises the following steps:
step S01: calculating the loop characteristics of the chlorine and hydrogen flows, and determining the reference values of the changes of the chlorine and hydrogen flows through a loop flow model;
step S02: calculating an optimum load change rate based on the reference value of the flow rate in step S01 in combination with the output characteristics of the controller;
step S03: calculating a load increment target value according to the load change rate in step S02, and performing variable load control according to the calculated change rate;
step S04: calculating an optimal control output value, and adjusting the set values of the chlorine flow and the hydrogen flow according to the optimal control output value when the actual ratio of the chlorine flow to the hydrogen flow measurement value deviates from the target ratio;
step S05: and calculating the locking characteristic of the regulating valve, and performing valve locking control according to the locking characteristic of the regulating valve when the chlorine flow and hydrogen flow loops are disturbed.
Further, the step S01 includes the following steps:
s101: chlorine gas flow loop in acquisition control system
Figure BDA0003346258050000021
Historical data of process values and operating values as historical values of chlorine flow loop process values
Figure BDA0003346258050000022
And historical values of operational values
Figure BDA0003346258050000023
Hydrogen flow loop in acquisition control system
Figure BDA0003346258050000024
Historical data of process values and operating values as historical values of hydrogen flow loop process values
Figure BDA0003346258050000025
And historical values of operational values
Figure BDA0003346258050000026
By pairs
Figure BDA0003346258050000027
Performing zero equalization treatment to obtain
Figure BDA0003346258050000028
Wherein Y is processed data, X is historical data,
Figure BDA0003346258050000029
is the mean value of;
s102: the operation values are divided into interval one [0, 20), interval two [20, 40), interval three [40, 60), and interval four [60, 100]According to
Figure BDA00033462580500000210
Selecting corresponding Y according to the value, performing model identification on the chlorine flow loop and the hydrogen flow loop to obtain open-loop models of the chlorine flow loop corresponding to the four intervals
Figure BDA00033462580500000211
And
Figure BDA00033462580500000212
wherein
Figure BDA00033462580500000213
Representing a model of the chlorine flow interval one to four, K in the formula1-K4、T1-T4、θ1-θ4Respectively representing the gain, the time constant and the pure lag time of one to four models of the chlorine flow interval; hydrogen flow loop open loop model
Figure BDA00033462580500000214
Figure BDA00033462580500000215
And
Figure BDA00033462580500000216
wherein
Figure BDA00033462580500000217
One to four models representing the hydrogen flow interval, K in the formula1-K4、T1-T4、θ1-θ4Respectively representing the gain, the time constant and the pure lag time of one to four models of the hydrogen flow interval;
s103: utilizing the parameter K in the chlorine flow loop open-loop model in S1021、T1、θ1、K2、T2、θ2、K3、T3、θ3、K4、T4、θ4And parameter K in the open loop model of the hydrogen flow loop1、T1、θ1、K2、T2、θ2、K3、T3、θ3、K4、T4、θ4According to the current operating values of the chlorine flow and hydrogen flow loops
Figure BDA00033462580500000218
Corresponding intervals divided in S102 are respectively selected corresponding open-loop model parameters, and the reference values of the chlorine flow and the hydrogen flow change rate are calculated
Figure BDA0003346258050000031
And
Figure BDA0003346258050000032
wherein the reference value delta1In (1)
Figure BDA0003346258050000033
According to the current operation value
Figure BDA0003346258050000034
Selecting the region [ K ]1、K2、K3、K4]A value of (1) in the same way
Figure BDA0003346258050000035
Selecting [ T ]1、T2、T3、T4]One value of,
Figure BDA00033462580500000331
Selecting [ theta ]1、θ2、θ3、θ4]Of value, reference value delta2The hydrogen flow loop open loop model parameters are selected by the same method.
Further, the step S02 includes the steps of,
s201: obtaining the reference value delta of the change rate of the chlorine flow and the hydrogen flow1And delta2Taking the minimum value as the reference value delta of the load change rate3；
S202: calculating the actual rate of load change Δ δ3γ, where γ represents the control loop output characteristic, γ is the maximum output value per cycle of the controller x the output cycle of the balance controller.
Further, the step S03 includes the steps of,
s301: obtaining the final target load instruction value of chlorine flow
Figure BDA0003346258050000036
Chlorine to hydrogen ratio
Figure BDA0003346258050000037
Chlorine flow loop in control system
Figure BDA0003346258050000038
The set value is the current load set value of chlorine flow
Figure BDA0003346258050000039
And controlling hydrogen flow loop in system
Figure BDA00033462580500000310
The set value is hydrogen flowCurrent load set point
Figure BDA00033462580500000311
S302: calculating the chlorine flow increasing and decreasing load if
Figure BDA00033462580500000312
If the chlorine flow increases the load, then
Figure BDA00033462580500000313
If it is not
Figure BDA00033462580500000314
When the chlorine flow is reduced, then
Figure BDA00033462580500000315
If it is not
Figure BDA00033462580500000316
When the chlorine flow enters a steady state, then
Figure BDA00033462580500000317
S303: calculating the hydrogen flow rate increase/decrease load if
Figure BDA00033462580500000318
If the hydrogen flow rate is increased, then
Figure BDA00033462580500000319
If it is not
Figure BDA00033462580500000320
If the hydrogen flow is reduced, then
Figure BDA00033462580500000321
If it is not
Figure BDA00033462580500000322
Figure BDA00033462580500000323
If the hydrogen flow enters a steady state, then
Figure BDA00033462580500000324
Further, the step S04 includes the steps of,
s401: obtaining
Figure BDA00033462580500000325
Chlorine to hydrogen ratio
Figure BDA00033462580500000326
Hydrogen flow loop in control system
Figure BDA00033462580500000327
Current process value
Figure BDA00033462580500000328
Chlorine flow loop in control system
Figure BDA00033462580500000329
Current process value
Figure BDA00033462580500000330
Calculating conversion load for converting hydrogen gas flow into chlorine gas flow
Figure BDA0003346258050000041
Calculating the conversion load from chlorine flow to hydrogen flow
Figure BDA0003346258050000042
S402: comparison
Figure BDA0003346258050000043
And
Figure BDA0003346258050000044
taking the minimum value and assigning the value to a chlorine flow loop in a control system
Figure BDA0003346258050000045
Setting a value;
s403: comparison
Figure BDA0003346258050000046
And
Figure BDA0003346258050000047
taking the maximum value and assigning the maximum value to a hydrogen flow loop in a control system
Figure BDA0003346258050000048
And (5) setting the value.
Further, the step S05 includes the steps of,
s501: obtaining the chlorine-hydrogen ratio
Figure BDA0003346258050000049
And real-time ratio
Figure BDA00033462580500000410
S502: according to the real-time ratio RCVTo target ratio
Figure BDA00033462580500000411
Size, increase and close of chlorine regulating valve XLOCK+ MVAnd hydrogen regulating valve locking block XLOCK- MVControl, i.e. when
Figure BDA00033462580500000412
The chlorine valve can no longer be increased and the hydrogen flow valve can no longer be decreased, where krRepresenting a preset ratio margin.
Compared with the prior art, the invention has the following advantages and beneficial effects:
(1) the chlorine flow loop characteristic and the hydrogen flow loop characteristic are obtained through a model identification method, and the loop characteristic is reflected more accurately.
(2) And determining the actual variable quantity of the load in each period according to the model parameters and the output characteristics of the controller, so that the variable load in the production process is more stable.
(3) The chlorine flow regulating valve and the hydrogen flow regulating valve are controlled by locking, so that the safety of the production process is ensured.
(4) Dynamically balancing the ratio of chlorine flow to hydrogen flow, and constantly ensuring that the ratio of chlorine to hydrogen meets the process requirements.
Drawings
FIG. 1 is a schematic flow chart of a chlorine-hydrogen ratio dynamic balance control method of the present invention;
FIG. 2 is a schematic flow diagram of a loop characteristic of a calculated chlorine hydrogen flow;
FIG. 3 is a schematic flow chart of calculating an optimal load change rate;
FIG. 4 is a schematic view of a variable load control flow;
FIG. 5 is a schematic flow chart of calculating an optimal control output value;
FIG. 6 is a flow chart illustrating a calculation of a damper lock-out characteristic.
Detailed Description
The following describes the embodiments of the present invention in detail with reference to the drawings, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, not all embodiments.
The invention provides a chlorine-hydrogen ratio value dynamic balance control method in a hydrogen chloride preparation process, which more accurately reflects loop characteristics, determines the actual load variation of each period according to model parameters and controller output characteristics to make the load variation of the production process more stable, controls a chlorine flow regulating valve and a hydrogen flow regulating valve in a locking manner to ensure the safety of the production process, and constantly ensures that the ratio of chlorine-hydrogen meets the process requirements through dynamically balancing the ratio of the chlorine flow to the hydrogen flow.
A chlorine-hydrogen ratio dynamic balance control method for a hydrogen chloride preparation process is shown in figure 1 and comprises the following steps:
step S01: calculating the loop characteristics of chlorine and hydrogen flow
When the production load is increased or decreased within a certain range of the chlorine pressure and the hydrogen pressure, the loop characteristics of the chlorine and hydrogen flow determine the stable characteristics of the instantaneous chlorine flow and the instantaneous hydrogen flow entering the synthesis furnace. Model identification is carried out on each loop to obtain a loop flow model, and reference values of changes of chlorine and hydrogen flows are determined according to the model, as shown in fig. 2, the specific method comprises the following steps:
s101: chlorine gas flow loop in acquisition control system
Figure BDA0003346258050000051
Historical data of process values and operating values as historical values of chlorine flow loop process values
Figure BDA0003346258050000052
And historical values of operational values
Figure BDA0003346258050000053
Hydrogen flow loop in acquisition control system
Figure BDA0003346258050000054
Historical data of process values and operating values as historical values of hydrogen flow loop process values
Figure BDA0003346258050000055
And historical values of operational values
Figure BDA0003346258050000056
By pairs
Figure BDA0003346258050000057
Performing zero equalization treatment to obtain
Figure BDA0003346258050000058
Wherein Y is processed data, X is historical data,
Figure BDA0003346258050000059
is the mean value of.
S102: in general, the operation value is 0 to 100, and the range is divided into four sections, namely a section one [0, 20 ], a section two [20, 40 ] and a section three [2 ]40, 60) and interval four [60, 100%]. According to
Figure BDA00033462580500000510
Figure BDA00033462580500000511
Selecting corresponding Y according to the value, and respectively carrying out model identification on the chlorine flow loop and the hydrogen flow loop by adopting an autoregressive moving average model (ARMAX) identification method to obtain open-loop models of the chlorine flow loop corresponding to four intervals
Figure BDA00033462580500000512
Figure BDA00033462580500000513
Figure BDA0003346258050000061
And
Figure BDA0003346258050000062
wherein
Figure BDA0003346258050000063
Representing a chlorine flow interval model, K1、T1And theta1The chlorine flow interval is expressed, namely model gain, time constant and pure lag time, and other parameters have the same meaning. Similarly, an open loop model of the hydrogen flow loop is obtained
Figure BDA0003346258050000064
Figure BDA0003346258050000065
Figure BDA0003346258050000066
Figure BDA0003346258050000067
And
Figure BDA0003346258050000068
wherein
Figure BDA0003346258050000069
Model representing hydrogen flow rate interval, K1、T1And theta1The model gain, time constant and pure lag time of the hydrogen flow interval are shown, and other parameters have the same meanings and other meanings.
S103: utilizing the chlorine flow loop open loop model parameter K in the step (2)1、T1、θ1、K2、T2、θ2、K3、T3、θ3、K4、T4、θ4And hydrogen flow loop open loop model parameter K1、T1、θ1、K2、T2、θ2、K3、T3、θ3、K4、T4、θ4According to the current operating values of the chlorine flow and hydrogen flow loops
Figure BDA00033462580500000610
Corresponding intervals divided in the step (2) are respectively selected corresponding open-loop model parameters, and finally, the reference values of the chlorine flow and the hydrogen flow change rate are calculated
Figure BDA00033462580500000611
And
Figure BDA00033462580500000612
wherein the reference value delta1In (1)
Figure BDA00033462580500000613
According to the current operation value
Figure BDA00033462580500000614
Selecting the region [ K ]1、K2、K3、K4]A value of (1) in the same way
Figure BDA00033462580500000615
Selecting [ T ]1、T2、T3、T4]One value of,
Figure BDA00033462580500000616
Selecting [ theta ]1、θ2、θ3、θ4]Of value, reference value delta2The hydrogen flow loop open loop model parameters are selected by the same method.
Step S02: calculating an optimal load change rate
Based on the reference value of the flow rate in step S01, in combination with the output characteristics of the actual controller, an optimal load rate is determined, as shown in fig. 3, the method includes:
s201: obtaining the reference value delta of the change rate of the chlorine flow and the hydrogen flow1And delta2Is going on delta1And delta2Comparing and obtaining the minimum value as the reference value delta of the load change rate3。
S202: calculating the actual rate of load change Δ δ3γ, where γ represents the control loop output characteristic, γ is the maximum output value per cycle of the controller x the output cycle of the balance controller.
Step S03: variable load control
According to the actual rate of load change in the step two, after receiving a load change instruction, automatically calculating a load increment target value of each control period, and enabling the load to change according to the calculated change rate, as shown in fig. 4, the specific method comprises the following steps:
s301: acquiring the final target load command value of chlorine flow input by an operator
Figure BDA00033462580500000617
Chlorine-hydrogen ratio value input by operator
Figure BDA0003346258050000071
Chlorine flow loop in control system
Figure BDA0003346258050000072
The set value is the current load set value of chlorine flow
Figure BDA0003346258050000073
And controlling hydrogen flow loop in system
Figure BDA0003346258050000074
The set value is the current load set value of the hydrogen flow
Figure BDA0003346258050000075
H2Following Cl by ratio2。
S302: and calculating the load of chlorine flow increase and decrease. If it is not
Figure BDA0003346258050000076
If the chlorine flow increases the load, then
Figure BDA0003346258050000077
If it is not
Figure BDA0003346258050000078
When the chlorine flow is reduced, then
Figure BDA0003346258050000079
If it is not
Figure BDA00033462580500000710
When the chlorine flow enters a steady state, then
Figure BDA00033462580500000711
S303: and calculating the hydrogen flow increase and decrease load. If it is not
Figure BDA00033462580500000712
If the hydrogen flow rate is increased, then
Figure BDA00033462580500000713
If it is not
Figure BDA00033462580500000714
If the hydrogen flow is reduced, then
Figure BDA00033462580500000715
If it is not
Figure BDA00033462580500000716
Figure BDA00033462580500000717
If the hydrogen flow enters a steady state, then
Figure BDA00033462580500000718
Step S04: and calculating an optimal control output value.
The controller outputs the result according to the third step, when the actual ratio of the measured values of the chlorine flow and the hydrogen flow deviates from the target ratio, the system will automatically adjust the set values of the chlorine flow and the hydrogen flow, as shown in fig. 5, the specific method includes:
s401: obtained by real-time calculation of S302 in "step S03
Figure BDA00033462580500000719
Chlorine-hydrogen ratio value input by operator
Figure BDA00033462580500000720
Hydrogen flow loop in control system
Figure BDA00033462580500000721
Current process value
Figure BDA00033462580500000722
Chlorine flow loop in control system
Figure BDA00033462580500000723
Current process value
Figure BDA00033462580500000724
Calculating conversion load for converting hydrogen gas flow into chlorine gas flow
Figure BDA00033462580500000725
Calculating the conversion load from chlorine flow to hydrogen flow
Figure BDA00033462580500000726
Figure BDA00033462580500000727
S402: comparison
Figure BDA00033462580500000728
And
Figure BDA00033462580500000729
taking the minimum value and assigning the value to a chlorine flow loop in a control system
Figure BDA00033462580500000730
And (5) setting the value.
S403: comparison
Figure BDA00033462580500000731
And
Figure BDA00033462580500000732
taking the maximum value and assigning the maximum value to a hydrogen flow loop in a control system
Figure BDA0003346258050000081
And (5) setting the value.
Step S05: calculating the locking characteristic of the regulating valve;
when the chlorine flow and hydrogen flow loops are disturbed, the flow can fluctuate greatly, and the excessive chlorine entering the synthesis furnace is prevented by adopting valve locking control, as shown in fig. 6, the specific method comprises the following steps:
s501: acquiring the chlorine-hydrogen ratio value input by an operator
Figure BDA0003346258050000082
And real-time ratio
Figure BDA0003346258050000083
S502: according to the real-time ratio RCVTo target ratio
Figure BDA0003346258050000084
Size, increase and close of chlorine regulating valve XLOCK+ MVAnd hydrogen regulating valve locking block XLOCK- MVControl, i.e. when
Figure BDA0003346258050000085
The chlorine valve can no longer be increased and the hydrogen flow valve can no longer be decreased, where krRepresenting a preset ratio margin.
Example 1
A chlorine-hydrogen ratio dynamic balance control method for preparing hydrogen chloride is carried out by the following steps,
step S01: calculating the loop characteristics of chlorine and hydrogen flow
Further, in step S01, the content specifically required to be calculated includes:
step S101: historical data of chlorine and hydrogen flow loop process values and operating values are obtained from a control system (such as DCS or PLC), and zero-mean processing is carried out on the data.
For example, historical chlorine flow loop process value data is obtained
Figure BDA0003346258050000086
And operation value history data
Figure BDA0003346258050000087
Figure BDA0003346258050000088
Figure BDA0003346258050000089
Figure BDA0003346258050000091
Computing
Figure BDA0003346258050000092
Mean value
Figure BDA0003346258050000093
Figure BDA0003346258050000094
Mean value
Figure BDA0003346258050000095
By passing
Figure BDA0003346258050000096
Calculating the data after the chlorine flow loop treatment:
Figure BDA0003346258050000097
Figure BDA0003346258050000098
the data of the hydrogen flow loop can be processed in the same way and will not be described in detail here.
Step S102: obtaining zero-averaged data
Figure BDA0003346258050000099
According to the operating value
Figure BDA00033462580500000910
And respectively carrying out model identification on the chlorine flow loop and the hydrogen flow loop in different intervals.
Within the range of 0 to 100, to
Figure BDA00033462580500000911
Carrying out interval division, wherein an interval I:
Figure BDA00033462580500000912
the interval two:
Figure BDA0003346258050000101
interval three:
Figure BDA0003346258050000102
interval four:
Figure BDA0003346258050000103
using model identification algorithms
Figure BDA0003346258050000104
The algorithm can be implemented on python, where na 1, nb 1, nc 1 and nk 1 are parameters set for identification, and a mathematical model of the interval one can be obtained
Figure BDA0003346258050000105
Correspondingly, by adopting the same method, a secondary mathematical model of the chlorine flow loop interval can be obtained by utilizing historical data
Figure BDA0003346258050000106
Three interval mathematical models:
Figure BDA0003346258050000107
interval four mathematical model
Figure BDA0003346258050000108
Accordingly, historical data is utilized for hydrogen flow circuit operating values
Figure BDA0003346258050000109
Four intervals are divided, interval one:
Figure BDA00033462580500001010
the interval two:
Figure BDA00033462580500001011
interval three:
Figure BDA00033462580500001012
interval four:
Figure BDA00033462580500001013
using a model recognition tool G ═ ARMAX (Y, [ na nb nc nk ]]) The model identification method is characterized in that Y represents data after zero equalization, na is 1, nb is 1, nc is 1 and nk is 1 which are parameters set during identification, and a mathematical model of a hydrogen flow loop interval is obtained
Figure BDA00033462580500001014
Interval two-mathematic model:
Figure BDA00033462580500001015
interval three-mathematical model
Figure BDA00033462580500001016
Interval four mathematical model
Figure BDA00033462580500001017
Step S103: according to the current operating value of the chlorine flow circuit
Figure BDA00033462580500001018
Current operating value of hydrogen flow loop
Figure BDA00033462580500001019
In the section, select "in step S102"
Figure BDA00033462580500001020
Figure BDA00033462580500001021
The open-loop model in (1) calculates the chlorine flow rate,A reference value of the hydrogen gas flow rate.
Such as the current chlorine flow loop operating value
Figure BDA00033462580500001022
Selecting
Figure BDA00033462580500001023
Hydrogen flow loop operating value
Figure BDA00033462580500001024
Selecting
Figure BDA00033462580500001025
According to
Figure BDA00033462580500001026
And
Figure BDA00033462580500001027
calculating a reference value delta of the chlorine flow rate change rate116.73/(61+9) ≈ 0.239 and a reference value δ for calculating a hydrogen flow rate change rate2＝16.11/(72+14)≈0.187。
Step S02: calculating an optimal variable load rate
Further, the content specifically required to be calculated in step S02 includes:
step S201: obtaining the reference value delta of the change rate of the chlorine flow and the hydrogen flow1About 0.239 and delta2≈0.187，δ1And delta2Comparing to obtain the minimum value as the reference value delta of the load change rate3＝δ20.187. Step S202: reference value delta of load change rate3When the output characteristic γ of the control loop is 0.187, the maximum output value of the controller per cycle × the output cycle of the balance controller is 0.8 × 1, 0.8, and finally the actual rate of the load change is obtained: delta-delta3×γ＝0.187×0.8≈0.1496。
Step S03: variable load control
Further, the content specifically required to be calculated in step S03 includes:
step S301: obtaining the final target load instruction value input by the operator
Figure BDA0003346258050000111
Operator input of hydrogen to hydrogen ratio
Figure BDA0003346258050000112
Chlorine flow loop in control system
Figure BDA0003346258050000113
The set value is the current load set value of chlorine flow
Figure BDA0003346258050000114
And controlling hydrogen flow loop in system
Figure BDA0003346258050000115
The set value is the current load set value of the hydrogen flow
Figure BDA0003346258050000116
Step S302: due to the fact that
Figure BDA0003346258050000117
Increase the load, therefore
Figure BDA0003346258050000118
Figure BDA0003346258050000119
Up to the point where
Figure BDA00033462580500001110
Until now.
Step S303: due to the fact that
Figure BDA00033462580500001111
Increase the load, therefore
Figure BDA00033462580500001112
Figure BDA00033462580500001113
Up to the point where
Figure BDA00033462580500001114
Figure BDA00033462580500001115
Step S04: calculating optimal control output
Step S401: obtaining real-time calculations in "step S302
Figure BDA00033462580500001116
Operator input of hydrogen to hydrogen ratio
Figure BDA00033462580500001117
Hydrogen flow loop in control system
Figure BDA00033462580500001118
Process value
Figure BDA00033462580500001119
Figure BDA00033462580500001120
Chlorine flow loop in control system
Figure BDA00033462580500001121
Current process value
Figure BDA00033462580500001122
The conversion load of hydrogen gas flow into chlorine gas flow
Figure BDA00033462580500001123
Load conversion of chlorine flow to hydrogen flow
Figure BDA00033462580500001124
Step S402: comparative chlorine gas
Figure BDA00033462580500001125
And converting the load
Figure BDA00033462580500001126
Taking the minimum value 40.1496, and assigning a value to a chlorine flow loop in the control system
Figure BDA00033462580500001127
And (5) setting the value.
Step S403: comparison of Hydrogen gas
Figure BDA0003346258050000121
And converting the load
Figure BDA0003346258050000122
Taking a maximum value 44.1646, and assigning a value to a hydrogen flow loop in the control system
Figure BDA0003346258050000123
And (5) setting the value.
Step S05: calculating loop regulator valve lockout characteristics
Further, the content specifically required to be calculated in step S05 includes:
step S501: obtaining an operator entered hydrogen to hydrogen ratio value
Figure BDA0003346258050000124
The ratio margin kr input by the operator is 0.95, and a hydrogen flow loop in the control system
Figure BDA0003346258050000125
Process value
Figure BDA0003346258050000126
Chlorine flow loop in control system
Figure BDA0003346258050000127
Current process value
Figure BDA0003346258050000128
Hydrogen flow loop in control system
Figure BDA0003346258050000129
Current operating value
Figure BDA00033462580500001210
Chlorine flow loop in control system
Figure BDA00033462580500001211
Current operating value
Figure BDA00033462580500001212
Figure BDA00033462580500001213
Hydrogen flow loop in control system
Figure BDA00033462580500001214
Lower limit value of current operation value
Figure BDA00033462580500001215
Chlorine flow loop in control system
Figure BDA00033462580500001216
Upper limit value of current operation value
Figure BDA00033462580500001217
Calculating a real-time ratio
Figure BDA00033462580500001218
Figure BDA00033462580500001219
Step S502: due to the fact that
Figure BDA00033462580500001220
Performing chlorine regulating valve locking
Figure BDA00033462580500001221
And hydrogen regulating valve reducing lock
Figure BDA00033462580500001222
When in use
Figure BDA00033462580500001223
When the circuit is in use, the locking control is reset,
Figure BDA00033462580500001224
although embodiments of the present invention have been shown and described above, it should be understood that the above embodiments are illustrative and not restrictive, and that those skilled in the art may make changes, modifications, substitutions and alterations to the above embodiments without departing from the scope of the present invention.

Claims (6)

1. A chlorine-hydrogen ratio dynamic balance control method in a hydrogen chloride preparation process is characterized by comprising the following steps:
step S01: calculating the loop characteristics of the chlorine and hydrogen flows, and determining the reference values of the changes of the chlorine and hydrogen flows through a loop flow model;
step S02: calculating an optimum load change rate based on the reference value of the flow rate in step S01 in combination with the output characteristics of the controller;
step S03: calculating a load increment target value according to the load change rate in step S02, and performing variable load control according to the calculated change rate;
step S04: calculating an optimal control output value, and adjusting the set values of the chlorine flow and the hydrogen flow according to the optimal control output value when the actual ratio of the chlorine flow to the hydrogen flow measurement value deviates from the target ratio;
step S05: and calculating the locking characteristic of the regulating valve, and performing valve locking control according to the locking characteristic of the regulating valve when the chlorine flow and hydrogen flow loops are disturbed.
2. The method for controlling dynamic balance of chlorine-hydrogen ratio in process of preparing hydrogen chloride according to claim 1, wherein said step S01 comprises the steps of:
s101: chlorine gas flow loop in acquisition control system
Figure FDA0003346258040000011
Historical data of process values and operating values as historical values of chlorine flow loop process values
Figure FDA0003346258040000012
And historical values of operational values
Figure FDA0003346258040000013
Hydrogen flow loop in acquisition control system
Figure FDA0003346258040000014
Historical data of process values and operating values as historical values of hydrogen flow loop process values
Figure FDA0003346258040000015
And historical values of operational values
Figure FDA0003346258040000016
By pairs
Figure FDA0003346258040000017
Performing zero equalization treatment to obtain
Figure FDA0003346258040000018
Wherein Y is processed data, X is historical data,
Figure FDA0003346258040000019
is the mean value of;
s102: the operation values are divided into interval one [0, 20), interval two [20, 40), interval three [40, 60), and interval four [60, 100]According to
Figure FDA00033462580400000110
Selecting corresponding Y according to the value, performing model identification on the chlorine flow loop and the hydrogen flow loop to obtain open-loop models of the chlorine flow loop corresponding to the four intervals
Figure FDA00033462580400000111
And
Figure FDA00033462580400000112
wherein
Figure FDA00033462580400000113
Representing a model of the chlorine flow interval one to four, K in the formula1-K4、T1-T4、θ1-θ4Respectively representing the gain, the time constant and the pure lag time of one to four models of the chlorine flow interval; hydrogen flow loop open loop model
Figure FDA0003346258040000021
Figure FDA0003346258040000022
Wherein
Figure FDA0003346258040000023
One to four models representing the hydrogen flow interval, K in the formula1-K4、T1-T4、θ1-θ4Respectively representing the gain, the time constant and the pure lag time of one to four models of the hydrogen flow interval;
s103: utilizing the parameter K in the chlorine flow loop open-loop model in S1021、T1、θ1、K2、T2、θ2、K3、T3、θ3、K4、T4、θ4And parameter K in the open loop model of the hydrogen flow loop1、T1、θ1、K2、T2、θ2、K3、T3、θ3、K4、T4、θ4According to the current operating values of the chlorine flow and hydrogen flow loops
Figure FDA0003346258040000024
Corresponding intervals divided in S102 are respectively selected corresponding open-loop model parameters, and the reference values of the chlorine flow and the hydrogen flow change rate are calculated
Figure FDA0003346258040000025
And
Figure FDA0003346258040000026
wherein the reference value delta1In (1)
Figure FDA0003346258040000027
According to the current operation value
Figure FDA0003346258040000028
Selecting the region [ K ]1、K2、K3、K4]A value of (1) in the same way
Figure FDA0003346258040000029
Selecting [ T ]1、T2、T3、T4]One value of,
Figure FDA00033462580400000210
Selecting [ theta ]1、θ2、θ3、θ4]Of value, reference value delta2The hydrogen flow loop open loop model parameters are selected by the same method.
3. The method for controlling dynamic balance of hydrogen chloride ratio in process of preparing hydrogen chloride according to claim 1, wherein said step S02 includes the steps of,
s201: obtaining references of chlorine flow and hydrogen flow change rateThe value delta1And delta2Taking the minimum value as the reference value delta of the load change rate3；
S202: calculating the actual rate of load change Δ δ3γ, where γ represents the control loop output characteristic, γ is the maximum output value per cycle of the controller x the output cycle of the balance controller.
4. The method for controlling dynamic balance of hydrogen chloride ratio in process of preparing hydrogen chloride according to claim 1, wherein said step S03 includes the steps of,
s301: obtaining the final target load instruction value of chlorine flow
Figure FDA00033462580400000211
Chlorine to hydrogen ratio
Figure FDA00033462580400000212
Chlorine flow loop in control system
Figure FDA00033462580400000213
The set value is the current load set value of chlorine flow
Figure FDA00033462580400000214
And controlling hydrogen flow loop in system
Figure FDA00033462580400000215
The set value is the current load set value of the hydrogen flow
Figure FDA00033462580400000216
S302: calculating the chlorine flow increasing and decreasing load if
Figure FDA00033462580400000217
If the chlorine flow increases the load, then
Figure FDA0003346258040000031
If it is not
Figure FDA0003346258040000032
When the chlorine flow is reduced, then
Figure FDA0003346258040000033
If it is not
Figure FDA0003346258040000034
When the chlorine flow enters a steady state, then
Figure FDA0003346258040000035
S303: calculating the hydrogen flow rate increase/decrease load if
Figure FDA0003346258040000036
If the hydrogen flow rate is increased, then
Figure FDA0003346258040000037
If it is not
Figure FDA0003346258040000038
If the hydrogen flow is reduced, then
Figure FDA0003346258040000039
If it is not
Figure FDA00033462580400000310
Figure FDA00033462580400000311
If the hydrogen flow enters a steady state, then
Figure FDA00033462580400000312
5. The method for controlling dynamic balance of hydrogen chloride ratio in process of preparing hydrogen chloride according to claim 1, wherein said step S04 includes the steps of,
s401: obtaining
Figure FDA00033462580400000313
Chlorine to hydrogen ratio
Figure FDA00033462580400000314
Hydrogen flow loop in control system
Figure FDA00033462580400000315
Current process value
Figure FDA00033462580400000316
Chlorine flow loop in control system
Figure FDA00033462580400000317
Current process value
Figure FDA00033462580400000318
Calculating conversion load for converting hydrogen gas flow into chlorine gas flow
Figure FDA00033462580400000319
Calculating the conversion load from chlorine flow to hydrogen flow
Figure FDA00033462580400000320
S402: comparison
Figure FDA00033462580400000321
And
Figure FDA00033462580400000322
taking the minimum value and assigning the value to a chlorine flow loop in a control system
Figure FDA00033462580400000323
Setting a value;
s403: comparison
Figure FDA00033462580400000324
And
Figure FDA00033462580400000325
taking the maximum value and assigning the maximum value to a hydrogen flow loop in a control system
Figure FDA00033462580400000326
And (5) setting the value.
6. The method for controlling dynamic balance of hydrogen chloride ratio in process of preparing hydrogen chloride according to claim 1, wherein said step S05 includes the steps of,
s501: obtaining the chlorine-hydrogen ratio
Figure FDA00033462580400000327
And real-time ratio
Figure FDA00033462580400000328
S502: according to the real-time ratio RCVTo target ratio
Figure FDA00033462580400000329
Size, increase and close of chlorine regulating valve XLOCK+ MVAnd hydrogen regulating valve locking block XLOCK- MVControl, i.e. when
Figure FDA00033462580400000330
The chlorine valve can no longer be increased and the hydrogen flow valve can no longer be decreased, where krRepresenting a preset ratio margin.
CN202111324006.0A 2021-11-10 2021-11-10 Dynamic balance control method for chlorine-hydrogen ratio in process of preparing hydrogen chloride Active CN114185369B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111324006.0A CN114185369B (en) 2021-11-10 2021-11-10 Dynamic balance control method for chlorine-hydrogen ratio in process of preparing hydrogen chloride

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111324006.0A CN114185369B (en) 2021-11-10 2021-11-10 Dynamic balance control method for chlorine-hydrogen ratio in process of preparing hydrogen chloride

Publications (2)

Publication Number Publication Date
CN114185369A true CN114185369A (en) 2022-03-15
CN114185369B CN114185369B (en) 2023-10-20

Family

ID=80540844

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111324006.0A Active CN114185369B (en) 2021-11-10 2021-11-10 Dynamic balance control method for chlorine-hydrogen ratio in process of preparing hydrogen chloride

Country Status (1)

Country Link
CN (1) CN114185369B (en)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110081561A1 (en) * 2009-05-29 2011-04-07 Majid Keshavarz Methods of producing hydrochloric acid from hydrogen gas and chlorine gas
CN105259940A (en) * 2015-11-24 2016-01-20 青海盐湖工业股份有限公司 Control system of hydrogen chloride synthesis purity
CN108715437A (en) * 2018-07-20 2018-10-30 唐山三友氯碱有限责任公司 The DCS control methods of load adjustment in hydrochloric acid production system and hydrogen chloride production
CN110589769A (en) * 2019-09-17 2019-12-20 德州实华化工有限公司 Automatic chlorine-hydrogen ratio control method and system for hydrogen chloride synthesis furnace and synthesis furnace
CN111503520A (en) * 2020-04-17 2020-08-07 浙江中智达科技有限公司 Air intake load control method, device, equipment and readable storage medium
CN112578745A (en) * 2020-09-28 2021-03-30 山东鲁泰化学有限公司 Intelligent control method for hydrogen chloride synthesis reaction process
CN113124318A (en) * 2021-04-20 2021-07-16 南通星球石墨股份有限公司 Automatic hydrogen and chlorine proportioning system for hydrochloric acid synthesis furnace and control method

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110081561A1 (en) * 2009-05-29 2011-04-07 Majid Keshavarz Methods of producing hydrochloric acid from hydrogen gas and chlorine gas
CN105259940A (en) * 2015-11-24 2016-01-20 青海盐湖工业股份有限公司 Control system of hydrogen chloride synthesis purity
CN108715437A (en) * 2018-07-20 2018-10-30 唐山三友氯碱有限责任公司 The DCS control methods of load adjustment in hydrochloric acid production system and hydrogen chloride production
CN110589769A (en) * 2019-09-17 2019-12-20 德州实华化工有限公司 Automatic chlorine-hydrogen ratio control method and system for hydrogen chloride synthesis furnace and synthesis furnace
CN111503520A (en) * 2020-04-17 2020-08-07 浙江中智达科技有限公司 Air intake load control method, device, equipment and readable storage medium
CN112578745A (en) * 2020-09-28 2021-03-30 山东鲁泰化学有限公司 Intelligent control method for hydrogen chloride synthesis reaction process
CN113124318A (en) * 2021-04-20 2021-07-16 南通星球石墨股份有限公司 Automatic hydrogen and chlorine proportioning system for hydrochloric acid synthesis furnace and control method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
胡丹宁;: "比值自动控制调节系统在氯化氢合成中的应用", 中国氯碱, no. 03, pages 37 - 41 *

Also Published As

Publication number Publication date
CN114185369B (en) 2023-10-20

Similar Documents

Publication Publication Date Title
JPH01109402A (en) Apparatus and method using adaptive gain scheduling algorism
CN104753439B (en) A kind of PID intelligent speed-regulating methods of motor
CN112965548B (en) Automatic control method for temperature of reaction kettle, and upper and lower limit controllers and system for opening degree of valve
CN108227500A (en) A kind of control method for coordinating and system of the quick peak regulation of fired power generating unit
CN114185369B (en) Dynamic balance control method for chlorine-hydrogen ratio in process of preparing hydrogen chloride
CN106054616B (en) The titanium strip coil continuous acid-washing looper height control method of fuzzy logic PID controller parameter
Tadeo et al. Control of neutralization processes by robust loop shaping
CN102873106B (en) Quick and precise elongation control method for temper mill
CN106468879A (en) A kind of liquid level flow nonlinear area control method
CN112578745B (en) Intelligent control method for hydrogen chloride synthesis reaction process
CN113325696A (en) Hybrid control method for combining single neuron PID and model prediction applied to crosslinked cable production equipment
CN116483003A (en) Control method, device and equipment for coal water slurry manufacturing process
CN103344115B (en) Novel sintering machine speed control method
CN105786055A (en) Control system and control method for ammonia-to-air ratio automatic set value in nitric acid production oxidation furnace
JP2001027903A (en) Automatic control method
CN108089442A (en) A kind of PI controller parameter automatic setting methods based on Predictive function control and fuzzy control
CN114428521A (en) Chlorine flow control method based on vortex shedding flowmeter
CN113300386A (en) Frequency controller design method and system based on alternating direction multiplier method
CN106698442A (en) Polysilicon coarse distillation control method
Tsamatsoulis et al. PID parameterization of cement kiln precalciner based on simplified modeling
CN113587120B (en) Control method of plasma ash melting furnace
CN117148707A (en) Method, device and equipment for controlling variable decoupling of mixed gas and storage medium
Higham A Different Approach for Self-Tuning in Process Controllers—The Case for Introducing an Expert System
CN114094637A (en) Active power adjusting method and device for speed regulator of hydroelectric generating set
CN114183692A (en) Automatic control method for four butterfly valves in gas mixing system

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