CN101458521B - Misoperation prevention method by DCS - Google Patents

Misoperation prevention method by DCS Download PDF

Info

Publication number
CN101458521B
CN101458521B CN2008102473682A CN200810247368A CN101458521B CN 101458521 B CN101458521 B CN 101458521B CN 2008102473682 A CN2008102473682 A CN 2008102473682A CN 200810247368 A CN200810247368 A CN 200810247368A CN 101458521 B CN101458521 B CN 101458521B
Authority
CN
China
Prior art keywords
value
amplitude
given
limit
time
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN2008102473682A
Other languages
Chinese (zh)
Other versions
CN101458521A (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.)
China National Offshore Oil Corp CNOOC
China BlueChemical Ltd
CNOOC Tianye Chemical Ltd
Original Assignee
China National Offshore Oil Corp CNOOC
China BlueChemical Ltd
CNOOC Tianye Chemical 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 China National Offshore Oil Corp CNOOC, China BlueChemical Ltd, CNOOC Tianye Chemical Ltd filed Critical China National Offshore Oil Corp CNOOC
Priority to CN2008102473682A priority Critical patent/CN101458521B/en
Publication of CN101458521A publication Critical patent/CN101458521A/en
Application granted granted Critical
Publication of CN101458521B publication Critical patent/CN101458521B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/02Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]

Landscapes

  • Feedback Control In General (AREA)

Abstract

The invention relates to a method of DCS preventing misoperation which includes steps as follows: firstly, determining a computing formula to boundary value SHk, SHk and limiting value SH, SL by using regression computing method and temporal analysis: recording initialize upper limit and down limit in DCS parameter as an initialize value before putting limit switch; secondly, judging amend given VS value whether in limiting range; thirdly, judging counter whether equal to pre-set counting value, namely judging whether delay arrival time; fourthly, judging boundary value SHk and SLk and amending SH and SL. The method can hold up most of misoperation by operator by using advanced control, intelligent operation judgement and misoperation holding up alarm prompt system, ensures production safety, and can enhances enterprise's economic efficiency by using long period production. The method provided by the invention can be used in various petrochemical operation produce.

Description

A kind of DCS prevents the method for maloperation
Technical field
The present invention relates to a kind of DCS (Distributed Control System (DCS)) method of operating, particularly about a kind of method that prevents the DCS maloperation of set-point (SV value) in proportion integration differentiation control (PID) loop of being used for of in advanced person control (APC) intellectualized technology field, revising.
Background technology
In Chemical Manufacture, widely-used DCS carries out the automatic control of production technology.In the DCS system, when the operative employee controls the SV value modification in (PID) loop to a proportion integration differentiation, because the size of SV value has determined the pid loop output controlling value size of PV value (PV is the pid loop measured value) amplitude amplitude just, do not exceed when the SV value under the situation of upper limit SH, lower limit SL scope, the DCS system can not do any judgement and carry out; When the SV value exceeds upper limit SH or lower limit SL, system will eject the alarm dialog box, and operation can be selected to continue to carry out or interrupt carrying out, and upper limit SH or lower limit SL are traditionally arranged to be technological limits parameter (as shown in Figure 6).Prior art is that the control method that high and low limit (being upper and lower limit) is carried out the edge amplitude limit is exported in the employing that skies company uses, so that do not exceed edge amplitude limit scope during operation, to protect handling safety.But in the high and low limited field of output, once revise given amplitude then without limits, so under many circumstances, because the operative employee's is careless, once revises and give the fixed step size amplitude excessive, under the effect of PID, can bring bigger overshoot to system, thereby security of system is brought very big danger, gently then cause system shutdown, heavy then can cause security incident.
Summary of the invention
At the problems referred to above, the purpose of this invention is to provide and a kind ofly can carry out intelligent decision given SV value step-length amplitude once, make the reliable DCS of security of system prevent the method for maloperation.
For achieving the above object, the present invention takes following technical scheme: a kind of DCS prevents the method for maloperation, and it may further comprise the steps: step 1, utilize regressing calculation method and temporal analysis to determine boundary value SH k, SL kComputing formula with amplitude limit value SH, SL: before dropping into the amplitude limit switch, initial setting upper limit magnitude and the Lower Limit Amplitude write down in the DCS parameter are initial set value, and according to technology with reference to the calculation of parameter value of revising the boundary: the boundary value SH of upper bound scope k=SH '-K pA; The boundary value SL of lower bound scope k=SL '+K pA, wherein A be technology with reference to parameter, K pIt is correction factor; At the amplitude limit switch is amplitude limit switch when dropping into for the first time: revising given maximum limit amplitude SH is SV+K/b 1-a 1Revising given irreducible minimum amplitude SL is SV-K/b 1-a 1 Step 2, judge to revise given SV value whether in limited range: when the described amplitude limit switch of input, if the given amplitude of given described SV value goes up the amplitude range value or during less than described amplitude range value down greater than described, the DCS system sends warning automatically and does not carry out, described SV value is that given preceding initial value is constant, original state is returned in the time-delay of not enabling counting device; If the given amplitude of given described SV value is less than described upward amplitude range value or greater than described amplitude range value down, in the time of in limited amplitude range, described SV value is this given SV value, start described counter time-delay; After described amplitude limit switch disconnects, then described upper and lower scope amplitude limit value is returned to described initial value; Step 3, judge that whether counter equals to set count value, promptly judges whether delay arrival time: when described Counter Value equals to set count value, just revise delay time then, enter boundary value and judge; When described Counter Value was not equal to the setting count value, just described correction time-delay was not then then returned original state, re-executes order, continues to wait for, equals to set count value up to described Counter Value; Step 4, judgement boundary value SH kAnd SL kAnd revise SH and SL: when described Counter Value equals to set count value, just revise delay time then, judge that whether upper boundary values is greater than SV+K/b earlier 1-a 1If described boundary value is greater than SV+K/b 1-a 1, then revise given amplitude maximal value and equal SV+K/b 1-a 1, enter lower border value and judge; If described upper boundary values is less than SV+K/b 1-a 1, the given amplitude maximal value of then described correction equals described upper boundary values, enters described lower border value and judges; If described lower border value is less than SV-K/b 1-a 1, then revise given amplitude minimum value and equal SV-K/b 1-a 1, described counter O reset; If described lower border value is greater than SV-K/b 1-a 1, the given amplitude minimum value of then described correction equals described lower border value, and original state is returned in described counter O reset.
Described amplitude limit switch is a soft switch that utilizes software to realize.
Described temporal analysis comprises first-order system and second-order system.
The described step size increments and the described maximum overshoot mathematical equation that utilize regressing calculation to obtain are: h Max=b 1Δ SV+a 1
H wherein MaxMaximum overshoot, Δ SV are step size increments, b 1And a 1Be constant; Utilize regressing calculation to obtain described step size increments and the mathematical equation of described adjusting time is:
T S=b 2·ΔSV+a 2
T wherein SBe adjusting time, b 2And a 2Be constant.
Under the prerequisite of described maximum overshoot smaller or equal to the maximum variable quantity of circuit controls value PV of described pid loop, the given described amplitude maximal value or the computing method of given described amplitude minimum value and the time interval and boundary value are: given amplitude maximal value SH=SV+K/b 1-a 1Given amplitude minimum value SL=SV-K/b 1-a 1Revise the time interval T of step size increments Δ SV Sp=b 2(K/b 1-a 1)+a 2The boundary value SH of upper bound scope k=SH '-K pA; The boundary value SL of lower bound scope k=SL '+K pA, wherein A be technology with reference to parameter, K pIt is correction factor.
The present invention is owing to take above technical scheme, it has the following advantages: 1, the present invention is owing to carry out intelligent decision the operative employee once being revised given SV value step-length amplitude, therefore super amplitude range then alarm point out and do not carry out, not super amplitude range is then only carried out in SL~SH amplitude range, and reaches upper boundary values SH when the SV value kOr lower border value SL kDuring scope, then can not operate execution up and down automatically.2, the present invention is equivalent to disturbance of pid loop owing to the given SV value of step, therefore the given SV value of step step-length amplitude is big more as can be known, then amplitude is big more in the pid loop amplitude process, and in the temporal analysis of the calculating employing automatic control system of carrying out control SL~SH amplitude range, the adjusting time T of the dynamic perfromance of data acquistion system by experiment sWith maximum overshoot h Max, determine step given step size increments Δ SV and maximum overshoot h according to the one-variable linear regression method MaxAnd adjusting time T sFuntcional relationship, thereby calculate SL and the upper and lower bounds value of SH SH kAnd SL k, so SH k, SL kValue can be done according to process system controlled variable A and on-the-fly modify.3, the present invention is because DCS adopts an amplitude limit switch (being soft switch) control operation, make fool proof input and withdraw from and obtained effective control, and the amplitude limit switch writes down former SL, SH automatically when not dropping into initial value is SH ' and SL ', drops into and withdraws from the value that original SL, SH are recovered in the back.4, the present invention is owing to after the change of SV Value Operations, postpone to regulate time T sAfter, can be once more to the operation of making amendment of SV value, and revise the value of SL, SH automatically, make the upper and lower limit amplitude range can follow the tracks of variation.In sum, the present invention has realized that intellectualized operation is judged and maloperation interception alarm, therefore most of maloperation of operative employee is blocked, and has ensured the security of producing, and can satisfy the requirement that long period is produced, and can improve economic benefit of enterprises.
The present invention can be widely used in the DCS application system of various petrochemical compleies productions.
Description of drawings
Fig. 1 is the limited synoptic diagram of opereating specification of the present invention
Fig. 2 is a first-order system mathematical model synoptic diagram of the present invention
Fig. 3 is the unit-step response curve synoptic diagram of first-order system of the present invention
Fig. 4 is a second-order system mathematical model synoptic diagram of the present invention
Fig. 5 is the unit-step response curve synoptic diagram of second-order system of the present invention
Fig. 6 is a prior art system opereating specification synoptic diagram
Embodiment
Below in conjunction with drawings and Examples the present invention is described in detail.
In the pid loop of DCS (Distributed Control System (DCS)), because when the given SV value of step function, loop output control PV value (PV is the loop measured value) an amplitude process can occur under the effect of pid parameter, be equivalent to the disturbance in loop, when given SV value step-length amplitude is big more, the amplitude of PV value is just big more, so the operative employee controls the amplitude of PV value by once revising given SV value amplitude.In order to improve the reliability of DCS, avoid maloperation, the present invention prevents the DCS maloperation by the amplitude of once revising given SV value is carried out intelligent decision.
As shown in Figure 1, the limited operation scope of given SV value comprises the boundary value SH of the upper and lower boundary scope of solid line part k, SL k, boundary value SH k, SL kCan do according to the process system controlled variable and on-the-fly modify, for example technology is with reference to the parameter A value.Apart from boundary value SH k, SL kUpper and lower edge dotted portion be respectively the initial value SH ' of system and the SL ' of upper and lower amplitude range value, apart from SH kAnd SL kBetween upper and lower ladder solid line be respectively amplitude range value SH and following amplitude range value SL, the ladder dotted line between amplitude range value SH and SL is pid loop modification set-point SV.What the present invention will realize is, in the DCS system, and when given SV value surpasses amplitude range SL~SH, auto-alarm-signal keying device prompt facility then.Specifically be divided into, when the SV value surpasses amplitude range SL~SH, then alarm prompting and not carrying out; When the SV value is no more than amplitude range SL~SH, then only in amplitude range SL~SH, carry out; When given SV value reaches upper boundary values SH kOr lower border value SL kDuring scope, then can not operate execution up or down automatically.
Because the amplitude of PV value is big or small relevant with the dynamic perfromance of SV value and pid loop, the state that finally can tend towards stability of the amplitude process in pid loop, hence one can see that, and the dynamic analysis of pid loop is relevant with the dynamic analysis of automatic control system.The dynamic analysis of automatic control system has temporal analysis, root-locus technique and frequency domain analysis, these three kinds of methods commonly used are analyzed the performance of linear control system, but diverse ways has the different characteristics and the scope of application, comparatively speaking, temporal analysis is a kind of method of directly in time domain system being analyzed, has intuitive and accurate advantage, and can provide the full detail of time response, be applicable to experimental technique, very easily the characterisitic parameter of controlled system.Therefore the present invention adopts the temporal analysis of automatic control system, the dynamic Characteristic Data of data acquistion system by experiment, and step response functions is the PV value, even h (t)=PV.Specify the analytic process of temporal analysis below:
One, typical original state analysis: the original state of stipulating system is a zero condition, that is:
c ( 0 - ) = c · ( 0 - ) = c · · ( 0 - ) = · · · = 0 .
Two, typical input signal: typical input signal is actual input approximate and abstract of numerous complicated, is meant the input signal form that often runs into according to system, more Utopian basic input functions in addition on mathematical description.The selection of typical case's input signal should make mathematical operation simple, is convenient to experimental verification again, and typical input signal commonly used has unit-step function:
1 ( t ) = 1 t > 0 0 t ≤ 0
Three, typical time response: the typical case is that original state is zero the output of system under typical input signal effect time response, it is made up of dynamic process and steady-state process two parts: dynamic process claims transient process or transient process again, the system that is meant is under typical input signal effect, and system's output is arrived the response process of end-state by original state; Steady-state process is meant system under the effect of typical input signal, when the time, t was tending towards infinity, and the form of expression of system's output quantity.
Four, the analysis of dynamic property and steady-state behaviour: the performance index of control system under typical input signal effect, form by dynamic property and steady-state behaviour two parts.
Dynamic property is meant describes stable system under the unit-step function effect, and dynamic process is the index of the changing condition of t in time, and the dynamic performance index of unit-step response has usually: t during delay d, refer to that response curve reaches half needed time of its final value for the first time; Rise time t r, refer to that response curve rises to the 90% needed time of final value from final value 10%, for the system that vibration is arranged, also may be defined as response and rise to the needed time of final value from zero first time, the rise time is a kind of tolerance of system response time; Time to peak t p, refer to that response surpasses final value and reaches the needed time of first peak value; Adjusting time t s, refer to that response reaches and remains on the needed time in final value ± 5% (or ± 2%); Maximum overshoot h Max, refer to the maximum deviation amount h (t that responds p) poor with final value h (∞).
Steady-state error is a kind of performance index of descriptive system steady-state behaviour.Measure or calculate under step function, ramp function and acceleration function effect, steady-state error is a kind of tolerance of system's control accuracy or disturbance rejection ability.When the time was tending towards infinity, the output quantity of system was not equal to definite function of input quantity or input quantity, and then there is steady-state error in system.
The present invention adopts first-order system and second-order system as input signal in the pid loop of DCS, and h (t)=PV can release the time-domain analysis of first-order system as follows by above-mentioned temporal analysis:
The mathematical model of first-order system: as shown in Figure 2, the open-loop transfer function of system is G ( s ) = 1 Ts = K s , The closed loop transfer function, of system is Φ ( s ) = K s + K = 1 Ts + 1 .
The response of first-order system: the present invention utilizes unit-step response, when input signal is unit step signal, i.e. and r (t)=1 (t), its Lars is transformed to R ( s ) = 1 s , Then have:
C ( s ) = Φ ( s ) R ( s ) = 1 s ( Ts + 1 ) = 1 s - 1 s + 1 / T
As shown in Figure 3, following formula is got the Lars inversion, can obtain the unit-step response of first-order system, its function is:
h(t)=L -1[C(s)]=1-e -t/T
According to the definition of dynamic performance index, can obtain dynamic performance index t time delay of first-order system d=0.69T, rise time t r=2.20T regulates time t s=3T (± 5%).For the unit-step response of first-order system, steady-state error e ss = lim t → ∞ e ( t ) = lim t → ∞ [ r ( t ) - c ( t ) ] = 0 , When being first-order system tracking step input signal, no steady-state error.
When in pid loop, adopting second-order system, can release time-domain analysis to second-order system by above-mentioned temporal analysis as input signal.
The mathematical model of typical case's second-order system: as shown in Figure 4, the open-loop transfer function of system is G ( s ) = ω n 2 s ( s + 2 ξ ω n ) = K s ( Ts + 1 ) , The closed loop transfer function, of system is: Φ ( s ) = ω n 2 s 2 + 2 ξ ω n s + ω n 2 , Wherein ξ is a damping ratio.
The unit-step response of second-order system: second order unit-step response function h (t) comprises steady-state component and dynamic component, and wherein steady-state component is 1, and dynamic component presents oscillatory extinction characteristic (as shown in Figure 5).The steady-state error e of second-order system Ss=0, when promptly typical second-order system was followed the tracks of step input signal, no steady-state error, system were astatic system.According to the definition of dynamic performance index, can obtain following two conclusions:
Conclusion one: damping ratio ξ is big more, and the overshoot of system is more little, and response steadily; Damping ratio ξ is more little, and the overshoot of system is big more, and the stationarity of response is poor more; When ξ=0, the response of system is h (t)=1-cos ω nT, h (t) is ω for frequency nContinuous oscillation, system can't enter balance duty, cisco unity malfunction; As ξ one timing, ω nBig more, the oscillation frequency ω of system dBig more, the stationarity of response is relatively poor.Therefore it is big to work as damping ratio ξ, frequencies omega nHour, the stationarity of system responses is good.
Conclusion two: regulate time t sComputing formula be approximate expression, in fact, when damping ratio ξ hour, the system responses speed of convergence was slow, regulated time t sLong; When damping ratio ξ was excessive, system responses was blunt, regulated time t sAlso longer.Therefore damping ratio ξ should get suitable numerical value, and the typical second-order system when damping ratio ξ=0.707 is called best second-order system, this moment overshoot h MaxBe 4.3%, regulate time t sBe 3/ ω n
Because given SV value is given by step function, in the pid loop of DCS, adopt first-order system and second-order system as input signal, can be similar to the step given step size increments Δ SV and the maximum overshoot h that think near given SV value MaxAnd adjusting time T SThere is linear relationship.Can determine step given step size increments Δ SV and maximum overshoot h by the one-variable linear regression method MaxAnd adjusting time T SFuntcional relationship, and determine computing method and the time interval T of given amplitude maximal value SH or given amplitude minimum value SL SpWith boundary value SH k, SL kExperimental method realizes that obtaining the regressing calculation data is in output loop, constantly to the given different range value of SV value, measures each maximum overshoot h MaxWith the adjusting time T S, draw some groups of data and carry out regressing calculation; In actual production process, realize obtaining the regressing calculation data and be data, analytic record is carried out in operative employee's the given operation of SV value, draw some groups of data and carry out regressing calculation by host computer opc server interface interchange DCS.Utilize the regressing calculation method can on-line measurement, and can remeasure calculating, make amendment calculating regression formula every a period of time.Utilize regressing calculation to calculate the given step size increments Δ SV that jumps and specifically comprise following step:
1) can determine to regulate time T according to the one-variable linear regression equation SWith the pass of step input signal step size increments Δ SV be:
If T S=f (Δ SV) regulates time T SFunction be the function of step input signal step size increments Δ SV, some groups of n some experimental data (Δ SV are arranged 1, T S1), (Δ SV 2, T S2) ..., (Δ SV n, T Sn).
Order T ^ S = a + bΔSV Be formula (1) that wherein a is a straight line intercept constant, b is the straight slope constant, regulates the estimated value of time
Figure G2008102473682D00071
With step size increments Δ SV be linear relationship, can determine the value of constant a and b according to least square method:
Order Q ( a , b ) = Σ i = 1 n [ T Si - ( a + b ( ΔSV ) i ) ] 2 , Then (a b) should satisfy at the limit place Q ∂ Q ∂ a = 0 ∂ Q ∂ b = 0 , That is:
Σ i = 1 n [ T Si - ( a + b ( ΔSV ) i ) ] = 0 Σ i = 1 n [ T Si - ( a + b ( ΔSV ) i ) ] ( ΔSV ) i = 0
Can get the system of linear equations of a, b by following formula, be designated as formula (2):
Σ i = 1 n T Si - na - b Σ i = 1 n ( ΔSV ) i = 0 Σ i = 1 n ( ΔSV ) i T Si - a Σ i = 1 n ( ΔSV ) i - b Σ i = 1 n ( ΔSV ) i 2 = 0
By above-mentioned normal equations can declare know Q (a b) has minimum point, separates normal equations and can obtain a and b one group and separate, be designated as (
Figure G2008102473682D00076
), can demonstrate,prove (
Figure G2008102473682D00077
) be Q (a, minimum point b), then
Figure G2008102473682D00078
With Formula be respectively formula (3) and formula (4).
Formula (3) is: a ^ = T S ‾ - b ^ ΔSV ‾ .
Formula (4) is: b ^ = Σ i = 1 n ( ΔSV ) i T Si - n ΔSVT S ‾ Σ i = 1 n ( ΔSV ) i 2 - n ( ΔSV ) 2 ‾ = Σ i = 1 n ( ( ΔSV ) i - ΔSV ‾ ) ( T Si - T S ‾ ) Σ i = 1 n ( ( ΔSV ) i - ΔSV ‾ ) 2 .
In the following formula ΔSV ‾ = Σ i = 1 n ( ΔSV ) i , T S ‾ = Σ i = 1 n T Si .
Formula (3) and (4) formula are determined
Figure G2008102473682D000714
Can obtain formula (5) in the substitution formula (1) is:
T ^ S = a ^ + b ^ ΔSV
Formula (5) is exactly step size increments Δ SV and adjusting time T SBetween the correlation experience formula, be called regression beeline equation.
2) to regulating time T SCarry out correlation test, have:
Residual sum of squares (RSS): Q = Σ i = 1 n ( T Si - T ^ Si ) 2 ;
Regression sum of square: U = Σ i = 1 n ( T ^ Si - T S ‾ ) 2 ;
Related coefficient: R = Σ i = 1 n ( T ^ Si - T s ‾ ) 2 Σ i = 1 n ( T Si - T S ‾ ) 2 ;
In the above-mentioned formula, when residual sum of squares (RSS) Q was big, then meter sensitivity increased, otherwise meter sensitivity diminishes.When coefficient R is big, then there is linear dependence; When R hour, then do not have linear dependence.
3) can determine maximum overshoot h according to the one-variable linear regression equation MaxWith the pass of step input signal step size increments Δ SV be:
If h Max=f (Δ SV), i.e. maximum overshoot h MaxFunction be the function of step input signal step size increments Δ SV, have n the point some groups of experimental data (Δ SV 1, h Max1), (Δ SV 2, h Max2) ..., (Δ SV n, h Maxn).
Order h ^ max = a + bΔSV Be formula (6) that wherein a is a straight line intercept constant, b is the straight slope constant, the estimated value of maximum overshoot
Figure G2008102473682D00082
With step size increments Δ SV be linear relationship, can determine the value of constant a and b according to least square method, solution procedure and above-mentioned definite adjusting time T SIdentical with the value method of finding the solution constant a and b in the step input signal increment Delta SV relation, therefore can obtain formula (7) and be:
h ^ max = a ^ + b ^ ΔSV
Formula (7) is exactly step size increments Δ SV and maximum overshoot h MaxBetween the correlation experience formula, be called regression beeline equation.
4) to maximum overshoot h MaxCarry out correlation test, have:
Residual sum of squares (RSS): Q = Σ i = 1 n ( h max i - h ^ max i ) 2
Regression sum of square: U = Σ i = 1 n ( h ^ max i - h max ‾ ) 2
Related coefficient: R = Σ i = 1 n ( h ^ max i - h max ‾ ) 2 Σ i = 1 n ( h max i - h max i ‾ ) 2
In the above-mentioned formula, when residual sum of squares (RSS) Q was big, then meter sensitivity increased, otherwise meter sensitivity diminishes.When coefficient R is big, then there is linear dependence; When R hour, then do not have linear dependence.
By above-mentioned 1)~4) the step analytical calculation, according to experimental data,, can obtain maximum overshoot h by regression beeline equation (7) MaxMathematical equation with step size increments Δ SV is designated as: h Max=b 1Δ SV+a 1By regression beeline equation (5), can adjusted time T SMathematical equation with step size increments Δ SV is designated as: T S=b 2Δ SV+a 2According to the requirement of process system, the maximum overshoot h of pid loop MaxShould≤K, K is the maximum variable quantity of the circuit controls value PV of technological requirement, and K value determined by the technology personnel, so can determine computing method and the time interval T of given amplitude limit maximal value SH or given amplitude limit minimum value SL SpWith boundary value SH k, SL k, then:
Given amplitude limit maximal value SH=SV+K/b 1-a 1
Given amplitude limit minimum value SL=SV-K/b 1-a 1
Revise the time interval T of step size increments Δ SV Sp=b 2(K/b 1-a 1)+a 2
The boundary value SH of upper bound scope k=SH '-K pA;
The boundary value SL of lower bound scope k=SL '+K pA,
Wherein A be technology with reference to parameter, K pBe correction factor, K pCan set according to actual conditions and experience.Because first-order system does not have overshoot h Max, can set the scope of initial set value SH and SL according to actual conditions, regulate time T SCan draw through the reality test according to actual conditions.
Comprehensive above explanation can obtain the following concrete operations step that amplitude that the present invention once revises given SV value is carried out intelligent decision:
Step 1, utilize the temporal analysis of regressing calculation method and automatic control system to determine to boundary value SH k, SL kComputing formula with amplitude limit value SH, SL: amplitude limit switch (for a soft switch that utilizes software to realize) is to drop into and the operation of withdrawing from anti-misoperation software in order to control.Before dropping into the amplitude limit switch, the initial setting upper limit magnitude SH, the Lower Limit Amplitude SL that write down in the DCS parameter are SH ' and SL ', and according to technology with reference to the calculation of parameter value of revising the boundary SH kAnd SL kBe respectively:
The boundary value SH of upper bound scope k=SH '-K pA;
The boundary value SL of lower bound scope k=SL '+K pA,
Wherein A be technology with reference to parameter, K pBe correction factor, and K pCan set according to actual conditions and experience.
At the amplitude limit switch is amplitude limit switch when dropping into for the first time:
Revising given maximum limit amplitude SH is SV+K/b 1-a 1
Revising given irreducible minimum amplitude SL is SV-K/b 1-a 1
Step 2, judge to revise set-point (SV value) whether in limited range: when input amplitude limit switch is the amplitude limit switch, if the given amplitude of given SV value is during greater than last amplitude range value SH or less than following amplitude range value SL, the DCS system sends warning automatically and does not carry out, the SV value is that given preceding initial value is constant, original state is returned in the time-delay of not enabling counting device; If the given amplitude of given SV value is less than last amplitude range value SH or greater than following amplitude range value SL, in the time of in limited amplitude range, the SV value is this given SV value, the time-delay of enabling counting device.
When the amplitude limit switch is after the amplitude limit switch disconnects, then upper and lower scope amplitude limit value SH and SL to be returned to initial value SH ' and SL '.
Step 3, judge that whether counter equals to set count value, promptly judges whether delay arrival time: when Counter Value equals to set count value, just revise delay time then, enter boundary value and judge; When Counter Value is not equal to the setting count value, just revise time-delay not then, then return original state, re-execute order, continue to wait for, equal to set count value up to Counter Value.
Step 4, judgement boundary value SH kAnd SL kAnd revise SH and SL: when Counter Value equals to set count value, just revise delay time then, judge upper boundary values SH earlier kWhether greater than SV+K/b 1-a 1If, upper boundary values SH kGreater than SV+K/b 1-a 1, then revise given amplitude maximal value SH=SV+K/b 1-a 1, enter lower border value SL kJudge; If upper boundary values SH kLess than SV+K/b 1-a 1, then revise given amplitude maximal value SH=SH k, enter lower border value SL kJudge.If lower border value SL kLess than SV-K/b 1-a 1, then revise given amplitude minimum value SL=SV-K/b 1-a 1, described counter O reset; If lower border value SL kGreater than SV-K/b 1-a 1, then revise given amplitude minimum value SL=SL k, original state is returned in described counter O reset.
By the various embodiments described above as can be seen, the present invention comes control operation by adopting an amplitude limit switch, automatically writing down former upper and lower amplitude range value SH and SL when the amplitude limit switch does not drop into is initial set value SH ' and SL ', amplitude limit switch input is withdrawn from the back and is recovered original amplitude range value SL, SH, the input of amplitude limit switch of the present invention and the setting of withdrawing from effective control data.And owing to after the change of SV Value Operations, postpone to regulate time T sAfter, can revise the value of SL, SH automatically, and the amplitude limit scope can be followed the tracks of variation.

Claims (6)

1. a DCS prevents the method for maloperation, and it may further comprise the steps:
Step 1, utilize regressing calculation method and temporal analysis to determine to upper boundary values SH k, lower border value SL kComputing formula with upper limit magnitude SH, Lower Limit Amplitude SL: before dropping into the amplitude limit switch, the initial setting upper limit magnitude SH and the Lower Limit Amplitude SL that write down in the DCS parameter are initial set value, and according to technology with reference to the calculation of parameter value of revising the boundary:
Upper boundary values SH k=SH '-K pA;
Lower border value SL k=SL '+K pA,
Wherein A be technology with reference to parameter, K pBe correction factor, SH ' is system's initial value of upper limit magnitude SH, and SL ' is system's initial value of Lower Limit Amplitude SL;
At the amplitude limit switch is amplitude limit switch when dropping into for the first time:
The maximum limit amplitude of revising given upper limit magnitude SH is SV+K/b 1-a 1
The irreducible minimum amplitude of revising given Lower Limit Amplitude SL is SV-K/b 1-a 1
Wherein SV is that step is given, and K is the maximum variable quantity of the circuit controls value PV of technological requirement, b 1, a 1Be constant;
Step 2, judge to revise the given SV value of step whether in limited range: when the described amplitude limit switch of input, if the given amplitude of the given SV value of given described step is during greater than the maximum limit amplitude of described given upper limit magnitude SH or less than the irreducible minimum amplitude of described given Lower Limit Amplitude SL, the DCS system sends warning automatically and does not carry out, the given SV value of described step is that given preceding initial value is constant, original state is returned in the time-delay of not enabling counting device; If the given amplitude of the given SV value of given described step is less than the maximum limit amplitude of described given upper limit magnitude SH or greater than the irreducible minimum amplitude of described given Lower Limit Amplitude SL, in the time of in limited amplitude range, the given SV value of described step is the given SV value of this given step, starts described counter time-delay;
After described amplitude limit switch disconnects, then the maximum limit amplitude of described given upper limit magnitude SH, the irreducible minimum amplitude of given Lower Limit Amplitude SL are returned to described initial value;
Step 3, judge that whether counter equals to set count value, promptly judges whether delay arrival time: when described Counter Value equals to set count value, just revise delay time then, enter boundary value and judge; When described Counter Value was not equal to the setting count value, just described correction time-delay was not then then returned original state, re-executes order, continues to wait for, equals to set count value up to described Counter Value;
Step 4, judgement upper boundary values SH kWith lower border value SL kAnd revise upper limit magnitude SH and Lower Limit Amplitude SL: when described Counter Value equals to set count value, just revise delay time then, judge upper boundary values SH earlier kWhether greater than SV+K/b 1-a 1If, described boundary value SH kGreater than SV+K/b 1-a 1, the maximum limit amplitude of then revising given upper limit magnitude SH equals SV+K/b 1-a 1, enter lower border value SL kJudge; If described upper boundary values SH kLess than SV+K/b 1-a 1, the maximum limit amplitude of the given upper limit magnitude SH of then described correction equals described upper boundary values SH k, enter described lower border value SL kJudge; If described lower border value SL kLess than SV-K/b 1-a 1, the irreducible minimum amplitude of then revising given Lower Limit Amplitude SL equals SV-K/b 1-a 1, described counter O reset; If described lower border value SL kGreater than SV-K/b 1-a 1, the irreducible minimum amplitude of the given Lower Limit Amplitude SL of then described correction equals described lower border value SL k, original state is returned in described counter O reset.
2. a kind of DCS as claimed in claim 1 prevents the method for maloperation, it is characterized in that: described amplitude limit switch is a soft switch that utilizes software to realize.
3. a kind of DCS as claimed in claim 1 prevents the method for maloperation, it is characterized in that: described temporal analysis comprises first-order system and second-order system.
4. a kind of DCS as claimed in claim 2 prevents the method for maloperation, it is characterized in that: described temporal analysis comprises first-order system and second-order system.
5. prevent the method for maloperation as claim 1 or 2 or 3 or 4 described a kind of DCS, it is characterized in that: the step size increments and the maximum overshoot mathematical equation that utilize regressing calculation to obtain are:
h max=b 1·ΔSV+a 1
H wherein MaxMaximum overshoot, Δ SV are step size increments, b 1And a 1Be constant; Utilize regressing calculation to obtain step size increments to be with the mathematical equation of regulating the time:
T S=b 2·ΔSV+a 2
T wherein SBe adjusting time, b 2And a 2Be constant.
6. a kind of DCS as claimed in claim 5 prevents the method for maloperation, it is characterized in that: under the prerequisite of described maximum overshoot smaller or equal to the maximum variable quantity of circuit controls value PV of pid loop, the computing method of the maximum limit amplitude of given described given upper limit magnitude SH, the irreducible minimum amplitude of given Lower Limit Amplitude SL and the time interval and boundary value are:
Maximum limit amplitude=SV+K/b of given upper limit magnitude SH 1-a 1
Irreducible minimum amplitude=SV-K/b of given Lower Limit Amplitude SL 1-a 1
Revise the time interval T of step size increments Δ SV Sp=b 2(K/b 1-a 1)+a 2
The boundary value SH of upper bound scope k=SH '-K pA;
The boundary value SL of lower bound scope k=SL '+K pA,
Wherein A be technology with reference to parameter, K pIt is correction factor.
CN2008102473682A 2008-12-30 2008-12-30 Misoperation prevention method by DCS Expired - Fee Related CN101458521B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2008102473682A CN101458521B (en) 2008-12-30 2008-12-30 Misoperation prevention method by DCS

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2008102473682A CN101458521B (en) 2008-12-30 2008-12-30 Misoperation prevention method by DCS

Publications (2)

Publication Number Publication Date
CN101458521A CN101458521A (en) 2009-06-17
CN101458521B true CN101458521B (en) 2010-12-29

Family

ID=40769441

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008102473682A Expired - Fee Related CN101458521B (en) 2008-12-30 2008-12-30 Misoperation prevention method by DCS

Country Status (1)

Country Link
CN (1) CN101458521B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102930219A (en) * 2012-10-31 2013-02-13 龚莉雅 Processing method and device for data setting
CN105388782B (en) * 2015-12-10 2018-06-08 天津电气科学研究院有限公司 A kind of algorithm for the sigmoid curve generator for being used to give link
CN106936801B (en) * 2015-12-31 2020-09-18 北京网御星云信息技术有限公司 Method and device for realizing message filtering
US10248097B2 (en) * 2016-05-09 2019-04-02 Yokogawa Electric Corporation Method and apparatus for optimizing process control systems

Also Published As

Publication number Publication date
CN101458521A (en) 2009-06-17

Similar Documents

Publication Publication Date Title
US20190187631A1 (en) Adaptive pid controller tuning via deep reinforcement learning
CN101458521B (en) Misoperation prevention method by DCS
US6631299B1 (en) System and method for self-tuning feedback control of a system
RU2501063C2 (en) Residue analysis apparatus and method for detecting system errors in aircraft system behaviour
CN111812967B (en) PID control parameter setting method based on stability margin and dynamic response index
CN108845491A (en) The wisdom PI composite control method of Correction for Large Dead Time System
CN104345640A (en) Motor servo system progressive tracking control method and system during input limitation
US7925362B2 (en) Autotuning method using integral of relay feedback response for extracting process information
Yerashunas et al. Control of lateral motion in moving webs
EP3575893A1 (en) Control parameter computation method and control parameter computation device
di Capaci et al. A revised technique of stiction compensation for control valves
Seki et al. Retuning oscillatory PID control loops based on plant operation data
Bertino et al. Experimental and analytical prescribed-time trajectory tracking control of a 7-DOF robot manipulator
CN101446829B (en) Method for preventing DCS misoperation
Huang Steady state and transient state identification in an industrial process
Papadoulis et al. A cautious self‐tuning controller for chemical processes
Lee et al. A controller design method for constructing a robust track-following system
JP5010717B2 (en) PID controller
Shi et al. Fmea-based control mechanism for embedded control software
Taylor et al. Loop-shaping controller design for nonlinear systems using the contoured robust controller bode (crcbode) plot
Visioli Assessment of tuning of PI controllers for self-regulating processes
Welz et al. Validation of a solution model for the optimization of a binary batch distillation column
Šarafín Input Shaping
CN117913853A (en) New energy power generation system frequency safety boundary prediction method based on LSTM
Campo et al. Decentralized Control System Design for a Heavy Oil Fractionator–

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20101229

Termination date: 20151230

EXPY Termination of patent right or utility model