A calculation method, system and medium for the thickness of any formation covering layer

Abstract

The present invention discloses a calculation method, system and medium for the thickness of the covering layer of any formation. The present invention includes establishing a three-dimensional geological model of the formation based on geological exploration data; importing the position coordinates of the location where the thickness of the formation covering layer needs to be calculated into the three-dimensional geological model, and calculating the projection distance from the position coordinates to the surface of the formation in the three-dimensional geological model as the thickness of the formation covering layer at the position coordinates. The present invention can quickly and accurately establish a three-dimensional geological model based on geological exploration data, with high precision and fast speed. Newly added borehole exploration points can be added in real time during the construction process, and the real-time model can be dynamically modified. The thickness of the formation at any unknown part can be calculated by interpolation analysis using the known exploration points. It has the advantages of fast calculation speed, high accuracy, intuitive graphics, accurate numerical values, and can avoid human calculation errors, which has great significance for the refined construction of projects such as rock foundation construction and pile foundation construction.

Description

Technical Field

[0001] The present invention relates to the field of geological engineering, and specifically relates to a method, system and medium for calculating the thickness of the covering layer of any formation, which can be applied to fields such as engineering exploration, pile foundation construction, and rock foundation construction. Background Art

[0002] During the construction of underground projects such as foundation treatment projects, geological exploration data is an important guiding document for underground design and construction, and its formation profile is one of its main contents. At present, domestic formation profiles usually adopt structural diagrams, and the data sources for their drawing mainly rely on actual drilling data. The thickness of the formation covering layer calculated by this method has high accuracy when the drilling spacing is small. However, when the drilling spacing is large, in order to ensure the accuracy of the data relied on by the design, the exploration boreholes are often encrypted, such as redrawing the formation profile, which is time-consuming and laborious. Moreover, the profile drawn by this means can only reflect the formation undulation changes between the drilling holes, and it is not convenient to calculate the formation conditions at positions outside the connection line of the exploration holes. Summary of the Invention

[0003] The technical problem to be solved by the present invention: Aiming at the above problems of the prior art, the present invention provides a method, system and medium for calculating the thickness of the covering layer of any formation. The present invention can quickly and accurately establish a three-dimensional geological model based on geological exploration data, with high precision and speed. Newly added drilling exploration points can be added in real time during the construction process, and the real-time model can be dynamically modified. The thickness of the formation at any unknown part is calculated by interpolation analysis using the known exploration points, which has the advantages of fast calculation speed, high accuracy, intuitive graphics, accurate numerical values, and can avoid human calculation errors, and has great significance for the refined construction of projects such as rock foundation construction and pile foundation construction.

[0004] In order to solve the above technical problems, the technical solution adopted by the present invention is as follows:

[0005] A method for calculating the thickness of the covering layer of any formation, comprising:

[0006] 1) Establish a three-dimensional geological model of the formation based on geological exploration data;

[0007] 2) Import the position coordinates where the thickness of the formation covering layer needs to be calculated into the three-dimensional geological model, and calculate the projection distance from the position coordinates to the formation surface in the three-dimensional geological model as the thickness of the formation covering layer at the position coordinates.

[0008] Optionally, step 1) includes:

[0009] 1.1) Extract the coordinates of each exploration hole and the surface elevation from the plane layout diagram of the geological exploration data. For each exploration hole, extract the elevation of each stratigraphic interface exposed by the exploration hole from the exploration hole columnar diagram of the geological exploration data, so that the coordinates of each exploration hole and the surface elevation correspond one by one with the elevation of each stratigraphic interface they expose, forming the point coordinates at the contact part between the exploration hole and each stratigraphic interface.

[0010] 1.2) Convert the point coordinates at the contact part between the exploration hole and each stratigraphic interface into three-dimensional discrete points. For the discrete points of each stratum in turn, interpolate the discrete points to form the topographic surface of the stratum surface.

[0011] 1.3) Extract the curved surface boundaries of each stratum based on the topographic surface of the stratum surface and stretch them to generate a solid. Cut the solid with each stratum surface as the boundary plane to form a three-dimensional geological model of the stratum.

[0012] Optionally, the step of interpolating the discrete points to form the topographic surface of the stratum surface in step 1.2) includes: first, combine the discrete points to generate an irregular triangular mesh, and then form a smooth surface from the triangular mesh to form the topographic surface of the stratum surface.

[0013] Optionally, the step of calculating the projection distance from the position coordinate to the stratum surface in the three-dimensional geological model in step 2) includes: first, project the position coordinate vertically onto each stratum surface in the three-dimensional geological model to obtain the projection points; then extract the elevations of each projection point and calculate the vertical distance between the projection points on each stratum surface, so as to obtain the thickness of the overburden layer of each stratum.

[0014] Optionally, after step 2), it further includes the steps of determining the distance from the pile foundation opening position to the top surface of the bearing stratum, and generating a pile foundation model using the surface point at the pile foundation opening position, the projection point on the top surface of the bearing stratum, and the pile diameter of the pile foundation opening.

[0015] Optionally, when generating the pile foundation model, it further includes the steps of exporting the three-dimensional coordinates of the surface point and the thickness of the surface point from the bearing stratum.

[0016] Optionally, establishing a three-dimensional geological model of the stratum based on the geological exploration data in step 1), and calculating the projection distance from the position coordinate to the stratum surface in the three-dimensional geological model as the thickness of the overburden layer of the position coordinate in step 2) are both implemented based on BIM software.

[0017] In addition, this embodiment also provides a calculation system for the thickness of the overburden layer of any stratum, including a microprocessor and a memory connected to each other. The microprocessor is programmed or configured to execute the steps of the calculation method for the thickness of the overburden layer of any stratum.

[0018] In addition, this embodiment also provides a calculation system for the thickness of the covering layer of any formation, including a microprocessor and a memory connected to each other. A computer program programmed or configured to execute the calculation method for the thickness of the covering layer of any formation is stored in the memory.

[0019] In addition, this embodiment also provides a computer-readable storage medium, in which a computer program programmed or configured to execute the calculation method for the thickness of the covering layer of any formation is stored.

[0020] Compared with the prior art, the present invention has the following advantages: The present invention includes establishing a three-dimensional geological model of the formation based on geological exploration data; importing the position coordinates where the thickness of the formation covering layer needs to be calculated into the three-dimensional geological model, and calculating the projection distance from the position coordinates to the formation surface in the three-dimensional geological model as the thickness of the formation covering layer at the position coordinates. The present invention can quickly and accurately establish a three-dimensional geological model based on geological exploration data, with high accuracy and speed. The whole process from importing data to generating the model takes about 10 - 30 minutes; the present invention can add newly added borehole exploration points in real time during the construction process to achieve dynamic modification of the three-dimensional geological model. Using the present invention, it is possible to extract the thickness of the covering layer of any formation below any point within the boundary of the geological exploration data. It has the advantages of fast calculation speed, high accuracy, intuitive graphics, accurate numerical values, and can avoid human calculation errors, which has great significance for the refined construction of projects such as rock foundation construction and pile foundation construction. BRIEF DESCRIPTION OF THE DRAWINGS

[0021] Figure 1 It is a schematic diagram of the basic process of the method in the embodiment of the present invention.

[0022] Figure 2 It is a node diagram for reading an excel file in the embodiment of the present invention.

[0023] Figure 3 It is a node diagram for generating point graphics based on Cartesian coordinates in an excel file in the embodiment of the present invention.

[0024] Figure 4 It is a node diagram for generating terrain and surface based on discrete points in the embodiment of the present invention.

[0025] Figure 5 It is to repeat Figure 4 nodes to establish a node diagram of a multi-layer formation surface model.

[0026] Figure 6 It is a node diagram for calculating the distance from the center point of the borehole opening position along the borehole direction to any formation surface in the embodiment of the present invention.

[0027] Figure 7 It is a node diagram for generating a borehole model according to the borehole position and diameter in the embodiment of the present invention.

[0028] Figure 8 It is the node diagram for creating a model parameter table for a single model in the embodiments of the present invention.

[0029] Figure 9 It is the node diagram for importing the model created in the embodiments of the present invention into BIM software.

[0030] Figure 10 It is the node diagram for recording information such as the calculated distance in the embodiments of the present invention through the model. Detailed implementation manners

[0031] As Figure 1 shown, the calculation method for the thickness of any formation overburden layer in this embodiment includes:

[0032] 1) Establish a three-dimensional geological model of the formation based on geological exploration data;

[0033] 2) Import the position coordinates where the formation overburden layer thickness needs to be calculated into the three-dimensional geological model, and calculate the projection distance from the position coordinates to the formation surface in the three-dimensional geological model as the formation overburden layer thickness of the position coordinates.

[0034] Refer to Figure 1 , as an optional implementation manner, before step 1) in this embodiment, there is also a step of sorting out and importing geological exploration data. The sorted geological exploration data in this embodiment includes:

[0035] The plan layout diagram of geological exploration data, including the coordinates of each exploration hole and the surface elevation information;

[0036] The boring log of geological exploration data, containing the elevation information of each formation interface exposed by the exploration hole.

[0037] In this embodiment, step 1) includes:

[0038] 1.1) Extract the coordinates of each exploration hole and the surface elevation according to the plan layout diagram of geological exploration data. For each exploration hole, respectively extract the elevation of each formation interface exposed by the exploration hole according to the boring log of geological exploration data, so that the coordinates of each exploration hole and the surface elevation correspond one by one with the elevation of each formation interface they expose, and form the point coordinates of the contact part between the exploration hole and each formation interface;

[0039] 1.2) Convert the point coordinates of the contact part between the exploration hole and each formation interface into three-dimensional discrete points. For the discrete points of each formation in turn, interpolate the discrete points to form the terrain surface of the formation surface;

[0040] 1.3) Based on the topographic surface of the stratum surface, the surface boundaries of each stratum are extracted and stretched to generate entities, and the entities are cut with the surfaces of each stratum as interfaces to form a three-dimensional geological model of the stratum.

[0041] In this embodiment, the step of interpolating discrete points to form a topographic surface of the stratum surface in step 1.2) includes: firstly, combining the discrete points to generate an irregular triangular mesh, and then forming a smooth surface from the triangular mesh to form a topographic surface of the stratum surface.

[0042] In this embodiment, the step of calculating the projection distance of the position coordinate to the stratum surface in the three-dimensional geological model in step 2) includes: first, projecting the position coordinate along the vertical direction onto each stratum surface in the three-dimensional geological model to obtain projection points; then extracting the elevation of each projection point, calculating the vertical distance between the projection points on each stratum surface, and thus obtaining the covering layer thickness of each stratum.

[0043] As an optional implementation, for the application of estimating the drilling depth of a pile foundation project, this embodiment also includes, after step 2), the step of determining the distance between the pile foundation opening position and the top surface of the bearing layer, and generating a pile foundation model using the surface point of the pile foundation opening position, the projection point on the top surface of the bearing layer, and the pile diameter of the pile foundation opening.

[0044] As an optional implementation, generating the pile foundation model in this embodiment further includes the step of deriving the three-dimensional coordinates of the surface point and the thickness of the surface point from the bearing layer.

[0045] In this embodiment, the three-dimensional geological model of the stratum is established according to the geological survey data in step 1), and the projection distance from the position coordinate to the stratum surface in the three-dimensional geological model is calculated as the stratum cover thickness of the position coordinate in step 2), which are both implemented based on BIM software. The following will take BIM software as an example to further explain in detail the calculation method of the thickness of any stratum cover in this embodiment, and its operation steps based on BIM software include:

[0046] S1. Import the geological survey data in Excel format into BIM software through the three nodes of BIM software file path (File Path), extract file (File.FromPath), and read Excel content (Excel.ReadFromFile) so that it can be converted into a three-dimensional geological model later. Figure 2 shown.

[0047] S2. Process the imported geological survey data in Excel format, remove redundant information such as stratum names, and retain coordinate data.

[0048] S3. Extract the coordinates of each exploration hole and the surface elevation from the layout plan of the geological exploration data. For each exploration hole, extract the elevation of each stratigraphic interface exposed by the exploration hole from the exploration hole columnar diagram of the geological exploration data, so that the coordinates of each exploration hole and the surface elevation correspond one by one with the elevation of each stratigraphic interface it exposes, forming the point coordinates at the contact part between the exploration hole and each stratigraphic interface. As Figure 3 shown, which is the content of step 1.3) above;

[0049] S4. Use the Point.ByCoordinates node of the BIM software to convert the point coordinates at the contact part between the exploration hole and each stratigraphic interface into three-dimensional discrete points; for the discrete points of each stratum in turn, first use the Topography.ByPoints node of the BIM software to combine the discrete points to generate an irregular triangular mesh, and then use the Topography.ToPolySurface-SpringNodes node of the BIM software to form a smooth surface (Surface surface) from the triangular mesh. As Figure 4 shown, thus forming the topographic surface of the stratum surface. If the calculation of the stratum thickness involves multiple strata, step S4 can be repeated to complete the modeling of the topographic surface of each stratum surface, and its nodes are as Figure 5 shown.

[0050] S5. After generating the three-dimensional geological model of the stratum, the position coordinates where the stratum overburden thickness needs to be calculated can be imported into the three-dimensional geological model. If the calculation of the overburden thickness involves multiple points, the three nodes of File Path, File.FromPath, and Excel.ReadFromFile in step S1 can be reused to batch import the position coordinates where the stratum overburden thickness needs to be calculated;

[0051] S6. After importing the position coordinates where the stratum overburden thickness needs to be calculated, first use the Point.Project node of the BIM software to project the position coordinates vertically (0, 0, -1) onto the surface (Surface surface) of each stratum in the three-dimensional geological model to obtain the projection points;

[0052] S7. Extract the elevation of each projection point, calculate the vertical distance between the projection points on the surface of each stratum, so as to obtain the overburden thickness of each stratum. The nodes are as Figure 6 shown.

[0053] After step 2) of this embodiment, it also includes the steps of determining the distance between the pile foundation opening position and the top surface of the bearing stratum, and generating a pile foundation model using the surface point of the pile foundation opening position, the projection point on the top surface of the bearing stratum, and the pile diameter of the pile foundation opening. In order to

[0054] As an alternative embodiment, when generating the pile foundation model in this embodiment, it further includes the steps of exporting the three-dimensional coordinates of the ground surface points and the thickness of the ground surface points from the bearing stratum. The overburden thickness calculation involved in this embodiment is a batch calculation at multiple positions. To ensure that the calculation results correspond one by one to the calculation positions, the calculation results can be recorded through the model for convenient later access to the results. According to requirements, use the Parameter.CreatProjectParameter node of the BIM software to create parameters for the model. The parameters created in this embodiment include the pile number, the xyz coordinates of the ground surface, and the thickness of the ground surface from the bearing stratum. The nodes are as Figure 8 shown.

[0055] In addition, this embodiment further includes importing the model generated in Dynamo in the BIM software into Revit in the BIM software for easy result display and access. The nodes are as Figure 9 shown.

[0056] In addition, this embodiment further includes writing the calculation results of the overburden thickness of the batch calculation into the model one by one in sequence. The nodes are as Figure 10 shown. It is also possible to directly export the calculation results of step S7 to Excel by selecting the WriteToExcel node of the BIM software.

[0057] It should be noted that the calculation method of the overburden thickness of any stratum in this embodiment does not depend on the BIM software, because establishing a three-dimensional geological model of the stratum based on the geological exploration data in step 1) and calculating the projection distance from the position coordinates to the stratum surface in the three-dimensional geological model as the overburden thickness of the position coordinates in step 2) can also be implemented using other three-dimensional software.

[0058] In summary, at present, domestic stratigraphic profiles usually adopt structural maps, and the source of their drawing data mainly depends on actual drilling data. The stratigraphic cover layer depth calculated by this method has high accuracy when the borehole spacing is small. However, when the borehole spacing is large, in order to ensure the accuracy of the data based on the design, the exploration boreholes are often encrypted, such as redrawing the stratigraphic profile, which is time-consuming and laborious. Some similar three-dimensional geological software can effectively solve the problem of repeatedly drawing stratigraphic profiles, but its output is limited to the profile on the line between the exploration points and the bar chart of the points on the line. The profile drawn by this means can only reflect the stratum fluctuations between the drilling holes and the drilling holes, and it is not convenient to calculate the stratigraphic conditions outside the exploration hole line. The calculation method of the thickness of any stratigraphic cover layer in this embodiment can quickly and accurately establish a three-dimensional geological model based on geological survey data, with high precision and high speed. Newly added drilling exploration points can be added in real time during the construction process, and the real-time model is dynamically modified. The stratigraphic thickness of any unknown part is calculated by interpolation analysis using known exploration points.

[0059] In addition, this embodiment also provides a system for calculating the thickness of an overburden layer of any stratum, comprising a microprocessor and a memory connected to each other, wherein the microprocessor is programmed or configured to execute the steps of the aforementioned method for calculating the thickness of an overburden layer of any stratum.

[0060] In addition, this embodiment also provides a system for calculating the thickness of an arbitrary stratum covering layer, comprising a microprocessor and a memory connected to each other, wherein the memory stores a computer program programmed or configured to execute the aforementioned method for calculating the thickness of an arbitrary stratum covering layer.

[0061] In addition, this embodiment also provides a computer-readable storage medium, which stores a computer program programmed or configured to execute the aforementioned method for calculating the thickness of any stratum covering layer.

[0062] Those skilled in the art should understand that the embodiments of the present application can be provided as methods, systems, or computer program products. Therefore, the present application can take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present application can take the form of a computer program product implemented on one or more computer-readable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code. The present application is a flowchart of the method, device (system), and computer program product according to the embodiment of the present application, and / or the processor executes instructions to generate instructions for implementing the flowchart. Figure 1 A process or multiple processes and / or boxes Figure 1a device for the functions specified in one or more boxes. These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, such that the instructions stored in the computer-readable memory produce a manufactured article including an instruction device that implements in the process Figure 1 one process or multiple processes and / or boxes Figure 1 a device for the functions specified in one or more boxes. These computer program instructions can also be loaded onto a computer or other programmable data processing device, such that a series of operating steps are performed on the computer or other programmable device to produce a computer-implemented process, so that the instructions executed on the computer or other programmable device provide for implementing in the process Figure 1 one process or multiple processes and / or boxes Figure 1 a device for the functions specified in one or more boxes.

[0063] The above is only the preferred embodiment of the present invention, and the protection scope of the present invention is not limited to the above embodiments. All technical solutions falling within the concept of the present invention belong to the protection scope of the present invention. It should be noted that for those of ordinary skill in the art in this technical field, several improvements and refinements made without departing from the principle of the present invention should also be regarded as within the protection scope of the present invention.

Claims

1. A calculation method for the thickness of any formation covering layer, characterized in that, Including: 1) Establish a three-dimensional geological model of the strata based on geological exploration data; 2) Import the position coordinates where the thickness of the stratum overburden needs to be calculated into the three-dimensional geological model, and calculate the projection distance from the position coordinates to the stratum surface in the three-dimensional geological model as the thickness of the stratum overburden at the position coordinates; The step of calculating the projection distance from the position coordinates to the stratum surface in the three-dimensional geological model includes: First, project the position coordinates vertically onto each stratum surface in the three-dimensional geological model to obtain projection points; then extract the elevations of each projection point, and calculate the vertical distances between the projection points on each stratum surface, so as to obtain the thickness of the overburden of each stratum; Step 1) includes: 1.1) Extract the coordinates of each exploration hole and the surface elevation from the plane layout diagram of the geological exploration data. For each exploration hole, extract the elevation of each stratum interface exposed by the exploration hole from the exploration hole columnar diagram of the geological exploration data, so that the coordinates of each exploration hole and the surface elevation correspond one by one to the elevation of each stratum interface exposed by it, forming the point coordinates of the contact part between the exploration hole and each stratum interface; 1.2) Convert the point coordinates of the contact part between the exploration hole and each stratum interface into three-dimensional discrete points. For the discrete points of each stratum in turn, interpolate the discrete points to form the topographic surface of the stratum surface; 1.3) Extract the curved surface boundaries of each stratum based on the topographic surface of the stratum surface and stretch them into entities, and cut the entity with each stratum surface as the boundary surface to form a three-dimensional geological model of the strata; The step of interpolating the discrete points to form the topographic surface of the stratum surface in step 1.2) includes: First, combine the discrete points to generate an irregular triangular mesh, and then form a smooth surface from the triangular mesh, so as to form the topographic surface of the stratum surface; After step 2), it also includes the step of determining the distance from the pile foundation opening position to the top surface of the bearing stratum, and generating a pile foundation model by using the surface point at the pile foundation opening position, the projection point on the top surface of the bearing stratum and the pile diameter of the pile foundation opening.

2. The calculation method of the thickness of any formation covering layer according to claim 1, characterized in that, When generating the pile foundation model, it also includes the step of exporting the three-dimensional coordinates of the surface point and the thickness of the surface point from the bearing stratum.

3. The calculation method of the thickness of any formation covering layer according to claim 1, characterized in that, Both establishing a three-dimensional geological model of the strata based on geological exploration data in step 1) and calculating the projection distance from the position coordinates to the stratum surface in the three-dimensional geological model as the thickness of the stratum overburden at the position coordinates in step 2) are implemented based on BIM software.

4. A calculation system for the thickness of any formation overburden, comprising a microprocessor and a memory connected to each other, characterized in that, The microprocessor is programmed or configured to execute the steps of the method for calculating the thickness of any stratum overburden described in any one of claims 1 to 3.

5. A calculation system for the thickness of any formation covering layer, comprising a microprocessor and a memory connected to each other, characterized in that, The memory stores a computer program that is programmed or configured to execute the method for calculating the thickness of any stratum overburden described in any one of claims 1 to 3 through the microprocessor.

6. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that is programmed or configured to execute the method for calculating the thickness of any stratum overburden described in any one of claims 1 to 3 through the microprocessor.

Images