Geological body model construction method, device and equipment and storage medium

Abstract

The application discloses a geological body model construction method, device and equipment and a storage medium. The method comprises the following steps: performing model fitting processing on an initial geological body model constructed by using an implicit modeling technology based on geological constraint points to obtain a transition model; and performing model fitting processing on the transition model based on geological constraint lines to obtain a geological body model with feature fitting. The geological constraint points are directional interpolation points for implicit surface interpolation, and the geological constraint points have a tangent constraint direction and a normal constraint direction. The geological constraint lines comprise non-boundary constraint lines and boundary constraint lines. The non-boundary constraint lines and the boundary constraint lines have a tangent constraint direction and a normal constraint direction. The non-boundary constraint lines are non-boundary interpolation lines for implicit surface interpolation, and the boundary constraint lines are boundary interpolation lines for implicit surface interpolation.

Description

Technical Field

[0001] This application relates to the field of geological body modeling, and in particular to a method, apparatus, equipment and storage medium for constructing geological body models. Background Technology

[0002] In the field of geological body modeling, geological sampling data is extremely useful and invaluable, and its accurate and effective participation in the modeling process should be ensured. For mines, detailed 3D modeling of geological bodies is a crucial guarantee for estimating mineral resource reserves. Implicit modeling techniques can interpolate geological body surface models using modeling constraints derived from geological feature transformations, thereby effectively improving the efficiency of geological body modeling and facilitating dynamic model updates. However, in practical applications, the geological body mesh model obtained through implicit modeling techniques is not easily and accurately fitted to the geological feature points and feature lines obtained from geological sampling data. This is mainly due to two factors: the implicit surface interpolation process and the implicit surface reconstruction process.

[0003] Implicit surface reconstruction primarily employs isosurface extraction. Generally, there are two main methods for ensuring the extracted isosurfaces fit the given feature points and feature lines in the model fitting process, employing pre-processing and post-processing approaches respectively. Pre-processing directly fits the feature points and feature lines during isosurface extraction. Post-processing, on the other hand, indirectly fits the resulting mesh model using feature points and feature lines after surface reconstruction. Pre-processing requires modification of the isosurface extraction process for specific methods, while post-processing can fit mesh models obtained by any isosurface extraction method, offering wider applicability. However, the relevant isosurface extraction methods often result in the geological body mesh model not accurately fitting to the geological feature points and feature lines obtained from geological sampling data, meaning the constructed geological body model lacks precision. Summary of the Invention

[0004] In view of this, embodiments of this application provide a method, apparatus, device and storage medium for constructing geological body models, aiming to improve the construction accuracy of geological body models.

[0005] The technical solution of this application embodiment is implemented as follows:

[0006] In a first aspect, embodiments of this application provide a method for constructing a geological body model, including:

[0007] The initial geological body model constructed using implicit modeling techniques is subjected to model fitting processing based on geological constraint points to obtain a transitional model;

[0008] The transition model is subjected to model fitting processing based on geological constraint lines to obtain a geological body model with feature fitting;

[0009] Wherein, the geological constraint point is a directional interpolation point used for implicit surface interpolation, and the geological constraint point has a tangential constraint direction and a normal constraint direction; the geological constraint line includes: a non-boundary constraint line and a boundary constraint line, both of which have a tangential constraint direction and a normal constraint direction, the non-boundary constraint line is a non-boundary interpolation line used for implicit surface interpolation, and the boundary constraint line is a boundary interpolation line used for implicit surface interpolation.

[0010] In some embodiments, the process of performing model fitting processing on the initial geological body model constructed using implicit modeling techniques based on geological constraint points to obtain a transitional model includes:

[0011] Obtain the constraint points of the initial geological body model within the set fitting range;

[0012] Calculate the distance from each constraint point to the mesh surface of the initial geological body model, and retain the constraint points whose distance is less than or equal to the maximum fitting distance as the geological constraint points;

[0013] The initial geological body model is subjected to model fitting processing by moving projection points to corresponding geological constraint points to obtain a transition model; wherein, the projection point is the point on the initial geological body model projected onto the corresponding geological constraint point, and the projection point was a mesh vertex of the initial geological body model before the movement.

[0014] In some embodiments, the process of performing model fitting based on geological constraint lines on the transition model to obtain a feature-fitted geological body model includes:

[0015] The transition model is subjected to model fitting processing based on the non-boundary constraint line and the boundary constraint line respectively to obtain a geological body model with feature fitting.

[0016] In some embodiments, the transition model is subjected to model fitting processing based on the boundary constraint lines, including:

[0017] Based on the boundary projection line projected onto the transition model surface by the boundary constraint line, a first projection neighborhood corresponding to the boundary constraint line is generated. The first projection neighborhood is the region where the set of topologically continuous triangular facets on the path where the boundary projection line is located is located.

[0018] Re-mesh the first projected neighborhood;

[0019] The remeshed triangular facets within the first projection neighborhood are merged with adjacent triangular facets, and the normals of the triangular facets are made uniform.

[0020] In some embodiments, the transition model is subjected to model fitting processing based on the non-boundary constraint lines, including:

[0021] Based on the non-boundary projection line projected onto the surface of the transition model, a second projection neighborhood corresponding to the non-boundary constraint line is generated. The second projection neighborhood is the region where the set of topologically continuous triangular facets on the path where the non-boundary projection line is located is located.

[0022] Remesh the second projection neighborhood;

[0023] Merge the remeshed triangular facets within the second projection neighborhood with adjacent triangular facets, and unify the normals of the triangular facets.

[0024] In some embodiments, the method further includes: trimming the model boundary of the feature-fitting geological body model to obtain a geological body model with matching boundaries.

[0025] In some embodiments, trimming the model boundary of the geological body model with the feature fit includes:

[0026] Along the direction of the boundary line of the geological body model that fits the features, start from the initial edge and traverse the edges of each boundary triangle. Based on the plane equation of the normal plane where each boundary triangle edge is located, determine whether the corresponding triangular facet is located within the model boundary.

[0027] Delete the triangular facets located outside the boundary of the model;

[0028] The initial edge is the edge formed by the selected starting point and the next adjacent point on the boundary line.

[0029] In some embodiments, determining whether a corresponding triangular facet lies within the model boundary based on the plane equation of the normal plane containing the sides of each boundary triangle includes:

[0030] If the coordinates of the three vertices of the triangular facet are substituted into the plane equation and the resulting function value is greater than or equal to zero, then the triangular facet is determined to be inside the model boundary; otherwise, the triangular facet is determined to be outside the model boundary.

[0031] Secondly, embodiments of this application provide a geological body model construction apparatus, comprising:

[0032] The first bonding module is used to bond the initial geological body model constructed using implicit modeling technology based on geological constraint points to obtain a transition model.

[0033] The second bonding module is used to perform model bonding processing on the transition model based on geological constraint lines to obtain a geological body model with feature bonding.

[0034] Wherein, the geological constraint point is a directional interpolation point used for implicit surface interpolation, and the geological constraint point has a tangential constraint direction and a normal constraint direction; the geological constraint line includes: a non-boundary constraint line and a boundary constraint line, both of which have a tangential constraint direction and a normal constraint direction, the non-boundary constraint line is a non-boundary interpolation line used for implicit surface interpolation, and the boundary constraint line is a boundary interpolation line used for implicit surface interpolation.

[0035] Thirdly, embodiments of this application provide a geological body model construction device, including: a processor and a memory for storing a computer program capable of running on the processor, wherein the processor, when running the computer program, executes the steps of the method described in the first aspect of embodiments of this application.

[0036] Fourthly, embodiments of this application provide a storage medium storing a computer program, which, when executed by a processor, implements the steps of the method described in the first aspect of embodiments of this application.

[0037] The technical solution provided in this application involves performing model fitting processing on an initial geological body model constructed using implicit modeling technology based on geological constraint points to obtain a transition model; and performing model fitting processing on the transition model based on geological constraint lines to obtain a feature-fitted geological body model. The geological constraint points are directional interpolation points used for implicit surface interpolation, and each geological constraint point has a tangential constraint direction and a normal constraint direction. The geological constraint lines include non-boundary constraint lines and boundary constraint lines, both of which have tangential and normal constraint directions. The non-boundary constraint lines are non-boundary interpolation lines used for implicit surface interpolation, and the boundary constraint lines are boundary interpolation lines used for implicit surface interpolation. Compared to directly using geological feature points and feature lines from geological interpretation for model fitting, the embodiments of this application perform mesh processing based on geological constraint points and lines derived from implicit modeling. This allows for precise fitting of model features, eliminating gaps between the model and features, improving the accuracy and reliability of geological body models, and enabling feature fitting for geological body models obtained by any isosurface extraction method, thus exhibiting broad applicability. For mines, this will help improve the reliability of mineral resource reserve estimation and mining design, solidifying the foundation for the digital and intelligent development of mines. Attached Figure Description

[0038] Figure 1 This is a flowchart illustrating the geological body model construction method of an embodiment of this application;

[0039] Figures 2A to 2C This is a schematic diagram of the implicit modeling constraints of geological bodies in an embodiment of this application, wherein, Figure 2A This is a schematic diagram of geological constraint points. Figure 2B This is a schematic diagram of non-boundary constraint lines. Figure 2C This is a schematic diagram of the boundary constraint lines;

[0040] Figure 3 This is a schematic diagram illustrating the principle of the constraint point in an application example of this application;

[0041] Figure 4 This is a schematic diagram illustrating the principle of remeshing based on the matching method of triangulation rules between two polylines in an application example of this application.

[0042] Figure 5 This is a schematic diagram illustrating the principle of conforming to boundary constraint lines in an application example of this application;

[0043] Figure 6 This is a schematic diagram illustrating the principle of fitting non-boundary constraint lines in an application example of this application;

[0044] Figure 7 This is a schematic diagram illustrating the principle of model boundary trimming for a geological body model with feature fitting in an application example of this application;

[0045] Figure 8 This is a schematic diagram illustrating the principle of conforming to geological constraint points in another application example of this application;

[0046] Figure 9 This is a schematic diagram illustrating the principle of conforming to geological constraint lines and trimming model boundaries in another application example of this application;

[0047] Figure 10 This is a schematic diagram of the geological body model construction device according to an embodiment of this application;

[0048] Figure 11 This is a schematic diagram of the geological body model construction device according to an embodiment of this application. Detailed Implementation

[0049] The present application will now be described in further detail with reference to the accompanying drawings and embodiments.

[0050] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs. The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the application.

[0051] In related technologies, implicit surface reconstruction mainly employs isosurface extraction methods. For example, taking the Marching Cubes (MC) algorithm, the MC algorithm only includes the intersection coordinates when calculating the intersection information between the isosurface and the voxel element, and constructs the geometric model within the voxel element based on the triangular facets connected by these intersections, without considering the normals of the intersections. This leads to a tendency for model feature loss during isosurface extraction. Due to the defects of such isosurface extraction algorithms, when the voxel element contains geometric model feature information (such as feature points, feature lines, etc.), the final constructed geometric model will lack this feature information, affecting the accuracy of the geological body model construction.

[0052] Based on this, in this embodiment of the application, based on the idea of ​​mesh model post-processing, taking the constraint data of implicit surfaces and the geological body mesh model generated by implicit modeling as the research object, a method for accurately fitting the features of the geological body model is proposed to solve the problem that the geological body mesh model obtained by implicit modeling technology is not easy to accurately fit the geological feature points and feature lines obtained based on geological sampling data.

[0053] like Figure 1 As shown in the figure, this application provides a method for constructing a geological body model, including:

[0054] Step 101: The initial geological body model constructed using implicit modeling technology is subjected to model fitting processing based on geological constraint points to obtain a transition model.

[0055] It should be noted that implicit modeling is a modeling technique well-suited for constructing interactive constraints, and the resulting model is relatively easy to update dynamically. This method transforms the geometric space constructed based on different geological data into an implicit function field using a distance function, expressing the 3D model as a mathematical function. The reconstructed surface is represented as the zero-level set of the implicit function. Because the amount of work required to create the model is relatively small, more time can be devoted to understanding the geological conditions and studying more complex details, such as faults, stratigraphic sequences, trends, and veins. Related implicit modeling techniques will not be elaborated upon here.

[0056] Here, geological constraint points are directional interpolation points used for implicit surface interpolation, and geological constraint points have tangential constraint directions and normal constraint directions.

[0057] It should be noted that geological constraint points are implicit function interpolation data required in implicit modeling of geological bodies, including but not limited to at least one of the following: borehole contact points, constraint line sampling points, and artificially added points.

[0058] For example, such as Figure 2AAs shown, geological constraint points can be displayed using a disk, where the direction of the disk's tilt indicates the tangential constraint direction of the geological constraint point. The inner and outer sides of the disk can be distinguished by different colors to indicate the normal constraint direction of the geological constraint point; for example, the inner side of the disk can be displayed in green, and the outer side in blue.

[0059] Step 102: Perform model fitting processing on the transition model based on geological constraint lines to obtain a geological body model with feature fitting.

[0060] Here, geological constraint lines include: non-boundary constraint lines and boundary constraint lines. Both non-boundary constraint lines and boundary constraint lines have tangential constraint directions and normal constraint directions. Non-boundary constraint lines are non-boundary interpolation lines used for implicit surface interpolation, and boundary constraint lines are boundary interpolation lines used for implicit surface interpolation.

[0061] For example, non-boundary constraint lines are as follows Figure 2B As shown, tangential and normal stripes can be used for display. The tangential constraint directions of non-boundary constraint lines include two directions: pointing inwards and outwards from the model. Boundary constraint lines are shown below. Figure 2C As shown, the tangential constraint direction of this boundary constraint line only includes the part pointing inwards from the model.

[0062] It should be noted that, compared to directly using geological feature points and feature lines from geological interpretation for model fitting, the embodiments of this application use geological constraint points and lines derived from implicit modeling for mesh processing. This allows for precise fitting of model features, eliminating gaps between the model and features, improving the accuracy and reliability of the geological body model, and enabling feature fitting for geological body models obtained by any isosurface extraction method, thus exhibiting broad applicability. For mines, this will help improve the reliability of mineral resource reserve estimation and mining design, solidifying the foundation for the digital and intelligent development of mines.

[0063] In some embodiments, the initial geological body model constructed using implicit modeling techniques is subjected to model fitting processing based on geological constraint points to obtain a transitional model, including:

[0064] Obtain the constraint points of the initial geological body model within the set fitting range;

[0065] Calculate the distance from each constraint point to the mesh surface of the initial geological body model, and retain the constraint points whose distance is less than or equal to the maximum fitting distance as geological constraint points;

[0066] The initial geological body model is fitted by moving the projection points to the corresponding geological constraint points to obtain the transition model. The projection points are the points on the initial geological body model projected onto the corresponding geological constraint points, and the projection points are the mesh vertices of the initial geological body model before the movement.

[0067] It should be noted that the above-mentioned fitting range can be understood as the distance tolerance from the constraint point to the mesh surface. The maximum fitting distance should be within the mesh surface accuracy value; for example, the maximum fitting distance can be half of the mesh table accuracy value.

[0068] In one application example, first obtain the constraint points P of the initial geological body model within the defined fitting range. i For {i = 1, 2, ..., n}, determine the maximum fitting distance d based on prior experience. max Half the mesh surface accuracy value (i.e., 1 / 2ε) mes Calculate the distance from the constraint point to the mesh surface. constraint point P i It should be projected onto the surface of the initial geological model M to obtain the distance constraint point P. i The nearest projection point Q i {i=1,2,…,n}. The distance from the constraint point to the projection point is calculated. By the maximum fitting distance d max The comparison is performed, and the constraint points within the maximum fitting distance are retained, while the constraint points outside the maximum fitting distance are discarded.

[0069] Based on the mesh processing method, to fit the initial geological body model M to the retained constraint points, this can be achieved by moving the projection points to their corresponding constraint point positions. Additionally, at projection point Q... i Before being moved, the projection point should be ensured to be a vertex of the initial geological body model M. Since there are three possible locations for the projection point: mesh vertex, mesh edge, and mesh interior, for projection points that are not mesh vertices, the triangular facet containing that point needs to be reconstructed using the projection point as the mesh vertex.

[0070] like Figure 3 As shown, after the constraint points are precisely fitted, the topological relationship of the triangular mesh on the model surface undergoes a significant change.

[0071] In some embodiments, the transition model is subjected to model fitting based on geological constraint lines to obtain a geological body model with feature fitting, including:

[0072] The transition model is subjected to model fitting processing based on non-boundary constraint lines and boundary constraint lines respectively to obtain a geological body model with feature fitting.

[0073] For example, geological constraint lines are generally derived from implicitly modeled geological feature lines, based on constraint line L. i Whether {i=1,2,…,n} lies on the boundary of the initial geological body model M can be used to divide the geological constraint line into boundary constraint lines L. b,i {i=1,2,…,n} and non-boundary constraint line L n,i {i=1,2,…,n}.

[0074] The process of conforming to geological constraint lines can be understood as remeshing the mesh matched by the geological constraint lines. Essentially, it's a process of triangulating two closed polylines. The key to remeshing is avoiding errors such as self-intersection in the remeshed triangular faces. Since the remeshing process of matching the mesh is a non-trivial problem, the quality of the mesh obtained based on implicit surface reconstruction should be improved as much as possible; otherwise, the reliability of the model conformation will be affected. Therefore, in some embodiments, remeshing can be performed based on the matching method of triangulation rules between two polylines. By matching the geometric trends of the two polylines, the quality of the triangulated mesh can be optimized. These triangulation rules refer to the pairing rules that the triangular faces formed by the vertices of the first and second polylines should satisfy, such as the minimum perimeter method and the minimum surface area method.

[0075] In one application example, such as Figure 4 As shown, the transition model is remeshed based on a matching method using triangulation rules between two polylines, including the following steps:

[0076] Step (1): Determine the initial matching point.

[0077] Determine the first polyline L f The second polyline L s The initial matching point (F) i＝0 {i=0,1,…,m-1} and S j＝0 ({j=0,1,…,n-1}). The initial matching point can be determined by searching for the nearest point between two polylines, such as... Figure 4 As shown in (b).

[0078] Step (2): Preprocess polylines.

[0079] Preprocess locations where the polyline's shape changes significantly. Sample the feature lines to ensure the point density of the two polylines is basically the same, then unify the point lists (F) of the two polylines. i and S j Adjust the direction of the polyline and the position where the polyline turns significantly to ensure that the shapes and trends of the two polylines are basically consistent, such as... Figure 4 As shown in (c).

[0080] Step (3): Construct vertex topology and mesh it.

[0081] Match the corresponding vertices of the two polylines according to the triangulation rules. Starting from the initial matching point (F... i＝0 and S j＝0 Starting with the corresponding triangulation rules, we construct the vertex topology between two polylines by taking two points on one polyline and one point on another polyline, such as... Figure 4 As shown in (d) and 4(e). Then, based on the vertex matching results, a triangular mesh topology diagram between the two polylines is constructed to obtain the remeshed mesh model, as shown in... Figure 4 As shown in (f).

[0082] In some embodiments, the transition model is subjected to model fitting processing based on boundary constraint lines, including:

[0083] Based on the boundary projection lines projected onto the surface of the transition model, a first projection neighborhood corresponding to the boundary constraint lines is generated. The first projection neighborhood is the region where the set of topologically continuous triangular facets on the path of the boundary projection lines is located.

[0084] Remesh the first projected neighborhood;

[0085] Merge the remeshed triangular facets in the first projection neighborhood with the adjacent triangular facets, and unify the normals of the triangular facets.

[0086] It should be noted that the boundary constraint line L b,i It serves to control the boundary of the geological body model. When fitting the boundary constraint line, it should be ensured that the geological body model is aligned with the boundary constraint line L. b,i There are no gaps between them, i.e., the boundary constraint line L b,i It should be precisely fitted; in other words, the model's boundary contours should maintain topological consistency with the actual boundary constraint lines.

[0087] In this embodiment, the remeshing of the first projection neighborhood adopts the aforementioned matching method based on the triangulation rules between two polylines to remesh the first projection neighborhood.

[0088] In one application example, such as Figure 5 As shown, the boundary constraint lines include:

[0089] Step (1): Determine the projection neighborhood of the boundary constraint line.

[0090] Boundary constraint line L b,i Projected onto the model surface, the boundary constraint line L is obtained. b,i Similar paths with the same number of points are denoted as boundary projection lines L′. b,i {i=1,2,…,n}. Based on the boundary projection line L′ b,i The location of the upper projection point, search boundary projection line L′b,i The set S of topologically continuous triangular faces along the path i {i=1,2,…,n}, let the set of triangular facets S i The area in question is defined as the boundary constraint line L. b,i The projected neighborhood Ω i {i=1,2,…,n}.

[0091] Step (2): Remesh the projected neighborhood of the boundary constraint lines.

[0092] Extracting the projected neighborhood Ω i At the boundary projection line L′ b,i The inner contour line is denoted as L. p,i And delete the set of triangular facets S i Using boundary constraint line L b,i The interpolation points on the grid are the grid vertices, in L b,i With L p,i Between the projected neighborhood Ω i Triangulation is performed using a remeshing method between two polylines to generate a new triangular facet S′. i {i=1,2,…,n}.

[0093] Step (3): Merge all facets of the model and unify the facet normals.

[0094] Projected neighborhood Ω i The triangular facet S′ with internal weight triangulation i Merge the untriangulated facets of the model to generate a new model M′, and unify the normals of the triangular facets of model M′.

[0095] In some embodiments, the transition model is subjected to model fitting processing based on non-boundary constraint lines, including:

[0096] Based on the non-boundary projection lines projected onto the surface of the transition model, a second projection neighborhood corresponding to the non-boundary constraint lines is generated. The second projection neighborhood is the region where the set of topologically continuous triangular facets on the path of the non-boundary projection lines is located.

[0097] Remesh the second projected neighborhood;

[0098] Merge the remeshed triangular facets in the second projection neighborhood with the adjacent triangular facets, and unify the normals of the triangular facets.

[0099] Due to the non-boundary constraint line L n,i It only affects the surface morphology of the geological body model. Within the tolerance range, as long as the geological body model fits the constraint points on the constraint line, the actual modeling requirements can be met. Based on this, when fitting the non-boundary constraint line L... n,iIn this case, the method of fitting constraint points can be used, without the need to eliminate non-boundary constraint lines L. n,i The existing gaps between the model and the boundary constraint line L. n,i The constraint points should first be projected onto the model surface to obtain the corresponding projection points. Then, based on the reconstructed mesh of the area where the projection points are located, the projection points can be moved to their corresponding constraint points.

[0100] Considering the principles of refined 3D modeling, and in order to improve the geometric quality and refinement of the model, it is required to implement the constraint lines L of the model. i Precise fit, meaning the fitted model M′ should be aligned with all constraint lines L i (including non-boundary constraint lines L) n,i and boundary constraint line L b,i There are no gaps between them, and the constraint line L i The constraint points and mesh edges on the model should maintain topological consistency with the mesh of the geological body model M′. Fitting non-boundary constraint lines can be considered as the general case of fitting boundary constraint lines.

[0101] In this embodiment, the remeshing of the second projection neighborhood is achieved by using the aforementioned matching method based on the triangulation rules between two polylines to remesh the second projection neighborhood.

[0102] In one application example, such as Figure 6 As shown, the steps for fitting non-boundary constraint lines are as follows:

[0103] Step (1): Determine the projection neighborhood of the non-boundary constraint line.

[0104] non-boundary constraint line L n,i Projected onto the model surface, we obtain the non-boundary constraint line L. n,i Paths with the same number of points are denoted as non-boundary projection lines L′. n,i {i=1,2,…,n}. Based on the non-boundary projection line L′ n,i The location of the upper projection point is used to search for the non-boundary projection line L′. n,i The set T of topologically continuous triangular faces along the path i {i=1,2,…,n}, let the set of triangular facets T i The region in question is defined as the non-boundary constraint line L. n,i The projection neighborhood R i {i=1,2,…,n}.

[0105] Step (2): Remesh the projected neighborhood of the non-boundary constraint lines.

[0106] Extracting the projected neighborhood R i At the non-boundary projection line L′ n,i The outlines on both sides are denoted as L.p,i,f and L p,i,s And delete the set of triangular facets T i Using non-boundary constraint line L n,i The constraint points on the grid are the grid vertices, respectively at L n,i With L p,i,f and L p,i,s Between the projected neighborhood R i Triangulation is performed using a remeshing method between two polylines to generate a new triangular facet T′. i {i=1,2,…,n}.

[0107] Step (3) Merge all facets of the model and unify the facet normals.

[0108] Project neighborhood R i The inner triangularized face T′ i Merge the untriangulated facets of the model to generate a new model M′, and unify the normals of the triangular facets of model M′.

[0109] Considering that if the boundary lines of a geological body model are only determined after the model is generated, the portion outside the model boundary can be trimmed through post-processing. Based on this, this application proposes a method for model boundary trimming to solve the problem of mismatched boundaries in some geological body models.

[0110] In some embodiments, the method further includes: trimming the model boundary of the feature-fitting geological body model to obtain a geological body model with matching boundaries.

[0111] It should be noted that before clipping the model boundary, in order to ensure the clipping effect, on the one hand, the boundary line should be closed and the direction of the boundary line should be consistent with the direction of the boundary constraint line; on the other hand, the model should be precisely fitted to the boundary constraint line.

[0112] For example, model boundary trimming is performed on a geological body model with feature fitting, including:

[0113] Along the direction of the boundary line of the geological body model with feature fit, start from the initial edge and traverse the edges of each boundary triangle. Based on the plane equation of the normal plane where each boundary triangle edge is located, determine whether the corresponding triangular facet is located within the model boundary.

[0114] Delete the triangular facets located outside the model boundary;

[0115] The initial edge is the edge formed by the selected starting point and the next adjacent point on the boundary line.

[0116] For example, determining whether a corresponding triangular facet lies within the model boundary based on the plane equation of the normal plane containing the sides of each boundary triangle includes:

[0117] If the coordinates of the three vertices of the triangular facet are substituted into the plane equation and the function value is greater than or equal to zero, then the triangular facet is determined to be inside the model boundary; otherwise, the triangular facet is determined to be outside the model boundary.

[0118] In one application example, such as Figure 7 As shown, the model boundary is trimmed for the geological body model with feature fitting, including:

[0119] Step (1): Select seed point and initial edge.

[0120] Randomly select a seed point on the boundary line l. According to the arrangement order of the points on the boundary line, the side of the triangle formed by the seed point and the next point is used as the initial side. Each side of the triangle forming the boundary line is denoted as the boundary triangle side l. i {i=1,2,…,n}. The direction of the boundary line can be directly determined based on the side of the tangential strip that is retained.

[0121] Step (2): Determine the grid inside and outside the boundary.

[0122] Traverse along the direction of boundary line l, starting from the initial edge, and traverse the boundary triangle edge l. i The cross product of the direction vector and the normal vector is shown in equation (1), thus obtaining the side l of the boundary triangle. i Normal vector of the plane in which it lies:

[0123] n p,i =l′ i ×n i (1)

[0124] Among them, l′ i {i = 1, 2, ..., n} is the side l of the boundary triangle i The direction vector, n i {i = 1, 2, ..., n} is the side l of the boundary triangle i The normal vector, n p,i {i = 1, 2, ..., n} is the side l of the boundary triangle i The normal vector of the plane in which it is located, and it points to the inside of the model boundary.

[0125] Then, n p,i =(a i ,b i ,c i ) and l i A little bit above p i =(x i ,y i ,z i Substituting into the plane point normal form equation, we get l i The equation of the plane of the normal plane α is (2).

[0126] f p (x)=a i (xx i )+b i (yy i )+c i (zz i (2)

[0127] If the coordinates of the three vertices of a triangular facet are substituted into the function value f of equation (2) p If all x are greater than or equal to 0, then the triangular facet is inside the boundary. Otherwise, the triangular facet is outside the boundary.

[0128] Step (3): Mark the grid inside and outside the boundary.

[0129] For a triangular mesh of a topological manifold, a boundary triangle edge l i There are two adjacent triangular faces, denoted as t. lef and t rig To facilitate manipulation of the triangular facets, a label is added to each triangular facet, denoted as t. lab . t lab =1 indicates that the triangular facet will be retained, t lab =0 indicates that the triangular facet will be deleted. This will form t lef and t rig Substitute the grid vertex coordinates into equation (2). If the triangular facet t lef or t rig Substituting the three vertices into the function value f of equation (2) p If all x are greater than or equal to 0, then the label of the triangular facet is set to 1(t). lab =1). Otherwise, set the label of the triangular facet to 0 (t). lab =0). Then, set the adjacent unlabeled triangles to the same label. Next, with label t lab The triangles with a value of 1 are the primary focus. A breadth-first search (BFS) algorithm is used to traverse all unlabeled triangles. This will be compared with t. lab The label of the triangle facet with a value of 1 is set to 1 for its adjacent triangle facets, and the label of the triangle facet with a value of 1 is set to 1. lab The label of the triangle facet with a value of 0 is set to 0 for the adjacent triangle facets.

[0130] Step (4): Delete the mesh outside the boundary.

[0131] All tags t lab Deleting triangles with a value of 0 allows for boundary trimming of the geological body model. Furthermore, a geological body model may have more than one boundary; when multiple boundaries exist, the above method can be used simultaneously to trim the model boundaries.

[0132] It is understood that the geological body model construction method of this application performs mesh processing on the implicitly modeled geological body model based on geological constraint points and geological constraint lines, which enables the model to accurately fit the geological feature points and feature lines, thereby achieving feature fitting of the model. In addition, different markings are made on the mesh patches on both sides of the model boundary constraint lines, and the mesh outside the boundary constraint lines is deleted according to different label values ​​to accurately trim the model boundary, thereby achieving accurate trimming of the model boundary and obtaining a geological body model with matching boundaries.

[0133] Figure 8 This diagram illustrates the principle of conforming to geological constraint points in an application example. Figure 9 A schematic diagram illustrating the principle of fitting geological constraint lines and trimming model boundaries in an application example is shown. Experiments demonstrate that the method of this embodiment is automatic and practical in improving the precision of geological body models, and can effectively solve the problem that geological body mesh models obtained through implicit modeling techniques are not easily and accurately fitted to geological feature points and feature lines obtained from geological sampling data.

[0134] To implement the method of the embodiments of this application, the embodiments of this application also provide a geological body model construction device, which is installed in the geological body model construction equipment, such as... Figure 10 As shown, the geological body model construction device includes: a first fitting module 1001 and a second fitting module 1002. The first fitting module 1001 is used to perform model fitting processing on the initial geological body model constructed using implicit modeling technology based on geological constraint points to obtain a transition model; the second fitting module 1002 is used to perform model fitting processing on the transition model based on geological constraint lines to obtain a feature-fitted geological body model; wherein, the geological constraint points are directional interpolation points used for implicit surface interpolation, and the geological constraint points have tangential constraint directions and normal constraint directions; the geological constraint lines include: non-boundary constraint lines and boundary constraint lines, both of which have tangential constraint directions and normal constraint directions, the non-boundary constraint lines are non-boundary interpolation lines used for implicit surface interpolation, and the boundary constraint lines are boundary interpolation lines used for implicit surface interpolation.

[0135] In some embodiments, the first bonding module 1001 is specifically used for:

[0136] Obtain the constraint points of the initial geological body model within the set fitting range;

[0137] Calculate the distance from each constraint point to the mesh surface of the initial geological body model, and retain the constraint points whose distance is less than or equal to the maximum fitting distance as the geological constraint points;

[0138] The initial geological body model is subjected to model fitting processing by moving projection points to corresponding geological constraint points to obtain a transition model; wherein, the projection point is the point on the initial geological body model projected onto the corresponding geological constraint point, and the projection point was a mesh vertex of the initial geological body model before the movement.

[0139] In some embodiments, the second bonding module 1002 is specifically used for:

[0140] The transition model is subjected to model fitting processing based on the non-boundary constraint line and the boundary constraint line respectively to obtain a geological body model with feature fitting.

[0141] In some embodiments, the second bonding module 1002 performs model bonding processing on the transition model based on the boundary constraint lines, including:

[0142] Based on the boundary projection line projected onto the transition model surface by the boundary constraint line, a first projection neighborhood corresponding to the boundary constraint line is generated. The first projection neighborhood is the region where the set of topologically continuous triangular facets on the path where the boundary projection line is located is located.

[0143] Re-mesh the first projected neighborhood;

[0144] The remeshed triangular facets within the first projection neighborhood are merged with adjacent triangular facets, and the normals of the triangular facets are made uniform.

[0145] In some embodiments, the second bonding module 1002 performs model bonding processing on the transition model based on the non-boundary constraint lines, including:

[0146] Based on the non-boundary projection line projected onto the surface of the transition model, a second projection neighborhood corresponding to the non-boundary constraint line is generated. The second projection neighborhood is the region where the set of topologically continuous triangular facets on the path where the non-boundary projection line is located is located.

[0147] Remesh the second projection neighborhood;

[0148] Merge the remeshed triangular facets within the second projection neighborhood with adjacent triangular facets, and unify the normals of the triangular facets.

[0149] In some embodiments, the geological body model construction device further includes a boundary trimming module 1003, used to trim the model boundary of the feature-fitting geological body model to obtain a boundary-fitting geological body model.

[0150] In some embodiments, the boundary trimming module 1003 is specifically used for:

[0151] Along the direction of the boundary line of the geological body model that fits the features, start from the initial edge and traverse the edges of each boundary triangle. Based on the plane equation of the normal plane where each boundary triangle edge is located, determine whether the corresponding triangular facet is located within the model boundary.

[0152] Delete the triangular facets located outside the boundary of the model;

[0153] The initial edge is the edge formed by the selected starting point and the next adjacent point on the boundary line.

[0154] In some embodiments, the boundary trimming module 1003 determines whether the corresponding triangular facet is located within the model boundary based on the plane equation of the normal plane containing the sides of each boundary triangle, including:

[0155] If the coordinates of the three vertices of the triangular facet are substituted into the plane equation and the resulting function value is greater than or equal to zero, then the triangular facet is determined to be inside the model boundary; otherwise, the triangular facet is determined to be outside the model boundary.

[0156] In practical applications, the first bonding module 1001, the second bonding module 1002, and the boundary trimming module 1003 can be implemented by the processor in the geological body model construction device. Of course, the processor needs to run the computer program in the memory to realize its functions.

[0157] It should be noted that the geological body model construction device provided in the above embodiments is only illustrated by the division of the above-described program modules when constructing geological body models. In practical applications, the above processing can be assigned to different program modules as needed, that is, the internal structure of the device can be divided into different program modules to complete all or part of the processing described above. In addition, the geological body model construction device and the geological body model construction method embodiments provided in the above embodiments belong to the same concept, and their specific implementation process can be found in the method embodiments, which will not be repeated here.

[0158] Based on the hardware implementation of the above program modules, and in order to implement the method of the embodiments of this application, the embodiments of this application also provide a geological body model construction device. Figure 11 This is only an exemplary structure of the geological body model building device, not the entire structure; implementation is possible as needed. Figure 11 The structure shown may be part or all of the structure.

[0159] like Figure 11As shown, the geological body model building device 1100 provided in this application embodiment includes: at least one processor 1101, a memory 1102, a user interface 1103, and at least one network interface 1104. The various components in the geological body model building device 1100 are coupled together through a bus system 1105. It can be understood that the bus system 1105 is used to realize the connection and communication between these components. In addition to a data bus, the bus system 1105 also includes a power bus, a control bus, and a status signal bus. However, for clarity, in... Figure 11 The general labeled all buses as Bus System 1105.

[0160] The user interface 1103 may include a monitor, keyboard, mouse, trackball, click wheel, buttons, touchpad, or touch screen.

[0161] The memory 1102 in this embodiment is used to store various types of data to support the operation of the geological body modeling device. Examples of such data include any computer program used to operate on the geological body modeling device.

[0162] The geological body model construction method disclosed in this application embodiment can be applied to or implemented by processor 1101. Processor 1101 may be an integrated circuit chip with signal processing capabilities. During implementation, each step of the geological body model construction method can be completed by integrated logic circuits in the hardware of processor 1101 or by instructions in software form. The processor 1101 can be a general-purpose processor, a digital signal processor (DSP), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. Processor 1101 can implement or execute the methods, steps, and logic block diagrams disclosed in the embodiments of this application. A general-purpose processor can be a microprocessor or any conventional processor. The steps of the method disclosed in the embodiments of this application can be directly manifested as execution by a hardware decoding processor, or execution by a combination of hardware and software modules in the decoding processor. The software modules can be located in a storage medium, specifically memory 1102. Processor 1101 reads information from memory 1102 and, in conjunction with its hardware, completes the steps of the geological body model construction method provided in the embodiments of this application.

[0163] In an exemplary embodiment, the geological body model building device 1100 may be implemented by one or more application-specific integrated circuits (ASICs), DSPs, programmable logic devices (PLDs), complex programmable logic devices (CPLDs), FPGAs, general-purpose processors, controllers, microcontrollers (MCUs), microprocessors, or other electronic components to perform the aforementioned method.

[0164] It is understood that memory 1102 can be volatile memory or non-volatile memory, or both. Non-volatile memory can be read-only memory (ROM), programmable read-only memory (PROM), erasable programmable read-only memory (EPROM), electrically erasable programmable read-only memory (EEPROM), ferromagnetic random access memory (FRAM), flash memory, magnetic surface memory, optical disc, or compact disc read-only memory (CD-ROM); magnetic surface memory can be disk storage or magnetic tape storage. Volatile memory can be random access memory (RAM), which is used as an external cache. By way of example, but not limitation, many forms of RAM are available, such as Static Random Access Memory (SRAM), Synchronous Static Random Access Memory (SSRAM), Dynamic Random Access Memory (DRAM), Synchronous Dynamic Random Access Memory (SDRAM), Double Data Rate Synchronous Dynamic Random Access Memory (DDRSDRAM), Enhanced Synchronous Dynamic Random Access Memory (ESDRAM), SyncLink Dynamic Random Access Memory (SLDRAM), and Direct Rambus Random Access Memory (DRRAM).The memories described in the embodiments of this application are intended to include, but are not limited to, these and any other suitable types of memories.

[0165] In an exemplary embodiment, this application also provides a storage medium, namely a computer storage medium, specifically a computer-readable storage medium, such as a memory 1102 that stores a computer program. This computer program can be executed by the processor 1101 of the geological model construction device to complete the steps described in the method of this application embodiment. The computer-readable storage medium can be a ROM, PROM, EPROM, EEPROM, Flash Memory, magnetic surface memory, optical disc, or CD-ROM, etc.

[0166] It should be noted that terms such as "first" and "second" are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence.

[0167] Furthermore, the technical solutions described in the embodiments of this application can be combined arbitrarily without conflict.

[0168] The above description is merely a specific embodiment of this application, but the scope of protection of this application is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in this application should be included within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. A method for constructing a geological body model, characterized in that, include: The initial geological body model constructed using implicit modeling techniques is subjected to model fitting processing based on geological constraint points to obtain a transitional model; The transition model is subjected to model fitting processing based on geological constraint lines to obtain a geological body model with feature fitting; Wherein, the geological constraint points are directional interpolation points used for implicit surface interpolation, and the geological constraint points have tangential constraint directions and normal constraint directions; the geological constraint lines include: non-boundary constraint lines and boundary constraint lines, and both the non-boundary constraint lines and the boundary constraint lines have tangential constraint directions and normal constraint directions. The non-boundary constraint lines are non-boundary interpolation lines used for implicit surface interpolation, and the boundary constraint lines are boundary interpolation lines used for implicit surface interpolation. The process of fitting the transition model based on geological constraint lines to obtain a geological body model with feature fitting includes: The transition model is subjected to model fitting processing based on the non-boundary constraint line and the boundary constraint line respectively to obtain a geological body model with feature fitting; The transition model is subjected to model fitting processing based on the boundary constraint lines, including: Based on the boundary projection line projected onto the transition model surface by the boundary constraint line, a first projection neighborhood corresponding to the boundary constraint line is generated. The first projection neighborhood is the region where the set of topologically continuous triangular facets on the path where the boundary projection line is located is located. Re-mesh the first projected neighborhood; Merge the remeshed triangular facets in the first projection neighborhood with the adjacent triangular facets, and unify the normals of the triangular facets; The transition model is subjected to model fitting processing based on the non-boundary constraint lines, including: Based on the non-boundary projection line projected onto the surface of the transition model, a second projection neighborhood corresponding to the non-boundary constraint line is generated. The second projection neighborhood is the region where the set of topologically continuous triangular facets on the path where the non-boundary projection line is located is located. Remesh the second projection neighborhood; Merge the remeshed triangular facets within the second projection neighborhood with adjacent triangular facets, and unify the normals of the triangular facets.

2. The method according to claim 1, characterized in that, The process of fitting the initial geological body model constructed using implicit modeling techniques to geological constraint points yields a transitional model, including: Obtain the constraint points of the initial geological body model within the set fitting range; Calculate the distance from each constraint point to the mesh surface of the initial geological body model, and retain the constraint points whose distance is less than or equal to the maximum fitting distance as the geological constraint points; The initial geological body model is subjected to model fitting processing by moving projection points to corresponding geological constraint points to obtain a transition model; wherein, the projection point is the point on the initial geological body model projected onto the corresponding geological constraint point, and the projection point was a mesh vertex of the initial geological body model before the movement.

3. The method according to any one of claims 1 to 2, characterized in that, The method further includes: The geological body model with the described feature fit is trimmed to obtain a geological body model with matching boundaries.

4. The method according to claim 3, characterized in that, The step of trimming the model boundary of the geological body model that fits the features includes: Along the direction of the boundary line of the geological body model that fits the features, start from the initial edge and traverse the edges of each boundary triangle. Based on the plane equation of the normal plane where each boundary triangle edge is located, determine whether the corresponding triangular facet is located within the model boundary. Delete the triangular facets located outside the boundary of the model; The initial edge is the edge formed by the selected starting point and the next adjacent point on the boundary line.

5. A geological body model construction device, characterized in that, include: The first bonding module is used to bond the initial geological body model constructed using implicit modeling technology based on geological constraint points to obtain a transition model. The second bonding module is used to perform model bonding processing on the transition model based on geological constraint lines to obtain a geological body model with feature bonding. Wherein, the geological constraint points are directional interpolation points used for implicit surface interpolation, and the geological constraint points have tangential constraint directions and normal constraint directions; the geological constraint lines include: non-boundary constraint lines and boundary constraint lines, and both the non-boundary constraint lines and the boundary constraint lines have tangential constraint directions and normal constraint directions. The non-boundary constraint lines are non-boundary interpolation lines used for implicit surface interpolation, and the boundary constraint lines are boundary interpolation lines used for implicit surface interpolation. The process of fitting the transition model based on geological constraint lines to obtain a geological body model with feature fitting includes: The transition model is subjected to model fitting processing based on the non-boundary constraint line and the boundary constraint line respectively to obtain a geological body model with feature fitting; The transition model is subjected to model fitting processing based on the boundary constraint lines, including: Based on the boundary projection line projected onto the transition model surface by the boundary constraint line, a first projection neighborhood corresponding to the boundary constraint line is generated. The first projection neighborhood is the region where the set of topologically continuous triangular facets on the path where the boundary projection line is located is located. Re-mesh the first projected neighborhood; Merge the remeshed triangular facets in the first projection neighborhood with the adjacent triangular facets, and unify the normals of the triangular facets; The transition model is subjected to model fitting processing based on the non-boundary constraint lines, including: Based on the non-boundary projection line projected onto the surface of the transition model, a second projection neighborhood corresponding to the non-boundary constraint line is generated. The second projection neighborhood is the region where the set of topologically continuous triangular facets on the path where the non-boundary projection line is located is located. Remesh the second projection neighborhood; Merge the remeshed triangular facets within the second projection neighborhood with adjacent triangular facets, and unify the normals of the triangular facets.

6. A geological body model construction device, characterized in that, include: A processor and memory for storing computer programs that can run on the processor, wherein, The processor, when running a computer program, performs the steps of the method according to any one of claims 1 to 4.

7. A storage medium storing a computer program, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 4.

Images