CN112229460B - Method and device for calculating mass flow rate of liquid pipeline leakage - Google Patents
Method and device for calculating mass flow rate of liquid pipeline leakage Download PDFInfo
- Publication number
- CN112229460B CN112229460B CN202010982849.9A CN202010982849A CN112229460B CN 112229460 B CN112229460 B CN 112229460B CN 202010982849 A CN202010982849 A CN 202010982849A CN 112229460 B CN112229460 B CN 112229460B
- Authority
- CN
- China
- Prior art keywords
- value
- initial
- average flow
- friction coefficient
- flow velocity
- 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
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/76—Devices for measuring mass flow of a fluid or a fluent solid material
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
- F17D—PIPE-LINE SYSTEMS; PIPE-LINES
- F17D5/00—Protection or supervision of installations
- F17D5/02—Preventing, monitoring, or locating loss
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M3/00—Investigating fluid-tightness of structures
- G01M3/02—Investigating fluid-tightness of structures by using fluid or vacuum
- G01M3/26—Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors
- G01M3/28—Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors for pipes, cables or tubes; for pipe joints or seals; for valves ; for welds
- G01M3/2807—Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors for pipes, cables or tubes; for pipe joints or seals; for valves ; for welds for pipes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Abstract
The embodiment of the application discloses a method and a device for calculating the mass flow rate of liquid pipeline leakage, which are used for improving the calculation efficiency of the mass flow rate of liquid pipeline leakage. The method comprises the following steps: acquiring a target energy formula; acquiring static parameters and real-time parameters; constructing a loss function about the leak according to the target energy, the static parameters and the real-time parameters; assigning the average flow velocity and the Vanning friction coefficient to an initial average flow velocity value and an initial Vanning friction coefficient value respectively, and calculating a loss value according to the initial average flow velocity value, the initial Vanning friction coefficient value and a loss function; judging whether the loss value converges or not, if so, determining the initial average flow velocity value as a target average flow velocity value; if not, updating the initial average flow velocity value and the initial van der waals friction coefficient value by a gradient descent method, and recalculating the loss value until the loss value converges, and then determining a target average flow velocity value; and calculating the mass flow rate according to the target average flow rate value.
Description
Technical Field
The embodiment of the application relates to the field of pipeline leakage detection, in particular to a method and a device for calculating the mass flow rate of liquid pipeline leakage.
Background
Liquid pipeline transportation is an important mode of cargo transportation, mainly for transporting liquids such as petroleum, seawater and the like. However, the pipeline itself must have a certain structural influence during transportation, for example: factors such as pipeline corrosion, natural damage, artificial damage, defects of the pipeline and the like can influence the safety problem of the pipeline, so that leakage occurs. When dangerous chemicals such as inflammable, explosive, severe toxicity and the like are conveyed by the liquid pipeline, once leakage occurs, if the leaked liquid pipeline is not treated as soon as possible, the leaked liquid pipeline can not only cause loss to economy of the liquid pipeline, but also cause major secondary disasters such as fire, explosion, poisoning, environmental pollution and the like.
The determination of the risk area of leakage of the liquid pipe has a direct relation to the leakage amount, and after the liquid pipe leaks, the mass flow rate of the leakage of the liquid pipe must be determined first to determine the repair scheme. Therefore, accurately calculating the mass flow rate is a necessary condition for repairing the leakage of the liquid pipeline, and can provide information for the result processing of accidents, avoid the occurrence of secondary accidents and reduce the economic and personnel safety loss generated in the process of remedying the leakage accidents.
In the conventional detection, the flow form of the liquid in the pipeline, the Reynolds number and the influence of the roughness coefficient of the net pipeline are considered, and the mass flow rate of the leakage part of the liquid pipeline is calculated according to the data, so that the detection in situ is performed manually by using a special instrument, the overall repair time is increased, and the calculation efficiency of the mass flow rate is reduced.
Disclosure of Invention
The embodiment of the application discloses a method and a device for calculating the mass flow rate of liquid pipeline leakage, which are used for improving the calculation efficiency of the mass flow rate of liquid pipeline leakage.
An embodiment of the present application provides a method for calculating a mass flow rate of a liquid pipeline leakage, including:
obtaining a target energy formula, wherein the target energy formula is a formula of liquid flow, and comprises two unknowns of average flow velocity and van der waals friction coefficient at a leakage position;
acquiring static parameters and real-time parameters, wherein the static parameters and the real-time parameters are respectively the structural parameters and the operation parameters of the liquid pipeline;
constructing a loss function for the leak based on the target energy, the static parameter, and the real-time parameter;
assigning an initial average flow velocity value and an initial van der waals friction coefficient value to the average flow velocity and the van der waals friction coefficient respectively, and calculating a loss value according to the initial average flow velocity value, the initial van der waals friction coefficient value and the loss function;
Judging whether the loss value converges or not, if so, determining that the initial average flow velocity value is a target average flow velocity value;
if not, updating the initial average flow velocity value and the initial van der waals friction coefficient value by a gradient descent method, and recalculating a loss value until the loss value converges, and then determining a target average flow velocity value;
and calculating the mass flow rate according to the target average flow rate value.
Optionally, the updating the initial average flow velocity value and the initial van der waals coefficient value by a gradient descent method includes:
respectively calculating a first gradient function and a second gradient function of the initial average flow velocity value and the initial van der waals friction coefficient value according to the loss function;
respectively calculating a first gradient value and a second gradient value of the initial average flow velocity value and the initial van der waals friction coefficient value according to the first gradient function and the second gradient function;
respectively calculating a first increment and a second increment of the initial average flow velocity value and the initial van der waals friction coefficient value according to the first gradient value and the second gradient value;
updating the initial average flow velocity value and the initial van der waals coefficient value according to the first increment and the second increment.
Optionally, the acquiring the target energy formula includes:
acquiring an initial energy formula, a friction coefficient term formula and a differential pressure loss term formula, wherein the initial energy formula comprises a friction coefficient term, the friction coefficient term formula comprises a differential pressure loss term, and the target energy formula comprises two unknowns of average flow velocity and van-to-n friction coefficient at a leakage position;
substituting the friction loss term formula into the friction coefficient term in the initial energy formula, and substituting the differential pressure loss term formula into the differential pressure loss term in the friction coefficient term formula to obtain a target energy formula.
Optionally, the acquiring the static parameter and the real-time parameter includes:
and acquiring static parameters and real-time parameters through the static database and the real-time database respectively.
Optionally, before the acquiring the target energy formula, the method further includes:
constructing a static database of the liquid pipeline, wherein the static database is used for storing structural parameters of the liquid pipeline;
acquiring the operation parameters of the liquid pipeline through a detection instrument;
and constructing a real-time database of the liquid pipeline according to the operation parameters.
A second aspect of the embodiments provides a device for calculating a mass flow rate of a liquid conduit leak, comprising:
The device comprises a first acquisition unit, a second acquisition unit and a third acquisition unit, wherein the first acquisition unit is used for acquiring a target energy formula, the target energy formula is a formula of liquid flow, and the target energy formula comprises two unknowns of average flow velocity and van-to-n friction coefficient at a leakage position;
the second acquisition unit is used for acquiring static parameters and real-time parameters, wherein the static parameters and the real-time parameters are respectively the structural parameters and the operation parameters of the liquid pipeline;
a first construction unit for constructing a loss function for the leak from the target energy, the static parameter and the real-time parameter;
the setting unit is used for respectively assigning the average flow velocity and the Vanning friction coefficient to an initial average flow velocity value and an initial Vanning friction coefficient value, and calculating a loss value according to the initial average flow velocity value, the initial Vanning friction coefficient value and the loss function;
a judging unit for judging whether the loss value converges or not;
a determining unit configured to determine the initial average flow velocity value as a target average flow velocity value when the judging unit confirms that the loss value converges;
a first updating unit, configured to update the initial average flow velocity value and the initial fanning friction coefficient value by a gradient descent method when the judging unit confirms that the loss value is not converged, and recalculate the loss value until the loss value is converged, and then determine a target average flow velocity value;
A first calculation unit for calculating a mass flow rate from the target average flow rate value.
Optionally, the first updating unit includes:
a second calculation module for calculating a first gradient function and a second gradient function of the initial average flow velocity value and the initial van der waals coefficient value according to the loss function, respectively;
a third calculation module for calculating a first gradient value and a second gradient value of the initial average flow velocity value and the initial fanning friction coefficient value according to the first gradient function and the second gradient function, respectively;
a fourth calculation module for calculating a first increment and a second increment of the initial average flow velocity value and the initial van der waals coefficient value according to the first gradient value and the second gradient value, respectively;
a second update module updates the initial average flow rate value and the initial van der Waals coefficient of friction value based on the first increment and the second increment.
Optionally, the first acquisition unit includes:
the third acquisition module is used for acquiring an initial energy formula, a friction coefficient term formula and a differential pressure loss term formula, wherein the initial energy formula comprises a friction coefficient term, the friction coefficient term formula comprises a differential pressure loss term, and the target energy formula comprises two unknowns of average flow velocity and van-to-n friction coefficient at a leakage position;
And a fourth obtaining module, configured to substitute the friction loss term formula into the friction coefficient term in the initial energy formula, and substitute the differential pressure loss term formula into the differential pressure loss term in the friction coefficient term formula, so as to obtain a target energy formula.
Optionally, the second obtaining unit includes:
and the fifth acquisition module is used for acquiring the static parameters and the real-time parameters through the static database and the real-time database respectively.
Optionally, the apparatus further comprises:
the second construction unit is used for constructing a static database of the liquid pipeline, and the static database is used for storing the structural parameters of the liquid pipeline;
a sixth acquisition unit for acquiring the operation parameters of the liquid pipeline through a detection instrument;
and the third construction unit is used for constructing a real-time database of the liquid pipeline according to the operation parameters.
A third aspect of the embodiments provides a device for calculating a mass flow rate of a liquid conduit leak, comprising:
the device comprises a processor, a memory, an input/output unit and a bus;
the processor is connected with the memory, the input/output unit and the bus;
the processor specifically performs the following operations:
Obtaining a target energy formula, wherein the target energy formula is a formula of liquid flow, and comprises two unknowns of average flow velocity and van der waals friction coefficient at a leakage position;
acquiring static parameters and real-time parameters, wherein the static parameters and the real-time parameters are respectively the structural parameters and the operation parameters of the liquid pipeline;
constructing a loss function for the leak based on the target energy, the static parameter, and the real-time parameter;
assigning an initial average flow velocity value and an initial van der waals friction coefficient value to the average flow velocity and the van der waals friction coefficient respectively, and calculating a loss value according to the initial average flow velocity value, the initial van der waals friction coefficient value and the loss function;
judging whether the loss value converges or not, if so, determining that the initial average flow velocity value is a target average flow velocity value;
if not, updating the initial average flow velocity value and the initial van der waals friction coefficient value by a gradient descent method, and recalculating a loss value until the loss value converges, and then determining a target average flow velocity value;
and calculating the mass flow rate according to the target average flow rate value.
Optionally, the processor is further configured to perform the operations of any of the alternatives in the first aspect.
A computer readable storage medium having stored thereon a program which when executed on a computer performs any of the alternative methods of the first aspect as well as the preceding first aspect.
From the above technical solutions, the embodiments of the present application have the following advantages:
and acquiring a formula (hereinafter referred to as target energy) of liquid flow of the liquid pipeline, and acquiring structural parameters of the liquid pipeline and operation parameters of the liquid pipeline during operation, so as to obtain a loss function of the liquid pipeline at a leakage position. Setting an initial value (an initial average flow velocity value and an initial Van Ning friction coefficient value) for the loss function, calculating a loss value according to the initial value and the loss function, judging whether the loss value is converged, if the loss value is not converged, updating the initial value through a gradient descent method, and calculating the loss value according to the updated initial value again until the loss value is converged. When the loss value converges, it can be confirmed that the current initial value is the optimal solution, i.e. closest to the actual situation, and finally the mass flow rate is calculated according to the optimal solution.
According to the method, the influence of the flowing form of the liquid in the pipeline, the Reynolds number and the rough coefficient of the net pipeline is not needed to be considered, the optimal solution can be calculated according to the loss function and the gradient descent algorithm only by obtaining the loss function of the pipeline according to the target energy type, the real-time parameters and the static parameters, the mass flow rate of the leakage part of the liquid pipeline is calculated through the optimal solution, the on-site detection is not needed to be carried out manually by using a special instrument, the overall repair time is reduced, and the calculation efficiency of the mass flow rate is improved.
Drawings
FIG. 1 is a flow chart of an embodiment of a method for calculating mass flow rate of liquid conduit leakage in an embodiment of the present application;
FIG. 2 is a flow chart of another embodiment of a method for calculating a mass flow rate of a liquid conduit leak in an embodiment of the present application;
FIG. 3 is a schematic diagram of an embodiment of a computing device for mass flow rate of liquid conduit leakage in an embodiment of the present application;
FIG. 4 is a schematic diagram of another embodiment of a computing device for mass flow rate of liquid conduit leakage in an embodiment of the present application;
FIG. 5 is a schematic diagram of another embodiment of a computing device for mass flow rate of liquid conduit leakage in an embodiment of the present application.
Detailed Description
In order to better understand the technical solutions of the present invention, the following description will clearly and completely describe the technical solutions of the embodiments of the present invention with reference to the drawings in the embodiments of the present invention, and it is apparent that the described embodiments are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, shall fall within the scope of the invention.
The embodiment of the application discloses a method and a device for calculating the mass flow rate of liquid pipeline leakage, which are used for improving the calculation efficiency of the mass flow rate of liquid pipeline leakage.
In this embodiment, the method for calculating the mass flow rate of the liquid pipeline leakage may be implemented in a system, or may be implemented in a terminal, and is not specifically limited. For convenience of description, embodiments of the present application will be described using a server as an example of an execution body.
Referring to fig. 1, an embodiment of a method for calculating a mass flow rate of a liquid pipeline leakage in an embodiment of the present application includes:
101. a target energy formula of a system, wherein the target energy formula is a formula of liquid flow and comprises two unknowns of average flow velocity and van-to-n friction coefficient at a leakage position;
the system obtains a target energy formula, which is a motion formula of the liquid, for calculating data of the liquid in the liquid pipe, and can calculate data related to the mass flow rate through the target energy formula, and further calculate the mass flow rate through the data.
Alternatively, in the present embodiment, the target energy formula is formula 1 as follows:
Δp—specify the pressure difference between the initial point and the pipe break in units: pa.
ρ—liquid density, unit: kg/m 3 。
Alpha, a dimensionless rate profile correction coefficient, for laminar flow, alpha takes 0.5; take 1.0 for turbulence.
g-gravitational acceleration, unit: kg/s 2 。
Δh—specify the height difference between the initial point and the pipe break, unit: m.
W s Shaft work, unit: pa.m.
A-pipeline sectional area, unit: m.
K-pressure differential loss caused by friction of the pipeline or pipeline fittings, and is dimensionless.
f-Vanning coefficient of friction, dimensionless.
L-length of pipe from initial point to pipe break, unit: m.
d-inner diameter of pipe, unit: m.
v 1 -designatingAverage instantaneous flow rate of starting point, unit: m/s;
v 2 average instantaneous flow rate at pipe break, unit: m/s.
102. The system acquires static parameters and real-time parameters, wherein the static parameters and the real-time parameters are respectively the structural parameters and the operation parameters of the liquid pipeline;
the system obtains relevant parameters in the target energy formula, including static parameters and real-time parameters.
The static parameters are structural parameters of the pipeline, and different detection modes are used for detecting the structural parameters on each liquid pipeline, wherein the structural parameters comprise fixed parameters such as the length, the diameter and the pipeline sectional area of the liquid pipeline. In this embodiment, only the static parameters required to be used in the formula are obtained from all the static parameters.
The static parameters may be detected in various ways, and may be obtained through machine measurement, manual measurement, or from a corresponding instruction for purchasing a liquid pipeline, which is not limited herein.
The real-time parameters are real-time dynamic data of liquid in the pipeline obtained by using different detection modes when the liquid pipeline starts to work, and the real-time parameters comprise: liquid flow rate, liquid pressure, liquid density, etc. at a specified point of the liquid conduit. From all the real-time parameters, only the real-time parameters needed to be used in the formula are obtained.
The real-time parameters may be detected in various ways, and may be obtained through machine measurement or manual measurement, which is not limited herein.
103. The system constructs a loss function about the leak according to the target energy, the static parameter and the real-time parameter;
the system acquires corresponding known parameters from all static parameters in real time according to the meaning of each parameter in the formula, substitutes the parameters into a target energy formula, and defines a loss function according to the energy formula.
Alternatively, in the present embodiment, the target energy formula (formula 1) is defined as a value of v by a predetermined rule 2 And f is the loss function of the unknown A number. Optionally, in this embodiment, the loss function defined by the target energy formula, the static parameter, and the parameter through a certain rule is the following formula 2:
the loss function is determined by f and v 2 The loss value of the loss function can be obtained, and the current van der waals friction coefficient and the average flow velocity at the leakage position of the liquid pipeline are determined according to the loss value.
104. The system respectively assigns the average flow velocity and the Van Nile friction coefficient to an initial average flow velocity value and an initial Van Nile friction coefficient value, and calculates a loss value according to the initial average flow velocity value, the initial Van Nile friction coefficient value and the loss function;
in order to calculate the accurate fanning friction coefficient and the average flow velocity at the leakage position of the liquid pipeline so as to enable the loss value to achieve convergence, according to the operation mode of the gradient descent method, an initial value needs to be defined for the fanning friction coefficient and the average flow velocity respectively, so that the loss value can be obtained according to the initial average flow velocity and the initial fanning friction coefficient.
Alternatively, in the present embodiment, the initial average flow velocity is set to be equal to the initial point v 1 Equal values, the initial van der waals coefficient of friction is assigned a random value greater than 0 and less than 1.
105. The system judges whether the loss value is converged or not; if yes, go to step 106; if not, go to step 107;
And substituting the initial average flow velocity and the initial Vanning friction coefficient into a loss function to obtain a loss value, and determining whether the loss value at the moment is converged by judging whether the loss value reaches a preset threshold value or not.
Optionally, in this embodiment, the preset threshold for determining that the loss value converges is 0, and the error is smaller than 0.00001.
106. When the system judges that the loss value is converged, the system determines that the initial average flow velocity value is a target average flow velocity value;
when the loss value reaches the preset threshold value, the current initial average flow rate and the initial fanning friction coefficient can be determined to be the optimal solution, namely the actual value which is the most consistent with the current leakage state of the liquid pipeline can be used as the basic parameter for calculating the mass flow rate.
107. When the system judges that the loss value is not converged, the system updates the initial average flow velocity value and the initial fanning friction coefficient value through a gradient descent method, and recalculates the loss value until the loss value is converged, and then the target average flow velocity value is determined;
and when the loss value does not reach the range of the preset threshold value, determining that the current initial average flow velocity value and the initial van der waals friction coefficient value are not close to the actual value, namely the optimal solution is not achieved. Only if the loss value calculated from the fanning friction coefficient value and the initial average flow rate value reaches a preset threshold value, it is determined that the calculated initial average flow rate value of the current leak point substantially reaches the actual initial average flow rate value. If the current calculated value does not accord with the calculated value, a new initial average flow velocity value and an initial van der waals friction coefficient value are needed to be reselected to calculate a loss value, namely the current initial average flow velocity and the initial van der waals friction coefficient value are needed to be updated under a certain rule, so that the loss value reaches a preset threshold value. Optionally, in this embodiment, the gradient descent method is used to iterate such that the current initial average flow rate and the initial fanning friction coefficient are updated. The detailed updating step is shown in the corresponding embodiment of fig. 2. This step is followed by step 104, where the loss value is recalculated.
108. The system calculates a mass flow rate from the target average flow rate value.
When the initial average flow rate and the initial fanning friction coefficient of the optimal solution are determined, selecting the average flow rate of the optimal solution to calculate the mass flow rate, and calculating the mass flow rate by using a formula 3:
wherein:
m—mass flow rate, unit: kg/s.
ρ—liquid density, unit: .
A-pipeline sectional area, unit: .
v 1 -mean instantaneous flow rate at a given starting point, unit: m/s;
v 2 average instantaneous flow rate at pipe break, unit: m/s.
In this embodiment, the loss function may be defined by defining a target energy formula, a static parameter and a real-time parameter, determining a loss value of the loss function, updating the initial average flow rate and the initial fanning friction coefficient by a gradient descent method, continuously obtaining an initial average flow rate at which the loss value reaches convergence, and calculating the mass flow rate by the initial average flow rate.
The mass flow rate can be calculated by acquiring basic data (average instantaneous flow rate of starting point, pressure difference shaft work of initial point and pipeline fracture) which can be measured in real time by the instrument and acquiring structural parameters in the liquid pipeline. The influence of the flowing form of the liquid in the pipeline, the Reynolds number and the roughness coefficient of the net pipeline does not need to be considered, so that the time for calculating the measurement data required by the mass flow rate is reduced, and the efficiency is improved.
In the above embodiments, the static parameters are obtained in various manners on the market, and in this application, a database is constructed first, so that the parameters can be obtained from the database more quickly. The manner in which the database is constructed is described in detail below.
Referring to fig. 2, another embodiment of a method for calculating a mass flow rate of a liquid pipeline leak in an embodiment of the present application includes:
201. the system constructs a static database of the liquid pipeline, wherein the static database is used for storing structural parameters of the liquid pipeline;
by recording data of each liquid pipeline, the structural parameters of all the liquid pipelines are stored in a static database.
Optionally, in this embodiment, in the process of purchasing the assembled pipes, data arrangement needs to be performed on parameters such as length, diameter, turning angle and the like of each liquid pipe, so as to ensure that detailed structural data of any liquid pipe exists in the database.
202. The system acquires the operation parameters of the liquid pipeline through a detection instrument;
by placing the instrument inside the liquid conduit, such that measurable data can be obtained while the liquid conduit is transporting liquid.
Alternatively, in this embodiment, the data of the liquid in the pipe during operation is measured by the instrument by installing a flow rate measuring device, a pressure measuring device, a liquid temperature measuring instrument, and the like in the pipe. Even if the liquid pipe leaks, the instrument performs measurement, and when abnormal data is measured, an abnormal alarm is generated.
203. The system builds a real-time database of the liquid pipeline according to the operation parameters;
and receiving real-time data fed back by the measuring instrument through the instrument for measuring in the pipeline, constructing a real-time database, and detecting the running state of the liquid pipeline at the moment. In case of leakage of the liquid pipe, the measurement frequency is quickened, so that the latest data are updated in the real-time database.
204. The system acquires an initial energy formula, a friction coefficient term formula and a differential pressure loss term formula, wherein the initial energy formula comprises a friction coefficient term, and the friction coefficient term formula comprises a differential pressure loss term;
the initial energy formula is obtained, and the energy formula selected in this embodiment is as follows formula 4:
Δp—specify the pressure difference between the initial point and the pipe break in units: pa.
ρ—liquid density, unit: kg/m 3 。
Δu—specify the average instantaneous flow rate difference of the liquid at the initial point and at the pipe break, in units: m/s.
Alpha, a dimensionless rate profile correction coefficient, for laminar flow, alpha takes 0.5; take 1.0 for turbulence.
g-gravitational acceleration, unit: kg/s 2 。
Δh—specify the height difference between the initial point and the pipe break, unit: m.
F-mechanical energy loss due to friction, is a function of the Van Ning coefficient of friction, unit: m.N/kg.
W s Shaft work, unit: pa.m.
A-pipeline sectional area, unit: m.
u—average instantaneous flow rate of liquid, unit: m/s.
In this embodiment, the friction coefficient term formula is as follows formula 5:
u—average instantaneous flow rate of liquid, unit: m/s.
K-pressure differential loss caused by friction of the pipeline or pipeline fittings, and is dimensionless.
n-total number of pipes and pipe fittings.
In this embodiment, the selected differential pressure loss term formula is as follows formula 6:
f-Vanning coefficient of friction, dimensionless.
L-length of pipe from initial point to pipe break, unit: m.
d-inner diameter of pipe, unit: m.
205. Substituting the friction loss term formula into the friction coefficient term in the initial energy formula by a system, and substituting the differential pressure loss term formula into the differential pressure loss term in the friction coefficient term formula to obtain a target energy formula;
substituting equation 5 and equation 6 into equation 4 yields the target energy equation (equation 1).
The system obtains a target energy formula, which is a motion formula of the liquid, for calculating data of the liquid in the liquid pipe, and can calculate data related to the mass flow rate through the target energy formula, and further calculate the mass flow rate through the data.
206. The system acquires static parameters and real-time parameters through a static database and a real-time database respectively;
the system provides the data required by the retrieval of the static database and the real-time database, and provides the data of each liquid pipeline designed by the database, so that formula data can be obtained more quickly in the calculation process.
207. The system constructs a loss function about the leak according to the target energy, the static parameter and the real-time parameter;
208. the system respectively assigns the average flow velocity and the Van Nile friction coefficient to an initial average flow velocity value and an initial Van Nile friction coefficient value, and calculates a loss value according to the initial average flow velocity value, the initial Van Nile friction coefficient value and the loss function;
209. the system judges whether the loss value is converged or not; if yes, go to step 210; if not, go to step 211;
210. when the system judges that the loss value is converged, the system determines that the initial average flow velocity value is a target average flow velocity value;
steps 207 and 210 in this embodiment are similar to steps 103 and 106 in the previous embodiment, and will not be repeated here.
211. When the system judges that the loss value is not converged, the system calculates a first gradient function and a second gradient function of the initial average flow velocity value and the initial van der waals friction coefficient value according to the loss function;
When the system judges that the loss value is not converged, the system firstly conducts derivation on the initial average flow velocity and the initial van der waals friction coefficient according to the loss function, and a first gradient function of the loss function relative to the initial average flow velocity and a second gradient function of the loss function relative to the initial van der waals friction coefficient value are obtained.
The first gradient function and the second gradient function are:
first gradient function
Second gradient function
212. The system calculates a first gradient value and a second gradient value of the initial average flow velocity value and the initial van der waals friction coefficient value according to the first gradient function and the second gradient function respectively;
the system calculates a gradient value from the gradient function.
Calculating a first gradient value of the initial average flow rate by a first gradient function, the first gradient value being expressed as:
first gradient value
Calculating a second gradient value of the initial average flow rate by a second gradient function, the second gradient value being expressed as:
second gradient value
n-nth iteration.
The gradient value is used to calculate a first increment and a second increment of the initial average flow rate corresponding to the initial van der waals coefficient of friction.
213. The system calculates a first increment and a second increment of the initial average flow velocity value and the initial van der waals friction coefficient value according to the first gradient value and the second gradient value respectively;
The system calculates a first increment and a second increment of the initial average flow velocity and the initial van der waals friction coefficient respectively according to the gradient value.
Calculating a first increment of the initial average flow rate from the first gradient value, the first increment being expressed as:
first increment
Calculating a second increment of the initial average flow rate from the second gradient value, the second increment being expressed as:
second increment
n-nth iteration.
η—learning rate or iteration step in gradient descent method.
The next group of initial average flow velocity and initial van ning friction coefficient can be calculated through increment.
214. A system updates the initial average flow velocity value and the initial van der waals coefficient value according to the first increment and the second increment;
the next group of initial average flow velocity and initial van ning friction coefficient can be calculated through increment.
v 2 (n+1)＝v 2 (n)+Δv 2 (n)
f(n+1)＝f(n)+Δf(n)
After updating the initial average flow rate and initial fanning friction coefficient, step 208 is performed to recalculate the loss value.
215. The system calculates a mass flow rate from the target average flow rate value.
Step 215 in this embodiment is similar to step 108 in the previous embodiment, and will not be described again here.
In this embodiment, the loss function may be defined by defining a target energy formula, a static parameter and a real-time parameter, determining a loss value of the loss function, updating the initial average flow rate and the initial fanning friction coefficient by a gradient descent method, continuously obtaining an initial average flow rate at which the loss value reaches convergence, and calculating the mass flow rate by the initial average flow rate.
The mass flow rate can be calculated by acquiring basic data (average instantaneous flow rate of starting point, pressure difference shaft work of initial point and pipeline fracture) which can be measured in real time by the instrument and acquiring structural parameters in the liquid pipeline. The influence of the flowing form of the liquid in the pipeline, the Reynolds number and the roughness coefficient of the net pipeline does not need to be considered, so that the time for calculating the measurement data required by the mass flow rate is reduced, and the efficiency is improved.
In addition, in the embodiment, a database construction mode is adopted, so that data can be extracted more quickly, the time for detecting the data in the repairing process is reduced, and the repairing efficiency is improved.
The method for calculating the mass flow rate of the liquid pipe leakage in the embodiment of the present application is described above, and the device for calculating the mass flow rate of the liquid pipe leakage in the embodiment of the present application will be described in detail below with reference to the accompanying drawings:
referring to fig. 3, an embodiment of a device for calculating a mass flow rate of a liquid conduit leak in an embodiment of the present application includes:
a first obtaining unit 301, configured to obtain a target energy formula, where the target energy formula is a formula of liquid flow, and the target energy formula includes two unknowns of an average flow velocity and a fanning friction coefficient at a leakage position;
A second obtaining unit 302, configured to obtain a static parameter and a real-time parameter, where the static parameter and the real-time parameter are a structural parameter and an operation parameter of the liquid pipeline respectively;
a first construction unit 303 for constructing a loss function for the leak from the target energy formula, the static parameter and the real-time parameter;
a setting unit 304, configured to assign an initial average flow velocity value and an initial fanning friction coefficient value to the average flow velocity and the fanning friction coefficient, respectively, and calculate a loss value according to the initial average flow velocity value, the initial fanning friction coefficient value and the loss function;
a judging unit 305 for judging whether the loss value converges;
a determining unit 306 configured to determine the initial average flow velocity value as a target average flow velocity value when the judging unit 305 confirms that the loss value converges;
a first updating unit 307, configured to update the initial average flow velocity value and the initial fanning friction coefficient value by a gradient descent method when the judging unit 305 confirms that the loss value is not converged, and recalculate the loss value until the loss value is converged, and then determine a target average flow velocity value;
A first calculation unit 308 for calculating a mass flow rate from the target average flow rate value.
Referring to fig. 4, another embodiment of a computing device for mass flow rate of liquid conduit leakage in an embodiment of the present application includes:
a second construction unit 401, configured to construct a static database of liquid pipes, where the static database is used to store structural parameters of the liquid pipes;
a sixth acquisition unit 402, configured to acquire an operation parameter of the liquid pipeline through a detection instrument;
a third construction unit 403 for constructing a real-time database of the liquid pipeline according to the operating parameters
A first obtaining unit 404, configured to obtain a target energy formula, where the target energy formula is a formula of liquid flow, and the target energy formula includes two unknowns of an average flow rate and a fanning friction coefficient at a leakage point;
in this embodiment, the first acquiring unit 404 includes a third acquiring module 4041 and a fourth acquiring module 4042.
The third obtaining module 4041 is configured to obtain an initial energy formula, a friction coefficient term formula, and a differential pressure loss term formula, where the initial energy formula includes a friction coefficient term, and the friction coefficient term formula includes a differential pressure loss term;
A fourth obtaining module 4042, configured to substitute the friction loss term formula into the friction coefficient term in the initial energy formula, and substitute the differential pressure loss term formula into the differential pressure loss term in the friction coefficient term formula, so as to obtain a target energy formula;
a second obtaining unit 405, configured to obtain a static parameter and a real-time parameter, where the static parameter and the real-time parameter are a structural parameter and an operation parameter of the liquid pipeline respectively;
in this embodiment, the second acquiring unit 405 includes a fifth acquiring module 4051.
The fifth obtaining module 4051 is configured to obtain the static parameter and the real-time parameter through the static database and the real-time database, respectively.
A first construction unit 406 for constructing a loss function for the leak from the target energy, the static parameters and the real-time parameters;
a setting unit 407, configured to assign the average flow velocity and the fanning friction coefficient to an initial average flow velocity value and an initial fanning friction coefficient value, respectively, and calculate a loss value according to the initial average flow velocity value, the initial fanning friction coefficient value, and the loss function;
a judging unit 408, configured to judge whether the loss value converges;
A determining unit 409 for determining the initial average flow velocity value as a target average flow velocity value when the judging unit 408 confirms that the loss value converges;
a first updating unit 410, configured to update the initial average flow velocity value and the initial fanning friction coefficient value by a gradient descent method when the judging unit 408 confirms that the loss value is not converged, and recalculate the loss value until the loss value is converged, and then determine a target average flow velocity value;
in the present embodiment, the first updating unit 410 includes a second computing module 4101, a third computing module 4102, a fourth computing module 4103, and a second updating module 4104.
A second calculation module 4101 for calculating a first gradient function and a second gradient function of the initial average flow velocity value and the initial fanning friction coefficient value, respectively, according to the loss function;
a third calculation module 4102 for calculating a first gradient value and a second gradient value of the initial average flow velocity value and the initial fanning friction coefficient value according to the first gradient function and the second gradient function, respectively;
a fourth calculation module 4103 for calculating a first increment and a second increment of the initial average flow velocity value and the initial fanning friction coefficient value according to the first gradient value and the second gradient value, respectively;
A second updating module 4104, configured to update the initial average flow velocity value and the initial fanning friction coefficient value according to the first increment and the second increment, and after the unit finishes executing, recalculate the loss value, that is, execute the step of setting the unit 407;
a first calculation unit 411 for calculating a mass flow rate from the optimal value.
Referring to fig. 5, another embodiment of a computing device for mass flow rate of liquid conduit leakage in an embodiment of the present application includes:
a processor 501, a memory 502, an input-output unit 503, and a bus 504;
the processor 501 is connected to the memory 502, the input/output unit 503, and the bus 504;
the processor 501 specifically performs the following operations:
obtaining a target energy formula, wherein the target energy formula is a formula of liquid flow, and comprises two unknowns of average flow velocity and van der waals friction coefficient at a leakage position;
acquiring static parameters and real-time parameters, wherein the static parameters and the real-time parameters are respectively the structural parameters and the operation parameters of the liquid pipeline;
constructing a loss function for the leak based on the target energy, the static parameter, and the real-time parameter;
assigning an initial average flow velocity value and an initial van der waals friction coefficient value to the average flow velocity and the van der waals friction coefficient respectively, and calculating a loss value according to the initial average flow velocity value, the initial van der waals friction coefficient value and the loss function;
Judging whether the loss value converges or not, if so, determining that the initial average flow velocity value is a target average flow velocity value;
if not, updating the initial average flow velocity value and the initial van der waals friction coefficient value by a gradient descent method, and recalculating a loss value until the loss value converges, and then determining a target average flow velocity value;
and calculating the mass flow rate according to the target average flow rate value.
In this embodiment, the functions of the processor 501 correspond to the steps in the embodiments shown in fig. 1 to 2, and are not described herein.
It will be clear to those skilled in the art that, for convenience and brevity of description, specific working procedures of the above-described systems, apparatuses and units may refer to corresponding procedures in the foregoing method embodiments, which are not repeated herein.
In the several embodiments provided in this application, it should be understood that the disclosed systems, apparatuses, and methods may be implemented in other ways. For example, the apparatus embodiments described above are merely illustrative, e.g., the division of the units is merely a logical function division, and there may be additional divisions when actually implemented, e.g., multiple units or components may be combined or integrated into another system, or some features may be omitted or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed with each other may be an indirect coupling or communication connection via some interfaces, devices or units, which may be in electrical, mechanical or other form.
The units described as separate units may or may not be physically separate, and units shown as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
In addition, each functional unit in each embodiment of the present application may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit. The integrated units may be implemented in hardware or in software functional units.
The integrated units, if implemented in the form of software functional units and sold or used as stand-alone products, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present application may be embodied in essence or a part contributing to the prior art or all or part of the technical solution in the form of a software product stored in a storage medium, including several instructions to cause a computer device (which may be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the methods described in the embodiments of the present application. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a read-only memory (ROM), a random access memory (RAM, random access memory), a magnetic disk, or an optical disk, or other various media capable of storing program codes.
Claims (10)
1. A method of calculating a mass flow rate of a liquid conduit leak, comprising:
obtaining a target energy formula, wherein the target energy formula is a formula of liquid flow, and comprises two unknowns of average flow velocity and van der waals friction coefficient at a leakage position;
acquiring static parameters and real-time parameters, wherein the static parameters and the real-time parameters are respectively the structural parameters and the operation parameters of the liquid pipeline;
constructing a loss function for the leak based on the target energy, the static parameter, and the real-time parameter;
assigning an initial average flow velocity value and an initial van der waals friction coefficient value to the average flow velocity and the van der waals friction coefficient respectively, and calculating a loss value according to the initial average flow velocity value, the initial van der waals friction coefficient value and the loss function;
judging whether the loss value converges or not, if so, determining that the initial average flow velocity value is a target average flow velocity value;
if not, updating the initial average flow velocity value and the initial van der waals friction coefficient value by a gradient descent method, and recalculating a loss value until the loss value converges, and then determining a target average flow velocity value;
and calculating the mass flow rate according to the target average flow rate value.
2. The method of calculating according to claim 1, wherein said updating said initial average flow velocity value and said initial van der waals coefficient of friction value by a gradient descent method comprises:
respectively calculating a first gradient function and a second gradient function of the initial average flow velocity value and the initial van der waals friction coefficient value according to the loss function;
respectively calculating a first gradient value and a second gradient value of the initial average flow velocity value and the initial van der waals friction coefficient value according to the first gradient function and the second gradient function;
respectively calculating a first increment and a second increment of the initial average flow velocity value and the initial van der waals friction coefficient value according to the first gradient value and the second gradient value;
updating the initial average flow velocity value and the initial van der waals coefficient value according to the first increment and the second increment.
3. The computing method of claim 1, wherein the obtaining the target energy formula comprises:
acquiring an initial energy formula, a friction coefficient term formula and a differential pressure loss term formula, wherein the initial energy formula comprises a friction coefficient term, the friction coefficient term formula comprises a differential pressure loss term, and the target energy formula comprises two unknowns of average flow velocity and van-to-n friction coefficient at a leakage position;
Substituting the friction coefficient term formula into the friction coefficient term in the initial energy formula, and substituting the differential pressure loss term formula into the differential pressure loss term in the friction coefficient term formula to obtain a target energy formula.
4. A computing method according to any one of claims 1 to 3, wherein the acquiring static parameters and real-time parameters comprises:
and acquiring static parameters and real-time parameters through the static database and the real-time database respectively.
5. The computing method of claim 4, wherein prior to the acquiring the target energy formula, the method further comprises:
constructing a static database of the liquid pipeline, wherein the static database is used for storing structural parameters of the liquid pipeline;
acquiring the operation parameters of the liquid pipeline through a detection instrument;
and constructing a real-time database of the liquid pipeline according to the operation parameters.
6. A computing device for mass flow rate of liquid conduit leakage, comprising:
the device comprises a first acquisition unit, a second acquisition unit and a third acquisition unit, wherein the first acquisition unit is used for acquiring a target energy formula, the target energy formula is a formula of liquid flow, and the target energy formula comprises two unknowns of average flow velocity and van-to-n friction coefficient at a leakage position;
The second acquisition unit is used for acquiring static parameters and real-time parameters, wherein the static parameters and the real-time parameters are respectively the structural parameters and the operation parameters of the liquid pipeline;
a first construction unit for constructing a loss function for the leak from the target energy, the static parameter and the real-time parameter;
the setting unit is used for respectively assigning the average flow velocity and the Vanning friction coefficient to an initial average flow velocity value and an initial Vanning friction coefficient value, and calculating a loss value according to the initial average flow velocity value, the initial Vanning friction coefficient value and the loss function;
a judging unit for judging whether the loss value converges or not;
a determining unit configured to determine the initial average flow velocity value as a target average flow velocity value when the judging unit confirms that the loss value converges;
a first updating unit, configured to update the initial average flow velocity value and the initial fanning friction coefficient value by a gradient descent method, and recalculate a loss value until the loss value converges, and then determine a target average flow velocity value;
a first calculation unit for calculating a mass flow rate from the target average flow rate value.
7. The computing device of claim 6, wherein the first updating unit comprises:
A second calculation module for calculating a first gradient function and a second gradient function of the initial average flow velocity value and the initial van der waals coefficient value according to the loss function, respectively;
a third calculation module for calculating a first gradient value and a second gradient value of the initial average flow velocity value and the initial fanning friction coefficient value according to the first gradient function and the second gradient function, respectively;
a fourth calculation module for calculating a first increment and a second increment of the initial average flow velocity value and the initial van der waals coefficient value according to the first gradient value and the second gradient value, respectively;
a second update module updates the initial average flow rate value and the initial van der Waals coefficient of friction value based on the first increment and the second increment.
8. The computing device of claim 6, wherein the first acquisition unit comprises:
the third acquisition module is used for acquiring an initial energy formula, a friction coefficient term formula and a differential pressure loss term formula, wherein the initial energy formula comprises a friction coefficient term, the friction coefficient term formula comprises a differential pressure loss term, and the target energy formula comprises two unknowns of average flow velocity and van-to-n friction coefficient at a leakage position;
And a fourth obtaining module, configured to substitute the friction coefficient term formula into the friction coefficient term in the initial energy formula, and substitute the differential pressure loss term formula into the differential pressure loss term in the friction coefficient term formula, so as to obtain a target energy formula.
9. The computing device of any of claims 6 to 8, wherein the second acquisition unit comprises:
and the fifth acquisition module is used for acquiring the static parameters and the real-time parameters through the static database and the real-time database respectively.
10. The computing device of claim 9, wherein the device further comprises:
the second construction unit is used for constructing a static database of the liquid pipeline, and the static database is used for storing the structural parameters of the liquid pipeline;
a sixth acquisition unit for acquiring the operation parameters of the liquid pipeline through a detection instrument;
and the third construction unit is used for constructing a real-time database of the liquid pipeline according to the operation parameters.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010982849.9A CN112229460B (en) | 2020-09-17 | 2020-09-17 | Method and device for calculating mass flow rate of liquid pipeline leakage |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010982849.9A CN112229460B (en) | 2020-09-17 | 2020-09-17 | Method and device for calculating mass flow rate of liquid pipeline leakage |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112229460A CN112229460A (en) | 2021-01-15 |
| CN112229460B true CN112229460B (en) | 2023-07-25 |
Family
ID=74108404
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010982849.9A Active CN112229460B (en) | 2020-09-17 | 2020-09-17 | Method and device for calculating mass flow rate of liquid pipeline leakage |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112229460B (en) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2018204723A1 (en) * | 2017-05-04 | 2018-11-08 | The Board Of Trustees Of The Leland Stanford Junior University | Leaky lamb wave flowmeter |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4726219A (en) * | 1986-02-13 | 1988-02-23 | Atlantic Richfield Company | Method and system for determining fluid pressures in wellbores and tubular conduits |
| US6970808B2 (en) * | 2004-04-29 | 2005-11-29 | Kingsley E. Abhulimen | Realtime computer assisted leak detection/location reporting and inventory loss monitoring system of pipeline network systems |
| CN102853261A (en) * | 2011-06-27 | 2013-01-02 | 国际商业机器公司 | Method and device for determining leakage amount of fluid in conveying pipeline |
| EP3112823A1 (en) * | 2015-07-03 | 2017-01-04 | Kamstrup A/S | System for monitoring a utility network |
| CN106650120B (en) * | 2016-12-27 | 2019-08-30 | 中国石油大学(北京) | A kind of liquid hydrocarbon pipe leakage method for measuring and calculating |
| CN107463737B (en) * | 2017-07-25 | 2019-12-06 | 中国石油大学(北京) | method and device for calculating leakage of liquid pipeline |
| CN108591836B (en) * | 2018-04-13 | 2020-06-26 | 中国石油大学(北京) | Method and device for detecting pipeline leakage |
-
2020
- 2020-09-17 CN CN202010982849.9A patent/CN112229460B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2018204723A1 (en) * | 2017-05-04 | 2018-11-08 | The Board Of Trustees Of The Leland Stanford Junior University | Leaky lamb wave flowmeter |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112229460A (en) | 2021-01-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hasan et al. | Probability assessment of burst limit state due to internal corrosion | |
| Teixeira et al. | Reliability analysis of pipelines with local corrosion defects under external pressure | |
| US5355324A (en) | Detecting leakage of fluid from a conduit | |
| CN112966378B (en) | Hydrogen leakage prediction method and system based on safety evaluation model | |
| Timashev et al. | Methods of assessing integrity of pipeline systems with different types of defects | |
| US10101194B2 (en) | System and method for identifying and recovering from a temporary sensor failure | |
| EP3594653B1 (en) | Diagnosis cost output device, diagnosis cost output method, and diagnosis cost output | |
| Witek et al. | Influence of gas transmission network failure on security of supply | |
| CN112229460B (en) | Method and device for calculating mass flow rate of liquid pipeline leakage | |
| WO2009022148A1 (en) | Method and apparatus for operating a pipeline | |
| CN114201840B (en) | Pipe section stress corrosion cracking risk identification method, device, equipment and medium | |
| CN117009801A (en) | Traffic prediction model training method, detection method, system and storage medium | |
| Zhang et al. | Reliability-Based Assessment of Cracked Pipelines Using Monte Carlo Simulation Technique With CorLAS™ | |
| Shakhmatov et al. | Safety evaluation of pipelines based on probabilistic approaches | |
| Kwestarz et al. | Method for leak detection and location for gas networks | |
| Riaz et al. | Volumetric flow and pressure gradient-based leak detection system for oil and gas pipelines | |
| CN105868536A (en) | Test data based health state assessment method for solid rocket engine | |
| CN102046955B (en) | Method and apparatus for monitoring for leaks in a stage ii fuel vapor recovery system | |
| Ergezinger et al. | Application of Noise Filtering Techniques for the Quantification of Uncertainty in Dent Strain Calculations | |
| KR102627849B1 (en) | Pipe prioriry inspection line selection system and method | |
| CN108345705B (en) | Method and device for evaluating pipeline leakage consequence influence area | |
| Paulin Jr et al. | Economic Impact of Current SIF, Flexibility, Inspection and Manufacturing Changes As They Relate to the B31 Piping Codes | |
| Godoy et al. | Plastic buckling of fixed-roof oil storage tanks under blast loads | |
| CN116698234A (en) | Pressure early warning method and system for pressure pipeline of power station boiler | |
| Adianto et al. | Demonstration of Limit States Design Method for Assessment of Corrosion and Crack Features |
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 |