CN115146513B - Simulation early warning method and system for pipeline leakage collapse - Google Patents
Simulation early warning method and system for pipeline leakage collapse Download PDFInfo
- Publication number
- CN115146513B CN115146513B CN202210886058.5A CN202210886058A CN115146513B CN 115146513 B CN115146513 B CN 115146513B CN 202210886058 A CN202210886058 A CN 202210886058A CN 115146513 B CN115146513 B CN 115146513B
- Authority
- CN
- China
- Prior art keywords
- road
- simulation
- collapse
- value
- medium fluid
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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 invention relates to a pipeline leakage collapse simulation early warning method and a pipeline leakage collapse simulation early warning system, which are characterized by firstly collecting simulation data of actual working conditions; secondly, establishing a discrete element-finite difference-particle flow dynamic and static coupling simulation model according to actual working condition simulation data; thirdly, performing a pipeline leakage collapse simulation test by using the simulation model to obtain a simulation test result; thirdly, obtaining a limit discrimination value according to the simulation test result; and finally, judging whether the actual working condition has the risk of road collapse based on the limit judgment value. The method can quantitatively represent the road collapse risk based on the actual working condition parameters, and has the advantages of comprehensive simulation elements, convenient parameter acquisition and accurate simulation result.
Description
Technical Field
The invention belongs to the technical field of road collapse test simulation, and particularly relates to a pipeline leakage collapse simulation early warning method and system.
Background
In the prior art, underground infrastructure is often in an overload running state, and the damage condition of the underground infrastructure is often caused because the construction standard of the existing underground infrastructure does not meet the current urban development requirement. The damage of underground pipelines is a typical case, namely, the underground pipelines are mostly cement concrete pipelines, the leakage phenomenon of pipeline interfaces is frequent under the long-time operation, and because the existing detection capability cannot meet the requirement of quick detection of the mileage diseases of the whole pipeline, all pipeline leakage points cannot be accurately positioned in a short time, the pipeline leakage erodes the soil body, the underground diseases are formed, and the road collapse is caused under the action of road surface load.
At present, the research in the aspect of numerical simulation in the field mainly adopts a discrete element method, the method can research the influence rule of pipeline leakage on a soil body stress field and a displacement field under the static or quasi-static condition of road surface load, but neglects the dynamic disturbance generated by the operation of a medium inside a pipeline on the peripheral soil body and the impact load caused to a leakage port, and the dynamic disturbance and dynamic scouring action of the pipeline often plays a promoting role in the formation of collapse hidden danger (especially the pipeline under pressure).
Therefore, it is necessary to explore a road collapse dynamic-static multiphase flow coupling numerical simulation algorithm.
Disclosure of Invention
In order to overcome the problems in the prior art, the invention provides a pipeline leakage collapse simulation early warning method and a pipeline leakage collapse simulation early warning system, which are used for solving the problems in the prior art.
A simulation early warning method for pipeline leakage collapse comprises the following steps:
s1, collecting simulation data of actual working conditions;
s2, establishing a discrete element-finite difference-particle flow dynamic and static coupling simulation model according to simulation data of actual working conditions;
s3, performing a pipeline leakage collapse simulation test by using the simulation model to obtain a simulation test result;
s4, obtaining a limit discrimination value according to the simulation test result;
and S5, judging whether the actual working condition has the risk of road collapse based on the limit judgment value.
In the above aspect and any possible implementation manner, there is further provided an implementation manner, in the step S1, simulation data of actual working conditions is collected through a field investigation or actual measurement method, where the simulation data includes modeling parameters and environmental variables;
the modeling parameters comprise the thickness, density, bulk modulus, shear modulus, internal friction angle, cohesion and uniaxial tensile strength of the road structure layer and the stratum;
the environmental variables include road load and flow rate of the media fluid.
The above-described aspect and any possible implementation manner further provide an implementation manner, and the step S2 includes:
s21, simulating the non-continuous characteristics of the road structure layer and the stratum, the near-homogeneous characteristics of the underground pipeline and the unsteady flow characteristics of the medium fluid;
s22, determining boundary conditions, transition areas, information link modes and transition area widths of the simulation model;
s23, establishing the discrete element-finite difference-particle flow dynamic and static coupling model according to the simulation data, the step S21 and the step S22.
The above aspects and any possible implementation manners further provide an implementation manner, wherein the simulation model is a rectangular space domain composed of 1 top surface and a bottom surface, and 4 side surfaces, and the boundary conditions are that the top surface applies a vehicle load, and the bottom surface and the 4 side surfaces apply a fixed displacement; the transition region comprises a road structure layer-stratum transition region, a stratum-pipeline transition region, a pipeline-medium fluid transition region and a stratum-medium fluid transition region.
In the aspect and any possible implementation manner described above, an implementation manner is further provided, where the road structure layer-stratum transition region information link manner is that a road structure layer is transmitted to a stratum in a one-way manner;
the stratum-pipeline transition domain information link mode is that stratum and pipeline are transmitted in two directions;
the formation-medium fluid transition region information link mode is that medium fluid is transmitted to the formation in a one-way mode.
The above aspects and any possible implementation manners further provide an implementation manner that the information transmitted by the medium fluid to the stratum in one direction is fluid impact load and impact vibration stress, and the fluid impact load is
In the above formula, σ a For fluid impact loads, ρ 0 Is the density of the medium fluid, v is the flow rate of the medium fluid;
the specific expression of the impact vibration stress is as follows:
in the above formula, σ b For shock vibration stress, ρ 1 Is the formation density, V is the formation volume, V 0 Is the wave velocity, L 2 Is a certain in the formationA point is spaced from the leak port.
The above-described aspect and any possible implementation manner further provide an implementation manner, where step S3 includes:
s31, flow velocity experiment of medium fluid: making the road surface load as a fixed value (a + b)/2, making the change interval of the medium fluid flow velocity as [ c, d ], and simulating the influence of the medium fluid flow velocity change on the road collapse by taking the unit flow velocity as a change gradient;
s32, road surface load experiment: the flow velocity of the medium fluid is set to be a constant value (c + d)/2, the change interval of the road surface load is set to be [ a, b ], and the influence of the road surface load change on the road collapse is simulated by taking the unit road surface load as the change gradient;
wherein a and b are respectively the minimum value and the maximum value of the road surface load, c and d are respectively the minimum value and the maximum value of the flow velocity of the medium fluid, and a, b, c and d are all obtained by the step S1.
In the aspect and any possible implementation manner described above, an implementation manner is further provided, where the limit criterion value in the step S4 includes a first limit criterion value and a second limit criterion value; the step S4 includes:
s41, if the road collapse occurs in the experimental processes of the step S31 and the step S32, outputting a first limit judgment value as a critical value of the flow rate of the medium fluid when the road collapse occurs in the step S31, and outputting a second limit judgment value as a critical value of the road surface load when the road collapse occurs in the step S32;
s42, if the road collapse does not occur in the experimental processes of the step S31 and the step S32, the output first limit judgment value and the output second limit judgment value are both-1;
s43, if the road collapse occurs in the step S31 and the road collapse does not occur in the step S32, outputting a first limit judgment value as a critical value of the flow rate of the medium fluid when the road collapse occurs in the step S31, and a second limit judgment value as-1;
and S44, if the road collapse does not occur in the step S31 and the road collapse occurs in the step S32, outputting a first limit judgment value of-1 and a second limit judgment value which is a critical value of the road surface load when the road collapse occurs in the step S32.
The above-described aspect and any possible implementation manner further provide an implementation manner, and the step S5 includes:
if the first limit judgment value and the second limit judgment value are both-1, judging that the actual working condition has no road collapse risk;
and if the first limit judgment value or the second limit judgment value is the other value, judging that the road collapse risk exists in the actual working condition.
The invention also provides a pipeline leakage collapse simulation early warning system, which comprises:
the acquisition module is used for acquiring analog simulation data in actual working conditions;
the model construction module is used for constructing a discrete element-finite difference-particle flow dynamic and static coupling model according to the simulation data;
the simulation module is used for carrying out pipeline leakage collapse simulation test by adopting the simulation model to obtain a simulation test result;
the result output module is used for outputting a limit discrimination value according to the simulation test result;
and the collapse early warning module is used for judging whether the actual working condition has the risk of road collapse according to the limit judgment value.
The invention has the advantages of
Compared with the prior art, the invention has the following beneficial effects:
according to the method and the system for simulating and early warning pipeline collapse caused by leakage, firstly, simulation data of actual working conditions are collected; secondly, establishing a discrete element-finite difference-particle flow dynamic and static coupling simulation model according to the simulation data of the actual working condition; thirdly, performing a pipeline leakage collapse simulation test by using the simulation model to obtain a simulation test result; thirdly, obtaining a limit discrimination value according to the simulation test result; and finally, judging whether the actual working condition has the risk of road collapse based on the limit judgment value. The method can quantitatively represent the road collapse risk based on the actual working condition parameters, and has the advantages of comprehensive simulation elements, convenient parameter acquisition and accurate simulation result.
Drawings
FIG. 1 is a schematic diagram of the process steps of the present invention.
Detailed Description
In order to better understand the technical solution of the present invention, the present disclosure includes but is not limited to the following detailed description, and similar techniques and methods should be considered as within the scope of the present invention. In order to make the technical problems, technical solutions and advantages of the present invention more apparent, the following detailed description is given with reference to the accompanying drawings and specific embodiments.
It should be understood that the described embodiments of the invention are only some embodiments of the invention, and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The terminology used in the embodiments of the invention is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used in the examples of the present invention and the appended claims, the singular forms "a," "an," and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise.
As shown in fig. 1, a simulation early warning method for pipeline leakage and collapse includes the following steps:
s1, collecting actual working condition simulation data through a field investigation or actual measurement method;
s2, establishing a discrete element-finite difference-particle flow dynamic and static coupling simulation model;
s3, carrying out a pipeline leakage collapse simulation test by adopting a control variable method;
s4, obtaining a limit discrimination value according to a simulation test result; (ii) a
And S5, judging whether the actual working condition has the risk of road collapse based on the output limit judgment value.
Preferably, in the present invention, the simulation data in step S1 includes modeling parameters and environmental variables; the method divides simulation data into modeling parameters and environment variables, is beneficial to distinguishing the influence degree of different parameters on the road collapse, the modeling parameters are basic physical and mechanical parameters which are difficult to change in a short time, the influence on the road collapse is small, and the environment variables can fluctuate violently in a short time and have important influence on the occurrence of road collapse accidents.
The modeling parameters comprise the thickness, density, volume modulus, shear modulus, internal friction angle, cohesion and uniaxial tensile strength of a road structure layer and a stratum, and the road structure layer is a combination of a surface layer and a base layer;
the modeling parameters further comprise the inner diameter, the outer diameter, the burial depth, the density, the bulk modulus, the cohesion, the internal friction angle, the shear modulus, the uniaxial tensile strength, the radius and the position of a leakage port and the density of a medium fluid;
the environment variables comprise road surface loads, the road surface load values are located in an interval [ a, b ], a is the minimum value of the road surface loads, and b is the maximum value of the road surface loads;
the environment variables also comprise the flow velocity of the medium fluid, the flow velocity value of the medium fluid is positioned in an interval [ c, d ], c is the minimum value of the flow velocity of the medium fluid, and d is the maximum value of the flow velocity of the medium fluid;
preferably, step S2 in the present invention includes:
s21, simulating the discontinuous characteristics of a road structure layer and a stratum by adopting a discrete element algorithm, simulating the near-homogeneous characteristics of an underground pipeline by adopting a finite difference algorithm, and simulating the unsteady flow characteristics of a medium fluid by adopting a particle flow algorithm.
And S22, determining the boundary condition, the transition domain, the information link mode and the transition domain width of the simulation model, wherein the step is favorable for improving the convergence speed of the model and the reliability of the simulation result by reasonably determining the boundary condition, the information link mode and the transition domain width of the simulation model.
The simulation model is a rectangular space domain consisting of 1 top surface, a bottom surface and 4 side surfaces, and the boundary conditions are that the top surface applies vehicle load and the bottom surface and the 4 side surfaces are fixedly displaced; the boundary condition refers to the initial condition of the calculation simulation model, and can be stress, displacement or other physical parameters.
The transition region comprises a road structure layer-stratum transition region, a stratum-pipeline transition region, a pipeline-medium fluid transition region and a stratum-medium fluid transition region.
The information link mode of the road structure layer-stratum transition region is that the road structure layer is transmitted to the stratum in a one-way mode; the transmission information of the road structure layer-stratum transition region is first stress and first deformation, and the unidirectional transmission of the stress and the deformation from the road structure layer to the stratum is favorable for feeding back the action state of the road structure layer subjected to vehicle-mounted on the road surface to the stratum in real time; the first stress and the first deformation are obtained through iteration of a discrete element algorithm of the road structure layer;
the information link mode of the road stratum-pipeline transition domain is bidirectional transmission of stratum and pipeline; the information transmitted by the stratum to the pipeline is second stress; the information fed back and transmitted to the stratum by the pipeline is second deformation, stress in the stratum directly acts on the surface of the pipeline in an actual working condition, the pipeline is likely to deform, meanwhile, the stress of the peripheral stratum is changed due to the deformation of the pipeline, a mutual feedback response relationship exists between the pipeline and the peripheral stratum, and the two-way transmission of the stratum and the pipeline can effectively represent the response relationship; the second stress is obtained through discrete element algorithm iteration of the stratum, and the second deformation is obtained through finite difference algorithm iteration of the pipeline;
the pipeline-medium fluid transition region does not generate information transmission;
the information link mode of the stratum-medium fluid transition region is that medium fluid is transmitted to the stratum in a one-way mode; the information transmitted to the stratum by the medium fluid is fluid impact load and impact vibration stress, the medium leakage in the flow pipeline can not only generate power impact action on the stratum, but also generate vibration waves when acting on the impact load of a leakage hole, the vibration waves are gradually reduced along with the increase of the distance, and the process that the medium fluid transmits the fluid impact load and the impact vibration stress to the stratum in a one-way mode can be accurately described; the fluid impact load and the impact vibration stress are obtained by embedding a specific expression of the fluid impact load and the impact vibration stress into a particle flow algorithm for calculation;
the specific expression of the fluid impact load is as follows:
in the above formula, σ a For fluid impact loads, ρ 0 Is the density of the medium fluid, v is the flow velocity of the medium fluid, p 0 And v are both known values;
the specific expression of the impact vibration stress is as follows:
in the above formula, σ b For shock vibration stress, ρ 1 Is the formation density, V is the formation volume, V 0 Taking 2480m/s, L as wave velocity 2 The distance between a certain point in the stratum and the leakage hole; rho 1 、V、v 0 And L 2 Is a known value;
the widths of the road structure layer-stratum transition region, the stratum-pipeline transition region, the pipeline-medium fluid transition region and the stratum-medium fluid transition region are all smaller than 1/100 of the minimum side length of the rectangular space region, different transition regions are arranged to enable simulation model parameters to be uniformly transmitted, iteration precision of a simulation model is improved, but the convergence time of the model is prolonged when the width of the transition region is larger than 1/100 of the minimum side length of the rectangular space region according to experimental statistics, the convergence time of the model is remarkably increased, and therefore the width critical value of the transition region is set to be 1/100 of the minimum side length of the rectangular space region.
S23, establishing a discrete element-finite difference-particle flow dynamic-static coupling model according to the analog simulation data obtained in the step 1, the analog simulation algorithm determined in the step S21, the boundary conditions, the transition domain information link mode and the transition domain width determined in the step S22;
the dynamic in the discrete element-finite difference-particle flow dynamic and static coupling model refers to the dynamic impact action of medium fluid on the stratum, and the static in the discrete element-finite difference-particle flow dynamic and static coupling model refers to the static mechanical action among the road structure layer, the stratum and the pipeline;
preferably, the step S3 in the present invention includes:
s31, flow velocity experiment of medium fluid: making the road surface load as a fixed value (a + b)/2, making the change interval of the medium fluid flow velocity as [ c, d ], and simulating the influence of the medium fluid flow velocity change on the road collapse by taking the unit flow velocity as a change gradient;
s32, road surface load experiment: the flow velocity of the medium fluid is set to be a constant value (c + d)/2, the change interval of the road surface load is set to be [ a, b ], and the influence of the road surface load change on the road collapse is simulated by taking the unit road surface load as the change gradient;
wherein a and b are respectively the minimum value and the maximum value of the road surface load, c and d are respectively the minimum value and the maximum value of the flow velocity of the medium fluid, and a, b, c and d are all obtained by the step S1. That is, when one variable changes, the other variable needs to be controlled to be a fixed value.
Preferably, in the present invention, the limit determination value in the step S4 includes a first limit determination value and a second limit determination value;
the step S4 includes:
s41, if the road collapses in the test processes of the step S31 and the step S32, outputting a first limit judgment value as a critical value of the flow rate of the medium fluid when the road collapses in the step S31, and outputting a second limit judgment value as a critical value of the road surface load when the road collapses in the step S32;
s42, if the road collapse does not occur in the test processes of the step S31 and the step S32, outputting a first limit judgment value and a second limit judgment value which are both-1;
s43, if the road collapse occurs in the step S31 and the road collapse does not occur in the step S32, outputting a first limit judgment value as a critical value of the flow rate of the medium fluid when the road collapse occurs in the step S31, and setting a second limit judgment value as-1;
s44, if the road collapse does not occur in the step S31 and the road collapse occurs in the step S32, outputting a first limit judgment value of-1 and a second limit judgment value which is a critical value of the road surface load when the road collapse occurs in the step S32;
the critical judgment criterion of the simulation model for the occurrence of the road collapse is that the maximum equivalent strain occurs at the bottom of a road structure layer; the maximum equivalent strain is obtained through iteration of a discrete element algorithm of a road structure layer; the first stress, the second stress, the first deformation and the second deformation before the maximum equivalent strain are obtained by inputting boundary conditions, a transition region information connection mode and a transition region thickness in calculation.
Preferably, step S5 in the present invention includes:
if the first limit judgment value and the second limit judgment value are both-1, judging that the actual working condition has no road collapse risk;
if at least one of the first limit discrimination value and the second limit discrimination value is not-1, namely the first limit discrimination value and the second limit discrimination value are not-1, the risk of road collapse in the actual working condition is judged under the condition that the first limit discrimination value and the second limit discrimination value are other values, the risk early warning value is 80% of the non-negative limit discrimination value, and the non-negative limit discrimination value is the maximum value of the first limit discrimination value and the second limit discrimination value. The safety redundancy of 20% is set, so that the rationality of the risk early warning value under the limit condition is improved.
The invention provides a discrete element-finite difference method-particle flow coupling numerical model establishing method according to discrete characteristics of a road structure layer and a stratum, near-uniform characteristics of a pipeline and particle characteristics of a medium fluid, completes the setting of a link mode and boundary conditions of a boundary transition domain, and realizes the establishment of a pipeline leakage collapse simulation model based on dynamic and static load coupling effects. Under the condition that 20% of safety redundancy is reserved, the invention determines the early warning value of the collapse risk under the actual working condition according to the calculation result, and the data of the required actual working condition are all the basic mechanics of a road structure body, a stratum, a pipeline and medium fluid.
Preferably, the invention also provides a pipeline leakage collapse simulation early warning system, which comprises:
the acquisition module is used for acquiring analog simulation data in actual working conditions;
the model construction module is used for constructing a discrete element-finite difference-particle flow dynamic and static coupling model according to the simulation data;
the simulation module is used for carrying out pipeline leakage collapse simulation test by adopting the simulation model to obtain a simulation test result;
the result output module is used for outputting a limit discrimination value according to the simulation test result;
and the collapse early warning module is used for judging whether the actual working condition has the risk of road collapse according to the limit judgment value. And parameters are acquired conveniently.
The foregoing description shows and describes several preferred embodiments of the present invention, but as aforementioned, it is to be understood that the invention is not limited to the forms disclosed herein, and is not to be construed as excluding other embodiments, and that the invention is capable of use in various other combinations, modifications, and environments and is capable of changes within the scope of the inventive concept as expressed herein, commensurate with the above teachings, or the skill or knowledge of the relevant art. And that modifications and variations may be effected by those skilled in the art without departing from the spirit and scope of the invention as defined by the appended claims.
Claims (8)
1. A simulation early warning method for pipeline leakage collapse is characterized by comprising the following steps: the method comprises the following steps:
s1, collecting simulation data of actual working conditions: collecting actual working condition simulation data by a field investigation or actual measurement method, wherein the simulation data comprises modeling parameters and environment variables;
the modeling parameters comprise the thickness, density, bulk modulus, shear modulus, internal friction angle, cohesion and uniaxial tensile strength of the road structure layer and the stratum;
the environmental variables comprise road surface load and flow rate of medium fluid;
s2, establishing a discrete element-finite difference-particle flow dynamic and static coupling simulation model according to actual working condition simulation data, wherein the discrete element-finite difference-particle flow dynamic and static coupling simulation model comprises the following steps:
s21, simulating discontinuous characteristics of the road structure layer and the stratum, simulating near-homogeneous characteristics of the underground pipeline and simulating unsteady flow characteristics of medium fluid;
s22, determining boundary conditions, transition areas, information link modes and transition area widths of the simulation model;
s23, establishing a discrete element-finite difference-particle flow dynamic and static coupling model according to the simulation data, the step S21 and the step S22;
s3, performing pipeline leakage collapse simulation test by using the simulation model to obtain a simulation test result;
s4, obtaining a limit discrimination value according to the simulation test result;
and S5, judging whether the actual working condition has the risk of road collapse based on the limit judgment value.
2. The early warning method according to claim 1, wherein the simulation model is a rectangular space domain consisting of 1 top surface, a bottom surface and 4 side surfaces, and the boundary conditions are that a vehicle load is applied to the top surface and a fixed displacement is applied to the bottom surface and the 4 side surfaces; the transition region comprises a road structure layer-stratum transition region, a stratum-pipeline transition region, a pipeline-medium fluid transition region and a stratum-medium fluid transition region.
3. The early warning method as claimed in claim 1, wherein the road structure layer-stratum transition region information link mode is one-way transmission from the road structure layer to the stratum;
the stratum-pipeline transition domain information link mode is that stratum and pipeline are transmitted in two directions;
the formation-medium fluid transition region information link mode is that medium fluid is transmitted to the formation in a one-way mode.
4. The early warning method according to claim 3, wherein the information transmitted by the medium fluid to the stratum in one direction is fluid impact load and impact vibration stress, and the specific expression of the fluid impact load is as follows:
in the above expression, σ a For fluid impact loading, ρ 0 Is the density of the medium fluid, v is the flow rate of the medium fluid;
the specific expression of the impact vibration stress is as follows:
in the above expression, σ b For shock vibration stress, ρ 1 Is the formation density, V is the formation volume, V 0 Is the wave velocity, L 2 Is the distance from a point in the formation to the leak-off port.
5. The warning method according to claim 1, wherein the step S3 comprises:
s31, medium fluid flow rate experiment: the road surface load is set to be a constant value (a + b)/2, the variation interval of the flow velocity of the medium fluid is set to be [ c, d ], and the influence of the variation of the flow velocity of the medium fluid on the road collapse is simulated by taking the unit flow velocity as the variation gradient;
s32, road surface load experiment: the flow velocity of the medium fluid is set to be a constant value (c + d)/2, the change interval of the road surface load is set to be [ a, b ], and the influence of the road surface load change on the road collapse is simulated by taking the unit road surface load as the change gradient;
wherein a and b are respectively the minimum value and the maximum value of the road surface load, c and d are respectively the minimum value and the maximum value of the flow velocity of the medium fluid, and a, b, c and d are all obtained by the step S1.
6. The warning method according to claim 1, wherein the limit discrimination value in the step S4 includes a first limit discrimination value and a second limit discrimination value;
the step S4 includes:
s41, if the road collapse occurs in the experimental processes of the step S31 and the step S32, outputting a first limit discrimination value as a critical value of the flow rate of the medium fluid when the road collapse occurs in the step S31, and outputting a second limit discrimination value as a critical value of the road surface load when the road collapse occurs in the step S32;
s42, if the road collapse does not occur in the experimental processes of the step S31 and the step S32, the output first limit judgment value and the output second limit judgment value are both-1;
s43, if the road collapse occurs in the step S31 and the road collapse does not occur in the step S32, outputting a first limit judgment value as a critical value of the flow rate of the medium fluid when the road collapse occurs in the step S31, and a second limit judgment value as-1;
and S44, if the road collapse does not occur in the step S31 and the road collapse occurs in the step S32, outputting a first limit judgment value of-1 and a second limit judgment value which is a critical value of the road surface load when the road collapse occurs in the step S32.
7. The warning method as claimed in claim 6, wherein the step S5 comprises:
if the first limit judgment value and the second limit judgment value are both-1, judging that the actual working condition has no road collapse risk;
and if the first limit judgment value or the second limit judgment value is the other value, judging that the road collapse risk exists in the actual working condition.
8. A pipeline leakage collapse simulation early warning system is characterized by comprising:
the acquisition module is used for acquiring analog simulation data in actual working conditions: collecting actual working condition simulation data by a field investigation or actual measurement method, wherein the simulation data comprises modeling parameters and environment variables;
the modeling parameters comprise the thickness, density, bulk modulus, shear modulus, internal friction angle, cohesion and uniaxial tensile strength of the road structure layer and the stratum;
the environmental variables comprise road surface load and flow rate of medium fluid;
the model construction module is used for constructing a discrete element-finite difference-particle flow dynamic and static coupling model according to analog simulation data, and comprises the following steps:
s21, simulating discontinuous characteristics of the road structure layer and the stratum, simulating near-homogeneous characteristics of the underground pipeline and simulating unsteady flow characteristics of medium fluid;
s22, determining boundary conditions, transition areas, information link modes and transition area widths of the simulation model;
s23, establishing a discrete element-finite difference-particle flow dynamic and static coupling model according to the simulation data, the step S21 and the step S22;
the simulation module is used for carrying out pipeline leakage collapse simulation test by adopting the model to obtain a simulation test result;
the result output module is used for outputting a limit discrimination value according to the simulation test result;
and the collapse early warning module is used for judging whether the actual working condition has the risk of road collapse according to the limit judgment value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210886058.5A CN115146513B (en) | 2022-07-26 | 2022-07-26 | Simulation early warning method and system for pipeline leakage collapse |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210886058.5A CN115146513B (en) | 2022-07-26 | 2022-07-26 | Simulation early warning method and system for pipeline leakage collapse |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115146513A CN115146513A (en) | 2022-10-04 |
| CN115146513B true CN115146513B (en) | 2023-04-07 |
Family
ID=83414527
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210886058.5A Active CN115146513B (en) | 2022-07-26 | 2022-07-26 | Simulation early warning method and system for pipeline leakage collapse |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115146513B (en) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114218651A (en) * | 2021-12-22 | 2022-03-22 | 中国地质大学（武汉） | Method and equipment for simulating particle discrete element in sandy soil foundation slab load test |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2017044074A1 (en) * | 2015-09-08 | 2017-03-16 | Halliburton Energy Services, Inc. | Domain-adaptive hydraulic fracture simulators and methods |
| CN106018736B (en) * | 2016-05-10 | 2018-02-02 | 北京工业大学 | Urban Underground pipeline seepage triggers the experimental rig of surface collapse |
| GB201711412D0 (en) * | 2016-12-30 | 2017-08-30 | Maxu Tech Inc | Early entry |
| CN109000158B (en) * | 2018-09-04 | 2020-05-08 | 深圳市市政设计研究院有限公司 | Pipeline leakage early warning system and method based on Internet of things and BIM and application thereof |
-
2022
- 2022-07-26 CN CN202210886058.5A patent/CN115146513B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114218651A (en) * | 2021-12-22 | 2022-03-22 | 中国地质大学（武汉） | Method and equipment for simulating particle discrete element in sandy soil foundation slab load test |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115146513A (en) | 2022-10-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| WO2020103838A1 (en) | Method and system for determining foundation bearing capacity during high-speed hydraulic tamping in real time | |
| CN113569312B (en) | Method for evaluating and pre-controlling safety margin of emptying road under urban underground construction scene | |
| CN106442937B (en) | A kind of novel sea shallow-layer soil strength variation detection system and its appraisal procedure | |
| CN105160138B (en) | Tunnel flow consolidates analysis of coupled system's stability method | |
| CN104568493A (en) | Method for rapidly identifying structure damage on basis of displacement time-history area under vehicle load | |
| CN104989456A (en) | Large-span underground construction surrounding rock excavation stability monitoring and early warning method | |
| CN114814165B (en) | Sandy seabed shear damage assessment method and assessment device under action of solitary waves | |
| CN102865952A (en) | Nondestructive testing method for working stress of concrete | |
| CN110220978A (en) | A kind of longitudinal wave guide monitoring method for tailings warehouse dam body dam break early warning | |
| CN103938660A (en) | Method for determining pile foundation bearing capacity after pile slipping | |
| CN115146513B (en) | Simulation early warning method and system for pipeline leakage collapse | |
| Wang et al. | A method for quantitative damage identification in a high-piled wharf based on modal strain energy residual variability | |
| US9995643B2 (en) | Detection of static tip resistance of a pile | |
| CN111859736B (en) | Service performance evaluation method under complex environment with multiple diseases of pipeline | |
| Abo-Youssef et al. | Numerical modelling of passive loaded pile group in multilayered soil | |
| CN111323192B (en) | Deep-buried pipeline vibration attenuation effect testing method based on dynamic compaction reinforcement | |
| CN105586974A (en) | Design method for foundation pit retaining structure | |
| CN109024717A (en) | SMW engineering method pile-type steel method for measuring stress | |
| CN111412885A (en) | Large deformation prediction method for extruded surrounding rock of large buried depth tunnel | |
| CN111428297A (en) | BP neural network-based pile foundation P-S curve determination method | |
| CN206177937U (en) | Ocean shallow soil bulk property detection system | |
| CN115455759A (en) | Dynamic monitoring and early warning method for mine side slope with weak interlayer | |
| CN103835278B (en) | Buried pipe is subject to the safety monitoring appraisal procedure that contiguous Construction Vibration affects under water | |
| Tao et al. | An anthropomorphic fuzzy model for the time-spatial assessment of sandstone seepage damage | |
| CN117074180B (en) | Method for measuring pressure change of underground soil body of building |
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 |