A multi-level feature analysis method for three-dimensional geological structural planes and related equipment

Abstract

The present invention provides a multi-level feature analysis method and related equipment for three-dimensional geological structural planes. According to the known point coordinates on the constructed geological body structural plane, the geological body structural plane is modeled to construct a geological trend surface, and the residual trend surface is extracted from the geological body structural plane; the radial basis function in the initial implicit function of the geological trend surface is determined and the linear coefficient term is solved to obtain the implicit function of the geological trend surface including the smoothing term constraint; the weight of the smoothing term constraint in the implicit function of the geological trend surface is adjusted to obtain the multi-level trend surface feature function of the geological structural plane and substitute the position coordinates of all points on the geological body structural plane to obtain the characteristic implicit function value corresponding to each point, and the multi-level feature analysis result of the geological body structural plane is obtained based on the characteristic implicit function values corresponding to all points, solving the problem of accurate fitting of the three-dimensional geological structural plane in the case of fewer known points.

Description

Technical Field

[0001] The present invention relates to the technical field of three-dimensional geological modeling, and particularly relates to a method for multi-level feature analysis of three-dimensional geological structural planes and related equipment. Background Art

[0002] With the continuous development of computer graphics technology and the continuous in-depth research on three-dimensional space data processing theory, information technology conditions are provided for establishing three-dimensional geological structural planes. By establishing three-dimensional geological structural planes, not only can the morphological characteristics of geological structures and the spatial relationships of structural elements be visually expressed, but also the analysis of the geometric morphology and spatial distribution characteristics of three-dimensional geological structural planes is of great significance for applications such as predicting geological disasters, optimizing engineering designs, simulating geological processes, predicting mineral resources, and oil and gas exploration.

[0003] For the modeling of three-dimensional geological structural planes, the currently widely used method is to construct three-dimensional geological structural planes by using traditional "outcrop points, outcrop lines + structural plane attitudes". This modeling method performs fuzzy processing or geological deduction processing on the basic source data of geological structural planes, and the established three-dimensional geological structural planes cannot comprehensively and truly reflect the development state of geological structural planes; for complex three-dimensional geological structural planes, in the case of fewer known points, it may be difficult to use implicit functions to fit their characteristic details by traditional methods, thus it is difficult to better simulate the actual geological interface, resulting in a low resolution of the geological structure model and making it difficult to observe its overall geometric morphology and spatial distribution characteristics. Summary of the Invention

[0004] The present invention provides a method for multi-level feature analysis of three-dimensional geological structural planes and related equipment, and its purpose is to solve the problem of accurate fitting of three-dimensional geological structural planes in the case of fewer known points.

[0005] To achieve the above purpose, the present invention provides a method for multi-level feature analysis of three-dimensional geological structural planes, including:

[0006] Step 1, based on the known point coordinates on the constructed geological body structural plane, model the geological body structural plane to obtain the implicit function of the geological body structural plane;

[0007] Step 2, construct a geological trend surface based on the implicit function of the geological body structural plane, and extract the structural plane with undulating characteristics in the geological body structural plane to obtain the remaining trend surface. The implicit function of the geological body structural plane is equal to the sum of the implicit functions of the geological trend surface and the remaining trend surface;

[0008] Step 3: Determine that the radial basis function in the initial implicit function of the geological trend surface is the surface radial basis function with the Green's function as the basis function. Solve the linear coefficient term of the initial implicit function of the geological trend surface according to the known point coordinates on the geological body structural plane and the surface radial basis function to obtain the implicit function of the geological trend surface. The implicit function of the geological trend surface includes a weighted smoothing term constraint;

[0009] Step 4: Adjust the weight of the smoothing term constraint in the implicit function of the geological trend surface to obtain the multi-level trend surface characteristic function of the geological structural plane;

[0010] Step 5: Substitute the position coordinates of all points on the geological body structural plane into each level of the trend surface characteristic function to obtain the characteristic implicit function value corresponding to each point, and obtain the multi-level characteristic analysis result of the geological body structural plane based on the characteristic implicit function values corresponding to all points.

[0011] Furthermore, Step 1 includes:

[0012] Extract the geological boundary from the obtained borehole data, profile data, and geological map data of the research area, and perform three-dimensional modeling of the geological body structural plane in the form of point coordinates to obtain the geological body structural plane of the research area;

[0013] Model the geological body structural plane according to the known point coordinates on the geological body structural plane and the Hermite radial basis function, and the implicit function of the geological body structural plane is:

[0014]

[0015] Wherein, represents the implicit function of the geological body structural plane, represents the number of known points on the geological body structural plane, is the linear coefficient term, , , represents the implicit function value of the initial spatial constraint point, represents the matrix form of the radial basis function, is the radial basis function, which satisfies , is the th known point on the geological body structural plane, represents any known point on the geological body structural plane.

[0016] Furthermore, Step 2 includes:

[0017] Construct a geological trend surface based on the implicit function of the geological body structural plane. The initial implicit function of the geological trend surface is:

[0018]

[0019] Among them, represents the linear coefficient term of the initial implicit function, , represents the radial basis function of the initial implicit function;

[0020] Extract the structural plane with undulating characteristics from the geological body structural plane to obtain the remaining trend surface;

[0021] The implicit function of the geological body structural plane is equal to the sum of the implicit function of the geological trend surface and the implicit function of the remaining trend surface, and the expression is:

[0022]

[0023] Among them, represents the implicit function of the remaining trend surface.

[0024] Furthermore, step 3 includes:

[0025] According to the principle of variational method, determine that the radial basis function in the initial implicit function of the geological trend surface is the surface radial basis function with the Green's function as the basis function:

[0026]

[0027] According to the known point coordinates on the geological body structural plane and the surface radial basis function, solve the linear coefficient term of the initial implicit function of the geological trend surface, and the analytical function of the linear coefficient term is:

[0028]

[0029] Based on the surface radial basis function and the analytical function of the linear coefficient term, obtain the implicit function of the geological trend surface as:

[0030]

[0031] Among them, represents the Green's function, represents the matrix form of the Green's function, represents the smoothing term weight, represents the unit vector.

[0032] Furthermore, step 4 includes:

[0033] Adjust the weight of the smoothing term constraint in the implicit function of the geological trend surface to obtain the multi-level trend surface characteristic function of the geological structural plane, including:

[0034]

[0035]

[0036]

[0037]

[0038] Among them, represents the order of the trend surface features.

[0039] Furthermore, based on the characteristic implicit function values corresponding to all points, the multi-level characteristic analysis results of the geological body structural plane are obtained, including:

[0040] Define the implicit function of the geological body structural plane to be equal to 0, and the relationship between the implicit function of the geological trend surface and the remaining trend surface is obtained as:

[0041]

[0042] According to the relationship between the implicit function of the geological trend surface and the remaining trend surface, the characteristic implicit function values corresponding to all points are converted into the implicit function values of the remaining trend surface. The implicit function values of the remaining trend surface are used to characterize the undulating characteristics of the geological body structural plane, and the implicit function values of the remaining trend surface are used as the multi-level characteristic analysis results of the geological body structural plane.

[0043] The present invention also provides a multi-level characteristic analysis device for a three-dimensional geological structural plane, including:

[0044] A modeling module, configured to model the geological body structural plane according to the known point coordinates on the constructed geological body structural plane, and obtain the implicit function of the geological body structural plane;

[0045] A construction module, configured to construct a geological trend surface based on the implicit function of the geological body structural plane, and extract the structural plane with undulating characteristics in the geological body structural plane to obtain the remaining trend surface. The implicit function of the geological body structural plane is equal to the sum of the implicit function of the geological trend surface and the implicit function of the remaining trend surface;

[0046] A solving module, configured to determine that the radial basis function in the initial implicit function of the geological trend surface is a surface radial basis function with the Green's function as the basis function, and solve the linear coefficient term of the initial implicit function of the geological trend surface according to the known point coordinates on the geological body structural plane and the surface radial basis function, so as to obtain the implicit function of the geological trend surface. The implicit function of the geological trend surface includes a weighted smoothing term constraint;

[0047] An adjustment module, configured to adjust the weight of the smoothing term constraint in the implicit function of the geological trend surface to obtain the multi-level trend surface characteristic function of the geological structural plane;

[0048] A substitution module, configured to substitute the position coordinates of all points on the geological body structural plane into each level of trend surface characteristic function to obtain the characteristic implicit function value corresponding to each point, and obtain the multi-level characteristic analysis results of the geological body structural plane based on the characteristic implicit function values corresponding to all points.

[0049] The present invention also provides a terminal device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, a multi-level feature analysis method for three-dimensional geological structural planes is implemented.

[0050] The present invention also provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements a multi-level feature analysis method for three-dimensional geological structural planes.

[0051] The above solution of the present invention has the following beneficial effects:

[0052] Based on the known point coordinates on the constructed geological body structural plane, the present invention models the geological body structural plane to obtain the implicit function of the geological body structural plane; constructs a geological trend surface based on the implicit function of the geological body structural plane, and extracts the structural plane with undulating characteristics in the geological body structural plane to obtain the remaining trend surface. The implicit function of the geological body structural plane is equal to the sum of the implicit functions of the geological trend surface and the remaining trend surface; determines that the radial basis function in the initial implicit function of the geological trend surface is a surface radial basis function with the Green's function as the basis function, and solves the linear coefficient term of the initial implicit function of the geological trend surface according to the known point coordinates on the geological body structural plane and the surface radial basis function to obtain the implicit function of the geological trend surface including the smooth term constraint; adjusts the weight of the smooth term constraint in the implicit function of the geological trend surface to obtain the multi-level trend surface characteristic function of the geological structural plane; substitutes the position coordinates of all points on the geological body structural plane into each level of the trend surface characteristic function to obtain the characteristic implicit function value corresponding to each point, and obtains the multi-level feature analysis result of the geological body structural plane based on the characteristic implicit function values corresponding to all points; compared with the prior art, the present invention constructs a trend surface of the geological body structural plane to achieve multi-level analysis of the characteristics of the geological structural plane, fundamentally solves the problem of accurate fitting of the three-dimensional geological structural plane in the case of fewer known points, and the trend surface still satisfies the known point position constraint and the smooth term constraint, and the obtained multi-level characteristic trend surface can observe the characteristics of the geological structural plane from multiple dimensions.

[0053] Other beneficial effects of the present invention will be described in detail in the subsequent specific implementation part. BRIEF DESCRIPTION OF THE DRAWINGS

[0054] Figure 1 It is a schematic flowchart of an embodiment of the present invention;

[0055] Figure 2 It is a schematic structural diagram of a multi-level feature analysis device in an embodiment of the present invention;

[0056] Figure 3 It is a schematic structural diagram of a terminal device in an embodiment of the present invention. DETAILED DESCRIPTION OF THE INVENTION

[0057] To make the technical problems, technical solutions and advantages to be solved by the present invention clearer, the following will be described in detail with reference to the accompanying drawings and specific embodiments. Obviously, the described embodiments are some embodiments of the present invention, rather than all embodiments. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0058] In the description of the present invention, it should be noted that the orientation or positional relationship indicated by the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", etc. is based on the orientation or positional relationship shown in the drawings. It is only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and thus cannot be construed as a limitation of the present invention. In addition, the terms "first", "second", "third" are only used for descriptive purposes and cannot be construed as indicating or implying relative importance.

[0059] In the description of the present invention, it should be noted that unless otherwise clearly specified and defined, the terms "installed", "connected", "connected" should be understood in a broad sense. For example, it can be a locking connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be directly connected or indirectly connected through an intermediate medium, and it can be the communication inside two elements. For those of ordinary skill in the art, the specific meanings of the above terms in the present invention can be understood according to specific situations.

[0060] In addition, the technical features involved in different embodiments of the present invention described below can be combined with each other as long as they do not conflict with each other.

[0061] The present invention provides a multi-level feature analysis method and related equipment for three-dimensional geological structural planes in view of existing problems.

[0062] As Figure 1 shown, an embodiment of the present invention provides a multi-level feature analysis method for three-dimensional geological structural planes, including:

[0063] Step 1, according to the known point coordinates on the constructed geological body structural plane, model the geological body structural plane to obtain the implicit function of the geological body structural plane;

[0064] Step 2, construct a geological trend surface based on the implicit function of the geological body structural plane, and extract the structural plane with undulating characteristics in the geological body structural plane to obtain the remaining trend surface. The implicit function of the geological body structural plane is equal to the sum of the implicit function of the geological trend surface and the implicit function of the remaining trend surface;

[0065] Step 3: Determine that the radial basis function in the initial implicit function of the geological trend surface is the surface radial basis function with the Green's function as the basis function. Solve the linear coefficient term of the initial implicit function of the geological trend surface according to the known point coordinates on the geological body structural plane and the surface radial basis function to obtain the implicit function of the geological trend surface, and the implicit function of the geological trend surface includes the weighted smoothing term constraint;

[0066] Step 4: Adjust the weight of the smoothing term constraint in the implicit function of the geological trend surface to obtain the multi-level trend surface characteristic function of the geological structural plane;

[0067] Step 5: Substitute the position coordinates of all points on the geological body structural plane into each level of the trend surface characteristic function to obtain the characteristic implicit function value corresponding to each point, and obtain the multi-level characteristic analysis result of the geological body structural plane based on the characteristic implicit function values corresponding to all points.

[0068] Specifically, Step 1 includes:

[0069] Extract the geological boundary from the obtained borehole data, profile data, and geological map data of the study area, and perform three-dimensional modeling of the geological body structural plane in the form of point coordinates to obtain the geological body structural plane of the study area;

[0070] In the three-dimensional geological space, define the field function as , to represent the scalar function at any point in the three-dimensional space. Model the geological body structural plane according to the known point coordinates on the geological body structural plane and the Hermite radial basis function, and obtain the implicit function of the geological body structural plane as:

[0071]

[0072] Among them, represents the implicit function of the geological body structural plane, represents the number of known points on the geological body structural plane, is the linear coefficient term, , , represents the implicit function value of the initial space constraint point, represents the matrix form of the radial basis function, is the radial basis function, which satisfies , is the th known point on the geological body structural plane, represents any one of the known points on the geological body structural plane;

[0073] The above implicit function satisfies:

[0074]

[0075] Based on the above relationship, the linear coefficient term in the implicit function can be solved. , and thus a continuous function expression of the initial structural plane model in the three-dimensional geological space is obtained.

[0076] It should be noted that borehole data, profile data, and geological map data are the most commonly used data in three-dimensional geological modeling.

[0077] Since a large number of constraint points are required to establish an implicit function for the structural plane of a geological body with very complex morphology and distribution characteristics, a relatively complex linear system is formed. On this basis, it is difficult to analyze the characteristics of the structural plane of this geological body. Therefore, step 2 includes:

[0078] Construct a geological trend surface based on the implicit function of the structural plane of the geological body. The initial implicit function of the geological trend surface is:

[0079]

[0080] Among them, represents the linear coefficient term of the initial implicit function, , represents the radial basis function of the initial implicit function;

[0081] Extract the structural plane with undulating characteristics in the structural plane of the geological body to obtain the residual trend surface;

[0082] The implicit function of the structural plane of the geological body is equal to the sum of the implicit function of the geological trend surface and the implicit function of the residual trend surface. The expression is:

[0083]

[0084] Among them, represents the implicit function of the residual trend surface.

[0085] From the above formula, the problem of analyzing the characteristics of the complex structural plane of the geological body can be transformed into the problem of solving the residual trend surface based on the implicit function of the trend surface. The characteristics of the structural plane of the geological body can be represented by the value of the residual trend surface. The numerical value of the residual trend surface can be directly characterized as the characteristic implicit function value of the point on the corresponding structural plane of the geological body on the geological trend surface, and shows a negative correlation with the implicit function value on the trend surface.

[0086] Specifically, step 3 includes:

[0087] According to the principle of variational method, it is determined that the radial basis function in the initial implicit function of the geological trend surface is the surface radial basis function with the Green's function as the basis function:

[0088]

[0089] To obtain the implicit function of the geological trend surface of the geological discontinuity surface mentioned above, it is also necessary to solve its linear coefficients. Although the geological trend surface is smoother than the geological discontinuity surface in space, it still needs to conform to the spatial position constraints of the known points, that is, at the positions of the known points in space, the implicit function value of the trend surface feature still needs to be equal to 0. Therefore, according to the coordinates of the known points on the geological discontinuity surface and the radial basis function of the surface, the linear coefficient term of the initial implicit function of the geological trend surface is solved, and the analytical function of the linear coefficient term is obtained as follows:

[0090]

[0091] Based on the radial basis function of the surface and the analytical function of the linear coefficient term, the implicit function of the geological trend surface is obtained as follows:

[0092]

[0093] Among them, represents the Green's function, represents the matrix form of the Green's function, , represents the smoothing term weight, represents the unit vector, avoiding the underdetermined solution process and realizing the smoothing characteristics of the characteristic trend surface.

[0094] In the embodiment of the present invention, in order to solve the geological trend surface of the geological structure, a weighted smoothing term constraint is added to the implicit function of the geological discontinuity surface to ensure that the geological trend surface has spatial smoothness on the premise of fitting the known points in space. Thus, the solution of the implicit function of the geological trend surface can be transformed into a variational problem, and the process is as follows:

[0095] First, the objective function of the variational problem is defined as follows:

[0096]

[0097] Among them, is the objective function of the variational problem, are the coordinates of the known points in space, represents the initial implicit function value of the geological trend surface in the geological discontinuity surface, is the smoothing term weight to control the degree of smoothing. When the smoothing term weight is lower, the generated geological trend surface is closer to the original discontinuity surface, is the Laplacian smoothing term, which can be represented as the Hessian operator, and its general form is as follows:

[0098]

[0099]

[0100] By solving the second-order partial derivatives, the distribution of the geological three-dimensional structural plane in space can be made smoother and continuous.

[0101] In the embodiment of the present invention, the Gauss-Newton algorithm is used to transform the non-linear constraints in the objective function into first-order linear constraint terms, that is, the objective function is transformed into solving a low-order linear problem:

[0102]

[0103] wherein, can be regarded as the step size term for the iterative update of the implicit function of the trend surface feature.

[0104] The process of transforming the non-linear constraints into first-order linear constraint terms by using the Gauss-Newton algorithm in the embodiment of the present invention is as follows:

[0105] Expand the above using the Gauss-Newton algorithm, and the following can be obtained:

[0106]

[0107] Here, since the iterative step size is small, the high-order infinitesimal terms can be discarded , and the further variational form of the above problem can be characterized as:

[0108]

[0109] For the above variational expression based on the updated step size smoothing term, the solution method of the regularization function with a smoothing term can be used to obtain the surface radial basis function with the Green's function as the basis function.

[0110] Specifically, step 4 includes:

[0111] Adjust the weight of the smoothing term constraint in the implicit function of the geological trend surface to obtain the multi-level trend surface characteristic function of the geological structural surface, including:

[0112]

[0113]

[0114]

[0115]

[0116] wherein, represents the series number of the trend surface feature. For the implicit functions of different-level trend surfaces, there are different linear coefficient terms, and the weights have different orders of magnitude at different levels.

[0117] Specifically, step 5 includes:

[0118] Substitute any point on the geological body structural plane into the characteristic function of each level of the trend surface to obtain the characteristic implicit function value corresponding to each point, and we can get:

[0119]

[0120]

[0121]

[0122]

[0123] Based on the characteristic implicit function values corresponding to all points, obtain the multi-level characteristic analysis results of the geological body structural plane, including:

[0124] Define the implicit function of the geological body structural plane to be equal to 0, and the relationship between the implicit function of the geological trend surface and the remaining trend surface is:

[0125]

[0126] According to the relationship between the implicit function of the geological trend surface and the remaining trend surface, convert the characteristic implicit function values corresponding to all points into the implicit function values of the remaining trend surface. The implicit function values of the remaining trend surface are used to characterize the undulation characteristics of the geological body structural plane, and the expression is:

[0127]

[0128]

[0129]

[0130]

[0131] Finally, take the implicit function values of the remaining trend surface as the multi-level characteristic analysis results of the geological body structural plane. The multi-level characteristic analysis results include the undulation characteristics, dip angle change characteristics, and convexity and concavity characteristics of the surface of the geological body structural plane.

[0132] In an embodiment of the present invention, based on the coordinates of known points on a geological structure surface, a model of the geological structure surface is constructed to obtain an implicit function of the geological structure surface; a geological trend surface is constructed based on the implicit function of the geological structure surface, and a structure surface with undulating characteristics is extracted from the geological structure surface to obtain a residual trend surface, where the implicit function of the geological structure surface is equal to the sum of the implicit functions of the geological trend surface and the residual trend surface; it is determined that the radial basis function in the initial implicit function of the geological trend surface is a surface radial basis function with a Green's function as the basis function, and based on the coordinates of known points on the geological structure surface and the surface radial basis function, the linear coefficient term of the initial implicit function of the geological trend surface is solved to obtain an implicit function of the geological trend surface including a smoothing term constraint; the weight of the smoothing term constraint in the implicit function of the geological trend surface is adjusted to obtain a multi-level trend surface characteristic function of the geological structure surface; the position coordinates of all points on the geological structure surface are substituted into each level of the trend surface characteristic function to obtain a characteristic implicit function value corresponding to each point, and a multi-level characteristic analysis result of the geological structure surface is obtained based on the characteristic implicit function values corresponding to all points; compared with the prior art, in the embodiment of the present invention, by constructing a trend surface of the geological structure surface, multi-level analysis of the characteristics of the geological structure surface is realized, fundamentally solving the problem of accurate fitting of a three-dimensional geological structure surface in the case of fewer known points, and the trend surface still satisfies the known point position constraint and the smoothing term constraint, and the obtained multi-level characteristic trend surface can observe the characteristics of the geological structure surface from multiple dimensions.

[0133] Corresponding to the multi-level characteristic analysis method of the three-dimensional geological structure surface described in the above embodiment, as Figure 2 shown, the present invention also provides a multi-level characteristic analysis device 100 for a three-dimensional geological structure surface. The multi-level characteristic analysis device 100 includes:

[0134] A modeling module 101, configured to construct a model of the geological structure surface based on the coordinates of known points on the constructed geological structure surface to obtain an implicit function of the geological structure surface;

[0135] A construction module 102, configured to construct a geological trend surface based on the implicit function of the geological structure surface, and extract a structure surface with undulating characteristics from the geological structure surface to obtain a residual trend surface, where the implicit function of the geological structure surface is equal to the sum of the implicit functions of the geological trend surface and the residual trend surface;

[0136] A solving module 103, configured to determine that the radial basis function in the initial implicit function of the geological trend surface is a surface radial basis function with a Green's function as the basis function, and based on the coordinates of known points on the geological structure surface and the surface radial basis function, solve the linear coefficient term of the initial implicit function of the geological trend surface to obtain an implicit function of the geological trend surface, where the implicit function of the geological trend surface includes a weighted smoothing term constraint;

[0137] An adjustment module 104 is configured to adjust the weight of the smoothing term constraint in the implicit function of the geological trend surface to obtain a multi-level trend surface feature function of the geological structural surface;

[0138] A substitution module 105 is configured to substitute the position coordinates of all points on the geological body structure surface into each level of the trend surface feature function to obtain the corresponding feature implicit function value for each point, and obtain a multi-level feature analysis result of the geological body structure surface based on the corresponding feature implicit function values of all points.

[0139] It should be noted that for the information interaction, execution process, etc. between the above-mentioned devices / units, since they are based on the same concept as the method embodiment of the present application, for their specific functions and the technical effects brought, please refer to the method embodiment part for details, and will not be elaborated here.

[0140] Those skilled in the art can clearly understand that for the convenience and simplicity of description, only the above-mentioned division of each functional unit and module is used as an example. In actual applications, the above-mentioned functions can be allocated to different functional units and modules according to needs, that is, the internal structure of the device is divided into different functional units or modules to complete all or part of the functions described above. Each functional unit and module in the embodiment can be integrated into a processing unit, or each unit can exist physically alone, or two or more units can be integrated into one unit. The above-mentioned integrated unit can be implemented in the form of hardware or in the form of a software functional unit. In addition, the specific names of each functional unit and module are only for the convenience of mutual distinction and do not limit the protection scope of the present application. The specific working process of the units and modules in the above system can refer to the corresponding process in the foregoing method embodiment and will not be elaborated here.

[0141] An embodiment of the present invention also provides a terminal device, such as Figure 3 shown. The terminal device D10 in this embodiment includes: at least one processor D100 ( Figure 3 only one processor is shown in the figure), a memory D101, and a computer program D102 stored in the memory D101 and executable on the at least one processor D100. When the processor D100 executes the computer program D102, the multi-level feature analysis method of the three-dimensional geological structural surface described above is implemented.

[0142] The terminal device D10 can be a computing device such as a desktop computer, a notebook, a palm computer, a server, a server cluster, and a cloud server. The terminal device may include, but is not limited to, a processor D100 and a memory D101. Those skilled in the art can understand, Figure 3This is only an example of the terminal device D10, which does not constitute a limitation on the terminal device D10. It may include more or fewer components than those shown in the figure, or combine certain components, or different components. For example, it may also include input / output devices, network access devices, etc.

[0143] The so-called processor D100 may be a central processing unit (CPU), and this processor D100 may also be other general-purpose processors, digital signal processors (DSPs), application specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor may be a microprocessor or any conventional processor, etc.

[0144] In some embodiments, the memory D101 may be an internal storage unit of the terminal device D10, such as the hard disk or memory of the terminal device D10. In other embodiments, the memory D101 may also be an external storage device of the terminal device D10, such as a plug-in hard disk, a smart media card (SMC), a secure digital (SD) card, a flash card, etc., equipped on the terminal device D10. Further, the memory D101 may also include both the internal storage unit and the external storage device of the terminal device D10. The memory D101 is used to store an operating system, application programs, a boot loader, data, and other programs, such as the program code of the computer program, etc. The memory D101 may also be used to temporarily store data that has been output or is to be output.

[0145] It should be noted that for the content such as information interaction and execution process between the above-mentioned devices / units, since it is based on the same concept as the method embodiment of this application, for its specific functions and the technical effects brought, reference may be specifically made to the method embodiment part, and details are not described herein again.

[0146] Those skilled in the art can clearly understand that, for the convenience and conciseness of description, only the division of the above-mentioned functional units and modules is used as an example. In actual applications, the above-mentioned functions can be allocated to different functional units and modules according to needs, that is, the internal structure of the device is divided into different functional units or modules to complete all or part of the functions described above. Each functional unit and module in the embodiments can be integrated into a processing unit, or each unit can exist physically alone, or two or more units can be integrated into one unit. The above-mentioned integrated unit can be implemented in the form of hardware or in the form of a software functional unit. In addition, the specific names of the functional units and modules are only for the convenience of mutual distinction and do not limit the protection scope of this application. The specific working processes of the units and modules in the above system can refer to the corresponding processes in the foregoing method embodiments and will not be elaborated here.

[0147] The present invention also provides a computer-readable storage medium storing a computer program, which when executed by a processor implements the method for multi-level feature analysis of three-dimensional geological structural planes.

[0148] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, to implement all or part of the processes in the above method embodiments of the present application, a computer program can be used to instruct the relevant hardware to complete. The computer program can be stored in a computer-readable storage medium. When the computer program is executed by a processor, the steps of the above method embodiments can be implemented. Among them, the computer program includes computer program code, and the computer program code can be in the form of source code, object code, executable file or some intermediate form. The computer-readable medium can at least include: any entity or device capable of carrying the computer program code to the construction device / terminal device, recording medium, computer memory, read-only memory (ROM, Read-Only Memory), random access memory (RAM, RandomAccess Memory), electrical carrier signal, telecommunication signal, and software distribution medium. For example, a USB flash drive, a mobile hard disk, a magnetic disk or an optical disc, etc.

[0149] The above is the preferred embodiment of the present invention. It should be noted that for those of ordinary skill in the art, without departing from the principle of the present invention, several improvements and refinements can be made, and these improvements and refinements should also be regarded as the protection scope of the present invention.

Claims

1. A method for multi-level feature analysis of three-dimensional geological structural planes, characterized in that, Including: Step 1: Model the geological structure plane according to the known point coordinates on the constructed geological structure plane to obtain the implicit function of the geological structure plane. Step 2: Construct a geological trend surface based on the implicit function of the geological structure plane, and extract the structure plane with undulating characteristics in the geological structure plane to obtain the remaining trend surface. The implicit function of the geological structure plane is equal to the sum of the implicit function of the geological trend surface and the implicit function of the remaining trend surface. Step 3: Determine that the radial basis function in the initial implicit function of the geological trend surface is a surface radial basis function with the Green's function as the basis function. According to the known point coordinates on the geological structure plane and the surface radial basis function, solve the linear coefficient term of the initial implicit function of the geological trend surface to obtain the implicit function of the geological trend surface, and the implicit function of the geological trend surface includes a weighted smoothing term constraint. Step 4: Adjust the weight of the smoothing term constraint in the implicit function of the geological trend surface to obtain the multi-level trend surface characteristic function of the geological structure plane. Step 5: Substitute the position coordinates of all points on the geological structure plane into each level of the trend surface characteristic function to obtain the characteristic implicit function value corresponding to each point, and obtain the multi-level characteristic analysis result of the geological structure plane based on the characteristic implicit function values corresponding to all points.

2. The multi-level feature analysis method of the three-dimensional geological structural plane according to claim 1, characterized in that The said Step 1 includes: Extract the geological boundary from the obtained borehole data, profile data and geological map data of the research area, and perform three-dimensional modeling of the geological structure plane in the form of point coordinates to obtain the geological structure plane of the research area. Model the geological structure plane according to the known point coordinates on the geological structure plane and the Hermite radial basis function, and the obtained implicit function of the geological structure plane is: Among them, represents the implicit function of the geological body structural plane, represents the number of known points on the geological body structural plane, is the linear coefficient term of the geological body structural plane, , , represents the implicit function value of the initial space constraint point, represents the matrix form of the radial basis function, is the radial basis function, which satisfies , is the th known point on the geological body structural plane, represents any known point on the geological body structural plane.

3. The multi-level feature analysis method of the three-dimensional geological structural plane according to claim 2, wherein The said Step 2 includes: Construct a geological trend surface based on the implicit function of the geological structure plane, and the initial implicit function of the geological trend surface is: Among them, represents the linear coefficient term of the initial implicit function, , represents the radial basis function of the initial implicit function; Extract the structure plane with undulating characteristics in the geological structure plane to obtain the remaining trend surface. The implicit function of the geological structure plane is equal to the sum of the implicit function of the geological trend surface and the implicit function of the remaining trend surface, and the expression is: Among them, represents the implicit function of the remaining trend surface.

4. The multi-level feature analysis method of the three-dimensional geological structural plane according to claim 3, characterized in that The said Step 3 includes: According to the principle of variational method, determine that the radial basis function in the initial implicit function of the geological trend surface is a surface radial basis function with the Green's function as the basis function as: According to the known point coordinates on the geological structure plane and the surface radial basis function, solve the linear coefficient term of the initial implicit function of the geological trend surface to obtain the analytical function of the linear coefficient term as: Based on the surface radial basis function and the analytical function of the linear coefficient term, obtain the implicit function of the geological trend surface as: Among them, represents the Green's function, represents the matrix form of the Green's function, represents the smooth term weight, represents the unit vector.

5. The multi-level feature analysis method of the three-dimensional geological structural plane according to claim 4, characterized in that The said Step 4 includes: Adjust the weight of the smoothing term constraint in the implicit function of the geological trend surface to obtain the multi-level trend surface characteristic function of the geological structure plane including: Among them, represents the order of the trend surface feature.

6. The multi-level feature analysis method of the three-dimensional geological structural plane according to claim 5, characterized in that, Based on the characteristic implicit function values corresponding to all points to obtain the multi-level characteristic analysis result of the geological structure plane, including: Define the implicit function of the geological structure plane to be equal to 0, and obtain the relationship between the implicit function of the geological trend surface and the remaining trend surface as: According to the relationship between the implicit function of the geological trend surface and the residual trend surface, convert the characteristic implicit function values corresponding to all points into the implicit function values of the residual trend surface. The implicit function values of the residual trend surface are used to characterize the undulation characteristics of the structural plane of the geological body, and use the implicit function values of the residual trend surface as the multi-level characteristic analysis results of the structural plane of the geological body.

7. A multi-level feature analysis device for three-dimensional geological structural planes, characterized in that, Including: A modeling module, configured to model the structural plane of the geological body according to the known point coordinates on the constructed structural plane of the geological body, and obtain the implicit function of the structural plane of the geological body; A construction module, configured to construct a geological trend surface based on the implicit function of the structural plane of the geological body, and extract the structural plane with undulation characteristics from the structural plane of the geological body to obtain a residual trend surface. The implicit function of the structural plane of the geological body is equal to the sum of the implicit function of the geological trend surface and the implicit function of the residual trend surface; A solution module, configured to determine that the radial basis function in the initial implicit function of the geological trend surface is a surface radial basis function with a Green's function as the basis function, and solve the linear coefficient term of the initial implicit function of the geological trend surface according to the known point coordinates on the structural plane of the geological body and the surface radial basis function, to obtain the implicit function of the geological trend surface. The implicit function of the geological trend surface includes a weighted smoothing term constraint; An adjustment module, configured to adjust the weight of the smoothing term constraint in the implicit function of the geological trend surface to obtain the multi-level trend surface characteristic function of the geological structural plane; A substitution module, configured to substitute the position coordinates of all points on the structural plane of the geological body into each level of the trend surface characteristic function to obtain the characteristic implicit function value corresponding to each point, and obtain the multi-level characteristic analysis result of the structural plane of the geological body based on the characteristic implicit function values corresponding to all points.

8. A terminal device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the multi-level characteristic analysis method of the three-dimensional geological structural plane according to any one of claims 1 to 6.

9. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by the processor, it implements the multi-level characteristic analysis method of the three-dimensional geological structural plane according to any one of claims 1 to 6.

Images