Three-dimensional geological model storage method and device based on geological contour line

Abstract

The invention relates to the technical field of storage, in particular to a geological contour line-based three-dimensional geological model storage method, which comprises the following steps of S1, establishing a three-dimensional geological body model; s2, according to each stratum, extracting upper surface data of the stratum; s3, extracting boundary data of the upper surface of a certain stratum; s4, geological contour line data are interpolated; and S5, intersection of the geological contour line and the upper surface boundary line is carried out to obtain local layer data. According to the method, the data storage efficiency is high, the geological contour line is a simplified and compressed expression of topographic data, the data volume can be greatly reduced, and compared with a direct storage complete three-dimensional grid model, the contour line data occupies less storage space; generation and data processing are easy, the contour lines are relatively simply extracted from geological exploration and measurement data, and the contour lines can be easily extracted through existing GIS tools and software; data format standardization and vectorization are facilitated, contour line data has a high standardization degree, and data exchange and sharing between different software and platforms are facilitated.

Description

Technical Field

[0001] The invention relates to the field of storage technology, in particular to a three-dimensional geological model storage method and device based on geological contour lines. Background Art

[0002] The technical background of the 3D geological model storage method is related to the combination of geology, computer science, and data storage and management technology. The research and application in this field are developing rapidly, mainly because of its important role in railway engineering construction, railway disaster prevention and mitigation, resource exploration, environmental management and urban planning. The background technologies include:

[0003] 1. Geological knowledge and application: The core of 3D geological models is to accurately represent underground structures, which requires deep geological knowledge, such as understanding of geological features such as strata, faults, and ore bodies. Geological models enable geologists to visualize and analyze geological data in 3D space, thereby making more accurate interpretations and predictions of geological structures and geological history.

[0004] 2. Computer science and information technology: This includes 3D modeling technology, graphics processing technology, computer vision, and machine learning. Computer science provides technical support for processing large amounts of geological data, performing complex calculations, and realizing data visualization. In particular, efficient algorithms and computing architectures are crucial in the construction and optimization of 3D models.

[0005] 3. Data storage and management technology: With the increase in the amount of geological survey data, how to effectively store, manage and retrieve these big data has become one of the key technologies. This involves database management, cloud storage solutions, and data compression and optimization technology to ensure the sustainability and accessibility of data.

[0006] 4. 3D visualization technology: 3D visualization technology enables complex geological data to be displayed intuitively and provides user interaction interface design to assist scientists and engineers in better understanding geological structures and making decision support.

[0007] 5. Geographic Information System (GIS) technology: GIS is an important tool for managing, analyzing and displaying geospatial data. In three-dimensional geological models, GIS technology can support spatial analysis of geological data and integration with other types of environmental and geographic data.

[0008] The combination of these technologies not only promotes the application of 3D geological models in scientific research, but also greatly promotes their use in practical fields such as engineering construction, resource exploration and environmental management. Through the comprehensive application of these technologies, geological risk assessment, resource assessment and environmental monitoring can be carried out more effectively, which is also of great practical significance for urban planning and disaster prevention.

[0009] However, due to the huge amount of 3D geological model data, direct storage will take up a lot of storage space and lead to inefficient data access and processing. Summary of the invention

[0010] The purpose of the present invention is to provide a three-dimensional geological model storage method and device based on geological contour lines to address the problems in the prior art that the amount of three-dimensional geological model data is huge, direct storage will occupy a large amount of storage space, and lead to low data access and processing efficiency.

[0011] In order to achieve the above object, the technical solution adopted by the present invention is:

[0012] In a first aspect, the present invention provides a method for storing a three-dimensional geological model based on geological contour lines, comprising the following steps:

[0013] S1. Establish a three-dimensional geological model;

[0014] S2, extracting the upper surface data of each stratum according to the stratum;

[0015] S3, extracting the upper surface boundary data of a certain stratum;

[0016] S4, interpolated geological contour data;

[0017] S5. The geological contour lines intersect with the upper surface boundary lines to obtain local layer data.

[0018] As a preferred technical solution of the present invention, in step S1, a three-dimensional geological model is constructed using boreholes, geological profiles, geological survey maps or surveying data.

[0019] As a preferred technical solution of the present invention, in step S2, the upper surface of each stratum is composed of a triangulated network, and the stratum surface data is exported in the following format:

[0020] The first row includes the number of vertices N and the number of triangles T;

[0021] The second to N+1 rows are vertex data areas, each row is a vertex;

[0022] The format of the vertices is as follows:

[0023] Vertex index, starting from 1, vertex X coordinate, vertex Y coordinate, vertex Z coordinate;

[0024] The N+2th line to the end of the file is the triangle data area, and each line is a triangle;

[0025] The format of the triangle is as follows:

[0026] Triangle number, vertex 1 index, vertex 2 index, vertex 3 index;

[0027] The vertices of the triangle are in counter-clockwise direction.

[0028] As a further preferred technical solution of the present invention, in step S3, the process of extracting the boundary includes: start → read the input file → read the number of vertices and triangles → parse the vertex data → parse the triangle data → construct the edge from the triangle → record the number of adjacent triangles for each edge → determine whether the number of adjacent triangles is 1 →

[0029] Yes, boundary edge → output boundary and vertex coordinates → end;

[0030] No, internal edge → end.

[0031] As a further preferred technical solution of the present invention, in step S4, S is sorted from small to large in the elevation direction, and S is a set of triangles.

[0032] As a further preferred technical solution of the present invention, step S4 includes the following steps:

[0033] S41. Sort S; Sort S in ascending order according to elevation. Since a triangle has three vertices and each vertex has a Z value, the Z value space occupied by the triangle is a range, namely:

[0034] TR i =[Z min , Z max ]

[0035] in:

[0036] Z min =Min(Z i1 , Z i2 , Z i3 )

[0037] Z max =Max(Z i1 , Z i2 , Z i3 )

[0038] Min is the minimum function, Max is the maximum function, Z i1 , Z i2 , Z i3 is the Z value of the three vertices of the i-th triangle; by calculating all triangles in S, we get the TR set. First, we calculate the Z value of the three vertices of the i-th triangle. min Sort from small to large. If there are more than two triangles, the Z min The same, then the Z of these triangles max Sort in ascending order, and finally get the sorted set TSeq;

[0039] S42, calculate the contour line starting surface and set the contour interval; the contour line starting surface is based on the Z value of the first TSeq min Calculate the starting surface Z of the contour line according to the contour drawing specification b :

[0040]

[0041] in, The symbol indicates rounding down, and h is the equal height interval;

[0042] S43. Calculate the plane equation Z of each contour line based on the starting surface and contour distance. i ;

[0043] Z i =i*h+Z b

[0044] Among them, i represents the contour line number, and the serial number starts from 0;

[0045] S44, interpolate geological contours; extract triangles T from TSeq in sequence j , j is T j The sequence number in TSeq is based on T j Query the contour plane equation within the range and get the set Z sub , Z sub Each equation in the triangle T j By finding the intersection, we finally get a series of points associated with each contour plane, and connect these points in sequence to get several contour lines.

[0046] As a preferred technical solution of the present invention, in step S1, GoCad, Petrel or EarthVision modeling software is used to build a three-dimensional geological model.

[0047] In a second aspect, the present invention further provides a three-dimensional geological model storage device based on geological contour lines, comprising:

[0048] Modeling module, used to build three-dimensional geological body models;

[0049] The first acquisition module is used to extract the upper surface data of each stratum in the three-dimensional geological body model;

[0050] The second acquisition module is used to extract the upper surface boundary data of a certain stratum;

[0051] Interpolation module, used to interpolate geological contour data in upper surface boundary data;

[0052] The third acquisition module is used to intersect the geological contour lines with the upper surface boundary lines to obtain local layer data.

[0053] In a third aspect, the present invention further provides an electronic device, comprising:

[0054] a memory having a computer program stored thereon;

[0055] A processor is used to execute the program in the memory to implement the three-dimensional geological model storage method based on geological contour lines as described in any one of the above.

[0056] In a fourth aspect, the present invention further provides a computer-readable storage medium having a computer program stored thereon, wherein when the program is executed by a processor, the three-dimensional geological model storage method based on geological contour lines as described in any one of the above items is implemented.

[0057] In summary, due to the adoption of the above technical solution, the beneficial effects of the present invention are:

[0058] The three-dimensional geological model storage method based on geological contour lines described in the present invention has the characteristics of high data storage efficiency, easy generation and data processing, and facilitation of data format standardization and vectorization. In terms of data storage efficiency, geological contour lines are a simplified and compressed expression of terrain data, which can greatly reduce the amount of data. Compared with directly storing a complete three-dimensional grid model, contour line data occupies less storage space. In terms of easy generation and data processing, it is relatively simple to extract contour lines from geological exploration and measurement data, and it can be easily implemented through existing GIS tools and software. In terms of facilitation of data format standardization and vectorization, contour line data has a high degree of standardization, which is convenient for data exchange and sharing between different software and platforms. BRIEF DESCRIPTION OF THE DRAWINGS

[0059] Figure 1 A schematic diagram for extracting the upper surface of a certain stratum;

[0060] Figure 2 It is a schematic diagram of the triangular network on the upper surface of the stratum;

[0061] Figure 3 This is a schematic diagram of the triangle data format;

[0062] Figure 4 Flowchart for extracting the boundary of a surface triangle mesh;

[0063] Figure 5 Schematic diagram for extracting the boundary of the surface triangle network;

[0064] Figure 6 For the set Z sub With triangle T j Schematic diagram of intersection;

[0065] Figure 7is a schematic diagram of interpolated contour lines;

[0066] Figure 8 Schematic diagram of the intersection of geological contour lines and upper surface boundary lines. DETAILED DESCRIPTION

[0067] The present invention is further described in detail below in conjunction with test examples and specific implementation methods. However, this should not be understood as the scope of the above subject matter of the present invention being limited to the following embodiments, and all technologies realized based on the content of the present invention belong to the scope of the present invention.

[0068] Unless otherwise specified, in the description of the specific embodiments of the present invention, the terms indicating the orientation or position relationship such as "up", "down", "left", "right", "center", "inside", "outside", etc. are all expressions based on the orientation or position relationship shown in the drawings, or are the orientation or position relationship when the invented product / equipment / device is usually used. These terms of orientation or position relationship are only for the convenience of describing the scheme of the present invention or simplifying the description in the specific embodiments, so as to facilitate the technicians to quickly understand the scheme, and do not indicate or imply that a specific device / component / element must have a specific orientation, or be constructed and operated in a specific position relationship, and therefore cannot be understood as a limitation on the present invention.

[0069] In addition, if the terms "horizontal", "vertical", "overhanging", "parallel" and the like appear, it does not mean that the corresponding devices / components / elements are required to be absolutely horizontal or vertical or overhanging or parallel, but may be slightly tilted or have deviations. For example, "horizontal" only means that its direction is more horizontal than "vertical", and does not mean that the structure must be completely horizontal, but may be slightly tilted. Alternatively, it can be simplified to mean that the corresponding devices / components / elements are set in directions such as "horizontal", "vertical", "overhanging", "parallel", etc., and can have an error / deviation of ±10% relative to the corresponding direction setting, more preferably an error / deviation within ±8%, more preferably an error / deviation within ±6%, more preferably an error / deviation within ±5%, and more preferably an error / deviation within ±4%. As long as the corresponding device / component / element is within the error / deviation range, it can still achieve its role in the scheme of the present invention.

[0070] In addition, the expressions “first”, “second”, “third”, etc., which appear in the terms, are merely used to distinguish the description of the same or similar components and should not be understood as emphasizing or implying the relative importance of specific components.

[0071] In addition, in the description of the embodiments of the present invention, "several", "plurality" and "a number" represent at least 2. It can be any number such as 2, 3, 4, 5, 6, 7, 8, 9, and even more than 9.

[0072] In addition, in the description of the technical solution of the present invention, unless otherwise clearly specified / defined / restricted, the terms "set", "install", "connect", "connected", "provided with", "laid", and "arranged" should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection, and can be welding, riveting, bolting, threading, and other commonly used connection means in the field. This connection can be a mechanical connection, an electrical connection, or a communication connection; it can be a direct connection, or an indirect connection through an intermediate medium, and it can be the internal connection of two elements.

[0073] In the related art, the data volume of the three-dimensional geological model is huge. Direct storage will occupy a lot of storage space and lead to low efficiency in data access and processing. Figures 1 to 8 To elaborate.

[0074] Example 1

[0075] like Figures 1 to 8 As shown, the three-dimensional geological model storage method based on geological contour lines of the present invention comprises the following steps:

[0076] Step 1: Modeling of three-dimensional geological bodies.

[0077] 3D geological models are constructed using drill holes, geological profiles, geological survey maps or mapping data. 3D geological models can be constructed using modeling software such as GoCad, Petrel or EarthVision.

[0078] 3D geological modeling technology is the basis for creating 3D geological models. 3D modeling technology involves methods for generating 3D models from 2D maps and geological data (such as drilling data, geological survey data, seismic data, and mapping data). These models can show the morphology and spatial relationships of underground structures, helping geologists and engineers to have a deeper understanding of geological conditions.

[0079] Step 2: Extract the upper surface data of each stratum.

[0080] like Figure 1 As shown, the strata are separated and the upper surface of the strata is extracted. The upper surface of the strata can be extracted using geological modeling software. The upper surface of each stratum is composed of a triangulated network, such as Figure 2 As shown in the figure, according to the triangle network rules, the stratum surface data is exported into the following format:

[0081] The first row includes the number of vertices N and the number of triangles T;

[0082] The second to N+1 rows are vertex data areas, each row is a vertex;

[0083] The format of the vertices is as follows:

[0084] Vertex index (starting from 1), vertex X coordinate, vertex Y coordinate, vertex Z coordinate;

[0085] The N+2th line to the end of the file is the triangle data area, and each line is a triangle;

[0086] The format of the triangle is as follows:

[0087] Triangle number, vertex 1 index, vertex 2 index, vertex 3 index;

[0088] The vertices of the triangle are arranged in a counterclockwise direction, such as Figure 3 As shown, Figure 3 The data format is as follows:

[0089] 4,2#Number of vertices N, number of triangles T

[0090] 1,10,10,0#Vertex data area

[0091] 2,0,0,0

[0092] 3,20,0,0

[0093] 4,10,-10,0

[0094] 1 1 2 3#Triangle data area

[0095] 2243.

[0096] Step 3: Extract the upper surface boundary data of a certain stratum.

[0097] Since the boundaries of the strata are discontinuous, there may be cracks or holes at the upper surface boundaries of the strata, such as Figure 1 Therefore, when generating geological contour lines, they need to be limited to the boundary. The process of extracting the boundary is as follows (see Figure 4 ):

[0098] Start → Read input file → Read the number of vertices and triangles → Parse vertex data → Parse triangle data → Build edges from triangles → Record the number of adjacent triangles for each edge → Determine whether the number of adjacent triangles is 1 → Yes (boundary edge → Output boundary and vertex coordinates → End) / No (internal edge → End).

[0099] The criterion for judging the boundary is whether the number of adjacent triangles of the edge is 1. If it is 1, it means that this edge is a boundary edge, otherwise it is an internal edge. After filtering out the boundary edges, they are connected into multiple closed lines, such as Figure 5 shown.

[0100] Step 4: Interpolate geological contour data.

[0101] from Figure 2 It can be seen that the upper surface of the stratum is composed of a triangular network, which is composed of a number of triangles (triangle set S). Therefore, when interpolating contour lines, the intersection relationship between the contour plane and the triangle is mainly calculated. In order to speed up the calculation, S is sorted from small to large in the elevation direction. The main steps are as follows:

[0102] 1. Sort S.

[0103] Sort S in ascending order by elevation. Since a triangle has three vertices, each vertex has a Z value, i.e., an elevation value. Therefore, the Z value space occupied by the triangle is a range, i.e.:

[0104] TR i =[Z min , Z max ]

[0105] in:

[0106] Z min =Min(Z i1 , Z i2 , Z i3 )

[0107] Z max =Max(Z i1 , Z i2 , Z i3 )

[0108] Min is the minimum function, Max is the maximum function, Z i1 , Z i2 , Z i3 are the Z values ​​of the three vertices of the i-th triangle.

[0109] By calculating all triangles in S, we get the TR set. First, we calculate the Z in the TR set. min Sort from small to large. If there are more than two triangles, the Z min The same, then the Z of these triangles max Sort them in ascending order, and finally get the sorted set TSeq.

[0110] 2. Calculate the starting surface of the contour lines and set the contour intervals.

[0111] The starting point of the contour line should be the Z value of the first TSeq value. min But due to Z min It is very likely to be a decimal, and it is not a multiple of the contour interval, which does not meet the contour mapping specifications. Therefore, it is necessary to adjust the starting surface Z of the contour line. b .

[0112]

[0113] in, The symbol represents rounding down, and h is the equal height interval.

[0114] 3. Calculate the plane equation Z of each contour line according to the starting surface and contour distance i .

[0115] Z i =i*h+Z b

[0116] Among them, i represents the number of contour lines, and the serial number starts from 0.

[0117] 4. Interpolate geological contours.

[0118] Take out triangles T from TSeq one by one j , j is T j The sequence number in TSeq. j Query the contour plane equation within the range and get the set Z sub . sub Each equation in the triangle T j Ask for intercourse, such as Figure 6 As shown, we finally get a series of points associated with each contour plane. By connecting these points in sequence, we can get several contour lines, such as Figure 7 shown.

[0119] Geological contour vectorization technology is an important geographic information system (GIS) technology, which is mainly used to convert three-dimensional geological models into two-dimensional representations, especially through contour lines. The key purpose of this technology is to solve the challenges of storage and processing of three-dimensional geological models, especially when the data volume is large and efficient processing is required.

[0120] The rapid 3D modeling technology of geological contour lines is an efficient technology for generating 3D geological models from 2D geological contour maps. This technology is very important in the field of geological science and engineering because it can quickly provide a 3D view of the terrain, thus helping professionals to make more accurate geological analysis and decisions.

[0121] Step 5: Find the intersection of geological contour lines and upper surface boundary lines.

[0122] The geological contour line is intersected with the upper surface boundary line to obtain the final result data of the local layer, such as Figure 8 shown.

[0123] The three-dimensional geological model storage method based on geological contour lines described in this embodiment has the characteristics of high data storage efficiency, easy generation and data processing, and facilitation of data format standardization and vectorization. In terms of data storage efficiency, geological contour lines are a simplified and compressed expression of terrain data, which can greatly reduce the amount of data. Compared with directly storing a complete three-dimensional grid model, contour line data occupies less storage space. In terms of easy generation and data processing, it is relatively simple to extract contour lines from geological exploration and measurement data, and it can be easily implemented through existing GIS tools and software. In terms of facilitation of data format standardization and vectorization, contour line data has a high degree of standardization, which facilitates data exchange and sharing between different software and platforms.

[0124] Example 2

[0125] The three-dimensional geological model storage device based on geological contour lines of the present invention comprises a modeling module, an interpolation module, a first acquisition module, a second acquisition module and a third acquisition module.

[0126] The modeling module is used to build a three-dimensional geological model;

[0127] The first acquisition module is used to extract the upper surface data of each stratum in the three-dimensional geological body model;

[0128] The second acquisition module is used to extract the upper surface boundary data of a certain stratum;

[0129] The interpolation module is used to interpolate geological contour data in the upper surface boundary data;

[0130] The third acquisition module is used to intersect the geological contour lines with the upper surface boundary lines to obtain local layer data.

[0131] The modeling module, the first acquisition module, the second acquisition module, the interpolation module and the third acquisition module are connected via a bus.

[0132] Example 3

[0133] An electronic device according to the present invention comprises:

[0134] a memory having a computer program stored thereon;

[0135] The processor is used to execute the program in the memory to implement the three-dimensional geological model storage method based on geological contour lines as described in Example 1.

[0136] As an optional solution of this embodiment, the electronic device may include: a processor, a memory, and the electronic device may also include one or more of a multimedia component, an input / output (I / O) interface, and a communication component.

[0137] The processor is used to control the overall operation of the electronic device to complete all or part of the steps in the above-mentioned three-dimensional geological model storage method based on geological contour lines.

[0138] The memory is used to store various types of data to support the operation of the electronic device, which data may include, for example, instructions for any application or method operating on the electronic device, and application-related data; the memory can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic memory, flash memory, magnetic disk or optical disk.

[0139] The multimedia component may include a screen and an audio component, wherein the screen may be, for example, a touch screen, and the audio component is used to output and / or input audio signals; for example, the audio component may include a microphone for receiving external audio signals, and the received audio signals may be further stored in a memory or sent via a communication component; the audio component also includes at least one speaker for outputting audio signals.

[0140] The I / O interface provides an interface between the processor and other interface modules, which may be a keyboard, a mouse, buttons, etc.; these buttons may be virtual buttons or physical buttons.

[0141] The communication component is used for wired or wireless communication between the electronic device and other devices; wireless communication, such as Wi-Fi, Bluetooth, Near Field Communication (NFC), 2G, 3G, 4G, 5G or 5.5G, or a combination of one or more of them, so the corresponding communication component may include: Wi-Fi module, Bluetooth module, Star Flash module, NFC module, mobile phone communication module.

[0142] As an optional solution of this embodiment, the electronic device can be implemented by one or more application specific integrated circuits (ASIC), digital signal processors (DSP), digital signal processing devices (DSPD), programmable logic devices (PLD), field programmable gate arrays (FPGA), controllers, microcontrollers, microprocessors or other electronic components to execute the above-mentioned three-dimensional geological model storage method based on geological contours.

[0143] In addition, the computer-readable storage medium provided in the embodiment of the present disclosure may be the above-mentioned memory including program instructions, and the above-mentioned program instructions may be executed by a processor of an electronic device to complete the above-mentioned three-dimensional geological model storage method based on geological contour lines.

[0144] Example 4

[0145] The computer-readable storage medium described in the present invention stores a computer program thereon, and is characterized in that when the program is executed by a processor, the three-dimensional geological model storage method based on geological contour lines as described in Example 1 is implemented.

[0146] Computer-readable storage media are used to store various types of data to support operations on the electronic device, which data may include, for example, instructions for any application or method operating on the electronic device, as well as application-related data; the memory may be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic memory, flash memory, magnetic disk or optical disk.

[0147] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should be included in the protection scope of the present invention.

Claims

1. A three-dimensional geological model storage method based on geological contour lines, It is characterized in that The following steps are involved: S1. Establish a three-dimensional geological model; S2, extracting the upper surface data of each stratum according to the stratum; S3, extracting the upper surface boundary data of a certain stratum; S4, interpolated geological contour data; S5. The geological contour lines intersect with the upper surface boundary lines to obtain local layer data.

2. The three-dimensional geological model storage method based on geological contour lines according to claim 1, It is characterized in that In step S1, a three-dimensional geological model is constructed using boreholes, geological profiles, geological survey maps or surveying data.

3. The three-dimensional geological model storage method based on geological contour lines according to claim 1, It is characterized in that In step S2, the upper surface of each stratum is composed of a triangulated network, and the stratum surface data is exported in the following format: The first row includes the number of vertices N and the number of triangles T; The second to N+1 rows are vertex data areas, each row is a vertex; The format of the vertices is as follows: Vertex index, starting from 1, vertex X coordinate, vertex Y coordinate, vertex Z coordinate; The N+2th line to the end of the file is the triangle data area, and each line is a triangle; The format of the triangle is as follows: Triangle number, vertex 1 index, vertex 2 index, vertex 3 index; The vertices of the triangle are in counter-clockwise direction.

4. The three-dimensional geological model storage method based on geological contour lines according to claim 3, It is characterized in that In step S3, the process of extracting the boundary includes: start → read the input file → read the number of vertices and triangles → parse the vertex data → parse the triangle data → build edges from triangles → record the number of adjacent triangles for each edge → determine whether the number of adjacent triangles is 1 → Yes, boundary edge → output boundary and vertex coordinates → end; No, internal edge → end.

5. The three-dimensional geological model storage method based on geological contour lines according to claim 3, It is characterized in that In step S4, S is sorted from small to large in the elevation direction, and S is a set of triangles.

6. The three-dimensional geological model storage method based on geological contour lines according to claim 5, It is characterized in that Step S4 includes the following steps: S41. Sort S; Sort S in ascending order according to elevation. Since a triangle has three vertices and each vertex has a Z value, the Z value space occupied by the triangle is a range, namely: TR i =[From min ,WITH max ] in: WITH min =Min(Z i1 ,WITH i2 ,WITH i3 ) WITH max =Max(Z i1 ,WITH i2 ,WITH i3 ) Min is the minimum function, Max is the maximum function, Z i1 , Z i2 , Z i3 is the Z value of the three vertices of the i-th triangle; by calculating all triangles in S, we get the TR set. First, we calculate the Z value of the three vertices of the i-th triangle. min Sort from small to large. If there are more than two triangles, the Z min The same, then the Z of these triangles max Sort in ascending order, and finally get the sorted set TSeq; S42, calculate the contour line starting surface and set the contour interval; the contour line starting surface is based on the Z value of the first TSeq min Calculate the starting point of the contour line according to the contour drawing specification. b : in, The symbol indicates rounding down, and h is the equal height interval; S43. Calculate the plane equation Z of each contour line based on the starting surface and contour distance. i ; Z i =i*h+Z b Among them, i represents the contour line number, and the serial number starts from 0; S44, interpolate geological contours; extract triangles T from TSeq in sequence j , j is T j The sequence number in TSeq is based on T j Query the contour plane equation within the range and get the set Z sub , Z sub Each equation in the triangle T j By finding the intersection, we finally get a series of points associated with each contour plane, and connect these points in sequence to get several contour lines.

7. The three-dimensional geological model storage method based on geological contour lines according to any one of claims 1 to 6, It is characterized in that In step S1, a three-dimensional geological model is constructed using GoCad, Petrel or EarthVision modeling software.

8. A three-dimensional geological model storage device based on geological contour lines, It is characterized in that include: Modeling module, used to build three-dimensional geological body models; The first acquisition module is used to extract the upper surface data of each stratum in the three-dimensional geological body model; The second acquisition module is used to extract the upper surface boundary data of a certain stratum; Interpolation module, used to interpolate geological contour data in upper surface boundary data; The third acquisition module is used to intersect the geological contour lines with the upper surface boundary lines to obtain local layer data.

9. An electronic device, It is characterized in that include: a memory having a computer program stored thereon; A processor is used to execute the program in the memory to implement the three-dimensional geological model storage method based on geological contour lines as described in any one of claims 1 to 7.

10. A computer-readable storage medium, It is characterized in that A computer program is stored thereon, characterized in that when the program is executed by a processor, the three-dimensional geological model storage method based on geological contour lines as described in any one of claims 1 to 7 is implemented.

Images