Method and device for realizing three-dimensional geological modeling using point layer thickness grid model

Abstract

The present invention relates to the fields of mathematical geology, geographic information and computer technology, and in particular to a method, device, equipment and computer storage medium for realizing three-dimensional geological modeling by using a point-layer thickness grid model. A point-layer thickness grid data model based on polygonal grids and generalized prismatic elements is used to express a three-dimensional geological model, accurately showing the three-dimensional spatial form of the stratum structure, geological structure and attribute characteristic distribution. A point-layer thickness interpolation calculation model is used to realize the generation and combination of stratum interface grids, which solves the problems of self-intersection of the top and bottom surfaces of thin stratum structure models and upper and lower penetration of multi-stratum structure models in the three-dimensional geological structure modeling method. A point-layer thickness division calculation model is used to realize the division and attribute interpolation of stratum generalized prismatic elements, which solves the problems of lack of stratum interface information in the attribute body model, non-overlapping of the stratum interfaces between the attribute body model and the structure model, and insufficient adaptability of the attribute body model to changes in stratum thickness in the three-dimensional geological attribute modeling method.

Description

Technical Field

[0001] The invention relates to the fields of mathematical geology, geographic information and computer technology, and in particular to a method, device, equipment and computer storage medium for realizing three-dimensional geological modeling by using a point layer thickness grid model. Background Art

[0002] The following problems often occur when constructing a three-dimensional geological model using existing technologies:

[0003] The top and bottom surfaces of the thin stratum structure model intersect: During the three-dimensional geological structure modeling process, due to the inadequacy of the stratum interface grid generation method, the top and bottom surfaces of the thin strata may intersect.

[0004] The structural model penetrates the upper and lower layers at the pinch-out position of the stratum: In the process of three-dimensional geological structure modeling, the adjacent strata are regarded as independent strata for modeling, and the stratum interfaces are not treated or treated improperly, especially the lack of proper treatment of the stratum missing areas in the stratum interfaces, resulting in the intersection of the bottom surface of the upper stratum and the top surface of the lower stratum, thus forming the interpenetration of the upper and lower strata.

[0005] The geological attribute volume model lacks stratigraphic interfaces, and its surface does not coincide with the surface of the structural model: In the process of three-dimensional geological attribute modeling, multiple strata are treated as a whole, and the overall outer contour surfaces of multiple strata are obtained at one time. Volume element segmentation and attribute interpolation are carried out based on the outer contour surfaces, ignoring the interface division between different strata, resulting in the lack of stratigraphic interface information in the generated attribute volume model; in the process of volume element segmentation, the top and bottom surfaces of the strata are not used to accurately constrain the top and bottom surfaces of the volume elements, resulting in the stratigraphic interface of the generated attribute volume model not coinciding with the stratigraphic interface of the structural model. The typical manifestation is that the surface of the attribute volume model does not coincide with the surface of the structural model.

[0006] The attribute body model is not able to adapt to the change of stratum thickness: In the process of 3D geological attribute modeling, a single resolution is used to generate the attribute body model. If the resolution is too low, the number of vertical subdivision layers of thin strata may be too small, or even the entire stratum may be subdivided into one layer, making it difficult to accurately reflect the vertical change characteristics of the attributes in the thin strata; if the resolution is set too high to process thin strata, the number of vertical subdivision layers of thick strata may be too large, thereby increasing the unnecessary amount of attribute space interpolation calculations, and may even cause the attribute interpolation algorithm calculation process to not converge. A single-resolution 3D geological attribute body model is difficult to meet the data processing requirements of both thin and thick strata.

[0007] In summary, how to construct a three-dimensional geological model with correct stratigraphic structure, accurate stratigraphic interface and adaptability to changes in stratigraphic thickness is a problem that needs to be solved at present. Summary of the invention

[0008] To this end, the technical problem to be solved by the present invention is to overcome the problems of incorrect stratum structure, inaccurate stratum interface and insufficient adaptability to stratum thickness changes in the prior art for constructing three-dimensional geological models.

[0009] In order to solve the technical problem, the present invention provides a method for realizing three-dimensional geological modeling, comprising:

[0010] The three-dimensional geological model is expressed by using a point-layer-thickness grid data model based on polygonal grids and generalized prismatic elements;

[0011] The point-by-point layer thickness interpolation calculation model is used to realize the formation interface grid generation and combination modeling, and construct a three-dimensional geological structure model;

[0012] The point-by-point layer thickness division calculation model is used to realize the generalized prism element segmentation and attribute interpolation of the stratum, and to construct a three-dimensional geological attribute body model.

[0013] Preferably, the three-dimensional geological model is expressed by using a point-layer thickness grid data model based on polygonal grids and generalized prismatic voxels, including:

[0014] Polygonal meshes are used to express the three-dimensional geological structure model: for the m strata stacked from top to bottom in the modeling area, m+1 polygonal meshes distributed from top to bottom are used to express the m+1 stratum interfaces in the strata.

[0015] The i-th and i+1-th polygonal meshes at any upper and lower positions represent the top and bottom surfaces of the i-th stratum, respectively.

[0016] All polygon meshes have the same number of vertices, the same number of polygons, the same number of edges, and the same way of vertices forming polygons;

[0017] The generalized prismatic volume element model is used to express the three-dimensional geological attribute volume model: for any layer L among the m layers i , using n layers of generalized prism elements stacked from top to bottom to express the space between the top and bottom surfaces of the stratum, all generalized prisms have the same number of edges, and each generalized prism element has 0 or one or more attribute values,

[0018] The top surface of each generalized prism in the uppermost voxel layer is related to the stratum L i The polygons in the top polygonal mesh overlap one by one, and the bottom surface of each generalized prism in the lowest voxel layer overlaps with the bottom surface of the stratum L. i The polygons in the bottom polygon mesh overlap one by one.

[0019] For any two adjacent voxel layers, the bottom surfaces of the generalized prisms in the upper voxel layer coincide with the top surfaces of the generalized prisms in the lower voxel layer.

[0020] Among them, m is an integer greater than 0, n is an integer greater than 1, the i-th and i+1-th polygonal meshes dynamically refer to two adjacent polygonal meshes in upper and lower positions being used, and the i-th stratum dynamically refers to the stratum being used.

[0021] Preferably, the strata are divided according to preset rules, and the preset rules include but are not limited to geological age, coal rock characteristics, rock mass characteristics, and geotechnical characteristics; the generalized prisms include conventional prisms and variant prisms, and the variant prisms include but are not limited to one or more of the following situations: the top surface is not parallel to the bottom surface, the top surface is not perpendicular to the edge, the bottom surface is not perpendicular to the edge, and the length of the edge is zero.

[0022] Preferably, the point-level layer thickness interpolation calculation model is used to realize the formation interface grid generation and combination modeling to construct a three-dimensional geological structure model, including:

[0023] For a three-dimensional geological structure model expressed by polygonal meshes, the three-dimensional coordinates of one or more sample points with stratigraphic boundary information on the top surface of the uppermost stratum and the polygonal mesh M of the top surface of the uppermost stratum are used. i t Any j-th vertex V j The plane coordinates of the vertex V are determined by using the spatial interpolation method. j The vertical coordinates of the polygonal mesh M are obtained i t The three-dimensional coordinates of each vertex;

[0024] According to any stratum L i The three-dimensional coordinates of one or more sample points with stratum boundary information on the top surface and the plane coordinate position of the stratum L at the above sample points i The thickness of the formation L i Top polygon mesh M i t Any j-th vertex V j The plane coordinates of the vertex V are determined by using the spatial interpolation method. j The plane coordinate position of the stratum L i The thickness of the polygonal mesh M is obtained i t The plane coordinate position of each vertex in the stratum L i Thickness;

[0025] Take the vertex V j The vertical coordinate of the vertex V minus or plus j The plane coordinate position of the stratum L i The thickness of the formation L i Bottom polygon mesh M ib The jth vertex V j The vertical coordinates of the polygonal mesh M are obtained i b The three-dimensional coordinates of each vertex in;

[0026] The three-dimensional coordinates include plane coordinates and vertical coordinates, and the jth vertex dynamically refers to the vertex being used in the polygonal mesh.

[0027] Preferably, the method of using a point-level layer thickness interpolation calculation model to realize formation interface grid generation and combination modeling to construct a three-dimensional geological structure model also includes:

[0028] For any stratum L in the 3D geological structure model expressed by polygonal meshes i , generate a polygon mesh M with the same vertex and polygon data as the top polygon mesh of the stratum Li i t , as the formation L i The top surface drawing data; if the formation L i If there is a missing layer in the stratum, then from the polygonal mesh M i t Delete all polygons whose vertices are located in the missing area of ​​the formation, and then remove the polygon mesh M i t Delete vertices that are not used by any polygon;

[0029] Generate a layer with L i The bottom polygon mesh has the same vertex and polygon data as the polygon mesh M i b , as the formation L i Bottom surface drawing data; if the stratum L i If there is a missing area in the stratum, then from the polygonal mesh M i b Delete all polygons whose vertices are located in the missing area of ​​the formation, and then remove the polygon mesh M i b Delete vertices that are not used by any polygon;

[0030] Extraction of formation L i Top and bottom drawing data polygon mesh M i t and M i b The boundary line in the expression layer L is generated according to the boundary line iA polygonal mesh of the surrounding interface; the boundary line is one or more broken lines without repeated vertices formed by connecting the open edges in the polygonal mesh according to the shared vertices between the head and the tail; the open edge is a polygonal edge used by only one polygon in the polygonal mesh.

[0031] Preferably, a point-level layer thickness division calculation model is used to implement stratum generalized prism element segmentation and attribute interpolation, and to construct a three-dimensional geological attribute body model, including:

[0032] For any stratum L in the 3D geological structure model expressed by polygonal meshes i , determine the number of voxel layers n that need to be divided;

[0033] In the stratum L i Top surface drawing data polygon mesh M i t Any k-th polygon and stratum L in i Bottom drawing data polygonal mesh M i b The spatial position between the k-th polygons in the , generates n-1 subdivided polygons arranged layer by layer from top to bottom;

[0034] The polygon mesh M i t The kth polygon in the mesh and the subdivided polygon at the top layer are used as the top and bottom surfaces of the generalized prism, respectively, to generate a generalized prism element; the mesh polygon M i b The kth polygon in the subdivided polygon and the subdivided polygon at the bottom layer are used as the top and bottom surfaces of the generalized prism respectively to generate a generalized prism element; each pair of two polygons in the subdivided polygons that are in the upper and lower adjacent positions are

[0035] The polygons are used as the top and bottom surfaces of the generalized prism to generate a generalized prism element;

[0036] Wherein, n is an integer greater than 1, and k dynamically refers to the polygons being used in the polygon mesh.

[0037] Preferably, the point layer thickness division calculation model is used to realize the generalized prism element segmentation and attribute interpolation of the stratum, and to construct a three-dimensional geological attribute body model, which also includes:

[0038] For any generalized prismatic element P in the three-dimensional geological attribute volume model expressed by the generalized prismatic element grid model k , according to the three-dimensional coordinates and attribute values ​​of the sample points with attribute information in the modeling area, and the voxel P k The three-dimensional coordinates of any point in the internal space or on the surface are determined by spatial interpolation calculation method to determine the volume element P kHas the attribute value.

[0039] The present invention also provides a device for implementing three-dimensional geological modeling, comprising:

[0040] A three-dimensional geological model data expression module is used to express a three-dimensional geological model using a point layer thickness grid data model based on polygonal grids and generalized prismatic elements;

[0041] The 3D geological structure modeling module is used to realize the grid generation and combination of stratum interface by using the point layer thickness interpolation calculation model to construct a 3D geological structure model;

[0042] The three-dimensional geological attribute modeling module is used to realize the generalized prism element segmentation and attribute interpolation of the stratum by adopting the point layer thickness division calculation model, and to construct a three-dimensional geological attribute body model.

[0043] The present invention also provides a device for implementing three-dimensional geological modeling, comprising:

[0044] Memory for storing computer programs;

[0045] A processor is used to implement the method steps for realizing three-dimensional geological modeling when executing the computer program.

[0046] The present invention also provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a computer, the steps of the method for realizing three-dimensional geological modeling are implemented.

[0047] The technical solution of the present invention has the following advantages over the prior art:

[0048] The method and device for implementing three-dimensional geological modeling by using a point-layer thickness grid model described in the present invention can correctly generate a three-dimensional geological structure model and / or a three-dimensional geological attribute body model in one go, and no human intervention is required during the entire model generation process. The core of this technology is to use a point-layer thickness grid data model based on polygonal grids and generalized prismatic elements to express a three-dimensional geological model, accurately displaying the three-dimensional spatial form of the stratigraphic structure, geological structure, and attribute characteristic distribution; to use a point-layer thickness interpolation calculation model to realize stratigraphic interface grid generation and combination into a model, to construct a three-dimensional geological structure model, to realize the constraint of stratigraphic thickness on the vertical coordinates of the top and / or bottom surfaces of the three-dimensional geological structure model, and to solve the problem of self-intersection of the top and bottom surfaces of the thin stratigraphic structure model; to realize that the bottom surface of the upper stratum in adjacent strata has the same number but different models as the top surface of the lower stratum, and to solve the problem of three-dimensional geological structure of multiple strata. The problem of layer penetration in the model is solved; the point layer thickness division calculation model is used to realize the generalized prism volume element segmentation and attribute interpolation of the stratum, and a three-dimensional geological attribute body model is constructed, so that the three-dimensional geological attribute body model contains accurate stratum interface information, and the problem of lack of stratum interface in the three-dimensional geological attribute body model and / or the stratum interface between the three-dimensional geological attribute body model and the three-dimensional geological structure model does not overlap is solved; different strata in the three-dimensional geological structure model are divided into volume element grid models with the same or different resolutions as needed, which solves the problem of insufficient adaptability of the single-resolution three-dimensional geological attribute body model to changes in stratum thickness.

[0049] In summary, the technology provided by the present invention effectively solves the problems of incorrect stratigraphic structure, inaccurate stratigraphic interface and insufficient adaptability to changes in stratigraphic thickness in constructing three-dimensional geological models in the prior art through a series of innovative methods and steps, realizes an automated multi-stratum, multi-resolution, high-precision three-dimensional geological modeling process, and effectively improves the accuracy of model result data, the efficiency of the model generation process and the convenience of human-computer interaction operations. BRIEF DESCRIPTION OF THE DRAWINGS

[0050] In order to make the content of the present invention more clearly understood, the present invention is further described in detail below according to specific embodiments of the present invention in conjunction with the accompanying drawings, wherein:

[0051] Figure 1 It is a flow chart of a method for realizing three-dimensional geological modeling provided by the present invention;

[0052] Figure 2 It is a schematic diagram of a method for expressing a three-dimensional geological model using a point-layer-thickness grid data model based on polygonal grids and generalized prismatic elements;

[0053] Figure 3 It is a schematic diagram of a method for realizing formation interface grid generation and combination modeling and constructing a three-dimensional geological structure model by using a point layer thickness interpolation calculation model;

[0054] Figure 4 This is a schematic diagram of a method for using a point-based layer thickness division calculation model to implement stratum generalized prism element segmentation and attribute interpolation and to construct a three-dimensional geological attribute body model. DETAILED DESCRIPTION

[0055] The core of the present invention is to provide a method, device, equipment and computer storage medium for realizing three-dimensional geological modeling using a point-layer thickness grid model, thereby realizing an automated multi-stratum, multi-resolution, high-precision three-dimensional geological modeling process, and effectively improving the accuracy of the model result data of three-dimensional geological modeling, the efficiency of the model generation process and the convenience of human-computer interaction.

[0056] In order to enable those skilled in the art to better understand the scheme of the present invention, the present invention is further described in detail below in conjunction with the accompanying drawings and specific implementation methods. Obviously, the described embodiments are only part of the embodiments of the present invention, rather than all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention.

[0057] Please refer to Figure 1. Figure 1 The following is a flow chart of a method for implementing three-dimensional geological modeling provided by the present invention; the specific operation steps are as follows:

[0058] S101: Use point-layer thickness grid data model based on polygonal grid and generalized prismatic element to express 3D geological model;

[0059] S102: Using the point layer thickness interpolation calculation model to realize the formation interface grid generation and combination modeling, and construct a three-dimensional geological structure model;

[0060] S103: A point-by-point layer thickness division calculation model is used to implement stratum generalized prism element segmentation and attribute interpolation, and a three-dimensional geological attribute body model is constructed.

[0061] Based on the above embodiments, this embodiment describes step S101 in detail:

[0062] Polygonal meshes are used to express the three-dimensional geological structure model: for the m strata stacked from top to bottom in the modeling area, (m+1) polygonal meshes distributed from top to bottom are used to express the (m+1) stratum interfaces in the above strata.

[0063] The i-th and (i+1)-th polygonal meshes at any upper and lower positions represent the top and bottom surfaces of the i-th stratum, respectively.

[0064] All polygon meshes have the same number of vertices, the same number of polygons, the same way of vertices forming polygons, and the same number of edges in all polygons in the meshes;

[0065] The generalized prismatic volume element model is used to express the three-dimensional geological attribute volume model: for any one of the m strata L i , n layers of generalized prism elements stacked from top to bottom are used to express the space between the top and bottom surfaces of the stratum, all generalized prisms have the same number of edges, and each generalized prism element may have one or more attribute values, or may not have an attribute value,

[0066] The top surface of each generalized prism in the uppermost voxel layer is related to the stratum L i The polygons in the top polygonal mesh overlap one by one, and the bottom surface of each generalized prism in the lowest voxel layer overlaps with the bottom surface of the stratum L. i The polygons in the bottom polygon mesh overlap one by one.

[0067] For any two adjacent voxel layers, the bottom surfaces of the generalized prisms in the upper voxel layer coincide with the top surfaces of the generalized prisms in the lower voxel layer.

[0068] The m is an integer greater than 0; the n is an integer greater than 1; the i-th and (i+1)-th polygonal meshes dynamically refer to two adjacent polygonal meshes in the upper and lower positions being used, and do not fixedly refer to two specific polygonal meshes; the i-th stratum dynamically refers to the stratum being used, and the sequence number is the same as the sequence number of the i-th polygonal mesh.

[0069] Based on the above embodiments, the strata are divided according to preset rules, and the preset rules include but are not limited to geological age, coal rock characteristics, rock mass characteristics, and geotechnical characteristics; the generalized prisms include conventional prisms and variant prisms, and the variant prisms include but are not limited to one or more of the following situations: the top surface is not parallel to the bottom surface, the top surface is not perpendicular to the edge, the bottom surface is not perpendicular to the edge, and the length of the edge is zero.

[0070] like Figure 2 As shown, step S101 is described with a specific example:

[0071] The three-dimensional geological structure model is expressed by quadrilateral meshes: the modeling area has three strata L1, L2, and L3 stacked from top to bottom, and four quadrilateral meshes M1 ({V1 i};i=1,2,3,...,12),M2({V2 i};i=1,2,3,...,12),M3({V3 i};i=1,2,3,...,12),M4({V4 i};i=1,2,3,...,12) respectively express the four stratigraphic interfaces existing in the above-mentioned strata, wherein {} represents a set consisting of one or more elements, i=1,2,3,... represents different elements formed by the placeholder elements in {} according to the different values ​​of i, and the subsequent expressions all have the above meanings.

[0072] M1 and M2 express the top and bottom surfaces of L1 respectively, M2 and M3 express the top and bottom surfaces of L2 respectively, M3 and M4 express the top and bottom surfaces of L3 respectively,

[0073] The number of vertices of M1, M2, M3, and M4 are all 12, the number of quadrilaterals is all 4, and the way the vertices form quadrilaterals is ({1,2,6,5};{2,3,7,6};{3,4,8,7};{5,6,10,9};{6,7,11,10};{7,8,12,11}), which means that the 12 vertices are arranged in sequence from 1 to 12, and the 4 vertices in each of the above curly brackets {} form a quadrilateral;

[0074] The three-dimensional geological attribute volume model is expressed by using a generalized tetrahedral prism grid model: Taking the stratum L1 as an example, three layers of generalized tetrahedral prism layers R1 are stacked from top to bottom. 1 ({V1 i};{T1 i};i=1,2,3,...,12}, R1 2 ({T1 i};{T2 i};i=1,2,3,...,12}), R1 3 ({T2 i};{V2 i};i=1,2,3,...,12}) expresses the space of the top and bottom surfaces of the stratum L1. The number of edges of all generalized tetrahedral prisms is 4. Each generalized tetrahedral prism element has an attribute value called porosity.

[0075] Take the generalized quadrangular prism P1 in the top voxel layer R11 1 ({V1 i};{T1 i};i=1,2,6,5) of the top quadrilateral Q1 p1t ({V1 i};i=1,2,6,5) as an example, Q1 p1t The quadrilateral Q1 in the top surface M1 of the layer L1 1 ({V1 i};i=1,2,6,5) overlap, with the lowest voxel layer R1 3 Generalized tetrahedral prism P1 in 3 ({V1 i};{T1i};i=1,2,6,5) of the bottom quadrilateral Q1 p3t ({V2 i};i=1,2,6,5) as an example, Q1 p3t The quadrilateral Q2 in the bottom surface M2 of the layer L1 1 ({V2 i};i=1,2,6,5) overlap,

[0076] The two adjacent voxel layers R1 1 and R1 2 For example, the upper voxel layer R1 1 Generalized tetrahedral prism P1 in 1 Bottom surface Q1 p1b ({T1 i};i=1,2,6,5) and the lower voxel layer R1 2 Generalized tetrahedral prism P1 in 2 Top surface Q1 P2t ({T1 i};i=1,2,6,5) overlap.

[0077] Based on the above embodiments, this embodiment describes step S102 in detail according to the top-down generation model:

[0078] For a three-dimensional geological structure model expressed by polygonal meshes, the three-dimensional coordinates (plane coordinates, vertical coordinates) of one or more sample points with stratigraphic boundary information on the top surface of the uppermost stratum and the polygonal mesh M of the top surface of the uppermost stratum are used. i t Any j-th vertex V j The plane coordinates of the vertex V are determined by using the spatial interpolation method. j The vertical coordinates of the polygonal mesh M are obtained i t The three-dimensional coordinates of each vertex (plane coordinates, vertical coordinates);

[0079] According to any stratum L i The three-dimensional coordinates (plane coordinates, vertical coordinates) of one or more sample points on the top surface with stratum boundary information and the stratum L at the plane coordinate position of the above sample points i The thickness of the formation L i Top polygon mesh M i t Any j-th vertex V j The plane coordinates of the vertex V are determined by using the spatial interpolation method. j The plane coordinate position of the stratum Li is the thickness, thereby obtaining the polygonal mesh M i tThe plane coordinate position of each vertex in the stratum L i Thickness;

[0080] Take the vertex V j The vertical coordinate of the vertex V minus or plus j The plane coordinate position of the stratum L i The thickness of the formation L i Bottom polygon mesh M i b The jth vertex V j The vertical coordinates of the polygonal mesh M are obtained i b The three-dimensional coordinates (plane coordinates, vertical coordinates) of each vertex in ;

[0081] The j-th vertex dynamically refers to the vertex being used in the polygonal mesh, and does not refer to a specific vertex in a fixed manner; the sample points may select all sample points or some of the sample points; the spatial interpolation calculation method includes but is not limited to inverse distance weighted (IDW) interpolation, various Kriging interpolations, and various trend interpolations;

[0082] For any stratum L in the 3D geological structure model expressed by polygonal meshes i , generating a layer with the formation L i The top polygon mesh has the same vertex and polygon data as the polygon mesh M i t , as the formation L i The top surface drawing data; if the formation L i If there is a missing area in the stratum, then from the polygonal mesh M i t Delete all polygons whose vertices are located in the missing area of ​​the formation, and then remove the polygon mesh M i t Delete the vertices that are not used by any polygons, and finally convert the polygon mesh M i t As the stratum L i The top surface drawing data of

[0083] Generate a layer with L i The bottom polygon mesh has the same vertex and polygon data as the polygon mesh M i b , as the formation L i Bottom surface drawing data; if the stratum L i If there is a missing area in the stratum, then from the polygonal mesh M i bDelete all polygons whose vertices are located in the missing area of ​​the formation, and then remove the polygon mesh M i b Delete the vertices that are not used by any polygons, and finally convert the polygon mesh M after the above processing i b As the stratum L i Bottom drawing data of ;

[0084] Extraction of formation L i Top and bottom drawing data polygon mesh M i t and M i b The boundary line in the expression layer L is generated according to the boundary line i A polygonal mesh of the surrounding interface; the boundary line is one or more broken lines without repeated vertices formed by connecting the open edges in the polygonal mesh according to the shared vertices between the head and the tail; the open edge is a polygonal edge used by only one polygon in the polygonal mesh.

[0085] like Figure 3 As shown, step S102 is described with a specific example:

[0086] For the three layers L1, L2, L3 and the four quadrilateral meshes M1, M2, M3, M4 in the modeling area, the vertex V1 in the quadrilateral mesh M1 on the top surface of the uppermost layer L1 is 6 For example, based on the four sample points with stratum boundary information located on the top surface of stratum L1 ({S1 i};i=1,2,6,7)'s three-dimensional coordinates ({(x s1 i ,y s1 i ,z s1 i )};i=1,2,6,7), and V1 6 The plane coordinates (x1 6 ,y1 6 ), using the inverse distance weighted (IDW) interpolation method to determine V1 6 The vertical coordinate (z1 6 ), thus obtaining V1 6 The three-dimensional coordinates (x1 6 ,y1 6 ,z1 6 ); Similarly, obtain the three-dimensional coordinates of all vertices in M1 ({(x1 i ,y1 i ,z1 i )};i=1,2,3,...,12);

[0087] Take the vertex V1 in the quadrilateral mesh M1 on the top surface of the layer L1 6 For example, based on the six sample points with stratum boundary information located on the top surface of stratum L1 ({S1 i};i=1,2,3,6,7,9)'s three-dimensional coordinates ({(x s1 i ,y s1 i ,z s1 i )}; i=1,2,3,6,7,9) and the plane coordinates of the above sample points ({(x s1 i ,y s1 i )};i=1,2,3,6,7,9) the thickness of the formation L1 ({d s1 i};i=1,2,3,6,7,9), and V1 6 The plane coordinates (x1 6 ,y1 6 ), using the polynomial trend surface interpolation method, determine the V1 6 The plane coordinates (x1 6 ,y1 6 ) The thickness of the layer L1 at position (d1 6 ); Similarly, obtain the plane coordinates of all vertices in the quadrilateral meshes M1, M2, and M3 ({(x j i ,y j i )};i=1,2,3,...,12;j=1,2,3) the thickness of the layers L1, L2, and L3 ({d j i};i=1,2,3,...,12;j=1,2,3);

[0088] Take the vertex V2 in the quadrilateral mesh M2 of the bottom surface of the layer L1 6 For example, take the vertex V1 in the quadrilateral mesh M1 on the top surface of the stratum L1 6 The vertical coordinate (z1 6 ), minus the value in V1 6 Plane coordinates (x1 6 ,y1 6 ) The thickness of the layer L1 at position (d1 6 ), and the obtained value (z1 6 -d1 6 ) as V2 6 The vertical coordinate (z2 6 ), thus obtaining V2 6 The three-dimensional coordinates (x2 6,y2 6 ,z2 6 ); Similarly, obtain the three-dimensional coordinates of all vertices in the quadrilateral meshes M2, M3, and M4 {(x j i ,y j i ,z j i );i=1,2,3,...,12;j=2,3,4};So far, the three-dimensional coordinates of all vertices in the quadrilateral meshes M1, M2, M3, and M4 have been obtained;

[0089] Taking the stratum L2 as an example, generate a quadrilateral mesh data M2 with the same vertices and quadrilaterals as the top surface quadrilateral mesh M2 of the stratum L2 t ({V2 i}; i = 1, 2, 3, ..., 12), as the top surface drawing data of layer L2; there is a layer missing area in layer L2, so from the top surface drawing data M2 t Delete the quadrilateral Q2 whose vertices are located in the missing area of ​​the formation 1 ({V2 i};i=3,4,8,7) and Q2 2 ({V2 i};i=7,8,12,11), then from the drawing data M2 t Delete the vertices not used by any quadrilateral ({V2 i};i=4,8,12), and finally the drawing data M2 after the above processing t ({V2 i}; i = 1, 2, 3, 5, 6, 7, 9, 10, 11) as the top surface drawing data of the stratum L2; similarly, the top surface drawing data M1 of the stratum L1 and L3 are obtained t 、M3 t ;

[0090] Taking the stratum L2 as an example, generate a quadrilateral mesh data M2 with the same vertices and quadrilaterals as the quadrilateral mesh M3 of the bottom surface of the stratum L2 b ({V3 i}; i = 1, 2, 3, ..., 12), as the bottom surface drawing data of layer L2; there is a layer missing area in layer L2, so the top surface drawing data M2 b Delete the quadrilateral Q3 whose vertices are located in the missing area of ​​the formation 1 ({V3 i};i=3,4,8,7) and Q3 2 ({V3 i};i=7,8,12,11), then from the drawing data M2 bDelete vertices not used by any quadrilateral ({V3 i};i=4,8,12), and finally the drawing data M2 after the above processing b ({V3 i}; i = 1, 2, 3, 5, 6, 7, 9, 10, 11) as the bottom surface drawing data of the stratum L2; similarly, the bottom surface drawing data M1 of the stratum L1 and L3 are obtained b 、M3 b ;

[0091] Taking stratum L2 as an example, extract L2 top surface drawing data M2 t Boundary line E2 t ({V2 i};i=1,2,3,5,7,9,10,11), bottom drawing data M2 b Boundary line E2 b ({V3 i}; i = 1, 2, 3, 5, 7, 9, 10, 11), generate the polygonal mesh M2 of the interface drawing data of the four surroundings of the stratum L2 according to the boundary line e ({V2 i},{V3 i}; i=1,2,3,5,7,9,10,11); with boundary line E2 t For example, E2 t is composed of quadrilateral mesh M2 t The 8 open edges in (V1 1 ,V1 2 )、(V2 2 ,V2 3 )、(V2 3 ,V2 7 )、(V2 7 ,V2 11 )、(V2 11 ,V2 10 )、(V2 10 ,V2 9 )、(V2 9 ,V2 5 )、(V2 5 ,V2 1 ) is a polyline without duplicate vertices connected by the shared vertices between the head and the tail; with open edges (V2 1 ,V2 2 ) as an example, (V2 1 ,V2 2 ) is the quadrilateral mesh M2 t There is only one quadrilateral Q({V2 i}; i = 1,2,6,5) quadrilateral edge used; similarly obtain the surrounding interface drawing data M2 of the strata L2 and L3 e 、M3 e At this point, the top surface, bottom surface, and surrounding interface mapping data of the strata L1, L2, and L3 have been obtained, thus constructing a three-dimensional geological structure model containing three strata.

[0092] Based on the above embodiments, this embodiment describes step S103 in detail:

[0093] For any stratum L in the 3D geological structure model expressed by polygonal meshes i , determine the number of voxel layers n that need to be divided;

[0094] In the stratum L i Top surface drawing data polygon mesh M i t Any k-th polygon and stratum L in i Bottom drawing data polygonal mesh M i b The spatial position between the k-th polygons in the , generates (n-1) subdivided polygons arranged layer by layer from top to bottom;

[0095] The polygon mesh M i t The kth polygon in the mesh and the subdivided polygon at the top layer are used as the top and bottom surfaces of the generalized prism, respectively, to generate a generalized prism element; the mesh polygon M i b The kth polygon in the subdivided polygon and the subdivided polygon at the bottom layer are respectively used as the top surface and the bottom surface of the generalized prism to generate a generalized prism element; each pair of two polygons in the subdivided polygons that are in the upper and lower adjacent positions are respectively used as the top surface and the bottom surface of the generalized prism to generate a generalized prism element;

[0096] The n is an integer greater than 1, and the k dynamically refers to a polygon being used in the polygonal mesh, rather than being fixed to a specific polygon;

[0097] For any generalized prismatic element P in the three-dimensional geological attribute volume model expressed by the generalized prismatic element grid model k , according to the three-dimensional coordinates (plane coordinates, vertical coordinates) and attribute values ​​of the sample points with attribute information in the modeling area, and the volume element P k The three-dimensional coordinates (plane coordinates, vertical coordinates) of any point in the internal space or on the surface are determined by spatial interpolation calculation method to determine the volume element P k The attribute value it has;

[0098] The sample points may be all sample points or some of the sample points, and the spatial interpolation calculation method may include but is not limited to inverse distance weighted (IDW) interpolation, various Kriging interpolations, and various trend interpolations.

[0099] like Figure 4 As shown, step S103 is described with a specific example:

[0100] For the three strata L1, L2, L3 and four quadrilateral meshes M1, M2, M3, M4 in the modeling area, taking stratum L1 as an example, it is determined that the number of volume element layers to be divided is 3;

[0101] Take the first quadrilateral Q1 of the quadrilateral mesh M1 on the top surface of layer L1 1 ({V1 i};i=1,2,6,5) and the first quadrilateral Q2 of the quadrilateral mesh M2 on the bottom of layer L1 1 ({V2 i}; i = 1, 2, 6, 5) as an example, generate two subdivided quadrilaterals Q1 arranged layer by layer from top to bottom p1 ({T1 i};i=1,2,6,5) and Q2 p1 ({T2 i};i=1,2,6,5);

[0102] The first quadrilateral Q1 in the quadrilateral mesh M1 1 and the top subdivided quadrilateral Q1 p1 The shapes are respectively used as the top and bottom surfaces of the generalized tetrahedron to generate a generalized tetrahedron volume element P1 1 ({V1 i},{T1 i};i=1,2,6,5); The first quadrilateral Q2 in the quadrilateral mesh M2 1 and the bottom subdivided quadrilateral Q2 p1 The shapes are respectively used as the top and bottom surfaces of the generalized tetrahedron to generate a generalized tetrahedron volume element P1 3 ({T1 i},{V2 i}; i = 1, 2, 6, 5); Subdivide each pair of two quadrilaterals in the upper and lower adjacent positions Q1 p1 、Q2 p1 The shapes are respectively used as the top and bottom surfaces of the generalized tetrahedron to generate a generalized tetrahedron volume element P1 2 ({T1 i},{T2 i};i=1,2,6,5); Similarly, each generalized tetrahedral voxel in the stratum L1 is generated; So far, all the generalized tetrahedral voxels in the stratum L1 have been generated, thus generating a generalized tetrahedral voxel grid model expressing the stratum L1;

[0103] The generalized tetrahedral element P1 in the generalized tetrahedral element grid model of the stratum L1 1 For example, 20 points out of 24 sample points with attribute information in the modeling area are selected ({S i};i=1,2,3,...,20), according to the three-dimensional coordinates of the above 20 points ({(x i ,y i ,z i )};i=1,2,3,...,20) and porosity property values ​​({a i};i=1,2,3,...,20), and P1 1 The three-dimensional coordinates (x c1 ,y c1 ,z c1 ), using the ordinary Kriging interpolation method, determine P1 1 Porosity property value (a c1 ); Similarly, determine the porosity attribute value of each generalized tetrahedral prism element in the stratum L1;

[0104] At this point, all the generalized tetrahedral volume elements in the stratum L1 and their attribute values ​​have been obtained, thereby generating the attribute volume model of the stratum L1; similarly, the attribute volume models of the strata L2 and L3 are generated, and the number of volume element layers of the stratum L2 is determined to be 2, and the number of volume element layers of the stratum L3 is determined to be 3, thereby constructing a three-dimensional geological attribute volume model containing 3 strata and 2 resolutions.

[0105] The embodiment of the present invention further provides a device for implementing three-dimensional geological modeling; the specific device may include:

[0106] A three-dimensional geological model data expression module is used to express a three-dimensional geological model using a point layer thickness grid data model based on polygonal grids and generalized prismatic elements;

[0107] The 3D geological structure modeling module is used to realize the grid generation and combination of stratum interface by using the point layer thickness interpolation calculation model to construct a 3D geological structure model;

[0108] The three-dimensional geological attribute modeling module is used to realize the generalized prism element segmentation and attribute interpolation of the stratum by adopting the point layer thickness division calculation model, and to construct a three-dimensional geological attribute body model.

[0109] The three-dimensional geological modeling device of this embodiment is used to implement the aforementioned method for implementing three-dimensional geological modeling. Therefore, the specific implementation method of the three-dimensional geological modeling device can be seen in the embodiment part of the method for implementing three-dimensional geological modeling in the previous text. For example, the three-dimensional geological model data expression module, the three-dimensional geological structure modeling module, and the three-dimensional geological attribute modeling module are respectively used to implement steps S101, S102 and S103 in the method for implementing three-dimensional geological modeling. Therefore, its specific implementation method can refer to the description of the corresponding embodiments of each part, which will not be repeated here.

[0110] A specific embodiment of the present invention further provides a device for implementing three-dimensional geological modeling, comprising: a memory for storing a computer program; and a processor for implementing the steps of the method for implementing three-dimensional geological modeling when executing the computer program.

[0111] A specific embodiment of the present invention further provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the steps of the method for implementing three-dimensional geological modeling are implemented.

[0112] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present application may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-LOM, optical storage, etc.) containing computer-usable program code.

[0113] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0114] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0115] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0116] Obviously, the above embodiments are merely examples for the purpose of clear explanation and are not intended to limit the implementation methods. For those skilled in the art, other different forms of changes or modifications can be made based on the above description. It is not necessary and impossible to list all the implementation methods here. The obvious changes or modifications derived therefrom are still within the scope of protection of the present invention.

Claims

1. A method for implementing three-dimensional geological modeling, characterized in that: include: The three-dimensional geological model is expressed by using a point-layer-thickness grid data model based on polygonal grids and generalized prismatic elements; Using the point-by-point layer thickness interpolation calculation model to realize the formation interface grid generation and combination modeling, and construct a three-dimensional geological structure model; and / or, The point-by-point layer thickness division calculation model is used to realize the generalized prism element segmentation and attribute interpolation of the stratum, and to construct a three-dimensional geological attribute body model; The three-dimensional geological model is expressed by using a point layer thickness grid data model based on polygonal grids and generalized prismatic elements, including: Polygonal meshes are used to express the three-dimensional geological structure model: for the m strata stacked from top to bottom in the modeling area, m+1 polygonal meshes distributed from top to bottom are used to express the m+1 stratum interfaces in the strata. The i-th and i+1-th polygonal meshes at any upper and lower positions represent the top and bottom surfaces of the i-th stratum, respectively. All polygon meshes have the same number of vertices, the same number of polygons, the same number of edges, and the same way of vertices forming polygons; The generalized prismatic volume element model is used to express the three-dimensional geological attribute volume model: for any layer L among the m layers i , using n layers of generalized prismatic elements stacked from top to bottom to express the space between the top and bottom surfaces of the stratum, All generalized prisms have the same number of edges, and each generalized prism element has 0 or one or more attribute values. The top surface of each generalized prism in the uppermost voxel layer is related to the stratum L i The polygons in the top polygonal mesh overlap one by one, and the bottom surface of each generalized prism in the lowest voxel layer overlaps with the bottom surface of the stratum L. i The polygons in the bottom polygon mesh overlap one by one. For any two adjacent voxel layers, the bottom surfaces of the generalized prisms in the upper voxel layer coincide with the top surfaces of the generalized prisms in the lower voxel layer. Wherein, m is an integer greater than 0, n is an integer greater than 1, the i-th and i+1-th polygonal meshes dynamically refer to two adjacent polygonal meshes being used, and the i-th stratum dynamically refers to the stratum being used; The method of using the point-level layer thickness interpolation calculation model to realize the formation interface grid generation and combination modeling and constructing the three-dimensional geological structure model includes: For a three-dimensional geological structure model expressed by polygonal meshes, the three-dimensional coordinates of one or more sample points with stratigraphic boundary information on the top surface of the uppermost stratum and the polygonal mesh M1 on the top surface of the uppermost stratum are used. t Any j-th vertex V j The plane coordinates of the vertex V are determined by using the spatial interpolation method. j The vertical coordinates of the polygonal mesh M1 are obtained t The three-dimensional coordinates of each vertex; According to any stratum L i The three-dimensional coordinates of one or more sample points with stratum boundary information on the top surface and the plane coordinate position of the stratum L at the above sample points i The thickness of the formation L i Top polygon mesh M i t Any j-th vertex V j The plane coordinates of the vertex V are determined by using the spatial interpolation method. j The plane coordinate position of the stratum L i The thickness of the polygonal mesh M is obtained i t The plane coordinate position of each vertex in the stratum L i Thickness; Take the vertex V j The vertical coordinate of the vertex V minus or plus j The plane coordinate position of the stratum L i The thickness of the formation L i Bottom polygon mesh M i b The jth vertex V j The vertical coordinates of the polygonal mesh M are obtained i b The three-dimensional coordinates of each vertex in; The three-dimensional coordinates include plane coordinates and vertical coordinates, and the jth vertex dynamically refers to the vertex being used in the polygonal mesh; The point-by-point layer thickness division calculation model is used to realize the generalized prism element segmentation and attribute interpolation of the stratum, and the three-dimensional geological attribute model is constructed, including: For any stratum L in the 3D geological structure model expressed by polygonal meshes i , determine the number of voxel layers n that need to be divided; In the stratum L i Top surface drawing data polygon mesh M i t Any k-th polygon and stratum L in i Bottom drawing data polygonal mesh M i b The spatial position between the k-th polygons in the , generates n-1 subdivided polygons arranged layer by layer from top to bottom; The polygon mesh M i t The kth polygon in the mesh and the subdivided polygon at the top layer are used as the top and bottom surfaces of the generalized prism, respectively, to generate a generalized prism element; the mesh polygon M i b The kth polygon in the subdivided polygon and the subdivided polygon at the bottom layer are respectively used as the top surface and the bottom surface of the generalized prism to generate a generalized prism element; each pair of two polygons in the subdivided polygons that are in the upper and lower adjacent positions are respectively used as the top surface and the bottom surface of the generalized prism to generate a generalized prism element; Wherein, n is an integer greater than 1, and k dynamically refers to the polygons being used in the polygon mesh.

2. The method for implementing three-dimensional geological modeling according to claim 1, characterized in that: The strata are divided according to preset rules, which include geological age, coal rock characteristics, rock mass characteristics, and rock and soil characteristics; the generalized prisms include conventional prisms and variant prisms, and the variant prisms include one or more of the following situations: the top surface is not parallel to the bottom surface, the top surface is not perpendicular to the edge, the bottom surface is not perpendicular to the edge, and the length of the edge is zero.

3. The method for implementing three-dimensional geological modeling according to claim 1, characterized in that: The point-level layer thickness interpolation calculation model is used to realize the formation interface grid generation and combination modeling, and the construction of the three-dimensional geological structure model also includes: For any stratum L in the 3D geological structure model expressed by polygonal meshes i , generating a layer with the formation L i The top polygon mesh has the same vertex and polygon data as the polygon mesh M i t , as the formation L i The top surface drawing data; if the stratum L i If there is a missing area in the stratum, then from the polygonal mesh M i t Delete all polygons whose vertices are located in the missing area of ​​the formation, and then remove the polygon mesh M i t Delete vertices that are not used by any polygon; Generate a layer with L i The bottom polygon mesh has the same vertex and polygon data as the polygon mesh M i b , as the formation L i Bottom surface drawing data; if the stratum L i If there is a missing area in the stratum, then from the polygonal mesh M i b Delete all polygons whose vertices are located in the missing area of ​​the formation, and then remove the polygon mesh M i b Delete vertices that are not used by any polygon; Extraction of formation L i Top and bottom drawing data polygon mesh M i t and M i b The boundary line in the expression layer L is generated according to the boundary line i A polygonal mesh of the surrounding interface; the boundary line is one or more broken lines without repeated vertices formed by connecting the open edges in the polygonal mesh according to the shared vertices between the head and the tail; the open edge is a polygonal edge used by only one polygon in the polygonal mesh.

4. The method for implementing three-dimensional geological modeling according to claim 1, characterized in that: The point layer thickness division calculation model is used to realize the generalized prism element segmentation and attribute interpolation of the stratum. The construction of the three-dimensional geological attribute body model also includes: For any generalized prismatic element P in the generalized prismatic element grid model k , according to the three-dimensional coordinates and attribute values ​​of the sample points with attribute information in the modeling area, and the voxel P k The three-dimensional coordinates of any point in the internal space or on the surface are determined by spatial interpolation calculation method to determine the volume element P k The attribute value it has; The three-dimensional coordinates include plane coordinates and vertical coordinates.

5. A device for implementing three-dimensional geological modeling, characterized in that: include: A three-dimensional geological model data expression module is used to express a three-dimensional geological model using a point layer thickness grid data model based on polygonal grids and generalized prismatic elements; A three-dimensional geological structure modeling module is used to generate and combine stratum interface grids using a point-by-point layer thickness interpolation calculation model to construct a three-dimensional geological structure model; and / or, The 3D geological attribute modeling module is used to implement the generalized prism element segmentation and attribute interpolation of the stratum using the point layer thickness division calculation model to construct a 3D geological attribute body model; The three-dimensional geological model is expressed by using a point layer thickness grid data model based on polygonal grids and generalized prismatic elements, including: Polygonal meshes are used to express the three-dimensional geological structure model: for the m strata stacked from top to bottom in the modeling area, m+1 polygonal meshes distributed from top to bottom are used to express the m+1 stratum interfaces in the strata. The i-th and i+1-th polygonal meshes at any upper and lower positions represent the top and bottom surfaces of the i-th stratum, respectively. All polygon meshes have the same number of vertices, the same number of polygons, the same number of edges, and the same way of vertices forming polygons; The generalized prismatic volume element model is used to express the three-dimensional geological attribute volume model: for any layer L among the m layers i , using n layers of generalized prismatic elements stacked from top to bottom to express the space between the top and bottom surfaces of the stratum, All generalized prisms have the same number of edges, and each generalized prism element has 0 or one or more attribute values. The top surface of each generalized prism in the uppermost voxel layer is related to the stratum L i The polygons in the top polygonal mesh overlap one by one, and the bottom surface of each generalized prism in the lowest voxel layer overlaps with the bottom surface of the stratum L. i The polygons in the bottom polygon mesh overlap one by one. For any two adjacent voxel layers, the bottom surfaces of the generalized prisms in the upper voxel layer coincide with the top surfaces of the generalized prisms in the lower voxel layer. Wherein, m is an integer greater than 0, n is an integer greater than 1, the i-th and i+1-th polygonal meshes dynamically refer to two adjacent polygonal meshes being used, and the i-th stratum dynamically refers to the stratum being used; The method of using the point-level layer thickness interpolation calculation model to realize the formation interface grid generation and combination modeling and constructing the three-dimensional geological structure model includes: For a three-dimensional geological structure model expressed by polygonal meshes, the three-dimensional coordinates of one or more sample points with stratigraphic boundary information on the top surface of the uppermost stratum and the polygonal mesh M1 on the top surface of the uppermost stratum are used. t Any j-th vertex V j The plane coordinates of the vertex V are determined by using the spatial interpolation method. j The vertical coordinates of the polygonal mesh M1 are obtained t The three-dimensional coordinates of each vertex; According to any stratum L i The three-dimensional coordinates of one or more sample points with stratum boundary information on the top surface and the plane coordinate position of the stratum L at the above sample points i The thickness of the formation L i Top polygon mesh M i t Any j-th vertex V j The plane coordinates of the vertex V are determined by using the spatial interpolation method. j The plane coordinate position of the stratum L i The thickness of the polygonal mesh M is obtained i t The plane coordinate position of each vertex in the stratum L i Thickness; Take the vertex V j The vertical coordinate of the vertex V minus or plus j The plane coordinate position of the stratum L i The thickness of the formation L i Bottom polygon mesh M i b The jth vertex V j The vertical coordinates of the polygonal mesh M are obtained i b The three-dimensional coordinates of each vertex in; The three-dimensional coordinates include plane coordinates and vertical coordinates, and the jth vertex dynamically refers to the vertex being used in the polygonal mesh; The point-by-point layer thickness division calculation model is used to realize the generalized prism element segmentation and attribute interpolation of the stratum, and the three-dimensional geological attribute model is constructed, including: For any stratum L in the 3D geological structure model expressed by polygonal meshes i , determine the number of voxel layers n that need to be divided; In the stratum L i Top surface drawing data polygon mesh M i t Any k-th polygon and stratum L in i Bottom drawing data polygonal mesh M i b The spatial position between the k-th polygons in the , generates n-1 subdivided polygons arranged layer by layer from top to bottom; The polygon mesh M i t The kth polygon in the mesh and the subdivided polygon at the top layer are used as the top and bottom surfaces of the generalized prism, respectively, to generate a generalized prism element; the mesh polygon M i b The kth polygon in the subdivided polygon and the subdivided polygon at the bottom layer are respectively used as the top surface and the bottom surface of the generalized prism to generate a generalized prism element; each pair of two polygons in the subdivided polygons that are in the upper and lower adjacent positions are respectively used as the top surface and the bottom surface of the generalized prism to generate a generalized prism element; Wherein, n is an integer greater than 1, and k dynamically refers to the polygons being used in the polygon mesh.

6. A device for implementing three-dimensional geological modeling, characterized in that: include: Memory for storing computer programs; A processor is used to implement the steps of a method for implementing three-dimensional geological modeling as described in any one of claims 1 to 4 when executing the computer program.

7. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, which, when executed by a user, implements the steps of a method for implementing three-dimensional geological modeling as claimed in any one of claims 1 to 4.

Images