A method for fast positioning and synthesizing three-dimensional numerical model based on two-dimensional profile topography

Abstract

The application relates to a kind of based on two-dimensional profile topography fast positioning synthesis three-dimensional numerical model method, it is related to geotechnical engineering numerical simulation technical field, including obtaining the two-dimensional profile of the region to be modeled, and pretreatment is carried out;From three-dimensional terrain, the plane position where each two-dimensional profile is located is determined, the azimuth angle and center point coordinate of two-dimensional profile are calculated;Extract three-dimensional terrain point cloud, then two-dimensional profile is processed one by one, to rotate and translate with rotation amount and best translation amount, synthesis three-dimensional stratum distribution model: for the internal stratum of three-dimensional stratum distribution model, increase stratum control point, improve the precision of three-dimensional stratum distribution model using distance reciprocal interpolation method, and using FLAC3D establishes three-dimensional numerical model;It can provide theoretical basis for slope excavation measure, disaster assessment determination, for complex working condition slope stability analysis and deformation mechanism.

Description

Technical Field

[0001] This invention relates to the field of numerical simulation technology in geotechnical engineering, specifically to a method for rapidly locating and synthesizing a three-dimensional numerical model based on two-dimensional profile terrain. Background Technology

[0002] Slope stability is a crucial aspect of mechanical analysis in hydraulic and geotechnical engineering. It serves as a vital basis for analyzing the stress characteristics of slopes under conditions such as earthquakes, rainfall, water level fluctuations, and excavation, exploring mechanical evolution, and formulating support schemes. However, slope stability analysis results rely heavily on a good numerical model. A good model often yields ideal results, while a model with poor mesh quality not only performs poorly under the same conditions but also has difficulty converging. Currently, engineering design units in my country often only provide three-dimensional slope topography, on-site slip surface outcrops, a few stratigraphic boreholes, or a few computational two-dimensional profiles when designing slopes. Establishing a refined three-dimensional numerical model requires detailed stratigraphic interfaces, and the more stratigraphic control points, the higher the stratigraphic control accuracy. This process is usually manually operated by modelers, who subjectively process stratigraphic information near the boundaries and then use bidirectional sweeping to form a stratigraphic model. This process is time-consuming, labor-intensive, highly subjective, and the resulting model accuracy is not high. Therefore, there is an urgent need for a method for rapidly locating and synthesizing three-dimensional numerical models based on two-dimensional profile topography to solve these problems. Summary of the Invention

[0003] The purpose of this invention is to provide a method for rapidly locating and synthesizing a three-dimensional numerical model based on two-dimensional profile terrain to solve the problems mentioned in the background art.

[0004] To solve the above-mentioned technical problems, the present invention adopts the following technical solution: a method for rapidly locating and synthesizing a three-dimensional numerical model based on two-dimensional profile terrain, comprising the following steps:

[0005] Obtain the two-dimensional profile of the region to be modeled and perform preprocessing;

[0006] Determine the planar position of each two-dimensional profile from the three-dimensional terrain, and calculate the azimuth and center point coordinates of the two-dimensional profile;

[0007] Extract the 3D terrain point cloud, and then process the 2D profiles one by one, including:

[0008] The azimuth angle of the corresponding two-dimensional profile is used as the rotation amount;

[0009] The optimal translation amount is determined by sampling along the two-dimensional profile line and minimizing the terrain elevation residual.

[0010] All two-dimensional profiles are rotated and translated with rotation and optimal translation to synthesize a three-dimensional stratigraphic distribution model.

[0011] Interpolation is performed on the internal strata of the three-dimensional stratigraphic distribution model, stratigraphic control points are added, and the accuracy of the three-dimensional stratigraphic distribution model is improved by using the inverse distance interpolation method. Finally, a three-dimensional numerical model is established using FLAC3D.

[0012] Preferably, the preprocessing includes: translating the two-dimensional profile to the AutoCAD origin while keeping the vertical elevation unchanged.

[0013] Preferably, the preprocessing further includes: translating the two-dimensional profile to make the elevation of typical feature locations consistent with the three-dimensional terrain, assigning the two-dimensional terrain lines to layer 0, and numbering the strata in order from top to bottom.

[0014] Preferably, the coordinates of two points (x1, y1) and (x2, y2) in the horizontal direction of each two-dimensional profile are obtained;

[0015] For the i-th two-dimensional profile, the azimuth angle alfai between the two-dimensional profile and the positive X-axis is calculated using the coordinates of points (x1, y1) and (x2, y2):

[0016] ax = x2 - x1;

[0017] ay = y2 - y1;

[0018] ;

[0019] Obtain the unit direction vector:

[0020] ax = ax / d;

[0021] ay = ay / d;

[0022] If ay is greater than 0, then the azimuth angle alfai = acos(ax) * 360 / 3.1415926;

[0023] If ay is less than 0, then the azimuth angle alfai = 360.0 - acos(ax) * 360 / 3.1415926;

[0024] The coordinates of the center point (xc0, yc0) are:

[0025] xc0 = (x1 + x2) / 2.0;

[0026] yc0=(y1+y2) / 2.0.

[0027] Preferably, a 3D terrain point cloud is extracted from an existing 3D terrain model as a reference for the real terrain; wherein processing each 2D profile specifically involves:

[0028] For the i-th two-dimensional profile, extract the points on its two-dimensional terrain lines;

[0029] Calculate the direction vector (ax, ay) of the two-dimensional profile, and take the first point (x1, y1) as the initial base point;

[0030] Starting from the initial base point, a candidate translation position (xk, yk) is generated along the azimuth direction of the two-dimensional profile line at a preset step size;

[0031] Translate all points on the 2D terrain line to the candidate position (xk, yk) and rotate them by an angle of alfai.

[0032] The true elevation is obtained by interpolating the planar coordinates of each translated and rotated point in the 3D terrain point cloud. Then, the difference between the elevation and the 2D profile elevation is calculated, squared, and summed to obtain the residual sum of squares for the candidate translation position.

[0033] Iterate through all candidate translation positions and record the position (xk_best, yk_best) that minimizes the sum of squared residuals.

[0034] The two-dimensional profile is translated with (xk_best, yk_best) as the base point and rotated by the angle alfai.

[0035] Preferably, the reciprocal distance interpolation method includes:

[0036] Suppose we want to find the spatial coordinates of (xi, yi) on the plane corresponding to the stratum on the ground. With (xi, yi) as the local coordinate base point, we divide the plane into ii intervals with the positive X-axis as 0° and the counterclockwise direction of the X-axis as 0°. In each interval, we search for a point that is closest to (xi, yi) and count its distance as dk (k=1,...,8).

[0037] Stop when there is a point (xm, ym, zm) whose horizontal distance is less than 1m, and zi directly takes the zm value of that point;

[0038] If no value exists with a distance less than 1m, then the value of zi is obtained using the sum of points zk and the weight (1.0 / dk), as shown in the following formula:

[0039] .

[0040] Preferably, the ground surface data encrypted using the inverse distance weighting method is exported, the portion outside the landslide boundary line is removed, and the contour lines within the stratum are retained as control stratum lines.

[0041] Preferably, a three-dimensional numerical model is established using three-dimensional software and exported as a macroscopic numerical model usable by FLAC3D; and PFC3D6.0 is used to transform the potential sliding body into a microscopic granular flow model for landslide mechanism simulation.

[0042] Beneficial effects: This invention utilizes the acquired two-dimensional profile, based on the two-dimensional topographic lines on the profile as the control point set and the profile position, to restore the two-dimensional profile to a spatial model through optimization search. It then densifies the number of stratigraphic control points to match the three-dimensional topography, thus facilitating the creation of a refined numerical model using commercial software for stability analysis. The method further transforms the region of interest into a micro-model for landslide mechanism and disaster analysis, providing a theoretical basis for slope excavation measures and disaster assessment. Furthermore, the numerical model generated using this method is fast, significantly improving research efficiency. Simultaneously, the micro-model employs a particle system servo method controlled by a measurement circle, reducing the problem of long computation time caused by large models and numerous internal particles, resulting in more realistic and reliable simulation results. Attached Figure Description

[0043] Figure 1 This is a two-dimensional cross-sectional view provided by the design institute in this invention;

[0044] Figure 2 This is a diagram showing the location of the three-dimensional terrain and two-dimensional profile that need to be adjusted in this invention;

[0045] Figure 3 This is a distribution map of the special stratigraphic control boundaries in this invention;

[0046] Figure 4 This is a typical two-dimensional stratigraphic control profile in this invention;

[0047] Figure 5 This is a schematic diagram illustrating the matching of two-dimensional profile terrain to a three-dimensional model in this invention;

[0048] Figure 6 This is a schematic diagram showing all two-dimensional cross-sections merged onto a three-dimensional terrain in this invention;

[0049] Figure 7 This is a schematic diagram of the bottom sliding surface synthesized using two-dimensional cross-sections in this invention;

[0050] Figure 8 This is a schematic diagram of the bottom sliding surface layer after interpolation in the present invention;

[0051] Figure 9 This is a schematic diagram of a continuous three-dimensional slope numerical model established using software in this invention;

[0052] Figure 10 This is a graph showing the calculation results of the continuous numerical model in this invention;

[0053] Figure 11 This is a schematic diagram of the continuous-discontinuous model formed based on the continuous model in this invention;

[0054] Figure 12 The figure shows the calculation results of the continuous-discontinuous numerical model in this invention. Detailed Implementation

[0055] To make the objectives and advantages of this invention clearer, the invention will be specifically described below with reference to embodiments. It should be understood that the following text is merely used to describe one or more specific embodiments of the invention and does not strictly limit the scope of protection specifically claimed by the invention.

[0056] A massive landslide deposit, with a volume approaching 50 million cubic meters, was reactivated by a sudden torrential rain. It only stabilized after a certain period following the rain, posing a key challenge for local hydropower construction. This necessitates systematic engineering geomechanical analysis, the establishment of a relatively continuous numerical model to evaluate slope stability, and the development of a refined continuous-discontinuous numerical model for landslide mechanism analysis. This invention employs a method based on rapid localization and synthesis of a three-dimensional numerical model using two-dimensional topographic profiles, as detailed below:

[0057] S1. Extract a section from the 3D topographic map of the area of ​​interest provided by the design institute, select the modeling range, and retain the 2D profile position lines and 3D topographic contour lines. The modeling uses a rectangular area parallel to the XOY coordinate axes. All elements to be retained are translated to near X=0.0, Y=0.0, and appropriately rotated to facilitate subsequent mechanical analysis. In this embodiment, the top view of the area of ​​interest provided by the design institute is shown below. Figure 1 The image shown is an AutoCAD top view of the three-dimensional terrain and two-dimensional profile location near the slope. It can be seen that the image contains a large amount of graphic element information. The top view is cleaned up by deleting unnecessary information and retaining the necessary information (three-dimensional terrain, two-dimensional profile location, and sliding body boundary line).

[0058] S11. For files in the geodetic coordinate system (with the Qingdao National Benchmark Point as the origin), ensure that the vertical (Z-axis) elevation remains unchanged and translate it to the vicinity of the AutoCAD origin to ensure that the mathematical operations on the coordinates in the numerical calculation do not exceed the upper and lower limits (-10^33~10^33) for real number types when calculated by commercial software.

[0059] S12. Design the transverse river direction as the positive X-axis of the modeling, select the modeling range, and based on the provided plan and the designed model placement direction, rotate all elements 90° around the AutoCAD origin, using the transverse river direction as the X-axis and the longitudinal river direction as the Y-axis. Delete all elements outside the modeling range; Figure 2 The coordinates of the lower left corner of the model after translation and rotation are xmin=0.0, ymin=0.0; the distribution map of the control boundary of special strata is as follows. Figure 3 As shown.

[0060] S13. Save the 3D terrain lines as a separate file, such as "3D Terrain Lines.dxf". Ensure that the file format is an R12 version dxf file (the minimum version of AutoCAD) so that it can be opened by different versions of AutoCAD.

[0061] S14. Based on the number of two-dimensional section position lines in the top view, find the coordinates of the two-dimensional section control points. Each two-dimensional section position line consists of two points, let these two points be (x1, y1) and (x2, y2).

[0062] S2. Translate each 2D profile (XOY coordinates) to ensure that the elevations of key points on the 2D model are consistent with those on the 3D model. Retain necessary topographic lines and stratigraphic lines that need to be considered, delete irrelevant primitives, and horizontally shift all primitives to around X=0.0. Set the topographic lines as layer 0, and represent other stratigraphic layers numerically from top to bottom. Specifically:

[0063] S21. Store each 2D profile as a separate file, and translate it according to typical feature positions on the 2D profile (such as excavation elevation, borehole elevation, elevation markings, etc.) to make the elevation consistent with the actual elevation. The horizontal 2D profiles should be located as close as possible to the AutoCAD origin. A total of 17 typical 2D profiles were provided for this landslide, such as... Figure 4 As shown.

[0064] S22. For each two-dimensional profile, the topographic lines of the two-dimensional profile are placed in layer 0. The stratigraphic lines to be considered are in layers 1, 2, 3, etc. The larger the number, the lower the stratigraphic level.

[0065] S23. Each two-dimensional profile is processed and then saved as a separate dxf file in R12 format.

[0066] S3. Extract the position of each two-dimensional profile from the three-dimensional terrain. Each two-dimensional profile only needs the coordinates of two points (x1, y1) and (x2, y2) in the horizontal direction. Calculate the angle alfai between these coordinates and the X-axis as the direction angle, and record the position of the center point (xc0, yc0) of the i-th two-dimensional profile. Specifically:

[0067] S31. For the i-th two-dimensional profile, calculate the angle (counterclockwise) between the two-dimensional profile and the positive X-axis using the two-dimensional profile control points (x1, y1) and (x2, y2), as follows:

[0068] ax = x2 - x1;

[0069] ay = y2 - y1;

[0070] ;

[0071] Obtain the unit direction vector:

[0072] ax = ax / d;

[0073] ay = ay / d;

[0074] If ay is greater than 0.0, then: alfai = acos(ax) * 360 / 3.1415926;

[0075] If ay is less than 0.0, then: alfai = 360.0 - acos(ax) * 360 / 3.1415926;

[0076] Step S32: Calculate the approximate center point coordinates (xc0, yc0) of the two-dimensional profile:

[0077] xc0 = (x1 + x2) / 2.0;

[0078] yc0=(y1+y2) / 2.0;

[0079] S4. Extract the 3D terrain point cloud, with NPP points, and store it in the value xyz0(NPP,3). Then, traverse N 2D profiles. For the i-th 2D profile, extract the terrain line (layer 0), with Ni0 points. Rotate the 2D profile line by an angle alfai, and calculate the average coordinate of the 2D terrain line, set as (xc1, yc1). Along the 2D profile line, calculate the interpolation of the terrain point (x, y) on the 3D terrain in steps of 1m. Calculate the coordinates of each point relative to the 2D terrain, take the sum of squared residuals of all terrain points as the consideration, and record the position of the planar model (xc2, yc2) when the residual is minimized. Rotate all strata on the 2D model by alfai, translate xc2 in the X direction, and yc2 in the Y direction. After traversing all 2D profiles, form a synthetic 3D strata distribution model. Specifically, this includes:

[0080] S41. Use programming languages ​​such as Fortran and Python to extract the coordinates of points on a 3D terrain line. Assume there are NPP points, which are stored in the array xyz0(NPP,3). The array NPP has 3 rows and 3 columns, where columns 1, 2, and 3 are the X, Y, and Z coordinates of the corresponding terrain points, respectively.

[0081] S42. Loop through N two-dimensional profiles; without loss of generality, for the i-th two-dimensional profile, extract its topographic lines (layer 0), and let the number of points be Ni0, store them in the array xy0(Ni0,2); traverse the Ni0 points, calculate the average value of its X coordinate xp, and set the average value of its Y coordinate yp to 0.0 (no need to calculate).

[0082] ; ;

[0083] Taking the first point of the two-dimensional terrain profile as the initial point, and then moving along the (ax, ay) direction, with each step being 1m, the total number of attempts is NUM = int(d), where int is the floor function; for any point in the total NUM attempts, let it be i, then the attempted position is:

[0084] (xk,yk)=(x1+(i-1)*ax,y1*(i-1)*ay);

[0085] Using (xp,yp) as the base point, translate the i-th two-dimensional profile topographic point to the (xk,yk) point;

[0086] That is, points on a two-dimensional terrain line are translated one by one using the following formula:

[0087] xy0(:,1)=xyz(:,1)+(xk-xp);

[0088] xy0(:,2)=xyz(:,2)+(yk-yp);

[0089] xy0(:,3)=xyz(:,3);

[0090] The colon (:) in the formula indicates line-by-line processing.

[0091] Using (xk, yk) as the base point, rotate counterclockwise by 12 degrees (which needs to be converted to radians) for Ni0 terrain points, where the formula for the i-th point is as follows:

[0092] xy0(ii,1)=xy0(ii,1)*cos(alfai)-xyz2(ii,2)*sin(alfai)+xc;

[0093] xy0(ii,2)=xy0(ii,2)*sin(alfai)+xyz2(ii,2)*cos(alfai)+yc;

[0094] xy0(ii,3)=xy0(ii,3);

[0095] In the formula, ii is the row number corresponding to the point, 1 is the X coordinate, and 2 is the Y coordinate.

[0096] Let the initial value of the sum of squared residuals of terrain be error = 10000000; loop through Ni0 terrain points, cancha = 0.0, arbitrarily select the i-th point, take its X and Y coordinates, interpolate to calculate the corresponding three-dimensional terrain value, denoted as Z2, and the terrain elevation translated from it is denoted as Z1 = xy0(ii,3), then cancha = cancha + (z1-z2)*(z1-z2), after traversing Ni0 terrain points, cancha is obtained; if the value of cancha is smaller than error, then error = cancha.

[0097] After traversing NUM positions and trying them, we will find the position with the smallest residual value (xk_best, yk_best);

[0098] S43. Using the above formula, translate and rotate all topographic lines, bottom lines, and other primitives in the two-dimensional profile by an angle of (xk_best, yk_best) to a 3D model; the two-dimensional profile terrain is matched to the three-dimensional model as follows. Figure 5 As shown, all two-dimensional profiles are merged into the three-dimensional terrain, as follows: Figure 6 As shown.

[0099] S44. Keep the X and Y coordinates of the points on the landslide boundary line unchanged, perform terrain interpolation on the Z coordinate, and replace it to prevent the boundary from being affected by errors.

[0100] S5. For the internal strata of the 3D stratigraphic distribution model, interpolation is performed by adding stratigraphic control points and using the inverse distance interpolation method to improve the accuracy of the 3D stratigraphic distribution model. Specifically, this includes:

[0101] S51. Using the closed landslide boundary line as the boundary, delete the topography or other strata within the boundary, retaining only the strata that need to be considered. Use the Kriging method in Surfer 8 for strata interpolation or the inverse distance weighting method for interpolation. The inverse distance weighting method includes:

[0102] Suppose we want to find the spatial coordinates of (xi, yi) on the plane corresponding to the stratum on the ground. With (xi, yi) as the local coordinate base point, and with the positive X-axis as 0°, we divide the space into ii (actually 8) intervals in a counterclockwise direction. In each interval, we search for a point that is closest to (xi, yi) and count its distance as dk (k=1,...,8).

[0103] Stop when there is a point (xm, ym, zm) whose horizontal distance is less than 1m, and zi directly takes the zm value of that store;

[0104] If no value exists with a distance less than 1m, the value of zi is obtained using the sum of points zk and the weight (1.0 / dk), as shown in the following formula:

[0105] ;

[0106] S52. Export the ground stratum data encrypted using Surfer interpolation or the inverse distance weighting method, remove the portion outside the slip body boundary line, and retain the contour lines within the stratum as control stratum lines. The two-dimensional profile of the bottom slip surface is synthesized using the two-dimensional profile as shown below. Figure 7 As shown, the bottom slip layer after intercalation is as follows: Figure 8 As shown.

[0107] S6. Use 3D software to create a 3D numerical model;

[0108] In one embodiment, a rhinoceros sieve is used to extrude the established strata to form a slope model; a mesh is generated and exported as a finite difference model for deformation stability analysis; key components such as the landslide body are transformed into a microscopic model for landslide mechanism analysis and disaster prediction analysis; the process includes the following steps:

[0109] S61. By setting separate layers for terrain and strata, and then using commercial software such as Rhino to manually process areas where the 2D profile and boundaries are not precise enough.

[0110] S62. Utilize Rhino's extrusion and other functions to build a three-dimensional numerical model and export it as a macroscopic numerical model usable in FLAC3D. A continuous three-dimensional slope numerical model built using the software is shown below. Figure 9 As shown in Figure 10, the calculation results of the continuous numerical model are as follows.

[0111] S63. Using PFC3D6.0, the potential sliding body is transformed into a micro-granular flow model to simulate the landslide mechanism.

[0112] Continuous-discontinuous models formed based on continuous models, such as Figure 11 As shown, the calculation results of the continuous-discontinuous numerical model Figure 12 As shown.

[0113] The embodiments of the present invention have been described in detail above with reference to the examples. However, the present invention is not limited to the above embodiments. For those skilled in the art, after learning the contents described in the present invention, several equivalent changes and substitutions can be made without departing from the principle of the present invention. These equivalent changes and substitutions should also be considered to fall within the protection scope of the present invention.

Claims

1. A method for rapidly locating and synthesizing a three-dimensional numerical model based on two-dimensional terrain profiles, characterized in that: Includes the following steps: Obtain a two-dimensional profile of the area to be modeled and perform preprocessing, including: translating the two-dimensional profile to the AutoCAD origin while keeping the vertical elevation unchanged; translating the two-dimensional profile so that the elevation of typical feature locations is consistent with the three-dimensional terrain; assigning the two-dimensional terrain lines to layer 0; and numbering the strata in order from top to bottom. Determine the planar position of each two-dimensional profile from the three-dimensional terrain, and calculate the azimuth and center point coordinates of the two-dimensional profile; Extract the 3D terrain point cloud, and then process the 2D profiles one by one, including: The azimuth angle of the corresponding two-dimensional profile is used as the rotation amount; The optimal translation amount is determined by sampling along the two-dimensional profile line and minimizing the terrain elevation residual. All two-dimensional profiles are rotated and translated with rotation and optimal translation to synthesize a three-dimensional stratigraphic distribution model. Interpolation is performed on the internal strata of the three-dimensional stratigraphic distribution model, stratigraphic control points are added, and the accuracy of the three-dimensional stratigraphic distribution model is improved by using the inverse distance interpolation method. Finally, a three-dimensional numerical model is established using FLAC3D.

2. The method for rapidly locating and synthesizing a three-dimensional numerical model based on two-dimensional profile terrain according to claim 1, characterized in that: Obtain the coordinates (x1, y1) and (x2, y2) of two points in the horizontal direction of each two-dimensional profile; For the i-th two-dimensional profile, the azimuth angle alfai between the two-dimensional profile and the positive X-axis is calculated using the coordinates of points (x1, y1) and (x2, y2): ax = x2 - x1; ay = y2 - y1; ； Obtain the unit direction vector: ax = ax / d; ay = ay / d; If ay is greater than 0, then the azimuth angle alfai = acos(ax) * 360 / 3.1415926; If ay is less than 0, then the azimuth angle alfai = 360.0 - acos(ax) * 360 / 3.1415926; The coordinates of the center point (xc0, yc0) are: xc0 = (x1 + x2) / 2.0; yc0=(y1+y2) / 2.

0.

3. The method for rapidly locating and synthesizing a three-dimensional numerical model based on two-dimensional profile terrain according to claim 2, characterized in that: Extract 3D terrain point clouds from existing 3D terrain models as a reference for real terrain; the specific processing of each 2D profile involves: For the i-th two-dimensional profile, extract the points on its two-dimensional terrain lines; Calculate the direction vector (ax, ay) of the two-dimensional profile, and take the first point (x1, y1) as the initial base point; Starting from the initial base point, a candidate translation position (xk, yk) is generated along the azimuth direction of the two-dimensional profile line at a preset step size; Translate all points on the 2D terrain line to the candidate position (xk, yk) and rotate them by an angle of alfai. The true elevation is obtained by interpolating the planar coordinates of each translated and rotated point in the 3D terrain point cloud. Then, the difference between the elevation and the 2D profile elevation is calculated, squared, and summed to obtain the residual sum of squares for the candidate translation position. Iterate through all candidate translation positions and record the position (xk_best, yk_best) that minimizes the sum of squared residuals. The two-dimensional profile is translated with (xk_best, yk_best) as the base point and rotated by the angle alfai.

4. The method for rapidly locating and synthesizing a three-dimensional numerical model based on two-dimensional profile terrain according to claim 3, characterized in that: The reciprocal distance interpolation method includes: Suppose we want to find the spatial coordinates of (xi, yi) on the plane corresponding to the stratum on the ground. With (xi, yi) as the local coordinate base point, we divide the plane into ii intervals with the positive X-axis as 0° and the counterclockwise direction of the X-axis as 0°. In each interval, we search for a point that is closest to (xi, yi) and count its distance as dk (k=1,...,8). Stop when there is a point (xm, ym, zm) whose horizontal distance is less than 1m, and zi directly takes the zm value of that point; If no value exists with a distance less than 1m, then the value of zi is obtained using the elevation zk of points ii and the weight (1.0 / dk), as shown in the following formula: 。 5. The method for rapidly locating and synthesizing a three-dimensional numerical model based on two-dimensional profile terrain according to claim 4, characterized in that: Export the ground stratum data encrypted using the inverse distance weighting method, remove the portion outside the landslide boundary line, and retain the contour lines within the stratum as control stratum lines.

6. The method for rapidly locating and synthesizing a three-dimensional numerical model based on two-dimensional terrain profiles according to claim 1, characterized in that: A three-dimensional numerical model was established using 3D software and exported as a macroscopic numerical model usable in FLAC3D. Then, using PFC3D6.0, the potential sliding body was transformed into a microscopic granular flow model to simulate the landslide mechanism.

Images