CN111487309A - Uranium ore geographical tracing method based on discriminant analysis - Google Patents
Uranium ore geographical tracing method based on discriminant analysis Download PDFInfo
- Publication number
- CN111487309A CN111487309A CN202010352005.6A CN202010352005A CN111487309A CN 111487309 A CN111487309 A CN 111487309A CN 202010352005 A CN202010352005 A CN 202010352005A CN 111487309 A CN111487309 A CN 111487309A
- Authority
- CN
- China
- Prior art keywords
- uranium ore
- sample
- geographical
- discriminant analysis
- method based
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- 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
- G01N27/626—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 using heat to ionise a gas
Landscapes
- Chemical & Material Sciences (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
The invention discloses a uranium ore geographical tracing method based on discriminant analysis. The method comprises the following steps: (1) collecting a uranium ore sample, and carrying out pretreatment to obtain a uranium ore sample solution; (2) measuring the content of rare earth elements in a uranium ore sample by adopting an inductively coupled plasma mass spectrometer; (3) analyzing the Ce isotope ratio in the uranium ore sample by adopting a thermal surface ionization mass spectrometer; (4) and (4) performing discriminant analysis on the content of the rare earth elements and the Ce isotope ratio, and judging the source area of the uranium ore sample according to the calculation result. The judgment method provided by the invention has the advantages that the judgment accuracy rate of the uranium ore source can reach 85%, the geographical tracing effect is ideal, and the method has important guiding significance for quickly realizing tracing of nuclear materials.
Description
Technical Field
The invention relates to the technical field of ore traceability, in particular to a uranium ore geographical traceability method based on discriminant analysis.
Background
In recent years, a comprehensive discipline involving multiple disciplines, the nuclear evidence collection, has been developed. The purpose of nuclear forensics is to provide characteristic quantities and process history information of intercepted nuclear materials as much as possible, provide evidences for tracing the origin, production process and transportation route of the materials, and provide technical support for proving the responsibility of exploring illegally owning the nuclear materials and radioactive materials. The tracing analysis of the geographic source of the nuclear material is the key point of the nuclear evidence research, and the uranium ore serving as the basic raw material of the nuclear material becomes an important analysis object of the nuclear evidence research.
In nuclear evidence research, a characteristic fingerprint that provides accurate and reliable geographic traceability information is generally referred to as a regional indicator. At present, the traditional uranium ore traceability analysis method mainly adopts trace element content or stable isotope composition as a region indicator for research, however, if the stable isotope composition is independently used as the region indicator, the defects of extremely limited contained information and low accuracy rate exist; if the trace elements are independently used as the regional indicators, the geographical tracing is usually performed by comparing distribution curves of the trace elements, and due to the fact that the trace elements are numerous in types and have large content differences, the geographical tracing is difficult to directly perform geographical tracing research by adopting the method, the subjectivity is strong, the quantification is difficult, and the error is large. In addition, because the chemical properties of the rare earth elements in the trace elements are relatively close and the correlation coefficient is relatively large, namely, the contained information is relatively overlapped, the problems of relatively large calculated amount and low analysis efficiency exist in the traditional method for performing nuclear material traceability analysis by taking the rare earth elements as a region indicator.
Disclosure of Invention
In view of this, the invention aims to provide an effective and high-accuracy uranium ore geographical tracing method based on discriminant analysis.
In order to achieve the above purpose, the technical scheme of the invention is as follows:
A uranium ore geographical tracing method based on discriminant analysis is characterized by comprising the following steps:
(1) Collecting m uranium ore samples, and carrying out pretreatment to obtain uranium ore sample solutions;
(2) Respectively measuring the content of rare earth elements in m uranium ore samples by adopting an inductively coupled plasma mass spectrometer;
(3) Respectively measuring Ce isotope ratio values in m uranium ore samples by adopting a thermal surface ionization mass spectrometer;
(4) And performing discriminant analysis according to the rare earth element content and the Ce isotope ratio of each uranium ore sample, and classifying the uranium ore samples according to the production place sources. Wherein, the step (4) specifically comprises the following steps:
Randomly taking 1 sample from m uranium ore samples as a sample to be distinguished, and dividing the rest (m-1) uranium ore samples into n types according to sources;
Step (4.2) determining the independent variables of the discriminant function and the discriminant function coefficients corresponding to the independent variables, and establishing corresponding discriminant functions for uranium ore samples of various source areas;
Step (4.3) substituting the rare earth element content and the Ce isotope ratio of the sample to be judged into each judgment function, and calculating to obtain a judgment function value of each judgment model;
And (4.4) comparing the discrimination function values, and determining the source corresponding to the maximum discrimination function value as the source of the sample to be discriminated.
Further, the step (4.2) specifically comprises the step (4.2.1): firstly, performing variable defect tolerance detection on the content of all rare earth elements and the Ce isotope ratio of the (m-1) uranium ore samples, eliminating the content of corresponding rare earth elements and/or the Ce isotope ratio with the tolerance less than 0.1, and taking the content of the rest rare earth elements and/or the Ce isotope ratio which are not eliminated as the independent variable of a discriminant function.
Further, the step (4.2) further comprises the step (4.2.2): the method for determining the coefficient of the discrimination function corresponding to each variable maximizes the ratio of the difference between various types and the difference in various types of the linear function composed of independent variables through variance analysis.
Further, the step (1) specifically comprises the following steps:
Step (1.1) grinding a uranium ore sample into powder with uniform particles by using a ball grinder;
Weighing a uranium ore powder sample and placing the uranium ore powder sample in a polytetrafluoroethylene beaker;
Digesting: sequentially adding HNO into a polytetrafluoroethylene beaker 3HF and HClO 4Heating the beaker on an electric heating plate for a certain time;
Step (1.4) after the sample solution is completely evaporated to dryness, utilizing concentrated HNO 3Dissolving;
Step (1.5) repeating operation step (1.4) for multiple times until the dissolution is complete;
Step (1.6) after the sample solution is evaporated to dryness again, HCl and H are used 3BO3And re-dissolving the mixed acid to obtain a uranium ore sample solution.
Further, in the step (1.3), HNO 3the concentration is 15 mol/L;
further, in the step (1.3), the concentration of HF is 20 mol/L;
Further, in the step (1.3), HClO 4the concentration is 12.4 mol/L;
Further, in the step (1.6), H 3BO3the concentration is 0.5 mol/L;
further, the inductively coupled plasma mass spectrometry in the step (2) is performed under the conditions that the sample introduction speed is 0.3rps, the sample is stable for 35 seconds before analysis, the He gas mode is adopted, the collection point number of the collection number per unit mass number is 3, the data collection repetition time is 3 times, the radio frequency power is 1550W, the carrier gas flow rate is 5.5m L/min, and the atomization chamber temperature is 2 ℃.
Further, the step (3) of measuring the Ce isotope ratio in the uranium ore sample by using a thermal surface ionization mass spectrometer specifically comprises the following steps:
Step (3.1) coating a sample solution on the center of a rhenium strip filament by using a micropipette;
Step (3.2) the filament after sample coating is loaded into an ion source, the filament is heated, a central Faraday cup is adopted to receive ions to be detected, and the ion current peak shape and the peak center are adjusted;
Step (3.3) adjusting the focusing parameters of the electron lens of the thermal surface ionization mass spectrometer, and heating the filament at a low speed to obtain a stable and strong filament 138Ce+And 140Ce+Ion current intensity;
Step (3.4) when the vacuum degree of the ion source reaches a certain value, obtaining 138Ce/140And obtaining the Ce isotope ratio.
further, the step (3.4) is specifically that the vacuum degree of the ion source is more than or equal to 1.1 × 10 -7mbar, is obtained 138Ce/140The ratio of Ce.
Compared with the prior art, the invention has the following advantages and beneficial effects:
1. In the aspect of selection of the traceability indicator, the invention adopts the content of the rare earth element and the Ce isotope ratio as the geographical traceability indicator of the uranium ore, and solves the problems that other common traceability indicator elements are easily influenced by the mining and metallurgy process or have poor synchronism and the like.
2. The method utilizes a discriminant analysis method to carry out geographic tracing research on the uranium ores, solves the problems that the existing visual comparison method is difficult to quantify and has obvious errors, can realize rapid and accurate nuclear geographic tracing, and provides a scientific method for establishing and perfecting a nuclear material tracing model.
3. The uranium ore geographical tracing method provided by the invention realizes accurate geographical tracing of the uranium ore by combining the rare earth element content and the Ce isotope composition, and can be used for tracing analysis of nuclear materials, radioactive materials and other related materials.
Detailed Description
The present invention will now be described more fully hereinafter with reference to the accompanying examples, in which some, but not all embodiments of the invention are shown. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
A uranium ore geographical tracing method based on discriminant analysis is characterized by comprising the following steps:
(1) Collecting m uranium ore samples from n source places, wherein the number of the uranium ore samples in each source place is not less than 3, and preprocessing the uranium ore samples to obtain uranium ore sample solutions;
(2) Respectively measuring the content of rare earth elements in m uranium ore samples by adopting an inductively coupled plasma mass spectrometer;
(3) Respectively measuring Ce isotope ratio values in m uranium ore samples by adopting a thermal surface ionization mass spectrometer;
(4) And performing discriminant analysis according to the rare earth element content and the Ce isotope ratio of each uranium ore sample, and judging the source area of the uranium ore sample according to the calculation result.
Further, the step (4) specifically comprises the following steps:
Randomly taking 1 sample from m uranium ore samples as a sample to be distinguished, and dividing the rest (m-1) uranium ore samples into n types according to sources;
Step (4.2) determining the independent variables of the discriminant function and the discriminant function coefficients corresponding to the independent variables, and establishing corresponding discriminant functions for uranium ore samples of various source areas;
Step (4.3) substituting the rare earth element content and the Ce isotope ratio of the sample to be judged into each judgment function, and calculating to obtain a judgment function value of each judgment model;
And (4.4) comparing the discrimination function values, and determining the source corresponding to the maximum discrimination function value as the source of the sample to be discriminated.
Further, step (4.2) comprises step (4.2.1): firstly, performing variable defect tolerance inspection on the content of all rare earth elements and the Ce isotope ratio of the (m-1) uranium ore samples, eliminating the content of corresponding rare earth elements and/or the Ce isotope ratio with the tolerance less than 0.1, and taking the content of the rest rare earth elements and/or the Ce isotope ratio which are not eliminated as the independent variable of a discriminant function.
Further, the step (4.2) further comprises the step (4.2.2): the method for determining the coefficient of the discrimination function corresponding to each variable maximizes the ratio of the difference between various types and the difference in various types of the linear function composed of independent variables through variance analysis. Further, the step (1) specifically comprises the following steps:
Step (1.1) grinding a uranium ore sample into a powder with uniform particles, the particle size being about 10 microns, by using a ball grinder;
Step (1.2) weighing a uranium ore powder sample with a certain mass by using a high-precision balance and placing the uranium ore powder sample into a polytetrafluoroethylene beaker;
Digesting: sequentially adding HNO into a polytetrafluoroethylene beaker 3HF and HClO 4Heating the beaker on an electric heating plate for a certain time;
Step (1.4) after the sample solution is completely evaporated to dryness, utilizing concentrated HNO 3Dissolving;
Step (1.5) repeating operation step (1.4) for multiple times until the dissolution is complete;
Step (1.6) after the sample solution is evaporated to dryness again, HCl and H are used 3BO3The mixed acid is dissolved again to eliminate fluorinated complexes possibly generated in the digestion process of the uranium ore sample, and the uranium ore sample solution is obtained.
Further, in the step (1.3), HNO 3the concentration is 15 mol/L;
further, in the step (1.3), the concentration of HF is 20 mol/L;
Further, in the step (1.3), HClO 4the concentration is 12.4 mol/L;
Further, in the step (1.6), H 3BO3the concentration is 0.5 mol/L;
further, the inductively coupled plasma mass spectrometry in the step (2) is performed under the conditions that the sample introduction speed is 0.3rps, the sample is stable for 35 seconds before analysis, the He gas mode is adopted, the collection point number of the collection number per unit mass number is 3, the data collection repetition time is 3 times, the radio frequency power is 1550W, the carrier gas flow rate is 5.5m L/min, and the atomization chamber temperature is 2 ℃.
Further, the step (3) of measuring the Ce isotope ratio in the uranium ore sample by using a thermal surface ionization mass spectrometer specifically comprises the following steps:
Step (3.1) coating the sample solution on the center of a rhenium strip filament by using a micropipette, which is beneficial to improving the ionization efficiency of the sample;
Step (3.2) the filament after sample coating is loaded into an ion source, the filament is heated, a central Faraday cup is adopted to receive ions to be detected, and the ion current peak shape and the peak center are adjusted;
Step (3.3) adjusting the focusing parameters of the electron lens of the thermal surface ionization mass spectrometer, and heating the filament at a low speed to obtain a stable and strong filament 138Ce+And 140Ce+Ion current intensity;
Step (3.4) when the vacuum degree of the ion source reaches a certain value, adopting a Faraday cup to receive the ion source simultaneously 138Ce+And 140Ce+Ion flow to obtain 138Ce/140The Ce ratio is obtained.
Further, the rhenium filament current rising speed in the step (3.2) is 200 mA/min-800 mA/min;
Further, the rhenium filament current rising speed in the step (3.3) is 30 mA/min-80 mA/min, the rhenium filament current rising speed is not too fast, otherwise, the too fast evaporation consumption of the sample is easily caused;
further, in the step (3.4), the vacuum degree of the ion source is more than or equal to 1.1 × 10 -7mbar, is obtained 138Ce/140The ratio of Ce;
Further, in the step (3.4), each measurement consists of 10-20 blocks, and each block contains 10-16 cycles;
Furthermore, the rhenium filament current increasing speed in the step (3.2) is 400 mA/min;
Further, the rhenium filament current rising speed in the step (3.3) is 60 mA/min;
Further, each measurement in step (3.4) consists of 16 blocks, each block containing 11 cycles.
The judgment method provided by the invention has the advantages that the judgment accuracy of the uranium ore source can reach 85%, the geographical tracing effect is ideal, and the method has important guiding significance for quickly realizing tracing of nuclear materials. The uranium ore geographical tracing method based on discrimination analysis can not only make full use of the existing sample information, but also can make use of the trace element content and stable isotope composition to perform accurate geographical tracing on the uranium ore.
Example 1
In this embodiment, a uranium ore geographical tracing method based on discriminant analysis specifically includes the following steps:
Collecting a uranium ore sample, and carrying out pretreatment to obtain a uranium ore sample solution;
The total number of uranium ore samples collected in the laboratory was 20, and the uranium ore samples were from three continents, north america, africa and oceania, and contained 10 samples of uranium ore produced in the united states, 7 samples of uranium ore produced in congo and 3 samples of uranium ore produced in australia.
Pre-treating a uranium ore sample:
(1.1) first grinding it into a uniform-grained powder having a grain size of about 10 μm using a ball mill to be digested into a lyso-sample;
(1.2) weighing a 500mg uranium ore powder sample by using a high-precision balance and placing the sample in a polytetrafluoroethylene beaker;
(1.3) 9m L HNO with a concentration of 15 mol/L was added in sequence 3HF at a concentration of 6m L of 20 mol/L and HClO at a concentration of 3m L of 12.4 mol/L 4Digesting, namely heating a polytetrafluoroethylene beaker on an electric hot plate to 210 ℃, and keeping the temperature for 24 hours;
(1.4) after the sample solution was completely evaporated to dryness, concentrated HNO of 10m L was used 3Dissolving;
(1.5) repeating the operation step (1.4) for a plurality of times until the dissolution is complete;
(1.6) after the sample solution was evaporated to dryness again, HCl of 2.5 mol/L concentration of 4m L and H of 0.5 mol/L concentration of 2m L were added 3BO3The mixed acid is redissolved at 90 ℃ to eliminate fluorinated complexes which may remain during digestion of the uranium ore sample.
And (4) completely digesting the uranium ore powder sample to obtain a sample solution.
And (2) measuring the content of rare earth elements in 20 uranium ore samples by adopting an inductively coupled plasma mass spectrometer:
the inductively coupled plasma mass spectrometer is used for detecting the content of L a, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb and L u of 14 rare earth elements in each uranium ore sample solution, and the inductively coupled plasma mass spectrometer is set under the conditions that the sample introduction speed is 0.3rps, the sample is stable for 35 seconds before analysis, the He gas mode is that the collection point number of the unit mass number is 3, the data collection repetition times is 3 times, the radio frequency power is 1550W, the carrier gas flow rate is 5.5m L/min, and the atomizing chamber temperature is 2 ℃ and the ICP-MS measurement result of the content of each rare earth element in 20 uranium ore samples is shown in table 1.
TABLE 1 measurement of the content (ppm) of each rare earth element and the Ce isotope ratio in uranium ore samples
TABLE 1 measurement of the ratio of the rare earth element content (unit: ppm) to the Ce isotope in uranium ore samples
And (3) measuring the Ce isotope ratio in the uranium ore sample by adopting a thermal surface ionization mass spectrometer:
in the sample coating process of the thermal ionization mass spectrum in the step (3.1), a rhenium band is generally adopted as a sample band, the length of the rhenium band is about 10mm, the width of the rhenium band is about 0.1mm, a micropipette is adopted to dropwise add the sample solution onto the rhenium band, and the sample solution is heated and evaporated to dryness under appropriate current.
Before the sample is coated, the blank rhenium belt is placed in vacuum of 10 -5Pa degassing apparatus, vacuum cleaning and degassing by heating current through rhenium strip, current for cleaning and degassing Typically 5.5A for a period of about 30 min. The program can effectively reduce the measurement background caused by hydrocarbon, reduce the interference of isobaric ions, improve the vacuum degree of the ion source during measurement and be beneficial to maintaining a good measurement environment. It is noted that to increase the ionic current strength, improve ionization efficiency and reduce isotope fractionation, micropipettes should be used to coat the sample solution in the center of the rhenium band as much as possible.
the specific sample coating conditions are that 1 mu L of sample solution is dripped on a rhenium band under the current of 0.9A, after the sample solution is dried by distillation, the current is added to 1.3A and kept for 2min, the current is added to 1.8A and kept for 30s, the current is added to 2A and kept for 2s, the current is slowly reduced to 1.3A and kept for 30s, the current is slowly reduced to 0A, and the rhenium band is taken down to be measured.
Step (3.2) the filament after sample coating is loaded into an ion source, the filament is heated to about 1000 ℃ at the speed of 200mA/min, and a central Faraday cup is adopted for receiving 140Ce+Adjusting ion flow peak shape and peak center;
Step (3.3) repeatedly adjusting the gathering parameters of the electron lens of the thermal surface ionization mass spectrometer, and optimizing the state of the electron lens of the thermal surface ionization mass spectrometer; the filament is heated at a low speed of 60mA/min to obtain the stable and strong filament 138Ce+And 140Ce+Ion current intensity;
step (3.4) when the vacuum degree of the ion source is better than 1.1 × 10 -7measurement of Ce at mbar +Ion flow, obtaining 138Ce/140The ratio of Ce. Each measurement consists of 16 blocks, each containing 11 cycles. The integration time and dead time for each cycle were set to 4.914s and 3s, respectively.
To correct for mass fractionation effects, to 136Ce/142The measured Ce isotope ratio is corrected by using an exponential law with Ce-0.01688 as a correction reference value to obtain the ratio of each uranium ore sample 138Ce/140The Ce ratio (hereinafter Ce is substituted) measurement results are shown in table 1.
(4) And (4) performing discriminant analysis on the content of the rare earth elements and the Ce isotope ratio to judge the source area of the uranium ore sample.
firstly, variable defect tolerance inspection is carried out on the rare earth element content and Ce isotope ratio of the 19 uranium ore samples, as can be seen from table 2, the tolerance of three variables of Ho, Er and L u is 0.001, the tolerance of Tm variable is 0 and is less than the minimum tolerance of 0.1, and therefore the samples need to be removed.
TABLE 2 variable Defect tolerance test
the contents of 10 rare earth elements (L a, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy and Yb) which are not removed after variable defect tolerance inspection and the Ce isotope ratio are used as independent variables of a discriminant function, and three types of U.S., Congo and Australia uranium ores are used as three groups of classification variables to establish various discriminant functions.
for example, when a No. 1 uranium ore sample from australia is taken as a sample to be discriminated from 20 uranium ore samples, the remaining 19 uranium ore samples are classified into 3 types by source, and the discrimination function coefficients are obtained by maximizing the ratio of the difference between classes to the difference within classes of the linear function of the independent variable composition by variance analysis, as shown in table 3, with the content of 10 rare earth elements (L a, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, and Yb) and the Ce isotope ratio as the discrimination function independent variable:
TABLE 3 discriminant function coefficients based on discriminant analysis method
Thus, the discriminant functions for each class are established as follows:
22.815 × L a +1.159 × Ce-21.625 × 0Pr +13.627 × Nd +18.965 × Sm-17.488 × Eu-167.850 × Gd-82.732 × Tb +92.719 × Dy-28.311 × Yb-310.948 × Ce-91682.02;
F (Congo) ═ 0.262 × L a +0.051 XCe +0.117 X0 Pr +0.007 XNd +0.157 XSm + 0.721 XEu-2.386 XGd-0.252 XTb +0.467 XDy-0.160 XYb-2.974 XCe-6.836;
F (Australia) ═ 0.023 × L a-0.022 × Ce-0.249 × 0Pr +0.106 × Nd-0.082 × Sm + 0.561 × Eu +0.545 × Gd +0.225 × Tb-0.479 × Dy +0.243 × Yb +3.460 × Ce-7.206.
Substituting the rare earth content and the Ce isotope ratio of the uranium ore sample to be distinguished into the distinguishing function for calculation, then comparing the three distinguishing function values, and if the F (American) value is the maximum, classifying the corresponding uranium ore sample into an American group; if the F (Congo) value is maximum, the corresponding uranium ore sample is classified into a Congo group; if the F (australia) value is the greatest, the corresponding uranium ore sample is classified into australia.
The judgment result obtained by adopting the judgment function is shown in a table 4, the upper half part of the table is the judgment information obtained by using a self-verification method, the judgment of all uranium ore samples of the American group, the Congo group and the Australian group can be known from the table to be correct, and the integral judgment accuracy is 100%; the lower half part of the table is discrimination information obtained by a cross validation method, and the discrimination results in the table show that except that 1 U.S. uranium ore sample is misjudged as a Congo group, 1 U.S. uranium ore sample is misjudged as a U.S. group, and 1 Australian uranium ore sample is misjudged as a U.S. group, the rest 17 uranium ore samples are all correctly discriminated, and the overall discrimination accuracy is 85%, so that the lower half part of the table is feasible and ideal in discrimination effect by adopting discrimination analysis to perform geographical tracing on the uranium ore samples based on the content of rare earth elements and the Ce isotope ratio information.
TABLE 4 discrimination results based on discriminant analysis method
Claims (9)
1. A uranium ore geographical tracing method based on discriminant analysis is characterized by comprising the following steps:
(1) Collecting m uranium ore samples from n source places, wherein the number of the uranium ore samples in each source place is not less than 3, and preprocessing the uranium ore samples to obtain uranium ore sample solutions;
(2) Respectively measuring the content of rare earth elements in m uranium ore samples by adopting an inductively coupled plasma mass spectrometer;
(3) Respectively measuring Ce isotope ratios in m uranium ore samples by adopting a thermal surface ionization mass spectrometer;
(4) Performing discriminant analysis according to the rare earth element content and the Ce isotope ratio of each uranium ore sample, and judging the source area of the uranium ore sample according to the calculation result;
Wherein, the step (4) specifically comprises the following steps:
Randomly taking 1 sample from m uranium ore samples as a sample to be distinguished, and classifying the rest (m-1) uranium ore samples into n types according to sources;
Step (4.2) determining the independent variables of the discriminant function and the discriminant function coefficients corresponding to the independent variables, and establishing corresponding discriminant functions for uranium ore samples of various source areas;
Step (4.3) substituting the rare earth element content and the Ce isotope ratio of the sample to be judged into each judgment function, and calculating to obtain the judgment function value of each judgment model;
And (4.4) comparing the discrimination function values, and determining the source corresponding to the maximum discrimination function value as the source of the sample to be discriminated.
2. The uranium ore geographical tracing method based on discriminant analysis according to claim 1, wherein: the step (4.2) specifically comprises the following steps of (4.2.1): firstly, performing variable defect tolerance inspection on the content of all rare earth elements and the Ce isotope ratio of the (m-1) uranium ore samples, eliminating the content of corresponding rare earth elements and/or the Ce isotope ratio with the tolerance less than 0.1, and taking the content of the rest rare earth elements and/or the Ce isotope ratio which are not eliminated as the independent variable of a discriminant function.
3. The uranium ore geographical tracing method based on discriminant analysis according to claim 1, wherein: step (4.2) further comprises step (4.2.2): the determination method of the discriminant function coefficient corresponding to each variable maximizes the ratio of the difference between each type and the difference in each type of the linear function composed of independent variables through variance analysis.
4. The uranium ore geographical tracing method based on discriminant analysis according to claim 1, wherein: the step (1) specifically comprises the following steps:
Step (1.1) grinding a uranium ore sample into powder with uniform particles by using a ball grinder;
Weighing a uranium ore powder sample and placing the uranium ore powder sample in a polytetrafluoroethylene beaker;
Digesting: sequentially adding HNO into a polytetrafluoroethylene beaker 3HF and HClO 4Heating the beaker on an electric heating plate for a certain time;
Step (1.4) after the sample solution is completely evaporated to dryness, utilizing concentrated HNO 3Dissolving;
Step (1.5) repeating operation step (1.4) for multiple times until the dissolution is complete;
Step (1.6) after the sample solution is evaporated to dryness again, HCl and H are used 3BO3And dissolving the mixed acid again to obtain a uranium ore sample solution.
5. The uranium ore geographical tracing method based on discriminant analysis according to claim 4, wherein: in the step (1.3), the HNO 3the concentration was 15 mol/L.
6. the uranium ore geographical tracing method based on discriminant analysis according to claim 4, wherein in step (1.3), the HF concentration is 20 mol/L.
7. The uranium ore geographical tracing method based on discriminant analysis according to claim 4, wherein: in the step (1.3), the HClO 4the concentration was 12.4 mol/L.
8. The uranium ore geographical tracing method based on discriminant analysis according to claim 4, wherein: in the step (1.3), the H 3BO3the concentration was 0.5 mol/L.
9. The uranium ore geographical tracing method based on discriminant analysis according to claim 1, wherein: the method is characterized in that the Ce isotope ratio in the uranium ore sample is measured by adopting a thermal surface ionization mass spectrometer in the step (3), and the method specifically comprises the following steps:
Step (3.1) coating a sample solution on the center of a rhenium strip filament by using a micropipette;
Step (3.2) the filament after sample coating is loaded into an ion source, the filament is heated, a central Faraday cup is adopted to receive ions to be detected, and the ion current peak shape and the peak center are adjusted;
Step (3.3) adjusting the focusing parameters of the electron lens of the thermal surface ionization mass spectrometer, and heating the filament at a low speed to obtain a stable and strong filament 138Ce+And 140Ce+Ion current intensity;
Step (3.4) when the vacuum degree of the ion source reaches a certain value, obtaining 138Ce/140The Ce ratio is obtained.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010352005.6A CN111487309B (en) | 2020-04-29 | 2020-04-29 | Uranium ore geographic tracing method based on discriminant analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010352005.6A CN111487309B (en) | 2020-04-29 | 2020-04-29 | Uranium ore geographic tracing method based on discriminant analysis |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111487309A true CN111487309A (en) | 2020-08-04 |
| CN111487309B CN111487309B (en) | 2023-05-09 |
Family
ID=71792261
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010352005.6A Active CN111487309B (en) | 2020-04-29 | 2020-04-29 | Uranium ore geographic tracing method based on discriminant analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111487309B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112180056A (en) * | 2020-10-27 | 2021-01-05 | 黄家岭 | Maotai-flavor liquor tracing method and system based on rare earth element detection |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2013061183A (en) * | 2011-09-12 | 2013-04-04 | Nisshin Flour Milling Inc | Origin discrimination method for australian wheat using trace element and isotopic ratio |
| CN105259160A (en) * | 2015-11-03 | 2016-01-20 | 中国农业科学院茶叶研究所 | West Lake Longjing tea production place identification method based on ionomics |
| CN105301089A (en) * | 2015-10-23 | 2016-02-03 | 宁波出入境检验检疫局检验检疫技术中心 | Method for measuring strontium isotope ratios of samples at different positions of strawberry tree and authenticating producing areas of waxberries |
| CN106248810A (en) * | 2016-06-08 | 2016-12-21 | 深圳出入境检验检疫局食品检验检疫技术中心 | A kind of wine place of production based on multielement and stable isotope source tracing method |
| CN107424003A (en) * | 2017-05-05 | 2017-12-01 | 浙江省农业科学院 | A kind of red bayberry place of production source tracing method based on isotope ratio rate and multielement |
| CN108982692A (en) * | 2018-07-27 | 2018-12-11 | 深圳出入境检验检疫局食品检验检疫技术中心 | The method that elemental analysis-stable isotope mass spectrum differentiates the milk powder place of production |
| CN109709200A (en) * | 2018-11-30 | 2019-05-03 | 宁波检验检疫科学技术研究院 | A method of the Rhizoma Atractylodis Macrocephalae place of production is identified based on Stable isotope ratio |
| CN109799051A (en) * | 2019-01-30 | 2019-05-24 | 航天科工防御技术研究试验中心 | A kind of processing method and system of multi coherent function spectrum |
-
2020
- 2020-04-29 CN CN202010352005.6A patent/CN111487309B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2013061183A (en) * | 2011-09-12 | 2013-04-04 | Nisshin Flour Milling Inc | Origin discrimination method for australian wheat using trace element and isotopic ratio |
| CN105301089A (en) * | 2015-10-23 | 2016-02-03 | 宁波出入境检验检疫局检验检疫技术中心 | Method for measuring strontium isotope ratios of samples at different positions of strawberry tree and authenticating producing areas of waxberries |
| CN105259160A (en) * | 2015-11-03 | 2016-01-20 | 中国农业科学院茶叶研究所 | West Lake Longjing tea production place identification method based on ionomics |
| CN106248810A (en) * | 2016-06-08 | 2016-12-21 | 深圳出入境检验检疫局食品检验检疫技术中心 | A kind of wine place of production based on multielement and stable isotope source tracing method |
| CN107424003A (en) * | 2017-05-05 | 2017-12-01 | 浙江省农业科学院 | A kind of red bayberry place of production source tracing method based on isotope ratio rate and multielement |
| CN108982692A (en) * | 2018-07-27 | 2018-12-11 | 深圳出入境检验检疫局食品检验检疫技术中心 | The method that elemental analysis-stable isotope mass spectrum differentiates the milk powder place of production |
| CN109709200A (en) * | 2018-11-30 | 2019-05-03 | 宁波检验检疫科学技术研究院 | A method of the Rhizoma Atractylodis Macrocephalae place of production is identified based on Stable isotope ratio |
| CN109799051A (en) * | 2019-01-30 | 2019-05-24 | 航天科工防御技术研究试验中心 | A kind of processing method and system of multi coherent function spectrum |
Non-Patent Citations (2)
| Title |
|---|
| ANJA MAHNE OPATI等: "Stable isotope ratio and elemental composition parameters in combination with discriminant analysis classification model to assign country of origin to commercial vegetables e A preliminary study", 《FOOD CONTROL》 * |
| XUEPENG SHAO等: "High-precision cerium isotope analysis by thermal ionization mass spectrometry using the Ce technique", 《JAAS》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112180056A (en) * | 2020-10-27 | 2021-01-05 | 黄家岭 | Maotai-flavor liquor tracing method and system based on rare earth element detection |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111487309B (en) | 2023-05-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109239176A (en) | A kind of yak meat place of production source tracing method based on multielement and stable isotope | |
| CN110687189A (en) | Method for distinguishing milk producing area tracing and identifying based on multiple elements and multiple indexes | |
| CN111272686B (en) | Hyperspectral detection method for iron grade of iron ore concentrate powder | |
| CN109030528A (en) | A kind of method that X-ray fluorescence spectra analyzes fluorine chlorine in smelting smoke dust | |
| CN111487309A (en) | Uranium ore geographical tracing method based on discriminant analysis | |
| CN103983715A (en) | TIMS measuring method for oxygen isotope ratio in uranium oxide | |
| CN111505101B (en) | Uranium ore producing area classification method based on principal component analysis | |
| CN104535598A (en) | Method for measuring isotopic abundance ratio of uranium in particles through accelerator mass spectrometry | |
| CN111505100B (en) | Uranium ore producing area classification method based on principal component-cluster analysis | |
| CN111505102B (en) | Uranium ore producing area classification method based on cluster analysis | |
| Gao et al. | An improved method for determination of isotope ratios in individual uranium particles by thermal ionization mass spectrometry | |
| CN110569581A (en) | Method for distinguishing production places of Chinese wolfberry based on multi-element combination random forest algorithm | |
| CN103983682A (en) | TIMS measuring method for oxygen isotopes | |
| CN110412116B (en) | Method for testing content of sulfur element and application thereof | |
| CN111539925B (en) | Quantitative lithium battery negative electrode lithium precipitation test method based on pixel analysis | |
| CN115015225A (en) | Method for measuring contents of phosphorus, iron and lithium in lithium iron phosphate for lithium ion battery | |
| CN110632162B (en) | Method for identifying wild astragalus and cultivated astragalus in ground | |
| CN113340975A (en) | Method for simultaneously determining 18 elements in copper ore | |
| CN113514586B (en) | Soybean origin tracing identification method based on combination of MALDI-TOF/TOF and multi-element analysis technology | |
| CN114034727B (en) | Rapid identification and quantitative detection method for niobium-rich minerals | |
| CN111220693A (en) | Hetian jade identification method for different producing areas based on multi-element content statistics | |
| CN114414653B (en) | Cherry origin tracing method based on kernel | |
| CN117907512B (en) | Sewage detection method constructed based on relation between solid phase extraction flow rate and internal standard selection | |
| CN111175368B (en) | Method for identifying addition of synthetic acetic acid in brewed vinegar | |
| CN117783261A (en) | Method for measuring element content in single-particle uranium mineral |
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 |