A method and system for implicitly modeling geological bodies by fusing stratigraphic constraints

Abstract

The application discloses a kind of implicit modeling methods and systems of geologic body fusing stratum constraint rule.It is determined that the application method includes the following steps: defining the topological relationship and implicit potential field function of geological unit in three-dimensional space;Based on geological data, each geological interface isogonic surface expression is generated using implicit function interpolation method, and the geological interface is sequentially spliced;According to the geological conditions, constraint definition is carried out;Based on the definition of geological age, the direction of geological interface and geological sequence are defined;Determine the lithology type at any position in space, and then use the multi-domain implicit surface reconstruction method to extract the geological body isosurface for segmentation and reconstruction, to obtain the topologically consistent three-dimensional geologic body model.The application method is based on generalized radial basis function interpolation method and PFM potential field method implicit modeling, meets a variety of geological rule constraints, uses the extended surface reconstruction method to automatically generate high-precision three-dimensional geologic body model, greatly improves the efficiency and reliability of geologic body implicit modeling.

Description

TECHNICAL FIELD

[0001] The application belongs to the field of digital signal modeling, and particularly relates to a geological body implicit modeling method and system fusing stratum constraint rules. BACKGROUND

[0002] Based on sparse finite geological data, constructing a high-precision three-dimensional geological model is an important challenge faced by transparent geology research. Emerging implicit modeling technology, by introducing potential field method interpolation theory, converts geological interface representation into scalar field function spatial interpolation and contour surface extraction problem, providing a new idea for solving many challenges in traditional explicit modeling, and showing unique advantages in geological rule construction, model dynamic updating and geological grid generation.

[0003] Existing spatial interpolation algorithms cannot effectively fuse stratum occurrence, tectonic evolution and other geological prior knowledge, resulting in insufficient geological rationality of the constructed model, and a suitable method is needed to automatically or semi-automatically define geological rule constraints. SUMMARY

[0004] In view of the deficiencies of the prior art, one of the purposes of the present application is to provide a geological body implicit modeling method fusing stratum constraint rules, which improves the efficiency and reliability of geological body implicit modeling.

[0005] The second purpose of the present application is to provide a system for implementing the geological body implicit modeling method fusing stratum constraint rules.

[0006] The present application provides a geological body implicit modeling method fusing stratum constraint rules, comprising the following steps:

[0007] S1. Defining the topological relationship of geological units in three-dimensional space, and defining an implicit potential field function for each geological unit;

[0008] S2. Based on geological data, using implicit function interpolation method to generate the equipotential surface expression of each geological interface, and based on combination rules to sequentially splice the geological interfaces to generate a layered structure represented by implicit functions;

[0009] S3. According to the geological conditions, defining the constraints of each geological interface; and based on the geological age, defining the direction and geological sequence of the geological interface;

[0010] S4. According to the geological body lithology evaluation rule, determining the lithology type at any position in space, and using the multi-domain implicit surface reconstruction method to extract the geological body contour surface for segmentation and reconstruction, to obtain a topologically consistent three-dimensional geological body model.

[0011] In step S1, for any point in the 3D geological space, the implicit function corresponding to the geological interface is a certain equipotential surface, satisfying a point set a zero or non-zero level surface; a potential field corresponds to a geological interface modeling; the geological interface is selected as a geological contact surface with relatively consistent characteristics; each geological interface is defined by different lithology types and is interpolated by a spatial interpolation method; the geological interface has directionality, and the lithology types on both sides are used to determine the types of each geological body;

[0012] define a geological sequence of geological interfaces, which corresponds to a set of geological interfaces defined in chronological order; under the influence of geological activity, older strata are always located below newer strata;

[0013] define a geological action relationship, which is a combination rule of potential fields of each geological interface in the geological sequence, used to correctly reflect the chronological order of geological events; each geological body is composed of geological interfaces according to geological rules.

[0014] Step S2 is specifically: sequentially splicing the geological interfaces based on the combination rule to obtain an implicit function of a geological body model characterized by multiple combination potential fields , which is expressed using the following formula: ; wherein, is the i-th sub-potential field, including an implicitly interpolated function and other combination potential fields; is calculated based on the combination operation rule between potential fields;

[0015] use intersection potential field combination constraints, union potential field combination constraints, and difference potential field combination constraints to represent different geological events and their corresponding geological action relationships.

[0016] The intersection potential field combination constraint is used to construct an intersection operation of implicit surfaces, and the intersection potential field combination constraint composed of the intersection of potential field A and potential field B is expressed using the following formula: ; wherein, is the implicit function of potential field A; is the implicit function of potential field B; is the larger value of and ; and

[0017] The intersection potential field combination constraint represents that the combined implicit surface is composed of the common area of two implicit surfaces.

[0018] The union potential field combination constraint is used to construct a union operation of implicit surfaces, and the union potential field combination constraint composed of the union of potential field A and potential field B is expressed using the following formula: ; wherein, is and a smaller value of the two;

[0019] The union potential field combination constraint represents that the combined implicit surface is composed of two implicit surfaces together.

[0020] The difference potential field combination constraint is used to construct the difference set operation of the implicit surface, and the difference potential field combination constraint composed of the difference set of the potential field A and the potential field B is expressed by the following formula: ; wherein, is the difference set of the potential field A and the potential field B; is the complement set of the potential field B;

[0021] The difference set constraint represents that the combined implicit surface of a certain implicit surface is limited outside another implicit surface.

[0022] The step S3 is specifically:

[0023] According to the geological conditions, define the geological interface constraint; the geological interface constraint includes four basic types of sedimentary constraint, erosion constraint, intrusive constraint and vein constraint; the geological interface constraint does not directly correspond to a specific geological structure type, but is used to generate an extension boundary representing the contact surface between different lithological units, and define the geological action characteristics such as "avoidance" or "cutting" between the extension boundaries;

[0024] Define the direction of the geological interface based on the geological age; each geological interface has a younger side (Younger Side) and an older side (Older Side), and the relative age between the two sides determines how they interact and affect the division of geological bodies; the normal direction of the geological interface is towards the younger side area, so the direction of the geological interface can also be called the younger side direction;

[0025] Define the geological sequence based on the geological age; the geological sequence defines a series of geological interfaces, their corresponding geological age order and the lithology types on both sides of the contact surface; the geological sequence relationship not only clarifies the time sequence of the strata, but also reflects important events in the geological history, such as deposition, erosion, volcanic activity, etc.

[0026] The sedimentary constraint is used in the sedimentary geological interface, which presents a sheet or layer structure; according to the principle of geology, the sedimentary geological action feature is that the new stratum is stacked on the older stratum and does not cut off the older stratum; the sedimentary interpolation constraint feature is that the lithology on the older side and the younger side of the sedimentary contact surface should contact or avoid the older lithology, which is specifically:

[0027] Taking n strata as an example, the geological age relationship from new to old is , assuming that stratum 1 corresponds to the sedimentary geological interface, if the target lithology area is located in the younger side direction of the geological interface, it should avoid the older stacked stratum, and the corresponding combined field implicit function of the lithology area is The constraint can be represented by potential field combination as follows:

[0028]

[0029] Wherein, is an implicit function obtained by geological interface interpolation through generalized radial basis function;

[0030] The sedimentary geological interface of stratum 1 obtained by interpolation is The geological interface satisfying geological rules obtained by combination constraint is ;

[0031] If the target lithology area is located in the direction of the older side of the geological interface, the older stacked stratum should be removed, and the implicit function of the corresponding combination field of the lithology area is The constraint can be represented by potential field combination as follows: The sedimentary geological interface of stratum 1 obtained by interpolation is The geological interface satisfying geological rules obtained by combination constraint is .

[0032] The erosion constraint is used in the erosion-type geological interface, which presents a sheet or layer structure; the characteristic of erosion is that the older side of the stratum is stacked on the older stratum and the younger side of the stratum cuts the older stratum; the characteristic of erosion interpolation constraint is that the lithology on the older side of the erosion contact surface should “contact” or “avoid” the older lithology, and the lithology on the younger side of the erosion contact surface can “cut” or “truncate” the older lithology, which is specifically as follows:

[0033] Taking n strata as an example, the geological age relationship from new to old is Suppose stratum 1j corresponds to the erosion-type geological interface, if the target lithology area is located in the direction of the younger side of the geological interface, the older stacked stratum can be eroded, and the implicit function of the corresponding combination field of the lithology area can be represented by potential field combination constraint as follows: ; ;

[0034] The erosion-type geological interface of stratum 1j obtained by interpolation is The geological interface satisfying geological rules obtained by combination constraint is: ;

[0035] If the target lithology area is located in the direction of the older side of the geological interface, the older stacked stratum should be removed, and the implicit function of the corresponding combination field of the lithology area is The constraint can be represented by potential field combination as follows: ;

[0036] The geological interface satisfying geological rules obtained by combination constraint is: .

[0037] The intrusive constraint is used in the intrusive geological interface, which has a smooth and complex shape. The intrusive geological action is characterized by the fact that the intrusive body can cross multiple different strata and fill the intrusive lithology type. The intrusive interpolation constraint is characterized by the fact that the younger side of the intrusive contact surface replaces the older lithology, and the older side of the intrusive contact surface directly "contacts" other lithology;

[0038] Suppose there are three strata a, b, c, and the geological age relationship from new to old is Assume that stratum b corresponds to the intrusive geological interface. According to the cutting and inserting law, the intrusive lithology area is located in the younger side direction of the geological interface, replacing the older stacked strata. The corresponding implicit function of the lithology area is represented by the potential field combination constraint as follows:

[0039]

[0040] The intrusive geological interface of stratum b obtained by interpolation is The geological interface that satisfies the geological rules obtained by combination constraint is: ;

[0041] The vein-shaped constraint is used in the vein-shaped geological interface, which has a thin layer shape along the fissure or crack filling. The vein-shaped geological interface contact surface is composed of a clear hanging wall surface (Hanging Wall Surface) and a foot wall surface (Foot Wall Surface), and the average surface (Mean Surface) between the hanging wall surface and the foot wall surface can control the thickness and range of the rock vein;

[0042] Suppose there are three strata d, e, f, and the geological age relationship from new to old is Assume that stratum e corresponds to the vein-shaped geological interface. According to the cutting and inserting law, the intrusive or cutting geological age is newer. The vein lithology area is located in the younger side direction of the geological interface, which can "replace" the older stacked strata. The corresponding implicit function of the lithology area is represented by the potential field combination constraint as follows:

[0043]

[0044] The implicit function representing the vein-shaped geological interface is obtained by the following potential field combination constraint:

[0045]

[0046] wherein, is the implicit function representing the roof of the vein-shaped geological interface; The implicit function of the vein-shaped geological interface is represented; the direction of the roof and floor geological interface is consistent with the direction of the vein-shaped geological interface;

[0047] The vein-shaped geological interface of the stratum e is obtained by interpolation The geological interface satisfying the geological rules is obtained by combining the constraints .

[0048] The direction of the geological interface is defined based on the geological age, and the implicit function of the geological interface is defined using the following formula :

[0049]

[0050] wherein x is an arbitrary position point in the geological space; the orientation of the normal constraint should be consistent with the orientation of the geological interface, and the normal direction of the geological interface is oriented to the young side area, so the direction of the geological interface can also be referred to as the young side direction.

[0051] The geological sequence is defined based on the geological age, and the list of lithology and the list of geological interfaces in the geological space are determined according to the geological sampling data (such as drilling, outcrop, etc.); in addition to the type of the created geological interface, such as sedimentary type, erosion type, intrusive type or vein-shaped type, and the young side direction, the corresponding geological age sequence and the young side lithology and the old side lithology of the contact surface should also be specified.

[0052] According to the direction of the geological interface and the direction of the potential field, the lithology on both sides of the geological interface is divided into young side lithology and old side lithology.

[0053] In step S4, the lithology type at an arbitrary position in space is determined according to the geological body lithology evaluation rule, which is specifically:

[0054] For an arbitrary given point in the 3D geological space, the following three evaluation rules are defined to determine the lithology type at the position:

[0055] Rule one: for multiple sedimentary type geological interface constraints, according to the stratigraphic sequence rule, the lithology type at an arbitrary given position is determined in order according to the geological interface from old to young in the geological age sequence, and according to the lithology type (sedimentary side) in the old side direction of the geological interface;

[0056] Rule two: for multiple erosion type geological interface constraints, remove the lithology on the old side direction of each geological interface that is eroded according to the erosion type constraint rule, and then treat it as a sedimentary type geological interface constraint, and determine the lithology type at an arbitrary given position in order according to the rule one method;

[0057] Rule three: for multiple intrusive or vein type geological interface constraints, according to the cutting through law, in the order of geological age from young to old, according to the lithology type on the young side of the geological interface (intrusive side), the lithology type at any given position is determined in turn;

[0058] According to the different lithology types in the whole geological space, the multi-label implicit function is defined to represent the geological body, specifically:

[0059] Suppose N scalar functions characterize the corresponding geological interface in the 3D geological space, the lithology type on the young side of each geological interface is represented as , and the lithology type on the old side of each geological interface is represented as ;

[0060] For any given position , all scalar values and the corresponding lithology label at this position are calculated, and then the final scalar value and the lithology type of the geological body are calculated according to the lithology evaluation rule;

[0061] Suppose , then satisfies the following formula:

[0062]

[0063] wherein, wherein, j is the serial number of a certain lithology label determined by the lithology evaluation rule;

[0064] Finally, through the following two evaluation strategies, the is calculated, specifically:

[0065] Evaluation strategy one: use the outer package of the model to limit the modeling range of each geological interface, and regard the outer package of the model as an implicit function to accelerate the calculation of the implicit function;

[0066] Evaluation strategy two: by adjusting the calculation order of the scalar values at the given position , avoid calculating all scalar values, the adjustment method of the calculation order is as follows:

[0067] First, the function value representing the intrusive or vein type geological interface is determined in the order of geological age from young to old; then, the function value representing the sedimentary or erosion type geological interface is determined in the order of geological age from old to young; finally, according to the erosion type constraint rule, it is determined in turn whether the geological interface is "eroded".

[0068] In step S4, the method of multi-domain implicit surface reconstruction extracts the geological body contour surface for segmentation and reconstruction, specifically:

[0069] By extending the Marching Tetrahedra method, using the vertex clustering algorithm to improve the quality of mesh generation and the surface tracking strategy to improve the efficiency of isosurface extraction, a moving tetrahedron method suitable for multi-domain surface reconstruction is realized, and smooth and ladder-shaped isosurfaces of multiple geological bodies are extracted;

[0070] In order to ensure the topological consistency of different geological interfaces, the triangular mesh storage mechanism is extended, and each triangular facet is assigned a specific label on both sides, indicating the lithology type it belongs to; the triangular network generated by the extended Marching Tetrahedra method is topologically consistent and has no ambiguity;

[0071] In terms of multi-domain tetrahedral subdivision, tetrahedral triangular subdivision configurations for different strata are designed, there are 256 different triangular facet configurations for the four vertices of a tetrahedron, and there are only 4 basic triangular facet configurations in space through rotational symmetry; for the facet configuration with multiple geological boundaries, auxiliary points (such as the incenter of a triangle and the incenter of a tetrahedron) are introduced, and new triangularization rules are proposed to ensure that the triangular facets of adjacent tetrahedrons are continuous in space and approximate the geological interface, thereby ensuring that the boundaries between different geological regions can be correctly extracted, and the topological structure of the model remains consistent.

[0072] The present application also provides a system for implementing the implicit modeling method of the geological body according to the fusion stratum constraint rule, which comprises a geological unit definition module, a geological interface splicing module, a constraint definition and direction sequence definition module, and a surface segmentation reconstruction modeling module.

[0073] The geological unit definition module defines the topological relationship of the geological unit in the three-dimensional space, defines an implicit potential field function for each geological unit, and uploads the data to the geological interface splicing module.

[0074] The geological interface splicing module generates the equipotential surface expression of each geological interface using the implicit function interpolation method according to the received data, and then splices the geological interfaces in order based on the combination rule to generate a layered structure represented by an implicit function, and uploads the data to the constraint definition and direction sequence definition module.

[0075] The constraint definition and direction sequence definition module defines the constraints for each geological interface according to the received data and geological conditions, and then defines the direction and geological sequence of the geological interface based on the geological age, and uploads the data to the surface segmentation reconstruction modeling module.

[0076] The curved surface segmentation reconstruction modeling module determines the lithology type at any position in space according to a lithology evaluation rule of the geological body, and then extracts the geological body isosurface for segmentation reconstruction by using a multi-domain implicit curved surface reconstruction method, so as to obtain a topologically consistent three-dimensional geological body model.

[0077] The application discloses a geological body implicit modeling method and system fusing stratum constraint rules, and is based on a generalized radial basis function interpolation method and a PFM potential field method implicit modeling, satisfies multiple geological rule constraints, and automatically generates a high-precision three-dimensional geological body model by using an extended curved surface reconstruction method, so that the efficiency and reliability of the geological body implicit modeling are greatly improved.

[0078] The method allows a geological engineer to define constraint relationships between geological interfaces based on prior geological cognition rules of contact relationships (integrated contact, unconformable contact, sedimentary contact and intrusive contact, etc.) between strata and rock bodies, and ensures that the constructed geological model conforms to geological rules. Under the constraint of the geological rules, the geological body implicit modeling method can enable the geological engineer to pay more attention to understanding and interpretation of geological phenomena, instead of spending more effort to maintain the topological consistency of the geological body model. BRIEF DESCRIPTION OF DRAWINGS

[0079] Figure 1 A flowchart of the method of the application is shown in the figure;

[0080] Figure 2 A constraint effect diagram of the interpolation modeling of the sedimentary type geological interface when the method of the application is implemented is shown in the figure; wherein, Figure 2 A is a drilling data diagram; Figure 2 B is an interpolation constraint effect diagram; Figure 2 C is a modeling effect diagram;

[0081] Figure 3 A constraint effect diagram of the interpolation modeling of the erosive type geological interface when the method of the application is implemented is shown in the figure; Figure 3 A is a drilling data diagram; Figure 3 B is an interpolation constraint effect diagram; Figure 3 C is a modeling effect diagram;

[0082] Figure 4 A constraint effect diagram of the interpolation modeling of the intrusive type geological interface when the method of the application is implemented is shown in the figure; Figure 4 A is a drilling data diagram; Figure 4 B is an interpolation constraint effect diagram; Figure 4 C is a modeling effect diagram;

[0083] Figure 5 A constraint effect diagram of the interpolation modeling of the vein type geological interface when the method of the application is implemented is shown in the figure; Figure 5 A is a drilling data diagram; Figure 5 B is an interpolation constraint effect diagram; Figure 5 C is a modeling effect diagram;

[0084] Figure 6 Fig. 1 is a configuration effect diagram of a multi-domain tetrahedral triangulation patch in geological surface reconstruction when the method of the present application is implemented; Figure 6 A and Figure 6 Fig. 2 is a configuration effect diagram of a multi-domain tetrahedral triangulation patch in different location conditions containing two strata; Figure 6 Fig. 3 is a configuration effect diagram of a multi-domain tetrahedral triangulation patch in three strata; Figure 6 Fig. 4 is a configuration effect diagram of a multi-domain tetrahedral triangulation patch in four strata;

[0085] Figure 7 Fig. 5 is a structural schematic diagram of the system of the present application;

[0086] Figure 8 Fig. 6 is an instance diagram of geological interface interpolation in the embodiment of the present application; Figure 8 Fig. 7a is a lithology distribution map of a borehole; Figure 8 Fig. 7b is a geological interface interpolation map;

[0087] Figure 9 Fig. 8 is an instance diagram of geological body 1 modeling under different strata constraint rules in the embodiment of the present application; Figure 9 Fig. 8a is a geological interface interpolation map; Figure 9 Fig. 8b is a geological body reconstruction map;

[0088] Figure 10 Fig. 9 is an instance diagram of geological body 2 modeling under different strata constraint rules in the embodiment of the present application; Figure 10 Fig. 9a is a geological interface interpolation map; Figure 10 Fig. 9b is a geological body reconstruction map;

[0089] Figure 11 Fig. 10 is an instance diagram of geological body implicit modeling in the embodiment of the present application; wherein, Figure 11 Fig. 10a is an interpolation constraint effect diagram; Figure 11 Fig. 10b is a geological interface effect diagram; Figure 11 Fig. 10c is a partial geological body modeling effect diagram; Figure 11 Fig. 10d is a complete geological body modeling effect diagram. DETAILED DESCRIPTION

[0090] The present application provides a geological body implicit modeling method fusing strata constraint rules, a flowchart thereof is shown in Fig. 1, which comprises the following steps: Figure 1

[0091] S1. Defining the topological relationship of geological units in a three-dimensional space, and defining an implicit potential field function for each geological unit;

[0092] In step S1, for any point in a 3D geological space, the geological interface corresponds to a certain equipotential surface of the implicit function , which satisfies​ a set of points a zero or non-zero iso-surface; one potential field corresponds to one geological interface modeling; the geological interface is selected as the geological contact surface with relatively consistent characteristics; each geological interface is defined by different lithology types and is interpolated by a spatial interpolation method; the geological interface has directionality, and the lithology types on both sides are used to determine the types of each geological body;

[0093] define a geological sequence of geological interfaces, which corresponds to a set of geological interfaces defined in chronological order; under the influence of geological activity, older strata are always located below newer strata;

[0094] define a geological action relationship, which is a combination rule of potential fields of each geological interface in the geological sequence, used to correctly reflect the chronological order of geological events; each geological body is composed of geological interfaces according to geological rules.

[0095] S2. Based on the geological data, an implicit function interpolation method is used to generate the potential surface expression of each geological interface, and then the geological interfaces are sequentially spliced based on the combination rule to generate a layered structure represented by implicit functions;

[0096] Step S2 is specifically: sequentially splicing the geological interfaces based on the combination rule to obtain a geological body model implicit function represented by multiple combination potential fields , which is represented by the following formula: ; wherein, is the i-th sub-potential field, which includes directly interpolated implicit functions and other combination potential fields; is calculated based on the combination operation rule between potential fields;

[0097] The intersection potential field combination constraint, the union potential field combination constraint, and the difference potential field combination constraint are used to represent different geological events and their corresponding geological action relationships.

[0098] The intersection potential field combination constraint is used to construct an intersection operation of implicit surfaces, and the intersection potential field combination constraint composed of the intersection of potential field A and potential field B is represented by the following formula: ; wherein, is the implicit function of potential field A; is the implicit function of potential field B; is the larger value of and ; and

[0099] The intersection potential field combination constraint represents that the combined implicit surface is composed of the common area of two implicit surfaces;

[0100] The union potential field combination constraint is used to construct a union operation of implicit surfaces, and the union potential field combination constraint composed of the union of potential field A and potential field B is expressed by the following formula: ; wherein, is a smaller value of and .

[0101] The union potential field combination constraint represents that the combined implicit surface is composed of two implicit surfaces;

[0102] The difference potential field combination constraint is used to construct a difference operation of implicit surfaces, and the difference potential field combination constraint composed of the difference set of potential field A and potential field B is expressed by the following formula: ; wherein, is the difference set of potential field A and potential field B; is the complement set of potential field B;

[0103] The difference set constraint represents that the combined implicit surface is limited outside another implicit surface.

[0104] S3. According to the geological conditions, the geological interfaces are defined by constraints, and then the geological interface directions and geological sequences are defined based on the geological ages;

[0105] The step S3 is specifically:

[0106] According to the geological conditions, the geological interface constraints are defined; the geological interface constraints include four basic types of sedimentary constraints, erosion constraints, intrusion constraints and vein constraints; the geological interface constraints do not directly correspond to specific geological structure types, but are used to generate extension boundaries representing the contact surfaces between different lithological units and define the geological action characteristics such as “avoidance” or “cutting” between the extension boundaries;

[0107] The geological interface directions are defined based on the geological ages; each geological interface has a younger side and an older side, and the relative ages between the two sides determine how they interact and affect the division of geological bodies; the normal direction of the geological interface is towards the younger side area, so the direction of the geological interface can also be called the younger side direction;

[0108] The geological sequences are defined based on the geological ages; the geological sequence defines a series of geological interfaces, their corresponding chronological order and the lithological types on both sides of the contact surface; the geological sequence relationship not only clarifies the time sequence of strata, but also reflects important events in geological history, such as deposition, erosion, volcanic activity, etc.

[0109] The deposition type constraint is used in a deposition type geological interface, and the deposition type geological interface presents a sheet or layer structure; according to the principle of geology, the characteristic of deposition geological action is that new strata are stacked on older strata and do not cut off the older strata; the deposition interpolation constraint has the characteristics that the lithology on the old side and the young side of the deposition contact surface should contact or avoid the older lithology, and specifically,

[0110] Taking n strata as an example, the geological age relationship from new to old is Suppose that stratum 1 corresponds to a deposition type geological interface, if the target lithology region is located in the young side direction of the geological interface, the older stacked strata should be avoided, and the corresponding combined field implicit function of the lithology region is The combined constraint of the potential field can be expressed as:

[0111]

[0112] Wherein, is an implicit function obtained by interpolation of the geological interface through a generalized radial basis function;

[0113] The deposition type geological interface of stratum 1 obtained by interpolation is The geological interface obtained by the combined constraint to meet the geological rules is ;

[0114] If the target lithology region is located in the old side direction of the geological interface, the older stacked strata should be removed, and the corresponding combined field implicit function of the lithology region is The combined constraint of the potential field can be expressed as: The deposition type geological interface of stratum 1 obtained by interpolation is The geological interface obtained by the combined constraint to meet the geological rules is .

[0115] The deposition type geological interface interpolation modeling effect is shown in Figure 2 From the modeling constraint effect, the deposition type and the erosion type geological interface with good layering are generally near-horizontal or inclined three-dimensional open surfaces; generally, there is only one contact point on a single borehole, and for the geological interface with good layering, a 2.5-dimensional method can be considered for interpolation.

[0116] The erosion type constraint is used in an erosion type geological interface, and the erosion type geological interface presents a sheet or layer structure; the characteristic of erosion geological action is that the old side of the stratum is stacked on the older stratum and the young side of the stratum cuts the older stratum; the erosion interpolation constraint has the characteristics that the lithology on the old side of the erosion contact surface should “contact” or “avoid” the older lithology, and the lithology on the young side of the erosion contact surface can “cut” or “cut off” the older lithology, and specifically,

[0117] Taking n strata as an example, the geological age relationship from new to old is , assuming that the stratum 1j is a corresponding erosion-type geological interface, if the target lithology region is located in the young side direction of the geological interface, the older stacked stratum can be eroded, and the corresponding combination field implicit function of the lithology region can be constrained by potential field combination is expressed as: ;

[0118] The erosion-type geological interface of the stratum 1j obtained by interpolation is The geological interface satisfying the geological rules obtained by combination constraint is: ;

[0119] If the target lithology region is located in the old side direction of the geological interface, the older stacked stratum is “removed”, and the corresponding combination field implicit function of the lithology region can be constrained by potential field combination and expressed as: ;

[0120] The geological interface satisfying the geological rules obtained by combination constraint is: ;

[0121] The erosion-type geological interface interpolation modeling effect is shown in Figure 3 . From the modeling constraint effect, the geological interfaces obtained by deposition-type and erosion-type interpolation are completely consistent, but the geological action modes of the geological body model constructed by the geological interfaces are different. The deposition-type geological interface cannot cut the older stratum, and the erosion-type geological interface can cut the older stratum.

[0122] The intrusion-type constraint is used in the intrusion-type geological interface, and the intrusion-type geological interface has a smooth and complex shape. The intrusion geological action feature is that the intrusion body can cross multiple different strata and fill the intrusion lithology type. The intrusion interpolation constraint feature is that the younger side lithology of the intrusion contact surface “replaces” the older lithology, and the older side of the intrusion contact surface directly “contacts” other lithology.

[0123] Suppose there are three strata a, b, and c, and the geological age relationship from new to old is , assuming that the stratum b is a corresponding intrusion-type geological interface, according to the cutting insertion law, the intrusion lithology region is located in the young side direction of the geological interface, and replaces the older stacked stratum, and the corresponding combination field implicit function of the lithology region is constrained by potential field combination and expressed as:

[0124]

[0125] The intrusion-type geological interface of the stratum b obtained by interpolation is The geological interface satisfying the geological rules obtained by combination constraint is: ;

[0126] The intrusion-type geological interface interpolation modeling effect is shown inFigure 4 The intruding type and the vein type geological interface are generally three-dimensional closed surfaces with large changes in spatial form in terms of modeling constraint effects; multiple contact points can exist on a single borehole, and a true three-dimensional method is generally used for interpolation.

[0127] The vein type constraint is used in vein type geological interfaces, which have thin layer-shaped appearances filled along fissures or cracks; the vein type geological interface contact surface is composed of a clear hanging wall surface and a foot wall surface, and the rock vein thickness and range can be controlled through the mean surface between the hanging wall surface and the foot wall surface;

[0128] Suppose there are three strata d, e, and f, and the geological age relationship from new to old is Suppose stratum e is a corresponding vein type geological interface; according to the cutting and penetrating law, the intruder or the cutter is younger in geological age; the vein type lithology area is located in the young side direction of the geological interface and can “replace” the older stacked strata, and the corresponding implicit function of the lithology area combination field is obtained through the potential field combination constraint is expressed as:

[0129]

[0130] Implicit function representing the vein type geological interface is obtained from the following potential field combination constraint:

[0131]

[0132] wherein, is an implicit function representing the roof of the vein type geological interface; is an implicit function representing the floor of the vein type geological interface; the direction of the roof and floor geological interface is consistent with the direction of the vein type geological interface;

[0133] The vein type geological interface of stratum e obtained through interpolation is The geological interface satisfying the geological rules obtained through combination constraint is: .

[0134] The interpolation modeling effect of the vein type geological interface is shown in Figure 5 From the modeling constraint effect, the vein type constraint belongs to a geometric shape constraint for interpolating the shape of a specific geological body, so as to better control the extension and extrapolation trend of interpolation. The modeling process is similar to the modeling of thin layer type geological bodies. For vein type geological interfaces with small changes in spatial form, a 2.5-dimensional method can be used for interpolation, and for vein type geological interfaces with large changes in spatial form, a true three-dimensional method can be used for interpolation.

[0135] The geological interface direction is defined based on geological age, and an implicit function representing the geological interface is defined using the following formula :

[0136]

[0137] wherein x is an arbitrary position point in the geological space; the normal constraint direction should be consistent with the direction of the geological interface, and the normal direction of the geological interface is directed to the young side region, and thus the direction of the geological interface can also be referred to as the young side direction.

[0138] The geological sequence defined based on the geological age is specifically: a list of all lithologies participating in modeling in the geological space and a list of geological interfaces are determined according to geological sampling data (such as drill holes, outcrops, etc.); in addition to the type of the created geological interface, such as sedimentary type, erosion type, intrusive type or vein type, and the young side direction, the corresponding geological age sequence and the young side lithology and the old side lithology of the contact surface should also be indicated.

[0139] According to the direction of the geological interface and the direction of the potential field, the lithologies on both sides of the geological interface are divided into young side lithology and old side lithology.

[0140] S4. According to the geological body lithology evaluation rule, the lithology type of an arbitrary position in space is determined, and then a multi-domain implicit surface reconstruction method is used to extract the geological body isosurface for segmentation and reconstruction to obtain a topologically consistent three-dimensional geological body model.

[0141] In step S4, the lithology type of an arbitrary position in space is determined according to the geological body lithology evaluation rule, which is specifically:

[0142] For an arbitrary given point in the 3D geological space, the following three evaluation rules are defined to determine the lithology type at the position:

[0143] Rule one: for multiple sedimentary type geological interface constraints, according to the stratigraphic sequence rule, according to the geological age sequence of the geological interface from old to young, and according to the lithology type (sedimentary side) in the young side direction of the geological interface, the lithology type of an arbitrary given position is determined in sequence;

[0144] Rule two: for multiple erosion type geological interface constraints, according to the erosion type constraint rule, remove the lithology on the old side direction of each geological interface that is eroded, and then treat it as a sedimentary type geological interface constraint, and determine the lithology type of an arbitrary given position in sequence according to rule one;

[0145] Rule three: for multiple intrusive or vein type geological interface constraints, according to the cutting and penetrating law, according to the geological age sequence of the geological interface from young to old, and according to the lithology type (intrusion side) in the young side direction of the geological interface, the lithology type of an arbitrary given position is determined in sequence;

[0146] According to different lithology types of the whole geological space, a multi-label implicit function is defined to represent the geological body, specifically:

[0147] Suppose N scalar functions represent the corresponding geological interfaces in the 3D geological space, and the lithology type of the young side of each geological interface is represented as , and the lithology type of the old side of each geological interface is represented as ;

[0148] For any given position , all scalar values and corresponding lithology labels of the position are calculated, and then the final scalar value and lithology type of the geological body are calculated according to the lithology evaluation rule;

[0149] Suppose , then satisfies the following formula:

[0150]

[0151] wherein, j is the serial number of a certain lithology label determined by the lithology evaluation rule;

[0152] Finally, the following two evaluation strategies are used to calculate , specifically:

[0153] Evaluation strategy one: use the outer package of the model to limit the modeling range of each geological interface, and consider the outer package of the model as an implicit function to accelerate the calculation of the implicit function;

[0154] Evaluation strategy two: by adjusting the calculation order of the scalar values of the given position , avoid calculating all scalar values, and the adjustment method of the calculation order is as follows:

[0155] First, determine the function values representing the intrusive or vein type geological interface in the order of geological age from young to old; then, determine the function values representing the sedimentary or erosional type geological interface in the order of geological age from old to young; finally, determine whether the geological interface is "eroded" according to the erosional constraint rule.

[0156] In step S4, the method of multi-domain implicit surface reconstruction extracts the geological body isosurface for segmentation and reconstruction, specifically:

[0157] By extending the Marching Tetrahedra method, using vertex clustering algorithm to improve the quality of mesh generation and surface tracking strategy to improve the efficiency of isosurface extraction, a moving tetrahedron method suitable for multi-domain surface reconstruction is realized, and smooth and ladder-shaped isosurfaces of multiple geological bodies are extracted;

[0158] In order to ensure the topological consistency of different geological interfaces, the triangular mesh storage mechanism is extended, and each triangular facet is assigned a specific label on both sides, indicating the lithology type it belongs to; The triangular network generated by the extended Marching Tetrahedra method is topologically consistent and has no ambiguity.

[0159] In terms of multi-domain tetrahedral subdivision, tetrahedral triangular subdivision configurations for different strata are designed, and there are 256 different triangular facet configurations for the four vertices of a tetrahedron. There are only 4 basic triangular facet configurations in space through rotational symmetry, as shown in Figure 6 For the facet configuration with multiple geological boundaries, auxiliary points (such as the incenter of a triangle and the incenter of a tetrahedron) are introduced, and new triangularization rules are proposed to ensure that the triangular facets of adjacent tetrahedrons are continuous in space and approximate the geological interface, so that the boundaries between different geological regions can be correctly extracted, and the topological structure of the model remains consistent.

[0160] The present application also provides a system for implementing the implicit modeling method of the fused stratum constraint rule, as shown in Figure 7 The structure diagram is shown, which includes a geological unit definition module, a geological interface splicing module, a constraint definition and direction sequence definition module, and a surface segmentation reconstruction modeling module.

[0161] The geological unit definition module defines the topological relationship of the geological unit in the three-dimensional space, defines the implicit potential field function for each geological unit, and uploads the data to the geological interface splicing module.

[0162] The geological interface splicing module generates the equipotential surface expression of each geological interface using the implicit function interpolation method according to the received data, and then splices the geological interfaces in order based on the combination rule to generate a layered structure represented by an implicit function, and uploads the data to the constraint definition and direction sequence definition module.

[0163] The constraint definition and direction sequence definition module defines the constraints for each geological interface according to the received data and geological conditions, and then defines the direction and geological sequence of the geological interface based on the geological age, and uploads the data to the surface segmentation reconstruction modeling module.

[0164] The surface segmentation and reconstruction modeling module determines the lithology type at any spatial location based on the lithology assessment rules of the geological body. Then, it uses the multi-domain implicit surface reconstruction method to extract the isosurface of the geological body for segmentation and reconstruction, resulting in a topologically consistent three-dimensional geological body model.

[0165] The method of the present invention will be further described below with reference to an embodiment:

[0166] This embodiment implements an implicit geological body modeling algorithm that integrates multi-potential field geological rules using the C++ programming language, and develops an automatic modeling system for 3D geological bodies oriented towards transparent geology on a 3D visualization platform. To test the modeling effect of geological bodies under different geological rule constraints, multiple geological modeling instances are constructed to verify the reliability of the implicit geological body modeling method based on generalized radial basis function interpolation.

[0167] To simplify the analysis of the multipotential field geological modeling process, four simplified geological modeling processes for strata and their corresponding geological rule construction methods were constructed. Figure 8 Four spatial interpolation examples of geological interfaces are presented. Based on the contact point data of different lithologies in the borehole, the surface level of each stratum is obtained by using the generalized radial basis function interpolation method.

[0168] Figure 9 and Figure 10 This demonstrates the modeling effects of geological bodies under different stratigraphic constraints. The geological interfaces are categorized as A, B, C, and D based on their geological age from youngest to oldest. Figure 9 The four stratigraphic constraint types are: A is sedimentary strata, B is sedimentary strata, C is sedimentary strata, and D is sedimentary strata. Figure 10 The four stratigraphic constraint types are: A is an erosional stratum, B is a sedimentary stratum, C is a sedimentary stratum, and D is a sedimentary stratum.

[0169] To further verify the transparent geological implicit modeling technology proposed in this invention, three-dimensional modeling analysis was conducted using the modeling of complex geological bodies in several different mining areas as examples, such as... Figure 11 As shown.

Claims

1. A method for implicitly modeling geological bodies by fusing stratigraphic constraints, characterized in that, The method comprises the following steps: S1. Defining the topological relationship of the geological units in a three-dimensional space, and defining an implicit potential field function for each geological unit; S2. Generating the equipotential surface expression of each geological interface by using an implicit function interpolation method based on geological data, and then sequentially splicing the geological interfaces based on combination rules to generate a layered structure represented by implicit functions; the implicit function interpolation method is to perform geological interface interpolation by using a generalized radial basis function; S3. Defining the constraints of each geological interface according to geological conditions; and then defining the direction and geological sequence of the geological interface based on the geological age; S4. Determining the lithology type at any position in space according to a geological body lithology evaluation rule, and then extracting the geological body isosurface by using a multi-domain implicit surface reconstruction method to perform segmentation and reconstruction, so as to obtain a three-dimensional geological body model with consistent topology; Step S3 is specifically: Defining geological interface constraints according to geological conditions; the geological interface constraints include four basic types of sedimentary constraints, erosion constraints, intrusive constraints and vein constraints; the geological interface constraints do not directly correspond to a specific geological structure type, but are used to generate an extension boundary representing the contact surface between different lithology units, and define the “avoidance” or “cutting” geological action characteristics between the extension boundaries; Defining the direction of the geological interface based on the geological age; each geological interface has a young side and an old side, and the relative age between the two sides determines how they interact and affect the division of the geological body; the normal direction of the geological interface is towards the young side area, so the direction of the geological interface can also be referred to as the young side direction; Defining the geological sequence based on the geological age; the geological sequence defines a series of geological interfaces, their corresponding geological age order and the lithology types on both sides of the contact surface; the geological sequence relationship not only clearly defines the time sequence of the strata, but also reflects important events in the geological history, including deposition, erosion and volcanic activity.

2. The method of implicitly modeling geobodies with fusion stratigraphic constraints rules according to claim 1, wherein, In step S1, for any point in the 3D geological space , the geological interface corresponds to a certain equipotential surface of the implicit function , satisfying The point set constitutes a zero or non-zero contour surface; one potential field corresponds to the modeling of one geological interface; the geological interface is selected as a geological contact surface with relatively consistent characteristics; each geological interface is defined by different lithology types and is obtained by interpolation using a spatial interpolation method; the geological interface has directionality, and the lithology types on both sides are used to determine the types of each geological body; Defining the geological sequence of the geological interface, the geological sequence corresponds to a set of geological interfaces defined in the order of geological time; without the influence of geological activities, the older stratum is always located below the newer stratum; Defining the geological action relationship, the geological action relationship is the potential field combination rule of each geological interface in the geological sequence, which is used to correctly reflect the chronological order of geological events; each geological body is composed of geological interfaces according to the geological rules.

3. The method of implicitly modeling geobodies with fusion stratigraphic constraints rules of claim 1, wherein, The step S2 specifically comprises: sequentially splicing the geological interfaces based on a combination rule to obtain an implicit function of a geological body model characterized based on multiple combination potential fields , and is expressed by the following formula: ; wherein, is the i-th sub-potential field, the sub-potential field comprising an implicitly function of direct interpolation and other combination potential fields; is calculated based on a combination operation rule between potential fields. The intersection potential field combination constraint, the union potential field combination constraint and the difference set potential field combination constraint are used to represent different geological events and their corresponding geological action relationships.

4. The method of implicitly modeling geobodies with strata-constrained rules of fusion according to claim 3, wherein, The intersection potential field combination constraint is used for constructing an intersection operation operation of an implicit surface, and the intersection potential field combination constraint composed of the intersection of the potential field A and the potential field B is expressed by the following formula: ; wherein, is an implicit function of the potential field A; is an implicit function of the potential field B; is the maximum value of and . The intersection potential field combination constraint represents that the combined implicit surface is composed of the common area of two implicit surfaces; The union potential field combination constraint is used for constructing the union operation of the implicit surface, and the union potential field combination constraint composed of the union of the potential field A and the potential field B is expressed by using the following formula: ; wherein, is the smaller value of and . The union potential field combination constraint represents that the combined implicit surface is composed of two implicit surfaces together; The difference set potential field combination constraint is used for constructing a difference set operation operation of an implicit surface, and the difference set potential field combination constraint composed of a difference set of the potential field A and the potential field B is expressed by using the following formula: ; wherein, is a difference set of the potential field A and the potential field B; is a complement set of the potential field B; The difference set constraint represents that the combined implicit surface of a certain implicit surface is limited outside another implicit surface.

5. The method of implicitly modeling geobodies with stratigraphic constraints of claim 1, wherein, The sedimentary constraint is used in a sedimentary geological interface, and the sedimentary geological interface presents a sheet or layer structure; according to the principle of geology, the sedimentary geological action feature is that the new stratum is stacked on the older stratum and does not cut off the older stratum; the sedimentary interpolation constraint feature is that the old side and the young side of the sedimentary contact surface contact or avoid the older lithology, specifically: Suppose there are n strata, and the geological age relationship from new to old is , assuming that stratum 1 is the corresponding sedimentary geological interface, if the target lithology area is located on the young side of the geological interface, avoiding the older stacked strata, the implicit function of the corresponding combination field of this lithology area can be expressed as: wherein, is an implicit function obtained by geological interface interpolation by generalized radial basis functions; The sedimentary geological interface of stratum 1 obtained by interpolation is The geological interface satisfying geological rules obtained by combination constraint is ; If the target lithology area is located in the old side of the geological interface, remove the older stacked strata, and the corresponding combination field implicit function of the lithology area is The constraint representation can be obtained by combining the potential field as: The geological interface that satisfies the geological rules obtained by combining the constraints is ; The erosion type constraint is used in an erosion type geological interface, which presents a sheet or layer structure; the erosion geological action feature is that the older stratum side is stacked on the older stratum and the younger stratum side cuts the older stratum; the erosion interpolation constraint feature is that the older side of the erosion contact surface is "contacted" or "avoided" by the older lithology, and the younger side of the erosion contact surface can "cut" or "truncate" the older lithology, specifically: Suppose there are n strata, and the geological age relationship from new to old is , assuming that stratum 1j is a corresponding erosion type geological interface, if the target lithology region is located in the young side direction of the geological interface, the older stacked strata can be eroded, and the corresponding combination field implicit function of the lithology region can be constrained by the potential field combination is expressed as: ; The interpolation erosion type geological interface of the stratum 1j is The geological interface satisfying the geological rules obtained by combining the constraints is: ; If the target lithology region is located in the old side direction of the geological interface, the corresponding combination field implicit function of the lithology region is obtained by "removing" the older stacked strata The constraint representation can be obtained by potential field combination as: ; The geological interface satisfying the geological rule by combining the constraints is: ; The intrusion type constraint is used in an intrusion type geological interface, which has a smooth and complex shape; the intrusion geological action feature is that the intrusion body can cross multiple different strata and fill the intrusion body lithology type; the intrusion interpolation constraint feature is that the younger side of the intrusion contact surface "replaces" the older lithology, and the older side of the intrusion contact surface directly "contacts" other lithology; Suppose there are three strata a, b, c, and the geological age relationship from new to old is Assume that stratum b is the corresponding intrusive geological interface. According to the cutting insertion law, the intrusive lithology area is located on the young side of the geological interface, replacing the older stacked strata. The corresponding combination field implicit function of this lithology area passes through the potential field combination constraint It is represented as: The interpolation of the intrusive geological interface of the stratum B is The geological interface satisfying the geological rules obtained by combining constraints is: ; The vein type constraint is used in a vein type geological interface, which has a thin layer shape filled along a fissure or crack; the vein type geological interface contact surface is composed of a clear hanging wall and footwall, and the vein thickness and range can be controlled by the average surface between the hanging wall and footwall; Suppose there are three strata d, e, f, and the geological age relationship from new to old is , assuming that stratum e is a corresponding vein type geological interface; according to the cutting through law, the intruder or the cutter is younger in geological age; the vein type lithology area is located in the young side direction of the geological interface, and can "replace" the older stacked strata, and the corresponding combination field implicit function of the lithology area passes through the potential field combination constraint is expressed as: Implicit function for characterizing a vein-type geological interface is obtained from the following potential field combination constraints: wherein, is an implicit function representing the top plate of the vein type geological interface; is an implicit function representing the bottom plate of the vein type geological interface; the direction of the top-bottom plate geological interface is consistent with the direction of the vein type geological interface; The geological interface of the vein type obtained by interpolation for the stratum e is The geological interface satisfying the geological rules obtained by combining constraints is .

6. The method of implicitly modeling geobodies with stratigraphic constraints of claim 5, wherein, The geological interface direction is defined based on the geological time, and an implicit function representing the geological interface is defined using the following formula : Wherein, x is an arbitrary position point in the geological space; the orientation of the normal constraint should be consistent with the orientation of the geological interface, and the normal direction of the geological interface is oriented to the younger side area, so the direction of the geological interface can also be called the younger side direction; The geological sequence is defined based on the geological age, specifically: determining a list of all lithologies participating in modeling and a list of geological interfaces in the geological space according to geological sampling data; the created geological interface, in addition to the specified deposition type, erosion type, intrusion type or vein type and younger side direction, also indicates the corresponding geological age sequence and the contact surface younger side lithology and older side lithology; According to the direction of the geological interface and the direction of the potential field, the lithology on both sides of the geological interface is divided into younger side lithology and older side lithology.

7. The method of implicitly modeling geobodies with stratigraphic constraints of claim 1, wherein, In step S4, the lithology type at an arbitrary position in space is determined according to the geological body lithology evaluation rule, specifically: For any given point in 3D geological space The following three evaluation rules are defined to determine the lithology type at that location: Rule one: for multiple deposition type geological interface constraints, according to the stratum sequence rule, the lithology type at an arbitrary given position is determined in order according to the geological age sequence of the geological interface from old to young and the lithology type in the direction of the older side of the geological interface; Rule two: for multiple erosion type geological interface constraints, remove the lithology in the direction of the older side of each geological interface that is eroded according to the erosion type constraint rule, and then treat it as a deposition type geological interface constraint, and determine the lithology type at an arbitrary given position in order according to rule one; Rule three: for multiple intrusion type or vein type geological interface constraints, according to the cut-through law, the lithology type at an arbitrary given position is determined in order according to the geological age sequence of the geological interface from young to old and the lithology type in the direction of the younger side of the geological interface; Then, according to different lithology types in the whole geological space, a multi-label implicit function is defined to represent the geological body, specifically: Assume that there are N scalar functions representing the respective geologic interfaces within the 3D geologic space, the young side lithology type of each geologic interface is represented as , and the old side lithology type of each geologic interface is represented as ; For any given location , all scalar values and corresponding lithology labels at that location are computed and the final geobody's scalar value and lithology type are derived according to the lithology evaluation rules Assume then satisfies the following equation: Wherein, j is the serial number of a certain type of lithology label determined by the lithology evaluation rule; Finally, the following two evaluation strategies are used to calculate , in particular: Evaluation strategy one: use the model outer package to limit the modeling range of each geological interface, and treat the model outer package as an implicit function to accelerate the calculation of the implicit function; Evaluation strategy two: Avoid computing all scalar values by adjusting the order of computation of the scalar values for a given position The adjustment method of the order of computation is as follows: First, the function value representing the intrusive or vein type geological interface is determined in the geological interface from young to old geological age sequence; then, the function value representing the sedimentary or erosion type geological interface is determined in the geological interface from old to young geological age sequence; finally, whether the geological interface is eroded is judged according to the erosion type constraint rule in turn.

8. The method of implicitly modeling geobodies with stratigraphic constraints of claim 1, wherein, In step S4, the method of multi-domain implicit surface reconstruction is used to extract the geological body isosurface for segmentation and reconstruction, specifically: By extending the Marching Tetrahedra moving tetrahedron method, using the vertex clustering algorithm to improve the grid generation quality and the surface tracking strategy to improve the isosurface extraction efficiency, a moving tetrahedron method suitable for multi-domain surface reconstruction is realized, and smooth and ladder-free isosurfaces of multiple geological bodies are extracted; In order to ensure the topological consistency of different geological interfaces, the triangular mesh storage mechanism is extended, and each triangular facet is assigned a specific label on both sides, indicating the lithology type it belongs to; the triangular network generated by the extended Marching Tetrahedra moving tetrahedron method is topologically consistent and has no ambiguity; In terms of multi-domain tetrahedral subdivision, tetrahedral triangular subdivision configurations for different strata are designed, there are 256 different triangular facet configurations for the four vertices of a tetrahedron, and there are only 4 basic triangular facet configurations in space through rotational symmetry; for the facet configuration with multiple geological boundaries, auxiliary points are introduced, and new triangularization rules are proposed to ensure that the triangular facets of adjacent tetrahedrons are continuous in space and approximate the geological interface, so that the boundaries between different geological regions can be correctly extracted, and the topological structure of the model remains consistent.

9. A system for implementing the method of implicitly modeling geological bodies subject to the rules of fusion of strata constraints according to any one of claims 1 to 8, characterized in that it comprises: The system comprises a geological unit definition module, a geological interface splicing module, a constraint definition and direction sequence definition module, and a surface segmentation and reconstruction modeling module. The geological unit definition module defines the topological relationship of the geological unit in the three-dimensional space, defines an implicit potential field function for each geological unit, and uploads the data to the geological interface splicing module. The geological interface splicing module generates the isosurface expression of each geological interface using the implicit function interpolation method based on the received data, and then splices the geological interfaces in order based on the combination rules to generate a layered structure represented by implicit functions, and uploads the data to the constraint definition and direction sequence definition module. The constraint definition and direction sequence definition module defines the constraints for each geological interface according to the received data and geological conditions, and then defines the direction and geological sequence of the geological interface based on the geological age, and uploads the data to the surface segmentation and reconstruction modeling module. The surface segmentation and reconstruction modeling module determines the lithology type at any position in space according to the geological body lithology evaluation rule, and then uses the method of multi-domain implicit surface reconstruction to extract the geological body isosurface for segmentation and reconstruction to obtain a topologically consistent three-dimensional geological body model.

Images