Multi-metal mineral exploration method and system

Abstract

The invention discloses a multi-metal mineral exploration method and system, and relates to the field of multi-metal mineral exploration, and the system comprises a data acquisition module, a structure analysis module, a geophysical exploration module and a mineral evaluation module. Comprising the following steps: collecting wrinkle information, fault information and geophysical prospecting information, carrying out wrinkle state and fault state analysis on the wrinkle information and the fault information to obtain a structural shape prospecting value of each region, comparing the structural shape prospecting value of each region with a set comparison threshold to obtain a judgment signal of each region, receiving a geophysical prospecting signal, and carrying out geophysical prospecting. And performing gravity anomaly, magnetic anomaly and electrical difference analysis on the high-possibility mineral product existence area to obtain a physical exploration shape value of each geophysical exploration area, and performing comparative analysis on the physical exploration shape value and a set reference standard interval to obtain a set of mineral product areas of different grades corresponding to different geophysical exploration areas, so that the multi-metal mineral product areas can be graded according to mineral product reserves. And a basis is provided for making multi-metal mineral exploration and mine development design.

Description

Technical Field

[0001] The present invention relates to the field of polymetallic mineral exploration, and in particular to a polymetallic mineral exploration method and system. Background Art

[0002] With the continuous development and growth of industry, the demand for polymetallic minerals such as copper, zinc, and lead continues to increase. Traditional polymetallic mineral exploration often relies on simple detection methods and manual field surveys, which are inefficient and cannot meet the continuously increasing demand for polymetallic minerals. At the same time, the accuracy of traditional methods cannot be guaranteed, so polymetallic mineral exploration methods and systems came into being.

[0003] In the process of exploration of polymetallic minerals, it is often impossible to conduct targeted exploration, which leads to increased exploration costs and waste of time. In addition, it is impossible to classify polymetallic mineral areas according to mineral reserves, and it is impossible to provide a basis for polymetallic mineral mining plans and mine designs. Therefore, multi-faceted targeted exploration and classification of mineral areas according to mineral reserves have become issues that need to be resolved.

[0004] In order to solve the above defects, a technical solution is now provided. Summary of the invention

[0005] In order to solve the technical problems raised by the above background technology, the present invention is proposed. The embodiment of the present invention provides a method and system for multi-metal mineral exploration.

[0006] The purpose of the present invention can be achieved through the following technical solutions:

[0007] In a first aspect, the present invention provides a method for prospecting for polymetallic minerals, comprising the following steps:

[0008] S1: Collect fold information, fault information and geophysical information;

[0009] S2: By analyzing the fold status and fault status of each fold information and fault information, the structural survey value of each area is obtained;

[0010] S3: Compare the wrinkle status value of each area with the set comparison threshold to obtain the judgment signal of each area;

[0011] S4: Receive geophysical exploration signals, analyze gravity anomalies, magnetic anomalies and electrical differences in areas with high probability of mineral deposits, and obtain geophysical exploration status values ​​of each geophysical exploration area;

[0012] S5: Receive the geophysical status values ​​of each geophysical exploration area, and compare and analyze them with the set reference standard interval to obtain different levels of mineral area sets corresponding to different geophysical exploration areas.

[0013] Furthermore, the analysis steps of the structural status values ​​of each region are as follows:

[0014] Obtain the dip angle of each rock layer in the region, perform average processing, and obtain the average dip angle pj of each rock layer in each region. Calculate the shape value zkz of each fold by combining with the fold amplitude value fz of each region, the complex shape value fx of each fold fissure and the connectivity value lt of each region.

[0015] The fold state survey value of each region and the fault state survey value of each region are multiplied by the correction factor coefficients d1 and d2 respectively, and then summed up to obtain the structural state survey value of each region.

[0016] Furthermore, the analysis steps of the complex shape value of the fold cracks in each region and the connectivity value of the cracks in each region are as follows:

[0017] The folded rock mass is scanned by three-dimensional scanning. According to the grayscale information and the spatial distribution of the point cloud data, the background point cloud is separated from the crack point cloud by cluster analysis. The center line of the fold crack is extracted by interpolation algorithm. The center line of the crack is placed in the three-dimensional coordinate system. The center line of the crack is completely covered by spheres of different radii, starting from radius R1, R1, R2, R3, R4, R5, R6, R7, R8, R9, R10, R11, R12, R13, R14, R15, R16, R17, R18, R19, R20, R21, R22, R23, R24, R25, R26, R27, R28, R29, R30, R31, R32, R33, R34, R35, R36, R37, R38, R39, R40, R 1 / 2 , R1 1 / 4 , R1 1 / 8 , R1 1 / 16 , ..., R1 1 / 2×N , record the number of spheres gs corresponding to each radius value, establish a two-dimensional coordinate system with logR as the horizontal axis and gs as the vertical axis, and perform linear fitting on the coordinate points in the coordinate system to obtain the slope of the fitting line and take the absolute value as the complex value of the fold cracks. Sum the complex values ​​of each fold crack and divide it by the number of fold cracks to obtain the average value of the complex values ​​of each fold crack and mark it as the complex value fx of the fold cracks in each region;

[0018] Put the center lines of each crack obtained in S4 into the three-dimensional coordinate system, and obtain the number of intersections between crack i and crack j in turn, mark them as the node number Kij, obtain the angle θij between the directions of crack i and crack j, and obtain the minimum spatial distance Dij between crack i and crack j, and analyze to obtain the connectivity value lt of the cracks in each area.

[0019] Furthermore, the analysis steps of the wrinkle amplitude value of each region are as follows:

[0020] The exploration area is divided into several areas. In each area, the position of the core layer of the same rock layer on both wings in the magnetic direction is obtained by using a magnetic positioning and orientation instrument. The angles θy and θe are obtained. θy and θe are the respective strikes of the same rock layer on both wings. The angles between the core layer and the core cylinder axis of the same rock layer on both wings are obtained, which are marked as the inclination angles αy and αe of the same rock layer on both wings, respectively. The core is obtained by vertical drilling. When the rock layer inclination is on the same side as the angle between the core cylinder axis, the inclination angles are δy=90°-αy and δe=90°-αe, respectively. When the rock layer inclination is on the opposite side to the angle between the core cylinder axis, the inclination angles are δy=90°+αy and δe=90°+αe, respectively.

[0021] The laser radar emits a laser beam in the survey area. After the laser beam interacts with the rock surface, it is received by the receiving device. The cloud data obtained is filtered to remove noise, and the fold feature points of the rock formation are extracted and converted into the engineering coordinate system to obtain the coordinate points of the fold core and the same rock formation position points on both wings. Then, the horizontal distance between the same rock formation position points on one wing and the same rock formation position points on the other wing is obtained by analysis, which are marked as Ly and Le respectively.

[0022] A three-dimensional rectangular coordinate system is established with the strike direction of the rock stratum as the x-axis, the direction perpendicular to the strike and horizontal as the y-axis, the vertical direction as the z-axis, and the fold core as the origin (0, 0, 0). The position points of the same rock stratum on both wings are determined as the first point coordinates (x1, y1, z1) and the second point coordinates (x2, y2, z2), respectively. The coordinate point values ​​of the first point coordinates and the second point coordinates are analyzed and obtained, and the amplitude value fd of the folded rock stratum is calculated. The amplitude values ​​of the folds at different cross-sectional positions are calculated, and the amplitude values ​​at different cross-sectional positions are summed and divided by the number of cross-sectional positions to obtain the average amplitude value of the folds at different cross-sectional positions, which is marked as the fold amplitude value fz in each region.

[0023] Furthermore, the analysis steps of the fault profile of each region are as follows:

[0024] By observing the mineral crystals of each fault rock under a transmission electron microscope in each region, the total length of the dislocation line within a certain region of the mineral crystal is obtained, and the dislocation density value wm of the fault rock mineral is obtained by analysis. The orientation difference between two adjacent scanning points of the mineral crystal of the fault rock is obtained through an electron microscope. When the orientation difference is greater than or equal to the set orientation difference threshold, a subgrain boundary is formed between the two adjacent scanning points, and the crystal orientation difference between adjacent scanning points inside each subgrain boundary is obtained, marked as each subgrain boundary angle value, and each subgrain boundary angle value is marked as each subgrain boundary angle value. The average value of the grain boundary angle is processed to obtain the average subgrain boundary angle, the number of subgrain boundaries and the average spacing of each adjacent subgrain boundary are obtained, the average subgrain boundary angle is multiplied by the number of subgrain boundaries to obtain the subgrain boundary product value, the subgrain boundary product value is divided by the average spacing of adjacent subgrain boundaries and multiplied by the correction factor coefficient to obtain the subgrain boundary shape value of the fault rock mineral, the fault displacement, the width of the fault fracture zone, the dislocation density value of the fault rock mineral and the subgrain boundary shape value of the fault rock mineral are obtained, and the fault shape exploration value of each area is obtained by analysis.

[0025] Furthermore, the analysis steps of the region determination signals are as follows:

[0026] The structural exploration value of each area is compared and analyzed with the comparison threshold XXT1. When the structural exploration value of each area is less than the comparison threshold XXT1, the corresponding area is a low-probability mineral existence area. When the structural exploration value of each area is greater than or equal to the comparison threshold XXT1, the corresponding area is a high-probability mineral existence area, and a geophysical exploration signal is generated.

[0027] Furthermore, the analysis steps of the geophysical exploration status values ​​of each geophysical exploration area are as follows:

[0028] The low resistivity ratio dzb, the minimum resistivity zd, the magnetic field intensity amplitude fq and the re-evaluation value zpz are calculated to obtain the geophysical exploration status wkz of each geophysical exploration area.

[0029] Furthermore, the analysis steps of the low resistivity ratio, minimum resistivity, magnetic field intensity amplitude value and re-evaluation value are as follows:

[0030] A square grid is used in the area with high possibility of mineral existence, in which the grid is set at a certain interval. The measured gravity value of each grid center is measured by a gravimeter. According to the geographical latitude and elevation value of the location of each grid center, the standard gravity value is obtained. The measured gravity value of each grid center is subtracted from the standard gravity value to obtain the gravity anomaly value of each grid center. With the position distance of each grid center as the horizontal axis and the gravity anomaly value as the vertical axis, a curve graph of the gravity anomaly value changing with the position distance is established. The two coordinate points corresponding to the gravity amplitude on both sides of each peak value dropping to half of the peak value are obtained, and the distance between the two coordinate points is taken and marked as the gravity anomaly half-width value. The average gravity anomaly half-width value pb of each peak value is taken, and the number of peaks fg and the average spacing of each peak pj are counted. The peak value of each peak is obtained and compared with the set large mineral peak value to obtain the number of peaks greater than the set large mineral peak value. The value is divided by the total number of peaks to obtain the large peak proportion value df, and the re-evaluation value zpz of each area is obtained by analysis;

[0031] Electrodes are set at a certain distance in the area with high possibility of mineral existence. The resistivity of each measuring point is obtained by a resistivity meter, and the resistivity is compared with the set resistivity threshold. When the resistivity is less than the set resistivity threshold, the measuring point is defined as a low-resistance abnormal point. The proportion of low-resistance abnormal points in the total number of measurements is counted and marked as the low-resistance proportion dzb. The lowest resistivity zd of the measuring point is obtained. The magnetometer is used to measure on a regular grid in the area with high possibility of mineral existence to obtain the magnetic field strength of each measuring point, and the maximum value of the magnetic field strength of each measuring point is taken and marked as the magnetic field strength amplitude value fq.

[0032] Furthermore, the analysis steps of the different geophysical exploration areas corresponding to the different grades of mineral area sets are as follows:

[0033] When the geophysical condition value of each geophysical exploration area is within the preset gradient reference standard interval XX0, the corresponding area is classified as the first-level large mineral area set W0; when the geophysical condition value of the geophysical exploration area is within the preset gradient reference standard interval XX1, the corresponding area is classified as the second-level medium-sized mineral area set W1; when the geophysical condition value of the geophysical exploration area is within the preset gradient reference standard interval XX2, the corresponding area is classified as the third-level small mineral area set W2; when the geophysical condition value of each geophysical exploration area is within the preset gradient reference standard interval XX3, the corresponding area is classified as the fourth-level non-ore abnormal area set W3, and the level set of each geophysical exploration area is divided and sent to the system display platform for display.

[0034] In a second aspect, the present invention provides a multi-metal mineral exploration system, characterized in that it comprises:

[0035] The data acquisition module is used to collect fold information, fault information and geophysical information, and send it to the structural analysis module and the geophysical exploration module;

[0036] The structural analysis module is used to receive the information of each fold and fault, analyze the fold state and fault state, obtain the structural state value of each area, and compare it with the set comparison threshold to obtain the judgment signal of each area;

[0037] The geophysical exploration module is used to receive geophysical exploration signals, analyze gravity anomalies, magnetic anomalies and electrical differences in areas with high probability of mineral deposits, and obtain geophysical exploration status values ​​of each geophysical exploration area;

[0038] The mineral assessment module is used to receive the geophysical conditions of each geophysical exploration area, and compare and analyze them with the set reference standard interval to obtain a set of mineral areas of different grades corresponding to different geophysical exploration areas.

[0039] Compared with the prior art, the present invention has the following beneficial effects:

[0040] 1. The present invention analyzes the fold state and fault state of fold information and fault information to obtain the average value of the rock layer inclination, the fold amplitude value, the complex shape value of the fold fissure, the connectivity value of the fissure and the fault state exploration value of each region, and then analyzes to obtain the structural state exploration value of each region, and compares it with the set comparison threshold to obtain the judgment signal of each region, receives the geophysical exploration signal, and analyzes the gravity anomaly, magnetic anomaly and electrical difference of the area with high possibility of mineral existence to obtain the geophysical exploration state value of each geophysical exploration area. It can first analyze the rock layer, folds and faults in each region, analyze the possibility of mineral existence area, and further conduct geophysical exploration analysis on the area with high possibility, thereby reducing the cost of exploration and reducing the waste of time.

[0041] 2. The present invention receives the geophysical status values ​​of each geophysical exploration area, and compares and analyzes them with the set reference standard interval to obtain a set of mineral areas of different grades corresponding to different geophysical exploration areas. It can classify polymetallic mineral areas according to mineral reserves, and provide a basis for formulating polymetallic mineral exploration and mine development design. BRIEF DESCRIPTION OF THE DRAWINGS

[0042] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the following briefly introduces the drawings required for describing the embodiments. The following drawings are not intentionally scaled according to the actual sizes, and the focus is on illustrating the main purpose of the present invention.

[0043] Figure 1 Schematic diagram of the method of the present invention

[0044] Figure 2 This is the overall system block diagram of the present invention. DETAILED DESCRIPTION

[0045] The technical solutions in the embodiments of the present invention will be described clearly and completely below in conjunction with the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work also fall within the scope of protection of the present invention.

[0046] like Figure 1 - Figure 2 As shown, a polymetallic mineral exploration system includes a data acquisition module, a structural analysis module, a geophysical exploration module and a mineral assessment module.

[0047] The data acquisition module is used to collect fold information, fault information and geophysical information, and send it to the structural analysis module and the geophysical exploration module.

[0048] The structural analysis module is used to receive the information of each fold and fault, analyze the fold state and fault state, obtain the structural state value of each area, and compare it with the set comparison threshold to obtain the judgment signal of each area, where the fold information includes rock formation information and fracture information. The specific analysis steps are as follows:

[0049] The specific steps for solving the wrinkle shape value are as follows:

[0050] Step 1: Divide the exploration area into several areas, and obtain the position of the core layer of the same rock layer on both wings in the magnetic direction in each area through the magnetic positioning and orientation instrument, and obtain the angles θy and θe, θy and θe are the respective strikes of the same rock layer on both wings, and obtain the angle between the core layer and the core cylinder axis of the same rock layer on both wings, which are marked as the inclination angles αy and αe of the position of the same rock layer on both wings, respectively. The core is obtained by vertical drilling. When the rock layer inclination is on the same side as the angle between the core cylinder axis, the inclination angles are δy=90°-αy and δe=90°-αe respectively, and when the rock layer inclination is on the opposite side of the angle between the core cylinder axis, the inclination angles are δy=90°+αy and δe=90°+αe respectively;

[0051] Step 2: Use laser radar to emit a laser beam in the survey area. After the laser beam interacts with the rock surface, it is received by the receiving device. The obtained cloud data is filtered to remove noise, and the fold feature points of the rock formation are extracted and converted into the engineering coordinate system to obtain the coordinate points of the fold core and the same rock formation position points on both wings. The abscissa point of the fold core is subtracted from the abscissa point of the same rock formation position point on one wing and squared to obtain the square value of the abscissa. The ordinate point of the fold core is subtracted from the ordinate point of the same rock formation position point on one wing and squared to obtain the square value of the ordinate. The square value of the abscissa and the square value of the ordinate are summed and the square root is taken to obtain the horizontal distance of the same rock formation position point on one wing. Repeat the above operation to obtain the horizontal distance of the same rock formation position point on the other wing, which are marked as Ly and Le respectively.

[0052] Step 3: Establish a three-dimensional rectangular coordinate system with the strike direction of the rock layer as the x-axis, the direction perpendicular to the strike and horizontal as the y-axis, and the vertical direction as the z-axis, with the fold core as the origin (0, 0, 0). Determine the position points of the same rock layer on both wings as the first point coordinates (x1, y1, z1) and the second point coordinates (x2, y2, z2). According to the formula and Get the coordinate point values ​​of the first point and the second point according to the set formula The amplitude value fd of the folded rock layer is obtained, where λ is the set correction factor coefficient, which is used to improve the accuracy of the calculation results. The specific value is determined by personnel in this professional field. The amplitude values ​​of the folds at different cross-sectional positions are obtained, and the amplitude values ​​at different cross-sectional positions are summed and divided by the number of cross-sectional positions to obtain the average amplitude value of the folds at different cross-sectional positions, and marked as the fold amplitude value fz of each region;

[0053] Step 4: Scan the folded rock mass through 3D scanning. According to the grayscale information and the spatial distribution of the point cloud data, the background point cloud is separated from the crack point cloud by cluster analysis. The center line of the fold crack is extracted by interpolation algorithm. The center line of the crack is placed in the 3D coordinate system. Spheres with different radii are used to completely cover the center line of the crack. Starting from radius R1, R1, R2 and R3 are selected in turn. 1 / 2 , R1 1 / 4 , R1 1 / 8 , R1 1 / 16 , ..., R1 1 / 2×N , record the number of spheres gs corresponding to each radius value, establish a two-dimensional coordinate system with logR as the horizontal axis and gs as the vertical axis, and perform linear fitting on the coordinate points in the coordinate system to obtain the slope of the fitting line and take the absolute value as the complex value of the fold cracks. Sum the complex values ​​of each fold crack and divide it by the number of fold cracks to obtain the average value of the complex values ​​of each fold crack and mark it as the complex value fx of the fold cracks in each region;

[0054] Step 5: Place the center lines of each crack obtained in step 4 into the three-dimensional coordinate system, obtain the number of intersections between crack i and crack j in turn, mark them as the node number Kij, obtain the angle θij between the directions of crack i and crack j, θij takes the value of [0, π], and obtain the minimum spatial distance Dij between crack i and crack j. When there is an intersection, Dij = 0, i = 1, 2, 3, ..., N, N is the maximum value of the crack number, according to the set formula The connectivity value lt of the cracks in each region is obtained, where to1, to2, to3, to4, to5, to6, to7 and to8 are all set weight factor coefficients to improve the accuracy of the calculation, to4>to1>to6, to7>to2, to8>to3>to5, and c is a natural constant with a specific value of 5.125;

[0055] Step 6: Obtain the inclination of each rock layer in the region, perform average processing, and obtain the average inclination of the rock layer in each region, marked as pj, and normalize it with the fold amplitude value fz of each region, the complex shape value fx of the fold fissure of each region, and the connectivity value lt of each region. Substitute it into the set formula zkz=f1×pj+f2×fz+f3×fx+f4×lt to calculate the shape value zkz of the fold in each region, where f1, f2, f3 and f 4 are the weight factor coefficients of the mean dip angle of each region, the fold amplitude value of each region, the complex shape value of each region's fold fissures and the connectivity value of each region's fissures, and the specific values ​​are 1.12, 2.13, 3.12 and 2.15, respectively. It should be noted that, as can be seen from the formula, the greater the dip angle of the rock formation, the greater the fold amplitude value, the greater the complex shape value of the fold fissures, the greater the connectivity value of the fissures, the greater the fold shape value, and the greater the possibility of the existence of polymetallic minerals;

[0056] By observing the mineral crystals of each fault rock under a transmission electron microscope in each region, the total length of the dislocation line within a certain region of the mineral crystal is obtained, marked as zl, where the area of ​​the certain range is S and the thickness is hd. According to the set formula wm=zl / (S×hd), the dislocation density value wm of the fault rock mineral is obtained. The orientation difference between two adjacent scanning points of the mineral crystal of the fault rock is obtained through an electron microscope. When the orientation difference is greater than or equal to the set orientation difference threshold, a subgrain boundary is formed between the two adjacent scanning points, and the crystal orientation difference between adjacent scanning points inside each subgrain boundary is obtained, marked as each The subgrain boundary angle value is obtained by processing the average value of each subgrain boundary angle value to obtain the mean value of the subgrain boundary angle, obtaining the number of subgrain boundaries and the average spacing of each adjacent subgrain boundary, multiplying the mean value of the subgrain boundary angle with the number of subgrain boundaries to obtain the subgrain boundary product value, dividing the subgrain boundary product value by the average spacing of adjacent subgrain boundaries, and multiplying it by the correction factor coefficient to obtain the subgrain boundary shape value of the fault rock mineral, obtaining the fault displacement, the width of the fault fracture zone, and the dislocation density value of the fault rock mineral and the subgrain boundary shape value of the fault rock mineral, performing weighted calculation, and multiplying it by the corresponding weight factor coefficient to obtain the fault shape exploration value of each area;

[0057] It should be noted that the fault displacement refers to the sum of the horizontal and vertical displacement distances through the same marker layer on both sides of the fault;

[0058] The fold state survey value of each region and the fault state survey value of each region are multiplied by the correction factor coefficients d1 and d2 respectively, and then summed to obtain the structural state survey value of each region;

[0059] Compare and analyze the structural exploration value of each area with the comparison threshold XXT1. When the structural exploration value of each area is less than the comparison threshold XXT1, the corresponding area is a low-possibility mineral existence area. When the structural exploration value of each area is greater than or equal to the comparison threshold XXT1, the corresponding area is a high-possibility mineral existence area, and a geophysical exploration signal is generated.

[0060] The geophysical exploration module is used to receive geophysical exploration signals, analyze gravity anomalies, magnetic anomalies and electrical differences in areas with high probability of mineral deposits, and obtain the geophysical status values ​​of each geophysical exploration area. The specific analysis steps are as follows:

[0061] A square grid is used in the area with high possibility of mineral existence, in which the grid is set at a certain interval. The measured gravity value of each grid center is measured by a gravimeter. According to the geographical latitude and elevation value of the location of each grid center, the standard gravity value is obtained. The measured gravity value of each grid center is subtracted from the standard gravity value to obtain the gravity anomaly value of each grid center. With the location distance of each grid center as the horizontal axis and the gravity anomaly value as the vertical axis, a curve graph of the gravity anomaly value versus location distance is established to obtain the corresponding gravity amplitude on both sides of each peak value when it drops to half of the peak value. The two coordinate points are taken, and the distance between the two coordinate points is taken, marked as the half-width value of the gravity anomaly, the average half-width value of the gravity anomaly pb of each peak is taken, and the number of peaks fg and the average spacing of each peak pj are counted, the peak value of each peak is obtained, and it is compared with the set large mineral peak value, and the number of peaks greater than the set large mineral peak value is obtained, and the value is divided by the total number of peaks to obtain the large peak proportion df, and the average gravity anomaly half-width value pb, the number of peaks fg, the average spacing of peaks pj and the large peak proportion df are normalized, according to the set formula The reassessment value zpz of each region is obtained, where XOX1, XOX2, XOX3 and XOX4 are the set weight factor coefficients of the average gravity anomaly half-width, the number of peaks, the average spacing between peaks and the proportion of large peaks, respectively. The specific values ​​are determined by professionals in this field, and e is a natural constant with a value of 2.718;

[0062] Electrodes are set at a certain distance in the area with high possibility of mineral existence, and the resistivity of each measuring point is obtained by a resistivity meter. The resistivity is compared with the set resistivity threshold. When the resistivity is less than the set resistivity threshold, the measuring point is regarded as a low resistivity abnormal point. The proportion of low resistivity abnormal points in the total number of measurements is counted and marked as the low resistivity proportion dzb. The lowest resistivity zd of the measuring point is obtained. The magnetometer is used to measure on a regular grid in the area with high possibility of mineral existence to obtain the magnetic field strength of each measuring point, and the maximum value of the magnetic field strength of each measuring point is taken and marked as the magnetic field strength amplitude value fq;

[0063] The low resistivity ratio dzb, the minimum resistivity zd, the magnetic field intensity amplitude fq and the re-evaluation value zpz are normalized according to the set formula.

[0064] The geophysical condition value wkz of each geophysical exploration area is calculated, where LP1, LP2, LP3 and LP4 are the preset weight factors of the proportion of low resistivity, the lowest resistivity, the amplitude value of magnetic field intensity and the re-evaluation value, and the values ​​are 1.122, 2.215, 3.214 and 3.21 respectively, and ρ is the preset correction factor, and the specific value is 0.651. It should be noted that when the resistivity meter is measured, the more metal minerals, the lower the resistivity value, and when the magnetometer is measured, the more metal minerals, the greater the magnetic field intensity;

[0065] The mineral assessment module is used to receive the geophysical conditions of each geophysical exploration area, and compare and analyze them with the set reference standard interval to obtain different levels of mineral area sets corresponding to different geophysical exploration areas. The specific analysis is as follows:

[0066] Set the reference standard intervals XX0, XX1, XX2, and XX3 for the geophysical condition values ​​of each geophysical exploration area, and substitute the geophysical condition values ​​of each geophysical exploration area into the preset gradient reference intervals XX0, XX1, XX2, and XX3 for comparative analysis, wherein the interval values ​​of XX0, XX1, XX2, and XX3 decrease in a gradient manner;

[0067] When the geophysical condition value of each geophysical exploration area is within the preset gradient reference standard interval XX0, the corresponding area is classified as a first-level large-scale mineral area set W0; when the geophysical condition value of the geophysical exploration area is within the preset gradient reference standard interval XX1, the corresponding area is classified as a second-level medium-sized mineral area set W1; when the geophysical condition value of the geophysical exploration area is within the preset gradient reference standard interval XX2, the corresponding area is classified as a third-level small-scale mineral area set W2; when the geophysical condition value of each geophysical exploration area is within the preset gradient reference standard interval XX3, the corresponding area is classified as a fourth-level non-ore abnormal area set W3, and the level set of each geophysical exploration area is divided and sent to the system display platform for display, providing a basis for formulating polymetallic mineral mining plans and conducting mine design;

[0068] A method for prospecting polymetallic minerals comprises the following steps:

[0069] S1: Collect fold information, fault information and geophysical information

[0070] S2: By analyzing the fold status and fault status of each fold information and fault information, the structural survey value of each area is obtained;

[0071] S3: Compare the structural status value of each area with the set comparison threshold to obtain the judgment signal of each area;

[0072] S4: Receive geophysical exploration signals, analyze gravity anomalies, magnetic anomalies and electrical differences in areas with high probability of mineral deposits, and obtain geophysical exploration status values ​​of each geophysical exploration area;

[0073] S5: Receive the geophysical status values ​​of each geophysical exploration area, and compare and analyze them with the set reference standard interval to obtain different levels of mineral area sets corresponding to different geophysical exploration areas.

[0074] The above is an explanation of the present invention and should not be considered as a limitation thereof. Although several exemplary embodiments of the present invention have been described, it will be readily appreciated by those skilled in the art that many modifications may be made to the exemplary embodiments without departing from the novel teachings and advantages of the present invention. Therefore, all such modifications are intended to be included within the scope of the present invention as defined in the claims. It should be understood that the above is an explanation of the present invention and should not be considered as being limited to the specific embodiments disclosed, and modifications to the disclosed embodiments and other embodiments are intended to be included within the scope of the appended claims. The present invention is defined by the claims and their equivalents.

Claims

1. A method for prospecting polymetallic minerals, characterized in that: The following steps are involved: S1: Collect fold information, fault information and geophysical information; S2: By analyzing the fold status and fault status of each fold information and fault information, the structural survey value of each area is obtained; S3: Compare the wrinkle status value of each area with the set comparison threshold to obtain the judgment signal of each area; S4: Receive geophysical exploration signals, analyze gravity anomalies, magnetic anomalies and electrical differences in areas with high probability of mineral deposits, and obtain geophysical exploration status values ​​of each geophysical exploration area; S5: Receive the geophysical status values ​​of each geophysical exploration area, and compare and analyze them with the set reference standard interval to obtain different levels of mineral area sets corresponding to different geophysical exploration areas.

2. A polymetallic mineral exploration method according to claim 1, characterized in that: The analysis steps of the structural status values ​​of each area are as follows: Obtain the inclination of each rock layer in the region, perform average processing, and obtain the average inclination of the rock layer in each region. Combined with the fold amplitude value of each region, the complex shape value of the fold fissure of each region, and the connectivity value of the fissure of each region, calculate and obtain the fold shape value of each region; The fold condition survey value of each region and the fault condition survey value of each region are multiplied by the corresponding correction factor coefficient, and then summed up to obtain the structural condition survey value of each region.

3. A polymetallic mineral exploration method according to claim 2, characterized in that: The analysis steps of the complex shape value of the fold cracks in each region and the connectivity value of the cracks in each region are as follows: The folded rock mass is scanned by three-dimensional scanning. According to the grayscale information and the spatial distribution of the point cloud data, the background point cloud is separated from the crack point cloud by cluster analysis. The center line of the fold crack is extracted by interpolation algorithm. The center line of the crack is placed in the three-dimensional coordinate system. Spheres of different radii are used to completely cover the center line of the crack. Starting from radius R1, R1, R1 1 / 2 , R1 1 / 4 , R1 1 / 8 , R1 1 / 16 , ..., R1 1 / 2×N , record the number of spheres gs corresponding to each radius value, establish a two-dimensional coordinate system with logR as the horizontal axis and gs as the vertical axis, and perform linear fitting on the coordinate points in the coordinate system to obtain the slope of the fitting line and take the absolute value as the complex value of the fold cracks. Sum the complex values ​​of each fold crack and divide it by the number of fold cracks to obtain the average value of the complex values ​​of each fold crack and mark it as the complex value of the fold cracks in each region; Put the center lines of each crack obtained in S4 into the three-dimensional coordinate system, and obtain the number of intersections between crack i and crack j in turn, mark them as the node number Kij, obtain the angle θij between the directions of crack i and crack j, and obtain the minimum spatial distance Dij between crack i and crack j, and analyze to obtain the connectivity value lt of the cracks in each area.

4. A polymetallic mineral exploration method according to claim 2, characterized in that: The analysis steps of the fold amplitude values ​​of each region are as follows: The exploration area is divided into several areas. In each area, the position of the core layer of the same rock layer on both wings in the magnetic direction is obtained by using a magnetic positioning and orientation instrument. The angles θy and θe are obtained. θy and θe are the respective strikes of the same rock layer on both wings. The angles between the core layer and the core cylinder axis of the same rock layer on both wings are obtained, which are marked as the inclination angles αy and αe of the same rock layer on both wings, respectively. The core is obtained by vertical drilling. When the rock layer inclination is on the same side as the angle between the core cylinder axis, the inclination angles are δy=90°-αy and δe=90°-αe, respectively. When the rock layer inclination is on the opposite side to the angle between the core cylinder axis, the inclination angles are δy=90°+αy and δe=90°+αe, respectively. The laser radar emits a laser beam in the survey area. After the laser beam interacts with the rock surface, it is received by the receiving device. The cloud data obtained is filtered to remove noise, and the fold feature points of the rock formation are extracted and converted into the engineering coordinate system to obtain the coordinate points of the fold core and the same rock formation position points on both wings. Then, the horizontal distance between the same rock formation position points on one wing and the same rock formation position points on the other wing is obtained by analysis, which are marked as Ly and Le respectively. A three-dimensional rectangular coordinate system is established with the strike direction of the rock stratum as the x-axis, the direction perpendicular to the strike and horizontal as the y-axis, the vertical direction as the z-axis, and the fold core as the origin (0, 0, 0). The position points of the same rock stratum on both wings are determined as the first point coordinates (x1, y1, z1) and the second point coordinates (x2, y2, z2), respectively. The coordinate point values ​​of the first point coordinates and the second point coordinates are analyzed and obtained, and the amplitude value fd of the folded rock stratum is calculated. The amplitude values ​​of the folds at different cross-sectional positions are calculated, and the amplitude values ​​at different cross-sectional positions are summed and divided by the number of cross-sectional positions to obtain the average amplitude value of the folds at different cross-sectional positions, which is marked as the fold amplitude value fz in each region.

5. A polymetallic mineral exploration method according to claim 2, characterized in that: The analysis steps of the fault profile survey of each area are as follows: By observing the mineral crystals of each fault rock under a transmission electron microscope in each region, the total length of the dislocation line within a certain region of the mineral crystal is obtained, and the dislocation density value wm of the fault rock mineral is obtained by analysis. The orientation difference between two adjacent scanning points of the mineral crystal of the fault rock is obtained through an electron microscope. When the orientation difference is greater than or equal to the set orientation difference threshold, a subgrain boundary is formed between the two adjacent scanning points, and the crystal orientation difference between adjacent scanning points inside each subgrain boundary is obtained, marked as each subgrain boundary angle value, and each subgrain boundary angle value is marked as each subgrain boundary angle value. The average value of the grain boundary angle is processed to obtain the average subgrain boundary angle, the number of subgrain boundaries and the average spacing of each adjacent subgrain boundary are obtained, the average subgrain boundary angle is multiplied by the number of subgrain boundaries to obtain the subgrain boundary product value, the subgrain boundary product value is divided by the average spacing of adjacent subgrain boundaries and multiplied by the correction factor coefficient to obtain the subgrain boundary shape value of the fault rock mineral, the fault displacement, the width of the fault fracture zone, the dislocation density value of the fault rock mineral and the subgrain boundary shape value of the fault rock mineral are obtained, and the fault shape exploration value of each area is obtained by analysis.

6. A polymetallic mineral exploration method according to claim 1, characterized in that: The analysis steps of the determination signals of each region are as follows: The structural exploration value of each area is compared and analyzed with the comparison threshold XXT1. When the structural exploration value of each area is greater than or equal to the comparison threshold XXT1, the corresponding area is a high-probability mineral existence area and a geophysical exploration signal is generated.

7. A polymetallic mineral exploration method according to claim 1, characterized in that: The specific analysis of the geophysical conditions of each geophysical exploration area is as follows: The proportion of low resistivity, minimum resistivity, magnetic field intensity amplitude and re-evaluation value are calculated to obtain the geophysical exploration status of each geophysical exploration area.

8. A method for prospecting for polymetallic minerals according to claim 7, characterized in that: The analysis steps of the low resistivity ratio, minimum resistivity, magnetic field intensity amplitude value and re-evaluation value are as follows: A square grid is used in the area with high possibility of mineral existence, in which the grid is set at a certain interval. The measured gravity value of each grid center is measured by a gravimeter. According to the geographical latitude and elevation value of the location of each grid center, the standard gravity value is obtained. The measured gravity value of each grid center is subtracted from the standard gravity value to obtain the gravity anomaly value of each grid center. With the position distance of each grid center as the horizontal axis and the gravity anomaly value as the vertical axis, a curve graph of the gravity anomaly value changing with the position distance is established. The two coordinate points corresponding to the gravity amplitude on both sides of each peak value dropping to half of the peak value are obtained, and the distance between the two coordinate points is taken and marked as the gravity anomaly half-width value. The average gravity anomaly half-width value pb of each peak value is taken, and the number of peaks fg and the average spacing of each peak pj are counted. The peak value of each peak is obtained and compared with the set large mineral peak value to obtain the number of peaks greater than the set large mineral peak value. The value is divided by the total number of peaks to obtain the large peak proportion value df, and the re-evaluation value zpz of each area is obtained by analysis; Electrodes are set at a certain distance in the area with high possibility of mineral existence. The resistivity of each measuring point is obtained by a resistivity meter, and the resistivity is compared with the set resistivity threshold. When the resistivity is less than the set resistivity threshold, the measuring point is defined as a low-resistance abnormal point. The proportion of low-resistance abnormal points in the total number of measurements is counted and marked as the low-resistance proportion dzb. The lowest resistivity zd of the measuring point is obtained. The magnetometer is used to measure on a regular grid in the area with high possibility of mineral existence to obtain the magnetic field strength of each measuring point, and the maximum value of the magnetic field strength of each measuring point is taken and marked as the magnetic field strength amplitude value fq.

9. A polymetallic mineral exploration method according to claim 1, characterized in that: The analysis steps of the different geophysical exploration areas corresponding to the different grades of mineral area sets are as follows: When the geophysical condition value of each geophysical exploration area is within the preset gradient reference standard interval XX0, the corresponding area is classified as the first-level large mineral area set W0; when the geophysical condition value of the geophysical exploration area is within the preset gradient reference standard interval XX1, the corresponding area is classified as the second-level medium-sized mineral area set W1; when the geophysical condition value of the geophysical exploration area is within the preset gradient reference standard interval XX2, the corresponding area is classified as the third-level small mineral area set W2; when the geophysical condition value of each geophysical exploration area is within the preset gradient reference standard interval XX3, the corresponding area is classified as the fourth-level non-ore abnormal area set W3, and the level set of each geophysical exploration area is divided and sent to the system display platform for display.

10. A polymetallic mineral exploration system, characterized in that A polymetallic mineral exploration method as described in any one of claims 1 to 9, comprising: The data acquisition module is used to collect fold information, fault information and geophysical information, and send it to the structural analysis module and the geophysical exploration module; The structural analysis module is used to receive fold information and fault information, analyze the fold state and fault state, obtain the structural state value of each area, and compare it with the set comparison threshold to obtain the judgment signal of each area; The geophysical exploration module is used to receive geophysical exploration signals, analyze gravity anomalies, magnetic anomalies and electrical differences in areas with high probability of mineral deposits, and obtain geophysical exploration status values ​​of each geophysical exploration area; The mineral assessment module is used to receive the geophysical conditions of each geophysical exploration area, and compare and analyze them with the set reference standard interval to obtain a set of mineral areas of different grades corresponding to different geophysical exploration areas.

Images