CN114925577A - Vulnerability analysis method and system for masonry building on creeping type landslide body - Google Patents
Vulnerability analysis method and system for masonry building on creeping type landslide body Download PDFInfo
- Publication number
- CN114925577A CN114925577A CN202210637742.XA CN202210637742A CN114925577A CN 114925577 A CN114925577 A CN 114925577A CN 202210637742 A CN202210637742 A CN 202210637742A CN 114925577 A CN114925577 A CN 114925577A
- Authority
- CN
- China
- Prior art keywords
- building
- landslide
- vulnerability
- real
- time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
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
- 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
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A50/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE in human health protection, e.g. against extreme weather
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
Abstract
The invention relates to a method and a system for analyzing vulnerability of a masonry building on a creeping type landslide body, wherein the method comprises the following steps: constructing a landslide and building coupling geological model according to the mapping data of the target landslide body; applying an inducing factor to the landslide and building coupled geological model; calculating the real-time inclination rate of a building and the real-time displacement of the surface of a sliding body where the building is located under a preset inducing factor by adopting a discrete unit numerical method based on continuous medium mechanics; determining a physical vulnerability curve of the building according to the real-time inclination rate of the building and the real-time displacement of the surface of a sliding body where the building is located; and evaluating the landslide risk according to the physical vulnerability curve of the building. The method can accurately evaluate the physical vulnerability of the building under the action of specific inducing factors, greatly improves the quantification and the reliability of the physical vulnerability evaluation of the masonry building on the creeping landslide body, and provides powerful basis for the risk quantitative evaluation and the accurate disaster reduction of landslide disasters.
Description
Technical Field
The invention relates to the technical field of geological disaster risk assessment, in particular to a method and a system for analyzing vulnerability of a masonry building on a creeping landslide body.
Background
Along with the current climate warm humidification aggravation, the earthquake activity is gradually intensified, the human engineering activity aggravation, the occurrence frequency and the scale of the loess area reviving landslide and the new landslide hazard are continuously increased, most of the landslide and the new landslide hazard belong to large or extra-large creeping landslides, the risk of continuous creeping deformation or local instability is high, villages and people on a slope body are mainly threatened, according to field investigation statistics, village buildings are mostly brick masonry structures, and the risk born by the village landslide hazard is very large.
The vulnerability assessment of the disaster-bearing body is a main technical bottleneck restricting the risk assessment and management of landslide disasters, but a widely accepted quantitative vulnerability assessment method of the disaster-bearing body is still lacked at present. Especially for landslide disaster research in the northwest loess areas of China, the research on the disaster mechanism and stability of the landslide disaster is paid more attention before, and the research on a method for quantitatively evaluating the vulnerability of typical landslides and building types is very deficient.
The physical vulnerability of the disaster-bearing body is generally under the action of natural conditions with specific strength, and the loss degree of the disaster-bearing body is from 0 (no damage) to 1 (complete damage). At present, the building vulnerability assessment method commonly used internationally is mainly divided into two types: a heuristic evaluation model and a statistical rule-based evaluation model. A heuristic evaluation method has an important guiding function on the analysis of vulnerability of a regional scale, but because the method usually ignores the dynamic damage action mechanism of a sliding body on a building, the evaluation precision and the applicability are poor. The evaluation method based on the statistical law obtains the vulnerability curve of the building in the disaster through nonlinear regression analysis on the basis of analyzing the disaster historical data, and the vulnerability curve has dependence and regional limitation on the disaster historical data, and the generalization capability of the model is poor.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention aims to provide a method and a system for analyzing the vulnerability of a masonry building on a creeping type landslide body.
In order to achieve the purpose, the invention provides the following scheme:
a method for analyzing vulnerability of masonry buildings on a creeping type landslide body comprises the following steps:
obtaining surveying and mapping data of a target landslide body; the surveying and mapping data comprises a slope body structure, a rock-soil body type, building distribution, a digital orthographic image and a digital elevation model;
constructing a landslide and building coupling geological model according to the mapping data of the target landslide body;
determining physical and mechanical parameters of rock and soil mass and a building masonry structure in a landslide and building coupling geological model;
applying an inducing factor into the landslide and building coupled geological model; the inducers include earthquakes and rainfall of varying intensity;
calculating the real-time inclination rate of the building and the real-time displacement of the surface of a sliding body where the building is located under a preset inducing factor by adopting a discrete unit numerical method based on continuous medium mechanics;
determining a physical vulnerability curve of the building according to the real-time inclination rate of the building and the real-time displacement of the surface of the sliding body where the building is located;
and evaluating the landslide risk according to the building physical vulnerability curve.
Preferably, the applying of the inducing factors to the landslide and building coupled geological model comprises:
acquiring a preset seismic waveform;
amplifying, low-pass filtering and integral operation processing are carried out on the preset seismic waveform to obtain a seismic velocity time-course curve;
and performing reference correction on the seismic motion speed time-course curve, and obtaining an external seismic load applied to the landslide and building coupling geological model by adopting an induction factor numerical conversion formula.
Preferably, the value conversion formula of the induction factor is as follows:
wherein σ n For normal stresses applied at viscous boundaries, σ s For tangential stresses applied at viscous boundaries, C p Is the wave velocity of P wave, C s Is the wave velocity, v, of the S wave n Is the normal velocity component on the model boundary, v s The tangential velocity component on the boundary of the model is shown, and rho is the density of the rock and soil mass of the sliding bed.
Preferably, the determining the physical vulnerability curve of the building according to the real-time inclination rate of the building and the real-time displacement of the sliding body surface where the building is located includes:
determining the damage value of the corresponding building according to the real-time tilt rate of the building;
and fitting by adopting an improved Weibull function according to the real-time displacement of the surface of the sliding body where the building is located and the data of the vulnerability value of the corresponding building to obtain a physical vulnerability curve of the building.
Preferably, the modified Weibull function is:
wherein V is physical vulnerability, D is landslide surface displacement, a represents a first fitting parameter, and b represents a second fitting parameter.
The invention also provides a system for analyzing the vulnerability of the masonry building on the creeping type landslide body, which comprises:
the surveying and mapping data acquisition module is used for acquiring surveying and mapping data of the target landslide body; the surveying and mapping data comprises a slope body structure, a rock-soil body type, building distribution, a digital orthographic image and a digital elevation model;
the geological model building module is used for building a landslide and building coupling geological model according to the surveying and mapping data of the target landslide body;
the mechanical parameter determining module is used for determining physical and mechanical parameters of rock and soil mass and building masonry structures in the landslide and building coupling geological model;
an induction factor applying module for applying an induction factor to the landslide and building coupled geological model; the inducers include earthquakes and rainfall of varying intensity;
the collaborative deformation calculation module is used for calculating the real-time inclination rate of the building and the real-time displacement of the surface of the sliding body where the building is located under the preset inducing factors by adopting a discrete unit numerical method based on continuous medium mechanics;
the vulnerability curve determining module is used for determining a physical vulnerability curve of the building according to the real-time inclination rate of the building and the real-time displacement of the surface of the sliding body where the building is located;
and the landslide risk evaluation module is used for evaluating the landslide risk according to the building physical vulnerability curve.
Preferably, the inducer applying module includes:
the induction factor acquisition unit is used for acquiring a preset seismic waveform;
the seismic waveform processing unit is used for amplifying the preset seismic waveform, performing low-pass filtering and performing integral operation processing to obtain a seismic velocity time-course curve;
and the seismic load determining unit is used for performing reference correction on the seismic velocity time-course curve and obtaining the external seismic load applied to the landslide and building coupling geological model by adopting an induction factor numerical value conversion formula.
Preferably, the value conversion formula of the induction factor is as follows:
wherein σ n For normal stresses imposed on viscous boundaries, σ s For tangential stresses applied at viscous boundaries, C p Is the wave velocity of P wave, C s Is the wave velocity, v, of the S wave n As model boundariesVelocity component of upper normal, v s The tangential velocity component on the boundary of the model is shown, and rho is the density of the rock and soil mass of the sliding bed.
Preferably, the vulnerability curve determination module includes:
the system comprises a vulnerability value determining unit, a building real-time gradient determining unit and a building real-time gradient determining unit, wherein the vulnerability value determining unit is used for determining the vulnerability value of a corresponding building according to the real-time gradient of the building;
and the fitting unit is used for fitting by adopting an improved Weibull function according to the real-time displacement of the surface of the sliding body where the building is located and the data of the vulnerability value of the corresponding building to obtain a physical vulnerability curve of the building.
Preferably, the modified Weibull function is:
wherein V is the physical vulnerability, D is the landslide surface displacement, a represents the first fitting parameter, and b represents the second fitting parameter.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
the method and the system for analyzing the vulnerability of the masonry building on the creeping type landslide body have the advantages that: compared with the prior art, the landslide and building coupling geological model is constructed by obtaining the mapping data of the target landslide body, then the induction factor is applied to the geological model, the real-time inclination rate of the building under the preset induction factor and the real-time displacement of the surface of the landslide body where the building is located are calculated, and then the physical vulnerability curve of the building is determined based on the real-time inclination rate and the real-time displacement, so that the physical vulnerability of the building can be accurately evaluated under the action of the specific induction factor, the quantification and the reliability of the physical vulnerability evaluation of the building on the creeping landslide body are greatly improved, and a powerful basis is provided for the disaster quantitative evaluation and the accurate disaster reduction of the landslide.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without inventive exercise.
FIG. 1 is a flow chart of a method for analyzing vulnerability of masonry buildings on a creeping type landslide mass according to an embodiment of the present invention;
FIG. 2 is a diagram illustrating physical and mechanical parameters of landslide rock-soil mass required for calculating dynamic action of landslide and buildings in an embodiment of the present invention;
FIG. 3 is a graph showing the calculated flexural tensile strength (f) of brickwork wall along mortar joints for landslide and building dynamics in an embodiment of the invention t ) And shear strength (f τ) values;
FIG. 4 illustrates the masonry structure building damage level and the quantization index limit in an embodiment of the present invention;
FIG. 5 is a diagram illustrating a relationship between a building tilt rate limit and a recommended vulnerability value in an embodiment of the present invention;
FIG. 6 is a coupled geological model of a landslide and masonry building in an embodiment provided by the invention;
FIG. 7 is a graph of horizontal (up) and vertical (down) seismic loads (velocity time course) to be input for the calculation of landslide building under the armed seismic condition in an embodiment provided by the present invention;
FIG. 8 is a cloud view of the A-D displacement of brick masonry buildings on a slope in an embodiment of the present invention;
FIG. 9 is a graph showing the ground surface displacement (left) of the slider where the building is located and the relationship between the ground surface displacement and the inclination rate of the building (right) in the embodiment provided by the present invention;
FIG. 10 is a graph of vulnerability of a building in a brickwork structure on a downhill slope in a armed seismic condition in an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the 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.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
Referring to fig. 1, a method for analyzing vulnerability of masonry building on creeping landslide body includes:
step 1: obtaining surveying and mapping data of a target landslide body; the surveying and mapping data comprises a slope body structure, a rock-soil body type, building distribution, a digital orthographic image and a digital elevation model;
step 2: constructing a landslide and building coupling geological model according to the mapping data of the target landslide body;
in practical application, the landslide coupling geological model can be found out by adjusting and drawing the slope body structure, the rock-soil body type and the building distribution of the landslide on site, the high-precision DOM (digital ortho-image) and DEM (digital elevation model) data of the landslide is obtained by surveying and mapping by using the unmanned aerial vehicle, and the landslide and building coupling geological model is established by selecting typical sections.
And 3, step 3: determining physical and mechanical parameters of rock and soil mass and a building masonry structure in a landslide and building coupling geological model;
furthermore, the method can obtain landslide rock and soil mass samples through drilling or pit testing, and then measure the basic physical and mechanical parameters of the rock and soil mass indoors, and mainly comprises the following steps: density, modulus of elasticity, poisson's ratio, cohesion and internal friction angle, as shown in figure 2. The values of the mechanical parameters of masonry structure buildings are mainly determined according to masonry structure design specifications (GB50003-2011), the damage of brick masonry buildings is generally damaged along the mortar joint section of masonry, so the values of the tensile strength and the shear strength are determined according to different mortar grades, and the value standard is shown in figure 3.
And 4, step 4: applying an inducing factor to the landslide and building coupled geological model; the inducers include earthquakes and rainfall of varying intensity;
the rainfall conditions are considered mainly by applying rainfall flow boundaries of different intensities (generally considering medium rain, heavy rain and extra heavy rain) on the model surface, and then carrying out saturated seepage and saturated region rock-soil body softening effect analysis. The consideration of earthquake intensity is mainly based on the Chinese earthquake peak acceleration region chart (GB 18306-2015). The method comprises the steps of firstly obtaining peak acceleration of a landslide site under three earthquake scenes of 50 years exceeding probability of 60% (frequently encountered earthquakes), 10% (fortifying earthquakes) and 2% (rarely encountered earthquakes), then selecting a typical earthquake reference waveform, carrying out filtering, amplification, correction and integration processing on earthquake waves by utilizing Seimosoft software to obtain earthquake motion speed time-course curves under the three earthquake scenes, and finally converting the earthquake motion speed time-course curves into a dynamic stress load form to be applied to a model bottom boundary.
The present invention takes the application of earthquake as an example, and explains the process of applying the inducing factors in the step 4:
acquiring a preset seismic waveform;
amplifying, low-pass filtering and integral operation processing are carried out on the preset seismic waveform to obtain a seismic velocity time-course curve;
and performing reference correction on the seismic motion speed time-course curve, and obtaining an external seismic load applied to the landslide and building coupling geological model by adopting an induction factor numerical conversion formula. Wherein, the numerical conversion formula of the induction factors is as follows:
wherein σ n For normal stresses imposed on viscous boundaries, σ s For tangential stresses applied at viscous boundaries, C p Is the wave velocity of P wave, C s Is the wave velocity, v, of the S wave n Is the normal velocity component on the model boundary, v s The tangential velocity component on the boundary of the model is shown, and rho is the density of the rock and soil mass of the sliding bed.
And 5: calculating the real-time inclination rate of a building and the real-time displacement of the surface of a sliding body where the building is located under a preset inducing factor by adopting a discrete unit numerical method based on continuous medium mechanics;
the invention adopts a continuous medium mechanics-based discrete unit method (CDEM) during calculation, effectively combines the advantages of finite elements and discrete elements, and can realize the analysis of the whole process of the starting course, the sliding and the movement of the landslide under the coupling action of internal and external power and the simulation of the dynamic destruction effect of the landslide deformation instability on the building. In the invention, the method mainly simulates the dynamic interaction of the sliding mass and the building under different rainfall and earthquake intensities, and discloses the cooperative deformation rule of the sliding mass and the brick masonry building. Furthermore, the method can calculate the real-time inclination rate L of the building by extracting the real-time displacement of the top and the bottom of the building, then extract the real-time displacement D of the surface of the sliding body where the building is located, and finally establish a correlation curve of the surface displacement of the sliding body and the inclination rate of the building by taking the surface displacement of the sliding body as an abscissa and the inclination rate of the building as an ordinate.
Step 6: determining a physical vulnerability curve of the building according to the real-time inclination rate of the building and the real-time displacement of the surface of the sliding body where the building is located;
further, step 6 comprises:
determining the damage value of the corresponding building according to the real-time tilt rate of the building;
in practical application, the invention refers to a structural damage grade dividing method commonly adopted in the seismic field, and the performance of the masonry building under the influence of landslide disasters is divided into: basic sound, mild failure, moderate failure, severe failure and failure were rated 5. Meanwhile, according to the field survey result, 5 damage levels corresponding to the limit value of the quantitative index taking the whole inclination rate of the building as vulnerability are determined (as shown in figures 4-5). Where substantially intact means that substantially no reinforcement or repair is required, slight to severe damage means that different degrees of reinforcement and repair are required, and failure means complete unavailability or that the repair reinforcement costs exceed the reconstruction costs.
And fitting by adopting an improved Weibull function according to the real-time displacement of the surface of the sliding body where the building is located and the data of the vulnerability of the corresponding building to obtain a physical vulnerability curve of the building. Wherein the improved Weibull function is:
wherein V is the physical vulnerability, D is the landslide surface displacement, a represents the first fitting parameter, and b represents the second fitting parameter.
Furthermore, according to the method, more than 10 groups of data of landslide surface displacement D and vulnerability V can be selected according to the correlation curve of the surface displacement of the sliding body and the inclination rate of the building and the corresponding relation between the limit value of the inclination rate of the building and the vulnerability recommended value, and finally, the improved Weibull function is used for establishing the physical vulnerability curve of the masonry structure.
And 7: and evaluating the landslide risk according to the physical vulnerability curve of the building.
The method comprises the steps of establishing a landslide and building coupled geological analysis model on the basis of a creeping landslide body structure, a disaster mode and building distribution, analyzing dynamic interaction and cooperative deformation rules of a creeping landslide body and a masonry building by using a discrete unit method (CDEM) based on continuous medium mechanics through inputting rainfall and earthquake motion inducing conditions with different probabilities, establishing a masonry building inclination rate as a vulnerability assessment index, and establishing a vulnerability curve of the masonry building based on an improved Weibull function by combining a performance level and a damage level of the masonry building, so that the physical vulnerability of the masonry building can be accurately assessed according to creeping landslide surface displacement under the action of specific inducing factors and strength, and a powerful basis is provided for quantitative assessment and accurate disaster reduction of landslide hazard risks. The method greatly improves the quantification and the reliability of physical vulnerability assessment of the masonry building on the creeping type landslide body, and has important significance for realizing the landslide risk quantitative assessment.
The method for analyzing the vulnerability of the masonry building on the creeping landslide body is further explained by taking a loess shale contact surface landslide as an example:
(1) landslide and building coupling geological model: according to the embodiment of the invention, the landslide of a contact surface of a loess mudstone is taken as a research object, and a vulnerability assessment curve of a masonry structure building on a landslide body under a fortification earthquake condition is established. The landslide is a typical loess-mudstone contact surface landslide in a loess area in the northwest, and through means of on-site fine investigation, unmanned aerial vehicle surveying and mapping and drilling, the landslide is disclosed as a landslide body structure which is formed by covering a fourth series of loess and underlying recent series of mudstone, the landslide body structure is a forward rock-soil mass mixed slope, a sliding belt develops along the loess-mudstone contact surface, and a generalized landslide and building coupling numerical analysis model is established in combination with the distribution of buildings on the landslide body (as shown in fig. 6).
(2) Determining input seismic motion intensity parameters: according to the Chinese seismic peak acceleration zoning map (GB18306-2015), the peak acceleration of the landslide is determined to be 0.3g under the seismic fortification condition (the exceeding probability is 10% in 50 years). In this embodiment, an earthquake waveform (horizontal PGA is 0.04g, and vertical PGA is 0.012g) monitored by an XX-year regan county earthquake Tianshui city amalgamation station is selected as a power input reference waveform. Firstly, using Seimosoft software to amplify the reference waveform by 7.5 times (making horizontal PGA 0.3g and vertical PGA 0.09g), then performing low-pass filtering and integral operation on the amplified seismic waves to obtain a seismic velocity time-course curve under a defense seismic scene, and finally performing reference correction on the velocity time-course curve (see fig. 7), and converting the seismic velocity time-course curve into a numerical value by using the following formula to analyze the applied external seismic load:
σ n ＝-2(ρC p )v n (formula 1)
σ s ＝-2(ρC s )v s (formula 2)
In the formula sigma n And σ s Normal and tangential stresses applied to the viscous boundary, C p And C s Wave velocities of P-and S-waves, respectively, where v n And v s The normal and tangential velocity components on the model boundary are respectively, and rho is the density of the rock and soil mass of the sliding bed.
(3) Determining physical and mechanical parameters of landslide rock-soil mass and masonry structure: the field is adopted loess original state appearance, acquires the mudstone sample through the probing, surveys the basic physics mechanics parameter of landslide ground body in the laboratory, mainly includes: the density, the elastic modulus, the Poisson ratio, the cohesive force and the internal friction angle are referred to the collected data and the research results of the same-row in the area, and the rock-soil body parameter value required by the numerical analysis is comprehensively determined (see figure 2). The strength parameter value of the masonry structure building is mainly determined according to masonry structure design specifications (GB50003-2011), the damage of the brick masonry building is generally damaged along the mortar joint section of the masonry, so the tensile strength and the shear strength of the brick masonry building are determined according to different mortar grades, the value standard is shown in figure 3, and the material parameters of the buildings (A, B, C and D) of the two-layer brick masonry structure on the slope body of the embodiment mainly adopt the mortar grade of M7.5 strength grade.
(4) Calculating the power action and the cooperative deformation of the sliding body and the building: according to the established landslide and building coupling numerical analysis model, a discrete unit numerical method (CDEM) based on continuous medium mechanics is adopted to calculate the cooperative deformation rule of the sliding body and the brick masonry building under the condition of a fortification earthquake (see figure 8). Free field boundaries are applied to the left side and the right side of the model, and viscous boundary conditions with normal seismic force and tangential seismic force are applied to the bottom of the model. Rayleigh damping is selected for dynamic calculation damping, and the center frequency f is determined by analyzing a velocity response spectrum min About 2Hz, the minimum critical damping ratio xi is selected min And 3 percent, the rigidity damping coefficient beta is calculated to be 2e-3, and the mass damping coefficient alpha is calculated to be 0.45. After the power calculation is finished, extracting the real-time inclination rate L of the building and the real-time displacement D of the surface of the sliding body where the building is located, establishing a correlation curve (see figure 9) of the surface displacement of the sliding body and the inclination rate of the building under logarithmic coordinates, and showing that the inclination rate of the building and the surface displacement of the sliding body form an obvious exponential function relationship.
(5) Building a vulnerability curve V ═ f (D) of the masonry building: according to a correlation curve (figure 9) of the surface displacement of the sliding body and the inclination rate of the building and a corresponding relation (figure 5) of the limit value of the inclination rate of the building and the vulnerability suggestion value, 15 groups of data of the surface displacement D of the sliding slope and the vulnerability V of the masonry are selected, and finally a physical vulnerability curve (figure 10) of the masonry structure is established according to an improved Weibull function. And the following vulnerability equation is obtained.
It should be noted that this curve is only a vulnerability curve of a brick masonry building (mortar strength M7.5) on a slider under specific seismic conditions (PGA ═ 0.3 g). By using the same idea, the building vulnerability assessment curve of buildings with different inducing factors (earthquake and rainfall conditions with different intensities) and different vulnerabilities (different building structures and different material intensities) can be constructed according to actual requirements.
The invention also provides a system for analyzing the vulnerability of the masonry building on the creeping type landslide body, which comprises:
the surveying and mapping data acquisition module is used for acquiring surveying and mapping data of the target landslide body; the surveying and mapping data comprises a slope body structure, a rock-soil body type, building distribution, a digital ortho-image and a digital elevation model;
the geological model building module is used for building a landslide and building coupling geological model according to the surveying and mapping data of the target landslide body;
the mechanical parameter determining module is used for determining physical and mechanical parameters of rock and soil mass and building masonry structures in the landslide and building coupling geological model;
an induction factor applying module for applying an induction factor to the landslide and building coupled geological model; the inducers include earthquakes and rainfall of varying intensity;
the collaborative deformation calculation module is used for calculating the real-time inclination rate of the building and the real-time displacement of the surface of a sliding body where the building is located under a preset inducing factor by adopting a discrete unit numerical method based on continuous medium mechanics;
the vulnerability curve determining module is used for determining a physical vulnerability curve of the building according to the real-time inclination rate of the building and the real-time displacement of the surface of the sliding body where the building is located;
and the landslide risk assessment module is used for assessing landslide risks according to the physical vulnerability curve of the building.
Preferably, the inducer applying module includes:
the induction factor acquisition unit is used for acquiring a preset seismic waveform;
the seismic waveform processing unit is used for amplifying the preset seismic waveform, performing low-pass filtering and performing integral operation processing to obtain a seismic velocity time-course curve;
and the seismic load determining unit is used for performing reference correction on the seismic velocity time-course curve and obtaining the external seismic load applied to the landslide and building coupling geological model by adopting an induction factor numerical conversion formula.
Preferably, the value conversion formula of the induction factor is as follows:
wherein σ n For normal stresses applied at viscous boundaries, σ s For tangential stresses applied at viscous boundaries, C p Is the wave velocity of P wave, C s Is the wave velocity, v, of the S wave n Is the normal velocity component on the model boundary, v s And the component is tangential velocity component on the model boundary, and rho is the density of the rock and soil mass of the sliding bed.
Preferably, the vulnerability profile determination module includes:
the system comprises a vulnerability value determining unit, a building real-time gradient determining unit and a building real-time gradient determining unit, wherein the vulnerability value determining unit is used for determining the vulnerability value of a corresponding building according to the real-time gradient of the building;
and the fitting unit is used for fitting by adopting an improved Weibull function according to the real-time displacement of the surface of the sliding body where the building is located and the data of the vulnerability value of the corresponding building to obtain a physical vulnerability curve of the building.
Preferably, the modified Weibull function is:
wherein V is physical vulnerability, D is landslide surface displacement, a represents a first fitting parameter, and b represents a second fitting parameter.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
according to the method, a landslide and building coupling geological model is constructed by obtaining the mapping data of a target landslide body, then an inducing factor is applied to the geological model, the real-time inclination rate of a building under the preset inducing factor and the real-time displacement of the surface of the landslide body where the building is located are calculated, and then a building physical vulnerability curve is determined based on the calculation, so that the physical vulnerability of the building can be accurately evaluated under the action of the specific inducing factor, the quantification and the reliability of the physical vulnerability evaluation of the masonry building on the creeping landslide body are greatly improved, and a powerful basis is provided for the risk quantitative evaluation and the accurate disaster reduction of landslide disasters.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other.
The principle and the embodiment of the present invention are explained by applying specific examples, and the above description of the embodiments is only used to help understanding the method and the core idea of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.
Claims (10)
1. A method for analyzing vulnerability of masonry buildings on a creeping type landslide body is characterized by comprising the following steps:
obtaining surveying and mapping data of a target landslide body; the surveying and mapping data comprises a slope body structure, a rock-soil body type, building distribution, a digital ortho-image and a digital elevation model;
constructing a landslide and building coupling geological model according to the surveying and mapping data of the target landslide body;
determining physical and mechanical parameters of rock and soil mass and a building masonry structure in the landslide and building coupling geological model;
applying an inducing factor into the landslide and building coupled geological model; the inducers include earthquakes and rainfall of varying intensity;
calculating the real-time inclination rate of the building and the real-time displacement of the surface of a sliding body where the building is located under a preset inducing factor by adopting a discrete unit numerical method based on continuous medium mechanics;
determining a physical vulnerability curve of the building according to the real-time inclination rate of the building and the real-time displacement of the surface of a sliding body where the building is located;
and evaluating the landslide risk according to the building physical vulnerability curve.
2. The method of claim 1, wherein the applying of the inducing factors to the landslide and building coupled geological model comprises:
acquiring a preset seismic waveform;
amplifying, low-pass filtering and integral operation processing are carried out on the preset seismic waveform to obtain a seismic velocity time-course curve;
and performing reference correction on the seismic velocity time-course curve, and obtaining an external seismic load applied to the landslide and building coupling geological model by adopting an induction factor numerical value conversion formula.
3. The method for analyzing the vulnerability of masonry buildings on creeping type landslide body according to claim 2, wherein the induction factor value conversion formula is as follows:
wherein σ n For normal stresses applied at viscous boundaries, σ s For tangential stresses applied at viscous boundaries, C p Is the wave velocity of P wave, C s Is the wave velocity, v, of the S wave n Is the normal velocity component on the model boundary, v s Is the tangential velocity component on the boundary of the model,rho is the density of the rock-soil body of the sliding bed.
4. The method for analyzing the vulnerability of masonry buildings on a creeping type sliding mass according to claim 1, wherein the determining the physical vulnerability curve of the building according to the real-time tilt rate of the building and the real-time displacement of the surface of the sliding mass on which the building is located comprises:
determining the damage value of the corresponding building according to the real-time tilt rate of the building;
and fitting by adopting an improved Weibull function according to the real-time displacement of the surface of the sliding body where the building is located and the data of the vulnerability value of the corresponding building to obtain a physical vulnerability curve of the building.
5. The method of claim 4, wherein the modified Weibull function is:
wherein V is the physical vulnerability, D is the landslide surface displacement, a represents the first fitting parameter, and b represents the second fitting parameter.
6. A system for analyzing vulnerability of masonry buildings on a creeping type landslide mass, comprising:
the surveying and mapping data acquisition module is used for acquiring surveying and mapping data of the target landslide body; the surveying and mapping data comprises a slope body structure, a rock-soil body type, building distribution, a digital ortho-image and a digital elevation model;
the geological model building module is used for building a landslide and building coupling geological model according to the mapping data of the target landslide body;
the mechanical parameter determining module is used for determining physical and mechanical parameters of rock and soil mass and building masonry structures in the landslide and building coupling geological model;
an induction factor applying module for applying an induction factor to the landslide and building coupled geological model; the inducers include earthquakes and rainfall of varying intensity;
the collaborative deformation calculation module is used for calculating the real-time inclination rate of the building and the real-time displacement of the surface of the sliding body where the building is located under the preset inducing factors by adopting a discrete unit numerical method based on continuous medium mechanics;
the vulnerability curve determining module is used for determining a physical vulnerability curve of the building according to the real-time inclination rate of the building and the real-time displacement of the surface of the sliding body where the building is located;
and the landslide risk assessment module is used for assessing landslide risks according to the physical vulnerability curve of the building.
7. The system of claim 6, wherein the inducer applying module comprises:
the induction factor acquisition unit is used for acquiring a preset seismic waveform;
the seismic waveform processing unit is used for amplifying the preset seismic waveform, performing low-pass filtering and performing integral operation processing to obtain a seismic velocity time-course curve;
and the seismic load determining unit is used for performing reference correction on the seismic velocity time-course curve and obtaining the external seismic load applied to the landslide and building coupling geological model by adopting an induction factor numerical value conversion formula.
8. The system of claim 7, wherein the inducer value conversion formula is:
wherein σ n For normal stresses imposed on viscous boundaries, σ s For applying tangential forces on viscous boundariesForce, C p Is the wave velocity of P wave, C s Is the wave velocity, v, of the S wave n Is the normal velocity component on the model boundary, v s The tangential velocity component on the boundary of the model is shown, and rho is the density of the rock and soil mass of the sliding bed.
9. The system of claim 6, wherein the vulnerability profile determination module comprises:
the system comprises a vulnerability value determining unit, a real-time gradient calculating unit and a real-time gradient calculating unit, wherein the vulnerability value determining unit is used for determining the vulnerability value of a corresponding building according to the real-time gradient of the building;
and the fitting unit is used for fitting by adopting an improved Weibull function according to the real-time displacement of the surface of the sliding body where the building is located and the data of the vulnerability value of the corresponding building to obtain a physical vulnerability curve of the building.
10. The system of claim 9, wherein the modified Weibull function is:
wherein V is physical vulnerability, D is landslide surface displacement, a represents a first fitting parameter, and b represents a second fitting parameter.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210637742.XA CN114925577B (en) | 2022-06-07 | 2022-06-07 | Vulnerability analysis method and system for masonry building on creeping type landslide body |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210637742.XA CN114925577B (en) | 2022-06-07 | 2022-06-07 | Vulnerability analysis method and system for masonry building on creeping type landslide body |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114925577A true CN114925577A (en) | 2022-08-19 |
| CN114925577B CN114925577B (en) | 2023-03-21 |
Family
ID=82812702
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210637742.XA Active CN114925577B (en) | 2022-06-07 | 2022-06-07 | Vulnerability analysis method and system for masonry building on creeping type landslide body |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114925577B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117852894A (en) * | 2024-03-06 | 2024-04-09 | 山东省地震工程研究院 | Rural masonry house bearing capacity assessment method and system |
| CN117951921A (en) * | 2024-03-26 | 2024-04-30 | 云南农业大学 | Slope earthquake vulnerability assessment method |
| CN117972437A (en) * | 2024-03-29 | 2024-05-03 | 四川省建筑设计研究院有限公司 | Regional building earthquake damage prediction method and system aiming at complex terrain geological conditions |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2008202994A (en) * | 2007-02-16 | 2008-09-04 | Public Works Research Institute | Method, device and program of estimating landslide surface shape |
| CN104063616A (en) * | 2014-07-04 | 2014-09-24 | 重庆大学 | Method for evaluating damage loss rate of buildings under action of landslide instability movement |
| CN108596518A (en) * | 2018-05-14 | 2018-09-28 | 中国路桥工程有限责任公司 | A kind of Highway Geological Disaster risk assessment method |
| CN109165424A (en) * | 2018-08-03 | 2019-01-08 | 四川理工学院 | A kind of landslide assessment of easy generation method based on domestic GF-1 satellite data |
| CN111090954A (en) * | 2019-12-17 | 2020-05-01 | 国家电网有限公司 | Geological disaster hidden danger identification and analysis, evaluation and treatment method |
| CN111191880A (en) * | 2019-12-13 | 2020-05-22 | 华能澜沧江水电股份有限公司 | Slope full life cycle safety management method based on digital mapping |
| CN111475978A (en) * | 2020-04-03 | 2020-07-31 | 中国地质科学院地质力学研究所 | Prediction method for damage engineering protection effect after high-position remote landslide |
| CN112241603A (en) * | 2020-10-28 | 2021-01-19 | 中国地质科学院地质力学研究所 | Numerical simulation method for high-order landslide impact scraping and underlayer converging process |
| CN112395789A (en) * | 2020-10-23 | 2021-02-23 | 马培峰 | Method for analyzing urban landslide deformation by coupling InSAR and numerical simulation |
-
2022
- 2022-06-07 CN CN202210637742.XA patent/CN114925577B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2008202994A (en) * | 2007-02-16 | 2008-09-04 | Public Works Research Institute | Method, device and program of estimating landslide surface shape |
| CN104063616A (en) * | 2014-07-04 | 2014-09-24 | 重庆大学 | Method for evaluating damage loss rate of buildings under action of landslide instability movement |
| CN108596518A (en) * | 2018-05-14 | 2018-09-28 | 中国路桥工程有限责任公司 | A kind of Highway Geological Disaster risk assessment method |
| CN109165424A (en) * | 2018-08-03 | 2019-01-08 | 四川理工学院 | A kind of landslide assessment of easy generation method based on domestic GF-1 satellite data |
| CN111191880A (en) * | 2019-12-13 | 2020-05-22 | 华能澜沧江水电股份有限公司 | Slope full life cycle safety management method based on digital mapping |
| CN111090954A (en) * | 2019-12-17 | 2020-05-01 | 国家电网有限公司 | Geological disaster hidden danger identification and analysis, evaluation and treatment method |
| CN111475978A (en) * | 2020-04-03 | 2020-07-31 | 中国地质科学院地质力学研究所 | Prediction method for damage engineering protection effect after high-position remote landslide |
| CN112395789A (en) * | 2020-10-23 | 2021-02-23 | 马培峰 | Method for analyzing urban landslide deformation by coupling InSAR and numerical simulation |
| CN112241603A (en) * | 2020-10-28 | 2021-01-19 | 中国地质科学院地质力学研究所 | Numerical simulation method for high-order landslide impact scraping and underlayer converging process |
Non-Patent Citations (1)
| Title |
|---|
| 车福东: "近远震作用下大型黄土—泥岩滑坡动力响应与稳定性研究", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117852894A (en) * | 2024-03-06 | 2024-04-09 | 山东省地震工程研究院 | Rural masonry house bearing capacity assessment method and system |
| CN117852894B (en) * | 2024-03-06 | 2024-05-24 | 山东省地震工程研究院 | Rural masonry house bearing capacity assessment method and system |
| CN117951921A (en) * | 2024-03-26 | 2024-04-30 | 云南农业大学 | Slope earthquake vulnerability assessment method |
| CN117951921B (en) * | 2024-03-26 | 2024-06-04 | 云南农业大学 | Slope earthquake vulnerability assessment method |
| CN117972437A (en) * | 2024-03-29 | 2024-05-03 | 四川省建筑设计研究院有限公司 | Regional building earthquake damage prediction method and system aiming at complex terrain geological conditions |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114925577B (en) | 2023-03-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114925577B (en) | Vulnerability analysis method and system for masonry building on creeping type landslide body | |
| Chieffo et al. | Comparative fragility methods for seismic assessment of masonry buildings located in Muccia (Italy) | |
| KR101642951B1 (en) | GIS-based real time earthquake prediction method | |
| CN113379105B (en) | Site earthquake liquefaction disaster rapid evaluation method based on artificial intelligence | |
| Ventura et al. | Structural dynamic properties of a reinforced concrete high-rise building during construction | |
| CN104989456B (en) | A kind of Large Span Underground engineering excavation surrounding rock stability monitoring and pre-alarming method | |
| CN111105600A (en) | Cutting slope stability dynamic monitoring and early warning system and method based on rainfall condition | |
| H. Sadeghi et al. | Seismic performance of historical vaulted adobe constructions: a numerical case study from Yazd, Iran | |
| Puppio et al. | Structural irregularity: The analysis of two reinforced concrete (rc) buildings | |
| Li | Comparison of empirical structural vulnerability rapid prediction models considering typical earthquakes | |
| Bouckovalas et al. | Stiff soil amplification effects in the 7 September 1999 Athens (Greece) earthquake | |
| CN111339602A (en) | Method for determining earthquake resistance evaluation index of underground vertical wall vault tunnel structure | |
| Kappos et al. | Analytical estimation of economic loss for buildings in the area struck by the 1999 Athens earthquake and comparison with statistical repair costs | |
| JP2002250027A (en) | Ground investigation method, construction method for base isolation building, and design method for building | |
| CN114925576A (en) | Earthquake landslide risk assessment method and system based on disaster dynamic process | |
| CN115577619A (en) | Method for predicting residual shear strength of liquefied soil based on machine learning | |
| CN115273410A (en) | Sudden landslide monitoring and early warning system based on big data | |
| Du et al. | Methodology for estimating human perception to tremors in high-rise buildings | |
| Biagini et al. | BLWT tests and design loads on the roof of the new Olympic stadium in Piraeus | |
| Murià-Vila et al. | Influence of soil structure interaction and subsoil consolidation in the seismic response of a building | |
| Gatti | Structural height, amplification and damages during the superficial earthquakes at Casamicciola, Ischia Island (2017), and Santa Venerina, Catania (2018), Italy | |
| CN117852894B (en) | Rural masonry house bearing capacity assessment method and system | |
| Kijewski et al. | Full-scale study of the behavior of tall buildings under winds | |
| Marzorati et al. | Detection of local site effects through the estimation of building damages | |
| CN115544794B (en) | Support steel frame design parameter selection method based on fatigue life prediction |
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 |