Caprock three-dimensional geological modeling method, system and storage medium

Abstract

This invention discloses a method, system, and storage medium for three-dimensional geological modeling of overburden layers, relating to the field of three-dimensional geological modeling of overburden layers. The method includes determining a standard profile based on the distribution of hydropower engineering exploration points and geotechnical engineering investigation information; the geotechnical engineering investigation information includes: borehole data, elevation information, and modeling range; extracting the minimum horizontal distance from the spatial points on the top surface of the overburden layer of the standard profile to the boundary line between the overburden layer and the surface bedrock, the relative elevation of the points to the surface bedrock, and the overburden layer depth based on the standard profile; constructing a polynomial model with the minimum horizontal distance from the spatial points on the top surface of the standard profile to the boundary line between the overburden layer and the surface bedrock and the relative elevation of the points to the surface bedrock as input, and the overburden layer depth as output; interpolating the nodes at the bottom surface of the deep overburden layer using the polynomial model to obtain the spatial information of the nodes at the bottom surface of the overburden layer, and establishing a three-dimensional geological model of the overburden layer. This invention can better reflect the true spatial distribution characteristics of the overburden layer.

Description

Technical Field

[0001] This invention relates to the field of three-dimensional geological modeling of overburden, and in particular to a method, system and storage medium for three-dimensional geological modeling of overburden. Background Technology

[0002] A three-dimensional geological model of the overburden layer is a model reflecting the spatial distribution characteristics of loose sediments on the ground. Overburden layers have a significant impact on engineering construction, ecological protection, and construction safety. Hydropower project overburden geological models are characterized by large areas, complex formations, and uneven spatial distribution. When establishing large-scale hydropower overburden three-dimensional geological models, two main methods are used: one is an automatic construction method based on existing borehole and trench data, and the other is a model reconstruction method based on standard geological profiles. However, limitations such as large modeling areas, limited existing borehole data, and significant topographic relief mean that using only discrete single-point borehole data leads to slow modeling speeds and models that violate geological laws. Therefore, establishing a three-dimensional geological model of hydropower overburden often requires data interpolation. During interpolation, the spatial location of the interpolation points and their relative position to the lithological boundary (the boundary between the surface bedrock and the overburden layer) must be considered.

[0003] Conventional model interpolation methods involve drawing constraint profiles based on existing borehole depth data of the overburden, connecting each constraint profile with the overburden boundary, and then linearly interpolating to generate the spatial positions of the remaining surface nodes.

[0004] However, conventional interpolation methods encounter the following problems: First, due to the large spatial scope of hydropower projects, the wide distribution of overburden, and the limited number of control boreholes, there are large unexplored areas where it is not feasible to draw constrained profiles. Second, when generating constrained profiles from discrete borehole data, only the spatial distribution of overburden in the profile strike direction is often considered, which differs significantly from the actual spatial distribution of overburden. The resulting 3D geological model of overburden often shows conflicting distribution trends in different directions, leading to abrupt changes that do not conform to geological laws. Third, the relative position of overburden thickness and overburden boundary, as well as the terrain features of the overburden spatial points, are not considered. In areas with steep slopes, the bottom elevation of the created 3D geological model of overburden may be higher than the surface. Fourth, simple linear interpolation methods have weak fitting capabilities, and the linear interpolation function cannot fully reflect the spatial distribution characteristics of overburden.

[0005] Therefore, based on the above problems, there is an urgent need to provide a three-dimensional geological modeling method or system for overburden that can better reflect the true spatial distribution characteristics of overburden. Summary of the Invention

[0006] The purpose of this invention is to provide a three-dimensional geological modeling method, system, and storage medium for overburden layers, which can better reflect the true spatial distribution characteristics of overburden layers.

[0007] To achieve the above objectives, the present invention provides the following solution:

[0008] A three-dimensional geological modeling method for overburden layers, comprising:

[0009] Standard profiles are determined based on the distribution of hydropower engineering exploration points and geotechnical engineering investigation information; the geotechnical engineering investigation information includes: borehole data, elevation information, and modeling scope;

[0010] The minimum horizontal distance from the spatial point on the top surface of the standard profile overburden to the boundary line of the bedrock overburden, the relative elevation of the point to the bedrock, and the depth of the overburden are extracted based on the standard profile. The relative elevation of the point to the bedrock is the relative elevation of the boundary point between the spatial point on the top surface of the profile overburden and the spatial point on the top surface of the standard profile overburden with the minimum horizontal distance from the boundary line of the bedrock overburden.

[0011] A polynomial model is constructed with the minimum horizontal distance from the spatial point on the top surface of the standard profile overburden to the boundary line of the bedrock overburden and the relative elevation of the point to the bedrock as input, and the overburden depth as output.

[0012] By interpolating the nodes at the bottom of the deep overburden layer using a polynomial model, spatial information of the nodes at the bottom of the overburden layer is obtained, and a three-dimensional geological model of the overburden layer is established.

[0013] Optionally, determining the standard profile based on the distribution of hydropower engineering exploration points and the amount of geotechnical engineering investigation information specifically includes:

[0014] AutoCAD was used to number the exploration points for hydropower projects and to lay out multiple standard profiles;

[0015] Number the standard sections of the layout;

[0016] Preprocessing of geotechnical engineering investigation information;

[0017] The three-dimensional topographic surface is determined based on the pre-processed geotechnical engineering investigation information;

[0018] The standard profile overburden layer is determined based on the three-dimensional topographic surface and the standard profile.

[0019] Optionally, the preprocessing of geotechnical engineering investigation information specifically includes:

[0020] Use the CGCS2000 coordinate system to unify the processing of geotechnical engineering investigation information;

[0021] Use AutoCAD to extract elevation points, primary curves, and calculation curves of the terrain within the modeling area, and delete outliers;

[0022] Convert the format of the borehole data.

[0023] Optionally, determining the three-dimensional topographic surface based on the preprocessed geotechnical engineering investigation information specifically includes:

[0024] Use ItasCAD to generate surfaces from the modeling area;

[0025] Use ItasCAD to generate the first elevation point file from the calculation curve and the primary curve data;

[0026] Use ItasCAD to generate a second elevation point file from the elevation points;

[0027] Determine the elevation point file based on the first elevation point file and the second elevation point file;

[0028] The 3D terrain surface is determined by applying fuzzy constraints to the surface generated within the modeling range based on the elevation point file.

[0029] Optionally, the step of extracting the minimum horizontal distance from the spatial point on the top surface of the standard profile overburden to the boundary line of the bedrock overburden, the relative elevation of the point to the bedrock, and the depth of the overburden based on the standard profile specifically includes:

[0030] Using formula Determine the minimum horizontal distance from the spatial point on the top surface of the standard profile overburden layer to the boundary line of the surface bedrock overburden layer;

[0031] Using the formula H = P z -P mz Determine the relative elevation of the point from the bedrock surface;

[0032] The depth of the overburden layer is determined using the formula D = Pz - Pbz;

[0033] Where Dmin represents the minimum horizontal distance from the top node of the overburden layer to the boundary line between the bedrock and the overburden layer; where Pmx i Let Pmy be the x-coordinate of the i-th node on the boundary line between the bedrock overburden and the surface. i Let Py be the y-coordinate of the i-th node on the boundary line between the bedrock and the surface overburden, n be the total number of boundary points constituting the boundary, Py be the ordinate of the spatial position of the node on the top surface of the overburden, Px be the abscissa of the spatial position of the node on the top surface of the overburden, H be the relative elevation of the point from the bedrock, and P be the y-coordinate of the node. z P represents the elevation of the node on the top surface of the overburden layer. mz denoted as , where is the elevation of the boundary node, D is the depth of the overburden layer, and Pbz is the elevation of the node at the bottom of the overburden layer.

[0034] Optionally, the polynomial model is:

[0035]

[0036] Where x is the input sample, and a = {a0, a1, ..., a...} m} represents the input sample coefficients, where m is determined by the maximum power of the polynomial. This is the predicted value from the regression model, and the predicted value is the cover depth.

[0037] A three-dimensional geological modeling system for overburden layers, comprising:

[0038] The standard profile determination module is used to determine the standard profile based on the distribution of hydropower engineering exploration points and geotechnical engineering investigation information; the geotechnical engineering investigation information includes: borehole data, elevation information, and modeling range.

[0039] The standard profile extraction module is used to extract the minimum horizontal distance between the spatial points on the top surface of the standard profile overburden and the boundary line of the bedrock overburden, the relative elevation of the points from the bedrock, and the depth of the overburden based on the standard profile. The relative elevation of the points from the bedrock is the relative elevation of the boundary point between the spatial points on the top surface of the profile overburden and the spatial points on the top surface of the standard profile overburden and the boundary line of the bedrock overburden.

[0040] The polynomial model building module is used to build a polynomial model with the minimum horizontal distance from the spatial point on the top surface of the standard profile overburden to the boundary line of the bedrock overburden and the relative elevation of the point to the bedrock as input, and the overburden depth as output.

[0041] The overburden 3D geological model construction module is used to interpolate the nodes of the bottom surface of the deep overburden using a polynomial model, obtain the spatial information of the nodes of the bottom surface of the overburden, and establish a 3D geological model of the overburden.

[0042] A storage medium storing computer program instructions, which, when executed by a processor, implement the aforementioned three-dimensional geological modeling method for overburden layers.

[0043] According to specific embodiments provided by the present invention, the present invention discloses the following technical effects:

[0044] The present invention provides a method, system, and storage medium for three-dimensional geological modeling of overburden layers. It uses the minimum horizontal distance between spatial points on the top surface of the overburden layer and the boundary line between the overburden layer and the bedrock surface, and the relative elevation of the points to the bedrock surface, to constrain interpolation. Considering the overall spatial distribution of the overburden layer, the constructed three-dimensional geological model of the overburden layer better conforms to geological laws. It uses a nonlinear fitting function instead of a linear fitting function, which can better fit the spatial distribution data of the overburden layer, making the three-dimensional geological model of the overburden layer more closely match the actual spatial distribution of the overburden layer. Furthermore, the use of a multinomial regression-based three-dimensional geological modeling method for overburden layers can obtain the same number of constraint points as the nodes of the overburden layer surface, greatly accelerating the overburden layer modeling speed. Attached Figure Description

[0045] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0046] Figure 1 This is a schematic diagram of a three-dimensional geological modeling method for overburden layers provided by the present invention. Detailed Implementation

[0047] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0048] The purpose of this invention is to provide a three-dimensional geological modeling method, system, and storage medium for overburden layers, which can better reflect the true spatial distribution characteristics of overburden layers.

[0049] To make the above-mentioned objects, features and advantages of the present invention more apparent and understandable, the present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0050] like Figure 1 As shown, the present invention provides a three-dimensional geological modeling method for overburden layers, comprising:

[0051] S101, determine the standard profile based on the distribution of hydropower engineering exploration points and geotechnical engineering investigation information; the geotechnical engineering investigation information includes: borehole data, elevation information and modeling range.

[0052] S101 specifically includes:

[0053] S1. AutoCAD is used to number the exploration points for hydropower projects and to lay out multiple standard profiles, with each standard profile passing through as many exploration points as possible.

[0054] S2, number the standard sections to be arranged.

[0055] S3 is used for preprocessing geotechnical engineering investigation information.

[0056] S3 specifically includes:

[0057] The CGCS2000 coordinate system is used to uniformly process geotechnical engineering investigation information, with elevation units in meters (m).

[0058] Use AutoCAD to extract elevation points, primary curves, and calculation curves of the terrain within the modeling area, and remove outliers. Save the resulting files as CSV, XLSX, and DAT formats.

[0059] The drilling data is converted into a file type that can be recognized and loaded by 3D modeling software and Python, typically in Excel spreadsheet format.

[0060] S4, determine the three-dimensional topographic surface based on the pre-processed geotechnical engineering investigation information.

[0061] S4 specifically includes:

[0062] Create a new project in ItasCAD (3D modeling software), import the modeling scope file (dxf format), and generate surfaces using closed rectangles.

[0063] Use ItasCAD to generate the first elevation point file from the calculation curve and the first curve data.

[0064] Use ItasCAD to generate a second elevation point file from the elevation points.

[0065] The elevation point file is determined based on the first elevation point file and the second elevation point file.

[0066] The 3D terrain surface is determined by applying fuzzy constraints to the surface generated within the modeling range based on the elevation point file.

[0067] S5, determine the standard profile overburden layer based on the three-dimensional topographic surface and standard profile.

[0068] S5 specifically includes:

[0069] Import the standard section lines into ItasCAD.

[0070] Import the drilling data into ItasCAD.

[0071] Using a standard profile to cut the cross section, a standard profile topographic line (dxf) with borehole depth projection is obtained.

[0072] Based on the depth information of the overburden revealed by the borehole, draw a standard profile of the bottom surface of the overburden.

[0073] S102, extract the minimum horizontal distance from the spatial point on the top surface of the standard profile overburden to the boundary line of the bedrock overburden, the relative elevation of the point to the bedrock, and the depth of the overburden based on the standard profile. The relative elevation of the point to the bedrock is the relative elevation of the boundary point between the spatial point on the top surface of the profile overburden and the spatial point on the top surface of the standard profile overburden and the minimum horizontal distance from the boundary line of the bedrock overburden.

[0074] S102 also includes:

[0075] Extract the nodes on the top surface of the overlay layer. Use ItasCAD to extract all nodes on the top surface of the overlay layer and save them as an Excel file;

[0076] Extract the bottom line nodes of the overburden layer. Import the standard section into ItasCAD and save the top and bottom line node data of the standard overburden layer section as an Excel file.

[0077] S102 specifically includes:

[0078] Using formula Determine the minimum horizontal distance from the spatial point at the bottom of the standard profile overburden layer to the boundary line of the surface bedrock overburden layer;

[0079] Using formula Determine the minimum horizontal distance from the spatial point on the top surface of the standard profile overburden layer to the boundary line of the surface bedrock overburden layer;

[0080] Using the formula H = P z -P mz Determine the relative elevation of the point from the bedrock surface;

[0081] The depth of the overburden layer is determined using the formula D = Pz - Pbz;

[0082] Where Dmin represents the minimum horizontal distance from the top node of the overburden layer to the boundary line between the bedrock and the overburden layer; where Pmx i Let Pmy be the x-coordinate of the i-th node on the boundary line between the bedrock overburden and the surface. i Let Py be the y-coordinate of the i-th node on the boundary line between the bedrock and the surface overburden, n be the total number of boundary points constituting the boundary, Py be the ordinate of the spatial position of the node on the top surface of the overburden, Px be the abscissa of the spatial position of the node on the top surface of the overburden, H be the relative elevation of the point from the bedrock, and P be the y-coordinate of the node. z P represents the elevation of the node on the top surface of the overburden layer. mz denoted as , where is the elevation of the boundary node, D is the depth of the overburden layer, and Pbz is the elevation of the node at the bottom of the overburden layer.

[0083] S103 constructs a polynomial model with the minimum horizontal distance from the spatial point on the top surface of the standard profile overburden to the boundary line of the bedrock overburden and the relative elevation of the point to the bedrock as inputs, and the overburden depth as output.

[0084] The polynomial model is:

[0085]

[0086] Where x is the input sample, and a = {a0, a1, ..., a...} m} represents the input sample coefficients, where m is determined by the maximum power of the polynomial. This is the predicted value from the regression model, and the predicted value is the cover depth.

[0087] The training of the multinomial model is as follows:

[0088] (1) Using D min The two features, H (minimum horizontal distance from the top node of the overburden layer to the boundary line of the bedrock overburden layer) and H (relative elevation of the point to the bedrock), are used as inputs to the multinomial regression model, and D (overburden layer depth) is used as the output of the multinomial regression model. By extracting standard profiles, more than 2,000 samples can generally be obtained.

[0089] (2) Divide the samples into training dataset and test dataset according to a ratio of 3:1.

[0090] Using the root mean square error as the model loss function, the model loss function is:

[0091]

[0092] Where J(a0, a1, ..., a n ) represents the sample coefficients a = {a0, a1, ..., a...} n The loss value of the model at time}, where m is the total number of samples and h is the total number of samples. a (x0, x1 , ..., x n y represents the model's predicted value (predicted overburden depth). i This represents the actual value (the actual overlay depth).

[0093] Adjusting the hyperparameters of the multinomial regression model. The degree (highest power) parameter of the multinomial in the multinomial regression model affects the model performance. Set different degrees to test the model, use the root mean square error (RMSE) for evaluation, and select the degree with the smallest RMSE as the optimal parameter for the model.

[0094] The model with the smallest root mean square error is selected to interpolate the bottom nodes of all unknown depth overburden layers to obtain the depth of all overburden layer nodes. The elevation of the bottom nodes of the overburden layer is obtained by subtracting the depth of the overburden layer nodes from the elevation of the overburden layer vertices and saving it as an Excel file.

[0095] S104. Using a polynomial model, interpolation of the nodes at the bottom of the deep overburden layer is performed to obtain the spatial information of the nodes at the bottom of the overburden layer, and a three-dimensional geological model of the overburden layer is established.

[0096] S104 specifically includes:

[0097] (1) Copy the top surface of the overlay in ItasCAD and name it the bottom surface of the overlay;

[0098] (2) Import the elevation file of the bottom node of the overburden layer into the 3D modeling software;

[0099] (3) Apply precise constraints to the bottom boundary of the cover layer, and apply fuzzy constraints to other nodes on the ground surface of the cover layer using the elevation points of the bottom surface of the cover layer.

[0100] (4) Use the spatial interpolation method in the 3D modeling software to obtain the bottom surface of the 3D geological model of the overburden layer;

[0101] (5) Merge the top and bottom surfaces of the three-dimensional geological model of the overburden to obtain the envelope of the three-dimensional geological model of the overburden.

[0102] This invention presents a 3D geological modeling method for overburden layers based on multinomial regression. This method is highly applicable and suitable for constructing 3D geological models of overburden layers for hydropower projects at various scales. It requires less drilling data than conventional 3D geological modeling methods and eliminates the need for drawing numerous constraint profiles, thus saving project costs. Furthermore, this invention uses nonlinear functions to fit the spatial distribution data of the overburden layer, resulting in a more accurate and reliable representation of its true spatial distribution compared to conventional methods, and providing stronger guidance for engineering projects. Finally, this invention offers fast modeling speed, low computer performance requirements, and improved work efficiency.

[0103] Corresponding to the above method, the present invention provides a three-dimensional geological modeling system for overburden layers, comprising:

[0104] The standard profile determination module is used to determine the standard profile based on the distribution of hydropower engineering exploration points and geotechnical engineering investigation information; the geotechnical engineering investigation information includes: borehole data, elevation information, and modeling range.

[0105] The standard profile extraction module is used to extract the minimum horizontal distance from the spatial point on the top surface of the standard profile overburden to the boundary line of the bedrock overburden, the relative elevation of the point to the bedrock, and the depth of the overburden based on the standard profile.

[0106] The polynomial model building module is used to construct a polynomial model with the minimum horizontal distance from the spatial point on the top surface of the standard profile overburden to the boundary line of the bedrock overburden and the relative elevation of the point to the bedrock as input, and the overburden depth as output.

[0107] The overburden 3D geological model construction module is used to interpolate the nodes of the bottom surface of the deep overburden using a polynomial model, obtain the spatial information of the nodes of the bottom surface of the overburden, and establish a 3D geological model of the overburden.

[0108] Based on the above description, the technical solution of the present invention, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods of the various embodiments of the present invention. The aforementioned computer storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory, random access memory, magnetic disks, or optical disks.

[0109] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from other embodiments. Similar or identical parts between embodiments can be referred to interchangeably. For the systems disclosed in the embodiments, since they correspond to the methods disclosed in the embodiments, the descriptions are relatively simple; relevant parts can be referred to the method section.

[0110] This document uses specific examples to illustrate the principles and implementation methods of the present invention. The descriptions of the above embodiments are only for the purpose of helping to understand the method and core ideas of the present invention. Furthermore, those skilled in the art will recognize that, based on the ideas of the present invention, there will be changes in the specific implementation methods and application scope. Therefore, the content of this specification should not be construed as a limitation of the present invention.

Claims

1. A method for three-dimensional geological modeling of overburden layers, characterized in that, include: Standard profiles are determined based on the distribution of hydropower engineering exploration points and geotechnical engineering investigation information. The geotechnical engineering investigation information includes: borehole data, elevation information, and modeling scope; The minimum horizontal distance from the spatial point on the top surface of the standard profile overburden to the boundary line of the bedrock overburden, the relative elevation of the point to the bedrock, and the depth of the overburden are extracted based on the standard profile. The relative elevation of the point to the bedrock is the relative elevation of the boundary point between the spatial point on the top surface of the profile overburden and the spatial point on the top surface of the standard profile overburden with the minimum horizontal distance from the boundary line of the bedrock overburden. A polynomial model is constructed with the minimum horizontal distance from the spatial point on the top surface of the standard profile overburden to the boundary line of the bedrock overburden and the relative elevation of the point to the bedrock as input, and the overburden depth as output. By interpolating the nodes at the bottom of the deep overburden layer using a polynomial model, spatial information of the nodes at the bottom of the overburden layer is obtained, and a three-dimensional geological model of the overburden layer is established. The extraction of the minimum horizontal distance from the spatial point on the top surface of the standard profile overburden layer to the boundary line of the bedrock overburden layer, the relative elevation of the point to the bedrock surface, and the depth of the overburden layer based on the standard profile specifically includes: Using the formula Dmin Determine the minimum horizontal distance from the spatial point on the top surface of the standard profile overburden layer to the boundary line of the surface bedrock overburden layer; Using formula Determine the relative elevation of the point from the bedrock surface; The depth of the overburden layer is determined using the formula D = Pz - Pbz; Where Dmin represents the minimum horizontal distance from the top node of the overburden layer to the boundary line between the bedrock and the overburden layer; Let x be the x-coordinate of the i-th node on the boundary line between the bedrock overburden and the surface. Let Py be the y-coordinate of the i-th node on the boundary line between the bedrock and the surface overburden, n be the total number of boundary points constituting the boundary, Px be the ordinate of the spatial position of the node on the top surface of the overburden, Px be the abscissa of the spatial position of the node on the top surface of the overburden, and H be the relative elevation of the point from the bedrock. The elevation of the nodes on the top surface of the cover layer. Where is the elevation of the boundary node, D is the depth of the overburden layer, and Pbz is the elevation of the node at the bottom of the overburden layer. The polynomial model is: ； Where x is the input sample, For the input sample coefficients, m is determined by the maximum power of the polynomial. This is the predicted value from the regression model, and the predicted value is the cover depth.

2. The three-dimensional geological modeling method for overburden layers according to claim 1, characterized in that, The determination of standard profiles based on the distribution of hydropower engineering exploration points and geotechnical engineering investigation information specifically includes: AutoCAD was used to number the exploration points for hydropower projects and to lay out multiple standard profiles; Number the standard sections of the layout; Preprocessing of geotechnical engineering investigation information; The three-dimensional topographic surface is determined based on the pre-processed geotechnical engineering investigation information; The standard profile overburden layer is determined based on the three-dimensional topographic surface and the standard profile.

3. The three-dimensional geological modeling method for overburden layers according to claim 2, characterized in that, The preprocessing of geotechnical engineering investigation information specifically includes: Use the CGCS2000 coordinate system to unify the processing of geotechnical engineering investigation information; Use AutoCAD to extract elevation points, primary curves, and calculation curves of the terrain within the modeling area, and delete outliers; Convert the format of the borehole data.

4. The three-dimensional geological modeling method for overburden layers according to claim 3, characterized in that, The process of determining the three-dimensional topographic surface based on the preprocessed geotechnical engineering investigation information specifically includes: Use ItasCAD to generate surfaces from the modeling area; Use ItasCAD to generate the first elevation point file from the calculation curve and the primary curve data; Use ItasCAD to generate a second elevation point file from the elevation point data; Determine the elevation point file based on the first elevation point file and the second elevation point file; The 3D terrain surface is determined by applying fuzzy constraints to the surface generated within the modeling range based on the elevation point file.

5. A three-dimensional geological modeling system for overburden layers, used to implement the three-dimensional geological modeling method for overburden layers as described in any one of claims 1-4, characterized in that, include: The standard profile determination module is used to determine the standard profile based on the distribution of hydropower engineering exploration points and geotechnical engineering investigation information. The geotechnical engineering investigation information includes: borehole data, elevation information, and modeling scope; The standard profile extraction module is used to extract the minimum horizontal distance between the spatial points on the top surface of the standard profile overburden and the boundary line of the bedrock overburden, the relative elevation of the points from the bedrock, and the depth of the overburden based on the standard profile. The relative elevation of the points from the bedrock is the relative elevation of the boundary point between the spatial points on the top surface of the profile overburden and the spatial points on the top surface of the standard profile overburden and the boundary line of the bedrock overburden. The polynomial model building module is used to build a polynomial model with the minimum horizontal distance from the spatial point on the top surface of the standard profile overburden to the boundary line of the bedrock overburden and the relative elevation of the point to the bedrock as input, and the overburden depth as output. The overburden 3D geological model construction module is used to interpolate the nodes of the bottom surface of the deep overburden using a polynomial model, obtain the spatial information of the nodes of the bottom surface of the overburden, and establish a 3D geological model of the overburden.

6. A storage medium storing computer program instructions thereon, characterized in that, When the computer program instructions are executed by the processor, a three-dimensional geological modeling method for overburden layers as described in any one of claims 1-4 is implemented.

Images