Investigation data driven three-dimensional space mineral substance migration path generation method

Abstract

The invention relates to an exploration data-driven three-dimensional space mineral migration path generation method, which comprises the following steps of: based on exploration data of a three-dimensional geological space, calculating distance and direction weighting of space points in a neighborhood of each voxel to obtain an optimal direction of each voxel; calculating a pairing point closest to each voxel in the optimal direction, and connecting the pairing points to form a migration path; according to the connection condition, the migration attribute of the space point is judged, and the migration attribute is a starting point, an end point or an intermediate point. The invention innovatively provides a migration path generation method oriented to deep three-dimensional space formation of minerals from a quantitative angle.

Description

Technical Field

[0001] The present invention relates to the field of geological mineral technology, and in particular to a method for generating a three-dimensional material-forming material migration path driven by exploration data. Background Art

[0002] Existing three-dimensional ore-forming material migration path simulation technologies mainly rely on numerical simulation and geomechanical models, which simulate the migration and enrichment paths of ore-forming materials in geological structures by coupling fluid dynamics and rock physical parameters. However, these technologies have the following problems and shortcomings: (1) The models are difficult to accurately simulate the interaction between complex strata, structures and ore-forming fluids; (2) The description of the non-stationarity of local geological strata and structures (such as the dynamic changes in the physical and chemical properties of ore-forming fluids) is not accurate enough, which affects the simulation accuracy; (3) The input parameters are heavily dependent on local geological exploration data, resulting in poor versatility of the methods. In addition, the high computational complexity and low computational efficiency of the models limit their application in large-scale mining area prediction. Summary of the Invention

[0003] In order to solve the problems existing in the above-mentioned prior art, the purpose of the present invention is to propose a method for generating migration paths of ore-forming materials in three-dimensional space driven by exploration data, taking the anisotropy of ore-controlling effects into consideration, and innovatively proposing a method for generating migration paths of ore-forming materials in deep three-dimensional space from a quantitative perspective.

[0004] To achieve the above object, the present invention provides the following solutions:

[0005] A method for generating 3D spatial ore-forming material migration paths driven by exploration data, comprising:

[0006] Based on the exploration data of three-dimensional geological space, the distance and direction of the spatial points in the neighborhood of each voxel are weighted and calculated to obtain the optimal direction of each voxel;

[0007] Calculate the paired points closest to each voxel in the optimal direction and connect them to form a migration path; determine the migration attribute of the spatial point based on the connection situation, wherein the migration attribute is a starting point, an end point, or an intermediate point.

[0008] Optionally, obtaining the optimal direction includes:

[0009] Based on the direction relationship model p1, a weighted distance and direction operation is performed on the spatial point:

[0010]

[0011] Among them, w ij is the weight, d ij is the distance from the voxel to the spatial point j, d maxis a parameter determined by the bandwidth, is the distance weight, α ij is the direction of the voxel, p i are the directions in the direction relationship model p1, For direction rights;

[0012] By comparing the performance of the spatial point in each direction in the direction relationship model p1, that is, the weight w ij , the weight w ij The direction closest to 1 is taken as the optimal direction.

[0013] Optionally, connecting the paired points includes:

[0014] On the basis of obtaining the optimal direction of each voxel, the paired point closest to the spatial point is obtained, and any spatial point is selected. If the spatial point has the next paired point, it is connected, otherwise it is not connected, until all spatial points are selected.

[0015] Optionally, judging the migration attribute of the spatial point according to the connection condition includes:

[0016] On the basis of the directional connection in the connection situation, the migration attribute of each element is determined according to the connection situation of each point.

[0017] Optionally, when the spatial point has only a corresponding post-matching point, the point attribute is determined to be a starting point.

[0018] Optionally, when the spatial point has a corresponding front matching point and a rear matching point, it is determined that the point attribute is an intermediate point.

[0019] Optionally, when the spatial point has only a corresponding preceding matching point, the point attribute is determined to be an end point.

[0020] The beneficial effects of the present invention are:

[0021] Based on geological exploration data, this paper takes into account the influence of three-dimensional spatial ore-controlling factors on spatial non-stationarity such as mineralization anisotropy, and introduces directional weights in GWR to characterize the directional characteristics of non-stationary relationships. From a quantitative perspective, it innovatively proposes a migration path generation method for deep three-dimensional ore-forming materials. BRIEF DESCRIPTION OF THE DRAWINGS

[0022] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0023] Figure 1 This is a flow chart of optimal direction calculation according to an embodiment of the present invention;

[0024] Figure 2 This is a flow chart of a directional connection algorithm according to an embodiment of the present invention;

[0025] Figure 3 A visualization diagram of regional ore-forming material migration paths according to an embodiment of the present invention;

[0026] Figure 4 This is a schematic diagram of a method for generating three-dimensional ore-forming material migration paths driven by exploration data in an embodiment of the present invention. DETAILED DESCRIPTION

[0027] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of the present invention.

[0028] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the present invention is further described in detail below with reference to the accompanying drawings and specific embodiments.

[0029] like Figure 4 As shown, this embodiment discloses a method for generating a three-dimensional space mineralization migration path driven by exploration data, including: based on the exploration data of the three-dimensional geological space, performing weighted operations on the distance and direction of the spatial points in the neighborhood of each element to obtain the optimal direction of each element; calculating the paired points closest to each element in the optimal direction, and connecting them to form a migration path; judging the migration attributes of the spatial points according to the connection situation, wherein the migration attributes are the starting point, end point or intermediate point.

[0030] Specifically, the technical solution of this embodiment is divided into three steps. First, the optimal direction of the spatial position point is calculated based on the exploration data, then the direction is connected based on the direction relationship model, and finally the point attributes are judged and classified to obtain the three-dimensional spatial mineralization migration path.

[0031] Furthermore, obtaining the optimal direction includes: performing a weighted operation of distance and direction on the spatial points based on the direction relationship model p1:

[0032]

[0033] Among them, w ij is the weight, d ij is the distance from the voxel to the spatial point j, d max is a parameter determined by the bandwidth, is the distance weight, α ij is the direction of the voxel, p i are the directions in the direction relationship model p1, For direction rights;

[0034] By comparing the performance of spatial points in each direction in the direction relationship model p1, that is, the weight w ij , the weight w ij The direction closest to 1 is taken as the optimal direction.

[0035] Specifically, the optimal direction is calculated:

[0036] Based on the exploration data of three-dimensional geological space, that is, the point data arranged in a cubic manner within the neighborhood of each voxel, the array CDirectedVertex(x c ,y c ,z c ,num,q j ,q x ) to represent the location of the ore-forming fluid in the three-dimensional geological space, x c 、y c 、z c Respectively represent the coordinates of the point on the x, y, and z axes, num represents the number of the point, and q j Indicates the angle between the direction of the point and the z-axis, q x Indicates the angle between the projection of the point on the xy plane and the north direction.

[0037] according to Figure 1 The process shown calculates the optimal direction of a spatial point. When the input spatial point p c After the spatial points in its neighborhood, the distance and direction weighted operation is performed based on the direction relationship model p1, and then the performance of each direction in the direction relationship model p1 is compared, that is, the weight w ij , take weight w ij The direction closest to 1 is the optimal direction.

[0038] Furthermore, connecting the paired points includes: obtaining the paired point closest to the spatial point based on the optimal direction of each voxel, selecting any spatial point, and connecting it if the spatial point has the next paired point, otherwise not connecting it, until all spatial points are selected.

[0039] Specifically, all the connectable spatial points and the spatial points in their neighborhood are connected. Figure 2 The directional connection algorithm shown in the figure starts by selecting any spatial point. If the spatial point has a next matching point, it is connected. Otherwise, it is not connected until all spatial points have been selected. Once all connectable points are connected, the migration path of the mineralization material in three-dimensional space is obtained.

[0040] Furthermore, judging the migration attribute of the spatial point according to the connection condition includes: judging the migration attribute of each voxel according to the connection condition of each point on the basis of the directional connection in the connection condition.

[0041] Specifically, based on whether there are front and back matching points, the point attributes are judged and classified, as shown in Table 1. When a spatial point only has a corresponding back matching point, the point attribute is judged to be a starting point. When a spatial point has corresponding front and back matching points, the point attribute is judged to be an intermediate point. When a spatial point only has a corresponding front matching point, the point attribute is judged to be an end point.

[0042] Table 1 Judgment of point migration attributes based on the front and back paired points

[0043]

[0044]

[0045] like Figure 3 As shown in the figure, the point-line image produced based on the 26-directional relationship model can better show the directional relationship between points in the local area, more clearly show the directional status of the area with complex point directions, and improve the accuracy of the simulation generation of the migration path of mineralizing materials in deep three-dimensional space.

[0046] Taking the Xiadian gold deposit in Jiaodong as an example, the optimal direction of the spatial position point is first calculated based on the exploration data of the deposit in the directional relationship model; then the paired points are calculated and the directions are connected; finally, the point attributes are judged and classified, and the migration path of the mineralizing materials in three-dimensional space is visualized.

[0047] Calculate the optimal direction: For any point in three-dimensional space, press Figure 1 The process shown in the figure calculates the optimal direction of the point. When the spatial point p c After the environmental data is calculated, a distance weighted calculation is performed, and a direction weighted calculation is performed in the direction p1 (direction relationship model). Finally, the performance of each direction in the direction relationship model is compared to obtain the optimal direction and its weighted residual square sum R 2 and the regression coefficient β i .

[0048] Directional connection: Calculate the pairing points of the cubically arranged spatial points according to their optimal direction, and connect them two by two through the directional connection algorithm to generate line files. After gradual connection, the migration direction and trend of the spatial position points are formed to express the migration path of mineralization materials in three-dimensional space.

[0049] The connection algorithm flow is as follows:

[0050] Step 1: Input point p c and its neighborhood;

[0051] Step 2: Determine whether there is a matching point. If so, connect it. Otherwise, repeat the first step and input the next point pi.

[0052] Step 3: Determine whether all points have been entered. If so, end. Otherwise, repeat the first step to enter the next point p. c .

[0053] Visual analysis: Classify point attributes into three types based on whether they have previous and subsequent paired points.

[0054] After the direction connection is completed, the generated line table is imported into the visualization software, and the image can reflect the overall directional trend of the mineralization migration. This method is used to show the dynamic effect of the ore-controlling factors on the mineralization of the Xiadian gold deposit. The direction connection and visualization of the point data of the 10*10*10 cubic distribution are performed, and the output results are as follows: Figure 3 As shown in the figure, it shows that for data points with large direction differences arranged in a cubic pattern, the direction connection method can more intuitively reflect the overall migration path of regional mineralization.

[0055] The embodiments described above are merely descriptions of preferred embodiments of the present invention and are not intended to limit the scope of the present invention. Without departing from the spirit of the present invention, various modifications and improvements made to the technical solutions of the present invention by persons skilled in the art should fall within the scope of protection defined by the claims of the present invention.

Claims

1. A method for generating a three-dimensional spatial ore-forming material migration path driven by exploration data, characterized in that: include: Based on the exploration data of three-dimensional geological space, the distance and direction of the spatial points in the neighborhood of each voxel are weighted and calculated to obtain the optimal direction of each voxel; Calculate the paired points closest to each voxel in the optimal direction and connect them to form a migration path; determine the migration attribute of the spatial point based on the connection situation, wherein the migration attribute is a starting point, an end point, or an intermediate point.

2. The method for generating a three-dimensional ore-forming material migration path driven by exploration data according to claim 1, characterized in that: Obtaining the optimal direction includes: Based on the direction relationship model p1, a weighted distance and direction operation is performed on the spatial point: Among them, w ij is the weight, d ij is the distance from the voxel to the spatial point j, d max is a parameter determined by the bandwidth, is the distance weight, α ij is the direction of the voxel, p i are the directions in the direction relationship model p1, For direction rights; By comparing the performance of the spatial point in each direction in the direction relationship model p1, that is, the weight w ij , the weight w ij The direction closest to 1 is taken as the optimal direction.

3. The method for generating a three-dimensional ore-forming material migration path driven by exploration data according to claim 1, characterized in that: Connecting paired points involves: On the basis of obtaining the optimal direction of each voxel, the paired point closest to the spatial point is obtained, and any spatial point is selected. If the spatial point has the next paired point, it is connected, otherwise it is not connected, until all spatial points are selected.

4. The method for generating a three-dimensional ore-forming material migration path driven by exploration data according to claim 1, characterized in that: Determining the migration attribute of the spatial point according to the connection condition includes: On the basis of the directional connection in the connection situation, the migration attribute of each element is determined according to the connection situation of each point.

5. The method for generating a three-dimensional ore-forming material migration path driven by exploration data according to claim 3, characterized in that: When the spatial point has only a corresponding post-matching point, the point attribute is determined to be a starting point.

6. The method for generating a three-dimensional ore-forming material migration path driven by exploration data according to claim 3, characterized in that: When the spatial point has a corresponding front matching point and a rear matching point, it is determined that the point attribute is an intermediate point.

7. The method for generating a three-dimensional ore-forming material migration path driven by exploration data according to claim 3, characterized in that: When the spatial point has only a corresponding preceding matching point, the point attribute is determined to be an end point.

Images