CN113177301B - Parameter optimization method for intelligent secondary water supply device for water hammer prevention - Google Patents

Parameter optimization method for intelligent secondary water supply device for water hammer prevention Download PDF

Info

Publication number
CN113177301B
CN113177301B CN202110394232.XA CN202110394232A CN113177301B CN 113177301 B CN113177301 B CN 113177301B CN 202110394232 A CN202110394232 A CN 202110394232A CN 113177301 B CN113177301 B CN 113177301B
Authority
CN
China
Prior art keywords
water
water tank
max
valve
value
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.)
Active
Application number
CN202110394232.XA
Other languages
Chinese (zh)
Other versions
CN113177301A (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.)
Fengtai Water Technology Hangzhou Co ltd
Original Assignee
Fengtai Water Technology Hangzhou Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fengtai Water Technology Hangzhou Co ltd filed Critical Fengtai Water Technology Hangzhou Co ltd
Priority to CN202110394232.XA priority Critical patent/CN113177301B/en
Publication of CN113177301A publication Critical patent/CN113177301A/en
Application granted granted Critical
Publication of CN113177301B publication Critical patent/CN113177301B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Electricity, gas or water supply
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/04—Constraint-based CAD
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A20/00—Water conservation; Efficient water supply; Efficient water use
    • Y02A20/152—Water filtration

Abstract

The invention discloses a parameter optimization method for a waterproof hammer intelligent secondary water supply device, which comprises the following steps: s1, establishing a water inlet model of the water tank; s2, collecting data and checking the water inlet model of the water tank; s3 determining an optimization variable, a constraint condition and an optimization target; s4 determines the optimization parameters by calculation. Above-mentioned technical scheme is through establishing the water tank model of intaking, and the cooperation optimizes the control effect of intaking of formula of relying on oneself liquid level control valve, can moisturizing as early as possible and guarantee the water tank regulation function, can reduce again and intake the sudden change and avoid, and the valve is opened or the valve of closing arouses the pipeline water hammer too fast to reduce the loss of valve, extension valve life.

Description

Parameter optimization method for intelligent secondary water supply device for water hammer prevention
Technical Field
The invention relates to the technical field of urban water supply, in particular to a parameter optimization method for an intelligent secondary water supply device for a water hammer.
Background
The water tank of the pump room in the district II is usually controlled by a ball float valve to feed water, has the problems of small water level control range, frequent action, water feeding vibration noise and the like, and is gradually replaced by a self-operated liquid level valve.
The data shows that in residential areas, residential water is concentrated in the morning and evening peaks, and different cells have different peak water consumption curves. If the self-operated liquid level control valve is put into operation according to fixed parameters, the self-operated liquid level control valve takes the difference in the aspects of water level control range, safe liquid level of a water tank, valve opening action amplitude and the like. Such as opening or closing the valve too quickly, causes a large water hammer, resulting in accelerated wear of the valve member.
Chinese patent document CN111637272A discloses a "self-operated liquid level adjusting device". Including the storage tank and receive the jar, through float governing valve intercommunication between storage tank and the receipt jar, the float governing valve includes flange one, float, valve body, flange two, sealed pad and balance pipe flange, flange one is through flange mouth and storage tank and receive jar all intercommunication, the valve body of flange one and the valve body of flange two link up, the float is installed in the valve body of flange two. According to the technical scheme, when the water consumption peak is not considered, the valve is opened or closed too fast, a large water hammer is caused, and the accelerated loss of the valve is caused.
Disclosure of Invention
The invention mainly solves the technical problem that the loss is accelerated due to the fact that a liquid level valve is switched too fast when the water consumption peak is not considered in the original technical scheme, and provides a parameter optimization method for an intelligent water supply device for preventing a water hammer.
The technical problem of the invention is mainly solved by the following technical scheme: the invention comprises the following steps:
s1, establishing a water inlet model of the water tank;
s2, collecting data and checking the water inlet model of the water tank;
s3 determining an optimization variable, a constraint condition and an optimization target;
s4 determines the optimization parameters by calculation.
Preferably, the step S1 specifically includes:
s1.1, setting water inlet model parameters of a water tank;
s1.2, calculating the representation of the adjusting action of the control valve in steady-state and transient hydraulic power;
s1.3, establishing a water inlet model of the water tank.
Preferably, the step S1.1 specifically includes: setting the deformation x of the spring K1, the height H of the water tank liquid level when water is in and out, and the height H of the water tank liquid level when the valve is closedmaxThe pressure of the cavity A is
FA＝PS＝ρghS
P is the pressure acting on the cavity A, S is the area of the diaphragm M, rho is the density of water, g is 9.8N/kg, and the elastic force of the spring K1 is
FK1＝Kx
Where K is the spring constant of spring K1
Balance equation of membrane in valve body:
pgh(t)S＝Kx(t) (1)
wherein t is time, h (t) is liquid level at time t, x (t) is spring deformation, the difference value between the current liquid level of the water tank and the highest liquid level influences the opening degree of the valve, and the larger the difference value is, the larger the opening degree is:
Figure BDA0003017961420000031
wherein
Figure BDA0003017961420000032
h(t)≤Hmax。
Preferably, the step S1.2 specifically includes: the control valve modulation is expressed in both steady-state and transient hydraulic calculations
Figure BDA0003017961420000033
Wherein QinΔ H and SinThe flow passing through the valve, the water head loss at the valve and the sectional area of a water inlet pipeline are respectively, and the delta H is equal to the water pressure of the pipe network and the height of the water tank, and is simplified and regarded as a constant.
Preferably, the step S1.3 specifically includes: cVThe control valve is a stop valve for the flow coefficient of the valve, and the relative flow coefficient is in direct proportion to the opening degree of the valve, namely, the following conditions are met:
Figure BDA0003017961420000034
wherein epsilon is a coefficient, and D is the inner diameter of the water inlet pipeline;
from conservation of mass:
Figure BDA0003017961420000035
wherein SWater tankIs the bottom area of the water tank C, Qout(t) is the outflow rate of the cistern, since QinIt is necessary to flow through a section of pipe network from the control valve to the top of the tank, so
Figure BDA0003017961420000036
However, since the difference between the two is small, the processing is simplified
Figure BDA0003017961420000037
Substitution
Figure BDA0003017961420000038
Therefore, the formula (2), (3), (4) and (5) are simplified to obtain the formula (6):
Figure BDA0003017961420000039
the formulas (1) and (6) form a water inlet model of the water tank, and the water inlet model is along with the water outlet flow QoutAnd (t) changing the liquid level of the water tank according to the water inlet model.
Preferably, the step S2 of checking the water inlet model of the water tank includes:
s2.1 Water consumption at the time of collecting the cell Qout(t), the inner diameter D of the water inlet pipeline of the community. Bottom area S of water tank CWater tankArea S of diaphragm M, height H of water level of water tank when valve is closedmaxHeight H of water tank liquid level when valve is fully openedminCalculating the sectional area S of the water inlet pipeline by using the elastic coefficient K of the spring K1inAnd head loss Δ H at the valve;
s2.2 Collection district Inlet pipe inflow Qin(t), calculating the inflow water flow rate;
and S2.3, checking the water inlet model of the water tank.
Preferably, the step S2.3 specifically includes:
s2.31 sets a value range of the coefficient ∈ [ epsilon min, epsilon max ], an initial value ∈ 0 ∈ min and a step length Δ ∈ (epsilon max-epsilon min)/100, and sets an optimization target that the actually measured liquid level value h (t) is most consistent with the simulation value hs (t), that is, Jh ∑ h (t) -hs (t) | is minimized;
s2.32 respectively substituting epsilon 0, epsilon 1, … … epsilon 100 and related data into an equation (6), carrying out 24-hour simulation, recording simulation data, namely hs (t) at each moment, and calculating to obtain objective function values Jh0 and Jh1 … … Jh 100;
s2.33 finds the minimum value Jhmin from Jh0 and Jh1 … … Jh100, the epsilon value at the moment is the fixed value of the model rate, and the water inlet model check is completed.
Preferably, the step S3 of determining the optimization variables, the constraints, and the optimization target specifically includes:
s3.1 when the diaphragm of the self-operated liquid level control valve is selected, the adjustable variables are the spring deformation x (t) and the height H of the liquid level of the water tank corresponding to the closing and opening valvesmaxAnd HminThe two quantities are related to the spring elastic coefficient K, and the optimized variable is the spring elastic coefficient K;
s3.2 constraint conditions: qin(t) not greater than MaxQ of maximum inflow of pipe networkinThe water level in the water tank can not be lower than the lowest level HminNamely:
0≤Qin(t)≤MaxQin (7)
Hmin≤h(t)≤Hmax (8)
s3.3 optimization goal: max (| Q) corresponding to K valueinReducing the values of' (t) | and Max (H-H (t)), adding weight to construct an optimized objective function:
J＝a*Max(|Qin′(t)|)+b*Max(H-h(t)) (9)
where a, b are weights such that J is minimized.
Preferably, the step 4 specifically comprises:
s4.1, determining the range, initial value and step length of the optimized variable K, and according to the maximum displacement of the valve rod, the lowest liquid level H of the water tank during water inlet and outletWater tank minDetermining Kmax(ii) a Minimum level H according to maximum displacement of valve stem 1/2Water tank minDetermining Kmin(ii) a Then there is Kmin≤K≤KmaxSetting K as initial Kmin, setting step length delta K as (Kmax-Kmin)/100;
s4.2, calculating an optimization target subentry index value under the current K value;
s4.3 if K<KmaxIf yes, returning to step S4.2;
s4.4, checking constraint conditions, drawing an optimization target subentry index curve under different K values, and aiming at Q under different K valuesin(t) h (t), checking whether the constraint condition is met, if not, rejecting the K value, and respectively drawing the corresponding Max (| Q) values under different K valuesin' (t) |), Max (H-H (t));
s4.5, calculating objective function values J under different K values, determining weights a and b, and firstly, calculating Max (| Q)in' (t) |, Max (H-H (t)) are respectively standardized; distributing 0.5 to the processed data weight factors, and calculating objective function values J under different K values; and searching Jmin, wherein the K value corresponding to the Jmin is the optimal solution.
Preferably, the step S4.2 specifically includes:
s4.21 takes 1 day with the maximum daily water consumption of the community as the basis, and assumes that 0 is the initial state, H (0) is Hmax，Qout(0) When K is 0 and K is Kmin, h' (0) is obtained by substituting formula (6);
s4.22h' (0) is multiplied by the time interval delta t to obtain h (1); then, Q is obtained from the formula (5)in(0)；
S4.23 relating h (1) to Q at the next momentout(1) The value is substituted for formula (6);
s4.24 repeating steps S4.21, S4.22 and S4.23, calculating all h' (t), h (t) and Q under the same K valuein(t)；
S4.25 calculating the corresponding | Q at each moment under the current K valuein' (t) | and H-H (t), and records Max (| Q)in' (t) |) and Max (H-H (t)).
The invention has the beneficial effects that: through establishing the water tank model of intaking, the cooperation optimizes the control effect of intaking of formula of relying on oneself liquid level control valve, can supply water as early as possible and guarantee the water tank regulation function, can reduce the sudden change of intaking again and avoid, and the valve is opened or the valve of closing arouses the pipeline water hammer too fast to reduce the loss of valve, extension valve life.
Drawings
FIG. 1 is a schematic view of the regulation of the water intake level in a cistern according to the present invention.
Fig. 2 is a flow chart of the present invention.
FIG. 3 shows the maximum daily water consumption Q of a certain district according to the present inventionoutCurve line.
FIG. 4 is a graph of the variation of the simulated liquid level and the actual liquid level of the water tank at the highest day according to the present invention.
FIG. 5 shows Max (| Q) values at different K values according to the present inventionin' (t) |) data map.
FIG. 6 is a graph of Max (H-H (t)) data at different K values according to the present invention.
Fig. 7 shows the variation of the function J for different values of K according to the invention.
FIG. 8 is a graph of tank level change for different K values of the present invention.
FIG. 9 shows Q at different K values according to the present inventionin' (t) variation curve.
Detailed Description
The technical scheme of the invention is further specifically described by the following embodiments and the accompanying drawings.
Example (b): the parameter optimization method for the intelligent two-water-supply device for the waterproof hammer in the embodiment is shown in fig. 2 and comprises the following steps:
s1, establishing a water inlet model of the water tank, specifically comprising:
s1.1, setting water tank water inlet model parameters, as shown in a water tank water inlet regulation schematic diagram of FIG. 1, setting a deformation x of a spring K1, and setting a water tank liquid level h (indicating the height from the bottom of a tank to the liquid level) when water enters or exits. Adjusting the set screw B can set the spring K1 value, the highest and lowest liquid level values, and the height of the water tank liquid level when the valve is closed is set to be Hmax。
According to the principle of a communicating vessel, the water tank C is communicated with the valve cavity A, and the force acting below the diaphragm M is FA+FAtmosphere (es)In which F isAPS ρ ghS, P is the pressure applied to the chamber a, S is the area of the diaphragm M, ρ is the density of water, and g is 9.8N/kg (usually, g is 10N/kg); according to Hooke's law FK1Kx, where K is a springK1, the pressure acting over the membrane M being FK1+FAtmosphere (es)(ii) a When the liquid level of the water tank reaches the maximum value, the pressure F of the cavity AAMaximum, spring force F of spring K1K1Cavity A pressure FAThe two forces are in a balanced state, and the main valve is in a closed state; fA+FAtmosphere (es)＝FK1+FAtmosphere (es)Elimination of FAtmosphere (es)，FA＝FK1. When the liquid level in the water tank drops, FABecome smaller, FK1The main valve is pushed to open until the lowest point of the set value of the liquid level of the water tank, and the main valve is fully opened. On the contrary, when the liquid level of the water tank is gradually increased, FAThe valve is gradually closed until the valve is completely closed.
Therefore, the balance equation of the membrane in the valve body can be obtained:
pgh(t)S＝Kx(t) (1)
wherein t is time, h (t), x (t) are liquid level and spring deformation at time t. Assuming the initial state (0 hours), the valve is closed, and the liquid level of the water tank is the highest liquid level Hmax. When the user's water progressively increases, the water tank liquid level reduces, and the valve aperture increases, and the water tank is intake and is increased. When the liquid level falls to the lowest liquid level HminWhen the water inlet valve is in a fully-opened state, the water inlet amount reaches the maximum value and is not increased any more. When the water consumption of the user is reduced, the liquid level of the water tank rises, the opening degree of the valve is reduced, and the water inlet of the water tank is reduced. When the liquid level reaches the highest liquid level HmaxWhen the valve is closed, the valve is closed.
Because the difference of the current liquid level of water tank and highest liquid level is influencing the opening degree size of valve, and the difference is big more, and the aperture is big more:
Figure BDA0003017961420000081
wherein
Figure BDA0003017961420000082
h(t)≤Hmax。
S1.2, calculating the representation of the regulating action of the control valve in steady-state and transient hydraulic power, and specifically comprising the following steps: the control valve modulation is expressed in both steady-state and transient hydraulic calculations
Figure BDA0003017961420000083
Wherein QinΔ H and SinThe flow passing through the valve, the water head loss at the valve and the sectional area of a water inlet pipeline are respectively, and the delta H is equal to the water pressure of the pipe network and the height of the water tank, and is simplified and regarded as a constant.
S1.3, establishing a water inlet model of the water tank, which specifically comprises the following steps: cVThe control valve is a stop valve for the flow coefficient of the valve, and the relative flow coefficient is in direct proportion to the opening degree of the valve, namely, the following conditions are met:
Figure BDA0003017961420000084
wherein epsilon is a coefficient, and D is the inner diameter of the water inlet pipeline;
from conservation of mass:
Figure BDA0003017961420000085
wherein SWater tankIs the bottom area of the water tank C, Qout(t) is the outflow rate of the cistern, since QinIt is necessary to flow through a section of pipe network from the control valve to the top of the tank, so
Figure BDA0003017961420000086
However, since the difference between the two is small, the processing is simplified
Figure BDA0003017961420000087
Substitution
Figure BDA0003017961420000088
Therefore, the simplified formulas (2), (3), (4) and (5) are simplified to obtain the formula (6):
Figure BDA0003017961420000089
the formulas (1) and (6) form a water inlet model of the water tank, and the water inlet model is along with the water outlet flow QoutAnd (t) changing the liquid level of the water tank according to the water inlet model, so that the change condition of the liquid level of the water tank of the second pump room in the community for 24 hours can be simulated. S2, collecting data and checking a water tank water inlet model, and the method specifically comprises the following steps:
s2.1 Water consumption at the time of collecting the cell Qout(t), the inner diameter D of the water inlet pipeline of the community. Bottom area S of tank CWater tankArea S of diaphragm M, height H of water tank level when valve is closedmaxHeight H of water tank level when valve is fully openedminCalculating the sectional area S of the water inlet pipeline by using the elastic coefficient K of the spring K1inAnd head loss Δ H at the valve;
s2.2 Collection district Inlet pipe inflow QinAnd (t) calculating the water inflow flow rate.
S2.3, the water inlet model of the water tank is checked, and the method specifically comprises the following steps:
s2.31 sets the value range of the coefficient ∈ min, ∈ max, the initial value ∈ 0 ∈ min, and the step length Δ ∈ max- ∈ min)/100, and sets the optimization target to make the actually measured liquid level value h (t) most match the analog value hs (t), that is, Jh ∑ h (t) -hs (t)) | minimize;
s2.32 respectively substituting epsilon 0, epsilon 1, … … epsilon 100 and related data into an equation (6), carrying out 24-hour simulation, recording simulation data, namely hs (t) at each moment, and calculating to obtain objective function values Jh0 and Jh1 … … Jh 100;
s2.33 finds the minimum value Jhmin from Jh0 and Jh1 … … Jh100, the epsilon value at the moment is the fixed value of the model rate, and the water inlet model check is completed.
S3 determining the optimization variables, constraints, and optimization objectives, specifically including:
s3.1 when the diaphragm of the self-operated liquid level control valve is selected, the adjustable variables are only the spring variable x (t) and the height H of the liquid level of the water tank corresponding to the closing and opening valvesmaxAnd Hmin. Both quantities are related to the spring constant K. The optimized variable is the spring rate K.
S3.2 constraint Condition：Qin(t) not greater than MaxQ of maximum inflow of pipe networkinThe water level in the water tank can not be lower than the lowest level HminNamely:
0≤Qin(t)≤MaxQin (7)
Hmin≤h(t)≤Hmax (8)
s3.3 optimization goal: in the face of water interference, water is supplemented to the water tank as soon as possible, and sudden change of water inlet (water hammer caused by too fast valve opening or closing) is reduced.
Namely Max (| Q) corresponding to K valuein' (t) |) and Max (H-H (t)) are both small, but these two goals are contradictory: max (| Q)inThe smaller the' (t) |), the larger the Max (H-H (t)). Adding weight, and constructing an optimization objective function:
J＝a*Max(|Qin′(t)|)+b*Max(H-h(t)) (9)
where a, b are weights such that J is minimized.
S4 determines the optimization parameters by calculation, which specifically includes:
s4.1, determining the range, initial value and step length of the optimized variable K, and according to the maximum displacement of the valve rod, the lowest liquid level H of the water tank during water inlet and outletWater tank minDetermining Kmax(ii) a Minimum level H based on maximum displacement of valve stem 1/2Water tank minDetermining Kmin(ii) a Then there is Kmin≤K≤KmaxSetting K as initial Kmin, setting step length delta K as (Kmax-Kmin)/100;
s4.2, calculating an optimization target subentry index value under the current K value, which specifically comprises the following steps:
s4.21 takes 1 day with the maximum daily water consumption of the community as the basis, and assumes that 0 is the initial state, H (0) is Hmax，Qout(0) When K is 0 and K is Kmin, h' (0) is obtained by substituting formula (6);
s4.22h' (0) is multiplied by the time interval delta t to obtain h (1); then, Q is obtained from the formula (5)in(0)；
S4.23 relating h (1) to Q at the next momentout(1) The value is substituted for formula (6);
s4.24 repeating steps S4.21, S4.22 and S4.23, calculating all h' (t), h (t) and Q under the same K valuein(t)；
S4.25 calculating the corresponding | Q at each moment under the current K valuein' (t) | and H-H (t), and records Max (| Q)in' (t) |) and Max (H-H (t)).
S4.3 if K<KmaxIf yes, returning to step S4.2;
s4.4, checking constraint conditions, drawing an optimization target subentry index curve under different K values, and aiming at Q under different K valuesin(t) h (t), checking whether the constraint condition is met, if not, rejecting the K value, and respectively drawing Max (| Q) corresponding to different K valuesin' (t) |), Max (H-H (t));
s4.5 calculating the value of the objective function J under different K values
The weights a, b are determined because of Max (| Q)in' (t) |), Max (H-H (t)) are not of the same order, so first Max (| Q)in' (t) |, Max (H-H (t)) are respectively standardized; considering that the importance of the two is basically the same, the weight factor of the processed data is distributed to 0.5, and objective function values J under different K values are calculated; and searching Jmin, wherein the K value corresponding to the Jmin is the optimal solution.
Example 2
The method of the present invention is described with a certain cell as an example.
The diameter of the water inlet pipe of the community is
Figure BDA0003017961420000112
The average water consumption of the residential quarter 11-12 months and days is 518m3Highest daily water consumption of 699m3According to the design requirement specification of the water tank, 20-25% of the daily maximum water consumption needs to be met, so that the design is carried out according to 700m3Considering that the effective volume of the water tank is 140m3。
Step 1: establishing water inlet model of water tank
Combining the data of the community, the pipe diameter of the water inlet pipe of the community is known to be
Figure BDA0003017961420000113
Rho is the density of water, the gravity acceleration g is 10N/kg, the inner diameter d of the water tank C is 9.5M, the flow speed is 1.4M/S, and the area S of the membrane M isWhen the square meter is 0.04, the water head loss Δ H at the valve can be obtained from the formula (10).
Δ H ═ pipe network water pressure-tank height (10)
The water pressure of the pipe network in the formula is approximately equal to 27 m; the height of the water tank is approximately equal to 3 m;
substituting to solve the problem to obtain delta H which is approximately equal to 24 m. The simplified objective function can be obtained by substituting the above quantities into formula (10):
Figure BDA0003017961420000111
unit: m/s.
Step 2: collecting relevant data and checking water tank water inlet model
Collecting and obtaining the maximum daily water consumption Q of the communityoutAnd the inner diameter D of the district water inlet pipeline is 100 mm. Bottom area S of water tank CWater tankSquare meter 70, square meter 0.4, water tank level H when valve is closedmax2m, the height H of the water tank when the valve is fully openedmin1.2m, 50000N/m elastic coefficient K of spring K1, and calculating sectional area S of water inlet pipelineinAnd head loss Δ H at the valve.
The highest daily water consumption Q of the communityoutThe data and curves are shown in table 1 and fig. 3.
Figure BDA0003017961420000121
TABLE 1
2) Collecting to obtain the flow rate of the water inlet pipe network of the community at 1.4m/s and the rest time periods QinThe flow is shown in table 2, and the water inlet flow of the water tank can reach 39.6m per hour when the valve is in a full-open state3。
Figure BDA0003017961420000122
Figure BDA0003017961420000131
TABLE 2
3) And checking the water inlet model.
Setting the value range of coefficient epsilon to be 0.5,2.5]Initial value ε 0 is equal to 0.5, and step size
Figure BDA0003017961420000132
Figure BDA0003017961420000133
Setting an optimization target that the actually measured liquid level value is most consistent with the analog value, namely Jh ═ h (t) -hs (t) and (t) are minimized;
respectively substituting epsilon 0, epsilon 1, … … epsilon 100 and related data into an equation (11), carrying out 24-hour simulation, recording simulation data (hs (t) at each moment), and calculating to obtain target function values Jh0 and Jh1 … … Jh 100;
and finding the minimum value Jhmin from Jh0 and Jh1 … … Jh100, wherein the epsilon is 1.02 which is the model rate fixed value, and completing the water inlet model check. The analog value of the model liquid level after checking is compared with the measured value, as shown in FIG. 4.
And 3, step 3: determining optimization variables, constraints and optimization objectives
1) When the self-operated liquid level control valve is selected, the adjustable variable is only the spring variable x (t), and the liquid level height H of the water tank corresponding to the closed and opened valvemaxAnd Hmin. These quantities are related to the spring constant K. The optimized variable is the spring rate K.
2) Constraint conditions are as follows: qin(t) not greater than MaxQ of maximum inflow of pipe networkinThe water level in the water tank can not be lower than the lowest level H for a long timeWater tank minNamely:
0≤Qin≤0.011
0.5≤h(t)≤2
3) optimizing the target: in the face of water interference, water is supplemented to the water tank as soon as possible, and sudden change of water inlet (water hammer caused by too fast valve opening or closing) is reduced.
Namely Max (| Q) corresponding to K valuein' (t) |) and Max (H-H (t)) are both smaller, but these two goals are contradictory: max (| Q)in′(t) |) the smaller, the larger the Max (H-H (t)). Adding weight, and constructing an optimization objective function:
J＝a*Max(|Qin′(t)|)+b*Max(H-h(t))#(9)
where a, b are weights such that J is minimized.
And 4, step 4: calculating optimization and determining optimization parameters
(1) Determining the range, initial value and step length of the optimized variable K
According to the maximum displacement of the valve rod, the lowest liquid level H of the water tank when water is fed and dischargedWater tank minK can be determined at 0.5 mmax81000N/m; maximum displacement, minimum H, of valve stem 1/2Water tank minK can be determined at 1 metermin44000N/m. Then K is equal to or more than 44000 and equal to or less than 81000N/m, K is set as initial value Kmin, step length
Figure BDA0003017961420000141
(2) Calculating the optimized target subentry index value under the current K value
The method is based on 1 day with the maximum daily water consumption of a community, and assumes that 0 is an initial state, and H (0) is Hmax，Qout(0) When K is 0 and K is Kmin, h' (0) is obtained by substituting equation (11);
h (1) can be obtained by multiplying h' (0) by the time interval delta t; then, Q can be obtained from the formula (5)in(0)；
③ h (1) and Q at the next momentout(1) The value continues to be substituted into formula (11);
fourthly, because of the water outlet flow Qout(t) changing in real time, repeating the third step, and calculating all h' (t), h (t) and Q under the same K valuein(t)；
Calculating corresponding | Q at each moment under current K valuein' (t) | and H-H (t), and records Max (| Q)in' (t) |) and Max (H-H (t)).
(3) If K<KmaxIf K is K + Δ K, return to step 4- (2)
(4) Checking constraint conditions, and drawing an optimized target subentry index curve under different K values
And (5) checking whether constraint conditions are met or not aiming at Qin (t) and h (t) under different K values, and if not, rejecting the K value.
Respectively drawing corresponding Max (| Q) under different K valuesin' (t) |), Max (H-H (t)), and a subentry index curve. As shown in fig. 5 and 6.
(5) And calculating the objective function value J at different K values.
The weights a, b are determined because of Max (| Q)in' (t) |), Max (H-H (t)) are not of the same order, so first Max (| Q)in' (t) |, Max (H-H (t)) are respectively standardized; considering that the importance of the two is basically the same, the weight factor of the processed data is distributed to 0.5, and objective function values J under different K values are calculated; as shown in fig. 7, the image is plotted, and Jmin-0.125 is found, where K corresponding to Jmin is the optimal solution, and K is 58800N/m.
And replacing the existing spring, and comparing the change rate of the water level of the water tank and the change rate of the water inlet speed under the condition that K is 58800N/m after replacement and K is 50000N/m before replacement, as shown in figures 8 and 9. It can be seen that when K is 58800N/m, although the water tank liquid level can be a little lower when the peak period water is used, the sudden change of the water inlet is reduced on the premise of meeting the water consumption of residents, the service life of the valve is prolonged, and therefore the effect is better.
The specific embodiments described herein are merely illustrative of the spirit of the invention. Various modifications or additions may be made to the described embodiments or alternatives may be employed by those skilled in the art without departing from the spirit or ambit of the invention as defined in the appended claims.
Although the terms tank inlet model, steady state and transient hydraulics are used more herein, the possibility of using other terms is not excluded. These terms are used merely to more conveniently describe and explain the nature of the present invention; they are to be construed as being without limitation to any additional limitations that may be imposed by the spirit of the present invention.

Claims (6)

1. A parameter optimization method for an intelligent two-water-supply device for a waterproof hammer is characterized by comprising the following steps:
s1, establishing a water inlet model of the water tank, which specifically comprises the following steps:
s1.1, water inlet model parameters of the water tank are set, and the method specifically comprises the following steps: setting the deformation x of the spring K1, the height H of the water tank liquid level when water is in and out, and the height H of the water tank liquid level when the valve is closedmaxThe pressure of the cavity A is
FA＝PS＝ρghS
P is the pressure acting on the cavity A, S is the area of the diaphragm M, rho is the density of water, g is 9.8N/kg, and the elastic force of the spring K1 is
FK1＝Kx
Wherein K is the elastic coefficient of the spring K1
Balance equation of membrane in valve body:
pgh(t)S＝Kx(t) (1)
where t is time, h (t) is the liquid level at time t, x (t) is the amount of spring deformation,
the difference of current liquid level of water tank and highest liquid level is influencing the opening degree size of valve, and the difference is big more, and the aperture is big more:
Figure FDA0003593840080000011
wherein
Figure FDA0003593840080000012
h(t)≤Hmax；
S1.2, calculating the representation of the adjusting action of the control valve in steady-state and transient hydraulic power;
s1.3, establishing a water inlet model of a water tank;
s2, collecting data and checking the water inlet model of the water tank;
s3, determining the optimization variables, constraints, and optimization objectives, specifically including:
s3.1 when the diaphragm of the self-operated liquid level control valve is selected, the adjustable variable is the spring deformation x (t), the closed valve and the corresponding liquid level height H of the water tank when the valve is fully openedmaxAnd HminThe two quantities are related to the spring elastic coefficient K, and the optimized variable is the spring elastic coefficient K;
s3.2 restraint strapA piece: qin(t) not greater than MaxQ of maximum inflow of pipe networkinThe water level in the water tank can not be lower than the lowest level HminNamely:
0≤Qin(t)≤MaxQin (7)
Hmin≤h(t)≤Hmax (8)
s3.3 optimization goal: max (| Q) corresponding to K valueinReducing the values of' (t) | and Max (H-H (t)), adding weight to construct an optimized objective function:
J＝a*Max(|Qin′(t)|)+b*Max(H-h(t)) (9)
wherein a, b are weights such that J is minimized;
s4 determines the optimization parameters by calculation, which specifically includes:
s4.1, determining the range, initial value and step length of the optimized variable K, and according to the maximum displacement of the valve rod, the lowest liquid level H of the water tank during water inlet and outletWater tank minDetermining Kmax(ii) a Minimum level H according to maximum displacement of valve stem 1/2Water tank minDetermining Kmin(ii) a Then there is Kmin≤K≤KmaxSetting K as initial Kmin, setting step length delta K as (Kmax-Kmin)/100;
s4.2, calculating an optimization target subentry index value under the current K value;
s4.3 if K<KmaxIf yes, returning to step S4.2;
s4.4, checking constraint conditions, drawing an optimization target subentry index curve under different K values, and aiming at Q under different K valuesin(t) h (t), checking whether the constraint condition is met, if not, rejecting the K value, and respectively drawing the corresponding Max (| Q) values under different K valuesin' (t) |), Max (H-H (t));
s4.5, calculating objective function values J under different K values, determining weights a and b, and firstly, calculating Max (| Q)in' (t) |, Max (H-H (t)) are respectively standardized; distributing 0.5 to the processed data weight factor, and calculating objective function values J under different K values; and searching Jmin, wherein the K value corresponding to the Jmin is the optimal solution.
2. The parameter optimization method for the intelligent secondary water supply device for the waterproof hammer according to claim 1, wherein the step S1.2 specifically comprises the following steps: the control valve modulation is expressed in both steady-state and transient hydraulic calculations
Figure FDA0003593840080000031
Wherein QinΔ H and SinRespectively the flow through the valve, the head loss at the valve and the cross-sectional area of the water inlet pipe, Δ H being the pipe network water pressure-water tank height, Δ H simplification being regarded as a constant, CVIs the flow coefficient of the valve.
3. The parameter optimization method for the intelligent secondary water supply device for the water hammer prevention according to claim 2, wherein the step S1.3 specifically comprises the following steps: cVThe control valve is a stop valve for the flow coefficient of the valve, and the relative flow coefficient is in direct proportion to the opening degree of the valve, namely, the following conditions are met:
Figure FDA0003593840080000032
wherein epsilon is a coefficient, and D is the inner diameter of the water inlet pipeline;
from conservation of mass:
Figure FDA0003593840080000033
wherein SWater tankIs the bottom area of the tank C, Qout(t) is the effluent flow rate due to QinIt is necessary to flow through a section of pipe network from the control valve to the top of the tank, so
Figure FDA0003593840080000034
However, since the difference between the two is small, the processing is simplified
Figure FDA0003593840080000035
Substitution
Figure FDA0003593840080000036
Therefore, the simplified formulas (2), (3), (4) and (5) are simplified to obtain the formula (6):
Figure FDA0003593840080000037
the formulas (1) and (6) form a water tank water inlet model, and the water inlet model is along with the water outlet flow QoutAnd (t) changing the liquid level of the water tank according to the water inlet model.
4. The parameter optimization method for the intelligent secondary water supply device for the waterproof hammer according to claim 1, wherein the step S2 of checking the water inlet model of the water tank comprises the following specific steps:
s2.1 collecting the effluent flow Qout(t) inner diameter D of water inlet pipeline of community and bottom area S of water tank CWater tankArea S of diaphragm M, height H of water level of water tank when valve is closedmaxHeight H of water tank liquid level when valve is fully openedminCalculating the sectional area S of the water inlet pipeline by using the elastic coefficient K of the spring K1inAnd head loss Δ H at the valve;
s2.2 Collection district Inlet pipe inflow Qin(t), calculating the inflow water flow rate;
s2.3, checking the water inlet model of the water tank.
5. The parameter optimization method for the intelligent secondary water supply device for the water hammer prevention according to claim 4, wherein the step S2.3 specifically comprises the following steps:
s2.31 sets the value range of the coefficient ∈ min, ∈ max, the initial value ∈ 0 ∈ min, and the step length Δ ∈ max- ∈ min)/100, and sets the optimization target to make the actually measured liquid level value h (t) most match the analog value hs (t), that is, Jh ∑ h (t) -hs (t)) | minimize;
s2.32 respectively substituting epsilon 0, epsilon 1, … … epsilon 100 and related data into an equation (6), carrying out 24-hour simulation, recording simulation data, namely hs (t) at each moment, and calculating to obtain objective function values Jh0 and Jh1 … … Jh 100; s2.33 finds the minimum value Jhmin from Jh0 and Jh1 … … Jh100, the epsilon value at the moment is the fixed value of the model rate, and the water inlet model check is completed.
6. The parameter optimization method for the intelligent secondary water supply device for the water hammer prevention according to claim 1, wherein the step S4.2 specifically comprises the following steps:
s4.21 takes 1 day with the maximum daily water consumption of the community as the basis, and assumes that 0 is the initial state, H (0) is Hmax，Qout(0) When K is 0 and K is Kmin, h' (0) is obtained by substituting formula (6);
s4.22h' (0) is multiplied by the time interval delta t to obtain h (1); then, Q is obtained from the formula (5)in(0)；
S4.23 comparing h (1) with Q at the next momentout(1) The value is substituted for formula (6);
s4.24 Steps S4.21, S4.22, S4.23 are repeated, and all of h' (t), h (t), and Q at the same K value are calculatedin(t)；
S4.25 calculating the corresponding | Q at each moment under the current K valuein' (t) | and H-H (t), and records Max (| Q)in' (t) |) and Max (H-H (t)).
CN202110394232.XA 2021-04-13 2021-04-13 Parameter optimization method for intelligent secondary water supply device for water hammer prevention Active CN113177301B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110394232.XA CN113177301B (en) 2021-04-13 2021-04-13 Parameter optimization method for intelligent secondary water supply device for water hammer prevention

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110394232.XA CN113177301B (en) 2021-04-13 2021-04-13 Parameter optimization method for intelligent secondary water supply device for water hammer prevention

Publications (2)

Publication Number Publication Date
CN113177301A CN113177301A (en) 2021-07-27
CN113177301B true CN113177301B (en) 2022-07-08

Family

ID=76923539

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110394232.XA Active CN113177301B (en) 2021-04-13 2021-04-13 Parameter optimization method for intelligent secondary water supply device for water hammer prevention

Country Status (1)

Country Link
CN (1) CN113177301B (en)

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100625071B1 (en) * 2002-12-19 2006-09-20 가부시키가이샤 후지킨 Method for closing fluid passage, and water hammerless closing device
CN101725740B (en) * 2009-12-31 2012-03-28 广东联塑科技实业有限公司 Impeller type water hammer-resistant mute check valve
CN108506628B (en) * 2018-04-18 2020-10-30 中国天辰工程有限公司 Safety control method for synchronous multistage valve positioning large-diameter pipeline system waterproof hammer
CN111932080A (en) * 2020-07-09 2020-11-13 上海威派格智慧水务股份有限公司 Early warning protection system and method applied to water service pipe network
CN111985166A (en) * 2020-07-31 2020-11-24 河海大学 Pipeline hydraulic transient simulation method and storage medium with implicit consideration of dynamic friction resistance
CN111853290B (en) * 2020-08-11 2021-12-10 中国水利水电科学研究院 Variable-flow-channel water hammer prevention air valve and design method
CN111985060B (en) * 2020-08-11 2023-09-29 中国水利水电科学研究院 Caliber calculation method for air valve waterproof hammer of pipeline water delivery system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Water hammer analysis when switching of parallel pumps based on contra-motion check valve;Yang Z 等;《Annals of Nuclear Energy》;20201231;第1-12页 *

Also Published As

Publication number Publication date
CN113177301A (en) 2021-07-27

Similar Documents

Publication Publication Date Title
CN111580382B (en) Unit-level heat supply adjusting method and system based on artificial intelligence
CN205314168U (en) Pump water supply installation in box no negative pressure silence pipe of intelligent regulation water tank water demand
CN113177301B (en) Parameter optimization method for intelligent secondary water supply device for water hammer prevention
CN208121879U (en) Automatically adjust the water tank water intake control system of flow of inlet water
CN108427459B (en) Automatic control method for piston type flow regulating valve of multi-water-source annular water supply pipe network
CN105068582B (en) DCS-based high-level water collection cooling tower water level maintaining control method
CN205013067U (en) Heat supply of backpressure machine, energy storage system
CN112064728A (en) Self-adaptive water supply system
CN102734636A (en) Static pressure continuous water drainage device
CN113836785B (en) Municipal area intelligent water supply system and artificial intelligent control optimization method thereof
CN115268518A (en) Secondary water supply inflow and main pipe network balancing method and system
CN2687458Y (en) Self-operated flow control valve
CN209227676U (en) A kind of intelligence is avoided the peak hour secondary water supply system
CN208474577U (en) Novel water supply decompressor
CN111754362A (en) Method for evaluating emptying capacity of multistage emptying dam
CN208718010U (en) Multifunction pressure tower
CN219827905U (en) Pressure reducing valve with built-in water pressure self-control pipeline
CN2225592Y (en) Steam trap regulating valve
CN214540536U (en) Self-adaptive water inlet control device for water tank
CN212375955U (en) Self-adaptive water supply system
CN202024071U (en) Water discharging device
CN109813399B (en) Self-adaptive constant-liquid-level flow stabilizer and operation method thereof
CN115268516A (en) Regulation and control method and system for improving utilization rate of secondary water supply tank
CN217153223U (en) Remote control ball cock with delay function of opening and close
CN217108767U (en) High-level water tank water hammer eliminating system

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 310016 second floor, building 15, No.260, No.6 street, economic and Technological Development Zone, Jianggan District, Hangzhou City, Zhejiang Province

Applicant after: Fengtai water technology (Hangzhou) Co.,Ltd.

Address before: 310016 second floor, building 15, No.260, No.6 street, economic and Technological Development Zone, Jianggan District, Hangzhou City, Zhejiang Province

Applicant before: HANGZHOU ZHONGZI FENGTAI ENVIRONMENT TECHNOLOGY Co.,Ltd.

GR01 Patent grant
GR01 Patent grant