Quantitative Analysis of Mineral Composition of Kyanite by X-ray Diffraction with Rietveld Refinement Method
-
摘要: 蓝晶石矿物定量分析通常采用的化学物相法分析流程繁琐,易受同质异象矿物和难熔矿物的干扰;采用X射线衍射内标法需要提纯矿物绘制标准曲线,但矿物包体的存在使得获取纯净单矿物成为难题;由于不同矿区蓝晶石矿物成分之间存在着差异性,以上两种方法都仅仅适用于同一矿区样品的矿物定量分析。为简化分析流程,提高测试效率,本文采用X射线衍射(XRD)全谱拟合法对蓝晶石国家二级标准物质和野外样品进行分析,并与RIR法、绝热法和K值法等衍射定量方法及化学分析结果进行了比较验证方法的可靠性。结果表明:全谱拟合法无需采用标准物质,也不用引入内标物,一次扫谱分析能获得样品中所有成分的信息,操作简单,能有效降低择优取向对衍射定量结果的影响;对于含量大于5%的矿物,组分分析结果的绝对误差小于1%,显著低于绝对误差允许限;采用Highscore和Jade两种衍射数据处理软件获得的定量结果吻合度高,目标矿物蓝晶石分析结果的双差均小于0.8%,相对偏差小于2.5%;蓝晶石的加标回收率在95.3%~101.0%之间,能实现不同矿区蓝晶石矿的快速定量分析。要点
(1) 采用X射线衍射全谱拟合法对蓝晶石矿标准物质和野外样品进行矿物成分定量分析。
(2) 全谱拟合法比RIR法、绝热法和K值法等衍射定量分析方法具有更高的准确度。
(3) 本方法适用于不同矿区蓝晶石矿的快速定量分析。
HIGHLIGHTS(1) Application of the Rietveld refinement method for quantification of mineral concentrations in kyanite was performed.
(2) The accuracy of the method was superior to the other methods of X-ray quantitative phase analysis, such as RIR method, adiabatic method, K-value method.
(3) The method can be used to quickly determine the contents of all minerals in kyanite from different mining area.
Abstract:BACKGROUNDThe content of minerals in kyanite can be analyzed by chemical phase method, but the analytical process is very tedious. Moreover, the existence of heteromorphism and refractory minerals can affect the accuracy of analysis results. An internal standard method by X-ray diffraction (XRD) for direct determination of kyanite requires the use of pure mineral to draw a standard curve, but it is extremely challenging to purify minerals due to the inclusion. The differences between the mineral components of kyanite in different mining areas makes the above two methods applicable only to the quantitative analysis of minerals in the same mining area.OBJECTIVESTo simplify analysis the process of mineral composition of kyanite and improve efficiency.METHODSThe content of minerals in kyanite was studied by X-ray diffraction. The X-ray diffraction Rietveld refinement method was used to analyze the secondary standard materials and field samples of kyanite, and the results were compared with the diffraction quantitative method and chemical analysis results including RIR method, adiabatic method and K-value method.RESULTSRietveld refinement method was simple and effective in correcting the diffraction intensity error caused by the preferred orientation. Accurate results can be obtained by Rietveld refinement method more than other methods of X-ray quantitative phase analysis, such as RIR method, adiabatic method, K-value method. The analytic absolute error of minerals with content greater than 5% was less 1%, and significantly lower than the allowable error. The results obtained by Highscore agree well with results measured by Jade, and the double difference of results was less than 0.8%. The relative standard deviation was less than 2.5%. The spiked recoveries of kyanite were 95.3%-101.0%.CONCLUSIONSX-ray diffraction with Rietveld refinement method is simple, and can be used to determine the content of all minerals in kyanite from different mining areas.-
Keywords:
- kyanite /
- phase composition /
- X-ray diffraction /
- Rietveld refinement method /
- preferred orientation
-
土壤中总有机碳(TOC)是指土壤中与有机质相关的碳元素总称,是土壤理化性质重要指标之一,主要用于指示土壤中有机质含量。研究发现,土壤中有机碳不仅在全球碳平衡中发挥着至关重要的作用,其含量与土壤全氮、阳离子交换量、pH值等土壤理化性质之间也具有相关性,而且还对有机污染物特别是多环芳烃(PAHs)在土壤中的迁移和吸附有重要影响[1-5]。土壤中总有机碳被认为是既能定量又能定性反映土壤信息的综合性参数,因此被广泛应用于土壤状况的测定与评估[3, 6]。而土壤中总有机碳标准样品作为配套土壤中总有机碳测定标准方法的实物标准,可有效保证土壤总有机碳监测数据的质量,对于客观评价土壤健康水平、优化土壤质量管理等工作具有重要意义[7-8]。
表 1列出了国内外主要的土壤中有机碳标准样品研究机构,以及所开发样品的主要信息。国外对土壤中有机碳的研究始于20世纪60年代。为满足对土壤有机碳存量等相关研究的需要,欧盟标准样品和测量研究所(IRMM)、美国RTC公司等国际标准样品研究机构分别针对黏土、砂土等不同质地的土壤,陆续开发出指标为总碳(TC)或总有机碳(TOC)的土壤标准样品。这些国外土壤标准样品中总有机碳的含量范围介于4.00~60.0mg/g之间(表 1)。多年来中国学者在土壤有机碳储量、存在形态、空间分布特征等领域的研究取得了十分丰硕的成果,但对土壤中有机碳标准样品的研究方面关注程度较为有限[9-13]。中国地质科学院地球物理地球化学勘查研究所(IGGE)研制的土壤中元素成分分析标准样品提供了土壤中总碳及有机碳的量值信息,可用于土壤肥力或土壤营养成分的评价等工作,但中国一直未曾推出专门用于土壤中有机碳分析检测的环境标准样品。2021年初,国家地质实验测试中心(NRCG)发布了土壤碳形态系列标准样品(GBW07974~GBW07977),量值指标包括总碳、有机碳以及重组有机碳。而上述中国土壤样品进行量值评定时,主要是使用传统的重铬酸钾容量法作为测定土壤中有机碳的分析方法,暂未采用当前国际上使用较多的燃烧氧化-非分散红外法。
表 1 国内外主要土壤标准样品中涉及总有机碳的信息Table 1. Information of TOC in major soil reference materials at home and abroad土壤标准样品研制机构 土壤标准样品编号 碳含量范围(mg/g) 土壤标准样品具体类型 总碳(TC) 总有机碳(TOC) 欧盟标准样品和测量研究所 IRMM-443-1~ IRMM-443-7 14.5~108.1 13.1~59.6 欧洲土壤 美国RTC公司 RTC-CCRM090~RTC-CCRM092 - 4.13~5.72 黏土、砂壤土、砂土 RTC-SQC014 - 6.51 国家地质实验测试中心 GBW07974~GBW07977 8.1~56.8 6.9~51.5 土壤
农用地土壤GBW07913~GBW07942 6.2~79.0 3.7~47.3 中国地质科学院地球物理地球化学勘查研究所 GBW07401~GBW07408 5.5~21.1 4.9~18.0 暗棕壤、栗钙土、黄棕壤、红壤等
主要平原土壤
盐碱土、棕漠土、灰钙土、滩涂沉积物等
土壤GBW07423~GBW07430 9.3~19.0 6.2~13.5 GBW07446~GBW07457 10.3~19.4 2.5~11.5 GBW07536~GBW07573 2.8~38.1 1.7~27.5 本文在梳理土壤中有机碳检测技术发展趋势的基础上,选择近些年逐渐流行的燃烧氧化-非分散红外法,以中国有机碳含量较高的农用地土壤为原材料,研制了土壤中总有机碳标准样品,并对样品量值及不确定度评定方法等进行了研究和讨论,为同类指标环境基体标准样品的研制提供参考。
1. 实验部分
1.1 分析仪器
总有机碳分析仪(Aurora 1030C + 1030S,美国OI Analytical公司);X射线荧光光谱仪(Axios PW4400,荷兰Malvern Panalytical公司);精密电子天平(AE240,瑞士Mettler Toledo公司)。
1.2 样品和主要试剂
葡萄糖基准试剂(CAS: 50-99-7,美国Sigma-Aldrich公司);土壤TOC质量控制样品[SQC014-100G,(5210±2280)mg/kg,美国Sigma-Aldrich公司];土壤成分分析标准样品(GBW07430,0.97%±0.12%);盐酸(MOS级,北京化学试剂研究所)。
1.3 候选物的采集与制备
考虑到土地功能、管理模式等人为因素对土壤中有机碳的巨大影响,选择受人为活动影响较少的农用地土壤作为土壤标准样品的候选物[14-16]。鉴于中国土壤标准样品中有机碳含量普遍偏低,在对中国土壤资源分布调研的基础上,确定以腐植质层深厚且钙积层较深的黑钙土为研究目标,并选定因山地气温低导致腐植质积累偏高的新疆地区为样品采集地。在新疆维吾尔自治区某牧场采集pH值为8.4的土壤原料样品,采样深度为0~20cm,共采集土壤样品约100kg。
除掉土壤中杂物等非土壤部分,在80℃低温条件下烘干样品,使样品水分含量低于1%。用陶瓷球磨机研磨样品,研磨后的样品过100目(150μm)不锈钢筛[17],样品粒度能够满足《土壤环境监测技术规范》(HJ/T 166—2004)对制样粒度的要求。收集筛下样品转移至混匀机中混匀8h,并采用X射线荧光光谱法(XRF)对混匀后样品进行均匀性初检。通过均匀性初检的土壤样品分装于40mL棕色玻璃瓶中,每瓶样品质量为30g,共计分装成品900瓶。以放射线60Co辐照灭菌,使土壤样品的菌落总数低于10CFU/g[18]。样品于室温下避光储存。
1.4 样品分析测试
1.4.1 仪器分析条件
检测模式为固体-总碳(Solids-TC),样品注射量0.200mL,系统压力20psi,排空时间15s,反应时间4min,检测时间3min,反应温度900℃,检测温度680℃,样品杯预处理时间1.5min,测量方式为标准曲线法。
1.4.2 检测样品的制备
在镂空底石英样品杯中填充少量石英棉,准确称取干燥后的待测土壤样品0.1g(精确到0.0001g)于石英样品杯中,缓慢滴加50%的盐酸约300μL,覆盖样品杯内全部样品。石英样品杯放入烘箱中120℃烘干2h,去除总无机碳(TIC)。样品妥善放置,待测。
2. 结果与讨论
2.1 测量方法的选择
土壤中总碳(TC)由总无机碳(TIC)和总有机碳(TOC)组成。土壤中总有机碳的测定原理是通过氧化的方式,将土壤中各种类型的有机碳(OC)转化为可以定量测定的化合物,利用碳元素之间的对应关系计算样品中有机碳含量。近年来,关于土壤有机碳测定方法的研究较多,按照测定过程来看,可分为土壤有机碳的氧化和氧化产物的测定两个步骤。几种常见的土壤中总有机碳的测定方法归纳于图 1。
根据土壤有机碳氧化方式的不同,可将氧化方法大致分为以下两种类型:一是使用强氧化剂作为氧化手段的化学氧化法,即湿法氧化法;二是使用干式燃烧作为氧化手段的高温氧化法,即燃烧法[3, 6, 10, 19]。湿法氧化法是指在土壤样品中加入过量的氧化剂,将样品中的有机碳全部氧化为二氧化碳,然后对氧化产物进行测定,进而换算出有机碳的含量。常用的湿法氧化法有重铬酸盐氧化法、过硫酸盐氧化法、过氧化物氧化法等[20-23]。燃烧法是将土壤样品置于燃烧炉中,使样品在高温下(通常超过900℃)的氧气流中燃烧,将样品中的碳元素全部氧化为二氧化碳。燃烧法既可先测定样品中总碳以及总无机碳,然后通过差减法间接计算出总有机碳含量;也可通过测定加酸并烘干去除总无机碳后的样品,直接获得样品中的总有机碳含量[21-22]。
对样品氧化后的产物(如:重铬酸根、二氧化碳等)进行测定,可采用滴定法、比色法、重量法、容量法、光度法、气相色谱法等多种不同技术[21-23]。手动滴定法是通过滴定过量的重铬酸根,然后换算出结果。比色法是测定溶液中反应产生的三价铬的吸光度,再间接确定指标含量。重量法是根据吸附装置在吸收二氧化碳后其质量发生的变化,从而换算出结果。容量法是通过滴定吸收二氧化碳后过量的氢氧化钡,进一步计算出有机碳含量。光度法是采用非色散红外(NDIR)检测器测定反应生成的二氧化碳的量,进而计算碳的含量。气相色谱法是使用配有热导检测器(TCD)的气相色谱仪测定二氧化碳的量,或将二氧化碳转化为甲烷后再用氢火焰离子化检测器(FID)进行测定[3, 6, 21-22]。
目前,测定土壤中有机碳含量最常用方法是重铬酸钾-外加热氧化法(湿法氧化)和总有机碳分析仪法(高温燃烧氧化法),两种氧化方法各有优缺点、适用范围及局限性[3, 10, 24]。重铬酸钾-外加热氧化法为测定土壤有机碳的传统方法,国际上采用该技术的标准方法为重铬酸钾氧化比色法《Soil quality—Determination of organic carbon by sulfochromic oxidation》(ISO 14235—1998),中国常用的标准方法为《区域地球化学样品分析方法第27部分:有机碳量测定重铬酸钾容量法》(DZ/T 0279.27—2016)。上述方法对于仪器设备要求少,成本较低,灵敏度高,不受土壤中碳酸盐的影响,用于碱性土壤的测定具有一定优势[3, 25-26]。但是,重铬酸钾-外加热氧化法存在氧化不完全的缺点,在测定高含量氯离子、亚铁离子、锰化合物的土壤时,测定结果会受到影响[27]。而且,该方法在实验中会消耗大量硫酸和重铬酸钾,产生大量废液,对环境不够友好[28-29]。相比之下,采用高温燃烧氧化法的有机碳分析仪具有更好的稳定性,对土壤中以各种形式存在的碳元素的氧化更彻底,实验测定结果的回收率较高,样品前处理工作量小,且不对环境产生污染。虽然燃烧氧化-非分散红外法的仪器昂贵且维护费用较高,但综合考虑其优势,以该方法检测土壤中有机碳的实验室越来越多[24, 28-29]。采用燃烧法的国际标准方法为《Soil quality—Determination of organic and total carbon after dry combustion (elementary analysis)》(ISO 10694—1995)和《Total organic carbon (TOC) in soil》(USEPA SW-846 Test Method 9060A)。中国的同类方法为《土壤有机碳的测定燃烧氧化-非分散红外法》(HJ 695—2014)。多篇文献曾报道关于燃烧氧化-非分散红外法与重铬酸钾容量法测定结果之间关系的研究,并对两种方法的特点进行了探讨,两种检测方法的优缺点汇总于表 2[3, 17, 25-26]。
表 2 土壤中总有机碳测定常用方法的比较Table 2. Comparison of common methods for determination of TOC in soil对比项目 燃烧氧化-非分散红外法 重铬酸钾容量法 优点 碳元素氧化完全,稳定性高,测定结果回收率高,样品前处理快捷,污染排放低 灵敏度高,操作简单,仪器设备费用低,使用成本低,更适合分析碱性土壤 缺点 仪器费用昂贵,使用维护费较高,长期使用仪器背景值易升高 存在氧化不完全的情况,易受氯离子、亚铁离子、锰化合物干扰,废液产生量大,不适合用于大量样品的快速测定 相应标准方法 《土壤有机碳的测定燃烧氧化-非分散红外法》 (HJ 695—2014) 《区域地球化学样品分析方法第27部分:有机碳量测定重铬酸钾容量法》(DZ/T 0279.27—2016) 鉴于中国尚未开发以燃烧氧化-非分散红外法用于定值的标准样品,且仪器法通常具有更好的稳定性,本研究采用燃烧氧化-非分散红外法(即总有机碳分析仪法)作为土壤中总有机碳标准样品均匀性检验、稳定性检验等分析测定的方法,采用重铬酸钾容量法和燃烧氧化-非分散红外法两种方法共同进行标准样品中总有机碳的协作定值。
2.2 均匀性检验
均匀性用来描述标准样品的空间分布特征,是表征样品量值指标一致性的重要参数。影响土壤样品均匀性的主要因素有:加工制备过程的混匀程度、研磨后粒径的均匀程度,以及分装过程中的潜在影响因素等。由于土壤样品为粉末颗粒,相比液体样品其流动性和分散性较差,因此对样品开展均匀性检验显得十分重要。
根据《标准样品工作导则(3) 标准样品定值的一般原则和统计方法》(GB/T 15000.3—2008)中的方法,采用分层随机抽样法,在土壤样品分装过程中的前、中、后期随机抽取10瓶样品,以燃烧氧化-非分散红外法检测标准样品中总有机碳指标的均匀性,每瓶样品平行测定2次。土壤中总有机碳标准样品均匀性检验数据详见表 3。由表中结果可见,样品均匀性实验的总均值为23.2mg/g,相对标准偏差(RSD)为1.4%。以单因素方差分析法对检测数据进行统计和评价[30]。统计量F的计算公式如下:
$ {F_{{\rm{cal}}}} = {\rm{M}}{{\rm{S}}_{\rm{b}}}/{\rm{M}}{{\rm{S}}_{\rm{w}}} = \frac{{{\rm{S}}{{\rm{S}}_{\rm{b}}}}}{{{v_{\rm{b}}}}}/\frac{{{\rm{S}}{{\rm{S}}_{\rm{w}}}}}{{{v_{\rm{w}}}}} $
(1) 表 3 土壤中总有机碳标准样品的均匀性检验结果Table 3. Homogeneity assessment results of TOC in soil reference material样品编号 总有机碳均匀性测定结果(mg/g) RSD(%) 平行测定1 平行测定2 均值 总均值 1 24.4 22.9 23.7 23.2 1.4 2 24.0 22.6 23.3 3 23.2 23.2 23.2 4 24.5 22.6 23.6 5 24.2 23.2 23.7 6 22.5 23.5 23.0 7 23.0 23.2 23.1 8 22.2 23.3 22.8 9 22.9 22.8 22.9 10 23.4 22.9 23.2 式(1)中:MSb为瓶间均方;MSw为瓶内均方;SSb为瓶间平方和;vb为瓶间自由度;SSw为瓶内平方和;vw为瓶内自由度。
F统计量=0.38<F0.05(9, 10)=3.02,说明在95%的置信区间内,样品瓶间及瓶内的有机碳含量无显著性差异,样品具有良好的均匀性。结合均匀性检验取样量以及测定所使用检测方法的情况,最终确定不影响本标准样品中总有机碳指标测量重复性的最小取样量为0.1g。
2.3 稳定性检验
稳定性是指在规定的时间间隔和环境条件下,标准样品的特性量值保持在规定范围内的性质[31]。在18个月稳定性检验期间,按照前密后疏的原则进行6次抽样及检测。每次随机抽取室温避光保存的土壤标准样品3瓶,采用燃烧氧化-非分散红外法进行测定。每瓶样品平行测定2次,以3瓶样品测定值的总平均值作为该时间点稳定性检验的检测值。土壤中总有机碳标准样品稳定性检验的研究结果详见表 4。由结果计算出样品稳定性实验测定的总有机碳总均值为23.4mg/g。
表 4 土壤中总有机碳标准样品的稳定性检验结果Table 4. Stability assessment results of TOC in soil reference material时间间隔(月) 总有机碳测定结果(mg/g) 时间间隔(月) 总有机碳测定结果(mg/g) 0 23.2 6 23.2 1 23.5 12 23.4 3 23.2 18 24.1 检测数据采用《标准样品工作导则(3) 标准样品定值的一般原则和统计方法》(GB/T 15000.3—2008)推荐的一元线性拟合模型进行统计检验,即假定特性量值变化与储存时间之间存在线性关系,通过考察斜率(b1)的显著性来评价样品的稳定性[32]。斜率b1的计算公式如下:
$ {b_1} = \sum\limits_{i = 1}^n {\left( {{X_i} - \bar X} \right)} \left( {{Y_i} - \bar Y} \right)/\sum\limits_{i = 1}^n {{{\left( {{X_i} - \bar X} \right)}^2}} $
(2) 斜率的标准偏差s(b1)的计算公式如下:
$ s\left( {{b_1}} \right) = \frac{{\sqrt {\frac{{\sum\limits_{i = 1}^n {{{\left( {{Y_i} - {b_0} - {b_1}{X_i}} \right)}^2}} }}{{n - 2}}} }}{{\sqrt {\sum\limits_{i = 1}^n {{{\left( {{X_i} - \bar X} \right)}^2}} } }} $
(3) 式(3)中:b0为截距,计算公式为b0=Y-b1X;n为稳定性检测次数。
采用以上公式计算,b1和s(b1)的结果分别为0.039和0.016。经查表t(0.95, n-2)=2.78,斜率|b1|小于t(1-α, n-2)×s(b1),表明在18个月的检测期内,在室温避光条件下,样品中有机碳含量随时间没有发生明显变化,即标准样品的稳定性良好。
2.4 定值分析与不确定度评定
2.4.1 定值方法
采用多家实验室协作定值的方式进行共同测定[18, 33-34]。邀请了中国在土壤样品检测方面能力较强的11家实验室,分别为:生态环境部环境发展中心环境标准样品研究所、天津市生态环境监测中心、国家地质实验测试中心、广西壮族自治区地质矿产测试研究中心、陕西省地质矿产实验研究所有限公司、河南省岩石矿物测试中心、沈阳矿产资源监督检测中心、长春矿产资源监督检测中心、成都矿产资源监督检测中心、武汉矿产资源监督检测中心、中国科学院新疆生态与地理研究所。分别向参加协作定值的每家实验室提供土壤样品6瓶,每家实验室反馈6个独立数据结果。参加定值的11家协作实验室共采用了2种分析方法:采用燃烧氧化-非分散红外法的实验室为5家,采用重铬酸钾容量法的实验室为7家,其中1家实验室同时采用了以上两种测量方法。各家实验室所采用的分析方法及测定结果汇总于表 5。
表 5 协作实验室所用分析方法和测定数据Table 5. Analysis methods and measurement data of the collaborative laboratories实验室代号 总有机碳分析方法 总有机碳协作测定数据(mg/g) RSD(%) S-1 S-2 S-3 S-4 S-5 S-6 均值 Lab-1 燃烧氧化-非分散红外法 23.7 23.3 23.2 23.6 23.7 23.0 23.4 1.3 Lab-2 燃烧氧化-非分散红外法 23.4 22.8 22.8 23.8 22.9 22.4 23.0 2.2 Lab-3 重铬酸钾容量法 26.7 26.4 26.4 26.3 26.5 26.4 26.5 0.6 Lab-4a 重铬酸钾容量法 25.4 25.6 26.1 26.0 26.0 25.4 25.8 1.3 Lab-4b 燃烧氧化-非分散红外法 25.9 25.9 25.5 25.7 25.7 26.1 25.8 0.9 Lab-5 重铬酸钾容量法 26.3 24.7 25.1 27.8 26.7 27.3 26.3 4.7 Lab-6 重铬酸钾容量法 25.8 26.3 26.5 25.9 26.3 26.0 26.1 1.1 Lab-7 重铬酸钾容量法 25.2 25.3 25.4 25.3 25.2 25.0 25.2 0.6 Lab-8 燃烧氧化-非分散红外法 21.1 21.1 21.5 21.1 21.3 21.3 21.2 0.8 Lab-9 重铬酸钾容量法 27.1 27.1 28.1 28.4 28.2 28.2 27.9 2.2 Lab-10 重铬酸钾容量法 26.1 26.0 25.4 26.3 24.9 26.0 25.8 2.1 Lab-11 燃烧氧化-非分散红外法 24.9 24.8 25.4 25.4 25.1 25.7 25.2 1.4 2.4.2 数据统计处理
在全部测定数据中,每家协作实验室提供的每个测定结果均满足溯源性要求。采用燃烧氧化-非分散红外法的5组数据的总均值为23.7mg/g,均值标准偏差为1.9mg/g;而采用重铬酸钾容量法的7组数据的总均值为26.2mg/g,均值标准偏差为0.9mg/g。上述两种方法检测结果的不确定度分量之间具有一致性,并且对样品特性值不确定度的贡献程度均较小。国内外多个课题组就燃烧氧化-非分散红外法与重铬酸钾容量法测定结果之间的关系进行了研究,所获得的结论是,在一定的含量范围内,上述两种方法的测定结果之间具有相关性[3, 25-26]。考虑到对于环境基体样品,单一分析方法可能存在基体影响或基体干扰的风险,两种方法测定结果的中位值将是更加稳健的统计量,用于评估标准样品的特性量值会更合适。
据此,将全部数据合并为一个总体进行技术审查,以Grubbs准则检验每家实验室内平行数据的一致性,以Dixon准则检验实验室间数据平均值的一致性,以Cochran准则检验实验室间标准偏差的一致性,最后采用偏态系数和峰态系数检验全部有效数据的正态分布[18, 33]。收回的12组72个检测结果全部通过统计检验,数据采用率为100%。当数据分布为正态分布或近似正态分布时,计算协作测定的总均值(
$\overline{\overline X} $ )、总均值标准偏差(s)。经统计分析,标准样品协作测定的总均值为25.2mg/g,均值标准偏差为1.8mg/g。采用燃烧氧化-非分散红外法定值的进口标准样品SQC014中TOC的标准值为(5.21±2.28)mg/g,实测值为4.45mg/g。采用重铬酸钾氧化-还原容量法定值的中国标准样品GBW07430中TOC的标准值为(9.7±1.2)mg/g,实测值为10.9mg/g,质控样品实测值均在不确定度范围内。2.4.3 标准值的认定和不确定度评估
经过对全部测定结果的统计处理,最终以多家实验室协作测定统计结果的总均值作为土壤总有机碳标准样品的标准值,该标准值为25.2mg/g。
本文土壤总有机碳标准样品的不确定度(uCRM)由测定不确定度分量(uchar)、瓶间均匀性不确定度分量(ubb)和长期稳定性不确定度分量(ults)组成。
测定不确定度分量(uchar)计算方法如下:
$ {u_{{\rm{char }}}} = s/\sqrt p $
(4) 式(4)中:s为总均值标准偏差;p为检测数据的组数。
经计算,测定不确定度(uchar)结果为0.52mg/g,相对测定不确定度为2.1%。
标准样品瓶间均匀性不确定度(ubb)是表征土壤样品均匀性的定量指标。而检测方法的重复性水平是评定瓶间均匀性不确定度的重要影响因素。当样品检测方法的重复性很好时,能够反映出瓶间很小的不均匀性,出现瓶间均方(MSb)明显大于瓶内均方(MSw)的情况,可用瓶间均匀性标准偏差来评估样品的瓶间均匀性不确定度,计算公式为:
$ {u_{{\rm{bb}}}} = \sqrt {\frac{{{\rm{M}}{{\rm{S}}_{\rm{b}}} - {\rm{M}}{{\rm{S}}_{\rm{w}}}}}{n}} $
(5) 当样品检测方法的重复性欠佳时,出现瓶间均方(MSb)接近或小于瓶内均方(MSw)的情况,采用检测方法的重复性标准偏差对瓶间均匀性标准偏差的贡献来评估样品的瓶间均匀性不确定度,计算公式为:
$ {u_{{\rm{bb}}}} = \sqrt {\frac{{{\rm{M}}{{\rm{S}}_{\rm{w}}}}}{n} \times \sqrt {\frac{2}{{{v_{\rm{w}}}}}} } $
(6) 根据方差检验结果,以公式(6)计算本标准样品的瓶间均匀性不确定度(ubb),其结果为0.36mg/g,相对瓶间均匀性不确定度为1.5%。
根据稳定性统计结果,计算标准样品的长期稳定性不确定度(ults)。对于没有显著稳定性变化趋势的标准样品,其长期稳定性不确定度计算公式为:
$ {u_{{\rm{lts}}}} = s\left( {{b_1}} \right) \times t $
(7) 式(7)中:s(b1)为斜率的标准偏差;t为稳定性检验的总时间(月)。
土壤总有机碳标准样品的稳定性检测共进行了6次,总时间累计18个月,即n=6,且t=18。经计算,长期稳定性不确定度(ults)为0.29mg/g,相对长期稳定性不确定度为1.2%。
土壤总有机碳标准样品的不确定度(uCRM)的计算公式为:
$ {u_{{\rm{CRM}}}} = \sqrt {u_{{\rm{char}}}^2 + u_{{\rm{bb}}}^2 + u_{{\rm{Its}}}^2} $
(8) 本文土壤总有机碳标准样品的不确定度(uCRM)为0.7mg/g,相对不确定度为2.8%。不确定度采用扩展不确定度(UCRM)表示,取包含因子k=2,并按照“只进不舍”原则对结果进行保留[18, 35]。则样品的扩展不确定度(UCRM)为1.4mg/g,相对扩展不确定度为5.6%。即土壤总有机碳标准样品的标准值和扩展不确定度为(25.2±1.4)mg/g。
2.4.4 标准值的溯源性
为保证样品标准值的质量,尽可能减小误差,本标准样品的研制主要采取了以下6项具体措施保证溯源性。①制作标准曲线的试剂均由基准试剂配制,可溯源至SI测量国际单位。②用于样品分析检测的方法均为国家或行业标准分析方法,且经实践证明是成熟、准确、可靠的方法,方法的准确度、精密度及检出限均满足标准样品分析要求。③实验所用计量器具、天平、测量仪器等均通过国家计量部门检定或校准,量值准确可靠,并可溯源至国家标准。④分析过程中,确保分析试剂的纯度,每次分析时均进行空白试验,扣减空白并进行背景校正。⑤参加本项目协作定值的实验室通过国家级计量认证认可,具备参与标准样品定值的能力、经验及质量保证体系。⑥各实验室均采用国家有证标准样品对测定全过程进行质量控制,确保分析结果准确可靠。
3. 结论
采集新疆农用地环境土壤作为原材料,经除杂、干燥、研磨混匀、分装灭菌等加工步骤,制备成土壤中总有机碳标准样品。采用燃烧氧化-非分散红外法分析测试并对数据结果进行统计检验,在最小取样量为0.1g时,该标准样品中总有机碳均匀性良好,在室温避光条件下的稳定性良好。组织多家具备能力和资质的实验室以燃烧氧化-非分散红外法和重铬酸钾容量法两种方法进行协作测定,在对检测结果进行量值评定后,确定了土壤中总有机碳标准样品的标准值和不确定度为(25.2±1.4)mg/g。该标准样品的研制程序完全符合GB/T 15000系列和HJ 173—2017等标准规范的要求。
研制的土壤中总有机碳标准样品为采用燃烧氧化-非分散红外法进行分析的标准样品。该样品中有机碳含量水平和不确定度范围等技术指标,能够满足中国土壤生态环境监测工作质量的要求,进一步完善了中国土壤环境标准样品体系,可广泛应用于土壤环境监测和土壤有机碳相关研究的质量控制、土壤有机碳分析仪的校准、土壤有机碳检测方法的评价,以及相关检测能力的考核等方面的工作。
-
表 1 样品全谱拟合分析结果和晶体结构参数
Table 1 Phase analysis results and crystallographic parameters obtained by Rietveld refinement of X-ray diffraction
样品编号 判定因子 卡片号 矿物名称 矿物含量(%) 比重(g/cm3) 晶胞体积(nm3) 晶体结构参数 a(nm) b(nm) c(nm) α(°) β(°) γ(°) 样品1 Rp=6.86%
Rwp=9.19%98-004-0799 石英 69.6 2.65 0.1130 0.4913 0.4913 0.5404 90.00 90.00 120.00 98-000-5699 蓝晶石 22.6 3.67 0.2932 0.5573 0.7119 0.7847 73.99 89.95 78.87 98-001-7444 叶腊石 3.9 2.73 0.2175 0.5286 0.5194 0.9341 79.49 85.41 59.61 98-002-8230 高岭石 0.6 2.56 0.1651 0.5143 0.5172 0.7403 84.14 75.20 60.15 98-004-5732 长石 0.5 2.59 0.3572 0.7181 0.7595 0.7899 113.23 103.90 103.54 98-006-1852 白云母 1.5 2.87 0.9307 0.5170 0.8946 2.0210 90.00 95.36 90.00 98-007-6735 金红石 1.2 4.25 0.0624 0.4592 0.4592 0.2961 90.00 90.00 90.00 样品2 Rp=5.67%
Rwp=8.03%98-004-5584 蓝晶石 60.6 3.67 0.2932 0.5573 0.7119 0.7847 73.99 89.96 78.89 98-005-3283 刚玉 29.6 3.99 0.2548 0.4759 0.4759 1.2990 90.00 90.00 120.00 98-002-3352 硬水铝石 5.6 3.32 0.1180 0.9423 0.2845 0.4402 90.00 90.00 90.00 98-001-1963 云母 0.5 2.83 0.7015 0.5171 0.5171 3.0291 90.00 90.00 120.00 98-001-7624 绿脆云母 1 3.25 0.4574 0.5232 0.9036 0.9842 90.00 100.55 90.00 98-001-7216 地开石 0.9 2.57 0.6564 0.5150 0.8899 1.4393 90.00 95.66 90.00 98-005-0903 重晶石 1.7 4.49 0.3451 0.8875 0.5475 0.7102 90.00 90.00 90.00 样品3 Rp=7.17%
Rwp=9.72%98-005-6873 石英 58.1 2.65 0.1129 0.4913 0.4913 0.5404 90.00 90.00 120.00 98-006-1891 蓝晶石 27.7 3.67 0.2931 0.5574 0.7115 0.7845 74.09 90.02 78.87 98-003-8360 铁铝榴石 1.8 4.31 1.5337 1.1532 1.1532 1.1532 90.00 90.00 90.00 98-007-7906 黑云母 3.4 3.06 0.4993 0.5361 0.9249 1.0246 90.00 100.64 90.00 98-001-1963 白云母 0.5 2.84 0.6998 0.5199 0.5199 2.9892 90.00 90.00 120.00 98-002-8573 透辉石 2.5 3.33 0.4366 0.9768 0.8908 0.5220 90.00 106.006 90.00 98-005-6870 长石 4.3 2.60 0.3353 0.7038 0.7648 0.7916 120.09 104.98 98.74 98-003-7562 磷灰石 1.7 3.40 0.5243 0.9380 0.9380 0.6882 90.00 90.00 120.00 样品4 Rp=7.61%
Rwp=10.45%98-004-5584 蓝晶石 95.0 3.67 0.2932 0.5574 0.7118 0.7846 73.99 89.96 78.89 98-004-0799 石英 2.0 2.64 0.1135 0.4912 0.4912 0.5431 90.00 90.00 120.00 98-004-5786 白云母 3.0 2.91 0.9054 0.5196 0.8610 2.0307 90.00 94.74 90.00 表 2 方法回收率
Table 2 Recovery tests of the method
样品编号 加入量(g) 回收量(g) 回收率(%) 样品1 0.9502 0.9505 100.0 样品3 0.9502 0.9355 98.4 样品2 0.9502 0.9058 95.3 样品4 0.9502 0.9608 101.0 表 3 Highscore和Jade衍射数据处理软件分析结果对比
Table 3 Comparison of analytical results determined by Highscore and Jade diffraction data treatment software
样品编号 Highscore软件 Jade软件 双差(%) 相对偏差(%) 成分 含量(%) 成分 含量(%) 样品1 低温石英 69.6 石英 69.7 0.1 0.19 蓝晶石 22.6 蓝晶石 23.1 0.5 2.32 叶腊石 3.9 叶腊石 2.4 1.5 46.4 高岭石 0.6 高岭石 1.0 0.4 42.9 长石 0.5 长石 0.8 0.3 53.1 云母 1.5 白云母 1.7 0.2 7.50 金红石 1.2 金红石 1.2 0 0.80 样品2 蓝晶石 60.6 蓝晶石 59.8 0.8 1.33 刚玉 29.6 刚玉 28.5 1.1 3.68 硬水铝石 5.6 硬水铝石 5.4 0.2 4.37 云母 0.5 1M型云母 1.6 - - 绿脆云母 1.0 2M型云母 2.5 - - 地开石 0.9 高岭石 1.2 - - 重晶石 1.7 重晶石 0.9 0.8 59.2 样品3 石英 58.1 石英 57.4 0.77 1.33 蓝晶石 27.7 蓝晶石 28.0 0.3 1.08 铁铝榴石 1.8 铁铝榴石 0.5 1.3 117 黑云母 3.4 黑云母 0.8 2.6 124 白云母 0.5 白云母 4.9 4.4 164 透辉石 2.5 透辉石 1.3 1.2 66.0 长石 4.3 长石 5.7 1.4 29.6 磷灰石 1.7 羟基磷灰石 1.4 0.3 17.5 样品4 蓝晶石 95.0 蓝晶石 95.4 0.4 0.37 石英 2.0 石英 1.9 0.1 2.05 白云母 3.0 白云母 2.7 0.3 10.5 表 4 样品化学成分分析误差统计
Table 4 Statistics of analytical error for the chemical components in samples
样品编号 参数 SiO2 Al2O3 Fe2O3 CaO MgO K2O Na2O TiO2 P2O5 MnO 样品1 含量(%) 78.76 14.69 2.41 0.29 0.05 0.26 0.06 1.24 0.25 0.01 偏差(%) 3.13 1.51 -2.36 - -0.05 -0.01 -0.06 0.01 - - 相对偏差(%) 3.98 10.25 -97.85 - -91.45 -4.41 -88.73 0.81 - - 样品2 含量(%) 23.29 71.53 0.2 0.32 - 0.16 0.16 0.064 0.22 - 偏差(%) -0.01 2.00 -0.19 0.05 - -0.10 - - - - 相对偏差(%) -0.03 2.79 -95.25 14.43 - -61.09 - - - - 样品3 含量(%) 74.56 18.46 2.42 1.19 0.32 1 0.2 0.22 0.82 0.038 偏差(%) 0.35 0.61 -0.68 0.24 0.57 -0.59 0.33 -0.14 -0.17 -0.01 相对偏差(%) 0.47 3.31 -28.14 19.76 178.84 -59.32 164.00 -63.65 -20.51 -37.70 样品4 含量(%) 37.9 59.73 0.48 0.03 < 0.24 < 0.04 0.03 0 偏差(%) 0.71 1.17 - - - 0.11 - - - - 相对偏差(%) 1.86 1.96 - - - 47.77 - - - - -
王梅英.蓝晶石的化学物相分析[J].非金属矿, 2006, 29(1):15-16. http://d.old.wanfangdata.com.cn/Periodical/fjsk200601007 Wang M Y.Chemical phase analysis of cyanite[J].Non-Metallic Mines, 2006, 29(1):15-16. http://d.old.wanfangdata.com.cn/Periodical/fjsk200601007
王梅英, 李鹏程, 李艳华, 等.蓝晶石矿中氟钠镁铝硅铁钛钾钙元素的X射线荧光光谱分析[J].岩矿测试, 2013, 32(6):909-914. doi: 10.3969/j.issn.0254-5357.2013.06.011 Wang M Y, Li P C, Li Y H, et al.Analysis of F, Na, Mg, Al, Si, Fe, Ti, K and Ca in cyanite ores by X-ray fluorescence spectrometry[J].Rock and Mineral Analysis, 2013, 32(6):909-914. doi: 10.3969/j.issn.0254-5357.2013.06.011
杜晓冉, 智红梅, 张金矿, 等.X射线衍射法定量分析蓝晶石样品[J].理化检验(化学分册), 2013, 49(4):402-404. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QKC20132013071700099060 Du X R, Zhi H M, Zhang J K, et al.Quantitative analysis of kyanite sample with X-ray diffractometry[J].Physical Testing and Chemical Analysis (Part B:Chemical Analysis), 2013, 49(4):402-404. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QKC20132013071700099060
邱贤荣, 齐砚勇, 唐志强.全谱拟合定量分析石灰石[J].分析科学学报, 2013, 29(1):146-148. http://d.old.wanfangdata.com.cn/Periodical/fxkxxb201301035 Qiu X R, Qi Y Y, Tang Z Q.Rietveld quantitative analysis of limestone[J].Journal of Analytical Science, 2013, 29(1):146-148. http://d.old.wanfangdata.com.cn/Periodical/fxkxxb201301035
冉敬, 杜谷, 王凤玉.X射线衍射全谱拟合法快速分析长石矿物含量[J].岩矿测试, 2017, 36(5):489-494. doi: 10.15898/j.cnki.11-2131/td.201610110154 Ran J, Du G, Wang F Y.Rapid analysis of feldspar by X-ray diffractometry Rietveld refinement method[J].Rock and Mineral Analysis, 2017, 36(5):489-494. doi: 10.15898/j.cnki.11-2131/td.201610110154
迟广成, 肖刚, 汪寅夫, 等.铁矿石矿物组分的X射线粉晶衍射半定量分析[J].冶金分析, 2015, 35(1):38-44. http://d.old.wanfangdata.com.cn/Periodical/yjfx201501007 Chi G C, Xiao G, Wang Y F, et al.Semi-quantitative analysis of the mineral components of iron ores by X-ray powder diffraction[J].Metallurgical Analysis, 2015, 35(1):38-44. http://d.old.wanfangdata.com.cn/Periodical/yjfx201501007
洪汉烈, 陈建军, 杨淑珍, 等.水泥熟料定定量分析的全谱拟合法[J].分析测试学报, 2001, 20(2):5-8. doi: 10.3969/j.issn.1004-4957.2001.02.002 Hong H L, Chen J J, Yang S Z, et al.Quantitative phase analysis of cement clinker by Rietveld full pattern fitting method[J].Journal of Instrumental Analysis, 2001, 20(2):5-8. doi: 10.3969/j.issn.1004-4957.2001.02.002
Gualtieri M L, Romagnoli M, Miselli P, et al.Full quantitative phase analysis of hydrated lime using the Rietveld method[J].Cement & Concrete Research, 2012, 42(9):1273-1279. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=70d31ed7a08e6334b2700a3083fec6dc
Santini T C.Application of the Rietveld refinement method for quantification of mineral concentrations in bauxite residues (alumina refining tailings)[J].International Journal of Mineral Processing, 2015, 139:1-10. doi: 10.1016/j.minpro.2015.04.004
Woodruff L, Cannon W F, Smith D B, et al.The distribu-tion of selected elements and minerals in soil of the conterminous United States[J].Journal of Geochemical Exploration, 2015, 154:49-60. doi: 10.1016/j.gexplo.2015.01.006
Gentili S, Comodi P, Bonadiman C, et al.Mass balance vs Rietveld refinement to determine the modal composition of ultramafic rocks:The case study of mantle peridotites from Northern Victoria Land (Antarctica)[J].Tectonophysics, 2015, 650:144-155. doi: 10.1016/j.tecto.2015.01.024
肖松.河南隐山蓝晶石矿物成因的研究[J].中国矿业, 2005, 14(2):47-49. doi: 10.3969/j.issn.1004-4051.2005.02.015 Xiao S.On mineral genesis of Yinsan kyanite in Henan[J].China Mineral Magazine, 2005, 14(2):47-49. doi: 10.3969/j.issn.1004-4051.2005.02.015
李锁成, 张松林, 火军昌, 等.萨尔哈布塔勒蓝晶石矿床地质及矿物学特征[J].矿产保护与利用, 2015(1):12-15. http://d.old.wanfangdata.com.cn/Periodical/kcbhyly201501003 Li S C, Zhang S L, Huo J C, et al.The process mineralogy and geological characters of Saerhabutale Kyanite deposit in Gansu Province[J].Conservation and Utilization of Mineral Resources, 2015(1):12-15. http://d.old.wanfangdata.com.cn/Periodical/kcbhyly201501003
万红波, 廖立兵.膨润土中蒙脱石物相的定量分析[J].硅酸盐学报, 2009, 37(12):2055-2060. doi: 10.3321/j.issn:0454-5648.2009.12.017 Wan H B, Liao L B.Quantitative phase analysis of montmorillonite in bentonite[J].Journal of the Chinese Ceramic Society, 2009, 37(12):2055-2060. doi: 10.3321/j.issn:0454-5648.2009.12.017
房俊卓, 张霞, 徐崇福.实验条件对X射线衍射物相定量分析结果的影响[J].岩矿测试, 2008, 27(1):60-62. doi: 10.3969/j.issn.0254-5357.2008.01.015 Fang J Z, Zhang X, Xu C F.Effect of experimental conditions on X-ray diffractometric quantitative phase analysis[J].Rock and Mineral Analysis, 2008, 27(1):60-62. doi: 10.3969/j.issn.0254-5357.2008.01.015
马礼敦.X射线粉末衍射的新起点——Rietveld全谱拟合[J].物理学进展, 1996, 16(2):251-271. doi: 10.3321/j.issn:1000-0542.1996.02.005 Ma L D.The new starting point of X-ray powder diffraction-Rietveld whole pattern fitting[J].Progress in Physics, 1996, 16(2):251-271. doi: 10.3321/j.issn:1000-0542.1996.02.005
林伟伟, 宋友桂.沉积物中X射线衍射物相定量分析中的两种方法对比研究[J].地球环境学报, 2017, 8(1):83-86. http://d.old.wanfangdata.com.cn/Periodical/dqhjxb201701010 Lin W W, Song Y G.A comparative study on X-ray diffraction mineral quantitative analysis of two methods in sediments[J].Journal of Earth Enviroment, 2017, 8(1):83-86. http://d.old.wanfangdata.com.cn/Periodical/dqhjxb201701010
伍月, 刘欣, 张波, 等.X射线粉晶衍射基体清洗法在矿物定量分析中的应用[J].地质与资源, 2017, 26(3):323-328. http://d.old.wanfangdata.com.cn/Periodical/gjsdz201703017 Wu Y, Liu X, Zhang B, et al.The application and research of X-ray powder diffraction matrix flushing method in quantitative analysis[J].Geology and Resources, 2017, 26(3):323-328. http://d.old.wanfangdata.com.cn/Periodical/gjsdz201703017
施倪承, 白文吉, 马喆生, 等.西藏罗布莎铬铁矿床中的金刚石包体X射线衍射研究[J].地质学报, 2002, 76(4):496-500. doi: 10.3321/j.issn:0001-5717.2002.04.008 Shi N C, Bai W J, Ma Z S, et al.A study of X-ray diffraction of diamond inclusions from Luobusha, Tibet[J].Acta Geologica Sinica, 2002, 76(4):496-500. doi: 10.3321/j.issn:0001-5717.2002.04.008
王德强, 王辅亚, 王冠鑫, 等.钾含量对白云母X射线衍射特征的影响[J].大地构造与成矿学, 2004, 28(1):69-73. doi: 10.3969/j.issn.1001-1552.2004.01.010 Wang D Q, Wang F Y, Wang G X, et al.XRD character of muscovites with different potassium contents[J].Geotectonica et Metallogenia, 2004, 28(1):69-73. doi: 10.3969/j.issn.1001-1552.2004.01.010
杜谷, 王坤阳, 冉敬, 等.红外光谱/扫描电镜等现代大型仪器岩石矿物鉴定技术及其应用[J].岩矿测试, 2014, 33(5):625-633. doi: 10.3969/j.issn.0254-5357.2014.05.003 Du G, Wang K Y, Ran J, et al.Application of IR/SEM and other modern instruments for mineral identification[J].Rock and Mineral Analysis, 2014, 33(5):625-633. doi: 10.3969/j.issn.0254-5357.2014.05.003
-
期刊类型引用(10)
1. 乔柱,王恒,苏明跃,赵秀荣,张萍萍,王永红. 质量吸收系数修正-X射线衍射外标法定量分析氧化铁皮中4种物相. 冶金分析. 2024(07): 55-61 . 百度学术
2. 李磊,张青,范慧,阳珊,遇祯. 土壤中高岭石的X射线衍射全谱拟合定量分析方法研究. 安徽地质. 2023(03): 282-287 . 百度学术
3. 张启燕,史维鑫,刘晓,回广骥,原春雨. 高光谱扫描在碳酸盐岩矿物组成分析中的应用. 岩矿测试. 2022(05): 815-825 . 本站查看
4. 赵彦钊,呼浩,杨崔月,胡智敏,李晋东,殷海荣. 锂电池正极材料烧成用匣钵物相的半定量分析. 陕西科技大学学报. 2021(05): 126-130+145 . 百度学术
5. 徐速,戴苏兰,石小平,陈大鹏,曲蔚,罗琴凤,夏玉梅. 碧石的宝石学分类及定义界定初探. 矿物岩石. 2021(04): 1-13 . 百度学术
6. 温利刚,贾木欣,王清,付强,赵建军. 自动矿物学:Ⅰ.技术进展与应用. 中国矿业. 2020(S1): 341-349 . 百度学术
7. 徐颖,邓利蓉,芦玉峰,左联,杜广报. 热处理对柯尔碱膨润土微观结构和物化性能的影响. 岩矿测试. 2019(03): 280-287 . 本站查看
8. 郑建国,于南树,刘永秀,胡善明,胡兆平. 利用磷尾矿合成新型缓释肥料聚磷酸钙镁的试验研究. 化工矿物与加工. 2019(07): 64-67+71 . 百度学术
9. 冉敬,郭创锋,杜谷,王凤玉. X射线衍射全谱拟合法分析蓝晶石的矿物含量. 岩矿测试. 2019(06): 660-667 . 本站查看
10. 陈乐,赵昕,游翠纯,朱肖峰. 金红石晶种的X射线衍射分析及其转化率测定. 分析测试技术与仪器. 2018(04): 236-240 . 百度学术
其他类型引用(4)