A water supply pipeline management system and method

Abstract

The present invention provides a water supply pipeline management system and method, comprising a leakage management module, a pipeline life assessment module, an operation and maintenance manpower estimation module, a water supply demand prediction module, and a water supply resource configuration module; the leakage management module is used to find out the section where the leaking pipeline is located and to give a reminder based on an electronic map, the pipeline life assessment module is used to estimate the remaining life of the pipeline, the operation and maintenance manpower estimation module is used to estimate the manpower required this year to ensure the normal operation of the built pipeline, the water supply demand prediction module is used to estimate the monthly water supply demand data, and the water supply resource configuration module is used to dynamically adjust the allocation of water supply resources. The system and method can improve the intelligence level of water supply pipeline management and optimize the water supply service quality of the water supply pipeline.

Description

Technical Field

[0001] The present invention relates to the field of software informatization, and in particular to a water supply pipeline management system and method. Background Art

[0002] Traditional water supply pipeline management system software is mainly based on basic information query, and there is still much room for improvement in the level of intelligence. It is necessary to introduce more refined and intelligent technologies and methods to optimize the water supply service quality of water supply pipelines. Summary of the invention

[0003] In order to improve the intelligent level of water supply pipeline management, the present invention designs a water supply pipeline management system and method.

[0004] The technical solution adopted by the present invention is a water supply pipeline management system: It includes leakage management module, pipeline life assessment module, operation and maintenance manpower estimation module, water supply demand forecasting module and water supply resource allocation module.

[0005] The leakage management module is used to find the section where the leaking pipe is located and to remind based on the electronic map, specifically in the following ways: A1, assume that the branches and intersections in the pipeline network are nodes of the pipeline, and the section between two adjacent nodes is a section of pipeline. Number each section of pipeline in the form of the longitude and latitude coordinates of the upstream node of this section of pipeline + the longitude and latitude coordinates of the downstream node, and the unit of the coordinate value is degree; for example, the coordinates of the upstream node of a section of pipeline are (30.54321N, 100.12345E), and the coordinates of the downstream node are (30.54320N, 100.12347E), then the number of this section of pipeline is: 30.54321N100.12345E30.54320N100.12347E.

[0006] A2, let the length of the pipe between the nodes be L, the diameter be D, and the friction coefficient be f ; Then the theoretical pressure loss H flowing through each section of the pipeline f = f ×L÷D×( v 2 )÷(2×g), where v is the flow velocity of the fluid in the pipe, and g is the acceleration due to gravity.

[0007] A3, the actual pressure loss H flowing through each section of the pipeline e =P in -P out , where P in The water pressure at the upstream node of this section of pipeline, P out The water pressure at the node downstream of this pipe segment.

[0008] A4, when H e >H f When the pipeline is detected, it is determined that there is a water leakage problem in this section of the pipeline; on the electronic map, pipelines in normal state are indicated in green, and pipelines with water leakage problems are indicated in red. When there is a red pipeline, a prompt sound is issued.

[0009] The pipeline life assessment module is used to estimate the remaining life of the pipeline in the following manner: B1, create a database of historical data of pipeline life t, average temperature of pipeline external environment T, and water pressure P, where the unit of t is unified as year, the unit of T is unified as degree Celsius, and the unit of P is unified as Pascal.

[0010] B2, construct pipeline life distribution function: ;in: ; β0 is the basic life parameter, which is determined by material, thickness and diameter. θ1 and θ2 are the influence coefficients of the average temperature of the external environment and water pressure on the pipeline life respectively.

[0011] B3, use maximum likelihood to estimate the values ​​of β0, θ1, θ2: ①Define the survival function: ; Probability density function: ; ② For the same sample data of β0 (t1, t1, …, t n ), the likelihood function is: ; ③ Take the log-likelihood function: ; Expand the log-likelihood function: ; ④ Randomly initialize the values ​​of β0, θ1, and θ2, and calculate their partial derivatives ∂β0, ∂θ1, and ∂θ2 respectively. 0 new =β 0 old +α×∂β 0 old ,θ 1 new =θ 1 old +α×∂θ 1 old ,θ 2 new =θ 2 old +α×∂θ 2 old The law is iterated, where β 0 new ,θ 1 new ,θ 2 new represents the newly generated value; β 0 old ,θ 1 old ,θ 2 oldrepresents the old value; α is the learning rate, which controls the step size of each update; until the value of the log-likelihood function reaches the maximum, the optimal values ​​of β0, θ1, and θ2 are obtained.

[0012] B4, evaluate the remaining life of the pipeline with the same β0: The β corresponding to the pipeline is calculated based on the average temperature and water pressure of the pipeline's external environment. The pipeline life distribution function is exponentially distributed, and the expected remaining life in years is approximately equal to the value of β. The value of β is the estimated remaining life of the pipeline.

[0013] The operation and maintenance manpower estimation module is used to estimate the manpower required this year to ensure the normal operation of the built pipelines. The specific method is as follows: Obtain the leakage status and remaining life of each pipe section, give priority to dealing with leaking pipes, and then replace pipes with nearly the same lifespan; sort pipes with a remaining lifespan less than (β+1) from least to most, and replace pipes with the lowest remaining lifespan first.

[0014] Suppose that there were m pipeline leaks in the previous year, and the average time spent on repairing leaks was K man-hours / time; suppose that it takes J0 man-hours / segment to replace pipelines on average, and the number of pipelines with a remaining life less than (β+1) at the beginning of the year is R; workers work h hours a day, then the number of pipeline operation and maintenance personnel to be prepared this year is N e =ceil [(m×K+R×J0)÷(365×h) ×(1+γ)], where γ is the proportion of reserve personnel and ceil() is the upward rounding function.

[0015] In actual implementation, γ is usually set to 10% to 20%.

[0016] The water supply demand forecasting module is used to estimate the monthly water supply demand data in the following way: C1, create a historical water supply database, including months and corresponding water supply volumes.

[0017] C2, create a prediction function according to the season, set the month x ,when x ≥12, converted to x =0; Suppose the spring prediction function is y spring = a spring × x + b spring , x =0, 1, 2; The summer prediction function is y summer = a summer × x + b summer , x =3, 4, 5; The autumn forecast function is y autumn = a autumn × x + b autumn , x =6,7,8; The winter prediction function is y winter = a winter × x + b winter , x =9, 10, 11; set up x 0. x 1. x 2. x 3. x 4. x 5. x 6. x 7. x 8. x 9. x 10 , x 11 The months they represent are 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, and 11, and the corresponding water supply is y 0. y 1. y 2. y 3. y 4. y 5. y 6. y 7. y 8. y 9. y 10 , y 11 ; a spring ={w spring ×[( x 0× y 0)+ ( x 1× y 1)+ ( x 2× y 2)]-( x 0+ x 1+ x 2)×( y 0+ y 1+ y 2)}÷[w spring ×( x 0 2 + x 12 + x 2 2 )-( x 0+ x 1+ x 2) 2 ]; w spring The number of sampling points for spring water supply; a summer ={w summer ×[( x 3× y 3)+ ( x 4× y 4)+ ( x 5× y 5)]-( x 3+ x 4+ x 5)×( y 3+ y 4+ y 5)}÷[w summer ×( x 3 2 + x 4 2 + x 5 2 )-( x 3+ x 4+ x 5) 2 ]; w summer The number of sampling points for summer water supply; a autumn ={w autumn ×[( x 6× y 6)+ ( x 7× y 7)+ ( x 8× y 8)]-( x 6+ x 7+ x 8)×( y 6+ y 7+ y 8)}÷[w autumn ×( x 6 2 + x 7 2 + x 8 2 )-( x 6+ x 7+ x 8) 2 ]; w autumn number of sampling points for fall water supply; a winter ={w winter ×[( x 9× y 9)+ ( x 10 × y 10 )+ ( x 11 × y 11 )]-( x 9+ x 10 + x 11 )×( y 9+ y 10 + y 11 )}÷[w winter ×( x 9 2 + x 10 2 + x 11 2 )-( x 9+ x 10 + x 11 ) 2 ]; w winter The number of sampling points for winter water supply; b spring =[( y 0+ y 1+ y 2)-a spring ×( x 0+ x 1+ x 2)]÷w spring ; b summer =[( y 3+ y 4+ y 5)-a summer ×( x 3+ x 4+ x 5)]÷w summer ; b autumn =[( y 6+ y 7+ y 8)-a autumn ×( x 6+ x 7+x 8)]÷w autumn ; b winter =[( y 9+ y 10 + y 11 )-a winter ×( x 9+ x 10 + x 11 )]÷w winter .

[0018] C3, yes x By assigning a month, you can get the demand corresponding to the predicted month.

[0019] The water supply resource configuration module is used to dynamically adjust the allocation of water supply resources, specifically in the following ways: D1, open multiple water supply points in the water supply pipeline network until the sum of the monthly opened water supply points in the water supply pipeline network is greater than the monthly demand predicted by C3, otherwise continue to increase the number of opened water supply points; D2, there are multiple demand points in the water supply pipeline network. Each water supply point stops supplying water after reaching its own capacity. At the same time, the number of demand points is also changing dynamically. Suppose there are U water supply points and V demand points at a certain moment. Suppose from the water supply point i To the demand point j The water distribution cost is C ij , set up water supply points i Assign to demand points j The amount of water is D ij , demand point j The quantity demanded is Q j ; Set the prerequisites: , and adjust the allocation of water supply resources based on this premise; Set adjustment goals: ; Min() is the minimization function, which is used as the goal to dynamically adjust the allocation of water supply resources.

[0020] A water supply pipeline management method: Step 1: Find the section where the leaking pipe is located and give a reminder based on the electronic map. The specific method is as follows: A1, assume that the branches and intersections in the pipeline network are nodes of the pipeline, and the section between two adjacent nodes is a section of pipeline. Each section of pipeline is numbered in the form of the longitude and latitude coordinates of the upstream node of this section of pipeline + the longitude and latitude coordinates of the downstream node, and the unit of the coordinate value is degree; for example, the coordinates of the upstream node of a section of pipeline are (30.54321N, 100.12345E), and the coordinates of the downstream node are (30.54320N, 100.12347E), then the number of this section of pipeline is: 30.54321N100.12345E30.54320N100.12347E.

[0021] A2, let the length of the pipe between the nodes be L, the diameter be D, and the friction coefficient be f ; Then the theoretical pressure loss H flowing through each section of the pipeline f = f ×L÷D×( v 2 )÷(2×g), where v is the flow velocity of the fluid in the pipe, and g is the acceleration due to gravity.

[0022] A3, the actual pressure loss H flowing through each section of the pipeline e =P in -P out , where P in The water pressure at the upstream node of this section of pipeline, P out The water pressure at the node downstream of this pipe segment.

[0023] A4, when H e >H f When the pipeline is detected, it is determined that there is a water leakage problem in this section of the pipeline; on the electronic map, pipelines in normal state are indicated in green, and pipelines with water leakage problems are indicated in red. When there is a red pipeline, a prompt sound is issued.

[0024] Step 2: Estimate the remaining life of the pipeline by: B1, create a database of historical data of pipeline life t, average temperature of pipeline external environment T, and water pressure P, where the unit of t is unified as year, the unit of T is unified as degree Celsius, and the unit of P is unified as Pascal.

[0025] B2, construct pipeline life distribution function: ;in: ; β0 is the basic life parameter, which is determined by material, thickness and diameter. θ1 and θ2 are the influence coefficients of the average temperature of the external environment and water pressure on the pipeline life respectively.

[0026] B3, use maximum likelihood to estimate the values ​​of β0, θ1, θ2: ①Define the survival function: ; Probability density function: ; ② For the same sample data of β0 (t1, t1, …, t n ), the likelihood function is: ; ③ Take the log-likelihood function: ; Expand the log-likelihood function: ; ④ Randomly initialize the values ​​of β0, θ1, and θ2, and calculate their partial derivatives ∂β0, ∂θ1, and ∂θ2 respectively. 0 new =β 0 old +α×∂β 0 old ,θ 1 new =θ 1 old +α×∂θ 1 old ,θ 2 new =θ 2 old +α×∂θ 2 old The law is iterated, where β 0 new ,θ 1 new ,θ 2 new represents the newly generated value; β 0 old ,θ 1 old ,θ 2 old represents the old value; α is the learning rate, which controls the step size of each update; until the value of the log-likelihood function reaches the maximum, the optimal values ​​of β0, θ1, and θ2 are obtained.

[0027] B4, evaluate the remaining life of the pipeline with the same β0: The β corresponding to the pipeline is calculated based on the average temperature and water pressure of the pipeline's external environment. The pipeline life distribution function is exponentially distributed, and the expected remaining life in years is approximately equal to the value of β. The value of β is the estimated remaining life of the pipeline.

[0028] Step 3: Estimate the manpower required this year to ensure the normal operation of the built pipelines. The specific method is as follows: Obtain the leakage status and remaining life of each pipe section, give priority to dealing with leaking pipes, and then replace pipes with nearly the same lifespan; sort pipes with a remaining lifespan less than (β+1) from least to most according to their remaining lifespan, and give priority to replacing pipes with the lowest remaining lifespan.

[0029] Suppose that there were m pipeline leaks in the previous year, and the average time spent on repairing leaks was K man-hours / time; suppose that it takes J0 man-hours / segment to replace pipelines on average, and the number of pipelines with a remaining life less than (β+1) at the beginning of the year is R; workers work h hours a day, then the number of pipeline operation and maintenance personnel to be prepared this year is N e =ceil [(m×K+R×J0)÷(365×h) ×(1+γ)], where γ is the proportion of reserve personnel and ceil() is the upward rounding function.

[0030] In actual implementation, γ is usually set to 10% to 20%.

[0031] Step 4: Estimate monthly water supply demand data by: C1, create a historical water supply database, including months and corresponding water supply volumes.

[0032] C2, create a prediction function according to the season, set the month x ,when x ≥12, converted to x =0; Suppose the spring prediction function is y spring = a spring × x + b spring , x =0, 1, 2; The summer prediction function is y summer = a summer × x + b summer , x =3, 4, 5; The autumn forecast function is y autumn = a autumn × x + b autumn , x =6,7,8; The winter prediction function is y winter = a winter × x + b winter , x =9, 10, 11; set up x 0. x 1. x 2. x 3. x 4. x 5. x 6. x 7. x 8.x 9. x 10 , x 11 The months they represent are 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, and 11, and the corresponding water supply is y 0. y 1. y 2. y 3. y 4. y 5. y 6. y 7. y 8. y 9. y 10 , y 11 ; a spring ={w spring ×[( x 0× y 0)+ ( x 1× y 1)+ ( x 2× y 2)]-( x 0+ x 1+ x 2)×( y 0+ y 1+ y 2)}÷[w spring ×( x 0 2 + x 1 2 + x 2 2 )-( x 0+ x 1+ x 2) 2 ]; w spring The number of sampling points for spring water supply; a summer ={w summer ×[( x 3× y 3)+ ( x 4× y 4)+ ( x 5× y 5)]-( x 3+ x 4+ x 5)×( y 3+ y 4+ y 5)}÷[w summer ×(x 3 2 + x 4 2 + x 5 2 )-( x 3+ x 4+ x 5) 2 ]; w summer The number of sampling points for summer water supply; a autumn ={w autumn ×[( x 6× y 6)+ ( x 7× y 7)+ ( x 8× y 8)]-( x 6+ x 7+ x 8)×( y 6+ y 7+ y 8)}÷[w autumn ×( x 6 2 + x 7 2 + x 8 2 )-( x 6+ x 7+ x 8) 2 ]; w autumn number of sampling points for fall water supply; a winter ={w winter ×[( x 9× y 9)+ ( x 10 × y 10 )+ ( x 11 × y 11 )]-( x 9+ x 10 + x 11 )×( y 9+ y 10 + y 11 )}÷[w winter ×( x 9 2 + x10 2 + x 11 2 )-( x 9+ x 10 + x 11 ) 2 ]; w winter The number of sampling points for winter water supply; b spring =[( y 0+ y 1+ y 2)-a spring ×( x 0+ x 1+ x 2)]÷w spring ; b summer =[( y 3+ y 4+ y 5)-a summer ×( x 3+ x 4+ x 5)]÷w summer ; b autumn =[( y 6+ y 7+ y 8)-a autumn ×( x 6+ x 7+ x 8)]÷w autumn ; b winter =[( y 9+ y 10 + y 11 )-a winter ×( x 9+ x 10 + x 11 )]÷w winter .

[0033] C3, yes x By assigning a month, you can get the demand corresponding to the predicted month.

[0034] Step 5: Dynamically adjust water supply resource allocation, specifically: D1, open multiple water supply points in the water supply pipeline network until the sum of the monthly opened water supply points in the water supply pipeline network is greater than the monthly demand predicted by C3, otherwise continue to increase the number of opened water supply points; D2, there are multiple demand points in the water supply pipeline network. Each water supply point stops supplying water after reaching its own capacity. At the same time, the number of demand points is also changing dynamically. Suppose there are U water supply points and V demand points at a certain moment. Suppose from the water supply point i To the demand point j The water distribution cost is C ij , set up water supply points i Assign to demand points j The amount of water is D ij , demand point j The quantity demanded is Q j ; Set the prerequisites: , and adjust the allocation of water supply resources based on this premise; Set adjustment goals: ; Min() is the minimization function, which is used as the goal to dynamically adjust the allocation of water supply resources.

[0035] The water supply pipeline management system and method of the present invention have the following advantages: (1) Find the location of the leaking pipeline and issue a reminder based on the electronic map, which is a reliable idea; (2) The formulas for estimating pipeline life, operation and maintenance manpower, and monthly water supply demand are ingenious and creative; (3) Innovative and practical method for dynamically adjusting water supply resource allocation.

[0036] Other features and advantages of the present invention will be set forth in the following description or may be learned by practicing the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] The drawings are only for the purpose of illustrating particular embodiments and are not to be considered limiting of the present invention. Like reference symbols denote like components throughout the drawings.

[0038] Figure 1 It is a flow chart of the method of the present invention. DETAILED DESCRIPTION

[0039] A water supply pipeline management system and method of the present invention will be further described in detail below in conjunction with the accompanying drawings and embodiments.

[0040] The technical solution adopted by the present invention is a water supply pipeline management system: It includes leakage management module, pipeline life assessment module, operation and maintenance manpower estimation module, water supply demand forecasting module and water supply resource allocation module.

[0041] The leakage management module is used to find the section where the leaking pipe is located and to remind based on the electronic map, specifically in the following ways: A1, assume that the branches and intersections in the pipeline network are nodes of the pipeline, and the section between two adjacent nodes is a section of pipeline. Each section of pipeline is numbered in the form of the longitude and latitude coordinates of the upstream node of this section of pipeline + the longitude and latitude coordinates of the downstream node, and the unit of the coordinate value is degree; for example, the coordinates of the upstream node of a section of pipeline are (30.54321N, 100.12345E), and the coordinates of the downstream node are (30.54320N, 100.12347E), then the number of this section of pipeline is: 30.54321N100.12345E30.54320N100.12347E.

[0042] A2, let the length of the pipe between the nodes be L, the diameter be D, and the friction coefficient be f ; Then the theoretical pressure loss H flowing through each section of the pipeline f = f ×L÷D×( v 2 )÷(2×g), where v is the flow velocity of the fluid in the pipe, and g is the acceleration due to gravity.

[0043] A3, the actual pressure loss H flowing through each section of the pipeline e =P in -P out , where P in The water pressure at the upstream node of this section of pipeline, P out The water pressure at the node downstream of this pipe segment.

[0044] A4, when H e >H f When the pipeline is detected, it is determined that there is a water leakage problem in this section of the pipeline; on the electronic map, pipelines in normal state are indicated in green, and pipelines with water leakage problems are indicated in red. When there is a red pipeline, a prompt sound is issued.

[0045] The pipeline life assessment module is used to estimate the remaining life of the pipeline in the following manner: B1, create a database of historical data of pipeline life t, average temperature of pipeline external environment T, and water pressure P, where the unit of t is unified as year, the unit of T is unified as degree Celsius, and the unit of P is unified as Pascal.

[0046] B2, construct pipeline life distribution function: ;in: ; β0 is the basic life parameter, which is determined by material, thickness and diameter. θ1 and θ2 are the influence coefficients of the average temperature of the external environment and water pressure on the pipeline life respectively.

[0047] B3, use maximum likelihood to estimate the values ​​of β0, θ1, θ2: ①Define the survival function: ; Probability density function: ; ② For the same sample data of β0 (t1, t1, …, t n ), the likelihood function is: ; ③ Take the log-likelihood function: ; Expand the log-likelihood function: ; ④ Randomly initialize the values ​​of β0, θ1, and θ2, and calculate their partial derivatives ∂β0, ∂θ1, and ∂θ2 respectively. 0 new =β 0 old +α×∂β 0 old ,θ 1 new =θ 1 old +α×∂θ 1 old ,θ 2 new =θ 2 old +α×∂θ 2 old The law is iterated, where β 0 new ,θ 1 new ,θ 2 new represents the newly generated value; β 0 old ,θ 1 old ,θ 2 old represents the old value; α is the learning rate, which controls the step size of each update; until the value of the log-likelihood function reaches the maximum, the optimal values ​​of β0, θ1, and θ2 are obtained.

[0048] B4, evaluate the remaining life of the pipeline with the same β0: The β corresponding to the pipeline is calculated based on the average temperature and water pressure of the pipeline's external environment. The pipeline life distribution function is exponentially distributed, and the expected remaining life in years is approximately equal to the value of β. The value of β is the estimated remaining life of the pipeline.

[0049] The operation and maintenance manpower estimation module is used to estimate the manpower required this year to ensure the normal operation of the built pipelines. The specific method is as follows: Obtain the leakage status and remaining life of each pipe section, give priority to dealing with leaking pipes, and then replace pipes with nearly the same lifespan; sort pipes with a remaining lifespan less than (β+1) from least to most, and replace pipes with the lowest remaining lifespan first.

[0050] Suppose that there were m pipeline leaks in the previous year, and the average time spent on repairing leaks was K man-hours / time; suppose that it takes J0 man-hours / segment to replace pipelines on average, and the number of pipelines with a remaining life less than (β+1) at the beginning of the year is R; workers work h hours a day, then the number of pipeline operation and maintenance personnel to be prepared this year is N e =ceil [(m×K+R×J0)÷(365×h) ×(1+γ)], where γ is the proportion of reserve personnel and ceil() is the upward rounding function.

[0051] In actual implementation, γ is usually set to 10% to 20%.

[0052] The water supply demand forecasting module is used to estimate the monthly water supply demand data in the following way: C1, create a historical water supply database, including months and corresponding water supply volumes.

[0053] C2, create a prediction function according to the season, set the month x ,when x ≥12, converted to x =0; Suppose the spring prediction function is y spring = a spring × x + b spring , x =0, 1, 2; The summer prediction function is y summer = a summer × x + b summer , x =3, 4, 5; The autumn forecast function is y autumn = a autumn × x + b autumn , x =6,7,8; The winter prediction function is y winter = a winter × x + b winter , x =9, 10, 11; set up x 0. x 1. x 2. x 3. x 4. x 5. x 6. x7. x 8. x 9. x 10 , x 11 The months they represent are 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, and 11, and the corresponding water supply is y 0. y 1. y 2. y 3. y 4. y 5. y 6. y 7. y 8. y 9. y 10 , y 11 ; a spring ={w spring ×[( x 0× y 0)+ ( x 1× y 1)+ ( x 2× y 2)]-( x 0+ x 1+ x 2)×( y 0+ y 1+ y 2)}÷[w spring ×( x 0 2 + x 1 2 + x 2 2 )-( x 0+ x 1+ x 2) 2 ]; w spring The number of sampling points for spring water supply; a summer ={w summer ×[( x 3× y 3)+ ( x 4× y 4)+ ( x 5× y 5)]-( x 3+ x 4+ x 5)×( y 3+ y 4+ y5)}÷[w summer ×( x 3 2 + x 4 2 + x 5 2 )-( x 3+ x 4+ x 5) 2 ]; w summer The number of sampling points for summer water supply; a autumn ={w autumn ×[( x 6× y 6)+ ( x 7× y 7)+ ( x 8× y 8)]-( x 6+ x 7+ x 8)×( y 6+ y 7+ y 8)}÷[w autumn ×( x 6 2 + x 7 2 + x 8 2 )-( x 6+ x 7+ x 8) 2 ]; w autumn number of sampling points for fall water supply; a winter ={w winter ×[( x 9× y 9)+ ( x 10 × y 10 )+ ( x 11 × y 11 )]-( x 9+ x 10 + x 11 )×( y 9+ y 10 + y 11 )}÷[w winter ×( x 92 + x 10 2 + x 11 2 )-( x 9+ x 10 + x 11 ) 2 ]; w winter The number of sampling points for winter water supply; b spring =[( y 0+ y 1+ y 2)-a spring ×( x 0+ x 1+ x 2)]÷w spring ; b summer =[( y 3+ y 4+ y 5)-a summer ×( x 3+ x 4+ x 5)]÷w summer ; b autumn =[( y 6+ y 7+ y 8)-a autumn ×( x 6+ x 7+ x 8)]÷w autumn ; b winter =[( y 9+ y 10 + y 11 )-a winter ×( x 9+ x 10 + x 11 )]÷w winter .

[0054] C3, yes x By assigning a month, you can get the demand corresponding to the predicted month.

[0055] The water supply resource configuration module is used to dynamically adjust the allocation of water supply resources, specifically in the following ways: D1, open multiple water supply points in the water supply pipeline network until the sum of the monthly opened water supply points in the water supply pipeline network is greater than the monthly demand predicted by C3, otherwise continue to increase the number of opened water supply points; D2, there are multiple demand points in the water supply pipeline network. Each water supply point stops supplying water after reaching its own capacity. At the same time, the number of demand points is also changing dynamically. Suppose there are U water supply points and V demand points at a certain moment. Suppose from the water supply point i To the demand point j The water distribution cost is C ij , set up water supply points i Assign to demand points j The amount of water is D ij , demand point j The quantity demanded is Q j ; Set the prerequisites: , and adjust the allocation of water supply resources based on this premise; Set adjustment goals: ; Min() is the minimization function, which is used as the goal to dynamically adjust the allocation of water supply resources.

[0056] like Figure 1 As shown, a water supply pipeline management method: Step 1: Find the section where the leaking pipe is located and give a reminder based on the electronic map. The specific method is as follows: A1, assume that the branches and intersections in the pipeline network are nodes of the pipeline, and the section between two adjacent nodes is a section of pipeline. Each section of pipeline is numbered in the form of the longitude and latitude coordinates of the upstream node of this section of pipeline + the longitude and latitude coordinates of the downstream node, and the unit of the coordinate value is degree; for example, the coordinates of the upstream node of a section of pipeline are (30.54321N, 100.12345E), and the coordinates of the downstream node are (30.54320N, 100.12347E), then the number of this section of pipeline is: 30.54321N100.12345E30.54320N100.12347E.

[0057] A2, let the length of the pipe between the nodes be L, the diameter be D, and the friction coefficient be f ; Then the theoretical pressure loss H flowing through each section of the pipeline f = f ×L÷D×( v 2 )÷(2×g), where v is the flow velocity of the fluid in the pipe, and g is the acceleration due to gravity.

[0058] A3, the actual pressure loss H flowing through each section of the pipeline e =Pin -P out , where P in The water pressure at the upstream node of this section of pipeline, P out The water pressure at the node downstream of this pipe segment.

[0059] A4, when H e >H f When the pipeline is detected, it is determined that there is a water leakage problem in this section of the pipeline; on the electronic map, pipelines in normal state are indicated in green, and pipelines with water leakage problems are indicated in red. When there is a red pipeline, a prompt sound is issued.

[0060] Step 2: Estimate the remaining life of the pipeline by: B1, create a database of historical data of pipeline life t, average temperature of pipeline external environment T, and water pressure P, where the unit of t is unified as year, the unit of T is unified as degree Celsius, and the unit of P is unified as Pascal.

[0061] B2, construct pipeline life distribution function: ;in: ; β0 is the basic life parameter, which is determined by material, thickness and diameter. θ1 and θ2 are the influence coefficients of the average temperature of the external environment and water pressure on the pipeline life respectively.

[0062] B3, use maximum likelihood to estimate the values ​​of β0, θ1, θ2: ①Define the survival function: ; Probability density function: ; ② For the same sample data of β0 (t1, t1, …, t n ), the likelihood function is: ; ③ Take the log-likelihood function: ; Expand the log-likelihood function: ; ④ Randomly initialize the values ​​of β0, θ1, and θ2, and calculate their partial derivatives ∂β0, ∂θ1, and ∂θ2 respectively. 0 new =β 0 old +α×∂β 0 old ,θ 1 new =θ 1 old +α×∂θ 1 old ,θ 2 new =θ 2 old +α×∂θ 2 old The law is iterated, where β 0 new ,θ1 new ,θ 2 new represents the newly generated value; β 0 old ,θ 1 old ,θ 2 old represents the old value; α is the learning rate, which controls the step size of each update; until the value of the log-likelihood function reaches the maximum, the optimal values ​​of β0, θ1, and θ2 are obtained.

[0063] B4, evaluate the remaining life of the pipeline with the same β0: The β corresponding to the pipeline is calculated based on the average temperature and water pressure of the pipeline's external environment. The pipeline life distribution function is exponentially distributed, and the expected remaining life in years is approximately equal to the value of β. The value of β is the estimated remaining life of the pipeline.

[0064] Step 3: Estimate the manpower required this year to ensure the normal operation of the built pipelines. The specific method is as follows: Obtain the leakage status and remaining life of each pipe section, give priority to dealing with leaking pipes, and then replace pipes with nearly the same lifespan; sort pipes with a remaining lifespan less than (β+1) from least to most, and replace pipes with the lowest remaining lifespan first.

[0065] Suppose that there were m pipeline leaks in the previous year, and the average time spent on repairing leaks was K man-hours / time; suppose that it takes J0 man-hours / segment to replace pipelines on average, and the number of pipelines with a remaining life less than (β+1) at the beginning of the year is R; workers work h hours a day, then the number of pipeline operation and maintenance personnel to be prepared this year is N e =ceil [(m×K+R×J0)÷(365×h) ×(1+γ)], where γ is the proportion of reserve personnel and ceil() is the upward rounding function.

[0066] In actual implementation, γ is usually set to 10% to 20%.

[0067] Step 4: Estimate monthly water supply demand data by: C1, create a historical water supply database, including months and corresponding water supply volumes.

[0068] C2, create a prediction function according to the season, set the month x ,when x ≥12, converted to x =0; Suppose the spring prediction function is y spring = a spring × x + b spring , x =0, 1, 2; The summer prediction function is y summer = asummer × x + b summer , x =3, 4, 5; The autumn forecast function is y autumn = a autumn × x + b autumn , x =6,7,8; The winter prediction function is y winter = a winter × x + b winter , x =9, 10, 11; set up x 0. x 1. x 2. x 3. x 4. x 5. x 6. x 7. x 8. x 9. x 10 , x 11 The months they represent are 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, and 11, and the corresponding water supply is y 0. y 1. y 2. y 3. y 4. y 5. y 6. y 7. y 8. y 9. y 10 , y 11 ; a spring ={w spring ×[( x 0× y 0)+ ( x 1× y 1)+ ( x 2× y 2)]-( x 0+ x 1+ x 2)×( y 0+ y 1+ y 2)}÷[wspring ×( x 0 2 + x 1 2 + x 2 2 )-( x 0+ x 1+ x 2) 2 ]; w spring The number of sampling points for spring water supply; a summer ={w summer ×[( x 3× y 3)+ ( x 4× y 4)+ ( x 5× y 5)]-( x 3+ x 4+ x 5)×( y 3+ y 4+ y 5)}÷[w summer ×( x 3 2 + x 4 2 + x 5 2 )-( x 3+ x 4+ x 5) 2 ]; w summer The number of sampling points for summer water supply; a autumn ={w autumn ×[( x 6× y 6)+ ( x 7× y 7)+ ( x 8× y 8)]-( x 6+ x 7+ x 8)×( y 6+ y 7+ y 8)}÷[w autumn ×( x 6 2 + x 7 2 + x 8 2 )-( x 6+ x 7+ x8) 2 ]; w autumn number of sampling points for fall water supply; a winter ={w winter ×[( x 9× y 9)+ ( x 10 × y 10 )+ ( x 11 × y 11 )]-( x 9+ x 10 + x 11 )×( y 9+ y 10 + y 11 )}÷[w winter ×( x 9 2 + x 10 2 + x 11 2 )-( x 9+ x 10 + x 11 ) 2 ]; w winter The number of sampling points for winter water supply; b spring =[( y 0+ y 1+ y 2)-a spring ×( x 0+ x 1+ x 2)]÷w spring ; b summer =[( y 3+ y 4+ y 5)-a summer ×( x 3+ x 4+ x 5)]÷w summer ; b autumn =[( y 6+ y 7+ y8)-a autumn ×( x 6+ x 7+ x 8)]÷w autumn ; b winter =[( y 9+ y 10 + y 11 )-a winter ×( x 9+ x 10 + x 11 )]÷w winter .

[0069] C3, yes x By assigning a month, you can get the demand corresponding to the predicted month.

[0070] Step 5: Dynamically adjust water supply resource allocation, specifically: D1, open multiple water supply points in the water supply pipeline network until the sum of the monthly opened water supply points in the water supply pipeline network is greater than the monthly demand predicted by C3, otherwise continue to increase the number of opened water supply points; D2, there are multiple demand points in the water supply pipeline network. Each water supply point stops supplying water after reaching its own capacity. At the same time, the number of demand points is also changing dynamically. Suppose there are U water supply points and V demand points at a certain moment. Suppose from the water supply point i To the demand point j The water distribution cost is C ij , set up water supply points i Assign to demand points j The amount of water is D ij , demand point j The quantity demanded is Q j ; Set the prerequisites: , and adjust the allocation of water supply resources based on this premise; Set adjustment goals: ; Min() is the minimization function, which is used as the goal to dynamically adjust the allocation of water supply resources.

[0071] Furthermore, in the field of water supply pipelines, the word "pipeline" is usually used when it is necessary to describe a pipeline with detailed attribute information such as size and material or to observe a part of the pipeline closely; in software systems, pipelines are usually only displayed with thin lines, so the word "pipeline" is usually used for macroscopic description; a pipeline network refers to a comprehensive mesh system composed of multiple pipelines.

[0072] The above description is only a preferred specific implementation manner of the present invention, but the protection scope of the present invention is not limited thereto. Any changes or substitutions that can be easily conceived by any technician familiar with the technical field within the technical scope disclosed by the present invention should be covered within the protection scope of the present invention.

Claims

1. A water supply pipeline management system, characterized in that: It includes leakage management module, pipeline life assessment module, operation and maintenance manpower estimation module, water supply demand prediction module, and water supply resource allocation module; The leakage management module is used to find the section where the leaking pipe is located and to remind based on the electronic map, specifically in the following ways: A1, assume that the branches and intersections in the pipeline network are nodes of the pipeline, and the section between two adjacent nodes is a section of pipeline. Number each section of pipeline in the form of the longitude and latitude coordinates of the upstream node of this section of pipeline + the longitude and latitude coordinates of the downstream node, and the unit of the coordinate value is degree; A2, let the length of the pipe between the nodes be L, the diameter be D, and the friction coefficient be f ; Then the theoretical pressure loss H flowing through each section of the pipeline f = f ×L÷D×( v 2 )÷(2×g), where v is the flow velocity of the fluid in the pipeline, g is the acceleration due to gravity; A3, the actual pressure loss H flowing through each section of the pipeline e =P in -P out , where P in The water pressure at the upstream node of this section of pipeline, P out The water pressure at the downstream node of this section of pipeline; A4, when H e >H f When the pipeline is detected, it is determined that there is a water leakage problem in this section of the pipeline; on the electronic map, normal pipelines are indicated in green, and pipelines with water leakage problems are indicated in red. When there is a red pipeline, a warning sound is emitted; The pipeline life assessment module is used to estimate the remaining life of the pipeline in the following manner: B1, create a database of historical data of pipeline life t, average temperature of pipeline external environment T, and water pressure P, where the unit of t is unified as year, the unit of T is unified as Celsius, and the unit of P is unified as Pascal; B2, construct pipeline life distribution function: ;in: ; β0 is the basic life parameter, which is determined by material, thickness and diameter. θ1 and θ2 are the influence coefficients of the average temperature and water pressure of the external environment on the life of the pipeline respectively. B3, use maximum likelihood to estimate the values ​​of β0, θ1, θ2: ①Define the survival function: ; Probability density function: ; ② For the same sample data of β0 (t1, t1, …, t n ), the likelihood function is: ； ③ Take the log-likelihood function: ； Expand the log-likelihood function: ； ④ Randomly initialize the values ​​of β0, θ1, and θ2, and calculate their partial derivatives ∂β0, ∂θ1, and ∂θ2 respectively. 0 new =β 0 old +α×∂β 0 old ,θ 1 new =θ 1 old +α×∂θ 1 old ,θ 2 new =θ 2 old +α×∂θ 2 old The law is iterated, where β 0 new ,θ 1 new ,θ 2 new represents the newly generated value; β 0 old ,θ 1 old ,θ 2 old represents the old value; α is the learning rate, which controls the step size of each update; until the log-likelihood function value reaches the maximum, the optimal values ​​of β0, θ1, and θ2 are obtained; B4, evaluate the remaining life of the pipeline with the same β0: According to the average temperature and water pressure of the external environment of the pipeline, calculate the corresponding β of the pipeline; the pipeline life distribution function is exponentially distributed, and the expected remaining life years are approximately equal to the value of β, which is the estimated remaining life of the pipeline; The operation and maintenance manpower estimation module is used to estimate the manpower required this year to ensure the normal operation of the built pipelines. The specific method is as follows: Obtain the leakage status and remaining life of each pipe section, and deal with leaking pipes first, followed by replacing pipes with a near-term life span; sort pipes with a remaining life span less than (β+1) from least to most, and replace pipes with a lower remaining life span first; Suppose that there were m pipeline leaks in the previous year, and the average time spent on repairing leaks was K man-hours / time; suppose that it takes J0 man-hours / segment to replace pipelines on average, and the number of pipelines with a remaining life less than (β+1) at the beginning of the year is R; workers work h hours a day, then the number of pipeline operation and maintenance personnel to be prepared this year is N e =ceil [(m×K+R×J0)÷(365×h) ×(1+γ)], where γ is the proportion of reserve personnel and ceil() is the upward rounding function; The water supply demand forecasting module is used to estimate the monthly water supply demand data in the following way: C1, create a historical water supply database, including months and corresponding water supply volumes; C2, create a prediction function according to the season, set the month x ,when x ≥12, converted to x =0; Suppose the spring prediction function is y spring = a spring × x + b spring , x =0, 1, 2; The summer prediction function is y summer = a summer × x + b summer , x =3, 4, 5; The autumn forecast function is y autumn = a autumn × x + b autumn , x =6,7,8; The winter prediction function is y winter = a winter × x + b winter , x =9, 10, 11; set up x 0. x 1. x 2. x 3. x 4. x 5. x 6. x 7. x 8. x 9. x 10 , x 11 The months they represent are 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, and 11, and the corresponding water supply is y 0. y 1. y 2. y 3. y 4. y 5. y 6. y 7. y 8. y 9. y 10 , y 11 ; a spring ={w spring ×[( x 0× y 0)+ ( x 1× y 1)+ ( x 2× y 2)]-( x 0+ x 1+ x 2)×( y 0+ y 1+ y 2)}÷[w spring ×( x 0 2 + x 1 2 + x 2 2 )-( x 0+ x 1+ x 2) 2 ]; w spring The number of sampling points for spring water supply; a summer ={w summer ×[( x 3× y 3)+ ( x 4× y 4)+ ( x 5× y 5)]-( x 3+ x 4+ x 5)×( y 3+ y 4+ y 5)}÷[w summer ×( x 3 2 + x 4 2 + x 5 2 )-( x 3+ x 4+ x 5) 2 ]; w summer The number of sampling points for summer water supply; a autumn ={w autumn ×[( x 6× y 6)+ ( x 7× y 7)+ ( x 8× y 8)]-( x 6+ x 7+ x 8)×( y 6+ y 7+ y 8)}÷[w autumn ×( x 6 2 + x 7 2 + x 8 2 )-( x 6+ x 7+ x 8) 2 ]; w autumn number of sampling points for fall water supply; a winter ={w winter ×[( x 9× y 9)+ ( x 10 × y 10 )+ ( x 11 × y 11 )]-( x 9+ x 10 + x 11 )×( y 9+ y 10 + y 11 )}÷[w winter ×( x 9 2 + x 10 2 + x 11 2 )-( x 9+ x 10 + x 11 ) 2 ]; w winter The number of sampling points for winter water supply; b spring =[( y 0+ y 1+ y 2)-a spring ×( x 0+ x 1+ x 2)]÷w spring ; b summer =[( y 3+ y 4+ y 5)-a summer ×( x 3+ x 4+ x 5)]÷w summer ; b autumn =[( y 6+ y 7+ y 8)-a autumn ×( x 6+ x 7+ x 8)]÷w autumn ; b winter =[( y 9+ y 10 + y 11 )-a winter ×( x 9+ x 10 + x 11 )]÷w winter ; C3, yes x Assign a month to get the demand corresponding to the predicted month; The water supply resource configuration module is used to dynamically adjust the allocation of water supply resources, specifically in the following ways: D1, open multiple water supply points in the water supply pipeline network until the sum of the monthly opened water supply points in the water supply pipeline network is greater than the monthly demand predicted by C3, otherwise continue to increase the number of opened water supply points; D2, there are multiple demand points in the water supply pipeline network. Each water supply point stops supplying water after reaching its own capacity. At the same time, the number of demand points is also changing dynamically. Suppose there are U water supply points and V demand points at a certain moment. Suppose from the water supply point i To the demand point j The water distribution cost is C ij , set up water supply points i Assign to demand points j The amount of water is D ij , demand point j The quantity demanded is Q j ; Set the prerequisites: , and adjust the allocation of water supply resources based on this premise; Set adjustment goals: ; Min() is the minimization function, which is used as the goal to dynamically adjust the allocation of water supply resources.

2. A water supply pipeline management method, characterized in that: Step 1: Find the section where the leaking pipe is located and give a reminder based on the electronic map. The specific method is as follows: A1, assume that the branches and intersections in the pipeline network are nodes of the pipeline, and the section between two adjacent nodes is a section of pipeline. Number each section of pipeline in the form of the longitude and latitude coordinates of the upstream node + the longitude and latitude coordinates of the downstream node of this section of pipeline, and the unit of the coordinate value is degree; A2, let the length of the pipe between the nodes be L, the diameter be D, and the friction coefficient be f ; Then the theoretical pressure loss H flowing through each section of the pipeline f = f ×L÷D×( v 2 )÷(2×g), where v is the flow velocity of the fluid in the pipeline, g is the acceleration due to gravity; A3, the actual pressure loss H flowing through each section of the pipeline e =P in -P out , where P in The water pressure at the upstream node of this section of pipeline, P out The water pressure at the downstream node of this section of pipeline; A4, when H e >H f When the pipeline is detected, it is determined that there is a water leakage problem in this section of the pipeline; on the electronic map, normal pipelines are indicated in green, and pipelines with water leakage problems are indicated in red. When there is a red pipeline, a warning sound is emitted; Step 2: Estimate the remaining life of the pipeline by: B1, create a database of historical data of pipeline life t, average temperature of pipeline external environment T, and water pressure P, where the unit of t is unified as year, the unit of T is unified as Celsius, and the unit of P is unified as Pascal; B2, construct pipeline life distribution function: ;in: ; β0 is the basic life parameter, which is determined by material, thickness and diameter. θ1 and θ2 are the influence coefficients of the average temperature and water pressure of the external environment on the life of the pipeline respectively. B3, use maximum likelihood to estimate the values ​​of β0, θ1, θ2: ①Define the survival function: ; Probability density function: ; ② For the same sample data of β0 (t1, t1, …, t n ), the likelihood function is: ； ③ Take the log-likelihood function: ； Expand the log-likelihood function: ； ④ Randomly initialize the values ​​of β0, θ1, and θ2, and calculate their partial derivatives ∂β0, ∂θ1, and ∂θ2 respectively. 0 new =β 0 old +α×∂β 0 old ,θ 1 new =θ 1 old +α×∂θ 1 old ,θ 2 new =θ 2 old +α×∂θ 2 old The law is iterated, where β 0 new ,θ 1 new ,θ 2 new represents the newly generated value; β 0 old ,θ 1 old ,θ 2 old represents the old value; α is the learning rate, which controls the step size of each update; until the log-likelihood function value reaches the maximum, the optimal values ​​of β0, θ1, and θ2 are obtained; B4, evaluate the remaining life of the pipeline with the same β0: According to the average temperature and water pressure of the external environment of the pipeline, calculate the corresponding β of the pipeline; the pipeline life distribution function is exponentially distributed, and the expected remaining life years are approximately equal to the value of β, which is the estimated remaining life of the pipeline; Step 3: Estimate the manpower required this year to ensure the normal operation of the built pipelines. The specific method is as follows: Obtain the leakage status and remaining life of each pipe section, and deal with leaking pipes first, followed by replacing pipes with a near-term life span; sort pipes with a remaining life span less than (β+1) from least to most, and replace pipes with a lower remaining life span first; Suppose that there were m pipeline leaks in the previous year, and the average time spent on repairing leaks was K man-hours / time; suppose that it takes J0 man-hours / segment to replace pipelines on average, and the number of pipelines with a remaining life less than (β+1) at the beginning of the year is R; workers work h hours a day, then the number of pipeline operation and maintenance personnel to be prepared this year is N e =ceil [(m×K+R×J0)÷(365×h) ×(1+γ)], where γ is the proportion of reserve personnel and ceil() is the upward rounding function; Step 4: Estimate monthly water supply demand data by: C1, create a historical water supply database, including months and corresponding water supply volumes; C2, create a prediction function according to the season, set the month x ,when x ≥12, converted to x =0; Suppose the spring prediction function is y spring = a spring × x + b spring , x =0, 1, 2; The summer prediction function is y summer = a summer × x + b summer , x =3, 4, 5; The autumn forecast function is y autumn = a autumn × x + b autumn , x =6,7,8; The winter prediction function is y winter = a winter × x + b winter , x =9, 10, 11; set up x 0. x 1. x 2. x 3. x 4. x 5. x 6. x 7. x 8. x 9. x 10 , x 11 The months they represent are 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, and 11, and the corresponding water supply is y 0. y 1. y 2. y 3. y 4. y 5. y 6. y 7. y 8. y 9. y 10 , y 11 ; a spring ={w spring ×[( x 0× y 0)+ ( x 1× y 1)+ ( x 2× y 2)]-( x 0+ x 1+ x 2)×( y 0+ y 1+ y 2)}÷[w spring ×( x 0 2 + x 1 2 + x 2 2 )-( x 0+ x 1+ x 2) 2 ]; w spring The number of sampling points for spring water supply; a summer ={w summer ×[( x 3× y 3)+ ( x 4× y 4)+ ( x 5× y 5)]-( x 3+ x 4+ x 5)×( y 3+ y 4+ y 5)}÷[w summer ×( x 3 2 + x 4 2 + x 5 2 )-( x 3+ x 4+ x 5) 2 ]; w summer The number of sampling points for summer water supply; a autumn ={w autumn ×[( x 6× y 6)+ ( x 7× y 7)+ ( x 8× y 8)]-( x 6+ x 7+ x 8)×( y 6+ y 7+ y 8)}÷[w autumn ×( x 6 2 + x 7 2 + x 8 2 )-( x 6+ x 7+ x 8) 2 ]; w autumn number of sampling points for fall water supply; a winter ={w winter ×[( x 9× y 9)+ ( x 10 × y 10 )+ ( x 11 × y 11 )]-( x 9+ x 10 + x 11 )×( y 9+ y 10 + y 11 )}÷[w winter ×( x 9 2 + x 10 2 + x 11 2 )-( x 9+ x 10 + x 11 ) 2 ]; w winter The number of sampling points for winter water supply; b spring =[( y 0+ y 1+ y 2)-a spring ×( x 0+ x 1+ x 2)]÷w spring ; b summer =[( y 3+ y 4+ y 5)-a summer ×( x 3+ x 4+ x 5)]÷w summer ; b autumn =[( y 6+ y 7+ y 8)-a autumn ×( x 6+ x 7+ x 8)]÷w autumn ; b winter =[( y 9+ y 10 + y 11 )-a winter ×( x 9+ x 10 + x 11 )]÷w winter ; C3, yes x Assign a month to get the demand corresponding to the predicted month; Step 5: Dynamically adjust water supply resource allocation, specifically: D1, open multiple water supply points in the water supply pipeline network until the sum of the monthly opened water supply points in the water supply pipeline network is greater than the monthly demand predicted by C3, otherwise continue to increase the number of opened water supply points; D2, there are multiple demand points in the water supply pipeline network. Each water supply point stops supplying water after reaching its own capacity. At the same time, the number of demand points is also changing dynamically. Suppose there are U water supply points and V demand points at a certain moment. Suppose from the water supply point i To the demand point j The water distribution cost is C ij , set up water supply points i Assign to demand points j The amount of water is D ij , demand point j The quantity demanded is Q j ; Set the prerequisites: , and adjust the allocation of water supply resources based on this premise; Set adjustment goals: ; Min() is the minimization function, which is used as the goal to dynamically adjust the allocation of water supply resources.

Images