CN107329167B - Detect method, system and the computer storage media of pipeline sediment distribution - Google Patents
Detect method, system and the computer storage media of pipeline sediment distribution Download PDFInfo
- Publication number
- CN107329167B CN107329167B CN201710431352.6A CN201710431352A CN107329167B CN 107329167 B CN107329167 B CN 107329167B CN 201710431352 A CN201710431352 A CN 201710431352A CN 107329167 B CN107329167 B CN 107329167B
- Authority
- CN
- China
- Prior art keywords
- wave
- test point
- deposit
- intensity
- pressurization
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/02—Generating seismic energy
- G01V1/133—Generating seismic energy using fluidic driving means, e.g. highly pressurised fluids; using implosion
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Pipeline Systems (AREA)
Abstract
The invention discloses a kind of method, system and computer storage medias for detecting pipeline sediment distribution.Method includes the following steps: obtaining the pressure value of the fluid at the test point being arranged in pipeline within a preset period of time since when incidence wave is inspired in pipeline；The mapping relations at time point and pressure value are generated according to the pressure value of acquisition；The arrival time of incidence wave and the arrival time of intensity, the arrival time for the pressurization wave being reflected at test point and intensity and the decompression wave being reflected at test point are obtained according to mapping relations；Length of the deposit apart from test point is determined according to the arrival time of incidence wave and pressurization wave；The length of deposit is determined according to the arrival time of pressurization wave and decompression wave；Length according to the intensity of incidence wave, the intensity and deposit of pressurization wave apart from test point determines gross blockage factor range；The length and gross blockage factor of length, deposit according to deposit apart from test point determine the distribution situation of deposit.
Description
Technical field
The present invention relates to pipe detection field more particularly to it is a kind of detect the method for pipeline sediment distribution, system and
Computer storage media.
Background technique
When the complex fluid of a variety of mixtures composition flows in pipeline, fluid will receive the conditions such as temperature, pressure
It influences, and then is easy to generate deposit on the inner wall of pipeline.Deposit may result in tube wall incrustation, in some instances it may even be possible to will appear
The accidents such as deposit blocking, solidifying pipe in managing, to can pipeline transportation be adversely affected and be endangered.
In order to avoid the generation of the accidents such as deposit blocking, solidifying pipe, need periodically to be purged the deposit in pipeline,
And the position of deposit and its deposit thickness in pipeline are needed to be grasped before removing to deposit.
In the prior art, exist and object detecting method is deposited to pipeline using pressure wave.This method is a certain usually using pipeline
The fluid pressure data of position or certain positions can determine the distribution situation of deposit along pipeline.Current pressure wave technology
The case where being only capable of the blocking of qualitative recognition local pipe, i.e., can only determine the blocking position and stemming length of tamper, but can not
The specific distribution situation of deposit in quantitative accurate detection pipeline.
Summary of the invention
In view of the above-mentioned problems, the purpose of the present invention is to provide it is a kind of detect the method for pipeline sediment distribution, system with
And computer storage media, for solving above-mentioned technical problem.
To achieve the goals above, a kind of method detecting pipeline sediment distribution provided by the invention, including following step
It is rapid:
Since when incidence wave is inspired in specified duct section, obtained in the duct section within a preset period of time
Test point at fluid pressure data；
Determine the mapping relations between the pressure value and corresponding sampling time in the fluid pressure data；
According to the mapping relations, obtain arrival time and intensity of the incidence wave at the test point, it is described enter
Arrival time and intensity and the incidence of the pressurization wave that ejected wave is formed in the duct section internal reflection at the test point
Arrival time of the decompression wave that wave is formed in the duct section internal reflection at the test point；
According to the arrival time of the incidence wave and the arrival time of the pressurization wave, the deposition in the duct section is determined
Length of the object apart from the test point；
According to the arrival time of the pressurization wave and the arrival time of the decompression wave, the length of the deposit is determined；
According to intensity at the test point of intensity of the incidence wave at the test point, the pressurization wave and
Length of the deposit apart from the test point determines the gross blockage factor range of the duct section；
The length of length, the deposit according to the deposit apart from the test point and the gross blockage factor model
It encloses, determines deposit distribution situation in the duct section.
In an optional embodiment, according to the arrival time of the incidence wave and it is described pressurization wave arrival time,
Determine length of the deposit apart from the test point in the duct section, comprising:
Obtain pipeline elasticity modulus, pipe diameter, pipe thickness, fluid density and the elasticity of fluid at the test point
Coefficient；
According to formulaDetermine the average velocity of wave in the duct section；
Wherein, aijFor the average velocity of wave of i-th section of j-th of spatial mesh size；fijFor the frictional resistance system of i-th section of j-th of spatial mesh size
Number；KijFor the elasticity of fluid coefficient of i-th section of j-th of spatial mesh size；ΔxijFor i-th section of j-th of spatial mesh size；ρ is that fluid is close
Degree；E is pipeline elasticity modulus；DiFor the pipe diameter of the i-th segment pipe；AiFor the cross-sectional area of the i-th segment pipe；E is thickness of pipe wall
Degree；C is constraint factor related with duct section constraint condition；
According to the arrival time and the average velocity of wave of the incidence wave and the pressurization wave determine the deposit away from
Length from the test point.
It is described according to the determination of the arrival time of the pressurization wave and the decompression wave in an optional embodiment
The length of deposit, comprising:
The deposit is determined according to the arrival time and the average velocity of wave of the pressurization wave and the decompression wave
Length.
In an optional embodiment, the intensity according to the incidence wave at the test point, the pressurization
Intensity and deposit length apart from the test point of the wave at the test point, determines the blocking of the duct section
Coefficient range, comprising:
Length according to the deposit apart from the test point, and determine that the incidence wave reaches institute by numerical simulation
State intensity when deposit proximal end；
The deposit proximal end is reached according to intensity of the incidence wave at the test point, and in conjunction with the incidence wave
When intensity determine the attenuation of the incidence wave；
Declining for the pressurization wave is determined by friction drag loss at position according to the attenuation of the incidence wave, pressure wave
Decrement；
According to the attenuation of the incidence wave, the attenuation of the pressurization wave, the incidence wave at the test point
The intensity of intensity and the pressurization wave at the test point, determines the gross blockage factor of the estimation of the duct section, and then really
Determine gross blockage factor range.
In an optional embodiment, the length according to the deposit apart from the test point, and pass through number
Value simulation determines the intensity when incidence wave reaches the deposit proximal end, comprising:
Pass through formula: Progress numerical simulation determines strong when the incidence wave reaches the deposit proximal end
Degree；
Wherein, x is the distance along duct section centerline direction；H is head；V is the mean flow rate of fluid in duct section；g
For acceleration of gravity；JsFor friction loss；F is the coefficient of friction resistance in Brunone model；kBIt is constant；A is water-hammer wave speed；α is pipe
Section inclination angle.
In an optional embodiment, the attenuation according to the incidence wave, pressure wave friction loss at
Determine the attenuation of the pressurization wave, comprising:
Pass through formula:Determine the attenuation of the pressurization wave；
Wherein, Δ hFAttenuation when reaching the deposit from the test point for the pressurization wave；ΔhfFor it is described enter
Ejected wave reaches attenuation when test point from the deposit；ΔHFDescribed sink is reached from the test point for the pressurization wave
Friction loss when at product object；ΔHfFriction loss when test point is reached from the deposit for the incidence wave.
It is described according to the attenuation of the incidence wave, the attenuation of the pressurization wave, institute in an optional embodiment
Intensity and pressurization wave intensity the test point at of the incidence wave at the test point is stated, determines the duct section
Gross blockage factor, comprising:
Pass through formula: Determine blocking
Coefficient range；
Wherein, k is local stoppages coefficient；D for no local stoppages pipeline section caliber；C is related with pipe-supporting mode
Parameter；HFFor the intensity of incidence wave；F1For incidence wave HFIntensity when reaching at deposit；f1For at deposit to test point
The intensity of the pressurization wave of reflection；HfFor intensity of the incidence wave at the test point；ΔhFIt is the pressurization wave from the inspection
Attenuation when measuring point reaches at the deposit；ΔhfDeclining when reaching test point from the deposit for the incidence wave
Decrement.
In an optional embodiment, the incidence wave is inspired by turning on or off pumping plant, or is passed through
The valve body positioned at pipe downstream end is closed to be inspired；
When the incidence wave is inspired by turning on or off pumping plant, the test point is located at the pumping plant；
When the incidence wave is inspired by the valve body that closing is located at pipe downstream end, the test point is located at described
At valve body.
In an optional embodiment, the step of determining the preset time period, comprising:
Determine the maximum allowable working pressure of duct section to be detected；
According to the maximum allowable working pressure of duct section to be detected, maximum allowable detection range is determined；
Preset time period is determined according to maximum allowable detection range.
In an optional embodiment, the relationship change curve of the mapping relations m- pressure when including.
Present invention also provides a kind of systems for detecting pipeline sediment distribution, comprising:
Pressure detecting module, for since when incidence wave is inspired in specified duct section, in preset time period
Fluid pressure data at the interior test point obtained in the duct section；
Graph generation module, for determining the mapping between the pressure value in the fluid pressure data and corresponding sampling time
Relationship；
Parameter acquisition module obtains arrival time of the incidence wave at the test point according to the mapping relations
And intensity, arrival time of the pressurization wave that the incidence wave is formed in the duct section internal reflection at the test point and strong
Arrival time of the decompression wave that degree and the incidence wave are formed in the duct section internal reflection at the test point；
First processing module, for determining according to the arrival time of the incidence wave and the arrival time of the pressurization wave
Length of the deposit apart from the test point in the duct section；
Second processing module, for determining the deposit according to the arrival time of the pressurization wave and the decompression wave
Length；
Third processing module, according to intensity of the incidence wave at the test point, the pressurization wave in the detection
The length of intensity and the deposit apart from the test point at point, determines the gross blockage factor range of the duct section；
Fourth processing module, the length of length, the deposit according to the deposit apart from the test point and
The gross blockage factor range determines deposit distribution situation in the duct section.
Present invention also provides a kind of computer storage medias, store computer program, the computer program thereon
It is performed the steps of when being executed by processor
Since when incidence wave is inspired in specified duct section, obtained in the duct section within a preset period of time
Test point at fluid pressure data；
Determine the mapping relations between the pressure value and corresponding sampling time in the fluid pressure data；
According to the mapping relations, obtain arrival time and intensity of the incidence wave at the test point, it is described enter
Arrival time and intensity and the incidence of the pressurization wave that ejected wave is formed in the duct section internal reflection at the test point
Arrival time of the decompression wave that wave is formed in the duct section internal reflection at the test point；
According to the arrival time of the incidence wave and the arrival time of the pressurization wave, the deposition in the duct section is determined
Length of the object apart from the test point；According to the arrival time of the pressurization wave and the arrival time of the decompression wave, institute is determined
State the length of deposit；
According to intensity at the test point of intensity of the incidence wave at the test point, the pressurization wave and
Length of the deposit apart from the test point determines the gross blockage factor range of the duct section；
The length of length, the deposit according to the deposit apart from the test point and the gross blockage factor model
It encloses, further more accurately determines deposit distribution situation in the duct section.
Method provided by the invention compared with the existing technology in pressure wave detection method, the position in addition to can determine that deposit
It sets outer, moreover it is possible to determine the length of deposit and the degree of blocking, have the characteristics of easy to operate, good economy performance.
Referring to following description and accompanying drawings, only certain exemplary embodiments of this invention is disclosed in detail, specifies original of the invention
Reason can be in a manner of adopted.It should be understood that embodiments of the present invention are not so limited in range.In appended power
In the range of the spirit and terms that benefit requires, embodiments of the present invention include many changes, modifications and are equal.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments of invention for those skilled in the art without any creative labor, can be with root
Other attached drawings are obtained according to these attached drawings.
Fig. 1 is the flow chart of the method for the detection pipeline sediment distribution that the embodiment of the present invention provides；
Fig. 2 is that the pipeline provided by the present invention based on pressure wave deposits in object detecting method, local stoppages section pressure wave
Reflect schematic diagram；
Fig. 3 is the schematic diagram that the pipeline provided by the present invention based on pressure wave deposits safe distance in object detecting method；
Fig. 4 is the time and pressure in the pipeline deposition object detecting method provided by the present invention based on pressure wave at test point
Relation curve between force value；
Fig. 5 is that pipeline provided by the present invention based on pressure wave deposits in object detecting method each pressure wave in the production of pipeline
Raw, transmitting and reflection schematic diagram；
Fig. 6 is the non-unified spatial mesh size grid dividing signal of the pipeline deposit provided by the present invention based on pressure wave
Figure；
Fig. 7 is that the pipeline provided by the present invention based on pressure wave deposits in object detecting method, by the deposit in immersed tube
It is equivalent to reducer pipe and its method of characteristic curves schematic diagram；
Fig. 8 is the ginseng matrix that the pipeline provided by the present invention based on pressure wave deposits local stoppages section in object detecting method
Number schematic diagram；
Fig. 9 is that the pipeline provided by the present invention based on pressure wave deposits pump-pipe-line system signal in object detecting method
Figure；
Figure 10 is that the pipeline provided by the present invention based on pressure wave deposits pipe end pressure in object detecting method embodiment 1
Power variation diagram；
Figure 11 is that the pipeline provided by the present invention based on pressure wave deposits moment closing end in object detecting method embodiment 2
Change of water level figure shown in the setting survey of end valve door；
Figure 12 be the pipeline provided by the present invention based on pressure wave deposit simulated in object detecting method embodiment 2 it is determining
To local stoppages section caliber and actual blockage caliber comparing result figure；
In the method for the detection pipeline sediment distribution that Figure 13 provides for the embodiment of the present invention in step S106 step by step
Flow chart；
Figure 14 is the structural schematic diagram of the system of detection pipeline sediment distribution provided by the present invention.
Specific embodiment
In order to make those skilled in the art better understand the technical solutions in the application, below in conjunction with the application reality
The attached drawing in example is applied, the technical scheme in the embodiment of the application is clearly and completely described, it is clear that described implementation
Example is merely a part but not all of the embodiments of the present application.Based on the embodiment in the application, this field is common
The application protection all should belong in technical staff's every other embodiment obtained without creative efforts
Range.
Shown in referring to Fig.1, which is a kind of flow chart of the method for detection pipeline sediment distribution of the invention.This method
It may comprise steps of:
S101: since when incidence wave is inspired in specified duct section, the pipe is obtained within a preset period of time
The fluid pressure data at test point in road section；
In this step, when the fluid in duct section is in a stable state, fluid pressure data in duct section,
It includes pressure value, will not generally occur more significantly to change, and is fluctuated in a metastable range.When in duct section
Somewhere when encountering blocking or closing, the fluid of trip disposed thereon because the effect continuation of inertia force is flowed downstream,
And then (extruding) is compressed to the fluid in downstream, so that the pressure moment for being located at downstream fluid increases.
After increasing the pressure moment of fluid, according to elastically-deformable principle in the prior art, fluid can be inspired pressure
Reeb.Pressure value variation is bigger, and the intensity of pressure wave is bigger, therefore can be by the way that pressure sensor, pressure is arranged in test point
Force snesor obtains the size of the pressure value at the point, can also detect that this presses the intensity of Reeb.
Referring to shown in Fig. 2, the different location in pipeline is can be set in test point.For example, pipeline transporter generally can be with
Including the pumping plant of upstream position is arranged in, the valve body of downstream position is arranged in and is connected to the pipeline of upstream and downstream.Pumping plant can be
Pressure needed for flowing provides in fluid pipeline, valve body also can control the flow velocity of fluids within pipes.
Pressure wave (incidence wave) can be generated by turning on or off pump, or be located at pipe downstream end by closing
Valve body generate.It is closing or is opening the moment of pumping plant and pressure wave can be inspired within the moment of valve body, pipeline,
Therefore test point can be at the position namely pumping plant import and export or valve body for inspiring pressure wave.
Wherein, it is also not routinely when pressure sensor continuously acquires the pressure value of fluid in the duct section at test point
It is obtained, but is obtained once with the interval preset time, carry out acquisition pressure for example, can be spaced 0.1 second, 1 second, 5 seconds etc.
Force value.Set time interval should reflect that the variation of pressure wave；And pressure data is decomposed with wavelet transformation, reconstruct, mistake
Filter high-frequency pressure noise.
In this step, the preset time that pressure sensor obtains pressure value can be born most according to pipe under test section
Big pressure value is determined.Specific method of determination are as follows: it is assumed that pipe under test maximum allowable working pressure H, in which:
H=H1+Hb1
H1For the pressure wave of excitation；Hb1For the to be measured section of pressure apart from pressure wave excitation point farthest；
Then, the distance between the position that maximum pressure value is generated on test point to duct section is determined according to maximum pressure value
(safe distance).Referring to shown in Fig. 3, by the point X for reaching maximum pressure value on the available duct section of determination2, the point and inspection
The distance between measuring point is then safe distance, then determines the velocity of wave in the pipe under test section, and then can determine that pressure wave reaches
Time at this, it may thereby determine that preset time.
Specifically, superpressure along pipeline in order to prevent, can take safety measures in back wave upstream (as opened pressure release
Valve).It should be guaranteed that the transmitting of taken security measures not parastic echo, that is, the position to be taken measures will be anti-
Ejected wave upstream.With the transmitting of back wave, upstream region (trapezoid area in Fig. 3) also changes.On the right side of trapezoid area
End valve is opened after corresponding time point, terminates measurement.
S102: the mapping relations between the pressure value and corresponding sampling time in the fluid pressure data are determined；
In this step, mapping relations can relationship change curve between pressure value and time.Specifically, pressure passes
Sensor is by the pressure value data transmission in the preset period to video generation device.The video generation device can be according to the number
Accordingly and the data corresponding time generates pressure value and the variation relation curve between the time.Referring to shown in Fig. 4, which is exactly
The pressure value and the relation curve of time generated according to the pressure value of acquisition.Wherein, lateral coordinates are sampling time, longitudinal coordinate
For pressure value.It can be seen that in t1When, pressure value has one to be mutated upward, shows incidence wave generation at this time, then pressure value
As the time constantly increases.During routinely increasing, once upward mutation (test point can be also reflected into
It is pressurized wave) and primary downward mutation (being reflected into the decompression wave of test point).
Of course, mapping relations in addition to being indicated in a manner of relationship change curve, can also pass through table, matrix etc.
Mode is indicated.Computer can determine the mapping relations according to preset condition, and then determine each parameter value.Example
Such as, the range that multiple mutation can be set in computer determines the pressure value when the variation of pressure value is greater than the mutational range
For the value of incidence wave or pressurization wave or decompression wave.
S103: according to the mapping relations, arrival time and intensity of the incidence wave at the test point, institute are obtained
State arrival time and intensity at the test point of pressurization wave that incidence wave is formed in the duct section internal reflection and described
Arrival time of the decompression wave that incidence wave is formed in the duct section internal reflection at the test point；
In this step, with continued reference to Fig. 4 and as shown in connection with fig. 5, incidence wave can be generated by closing valve, work as incidence
While wave generates, test point is located at the position of valve, and the pressure sensor at test point detects the prominent of a pressure value in turn
Become.Referring to shown in Fig. 4, t1This is the time generated for the incidence wave to point, and can obtain incidence wave by the corresponding coordinate of the point
Intensity when generation.After incidence wave generates, incidence wave can upstream direction propagation.Wherein, F indicates that incidence wave, f indicate reflection
Wave, F1To F5It is propagated for incidence wave upstream direction and further through the reflected schematic diagram of deposit.
Simultaneously, referring to Figure 5, after the valve of pipe downstream is closed suddenly, the fluid of upstream can be in pumping plant
Continue to move downward under effect, build the pressure in pipeline section (C after resulting in blockage a little1To between test point), raw transient state positive pressure wave f1(increase
Press wave)；And pipeline section (the C after plugging point2To between pumping plant) due to losing segment fluid flow supply, and transient pressure drop is generated, phase
When in generation transient negative pressure f2(decompression wave).Decompression wave f2By C1Reflection can occur when point again and transmit, wherein at valve after
Resuming the decompression wave passed is F5。
After being pressurized wave generation, can downstream it be propagated along the direction of pipeline.Therefore, figure 4, it is seen that in t2When
It carves, pressurization wave reaches at test point, and the pressure at test point can generate mutation because being pressurized the appearance of wave.It therefore, can be from change
Change the intensity of arrival time and the wave that pressurization wave is obtained on relation curve.Similarly, after pressurization wave reaches test point, decompression wave
It can be in t3Moment reaches test point, the arrival time of available decompression wave from Fig. 4.
S104: it according to the arrival time of the incidence wave and the arrival time of the pressurization wave, determines in the duct section
Length of the deposit apart from the test point；
In this step, according to the propagation principle of wave, the generation time of incidence wave is identical with the pressurization generation time of wave, produces
After life, incidence wave upstream direction is propagated, and pressurization wave direction downstream direction is propagated.It is assumed that the spread speed of wave in a fluid is identical,
When being pressurized wave arrival test point, can wait is all incidence wave in C1Place reflects and reaches test point (at valve), therefore
It can be according to the C for determining deposit in the way of sound ranging in the prior art1Distance x of the point apart from test point.Wherein, public
Formula are as follows:
A is average velocity of wave, unit m/s；t1For the arrival time of incidence wave, unit s；t2For reflected pressurization
Wave arrival time, unit s.
Wherein, in this step, average velocity of wave a can be provided with the springform of the duct section of the test point by obtaining
Amount, pipe diameter, pipe thickness, in pipeline in the density of fluid and duct section fluid coefficient of elasticity, and according to formula:
The average velocity of wave of each spatial mesh size is calculated, then ask the average value of all step-length velocities of wave obtain along average wave
Fast a.
In formula (2), aijFor the average velocity of wave of i-th section of j-th of spatial mesh size, unit m/s；fijFor i-th section j-th
The coefficient of friction resistance of spatial mesh size；KijFor the elasticity of fluid coefficient of i-th section of j-th of spatial mesh size, unit Pa；ΔxijIt is i-th section
J-th of spatial mesh size, unit m；ρ is fluid density, unit kg/m3；E is tubing elasticity modulus, unit Pa；DiIt is i-th
Segment pipe diameter, unit m；AiFor the i-th segment pipe cross-sectional area, unit m2；E is pipe thickness, unit m；C be and pipe
The related constraint factor of road constraint condition.
S105: according to the arrival time of the pressurization wave and the arrival time of the decompression wave, the deposit is determined
Length；
In this step, reference Fig. 4 is shown and above content is it is found that pressurization wave is C1Point is issued to test point.Decompression
Wave is C2Point is issued to test point.The arrival time interval of pressurization wave and decompression wave is exactly the time (C that wave is propagated in deposit1
To C2Point), therefore can be according to the length L for determining deposit using the formula of sound ranging in the prior art；Formula are as follows:
Wherein, a is average velocity of wave, unit m/s；t2The time of test point, unit are reached for reflected pressurization wave
For s；t3The time of test point, unit s are reached for reflected decompression wave.
S106: strong at the test point according to intensity of the incidence wave at the test point, the pressurization wave
The length of degree and the deposit apart from the test point, determines the gross blockage factor range of the duct section；
Referring to Fig.1 shown in 3, step S106 further include it is following step by step:
S1061: the length according to the deposit apart from the test point determines that incidence wave reaches institute by numerical simulation
State intensity when deposit proximal end；
In this step, because that can not determine the location point of deposit, therefore can not before determining deposit specific location
Pressure sensor can be provided at deposit goes detection incidence wave to reach deposit proximal end (C1) intensity, we so need
It is calculated by simulation come intensity when determining that incidence wave reaches this.Wherein, referring to Figure 5, when incidence wave is by closing valve
When door generates, deposit proximal end is in the opposite both ends of deposit close to one end of valve.
Specifically, in the prior art, governing equation of the fluid in the duct under transient condition includes:
In formula (4)-(7), x is the distance along pipe centerline direction, unit m；T is time, unit s；H is water
Head, unit m；V is the mean flow rate of fluids within pipes, unit m/s；G is acceleration of gravity, unit m2/s；JsFor frictional resistance
Loss, unit m；F is the coefficient of friction resistance in Brunone model；kBIt is the empirical between 0.01 and 0.03；A is water
Hit velocity of wave, unit m/s；α is pipeline section inclination angle, and unit is degree.
Formula (4)-(7) constitute a partial differential equation group, which can be in the feelings for the position for determining deposit
Under condition, the situation of change of pressure wave (including incidence wave) is determined, therefore can length according to deposit apart from the test point
Intensity when determining that incidence wave reaches the deposit proximal end by numerical simulation.
In addition, the accuracy in order to guarantee numerical simulation.In actual calculating process, the numerical computation method considers
Pipe roughness caused by viscosity variation, deposit changes along institute trandfer fluid, describes series connection internal edges with can be convenient
Boundary's condition (shown in Fig. 6 to Fig. 8).In figure, in order to distinguish the duct section of upstream and downstream, using the representation method of double subscript, wherein
One subscript indicates that the number of pipeline section, second subscript indicate that the spatial mesh size of pipeline section section number or the cross-sectional right side is compiled
Number, subscript P indicates the variable at t+ Δ t time point.Wherein, spatial mesh size refers to the segment space in duct section along its length.
Using the method for characteristic curves, the 1st node of N+1 node and downstream pipeline section for upstream pipeline section has:
C+:HPi,N+1=CPi,N-Bi,NQPi,N+1 (8)
C-:HPi+1,1=CMi+1,1+Bi+1,1QPi+1,1 (9)
CPi,N=Hi,N+Bi,NQi,N-Ri,NQi,N|Qi,N| (10)
CMi+1,1=Hi+1,2-Bi+1,1Qi+1,2+Ri+1,1Qi+1,2|Qi+1,2| (11)
Wherein, formula (8)-(13) are for discrete partial differential equation group (4)-(7), to determine that incidence wave reaches the deposit
Intensity F when proximal end1。
S1062: intensity when reaching the deposit proximal end according to the intensity of the incidence wave and the incidence wave is true
The attenuation of the fixed incidence wave；
In this step, referring to the content of above-mentioned S1061, formula can be passed through:
ΔhF=HF-F1 (14)
Determine incidence wave attenuation Δ hF, wherein HFIndicate intensity of the incidence wave at test point, F1For incidence wave arrival
Deposit proximal end (the C1) when intensity.
S1063: according to the attenuation Δ h of the incidence waveF, incident intensity of wave HFThe value H of end is reached with back wavef
Determine gross blockage factor range；
In this step, the case where closing valve for pipe end can obtain gross blockage factor model by solving following equation
It encloses:
Wherein, k is local stoppages coefficient, dimensionless；Caliber of the D for no local stoppages pipeline section, unit m；C is and pipeline
The related parameter of support pattern, dimensionless；HFFor the intensity of incidence wave, unit Pa；HfFor pressurization wave reach test point when it is strong
Degree, unit Pa；ΔhFFor incidence wave HFDecaying when being reached from deposit from test point, unit Pa；ΔhfFor pressurization wave from
Decaying when test point, unit Pa are reached at deposit.
In formula (15)-(17), there are two unknown quantity Δ hfAnd k, can not direct solution go out k.Below by Δ hfIt is indicated with k
Out.Assuming that the decaying of pressure wave and pressure wave, available Δ h proportional by the friction loss at positionfBetween k
Functional relation:
Wherein, Δ HfWith Δ HFRespectively pressure wave by the friction loss at position.
Assuming that the position and flow distribution after incidence wave, back wave wave crest warp is linear after end moment closes valve.It can obtain
Flow distribution (Fig. 9) after to incidence wave wave crest:
Wherein, Q0Flow when for pipeline steady-state operation.
Pressure after back wave wave crest is distributed (Fig. 9):
Wherein, r12It (k) is incidence wave in C1The reflection coefficient at place is the function of k.It is calculated with following formula:
The distance between end to blocking section proximal end is divided into n sections, every segment length is Δ x=l/n, is calculated separately each
Then friction loss on Δ x is accumulated by Δ HfWith Δ HF.It is expressed with following formula
In formula, Qil、QilThe respectively flow of pressure wave left and right side.Have for incidence wave: Qil=Q0, Qir=Q1(i
Δx)；Have for back wave: Qil=Q2(i Δ x), Qir=Q1(iΔx)。
Due to Δ HfDependent on gross blockage factor k, equation (15) intercouples with (18), can not direct solution, use is following
Searching algorithm solves it in (0,1) section, and solution procedure is as follows:
(1) a gross blockage factor initial value k=k is set0；
(2) Δ h is obtained with formula (18)-formula (22) solutionf；
(3) the Δ h for solving previous stepfIt substitutes into equation (15) (16), solution obtains k=k1；
(4) if | k0-k1| < δ (precision that δ is setting), exits calculating；Otherwise k=k+ Δ k is updated, (2) is returned and carries out
It calculates, it is final to determine a gross blockage factor range.
Shown in referring to Fig.1, embodiment provided by the invention is further comprised the steps of:
S107: the length and the blocking of length, the deposit according to the deposit apart from the test point
Coefficient range determines deposit distribution situation in the duct section；
In this step: calculated result and reality in order to minimize the method for detecting pipeline sediment distribution of the invention
Error between actual value, position of one or more deposit apart from test point on determining pipeline, the length of deposit and
On the basis of gross blockage factor range, pass through formula
The fine solution of further progress, the accurate distribution situation for determining deposit.
In formula (23), OF is objective function；DT is total time step；H is the piezometric head of simulation, m；H* is measurement
Piezometric head out, m；nHIt is the quantity of measurement point；BLOC is the matrix of variables to be recognized.
Wherein, there are in the duct section of n deposit section, sediment distribution situation uses the matrix B LOC of n × 3
To indicate (as shown in Figure 8):
Every a line represents a deposit section (duct section) in formula (24), and i-th section of deposit section is retouched using 3 parameters
It states: distance x of the deposit apart from test pointi, deposit length li, gross blockage factor ki。
Specifically, the result that formula (1), (3) and (15) determines all is the estimated value of three blocking parameters.In order to more quasi-
It really determines that the distribution situation of deposit, applying equation (23) are solved, but needs to provide it with asking for a blocking parameter
Solve range.Therefore, the value range of blocking position and stemming length is estimated according to formula (1) and (3), the value of gross blockage factor range
Range is estimated according to formula (15).After the value range for determining three blocking parameters, the blocking that multiple groups determine can be joined
Number is updated in matrix (24), is then calculated again by formula (23).
Although procedures described above process includes the multiple operations reached with particular order, it should however be appreciated that understand,
These processes may include more or fewer operations, these operations can be executed sequentially or be executed parallel (such as using parallel
Processor or multi-thread environment)
Embodiment 1
In order to verify the feasibility for the method that embodiments herein provides, an experimental system is set, which shows
It is intended to as shown in Figure 9.The station spacing of the horizontal pipe is 100km, internal diameter of the pipeline 443mm, and the absolute equivalent roughness of pipeline is
0.01mm, elasticity modulus are 207 × 103MPa.The density of Experimental Flowing Object is 842.1kg/m3, initial station pump lift is 30m.Pipeline
On be provided with regulating valve.Initial station intermediate station has control valve.Downstream boundary is valve and fixed water level tank, in fixed water level tank
Pipe end valve generates a constant pressure 0.4MPa.Assuming that play its maximum under the premise of meeting conditions above defeated for the system
Send ability.
Assuming that having a length at the second distance between sites intermediate station 45km is the blocking section of 2km, gross blockage factor 0.7.
End valve moment closes valve, and simulation obtains end and blocking endpoint (C in Fig. 51Point) at pressure change.It is passed now by pressure
Sensor detects the variation of the pressure value at valve and the variation production Methods curve (solid line in Figure 10) according to pressure value, it is desirable that office
Portion blocks parameter.
The method according to the invention, according to the available H of the relation curve of pressure value and timeF、Hf、t1、t2、t3.Velocity of wave a
It can be calculated by formula (2).First choice calculates position and the length of local stoppages according to formula (1) and (3).It is stifled estimating
After filling in position, Δ hFIt can be estimated by numerical simulation, then pass through formula: Δ hf=(3-4ln2) Δ hFBy Δ hfIt determines
Out.
The above given data is substituted into formula (16) available local stoppages coefficient.
Estimation method inverse obtains local stoppages parameter x=30891.36m, l=2007.52m, k=0.65, it is opposite with
The numerical error that actual experiment determines is respectively 0.3%, 0.3%, 7.1%.
Embodiment 2
It is identical with pump-pipe-line system in example 1, it is assumed that there are a local stoppages section between second station, parameter is such as
Shown in table 3, it is assumed that moment closes end valve, is arranged two measurement points (end and at the 10km of end), it is available such as
Change of water level figure shown in Figure 11.It is now assumed that having obtained change of water level figure as shown in figure 11 by measurement, it is desirable that 7 blockings
Gross blockage factor (the respectively k of section1、k2、k3、k4、k5、k6、k7)。
Firstly, being programmed with the fluid line transient flow calculation procedure of local stoppages according to method described in patent specification；
Then, the optimal objective function value function of discrete form is write out；And by pressure waveform in observation Figure 11, following constraint item is added
Part: k1﹤ k2﹤ k3﹤ k4, and k4﹥ k5﹥ k6﹥ k7, it is set as [0.6,1] by estimation by bound is solved, then the blocking parameter is sharp
It is solved with formula (23), and is solved twice with genetic algorithm, solving maximum algebra was 100 generations, and average used time 14.58h is obtained
The caliber (table 1, Figure 12) of 7 continuous and locals blocking sections, and then it can be concluded that in pipeline deposit distribution situation, and can be with
Find out the error between calculated value and practical caliber in reasonable range.
Table 1
| Position (km) | 45-47 | 47-49 | 49-51 | 51-53 | 53-55 | 55-57 | 57-59 |
| Practical caliber (mm) | 420.66 | 376.38 | 332.10 | 287.82 | 332.10 | 376.38 | 420.66 |
| It calculates caliber (mm) | 386.14 | 360.74 | 328.38 | 295.98 | 340.07 | 380.91 | 419.50 |
From above-mentioned multiple embodiments as can be seen that method provided by the invention compared with the existing technology in pressure wave
Detection method, other than the position that can determine that deposit, moreover it is possible to determine the length of deposit and the degree of blocking, have operation
Simply, the characteristics of good economy performance.
Referring to Fig.1 shown in 4, the present invention also provides a kind of system for detecting pipeline sediment distribution, which be can wrap
It includes:
Pressure detecting module 111, for since when incidence wave is inspired in specified duct section, in preset time
The fluid pressure data at the test point in the duct section is obtained in section；Mapping relations generation module 112, described in determining
Pressure value in fluid pressure data and the mapping relations between the corresponding sampling time；Parameter acquisition module 113, for according to
Mapping relations obtain arrival time and intensity of the incidence wave at the test point, and the incidence wave is in the duct section
Arrival time and intensity and the incidence wave of the pressurization wave that internal reflection is formed at the test point are in the duct section
Arrival time of the decompression wave reflected to form at the test point；First processing module 114, for according to the incidence wave
Arrival time and the arrival time of the pressurization wave, determine length of the deposit apart from the test point in the duct section；
Second processing module 115, for according to the arrival time of the pressurization wave and the arrival time of the decompression wave, it to be described heavy to determine
The length of product object；Third processing module 116, for according to intensity of the incidence wave at the test point, the pressurization wave
The length of intensity and the deposit apart from the test point at the test point, determines the blocking system of the duct section
Number range；Fourth processing module 117, the length for length, the deposit according to the deposit apart from the test point
Degree and the gross blockage factor range, determine deposit distribution situation in the duct section.
Disclosed herein as well is a kind of computer storage medias, store computer program, the computer program thereon
It is performed the steps of when being executed by processor
Since when incidence wave is inspired in specified duct section, obtained in the duct section within a preset period of time
Test point at fluid pressure data；Determine the mapping between the pressure value and corresponding sampling time in the fluid pressure data
Relationship；According to the mapping relations, arrival time and intensity of the incidence wave at the test point, the incidence wave are obtained
Exist in arrival time and intensity and the incidence wave of the pressurization wave that the duct section internal reflection is formed at the test point
Arrival time of the decompression wave that the duct section internal reflection is formed at the test point；According to the arrival time of the incidence wave
And the arrival time of the pressurization wave, determine length of the deposit apart from the test point in the duct section；According to described
It is pressurized the arrival time of wave and the arrival time of the decompression wave, determines the length of the deposit；Existed according to the incidence wave
The intensity and the deposit of intensity, the pressurization wave at the test point at the test point is apart from the test point
Length, determine the gross blockage factor range of the duct section；It is length according to the deposit apart from the test point, described heavy
The length and the gross blockage factor range of product object, determine deposit distribution situation in the duct section.
For convenience of description, it is divided into various modules when description apparatus above with function to describe respectively.Certainly, implementing this
The function of each module can be realized in the same or multiple software and or hardware when application.
It should be understood by those skilled in the art that, the embodiment of the present invention can provide as method, system or computer program
Product.Therefore, complete hardware embodiment, complete software embodiment or reality combining software and hardware aspects can be used in the present invention
Apply the form of example.Moreover, it wherein includes the computer of computer usable program code that the present invention, which can be used in one or more,
The computer program implemented in usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.) produces
The form of product.
The present invention be referring to according to the method for the embodiment of the present invention, the process of equipment (system) and computer program product
Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions
The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs
Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real
The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates,
Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or
The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one
The step of function of being specified in a box or multiple boxes.
In a typical configuration, calculating equipment includes one or more processors (CPU), input/output interface, net
Network interface and memory.
Memory may include the non-volatile memory in computer-readable medium, random access memory (RAM) and/or
The forms such as Nonvolatile memory, such as read-only memory (ROM) or flash memory (flash RAM).Memory is computer-readable medium
Example.
Computer-readable medium includes permanent and non-permanent, removable and non-removable media can be by any method
Or technology come realize information store.Information can be computer readable instructions, data structure, the module of program or other data.
The example of the storage medium of computer includes, but are not limited to phase change memory (PRAM), static random access memory (SRAM), moves
State random access memory (DRAM), other kinds of random access memory (RAM), read-only memory (ROM), electric erasable
Programmable read only memory (EEPROM), flash memory or other memory techniques, read-only disc read only memory (CD-ROM) (CD-ROM),
Digital versatile disc (DVD) or other optical storage, magnetic cassettes, tape magnetic disk storage or other magnetic storage devices
Or any other non-transmission medium, can be used for storage can be accessed by a computing device information.As defined in this article, it calculates
Machine readable medium does not include temporary computer readable media (transitory media), such as the data-signal and carrier wave of modulation.
It should also be noted that, the terms "include", "comprise" or its any other variant are intended to nonexcludability
It include so that the process, method, commodity or the equipment that include a series of elements not only include those elements, but also to wrap
Include other elements that are not explicitly listed, or further include for this process, method, commodity or equipment intrinsic want
Element.In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that including described want
There is also other identical elements in the process, method of element, commodity or equipment.
It will be understood by those skilled in the art that embodiments herein can provide as method, system or computer program product.
Therefore, complete hardware embodiment, complete software embodiment or embodiment combining software and hardware aspects can be used in the application
Form.It is deposited moreover, the application can be used to can be used in the computer that one or more wherein includes computer usable program code
The shape for the computer program product implemented on storage media (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.)
Formula.
The application can describe in the general context of computer-executable instructions executed by a computer, such as program
Module.Generally, program module includes routines performing specific tasks or implementing specific abstract data types, programs, objects, group
Part, data structure etc..The application can also be practiced in a distributed computing environment, in these distributed computing environments, by
Task is executed by the connected remote processing devices of communication network.In a distributed computing environment, program module can be with
In the local and remote computer storage media including storage equipment.
All the embodiments in this specification are described in a progressive manner, same and similar portion between each embodiment
Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.Especially for system reality
For applying example, since it is substantially similar to the method embodiment, so being described relatively simple, related place is referring to embodiment of the method
Part explanation.
The above description is only an example of the present application, is not intended to limit this application.For those skilled in the art
For, various changes and changes are possible in this application.All any modifications made within the spirit and principles of the present application are equal
Replacement, improvement etc., should be included within the scope of the claims of this application.
Claims (12)
1. a method of detection pipeline sediment distribution, which comprises the following steps:
Since when incidence wave is inspired in specified duct section, the inspection in the duct section is obtained within a preset period of time
Fluid pressure data at measuring point；
Determine the mapping relations between the pressure value and corresponding sampling time in the fluid pressure data；
According to the mapping relations, arrival time and intensity of the incidence wave at the test point, the incidence wave are obtained
Exist in arrival time and intensity and the incidence wave of the pressurization wave that the duct section internal reflection is formed at the test point
Arrival time of the decompression wave that the duct section internal reflection is formed at the test point；
According to the arrival time of the incidence wave and the arrival time of the pressurization wave, determine deposit in the duct section away from
Length from the test point；
According to the arrival time of the pressurization wave and the arrival time of the decompression wave, the length of the deposit is determined；
According to intensity at the test point of intensity of the incidence wave at the test point, the pressurization wave and described
Length of the deposit apart from the test point determines the gross blockage factor range of the duct section；
The length and the gross blockage factor range of length, the deposit according to the deposit apart from the test point,
Determine deposit distribution situation in the duct section；Wherein,
The intensity according to the incidence wave at the test point, the intensity of the pressurization wave at the test point and
Length of the deposit apart from the test point determines the gross blockage factor range of the duct section, comprising:
The case where closing valve for pipe end can obtain gross blockage factor range by solving following equation:
Wherein, k is local stoppages coefficient；F (k) is the function of local stoppages coefficient k；D for no local stoppages pipeline section caliber；c
For parameter related with pipe-supporting mode；HFFor the intensity of incidence wave；HfIntensity when test point is reached for pressurization wave；ΔhF
For incidence wave HFDecaying when being reached from deposit from test point；ΔhfDeclining when reaching test point from deposit for pressurization wave
Subtract；E is tubing elasticity modulus；E is pipe thickness.
2. the method for detection pipeline sediment distribution according to claim 1, which is characterized in that
According to the arrival time of the incidence wave and the arrival time of the pressurization wave, determine deposit in the duct section away from
Length from the test point, comprising:
Obtain pipeline elasticity modulus, pipe diameter, pipe thickness, fluid density and the elasticity of fluid system at the test point
Number；
According to formulaDetermine the average velocity of wave in the duct section；
Wherein, aijFor the average velocity of wave of i-th section of j-th of spatial mesh size；KijFor the elasticity of fluid system of i-th section of j-th of spatial mesh size
Number；ρ is fluid density；E is pipeline elasticity modulus；DiFor the pipe diameter of the i-th segment pipe；E is pipe thickness；C is and pipeline
The related constraint factor of section constraint condition；According to the arrival time and the average velocity of wave of the incidence wave and the pressurization wave
Determine length of the deposit apart from the test point.
3. the method for detection pipeline sediment distribution according to claim 2, which is characterized in that described according to the pressurization
The arrival time of wave and the decompression wave determines the length of the deposit, comprising:
The length of the deposit is determined according to the arrival time and the average velocity of wave of the pressurization wave and the decompression wave.
4. the method for detection pipeline sediment distribution according to claim 1, which is characterized in that described according to the incidence
Intensity, pressurization wave intensity and the deposit the test point at of the wave at the test point is apart from the inspection
The length of measuring point determines the gross blockage factor range of the duct section, comprising:
Length according to the deposit apart from the test point, and determine that the incidence wave reaches described sink by numerical simulation
Intensity when product object proximal end；
According to intensity of the incidence wave at the test point, and when reaching the deposit proximal end in conjunction with the incidence wave
Intensity determines the attenuation of the incidence wave；
The decaying of the pressurization wave is determined by friction drag loss at position according to the attenuation of the incidence wave, pressure wave
Amount；
According to the intensity of the attenuation of the incidence wave, the attenuation, the incidence wave of the pressurization wave at the test point
And intensity of the pressurization wave at the test point, determine the gross blockage factor of the estimation of the duct section, and then determine stifled
Fill in coefficient range.
5. the method for detection pipeline sediment distribution according to claim 4, which is characterized in that described according to the deposition
Length of the object apart from the test point, and it is strong when determining that the incidence wave reaches the deposit proximal end by numerical simulation
Degree, comprising:
Pass through formula: Progress numerical simulation determines strong when the incidence wave reaches the deposit proximal end
Degree；
Wherein, x is the distance along duct section centerline direction；H is head；V is the mean flow rate of fluid in duct section；G attaches most importance to
Power acceleration；JsFor friction loss；F is the coefficient of friction resistance in Brunone model；kBIt is constant；A is water-hammer wave speed；α inclines for pipeline section
Oblique angle；HxIt adjusts the distance the local derviation of x for head；HtIt is head to the local derviation of time t；VxIt adjusts the distance the local derviation of x for flow velocity；VtFor stream
Local derviation of the speed to time t；D for no local stoppages pipeline section caliber；fSFor the stable state coefficient of friction resistance.
6. the method for detection pipeline sediment distribution according to claim 4, which is characterized in that described according to the incidence
The attenuation of wave, pressure wave friction loss at determine the attenuation of the pressurization wave, comprising:
Pass through formula:Determine the attenuation of the pressurization wave；
Wherein, Δ hFAttenuation when reaching the deposit from the test point for the pressurization wave；ΔhfFor the incidence wave
Attenuation when test point is reached from the deposit；ΔHFThe deposit is reached from the test point for the pressurization wave
Friction loss when place；ΔHfFriction loss when test point is reached from the deposit for the incidence wave.
7. the method for detection pipeline sediment distribution according to claim 4, which is characterized in that described according to the incidence
The intensity and the pressurization wave of the attenuation of wave, the attenuation, the incidence wave of the pressurization wave at the test point exist
Intensity at the test point determines the gross blockage factor of the duct section, comprising:
Pass through formula: Determine gross blockage factor
Range；
Wherein, k is local stoppages coefficient；D for no local stoppages pipeline section caliber；C is parameter related with pipe-supporting mode；
HFFor the intensity of incidence wave；HfFor intensity of the incidence wave at the test point；ΔhFIt is the pressurization wave from the detection
Attenuation when point reaches at the deposit；ΔhfDecaying when test point is reached from the deposit for the incidence wave
Amount；C is parameter related with pipe-supporting mode；E is tubing elasticity modulus；E is pipe thickness.
8. the method for detection pipeline sediment distribution according to claim 1, which is characterized in that the incidence wave is by beating
On or off is closed pumping plant and is inspired, or is inspired by the valve body that closing is located at pipe downstream end；
When the incidence wave is inspired by turning on or off pumping plant, the test point is located at the pumping plant；
When the incidence wave is inspired by the valve body that closing is located at pipe downstream end, the test point is located at the valve body
Place.
9. the method for detection pipeline sediment distribution according to claim 1, which is characterized in that determine the preset time
The step of section, comprising:
Determine the maximum allowable working pressure of duct section to be detected；
According to the maximum allowable working pressure of duct section to be detected, maximum allowable detection range is determined；
Preset time period is determined according to maximum allowable detection range.
10. the method for detection pipeline sediment distribution according to claim 1, which is characterized in that the mapping relations packet
The relationship change curve of m- pressure when including.
11. a kind of system for detecting pipeline sediment distribution characterized by comprising
Pressure detecting module, for being obtained since when incidence wave is inspired in specified duct section within a preset period of time
Take the fluid pressure data at the test point in the duct section；
Graph generation module, for determining that the mapping between the pressure value in the fluid pressure data and corresponding sampling time is closed
System；
Parameter acquisition module obtains arrival time of the incidence wave at the test point and strong according to the mapping relations
Degree, arrival time and intensity of the pressurization wave that the incidence wave is formed in the duct section internal reflection at the test point, with
And arrival time of the decompression wave that is formed in the duct section internal reflection of the incidence wave at the test point；
First processing module, described in determining according to the arrival time of the incidence wave and the arrival time of the pressurization wave
Length of the deposit apart from the test point in duct section；
Second processing module, for determining the length of the deposit according to the arrival time of the pressurization wave and the decompression wave
Degree；
Third processing module, according to intensity of the incidence wave at the test point, the pressurization wave at the test point
Length apart from the test point of intensity and the deposit, determine the gross blockage factor range of the duct section；
Fourth processing module, the length of length, the deposit according to the deposit apart from the test point and described
Gross blockage factor range determines deposit distribution situation in the duct section；Wherein,
The intensity according to the incidence wave at the test point, the intensity of the pressurization wave at the test point and
Length of the deposit apart from the test point determines the gross blockage factor range of the duct section, comprising:
The case where closing valve for pipe end can obtain gross blockage factor range by solving following equation:
Wherein, k is local stoppages coefficient；F (k) is the function of local stoppages coefficient k；D for no local stoppages pipeline section caliber；c
For parameter related with pipe-supporting mode；HFFor the intensity of incidence wave；HfIntensity when test point is reached for pressurization wave；ΔhF
For incidence wave HFDecaying when being reached from deposit from test point；ΔhfDeclining when reaching test point from deposit for pressurization wave
Subtract；E is tubing elasticity modulus；E is pipe thickness.
12. a kind of computer storage media, stores computer program thereon, which is characterized in that the computer program is located
Reason device performs the steps of when executing
Since when incidence wave is inspired in specified duct section, the inspection in the duct section is obtained within a preset period of time
Fluid pressure data at measuring point；
Determine the mapping relations between the pressure value and corresponding sampling time in the fluid pressure data；
According to the mapping relations, arrival time and intensity of the incidence wave at the test point, the incidence wave are obtained
Exist in arrival time and intensity and the incidence wave of the pressurization wave that the duct section internal reflection is formed at the test point
Arrival time of the decompression wave that the duct section internal reflection is formed at the test point；
According to the arrival time of the incidence wave and the arrival time of the pressurization wave, determine deposit in the duct section away from
Length from the test point；According to the arrival time of the pressurization wave and the arrival time of the decompression wave, it is described heavy to determine
The length of product object；
According to intensity at the test point of intensity of the incidence wave at the test point, the pressurization wave and described
Length of the deposit apart from the test point determines the gross blockage factor range of the duct section；
The length and the gross blockage factor range of length, the deposit according to the deposit apart from the test point,
Further more accurately determine deposit distribution situation in the duct section；Wherein,
The intensity according to the incidence wave at the test point, the intensity of the pressurization wave at the test point and
Length of the deposit apart from the test point determines the gross blockage factor range of the duct section, comprising:
The case where closing valve for pipe end can obtain gross blockage factor range by solving following equation:
Wherein, k is local stoppages coefficient；F (k) is the function of local stoppages coefficient k；D for no local stoppages pipeline section caliber；c
For parameter related with pipe-supporting mode；HFFor the intensity of incidence wave；HfIntensity when test point is reached for pressurization wave；ΔhF
For incidence wave HFDecaying when being reached from deposit from test point；ΔhfDeclining when reaching test point from deposit for pressurization wave
Subtract；E is tubing elasticity modulus；E is pipe thickness.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710431352.6A CN107329167B (en) | 2017-06-09 | 2017-06-09 | Detect method, system and the computer storage media of pipeline sediment distribution |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710431352.6A CN107329167B (en) | 2017-06-09 | 2017-06-09 | Detect method, system and the computer storage media of pipeline sediment distribution |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107329167A CN107329167A (en) | 2017-11-07 |
| CN107329167B true CN107329167B (en) | 2019-05-07 |
Family
ID=60194668
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710431352.6A Active CN107329167B (en) | 2017-06-09 | 2017-06-09 | Detect method, system and the computer storage media of pipeline sediment distribution |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107329167B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109780447B (en) * | 2018-12-28 | 2020-11-06 | 大连理工大学 | Method for detecting blockage condition in pressure pipeline by using pulse pressure wave |
| CN112647929B (en) * | 2019-10-11 | 2024-05-14 | 中国石油化工股份有限公司 | Experimental device for be used for detecting pit shaft deposit |
| CN113970292B (en) * | 2021-10-09 | 2023-10-31 | 同济大学 | Method for identifying object surface attachment area |
| CN116446829A (en) * | 2023-06-16 | 2023-07-18 | 什邡慧丰采油机械有限责任公司 | Wellhead safety automatic control integrated system based on flow |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS59171877A (en) * | 1983-03-19 | 1984-09-28 | Jgc Corp | Method for detecting blocked position of pipeline |
| CN202040552U (en) * | 2010-08-12 | 2011-11-16 | 北京知容寓远软件科技有限公司 | System for monitoring the leakage or obstruction of pipeline base on low frequency sound wave |
| CN203162549U (en) * | 2012-11-26 | 2013-08-28 | 昆明理工大学 | Device for monitoring blockage of ore slurry pipeline |
| CN104034800A (en) * | 2014-03-25 | 2014-09-10 | 李光 | Assessment method and system for hydraulic detection of conveying pipeline and for state of carrier fluid pipeline |
| CN104930355A (en) * | 2015-06-09 | 2015-09-23 | 段焕丰 | Online nondestructive detection method and device applicable to urban water supply pipeline system |
| CN104989958A (en) * | 2015-07-30 | 2015-10-21 | 云南大红山管道有限公司 | Positioning method for blocking position of ore slurry conveying pipeline |
| CN105135216A (en) * | 2015-07-30 | 2015-12-09 | 云南大红山管道有限公司 | Method for locating leakage position of ore pulp conveying pipe |
| CN106090621A (en) * | 2016-06-02 | 2016-11-09 | 东北电力大学 | A kind of water supply network leakage, plugging fault diagnosis and localization method analyzed based on pressure signal |
-
2017
- 2017-06-09 CN CN201710431352.6A patent/CN107329167B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS59171877A (en) * | 1983-03-19 | 1984-09-28 | Jgc Corp | Method for detecting blocked position of pipeline |
| CN202040552U (en) * | 2010-08-12 | 2011-11-16 | 北京知容寓远软件科技有限公司 | System for monitoring the leakage or obstruction of pipeline base on low frequency sound wave |
| CN203162549U (en) * | 2012-11-26 | 2013-08-28 | 昆明理工大学 | Device for monitoring blockage of ore slurry pipeline |
| CN104034800A (en) * | 2014-03-25 | 2014-09-10 | 李光 | Assessment method and system for hydraulic detection of conveying pipeline and for state of carrier fluid pipeline |
| CN104930355A (en) * | 2015-06-09 | 2015-09-23 | 段焕丰 | Online nondestructive detection method and device applicable to urban water supply pipeline system |
| CN104989958A (en) * | 2015-07-30 | 2015-10-21 | 云南大红山管道有限公司 | Positioning method for blocking position of ore slurry conveying pipeline |
| CN105135216A (en) * | 2015-07-30 | 2015-12-09 | 云南大红山管道有限公司 | Method for locating leakage position of ore pulp conveying pipe |
| CN106090621A (en) * | 2016-06-02 | 2016-11-09 | 东北电力大学 | A kind of water supply network leakage, plugging fault diagnosis and localization method analyzed based on pressure signal |
Non-Patent Citations (5)
| Title |
|---|
| 压力传感器在给水管道阻力系数测试中的应用;魏希柱 等;《传感器技术》;20011231;第20卷(第4期);第36-38页 |
| 基于压力波法的管道堵塞检测技术;刘恩斌 等;《天然气工业》;20060430;第26卷(第4期);第112-114页 |
| 天然气管道堵塞检测理论方法研究;邓志彬;《中国博士学位论文全文数据库 工程科技I辑》;20170515(第05(2017)期);第B019-33页 |
| 煤层气输气管道内沉积物定位检测方法研究;周东平 等;《煤气与热力》;20100430;第30卷(第4期);第27-29页 |
| 输气管道清管堵塞定位技术研究;张红兵 等;《天然气技术》;20080229;第2卷(第1期);第42-43页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107329167A (en) | 2017-11-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107329167B (en) | Detect method, system and the computer storage media of pipeline sediment distribution | |
| CN109344436B (en) | Online simulation method for large complex natural gas pipe network system | |
| US11592351B2 (en) | Urban non-metallic pipeline leakage location method | |
| CN103814277B (en) | Measurement of rate of flow method and device | |
| CN109578816B (en) | Water supply pipe network pipe explosion detection method and device, control method and device | |
| CN103353908B (en) | A kind of pipe resistance coefficient Method for Accurate Calculation based on numerical computations | |
| CN108662442B (en) | The localization method and device of pipe leakage | |
| Nikpour et al. | Experimental and numerical simulation of water hammer | |
| CN107524921B (en) | Products Batch Transportation Pipeline batching interface tracking and device | |
| CN108197070A (en) | Natural gas not exclusively blocks pipeline method for numerical simulation | |
| CN106503396A (en) | The multidimensional Hydraulic Power System transition analogy method coupled with finite volume method based on finite difference calculus | |
| US20150143920A1 (en) | Method of assessing and condition monitoring of fluid conduits and apparatus therefor | |
| CN111102476B (en) | Pipe burst searching method based on three-dimensional underground pipeline | |
| CN113348349A (en) | Method and system for analyzing pipe condition | |
| CN110348112B (en) | Self-adaptive constraint method for calculating maximum flow rate of natural gas emptying pipeline | |
| CA3086168C (en) | Method and system for detecting and quantifying irregularities in a fluidic channel | |
| CN111125867B (en) | Method for establishing and calculating real-time transient model of chemical production pipeline based on chaotic particle swarm | |
| Jiang et al. | Study on pressure transients in low pressure water-hydraulic pipelines | |
| Vítkovský et al. | Minimization algorithms and experimental inverse transient leak detection | |
| CN108490804A (en) | Pipeline minute leakage detectability appraisal procedure under a kind of complex environment | |
| Dupont | Galerkin methods for modeling gas pipelines | |
| CN103353319A (en) | Moisture flow measuring method based on straight through type gas ultrasonic flowmeter | |
| Huang et al. | An Efficient Lagrangian-Based Approach for Transient Modelling and Analysis in Urban Water Distribution Networks:(055) | |
| Abdol-Hamid et al. | The effect of leakage on railroad brake pipe steady state behavior | |
| CN110705185A (en) | Method for predicting pipeline air hammer |
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 |