CN112462406B - Deep uranium mineralization radioactivity and deep penetration geochemistry combination identification method - Google Patents
Deep uranium mineralization radioactivity and deep penetration geochemistry combination identification method Download PDFInfo
- Publication number
- CN112462406B CN112462406B CN202011253447.1A CN202011253447A CN112462406B CN 112462406 B CN112462406 B CN 112462406B CN 202011253447 A CN202011253447 A CN 202011253447A CN 112462406 B CN112462406 B CN 112462406B
- Authority
- CN
- China
- Prior art keywords
- uranium
- deep
- data
- soil
- variable
- 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
- 229910052770 Uranium Inorganic materials 0.000 title claims abstract description 95
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 title claims abstract description 95
- 238000000034 method Methods 0.000 title claims abstract description 32
- 230000033558 biomineral tissue development Effects 0.000 title claims abstract description 28
- 230000035515 penetration Effects 0.000 title claims abstract description 26
- 239000002689 soil Substances 0.000 claims abstract description 40
- 230000000694 effects Effects 0.000 claims abstract description 39
- 230000009467 reduction Effects 0.000 claims abstract description 19
- 238000012545 processing Methods 0.000 claims abstract description 13
- 238000001914 filtration Methods 0.000 claims abstract description 10
- 230000002349 favourable effect Effects 0.000 claims abstract description 6
- 238000007906 compression Methods 0.000 claims abstract description 5
- 230000006835 compression Effects 0.000 claims abstract description 5
- 239000003153 chemical reaction reagent Substances 0.000 claims description 17
- 238000001228 spectrum Methods 0.000 claims description 16
- 230000009466 transformation Effects 0.000 claims description 16
- CIWBSHSKHKDKBQ-JLAZNSOCSA-N Ascorbic acid Chemical compound OC[C@H](O)[C@H]1OC(=O)C(O)=C1O CIWBSHSKHKDKBQ-JLAZNSOCSA-N 0.000 claims description 14
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 claims description 11
- 229910052802 copper Inorganic materials 0.000 claims description 11
- 239000010949 copper Substances 0.000 claims description 11
- 239000007788 liquid Substances 0.000 claims description 11
- VEXZGXHMUGYJMC-UHFFFAOYSA-N Hydrochloric acid Chemical compound Cl VEXZGXHMUGYJMC-UHFFFAOYSA-N 0.000 claims description 10
- 238000005259 measurement Methods 0.000 claims description 9
- 229960005070 ascorbic acid Drugs 0.000 claims description 7
- 235000010323 ascorbic acid Nutrition 0.000 claims description 7
- 239000011668 ascorbic acid Substances 0.000 claims description 7
- 238000001095 inductively coupled plasma mass spectrometry Methods 0.000 claims description 7
- KCXVZYZYPLLWCC-UHFFFAOYSA-N EDTA Chemical compound OC(=O)CN(CC(O)=O)CCN(CC(O)=O)CC(O)=O KCXVZYZYPLLWCC-UHFFFAOYSA-N 0.000 claims description 5
- FQENQNTWSFEDLI-UHFFFAOYSA-J sodium diphosphate Chemical compound [Na+].[Na+].[Na+].[Na+].[O-]P([O-])(=O)OP([O-])([O-])=O FQENQNTWSFEDLI-UHFFFAOYSA-J 0.000 claims description 5
- 229940048086 sodium pyrophosphate Drugs 0.000 claims description 5
- 235000019818 tetrasodium diphosphate Nutrition 0.000 claims description 5
- 239000001577 tetrasodium phosphonato phosphate Substances 0.000 claims description 5
- YWYZEGXAUVWDED-UHFFFAOYSA-N triammonium citrate Chemical compound [NH4+].[NH4+].[NH4+].[O-]C(=O)CC(O)(CC([O-])=O)C([O-])=O YWYZEGXAUVWDED-UHFFFAOYSA-N 0.000 claims description 5
- 238000010606 normalization Methods 0.000 claims description 4
- 230000010355 oscillation Effects 0.000 claims description 4
- 239000008367 deionised water Substances 0.000 claims description 3
- 229910021641 deionized water Inorganic materials 0.000 claims description 3
- 238000002156 mixing Methods 0.000 claims description 3
- 238000004080 punching Methods 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 238000007873 sieving Methods 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 3
- 238000005406 washing Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 238000001035 drying Methods 0.000 claims description 2
- 238000009472 formulation Methods 0.000 claims description 2
- 239000000203 mixture Substances 0.000 claims description 2
- 238000000691 measurement method Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000012512 characterization method Methods 0.000 description 3
- 238000013508 migration Methods 0.000 description 3
- 230000005012 migration Effects 0.000 description 3
- 230000002285 radioactive effect Effects 0.000 description 3
- 229910052704 radon Inorganic materials 0.000 description 3
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000000904 thermoluminescence Methods 0.000 description 2
- 150000001224 Uranium Chemical class 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 229910052797 bismuth Inorganic materials 0.000 description 1
- JCXGWMGPZLAOME-UHFFFAOYSA-N bismuth atom Chemical compound [Bi] JCXGWMGPZLAOME-UHFFFAOYSA-N 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000011068 loading method Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 229910052699 polonium Inorganic materials 0.000 description 1
- HZEBHPIOVYHPMT-UHFFFAOYSA-N polonium atom Chemical compound [Po] HZEBHPIOVYHPMT-UHFFFAOYSA-N 0.000 description 1
- 229910052705 radium Inorganic materials 0.000 description 1
- HCWPIIXVSYCSAN-UHFFFAOYSA-N radium atom Chemical compound [Ra] HCWPIIXVSYCSAN-UHFFFAOYSA-N 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/167—Measuring radioactive content of objects, e.g. contamination
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/62—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating the ionisation of gases, e.g. aerosols; by investigating electric discharges, e.g. emission of cathode
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/02—Dosimeters
- G01T1/10—Luminescent dosimeters
- G01T1/11—Thermo-luminescent dosimeters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/17—Circuit arrangements not adapted to a particular type of detector
- G01T1/178—Circuit arrangements not adapted to a particular type of detector for measuring specific activity in the presence of other radioactive substances, e.g. natural, in the air or in liquids such as rain water
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/147—Discrete orthonormal transforms, e.g. discrete cosine transform, discrete sine transform, and variations therefrom, e.g. modified discrete cosine transform, integer transforms approximating the discrete cosine transform
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/30—Prediction of properties of chemical compounds, compositions or mixtures
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Chemical & Material Sciences (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computing Systems (AREA)
- Analytical Chemistry (AREA)
- Discrete Mathematics (AREA)
- Biochemistry (AREA)
- Crystallography & Structural Chemistry (AREA)
- Electrochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Algebra (AREA)
- Immunology (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Measurement Of Radiation (AREA)
Abstract
The invention discloses a deep uranium mineralization radioactivity and deep penetration geochemistry combination identification method, which comprises the following steps: determining a coordinate value of a kilometer network; measuring the content of active uranium in soil to be measured, the specific activity of Po-210 and the natural thermoluminescent dosage; the kilometer net coordinate values, the measured active uranium content of the corresponding soil to be measured, the specific activity of Po-210 and the natural pyroelectric light dose are arranged to form two-dimensional space scattered point data; carrying out noise reduction treatment on the active uranium content, the Po-210 specific activity and the natural pyroelectric dose of the soil to be detected to obtain a variable value after noise reduction; normalized compression processing is carried out on the data after the noise reduction processing, and calculation processing is carried out to form gridding dataAnd then carrying out fractal filtering treatment, wherein a data range with a gridding raster data value larger than or equal to 1 after the fractal treatment is the favorable remote scenic region for deep uranium mineralization. The method can effectively avoid the interference of the ground surface background value and deep uranium oresPoor ore finding effect.
Description
Technical Field
The invention belongs to the field of uranium mine exploration, and particularly relates to a deep uranium mineralization radioactivity and deep penetration geochemistry combination identification method.
Background
The uranium element in the deep uranium ore body migrates to the earth surface along with the earth gas under the condition of the difference of the earth pressure, the temperature and the self content, enrichment occurs under the conditions of surface geochemistry barrier and the like, uranium continuously decays in the migration and enrichment process to generate radionuclides such as uranium series daughter radium, radon, polonium, bismuth, lead and the like, in a long-term geological environment, the uranium and the decayed series daughter are in a radioactive equilibrium state, the radioactive dose generated by the decay of the uranium and the daughter can be totally detected through soil natural thermoluminescence measurement, the result of the soil natural thermoluminescence measurement reflects the total quantity of the uranium and the daughter, the daughter of the uranium is the only gas molecule of the uranium decay series, the strong migration characteristic is realized, the radon half-life is short, the migration distance is short in theory, the decay half-life of the daughter of the uranium is relatively long, the daughter Po-210 is an indicator of good deep uranium mineralization, and the specific activity is obtained through soil Po-210 measurement. Meanwhile, a part of nano-active uranium which is not decayed migrates to the surface along the gaps and cracks, and can be obtained through metal multiple dynamic measurement.
However, in the analysis of the actual measurement data, whether it is radiometric or active uranium obtained by deep penetration geochemistry, the result of the indication is not obvious, but is disturbed by the ground background value.
Therefore, development of a new measuring method for deep uranium ores is needed to eliminate surface interference, enhance characterization of deep mineral information and achieve the purpose of deep mining.
Disclosure of Invention
The invention aims to provide a combination identification method of deep uranium mineralization radioactivity and deep penetration geochemistry, which effectively overcomes the defects that the characterization geochemistry information of the deep uranium ores is weak and the exploration mark of a single method is not obvious, and effectively avoids the problems that active uranium obtained by a radioactivity measurement method and a deep penetration geochemistry method is interfered by ground surface background values and the exploration effect of the deep uranium ores is poor.
The technical scheme for realizing the purpose of the invention comprises the following steps: a method of deep uranium mineralization radioactivity and deep penetration geochemistry combination identification, the method comprising the steps of:
step 1, designing deep uranium mine measuring points, determining GPS coordinate points, carrying out site fixed points, punching a plurality of holes at different positions in a range near the site fixed points, and determining kilometer network coordinate values (X, Y) at different depths of each hole;
step 2, respectively collecting fresh soil at the kilometer network coordinate values (X, Y) determined in the step 1, airing, fully mixing, and sieving to obtain soil to be measured for later use;
step 3, preparing a formula reagent of the active uranium, and measuring the content of the active uranium in the soil to be measured through ICP-MS;
step 4, preparing a formula reagent of the Po-210, and measuring the specific activity of the Po-210 of the soil to be measured through a low-background alpha measuring instrument;
step 5, measuring the natural thermoluminescent dosage of the soil to be measured through a thermoluminescent instrument;
step 6, the kilometer network coordinate values (X, Y) determined in the step 1 are arranged with the corresponding active uranium content, po-210 specific activity and natural thermoluminescent dose of the soil to be detected, which are measured in the step 3-5, so as to form two-dimensional space scattered point data;
step 7, carrying out noise reduction treatment on the active uranium content, the Po-210 specific activity and the natural pyroelectric light dose of the soil to be detected which is finished in the step 6, and obtaining a variable value after noise reduction;
and 8, carrying out normalized compression treatment on the data subjected to the noise reduction treatment in the step 7, processing the data into data between 0 and 1, and respectively marking new active uranium content, po-210 specific activity and natural thermoluminescent dosage values as: c (C) U 、C Po 、C TL ；
Step 9, the new active uranium content C formed in the step 8 U Po-210 specific activity C Po Natural thermoluminescent dosage C TL Numerical calculation to form a ratio
Step 10, the ratio formed in the step 9The calculated value of (1) adopts Keli Jin Wangge to interpolate to form gridding data +.>
Step 11, gridding the data formed in step 10And (3) performing fractal filtering processing, drawing a plane contour map, and obtaining a data range with a gridding grid data value of more than or equal to 1 after the fractal processing, namely the favorable remote scenic region for deep uranium mineralization.
Further, the step 7 specifically includes:
step 7.1, converting the active uranium content, po-210 specific activity and natural thermoluminescent dosage data of the soil to be measured which is finished in the step 6 from a two-dimensional space domain into a discrete sine domain by utilizing discrete sine transformation to obtain spectrum data in the discrete sine domain
Step 7.2, in the discrete sine domain, obtaining the average value of 5 times of the spectrum data in the step 7.1, rounding, and determining the spectrum threshold E n ；
Step 7.3, spectral threshold E determined according to step 7.2 n For variable data in discrete sine domainFiltering to make the frequency spectrum in sine domain larger than that in sine domainThreshold E n Is->Assigning zero;
and 7.4, performing sinusoidal inverse transformation according to the variable data subjected to the filtering processing in the step 7.3, and converting the variable data into a two-dimensional space domain from a discrete sinusoidal domain to obtain a variable value subjected to noise reduction.
Further, the discrete sine transformation formula in the step 7.1 is as follows:
wherein j=1, 2,3, respectively represent the active uranium content, po-210 specific activity and natural pyroelectric dose; i is different data of the j variable; n is the number of data of the j variable;the ith original data of the jth variable in the two-dimensional space domain; />Discrete sine transform values of the ith data for the jth variable; a (i) is a transform coefficient.
Further, the calculation formula of the transformation coefficient a (i) is as follows:
further, the formula of the discrete sine inverse transformation in the step 7.4 is as follows:
wherein C' ij The value of the ith measurement point representing the jth variable is subjected to noise reduction.
Further, the formula of the normalization processing in the step 8 is:
wherein C ij The value of the ith measuring point of the jth variable is normalized and compressed, C' j，Min Representing the minimum value of the j-th variable, C' j，Max Represents the maximum value of the j-th variable.
Further, in the step 9The calculation formula of (2) is as follows:
wherein C is U，i 、C Po，i 、G TL，i Respectively represent the content C of active uranium U Po-210 specific activity C Po Natural thermoluminescent dosage C TL The i-th value of the number,representing the corresponding ratio.
Further, the step 3 specifically includes: preparing a formula reagent of the active uranium, injecting the formula reagent of the active uranium into soil to be detected to form mixed liquid, vibrating the mixed liquid, performing solid-liquid separation, and performing ICP-MS (inductively coupled plasma-mass spectrometry) measurement on the separated liquid to obtain the content of the active uranium.
Further, the formula reagent of the active uranium in the step 3 consists of 40g/L ammonium citrate solution, 30g/L sodium pyrophosphate solution and 30g/L ethylenediamine tetraacetic acid solution, wherein the volume ratio of the ammonium citrate solution to the sodium pyrophosphate solution to the ethylenediamine tetraacetic acid solution is 2:1:1.
further, the step 4 specifically includes: sequentially placing soil to be measured, a formula reagent of Po-210 and a copper sheet into a beaker for oscillation, taking out the copper sheet and deionized water after oscillation, washing the copper sheet, drying and standing, and measuring by adopting a low-background alpha measuring instrument to obtain the specific activity of Po-210.
Further, in the step 4, the formulation reagent of Po-210 consists of 2mol/L hydrochloric acid and ascorbic acid, wherein the volume mass ratio of the hydrochloric acid to the ascorbic acid is 10:1, and the mass ratio of the ascorbic acid to the soil to be detected is 2:1.
The beneficial technical effects of the invention are as follows:
1. the invention relates to a combination identification method of deep uranium mineralization radioactivity and deep penetration geochemistry, which comprises the steps of obtaining deep nano-to micron-level deep penetration geochemistry uranium content by adopting a prepared extraction reagent, obtaining deep uranium decay daughter Po-210 specific activity by adopting a radon daughter Po-210 measurement method, obtaining radioactive natural pyroelectric dosage by adopting a natural pyroelectric measurement method, compressing three data to be between 0 and 1, and dividing a favorable prediction area of uranium ore mineralization by utilizing a data range with a gridding raster data value of more than or equal to 1 as a favorable far-field area of deep uranium mineralization;
2. the deep uranium mineralization radioactivity and deep penetration geochemistry combination identification method effectively overcomes the defects that the deep uranium ore characterization geochemistry information is weak and a single method is not obvious in prospecting mark;
3. the combination identification method for deep uranium mineralization radioactivity and deep penetration geochemistry effectively avoids the problems that active uranium obtained by a radioactivity measurement method and a deep penetration geochemistry method is interfered by ground surface background values and the prospecting effect of deep uranium ores is poor.
Drawings
Fig. 1 is an abnormal synthesis diagram obtained by a combination of deep uranium mineralization radioactivity and deep penetration geochemistry identification method according to the present invention.
Detailed Description
The invention is described in further detail below with reference to the drawings and examples.
The invention discloses a deep uranium mineralization radioactivity and deep penetration geochemistry combination identification method, which comprises the following steps:
step 1, selecting a certain mining area and a peripheral research area, and designing measuring points according to the point-line distance of 100m multiplied by 500m to form GPS coordinate points; according to the determined GPS coordinate points, performing site fixed points and marking by using a GPS positioning instrument, punching 3 holes at different positions of a 1m range of a square circle at the position of the designed coordinate point of each measuring point, respectively taking different depths of 40cm, 50cm and 60cm from each hole, and determining coordinate values (X, Y) of a kilometer network;
step 2, respectively collecting 300g of fresh soil at kilometer network coordinate values (X, Y) determined in the step 1, loading a total of 2.7kg into a sample cloth bag, marking the sample cloth bag, airing a soil sample, fully mixing, sieving the soil sample with a 200-mesh sieve, and taking 2g, 5g and 1g of each measuring point soil sample as soil to be measured for subsequent analysis and test of active uranium content, po-210 specific activity and natural pyroelectric dose;
step 3, selecting 40g/L ammonium citrate solution, 30g/L sodium pyrophosphate solution and 30g/L ethylenediamine tetraacetic acid solution according to the volume ratio of 2:1:1 preparing a formula reagent for forming active uranium, putting 1g of soil to be detected, which is obtained in the step 2, into a test tube, injecting 30mL of the formula reagent to form mixed liquid fused with the soil to be detected, fixing the test tube on a shaking table, shaking for 4 hours, placing the test tube on a centrifuge for solid-liquid separation, taking out the liquid, and measuring by ICP-MS to obtain the content (unit: ng/g) of active uranium;
step 4, putting 5g of the soil to be detected, which is obtained in the step 2, into a beaker, sequentially injecting 100mL of 2mol/L hydrochloric acid, 10g of ascorbic acid and 1.5cm diameter round copper flakes into the beaker, then fixing the beaker on a constant-temperature shaking table at 85 ℃ for shaking for 2 hours, taking out the copper flakes by using plastic tweezers, washing the copper flakes with deionized water, airing the copper flakes for 4 hours, and measuring the copper flakes by using a special low-background alpha measuring instrument to obtain the specific activity (unit: bq/kg);
step 5, measuring 2g of the soil to be measured obtained in the step 2 on a thermoluminescent instrument to obtain natural thermoluminescent dose (unit: mu Gy);
step 6, the kilometer network coordinate values (X, Y) determined in the step 1, the corresponding active uranium content, po-210 specific activity and natural thermoluminescent dosage of the soil to be detected, which are measured in the step 3-5, are arranged into EXCEL format files to form two-dimensional space scattered point data;
step 7, carrying out noise reduction treatment on the active uranium content, the Po-210 specific activity and the natural pyroelectric light dose of the soil to be detected which is finished in the step 6, and obtaining a variable value after noise reduction;
step 7.1, converting the active uranium content, po-210 specific activity and natural thermoluminescent dosage data of the soil to be measured which is finished in the step 6 from a two-dimensional space domain into a discrete sine domain by utilizing discrete sine transformation to obtain spectrum data in the discrete sine domainThe formula of the discrete sine transformation is:
wherein j=1, 2,3, respectively represent the active uranium content, po-210 specific activity and natural pyroelectric dose; i is different data of the j variable; n is the number of data of the j variable;the ith original data of the jth variable in the two-dimensional space domain; />Discrete sine transform values of the ith data for the jth variable; a (i) is a transformation coefficient, and a calculation formula is as follows:
wherein N is the number of data for the j variable;
step 7.2, in the discrete sine domain, obtaining the average value of 5 times of the spectrum data in the step 7.1, rounding, and determining the spectrum threshold E n ；
Step 7.3, spectral threshold E determined according to step 7.2 n For variable data in discrete sine domainFiltering to obtain a sinusoidal spectrum with a frequency spectrum greater than the frequency spectrum threshold E n Is>Assigning zero;
and 7.4, performing sinusoidal inverse transformation according to the variable data subjected to the filtering processing in the step 7.3, and converting the variable data into a two-dimensional space domain from a discrete sinusoidal domain to obtain a variable value subjected to noise reduction, wherein the formula of the discrete sinusoidal inverse transformation is as follows:
wherein j=1, 2,3, respectively represent the active uranium content, po-210 specific activity and natural pyroelectric dose; i is different data of the j variable; n is the number of data of the j variable;discrete sine transform values of the ith data for the jth variable; c'. ij The value of the ith measuring point of the jth variable is subjected to noise reduction treatment; a (i) is a transformation coefficient, and a calculation formula is as follows:
wherein N is the number of data for the j variable;
and 8, carrying out normalized compression treatment on the active uranium content, the Po-210 specific activity and the natural pyroelectric dosage value subjected to noise reduction treatment in the step 7 according to a normalization treatment formula, and processing the normalized compression treatment to obtain data between 0 and 1, wherein the new active uranium content, the Po-210 specific activity and the natural pyroelectric dosage value are respectively recorded as follows: c (C) U 、C Po 、C TL The formula of the normalization process is:
wherein j=1, 2,3, respectively represent the active uranium content, po-210 specific activity and natural pyroelectric dose; i is different data of the j variable; c'. ij The value of the ith measuring point of the jth variable is subjected to noise reduction treatment; c' ij The value of the ith measuring point of the jth variable is normalized and compressed, C' j，Min Representing the minimum value of the j-th variable, C' j，Max Represents the maximum value of the j-th variable;
step 9, normalizing the active uranium content C obtained in the step 8 U Po-210 specific activity C Po Natural thermoluminescent dosage C TL The numerical value of (2) is calculated to form a ratio according to the following formula
Wherein C is U，i 、C Po，i 、C TL，i Respectively represent the content C of active uranium U Po-210 specific activity C Po Natural thermoluminescent dosage C TL Is a function of the i-th value of (c),representing the corresponding ratio;
step 10, the ratio formed in the step 9The calculated value of (1) adopts Keli Jin Wangge to interpolate to form gridding data +.>
Step 11, gridding the data formed in step 10Fractal filtering and paintingAnd (3) preparing a plane contour map, wherein a data range with a gridding raster data value being more than or equal to 1 after fractal processing is the favorable remote scenic region for deep uranium mineralization.
The abnormal comprehensive diagram obtained by the combination identification method of deep uranium mineralization radioactivity and deep penetration geochemistry is shown in figure 1.
The present invention has been described in detail with reference to the drawings and the embodiments, but the present invention is not limited to the embodiments described above, and various changes can be made within the knowledge of those skilled in the art without departing from the spirit of the present invention. The invention may be practiced otherwise than as specifically described.
Claims (10)
1. A method for identifying a combination of deep uranium mineralization radioactivity and deep penetration geochemistry, the method comprising the steps of:
step 1, designing deep uranium mine measuring points, determining GPS coordinate points, carrying out site fixed points, punching a plurality of holes at different positions in a range near the site fixed points, taking a plurality of different depths from each hole, and determining kilometer network coordinate values;
step 2, respectively collecting fresh soil at the kilometer network coordinate values determined in the step 1, airing, fully mixing, and sieving to obtain soil to be measured for later use;
step 3, preparing a formula reagent of the active uranium, and measuring the content of the active uranium in the soil to be measured through ICP-MS;
step 4, preparing a formula reagent of the Po-210, and measuring the specific activity of the Po-210 of the soil to be measured through a low-background alpha measuring instrument;
step 5, measuring the natural thermoluminescent dosage of the soil to be measured through a thermoluminescent instrument;
step 6, finishing the kilometer network coordinate value determined in the step 1, the active uranium content, the Po-210 specific activity and the natural pyroelectric light dose of the corresponding soil to be detected, which are determined in the step 3-5, to form two-dimensional space scattered point data;
step 7, carrying out noise reduction treatment on the active uranium content, the Po-210 specific activity and the natural pyroelectric light dose of the soil to be detected which is finished in the step 6, and obtaining a variable value after noise reduction;
and 8, carrying out normalized compression treatment on the data subjected to the noise reduction treatment in the step 7, processing the data into data between 0 and 1, and respectively marking new active uranium content, po-210 specific activity and natural thermoluminescent dosage values as: c (C) U 、C Po 、C TL ；
Step 9, the new active uranium content C formed in the step 8 U Po-210 specific activity C Po Natural thermoluminescent dosage C TL Numerical calculation to form a ratio
Step 10, the ratio formed in the step 9The calculated value of (1) adopts Keli Jin Wangge to interpolate to form gridding data
Step 11, gridding the data formed in step 10Performing fractal filtering treatment, and drawing a planar contour map, wherein a data range with a gridding grid data value being greater than or equal to 1 after the fractal treatment is a favorable remote scenic region for deep uranium mineralization;
in the step 9The calculation formula of (2) is as follows:
wherein C is U，i 、C Po，i 、C TL，i Respectively represent the content C of active uranium U 、Po-210 specific activity C Po Natural thermoluminescent dosage C TL The i-th value of the number,representing the corresponding ratio.
2. The method for identifying a combination of deep uranium mineralization radioactivity and deep penetration geochemistry according to claim 1, wherein the step 7 specifically includes:
step 7.1, converting the active uranium content, po-210 specific activity and natural thermoluminescent dosage data of the soil to be measured which is finished in the step 6 from a two-dimensional space domain into a discrete sine domain by utilizing discrete sine transformation to obtain spectrum data in the discrete sine domain
Step 7.2, in the discrete sine domain, obtaining the average value of 5 times of the spectrum data in the step 7.1, rounding, and determining the spectrum threshold E n ；
Step 7.3, spectral threshold E determined according to step 7.2 n For variable data in discrete sine domainFiltering to obtain a sinusoidal spectrum with a frequency spectrum greater than the frequency spectrum threshold E n Is>Assigning zero;
and 7.4, performing sinusoidal inverse transformation according to the variable data subjected to the filtering processing in the step 7.3, and converting the variable data into a two-dimensional space domain from a discrete sinusoidal domain to obtain a variable value subjected to noise reduction.
3. The method for identifying deep uranium mineralization radioactivity and deep penetration geochemistry according to claim 2, wherein the formula of the discrete sine transformation in the step 7.1 is as follows:
wherein j=1, 2,3, respectively represent the active uranium content, po-210 specific activity and natural pyroelectric dose; i is different data of the j variable; n is the number of data of the j variable;the ith original data of the jth variable in the two-dimensional space domain; />Discrete sine transform values of the ith data for the jth variable; a (i) is a transform coefficient.
4. A combination deep uranium mineralization and deep penetration geochemical identification method according to claim 3, wherein the transformation coefficient a (i) is calculated according to the formula:
5. the method for identifying deep uranium mineralization radioactivity and deep penetration geochemical combination according to claim 4, wherein the formula of discrete sine inverse transformation in the step 7.4 is as follows:
wherein C' ij The value of the ith measurement point representing the jth variable is subjected to noise reduction.
6. The method for identifying a combination of deep uranium mineralization radioactivity and deep penetration geochemistry according to claim 5, wherein the formula of normalization in step 8 is as follows:
wherein C ij The value of the ith measuring point of the jth variable is normalized and compressed, C' j，Min Representing the minimum value of the j-th variable, C' j，Max Represents the maximum value of the j-th variable.
7. The method for identifying deep uranium mineralization radioactivity and deep penetration geochemistry according to claim 1, wherein the step 3 specifically includes: preparing a formula reagent of the active uranium, injecting the formula reagent of the active uranium into soil to be detected to form mixed liquid, vibrating the mixed liquid, performing solid-liquid separation, and performing ICP-MS (inductively coupled plasma-mass spectrometry) measurement on the separated liquid to obtain the content of the active uranium.
8. The method for identifying deep uranium mineralized radioactivity and deep penetration geochemical combination according to claim 7, wherein the formula reagent of the active uranium in the step 3 is composed of 40g/L ammonium citrate solution, 30g/L sodium pyrophosphate solution and 30g/L ethylenediamine tetraacetic acid solution, and the volume ratio of the ammonium citrate solution, the sodium pyrophosphate solution and the ethylenediamine tetraacetic acid solution is 2:1:1.
9. the method for identifying a combination of deep uranium mineralization radioactivity and deep penetration geochemistry according to claim 1, wherein the step 4 specifically includes: sequentially placing soil to be measured, a formula reagent of Po-210 and a copper sheet into a beaker for oscillation, taking out the copper sheet and deionized water after oscillation, washing the copper sheet, drying and standing, and measuring by adopting a low-background alpha measuring instrument to obtain the specific activity of Po-210.
10. The method for identifying deep uranium mineralization radioactivity and deep penetration geochemistry according to claim 9, wherein in the step 4, the formulation reagent of Po-210 consists of 2mol/L hydrochloric acid and ascorbic acid, the volume-mass ratio of hydrochloric acid to ascorbic acid is 10:1, and the mass ratio of ascorbic acid to soil to be detected is 2:1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011253447.1A CN112462406B (en) | 2020-11-11 | 2020-11-11 | Deep uranium mineralization radioactivity and deep penetration geochemistry combination identification method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011253447.1A CN112462406B (en) | 2020-11-11 | 2020-11-11 | Deep uranium mineralization radioactivity and deep penetration geochemistry combination identification method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112462406A CN112462406A (en) | 2021-03-09 |
| CN112462406B true CN112462406B (en) | 2024-03-19 |
Family
ID=74826579
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011253447.1A Active CN112462406B (en) | 2020-11-11 | 2020-11-11 | Deep uranium mineralization radioactivity and deep penetration geochemistry combination identification method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112462406B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5235190A (en) * | 1991-01-29 | 1993-08-10 | Gemini Research, Inc. | Continuous air monitoring system |
| CN101620274A (en) * | 2009-08-12 | 2010-01-06 | 成都理工大学 | Alpha energy spectrum measuring method of radon in soil based on static diffusion and electrostatic adsorption principle of radon in soil |
| CN102478674A (en) * | 2010-11-25 | 2012-05-30 | 核工业北京地质研究院 | Method for combined tracing exploration on hydrothermal uranium ore through geochemical elements |
| CN105738937A (en) * | 2016-04-29 | 2016-07-06 | 成都理工大学 | Static alpha-cup multifunctional emanometer |
| CN110908010A (en) * | 2019-10-25 | 2020-03-24 | 东华理工大学 | Effective geophysical method for finding shallow sandstone type uranium deposit within 800 meters |
| CN111142167A (en) * | 2019-11-06 | 2020-05-12 | 中国辐射防护研究院 | Method for detecting position of spontaneous combustion source of coal by using polonium-210 method |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE602004009184T2 (en) * | 2003-04-30 | 2008-07-03 | Ramot At Tel Aviv University Ltd. | METHOD AND DEVICE FOR RADIOTHERAPY |
| US7973287B2 (en) * | 2007-08-14 | 2011-07-05 | Nucsafe, Inc | Detection of alpha-particle emitting isotopes |
-
2020
- 2020-11-11 CN CN202011253447.1A patent/CN112462406B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5235190A (en) * | 1991-01-29 | 1993-08-10 | Gemini Research, Inc. | Continuous air monitoring system |
| CN101620274A (en) * | 2009-08-12 | 2010-01-06 | 成都理工大学 | Alpha energy spectrum measuring method of radon in soil based on static diffusion and electrostatic adsorption principle of radon in soil |
| CN102478674A (en) * | 2010-11-25 | 2012-05-30 | 核工业北京地质研究院 | Method for combined tracing exploration on hydrothermal uranium ore through geochemical elements |
| CN105738937A (en) * | 2016-04-29 | 2016-07-06 | 成都理工大学 | Static alpha-cup multifunctional emanometer |
| CN110908010A (en) * | 2019-10-25 | 2020-03-24 | 东华理工大学 | Effective geophysical method for finding shallow sandstone type uranium deposit within 800 meters |
| CN111142167A (en) * | 2019-11-06 | 2020-05-12 | 中国辐射防护研究院 | Method for detecting position of spontaneous combustion source of coal by using polonium-210 method |
Non-Patent Citations (4)
| Title |
|---|
| 二连盆地第四系沉积物天然热释光特征及其找矿意义;王南萍, 侯胜利, 刘海生, 李怀渊, 程业勋;物探与化探(04);全文 * |
| 分量化探法在铀资源勘查中的应用;葛祥坤;尹金双;范光;艾永亮;张建锋;金贵善;张彦辉;张良圣;;铀矿地质(01);全文 * |
| 学堂坳地区铀矿物探方法找矿应用效果探讨;陈文财;;西部资源(06);全文 * |
| 鹿井地区铀矿化放射性信息提取与应用;符志军;杨亚新;罗齐彬;王帅帅;洪昆;李星阳;;江西科学(02);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112462406A (en) | 2021-03-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103824133B (en) | A kind of granite type U-ore field prospective area Comprehensive Prediction Method | |
| CN111679342B (en) | Method for quickly searching basin hidden sandstone type uranium ores | |
| CN203376483U (en) | Fine gamma ray full-spectrum fast radon-measurement apparatus | |
| Duggal et al. | Measurement of soil-gas radon in some areas of northern Rajasthan, India | |
| CN107144868A (en) | A kind of method that live multi-parameter fast investigation evaluates coalfield radioactive radiation environment | |
| CN103076292B (en) | System and method for quantitatively evaluating concentration of rare earth in water body based on ground reflection spectrum | |
| CN103678778A (en) | Method for radioactive geophysical and geochemical exploration information integration | |
| Oliver et al. | A geostatistical investigation of the spatial variation of radon in soil | |
| Ramola et al. | Radon occurrence in soil–gas and groundwater around an active landslide | |
| Chiozzi et al. | Ground radiometric survey of U, Th and K on the Lipari Island, Italy | |
| Buccianti et al. | Natural radioactivity levels (K, Th, U and Rn) in the Cecita Lake area (Sila Massif, Calabria, Southern Italy): An attempt to discover correlations with soil features on a statistical base | |
| Gaafar et al. | Analysis of geological and geophysical datasets for radioelement exploration in Kab Amiri area, Central Eastern Desert, Egypt | |
| CN112462406B (en) | Deep uranium mineralization radioactivity and deep penetration geochemistry combination identification method | |
| CN107664765A (en) | A kind of radon daughter board accumulates absorption type Soil Radon rapid measurement device | |
| CN107976718B (en) | A kind of deep sandstone type uranium mineralization direct information exploration method | |
| CN110095161A (en) | A kind of acquisition system and its method for geological environment exploration | |
| Morgenstern et al. | 32Si in limestone aquifers | |
| CN111983721A (en) | Sandstone uranium ore mud-sand-mud geological structure identification method and system | |
| Oliver et al. | Investigating the spatial variation of radon in soil geostatistically | |
| Robinson et al. | Portable gamma ray spectrometry for archaeological prospection: A preliminary investigation at Silchester Roman Town | |
| CN113406711B (en) | Method for circling uranium-bearing ore-forming beneficial zone based on radon-active uranium comprehensive index | |
| CN113640498B (en) | Effective combination method for verifying navigation weak information of coverage area | |
| CN109901211A (en) | A kind of fine gamma based on lanthanum bromide crystal is composed entirely surveys radon method | |
| Thomas | Field Tests of a UAV-compatible Spectrometer to Evaluate Its Suitability for Detailed Soil Radon Potential Mapping | |
| Asfahani et al. | Identifying spectrometric signatures of phosphate deposits and enclosing sediments in Al-Awabed area, Northern Palmyrides, Central Syria, by the use of statistical factor analysis |
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 |