CN116206079B - Moving tetrahedron-based geologic body modeling method and related equipment - Google Patents
Moving tetrahedron-based geologic body modeling method and related equipment Download PDFInfo
- Publication number
- CN116206079B CN116206079B CN202310499952.1A CN202310499952A CN116206079B CN 116206079 B CN116206079 B CN 116206079B CN 202310499952 A CN202310499952 A CN 202310499952A CN 116206079 B CN116206079 B CN 116206079B
- Authority
- CN
- China
- Prior art keywords
- tetrahedron
- implicit
- points
- elements
- target area
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- Image Generation (AREA)
Abstract
The invention provides a moving tetrahedron-based geologic body modeling method and related equipment, wherein the moving tetrahedron-based geologic body modeling method comprises the following steps: acquiring drilling exploration data of a target area to obtain boundary point information data of the surface of a geologic body; constructing an implicit curved surface function model based on the boundary point information data, and solving the implicit curved surface function model to obtain a geological implicit curved surface function model; constructing a plurality of cube elements according to geological data of a target area and preset grid resolution, and dividing each cube element to obtain a plurality of tetrahedron elements; substituting coordinate values of four points on the tetrahedron element into a geological implicit surface function model for calculation according to each tetrahedron element to obtain implicit surface function values corresponding to the four points on all the tetrahedron elements; cutting all tetrahedral elements according to the function values of each implicit curved surface to obtain a geological three-dimensional entity model of the target area; the modeling precision and the matching degree of the geologic body are improved.
Description
Technical Field
The invention relates to the technical field of three-dimensional geological modeling, in particular to a geological modeling method based on a mobile tetrahedron and related equipment.
Background
At present, three-dimensional geologic modeling has become a popular research object in the field of geologic modeling, and two-dimensional models cannot conveniently and clearly show geologic conditions, but compared with the traditional information display method, the three-dimensional geologic model can more clearly and intuitively represent various environmental conditions of geology. The three-dimensional geologic model can characterize the spatial characteristics of the geologic body from different angles, and is helpful for researchers to analyze and obtain more geologic information.
Implicit surface functions can be used to characterize geological surfaces, and there are two main methods for visualizing implicit surface functions: surface painting and volume painting methods. The mobile tetrahedron algorithm belongs to a first visualization method, namely a surface drawing method, and the surface of the geologic body is characterized by extracting an isosurface in an implicit surface function. In the case of visualizing an implicit surface function model, common basic units include cubes, hexahedrons, tetrahedrons, and the like, and the basic units are also called voxels or primitive elements for representing a filling space. Each surface in hexahedron and tetrahedron usually keeps parallel state with coordinate axis, and the geological model that builds under this condition, the boundary is comparatively coarse, and is not smooth enough, can't be with the accurate characterization of geological boundary.
Disclosure of Invention
The invention provides a moving tetrahedron-based geologic body modeling method and related equipment, and aims to improve accuracy and matching degree of geologic body modeling.
In order to achieve the above object, the present invention provides a moving tetrahedron-based geologic body modeling method, comprising:
step 2, constructing an implicit surface function model based on the boundary point information data, and solving the implicit surface function model to obtain a geological implicit surface function model;
step 3, constructing a plurality of cube elements according to geological data of a target area and preset grid resolution, and dividing each cube element to obtain a plurality of tetrahedron elements;
step 4, substituting coordinate values of four points on the tetrahedron element into a geological implicit surface function model for calculation according to each tetrahedron element to obtain implicit surface function values corresponding to the four points on all the tetrahedron elements;
and 5, cutting all tetrahedral elements according to the function values of each implicit curved surface to obtain a three-dimensional solid model of the geologic body of the target area.
Further, step 1 includes:
acquiring borehole exploration data of a target area;
carrying out combined sample analysis on the drilling exploration data to obtain lithology samples of the target area;
and processing adjacent different lithology samples in the target area to obtain boundary point information data of the surface of the geologic body.
Further, before step 5, the method further comprises:
according to the relation between the implicit surface function values corresponding to the four points on all tetrahedron elements and the preset grid size, solving the four points on all tetrahedron elements by a central difference method to obtain gradient values corresponding to the points;
and taking the gradient value corresponding to each point position as the normal vector of four point positions on all tetrahedron elements.
Further, the expression of the center difference method is:
solving four points on all tetrahedron elements by a central difference method to obtain gradient values corresponding to the points as follows:
wherein the gradient valueNormal vector as four points on all tetrahedrons +.>,/>For the grid serial number> The length of the grid in the three directions X, Y, Z, respectively,>is a solving expression of the implicit surface function value.
Further, before step 5, the method further comprises:
according to the implicit surface function values corresponding to the four points on all tetrahedron elements, respectively extracting equivalent points on all sides of each tetrahedron element to obtain an equivalent point set;
and according to the equivalence point set, triangular surface patch drawing is sequentially carried out in each tetrahedron element, and an equivalence surface corresponding to the geologic body of the target area is obtained.
Further, step 5 includes:
defining tetrahedral elements with constraint values of four points being greater than zero as absolute internal tetrahedrons,
defining tetrahedral elements with constraint values smaller than or equal to zero in at least one of the four points as boundary tetrahedrons;
taking the isosurface as an interface, and carrying out tetrahedron cutting on the boundary tetrahedron to obtain a part of the boundary tetrahedron positioned in the implicit curved surface;
the part of the absolute internal tetrahedron and the boundary tetrahedron positioned in the implicit curved surface forms a geological body three-dimensional solid model of the target area.
Further, after step 5, the method further comprises:
and calculating the three-dimensional solid model of the geologic body by using an Euler tetrahedral formula to obtain the volume of the geologic body.
The invention also provides a geological body modeling device based on the mobile tetrahedron, which comprises:
the acquisition module is used for acquiring drilling exploration data of the target area to obtain boundary point information data of the surface of the geologic body;
the construction module is used for constructing an implicit curved surface function model based on the boundary point information data, and solving the implicit curved surface function model to obtain a geological implicit curved surface function model;
the dividing module is used for constructing a plurality of cube elements according to geological data of the target area and preset grid resolution, and dividing each cube element to obtain a plurality of tetrahedron elements;
the calculation module is used for substituting coordinate values of four points on the tetrahedron element into the geological implicit surface function model for calculation according to each tetrahedron element to obtain implicit surface function values corresponding to the four points on all the tetrahedron elements;
and the cutting module is used for cutting all tetrahedron elements according to the function values of each implicit curved surface to obtain a geological three-dimensional entity model of the target area.
The invention also provides a computer readable storage medium, wherein the computer readable storage medium stores a computer program, and the computer program realizes a moving tetrahedron-based geologic body modeling method when being executed by a processor.
The invention also provides a terminal device, which comprises a memory, a processor and a computer program stored in the memory and capable of running on the processor, wherein the processor realizes the modeling method of the geologic body based on the mobile tetrahedron when executing the computer program.
The scheme of the invention has the following beneficial effects:
according to the method, boundary point information data of the surface of the geologic body is obtained by obtaining drilling exploration data of a target area; constructing an implicit curved surface function model based on the boundary point information data, and solving the implicit curved surface function model to obtain a geological implicit curved surface function model; constructing a plurality of cube elements according to geological data of a target area and preset grid resolution, and dividing each cube element to obtain a plurality of tetrahedron elements; substituting coordinate values of four points on the tetrahedron element into a geological implicit surface function model for calculation according to each tetrahedron element to obtain implicit surface function values corresponding to the four points on all the tetrahedron elements; cutting all tetrahedral elements according to the function values of each implicit curved surface to obtain a geological three-dimensional entity model of the target area; the method solves the problems of lower accuracy of the geologic body and lower matching degree with the geologic body in the traditional method, establishes a geologic entity model which is more matched with the actual geologic body, and improves modeling accuracy and matching degree of the geologic body.
Other advantageous effects of the present invention will be described in detail in the detailed description section which follows.
Drawings
FIG. 1 is a schematic flow chart of an embodiment of the present invention;
fig. 2 (a) is a schematic diagram of a first contour point extraction case according to an embodiment of the present invention; fig. 2 (b) is a schematic diagram of a second contour point extraction case according to an embodiment of the present invention; fig. 2 (c) is a schematic diagram of a third contour point extraction case according to an embodiment of the present invention; fig. 2 (d) is a schematic diagram of a fourth contour point extraction case according to an embodiment of the present invention; fig. 2 (e) is a schematic diagram of a fifth contour point extraction case according to an embodiment of the present invention; fig. 2 (f) is a schematic diagram of a sixth contour point extraction case according to an embodiment of the present invention; fig. 2 (g) is a schematic diagram of a seventh contour point extraction case according to an embodiment of the present invention; fig. 2 (h) is a schematic diagram of an eighth contour point extraction case according to an embodiment of the present invention; fig. 2 (i) is a schematic diagram of a ninth contour point extraction case according to an embodiment of the present invention; fig. 2 (j) is a schematic diagram of a tenth contour point extraction case according to an embodiment of the present invention; fig. 2 (k) is a schematic diagram of an eleventh contour point extraction case according to an embodiment of the present invention; fig. 2 (l) is a schematic diagram of a twelfth contour point extraction case according to an embodiment of the present invention; FIG. 2 (m) is a schematic diagram of a thirteenth contour point extraction case according to an embodiment of the present invention; fig. 2 (n) is a schematic diagram of a fourteenth contour point extraction case according to an embodiment of the present invention; FIG. 2 (o) is a schematic diagram of a fifteenth contour point extraction case according to an embodiment of the present invention; fig. 2 (p) is a schematic diagram of a sixteenth contour point extraction case according to an embodiment of the present invention;
FIG. 3 (a) is a schematic diagram of a first specific cutting case based on a mobile tetrahedral algorithm according to an embodiment of the present invention; FIG. 3 (b) is a schematic diagram of a second specific cutting case based on a mobile tetrahedral algorithm according to an embodiment of the present invention; FIG. 3 (c) is a schematic diagram of a third specific cutting case based on a mobile tetrahedral algorithm according to an embodiment of the present invention; FIG. 3 (d) is a schematic diagram of a fourth specific cutting case based on a mobile tetrahedral algorithm according to an embodiment of the present invention; FIG. 3 (e) is a schematic diagram of a fifth specific cutting case based on a mobile tetrahedral algorithm according to the embodiment of the present invention; FIG. 3 (f) is a schematic diagram of a sixth specific cutting case based on a mobile tetrahedral algorithm according to the embodiment of the present invention; FIG. 3 (g) is a schematic diagram of a seventh specific cutting case based on a mobile tetrahedral algorithm according to an embodiment of the present invention; fig. 3 (h) is a schematic diagram of an eighth specific cutting case based on a mobile tetrahedral algorithm according to an embodiment of the present invention; FIG. 3 (i) is a schematic diagram of a ninth specific cutting case based on a mobile tetrahedral algorithm according to an embodiment of the present invention; fig. 3 (j) is a schematic diagram of a tenth specific cutting case based on a mobile tetrahedral algorithm according to an embodiment of the present invention; FIG. 3 (k) is a schematic diagram of an eleventh specific cutting case based on a mobile tetrahedral algorithm according to an embodiment of the present invention; FIG. 3 (l) is a schematic diagram of a twelfth specific cutting case based on a mobile tetrahedral algorithm according to an embodiment of the present invention; FIG. 3 (m) is a schematic diagram of a thirteenth specific cutting case based on a mobile tetrahedral algorithm according to the embodiment of the present invention; FIG. 3 (n) is a schematic diagram of a fourteenth specific cutting case based on a mobile tetrahedral algorithm according to an embodiment of the present invention; FIG. 3 (o) is a schematic diagram of a fifteenth specific cutting case based on a mobile tetrahedral algorithm according to an embodiment of the present invention; FIG. 3 (p) is a schematic diagram of a sixteenth specific cutting case based on a mobile tetrahedral algorithm according to an embodiment of the present invention;
fig. 4 is a cross-sectional view of a three-dimensional geologic body solid model constructed in accordance with an embodiment of the invention.
Detailed Description
In order to make the technical problems, technical solutions and advantages to be solved more apparent, the following detailed description will be given with reference to the accompanying drawings and specific embodiments. It will be apparent that the described embodiments are some, but not all, embodiments of the invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
In the description of the present invention, it should be noted that the directions or positional relationships indicated by the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", etc. are based on the directions or positional relationships shown in the drawings, are merely for convenience of describing the present invention and simplifying the description, and do not indicate or imply that the devices or elements referred to must have a specific orientation, be configured and operated in a specific orientation, and thus should not be construed as limiting the present invention. Furthermore, the terms "first," "second," and "third" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance.
In the description of the present invention, it should be noted that, unless explicitly stated and limited otherwise, the terms "mounted," "connected," and "connected" are to be construed broadly, and may be, for example, a locked connection, a removable connection, or an integral connection; can be mechanically or electrically connected; can be directly connected or indirectly connected through an intermediate medium, and can be communication between two elements. The specific meaning of the above terms in the present invention will be understood in specific cases by those of ordinary skill in the art.
In addition, the technical features of the different embodiments of the present invention described below may be combined with each other as long as they do not collide with each other.
The invention provides a geologic body modeling method based on a mobile tetrahedron and related equipment aiming at the existing problems.
As shown in fig. 1, an embodiment of the present invention provides a moving tetrahedron-based geologic body modeling method, including:
step 2, constructing an implicit surface function model based on the boundary point information data, and solving the implicit surface function model to obtain a geological implicit surface function model;
step 3, constructing a plurality of cube elements according to geological data of a target area and preset grid resolution, and dividing each cube element to obtain a plurality of tetrahedron elements;
step 4, substituting coordinate values of four points on the tetrahedron element into a geological implicit surface function model for calculation according to each tetrahedron element to obtain implicit surface function values corresponding to the four points on all the tetrahedron elements;
and 5, cutting all tetrahedral elements according to the function values of each implicit curved surface to obtain a three-dimensional solid model of the geologic body of the target area.
Specifically, step 1 includes:
acquiring borehole exploration data of a target area;
carrying out combined sample analysis on the drilling exploration data to obtain lithology samples of the target area;
and processing adjacent different lithology samples in the target area to obtain boundary point information data of the surface of the geologic body.
In the embodiment of the invention, when drilling and exploration are carried out on a target area, data in the exploration process, such as drilling and exploration data of azimuth angle, dip angle, layer bottom elevation, layer bottom depth, layering thickness and the like of a drilling section, are required to be recorded; after the drilling exploration engineering is completed, the data are processed and arranged into corresponding column diagrams, and analysis can be performed according to the drilling column diagrams to obtain the specific position of the drilling, the inclinometry information of the drilling, the sampling condition and the like. Four relevant tables can be obtained by sorting the data: a borehole perforation table, a borehole inclinometry data table, a borehole lithology table, and a borehole sample analysis table. And carrying out combined sample analysis on the data to obtain a combined sample table of the stratum. Through the series of sorting, analyzing and cataloguing operations, lithology and other relevant information of a certain section of drilling can be obtained rapidly. The boundary point information data can be obtained by processing the adjacent different lithology sample segments, and the boundary point information data can be used as initial data.
Specifically, step 2, constructing an implicit surface function model based on boundary point information data, and solving the implicit surface function model to obtain a geological implicit surface function model, which comprises the following steps:
after preprocessing original drilling exploration data, information such as coordinates of a boundary point set and normal vectors on the surface of a geologic body can be obtained, the data information is used as initial input data for solving an implicit curved surface function, a proper implicit curved surface reconstruction method is selected to construct a geologic body implicit curved surface function model, and the implicit curved surface reconstruction method used in the embodiment of the invention is a radial basis implicit curved surface reconstruction method based on a first order polynomial, an implicit curved surface function equation is constructed, and the implicit curved surface function model of the surface of the geologic body is obtained by solving the implicit curved surface function equation.
Specifically, step 3, constructing a plurality of cube elements according to geological data of a target area and a preset grid resolution, and dividing each cube element to obtain a plurality of tetrahedron elements, including:
and analyzing geological data of the target area, constructing cube elements according to the set grid resolution, and then subdividing cubes one by one to divide each cube element into five tetrahedron elements.
Specifically, step 4, respectively substituting coordinate values of four points on the tetrahedron element into the geological implicit surface function model for calculation to obtain implicit surface function values corresponding to the four points on all the tetrahedron element, includes:
in the embodiment of the invention, coordinate values of four points on all tetrahedron elements are substituted into an implicit surface function model of the surface of the geological body respectively to obtain implicit surface function values corresponding to the four points on all tetrahedron elements; the implicit surface function value of each point represents the distance between the point and the implicit surface, the implicit surface function value is regular and represents that the point is positioned inside the implicit surface, the implicit surface function value is negative and represents that the point is positioned outside the implicit surface, and if the implicit surface function value is zero, the point is positioned on the implicit surface.
Specifically, before step 5, the method further comprises:
according to the relation between the implicit surface function values corresponding to the four points on all tetrahedron elements and the preset grid size, solving the four points on all tetrahedron elements by a central difference method to obtain gradient values corresponding to the points;
and taking the gradient value corresponding to each point position as the normal vector of four point positions on all tetrahedron elements.
In the embodiment of the invention, after the implicit surface function values corresponding to the four points on all tetrahedron elements are obtained, the gradient values of the points are solved by using a central difference method in combination with the set geobody grid size, and the gradient values corresponding to the points are used as normal vectors of the four points on all tetrahedron elements, so that the built model is more continuous and smooth.
Specifically, the expression of the center difference method is:
solving four points on all tetrahedron elements by a central difference method to obtain gradient values corresponding to the points as follows:
wherein the gradient valueNormal vector as four points on all tetrahedrons +.>,/>For the grid serial number> The length of the grid in the three directions X, Y, Z, respectively,>is a solving expression of the implicit surface function value.
Specifically, before step 5, the method further comprises:
according to the implicit surface function values corresponding to the four points on all tetrahedron elements, respectively extracting equivalent points on all sides of each tetrahedron element to obtain an equivalent point set;
according to the implicit surface function value of each point on the tetrahedron element, the intersection condition of the implicit surface model and the tetrahedron element can be obtained by analysis, in the embodiment of the invention, taking one tetrahedron element as an example, four point positions on the tetrahedron element are defined as p 1 、p 2 、p 3 、p 4 Sixteen conditions exist in relation between the implicit surface function model and tetrahedron elements, as shown in fig. 2, if the implicit function values of the four points in (a) are all greater than zero, the tetrahedron is located in the implicit surface, and no equivalent point exists; (b) Middle p 1 Is smaller than zero, p 2 、p 3 、p 4 If the implicit function value of (2) is greater than zero, extracting p respectively 1 p 2 、p 1 p 3 、p 1 p 4 The equivalent points c, a and b are used as triangle vertexes to form an equivalent surface; (c) Middle p 2 Is smaller than zero, p 1 、p 3 、p 4 If the implicit function value of (2) is greater than zero, extracting p respectively 1 p 2 、p 2 p 3 、p 3 p 4 The equivalent points c, a and b are used as triangle vertexes to form an equivalent surface; (d) Middle p 3 Is smaller than zero, p 1 、p 2 、p 4 If the implicit function value of (2) is greater than zero, extracting p respectively 1 p 3 、p 2 p 3 、p 3 p 4 The equivalent points a, b and c are used as triangle vertexes to form an equivalent surface; (e) Middle p 4 Is smaller than zero, p 1 、p 2 、p 3 If the implicit function value of (2) is greater than zero, extracting p respectively 1 p 4 、p 2 p 3 、p 3 p 4 The equivalent points c, b and a are used as triangle vertexes to form an equivalent surface; (f) Middle p 1 、p 2 Is smaller than zero, p 3 、p 4 If the implicit function value of (2) is greater than zero, extracting p respectively 1 p 3 、p 2 p 3 、p 2 p 4 、p 1 p 4 The upper equivalence point a, b, c, d takes a, b, c, d as a quadrilateral vertex to form an equivalence surface; (g) Middle p 1 、p 3 Is smaller than zero, p 2 、p 4 If the implicit function value of (2) is greater than zero, extracting p respectively 1 p 2 、p 2 p 3 、p 3 p 4 、p 1 p 4 The upper equivalence point a, b, c, d takes a, b, c, d as a quadrilateral vertex to form an equivalence surface; (h) Middle p 1 、p 4 Is smaller than zero, p 2 、p 3 If the implicit function value of (2) is greater than zero, extracting p respectively 1 p 2 、p 1 p 3 、p 3 p 4 、p 1 p 4 The upper equivalence point a, b, c, d takes a, b, c, d as a quadrilateral vertex to form an equivalence surface; (i) P is p 2 、p 3 Is smaller than zero, p 1 、p 4 If the implicit function value of (2) is greater than zero, extracting p respectively 1 p 3 、p 1 p 2 、p 2 p 4 、p 1 p 2 The upper equivalence point a, b, c, d takes a, b, c, d as a quadrilateral vertex to form an equivalence surface; (j) Middle p 2 、p 4 Is smaller than zero, p 1 、p 3 If the implicit function value of (2) is greater than zero, extracting p respectively 1 p 2 、p 2 p 3 、p 3 p 4 、p 1 p 4 The upper equivalence point a, b, c, d takes a, b, c, d as a quadrilateral vertex to form an equivalence surface; (k) Middle p 3 、p 4 Is smaller than zero, p 1 、p 2 If the implicit function value of (2) is greater than zero, extracting p respectively 1 p 3 、p 2 p 3 、p 2 p 4 、p 1 p 4 The upper equivalence point a, b, c, d takes a, b, c, d as a quadrilateral vertex to form an equivalence surface; (l) Middle p 1 、p 2 、p 3 Is smaller than zero, p 4 If the implicit function value of (2) is greater than zero, extracting p respectively 1 p 4 、p 2 p 4 、p 3 p 4 The equivalent points c, b and a are used as triangle vertexes to form an equivalent surface; p in (m) 1 、p 2 、p 4 Is less than zero p 3 If the implicit function value of (2) is greater than zero, extracting p respectively 1 p 3 、p 2 p 3 、p 3 p 4 The equivalent points a, b and c are used as triangle vertexes to form an equivalent surface; p in (n) 1 、p 3 、p 4 Is smaller than zero, p 2 If the implicit function value of (2) is greater than zero, extracting p respectively 1 p 2 、p 2 p 3 、p 2 p 4 The equivalent points c, a and b are used as triangle vertexes to form an equivalent surface; p in (o) 2 、p 3 、p 4 Is smaller than zero, p 1 If the implicit function value of (2) is greater than zero, extracting p respectively 1 p 2 、p 1 p 3 、p 1 p 4 The equivalent points c, a and b are used as triangle vertexes to form an equivalent surface; the implicit function values of the four points in (p) are all smaller than zero, and the tetrahedron is positioned outside the implicit curved surface without equivalent points; and extracting the equivalent points on each side of the tetrahedron element, traversing all the tetrahedron elements in sequence, and extracting the equivalent points to obtain a final equivalent point set. In the process of extracting the equivalent points, a linear interpolation method is used, and under the condition of obtaining the implicit function value of the equivalent points to be extracted, the position coordinates on two points of the tetrahedron element edge are usedX、Y、ZLinear interpolation is carried out on the normal vector and the corresponding implicit surface function value to obtain an equivalent pointX、Y、ZThe coordinates and normal vector of the sample,the linear interpolation formula is:
wherein, the liquid crystal display device comprises a liquid crystal display device,coordinate value of equivalent point->、/>Coordinate values of two endpoints of a straight line where the equivalent point is located, < + >>For extracting the threshold value of the isosurface +.>A solving expression for the implicit surface function;
and according to the equivalence point set, triangular surface patch drawing is sequentially carried out in each tetrahedron element, and an equivalence surface corresponding to the geologic body of the target area is obtained.
Specifically, step 5 includes:
in the embodiment of the invention, tetrahedron cutting and extraction of tetrahedrons inside the implicit curved surface are carried out, and when the implicit curved surface is judged to pass through the tetrahedron elements, the tetrahedron elements are cut on the basis of the extraction of the isosurface, so that basic units such as tetrahedrons, triangular prisms and the like are obtained. The tetrahedron element is cut in sixteen cases, the cutting mode is shown in fig. 3, sixteen cutting diagrams from (a) - (p) in fig. 3 are in one-to-one correspondence with sixteen equivalent point extraction cases from (a) - (p) in fig. 2.
Carrying out tetrahedron cutting and internal tetrahedron extraction according to the implicit surface function values of four points of the regular tetrahedron; firstly, defining tetrahedral elements with constraint values larger than zero in four points as absolute internal tetrahedrons, and defining tetrahedral elements with constraint values smaller than or equal to zero in at least one of the four points as boundary tetrahedrons; taking the isosurface as an interface, and carrying out tetrahedron cutting on the boundary tetrahedron to obtain a part of the boundary tetrahedron positioned outside the implicit curved surface and a part positioned inside the implicit curved surface; and (4) reserving the part of the boundary tetrahedron positioned in the implicit curved surface, completing the segmentation of the boundary tetrahedron, and then judging the tetrahedron with four points positioned in the implicit curved surface as an absolute tetrahedron, wherein the part of the absolute inner tetrahedron positioned in the implicit curved surface and the boundary tetrahedron form a geological body three-dimensional solid model of the target area, and the cross section of the geological body three-dimensional solid model is shown in figure 4.
Specifically, after step 5, the method further comprises:
and calculating the three-dimensional solid model of the geologic body by using an Euler tetrahedral formula to obtain the volume of the geologic body.
In the embodiment of the invention, the volume of the geologic body is as follows:
the constructed three-dimensional solid model of the geologic body is composed of tetrahedrons, triangular prisms and other elements;
for tetrahedral elements, calculating the volume of the tetrahedral using the euler tetrahedral formula;
for triangular prism elements, the triangular prism elements are firstly divided into three tetrahedron elements, and then the volume of the triangular prism is calculated by using Euler tetrahedron formula.
The Euler tetrahedron volume calculation formula is used for calculating according to the length remembering of six sides of the tetrahedron element to obtain the volume of the tetrahedron; defining the origin as O, the coordinates of the points A, B, C are respectively,/>,/>The lengths of six edges of the tetrahedron O-ABC are respectively 1 and m, n, p, q, r, and the volume calculation formula of the tetrahedron is as follows:
and calculating the volumes of all tetrahedrons and triangular prisms according to the formula, and summing to finally obtain the volume of the geologic body.
According to the embodiment of the invention, the boundary point information data of the surface of the geologic body is obtained by obtaining the drilling exploration data of the target area; constructing an implicit curved surface function model based on the boundary point information data, and solving the implicit curved surface function model to obtain a geological implicit curved surface function model; constructing a plurality of cube elements according to geological data of a target area and preset grid resolution, and dividing each cube element to obtain a plurality of tetrahedron elements; substituting coordinate values of four points on the tetrahedron element into a geological implicit surface function model for calculation according to each tetrahedron element to obtain implicit surface function values corresponding to the four points on all the tetrahedron elements; cutting all tetrahedral elements according to the function values of each implicit curved surface to obtain a geological three-dimensional entity model of the target area; the method solves the problems of lower accuracy of the geologic body and lower matching degree with the geologic body in the traditional method, establishes a geologic entity model which is more matched with the actual geologic body, and improves modeling accuracy and matching degree of the geologic body.
The embodiment of the invention also provides a geologic body modeling device based on the mobile tetrahedron, which comprises:
the acquisition module is used for acquiring drilling exploration data of the target area to obtain boundary point information data of the surface of the geologic body;
the construction module is used for constructing an implicit curved surface function model based on the boundary point information data, and solving the implicit curved surface function model to obtain a geological implicit curved surface function model;
the dividing module is used for constructing a plurality of cube elements according to geological data of the target area and preset grid resolution, and dividing each cube element to obtain a plurality of tetrahedron elements;
the calculation module is used for substituting coordinate values of four points on the tetrahedron element into the geological implicit surface function model for calculation according to each tetrahedron element to obtain implicit surface function values corresponding to the four points on all the tetrahedron elements;
and the cutting module is used for cutting all tetrahedron elements according to the function values of each implicit curved surface to obtain a geological three-dimensional entity model of the target area.
It should be noted that, because the content of information interaction and execution process between the above devices/units is based on the same concept as the method embodiment of the present invention, specific functions and technical effects thereof may be referred to in the method embodiment section, and will not be described herein.
It will be apparent to those skilled in the art that, for convenience and brevity of description, only the above-described division of the functional units and modules is illustrated, and in practical application, the above-described functional distribution may be performed by different functional units and modules according to needs, i.e. the internal structure of the apparatus is divided into different functional units or modules to perform all or part of the above-described functions. The functional units and modules in the embodiment may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit, where the integrated units may be implemented in a form of hardware or a form of a software functional unit. In addition, specific names of the functional units and modules are only for convenience of distinguishing from each other, and are not used for limiting the protection scope of the embodiments of the present invention. The specific working process of the units and modules in the above system may refer to the corresponding process in the foregoing method embodiment, which is not described herein again.
The embodiment of the invention also provides a computer readable storage medium, wherein the computer readable storage medium stores a computer program, and the computer program realizes a moving tetrahedron-based geologic body modeling method when being executed by a processor.
The integrated units, if implemented in the form of software functional units and sold or used as stand-alone products, may be stored in a computer readable storage medium. Based on such understanding, the implementation of all or part of the flow of the method of the foregoing embodiments of the present invention may be accomplished by a computer program to instruct related hardware, where the computer program may be stored in a computer readable storage medium, and the computer program may implement the steps of each of the foregoing method embodiments when executed by a processor. Wherein the computer program comprises computer program code which may be in source code form, object code form, executable file or some intermediate form etc. The computer readable medium may include at least: any entity or device capable of carrying computer program code to construct an apparatus/terminal equipment, recording medium, computer Memory, read-Only Memory (ROM), random access Memory (RAM, random Access Memory), electrical carrier signals, telecommunications signals, and software distribution media. Such as a U-disk, removable hard disk, magnetic or optical disk, etc. In some jurisdictions, computer readable media may not be electrical carrier signals and telecommunications signals in accordance with legislation and patent practice.
The embodiment of the invention also provides a terminal device which comprises a memory, a processor and a computer program stored in the memory and capable of running on the processor, wherein the processor realizes the moving tetrahedron-based geologic body modeling method when executing the computer program.
It should be noted that the terminal device may be a mobile phone, a tablet computer, a notebook computer, an Ultra mobile personal computer (UMPC, ultra-mobile Personal Computer), a netbook, a personal digital assistant (PDA, personal Digital Assistant), or the like, and the terminal device may be a station (ST, stand) in a WLAN, for example, a cellular phone, a cordless phone, a session initiation protocol (SIP, session Initiation Protocol) phone, a wireless local loop (WLL, wireless Local Loop) station, a personal digital processing (PDA, personal Digital Assistant) device, a handheld device having a wireless communication function, a computing device, or other processing device connected to a wireless modem, a computer, a laptop computer, a handheld communication device, a handheld computing device, a satellite wireless device, or the like. The embodiment of the invention does not limit the specific type of the terminal equipment.
The processor may be a central processing unit (CPU, central Processing Unit), but may also be other general purpose processors, digital signal processors (DSP, digital Signal Processor), application specific integrated circuits (ASIC, application Specific Integrated Circuit), off-the-shelf programmable gate arrays (FPGA, field-Programmable Gate Array) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, or the like. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The memory may in some embodiments be an internal storage unit of the terminal device, such as a hard disk or a memory of the terminal device. The memory may in other embodiments also be an external storage device of the terminal device, such as a plug-in hard disk provided on the terminal device, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card), etc. Further, the memory may also include both an internal storage unit and an external storage device of the terminal device. The memory is used to store an operating system, application programs, boot loader (BootLoader), data, and other programs, etc., such as program code for the computer program, etc. The memory may also be used to temporarily store data that has been output or is to be output.
It should be noted that, because the content of information interaction and execution process between the above devices/units is based on the same concept as the method embodiment in the embodiment of the present invention, specific functions and technical effects thereof may be referred to in the method embodiment section, and will not be described herein.
While the foregoing is directed to the preferred embodiments of the present invention, it will be appreciated by those skilled in the art that various modifications and adaptations can be made without departing from the principles of the present invention, and such modifications and adaptations are intended to be comprehended within the scope of the present invention.
Claims (10)
1. A mobile tetrahedron-based geologic body modeling method, comprising:
step 1, acquiring drilling exploration data of a target area to obtain boundary point information data of the surface of a geologic body;
step 2, an implicit surface function model is built based on the boundary point information data, and the implicit surface function model is solved to obtain a geological implicit surface function model;
step 3, constructing a plurality of cube elements according to the geological data of the target area and the preset grid resolution, and dividing each cube element to obtain a plurality of tetrahedron elements;
step 4, substituting coordinate values of four points on each tetrahedron element into the implicit curved surface function model of the geologic body for calculation to obtain implicit curved surface function values corresponding to the four points on all tetrahedron elements;
and 5, cutting all tetrahedral elements according to the function values of each implicit curved surface to obtain a geological three-dimensional entity model of the target area.
2. The mobile tetrahedron-based geobody modeling method of claim 1, wherein step 1 comprises:
acquiring borehole exploration data of a target area;
carrying out combined sample analysis on the drilling exploration data to obtain lithology samples of a target area;
and processing different adjacent lithology samples in the target area to obtain boundary point information data of the surface of the geologic body.
3. The mobile tetrahedron-based geobody modeling method of claim 1, further comprising, prior to said step 5:
according to the relation between the implicit surface function values corresponding to the four points on all tetrahedron elements and the preset grid size, solving the four points on all tetrahedron elements by a central difference method to obtain gradient values corresponding to the points;
and taking the gradient value corresponding to each point position as the normal vector of four point positions on all tetrahedron elements.
4. A mobile tetrahedron-based geobody modeling method according to claim 3, wherein the expression of the central difference method is:
solving four points on all tetrahedron elements by a central difference method to obtain gradient values corresponding to the points as follows:
5. A mobile tetrahedron-based geobody modeling method according to claim 3, further comprising, prior to said step 5:
according to implicit surface function values corresponding to four points on all tetrahedron elements, respectively extracting equivalent points on all sides of each tetrahedron element to obtain an equivalent point set;
and according to the contour point set, triangular surface patch drawing is sequentially carried out in each tetrahedron element, and a contour surface corresponding to the geologic body of the target area is obtained.
6. The mobile tetrahedron-based geobody modeling method of claim 5, wherein said step 5 comprises:
defining tetrahedral elements with constraint values of four points being greater than zero as absolute internal tetrahedrons,
defining tetrahedral elements with constraint values smaller than or equal to zero in at least one of the four points as boundary tetrahedrons;
cutting tetrahedrons on the boundary tetrahedrons by taking the isosurface as an interface to obtain a part of the boundary tetrahedrons positioned in the implicit curved surface;
and the absolute internal tetrahedron and the part of the boundary tetrahedron, which is positioned in the implicit curved surface, form a geological body three-dimensional entity model of the target area.
7. The mobile tetrahedron-based geobody modeling method of claim 6, further comprising, after step 5:
and calculating the three-dimensional solid model of the geologic body by using an Euler tetrahedral formula to obtain the volume of the geologic body.
8. A mobile tetrahedron-based geologic body modeling apparatus, comprising:
the acquisition module is used for acquiring drilling exploration data of the target area to obtain boundary point information data of the surface of the geologic body;
the construction module is used for constructing an implicit curved surface function model based on the boundary point information data, and solving the implicit curved surface function model to obtain a geological implicit curved surface function model;
the dividing module is used for constructing a plurality of cube elements according to the geological data of the target area and the preset grid resolution, and dividing each cube element to obtain a plurality of tetrahedron elements;
the calculation module is used for substituting coordinate values of four points on each tetrahedron element into the geological implicit surface function model for calculation to obtain implicit surface function values corresponding to the four points on all the tetrahedron elements;
and the cutting module is used for cutting all tetrahedral elements according to the function values of each implicit curved surface to obtain a geological three-dimensional entity model of the target area.
9. A computer readable storage medium storing a computer program, wherein the computer program when executed by a processor implements the mobile tetrahedron based geobody modeling method according to any one of claims 1 to 7.
10. A terminal device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that the processor implements the mobile tetrahedron based geobody modeling method according to any one of claims 1 to 7 when the computer program is executed.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310499952.1A CN116206079B (en) | 2023-05-06 | 2023-05-06 | Moving tetrahedron-based geologic body modeling method and related equipment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310499952.1A CN116206079B (en) | 2023-05-06 | 2023-05-06 | Moving tetrahedron-based geologic body modeling method and related equipment |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116206079A CN116206079A (en) | 2023-06-02 |
CN116206079B true CN116206079B (en) | 2023-06-30 |
Family
ID=86517751
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310499952.1A Active CN116206079B (en) | 2023-05-06 | 2023-05-06 | Moving tetrahedron-based geologic body modeling method and related equipment |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116206079B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116630568B (en) * | 2023-07-24 | 2023-11-03 | 航天宏图信息技术股份有限公司 | Geological three-dimensional model generation method and device based on point, line and surface structures |
CN116977588B (en) * | 2023-09-21 | 2023-12-19 | 煤炭科学研究总院有限公司 | Cutting curve generation method and device based on voxel grid geological model |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5390035A (en) * | 1992-12-23 | 1995-02-14 | International Business Machines Corporation | Method and means for tetrahedron/octahedron packing and tetrahedron extraction for function approximation |
EP1024436A2 (en) * | 1999-01-29 | 2000-08-02 | Mitsubishi Denki Kabushiki Kaisha | Method for generating graphic objects represented as surface elements |
CN102194253A (en) * | 2011-06-15 | 2011-09-21 | 北京航空航天大学 | Method for generating tetrahedron gridding for three-dimensional geological structure |
CN102254349A (en) * | 2011-06-30 | 2011-11-23 | 华东师范大学 | Method for constructing three-dimensional physical model of sedimentary stratum system by using drilling data |
CN111986325A (en) * | 2020-08-25 | 2020-11-24 | 东北大学 | Implicit sectioning method for three-dimensional geologic body model |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9902114B2 (en) * | 2014-01-09 | 2018-02-27 | Siemens Product Lifecycle Management Software Inc. | Method for creating three dimensional lattice structures in computer-aided design models for additive manufacturing |
FR3016983A1 (en) * | 2014-01-28 | 2015-07-31 | IFP Energies Nouvelles | METHOD FOR CONSTRUCTING A VOLUME MESH FOR MODELING GEOLOGICAL STRUCTURES |
WO2016120225A1 (en) * | 2015-01-28 | 2016-08-04 | Koninklijke Philips N.V. | Finite element modeling of anatomical structure |
-
2023
- 2023-05-06 CN CN202310499952.1A patent/CN116206079B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5390035A (en) * | 1992-12-23 | 1995-02-14 | International Business Machines Corporation | Method and means for tetrahedron/octahedron packing and tetrahedron extraction for function approximation |
EP1024436A2 (en) * | 1999-01-29 | 2000-08-02 | Mitsubishi Denki Kabushiki Kaisha | Method for generating graphic objects represented as surface elements |
CN102194253A (en) * | 2011-06-15 | 2011-09-21 | 北京航空航天大学 | Method for generating tetrahedron gridding for three-dimensional geological structure |
CN102254349A (en) * | 2011-06-30 | 2011-11-23 | 华东师范大学 | Method for constructing three-dimensional physical model of sedimentary stratum system by using drilling data |
CN111986325A (en) * | 2020-08-25 | 2020-11-24 | 东北大学 | Implicit sectioning method for three-dimensional geologic body model |
Non-Patent Citations (1)
Title |
---|
移动立方体算法的地质体三维空间形态模拟;邹艳红;何建春;;测绘学报(06);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN116206079A (en) | 2023-06-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116206079B (en) | Moving tetrahedron-based geologic body modeling method and related equipment | |
Rottensteiner et al. | Automated delineation of roof planes from lidar data | |
CN109753707B (en) | Method for extracting stratum boundary line by utilizing exploration line profile to develop three-dimensional modeling | |
CN104966317B (en) | A kind of three-dimensional method for automatic modeling based on ore body contour line | |
CN102999937A (en) | Curved planar reconstruction method for cardiac scattered-point cloud data | |
CN111932669A (en) | Deformation monitoring method based on slope rock mass characteristic object | |
CN116030218A (en) | Tetrahedral mesh dividing method, device, system and storage medium | |
CN115861527A (en) | Method and device for constructing live-action three-dimensional model, electronic equipment and storage medium | |
CN115457212A (en) | Tree image processing method and device, terminal equipment and storage medium | |
CN110675498B (en) | Long and large strip-shaped three-dimensional geological environment modeling method and device and terminal equipment | |
CN115546116A (en) | Method and system for extracting and calculating spacing of discontinuous surface of fully-covered rock mass | |
Chaudhuri et al. | Scalable computation of distributions from large scale data sets | |
CN109163708B (en) | Optimization method, device and medium for unmanned aerial vehicle large-scale topographic map house line | |
CN114785998A (en) | Point cloud compression method and device, electronic equipment and storage medium | |
CN109191484B (en) | A method of the rapidly extracting planar chip from airborne laser radar point cloud | |
CN110246223A (en) | A kind of tunnel modeling method and device | |
CN116012613B (en) | Method and system for measuring and calculating earthwork variation of strip mine based on laser point cloud | |
CN115792900A (en) | Deformation detection method and system for geological sensitive area of power transmission channel | |
CN116452764A (en) | Airborne laser radar point cloud hierarchical rarefaction method integrating terrain information | |
CN112907649B (en) | Method, device, equipment and storage medium for measuring volume of soil and stone in pipe ditch | |
CN115481268A (en) | Parallel line wall data identification method, device, equipment and storage medium | |
Sud et al. | Interactive computation of discrete generalized voronoi diagrams using range culling | |
Zhang et al. | Extraction, quantification and visualization of protein pockets | |
CN115328373B (en) | Deterministic coefficient calculation method and system | |
CN115034151B (en) | Fault vertical dominant migration channel analysis method and system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |