Geological interface fitting method based on space rotation, medium, equipment and product

Abstract

The invention provides a geological interface fitting method based on spatial rotation, and relates to the technical field of geological information, and the method comprises the steps: obtaining stratum bottom boundary data of different years in a three-dimensional geological model; identifying boundary lines with a multi-valued function relationship from the boundary line data; performing coordinate rotation transformation according to the inclination direction and the inclination angle of the boundary with the multi-valued function relationship in the space; if the boundary line after rotation still has the multi-valued function relationship, adjusting the rotation axis and / or the angle until the boundary line after rotation is a single value or reaches the set maximum number of iterations; in the final flattened coordinate system, fitting data points of the final boundary line, generating a geological interface curved surface and optimizing the geological interface curved surface; and reversely rotating the optimized curved surface and the final boundary from the final flattened coordinate system to the original coordinate system to serve as a shear surface, and segmenting three-dimensional entity models of stratums of different years. According to the method, the modeling problem of the geological structure with the multi-valued function relation can be efficiently solved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of geological information, and in particular to a geological interface fitting method based on spatial rotation, a medium, an equipment and a product. BACKGROUND

[0002] In the process of three-dimensional geological structure modeling, the conventional spatial interpolation method is based on the single value assumption, that is, the target surface is a single value function on the projection plane (usually XY plane) by default or requirement, that is, for a given (X, Y) coordinate, there is only one Z value. For the stratum boundary with complex folds or faults, it may present a multi-value function relationship on the horizontal plane, that is, there are multiple elevation values corresponding to at least one plane position. When there are multi-value points (for example, drilling through the same rock layer twice of an inverted fold at the same position on the ground surface), these points will be regarded as two different Z values at the same (X, Y) position by the interpolation algorithm. This violates the single value premise of the algorithm, and usually leads to: algorithm error, crash; generating distorted, self-intersecting, unreasonable surfaces; only one Z value (such as the average value, the maximum value, the minimum value) can be forced to be selected, and the key structure information is lost. This multi-value relationship leads to the inability to directly generate a geological interface surface by using the conventional interpolation method.

[0003] The existing inverted fold modeling method (such as the inverted fold three-dimensional geological modeling controlled by "keel line") can handle simple folded stratum conditions, but it is difficult to work when facing complex geological structures with multiple intersecting profile stratum lines. By using surface discretization and partition modeling, a single, globally continuous mathematical surface is abandoned, and instead, a discrete triangular mesh or regular mesh is used to represent the interface, and the area containing multi-value points is divided into different "branches" or "patches". The surface discretization and partition modeling highly depend on the experience and subjective judgment of the modeler, and have high artificial cost, relatively low efficiency and limited automation. The core idea of implicit modeling is that the surface Z=f(X, Y) is not directly represented, but a three-dimensional scalar field function F(X, Y, Z) is defined. The target geological interface is represented as a point set satisfying a certain equal value condition (usually F(X, Y, Z)=0). This method allows one (X, Y) position to correspond to multiple Z values. The surface shape generated by implicit modeling is sometimes difficult to directly correspond to the intuitive geological interpretation of geologists, and additional visualization tools are needed, and the parameter tuning is complex, the computational overhead is large, and the shape control of the surface at a specific position (such as a fault) is not accurate enough. Therefore, there is an urgent need for a three-dimensional fitting method capable of processing multi-value function relationship geological interface to overcome the limitations of the prior art and improve the accuracy and efficiency of complex geological structure model construction. SUMMARY

[0004] The application aims to solve the problems of low accuracy and efficiency and limited automation of multi-value function geological interface modeling, and provides a geological interface fitting method based on space rotation, comprising the following steps:

[0005] S1, obtaining stratum bottom boundary data of different ages in a three-dimensional geological model from geological exploration data of a target area;

[0006] S2, identifying stratum bottom boundary data with multi-value function relationship from the stratum bottom boundary data, wherein the multi-value function relationship refers to that one plane coordinate point in a three-dimensional space corresponds to multiple elevation values of a stratum boundary;

[0007] S3, determining the tilt direction and tilt angle of the stratum bottom boundary data with multi-value function relationship in space, selecting a rotation axis according to the tilt direction, and calculating an initial rotation angle value according to the tilt angle;

[0008] S4, performing coordinate transformation on the stratum bottom boundary data with multi-value function relationship around the rotation axis according to the initial rotation angle value to obtain a rotated boundary;

[0009] S5, judging whether the rotated boundary still has multi-value function relationship, if yes, adjusting the rotation axis and / or angle according to the residual tilt condition, repeating steps S4 and S5 until the rotated boundary does not have multi-value function relationship or reaches a set maximum iteration number, and obtaining a final boundary in a final flattened coordinate system;

[0010] S6, fitting data points of the final boundary in the final flattened coordinate system to generate a geological interface surface, and optimizing the geological interface surface;

[0011] S7, inversely rotating the optimized geological interface surface and the final boundary from the final flattened coordinate system back to an original coordinate system;

[0012] S8, taking the geological interface surface restored to the original coordinate system as a shear surface to shear the three-dimensional geological model to segment a three-dimensional entity model of strata of different ages.

[0013] Further, the selection of the rotation axis according to the tilt direction is as follows: if the boundary tilts northward or southward, rotating around the X axis is selected; if the boundary tilts eastward or westward, rotating around the Y axis is selected; if the tilt direction of the boundary does not belong to any of the northward, southward, eastward and westward, rotating around the X axis and the Y axis in sequence is selected.

[0014] Further, the tilt angle is an angle of the stratum bottom boundary data with multi-value function relationship relative to a horizontal plane.

[0015] Further,

[0016] The coordinate transformation relationship when rotating around the X axis by an angle a is as follows:

[0017] [x′ y′ z′ 1]=[x y z 1]·R x (α)

[0018]

[0019] The coordinate transformation when rotating the Y axis by an angle β is as follows:

[0020] [x″ y″ z″ 1]=[x y z 1]·R y (β)

[0021]

[0022] Wherein, (x, y, z) represents the coordinates before rotation, (x', y', z') represents the coordinates after rotating the X axis by an angle α, (x'', y'', z'') represents the coordinates after rotating the Y axis by an angle β, α and β are rotation angles, R x (α) represents the rotation matrix of rotating the X axis by an angle α, R y (β) represents the rotation matrix of rotating the Y axis by an angle β.

[0023] Further, the Kriging interpolation algorithm is used to fit the data points of the final boundary, and the discrete point data of the final boundary is spatially interpolated into a continuous elevation surface, and an irregular triangular mesh geological interface surface is generated through triangulation.

[0024] Further, after the three-dimensional entity model of different ages of strata is segmented, the three-dimensional entity model of different ages of strata is verified for closure, and open boundaries generated by cutting are checked and repaired.

[0025] The application further provides a computer readable storage medium, the computer readable storage medium stores a computer program, and the computer program is executed by a processor to realize the above-mentioned geological interface fitting method based on spatial rotation.

[0026] The application further provides an electronic device, including a processor and a memory, the processor and the memory are connected with each other, wherein the memory is used for storing a computer program, the computer program includes computer readable instructions, the processor is configured to call the computer readable instructions, and the above-mentioned geological interface fitting method based on spatial rotation is executed.

[0027] The application further provides a computer program product, including computer program / instructions, which are executed by a processor to realize the steps of the above-mentioned geological interface fitting method based on spatial rotation.

[0028] The technical scheme provided by the application has the beneficial effects that:

[0029] The present application can generate the corresponding triangular mesh surface by using interpolation algorithm by rotating the stratum bottom boundary of complex geological interface in the spatial coordinate system around X axis or Y axis by a proper angle to project it as a single-valued function relationship, and then restoring the surface to the original coordinate system and cutting the initial geological body model to extract the three-dimensional entity of strata of different ages layer by layer. The method of the present application can efficiently process the modeling problem of complex geological structure with multiple-valued function relationship, overcome the bottleneck of being unable to fit a surface due to at least one plane coordinate corresponding to multiple elevation values, greatly improve the efficiency and accuracy of three-dimensional geological model construction, and generate a more refined model that fits the real geological conditions with high automation. BRIEF DESCRIPTION OF DRAWINGS

[0030] Figure 1 is a flowchart of a geological interface fitting method based on spatial rotation according to an embodiment of the present application;

[0031] Figure 2 is a schematic diagram of multiple sets of stratum bottom boundary lines of different ages of complex geological structure in a series of geological profiles according to an embodiment of the present application;

[0032] Figure 3 is a comparison diagram of a stratum bottom boundary line of a certain age before and after being processed and rotated and interpolated to generate a shear surface according to an embodiment of the present application;

[0033] Figure 4 is a refined three-dimensional geological model effect diagram of complex geological structure according to an embodiment of the present application;

[0034] Figure 5 is an internal structure fence diagram of a three-dimensional geological structure model according to an embodiment of the present application;

[0035] Figure 6 is a block diagram of an electronic device in an exemplary embodiment according to an embodiment of the present application. DETAILED DESCRIPTION

[0036] To make the objects, technical solutions and advantages of the present application clearer, the embodiments of the present application will be further described below with reference to the drawings.

[0037] The flowchart of the geological interface fitting method based on spatial rotation according to an embodiment of the present application is as shown in Figure 1 , which specifically includes the following steps:

[0038] S1, obtain geological exploration data of a target area, construct a three-dimensional geological model of the target area based on the geological exploration data, and obtain stratum bottom boundary line data of different ages in the three-dimensional geological model from a series of geological profiles of the three-dimensional geological model. The stratum bottom boundary lines of different ages are sorted and classified respectively, for example, stored in independent layers and marked with numbers for subsequent processing.

[0039] S2, identify the existence of a multi-value function relationship from the bottom boundary data of the formation, the multi-value function relationship refers to a plane coordinate point corresponding to multiple elevation values in a three-dimensional space, which usually occurs in complex areas with developed folds or faults.

[0040] S3, determine the tilt direction and tilt angle of the bottom boundary of the formation with multi-value function relationship in space, select the rotation axis according to the tilt direction, and calculate the initial rotation angle value according to the tilt angle.

[0041] Determine the tilt direction of the bottom boundary of the formation in space, that is, determine the main direction of the boundary. According to the tilt direction of the boundary, the rotation axis is selected: if the boundary tilts north or south as a whole, the X-axis rotation is selected; if it tilts east or west, the Y-axis rotation is selected; if the tilt direction of the boundary does not belong to any of the north, south, east, and west (such as the southeast direction), it is recorded that the X-axis and Y-axis need to be rotated in turn. The positive direction of rotation follows the right-hand rule: for example, from the positive direction of the X-axis to the origin, the boundary rotates counterclockwise to the positive angle of the X-axis.

[0042] The determination of the rotation angle is based on the spatial distribution characteristics of the bottom boundary of the formation, and the average tilt angle of the boundary relative to the horizontal plane is usually taken as the initial value of the rotation angle.

[0043] S4, the bottom boundary of the formation with multi-value function relationship is transformed around the rotation axis according to the initial rotation angle value, and the rotated boundary is obtained.

[0044] The coordinate transformation relationship when rotating α angle around the X-axis is as follows:

[0045] [x′ y′ z′ 1]=[x y z 1]·R x (α)

[0046]

[0047] The coordinate transformation when rotating β angle around the Y-axis is as follows:

[0048] [x″ y″ z″ 1]=[x y z 1]·R y (β)

[0049]

[0050] Where (x, y, z) represents the coordinates before rotation, (x′, y′, z′) represents the coordinates after rotating α angle around the X-axis, (x″, y″, z″) represents the coordinates after rotating β angle around the Y-axis, α and β are rotation angles, R x (α) represents the rotation matrix of rotating α angle around the X-axis, and Ry (b) denotes a rotation matrix that rotates by an angle β about the Y axis.

[0051] S5, judging whether the rotated boundary line still has a multi-value function relationship, if yes, adjusting the rotation axis and / or angle according to the residual inclination, repeating steps S4 and S5 until the rotated boundary line has no multi-value function relationship (any position corresponds to a unique elevation value) or reaches a set maximum iteration number, obtaining a final boundary line in a final flattened coordinate system.

[0052] S6, in the final flattened coordinate system, using a Kriging interpolation algorithm to fit the data points of the final boundary line, spatially interpolating the discrete point data of the final boundary line into a continuous elevation surface, and generating an irregular triangular mesh geological interface surface through triangulation, checking the triangular mesh topology on the surface, removing or adjusting abnormal triangular elements such as long and narrow triangles or folded triangular patches. Smooth the elevations of the surface nodes to eliminate sharp changes. After iterative optimization, the surface is smooth and continuous without error topology.

[0053] S7, reversing the optimized geological interface surface and the final boundary line from the final flattened coordinate system back to the original coordinate system. If the boundary line has undergone multiple rotations, apply the inverse angle rotation transformation in the reverse order. For example, if the boundary line has been rotated by an angle β about the Y axis and then by an angle a about the X axis, then restore it by first rotating by an angle -a about the X axis and then by an angle -β about the Y axis. After inverse transformation, the surface and the boundary line can be accurately restored to the correct position in the initial three-dimensional space.

[0054] S8, taking the geological interface surface restored to the original coordinate system as a shear surface, shearing the three-dimensional geological model to segment the three-dimensional solid model of the stratum of different ages. The three-dimensional geological model is the overall three-dimensional geological model of the study area in S1 (e.g., a simple block containing all strata). The initial model is truncated from the corresponding interface using the shear surface to segment the three-dimensional solid model (geological body) of the stratum of that age. The resulting geological body is labeled and removed from the initial model, and the remaining part is retained for subsequent cutting.

[0055] Closedness test: perform a closedness test on the three-dimensional solid model of the stratum of that age obtained by cutting. If some boundaries are found to have open topology (hanging edges) due to cutting, perform hanging edge processing by adding boundary surfaces or extending the surface, etc., to close the open boundaries and ensure that the stratum solid becomes a closed three-dimensional body. After completing the closedness test, the complete three-dimensional model of the stratum of that age is obtained.

[0056] Figure 2For the series of geological profiles of the embodiment of the present application, the complex geological structure of the multiple sets of different age stratum bottom boundary lines is shown, and six sets of different age stratum bottom boundary lines are displayed, each of which is stored in an independent layer and is sequentially marked as 1-6 (from the upper stratum to the lower stratum). According to the data sequence shown, the steps S2-S8 of selecting each age stratum bottom boundary line are repeated layer by layer. The contrast of the unprocessed and rotated stratum bottom boundary lines of the embodiment of the present application is shown in the figure Figure 2 , where the triangular mesh surface structure directly constructed from the unrotated boundary line is disordered and staggered, and the interpolation of the stratum bottom boundary line after the appropriate rotation around the X-axis and Y-axis can generate a regular shear surface. Finally, all the age strata are separated from the initial model, and the fine three-dimensional geological model effect diagram of the complex geological structure is shown in Figure 3 . Figure 4 , which shows the three-dimensional geological model of the geological body obtained by layering according to the geological age sequence by using the method of the present application. Figure 4 Figure 5 is a grid diagram of the internal structure of the three-dimensional geological structure model, which shows the internal structure details of the obtained three-dimensional model and the spatial distribution and mutual relationship of the age stratum interfaces in the model. The method of the present application can effectively meet the three-dimensional modeling needs of complex geological structures and improve the modeling efficiency and accuracy of the results.

[0057] In an exemplary embodiment, a computer readable storage medium is included, and the computer readable storage medium stores a computer program. When the computer program is executed by a processor, the above-mentioned geological interface fitting method based on spatial rotation is realized.

[0058] Please refer to Figure 6 , in an exemplary embodiment, an electronic device is also included, which comprises at least one processor, at least one memory, and at least one communication bus.

[0059] , the memory stores a computer program, and the computer program includes computer readable instructions. The processor calls the computer readable instructions stored in the memory through the communication bus and executes the above-mentioned geological interface fitting method based on spatial rotation.

[0060] In an exemplary embodiment, a computer program product is proposed, which includes a computer program / instruction. When the computer program / instruction is executed by a processor, the steps of the above-mentioned geological interface fitting method based on spatial rotation are realized.

[0061] ​The foregoing description of the disclosed embodiments enables a person skilled in the art to make or use the application. Modifications of these embodiments will occur to persons of skill in the art, and that the appended claims are intended to cover all such modifications that do not depart from the true spirit and scope of the application. Therefore, the application is not limited to the embodiments shown but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A geological interface fitting method based on spatial rotation, characterized in that: The following steps are involved: S1. Obtaining the bottom boundary data of strata of different ages in the three-dimensional geological model from the geological exploration data of the target area; S2. Identifying a stratum bottom boundary having a multi-valued function relationship from the stratum bottom boundary data, wherein the multi-valued function relationship means that a plane coordinate point of a stratum boundary in three-dimensional space corresponds to multiple elevation values; S3, determining the inclination direction and inclination angle of the bottom boundary of the stratum having a multi-valued function relationship in space, selecting a rotation axis according to the inclination direction, and calculating an initial rotation angle value according to the inclination angle; S4, performing coordinate transformation on the bottom boundary of the stratum having a multi-valued function relationship around the rotation axis according to the initial rotation angle value to obtain a rotated boundary; S5. Determine whether the rotated boundary still has a multi-valued function relationship. If so, adjust the rotation axis and / or angle according to the residual tilt. Repeat steps S4 and S5 until the rotated boundary no longer has a multi-valued function relationship or the maximum number of iterations is reached, thereby obtaining the final boundary in the final flattened coordinate system. S6. Fitting the data points of the final boundary in the final flattened coordinate system to generate a geological interface surface, and optimizing the geological interface surface; S7, reversely rotating the optimized geological interface surface and the final boundary from the final flattened coordinate system back to the original coordinate system; S8. Using the geological interface surface restored to the original coordinate system as a shearing surface, shearing the three-dimensional geological model, and segmenting the three-dimensional solid models of strata of different ages.

2. A geological interface fitting method based on spatial rotation according to claim 1, characterized in that: The rotation axis is selected according to the tilt direction: if the boundary line tilts to the north or south, choose to rotate around the X axis; if it tilts to the east or west, choose to rotate around the Y axis; if the tilt direction of the boundary line is not to the north, south, east, or west, choose to rotate around the X axis and Y axis respectively.

3. The geological interface fitting method based on spatial rotation according to claim 1, characterized in that: The inclination angle is the inclination angle of the bottom boundary of the stratum with respect to the horizontal plane having a multi-valued function relationship.

4. The geological interface fitting method based on spatial rotation according to claim 1, characterized in that: The coordinate transformation relationship when rotating around the X axis by an angle α is as follows: [x′ y′ z′ 1]=[xyz 1]·R x (α) The coordinate transformation when rotating around the Y axis by an angle of β is as follows: [x″ y″ z″ 1]=[xyz 1]·R y (β) Where (x, y, z) represents the coordinates before rotation, (x′, y′, z′) represents the coordinates after rotation around the X axis by an angle α, and (x″, y″, z″) represents the coordinates after rotation around the Y axis by an angle β. α and β are the rotation angles, and R x (α) represents the rotation matrix of the angle α around the X axis, R y (β) represents the rotation matrix for the angle β around the Y axis.

5. The geological interface fitting method based on spatial rotation according to claim 1, characterized in that: The Kriging interpolation algorithm is used to fit the data points of the final boundary, and the discrete point data of the final boundary are spatially interpolated into a continuous elevation surface, and the irregular triangulated geological interface surface is generated through triangulation.

6. The geological interface fitting method based on spatial rotation according to claim 1, characterized in that: After segmenting the 3D solid models of strata of different ages, the closure of the 3D solid models of strata of different ages is verified, and the open boundaries caused by cutting are checked and repaired.

7. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the method according to any one of claims 1 to 6 is implemented.

8. An electronic device, characterized in that: The method comprises a processor and a memory, wherein the processor and the memory are interconnected, wherein the memory is used to store a computer program, the computer program includes computer-readable instructions, and the processor is configured to call the computer-readable instructions to execute the method according to any one of claims 1 to 6.

9. A computer program product comprising a computer program / instructions, characterized in that When the computer program / instructions are executed by a processor, the steps of the method according to any one of claims 1 to 6 are implemented.

Images