Water supply network node water demand initialization method based on dynamic population distribution

Abstract

The invention belongs to the technical field of hydraulic model correction, and aims to solve the problems that the current bottom-up method is tedious, time-consuming and difficult to reflect real-time water consumption change and the node water demand obtained by the top-down method is not consistent with the actual condition. The invention provides a dynamic population distribution-based water supply network node water demand initialization method, which comprises the following steps of: calculating the population number in each online metering water meter service area and the population number in each grid unit in a water supply range at the current moment based on a mobile phone signaling data platform; and calculating the population number of off-line metering water in each water supply grid unit at the current moment, and then initializing the node water demand of the water supply network at the current moment based on the on-line metering water data at the current moment and the population number of off-line metering water in each node service area. According to the method, the population dynamics and the influence of the population dynamics on the water demand can be accurately reflected, so that the node water demand initialization in the high population flow area is more in line with the actual situation.

Description

Technical Field

[0001] The invention belongs to the technical field of hydraulic model correction, and in particular relates to a method for initializing water demand of a water supply network node based on dynamic population distribution. Background Art

[0002] The water demand at the nodes of the water supply network is the driving force for the flow of water in the pipeline. When the hydraulic model of the water supply network is used to simulate the network operation parameters (such as node pressure and pipeline flow), the accuracy of the simulation depends on the accurate allocation of the model node water demand. If the node water demand is not initialized properly, it will lead to a significant deviation between the model simulation value and the actual observation value, thus affecting the optimal design, operation and maintenance decision of the water supply network based on the hydraulic model. Therefore, the reasonable initialization of the node water demand is the key to achieve high-precision simulation of the hydraulic model.

[0003] Existing node water demand initialization methods are mainly divided into two categories: bottom-up initialization methods and top-down initialization methods, such as Figure 1 As shown in the figure. The bottom-up node water demand initialization method refers to starting from the user's meter reading data, assigning a geographic coordinate to each user's water meter and linking it to the nearest hydraulic model node: (1) Revenue water allocation: For online metering users, the real-time water consumption is directly allocated to the corresponding node. For other users (customers with monthly meter reading or IC card), the water consumption data of the past month is counted, and the daily average water consumption is calculated. Then, it is multiplied by the water consumption coefficient at each moment to obtain the water consumption at each moment, and finally allocated to the corresponding node; (2) Non-revenue water (total water supply - revenue water) allocation: It is allocated to the corresponding node in proportion to the service area of ​​each node or the length of the pipeline adjacent to each node.

[0004] The top-down node demand initialization method is the exact opposite of the bottom-up initialization method. Instead of starting with the billing records of individual customers, this method starts from the water source (the "upper layer") and then systematically allocates it downward to each node. Online metered water consumption is first allocated to the nearest node, and then the remaining non-online metered water consumption is allocated in proportion to the area of ​​each node's service area or the historical statistical population in the service area. Thiessen polygons are a representative method for determining node service areas. It divides the entire water supply area into multiple polygons, each centered on a node, ensuring that any user within the polygon is closer to its node than any other node.

[0005] However, the existing node water demand initialization method has many defects, as follows:

[0006] (1) Bottom-up node water demand initialization method: ① This method requires detailed linking of the customer’s billing records with the model nodes at the beginning. This process is very cumbersome and time-consuming, and subsequent data maintenance and updates still require continuous investment, especially when user information or service areas change. Maintenance is difficult; ② This method relies on historical data or average water demand calculations, which may not accurately reflect real-time water use changes, resulting in inconsistency with actual conditions; ③ For small users, water demand calculations may be affected by greater random fluctuations, thereby increasing uncertainty.

[0007] (2) Top-down node water demand initialization method: ① This method allocates water demand in proportion based on the node service area and historical population data, ignoring the actual differences between different areas (such as residential and commercial). In addition, outdated demographic data (derived from the census) fail to reflect the real-time changes in population mobility and water demand in a timely manner, resulting in a significant difference between the initialized node water demand and the actual value, affecting the simulation accuracy of the hydraulic model. ② In areas with high population mobility, this method cannot accurately reflect the water demand at different times and fails to adjust the water demand allocation in a timely manner, resulting in the initialization of the node water demand inconsistent with the actual water demand. ③ The actual structure of the water supply network and the flow direction of piped water are not considered, which easily leads to a mismatch between the allocated node water demand and the water supply capacity of the network. Summary of the invention

[0008] The purpose of the present invention is to overcome the deficiencies of the prior art and to provide a method for initializing water demand at a water supply network node based on dynamic population distribution.

[0009] To achieve the above object, the technical solution adopted by the present invention is:

[0010] A method for initializing water demand of water supply network nodes based on dynamic population distribution comprises the following steps:

[0011] S1. Calculate the population in the service area of ​​each online water meter and the population in each grid unit at the current moment;

[0012] S2. Based on step S1, calculate the non-online metered water service area in each grid unit, and calculate the population of the non-online metered water service area in each grid unit at the current moment;

[0013] S3. Based on step S2 and the online water metering data at the current moment, initialize the water demand of the water supply network node at the current moment.

[0014] Preferably, step S1 comprises the following steps:

[0015] S11. According to the service area of ​​each online metering water meter, the service area of ​​all online metering water meters in the water supply network is obtained, that is, the entire water supply area. The expression of the entire water supply area is:

[0016] Γ={γ1,γ2,...,γ m}

[0017] In the formula, γ i is the service area of ​​the online metering water meter, i=1, 2, ..., m, Γ is the set of service areas of all online metering water meters in the water supply network;

[0018] S12, dividing the entire water supply area obtained in step S11 into grid units to obtain multiple grid units, and the set expression of all grid units is:

[0019] G = {g1, g2, ..., g n}

[0020] In the formula, g j is the grid unit divided into the entire water supply area, j = 1, 2, ..., n, the size of the grid unit is determined by the accuracy of the mobile phone signaling data platform, and G is the set of all grid units;

[0021] S13. The entire water supply area and multiple grid units are input as input parameters to the mobile signaling data platform to calculate the population in the service area of ​​each online water meter and the population in each grid unit at the current moment.

[0022] Preferably, step S2 comprises the following steps:

[0023] S21. Intersect any online water meter service area with the set of all grid cells to obtain at least one sub-area. The expression of the sub-area is:

[0024]

[0025] In the formula, γ i ∩g j Represents γ i With g j Intersection operation; γ i is the service area of ​​the online water meter, g j is the grid unit divided into the entire water supply area, Λ is γ i With g j The set of sub-regions generated by the intersection operation; K is γ i The number of sub-regions generated after the set G of all grid cells is intersected and divided; λ ik For γ i Sub-regions after segmentation;

[0026] S22, according to the population number in the service area of ​​each online water meter obtained in step S1 at the current moment, combined with the sub-area obtained in step S21, the population number in the sub-area at the current moment is calculated, and the calculation formula is as follows:

[0027]

[0028] In the formula, is the sub-region λ ik The current population of the area; The current online water meter service area γ i The population within is the sub-region λ ik area; For γ i the area of ​​the region;

[0029] S23. According to step S21, all sub-areas in any grid unit are removed to obtain the non-online metered water service area in the grid unit, which is expressed as:

[0030] g′ j =g j \Λ,g′ j ∈G′

[0031] In the formula, g′ j For non-online water metering service area, g j \Λ is g j The difference between Λ and G′ is g′ j (j=1, 2, ..., n);

[0032] S24, according to the population number in each grid unit obtained in step S1, the population number in the sub-area at the current moment obtained in step S22 and the non-online metered water service area obtained in step S23, the population number in the non-online metered water service area in each grid unit at the current moment is obtained, and the calculation formula is as follows:

[0033]

[0034] In the formula, is the grid unit g j The population within g j The current population of all sub-areas within; is g′ j The population in the region, g′ j This is a non-online water metering service area.

[0035] Preferably, step S3 comprises the following steps:

[0036] S31, taking each water supply network node as the center, draw a Thiessen polygon Z = {z1, z2, ..., z N};

[0037] S32, performing an intersection operation on the non-online water metering service area in all grid cells obtained in step S2 and the Thiessen polygon obtained in step S31 to obtain a set of segmented sub-areas, expressed as:

[0038]

[0039] In the formula, s lr is g′ j A single sub-region after being divided by Z, g′ j is the non-online water metering service area, G′ is g′ j A collection of components;

[0040] S33. Based on step S32, the number of people using non-online metered water in the service area of ​​each water supply network node at the current moment is calculated, and the calculation formula is:

[0041]

[0042] In the formula, s lr is g′ j A single sub-region after being divided by Z; For sub-region s lr area; is the region g′ j area;

[0043] S34, based on step S33 and the online water metering data at the current moment, initialize the water demand of the water supply network node at the current moment, and obtain the initial value of the water demand of the water supply network node. The calculation formula is:

[0044]

[0045] Where, d j is the initial value of water demand at the water supply network node j at the current moment; d non is the non-metered water demand at node j at the current moment; d tol is the total non-metered water demand of the water supply network at the current moment; P tol is the total population of the entire water supply area at the current moment; ∑d online is the total water demand of all online water meters corresponding to node j at the current moment, P j It is the number of people using non-online metered water in the service area of ​​each water supply network node at the current moment.

[0046] Compared with the prior art, the present invention has the following beneficial effects:

[0047] (1) The present invention calculates the population in the service area of ​​each online water meter and the population in each grid unit within the water supply range at the current moment through a real-time mobile signaling data platform. Since the mobile signaling data platform has a wide coverage and high user penetration rate, the present invention improves the accuracy of population mobility analysis, can accurately reflect population dynamics and its impact on water demand, and makes the water demand initialization of nodes in high population mobility areas more in line with actual conditions;

[0048] (2) The method of the present invention performs well in urban environments with large population mobility, has strong adaptability, is suitable for complex urban water supply systems, and has strong robustness; specifically: the water demand of the pipe network node is positively correlated with the number of people within its service area, that is, the larger the population, the greater the water demand. In an urban environment with significant temporal and spatial population mobility, the number of people in the service area of ​​each node will change significantly over time. There is often a large deviation between the traditional water demand initialization method and the actual water demand. The method proposed in the present invention can accurately determine the dynamic population of each node service area at each moment. Compared with traditional methods, it is more suitable for complex urban water supply systems with strong population mobility;

[0049] (3) The method provided by the present invention is seamlessly integrated with the existing geographic information system and water supply network data platform, thereby enhancing the wide application range. BRIEF DESCRIPTION OF THE DRAWINGS

[0050] Figure 1 Initialize flow charts for existing top-down and bottom-up node water demands;

[0051] Figure 2 A flow chart of a method for initializing water demand of a water supply network node based on dynamic population distribution provided by an embodiment of the present invention;

[0052] Figure 3 An example diagram for calculating the population in the service area of ​​each online water meter and the population in each grid cell within the water supply area at the current moment;

[0053] in, Figure 3 (a) is the topological structure and water supply range of an example water supply network; Figure 3 (b) water supply service area and grid unit; Figure 3 (c) is a mobile phone signaling data platform; Figure 3 (d) is the population of the current sample area;

[0054] Figure 4 This is an example diagram for calculating the number of people using non-online water metering in each grid cell at the current moment;

[0055] in, Figure 4 (a) is the water supply grid unit and the population of each grid unit at the current moment; Figure 4 (b) is the number of sub-regions and the population using online metered water in each sub-region at the current moment; Figure 4 (c) is the number of people using non-online water metering at the current moment;

[0056] Figure 5 An example diagram for calculating the water demand of the water supply network node at the initialization moment;

[0057] in, Figure 5 (a) is a schematic diagram of the Thiessen polygon and the non-online metered water service area in all grid cells; Figure 5 (b) is the area map of the sub-region obtained by segmentation; Figure 5 (c) is an example diagram for calculating the number of people using non-online water metering within the service area of ​​node j1 at the current moment. DETAILED DESCRIPTION

[0058] The following is a description of the embodiments of the present invention. Figures 1 to 5 , the technical solutions in the embodiments of the present invention are clearly and completely described. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of them. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0059] The embodiment of the present invention provides a method for initializing water demand of a water supply network node based on dynamic population distribution, which specifically includes the following steps:

[0060] S1, calculating the population in the service area of ​​each online water meter and the population in each grid unit at the current moment, including the following steps:

[0061] S11. According to the service area of ​​each online metering water meter, the service area of ​​all online metering water meters in the water supply network is obtained, that is, the entire water supply area. The expression of the entire water supply area is:

[0062] Γ={γ1,γ2,...,γ m} Formula 1-1

[0063] In the formula, γ i is the service area of ​​the online metering water meter, i = 1, 2, ..., m, m is the total number of online metering water meters, γ i Usually the plane outline of the building where the water meter is located is taken as its service area, or it is determined by the user according to the metering area. Γ is the set of service areas of all online metering water meters in the water supply network;

[0064] S12, dividing the entire water supply area obtained in step S11 into grid units to obtain multiple grid units, and the set expression of all grid units is:

[0065] G = {g1, g2, ... g n} Formula 1-2

[0066] In the formula, g j is the grid unit divided into the entire water supply area, j = 1, 2, ..., n, n is the total number of divided grid units, and the size of the grid unit is determined by the accuracy of the mobile phone signaling data platform provided by the operator. For example, if the platform provides a real-time population of 250m×250m, the grid unit size is set to 250m×250m; G is the set of all grid units;

[0067] S13, the entire water supply area Γ and multiple grid units G are input as input parameters to the mobile signaling data platform provided by the operator, and the population in the service area of ​​each online water meter and the population in each grid unit at the current moment are calculated. The mobile signaling data platform used in the embodiment of the present invention is a prior art, specifically a mobile signaling data platform provided by China Mobile, China Unicom or a telecommunications company, which provides paid services to users, and users can pay to obtain the regional population at each moment.

[0068] S2. Based on step S1, the non-online metered water service area in each grid unit is calculated, and the population of the non-online metered water service area in each grid unit at the current moment is calculated, including the following steps:

[0069] S21, any online water meter service area γ i (i=1, 2, ..., m) performs an intersection operation with the set G of all grid cells to obtain at least one sub-region. Specifically, when γ i completely within a grid cell g j , j = 1, 2, ..., n, the segmentation results in one sub-region, and when γ i When it intersects with the boundaries of some grid cells, it is divided into multiple sub-regions, and the expression of the sub-region is:

[0070]

[0071] In the formula, γ i ∩g j Represents γ i With g j Intersection operation; γ i is the service area of ​​the online water meter, g j is the grid unit divided into the entire water supply area, Λ is γ i With gj The set of sub-regions generated by the intersection operation; K is γ i The number of sub-regions generated after the set G of all grid cells is intersected and divided; λ ik (k=1,2,...,K) is γ i The sub-region after segmentation;

[0072] S22, according to the population number in the service area of ​​each online water meter obtained in step S1 at the current moment, combined with the sub-area obtained in step S21, the population number in the sub-area at the current moment is calculated, and the calculation formula is as follows:

[0073]

[0074] In the formula, is the sub-region λ ik The current population of the area; The current online water meter service area γ i The population within is the sub-region λ ik area; For γ i the area of ​​the region;

[0075] S23, according to step S21, any grid unit g j All sub-areas in (j=1, 2, ..., n) are removed to obtain the non-online water metering service area g′ in the grid unit. j , g′ j All users are non-online water users, g′ j The expression is:

[0076] g′ j =g j \Λ,g′ j ∈G′ Formula 1-6

[0077] In the formula, g′ j For non-online water metering service area, g j \Λ is g j The difference between Λ and G′ is g′ j (j=1, 2, ..., n);

[0078] S24, according to the population number in each grid unit obtained in step S1, the population number in the sub-area at the current moment obtained in step S22 and the non-online metered water service area obtained in step S23, the population number in the non-online metered water service area in each grid unit at the current moment is obtained, and the calculation formula is as follows:

[0079]

[0080] In the formula, is the grid unit g j The population within g j The current population of all sub-areas within; is g′ j The population in the region, g′ j It is a non-online metered water service area (all users are non-online metered water users).

[0081] S3, based on step S2 and the online water metering data at the current moment, initializing the water demand of the water supply network node at the current moment, including the following steps:

[0082] S31, taking each water supply network node j (j=1, 2, ..., N, N is the number of water supply network nodes) as the center, draw a Thiessen polygon Z={z1, z2, ..., z N}, and adjust the boundary of each Thiessen polygon according to the actual situation of the water supply network (water supply area, administrative boundary, etc.), and the Thiessen polygon corresponding to each node after adjustment is the water supply service area of ​​the node;

[0083] S32, performing an intersection operation on the non-online water metering service area G′ in all grid cells obtained in step S2 and the Thiessen polygon Z obtained in step S31 to obtain a segmented sub-area set S, where S is expressed as:

[0084]

[0085] In the formula, s lr is g′ j A single sub-region after being divided by Z, g′ j is the non-online water metering service area, G′ is g′ j A collection of components;

[0086] S33. Based on step S32, the number of people using non-online metered water in the service area of ​​each water supply network node at the current moment is calculated, and the calculation formula is:

[0087]

[0088] In the formula, s lr is g′ j A single sub-region after being divided by Z; For sub-region s lr area; is the region g′ j area;

[0089] S34, based on step S33 and the online water metering data at the current moment, initialize the water demand of the water supply network node at the current moment, and obtain the initial value of the water demand of the water supply network node. The calculation formula is:

[0090]

[0091] Where, d j is the initial value of water demand at the water supply network node j (j=1, 2, ..., N) at the current moment; d non is the non-metered water demand at node j at the current moment; d tol P is the total non-metered water demand of the water supply network at the current moment (total water supply - total online metered water); tol is the total population of the entire water supply area at the current moment; ∑d online is the total water demand of all online water meters corresponding to node j at the current moment, P j It is the number of people using non-online metered water in the service area of ​​each water supply network node at the current moment.

[0092] The following takes a simple water supply network as an example to further illustrate the specific implementation process of the method for initializing water demand of water supply network nodes based on dynamic population distribution provided by the embodiment of the present invention.

[0093] like Figure 3 As shown in (a), the water supply network consists of 10 nodes and 9 pipes, and there is 1 online water meter in the water supply area.

[0094] 1. Calculate the population in the service area of ​​each online water meter and the population in each grid unit within the water supply range at the current moment:

[0095] (1) Take the outline of the building where the online water meter is located as its water supply service area ( Figure 3 (b)), then Γ={γ1}; Since the grid accuracy of the mobile signaling data platform provided by the mobile operator is 250m×250m, the entire water supply area is divided into 8 grid units according to the grid accuracy ( Figure 3 b), that is, G = {g1, g2, ..., g8}.

[0096] (2) Figure 3 (c) and Figure 3 As shown in (d), Γ and G are used as input parameters and input into the mobile phone signaling data platform provided by the mobile operator to obtain the population of each region at the current moment, that is,

[0097] 2. Calculate the number of people using non-online water metering in each grid cell at the current moment:

[0098] (1) The intersection operation between the online water meter service area γ1 and the grid unit set G ( Figure 4 (a)) is divided into two sub-regions λ 11 and λ 12 ( Figure 4 (b)), that is, Λ={λ 11 ,λ 12}, Their areas are and According to the proportion of the sub-region to the area of ​​γ1, the number of people using online water metering in each sub-region at the current moment is and

[0099] (2) For the online water meter service area γ1, only g2 and g3 intersect with γ1, and the outlines and populations of the remaining grid cells remain unchanged. The grid cells after removing γ1 from g2 and g3 and updating are as follows: Figure 4 As shown in (c), the number of people using non-online water metering at the current moment is and

[0100] 3. Initialize the water demand of the water supply network node at the current moment:

[0101] (1) With each water supply network node j (j = 1, 2, ..., 10) as the center, draw a Thiessen polygon Z = {z1, z2, ..., z 10},like Figure 5 (a) shown.

[0102] (2) G′ and Z are intersected, and Z is divided into multiple sub-regions. The area of ​​each sub-region is as follows: Figure 5 (b) As shown in Figure 2, taking node j1 as an example, its service area z1 is divided into s 11 and 51 There are two sub-regions, whose areas are and If Figure 5 As shown in (C), the number of people using non-online water metering in the service area of ​​node j1 at the current moment is At the current moment, the number of people using non-online water metering in the service areas of all other nodes is The total population of non-online water metering users in the water supply network at the current moment is

[0103] (3) Based on the above-calculated population of non-online water metering users and the current online water metering data, the node water demand at the current moment is initialized. The water volume of the online water meter at the current moment is 3.6 L / s, and the total water supply of the example pipe network at the current moment is 67.6 L / s. Therefore, the non-online water demand at the current moment is d tol =67.6-3.6=64L / s. Taking node j4 as an example, the initial value of the water demand of node j4 at the current moment is The initial values ​​of water demand of the remaining nodes at the current moment are:

[0104] Although the embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that various changes, modifications, substitutions and variations may be made to the embodiments without departing from the principles and spirit of the present invention, and that the scope of the present invention is defined by the claims and their equivalents.

Claims

1. A method for initializing water demand of water supply network nodes based on dynamic population distribution, characterized in that: The following steps are involved: S1. Calculate the population in the service area of ​​each online water meter and the population in each grid unit at the current moment; S2. Based on step S1, calculate the non-online metered water service area in each grid unit, and calculate the population of the non-online metered water service area in each grid unit at the current moment; S3. Based on step S2 and the online water metering data at the current moment, initialize the water demand of the water supply network node at the current moment.

2. The method for initializing water demand of water supply network nodes based on dynamic population distribution according to claim 1, characterized in that: Step S1 includes the following steps: S11. According to the service area of ​​each online metering water meter, the service area of ​​all online metering water meters in the water supply network is obtained, that is, the entire water supply area. The expression of the entire water supply area is: Γ={γ1,γ2,...,γ m } In the formula, γ i is the service area of ​​the online metering water meter, i=1, 2, ..., m, Γ is the set of service areas of all online metering water meters in the water supply network; S12, dividing the entire water supply area obtained in step S11 into grid units to obtain multiple grid units, and the set expression of all grid units is: G={g1,g2,...,g n } In the formula, g j The grid cells divided into the entire water supply area, j = 1, 2, ..., n, the size of the grid cell is determined by the accuracy of the mobile phone signaling data platform, and G is the set of all grid cells; S13. The entire water supply area and multiple grid units are input as input parameters to the mobile signaling data platform to calculate the population in the service area of ​​each online water meter and the population in each grid unit at the current moment.

3. The method for initializing water demand of water supply network nodes based on dynamic population distribution according to claim 2, characterized in that: Step S2 includes the following steps: S21. Intersect any online water meter service area with the set of all grid cells to obtain at least one sub-area. The expression of the sub-area is: In the formula, γ i ∩g j Represents γ i With g j Intersection operation; γ i is the service area of ​​the online water meter, g j is the grid unit divided into the entire water supply area, Λ is γ i With g j The set of sub-regions generated by the intersection operation; K is γ i The number of sub-regions generated after the set G of all grid cells is intersected and divided; λ ik is the sub-region after segmentation; S22, according to the population number in the service area of ​​each online water meter obtained in step S1 at the current moment, combined with the sub-area obtained in step S21, the population number in the sub-area at the current moment is calculated, and the calculation formula is as follows: In the formula, is the sub-region λ ik The current population of the area; The current online water meter service area γ i The population within is the sub-region λ ik area; For γ i the area of ​​the region; S23. According to step S21, all sub-areas in any grid unit are removed to obtain the non-online metered water service area in the grid unit, which is expressed as: g′ j =g j \Λ,g′ j ∈G′ In the formula, g′ j For non-online water metering service area, g i \Λ is g j The difference between Λ and G′ is g′ j A collection of components; S24, according to the population number in each grid unit obtained in step S1, the population number in the sub-area at the current moment obtained in step S22 and the non-online metered water service area obtained in step S23, the population number in the non-online metered water service area in each grid unit at the current moment is obtained, and the calculation formula is as follows: In the formula, is the grid unit g j The population within g j The current population of all sub-areas within; is g′ j The population in the region, g′ j This is a non-online water metering service area.

4. The method for initializing water demand of water supply network nodes based on dynamic population distribution according to claim 3 is characterized in that: Step S3 includes the following steps: S31, taking each water supply network node as the center, draw a Thiessen polygon Z = {z1, z2, ..., z N }; S32, performing an intersection operation on the non-online water metering service area in all grid cells obtained in step S2 and the Thiessen polygon obtained in step S31 to obtain a set of segmented sub-areas, expressed as: In the formula, s lr is g′ j A single sub-region after being divided by Z, g′ j is the non-online water metering service area, G′ is g′ j A collection of components; S33. Based on step S32, the number of people using non-online metered water in the service area of ​​each water supply network node at the current moment is calculated, and the calculation formula is: In the formula, s lr is g′ j A single sub-region after being divided by Z; For sub-region s lr area; is the region g′ j area; S34, based on step S33 and the online water metering data at the current moment, initialize the water demand of the water supply network node at the current moment, and obtain the initial value of the water demand of the water supply network node. The calculation formula is: Where, d j is the initial value of water demand at the water supply network node j at the current moment; d non is the non-metered water demand at the current time node; d tol is the total non-metered water demand of the water supply network at the current moment; P tol is the total population of the entire water supply area at the current moment; ∑d online is the total water demand of all online water meters corresponding to node j at the current moment, P j It is the number of people using non-online metered water in the service area of ​​each water supply network node at the current moment.

Images