This patent application claims the benefit and priority of Chinese Patent Application No. 202111586366.8, filed with the China National Intellectual Property Administration on Dec. 23, 2021, the disclosure of which is incorporated by reference herein in its entirety as part of the present application.
The present disclosure relates to the technical field of ore deposit drilling information analysis methods.
How to interpret and present ore deposit exploration information and guide the exploration, development and comprehensive utilization of ore resources is an important research content of ore deposit science. Compared with the basic geological survey, geochemical exploration, geophysical exploration, remote sensing, and prospecting projects carried out on the surface and shallow underground, the exploration borehole information can directly and objectively reflect geological information such as underground lithology, occurrence, structure, construction, mineralization, alteration, hidden explosion and crushing within the borehole range, providing direct, real and reliable data support for geological exploration. In addition, in recent years, with the development of analysis technology, the element analysis data of drilling core have been enriched.
In the prior art, the drilling exploration data is processed and presented mainly through the geological profile and geochemical interpolation map of the exploration line, which usually follows the processing sequence of the single project borehole histogram, exploration line profile, and sectional plan. This method has the following obvious defects:
(1) During construction of geological/geochemical profiles of exploration lines, due to the subjectivity of manual linking and the distortion or over fitting of interpolation methods, they often deviate from the geological reality. (2) The boundaries of different types of geological and geochemical maps are quite different, so it is often difficult to carry out spatial comprehensive matching analysis, and only rough superposition analysis can be carried out. (3) Due to the lack of vertical constraint between the two boreholes, the commonly used processing sequence of the profile has congenital errors, and the errors generated in the profile drawing phase will further accumulate during drawing of the plan, resulting in distortion.
On the other hand, in recent years, researchers have gradually realized that the geological system, especially the system of magma, structure and fluid interaction, is an open, complex and fuzzy nonlinear system with self-organized criticality. In addition, endogenetic deposits, especially fluid related deposits, often form complex superposition of multiple systems from different geological historical periods. It is difficult to interpret it as a whole under the existing technical conditions. Both empirical inference and interpolation based on existing data are likely to be misleading.
An objective of the present disclosure is to provide an improved ore deposit drilling information processing method and device, which can use computer technology to visually analyze and present sectional and layered geological borehole information. The processing method can solve the problem that it is difficult to perform comprehensive comparative analysis of multivariate geological information and geochemistry for the map results in the existing methods, and can further solve the problems of systematic errors caused by the lack of vertical constraints between boreholes, human factors having great influence on the delineation of geological body boundaries of exploration profiles, as well as the obvious interpolation distortion in the drawing of geochemical maps in the traditional method for obtaining the sectional plan through the exploration line profile.
A finite element-based ore deposit drilling information processing method includes:
obtaining borehole information, sampled sample information, and geological information corresponding to a borehole of a to-be-analyzed ore deposit, where the sampled sample information includes position information and chemical analysis information of samples sampled through the borehole;
based on three standard values of an associated grade, a cut-off grade, and an industrial grade specified in an industrial standard of mineral geological exploration specifications, performing grade classification of ore deposits corresponding to the samples at different positions according to the sampled sample information;
according to the borehole information, the sampled sample information, and the geological information corresponding to the borehole, performing vertical section analysis of the to-be-analyzed ore deposit, where the analysis includes: determining a spacing and number of vertical sections, calculating mineralization, hidden explosion, and alteration parameters of each of the sections, and classifying the calculated parameters;
obtaining a plane projection map of exploration borehole control points in different vertical sections, and performing map zoning according to a control range of each exploration point to form a sectional finite element zoning map;
according to results of the vertical section analysis, assigning different parameters to each finite element in the sectional finite element zoning map of different vertical sections to form a sectional finite element plan of each parameter;
performing three-dimensional (3D) superposition on the sectional finite element plan of each parameter according to spatial positions, and during the superposition, combining parameters at a same position to obtain a 3D composite data graph; and
scanning and analyzing the 3D composite data graph on a vertical profile to obtain a corresponding finite element profile containing composite data.
According to some implementations of the present disclosure, the borehole information includes position coordinates, an elevation, an azimuth angle, a dip angle and a depth of the borehole.
According to some implementations of the present disclosure, the sampled sample information includes starting position coordinates of sampling, ending position coordinates of sampling, and chemical analysis items and chemical analysis result data of the samples.
According to some implementations of the present disclosure, the geological information includes a name of rock and ore corresponding to the borehole, geological description, an alteration type, and a hidden explosion and crushing degree.
According to some implementations of the present disclosure, the grades include: grade A: greater than or equal to associated grade and less than cut-off grade; grade B: greater than or equal to cut-off grade and less than industrial grade; and grade C: greater than or equal to industrial grade.
According to some implementations of the present disclosure, the vertical section analysis specifically includes:
setting a vertical spacing of the section analysis;
according to the borehole information, obtaining an elevation H(max) of a maximum opening and an elevation H(min) of a minimum final hole of the borehole, and determining a number of the sections [N]+1 in combination with the vertical spacing
and
according to the geological information corresponding to the borehole, calculating a mineralization thickness percentage parameter, a thickness percentage parameter of a hidden explosion section, and a thickness percentage parameter of an altered section of the boreholes in each vertical section with different grades, and classifying the percentage parameters into several levels from high to low,
where the percentage parameters are calculated as follows:
mineralization thickness percentage=(mineralization thickness in vertical section/thickness of sample section)×100%;
thickness percentage of hidden explosion section=(length of hidden explosion section in vertical section/thickness of sample section)×100%; and
thickness percentage of altered section=(length of altered section in vertical section/thickness of sample section)×100%.
According to some implementations of the present disclosure, the vertical spacing is 20-25 m.
According to some implementations of the present disclosure, the levels include the following five levels: 100-80% for level I, 80-50% for level II, 50-30% for level III, 30-10% for level IV, and 10-0% for level V.
According to some implementations of the present disclosure, a process of obtaining the sectional finite element zoning map specifically includes:
in any vertical section, projecting the exploration borehole control points on the plan;
performing finite extrapolation of outermost control points and connecting an envelope of an extrapolated circle to determine a boundary of the finite element zoning map;
according to a principle of adjacent points sharing a same edge, dividing a boundary line of the finite elements along a midpoint between adjacent projects, and dividing the finite elements accordingly to form the sectional finite element zoning map;
according to a regional geological structure, a degree of engineering deflection, and distribution of discontinuities, and in combination with a geological expert knowledge base, adjusting boundary division of the finite elements; and
assigning the results of the vertical section analysis to each finite element class by class to form the sectional finite element plan of each parameter.
According to some implementations of the present disclosure, an extrapolated spacing is ½ of an average spacing of mining projects.
According to some implementations of the present disclosure, in the sectional finite element plan, data compositions at different points are as follows:
for any point Dm,
if the point is located in the finite element, there is:
D
m
={X
m
,Y
m
,P,α
p,βp,γp . . . }; and
if the point is located on a boundary of the finite element, there is:
D
m
={X
m
,Y
m
,H
n
,B},
where Xm and Ym represent plane coordinates of the point, P represents the finite element that the point belongs to on the plane, αp, βp, γp . . . represents various to-be-analyzed attributes to which the finite element is assigned, B represents a boundary point nature of the point, and Hn represents an elevation of a section where the point is located.
According to some implementations of the present disclosure, a process of obtaining the finite element profile containing composite data specifically includes:
on each horizontal plane of the 3D composite data graph, arranging a parallel vertical profile with a horizontal spacing of
when the vertical profile intersects a sectional plane of the sectional finite element zoning map, assigning an intersection point according to the combination of the parameters at the same position;
based on boundary points assigned in the sectional finite element plan, dividing an intersection line of the vertical profile and the sectional plane, extending each boundary point on the intersection line by ½
based on attribute points on the intersection line assigned with (P, αp, βp, γp), assigning the profile finite element, and constructing the analytic finite element profile and a corresponding database thereof.
According to the above information processing method, an analysis method can be further obtained, including: performing superposition analysis on the sectional finite element plan and/or analytical profile of the parameters.
According to some implementations of the present disclosure, the analysis method further includes: determining a spatial correlation of multiple objects according to results of the superposition analysis, delineating a potential metallogenic prospect zone in combination with a development trend of a mineralization grade, a hidden explosion grade, and an alteration grade and a geological law, and predicting and evaluating a resource potential of surrounding and deep ore bodies of a mining area.
According to the above information processing method and/or analysis method, an ore deposit drilling information analysis device can be further obtained, including a storage medium and an arithmetic unit (ALU). The storage medium stores program and/or structural data for implementing any one of the above information processing methods and/or the analysis method.
The present disclosure has the following beneficial effects:
Through the processing or analysis method of the present disclosure, a complete set of sectional planar projection analytic plan and hierarchical planar scanning analytic profile can be obtained, and the following can be realized:
(1) The geological and geochemical information of borehole exploration can be regularly processed to reduce the subjectivity of human interpretation and the distortion of interpolation analysis and more objectively and truly reflect underground stratum lithology, structure, mineralization, alteration, hidden explosion and crushing in the borehole range, and the maps presented are more consistent with the geological reality.
(2) By constructing the finite element only related to the spatial distribution of borehole control, and assigning the multivariate information to the finite element with standardized shape, it is convenient for geological and geochemical multidimensional analysis for the complex system with multi-phase superposition.
(3) Through the sectional finite element plan at different depths, spatial 3D co-position superposition, and scanning analytic finite element profile, and by making use of the more uniform distribution of the borehole control points on the plane than those on the profile, errors caused by the lack of vertical constraints between boreholes are reduced in the traditional interpretation method of the single project borehole histogram, exploration line profile, and sectional plan.
The present disclosure systematically solves the problems of few contents, single use, and poor metallogenic prediction effect of the current traditional geological borehole information expression method, which has excellent practical value and indicating significance for accurate delineation of prospecting target areas, accurate mining, dressing and smelting, and structural metallogenic research.
The present disclosure is described in detail below in combination with the embodiments and drawings, but it should be understood that the embodiments and drawings are only used to describe the present disclosure in an exemplary manner, and cannot constitute any limitation on the scope of protection of the present disclosure. All reasonable transformations and combinations included within the purpose of the present disclosure fall into the scope of protection of the present disclosure.
Some specific implementations of the present disclosure are described by taking the example of delineating favorable ore forming areas based on borehole mineralization, hidden explosion (blasting phenomena related to high-pressure fluid occurring in underground strata, usually related to fluid activity and mineralization) and alteration information.
Referring to the specification drawings
Step 1 Taking a geological work project or an ore deposit as a unit, map data, text data, and table data related to geological borehole data are comprehensively collected, and the collected geological borehole data is systematically checked and sorted out.
The collected data mainly includes: geological work report, engineering layout diagram, borehole sampling grade (basic) analysis data sheet, borehole engineering parameters (such as position parameters and inclinometer parameters), borehole histogram, and exploration line profile.
The Inspection May be as Follows.
(1) Integrity check: whether the geological borehole coordinates are complete, whether the borehole data amount is sufficient, whether sections of the single project are continuous, and whether the borehole sampling analysis data is complete.
(2) Validity check: whether the basic chemical analysis data of the borehole has undetected and incorrect numerical data, and whether the logical relationship between the analysis data items is correct.
(3) Normative check: whether the name, geological description, alteration type and crushing data of rock and ore are standardized and scientific.
Step 2 After step 1, the basic information, chemical analysis information and cataloging information of the borehole are sorted out.
The basic information of the borehole can include borehole construction parameters, such as the position coordinates X and Y of the hole, and the elevation H of the hole, and inclinometer parameters, such as the azimuth angle, dip angle and depth of the borehole.
The basic chemical analysis information of the borehole can include the sampling position in the borehole, such as the starting position and the ending position of the sampling, and analysis item names and analysis results.
The cataloging information can include the name of rock and ore, geological description, alteration type and the hidden explosion and crushing degree.
The above basic information can be further filled in and sorted out through the following basic borehole information table (Table 1), basic borehole chemical analysis information table (Table 2), and borehole cataloging information table (Table 3):
Step 3 After step 2, based on three standard values of an associated grade, a cut-off grade, and an industrial grade specified in an industrial standard of mineral geological exploration specifications, basic chemical analysis data of the borehole is classified into three grades.
More specifically, the three grades include: grade A: greater than or equal to associated grade and less than cut-off grade, grade B: greater than or equal to cut-off grade and less than industrial grade, and grade C: greater than or equal to industrial grade.
Step 4 A vertical spacing B of the section analysis and a number of sections are determined, and classification parameters of mineralization, hidden explosion, and alteration of the sections are calculated.
More specifically, the step can include the following sub-steps.
(1) A vertical spacing
(2) According to an elevation H(max) of a maximum opening and an elevation H(min) of a minimum final hole of the borehole and the vertical spacing
According to the borehole cataloging information and the three grades classified in step 3, the mineralization thickness percentage, the thickness percentage of the hidden explosion section, and the thickness percentage of the altered section of borehole grades (grade A or grade B or grade C) of each section are calculated, and the percentage data is classified into several levels from level I to level V from high to low, and filled in the following table.
Specific Parameters are Calculated as Follows:
mineralization thickness percentage=(mineralization thickness in section/thickness of sample section)×100%,
thickness percentage of hidden explosion section=(length of hidden explosion section within section/thickness of sample section)×100%, and
thickness percentage of altered section=(length of altered section within section/thickness of sample section)×100%.
The Classification can be Set as Follows:
100-80% for level I, 80-60% for level II, 60-40% for level III, 40-10% for level IV, and 10-0% for level V.
Step 5 On the exploration control plan with control point positions, according to the principle of control area of adjacent points sharing a same edge, the map reflecting the zoning of control area of each exploration point is determined, the sectional finite element zoning map is obtained, and values are assigned for each finite element in the zoning map according to the results obtained from sectional quantitative statistics in step 4 to form the sectional finite element plan (hereinafter referred to as the sectional plan).
The Step can More Specifically Include the Following Sub-Steps.
(1) The exploration borehole control points are projected on the sectional plans.
The exploration borehole drilled downward can be divided into inclined and vertical holes. The control point coordinates Xi and Yi of the vertical borehole on the sectional plan are the opening position coordinates X and Y in step 2. For the control point coordinates of the inclined hole on the sectional plan, according to the different borehole construction parameters, inclinometer parameters and elevation of each sectional plan in step 2, the control point coordinates Xi and Yi on the sectional plan can be corrected and assigned by depth section using the full-angle half-distance method with the opening position coordinates X and Y as the known starting points, as shown in
The full-angle half-distance method can be more specifically as follows:
where Xi and Yi are control point coordinates of the inclined hole on each sectional plan (i=0, 1, 2, . . . , N), Li is a length from the control point on each sectional plan to the opening positions X and Y, Zi is a dip angle of the control point on each sectional plan, Mi is an azimuth angle of the control point on each sectional plan, and dM is a magnetic declination. Li, Zi, and Mi can be checked in the inclinometer parameters in step 2. dM can be checked in the 1/50,000 topographic map of the borehole area.
(2) Finite extrapolation of outermost control points is performed and an envelope of an extrapolated circle is connected to determine a boundary of the finite element zoning map. The extrapolated spacing can be ½ of an average spacing of mining projects.
(3) According to a principle of adjacent points sharing a same edge, a boundary line of the finite elements is divided along a midpoint between adjacent projects, and the finite elements are divided accordingly to form the sectional finite element zoning map. The minimum scale of the sectional finite element zoning map can generally be 1/25,000.
(4) According to a regional geological structure, a degree of engineering deflection, and distribution of discontinuities, and in combination with a geological expert knowledge base, boundary division of the finite elements are moderately adjusted.
(5) The results in the sectional statistical table are assigned to each finite element class by class to form the sectional finite element plan and corresponding database.
Furthermore, since the division of the finite element in the sectional plan is only related to the position of the borehole control point, and is not affected by other variables, the finite element can contain multiple groups of information. Taking a point Dm on the sectional plan as an example, the structure of the database is described as follows.
For any point Dm inside the finite element, there is:
D
m
={X
m
,Y
m
,P,α
p,βp,γp . . . }
where Xm and Ym represent plane coordinates of the point, P represents the finite element that the point belongs to on the plane, which can be regarded as a fixed value after the division of the sectional finite element, and αp, βp, γp . . . represents various to-be-analyzed attributes to which the finite element is assigned.
For any point on a boundary of the finite element, there is:
D
m
={X
m
,Y
m
,H
m
,B},
where B represents a boundary point nature of the point, and Hn is an elevation of a section where the point is located.
According to the obtained sectional finite element plan and database, multiple types of geological information can be analyzed and presented separately or comprehensively, for example, the finite element plan of sectional and graded mineralization thickness, the finite element plan of sectional hidden explosion section thickness-alteration mineralization superposition, and the sectional mineralization thickness-hidden explosion section thickness-alteration section thickness superposition finite element plan.
Step 6 The sectional plans of various data are input into the 3D spatial database according to their spatial positions to form a composite database with co-position superposition.
More specifically, taking a point Dm in the n-th section from top to bottom of the 3D spatial database as an example, the structure of the database is illustrated. If the point is located inside a finite element, then:
D
m
={X
m
,Y
m
,H
n
,P,α
p,βp,γp . . . }
where Hn is the elevation of a section where the point is located (subject to the midpoint of the vertical spacing of the section), which satisfies:
H
n
=H(max)−
where H(max)≥Hn≥H(min)
If the point is located on a boundary of the finite element, there is:
D
m
={X
m
,Y
m
,H
n
,B}.
Step 7 The corresponding scanning analytic finite element profile (hereinafter referred to as analytic profile) is constructed using geographic information system (GIS) software such as Mapgis.
More specifically, the step can include the following sub-steps.
(1) A parallel vertical profile with a horizontal spacing of
(2) When the vertical profile intersects the sectional plane, it is assigned based on the co-position superposition database established in step 6. Taking a vertical profile Ll(Xl, Yl, Hl) as an example, when it coincides with the point Dm in the n-th section mentioned in step 6, that is, when (Xl=Xm, Yl=Ym, Hl=Hl), the coincidence point on the profile is assigned with (P, αp, βp, γp) or (B).
(3) Based on boundary points assigned in the previous step (that is, all points assigned with (B), an intersection line of the vertical profile and the sectional plane is divided. Thereafter, each boundary point on the intersection line is extended by ½
(4) Based on attribute points on the intersection line assigned with (P, αp, βp, γp), the profile finite element is assigned, and the analytic profile database is constructed.
Step 8 Superposition analysis of sectional plans or analytic profiles of different types of data is performed using multiple areal spatial analysis and superposition function of GIS software such as Mapgis. A spatial correlation of multiple objects is determined. A potential metallogenic prospect zone is delineated in combination with a development trend of a mineralization grade, a hidden explosion grade, and an alteration grade and a general geological law, and a resource potential of surrounding and deep ore bodies of a mining area is predicted and evaluated.
In the above process, the scanning and analysis of the finite element profile in steps 5 to 7 is shown in
According to the above specific implementations, taking a copper molybdenum deposit as an example, data related to geological borehole is comprehensively collected, and the collected geological borehole data is systematically checked and sorted out. A total of 78 boreholes is drilled in the ore deposit, with a total drill footage of 25,427.72 m. A total of 2,796 basic chemical analysis samples of lead is collected. A total of 13,074.93 m is drilled in the sampling section, accounting for 51.42% of the total drill footage.
In order to finely characterize the single element mineralization classification of Cu and Mo, based on three standard values of an associated grade, a cut-off grade, and an industrial grade specified in an industrial standard of mineral geological exploration specifications, basic chemical analysis data of the two mineralized elements of Cu and Mo in the borehole is divided into three grades from high to low (Table 5).
Through the processing of borehole engineering parameters and basic analysis data, the vertical spacing of section analysis is determined as 25 m. 78 exploration borehole control points are projected on each sectional plan with a spacing of 25 m. The mineralization classification parameters of different sections are calculated, and the mineralization thickness percentages are divided into level I to level V from high to low (Table 6).
On the exploration control plan with control point positions, according to the principle of control area of adjacent points sharing a same edge, the map reflecting the zoning of control area of each exploration point is determined, the sectional finite element zoning map is obtained, and values are assigned for each finite element in the zoning map according to the results obtained from sectional statistics to obtain the grade A mineralization sectional plan of Cu and Mo elements as shown in
The sectional plans of various data are input into the 3D spatial database according to their spatial positions to form a composite database with co-position superposition. The analytical profiles of grade A mineralization of Cu and Mo elements shown in
Superposition analysis of sectional plans or analytic profiles is performed using multiple areal spatial analysis and superposition function of GIS software such as Mapgis.
The sectional plan of
The analytic profiles of
The above embodiments merely describe preferred implementations of the present disclosure, and the scope of protection of the present disclosure is not limited to the above embodiments. All technical solutions under the idea of the present disclosure belong to the scope of protection of the present disclosure. It should be noted that those of ordinary skill in the art can make several improvements and modifications without departing from the principles of the present disclosure. These improvements and modifications should also be considered as falling within the scope of protection of the present disclosure.
Number | Date | Country | Kind |
---|---|---|---|
202111586366.8 | Dec 2021 | CN | national |