A 3D geological modeling method for overburden strata with arbitrary constraints

Abstract

The present invention relates to a three-dimensional geological modeling method for a covering layer stratum to which arbitrary constraints can be added, the method comprising: obtaining a reference surface model of a site covering layer, and determining the outcrop line type and boundary range of the site covering layer. Based on the outcrop line type and boundary range of the site covering layer, the geological information of the site covering layer is selected, and the geological information includes geological points and boreholes of the site covering layer. Auxiliary lines are introduced and outlined into the site covering layer, and the corresponding positions of the auxiliary lines are given inferred thicknesses, and the calculation direction of the thickness is defined to infer the thickness distribution of the covering layer at a typical position. At the same time, the exploration points used in the process of modeling the site covering layer are selected as precise constraints, and the iterative calculation steps of the discrete smooth interpolation approximation operation are set, and the covering layer model is generated by using geometric interpolation and iterative operations. The problem of interlacing the surface and unreasonable thickness restoration during the modeling of the covering layer is solved, and the low efficiency of three-dimensional geological modeling caused by repeated adjustments and modifications is avoided.

Description

Technical Field

[0001] The invention relates to the technical field of geological modeling, and in particular to a three-dimensional geological modeling method for an overburden stratum capable of adding arbitrary constraints. Background Art

[0002] Traditional two-dimensional geological models are difficult to accurately reflect the complex underground geological conditions, while three-dimensional geological models can more comprehensively display underground geological structures, stratigraphic distribution, fault zones and other information, which can help engineering designers more accurately assess geological risks and underground engineering conditions, provide a reliable basis for geotechnical engineering design and construction, and improve the safety and efficiency of the project. As an important part of three-dimensional geological modeling, overburden modeling can quickly and accurately establish a site overburden stratigraphic model, assisting engineers to more accurately understand the properties, thickness, distribution and other information of the overburden, which can be used in various aspects such as geological disaster assessment, resource development and environmental protection, and engineering construction plan optimization. However, the overburden has the characteristics of a wide distribution area, multiple types, irregular boundaries, and little constraint data.

[0003] The objects of overburden modeling include overburden with outcrops of bedrock materials (such as residual layers) and overburden with outcrops of Quaternary materials (such as landslide deposits). A common problem in the modeling process is the lack of depth control information. A large amount of manual intervention is required randomly according to the rationality of the results during the modeling process, and it is difficult to perform batch processing like soil layer modeling. Therefore, there is an urgent need for a modeling method that can add arbitrary constraints to assist in establishing overburden strata.

[0004] At present, a method, system and electronic device for correcting dense surface grids under sparse constraint conditions are disclosed in the prior art. The method draws a closed line based on geological outcrop points and the estimated stratum thickness distribution range and projects it onto the terrain surface grid, then adds the estimated thickness constraint distribution data to the dense surface grid, and iterates and corrects the dense surface grid through geometric interpolation, so that the exposed boundary at the dense surface grid is consistent with the line connecting the outcrop points. However, the method is difficult to control the bottom plate morphology of the overburden layer generated when solving the geological surface modeling of the overburden layer. Since the isopach lines are the estimated results, the rationality of the bottom plate morphology may be affected if the estimation is inappropriate. In addition, the technicians cannot add auxiliary lines at any position of the bottom plate, and thus cannot use thickness data such as geological profile lines as thickness constraints. The terrain surface grid generated by the method after the thickness distribution interpolation and geometric iteration of the overburden layer is uneven, and the thickness distribution between the terrain surfaces cannot well meet the estimation of the geologists. After the modeling is completed, it is generally necessary to check the morphology of each overburden layer bottom plate, so a large amount of manual intervention is required for adjustment. Summary of the invention

[0005] Based on this, it is necessary to provide a three-dimensional geological modeling method for overburden strata to which arbitrary constraints can be added in order to address the above technical issues.

[0006] The present invention provides a three-dimensional geological modeling method for overburden strata capable of adding arbitrary constraints, the method comprising:

[0007] Obtaining a reference surface model of a site cover layer, and determining the outcrop line type and boundary range of the site cover layer;

[0008] Selecting geological information of the site covering layer based on the outcrop line type and boundary range of the site covering layer, the geological information including geological points and boreholes of the site covering layer;

[0009] Introducing and outlining auxiliary lines into the site covering layer, assigning estimated thicknesses to positions corresponding to the auxiliary lines, defining calculation directions of the thicknesses, and estimating thickness distribution of the covering layer at typical positions;

[0010] Set the number of iterative calculation steps for discrete smooth interpolation approximation operation, check the exploration points used in the site cover modeling process as precise constraints, and generate the cover model using geometric interpolation and iterative operations.

[0011] Furthermore, the obtaining of a reference surface model of the site covering layer and determining the outcrop line type and boundary range of the site covering layer includes:

[0012] Acquire a ground surface model of the site cover layer, and determine the outcrop line type and boundary range of the site cover layer to import the geological outcrop line;

[0013] Drawing a closed line according to the geological outcrop points to form the geological outcrop line, and performing encryption and smoothing processing on the geological outcrop line;

[0014] The processed geological outcrop points and geological outcrop lines are projected onto the ground surface, and all the geological outcrop points are defined as the precise constraints.

[0015] Furthermore, the method of obtaining a reference surface model of the site covering layer and determining the outcrop line type and boundary range of the site covering layer further includes:

[0016] The ground surface model is cut and trimmed according to the projected geological outcrop points and geological outcrop lines to delete the initial mesh surface formed by the geological outcrop points and the parts other than the closed lines in the geological outcrop lines;

[0017] The ground surface grid within the closed geological outcrop line is backed up, and the ground surface grid is the initial grid surface used to adjust the bottom plate surface morphology of the site covering layer.

[0018] Furthermore, the method further comprises:

[0019] When the ground surface grid inside the closed geological outcrop line is a cover layer outcrop, the initial grid surface inside the closed line is retained as the initial grid surface of the bottom plate of the site cover layer;

[0020] When the covering layer type corresponding to the ground surface grid inside the closed geological outcrop line is bedrock outcrop, the ground surface grid outside the closed geological outcrop line is used as the initial grid surface of the bottom plate surface of the site covering layer.

[0021] Furthermore, based on the outcrop line type and boundary range of the site covering layer, the geological information of the site covering layer is selected, and the geological information includes geological points and boreholes of the site covering layer, including:

[0022] Constraints are set on the initial grid surface of the site cover bottom surface morphology, the constraints include the precisely constrained cover outcrop line and the geological information, the geological information also includes the cover bottom position depth value and the geological boundary points determined by the geological points.

[0023] Furthermore, the method of introducing and outlining auxiliary lines into the site covering layer, assigning estimated thicknesses to positions corresponding to the auxiliary lines, and defining calculation directions of the thicknesses to estimate thickness distribution of the covering layer at typical positions includes:

[0024] Import the existing geological profile bottom line, according to the thickness of the site cover layer and the development law of the cover layer revealed by drilling, and introduce the auxiliary line by defining the thickness range and representative position of the cover layer to define the estimated thickness value of the cover layer;

[0025] The information of the depth of the site covering layer relative to the reference surface is obtained as a constraint condition for controlling the bottom surface morphology of the covering layer, and the inferred thickness value difference is distributed to each node of the ground surface grid within the modeling range through spatial data interpolation technology.

[0026] Furthermore, the iterative calculation steps of the discrete smooth interpolation approximation operation are set, and the exploration points used in the site cover modeling process are checked as precise constraints, and the cover model is generated by using geometric interpolation and iterative operation, including:

[0027] The estimated thickness value corresponding to the local area of ​​the site cover layer is extended to the entire space range for interpolation operation by a discrete smooth interpolation algorithm, so as to generalize the sparse constraint condition to each node on the initial grid surface;

[0028] Verify whether a given node attribute meets known constraint conditions based on the grid nodes on the initial grid surface, and solve the attribute value at the given node position point in a recursive manner through a linear equation.

[0029] The present invention also provides a three-dimensional geological modeling device for overburden strata capable of adding arbitrary constraints, the device comprising:

[0030] A covering layer determination module is used to obtain a reference surface model of the site covering layer and determine the outcrop line type and boundary range of the site covering layer;

[0031] An information selection module, used for selecting geological information of the site covering layer based on the outcrop line type and boundary range of the site covering layer, wherein the geological information includes geological points and boreholes of the site covering layer;

[0032] An auxiliary inference module is used to introduce and outline auxiliary lines into the site covering layer, assign inferred thickness to the positions corresponding to the auxiliary lines, and define the calculation direction of the thickness to infer the thickness distribution of the covering layer at a typical position;

[0033] The model generation module is used to set the iterative calculation steps of the discrete smooth interpolation approximation operation, and to check the exploration points used in the site cover modeling process as precise constraints, and to generate the cover model using geometric interpolation and iterative operations.

[0034] The present invention also provides an electronic device, comprising a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, it implements a three-dimensional geological modeling method for overburden strata to which arbitrary constraints can be added as described in any one of the above.

[0035] The present invention also provides a computer storage medium storing a computer program, wherein when the computer program is executed by a processor, the computer program implements a three-dimensional geological modeling method for overburden strata to which arbitrary constraints can be added as described in any one of the above.

[0036] The present invention also provides a computer program product, comprising a computer program, wherein when the computer program is executed by a processor, the computer program implements any of the above-mentioned three-dimensional geological modeling methods for overburden strata to which arbitrary constraints can be added.

[0037] The above-mentioned three-dimensional geological modeling method of the cover layer stratum to which arbitrary constraints can be added obtains the reference surface model of the site cover layer and determines the outcrop line type and boundary range of the site cover layer. Then, based on the outcrop line type and boundary range of the site cover layer, the geological information of the site cover layer is selected, and the geological information includes geological points and boreholes of the site cover layer. Subsequently, auxiliary lines are introduced and outlined into the site cover layer, and the estimated thickness is assigned to the corresponding position of the auxiliary line, and the calculation direction of the thickness is defined to infer the thickness distribution of the cover layer at a typical position. Finally, the iterative calculation steps of the discrete smooth interpolation approximation operation are set, and the exploration points used in the site cover layer modeling process are checked as precise constraints, and the cover layer model is generated using geometric interpolation and iterative operations. This method uses geometric information, namely the results of surface outcrop surveys and exploration results, to provide a fast and effective way to integrate various geological surveying and mapping data such as drilling, and supports the addition of arbitrary constraints to model the overburden. The modeling process follows a series of set basic principles, which fundamentally guarantees the geological rationality of the model and solves the problems of interlacing the surface and unreasonable thickness restoration when modeling the overburden, while avoiding the pain points of low efficiency of three-dimensional geological modeling due to repeated adjustments and modifications. BRIEF DESCRIPTION OF THE DRAWINGS

[0038] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0039] Figure 1 One of the flow diagrams of a method for three-dimensional geological modeling of overburden strata to which arbitrary constraints can be added provided by the present invention;

[0040] Figure 2 A second flow chart of a method for three-dimensional geological modeling of overburden strata to which arbitrary constraints can be added provided by the present invention;

[0041] Figure 3 A third flow chart of a method for three-dimensional geological modeling of overburden strata capable of adding arbitrary constraints provided by the present invention;

[0042] Figure 4 A fourth flow chart of a method for three-dimensional geological modeling of overburden strata capable of adding arbitrary constraints provided by the present invention;

[0043] Figure 5 A fifth flow chart of a method for three-dimensional geological modeling of overburden strata capable of adding arbitrary constraints provided by the present invention;

[0044] Figure 6A sixth flow chart of a method for three-dimensional geological modeling of overburden strata capable of adding arbitrary constraints provided by the present invention;

[0045] Figure 7 A schematic diagram of the overall process of three-dimensional geological modeling of overburden layers in a three-dimensional geological modeling method of overburden layers capable of adding arbitrary constraints in a specific embodiment provided by the present invention;

[0046] Figure 8 A schematic diagram of the structure of a three-dimensional geological modeling device for overburden strata to which arbitrary constraints can be added provided by the present invention;

[0047] Fig. 9 This is a diagram of the internal structure of the computer device provided by the present invention. DETAILED DESCRIPTION

[0048] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0049] Combine the following Figure 1-Figure 9 The present invention describes a three-dimensional geological modeling method for overburden strata capable of adding arbitrary constraints.

[0050] like Figure 1 As shown, in one embodiment, a method for three-dimensional geological modeling of overburden strata capable of adding arbitrary constraints comprises the following steps:

[0051] Step S110, obtaining a reference surface model of the site covering layer, and determining the outcrop line type and boundary range of the site covering layer.

[0052] Step S120 , based on the outcrop line type and boundary range of the site covering layer, geological information of the site covering layer is selected, the geological information including geological points and boreholes of the site covering layer.

[0053] Specifically, constraints are set for the initial grid surface of the determined overburden bottom surface morphology, the constraints include precisely constrained overburden outcrop lines and geological information, the geological information also includes the overburden bottom position depth value and geological boundary points determined by geological points.

[0054] Step S130, introducing and outlining auxiliary lines into the site covering layer, assigning estimated thicknesses to positions corresponding to the auxiliary lines, defining thickness calculation directions, and estimating thickness distribution of the covering layer at typical positions.

[0055] Step S140, setting the iterative calculation steps of the discrete smooth interpolation approximation operation, and checking the exploration points used in the site cover modeling process as precise constraints, and generating the cover model by using geometric interpolation and iterative operation.

[0056] The above-mentioned three-dimensional geological modeling method of the cover layer strata to which arbitrary constraints can be added obtains the reference surface model of the site cover layer and determines the outcrop line type and boundary range of the site cover layer. Based on the outcrop line type and boundary range of the site cover layer, the geological information of the site cover layer is selected, and the geological information includes geological points and boreholes of the site cover layer. Subsequently, auxiliary lines are introduced and outlined into the site cover layer, and the estimated thickness is assigned to the corresponding position of the auxiliary line, and the calculation direction of the thickness is defined to infer the thickness distribution of the cover layer at a typical position. Finally, the iterative calculation steps of the discrete smooth interpolation approximation operation are set, and the exploration points used in the modeling process of the site cover layer are checked as precise constraints, and the cover layer model is generated using geometric interpolation and iterative operations. This method uses geometric information, namely the results of surface outcrop surveys and exploration results, to provide a fast and effective way to integrate various geological surveying and mapping data such as drilling, and supports the addition of arbitrary constraints to model the overburden. The modeling process follows a series of set basic principles, which fundamentally guarantees the geological rationality of the model and solves the problems of interlacing the surface and unreasonable thickness restoration when modeling the overburden, while avoiding the pain points of low efficiency of three-dimensional geological modeling due to repeated adjustments and modifications.

[0057] like Figure 2 As shown, in some embodiments, the present invention provides a three-dimensional geological modeling method for covering layer strata to which arbitrary constraints can be added, obtains a reference surface model of the site covering layer, and determines the outcrop line type and boundary range of the site covering layer, and specifically includes the following steps:

[0058] Step S111, obtaining a ground surface model of the site cover layer, and determining the outcrop line type and boundary range of the site cover layer to import the geological outcrop line.

[0059] Step S112, drawing closed lines according to geological outcrop points to form geological outcrop lines, and performing encryption and smoothing processing on the geological outcrop lines.

[0060] Step S113, projecting the processed geological outcrop points and geological outcrop lines onto the ground surface, and defining all geological outcrop points as precise constraints.

[0061] like Figure 3 As shown, in some embodiments, the present invention provides a method for three-dimensional geological modeling of overburden strata to which arbitrary constraints can be added, further comprising the following steps:

[0062] Step S310, cutting and trimming the ground surface model according to the projected geological outcrop points and geological outcrop lines, so as to delete the initial mesh surface formed by the geological outcrop points and the parts of the geological outcrop lines other than the closed lines.

[0063] Step S320, backing up the ground surface grid within the closed geological outcrop line, the ground surface grid is the initial grid surface used to adjust the bottom plate surface morphology of the site cover layer.

[0064] like Figure 4 As shown, in some embodiments, the present invention provides a method for three-dimensional geological modeling of overburden strata to which arbitrary constraints can be added, further comprising the following steps:

[0065] Step S410: When the ground surface grid inside the closed geological outcrop line is a cover layer outcrop, the initial grid surface inside the closed line is retained as the initial grid surface of the site cover layer bottom surface morphology.

[0066] Step S420: When the type of the covering layer corresponding to the ground surface grid inside the closed geological outcrop line is bedrock outcrop, the ground surface grid outside the closed geological outcrop line is used as the initial grid surface of the bottom plate surface of the site covering layer.

[0067] like Figure 5 As shown, in some embodiments, the present invention provides a three-dimensional geological modeling method for covering layer strata with arbitrary constraints, which introduces and outlines auxiliary lines into the site covering layer, assigns inferred thickness to the corresponding positions of the auxiliary lines, and defines the calculation direction of the thickness to infer the thickness distribution of the covering layer at a typical position, specifically including the following steps:

[0068] Step S131, import the existing geological profile bottom line, according to the thickness of the site cover layer revealed by drilling and the law of cover layer development, and introduce auxiliary lines by defining the cover layer thickness range and representative positions to define the estimated thickness value of the cover layer.

[0069] Step S132, obtaining information on the depth of the site cover relative to the reference surface as a constraint condition for controlling the bottom surface morphology of the cover, and distributing the estimated thickness value difference to each node of the ground surface grid within the modeling range through spatial data interpolation technology.

[0070] like Figure 6 As shown, in some embodiments, the present invention provides a three-dimensional geological modeling method for cover layer strata to which arbitrary constraints can be added, sets the iterative calculation steps of discrete smooth interpolation approximation operation, and selects the exploration points used in the process of site cover layer modeling as precise constraints, and generates a cover layer model by using geometric interpolation and iterative operation, which specifically includes the following steps:

[0071] Step S141, using a discrete smooth interpolation algorithm, the estimated thickness value corresponding to the local site coverage layer is expanded to the entire space range for interpolation operation, so as to generalize the sparse constraint condition to each node on the initial grid surface.

[0072] Step S142, verifying whether the given node attribute meets the known constraint conditions according to the grid nodes on the initial grid surface, and solving the attribute value at the given node position point in a recursive manner through a linear equation.

[0073] Combination Figure 7 As shown, in a specific embodiment, the present invention provides a three-dimensional geological modeling method for overburden strata to which arbitrary constraints can be added. First, a ground surface model of the site overburden is obtained, the outcrop line type and boundary range of the overburden are determined, and geological outcrop lines are imported, or closed lines are drawn according to geological outcrop points to form geological outcrop lines, the geological outcrop lines are encrypted and smoothed, and then the processed geological outcrop points and geological outcrop lines are projected onto the ground surface, and all geological outcrop points are defined as precise constraints.

[0074] Secondly, the computer program backs up the ground surface grid within the closed geological outcrop line in the background, and the ground surface grid will be used as the initial grid surface for adjusting the bottom plate shape of the overburden layer. The grid surface is the initial dense grid surface formed by the computer intersecting and cutting the terrain surface model according to the projected geological outcrop line and deleting the part outside the closed line.

[0075] In addition, the dense mesh surface needs to be processed separately according to the type of cover layer; when the mesh inside the closed geological outcrop line represents the cover layer outcrop, the dense mesh inside the closed line is retained as the initial dense mesh surface of the cover layer bottom surface morphology. When the cover layer type of the mesh inside the closed geological outcrop line is bedrock outcrop, the ground surface mesh outside the geological outcrop closed line is used as the initial dense mesh surface of the cover layer bottom surface morphology.

[0076] Thirdly, various constraints are set for the determined initial dense grid of the overburden bottom surface, including the outcrop line of the overburden and exploration data as precise constraints. The exploration data include the depth value of the overburden bottom position revealed by drilling holes and flat holes, and the geological boundary points determined by geological points, so that 100% of the triangulated network nodes constituting the overburden interface pass through the point positions under precise constraints, thereby more accurately controlling the morphology of the overburden bottom.

[0077] After that, the existing inferred data (bottom line of geological profile) is imported, or according to the thickness of the overburden revealed by drilling and the development law of the overburden, the inferred thickness value of the overburden is defined by defining the overburden thickness range and introducing auxiliary lines (bottom line of geological profile) at representative positions, and the information of the depth of the overburden relative to the reference surface (surface) is obtained as a constraint condition for controlling the morphology of the overburden floor. The former ensures that the thickness inferred value is randomly distributed within the overburden modeling range by setting local inferred thickness points, and the latter assigns the thickness inferred value to the local position of the overburden floor by importing or drawing auxiliary lines. The thickness value is interpolated and distributed to each grid node within the modeling range through spatial data interpolation technology. Before interpolation, the thickness calculation direction of the overburden floor relative to the reference surface needs to be defined, and the interpolation operation is generally performed in the vertical direction by default.

[0078] It should be noted that when the bottom boundary of the geological profile is set as a boundary constraint, the bottom plate morphology of the overburden layer can be controlled more accurately and an overburden layer model that conforms to geological laws can be established. If the bottom boundary of the geological profile is not introduced as a constraint, the dense grid of the overburden layer needs to be locally adjusted and corrected, and constraints need to be added multiple times for interpolation iteration to achieve the expected shape.

[0079] In this embodiment, the spatial interpolation technology used in the interpolation of the overburden modeling, namely DSI (discrete smooth interpolation algorithm), can fit and construct complex geological models in real time according to various constraints. The geological rationality of the overburden thickness at a typical location is achieved by introducing auxiliary lines, assigning the estimated thickness to the auxiliary line parts, and using the DSI interpolation algorithm to extend the local estimated values ​​to the entire space range for interpolation operations by carrying the relevant attributes of the local estimated thickness points on the auxiliary lines. The algorithm establishes mutually connected grids for a discretized natural body model, and generalizes the sparse constraints to each node of the dense surface grid, so that each node of the initial dense surface grid can obtain the thickness constraints relative to the terrain surface, and verify whether the given node attributes meet the known constraints according to the network nodes, and then the attribute values ​​at the location points can be solved recursively through linear equations, and the interpolation accuracy of the created surface can be adjusted by setting the number of recursive steps.

[0080] Finally, according to the geological outcrop line, precisely constrained exploration data, local estimated thickness values ​​of typical parts and the bottom boundary of the geological profile that controls the bottom plate morphology of the overburden, the iterative calculation steps of the DSI (discrete smooth interpolation) approximation operation are set to achieve adjustment and correction of the bottom plate of the overburden and generate a three-dimensional model of the overburden.

[0081] It should be noted that the overburden is a common geological body with a semi-closed surface as its boundary. The morphology of the overburden bottom plate often has certain similarities with the surface. When modeling the overburden bottom plate, it is often necessary to refer to the surface morphology. Therefore, it is necessary to import the overburden outcrop line and terrain surface model.

[0082] In this embodiment, the geological modeling rarely uses the method of connecting known points to construct lines or surfaces. Instead, the known points are regarded as constraints, and the role of constraints is reflected by continuously approaching or completely passing through. Since the known points are generally sparse, the control of the shape and rationality between the known points is often reflected through constraints. In many cases, constraints are the most common and effective means to achieve geological speculation. The constraints here refer to precise constraints and fuzzy constraints.

[0083] In this embodiment, the outcrop line of the overburden will be used as a precise constraint, that is, the shape and position of the outcrop line in three-dimensional space will not change with the change of other constraints (fuzzy constraints generated by inferred data and DSI interpolation). At the same time, the exploration data will also be used as a precise constraint, including the depth value of the overburden floor position revealed by drilling and flat holes, and the geological boundary points determined by geological points, so that the triangulated network nodes constituting the overburden interface pass through the point positions under the precise constraints 100%, and the shape of the overburden floor is controlled more accurately. Then, the existing inferred data (bottom boundary line of the geological profile) is imported, or according to the thickness of the overburden revealed by drilling and the development law of the overburden, auxiliary lines are drawn on the terrain surface model, and the thickness inferred value is assigned. The computer projects the nodes on these auxiliary lines vertically on the terrain surface in the background, and uses the thickness inferred value as the vertical downward depth, so that the triangulated network nodes constituting the overburden interface are close to the nodes after the downward movement, that is, as the fuzzy constraint condition for controlling the shape of the overburden floor, the shape of the overburden floor is controlled.

[0084] In this embodiment, the terrain surface model is copied, and then the copied model is cut with the outcrop line to delete unnecessary parts (when the outcrop line is a landslide accumulation, the terrain surface grid within the outcrop line boundary is used as the modeling range of the cover layer, and the terrain surface grid outside the outcrop line boundary is the unnecessary part. When the outcrop line is a bedrock outcrop, the terrain surface grid outside the outcrop line boundary range is used as the modeling range of the cover layer, and the terrain surface grid within the outcrop line boundary range is the unnecessary part. When the modeling range extends beyond the outcrop line boundary range, the extended segment uses the terrain surface boundary as the outer boundary of the cover layer modeling range, and the terrain surface grid within the outcrop line boundary range is the unnecessary part.), and the retained part is subjected to the overburden interface morphology correction using the overburden interface thickness exposed by the drill hole and the flat hole and the inferred thickness data as constraints.

[0085] In this embodiment, constraint data is added to the overburden floor model, namely, exploration data (overburden interface in boreholes and flat holes) that reveals the position of the overburden floor, and real and reliable depth data from the terrain surface obtained through exploration, that is, overburden thickness data. The exploration interface provides precise constraints for the overburden (precise constraints refer to fixing the spatial position of the overburden interface revealed by the boreholes and flat holes and preventing it from changing position with subsequent constraint adjustments). These data can be used in the modeling process to more accurately control the morphology of the overburden floor.

[0086] In this embodiment, arbitrary constraints are added: the estimated thickness value of the overburden is defined by defining the overburden thickness range and introducing auxiliary lines (bottom boundary of geological profile) at representative positions, and obtaining information on the depth of the overburden relative to the reference surface (surface) as a constraint to control the morphology of the overburden floor. The former ensures that the estimated thickness value is randomly distributed within the overburden modeling range by setting local estimated thickness points, and the latter assigns the estimated thickness value to the local position of the overburden floor by importing or drawing auxiliary lines. The thickness value is interpolated and distributed to each grid node within the modeling range through spatial data interpolation technology. Before interpolation, the thickness calculation direction of the overburden floor relative to the reference surface needs to be defined, and the interpolation operation is generally performed in the vertical direction by default.

[0087] It should be noted that the basic idea of ​​adjusting and correcting the bottom plate morphology here is the DSI interpolation operation after constraint processing, that is, the joint application of constraint data and DSI discrete smooth interpolation algorithm.

[0088] There are two types of constraints that geologists often use during modeling: vector constraints and thickness constraints. Vector constraints use the occurrence of the exploration point as an element to "guess" the coordinates of unknown points, correctly reflecting the undulating shape; thickness constraints are mostly used to agree on the relationship between the upper and lower interfaces of the same geological body. The thickness is not less than zero (no interpenetration) is the most basic constraint condition (pinch-out when it is zero), which effectively avoids the interpenetration phenomenon that may occur in the mathematical calculation process. When mathematical language is used to describe geological laws and converted into constraints for the interpolation process, the reliability of the results and work efficiency will be greatly improved.

[0089] In this embodiment, in addition to the constraints, the spatial interpolation technology used in the overburden modeling, namely DSI (discrete smooth interpolation algorithm), can fit and construct complex geological models in real time according to various constraints. The geological rationality of the overburden thickness at a typical location is achieved by introducing auxiliary lines, assigning the estimated thickness to the auxiliary line parts, and using the DSI interpolation algorithm to extend the local estimated values ​​to the entire space range for interpolation operations by carrying the relevant attributes of the local estimated thickness points on the auxiliary lines. The algorithm establishes a mutually connected grid for a discretized natural body model, and generalizes the sparse constraints to each node of the dense surface grid, so that each node of the initial dense surface grid can obtain the thickness constraint relative to the terrain surface, and verify whether the given node attribute meets the known constraint according to the network node, and then the attribute value at the position point can be solved recursively through the linear equation, and the interpolation accuracy of the created surface can be adjusted by setting the number of recursive steps.

[0090] Based on the geological outcrop line, precisely constrained exploration data, local estimated thickness values ​​of typical parts and the bottom boundary of the geological section to control the bottom plate morphology of the overburden, the DSI interpolation algorithm is used to perform interpolation operations to adjust and correct the bottom plate of the overburden and generate a three-dimensional model of the overburden.

[0091] The following is a description of a three-dimensional geological modeling device for an overburden layer stratum to which arbitrary constraints can be added provided by the present invention. The three-dimensional geological modeling device for an overburden layer stratum to which arbitrary constraints can be added described below and the three-dimensional geological modeling method for an overburden layer stratum to which arbitrary constraints can be added described above can refer to each other.

[0092] like Figure 8 As shown, in one embodiment, a three-dimensional geological modeling device for overburden strata to which arbitrary constraints can be added includes an overburden determination module 810 , an information selection module 820 , an auxiliary inference module 830 and a model generation module 840 .

[0093] The cover layer determination module 810 is used to obtain a reference surface model of the site cover layer and determine the outcrop line type and boundary range of the site cover layer.

[0094] The information selection module 820 is used to select geological information of the site covering layer based on the outcrop line type and boundary range of the site covering layer, and the geological information includes geological points and boreholes of the site covering layer.

[0095] The auxiliary inference module 830 is used to introduce and outline auxiliary lines into the site covering layer, assign inferred thickness to the positions corresponding to the auxiliary lines, and define the calculation direction of the thickness to infer the thickness distribution of the covering layer at typical positions.

[0096] The model generation module 840 is used to set the iterative calculation steps of the discrete smooth interpolation approximation operation, and to select the exploration points used in the site cover modeling process as precise constraints, and to generate the cover model using geometric interpolation and iterative operations.

[0097] In this embodiment, the present invention provides a three-dimensional geological modeling device for overburden strata capable of adding arbitrary constraints, and the overburden determination module is specifically used for:

[0098] Obtain the ground surface model of the site cover, and determine the outcrop line type and boundary range of the site cover to import the geological outcrop line.

[0099] According to the geological outcrop points, closed lines are drawn to form geological outcrop lines, and the geological outcrop lines are encrypted and smoothed.

[0100] The processed geological outcrop points and geological outcrop lines are projected onto the ground surface, and all geological outcrop points are defined as precise constraints.

[0101] In this embodiment, the present invention provides a three-dimensional geological modeling device for overburden strata capable of adding arbitrary constraints, further comprising a grid surface clipping module for:

[0102] The ground surface model is intersected and trimmed according to the projected geological outcrop points and geological outcrop lines to delete the initial mesh surface formed by the parts outside the closed lines of the geological outcrop points and geological outcrop lines.

[0103] The ground surface grid within the closed geological outcrop line is backed up, and the ground surface grid is the initial grid surface used to adjust the bottom plate surface shape of the site cover layer.

[0104] In this embodiment, the present invention provides a three-dimensional geological modeling device for covering layer strata to which arbitrary constraints can be added, and further includes a covering surface type screening module for:

[0105] When the ground surface grid inside the closed geological outcrop line is a cover layer outcrop, the initial grid surface inside the closed line is retained as the initial grid surface of the site cover layer bottom plate surface morphology.

[0106] When the cover type corresponding to the ground surface grid inside the closed geological outcrop line is bedrock outcrop, the ground surface grid outside the closed geological outcrop line is used as the initial grid surface of the bottom plate surface of the site cover.

[0107] In this embodiment, the present invention provides a three-dimensional geological modeling device for overburden strata capable of adding arbitrary constraints, and the information selection module is specifically used for:

[0108] Constraints are set on the initial grid surface of the overburden bottom surface morphology. The constraints include precisely constrained overburden outcrop lines and geological information. The geological information also includes the overburden bottom position depth value and geological boundary points determined by geological points.

[0109] In this embodiment, the present invention provides a three-dimensional geological modeling device for overburden strata capable of adding arbitrary constraints, wherein the auxiliary inference module is specifically used for:

[0110] Import the existing geological profile bottom boundary, based on the thickness of the site cover layer revealed by drilling and the law of cover layer development, and introduce auxiliary lines by defining the cover layer thickness range and representative positions to define the estimated thickness value of the cover layer.

[0111] The information on the depth of the site cover relative to the reference surface is obtained as a constraint condition for controlling the bottom surface morphology of the cover layer, and the inferred thickness value difference is distributed to each node of the ground surface grid within the modeling range through spatial data interpolation technology.

[0112] In this embodiment, the present invention provides a three-dimensional geological modeling device for overburden strata capable of adding arbitrary constraints, and the model generation module is specifically used for:

[0113] The inferred thickness value corresponding to the local site coverage layer is extended to the entire space for interpolation operation through the discrete smooth interpolation algorithm, so as to generalize the sparse constraint condition to each node on the initial grid surface.

[0114] Verify whether the given node attributes meet the known constraints based on the grid nodes on the initial grid surface, and solve the attribute value at the given node position point recursively through linear equations.

[0115] Fig. 9 An example of a physical structure diagram of an electronic device is shown. The electronic device may be a smart terminal, and its internal structure diagram may be as follows: Fig. 9 As shown. The electronic device includes a processor, a memory and a network interface connected through a system bus. Among them, the processor of the electronic device is used to provide computing and control capabilities. The memory of the electronic device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The network interface of the electronic device is used to communicate with an external terminal through a network connection. When the computer program is executed by the processor, a three-dimensional geological modeling method of an overburden stratum to which arbitrary constraints can be added is implemented, and the method includes:

[0116] Obtain the reference surface model of the site cover layer and determine the outcrop line type and boundary range of the site cover layer;

[0117] Based on the outcrop line type and boundary range of the site cover layer, the geological information of the site cover layer is selected, and the geological information includes geological points and drill holes of the site cover layer;

[0118] Introduce and draw auxiliary lines into the site cover layer, assign estimated thickness to the corresponding positions of the auxiliary lines, define the calculation direction of the thickness, and estimate the thickness distribution of the cover layer at typical locations;

[0119] Set the number of iterative calculation steps for discrete smooth interpolation approximation operations, check the exploration points used in the site cover modeling process as precise constraints, and use geometric interpolation and iterative operations to generate the cover model.

[0120] Those skilled in the art will understand that Fig. 9 The structure shown in the figure is only a block diagram of a part of the structure related to the scheme of the present invention, and does not constitute a limitation on the electronic device to which the scheme of the present invention is applied. The specific electronic device may include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components.

[0121] On the other hand, the present invention also provides a computer storage medium storing a computer program, which, when executed by a processor, implements a three-dimensional geological modeling method for overburden strata capable of adding arbitrary constraints, the method comprising:

[0122] Obtain the reference surface model of the site cover layer and determine the outcrop line type and boundary range of the site cover layer;

[0123] Based on the outcrop line type and boundary range of the site cover layer, the geological information of the site cover layer is selected, and the geological information includes geological points and drill holes of the site cover layer;

[0124] Introduce and draw auxiliary lines into the site cover layer, assign estimated thickness to the corresponding positions of the auxiliary lines, define the calculation direction of the thickness, and estimate the thickness distribution of the cover layer at typical locations;

[0125] Set the number of iterative calculation steps for discrete smooth interpolation approximation operations, check the exploration points used in the site cover modeling process as precise constraints, and use geometric interpolation and iterative operations to generate the cover model.

[0126] In another aspect, a computer program product or a computer program is provided, the computer program product or the computer program comprising computer instructions, the computer instructions being stored in a computer-readable storage medium. A processor of an electronic device reads the computer instructions from the computer-readable storage medium, and when the processor executes the computer instructions, a method for three-dimensional geological modeling of an overburden stratum capable of adding arbitrary constraints is implemented, the method comprising:

[0127] Obtain the reference surface model of the site cover layer and determine the outcrop line type and boundary range of the site cover layer;

[0128] Based on the outcrop line type and boundary range of the site cover layer, the geological information of the site cover layer is selected, and the geological information includes geological points and drill holes of the site cover layer;

[0129] Introduce and draw auxiliary lines into the site cover layer, assign estimated thickness to the corresponding positions of the auxiliary lines, define the calculation direction of the thickness, and estimate the thickness distribution of the cover layer at typical locations;

[0130] Set the number of iterative calculation steps for discrete smooth interpolation approximation operations, check the exploration points used in the site cover modeling process as precise constraints, and use geometric interpolation and iterative operations to generate the cover model.

[0131] Those of ordinary skill in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be implemented by instructing related hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided by the present invention may include non-volatile and / or volatile memory. Non-volatile memory may include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM) or flash memory. Volatile memory may include random access memory (RAM) or external cache memory.

[0132] By way of illustration and not limitation, RAM is available in many forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.

[0133] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0134] The above-mentioned embodiments only express several implementation methods of the present invention, and the description thereof is relatively specific and detailed, but it cannot be understood as limiting the scope of the patent of the present invention. It should be pointed out that, for ordinary technicians in this field, several variations and improvements can be made without departing from the concept of the present invention, which all belong to the protection scope of the present invention. Therefore, the protection scope of the patent of the present invention shall be subject to the attached claims.

Claims

1. A three-dimensional geological modeling method for overburden strata capable of adding arbitrary constraints, characterized in that: The method comprises: Obtaining a reference surface model of a site cover layer, and determining the outcrop line type and boundary range of the site cover layer; Selecting geological information of the site covering layer based on the outcrop line type and boundary range of the site covering layer, the geological information including geological points and boreholes of the site covering layer; Introducing and outlining auxiliary lines into the site covering layer, assigning estimated thicknesses to positions corresponding to the auxiliary lines, defining calculation directions of the thicknesses, and estimating thickness distribution of the covering layer at typical positions; Set the number of iterative calculation steps for discrete smooth interpolation approximation operation, check the exploration points used in the site cover modeling process as precise constraints, and use geometric interpolation and iterative operation to generate the cover model; The selecting of geological information of the site covering layer based on the outcrop line type and boundary range of the site covering layer, wherein the geological information includes geological points and boreholes of the site covering layer, comprises: Setting constraints on the initial grid surface of the site cover bottom plate surface morphology, the constraints include the precisely constrained cover outcrop line and the geological information, the geological information also includes the cover bottom plate position depth value and the geological boundary point determined by the geological point; The step of introducing and outlining auxiliary lines into the site covering layer, assigning estimated thicknesses to positions corresponding to the auxiliary lines, and defining a calculation direction of the thicknesses to estimate the thickness distribution of the covering layer at a typical position includes: Import the existing geological profile bottom line, according to the thickness of the site cover layer and the development law of the cover layer revealed by drilling, and introduce the auxiliary line by defining the thickness range and representative position of the cover layer to define the estimated thickness value of the cover layer; The information of the depth of the site covering layer relative to the reference surface is obtained as a constraint condition for controlling the bottom surface morphology of the covering layer, and the inferred thickness value difference is distributed to each node of the ground surface grid within the modeling range through spatial data interpolation technology.

2. The three-dimensional geological modeling method of overburden strata capable of adding arbitrary constraints according to claim 1, characterized in that: The step of obtaining a reference surface model of a site covering layer and determining an outcrop line type and a boundary range of the site covering layer includes: Acquire a ground surface model of the site cover layer, and determine the outcrop line type and boundary range of the site cover layer to import the geological outcrop line; Drawing a closed line according to the geological outcrop points to form the geological outcrop line, and performing encryption and smoothing processing on the geological outcrop line; The processed geological outcrop points and geological outcrop lines are projected onto the ground surface, and all the geological outcrop points are defined as the precise constraints.

3. The three-dimensional geological modeling method of overburden strata capable of adding arbitrary constraints according to claim 2, characterized in that: The step of obtaining a reference surface model of the site covering layer and determining the outcrop line type and boundary range of the site covering layer further includes: The ground surface model is cut and trimmed according to the projected geological outcrop points and geological outcrop lines to delete the initial mesh surface formed by the geological outcrop points and the parts other than the closed lines in the geological outcrop lines; The ground surface grid within the closed geological outcrop line is backed up, and the ground surface grid is the initial grid surface used to adjust the bottom plate surface morphology of the site covering layer.

4. The method for three-dimensional geological modeling of overburden strata capable of adding arbitrary constraints according to claim 3, characterized in that: The method further comprises: When the ground surface grid inside the closed geological outcrop line is a cover layer outcrop, the initial grid surface inside the closed line is retained as the initial grid surface of the bottom plate of the site cover layer; When the covering layer type corresponding to the ground surface grid inside the closed geological outcrop line is bedrock outcrop, the ground surface grid outside the closed geological outcrop line is used as the initial grid surface of the bottom plate surface of the site covering layer.

5. The three-dimensional geological modeling method of overburden strata capable of adding arbitrary constraints according to claim 1, characterized in that: The step of setting the iterative calculation steps of the discrete smooth interpolation approximation operation, and checking the exploration points used in the site cover modeling process as precise constraints, and using geometric interpolation and iterative operations to generate a cover model, includes: The estimated thickness value corresponding to the local area of ​​the site cover layer is extended to the entire space range for interpolation operation by a discrete smooth interpolation algorithm, so as to generalize the sparse constraint condition to each node on the initial grid surface; Verify whether a given node attribute meets known constraint conditions based on the grid nodes on the initial grid surface, and solve the attribute value at the given node position point in a recursive manner through a linear equation.

Images