Method and device for evaluating safety of buildings in land subsidence area

Abstract

The invention relates to the technical field of engineering geology research, in particular to a surface subsidence area building safety assessment method and device, and the method comprises the steps: constructing a three-dimensional calculation model of a stratum and a building of a research area, and the three-dimensional calculation model comprises a stratum model and a building model; determining a target settlement rate of each node on the top surface of the stratum model by using ground synthetic aperture radar interference monitoring data in a target time period of the research area; applying the target settlement rate to the top surface of the stratum model step by step by utilizing an FLAC3D numerical simulation technology so as to determine a settlement displacement field of the three-dimensional calculation model; and on the basis of a final settlement field determined by the settlement displacement field, extracting mechanical response characteristics of a target building in the building model so as to generate a building safety assessment result of the research area in combination with a target safety specification. Therefore, the problems that the land subsidence of the research area is difficult to accurately obtain, the calculation efficiency is low, and the efficiency and the accuracy of building and structure safety assessment are reduced in the prior art are solved.

Description

Technical Field

[0001] The present application relates to the technical field of engineering geological research, and in particular to a method and device for assessing the safety of buildings and structures in land subsidence areas. Background Art

[0002] In related technologies, numerical simulation methods can reveal the development laws of ground subsidence induced by underground mining and its impact on surface engineering structures by establishing a three-dimensional calculation model and combining effective calculation parameters, constitutive models and boundary conditions, thereby evaluating the safety of structures.

[0003] However, using numerical simulation methods to study the development of land subsidence and its impact on engineering structures still presents the challenge of obtaining the physical and mechanical parameters of the rock and soil (such as compression modulus, Poisson's ratio, cohesion, and friction angle). Therefore, to obtain accurate physical and mechanical parameters and make numerical simulation results more realistic, researchers typically use parameter inversion to adjust the calculation parameters.

[0004] However, parameter inversion requires running the numerical model multiple times to iteratively optimize the parameters, and a single calculation of a complex numerical model may take hours or even days, resulting in extremely low inversion efficiency; secondly, different parameter combinations may lead to similar model output results, resulting in the inversion parameters being unable to be uniquely determined, and the noise, sparsity or insufficient temporal and spatial coverage of the observation data will lead to unstable or biased inversion results.

[0005] Therefore, it is still difficult to accurately obtain the ground subsidence in the study area using the above numerical simulation method, and the calculation efficiency is low, which makes it difficult to obtain the deformation and failure characteristics of surface engineering structures in the ground subsidence area, and thus difficult to conduct safety evaluation of surface engineering structures in the ground subsidence area, which needs to be solved urgently. Summary of the Invention

[0006] This application is based on the following problems and understandings made by the inventors:

[0007] Ground subsidence induced by underground mining is a widely recognized geological hazard, characterized by its wide impact, long duration, slow development, complex mechanisms, and significant difficulty in prevention and control. Ground subsidence can weaken the bearing capacity of foundations, causing deformation and damage to surface engineering structures, posing a serious threat to the surrounding environment and human life and property. Therefore, accurately understanding the deformation and damage characteristics and patterns of surface engineering structures in subsidence-prone areas is crucial for scientifically assessing engineering safety risks and developing targeted prevention and control measures.

[0008] With the rapid development of computer technology, numerical simulation methods have been widely used to study these issues due to their economical, convenient, and reliable nature. In related technologies, numerical simulation methods, by establishing three-dimensional computational models and combining effective computational parameters, constitutive models, and boundary conditions, can reveal the development patterns of ground subsidence induced by underground mining and its impact on surface engineering structures, thereby assessing the safety of these structures.

[0009] However, there are still many problems in using numerical simulation methods to study the development laws of ground subsidence and its impact on engineering structures. The main problems are: the physical and mechanical parameters of rock and soil (compression modulus, Poisson's ratio, cohesion and friction angle, etc.) are difficult to obtain and usually difficult to measure directly, especially under large-scale or complex geological conditions. Generally, drilling is used to obtain soil or rock samples to measure physical and mechanical parameters, but due to the existence of size effects, the parameters obtained are not equal to the actual parameters of the rock and soil. At the same time, rock and soil have spatial variability, and physical and mechanical parameters may vary significantly with spatial position. Local measurements alone cannot represent the overall characteristics of the rock and soil.

[0010] In order to obtain accurate physical and mechanical parameters and make the numerical simulation results more realistic, researchers usually use parameter inversion methods to adjust the calculation parameters. The specific process is as follows:

[0011] ① Determine the parameters to be inverted and determine a reasonable range of parameters to avoid inversion results that are out of touch with reality; ② Collect ground subsidence monitoring data in the study area, such as leveling or GNSS (Global Navigation Satellite System) monitoring data; ③ Design a suitable orthogonal table based on the number and levels of the parameters to be inverted, so that each level of each parameter is evenly matched with the levels of other parameters to cover representative combinations in the parameter space; ④ Assign different calculation parameters to the calculation model, start simulation and record the simulation results; ⑤ Compare the simulated ground subsidence with the monitoring data to screen out the appropriate parameter combination.

[0012] However, parameter inversion methods still face the following challenges: ① Parameter inversion requires multiple runs of the numerical model to iteratively optimize parameters. However, a single computation of a complex numerical model can take hours or even days, resulting in extremely low inversion efficiency. ② Different parameter combinations can lead to similar model outputs, making it impossible to uniquely determine the inverted parameters. ③ Noise, sparsity, or insufficient spatial and temporal coverage of the observational data can lead to unstable or biased inversion results.

[0013] Therefore, it is still difficult to accurately obtain the ground subsidence in the study area using the above simulation method, and the calculation efficiency is low, which makes it difficult to obtain the deformation and failure characteristics of surface engineering structures in the ground subsidence area, and thus difficult to conduct safety evaluation, which needs to be solved urgently.

[0014] The present application provides a method and device for safety assessment of buildings and structures in land subsidence areas, in order to solve the problems in related technologies such as difficulty in accurately obtaining the land subsidence in the study area, low calculation efficiency, and reduced efficiency and accuracy of building safety assessment.

[0015] A first aspect of the present application provides a method for safety assessment of buildings and structures in a ground subsidence area, comprising the following steps: constructing a target three-dimensional computational model of strata and buildings in a target study area, wherein the target three-dimensional computational model includes a stratum model and a building model; determining a target settlement rate for each node on the top surface of the stratum model during the target time period using ground synthetic aperture radar interferometry monitoring data of the target study area; applying the target settlement rate of each node to the top surface of the stratum model in steps using target FLAC3D numerical simulation technology to determine a target settlement displacement field of the target three-dimensional computational model; extracting the mechanical response characteristics of the target building in the building model of the target study area based on a final settlement field determined by the target settlement displacement field of the target three-dimensional computational model, so as to generate a safety assessment result of the building in the target study area using the mechanical response characteristics of the target building and the target safety specification.

[0016] Optionally, in one embodiment of the present application, constructing a target three-dimensional computational model of the strata and structures in the target study area includes: determining the target study area and establishing a target three-dimensional computational model of the strata and structures in the target study area; meshing the target three-dimensional computational model and setting the initial geostress field and initial pore water pressure field of the target three-dimensional computational model.

[0017] Optionally, in one embodiment of the present application, the target settlement rate of each node on the top surface of the formation model during the target time period is determined by using the ground synthetic aperture radar interferometry monitoring data of the target study area during the target time period, including: obtaining the surface settlement value of the target study area during the target time period using the ground synthetic aperture radar interferometry monitoring technology; interpolating the ground settlement amount of each node on the top surface of the formation model using the surface settlement value to determine the target settlement value of each node on the top surface of the formation model during the target time period, so as to obtain the target settlement rate based on the target settlement value.

[0018] Optionally, in one embodiment of the present application, the target FLAC3D numerical simulation technology is used to apply the target settlement rate of each node step by step to the top surface of the formation model to determine the target settlement displacement field of the target three-dimensional calculation model, including: based on the target FLAC3D numerical simulation technology, using the surface settlement value of the target study area within the target time period to determine the target node rate loading steps of the formation model; applying the target settlement rate of each node step by step to the nodes on the top surface of the formation model, wherein no conditions are applied to the internal nodes of the formation model; after one iterative calculation of the formation model, the node rate of the top surface of the formation model is cleared, and the surface displacement of the formation model is fixed, and the formation model is iterated multiple times using the target node rate loading steps until the calculation of all time steps is completed, so that the formation model is iterated to equilibrium to form the target settlement displacement field of the target three-dimensional calculation model.

[0019] Optionally, in one embodiment of the present application, the use of the mechanical response characteristics of the target structure and the target safety specifications to generate the building safety assessment results of the target study area includes: analyzing the mechanical response characteristics of the target structure to determine the deformation characteristics and stress characteristics of the target structure; and using the deformation characteristics, the stress characteristics and the target safety specifications of the target structure to evaluate the safety of the buildings in the target study area to generate the building safety assessment results.

[0020] According to a second aspect of the present application, an embodiment provides a device for assessing the safety of buildings and structures in a ground subsidence area, comprising: a construction module for constructing a target three-dimensional computational model of the strata and buildings in a target study area, wherein the target three-dimensional computational model includes a stratum model and a building model; an acquisition module for determining the target settlement rate of each node on the top surface of the stratum model during the target time period using ground synthetic aperture radar interferometry monitoring data of the target study area; a determination module for applying the target settlement rate of each node to the top surface of the stratum model in steps using target FLAC3D numerical simulation technology to determine the target settlement displacement field of the target three-dimensional computational model; and an assessment module for extracting the mechanical response characteristics of the target building in the building model of the target study area based on the final settlement field determined by the target settlement displacement field of the target three-dimensional computational model, so as to generate a safety assessment result of the building in the target study area using the mechanical response characteristics of the target building and the target safety specification.

[0021] Optionally, in one embodiment of the present application, the construction module includes: a determination unit, used to determine the target study area and establish a target three-dimensional computational model of the strata and structures in the target study area; a construction unit, used to grid the target three-dimensional computational model and set the initial geostress field and initial pore water pressure field of the target three-dimensional computational model.

[0022] Optionally, in one embodiment of the present application, the acquisition module includes: an acquisition unit, used to use the ground synthetic aperture radar interferometry monitoring technology to obtain the surface settlement value of the target study area in the target time period; a first determination unit, used to use the surface settlement value to interpolate the ground settlement amount of each node on the top surface of the formation model to determine the target settlement value of each node on the top surface of the formation model in the target time period, so as to obtain the target settlement rate based on the target settlement value.

[0023] Optionally, in one embodiment of the present application, the determination module includes: a second determination unit, used to determine the target node rate loading steps of the formation model based on the target FLAC3D numerical simulation technology and the surface settlement value within the target time period of the target study area; an application unit, used to apply the target settlement rate of each node to the top surface node of the formation model in steps, wherein no conditions are applied to the internal nodes of the formation model; a processing unit, used to clear the node rate of the top surface of the formation model after one iterative calculation of the formation model, fix the surface displacement of the formation model, and iterate the formation model multiple times using the target node rate loading steps until the calculation of all time steps is completed, so that the formation model is iterated to equilibrium to form the target settlement displacement field of the target three-dimensional calculation model.

[0024] Optionally, in one embodiment of the present application, the evaluation module includes: an analysis unit for analyzing the mechanical response characteristics of the target building structure to determine the deformation characteristics and stress characteristics of the target building structure; an evaluation unit for evaluating the safety of the building structure in the target study area using the deformation characteristics of the target building structure, the stress characteristics and target safety specifications to generate the building structure safety evaluation result.

[0025] The third aspect of the present application provides an electronic device, comprising: a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor executes the program to implement a method for assessing the safety of buildings and structures in a ground subsidence area as described in the above embodiment.

[0026] A fourth aspect of the present application provides a computer-readable storage medium, which stores a computer program. When the program is executed by a processor, it implements the above-mentioned method for assessing the safety of buildings and structures in a land subsidence area.

[0027] The fifth embodiment of the present application provides a computer program product, including a computer program, which, when executed, is used to implement the above-mentioned method for assessing the safety of buildings and structures in land subsidence areas.

[0028] The embodiment of the present application can construct a three-dimensional computational model including a stratum model and a building model of the study area, and use the ground InSAR (synthetic aperture radar interferometry) monitoring data of the study area within a certain period of time to determine the target settlement rate of each node on the top surface of the stratum model. Then, using the FLAC3D numerical simulation technology, the target settlement rate is applied step by step to the top surface of the stratum model to determine the settlement displacement field of the three-dimensional computational model, and then the mechanical response characteristics of the target building in the building model are extracted to generate the building safety assessment results of the study area in combination with the target safety specifications, effectively improving the efficiency and accuracy of the building safety assessment. Thus, it solves the problems in the related art that it is difficult to accurately obtain the ground settlement of the study area, and the calculation efficiency is low, which reduces the efficiency and accuracy of the building safety assessment.

[0029] Additional aspects and advantages of the present application will be given in part in the description below, and in part will become apparent from the description below, or will be learned through practice of the present application. BRIEF DESCRIPTION OF THE DRAWINGS

[0030] The above and / or additional aspects and advantages of the present application will become apparent and easily understood from the following description of the embodiments in conjunction with the accompanying drawings, in which:

[0031] Figure 1 A flow chart of a method for safety assessment of buildings and structures in land subsidence areas provided according to an embodiment of the present application;

[0032] Figure 2 A schematic diagram of a stratum model and a railway model of a study area according to a specific embodiment of the present application;

[0033] Figure 3 This is a schematic diagram of setting the node rate on the top surface of the formation model according to a specific embodiment of the present application;

[0034] Figure 4 A schematic diagram of a settlement field of a three-dimensional calculation model in 2016 according to a specific embodiment of the present application;

[0035] Figure 5 A schematic diagram of a settlement field of a three-dimensional calculation model in 2019 according to a specific embodiment of the present application;

[0036] Figure 6 A schematic diagram of a settlement field of a three-dimensional calculation model in 2023 according to a specific embodiment of the present application;

[0037] Figure 7 This is a schematic diagram of the settlement and deformation of the railway subgrade along the direction in 2023 according to a specific embodiment of the present application;

[0038] Figure 8 This is a schematic diagram of the stress distribution along the railway subgrade in 2023 according to a specific embodiment of the present application;

[0039] Figure 9 A schematic diagram of a curve of uneven settlement of the railway subgrade along the strike direction according to a specific embodiment of the present application;

[0040] Figure 10 This is a flow chart of a method for safety assessment of buildings and structures in land subsidence areas according to a specific embodiment of the present application;

[0041] Figure 11 This is a schematic structural diagram of a device for assessing the safety of buildings and structures in land subsidence areas according to an embodiment of the present application;

[0042] Figure 12 A schematic diagram of the structure of an electronic device provided according to an embodiment of the present application. DETAILED DESCRIPTION

[0043] The following describes in detail embodiments of the present application, examples of which are shown in the accompanying drawings, wherein the same or similar reference numerals throughout represent the same or similar elements or elements having the same or similar functions. The embodiments described below with reference to the accompanying drawings are exemplary and are intended to be used to explain the present application, and should not be construed as limiting the present application.

[0044] The following describes a method and device for assessing the safety of buildings and structures in a ground subsidence area according to an embodiment of the present application with reference to the accompanying drawings. In view of the fact that it is difficult to accurately obtain the ground subsidence of the study area in the related art mentioned in the background technology center, and the calculation efficiency is low, which reduces the efficiency and accuracy of the safety assessment of buildings and structures, the present application provides a method for assessing the safety of buildings and structures in a ground subsidence area. In this method, a three-dimensional calculation model including a stratum model and a building model of the study area can be constructed. The ground InSAR monitoring data within a certain period of time in the study area is used to determine the target settlement rate of each node on the top surface of the stratum model. Then, the FLAC3D numerical simulation technology is used to apply the target settlement rate step by step to the top surface of the stratum model to determine the settlement displacement field of the three-dimensional calculation model, and then the mechanical response characteristics of the target building in the building model are extracted to generate the safety assessment results of the building and structures in the study area in combination with the target safety specifications, effectively improving the efficiency and accuracy of the safety assessment of buildings and structures. Thus, the problem that it is difficult to accurately obtain the ground subsidence of the study area in the related art, and the calculation efficiency is low, which reduces the efficiency and accuracy of the safety assessment of buildings and structures is solved.

[0045] Specifically, Figure 1 A flow chart of a method for safety assessment of buildings and structures in land subsidence areas provided in an embodiment of the present application.

[0046] like Figure 1 As shown, a method for safety assessment of buildings and structures in land subsidence areas includes the following steps:

[0047] In step S101 , a target three-dimensional computational model of the strata and buildings in the target study area is constructed, wherein the target three-dimensional computational model includes a stratum model and a building model.

[0048] In the embodiment of the present application, the target study area is an area for conducting safety assessment of buildings and structures in a land subsidence area.

[0049] It can be understood that the embodiments of the present application can construct a three-dimensional computational model of the strata and structures in the study area, wherein the three-dimensional computational model includes a stratum model and a structure model. For example, the present application can use FLAC3D technology to establish a stratum model and a structure model, and mesh the stratum model and the structure model in the following steps, so as to accurately simulate the interaction between structures and strata under complex geological conditions and effectively evaluate the stability and safety of the structures.

[0050] Among them, in one embodiment of the present application, a target three-dimensional computational model of the strata and structures in the target study area is constructed, including: determining the target study area, and establishing a target three-dimensional computational model of the strata and structures in the target study area; meshing the target three-dimensional computational model, and setting the initial geostress field and initial pore water pressure field of the target three-dimensional computational model.

[0051] In the actual implementation process, the embodiment of the present application can select the research scope and establish a three-dimensional calculation model based on the engineering geological survey data, including the stratum calculation model and the building calculation model, and further, the stratum calculation model and the building calculation model are meshed. For example, Figure 2 As shown, this application takes the safety assessment of railway subgrade in a land subsidence active area as an example, and uses FLAC3D technology to establish a stratum calculation model and a railway calculation model. The stratum model size is 300m*300m*100m, and the stratum is divided into ten layers. The layer thickness, stratum lithology and calculation parameters are shown in Table 1. Table 1 is a table of stratum thickness and mechanical parameters of the calculation model. The specific Table 1 is as follows:

[0052] Table 1

[0053]

[0054] Among them, the constitutive model of each soil layer is uniformly set to the classic Mohr-Coulomb model, which can describe the stress-strain characteristics of most rock and soil masses.

[0055] In addition, the railway subgrade structure was simulated and analyzed using solid elements. According to the specifications, the subgrade surface width was set to 21.6m, the top surface width was set to 13.5m, the bottom surface width was set to 21.6m, the subgrade thickness was set to 2.7m, and the slope was set to 1:1.5. The calculation parameters of the railway subgrade are shown in Table 2, which is a table of the physical and mechanical parameters of the railway subgrade unit. The specific table 2 is as follows:

[0056] Table 2

[0057] Roadbed structure model E(MPa) μ Thickness (m) elastic 250 0.25 2.7

[0058] The constitutive model of the railway subgrade is an isotropic elastic model. The railway calculation model uses triangular prism elements to perform meshing in the following steps. After meshing, the number of nodes and the number of elements of the railway calculation model are 85,030 and 158,014, respectively.

[0059] Furthermore, in order to effectively study the impact of ground subsidence on engineering structures, it is first necessary to initialize the ground stress field of the three-dimensional calculation model. Among them, this application adopts a staged elastic-plastic solution method to simulate the initial ground stress field. During the calculation process, a fixed constraint is applied to the lower boundary of the three-dimensional calculation model, that is, no displacement will occur in the x, y, and z directions; the front, back, left, and right boundaries of the three-dimensional calculation model only constrain the normal displacement, and the upper boundary is the ground surface, which is not constrained for the time being; the gravity acceleration is set to -9.8m / s 2 The specific implementation method is as follows:

[0060] First, the constitutive model of the geological material is set as the Mohr-Coulomb model, and its cohesion and tensile strength are set to 1000 MPa. Under the action of gravity alone, iterative calculations are performed until the ratio of the maximum unbalanced force to the typical internal force is less than 10. -5 Then, the material cohesion and tensile strength are reset to the initial setting value of the soil layer, and the iterative calculation is repeated until the ratio of the maximum unbalanced force to the typical internal force is less than 10. -5 , the three-dimensional calculation model stress field is initialized.

[0061] Furthermore, to better simulate the ground subsidence displacement field, it is necessary to set the initial pore water pressure field for the 3D computational model. Based on hydrogeological survey data, the initial water level of the stratum model is set using the water table command. This also sets the hydraulic head boundary conditions for the stratum model, automatically generating a pore water pressure field within the stratum model that is linearly distributed according to the water level. The unit permeability coefficients are set according to the soil properties of each stratum in the stratum model. The mechanics calculation mode for the stratum model is disabled, and the fluid calculation mode is enabled. Balance calculations are performed to the set unbalanced flow ratio to generate the initial seepage field.

[0062] During the simulation process, the groundwater level remains constant for a given time period. When entering the next time period, the groundwater level needs to be readjusted based on the specific groundwater level data for that time period, thereby generating a new seepage field. This method accounts for changes in the groundwater level, making the evolution of the settlement field in the stratigraphic model used in this application more realistic.

[0063] In step S102, the target subsidence rate of each node on the top surface of the stratum model in the target time period is determined using the ground synthetic aperture radar interferometry monitoring data of the target study area in the target time period.

[0064] In the embodiment of the present application, the target time period may be several years, several months, several weeks, etc., which is specifically set by relevant technical personnel and is not specifically limited here.

[0065] It can be understood that the embodiment of the present application can use InSAR (synthetic aperture radar interferometry) technology to obtain ground monitoring data for a certain time period within the study scope, that is, surface settlement values, so as to determine the ground settlement values of each node on the top surface of the formation model within a certain time period within the study scope, and then obtain the settlement rate of each node on the top surface of the formation model based on the ground settlement values, thereby effectively improving the ability to capture the dynamic change response of the formation and the analysis accuracy.

[0066] Among them, in one embodiment of the present application, the target settlement rate of each node on the top surface of the formation model in the target time period is determined by using the ground synthetic aperture radar interferometry monitoring data of the target study area in the target time period, including: using the ground synthetic aperture radar interferometry monitoring technology to obtain the surface settlement value of the target study area in the target time period; using the surface settlement value to interpolate the ground settlement amount of each node on the top surface of the formation model to determine the target settlement value of each node on the top surface of the formation model in the target time period, so as to obtain the target settlement rate according to the target settlement value.

[0067] For example, the embodiment of the present application adopts the surface settlement values of a certain area in the North China Plain from 2016 to 2023, and uses the Kriging interpolation method to interpolate the ground settlement amounts of each node on the top surface of the stratum model, thereby determining the ground settlement values of each node on the top surface of the stratum model, so as to achieve the purpose of finely restoring the surface settlement displacement field of the stratum model. While ensuring the accuracy of the ground settlement, the simulation time is greatly shortened and the complexity of the calculation is effectively reduced.

[0068] In FLAC3D, displacement is applied in the form of node settlement rate. Therefore, the ground settlement of each node within a certain time period is expressed in the form of node settlement rate Δ, with the unit of m / s.

[0069] In step S103, the target settlement rate of each node is applied step by step to the top surface of the formation model using the target FLAC3D numerical simulation technology to determine the target settlement displacement field of the target three-dimensional calculation model.

[0070] It can be understood that the embodiment of the present application can use the target FLAC3D numerical simulation technology to evenly divide the settlement rate of each node within a certain time period into multiple steps, and apply them step by step to the top surface of the formation model. For example, the ground InSAR monitoring data distribution from 2016 to 2023 can be applied to the surface of the formation model to determine the settlement displacement field of the three-dimensional calculation model, and then the ground settlement field between 2016 and 2023 can be formed, which effectively ensures the accuracy and reliability of the calculation.

[0071] In one embodiment of the present application, the target FLAC3D numerical simulation technology is used to apply the target settlement rate of each node to the top surface of the formation model in steps to determine the target settlement displacement field of the target three-dimensional calculation model, including: based on the target FLAC3D numerical simulation technology, the surface settlement value within the target time period of the target study area is used to determine the target node rate loading steps of the formation model; the target settlement rate of each node is applied step by step to the nodes on the top surface of the formation model, wherein no conditions are applied to the nodes inside the formation model; after one iterative calculation of the formation model, the node rate of the top surface of the formation model is cleared, and the surface displacement of the formation model is fixed, and the formation model is iterated multiple times using the target node rate loading steps until the calculation of all time steps is completed, so that the formation model is iterated to equilibrium to form the target settlement displacement field of the target three-dimensional calculation model.

[0072] As a possible implementation method, the embodiment of the present application can evenly divide the node settlement rate Δ within a certain time period into multiple steps based on the target FLAC3D numerical simulation technology and apply it to the surface of the formation model. Figure 3 As shown in Figure 1, the number of loading steps N for the node settlement rate Δ is determined. First, the node rates of the top surface, i.e., the settlement rate of each node on the top surface Δ1 = Δ / N, are applied to the nodes on the top surface of the formation model, and no other conditions are applied to the internal nodes. Then, the formation model is iterated for one step. Finally, the node rates of the top surface of the formation model are reset to zero and the surface displacement is fixed. Based on the number of loading steps N for the node settlement rate Δ, the formation model is iterated to equilibrium. The top surface displacement field is transmitted to the interior of the formation model, and the internal nodes respond and generate displacements, forming the settlement displacement field of the formation model. For details, see Figure 3 .

[0073] Then, after the previous calculation is completed, the node settlement rate Δ2 = Δ / N is continued to be applied to the top of the formation model, and no other conditions are applied to the internal nodes; then, the formation model is iterated for one step, and then the rates of all nodes in the formation model are reset to zero, and the displacement of the top surface of the formation model is fixed again, so that the formation model is iterated to equilibrium and a new settlement displacement field of the formation model is formed.

[0074] Repeat the above iterative process and apply the InSAR monitoring data from 2016 to 2023 to the model surface to form the following Figure 4 、 Figure 5 and Figure 6 The ground subsidence field between 2016 and 2023 is shown, which can accurately simulate the dynamic changes of the strata over time and ensure the accuracy and stability of the subsidence displacement field.

[0075] In other words, the present embodiment can apply InSAR monitoring data as displacement boundary conditions on top of the computational model. This allows for rapid and accurate calculation of the ground subsidence field in the study area, avoiding the complex and time-consuming process of ground subsidence simulation and significantly simplifying the simulation. Furthermore, by acquiring InSAR monitoring data from different study areas and time periods, the development patterns of ground subsidence in different regions and time periods, as well as its impact on engineering structures, can be studied, significantly improving the applicability of the building safety assessment method.

[0076] In step S104, based on the final settlement field determined by the target settlement displacement site of the target three-dimensional calculation model, the mechanical response characteristics of the target building structure in the building model of the target study area are extracted to generate the building safety assessment results of the target study area using the mechanical response characteristics of the target building structure and the target safety specifications.

[0077] In the embodiment of the present application, the target building is a building in a ground subsidence area, and the present application uses a railway subgrade as an example for explanation.

[0078] It can be understood that the embodiment of the present application can extract the mechanical response characteristics of the buildings and structures in the building model based on the final settlement field determined by the settlement displacement field of the three-dimensional calculation model. For example, in the process of calculating the ground settlement field in the embodiment of the present application, the railway model on the surface synchronously generates a calculation response. Therefore, after the ground settlement field calculation in the above steps is completed, the mechanical response characteristics of the railway structure can be directly obtained, such as settlement deformation characteristics, stress distribution characteristics, etc., so that the mechanical response characteristics of the buildings and structures in the following steps can be combined with certain safety specifications to establish a safety evaluation standard for the buildings and structures, and then evaluate the safety of the buildings and structures in the study area, and obtain the safety evaluation results of the buildings and structures, which effectively improves the efficiency and accuracy of the safety evaluation of buildings and structures in the ground settlement area, helps to accurately grasp the deformation and destruction characteristics and laws of engineering buildings in the settlement activity area, and is of great significance for scientifically evaluating engineering safety risks and formulating targeted prevention and control measures.

[0079] Optionally, in one embodiment of the present application, the mechanical response characteristics of the target structure and the target safety specifications are used to generate a building safety assessment result for the target study area, including: analyzing the mechanical response characteristics of the target structure to determine the deformation characteristics and stress characteristics of the target structure; and evaluating the safety of the buildings in the target study area using the deformation characteristics, stress characteristics and target safety specifications of the target structure to generate a building safety assessment result.

[0080] For example, the embodiment of the present application can analyze the mechanical response characteristics of the railway structure to obtain the following Figure 7 The settlement deformation characteristics of the railway structure shown in Figure 8 The stress distribution characteristics shown are then obtained. After obtaining the deformation characteristics of the building structure, that is, the railway structure, the safety of the railway structure is evaluated in combination with safety regulations. The specific method is as follows:

[0081] Taking ballastless track as an example, my country's standards for limiting subgrade settlement after ballastless track construction clearly specify the track deformation limit. Table 3 is a standard table of subgrade settlement limits after ballastless track construction. The specific details of Table 3 are as follows:

[0082] Table 3

[0083] Post-construction settlement is generally allowed Subgrade length ≥ 20m allows post-construction settlement Uneven settlement Differential settlement ≤15mm ≤30mm 20mm / 20m ≤5mm

[0084] If all items are less than the allowable values in Table 3, the railway track is considered to be basically stable. The specific deformation standards are as follows:

[0085] In order to reasonably analyze the impact of ground subsidence on railway subgrade, this application extracts the Z-direction deformation values of 50 points along the railway subgrade surface and calculates the inclined deformation along the subgrade. This value can represent the degree of uneven settlement. The calculation formula of the inclined deformation i is as follows:

[0086]

[0087] Where Δdisz is the difference in Z-direction displacement between the extraction points, and Δl is the distance between the two extraction points.

[0088] For example, Figure 9 As shown in Figure 2, there is obvious uneven settlement of the railway subgrade surface, where the left half of the model shows settlement and the right half shows uplift. Figure 8 The results show that the differential settlement of the railway subgrade gradually increased over time between 2016 and 2022, but by 2023, the differential settlement had significantly decreased, indicating a trend toward uniform settlement. From 2016 to 2023, the differential settlement of the top and bottom surfaces of the railway subgrade remained below 0.8 mm / m, remaining within the limits specified for ballastless track subgrades, indicating that the track remains generally stable.

[0089] For example, if Figure 10 As shown, the working principle of the embodiment of the present application is described in detail below with a specific embodiment.

[0090] Step S1001: Establish a three-dimensional computational model of the strata and structures. That is, the embodiment of the present application can construct a target three-dimensional computational model of the strata and structures in the study area, wherein the target three-dimensional computational model includes a stratum model and a structure model.

[0091] Step S1002: Calculate the initial geostress field and the pore water pressure field. That is, the embodiment of the present application can set the initial geostress field and the initial pore water pressure field of the target three-dimensional calculation model.

[0092] Step S1003: Use ground synthetic aperture radar interferometry monitoring data to obtain the settlement rate of the top surface nodes of the formation model. That is, the embodiment of the present application can use ground synthetic aperture radar interferometry monitoring data within a certain period of time in the study area to determine the settlement rate of the top surface nodes of the formation model in the target three-dimensional calculation model.

[0093] Step S1004: Apply the node settlement rate in steps to the surface of the formation model to form a settlement displacement field. That is, the embodiment of the present application can use the target FLAC3D numerical simulation technology to apply the node settlement rate in steps to the surface of the formation model to form the settlement displacement field of the target three-dimensional calculation model.

[0094] Step S1005: Repeat steps S1003-S1004 to obtain the final ground subsidence field. That is, the embodiment of the present application can use the node rate loading steps to iterate the stratum model multiple times until the calculation of all time steps is completed, so that the stratum model is iterated to equilibrium to form the final ground subsidence field of the target three-dimensional calculation model.

[0095] Step S1006: Extract the mechanical response characteristics of the building and structure and analyze them. That is, the embodiment of the present application can extract the mechanical response characteristics of the building and structure in the building model of the study area, such as deformation characteristics and stress characteristics, based on the final settlement field determined by the settlement displacement field of the target three-dimensional calculation model.

[0096] Step S1007: Evaluate the safety of buildings and structures in combination with safety regulations. That is, the embodiment of the present application can use the deformation characteristics, stress characteristics and safety regulations of buildings and structures to evaluate the safety of buildings and structures in the study area to generate building safety assessment results, effectively improving the efficiency and accuracy of building safety assessments in ground subsidence areas.

[0097] According to a method for assessing the safety of buildings and structures in a ground subsidence area proposed in an embodiment of the present application, a three-dimensional computational model including a stratum model and a building model of the study area can be constructed. The target settlement rate of each node on the top surface of the stratum model is determined using ground InSAR monitoring data within a certain period of time in the study area. Then, the target settlement rate is applied step by step to the top surface of the stratum model using FLAC3D numerical simulation technology to determine the settlement displacement field of the three-dimensional computational model, and then the mechanical response characteristics of the target building in the building model are extracted to generate the safety assessment results of the buildings and structures in the study area in combination with the target safety specifications, thereby effectively improving the efficiency and accuracy of the safety assessment of the buildings and structures. Thus, the problems in the related art of being difficult to accurately obtain the ground subsidence of the study area, and the low computational efficiency, which reduces the efficiency and accuracy of the safety assessment of the buildings and structures, are solved.

[0098] Next, a device for assessing the safety of buildings and structures in land subsidence areas according to an embodiment of the present application will be described with reference to the accompanying drawings.

[0099] Figure 11 It is a block diagram of a device for assessing the safety of buildings and structures in land subsidence areas according to an embodiment of the present application.

[0100] like Figure 11 As shown, a device 10 for assessing the safety of buildings and structures in land subsidence areas includes: a construction module 100 , an acquisition module 200 , a determination module 300 and an assessment module 400 .

[0101] Specifically, the construction module 100 is used to construct a target three-dimensional computational model of the strata and buildings in the target research area, wherein the target three-dimensional computational model includes a stratum model and a building model.

[0102] The acquisition module 200 is used to determine the target subsidence rate of each node on the top surface of the stratum model in the target time period by using the ground synthetic aperture radar interferometry monitoring data of the target study area in the target time period.

[0103] The determination module 300 is used to apply the target settlement rate of each node to the top surface of the formation model in steps using the target FLAC3D numerical simulation technology to determine the target settlement displacement field of the target three-dimensional calculation model.

[0104] The evaluation module 400 is used to extract the mechanical response characteristics of the target buildings in the building model of the target study area based on the final settlement field determined by the target settlement displacement site of the target three-dimensional calculation model, so as to generate the safety evaluation results of the buildings in the target study area using the mechanical response characteristics of the target buildings and the target safety specifications.

[0105] Optionally, in one embodiment of the present application, the construction module 100 includes: a determination unit and a construction unit.

[0106] The determination unit is used to determine the target study area and establish a target three-dimensional calculation model of the strata and structures in the target study area.

[0107] The construction unit is used to mesh the target three-dimensional calculation model and set the initial ground stress field and initial pore water pressure field of the target three-dimensional calculation model.

[0108] Optionally, in one embodiment of the present application, the acquisition module 200 includes: an acquisition unit and a first determination unit.

[0109] The acquisition unit is used to obtain the surface subsidence value of the target study area during the target time period using ground synthetic aperture radar interferometry monitoring technology.

[0110] The first determining unit is used to interpolate the ground settlement amount of each node on the top surface of the stratum model using the surface settlement value to determine the target settlement value of each node on the top surface of the stratum model in the target time period, so as to obtain the target settlement rate according to the target settlement value.

[0111] Optionally, in one embodiment of the present application, the determination module 300 includes: a second determination unit, an applying unit, and a processing unit.

[0112] The second determining unit is used to determine the target node rate loading steps of the stratum model based on the target FLAC3D numerical simulation technology and the surface settlement value within the target time period of the target study area.

[0113] The applying unit is used to apply the target settlement rate of each node to the top surface nodes of the formation model in steps, wherein no conditions are applied to the internal nodes of the formation model.

[0114] The processing unit is used to clear the node rate of the top surface of the formation model to zero after one step of iterative calculation of the formation model, fix the surface displacement of the formation model, and iterate the formation model multiple times using the target node rate loading step number until the calculation of all time steps is completed, so that the formation model is iterated to equilibrium to form the target settlement displacement field of the target three-dimensional calculation model.

[0115] Optionally, in one embodiment of the present application, the evaluation module 400 includes: an analysis unit and an evaluation unit.

[0116] The analysis unit is used to analyze the mechanical response characteristics of the target building structure to determine the deformation characteristics and stress characteristics of the target building structure.

[0117] The evaluation unit is used to evaluate the safety of buildings and structures in the target study area by using the deformation characteristics, stress characteristics and target safety specifications of the target buildings and structures to generate building safety evaluation results.

[0118] It should be noted that the above explanation of the embodiment of a method for safety assessment of buildings and structures in a land subsidence area is also applicable to the device for safety assessment of buildings and structures in a land subsidence area of this embodiment, which will not be repeated here.

[0119] According to the device for assessing the safety of buildings and structures in a ground subsidence area proposed in an embodiment of the present application, a three-dimensional computational model including a stratum model and a building model of the study area can be constructed. The target settlement rate of each node on the top surface of the stratum model is determined by using the ground InSAR monitoring data within a certain period of time in the study area. Then, the target settlement rate is applied step by step to the top surface of the stratum model using the FLAC3D numerical simulation technology to determine the settlement displacement field of the three-dimensional computational model, and then the mechanical response characteristics of the target building in the building model are extracted to generate the safety assessment results of the buildings and structures in the study area in combination with the target safety specifications, thereby effectively improving the efficiency and accuracy of the safety assessment of the buildings and structures. Thus, the problems in the related art of being difficult to accurately obtain the ground subsidence of the study area, and the low computational efficiency, which reduces the efficiency and accuracy of the safety assessment of the buildings and structures, are solved.

[0120] Figure 12 This is a schematic diagram of the structure of an electronic device provided in an embodiment of the present application. The electronic device may include:

[0121] A memory 1201 , a processor 1202 , and a computer program stored in the memory 1201 and executable on the processor 1202 .

[0122] When the processor 1202 executes the program, a method for assessing the safety of buildings and structures in land subsidence areas provided in the above embodiment is implemented.

[0123] Furthermore, the electronic device further includes:

[0124] The communication interface 1203 is used for communication between the memory 1201 and the processor 1202 .

[0125] The memory 1201 is used to store computer programs that can be run on the processor 1202 .

[0126] The memory 1201 may include a high-speed RAM memory, and may also include a non-volatile memory (non-volatile memory), such as at least one disk memory.

[0127] If the memory 1201, processor 1202, and communication interface 1203 are implemented independently, the communication interface 1203, memory 1201, and processor 1202 can be interconnected via a bus and communicate with each other. The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, or an Extended Industry Standard Architecture (EISA) bus. Buses can be divided into address buses, data buses, control buses, etc. For ease of representation, Figure 12 Only one thick line is used in the diagram, but this does not mean that there is only one bus or one type of bus.

[0128] Optionally, in a specific implementation, if the memory 1201, the processor 1202 and the communication interface 1203 are integrated on a chip, the memory 1201, the processor 1202 and the communication interface 1203 can communicate with each other through an internal interface.

[0129] The processor 1202 may be a central processing unit (CPU), an application specific integrated circuit (ASIC), or one or more integrated circuits configured to implement the embodiments of the present application.

[0130] This embodiment also provides a computer-readable storage medium having a computer program stored thereon. When the program is executed by a processor, the method for safety assessment of buildings and structures in land subsidence areas as described above is implemented.

[0131] This embodiment also provides a computer program product, including a computer program. When the computer program is executed, it is used to implement the above-mentioned method for assessing the safety of buildings and structures in land subsidence areas.

[0132] In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described can be combined in any one or N embodiments or examples in a suitable manner. In addition, those skilled in the art can combine and combine different embodiments or examples described in this specification and features of different embodiments or examples without contradiction.

[0133] Furthermore, the terms "first" and "second" are used for descriptive purposes only and should not be understood to indicate or imply relative importance or implicitly specify the number of technical features indicated. Thus, a feature specified as "first" or "second" may explicitly or implicitly include at least one such feature. In the description of this application, "N" means at least two, for example, two, three, etc., unless otherwise specifically defined.

[0134] Any process or method description in a flowchart or otherwise described herein may be understood to represent a module, fragment or portion of code comprising one or N executable instructions for implementing a custom logical function or process step, and the scope of the preferred embodiments of the present application includes alternative implementations in which functions may be performed in a different order than shown or discussed, including performing functions in a substantially simultaneous manner or in a reverse order depending on the functions involved, which should be understood by those skilled in the art to which the embodiments of the present application pertain.

[0135] The logic and / or steps represented in the flowcharts or otherwise described herein, for example, can be considered as a sequenced list of executable instructions for implementing the logical functions, and can be embodied in any computer-readable medium for use by, or in conjunction with, an instruction execution system, apparatus, or device (e.g., a computer-based system, a system including a processor, or other system that can fetch and execute instructions from an instruction execution system, apparatus, or device). For purposes of this specification, a "computer-readable medium" can be any device that can contain, store, communicate, propagate, or transport a program for use by, or in conjunction with, an instruction execution system, apparatus, or device. More specific examples (a non-exhaustive list) of computer-readable media include the following: an electrical connection with one or N wires (electronic devices), a portable computer disk cartridge (magnetic device), a random access memory (RAM), a read-only memory (ROM), an erasable and programmable read-only memory (EPROM or flash memory), a fiber optic device, and a portable compact disc read-only memory (CDROM). In addition, the computer-readable medium may even be paper or other suitable medium on which the program is printed, since the program can be obtained electronically by optically scanning the paper or other medium and then editing, interpreting or processing it in other suitable ways as necessary, and then storing it in a computer memory.

[0136] It should be understood that various parts of the present application can be implemented using hardware, software, firmware, or a combination thereof. In the above embodiment, the N steps or methods can be implemented using software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented using hardware, as in another embodiment, it can be implemented using any one or a combination of the following technologies known in the art: a discrete logic circuit having a logic gate circuit for implementing a logic function on a data signal, an application-specific integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0137] Those skilled in the art will understand that all or part of the steps in the method of the above embodiment can be completed by instructing related hardware through a program, and the program can be stored in a computer-readable storage medium. When the program is executed, it includes one or a combination of the steps of the method embodiment.

[0138] In addition, the functional units in the various embodiments of the present application may be integrated into a processing module, or each unit may exist physically separately, or two or more units may be integrated into a module. The above-mentioned integrated module may be implemented in the form of hardware or in the form of a software functional module. If the integrated module is implemented in the form of a software functional module and sold or used as an independent product, it may also be stored in a computer-readable storage medium.

[0139] The storage medium mentioned above may be a read-only memory, a magnetic disk, or an optical disk, etc. Although the embodiments of the present application have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting the present application. Persons skilled in the art may make changes, modifications, substitutions, and variations to the above embodiments within the scope of the present application.

Claims

1. A method for safety assessment of buildings and structures in land subsidence areas, characterized in that: The following steps are involved: Constructing a target three-dimensional computational model of the strata and structures in the target study area, wherein the target three-dimensional computational model includes a stratum model and a structure model; Determining a target subsidence rate for each node on the top surface of the stratum model during the target time period using ground synthetic aperture radar interferometry monitoring data of the target study area during the target time period; Using the target FLAC3D numerical simulation technology, the target settlement rate of each node is applied step by step to the top surface of the formation model to determine the target settlement displacement field of the target three-dimensional calculation model; Based on the final settlement field determined by the target settlement displacement site of the target three-dimensional calculation model, the mechanical response characteristics of the target building in the building model of the target study area are extracted to generate the building safety assessment results of the target study area using the mechanical response characteristics of the target building and the target safety specifications.

2. The method according to claim 1, characterized in that The target three-dimensional computational model of the strata and structures in the target study area is constructed, including: Determine the target study area and establish a target three-dimensional computational model of the strata and structures in the target study area; The target three-dimensional calculation model is meshed, and an initial geostress field and an initial pore water pressure field of the target three-dimensional calculation model are set.

3. The method according to claim 1, characterized in that Determining the target subsidence rate of each node on the top surface of the stratum model during the target time period by using ground synthetic aperture radar interferometry monitoring data of the target study area during the target time period includes: Using the ground synthetic aperture radar interferometry monitoring technology to obtain the surface subsidence value of the target study area during the target time period; The surface settlement value is used to interpolate the ground settlement amount of each node on the top surface of the stratum model to determine the target settlement value of each node on the top surface of the stratum model in the target time period, so as to obtain the target settlement rate according to the target settlement value.

4. The method according to claim 3, characterized in that The target settlement rate of each node is applied step by step to the top surface of the formation model using the target FLAC3D numerical simulation technology to determine the target settlement displacement field of the target three-dimensional calculation model, including: Based on the target FLAC3D numerical simulation technology, the target node rate loading steps of the formation model are determined using the surface settlement value of the target study area within the target time period; Applying the target settlement rate of each node to the top surface nodes of the formation model in steps, wherein no conditions are applied to the internal nodes of the formation model; After one step of iterative calculation of the formation model, the node rate of the top surface of the formation model is cleared, and the surface displacement of the formation model is fixed. The formation model is iterated multiple times using the target node rate loading step number until the calculation of all time steps is completed, so that the formation model is iterated to equilibrium to form the target settlement displacement field of the target three-dimensional calculation model.

5. The method according to claim 1, wherein The generating of the building safety assessment results of the target study area by using the mechanical response characteristics of the target building and structure and the target safety specifications includes: Analyzing the mechanical response characteristics of the target building structure to determine the deformation characteristics and stress characteristics of the target building structure; The safety of the buildings and structures in the target study area is evaluated using the deformation characteristics, the stress characteristics and the target safety specifications of the target buildings and structures to generate the building and structure safety evaluation results.

6. A device for assessing the safety of buildings and structures in land subsidence areas, characterized in that: include: A construction module is used to construct a target three-dimensional computational model of the strata and buildings in the target study area, wherein the target three-dimensional computational model includes a stratum model and a building model; an acquisition module, configured to determine a target subsidence rate for each node on the top surface of the stratum model during a target period of time using ground synthetic aperture radar interferometry monitoring data of the target study area during a target period of time; a determination module for applying the target settlement rate of each node to the top surface of the formation model in steps using a target FLAC3D numerical simulation technology to determine a target settlement displacement field of the target three-dimensional computational model; An evaluation module is used to extract the mechanical response characteristics of the target building structure in the building model of the target study area based on the final settlement field determined by the target settlement displacement field of the target three-dimensional calculation model, so as to generate the building safety evaluation results of the target study area by using the mechanical response characteristics of the target building structure and the target safety specifications.

7. The device according to claim 6, characterized in that The building blocks include: a determination unit, configured to determine the target study area and establish a target three-dimensional computational model of the strata and structures in the target study area; The construction unit is used to perform grid division on the target three-dimensional calculation model and set the initial ground stress field and initial pore water pressure field of the target three-dimensional calculation model.

8. The device according to claim 6, characterized in that The acquisition module includes: an acquiring unit, configured to acquire the surface subsidence value of the target study area during the target time period by using the ground synthetic aperture radar interferometry monitoring technology; The first determination unit is used to interpolate the ground settlement amount of each node on the top surface of the stratum model using the surface settlement value to determine the target settlement value of each node on the top surface of the stratum model in the target time period, so as to obtain the target settlement rate according to the target settlement value.

9. The device according to claim 8, characterized in that The determination module includes: A second determining unit is configured to determine a target node rate loading step number of the stratum model using the surface settlement value of the target study area within the target time period based on a target FLAC3D numerical simulation technology; an applying unit, configured to apply the target settlement rate of each node to the nodes on the top surface of the formation model in steps, wherein no conditions are applied to the nodes inside the formation model; The processing unit is used to clear the node rate of the top surface of the formation model to zero after one step of iterative calculation of the formation model, fix the surface displacement of the formation model, and iterate the formation model multiple times using the target node rate loading step number until the calculation of all time steps is completed, so that the formation model is iterated to equilibrium to form the target settlement displacement field of the target three-dimensional calculation model.

10. The device according to claim 6, characterized in that The evaluation module includes: An analysis unit, configured to analyze the mechanical response characteristics of the target building or structure to determine the deformation characteristics and stress characteristics of the target building or structure; An evaluation unit is used to evaluate the safety of the buildings and structures in the target study area by using the deformation characteristics, the stress characteristics and the target safety specifications of the target buildings and structures to generate the building and structure safety evaluation result.

11. An electronic device, characterized in that: include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement a method for safety assessment of buildings and structures in a land subsidence area as described in any one of claims 1 to 5.

12. A computer-readable storage medium having a computer program stored thereon, characterized in that: The program is executed by a processor to implement a method for safety assessment of buildings and structures in land subsidence areas as described in any one of claims 1 to 5.

13. A computer program product comprising a computer program, characterized in that The computer program is executed by a processor to implement a method for safety assessment of buildings and structures in land subsidence areas as described in any one of claims 1 to 5.

Images