CN113435040A - Explosion diameter inversion method based on transient flow - Google Patents

Explosion diameter inversion method based on transient flow Download PDF

Info

Publication number
CN113435040A
CN113435040A CN202110715140.7A CN202110715140A CN113435040A CN 113435040 A CN113435040 A CN 113435040A CN 202110715140 A CN202110715140 A CN 202110715140A CN 113435040 A CN113435040 A CN 113435040A
Authority
CN
China
Prior art keywords
diameter
signal
pipe
path
detonation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110715140.7A
Other languages
Chinese (zh)
Other versions
CN113435040B (en
Inventor
郑国磊
田一梅
彭森
李政翾
程景
梁建文
何楠
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN202110715140.7A priority Critical patent/CN113435040B/en
Publication of CN113435040A publication Critical patent/CN113435040A/en
Application granted granted Critical
Publication of CN113435040B publication Critical patent/CN113435040B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • E—FIXED CONSTRUCTIONS
    • E03—WATER SUPPLY; SEWERAGE
    • E03B—INSTALLATIONS OR METHODS FOR OBTAINING, COLLECTING, OR DISTRIBUTING WATER
    • E03B7/00—Water main or service pipe systems
    • E—FIXED CONSTRUCTIONS
    • E03—WATER SUPPLY; SEWERAGE
    • E03B—INSTALLATIONS OR METHODS FOR OBTAINING, COLLECTING, OR DISTRIBUTING WATER
    • E03B7/00—Water main or service pipe systems
    • E03B7/02—Public or like main pipe systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Hydrology & Water Resources (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Measuring Fluid Pressure (AREA)
  • Pipeline Systems (AREA)

Abstract

The invention relates to a method for realizing the inversion of the diameter of a burst pipe, in order to realize the location of the burst pipe in a pipe network, a high-frequency SCADA system is applied to monitor a pressure abnormal signal in the pipe network, and on the basis of the location of the burst pipe, a total monitoring signal is decomposed into branch monitoring signals from each propagation path according to the propagation characteristics of transient flow; inverting the propagation path and characteristics of the transient flow, namely taking the reciprocal of the local loss and the frictional resistance loss in the original propagation process to obtain gain, inverting the excitation signals from each path, and superposing the excitation signals to form a total inversion excitation signal; and according to the relation between the detonation excitation signal and the detonation diameter, the excitation signal takes the maximum amplitude, the detonation diameter takes the maximum value, and the detonation diameter is calculated. The comparison of the diameter of the detonation calculated by the method with the actual diameter of the detonation and the instantaneous flow attenuation method proves that the positioning precision of the method is higher, and the method can effectively solve the problem of the diameter inversion of the detonation.

Description

Explosion diameter inversion method based on transient flow
Technical Field
The invention relates to a method for inverting the diameter of a burst pipe, belongs to the field of engineering, and particularly relates to a method for inverting the diameter of the burst pipe based on transient flow.
Background
Water supply networks are an important component of the urban infrastructure, called the life line of the city. When the pipe burst happens, a decision is needed to be provided for pipe burst repair according to pipe burst information. In recent years, an SCADA system is widely introduced into a water supply network to monitor flow and pressure in real time, and a plurality of methods for positioning a burst pipe by using monitoring data cannot infer the burst pipe area, so that certain difficulty is brought to burst pipe repair decision. Therefore, the inversion of the caliber of the pipe burst is of great significance.
In recent years, high-frequency pressure monitoring in an SCADA system is introduced into a water supply pipe network, and the transient flow monitored by the system can realize accurate positioning of a complex pipe network, but pipe burst caliber information is difficult to infer. Therefore, it is necessary to research an inversion method of the pipe burst caliber to provide pipe burst repair decision.
Transient flow monitored by the SCADA system is formed by overlapping transient flow excited by pipe explosion and propagated to a monitoring point through multiple paths, and local loss and friction loss exist in the propagation process. On the basis of pipe burst positioning, the propagation path is reversed, two losses are converted into two gains, transient flow is propagated from a monitoring point to a pipe burst point through multiple paths and then overlapped to form an excitation signal, and the pipe burst caliber is calculated according to the corresponding relation between the excitation signal and the pipe burst caliber.
Disclosure of Invention
Aiming at the problem that the caliber of the pipe burst in the water supply network is difficult to invert, the invention provides a pipe burst diameter inversion method based on transient flow on the basis of pipe burst positioning, and the inversion of the caliber of the pipe burst in the water supply network is effectively realized. Therefore, the technical scheme of the invention is as follows:
a burst pipe diameter inversion method based on transient flow comprises the following steps:
s1, decomposing the total monitoring signal into branch monitoring signals from each propagation path according to the propagation characteristics of the transient flow;
s2, reversing the propagation path and the characteristics of the transient flow, and inverting the excitation signal;
and S3, calculating the diameter of the detonator according to the relationship between the detonator excitation signal and the diameter of the detonator.
Further, in step S1, the specific method of the total monitoring signal decomposition is:
according to the propagation theory of the transient flow, the propagating transient flow signals in each propagation path from the pipe explosion point to the monitoring point, called branch analog signals, are calculated and are superposed to form a total analog signal, called a total analog signal, and the amplitude of the monitoring signal from each propagation path can be expressed as:
Figure BDA0003134881970000011
in the formula (1), k is the propagation path number of the transient stream, Ek(t) the amplitude of the branched analogue signal at time t, M (t) the amplitude of the total monitored signal at time t, εk(t) is the ratio of the monitor signal in the kth path at time t.
Further, in step S2, the specific method for inverting the excitation signal is as follows:
and (3) reversely transmitting the signal in each path from the monitoring point to the tube explosion point according to the original path, taking the reciprocal of the local loss and the friction loss in the original transmission process, and changing the reciprocal into gain, wherein the amplitude expression of the excitation signal from each path is as follows:
Figure BDA0003134881970000021
in the formula (2), TkIs the propagation time of the transient stream in the kth path, AkiTransmission or reflection coefficient of i-th node, RkjFriction factor per unit length of j-th section of pipe, e-RkjFrictional drag loss per unit length of j-th section of pipe, LkjIs the length of the j-th section of pipe.
And (3) superposing the excitation signals from each path, wherein the superposed signals are total excitation signals, and the amplitude expression of the superposed signals is as follows:
Figure BDA0003134881970000022
further, in step S3, the specific method of calculating the diameter of the squib is:
regarding the tube explosion surface as a circle, according to the relation between the excitation signal and the tube explosion diameter, the excitation signal takes the maximum amplitude, the tube explosion diameter takes the maximum value, and then the expression is:
Figure BDA0003134881970000023
in the formula (4), ImaxMu is the burst coefficient and D is the burst diameter for the maximum amplitude of the excitation signal.
The diameter of the detonation tube can be calculated by the formula (4).
Drawings
FIG. 1 is a schematic diagram of an experimental pipe network;
FIG. 2 is a graph of a total monitoring signal;
FIG. 3 is a graph of the inverted total excitation signal;
FIG. 4 is a flow chart of the present invention.
Detailed Description
The invention is further illustrated with reference to the following figures and specific examples, which are not intended to limit the invention in any way.
The method comprises the following steps of realizing pipe explosion diameter inversion in an experimental pipe network, and specifically comprises the following operation steps:
s1, decomposition of monitoring signal
The experimental equipment mainly comprises: the water pump, water tank, high frequency pressure sensor (3), DN100, DN50, DN25 steel pipe, DN50PVC pipe, water valve, each pipe section makes up the pipe network that has 2 rings, and high frequency pressure sensor acquisition frequency is 10000 Hz. Each node in the pipe network is numbered with a letter "B" and a number, for example: b1, B5 and B10. Each pipe section number is numbered with a letter "T" and the number of the node at both ends of the pipe section, for example: t2-3, T10-13 (FIG. 1). In an experimental pipe network, a water valve is quickly opened to manufacture a primary pipe explosion event, pressure signals (shown in figure 2) before and after the event are monitored, and a pipe explosion point is positioned on a T9-10 pipe section and 2.3m away from a B9 node by applying a pipe explosion positioning method based on transient flow.
According to the propagation theory of the transient flow, the propagating transient flow signals in each propagation path from the pipe explosion point to the monitoring point, called branch analog signals, are calculated and are superposed to form a total analog signal, called a total analog signal, and the amplitude of the monitoring signal from each propagation path can be expressed as:
Figure BDA0003134881970000031
in the formula (1), k is the propagation path number of the transient stream, Ek(t) the amplitude of the branched analogue signal at time t, M (t) the amplitude of the total monitored signal at time t, εk(t) is the ratio of the monitor signal in the kth path at time t.
S2 inversion of excitation signal
And (3) reversely transmitting the signal in each path from the monitoring point to the tube explosion point according to the original path, taking the reciprocal of the local loss and the friction loss in the original transmission process, and changing the reciprocal into gain, wherein the amplitude expression of the excitation signal from each path is as follows:
Figure BDA0003134881970000032
in the formula (2), TkIs the propagation time of the transient stream in the kth path, AkiTransmission or reflection coefficient of i-th node, RkjFriction factor per unit length of j-th section of pipe, e-RkjFrictional drag loss per unit length of j-th section of pipe, LkjIs the length of the j-th section of pipe.
And (3) superposing the excitation signals from each path, wherein the superposed signals are total excitation signals, and the amplitude expression of the superposed signals is as follows:
Figure BDA0003134881970000033
the monitoring signal from each propagation path obtained in step S1 is propagated in reverse direction, and then the excitation signal from each path is obtained by equation (2), and then the total excitation signal is obtained by equation (3) (fig. 3).
S3, calculating the diameter of the pipe explosion
Regarding the tube explosion surface as a circle, according to the relation between the excitation signal and the tube explosion diameter, the excitation signal takes the maximum amplitude, the tube explosion diameter takes the maximum value, and then the expression is:
Figure BDA0003134881970000034
in the formula (4), ImaxMu is the burst coefficient and D is the burst diameter for the maximum amplitude of the excitation signal.
The maximum amplitude (1.617MPa) in FIG. 3 was taken, and the squib diameter of 32.92mm was calculated by equation (4).
Table 1: comparison table of diameter inversion results of pipe explosion by using transient current attenuation method
Figure BDA0003134881970000041
As can be seen from Table 1, the pipe burst diameter inversion of the invention has an absolute error of 0.92mm and a relative error of 2.88%; compared with a transient flow attenuation method, the absolute value of the inversion accuracy of the diameter of the burst pipe is improved by 4.52mm, and the relative value is improved by 14.12%; in conclusion, the absolute error and the relative error of the reversal of the pipe burst diameter are small, and the problem that the pipe burst diameter in a pipe network is difficult to reverse is solved.

Claims (4)

1. The method for inverting the diameter of the pipe burst based on the transient flow is characterized by comprising the following steps of:
s1, decomposing the total monitoring signal into branch monitoring signals from each propagation path according to the propagation characteristics of the transient flow;
s2, reversing the propagation path and the propagation characteristics of the transient flow, and inverting the excitation signal;
and S3, calculating the diameter of the detonator according to the relationship between the detonator excitation signal and the diameter of the detonator.
2. The transient flow-based method for retrieving the detonation diameter of the pipe in the step S1 is characterized in that the specific method for decomposing the total monitoring signal is as follows:
according to the propagation theory of the transient flow, calculating the propagating transient flow signals in each propagation path from the pipe explosion point to the monitoring point, called branch analog signals, and superposing the branch analog signals to form a total analog signal, called a total analog signal, wherein the amplitude of the monitoring signal from each propagation path can be expressed as:
Figure FDA0003134881960000011
in the formula (1), k is the propagation path number of the transient stream, Ek(t) the amplitude of the branched analogue signal at time t, M (t) the amplitude of the total monitored signal at time t, εk(t) is the ratio of the monitor signal in the kth path at time t.
3. The transient flow-based method for inverting the diameter of the pipe burst according to claim 1, wherein in step S2, the specific method for inverting the excitation signal is as follows:
and (3) reversely transmitting the signal in each path from the monitoring point to the tube explosion point according to the original path, taking the reciprocal of the local loss and the friction loss in the original transmission process, and changing the reciprocal into gain, wherein the amplitude expression of the excitation signal from each path is as follows:
Figure FDA0003134881960000012
in the formula (2), TkIs the propagation time of the transient stream in the kth path, AkiTransmission or reflection coefficient of i-th node, RkjFriction factor per unit length of j-th section of pipe, e-RkjFrictional drag loss per unit length of j-th section of pipe, LkjIs the length of j-th section of pipeAnd superposing the excitation signals from each path, wherein the superposed signals are total excitation signals, and the amplitude expression of the superposed signals is as follows:
Figure FDA0003134881960000013
4. the transient flow-based method for inverting the detonation diameter of claim 1, wherein in step S3, the specific method for calculating the detonation diameter is as follows:
regarding the tube explosion surface as a circle, according to the relation between the excitation signal and the tube explosion diameter, the excitation signal takes the maximum amplitude, the tube explosion diameter takes the maximum value, and then the expression is:
Figure FDA0003134881960000014
in the formula (4), ImaxMu is the burst coefficient and D is the burst diameter for the maximum amplitude of the excitation signal.
CN202110715140.7A 2021-06-26 2021-06-26 Explosion diameter inversion method based on transient flow Active CN113435040B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110715140.7A CN113435040B (en) 2021-06-26 2021-06-26 Explosion diameter inversion method based on transient flow

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110715140.7A CN113435040B (en) 2021-06-26 2021-06-26 Explosion diameter inversion method based on transient flow

Publications (2)

Publication Number Publication Date
CN113435040A true CN113435040A (en) 2021-09-24
CN113435040B CN113435040B (en) 2022-09-23

Family

ID=77754725

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110715140.7A Active CN113435040B (en) 2021-06-26 2021-06-26 Explosion diameter inversion method based on transient flow

Country Status (1)

Country Link
CN (1) CN113435040B (en)

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104359979A (en) * 2014-11-14 2015-02-18 西安交通大学 Method for detecting interlamination transverse cracks of carbon steel/aluminum explosion composite pipe
CN104615576A (en) * 2015-03-02 2015-05-13 中国人民解放军国防科学技术大学 CPU+GPU processor-oriented hybrid granularity consistency maintenance method
CN105787174A (en) * 2016-02-25 2016-07-20 武汉大学 High-rockfill-dam transient-rheological-parameter inversion method based on response surface method
CN107045017A (en) * 2017-04-25 2017-08-15 南京信息工程大学 Crack In Thin Plate depth analysis method based on ultrasonic Lamb waves and time-reversal theory
CN107218516A (en) * 2017-07-19 2017-09-29 中国水利水电科学研究院 A kind of water delivery in pipeline system multiple spot minute leakage detection means and method
CN108332059A (en) * 2018-01-16 2018-07-27 浙江大学 Serve the pressure tap optimization placement method of water supply network booster monitoring
CN108460194A (en) * 2018-02-08 2018-08-28 天津大学 A kind of energisation mode optimization method in pipeline transient flow minute leakage detection
CN110081319A (en) * 2019-05-29 2019-08-02 中国水利水电科学研究院 A kind of design method of controllable transient flow leak detection system
CN110410685A (en) * 2019-07-11 2019-11-05 西北工业大学 A kind of underground pipe network leakage orienting system and method based on time domain reflection technology
CN110762398A (en) * 2019-11-05 2020-02-07 中国石油大学(华东) Oil and gas pipeline leakage detection method and system based on excitation response
CN112393125A (en) * 2020-12-08 2021-02-23 哈尔滨石油学院 Gas pipe network leakage detection system and method based on inverse problem analysis
CN112797860A (en) * 2020-12-30 2021-05-14 西南交通大学 Construction method for preventing and treating rock burst of high-ground-temperature tunnel
CN112989535A (en) * 2021-03-09 2021-06-18 昆明理工大学 Water supply network pressure monitoring point optimal layout method based on pipe burst detection benefit

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104359979A (en) * 2014-11-14 2015-02-18 西安交通大学 Method for detecting interlamination transverse cracks of carbon steel/aluminum explosion composite pipe
CN104615576A (en) * 2015-03-02 2015-05-13 中国人民解放军国防科学技术大学 CPU+GPU processor-oriented hybrid granularity consistency maintenance method
CN105787174A (en) * 2016-02-25 2016-07-20 武汉大学 High-rockfill-dam transient-rheological-parameter inversion method based on response surface method
CN107045017A (en) * 2017-04-25 2017-08-15 南京信息工程大学 Crack In Thin Plate depth analysis method based on ultrasonic Lamb waves and time-reversal theory
CN107218516A (en) * 2017-07-19 2017-09-29 中国水利水电科学研究院 A kind of water delivery in pipeline system multiple spot minute leakage detection means and method
CN108332059A (en) * 2018-01-16 2018-07-27 浙江大学 Serve the pressure tap optimization placement method of water supply network booster monitoring
CN108460194A (en) * 2018-02-08 2018-08-28 天津大学 A kind of energisation mode optimization method in pipeline transient flow minute leakage detection
CN110081319A (en) * 2019-05-29 2019-08-02 中国水利水电科学研究院 A kind of design method of controllable transient flow leak detection system
CN110410685A (en) * 2019-07-11 2019-11-05 西北工业大学 A kind of underground pipe network leakage orienting system and method based on time domain reflection technology
CN110762398A (en) * 2019-11-05 2020-02-07 中国石油大学(华东) Oil and gas pipeline leakage detection method and system based on excitation response
CN112393125A (en) * 2020-12-08 2021-02-23 哈尔滨石油学院 Gas pipe network leakage detection system and method based on inverse problem analysis
CN112797860A (en) * 2020-12-30 2021-05-14 西南交通大学 Construction method for preventing and treating rock burst of high-ground-temperature tunnel
CN112989535A (en) * 2021-03-09 2021-06-18 昆明理工大学 Water supply network pressure monitoring point optimal layout method based on pipe burst detection benefit

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
康宪芝等: "《基于瞬变流法的管道泄漏定位研究》", 《计算机工程与应用》 *
康宪芝等: "《基于瞬变流法的管道泄漏定位研究》", 《计算机工程与应用》, 31 December 2015 (2015-12-31) *
郭瑾: "《基于水力瞬变分析的输油管线漏失检测方法研究》", 《万方数据》 *
郭瑾: "《基于水力瞬变分析的输油管线漏失检测方法研究》", 《万方数据》, 5 May 2016 (2016-05-05) *

Also Published As

Publication number Publication date
CN113435040B (en) 2022-09-23

Similar Documents

Publication Publication Date Title
CN109114424B (en) Method and device for determining position of pipeline leakage point
CN104180166A (en) Pipeline leakage detection method based on pipeline pressure data
CN113435040B (en) Explosion diameter inversion method based on transient flow
CN103194758B (en) Forming device for gas production pipeline corrosion inhibitor slug pre-film
CN108804537B (en) Method for internally detecting data alignment
CN111310323B (en) Deep sea pipeline comprehensive risk index evaluation method based on MATLAB software
CN113221354A (en) Pipeline bending deformation fitting algorithm
CN113446521B (en) Tube burst positioning method based on transient flow
CN110688782B (en) Long-distance superconducting energy pipeline simulation method considering axial and radial field distribution
CN112668182A (en) Analysis method for failure of natural gas long-distance pipeline
CN103245222A (en) Pipeline configuration structure of triple-inlet air cooler system for hydrogenation reaction effluents
CN208666774U (en) It is a kind of for transporting the fixture of heavy-caliber plastic pipe and composite coiled tubing
CN111578147A (en) Automatic fault diagnosis method for coal mine underground coal seam gas extraction pipeline
CN109307158B (en) Method and device for determining pipeline leakage
Tabella et al. Data fusion for subsea oil spill detection through wireless sensor networks
CN105013823A (en) Processing technology of production method of seamless steel tubes for high-temperature and high-pressure pipes
CN114971076A (en) Multi-objective optimal arrangement method for monitoring points of water supply network
CN205877444U (en) Multi -functional non - metal composite pipe
CN101004231A (en) Polyethylene wearproof composite tube with steel tightened liner, and preparation method
Liu et al. Research on detecting and locating tubing leakage of offshore gas wells based on acoustic method
CN209413865U (en) It is a kind of can fixed data acquisition optical fiber box cupling
CN208487297U (en) A kind of heavy caliber high flexibility intelligent pipeline
CN105160062A (en) Hydraulic check method for reverse return pipe network
Zhang et al. Application of the cable laying coiled tubing in electric submersible pump
CN203950586U (en) A kind of fluoroplastic polymer insulation high-temperature-resistant cable

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant