CN113075749A - Method and system for locating favorable ore-forming space of sandstone-type uranium ore in anticline development area - Google Patents
Method and system for locating favorable ore-forming space of sandstone-type uranium ore in anticline development area Download PDFInfo
- Publication number
- CN113075749A CN113075749A CN202110354196.4A CN202110354196A CN113075749A CN 113075749 A CN113075749 A CN 113075749A CN 202110354196 A CN202110354196 A CN 202110354196A CN 113075749 A CN113075749 A CN 113075749A
- Authority
- CN
- China
- Prior art keywords
- favorable
- anticline
- area
- space
- ore
- 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.)
- Granted
Links
- 230000002349 favourable effect Effects 0.000 title claims abstract description 148
- 238000011161 development Methods 0.000 title claims abstract description 45
- 229910052770 Uranium Inorganic materials 0.000 title claims abstract description 42
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 title claims abstract description 42
- 238000000034 method Methods 0.000 title claims abstract description 37
- 239000004576 sand Substances 0.000 claims abstract description 50
- 230000018109 developmental process Effects 0.000 claims abstract description 44
- 230000033558 biomineral tissue development Effects 0.000 claims abstract description 39
- 230000007704 transition Effects 0.000 claims abstract description 29
- 230000007480 spreading Effects 0.000 claims abstract description 12
- 230000000630 rising effect Effects 0.000 claims abstract description 5
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 45
- 239000011707 mineral Substances 0.000 claims description 45
- 229910052704 radon Inorganic materials 0.000 claims description 44
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 claims description 44
- 230000002159 abnormal effect Effects 0.000 claims description 33
- 238000005259 measurement Methods 0.000 claims description 24
- 239000002689 soil Substances 0.000 claims description 24
- 238000010586 diagram Methods 0.000 claims description 13
- 230000005856 abnormality Effects 0.000 claims description 6
- 238000011160 research Methods 0.000 claims description 6
- 239000000126 substance Substances 0.000 claims description 6
- 229910001727 uranium mineral Inorganic materials 0.000 claims description 6
- 230000000994 depressogenic effect Effects 0.000 claims description 4
- 230000005389 magnetism Effects 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 description 8
- 238000012937 correction Methods 0.000 description 8
- 230000008859 change Effects 0.000 description 7
- 150000001875 compounds Chemical class 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 4
- 230000005358 geomagnetic field Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000005553 drilling Methods 0.000 description 3
- 238000005065 mining Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241000251468 Actinopterygii Species 0.000 description 1
- 241001062954 Clinopodium Species 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Images
Classifications
-
- 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
Abstract
The invention relates to a method and a system for positioning a favorable ore-forming space of sandstone-type uranium ore in an anticline development area. Determining a space distribution form of favorable synthetic ore sand bodies, a redox transition zone position of a depression coverage area and a favorable synthetic ore area above an anticline rising area; determining a mineral-forming space favorable for the depression coverage area according to the position of the redox transition zone of the depression coverage area and the space spreading form of the mineral-forming sand body; determining the favorable ore-forming space of the anticline uplifting area according to the favorable ore-forming area above the anticline uplifting area and the space distribution form of the favorable ore-forming sand body; and determining that the uranium mineralization in the anticline development area is favorable for the mineralization space according to the depression coverage area favorable for the mineralization space and the anticline uplift area favorable for the mineralization space. The invention can accurately position the space position of the synthetic ore body near the anticline raised area.
Description
Technical Field
The invention relates to the field of uranium ore geological exploration, in particular to a method and a system for locating a favorable ore-forming space of sandstone-type uranium ore in an anticline development area.
Background
With the rapid development of the nuclear power industry in China, the demand of uranium resources is increasing. Sandstone-type uranium ores are one of main ore finding types in China in recent years, and the forming process of the sandstone-type uranium ores is influenced by various aspects such as structures, stratums, hydrological environments, alteration environments and the like. The anticline structure plays a vital role in the formation, enrichment, preservation and even destruction of sandstone-type uranium ores. With the gradual advance of sandstone-type uranium deposit exploration in China, a plurality of sandstone-type uranium deposit is found in a fold development area of a large number of sedimentary basins, and the fact that the uranium mineralization is closely related to the anticline is verified. However, the distribution rule of main ore control factors cannot be determined at present, particularly, the spatial position of an ore forming sand body beneficial to the vicinity of a anticline uplift area is difficult to determine, and a systematic ore forming space positioning method is lacked, so that the uranium ore exploration work of the type is difficult to make a major breakthrough. Therefore, a method and a system for positioning the favorable mineralization space of sandstone-type uranium ores in an anticline development area are needed to provide a basis for rapid space positioning of the favorable uranium ore areas in a fold development area.
Disclosure of Invention
The invention aims to provide a method and a system for positioning a favorable ore-forming space of sandstone-type uranium ore in an anticline development area, which can accurately position the space position of a favorable ore-forming sand body near an anticline uplift area.
In order to achieve the purpose, the invention provides the following scheme:
a method for locating favorable ore-forming space of sandstone-type uranium ore in an anticline development area comprises the following steps:
acquiring the stratum relief condition of a anticline development area, and dividing a anticline uplifting area and a depression coverage area according to the stratum relief condition of the research area;
determining resistivity characteristics according to the resistivity logging curve;
inverting a resistivity section equivalent graph according to an audio magnetotelluric sounding method;
determining a space distribution form favorable for mineral sand bodies according to the resistivity section equivalent diagram and the resistivity characteristics;
acquiring residual magnetic anomaly data;
determining the position of a redox transition zone of a depression coverage area according to the residual magnetic anomaly data;
acquiring soil radon concentration data;
determining a favorable mineralization area above the anticline upheaval area according to the radon concentration data of the soil;
determining a mineral-forming space of the depression coverage area according to the position of the redox transition zone of the depression coverage area and the space spreading form of the mineral-forming sand body;
determining favorable ore-forming space of the anticline uplifting area according to the favorable ore-forming area above the anticline uplifting area and the space distribution form of the favorable ore-forming sand body;
and determining the favorable ore-forming space of uranium ore-forming in the anticline development area according to the favorable ore-forming space of the depression coverage area and the favorable ore-forming space of the anticline upheaval area.
Optionally, the determining, according to the resistivity section equivalent diagram and the resistivity characteristic, a spatial distribution form favorable for mineral sand formation specifically includes:
determining the range of each section favorable for mineral sand bodies according to the resistivity section equivalent diagram and the resistivity characteristics;
and connecting the ranges of the sections favorable for the mineral sand bodies along the shortest path to obtain the space distribution form favorable for the mineral sand bodies.
Optionally, the obtaining of the remaining magnetic anomaly data specifically includes:
using formula according to the abnormal data of polarized magnetic field
And Δ TRemainder of＝ΔTPolarization plate-ΔTRegion(s)Obtaining residual magnetic anomaly data;
wherein, Delta TRegion(s)Regional magnetic anomalies, Sum Δ T, representing survey pointsPolarization plate(0, R) represents the sum of all measured polarized magnetic anomalies within a range of R from the measured point radius; n represents the number of measuring points in the range, R represents the distance of the window participating in calculation, and is at least 3 times of the distance between the measuring lines, and delta TRemainder ofRepresenting residual magnetic anomalies at the measurement points.
Optionally, the determining the position of the redox transition zone of the depression coverage area according to the residual magnetic anomaly data specifically includes:
adopting a formula according to the depression coverage area and the residual magnetism abnormal data
Determining the position of a redox transition zone of a down-cast coverage area;
wherein the content of the first and second substances,
representing the average value of residual abnormality of 20 percent measuring points with the maximum residual magnetic abnormal value of all measuring points,
representing the average value of residual abnormality of 20 percent of measuring points with the minimum residual magnetic abnormal value of all measuring points,
mean value, Δ T, representing residual magnetic anomalies of the measurement pointRemainder ofRepresenting residual magnetic anomalies at the measurement points.
Optionally, the determining, according to the radon concentration data of the soil, a favorable mineralization region above the anticline upheaval region specifically includes:
adopting a formula according to the soil radon concentration data
And
determining the abnormal lower limit value of radon concentration;
determining a favorable mineralization area above the anticline uplift area according to the anticline uplift area and the radon concentration abnormal lower limit value;
wherein the content of the first and second substances,
for smoothingThe average value of the data is shown, sigma is the standard deviation of the data after smoothing, A (R) is the abnormal lower limit value of radon concentration, and m and n are the total number of measuring lines and the number of measuring lines on each measuring line.
Optionally, the determining a favorable mineral forming space of the depressed coverage area according to the redox transition zone position of the depressed coverage area and the space spreading form of the favorable mineral sand body specifically includes:
intersecting a vertical projection zone of the redox transition zone position of the depression coverage area and the spatial distribution form of the mineral-mineral.
Optionally, the determining the favorable ore-forming space of the anticline humping area according to the favorable ore-forming area above the anticline humping area and the space distribution form of the favorable ore-forming sand body specifically includes:
and intersecting the vertical projection zone of the favorable ore-forming area above the anticline elevation area with the space distribution form of the favorable ore-forming sand body to obtain a superposed second space body, wherein the second space body is the favorable ore-forming space of the anticline elevation area.
Optionally, determining a favorable mineralization space of an anticline development area according to the favorable mineralization space of the depressed coverage area and the favorable mineralization space of the anticline upheaval rising area specifically includes:
and performing union treatment on the favorable mineral forming space of the depression coverage area and the favorable mineral forming space of the anticline upheaval area, and determining the favorable mineral forming space of the anticline developmental area.
A favorable ore-forming space positioning system for sandstone-type uranium ores in an anticline development area comprises:
the anticline development region dividing module is used for acquiring the relief condition of the anticline development region and dividing the anticline uplift region and the down-depression coverage region according to the relief condition of the study region;
the resistivity characteristic determining module is used for determining the resistivity characteristic according to the resistivity logging curve;
the resistivity section equivalent graph inversion module is used for inverting the resistivity section equivalent graph according to an audio magnetotelluric sounding method;
the space spreading form determination module is used for determining the space spreading form of the mineral sand bodies according to the resistivity section equivalent graph and the resistivity characteristics;
the residual magnetic anomaly data acquisition module is used for acquiring residual magnetic anomaly data;
a redox transition zone position determining module of the depression coverage area, configured to determine the redox transition zone position of the depression coverage area according to the residual magnetic anomaly data;
the soil radon concentration data acquisition module is used for acquiring the soil radon concentration data;
the favorable mineralization area determination module is used for determining the favorable mineralization area above the anticline upheaval area according to the soil radon concentration data;
the favorable mineral-forming space determination module of the depression coverage area is used for determining a favorable mineral-forming space of the depression coverage area according to the redox transition zone position of the depression coverage area and the space spreading form of the favorable mineral-forming sand body;
the favorable ore-forming space determination module of the anticline humping area is used for determining the favorable ore-forming space of the anticline humping area according to the favorable ore-forming area above the anticline humping area and the space distribution form of the favorable ore-forming sand body;
and the anticline development area uranium mineral favorable mineral forming space determining module is used for determining the anticline development area uranium mineral favorable mineral forming space according to the depression covering area favorable mineral forming space and the anticline rising area favorable mineral forming space.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
the method adopts a magnetic method, an audio geoelectromagnetic method, soil radon gas measurement and other comprehensive geophysical detection methods, can realize the rapid positioning of favorable areas for ore formation of the sandstone-type uranium ores in the anticline development area, reduce the ore search area, reduce unnecessary drilling work, accelerate the ore search period and realize the purposes of economic ore search and efficient ore search of the sandstone-type uranium ores in the anticline development area.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without inventive exercise.
FIG. 1 is a flow chart of a favorable ore-forming space positioning method for sandstone-type uranium ores in an anticline development area according to the invention;
fig. 2 is a structural diagram of a beneficial ore-formation space positioning system for sandstone-type uranium ore in an anticline development area.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention aims to provide a method and a system for positioning a favorable ore-forming space of sandstone-type uranium ore in an anticline development area, which can accurately position the space position of a favorable ore-forming sand body near an anticline uplift area.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
Fig. 1 is a flow chart of a favorable ore-forming space positioning method for sandstone-type uranium ores in an anticline development area. As shown in fig. 1, a method for locating a favorable mineralization space of sandstone-type uranium ore in an anticline development area includes:
step 101: and acquiring the stratum relief condition of a anticline development area, and dividing the anticline uplifting area and the depression coverage area according to the stratum relief condition of the research area.
Step 1011: and determining the working network degree of high-precision magnetic measurement according to geological features, wherein the line distance is not more than 1000 meters, the point distance is not more than 50 meters, the direction is vertical to the trend of the geologic body or the boundary line of the stratum, and the measurement area contains different exposed stratums as much as possible.
Step 1012: the used high-precision magnetometer is required to have stable performance, the resolution is better than 0.1nT, the repeatability precision is better than 0.2nT, the observation mean square error is better than 2nT, and the magnetic field data is collected and the coordinates and the elevation of the measuring point are stored at the same time.
Step 1013: and performing daily change correction and height correction on the obtained magnetic total field data by using the following formula, and removing a normal geomagnetic field to obtain the magnetic anomaly.
TImprovement of＝TMeasuring-(TDay(s)-T0)
ΔT＝TImprovement of-THeight of-TIs just
In the formula (I), the compound is shown in the specification,
Timprovement ofThe magnetic observed value after the daily variation correction is expressed in nT.
TMeasuringThe magnetic observations at the test points are in nT.
TDay(s)And the magnetic observed value of the day-to-day station is nT.
T0Is the basic magnetic field value of the day-changing station and has the unit of nT.
THeight ofThe unit is nT for a high correction of the magnetic field at the measuring point.
R is the earth mean radius, specifically 6371200 m.
H (I) is elevation of measurement point, and H (o) is elevation of base point in m.
TIs justAnd (4) obtaining the normal geomagnetic field value of the measured point by inquiring an international geomagnetic reference model IGRF, wherein the unit is nT.
And delta T is a magnetic abnormal value of a measuring point, and the unit is nT.
Step 1014: importing the magnetic anomaly data into Oasis Montaj software, and solving the polarized magnetic anomaly delta T by utilizing the polarized processing function in the Magmap modulePolarization plate. The selected local magnetic dip angle and magnetic declination parameter is obtained through international geomagnetic reference model IGRF inquiry.
Step 1015: and combining geological data, defining the region with stable and continuous polarimetric magnetic anomaly map, no singular point and abnormal amplitude change higher than 2/3 as a anticline uplifting region, and remaining as a down-depression covering region.
Step 102: and determining the resistivity characteristics according to the resistivity logging curve.
Step 103: and inverting the resistivity section equivalent graph according to an audio magnetotelluric sounding method.
And step 1031, the direction of the audio magnetotelluric measuring line is required to be consistent with that of the high-precision magnetic measurement in step 101, the number and the distance of the measuring lines are determined according to actual conditions, the point distance is generally less than or equal to 200 meters, and the measuring lines pass through exposed strata in the early age as far as possible.
Step 1032 is that the measured point data obtained in the step 1031 are imported into EMAGE-2D two-dimensional inversion software, Both and the rounding coefficient are selected as 3 in an inversion mode, the minimum error of the data is respectively 5% in a TM mode and 10% in a TE mode, and the resistivity data are obtained through automatic inversion calculation.
And 1033, loading the resistivity data of each section in the step 1032 by surfer software to obtain an audio magnetotelluric inversion resistivity section equivalent diagram.
Step 104: according to the resistivity section equivalent graph and the resistivity characteristics, determining a space distribution form favorable for mineral sand bodies, which specifically comprises the following steps:
and determining the range of each section favorable for mineral sand bodies according to the resistivity section equivalent graph and the resistivity characteristics. Specifically, the resistivity change range of the low-resistance region is determined according to the resistivity characteristics in the step 102, and an area with the thickness larger than 30m and stable resistivity change in the horizontal direction is selected from the resistivity section equivalent diagram, so that the range of each section which is favorable for forming ore sand bodies is determined.
Connecting the ranges of the sections which are beneficial to the mineral sand bodies along the shortest path to obtain the space distribution form which is beneficial to the mineral sand bodies and is marked as S (E).
Step 105: acquiring residual magnetic anomaly data, specifically comprising:
using formula according to the abnormal data of polarized magnetic field
And Δ TRemainder of＝ΔTPolarization plate-ΔTRegion(s)Obtaining residual magnetic anomaly data;
wherein, Delta TRegion(s)The regional magnetic anomaly of the measuring point is represented in unit nT; sum Δ TPolarization plate(0, R) represents the sum of all measured polarized magnetic anomalies within a range of R from the measured point radius; n represents the number of measuring points in the range, R represents the distance of the window participating in calculation, and is at least 3 times of the distance between the measuring lines, and delta TRemainder ofIndicates the residual magnetic anomaly at the measurement point in nT.
Step 106: determining the position of a redox transition zone of a down-depression coverage area according to the residual magnetic anomaly data, which specifically comprises the following steps:
adopting a formula according to the depression coverage area and the residual magnetism abnormal data
Determining the position of a redox transition zone of a down-cast coverage area;
wherein the content of the first and second substances,
representing the average value of residual abnormality of 20 percent measuring points with the maximum residual magnetic abnormal value of all measuring points,
representing the average value of residual abnormality of 20 percent of measuring points with the minimum residual magnetic abnormal value of all measuring points,
mean value, Δ T, representing residual magnetic anomalies of the measurement pointRemainder ofRepresenting residual magnetic anomalies at the measurement points.
Step 107: and acquiring the soil radon concentration data.
Step 1071: covering the anticline raised area defined in the step 101 in the working range of soil radon gas, wherein the line distance is not more than 200m, the point distance is not more than 50 m, the measurement time of each measuring point is not less than 10 minutes, and obtaining the radon concentration R of the measuring pointMeasuring(i,j)。
1072: radon concentration data R of each measuring pointMeasuringAnd (i, j) importing Oasis Montaj software, and obtaining smoothed data R (i, j) by using low-pass filtering processing in a Magmap module.
Step 108: according to soil radon concentration data, determining a favorable mineralization area above a anticline humping area, and specifically comprising the following steps:
step 1081: adopting a formula according to the soil radon concentration data
And
determining the abnormal lower limit value of radon concentration; wherein the content of the first and second substances,
the average value of the smoothed data is represented by Bq/m3(ii) a σ is the standard deviation of the smoothed data in Bq/m3(ii) a A (R) is the abnormal lower limit value of radon concentration and the unit is Bq/m3(ii) a m and n are the total number of measuring lines and the number of measuring points on each measuring line.
Step 1082: and determining a favorable mineralization area above the anticline uplifting area according to the anticline uplifting area and the radon concentration abnormal lower limit value. And loading the abnormal lower limit value of the radon concentration by surfer software to obtain a radon concentration contour map, and defining the area which is positioned in the range of the anticline humped area defined in the step 101 and has the radon concentration value higher than A (R) as a favorable mining area S (R).
Step 109: according to the redox transition zone position of the depression coverage area and the spatial spreading form of the mineral-:
intersecting a vertical projection zone of the redox transition zone position of the depression coverage area and the spatial distribution form of the mineral-mineral.
Step 110: according to the favorable ore-forming area above the anticline humping area and the space distribution form of the favorable ore-forming sand body, determining the favorable ore-forming space of the anticline humping area, and specifically comprising the following steps:
and intersecting the vertical projection zone of the favorable ore-forming area above the anticline elevation area with the space distribution form of the favorable ore-forming sand body to obtain a superposed second space body, wherein the second space body is the favorable ore-forming space of the anticline elevation area.
Step 111: according to the favorable mineralization space of the depression coverage area and the favorable mineralization space of the anticline upbeat area, determining the favorable mineralization space of the uranium in the anticline development area, which specifically comprises the following steps:
and performing union treatment on the favorable mineral forming space of the depression coverage area and the favorable mineral forming space of the anticline upheaval area, and determining the favorable mineral forming space of the anticline developmental area.
Fig. 2 is a structural diagram of a beneficial ore-formation space positioning system for sandstone-type uranium ore in an anticline development area. As shown in fig. 2, a system for locating a favorable mineralization space of a sandstone-type uranium ore in an anticline developmental zone includes:
the anticline development region dividing module 201 is used for acquiring the relief condition of the anticline development region and dividing the anticline uplift region and the depression coverage region according to the relief condition of the study region;
a resistivity characteristic determination module 202 for determining a resistivity characteristic according to the resistivity log;
the resistivity section equivalent graph inversion module 203 is used for inverting the resistivity section equivalent graph according to the audio magnetotelluric sounding method;
the spatial distribution form determination module 204 is used for determining the spatial distribution form of the favorable ore-forming sand bodies according to the resistivity section equivalent diagram and the resistivity characteristics;
a residual magnetic anomaly data acquisition module 205 configured to acquire residual magnetic anomaly data;
a redox transition zone position determining module 206 of the depression coverage area, configured to determine a redox transition zone position of the depression coverage area according to the residual magnetic anomaly data;
a soil radon concentration data acquisition module 207 for acquiring the soil radon concentration data;
a favorable-mineralization-region determining module 208 above the anticline upheaval region, configured to determine a favorable-mineralization region above the anticline upheaval region according to soil radon concentration data;
a mineral-favorable space determination module 209 for determining a mineral-favorable space of the depression coverage area according to the redox transition zone position of the depression coverage area and the spatial spread morphology of the mineral-favorable sand;
the favorable ore-forming space determination module 210 for determining the favorable ore-forming space of the anticline upheaval area according to the favorable ore-forming area above the anticline upheaval area and the space distribution form of the favorable ore-forming sand body;
and the anticline development area uranium mineral favorable mineral forming space determining module 211 is configured to determine a anticline development area uranium mineral favorable mineral forming space according to the depression coverage area favorable mineral forming space and the anticline rising area favorable mineral forming space.
Example 1:
taking the embodiment of the fish card area of the chada basin as an example, the method for positioning the favorable ore forming space of the sandstone-type uranium ore in the anticline development area sequentially comprises the following steps:
step 1: and (3) rapidly determining the stratum relief condition of the anticline development area by using high-precision magnetic measurement, and dividing the anticline uplifting area and the depression coverage area in the research area.
The step 1 comprises the following steps:
step 1.1, determining the working network degree of high-precision magnetic measurement according to geological features, wherein the line distance is 500 meters, the point distance is 50 meters, the direction is vertical to the trend of a geological body or a stratum boundary line, and a measurement area contains different exposed stratums as much as possible.
The high-precision magnetometer used in the step 1.2 has stable performance, the resolution is 0.01nT, and the repeatability precision is 0.1 nT. And observing the mean square error of 1.28nT, and acquiring magnetic field data and simultaneously storing the coordinates and the elevation of the measuring points.
And 1.3, carrying out daily variation correction and height correction on the obtained magnetic total field data by using the following formula, and removing a normal geomagnetic field to obtain magnetic anomaly.
TImprovement of＝TMeasuring-(TDay(s)-T0)
ΔT＝TImprovement of-THeight of-TIs just
In the formula (I), the compound is shown in the specification,
Timprovement ofThe magnetic observed value after the daily variation correction is expressed in nT.
TMeasuringThe magnetic observations at the test points are in nT.
TDay(s)And the magnetic observed value of the day-to-day station is nT.
T0Is the basic magnetic field value of the day-changing station and has the unit of nT.
THeight ofThe unit is nT for a high correction of the magnetic field at the measuring point.
R is the earth mean radius, specifically 6371200 m.
H (I) is elevation of measurement point, and H (o) is elevation of base point in m.
TIs justAnd (4) obtaining the normal geomagnetic field value of the measured point by inquiring an international geomagnetic reference model IGRF, wherein the unit is nT.
And delta T is a magnetic abnormal value of a measuring point, and the unit is nT.
Step 1.4, importing the magnetic anomaly data into Oasis Montaj software, and utilizing a Magmap module to carry out neutralizationPole processing function to solve pole magnetic anomaly delta TPolarization plate. And inquiring the local magnetic dip angle and the magnetic declination angle through an international geomagnetic reference model IGRF.
And step 1.5, combining geological data, and defining a stable, continuous and singularity-free region with the abnormal amplitude change higher than 2/3 of the polarized magnetic abnormal graph as a anticline uplifting region.
Step 2: and collecting drilling data of the research area, and determining the position and the resistivity characteristics of the sand body according to the resistivity logging curve.
And (3) collecting the information of the drilled hole ZKY-1 in the anticline uplift region defined in the step (1), and determining that the lithology of the low-resistance region is sandstone or sand-shale interbed according to the resistivity curve, wherein the specific resistivity variation range is 0-20 omega m.
And step 3: the space spreading form S (E) beneficial to forming ore sand bodies is defined by using an audio magnetotelluric sounding method.
The step 3 comprises the following steps:
and 3.1, the direction of the audio magnetotelluric survey line is required to be consistent with that of the high-precision magnetic survey in the step 1, the distance between two survey lines is 100 meters, and the survey lines pass through exposed ancient and Jurassic stratums.
And 3.2, importing the measured point data obtained in the step 3.1 into EMAGE-2D two-dimensional inversion software, selecting Both and rounding coefficients as 3 in an inversion mode, and automatically performing inversion calculation to obtain resistivity data by respectively adopting a TM mode of 5% and a TE mode of 10% for minimum data errors.
And 3.3, loading the resistivity data of each section in the step 3.2 by using surfer software to obtain an audio magnetotelluric inversion resistivity section equivalent diagram. And selecting an area with the resistivity range of 0-20 omega m, the thickness of more than 30m and stable resistivity change in the horizontal direction in the resistivity section equivalent diagram, and delineating the range of each section favorable for ore sand formation.
And 3.4, connecting the range which is circled on each section in the step 3.3 and is beneficial to mineral sand body along the shortest path to obtain the space distribution form S (E) of the mineral sand body.
And 4, step 4: and determining the redox transition position S (M) of the depression area by using the residual magnetic anomaly.
The step 4 comprises the following steps:
and 4.1, solving the residual magnetic anomaly. Polarizing magnetic anomaly delta T obtained in step 1Polarization plateThe calculation was processed as follows to obtain residual magnetic anomalies.
ΔTRemainder of＝ΔTPolarization plate-ΔTRegion(s)
In the formula (I), the compound is shown in the specification,
ΔTregion(s)The area representing the measurement point is magnetically abnormal, in nT.
SumΔTPolarization plate(0, R) represents the sum of all measured polarized magnetic anomalies within a range of R from the measured point radius; n represents the number of the measurement points in the range. And R is the window distance, and is 1500 meters.
ΔTRemainder ofIndicates the residual magnetic anomaly at the measurement point in nT.
And 4.2, loading the obtained residual magnetic anomaly constant by surfer software to obtain a residual anomaly map. A region which is located within the range of the circumscribed depression footprint of step 1 and in which the residual magnetic anomaly satisfies the following formula is simultaneously circumscribed as a redox transition zone position s (m).
In the formula (I), the compound is shown in the specification,
And 5: determining the area S (R) above the anticline uplifting area by using the radon concentration of soil, which is favorable for mineralization.
The step 5 comprises the following steps:
step 5.1 the working range of soil radon gas is required to cover the back slope range defined in step 1, the line distance is 200 meters, the point distance is 50 meters, the measurement time of each measuring point is 12 minutes, and the radon concentration R of the measuring point is obtainedMeasuring(i,j)。
Step 5.2, radon concentration data R of each measuring pointMeasuringAnd (i, j) importing Oasis Montaj software, and obtaining smoothed data R (i, j) by using low-pass filtering processing in a Magmap module.
And 5.3, calculating the radon abnormal delineation lower limit by using the following formula.
In the formula (I), the compound is shown in the specification,
the average value of the smoothed data is represented by Bq/m3；
σ is the standard deviation of the smoothed data in Bq/m3；
A (R) is the abnormal lower limit value of radon concentration and the unit is Bq/m3；
m and n are the total number of measuring lines and the number of measuring points on each measuring line.
And 5.4, loading the data by using surfer software to obtain a radon concentration contour map. And (3) the area which is positioned in the range of the dorsal clinopodium area determined in the step (1) and has the radon concentration value higher than A (R) is determined as a favorable mining area S (R).
Step 6: the formed ore is defined to be beneficial to the formed ore space.
Intersecting the depression coverage area vantage region S (M) defined in step 4 with the sand vantage region S (E) defined in step 3 by audio frequency earth electromagnetic sounding, the coinciding space being a depression coverage area mineral vantage space S (EM); intersecting the vertical projection band of the anticline upper favorable area S (R) defined in the step 5 with the favorable sand body area S (E) defined in the step 3 by the audio-frequency earth electromagnetic sounding, wherein the superposed space body is the favorable mining space S (ER) of the area.
S (EM) U.S (ER) is an ore-forming space favorable for uranium ore-forming in the anticline development area.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. For the system disclosed by the embodiment, the description is relatively simple because the system corresponds to the method disclosed by the embodiment, and the relevant points can be referred to the method part for description.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.
Claims (9)
1. A method for locating favorable ore-forming space of sandstone-type uranium ore in an anticline development area is characterized by comprising the following steps:
acquiring the stratum relief condition of a anticline development area, and dividing a anticline uplifting area and a depression coverage area according to the stratum relief condition of the research area;
determining resistivity characteristics according to the resistivity logging curve;
inverting a resistivity section equivalent graph according to an audio magnetotelluric sounding method;
determining a space distribution form favorable for mineral sand bodies according to the resistivity section equivalent diagram and the resistivity characteristics;
acquiring residual magnetic anomaly data;
determining the position of a redox transition zone of a depression coverage area according to the residual magnetic anomaly data;
acquiring soil radon concentration data;
determining a favorable mineralization area above the anticline upheaval area according to the radon concentration data of the soil;
determining a mineral-forming space of the depression coverage area according to the position of the redox transition zone of the depression coverage area and the space spreading form of the mineral-forming sand body;
determining favorable ore-forming space of the anticline uplifting area according to the favorable ore-forming area above the anticline uplifting area and the space distribution form of the favorable ore-forming sand body;
and determining the favorable ore-forming space of uranium ore-forming in the anticline development area according to the favorable ore-forming space of the depression coverage area and the favorable ore-forming space of the anticline upheaval area.
2. The method for positioning the favorable mineralization space of the sandstone-type uranium ore in the anticline developmental zone according to claim 1, wherein the method for determining the space distribution form of the favorable mineralization sand body according to the resistivity section equivalent diagram and the resistivity characteristics specifically comprises:
determining the range of each section favorable for mineral sand bodies according to the resistivity section equivalent diagram and the resistivity characteristics;
and connecting the ranges of the sections favorable for the mineral sand bodies along the shortest path to obtain the space distribution form favorable for the mineral sand bodies.
3. The method for positioning favorable mineralization space of the sandstone-type uranium ore in the anticline developmental zone according to claim 1, wherein the obtaining of the residual magnetic anomaly data specifically comprises:
according to the difference of polarized magnetic polesConstant data using formulas
And Δ TRemainder of＝ΔTPolarization plate-ΔTRegion(s)Obtaining residual magnetic anomaly data;
wherein, Delta TRegion(s)Regional magnetic anomalies, Sum Δ T, representing survey pointsPolarization plate(0, R) represents the sum of all measured polarized magnetic anomalies within a range of R from the measured point radius; n represents the number of measuring points in the range, R represents the distance of the window participating in calculation, and is at least 3 times of the distance between the measuring lines, and delta TRemainder ofRepresenting residual magnetic anomalies at the measurement points.
4. The method for positioning favorable mineralization space of a sandstone-type uranium ore in an anticline developmental zone according to claim 1, wherein the determining a redox transition position of a depression coverage area according to the residual magnetic anomaly data specifically comprises:
adopting a formula according to the depression coverage area and the residual magnetism abnormal data
Determining the position of a redox transition zone of a down-cast coverage area;
wherein the content of the first and second substances,
representing the average value of residual abnormality of 20 percent measuring points with the maximum residual magnetic abnormal value of all measuring points,
representing the average value of residual abnormality of 20 percent of measuring points with the minimum residual magnetic abnormal value of all measuring points,
mean value, Δ T, representing residual magnetic anomalies of the measurement pointRemainder ofRepresenting residual magnetic anomalies at the measurement points.
5. The method for positioning favorable mineralization space of a sandstone-type uranium ore in an anticline developmental zone according to claim 1, wherein the step of determining the favorable mineralization area above the anticline uplift zone according to soil radon concentration data comprises the following specific steps:
adopting a formula according to the soil radon concentration data
And
determining the abnormal lower limit value of radon concentration;
determining a favorable mineralization area above the anticline uplift area according to the anticline uplift area and the radon concentration abnormal lower limit value;
wherein the content of the first and second substances,
the average value of the smoothed data is sigma, the standard deviation of the smoothed data is A (R), the abnormal lower limit value of the flat radon concentration is A (R), and m and n are the total number of measuring lines and the number of measuring points on each measuring line.
6. The method for positioning a favored mineralization space of a sandstone-type uranium deposit in an anticline developmental zone according to claim 1, wherein the determining the favored mineralization space of the depression coverage zone according to the position of the redox transition zone of the depression coverage zone and the spatial spread morphology of the favored mineral sand comprises:
intersecting a vertical projection zone of the redox transition zone position of the depression coverage area and the spatial distribution form of the mineral-mineral.
7. The method for positioning the favorable mineralization space of the sandstone-type uranium ore in the anticline developmental zone according to claim 1, wherein the method for determining the favorable mineralization space in the anticline uplift zone according to the favorable mineralization area above the anticline uplift zone and the space distribution form of the favorable mineralization sand body comprises the following steps:
and intersecting the vertical projection zone of the favorable ore-forming area above the anticline elevation area with the space distribution form of the favorable ore-forming sand body to obtain a superposed second space body, wherein the second space body is the favorable ore-forming space of the anticline elevation area.
8. The method according to claim 1, wherein the determining the favorable ore-forming space for uranium mineralization in the anticline development area according to the favorable ore-forming space for the depressed covering area and the favorable ore-forming space for the anticline upheaval rise area specifically comprises:
and performing union treatment on the favorable mineral forming space of the depression coverage area and the favorable mineral forming space of the anticline upheaval area, and determining the favorable mineral forming space of the anticline developmental area.
9. A system for locating an ore-forming space favorable for sandstone-type uranium ores in an anticline development area is characterized by comprising:
the anticline development region dividing module is used for acquiring the relief condition of the anticline development region and dividing the anticline uplift region and the down-depression coverage region according to the relief condition of the study region;
the resistivity characteristic determining module is used for determining the resistivity characteristic according to the resistivity logging curve;
the resistivity section equivalent graph inversion module is used for inverting the resistivity section equivalent graph according to an audio magnetotelluric sounding method;
the space spreading form determination module is used for determining the space spreading form of the mineral sand bodies according to the resistivity section equivalent graph and the resistivity characteristics;
the residual magnetic anomaly data acquisition module is used for acquiring residual magnetic anomaly data;
a redox transition zone position determining module of the depression coverage area, configured to determine the redox transition zone position of the depression coverage area according to the residual magnetic anomaly data;
the soil radon concentration data acquisition module is used for acquiring the soil radon concentration data;
the favorable mineralization area determination module is used for determining the favorable mineralization area above the anticline upheaval area according to the soil radon concentration data;
the favorable mineral-forming space determination module of the depression coverage area is used for determining a favorable mineral-forming space of the depression coverage area according to the redox transition zone position of the depression coverage area and the space spreading form of the favorable mineral-forming sand body;
the favorable ore-forming space determination module of the anticline humping area is used for determining the favorable ore-forming space of the anticline humping area according to the favorable ore-forming area above the anticline humping area and the space distribution form of the favorable ore-forming sand body;
and the anticline development area uranium mineral favorable mineral forming space determining module is used for determining the anticline development area uranium mineral favorable mineral forming space according to the depression covering area favorable mineral forming space and the anticline rising area favorable mineral forming space.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110354196.4A CN113075749B (en) | 2021-04-01 | 2021-04-01 | Method and system for locating favorable ore-forming space of sandstone uranium deposit in anticline development area |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110354196.4A CN113075749B (en) | 2021-04-01 | 2021-04-01 | Method and system for locating favorable ore-forming space of sandstone uranium deposit in anticline development area |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113075749A true CN113075749A (en) | 2021-07-06 |
| CN113075749B CN113075749B (en) | 2024-03-22 |
Family
ID=76614598
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110354196.4A Active CN113075749B (en) | 2021-04-01 | 2021-04-01 | Method and system for locating favorable ore-forming space of sandstone uranium deposit in anticline development area |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113075749B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2668387A1 (en) * | 2001-04-24 | 2002-10-31 | Shell Canada Limited | In situ recovery from a tar sands formation |
| CN103837908A (en) * | 2014-03-05 | 2014-06-04 | 核工业北京地质研究院 | Rapid prospecting positioning method applicable to hidden sandstone-type uranium mine |
| WO2015052590A2 (en) * | 2013-10-10 | 2015-04-16 | Opto Trace (Suzhou) Technologies, Inc. | System and method for detecting crude oil or gas underground using light scattering spectral analyses |
| CN105388531A (en) * | 2015-10-19 | 2016-03-09 | 成都理工大学 | Lithology identification method based on support vector regression machine and Kernel Fisher discriminant analysis |
| CN106646660A (en) * | 2015-10-30 | 2017-05-10 | 核工业北京地质研究院 | Sandstone type uranium mine integrated geophysical exploration method |
| CN107479100A (en) * | 2017-08-24 | 2017-12-15 | 山东省物化探勘查院 | A kind of oblique deep concealed iron ore method of prospecting of steep dip |
| CN110908010A (en) * | 2019-10-25 | 2020-03-24 | 东华理工大学 | Effective geophysical method for finding shallow sandstone type uranium deposit within 800 meters |
-
2021
- 2021-04-01 CN CN202110354196.4A patent/CN113075749B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2668387A1 (en) * | 2001-04-24 | 2002-10-31 | Shell Canada Limited | In situ recovery from a tar sands formation |
| WO2015052590A2 (en) * | 2013-10-10 | 2015-04-16 | Opto Trace (Suzhou) Technologies, Inc. | System and method for detecting crude oil or gas underground using light scattering spectral analyses |
| CN103837908A (en) * | 2014-03-05 | 2014-06-04 | 核工业北京地质研究院 | Rapid prospecting positioning method applicable to hidden sandstone-type uranium mine |
| CN105388531A (en) * | 2015-10-19 | 2016-03-09 | 成都理工大学 | Lithology identification method based on support vector regression machine and Kernel Fisher discriminant analysis |
| CN106646660A (en) * | 2015-10-30 | 2017-05-10 | 核工业北京地质研究院 | Sandstone type uranium mine integrated geophysical exploration method |
| CN107479100A (en) * | 2017-08-24 | 2017-12-15 | 山东省物化探勘查院 | A kind of oblique deep concealed iron ore method of prospecting of steep dip |
| CN110908010A (en) * | 2019-10-25 | 2020-03-24 | 东华理工大学 | Effective geophysical method for finding shallow sandstone type uranium deposit within 800 meters |
Non-Patent Citations (3)
| Title |
|---|
| 廉康;赵兴齐;王继斌;张淼;: "柴北缘上干柴沟组砂体特征及铀成矿作用研究", 东华理工大学学报(自然科学版), no. 04 * |
| 王会波;陈载林;: "剩余磁异常的计算与深部磁性盲矿层的预测", 中国锰业, no. 06 * |
| 黄国龙, 陈贤春: "柴达木盆地西北部可地浸砂岩型铀矿成矿条件及找矿方向", 铀矿地质, no. 06 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113075749B (en) | 2024-03-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Grayver et al. | 3D inversion and resolution analysis of land-based CSEM data from the Ketzin CO 2 storage formation | |
| Wang et al. | A new edge recognition technology based on the normalized vertical derivative of the total horizontal derivative for potential field data | |
| CN109425896A (en) | Dolomite oil and gas reservoir distribution forecasting method and device | |
| CN107329171A (en) | Depth Domain reservoir seismic inversion method and device | |
| Baumann-Wilke et al. | P-wave traveltime tomography for a seismic characterization of black shales at shallow depth on Bornholm, Denmark | |
| CN105510993A (en) | Foreland basin deep buried and compressed type complex gypsum-salt rock identification and distribution prediction method | |
| US20120296618A1 (en) | Multiscale Geologic Modeling of a Clastic Meander Belt Including Asymmetry Using Multi-Point Statistics | |
| CN105738952B (en) | A kind of horizontal well region reservoir rock phase modeling method | |
| CN106405664B (en) | A kind of magnetic anomaly normalizing pole method | |
| CN106650192B (en) | A kind of Volcanic Type Uranium Deposits magnetic interface inversion method | |
| CN105158796A (en) | Method and device for determining TOC content | |
| CN109839670A (en) | A kind of hydrothermal-type uranium deposit substrate interface inversion method | |
| US10451769B2 (en) | Method for determining petrophysical properties from logging measurements | |
| Andrieux et al. | The so-called static corrections in magnetotelluric measurements | |
| Zhu et al. | Normalized vertical derivatives in the edge enhancement of maximum-edge-recognition methods in potential fields | |
| Whitelaw et al. | High-resolution gravity study of the Gray Fossil Site | |
| Witter et al. | 2D and 3D Potential field mapping and modelling at the Fallon FORGE site, Nevada, USA | |
| CN113075749A (en) | Method and system for locating favorable ore-forming space of sandstone-type uranium ore in anticline development area | |
| Obiora et al. | Interpretation of aeromagnetic data of Idah area in north central Nigeria using combined methods | |
| Szarka | Geophysical mapping by stationary electric and magnetic field components: a combination of potential gradient mapping and magnetometric resistivity (MMR) methods | |
| CN105242317B (en) | A kind of determination method and device of velocity of longitudinal wave | |
| Maus et al. | Variogram analysis of helicopter magnetic data to identify paleochannels of the Omaruru River, Namibia | |
| Akingboye et al. | Development of Novel Velocity-Resistivity Geostatistical Relationships for Tropical Granitic Environments | |
| WO2012060888A1 (en) | System and method for providing a physical property model | |
| Rabeh et al. | Interpretation case study of the Sahl El Qaa area, southern Sinai Peninsula, Egypt |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |