CN116863095B - Implementation method and device for large-scale high-precision three-dimensional geological implicit modeling - Google Patents
Implementation method and device for large-scale high-precision three-dimensional geological implicit modeling Download PDFInfo
- Publication number
- CN116863095B CN116863095B CN202311118560.2A CN202311118560A CN116863095B CN 116863095 B CN116863095 B CN 116863095B CN 202311118560 A CN202311118560 A CN 202311118560A CN 116863095 B CN116863095 B CN 116863095B
- Authority
- CN
- China
- Prior art keywords
- sub
- voxel
- lithology
- block
- blocks
- 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
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000000638 solvent extraction Methods 0.000 claims abstract description 14
- 238000012545 processing Methods 0.000 claims description 23
- 238000004590 computer program Methods 0.000 claims description 7
- 230000008859 change Effects 0.000 abstract description 4
- 230000015572 biosynthetic process Effects 0.000 description 16
- 238000005755 formation reaction Methods 0.000 description 16
- 239000003245 coal Substances 0.000 description 13
- 238000010586 diagram Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000005192 partition Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
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
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T1/00—General purpose image data processing
- G06T1/20—Processor architectures; Processor configuration, e.g. pipelining
-
- 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/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Software Systems (AREA)
- Business, Economics & Management (AREA)
- Computer Graphics (AREA)
- Mining & Mineral Resources (AREA)
- Economics (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- Health & Medical Sciences (AREA)
- Agronomy & Crop Science (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Image Generation (AREA)
Abstract
The application provides a method and a device for realizing large-scale high-precision three-dimensional geological implicit modeling, wherein the method comprises the following steps: voxelization is carried out on the space range where the modeling data source is located, and a voxel grid model is obtained; partitioning the voxel grid model to obtain a plurality of first sub-blocks which are partially overlapped; determining a lithology value of each voxel in each first sub-block; determining an overlapping region and a non-overlapping region of adjacent first sub-blocks in the plurality of first sub-blocks; and determining the lithology value of the first voxel according to the mean value of the lithology values of the voxels at corresponding positions in a plurality of first sub-blocks corresponding to the overlapping region of the first voxel in the voxel grid model, and finally obtaining the three-dimensional geological model expressed by the voxel grid. Therefore, smoothness of lithologic value change trend of adjacent first adjacent voxels is enhanced, and accuracy of the constructed three-dimensional geological model is improved.
Description
Technical Field
The application relates to the technical field of mines and computers, in particular to a method and a device for realizing large-scale high-precision three-dimensional geological implicit modeling.
Background
In recent years, with the development of the application of the technology of digital mines, a three-dimensional geological model is an important component and a data base of the construction of the digital mines, and the space analysis calculation and the comprehensive management function of the digital mines are all carried out based on the three-dimensional geological model.
However, due to the large scale characteristics of the geologic volume and the limitations of geologic condition and exploration techniques, it is impossible for geologic survey measurements to obtain complete and regular geologic volume data, so that there is a large difference between a three-dimensional geologic model constructed based on sparse sampled data and a real geologic volume structure. Thus, there is a need for a method of accurately constructing three-dimensional geologic models.
Disclosure of Invention
The application provides a large-scale high-precision three-dimensional geological implicit modeling method and device. The specific scheme is as follows:
in one aspect, the embodiment of the application provides a large-scale high-precision three-dimensional geological implicit modeling method, which comprises the following steps:
voxelization is carried out on the space range where the modeling data source is located, and a voxel grid model is obtained;
partitioning the voxel grid model to obtain a plurality of first sub-blocks, wherein adjacent first sub-blocks are partially overlapped;
determining a lithology value of each voxel in each first sub-block;
determining an overlapping region and a non-overlapping region of adjacent first sub-blocks in the plurality of first sub-blocks;
determining the lithology value of a first voxel in an overlapping region according to the mean value of lithology values of the voxels in corresponding positions in a plurality of first sub-blocks corresponding to the overlapping region of the first voxel in the voxel grid model;
and determining the lithology value of the voxel at the corresponding position in the first sub-block corresponding to the non-overlapping area to which the second voxel belongs as the lithology value of the second voxel in the non-overlapping area, and finally obtaining the three-dimensional geological model expressed by the voxel grid.
In another aspect, the embodiment of the application provides a device for implementing large-scale high-precision three-dimensional geological implicit modeling, which comprises:
the voxelization module is used for voxelization of the space range where the modeling data source is positioned to obtain a voxel grid model;
the partitioning module is used for partitioning the voxel grid model to obtain a plurality of first sub-blocks, wherein adjacent first sub-blocks are partially overlapped;
the first determining module is used for determining the lithology value of each voxel in each first sub-block;
a second determining module, configured to determine an overlapping region and a non-overlapping region of adjacent first sub-blocks in the plurality of first sub-blocks;
the third determining module is used for determining the lithology value of the first voxel according to the mean value of the lithology values of the voxels at the corresponding positions in a plurality of first sub-blocks corresponding to the overlapping area of the first voxel in the voxel grid model;
and the third determining module is used for determining the lithology value of the second voxel in the corresponding first sub-block of the non-overlapping area corresponding to the second voxel in the voxel grid model as the lithology value of the second voxel, and finally obtaining the three-dimensional geological model expressed by the voxel grid.
Another embodiment of the present application provides a computer device, including a processor and a memory;
wherein the processor runs a program corresponding to the executable program code by reading the executable program code stored in the memory for implementing the method of the embodiment of the above aspect.
Another aspect of the present application provides a non-transitory computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements a method as in the above aspect.
In a further aspect, the application proposes a computer program product comprising a computer program which, when executed by a processor, implements the steps of the method of the above embodiment.
Additional aspects and advantages of the application will be set forth in part in the description which follows and, in part, will be obvious from the description, or may be learned by practice of the application.
Drawings
The foregoing and/or additional aspects and advantages of the application will become apparent and readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings, in which:
FIG. 1 (a) is a schematic flow chart of a method for implementing large-scale high-precision three-dimensional geological implicit modeling according to an embodiment of the present application;
FIG. 1 (b) is a schematic flow chart of another implementation method of the large-scale high-precision three-dimensional geological implicit modeling provided by the embodiment of the application;
FIG. 2 is a schematic diagram of a modeling data source provided by the present application;
FIG. 3 is a schematic flow chart of another implementation method of large-scale high-precision three-dimensional geological implicit modeling according to an embodiment of the present application;
FIG. 4 (a) is a schematic diagram of a voxel grid model according to an embodiment of the present application;
FIG. 4 (b) is a schematic diagram of a second sub-block according to an embodiment of the present application;
FIG. 4 (c) is a schematic view of a second sub-block perpendicular to the y-axis according to an embodiment of the present application;
FIG. 4 (d) is a schematic view of a first sub-block perpendicular to the y-axis according to an embodiment of the present application;
FIG. 5 is a schematic flow chart of another implementation method of the large-scale high-precision three-dimensional geological implicit modeling provided by the embodiment of the application;
FIG. 6 is a schematic diagram of all lithology feature points in a first sub-block according to the present application;
fig. 7 is a schematic structural diagram of a device for implementing large-scale high-precision three-dimensional geological implicit modeling according to an embodiment of the present application.
Detailed Description
Embodiments of the present application are described in detail below, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to like or similar elements or elements having like or similar functions throughout. The embodiments described below by referring to the drawings are illustrative and intended to explain the present application and should not be construed as limiting the application.
The method and the device for realizing the large-range high-precision three-dimensional geological implicit modeling are described below with reference to the accompanying drawings.
Fig. 1 (a) is a schematic flow chart of a method for implementing large-scale high-precision three-dimensional geological implicit modeling according to an embodiment of the present application.
As shown in fig. 1 (a), the implementation method of the large-scale high-precision three-dimensional geological implicit modeling comprises the following steps:
and step 101, voxelizing the space range where the modeling data source is located to obtain a voxel grid model.
In the application, the data of the stratum demarcation points of the coal seam, the stratum above the coal seam, the stratum below the coal seam and the like of the working surface can be extracted from the geological data of the working surface in advance. And storing the formation demarcation point data as a modeling data source in the system. The geological data may include, but is not limited to, geological borehole data, geological profile data, and the like.
For example, a discretization may be performed on a geological section line in the face geological section data to obtain a discretized formation demarcation point. The geologic section lines may include, among other things, upper strata top surface boundaries, coal seam top surface boundaries, lower strata bottom surface boundaries, and the like.
In the present application, the modeling data source may include one or more of the following:
a formation demarcation point in the geological borehole;
a string of formation demarcation points formed by discretizing the formation demarcation line;
a slice of formation demarcation point formed by discretizing the fault surface.
For ease of understanding, fig. 2 is a schematic diagram of a modeling data source provided by the present application, where there are 3 strata, namely, a coal seam, an upper strata of the coal seam and a lower strata of the coal seam, but the strata that can be processed by the present application are not limited to 3 strata. The coordinate system XYZ shown in fig. 2 is a geodetic coordinate system. As shown in fig. 2, taking an example of creating a three-dimensional geological model of a coal seam and its top and bottom plates, geological borehole data may include: an upper formation top surface demarcation point, a coal seam top surface demarcation point, a lower formation top surface demarcation point, and a lower formation bottom surface demarcation point; the geologic profile data may include boundaries between adjacent formations: an upper stratum boundary line, a coal seam top boundary line, a lower stratum top boundary line and a lower stratum bottom boundary line. And a coal bed is arranged between the boundary line of the top surface of the coal bed and the boundary line of the top surface of the lower stratum. A three-dimensional geologic model may be created that contains multiple strata including rock formations and/or coal seams in a similar manner.
According to the voxel generation method, the voxel generation can be carried out on the space range where the modeling data source is located according to the obtained voxel generation parameters, and a voxel grid model is obtained. Wherein the voxelization parameters may include, but are not limited to, column number, row number, layer number, and the like. The length, width, height, column number, line number and layer number of the space range determine the dimension of the voxel grid model obtained by voxelization.
It should be noted that, the voxelized parameters may be obtained according to a configuration operation of a user, may be default, or may be obtained by other manners, which is not limited in the present application.
And 102, partitioning the voxel grid model to obtain a plurality of first sub-blocks, wherein adjacent first sub-blocks are partially overlapped.
In the application, in order to improve the modeling efficiency of the model, the voxel grid model can be segmented, so that the voxel grid model is divided into a plurality of first sub-blocks.
As a possible implementation manner, the voxel grid model may be segmented based on a preset size of the first sub-block and a preset sliding step length in each direction of the coordinate system, so as to obtain a plurality of first sub-blocks. The preset sliding step length of each direction of the coordinate system is smaller than or equal to the length of the first sub-block in each direction. The sliding step is the number of voxels per translation. When the preset sliding step length in each direction of the coordinate system is smaller than the length of the first sub-block in each direction, an overlapping part exists between the adjacent first sub-blocks. When the preset sliding step length of each direction of the coordinate system is equal to the length of the first sub-block in each direction, no overlapping part exists between the adjacent first sub-blocks.
Step 103, determining a lithology value of each voxel in each first sub-block.
In the application, the first sub-block and the corresponding voxel of each stratum demarcation point can be determined based on the coordinates of each stratum demarcation point in the modeling data source. And determining the voxel corresponding to the stratum demarcation point in each first sub-block as a lithology characteristic point.
In addition, the modeling data source comprises the lithology value of the stratum demarcation point and the space extension trend information of the lithology value besides the position information of the stratum demarcation point. Therefore, the position information, lithology value and space extension trend information of each stratum demarcation point are determined as the attribute information of lithology characteristic points corresponding to each stratum demarcation point. Wherein the lithology value is used to indicate a formation above a formation demarcation point, such as the numerical sequence numbers of the above-described formation in the formation throughout the modeling range, for example: 0,1,2,3, etc. Wherein, the top surface lithology value of the stratum at the uppermost layer in the whole modeling range is generally 0, the top surface lithology value of the next layer (namely the bottom surface of the layer above the current layer) is generally 1, the lithology value of the next layer is generally 2, and the like. The lithology values are not limited to 0,1,2, and 3, and the lithology value of the lower layer may be set to be larger than the lithology value of the upper layer.
It will be appreciated that the lithology values of the lithology feature points in the first sub-block are known and the lithology values of the other voxels in the first sub-block, except for the lithology feature points, are unknown.
In the application, the closed triangular net corresponding to each first sub-block can be determined according to all lithology characteristic points in each first sub-block. And then, respectively carrying out lithology interpolation processing according to lithology characteristic points in the triangular network corresponding to each first sub-block so as to obtain lithology values of each voxel in the first sub-block. Therefore, lithology interpolation processing is carried out based on lithology characteristic points in the triangular network, the lithology values of unreliable voxels are removed, and the lithology values of the reliable voxels are reserved. Thereby improving the accuracy of determining the lithology value of the voxel and further improving the accuracy of constructing the three-dimensional geological model.
The lithology interpolation processing is performed according to lithology characteristic points in the triangular network corresponding to each first sub-block, so as to obtain lithology values of each voxel in the first sub-block, which may be: for each first sub-block, interpolation can be performed according to the lithology value and the relief trend information of each lithology characteristic point in the triangular network corresponding to the first sub-block, so as to obtain lithology values of other body elements in the triangular network. And determining a non-lithology characteristic point non-lithology value outside the triangular network corresponding to the first sub-block, namely that the non-lithology characteristic point position outside the triangular network corresponding to the first sub-block does not correspond to any stratum. Wherein the non-lithologic feature points are voxels that do not correspond to any formation demarcation point.
In addition, the interpolation may be performed using a spatial interpolation method, which may include, but is not limited to, discrete smooth interpolation, finite difference interpolation, k-space (or kriging) interpolation, and the like.
It can be understood that the lithology interpolation processing is performed based on lithology characteristic points in the triangular net, the lithology value of the unreliable voxel is removed, and the lithology value of the reliable voxel is reserved, so that the accuracy of determining the lithology value of the voxel is improved, and the accuracy of the constructed three-dimensional geological model is further improved.
Optionally, as shown in fig. 1 (b), 1 or more lithology interpolation requests may also be sent to multiple computing devices distributed in the network environment at the same time, where the lithology interpolation requests include one or more first sub-blocks and modeling data sources. After receiving the lithology interpolation request, the computing device may invoke available resources of a central processing unit (Central Processing Unit, CPU) and/or resources of a graphics processor (Graphics Processing Unit, GPU) based on the modeling data source, perform interpolation processing on the first sub-block, and obtain a lithology value of each voxel in the first sub-block. Thereafter, the computing device may generate and send feedback information to the device issuing the lithology interpolation request based on the lithology value of each voxel in the first sub-block. Thus, the lithology value of each voxel in each first sub-block may be determined.
The computing device may run a lithologic interpolation computation service program in real-time to monitor and respond to lithologic interpolation requests. The interpolation process may be: and determining a closed triangular net corresponding to each first sub-block according to all lithology characteristic points in each first sub-block. And then, respectively carrying out lithology interpolation processing according to lithology characteristic points in the triangular network corresponding to each first sub-block so as to obtain lithology values of each voxel in the first sub-block.
Step 104, determining an overlapping region and a non-overlapping region of adjacent first sub-blocks in the plurality of first sub-blocks.
In the application, the positions of all the voxels in the adjacent first sub-blocks can be matched, and the overlapping area and the non-overlapping area of the adjacent first sub-blocks can be determined.
Step 105, determining the lithology value of the first voxel in the overlapping region according to the mean value of the lithology values of the voxels in the corresponding positions in the plurality of first sub-blocks corresponding to the overlapping region of the first voxel.
For example, a certain overlapping area is generated by partially overlapping the first sub-block 1 and the first sub-block 2, and the lithology value of a certain first voxel in the overlapping area may be a mean value of the lithology value of the voxel corresponding to the first voxel in the first sub-block 1 and the lithology value of the voxel corresponding to the first voxel in the first sub-block 2.
And 106, determining the lithology value of the voxel at the corresponding position in the first sub-block corresponding to the non-overlapping area to which the second voxel belongs as the lithology value of the second voxel in the non-overlapping area, and finally obtaining the three-dimensional geological model expressed by the voxel grid.
Optionally, in any one of the first sub-blocks corresponding to the first sub-blocks in the overlapping area to which the first voxel belongs, under the condition that the lithology value of the corresponding position voxel of the first voxel is unknown, determining the lithology value of the first voxel according to the average value of lithology values of the corresponding position voxels of the first voxel in other first sub-blocks except for any one of the first sub-blocks corresponding to the overlapping area to which the first voxel belongs.
For example, if a certain overlapping area is generated by partially overlapping the first sub-block 1, the first sub-block 2 and the first sub-block 3, and the lithology value of the voxel corresponding to the position of a certain first voxel in the overlapping area in the first sub-block 1 is unknown, the mean value of the lithology value of the voxel corresponding to the position of the first voxel in the first sub-block 2 and the lithology value of the voxel corresponding to the position of the first voxel in the first sub-block 3 can be determined as the lithology value of the first voxel.
According to the method, a space range where a modeling data source is located can be subjected to voxelization to obtain a voxel grid model, the voxel grid model is segmented to obtain a plurality of first sub-blocks with overlapping portions, then the lithology value of each voxel in each first sub-block is determined, the overlapping area and the non-overlapping area of the adjacent first sub-blocks in the plurality of first sub-blocks are determined, then the lithology value of a first voxel is determined according to the mean value of the lithology values of the voxels in the corresponding positions in the plurality of first sub-blocks corresponding to the overlapping area of the first voxel in the voxel grid model, and the lithology value of the voxel in the corresponding positions is determined as the lithology value of the second voxel in the first sub-block corresponding to the non-overlapping area of the second voxel in the voxel grid model, so that the three-dimensional geological model expressed by the voxel grid is finally obtained. Therefore, based on the lithology value corresponding to the redundant part voxel of the adjacent first sub-block, the lithology value of the voxel of the redundant part at the corresponding position of the voxel grid model is determined, the smoothness of the lithology value change trend of the adjacent first adjacent voxel is enhanced, and the accuracy of the constructed three-dimensional geological model is improved.
FIG. 3 is a schematic flow chart of another implementation method of large-scale high-precision three-dimensional geological implicit modeling according to an embodiment of the present application.
As shown in FIG. 3, the implementation method of the large-range high-precision three-dimensional geological implicit modeling comprises the following steps:
step 301, voxelizing the space range where the modeling data source is located, and obtaining a voxel grid model.
In the present application, the specific implementation process of step 301 may be referred to the detailed description of any embodiment of the present application, and will not be repeated here.
And step 302, partitioning the voxel grid model to obtain a plurality of second sub-blocks without overlapping parts.
In the present application, a specific implementation process of partitioning a voxel grid model to obtain a plurality of second sub-blocks without overlapping portions may refer to a detailed description of partitioning a voxel grid model to obtain a plurality of first sub-blocks without overlapping portions in any embodiment of the present application, which is not described herein.
Step 303, obtaining the expansion direction and the expansion number corresponding to the expansion direction corresponding to each second sub-block.
In the application, the expansion direction corresponding to each second sub-block and the expansion quantity corresponding to the expansion direction can be preset in the system. The expansion direction corresponding to each second sub-block may be one or more. The expansion direction and the expansion number corresponding to the expansion direction of each second sub-block may be the same or different.
Step 304, adding an expansion number of voxels to the boundary of each second sub-block in the expansion direction to obtain a first sub-block corresponding to each second sub-block.
Fig. 4 (a) is a schematic diagram of a voxel grid model in which one cube is one voxel. Fig. 4 (b) is a schematic diagram of a second sub-block for partition determination of the voxel grid model of fig. 4 (a). Fig. 4 (c) is a schematic view of the second sub-block of fig. 4 (b) perpendicular to the y-axis. A schematic view of a plane perpendicular to the y-axis of a first sub-block obtained by expanding each second sub-block by one voxel in the positive x-axis direction is shown in fig. 4 (d). The voxels of the same shaded portion in 4 (d) are overlapping portion voxels.
It can be understood that, in the first sub-block obtained after expanding the second sub-block, an overlapping portion exists between adjacent first sub-blocks, so that the first sub-block has partial redundancy. And determining the lithology value of the redundant partial voxel at the corresponding position of the voxel grid model based on lithology values corresponding to the redundant partial voxels of the plurality of first sub-blocks, thereby being beneficial to improving the accuracy of the lithology value of the voxel at the boundary of the sub-blocks and enhancing the smoothness of the lithology value change trend between adjacent sub-blocks, and further improving the accuracy of the constructed three-dimensional geological model.
And 305, determining a closed triangular network corresponding to each first sub-block according to all lithology characteristic points in each first sub-block.
And 306, performing lithologic interpolation processing according to lithologic feature points in the triangular network corresponding to each first sub-block to obtain lithologic values of each voxel in the first sub-block.
In the present application, the specific implementation process of step 305 to step 306 can be referred to the detailed description of any embodiment of the present application, and will not be repeated here.
Step 307, determining a three-dimensional geological model expressed by the voxel grid based on the lithology values of the voxels of each first sub-block.
In the application, under the condition that no overlapping part exists between the adjacent first sub-blocks, the lithology values of the voxels in the plurality of first sub-blocks can be assigned to the voxels at the corresponding positions in the voxel grid model to obtain the three-dimensional geological model expressed by the voxel grid.
Optionally, when there is an overlapping portion between the adjacent first sub-blocks, an overlapping area and a non-overlapping area of the adjacent first sub-blocks in the plurality of first sub-blocks may be determined, and for a certain first voxel located in the overlapping area in the voxel grid model, a mean value of lithology values of the corresponding voxel in the corresponding location in the plurality of first sub-blocks corresponding to the overlapping area to which the first voxel belongs may be determined as the lithology value of the first voxel. And determining the lithology value of the second voxel at the corresponding position in the first sub-block corresponding to the non-overlapping area to which the second voxel belongs as the lithology value of the second voxel aiming at the second voxel in the non-overlapping area in the voxel grid model.
Optionally, in any one of the first sub-blocks corresponding to the first sub-blocks in the overlapping area to which the first voxel belongs, under the condition that the lithology value of the corresponding position voxel of the first voxel is unknown, determining the lithology value of the first voxel according to the average value of lithology values of the corresponding position voxels of the first voxel in other first sub-blocks except for any one of the first sub-blocks corresponding to the overlapping area to which the first voxel belongs.
In the application, after voxel forming is carried out on a space range where a modeling data source is located to obtain a voxel grid model, the voxel grid model can be segmented to obtain a plurality of second sub-blocks without overlapped parts, the expansion direction corresponding to each second sub-block and the expansion quantity corresponding to the expansion direction are obtained, so that the expansion quantity of voxels in the expansion direction is increased for the boundary of each second sub-block to obtain first sub-blocks corresponding to each second sub-block, then, a closed triangular net corresponding to each first sub-block is determined according to all lithology characteristic points in each first sub-block, lithology interpolation processing is carried out according to lithology characteristic points in the triangular net corresponding to each first sub-block respectively, so as to obtain the lithology value of each voxel in the first sub-block, and then, the three-dimensional geological model expressed by the voxel grid is determined based on the lithology value of the voxels of each first sub-block. Thereby improving the accuracy of determining the lithology value of the voxel and further improving the accuracy of constructing the three-dimensional geological model.
FIG. 5 is a schematic flow chart of another implementation method of large-scale high-precision three-dimensional geological implicit modeling according to an embodiment of the present application.
As shown in FIG. 5, the implementation method of the large-range high-precision three-dimensional geological implicit modeling comprises the following steps:
and step 501, voxelizing the space range where the modeling data source is located to obtain a voxel grid model.
And step 502, partitioning the voxel grid model to obtain a plurality of first sub-blocks, wherein adjacent first sub-blocks are partially overlapped.
In the present application, the specific implementation process of step 501 to step 502 can be referred to the detailed description of any embodiment of the present application, and will not be repeated here.
At step 503, all lithology feature points in the first sub-block are projected to the target coordinate plane.
The target coordinate plane may be a plane in which an x-axis and a y-axis are located, or a plane in which an x-axis and a z-axis are located, or a plane in which a y-axis and a z-axis are located, which is not limited in the present disclosure.
Step 504, determining a first plane convex hull according to projection points of all lithology characteristic points in the target coordinate plane.
In the application, the projection points of all lithology characteristic points in the target coordinate plane can be processed based on a preset arbitrary convex hull algorithm and the like to generate the first plane convex hull.
It will be appreciated that the first planar convex hull is polygonal and that the projected points of all lithology feature points in the target coordinate plane are located inside or on the sides of the first planar convex hull.
And 505, determining the minimum coordinate value and the maximum coordinate value of all lithology characteristic points on a target coordinate axis, wherein the target coordinate axis is perpendicular to the target coordinate plane.
It will be appreciated that when the target coordinate plane is the plane in which the x-axis and y-axis lie, the target coordinate axis is the z-axis. When the target coordinate plane is the plane in which the x-axis and the z-axis are located, the target coordinate axis is the y-axis. When the target coordinate plane is the plane in which the y-axis and the z-axis are located, the target coordinate axis is the x-axis.
As shown in fig. 6, fig. 6 is a schematic view of all lithology feature points in the first sub-block, each cube corresponding to one lithology feature point. Assuming that the target coordinate plane is the plane in which the x-axis and the y-axis are located, the coordinate of the H1 lithology characteristic point on the z-axis of the target coordinate axis is the minimum coordinate value, and the coordinate of the H2 lithology characteristic point on the z-axis of the target coordinate axis is the maximum coordinate value.
Step 506, generating a second planar convex hull and a third planar convex hull at the minimum coordinate value and the maximum coordinate value respectively according to the first planar convex hull.
In the application, the first plane convex hull can be translated to the minimum coordinate value to generate the second plane convex hull. And translating the first plane convex hull to the maximum coordinate value to generate a third plane convex hull.
And 507, generating a closed triangular net according to the second plane convex hull and the third plane convex hull.
In the present application, as shown in fig. 6, the vertices of the second planar convex hull and the vertices of the third planar convex hull may be triangulated to generate a self-enclosed triangulation network. Thereby improving the accuracy of the generated triangle net.
Alternatively, the three-dimensional external boundary point may be determined from all lithology feature points in the first sub-block based on a preset convex hull algorithm. And then triangulating the external boundary points to generate a closed triangular net (namely, a three-dimensional convex hull containing all lithology characteristic points in the first sub-block). Thereby improving the accuracy of the generated triangle net.
And step 508, performing lithology interpolation processing according to lithology characteristic points in the triangular network corresponding to each first sub-block so as to obtain lithology values of each voxel in the first sub-block.
Step 509, determining a three-dimensional geological model expressed by the voxel grid based on the lithology values of the voxels of each first sub-block.
In the present application, the specific implementation process of step 508 to step 509 can be referred to the detailed description of any embodiment of the present application, and will not be repeated here.
In the application, all lithology characteristic points in a first sub-block can be projected to a target coordinate plane, a first plane convex hull is determined according to the projection points of all lithology characteristic points in the target coordinate plane, then the minimum coordinate value and the maximum coordinate value of all lithology characteristic points on the target coordinate axis vertical to the target coordinate plane are determined, so that a second plane convex hull and a third plane convex hull are respectively generated at the positions of the minimum coordinate value and the maximum coordinate value according to the first plane convex hull, and then a closed triangular net is generated according to the second plane convex hull and the third plane convex hull. Thereby improving the accuracy of generating the triangle network and further improving the accuracy of the three-dimensional address model.
In order to achieve the above embodiment, the embodiment of the application also provides a device for achieving large-scale high-precision three-dimensional geological implicit modeling. Fig. 7 is a schematic structural diagram of a device for implementing large-scale high-precision three-dimensional geological implicit modeling according to an embodiment of the present application.
As shown in fig. 7, the apparatus 700 includes:
the voxelization module 710 is configured to voxelize a spatial range where the modeling data source is located to obtain a voxel grid model;
the partitioning module 720 is configured to partition the voxel grid model to obtain a plurality of first sub-blocks, where adjacent first sub-blocks are partially overlapped;
a first determining module 730, configured to determine a lithology value of each voxel in each first sub-block;
a second determining module 740, configured to determine an overlapping region and a non-overlapping region of adjacent first sub-blocks in the plurality of first sub-blocks;
a third determining module 750, configured to determine, for a first voxel located in the overlapping region in the voxel grid model, a lithology value of the first voxel according to a mean value of lithology values of the voxels at corresponding positions in a plurality of first sub-blocks corresponding to the overlapping region to which the first voxel belongs;
the third determining module 750 is configured to determine, for a second voxel located in the non-overlapping region in the voxel grid model, a lithology value of the second voxel at a corresponding position in a first sub-block corresponding to the non-overlapping region to which the second voxel belongs, as a lithology value of the second voxel, and finally obtain a three-dimensional geological model expressed by the voxel grid.
In one possible implementation manner of the embodiment of the present application, the blocking module 720 is configured to:
partitioning the voxel grid model to obtain a plurality of second sub-blocks without overlapping parts;
acquiring the expansion direction corresponding to each second sub-block and the expansion quantity corresponding to the expansion direction;
and adding an expansion number of voxels to the boundary of each second sub-block in the expansion direction to acquire a first sub-block corresponding to each second sub-block.
In one possible implementation manner of the embodiment of the present application, the first determining module 730 is configured to:
determining a closed triangular net corresponding to each first sub-block according to all lithology characteristic points in each first sub-block;
and respectively carrying out lithology interpolation processing according to lithology characteristic points in the triangular network corresponding to each first sub-block so as to obtain lithology values of each voxel in the first sub-block.
In one possible implementation manner of the embodiment of the present application, the first determining module 730 is configured to:
projecting all lithology characteristic points in the first sub-block to a target coordinate plane;
determining a first plane convex hull according to projection points of all lithology characteristic points in a target coordinate plane;
determining the minimum coordinate value and the maximum coordinate value of all lithology characteristic points on a target coordinate axis, wherein the target coordinate axis is perpendicular to a target coordinate plane;
generating a second plane convex hull and a third plane convex hull at the positions of the minimum coordinate value and the maximum coordinate value respectively according to the first plane convex hull;
and generating a closed triangular net according to the second plane convex hull and the third plane convex hull.
In one possible implementation manner of the embodiment of the present application, the first determining module 730 is configured to:
determining three-dimensional external boundary points from all lithology characteristic points in the first sub-block;
triangulation is performed on the external boundary points to generate a self-closing triangulation network.
In one possible implementation manner of the embodiment of the present application, the first determining module 730 is configured to:
1 or more lithologic interpolation requests are respectively sent to a plurality of computing devices distributed in a network environment at one time, the computing devices call respective CPU (central processing unit) resources and/or GPU (graphics processing unit) resources, interpolation processing is carried out on each first sub-block, and feedback information is sent to the device sending the lithologic interpolation requests after the processing is finished;
wherein the lithology interpolation request comprises first sub-blocks and modeling data sources; the lithology interpolation calculation service program is operated in the calculation equipment in real time, and lithology interpolation requests are monitored and responded; the feedback information contains lithology values of each voxel in each first sub-block.
In one possible implementation manner of the embodiment of the present application, the second determining module 750 is configured to:
and determining the lithology value of the first voxel according to the average value of lithology values of the corresponding position voxels in the first sub-blocks except any first sub-block corresponding to the overlapping region of the first voxel under the condition that the lithology value of the corresponding position voxel of the first voxel is unknown in any first sub-block of the plurality of first sub-blocks corresponding to the overlapping region of the first voxel.
The explanation of the embodiment of the implementation method of the implicit modeling of the three-dimensional geology with high precision in a large range is also applicable to the implementation device of the implicit modeling of the three-dimensional geology with high precision in a large range in this embodiment, so that the explanation is omitted here.
According to the method, a space range where a modeling data source is located can be subjected to voxelization to obtain a voxel grid model, the voxel grid model is segmented to obtain a plurality of first sub-blocks with overlapping portions, then the lithology value of each voxel in each first sub-block is determined, the overlapping area and the non-overlapping area of the adjacent first sub-blocks in the plurality of first sub-blocks are determined, then the lithology value of a first voxel is determined according to the mean value of the lithology values of the voxels in the corresponding positions in the plurality of first sub-blocks corresponding to the overlapping area of the first voxel in the voxel grid model, and the lithology value of the voxel in the corresponding positions is determined as the lithology value of the second voxel in the first sub-block corresponding to the non-overlapping area of the second voxel in the voxel grid model, so that the three-dimensional geological model expressed by the voxel grid is finally obtained. Therefore, based on the lithology value corresponding to the redundant part voxel of the adjacent first sub-block, the lithology value of the voxel of the redundant part at the corresponding position of the voxel grid model is determined, the smoothness of the lithology value change trend of the adjacent first adjacent voxel is enhanced, and the accuracy of the constructed three-dimensional geological model is improved.
In order to implement the above embodiments, the embodiments of the present application further provide a computer device, including a processor and a memory;
wherein the processor runs a program corresponding to the executable program code by reading the executable program code stored in the memory, for implementing the method for implementing the extensive high-precision three-dimensional geological implicit modeling as described in the above embodiments.
In order to implement the above embodiment, the embodiment of the present application further proposes a non-transitory computer-readable storage medium, on which a computer program is stored, which when executed by a processor implements the method for implementing high-precision three-dimensional geological implicit modeling in a large range as described in the above embodiment.
In order to implement the above embodiments, the embodiments of the present application also propose a computer program product comprising a computer program which, when being executed by a processor, implements the steps of the method described in the above embodiments.
In the description of this specification, the terms "first," "second," and the like are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defining "a first" or "a second" may explicitly or implicitly include at least one such feature. In the description of the present application, the meaning of "plurality" means at least two, for example, two, three, etc., unless specifically defined otherwise.
While embodiments of the present application have been shown and described above, it will be understood that the above embodiments are illustrative and not to be construed as limiting the application, and that variations, modifications, alternatives and variations may be made to the above embodiments by one of ordinary skill in the art within the scope of the application.
Claims (7)
1. The implementation method of the three-dimensional geological implicit modeling is characterized by comprising the following steps:
voxelization is carried out on the space range where the modeling data source is located, and a voxel grid model is obtained;
partitioning the voxel grid model to obtain a plurality of first sub-blocks, wherein adjacent first sub-blocks are partially overlapped;
determining a closed triangular net corresponding to each first sub-block according to all lithology characteristic points in each first sub-block; respectively carrying out lithology interpolation processing according to lithology characteristic points in the triangular network corresponding to each first sub-block so as to obtain lithology values of each voxel in the first sub-block; the determining the closed triangular net corresponding to each first sub-block according to all lithology characteristic points in each first sub-block comprises the following steps: projecting all lithology characteristic points in the first sub-block to a target coordinate plane; determining a first plane convex hull according to projection points of all lithology characteristic points in the target coordinate plane; determining the minimum coordinate value and the maximum coordinate value of all lithology characteristic points on a target coordinate axis, wherein the target coordinate axis is perpendicular to the target coordinate plane; generating a second plane convex hull and a third plane convex hull at the positions of the minimum coordinate value and the maximum coordinate value respectively according to the first plane convex hull; generating a closed triangular net according to the second plane convex hull and the third plane convex hull; or determining three-dimensional external boundary points from all lithology characteristic points in the first sub-block; triangulating the external boundary points to generate a self-closed triangulation network;
determining an overlapping region and a non-overlapping region of adjacent first sub-blocks in the plurality of first sub-blocks;
determining lithology values of first voxels positioned in the overlapping region in the voxel grid model according to lithology value average values of the corresponding position voxels in a plurality of first sub-blocks corresponding to the overlapping region to which the first voxels belong;
and determining the lithology value of the corresponding position voxel in a first sub-block corresponding to the non-overlapping area to which the second voxel belongs as the lithology value of the second voxel aiming at the second voxel in the non-overlapping area in the voxel grid model, and finally obtaining the three-dimensional geological model expressed by the voxel grid.
2. The method of claim 1, wherein the partitioning the voxel grid model to obtain a plurality of first sub-blocks comprises:
partitioning the voxel grid model to obtain a plurality of second sub-blocks without overlapping parts;
acquiring the expansion direction corresponding to each second sub-block and the expansion quantity corresponding to the expansion direction;
and adding the volume elements of the expansion quantity to the boundary of each second sub-block in the expansion direction so as to obtain a first sub-block corresponding to each second sub-block.
3. The method of claim 1, wherein said determining the lithology value of each voxel in each of said first sub-blocks comprises:
1 or more lithologic interpolation requests are respectively sent to a plurality of computing devices distributed in a network environment at one time, the computing devices call respective CPU (central processing unit) resources and/or GPU (graphics processing unit) resources, interpolation processing is carried out on each first sub-block, and feedback information is sent to the device sending the lithologic interpolation requests after the processing is finished;
wherein the lithology interpolation request comprises the first sub-block and the modeling data source; the computing equipment runs a lithology interpolation computing service program in real time, and monitors and responds to lithology interpolation requests; the feedback information comprises lithology values of each voxel in each first sub-block.
4. The method of claim 1, wherein determining the lithology value of the first voxel from the mean of lithology values of the corresponding location voxels in the plurality of first sub-blocks corresponding to the overlapping region to which the first voxel belongs comprises:
and determining the lithology value of the first voxel according to the average value of lithology values of the corresponding position voxels in other first sub-blocks except the first sub-block corresponding to the overlapping region of the first voxel under the condition that the lithology value of the corresponding position voxel of the first voxel is unknown in any one of the first sub-blocks of the plurality of first sub-blocks corresponding to the overlapping region of the first voxel.
5. A three-dimensional geological implicit modeling implementation device, comprising:
the voxelization module is used for voxelization of the space range where the modeling data source is positioned to obtain a voxel grid model;
the block dividing module is used for dividing the voxel grid model into a plurality of first sub-blocks, wherein the adjacent first sub-blocks are partially overlapped;
the first determining module is used for determining a closed triangular net corresponding to each first sub-block according to all lithology characteristic points in each first sub-block; respectively carrying out lithology interpolation processing according to lithology characteristic points in the triangular network corresponding to each first sub-block so as to obtain lithology values of each voxel in the first sub-block; the determining the closed triangular net corresponding to each first sub-block according to all lithology characteristic points in each first sub-block comprises the following steps: projecting all lithology characteristic points in the first sub-block to a target coordinate plane; determining a first plane convex hull according to projection points of all lithology characteristic points in the target coordinate plane; determining the minimum coordinate value and the maximum coordinate value of all lithology characteristic points on a target coordinate axis, wherein the target coordinate axis is perpendicular to the target coordinate plane; generating a second plane convex hull and a third plane convex hull at the positions of the minimum coordinate value and the maximum coordinate value respectively according to the first plane convex hull; generating a closed triangular net according to the second plane convex hull and the third plane convex hull; or determining three-dimensional external boundary points from all lithology characteristic points in the first sub-block; triangulating the external boundary points to generate a self-closed triangulation network;
a second determining module, configured to determine an overlapping region and a non-overlapping region of adjacent first sub-blocks in the plurality of first sub-blocks;
the third determining module is used for determining lithology value of the first voxel in the overlapping region according to lithology value mean value of the corresponding position voxel in a plurality of first sub-blocks corresponding to the overlapping region of the first voxel in the voxel grid model;
the third determining module is configured to determine, for a second voxel located in the non-overlapping region in the voxel grid model, a lithology value of the second voxel corresponding to a position voxel in a first sub-block corresponding to the non-overlapping region to which the second voxel belongs, as a lithology value of the second voxel, and finally obtain a three-dimensional geological model expressed by the voxel grid.
6. A computer device comprising a processor and a memory;
wherein the processor runs a program corresponding to executable program code stored in the memory by reading the executable program code for implementing the method according to any one of claims 1-4.
7. A non-transitory computer readable storage medium, on which a computer program is stored, characterized in that the program, when executed by a processor, implements the method according to any one of claims 1-4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311118560.2A CN116863095B (en) | 2023-09-01 | 2023-09-01 | Implementation method and device for large-scale high-precision three-dimensional geological implicit modeling |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311118560.2A CN116863095B (en) | 2023-09-01 | 2023-09-01 | Implementation method and device for large-scale high-precision three-dimensional geological implicit modeling |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116863095A CN116863095A (en) | 2023-10-10 |
| CN116863095B true CN116863095B (en) | 2023-11-07 |
Family
ID=88230762
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311118560.2A Active CN116863095B (en) | 2023-09-01 | 2023-09-01 | Implementation method and device for large-scale high-precision three-dimensional geological implicit modeling |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116863095B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113689559A (en) * | 2021-08-12 | 2021-11-23 | 贵州省地质调查院(贵州省地质矿产勘查开发局地质科学研究所) | Grid-shaped graph cutting geological section wiring method for improving geological three-dimensional model content expression precision |
| CN114255188A (en) * | 2021-12-23 | 2022-03-29 | 中国矿业大学(北京) | Boundary linear smoothing method and device for drawing three-dimensional geological lithologic grid model |
| CN116433864A (en) * | 2023-06-14 | 2023-07-14 | 煤炭科学研究总院有限公司 | Coal face three-dimensional geological implicit modeling method and device |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7636610B2 (en) * | 2006-07-19 | 2009-12-22 | Envisiontec Gmbh | Method and device for producing a three-dimensional object, and computer and data carrier useful therefor |
-
2023
- 2023-09-01 CN CN202311118560.2A patent/CN116863095B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113689559A (en) * | 2021-08-12 | 2021-11-23 | 贵州省地质调查院(贵州省地质矿产勘查开发局地质科学研究所) | Grid-shaped graph cutting geological section wiring method for improving geological three-dimensional model content expression precision |
| CN114255188A (en) * | 2021-12-23 | 2022-03-29 | 中国矿业大学(北京) | Boundary linear smoothing method and device for drawing three-dimensional geological lithologic grid model |
| CN116433864A (en) * | 2023-06-14 | 2023-07-14 | 煤炭科学研究总院有限公司 | Coal face three-dimensional geological implicit modeling method and device |
Non-Patent Citations (2)
| Title |
|---|
| 基于HDF5的煤矿地质三维层叠网格模型分布式存储研究;郭军;《工矿自动化》;第49卷(第1期);第153-161页 * |
| 煤矿区域措施钻孔三维反演系统开发;岳俊;《煤矿机械》;第44卷(第6期);第209-212页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116863095A (en) | 2023-10-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| USRE49507E1 (en) | Faulted geological structures having unconformities | |
| US8463586B2 (en) | Machine, program product, and computer-implemented method to simulate reservoirs as 2.5D unstructured grids | |
| CN110529077B (en) | Modeling cross-faults and complex wellbores in reservoir simulation | |
| US8965745B2 (en) | Grid from depositional space | |
| US10571601B2 (en) | Geologic fault representation for flow simulation | |
| US9792726B2 (en) | Determining and displaying reserve estimates for a reservoir | |
| EP2869096A2 (en) | Systems and methods of multi-scale meshing for geologic time modeling | |
| CN106981093A (en) | A kind of three-dimensional formation parallel modeling method of subregion constraint coupling | |
| CN104884974A (en) | Systems and methods for 3d seismic data depth conversion utilizing artificial neural networks | |
| US10387583B2 (en) | Rotations from gradient directions | |
| KR20180056086A (en) | Constructing Method of Three Dimension Visualization Model for Tunnel Information and System of The Same | |
| CN102057368A (en) | Distribution of properties in a 3D volumetric model using a maximum continuity field | |
| CN116433864B (en) | Coal face three-dimensional geological implicit modeling method and device | |
| KR20180069547A (en) | Constructing Method of Three Dimension Model Using Line Information Based on Two Dimension Drawing and System of The Same | |
| WO2016001697A1 (en) | Systems and methods for geologic surface reconstruction using implicit functions | |
| EP3455459A2 (en) | Three-dimensional fracture abundance evaluation of subsurface formations | |
| AU2011271175B2 (en) | Machine, program product, and computer-implemented method to simulate reservoirs as 2.5D unstructured grids | |
| CN116863095B (en) | Implementation method and device for large-scale high-precision three-dimensional geological implicit modeling |
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 |