CN115238553B - Method and system for dividing dangerous area of buried pipeline leakage erosion - Google Patents
Method and system for dividing dangerous area of buried pipeline leakage erosion Download PDFInfo
- Publication number
- CN115238553B CN115238553B CN202210893960.XA CN202210893960A CN115238553B CN 115238553 B CN115238553 B CN 115238553B CN 202210893960 A CN202210893960 A CN 202210893960A CN 115238553 B CN115238553 B CN 115238553B
- Authority
- CN
- China
- Prior art keywords
- particle
- unit
- area
- simulation model
- leakage
- 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- 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/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
Abstract
The invention relates to a method and a system for dividing a dangerous area of buried pipeline leakage erosion, wherein the method comprises the following steps: establishing a geological generalization model; converting the generalized model into a simulation model; monitoring each unit grid in the simulation model and tracking the stress state of each unit grid; when the stress state meets the maximum elastic limit condition, deleting the unit grid, and simultaneously creating discrete particle units; continuously repeating the deleting and the creating so that the discrete unit particles are contacted with each other to form a particle cluster; monitoring a simulation model comprising the particle clusters, and analyzing a potential instability process; and partitioning the dangerous area of the leakage erosion of the buried pipeline according to the influence range of the simulation model after collapse to obtain areas with different dangerous degrees. The method makes up the one-sidedness of the structure deduction of the traditional single model, improves the reliability of the disaster-recovery development structural characteristics of the leakage pipeline, and has the advantages of simple operation, strong engineering applicability and high calculation precision.
Description
Technical Field
The invention belongs to the technical field of road collapse monitoring, and particularly relates to a method and a system for dividing a dangerous area of buried pipeline leakage erosion.
Background
The accessible road construction is an important junction for realizing convenience of urban resident life, external transportation and movement, and road bed water and soil loss caused by underground pipeline leakage easily causes road settlement and collapse, and seriously affects normal service of roads and safe trip of residents. By means of the simulation model method of the buried pipeline leakage erosion collapse, the relation between pipeline leakage and the dynamic development of the underground hidden danger of the road is deduced, the leakage collapse mechanism of the buried pipeline is disclosed, and the method is one of important means for solving disaster accidents of the road caused by the leaked pipeline.
The calculation of the road collapse simulation model at the present stage is mainly carried out by the existing finite element, finite difference and discrete element software, and the physical and mechanical characteristics of the model in the leakage collapse evolution process are analyzed. In the buried pipeline leakage erosion collapse model, leakage water is a continuous unit structure, and a soil body structure is divided into two stages, wherein the soil body structure is a continuous unit structure when not subjected to erosion action, and is converted into a discrete unit structure under the action of buoyancy, gravity and erosion force when subjected to leakage water erosion. The numerical analysis models adopted in the current stage are single models, the conversion of unit structures in the calculation process cannot be realized by adopting the coupling effect of the two models, so that the analysis rule of the model in the whole leakage erosion process in the current stage has certain one-sidedness for solving the collapse mechanism of the underground leakage pipeline of the road, and the misjudgment of the safety state of the road is easily caused due to the adaptability problem of the model structure, and safety accidents are caused.
Aiming at the defects of the existing stage model, a buried pipeline leakage erosion finite element conversion discrete particle collapse model method needs to be developed urgently, the characteristics of the whole process of buried pipeline leakage collapse are reflected by the aid of the model method, and the danger subareas of the road are fed back reasonably, objectively and accurately. Therefore, the simulation model is established through the limited grids, the unit grids are deleted and the particle clusters are synchronously created under the maximum elastic limit condition, the potential instability process analysis of the model structure is realized, the leakage erosion danger subarea of the buried pipeline is further determined, and the method has important significance for comprehensively mastering the whole process rule of the leakage-induced collapse evolution of the buried pipeline and ensuring the urban road updating speed and the safe service.
Disclosure of Invention
In order to overcome the problems in the prior art, the invention provides a method and a system for dividing a dangerous area of underground pipeline leakage erosion, which are used for solving the problems in the prior art.
A method for dividing a dangerous area of buried pipeline leakage erosion comprises the following steps:
s1, establishing a geological generalized model of leakage erosion of the buried pipeline;
s2, converting the geological generalized model into a simulation model of leakage erosion of the buried pipeline comprising a plurality of unit grids;
s3, performing mechanical calculation on each unit grid, and monitoring and tracking the stress state of each unit grid;
s4, when the current stress of the unit grid is larger than or equal to the maximum tensile stress, deleting the unit grid, and creating discrete unit particles at the deleted unit grid;
s5, continuously and repeatedly deleting the unit grids and creating discrete unit particles to enable the discrete unit particles to be mutually contacted to form a particle cluster;
s6, setting a measuring point in a simulation model comprising the particle cluster and the unit grid, and recording the stress and displacement increment change at the measuring point;
s7, judging whether a collapse structure of the simulation model is formed or not according to the stress and displacement increment states of the measuring points;
and S8, partitioning the dangerous area formed by collapse according to the collapse position of the simulation model, the cavity position and the influence zone position to obtain areas with different danger degrees.
The above-mentioned aspects and any possible implementation manners further provide an implementation manner, and the S1 includes determining physical and mechanical parameters of a stratigraphic structure of the leakage road, pipeline type parameters and a leakage form according to the working condition characteristics of the buried pipeline leakage road, and establishing a geological generalized model of the buried pipeline leakage erosion according to the parameters.
The above-mentioned aspects and any possible implementation manners further provide an implementation manner, and the S2 includes establishing the simulation model of the underground pipeline leakage erosion by using a finite element method, and obtaining a plurality of unit grids by performing discretization by using finite grids.
The above aspects, and any possible implementations, further provide an implementation,
the current stress of the unit grid satisfies the following relational expression:
D≥F+(G-F L )cosδ-F z
wherein F represents the erosion force of the leakage water to the unit grid, G represents the gravity of the unit grid, and F L Showing the buoyancy of the cell grid, F z The motion resistance of the unit grids is represented, D represents the current stress of the unit grids, delta represents the included angle between the hydraulic etching direction and the horizontal plane, and the right side of the inequality integrally represents the maximum tensile stress.
The above aspects and any possible implementations further provide an implementation in which the particle clusters have a number and are randomly distributed, wherein:
(1) The base radius of each unit particle is R, and is determined by the following relation:
the method comprises the following steps that S represents the area of a certain unit grid, N represents the number of the unit grids, N is a positive integer, 2D represents a two-dimensional plane simulation model, and 3D represents a three-dimensional simulation model;
(2) Generating boundary conditions of the particle cluster, the distance L between the sphere center of the ith particle and the jth boundary of the adjacent unit grid ij And the distance D between the center of the ith particle and the center of the kth particle ik And satisfies the following conditions:
wherein L is ij Denotes the distance, D, between the center of the ith particle and the jth boundary of the adjacent cell grid ik Denotes the distance between the spherical center of the ith particle and the spherical center of the kth particle, delta denotes the space margin, i, j, k are positive integers greater than 0, R i ,R k Respectively representing the radius of the ith particle and the kth particle;
(3) And after the ith particle is created, completely assigning the mass, the speed, the contact rigidity, the contact stress and the material parameters of the deleted unit grid to the particle.
The above-described aspect and any possible implementation further provide an implementation, where the displacement increment is related by:
wherein H u+1 Represents the amount of vertical displacement of the u +1 th measurement point, H u Represents the vertical displacement amount, P, of the u-th measurement point u+1 Represents the horizontal distance, P, from the leak point of the u +1 th measuring point u And the horizontal distance between the u-th measuring point and the leakage point is shown, and u is a positive integer greater than 0.
The above-described aspects and any possible implementations further provide an implementation in which the risk zones include a high risk zone, a medium risk zone, a low risk zone, and a no risk zone.
According to the above aspects and any possible implementation manner, an implementation manner is further provided, under the action of leakage erosion of the buried pipeline, the high-risk area is a collapse hole area formed below a road surface after the instability of a road structure, the area is easy to collapse, and high potential safety hazards exist;
the middle risk area is a cavity area formed below the road surface after the road structure is unstable, the area is easy to collapse, and higher potential safety hazards exist;
the low-risk area is a collapse cavity formed below the road surface after the road structure is unstable, and the area has the possibility of collapse and has certain potential safety hazard;
the risk-free area is an area without influence after the instability of the road structure, the road bearing capacity of the area meets the requirement, and the safety is high.
The above-described aspects and any possible implementation further provide an implementation that indicates the formation of a collapsed structure of the simulation model when the incremental states of stress and displacement of the measuring point are both zero.
The invention also provides a system for dividing the dangerous area of the buried pipeline leakage erosion, which comprises:
the model building module is used for building a geological generalized model of the leakage erosion of the buried pipeline;
the model conversion module is used for converting the geological generalized model into a simulation model of the leakage erosion of the buried pipeline comprising a plurality of unit grids;
the calculation and monitoring module is used for performing mechanical calculation on each unit grid, monitoring and tracking the stress state of the soil body structure in each unit grid;
the particle creating module is used for deleting the unit grid when the current stress of the unit grid is greater than or equal to the maximum tensile stress value, and creating discrete unit particles at the deleted unit grid;
the particle cluster creating module continuously and repeatedly deletes the unit grids and creates the discrete unit particles so that the discrete unit particles are mutually contacted to form a particle cluster;
the device comprises a determining module, a calculating module and a calculating module, wherein the determining module is used for establishing measuring points in a simulation model comprising particle clusters and unit grids and recording the stress and displacement increment changes at the measuring points;
the judging module is used for judging whether the collapse of the simulation model is formed or not according to the stress and displacement increment states of the measuring points;
and the dangerous area dividing module is used for dividing the dangerous area formed by collapse according to the collapse position of the simulation model, the cavity position and the position of the influence zone to obtain areas with different dangerous degrees.
The invention has the advantages of
Compared with the prior art, the invention has the following beneficial effects:
due to the adoption of the technical scheme, compared with the prior art, the invention has the following advantages:
(1) According to the method, a simulation model is established through a limited grid, under the condition of the maximum elastic limit, unit grids are deleted and discrete unit particles are synchronously created, deformation occurs under the action of water leakage of pipelines, displacement increment and stress change occur in measuring points arranged in the simulation model, the road collapses until the displacement increment and the stress below a road surface structure are zero at the same time, and after the road collapses, dangerous areas are divided, so that the process monitoring of the instability development of a model structure is realized, the one-sidedness of structure deduction of a traditional single model is compensated, and the reliability of the structural characteristics of the leakage pipeline disaster-pregnancy development is improved;
(2) According to the method for creating the discrete particle model by deleting the unit grids of the buried pipeline leakage erosion simulation model, discontinuous deformation and stress data of measuring points in the direction perpendicular to the measuring line are analyzed through the mutual contact action of particle clusters, the discontinuous deformation is the characteristic of a particle structure, particles are discontinuously deformed under the action of seepage, and the stress data is data monitored by the measuring points and can be monitored through the set measuring points, so that the discontinuous deformation monitoring of the continuous simulation model is realized, and the method has the advantages of simplicity in operation, strong engineering applicability and high calculation precision;
(3) According to the method, the simulation model is adopted to deduce the influence range of the road after collapse, the division of the dangerous area of the buried pipeline leakage erosion is realized, the interference of subjective factors is avoided, and the precision and the reliability of dangerous subareas are improved.
Drawings
FIG. 1 is a flow chart of a method implementation of the present invention;
FIG. 2 is a schematic diagram of a buried pipeline seepage erosion model of the present invention;
FIG. 3 is a schematic view of a simulation model of buried pipeline leakage erosion according to the present invention;
FIG. 4 is a schematic flow chart of the invention for deleting a mesh of cells to create discrete cell particles;
FIG. 5 is a schematic illustration of a particle cluster of the present invention;
FIG. 6 is a schematic view of an incremental displacement survey line arrangement of particle clusters in accordance with the present invention;
FIG. 7 is a front view of a section of the present invention at risk of buried pipeline seepage erosion;
fig. 8 is a plan view of the region of the risk of underground piping leakage etching according to the present invention.
In the figure: 1-pavement, 2-leakage pipeline, 3-roadbed stratum, 4-soil structure, 5-unit grid and 6-particle cluster.
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. 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 method for dividing a dangerous area of buried pipeline leakage erosion comprises the following steps:
A. determining road surface, stratum structure and pipeline parameters based on real working conditions, establishing a geological generalized model of the buried pipeline leakage erosion, determining road structure parameters of the buried pipeline leakage erosion according to the actual road stratum structure, pipeline distribution and road size, and establishing a geological generalized model of the buried pipeline leakage erosion according to the determined road structure parameters; the established geological generalized model is a structural model consistent with the road parameters of the real working conditions, the real road parameters are provided for establishing the simulation model, and the improved geological generalized model parameters are beneficial to improving the accuracy and the construction efficiency of the subsequent simulation model;
B. establishing a simulation model of the buried pipeline leakage erosion through finite elements based on the pavement structure, the stratum structure, the pipeline parameters, the boundary conditions and the like of a geological generalized model of the buried pipeline leakage erosion, wherein the simulation model comprises a plurality of unit grids, and the step realizes the accurate and rapid establishment of the simulation model of the buried pipeline leakage erosion;
C. performing mechanical calculation on each unit grid, monitoring and tracking the stress state of each unit grid, realizing self-monitoring of stress conditions in the simulation model, and accurately monitoring any grid unit in the simulation model;
D. when the current stress is greater than or equal to the comprehensive stress value, deleting the unit grid where the current stress is located, and simultaneously creating a particle cluster at the deleted unit grid to realize the creation of discrete unit particles and the construction of an efficient substitute unit structure;
E. continuously and repeatedly deleting the unit grid and creating discrete unit particles, so that the discrete unit particles are mutually contacted to form a particle cluster;
F. setting a vertical measuring line for a simulation model comprising particle clusters and unit grids, setting a measuring point at an interval of s distance in the direction of the measuring line, and recording displacement and stress data of each measuring point in the direction of the measuring line; judging the collapse state of the simulation model of the leakage pipeline erosion according to the displacement increment of the measuring point and the change of the stress data under the road surface;
G. and analyzing the displacement increment and the stress data, determining the process that the simulation model deforms due to stress until the simulation model collapses under the condition of buried pipeline leakage erosion, partitioning the dangerous area formed by collapsing to obtain areas with different dangerous degrees, wherein the areas with different dangerous degrees correspond to the areas with different dangerous degrees in the actual road collapse working condition, and thus realizing accurate judgment of the dangerous area of the actual road collapse working condition.
Preferably, as shown in fig. 2, when a geological modeling model of the buried pipeline leakage erosion is established, firstly, according to the characteristic of the working condition of a typical buried pipeline leakage road, then, the physical and mechanical parameters of the stratum structure, the pipeline type parameters and the leakage form are determined, and finally, the geological modeling model of the buried pipeline leakage erosion is established based on the above situation.
Preferably, as shown in fig. 3, the method for establishing a simulation model by using a finite grid in the present invention is to establish a simulation model of buried pipeline leakage erosion by a wire element method; the finite element mesh is discretized by using a thicker finite element mesh.
Preferably, the stress state of the soil structure is monitored and tracked, the stress characteristics of the soil structure are monitored through the serial number ID of each unit grid of the simulation model, each unit grid in the simulation model has one serial number ID, and the positions of the unit grids can be quickly found through the IDs, so that the current stress data of the unit grids can be directionally monitored through the serial number IDs of the unit grids, and particle units can be further created;
preferably, as shown in fig. 4 and 5, firstly, when the soil structure of the simulation model of the present invention is eroded and the unit grid reaches a value greater than or equal to the maximum tensile stress value; secondly, deleting the unit grids reaching the maximum tensile stress condition; finally, the unit grid is converted into the interior discrete unit particles.
Preferably, in the invention, when the unit grid meets the maximum tensile stress condition, the unit grid is deleted and the discrete particles are created, otherwise, the mechanical calculation is continued.
The simulation model is built according to the real soil structure and has the same parameters and attributes with the soil under the actual working condition, so that the unit grids are real in the simulation model, the road structure reflecting the real working condition is built through the model, the stress, the deformation and the like of the road structure are calculated, and the disaster recovery characteristics of the real road are deduced.
Preferably, the creating of the randomly distributed discrete unit particles in the present invention is to generate a certain number of randomly distributed particle clusters in the deleted unit grid region, and the specific method is as follows:
(1) Determining the base radius R of the particles, obtained by:
wherein S represents the area of each unit grid, N represents the number of the unit grids, N is more than or equal to 1,2D represents a two-dimensional plane simulation model, 3D represents a three-dimensional simulation model, each particle is different in size, the radius of each particle is not completely equal to the basic radius, and the radius of each particle is 0.3-1.5 times of the basic radius;
(2) Generating boundary conditions of the particle cluster, the distance L between the sphere center of the ith particle and the jth boundary of the adjacent unit grid ij And the distance D between the center of the ith particle and the center of the kth particle ik And satisfies the following conditions:
wherein L is ij Denotes the distance, D, between the center of the ith particle and the jth boundary of the adjoining cell grid ik Denotes the distance between the sphere center of the ith particle and the sphere center of the kth particle, delta denotes the space margin, i, j, k are positive integers greater than 0, R i ,R k Respectively representing the radius of the ith particle and the kth particle;
(3) After the ith particle is created, completely assigning the mass, the speed, the contact rigidity, the contact stress and the material parameters of the deleted unit grid to the particle;
(4) And repeating the steps, continuously generating discrete unit particles in the area of the deleted unit grid, and enabling the particles to be mutually contacted under the action of pipeline seepage flow to form a particle cluster.
Preferably, as shown in fig. 5, continuous unit grids in the simulation model are converted into discrete unit particles, the discrete unit particles interact to form particle clusters, so that the structure of the simulation model can better reflect the real working conditions, the displacement increment is the displacement increment of the measuring points arranged on the measuring line, the measuring points do not change due to any unit change and are fixed points in the simulation model, the displacement increment and the stress change of the measuring points in the simulation model are monitored, whether the simulation model forms collapse or not is reflected through the displacement increment and the stress value of the measuring points, specifically, for the simulation model, the measuring line is made through a vertical road surface, one monitoring point is arranged at certain intervals along the measuring line direction, the displacement and the stress data of each measuring point in the measuring line direction are recorded, and data support is provided for judging the forming process of collapse by monitoring the stress and the displacement data of the measuring points. The displacement increment is expressed as follows:
wherein H u+1 Represents the vertical displacement of the u +1 measuring point of the simulation model, H u Represents the vertical displacement, P, of the u-th measuring point of the simulation model u+1 Represents the horizontal distance between the u +1 measuring point and the leakage point of the simulation model, P u And the horizontal distance between the u-th measuring point of the model and the leakage point is represented, and u is a positive integer greater than 0.
According to the invention, under the hydraulic action of the leakage pipeline, after the simulation model is converted into the particle cluster after the unit grid is deleted, the area range is enlarged, so that the stratum overlying the leakage pipeline collapses under the self-weight action, the stratum is easy to deform, and whether the simulation model collapses or not is judged by taking the displacement increment and the stress below the road surface as zero reflection through the displacement increment and the stress data of the measuring points in the measuring line direction. The displacement increment is 0, which indicates that the calculation of the simulation model has reached balance, the calculation reaches the precision requirement, and the whole simulation process is stopped; and when the stress below the road surface is 0, indicating that the roadbed soil in the simulation model is eroded to form a cavity below the road surface, radiating the cavity to the periphery by taking the stress below the road surface as a middle point according to different boundary conditions of the road structure after collapse, and dividing dangerous areas by taking the collapsed area as a high-risk area, the area with the cavity as a middle-risk area, the area with influence as a low-risk area and the area without influence as a no-risk area.
Preferably, as shown in fig. 6, 7 and 8, the invention calculates according to the model value, and the calculation is that the simulation model obtains the data result of the required stress and displacement through the calculation of the model structure, and the influence range of the collapse is presented. The collapse influence range formed after the method is finished is characterized in that a circle with the diameter of CD is a collapsed hole area, an annular area formed by BC and DE is a cavity area, an annular area formed by AB and EF is an etching collapse influence area, an area formed by OA and FG is a compact stratum structure, the leakage of the buried pipeline is subjected to dangerous zoning, the specific zoning is that the collapse result formed after the calculation of numerical values is finished, the stress below the road surface is taken as a middle point, the area with collapse is taken as a high risk area, the area with cavities is taken as an intermediate risk area, the area with influence is taken as a low risk area, and the area without influence is taken as a no risk area, and the structure formed after the calculation of a simulation model is finished is shown in figure 7. Wherein the zones are illustrated below:
high risk zone: the method is characterized in that a collapse opening area is formed below a road surface after the road structure is unstable under the action of leakage erosion of an underground pipeline, and the area is a circular area shown by the diameter CD and is easy to collapse under the action of slight road surface load, so that high potential safety hazards exist;
area of stroke risk: the method is characterized in that under the action of leakage erosion of the buried pipeline, a cavity region is formed below the road surface after the road structure is unstable, for example, an annular region formed by BC and DE is easy to collapse under the action of common road surface load, and higher potential safety hazard exists;
low risk zone: the method comprises the steps that a collapse cavity is formed below a road surface after a road structure is unstable under the action of buried pipeline leakage erosion, a perpendicular line is formed along a boundary line inflection point of the cavity in an area below the leakage pipeline and perpendicular to a line connecting the leakage point, and an erosion collapse influence area, such as an annular area formed by AB and EF, is determined, and has the possibility of collapse under the action of high-strength road surface load, so that certain potential safety hazards exist;
no risk zone: the method shows that under the action of leakage erosion of the buried pipeline, no influence area is generated after the road structure is unstable, such as areas formed by OA and FG, the road bearing capacity of the area meets the requirement of a compact road structure, and the safety is high.
According to the method, the range of the road collapse influence zone is determined according to the collapse result of the leakage pipeline simulation model, different danger levels are divided, a basis is provided for judging the range of the actual leakage influence zone of the buried pipeline, the danger zone of the buried pipeline in the actual road condition can be judged quickly and accurately, and an important basis is provided for the safe and timely disposal of the road.
Preferably, the present invention also provides a system for dividing a hazardous area of buried pipeline leakage erosion, comprising:
the model building module is used for building a geological generalized model of the leakage erosion of the buried pipeline;
the model conversion module is used for converting the geological generalized model into a simulation model of the leakage erosion of the underground pipeline comprising a plurality of unit grids;
the calculation and monitoring module is used for performing mechanical calculation on each unit grid, monitoring and tracking the stress state of the soil body structure in each unit grid;
the particle creating module deletes one unit grid when the stress state of the soil body structure of the unit grid meets the maximum elastic limit condition, and creates discrete unit particles at the deleted unit grid;
the particle cluster creating module continuously and repeatedly deletes the unit grids and creates the discrete unit particles so that the discrete unit particles are mutually contacted to form a particle cluster;
the determining module is used for acquiring displacement increment and stress of the simulation model including the particle clusters under the erosion of the leakage pipeline and determining the destabilization development process of the simulation model;
and the dangerous area dividing module is used for analyzing the displacement increment and the stress data, determining the process that the simulation model collapses under the condition of leakage erosion of the buried pipeline, and dividing the dangerous area formed by collapse to obtain areas with different dangerous degrees.
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 (9)
1. A method for dividing a dangerous area of buried pipeline leakage erosion is characterized by comprising the following steps:
s1, establishing a geological generalized model of the leakage erosion of the buried pipeline;
s2, converting the geological generalized model into a simulation model of leakage etching of the underground pipeline comprising a plurality of unit grids;
s3, performing mechanical calculation on each unit grid, and monitoring and tracking the stress state of each unit grid;
s4, when the current stress of the unit grid is larger than or equal to the maximum tensile stress, deleting the unit grid, and creating discrete unit particles at the deleted unit grid;
s5, continuously and repeatedly deleting the unit grids and creating discrete unit particles to enable the discrete unit particles to be mutually contacted to form a particle cluster;
s6, setting measuring points in a simulation model comprising the particle clusters and the unit grids, and recording the stress and displacement increment changes at the measuring points;
s7, judging whether a collapse structure of the simulation model is formed or not according to the stress and displacement increment states of the measuring points;
s8, partitioning a dangerous area formed by collapse according to the collapse position, the cavity position and the influence zone position of the simulation model to obtain areas with different danger degrees;
wherein the particle clusters have a number and are randomly distributed, wherein:
(1) The base radius of each unit particle is R, which is determined by the following relation:
the method comprises the following steps that S represents the area of a certain unit grid, N represents the number of the unit grids, N is a positive integer, 2D represents a two-dimensional plane simulation model, and 3D represents a three-dimensional simulation model;
(2) Generating a boundary condition of a particle cluster, a distance L between a sphere center of an ith particle and a jth boundary of an adjacent unit grid ij And the distance D between the center of the ith particle and the center of the kth particle ik And satisfies the following conditions:
wherein L is ij Denotes the distance, D, between the center of the ith particle and the jth boundary of the adjoining cell grid ik Denotes the distance between the sphere center of the ith particle and the sphere center of the kth particle, delta denotes the space margin, i, j, k are positive integers greater than 0, R i ,R k Respectively representing the radius of the ith particle and the kth particle;
(3) And after the ith particle is created, completely assigning the mass, the speed, the contact rigidity, the contact stress and the material parameters of the deleted unit grid to the particle.
2. The method for dividing the dangerous area of the buried pipeline leakage erosion according to claim 1, wherein the step S1 comprises determining the physical and mechanical parameters of the stratum structure of the leakage road, the pipeline type parameters and the leakage form according to the working condition characteristics of the buried pipeline leakage road, and establishing a geological modeling of the buried pipeline leakage erosion according to the parameters.
3. The method for dividing the dangerous area of the underground pipeline leakage erosion according to claim 2, wherein the step S2 comprises establishing the underground pipeline leakage erosion simulation model by using a finite element method, and obtaining a plurality of unit grids by using finite grids for dispersion.
4. The method for dividing a dangerous area for underground pipeline leakage erosion according to claim 1,
the current stress of the unit grid satisfies the following relation:
D≥F+(G-F L )cosδ-F z
wherein F represents the erosion force of the leakage water to the unit grid, G represents the gravity of the unit grid, and F L Showing the buoyancy of the cell grid, F z The motion resistance of the unit grid is represented, D represents the current stress of the unit grid, delta represents the included angle between the hydraulic etching direction and the horizontal plane, and the right integer of the inequalityThe maximum tensile stress is indicated.
5. The method of dividing a hazardous area for buried pipeline leak erosion according to claim 1, wherein said displacement increment is expressed by the following relation:
wherein H u+1 Represents the amount of vertical displacement of the u +1 th measurement point, H u Represents the vertical displacement amount, P, of the u-th measurement point u+1 Represents the horizontal distance, P, from the leak point of the u +1 th measuring point u And the horizontal distance between the u-th measuring point and the leakage point is shown, and u is a positive integer greater than 0.
6. The method for dividing a dangerous area for underground pipeline leakage attack according to claim 1, wherein the dangerous area comprises a high-risk area, a medium-risk area, a low-risk area and a no-risk area.
7. The method for dividing the dangerous area of the buried pipeline leakage erosion according to claim 6, wherein under the action of the buried pipeline leakage erosion, the high-risk area is a cave-in opening area formed below a road surface after the instability of a road structure, the area is easy to collapse, and high potential safety hazards exist;
the middle danger zone is a cavity zone formed below the road surface after the road structure is unstable, the zone is easy to collapse, and higher potential safety hazard exists;
the low-risk area is a collapse cavity formed below the road surface after the road structure is unstable, and the area has the possibility of collapse and has certain potential safety hazard;
the risk-free area is an area without influence after the instability of the road structure, the road bearing capacity of the area meets the requirement, and the safety is high.
8. The method for dividing the dangerous area for the leakage erosion of the buried pipeline according to claim 1, wherein when the stress and displacement increment states of the measuring points are both zero, the collapse structure of the simulation model is formed.
9. A system for dividing a hazardous area for underground pipeline leakage erosion, comprising:
the model building module is used for building a geological generalized model of the leakage erosion of the buried pipeline;
the model conversion module is used for converting the geological generalized model into a simulation model of the leakage erosion of the underground pipeline comprising a plurality of unit grids;
the calculation and monitoring module is used for performing mechanical calculation on each unit grid, monitoring and tracking the stress state of the soil body structure in each unit grid;
the particle creating module is used for deleting the unit grid when the current stress of the unit grid is greater than or equal to the maximum tensile stress value, and creating discrete unit particles at the deleted unit grid;
the particle cluster creating module continuously and repeatedly deletes the unit grids and creates the discrete unit particles so that the discrete unit particles are mutually contacted to form a particle cluster;
the device comprises a determining module, a calculating module and a calculating module, wherein the determining module is used for establishing measuring points in a simulation model comprising particle clusters and unit grids and recording the stress and displacement increment changes at the measuring points;
the judging module is used for judging whether the collapse of the simulation model is formed or not according to the stress and displacement increment states of the measuring points;
the danger area dividing module is used for dividing the danger area formed by collapse according to the collapse position of the simulation model, the position of the cavity and the position of the influence zone to obtain areas with different danger degrees,
wherein the particle clusters have a number and are randomly distributed, wherein:
(1) The base radius of each unit particle is R, and is determined by the following relation:
the method comprises the following steps that S represents the area of a certain unit grid, N represents the number of the unit grids, N is a positive integer, 2D represents a two-dimensional plane simulation model, and 3D represents a three-dimensional simulation model;
(2) Generating boundary conditions of the particle cluster, the distance L between the sphere center of the ith particle and the jth boundary of the adjacent unit grid ij And the distance D between the center of the ith particle and the center of the kth particle ik And satisfies the following conditions:
wherein L is ij Denotes the distance, D, between the center of the ith particle and the jth boundary of the adjacent cell grid ik Denotes the distance between the sphere center of the ith particle and the sphere center of the kth particle, delta denotes the space margin, i, j, k are positive integers greater than 0, R i ,R k Respectively representing the radius of the ith particle and the kth particle;
(3) After the ith particle is created, the mass, the speed, the contact rigidity, the contact stress and the material parameters of the deleted unit grid are completely assigned to the particle.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210893960.XA CN115238553B (en) | 2022-07-27 | 2022-07-27 | Method and system for dividing dangerous area of buried pipeline leakage erosion |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210893960.XA CN115238553B (en) | 2022-07-27 | 2022-07-27 | Method and system for dividing dangerous area of buried pipeline leakage erosion |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115238553A CN115238553A (en) | 2022-10-25 |
| CN115238553B true CN115238553B (en) | 2023-03-10 |
Family
ID=83677274
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210893960.XA Active CN115238553B (en) | 2022-07-27 | 2022-07-27 | Method and system for dividing dangerous area of buried pipeline leakage erosion |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115238553B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116341295B (en) * | 2023-05-30 | 2023-07-25 | 交通运输部公路科学研究所 | Road stability assessment method and system based on underground pipeline |
| CN117236706B (en) * | 2023-11-16 | 2024-02-20 | 海纳云物联科技有限公司 | Evaluation method, device, equipment and storage medium for leakage risk of heating power pipeline |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110487697A (en) * | 2019-07-29 | 2019-11-22 | 北京科技大学 | Infuse supercritical carbon dioxide coal petrography mechanical property testing and fracturing experiments device |
| CN111122348A (en) * | 2019-11-25 | 2020-05-08 | 温州大学 | Computing method for predicting rheological strength of ultra-soft soil with sea area |
| CN111862323A (en) * | 2020-07-08 | 2020-10-30 | 北京科技大学 | Multi-element pregnant disaster digital twin intelligent perception identification early warning system and method |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8081796B2 (en) * | 2008-11-24 | 2011-12-20 | Ingrain, Inc. | Method for determining properties of fractured rock formations using computer tomograpic images thereof |
| CN107978138A (en) * | 2017-12-05 | 2018-05-01 | 南瑞集团有限公司 | A kind of disaster monitoring method for early warning based on mountain torrents dynamical evolution simulation model |
| CN114492140A (en) * | 2022-02-18 | 2022-05-13 | 中国科学院武汉岩土力学研究所 | Method for dividing dangerous partitions in excavation of high and steep joint slope |
-
2022
- 2022-07-27 CN CN202210893960.XA patent/CN115238553B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110487697A (en) * | 2019-07-29 | 2019-11-22 | 北京科技大学 | Infuse supercritical carbon dioxide coal petrography mechanical property testing and fracturing experiments device |
| CN111122348A (en) * | 2019-11-25 | 2020-05-08 | 温州大学 | Computing method for predicting rheological strength of ultra-soft soil with sea area |
| CN111862323A (en) * | 2020-07-08 | 2020-10-30 | 北京科技大学 | Multi-element pregnant disaster digital twin intelligent perception identification early warning system and method |
Non-Patent Citations (6)
| Title |
|---|
| Determination of microcystins in environmental water samples with ionic liquid magnetic graphene;Xiaoyan Liu等;《Ecotoxicology and Environmental Safety》;20191231;20-26 * |
| Mining fault association rules in the perception layer of electric power sensor network based on improved Eclat;Yongsheng Zhao;《IWCMC》;20211231;418-423 * |
| 基于临界应力指数法巷道冲击地压危险性评价;王爱文等;《煤炭学报》;20200515(第05期);76-84 * |
| 复杂地质构造对热采水平井弯管段套管应力影响的数值分析;吕祥锋等;《武汉工业学院学报》;20090915(第03期);63-67 * |
| 气固颗粒系统模拟的研究进展;苏军伟等;《化学反应工程与工艺》;20160625(第03期);71-86 * |
| 盾构隧道洞周土压力颗粒流数值分析;邢心魁等;《重庆交通大学学报(自然科学版)》;20151231(第05期);33-36+52 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115238553A (en) | 2022-10-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN115238553B (en) | Method and system for dividing dangerous area of buried pipeline leakage erosion | |
| Nikolić et al. | Overview of the numerical methods for the modelling of rock mechanics problems | |
| Katzenbach et al. | Assessing settlement of high‐rise structures by 3D simulations | |
| CN102368277B (en) | Building method of load-structure model in view of arch effect of tunnel stress | |
| CN114297864B (en) | Cracked loose rock mass slope stability analysis method controlled by steep and gentle dip angles | |
| CN115167212B (en) | Dynamic construction control system and method for foundation pit based on monitoring platform | |
| CN104989456A (en) | Large-span underground construction surrounding rock excavation stability monitoring and early warning method | |
| CN108842821B (en) | Calculation method for reasonable buried depth of submarine tunnel constructed by drilling and blasting method | |
| CN113700053B (en) | BIM-based construction process foundation pit deformation monitoring and early warning method and system | |
| CN111400949A (en) | Rock mass underground structure three-dimensional model construction and seismic response analysis method and system | |
| CN112199758A (en) | Foundation pit engineering design and dynamic risk analysis method and system based on BIM technology | |
| Wang et al. | Three-dimensional flood routing of a dam break based on a high-precision digital model of a dense urban area | |
| Li et al. | A meshed kinematical approach for 3D slope stability analysis | |
| Zhang et al. | Wireless monitoring–based real-time analysis and early-warning safety system for deep and large underground caverns | |
| Chen et al. | Spatial estimation of material parameters and refined finite-element analysis of rockfill dam based on construction digitization | |
| CN113221462A (en) | Foundation pit soil parameter inversion method combining PSO-DE algorithm and neural network | |
| Huang et al. | Connectivity evaluation for three‐dimensional fracture network in support‐based model: A case study in the Ordos Basin, China | |
| Ai et al. | The use of shape accel array for deformation monitoring and parameter inversion of a 300 m ultrahigh rockfill dam | |
| Vahdati et al. | Inverse hardening soil parameter identification of an earth and rockfill dam by genetic algorithm optimization | |
| CN114492140A (en) | Method for dividing dangerous partitions in excavation of high and steep joint slope | |
| CN113221399A (en) | Decision-making method for railway tunnel fabricated lining component division scheme | |
| CN113435075A (en) | Computer method for simulating load response caused by landslide disaster impact on seabed energy pipeline | |
| CN111046468A (en) | Underground comprehensive pipe gallery anti-seismic design method | |
| CN206627073U (en) | The quick monitoring and forecasting device of underground engineering construction based on 3D laser scannings | |
| Dmitriev et al. | Development and approbation of the mathematical multi-model basis of the monitoring system for the “heterogeneous filtering base–structures” systems of the responsible building objects (on the example of Zagorskaya PSPP) |
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 |