Preparation of Ion-adsorption Type REE Monitoring Samples
-
摘要: 离子吸附型稀土矿是中国极其重要、世界罕见的矿床类型，是中国的优势矿产资源。目前现有的稀土标样全部为稀土矿石成分分析标样，稀土元素以氧化物形式稳定存在，无法淋滤浸取，不能对离子吸附型稀土淋滤过程进行监控。为进一步满足离子吸附型稀土资源勘查和评价需要，本文按照导则《标准物质定值的通用原则及统计学原理》（JJF 1343—2012），制备了三种岩性共9个离子吸附型稀土监控样。样品采自南岭地区典型富轻稀土（HREE）离子吸附型稀土风化壳，经干燥、球磨至200目后混合机混匀，以硫酸铵淋滤、电感耦合等离子体质谱法（ICP-MS）测试离子相稀土元素含量，检验样品均匀性，结果表明F值小于临界值F0.05（19，20），样品均匀性良好。在两年内对样品进行4次稳定性检验，在95%置信度时│β1│ <t0.05×s（β1），未发现明显不稳定变化，表明稳定性良好。由8家技术权威的实验室进行协作定值，经过统计计算给出各离子相稀土元素含量的加权平均值和扩展不确定度，定值结果涵盖除Sc以外的15种稀土元素。该系列监控样的研制能够为离子吸附型稀土矿产资源评价和有效利用等工作提供计量支撑。要点
(1) 制备了9个离子吸附型稀土监控样，定值成分达15种。
(2) 监控样的均匀性、稳定性及定值均满足JJF 1343—2012的要求。
(3) 监控样的制备为完善稀土标准物质体系奠定了基础。
HIGHLIGHTS(1) Preparation of 9 ion-adsorption type REE monitoring samples and 15 elements were characterized as certified values.
(2) The homogeneity, stability and certified values met the requirements of JJF 1343—2012.
(3) The preparation of monitoring samples layed a foundation for improving the rare earth reference materials system.
Abstract:OBJECTIVESIon-adsorption type rare earth ore is extremely important in China and rare in the world. It is China's dominant mineral resource. At present, all the rare earth elements in the existing REE reference materials are oxide minerals, not ion-adsorption type, and cannot be exchanged with strong electrolytes. These reference materials cannot be used to monitor the leaching process of ion-adsorption type REE.METHODSThe samples were collected from the typical ion-adsorption type REE weathering crust in Jiangxi Province. After being dried and ball-milled to 200 mesh, they were mixed for 1.5h. The contents of rare earth elements in the ionic phase were analyzed to test the homogeneity by ammonium sulfate leaching and inductively coupled plasma-mass spectrometry (ICP-MS).RESULTSResults showed that the F values for the variance test were less than the threshold, which indicated a good homogeneity. The stability of monitoring samples was tested 4 times in two years, and no statistically significant changes were observed, indicating good stability of the samples. An accurate and reliable leaching analytical method was used in 8 technologically significant laboratories to finalize the contents of ionic phase rare earth elements. The determined values of nine monitoring samples were given by statistical calculation, including weighted average value and total uncertainty results which were from uncertainties in certified values, between-bottle homogeneity and long-term stability. These nine monitoring samples have certified values of 15 ionic phase rare earth elements except Sc.CONCLUSIONSThe reference materials can be used to monitor the leaching process and provide metrological support for the analysis of rare earth elements in ionic phase during the evaluation and effective utilization of REE mineral resources.
-
Keywords:
- ion-adsorption type rare earth elements /
- leaching /
- monitoring samples /
- certified value /
-
-
释光测年技术最早可追溯至20世纪60年代，Daniels等首次提出天然矿物对电离辐射的发光响应可应用在加热(温度高于500℃)的古代陶器测年中，称为热释光测年(thermoluminescence dating)[1-2]；70年代，Aitken等发现热释光的电子信号还可以被光激发[3]，由此将热释光测年应用在非受热的沉积物中[4]。1985年Huntley等[5]提出释光技术应用在沉积物测年领域，使用绿光照射矿物颗粒激发冷光，并将此命名为“optical dating”。随后近十年发展为光释光测年(optically stimulated luminescence dating, OSL)，对沉积物最后一次曝光(或受热)事件的年代测定，测年矿物主要为石英和钾长石，而锆石、方解石、燧石虽有释光特征却并不常用。石英释光信号晒退快但饱和剂量低，目前比较成熟的石英释光测年可用于100~200ka以来第四纪沉积物的年龄测定；钾长石晒退慢，但红外释光(IRSL)信号饱和剂量更高，目前比较成熟的钾长石释光测年可用于250~300ka以来全球不同区域沉积物的年龄测定。这些都远远大于常用的放射性碳-14测年40~50ka的测年上限，因此被广泛应用在考古学、人类学、地质考古学、地貌学、第四纪地质环境等方面。
传统石英及钾长石测年主要使用每个测片上均为几百至上千颗粒的单片测年技术[6-7]，后来随着释光测年技术进步，发展了可以开展测试每个矿物颗粒释光信号的单颗粒释光测年技术[8]。早期，Lamothe等[8]对14C年龄为9.5~10ka的海沙进行长石单颗粒释光测试，认为样品不完全晒退和长石的异常衰退很大程度上影响了等效剂量(equivalent dose, De)的可靠测量，首次证明了单颗粒释光测年的可行性。Murray等[9]对澳大利亚六条河流的漫滩相沉积物进行释光测年时发现一小部分晒退不好的石英颗粒最终污染了整个样品，残余剂量约0.02Gy导致年龄被高估几百年，认为应该对样品的每个颗粒单独测试来筛除被污染样品。Murray和Roberts[10]将澳大利亚南部风成沉积物石英颗粒每一粒单独制成单片进行单颗粒测年，所得结果与多颗粒单片再生剂量法、多颗粒单片附加剂量法、14C和热释光年龄吻合，并系统分析了单颗粒结果离散的原因。Roberts等[11]对澳大利亚Jinmium岩画遗址的样品进行单颗粒研究，区分出未晒退或晒退不完全的样品，认为是由于岩石风化后颗粒进入新的地层中导致。Olley等[12]对深海沉积物样品进行单颗粒测年，与14C进行比较，结果显示光释光测年在1.5~50ka范围内与14C结果符合，提出单颗粒方法可以有效地应用于晒退不均匀的样品。随后二十年，单颗粒释光测年技术不断发展[9-10, 13]，极大程度上提高了光释光测年的准确度。释光测年最后结果取自95%置信区间即中间值的2倍标准差(2σ)之内，对于地质考古样品单颗粒De测量结果通常较为分散，主要可能是由于以下原因导致：①样品本身的不均一性，包括样品热稳定性不同[14-15]，颗粒之间被晒退能力不同[16-17]，导致最终De值分散。②在样品埋藏时期受到不同程度的周边环境辐射，尤其是β辐射变化较大，如硅质骨料和周边长石颗粒中含有的40K可能会导致石英颗粒成为高辐射的热点，而土壤中的碳酸盐和贝壳则可能会导致石英低辐射区，辐照的不同导致矿物颗粒最后储存的能量不同，最后De值表现分散。③沉积后的干扰，如沉积地层中埋藏的碎石、风化的基岩崩塌掉落的颗粒、生物活动、人类活动导致地层扰动，最终导致年轻颗粒进入较老的沉积物中，或较老的颗粒进入年轻沉积物中。④沉积物最后一次埋藏前为不完全晒退或者晒退不均一，不同沉积类型的沉积物(风成沉积物、河流相与湖相沉积物、重力堆积物与洪积物、洞穴堆积物)晒退程度相差较大[18-19]。这些扰动、侵入，不同材料的辐射等原因，最终都会传递到De上，导致在地质考古点中光释光样品的De值更为分散且分散特征多样化[20]。对于这些De值比较分散的样品进行光释光年代测定时，单颗粒技术尤为必要。目前单颗粒释光测年方法被国际同行应用在地质考古领域取得大量高水平成果[11, 13, 21-26]。
在我国，石英测年及钾长石pIRIR测年已广泛应用于西北干旱区及青藏高原地区沉积物的可靠年龄测量，石英释光测年可用于50~60ka以来沉积物样品年龄测定，钾长石pIRIR测年已成功用于我国西北干旱区沙漠、黄土、湖泊沉积250ka以来样品可靠年龄的测定[27-29]。但是在我国多为使用多颗粒的释光测年，目前单颗粒释光测年应用较少，未来应用前景广泛。
1. 光释光测年主要原理
光释光测年理论基础是能级理论(energy-level)，自然界中一些矿物(石英、长石等)普遍存在一种现象，当它们暴露在α、β、γ辐射下时，可以在晶体结构内储存一小部分辐射传递过来的能量，这个能量在见光或者受热后归零并释放出冷光，当再次被埋藏时又可以继续储存能量，矿物的这种储存能力像充电电池一样可以反复储存—释放—储存。被储存的能量表示为公式(1)，释光测年应用中通常默认
$ {\dot D(t)}$ 为常量，则年龄A函数可简化为公式(2)，通过测定矿物储存的能量和能量吸收速度来确定沉积物最后一次曝光(或受热)时的年代(图 1)[30]。![]()
$ {D_{\rm{e}}} = \int\limits_0^A {\dot D(t){\rm{d}}t} $
(1) 式中：
$ {\dot D(t)}$ 指环境剂量率，是时间t的函数。如果$ {\dot D(t)}$ 是常量，则年龄A函数可简化为De/$ {\dot D(t)}$ 。$ A = {D_{\rm{e}}}/\dot D(t) $
(2) 式中：De指等效剂量，是电子陷阱自上一次见光归零之后埋藏至今储存的能量总和，单位是Gy(1Gy=1J/kg)；
$ {\dot D(t)}$ 指环境剂量率，又称年剂量率，单位是Gy/a。如图 1所示，光释光测年原理可分为以下四个步骤理解。
(a) 电离：自然界中矿物被埋藏时，受到周围环境中的放射性物质(主要是U、Th、K)衰变释放产生的能量和宇宙射线作用，价带中的晶体电子发生电离，产生自由电子-空穴对，自由电子被晶格缺陷产生的电子陷阱捕获，空穴则位于能量较低的释光中心。
(b) 储存：随着埋藏时间增加，被捕获的电子越来越多，但是导带下被捕获电子的寿命取决于捕获位置深度，对于测年而言，要求捕获电子足够深(寿命足够长)，以至于在测年范围内的泄露数可忽略。
(c) 释放：在接受阳光照射或受热后，在光敏陷阱中的电子进入导带，随后被其他陷阱捕获或者与释光中心结合，并释放出冷光，照射一段时间后，光敏陷阱清空。
(d) 实验：样品采集到实验室时，使用激光照射样品，将地质历史埋藏时期被捕获的电子释放，再对其加载不同剂量的β剂量，分别测试每种剂量下的电子释放数，建立样品生长曲线，再将自然信号对应在生长曲线上得到相应的β剂量数，最后除以年剂量率得到年代。
前人对如何计算年剂量率已有详细研究[30-35]，单颗粒释光测年主要是对等效剂量的一种测量方式，因此本文对年剂量率仅作简要阐述。
$ {\dot D(t)}$ 是指每个颗粒每年吸收的α、β、γ剂量和宇宙射线剂量的总和。年剂量率主要来自于样品本身和周围的放射性元素(主要是238U、235U、232Th以及衰变链中子体，40K)衰变产生；年剂量率少量由宇宙射线贡献，由样品的海拔、纬度、经度和埋藏深度决定[31]。α射线射程短，仅百分之几毫米(~25μm)，仅作用在矿物颗粒外表面，对样品进行前处理时使用氢氟酸将受到α射线影响的外表面去除，因此计算年剂量率时主要计算β、γ剂量率和宇宙射线的微量影响，但是对于长石颗粒还需要考虑颗粒本身内部40K和87Rb的辐射影响。De的测量目前主要使用单片再生剂量法(single aliquot regenerative-dose method, SAR)[36-37]和简单多片再生法(simplified multiple aliquot regenerative-dose protocol, MAR)[7]。但是，这些方法的De测试结果是几百至上千颗粒的平均状态，无法对样品单个颗粒[8-9, 38-39]进行分析。单颗粒释光测年技术通过测量单个石英或长石颗粒，分析不同颗粒之间的差异性，基于数理统计分析获得样品可靠的De结果[10-11, 40-43]。
2. 单颗粒释光测年技术实验流程
2.1 单颗粒释光测年样品前处理
一般地质考古的光释光样品在采集时选择不透明钢管，或者在弱红光灯(波长约630nm)[44]的条件下采集黑色塑料袋袋装样品。样品运至实验室之后，整个实验在弱红光灯条件下进行操作，将样品进行湿筛或者干筛得到180~250μm样品，加入10%盐酸去除碳酸盐部分、30%双氧水去除有机质[45]，清洗样品以后进行二筛得到180~212μm或180~250μm样品，烘干之后用重液(多钨酸钠或多钨酸锂)进行矿物分离。分离石英的重液组合为2.75、2.62g/cm3，分离钾长石的重液组合为2.75、2.62、2.58g/cm3[32]。对分离出来的石英、钾长石进行氢氟酸溶蚀，溶蚀的目的主要是为了溶蚀样品外层，去除外部α射线对样品年剂量率的影响，同时由于长石更容易被溶蚀，对于石英样品用氢氟酸溶蚀可以除去掺杂的长石。溶蚀石英选择40%氢氟酸反应45min或45%氢氟酸反应40min；溶蚀长石选择10%氢氟酸反应40min(少量长石则为10min)。之后加入足量10%盐酸，去除反应产生的氟化物。最后用蒸馏水洗净之后烘干，再次干筛出180~212μm或180~250μm石英或钾长石矿物，完成样品前处理。
2.2 De测量程序
Murray等[10]开始进行单颗粒尝试时是将每个粗颗粒单独粘在一个不锈钢测片上，每个测片独立测试。这种方法非常耗时，并且之前的多颗粒仪器使用的是卤钨灯或发光二极管，发射40mW/cm2的绿光或蓝光照射石英颗粒，这种强度下只有非常亮的颗粒才能发出可探测到的冷光，因此极大地限制了单颗粒测试。随着技术的进步，目前Automated Riso luminescence reader这款设备可以开展单颗粒样品的全自动批量测量，使用10mW绿光激光(波长532nm)照射石英，使用150mW红外激光(波长830nm)照射长石，可以准确地对每一个颗粒进行激发。绿光激光照射在样品上能量可达40mW/cm2，聚焦20μm的直径范围内，红外激光照射到样品上最强可达到400mW/cm2，考虑到稳定性使用时一般选择最高不超过200mW/cm2[46]。单颗粒测片为镀金铝片，每个片子有10×10个空洞，每个空洞直径为300μm，深度300μm，测试的单颗粒最佳粒径为180~250μm，每个样片可装载100颗样品。另外，在一个测片上的所有颗粒可以同时进行90Sr β源辐照或X射线辐照，极大地提高了测量效率。
石英单颗粒释光测试流程主要由SAR发展而来，利用IRSL信号检测长石污染，石英110℃TL峰校正信号，小剂量的实验剂量进行灵敏度校正，采用绿光刺激下125℃时释放出来的信号，测试最接近石英OSL信号的325℃ TL峰值信号[6, 47-50]。钾长石单颗粒释光测试考虑到长石存在由电子-空穴对之间的隧道复合而引起的异常衰退[30, 51-53]，目前常用两步加热法[54-56]和多步加热法[57]。
2.3 样品纯度筛选
对于石英单颗粒测年，为了检验提纯获得的石英颗粒是否受到发光信号更强的长石信号的影响，最常用的方法是开展OSL IR depletion ratio检测，在De测试程序最后使用红外光刺激，计算IR刺激之后的OSL信号/IR刺激之前的OSL信号比值，如果这个比值近似1则表示为纯净石英[58]。此外，还可以用不同温度下的OSL信号衰减曲线(decay curve)。对于石英，随着温度的增加衰减更快，但是要考虑高温下的热淬火(thermal quenching)；而长石的衰减曲线则没有明显的区别[58]。后一种方法比较耗时耗样品，前一种方法可以直接加在De测量的最后一步，并不需要准备额外的样品，测试时间更短，因此被广泛应用。
2.4 De测量可靠性分析主要标准
单颗粒释光测年中，为了获得准确的De，对每个颗粒进行测试后需要通过严格的条件筛选出合适的样品颗粒来计算De[6, 59-60]。目前De测量可靠性分析主要有以下5条标准。
(1) 单颗粒测年需要选择足够亮的颗粒，筛去不发光的杂质颗粒、亮度不够或不发光的石英或钾长石，利用自然剂量下加实验剂量的实验信号(Tn)将不亮的颗粒筛去，条件如下：①Tn是否大于背景值的3倍标准偏差(3σ)，舍去Tn小于3σ的颗粒[59]；②Tn的相对标准差，一般筛去Tn大于20%~30%的颗粒。
(2) 利用Tn选择出快组分多的颗粒。样品受到光刺激之后释放出来的信号可以分为快组分、中组分、慢组分[30]，释光测年主要采用快组分[61]，因为快组分信号在百万年尺度范围内更稳定[14]，同时在自然界中也最容易晒退[62]，这里主要是指石英颗粒，长石颗粒多考虑异常衰退问题。快组分的贡献程度将很大程度上决定等效剂量计算的精度[50]。Ballarini等[63]在对小于300年的年轻沉积物进行单颗粒测试时提出使用前背景法(EBG)计算快组分，即背景信号采用紧邻初始信号积分区间之后的那部分，这一方法相比传统的背景信号选择方法(选择信号积分区间的最后一段的信号平均值)，最终将中组分的贡献减少83%，慢组分的贡献减少99%。这一方法虽然更大程度减弱了中、慢组分的影响，但是同时也减少了用于计算等效剂量的净信号，降低了De分析的精度。Cunningham等[64]对4个年轻样品进行了单片法测试，进一步发展了EBG方法，得到结论是对于假定由快速分量控制的信号，采用大约是初始信号长度2.5倍的前背景值计算，不会影响等效剂量的精度。另外，也可以使用快组分比例(fast ratio)，即快组分减去背景值除以中组分减去背景值[61, 65]。在等效剂量计算的筛选中可以用fast ratio值来区分快速和非快速主导的CW-OSL信号[65]。
(3) Recycling Ratio。测试程序最后一步为重复测量一次相同的实验剂量，前后两次结果的比值为recycling ratio。样品受到反复加热、反复给实验剂量后灵敏度发生变化，通过recycling ratio来选择灵敏度变化不大的样品。Murray等[6]认为recycling ratio在1.0±0.1范围内可接受。也有学者如Jacobs等[60]认为recycling ratio必须要在1.0±0.1范围内且大于1.0才可接受。
(4) Recuperation。测试程序中多次给实验剂量之后还需要测一个实验剂量为0的信号，将这一步的校正后结果与自然信号或者给定的最大的实验剂量的灵敏度校正后结果的比值为recuperation。Murray等[6]建议这个比值选择不应超过5%。
(5) 曲线拟合程度。单颗粒测试之后每个颗粒可以采用不同的曲线拟合，通常在单指数饱和、指数加线性、双指数饱和曲线中选择最佳的曲线，剔除曲线拟合度较低或者通过外延得到自然信号值的颗粒。
3. 单颗粒De测量结果的主要分析模型
技术和方法的发展极大地提高了大量单颗粒数据的计算效率。对于一个理想样品，单颗粒之间的De值所观察到的离散，应该与仪器测量误差一致。如本文引言所述，地质考古样品De值分布的离散程度通常大于仅从统计学角度考虑的结果，这些过度分散(over dispersion, OD)是指不能仅用测量不确定性来解释，而是由于实验的其他原因或者颗粒之间自然信号的变化引起[60]，如仪器测量的重现性[66-71]。已有研究表明，信号较好的亮颗粒OD值可能低于1%，而地质考古样品由于其特殊性，OD值经常高达20%，甚至有些暗颗粒OD值可能超过50%[60, 72-73]。面对大量分散的De数据，目前最常用的方法是使用Radio plot图解来提取有效信息，它可以显示每个颗粒的De值及其测量不确定性[74-76]。Radio plot图(图 2)从左侧(标准化估计值)轴的原点经过数据点直线延伸到右侧的Y轴，交点为De值，数据点垂直延伸到底部X轴为该数值的测量不确定性，通常上列用百分数表示相对误差，下列用其倒数表示精度。Radio plot图解借鉴裂变径迹的统计方法，可以统计传统的加权直方图(概率密度函数)无法统计的一些精度较低、跨越多个分区的De值[79]。
![]()
图 2 从考古遗址提取的单一石英颗粒De值Radio plot图解(a)南非Blombos洞穴考古遗址中晒退较好的样品，使用中心年龄模式得到等效剂量值[77]。
(b)南非Blombos洞穴考古遗址中晒退不均匀的样品，使用有限混合模式得到等效剂量值[77]。
(c)澳大利亚最早的人类遗址点Lake Mungo中二次晒退样品，使用最大年龄模式得到等效剂量值[78]。
(d)澳大利亚Jinmium岩画遗址点不完全晒退样品，使用最小年龄模式得到等效剂量值[11]。Figure 2. Radio plots of the De values obtained from single grains of quartz extracted from archeological sites单颗粒实验获得的所有颗粒De信息，通常使用以下5种模式开展De的分析及可靠年龄的计算。
(1) 普通年龄模式(common age model)。这一模式主要用于在沉积前样品晒退好，并且在埋藏过程中受到均一的外界(和内部)辐照，后期也没有受到任何干扰的这一类样品。该模式应用于非对数化的De值估计，因此使用这一模式计算时，随着单个颗粒的样本数量增加，De的标准误差不会发生明显变化，这也表明普通年龄模式适用于样品离散度非常低的情况，然而自然界中这类样品极少，因此普通年龄模式适用范围很窄[76, 80]。
(2) 中心年龄模式(central age model)。如普通年龄模式中所述，自然界中几乎所有样品的单颗粒数据之间都会有离散度(OD)，不适合普通年龄模式，这时就需要选择中心年龄模式。这一模式会将每一个颗粒的数据都放入最终De值的计算中，当中心年龄模式的标准偏差为0时，则和普通年龄模式得到的De值相同[76, 78]。对于晒退较好的地质考古样品，OD值一般小于20%[12, 17, 72]，De值相对集中分布在一个组分中，通常使用中心年龄模式计算结果。
Jacobs等[59, 77]对南非Blombos洞穴进行了详细的释光测年，中石器M1期地层ZB4样品大部分De值落在2σ之内，单颗粒中心年龄模式得到De值为79.5±1.5Gy(图 2a)，离散值为13.5%，与5mm大测片SAR和2mm小测片SAR得到的De值83.0±1.8Gy、79.3±1.0Gy相吻合。Olley等[78]对澳大利亚最早的人类遗址点Lake Mungo进行单颗粒测年，中心年龄模式得到De离散值为19%，最终年龄为37.0±1.6ka，与之前使用澳大利亚滑移法得到的38.2±1.3ka、37.8±1.9ka[21]匹配。
(3) 最大年龄模式(maximum age model)。这一模式主要适用于最后一次沉积前已经彻底晒退并且在随后的沉积过程中有年轻颗粒混入，或者沉积前晒退很好而沉积后有部分样品又再次被晒退的样品。在这些情况下，在最后一次沉积前晒退很好并且之后一直被埋藏的颗粒是可得到正确De值的颗粒。最大年龄模式应用较窄，如应用在澳大利亚裸露岩石上的岩画测年，现代阳光的照射能将岩石表面几毫米处的颗粒完全或部分晒退，因此我们计算岩画年龄使用最大年龄模式[81-82]；再如Olley等[78]对澳大利亚最早的人类遗址点Lake Mungo样品进行了单颗粒释光测年，由于该样品是用半透明树脂保存30年的沙丘砂，利用LM-OSL信号排除了保存过程中再次晒退的颗粒，使用最大年龄模式计算得到年龄值为41± 4ka(图 2c)。与Bowler等[21]对同一沉积单元获得的OSL年龄结果一致，证明了最大年龄模式的可靠性。
(4) 最小年龄模式(minimum age model)。这一模式与最大年龄模式相反，主要适用于样品在最后一次埋藏前有不完全晒退[76]，或发生了沉积后的混合，尤其是含水沉积物、土壤和考古沉积物等可能是由混合年代沉积物组成[11, 13, 17-18, 40, 77, 83-88]。不完全晒退沉积物的De值因为残余信号，因此系数γ比晒退完全的沉积物的De值大，因此该年龄模式通过γ来确定最后计算的De值被截取的那一段正态分布。当然，除了测量误差和不均匀晒退、沉积混合等造成De值的扩散，对于完全晒退的样品颗粒同样可能存在过度分散，因此该模式还需要每一个De值的相对标准误差上增加一个σb值(OD)。通常可以从同一剖面的相同矿物中晒退较好的样品中获得σb值，理想情况下是从同等年龄样品中获得，这样能更好地考虑时间等相关因素[89]。但是，如果缺乏此样品，单颗粒石英样品σb值可以使用20%、30%来进行敏感性实验[90]。
Roberts等[11, 13, 76]对澳大利亚Jinmium岩画遗址点进行了详细的单颗粒释光测年，认为成岩地层中存在大量不完全晒退颗粒，而之前的TL年龄未排除这部分颗粒，得到的116~176ka人类最早的岩石艺术时代和50~75ka人类最早走出非洲的时代不准确，光释光单颗粒测年使用最小年龄模式，排除掉不完全晒退的颗粒，得到这些岩石雕刻的年龄不早于10ka(图 2d)[11]，与14C年龄吻合，而使用光释光多颗粒方法以及单颗粒的中心年龄模式得到的结果都偏老。
(5) 有限混合模式(finite mixture model)。这一模式主要适用于不适合前几种模式的其他样品，比如与最小年龄模式相关的样品中混入了更年轻的样品或有些颗粒受到低于平均水平的辐射(图 2b)[77]，再如与最大年龄模式相关的样品中有些颗粒并未完全晒退，尤其是地质考古样品，人类活动、动植物的干扰可能导致不同年龄单元的颗粒混合，采样时样品的污染，地层倾斜导致采样时样品跨越多个沉积单元，样品本身的不均匀晒退等情况[77, 86, 88, 91-94]，这些情况可以考虑使用有限混合模式[77, 95-96]。
Roberts等对澳大利亚Jinmium岩画考古点[11, 13]和Malakunanja考古点[97]样品进行了详细的单颗粒测年，随后Roberts等对澳大利亚Cooloola地区风成砂进行完全晒退后分别辐照5、10、20Gy，测试之后利用有限混合模式能很好地将这三类颗粒分开，验证了有限混合模式的合理性[95]。Roberts等[81]和Galbraith[96]进一步发展优化了有限混合模式。Jacobs等[77, 98]将有限混合模式运用于南非Blomobos洞穴考古遗址中，将之前得到的中石器时代最老沙丘年龄70ka提早到143.2±5.5ka。
该模式可以将各个独立单元区分开来，应用有限混合模式时需要确定好分几个组分，然后根据拟合得到的最大似然估计来确定结果是否合适，OD值确定方法同最小年龄模式。
4. 单颗粒测试所需样品数量
到底单颗粒测年需要测量多少样品才能得到满意结果？这主要取决于样品充分晒退的数量、真实De值的分散情况，以及年代要求的精确度等等。无论使用哪种模型，样本量越大越好，这样可以包含来自各个混合成分足够多的颗粒，尤其是当完全晒退比例较少时，或者OD值过大时，则可能需要相当大的样本量。Rodnight[99]提出对于部分晒退完全的样品，至少需要50个独立的De值。而对于只有少量完全晒退、高度分散的样品而言，则需要更多的颗粒；对于晒退完全且后期没有被污染的样品，则只需要少量De值。
5. 问题与展望
光释光测年是一种重要的第四纪测年手段，而单颗粒技术的发展进一步提高了测年精度，尤其是对于晒退不充分等原因造成De分散的样品，提供了新的测年方法。然而单颗粒技术仍然存在很多有待进一步开展深入研究的方面，比如等效剂量值超分散的原因，实验表明一些完全晒退后辐照相同实验剂量的样品，也存在等效剂量值分散的现象[13, 72, 77, 95]，过度分散将会导致最终单颗粒释光测年结果精度下降。可针对同一样品单个颗粒释光行为的根本区别开展进一步研究，改良测量和分析程序，更好地理解等效剂量值的含义、随机不确定性的来源。针对具有不同物源和沉积过程的地质考古样品，改进测量程序，发展新的数据分析方法，更好地理解测量获得的De值的含义，对于单颗粒释光测年的广泛应用具有重要意义。
此外，单颗粒释光测年结果尽可能与石英及钾长石释光单片法、多片法，以及其他可靠的独立年龄进行对比分析，以及多种测年方法交叉验证依然是建立第四纪地质考古样品可靠年龄标尺的主要途径。
-
表 1 监控样均匀性检验结果
Table 1 Analytical results of the homogeneity test for monitoring samples
样品编号 参数 La Ce Pr Nd Sm Eu Gd Tb Dy Ho Er Tm Yb Lu Y HS-1 X(×10-6) 462.2 4.8 97.5 383.1 69.3 6.2 59.5 9.7 45.2 7.5 17.4 2.1 13.2 1.6 275.4 Q1 150.1 0.007 111.5 368 22.8 0.1 17.2 0.3 13.6 0.5 0.6 0.006 0.4 0.003 20.1 Q2 232.5 0.07 46.1 186.5 34.2 0.09 10.1 0.1 5.9 0.4 2 0.03 0.8 0.004 109.6 F 0.51 0.08 1.91 1.56 0.53 0.87 1.35 2.01 1.82 0.93 0.22 0.17 0.37 0.52 0.14 RSD(%) 1.08 0.73 4.42 2.04 2.81 2.11 2.84 2.13 3.33 3.92 1.76 1.44 1.90 1.32 0.66 HS-2 X(×10-6) 851.6 46.6 226.7 849.4 108.4 6.7 54 6.4 24.6 4.2 12.5 1.5 7.5 0.96 94.4 Q1 329.9 26.4 266.7 1697.9 12.7 0.05 31.4 0.03 0.2 0.05 0.3 0.002 0.2 0.0009 38.4 Q2 1523.1 15.9 105.3 1425.2 19.3 0.03 12.9 0.02 0.4 0.02 0.1 0.005 0.3 0.004 19.4 F 0.17 1.31 2.00 0.94 0.52 1.80 1.93 1.06 0.42 1.92 1.70 0.26 0.55 0.16 1.56 RSD(%) 0.87 4.50 2.94 1.98 1.34 2.52 4.24 1.85 0.77 4.43 4.44 4.84 2.44 1.28 3.50 HS-3 X(×10-6) 920.9 28.9 262.2 804.8 112.4 7 41.6 4.6 17.7 2.8 8.8 0.8 5.6 0.7 59.3 Q1 2672.5 3.2 544.5 668.1 4.5 0.09 37.2 0.6 1.6 0.05 0.001 0.001 0.3 0.002 28.5 Q2 1225.1 2.1 485.8 456.4 29.3 0.08 33.5 0.2 4.8 0.02 0.5 0.003 0.2 0.0007 47.8 F 1.72 1.20 0.88 1.16 0.12 0.90 0.88 2.06 0.26 1.96 0.00 0.36 1.33 2.03 0.47 RSD(%) 2.29 2.54 3.63 3.80 0.77 1.79 4.98 4.81 2.90 4.22 0.18 1.69 4.21 2.57 3.68 BZ-1 X(×10-6) 101.9 9.4 16.8 56.6 10.2 2.1 13.3 2.1 11.4 2.1 5.5 0.5 3.3 0.4 64.6 Q1 56.7 0.3 2.76 17.8 0.3 0.02 0.8 0.03 0.09 0.0001 0.07 0.0002 0.1 0.002 15.5 Q2 96.4 0.004 2.16 31.4 2.2 0.03 0.8 0.04 0.06 0.0005 0.1 0.0009 0.1 0.0008 8.5 F 0.46 1.68 1.01 0.45 0.10 0.39 0.80 0.60 1.28 0.21 0.56 0.20 0.75 1.57 1.44 RSD(%) 3.01 4.40 4.04 3.05 2.18 2.57 2.77 3.54 2.52 2.53 1.95 4.40 4.14 4.33 4.19 BZ-2 X(×10-6) 158.4 5.5 29.4 104.2 17, 1 3.2 19.7 2.5 14.4 2.8 7.4 0.9 4.8 0.6 71.8 Q1 81.9 0.4 3.71 24.5 0.7 0.003 3.4 0.01 0.02 0.04 0.04 0.0006 0.01 0.001 1 Q2 33.6 0.2 1.52 9.1 1.92 0.02 1.5 0.005 0.4 0.02 0.06 0.002 0.02 0.0008 5.7 F 1.93 1.57 1.93 2.13 0.29 0.18 1.74 1.80 0.06 1.64 0.57 0.31 0.58 1.32 0.14 RSD(%) 2.99 3.98 2.68 1.94 2.00 0.77 4.97 2.33 0.45 4.87 1.13 4.63 0.90 2.68 0.57 BZ-3 X(×10-6) 270.8 13.3 57.1 206 29.6 4.3 23.3 2.3 11.1 2.1 5.6 0.6 3.6 0.4 53 Q1 253.5 0.01 91.6 16.7 6.4 0.02 1 0.02 0.2 0.1 0.7 0.003 0.1 0.0006 10.9 Q2 203.2 0.3 54.7 46.7 4.1 0.01 0.4 0.008 0.3 0.007 0.3 0.003 0.1 0.001 5.6 F 0.98 0.03 1.32 0.28 1.25 1.58 2.03 2.15 0.47 1.37 2.04 0.97 0.65 0.44 1.54 RSD(%) 2.40 0.33 4.00 0.81 3.50 1.68 3.35 4.69 1.42 4.12 4.49 4.27 3.89 2.93 3.35 CJ-1 X(×10-6) 44.5 1.4 16.7 90.7 57.2 0.4 60.2 9.9 52.8 9.6 30.1 4.4 27.3 3.9 237.7 Q1 1.8 0.002 0.33 13.5 114.2 0.001 9.6 0.5 27.5 0.1 8.6 0.09 2.3 0.01 230.3 Q2 3.7 0.002 0.54 22.7 78.1 0.0003 26.3 0.2 10.9 0.1 7.3 0.03 1.3 0.04 924.2 F 0.40 0.69 0.49 0.47 1.15 2.12 0.29 1.98 1.99 0.97 0.93 1.98 1.46 0.27 0.20 RSD(%) 1.26 1.28 1.41 1.65 7.63 4.29 2.11 4.58 4.05 2.43 3.99 4.44 2.28 1.22 2.61 CJ-2 X(×10-6) 425.7 15.2 101.8 336.8 71.2 5.6 46.3 6.5 27.3 4.5 9.7 1.1 5.4 0.7 102.3 Q1 400.4 0.2 134.5 1309.4 18.7 0.06 28 0.08 0.2 0.04 0.4 0.002 0.05 0.002 213.6 Q2 1055.3 0.1 122.5 753.8 24.2 0.03 11.4 0.03 3.6 0.02 0.3 0.001 0.02 0.001 170.3 F 0.30 1.88 0.87 1.37 0.61 1.36 1.94 2.12 0.05 1.88 1.13 1.67 1.72 1.19 0.99 RSD(%) 1.92 2.99 4.65 4.38 2.48 3.72 4.86 3.79 0.74 4.46 2.77 4.83 4.04 2.82 4.83 CJ-3 X(×10-6) 342.8 12.4 78.8 249.7 55.8 6.9 56.7 10 64.4 11.7 33.4 3.8 24.4 3.4 315.3 Q1 221.9 0.5 68.9 48.9 14.2 0.07 17.5 0.2 6.8 0.7 1.1 0.1 1.4 0.01 529.5 Q2 118.8 0.2 93.2 87.3 8.7 0.12 7.8 0.1 5.7 0.3 3.2 0.08 1.6 0.04 338.5 F 1.47 1.85 0.58 0.44 1.29 0.48 1.75 1.87 0.95 1.70 0.28 1.07 0.67 0.30 1.23 RSD(%) 1.77 2.45 4.30 1.14 2.76 1.60 4.12 4.12 4.11 2.88 1.29 4.23 1.97 1.44 2.98 注：X为浓度平均值；Q1为组间差方和；Q2为组内差方和；F为统计量；RSD为相对标准偏差。 
下载: 导出CSV
表 2 监控样稳定性检验结果
Table 2 Analytical results of the long-term stability test for monitoring samples
样品编号 参数 La Ce Pr Nd Sm Eu Gd Tb Dy Ho Er Tm Yb Lu Y HS-1 X(×10-6) 516.5 5.4 99.9 354.8 74.4 6.8 66.6 9.0 49.3 7.9 18.8 2.0 11.2 1.5 222.5 β1 -7.902 -0.063 -1.269 -4.616 -0.738 -0.051 -1.079 -0.074 -0.256 -0.099 -0.197 -0.016 -0.148 -0.006 -3.143 β0 611.273 6.131 115.115 410.186 83.214 7.422 79.546 9.920 52.329 9.076 21.199 2.218 12.965 1.529 260.182 s(β1) 3.450 0.093 0.258 0.878 0.163 0.022 0.113 0.029 0.192 0.011 0.077 0.003 0.021 0.001 1.543 t0.05×s(β1) 43.810 1.181 3.281 11.149 2.070 0.281 1.439 0.369 2.440 0.144 0.975 0.040 0.272 0.016 19.596 HS-2 X(×10-6) 951.5 42.5 230.1 852.4 108.0 6.4 49.5 4.9 27.6 4.0 13.1 1.2 7.3 1.0 95.6 β1 -8.136 -0.350 -1.297 -3.568 -0.991 -0.014 -0.641 -0.021 -0.352 -0.002 -0.246 0.000 -0.002 0.001 -0.613 β0 1049.118 46.672 245.682 895.204 119.924 6.618 57.181 5.163 31.821 4.029 16.090 1.187 7.293 0.939 102.917 s(β1) 10.139 0.328 1.344 4.581 0.572 0.044 0.430 0.002 0.212 0.020 0.146 0.004 0.028 0.003 0.422 t0.05×s(β1) 128.763 4.171 17.064 58.184 7.261 0.564 5.461 0.027 2.692 0.260 1.853 0.051 0.351 0.036 5.364 HS-3 X(×10-6) 1003.2 29.5 223.6 807.1 113.1 6.7 42.6 3.5 19.9 2.8 9.7 0.8 5.1 0.7 61.9 β1 -5.585 -0.128 -0.915 -2.283 -0.680 -0.004 -0.314 -0.009 -0.250 0.004 -0.141 0.003 0.017 0.002 -0.287 β0 1070.210 31.039 234.605 834.469 121.292 6.748 46.407 3.629 22.891 2.707 11.371 0.808 4.942 0.635 65.346 s(β1) 10.795 0.212 1.274 4.811 0.564 0.037 0.257 0.017 0.163 0.010 0.184 0.004 0.022 0.001 0.322 t0.05× s(β1) 137.096 2.690 16.183 61.104 7.164 0.467 3.266 0.219 2.070 0.132 2.342 0.054 0.276 0.015 4.084 BZ-1 X(×10-6) 94.6 9.7 16.7 57.1 10.5 2.1 11.8 1.9 11.2 2.1 5.4 0.6 3.0 0.4 70.9 β1 -0.715 -0.019 -0.074 -0.199 -0.042 -0.007 -0.052 -0.002 -0.028 -0.002 -0.009 0.000 -0.002 0.000 -0.323 β0 103.195 9.978 17.554 59.493 11.049 2.228 12.386 1.908 11.508 2.094 5.480 0.582 3.054 0.371 74.769 s(β1) 0.841 0.048 0.133 0.376 0.062 0.010 0.074 0.007 0.063 0.007 0.028 0.002 0.008 0.000 0.467 t0.05× s(β1) 10.675 0.613 1.693 4.770 0.793 0.126 0.936 0.094 0.798 0.091 0.357 0.026 0.106 0.004 5.931 BZ-2 X(×10-6) 174.1 6.3 31.6 99.4 17.3 3.3 16.6 2.3 14.3 2.5 6.9 0.8 4.4 0.5 81.5 β1 -1.698 -0.025 -0.324 -1.030 -0.209 -0.020 -0.188 -0.007 -0.147 -0.005 -0.044 -0.003 -0.010 0.000 -0.565 β0 194.453 6.578 35.461 111.716 19.855 3.508 18.866 2.402 16.046 2.612 7.429 0.838 4.561 0.548 88.259 s(β1) 1.858 0.011 0.408 0.570 0.071 0.019 0.079 0.014 0.056 0.012 0.044 0.003 0.019 0.003 0.577 t0.05×s(β1) 23.594 0.146 5.187 7.235 0.904 0.243 1.001 0.180 0.714 0.147 0.563 0.042 0.236 0.034 7.326 BZ-3 X(×10-6) 278.5 14.7 57.4 195.3 29.4 4.3 18.1 2.0 11.2 2.0 5.7 0.6 3.5 0.4 49.5 β1 -3.127 -0.089 -0.536 -1.787 -0.357 -0.035 -0.248 -0.005 -0.114 -0.009 -0.065 -0.003 -0.017 -0.001 -0.602 β0 316.050 15.804 63.866 216.697 33.727 4.730 21.030 2.015 12.580 2.080 6.513 0.657 3.752 0.439 56.739 s(β1) 2.035 0.154 0.470 1.454 0.205 0.037 0.164 0.025 0.124 0.019 0.086 0.006 0.035 0.003 0.055 t0.05× s(β1) 25.847 1.962 5.968 18.461 2.608 0.469 2.083 0.323 1.581 0.240 1.086 0.071 0.450 0.040 0.703 CJ-1 X(×10-6) 49.9 1.8 17.3 80.3 52.7 0.4 63.8 9.0 49.3 9.4 27.1 3.9 24.4 3.6 267.8 β1 -0.572 0.012 -0.101 -0.570 -0.458 -0.002 -0.491 -0.035 -0.224 -0.036 -0.062 -0.010 -0.060 -0.005 -1.816 β0 56.761 1.603 18.555 87.150 58.201 0.413 69.727 9.384 52.005 9.868 27.808 4.047 25.131 3.708 289.626 s(β1) 0.534 0.002 0.127 0.577 0.318 0.002 0.513 0.048 0.320 0.030 0.157 0.018 0.118 0.014 1.871 t0.05×s(β1) 6.787 0.029 1.608 7.324 4.035 0.029 6.516 0.607 4.070 0.383 1.995 0.230 1.497 0.174 23.760 CJ-2 X(×10-6) 472.7 16.9 96.4 336.7 61.9 5.3 49.9 5.8 28.2 4.3 10.0 0.9 5.3 0.7 106.4 β1 -4.891 -0.155 -0.918 -2.824 -0.690 -0.034 -0.537 -0.046 -0.279 -0.034 -0.107 -0.004 -0.025 -0.001 -0.670 β0 531.381 18.774 107.425 370.632 70.190 5.705 56.393 6.394 31.566 4.713 11.255 0.972 5.648 0.722 114.396 s(β1) 5.422 0.174 0.593 2.287 0.276 0.043 0.273 0.018 0.217 0.019 0.113 0.006 0.035 0.005 0.118 t0.05×s(β1) 68.865 2.208 7.529 29.049 3.510 0.545 3.470 0.230 2.760 0.236 1.441 0.081 0.440 0.059 1.501 CJ-3 X(×10-6) 382.8 13.8 75.1 261.7 55.8 7.1 60.4 9.7 57.3 11.0 30.8 4.1 25.0 3.4 356.6 β1 -3.982 -0.069 -0.609 -1.678 -0.588 -0.043 -0.750 -0.085 -0.495 -0.096 -0.282 -0.030 -0.191 -0.006 -3.210 β0 430.630 14.605 82.402 281.826 62.839 7.563 69.410 10.676 63.229 12.143 34.213 4.450 27.321 3.466 395.093 s(β1) 4.486 0.123 0.507 1.730 0.329 0.056 0.421 0.053 0.359 0.037 0.221 0.015 0.140 0.023 2.069 t0.05×s(β1) 56.978 1.568 6.435 21.973 4.172 0.708 5.352 0.675 4.557 0.471 2.803 0.191 1.783 0.288 26.276 注：x为浓度平均值；β1和β0为回归系数；s(β1)为β1的标准偏差；t0.05为自由度0.05下的t因子。
下载: 导出CSV
表 3 监控样定值结果和不确定度
Table 3 Certified values and uncertainty for monitoring samples
定值元素 $\overline{\overline X} = \pm {U_{{\rm{CRM}}}} $ (×10-6)HS-1 HS-2 HS-3 BZ-1 BZ-2 BZ-3 CJ-1 CJ-2 CJ-3 La 453.9±33.2 507.3±47.3 542.3±48.3 79.3±8.2 152.8±15.1 242.3±24.5 35.7±5.1 429.5±17.6 341.4±28.7 Ce 4.4±0.5 29.2±3.9 28.0±1.8 9.6±0.3 6.2±0.1 13.9±0.2 1.5±0.1 16.3±0.5 10.7±1.5 Pr 83.0±14.6 260.1±24.8 245.8±12.7 12.8±1.8 27.2±2.9 48.5±8.3 15.8±2.0 65.9±4.2 68.4±11.3 Nd 321.4±51.8 754.2±36 720.0±12 57.6±5.5 93.4±13.0 206.1±3.7 75.2±10.9 308.5±35.2 213.9±32.9 Sm 50.8±9.6 90.8±12.8 114.3±24.5 10.5±1.1 16.5±1.8 23.7±4.1 51.0±4.9 53.6±6.6 45.7±6.9 Eu 6.7±0.8 6.9±0.6 6.7±1.0 1.9±0.1 2.7±0.2 4.1±0.4 0.4±0.02 5.4±0.2 6.6±0.7 Gd 64.1±16.6 48.7±7.3 76.1±9.1 9.8±1.1 15.2±1.9 20.3±2.7 57.7±7.0 45.6±6.4 59.1±10.1 Tb 6.8±1.6 6.1±0.4 4.7±1.1 1.9±0.1 2.2±0.3 2.3±0.1 8.7±1.0 5.1±0.8 7.8±1.3 Dy 41.2±8.0 22.6±3.7 19.4±2.9 9.3±1.5 14.1±0.4 11.1±0.2 53.5±5.4 28.4±3.0 57.8±7.4 Ho 5.1±0.9 4.2±0.3 3.1±0.3 2.0±0.1 2.3±0.3 1.9±0.1 7.5±1.2 4.3±0.4 11.3±1.4 Er 12.6±1.9 9.7±1.4 8.7±1.7 5.1±0.3 7.2±0.4 5.1±0.6 21.8±4.0 9.3±0.7 29.4±4.3 Tm 1.9±0.4 1.3±0.1 0.7±0.1 0.7±0.01 0.8±0.1 0.4±0.1 3.3±0.6 1.1±0.1 3.8±0.7 Yb 7.3±1.5 6.5±1.4 5.1±1.1 2.9±0.3 4.6±0.4 3.2±0.3 21.7±4.5 5.5±0.7 22.4±4.7 Lu 1.6±0.1 1.3±0.1 0.5±0.1 0.5±0.03 0.5±0.07 0.3±0.06 3.9±0.5 0.7±0.1 3.5±0.4 Y 187.1±35.4 94.3±15.4 52.8±8.0 62.1±4.2 70.9±8.6 43.9±7.5 182.4±34.7 96.6±11.1 316.7±49.6
下载: 导出CSV
-
Xiao Y F, Feng Z Y, Hu G H, et al.Reduction leaching of rare earth from ion-adsorption type rare earths ore with ferrous sulfate[J].Journal of Rare Earths, 2016, 34(9):917-923.
Georgiana A M, Vladimiros G P.Leaching of lanthanides from various weathered elution deposited ores[J].Canadian Metallurgical Quarterly, 2013, 52(3):257-264. doi: 10.1179/1879139513Y.0000000060
赵芝, 王登红, 王成辉, 等.离子吸附型稀土找矿及研究新进展[J].地质学报, 2019, 93(6):1454-1465. Zhao Z, Wang D H, Wang C H, et al.Progress in prospecting and research of ion-adsorption type REE deposits[J].Acta Geologica Sinica, 2019, 93(6):1454-1465.
程丽娅.离子型稀土矿中离子稀土的ICP-AES测定方法研究[J].安徽地质, 2017, 27(2):147-149, 160. http://www.cnki.com.cn/Article/CJFDTOTAL-AHDZ201702017.htm Cheng L Y.Study on determination of rare earth ions in ion-adsorbed rare earth mineral using ICP-MS[J].Geology of Anhui, 2017, 27(2):147-149, 160. http://www.cnki.com.cn/Article/CJFDTOTAL-AHDZ201702017.htm
肖燕飞, 黄莉, 李明来, 等.ICP-AES法测定离子吸附型稀土矿镁盐体系稀土浸出液中稀土与非稀土杂质[J].稀有金属, 2017, 41(4):390-397. http://d.wanfangdata.com.cn/Periodical/xyjs201704009 Xiao Y F, Huang L, Li M L, et al.Determination of rare earth elements and non-rare earth elements in leaching solution with magnesium salt system of ion-adsorption type rare earth ore by ICP-AES[J].Chinese Journal of Rare Metals, 2017, 41(4):390-397. http://d.wanfangdata.com.cn/Periodical/xyjs201704009
Yang X J, Lin A J, Li X L, et al.China's ion-adsorption rare earth resources, mining consequences and preservation[J].Environmental Development, 2013, 8:131-136.
Georgiana A M, Vladimiros G P.Recovery of rare earth elements adsorbed on clay minerals:Ⅱ.Leaching with ammonium sulfate[J].Hydrometallurgy, 2013, 131:158-166. http://www.sciencedirect.com/science/article/pii/S0304386X1200240X
Xiao Y F, Gao G H, Huang L, et al.A discussion on the leaching process of the ion-adsorption type rare earth ore with the electrical double layer model[J].Minerals Engineering, 2018, 120:35-43. http://www.sciencedirect.com/science/article/pii/S0892687518300876
Sanjukta A K, Shailaja P P, Sangita D K.Determination of rare earth elements in Indian kimberlite using inductively coupled plasma mass spectrometer (ICP-MS)[J].Journal of Radioanalytical and Nuclear Chemistry, 2012, 294:419-424. doi: 10.1007/s10967-012-1704-3
代小吕, 赵金宝, 贺颖婷, 等.电感耦合等离子体质谱(ICP-MS)法测定离子吸附型稀土矿中的浸出稀土元素[J].中国无机分析化学, 2015, 5(4):35-40. http://www.cqvip.com/QK/60072X/201504/666946629.html Dai X L, Zhao J B, He Y T, et al.Determination of leaching rare earth elements in ion adsorption type rare earth ore by ICP-MS[J].Chinese Journal of Inorganic Analytical Chemistry, 2015, 5(4):35-40. http://www.cqvip.com/QK/60072X/201504/666946629.html
Raunt N M, Huang L S, Lin K C, et al.Uncertainty propagation through correction methodology for the determination of rare earth elements by quadrupole based inductively coupled plasma mass spectrometry[J].Analytica Chimica Acta, 2005, 530(1):91-103.
刘贵磊, 许俊玉, 温宏利, 等.动态反应池-电感耦合等离子体质谱法精确测定配分差异显著的重稀土元素[J].桂林理工大学学报, 2016, 36(1):176-183. Liu G L, Xu J Y, Wen H L, et al.Determination of heavy rare earth elements of special rare earth ores by inductively coupled plasma mass spectrometry with a dynamic reaction cell[J].Journal of Guilin University of Technology, 2016, 36(1):176-183.
Botha A, Ellison S, Linsinger T, et al.Outline for the revision of ISO Guide 35[J].Accreditation and Quality Assurance, 2013, 18:115-118. doi: 10.1007/s00769-012-0940-0
Linsinger J T P, Pauwels J, van der Veen A M H, et al.Homogeneity and stability of reference materials[J].Accreditation and Quality Assurance, 2001, 6(1):20-25. doi: 10.1007/s007690000261
王晓红, 王毅民, 高玉淑, 等.地质标准物质均匀性检验方法评价与探讨[J].岩矿测试, 2010, 29(6):735-741. doi: 10.15898/j.cnki.11-2131/td.201907310115 Wang X H, Wang Y M, Gao Y S, et al.A view on homogeneity testing techniques for geochemical reference material in China[J].Rock and Mineral Analysis, 2010, 29(6):735-741. doi: 10.15898/j.cnki.11-2131/td.201907310115
van der Veen A M H, Linsinger P T, Pauwels J.Uncer-tainty calculation in the certification of reference materials.2.Homogeneity study[J].Accreditation and Quality Assurance, 2001, 6(1):26-30. doi: 10.1007/s007690000292
杜烨, 王春龙, 刘俊保, 等.GBW(E)130573a标准物质复制均匀性与稳定性评价[J].合成材料老化与应用, 2020, 49(2):44-46, 127. Du Y, Wang C L, Liu J B, et al.Evaluation of uniformity and stability of reference material of GBW(E)130573a[J].Synthetic Materials Aging and Application, 2020, 49(2):44-46, 127.
白玉洁, 李红亮, 李微微, 等.总有机碳分析仪校准用标准物质的研制[J].工业计量, 2020, 30(2):80-82. Bai Y J, Li H L, Li W W, et al.Development of certified reference material for calibration of total organic carbon analyser[J].Industrial Measurement, 2020, 30(2):80-82.
曾美云, 刘金, 邵鑫, 等.磷矿石化学成分分析标准物质研制[J].岩矿测试, 2017, 36(6):633-640. doi: 10.15898/j.cnki.11-2131/td.201907310115 Zeng M Y, Liu J, Shao X, et al.Preparation of phosphate ore reference materials for chemical composition analysis[J].Rock and Mineral Analysis, 2017, 36(6):633-640. doi: 10.15898/j.cnki.11-2131/td.201907310115
杨理勤, 陈占生, 谢璐, 等.卡林型金矿砷成分分析标准物质研制[J].岩矿测试, 2018, 37(2):209-216. doi: 10.15898/j.cnki.11-2131/td.201907310115 Yang L Q, Chen Z S, Xie L, et al.Preparation of gold and arsenic certified reference materials for chemical composition analysis in Carlin-type gold deposits[J].Rock and Mineral Analysis, 2018, 37(2):209-216. doi: 10.15898/j.cnki.11-2131/td.201907310115
金浩, 韩永志.标准物质及其应用技术[M].北京:中国标准出版社, 2003:34-35. Jin H, Han Y Z.Standard reference material and its application technology[M].Beijing:China Standards Press, 2003:34-35.
宋丽华, 郝原方, 杨柳, 等.地质标准物质的研制方法[J].地质与资源, 2013(5):419-421. Song L H, Hao Y F, Yang L, et al.Preparation method of geochemical reference materials[J].Geology and Resources, 2013(5):419-421.
辛文彩, 夏宁, 徐磊, 等.长江三角洲沉积物标准物质研制[J].岩矿测试, 2017, 36(4):388-395. doi: 10.15898/j.cnki.11-2131/td.201907310115 Xin W C, Xia N, Xu L, et al.Preparation of Yangtze River Delta sediment reference materials[J].Rock and Mineral Analysis, 2017, 36(4):388-395. doi: 10.15898/j.cnki.11-2131/td.201907310115
田衎, 杨珺, 孙自杰, 等.矿区污染场地土壤重金属元素分析标准样品的研制[J].岩矿测试, 2017, 36(1):82-88. doi: 10.15898/j.cnki.11-2131/td.201907310115 Tian K, Yang J, Sun Z J, et al.Preparation of soil certified reference materials for heavy metals in contaminated sites[J].Rock and Mineral Analysis, 2017, 36(1):82-88. doi: 10.15898/j.cnki.11-2131/td.201907310115
高捷, 盛成, 卓尚军.X射线荧光光谱分析用的含铁尘泥标准样品的研制[J].冶金分析, 2015, 35(2):74-78. Gao J, Sheng C, Zhuo S J.Development of ferric-containing dust and sludge standard samples used for X-ray fluorescence spectrometric analysis[J].Metallurgical Analysis, 2015, 35(2):74-78.
汪斌, 卢晓华, 王茜.质量控制图在标准物质稳定性评估中的应用[J].化学试剂, 2019, 41(5):475-477. Wang B, Lu X H, Wang Q.Application of control chart for assessment of stability of reference materials[J].Huaxue Shiji, 2019, 41(5):475-477.
李津, 马健雄, 闫斌, 等.黑色页岩铁同位素标准物质的研制[J].地球学报, 2020, 41(5):623-629. http://www.cqvip.com/QK/98325A/202005/7102942635.html Li J, Ma J X, Yan B, et al.The preparation of reference material for Fe isotope measurement of black shale samples[J].Acta Geoscientica Sinica, 2020, 41(5):623-629. http://www.cqvip.com/QK/98325A/202005/7102942635.html
胡德龙, 张雯, 陈家颖.二氧化碳中一氧化氮气体标准物质研制[J].化学分析计量, 2020, 29(2):8-11. Hu D L, Zhang W, Chen J Y. Preparation for gas reference material of nitric oxide in carbon dioxide[J].Chemical Analysis and Meterage, 2020, 29(2):8-11.
刘妹, 顾铁新, 潘含江, 等.泛滥平原沉积物标准物质研制[J].岩矿测试, 2018, 37(5):558-571. doi: 10.15898/j.cnki.11-2131/td.201907310115 Liu M, Gu T X, Pan H J, et al.Preparation of seven certified reference materials for floodplain sediments[J].Rock and Mineral Analysis, 2018, 37(5):558-571. doi: 10.15898/j.cnki.11-2131/td.201907310115
张磊, 周伟, 朱云, 等.硫酸铵溶液淋滤-电感耦合等离子体质谱测定离子相稀土分量的方法优化[J].岩矿测试, 2018, 37(5):518-525. doi: 10.15898/j.cnki.11-2131/td.201907310115 Zhang L, Zhou W, Zhu Y, et al.An optimized method for determination of ionic-phase rare earth elements by ICP-MS using ammonium sulfate leaching[J].Rock and Mineral Analysis, 2018, 37(5):518-525. doi: 10.15898/j.cnki.11-2131/td.201907310115
-
期刊类型引用(30)
1. 范晨子，孙冬阳，赵令浩，袁继海，胡明月，赵明. 激光剥蚀电感耦合等离子体质谱法微区原位定量分析锂铍矿物化学成分. 岩矿测试. 2024(01): 87-100 . 
本站查看
2. 周伟，张嘉升，祁晓鹏，徐磊，杨杰. X射线衍射和TIMA研究陕南镇巴地区富锂黏土岩的矿物组成及锂的赋存状态. 岩矿测试. 2024(01): 76-86 . 本站查看
3. 苑东洋，石威，刘德民，黄天统，徐龙生，陈超，孔令耀. 鄂北耿集-板桥地区新元古代变火山岩系地层归属——来自锆石U-Pb测年、Hf同位素的证据. 吉林大学学报(地球科学版). 2024(02): 516-530 . 百度学术
4. 叶健华，方洁，钟石玉，杨成，张凯玲，周梦阁. 桐柏—大别造山带中部娘娘顶花岗岩成因——锆石U-Pb年龄、岩石地球化学及Sr-Nd同位素证据. 资源环境与工程. 2024(04): 381-392 . 百度学术
5. 李国占，涂家润，崔玉荣，毕君辉，郝爽，刘文刚，周红英. NEPTUNE型质谱仪等离子体源的维护与故障分析. 华北地质. 2024(04): 98-105 . 百度学术
6. 赵峰，王占明，王挺，龙绍军. 激光剥蚀电感耦合等离子体质谱(LA-ICP-MS)法测定涂层氧化锆颗粒铌层中24种微量杂质元素. 中国无机分析化学. 2023(04): 374-381 . 百度学术
7. 胡子奇，张德贤，刘磊. 束斑直径和能量密度对锆石U-Pb定年准确度的影响研究. 现代地质. 2023(03): 722-732 . 百度学术
8. 童铄云，潘诗洋，胡伟康，刘爽，董学林，鲁力. 激光剥蚀-电感耦合等离子体质谱(LA-ICP-MS)法测定底栖有孔虫中的元素/钙比值. 中国无机分析化学. 2023(09): 967-974 . 百度学术
9. 刘应冬，承道明. 青海大场金矿矿集区成矿规律与找矿潜力. 矿产综合利用. 2022(01): 7-14 . 百度学术
10. 钟石玉，徐浩，汤旋，胡太平，杨成，王小敏，阮业东，黄丽. 黄梅县婆婆寨岩体锆石U-Pb年代学及岩石地球化学特征研究. 资源环境与工程. 2022(04): 453-462 . 百度学术
11. 谭细娟，郭超，凤永刚，周义，梁婷. 激光剥蚀系统气体流速变化对LA-ICP-MS锆石U-Pb定年精度的影响. 岩矿测试. 2022(04): 554-563 . 本站查看
12. 王森，张拴宏，张琪琪，梁霞，孔令昊，胡国辉，王开，王宏宇，蔡瑜杭. LA-ICP-MS锆石微区U-Pb定年方法及不同束斑直径对年龄结果的影响作用探讨. 地质力学学报. 2022(04): 642-652 . 百度学术
13. 陈玺贇，吕开来，王平，黄湘通，孔兴功. 大样本量(large-n)碎屑锆石U-Pb年代学分析技术研究进展. 岩矿测试. 2022(06): 920-934 . 本站查看
14. 孔令耀，杨成，李琳静，周峰，郭盼，陈超. 扬子北缘南华系武当群锆石U-Pb年龄：对扬子北缘南华纪裂解模式的启示. 地质通报. 2022(12): 2224-2243 . 百度学术
15. 崔玉荣，涂家润，李国占，周红英，耿建珍. LA-ICP-MS榍石U-Pb定年方法. 华北地质. 2022(04): 53-59 . 百度学术
16. 钟石玉，杨成，李书涛，石先滨，胡太平，秦志军. 大悟县娘娘顶钨铍矿区隐伏中性岩体锆石U-Pb年龄、地球化学特征及地质意义. 资源环境与工程. 2021(01): 20-29+52 . 百度学术
17. 黄亮，吕勃烨，王晓林，丛峰，黄晓明，熊波，宋冬虎. 滇西漕涧地区两类晚奥陶世强过铝花岗岩锆石U-Pb年龄、地球化学和Hf同位素特征及对泛非运动的指示. 地质通报. 2021(07): 1033-1046 . 百度学术
18. 郭盼，陈超，吴波，陈小龙，孔令耀，刘成新. 西大别～2.0 Ga变质花岗岩的发现及其Hf同位素特征. 中国地质. 2021(04): 1267-1279 . 百度学术
19. 靳梦琪，李艳广，王鹏，汪双双，黎卫亮. 榍石LA-ICP-MS U-Pb定年中元素分馏的影响及校正研究. 岩矿测试. 2020(02): 274-284 . 本站查看
20. 李阳，邹灏，刘行，蒋修未，李蝶. SILLS软件在单个萤石流体包裹体LA-ICP-MS微量元素分析数据处理中的应用. 岩矿测试. 2020(02): 300-310 . 本站查看
21. 陈超，苑金玲，郭盼，黎方云，孔令耀，杨金香，毛新武. 扬子陆块～2.0 Ga的区域变质事件对南北黄陵古元古代差异演化的启示. 中国地质. 2020(04): 899-913 . 百度学术
22. 万秋，李延河，王利民，段超，施珂，谭德兴. 北淮阳晓天火山岩盆地片麻状花岗岩成岩年代学及地球化学特征. 岩矿测试. 2020(04): 620-630 . 本站查看
23. 黄亮，吴嘉林，王晓林，熊波，浦涛，朱勋早，方雄. 滇西漕涧地区崇山造山带伟晶岩的发现与地质意义. 矿物学报. 2020(05): 584-594 . 百度学术
24. 钟石玉，杨成，李书涛，石先滨，胡太平，秦志军. 蕲春太阳脑岩体锆石U-Pb年代学及岩石地球化学特征研究. 资源环境与工程. 2020(S2): 8-15+77 . 百度学术
25. 谢纪海，胡正祥，毛新武，孔令耀，杨青雄，杨成，郭盼. 鄂北大洪山晋宁期MORB-like玄武岩的识别与洋内俯冲作用. 中国地质. 2019(06): 1496-1511 . 百度学术
26. 魏均启，桂博艺，朱丹，王芳，鲁力，潘诗洋. 激光剥蚀电感耦合等离子体质谱联用工作参数优化及在黄铁矿原位分析中的应用. 资源环境与工程. 2019(04): 587-592 . 百度学术
27. 王辉，汪方跃，关炳庭，盛兆秋. 激光能量密度对LA-ICP-MS分析数据质量的影响研究. 岩矿测试. 2019(06): 609-619 . 本站查看
28. 陈超，苑金玲，孔令耀，叶竹君，杨青雄，杨成，周峰. 扬子北缘大洪山地区早古生代基性岩脉的厘定及其地质意义. 地球科学. 2018(07): 2370-2388 . 百度学术
29. 周向科，毛冰，祁小军，谭锴，万海洋. 金沙江-红河成矿带北衙金多金属矿床成岩成矿时代研究. 中国矿业. 2018(S1): 203-210 . 百度学术
30. 黄新鹏. 福建霞浦大湾钼铍矿区碱长花岗岩LA-ICP-MS锆石U-Pb测年研究. 岩矿测试. 2018(05): 572-579 . 本站查看
其他类型引用(3)
计量
- 文章访问数: 3492
- HTML全文浏览量: 1312
- PDF下载量: 48
- 被引次数: 33
京公网安备 11010202008159号