CN108490989A - A kind of flow control system and efficient control method - Google Patents
A kind of flow control system and efficient control method Download PDFInfo
- Publication number
- CN108490989A CN108490989A CN201810076141.XA CN201810076141A CN108490989A CN 108490989 A CN108490989 A CN 108490989A CN 201810076141 A CN201810076141 A CN 201810076141A CN 108490989 A CN108490989 A CN 108490989A
- Authority
- CN
- China
- Prior art keywords
- flow
- control system
- value
- pump
- flow control
- 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.)
- Withdrawn
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D7/00—Control of flow
- G05D7/06—Control of flow characterised by the use of electric means
- G05D7/0617—Control of flow characterised by the use of electric means specially adapted for fluid materials
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D7/00—Control of flow
- G05D7/06—Control of flow characterised by the use of electric means
- G05D7/0617—Control of flow characterised by the use of electric means specially adapted for fluid materials
- G05D7/0629—Control of flow characterised by the use of electric means specially adapted for fluid materials characterised by the type of regulator means
- G05D7/0676—Control of flow characterised by the use of electric means specially adapted for fluid materials characterised by the type of regulator means by action on flow sources
Abstract
The present invention provides a kind of flow control system and efficient control method.Its slope formula is obtained first.By carrying out fitting a straight line to M point under Relative steady-state state, and the slope K of the straight line is obtained, to obtain pressure P；Secondly, the intersection point of the heat-capacity curve and the head curve and corresponding high efficiency range boundary by every pump of Q H Plane-points (Q, P) is obtained；Finally obtain Effec-tive Function performance parameter, and then the pump work for selecting performance parameter optimal, it is ensured that flow control system is in Effec-tive Function.The Effec-tive Function control of flow control system can be realized without press detection sensor and auxiliary circuit by the present invention, eliminate the installation and debugging required time and cost of pressure sensor and additional processing circuitry, so that system structure is simpler, system cost is lower, operational efficiency higher, system lifetim and reliability are improved, reliable guarantee is provided for the safe efficient operation of flow control system.
Description
Technical field
The invention belongs to process control fields, and in particular to a kind of flow control system and efficient control method realize stream
The Effec-tive Function of amount control system.
Background technology
Flow-rate adjustment control has a wide range of applications in fields such as chemical industry, food, medicine, water supply.Early stage flow control master
To be opened by the aperture of adjusting control valve and output flow be adjusted, but there are high energy consumption, adjustable range is little the deficiencies of.
The current main scheme using governor impeller realizes that the adjusting of output flow, principle are mainly flowed by detecting output flow with setting
Deviation between amount, and feedback compensation control algorithm is carried out to the deviation, and then the output frequency of frequency converter is adjusted, change pump
Rotating speed realizes the stability contorting of output flow.However, delivered in control theory and application periodical by Zhang Chenghui etc.《Frequency conversion
Governing feedwater pump station efficiency-optimized control strategy》One text is it is found that pump there are one by heat-capacity curve, similar operating condition parabola
The Effec-tive Function section of composition.Pump operation can realize Effec-tive Function in this section；Otherwise, the way and service life will significantly
It reduces.In process control fields such as chemical industry, food, medicine, pump is widely used so that efficiency all can even only improving 1%
To energy saving and environmentally friendly bring huge interests.
However, flow control system operating status is not invariable, output flow and lift have in time
There is time variation, it cannot be guaranteed that pump operates between high efficient area always.On the one hand, in chemical industry, food, medicine, the fields such as supply water, by
It is run in the duct for a long time in the liquid of conveying, may have dirt deposition, entire pipeline effective sectional area is caused to become smaller,
Pipe resistance characteristic is deteriorated, and in the case where setting traffic conditions, pump discharge and ductwork pressure increased dramatically, and pump operation state is caused to become
Change, high efficiency range may be deviateed；On the other hand, since flow control system liquid matter may be transported in different containers, cause
Its lift changes, and pump operation state is caused to change, and may deviate high efficiency range；In another aspect, flow control system
Output flow be limited by the working condition of real system or the difference of working stage and different, cause pump operation state to occur
Variation, may deviate high efficiency range；In non-efficient section, operation can lead to the efficiency of frequency conversion flow control system to pump for a long time
It reduces, even results in overload/low-frequency operation of frequency converter and pump, increase the failure risk of frequency conversion flow control system.In order to
Ensure the Effec-tive Function of flow control system, must just carry out effectively optimizing control to pump.
Invention content
It is an object of the present invention to overcome the above deficiencies, proposes a kind of flow control simple in structure, applicability is good
System and efficient control method.
The present invention provides a kind of flow control system and efficient control method, and its step are as follows：
1) pressure value P of the flow control system in stable state and t ∈ [0, T are establishedd] changes in flow rate amount Δ q1(t) pass
It is formula：Wherein；P is ductwork pressure value, and F is frequency converter output frequency, and Q is liquid in-out
Body flow, T are environment temperature, TbFor air pressure tank rated temperature, VbFor air pressure tank gas chamber nominal volume, PbFor air pressure tank gas chamber volume
Constant-pressure, t are time variable, TdFor pre-defined observation interval, Δ F is frequency disturbance increment；
2) by obtaining slope in the relational expression of foundation in step 1)
3) with sampling period TsFlow control system flow value and frequency converter output frequency are sampled for interval, and obtained
It is sampling number to take flow value q (k) and output frequency f (k), wherein k；
4) t=kT is found out by pid control algorithmsThe output frequency value f (k) of moment frequency converter=f (k-1)+Kp[e(k)-e
(k-1)]+Kie(k)+Kd[e(k)-2e(k-1)+e(k-2)]；
Wherein, e (k-1), f (k-1) are respectively t=(k-1) TsThe flow error at moment and the output frequency of frequency converter；e
(k-2) it is t=(k-2) TsThe flow error at moment, Kp、KiAnd KdProportionality coefficient in respectively preset pid algorithm,
Integral coefficient and differential coefficient, flow error e (k)=Qset-q(k)；Wherein, e (i) |＜=0 i=0；QsetTo set output flow
Value；The flow value that q (k) is sampling number when being k, the output frequency value of f (k) is sampling number when being k frequency converter；f(i)|＜=0 i
=0；
5) and according to the flow value q (k) and output frequency f (k) sampled, the flow value number being made of N number of element is established
Group { q (i) } and frequency converter output frequency array { f (i) }, wherein i=k-N+1, k-N+2 ... and k }, N is to preset
Be more than 1 positive integer, q (i) |＜=0 i=0, f (i) |＜=0 i=0；
6) average value of flow value array { q (i) } is obtainedAnd judge whether flow control system is located
In metastable state；
7) when determining that flow control system is in metastable state, then judgeWhether at
Vertical, wherein θ is setting positive value, QsetTo set output flow；
8) in determinationWhen invalid, then it is considered as the pump M of work at presentjH-MaxToo
It is small, j=1,2,3, enable control switch Sj=0, control MjIt is out of service；Meanwhile enabling control switch Sj+1=1, control lift is big by one
The pump M of gradej+1Work, j+1<=3, update k=k+1；It is sampled next time, and output flow value and frequency converter is marked to export
The sampled value of frequency is q (k) and f (k), repeats above step；
9) in determinationWhen establishment, then the average value of frequency converter output frequency is obtained
10) t=0 is denoted as with blaze at this time, gives output frequency one smaller disturbance quantity Δ F, f (mTs)=F+ Δs F；
11) judge whether m ＞ M are true, wherein
12) if m ＞ M are invalid, in t=mTsMoment, sample streams magnitude q (m)；It obtainsAnd sentence
It is disconnectedWherein α is setting positive value；If so, then update m=m+1, and re-start m ＞ M judgement, if not at
It is vertical, then update k=k+1；It is sampled next time, and it is q (k) to mark output flow value and the sampled value of frequency converter output frequency
With f (k), above step is repeated；IfIt sets up, then to M point (Ts,Δq(1)),(2Ts,Δq(2)),…,
(MTs, Δ q (M)) and fitting a straight line is carried out, and obtain the slope K of the straight line；
13) pass through acquisition in step 2)Estimate pressure P；
14) have the characteristics that translate according to point r (Q, P) and governor impeller its lift characteristic, obtain the pump of passing point r (Q, P)
MjQ-H heat-capacity curvesAnd obtain pump MjQ-H heat-capacity curvesCorresponding frequency
Fj；
15) pump M is obtainedjQ-H heat-capacity curvesTo similar parabolaIntersection point ajAnd bj, j=1,2,
3；
16) it calculates and meetsThe pump M of conditionjEffec-tive Function performance parameter ψj, ψjFor nonnegative number, ψjFor
Intersection point aj、bjWith r (Q, P) and corresponding high efficiency range AjBjCjDjFunction；
17) ψ is soughti=max { ψ1,ψ2,ψ3Corresponding i, i=1,2,3, controller is by corresponding switch Si(t)=1, Su
=0, u=1,2,3 ∩ u ≠ i, and update k=k+1；It is sampled next time, and output flow value and frequency converter is marked to export frequency
The sampled value of rate is q (k) and f (k), repeats above step.
Include the following steps in step 1)：
1. establishing flow control system water pump output power equation：
Wherein：η is the efficiency of pump, i.e., the ratio between motor effective power and shaft power, s are revolutional slip, R1,R2,X1σ,
X2σ,m1,For the intrinsic parameter of pump motor；
2. carrying out small signal disturbance to the equation of step 1., obtains relational expression and be reduced to
QΔp(t)+PΔq1(t)+Δq1(t) Δ p (t)=k'(2F Δ F+ Δs F2), wherein：q1(t)=Q+ Δs q1(t), f
(t)=F+ Δs F, p (t)=P+ Δ p (t), k'=η k/ ρ,
3. in t ∈ [0, Td] obtain system small-signal model equation：
4. obtaining in t ∈ [0, Td], the volume change of air pressure tank liquid chamber；
And thus obtain t ∈ [0, Td] when liquid chamber volume
Chamber volumeAnd it is obtained according to equation for ideal gases
The air chamber pressure variable quantity of air pressure tankAnd thus obtain ductwork pressure variable quantityAnd in pa(0)=P when obtain
5. 3. 4. being obtained with step according to stepAnd it finally obtains
The efficient region ABCD is rated frequency fNHeat-capacity curve HN, low-limit frequency fminLift characteristic it is bent
Line Hmin, similar operating condition parabola li1, similar operating condition parabola li2The fan annular region surrounded.
The average value of flow value array { q (i) } is obtained in step 6)And it solvesJudge whether to meet：σq≤εq, wherein：εqTo set positive value, if satisfied, then thinking flow control
System processed is in stable state.
ψjRegion A is arrived for r (Q, P)jBjCjDjGeometric center point distance derivative, wherein j=1,2,3.
ψjWhether it is in A for r (Q, P)jBjCjDjIn region and r (Q, P) arrives AjBjCjDjThe weighting letter of the distance of boundary curve
Number, wherein j=1,2,3.
ψjWhether it is in A for r (Q, P)jBjCjDjIn region and r (Q, P) arrives aj、bjThe weighting function of distance, wherein j=1,
2,3.
The present invention has the advantages that：
One, flow control system of the present invention and efficient control method have system pressure on-line checking, without pressure
Force snesor has saved system installation and debugging required time and cost so that system structure is simpler, and system cost is lower；
Two, pressure online test method of the present invention is simple with algorithm, detection speed is fast, highly practical and reliable
The advantages that property is high；
Three, flow control system of the present invention and efficient control method can obtain output flow Q and lift (or pressure
Power) on the basis of the data such as P, the heat-capacity curve according to variable-frequency control pump has translation feature, obtains by Q-H planes
The heat-capacity curve of every pump of point (Q, P) and the intersection point of the head curve and corresponding high efficiency range boundary.It is pumped according to every
High efficiency range obtain Effec-tive Function performance parameter with position relationship between (Q, P) and intersection point, and then it is former based on efficiency optimization
Then intelligence is switched to model pump M appropriatei(i=1,2,3) works, it is ensured that system high efficiency is run, to significantly improve flow control
The working efficiency of system processed；
Four, flow control system of the present invention and efficient control method are widely portable to using variable-frequency control pump
Field of flow control has extensive versatility.This is because pressure P meets formula
The formula is by parameter, Δ q (t), Q, Δ F, F, Tb、Pb、Vb, T and t determine pressure P.Wherein：Δq(t)、Q、ΔF、F、Tb、Pb、
Vb, T and t flow when flow deviates the undulate quantity of opposite stationary value, opposite stable operation when being expressed as the F disturbance operations of frequency Δ
Nominal temperature, air pressure tank are specified when the specified operation of frequency converter output frequency, air pressure tank when amount, frequency disturbance increment, stable operation
Nominal volume, current environmental temperature and time variable when the specified operation of nominal pressure, air pressure tank when operation.Calculation of pressure formula with
Motor design parameter is unrelated, has extensive versatility.
Description of the drawings
Fig. 1 is the structure diagram of flow control system；
Fig. 2 is flow control system lift-pipe resistance characteristic figure.
Fig. 3 is governor impeller Effec-tive Function area schematic diagram.
Fig. 4 is pump operation interval diagram
Fig. 5 operating points (Q, P) and M1High efficiency range position relationship schematic diagram
Fig. 6 operating points (Q, P) and M2High efficiency range position relationship schematic diagram
Fig. 7 operating points (Q, P) and M3High efficiency range position relationship schematic diagram
Specific implementation mode
Embodiments of the present invention is further illustrated below in conjunction with the accompanying drawings：
The present invention provides a kind of flow control system and efficient control methods, mainly establish the number of flow control system
Learn model, and according to the Q-H lift characteristics of the mathematical model of foundation and pump and similar operating condition parabolic surround Effec-tive Function region to
Efficient control method is gone out.Flow control system mathematical model to establish process as follows：
Flow control system schematic diagram is as shown in Figure 1, include mainly fluid supply 1, check valve 2, M1、M2、M3For output flow
The different pump of the identical lift of range, corresponding H-Max are respectivelyWith(wherein：
)、M1Control switch S1、M2Control switch S2、M3Control switch S3, flowmeter 3, air pressure tank 4, frequency conversion and controller 5, temperature sensing
Device 6 etc..Overstriking line indicates that power cord, arrow direction indicate power transfer direction in Fig. 1.Fluid supply 1 is mainly liquid medium, can
To be water, oil, chemical solution or other liquid；2 major function of check valve is to prevent liquid from flowing backwards；Pump Mi(i=1,2,3) pass through
Liquid in fluid supply is transported to pipeline by impeller blade high speed rotation；Switch Si(i=1,2,3) M is controllediWhether operation；Flow
Meter 3 is for detecting pump discharge flow；Air pressure tank 4 is mainly the function of stablizing ductwork pressure；Frequency conversion and controller 5 mainly realize phase
The input of related parameter, the display of operating status and the operation of system control program adjust revolution speed, realize pump output flow control
System；Temperature sensor 6 is used for detecting system Current Temperatures.
Variable declaration is as follows：Q (t) is pump discharge flow；q2(t) it is air pressure tank rate of discharge；P (t) is the pressure of pipe network
Value；F (t) is frequency converter output frequency；Air pressure tank chamber volume is v1(t)；Air pressure tank air chamber pressure pa(t), air pressure tank liquid chamber body
Product is v2(t), air pressure tank total volume is Vz, air pressure tank rated pressure value Pb, air pressure tank gas chamber nominal volume Vb, air pressure tank is specified
Temperature Tb, environment temperature is T (t), and t is time variable, and ρ is fluid density.
When flow control system Relative steady-state：Ductwork pressure value is P, and frequency converter output frequency is F, passes in and out fluid flow
For Q, environment temperature T, air pressure tank chamber volume is V1, liquid chamber volume is V2, the unit of above-mentioned all amounts is international unit.
Define the t=0 moment be system with the last moment of frequency F stable operations, that is, exist：
Assuming that (0, Td] running frequency that pumps in the time is：F (t)=F+ Δ F, Δ F is frequency disturbance increment, usual feelings
Under condition | Δ F | ＜ ＜ F；TdFor pre-defined observation interval, for the time value more than 0, according to flow control system
Energy index is different and artificially determines；Then pressure value is p (t)=P+ Δ p (t), and Δ p (t) is pressure oscillation value caused by Δ F；Pump
Rate of discharge is q (t)=Q+ Δ q (t), and Δ q (t) is pump discharge flow undulating value caused by Δ F；Air pressure tank rate of discharge is
q2(t)=Q+ Δs q2(t), Δ q2(t) it is air pressure tank rate of discharge undulating value caused by Δ F；It is controlled by motor frequency conversion it is found that pumping
The relationship of output power be：
Wherein：ρ × the q (t) on the equation left side × p (t) is the shaft power of pump；η is the efficiency of pump；For the output power of motor；S is revolutional slip； R1,R2,X1σ,X2σ,m1,For pump motor
Intrinsic parameter；
Since pump motor uses variable frequency regulating speed control, so s is held essentially constant.It enables：
K is only related with motor structural parameters itself, unrelated with flow, pressure.So formula (1) can be reduced to：
Q (t) p (t)=k η f (t)2/ρ (3)
Enable k'=η k/ ρ.Then in t=0, have：
QP=k'F2 (4)
In t ∈ (0, Td], q (t)=Q+ Δ q (t), f (t)=F+ Δs F and p (t)=P+ Δ p (t) are substituted into formula (4)：
(Q+ Δ q (t)) (P+ Δ p (t))=k'(F+ Δ F)2 (5)
It is unfolded (5), and arranges：
PQ+Q Δs p (t)+P Δ q (t)+Δ q (t) Δ p (t)=k'(F2+2FΔF+ΔF2) (6)
(4) substitution (6) can be obtained：
Q Δs p (t)+P Δ q (t)+Δ q (t) Δ p (t)=k'(2F Δ F+ Δs F2) (7)
Since there are the big inertia damping links of air pressure tank, then in t ∈ (0, Td] changes in flow rate amount Δ q (t) draws in the short time
Pressure variety Δ p (t) very littles risen meet：
| Δ p (t) | ＜ ＜ P (8)
It is obtained so arranging (7)：
Q Δs p (t)+P Δ q (t)=k'(2F Δ F+ Δs F2) (9)
By formula (9) divided by (4) and consider | Δ F | ＜ ＜ F can be obtained：
Due in t ∈ (0, Td] have | Δ p (t) |《P, i.e. ductwork pressure are kept approximately constant, and are not changed in pipe resistance characteristic
In the case of, the rate of discharge variation delta q of air pressure tank2(t) 0 ≈, i.e. q2(t)≈Q.Have according to air pressure tank kinetics equation：In t
∈(0,Td], the volume change of air pressure tank liquid chamber is：
So t ∈ (0, Td] liquid chamber volume is：
Because V is remained unchanged, thus chamber volume is：
In t ∈ (0, Td] in the time, environment temperature remains unchanged, then from equation for ideal gases：
(13) are substituted into (14) and are arranged：
Enable Δ pa(t)=pa(t)-pa(0) it is air pressure tank air chamber pressure variable quantity, then：
According to hydraulic principle it is found that ductwork pressure variable quantity is：
By pa(0)=P substitutes into formula (17), can obtain：
Simultaneous (18) and (10) simultaneously arrange：
It enables：Then have：Y'(t)=Δ q (t), thus have：Y (0)=0, can to formula (19) arrangement
：
The differential equation (20) is arranged and considers that Δ F ＜ ＜ F, 2 × Δ F ＜ ＜ F can be obtained：
Solving (21) can obtain：
It willIt substitutes into formula (22) and arranges：
In t ∈ (0, Td], due to | Δ F |《F and | Δ p (t) |《P, according to (5) it is found that Δ q (t) ＜ ＜ Q, so having：
Below for Δ F withSyntactics discuss analysis：As Δ F ＞ 0, due to f (t)
=F+ Δ F ＞ F, thus q (t)=Q+ Δ q (t) ＞ Q, so there is Δ q (t) ＞ 0；Similarly, as Δ F ＜ 0, due to f (t)=F+
Δ F ＜ F, thus q (t)=Q+ Δ q (t) ＜ Q, so there is Δ q (t) ＜ 0；So：Δ F and Δ q (t) jack per lines, that is, Δ F with
Y (t) jack per lines.So having：
Again due in t ∈ (0, Td], the right end of formula (24) meets：Qt ＞ 0, so having：
So formula (24) can arrange：
(27) are solved equation to obtain：
Again because of Δ q (t)=y'(t), have：
Because of air pressure tank No leakage, then from equation for ideal gases：
Simultaneous formula (29) and (30), and arrange：
Due to t ∈ (0, Td], if to TdSelection meets inequality：
Then (31) are carried out Taylor series expansion and arranged to obtain：
Due toSo expression formula (33) can be approximately：
Due to parameter, Δ q (t), Q, F, Δ F, Pb、Vb、Tb, T and t be observable quantity and known quantity, thus pass through acquisition
The value of changes in flow rate amount Δ q (t) can values of pressure p of the on-line measurement outflow control system in stable state size.
The pressure P at flow control system arbitrary Relative steady-state moment can be found out according to formula (34).At the same time, flow control
System output flow Q processed can be obtained by flow sensor, and then get flow control system in the operating point of Q-H planes.
Fig. 3 show governor impeller Effec-tive Function area schematic diagram, and the Effec-tive Function section of pump is rated frequency fNLift it is special
Linearity curve HN, low-limit frequency fminHeat-capacity curve Hmin, similar operating condition parabola li1, similar operating condition parabola li2It surrounds
Fan annular region ABCD.If pump is in region ABCD in the characteristic operating points Q-H, pump is in Effec-tive Function；Conversely,
Pump is in non-efficient operating status.
Since flow control system realizes flow-rate adjustment control mode, thus different running frequency situations using frequency control
The heat-capacity curve of lower pump has translation feature.Pump operation section distribution situation is described in detail with reference to Fig. 3.
(1) flow control system output flow is Q1：
Assuming that when the running frequency of front pump is f1, then the heat-capacity curve pumped is H1, flow Q1Corresponding operating point
Pressure value be P1.From the figure 3, it may be seen that when front pump is in efficient region ABCD.If certain moment flow control system is because of other
When factor causes pressure to reduce (for example, fluid pipeline change causes pipe resistance to reduce, and liquid enters reaction groove tank of low lift etc.),
Then maintaining output flow Q1In the case of constant, the running frequency of pump must be reduced, it is assumed that the running frequency of pump is at this time
f2, heat-capacity curve is switched to H2.From the figure 3, it may be seen that characteristic curve H2Middle flow is Q1The pressure value of corresponding operating point is
P2, at this time pump operation point be not in efficient region ABCD, the inefficiency of pump, fever is serious.
(2) flow control system output flow is by Q1It is adjusted to Q2
Assuming that when the running frequency of front pump is f1, then the heat-capacity curve pumped is H1, flow Q1Corresponding operating point
Pressure value be P1.If certain moment flow control system setting output flow increases to Q2, then not due to the pipe resistance characteristic of system
Become, increasing output flow necessarily causes pipe resistance to increase, and must improve the running frequency of pump, it is assumed that the running frequency of pump is at this time
f3, heat-capacity curve is switched to H3.From the figure 3, it may be seen that characteristic curve H3Middle flow is Q2The pressure value of corresponding operating point is
P3, at this time pump operation point be not in efficient region ABCD, the inefficiency of pump, fever is serious.Similarly, it is assumed that flow control system
System is not currently in high efficiency range operation, then due to the change of external operating mode, system operating point migrates, it is possible in height
Imitate section operation.
By above-mentioned analysis it is found that the traffic coverage of flow control system pump is not to be constantly in high efficient district, with
The variation of output flow and system pipes resistance and change, in order to realize the safe and reliable operation of flow control system, then have to pair
Flow control system is efficiently controlled.
Using the flow control system singly pumped, due to output flow can with the difference of technique flow or reaction process and
Variation, pipe resistance characteristic with variation can also lead to that Effec-tive Function section can not possibly be constantly in.Output stream may be used thus
Amount range is identical, and lift is divided into three large, medium and small pump redundancy of effort modes, and control strategy is dynamic according to the operating mode of system
The optimal pump work of state efficiency of selection, and in addition two pumps are then in stand-by state.On the one hand, the efficiency of system is improved；Separately
On the one hand, the reliability of system is improved.
The present invention provides a kind of flow control system and efficient control methods, include the following steps：
(1) with sampling period TsFlow control system output flow value and frequency converter output frequency are adopted for interval
First time sampled value is labeled as q (1) and f (1) by sample；Label present sample number is k；
Define flow error e (k)=Qset-q(k)；Wherein, e (i) |＜=0 i=0；QsetTo set output flow value；q(k)
Flow value when for sampling number being k, the output frequency value of f (k) is sampling number when being k frequency converter；f(i)|＜=0 i=0；
Enable k=1；
(2) t=kT is found out by pid control algorithmsThe output frequency value f (k) of moment frequency converter=f (k-1)+Kp[e(k)-
e(k-1)]+Kie(k)+Kd[e(k)-2e(k-1)+e(k-2)]；
Wherein, e (k-1), f (k-1) are respectively t=(k-1) TsThe flow error at moment and the output frequency of frequency converter；e
(k-2) it is t=(k-2) TsThe flow error at moment；
Kp、KiAnd KdProportionality coefficient, integral coefficient and differential coefficient in respectively preset pid algorithm；
More new variables enables e (k-2)=e (k-1), e (k-1)=e (k), f (k-1)=f (k)；
(3) the flow value array { q (i) } being made of N number of element and frequency converter output frequency array { f (i) } are established,
Wherein i={ k-N+1, k-N+2 ... k }, N be it is preset be more than 1 positive integer, k is present sample number；q(i)
|＜=0 i=0, f (i) |＜=0 i=0；
(4) judge whether flow control system is in metastable state, the definition of metastable state is：Calculate { q
(i) } average valueAnd it solvesJudge whether to meet：σq≤εq,
In：εqTo set positive value, can be set according to real system, for example 0.05 or 0.1 can be taken.If it is satisfied, then thinking
Flow control system is in metastable state, enters step (5)；Otherwise, flow control system plays pendulum, and is transferred to
Step (20)；
(5) judgeWhether it is true (wherein θ be setting positive value, can be carried out according to real system
Setting, for example 0.01 or 0.03 can be taken；QsetTo set output flow).If it is satisfied, the pump M that explanation is currentj(j=1,
2,3) flow control requirement can be met, be transferred to step (7)；Otherwise, illustrate the pump M of work at presentj(j=1,2,3) H-MaxIt is too small, output flow requirement cannot be met, entered step (6).
(6) control switch S is enabledj=0, that is, control MjIt is out of service；Meanwhile enabling control switch Sj+1=1, that is, it is big to control lift
The pump M of level-onej+1(j+1<=3) it works, is transferred to step (20).If work at present MjIt has been lift maximum pump, then has illustrated the stream
There are problems for the selection of amount control system pump, cannot be satisfied system requirements, and control method of the invention is not applicable in this case.
(7) average value of frequency converter output frequency is solved
(8) t=0 is denoted as with blaze at this time, gives output frequency one smaller disturbance quantity Δ F, i.e. f (mTs)=F+ Δs F；
(9) m=1 is enabled；
(10) judge whether m ＞ M are true.It sets up, then enters step (13)；Otherwise, in t=mTsMoment, sampling flow
Value is denoted as q (m)；It obtains
(11) judge(α is setting positive value, can be set according to real system, for example can take 0.01
Or it is 0.1) whether true.It is invalid, it is transferred to step (20)；Otherwise, it enters step (12)；
(12) more new variables：Enable m=m+1；Return to step (10).
(13) to M point (Ts,Δq(1)),(2Ts,Δq(2)),…,(MTs, Δ q (M)) and fitting a straight line is carried out, and obtain
The slope K of the straight line.
(14) according to formulaEstimate pressure P.
(15) have the characteristics that translate according to point r (Q, P) and governor impeller its lift characteristic, obtain the pump of passing point r (Q, P)
MjThe Q-H heat-capacity curves of (j=1,2,3), are denoted as
(16) heat-capacity curve is solvedCorresponding frequency Fj。
(17) it solvesTo similar parabolaIntersection point ajAnd bj(j=1,2,3).
(18) it calculates and meetsThe pump M of conditionjEffec-tive Function performance parameter ψj。ψjFor nonnegative number, and
ψjBigger, then performance is better.Wherein：ψjFor parameter aj、bjWith r (Q, P) and corresponding high efficiency range AjBjCjDjFunction.For example, meter
It calculates r (Q, P) and arrives region AjBjCjDjGeometric center point distance derivative, or calculate r (Q, P) whether be in AjBjCjDjArea
In domain and r (Q, P) arrives AjBjCjDjThe weighting function of the distance of boundary curve, or to calculate whether r (Q, P) is in
AjBjCjDjIn region and r (Q, P) arrives aj、bjThe weighting function (j=1,2,3) of distance.
(19) ψ is soughti=max { ψ1,ψ2,ψ3Corresponding i (i=1,2,3).Controller is by corresponding switch Si(t)=1,
Su=0 (u=1,2,3 ∩ u ≠ i), to select the pump work of suitable lift, improve the efficiency of system, and enter step (20).
(20) k=k+1 is enabled；It is sampled next time, and marks the sampled value of output flow value and frequency converter output frequency
For q (k) and f (k)；Return to step (2).
Embodiment is not construed as limitation of the present invention, any spiritual improvements introduced based on the present invention, all Ying Ben
Within the protection domain of invention.
Claims (7)
1. a kind of flow control system and efficient control method, it is characterised in that：Its step are as follows：
1) pressure value P of the flow control system in stable state and t ∈ [0, T are establishedd] changes in flow rate amount Δ q1(t) relational expression：Wherein；P is ductwork pressure value, and F is frequency converter output frequency, and Q is disengaging liquid flow
Amount, T is environment temperature, TbFor air pressure tank rated temperature, VbFor air pressure tank gas chamber nominal volume, PbFor the specified pressure of air pressure tank gas chamber
Power, t are time variable, TdFor pre-defined observation interval, Δ F is frequency disturbance increment；
2) by obtaining slope in the relational expression of foundation in step 1)
3) with sampling period TsFlow control system flow value and frequency converter output frequency are sampled for interval, and obtain stream
Magnitude q (k) and output frequency f (k), wherein k are sampling number；
4) t=kT is found out by pid control algorithmsThe output frequency value f (k) of moment frequency converter=f (k-1)+Kp[e(k)-e(k-
1)]+Kie(k)+Kd[e(k)-2e(k-1)+e(k-2)]；
Wherein, e (k-1), f (k-1) are respectively t=(k-1) TsThe flow error at moment and the output frequency of frequency converter；e(k-2)
For t=(k-2) TsThe flow error at moment, Kp、KiAnd KdProportionality coefficient, integration system in respectively preset pid algorithm
Number and differential coefficient, flow error e (k)=Qset-q(k)；Wherein, e (i) |＜=0 i=0；QsetTo set output flow value；q(k)
Flow value when for sampling number being k, the output frequency value of f (k) is sampling number when being k frequency converter；f(i)|＜=0 i=0；
5) and according to the flow value q (k) and output frequency f (k) sampled, the flow value array { q being made of N number of element is established
(i) } and frequency converter output frequency array { f (i) }, wherein i={ k-N+1, k-N+2 ... k }, N is preset is more than
1 positive integer, q (i) |＜=0 i=0, f (i) |＜=0 i=0；
6) average value of flow value array { q (i) } is obtainedAnd judge whether flow control system is in phase
To stable state；
7) when determining that flow control system is in metastable state, then judgeIt is whether true,
Middle θ is setting positive value, QsetTo set output flow；
8) in determinationWhen invalid, then it is considered as the pump M of work at presentjH-MaxIt is too small, j=
1,2,3, enable control switch Sj=0, control MjIt is out of service；Meanwhile enabling control switch Sj+1=1, control the pump of the big level-one of lift
Mj+1Work, j+1<=3, update k=k+1；It is sampled next time, and marks output flow value and frequency converter output frequency
Sampled value is q (k) and f (k), repeats above step；
9) in determinationWhen establishment, then the average value of frequency converter output frequency is obtained
10) t=0 is denoted as with blaze at this time, gives output frequency one smaller disturbance quantity Δ F, f (mTs)=F+ Δs F；
11) judge whether m ＞ M are true, wherein
12) if m ＞ M are invalid, in t=mTsMoment, sample streams magnitude q (m)；It obtainsAnd judgeWherein α is setting positive value；If so, m=m+1 is then updated, and re-starts m ＞ M judgements, if not,
Then update k=k+1；It is sampled next time, and it is q (k) and f to mark output flow value and the sampled value of frequency converter output frequency
(k), above step is repeated；IfIt sets up, then to M point (Ts,Δq(1)),(2Ts,Δq(2)),…,(MTs,
Δ q (M)) fitting a straight line is carried out, and obtain the slope K of the straight line；
13) pass through acquisition in step 2)Estimate pressure P；
14) have the characteristics that translate according to point r (Q, P) and governor impeller its lift characteristic, obtain the pump M of passing point r (Q, P)jQ-
H heat-capacity curvesJ=1,2,3；And obtain pump MjQ-H heat-capacity curvesCorresponding frequency Fj；
15) pump M is obtainedjQ-H heat-capacity curvesTo similar parabolaIntersection point ajAnd bj, j=1,2,3；
16) it calculates and meetsThe pump M of conditionjEffec-tive Function performance parameter ψj, ψjFor nonnegative number, ψjFor intersection point
aj、bjWith r (Q, P) and corresponding high efficiency range AjBjCjDjFunction；
17) ψ is soughti=max { ψ1,ψ2,ψ3Corresponding i, i=1,2,3, controller is by corresponding switch Si(t)=1, Su=0, u
=1,2,3 ∩ u ≠ i, and update k=k+1；It is sampled next time, and marks output flow value and frequency converter output frequency
Sampled value is q (k) and f (k), repeats above step.
2. a kind of flow control system according to claim 1 and efficient control method, it is characterised in that：Packet in step 1)
Include following steps：
1. establishing flow control system water pump output power equation：
Wherein：η is the efficiency of pump, i.e., the ratio between motor effective power and shaft power, s are revolutional slip, R1,R2,X1σ,X2σ,m1,For the intrinsic parameter of pump motor；
2. carrying out small signal disturbance to the equation of step 1., obtains relational expression and be reduced to Q Δs p (t)+P Δs q1(t)+Δq1(t)Δp
(t)=k'(2F Δ F+ Δs F2), wherein：q1(t)=Q+ Δs q1(t), f (t)=F+ Δs F, p (t)=P+ Δ p (t), k'=η k/ ρ,
3. in t ∈ [0, Td] obtain system small-signal model equation：
4. obtaining in t ∈ [0, Td], the volume change of air pressure tank liquid chamber；
And thus obtain t ∈ [0, Td] when liquid chamber volume
Chamber volumeAnd the air chamber pressure variable quantity of air pressure tank is obtained according to equation for ideal gasesAnd thus obtain ductwork pressure variable quantityAnd
In pa(0)=P when obtain
5. 3. 4. being obtained with step according to stepAnd it finally obtains
3. a kind of flow control system according to claim 1 and efficient control method, it is characterised in that：The high efficient district
Domain ABCD is rated frequency fNHeat-capacity curve HN, low-limit frequency fminHeat-capacity curve Hmin, similar operating condition parabola
li1, similar operating condition parabola li2The fan annular region surrounded.
4. a kind of flow control system according to claim 1 and efficient control method, it is characterised in that：It is obtained in step 6)
Take the average value of flow value array { q (i) }And it solvesJudge whether full
Foot：σq≤εq, wherein：εqTo set positive value, if satisfied, then thinking that flow control system is in stable state.
5. a kind of flow control system according to claim 1 and efficient control method, it is characterised in that：ψjFor r (Q, P)
To region AjBjCjDjGeometric center point distance derivative, wherein j=1,2,3.
6. a kind of flow control system according to claim 1 and efficient control method, it is characterised in that：ψjFor r (Q, P)
Whether A is injBjCjDjIn region and r (Q, P) arrives AjBjCjDjThe weighting function of the distance of boundary curve, wherein j=1,2,3.
7. a kind of flow control system according to claim 1 and efficient control method, it is characterised in that：ψjFor r (Q, P)
Whether A is injBjCjDjIn region and r (Q, P) arrives aj、bjThe weighting function of distance, wherein j=1,2,3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810076141.XA CN108490989A (en) | 2018-01-26 | 2018-01-26 | A kind of flow control system and efficient control method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810076141.XA CN108490989A (en) | 2018-01-26 | 2018-01-26 | A kind of flow control system and efficient control method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108490989A true CN108490989A (en) | 2018-09-04 |
Family
ID=63343736
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810076141.XA Withdrawn CN108490989A (en) | 2018-01-26 | 2018-01-26 | A kind of flow control system and efficient control method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108490989A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110953169A (en) * | 2019-12-03 | 2020-04-03 | 温州大学 | Control method of parallel variable-frequency constant-voltage control system |
| CN114357664A (en) * | 2022-03-15 | 2022-04-15 | 河北建投水务投资有限公司 | Modeling method and system for mathematical model of variable-frequency speed-regulating water pump |
-
2018
- 2018-01-26 CN CN201810076141.XA patent/CN108490989A/en not_active Withdrawn
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110953169A (en) * | 2019-12-03 | 2020-04-03 | 温州大学 | Control method of parallel variable-frequency constant-voltage control system |
| CN114357664A (en) * | 2022-03-15 | 2022-04-15 | 河北建投水务投资有限公司 | Modeling method and system for mathematical model of variable-frequency speed-regulating water pump |
| CN114357664B (en) * | 2022-03-15 | 2022-07-01 | 河北建投水务投资有限公司 | Modeling method and system for mathematical model of variable-frequency speed-regulating water pump |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108287571B (en) | Method for judging running interval of flow control system pump | |
| CN108415467A (en) | A kind of flow control system and Optimal Control Strategy | |
| CN106014999B (en) | A kind of water-ring vacuum pump assembly Effec-tive Function control method and control system | |
| CN108445921A (en) | A kind of section recognition methods of flow control system pump operation | |
| CN105302984B (en) | A kind of earth source heat pump unit modeling and simulating method | |
| CN202165326U (en) | Industry centrifugal compressor surge-proof control device | |
| US20080264086A1 (en) | Method for improving efficiency in heating and cooling systems | |
| CN108490989A (en) | A kind of flow control system and efficient control method | |
| CN108415466A (en) | A kind of flow control system pump operation method of interval estimation | |
| CN107131131A (en) | A kind of flow equilibrium method of water pump serial-connection system | |
| CN108490987A (en) | A kind of flow control system pump operation interval prediction method | |
| CN108490988A (en) | A kind of flow control system pump operation section method of discrimination | |
| CN102200120A (en) | Precomputation method of frequency conversion transformation working condition points of condensate pump for thermoelectric generating set | |
| CN110736129A (en) | Intelligent balance regulation and control system and method for urban heat supply pipe networks | |
| Sapaev et al. | Research energy and resource saving operating modes of the pump unit | |
| Zhang et al. | CFD research on runaway transient of pumped storage power station caused by pumping power failure | |
| CN105543443B (en) | A method of being used for coal gas of converter pressuring machine antisurge control system | |
| CN101871225B (en) | Water supply system with independently controlled temperature, flow and pressure | |
| CN108268060B (en) | A kind of flow control methods | |
| CN100561043C (en) | A kind of boiler feedwater flow control and pressure compensating system | |
| CN108415468A (en) | A kind of flow control system and optimal control method | |
| CN108490990A (en) | A kind of flow control system and efficient control strategy | |
| CN108319137A (en) | A kind of flow control system pressure prediction method | |
| CN108333932A (en) | A kind of flow control system pressure method of estimation | |
| CN108426668A (en) | A kind of flow control system pressure measurement method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20180904 |
|
| WW01 | Invention patent application withdrawn after publication |