Three-dimensional geological model adaptive correction method based on settlement monitoring data

Abstract

The invention discloses a three-dimensional geological model adaptive correction method based on settlement monitoring data. The method comprises the following steps: S1, acquiring and sorting site geological survey data; s2, plane information of terrain points and geological drill holes and a boundary range used for later calculation of settlement are obtained, Denaulay triangulation is carried out in the range, and interpolation is carried out to form a plane; s3, constructing a three-dimensional geologic model; s4, current-period settlement monitoring data of all ground surface settlement monitoring points in the field settlement area range are obtained and sorted; s5, calculating ground surface settlement and generating a settlement cloud picture; and S6, based on the stratum distribution condition in the three-dimensional geologic model, the compression modulus of the lithology of each soil layer and the depth of the soil layer, calculating the settlement amount of each node on the bottom surface of each stratum. According to the method, the three-dimensional geological model is dynamically optimized through real-time data feedback of the ground surface settlement monitoring data, so that more accurate geological analysis and prediction are expected to be achieved.

Description

Technical Field

[0001] The invention relates to the technical field of construction design, and in particular to a three-dimensional geological model adaptive correction method based on settlement monitoring data. Background Art

[0002] Settlement monitoring in geological engineering is an important part of ensuring structural safety and environmental stability. By monitoring the vertical displacement of the surface or structures, the changes in geological conditions and the stability of engineering structures can be evaluated. The three-dimensional geological model is a method of digitally expressing underground geological structures. It can display the spatial distribution, morphological characteristics and their interrelationships of geological bodies in a three-dimensional visual form. During the construction and design phase of an engineering project, engineering geological information is often obtained through drilling and other means. The layered summation method is generally used to analyze and calculate surface settlement based on indoor test parameters or in-situ test parameters of geotechnical surveys. On the one hand, the compression modulus of each stratum is generally obtained by indoor tests, which is very different from the compression modulus of geotechnical in-situ tests. At the same time, the spatial variability characteristics of geotechnical parameters are not considered. On the other hand, the three-dimensional geological model generated based on geological drilling data does not take into account the treatment processes such as foundation reinforcement during the construction phase, and cannot achieve dynamic adjustment of the three-dimensional geological model during the construction phase. During the construction phase, project personnel still need to rely on their own engineering experience and subsequent settlement monitoring to judge the settlement of the site.

[0003] Although the compression modulus of geotechnical materials is very different from the in-situ parameters of geotechnical materials, it is generally believed that the ratio of the compression modulus between geotechnical materials is equivalent to the ratio of the in-situ compression modulus of geotechnical materials. By using this characteristic of geotechnical material parameters and combining the on-site monitoring settlement data, the settlement of each point on the bottom of each stratum can be estimated, and the adaptive correction of the three-dimensional geological model can be achieved based on the drilling data in the early stage of the project.

[0004] Therefore, a method for adaptively correcting a three-dimensional geological model based on settlement monitoring data is provided. Summary of the invention

[0005] The purpose of the present invention is to overcome the existing defects and provide a three-dimensional geological model adaptive correction method based on settlement monitoring data. Through real-time data feedback of surface settlement monitoring data, the three-dimensional geological model is dynamically optimized to achieve more accurate geological analysis and prediction.

[0006] The technical solution to achieve the above purpose is:

[0007] A three-dimensional geological model adaptive correction method based on settlement monitoring data, comprising:

[0008] Step S1, obtaining and collating site geological survey data;

[0009] Step S2, obtaining the plane information of the topographic points and geological boreholes and the boundary range used for the later calculation of the settlement, performing Denaulay triangulation within the range, and interpolating it into a surface;

[0010] Step S3, constructing a three-dimensional geological model;

[0011] Step S4, obtaining and collating the current settlement monitoring data of all surface settlement monitoring points within the settlement area of ​​the site;

[0012] Step S5, calculating the surface settlement and generating a settlement cloud map;

[0013] Step S6, calculating the settlement of each node at the bottom of each stratum based on the stratum distribution in the three-dimensional geological model, the compression modulus of each soil layer lithology and the depth of the soil layer;

[0014] Step S7, performing shelling processing on the three-dimensional geological model with the settlement amount of each point calculated, and creating a vertex buffer to store the vertex data of the triangulated mesh, wherein the vertex data of the triangulated mesh includes the settlement attribute values ​​of each vertex obtained by the settlement calculation in steps S5 and S6 and the current surface settlement monitoring time;

[0015] Step S8, using WebGL technology to dynamically display the calculated structure, thereby achieving adaptive correction and dynamic display of three-dimensional geology.

[0016] Preferably, in step S1, the geological survey data includes compression models of various rock and soil layers, geological drilling data, and monitoring point coordinates and elevation information.

[0017] Preferably, in step S2, two geometric parameters are defined:

[0018] The length scale is used to assign a geometric parameter to a grid point, which is defined as the average distance from the point to the two adjacent grid points on the boundary. The length scale of the internal grid points is obtained by interpolating the length scale of the grid points on the boundary using the following inverse principle. Let this point be Q, then the length scale of this point is L(Q);

[0019] Dimensionless radius of the circumscribed circle of a triangle. Let any Delaunay triangle be Δk, the radius of its circumscribed circle be r(k), and the length scale of the center of its circumscribed circle be L(k). Then the dimensionless radius of its circumscribed circle R(k) is:

[0020]

[0021] Preferably, in step S2, the step of automatically adding points and triangulating the settlement range according to the length scale and the dimensionless radius of the triangle circumscribed circle includes:

[0022] Step S21, applying the reciprocal principle, calculating the length scale L(k) of the center of the circumscribed circle of all generated Delaunay triangles and the radius r(k) of the circumscribed circle;

[0023] Step S22, calculating the dimensionless radius R(k) of the circumscribed circle of all triangles;

[0024] Step S23, sorting all existing triangles according to the size of the dimensionless radius R(k) of the circumscribed circle, and the triangle with the largest dimensionless radius R(k) is at the front of the sequence;

[0025] Step S24, adding a new inner point P at the center of the circumscribed circle of the triangle that does not satisfy the circumscribed circle dimensionless radius R(k)<1 and the triangle whose circumscribed circle dimensionless radius R(k) is the largest;

[0026] Step S25, eliminating the common edges of the triangles whose circumscribed circles contain point P to form a Delaunay cavity, connecting point P with related grid points, and generating a new set of Delaunay triangles;

[0027] Step S26, looping the steps S24 to S25 until the dimensionless radius R(k) of the circumscribed circles of all triangles is less than 1;

[0028] Step S27, if the side length of the triangle exceeds the set maximum side length of the triangle, the side is divided equally, and the equally divided points are connected with the corresponding triangle vertices to form a new triangle, wherein the number of equally divided triangles is related to the multiple of the side length exceeding the set maximum side length of the triangle;

[0029] Step S28, for the newly generated triangle, if the dimensionless radius of the circumscribed circle of the triangle is R(k)<1, loop the steps S24-S25 until the dimensionless radius of the circumscribed circle of all triangles is R(k)<1, and finally obtain a planar triangulated network.

[0030] Preferably, step S3 comprises:

[0031] Step S31, using the location information of the geological survey borehole and the geological stratification information of the borehole to establish a reasonable standard stratum;

[0032] Step S32, selecting a suitable interpolation algorithm to construct corresponding stratum levels according to the scope of the subsidence area and the obtained plane triangulation;

[0033] Step S33, closing each stratigraphic layer to form a final three-dimensional geological model.

[0034] Preferably, step S5 comprises:

[0035] The first law of geoscience is used to calculate the settlement value of each node of the triangulated surface, and the specific settlement value of the adjacent unit of a node is calculated based on the extracted layer node information;

[0036] According to the first law of geoscience, the estimated settlement value of a node p(x0,y0,z0) on the surface is S', which is the weighted sum of the settlement of some known points. The formula is:

[0037]

[0038] Where n is the number of nodes in the neighborhood of point p(x0,y0,z0), S i represents the settlement value of the ith node in the neighborhood, w i Represents the weight coefficient of the i-th node;

[0039] To obtain accurate results, w i The difference E between the estimated settlement value S' and the true value S0 at point p (x0, y0, z0) should be minimized, then the following formula is satisfied:

[0040] E(S-S0)=0;

[0041] At the same time, the minimum value of the variance of the difference between the estimated value S' of the node settlement and the true value S0 is calculated, that is:

[0042] min(σ 2 ) = min(Var(S'-S0));

[0043] The settlement value S = S (x, y, z) at any point in space is composed of the regional average value c and the random deviation R (x, y, z), where the variance of the deviation is a constant, then:

[0044] S(x,y,z)=E[S(x,y,z)+R(x,y,z)]=c+R(x,y,z);

[0045] By unbiased constraints:

[0046]

[0047] And optimizing the objective function, we can get the following matrix equation:

[0048]

[0049] For S i =S(x i ,y i ,z i ), S j =S(x j ,y j ,z j (Its semivariance is:

[0050]

[0051] In the formula, γ ij It is expressed as the semivariance between the i-th and j-th neighboring points, i.e., λ is a constant, where γ ij (j≠0) is known, γ ij (j=0) is unknown, and with the help of a suitable fitting function γ=f(d), γ ij (j=0) is solved and substituted into the following formula:

[0052]

[0053] Find w1,…,w n , and put it into the following formula to get the final estimate of S':

[0054]

[0055] Finally, the settlement value of each node on the surface is interpolated through the algorithm and a settlement cloud map is generated.

[0056] Preferably, in step S6, the settlement value of each triangulated corner point in the model is calculated according to the relationship between the settlement amount of different soil layers and the depth and compression modulus, as well as the total settlement, so as to realize the adaptive deformation calculation of the model, that is, assuming that the depth range of the geological model is the settlement deformation range, the calculation formula for the total settlement s is obtained as follows:

[0057]

[0058] Where a is the settlement coefficient, a constant related to the settlement empirical coefficient and the additional surface load, and A i is the integral value of the additional stress coefficient of the i-th layer of soil along the soil layer thickness, which can be obtained according to the formula of soil layer thickness and additional stress coefficient, E si is the compression modulus of the i-th soil layer;

[0059] According to the above calculation formula of total settlement s, the compression deformation s of the i-th layer can be derived i Calculation formula:

[0060]

[0061] Based on the settlement monitoring data, the settlement value at each location on the surface has been calculated as the total settlement s at that point in the three-dimensional geological model;

[0062] In the three-dimensional stratum model, the settlement of all points in the vertical direction of the point is calculated according to the above formula. The settlement value of each layer is s i Further obtain the settlement of all points.

[0063] Preferably, the step S8 comprises:

[0064] Step S81, write vertex shader and fragment shader, define color rendering logic, convert vertex settlement attribute values ​​into RGB arrays, and use WebGL to render and display the model;

[0065] Step S82, through the surface settlement monitoring time, retrieve the settlement value of each vertex of the three-dimensional geological model, adjust the Z value coordinate of each vertex, realize the adaptive correction deformation of the three-dimensional geological model based on the settlement monitoring data, and dynamically display the settlement changes of the three-dimensional geological model.

[0066] The beneficial effects of the present invention are as follows: the present invention uses the first law of geology to perform interpolation calculation on on-site surface settlement monitoring data to obtain a surface settlement cloud map; uses terrain points and monitoring points as sample data to perform Denaulay triangulation on a three-dimensional grid surface within the calculation settlement range, and further constructs a three-dimensional geological model by generating a ground surface based on interpolation of the triangulated grid surface; combines the three-dimensional geological model with geological survey data, calculates settlement values ​​of each node of the three-dimensional geological model, and realizes adaptive correction and dynamic display of three-dimensional geology; the present invention dynamically optimizes the three-dimensional model through real-time data feedback of surface settlement monitoring data, so as to achieve more accurate geological analysis and prediction. BRIEF DESCRIPTION OF THE DRAWINGS

[0067] Figure 1 The process of the method for adaptively correcting a three-dimensional geological model based on settlement monitoring data of the present invention is as follows;

[0068] Figure 2 It is a specific flow chart of obtaining the plane information of topographic points and geological boreholes and the boundary range for later calculation of settlement in the present invention, and performing Denaulay triangulation within the range;

[0069] Figure 3 It is a specific flow chart of constructing a three-dimensional geological model in the present invention;

[0070] Figure 4 It is a specific flow chart of dynamically displaying the computing structure using WebGL technology in the present invention;

[0071] Figure 5 is a schematic diagram of the dimensionless radius of the circumscribed circle of a triangle in the present invention;

[0072] Figure 6 is a schematic diagram of the grid distribution and stratum distribution of the three-dimensional geological model in an embodiment of the present invention;

[0073] Figure 7 is a planar position diagram of each monitoring point in space in an embodiment of the present invention;

[0074] Figure 8 is a surface settlement cloud map in an embodiment of the present invention;

[0075] Fig. 9 It is a dynamic display process diagram in an embodiment of the present invention. DETAILED DESCRIPTION

[0076] The technical solution of the present invention will be described clearly and completely below in conjunction with the accompanying drawings. In the description of the present invention, it should be noted that the terms "center", "up", "down", "left", "right", "vertical", "horizontal", "inside", "outside" and the like indicate directions or positional relationships based on the directions or positional relationships shown in the accompanying drawings, which are only for the convenience of describing the present invention and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific direction, be constructed and operated in a specific direction, and therefore cannot be understood as a limitation on the present invention. In addition, the terms "first", "second", and "third" are used for descriptive purposes only and cannot be understood as indicating or implying relative importance.

[0077] The present invention will be further described below in conjunction with the accompanying drawings.

[0078] like Figure 1 As shown, a three-dimensional geological model adaptive correction method based on settlement monitoring data includes:

[0079] Step S1, obtaining and arranging site geological survey data.

[0080] In the embodiment, the geological survey data includes compression models of various rock and soil layers, geological drilling data, and monitoring point coordinates and elevation information.

[0081] Step S2, obtaining the plane information of the topographic points and geological boreholes and the boundary range used for the later calculation of settlement, performing Denaulay triangulation within the range, and interpolating it into a surface.

[0082] In the embodiment, two geometric parameters are defined:

[0083] The length scale is used to assign a geometric parameter to a grid point, which is defined as the average distance from the point to the two adjacent grid points on the boundary. The length scale of the internal grid points is obtained by interpolating the length scale of the grid points on the boundary using the following inverse principle. Let this point be Q, then the length scale of this point is L(Q);

[0084] Assume that point Q is to be inserted into the next Figure 5 A point in the middle triangle 145, then the length scale L(Q) of Q is:

[0085]

[0086] Where L(1), L(4), L(5) are the length scales of the boundary grid points; l1, l4, l5 are the distances from point Q to the three vertices;

[0087] Dimensionless radius of the circumscribed circle of a triangle. Let any Delaunay triangle be Δk, the radius of its circumscribed circle be r(k), and the length scale of the center of its circumscribed circle be L(k). Then the dimensionless radius of its circumscribed circle R(k) is:

[0088]

[0089] like Figure 2 As shown, the steps of automatically adding points and triangulating the settlement range according to the length scale and the dimensionless radius of the triangle circumscribed circle include:

[0090] Step S21, applying the reciprocal principle, calculates the length scale L(k) of the center of the circumscribed circle of all generated Delaunay triangles and the radius r(k) of the circumscribed circle.

[0091] Step S22, calculating the dimensionless radius R(k) of the circumscribed circles of all triangles.

[0092] Step S23, sort all existing triangles according to the size of the dimensionless radius R(k0) of the circumscribed circle, and the triangle with the largest dimensionless radius R(k0) of the circumscribed circle is at the front of the sequence.

[0093] Step S24, adding a new inner point P at the center of the circumscribed circle of the triangle that does not satisfy the circumscribed circle dimensionless radius R(k)<1 and whose circumscribed circle dimensionless radius R(k) is the largest.

[0094] Step S25, eliminating the common edges of the triangles whose circumscribed circles contain point P to form a Delaunay cavity, connecting point P with related grid points, and generating a new set of Delaunay triangles.

[0095] Step S26, looping the steps S24 to S25 until the dimensionless radius R(k) of the circumscribed circles of all triangles is less than 1.

[0096] Step S27, for a triangle whose side length exceeds the set maximum side length of a triangle, the side is divided into equal parts, and the dividing points are connected to the corresponding triangle vertices to form a new triangle, wherein the number of equal parts is related to the multiple of the side length exceeding the set maximum side length of the triangle.

[0097] Step S28, for the newly generated triangle, if the dimensionless radius of the circumscribed circle of the triangle is R(k)<1, loop the steps S24-S25 until the dimensionless radius of the circumscribed circle of all triangles is R(k)<1, and finally obtain a planar triangulated network.

[0098] Step S3, constructing a three-dimensional geological model.

[0099] like Figure 3 As shown, step S3 includes:

[0100] Step S31, using the location information of geological survey boreholes and the geological stratification information of the boreholes to establish a reasonable standard stratum.

[0101] Step S32, selecting a suitable interpolation algorithm to construct corresponding stratigraphic levels according to the scope of the subsidence area and the obtained plane triangulation.

[0102] Step S33, closing each stratigraphic layer to form a final three-dimensional geological model.

[0103] Step S4, obtaining and collating the current settlement monitoring data of all surface settlement monitoring points within the site settlement area.

[0104] Step S5, calculating the surface settlement and generating a settlement cloud map.

[0105] In the embodiment, step S5 includes:

[0106] The first law of geoscience is used to calculate the settlement value of each node of the triangulated surface, and the specific settlement value of the adjacent unit of a node is calculated based on the extracted layer node information;

[0107] According to the first law of geoscience, the estimated settlement value of a node p(x0,y0,z0) on the surface is S', which is the weighted sum of the settlement of some known points. The formula is:

[0108]

[0109] Where n is the number of nodes in the neighborhood of point p(x0,y0,z0), S i represents the settlement value of the ith node in the neighborhood, w i Represents the weight coefficient of the i-th node;

[0110] To obtain accurate results, w i The difference E between the estimated settlement value S' and the true value S0 at point p (x0, y0, z0) should be minimized, then the following formula is satisfied:

[0111] E(S-S0)=0;

[0112] At the same time, the minimum value of the variance of the difference between the estimated value S' of the node settlement and the true value S0 is calculated, that is:

[0113] min(σ 2 ) = min(Var(S'-S0));

[0114] The settlement value S = S (x, y, z) at any point in space is composed of the regional average value c and the random deviation R (x, y, z), where the variance of the deviation is a constant, then:

[0115] S(x,y,z)=E[S(x,y,z)+R(x,y,z)]=c+R(x,y,z);

[0116] By unbiased constraints:

[0117]

[0118] And optimizing the objective function, we can get the following matrix equation:

[0119]

[0120] For S i =S(x i ,y i ,z i ), S j =S(x j ,y j ,z j (Its semivariance is:

[0121]

[0122] In the formula, γ ij It is expressed as the semivariance between the i-th and j-th neighboring points, i.e., λ is a constant, where γ ij (j≠0) is known, γ ij (j=0) is unknown, and with the help of a suitable fitting function γ=f(d), γ ij (j=0) is solved and substituted into the following formula:

[0123]

[0124] Find w1,…,w n , and put it into the following formula to get the final estimate of S':

[0125]

[0126] Finally, the settlement value of each node on the surface is interpolated through the algorithm and a settlement cloud map is generated.

[0127] Step S6, based on the stratum distribution in the three-dimensional geological model, the compression modulus of each soil layer lithology and the depth of the soil layer, the settlement of each node at the bottom of each stratum is calculated.

[0128] In the embodiment, according to the relationship between the settlement amount of different soil layers and the depth and compression modulus, as well as the total settlement, the settlement value of each triangulated corner point in the model is calculated, so as to realize the adaptive deformation calculation of the model, that is, assuming that the depth range of the geological model is the settlement deformation range, the calculation formula of the total settlement s is obtained as follows:

[0129]

[0130] Where a is the settlement coefficient, a constant related to the settlement empirical coefficient and the additional surface load, and A i is the integral value of the additional stress coefficient of the i-th layer of soil along the soil layer thickness, which can be obtained according to the formula of soil layer thickness and additional stress coefficient, E si is the compression modulus of the i-th soil layer;

[0131] According to the above calculation formula of total settlement s, the compression deformation s of the i-th layer can be derived i Calculation formula:

[0132]

[0133] Based on the settlement monitoring data, the settlement value at each location on the surface has been calculated as the total settlement s at that point in the three-dimensional geological model;

[0134] In the three-dimensional stratum model, the settlement of all points in the vertical direction of the point is calculated according to the above formula. The settlement value of each layer is s i Further obtain the settlement of all points.

[0135] Step S7, shelling the three-dimensional geological model with the settlement amount of each point calculated, and creating a vertex buffer to store the vertex data of the triangulated mesh, wherein the vertex data of the triangulated mesh includes the settlement attribute values ​​of each vertex obtained by the settlement calculation in steps S5 and S6 and the current surface settlement monitoring time.

[0136] Step S8, using WebGL technology to dynamically display the calculated structure, thereby achieving adaptive correction and dynamic display of three-dimensional geology.

[0137] like Figure 4 As shown, step S8 includes:

[0138] Step S81, write vertex shader and fragment shader, define color rendering logic, convert vertex settlement attribute values ​​into RGB arrays, and use WebGL to render and display the model.

[0139] Step S82, through the surface settlement monitoring time, retrieve the settlement value of each vertex of the three-dimensional geological model, adjust the Z value coordinate of each vertex, realize the adaptive correction deformation of the three-dimensional geological model based on the settlement monitoring data, and dynamically display the settlement changes of the three-dimensional geological model.

[0140] The specific implementation method is described in detail according to the following embodiments:

[0141] A site with an area of ​​5652.18 m2 was selected in the region. The geotechnical survey data (compression modulus of each rock and soil layer), coordinate information of each monitoring point, and surface settlement information are shown in Tables 1 and 2.

[0142] Geotechnical Name Compression modulus / MPa D-Medium sand 35 D-Silty sand 15 R-Sandy silt-Clayey silt 20 R-Clay 10

[0143] Table 1

[0144] Subsidence monitoring points X coordinate Y coordinate Z coordinate Sedimentation value / m CJD01 381240.00 3539305.00 36.67 1.0 CJD02 381195.00 3539355.00 35.69 1.2 CJD03 381225.00 3539330.00 36.27 1.0 CJD04 381190.00 3539300.00 39.05 1.1 CJD05 381245.00 3539360.00 34.87 1.3

[0145] Table 2

[0146] According to the size of the site, the location information of the exploration boreholes and the geological stratification information, a three-dimensional geological model with tetrahedrons as units is established using triangulation technology. The grid distribution and stratum distribution of the model are shown in the figure below. Figure 6 As shown in the figure, the plane position diagram of each monitoring point in space is as follows Figure 7 shown.

[0147] According to the monitoring point coordinate information and the settlement data of the monitoring points provided in Table 2, combined with the terrain model extracted from the three-dimensional geological model, the basic data and grid information of the settlement cloud map can be formed. Relying on the data of the five known monitoring points, according to the interpolation calculation process of step S5, the settlement value of each node on the surface can be known, and the surface settlement can be displayed in the form of a cloud map. The settlement cloud map is as follows: Figure 8 shown.

[0148] According to the surface settlement cloud map calculated in step S5, the settlement value of each position on the surface is the total settlement value of the point; for example, the coordinates of the settlement monitoring point CJD01 are (381240.00, 3539305.00, 36.67), and its settlement value is 1.0m. Then the settlement value of each stratum layer at the corresponding position in the vertical direction of the position (381240.00, 3539305.00) can be obtained by Similarly, the settlement values ​​at all locations in the three-dimensional geological model can be calculated and dynamically displayed. The dynamic display process is as follows: Fig. 9 shown.

[0149] The above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit the same. Although the present invention has been described in detail with reference to the above embodiments, those skilled in the art should understand that the technical solutions described in the above embodiments may still be modified, or some or all of the technical features may be replaced by equivalents. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the scope of the technical solutions of the embodiments of the present invention.

Claims

1. A three-dimensional geological model adaptive correction method based on settlement monitoring data, characterized in that: include: Step S1, obtaining and collating site geological survey data; Step S2, obtaining the plane information of the topographic points and geological boreholes and the boundary range used for the later calculation of settlement, performing Denaulay triangulation within the range, and interpolating it into a surface; Step S3, constructing a three-dimensional geological model; Step S4, obtaining and collating the current settlement monitoring data of all surface settlement monitoring points within the settlement area of ​​the site; Step S5, calculating the surface settlement and generating a settlement cloud map; Step S6, calculating the settlement of each node at the bottom of each stratum based on the stratum distribution in the three-dimensional geological model, the compression modulus of each soil layer lithology and the depth of the soil layer; Step S7, performing shelling processing on the three-dimensional geological model with the settlement amount of each point calculated, and creating a vertex buffer to store the vertex data of the triangulated mesh, wherein the vertex data of the triangulated mesh includes the settlement attribute values ​​of each vertex obtained by the settlement calculation in steps S5 and S6 and the current surface settlement monitoring time; Step S8, using WebGL technology to dynamically display the calculated structure, thereby achieving adaptive correction and dynamic display of three-dimensional geology.

2. A three-dimensional geological model adaptive correction method based on settlement monitoring data according to claim 1, characterized in that: In step S1, the geological survey data includes compression models of various rock and soil layers, geological drilling data, and monitoring point coordinates and elevation information.

3. The method for adaptively correcting a three-dimensional geological model based on settlement monitoring data according to claim 1, characterized in that: In step S2, two geometric parameters are defined: The length scale is used to assign a geometric parameter to a grid point, which is defined as the average of the distance from the point to the two adjacent grid points on the boundary. The length scale of the internal grid points is obtained by interpolating the length scale of the grid points on the boundary using the following inverse principle. Let this point be Q, then the length scale of this point is L(Q); Dimensionless radius of the circumscribed circle of a triangle. Let any Delaunay triangle be Δk, the radius of its circumscribed circle be r(k), and the length scale of the center of its circumscribed circle be L(k). Then the dimensionless radius of its circumscribed circle R(k) is:

4. A three-dimensional geological model adaptive correction method based on settlement monitoring data according to claim 3, characterized in that: In step S2, the step of automatically adding points and triangulating the settlement range according to the length scale and the dimensionless radius of the triangle circumscribed circle includes: Step S21, applying the reciprocal principle, calculating the length scale L(k) of the center of the circumscribed circle of all generated Delaunay triangles and the radius r(k) of the circumscribed circle; Step S22, calculating the dimensionless radius R(k) of the circumscribed circle of all triangles; Step S23, sorting all existing triangles according to the size of the dimensionless radius R(k) of the circumscribed circle, and the triangle with the largest dimensionless radius R(k) is at the front of the sequence; Step S24, adding a new inner point P at the center of the circumscribed circle of the triangle that does not satisfy the circumscribed circle dimensionless radius R(k)<1 and the triangle whose circumscribed circle dimensionless radius R(k) is the largest; Step S25, eliminating the common edges of the triangles whose circumscribed circles contain point P to form a Delaunay cavity, connecting point P with related grid points, and generating a new set of Delaunay triangles; Step S26, looping the steps S24 to S25 until the dimensionless radius R(k) of the circumscribed circles of all triangles is less than 1; Step S27, if the side length of the triangle exceeds the set maximum side length of the triangle, the side is divided equally, and the equally divided points are connected with the corresponding triangle vertices to form a new triangle, wherein the number of equally divided triangles is related to the multiple of the side length exceeding the set maximum side length of the triangle; Step S28, for the newly generated triangle, if the dimensionless radius of the circumscribed circle of the triangle is R(k)<1, loop the steps S24-S25 until the dimensionless radius of the circumscribed circle of all triangles is R(k)<1, and finally obtain a planar triangulated network.

5. The method for adaptively correcting a three-dimensional geological model based on settlement monitoring data according to claim 1, characterized in that: The step S3 comprises: Step S31, using the location information of the geological survey borehole and the geological stratification information of the borehole to establish a reasonable standard stratum; Step S32, selecting a suitable interpolation algorithm to construct corresponding stratum levels according to the scope of the subsidence area and the obtained plane triangulation; Step S33, closing each stratigraphic layer to form a final three-dimensional geological model.

6. The method for adaptively correcting a three-dimensional geological model based on settlement monitoring data according to claim 1, characterized in that: The step S5 comprises: The first law of geoscience is used to calculate the settlement value of each node of the triangulated surface, and the specific settlement value of the adjacent unit of a node is calculated based on the extracted layer node information; According to the first law of geoscience, the estimated settlement value of a node p(x0,y0,z0) on the surface is S', which is the weighted sum of the settlement of some known points. The formula is: Where n is the number of nodes in the neighborhood of point p(x0,y0,z0), S i represents the settlement value of the ith node in the neighborhood, w i Represents the weight coefficient of the i-th node; To obtain accurate results, w i The difference E between the estimated settlement value S' and the true value S0 at point p (x0, y0, z0) should be minimized, then the following formula is satisfied: E(S-S0)=0; At the same time, the minimum value of the variance of the difference between the estimated value S' of the node settlement and the true value S0 is calculated, that is: min(σ 2 )=min(Var(S'-S0)); The settlement value S = S (x, y, z) at any point in space is composed of the regional average value c and the random deviation R (x, y, z), where the variance of the deviation is a constant, then: S(x,y,z)=E[S(x,y,z)+R(x,y,z)]=c+R(x,y,z); By unbiased constraints: And optimizing the objective function, we can get the following matrix equation: For S i =S(x i ,y i ,z i ), S j =S(x j ,y j ,z j ) Its semivariance is: In the formula, γ ij It is expressed as the semivariance between the i-th and j-th neighboring points, i.e., λ is a constant, where γ ij (j≠0) is known, γ ij (j=0) is unknown, and with the help of a suitable fitting function γ=f(d), γ ij (j=0) is solved and substituted into the following formula: Find w1,…,w n , and put it into the following formula to get the final estimate of S': Finally, the settlement value of each node on the surface is interpolated through the algorithm and a settlement cloud map is generated.

7. The method for adaptively correcting a three-dimensional geological model based on settlement monitoring data according to claim 1, characterized in that: In step S6, the settlement value of each triangulated corner point in the model is calculated according to the relationship between the settlement amount of different soil layers and the depth and compression modulus, as well as the total settlement, so as to realize the adaptive deformation calculation of the model, that is, assuming that the depth range of the geological model is the settlement deformation range, the calculation formula of the total settlement s is obtained as follows: Where a is the settlement coefficient, a constant related to the settlement empirical coefficient and the additional surface load, and A i is the integral value of the additional stress coefficient of the i-th layer of soil along the soil layer thickness, which can be obtained according to the formula of soil layer thickness and additional stress coefficient, E si is the compression modulus of the i-th soil layer; According to the above calculation formula of total settlement s, the compression deformation s of the i-th layer can be derived i Calculation formula: Based on the settlement monitoring data, the settlement value at each location on the surface has been calculated as the total settlement s at that point in the three-dimensional geological model; In the three-dimensional stratum model, the settlement of all points in the vertical direction of the point is calculated according to the above formula. The settlement value of each layer is s i Further obtain the settlement of all points.

8. The method for adaptively correcting a three-dimensional geological model based on settlement monitoring data according to claim 1, characterized in that: The step S8 comprises: Step S81, write vertex shader and fragment shader, define color rendering logic, convert vertex settlement attribute values ​​into RGB arrays, and use WebGL to render and display the model; Step S82, through the surface settlement monitoring time, retrieve the settlement value of each vertex of the three-dimensional geological model, adjust the Z value coordinate of each vertex, realize the adaptive correction deformation of the three-dimensional geological model based on the settlement monitoring data, and dynamically display the settlement changes of the three-dimensional geological model.

Images