CN118505921A - Model three-dimensional visualization method and system for geological archives - Google Patents
Model three-dimensional visualization method and system for geological archives Download PDFInfo
- Publication number
- CN118505921A CN118505921A CN202410955946.7A CN202410955946A CN118505921A CN 118505921 A CN118505921 A CN 118505921A CN 202410955946 A CN202410955946 A CN 202410955946A CN 118505921 A CN118505921 A CN 118505921A
- Authority
- CN
- China
- Prior art keywords
- model
- contour line
- geological
- dimensional
- dimensional model
- 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.)
- Granted
Links
- 238000007794 visualization technique Methods 0.000 title claims abstract description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 68
- 238000005553 drilling Methods 0.000 claims abstract description 41
- 238000013519 translation Methods 0.000 claims abstract description 32
- 238000000034 method Methods 0.000 claims abstract description 26
- 238000012800 visualization Methods 0.000 claims abstract description 14
- 230000000007 visual effect Effects 0.000 claims abstract description 7
- 239000013598 vector Substances 0.000 claims description 17
- 239000011159 matrix material Substances 0.000 claims description 11
- 230000015572 biosynthetic process Effects 0.000 claims description 10
- 238000009499 grossing Methods 0.000 claims description 8
- 238000007781 pre-processing Methods 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 4
- 239000006185 dispersion Substances 0.000 claims description 4
- 238000006073 displacement reaction Methods 0.000 claims description 4
- 238000013507 mapping Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 239000002245 particle Substances 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 2
- 238000005755 formation reaction Methods 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 239000013049 sediment Substances 0.000 description 2
- 238000012876 topography Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 239000003643 water by type Substances 0.000 description 1
Landscapes
- Processing Or Creating Images (AREA)
Abstract
The invention discloses a model three-dimensional visualization method and a model three-dimensional visualization system for geological files, and relates to the technical field of geological three-dimensional modeling, wherein the method comprises the steps of obtaining drilling data of a region to be modeled, and generating geological profile data corresponding to topographic features according to a sub-drilling dataset; based on the sub-drilling data and the geological profile data, generating a local geological model corresponding to each topographic feature, setting a water area three-dimensional model as a reference model, and extracting a first boundary contour line of the water area three-dimensional model and a second contour line of the intersection of the outer surface of the plain three-dimensional model and a fracture surface; matching the first contour line and the second contour line by a curve matching method to obtain the corresponding relation between the first contour line and the second contour line; and obtaining rotation and translation parameters of the plain three-dimensional model based on the corresponding relation between the first contour line and the second contour line, and splicing the water area three-dimensional model and the plain three-dimensional model through the rotation and translation parameters to generate a three-dimensional visual model of the area to be modeled.
Description
Technical Field
The invention relates to the technical field of geological three-dimensional modeling, in particular to a model three-dimensional visualization method and a model three-dimensional visualization system for geological files.
Background
With the continued development of geological exploration technology, a large amount of geological data is acquired and stored. However, the conventional two-dimensional graphic representation cannot fully display the three-dimensional spatial characteristics of the geological data, so that the interpretation of the geological information has a great limitation. Therefore, three-dimensional visualization of geological data is needed to be realized, accuracy and efficiency of geological information interpretation are improved, three-dimensional geological modeling mainly refers to a process of generating a three-dimensional model by using a computer image forming technology, and the constructed three-dimensional geological model can enable geological staff to clearly observe ore distribution, reserves and structures of underground geological bodies, so that subsequent research can be facilitated.
In the prior art, a plurality of local geologic models corresponding to the topographic features are generated in real time through drilling data, and then the local geologic models are spliced to generate a target three-dimensional geologic model.
However, in the actual modeling process, it may be difficult to accurately reflect the true morphology and boundaries of the formation due to the complexity and uncertainty of the geological phenomenon. This may result in discontinuities or inconsistencies in geometry, topology or spatial locations between the local geologic models during the stitching process.
Disclosure of Invention
In order to solve the problems in the prior art, the invention provides a model three-dimensional visualization method and a system for geological files.
The technical scheme adopted by the invention is as follows:
The first aspect of the application provides a model three-dimensional visualization method of a geological archive, which comprises the following steps: obtaining drilling data of an area to be molded, wherein the area to be molded comprises a plurality of topographic features, the drilling data comprises a plurality of sub-drilling data sets, and each sub-drilling data set corresponds to each topographic feature; generating geological profile data corresponding to the topographic features according to the sub-drilling data set; generating a local geological model corresponding to each of the topographical features based on the sub-borehole data and the geological profile data, wherein the topographical features comprise a water area and a plains, and the local geological model comprises a water area three-dimensional model and a plains three-dimensional model; setting a water area three-dimensional model as a reference model, setting a surface of the water area three-dimensional model intersecting with a plain three-dimensional model as a fracture surface, and extracting a first boundary contour line of the water area three-dimensional model and a second contour line of the plain three-dimensional model intersecting with the fracture surface; matching the first contour line and the second contour line by a curve matching method to obtain the corresponding relation between the first contour line and the second contour line; obtaining rotation and translation parameters of the plain three-dimensional model based on the corresponding relation between the first contour line and the second contour line; and (3) splicing the water area three-dimensional model and the plain three-dimensional model through rotation and translation parameters, and generating a three-dimensional visual model of the area to be modeled.
Preferably, generating geological profile data corresponding to the topographical features from the sub-borehole dataset includes:
Screening feature data related to specific topographic features from the sub-borehole data set, wherein the feature data comprises coordinates, depth and formation information of the borehole;
and determining the position and the direction of the geological profile to be generated according to the topographic features and modeling requirements, and processing the feature data by using an interpolation algorithm to generate geological profile data corresponding to the continuous topographic features.
Preferably, generating a local geologic model corresponding to each of the topographical features based on the sub-borehole data and the geologic profile data comprises:
constructing a geological point data set corresponding to the topographic features based on the sub-borehole data and the geological profile data;
Generating a geological curved surface model dataset corresponding to each geological layer in the topographic features according to the ground particle dataset; the geological curved surface model data set is used for representing three-dimensional space distribution of each geological layer;
Extracting curved surface contour line data of each geological layer;
Taking one geological layer as a reference geological layer, and generating a side elevation data set corresponding to a topographic feature based on curved surface contour line data corresponding to the reference geological layer and the curved surface contour line data corresponding to the adjacent geological layer of the reference geological layer;
And combining the geological curved surface model data set and the side elevation data set to generate a local three-dimensional model corresponding to the topographic features.
Preferably, setting the intersecting surface of the three-dimensional model of the water area and the three-dimensional model of the plain as the fracture surface includes the following:
Performing gridding treatment on the water area three-dimensional model and the plain three-dimensional model;
Traversing all grid surfaces of the water area three-dimensional model and the plain three-dimensional model, and calculating the minimum distance between the two model grid surfaces;
When the minimum distance between the two grid surfaces is smaller than or equal to a preset threshold value, determining that the model areas represented by the two grid surfaces are contacted, and marking the two grid surfaces and the associated vertexes and edges as components of potential fracture surfaces;
And carrying out topology analysis on the marked grid surface, finding out a continuous contact surface network, and confirming that the intersecting surface of the water area three-dimensional model and the plain three-dimensional model is a fracture surface.
Preferably, the matching of the first contour line and the second contour line by the curve matching method, the obtaining of the corresponding relation between the first contour line and the second contour line includes the following contents:
Preprocessing the extracted first contour line and second contour line, wherein the preprocessing comprises removing noise, smoothing curve, normalizing coordinates and calculating characteristic points;
Calculating a feature description factor for each feature point; wherein the feature descriptors include curvature, arc length, normal vector, local shape index;
And based on the feature description factors, obtaining the corresponding relation between the first contour line and the second contour line by adopting a curve matching algorithm, wherein the corresponding relation comprises mapping one point on the first contour line to a corresponding point on the second contour line.
Preferably, the rotation and translation parameters of the plain three-dimensional model obtained based on the correspondence between the first contour line and the second contour line include the following:
Selecting a plurality of corresponding point pairs from the corresponding relation between the first contour line and the second contour line as a basis for calculating rotation and translation parameters; each corresponding point pair comprises two points which come from different local three-dimensional models and have corresponding relations;
respectively extracting three-dimensional coordinates of corresponding point pairs on the water area three-dimensional model and the plain three-dimensional model;
calculating the center point of the corresponding point pair set based on the three-dimensional coordinates, and taking the center point as a reference point for rotation and translation operation;
translating coordinates of all points in the corresponding point pair set into a new coordinate system taking the center point as an origin based on the center point;
calculating covariance matrixes of the corresponding point pairs under a new coordinate system;
Performing eigenvalue decomposition on the covariance matrix to obtain three eigenvalues and corresponding eigenvectors thereof; the characteristic vector represents the directions of x, y and z axes in the original coordinate system in the new coordinate system, and the characteristic value represents the dispersion degree of the corresponding axes in the new coordinate system;
According to the feature vector, calculating a rotation angle of the plain three-dimensional model relative to the water area three-dimensional model, wherein cosine and sine values of the rotation angle form elements of a rotation matrix;
And constructing a rotation matrix of the plain three-dimensional model relative to the water area three-dimensional model, and calculating displacement vectors between the rotated plain three-dimensional model reference points and corresponding points of the water area three-dimensional model.
Preferably, the splicing of the water area three-dimensional model and the plain three-dimensional model is completed through rotation and translation parameters, and the three-dimensional visualized model of the area to be modeled is generated and then further comprises the following contents:
and performing smoothing operation on the model splicing boundary and the adjacent areas thereof, and filling the data holes by adopting an interpolation algorithm.
The second aspect of the present application provides a three-dimensional model visualization system for a geological archive, which is applied to the three-dimensional model visualization method for a geological archive, and comprises the following steps:
the data acquisition module is used for acquiring drilling data of an area to be molded, wherein the area to be molded comprises a plurality of topographic features, the drilling data comprises a plurality of sub-drilling data sets, and each sub-drilling data set corresponds to each topographic feature;
The geological profile generation module is used for generating geological profile data corresponding to the topographic features according to the sub-drilling data set;
A local geologic model construction module for generating a local geologic model corresponding to each of the topographical features based on the sub-borehole data and the geologic profile data, wherein the plurality of topographical features includes a water area and a plains, the local geologic model includes a water area three-dimensional model and a plains three-dimensional model;
the model splicing module is used for setting a water area three-dimensional model as a reference model, setting a surface of the water area three-dimensional model intersected with a plain three-dimensional model as a fracture surface, and extracting a first boundary contour line of the water area three-dimensional model and a second contour line of the plain three-dimensional model intersected with the fracture surface; matching the first contour line and the second contour line by a curve matching method to obtain the corresponding relation between the first contour line and the second contour line; obtaining rotation and translation parameters of the plain three-dimensional model based on the corresponding relation between the first contour line and the second contour line; and (3) splicing the water area three-dimensional model and the plain three-dimensional model through rotation and translation parameters, and generating a three-dimensional visual model of the area to be modeled.
The invention has the beneficial effects that at least one of the following is adopted: and (3) by acquiring drilling data of the region to be modeled and combining geological profile data, generating a local geological model accurately corresponding to the topographic features, such as a water area three-dimensional model and a plain three-dimensional model. The spatial distribution, morphological characteristics and internal structure of the geologic body are reflected more truly and finely.
The method comprises the steps of setting a reference model to identify a fracture surface, extracting a contour line, determining a corresponding relation of the contour line by using a curve matching method, and further calculating rotation and translation parameters of a plain three-dimensional model, so that accurate model splicing is realized. The process can effectively eliminate the discontinuity on the geometric shape, the topological structure or the spatial position of the model boundary, and improves the integrity and the consistency of the whole model.
Drawings
Fig. 1 is a flowchart of a method according to a first embodiment of the present invention.
Detailed Description
Embodiments of the present invention will be described in detail below with reference to the accompanying drawings.
Example 1
As shown in fig. 1, a three-dimensional visualization method for a model of a geological archive includes the following steps:
s101, acquiring drilling data of a region to be molded.
The area to be molded comprises a plurality of topographic features, the drilling data comprises a plurality of sub-drilling data sets, and each sub-drilling data set corresponds to each topographic feature.
The area to be modeled may be a geological area including lakes, rivers and surrounding plains. The region has been subjected to multiple borehole samples, resulting in a rich borehole data set. These data sets are stored in categories according to topographical features (waters and plains), such as "lake borehole data sets", "river borehole data sets" and "plains borehole data sets". Each dataset contains position coordinates (latitude and longitude) of the respective borehole, borehole depth, and formation information sampled at different depths (e.g., rock type, thickness, mine content, etc.). The data are collected through geological investigation equipment and are recorded into a geological archive database through professional arrangement, so that basic data support is provided for subsequent modeling work.
S102, generating geological profile data corresponding to the topographic features according to the sub-drilling data set.
In one possible embodiment, generating geological profile data corresponding to the topographical features from the sub-borehole dataset includes: feature data associated with the particular topographical feature is screened from the sub-borehole data set, wherein the feature data includes borehole coordinates, depth, formation information.
And determining the position and the direction of the geological profile to be generated according to the topographic features and modeling requirements, and processing the feature data by using an interpolation algorithm to generate geological profile data corresponding to the continuous topographic features.
Referring to the lake borehole data set, feature data related to the lake topography features, including coordinates, depths and formation information of each borehole, are first screened from the data set. According to the shape, the size and modeling requirements of the lake (for example, the distribution of the lake bottom sediment needs to be displayed), a drill hole is determined to be selected at certain intervals in the long axis direction of the lake as a section starting point, and a plurality of geological section lines are generated. Then, a Kriging interpolation algorithm commonly used in the prior art is applied to the drilling data on the selected profile line to convert discrete drilling data points into continuous geological profile data, and the stratum structure and the material distribution condition of the lake bottom are clearly displayed. Similarly, the river and plain borehole data sets are processed accordingly to generate corresponding geologic profile data.
And S103, generating a local geological model corresponding to each terrain feature based on the sub-drilling data and the geological profile data.
Wherein the plurality of topographical features comprises a body of water and a plains, and the local geologic model comprises a body of water three-dimensional model and a plains three-dimensional model.
In one possible implementation, generating a local geologic model corresponding to each of the topographical features based on the sub-borehole data and the geological profile data includes: and constructing a geological point data set corresponding to the topographic features based on the sub-borehole data and the geological profile data. Generating a geological curved surface model dataset corresponding to each geological layer in the topographic features according to the ground particle dataset; the geological curved surface model data set is used for representing three-dimensional space distribution of each geological layer. And extracting curved surface contour line data of each geological layer, taking one geological layer as a reference geological layer, generating a side elevation data set corresponding to the topographic features based on the curved surface contour line data corresponding to the reference geological layer and the curved surface contour line data corresponding to the adjacent geological layer of the reference geological layer, and generating a local three-dimensional model corresponding to the topographic features by combining the geological curved surface model data set and the side elevation data set.
Referring to the above, for a water area (lake, river) section, a ground point data set including the formation interface point coordinates corresponding to each drilling position is first constructed based on its drilling data and the generated geological profile data. These points are then used to generate three-dimensional curved surface models of various geologic formations, such as sedimentary formations of the lake bottom, base strata, etc. And (3) extracting contour line data of each curved surface, selecting one layer (such as a lake bottom sediment layer) as a reference layer, and generating a lake side elevation data set by combining the contour line data of the reference layer and adjacent layers thereof. And finally, combining the geological curved surface model data set and the side elevation data set to generate a complete lake three-dimensional model.
Similarly, the plain drilling data set is similarly processed, a geological point data set comprising a soil layer, a sand layer, a base rock layer and the like is constructed, three-dimensional curved surface models of all layers are generated, contour lines are extracted, a side elevation data set is built, and finally the plain three-dimensional model is synthesized.
S104, setting the water area three-dimensional model as a reference model, setting the intersecting surface of the water area three-dimensional model and the plain three-dimensional model as a fracture surface, and extracting a first boundary contour line of the water area three-dimensional model and a second contour line of the outer surface of the plain three-dimensional model intersecting the fracture surface.
In one possible embodiment, the setting of the intersection surface of the three-dimensional model of the water area and the three-dimensional model of the plain as the fracture surface includes the following: and performing gridding treatment on the water area three-dimensional model and the plain three-dimensional model. Traversing all grid surfaces of the water area three-dimensional model and the plain three-dimensional model, and calculating the minimum distance between the two model grid surfaces. And when the minimum distance between the two grid surfaces is smaller than or equal to a preset threshold value, determining that the model areas represented by the two grid surfaces are contacted, and marking the two grid surfaces and the associated vertexes and edges as components of the potential fracture surfaces. And carrying out topology analysis on the marked grid surface, finding out a continuous contact surface network, and confirming that the intersecting surface of the water area three-dimensional model and the plain three-dimensional model is a fracture surface.
By reference, the three-dimensional model of the lake which is already constructed is set as a reference model. Then, the lake and the plain three-dimensional model are subjected to gridding treatment, and the lake and the plain three-dimensional model are divided into a plurality of small grid surfaces. And traversing all the grid surfaces, and calculating the minimum distance between the lake model grid surface and the plain model grid surface. When the distance between the two mesh surfaces is less than or equal to a set threshold (e.g., a reasonable margin of error in the thickness of the geologic formation), then the two models are considered to be in contact there, marking the mesh surfaces and their associated vertices and edges as components of a potentially fractured surface. By performing topology analysis on the marked grid surfaces, a continuous contact surface network is identified, and the surface where the lake model and the plain model actually intersect is confirmed to be a fracture surface. At the same time, a first boundary contour line of the lake model (namely, a lake shoreline connected with the plain) and a second contour line of the outer surface of the plain model intersecting with the fracture surface (namely, projection of the lake shoreline on the plain model) are extracted.
It should be noted that, the "first boundary contour line" herein refers to an outer boundary of the lake model, that is, an edge line where the lake contacts with the surrounding land (plain), and generally represents a closed curve, which may be referred to as a "lake shoreline". The line is the projection boundary of the lake model in the direction vertical to the ground, and the actual shape and boundary position of the lake are intuitively reflected. This line is extracted to define the contact boundary when the lake model is docked with the plain model.
A second contour line (i.e. the projection of the shoreline on the plain model) of the intersection of the outer surface of the plain model and the fracture surface
The "second contour line" refers to a boundary line corresponding to the intersection of the plain model and the lake model. Since the lake model and the plain model are independently generated based on the respective geological data at the time of constructing the model, they are not precisely aligned in the initial state. In order to achieve seamless splicing between the two, a part corresponding to the lake shoreline of the lake model on the plain model needs to be found.
The "fracture plane" is determined in the previous step and identifies the interface where the lake model and the plain model actually contact or cross each other in three-dimensional space. On a plain model, this interface appears as a curve formed by the intersection of the model outer surface with the fracture surface. Since this curve corresponds to the topography of the lake edge on the plains model, it can be understood as the projection of the lake shoreline on the plains model, i.e. its geometrical mapping on the surface of the plains model as seen from the lake shoreline perspective of the lake model.
S105, matching the first contour line and the second contour line through a curve matching method to obtain the corresponding relation between the first contour line and the second contour line.
In a possible implementation manner, matching the first contour line and the second contour line through a curve matching method to obtain the corresponding relation between the first contour line and the second contour line comprises the following steps of preprocessing the extracted first contour line and the extracted second contour line, wherein the preprocessing comprises removing noise, smoothing curves, normalizing coordinates and calculating feature points, and calculating feature description factors for each feature point; the feature description factors comprise curvature, arc length, normal vector and local shape index, and a curve matching algorithm is adopted to obtain the corresponding relation between the first contour line and the second contour line based on the feature description factors, wherein the corresponding relation comprises mapping of one point on the first contour line to a corresponding point on the second contour line.
By way of reference, the first contour line (lakeshoreline) and the second contour line (lakeshoreline projection) extracted are preprocessed, including noise data filtering, smoothing the curve, unifying the coordinate system, and calculating key feature points (e.g., turning points, peak points). Feature descriptors, such as curvature, arc length, normal vector, and local shape index, are then calculated for each feature point. And (3) using curve matching algorithms such as a Gaussian-Markov optical flow method and the like commonly used in the prior art, and finding out the corresponding relation between the lake shoreline and the lake shoreline projection according to the characteristic description factors, namely determining the corresponding point of each point on the lake shoreline projection.
And S106, obtaining rotation and translation parameters of the plain three-dimensional model based on the corresponding relation between the first contour line and the second contour line.
In one possible implementation, the rotation and translation parameters of the plain three-dimensional model obtained based on the correspondence between the first contour line and the second contour line include the following: selecting a plurality of corresponding point pairs from the corresponding relation between the first contour line and the second contour line as a basis for calculating rotation and translation parameters; wherein each corresponding point pair comprises two points which come from different local three-dimensional models and have corresponding relations. Respectively extracting three-dimensional coordinates of corresponding point pairs on a water area three-dimensional model and an original three-dimensional model, calculating central points of a corresponding point pair set based on the three-dimensional coordinates, taking the central points as reference points for rotation and translation operation, translating coordinates of all points in the corresponding point pair set into a new coordinate system taking the central points as an original point based on the central points, calculating covariance matrixes of the corresponding point pairs under the new coordinate system, and carrying out eigenvalue decomposition on the covariance matrixes to obtain three eigenvalues and eigenvectors corresponding to the three eigenvalues; the characteristic vector represents the directions of x, y and z axes in an original coordinate system in a new coordinate system, the characteristic value represents the dispersion degree of corresponding axes in the new coordinate system, the rotation angle of the plain three-dimensional model relative to the water area three-dimensional model is calculated according to the characteristic vector, cosine and sine values of the rotation angle form elements of a rotation matrix, the rotation matrix of the plain three-dimensional model relative to the water area three-dimensional model is constructed, and the displacement vector between a reference point of the plain three-dimensional model and a corresponding point of the water area three-dimensional model after rotation is calculated.
By reference, a plurality of corresponding point pairs are selected from the obtained corresponding relation between the first contour line and the second contour line and serve as the basis for calculating rotation and translation parameters. And extracting three-dimensional coordinates of each pair of corresponding points on the lake model and the plain model respectively. The center points of these corresponding pairs of points are calculated as reference points for rotation and translation operations. The coordinates of all corresponding point pairs are translated into a new coordinate system with the reference point as the origin. And under the new coordinate system, calculating a covariance matrix of the corresponding point pair set, and decomposing the characteristic values to obtain three characteristic values and corresponding characteristic vectors. These eigenvectors represent the directions of the x, y, z axes in the original coordinate system in the new coordinate system, and the eigenvalue size represents the dispersion degree of the corresponding axes in the new coordinate system. Based on these feature vectors, the rotation angle (Euler angle) of the plain three-dimensional model relative to the lake model is calculated, and a rotation matrix is constructed based on this. And meanwhile, calculating a displacement vector between the rotated plain model reference point and the lake model corresponding point to serve as a translation parameter.
And S107, splicing the water area three-dimensional model and the plain three-dimensional model through rotation and translation parameters, and generating a three-dimensional visual model of the area to be modeled.
In one possible implementation manner, the splicing of the water area three-dimensional model and the plain three-dimensional model is completed through rotation and translation parameters, and the three-dimensional visualization model of the area to be modeled is generated by the following steps: and performing smoothing operation on the model splicing boundary and the adjacent areas thereof, and filling the data holes by adopting an interpolation algorithm.
By reference, the plain model rotation matrix and translation vector obtained in the step S106 are utilized to perform corresponding rotation and translation operations on the plain three-dimensional model, so that the plain three-dimensional model is accurately aligned with the lake model. And (3) applying smoothing algorithms such as high-order spline interpolation and the like to the spliced boundary and the adjacent region of the model, eliminating possible boundary discontinuity and ensuring smooth transition of the surface of the model. And filling data into a cavity area possibly with data loss in the model by adopting methods such as Kriging interpolation and the like to ensure that the whole model has no obvious blank, thereby forming a complete and coherent three-dimensional visualization model of the area to be modeled.
Embodiment two:
A model three-dimensional visualization system of a geological archive is applied to the model three-dimensional visualization method of a geological archive, which comprises the following steps:
the data acquisition module is used for acquiring drilling data of an area to be molded, wherein the area to be molded comprises a plurality of topographic features, the drilling data comprises a plurality of sub-drilling data sets, and each sub-drilling data set corresponds to each topographic feature;
The geological profile generation module is used for generating geological profile data corresponding to the topographic features according to the sub-drilling data set;
A local geologic model construction module for generating a local geologic model corresponding to each of the topographical features based on the sub-borehole data and the geologic profile data, wherein the plurality of topographical features includes a water area and a plains, the local geologic model includes a water area three-dimensional model and a plains three-dimensional model;
the model splicing module is used for setting a water area three-dimensional model as a reference model, setting a surface of the water area three-dimensional model intersected with a plain three-dimensional model as a fracture surface, and extracting a first boundary contour line of the water area three-dimensional model and a second contour line of the plain three-dimensional model intersected with the fracture surface; matching the first contour line and the second contour line by a curve matching method to obtain the corresponding relation between the first contour line and the second contour line; obtaining rotation and translation parameters of the plain three-dimensional model based on the corresponding relation between the first contour line and the second contour line; and (3) splicing the water area three-dimensional model and the plain three-dimensional model through rotation and translation parameters, and generating a three-dimensional visual model of the area to be modeled.
The foregoing examples merely illustrate specific embodiments of the invention, which are described in greater detail and are not to be construed as limiting the scope of the invention. It should be noted that it will be apparent to those skilled in the art that several variations and modifications can be made without departing from the spirit of the invention, which are all within the scope of the invention.
Claims (8)
1. The three-dimensional visualization method for the model of the geological archive is characterized by comprising the following steps of:
Obtaining drilling data of an area to be molded, wherein the area to be molded comprises a plurality of topographic features, the drilling data comprises a plurality of sub-drilling data sets, and each sub-drilling data set corresponds to each topographic feature;
generating geological profile data corresponding to the topographic features according to the sub-drilling data set;
generating a local geological model corresponding to each of the topographical features based on the sub-borehole data and the geological profile data, wherein the topographical features comprise a water area and a plains, and the local geological model comprises a water area three-dimensional model and a plains three-dimensional model;
Setting a water area three-dimensional model as a reference model, setting a surface of the water area three-dimensional model intersecting with a plain three-dimensional model as a fracture surface, and extracting a first boundary contour line of the water area three-dimensional model and a second contour line of the plain three-dimensional model intersecting with the fracture surface;
Matching the first contour line and the second contour line by a curve matching method to obtain the corresponding relation between the first contour line and the second contour line;
Obtaining rotation and translation parameters of the plain three-dimensional model based on the corresponding relation between the first contour line and the second contour line;
and (3) splicing the water area three-dimensional model and the plain three-dimensional model through rotation and translation parameters, and generating a three-dimensional visual model of the area to be modeled.
2. A method for three-dimensional visualization of a model of a geological archive as claimed in claim 1, wherein: generating geological profile data corresponding to the topographical features from the sub-borehole dataset includes:
Screening feature data related to specific topographic features from the sub-borehole data set, wherein the feature data comprises coordinates, depth and formation information of the borehole;
and determining the position and the direction of the geological profile to be generated according to the topographic features and modeling requirements, and processing the feature data by using an interpolation algorithm to generate geological profile data corresponding to the continuous topographic features.
3. A method for three-dimensional visualization of a model of a geological archive as claimed in claim 1, wherein: based on the sub-borehole data and the geologic profile data, generating a local geologic model corresponding to each of the topographical features includes:
constructing a geological point data set corresponding to the topographic features based on the sub-borehole data and the geological profile data;
Generating a geological curved surface model dataset corresponding to each geological layer in the topographic features according to the ground particle dataset; the geological curved surface model data set is used for representing three-dimensional space distribution of each geological layer;
Extracting curved surface contour line data of each geological layer;
Taking one geological layer as a reference geological layer, and generating a side elevation data set corresponding to a topographic feature based on curved surface contour line data corresponding to the reference geological layer and the curved surface contour line data corresponding to the adjacent geological layer of the reference geological layer;
And combining the geological curved surface model data set and the side elevation data set to generate a local three-dimensional model corresponding to the topographic features.
4. A method for three-dimensional visualization of a model of a geological archive as claimed in claim 1, wherein: setting the intersecting surface of the water area three-dimensional model and the plain three-dimensional model as a fracture surface comprises the following contents:
Performing gridding treatment on the water area three-dimensional model and the plain three-dimensional model;
Traversing all grid surfaces of the water area three-dimensional model and the plain three-dimensional model, and calculating the minimum distance between the two model grid surfaces;
When the minimum distance between the two grid surfaces is smaller than or equal to a preset threshold value, determining that the model areas represented by the two grid surfaces are contacted, and marking the two grid surfaces and the associated vertexes and edges as components of potential fracture surfaces;
And carrying out topology analysis on the marked grid surface, finding out a continuous contact surface network, and confirming that the intersecting surface of the water area three-dimensional model and the plain three-dimensional model is a fracture surface.
5. A method of three-dimensional visualization of a model of a geological archive as defined in claim 4, wherein: matching the first contour line and the second contour line by a curve matching method, and obtaining the corresponding relation between the first contour line and the second contour line comprises the following contents:
Preprocessing the extracted first contour line and second contour line, wherein the preprocessing comprises removing noise, smoothing curve, normalizing coordinates and calculating characteristic points;
Calculating a feature description factor for each feature point; wherein the feature descriptors include curvature, arc length, normal vector, local shape index;
And based on the feature description factors, obtaining the corresponding relation between the first contour line and the second contour line by adopting a curve matching algorithm, wherein the corresponding relation comprises mapping one point on the first contour line to a corresponding point on the second contour line.
6. A method for three-dimensional visualization of a model of a geological archive as defined in claim 5, wherein: the rotation and translation parameters of the plain three-dimensional model obtained based on the corresponding relation between the first contour line and the second contour line comprise the following contents:
Selecting a plurality of corresponding point pairs from the corresponding relation between the first contour line and the second contour line as a basis for calculating rotation and translation parameters; each corresponding point pair comprises two points which come from different local three-dimensional models and have corresponding relations;
respectively extracting three-dimensional coordinates of corresponding point pairs on the water area three-dimensional model and the plain three-dimensional model;
calculating the center point of the corresponding point pair set based on the three-dimensional coordinates, and taking the center point as a reference point for rotation and translation operation;
translating coordinates of all points in the corresponding point pair set into a new coordinate system taking the center point as an origin based on the center point;
calculating covariance matrixes of the corresponding point pairs under a new coordinate system;
Performing eigenvalue decomposition on the covariance matrix to obtain three eigenvalues and corresponding eigenvectors thereof; the characteristic vector represents the directions of x, y and z axes in the original coordinate system in the new coordinate system, and the characteristic value represents the dispersion degree of the corresponding axes in the new coordinate system;
According to the feature vector, calculating a rotation angle of the plain three-dimensional model relative to the water area three-dimensional model, wherein cosine and sine values of the rotation angle form elements of a rotation matrix;
And constructing a rotation matrix of the plain three-dimensional model relative to the water area three-dimensional model, and calculating displacement vectors between the rotated plain three-dimensional model reference points and corresponding points of the water area three-dimensional model.
7. A method of three-dimensional visualization of a model of a geological archive as defined in claim 6, wherein: the splicing of the water area three-dimensional model and the plain three-dimensional model is completed through rotation and translation parameters, and the three-dimensional visualized model of the area to be modeled is generated and then comprises the following contents:
and performing smoothing operation on the model splicing boundary and the adjacent areas thereof, and filling the data holes by adopting an interpolation algorithm.
8. A model three-dimensional visualization system of a geological archive, applied to a model three-dimensional visualization method of a geological archive according to any one of claims 1 to 7, comprising:
the data acquisition module is used for acquiring drilling data of an area to be molded, wherein the area to be molded comprises a plurality of topographic features, the drilling data comprises a plurality of sub-drilling data sets, and each sub-drilling data set corresponds to each topographic feature;
The geological profile generation module is used for generating geological profile data corresponding to the topographic features according to the sub-drilling data set;
A local geologic model construction module for generating a local geologic model corresponding to each of the topographical features based on the sub-borehole data and the geologic profile data, wherein the plurality of topographical features includes a water area and a plains, the local geologic model includes a water area three-dimensional model and a plains three-dimensional model;
the model splicing module is used for setting a water area three-dimensional model as a reference model, setting a surface of the water area three-dimensional model intersected with a plain three-dimensional model as a fracture surface, and extracting a first boundary contour line of the water area three-dimensional model and a second contour line of the plain three-dimensional model intersected with the fracture surface; matching the first contour line and the second contour line by a curve matching method to obtain the corresponding relation between the first contour line and the second contour line; obtaining rotation and translation parameters of the plain three-dimensional model based on the corresponding relation between the first contour line and the second contour line; and (3) splicing the water area three-dimensional model and the plain three-dimensional model through rotation and translation parameters, and generating a three-dimensional visual model of the area to be modeled.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410955946.7A CN118505921B (en) | 2024-07-17 | 2024-07-17 | Model three-dimensional visualization method and system for geological archives |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410955946.7A CN118505921B (en) | 2024-07-17 | 2024-07-17 | Model three-dimensional visualization method and system for geological archives |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN118505921A true CN118505921A (en) | 2024-08-16 |
| CN118505921B CN118505921B (en) | 2024-10-11 |
Family
ID=92238617
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410955946.7A Active CN118505921B (en) | 2024-07-17 | 2024-07-17 | Model three-dimensional visualization method and system for geological archives |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118505921B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111968231A (en) * | 2020-08-14 | 2020-11-20 | 成都理工大学 | Three-dimensional stratum modeling method based on geological map section |
| CN112381937A (en) * | 2020-11-03 | 2021-02-19 | 广东省地质调查院 | Multi-source geological data coupling modeling method based on drilling and complex geological profile |
| CN116630568A (en) * | 2023-07-24 | 2023-08-22 | 航天宏图信息技术股份有限公司 | Geological three-dimensional model generation method and device based on point, line and surface structures |
| WO2023185734A1 (en) * | 2022-03-28 | 2023-10-05 | 中国矿业大学 | Mine digital model construction method based on three-dimensional point cloud database |
| CN117456122A (en) * | 2023-12-22 | 2024-01-26 | 江西省地质博物馆 | Three-dimensional geological model generation method and system |
-
2024
- 2024-07-17 CN CN202410955946.7A patent/CN118505921B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111968231A (en) * | 2020-08-14 | 2020-11-20 | 成都理工大学 | Three-dimensional stratum modeling method based on geological map section |
| CN112381937A (en) * | 2020-11-03 | 2021-02-19 | 广东省地质调查院 | Multi-source geological data coupling modeling method based on drilling and complex geological profile |
| WO2023185734A1 (en) * | 2022-03-28 | 2023-10-05 | 中国矿业大学 | Mine digital model construction method based on three-dimensional point cloud database |
| CN116630568A (en) * | 2023-07-24 | 2023-08-22 | 航天宏图信息技术股份有限公司 | Geological three-dimensional model generation method and device based on point, line and surface structures |
| CN117456122A (en) * | 2023-12-22 | 2024-01-26 | 江西省地质博物馆 | Three-dimensional geological model generation method and system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN118505921B (en) | 2024-10-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106934860B (en) | Three-dimensional geological modeling method based on T-spline | |
| Guo et al. | Three-dimensional geological modeling and spatial analysis from geotechnical borehole data using an implicit surface and marching tetrahedra algorithm | |
| Caumon et al. | Building and editing a sealed geological model | |
| Lemon et al. | Building solid models from boreholes and user-defined cross-sections | |
| CN106844871A (en) | Three-dimensional geotechnical engineering investigation information model construction method based on BIM | |
| CN108986213B (en) | Three-dimensional stratum modeling method based on stacking technology | |
| Gold et al. | A spatial data structure integrating GIS and simulation in a marine environment | |
| CN111612903A (en) | Geological data visualization method based on mixed data model | |
| Zheng et al. | Three-dimensional unstructured mesh generation: Part 1. Fundamental aspects of triangulation and point creation | |
| CN118298117A (en) | Three-dimensional multi-level geological model construction method | |
| CN113628339A (en) | Broken layer slice extraction method based on local reservation projection | |
| CN118505921B (en) | Model three-dimensional visualization method and system for geological archives | |
| de Oliveira Miranda et al. | Finite element mesh generation for subsurface simulation models | |
| RU2484427C1 (en) | Method for cartographic display of two-dimensional distributions given in digital form | |
| Jones | Solid modelling of earth masses for applications in geotechnical engineering | |
| Lee | Automatic metric 3D surface mesh generation using subdivision surface geometrical model. Part 1: Construction of underlying geometrical model | |
| CN113658333A (en) | Geologic body modeling method based on isosurface extraction | |
| CN114494637B (en) | Sandstone three-dimensional real model reconstruction method based on structural matrix | |
| EP3262444A1 (en) | Method for parameterizing a 3d domain with discontinuities | |
| Lee | Automatic metric 3D surface mesh generation using subdivision surface geometrical model. Part 2: Mesh generation algorithm and examples | |
| Wang et al. | Representing the geological body heterogeneous property field using the quadratic generalized tri-prism volume function model (QGTPVF) | |
| Zhong et al. | NURBS reconstruction of digital terrain for hydropower engineering based on TIN model | |
| Matuk et al. | Skeleton based contour line generalization | |
| Miao et al. | Automatic generation method of geological cross-sections in dredging engineering based on 3D solid NURBS models | |
| Li et al. | Alternative 3D modeling approaches based on complex multi-source geological data interpretation |
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 |