Determination of Aluminum Oxide in Rare Earth Ore by Inductively Coupled Plasma-Atomic Emission Spectrometry
-
摘要: 铝是稀土矿石常检的杂质元素，目前采用电感耦合等离子体发射光谱法(ICP-AES)基体匹配校正模式测定矿石中的铝量，对基体中铝的空白及基体量有一定的要求；且稀土矿石的品种很多，铝在其中的存在形式也较复杂，简单的酸溶很难彻底地将铝转入溶液中进行准确测定。本文采用碱熔法处理样品，滤液酸化后用ICP-AES测定稀土矿石中的Al2O3。通过共存元素干扰实验发现样品中的稀土元素及钍对铝的测定产生严重的光谱干扰，提出预先以碱分离除去消除干扰，有效地降低了检测下限；采用基体校正模式，消除了盐分对测定的影响。方法检出限为0.021~0.035 mg/g，测定范围为0.50%~8.00%，精密度(RSD) < 7.1%。对不同含量的样品进行分析，测定值与化学容量法测定结果基本相符。该方法具有测定含量范围宽、分析速度快、结果准确等优点。
-
关键词:
- 稀土矿石 /
- 三氧化二铝 /
- 碱熔 /
- 电感耦合等离子体发射光谱法
Abstract: The element Al is determined frequently in rare earth ore since Al is a common impurity element in this ore. At present, the Inductively Coupled Plasma-Atomic Emission Spectrometric (ICP-AES) method, which has certain requirements of the Al blank and the matrix concentration, is used to determine the content of Al in ore with a matrix matching calibration. Because Al has complicated patterns in many kinds of rare earth ores, it is difficult to dissolve Al completely into solution for determination by a routine acid-leaching method. In this paper, a new method is described, whereby the sample is digested by alkali fusion and the filtrate acidulated before the content Al in rare earth ore is determined by ICP-AES. The effect of spectrum interferences of REEs and Th on the Al was studied. The detection limit was reduced effectively by eliminating the interferences of coexisting elements. The matrix correction mode was undertaken to eliminate the salt influence. The determination limit of the method was from 0.021 mg/g to 0.035 mg/g and the determination content range was from 0.50 % to 8.00%. The RSD was less than 7.1%. The result of the samples with the presented method is consistent with the volumetric method. The method is fast with a wide measuring range and extremely high precision. -
海相碳酸盐碳同位素(δ13C)组成可以有效地恢复地质历史时期海水的碳同位素组成，在古海洋碳循环[1, 2, 3]、古气候[4, 5]和地层对比[6, 7, 8, 9, 10]中发挥着极为重要的作用。由于碳酸盐氧同位素(δ18O)容易受到后期水岩交换和成岩作用的影响，因此不能够准确代表古海水的原始信息[11]，但其变化趋势在一定程度上可以指示古气候的变化情况。例如，一些碳酸盐壳体的生物化石可以有效地恢复古水温的变化情况[12, 13]。因此，建立可靠的碳酸盐碳、氧同位素分析方法对上述研究领域非常重要。
碳酸盐碳氧同位素分析中的制样过程是将碳酸根以二氧化碳的形式提取出来，利用气体稳定同位素质谱仪对碳氧同位素组成进行测定。不同的制样方法对样品性质、样品量的要求不同，分析结果的精度也各有差异。因此，制样方法在碳酸盐碳氧同位素的分析测试中十分重要。目前碳酸盐碳、氧同位素的制样方法包括：元素分析仪燃烧法[14, 15]、激光熔蚀法[16, 17, 18]和磷酸分解法。由于前两种方法存在一定的局限性，当前测定碳酸盐δ13C、δ18O的制样方法以磷酸分解法为主。磷酸分解制样法可分为：传统离线法制备CO2，经双路进样系统进入气体稳定同位素质谱仪分析测定[19]和在线制备CO2与连续流气体稳定同位素质谱仪联用[20](如GasBench-MAT253、MultiPrep-IsoPrimeⅡ)。传统的磷酸分解法制样耗时、样品和磷酸消耗量大[20]，易引入人为误差[21]。因此，不能满足当前研究人员对微区、微量碳氧同位素的分析测试需求。GasBenchⅡ技术需样量小(100~400 μg)，可以自动完成排空、加酸、气体采集、分离等步骤，之后用质谱仪进行碳氧同位素测定。利用微钻对样品进行原位采样，可以实现碳酸盐岩微区、微量碳氧同位素的分析测试。目前，GasBenchⅡ-MAT253连续流质谱仪联用(GasBenchⅡ-IRMS)测定碳酸盐岩δ13C、δ18O同位素组成已被广泛应用。
已有研究认为，在GasBenchⅡ-IRMS分析中，不同的实验条件直接影响着分析δ13C、δ18O同位素测试结果的准确性。例如，杜广鹏等[22]和Paul等[23]研究了顶空进样瓶的排空时间对δ13C、δ18O分析精度和准确性的影响，发现当排空时间不足时，空气会对样品分析结果造成干扰。陶成等[21]分析了色谱分离温度对测试结果的影响，发现色谱柱柱温在低温范围(40℃，65℃)对分析结果影响不大，当柱温为120℃时，空气峰与CO2峰发生重叠，影响分析结果。Breitenbach等[24]和Paul等[25]分别考察了样品质量、反应温度和反应时间对δ13C、δ18O同位素测试结果的影响，研究表明，样品质量过小会造成分析精度下降；反应温度过低会导致酸泵进样针堵塞，同时也直接影响平衡时间的长短；反应时间过短会导致同位素不能达到平衡状态而影响分析结果。但是，这些研究都是针对单一或者某两个实验条件对同位素分析结果的影响，并未对该联用方法所涉及的实验条件进行逐一考察和优化，因此不能为实验条件的选择提供全面的参考信息。本文在前人的工作基础上，利用GasBenchⅡ-MAT253联用方法，选取两个国家标准物质和一个实验室石笋标样作为研究对象，综合考察顶空瓶排空时间、反应温度、反应时间和色谱分离温度，以及同位素结果标准归一化方法对碳酸盐岩δ13C、δ18O分析结果的影响，拟为今后优化和选取合适的分析条件、获得可靠的碳酸盐岩δ13C和δ18O结果提供一定的科学参考。
1. 实验方法
1.1 仪器和主要试剂
GasBenchⅡ前处理装置(CTC Analytics，瑞士)：GC Combi PAL自动进样器，恒温样品盘，Poraplot Q色谱柱(30 m×0.32 mm)，酸泵，硼硅酸盐反应瓶(12 mL)，进样针和吹气针。
Thermo FisherTM MAT253气体稳定同位素质谱仪：高灵敏电子轰击离子源，ISODAT NT 2.0控制软件。
高纯He气(纯度优于99.999%)，CO2参考气(纯度优于99.995%)，正磷酸99%分析级晶体EMSURE(Merck KGaA，德国)。
1.2 分析样品
本文选取不同碳酸盐标准物质作为研究对象进行分析测试，包括GBW 04416 (TB-1)、GBW 04417 (TB-2) 为国家标准物质，石笋标样为实验室内部标准样品，均具有不同的碳氧同位素值。
1.3 分析方法
称取适量样品于反应瓶中，用带硅胶隔垫的盖子密封，通过排空针按顺序依次排空反应瓶中的空气，消除其对δ13C、δ18O分析结果的影响。经排空处理后，通过酸泵和酸针向反应瓶中加入7滴磷酸(足量)，在设定温度下反应一定时间，磷酸与碳酸盐反应释放出的CO2在高纯氦气的带动下在线除水，通过色谱柱与其他杂质气体分离；分离后的CO2直接进入MAT 253稳定同位素质谱仪进行δ13C、δ18O测定。分析方法详见文献[20, 21, 22]。
2. 实验条件的考察与优化
2.1 排空时间
采用GasBenchⅡ在线制样时，利用顶空排气针向反应瓶中注入高纯氦气，排空反应瓶中的空气，从而消除空气中CO2等气体对测试结果的影响。排空时间的长短是影响分析测试结果的重要因素：时间过短，反应瓶中残余的空气会对测试结果产生干扰[22-23]；时间过长，一定程度上浪费机时和载气。
按1.2节的分析流程，挑取适量(约300 μg；保证样品m/z=44的信号强度大于2000 mV，否则测量内精度变差[24])国家标准物质和实验室标样(按TB-1—TB-2—石笋标样—TB-1—TB-2交叉顺序放入恒温样品盘中，共8组)，以2 min为时间间隔设定排空时间依次为1、3、5、7、9、11、12、13 min，反应温度为72℃，反应时间为60 min，色谱分离温度为60℃。
实验结果表明(表 1)：当排空时间为1~7 min时，在样品峰之前出现Paul等[23]报道的空气次级峰，次级峰的信号强度随排空时间的增加逐渐减小。Paul等[23]指明有空气次级峰存在的情况下，测试结果均不可靠。因此，该结果表明仪器使用手册推荐的5 min排空时间[26]不能完全消除空气干扰。当排空时间大于9 min时，无明显空气次级峰，同位素值较稳定。当排空时间增加至12 min和13 min时，测试精度进一步提高。综合考虑分析效率和氦气的消耗量，选取12 min作为排空时间。
表 1 不同排空时间下石笋标准的碳氧同位素组成Table 1. δ13C>and δ18Ocompositions of the stalagmite standard on varied flushing times排空时间
(min)δ13CVPDB
(‰)SD
(1σ, n=4)δ18OVPDB
(‰)SD
(1σ, n=4)1~7 存在空气次级峰 存在空气次级峰 9 -10.740 0.038 -8.567 0.010 11 -10.784 0.026 -8.480 0.010 12 -10.774 0.012 -8.563 0.004 13 -10.798 0.003 -8.574 0.010 2.2 反应温度
反应温度是磷酸法制备CO2气体的一个重要条件，磷酸与碳酸盐反应速率随温度的升高而加快[25]，由于同位素分馏系数与温度密切相关[27]，反应过程中的温度直接决定着分析结果的准确性[28]。按2.1节实验方法，设定排空时间为12 min，反应时间为60 min，色谱分离温度为60℃。Paul等[25]通过研究指出碳酸盐与磷酸反应温度应控制在50~90℃之内，由于温度较低磷酸易凝固堵塞酸针，所以选择的最低温度为50℃，铝制样品盘可调的最高温度为80℃。综上，本文考察反应温度分别为50℃、60℃、72℃和80℃对同位素测试结果的影响。
结果表明(表 2)：反应温度为50℃、60℃和72℃三组实验中δ13C测试精度略优于反应温度为80℃的结果，该现象与Paul等[25]报道的结果一致。由此表明50~72℃范围的温度均为合适的反应温度。此外，实验过程中还发现80℃可以造成反应瓶的密封性变差。在保证测试精度相当的情况下，选择较高的反应温度可以缩短反应时间，因此选择72℃作为反应温度。
表 2 不同反应温度下石笋标准的碳氧同位素组成Table 2. δ13C and δ18O compositions of the stalagmite standard on varied temperatures反应温度
(℃)δ13CVPDB
(‰)SD
(1σ, n=4)δ18OVPDB
(‰)SD
(1σ, n=4)50 -10.785 0.031 -8.588 0.024 60 -10.726 0.050 -8.629 0.013 72 -10.758 0.039 -8.550 0.027 80 -10.811 0.065 -8.577 0.022 2.3 反应时间
磷酸与碳酸盐在一定温度下需反应一定时间才能反应完全，气相-液相实现同位素平衡也需要一定时间[25]，同位素达到平衡后才能获取准确的测试结果。按2.1节实验方法，设定排空时间为12 min，反应温度为72℃，色谱分离温度为60℃。在常温下(25℃)，碳酸盐和磷酸需反应24 h才能达到同位素平衡。当反应温度升高至50℃以上时，反应时间可以缩短至1 h。参考文献[25]，本文考察的反应时间范围为30~90 min，由于测量一个样品需15 min，故设定反应时间的增幅为15 min。
结果表明(表 3)，反应时间在30~75 min范围内变化，同位素结果差异较小，分析精度较稳定，这与Paul等[25]的实验结论相吻合；反应时间为90 min时，δ13C偏离其他四组实验结果约-0.2‰，而氧同位素无明显差异，Tu等[29]也报道了类似的实验结果。不同碳酸盐矿物与磷酸的反应速率不同，郑永飞等[27]提出了方解石、石灰石和菱铁矿等的反应温度和所需时间，Paul等[25]通过实验发现碳酸盐与磷酸反应不彻底会导致结果发生偏离。综合本文选取石笋标样的物质组成以及气-液两相同位素的平衡时间，选取仪器厂商推荐的60 min作为反应时间。
表 3 不同反应时间下石笋标准的碳氧同位素组成Table 3. δ13C and δ18O compositions of the stalagmite standard on varied reaction times反应时间
(min)δ13CVPDB
(‰)SD
(1σ, n=4)δ18OVPDB
(‰)SD
(1σ, n=4)30 -10.776 0.009 -8.562 0.036 45 -10.749 0.046 -8.512 0.006 60 -10.765 0.036 -8.524 0.009 75 -10.760 0.011 -8.533 0.044 90 -10.944 0.036 -8.522 0.027 2.4 色谱柱分离温度
色谱柱温度直接影响分离效能[30]，氦气黏度因色谱分离温度的变化而改变，进而影响流动相的分离效果，流动相的分离差异是导致同位素分馏的重要因素。按2.1中实验方法，设定排空时间为12 min，反应温度为72℃，反应时间为60 min，考察色谱分离温度为45℃[22]、60℃[31]和70℃[20, 32]对测试结果的影响。
结果表明(表 4)，考察的色谱分离温度均能保证空气干扰峰(N2O)与CO2峰完全分离。色谱分离温度为45℃和60℃时，测试结果基本一致，45℃条件下分析精度相对较差。并且，色谱分离温度为45℃时，温控显示有±1℃的波动，可能是造成精度变差的原因。分离温度为70℃时，δ13C分析结果发生一定偏离。因此，本文选取的色谱分离温度为60℃。
表 4 不同色谱分离温度下石笋标准的碳氧同位素组成Table 4. δ13Cand δ18Ocompositions of the stalagmite standard on varied chromatographic separation temperatures色谱温度
(℃)δ13CVPDB
(‰)SD
(1σ, n=4)δ18OVPDB
(‰)SD
(1σ, n=4)45 -10.823 0.050 -8.488 0.053 60 -10.811 0.024 -8.485 0.036 70 -10.760 0.024 -8.584 0.047 3. 同位素数据的处理方法
通过以上实验条件的优化与选择，为获得碳酸盐碳氧同位素可靠的分析结果提供了重要前提。同位素数据归一化处理方法也会对最终的同位素分析结果产生重要的影响[25]。因此，在确定最优实验条件的前提下，选择合适的同位素数据归一化方法是对分析结果的进一步保证。同位素数据的归一化(标准化)指将仪器计算的同位素原始数据(分析物相对参考气的同位素值)转化为相对国际标准的同位素标准值[33, 34, 35]，以便进行同位素结果对比。Paul等[34]认为样品分析结果的准确性在一定程度上取决于采用的归一化方法。Finnigan GasBench操作说明书[26]推荐使用单一标准物质校准，但我们在工作中发现，当实际样品同位素组成偏离标准物质较大时，会造成较大的校准误差[33, 34]。本文在获得近4000件实际样品数据的基础上，对比了单一标准物质和双标准物质线性校准(即标准物质的推荐值与测量值的线性方程)的计算结果，发现双标准物质校准偏差小于单一标准物质校准偏差[33, 34]。
4. 结语
本文系统地考察了GasBenchⅡ-IRMS联用方法中排空时间、反应温度、反应时间和色谱分离温度对碳酸盐碳氧同位素分析结果的影响，发现排空时间过短反应瓶中残留的空气会对分析结果造成干扰；反应温度过低会导致酸针堵塞，过高会造成反应瓶瓶盖漏气；反应时间主要受反应温度控制；色谱分离温度对分析结果有一定的影响；通过优化分别选取12 min、72℃、60 min和60℃作为排空时间、反应温度、反应时间和色谱分离温度。与现有的报道相比，本研究更加综合全面地评估了各实验条件对碳酸盐碳氧同位素分析结果的影响。此外，在获得大量分析测试结果的基础上，对同位素数据的归一化方法进行比较，发现双标准物质较准结果明显优于单标准物质。
本研究为GasBenchⅡ-IRMS联用技术中实验条件的选取提供了一定的参考，保证碳氧同位素测试结果的可靠性和准确性。但是，本文主要对标准物质进行系统研究，未对不同性质的实际样品进行考察。而实际样品成分复杂且不均一，因此在开展实际工作中要根据样品的性质进一步对实验条件进行考察。
-
表 1 熔剂的选择
Table 1 Choice of alkali dissolution flux
稀土矿石
样品种类氢氧化钠+过氧化钠 碳酸钠+硼酸
(质量比2∶1)包头矿 清亮 清亮 四川矿 清亮 有不溶物(少量) 独居石 清亮 有不溶物(大量) 
下载: 导出CSV
表 2 熔剂配比的影响
Table 2 Effect of flux ratio
氢氧化钠+
过氧化钠配比空白测定值
ρ(Al)/(μg·mL-1)w(Al2O3)/% 5 g+1 g 0.007 2.03 5 g+2 g 0.007 2.04 4 g+1 g 0.007 2.09 4 g+2 g 0.007 2.08
下载: 导出CSV
表 3 过滤条件选择
Table 3 Choice of filter conditions
w(Al2O3)/% 稀土矿石
样品种类碱分离，过滤洗涤测定 直接定容，
干过滤滤液中
Al2O3的含量滤液中
Al2O3的含量残渣中
Al2O3的含量包头矿 0.40 0.02 0.45 四川矿 1. 47 0.07 1.62 独居石 4.12 0.05 4.34 包头矿标准样品 1.90 0.10 2.09
下载: 导出CSV
表 4 放置时间的选择
Table 4 Choice of standing time
w(Al2O3)/% 放置时间
t/h包头稀土矿标准样品
Al2O3测定值0 2.09 0.25 2.07 0.50 2.09 0.75 2.00 1.0 1.97 2.0 1.81 3.0 1.69 24.0 1.50
下载: 导出CSV
表 5 共存元素的干扰实验
Table 5 Interference test of co-existing ions
分析波长
λ/nm干扰元素限量/% Ca Mg Fe Mn Th P Ni La Ce Pr Nd Sm Eu Gd Tb Dy Ho Er Tm Yb Lu Y 308.2 - - 40 32 0.61 - 89 - - 0 36 - 50 0.4 1.5 - - 0 21 50 - - 309.2 - - 40 56 - - 89 - - - 0 16 12 25 36 36 6 32 - - 13 36 237.3 20 25 - 58 - - - - - 2 48 5 36 - - - 14 - 23 - - 45 394.4 26 80 10 64 - - 23 3 0 0.5 4 3 10 6 - - 13 - - - - - 396.1 40 - 70 50 - - - - 17 2 3 - - 21 25 25 16 - - - - -
下载: 导出CSV
表 6 基体效应
Table 6 Matrix effect
稀土矿石
样品种类基体校正因子 盐分稀释5倍 盐分稀释10倍 盐分稀释20倍 包头矿 0.78 0.91 0.95 四川矿 0.82 0.92 0.95 独居石 0.85 0.91 0.96 铁矿石 0.81 0.91 0.95
下载: 导出CSV
表 7 方法测定下限
Table 7 Detection limits of the method
分析波长
λ/nm标准偏差
s/(mg·g-1)方法检出限
3s/(mg·g-1)测定下限
30s/(mg·g-1)308.2 0.012 0.035 0.35 309.2 0.0096 0.029 0.29 237.3 0.011 0.034 0.34 396.1 0.0070 0.021 0.21
下载: 导出CSV
表 8 方法精密度
Table 8 Precision tests of the method
稀土矿石
样品种类w(Al2O3)/% RSD/% 11次测定值 平均值 包头矿 0.51 0.48 0.46 0.44 0.42 0.49 0.457 7.1 0.47 0.48 0.45 0.42 0.41 四川矿
(加标)1.67 1.62 1.71 1.74 1.79 1.59 1.680 4.4 1.59 1.69 1.60 1.69 1.79 独居石 4.35 4.34 4.38 4.24 4.24 4.35 4.310 1.4 4.30 4.29 4.33 4.21 4.38
下载: 导出CSV
表 9 方法对照实验
Table 9 Comparison of analytical results of Al2O3 by analytical methods
稀土矿石
样品种类w(Al2O3)/% 相对偏差/% ICP-AES法 化学容量法 独居石 4.310 4.280 0.70 包头矿 0.457 0.472 -3.18 四川矿(加标) 1.680 1.697 1.00
下载: 导出CSV
-
冯静.稀土矿石成分分析标准物质的研制[J].化学分析计量,2005,14(4): 1-3,27. http://www.cnki.com.cn/Article/CJFDTOTAL-HXFJ200504000.htm 孙肃,白立忠,于化琴.降低氯化钕溶液中铝离子的工艺研究[J].无机盐工业,2007,39(11): 34-35. doi: 10.3969/j.issn.1006-4990.2007.11.011 程明焱,刘和连,吴伟明,罗飞扬,孙仙源,李安运,陈金清.稀土分析检测方法标准述评[J].有色金属科学与工程,2012,3(4): 108-114. http://www.cnki.com.cn/Article/CJFDTOTAL-JXYS201204021.htm 黎香荣,陈永欣,吕泽娥,罗明贵,谢毓群,刘顺琼,阮贵武.乳化剂增敏铬天青S分光光度法测定铜精矿中的铝[J].分析化学,2009,37(Z1): B142. http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGHY200910003240.htm 熊维巧,张霞铝.试剂示差光度法测定砂岩矿中的三氧化二铝[J].中国非金属矿工业导刊,2005(6): 40-42. http://www.cnki.com.cn/Article/CJFDTOTAL-LGFK200506012.htm 刘建华.分光光度法测定硅铝矿渣中铝的含量[J].武汉化工学院学报,2006,28(2): 20-21. http://www.cnki.com.cn/Article/CJFDTOTAL-WHHG200602005.htm 赵树宝.EDTA络合滴定法连续测定铁矿石中铝铅锌[J].冶金分析,2011,31(11): 66-69. doi: 10.3969/j.issn.1000-7571.2011.11.014 孟亚东,孙洛新,傅晓强.氟盐取代-EDTA滴定法测定铝土矿中铝量[J].岩矿测试,2008,27(6): 475-476. http://www.cnki.com.cn/Article/CJFDTOTAL-YKCS200806016.htm 孙哲平.高铝岩矿中Fe、A1、Ti的EDTA滴定[J].现代科学仪器,2005(5): 71-73. 殷凤玲.铝镁尖晶石中铝、钛、钙、镁的测定[J].矿业快报,2008(12): 110-111. http://www.cnki.com.cn/Article/CJFDTOTAL-KYKB200812041.htm 周尚元.肖伟.铝镁碳砖中三氧二铝的测定[J].湖南冶金,2003,31(1): 44-45. 于永生,王艳蕊,王景霞.珍珠岩矿中Si、Al、Ti含量的测定[J].信阳师范学院学报:自然科学版,2010,23(2): 278-280. http://www.cnki.com.cn/Article/CJFDTOTAL-XYSK201002031.htm 皮业华,光红琼.王啸群.炉渣中铝的测定[J].资源环境与工程,2007,21(5): 613-614. http://www.cnki.com.cn/Article/CJFDTOTAL-HBDK200705026.htm 胡顺峰,王霞,郭合颜,金伟.电感耦合等离子体发射光谱法测定红土镍矿石中镍铬镁铝钴[J].岩矿测试,2011,30(4): 465-468. http://www.cnki.com.cn/Article/CJFDTOTAL-YKCS201104019.htm 张超,李享.电感耦合等离子体发射光谱法测定镍矿石中镍铝磷镁钙[J].岩矿测试,2011,29(4): 473-476. http://www.cnki.com.cn/Article/CJFDTOTAL-YKCS201104021.htm 吕新明,贺国庆,赵晶晶.电感耦合等离子体发射光谱法测定铅精矿中锌、铜、铝、镁含量[J].分析仪器,2010(1): 43. http://www.cnki.com.cn/Article/CJFDTOTAL-FXYQ201001015.htm 孙喜顺,王彦茹,阎雪.电感耦合等离子体原子发射光谱法测定钒钛铁精矿中的钒钛铝镁锰[J].冶金分析,2011,31(8): 79-82. http://www.cnki.com.cn/Article/CJFDTOTAL-YJFX201108021.htm 刘稚,丁仕兵,闵国华.电感耦合等离子体原子发射光谱法测定镍矿中镍铝铬镁钴[J].冶金分析,2008, 28(Z1): 621-623. 金献忠,谢健梅,梁帆,朱丽辉,陈建国.碱熔融-电感耦合等离子体原子发射光谱法测定铬矿石中铬铝铁镁硅[J].冶金分析,2010,30(1): 29-33. http://www.cnki.com.cn/Article/CJFDTOTAL-YJFX201001006.htm
计量
- 文章访问数: 1150
- HTML全文浏览量: 398
- PDF下载量: 15
京公网安备 11010202008159号