A polymetallic mineral exploration method and system

Abstract

The present invention discloses a polymetallic mineral exploration method and system, which relates to the field of polymetallic mineral exploration and includes a data acquisition module, a structural analysis module, a geophysical exploration module, and a mineral assessment module. The method includes the following steps: collecting fold information, fault information, and geophysical information; analyzing the fold state and fault state of each fold information and fault information to obtain a structural state value for each region; comparing the fold state value of each region with a set comparison threshold to obtain a determination signal for each region; receiving geophysical exploration signals; analyzing gravity anomalies, magnetic anomalies, and electrical differences in regions with high probability of mineral presence to obtain geophysical state values ​​for each geophysical exploration region; and comparing and analyzing the geophysical state values ​​with a set reference standard interval to obtain a set of mineral regions of different grades corresponding to different geophysical exploration regions. The method can classify polymetallic mineral regions according to mineral reserves, providing a basis for formulating polymetallic mineral exploration and mine development designs.

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. These inefficiencies cannot meet the growing demand for polymetallic minerals, and traditional methods cannot guarantee accuracy. This has led to the emergence of polymetallic mineral exploration methods and systems.

[0003] During the exploration process 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 solved.

[0004] In order to solve the above-mentioned 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 embodiments of the present invention provide a method and system for multi-metallic 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 polymetallic minerals, comprising the following steps:

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

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

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

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

[0012] S5: Receive the geophysical condition values ​​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.

[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. Combined with the fold amplitude value fz of each region, the complex shape value fx of each region's fold fractures, and the connectivity value lt of each region's fractures, calculate the fold shape value zkz of each region.

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

[0016] Furthermore, the analysis steps of the complex shape value of the folds and 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 and the crack point cloud are separated 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, R29, R29, R21, R2 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 within 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] Place the center lines of each crack obtained in S4 into a three-dimensional coordinate system, and obtain the number of intersections between crack i and crack j in turn, marked 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 cylindrical axis of the same rock layer on both wings are obtained and 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 dips to the same side as the angle between the core cylindrical axis, the dip angles are δy = 90° - αy and δe = 90° - αe, respectively. When the rock layer dips to the opposite side of the angle between the core cylindrical axis, the dip 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 resulting cloud data is filtered to remove noise, and the characteristic points of the rock folds are extracted and converted into the engineering coordinate system. The coordinates of the fold core and the same rock stratum positions on both flanks are obtained. The horizontal distances between the same rock stratum positions on one flank and the same rock stratum positions on the other flank are then analyzed and 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, and the vertical direction as the z-axis, with 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 obtained by analysis, and the amplitude value fd of the fold rock stratum is obtained by calculation. The amplitude values ​​of the folds at different cross-sectional positions are calculated, and the amplitude values ​​of the folds 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 of 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 by electron microscopy. 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 of the 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 value is processed to obtain the average subgrain boundary angle value, the number of subgrain boundaries and the average spacing of each adjacent subgrain boundary are obtained, the average subgrain boundary angle value and the number of subgrain boundaries are multiplied 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 fault fracture zone width, 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 determination signals of each region 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 steps for analyzing the geophysical conditions of each geophysical exploration area are as follows:

[0028] The low resistivity ratio dzb, the lowest resistivity zd, the magnetic field intensity amplitude fq and the reassessment value zpz are calculated to obtain the geophysical condition value wkz of each geophysical exploration area.

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

[0030] A square grid is used in areas with high probability of mineral deposits, in which the grids are set at a certain spacing. The measured gravity value of each grid center is measured by a gravimeter. According to the geographical latitude and elevation values ​​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 versus 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 between each peak value 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. The reassessment 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. 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 anomaly point. The proportion of low-resistance anomaly 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. The maximum magnetic field strength of each measuring point is taken and marked as the magnetic field strength amplitude value fq.

[0032] Furthermore, the analysis steps for the different geophysical exploration areas corresponding to the different grades of mineral regions 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 grade 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 polymetallic mineral exploration system, characterized by comprising:

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

[0036] The structural analysis module is used to receive the information of each fold and fault, analyze the fold status and fault status, obtain the structural status 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 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 status and fault status of fold information and fault information to obtain the average rock layer dip angle, fold amplitude value, complex shape value of fold fissures, connectivity value of fissures and fault condition value of each region, and then analyzes and obtains the structural condition 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 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 conduct further 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 conditions 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, providing 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 to the actual size, 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 following is a clear and complete description of the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the embodiments described 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 making creative efforts 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 geophysical exploration module.

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

[0049] The specific steps for solving the fold condition value are as follows:

[0050] Step 1: Divide the exploration area into several areas. In each area, use a magnetic positioning and orientation instrument to obtain the position of the core layer of the same rock layer on both wings in the magnetic direction. The obtained angles θy and θe, θy and θe are the respective strikes of the same rock layer on both wings, and the angles between the core layer and the core cylindrical axis of the same rock layer on both wings are obtained, 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 dips to the same side as the angle between the rock layer and the core cylindrical axis, the dip angles are δy = 90° - αy and δe = 90° - αe, respectively. When the rock layer dips to the opposite side of the angle between the rock layer and the core cylindrical axis, the dip angles are δy = 90° + αy and δe = 90° + αe, respectively.

[0051] Step 2: Use a laser radar to emit a laser beam in the survey area. After the laser beam interacts with the rock surface, it is received by a receiving device. The obtained cloud data is filtered to remove noise, and the characteristic points of the rock fold are extracted and converted into the engineering coordinate system to obtain the coordinates of the fold core and the same rock stratum position points on both wings. The abscissa of the fold core is subtracted from the abscissa of the same rock stratum position point on one flank and squared to obtain the square value of the abscissa. The ordinate of the fold core is subtracted from the ordinate of the same rock stratum position point on one flank and squared to obtain the square value of the ordinate. The square values ​​of the abscissa and ordinate are summed and the square root is taken to obtain the horizontal distance between the points in the same rock stratum on one flank. Repeat the above operation to obtain the horizontal distance between the points in the same rock stratum on the other flank, 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 Obtain the amplitude value fd of the folded rock layer, where λ is the set correction factor coefficient used to improve the accuracy of the calculation results. The specific value is determined by professionals in this field. Obtain the amplitude value of the fold at different cross-sectional positions and calculate its amplitude value. Sum the amplitude values ​​at different cross-sectional positions and divide them by the number of cross-sectional positions to obtain the average amplitude value of the fold at different cross-sectional positions, and mark it 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 and the crack point cloud are separated 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. The center line of the crack is completely covered by spheres of different radii, starting from radius R1, and 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 within 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, and obtain the number of intersections between crack i and crack j in turn, marked 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 Obtain the connectivity value lt of each region's cracks, where to1, to2, to3, to4, to5, to6, to7, and to8 are all set weight factor coefficients used 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 dip angle of each rock layer in the region, perform average processing, and obtain the average dip angle of each 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 each region's fold cracks, and the connectivity value lt of each region's cracks. Substitute it into the set formula zkz=f1×pj+f2×fz+f3×fx+f4×lt to calculate the fold shape value zkz of each region, where f1, f2, f3 and f 4 are the weight factor coefficients for the mean dip angle of each region, the fold amplitude value of each region, the complex shape value of each region's folds and fissures, and the connectivity value of each region's fissures. 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 layer, the greater the fold amplitude value, the greater the complex shape value of the folds and 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 in 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 by electron microscopy. 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 of the 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 subgrain boundary angle value, obtaining the number of subgrain boundaries and the average spacing of each adjacent subgrain boundary, multiplying the mean subgrain boundary angle value and 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 and the width of the fault fracture zone, and performing weighted calculation with the dislocation density value and the subgrain boundary shape value of the fault rock mineral, and multiplying it by the corresponding weight factor coefficient to obtain the fault shape survey 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] Multiply the fold and fault profile values ​​of each region by the correction factor coefficients d1 and d2 respectively, and then sum them up to obtain the structural profile 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-probability mineral deposit 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 deposit 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 geophysical status values ​​for each geophysical exploration area. The specific analysis steps are as follows:

[0061] A square grid is used in areas with high probability of mineral deposits, where the grids are 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 versus position 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 compared with the set large mineral peak value, 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, 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, and the formula set is used. Get the reassessment value zpz of each region, where XOX1, XOX2, XOX3 and XOX4 are the weight factor coefficients of the mean 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. 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 deposits. 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 marked as a low-resistance anomaly point. The proportion of low-resistance anomaly points in the total number of measurements is counted and marked as the low-resistance ratio 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 deposits to obtain the magnetic field strength of each measuring point. The maximum magnetic field strength of each measuring point is taken and marked as the magnetic field strength amplitude value fq.

[0063] Normalize the low resistivity ratio dzb, the lowest resistivity zd, the magnetic field intensity amplitude fq and the reassessment value zpz 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 low resistivity ratio, the lowest resistivity, the magnetic field intensity amplitude value and the reassessment value, respectively, with values ​​of 1.122, 2.215, 3.214 and 3.21, and ρ is the preset correction factor, with a specific value of 0.651. It should be noted that when the resistivity meter is measured, the more metal minerals, the lower the resistivity value; 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 a set of mineral areas of different grades 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. Among them, 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 into 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 into 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 into 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 into the fourth-level no-ore abnormal area set W3, and the grade 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 polymetallic mineral exploration method comprises the following steps:

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

[0070] S2: By analyzing the fold and fault status of each fold 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 potential for mineral deposits, and obtain geophysical status values ​​for each geophysical exploration area;

[0073] S5: Receive the geophysical condition values ​​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.

[0074] The above is an illustration of the present invention and should not be considered as limiting thereof. Although several exemplary embodiments of the present invention have been described, it will be readily understood 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 by the claims. It should be understood that the above is an illustration 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 and fault status of each fold and fault information, the structural survey value of each area is obtained; The analysis steps of the structural survey values ​​of each area are as follows: Obtain the dip angle of each rock layer in the region, perform average processing, and obtain the average dip angle of each rock layer in each region. Combined with the fold amplitude value of each region, the complex shape value of each region's fold fractures, and the connectivity value of each region's fractures, the fold shape value of each region is calculated; The analysis steps of the wrinkle amplitude values ​​in 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, and 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 cylindrical axis of the same rock layer on both wings are obtained, and are marked as the position 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 cylindrical axis, the inclination angles are δy=90°-αy and δe=90°-αe, respectively. When the rock layer inclination is on the opposite side of the angle between the core cylindrical 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 resulting cloud data is filtered to remove noise, and the characteristic points of the rock folds are extracted and converted into the engineering coordinate system. The coordinates of the fold core and the same rock stratum positions on both flanks are obtained. The horizontal distances between the same rock stratum positions on one flank and the same rock stratum positions on the other flank are then analyzed and marked as Ly and Le, respectively. Establish a three-dimensional rectangular coordinate system with the rock formation strike direction 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). Determine the position points of the same rock formation 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 coordinate and the second point coordinate according to the set formula Get the amplitude value fd of the folded rock layer, then repeat the above calculation process to obtain the amplitude value of the fold at different cross-sectional positions, sum the amplitude values ​​at different cross-sectional positions and divide them by the number of cross-sectional positions to get the average amplitude value of the fold at different cross-sectional positions, and mark it as the fold amplitude value fz of each region; Multiply the fold and fault profile values ​​of each region by the corresponding correction factor coefficients, and then sum them up to obtain the structural profile value of each region; S3: Compare the structural 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 potential for mineral deposits, and obtain geophysical status values ​​for each geophysical exploration area; S5: Receive the geophysical condition values ​​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.

2. A polymetallic mineral exploration method according to claim 1, characterized in that: The analysis steps of the complex shape value of the folds and 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 and the crack point cloud are separated 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, and R2 are selected in turn. 1 / 2 , R1 1 / 4 , R1 1 / 8 , R1 1 / 16 , ..., R1 1 / 2×N , N is the maximum value of the crack number, record the number of spheres gs corresponding to each radius value, use logR as the horizontal axis and gs as the vertical axis to establish a two-dimensional coordinate system, and perform linear fitting on the coordinate points in the two-dimensional coordinate system to obtain the slope of the fitting line and take its absolute value as the complex value of the fold crack. 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 crack in each region. The center lines of each fracture extracted by the interpolation algorithm are placed in a three-dimensional coordinate system. The number of intersection points between fracture i and fracture j is obtained in turn, marked as the node number Kij, the angle θij between the directions of fracture i and fracture j is obtained, and the minimum spatial distance Dij between fracture i and fracture j is obtained. The connectivity value lt of the fractures in each region is obtained by analysis.

3. A polymetallic mineral exploration method according to claim 1, characterized in that: The analysis steps of the fault profile survey of each region are as follows: The mineral crystals of the fault rocks in each region were observed by transmission electron microscopy, and the total length of the dislocation lines in a certain region of the mineral crystals was obtained. The dislocation density value wm of the fault rock minerals was obtained by analysis. The orientation difference of two adjacent scanning points of the mineral crystals of the fault rock was obtained by electron microscopy. When the orientation difference was greater than or equal to the set orientation difference threshold, a subgrain boundary was formed between the two adjacent scanning points. The crystal orientation difference of the adjacent scanning points inside each subgrain boundary was obtained and marked as the angle value of each subgrain boundary. The average value of the subgrain 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 and the number of subgrain boundaries are multiplied 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 fault fracture zone width, 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.

4. 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.

5. 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 reassessment value are calculated to obtain the geophysical condition value of each geophysical exploration area.

6. A polymetallic mineral exploration method according to claim 5, characterized in that: The analysis steps for the low resistivity ratio, minimum resistivity, magnetic field intensity amplitude value and reassessment value are as follows: A square grid is used in areas with high probability of mineral deposits, in which the grids are set at a certain spacing. The measured gravity value of each grid center is measured by a gravimeter. According to the geographical latitude and elevation values ​​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 versus 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 between each peak value 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. The reassessment 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 compared with the set resistivity threshold. When the resistivity is less than the set resistivity threshold, the measuring point is determined to be a low-resistance anomaly point. The proportion of low-resistance anomaly 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. The maximum magnetic field strength of each measuring point is taken and marked as the magnetic field strength amplitude value fq.

7. A polymetallic mineral exploration method according to claim 1, characterized in that: The analysis steps for the different geophysical exploration areas corresponding to the different grades of mineral areas 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 grade set division of each geophysical exploration area is sent to the system display platform for display.

8. A polymetallic mineral exploration system, characterized in that A polymetallic mineral exploration method according to any one of claims 1 to 7, comprising: Data acquisition module, used to collect fold information, fault information and geophysical information, and send it to the structural analysis module and geophysical exploration module; The structural analysis module is used to receive fold information and fault information, analyze the fold status and fault status, obtain the structural status 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 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