Geological occurrence difference algorithm and device based on Delaunay triangulation

Abstract

The invention relates to the technical field of geology, in particular to a geologic occurrence difference algorithm based on Delaunay triangulation, which comprises the following steps of: 1, establishing a normalized geologic occurrence table which comprises a number, N, E, an occurrence direction D and an occurrence angle P; 2, constructing a Delaunay triangle according to N and E in the normalized geological occurrence table, establishing a Delaunay triangle table, and storing in a manner of triangle number, vertex 1, vertex 2 and vertex 3; 3, extracting edges of a Delaunay triangle based on the Delaunay triangle table, and establishing an edge information table which comprises an edge number, a vertex 1, a vertex 2 and a subordinate triangle; 4, a plurality of points P1, P2,..., Pn are obtained through intersection with the geological boundary line DZ1, the points are stored in an intersection point information table, and the intersection point information table comprises intersection point numbers, geological boundary numbers, edge numbers, N and E; and 5, solving the direction and angle of the occurrence of the intersection point, and interpolating the occurrence of any point of the geological boundary line.

Description

Technical Field

[0001] The invention relates to the field of geological technology, in particular to a geological occurrence difference algorithm and device based on Delaunay triangulation. Background Art

[0002] Three-dimensional geological modeling is a powerful tool for quantitatively studying underground geological information. It is widely used to display and analyze underground geological structures. High-precision three-dimensional geological models are not only an important guarantee for geological applications such as spatial analysis, numerical simulation, and resource calculation, but also an important means to identify and deepen the understanding of underground geological phenomena.

[0003] Geological occurrence data is an important basis for 3D geological modeling. Through accurate occurrence measurement and analysis, the accuracy and reliability of geological models can be effectively improved, providing a scientific basis for geological research and practical applications. Geological occurrence can determine the structure, layer and structural surface of the stratum, thereby helping geologists establish stratum boundaries and construct 3D geological models. The main methods of using occurrence modeling include:

[0004] 1. Spline numerical interpolation is an advanced method commonly used to construct three-dimensional geological models. The basic idea is to interpolate geological data through spline curves (such as cubic splines) to generate smooth and continuous geological layers and structures. The spline interpolation method based on geological occurrence specifically uses the geological occurrence as the control constraint value of the spline curve. When generating the spline curve, these occurrence data are used to control the strike and dip of the stratum. This method can more accurately reflect the actual spatial distribution and structural morphology of the geological body and provide a high-precision three-dimensional geological model.

[0005] 2. Kriging interpolation is a statistical interpolation method commonly used in geological modeling. It is based on the theory of regionalized variables and uses spatial autocorrelation to estimate the values ​​of unknown points to generate a smooth and accurate geological model. The Kriging interpolation method requires data preparation, preprocessing and analysis to find the distance and difference between data points, analyze the spatial autocorrelation of the data, determine the variability and related distance of the data, and construct a Kriging equation system. The weight of each known point to the interpolation point is calculated using the variogram model, and finally the Kriging equation system is solved to obtain the estimated value and variance of the interpolation point. In the Kriging interpolation method, the geological occurrence is used as a known point in the variogram model and is assigned a corresponding weight. In order to obtain a more accurate model, drilling data, geological profile data, etc. can also be used as known data to modify the model.

[0006] 3. Numerical simulation is a method of constructing a three-dimensional geological model by simulating geological processes and phenomena by computer. This method can simulate complex geological structures and evolution processes and generate high-precision geological models. The numerical simulation method uses the physical parameters of each geological unit, such as density, elastic modulus, permeability, etc., to establish the boundary conditions of the geological model, such as stress boundary, fluid boundary, etc. The stress field, deformation field and heat conduction of the geological body are simulated by the finite element analysis method, and the geological model is discretized into finite elements or finite difference grids for numerical calculation. The model state is gradually updated by iteratively solving the simulation equations.

[0007] However, the existing geological occurrence difference algorithm is difficult to generate the occurrence of any point on the geological boundary line. Summary of the invention

[0008] The purpose of the present invention is to provide a geological occurrence difference algorithm and device based on Delaunay triangulation to solve the problem that the existing geological occurrence difference algorithm in the prior art is difficult to generate the occurrence of any point on the geological boundary line.

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

[0010] In a first aspect, the present invention provides a geological occurrence difference algorithm based on Delaunay triangulation, comprising the following steps:

[0011] S1. Establish a normalized geological occurrence table, which includes number, N, E and rock formation occurrence, where number is the occurrence number and is unique, N represents the north coordinate in the engineering independent coordinate system, E represents the east coordinate in the engineering independent coordinate system, and rock formation occurrence describes the occurrence information of the occurrence point, including the direction D and the angle P of the occurrence;

[0012] S2. Construct a Delaunay triangle according to N and E in the normalized geological occurrence table, and establish a Delaunay triangle table, which includes a triangle number, vertex 1, vertex 2, and vertex 3. The triangle number is unique, each row represents a triangle, and vertex 1, vertex 2, and vertex 3 are arranged counterclockwise. The columns of vertex 1, vertex 2, and vertex 3 store the numbers in the digital geological occurrence table;

[0013] S3, extracting the edges of the Delaunay triangle based on the Delaunay triangle table, and establishing an edge information table, which includes the edge number, vertex 1, vertex 2, and the belonging triangle;

[0014] S4. The edge set in the edge information table intersects with the geological boundary line DZ1 to obtain several points P1, P2, ..., P n, store these points in the intersection information table, which contains the intersection number, geological boundary number, edge number, N and E, and each row represents an intersection;

[0015] S5. Find the direction of the occurrence of the intersection point;

[0016] S6. Find the angle of the occurrence of the intersection point;

[0017] S7. Interpolate the occurrence of any point on the geological boundary line.

[0018] Among them, Delaunay triangulation is a technique widely used in computational geometry and computer graphics. Its purpose is to split a set of points into several non-overlapping triangles to meet certain optimization characteristics. For a given set of plane points, Delaunay triangulation constructs a triangle network so that no point is within the circumscribed circle of any triangle. This triangulation has the following good properties:

[0019] Maximize the minimum angle: Delaunay triangulation tends to generate more uniform triangles and avoid very sharp triangles. It maximizes the minimum internal angle of all triangles.

[0020] Empty circle property: In Delaunay triangulation, for any triangle, its circumscribed circle does not contain any other points.

[0021] Uniqueness: For a set of points with no collinear points, the Delaunay triangulation is unique. However, if there are multiple collinear points, the triangulation may not be unique.

[0022] As a preferred technical solution of the present invention, in step S1, a geological occurrence information table is established based on the information of the occurrence points, and then normalized.

[0023] As a further preferred technical solution of the present invention, for D, the north direction N is defined as 0°, and it is digitized in the counterclockwise direction, and its value range is [0°, 360°). P takes the vector line D as the reference, the horizontal plane on the right hand of D as the reference plane 0°, and the horizontal plane on the left hand of D is 180°, and its value range is [0°, 180°].

[0024] As a preferred technical solution of the present invention, in step S2, a Delaunay triangle is constructed by using GIS software.

[0025] As a preferred technical solution of the present invention, in step S3, for the common edge of two triangles, one piece of edge information is stored in the edge information table, and the number of subordinate triangles is 2.

[0026] As a preferred technical solution of the present invention, in step S5, assuming that the geological boundary line DZ1 and the triangle side E1 intersect at point P1, and the two vertices of the side E1 are V1 and V2, the calculation of the orientation direction of P1 includes the following steps:

[0027] S51, calculate the weights of V1 and V2,

[0028]

[0029] Where ‖P1V2‖ refers to the distance between point P1 and point V2, and the other terms in the formula are the same;

[0030] S52. Use the following formula to process the directions of V1 and V2.

[0031] f(x)=sin(x / 2)

[0032] but:

[0033]

[0034] Among them, D V1 is the occurrence direction of V1, D V2 is the occurrence direction of V2,

[0035] The interpolated data is:

[0036] D p1 =ω1μ1+ω2μ2;

[0037] S53, calculate the direction D of P1 occurrence based on Dp1,

[0038] w=2*arcsin(D p1 )

[0039] but:

[0040]

[0041] As a further preferred technical solution of the present invention, in step S6, the angle P of the occurrence of P1 is calculated using the following formula:

[0042] P=ω1P v1 +ω2P v2

[0043] Among them, P V1 is the occurrence direction of V1, P V2 is the occurrence direction of V2.

[0044] As a further preferred technical solution of the present invention, in step S7, for the edge intersecting with the modeling boundary B, since there is no generation point at one end, P1P2 is interpolated on the internal line segment of the Delaunay triangle, P1 and P2 are used as known generation points, and then steps S5 and S6 are used to interpolate the internal points.

[0045] In a second aspect, the present invention further provides a geological occurrence difference device based on Delaunay triangulation, comprising:

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

[0047] A processor is used to execute the program in the memory to implement the geological occurrence difference algorithm based on Delaunay triangulation as described in any one of the above items.

[0048] In a third aspect, the present invention further provides a computer-readable storage medium having a computer program stored thereon, characterized in that when the program is executed by a processor, the geological occurrence difference algorithm based on Delaunay triangulation as described in any one of the above items is implemented.

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

[0050] The present invention discloses a geological occurrence difference algorithm based on Delaunay triangulation, which uses geological occurrence points to construct Delaunay triangles, intersects geological boundary lines with Delaunay triangles to obtain intersection points of boundary lines, and the occurrence of the intersection points is calculated by inverse distance weighted interpolation. Each intersection point generates two values, one is the direction of the occurrence point, and the other is the angle of the occurrence point. The occurrences of other points on the boundary line are also calculated by inverse distance weighted interpolation. Finally, the occurrence can be obtained at any point on the boundary line. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] Figure 1 It is a schematic diagram of the D direction;

[0052] Figure 2 It is a schematic diagram of the angle of occurrence;

[0053] Figure 3 A schematic diagram of constructing a Delaunay triangle based on coordinates;

[0054] Figure 4 This is a schematic diagram of the intersection between the geological boundary line and the Delaunay edge;

[0055] Figure 5 A schematic diagram for solving the occurrence of the intersection point;

[0056] Figure 6Schematic diagram of the interpolated boundary line occurrence point. DETAILED DESCRIPTION

[0057] 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.

[0058] 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.

[0059] 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.

[0060] 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.

[0061] 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.

[0062] 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.

[0063] In the related art, the existing geological occurrence difference algorithm is difficult to generate the occurrence of any point on the geological boundary line. Figures 1 to 6 To elaborate.

[0064] Example 1

[0065] like Figures 1 to 6 As shown, the geological occurrence difference algorithm based on Delaunay triangulation described in the present invention includes the following steps:

[0066] Step 1: Data preparation.

[0067] Fill in the information of the occurrence point as shown in Table 1.

[0068] Table 1. Geological occurrence information table

[0069] serial number N (North coordinate) E(East coordinate) Rock formation Rock Formation Description <![CDATA[V1]]> 3287416.42 499949.90 N45E / 40SE <![CDATA[V2]]> 3287476.87 499800.90 S13W / 85SE …… …… …… …… ……

[0070] The number is the occurrence number, which must be unique. N (North coordinate) represents the north coordinate in the engineering independent coordinate system, and E (East coordinate) represents the east coordinate in the engineering independent coordinate system. The rock formation occurrence describes the occurrence information of the occurrence point. The direction D is before the “ / ” sign, and the angle P is after the “ / ” sign. The rock formation description is a text description of the occurrence point.

[0071] Step 2: Normalize the occurrence information.

[0072] 1. Direction of occurrence D.

[0073] D has four characters, with letters on both sides and a character in the middle. For example, N60E means 60° east of north. The first letter of D can only be N or S, indicating the north or south direction, and the last character must be E or W, indicating the east or west direction. The number in the middle indicates how many degrees west (or east) the north (or south) is. Figure 1 shown.

[0074] 2. Angle P of occurrence.

[0075] P has four characters. The first two characters are numbers, indicating the angle of the occurrence cut, and the last two characters are characters, indicating the direction of the cut. For example, 40SE means that the occurrence cuts 40°, which is perpendicular to direction D and faces southeast. Figure 2 As shown, the angle of occurrence is the angle between plane P and the horizontal plane, and the direction is southeast.

[0076] 3. Digitalization of occurrence status.

[0077] In order to facilitate the interpolation of D and P, D and P need to be digitized. For D, the north direction N is defined as 0° and digitized in the counterclockwise direction. Its value range is [0°, 360°). P takes the vector line D as the reference, and the horizontal plane on the right hand of D is the reference plane (0°), to the horizontal plane on the left hand of D (180°), and its value range is [0°, 180°].

[0078] The digitized occurrence information is shown in Table 2.

[0079] Table 2. Normalized geological occurrence table

[0080] serial number N (North coordinate) E(East coordinate) D(Direction) P(angle) Rock Formation Description <![CDATA[V1]]> 3287416.42 499949.90 45 40 <![CDATA[V2]]> 3287476.87 499800.90 193 95 …… …… …… …… ……

[0081] Step 3: Construct a Delaunay triangle.

[0082] According to N and E in Table 2, Delaunay triangles can be constructed using professional GIS software (such as arcgis), or a program can be written to construct Delaunay triangles, such as Figure 3 The constructed triangles are stored in Table 3.

[0083] Table 3. Delaunay triangle table

[0084] Triangle number Vertex 1 Vertex 2 Vertex 3 describe <![CDATA[T1]]> <![CDATA[V1]]> <![CDATA[V2]]> <![CDATA[V3]]> <![CDATA[T2]]> <![CDATA[V2]]> <![CDATA[V3]]> <![CDATA[V4]]> …… …… …… …… ……

[0085] In Table 3, the triangle numbers are unique, each row represents a triangle, and vertices 1, 2, and 3 are arranged counterclockwise. Vertex 1, 2, and 3 columns store the numbers in Table 2. By querying Table 2, complete information of each occurrence point can be obtained.

[0086] Step 4: Construct side information.

[0087] In the process of finding the intersection between the Delaunay triangle and the geological boundary, it is necessary to extract the edges of the Delaunay triangle. The edges of the Delaunay triangle are extracted according to Table 3 and stored in Table 4.

[0088] Table 4. Edge information table

[0089] Edge Number Vertex 1 Vertex 2 Belongs to the triangle describe <![CDATA[E1]]> <![CDATA[V1]]> <![CDATA[V2]]> <![CDATA[T1]]> <![CDATA[E2]]> <![CDATA[V2]]> <![CDATA[V3]]> <![CDATA[T1,T2]]> …… …… …… …… ……

[0090] In particular, for the common edge of two triangles, only one piece of edge information needs to be stored in Table 4, and it can belong to two triangles.

[0091] Step 5: Find the intersection of the geological boundary line and the Delaunay triangle edge.

[0092] like Figure 4 As shown, take a geological boundary line DZ1, and intersect DZ1 with the edge set in Table 4 to obtain several points P1, P2, ..., P n , and store these points in Table 5.

[0093] Table 5. Intersection information table

[0094] Intersection number Geological Boundary Number Edge Number N (North coordinate) E(East coordinate) describe <![CDATA[P1]]> <![CDATA[DZ1]]> <![CDATA[E1]]> 3287412.42 499943.91 <![CDATA[P2]]> <![CDATA[DZ1]]> <![CDATA[E2]]> 3287473.87 499801.91 …… …… …… …… ……

[0095] Each row in Table 5 represents an intersection point, and N and E are the coordinates of the intersection point.

[0096] Step 6: Find the direction of the intersection.

[0097] like Figure 5 As shown, assuming that DZ1 and E1 intersect at point P1, and the two vertices of edge E1 are V1 and V2, the direction of the occurrence of P1 can be calculated by the following steps:

[0098] First, calculate the weights of V1 and V2,

[0099]

[0100] Among them, ||P1V2|| refers to the distance between point P1 and point V2, and the other terms in the formula are the same. Secondly, the direction of the occurrence of V1 and V2 is processed using the following formula:

[0101] f(x)=sin(x / 2)

[0102] but:

[0103]

[0104] Among them, D v1 is the occurrence direction of V1, D v2 is the strike direction of V2, then the interpolation data is:

[0105] D p1 =ω1μ1+ω2μ2.

[0106] Finally, the direction D of P1 occurrence is calculated based on Dp1.

[0107] w=2*arcsin(D p1 )

[0108] but:

[0109]

[0110] Step 7: Find the angle of the intersection.

[0111] The angle P of P1 is calculated using the following formula.

[0112] P=ω1P v1 +ω2P v2

[0113] Among them, P v1 is the occurrence angle of V1, P v2 is the occurrence angle of V2.

[0114] Step 8: Interpolate the occurrence of geological boundaries.

[0115] from Figure 6 It can be seen that DZ1 is divided into several parts by the Delaunay triangle and the modeling boundary B, and the intersection with B is B1. When interpolating DZ1, it is necessary to divide DZ1 into several parts and interpolate them separately. For the interpolation of the internal line segments of the Delaunay triangle (such as: P1P2), P1 and P2 can be used as known occurrence points, and then steps six and seven are used to interpolate the internal points.

[0116] For the edge that intersects the boundary, since there is no occurrence point at one end, its occurrence needs to be inferred through P1P2:

[0117] Calculate the weights of P1 and P2,

[0118]

[0119] Among them, X1 is a point in the line segment B1P1, and then steps 6 and 7 are used to interpolate the internal points.

[0120] The geological occurrence difference algorithm based on Delaunay triangulation described in this embodiment uses geological occurrence points to construct Delaunay triangles, and intersects the geological boundary line with the Delaunay triangle to obtain the intersection of the boundary line. The occurrence of the intersection point is calculated by the inverse distance weighted interpolation method. Each intersection generates two values, one is the direction of the occurrence point, and the other is the angle of the occurrence point. The occurrence of other points on the boundary line is also calculated by the inverse distance weighted interpolation method. Finally, the occurrence can be obtained at any point on the boundary line.

[0121] Example 2

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

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

[0124] A processor is used to execute the program in the memory to implement the geological occurrence difference algorithm based on Delaunay triangulation as described in Example 1.

[0125] 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.

[0126] 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 geological occurrence difference algorithm based on Delaunay triangulation.

[0127] 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.

[0128] 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.

[0129] 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.

[0130] 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.

[0131] 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 geological occurrence difference algorithm based on Delaunay triangulation.

[0132] 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 geological occurrence difference algorithm based on Delaunay triangulation.

[0133] Example 3

[0134] The computer-readable storage medium described in the present invention stores a computer program thereon, characterized in that when the program is executed by a processor, the geological occurrence difference algorithm based on Delaunay triangulation as described in Example 1 is implemented.

[0135] 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.

[0136] 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 geological occurrence difference algorithm based on Delaunay triangulation, characterized in that: The following steps are involved: S1. Establish a normalized geological occurrence table, which contains number, N, E and rock formation occurrence, where number is the occurrence number and is unique, N represents the north coordinate in the engineering independent coordinate system, E represents the east coordinate in the engineering independent coordinate system, and rock formation occurrence describes the occurrence information of the occurrence point, including the direction D and the angle P of the occurrence; S2. Construct a Delaunay triangle according to N and E in the normalized geological occurrence table, and establish a Delaunay triangle table, which contains triangle numbers, vertex 1, vertex 2, and vertex 3. The triangle numbers are unique, and each row represents a triangle. Vertex 1, vertex 2, and vertex 3 are arranged counterclockwise, and the columns of vertex 1, vertex 2, and vertex 3 store the numbers in the digital geological occurrence table; S3, extracting the edges of the Delaunay triangle based on the Delaunay triangle table, and establishing an edge information table, which includes the edge number, vertex 1, vertex 2, and the belonging triangle; S4. The edge set in the edge information table intersects with the geological boundary line DZ1 to obtain several points P1, P2, ..., P n , store these points in the intersection information table, which contains the intersection number, geological boundary number, edge number, N and E, and each row represents an intersection; S5. Find the direction of the occurrence of the intersection point; S6. Find the angle of the occurrence of the intersection point; S7. Interpolate the occurrence of any point on the geological boundary line.

2. The geological occurrence difference algorithm based on Delaunay triangulation according to claim 1 is characterized in that: In step S1, a geological occurrence information table is established based on the information of the occurrence points, and then normalized.

3. The geological occurrence difference algorithm based on Delaunay triangulation according to claim 2 is characterized in that: For D, the north direction N is defined as 0°, and it is digitized in the counterclockwise direction, and its value range is [0°, 360°). P takes the vector line D as the reference, the horizontal plane on the right hand of D as the reference plane 0°, and the horizontal plane on the left hand of D is 180°, and its value range is [0°, 180°].

4. The geological occurrence difference algorithm based on Delaunay triangulation according to claim 1 is characterized in that: In step S2, a Delaunay triangle is constructed using GIS software.

5. The geological occurrence difference algorithm based on Delaunay triangulation according to claim 1 is characterized in that: In step S3, for the common edge of two triangles, one piece of edge information is stored in the edge information table, and the number of the belonging triangles is 2.

6. The geological occurrence difference algorithm based on Delaunay triangulation according to any one of claims 1 to 5, characterized in that: In step S5, assuming that the geological boundary line DZ1 and the triangle side E1 intersect at point P1, and the two vertices of the side E1 are V1 and V2, the calculation of the orientation direction of P1 includes the following steps: S51, calculate the weights of V1 and V2, Where ‖P1V2‖ refers to the distance between point P1 and point V2, and the other terms in the formula are the same; S52. Use the following formula to process the directions of V1 and V2. f(x)=sin(x / 2) but: Among them, D V1 is the occurrence direction of V1, D V2 is the occurrence direction of V2, The interpolated data is: D p1 =ω1μ1+ω2μ2; S53, according to D p1 Calculate the direction D of P1 occurrence, w=2*arcsin(D p1 )but:

7. The geological occurrence difference algorithm based on Delaunay triangulation according to claim 6 is characterized in that: In step S6, the angle P of the P1 occurrence is calculated using the following formula: P=ω1P v1 +ω2P v2 Among them, P V1 is the occurrence direction of V1, P V2 is the occurrence direction of V2.

8. The geological occurrence difference algorithm based on Delaunay triangulation according to claim 7 is characterized in that: In step S7, for the edge intersecting with the modeling boundary B, since one end has no origin point, the internal line segment of the Delaunay triangle is interpolated P1P2, P1 and P2 are used as known origin points, and then steps S5 and S6 are used to interpolate the internal points.

9. A geological occurrence difference device based on Delaunay triangulation, 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 geological occurrence difference algorithm based on Delaunay triangulation as described in any one of claims 1 to 8.

10. A computer-readable storage medium, characterized in that: A computer program is stored thereon, characterized in that when the program is executed by a processor, the geological occurrence difference algorithm based on Delaunay triangulation as described in any one of claims 1 to 8 is implemented.

Images