Method and system for abnormal monitoring of pipeline network flow based on ultrasonic flowmeter

Abstract

The present invention provides a method and system for monitoring abnormal pipeline flow based on an ultrasonic flowmeter, which relates to the technical field of flow monitoring. The method includes: acquiring first flow data collected by a first ultrasonic flowmeter at the pipeline inlet; acquiring second flow data collected by a second ultrasonic flowmeter at the pipeline outlet; determining the association relationship between the pipeline inlet and the pipeline outlet according to the design parameters of the pipe network; determining the association type between the pipeline inlet and the pipeline outlet with an association relationship; determining the flow response duration according to the first flow data, the second flow data, and the design parameters; determining whether there is an abnormality in the pipeline between the pipeline inlet and the pipeline outlet with an association relationship according to the flow response duration, the design parameters, and the association type. According to the present invention, pipeline abnormalities can be discovered in a timely manner, the health status of the pipe network can be effectively supervised, and potential safety hazards can be reduced.

Description

Technical Field

[0001] The present invention relates to the technical field of flow monitoring, and particularly to a method and system for monitoring abnormal pipeline network flow based on an ultrasonic flowmeter. Background Art

[0002] In the related art, a pipeline network can be composed of multiple pipelines, and the connection relationships between the pipelines are complex. Therefore, in the monitoring and management of the pipeline network, it is difficult to determine whether the flow rates of the pipelines in the pipeline network are abnormal, and it is also difficult to determine the type of abnormality. For example, it is difficult to determine whether there is congestion or leakage in the pipeline. As a result, it is difficult to effectively supervise the health status of the pipeline network, which may lead to potential safety hazards when there is a large flow rate. For example, during heavy rain, the drainage pipeline network may have potential safety hazards due to congestion or leakage. Summary of the Invention

[0003] An embodiment of the present invention provides a method and system for monitoring abnormal pipeline network flow based on an ultrasonic flowmeter, which can timely detect pipeline abnormalities, effectively supervise the health status of the pipeline network, and reduce potential safety hazards.

[0004] According to the first aspect of the embodiments of the present invention, there is provided a method for monitoring abnormal pipeline network flow based on an ultrasonic flowmeter, including:

[0005] Obtaining first flow rate data collected by a first ultrasonic flowmeter disposed at the inlets of multiple pipelines in the pipeline network at multiple moments within a preset time period;

[0006] Obtaining second flow rate data collected by a second ultrasonic flowmeter disposed at the outlets of multiple pipelines in the pipeline network at multiple moments within a preset time period;

[0007] Determining the association relationship between each pipeline inlet and each pipeline outlet according to the design parameters of the pipeline network;

[0008] Determining the association type between the pipeline inlet and the pipeline outlet with an existing association relationship;

[0009] Determining the flow response duration between the pipeline inlet and the pipeline outlet according to the first flow rate data of the pipeline inlet with an existing association relationship, the second flow rate data of the pipeline outlet, and the design parameters of the pipeline network;

[0010] Determining whether there is an abnormality in the pipeline between the pipeline inlet and the pipeline outlet with an existing association relationship according to the flow response duration, the design parameters of the pipeline network, and the association type.

[0011] According to an embodiment of the present invention, determining the flow response duration between the pipeline inlet and the pipeline outlet according to the first flow rate data of the pipeline inlet with an existing association relationship, the second flow rate data of the pipeline outlet, and the design parameters of the pipeline network includes:

[0012] Obtain a first flow vector of the first flow data of the a-th pipeline inlet having an associated relationship and a second flow vector of the second flow data of the b-th pipeline outlet, where a is a positive integer and a is less than or equal to the total number of pipeline inlets, and b is a positive integer and b is less than or equal to the total number of pipeline outlets;

[0013] Determine the discrete shift response duration according to the first flow vector, the second flow vector, and the design parameters of the pipe network;

[0014] Obtain a first flow function of the first flow data according to the first flow data of the a-th pipeline inlet and the design parameters of the pipe network;

[0015] Obtain a second flow function of the second flow data according to the second flow data of the b-th pipeline outlet and the design parameters of the pipe network;

[0016] Obtain the continuous shift response duration according to the first flow function and the second flow function;

[0017] Obtain the flow response duration between the pipeline inlet and the pipeline outlet according to the discrete shift response duration and the continuous shift response duration.

[0018] According to an embodiment of the present invention, determining the discrete shift response duration according to the first flow vector, the second flow vector, and the design parameters of the pipe network includes: according to the formula

[0019]

[0020] Obtain the maximum discrete similarity , where is the forward displacement matrix at the s-th step

[0021]

[0022] is the negative displacement matrix at the s-th step

[0023]

[0024] Indicates that in the case where the number of steps s is the k-th step, = 1, otherwise = 0, both s and k are positive integers less than or equal to n, and n is the number of moments within a preset time period, is the second flow data at the i-th moment of the b-th pipeline outlet, is the cross-sectional area of the b-th pipeline outlet, is the first flow data at the i-th moment of the a-th pipeline inlet, is the cross-sectional area of the a-th pipe inlet, and i is a positive integer less than or equal to n;

[0025] Determine the number of steps corresponding to the maximum discrete similarity;

[0026] According to the formula

[0027]

[0028] Determine the discrete shift response duration , where is the number of steps corresponding to the maximum discrete similarity, is the interval duration between each moment within a preset time period.

[0029] According to an embodiment of the present invention, obtaining the continuous shift response duration according to the first flow function and the second flow function includes:

[0030] According to the formula

[0031] Obtain the minimum continuous deviation , where is the start time of the preset time period, is the end time of the preset time period, t is the time within the preset time period, is the time offset, is the first flow function, is the second flow function;

[0032] Determine the time offset corresponding to the minimum continuous deviation;

[0033] According to the formula

[0034] Determine the continuous shift response duration , where is the time offset corresponding to the minimum continuous deviation.

[0035] According to an embodiment of the present invention, obtaining the flow response duration between the pipe inlet and the pipe outlet according to the discrete shift response duration and the continuous shift response duration includes:

[0036] According to the formula

[0037]

[0038] Obtain the flow response duration , where and are weight parameters and satisfy the formula

[0039] 。

[0040] According to an embodiment of the present invention, the association type includes any one of the following:

[0041] One pipe inlet is only associated with one pipe outlet;

[0042] Two or more pipe inlets are associated with one pipe outlet;

[0043] One pipe inlet is associated with two or more pipe outlets;

[0044] Two or more pipe inlets are associated with two or more pipe outlets.

[0045] According to an embodiment of the present invention, according to the flow response duration, the design parameters of the pipe network, and the association type, determining whether there is an abnormality in the pipe between the pipe inlet and the pipe outlet with an association relationship includes:

[0046] According to the design parameters of the pipe network and the association type, determining the theoretical response duration of the pipe between the pipe inlet and the pipe outlet with an association relationship;

[0047] When the deviation amplitude between the flow response duration and the theoretical response duration is greater than or equal to the amplitude threshold, it is determined that there is an abnormality in the pipe between the pipe inlet and the pipe outlet with an association relationship.

[0048] According to the second aspect of the embodiment of the present invention, there is provided a pipe network flow abnormality monitoring system based on an ultrasonic flowmeter, including:

[0049] A first flow data module, configured to obtain first flow data collected by a first ultrasonic flowmeter provided at a plurality of pipe inlets in the pipe network at a plurality of moments within a preset time period;

[0050] A second flow data module, configured to obtain second flow data collected by a second ultrasonic flowmeter provided at a plurality of pipe outlets in the pipe network at a plurality of moments within a preset time period;

[0051] An association relationship module, configured to determine the association relationship between each pipe inlet and each pipe outlet according to the design parameters of the pipe network;

[0052] An association type module, configured to determine the association type between the pipe inlet and the pipe outlet with an association relationship;

[0053] A flow response duration module, configured to determine the flow response duration between the pipe inlet and the pipe outlet according to the first flow data of the pipe inlet with an association relationship, the second flow data of the pipe outlet, and the design parameters of the pipe network;

[0054] Anomaly judgment module, configured to determine whether there is an anomaly in the pipeline between the pipeline inlet and the pipeline outlet with an associated relationship according to the flow response duration, the design parameters of the pipeline network, and the association type.

[0055] According to an embodiment of the present invention, the flow response duration module is further configured to:

[0056] Obtain a first flow vector of the first flow data of the a-th pipeline inlet with an associated relationship and a second flow vector of the second flow data of the b-th pipeline outlet, where a is a positive integer and a is less than or equal to the total number of pipeline inlets, and b is a positive integer and b is less than or equal to the total number of pipeline outlets;

[0057] Determine the discrete shift response duration according to the first flow vector, the second flow vector, and the design parameters of the pipeline network;

[0058] Obtain a first flow function of the first flow data according to the first flow data of the a-th pipeline inlet and the design parameters of the pipeline network;

[0059] Obtain a second flow function of the second flow data according to the second flow data of the b-th pipeline outlet and the design parameters of the pipeline network;

[0060] Obtain the continuous shift response duration according to the first flow function and the second flow function;

[0061] Obtain the flow response duration between the pipeline inlet and the pipeline outlet according to the discrete shift response duration and the continuous shift response duration.

[0062] According to an embodiment of the present invention, the flow response duration module is further configured to:

[0063] According to the formula

[0064]

[0065] Obtain the maximum discrete similarity , where is the forward displacement matrix at the s-th step

[0066]

[0067] is the negative displacement matrix at the s-th step

[0068]

[0069] Indicates that in the case where the number of steps s is the k-th step, = 1, otherwise = 0, where both s and k are positive integers less than or equal to n, and n is the number of moments within a preset time period. is the second flow rate data at the i-th moment of the b-th pipe outlet. is the cross-sectional area of the b-th pipe outlet. is the first flow rate data at the i-th moment of the a-th pipe inlet. is the cross-sectional area of the a-th pipe inlet, and i is a positive integer less than or equal to n.

[0070] Determine the number of steps corresponding to the maximum discrete similarity.

[0071] According to the formula

[0072]

[0073] Determine the discrete shift response duration , where is the number of steps corresponding to the maximum discrete similarity. is the interval duration between each moment within the preset time period.

[0074] According to an embodiment of the present invention, the flow response duration module is further configured to:

[0075] According to the formula

[0076]

[0077] Obtain the minimum continuous deviation , where is the start moment of the preset time period. is the end moment of the preset time period, t is a moment within the preset time period. is the time offset. is the first flow rate function. is the second flow rate function.

[0078] Determine the time offset corresponding to the minimum continuous deviation.

[0079] According to the formula

[0080]

[0081] Determine the continuous shift response duration , where is the time offset corresponding to the minimum continuous deviation.

[0082] According to an embodiment of the present invention, the flow response duration module is further configured to:

[0083] According to the formula

[0084]

[0085] Obtain the flow response duration , wherein, and are weight parameters and satisfy the formula

[0086] .

[0087] According to an embodiment of the present invention, the association type includes any one of the following:

[0088] One pipeline inlet is only associated with one pipeline outlet;

[0089] Two or more pipeline inlets are associated with one pipeline outlet;

[0090] One pipeline inlet is associated with two or more pipeline outlets;

[0091] Two or more pipeline inlets are associated with two or more pipeline outlets.

[0092] According to an embodiment of the present invention, the abnormality judgment module is further configured to:

[0093] Determine the theoretical response duration of the pipeline between the pipeline inlet and the pipeline outlet with an association relationship according to the design parameters of the pipe network and the association type;

[0094] When the deviation amplitude between the flow response duration and the theoretical response duration is greater than or equal to the amplitude threshold, determine that there is an abnormality in the pipeline between the pipeline inlet and the pipeline outlet with an association relationship.

[0095] According to the third aspect of the embodiments of the present invention, there is provided a pipe network flow abnormality monitoring device based on an ultrasonic flowmeter, including: a processor; a memory for storing instructions executable by the processor; wherein, the processor is configured to call the instructions stored in the memory to execute the pipe network flow abnormality monitoring method based on the ultrasonic flowmeter.

[0096] According to the fourth aspect of the embodiments of the present invention, there is provided a computer-readable storage medium, on which computer program instructions are stored, and when the computer program instructions are executed by a processor, the pipe network flow abnormality monitoring method based on the ultrasonic flowmeter is implemented.

[0097] The method for monitoring abnormal pipeline network flow based on an ultrasonic flowmeter according to an embodiment of the present invention can obtain real-time flow data at multiple pipeline inlets and outlets in the pipeline network through the ultrasonic flowmeter, improving the real-time performance of flow monitoring. Based on the flows at the pipeline inlet and outlet with an associated relationship, the flow response duration between the pipeline inlet and outlet can be determined, thereby determining whether there is an abnormality in the pipeline between the pipeline inlet and outlet, so that pipeline abnormalities can be detected in a timely manner, effectively supervising the health status of the pipeline network and reducing potential safety hazards. When solving for the flow response duration, a flow response duration with both duration and accuracy in terms of a reference discrete shift response duration and a continuous shift response duration can be obtained, improving the accuracy and objectivity of the flow response duration. When solving for the discrete shift response duration, to address the possible difference in cross-sectional areas at the pipeline inlet and outlet, the cross-sectional area factor can be eliminated through the design parameters of the pipeline network, making the flow data at the pipeline inlet and the flow data at the pipeline outlet comparable and improving the accuracy of subsequent calculations. To address the problem of the time difference in the changes in flow data at the inlet and outlet caused by the same liquid flow, the maximum discrete similarity can be found through shifting, and the displacement speed can be increased through two displacement matrices, namely a positive displacement matrix and a negative displacement matrix, improving the efficiency of finding the maximum discrete similarity. Furthermore, based on the number of steps of displacement when the maximum discrete similarity is obtained, the discrete shift response duration can be determined, thereby obtaining the time difference in the changes in flow data at the inlet and outlet caused by the same liquid flow and improving the accuracy of calculating the time difference. When solving for the continuous shift response duration, by offsetting the first flow function and the second flow function, the minimum deviation of the integral value of the two functions after translation is found, and then the time offset corresponding to the minimum continuous deviation and the continuous shift response duration are determined. At the same time, offsetting the two flow functions can improve the efficiency of finding the minimum deviation of the integral value and improve the accuracy of the continuous shift response duration.

[0098] It should be understood that the above general description and the following detailed description are merely exemplary and explanatory, and do not limit the present invention. Other features and aspects of the present invention will become clearer based on the following detailed description of exemplary embodiments with reference to the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0099] To more clearly illustrate the technical solutions in the embodiments of the present invention or in the prior art, the following will briefly introduce the drawings required for use in the description of the embodiments or the prior art. Obviously, the following drawings are only some embodiments of the present invention. For those of ordinary skill in the art, other embodiments can be obtained based on these drawings without creative efforts.

[0100] Figure 1 Exemplarily shows a schematic flowchart of a method for monitoring abnormal pipeline network flow based on an ultrasonic flowmeter according to an embodiment of the present invention;

[0101] Figure 2 Exemplarily shown is a schematic diagram of a pipeline network flow anomaly monitoring system based on an ultrasonic flowmeter according to an embodiment of the present invention. Detailed implementation manners

[0102] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Apparently, the described embodiments are only some of the embodiments of the present invention, rather than all of the embodiments. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0103] The technical solutions of the present invention will be described in detail below with specific embodiments. These specific embodiments may be combined with each other, and the same or similar concepts or processes may not be repeated in some embodiments.

[0104] Figure 1 Exemplarily shown is a schematic flowchart of a pipeline network flow anomaly monitoring method based on an ultrasonic flowmeter according to an embodiment of the present invention, as Figure 1 shown, the method includes:

[0105] Step S101, obtaining first flow data collected by a first ultrasonic flowmeter provided at multiple pipeline inlets in the pipeline network at multiple moments within a preset time period;

[0106] Step S102, obtaining second flow data collected by a second ultrasonic flowmeter provided at multiple pipeline outlets in the pipeline network at multiple moments within a preset time period;

[0107] Step S103, determining the association relationship between each pipeline inlet and each pipeline outlet according to the design parameters of the pipeline network;

[0108] Step S104, determining the association type between a pipeline inlet and a pipeline outlet with an association relationship;

[0109] Step S105, determining the flow response duration between a pipeline inlet and a pipeline outlet according to the first flow data of the pipeline inlet with an association relationship, the second flow data of the pipeline outlet, and the design parameters of the pipeline network;

[0110] Step S106, determining whether there is an anomaly in the pipeline between a pipeline inlet and a pipeline outlet with an association relationship according to the flow response duration, the design parameters of the pipeline network, and the association type.

[0111] The abnormal monitoring method of pipeline network flow based on ultrasonic flowmeter according to an embodiment of the present invention can obtain real-time flow data at multiple pipeline inlets and outlets in the pipeline network through the ultrasonic flowmeter, improving the real-time performance of flow monitoring. And based on the flow rates at the pipeline inlets and outlets with associated relationships, the flow response duration between the pipeline inlet and the pipeline outlet can be determined, so as to determine whether there is an abnormality in the pipeline between the pipeline inlet and the pipeline outlet, so that pipeline abnormalities can be discovered in time, the health status of the pipeline network can be effectively supervised, and potential safety hazards can be reduced.

[0112] According to an embodiment of the present invention, the pipeline network may include multiple pipelines, which are intricate and may have situations where multiple branch pipelines converge into the same pipeline, or a pipeline may branch out into multiple branch pipelines. Of course, it is also possible that both of the above situations exist. Therefore, it is complex to judge whether the flow in the pipeline network is abnormal, and it is also complex to judge whether a certain measured flow data indicates that there are faults such as blockage or leakage in a certain or certain pipelines.

[0113] According to an embodiment of the present invention, in this case, first ultrasonic flowmeters may be respectively arranged at multiple pipeline inlets in the pipeline network, and second ultrasonic flowmeters may be respectively arranged at multiple pipeline outlets in the pipeline network, so as to judge whether the flow of a certain or certain pipelines is abnormal based on the flow data measured by the first ultrasonic flowmeters and the flow data measured by the second ultrasonic flowmeters, and whether there are faults such as blockage or leakage in the pipeline can be predicted.

[0114] According to an embodiment of the present invention, in step S101, the first flow data of the first ultrasonic flowmeter may be collected at multiple moments within a preset time period. The preset time period may be a time period with a duration of one hour, three hours, one day, etc. The present invention does not limit the duration of the preset time period. The preset time period may include multiple moments. For example, each minute may be used as a moment, and every 3 minutes may be used as a moment, etc. The present invention does not limit the time interval between moments.

[0115] According to an embodiment of the present invention, the first flow data of the first ultrasonic flowmeter at multiple pipeline inlets with non-zero flow rates may be collected within a preset time period. For example, there are 100 pipeline inlets of the sewer pipelines in a certain area, and 60 pipeline inlets with non-zero flow rates within the preset time period. Then, the first flow data of the first ultrasonic flowmeters at these 60 pipeline inlets within the preset time period can be obtained. The first flow data collected by the other 40 first ultrasonic flowmeters within the preset time period are all 0, and these all-zero data can be discarded.

[0116] According to an embodiment of the present invention, in step S102, the second flow data collected by the second ultrasonic flowmeters arranged at multiple pipeline outlets in the pipeline network at multiple moments may be obtained within a preset time period.

[0117] According to an embodiment of the present invention, in step S103, design parameters of the pipe network can be obtained. The design parameters may include size parameters of each pipeline in the pipe network, such as cross-sectional area, length and other parameters, and may also include the association relationship between each inlet and outlet. In the example, it can be determined which or which pipeline outlets the flow at a certain pipeline inlet finally flows to.

[0118] According to an embodiment of the present invention, the second flow data in step S102 can also be filtered based on the association relationship. For example, the second flow data of the pipeline outlets associated with the pipeline inlets with non-zero flow within a preset time period can be filtered out, and the second flow data at other pipeline outlets can be discarded.

[0119] According to an embodiment of the present invention, in step S104, the association type between the pipeline inlet and the pipeline outlet with an association relationship can be determined. For example, if a pipeline corresponding to a certain pipeline inlet can be divided into multiple pipelines, then the pipeline inlet can correspond to the pipeline outlets of the divided multiple pipelines. For another example, if the pipelines corresponding to multiple pipeline inlets merge into one pipeline, then the multiple pipeline inlets correspond to the pipeline outlet of the merged pipeline.

[0120] According to an embodiment of the present invention, the association type includes any one of the following: one pipeline inlet is only associated with one pipeline outlet; two or more pipeline inlets are associated with one pipeline outlet; one pipeline inlet is associated with two or more pipeline outlets; two or more pipeline inlets are associated with two or more pipeline outlets.

[0121] According to an embodiment of the present invention, the association type in which one pipeline inlet is only associated with one pipeline outlet means that the pipeline inlet and the pipeline outlet are only connected by one pipeline, and the flow at the pipeline inlet directly flows to the pipeline outlet, and there is no situation where other pipelines are merged into this pipeline or this pipeline merges into other pipelines. The cross-sectional area of this pipeline can remain unchanged or can also change, and the present invention does not limit this.

[0122] According to an embodiment of the present invention, the association type in which two or more pipeline inlets are associated with one pipeline outlet means that the pipelines at two or more pipeline inlets merge into one pipeline. Therefore, two or more pipeline inlets are associated with the pipeline outlet of the merged one pipeline.

[0123] According to an embodiment of the present invention, the association type in which one pipeline inlet is associated with two or more pipeline outlets means that the pipeline at one pipeline inlet is divided into multiple branch pipelines. Therefore, this pipeline inlet is associated with the pipeline outlets of the multiple branch pipelines.

[0124] According to an embodiment of the present invention, the association type in which two or more pipe inlets are associated with two or more pipe outlets means that the pipes at the two or more pipe inlets merge and branch twice or more times to reach the two or more pipe outlets. Therefore, the above two or more pipe inlets are associated with the above two or more pipe outlets.

[0125] According to an embodiment of the present invention, after determining the above association relationship and association type, the flow response duration between the pipe inlet and the pipe outlet with the association relationship can be determined, that is, the duration for the same liquid flow to flow from the pipe inlet to the pipe outlet. Furthermore, this duration can be used as a basis for determining whether there is an abnormality in the pipe between the pipe inlet and the pipe outlet.

[0126] According to an embodiment of the present invention, in step S105, the flow response duration between the pipe inlet and the pipe outlet can be determined based on the first flow data of the pipe inlet with the association relationship, the second flow data of the pipe outlet, and the design parameters of the pipe network (such as cross-sectional area, etc.). Step S105 may include: obtaining a first flow vector of the first flow data of the a-th pipe inlet with the association relationship and a second flow vector of the second flow data of the b-th pipe outlet, where a is a positive integer and a is less than or equal to the total number of pipe inlets, and b is a positive integer and b is less than or equal to the total number of pipe outlets; determining the discrete shift response duration according to the first flow vector, the second flow vector, and the design parameters of the pipe network; obtaining a first flow function of the first flow data according to the first flow data of the a-th pipe inlet and the design parameters of the pipe network; obtaining a second flow function of the second flow data according to the second flow data of the b-th pipe outlet and the design parameters of the pipe network; obtaining the continuous shift response duration according to the first flow function and the second flow function; and obtaining the flow response duration between the pipe inlet and the pipe outlet according to the discrete shift response duration and the continuous shift response duration.

[0127] According to an embodiment of the present invention, there is an association relationship between the a-th pipe inlet and the b-th pipe outlet. At multiple moments within a preset time period, the first ultrasonic flowmeter at the a-th pipe inlet can obtain multiple first flow data, and each first flow data can be used as a vector element to obtain a high-dimensional first flow vector. Similarly, at multiple moments within a preset time period, the second ultrasonic flowmeter at the b-th pipe outlet can obtain multiple second flow data, and each second flow data can be used as a vector element to obtain a high-dimensional second flow vector. The dimensions of the first flow vector and the second flow vector are the same and are both equal to the number of moments within the preset time period.

[0128] According to an embodiment of the present invention, the elements in the above first flow vector and second flow vector are flow rate data at discrete moments. Therefore, the discrete shift response duration can be obtained from the perspective of discrete data, that is, the duration for the same liquid flow to flow from the a-th pipe inlet to the b-th pipe outlet determined based on the above discrete data.

[0129] According to an embodiment of the present invention, determining the discrete shift response duration according to the first flow vector, the second flow vector and the design parameters of the pipe network includes:

[0130] Obtain the maximum discrete similarity according to formula (1) ,

[0131] (1)

[0132] Where, is the forward displacement matrix at the s-th step, as shown in the following formula (2),

[0133] (2)

[0134] is the negative displacement matrix at the s-th step, as shown in the following formula (3),

[0135] (3)

[0136] Indicates that in the case where the number of steps s is the k-th step, = 1, otherwise = 0, both s and k are positive integers less than or equal to n, and n is the number of moments within a preset time period, is the second flow rate data at the i-th moment of the b-th pipe outlet, is the cross-sectional area of the b-th pipe outlet, is the first flow rate data at the i-th moment of the a-th pipe inlet, is the cross-sectional area of the a-th pipe inlet, and i is a positive integer less than or equal to n;

[0137] Determine the number of steps corresponding to the maximum discrete similarity;

[0138] Determine the discrete shift response duration according to formula (4) ,

[0139] (4)

[0140] Where, is the number of steps corresponding to the maximum discrete similarity, is the interval duration between each moment within the preset time period.

[0141] According to an embodiment of the present invention, in formula (1), since the cross-sectional areas of the a-th pipe inlet and the b-th pipe outlet may be different, in order to make the first flow rate vector and the second flow rate vector comparable, each element of the first flow rate vector can be divided by the cross-sectional area of the a-th pipe inlet to obtain the first flow rate vector after removing the area factor. Similarly, each element of the second flow rate vector is divided by the cross-sectional area of the b-th pipe outlet to obtain the second flow rate vector after removing the area factor. .

[0142] According to an embodiment of the present invention, although the first flow rate vector and the second flow rate vector exclude the cross-sectional area factor, the two obtained vectors are still not similar. That is, since there is a certain distance between the inlet and the outlet of the pipe, and it takes a certain amount of time for the liquid flow to pass through this distance. Therefore, there is a certain time difference when the same liquid flow passes through the a-th pipe inlet and the b-th pipe outlet. For example, when the liquid flow passes through the a-th pipe inlet, the first flow rate data increases. When this liquid flow passes through the b-th pipe outlet, the second flow rate data increases. However, at this time, the first flow rate data has decreased, resulting in a relatively low correlation between the first flow rate data and the second flow rate data at the same moment, but a relatively high correlation between the first flow rate data and the second flow rate data at different moments.

[0143] Based on the above analysis, according to an embodiment of the present invention, if the similarity between the first flow rate vector after removing the area factor and the second flow rate vector after removing the area factor is directly determined according to the corresponding elements, the accuracy of the obtained similarity is not high, and the above time difference factor cannot be considered, and the above time difference cannot be solved. Therefore, the forward shift matrix and the negative shift matrix can be used to respectively perform forward shift on the first flow rate vector after removing the area factor and perform negative shift on the second flow rate vector after removing the area factor. And after each displacement, the cosine similarity between the first flow rate vector after removing the area factor after displacement and the second flow rate vector after removing the area factor after displacement is solved. After obtaining the cosine similarities of multiple displacements, the maximum value of the obtained cosine similarities is determined as the maximum discrete similarity.

[0144] According to an embodiment of the present invention, during the shifting process, since the same liquid flow passes through the a-th pipe inlet first, the first flow rate data at the a-th pipe inlet changes first. After a certain period of time, the second flow rate data at the b-th pipe outlet will change. Therefore, in order to align the first flow rate data and the second flow rate data changed by the same liquid flow during the process of solving the similarity, the first flow rate data can be shifted backward, and after shifting, the data before the first bit is filled with 0. Or the second flow rate data can be shifted forward, and after shifting, the data after the last bit is filled with 0. Of course, in order to find the maximum value of the cosine similarity faster and solve the time difference for the same liquid flow to flow from the a-th pipe inlet to the b-th pipe outlet faster, the above backward and forward shifts can be performed simultaneously to improve the solving efficiency.

[0145] In the example, the elements of the first row of the forward displacement matrix are: the element in the first row and the first column is , which means that when the number of displacement steps s is 0, that is, when no displacement is performed, this element is 1, otherwise it is 0, and the other elements in the first row are all 0. The elements of the second row of the forward displacement matrix are: the element in the second row and the first column is , which means that when the number of displacement steps s is 1, this element is 1, otherwise it is 0, and the element in the second row and the second column is , and the rest of the elements are 0. The elements of the third row of the forward displacement matrix are: the element in the third row and the first column is , which means that when the number of displacement steps s is 2, this element is 1, otherwise it is 0, the element in the third row and the second column is , the element in the third row and the third column is , and the rest of the elements are 0. And so on, the elements of the n-th row of the forward displacement matrix are: the element in the n-th row and the first column is , which means that when the number of displacement steps s is n - 1, this element is 1, otherwise it is 0, the element in the n-th row and the second column is , which means that when the number of displacement steps s is n - 2, this element is 1, otherwise it is 0... the element in the n-th row and the n-th column is .

[0146] In the example, multiply the first flow rate vector after removing the area factor by this forward displacement matrix. When the number of displacement steps s is 0, this vector remains unchanged. When the number of displacement steps s is 1, we get , when the number of displacement steps s is 2, we get ... when the number of displacement steps s is n - 1, we get .

[0147] In the example, the elements of the first row of the negative displacement matrix are: the element in the first row and the first column is , which means that when the number of displacement steps s is 0, that is, when no displacement is performed, this element is 1, otherwise it is 0, and the element in the first row and the second column is which means that when the number of steps s of the displacement is 1, this element is 1, otherwise it is 0. The element in the first row and the third column is which means that when the number of steps s of the displacement is 2, this element is 1, otherwise it is 0... The element in the first row and the (n - 1)th column is which means that when the number of steps s of the displacement is n - 2, this element is 1, otherwise it is 0. The element in the first row and the nth column is which means that when the number of steps s of the displacement is n - 1, this element is 1, otherwise it is 0. The elements in the second row of the negative displacement matrix are: the element in the second row and the first column is 0, the element in the second row and the second column is which means that when the number of steps s of the displacement is 0, this element is 1, otherwise it is 0. The element in the second row and the third column is which means that when the number of steps s of the displacement is 1, this element is 1, otherwise it is 0... The element in the second row and the nth column is which means that when the number of steps s of the displacement is n - 2, this element is 1, otherwise it is 0. By analogy, the elements in the nth row of the negative displacement matrix are: the element in the nth row and the first column is 0, the element in the nth row and the second column is 0... the element in the nth row and the (n - 1)th column is 0, and the element in the nth row and the nth column is which means that when the number of steps s of the displacement is 0, this element is 1, otherwise it is 0.

[0148] In the example, the negative displacement matrix is used to multiply the second flow vector after removing the area factor When the number of steps s of the displacement is 0, this vector remains unchanged. When the number of steps s of the displacement is 1, we get When the number of steps s of the displacement is 2, we get ... When the number of steps s of the displacement is n - 1, we get .

[0149] According to an embodiment of the present invention, the cosine similarity of the vectors obtained after each step of displacement above can be solved. The maximum value of the cosine similarity is the maximum discrete similarity, and when obtaining the maximum value of the cosine similarity, the number of steps of the displacement is the number of steps corresponding to the maximum discrete similarity.

[0150] According to an embodiment of the present invention, in order to improve the search efficiency of the maximum value of the cosine similarity, the forward displacement matrix and the negative displacement matrix are used simultaneously to displace the first flow vector after removing the area factor and the second flow vector after removing the area factor respectively. In each displacement step, the first flow vector after removing the area factor is displaced backward by one step, and the second flow vector after removing the area factor is displaced forward by one step, so the relative displacement between the two is two steps. Moreover, the maximum discrete similarity can represent the similarity of the flow data changes caused by the same liquid flow at the inlet of the a-th pipe and the flow data changes caused at the outlet of the b-th pipe. Therefore, the number of steps corresponding to the maximum discrete similarity is the number of steps when the first flow vector after removing the area factor and the second flow vector after removing the area factor are displaced to place the flow data changes caused by the same liquid flow at the same position. This number of steps is half of the number of steps difference between the flow data changes caused by the same liquid flow at the inlet of the a-th pipe and the flow data changes caused at the outlet of the b-th pipe. Therefore, in formula (4), 2 times of represents the number of steps difference between the flow data changes caused by the same liquid flow at the inlet of the a-th pipe and the flow data changes caused at the outlet of the b-th pipe. By multiplying the number of steps difference by the time interval between each moment, the discrete shift response duration representing the time difference between the flow data changes caused by the same liquid flow at the inlet of the a-th pipe and the flow data changes caused at the outlet of the b-th pipe can be obtained.

[0151] In this way, for the problem that the cross-sectional areas at the pipe inlet and the pipe outlet may be different, the area factor can be excluded through the design parameters of the pipe network, making the flow data at the pipe inlet and the flow data at the pipe outlet comparable and improving the accuracy of subsequent calculations. And for the problem that there is a time difference in the flow data changes at the inlet and the outlet caused by the same liquid flow, the maximum discrete similarity can be found through displacement, and the displacement speed can be increased by using two displacement matrices, namely the forward displacement matrix and the negative displacement matrix, to improve the efficiency of finding the maximum discrete similarity. Furthermore, based on the number of steps of displacement when the maximum discrete similarity is obtained, the discrete shift response duration is determined, so as to obtain the time difference between the flow data changes at the inlet and the outlet caused by the same liquid flow, improving the accuracy of calculating the time difference.

[0152] According to an embodiment of the present invention, the discrete shift response duration is obtained based on the discrete flow data at each moment above. Below, the continuous shift response duration can also be obtained based on the continuous flow data, providing reference data for calculating the flow response duration and improving the accuracy and objectivity of the flow response duration.

[0153] According to an embodiment of the present invention, continuous flow rate data can be obtained first. For example, the first flow rate data at the inlet of the a-th pipeline can be fitted, and the second flow rate data at the outlet of the b-th pipeline can be fitted. Moreover, in order to exclude the influence caused by different cross-sectional areas, the factor of the cross-sectional area can be directly excluded during fitting, that is, the ratio of the first flow rate data to the cross-sectional area of the inlet of the a-th pipeline can be fitted to obtain the first flow rate function. Similarly, the ratio of the second flow rate data to the cross-sectional area of the outlet of the b-th pipeline can be fitted to obtain the second flow rate function.

[0154] According to an embodiment of the present invention, based on the first flow rate function and the second flow rate function, a continuous shift response duration is obtained, including: obtaining a minimum continuous deviation according to formula (5) ,

[0155] (5)

[0156] wherein, is the start time of a preset time period, is the end time of the preset time period, t is the time within the preset time period, is the time offset, is the first flow rate function, is the second flow rate function;

[0157] Determine the time offset corresponding to the minimum continuous deviation;

[0158] According to formula (6), determine the continuous shift response duration , (6)

[0159] wherein, is the time offset corresponding to the minimum continuous deviation.

[0160] According to an embodiment of the present invention, the first flow rate function is a relationship function between the first flow rate data and time, and the second flow rate function is a relationship function between the second flow rate data and time. As described above, the change time of the first flow rate data caused by the same liquid flow is different from the change time of the second flow rate data. Therefore, the waveforms of the two functions at the same time are not similar, and it is meaningless to directly subtract them or subtract the integral values. Therefore, and The integral values of two functions are subtracted and compared. That is, the function after the first flow function is shifted to the right by a time offset is compared with the function after the second flow function is shifted to the left by the time offset, and the integral intervals of the two shifted functions are adjusted synchronously. During the shifting process, the minimum difference in the integral values of the shifted functions can be solved as the minimum continuous deviation, that is, the deviation between the waveforms of the two functions after aligning the change times of the first flow data and the second flow data caused by the same liquid flow, and the time offset when the deviation between the integral values of the two functions reaches the minimum can be determined, that is, the time offset corresponding to the minimum continuous deviation.

[0161] According to an embodiment of the present invention, since the time offset corresponding to the minimum continuous deviation is the time offset for shifting the above two functions simultaneously, the time difference between the change time of the first flow data and the change time of the second flow data caused by the same liquid flow is twice this time offset. Therefore, in formula (6), the continuous shift response duration is equal to twice the time offset corresponding to the minimum continuous deviation.

[0162] In this way, for the problem that the cross-sectional areas at the pipe inlet and the pipe outlet may be different, the cross-sectional area factor can be eliminated through the design parameters of the pipe network, making the flow data at the pipe inlet and the flow data at the pipe outlet comparable and improving the accuracy of subsequent calculations. And for the problem that there is a time difference in the changes of the flow data at the inlet and the outlet caused by the same liquid flow, by shifting the first flow function and the second flow function, the minimum deviation of the integral values of the two shifted functions can be found, and then the time offset corresponding to the minimum continuous deviation and the continuous shift response duration can be determined. Shifting the two flow functions simultaneously can improve the efficiency of finding the minimum deviation of the integral values, thereby obtaining the time difference in the changes of the flow data at the inlet and the outlet caused by the same liquid flow and improving the accuracy of calculating the time difference.

[0163] According to an embodiment of the present invention, the discrete shift response duration and the continuous shift response duration are respectively obtained based on the discrete flow data and the continuous flow data, and the flow response duration can be determined based on the two durations. According to the discrete shift response duration and the continuous shift response duration, obtaining the flow response duration between the pipe inlet and the pipe outlet includes:

[0164] Obtain the flow response duration according to formula (7) ,

[0165] (7)

[0166] Wherein, and is a weight parameter and satisfies formula (8).

[0167] (8)

[0168] According to an embodiment of the present invention, based on formula (7), the discrete shift response duration and the continuous shift response duration can be weighted and summed to obtain the flow response duration. The weights are respectively and . Among them, the maximum discrete similarity can represent the similarity between the first flow data after removing the cross-sectional area factor and the second flow data after removing the cross-sectional area factor caused by the same liquid flow, and can also represent the accuracy of the discrete shift response duration. Therefore, the maximum discrete similarity can be used as one of the coefficients in the weight of the discrete shift response duration. The minimum continuous deviation can represent the minimum value of the integral value difference between the first flow function and the second flow function caused by the same liquid flow, and is inversely related to the accuracy of the continuous shift response duration. Therefore, the minimum continuous deviation can be normalized by times the arctangent function, and a parameter positively correlated with the accuracy of the continuous shift response duration can be obtained by subtracting the normalization result from 1. This parameter is used as one of the coefficients in the weight of the continuous shift response duration. Further, the preset weight parameters and can also be used to make the sum of the weights of the discrete shift response duration and the continuous shift response duration equal to 1, that is, satisfy formula (8). Based on the above weighted summation process, the flow response duration considering the duration and accuracy of the reference discrete shift response duration and the continuous shift response duration can be obtained, improving the accuracy and objectivity of the flow response duration.

[0169] According to an embodiment of the present invention, in step 106, based on the flow response duration obtained above, as well as the design parameters and association types of the pipe network, it can be determined whether there is an abnormality in the pipeline between the pipe inlet and the pipe outlet with an associated relationship. For example, it can be determined whether there is an abnormality in the pipeline between the a-th pipe inlet and the b-th pipe outlet.

[0170] According to an embodiment of the present invention, step S106 may include: determining the theoretical response duration of the pipeline between the pipe inlet and the pipe outlet with an associated relationship according to the design parameters and the association type of the pipe network; determining that there is an abnormality in the pipeline between the pipe inlet and the pipe outlet with an associated relationship when the deviation amplitude between the flow response duration and the theoretical response duration is greater than or equal to the amplitude threshold.

[0171] According to an embodiment of the present invention, the theoretical response duration of the pipeline between the pipe inlet and the pipe outlet with an associated relationship can be determined based on the design parameters and association type of the pipe network.

[0172] In the example, the length and cross-sectional area of each pipeline can be determined based on the design parameters of the pipeline network. For a pipeline with a longer length, since the liquid flow needs to pass through a longer distance, its theoretical response time is also longer. At the same flow rate, for a pipeline with a smaller cross-sectional area, the theoretical response time is smaller because the flow rate represents the volume of the liquid flow passing through. If the cross-sectional area is smaller, the flow velocity will be faster, resulting in a shorter theoretical response time.

[0173] In the example, the theoretical response time may be related to the association type. For example, if the association type between the pipeline inlet and the pipeline outlet with an association relationship is a type where one pipeline inlet is only associated with one pipeline outlet, then the theoretical response time is only related to the design parameters of the pipeline network described above.

[0174] In the example, if the association type between the pipeline inlet and the pipeline outlet with an association relationship is a relationship where two or more pipeline inlets are associated with one pipeline outlet, then the pipeline outlet is affected by the liquid flows from multiple pipeline inlets. Considering the resistance of the liquid when multiple pipelines converge, the theoretical response time between the pipeline inlet and the pipeline outlet is longer than the case where one pipeline inlet is only associated with one pipeline outlet.

[0175] In the example, if the association type between the pipeline inlet and the pipeline outlet with an association relationship is a relationship where one pipeline inlet is associated with two or more pipeline outlets, then the liquid flow flowing into the pipeline inlet is divided into flows out from multiple pipeline outlets. In this case, considering that the diversion effect of multiple pipeline outlets reduces the water pressure of the main pipeline corresponding to the pipeline inlet, therefore, the theoretical response time between the pipeline inlet and the pipeline outlet is longer than the case where one pipeline inlet is only associated with one pipeline outlet.

[0176] In the example, if the association type between the pipeline inlet and the pipeline outlet with an association relationship is a relationship where two or more pipeline inlets are associated with two or more pipeline outlets, then there are two or more cases of pipeline convergence and pipeline branching between the pipeline inlet and the pipeline outlet. Therefore, the theoretical response time between the pipeline inlet and the pipeline outlet is longer than the case where one pipeline inlet is only associated with one pipeline outlet.

[0177] According to an embodiment of the present invention, considering the above factors, the theoretical response time can be obtained by conducting tests when no faults occur in the pipeline, or the theoretical response time can be obtained by simulating the liquid flow based on the above factors. The present invention does not limit the method for obtaining the theoretical response time.

[0178] According to an embodiment of the present invention, after obtaining the theoretical response duration, it is possible to determine whether the deviation amplitude between the flow response duration and the theoretical response duration is greater than or equal to the amplitude threshold. For example, the absolute value of the difference between the flow response duration and the theoretical response duration can be determined and divided by the theoretical response duration to obtain the deviation amplitude. If the deviation amplitude is too large, there is an abnormality in the pipeline between the pipeline inlet and the pipeline outlet with an associated relationship. For example, there is a pipeline blockage fault, resulting in a reduction in the cross-sectional area at the pipeline blockage, an increase in the flow velocity, and a significant shortening of the flow response duration compared to the theoretical response duration. Or, there is a pipeline leakage fault, resulting in a decrease in the flow rate and a decrease in the flow velocity, making the flow response duration significantly longer than the theoretical response duration, etc. The present invention does not limit the basis for determining abnormalities and the types of abnormalities.

[0179] For the method for monitoring abnormal pipeline network flow based on an ultrasonic flowmeter according to an embodiment of the present invention, real-time flow data of multiple pipeline inlets and outlets in the pipeline network can be obtained through the ultrasonic flowmeter, improving the real-time performance of flow monitoring. And based on the flow rates of the pipeline inlet and the pipeline outlet with an associated relationship, the flow response duration between the pipeline inlet and the pipeline outlet can be determined, thereby determining whether there is an abnormality in the pipeline between the pipeline inlet and the pipeline outlet, so that pipeline abnormalities can be detected in a timely manner, effectively supervising the health status of the pipeline network and reducing potential safety hazards. When solving the flow response duration, a flow response duration with both duration and accuracy in terms of the reference discrete shift response duration and the continuous shift response duration can be obtained, improving the accuracy and objectivity of the flow response duration. When solving the discrete shift response duration, for the problem that the cross-sectional areas at the pipeline inlet and the pipeline outlet may be different, the cross-sectional area factor can be excluded through the design parameters of the pipeline network, making the flow data at the pipeline inlet and the flow data at the pipeline outlet comparable and improving the accuracy of subsequent calculations. And for the problem that there is a time difference in the changes of the flow data at the inlet and the outlet caused by the same liquid flow, the maximum discrete similarity can be found through shifting, and the displacement speed can be increased through two displacement matrices, namely the positive displacement matrix and the negative displacement matrix, improving the efficiency of finding the maximum discrete similarity. Furthermore, based on the number of steps of displacement when obtaining the maximum discrete similarity, the discrete shift response duration can be determined, thereby obtaining the time difference in the changes of the flow data at the inlet and the outlet caused by the same liquid flow and improving the accuracy of calculating the time difference. When solving the continuous shift response duration, by offsetting the first flow function and the second flow function, the minimum deviation of the integral value of the two functions after translation is found, and then the time offset corresponding to the minimum continuous deviation and the continuous shift response duration are determined. At the same time, offsetting the two flow functions can improve the efficiency of finding the minimum deviation of the integral value and improve the accuracy of the continuous shift response duration.

[0180] Figure 2 Exemplarily shown is a schematic diagram of a pipeline network flow anomaly monitoring system based on an ultrasonic flowmeter according to an embodiment of the present invention, asFigure 2 As shown, the system includes:

[0181] A first flow data module 101, configured to obtain first flow data collected by first ultrasonic flowmeters disposed at multiple pipeline inlets in a pipeline network at multiple moments within a preset time period;

[0182] A second flow data module 102, configured to obtain second flow data collected by second ultrasonic flowmeters disposed at multiple pipeline outlets in the pipeline network at multiple moments within the preset time period;

[0183] An association relationship module 103, configured to determine an association relationship between each pipeline inlet and each pipeline outlet according to design parameters of the pipeline network;

[0184] An association type module 104, configured to determine an association type between a pipeline inlet and a pipeline outlet having an association relationship;

[0185] A flow response duration module 105, configured to determine a flow response duration between a pipeline inlet and a pipeline outlet according to first flow data of a pipeline inlet having an association relationship, second flow data of a pipeline outlet, and design parameters of the pipeline network;

[0186] An abnormality determination module 106, configured to determine whether there is an abnormality in a pipeline between a pipeline inlet and a pipeline outlet having an association relationship according to the flow response duration, design parameters of the pipeline network, and the association type.

[0187] According to an embodiment of the present invention, the flow response duration module is further configured to:

[0188] Obtain a first flow vector of first flow data of the a-th pipeline inlet having an association relationship and a second flow vector of second flow data of the b-th pipeline outlet, where a is a positive integer and a is less than or equal to the total number of pipeline inlets, and b is a positive integer and b is less than or equal to the total number of pipeline outlets;

[0189] Determine a discrete shift response duration according to the first flow vector, the second flow vector, and design parameters of the pipeline network;

[0190] Obtain a first flow function of the first flow data according to the first flow data of the a-th pipeline inlet and design parameters of the pipeline network;

[0191] Obtain a second flow function of the second flow data according to the second flow data of the b-th pipeline outlet and design parameters of the pipeline network;

[0192] Obtain a continuous shift response duration according to the first flow function and the second flow function;

[0193] Obtain the flow response duration between the pipeline inlet and the pipeline outlet according to the discrete shift response duration and the continuous shift response duration.

[0194] According to an embodiment of the present invention, the flow response duration module is further configured to:

[0195] According to the formula

[0196]

[0197] Obtain the maximum discrete similarity , where is the forward displacement matrix at the s-th step

[0198]

[0199] is the negative displacement matrix at the s-th step

[0200]

[0201] indicates that when the number of steps s is the k-th step, = 1, otherwise = 0, both s and k are positive integers less than or equal to n, and n is the number of moments within a preset time period, is the second flow data at the i-th moment of the b-th pipeline outlet, is the cross-sectional area of the b-th pipeline outlet, is the first flow data at the i-th moment of the a-th pipeline inlet, is the cross-sectional area of the a-th pipeline inlet, and i is a positive integer less than or equal to n;

[0202] Determine the number of steps corresponding to the maximum discrete similarity;

[0203] According to the formula

[0204]

[0205] Determine the discrete shift response duration , where is the number of steps corresponding to the maximum discrete similarity, is the interval duration between each moment within a preset time period.

[0206] According to an embodiment of the present invention, the flow response duration module is further configured to:

[0207] According to the formula

[0208]

[0209] Obtain the minimum continuous deviation , where is the start time of the preset time period, is the end time of the preset time period, t is the time within the preset time period, is the time offset, is the first flow function, is the second flow function;

[0210] Determine the time offset corresponding to the minimum continuous deviation;

[0211] According to the formula

[0212]

[0213] Determine the continuous shift response duration , where is the time offset corresponding to the minimum continuous deviation.

[0214] According to an embodiment of the present invention, the flow response duration module is further configured to:

[0215] According to the formula

[0216]

[0217] Obtain the flow response duration , where and are weight parameters and satisfy the formula

[0218] .

[0219] According to an embodiment of the present invention, the association type includes any one of the following:

[0220] One pipeline inlet is only associated with one pipeline outlet;

[0221] Two or more pipeline inlets are associated with one pipeline outlet;

[0222] One pipeline inlet is associated with two or more pipeline outlets;

[0223] Two or more pipeline inlets are associated with two or more pipeline outlets.

[0224] According to an embodiment of the present invention, the abnormal judgment module is further configured to:

[0225] According to the design parameters of the pipe network and the association type, determine the theoretical response duration of the pipeline between the pipeline inlet and the pipeline outlet with an association relationship;

[0226] When the deviation amplitude between the flow response duration and the theoretical response duration is greater than or equal to the amplitude threshold, it is determined that there is an abnormality in the pipeline between the pipeline inlet and the pipeline outlet where the correlation exists.

[0227] According to an embodiment of the present invention, there is provided a pipeline network flow anomaly monitoring device based on an ultrasonic flowmeter, including: a processor; a memory for storing instructions executable by the processor; wherein, the processor is configured to call the instructions stored in the memory to execute the pipeline network flow anomaly monitoring method based on the ultrasonic flowmeter.

[0228] According to an embodiment of the present invention, there is provided a computer-readable storage medium, on which computer program instructions are stored, and when the computer program instructions are executed by a processor, the pipeline network flow anomaly monitoring method based on the ultrasonic flowmeter is implemented.

[0229] The present invention can be a method, a device, a system, and / or a computer program product. The computer program product can include a computer-readable storage medium, on which computer-readable program instructions for executing various aspects of the present invention are loaded.

[0230] Those skilled in the art should understand that the embodiments of the present invention described above and shown in the drawings are only examples and do not limit the present invention. The object of the present invention has been fully and effectively achieved. The function and structural principle of the present invention have been demonstrated and explained in the embodiments, and without departing from the principle, the embodiments of the present invention can have any deformation or modification.

[0231] Finally, it should be noted that: the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that: they can still modify the technical solutions recorded in the foregoing embodiments, or perform equivalent replacements on some or all of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for monitoring abnormal pipeline network flow based on an ultrasonic flowmeter, characterized in that, Including: Obtaining first flow rate data collected by a plurality of first ultrasonic flow meters provided at the inlets of a plurality of pipelines in the pipe network at multiple moments within a preset time period; Obtaining second flow rate data collected by a plurality of second ultrasonic flow meters provided at the outlets of a plurality of pipelines in the pipe network at multiple moments within a preset time period; Determining the association relationship between each pipeline inlet and each pipeline outlet according to the design parameters of the pipe network; Determining the association type between the pipeline inlet and the pipeline outlet with an existing association relationship; Determining the flow response duration between the pipeline inlet and the pipeline outlet according to the first flow rate data of the pipeline inlet with an existing association relationship, the second flow rate data of the pipeline outlet, and the design parameters of the pipe network; Determining whether there is an abnormality in the pipeline between the pipeline inlet and the pipeline outlet with an existing association relationship according to the flow response duration, the design parameters of the pipe network, and the association type; Determining the flow response duration between the pipeline inlet and the pipeline outlet according to the first flow rate data of the pipeline inlet with an existing association relationship, the second flow rate data of the pipeline outlet, and the design parameters of the pipe network, including: Obtaining a first flow rate vector of the first flow rate data of the a-th pipeline inlet with an existing association relationship and a second flow rate vector of the second flow rate data of the b-th pipeline outlet, where a is a positive integer and a is less than or equal to the total number of pipeline inlets, and b is a positive integer and b is less than or equal to the total number of pipeline outlets; Determining the discrete shift response duration according to the first flow rate vector, the second flow rate vector, and the design parameters of the pipe network; Obtaining a first flow rate function of the first flow rate data according to the first flow rate data of the a-th pipeline inlet and the design parameters of the pipe network; Obtaining a second flow rate function of the second flow rate data according to the second flow rate data of the b-th pipeline outlet and the design parameters of the pipe network; Obtaining the continuous shift response duration according to the first flow rate function and the second flow rate function; Obtaining the flow response duration between the pipeline inlet and the pipeline outlet according to the discrete shift response duration and the continuous shift response duration.

2. The abnormal pipeline flow monitoring method based on an ultrasonic flowmeter according to claim 1, wherein Determining the discrete shift response duration according to the first flow rate vector, the second flow rate vector, and the design parameters of the pipe network, including: According to the formula Obtain the maximum discrete similarity S dmax , where is the forward displacement matrix at the s-th step is the negative displacement matrix for the s-th step T if(s＝k，1，0) represents that when the number of steps s is the k-th step, T if(s＝k,1,0) = 1, otherwise T if(s＝k,1,0) = 0, both s and k are positive integers less than or equal to n, and n is the number of moments within a preset time period. is the second flow rate data at the i-th moment of the b-th pipe outlet, S b is the cross-sectional area of the b-th pipe outlet. is the first flow rate data at the i-th moment of the a-th pipe inlet, S a is the cross-sectional area of the a-th pipe inlet, and i is a positive integer less than or equal to n. Determining the number of steps corresponding to the maximum discrete similarity; According to the formula t dA = 2k max Δt Determine the discrete shift response duration t dA , where k max is the number of steps corresponding to the maximum discrete similarity, and Δt is the interval duration between each moment within a preset time period.

3. The method for monitoring abnormal pipeline network flow based on an ultrasonic flowmeter according to claim 2, wherein, Obtaining the continuous shift response duration according to the first flow rate function and the second flow rate function, including: According to the formula Obtain the minimum continuous deviation E cmin , where t0 is the start time of the preset time period, t n is the end time of the preset time period, t is the time within the preset time period, t s is the time offset, f1 is the first flow function, and f2 is the second flow function; Determining the time offset corresponding to the minimum continuous deviation; According to the formula t cA = 2t smin Determine the continuous shift response duration t cA , where t smin is the time offset corresponding to the minimum continuous deviation.

4. The abnormal pipeline network flow monitoring method based on an ultrasonic flowmeter according to claim 3, characterized in that Obtaining the flow response duration between the pipeline inlet and the pipeline outlet according to the discrete shift response duration and the continuous shift response duration, including: According to the formula Obtain the flow response duration t A , where α and β are weight parameters and satisfy the formula 5. The abnormal pipeline flow monitoring method based on an ultrasonic flowmeter according to claim 1, wherein, The association type includes any one of the following: One pipeline inlet is only associated with one pipeline outlet; Two or more pipeline inlets are associated with one pipeline outlet; One pipeline inlet is associated with two or more pipeline outlets; Two or more pipeline inlets are associated with two or more pipeline outlets.

6. The abnormal pipeline flow monitoring method based on an ultrasonic flowmeter according to claim 1, characterized in that Determining whether there is an abnormality in the pipeline between the pipeline inlet and the pipeline outlet with an existing association relationship according to the flow response duration, the design parameters of the pipe network, and the association type, including: Determine the theoretical response duration of the pipeline between the pipeline inlet and the pipeline outlet with an associated relationship according to the design parameters of the pipeline network and the associated type; When the deviation amplitude between the flow response duration and the theoretical response duration is greater than or equal to the amplitude threshold, determine that there is an abnormality in the pipeline between the pipeline inlet and the pipeline outlet with an associated relationship.

7. A pipeline network flow anomaly monitoring system based on an ultrasonic flowmeter, for implementing the method according to any one of claims 1-6, characterized in that, Comprising: A first flow data module, configured to obtain first flow data collected by first ultrasonic flowmeters provided at multiple pipeline inlets in the pipeline network at multiple moments within a preset time period; A second flow data module, configured to obtain second flow data collected by second ultrasonic flowmeters provided at multiple pipeline outlets in the pipeline network at multiple moments within a preset time period; An association relationship module, configured to determine the association relationship between each pipeline inlet and each pipeline outlet according to the design parameters of the pipeline network; An association type module, configured to determine the association type between the pipeline inlet and the pipeline outlet with an associated relationship; A flow response duration module, configured to determine the flow response duration between the pipeline inlet and the pipeline outlet according to the first flow data of the pipeline inlet with an associated relationship, the second flow data of the pipeline outlet, and the design parameters of the pipeline network; An abnormality determination module, configured to determine whether there is an abnormality in the pipeline between the pipeline inlet and the pipeline outlet with an associated relationship according to the flow response duration, the design parameters of the pipeline network, and the associated type.

8. A pipeline network flow anomaly monitoring device based on an ultrasonic flowmeter, characterized in that, Comprising: A processor; A memory for storing instructions executable by the processor; wherein, the processor is configured to call the instructions stored in the memory to execute the method according to any one of claims 1-6.

9. A computer-readable storage medium, characterized in that, A computer program instruction is stored thereon, and when the computer program instruction is executed by the processor, the method according to any one of claims 1-6 is implemented.

Images