CN114998367A

CN114998367A - Plane partitioning method and device for karst development degree and electronic equipment

Info

Publication number: CN114998367A
Application number: CN202210628769.2A
Authority: CN
Inventors: 彭丁茂; 桂炎德; 施兹国; 陈侃福; 舒宗运; 徐国锋; 崔永兴
Original assignee: Zhejiang Shuzhijiaoyuan Technology Co Ltd
Current assignee: Zhejiang Shuzhijiaoyuan Technology Co Ltd
Priority date: 2022-06-06
Filing date: 2022-06-06
Publication date: 2022-09-02
Anticipated expiration: 2042-06-06
Also published as: CN114998367B

Abstract

The application provides a plane partitioning method for karst development degree, which comprises the following steps: determining a drilling number accumulated value according to the distribution characteristics of the karst section limestone accumulated height and the characteristics of the limestone top surface elevation slope change in the drilling karst single index relation graph, and determining a first vertical line and a second vertical line according to the drilling number accumulated value and the drilling number of the set ratio in the karst drilling index scatter diagram; determining a first datum plane and a second datum plane according to the first perpendicular line, the second perpendicular line and an index value in the karst drilling index scatter diagram, wherein the first datum plane and the second datum plane are partitioned elevation lines of the karst of a target area on a geological vertical section; the target plane is partitioned by the first reference plane and the second reference plane. According to the method and the device, the karst plane is partitioned by determining the partition elevation line of the karst of the target area on the geological longitudinal section through the elevation line so as to determine the karst development distribution condition on the plane.

Description

Plane partitioning method and device for karst development degree and electronic equipment

Technical Field

The application relates to the field of engineering geology, in particular to a plane partitioning method and device for karst development degree and electronic equipment.

Background

Karst is a general and complex unfavorable geological phenomenon, has strong heterogeneity, and can cause geological disasters under certain conditions, thereby seriously threatening engineering safety. Therefore, before the construction of the target area, the karst distribution of the target area needs to be determined in advance according to the development degree of the karst, and then engineering design is carried out. Generally, the development degree of the karst needs to be determined according to the geological zoning condition of the karst. However, the karst is difficult to partition geologically by the current exploration means and analysis method, and the development degree of the karst cannot be intuitively understood.

Disclosure of Invention

In view of the above, an object of the embodiments of the present application is to provide a method, an apparatus, an electronic device and a readable storage medium for plane partitioning of a karst development degree. The karst geology can be partitioned, and then the karst development condition on the plane can be determined.

In a first aspect, an embodiment of the present application provides a planar partitioning method for a karst development degree, including: determining a drilling number accumulated value according to distribution characteristics of karst section limestone accumulated height and elevation slope change characteristics of limestone top surface in a drilling karst single index relational graph, wherein the drilling karst single index relational graph is a relational graph in which a plurality of single index values of drilling holes are projected in a first target coordinate system according to a set sequence; determining a first perpendicular line and a second perpendicular line according to the accumulated value of the number of the drill holes and the number of the drill holes with a set ratio in a karst drill hole index scatter diagram, wherein the karst drill hole index scatter diagram is a relational diagram which takes the elevation of a karst cave bottom plate as a reference according to a set sequence and projects all index values of the drill holes corresponding to the elevation of the karst cave bottom plate in a second target coordinate system, and the first perpendicular line and the second perpendicular line are used for dividing the karst drill hole index scatter diagram into a plurality of areas according to the distribution situation of all index data of the drill holes; determining a first datum plane and a second datum plane according to the first perpendicular line, the second perpendicular line and an index numerical value in a karst drilling index scatter diagram, wherein the first datum plane and the second datum plane are partition elevation lines for partitioning a target area karst on a geological vertical section; and partitioning the target plane through the first reference plane and the second reference plane.

In the implementation process, the acquired index values of the drill holes are processed in the coordinate system, so that different types of index relation graphs are obtained. Corresponding processing values can be obtained based on different index relation graphs, and the shallow erosion reference surface of the valley region and the shallow erosion reference surface of the mountain front slope region can be determined through the processing values, so that the karst plane can be partitioned, and the karst development condition on the plane can be determined.

With reference to the first aspect, an embodiment of the present application provides a first possible implementation manner of the first aspect, where: the determining of the first vertical line and the second vertical line according to the accumulated value of the number of drill holes and the number of drill holes with the set ratio in the karst drill hole index scatter diagram comprises: determining a first vertical line according to the accumulated value of the drilling number; determining the number of drill holes between the first set value of the karst cave bottom plate elevation and the first vertical line, wherein the drill holes meet the set ratio; and determining a second perpendicular line according to the number of the drill holes meeting the set ratio.

In the implementation process, the first perpendicular line and the second perpendicular line are determined in the karst drilling index scatter diagram, each index of a drill hole in the karst drilling index scatter diagram can be partitioned, the actual distribution situation of each index of the drill hole is convenient to analyze, and the accuracy of plane partitioning of the karst development degree is improved.

With reference to the first possible implementation manner of the first aspect, an embodiment of the present application provides a second possible implementation manner of the first aspect, where: the index values include: cumulative height of karst cave, karst cave bottom plate elevation and drill way elevation, according to index value in first perpendicular, second perpendicular and the karst drilling index scatter diagram confirms first reference surface and second reference surface, includes: determining a first datum plane according to the first vertical line, the accumulated height of the karst cave and the elevation of the bottom plate of the karst cave; a second datum is determined from the aperture elevation.

In the implementation process, the first reference surface and the second reference surface are determined through a plurality of index values according to the drilling, and the actual distribution condition of the drilling which is more fit for the determination of the first reference surface and the second reference surface can be ensured because a plurality of indexes of the drilling are actual index data of the drilling, so that the accuracy of the partition is improved.

With reference to the second possible implementation manner of the first aspect, an embodiment of the present application provides a third possible implementation manner of the first aspect, where determining a first reference plane according to the first vertical line, the second vertical line, the cumulative height of the cavern, and the elevation of the cavern floor includes: determining a karst cave accumulated height concentration area according to the index dispersion degree of the karst cave accumulated height; determining a set value of the accumulated karst cave height in the accumulated karst cave height concentration area between the first vertical line and the second vertical line; and determining the first datum plane according to the second set value of the karst cave bottom plate elevation and the accumulated height set value of the karst cave between the first vertical line and the second vertical line.

In the implementation process, the first datum plane is determined according to the elevation of the karst cave base plate of the drilled hole and the accumulated height of the karst cave, and the actual distribution condition of the drilled hole can be ensured to be more fit for the determination of the first datum plane because the elevation of the karst cave base plate of the drilled hole and the accumulated height of the karst cave are actual index data of the drilled hole, so that the accuracy of the partition is improved.

With reference to the third possible implementation manner of the first aspect, the present example provides a fourth possible implementation manner of the first aspect, where the determining a second reference plane according to the aperture elevation includes: determining the target orifice elevation with the numerical value change of the adjacent orifice elevation larger than a set change value; determining the adjacent orifice elevations of the target orifice elevations, wherein the numerical change of the adjacent orifice elevations is larger than a set value for the first time under a certain sequence; and determining a second datum plane according to the aperture elevation with smaller value in the adjacent aperture elevations with the numerical value change being larger than the set value for the first time.

In the implementation process, the first datum plane is determined according to the numerical value change of the elevation of the hole opening of the drilled hole, and the numerical value change of the elevation of the hole opening of the drilled hole is the actual index data change of the drilled hole, so that the actual distribution condition of the drilled hole which is more fit for the determination of the second datum plane can be guaranteed, and the accuracy of the subarea is improved.

With reference to the fourth possible implementation manner of the first aspect, this application provides a fifth possible implementation manner of the first aspect, where the partitioning a target plane by the first reference plane and the second reference plane includes: projecting the first datum plane and the second datum plane into a geological longitudinal section map so as to divide the geological longitudinal section map into three regions, wherein the geological longitudinal section map is a relational map of all drill holes of the target region projected in a third target coordinate system according to the drill hole distance; and matching the drilling positions in the geological longitudinal section map with the drilling positions in the litholytic geological plane map so as to partition the target plane.

In the implementation process, the drilling of the geological longitudinal section diagram is partitioned according to the actual situation of the elevation through the relation between the elevation of the drilling index numerical value and the elevation of the datum plane, and the partitioning of the drilling is performed through the actual index numerical value of the drilling, so that the partitioning of the geological longitudinal section diagram is more fit for the actual situation. And corresponding drilling positions of the classified drilling holes in the longitudinal geological section map in the karst geological plan map, so that the drilling holes are subjected to plane partitioning, and then the karst plane is partitioned to determine the development condition of the karst.

With reference to the fifth possible implementation manner of the first aspect, an embodiment of the present application provides a sixth possible implementation manner of the first aspect, where the determining a drilling number accumulated value according to distribution characteristics of karst section limestone accumulated heights and limestone top surface elevation slope change characteristics in a drilling karst single index relational graph includes: determining a target development section according to the distribution characteristics of the karst section limestone accumulated height; determining discrete drilling holes according to the elevation slope change characteristics of the limestone top surface; and determining the accumulated value of the drilling hole number through the drilling hole number of the target development segment and the discrete drilling holes.

In the implementation process, the accumulated value of the number of the drilled holes is determined according to the two indexes of the accumulated height of the limestone of the karst section and the elevation of the top surface of the limestone of the drilled holes, so that the accuracy of the accumulated value of the number of the drilled holes is ensured, and the accuracy of the zoning is improved.

In a second aspect, an embodiment of the present application further provides a planar partitioning apparatus for a development degree of a karst, including: a first determination module: the device is used for determining the accumulated value of the number of drill holes according to the distribution characteristic of the karst section limestone accumulated height and the elevation slope change characteristic of the top surface of limestone in the drill hole karst single index relation graph, and the drill hole karst single index relation graph is a relation graph in which a plurality of single index values of the drill holes are projected in a first target coordinate system according to a set sequence; a second determination module: the karst drilling index scatter diagram is a relational diagram which takes the elevation of a karst cave bottom plate as a reference according to a set sequence and projects all index values of the drill holes into a second target coordinate system, and the first perpendicular line and the second perpendicular line are used for dividing the karst drilling index scatter diagram into a plurality of areas according to the distribution condition of all index data of the drill holes; a third determination module: the first datum plane and the second datum plane are determined according to the first perpendicular line, the second perpendicular line and index values in the karst drilling index scatter diagram, and are partition elevation lines for partitioning the karst of the target area on the geological vertical section; a partitioning module: the first reference surface and the second reference surface are used for partitioning a target plane.

In a third aspect, an embodiment of the present application further provides an electronic device, including: a processor, a memory storing machine-readable instructions executable by the processor, the machine-readable instructions, when executed by the processor, performing the steps of the method of the first aspect described above, or any possible implementation of the first aspect, when the electronic device is run.

In a fourth aspect, the present application further provides a computer-readable storage medium, on which a computer program is stored, where the computer program, when executed by a processor, performs the steps of the above-mentioned method for plane partitioning of the degree of development of a karst, or any possible implementation manner of the first aspect.

In order to make the aforementioned objects, features and advantages of the present application more comprehensible, embodiments accompanied with figures are described in detail below.

Drawings

In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are required to be used in the embodiments will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present application and therefore should not be considered as limiting the scope, and for those skilled in the art, other related drawings can be obtained from the drawings without inventive effort.

FIG. 1 is a flow chart of a planar zoning method of karst development degree provided by the embodiment of the application;

FIG. 2 is a borehole karst single-index relationship diagram provided in an embodiment of the present application;

fig. 3 is a karst borehole index scattergram provided in an embodiment of the present application;

FIG. 4 is a flowchart of the planar zoning method step 202 of the development degree of the karst provided by the embodiment of the application;

FIG. 5 is a longitudinal geological section provided by an embodiment of the present application;

FIG. 6 is a litholytic geological plan provided by an embodiment of the present application;

FIG. 7 is a functional block diagram of a planar partitioning apparatus for the development degree of a karst according to an embodiment of the present disclosure;

fig. 8 is a block diagram of an electronic device according to an embodiment of the present application.

Detailed Description

The technical solution in the embodiments of the present application will be described below with reference to the drawings in the embodiments of the present application.

It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, it need not be further defined or explained in subsequent figures. Meanwhile, in the description of the present application, the terms "first", "second", and the like are used only for distinguishing the description, and are not to be construed as indicating or implying relative importance.

In engineering construction, the existence of karst can seriously affect the safety of engineering. Therefore, the real development condition of the karst in the geology is generally considered in the engineering design and construction process. For example, when the karst is mainly a shallow karst groove, if a pile foundation is adopted, not only the karst condition of the pile end needs to be considered, but also the influence of the slope of the karst groove on the stability of the pile foundation needs to be analyzed; when the karst is mainly a deep erosion cave or pipeline, the position of the bottom plate of the karst cave is mainly considered. The exposed karst is generally visual and clear in shape, such as hilly depression, karst funnels, caves and the like. However, current exploration means and analysis methods are difficult to characterize the development of the karst.

In the investigation and design process of each stage such as feasibility research, preliminary design and construction drawing of line engineering, drilling is the most extensive and commonly used means for karst investigation application, and abundant drilling data are accumulated in the engineering construction of karst areas. However, the application analysis of the current karst region drilling data is quite limited, and the application analysis is mainly used for marking the karst cave position in a geological profile and counting the probability problems of karst development such as the hole encountering rate and the karst rate of a drill hole.

In view of this, the inventor of the present application proposes a planar partitioning method for a karst development degree, in which index values of drill holes are projected in a constructed specific coordinate system and subjected to a series of analyses, so as to partition the drill holes, thereby implementing partitioning on a karst geology and further determining a karst development condition.

To facilitate understanding of the present embodiment, the following describes in detail the implementation process of the planar partitioning method of the development degree of the karst through several embodiments.

For ease of understanding, some terms used in the examples of this application are first explained:

elevation: the distance in the vertical direction from a point to the reference sea level, i.e. the absolute elevation.

Height: the ground or reference surface is facing up to a distance in the vertical direction, i.e. the relative height.

Cumulative height of karst cave: the drilling reveals the sum of the heights of all the caverns revealed within the depth range.

Height of limestone in karst section: and in the depth range of the exposed drill hole, removing the sum of the heights of all the remaining limestone sections after the accumulated height of the karst cave is removed.

Please refer to fig. 1, which is a flowchart illustrating a planar partitioning method for karst development degree according to an embodiment of the present disclosure. The specific process shown in FIG. 1 will be described in detail below.

Step 201, determining a drilling number accumulated value according to distribution characteristics of karst section limestone accumulated height and elevation slope change characteristics of a limestone top surface in a drilling karst single index relation graph.

The drill hole karst single-index relational graph is a relational graph in which a plurality of single-index values of the drill hole are projected in a first target coordinate system according to a set sequence.

As shown in fig. 2, the abscissa of the first target coordinate system may be the number of drilled holes, the primary ordinate may be the elevation, and the secondary ordinate may be the height. The number of holes is here the number of holes drilled. Illustratively, a single index value corresponding to an abscissa of 5 indicates that this is the single index value for the 5 th borehole. It is understood that the number of boreholes is intuitively reflected according to the value of the abscissa, for example, 5 on the abscissa indicates 5 on the number of boreholes between the coordinate point and the origin of coordinates, and 10 on the abscissa indicates 10 on the number of boreholes between the coordinate point and the origin of coordinates.

It is understood that the single index value of the drilling hole may include one or more of the index values of the hole opening elevation, the limestone top surface elevation, the karst cave top plate elevation, the karst cave bottom plate elevation, the cumulative karst cave height and the cumulative karst section limestone height. And each index of the single index values of the multiple drill holes is not related and is projected in the first target coordinate system according to a set sequence. The setting order may be in the order of the index value from small to large, or in the order of the index value from large to small. The setting sequence can be adjusted correspondingly according to the actual situation, and the application is not limited.

Alternatively, the distribution characteristic of the karst section limestone accumulation height may be a characteristic of numerical change of the karst section limestone accumulation height, or may be a characteristic of data distribution of the karst section limestone accumulation height. The slope change characteristic of the elevation of the limestone top surface can be the slope of the numerical change of the elevation of the limestone top surface, and can also be the distribution condition of the elevation of the limestone top surface and the like.

The accumulated value of the number of the drill holes can be the accumulated value of the number of the drill holes meeting the distribution characteristic condition of the accumulated height of the limestone in the karst section and the number of the drill holes meeting the characteristic condition of the change of the elevation slope of the top surface of the limestone.

And 202, determining a first perpendicular line and a second perpendicular line according to the accumulated value of the number of drill holes and the number of drill holes with a set ratio in the karst drill hole index scatter diagram.

The karst drilling hole index scatter diagram is a relational diagram which takes the elevation of the karst cave bottom plate as a reference according to a set sequence and projects all index values of the drill holes corresponding to the elevation of the karst cave bottom plate in a second target coordinate system. Namely, after the elevations of the karst cave bottom plate are arranged according to a set sequence, the rest index values of all the drill holes are correspondingly projected into a relational graph in a second target coordinate system. All the index values of the borehole are correlated.

As shown in fig. 3, the abscissa of the second target coordinate system may be the number of drilled holes, the primary ordinate may be the elevation, and the secondary ordinate may be the height. It is to be understood that the first target coordinate system and the second target coordinate system of the embodiment of the present application may be the same coordinate system, or may be different coordinate systems.

If the first target coordinate system and the second target coordinate system are the same coordinate system, after the accumulated value of the drilling number is determined in the drilling karst single index relational graph, all index values of the drilling holes corresponding to the elevation of the karst cave bottom plate are projected in the coordinate system by taking the elevation of the karst cave bottom plate as a reference according to a set sequence in the coordinate system so as to establish a karst drilling index scatter diagram.

If the first target coordinate system and the second target coordinate system are different coordinate systems, after the accumulated value of the number of drill holes is determined in the drill hole karst single index relational graph, all index values of the drill holes corresponding to the elevation of the karst cave bottom plate are projected in the coordinate system by taking the elevation of the karst cave bottom plate as a reference according to a set sequence in the second target coordinate system so as to establish a karst drill hole index scatter diagram. Or the drill hole karst single index relation graph and the karst drill hole index scatter diagram are simultaneously established, and after the accumulated value of the number of drill holes is determined in the drill hole karst single index relation graph, the first perpendicular line and the second perpendicular line are directly determined in the karst drill hole index scatter diagram.

The first perpendicular and the second perpendicular are used for dividing the karst borehole index scattergram into a plurality of areas according to the distribution of index data of all boreholes. The first vertical line and the second vertical line are two straight lines perpendicular to the abscissa of the second target coordinate system, the first vertical line can be determined by the accumulated value of the number of drilled holes, and the second vertical line can be determined by the number of drilled holes with the elevation setting ratio of the karst cave bottom plate.

And step 203, determining a first datum plane and a second datum plane according to the first perpendicular line, the second perpendicular line and the index numerical value in the karst drilling index scatter diagram.

The index value can be a value related to the elevation index of the karst cave floor. The index value may include: the total height of the karst cave, the elevation of the bottom plate of the karst cave, the elevation of the hole opening and other numerical values. The first reference surface and the second reference surface are partitioned elevation lines of the karst of the target area partitioned on the geological profile. The first reference plane and the second reference plane may each be defined by two straight lines perpendicular to the ordinate of the second target coordinate system. It is to be understood that the first reference surface may be a shallow erosion reference surface of the valley region and the second reference surface may be a shallow erosion reference surface of the hill front slope region.

And step 204, partitioning the target plane through the first reference plane and the second reference plane.

It will be appreciated that the geological profile of the target region may be divided into three regions by partitioning the geological profile by the first and second datum.

In the implementation process, the acquired index values of the drill holes are processed and analyzed in the coordinate system, so that different types of index relation graphs are obtained. Corresponding processing values can be obtained based on different index relation graphs, and the shallow erosion reference surface of the valley region and the shallow erosion reference surface of the mountain front slope region can be determined through the processing values, so that the karst plane can be partitioned, and the karst development condition can be further determined.

In one possible implementation, step 201 includes: determining a target development section according to distribution characteristics of the karst section limestone accumulated height, determining discrete drilling holes according to elevation slope change characteristics of the limestone top surface, and determining the accumulated value of the drilling holes according to the number of the drilling holes of the target development section and the discrete drilling holes.

Understandably, two straight lines of weak development and strong development can be fitted according to the distribution characteristics of the cumulative heights of the limestone in the karst section. Wherein, the slope of the karst section limestone accumulated height change straight line is larger, which is a strong development section. The slope of the cumulative height change straight line of the limestone at the karst section is smaller, and the karst section is a weak development section.

The target development section can be a strong development section of the cumulative height of the limestone of the karst section, and the number of drill holes of the strong development section is counted after the strong development section of the cumulative height of the limestone of the karst section is determined.

And if the elevation of part of the limestone top surface has the obvious phenomenon of changing according to the slope of the elevation of the limestone top surface which is not developed according to the karst, the drill holes corresponding to the elevation of the limestone top surface which does not change according to the regular slope are discrete drill holes, and the number of the discrete drill holes is counted. The discrete number of boreholes reflects the subsurface development overboard situation and the strength of the karst development.

In one embodiment, if there is no large difference between the change in the elevation of the limestone top surface and the change in the slope of the elevation of the limestone top surface, the discrete number of boreholes need not be counted.

After the number of drill holes in the strong development section of the karst section limestone accumulated height and the discrete drill holes in the elevation of the top surface of the limestone are determined, the number of the drill holes and the discrete drill holes are accumulated to obtain an accumulated value of the number of the drill holes.

Taking fig. 2 as an example, it is shown that if the slope of the right-side straight line of the cumulative height of the limestone in the karst section is larger, the right-side straight line is determined to be a strong development section, and the number of the drill holes in the section is counted to be 22. According to the distribution characteristics of the elevation of the top surface of the limestone in the graph, if 2 discrete drill holes are obviously formed on the leftmost side of the elevation distribution of the top surface of the limestone, the number of the discrete drill holes in the elevation of the top surface of the limestone is determined to be 2. And accumulating the number of the drill holes and the discrete drill holes in the elevation of the limestone top surface to obtain 24 accumulated values of the number of the drill holes.

In the implementation process, the accumulated value of the number of the drill holes is determined according to two indexes of the accumulated height of the limestone of the karst section of the drill holes and the elevation of the top surface of the limestone, so that the accuracy of the accumulated value of the number of the drill holes is ensured, and the accuracy of zoning is improved.

In one possible implementation, step 202, as shown in FIG. 4, includes steps 2021-2023.

Step 2021, determine a first vertical line from the accumulated value of the number of boreholes.

Here, step 2021 specifically is: and after the accumulated value of the number of the drill holes is determined, determining a perpendicular line perpendicular to the abscissa of the second target coordinate system in the karst drill hole index scatter diagram according to the accumulated value of the number of the drill holes, wherein the perpendicular line of the accumulated value of the number of the drill holes as the abscissa is the first perpendicular line.

Step 2022, determining the number of the drilling holes between the first set value of the elevation of the karst cave bottom plate and the first perpendicular line, wherein the drilling holes satisfy the set ratio.

Optionally, the first set value may be a maximum value of a karst cave floor elevation, may also be a minimum value of the karst cave floor elevation, and may also be any other intermediate value. The set ratio can be directly obtained from the outside or obtained by specific calculation. The set ratio can be 90%, 80%, 50%, 30%, etc., and the specific value of the set ratio can be adjusted according to the actual situation, which is not limited in the present application.

For example, if the first set value is the maximum elevation value of the karst cave floor, step 2022 specifically includes: and in the karst drilling index scatter diagram, determining all drilling numbers from the maximum elevation value of the karst cave bottom plate to the first vertical line, multiplying the drilling numbers by a set ratio on the basis of all the drilling numbers, and determining the drilling numbers meeting the set ratio. For example, if the set ratio is 90%, and it is determined that all the drill holes between the maximum elevation of the karst cave floor and the first vertical line are 100, the number of drill holes satisfying the set ratio is 90. If the set ratio is 50%, and all the drilling holes between the maximum elevation value of the karst cave bottom plate and the first vertical line are determined to be 100, the drilling holes meeting the set ratio are 50.

Optionally, step 2022 may further include: and determining the number of the drilling holes which meet a set ratio from the first set value of the elevation of the limestone top surface to the first vertical line, or determining the number of the drilling holes which meet the set ratio from the first set value of the elevation of the drilling air interface to the first vertical line, and the like.

Step 2023, determining a second perpendicular line according to the number of drilled holes satisfying the set ratio.

Specifically, step 2023 includes: counting the elevation of the karst cave base plate from the first set value of the elevation of the karst cave base plate to the first vertical line, and determining the elevation of the karst cave base plate close to the first vertical line in all the elevations of the karst cave base plate meeting the drilling number. And determining a second perpendicular line according to the elevation of the karst cave bottom plate close to the first perpendicular line in all the elevations of the karst cave bottom plates meeting the drilling number. For example, if the first set value is the maximum value, after the number of drilled holes meeting the set ratio is determined, counting towards the first vertical line from the maximum value of the height of the karst cave floor, and when the last value of the number of drilled holes meeting the set ratio is counted, determining the last value as the height of the karst cave floor close to the first vertical line in all the heights of the karst cave floor meeting the number of drilled holes. Taking fig. 3 as an example, if the maximum elevation value of the karst cave floor is 285m, the corresponding drilling number is 87, and the drilling number corresponding to the first vertical line is 24, it can be determined that all the drilling numbers from the maximum elevation value of the karst cave floor to the first vertical line are 63. If the set ratio is 90%, the number of drilled holes that satisfy the set ratio is 57 accordingly. And counting the maximum elevation value of the karst cave base plate to the first vertical line, wherein the 57 th drilling hole corresponds to the karst cave base plate with the elevation of 231m, and determining that the 231m is the karst cave base plate elevation close to the first vertical line in all the karst cave base plate elevations meeting the drilling hole number.

Understandably, after determining the karst cave bottom plate elevation close to the first perpendicular line in all the karst cave bottom plate elevations meeting the drilling number, the second perpendicular line can be directly determined according to the karst cave bottom plate elevation. For example, if the height of the karst cave floor close to the first vertical line among all the heights of the karst cave floor determined to satisfy the drilling number is 231m, the abscissa corresponding to 231m is determined to be 31. At this time, a perpendicular line perpendicular to the abscissa axis, which is the second perpendicular line, is determined at the point of the abscissa axis 31. And if the karst cave bottom plate elevation close to the first vertical line in all the karst cave bottom plate elevations meeting the drilling number is determined to be 239m, determining that the abscissa corresponding to the 239m is 36. At this time, a perpendicular line perpendicular to the abscissa axis, which is the second perpendicular line, is determined at the point of the abscissa of 36. In the implementation process, the first perpendicular line and the second perpendicular line are determined in the karst drilling index scatter diagram, each index of a drill hole in the karst drilling index scatter diagram can be partitioned, the actual distribution situation of each index of the drill hole is convenient to analyze, and the accuracy of plane partitioning of the karst development degree is improved.

In one possible implementation, step 203 includes: and determining a first datum plane according to the first perpendicular line, the second perpendicular line, the accumulated height of the karst cave and the elevation of the bottom plate of the karst cave, and determining a second datum plane according to the elevation of the hole opening.

The elevation of the cavern floor here means the distance of the lowermost cavern floor to the reference sea level in the vertical direction. The orifice elevation represents the distance of the orifice from the reference sea level in the vertical direction.

It is understood that embodiments of the present application may also include more reference surfaces in addition to the first and second reference surfaces.

In a possible implementation manner, step 203 specifically includes: determining a karst cave accumulated height concentration area according to the index dispersion degree of the karst cave accumulated height, determining a karst cave accumulated height set value in the karst cave accumulated height concentration area between the first vertical line and the second vertical line, and determining a first datum plane according to a second karst cave bottom plate elevation set value and a karst cave accumulated height set value between the first vertical line and the second vertical line.

The cumulative height of the karst cave is the cumulative sum of the heights of a plurality of karst caves.

It is understood that the karst cave accumulated height setting value can be a maximum karst cave accumulated height value or a minimum karst cave accumulated height value. The second set value of the karst cave bottom plate elevation can be the maximum value of the karst cave bottom plate elevation or the minimum value of the karst cave bottom plate elevation. The above setting value can be adjusted according to actual conditions, and the application is not particularly limited.

The karst cave accumulated height concentration area is an area where indexes of the karst cave accumulated height are dispersed and concentrated.

If the set value of the cumulative height of the karst cave is the minimum value of the cumulative height of the karst cave, and the second set value of the elevation of the bottom plate of the karst cave is the minimum value of the elevation of the bottom plate of the karst cave, after the cumulative height concentration area of the karst cave is determined, the number of drill holes meeting the minimum value of the cumulative height of the karst cave in the cumulative height concentration area of the karst cave between the first vertical line and the second vertical line is further determined, and the elevation value corresponding to the first datum plane is determined according to the minimum value of the elevation of the bottom plate of the karst cave corresponding to the number of drill holes. Taking fig. 3 as an example, in the diagram, the number of drill holes meeting the minimum value of the cumulative height of the karst cave in the karst cave cumulative height concentration region between the first vertical line and the second vertical line is 26, the minimum value of the elevation of the karst cave bottom plate corresponding to the number of drill holes is 222.6m, and the elevation of the corresponding first datum plane is 222.6 m.

In a possible implementation manner, step 203 specifically includes: determining the target orifice elevation with the numerical change of the adjacent orifice elevations larger than the set change value, determining the adjacent orifice elevation with the numerical change larger than the set value for the first time in the target orifice elevations in a certain sequence, and determining a second reference surface according to the orifice elevation with the smaller numerical value in the adjacent orifice elevations with the numerical change larger than the set value for the first time.

The set change value here may be an elevation difference that sets the elevation of adjacent apertures, may be a set adjacent aperture elevation change rate, or the like. The target orifice elevation may include one or more orifice elevations. The above-mentioned certain sequence may be from large to small, or from small to large.

For example, in fig. 3, before the number of boreholes is 65, the numerical value of the aperture elevation changes within a certain range (the certain range may be regarded as a set value), and when the number of boreholes is 65, 69, 74, 77, 83, 87, the numerical value of the aperture elevation changes suddenly (the sudden change is regarded as exceeding the set value), and the sudden aperture elevations may be regarded as target aperture elevations. Further, it is determined that a first numerical variation of the target aperture elevations is larger than the set aperture elevation, and if the adjacent aperture elevation with the first numerical variation larger than the set aperture elevation is an aperture elevation with an elevation of 251.8m, the corresponding second datum plane elevation is 251.8 m.

In one possible implementation, step 204 includes: and projecting the first datum plane and the second datum plane into a geological longitudinal section map so as to divide the geological longitudinal section map into three areas, and matching the drilling positions in the geological longitudinal section map with the drilling positions in the litholytic geological plane map so as to partition the target plane.

The geological profile here is a plot of all boreholes of the target area as projected in the third target coordinate system in terms of borehole distance. The abscissa of the third target coordinate system may be the borehole distance and the primary ordinate may be elevation. The borehole distance here is the actual distance between the individual boreholes. In the geological longitudinal section map, the ordinate values corresponding to the same abscissa are the index values of the same borehole.

After the first datum plane and the second datum plane are determined, elevation values corresponding to the first datum plane and the second datum plane are projected on the longitudinal geological section map, and the longitudinal geological section map can be divided into three areas by the first datum plane and the second datum plane. For example, as shown in fig. 5, if the determined elevations of the first reference surface and the second reference surface are 222.6m and 251.8m respectively, two first reference surface elevation lines and two second reference surface elevation lines perpendicular to the ordinate are made in the geological longitudinal section map with the elevations of 222.6m and 251.8m respectively. The elevation line of the first datum plane and the elevation line of the second datum plane divide the geological longitudinal section map into three areas.

It is understood that the above karst partition can be classified into karst class a, karst class B1, and karst class B2. It can be further determined that the area to which the drill hole with the elevation value larger than the second datum plane elevation line belongs to the karst type B2, the area to which the drill hole with the elevation value smaller than the second datum plane elevation line and larger than the first datum plane elevation line belongs to the karst type B1, and the area to which the drill hole with the elevation value smaller than the first datum plane elevation line belongs to the karst type A.

After the geological longitudinal section map is divided into three areas, the drilling hole positions in the geological longitudinal section map are matched with the drilling hole positions in the litholytic geological plane map respectively so as to determine the partition condition of the target plane. The litholytic geological plan is a litholytic geological actual plan. For example, after the longitudinal geological section map is partitioned, the classification condition of each drill hole can be respectively determined, and accordingly, the litholytic geological plan map can be partitioned according to the classification condition of each drill hole. Taking fig. 5 as an example, in fig. 5, the drill holes with the drill hole distance of 0-350m belong to category B1, the drill holes with the drill hole distance of 350-. Correspondingly, all the drill holes are in one-to-one correspondence with the positions of the drill holes in the litholytic geological plan according to the classification. As shown in FIG. 6, the drill holes with the drilling distance of 0-350m are classified into B1 category, the drilling distance of 350-850m is classified into A category and B1 category according to actual conditions, and the drilling distance of 850-1200m is classified into B2 category.

In the implementation process, the drilling of the geological longitudinal section diagram is partitioned according to the actual condition of elevation through the relation between the elevation of the drilling index numerical value and the elevation of the datum plane, the partitioning of the drilling is performed through the actual index numerical value of the drilling, and the fact that the partitioning of the geological longitudinal section diagram is more fit with the actual condition is guaranteed. And corresponding drilling positions of the drilling holes classified in the geological longitudinal section diagram in the karst geological plan diagram are corresponded, so that the planar partitioning of the drilling holes is realized, and then the planar partitioning of the karst is realized, so that the development condition of the karst is determined.

Based on the same application concept, a planar partitioning device of the karst development degree corresponding to the planar partitioning method of the karst development degree is further provided in the embodiment of the present application, and as the principle of solving the problem of the device in the embodiment of the present application is similar to the planar partitioning method of the karst development degree in the embodiment of the present application, the implementation of the device in the embodiment of the present application may refer to the description in the embodiment of the method, and repeated parts are not described again.

Please refer to fig. 7, which is a functional block diagram of a planar partitioning apparatus for karst development degree according to an embodiment of the present disclosure. The modules in the planar partitioning apparatus for the development degree of the karst in the present embodiment are used for executing the steps in the above method embodiments. The plane partitioning device for the development degree of the karst comprises a first determining module 301, a second determining module 302, a third determining module 303 and a partitioning module 304; wherein the content of the first and second substances,

the first determining module 301 is configured to determine a drilling hole number accumulated value according to distribution characteristics of karst section limestone accumulated heights and limestone top surface elevation slope change characteristics in a drilling hole karst single index relational graph, where the drilling hole karst single index relational graph is a relational graph in which a plurality of single index values of drilling holes are projected in a first target coordinate system according to a set order.

The second determining module 302 is configured to determine a first perpendicular and a second perpendicular according to the cumulative value of the number of drill holes and the number of drill holes in a karst drill hole index scattergram, where the karst drill hole index scattergram is a relational graph obtained by projecting all index values of drill holes in a second target coordinate system based on a set sequence of karst cave floor elevations, and the first perpendicular and the second perpendicular are used to divide the karst drill hole index scattergram into a plurality of areas according to distribution of all index data of the drill holes.

The third determination module 303 is configured to determine a first reference plane and a second reference plane according to the first perpendicular line, the second perpendicular line, and the index value in the karst drilling index scattergram, where the first reference plane and the second reference plane are partition elevation lines on the geological profile for the target area karst.

A partitioning module 304, configured to partition the target plane by the first reference plane and the second reference plane.

In a possible implementation, the second determining module 302 is further configured to determine the first perpendicular line according to the accumulated value of the drilling numbers; determining the number of drilling holes between the first set value of the karst cave bottom plate elevation and the first vertical line, wherein the drilling holes meet the set ratio; counting the elevation of the karst cave bottom plate from the first set value of the elevation of the karst cave bottom plate to the first vertical line, and determining the elevation of the karst cave bottom plate close to the first vertical line in all the elevations of the karst cave bottom plate meeting the drilling number; and determining a second perpendicular line according to the elevation of the karst cave bottom plate close to the first perpendicular line in all the elevations of the karst cave bottom plates meeting the drilling number.

In a possible implementation manner, the third determining module 303 is further configured to determine a first reference plane according to the first vertical line, the cumulative height of the karst cave, and the elevation of the karst cave bottom plate; a second datum is determined from the aperture elevation.

In a possible implementation manner, the third determining module 303 is specifically configured to determine a karst cave accumulated height concentration area according to the karst cave accumulated height; determining a set value of the accumulated karst cave height in the accumulated karst cave height concentration area between the first vertical line and the second vertical line; and determining the first datum plane according to the second set value of the karst cave bottom plate elevation and the accumulated height set value of the karst cave between the first vertical line and the second vertical line.

In a possible implementation, the third determining module 303 is specifically configured to determine a target aperture elevation where a numerical change of adjacent aperture elevations is greater than a set change value; determining the adjacent orifice elevations of the target orifice elevations, wherein the numerical change of the adjacent orifice elevations is larger than a set value for the first time under a certain sequence; and determining a second datum plane according to the hole opening elevation with smaller value in the adjacent hole opening elevations with the numerical value change larger than the set value for the first time.

In a possible implementation, the partitioning module 304 is further configured to project the first reference plane and the second reference plane into a geological longitudinal section map, so as to divide the geological longitudinal section map into three regions, where the geological longitudinal section map is a relational map of all drill holes of the target region projected in a third target coordinate system according to the drill hole distance; and matching the drilling positions in the geological longitudinal section map with the drilling positions in the litholytic geological plane map so as to partition the target plane.

In a possible implementation manner, the first determining module 301 is further configured to determine a target development segment according to a distribution characteristic of the cumulative height of the limestone of the karst segment; determining discrete drilling holes according to the elevation slope change characteristics of the limestone top surface; and determining the accumulated value of the drilling hole number through the drilling hole number of the target development segment and the discrete drilling holes.

To facilitate understanding of the present embodiment, an electronic device for performing planar partitioning of the degree of development of a karst as disclosed in the embodiments of the present application will be described in detail below.

Fig. 8 is a block diagram of an electronic device. The electronic device 100 may include a memory 111, a memory controller 112, a processor 113, a peripheral interface 114, an input-output unit 115, and a display unit 116. It will be understood by those skilled in the art that the structure shown in fig. 8 is merely an illustration and is not intended to limit the structure of the electronic device 100. For example, electronic device 100 may also include more or fewer components than shown in FIG. 8, or have a different configuration than shown in FIG. 1.

The above-mentioned elements of the memory 111, the memory controller 112, the processor 113, the peripheral interface 114, the input/output unit 115 and the display unit 116 are electrically connected to each other directly or indirectly, so as to implement data transmission or interaction. For example, the components may be electrically connected to each other via one or more communication buses or signal lines. The processor 113 is used to execute the executable modules stored in the memory.

The Memory 111 may be, but is not limited to, a Random Access Memory (RAM), a Read Only Memory (ROM), a Programmable Read-Only Memory (PROM), an Erasable Read-Only Memory (EPROM), an electrically Erasable Read-Only Memory (EEPROM), and the like. The memory 111 is configured to store a program, and the processor 113 executes the program after receiving an execution instruction, and the method executed by the electronic device 100 defined by the process disclosed in any embodiment of the present application may be applied to the processor 113, or implemented by the processor 113.

It is understood that the memory 111 may also be used to store index values of the first target coordinate system, the second target coordinate system, the third target coordinate system, and the drill holes, such as an elevation of an opening, an elevation of a top surface of the limestone, an elevation of a top plate of the karst cave, an elevation of a bottom plate of the karst cave, an accumulated height of the karst cave, and an accumulated height of the limestone of the karst section.

The processor 113 may be an integrated circuit chip having signal processing capability. The Processor 113 may be a general-purpose Processor, and includes a Central Processing Unit (CPU), a Network Processor (NP), and the like; the Integrated Circuit may also be a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA) or other programmable logic device, a discrete gate or transistor logic device, or a discrete hardware component. The various methods, steps, and logic blocks disclosed in the embodiments of the present application may be implemented or performed. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.

The peripheral interface 114 couples various input/output devices to the processor 113 and memory 111. In some embodiments, the peripheral interface 114, the processor 113, and the memory controller 112 may be implemented in a single chip. In other examples, they may be implemented separately from the individual chips.

The input/output unit 115 is used to provide input data to the user. The input/output unit 115 may be, but is not limited to, a mouse, a keyboard, and the like.

It is understood that the input/output unit 115 may be configured to input index values of the drilling hole, such as an elevation of an orifice, an elevation of a limestone top surface, an elevation of a karst cave top plate, an elevation of a karst cave bottom plate, an accumulated height of the karst cave, and an accumulated height of the limestone in the karst section, and determination values of the method, such as a first set value, a second set value, and a set variation value.

The display unit 116 provides an interactive interface (e.g., a user operation interface) between the electronic device 100 and the user or is used for displaying image data to the user for reference. In this embodiment, the display unit may be a liquid crystal display or a touch display. In the case of a touch display, the display can be a capacitive touch screen or a resistive touch screen, which supports single-point and multi-point touch operations. The support of single-point and multi-point touch operations means that the touch display can sense touch operations simultaneously generated from one or more positions on the touch display, and the sensed touch operations are sent to the processor for calculation and processing.

It is understood that the display unit 116 may be used to display a karst single index relational map, a karst borehole index scatter plot, a geological profile, a karst geological plan, and the like.

The electronic device 100 in this embodiment may be configured to perform each step in each method provided in this embodiment.

Furthermore, an embodiment of the present application further provides a computer-readable storage medium, where a computer program is stored on the computer-readable storage medium, and when being executed by a processor, the computer program performs the steps of the planar partitioning method for the development degree of the karst described in the above method embodiment.

The computer program product of the planar partition method for the development degree of a karst provided in the embodiment of the present application includes a computer-readable storage medium storing a program code, where instructions included in the program code may be used to execute the steps of the planar partition method for the development degree of a karst described in the above method embodiment, which may be specifically referred to in the above method embodiment and are not described herein again.

In the embodiments provided in the present application, it should be understood that the disclosed apparatus and method can be implemented in other ways. The apparatus embodiments described above are merely illustrative, and for example, the flowchart and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of apparatus, methods and computer program products according to various embodiments of the present application. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.

In addition, functional modules in the embodiments of the present application may be integrated together to form an independent part, or each module may exist separately, or two or more modules may be integrated to form an independent part.

The functions, if implemented in the form of software functional modules and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solutions of the present application or portions thereof that substantially contribute to the prior art may be embodied in the form of a software product, which is stored in a storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the methods described in the embodiments of the present application. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes. It is noted that, herein, relational terms such as first and second, and the like may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Also, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising … …" does not exclude the presence of other identical elements in a process, method, article, or apparatus that comprises the element.

The above description is only a preferred embodiment of the present application and is not intended to limit the present application, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present application shall be included in the protection scope of the present application. It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, it need not be further defined and explained in subsequent figures.

The above description is only for the specific embodiments of the present application, but the scope of the present application is not limited thereto, and any person skilled in the art can easily conceive of the changes or substitutions within the technical scope of the present application, and shall be covered by the scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.

Claims

1. A planar partitioning method for the development degree of a karst is characterized by comprising the following steps:

determining a drilling number accumulated value according to distribution characteristics of karst section limestone accumulated height and elevation slope change characteristics of limestone top surface in a drilling karst single index relational graph, wherein the drilling karst single index relational graph is a relational graph in which a plurality of single index values of drilling holes are projected in a first target coordinate system according to a set sequence;

determining a first perpendicular line and a second perpendicular line according to the accumulated value of the number of drill holes and the number of drill holes with a set ratio in a karst drill hole index scatter diagram, wherein the karst drill hole index scatter diagram is a relational diagram obtained by projecting all index values of drill holes corresponding to the elevation of a karst cave bottom plate in a second target coordinate system on the basis of the elevation of the karst cave bottom plate according to a set sequence, and the first perpendicular line and the second perpendicular line are used for dividing the karst drill hole index scatter diagram into a plurality of areas according to the distribution situation of all index data of the drill holes;

determining a first datum plane and a second datum plane according to the first perpendicular line, the second perpendicular line and an index numerical value in a karst drilling index scatter diagram, wherein the first datum plane and the second datum plane are partition elevation lines for partitioning a target area karst on a geological vertical section;

and partitioning the target plane through the first reference plane and the second reference plane.

2. The method of claim 1, wherein determining the first vertical and the second vertical based on the accumulated number of boreholes and a set ratio of the number of boreholes in the karst borehole index scattergram comprises:

determining a first vertical line according to the accumulated value of the drilling number;

determining the number of drilling holes between the first set value of the karst cave bottom plate elevation and the first vertical line, wherein the drilling holes meet the set ratio;

and determining a second perpendicular line according to the number of the drill holes meeting the set ratio.

3. The method of claim 1, wherein the index value comprises: the accumulated height of the karst cave, the elevation of a bottom plate of the karst cave and the elevation of an orifice of the karst cave, and determining a first datum plane and a second datum plane according to index values in a first perpendicular line, a second perpendicular line and a karst drilling index scatter diagram, wherein the method comprises the following steps:

determining a first datum plane according to the first perpendicular line, the second perpendicular line, the accumulated height of the karst cave and the elevation of the bottom plate of the karst cave;

a second datum is determined from the aperture elevation.

4. The method of claim 3, wherein determining a first datum level from the first vertical line, the second vertical line, the cumulative height of the cavern, and the elevation of the cavern floor comprises:

determining a karst cave accumulated height concentration area according to the index dispersion degree of the karst cave accumulated height;

determining a set value of the accumulated karst cave height in the accumulated karst cave height concentration area between the first vertical line and the second vertical line;

and determining the first datum plane according to the second set value of the karst cave bottom plate elevation and the set value of the karst cave accumulated height between the first vertical line and the second vertical line.

5. The method of claim 3, wherein determining a second reference plane from the aperture elevation comprises:

determining the target orifice elevation with the numerical value change of the adjacent orifice elevation larger than a set change value;

determining the adjacent orifice elevations of the target orifice elevations, wherein the numerical change of the adjacent orifice elevations is larger than a set value for the first time under a certain sequence;

and determining a second datum plane according to the hole opening elevation with smaller value in the adjacent hole opening elevations with the numerical value change larger than the set value for the first time.

6. The method of claim 1, wherein the partitioning a target plane by the first reference plane and the second reference plane comprises:

projecting the first datum plane and the second datum plane into a geological longitudinal section map so as to divide the geological longitudinal section map into three regions, wherein the geological longitudinal section map is a relational map of all drill holes of the target region projected in a third target coordinate system according to the drill hole distance;

and matching the positions of the drill holes in the geological longitudinal section map with the positions of the drill holes in the karst geological plane map so as to partition the target plane.

7. The method of claim 1, wherein determining the cumulative number of boreholes from the distribution characteristics of the cumulative height of limestone at the karst section and the elevation slope change characteristics of the limestone top surface in the karst single index relational graph of the boreholes comprises:

determining a target development section according to the distribution characteristics of the karst section limestone accumulated height;

determining discrete drilling holes according to the elevation slope change characteristics of the limestone top surface;

and determining the accumulated value of the drilling hole number through the drilling hole number of the target development segment and the discrete drilling holes.

8. A planar zoning apparatus for the degree of development of a karst, comprising:

a first determination module: the device is used for determining the accumulated value of the number of drill holes according to the distribution characteristic of the karst section limestone accumulated height and the elevation slope change characteristic of the top surface of limestone in the drill hole karst single index relation graph, and the drill hole karst single index relation graph is a relation graph in which a plurality of single index values of the drill holes are projected in a first target coordinate system according to a set sequence;

a second determination module: the karst drilling index scattergram is a relational graph which takes the elevation of a karst cave bottom plate as a reference according to a set sequence and projects all index values of the drill holes in a second target coordinate system, and the first vertical line and the second vertical line are used for dividing the karst drilling index scattergram into a plurality of areas according to the distribution condition of all index data of the drill holes;

a third determining module: the first datum plane and the second datum plane are determined according to the first perpendicular line, the second perpendicular line and index values in the karst drilling index scatter diagram, and are partition elevation lines for partitioning the karst of the target area on the geological vertical section;

a partitioning module: the first reference surface and the second reference surface are used for partitioning a target plane.

9. An electronic device, comprising: a processor, a memory storing machine-readable instructions executable by the processor, the machine-readable instructions when executed by the processor performing the steps of the method of any of claims 1 to 7 when the electronic device is run.

10. A computer-readable storage medium, having stored thereon a computer program which, when being executed by a processor, is adapted to carry out the steps of the method according to any one of claims 1 to 7.