CN117310826A - Geophysical method for delineating favorable target area of ore formation and prospecting - Google Patents
Geophysical method for delineating favorable target area of ore formation and prospecting Download PDFInfo
- Publication number
- CN117310826A CN117310826A CN202311272615.5A CN202311272615A CN117310826A CN 117310826 A CN117310826 A CN 117310826A CN 202311272615 A CN202311272615 A CN 202311272615A CN 117310826 A CN117310826 A CN 117310826A
- Authority
- CN
- China
- Prior art keywords
- ore
- target area
- geophysical
- prospecting
- aeromagnetic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 64
- 230000002349 favourable effect Effects 0.000 title claims abstract description 33
- 230000015572 biosynthetic process Effects 0.000 title claims description 32
- 238000005259 measurement Methods 0.000 claims abstract description 65
- 230000002159 abnormal effect Effects 0.000 claims abstract description 26
- 238000012795 verification Methods 0.000 claims abstract description 18
- 238000012545 processing Methods 0.000 claims abstract description 12
- 238000005553 drilling Methods 0.000 claims abstract description 11
- 238000011835 investigation Methods 0.000 claims abstract description 9
- 230000010287 polarization Effects 0.000 claims description 17
- 238000012937 correction Methods 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 7
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 6
- 239000011707 mineral Substances 0.000 claims description 6
- 238000005065 mining Methods 0.000 claims description 5
- 230000007480 spreading Effects 0.000 claims description 4
- 238000003892 spreading Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 239000011435 rock Substances 0.000 abstract description 21
- 238000011160 research Methods 0.000 abstract description 11
- 238000001514 detection method Methods 0.000 abstract description 2
- 230000005856 abnormality Effects 0.000 description 20
- 230000005284 excitation Effects 0.000 description 14
- 230000009286 beneficial effect Effects 0.000 description 10
- 230000033558 biomineral tissue development Effects 0.000 description 10
- WFKWXMTUELFFGS-UHFFFAOYSA-N tungsten Chemical compound [W] WFKWXMTUELFFGS-UHFFFAOYSA-N 0.000 description 10
- 229910052721 tungsten Inorganic materials 0.000 description 10
- 239000010937 tungsten Substances 0.000 description 10
- 238000010276 construction Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 4
- 230000008901 benefit Effects 0.000 description 4
- 239000010438 granite Substances 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- JQJCSZOEVBFDKO-UHFFFAOYSA-N lead zinc Chemical compound [Zn].[Pb] JQJCSZOEVBFDKO-UHFFFAOYSA-N 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 102100033538 Clusterin-associated protein 1 Human genes 0.000 description 2
- 240000009087 Crescentia cujete Species 0.000 description 2
- 235000005983 Crescentia cujete Nutrition 0.000 description 2
- 101000945057 Danio rerio Clusterin-associated protein 1 homolog Proteins 0.000 description 2
- 101000945060 Homo sapiens Clusterin-associated protein 1 Proteins 0.000 description 2
- 235000009797 Lagenaria vulgaris Nutrition 0.000 description 2
- ZOKXTWBITQBERF-UHFFFAOYSA-N Molybdenum Chemical compound [Mo] ZOKXTWBITQBERF-UHFFFAOYSA-N 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 238000005520 cutting process Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 229910052742 iron Inorganic materials 0.000 description 2
- SZVJSHCCFOBDDC-UHFFFAOYSA-N iron(II,III) oxide Inorganic materials O=[Fe]O[Fe]O[Fe]=O SZVJSHCCFOBDDC-UHFFFAOYSA-N 0.000 description 2
- 229910052750 molybdenum Inorganic materials 0.000 description 2
- 239000011733 molybdenum Substances 0.000 description 2
- 230000035772 mutation Effects 0.000 description 2
- 229910052720 vanadium Inorganic materials 0.000 description 2
- LEONUFNNVUYDNQ-UHFFFAOYSA-N vanadium atom Chemical compound [V] LEONUFNNVUYDNQ-UHFFFAOYSA-N 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 240000001592 Amaranthus caudatus Species 0.000 description 1
- 235000009328 Amaranthus caudatus Nutrition 0.000 description 1
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- CWYNVVGOOAEACU-UHFFFAOYSA-N Fe2+ Chemical compound [Fe+2] CWYNVVGOOAEACU-UHFFFAOYSA-N 0.000 description 1
- 241000923606 Schistes Species 0.000 description 1
- 235000012735 amaranth Nutrition 0.000 description 1
- 239000004178 amaranth Substances 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 229910001919 chlorite Inorganic materials 0.000 description 1
- 229910052619 chlorite group Inorganic materials 0.000 description 1
- QBWCMBCROVPCKQ-UHFFFAOYSA-N chlorous acid Chemical compound OCl=O QBWCMBCROVPCKQ-UHFFFAOYSA-N 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- -1 gold copper lead zinc chromium iron tungsten rare earth Chemical class 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000005389 magnetism Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 229910052683 pyrite Inorganic materials 0.000 description 1
- 239000011028 pyrite Substances 0.000 description 1
- NIFIFKQPDTWWGU-UHFFFAOYSA-N pyrite Chemical compound [Fe+2].[S-][S-] NIFIFKQPDTWWGU-UHFFFAOYSA-N 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- 229910052761 rare earth metal Inorganic materials 0.000 description 1
- 229910052701 rubidium Inorganic materials 0.000 description 1
- IGLNJRXAVVLDKE-UHFFFAOYSA-N rubidium atom Chemical compound [Rb] IGLNJRXAVVLDKE-UHFFFAOYSA-N 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/15—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for use during transport, e.g. by a person, vehicle or boat
- G01V3/16—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for use during transport, e.g. by a person, vehicle or boat specially adapted for use from aircraft
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Aviation & Aerospace Engineering (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a geophysical method for delineating an ore-forming favorable target area and finding ores, which comprises the following steps: acquiring aeromagnetic measurement data of a first target area; processing the aeromagnetic measurement data to obtain an aeromagnetic measurement Δt contour map; setting an ore-beneficial target area in an inner circle of an aeromagnetic measurement delta T contour map; carrying out ground investigation and geophysical prospecting accurate detection section verification in an ore-forming favorable target area, and delineating an abnormal area of ore; establishing a geological-geophysical model by multi-parameter mutual constraint inversion of the combined geophysical prospecting measurement; drilling verification is carried out to determine the position of the mine to be found, and thus mine finding is completed. The invention can collect a large amount of rich magnetic field data by utilizing a large proportion of low-altitude aeromagnetic, and can divide the breaking structure of a research area by utilizing the data, define basic and super basic rock boundaries and define an ore target area by combining an ore forming theory of the research area.
Description
Technical Field
The invention belongs to the field of ore-forming searching, and particularly relates to a geophysical method for defining an advantageous target area of ore-forming and searching.
Background
At present, conventional means for finding mineralized bodies are geological, geochemical and remote sensing technologies, but the prospecting methods are only suitable for shallow surface coverage areas, the loess coverage area of China exceeds 30 ten thousand square kilometers, loess coverage rock (ore) is less exposed, surface exploration technologies such as geology and geochemistry are difficult to play a role, and the mineral resource exploration degree is low overall. Along with the deep mineral exploration degree of China, how to rapidly and effectively outline the target area of the beneficial area of the mine in the loess coverage area is a main problem at present.
Disclosure of Invention
The invention aims to provide a geophysical method for defining an advantageous target area of ore and finding ores so as to overcome the technical defects.
In order to solve the technical problems, the invention provides a geophysical method for delineating an advantageous target area of a mine and finding the mine, which comprises the following steps:
acquiring aeromagnetic measurement data of a first target area;
processing the aeromagnetic measurement data to obtain an aeromagnetic measurement Δt contour map;
setting an ore-beneficial target area in an inner circle of an aeromagnetic measurement delta T contour map;
carrying out ground investigation and geophysical prospecting accurate detection section verification in an ore-forming favorable target area, and delineating an abnormal area of ore;
establishing a geological-geophysical model by multi-parameter mutual constraint inversion of the combined geophysical prospecting measurement;
drilling verification is carried out to determine the position of the mine to be found, and thus mine finding is completed.
Acquiring aeromagnetic measurement data of a first target region, comprising:
taking a mining remote scenic spot defined based on Chinese basic geological survey results as a first target area;
carrying out low-altitude aeromagnetic measurement on a first target area, and flying along the relief of the terrain, wherein the flying height is within 160 m;
and recording aeromagnetic measurement data, wherein the aeromagnetic measurement data at least comprises a total magnetic field T in the horizontal gradient direction.
Processing the aeromagnetic measurement data to obtain an aeromagnetic measurement Δt contour map, comprising:
performing daily change correction, altitude correction and gradient correction on the total magnetic field T in the horizontal gradient direction to obtain a magnetic anomaly delta T before non-leveling;
and (3) carrying out trend approximation leveling and low-pass filtering treatment on the magnetic anomaly DeltaT before non-leveling, finally obtaining DeltaT magnetic anomaly data, meshing the DeltaT magnetic anomaly data, generating GRD grid data, and drawing an aeromagnetic measurement DeltaT contour map.
An ore-bearing beneficial target zone is defined in an aeromagnetic survey DeltaT contour map, comprising:
and (3) sequentially carrying out magnetic pole formation, upward extension and vertical second derivative comprehensive treatment on the delta T magnetic anomaly data, then encircling the boundary of a fixed magnetic body in an aeromagnetic measurement delta T contour map by taking a vertical second derivative 0 value line as a boundary, and combining with an ore formation theoretical model of the area to obtain an ore formation favorable target area.
If the number of the ore-forming favorable target areas is at least 2, selecting an ore-forming potential target area, wherein the selecting steps are as follows:
after carrying out magnetic pole processing on delta T magnetic anomaly data, defining a presumed fracture zone in a delta T contour map of the magnetic measurement according to the judgment basis of the magnetic fracture, or carrying out magnetic pole precedent on delta T magnetic anomaly data and then defining a presumed fracture zone in the delta T contour map of the magnetic measurement by using horizontal first derivative strip anomaly;
the ore-forming favorable target area on the presumed fracture zone is taken as an ore-forming potential target area.
Performing ground investigation and geophysical prospecting profile verification in an ore-forming beneficial target area, and delineating an abnormal region of ore, comprising: and (3) checking an ore-forming favorable target area or an ore-forming potential target area by adopting a comprehensive geophysical prospecting method, determining the abnormal characteristics of the ore according to the found ore body, wherein the area with the abnormal characteristics of the ore is an abnormal area of the ore.
Multi-parameter, mutually constrained inversion of joint geophysical prospecting measurements to build a geologic-geophysical model comprising: the comprehensive geophysical prospecting method comprises a magnetic method, a gravitational method, a natural potential method, a charging method, a frequency spectrum induced polarization method, an induced polarization center ladder method and an induced polarization sounding method, at least two comprehensive geophysical prospecting methods are selected to form a comprehensive geophysical refined survey section, and a geological-geophysical model is established by multi-parameter mutual constraint inversion related to the comprehensive geophysical refined survey section.
Drilling verification to determine the location of the prospected body, completing prospecting, comprising: according to the geological-geophysical model, the tendency, the dip angle, the scale and the spatial spreading characteristics of mineralized bodies are ascertained, and drilling verification is scientifically arranged.
The beneficial effects of the invention are as follows:
aiming at the areas where the geological and chemical exploration of the loess coverage area is difficult to develop, the invention can acquire a large amount of rich magnetic field data by utilizing a large proportion of low-altitude aeromagnets, can divide the fracture structure of a research area by utilizing the data, and can define the boundaries of basic and super basic rock bodies and define an ore target area by combining the ore theory of the research area. And carrying out ground electric method, magnetic method and gravity in the target area for key verification, adopting multi-method joint constraint inversion, and finally establishing a true and reliable geological-geophysical model, scientifically and effectively guiding drilling arrangement, and achieving the purpose of fine prospecting.
In order to make the above-mentioned objects of the present invention more comprehensible, preferred embodiments accompanied with figures are described in detail below.
Drawings
FIG. 1 is a flow chart of a geophysical method for delineating an advantageous target area for mine production.
Fig. 2 is a contour plot of the aeromagnetic measurement Δt.
Fig. 3 is a graph of the post-abnormal extension 3D surface of the aeromagnetic survey Δt.
Fig. 4 is a diagram of a presumed fracture structure and a diagram of an ore-forming advantageous target area b.
FIG. 5 is a plot of the section of the ladder 2 line in an excitation and a quantitative inversion chart.
FIG. 6 is a ladder-in-excitation measurement ρ s 、η s A graph;
fig. 7 is a 7-line comprehensive cross-sectional view of a number 1 mine-forming beneficial target area.
Detailed Description
Further advantages and effects of the present invention will become apparent to those skilled in the art from the disclosure of the present specification, by describing the embodiments of the present invention with specific examples.
In the present invention, the upper, lower, left and right directions in the drawings are regarded as the upper, lower, left and right directions of the geophysical method for defining the favorable target area for ore and prospecting described in the present specification.
The exemplary embodiments of the present invention will now be described with reference to the accompanying drawings, however, the present invention may be embodied in many different forms and is not limited to the examples described herein, which are provided to fully and completely disclose the present invention and fully convey the scope of the invention to those skilled in the art. The terminology used in the exemplary embodiments illustrated in the accompanying drawings is not intended to be limiting of the invention. In the drawings, like elements/components are referred to by like reference numerals.
Unless otherwise indicated, terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art. In addition, it will be understood that terms defined in commonly used dictionaries should be interpreted as having a meaning that is consistent with their meaning in the context of the relevant art and will not be interpreted in an idealized or overly formal sense.
The following magnetic anomalies are aeromagnetic anomalies.
The present embodiments relate to a geophysical method for delineating an ore-forming beneficial target area and for prospecting, comprising:
step 100, obtaining aeromagnetic measurement data of a first target region.
The implementation mode adopts a power glider aeromagnetic measurement system to carry out low-altitude aeromagnetic measurement work, and the specific steps include:
step 101, taking a prospecting remote scenic spot defined based on Chinese basic geological survey results as a first target area.
The Chinese basic geological survey refers to the completed research results obtained by 1:10 ten thousand, 1:5 ten thousand regional geological survey and mineral geological survey, the regional ore-forming geological conditions are combed on the basis of the Chinese basic geological survey results, and the regional ore-forming geological conditions are preferably 50-100 km in an ore collection area or an ore-forming remote scenic area 2 The investigation area is used as a prospecting unit.
Here, the first target area satisfies the above conditions, and also satisfies the areas where rock (ore) is less exposed, such as the loess coverage area or the desert area.
102, carrying out low-altitude aeromagnetic measurement on a first target area, flying along the relief of the terrain, wherein the flying height is within 160 m.
The low-altitude aeromagnetic measurement is carried out, the working scale is selected to be 1:5000-1:10000, the measuring line distance is 50-100 meters, the cutting line distance is 1000-2000 meters, the point distance is 1-2 meters, the flying ground clearance is controlled to be 30-150 meters, and the measuring line direction is perpendicular to the trend of main stratum in the area.
Step 103, recording aeromagnetic measurement data, wherein the aeromagnetic measurement data at least comprises a total magnetic field T in a horizontal gradient direction.
During low-altitude aeromagnetic measurement, the direction of the survey line is perpendicular to the trend of the main stratum in the first target area, the total magnetic field in the first target area is measured, and finally the total magnetic field T in the horizontal gradient direction is obtained.
Step 200, processing the aeromagnetic measurement data to obtain an aeromagnetic measurement Δt contour map, specifically as follows:
performing daily change correction, altitude correction and gradient correction on a total magnetic field T in a horizontal gradient direction to obtain a magnetic anomaly delta T before non-leveling, wherein most of aeromagnetic anomalies have strip-shaped false anomalies due to system noise, steering difference and the like, performing trend approximation leveling (trend leveling for short) and low-pass filtering treatment on the magnetic anomaly delta T before non-leveling to finally obtain delta T magnetic anomaly data, and performing gridding on the processed total magnetic field in surfer8.0 software by using a Kriging method to generate GRD grid data to obtain a GRD file, and then drawing an aeromagnetic measurement delta T contour map.
It should be noted that the surfer8.0 software is used to grid Δt magnetic anomaly data, so that other software may be used to implement data meshing, and is not limited to surfer8.0 software.
In this embodiment, a plan view of a Δt contour line of 1:1 magneto-aero measurement is drawn (see fig. 2), to further illustrate the advantages of the present magneto-aero measurement, the 1:1 magneto-aero measurement result is compared with the previous 1:5 magneto-aero measurement result, the present 1:1 magneto-aero Δt contour line is basically consistent with the 1:5 magneto-aero Δt contour line in large form, but the present magneto-aero measurement is more detailed in terms of magnetic field details, the abnormal form is more complete, and part of the abnormal form is subdivided into a plurality of peak abnormal components, so that the internal characteristics of the underground magnetic body are more clearly reflected.
Through this time aeromagnetic measurement, provide more magnetic field information for geological prospecting.
Step 300, defining an ore-forming favorable target area in the aeromagnetic measurement DeltaT contour map, specifically, performing structural inference based on a specific magnetic field form, and simultaneously defining the ore-forming favorable target area based on an in-area ore-forming theory, wherein the method comprises the following steps of:
and (3) sequentially carrying out magnetic pole formation, upward extension and vertical second derivative comprehensive treatment on the delta T magnetic anomaly data, then encircling the boundary of a fixed magnetic body in an aeromagnetic measurement delta T contour map by taking a vertical second derivative 0 value line as a boundary, and combining with an ore formation theoretical model of the area to obtain an ore formation favorable target area.
And (3) carrying out pole formation processing (pole formation parameters are survey area magnetic declination and magnetic inclination angle), carrying out upward extension processing on the basis of pole formation magnetic anomaly delta T, wherein extension height can be set between 50 and 2000 meters, referring to FIG. 3, FIG. 3 is a post-pole formation 3D surface map of the pole formation magnetic measurement delta T, which is an original contour map, an upward extension 100mm contour map, an upward extension 200mm contour map and an upward extension 500mm contour map in sequence from bottom to top, then carrying out vertical second derivative processing on each extended magnetic anomaly delta T, and demarcating the boundary of a main magnetic body at different deep parts by taking a vertical second derivative 0 value line as a boundary, and demarcating an ore beneficial target area according to an ore formation theory of a research area.
If the number of the ore-forming favorable target areas is at least 2, selecting an ore-forming potential target area, wherein the selecting steps are as follows:
after carrying out magnetic pole processing on delta T magnetic anomaly data, defining a presumed fracture zone in a delta T contour map of the magnetic measurement according to the judgment basis of the magnetic fracture, or carrying out magnetic pole precedent on delta T magnetic anomaly data and then defining a presumed fracture zone in the delta T contour map of the magnetic measurement by using horizontal first derivative strip anomaly; the ore-forming favorable target area on the presumed fracture zone is taken as an ore-forming potential target area.
Specifically, because the construction movement of the position of the research area is frequent, and the area is provided with more folds and cracks under the construction action, the construction estimation can be performed by utilizing a specific magnetic field form, and common distinguishing marks are as follows: gradient bands, different magnetic field boundaries, linear anomaly bands, anomaly mutation bands, beaded anomaly bands, anomaly stagger bands, wild goose-shaped anomaly bands, radial anomaly bands, etc.
The delineation of an ore-bearing beneficial target zone within the aeromagnetic survey Δt contour plot is described in detail below:
referring to fig. 4, in this aeromagnetic process, by performing upward continuation and vertical second derivative treatment on Δt magnetic anomalies of the magnetic poles, four fracture zones are defined in total in a region according to aeromagnetic fracture discrimination, and the positions and directions of the fracture zones are shown in fig. 4a (presumed fracture structure diagram), wherein the F2 fracture zone is near a large pier village and has been verified to be completely coincident, and tungsten ore points are found at the fracture zones, and ore bodies are commonly controlled by strata and fractures.
Performing pole formation treatment on the delta T magnetic anomaly, performing extension treatment of different heights upwards on the basis of the pole formation, performing vertical second derivative treatment on the extended magnetic anomaly respectively, delineating a granite boundary by taking a 0-value line as a boundary, and referring to the existing research results, wherein the ore formation in the area is closely related to the rock mass and is mainly in a concave space part of the rock mass.
In addition, the northwest fracture zone has the functions of guiding and controlling ore formation, so the characteristic of rock mass, fracture zone and ground checking condition are synthesized for the time, the boundary shape of the rock mass which is shaped like a calabash and is shaped like a neck is defined as an ore formation favorable zone, the regions are defined as favorable target regions of ore formation at 3 positions, see figure 4b, and are respectively numbered 1, 2 and 3, namely 3 positions in the circle in figure 4b, wherein tungsten mineralization clues are found in the checking process of the target region 1, and the later stage is subjected to key checking.
In other words, target 1 is located on or contains the fracture zone, and also belongs to the ore-forming favorable target delimited by the vertical second derivative 0 value line, and is therefore determined to be the ore-forming potential target.
It should be noted that, the target area No. 3 may also be in line with the target area which is located on the fracture zone and belongs to the delineated mine-forming favorable target area, and whether the target area No. 3 is the mine-forming potential target area may be verified in the same manner, which is not verified one by one in the embodiment.
Step 400, performing ground investigation and geophysical prospecting profile verification in an ore-forming beneficial target area, and delineating an abnormal region of ore, including: and (3) checking an ore-forming favorable target area or an ore-forming potential target area by adopting a comprehensive geophysical prospecting method, determining the abnormal characteristics of the ore according to the found ore body, wherein the area with the abnormal characteristics of the ore is an abnormal area of the ore.
The method comprises the steps of searching various mineralization clues related to ore formation through ground investigation on the well-defined ore formation target areas (No. 1, no. 2 and No. 3 target areas in fig. 4 b) or the ore formation potential target area (No. 1 target area in fig. 4 b), screening out the ore formation target area with large ore formation potential, for example, arranging 1:5000-1:10000 in-plane area measurement, controlling the measuring line distance to be 50-100 meters, controlling the point distance to be 10-40 meters, enabling the measuring line direction to be perpendicular to the stratum trend, and defining the plane range and trend characteristics of a mineralization body.
For a specific manner, reference can be made to example 1.
Step 500, establishing a geological-geophysical model by multi-parameter mutual constraint inversion of combined geophysical prospecting measurement.
The comprehensive geophysical prospecting method comprises a magnetic method, a gravitational method, a natural potential method, a charging method, a frequency spectrum induced polarization method, an induced polarization center ladder method and an induced polarization sounding method, at least two comprehensive geophysical prospecting methods are selected to form a comprehensive geophysical refined survey section, and a geological-geophysical model is established by multi-parameter mutual constraint inversion related to the comprehensive geophysical refined survey section.
Step 600, drilling verification to determine the location of the retrieved ore body, and completing the retrieving.
According to the geological-geophysical model, the tendency, the dip angle, the scale and the spatial spreading characteristics of mineralized bodies are ascertained, and drilling verification is scientifically arranged.
Example 1:
the research area is located in a Huajialian area of the fixed western city of the loess covered area, and the area is basically formed by a large river basin, a long and slow undulating mountain and loess beam and a long and thin gully, and loess is widely distributed. The earth structure belongs to the arc of the Saint Jobi magma in the eastern section of the Qilin mountain area and the arc of the island in the south of the corridor. The ore-forming belt belongs to an Algold-Qilin ore-forming belt of North Qilin gold copper lead zinc chromium iron tungsten rare earth. The ore forming geological conditions of lead zinc, tungsten, gold and copper in the area are good, and the flood palm lead zinc ore, high-family gorge vanadium titano-magnetite, high-family gorge molybdenum ore and the like are found on an ore forming belt, and the ore deposit type is mainly composed of sea-phase volcanic lead zinc ore, magma rock heteromorphic vanadium titano-magnetite, porphyry molybdenum ore and the like.
The area is mostly covered by fourth loess, only a small amount of the grains are exposed, and most of the rocks are non-magnetic or weakly magnetic, wherein only a small amount of the rocks are strong in magnetism, mainly medium-base rocks and old metamorphic rocks, such as maculous secondary granite, amphiclastic rock, amaranth secondary granite, and the magnetic susceptibility is in the range of (11313-55262) x 10 －6 Between SI; the magnetic property of sedimentary rock is very weak or non-magnetic, and the magnetic susceptibility is 625 multiplied by 10 －6S I is within.
The difference of the electrical characteristics of the rock is obvious, the pyrite-based breccia is a low-resistance high-polarization characteristic, the limonite-based chlorite-based sericite is a low-resistance medium-polarization characteristic, the tungsten ore is a low-resistance high-polarization characteristic, the other rock has lower polarization rate, and the difference of the electric resistance is larger.
The low-resistance high-polarization abnormality is used as a prospecting mark, namely the mineral abnormality is characterized by low-resistance high-polarization abnormality, then the low-resistance high-polarization abnormality is further divided into three types according to ground investigation, the first type 1 is an abnormality with obvious mineralization bodies, the second type 1 is an abnormality with yellow iron and brown iron mineralization, the second type 2 is an abnormality with obvious excitation abnormality, but after covering, the property abnormality cannot be judged, the third type abnormality is an excitation abnormality, but the abnormality can be determined to be caused by non-mineral abnormality.
100km of working area investigation area 2 Firstly, a large scale low-altitude aeromagnetic measurement is carried out by adopting the method flow of fig. 1, after daily change correction, altitude correction, gradient correction and cutting line leveling treatment are carried out on collected data, the processed total magnetic field is gridded by adopting a Kriging method in surfer8.0 software to obtain a grd file, and fig. 2 is drawn.
And carrying out upward continuation treatment on the delta T magnetic anomalies, wherein the maximum values of the magnetic anomalies after continuation are shown as follows:
as can be seen from the above table, the decay of C-74-14, C-18-4 and C-18-2, C-83-119, C-18-5 and C-18-6 is faster in the upward extension process, and the decay of C-18-2, C-18-6 is slower, which shows that the field sources causing the magnetic anomalies of C-74-14 and C-18-2 are shallower, and the field sources causing the magnetic anomalies of C-18-2, C-83-119, C-18-5 and C-18-6 are deeper and are homologous magnetic bodies.
Because the construction movement of the position of the research area is frequent, and the area is provided with more folds and cracks under the construction action, the construction estimation can be carried out by utilizing the specific magnetic field form, and the common distinguishing marks are as follows: gradient bands, different magnetic field boundaries, linear anomaly bands, anomaly mutation bands, beaded anomaly bands, anomaly stagger bands, wild goose-shaped anomaly bands, radial anomaly bands, etc. In this aeromagnetic process, by performing upward continuation and vertical second derivative treatment on the delta T (pole) magnetic anomaly, four fracture zones, namely F1, F2, F3 and F4, are defined in the area according to the aeromagnetic fracture discrimination basis, and refer to fig. 4a, wherein the F2 fracture zone is verified to exist near the great pier village, the position and the trend are completely coincident, and a tungsten mine point is found at the position, and the ore body is controlled by the stratum and the fracture.
Performing chemical pole treatment on the aeromagnetic delta T contour line in the research area, performing upward extension on the basis of the chemical pole, further performing vertical second derivative treatment, delineating a granite boundary by taking a 0 value line as a boundary, and referring to the existing ore forming theoretical system, wherein ore forming in the area is closely related to a rock body and is mainly in a concave space part of the rock body. In addition, the northwest fracture zone has the functions of guiding and controlling ore formation, so the characteristic of rock mass, fracture zone and ground checking condition are synthesized for the time, the boundary shape of the rock mass which is shaped like a calabash and is shaped like a neck is defined as an ore formation favorable zone, the regions are defined as favorable target regions for ore formation at 3 positions, and the favorable target regions are numbered as No. 1, no. 2 and No. 3 respectively, wherein tungsten mineralization clues are found in the target region No. 1 in the checking process, key checking is carried out in the later stage, and the inferred graph of the ore formation target region in the region is shown in fig. 4b.
For the ore-forming target area defined in this time, a target area 1 is preferred, 4 ladder sections in the excitation are arranged according to the method flow shown in fig. 1, and referring to fig. 4b, the line numbers are 2 lines, 4 lines, 6 lines and 7 lines respectively, the section length is 1200 meters, and the section direction is 110 degrees. According to the characteristic of the abnormal excitation, 1 low-resistance high-polarization abnormal zone in the north east direction is defined, the number is JD1, the abnormal excitation curve is shown in figure 6, the change of 4 section polarization values is 0.11-6.4%, the change of apparent resistivity is 8.6-1399Ω & OMEGA M, the abnormal excitation is obvious in the figure, and the background value is stable.
FIG. 6 is a ladder-in-excitation measurement ρ s 、η s Graph with ordinate representing apparent resistivity ρ s And the visual polarization rate eta s The abscissa represents the point number, where ρ s Is apparent resistivity in m.m.m.eta.m s Is the visual polarizability, and FIG. 7 is a 7-line comprehensive cross-sectional view of a No. 1 mining beneficial target area, specifically defined as follows:
the tungsten mineralized green mud quartz schist is collected at the excitation abnormality position near the No. 55 line 4, the polarization value is measured to be 5.73%, the resistivity is 155 omega M, the excitation abnormality can be caused, and the excitation abnormality is highly matched with the tungsten mineralization point position, so that JD1 is presumed to be a mining abnormality. The 4 profile excitation abnormal curves are similar in morphology and are all characterized by low-resistance high-polarization abnormality, stable in northeast trend, and a tungsten mineralization zone JD1 is defined by taking the polarization rate of 2.3% as an abnormal boundary, the mineralization bandwidth is 150 meters, the length is 800 meters, the trend is 50 degrees in the northeast, and the mineralization body tendency to overstock can be judged according to the characteristic of slow east of the low-resistance high-polarization curve.
The 2-line profile curve has regular shape and smooth curve, and the 2-line profile curve is semi-quantitatively explained by using approximate formulas h (about 0.5 q) and h (about 0.6 m), wherein h is the top buried depth of the pulse-like polarizer, q is the chord of the half extremum point of the profile curve, and m is the chord tangent distance passing through the inflection point of the eta curve. Fig. 5 is a graph of a section of a line of a ladder 2 in an excitation and a quantitative inversion chart, and as can be seen from fig. 5, q=130 meters, the top buried depth h of the excitation is calculated to be about 65 meters, m=120, and the top buried depth h of the excitation is calculated to be about 72 meters.
In order to highlight the advantage of the method, on the marked abnormal zone of the mineralized body JD1, please refer to fig. 7, the electric sounding work is arranged between the 45-69 # points of the 7-line electric-excitation middle ladder section which is completely covered by the earth surface, the electric characteristics of the sections of the 45-69 # points are obtained through constraint joint inversion, the space spreading characteristics of the mineralized body in the deep part are further ascertained, the polarization rate of the electric sounding is shown in fig. 7, the polarization rate value change is 0.3-4.5%, the electric sounding abnormality is obvious high polarization abnormality in the deep part, a certain extension is arranged in the deep part, the mineralized body range is defined by taking 2.3% as a boundary, the width between the 50-65 # points of the earth surface is about 150 meters, and the underground burial depth is about 370 meters between-150-520 meters.
According to the geological-geophysical model established by the above induced-polarization middle ladder measurement and induced-polarization sounding measurement, drilling verification is carried out at the point of the No. 7 line 51, and 3 layers of tungsten ore bodies and 1 layer of rubidium ore bodies are found, so that a good ore finding effect is achieved.
The invention adopts the low-altitude aeromagnetic rapid sweeping, the abnormal is defined, and then the geophysical prospecting comprehensive method for the electric profile verification is high in prospecting efficiency, low in cost and obvious in effect in the loess coverage area, has obvious advantages, and can provide references for other similar loess coverage areas.
It will be understood by those of ordinary skill in the art that the foregoing embodiments are specific examples of carrying out the invention and that various changes in form and details may be made therein without departing from the spirit and scope of the invention.
Claims (8)
1. A geophysical method for delineating an advantageous target area of a mine and for prospecting comprising:
acquiring aeromagnetic measurement data of a first target area;
processing the aeromagnetic measurement data to obtain an aeromagnetic measurement Δt contour map;
setting an ore-beneficial target area in the inner circle of the aeromagnetic measurement delta T contour map;
performing ground investigation and geophysical prospecting profile verification in the ore-forming favorable target area, and delineating an abnormal area of ore;
establishing a geological-geophysical model by multi-parameter mutual constraint inversion of the combined geophysical prospecting measurement;
drilling verification is carried out to determine the position of the mine to be found, and thus mine finding is completed.
2. The geophysical method for delineating an ore-bearing favorable target area and for prospecting as set forth in claim 1 wherein obtaining aeromagnetic measurement data for the first target area comprises:
taking a mining remote scenic spot defined based on Chinese basic geological survey results as a first target area;
performing low-altitude aeromagnetic measurement on the first target area, and flying along the relief of the terrain, wherein the flying height is within 160 m;
and recording aeromagnetic measurement data, wherein the aeromagnetic measurement data at least comprises a total magnetic field T in the horizontal gradient direction.
3. A geophysical method for delineating an ore-bearing target area and for prospecting as claimed in claim 2, wherein processing the aeromagnetic survey data to obtain an aeromagnetic survey Δt contour map comprises:
performing daily change correction, altitude correction and gradient correction on the total magnetic field T in the horizontal gradient direction to obtain a magnetic anomaly delta T before non-leveling;
and carrying out trend approximation leveling and low-pass filtering treatment on the magnetic anomaly DeltaT before non-leveling to finally obtain DeltaT magnetic anomaly data, gridding the DeltaT magnetic anomaly data to generate GRD grid data, and drawing an aeromagnetic measurement DeltaT contour map.
4. A geophysical method for delineating and locating an ore-bearing target as set forth in claim 3 wherein delineating an ore-bearing target within the aeromagnetic survey Δt contour map comprises:
and sequentially carrying out magnetic pole formation, upward continuation and vertical second derivative comprehensive treatment on the delta T magnetic anomaly data, then, defining the boundary of a magnetic body in the aeromagnetic measurement delta T contour map by taking a vertical second derivative 0 value line as a boundary, and combining an ore formation theoretical model of the area to obtain the ore formation favorable target area.
5. The geophysical method for delineating and prospecting ore forming favorable targets as recited in claim 4, wherein if there are at least 2 ore forming favorable targets, then selecting an ore forming potential target is performed as follows:
after carrying out magnetic pole processing on the delta T magnetic anomaly data, defining a presumed fracture zone in the delta T contour map of the magnetic measurement according to the judgment of the magnetic fracture, or carrying out magnetic pole precedent on the delta T magnetic anomaly data and then defining a presumed fracture zone in the delta T contour map of the magnetic measurement according to horizontal first derivative strip-shaped anomalies;
the ore-forming favorable target area on the presumptive fracture zone is taken as an ore-forming potential target area.
6. The geophysical method for delineating an ore-forming favorable target area and for prospecting as set forth in claim 5, wherein performing a ground survey and a geophysical prospecting profile verification within the ore-forming favorable target area and delineating an abnormal region of mineral matter comprises: and (3) checking the ore-forming favorable target area or the ore-forming potential target area by adopting a comprehensive geophysical prospecting method, and determining the abnormal characteristics of the ore according to the found ore body, wherein the area with the abnormal characteristics of the ore is an abnormal area of the ore.
7. The geophysical method for delineating an ore-bearing target area and for prospecting as set forth in claim 6 wherein the multiparameter, mutually constrained inversion of the joint geophysical prospecting measurements creates a geologic-geophysical model comprising: the comprehensive geophysical prospecting method comprises a magnetic method, a gravitational method, a natural potential method, a charging method, a frequency spectrum induced polarization method, an induced polarization center ladder method and an induced polarization depth measuring method, at least two comprehensive geophysical prospecting methods are selected to form a comprehensive geophysical refined survey section, and a geological-geophysical model is established by combining multi-parameter mutual constraint inversion related to the comprehensive geophysical refined survey section.
8. Geophysical method for delineating a favorable target area for mining and prospecting according to claim 1 or 7, wherein the drilling verification to determine the location of the prospected body, the prospecting is completed, comprises: according to the geological-geophysical model, the tendency, the dip angle, the scale and the spatial spreading characteristics of mineralized bodies are ascertained, and drilling verification is scientifically arranged.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311272615.5A CN117310826A (en) | 2023-09-28 | 2023-09-28 | Geophysical method for delineating favorable target area of ore formation and prospecting |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311272615.5A CN117310826A (en) | 2023-09-28 | 2023-09-28 | Geophysical method for delineating favorable target area of ore formation and prospecting |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117310826A true CN117310826A (en) | 2023-12-29 |
Family
ID=89242091
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311272615.5A Pending CN117310826A (en) | 2023-09-28 | 2023-09-28 | Geophysical method for delineating favorable target area of ore formation and prospecting |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117310826A (en) |
-
2023
- 2023-09-28 CN CN202311272615.5A patent/CN117310826A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109828316B (en) | Calcium-bonded rock type uranium ore prospecting method | |
| CN111045114B (en) | Method for identifying and positioning favorable sand bodies of basalt coverage area sandstone-type uranium deposit mineralization | |
| CN112965141B (en) | Delineating method of ore-forming favorable section of uranium polymetallic ore | |
| CN105510993A (en) | Foreland basin deep buried and compressed type complex gypsum-salt rock identification and distribution prediction method | |
| CN112711078B (en) | Identification method for sedimentary basin deep favorable sandstone type uranium ore formation sand body | |
| CN111045091B (en) | Recognition and positioning method for ancient river channel covered by basalt | |
| CN106772676A (en) | A kind of method of searching heavy-cover area skarn type mineral metal products | |
| CN111967631A (en) | Method for predicting prospect area of early mineralization of sandstone-type uranium deposit exploration | |
| Boubaya | Combining resistivity and aeromagnetic geophysical surveys for groundwater exploration in the Maghnia Plain of Algeria | |
| Rabeh et al. | Structural control of hydrogeological aquifers in the Bahariya Oasis, Western Desert, Egypt | |
| Abdel-Raouf et al. | Implementation of magnetic, gravity and resistivity data in identifying groundwater occurrences in El Qaa Plain area, Southern Sinai, Egypt | |
| Yan et al. | A 3D Geological Model Constrained by Gravity and Magnetic Inversion and its Exploration Implications for the World‐class Zhuxi Tungsten Deposit, South China | |
| CN106772677A (en) | A kind of method for finding area of coverage skarn type Eastern Pamir | |
| CN111045111A (en) | Method suitable for recognizing comprehensive geophysical target area of ground-leaching sandstone-type uranium-bearing basin | |
| CN114814978B (en) | Granite area tungsten-tin ore exploration method based on multiple depth scales | |
| Okwoli et al. | Ground magnetic and electrical resistivity mapping for basement structures over charnokitic terrain in Ado-Ekiti area, Southwestern Nigeria | |
| Oyedele et al. | Groundwater potential evaluation using surface geophysics at Oru-Imope, South-Western Nigeria | |
| CN117310826A (en) | Geophysical method for delineating favorable target area of ore formation and prospecting | |
| Yusoh et al. | Optimization of Archeological anomalies using GIS method for magnetic and resistivity study at Sungai Batu, Lembah Bujang, Kedah (Malaysia) | |
| Anbazhagan et al. | Electrical resistivity survey and shallow subsurface geological study in hard rock terrain, Southern India | |
| Pospeeva et al. | Deep structure of the junction zones of the Chuya tectonic depression and its mountainous frame in Gorny Altai according to magnetotelluric studies | |
| Mendoza et al. | GPR and Magnetic Techniques to Locate Ancient Mining Galleries (Linares, Southeast Spain) | |
| Majzoub | Characterization and delineation of karst geohazards along RM652 using electrical resistivity tomography, Culberson County, Texas | |
| Karriqi et al. | Combination of resistivity “real section” with quantitative interpretation of vertical electrical soundings | |
| Mathieson et al. | The National Museums of Scotland Saqqara survey project, earth sciences 1990–1998 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |