Method and related device for complex formation cave volume calculation

Abstract

The embodiment of the application provides a complex stratum cave volume calculation method and related equipment, and belongs to the technical field of geological survey. The method comprises the following steps: obtaining a geological profile and a drilling columnar graph; performing information extraction on the drilling columnar graph to obtain drilling information of each drilling point; performing structured processing on the drilling information of all drilling point positions to obtain an original stratum distribution table; performing elevation supplement on the original stratum distribution table according to the geological profile to obtain a standard stratum distribution table; based on the standard stratum distribution table, performing data splitting on each stratum as a derived unit to obtain a plurality of stratum derived tables; performing three-dimensional model reconstruction on each stratum derived table to obtain a stratum interface entity model of each stratum; performing target detection on all stratum interface entity models according to a preset cave label to obtain a target cave model; and performing face domain attribute analysis on the target cave model to obtain a target cave volume. The application can improve the accuracy of cave volume calculation.

Description

Technical Field

[0001] This application relates to the field of geological exploration technology, and in particular to a method and related equipment for calculating the volume of karst caves in complex strata. Background Technology

[0002] Currently, the mainstream methods for calculating the volume of karst caves include empirical estimation, two-dimensional drawing calculation, and equipment measurement.

[0003] Empirical estimation methods are susceptible to human error and often contain significant errors. These errors are particularly pronounced in complex geological formations and may lead to underestimation or overestimation of the volume of karst caves, thereby affecting the accuracy of engineering design and the safety of construction.

[0004] While two-dimensional CAD drawings can provide some geological information, their limitations become increasingly apparent when dealing with complex strata. Furthermore, two-dimensional drawings struggle to visually represent the spatial morphology and distribution of caves, failing to meet the demands of modern engineering construction for three-dimensional visualization and precise modeling.

[0005] The method for measuring the volume of karst caves based on the ideal gas law has certain effectiveness in engineering applications in urban karst sites, but it has disadvantages such as high equipment cost, high equipment risk, low measurement efficiency, and great influence from the site environment and temperature, making it unsuitable for calculating the volume of karst cave entities on construction sites.

[0006] Furthermore, the above methods are difficult to achieve efficient data sharing and updating, and also present many inconveniences in collaborative work, making it difficult to adapt to the trend of the construction industry's transformation towards intelligent design and management. Summary of the Invention

[0007] The main objective of this application is to propose a method and related equipment for calculating the volume of karst caves in complex strata, aiming to achieve efficient conversion from geological data to three-dimensional geological models and ensure the accuracy of karst cave volume calculation.

[0008] To achieve the above objectives, one aspect of this application proposes a method for calculating the volume of karst caves in complex strata, the method comprising:

[0009] Acquire multi-source geological survey data, wherein the multi-source geological survey data includes geological profile maps and borehole columnar sections of multiple borehole locations;

[0010] Information is extracted from the borehole columnar section to obtain borehole information for each borehole point, wherein the borehole information includes coordinate information, number information, and elevation information;

[0011] The borehole information from all borehole locations is structured to obtain the original stratigraphic distribution table;

[0012] Based on the geological profile, the missing strata in the original stratigraphic distribution table are supplemented by elevation to obtain a standard stratigraphic distribution table.

[0013] Based on the standard stratigraphic distribution table, data is split into several stratigraphic derivation tables by taking each stratigraphic unit as the derivation unit.

[0014] A three-dimensional model is reconstructed based on each of the aforementioned stratigraphic derivation tables to obtain a stratigraphic interface entity model for each stratigraphic layer.

[0015] Based on the preset cave tags, target detection is performed on all the aforementioned strata interface entity models to obtain the target cave model;

[0016] The surface attribute analysis of the target cave model is performed to obtain the volume of the target cave.

[0017] In some embodiments, the coordinate information includes abscissa and ordinate, the numbering information includes borehole number and formation number, and the elevation information includes borehole elevation and bottom elevation of all formations. The step of performing structured processing based on the borehole information of all borehole locations to obtain the original formation distribution table includes the following steps:

[0018] Using the borehole number as the first column, the horizontal axis as the second column, and the vertical axis as the third column, the borehole number and its corresponding horizontal and vertical axes are entered sequentially according to the left-right arrangement order of the borehole number in the borehole bar chart to obtain the original borehole table. The vertical arrangement order of the borehole number in the original borehole table is consistent with the left-right arrangement order of the borehole number in the borehole bar chart.

[0019] Based on the vertical arrangement of the strata numbers in the borehole columnar section and the horizontal arrangement of the bottom elevations, the elevation information is entered sequentially line by line into the original borehole table to obtain the original strata distribution table.

[0020] In some embodiments, the step of supplementing the missing strata in the original stratigraphic distribution table with elevation based on the geological profile to obtain a standard stratigraphic distribution table includes the following steps:

[0021] The original stratigraphic distribution table is used to query the stratigraphic layers with null bottom elevation values ​​to obtain the missing stratigraphic layers and the associated stratigraphic layers adjacent to the missing stratigraphic layers. The associated stratigraphic layers include local overlying strata and local underlying strata.

[0022] The missing type is determined based on the local overlying strata and the local underlying strata to obtain the missing result, wherein the missing type includes single-layer missing and multi-layer missing;

[0023] When the missing result is a single-layer missing, the local overlying stratum is determined as the reference stratum;

[0024] When the missing result is multi-layer missing, the stratigraphic sequence between the associated strata is reconstructed according to the geological profile, the global overlying strata are determined based on the reconstructed stratigraphic sequence, and the global overlying strata are determined as the reference strata.

[0025] The missing strata are sequentially inherited and supplemented based on the bottom elevation of the reference strata to obtain a standard stratum distribution table.

[0026] In some embodiments, determining the missing type based on the local overlying strata and the local underlying strata to obtain the missing result includes the following steps:

[0027] Query the number of strata located between the local overlying strata and the local underlying strata in the original stratigraphic distribution table;

[0028] The missing type is determined based on the number of strata. If the number of strata is equal to one, the missing result is confirmed to be a single-layer missing.

[0029] If the number of strata is greater than one, the missing result is confirmed to be a multi-layer missing result.

[0030] In some embodiments, the step of sequentially inheriting and supplementing the missing strata according to the bottom elevation of the reference strata to obtain a standard stratigraphic distribution table includes the following steps:

[0031] The bottom elevation of the reference stratum is copied as the bottom elevation of the missing stratum to obtain supplementary elevation data;

[0032] The supplementary elevation data is labeled, and the labeled supplementary elevation data is filled into the original stratigraphic distribution table to obtain the standard stratigraphic distribution table.

[0033] In some embodiments, the step of splitting data based on the standard stratigraphic distribution table, with each stratigraphic unit as the derivation unit, to obtain several stratigraphic derivation tables includes the following steps:

[0034] The list is initialized using each column of the standard stratigraphic distribution table as a derivation unit, resulting in several original derivation tables;

[0035] Based on the standard stratigraphic distribution table, the original derived table is mapped according to the derived unit to obtain a stratigraphic derived table. Each row of data in the stratigraphic derived table corresponds to a borehole location, and each row of data includes the horizontal coordinate, the vertical coordinate, and the bottom elevation of the derived unit.

[0036] In some embodiments, the step of performing target detection on all the geological interface entity models according to preset cave tags to obtain the target cave model includes the following steps:

[0037] Target matching is performed based on the preset cave label and the stratum number to obtain a target number that matches the preset cave label;

[0038] The target cave model is obtained by isolating the geological interface entity model corresponding to the target number from all the geological interface entity models.

[0039] To achieve the above objectives, another aspect of this application proposes a system for calculating the volume of karst caves in complex strata, the system comprising:

[0040] The data acquisition module is used to acquire multi-source geological survey data, which includes geological profile maps and borehole columnar sections of multiple borehole locations.

[0041] The information extraction module is used to extract information from the borehole columnar section to obtain borehole information for each borehole point, wherein the borehole information includes coordinate information, number information and elevation information;

[0042] The structuring module is used to perform structuring processing on the borehole information of all borehole locations to obtain the original stratigraphic distribution table;

[0043] The data supplementation module is used to supplement the elevation of the missing strata in the original stratigraphic distribution table based on the geological profile map, so as to obtain a standard stratigraphic distribution table.

[0044] The data splitting module is used to split the data based on the standard stratigraphic distribution table, with each stratigraphic unit as the derivation unit, to obtain several stratigraphic derivation tables.

[0045] The model reconstruction module is used to reconstruct a three-dimensional model based on each of the stratigraphic derivation tables to obtain the stratigraphic interface entity model of each stratigraphic layer.

[0046] The target detection module is used to perform target detection on all the geological interface entity models according to the preset cave labels to obtain the target cave model;

[0047] The volume calculation module is used to perform surface attribute analysis on the target cave model to obtain the volume of the target cave.

[0048] To achieve the above objectives, another aspect of this application provides an electronic device, which includes a memory and a processor. The memory stores a computer program, and the processor executes the computer program to implement the above-described method.

[0049] To achieve the above objectives, another aspect of the embodiments of this application proposes a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method.

[0050] To achieve the above objectives, another aspect of this application provides a computer program product, including a computer program that, when executed by a processor, implements the above-described method.

[0051] The embodiments of this application include at least the following beneficial effects: This application provides a method and related equipment for calculating the volume of karst caves in complex strata. This scheme acquires multi-source geological survey data, including geological profile maps and borehole columnar sections of multiple borehole locations; extracts information from the borehole columnar sections to obtain borehole information for each borehole location, including coordinate information, number information, and elevation information; performs structured processing based on the borehole information of all borehole locations to obtain an original stratigraphic distribution table; supplements the missing strata in the original stratigraphic distribution table with elevations based on the geological profile maps to obtain a standard stratigraphic distribution table; based on the standard stratigraphic distribution table, splits the data with each stratum as a derivative unit to obtain several stratigraphic derivative tables; and reconstructs a three-dimensional model based on each stratigraphic derivative table to obtain a stratigraphic interface entity model for each stratum, thus achieving efficient conversion from geological data to a three-dimensional geological model. Based on preset cave tags, target detection is performed on all geological interface entity models to obtain target cave models; surface attribute analysis is performed on the target cave models to obtain the target cave volume, which can ensure the accuracy of cave volume calculation. Attached Figure Description

[0052] Figure 1 This is a flowchart of the method for calculating the volume of karst caves in complex strata provided in the embodiments of this application;

[0053] Figure 2 This is a schematic diagram of the detailed survey drawings provided in the embodiments of this application;

[0054] Figure 3 This is a schematic diagram of the borehole columnar section provided in the embodiments of this application;

[0055] Figure 4 This is a schematic diagram of the original stratigraphic distribution table provided in the embodiments of this application;

[0056] Figure 5 This is a schematic diagram of the standard stratigraphic distribution table provided in the embodiments of this application;

[0057] Figure 6 This is a schematic diagram of multiple stratigraphic derivation tables provided in the embodiments of this application;

[0058] Figure 7This is a schematic diagram of the target karst cave volume calculation results provided in the embodiments of this application;

[0059] Figure 8 This is a schematic diagram of the 0-2D volume measurement of a karst cave provided in an embodiment of this application;

[0060] Figure 9 This is a schematic diagram of the original borehole table provided in an embodiment of this application;

[0061] Figure 10 This is a schematic diagram of the original stratigraphic distribution table provided in another embodiment of this application;

[0062] Figure 11 This is a schematic diagram showing the missing original stratigraphic distribution table data provided in the embodiments of this application;

[0063] Figure 12 This is a schematic diagram illustrating the sequential inheritance supplement provided in the embodiments of this application;

[0064] Figure 13 This is a schematic diagram of a geological profile provided in an embodiment of this application;

[0065] Figure 14 This is a schematic diagram showing the missing original stratigraphic distribution table data provided in another embodiment of this application;

[0066] Figure 15 This is a schematic diagram of the stratigraphic derivation table provided in the embodiments of this application;

[0067] Figure 16 This is a schematic diagram of the stratigraphic color information provided in the embodiments of this application;

[0068] Figure 17 This is a schematic diagram of surface generation provided in an embodiment of this application;

[0069] Figure 18 This is a schematic diagram of the imported stratigraphic derivation table provided in an embodiment of this application;

[0070] Figure 19 This is a schematic diagram of the generated orifice surface layer provided in an embodiment of this application;

[0071] Figure 20 This is a schematic diagram of the generation of 1-plain fill 1-2 surfaces provided in the embodiments of this application;

[0072] Figure 21 This is a schematic diagram illustrating the generation of solids from a surface according to an embodiment of this application;

[0073] Figure 22 This is a schematic diagram of the 1-plain fill soil 1-2 stratum entity interface model provided in the embodiments of this application;

[0074] Figure 23This is a schematic diagram of the isolated target karst cave model provided in the embodiments of this application. Detailed Implementation

[0075] To make the objectives, technical solutions, and advantages of this application clearer, the following detailed description is provided in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of this application and are not intended to limit it. In the following description, when referring to the accompanying drawings, unless otherwise indicated, the same numbers in different drawings represent the same or similar elements. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with those of this application; they are merely examples of systems and methods consistent with some aspects of the embodiments of this application as detailed in the appended claims.

[0076] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs. The terminology used herein is for the purpose of describing embodiments of this application only and is not intended to limit this application.

[0077] The method for calculating the volume of karst caves in complex strata provided in this application relates to the field of geological exploration. This method can be applied to a terminal, a server, or software running on either a terminal or a server. In some embodiments, the terminal can be a smartphone, tablet, laptop, desktop computer, smart speaker, smartwatch, or vehicle-mounted terminal, but is not limited to these. The server can be configured as an independent physical server, a server cluster or distributed system composed of multiple physical servers, or a cloud server providing basic cloud computing services such as cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communication, middleware services, domain name services, security services, CDN, and big data and artificial intelligence platforms. The server can also be a node server in a blockchain network. The software can be an application that implements the method for calculating the volume of karst caves in complex strata, but is not limited to the above forms.

[0078] This application can be used in a wide variety of general-purpose or special-purpose computer system environments or configurations. Examples include: personal computers, server computers, handheld or portable devices, tablet devices, multiprocessor systems, microprocessor-based systems, set-top boxes, programmable consumer electronics devices, network PCs, minicomputers, mainframe computers, and distributed computing environments including any of the above systems or devices. This application can be described in the general context of computer-executable instructions executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, etc., that perform specific tasks or implement specific abstract data types. This application can also be practiced in distributed computing environments where tasks are performed by remote processing devices connected via a communication network. In distributed computing environments, program modules can reside in local and remote computer storage media, including storage devices.

[0079] Figure 1 This is an optional flowchart of the method for calculating the volume of karst caves in complex strata provided in the embodiments of this application. Figure 1 The method may include, but is not limited to, steps S101 to S108.

[0080] Step S101: Obtain multi-source geological survey data, which includes geological profile maps and borehole columnar sections of multiple borehole locations.

[0081] Step S102: Extract information from the borehole columnar section to obtain borehole information for each borehole location. The borehole information includes coordinate information, number information, and elevation information.

[0082] Step S103: Perform structured processing based on the borehole information of all borehole locations to obtain the original stratigraphic distribution table.

[0083] Step S104: Based on the geological profile, the missing strata in the original stratigraphic distribution table are supplemented with elevation data to obtain a standard stratigraphic distribution table.

[0084] Step S105: Based on the standard stratigraphic distribution table, the data is split into several stratigraphic derivation tables, with each stratigraphic unit as the derivation unit.

[0085] Step S106: Reconstruct the three-dimensional model based on the stratum derivation table to obtain the stratum interface entity model of each stratum.

[0086] Step S107: Target detection is performed on all geological interface entity models according to the preset cave labels to obtain the target cave model.

[0087] Step S108: Perform surface attribute analysis on the target cave model to obtain the target cave volume.

[0088] In this embodiment, the geological data consists of CAD drawings provided by the surveying unit, mainly including borehole data such as advanced drilling, detailed exploration, and edge probing. Advanced drilling is an exploration of the pile location conducted before pile foundation construction, primarily using drilling methods to obtain rock core samples from the bearing stratum at the pile tip. Analysis of the rock cores determines the rock properties and integrity. Detailed exploration, or detailed investigation, is a more in-depth geological investigation based on the preliminary exploration. It involves the comprehensive application of various methods to fully understand the site's geological conditions, such as... Figure 2 As shown, geological data from explorations are typically stored in the form of borehole plan locations, geological columnar sections, and geological profiles.

[0089] Specifically, multi-source geological survey data mainly uses geological profiles from CAD drawings that reflect the continuous spatial variation of strata and borehole columnar sections from various discrete borehole locations. The model is built by extracting elevation points from the borehole columnar sections, and later, the geological profiles are combined for auxiliary modeling to achieve efficient conversion from geological data to three-dimensional geological models.

[0090] Understandably, each borehole column chart independently records the vertical stratigraphic conditions at the location of a borehole point. By obtaining borehole column charts from multiple discrete borehole points across the entire survey area, a geological profile can be drawn that shows the actual (or inferred) cross-section of the geological structure at the surface or within a certain depth, providing constraints on the spatial relationships of the strata.

[0091] Optionally, to ensure the integrity and accuracy of the three-dimensional geological model, all borehole information provided by the exploration unit can be used.

[0092] After acquiring multi-source geological survey data, artificial intelligence technology can be used to extract coordinate information for determining planar location, numbering information for identification and tracing, and elevation information for constructing three-dimensional space from borehole columnar diagrams through image recognition. The coordinate information includes the borehole's x-coordinate (X) and y-coordinate (Y) coordinates, the numbering information includes the borehole number and stratum number, and the elevation information includes the borehole opening elevation and the bottom elevation of all strata.

[0093] For example, such as Figure 3As shown, the CAD drawings are exported as PDFs and uploaded to a pre-trained image detection neural network model. The model automatically extracts information such as "borehole number", "borehole elevation", "coordinates", "stratum number", and "bottom elevation" from the borehole columnar section. For example, the borehole with "borehole number" "DK02" has "coordinates" of "X=233836.43, Y=35198.85" and "borehole diameter" of "7.62". The bottom elevation of the stratum with "stratum number" "plain fill 1-2" is "5.62" and the bottom elevation of the stratum with "stratum number" "silty fine sand 2-2" is "3.22", and so on.

[0094] It should be noted that the types of information extracted and the format of the output information by the image detection neural network model can be adjusted according to actual needs. This application is only exemplary and does not impose any specific limitations.

[0095] The borehole information from all borehole locations is structured and unified to construct a two-dimensional matrix. The borehole number is used as the row index, and the stratigraphic number of each borehole location is used as the column index. The borehole information extracted in step S102 is then entered into the corresponding positions in the two-dimensional matrix, resulting in an original stratigraphic distribution table, as shown below. Figure 4 As shown, the borehole information corresponding to borehole number DK02 will be sequentially filled into the original stratigraphic distribution table, for reference. Figure 3 It can be seen that borehole DK02 has 9 layers, from top to bottom: borehole opening, plain fill 1-2, silty fine sand 2-2, ..., slightly weathered carbonaceous limestone 9C-1. Therefore, according to the borehole number DK02 and the corresponding stratum number, the elevation information is filled into the corresponding position.

[0096] Because the underground strata themselves have faults and the boreholes are sparsely distributed during the survey, not every borehole can collect data from all strata. Therefore, as borehole information from different boreholes is entered, the number of columns in the two-dimensional matrix will increase with the newly appearing strata until all borehole columnar sections are traversed, all collected borehole information is recorded, and an original stratum distribution table is obtained.

[0097] Furthermore, in the original stratigraphic distribution table, due to the differences in strata at different borehole locations, some borehole locations will have elevation information for that stratum, while others will lack such information. This results in a large number of strata with blank bottom elevations in the original stratigraphic distribution table. To address this issue, this embodiment employs data interpolation. For strata with unknown data due to missing information, elevation information is supplemented based on known elevation data. Different data supplementation methods are used depending on the nature of the missing information. For example, when the missing stratum in a borehole is located between two consecutive... When there are strata with available data, the sequence of the missing strata is clear, and the supplementary elevation information can directly inherit the sequence of the strata above the missing strata. However, if there are multiple different strata between the upper and lower strata of the missing strata, there are multiple possible insertion positions for the missing strata. In this case, spatial constraints need to be applied based on the geological profile, and a comprehensive judgment needs to be made using stratum development patterns, pinch-out point locations, and the sequence of surrounding boreholes to determine the actual insertion position of the missing strata at the current borehole location. Then, data supplementation is performed based on the determined sequence, ultimately completing all the missing data in the original stratigraphic distribution table to obtain a complete table. Figure 5 The complete standard stratigraphic distribution table is shown.

[0098] Next, using each column of the standard stratigraphic distribution table as a derivation unit, a separate planar coordinate table containing only the plane coordinates of all borehole locations and the bottom elevation of the stratigraphic layer at each borehole location is generated for each stratigraphic layer, resulting in a stratigraphic derivation table for each layer. By splitting the data of each column of stratigraphic layers, several stratigraphic derivation tables with the same number of layers are obtained, such as... Figure 6 As shown.

[0099] After obtaining the stratigraphic derivation table, the discrete two-dimensional data stored in each stratigraphic derivation table are transformed into a three-dimensional model that can measure volume. The three-dimensional model reconstruction process is mainly divided into two stages: surface generation and solid generation. First, using the stratigraphic derivation table as the input data for reconstruction, a corresponding surface is created for each stratigraphic layer using the idea of ​​points forming surfaces. By transforming multiple discrete elevation data recorded in the stratigraphic derivation table into continuous surfaces, stratigraphic interfaces representing the spatial interface morphology at the bottom of the stratigraphic layer are obtained. Then, based on the idea of ​​merging two stratigraphic interfaces into a single entity, two adjacent stratigraphic interfaces are closed to generate a solid located between the surfaces, thus obtaining the stratigraphic interface solid model of each stratigraphic layer. Multiple stratigraphic interface solid models can be tightly stacked to reconstruct the spatial distribution of the subsurface layers.

[0100] Finally, based on the preset cave tags, the research subject, namely the cave, is separated from all the stratigraphic interface entity models. It should be noted that in this embodiment, each cave has been defined as an independent stratum. Based on the preset cave tags, semantic query and filtering are performed. By comparing whether the stratum number of each stratigraphic interface entity model is consistent with the preset cave tag, the target cave model is separated from all the stratigraphic interface entity models. Then, the volume of the closed target cave model is calculated by using the surface attribute analysis function in the 3D modeling software. The volume of each cave is output simultaneously, as well as the total target cave volume after summing the volumes of multiple caves.

[0101] For example, refer to Figure 7 Using 26 pre-set cave tags, target detection is performed in 95 geological interface entity models. Each time a geological interface model matching a cave tag is detected, the geological interface model of that cave is isolated.

[0102] like Figure 8 As shown, the volume of a single isolated cave (0-2D) is obtained. The volume is calculated using the massprop function, yielding a volume of 52864.826 m³. 3 By summing the volumes of all the caves, the target cave volume is obtained as 1620969.528 m³. 3 .

[0103] In some embodiments, step S103 may include, but is not limited to, steps S201 to S202.

[0104] Step S201: Using the borehole number as the first column, the horizontal axis as the second column, and the vertical axis as the third column, according to the left-right arrangement order of the borehole numbers in the borehole bar chart, enter the borehole numbers and their corresponding horizontal and vertical axes in sequence to obtain the original borehole table. The vertical arrangement order of the borehole numbers in the original borehole table is consistent with the left-right arrangement order of the borehole numbers in the borehole bar chart.

[0105] Step S202: Based on the vertical arrangement of the strata numbers in the borehole columnar section and the order of the bottom elevations, enter the elevation information line by line in the original borehole table to obtain the original strata distribution table.

[0106] In this embodiment, refer to Figure 9The original stratigraphic distribution table is an Excel spreadsheet created based on borehole information extracted from borehole columnar sections. The first column (column A) corresponds to the borehole number in the drawing, and the second column (column B) and third column (column C) represent the X and Y coordinates of the borehole locations, respectively. To ensure the systematic nature of the data entry, the borehole numbering order from left to right in the multiple borehole columnar sections in the drawing is used as the baseline order for data entry. The borehole numbers are filled into the first column sequentially from top to bottom. Then, the x and y coordinates corresponding to each borehole number are entered into the second and third columns of the same row, respectively, resulting in a sequence reflecting the planar distribution of boreholes. This establishes a basic data framework with spatial location indexed by borehole numbers for subsequent steps, resulting in the original borehole table.

[0107] For example, assuming the borehole column chart in the drawing is numbered DK02, DK04, DK05, DK06... from left to right, then in this order, enter DK02 in the first row of the borehole numbers, DK04 in the second row, and so on, until all borehole numbers have been entered. When entering the borehole numbers, simultaneously enter the X and Y coordinates in the last two columns of the row containing the borehole number.

[0108] It should be noted that the borehole numbering order from top to bottom in the original borehole table must be consistent with the borehole numbering order from left to right in the borehole column diagram on the drawing.

[0109] After the original borehole table was constructed, the elevation information of each stratum was gradually filled in to form the original stratigraphic distribution table.

[0110] Specifically, using the borehole number as an index, the stratigraphic number and description of each stratum are read from top to bottom, starting from the borehole opening. For example, according to... Figure 3 The borehole columnar section shown is arranged in the following order: borehole opening, plain fill 1-2, silty fine sand 2-2, ..., slightly weathered carbonaceous limestone 9C-1. This order is used as the stratigraphic sequence of the bottom layer. Starting from 0, each stratum is assigned a stratigraphic sequence number. For example, the stratigraphic sequence number of the borehole opening is 0, and the stratigraphic sequence number of the plain fill is 1. Starting from the fourth column in the header of the original borehole table, the elevation of each stratum at the bottom of the borehole number is entered row by row in ascending order of stratigraphic sequence number.

[0111] For example, such as Figure 4 As shown, using borehole number DK02 as an index, the borehole elevation of 7.62 is entered in the fourth column of the first row, and the bottom elevation of fill soil 1-2 of 5.62 is entered in the fifth column, and so on. Because the stratigraphic sequence is defined in a top-down order, the elevation information in the same row decreases from left to right.

[0112] It should be noted that the strata detected at different borehole locations vary. Therefore, there may be instances where some strata are present at one borehole location but not at another. Figure 10 As shown, at this point, there will be strata that are not present in DK02 but are present in DK04, such as "moderately weathered carbonaceous limestone 8C-1", or strata that are present in DK02 but not in DK04, such as "slightly weathered carbonaceous limestone 9C-1". In this case, it is necessary to supplement the missing strata in the list header determined by DK02. However, since these two strata do not appear in the same borehole location, it is not possible to directly determine their stratigraphic sequence through their vertical relationship. Therefore, we will not explore the stratigraphic sequence relationship between the two at this time, but instead insert the missing strata between the adjacent strata based on their respective adjacent strata at the borehole location.

[0113] Understandably, during insertion, it's still necessary to ensure that the bottom elevations of all strata decrease from left to right in each row after insertion. For example... Figure 11 As shown, borehole DK01 recorded four strata, borehole DK02 recorded six strata, and borehole DK03 recorded six strata. Based on the decreasing order of strata 1-2, 1-1, and 2-1 recorded in DK02, it can be determined that strata 1-2 with a bottom elevation of 6.87 recorded in DK02 should be inserted between columns E and G. However, based on strata 1-2 with a bottom elevation of 4.70 and strata 2-1 with a bottom elevation of 4.91 recorded in DK03, it can be seen that the bottom elevation of strata 1-2 recorded in DK02 is lower than that of strata 2-1. Therefore, the stratigraphic sequence of this stratum is different from that of strata 1-2 recorded in DK02. Thus, a new column for strata 1-2 should be added after column G to ensure the decreasing order of the data.

[0114] By preserving the blanks of missing data, the uncertainty of stratigraphy caused by data gaps is reflected, providing real data support for subsequent standardization and data supplementation of stratigraphic sequences. The decreasing nature of the data reflects the faulting phenomenon of the strata.

[0115] In some embodiments, step S104 may include, but is not limited to, steps S301 to S305.

[0116] Step S301: Query the original stratigraphic distribution table for stratigraphic layers with null bottom elevation values ​​to obtain the missing stratigraphic layers and the associated stratigraphic layers adjacent to the missing stratigraphic layers. The associated stratigraphic layers include local overlying strata and local underlying strata.

[0117] Step S302: Determine the missing type based on the local overlying strata and the local underlying strata to obtain the missing result. The missing type includes single-layer missing and multi-layer missing.

[0118] Step S303: When the missing result is a single-layer missing, the local overlying strata are determined as the reference strata.

[0119] Step S304: When the missing result is multi-layer missing, the stratigraphic sequence between the associated strata is reconstructed according to the geological profile. The global overlying strata are determined based on the reconstructed stratigraphic sequence and the global overlying strata are determined as the reference strata.

[0120] Step S305: Based on the bottom elevation of the reference strata, the missing strata are sequentially inherited and supplemented to obtain the standard stratum distribution table.

[0121] In this embodiment, the original stratigraphic distribution table is scanned to identify all cells with blank data. Each blank cell represents a specific missing stratigraphic layer, meaning that the stratigraphic layer corresponding to the blank cell was not recorded at the corresponding borehole location. The stratigraphic layer that is closest to the missing stratigraphic layer and has bottom elevation data records is searched in the data of the obtained borehole location to be identified as the associated stratigraphic layer. The local overlying stratigraphic layer is the stratigraphic layer that covers the missing stratigraphic layer and is represented in the original stratigraphic distribution table as the first stratigraphic layer with data to the left of the missing stratigraphic layer. The local underlying stratigraphic layer is the stratigraphic layer that is deposited below the missing stratigraphic layer and is represented in the original stratigraphic distribution table as the first stratigraphic layer with data to the right of the missing stratigraphic layer.

[0122] The type of missing strata is determined by counting the total number of strata between the local overlying strata and the local underlying strata. A single missing strata indicates that there is only one missing stratum between the local overlying strata and the local underlying strata, while a multi-layer missing strata indicates that there are other strata between the local overlying strata and the local underlying strata besides the missing stratum.

[0123] When the missing result is a single layer missing, the sequence of the missing stratum is uniquely determined. Its sequence is located after the local overlying stratum and before the local underlying stratum. Therefore, the local overlying stratum can be directly used as the reference stratum for subsequent data supplementation.

[0124] When the missing result is multi-layered, a single borehole data cannot resolve the possibility of multiple solutions. The missing strata may be adjacent to the local overlying strata, or there may be several missing strata separated from the local overlying strata. Therefore, it is necessary to combine the geological profile map to see the actual insertion location of the missing strata.

[0125] For example, Figure 13 The geological profile provided in this embodiment is referenced. Figure 14Taking the inserted slightly weathered 9C-1 stratum as an example, in the borehole columnar section of the drawing, 9C-1 is located between the green-marked 2-2 (local overlying stratum) and 0-1A (local underlying stratum). However, since there are already strata 3-2 / 4N-2 / 3-1 in strata 2-2 and 0-1A, there are many possibilities for inserting 9C-1, such as between 2-2 and 3-2, between 3-2 and 4N-2, between 4N-2 and 3-1, or between 3-1 and 0-1A. Different insertion positions will affect the supplementary elevation when supplementing data. Therefore, it is necessary to first determine the actual global overlying stratum of 9C-1 as 3-1 by combining the geological profile map, and then use the global overlying stratum as the reference stratum for the missing strata in the case of multiple missing layers.

[0126] After determining the reference stratigraphic layers for each missing stratum, a sequential inheritance supplementation process is performed uniformly. The bottom elevation of the reference layer is used to fill in the blank cells in the missing strata, ultimately resulting in a complete and ordered standard stratigraphic distribution table without null values, such as... Figure 12 As shown, assuming the reference stratum for the missing stratum 8C-1 is 5C-18, the bottom elevation of 5C-18 is copied as the bottom elevation of 8C-1.

[0127] In some embodiments, step S302 may include, but is not limited to, steps S401 to S403.

[0128] Step S401: Query the number of strata located between the local overlying strata and the local underlying strata in the original stratigraphic distribution table.

[0129] Step S402: Determine the missing type based on the number of strata. If the number of strata is equal to one, then the missing result is confirmed as a single-layer missing.

[0130] Step S403: If the number of layers is greater than one, the missing result is confirmed to be multi-layered missing.

[0131] In this embodiment, after identifying the missing strata, the missing scale caused by the missing strata is quantified to determine the missing type.

[0132] Specifically, for the missing strata in the current analysis, by querying the number of strata from the local overlying strata to the local underlying strata, if the number of strata is exactly one, it indicates that the missing strata is uniquely sandwiched between the local overlying strata and the local underlying strata. For example... Figure 11 In column F, 1-2 has only layer 1-2 between its local overlying stratum 1-1 and its local underlying stratum 2-1. Therefore, its stratigraphic sequence can be determined to be one layer after 1-1.

[0133] If the number of strata is greater than one, it indicates that there may be a gap between the overlying and underlying strata, consisting of multiple strata with unclear internal sequences. In this case, due to the lack of data, it is impossible to accurately determine the spatial order of these strata. Therefore, the current result of the missing data is confirmed as multi-layered missing data. Multi-layered missing data signifies that data supplementation requires initiating a more advanced processing procedure, such as calling geological profiles and using spatial constraints to reconstruct the sequence stratigraphy.

[0134] In some embodiments, step S305 may include, but is not limited to, steps S501 to S502.

[0135] Step S501: Copy the bottom elevation of the reference stratum as the bottom elevation of the missing stratum to obtain supplementary elevation data.

[0136] Step S502: Label the supplementary elevation data and fill the labeled supplementary elevation data into the original stratigraphic distribution table to obtain the standard stratigraphic distribution table.

[0137] In this embodiment, sequential inheritance supplementation is an operation that uses a reference stratum above the missing stratum to perform amplitude calculations. This operation is based on a clear geological logic. When the stratigraphic sequence is clear and there are no multiple solutions, the missing stratum can be regarded as having the same bottom elevation as the overlying stratum that is continuous in the sequence. By copying the bottom elevation of the reference stratum as supplementary elevation data, it can be ensured that the missing stratum is smoothly connected in subsequent geological 3D modeling, thereby providing a geometric basis for volume calculation.

[0138] Furthermore, to distinguish which data in the standard stratigraphic distribution table were obtained through field surveys and which data were obtained through data supplementation in this application, data annotation can be used to annotate the supplementary elevation data.

[0139] Preferably, the supplementary elevation data can be highlighted in red for visual annotation. The annotated supplementary elevation data is then filled into the original stratigraphic distribution table to fill in the missing strata, ultimately forming a standard stratigraphic distribution table where all cells have values ​​and the data source is clearly defined, such as... Figure 5 As shown.

[0140] In some embodiments, step S105 may include, but is not limited to, steps S601 to S602.

[0141] Step S601: Initialize the list by using each column of the standard stratigraphic distribution table as the derivation unit to obtain several original derivation tables.

[0142] Step S602: Based on the standard stratigraphic distribution table, the original derived table is mapped according to the derived unit to obtain the stratigraphic derived table. Each row of data in the stratigraphic derived table corresponds to a borehole location, and each row of data includes the horizontal coordinate, the vertical coordinate, and the bottom elevation of the derived unit.

[0143] In this embodiment, after obtaining the standard stratigraphic distribution table, in order to be suitable for 3D modeling software, the standard stratigraphic distribution table is deconstructed into a series of original derived tables.

[0144] Specifically, the first step is to initialize the list to create a data container, using the stratigraphy of each column in the standard stratigraphic distribution table as the derived unit, for example... Figure 5 The stratum with the soil layer number (i.e., the sequence number) of plain fill soil 1-2 is a derived unit, and the slightly weathered carbonaceous limestone 0-2C with the soil layer number 87 is a derived unit. A corresponding original derivation table is generated for each of the 96 strata with soil layer numbers from 0 to 95.

[0145] Next, using the derived unit as an index, a scan and data extraction are performed in the standard stratigraphic distribution table. The x-coordinate, y-coordinate, and bottom elevation of the corresponding stratum of the derived unit below the borehole point in each row of the standard stratigraphic distribution table are written into the original derived table, such as... Figure 15 As shown, by converting the abscissa and ordinate of the borehole locations, as well as the bottom elevation, into three-dimensional information of the strata, and aggregating the discrete points from different borehole locations, a three-dimensional point set containing all survey points of the strata can be obtained, providing core parameters for subsequent point-to-surface modeling methods.

[0146] In some embodiments, step S107 may include, but is not limited to, steps S701 to S702.

[0147] Step S701: Target matching is performed based on the preset cave label and stratum number to obtain the target number that matches the preset cave label.

[0148] Step S702: Isolate the stratum interface entity model corresponding to the target number from all stratum interface entity models to obtain the target karst cave model.

[0149] In this embodiment, the preset cave label is a specific code for identifying the geology of the cave, such as 0-2D. By semantically matching the preset cave label and the stratum number, it is compared one by one to see if there is a stratum number that matches the preset cave label. If there is a stratum number that matches the preset cave label, then the stratum number is determined as the target number.

[0150] It is understood that the form and number of preset cave labels can be set according to actual needs. This application is only exemplary and does not impose any specific restrictions.

[0151] After obtaining the target number, the associated stratigraphic interface entity model is isolated from all other stratigraphic interface entity models based on the target number. This isolation typically involves hiding or turning off all stratigraphic interface entity models that are not associated with the target number, retaining and highlighting only the stratigraphic interface entity models selected by the target number. By isolating these entities, the karst cave entity, which is part of the geology, is independently extracted, forming a target karst cave model containing only the cave and unaffected by other strata. The volume of this target karst cave model is then calculated, yielding the karst cave volume.

[0152] The following is a detailed description and explanation of the solutions in the embodiments of the present invention, using specific application examples:

[0153] With the acceleration of urbanization, underground space development has become an important way to expand urban space and enhance urban functions. However, underground space development faces complex and diverse geological environments, among which the problem of karst caves in complex strata is particularly prominent. The existence of karst caves brings many challenges to engineering construction, and accurate calculation of karst cave volume is crucial for engineering design, construction safety, and the assessment of grouting volume for karst cave treatment. Precise data on karst cave volume can help engineers assess geological risks, optimize design schemes, and provide a basis for project cost accounting.

[0154] In view of this, this application utilizes Civil 3D's ability to process complex stratigraphic data and proposes a method for calculating the volume of karst caves in complex stratigraphic formations based on Civil 3D three-dimensional geological modeling. This method achieves efficient conversion from geological data to a three-dimensional geological model by extracting CAD drawing data provided by the exploration unit, creating detailed data tables, generating the model, and performing volume calculations. This ensures the accuracy of the karst cave volume calculation and provides a new technical approach and practical reference for calculating the volume of karst caves in similar complex stratigraphic engineering projects.

[0155] Specifically, accurate borehole data is crucial for building reliable geological models, and Civil 3D can make full use of this data to generate accurate three-dimensional geological models, providing a model carrier for calculating the volume of karst caves in complex soil layers.

[0156] First, identify the drawings and data provided by the surveying unit, extract borehole information based on the borehole columnar section and geological profile of all strata and karst cave borehole locations, and then create the original stratigraphic distribution table in an Excel spreadsheet.

[0157] Optionally, when encountering a new stratum during the data entry process, the corresponding stratum column number should be added to the table in a timely manner. Alternatively, a borehole location with a relatively complete stratum can be selected before data entry, and the stratum number and information can be pre-entered into the table.

[0158] Next, considering the possible data omissions or anomalies in the drawings, reasonable data supplementation or correction measures were taken in a timely manner to obtain a standard stratigraphic distribution table, thereby providing reliable data support for subsequent geological model construction and cave volume calculation.

[0159] It should be noted that during the data supplementation process, if the insertion location is difficult to determine, for example... Figure 14 The slightly weathered carbonaceous limestone 9C-1 shown has multiple possible insertion locations. In such cases, it is necessary to refer to the geological profile to determine the accurate insertion location of the strata. The auxiliary analysis of the geological profile can ensure the accuracy of the stratigraphic data and the rationality of the geological model.

[0160] The purpose of collecting tabular data in the above steps is to establish a unified stratigraphic sequence and comprehensively sort the strata, serving as the data basis for the generation of subsequent geological models and the calculation of cave volumes. Without establishing a unified stratigraphic sequence, it is impossible to generate a usable three-dimensional geological model from complex strata.

[0161] Since Civil 3D software only supports comma-separated values ​​(CSV) format data files, each stratum needs to be stored separately as a CSV file, using plain text to form a CSV table, where each row represents a record and fields are separated by commas.

[0162] Specifically, based on the standard stratigraphic distribution table, data is split into derived units for each stratigraphic layer. Each table needs to enter coordinate data in three columns: X, Y, and Z. The X and Y coordinates are the same, while the Z coordinate (elevation) is different. The tables are named according to the rule of "sequence number-stratigraphic description-stratigraphic number", such as "1-plain fill 1-2".

[0163] Furthermore, in Civil 3D software, constructing reasonable and effective layer information plays a crucial role in ensuring the clear presentation and orderly arrangement of the geological model, preparing for subsequent calculation of cave volume. The layer information should correspond to the information of each stratum.

[0164] For example, you can refer to Figure 5 All integrated stratigraphic information was used to determine the soil layer categories within the project in a newly created table using paste, transpose, and filter commands, ultimately resulting in 26 soil categories. Each soil category was assigned a corresponding color to distinguish the properties of different strata. Civil 3D software was then opened, and corresponding layer information was created in the layer properties. Based on this layer information, a soil layer and cave information table was created for subsequent volumetric data entry and statistics. Figure 16 As shown.

[0165] Specifically, layer generation can be referenced. Figure 17 First, in Civil 3D software, right-click on the surface and create a surface. Rename it to "0-Orifice High Level," the same name as the first CSV table. Then, as follows... Figure 18 As shown, in the newly created surface, select "Point File," right-click to add, and click the "+" icon in the upper right corner of the interface to import the corresponding stratum's CSV table. Click "OK" to generate the first wellhead surface level. Figure 19 As shown.

[0166] Next, repeat the process of creating surfaces, point files, and generating layers, importing the CSV table of 1-fill soil 1-2 to generate as shown. Figure 20 The second side shown, then select 1-plain fill 1-2 side, right-click and select, in Figure 21 In the interface shown, select "Generate Solid from Surface", choose the appropriate layers and colors, and then generate the terrain solid, as shown. Figure 22 The solid model of the 1-plain fill soil 1-2 stratum interface is shown.

[0167] Repeat the above steps until all 95 strata corresponding to the strata interface entity models are generated. In the layer properties, select all strata and soil layers, right-click and isolate the selected layers. This will display all 95 entities as shown below. Figure 23 As shown.

[0168] Taking the 0-2D cave layer as an example, by isolating the 0-2D cave layer and using the massprop function to calculate the soil volume, the result can be obtained. Figure 8 The volume calculation results shown are used to fill the volume of the cave into a large table. By isolating and calculating the volume of each cave layer, the following results are obtained: Figure 7 The large table shown is used to summarize the volumes recorded in the table, thus completing the volume calculation of the target cave.

[0169] This application also provides a system for calculating the volume of karst caves in complex strata, which can implement the above-mentioned method. The system includes:

[0170] The data acquisition module is used to acquire multi-source geological survey data, which includes geological profile maps and borehole columnar sections from multiple borehole locations.

[0171] The information extraction module is used to extract information from the borehole columnar section to obtain borehole information for each borehole point. The borehole information includes coordinate information, number information, and elevation information.

[0172] The structuring module is used to perform structuring processing on the borehole information of all borehole locations to obtain the original stratigraphic distribution table.

[0173] The data supplementation module is used to supplement the elevation of missing strata in the original stratigraphic distribution table based on the geological profile map, so as to obtain a standard stratigraphic distribution table.

[0174] The data splitting module is used to split data based on the standard stratigraphic distribution table, with each stratigraphic unit as the derivation unit, to obtain several stratigraphic derivation tables.

[0175] The model reconstruction module is used to reconstruct the three-dimensional model based on each stratigraphic derivation table, and obtain the stratigraphic interface entity model of each stratigraphic layer.

[0176] The target detection module is used to perform target detection on all geological interface entity models based on preset cave labels to obtain the target cave model.

[0177] The volume calculation module is used to perform surface attribute analysis on the target cave model to obtain the volume of the target cave.

[0178] It is understood that the content of the above method embodiments is applicable to this system embodiment. The specific functions implemented in this system embodiment are the same as those in the above method embodiments, and the beneficial effects achieved are also the same as those achieved in the above method embodiments.

[0179] This application also provides an electronic device, which includes a memory and a processor. The memory stores a computer program, and the processor executes the computer program to implement the above-described method. This electronic device can be any smart terminal, including tablet computers, in-vehicle computers, etc.

[0180] It is understood that the content of the above method embodiments is applicable to this device embodiment. The specific functions implemented by this device embodiment are the same as those of the above method embodiments, and the beneficial effects achieved are also the same as those achieved by the above method embodiments.

[0181] This application also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method.

[0182] It is understood that the content of the above method embodiments is applicable to this storage medium embodiment. The specific functions implemented in this storage medium embodiment are the same as those in the above method embodiments, and the beneficial effects achieved are also the same as those achieved in the above method embodiments.

[0183] This application also provides a computer program product, including a computer program that, when executed by a processor, implements the above-described method.

[0184] It is understood that the content of the above method embodiments is applicable to the embodiments of this program product. The specific functions implemented by the embodiments of this program product are the same as those of the above method embodiments, and the beneficial effects achieved are also the same as those achieved by the above method embodiments.

[0185] Memory, as a non-transitory computer-readable storage medium, can be used to store non-transitory software programs and non-transitory computer-executable programs. Furthermore, memory may include high-speed random access memory, and may also include non-transitory memory, such as at least one disk storage device, flash memory device, or other non-transitory solid-state storage device. In some embodiments, memory may optionally include memory remotely located relative to the processor, and these remote memories can be connected to the processor via a network. Examples of such networks include, but are not limited to, the Internet, intranets, local area networks, mobile communication networks, and combinations thereof.

[0186] The method and related equipment for calculating the volume of karst caves in complex strata provided in this application acquire multi-source geological survey data, including geological profiles and borehole columnar sections of multiple borehole locations. Information is extracted from the borehole columnar sections to obtain borehole information for each borehole location, including coordinate, number, and elevation information. The borehole information from all locations is then structured to obtain an original stratigraphic distribution table. The missing strata in the original stratigraphic distribution table are supplemented with elevation data based on the geological profiles to obtain a standard stratigraphic distribution table. Based on the standard stratigraphic distribution table, data is split using each stratum as a derivative unit to obtain several stratigraphic derivative tables. A three-dimensional model is reconstructed based on each stratigraphic derivative table to obtain a stratigraphic interface entity model for each stratum. This enables efficient conversion from geological data to a three-dimensional geological model. Based on preset cave tags, target detection is performed on all geological interface entity models to obtain target cave models; surface attribute analysis is performed on the target cave models to obtain the target cave volume, which can ensure the accuracy of cave volume calculation.

[0187] The embodiments described in this application are for the purpose of more clearly illustrating the technical solutions of the embodiments of this application, and do not constitute a limitation on the technical solutions provided by the embodiments of this application. As those skilled in the art will know, with the evolution of technology and the emergence of new application scenarios, the technical solutions provided by the embodiments of this application are also applicable to similar technical problems.

[0188] Those skilled in the art will understand that the technical solutions shown in the figures do not constitute a limitation on the embodiments of this application, and may include more or fewer steps than shown, or combine certain steps, or different steps.

[0189] The system embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate; that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs.

[0190] Those skilled in the art will understand that all or some of the steps in the methods disclosed above, as well as the functional modules / units in the systems and devices, can be implemented as software, firmware, hardware, or suitable combinations thereof.

[0191] The terms “first,” “second,” “third,” “fourth,” etc. (if present) in the specification and accompanying drawings of this application are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that the embodiments of this application described herein can be implemented in orders other than those illustrated or described herein. Furthermore, the terms “comprising” and “having,” and any variations thereof, are intended to cover non-exclusive inclusion; for example, a process, method, system, product, or apparatus that comprises a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or apparatus.

[0192] It should be understood that in this application, "at least one (item)" means one or more, and "more than" means two or more. "And / or" is used to describe the relationship between related objects, indicating that three relationships can exist. For example, "A and / or B" can represent three cases: only A exists, only B exists, and both A and B exist simultaneously, where A and B can be singular or plural. The character " / " generally indicates that the preceding and following related objects are in an "or" relationship. "At least one (item) of the following" or similar expressions refer to any combination of these items, including any combination of single or plural items. For example, at least one (item) of a, b, or c can represent: a, b, c, "a and b", "a and c", "b and c", or "a and b and c", where a, b, and c can be single or multiple.

[0193] In the embodiments provided in this application, it should be understood that the disclosed systems and methods can be implemented in other ways. For example, the system embodiments described above are merely illustrative; for instance, the division of the units described above is only a logical functional division, and in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the coupling or direct coupling or communication connection shown or discussed may be through some interfaces; the indirect coupling or communication connection between systems or units may be electrical, mechanical, or other forms.

[0194] The units described above as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of this embodiment according to actual needs.

[0195] Furthermore, the functional units in the various embodiments of this application can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. The integrated unit can be implemented in hardware or as a software functional unit.

[0196] If the integrated unit is implemented as a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, or all or part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes multiple instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods of the various embodiments of this application. The aforementioned storage medium includes various media capable of storing programs, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0197] The preferred embodiments of the present application have been described above with reference to the accompanying drawings, but this does not limit the scope of the claims of the present application. Any modifications, equivalent substitutions, and improvements made by those skilled in the art without departing from the scope and substance of the embodiments of the present application shall be within the scope of the claims of the present application.

Claims

1. A method for calculating the volume of karst caves in complex strata, characterized in that, The method includes the following steps: Acquire multi-source geological survey data, wherein the multi-source geological survey data includes geological profile maps and borehole columnar sections of multiple borehole locations; Information is extracted from the borehole columnar section to obtain borehole information for each borehole point, wherein the borehole information includes coordinate information, number information, and elevation information; The borehole information from all borehole locations is structured to obtain the original stratigraphic distribution table; Based on the geological profile, the missing strata in the original stratigraphic distribution table are supplemented by elevation to obtain a standard stratigraphic distribution table. Based on the standard stratigraphic distribution table, data is split into several stratigraphic derivation tables by taking each stratigraphic unit as the derivation unit. A three-dimensional model is reconstructed based on each of the aforementioned stratigraphic derivation tables to obtain a stratigraphic interface entity model for each stratigraphic layer. Based on the preset cave tags, target detection is performed on all the aforementioned strata interface entity models to obtain the target cave model; The surface attribute analysis of the target cave model is performed to obtain the volume of the target cave.

2. The method according to claim 1, characterized in that, The coordinate information includes abscissa and ordinate; the numbering information includes borehole number and stratum number; the elevation information includes borehole elevation and bottom elevation of all strata; the step of performing structured processing based on the borehole information of all borehole locations to obtain the original stratum distribution table includes the following steps: Using the borehole number as the first column, the horizontal axis as the second column, and the vertical axis as the third column, the borehole number and its corresponding horizontal and vertical axes are entered sequentially according to the left-right arrangement order of the borehole number in the borehole bar chart to obtain the original borehole table. The vertical arrangement order of the borehole number in the original borehole table is consistent with the left-right arrangement order of the borehole number in the borehole bar chart. Based on the vertical arrangement of the strata numbers in the borehole columnar section and the horizontal arrangement of the bottom elevations, the elevation information is entered sequentially line by line into the original borehole table to obtain the original strata distribution table.

3. The method according to claim 2, characterized in that, The process of supplementing the missing strata in the original stratigraphic distribution table with elevation data based on the geological profile to obtain a standard stratigraphic distribution table includes the following steps: The original stratigraphic distribution table is used to query the stratigraphic layers with null bottom elevation values ​​to obtain the missing stratigraphic layers and the associated stratigraphic layers adjacent to the missing stratigraphic layers. The associated stratigraphic layers include local overlying strata and local underlying strata. The missing type is determined based on the local overlying strata and the local underlying strata to obtain the missing result, wherein the missing type includes single-layer missing and multi-layer missing; When the missing result is a single-layer missing, the local overlying stratum is determined as the reference stratum; When the missing result is multi-layer missing, the stratigraphic sequence between the associated strata is reconstructed according to the geological profile, the global overlying strata are determined based on the reconstructed stratigraphic sequence, and the global overlying strata are determined as the reference strata. The missing strata are sequentially inherited and supplemented based on the bottom elevation of the reference strata to obtain a standard stratum distribution table.

4. The method according to claim 3, characterized in that, The process of determining the missing type based on the local overlying strata and the local underlying strata to obtain the missing result includes the following steps: Query the number of strata located between the local overlying strata and the local underlying strata in the original stratigraphic distribution table; The missing type is determined based on the number of strata. If the number of strata is equal to one, the missing result is confirmed to be a single-layer missing. If the number of strata is greater than one, the missing result is confirmed to be a multi-layer missing result.

5. The method according to claim 3, characterized in that, The step of sequentially inheriting and supplementing the missing strata according to the bottom elevation of the reference strata to obtain a standard stratigraphic distribution table includes the following steps: The bottom elevation of the reference stratum is copied as the bottom elevation of the missing stratum to obtain supplementary elevation data; The supplementary elevation data is labeled, and the labeled supplementary elevation data is filled into the original stratigraphic distribution table to obtain the standard stratigraphic distribution table.

6. The method according to claim 2, characterized in that, The process of splitting data based on the standard stratigraphic distribution table, with each stratigraphic unit as the derivation unit, to obtain several stratigraphic derivation tables includes the following steps: The list is initialized using each column of the standard stratigraphic distribution table as a derivation unit, resulting in several original derivation tables; Based on the standard stratigraphic distribution table, the original derived table is mapped according to the derived unit to obtain a stratigraphic derived table. Each row of data in the stratigraphic derived table corresponds to a borehole location, and each row of data includes the horizontal coordinate, the vertical coordinate, and the bottom elevation of the derived unit.

7. The method according to claim 1, characterized in that, The step of performing target detection on all the geological interface entity models based on preset cave tags to obtain the target cave model includes the following steps: Target matching is performed based on the preset cave label and the stratum number to obtain a target number that matches the preset cave label; The target cave model is obtained by isolating the geological interface entity model corresponding to the target number from all the geological interface entity models.

8. A system for calculating the volume of karst caves in complex strata, characterized in that, The system includes: The data acquisition module is used to acquire multi-source geological survey data, which includes geological profile maps and borehole columnar sections of multiple borehole locations. The information extraction module is used to extract information from the borehole columnar section to obtain borehole information for each borehole point, wherein the borehole information includes coordinate information, number information and elevation information; The structuring module is used to perform structuring processing on the borehole information of all borehole locations to obtain the original stratigraphic distribution table; The data supplementation module is used to supplement the elevation of the missing strata in the original stratigraphic distribution table based on the geological profile map, so as to obtain a standard stratigraphic distribution table. The data splitting module is used to split the data based on the standard stratigraphic distribution table, with each stratigraphic unit as the derivation unit, to obtain several stratigraphic derivation tables. The model reconstruction module is used to reconstruct a three-dimensional model based on each of the stratigraphic derivation tables to obtain the stratigraphic interface entity model of each stratigraphic layer. The target detection module is used to perform target detection on all the geological interface entity models according to the preset cave labels to obtain the target cave model; The volume calculation module is used to perform surface attribute analysis on the target cave model to obtain the volume of the target cave.

9. An electronic device, characterized in that, The electronic device includes a memory and a processor, the memory storing a computer program, and the processor executing the computer program to implement the method according to any one of claims 1 to 7.

10. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by a processor, it implements the method of any one of claims 1 to 7.

Images