Method for dividing ore deposit two-dimensional plane finite element by taking exploration project as center

Abstract

The invention provides a method for dividing a mineral deposit two-dimensional plane finite element by taking an exploration project as a center. The method comprises the following steps: acquiring exploration project data; selecting a plane of a to-be-drawn two-dimensional map from the target area as a section plane, and determining an intersection point with the section plane based on the position of the exploration engineering observation / sampling point as a range center point; vertically extending the height of the section by taking the plane of the section as a reference to form a space which is called as the section; giving geological information collected by observation / sampling points in the section range to each path center point; forming a unit plane by using a trap curve / broken line by taking each path center point as a center, and constructing a path center area by taking the unit plane as a reference; and determining a two-dimensional geological distribution diagram of the ore deposit formed by the boundary of each range center area. According to the method, the mineralization system with high inhomogeneity is divided into mutually independent units on the plane, so that the organic defect of a geostatistical method can be avoided, and the actual distribution condition of ore deposit resources can be reflected more accurately.

Description

Technical Field

[0001] The present invention relates to the technical field of mineral deposit resource exploration, and in particular to a method for dividing a mineral deposit two-dimensional plane finite element with exploration engineering as the center. Background Art

[0002] Mineral deposits, especially endogenous deposits, are often controlled by multiple ore-controlling factors, with complex occurrences and often containing multiple superpositions of ore-forming elements. How to identify and objectively reflect the appearance of mineralization, so as to guide exploration and development, is an important topic that has long been of concern to academia and industry.

[0003] Exploration projects represented by drilling can provide a variety of information including lithology, occurrence, structure, structure, mineralization, alteration, and crypto-explosion and fragmentation. Especially with the advancement and popularization of testing and analysis instruments, a large amount of elemental analysis data has been accumulated, allowing mineral exploration to progress from the qualitative stage to the quantitative and semi-quantitative stage.

[0004] Although the three-dimensional modeling of mineral deposits has made great progress, two-dimensional maps are still the most direct carrier for presenting the appearance of mineral deposits, and are also the main guiding maps used in the front line of exploration and development. However, the projection points of exploration projects on the plane will inevitably present discrete forms. How to use discrete exploration project information to depict the characteristics of mineral deposits on a specific plane as objectively as possible has become an urgent problem to be solved.

[0005] Previous solutions can be divided into two categories. One type of solution regards the geological and mineralization information of the entire ore deposit as a unified evolutionary set, and believes that it has an internally related spatial structure. This type of method usually uses geostatistical methods to perform two-dimensional or three-dimensional interpolation operations based on discrete data to draw relatively smooth graphics to show the characteristics of mineralization. For example, "Adaptive Fine Mining Method Based on Two-Dimensional Sections (CN201010234694)" uses "Boolean operations" to perform three-dimensional geological modeling, thereby achieving limited cutting of two-dimensional planes; Wei Hongxing of Beijing Nonferrous Metallurgical Design and Research Institute described the method of using Kriging method to calculate the plane distribution characteristics of ore bodies and conduct exploration and mining comparison in the article "Application and Practice of Kriging Method in Reserve Calculation of Shandak Copper Deposit".

[0006] The above methods have played a certain role, but they often cause problems or even mislead in actual application. For example, in exploration practice, boreholes without ore are often drilled in the delineated high anomaly areas; in the mining process, the prediction results of the above maps are often inconsistent with the actual mining situation. The root of the problem lies in the basic premise of geostatistics that "various geological characteristics related to geological bodies can be characterized by the spatial distribution characteristics of regionalized variables", that is, "the geological characteristics of ore deposits can be regarded as a unified system with space as the independent variable", which does not conform to geological reality.

[0007] In recent years, geological research has shown that the mineralization of endogenous deposits is a system of magma, structure, and fluid complex action, presenting an open, complex, fuzzy, nonlinear system with self-organized criticality. In other words, although the mineralization characteristics are constrained by underground space elements, they do not change regularly with space as the independent variable, but show a high degree of mutation and heterogeneity. After processing according to traditional geostatistical methods, on the one hand, smoothing makes local outliers inconspicuous, and on the other hand, false anomalies may appear in the neighboring areas of high outlier areas. In addition, geological and chronological studies have shown that the same spatial location often experiences more than several superimposed mineralizations. Therefore, the information of a mineral deposit at a certain spatial location should be a collection formed by multiple actions, rather than a unified system formed at a single time. The geostatistical method may forcibly associate geological information formed by different actions in different periods, thereby drawing misleading conclusions. In short, due to the contradiction between basic assumptions and geological facts, traditional geostatistical methods have organic defects.

[0008] In view of this, people have proposed the second type of idea to characterize the plane characteristics of the ore deposit, that is, to divide the two-dimensional space into independent finite elements and assign independent values ​​to characterize the appearance of the ore deposit. For example, in the article "Research on Simulation of Self-Organized Criticality of Tectonic-Fluid-Ore-Forming System", Xiong Yihui divided the two-dimensional space into three forms: hexagons, squares, and triangles, and constructed a "cellular automaton" to simulate the evolution of the ore-forming system. However, this division method does not take into account the actual distribution of the exploration project and can only be used for simulation calculations. "Mineral Deposit Drilling Information Processing and Analysis Method and Device Based on Finite Elements (CN114266179B)" proposes "On the exploration control plane map with control point positions, according to the principle of common edges of adjacent point control areas, determine the map reflecting the control range zoning of each exploration point (segmented finite element zoning map)". However, this scheme does not propose a systematic planar unit division method.

[0009] In summary, of the two existing approaches to using discrete exploration projects to characterize the planar characteristics of mineral deposits, the geostatistical approach has organic defects. The finite element division approach can more objectively characterize the planar characteristics of mineral deposits, but there is currently a lack of a systematic division method, and it has not been combined with exploration standards, so there are shortcomings in practical application. Summary of the invention

[0010] In order to overcome the shortcomings of the prior art, the object of the present invention is to provide a method for dividing a two-dimensional plane finite element of a mineral deposit with exploration engineering as the center.

[0011] To achieve the above objectives, the present invention provides the following solutions:

[0012] A method for dividing a two-dimensional plane finite element of a mineral deposit centered on exploration engineering, comprising:

[0013] Acquire exploration engineering data in the target area; the exploration engineering data at least includes the location of exploration engineering observation / sampling points and geological information collected at each observation / sampling point;

[0014] The plane where the two-dimensional map is to be drawn is selected from the target area as the section plane, and the intersection with the section plane is determined based on the location of the survey engineering observation / sampling point as the center point;

[0015] Taking the bid section plane as the reference, vertically extending the bid section height, the space formed is called the bid section;

[0016] Assign the geological information collected from the observation / sampling points within the bidding section to each center point to form a bidding section plane database;

[0017] A unit plane is formed by using a closed curve / broken line with each process center point as the center, and a process center area is constructed based on the unit plane;

[0018] The boundaries of each core area are determined to form a two-dimensional geological distribution map of the deposit.

[0019] Preferably, the geological information collected from the observation / sampling points within the bidding section is assigned to each center point to form a bidding section plane database, including:

[0020] Based on the footage information of the observation / sampling point, determine whether there is an observation / sampling point within the bidding section;

[0021] If there is no observation / sampling point within the bid section, the value is assigned as uncontrolled;

[0022] If there is one observation / sampling point within the bidding section, the geological information collected at the observation / sampling point will be assigned to the center point;

[0023] If the bid section includes multiple observation points, the mode value of the geological information recorded in the observation is taken and assigned to the center point. If the bid section involves multiple sampling points, the average value of the analysis results of the sampling points is taken and assigned to the center point.

[0024] Preferably, the process of forming a unit plane using closed curves / broken lines with each center point as the center and constructing the center area with the unit plane as the reference includes:

[0025] Screen out target mineralized elements that meet industrial utilization standards and their corresponding exploration project spacings; the exploration project spacings include: exploration project spacings for proven resources, exploration project spacings for controlled resources, and exploration project spacings for inferred resources;

[0026] Taking the center point as the center and expanding outward with a radius of half the survey project spacing to form a perfect circle area, the perfect circle area is taken as the initial center point area, and the boundary of the initial center point area is taken as the first-level enclosure line;

[0027] Determine whether the distance between the center points is smaller than the distance between the exploration projects. When the distance between the center points is smaller than the distance between the exploration projects, the initial center area generated based on the two center points must have two intersection points, and the connecting line of the two intersection points is used as the secondary enclosure line.

[0028] A new process center area is constructed based on the primary enclosing line and the secondary enclosing line.

[0029] Preferably, the constructing a new process center area based on the primary enclosing line and the secondary enclosing line includes:

[0030] Search the distance of the exploration project spacing along 360° from the center point outward. If there is a secondary enclosure line generated by the initial center area of ​​the center point in the search direction, the secondary enclosure line is used as the boundary of the new center area. If there are multiple secondary enclosure lines that meet the conditions, the center point closest to the corresponding center point is used as the boundary of the new center area. If there is no secondary enclosure line that meets the conditions, the primary enclosure line is used as the boundary of the new center area, and the area enclosed by each boundary is used as the scope of the new center area.

[0031] Assign the section plane data in the corresponding process center point to the entire new process center area;

[0032] The new Chengxin areas with equal contents of target mineralized elements are merged to complete the grading, forming a merged Chengxin area, and the outer boundary of the merged Chengxin area is used as the envelope line of the Chengxin area.

[0033] Preferably, the step of determining the boundary of each core area to form a two-dimensional geological distribution map of the ore deposit comprises:

[0034] When the distance between the Chengxin points in the merged Chengxin area is less than the distance between the proven exploration projects, a connecting line is established between the two points to form a triangle, and the triangles are merged to form a merged polygon, and the outer boundary of the merged polygon is used as the boundary of the Chengxin point area.

[0035] According to the specific embodiments provided by the present invention, the present invention discloses the following technical effects:

[0036] The present invention provides a method for dividing a mineral deposit into two-dimensional plane finite units with exploration engineering as the center. Compared with the prior art, the present invention can circumvent the organic defects of the geological statistical method and more accurately reflect the actual distribution of mineral deposit resources by dividing a mineralization system with high heterogeneity into independent units on the plane. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0038] Figure 1 A flow chart of a method for dividing two-dimensional plane finite elements of a mineral deposit centered on exploration engineering provided by the present invention;

[0039] Figure 2 A flowchart for constructing a process center area provided by the present invention;

[0040] Figure 3 A schematic diagram of the engineering space analysis and the construction of the center point distribution map provided by the present invention;

[0041] Figure 4 A schematic diagram of the boundary of a certain level of the process center point area and the envelope surface of the process center area provided by the present invention;

[0042] Figure 5 The results of the construction of the Chengxin area of ​​a copper polymetallic deposit provided by the present invention; (a) the results of the construction of the Chengxin point; (b) the results of the construction of the Chengxin point area; (c) the results of the construction of the Chengxin area; (d) the completed construction of the Chengxin area. La, Lb, Lc-the envelope of the Chengxin area at the level of proven, controlled, and inferred; BLa, BLb, BLc-the boundaries of the Chengxin area at the level of proven, controlled, and inferred. DETAILED DESCRIPTION

[0043] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. 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 creative work are within the scope of protection of the present invention.

[0044] Reference to "embodiments" herein means that a particular feature, structure, or characteristic described in conjunction with the embodiments may be included in at least one embodiment of the present application. The appearance of the phrase in various locations in the specification does not necessarily refer to the same embodiment, nor is it an independent or alternative embodiment that is mutually exclusive with other embodiments. It is explicitly and implicitly understood by those skilled in the art that the embodiments described herein may be combined with other embodiments.

[0045] The terms "first", "second", "third" and "fourth" in the specification and claims of the present application and the drawings are used to distinguish different objects rather than to describe a specific order. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions. For example, a series of steps, processes, methods, etc. are not limited to the listed steps, but may optionally include steps that are not listed, or may optionally include other step elements inherent to these processes, methods, products or devices.

[0046] In order to make the above-mentioned objectives, 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.

[0047] See also Figure 1-4 , a method for dividing the two-dimensional plane finite element of the ore deposit centered on the exploration project, including:

[0048] Step 1: Acquire exploration engineering data in the target area; the exploration engineering data at least includes the location of exploration engineering observation / sampling points and geological information collected at each observation / sampling point;

[0049] In practical applications, the present invention needs to fully collect exploration projects such as drilling within the working area, and organize the geospatial information of the project, including but not limited to the opening position, azimuth, inclination, total footage (depth), etc. Comprehensively organize the geological information obtained by the exploration project, including but not limited to the observation / sampling position, lithology, structural occurrence, deformation characteristics, element anomalies, mineralization grade, small weight, mechanical properties, alteration type, and degree of cryptoexplosion and fragmentation, to form an engineering space database (DBS: Database-space). The engineering space database includes engineering geospatial information, construction parameter information, and geological information.

[0050]

[0051] Among them, X n ,Y n ,A n ,Dip n ,Dept n They are respectively the X-axis coordinate, Y-axis coordinate, azimuth, inclination, and total footage of project n; Dep n1 is the footage of observation / sampling point n1, α n1 , β n1 Etc. are the corresponding geological information, and so on from n1, n2 to nm points.

[0052] Step 2: Select the plane where the two-dimensional map is to be drawn from the target area as the section plane, and determine the intersection with the section plane based on the location of the survey engineering observation / sampling point as the center point;

[0053] See also Figure 3 In the implementation of the present invention, it is necessary to select the plane on which the two-dimensional map is to be drawn as the section plane, conduct spatial analysis on the survey projects one by one, and determine the intersection point with the section plane, namely the engineering center point (ECP). After traversing all the survey projects, a projection map of the engineering center point (engineering center point) is formed.

[0054] Step 3: Taking the bid section plane as the reference, vertically extend the bid section height, and the space formed is called the bid section;

[0055] Taking into account the complexity of geological observation and sampling, the present invention can set a certain elevation interval above and below the section elevation, i.e., the section height (h), as the data source range for plane mapping according to actual exploration / mining needs or average sampling / observation elevation intervals.

[0056] Step 4: Assign the geological information collected from the observation / sampling points within the bidding section to each center point to form a bidding section plane database;

[0057] Furthermore, step 4 includes:

[0058] Based on the footage information of the observation / sampling point, determine whether there is an observation / sampling point within the bidding section;

[0059] If there is no observation / sampling point within the bid section, the value is assigned as uncontrolled;

[0060] If there is one observation / sampling point within the bidding section, the geological information collected at the observation / sampling point will be assigned to the center point;

[0061] If the bid section includes multiple observation points, the mode value of the geological information recorded in the observation is taken and assigned to the center point. If the bid section involves multiple sampling points, the average value of the analysis results of the sampling points is taken and assigned to the center point.

[0062] In the embodiment of the present invention, the footage information (Dep n1 -Dep nm ) to perform spatial coincidence analysis. If there is no observation / sampling point within the range h, the value is assigned to uncontrolled (null); if there is one observation / sampling point within the range h, the information is assigned to the center point; if there are multiple observation points within the range h, the mode value of the observation record information is taken; if multiple sampling points are involved, the average value of the results is taken and assigned to the center point. The obtained multi-dimensional information is assigned to each center point to form a section plane database (DBP: Database-plane), and the basic element structure of the database is as follows.

[0063] DBP={X k,Y k ,α k ,β k …}

[0064] Among them, X k , Y k is the spatial coordinate of the center point k, α k , β k etc. are the corresponding geological information.

[0065] Step 5: Using closed curves / broken lines to form a unit plane with each center point as the center, and constructing the center area based on the unit plane;

[0066] Furthermore, step 5 includes:

[0067] Screen out target mineralized elements that meet industrial utilization standards and their corresponding exploration project spacings; the exploration project spacings include: exploration project spacings for proven resources, exploration project spacings for controlled resources, and exploration project spacings for inferred resources;

[0068] The present invention needs to examine the enrichment of available elements based on basic analysis data. The mineralized elements that meet industrial utilization standards are screened out, and the exploration project spacing requirements (along the ore body trend direction) of different element mineralization types (ore species) are comprehensively compared. The exploration project spacing (2l a ), control level resource exploration project spacing (2l b ) and inferred resource exploration project spacing (2l c ) is used as the basic indicator for two-dimensional finite element division.

[0069] See also Figure 4 , Chengxinyuan, that is, FEC-FiniteElementbyEngineering Center. In space, it is a unit composed of closed curves / broken lines with ECP as the center. On a two-dimensional map, Chengxinyuan is expressed as a three-dimensional block with ECP as the center and a height of h. The intersection of Chengxinyuan and the section plane is the Chengxin area, that is, ECS: Engineering Centered Sector. Conversely, Chengxinyuan can also be defined as a three-dimensional unit formed by vertically extending h based on the Chengxin area. The basic steps for constructing the Chengxin area are as follows:

[0070] ① Taking the center point as the center, with half of the survey project spacing (l a , l b , l c ) is the radius to expand outward to form a perfect circular area, the perfect circular area is taken as the initial process center area, and the boundary of the initial process center area is taken as the first-level enclosure line;

[0071] ② Determine whether the distance between the center points (d-distance) is less than the distance between the survey projects. x ), then the initial center area generated based on the two center points must have two intersection points, and the connecting line of the two intersection points is used as the secondary enclosure line (En 2 -Envelope 2 );

[0072] ③ Search the distance of the exploration project spacing along 360° from the center point outward. If there is a secondary enclosure line generated by the initial center area of ​​the center point in the search direction, the secondary enclosure line is used as the boundary of the new center area; if there are multiple secondary enclosure lines that meet the conditions, the center point closest to the corresponding center point is used as the boundary of the new center area. If there is no secondary enclosure line that meets the conditions, the primary enclosure line is used as the boundary of the new center area, and the area enclosed by each boundary is used as the scope of the new center area;

[0073] ④ Assign the section plane data in the corresponding process center point to the entire new process center area;

[0074] ⑤ Merge the new Chengxin areas with equal content of target mineralized elements to complete the classification, forming a merged Chengxin area, and use the outer boundary of the merged Chengxin area as the envelope line of the Chengxin area.

[0075] Step 6: Determine the boundaries of each core area to form a two-dimensional geological distribution map of the deposit.

[0076] See also Figure 2 In step 6, when the distance between the center points in the merged center area is less than the distance between the proven exploration projects, a connecting line is established between the two points to form a triangle, and the triangles are merged to form a merged polygon. The outer boundary of the merged polygon is used as the boundary of the center point area (BL a ,BL b ,BL c ).

[0077] See also Figure 5 The following is a further explanation of a method of dividing a two-dimensional plane finite element of a ore deposit centered on an exploration project, taking a certain ore section of a copper polymetallic ore deposit as an example:

[0078] (1) Collect data and build a spatial database: The exploration data of the ore deposit was comprehensively sorted out, and the exploration data of 58 boreholes were collected. Based on this, the engineering spatial database DBS was constructed, and its basic structure is as follows:

[0079]

[0080] Among them, X n ,Yn ,A n ,Dip n ,Dept n They are respectively the X-axis coordinate, Y-axis coordinate, azimuth, inclination, and total footage of project n; Dep n1 is the footage of observation / sampling point n1, Cu n1 , Au n1 Etc. are geological information such as copper and gold grades, and so on from n1, n2 to nm points.

[0081] (2) Selection of bidding sections and mapping of the center point: The bidding section at an altitude of 3,500 meters was selected, and an intersection analysis was conducted with the project based on the spatial database to form a center point map of the 3,500-meter bidding section.

[0082] (3) Construction of bid section plane database: Based on the actual survey, especially the observation / sampling spacing and mining design indicators, 50m is determined as the bid section height (h). Spatial analysis is performed with the observation / sampling points in the engineering space database to construct the bid section plane database DBP:

[0083] DBP={X k ,Y k ,Cu k ,Au k …}

[0084] Where X k , Y k is the spatial coordinate of the center point k, Cu k , Au k Etc. are geological information such as copper grade and gold grade corresponding to the center point.

[0085] (4) Determination of two-dimensional unit division index: Based on the exploration industry standard of copper ore, the main mineral species, the interval between exploration projects at the proven level (2l a =80m), control level exploration project spacing (2l b =120m) and inferred level exploration project spacing (2l c =200m).

[0086] (5) Construction of the core area:

[0087] ① Based on the center point, a =40m,l b =60m,l c =100m as radius, construct the initial process center area and the first-level enclosure line;

[0088] ② Carry out spatial analysis on the primary enclosure lines of each level of control to form secondary enclosure lines;

[0089] ③Search for enclosing lines from each center point outward to form a center area;

[0090] ④ According to the bid section plane data DBP, assign the data information corresponding to the center point to the entire center area;

[0091] ⑤According to the outer enclosing lines of the boundary core area at different levels, the core area envelope line (L a ,L b ,L c ).

[0092] (6) Establishment of benchmark boundaries for each level of the project center area

[0093] According to the survey project spacing (2l a =80m, 2l b =120m, 2l c =200m) to screen and establish a connection grid. Taking the proven engineering spacing as an example, when the distance between two points is less than the proven engineering spacing (d<2l a ), a connecting line is created between the two points.

[0094] Based on the connection lines between the center points, construct the center point triangulation network of different levels. After merging all triangles, obtain the merged polygon and establish the center point area boundary (BL a ,BL b ,BL c ).

[0095] According to the specific embodiments of the present invention, the present invention discloses the following effects:

[0096] (1) Avoid the organic defects of geostatistical methods and divide the highly heterogeneous mineralized system into independent units on the plane, so as to present the two-dimensional geological and mineralized characteristics as completely and objectively as possible; (2) Make substantial innovations to the existing plane division methods and establish a detailed plane unit division method centered on exploration engineering; (3) Combine the objective characterization of ore deposit geological information with the evaluation of engineering control degree to facilitate industrial applications such as exploration and mining.

[0097] In this specification, each embodiment is described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the embodiments can be referred to each other. For the method disclosed in the embodiment, since it corresponds to the device disclosed in the embodiment, the description is relatively simple, and the relevant parts can be referred to the device part description.

[0098] The principles and implementation methods of the present invention are described in this article using specific examples. The description of the above embodiments is only used to help understand the method and core idea of ​​the present invention. At the same time, for those skilled in the art, according to the idea of ​​the present invention, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as limiting the present invention.

Claims

1. A method for dividing a two-dimensional plane finite element of a mineral deposit centered on an exploration project, characterized in that: include: Acquire exploration engineering data in the target area; the exploration engineering data at least includes the location of exploration engineering observation / sampling points and geological information collected at each observation / sampling point; The plane where the two-dimensional map is to be drawn is selected from the target area as the section plane, and the intersection with the section plane is determined based on the location of the survey engineering observation / sampling point as the center point; Taking the bid section plane as the reference, vertically extending the bid section height, the space formed is called the bid section; Assign the geological information collected from the observation / sampling points within the bidding section to each center point to form a bidding section plane database; A unit plane is formed by using a closed curve / broken line with each process center point as the center, and a process center area is constructed with the unit plane as the reference; The boundaries of each core area are determined to form a two-dimensional geological distribution map of the deposit.

2. A method for dividing a two-dimensional plane finite element of a mineral deposit centered on an exploration project as claimed in claim 1, characterized in that: The geological information collected from the observation / sampling points within the bidding section is assigned to each project center point to form a bidding section plane database, including: Based on the footage information of the observation / sampling point, determine whether there is an observation / sampling point within the bidding section; If there is no observation / sampling point within the bid section, the value is assigned as uncontrolled; If there is one observation / sampling point within the bidding section, the geological information collected at the observation / sampling point will be assigned to the center point; If the bid section includes multiple observation points, the mode value of the geological information recorded in the observation is taken and assigned to the center point. If the bid section involves multiple sampling points, the average value of the analysis results of the sampling points is taken and assigned to the center point.

3. A method for dividing a two-dimensional plane finite element of a mineral deposit centered on an exploration project as claimed in claim 2, characterized in that: The method of forming a unit plane using closed curves / broken lines with each center point as the center and constructing the center area with the unit plane as the reference includes: Screen out target mineralized elements that meet industrial utilization standards and their corresponding exploration project spacings; the exploration project spacings include: exploration project spacings for proven resources, exploration project spacings for controlled resources, and exploration project spacings for inferred resources; Taking the center point as the center and expanding outward with a radius of half the survey project spacing to form a perfect circle area, the perfect circle area is taken as the initial center point area, and the boundary of the initial center point area is taken as the first-level enclosure line; Determine whether the distance between the center points is smaller than the distance between the exploration projects. When the distance between the center points is smaller than the distance between the exploration projects, the initial center area generated based on the two center points must have two intersection points, and the connecting line of the two intersection points is used as the secondary enclosure line. A new process center area is constructed based on the primary enclosing line and the secondary enclosing line.

4. A method for dividing a two-dimensional plane finite element of a mineral deposit centered on an exploration project as claimed in claim 3, characterized in that: The constructing a new process center area based on the primary enclosing line and the secondary enclosing line includes: Search the distance of the exploration project spacing along 360° from the center point outward. If there is a secondary enclosure line generated by the initial center area of ​​the center point in the search direction, the secondary enclosure line is used as the boundary of the new center area. If there are multiple secondary enclosure lines that meet the conditions, the center point closest to the corresponding center point is used as the boundary of the new center area. If there is no secondary enclosure line that meets the conditions, the primary enclosure line is used as the boundary of the new center area, and the area enclosed by each boundary is used as the scope of the new center area. Assign the section plane data in the corresponding process center point to the entire new process center area; The new Chengxin areas with equal contents of target mineralized elements are merged to complete the grading, forming a merged Chengxin area, and the outer boundary of the merged Chengxin area is used as the envelope line of the Chengxin area.

5. A method for dividing a two-dimensional plane finite element of a mineral deposit centered on an exploration project as claimed in claim 4, characterized in that: The two-dimensional geological distribution map of the ore deposit formed by determining the boundary of each core area includes: When the distance between the Chengxin points in the merged Chengxin area is less than the distance between the proven exploration projects, a connecting line is established between the two points to form a triangle, and the triangles are merged to form a merged polygon, and the outer boundary of the merged polygon is used as the boundary of the Chengxin point area.

Images