A method and device for evaluating the safety of a building in a land subsidence area

Abstract

The application relates to the technical field of engineering geology research, in particular to a ground subsidence area building structure safety evaluation method and device, which comprises the following steps: constructing a three-dimensional calculation model of strata and building structures in a research area, including a strata model and a building structure model; determining the target subsidence rate of each node on the top surface of the strata model by using ground synthetic aperture radar interferometric monitoring data of a target time period of the research area; applying the target subsidence rate step by step to the top surface of the strata model by using FLAC3D numerical simulation technology to determine the subsidence displacement field of the three-dimensional calculation model; and extracting the mechanical response characteristics of a target building structure in the building structure model based on the final subsidence field determined by the subsidence displacement field to generate the building structure safety evaluation result of the research area in combination with a target safety specification. Therefore, the problems that the ground subsidence of the research area cannot be accurately obtained in the related art, the calculation efficiency is low, and the efficiency and accuracy of building structure safety evaluation are reduced are solved.

Description

Technical Field

[0001] This application relates to the field of engineering geological research technology, and in particular to a method and apparatus for safety assessment of buildings and structures in land subsidence areas. Background Technology

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

[0003] However, using numerical simulation methods to study the development patterns of land subsidence and its impact on engineering structures still faces the challenge of obtaining the physical and mechanical parameters of the soil and rock mass (compression modulus, Poisson's ratio, cohesion, and friction angle, etc.). Therefore, in order to obtain accurate physical and mechanical parameters and make the numerical simulation results more consistent with reality, researchers usually use parameter inversion methods to adjust the calculation parameters.

[0004] However, parameter inversion requires running the numerical model multiple times to iteratively optimize the parameters. A single calculation of a complex numerical model can take hours or even days, resulting in extremely low inversion efficiency. Secondly, different parameter combinations may lead to similar model output results, making it impossible to uniquely determine the parameters obtained from the inversion. Furthermore, noise, sparsity, or insufficient spatiotemporal coverage of the observation data can lead to unstable or biased inversion results.

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

[0006] This application is based on the inventor's understanding and insights into the following issues:

[0007] Ground subsidence induced by underground mining is a major geological hazard issue, characterized by its wide impact, long duration, slow development, complex formation mechanisms, and difficulty in prevention and control. Ground subsidence weakens the bearing capacity of the foundation, leading to 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 failure characteristics and patterns of surface engineering structures in subsidence activity 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 the aforementioned problems due to their economy, convenience, and reliability. Among related technologies, numerical simulation methods, by establishing a three-dimensional computational model 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, many challenges remain in using numerical simulation methods to study the development patterns of land subsidence and its impact on engineering structures. These challenges primarily stem from the difficulty in obtaining the physical and mechanical parameters of soil and rock masses (compression modulus, Poisson's ratio, cohesion, and friction angle, etc.), which are often difficult to measure directly, especially under large-scale or complex geological conditions. Generally, drilling is used to obtain soil or rock samples to determine these parameters, but due to size effects, the obtained parameters do not necessarily match the actual parameters of the soil and rock mass. Furthermore, soil and rock masses exhibit spatial variability; their physical and mechanical parameters may vary significantly with spatial location, and local measurements alone cannot represent the overall characteristics of the soil and rock mass.

[0010] To obtain accurate physical and mechanical parameters and make numerical simulation results more consistent with reality, researchers typically use parameter inversion methods to adjust the calculation parameters. The specific process is as follows:

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

[0012] However, parameter inversion methods still have the following problems: ① Parameter inversion requires running the numerical model multiple times to iteratively optimize the parameters, and a single calculation of a complex numerical model can take hours or even days, resulting in extremely low inversion efficiency. ② Different parameter combinations may lead to similar model output results, making it impossible to uniquely determine the parameters obtained from the inversion. ③ Noise, sparsity, or insufficient spatiotemporal coverage of the observation data can lead to unstable or biased inversion results.

[0013] Therefore, the above simulation method is still difficult to accurately obtain the ground subsidence in the study area, and the calculation efficiency is low. This 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 assessments. This problem urgently needs to be solved.

[0014] This application provides a method and apparatus for safety assessment of buildings and structures in ground settlement zones, in order to solve the problems in related technologies such as difficulty in accurately obtaining ground settlement in the study area and low computational efficiency, which reduce the efficiency and accuracy of building and structure safety assessment.

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

[0016] Optionally, in one embodiment of this application, constructing a target three-dimensional computational model of the strata and structures of the target study area includes: determining the target study area and establishing a target three-dimensional computational model of the strata and structures of 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 this application, determining the target subsidence rate of each node on the top surface of the stratigraphic model during the target time period using ground-based synthetic aperture radar interferometry monitoring data of the target study area includes: obtaining the surface subsidence value of the target study area during the target time period using the ground-based synthetic aperture radar interferometry monitoring technology; interpolating the surface subsidence amount of each node on the top surface of the stratigraphic model using the surface subsidence value to determine the target subsidence value of each node on the top surface of the stratigraphic model during the target time period, so as to obtain the target subsidence rate based on the target subsidence value.

[0018] Optionally, in one embodiment of this application, the step of applying the target settlement rate of each node to the top surface of the geological model step by step using the target FLAC3D numerical simulation technology to determine the target settlement displacement field of the target three-dimensional computational model includes: determining the target node rate loading step number of the geological model based on the target FLAC3D numerical simulation technology using the surface settlement value of the target study area within the target time period; applying the target settlement rate of each node step by step to the nodes on the top surface of the geological model, wherein no conditions are applied to the nodes inside the geological model; after the geological model has been iterated for one step, the node rate on the top surface of the geological model is cleared to zero, and the surface displacement of the geological model is fixed, and the geological 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 geological model iterates to equilibrium to form the target settlement displacement field of the target three-dimensional computational model.

[0019] Optionally, in one embodiment of this application, generating a building safety assessment result for the target study area using the mechanical response characteristics of the target building and the target safety specifications includes: analyzing the mechanical response characteristics of the target building to determine the deformation and stress characteristics of the target building; and assessing the building safety of the target study area using the deformation characteristics, the stress characteristics, and the target safety specifications to generate the building safety assessment result.

[0020] A second aspect of this application provides a device for assessing the safety of buildings and structures in a land subsidence area, comprising: a construction module for constructing a target three-dimensional computational model of the strata and buildings and structures in a target study area, wherein the target three-dimensional computational model includes a strata model and a building and structure model; an acquisition module for determining the target settlement rate of each node on the top surface of the strata model during the target time period using ground synthetic aperture radar interferometry monitoring data of the target study area during the target time period; a determination module for applying the target settlement rate of each node stepwise to the top surface of the strata model 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 buildings and structures in the building and structure 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 buildings and structures in the target study area using the mechanical response characteristics of the target buildings and the target safety specifications.

[0021] Optionally, in one embodiment of this application, the construction module includes: a determining 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; and a construction unit, used to mesh 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 this application, the acquisition module includes: an acquisition unit, configured to acquire the surface subsidence value of the target study area during the target time period using the ground synthetic aperture radar interferometric monitoring technology; and a first determination unit, configured to interpolate the surface subsidence amount of each node on the top surface of the stratigraphic model using the surface subsidence value to determine the target subsidence value of each node on the top surface of the stratigraphic model during the target time period, so as to acquire the target subsidence rate based on the target subsidence value.

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

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

[0025] A third aspect of this application provides an electronic device, including: 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 assessing the safety of buildings and structures in a land subsidence zone as described in the above embodiments.

[0026] A fourth aspect of this application provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method for assessing the safety of buildings and structures in a land subsidence zone.

[0027] A fifth aspect of this application provides a computer program product, including a computer program that, when executed, is used to implement the above-described method for assessing the safety of buildings and structures in a ground subsidence area.

[0028] This application's embodiments can construct a three-dimensional computational model of the study area, including a geological model and a building / structure model. Using ground-based InSAR (Synthetic Aperture Radar Interferometry) monitoring data of the study area over a certain period, the target settlement rate of each node on the top surface of the geological model is determined. Then, using FLAC3D numerical simulation technology, the target settlement rate is applied step-by-step to the top surface of the geological model to determine the settlement displacement field of the three-dimensional computational model. Furthermore, the mechanical response characteristics of the target buildings / structures in the building / structure model are extracted. Combined with the target safety specifications, the safety assessment results of the buildings / structures in the study area are generated, effectively improving the efficiency and accuracy of building / structure safety assessment. This solves the problems in related technologies, such as the difficulty in accurately obtaining ground settlement in the study area and low computational efficiency, which reduces the efficiency and accuracy of building / structure safety assessment.

[0029] Additional aspects and advantages of this application will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of this application. Attached Figure Description

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

[0031] Figure 1 This is a flowchart of a method for assessing the safety of buildings and structures in a ground settlement zone according to an embodiment of this application;

[0032] Figure 2 This is a schematic diagram of the stratigraphic model and railway model of the study area for a specific embodiment of this application;

[0033] Figure 3 A schematic diagram showing the node rate settings on the top surface of a formation model according to a specific embodiment of this application;

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

[0035] Figure 5 This is a schematic diagram of the settlement field of a 2019 three-dimensional computational model, which is a specific embodiment of this application.

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

[0037] Figure 7 This is a schematic diagram of the settlement deformation of a railway subgrade along its direction in 2023, according to a specific embodiment of this application.

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

[0039] Figure 9 This is a schematic diagram of the uneven settlement curve of the railway subgrade surface along the strike, according to a specific embodiment of this application.

[0040] Figure 10 A flowchart illustrating a method for safety assessment of buildings in land subsidence zones, as a specific embodiment of this application;

[0041] Figure 11 This is a structural schematic diagram of a safety assessment device for buildings in ground settlement zones provided in accordance with an embodiment of this application;

[0042] Figure 12 This is a schematic diagram of the structure of an electronic device provided according to an embodiment of this application. Detailed Implementation

[0043] The embodiments of this application are described in detail below. Examples of these embodiments are shown in the accompanying drawings, wherein the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The embodiments described below with reference to the accompanying drawings are exemplary and intended to explain this application, and should not be construed as limiting this application.

[0044] The following description, with reference to the accompanying drawings, illustrates a method and apparatus for assessing the safety of buildings and structures in a ground settlement zone according to an embodiment of this application. Addressing the problems mentioned in the background art, such as the difficulty in accurately obtaining ground settlement in the study area and low computational efficiency, which reduces the efficiency and accuracy of building and structure safety assessments, this application provides a method for assessing the safety of buildings and structures in a ground settlement zone. In this method, a three-dimensional computational model of the study area, including a geological model and a building / structure model, can be constructed. Using ground-based InSAR monitoring data of the study area over a certain period, the target settlement rate of each node on the top surface of the geological model is determined. Then, using FLAC3D numerical simulation technology, the target settlement rate is applied step-by-step to the top surface of the geological model to determine the settlement displacement field of the three-dimensional computational model. Furthermore, the mechanical response characteristics of the target buildings and structures in the building / structure model are extracted, and combined with the target safety specifications, the safety assessment results of the buildings and structures in the study area are generated, effectively improving the efficiency and accuracy of building and structure safety assessments. Thus, the problems of difficulty in accurately obtaining ground settlement in the study area and low computational efficiency in related technologies, which reduce the efficiency and accuracy of building and structure safety assessments, are solved.

[0045] Specifically, Figure 1 This is a flowchart illustrating a method for assessing the safety of buildings in a ground subsidence zone, as provided in an embodiment of this application.

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

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

[0048] In this embodiment of the application, the target study area is the area where the safety assessment of buildings and structures in the ground subsidence zone is carried out.

[0049] It is understood that the embodiments of this application can construct a three-dimensional computational model of the strata and buildings in the study area. The three-dimensional computational model includes a strata model and a building model. For example, this application can use FLAC3D technology to establish a strata model and a building model, and perform mesh subdivision of the strata model and the building model in the following steps, so as to accurately simulate the interaction between buildings and strata under complex geological conditions and effectively evaluate the stability and safety of buildings.

[0050] In one embodiment of this 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.

[0051] In practical implementation, this application embodiment can select the research scope and establish a three-dimensional calculation model based on engineering geological survey data, including a stratigraphic calculation model and a building / structure calculation model. Furthermore, the stratigraphic calculation model and the building / structure calculation model are meshed. For example, such as... Figure 2 As shown, this application takes the safety evaluation of railway subgrade in a land subsidence area as an example. A geological calculation model and a railway calculation model are established using FLAC3D technology. The geological model has dimensions of 300m*300m*100m and consists of ten layers. The layer thickness, lithology, and calculation parameters are shown in Table 1. Table 1 is a table of geological thickness and mechanical parameters for the calculation model. The specific details of Table 1 are as follows:

[0052] Table 1

[0053]

[0054] The constitutive model for each soil layer is uniformly set to the classic Mohr-Coulomb model, which can describe the stress-strain characteristics of most soil and rock 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 to 13.5m, the bottom surface width to 21.6m, the subgrade thickness to 2.7m, and the slope to 1:1.5. The calculation parameters for the railway subgrade are shown in Table 2, which is a table of physical and mechanical parameters for railway subgrade elements. The specific details of Table 2 are 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 computational model is meshed using triangular prism elements in the following steps. The number of nodes in the meshed railway computational model is 85030 and the number of elements is 158014.

[0059] Furthermore, in order to effectively study the impact of ground settlement on engineering structures, it is first necessary to initialize the geostress field of the three-dimensional calculation model. This application employs a staged elastoplastic solution method to simulate the initial geostress field. During the calculation, the lower boundary of the three-dimensional calculation model is subject to fixed constraints, meaning there will be no displacement in the x, y, and z directions. The front, back, left, and right boundaries of the three-dimensional calculation model are only constrained for normal displacement, while the upper boundary, being the ground surface, is not constrained at this stage. The gravitational acceleration is set to -9.8 m / s². 2 The specific implementation method is as follows:

[0060] First, the constitutive model for the geological material is set as the Mohr-Coulomb model, with its cohesion and tensile strength set at 1000 MPa. Iterative calculations are performed considering only gravity 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 values ​​set for the soil layer, and the calculation is iterated again until the ratio of the maximum unbalanced force to the typical internal force is less than 10. -5 The initialization of the geostress field in the three-dimensional calculation model is complete.

[0061] In addition, to better simulate the ground settlement and displacement field, it is necessary to set the initial pore water pressure field of the three-dimensional calculation model. Based on hydrogeological survey data, the initial water level height of the formation model is set using the water table command, and the water head boundary conditions of the formation model are set so that the pore water pressure field is automatically generated within the formation model according to the linear distribution of the water level height. The unit permeability coefficient is set according to the soil properties of each stratum in the formation model, the mechanical calculation mode of the formation model is turned off, the fluid calculation mode of the formation model is turned on, and equilibrium calculations are performed to the set unbalanced flow ratio, thereby generating the initial seepage field.

[0062] During the simulation, the groundwater level remains constant for a given time period. However, as the simulation progresses to the next time period, the groundwater level needs to be readjusted based on the specific data for that period, thus generating a new seepage field. This method considers changes in the groundwater level, making the evolution of the settlement field in the geological model of this application more consistent with reality.

[0063] In step S102, ground synthetic aperture radar interferometry monitoring data of the target study area during the target time period are used to determine the target subsidence rate of each node on the top surface of the stratigraphic model during the target time period.

[0064] In the embodiments of this application, the target time period can be several years, several months, several weeks, etc., which can be set by relevant technical personnel and are not specifically limited here.

[0065] It is understood that the embodiments of this application can utilize InSAR (Synthetic Aperture Radar Interferometry) technology to obtain ground monitoring data, i.e., surface subsidence values, within a certain time period within the study area. This allows for the determination of the ground subsidence values ​​at various nodes on the top surface of the stratum model within a certain time period within the study area. Furthermore, the settlement rate at each node on the top surface of the stratum model can be obtained based on the ground subsidence values, effectively improving the ability to capture and analyze dynamic changes in the stratum.

[0066] In one embodiment of this application, the target subsidence rate of each node on the top surface of the stratigraphic model during the target time period is determined using ground synthetic aperture radar interferometry monitoring data of the target study area during the target time period. This includes: obtaining the surface subsidence value of the target study area during the target time period using ground synthetic aperture radar interferometry monitoring technology; interpolating the surface subsidence amount of each node on the top surface of the stratigraphic model using the surface subsidence value to determine the target subsidence value of each node on the top surface of the stratigraphic model during the target time period, so as to obtain the target subsidence rate based on the target subsidence value.

[0067] For example, in this embodiment, the surface subsidence values ​​of a certain region in the North China Plain from 2016 to 2023 are used. The Kriging interpolation method is used to interpolate the ground subsidence of each node on the top surface of the stratigraphic model, thereby determining the ground subsidence values ​​of each node on the top surface of the stratigraphic model. This achieves the purpose of refining the surface subsidence displacement field of the stratigraphic model. While ensuring the accuracy of ground subsidence, the simulation time is greatly shortened and the computational complexity 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 represented by the node settlement rate Δ, with the unit being 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 is understood that the embodiments of this application can utilize the target FLAC3D numerical simulation technology to uniformly 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 stratum model. For example, the ground InSAR monitoring data from 2016 to 2023 can be distributed and applied to the surface of the stratum model to determine the settlement displacement field of the three-dimensional calculation model, thereby forming the ground settlement field from 2016 to 2023, effectively ensuring the accuracy and reliability of the calculation.

[0071] In one embodiment of this application, the target settlement rate of each node is applied stepwise to the top surface of the geological model using the target FLAC3D numerical simulation technology to determine the target settlement displacement field of the target three-dimensional computational model. This includes: determining the target node rate loading step number of the geological model based on the target FLAC3D numerical simulation technology using the surface settlement value of the target study area within the target time period; applying the target settlement rate of each node stepwise to the nodes on the top surface of the geological model, wherein no conditions are applied to the nodes inside the geological model; after one iteration of the geological model calculation, the node rate on the top surface of the geological model is cleared to zero, and the surface displacement of the geological model is fixed. The geological model is then iterated multiple times using the target node rate loading step number until all time steps are calculated, so that the geological model iterates to equilibrium to form the target settlement displacement field of the target three-dimensional computational model.

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

[0073] Next, after the previous calculation is completed, the nodal settlement rate Δ2 = Δ / N is applied to the top of the stratum model, and no other conditions are applied to the internal nodes. Then, the stratum model is iterated one more time, and then all the nodal rates of the stratum model are cleared to zero. The displacement of the top surface of the stratum model is fixed again, and the stratum model is iterated to equilibrium, forming a new settlement displacement field of the stratum model.

[0074] Repeating the above iterative process, applying InSAR monitoring data from 2016 to 2023 to the model surface, can produce results such as... Figure 4 , Figure 5 and Figure 6 The ground subsidence field shown is from 2016 to 2023, which can accurately simulate the dynamic changes of the strata over time, ensuring the accuracy and stability of the subsidence displacement field.

[0075] In other words, the embodiments of this application can apply InSAR monitoring data as displacement boundary conditions to the top of the calculation model. Through calculation, the ground settlement field of the study area can be obtained quickly and accurately, avoiding the complex and time-consuming ground settlement simulation process and greatly simplifying the simulation complexity. At the same time, by acquiring InSAR monitoring data from different study areas and time periods, the development patterns of ground settlement in different areas and time periods and their impact on engineering structures can be studied, greatly 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 location in the target three-dimensional calculation model, the mechanical response characteristics of the target buildings in the building model of the target study area are extracted, so as to generate the building safety assessment results of the target study area using the mechanical response characteristics of the target buildings and the target safety specifications.

[0077] In the embodiments of this application, the target building is a building in the ground settlement zone, and the railway subgrade is used as an example in this application.

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

[0079] Optionally, in one embodiment of this application, the safety assessment results of buildings in the target study area are generated by utilizing the mechanical response characteristics of the target building and the target safety specifications. This includes: analyzing the mechanical response characteristics of the target building to determine the deformation and stress characteristics of the target building; and using the deformation and stress characteristics of the target building and the target safety specifications to assess the safety of buildings in the target study area, thereby generating the safety assessment results of buildings.

[0080] For example, embodiments of this application can analyze the mechanical response characteristics of railway structures to obtain, for instance, Figure 7 The settlement and deformation characteristics of the railway structure shown are Figure 8 The stress distribution characteristics are shown. Then, after obtaining the deformation characteristics of the structure, i.e., the railway structure, the safety of the railway structure is evaluated in conjunction with safety specifications. The specific method is as follows:

[0081] Taking ballastless track as an example, my country's post-construction subgrade settlement limit standard for ballastless track clearly specifies the track deformation limit. Table 3 is the standard table of post-construction subgrade settlement limit values ​​for ballastless track, as detailed in Table 3 below:

[0082] Table 3

[0083] Post-construction settlement is generally permissible. Post-construction settlement is allowed for roadbed lengths ≥20m. Uneven settlement Differential Settlement of Misaligned Platforms ≤15mm ≤30mm 20mm / 20m ≤5mm

[0084] If all values ​​are less than the allowable values ​​in Table 3, the railway track is considered basically stable. Specific deformation standards are as follows:

[0085] To reasonably analyze the impact of ground settlement on railway subgrade, this application extracts the Z-direction deformation values ​​at 50 points along the railway subgrade surface and calculates the tilt deformation along the direction of the subgrade. This value can characterize the degree of uneven settlement. The formula for calculating the tilt deformation i is as follows:

[0086]

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

[0088] For example, such as Figure 9 As shown, there is significant uneven settlement on the railway subgrade surface, with the left half of the model representing settlement and the right half representing bulging; [The text abruptly ends here, likely due to an incomplete translation or missing information.] Figure 8 It can be seen that between 2016 and 2022, the uneven settlement of the railway subgrade surface gradually increased over time. However, in 2023, the degree of uneven settlement of the railway decreased significantly, showing a trend of uniform settlement. Between 2016 and 2023, the uneven settlement of the top and bottom surfaces of the railway subgrade was less than 0.8 mm / m, which is still within the limit required by the ballastless track subgrade specifications, indicating that the railway track is basically stable.

[0089] For example, such as Figure 10 As shown, the working principle of the embodiments of this application will be described in detail below with a specific example.

[0090] Step S1001: Establish a three-dimensional calculation model of the strata and buildings. In other words, the embodiments of this application can construct a target three-dimensional calculation model of the strata and buildings in the study area, wherein the target three-dimensional calculation model includes a strata model and a building model.

[0091] Step S1002: Calculate the initial geostress field and pore water pressure field. In other words, the embodiments of this application can set the initial geostress field and initial pore water pressure field of the target three-dimensional calculation model.

[0092] Step S1003: Using ground-based synthetic aperture radar interferometric monitoring data, the settlement rate of the top surface node of the stratum model is obtained. In other words, the embodiments of this application can use ground-based synthetic aperture radar interferometric monitoring data within a certain time period of the study area to determine the settlement rate of the top surface node of the stratum model in the target three-dimensional calculation model.

[0093] Step S1004: Apply the nodal settlement rate stepwise to the surface of the formation model to form a settlement displacement field. In other words, the embodiments of this application can use the target FLAC3D numerical simulation technology to apply the nodal settlement rate stepwise to the surface of the formation model to form a settlement displacement field of the target three-dimensional calculation model.

[0094] Step S1005: Repeat steps S1003-S1004 to obtain the final ground subsidence field. In other words, the embodiments of this application can use the node rate loading step to iterate the stratum model multiple times until all time steps are calculated, so that the stratum model iterates to equilibrium to form the final ground subsidence field of the target three-dimensional calculation model.

[0095] Step S1006: Extract and analyze the mechanical response characteristics of the building structure. In other words, the embodiments of this application can extract the mechanical response characteristics of the building structure in the study area model, such as deformation characteristics and stress characteristics, based on the final settlement field determined by the settlement displacement location of the target three-dimensional calculation model.

[0096] Step S1007: Combine safety standards to evaluate the safety of buildings and structures. In other words, the embodiments of this application can use the deformation characteristics, stress characteristics and safety standards of buildings and structures to evaluate the safety of buildings and structures in the study area, so as to generate building and structure safety evaluation results, effectively improving the efficiency and accuracy of building and structure safety evaluation in the ground settlement area.

[0097] A method for assessing the safety of buildings and structures in a ground subsidence zone, as proposed in this application, can construct a three-dimensional computational model of the study area, including a geological model and a building / structure model. Using ground-based InSAR monitoring data of the study area over a certain period, the target settlement rate of each node on the top surface of the geological model is determined. Then, using FLAC3D numerical simulation technology, the target settlement rate is applied step-by-step to the top surface of the geological model to determine the settlement displacement field of the three-dimensional computational model. Furthermore, the mechanical response characteristics of the target buildings and structures in the building / structure model are extracted. Combined with the target safety standards, the method generates the building / structure safety assessment results for the study area, effectively improving the efficiency and accuracy of building / structure safety assessment. This solves the problems in related technologies, such as the difficulty in accurately obtaining ground subsidence in the study area and low computational efficiency, which reduces the efficiency and accuracy of building / structure safety assessment.

[0098] Next, referring to the accompanying drawings, a safety assessment device for buildings in ground subsidence areas is described according to an embodiment of this application.

[0099] Figure 11 This is a block diagram of a safety assessment device for buildings in ground subsidence areas according to an embodiment of this application.

[0100] like Figure 11 As shown, a safety assessment device 10 for buildings in a ground subsidence area includes: a construction module 100, an acquisition module 200, a determination module 300, and an assessment module 400.

[0101] Specifically, module 100 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 strata 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 stratigraphic model during the target time period by using ground synthetic aperture radar interferometric monitoring data of the target study area during 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 stratum model step by step using the target FLAC3D numerical simulation technology in order to determine the target settlement displacement field of the target three-dimensional calculation model.

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

[0105] Optionally, in one embodiment of this application, the construction module 100 includes a determining unit and a construction unit.

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

[0107] The building blocks are used to mesh the target 3D computational model and set the initial geostress field and initial pore water pressure field of the target 3D computational model.

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

[0109] The acquisition unit is used to acquire the surface subsidence value of the target study area for the target time period using ground-based synthetic aperture radar interferometric monitoring technology.

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

[0111] Optionally, in one embodiment of this application, the determining module 300 includes: a second determining 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 stratigraphic model based on the target FLAC3D numerical simulation technology and the surface subsidence value within the target time period of the target study area.

[0113] An application unit is used to apply the target settlement rate of each node stepwise to the nodes on the top surface of the formation model, wherein no conditions are applied to the nodes inside 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 iteration of the formation model calculation, fix the surface displacement of the formation model, and use the target node rate to load the number of steps to perform multiple iterations on the formation model until all time steps are calculated, so that the formation model iterates to equilibrium to form the target settlement displacement field of the target three-dimensional calculation model.

[0115] Optionally, in one embodiment of this 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 structure in order to determine its deformation and stress characteristics.

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

[0118] It should be noted that the foregoing explanation of an embodiment of a method for assessing the safety of buildings and structures in a land subsidence area also applies to an embodiment of a device for assessing the safety of buildings and structures in a land subsidence area, and will not be repeated here.

[0119] According to an embodiment of this application, a device for assessing the safety of buildings and structures in a ground subsidence zone can construct a three-dimensional computational model of the study area, including a geological model and a building / structure model. Using ground-based InSAR monitoring data of the study area over a certain period, the target settlement rate of each node on the top surface of the geological model is determined. Then, using FLAC3D numerical simulation technology, the target settlement rate is applied step-by-step to the top surface of the geological model to determine the settlement displacement field of the three-dimensional computational model. Furthermore, the mechanical response characteristics of the target buildings and structures in the building / structure model are extracted. Combined with the target safety standards, the device generates the building / structure safety assessment results for the study area, effectively improving the efficiency and accuracy of building / structure safety assessment. This solves the problems in related technologies, such as the difficulty in accurately obtaining ground subsidence in the study area and low computational efficiency, which reduce the efficiency and accuracy of building / structure safety assessment.

[0120] Figure 12 A schematic diagram of the structure of an electronic device provided in an embodiment of this application. The electronic device may include:

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

[0122] When the processor 1202 executes the program, it implements a method for assessing the safety of buildings and structures in the ground subsidence area provided in the above embodiments.

[0123] Furthermore, electronic devices also include:

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

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

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

[0127] If the memory 1201, processor 1202, and communication interface 1203 are implemented independently, then the communication interface 1203, memory 1201, and processor 1202 can be interconnected via a bus to complete communication between them. The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, or an Extended Industry Standard Architecture (EISA) bus, etc. The bus can be divided into address bus, data bus, control bus, etc. For ease of representation, Figure 12 The bus is represented by a single thick line, 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, processor 1202, and communication interface 1203 are integrated on a single chip, then the memory 1201, processor 1202, and 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 this application.

[0130] This embodiment also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method for assessing the safety of buildings and structures in a land subsidence zone.

[0131] This embodiment also provides a computer program product, including a computer program, which, when executed, is used to implement the above-described method for assessing the safety of buildings and structures in a ground settlement zone.

[0132] In the description of this specification, the references to terms such as "one embodiment," "some embodiments," "example," "specific example," or "some examples," etc., indicate that a specific feature, structure, material, or characteristic described in connection with that embodiment or example is included in at least one embodiment or example of this application. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in one or more embodiments or examples. Moreover, without contradiction, those skilled in the art can combine and integrate the different embodiments or examples described in this specification, as well as the features of different embodiments or examples.

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

[0134] Any process or method described in the flowchart or otherwise herein can be understood as representing a module, segment, or portion of code comprising one or N executable instructions for implementing custom logic functions or processes, and the scope of the preferred embodiments of this application includes additional implementations in which functions may be performed not in the order shown or discussed, including substantially simultaneously or in reverse order depending on the functions involved, as should be understood by those skilled in the art to which embodiments of this application pertain.

[0135] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as a ordered list of executable instructions for implementing 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 (such as a computer-based system, a processor-included system, or other system that can fetch and execute instructions from, an instruction execution system, apparatus, or device). For the purposes of this specification, "computer-readable medium" can be any means that can contain, store, communicate, propagate, or transmit programs 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: an electrical connection having one or more wires (electronic device), a portable computer disk drive (magnetic device), random access memory (RAM), read-only memory (ROM), erasable and editable read-only memory (EPROM or flash memory), fiber optic devices, and portable optical disc read-only memory (CDROM). Alternatively, the computer-readable medium may be paper or other suitable media on which the program can be printed, since the program can be obtained electronically by optically scanning the paper or other medium, followed by editing, interpreting, or otherwise processing as necessary, and then stored in a computer memory.

[0136] It should be understood that the various parts of this application can be implemented using hardware, software, firmware, or a combination thereof. In the above embodiments, the N steps or methods can be implemented using software or firmware stored in memory and executed by a suitable instruction execution system. For example, if implemented in hardware, as in another embodiment, it can be implemented using any one or a combination of the following techniques known in the art: discrete logic circuits having logic gates for implementing logical functions on data signals, application-specific integrated circuits (ASICs) having suitable combinational logic gates, programmable gate arrays (PGAs), field-programmable gate arrays (FPGAs), etc.

[0137] Those skilled in the art will understand that all or part of the steps of the methods in the above embodiments can be implemented by a program instructing related hardware. The program can be stored in a computer-readable storage medium, and when executed, the program includes one or a combination of the steps of the method embodiments.

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

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

Claims

1. A method for evaluating safety of a building structure in a land subsidence area, characterized by, Includes the following steps: Construct a target three-dimensional computational model of the stratigraphy and buildings in the target study area, wherein the target three-dimensional computational model includes a stratigraphic model and a building model; Using ground-based synthetic aperture radar interferometric monitoring data of the target study area during the target time period, the target subsidence rate of each node on the top surface of the stratigraphic model during the target time period is determined. Using the target FLAC3D numerical simulation technology, the target settlement rate of each node is applied stepwise to the top surface of the formation model to determine the target settlement displacement field of the target three-dimensional computational model. This step of applying the target settlement rate of each node stepwise to the top surface of the formation model to determine the target settlement displacement field of the target three-dimensional computational model includes: Based on the target FLAC3D numerical simulation technology, the target node rate loading steps of the stratigraphic model are determined using the surface subsidence value of the target study area within the target time period. The target subsidence rate of each node is applied step by step to the top node of the stratigraphic model, wherein no conditions are applied to the nodes inside the stratigraphic model. After one iteration of the stratigraphic model, the node rate of the top surface of the stratigraphic model is cleared to zero, and the surface displacement of the stratigraphic model is fixed. The stratigraphic model is iterated multiple times using the target node rate loading steps until the calculation of all time steps is completed, so that the stratigraphic model iterates to equilibrium to form the target subsidence displacement field of the target three-dimensional calculation model. Based on the final settlement field determined by the target settlement displacement location in the target three-dimensional calculation model, the mechanical response characteristics of the target buildings in the building model of the target study area are extracted, so as to generate the building safety assessment results of the target study area using the mechanical response characteristics of the target buildings and the target safety specifications.

2. The method of claim 1, wherein, The target three-dimensional computational model of the stratigraphy and structures in the target study area includes: The target study area is determined, and a target three-dimensional computational model of the strata and structures in the target study area is established. The target three-dimensional computational model is meshed, and the initial geostress field and initial pore water pressure field of the target three-dimensional computational model are set.

3. The method of claim 1, wherein, The step of determining the target subsidence rate of each node on the top surface of the stratigraphic model during the target time period using ground-based synthetic aperture radar interferometric monitoring data of the target study area includes: The ground-based synthetic aperture radar interferometric monitoring technology is used to obtain the surface subsidence value of the target study area during the target time period; The ground settlement of each node on the top surface of the stratigraphic model is interpolated using the surface settlement value to determine the target settlement value of each node on the top surface of the stratigraphic model for the target time period, so as to obtain the target settlement rate based on the target settlement value.

4. The method of claim 1, wherein, The process of generating safety assessment results for buildings and structures in the target study area using the mechanical response characteristics of the target buildings and structures and the target safety standards includes: Analyze the mechanical response characteristics of the target building to determine its deformation and stress characteristics; The safety of buildings in the target study area is assessed using the deformation characteristics, stress characteristics, and target safety standards of the target buildings to generate the safety assessment results.

5. A ground subsidence zone construction structure safety evaluation device characterized by comprising: 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 strata model and a building model; The acquisition module is used to determine the target subsidence rate of each node on the top surface of the stratigraphic model during the target time period using ground synthetic aperture radar interferometric monitoring data of the target study area during the target time period. The determination module is used to apply the target settlement rate of each node stepwise 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. The step of applying the target settlement rate of each node stepwise 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 includes: Based on the target FLAC3D numerical simulation technology, the target node rate loading steps of the stratigraphic model are determined using the surface subsidence value of the target study area within the target time period. The target subsidence rate of each node is applied step by step to the top node of the stratigraphic model, wherein no conditions are applied to the nodes inside the stratigraphic model. After one iteration of the stratigraphic model, the node rate of the top surface of the stratigraphic model is cleared to zero, and the surface displacement of the stratigraphic model is fixed. The stratigraphic model is iterated multiple times using the target node rate loading steps until the calculation of all time steps is completed, so that the stratigraphic model iterates to equilibrium to form the target subsidence displacement field of the target three-dimensional calculation model. The evaluation module 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 location of the target three-dimensional calculation model, so as to generate the building safety evaluation results of the target study area using the mechanical response characteristics of the target buildings and the target safety specifications.

6. The apparatus of claim 5, wherein, The building module includes: A 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; The construction unit is used to mesh 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.

7. The apparatus of claim 5, wherein, The acquisition module includes: The acquisition unit is used to acquire the surface subsidence value of the target study area during the target time period using the ground synthetic aperture radar interferometric monitoring technology; The first determining unit is used to interpolate the ground settlement of each node on the top surface of the stratum model using the surface settlement value, so as to determine the target settlement value of each node on the top surface of the stratum model during the target time period, and to obtain the target settlement rate based on the target settlement value.

8. The apparatus of claim 5, wherein, The evaluation module includes: The analysis unit is used to analyze the mechanical response characteristics of the target structure in order to determine the deformation and stress characteristics of the target structure. An evaluation unit is used to evaluate the safety of buildings in the target study area using the deformation characteristics, stress characteristics, and target safety standards of the target buildings, so as to generate the safety evaluation results of the buildings.

9. An electronic device, comprising: include: The system includes 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 assessing the safety of buildings and structures in a land subsidence zone as described in any one of claims 1-4.

10. A computer-readable storage medium having stored thereon a computer program, characterized in that, The program is executed by the processor to implement a method for assessing the safety of buildings and structures in a ground subsidence zone as described in any one of claims 1-4.

11. A computer program product comprising a computer program, characterized in that, The computer program is executed by a processor to implement a method for assessing the safety of buildings and structures in a land subsidence zone as described in any one of claims 1-4.

Images