Multi-level feature analysis method for three-dimensional geological structural surface and related equipment

Abstract

The invention provides a multilevel feature analysis method of a three-dimensional geologic structural plane and related equipment, and the method comprises the steps: carrying out the modeling of a geologic body structural plane according to the coordinates of a known point on the constructed geologic body structural plane, so as to construct a geologic trend plane, and extracting a remaining trend plane from the geologic body structural plane; determining a radial basis function in the initial implicit function of the geological trend surface and solving a linear coefficient item to obtain an implicit function of the geological trend surface including smooth item constraint; adjusting a weight constrained by a smooth item in the implicit function of the geological trend surface to obtain a multistage trend surface feature function of the geological structural surface, and substituting the multistage trend surface feature function into position coordinates of all points on the geological body structural surface to obtain a feature implicit function value corresponding to each point; and based on the feature implicit function values corresponding to all the points, a multi-level feature analysis result of the geologic body structural plane is obtained, and the problem of accurate fitting of the three-dimensional geologic structural plane under the condition that the known points are few is solved.

Description

Technical Field

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

[0002] With the continuous development of computer graphics technology and the deepening of theoretical research on three-dimensional spatial data processing, information technology conditions have been provided for the establishment of three-dimensional geological structural surfaces. By establishing three-dimensional geological structural surfaces, not only can the morphological characteristics of geological structures and the spatial relationship of structural elements be expressed vividly and intuitively, but the analysis of the geometric morphology and spatial distribution characteristics of three-dimensional geological structural surfaces is of great significance for the prediction of geological disasters, optimization of engineering design, simulation of geological processes, prediction of mineral resources, oil and gas exploration and other applications.

[0003] The currently widely used method for modeling three-dimensional geological structural surfaces is to use the traditional "outcrop point, outcrop line + structural surface attitude" to construct a three-dimensional geological structural surface. This modeling method fuzzies or performs geological deduction on the basic source data of the geological structural surface. The established three-dimensional geological structural surface cannot fully and truly reflect the development status of the geological structural surface. For geological three-dimensional structural surfaces with complex morphology, traditional methods may find it difficult to use implicit functions to fit their characteristic details when there are few known points, making it difficult to simulate the actual geological interface well, resulting in a low resolution of the geological structural model and difficulty in observing its overall geometric morphology and spatial distribution characteristics. Summary of the invention

[0004] The present invention provides a multi-level feature analysis method for a three-dimensional geological structure surface and related equipment, the purpose of which is to solve the problem of accurate fitting of the three-dimensional geological structure surface when there are fewer known points.

[0005] In order to achieve the above object, the present invention provides a multi-level feature analysis method for a three-dimensional geological structure surface, comprising: Step 1, modeling the geological body structural surface according to the coordinates of the known points on the constructed geological body structural surface, and obtaining the implicit function of the geological body structural surface; Step 2, constructing a geological trend surface based on the implicit function of the geological body structural surface, and extracting the structural surface with undulating characteristics from the geological body structural surface to obtain a residual trend surface, wherein the implicit function of the geological body structural surface is equal to the sum of the implicit function of the geological trend surface and the implicit function of the residual trend surface; Step 3, determining that the radial basis function in the initial implicit function of the geological trend surface is a surface radial basis function with Green's function as the basis function, and solving 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 surface and the surface radial basis function, to obtain the implicit function of the geological trend surface, wherein the implicit function of the geological trend surface includes a weighted smooth term constraint; Step 4, adjusting the weight of the smoothing term constraint in the implicit function of the geological trend surface to obtain a multi-level trend surface characteristic function of the geological structure surface; Step 5, substitute the position coordinates of all points on the geological body structural surface into the characteristic function of each level trend surface, obtain the characteristic implicit function value corresponding to each point, and obtain the multi-level characteristic analysis results of the geological body structural surface based on the characteristic implicit function values ​​corresponding to all points.

[0006] More specifically, step 1 includes: Extract geological boundaries from the acquired drilling data, profile data and geological map data of the study area, and perform three-dimensional modeling of the geological body structural surface in the form of point coordinates to obtain the geological body structural surface of the study area; The geological body structural surface is modeled according to the known point coordinates on the geological body structural surface and the Hermite radial basis function, and the implicit function of the geological body structural surface is obtained as follows: in, The implicit function representing the structural surface of the geological body, Represents the number of known points on the structural surface of the geological body, is the linear coefficient term, , , represents the implicit function value of the initial spatial constraint point, Represents the matrix form of radial basis function, is a radial basis function that satisfies , The first A known point, Represents any known point on the structural surface of a geological body.

[0007] More specifically, step 2 includes: The geological trend surface is constructed based on the implicit function of the geological body structural surface. The initial implicit function of the geological trend surface is: in, represents the linear coefficient term of the initial implicit function, , The radial basis function representing the initial implicit function; Extract the structural surface with undulating characteristics from the structural surface of the geological body to obtain the residual trend surface; The implicit function of the geological structure surface is equal to the sum of the implicit function of the geological trend surface and the implicit function of the residual trend surface, and the expression is: in, Implicit function representing the residual trend surface.

[0008] More specifically, step 3 includes: According to the variational method, the radial basis function in the initial implicit function of the geological trend surface is determined as the surface radial basis function with Green's function as the basis function: According to the known point coordinates on the geological body structure surface and the radial basis function of the curved 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: Based on the surface radial basis function and the analytical function of the linear coefficient term, the implicit function of the geological trend surface is obtained as follows: in, represents the Green's function, represents the matrix form of Green's function, represents the weight of the smoothing term, represents a unit vector.

[0009] More specifically, step 4 includes: 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 characteristic function of the geological structure surface, including: in, The level that represents the characteristics of the trend surface.

[0010] Furthermore, based on the characteristic implicit function values ​​corresponding to all points, the multi-level characteristic analysis results of the geological structure surface are obtained, including: Define the implicit function of the geological structure surface to be equal to 0, and the relationship between the implicit function of the geological trend surface and the residual trend surface is: According to the relationship between the implicit function of the geological trend surface and the residual trend surface, the characteristic implicit function values ​​corresponding to all points are converted into the implicit function values ​​of the residual trend surface. The implicit function values ​​of the residual trend surface are used to characterize the undulating characteristics of the geological body structural surface, and the implicit function values ​​of the residual trend surface are used as the multi-level characteristic analysis results of the geological body structural surface.

[0011] The present invention also provides a multi-level feature analysis device for a three-dimensional geological structure surface, comprising: A modeling module is used to model the geological body structural surface according to the coordinates of known points on the constructed geological body structural surface to obtain the implicit function of the geological body structural surface; A construction module is used to construct a geological trend surface based on the implicit function of the geological body structural surface, and extract the structural surface with undulating characteristics from the geological body structural surface to obtain the residual trend surface, wherein the implicit function of the geological body structural surface 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 is used to determine that the radial basis function in the initial implicit function of the geological trend surface is a surface radial basis function with 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 structure surface and the surface radial basis function to obtain the implicit function of the geological trend surface, wherein the implicit function of the geological trend surface includes a weighted smooth term constraint; An adjustment module is used 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 characteristic function of the geological structure surface; The substitution module is used to substitute the position coordinates of all points on the geological body structural surface into the characteristic function of each level of trend surface, obtain the characteristic implicit function value corresponding to each point, and obtain the multi-level characteristic analysis results of the geological body structural surface based on the characteristic implicit function values ​​corresponding to all points.

[0012] 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 a three-dimensional geological structure surface is implemented.

[0013] The present invention also provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, a multi-level feature analysis method for a three-dimensional geological structure surface is implemented.

[0014] The above scheme of the present invention has the following beneficial effects: The invention models the geological body structural surface according to the coordinates of known points on the constructed geological body structural surface to obtain the implicit function of the geological body structural surface; constructs a geological trend surface based on the implicit function of the geological body structural surface, extracts the structural surface with undulating characteristics in the geological body structural surface to obtain the residual trend surface, and the implicit function of the geological body structural surface is equal to the sum of the implicit function of the geological trend surface and the implicit function of the residual trend surface; determines that the radial basis function in the initial implicit function of the geological trend surface is a curved surface radial basis function with 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 coordinates of known points on the geological body structural surface and the curved surface radial basis function to obtain the implicit function of the geological trend surface including the smooth term constraint; and calculate ... The weight of the smoothing term constraint in the implicit function of the potential 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 body structural surface are substituted into the characteristic function of each level of the trend surface to obtain the characteristic implicit function value corresponding to each point, and the multi-level characteristic analysis result of the geological body structural surface is obtained based on the characteristic implicit function values ​​corresponding to all points; compared with the prior art, the present invention realizes multi-level analysis of the characteristics of the geological structure surface by constructing the trend surface of the geological body structural surface, fundamentally solves the problem of accurate fitting of the three-dimensional geological structure surface when there are fewer known points, and the trend surface still satisfies the known point position constraints and the smoothing term constraints, and the obtained multi-level characteristic trend surface can observe the characteristics of the geological structure surface from multiple dimensions.

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

[0016] Figure 1 It is a schematic diagram of a flow chart of an embodiment of the present invention; Figure 2 Schematic diagram of the structure of a multi-stage feature analysis device in an embodiment of the present invention; Figure 3 Schematic diagram of the structure of a terminal device in an embodiment of the present invention. DETAILED DESCRIPTION

[0017] In order 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 part of the embodiments of the present invention, rather than all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0018] In the description of the present invention, it should be noted that the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", etc., indicating the orientation or positional relationship, are based on the orientation or positional relationship shown in the drawings, and are 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 therefore cannot be understood as limiting the present invention. In addition, the terms "first", "second", and "third" are used for descriptive purposes only, and cannot be understood as indicating or implying relative importance.

[0019] In the description of the present invention, it should be noted that, unless otherwise clearly specified and limited, the terms "installation", "connection" and "connection" 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 a direct connection or an indirect connection through an intermediate medium, or it can be the internal communication of two components. For ordinary technicians in this field, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

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

[0021] In view of the existing problems, the present invention provides a multi-level characteristic analysis method of a three-dimensional geological structure surface and related equipment.

[0022] like Figure 1 As shown, an embodiment of the present invention provides a multi-level feature analysis method for a three-dimensional geological structure surface, comprising: Step 1, modeling the geological body structural surface according to the coordinates of the known points on the constructed geological body structural surface, and obtaining the implicit function of the geological body structural surface; Step 2, constructing a geological trend surface based on the implicit function of the geological body structural surface, and extracting the structural surface with undulating characteristics from the geological body structural surface to obtain a residual trend surface, wherein the implicit function of the geological body structural surface is equal to the sum of the implicit function of the geological trend surface and the implicit function of the residual trend surface; Step 3, determining that the radial basis function in the initial implicit function of the geological trend surface is a surface radial basis function with Green's function as the basis function, and solving 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 surface and the surface radial basis function, to obtain the implicit function of the geological trend surface, wherein the implicit function of the geological trend surface includes a weighted smooth term constraint; Step 4, adjusting the weight of the smoothing term constraint in the implicit function of the geological trend surface to obtain a multi-level trend surface characteristic function of the geological structure surface; Step 5, substitute the position coordinates of all points on the geological body structural surface into the characteristic function of each level trend surface, obtain the characteristic implicit function value corresponding to each point, and obtain the multi-level characteristic analysis results of the geological body structural surface based on the characteristic implicit function values ​​corresponding to all points.

[0023] Specifically, step 1 includes: Extract geological boundaries from the acquired drilling data, profile data and geological map data of the study area, and perform three-dimensional modeling of the geological body structural surface in the form of point coordinates to obtain the geological body structural surface of the study area; In three-dimensional geological space, the field function is defined as , to represent any point in three-dimensional space Scalar function at , the geological body structural surface is modeled according to the known point coordinates on the geological body structural surface and the Hermite radial basis function, and the implicit function of the geological body structural surface is obtained as follows: in, The implicit function representing the structural surface of the geological body, Represents the number of known points on the structural surface of the geological body, is the linear coefficient term, , , represents the implicit function value of the initial spatial constraint point, Represents the matrix form of radial basis function, is a radial basis function that satisfies , The first A known point, Represents any known point on the structural surface of a geological body; The above implicit function satisfies: Based on the above relationship, the linear coefficient term in the implicit function can be solved , thereby obtaining a continuous function expression of the initial structural surface model in three-dimensional geological space.

[0024] It should be noted that borehole data, profile data and geological map data are the most commonly used data in 3D geological modeling.

[0025] Since the geological structure surface with very complex morphology and distribution characteristics requires a large number of constraint points to establish implicit functions, thus forming a relatively complex linear system, on this basis, it is very difficult to analyze the characteristics of the geological structure surface, so step 2 includes: The geological trend surface is constructed based on the implicit function of the geological body structural surface. The initial implicit function of the geological trend surface is: in, represents the linear coefficient term of the initial implicit function, , The radial basis function representing the initial implicit function; Extract the structural surface with undulating characteristics from the structural surface of the geological body to obtain the residual trend surface; The implicit function of the geological structure surface is equal to the sum of the implicit function of the geological trend surface and the implicit function of the residual trend surface, and the expression is: in, Implicit function representing the residual trend surface.

[0026] The above formula can transform the complex geological body structural surface characteristic analysis problem into the residual trend surface problem based on the implicit function of the trend surface. The geological body structural surface characteristics can be represented by the value of the residual trend surface. The residual trend surface value can be directly represented as the characteristic implicit function value of the point on the corresponding geological body structural surface on the geological trend surface, and it is negatively correlated with the implicit function value on the trend surface.

[0027] Specifically, step 3 includes: According to the variational method, the radial basis function in the initial implicit function of the geological trend surface is determined as the surface radial basis function with Green's function as the basis function: In order to obtain the implicit function of the geological trend surface of the above-mentioned geological body structural surface, its linear coefficient needs to be solved. Although the geological trend surface is smoother than the geological structure surface in space, it still needs to meet the spatial position constraints of the known points, that is, at the position of the known points in space, the characteristic implicit function value of the trend surface still needs to be equal to 0. Therefore, according to the coordinates of the known points on the geological body structural 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: Based on the surface radial basis function and the analytical function of the linear coefficient term, the implicit function of the geological trend surface is obtained as follows: in, represents the Green's function, represents the matrix form of Green's function, , represents the weight of the smoothing term, Represents a unit vector, avoids the underdetermined solution process, and realizes the smoothness of the characteristic trend surface.

[0028] In order to solve the geological trend surface of the geological body structure, the embodiment of the present invention adds a weighted smooth term constraint to the implicit function of the geological body structure surface to ensure that the geological trend surface has spatial smoothness under the premise of fitting the known points in space. Therefore, the implicit function solution of the geological trend surface can be transformed into a variational problem, and the process is as follows: First, define the objective function of the variational problem as: in, is the objective function of the variational problem, are the coordinates of a known point in space, Characterize the initial implicit function value of the geological trend surface in the geological body structure surface, is the weight of the smoothing term to control the degree of smoothing. When the weight of the smoothing term is lower, the generated geological trend surface is closer to the original structural surface. is the Laplace smoothing term, which can be characterized as the Hessian operator, and its general form is: By solving the second-order partial derivatives, the distribution of geological three-dimensional structural surfaces in space can be made smoother and more continuous; The embodiment of the present invention uses the Gauss-Newton algorithm to convert the nonlinear constraints in the objective function into Transformed into a first-order linear constraint term, that is, the objective function is transformed into solving a low-order linear problem: in, It can be regarded as the step length term for iterative update of the characteristic implicit function of the trend surface.

[0029] The process of converting nonlinear constraints into first-order linear constraints using the Gauss-Newton algorithm in the embodiment of the present invention is as follows: The above Expanding the Gauss-Newton algorithm, we can obtain:

[0030] Here, due to the small iteration step size, the higher-order infinitesimal terms can be discarded , so the further variational form of the above problem can be characterized as: For the above variational expression based on the update step smoothing term, a regularized function solving method with a smoothing term can be used to obtain a surface radial basis function with the Green's function as the basis function.

[0031] Specifically, step 4 includes: 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 characteristic function of the geological structure surface, including: in, The series representing the trend surface characteristics has different linear coefficients for different levels of trend surface implicit functions, and the weights have different orders of magnitude at different levels.

[0032] Specifically, step 5 includes: Any point on the geological structure surface Substituting into the characteristic function of each level trend surface, we can get the characteristic implicit function value corresponding to each point, and we can get: Based on the characteristic implicit function values ​​corresponding to all points, the multi-level characteristic analysis results of the geological structure surface are obtained, including: Define the implicit function of the geological structure surface to be equal to 0, and the relationship between the implicit function of the geological trend surface and the residual trend surface is: According to the relationship between the implicit function of the geological trend surface and the residual trend surface, the characteristic implicit function values ​​corresponding to all points are converted into the implicit function values ​​of the residual trend surface. The implicit function value of the residual trend surface is used to characterize the undulating characteristics of the geological body structural surface. The expression is: Finally, the implicit function value of the remaining trend surface is used as the multi-level characteristic analysis result of the geological body structural surface. The multi-level characteristic analysis result includes the undulating characteristics, inclination change characteristics, and surface concave-convex characteristics of the geological body structural surface.

[0033] The embodiment of the present invention models the geological body structural surface according to the known point coordinates on the constructed geological body structural surface to obtain the implicit function of the geological body structural surface; constructs a geological trend surface based on the implicit function of the geological body structural surface, extracts the structural surface with undulating characteristics in the geological body structural surface to obtain the residual trend surface, and the implicit function of the geological body structural surface is equal to the sum of the implicit function of the geological trend surface and the implicit function of the residual trend surface; determines that the radial basis function in the initial implicit function of the geological trend surface is a curved surface radial basis function with 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 surface and the curved surface radial basis function to obtain the implicit function of the geological trend surface including the smooth term constraint; The weight of the smoothing term constraint in the implicit function of the potential 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 body structural surface are substituted into the characteristic function of each level of the trend surface to obtain the characteristic implicit function value corresponding to each point, and the multi-level characteristic analysis result of the geological body structural surface is obtained based on the characteristic implicit function values ​​corresponding to all points; compared with the prior art, the embodiment of the present invention realizes multi-level analysis of the characteristics of the geological structure surface by constructing the trend surface of the geological body structural surface, fundamentally solves the problem of accurate fitting of the three-dimensional geological structure surface when there are fewer known points, and the trend surface still satisfies the known point position constraints and the smoothing term constraints, and the obtained multi-level characteristic trend surface can observe the characteristics of the geological structure surface from multiple dimensions.

[0034] Corresponding to the multi-level feature analysis method of the three-dimensional geological structure surface described in the above embodiment, Figure 2 As shown, the present invention also provides a multi-level feature analysis device 100 for a three-dimensional geological structure surface, the multi-level feature analysis device 100 comprising: A modeling module 101 is used to model the geological body structural surface according to the coordinates of known points on the constructed geological body structural surface to obtain an implicit function of the geological body structural surface; A construction module 102 is used to construct a geological trend surface based on the implicit function of the geological body structural surface, and extract the structural surface with undulating characteristics from the geological body structural surface to obtain a residual trend surface, wherein the implicit function of the geological body structural surface 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 103 is used to determine that the radial basis function in the initial implicit function of the geological trend surface is a surface radial basis function with 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 structure surface and the surface radial basis function to obtain the implicit function of the geological trend surface, wherein the implicit function of the geological trend surface includes a weighted smooth term constraint; An adjustment module 104 is used 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 characteristic function of the geological structure surface; The substitution module 105 is used to substitute the position coordinates of all points on the geological body structural surface into each level 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 surface based on the characteristic implicit function values ​​corresponding to all points.

[0035] It should be noted that the information interaction, execution process, etc. between the above-mentioned devices / units are based on the same concept as the method embodiment of the present application. Their specific functions and technical effects can be found in the method embodiment part and will not be repeated here.

[0036] The technicians in the relevant field can clearly understand that for the convenience and simplicity of description, only the division of the above-mentioned functional units and modules is used as an example for illustration. In practical applications, the above-mentioned function allocation can be completed by different functional units and modules as needed, that is, the internal structure of the device can be divided into different functional units or modules to complete all or part of the functions described above. The functional units and modules in the embodiment can be integrated in a processing unit, or each unit can exist physically separately, or two or more units can be integrated in one unit. The above-mentioned integrated unit can be implemented in the form of hardware or in the form of software functional units. In addition, the specific names of the functional units and modules are only for the convenience of distinguishing each other, and are not used to limit the scope of protection of this application. The specific working process of the units and modules in the above-mentioned system can refer to the corresponding process in the aforementioned method embodiment, which will not be repeated here.

[0037] The embodiment of the present invention also provides a terminal device, such as Figure 3 As shown, the terminal device D10 of 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 above-mentioned multi-level feature analysis method of the three-dimensional geological structure surface is implemented.

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

[0039] The processor D100 may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or any conventional processor, etc.

[0040] In some embodiments, the memory D101 may be an internal storage unit of the terminal device D10, such as a 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 memory card (SMC, SmartMedia Card), a secure digital (SD, Secure Digital) card, a flash card (Flash Card), etc. equipped on the terminal device D10. Further, the memory D101 may also include both an internal storage unit of the terminal device D10 and an external storage device. The memory D101 is used to store an operating system, an application program, a boot loader (BootLoader), 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.

[0041] It should be noted that the information interaction, execution process, etc. between the above-mentioned devices / units are based on the same concept as the method embodiment of the present application. Their specific functions and technical effects can be found in the method embodiment part and will not be repeated here.

[0042] The technicians in the relevant field can clearly understand that for the convenience and simplicity of description, only the division of the above-mentioned functional units and modules is used as an example for illustration. In practical applications, the above-mentioned function allocation can be completed by different functional units and modules as needed, that is, the internal structure of the device can be divided into different functional units or modules to complete all or part of the functions described above. The functional units and modules in the embodiment can be integrated in a processing unit, or each unit can exist physically separately, or two or more units can be integrated in one unit. The above-mentioned integrated unit can be implemented in the form of hardware or in the form of software functional units. In addition, the specific names of the functional units and modules are only for the convenience of distinguishing each other, and are not used to limit the scope of protection of this application. The specific working process of the units and modules in the above-mentioned system can refer to the corresponding process in the aforementioned method embodiment, which will not be repeated here.

[0043] The present invention also provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, a multi-level feature analysis method for a three-dimensional geological structure surface is implemented.

[0044] 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, the present application implements all or part of the processes in the above-mentioned embodiment method, which can be completed by instructing the relevant hardware through a computer program, and the computer program can be stored in a computer-readable storage medium, and the computer program can implement the steps of the above-mentioned various method embodiments when executed by the processor. Among them, the computer program includes computer program code, and the computer program code can be in source code form, object code form, executable file or some intermediate form. The computer-readable medium may at least include: any entity or device that can carry 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), electric carrier signal, telecommunication signal and software distribution medium. For example, a USB flash drive, a mobile hard disk, a disk or an optical disk.

[0045] The above is a preferred embodiment of the present invention. It should be pointed out that for ordinary technicians in this technical field, several improvements and modifications can be made without departing from the principles of the present invention. These improvements and modifications should also be regarded as the scope of protection of the present invention.

Claims

1. A multi-level feature analysis method for a three-dimensional geological structure surface, characterized in that: include: Step 1, modeling the geological body structural surface according to the coordinates of known points on the constructed geological body structural surface to obtain the implicit function of the geological body structural surface; Step 2, constructing a geological trend surface based on the implicit function of the geological body structural surface, and extracting a structural surface with undulating characteristics from the geological body structural surface to obtain a residual trend surface, wherein the implicit function of the geological body structural surface is equal to the sum of the implicit function of the geological trend surface and the implicit function of the residual trend surface; Step 3, determining that the radial basis function in the initial implicit function of the geological trend surface is a surface radial basis function with Green's function as the basis function, and solving 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 surface and the surface radial basis function to obtain the implicit function of the geological trend surface, wherein the implicit function of the geological trend surface includes a weighted smooth term constraint; Step 4, adjusting the weight of the smoothing term constraint in the implicit function of the geological trend surface to obtain a multi-level trend surface characteristic function of the geological structure surface; Step 5, substitute the position coordinates of all points on the geological body structural surface into each level 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 surface based on the characteristic implicit function values ​​corresponding to all points.

2. The multi-level feature analysis method of a three-dimensional geological structure surface according to claim 1, characterized in that: The step 1 comprises: Extract geological boundaries from the acquired drilling data, profile data and geological map data of the study area, and perform three-dimensional modeling of the geological body structural surface in the form of point coordinates to obtain the geological body structural surface of the study area; The geological body structural surface is modeled according to the known point coordinates on the geological body structural surface and the Hermite radial basis function, and the implicit function of the geological body structural surface is obtained as follows: in, The implicit function representing the structural surface of the geological body, Represents the number of known points on the structural surface of the geological body, is the linear coefficient term of the geological body structural surface, , , represents the implicit function value of the initial spatial constraint point, Represents the matrix form of radial basis function, is a radial basis function that satisfies , The first A known point, Represents any known point on the structural surface of a geological body.

3. The multi-level feature analysis method of a three-dimensional geological structure surface according to claim 2, characterized in that: The step 2 comprises: A geological trend surface is constructed based on the implicit function of the geological body structural surface. The initial implicit function of the geological trend surface is: in, represents the linear coefficient term of the initial implicit function, , The radial basis function representing the initial implicit function; Extracting structural surfaces with undulating characteristics from the structural surfaces of the geological body to obtain residual trend surfaces; The implicit function of the geological body structural surface is equal to the sum of the implicit function of the geological trend surface and the implicit function of the residual trend surface, and the expression is: in, Implicit function representing the residual trend surface.

4. The multi-level feature analysis method of a three-dimensional geological structure surface according to claim 3, characterized in that: The step 3 comprises: According to the variational method principle, the radial basis function in the initial implicit function of the geological trend surface is determined to be a surface radial basis function with Green's function as the basis function: According to the known point coordinates on the geological body structural surface and the surface radial basis function, 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: Based on the surface radial basis function and the analytical function of the linear coefficient term, the implicit function of the geological trend surface is obtained as follows: in, represents the Green's function, represents the matrix form of Green's function, represents the weight of the smoothing term, represents a unit vector.

5. The multi-level feature analysis method of a three-dimensional geological structure surface according to claim 4, characterized in that: The step 4 comprises: 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 characteristic function of the geological structure surface, including: in, The level that represents the characteristics of the trend surface.

6. The multi-level feature analysis method of a three-dimensional geological structure surface according to claim 5, characterized in that: Based on the characteristic implicit function values ​​corresponding to all points, the multi-level characteristic analysis results of the geological body structural surface are obtained, including: The implicit function of the geological body structural surface is defined as 0, and the relationship between the implicit function of the geological trend surface and the residual trend surface is obtained as follows: According to the relationship between the implicit function of the geological trend surface and the residual trend surface, the characteristic implicit function values ​​corresponding to all points are converted into implicit function values ​​of the residual trend surface. The implicit function values ​​of the residual trend surface are used to characterize the undulating characteristics of the geological body structural surface. The implicit function values ​​of the residual trend surface are used as the multi-level characteristic analysis results of the geological body structural surface.

7. A multi-level feature analysis device for a three-dimensional geological structure surface, characterized in that: include: A modeling module, used for modeling the geological body structural surface according to the coordinates of known points on the constructed geological body structural surface, and obtaining the implicit function of the geological body structural surface; A construction module, used for constructing a geological trend surface based on the implicit function of the geological body structural surface, and extracting a structural surface with undulating characteristics from the geological body structural surface to obtain a residual trend surface, wherein the implicit function of the geological body structural surface 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, used for determining that the radial basis function in the initial implicit function of the geological trend surface is a surface radial basis function with Green's function as the basis function, and solving 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 surface and the surface radial basis function to obtain the implicit function of the geological trend surface, wherein the implicit function of the geological trend surface includes a weighted smooth term constraint; An adjustment module, used for adjusting the weight of the smoothing term constraint in the implicit function of the geological trend surface to obtain a multi-level trend surface characteristic function of the geological structure surface; The substitution module is used to substitute the position coordinates of all points on the geological body structural surface into the characteristic function of each level trend surface to obtain the characteristic implicit function value corresponding to each point, and obtain the multi-level characteristic analysis results of the geological body structural surface 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, the multi-level feature analysis method for a three-dimensional geological structure surface as described in any one of claims 1 to 6 is implemented.

9. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the multi-level feature analysis method for a three-dimensional geological structure surface as described in any one of claims 1 to 6 is implemented.

Images