CN106896721A - A kind of binary distillation column centerized fusion method and system - Google Patents
A kind of binary distillation column centerized fusion method and system Download PDFInfo
- Publication number
- CN106896721A CN106896721A CN201710163840.3A CN201710163840A CN106896721A CN 106896721 A CN106896721 A CN 106896721A CN 201710163840 A CN201710163840 A CN 201710163840A CN 106896721 A CN106896721 A CN 106896721A
- Authority
- CN
- China
- Prior art keywords
- parameter
- inversion model
- matrix
- frequency characteristic
- omega
- 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
Links
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
Landscapes
- 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)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
The present invention discloses a kind of binary distillation column centerized fusion method and system.By recognizing the channel transfer function between each controlled variable and control variables, the model of binary distillation column is set up.Binary distillation column model to obtaining is asked for and analyzes the frequency domain characteristic of its inversion model to determine the filter parameter in centralized PID controller, and the regulation parameter of PID controller is determined further according to desired dynamic property and robustness.Control signal is adjusted using the centralized PID controller for calculating so that two-component fractionation column overhead reaches requirement with bottom of towe light component with the purity of heavy constituent.This method thinking is simple, it is easy to accomplish, effective control can be implemented to the binary distillation column that there is coupling.
Description
Technical field
The invention belongs to the automation field of petrochemical industry rectifying column, particularly binary distillation column centerized fusion method
And system.
Background technology
Rectifying is the mass transport process being most widely used in petrochemical industry, and purification and the recovery of about 90% product are to pass through
Rectifying is completed.The mechanism that rectifying is separate mainly uses the volatility difference of different component in mixture to be separated.Mixture
By after rectifying column, the product produced from tower top or bottom of towe meets certain purity requirement.Binary distillation column is mainly used to separate
Containing two kinds of mixtures of component.Binary distillation column is a multivariable process, and internal mechanism is complicated, and dynamic response is slow and controls
Coupling between variable and controlled variable is more serious, therefore requirement higher is proposed to control program.In rectifying column
Control program in, relatively simple is that the distillation process of multivariable is split into multiple single loops by the method for variable pairing
It is controlled, this method is more effective when single-ended component is controlled, but it is not carried out to the coupling between variable
Effective compensation, and the criterion of pairing is often according to static Relative increasing rate (RGA), therefore when carrying out two ends to tower
It is necessary to carry out comprehensive analysis to process during component control.Uneoupled control is a good method coupled between eliminating variable, but existing
Uneoupled control in, though static decoupling easily realize, equally face ignore process dynamics information problem, and ideal decoupling or
Process is difficult to ask for again after its decouplers of dynamic decoupling method such as simple decoupling or decoupling, and this all brings difficulty to design controller
Degree.Now about 90% controller is still PID controller in the industry, how to design effective PID controller, it is controlled well
Binary distillation column processed is a good problem to study.
The content of the invention
In view of this, the main object of the present invention is to provide a kind of binary distillation column centerized fusion method.Including following
Step:
S1, each input of measurement binary distillation column and the numerical value of each output, and will measure the value pick-up for obtaining is electricity
Signal;
S2, the electric signal that S1 steps are obtained is stored, by after storage time t, electric signal being sent into system and being distinguished
Know module;
S3, the transmission function that each control passage is recognized using linear least squares method algorithm, set up the biography of binary distillation column
Delivery function matrix model;
S4, the transfer function model obtained according to S3 steps, are calculated the frequency characteristic of normalization inversion model, Ran Houtong
Cross analysis inversion model frequency characteristic design wave filter.
The wave filter that S5, the transfer function model obtained according to S3 steps and S4 steps are obtained, according to given robustness
And dynamic performance requirements, PID controller parameter is calculated, and determine final centralized PID controller;
S6, adjust input using the centralized PID controller that is obtained by S5 steps.
The storage time t is from the step signal moment for recognizing each passage is added, to be reached to output new steady
Time period untill state value.
The detailed process of S3 is:The output Y of amplitude h, sampling period Ts and acquisition target according to the step signal
K (), is calculated the parameter in the transmission function of each passage according to the following formula
θ=(PT·P)-1PTZ
Wherein P and Z are respectively the coefficient matrixes being made up of h, Ts, Y (k), and θ is the transmission function for identifying each branch road
Parameter.
S4's concretely comprises the following steps:
S401, the steady-state gain inverse matrix K for asking for process.
S402, according to process transfer function matrix obtain respectively and export 1 corresponding all passages in minimum time lag τ1
And with output 2 corresponding all passages in minimum time lag τ2.。
S403, will with output 1 corresponding all passages time lag part reduce τ1, with 2 corresponding all passages of output
Time lag part reduces τ2Afterwards, the model G for designing wave filter is obtained0。
S404, selected peak frequency ωmax, in frequency range [0, ωmax] on draw normalized inversion model each
ElementFrequency characteristic (Nyquist figures).Peak frequency ωmaxMeet all elements of normalization inversion model
Frequency characteristic all in the RHP of complex plane.The frequency characteristic of normalized inversion model element is obtained as the following formula
Wherein, j is imaginary unit, and ω is frequency, and K is the steady-state gain inverse matrix of process,It is inversion model frequency
Feature matrix,It is normalization inversion model frequency characteristic matrix, subscript i and k represent the i-th row and the k row of matrix respectively
Element.
S405, the frequency characteristic to the element of normalization inversion model are analyzed to determine controller parameter by column
γkValue, subscript k represents parameter γkThe kth column element of correspondence normalization inversion model, parameter γkSelection standard be selection one
Individual smaller value causes the frequency characteristic through the kth row all elements of the normalization inversion model after overcompensation all in stable region
Domain, the kth column element for normalizing inversion model is compensated according to the following formula
Wherein, j is imaginary unit, and ω is frequency, γkIt is corresponding controller parameter,It is normalization inversion model
Frequency characteristic matrix,It is the normalization inversion model frequency characteristic matrix after compensation, subscript i and k represent matrix respectively
The element of the i-th row and kth row,
S406, using the normalization inversion model after the compensation for obtaining, taken on the frequency characteristic of each of which element N number of
Point, the corresponding frequency of these points is respectively ω1, ω2... ωN, wherein, ω1=0.The method of compound curve fitting is recycled to obtain
To the transmission function of corresponding wave filter
Wherein, Qik(s) be by compensating after normalization inversion model the i-th row kth column element fitting obtain wave filter,
AikAnd BikIt is the corresponding parameter value of wave filter, KikIt is the i-th row kth column element of the steady-state gain inverse matrix of process, s is general to draw
Laplacian operater, parameter is asked for according to following compound curve approximating method:
Wherein, Qik(jωm) it is the i-th row kth column element, take ω in frequency on frequency characteristicmWhen corresponding plural number, |
Qik(jωm) | it is to ask for the plural modulus value, Re (Qik(jωm)) it is to ask for the real, Im (Qik(jωm)) it is to ask for
The plural imaginary part, N is the number of the Frequency point chosen, and a, b, c, d are intermediate variable.
The detailed process of S5 is
The value of S501, given robustness index Ms, controller regulation parameter λ1And λ2Ask for according to the following formula
S502, the sub- PID controller parameter of calculating
Wherein, KpIt is proportional gain, KIIt is storage gain, KdIt is the differential gain, TfIt is filter time constant, these ginsengs
Several specific acquiring methods is as follows:
The centralized PID controller that can obtain final band filter is:
Wherein, C is total centralized PID controller, C1And C2It is the above-mentioned sub- PID controller asked for, QikTo ask before
The wave filter for taking.
Second aspect, the present invention provides a kind of binary distillation column Centralized Control System, including the measurement being sequentially connected becomes
Device and data storage and output unit are sent, being also linked in sequence after data storage and output unit has System Discrimination unit, parameter
Analytic unit and controller unit, wherein
Measuring transducer, for gathering input signal and output signal;
Data storage and output unit, the data that measuring transducer transmission comes are stored, and export to system
Identification unit;
System Discrimination unit, the transmission function of each passage is identified using identification algorithm, constitutes transfer function matrix;
Parameter analysis unit, the transfer function matrix to process is analyzed, and designs centralized wave filter, and using fixed
Robustness index calculates PID controller parameter;
Controller unit, centralized PID controller is constituted using the filter parameter and PID controller parameter for obtaining, from
And output control signal.
According to the controller that this method is designed, coupling can be effectively reduced for the influence that system is exported, so that
Effective lifting tower top and the product quality of bottom of towe extraction.And controller architecture is simple, it is easy to accomplish, controller parameter can be straight
Obtain and take, it is to avoid the process rule of thumb adjusted.
Brief description of the drawings
Fig. 1 is the flow chart of embodiment of the present invention binary distillation column centerized fusion method.
Fig. 2 is the structure chart of embodiment of the present invention binary distillation column Centralized Control System.
Fig. 3 is frequency characteristic of each element of embodiment of the present invention normalization inversion model in regulation frequency range.
Fig. 4 is the judgement of stability schematic diagram of each element of embodiment of the present invention normalization inversion model.
Fig. 5 is frequency characteristic of each element of the normalization inversion model after embodiment of the present invention compensation in regulation frequency range
Curve.
Specific embodiment
The present invention is described in detail with specific embodiment below in conjunction with the accompanying drawings.
So that component is for the binary distillation column of first alcohol and water as an example, the input (control variables) that the process is chosen is returned for tower top
The steam flow of flow and bottom of towe, output (controlled variable) selection is the molar fraction of tower top methyl alcohol component and rubbing for bottom of towe methyl alcohol
That fraction.There is stronger coupling between the input and output of process, control targe is the methyl alcohol in the product of tower top and bottom of towe
Molar fraction be intended to reach the requirement of setting.Binary distillation column centerized fusion embodiment of the method is as shown in figure 1, including following
Step:
Step 1, the input end measuring overhead reflux value from process, tower bottom steam flow value, from the measurement of the output of process end
The Methanol Molar fraction of tower top and bottom of towe, is electric signal by flow value and molar fraction pick-up;
Electric signal in step 2, storing step 1, until after the quantity required for meeting each passage of identification, electric signal being passed
Deliver to System Discrimination module;
Step 3, the transmission function that each passage is recognized according to the linear least squares method algorithm of setting, and the transmission for picking out
Function constitutes transfer function matrix;
Recognized, the transfer function matrix for obtaining is:
Wherein, y1And y2The molar fraction of the methyl alcohol respectively in tower top and tower bottom product, R and S are respectively overhead reflux amount
Steam flow is heated with bottom of towe.
Step 4, the process model obtained using identification, the centralized wave filter of analysis design.
(1) Steady-state process gain inverse matrix K can be obtained first
(2) the minimum time lag related to 1 (molar fraction of methyl alcohol in the product of overhead extraction) of output is 1, with output 1
(molar fraction of methyl alcohol in the product of bottom of towe extraction) related minimum time lag is 3, i.e. τ1=1, τ2=3.
(3) the time lag part with 1 corresponding all passages of output is reduced into τ1, with output 2 corresponding all passages when
Stagnant part reduces τ2, obtain the model G for designing wave filter0。
(4) peak frequency ω is chosenmax=10-2, then [0, ωmax] on, according to G0Draw the every of normalization inversion model
The frequency characteristic of individual element, as shown in Figure 3.
(5) stability of analytical element is so that it is determined that controller parameter γk, the stability of element judged according to Fig. 4.
In Fig. 4, be divided into for the RHP of complex plane by real axis and the dotted line (straight line that the permanent plural number for 1 of real part is constituted) perpendicular to real axis
4 regions, if the frequency characteristic of normalization inversion model element is in region II or region III, element be it is stable, such as
Fruit is in region I or region IV, and element is unstable, therefore normalization inversion model is compensated so that the normalization after compensation
The frequency characteristic of all elements of inversion model is all in region II or region III.So wave filter is by using the normalizing after compensation
Change inversion model to ask for.Unnecessary dynamic will certainly be introduced after due to introducing compensator, therefore is to use up when controller is designed
Possible this partial dynamic of counteracting, so compensator parameter γkOccur in the controller, as controller parameter.Parameter γk
Selection standard be selection one smaller value so that through after overcompensation normalization inversion model kth row all elements frequency
Characteristic curve all in stability region, the element of the former normalization permutation of inversion model one will be compensated due to same compensator, for example, containing
There is parameter γ1Compensator will act on two elements of first row, contain parameter γ2Compensator will act on secondary series
Two elements, therefore γ1Will be selected to be can make two smaller values of element Simultaneous Stabilization of inversion model first row, γ after compensation2Will choosing
It is selected as that two smaller values of element Simultaneous Stabilization of inversion model secondary series after compensation can be made.In this example, it is computed, optional γ1=2,
γ2=10.It can thus be concluded that after compensation the frequency characteristic of element as shown in figure 5, as shown in Figure 5, it is compensated after, all of unit
The frequency characteristic of element is all in stability region.Therefore filtering will be fitted using the frequency characteristic of the element after compensation
Device.
(6) 100 points are taken on the frequency characteristic of element one by one after each compensation, these corresponding frequencies minute of point
Wei not ω1,ω2,……ω100, ω1=0.By taking the element of the first row first row as an example, the filter parameter asked for is:
According to this method, can obtain whole electric-wave filter matrix is:
The wave filter of step 5, the process transfer function matrix based on acquisition and design, determines PID controller parameter.
(1) robustness index Ms=1.6 is given, by the minimum time lag τ related to each output to obtain1=1, τ2=
3, the another set parameter for determining controller can be obtained:
(2) using the parameter lambda obtained1And λ2, and the γ obtained by step 41And γ2, it may be determined that final PID ginsengs
Number, bringing formula into can obtain:
Thus obtaining final centralized PID controller is:
The heating quantity of steam of step 6, the overhead reflux amount using the PID controller regulating system obtained and bottom of towe.
The present embodiment is binary distillation column Centralized Control System, and its structure is as shown in Figure 2.
The present embodiment also includes in addition to comprising a binary distillation column:Measuring transducer, for measuring overhead reflux
Amount, bottom of towe heating quantity of steam and tower top and the Methanol Molar fraction of bottom of towe extraction, and be electric signal by above parameter pick-up;
Data storage and output unit, the signal of transmitter of testing oneself in the future are stored, until meeting each passage of identification
Required data volume, then output a signal to System Discrimination unit;
System Discrimination unit, for by calculating, picking out the transmission function of each branch road;
Parameter analysis unit, the transfer function matrix to process is analyzed, and designs centralized wave filter, and utilization is given
Determine robustness index and calculate PID controller parameter;
Controller unit, centralized PID controller is constituted using the filter parameter and PID controller parameter for obtaining, from
And output control signal.
Above-described embodiment is only the displaying to spirit of the present invention, and the present invention is not limited to the present embodiment.
Claims (6)
1. a kind of binary distillation column centerized fusion method, it is characterised in that comprise the following steps:
S1, each input of measurement binary distillation column and the numerical value of each output, and will measure the value pick-up for obtaining is electric signal;
S2, the electric signal that S1 steps are obtained is stored, by after storage time t, electric signal being sent into System Discrimination mould
Block;
S3, the transmission function that each control passage is recognized using linear least squares method algorithm, set up the transmission letter of binary distillation column
Matrix number model;
S4, the transfer function model obtained according to S3 steps, are calculated the frequency characteristic of normalization inversion model, then by dividing
Analysis inversion model frequency characteristic design wave filter,
The frequency characteristic of the normalization inversion model is obtained as the following formula
Wherein, j is imaginary unit, and ω is frequency, and K is the steady-state gain inverse matrix of process,It is inversion model frequency characteristic
Matrix,It is normalization inversion model frequency characteristic matrix, subscript i and k represent the unit of the i-th row and the kth row of matrix respectively
Element.
The wave filter that S5, the transfer function model obtained according to S3 steps and S4 steps are obtained, according to given robustness and dynamic
State performance requirement, calculates PID controller parameter, and determine final centralized PID controller;
S6, adjust input using the centralized PID controller that is obtained by S5 steps.
2. binary distillation column centerized fusion method according to claim 1, it is characterised in that the storage time t is
From the step signal moment for recognizing each passage is added, the time period new steady-state value is reached to output.
3. binary distillation column centerized fusion method according to claim 1, it is characterised in that the specific steps of the S3
For:Output Y (k) of amplitude h, sampling period Ts and acquisition target according to the step signal, is calculated according to the following formula
Parameter in the transmission function of each passage
θ=(PT·P)-1PTZ
Wherein P and Z are respectively the coefficient matrixes being made up of h, Ts, Y (k), and θ is the ginseng for identifying the transmission function of each branch road
Number.
4. binary distillation column centerized fusion method according to claim 1, it is characterised in that the detailed process of S4 is:
S401, the steady-state gain inverse matrix K for asking for process,
S402, according to process transfer function matrix obtain respectively and export 1 corresponding all passages in minimum time lag τ1And with it is defeated
The minimum time lag τ gone out in 2 corresponding all passages2.,
S403, will with output 1 corresponding all passages time lag part reduce τ1, the time lag portion with 2 corresponding all passages of output
Divide and reduce τ2Afterwards, the model G for designing wave filter is obtained0,
S404, selected peak frequency ωmax, in frequency range [0, ωmax] on draw each element of normalized inversion modelFrequency characteristic, i.e. Nyquist figure, peak frequency ωmaxMeet all elements of normalization inversion model
Frequency characteristic all in the RHP of complex plane, obtain as the following formula by the frequency characteristic of normalized inversion model element
Wherein, j is imaginary unit, and ω is frequency, and K is the steady-state gain inverse matrix of process,It is inversion model frequency characteristic
Matrix,It is normalization inversion model frequency characteristic matrix, subscript i and k represent the unit of the i-th row and the kth row of matrix respectively
Element,
S405, the frequency characteristic to the element of normalization inversion model are analyzed to determine controller parameter γ by columnk's
Value, subscript k represents parameter γkThe kth column element of correspondence normalization inversion model, parameter γkSelection standard be selection one it is smaller
Value so that through the frequency characteristic of the kth row all elements for normalizing inversion model after overcompensation all in stability region, normalizing
The kth column element for changing inversion model is compensated according to the following formula
Wherein, j is imaginary unit, and ω is frequency, γkIt is corresponding controller parameter,For normalization inversion model frequency is special
Property matrix,Be the normalization inversion model frequency characteristic matrix after compensation, subscript i and k represent respectively matrix the i-th row and
The element of kth row,
S406, using the normalization inversion model after the compensation for obtaining, N number of point is taken on the frequency characteristic of each of which element,
The corresponding frequency of these points is respectively ω1, ω2... ωN, wherein, ω1=0.The method of compound curve fitting is recycled to obtain right
The transmission function of the wave filter answered
Wherein, Qik(s) be by compensating after normalization inversion model the i-th row kth column element fitting obtain wave filter, AikWith
BikIt is the corresponding parameter value of wave filter, KikIt is the i-th row kth column element of the steady-state gain inverse matrix of process, s is Laplce
Operator, parameter is asked for according to following compound curve approximating method:
Wherein, Qik(jωm) it is the i-th row kth column element, take ω in frequency on frequency characteristicmWhen corresponding plural number, | Qik(j
ωm) | it is to ask for the plural modulus value, Re (Qik(jωm)) it is to ask for the real, Im (Qik(jωm)) it is to ask for this to answer
Several imaginary parts, N is the number of the Frequency point chosen, and a, b, c, d are intermediate variable.
5. binary distillation column centerized fusion method according to claim 1, it is characterised in that the detailed process of S5 is
The value of S501, given robustness index Ms, controller regulation parameter λ1And λ2Ask for according to the following formula
S502, the sub- PID controller parameter of calculating
Wherein, KpIt is proportional gain, KIIt is storage gain, KdIt is the differential gain, TfIt is filter time constant, the tool of these parameters
Body acquiring method is as follows:
The centralized PID controller that can obtain final band filter is:
Wherein, C is total centralized PID controller, C1And C2It is the above-mentioned sub- PID controller asked for, QikIt is the filter asked for before
Ripple device.
6. a kind of binary distillation column Centralized Control System, it is characterised in that including the measuring transducer being sequentially connected and data
Also be linked in sequence after storage and output unit, data storage and output unit have System Discrimination unit, parameter analysis unit and
Controller unit, wherein
Measuring transducer, for gathering input signal and output signal;
Data storage and output unit, the data that measuring transducer transmission comes are stored, and export to System Discrimination
Unit;
System Discrimination unit, the transmission function of each passage is identified using identification algorithm, constitutes transfer function matrix;
Parameter analysis unit, the transfer function matrix to process is analyzed, and designs centralized wave filter, and using given Shandong
Rod index calculates PID controller parameter;
Controller unit, constitutes centralized PID controller, so that defeated using the filter parameter and PID controller parameter for obtaining
Go out control signal.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710163840.3A CN106896721B (en) | 2017-03-18 | 2017-03-18 | A kind of binary distillation column centerized fusion method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710163840.3A CN106896721B (en) | 2017-03-18 | 2017-03-18 | A kind of binary distillation column centerized fusion method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106896721A true CN106896721A (en) | 2017-06-27 |
CN106896721B CN106896721B (en) | 2019-11-15 |
Family
ID=59193664
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710163840.3A Active CN106896721B (en) | 2017-03-18 | 2017-03-18 | A kind of binary distillation column centerized fusion method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106896721B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112213943A (en) * | 2019-07-09 | 2021-01-12 | 北京化工大学 | Distillation tower inverted decoupling internal model control method based on equi-fractional order Butterworth filter |
CN114995155A (en) * | 2022-06-21 | 2022-09-02 | 南通大学 | Robust decoupling control system and control method for high-purity rectification process |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000257824A (en) * | 1999-03-10 | 2000-09-22 | Mitsubishi Electric Corp | Method and device for controlling combustion facility |
CN1762525A (en) * | 2005-10-14 | 2006-04-26 | 清华大学 | Rectification tower automatic control and optimization method |
CN1962014A (en) * | 2006-10-30 | 2007-05-16 | 浙江大学 | High-purity distillation general model control system and method |
CN104834217A (en) * | 2015-04-27 | 2015-08-12 | 北京化工大学 | Binary rectifying tower anti-saturation control analysis system |
CN104950725A (en) * | 2015-07-31 | 2015-09-30 | 北京化工大学 | Distributed control system for binary distillation column |
US20160169853A1 (en) * | 2014-12-15 | 2016-06-16 | Exxonmobil Research And Engineering Company | Method for determining suitability of marine fuels |
-
2017
- 2017-03-18 CN CN201710163840.3A patent/CN106896721B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000257824A (en) * | 1999-03-10 | 2000-09-22 | Mitsubishi Electric Corp | Method and device for controlling combustion facility |
CN1762525A (en) * | 2005-10-14 | 2006-04-26 | 清华大学 | Rectification tower automatic control and optimization method |
CN1962014A (en) * | 2006-10-30 | 2007-05-16 | 浙江大学 | High-purity distillation general model control system and method |
US20160169853A1 (en) * | 2014-12-15 | 2016-06-16 | Exxonmobil Research And Engineering Company | Method for determining suitability of marine fuels |
CN104834217A (en) * | 2015-04-27 | 2015-08-12 | 北京化工大学 | Binary rectifying tower anti-saturation control analysis system |
CN104950725A (en) * | 2015-07-31 | 2015-09-30 | 北京化工大学 | Distributed control system for binary distillation column |
Non-Patent Citations (5)
Title |
---|
DAZI LI: "Applications of an IMC based PID Controller tuning strategy in atmospheric and vacuum distillation units", 《NONLINEAR ANALYSIS: REAL WORLD APPLICATIONS》 * |
QIBING JIN: "Iterative Method for Frequency Domain Identification of Continuous Processes With Delay Time", 《THE CANADIAN JOURNAL OF CHEMICAL ENGINEERING》 * |
W.R.CLUETT: "Development of quality bounds for time and frequency domain models: application to the Shell distillation column", 《JOURNAL OF PROCESS CONTROL》 * |
张志强: "基于频域特性的炉温PID控制参数优化", 《冶金自动化》 * |
朱志强等: "PID 控制器的频域特性与无模型参数调节", 《控制与决策》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112213943A (en) * | 2019-07-09 | 2021-01-12 | 北京化工大学 | Distillation tower inverted decoupling internal model control method based on equi-fractional order Butterworth filter |
CN114995155A (en) * | 2022-06-21 | 2022-09-02 | 南通大学 | Robust decoupling control system and control method for high-purity rectification process |
Also Published As
Publication number | Publication date |
---|---|
CN106896721B (en) | 2019-11-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109597449A (en) | A kind of ultrasonic wave separating apparatus temprature control method neural network based | |
CN110531612A (en) | A kind of parameter tuning method of Fractional Order PID Controller | |
CN105045233B (en) | The Optimization Design of PID controller based on time metric in Power Plant Thermal system | |
CN105547325B (en) | A kind of optical fiber based on K mean cluster is used to a group temperature model coefficient and determines method | |
CN109870899A (en) | A kind of photoelectric follow-up control method based on extended state observer | |
CN106896721A (en) | A kind of binary distillation column centerized fusion method and system | |
CN106168815A (en) | A kind of acid liquor temperature control system based on Neural network PID and method | |
CN111812967A (en) | PID control parameter setting method based on stability margin and dynamic response index | |
CN106227964B (en) | Nonlinear systems with hysteresis parameter identification method based on extended state observer | |
CN104834217B (en) | Binary distillation column anti-saturation controls analysis system | |
DE102016117944A1 (en) | Servo control system with function for measuring a property of a learning controller | |
CN110378011A (en) | A kind of metal bellows hydraulic bulging process robust design method | |
DE102007021323B4 (en) | The thermal analysis apparatus | |
CN104865944B (en) | Gas separation unit control system performance estimating method based on PCA LSSVM | |
EP2884244B1 (en) | Method for operating a coriolis mass flowmeter | |
CN110034715B (en) | Voice coil motor motion control method based on disturbance estimation and related equipment | |
CN104898587B (en) | Industrial process modeling System and method for based on parallel diffused intelligent search algorithm | |
CN104793496B (en) | Two inputs two export the decoupling and controlling system of polymer reactor | |
CN108345574A (en) | Related dual data stream abnormality detection and modified method | |
CN108170637A (en) | A kind of transfer function model discrimination method with derivative characteristic process | |
CN113244647A (en) | Rectifying tower VOFFLC temperature control method and system based on matrix decoupling | |
CN100441552C (en) | Process for advanced controlling rectifying apparatus of butadiene | |
CN106950835A (en) | Succinct Robust Two-Degree-of-Freedom proportional plus integral control method | |
CN114706311B (en) | Control method of multivariable control system | |
CN109992887B (en) | Anti-interference control method and system for binary distillation tower |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |