CN100510590C - General model connecting system of air separation tower and method thereof - Google Patents

Abstract

A general mode control system of space dividing tower comprises an on-site intelligent device connected with the space dividing tower, a data storage device for storing history data, and a superior machine, wherein the intelligent device, data storage device and superior are serially connected, the superior machine is general mode controller comprising a component controller and a general mode controller. The component controller comprises a checking module, an I/O module, and a component analyzing module. The general mode controller comprises an I/O module, a secondary curvature simulating module, a general mode controlling calculating module and a control output module. The invention also provides a relative control method. The invention can be applicable for the dynamic, non-liner and couple properties of space dividing operation, to obtain better control effect.

Description

The general model control system of air separation column and method

(1) technical field

The present invention relates to the control system and the method design field of air separation column, especially, relate to a kind of general model control system and method for air separation column.

(2) background technology

Air separation unit separates the sky branch exactly, and obtains the device of high-purity industrial gasses such as oxygen, nitrogen, argon.It is the supportive unit operations of numerous industries that concern the life-blood of the national economy, as chemical industry, metallurgy, electronics, the energy, Aero-Space, food and drink etc., belong to national substantial equipment, its development scale and technology status are to weigh the industry of a country and an important symbol of development in science and technology level.

Empty branch operation is one and relates to low temperature, many equipment, long flow process, complicated operation, the exigent complex process of safety in production.In the production, the purity of oxygen, nitrogen, argon product often requires up to more than 99%, belong to high-purity distillation control problem, stationarity to the air separation column operation requires very high, and the high-purity distillation process is because the coupling between dynamic characteristic, strong nonlinear and the loop of the complexity that it showed, and traditional is difficult to it is controlled effect preferably as Linear Control schemes such as PID.

(3) summary of the invention

In order to overcome coupling between empty dynamic characteristic, strong nonlinear and the loop of dividing operation of can not adapting to of existing air separation column control scheme, can not to obtain the deficiency of good control effect, the invention provides a kind of coupling that can adapt between empty dynamic characteristic, strong nonlinear and the loop of dividing operation, can access the general model control system and the method for the air separation column of good control effect.

The technical solution adopted for the present invention to solve the technical problems is:

A kind of general model control system of air separation column, comprise and the direct-connected field intelligent instrument of air separation column, the data storage device that is used for storing history data and host computer, intelligence instrument, data storage device and host computer link to each other successively, described host computer is the universal model controller, described universal model controller comprises component deduction control section and universal model control section, described component infers that control section comprises: the instrumentation module, comprise detector unit and pressure detecting element, be used to detect the temperature and pressure of the last tower of air separation column; The I/O component module is used for the transmission between controller inside and controller and DCS of the signal of telecommunication and data-signal, and the component inference module is used for inferring component according to detecting the temperature and the pressure data that obtain, and its formula is (1), (2):

${\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_1=\frac{\mathrm{\&alpha;}}{\mathrm{\&alpha;}-1}-\frac{{10}^{(a-\frac{\mathrm{<figure-callout\; id="1"\; label="b\; T"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">b</figure-callout>}}{T_{}})}}{(\mathrm{\&alpha;}-1)P}---\left(1\right)$

$\mathrm{Xn}=\frac{\mathrm{P\&alpha;}}{{10}^{(a-\frac{T_n+c}{b})}(\mathrm{\&alpha;}-1)}-\frac{1}{\mathrm{\&alpha;}-1}---\left(2\right)$

Wherein, Y ₁Be the component of nitrogen in the nitrogen product in the air separation column, Xn is the component of nitrogen in the liquid oxygen product, and P is last tower pressure, T ₁, T _nBe respectively column overhead, column bottom temperature, α is a relative volatility, and a, b, c are the Peter Antonie constant;

Described universal model control section comprises: the I/O component module is used for the inside of universal model controller and the signal of telecommunication between controller and the DCS, the transmission of data-signal; Conic section fits module, and the historical data that is used for that the data storage device is obtained is carried out conic section and fitted, and obtains the nitrogen component Y in the product nitrogen gas ₁About the following quadratic equation of the capacity of returns R of the supreme tower of tower liquid nitrogen liquid oxygen: Y ₁(k)=a ₁* R (k) ², and the nitrogen component Xn in the product liquid oxygen is about the quadratic equation of the output flow L of liquid oxygen product: X _n(k)=a ₂* L (k) ²

Universal model control computing module is used for the Y that obtains by the component inference section ₁With the Xn value, calculate the capacity of returns R of tower liquid nitrogen liquid air under the current control variables and the value that goes up the output flow L of tower liquid oxygen product, its formula is (3), (4), (5):

$\mathrm{\&Delta;R}\left(k\right)=\frac{1}{2\sqrt{a_1}}\&times;\frac{(1-T\&times;K_1)\&times;{\stackrel{.}{Y}}_1\left(k\right)+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\&times;T\&times;({\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_{1\mathrm{set}}-Y_1\left(k\right))}{\sqrt{Y_1\left(k\right)+T(K_1({\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_{1\mathrm{set}}-Y_1\left(k\right))+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\mathrm{\&Sigma;}({\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_{1\mathrm{set}}-Y_1\left(k\right))}}---\left(3\right)$

$\mathrm{\&Delta;L}\left(k\right)=\frac{1}{2\sqrt{{\mathrm{<figure-callout\; id="2"\; label="a"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">a</figure-callout>}}_2}}\&times;\frac{(1-T\&times;K_1)\&times;{\stackrel{.}{X}}_n\left(k\right)+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\&times;T\&times;(X_{\mathrm{nset}}-X_n\left(k\right))}{\sqrt{X_n\left(k\right)+T(K_1(X_{\mathrm{nset}}-X_n\left(k\right))+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\mathrm{\&Sigma;}(X_{\mathrm{nset}}-X_n\left(k\right))}}---\left(4\right)$

R(k)＝R(k-1)+ΔR(k)

(5)

L(k)＝L(k-1)+ΔL(k)

Wherein, T is the sampling period, K ₁, K ₂Be controller adjustable parameter, Y _1set, X _NsetBe respectively Y ₁With X _nSetting value; The control output module is used for the R (k) that will calculate, and the data-signal of L (k) outputs to air separation column.

As preferred a kind of scheme: described general model control system also comprises the DCS system, and described DCS system is made of data-interface, control station and historical data base, and described data storage device is the historical data base of DCS system.

As preferred another kind of scheme: described field intelligent instrument, DCS system, universal model controller connect successively by fieldbus

As preferred another scheme: described universal model controller also comprises human-computer interface module, is used for the control variables R (k) that will calculate, the value of L (k), and detect the Y that obtains ₁, X _nValue on the man-machine interface of controller, show.

The control method that the general model control system of the described air separation column of a kind of usefulness is realized, described control method may further comprise the steps:

(1) determines the bi-component setting value Y of air separation column _1set, X _Nset, and sampling period T;

(2) after doing step test around the setting value, from data storage device, obtain historical data, fit module by conic section and obtain, the nitrogen component Y in the nitrogen ₁About the following quadratic equation of the capacity of returns R of the supreme tower of tower liquid nitrogen liquid oxygen: Y ₁(k)=a ₁* R (k) ², and the nitrogen component X in the product liquid oxygen _nQuadratic equation about the output flow L of liquid oxygen product: X _n(k)=a ₂* L (k) ²

(3) each sampling instant KT infers component according to detecting the temperature and the pressure data that obtain, and its formula is (1), (2):

${\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_1=\frac{\mathrm{\&alpha;}}{\mathrm{\&alpha;}-1}-\frac{{10}^{(a-\frac{\mathrm{<figure-callout\; id="1"\; label="b\; T"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">b</figure-callout>}}{T_{}})}}{(\mathrm{\&alpha;}-1)P}---\left(1\right)$

$\mathrm{Xn}=\frac{\mathrm{P\&alpha;}}{{10}^{(a-\frac{T_n+c}{b})}(\mathrm{\&alpha;}-1)}-\frac{1}{\mathrm{\&alpha;}-1}---\left(2\right)$

Wherein, Y ₁Be the component of nitrogen in the nitrogen product in the air separation column, Xn is the component of nitrogen in the liquid oxygen product, and P is last tower pressure, T ₁, T _nBe respectively column overhead, column bottom temperature, α is a relative volatility, and a, b, c are the Peter Antonie constant;

(4) controller reads Y from data storage device ₁With the value of Xn as input, the s operation control variable R, the value of L, its formula is (3), (4), (5)::

$\mathrm{\&Delta;R}\left(k\right)=\frac{1}{2\sqrt{a_1}}\&times;\frac{(1-T\&times;K_1)\&times;{\stackrel{.}{Y}}_1\left(k\right)+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\&times;T\&times;({\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_{1\mathrm{set}}-Y_1\left(k\right))}{\sqrt{Y_1\left(k\right)+T(K_1({\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_{1\mathrm{set}}-Y_1\left(k\right))+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\mathrm{\&Sigma;}({\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_{1\mathrm{set}}-Y_1\left(k\right))}}---\left(3\right)$

$\mathrm{\&Delta;L}\left(k\right)=\frac{1}{2\sqrt{{\mathrm{<figure-callout\; id="2"\; label="a"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">a</figure-callout>}}_2}}\&times;\frac{(1-T\&times;K_1)\&times;{\stackrel{.}{X}}_n\left(k\right)+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\&times;T\&times;(X_{\mathrm{nset}}-X_n\left(k\right))}{\sqrt{X_n\left(k\right)+T(K_1(X_{\mathrm{nset}}-X_n\left(k\right))+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\mathrm{\&Sigma;}(X_{\mathrm{nset}}-X_n\left(k\right))}}---\left(4\right)$

R(k)＝R(k-1)+ΔR(k)

(5)

L(k)＝L(k-1)+ΔL(k)

Wherein, T is the sampling period, K ₁, K ₂Be controller adjustable parameter, Y _1set, X _NsetBe respectively Y ₁Setting value with Xn;

(5) with R (k), the data-signal of L (k) returns to air separation column.

As preferred a kind of scheme: described control method also comprises: (6), fall into a trap in described (4) and to have calculated control variables R (k), and the value of L (k), and with it and detect the Y that obtains ₁, Xn value on the man-machine interface of controller, show.

As preferred another kind of scheme: described data storage device is the historical data base of DCS system, described DCS system is made of data-interface, control station and historical data base, in described (6), data are passed to the DCS system, and at the control station procedure for displaying state of DCS.

Technical conceive of the present invention is: the component Y that adopts nitrogen in the nitrogen product ₁, the component Xn of nitrogen is a controlled variable in the liquid oxygen product, the flow of the liquid oxygen of the supreme tower of following tower and the capacity of returns of liquid air, liquid oxygen product is the control corresponding variable.

Infer control section, be used to solve the difficult problem that the industry spot product component can not directly be measured,, can eliminate greatly and measure hysteresis and have stronger reliability relatively with respect to the chromatographic way of online applicable industry.The universal model control section is used to use the universal model control algolithm to obtain the value of real-time control variables.

General model control system has effectively solved air separation column in high-purity problem of strong nonlinearity down, realized quiet run to bi-component control at the bottom of the high-purity distillation process cat head tower, and had zero deflection characteristic faster, more traditional PID control system is significantly improved on dynamic property, so very large application prospect is arranged.

Select suitable universal model controller parameter K1 and K2.Accompanying drawing 4 is choosing with reference to figure of universal model controller parameter K1, K2, and it is a system:

$\frac{x}{x^{*}}=\frac{2\mathrm{\&tau;\&xi;s}+1}{{\mathrm{\&tau;}}^2{\mathrm{<figure-callout\; id="2"\; label="s"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">s</figure-callout>}}^2+2\mathrm{\&tau;\&xi;s}+1}$

Be similar to the standard step response curve of the x/x* of classical second-order system about t/ τ.

It is as follows specifically to choose step:

1) from accompanying drawing 4, selects suitable ξ according to needed desirable dynamic characteristic.

2) by selecting the dynamic performance index under this ξ to calculate the τ value.

3) according to formula $\mathrm{\&tau;}=\frac{1}{\sqrt{{\mathrm{<figure-callout\; id="2"\; label="k"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">k</figure-callout>}}_2}},$ $\mathrm{\&xi;}=\frac{{\mathrm{<figure-callout\; id="1"\; label="k"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">k</figure-callout>}}_1}{2\sqrt{{\mathrm{<figure-callout\; id="2"\; label="k"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">k</figure-callout>}}_2}}$ The anti-calculating K of separating ₁, K ₂

Beneficial effect of the present invention mainly shows: 1, can adapt to the empty coupling that divides between dynamic characteristic, strong nonlinear and the loop of operating, realize the quiet run to bi-component control at the bottom of the high-purity distillation process cat head tower; 2, can access good control effect.

(4) description of drawings

Fig. 1 is the hardware connection layout of the general model control system of air separation column proposed by the invention.

Fig. 2 is the theory diagram of universal model controller of the present invention.

Fig. 3 is the on-the-spot connection layout of the general model control system of air separation column proposed by the invention.

Fig. 4 is that universal model controller parameter K1, K2 choose with reference to figure.

(5) specific embodiment

Below in conjunction with accompanying drawing the present invention is further described.The embodiment of the invention is used for the present invention that explains, rather than limits the invention, and in the protection domain of spirit of the present invention and claim, any modification and change to the present invention makes all fall into protection scope of the present invention.

Embodiment 1

With reference to Fig. 1～Fig. 4, a kind of general model control system of air separation column, comprise and air separation column 1 direct-connected field intelligent instrument 2, the data storage device and the host computer 6 that are used for storing history data, intelligence instrument 2, data storage device and host computer 6 link to each other successively, described host computer 6 is the universal model controller, described universal model controller comprises component deduction control section and universal model control section, described component infers that control section comprises: instrumentation module 7, comprise detector unit and pressure detecting element, be used to detect the temperature and pressure of the last tower of air separation column; I/O component module 9 is used for the transmission between controller inside and controller and DCS of the signal of telecommunication and data-signal, and component inference module 10 is used for inferring component according to detecting the temperature and the pressure data that obtain, and its formula is (1), (2):

${\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_1=\frac{\mathrm{\&alpha;}}{\mathrm{\&alpha;}-1}-\frac{{10}^{(a-\frac{\mathrm{<figure-callout\; id="1"\; label="b\; T"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">b</figure-callout>}}{T_{}})}}{(\mathrm{\&alpha;}-1)P}---\left(1\right)$

$\mathrm{Xn}=\frac{\mathrm{P\&alpha;}}{{10}^{(a-\frac{T_n+c}{b})}(\mathrm{\&alpha;}-1)}-\frac{1}{\mathrm{\&alpha;}-1}---\left(2\right)$

Wherein, Y ₁Be the component of nitrogen in the nitrogen product in the air separation column, Xn is the component of nitrogen in the liquid oxygen product, and P is last tower pressure, T ₁, T _nBe respectively column overhead, column bottom temperature, α is a relative volatility, and a, b, c are the Peter Antonie constant;

Described universal model control section 11 comprises: the I/O component module is used for the inside of universal model controller and the signal of telecommunication between controller and the DCS, the transmission of data-signal; Conic section fits module, and the historical data that is used for that the data storage device is obtained is carried out conic section and fitted, and obtains the nitrogen component Y in the product nitrogen gas ₁About the following quadratic equation of the capacity of returns R of the supreme tower of tower liquid nitrogen liquid oxygen: Y ₁(k)=a ₁* R (k) ², and the nitrogen component Xn in the product liquid oxygen is about the quadratic equation of the output flow L of liquid oxygen product: X _n(k)=a ₂* L (k) ²Universal model control computing module is used for the Y that obtains by the component inference section ₁With the Xn value, calculate the capacity of returns R of tower liquid nitrogen liquid air under the current control variables and the value that goes up the output flow L of tower liquid oxygen product, its formula is (3), (4), (5):

$\mathrm{\&Delta;R}\left(k\right)=\frac{1}{2\sqrt{a_1}}\&times;\frac{(1-T\&times;K_1)\&times;{\stackrel{.}{Y}}_1\left(k\right)+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\&times;T\&times;({\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_{1\mathrm{set}}-Y_1\left(k\right))}{\sqrt{Y_1\left(k\right)+T(K_1({\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_{1\mathrm{set}}-Y_1\left(k\right))+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\mathrm{\&Sigma;}({\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_{1\mathrm{set}}-Y_1\left(k\right))}}---\left(3\right)$

$\mathrm{\&Delta;L}\left(k\right)=\frac{1}{2\sqrt{a_2}}\&times;\frac{(1-T\&times;K_1)\&times;{\stackrel{.}{X}}_n\left(k\right)+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\&times;T\&times;(X_{\mathrm{nset}}-X_n\left(k\right))}{\sqrt{X_n\left(k\right)+T(K_1(X_{\mathrm{nset}}-X_n\left(k\right))+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\mathrm{\&Sigma;}(X_{\mathrm{nset}}-X_n\left(k\right))}}---\left(4\right)$

R(k)＝R(k-1)+ΔR(k)

(5)

L(k)＝L(k-1)+ΔL(k)

Wherein, T is the sampling period, K ₁, K ₂Be controller adjustable parameter, Y _1set, X _NsetBe respectively Y ₁Setting value with Xn; The control output module is used for the R (k) that will calculate, and the data-signal of L (k) outputs to air separation column.

Described general model control system also comprises DCS system 12, and described DCS system 12 is made of data-interface 3, control station 4 and historical data base 5, and described data storage device is the historical data base 5 of DCS system.Described field intelligent instrument 2, DCS system, universal model controller 6 connect successively by fieldbus.

With reference to Fig. 1, the general model control system of the air separation column of present embodiment, comprise the field intelligent instrument 2, DCS system and the universal model controller 6 that link to each other with on-the-spot air separation column 1, described DCS system is made of data-interface 3, control station 4 and historical data base 5; On-the-spot air separation column object 1, intelligence instrument 2, DCS system, universal model controller 6 connect successively by fieldbus.

The general model control system hardware structure diagram of the air separation column of present embodiment as shown in Figure 1, the core of described general model control system is a universal model controller 6, comprises in addition: field intelligent instrument 2, DCS system and fieldbus.On-the-spot air separation column 1, intelligence instrument 2, DCS system, universal model controller 6 link to each other successively by fieldbus, and uploading of information of realization assigned.General model control system in time obtains the value of the control variables of current time by industry spot data detected and that extract from historical data base 5, and returns to first floor system, in time system is dynamically made a response.

The theory diagram of the universal model controller of the air separation column of present embodiment as shown in Figure 2, the universal model controller of described air separation column comprises:

Infer control section, be used to solve the difficult problem that the industry spot product component can not directly be measured,, can eliminate greatly and measure hysteresis and have stronger reliability relatively with respect to the chromatographic way of online applicable industry.

1) the instrumentation module 7: comprise detector unit, can adopt the thermojunction type temperature transmitter, and pressure detecting element, can adopt the pressure resistance type transmitter.

2) the I/O component module 9: be used for the transmission between controller inside and controller and DCS of the signal of telecommunication and data-signal.

3) the component inference module 10: be used for inferring component according to detecting the temperature and the pressure data that obtain.Its formula is (1), (2):

${\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_1=\frac{\mathrm{\&alpha;}}{\mathrm{\&alpha;}-1}-\frac{{10}^{(a-\frac{\mathrm{<figure-callout\; id="1"\; label="b\; T"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">b</figure-callout>}}{T_{}})}}{(\mathrm{\&alpha;}-1)P}---\left(1\right)$

$\mathrm{Xn}=\frac{\mathrm{P\&alpha;}}{{10}^{(a-\frac{T_n+c}{b})}(\mathrm{\&alpha;}-1)}-\frac{1}{\mathrm{\&alpha;}-1}---\left(2\right)$

Wherein P is last tower pressure, T ₁, T _nBe respectively column overhead, column bottom temperature, α is a relative volatility, and a, b, c are the Peter Antonie constant;

Universal model control section 11 is used to use the universal model control algolithm to obtain the value of real-time control variables.

1) I/O element: be used for the inside of universal model control and the signal of telecommunication between controller and the DCS, the transmission of data-signal.

2) conic section fits module, is used for that test or the data that obtain of DCS historical data base are carried out conic section and fits, thereby obtain nitrogen component Y in the product nitrogen gas ₁About the following quadratic equation of the capacity of returns R of the supreme tower of tower liquid nitrogen liquid oxygen: Y ₁(k)=a ₁* R (k) ², and the nitrogen component Xn in the product liquid oxygen is about the quadratic equation of the output flow L of liquid oxygen product: X _n(k)=a ₂* L (k) ²

3) universal model control computing module is used for the Y that obtains by the component inference section ₁With the Xn value, calculate the capacity of returns R of tower liquid nitrogen liquid air under the current control variables and the value that goes up the output flow L of tower liquid oxygen product, its formula is (3), (4), (5):

$\mathrm{\&Delta;R}\left(k\right)=\frac{1}{2\sqrt{a_1}}\&times;\frac{(1-T\&times;K_1)\&times;{\stackrel{.}{Y}}_1\left(k\right)+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\&times;T\&times;({\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_{1\mathrm{set}}-Y_1\left(k\right))}{\sqrt{Y_1\left(k\right)+T(K_1({\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_{1\mathrm{set}}-Y_1\left(k\right))+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\mathrm{\&Sigma;}({\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_{1\mathrm{set}}-Y_1\left(k\right))}}---\left(3\right)$

$\mathrm{\&Delta;L}\left(k\right)=\frac{1}{2\sqrt{a_2}}\&times;\frac{(1-T\&times;K_1)\&times;{\stackrel{.}{X}}_n\left(k\right)+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\&times;T\&times;(X_{\mathrm{nset}}-X_n\left(k\right))}{\sqrt{X_n\left(k\right)+T(K_1(X_{\mathrm{nset}}-X_n\left(k\right))+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\mathrm{\&Sigma;}(X_{\mathrm{nset}}-X_n\left(k\right))}}---\left(4\right)$

R(k)＝R(k-1)+ΔR(k)

(5)

L(k)＝L(k-1)+ΔL(k)

Wherein, T is the sampling period, K ₁, K ₂Be controller adjustable parameter, Y _1set, X _NsetBe respectively Y ₁Setting value with Xn.

The universal model controller of described air separation column also comprises human-computer interface module 8, is used for the demonstration of historical data and system's current state, and the operation of control system parameter selection etc.

The on-the-spot connection layout of the general model control system of the air separation column of present embodiment as shown in Figure 3, system adopts the component Y that goes up nitrogen in the tower 14 top nitrogen products ₁, the component Xn that goes up nitrogen in the tower 1 bottom liquid oxygen product is controlled variable, the flow of the liquid oxygen of following tower 15 supreme towers 14 and the capacity of returns of liquid air, liquid oxygen product is the control corresponding variable.Connect a detector unit TT and pressure detecting element PT respectively and be delivered to upper system at the bottom of the last tower 14 cat head towers, the universal model controller by the data computation current time of on-the-spot and historical data base the control variables value and pass to down layer system, the scene changes the value of control variables by the change valve opening by flow controller FC.

The universal model control method of described air separation column realizes according to following steps:

1, system initialization

(1) in universal model controller 6, sets the bi-component setting value Y of air separation column _1set, X _Nset, sampling period T and sampling period among the DCS is set.

(2) after doing the open loop step test around the setting value, from DCS historical data base 5, obtain historical data (or directly from historical data base 5, obtaining), fit module by conic section and obtain, the nitrogen component Y in the nitrogen ₁About the following quadratic equation of the capacity of returns R of the supreme tower of tower liquid nitrogen liquid oxygen: Y ₁(k)=a ₁* R (k) ², and the nitrogen component Xn in the product liquid oxygen is about the quadratic equation of the output flow L of liquid oxygen product: X _n(k)=a ₂* L (k) ²

(3) select suitable universal model controller 6 parameter K 1 and K2.

Accompanying drawing 4 is choosing with reference to figure of universal model controller parameter K1, K2, and it is a system:

$\frac{x}{x^{*}}=\frac{2\mathrm{\&tau;\&xi;s}+1}{{\mathrm{\&tau;}}^2{\mathrm{<figure-callout\; id="2"\; label="s"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">s</figure-callout>}}^2+2\mathrm{\&tau;\&xi;s}+1}$

Be similar to the standard step response curve of the x/x* of classical second-order system about t/ τ.

It is as follows specifically to choose step:

1) from accompanying drawing 4, selects suitable ξ according to needed desirable dynamic characteristic.

2) by selecting the dynamic performance index under this ξ to calculate the τ value.

3) according to formula $\mathrm{\&tau;}=\frac{1}{\sqrt{{\mathrm{<figure-callout\; id="2"\; label="k"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">k</figure-callout>}}_2}},$ $\mathrm{\&xi;}=\frac{{\mathrm{<figure-callout\; id="1"\; label="k"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">k</figure-callout>}}_1}{2\sqrt{{\mathrm{<figure-callout\; id="2"\; label="k"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">k</figure-callout>}}_2}}$ The anti-calculating K of separating ₁, K ₂

2, system puts into operation.

(1) each DCS sampling instant, intelligence instrument 2 detects temperature, the pressure data of on-the-spot air separation column 1 and is sent in the DCS historical data base 5;

(2) each controller sampling instant, universal model controller 6 reads the temperature and pressure data from DCS historical data base 5, calculate current time controlled variable Y by component inference module 10 ₁, the value of Xn, its formula is (1), (2);

${\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_1=\frac{\mathrm{\&alpha;}}{\mathrm{\&alpha;}-1}-\frac{{10}^{(a-\frac{\mathrm{<figure-callout\; id="1"\; label="b\; T"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">b</figure-callout>}}{T_{}})}}{(\mathrm{\&alpha;}-1)P}---\left(1\right)$

$\mathrm{Xn}=\frac{\mathrm{P\&alpha;}}{{10}^{(a-\frac{T_n+c}{b})}(\mathrm{\&alpha;}-1)}-\frac{1}{\mathrm{\&alpha;}-1}---\left(2\right)$

Wherein, Y ₁Be the component of nitrogen in the nitrogen product in the air separation column, Xn is the component of nitrogen in the liquid oxygen product, and P is last tower pressure, T ₁, T _nBe respectively column overhead, column bottom temperature, α is a relative volatility, and a, b, c are the Peter Antonie constant;

(3) Y by obtaining from component inference module 10 ₁, Xn value, the computing by universal model control module 11 obtains the control variables R of current time, the value of L, its formula is (3), (4), (5):

$\mathrm{\&Delta;R}\left(k\right)=\frac{1}{2\sqrt{a_1}}\&times;\frac{(1-T\&times;K_1)\&times;{\stackrel{.}{Y}}_1\left(k\right)+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\&times;T\&times;({\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_{1\mathrm{set}}-Y_1\left(k\right))}{\sqrt{Y_1\left(k\right)+T(K_1({\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_{1\mathrm{set}}-Y_1\left(k\right))+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\mathrm{\&Sigma;}({\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_{1\mathrm{set}}-Y_1\left(k\right))}}---\left(3\right)$

$\mathrm{\&Delta;L}\left(k\right)=\frac{1}{2\sqrt{a_2}}\&times;\frac{(1-T\&times;K_1)\&times;{\stackrel{.}{X}}_n\left(k\right)+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\&times;T\&times;(X_{\mathrm{nset}}-X_n\left(k\right))}{\sqrt{X_n\left(k\right)+T(K_1(X_{\mathrm{nset}}-X_n\left(k\right))+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\mathrm{\&Sigma;}(X_{\mathrm{nset}}-X_n\left(k\right))}}---\left(4\right)$

R(k)＝R(k-1)+ΔR(k)

(5)

L(k)＝L(k-1)+ΔL(k)

Wherein, T is the sampling period, K ₁, K ₂Be controller adjustable parameter, Y _1set, X _NsetBe respectively Y ₁Setting value with Xn;

(4) with R (k), the data-signal of L (k) returns to the DCS system, and acts on air separation column.

(5) result is delivered on the display module of each level system and show, make things convenient for the engineer in time process dynamically to be reacted and operate, comprise human-computer interface module 8, DCS system active station 4 and the operator station of universal model controller.

Embodiment 2

With reference to Fig. 1～Fig. 4, the control method that the general model control system of the described air separation column of a kind of usefulness is realized, described control method may further comprise the steps:

(1) determines the bi-component setting value Y of air separation column _1set, X _Nset, and sampling period T;

(2) after doing step test around the setting value, from data storage device, obtain historical data, fit module by conic section and obtain, the nitrogen component Y in the nitrogen ₁About the following quadratic equation of the capacity of returns R of the supreme tower of tower liquid nitrogen liquid oxygen: Y ₁(k)=a ₁* R (k) ², and the nitrogen component Xn in the product liquid oxygen is about the quadratic equation of the output flow L of liquid oxygen product: X _n(k)=a ₂* L (k) ²

(3) each sampling instant KT infers component according to detecting the temperature and the pressure data that obtain, and its formula is (1), (2):

${\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_1=\frac{\mathrm{\&alpha;}}{\mathrm{\&alpha;}-1}-\frac{{10}^{(a-\frac{\mathrm{<figure-callout\; id="1"\; label="b\; T"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">b</figure-callout>}}{T_{}})}}{(\mathrm{\&alpha;}-1)P}---\left(1\right)$

$\mathrm{Xn}=\frac{\mathrm{P\&alpha;}}{{10}^{(a-\frac{T_n+c}{b})}(\mathrm{\&alpha;}-1)}-\frac{1}{\mathrm{\&alpha;}-1}---\left(2\right)$

Wherein, Y ₁Be the component of nitrogen in the nitrogen product in the air separation column, Xn is the component of nitrogen in the liquid oxygen product, and P is last tower pressure, T ₁, T _nBe respectively tower, cat head column bottom temperature, α is a relative volatility, and a, b, c are the Peter Antonie constant;

(4) controller reads Y from data storage device ₁With the value of Xn as input, the s operation control variable R, the value of L, its formula is (3), (4), (5)::

$\mathrm{\&Delta;R}\left(k\right)=\frac{1}{2\sqrt{a_1}}\&times;\frac{(1-T\&times;K_1)\&times;{\stackrel{.}{Y}}_1\left(k\right)+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\&times;T\&times;({\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_{1\mathrm{set}}-Y_1\left(k\right))}{\sqrt{Y_1\left(k\right)+T(K_1({\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_{1\mathrm{set}}-Y_1\left(k\right))+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\mathrm{\&Sigma;}({\mathrm{<figure-callout\; id="1"\; label="Y"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">Y</figure-callout>}}_{1\mathrm{set}}-Y_1\left(k\right))}}---\left(3\right)$

$\mathrm{\&Delta;L}\left(k\right)=\frac{1}{2\sqrt{a_2}}\&times;\frac{(1-T\&times;K_1)\&times;{\stackrel{.}{X}}_n\left(k\right)+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\&times;T\&times;(X_{\mathrm{nset}}-X_n\left(k\right))}{\sqrt{X_n\left(k\right)+T(K_1(X_{\mathrm{nset}}-X_n\left(k\right))+{\mathrm{<figure-callout\; id="2"\; label="K"\; filenames="C200610155486E00191.png,C200610155486E00202.png"\; state="\{\{state\}\}">K</figure-callout>}}_2\mathrm{\&Sigma;}(X_{\mathrm{nset}}-X_n\left(k\right))}}---\left(4\right)$

R(k)＝R(k-1)+ΔR(k)

(5)

L(k)＝L(k-1)+ΔL(k)

Wherein, T is the sampling period, K ₁, K ₂Be controller adjustable parameter, Y _1set, X _NsetBe respectively Y ₁Setting value with Xn;

(5) with R (k), the data-signal of L (k) returns to air separation column.

Described control method also comprises: (6), fall into a trap in described (4) and to have calculated control variables R (k), and the value of L (k), and with it and detect the Y that obtains ₁, Xn value on the man-machine interface of controller, show.Described data storage device is the historical data base of DCS system 12, and described DCS system is made of data-interface 3, control station 4 and historical data base 5, in described (6), data is passed to the DCS system, and at the control station procedure for displaying state of DCS.

Claims (7)

1, a kind of general model control system of air separation column, comprise and the direct-connected field intelligent instrument of air separation column, the data storage device that is used for storing history data and host computer, intelligence instrument, data storage device and host computer link to each other successively, it is characterized in that: described general model control system also comprises the DCS system, described host computer is the universal model controller, described universal model controller comprises component deduction control section and universal model control section

Described component infers that control section comprises:

The instrumentation module comprises detector unit and pressure detecting element, is used to detect the temperature and pressure of the last tower of air separation column;

The I/O component module is used for the transmission between controller inside and controller and described DCS system of the signal of telecommunication and data-signal,

The component inference module is used for inferring component according to detecting the temperature and the pressure data that obtain, and its formula is (1), (2):

$Y_1=\frac{\mathrm{\&alpha;}}{\mathrm{\&alpha;}-1}-\frac{{10}^{(a-\frac{b}{T_1+c})}}{(\mathrm{\&alpha;}-1)P}---\left(1\right)$

$\mathrm{Xn}=\frac{\mathrm{P\&alpha;}}{{10}^{(a-\frac{T_n+c}{b})}(\mathrm{\&alpha;}-1)}-\frac{1}{\mathrm{\&alpha;}-1}---\left(2\right)$

Wherein, Y ₁Be the component of nitrogen in the nitrogen product in the air separation column, Xn is the component of nitrogen in the liquid oxygen product, and P is last tower pressure, T ₁, T _nBe respectively column overhead and column bottom temperature, α is a relative volatility, and a, b, c are the Peter Antonie constant;

Described universal model control section comprises:

The I/O component module is used for the signal of telecommunication between the inside of universal model controller and controller and the described DCS system, the transmission of data-signal;

Conic section fits module, and the historical data that is used for that the data storage device is obtained is carried out conic section and fitted, and obtains the nitrogen component Y in the product nitrogen gas ₁About the following quadratic equation of the capacity of returns R of the supreme tower of tower liquid nitrogen liquid oxygen: Y ₁(k)=a ₁* R (k) ², and the nitrogen component Xn in the product liquid oxygen is about the quadratic equation of the output flow L of liquid oxygen product: X _n(k)=a ₂* L (k) ²

Universal model control computing module is used for the Y that obtains by the component inference section ₁With the Xn value, calculate the capacity of returns R of tower liquid nitrogen liquid air under the current control variables and the value that goes up the output flow L of tower liquid oxygen product, its formula is (3), (4), (5):

$\mathrm{\&Delta;R}\left(k\right)=\frac{1}{2\sqrt{a_1}}\&times;\frac{(1-T\&times;K_1)\&times;{\stackrel{.}{Y}}_1\left(k\right)+K_2\&times;T\&times;(Y_{1\mathrm{set}}-Y_1\left(k\right))}{\sqrt{Y_1\left(k\right)+T(K_1(Y_{1\mathrm{set}}-Y_1\left(k\right))+K_2\mathrm{\&Sigma;}(Y_{1\mathrm{set}}-Y_1\left(k\right))}}---\left(3\right)$

$\mathrm{\&Delta;L}\left(k\right)=\frac{1}{2\sqrt{a_2}}\&times;\frac{(1-T\&times;K_1)\&times;{\stackrel{.}{X}}_n\left(k\right)+K_2\&times;T\&times;(X_{\mathrm{nset}}-X_n\left(k\right))}{\sqrt{X_n\left(k\right)+T(K_1(X_{\mathrm{nset}}-X_n\left(k\right))+K_2\mathrm{\&Sigma;}(X_{\mathrm{nset}}-X_n\left(k\right))}}---\left(4\right)$

R(k)＝R(k-1)+ΔR(k)

(5)

L(k)＝L(k-1)+ΔL(k)

Wherein, T is the sampling period, K ₁, K ₂Be controller adjustable parameter, Y _1set, X _NsetBe respectively Y ₁Setting value with Xn;

The control output module is used for the R (k) that will calculate, and the data-signal of L (k) outputs to air separation column.

2, the general model control system of air separation column as claimed in claim 1 is characterized in that: described DCS system is made of data-interface, control station and historical data base, and described data storage device is the historical data base of DCS system.

3, the general model control system of air separation column as claimed in claim 2 is characterized in that: described field intelligent instrument, DCS system, universal model controller connect successively by fieldbus.

4, as the general model control system of the described air separation column of one of claim 1～3, it is characterized in that: described universal model controller also comprises human-computer interface module, is used for the control variables R (k) that will calculate, the value of L (k), and detect the Y that obtains ₁, Xn value on the man-machine interface of controller, show.

5, the control method of the general model control system of a kind of usefulness air separation column as claimed in claim 1 realization, it is characterized in that: described control method may further comprise the steps:

1) determines the bi-component setting value Y of air separation column _1set, X _Nset, and sampling period T;

2) after doing step test around the setting value, from data storage device, obtain historical data, fit module by conic section and obtain, the nitrogen component Y in the nitrogen ₁About the following quadratic equation of the capacity of returns R of the supreme tower of tower liquid nitrogen liquid oxygen: Y ₁(k)=a ₁* R (k) ², and the nitrogen component Xn in the product liquid oxygen is about the quadratic equation of the output flow L of liquid oxygen product: X _n(k)=a ₂* L (k) ²

3) each sampling instant KT infers component according to detecting the temperature and the pressure data that obtain, and its formula is (1), (2):

$Y_1=\frac{\mathrm{\&alpha;}}{\mathrm{\&alpha;}-1}-\frac{{10}^{(a-\frac{b}{T_1+c})}}{(\mathrm{\&alpha;}-1)P}---\left(1\right)$

$\mathrm{Xn}=\frac{\mathrm{P\&alpha;}}{{10}^{(a-\frac{T_n+c}{b})}(\mathrm{\&alpha;}-1)}-\frac{1}{\mathrm{\&alpha;}-1}---\left(2\right)$

Wherein, Y ₁Be the component of nitrogen in the nitrogen product in the air separation column, Xn is the component of nitrogen in the liquid oxygen product, and P is last tower pressure, T ₁, T _nBe respectively column overhead, column bottom temperature, α is a relative volatility, and a, b, c are the Peter Antonie constant;

4) controller reads Y from data storage device ₁With the value of Xn as input, the s operation control variable R, the value of L, its formula is (3), (4), (5)::

$\mathrm{\&Delta;R}\left(k\right)=\frac{1}{2\sqrt{a_1}}\&times;\frac{(1-T\&times;K_1)\&times;{\stackrel{.}{Y}}_1\left(k\right)+K_2\&times;T\&times;(Y_{1\mathrm{set}}-Y_1\left(k\right))}{\sqrt{Y_1\left(k\right)+T(K_1(Y_{1\mathrm{set}}-Y_1\left(k\right))+K_2\mathrm{\&Sigma;}(Y_{1\mathrm{set}}-Y_1\left(k\right))}}---\left(3\right)$

$\mathrm{\&Delta;L}\left(k\right)=\frac{1}{2\sqrt{a_2}}\&times;\frac{(1-T\&times;K_1)\&times;{\stackrel{.}{X}}_n\left(k\right)+K_2\&times;T\&times;(X_{\mathrm{nset}}-X_n\left(k\right))}{\sqrt{X_n\left(k\right)+T(K_1(X_{\mathrm{nset}}-X_n\left(k\right))+K_2\mathrm{\&Sigma;}(X_{\mathrm{nset}}-X_n\left(k\right))}}---\left(4\right)$

R(k)＝R(k-1)+ΔR(k)

(5)

L(k)＝L(k-1)+ΔL(k)

Wherein, T is the sampling period, K ₁, K ₂Be controller adjustable parameter, Y _1set, X _NsetBe respectively Y ₁Setting value with Xn;

5) with R (k), the data-signal of L (k) returns to air separation column.

6, control method as claimed in claim 5 is characterized in that: described control method also comprises:

6), fall into a trap in described step 4) and to have calculated control variables R (k), the value of L (k), and with it and detect the Y that obtains ₁, Xn value on the man-machine interface of controller, show.

7, control method as claimed in claim 6, it is characterized in that: described data storage device is the historical data base of DCS system, described DCS system is made of data-interface, control station and historical data base, in described step 6), data are passed to the DCS system, and at the control station procedure for displaying state of DCS.

Images