CN113075749B - Method and system for locating favorable ore-forming space of sandstone uranium deposit in anticline development area - Google Patents
Method and system for locating favorable ore-forming space of sandstone uranium deposit in anticline development area Download PDFInfo
- Publication number
- CN113075749B CN113075749B CN202110354196.4A CN202110354196A CN113075749B CN 113075749 B CN113075749 B CN 113075749B CN 202110354196 A CN202110354196 A CN 202110354196A CN 113075749 B CN113075749 B CN 113075749B
- Authority
- CN
- China
- Prior art keywords
- anticline
- ore
- space
- favorable
- area
- 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.)
- Active
Links
- 230000002349 favourable effect Effects 0.000 title claims abstract description 141
- 238000011161 development Methods 0.000 title claims abstract description 54
- 238000000034 method Methods 0.000 title claims abstract description 52
- 229910052770 Uranium Inorganic materials 0.000 title claims abstract description 47
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 title claims abstract description 47
- 239000004576 sand Substances 0.000 claims abstract description 52
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 32
- 230000007704 transition Effects 0.000 claims abstract description 29
- 238000009826 distribution Methods 0.000 claims abstract description 23
- 230000007480 spreading Effects 0.000 claims abstract description 17
- 238000003892 spreading Methods 0.000 claims abstract description 17
- 230000018109 developmental process Effects 0.000 claims description 51
- 229910052704 radon Inorganic materials 0.000 claims description 45
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 claims description 45
- 230000009286 beneficial effect Effects 0.000 claims description 33
- 238000005259 measurement Methods 0.000 claims description 27
- 230000002159 abnormal effect Effects 0.000 claims description 26
- 239000002689 soil Substances 0.000 claims description 24
- 239000000126 substance Substances 0.000 claims description 22
- 230000033558 biomineral tissue development Effects 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 7
- 238000005266 casting Methods 0.000 claims description 6
- 230000005856 abnormality Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 description 8
- 238000012937 correction Methods 0.000 description 8
- 230000006872 improvement Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000011160 research Methods 0.000 description 6
- 230000005358 geomagnetic field Effects 0.000 description 4
- 238000005553 drilling Methods 0.000 description 3
- 210000000746 body region Anatomy 0.000 description 2
- 230000000994 depressogenic effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000010276 construction Methods 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
- 230000008676 import Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
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
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to a method and a system for locating favorable ore-forming space of sandstone-type uranium ores in anticline development areas. Determining a spatial distribution form of favorable ore bodies, a redox transition zone position of a depression coverage area and a favorable ore formation area above a anticline uplift area; determining favorable ore 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 favorable ore forming sand bodies; determining the favorable ore forming space of the anticline bulge region according to the favorable ore forming region above the anticline bulge region and the space spreading form of favorable ore forming sand bodies; and determining the uranium ore-forming favorable ore-forming space of 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 uplift area. The invention can accurately position the space position of the ore forming sand body near the anticline uplift 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 favorable ore formation space of sandstone-type uranium ores in a anticline development area.
Background
Along with the rapid development of the nuclear power industry in China, the uranium resource demand is larger and larger. Sandstone-type uranium ores are one of the main prospecting types in recent years in China, and the formation process is influenced by various aspects such as construction, stratum, hydrologic environment, alteration environment and the like. Wherein the anticline structure plays a vital role in the formation, enrichment, preservation and even destruction of sandstone-type uranium ores. Along with the gradual advancement of sandstone-type uranium deposit investigation in China, a plurality of sandstone-type uranium deposits are found in a large number of deposit basin fold development areas, and the close relationship between uranium deposit and anticline is verified. However, the distribution rule of main ore control factors cannot be clarified at present, particularly, the space position of the favorable ore-forming sand body near the anticline uplift area is difficult to determine, and a systematic ore-forming space positioning method is lacked, so that the type of uranium ore exploration work is difficult to obtain a great breakthrough. Therefore, a method and a system for locating the favorable ore formation space of sandstone uranium ores in a anticline development area are needed to provide a basis for the rapid spatial location of the favorable uranium ore formation area in the 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 uranium ores in a anticline development area, which can accurately position the space position of favorable ore-forming sand bodies near the anticline bulge area.
In order to achieve the above object, the present invention provides the following solutions:
the method for locating the favorable ore-forming space of sandstone-type uranium deposit in anticline development area comprises the following steps:
acquiring stratum relief conditions of a anticline development area, and dividing anticline bulge areas and depression coverage areas according to the stratum relief conditions of the research area;
determining resistivity characteristics according to the resistivity log;
inverting the resistivity section contour map according to an audio magnetotelluric sounding method;
determining the space spreading form of the favorable ore forming sand body according to the resistivity section contour map and the resistivity characteristics;
obtaining residual magnetic anomaly data;
determining the redox transition zone position of the depression coverage area according to the residual magnetic anomaly data;
acquiring the radon concentration data of the soil;
determining an advantageous ore forming area above the anticline uplift area according to the soil radon concentration data;
determining a favorable ore-forming space of the depression coverage area according to the redox transition zone position of the depression coverage area and the space distribution form of the favorable ore-forming sand body;
determining a favorable ore forming space of the anticline uplift region according to the favorable ore forming region above the anticline uplift region and the space distribution form of the favorable ore forming sand body;
and determining the beneficial ore-forming space of uranium in the anticline development area according to the beneficial ore-forming space of the depression coverage area and the beneficial ore-forming space of the anticline uplift area.
Optionally, the determining the space distribution form of the favorable ore-forming sand body according to the resistivity section contour map and the resistivity characteristics specifically comprises the following steps:
determining the range of each section favorable for the ore formation sand body according to the resistivity section contour map and the resistivity characteristics;
and connecting the ranges of the favorable ore-forming sand bodies of the sections along the shortest paths to obtain the space distribution form favorable for the ore-forming sand bodies.
Optionally, the acquiring the residual magnetic anomaly data specifically includes:
adopting a formula according to the abnormal data of the chemical pole magnetismAnd DeltaT Residual of ＝ΔT Chemical electrode -ΔT Region(s) Obtaining residual magnetic anomaly data;
wherein DeltaT Region(s) SumDeltaT representing regional magnetic anomalies at the measurement site Chemical electrode (0, R) represents the sum of all the pole magnetic anomalies of the measuring point within the range of the radius R of the measuring point; n represents the number of measuring points in the range, R represents the window distance involved in calculation, which is at least 3 times of the distance between the measuring lines, and delta T Residual of Indicating the residual magnetic anomalies at the measurement point.
Optionally, determining the redox transition zone position 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 magnetic anomaly dataDetermining a redox transition zone location of the depression footprint;
wherein,indicating the residual anomaly average value of 20% of the measuring points with the maximum residual magnetic anomaly value of all the measuring points, ++>The residual abnormal average value of 20% measuring points with the minimum residual magnetic abnormal value of all measuring points is represented,mean value of residual magnetic anomalies of measuring points, delta T Residual of Indicating the residual magnetic anomalies at the measurement point.
Optionally, the determining the favorable ore forming area above the anticline uplift area according to the soil radon concentration data specifically comprises the following steps:
adopting a formula according to the soil radon concentration dataAnd->Determining an abnormal radon concentration lower limit value;
determining an advantageous ore formation area above the anticline bulge area according to the anticline bulge area and the radon concentration abnormal lower limit value;
wherein,for the average value of the smoothed data, sigma is the standard deviation of the smoothed data, A (R) is the abnormal lower limit value of radon concentration, and m and n are the total line number and the number of measurement points on each line.
Optionally, the determining the depression coverage area favorable ore forming space according to the redox transition zone position of the depression coverage area and the spatial distribution form of the favorable ore forming sand body specifically comprises:
and intersecting the vertical casting belt at the redox transition belt position of the depression coverage area with the space spreading form of the favorable ore formation sand body to obtain a superposed first space body, wherein the first space body is the favorable ore formation space of the depression coverage area.
Optionally, the determining the beneficial ore forming space of the anticline uplift area according to the beneficial ore forming area above the anticline uplift area and the spatial distribution form of the beneficial ore forming sand body specifically comprises:
and intersecting the vertical casting belt of the favorable ore forming area above the anticline uplift area with the space spreading form of the favorable ore forming sand body to obtain a superposed second space body, wherein the second space body is a favorable ore forming space of the anticline uplift area.
Optionally, the determining the anticline development area uranium ore-forming beneficial ore-forming space according to the depression coverage area beneficial ore-forming space and the anticline uplift area beneficial ore-forming space specifically includes:
and carrying out union processing on the favorable ore forming space of the depression coverage area and the favorable ore forming space of the anticline uplift area, and determining the favorable ore forming space of uranium ore forming in the anticline development area.
A system for locating a favorable ore-forming space of sandstone-type uranium deposit in a anticline development region, comprising:
the anticline development area dividing module is used for acquiring the anticline development area stratum relief condition and dividing anticline bulge areas and depression coverage areas according to the research area stratum relief condition;
the resistivity characteristic determining module is used for determining resistivity characteristics according to the resistivity logging curve;
the resistivity section contour map inversion module is used for inverting the resistivity section contour map according to an audio magnetotelluric sounding method;
the space spreading form determining module is used for determining the space spreading form of the favorable ore-forming sand body according to the resistivity section contour map and the resistivity characteristics;
the residual magnetic abnormal data acquisition module is used for acquiring residual magnetic abnormal data;
the redox transition zone position determining module is used for determining 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 beneficial ore forming area determining module is used for determining the beneficial ore forming area above the anticline bulge area according to the soil radon concentration data;
a depression coverage area favorable ore forming space determining module, configured to determine a depression coverage area favorable ore forming space according to a redox transition zone position of the depression coverage area and a spatial distribution form of the favorable ore forming sand body;
the anticline uplift region favorable ore forming space determining module is used for determining the anticline uplift region favorable ore forming space according to the favorable ore forming region above the anticline uplift region and the space distribution form of the favorable ore forming sand body;
and the determination module is used for determining the beneficial ore forming space of the uranium ore forming in the anticline development zone according to the beneficial ore forming space of the depression coverage area and the beneficial ore forming space of the anticline uplift zone.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
the comprehensive geophysical detection methods such as a magnetic method, an audio magnetotelluric method and soil radon measurement are adopted, so that the rapid positioning of the favorable region for ore formation of sandstone-type uranium ores in the anticline development region can be realized, the ore-finding region is reduced, unnecessary drilling work is reduced, the ore-finding period is accelerated, and the purposes of economic ore finding and efficient ore finding of the sandstone-type uranium ores in the anticline development region are realized.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions of the prior art, the drawings that are needed in the embodiments will be briefly described below, it being obvious that the drawings in the following description are only some embodiments of the present invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
Fig. 1 is a flow chart of a method for locating a favorable ore-forming space of sandstone-type uranium deposit in a anticline development area according to the present invention;
fig. 2 is a block diagram of a favorable ore-forming space positioning system of sandstone-type uranium ores in anticline development areas.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
The invention aims to provide a method and a system for positioning a favorable ore-forming space of sandstone uranium ores in a anticline development area, which can accurately position the space position of favorable ore-forming sand bodies near the anticline bulge area.
In order that the above-recited objects, features and advantages of the present invention will become more readily apparent, a more particular description of the invention will be rendered by reference to the appended drawings and appended detailed description.
Fig. 1 is a flow chart of a method for locating a favorable ore-forming space of sandstone-type uranium deposit in a anticline development area. As shown in fig. 1, a method for locating a favorable ore-forming space of sandstone-type uranium deposit in anticline development area comprises the following steps:
step 101: and acquiring the stratum fluctuation condition of the anticline development area, and dividing the anticline uplift area and the depression coverage area according to the stratum fluctuation condition of the research area.
Step 1011: and determining the working network of the high-precision magnetic measurement according to the geological features, wherein the line distance is not more than 1000 meters, the point distance is not more than 50 meters, the direction is perpendicular to the trend of the geological body or the stratum boundary line, and the measuring area comprises different exposed stratum as much as possible.
Step 1012: the high-precision magnetometer used has stable performance, resolution superior to 0.1nT, repeatability precision superior to 0.2nT, observation mean square error superior to 2nT, and magnetic field data acquisition and coordinate and elevation storage of the measuring point.
Step 1013: and performing daily change correction and altitude correction on the obtained magnetic total field data by using the following method, and removing the normal geomagnetic field to obtain magnetic anomalies.
T Improvement of ＝T Measuring -(T Day of the day -T 0 )
ΔT＝T Improvement of -T High height -T Positive direction
In the method, in the process of the invention,
T improvement of The unit is nT for the magnetic force observation value after the daily change correction.
T Measuring The unit is nT, which is the observed value of magnetic force of the measuring point.
T Day of the day Is the magnetometric observation value of the solar power station, and the unit is nT.
T 0 The basic magnetic field value of the solar power station is given as nT.
T High height For height correction of the magnetic field of the measuring point, the unit is nT.
R is the average radius of the earth, and is specifically 6371200m.
H (I) is the elevation of the measuring point, H (o) is the elevation of the base point, and the unit is m.
T Positive direction And obtaining the normal geomagnetic field value of the measuring point through the IGRF query of the international geomagnetic reference model, wherein the unit is nT.
Delta T is the magnetic anomaly value of the measuring point, and the unit is nT.
Step 1014: the magnetic anomaly data are imported into Oasis Montaj software, and the polarized magnetic anomaly delta T is obtained by utilizing the polarized processing function in the Magmap module Chemical electrode . The selected local geomagnetic inclination angle and declination angle parameters are obtained through the IGRF query of an international geomagnetic reference model.
Step 1015: and (3) combining geological data, and delineating a region which has stable, continuous and no singular point and has an abnormal amplitude change higher than 2/3 of the maximum magnetic anomaly map as a anticline bulge region and the rest as a depression coverage region.
Step 102: from the resistivity log, a resistivity characteristic is determined.
Step 102 determines the horizon and resistivity characteristics of the beneficial sand body from the resistivity log by collecting the study area borehole data. The selected borehole should be located in the anticline bump area outlined in step 101, and the lithology revealed by the borehole is rich in more sandstone or sand-mud interbedded. And determining the lithology of the low-resistance area and the change range of specific resistivity according to the resistivity curve.
Step 103: inverting the resistivity section contour map according to an audio magnetotelluric sounding method.
The audio magnetotelluric line direction in step 1031 needs to be consistent with the high-precision magnetotelluric line in step 101, the number and the distance between the line are determined according to practical conditions, the point distance is generally less than or equal to 200 meters, and the line passes through the exposed stratum as early as possible.
Step 1032 is to import the measurement point data obtained in step 1031 into EMAGE-2D two-dimensional inversion software, select the inversion mode to be Both and the smooth coefficient to be 3, and the minimum error of the data is 5% of the TM mode and 10% of the TE mode respectively, and the resistivity data is obtained by automatic inversion calculation.
Step 1033, loading the section resistivity data of step 1032 by using surfer software to obtain an audio magnetotelluric inversion resistivity section contour map.
Step 104: and determining the space distribution form of the favorable ore-forming sand body according to the resistivity section contour map and the resistivity characteristics, wherein the method specifically comprises the following steps of:
and determining the range of favorable ore formation of each section according to the resistivity section contour map and the resistivity characteristics. Specifically, the resistivity variation range of the low-resistance area is determined according to the resistivity characteristics in the step 102, a region with the thickness larger than 30m and stable resistivity variation in the horizontal direction is selected from the contour map of the resistivity section, and the range of favorable ore formation of each section is defined.
And connecting the ranges of the favorable ore-forming sand bodies of the sections along the shortest paths to obtain the space distribution form of the favorable ore-forming sand bodies, which is marked as S (E).
Step 105: the method for acquiring the residual magnetic anomaly data specifically comprises the following steps:
adopting a formula according to the abnormal data of the chemical pole magnetismAnd DeltaT Residual of ＝ΔT Chemical electrode -ΔT Region(s) Obtaining residual magnetic anomaly data;
wherein DeltaT Region(s) The magnetic abnormality of the region of the measuring point is represented by a unit nT; sum delta T Chemical electrode (0, R) represents the sum of all the pole magnetic anomalies of the measuring point within the range of the radius R of the measuring point; n represents the number of measuring points in the range, R represents the window distance involved in calculation, which is at least 3 times of the distance between the measuring lines, and delta T Residual of The residual magnetic anomalies at the measurement point are expressed in nT.
Step 106: determining a redox transition zone position of the depression coverage area according to the residual magnetic anomaly data, wherein the method specifically comprises the following steps of:
adopting a formula according to the depression coverage area and the residual magnetic anomaly dataDetermining a redox transition zone location of the depression footprint;
wherein,indicating the residual anomaly average value of 20% of the measuring points with the maximum residual magnetic anomaly value of all the measuring points, ++>The residual abnormal average value of 20% measuring points with the minimum residual magnetic abnormal value of all measuring points is represented,mean value of residual magnetic anomalies of measuring points, delta T Residual of Indicating the residual magnetic anomalies at the measurement point.
Step 107: and acquiring the radon concentration data of the soil.
Step 1071: the soil radon gas working range needs to cover the anticline uplift area outlined in the step 101, the line distance is not more than 200 meters, the point distance is not more than 50 meters, the measurement time of each measuring point is not less than 10 minutes, and the measuring points are obtainedRadon concentration R Measuring (i,j)。
Step 1072: radon concentration data R of each measuring point Measuring (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 the radon concentration data of the soil, the method for determining the favorable ore forming area above the anticline uplift area comprises the following steps:
step 1081: adopting a formula according to the soil radon concentration data And->Determining an abnormal radon concentration lower limit value; wherein (1)>The unit is Bq/m, which is the average value of the smoothed data 3 The method comprises the steps of carrying out a first treatment on the surface of the Sigma is the standard deviation of the smoothed data in Bq/m 3 The method comprises the steps of carrying out a first treatment on the surface of the A (R) is radon concentration anomaly lower limit value, and the unit is Bq/m 3 The method comprises the steps of carrying out a first treatment on the surface of the m and n are the total number of the measuring lines and the number of the measuring points on each measuring line.
Step 1082: and determining an advantageous ore forming area above the anticline bulge area according to the anticline bulge area and the radon concentration abnormal lower limit value. Loading the radon concentration anomaly lower limit value with surfer software to obtain a radon concentration contour map, and delineating a region which is positioned in the range of the anticline bulge region delineated in the step 101 and has a radon concentration value higher than A (R) as an advantageous mining area S (R).
Step 109: determining a favorable ore-forming space of the depression footprint from the redox transition zone position of the favorable ore-forming sand body and the spatial distribution morphology of the favorable ore-forming sand body, specifically comprising:
and intersecting the vertical casting belt at the redox transition belt position of the depression coverage area with the space spreading form of the favorable ore formation sand body to obtain a superposed first space body, wherein the first space body is the favorable ore formation space of the depression coverage area.
Step 110: determining the favorable ore forming space of the anticline uplift region according to the favorable ore forming region above the anticline uplift region and the space distribution form of the favorable ore forming sand body, and specifically comprising:
and intersecting the vertical casting belt of the favorable ore forming area above the anticline uplift area with the space spreading form of the favorable ore forming sand body to obtain a superposed second space body, wherein the second space body is a favorable ore forming space of the anticline uplift area.
Step 111: determining a back-inclined development zone uranium ore-forming favorable ore-forming space according to the depression coverage favorable ore-forming space and the back-inclined bulge zone favorable ore-forming space, specifically comprising:
and carrying out union processing on the favorable ore forming space of the depression coverage area and the favorable ore forming space of the anticline uplift area, and determining the favorable ore forming space of uranium ore forming in the anticline development area.
Fig. 2 is a block diagram of a favorable ore-forming space positioning system of sandstone-type uranium ores in anticline development areas. As shown in fig. 2, a system for locating a favorable ore-forming space of sandstone-type uranium deposit in a anticline development area comprises:
the anticline development area dividing module 201 is configured to obtain anticline development area stratum relief conditions, and divide anticline bulge areas and depression coverage areas according to the research area stratum relief conditions;
a resistivity characteristic determination module 202 for determining a resistivity characteristic from the resistivity log;
the resistivity section contour map inversion module 203 is configured to invert the resistivity section contour map according to an audio magnetotelluric sounding method;
a spatial spreading morphology determination module 204 for determining a spatial spreading morphology of the favorable ore body based on the resistivity profile contour map and the resistivity features;
the residual magnetic anomaly data acquisition module 205 is configured to acquire residual magnetic anomaly data;
a redox transition zone position determination module 206 for determining a redox transition zone position of the depressed coverage zone from the residual magnetic anomaly data;
a soil radon concentration data acquisition module 207 for acquiring the soil radon concentration data;
an advantageous mineralisation zone determination module 208 above the anticline raised zone for determining an advantageous mineralisation zone above the anticline raised zone from soil radon concentration data;
a depression footprint favorable mineralisation space determination module 209 for determining a depression footprint favorable mineralisation space from a redox transition zone position of the depression footprint and a spatial spread morphology of the favorable mineralisation sand;
a anticline bump area favorable ore forming space determining module 210, configured to determine a anticline bump area favorable ore forming space according to a favorable ore forming area above the anticline bump area and a spatial distribution form of the favorable ore forming sand body;
the beneficial mineralization space determining module 211 for uranium mineralization in the anticline development area is used for determining the beneficial mineralization space for uranium mineralization in the anticline development area according to the beneficial mineralization space in the depression coverage area and the beneficial mineralization space in the anticline uplift area.
Example 1:
taking a firewood basin fish-calorie region as an example, the positioning method of the sandstone-type uranium deposit favorable to the ore formation space in the anticline development region provided by the invention sequentially comprises the following steps:
step 1: and rapidly determining stratum fluctuation conditions of the anticline development area by utilizing high-precision magnetic measurement, and dividing anticline uplift areas and depression coverage areas in the research area.
The step 1 comprises the following steps:
step 1.1, determining high-precision magnetic measurement working network according to geological features, wherein the line distance is 500 meters, the point distance is 50 meters, the direction is perpendicular to the trend of the geologic body or the stratum boundary line, and the measurement area comprises different exposed stratum 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.1nT. And observing the mean square error of 1.28nT, collecting magnetic field data and simultaneously storing coordinates and elevation of the measuring point.
And 1.3, performing daily change correction and altitude correction on the obtained magnetic total field data by using the following formula, and removing the normal geomagnetic field to obtain magnetic anomalies.
T Improvement of ＝T Measuring -(T Day of the day -T 0 )
ΔT＝T Improvement of -T High height -T Positive direction
In the method, in the process of the invention,
T improvement of The unit is nT for the magnetic force observation value after the daily change correction.
T Measuring The unit is nT, which is the observed value of magnetic force of the measuring point.
T Day of the day Is the magnetometric observation value of the solar power station, and the unit is nT.
T 0 The basic magnetic field value of the solar power station is given as nT.
T High height For height correction of the magnetic field of the measuring point, the unit is nT.
R is the average radius of the earth, and is specifically 6371200m.
H (I) is the elevation of the measuring point, H (o) is the elevation of the base point, and the unit is m.
T Positive direction And obtaining the normal geomagnetic field value of the measuring point through the IGRF query of the international geomagnetic reference model, wherein the unit is nT.
Delta T is the magnetic anomaly value of the measuring point, and the unit is nT.
Step 1.4, the magnetic anomaly data are imported into Oasis Montaj software, and the polarized magnetic anomaly delta T is obtained by utilizing the polarized processing function in the Magmap module Chemical electrode . And inquiring the local geomagnetic inclination angle and the declination angle through an international geomagnetic reference model IGRF.
And 1.5, combining geological data, and delineating a region which has stable, continuous and no singular point and is higher than 2/3 abnormal amplitude change of the polarized magnetic anomaly map as a anticline bulge region.
Step 2: and collecting drilling data of a research area, and determining the horizon and the resistivity characteristics of the favorable sand body according to the resistivity logging curve.
And (3) collecting drilling ZKY-1 data in the anticline uplift area outlined in the step (1), and determining that the lithology of the low-resistance area is sandstone or sand-shale interbedded according to the resistivity curve, wherein the specific resistivity change range is 0-20 omega m.
Step 3: and (3) defining the space spreading form S (E) of the favorable ore forming sand body by utilizing an audio magnetotelluric sounding method.
The step 3 comprises the following steps:
and 3.1, the direction of the audio magnetotelluric line is required to be consistent with that of the high-precision magnetotelluric line in the step 1, the two measuring lines are separated by 100 meters, and the measuring lines pass through exposed ancient and dwarf stratum.
And 3.2, importing the measurement point data obtained in the step 3.1 into EMAGE-2D two-dimensional inversion software, selecting a Both inversion mode and a round coefficient to be 3, respectively adopting a TM mode and a TE mode for the minimum error of the data to be 5% and 10% respectively, and automatically inverting and calculating to obtain resistivity data.
And 3.3, loading the section resistivity data of the step 3.2 by using surfer software to obtain an audio magnetotelluric inversion resistivity section contour map. And selecting a region with the resistivity ranging from 0 to 20 omega m and the thickness being more than 30m and the resistivity variation in the horizontal direction being stable from the resistivity section contour map, and defining the range of favorable ore forming sand bodies of each section.
And 3.4, connecting the range of the favorable ore-forming sand bodies defined on each section of the step 3.3 along the shortest path to obtain the space distribution form S (E) of the favorable ore-forming sand bodies.
Step 4: the remaining magnetic anomalies are used to determine the redox transition zone position S (M) of the depression coverage area.
The step 4 comprises the following steps:
and 4.1, obtaining residual magnetic anomalies. The magnetic anomaly delta T of the chemical pole obtained in the step 1 is processed Chemical electrode The calculation is processed as follows to obtain the residual magnetic anomaly.
ΔT Residual of ＝ΔT Chemical electrode -ΔT Region(s)
In the method, in the process of the invention,
ΔT region(s) The area magnetic anomaly of the measurement point is represented by the unit nT.
SumΔT Chemical electrode (0, R) represents the sum of all the pole magnetic anomalies of the measuring point within the range of the radius R of the measuring point; n represents the number of measuring points in the range. R is window distance, taking 1500 meters.
ΔT Residual of The residual magnetic anomalies at the measurement point are expressed in nT.
And 4.2, loading the obtained residual magnetic anomaly data by using surfer software to obtain a residual anomaly map. The region that lies within the range of the step 1 depressed coverage region and where the residual magnetic anomaly satisfies the following formula is satisfied while the region is delineated as a redox transition zone position S (M).
In the method, in the process of the invention,
the residual anomaly average value of the 20% measuring points with the maximum residual magnetic anomaly value of all measuring points is represented by the unit nT.
The residual anomaly average value of 20% of the measurement points with the minimum residual magnetic anomaly value of all the measurement points is represented by the unit nT.
The average value of residual magnetic anomalies at the measurement point is expressed in nT.
Step 5: and determining an area S (R) which is favorable for ore formation and above the anticline uplift area by utilizing the radon concentration of the soil.
The step 5 comprises the following steps:
step 5.1 soil radon gas working Range requires covering step1, a line distance is 200m, a point distance is 50 m, and the measurement time of each measuring point is 12 minutes, so as to obtain radon concentration R of the measuring point Measuring (i,j)。
Step 5.2, radon concentration data R of each measuring point Measuring (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 a radon anomaly delineation lower limit by using the following formula.
In the method, in the process of the invention,the unit is Bq/m, which is the average value of the smoothed data 3 ；
Sigma is the standard deviation of the smoothed data in Bq/m 3 ；
A (R) is radon concentration anomaly lower limit value, and the unit is Bq/m 3 ；
m and n are the total number of the measuring lines and the number of the measuring points on each measuring line.
And 5.4, loading the data by using surfer software to obtain a radon concentration contour map. The area which is located in the range of the anticline bulge area encircled in the step 1 and has radon concentration value higher than A (R) is encircled as an advantageous mining area S (R).
Step 6: the delineating the ore is beneficial to the ore-forming space.
Intersecting the vertical projection belt of the favorable region S (M) of the depression coverage area outlined in the step 4 with the favorable sand body region S (E) of the audio magnetotelluric sounding outlined in the step 3, wherein the overlapped space body is the favorable mineralization space S (EM) of the depression coverage area; and (3) intersecting the perpendicular projection belt of the back inclined upper favorable region S (R) outlined in the step (5) with the favorable sand body region S (E) outlined in the step (3) by the audio magnetotelluric sounding, wherein the overlapped space body is the favorable ore forming space S (ER) of the region.
S (EM) U S (ER) is the beneficial ore space of uranium ore in the anticline development area.
In the present specification, each embodiment is described in a progressive manner, and each embodiment is mainly described in a different point from other embodiments, and identical and similar parts between the embodiments are all enough to refer to each other. For the system disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant points refer to the description of the method section.
The principles and embodiments of the present invention have been described herein with reference to specific examples, the description of which is intended only to assist in understanding the methods of the present invention and the core ideas thereof; also, it is within the scope of the present invention to be modified by those of ordinary skill in the art in light of the present teachings. In view of the foregoing, this description should not be construed as limiting the invention.
Claims (7)
1. A method for locating a favorable ore-forming space of sandstone-type uranium deposit in a anticline development area, which is characterized by comprising the following steps:
acquiring stratum fluctuation conditions of a anticline development area, and dividing anticline uplift areas and depression coverage areas according to the stratum fluctuation conditions of the anticline development area;
determining resistivity characteristics according to the resistivity log;
inverting the resistivity section contour map according to an audio magnetotelluric sounding method;
determining the space spreading form of the favorable ore forming sand body according to the resistivity section contour map and the resistivity characteristics;
the method for acquiring the residual magnetic anomaly data specifically comprises the following steps:
adopting a formula according to the abnormal data of the chemical pole magnetismAnd DeltaT Residual of ＝ΔT Chemical electrode -ΔT Region(s) Obtaining residual magnetic anomaly data;
wherein DeltaT Region(s) The magnetic abnormality of the region of the measuring point is represented; sum delta T Chemical electrode (0, R) represents the sum of all the pole magnetic anomalies of the measuring point within the range of the radius R of the measuring point; n represents the number of measuring points in the range, R represents the window distance involved in calculation, which is at least 3 times of the distance between the measuring lines, and delta T Residual of Indicating residual magnetic anomalies of the measurement point; delta T Chemical electrode Representing the magnetic anomaly of the chemical pole;
determining a redox transition zone position of the depression coverage area according to the residual magnetic anomaly data, wherein the method specifically comprises the following steps of:
adopting a formula according to the depression coverage area and the residual magnetic anomaly dataDetermining a redox transition zone location of the depression footprint;
wherein,the residual abnormal average value of the 20% measuring points with the maximum residual magnetic abnormal values of all measuring points is represented,residual anomaly average value of 20% measuring points with minimum residual magnetic anomaly value of all measuring points is represented by +.>Mean value of residual magnetic anomalies of measuring points, delta T Residual of Indicating residual magnetic anomalies of the measurement point;
acquiring soil radon concentration data of the anticline uplift area;
determining an advantageous ore forming area above the anticline uplift area according to the soil radon concentration data;
determining a favorable ore-forming space of the depression coverage area according to the redox transition zone position of the depression coverage area and the space distribution form of the favorable ore-forming sand body;
determining a favorable ore forming space of the anticline uplift region according to the favorable ore forming region above the anticline uplift region and the space distribution form of the favorable ore forming sand body;
and determining the beneficial ore-forming space of uranium in the anticline development area according to the beneficial ore-forming space of the depression coverage area and the beneficial ore-forming space of the anticline uplift area.
2. The method for locating a favorable ore formation space of a sandstone uranium deposit in a anticline development area according to claim 1, wherein the determining a space-wise distribution form of favorable ore formation sand according to the resistivity section contour map and the resistivity characteristics specifically includes:
determining the range of each section favorable for the ore formation sand body according to the resistivity section contour map and the resistivity characteristics;
and connecting the ranges of the favorable ore-forming sand bodies of the sections along the shortest paths to obtain the space distribution form favorable for the ore-forming sand bodies.
3. The method for locating the favorable ore formation space of sandstone uranium ores in anticline development area according to claim 1, wherein the determining the favorable ore formation area above the anticline bulge area according to soil radon concentration data specifically comprises the following steps:
adopting a formula according to the soil radon concentration dataAnddetermining an abnormal radon concentration lower limit value;
determining an advantageous ore formation area above the anticline bulge area according to the anticline bulge area and the radon concentration abnormal lower limit value;
wherein,for the average value of the smoothed data, sigma is the standard deviation of the smoothed data, A (R) is the abnormal lower limit value of radon concentration, m and n are the total line number and the number of measuring points on each measuring line, R Measuring (i, j) represents the radon concentration data of the measurement point, R (i, j) represents the radon concentration data R of each measurement point Measuring (i, j) smoothed data, (i, j) respectively represent the line number and the point number of the measurement point.
4. A method of locating a favorable ore formation space of a anticline development zone sandstone uranium deposit according to claim 1, wherein the determining a favorable ore formation space of a depression coverage zone from redox transition zone locations of the depression coverage zone and a spatially spread morphology of the favorable ore formation sand body, specifically comprises:
and intersecting the vertical casting belt at the redox transition belt position of the depression coverage area with the space spreading form of the favorable ore formation sand body to obtain a superposed first space body, wherein the first space body is the favorable ore formation space of the depression coverage area.
5. The method for locating the favorable ore formation space of the sandstone uranium deposit in the anticline development area according to claim 1, wherein the determining the favorable ore formation space of the anticline uplift area according to the favorable ore formation area above the anticline uplift area and the space distribution form of the favorable ore formation sand body specifically comprises:
and intersecting the vertical casting belt of the favorable ore forming area above the anticline uplift area with the space spreading form of the favorable ore forming sand body to obtain a superposed second space body, wherein the second space body is a favorable ore forming space of the anticline uplift area.
6. A method of locating a beneficial mineralization space of a sandstone uranium deposit in a anticline development zone according to claim 1, wherein the determining the beneficial mineralization space of the anticline development zone based on the beneficial mineralization space of the depression coverage zone and the beneficial mineralization space of the anticline bulge zone specifically comprises:
and carrying out union processing on the favorable ore forming space of the depression coverage area and the favorable ore forming space of the anticline uplift area, and determining the favorable ore forming space of uranium ore forming in the anticline development area.
7. A system for locating a favorable ore-forming space of sandstone-type uranium deposit in a anticline development region, comprising:
the anticline development area dividing module is used for acquiring the anticline development area stratum relief condition and dividing anticline bulge areas and depression coverage areas according to the anticline development area stratum relief condition;
the resistivity characteristic determining module is used for determining resistivity characteristics according to the resistivity logging curve;
the resistivity section contour map inversion module is used for inverting the resistivity section contour map according to an audio magnetotelluric sounding method;
the space spreading form determining module is used for determining the space spreading form of the favorable ore-forming sand body according to the resistivity section contour map and the resistivity characteristics;
the residual magnetic abnormal data acquisition module is used for acquiring residual magnetic abnormal data and specifically comprises the following steps:
adopting a formula according to the abnormal data of the chemical pole magnetismAnd DeltaT Residual of ＝ΔT Chemical electrode -ΔT Region(s) Obtaining residual magnetic anomaly data;
wherein DeltaT Region(s) The magnetic abnormality of the region of the measuring point is represented; sum delta T Chemical electrode (0, R) represents the sum of all the pole magnetic anomalies of the measuring point within the range of the radius R of the measuring point; n represents the number of measuring points in the range, R represents the window distance involved in calculation, which is at least 3 times of the distance between the measuring lines, and delta T Residual of Indicating residual magnetic anomalies of the measurement point; delta T Chemical electrode Representing the magnetic anomaly of the chemical pole;
the redox transition zone position determining module of the depression coverage area is used for determining the redox transition zone position of the depression coverage area according to the residual magnetic anomaly data, and specifically comprises the following steps:
adopting a formula according to the depression coverage area and the residual magnetic anomaly dataDetermining a redox transition zone location of the depression footprint;
wherein,the residual abnormal average value of the 20% measuring points with the maximum residual magnetic abnormal values of all measuring points is represented,residual anomaly average value of 20% measuring points with minimum residual magnetic anomaly value of all measuring points is represented by +.>Mean value of residual magnetic anomalies of measuring points, delta T Residual of Indicating residual magnetic anomalies of the measurement point;
the soil radon concentration data acquisition module is used for acquiring the soil radon concentration data of the anticline uplift area;
the beneficial ore forming area determining module is used for determining the beneficial ore forming area above the anticline bulge area according to the soil radon concentration data;
a depression coverage area favorable ore forming space determining module, configured to determine a depression coverage area favorable ore forming space according to a redox transition zone position of the depression coverage area and a spatial distribution form of the favorable ore forming sand body;
the anticline uplift region favorable ore forming space determining module is used for determining the anticline uplift region favorable ore forming space according to the favorable ore forming region above the anticline uplift region and the space distribution form of the favorable ore forming sand body;
and the determination module is used for determining the beneficial ore forming space of the uranium ore forming in the anticline development zone according to the beneficial ore forming space of the depression coverage area and the beneficial ore forming space of the anticline uplift zone.
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 CN113075749A (en) | 2021-07-06 |
| CN113075749B true 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) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114609683B (en) * | 2022-03-15 | 2024-07-12 | 核工业北京地质研究院 | Method and system for determining beneficial region of sandstone uranium ore sand body |
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 |
|---|
| 剩余磁异常的计算与深部磁性盲矿层的预测;王会波;陈载林;;中国锰业(06);全文 * |
| 柴北缘上干柴沟组砂体特征及铀成矿作用研究;廉康;赵兴齐;王继斌;张淼;;东华理工大学学报(自然科学版)(04);全文 * |
| 柴达木盆地西北部可地浸砂岩型铀矿成矿条件及找矿方向;黄国龙, 陈贤春;铀矿地质(06);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113075749A (en) | 2021-07-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Dannemiller et al. | A new method for determination of magnetization direction | |
| CN111856599B (en) | Magnetic measurement data equivalent source pole and type conversion method based on PDE | |
| CN109581513B (en) | Sandstone-type uranium ore mineralization target area delineating method based on multiple spatial scales | |
| CN106650192B (en) | A kind of Volcanic Type Uranium Deposits magnetic interface inversion method | |
| Weller et al. | Three-dimensional inversion of induced polarization data from simulated waste | |
| Ibraheem et al. | Subsurface investigation of the Neogene Mygdonian Basin, Greece using magnetic data | |
| CN106405664A (en) | Magnetic anomaly RTP (Reduction to the Pole) method | |
| CN109884710A (en) | For the micro logging chromatography imaging method of excitation well depth design | |
| CN102236108A (en) | Three-dimensional terrain correcting method for magnetic surface | |
| CN107784159A (en) | A kind of determination method of reservoir resistivity anisotropy coefficient | |
| CN109839670A (en) | A kind of hydrothermal-type uranium deposit substrate interface inversion method | |
| CN105911592B (en) | A kind of real three dimensional seismic data long wavelength's static correcting method for boring constraint | |
| CN113075749B (en) | Method and system for locating favorable ore-forming space of sandstone uranium deposit in anticline development area | |
| Fedorova et al. | Numerical Modeling of the Sources of Magnetic Anomalies in the South Urals Earth’s Crust1 | |
| CN108008456B (en) | A method of delineation mesothermal gold deposits deep three-dimensional emphasis U metallogeny Favourable Target Areas | |
| CN111859251B (en) | Magnetic measurement data equivalent source extension and extension method based on PDE | |
| Muthamilselvan et al. | Spatial confirmation of major lineament and groundwater exploration using ground magnetic method near mecheri village, salem district of Tamil nadu, India | |
| Aali et al. | Geophysical Study to Identify Iron Mineralization Anomalies Using Terrestrial Magnetometry in the Chak-Chak Exploration Area, Iran | |
| CN107797148B (en) | A kind of aeromagnetic anomaly field separation method and system based on three-dimensional geological modeling | |
| Obiora et al. | Interpretation of aeromagnetic data of Idah area in north central Nigeria using combined methods | |
| CN113536693B (en) | Aviation-ground-well magnetic anomaly data joint inversion method based on well rock physical property constraint | |
| CN105259577B (en) | A kind of method and device for the angle information for determining bed boundary | |
| CN105242317B (en) | A kind of determination method and device of velocity of longitudinal wave | |
| CN110618455B (en) | Quantitative evaluation method for sedimentary basin structure uplift | |
| 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 |