-
纳米材料是指在三维空间尺度中，至少有一维的尺度在1—100 nm之间的材料[1]。随着纳米技术的不断发展，各种各样的纳米材料已经广泛应用于生物医药、航空航天、环境保护等很多领域，并发挥着巨大的作用[2-5]。纳米二氧化钛，又称纳米钛白粉，是目前使用较为广泛的无机纳米材料之一[6]。由于纳米二氧化钛具有良好的光催化特性以及紫外线屏蔽能力，其被广泛地应用于化妆品、涂料、油漆和光催化剂等产品中[7-9]。在纳米二氧化钛广泛应用的同时，其对生态环境的负面影响也受到了研究人员的关注[10]。在大量生产、运输、使用和废弃过程中，纳米二氧化钛颗粒会被不可避免地释放到环境中去[11-13]。当纳米二氧化钛被释放到环境介质中时，巨大的比表面积以及超强吸附性能使得它极易吸附环境介质中存在的有毒物质，如持久性有机污染物、重金属等，它可以充当这些有毒有害物质在环境介质中转移的载体或流动相，促进这些污染物在环境介质中的迁移，从而导致污染距离的增大和污染范围的扩大，最终会威胁生态环境安全和人体健康[12,14]。迁移过程是纳米二氧化钛环境地球化学行为中的重要步骤，研究纳米二氧化钛在自然界多孔介质中的迁移行为是揭示其在“固-液”和“固-液-气”界面迁移机制及对其在实际环境介质中的迁移行为做出正确评估的关键。因此，有必要进行纳米二氧化钛在环境介质中的迁移的研究。
近年来，国内外研究工作者关于探究纳米二氧化钛在多孔介质中的迁移规律做了大量的研究。研究表明，影响纳米二氧化钛在多孔介质中迁移行为的因素有很多，如纳米颗粒粒径、浓度、表面性质等纳米材料的性质[15-17]，介质种类、孔隙结构、粒径分布等多孔介质的性质[18-19]，以及pH 值、离子强度、离子组成、表面活性剂等背景溶液的性质[12,19-24]都会影响纳米二氧化钛在多孔介质中的迁移。目前，虽然有关纳米二氧化钛在多孔介质中迁移行为的研究较多，但这些研究大多限于均质的石英砂介质[25]。相比于石英砂，双孔隙介质更为复杂，与真实土壤介质更为接近。双孔隙介质这一概念最早由PIRSON[26]提出用于解释裂缝性油藏原油采收率低的问题，他表示裂缝性储集岩的孔隙度可以认为是由占储集层大部分孔隙空间的颗粒间孔隙度和占小部分孔隙空间的裂缝和裂隙空间组成的两个以上的孔隙-渗透系统组成。之后，PHILIP[27]提出了另一种双孔隙介质类型——聚集质介质（aggregated media）的概念，他表示这种类型介质的孔隙空间由两种不同类型的孔隙度组成：大孔隙 (颗粒间孔隙度)和微孔隙率(颗粒内孔隙率)，并表示聚集体介质构成了聚集成岩土壤的真实表现。纳米二氧化钛在相同环境条件下在石英砂和双孔介质中的迁移规律是否一致，是有待研究的一个问题。目前，关于纳米二氧化钛在双孔隙介质中迁移行为的研究较少。因此，系统研究纳米二氧化钛在双孔隙介质中的迁移过程，对于理解和评估纳米二氧化钛在实际土壤及地下水中的环境风险是非常有必要的。
本研究选取石英砂和双孔隙介质两种多孔介质，探究pH、离子强度（IS）和十二烷基苯磺酸钠（SDBS）的3种环境因素对纳米二氧化钛迁移行为的影响，并通过HYDRUS-1D软件模拟实验结果，预测纳米二氧化钛在两种多孔介质中的迁移机理，为评价纳米二氧化钛的环境风险提供更为全面的理论依据。
-
实验试剂：纳米二氧化钛（99.8%，锐钛矿，亲水），南京先丰纳米材料科技有限公司；氯化钠（NaCl），AR，天津市科密欧化学试剂有限公司；盐酸（HCl，AR）；氢氧化钠（NaOH，AR），天津市天力化学试剂有限公司；十二烷基苯磺酸钠（SDBS C18H29NaO3S，CP），天津市科密欧化学试剂有限公司；石英砂 ( 98% SiO2 ; Fe, Al < 2 %)；双孔隙介质（80% SiO2 ; 18% Ca ; Mg, K, Na < 1 %）；实验用水为超纯水。
实验仪器：TU-1810型紫外可见分光光度计；ISO 9001型蠕动泵；HZT-A2000型电子天平；101 型电热恒温鼓风干燥箱；KQ-50型超声波清洗器；JY92型超声波细胞粉碎机；ME204分析天平。
-
用分析天平称取1 g纳米二氧化钛于500 mL烧杯中，加去离子水后，用超声波细胞粉碎机，325 W功率下超声2 min，倒入1 L的容量瓶中定容，配成1 g·L−1的纳米二氧化钛溶液（储备液A）并用紫外分光光度计测量A0。从储备液A中取出25 mL于1 L容量瓶中，加去离子水稀释定容，配成25 mg·L−1的纳米二氧化钛溶液（储备液B）。
-
用稀释后的盐酸和氢氧化钠溶液调节去离子水的pH，分别配制pH = 2（结合文献[28]，为了探究低pH环境因素对纳米二氧化钛在实验介质中迁移行为的影响）、pH = 6.9和pH = 12（由于石英砂柱实验中pH=12条件下的流出液浑浊，所以纳米二氧化钛数据无法用分光光度法测量，故此文对其不再做讨论）的背景溶液。
-
用电子天平称取58.44 g NaCl晶体于1 L容量瓶中，加去离子水至刻度线定容，配成1 g·L−1的NaCl溶液备用。用移液管移取配置好的NaCl溶液，稀释成1 mmol·L−1、50 mmol·L−1的NaCl背景溶液和50 mmol·L−1的NaCl示踪剂溶液。
-
SDBS的临界胶束浓度（cmc）为1.2 mmol·L−1。首先，称取2.0909 g SDBS固体粉末于500 mL容量瓶内，加去离子水至刻度线定容，配制12 mmol·L−1的SDBS储备液；分别移取0.05 mL、0.5 mL SDBS储备液和1.25 mL纳米二氧化钛储备液A于50 mL容量瓶定容，配制成含SDBS 0.01 cmc、0.1 cmc的纳米二氧化钛悬浊液；取1.25 mL纳米二氧化钛储备液A于50 mL容量瓶定容配制成含0 cmc SDBS的纳米二氧化钛溶液作参照实验。
-
首先对石英砂反复冲洗；第二步用pH = 3—4的盐酸浸泡24 h之后，用去离子水反复冲洗；第三步用pH = 10—11的氢氧化钠浸泡清洗，继续用去离子水清洗至pH呈中性，以除去挥发性有机物；双孔隙介质用去离子水彻底清洗消除其表面的微粒，两种介质都放在烘箱中在105 ℃烘干，放在带盖的桶里储藏，以待实验[29]。
-
采用激光粒度仪（Malvern 2000）来确定两种多孔介质的粒径分布；使用注汞式孔隙度仪（AutoPoreIV9500V1.07）来确定两种多孔介质的孔径分布。利用Philip XL-30环境扫描电子显微镜（ESEM）对多孔介质进行表征。
-
首先将多孔介质颗粒碾碎成接近胶体的细小颗粒，然后将小颗粒悬浮在离子强度为0.1 mmol·L−1 NaCl的溶液中，并在小瓶中剧烈摇动。然后收集一份上清液，采用纳米粒度电位仪Malvern ZS3000 HAS测定两种多孔介质的表面Zeta电位，结果表明，石英砂的表面Zeta电位为 (–39.6 ± 1.8) mV；双孔隙介质的表面Zeta电位为 (–12.5 ± 1.8) mV。
-
装柱：迁移实验使用高度15 cm，直径3.5 cm的玻璃高压层析柱，空柱质量为399.49 g，空柱体积为144.32 cm3。在柱下端放置一张尼龙网以防止多孔介质堵塞出口，将已经洗净、烘干的多孔介质缓慢倒入柱中，轻敲柱壁夯实，重复上述操作至满柱。称取装满多孔介质的柱重，计算出柱中多孔介质质量[30]。
洗柱：连接柱实验装置，从左至右分别为背景溶液、蠕动泵、层析柱、比色皿。调整蠕动泵转速在10—11，恒定流速在3.3—3.7（mL∙min−1）范围，自下而上泵入背景溶液饱和柱子，饱和时间2 h，饱和后拆下柱子称取其质量，计算柱体中含水质量，进而求出水的体积，也即柱的孔隙体积，最后由孔隙体积算出多孔介质的孔隙率。实验测得石英砂柱的孔隙率在0.4—0.5范围内，而双孔隙介质柱的孔隙率在0.8—0.85的范围内。
进样：待柱子饱和后，用移液管移取15 mL的示踪剂或纳米二氧化钛悬浊液于小烧杯中，开始进样的同时用秒表计时，以进样时背景溶液的体积为初始体积，与此同时收集流出液在比色皿中；进样结束时，记录进样时间，继续泵入背景溶液。记录每一次接样时的时间和背景溶液体积，直到收集液中检测不到纳米二氧化钛。
测量：首先，预热紫外-可见分光光度计，之后对光度计进行暗电流矫正；然后在200—400 nm波长范围内，对纳米二氧化钛悬浊液样品进行光谱扫描，得出该浓度悬浊液的紫外-可见扫描光谱图，确定最大波长；在最大波长下，以去离子水为基准线，先测量纳米悬浊液的初始吸光度A1。进样后以收集的第一个样为基准，测量收集流出液样品的吸光度记录数据，利用吸光度与浓度关系计算出相应纳米二氧化钛颗粒浓度。
流出曲线：以孔隙体积数（V/V0）为横坐标，纳米二氧化钛流出浓度与初始浓度比值（C/C0）为纵坐标绘制纳米二氧化钛流出曲线图（BTC）。
-
对纳米二氧化钛在多孔介质中的穿透曲线采用动力学截留模型进行拟合，相应的质量平衡用一维对流－弥散方程来描述：
式中，θ 为水饱和度；C为水相中纳米二氧化钛的浓度；t为时间；ρ为装柱密度；x为纵向坐标；D为水力扩散系数；q为达西流速；s为截留位的固相颗粒浓度。
物质在固液间转移的质量平衡可由下式来表示：
式中，katt和kd分别是吸附和解吸系数；s为固相颗粒浓度；ψ是沉积作用有关的一个功能函数，其动力学可表示为下列公式：
式中，smax为截留位上的最大固相颗粒浓度。dc为介质颗粒的直径；x为与柱子深度有关的截留发生的起点坐标；β为控制颗粒空间分布的一个经验参数。Bradford等发现对于与柱子深度有关的截留作用取0.432最优[15]。
-
颗粒之间的相互作用力一般采用DLVO 理论来描述，该理论将颗粒间的总势能Φtot表示为范德华力ΦVDW和静电作用力ΦEL之和：
通常情况下，总势能Φtot、范德华力ΦVDW和静电作用力ΦEL均除以玻尔兹曼常数（k = 1.38 ×10−23 J·K−1）和绝对温度（T），即用无量纲形式来表示，故Φ以kT为度量。
粒子之间的静电作用可用下式来表示：
式中，ε = 78.5，为水的相对介电常数；ε0 = 8.854×10−12 C2·N−1·m2为真空绝对介电常数；rp是二氧化钛的平均粒径，ψ01和ψ02分别是二氧化钛和多孔介质的表面Zeta电势；h是二氧化钛和多孔介质颗粒表面之间的距离，κ是德拜参数。
粒子间的范德华作用力可由下式表示：
式中，λ为特征波长，一般取100 nm；A = 1×10−20为Hamaker常数。
-
实验用石英砂的密度为2.65 g·cm−3。由图1a可知，石英砂的粒径分布在0.16—0.79 mm之间，其平均粒径为0.36 mm；由图1b可知，石英砂的孔径分布为5—200 μm，孔径分布图只有1个峰，且峰值为55 μm。
实验用双孔隙介质的密度为2.49 g·cm−3。由图1a可知，双孔隙介质的粒径在0.4—5.0 mm之间，平均粒径为1.5 mm；由图1b可以看出，双孔隙介质孔隙分布复杂，孔隙分布图有两个孔径范围，1个在0.005—5 μm，峰值为0.035 μm；1个是孔隙范围在5—360 μm，有两个峰值，15 μm和200 μm。
从孔径分布图可以看出，相较于石英砂的单峰孔径分布图，双孔隙介质的孔径分布图，呈多范围、多峰值的特点，说明双孔隙介质内部孔道更为发达，内部结构更为复杂。图1c、图1d分别是石英砂和双孔隙介质的扫描电镜图。从图像中可以清晰地看到，石英砂颗粒表面均匀，颗粒内部无孔隙；而双孔隙介质表面粗糙，颗粒内部有很多大小不一的小孔。
-
采用50 mmol·L−1的NaCl作为示踪剂，分别对石英砂和双孔隙介质两种不同多孔介质中的水流特征进行表征。示踪剂在石英砂和多孔介质中的实验条件如表1所示，石英砂在柱子里的堆积密度为1.42 g·cm−3，而双孔隙介质的堆积密度是0.44 g·cm−3. 石英砂在柱子中的孔隙度为0.47，而双孔隙介质在柱子中的孔隙度为0.83。相比于石英砂，双孔隙介质的堆积密度很小，孔隙体积很大，这与前文介质表征时所表明的双孔隙介质内部由复杂的孔隙结构密切相关。
图2是示踪剂在两种多孔介质中的穿透曲线。由图2可以发现，示踪剂穿过石英砂的时候，穿透曲线比较对称，这说明水流经过石英砂的时候流动状态比较均匀；而示踪剂穿过双孔隙介质的时候，穿透曲线是不对称的，有明显拖尾现象，这说明水流经过双孔隙介质的时候，由于双孔隙介质内部孔隙复杂不均，导致流动状态是不均匀的。表2是拟合参数表，示踪剂实验中的实验现象与HYDRUS模拟软件得出的拟合参数相一致：示踪剂穿过装有石英砂的柱中，水流扩散系数是0.19 cm，而穿过装有双孔隙介质的柱中，水流扩散系数为1.95 cm。
-
图3a是不同pH条件下，实验所得纳米二氧化钛在饱和石英砂中的穿透曲线及模拟曲线。从图3可以看出，不同pH条件下，纳米二氧化钛在石英砂中的穿透曲线都有峰，表明两种条件下均有二氧化钛从饱和石英砂中流出；然而，相比于pH = 6.9时的穿透曲线，pH = 2时的穿透曲线峰值很低，与横坐标轴围成的图形较小，表明pH = 2时饱和石英砂中纳米二氧化钛流出较少。
由实验条件表表3可知，在饱和石英砂介质中，pH = 2时，纳米二氧化钛的回收率为1.7%；而在pH = 6.9时，回收率为30%。
通过HYDRUS-1D软件，对实验进行了多次模拟，得到实验模拟结果参数表（表4）。由表4可知，参数katt的值在pH = 2时远大于pH = 6.9时的值。结合两点动力学理论分析，在两点动力学理论中，katt越大说明介质表面对胶体颗粒的吸附能力越强，因此在实验条件下，pH = 2时石英砂对纳米二氧化钛的吸附能力要远大于pH = 6.9时的吸附能力，所以导致纳米二氧化钛大量滞留不易流出。
图4a为不同pH条件下石英砂与纳米二氧化钛之间的能势图。结合DLVO理论分析，由图4可知pH = 2时，纳米二氧化钛颗粒与石英砂表面之间的总势能Φtot小于0，且随着纳米二氧化钛颗粒与石英砂表面的距离增加，总势能Φtot的值不断趋近于0 KT。Φtot小于0，表明石英砂表面与纳米二氧化钛颗粒之间成吸附作用，与模拟结果相一致。相反地，pH = 6.9时，在纳米二氧化钛颗粒距离石英砂介质表面0—20 nm范围内，能势曲线存在一个峰，峰值为200 KT，表明纳米二氧化钛与石英砂之间表现为静电斥力作用，能垒较大，导致纳米二氧化钛不易吸附在石英砂表面，从而导致pH = 6.9时纳米二氧化钛从饱和石英砂中的流出能力高于pH = 2。据文献表述[21]，纳米二氧化钛的等电点为4.4—6.2，当pH小于等电点时，未被包被的二氧化钛会在带负电的介质中大量滞留，实验结论与其相一致。
图3b为不同pH条件下，纳米二氧化钛在双孔隙介质中的穿透曲线图，与石英砂相似，pH = 2时纳米二氧化钛在双孔隙介质中的穿透曲线与坐标轴围成的图形要小于pH = 6.9时的穿透曲线与坐标轴围成的图形。由表3可知，与在石英砂中相似，pH = 2时纳米二氧化钛在双孔隙介质中的回收率为2.6%，远小于pH = 6.9时的50.8%。由表4可知，与在石英砂中模拟结果相似，pH = 2时的katt为19.4，远大于pH = 6.9时的0.753。图4b是纳米二氧化钛与双孔隙介质的能势图，与在石英砂中情况类似，pH = 2时的纳米二氧化钛与双孔隙介质表面的总势能Φtot小于0 kT，使得纳米二氧化钛与双孔隙介质表面成吸附作用；而pH = 6.9时的Φtot则大于0 kT，且存在一个能量壁垒，使得纳米二氧化钛与双孔隙介质表面成静电斥力作用，导致pH = 2时纳米二氧化钛在双孔隙介质中穿透能力小于pH = 6.9时。
结论虽相似，然细节有不同。相比于纳米二氧化钛在石英砂中的穿透曲线，所有pH条件下，纳米二氧化钛在双孔隙介质中的穿透曲线对称度较低且有明显的拖尾现象；相比于饱和石英砂中纳米二氧化钛的迁移回收率，饱和双孔隙介质中纳米二氧化钛的迁移回收率较大。从表4参数对比中可以看出，相同pH条件对比下，饱和石英砂中的katt值均大于饱和双孔隙介质中的katt值，表示实验中相同pH条件下石英砂对纳米化钛的吸附能力要大于双孔隙介质，造成纳米二氧化钛在饱和双孔隙介质中的迁移能力较强，回收率较大。从能势图对比可以看出，pH = 2的条件下，纳米二氧化钛与双孔隙介质的初级能阱绝对值为60 kT，小于纳米二氧化钛与石英砂的初级能阱的100 KT，表示pH = 2时，纳米二氧化钛与石英砂的吸附作用大于纳米二氧化钛与双孔隙介质的吸附作用，与前述模拟结果相符；然而pH = 6.9的条件下，纳米二氧化钛与石英砂表面的能量势垒峰值为200 kT，大于与双孔隙介质的能量壁垒峰值20 kT，表示pH = 6.9条件下，纳米二氧化钛与石英砂表面的静电斥力作用要大于与双孔隙介质表面的静电斥力作用，以此推测，纳米二氧化钛在石英砂中的迁移能力要大于在饱和双孔隙介质中的迁移能力，与前述结论不符。纳米颗粒的团聚作用是造成pH = 6.9条件下与结论不符的原因，由于静电斥力作用，纳米二氧化钛颗粒会发生团聚作用，导致纳米二氧化钛颗粒粒径增大，静电斥力作用越大，团聚作用越剧烈。纳米二氧化钛与石英砂表面的静电斥力作用远大于与双孔隙介质间的静电斥力作用，导致纳米二氧化钛在石英砂中的团聚作用较强，从而迁移能力减弱；另一方面，有文献报道，纳米颗粒在孔径比较大的多孔介质中的回收率要大于孔径小的回收率[31]，由介质表征可知双孔隙介质的孔径要大于石英砂介质的孔径，从而，相比之下，双孔隙介质中本就较弱的团聚作用的影响也不明显。
综上所述，pH对纳米二氧化钛在双孔隙介质中迁移行为的影响类似于纳米二氧化钛子在饱和石英砂中的迁移行为的影响，然而由于双孔隙介质自身孔隙较大的特性，相同pH条件下的纳米二氧化钛在双孔隙介质中的回收率要大于纳米二氧化钛在石英砂中的回收率。
-
图5a是在不同离子强度下，纳米二氧化钛在石英砂中的穿透曲线。可以看出，纳米二氧化钛均可穿透饱和石英砂柱，且低离子强度(1 mmol∙L−1)流出曲线在上，高离子强度(50 mmol∙L−1)流出曲线在下。由表3可知高离子强度下纳米二氧化钛的回收率为64.2%，小于低离子强度下的67.6%。由表4可知，高离子强度的katt参数为8.73，而低离子强度的katt值是0.316。高离子强度的katt值比低离子强度大了一个数量级，结合两点动力学理论可知高离子强度下石英砂对纳米二氧化钛的吸附能力较强，因而造成纳米二氧化钛的不易迁移，回收率较低。
同时，图6a是石英砂表面与纳米颗粒间的能势图，高离子强度条件下的能势图在下，总势能值始终为负值；低离子强度条件下的势能图在上，存在能垒，总势能值始终为正值；结合DLVO理论可知，高离子强度下石英砂表面与纳米二氧化钛颗粒之间呈吸附作用，与前述结论一致。
造成这一现象的原因是，外来离子会对纳米颗粒表面的双电层产生影响，高离子条件下，纳米二氧化钛双电层被压缩，石英砂与纳米二氧化钛颗粒之间的静电斥力作用减弱，吸附作用增强，从而抑制了纳米二氧化钛的迁移，造成高离子强度下纳米二氧化钛的回收率较低。
在饱和双孔隙介质中，不同离子强度下，纳米二氧化钛的穿透曲线如图5b。相比于饱和石英砂，不同离子条件下双孔隙介质中纳米二氧化钛的穿透曲线有明显的拖尾现象，并且，双孔隙中纳米二氧化钛的初始流出时间要早于石英砂中纳米二氧化钛的流出时间。由表3可知，高离子强度下纳米二氧化钛的回收率为42.4%，低离子强度条件下的回收率为47.1%。与饱和石英砂条件下类似，饱和双孔隙介质中，高离子强度条件下纳米二氧化钛的回收率小于低离子强度条件下的回收率。然而，与石英砂相比，在所有离子强度下，双孔隙介质条件下的纳米二氧化钛的回收率要小于石英砂中纳米二氧化钛的回收率。
由表4可知，与石英砂条件下类似，双孔隙介质中，高离子强度下的katt值大于低离子强度下的katt值，表明高离子强度下双孔隙介质对纳米二氧化钛颗粒的吸附作用大于低离子强度的吸附作用从而与穿透曲线结果相符。但是，相较于石英砂低离子强度条件下的katt值，双孔隙介质低离子强度条件下的katt值较大，但相差不大；同时，相比于石英砂高离子强度条件下的katt值，双孔隙介质高离子强度下的katt值较小。以此结合两点动力学推测，不同离子条件下应是石英砂条件下的回收率较低或者两者差别不大。但是，相同离子条件下，石英砂中的kdet值远大于双孔隙介质中的kdet值，结合两点动力学理论，kdet越大表示纳米颗粒的脱附能力越强。上述表明，不同离子条件下，纳米二氧化钛在石英砂中的脱附能力远大于在双孔隙介质中的脱附能力，从而使得其在石英砂中的回收率远大于在双孔隙介质中的回收率；而此时，纳米二氧化钛在不同介质对的吸附能力大小对纳米颗粒的迁移影响甚微。
图6b是双孔隙介质与纳米二氧化钛颗粒间的能势图，与石英砂中相似，低离子强度条件下，双孔隙介质表面对纳米二氧化钛呈静电斥力作用，且存在能量壁垒；而相反，高离子强度下的能势曲线在下，且小于0 kT，表示此时双孔隙介质对纳米二氧化钛成吸附作用。但是，相比于石英砂，双孔隙介质低离子强度下最大势能绝对值较小，表示此时双孔隙介质对纳米二氧化钛的静电斥力作用小于石英砂对纳米二氧化钛的静电斥力作用，导致相比之下，纳米二氧化钛更易在双孔隙介质中滞留，从而使得低离子强度下双孔隙介质中的回收率较低。纳米颗粒与介质表面同一距离下，双孔隙介质高离子强度下的总势能绝对值与石英砂中相似，表明此时两者对纳米二氧化钛颗粒的吸附能力相同，结合实验模拟分析，在吸附能力相近的情况下，纳米二氧化钛在石英砂中的脱附能力较大，故而纳米二氧化钛在石英砂中的迁移能力较强，回收率高。
综上所述，离子强度对两种介质中纳米二氧化钛迁移行为的影响是相似的，都是高离子强度抑制，低离子强度促进纳米二氧化钛的迁移。但是，由于纳米二氧化钛在两种介质中的吸附和脱附能力的不同，导致相同离子强度条件下，双孔隙介质表现为抑制纳米二氧化钛的迁移。
-
图7a是纳米二氧化钛在饱和石英砂中的穿透曲线。
由图7可知，实验条件下不同浓度的SDBS、纳米二氧化钛均能从石英砂中穿透，并且从3种SDBS条件下纳米二氧化钛均是从0.5个孔隙体积左右开始流出，之后在1.25孔隙体积流出浓度达到峰值。在开始的1.25个孔隙体积数内相对流出浓度从高到低分别对应SDBS的含量为 0.1 cmc、0.01 cmc、0 cmc，但是在1.25个孔隙体积数后，从高到低的顺序变为了0.01 cmc、0.1 cmc、0 cmc。介质表征可知石英砂的表面光滑，所以其对SDBS的吸附能力较弱，所以在最初1.25个孔隙体积数时，表面活性剂对的纳米二氧化钛流出曲线的作用主要体现在纳米二氧化钛和介质之间的静电排斥力，体现到曲线上就是1 cmc的体系中纳米二氧化钛的相对流出浓度较高。
由表3可知，随着SDBS含量的增加，纳米二氧化钛的回收率也越来越大，表明随着SDBS的增加，纳米二氧化钛在石英砂中的迁移能力越来越强。SDBS是一种典型的阴离子表面活性剂，它可以吸附于纳米二氧化钛表面或者吸附于石英砂表面，从而改变石英砂与二氧化钛之间的空间结构及相互作用力，从而影响纳米二氧化钛在介质中的迁移。从表4可以看出，随着SDBS浓度的增大，对应的katt值不断减小，结合两点动力学理论，katt值减小说明纳米二氧化钛的吸附能力减弱，从而其迁移能力不断增强，导致回收率不断增加。Godinez [23]研究表明，SDBS的存在和浓度的增加增强了纳米二氧化钛在饱和多孔介质中的迁移，表面活性剂对空间斥力的影响使纳米二氧化钛保持稳定或分散，从而增加了迁移率。
图7b是双孔隙介质中不同SDBS存在条件下的纳米二氧化钛的穿透曲线图。由图7可知，相比于石英砂的流出曲线，同一SDBS浓度下，双孔隙介质中纳米二氧化钛的迁移曲线具有较为明显的拖尾现象，同时随着表面活性剂浓度的升高拖尾越严重。由图7可知，双孔隙介质中，纳米二氧化钛从0.1孔隙体积开始流出，比石英砂的0.5孔隙体积流出较快。在开始的0.5个孔隙体积数内，纳米二氧化钛在双孔隙介质柱的相对流出浓度从高到低分别对应的含SDBS浓度为 0.01 cmc、0 cmc、0.1 cmc，在0.5个孔隙体积数后，从高到低的顺序变为了0.1cmc、0.01cmc、0 cmc。这是由于在最初的0.5个孔隙体积数内，SDBS的桥键作用起了主导优势，SDBS浓度越高，桥键作用越大，因此0.1 cmc体系中纳米二氧化钛的相对流出浓度最低；随着时间的推移，淋溶液中阴离子表面活性剂SDBS占据了介质表面的正电荷位点，随着SDBS在介质上的吸附量逐渐增加，纳米二氧化钛和介质之间的静电排斥力增加且随着SDBS浓度的增加更加明显。因此在0.5个孔隙体积数后，静电排斥力的优势更明显，使含高浓度SDBS体系中的纳米二氧化钛更不容易被多孔介质截留，所以最后0.1 cmc SDBS的体系中纳米二氧化钛的相对流出浓度最高。由表3可知与石英砂中相似，纳米二氧化钛的回收率随着SDBS浓度的增加而增加。然而，相比于石英砂中纳米二氧化钛的流出率，同一SDBS浓度下，双孔隙介质中纳米二氧化钛的回收率较小。从表4可以看出随着SDBS的浓度增加，katt值逐渐降低，表示纳米二氧化钛的吸附能力逐渐降低。
结合表3、表4分析，SDBS存在条件下，双孔隙介质实验参数模拟结果的katt值与石英砂实验参数模拟结果相比差距不大，结合两点里动力学分析，此时两种介质表面对纳米二氧化钛颗粒的吸附能力相当；然而相对应的双孔隙介质实验参数模拟的的kdet值要比石英砂的小几个数量级，表明此时虽然吸附能力相当，但是纳米二氧化钛颗粒在石英砂中的脱附能力要远大于在双孔隙介质中的脱附能力，从而导致同一SDBS浓度下，纳米二氧化钛在石英砂中的滞留较少，回收率较高。
综上所述，无论多孔介质是均质的石英砂还是非均质的双孔隙介质，SDBS的加入，都会大幅度提高纳米二氧化钛的迁移能力。然而，相比于在石英砂中而言，纳米二氧化钛在双孔隙介质中的回收率较小。
-
(1) pH对纳米二氧化钛在两种多孔介质中迁移行为的影响与纳米二氧化钛等电点有关，当pH低于等电点时，纳米二氧化钛会在石英砂和双孔隙中大量滞留。
(2) 由于双孔隙介质对纳米二氧化钛的吸附能力较弱以及本身孔隙较多的综合原因，使得相同pH条件下，纳米二氧化钛在双孔隙介质中比在细砂中的迁移能力较强，回收率较高。
(3) IS对纳米二氧化钛在两种介质中迁移行为的影响相似，均为：高离子强度抑制、低离子强度促进纳米二氧化钛的迁移。
(4) 由于纳米二氧化钛在双孔隙介质中比在细砂介质中的脱附能力弱，导致相比于细砂介质，在相同离子强度条件下，双孔隙介质表现为抑制纳米二氧化钛的迁移。
(5) SDBS的存在促进了纳米二氧化钛在两种多孔介质中的迁移，且随着SDBS浓度的增加，纳米二氧化钛的迁移能力增强，回收率增加。
(6) 由于纳米二氧化钛颗粒在石英砂中的脱附能力要远大于在双孔隙介质中的脱附能力，从而导致相较于细砂介质，在同一SDBS浓度下，纳米二氧化钛在双孔隙介质中的滞留较多，回收率较低。
不同环境因素对纳米二氧化钛在双孔隙介质中迁移行为的影响
TiO2 nanoparticles transport in aggregated porous media with intra- and interaggregate porosities: Effects of pH, ionic strength and sodium dodecylbenzene sulfonate
-
摘要: 选取石英砂和双孔隙介质两种多孔介质，通过柱实验探究pH、离子强度（IS）和阴离子表面活性剂十二烷基苯磺酸钠（SDBS）的3种环境因素对纳米二氧化钛在饱和多孔介质中迁移行为的影响。通过HYDRUS-1D用两点动力学模型对实验结果进行模拟，并结合DLVO理论进行分析。结果表明，pH对纳米二氧化钛在两种多孔介质中迁移行为的影响与纳米二氧化钛等电点有关，当pH低于等电点时，纳米二氧化钛会在石英砂和双孔隙中大量滞留；由于双孔隙介质孔隙较多，且对纳米二氧化钛的吸附能力较弱，因而在相同pH下，纳米二氧化钛在双孔隙介质中比在细砂中的迁移能力较强，回收率较高；IS对纳米二氧化钛在两种介质中迁移行为的影响相似，均为高离子强度抑制、低离子强度促进纳米二氧化钛的迁移；由于纳米二氧化钛在双孔隙介质中比在细砂介质中的脱附能力弱，导致相比于细砂介质，在相同离子强度条件下，双孔隙介质表现为抑制纳米二氧化钛的迁移；SDBS的存在促进了纳米二氧化钛在两种多孔介质中的迁移，且随着SDBS浓度的增加，纳米二氧化钛的迁移能力增强，回收率增加；由于纳米二氧化钛颗粒在石英砂中的脱附能力要远大于在双孔隙介质中的脱附能力，从而导致相较于细砂介质，在同一SDBS浓度下，纳米二氧化钛在双孔隙介质中的滞留较多，回收率较低。Abstract: The influences of pH, ionic strength (IS) and sodium dodecylbenzene sulfonate (SDBS) on TiO2 nanoparticles transport behavior in saturated porous media were investigated by column experiment, and by selecting two kinds of porous media, i.e. quartz sand and aggregated porous media with intra- and interaggregate porosities (dual-porosity media). Experimental results were simulated using two-site model by HYDRUS-1D code and analyzed theoretically according to DLVO theory. The results showed that the effect of pH on the transport of TiO2 nanoparticles in both porous media was related to the isoelectric point of TiO2 nanoparticles. When pH was lower than the isoelectric point, TiO2 nanoparticles would be largely retained in both media. Due to the large number of pores and the weak adsorption capacity of the dual-porosity media for TiO2 nanoparticles and the large number of pores, the TiO2 nanoparticles in the dual-porosity media had a stronger transport capacity and a higher recovery rate under the same pH, compared to fine sand. Given to the effect of IS, high IS inhibited the transport of TiO2 nanoparticles, and low IS promoted the transport of TiO2 nanoparticles for both porous media. However, the transport of TiO2 nanoparticles was inhibited in the dual-porosity media under the same ionic strength compared to fine sand, because the desorption ability of TiO2 nanoparticles from the dual-porosity media was weak. The presence of SDBS promoted the transport of TiO2 nanoparticles in both porous media. Furthermore, with the increase of SDBS concentration, the transport ability of TiO2 nanoparticles was enhanced, and the recovery rate increased correspondingly. Because the desorption ability of TiO2 nanoparticles in fine sand was much greater than that in dual-porosity media, TiO2 nanoparticles remained more in dual-porosity media with the same SDBS concentration, and the recovery rate was correspondingly lower.
-
Key words:
- TiO2 nanoparticle /
- saturated porous media /
- transport
-
-
图 1 两种多孔介质的粒径分布(a)，孔径分布(b)和扫描电镜图（c:石英砂；d:双孔隙介质）
Figure 1. Particle size distribution (a), pore size distribution (b) and Scanning electron microscopy (c: quartz sand; d: dual-porosity media)
图 2 示踪剂在石英砂和双孔隙介质中穿透曲线和拟合曲线
Figure 2. Penetration curves and fitting curves of the tracer in quartz sand and dual-porosity media
图 3 不同pH下纳米二氧化钛在石英砂(a)和双孔隙介质(b)中的穿透曲线
Figure 3. Penetration curves of TiO2 nanoparticles in quartz sand (a) and dual-porosity media (b) at different pH
图 4 pH对纳米颗粒与石英砂(a)和双孔隙介质(b)间势能的影响
Figure 4. Influence of pH on potential energy between nanoparticles and quartz sand (a) and dual-porosity media (b)
图 5 不同IS下纳米二氧化钛在石英砂(a)和双孔隙介质(b)中的穿透曲线
Figure 5. Penetration curves of TiO2 nanoparticles in quartz sand (a) and dual-porosity media (b) under different IS
图 6 IS对纳米颗粒与石英砂(a)和双孔隙介质(b)间势能的影响
Figure 6. Effect of IS on potential energy between nanoparticles and quartz sand (a) and dual-porosity media (b)
图 7 不同浓度SDBS下纳米二氧化钛在石英砂(a)和双孔隙介质(b)中的穿透曲线
Figure 7. Penetration curves of TiO2 nanoparticles in quartz sand (a) and dual-porosity media (b) at different concentrations of SDBS
表 1 示踪剂实验条件表
Table 1. Tracer experimental conditions
多孔介质
Porous media介质堆积密度/（g·cm−3）
Bulk density孔隙度Porosity 流速/（cm·min−1）
Darcy velocity注入时间/min
Injection time回收率/%
Recovery石英砂
Quartz sand1.42 0.47 0.37 4.28 99.4 双孔隙介质
Dual-porosity media0.44 0.83 0.36 4.32 99.9 
下载: 导出CSV
表 2 示踪剂实验参数模拟结果表
Table 2. Simulation results of tracer experimental parameters
多孔介质
Porous media扩散系数/cm
Diffusion coefficient扩散速率/min−1
Diffusion rateR2 石英砂
Quartz sand0.19 1.03×10−4 0.9961 双孔隙介质
Dual-porosity media1.95 3.4×10−1 0.9973
下载: 导出CSV
表 3 柱实验条件表
Table 3. Column experimental conditions
多孔介质
Porous media环境因素
Environmental conditions装柱密度/（g·cm−3）
Bulk density孔隙度/%
Porosity达西流速/（cm·min−1 ）
Darcy velocity注入时间/min
Injection time回收率/%
Recovery滞留率/%
Detention石英砂
Quartz sandpH 2 1.33 0.5 0.37 4.28 1.7 98.3 6.9 1.35 0.49 0.35 4.4 30 70 IS 1 1.42 0.47 0.38 4.12 67.6 32.4 50 1.49 0.44 0.36 4.33 64.2 35.8 SDBS 0 1.35 0.49 0.35 4.4 30 70 0.01 1.4 0.47 0.43 3.97 74.2 25.8 0.1 1.41 0.47 0.39 4.32 97.5 2.5 双孔隙介质
Dual-porosity mediapH 2 0.42 0.83 0.32 4.82 2.6 97.4 6.9 0.43 0.83 0.36 4.32 50.8 49.2 IS 1 0.45 0.82 0.38 4.12 47.1 52.9 50 0.46 0.82 0.35 4.55 42.4 57.6 SDBS 0 0.43 0.83 0.36 4.32 50.8 49.2 0.01 0.44 0.82 0.36 4.37 67 33 0.1 0.45 0.82 0.35 4.48 78.1 21.9
下载: 导出CSV
表 4 柱实验模拟结果参数表
Table 4. Parameters of column experiment simulation results
多孔介质
Porous media环境因素
Environmental conditions模拟结果
Simulation resultskatt kdet R2 石英砂
Quartz sandpH 2 1.31×102 1.23×10−2 0.7872 6.9 1.28 3.95×10−2 0.7493 IS 1 0.316 1.34 0.9457 50 8.73 1.42 0.8506 SDBS 0 1.28 3.95×10−2 0.7493 0.01 0.228 0.194 0.9685 石英砂
Quartz sandSDBS 0.1 0.122 7.40×10−3 0.9939 双孔隙介质
Dual-porosity mediapH 2 19.4 3.50×10−4 0.9515 6.9 0.753 2.36×10−6 0.9467 IS 1 0.999 1.41×10−6 0.9828 50 1.77 3.92×10−5 0.9609 SDBS 0 0.753 2.36×10−6 0.9467 0.01 0.319 5.37×10−5 0.9647 0.1 0.261 1.39×10−5 0.965
下载: 导出CSV
-
[1] 沈钟, 赵振国, 康万利. 胶体与表面化学[M]. 4版. 北京: 化学工业出版社, 2012: 329. SHEN Z, ZHAO Z G, KANG W L. Colloid and surface chemistry[M]. Beijing: Chemical Industry Press, 2012: 329(in Chinese).
[2] 于荣丽, 孙铁珩, 胡晓钧. 纳米材料的生物负效应研究进展 [J]. 安全与环境学报, 2009, 9(2): 121-124. YU R L, SUN T H, HU X J. Research progress on the negative effects of nanomaterials on biological systems [J]. Journal of Safety and Environment, 2009, 9(2): 121-124(in Chinese).
[3] KHAN S, NAUSHAD M, AL-GHEETHI A, et al. Engineered nanoparticles for removal of pollutants from wastewater: Current status and future prospects of nanotechnology for remediation strategies [J]. Journal of Environmental Chemical Engineering, 2021, 9(5): 106160. doi: 10.1016/j.jece.2021.106160 [4] 李冉, 吴霜, 付鹏辉. 浅谈纳米技术在生活生产中的应用 [J]. 信息记录材料, 2019, 20(5): 14-15. LI R, WU S, FU P H. Application of nanotechnology in life production [J]. Information Recording Materials, 2019, 20(5): 14-15(in Chinese).
[5] 冯辰昀, 李旭东, 郑妤婕, 等. 纳米材料的毒理学研究进展[J/OL]. 中国科学: 化学: 2021, 51,doi: 10.1360/SSC-2021-0138. FENG C Y, LI X D, ZHENG Y J, et al. Recent progress in nanotoxicology of nanomaterials[J/OL]. Sci Sin Chim, 2021, 51,doi: 10.1360/SSC-2021-0138 (in Chinese).
[6] CHEN X B, MAO S S. Titanium dioxide nanomaterials: Synthesis, properties, modifications, and applications [J]. Chemical Reviews, 2007, 107(7): 2891-2959. doi: 10.1021/cr0500535 [7] 刁润丽, 赵世伟. 纳米二氧化钛的应用研究进展 [J]. 山西化工, 2021, 41(3): 25-26,31. DIAO R L, ZHAO S W. Application research progress of Nano-titanium dioxide [J]. Shanxi Chemical Industry, 2021, 41(3): 25-26,31(in Chinese).
[8] 张海丰, 张鹏宇, 赵贵龙, 等. 纳米二氧化钛的制备及其应用研究进展 [J]. 东北电力大学学报, 2014, 34(2): 52-56. doi: 10.3969/j.issn.1005-2992.2014.02.011 ZHANG H F, ZHANG P Y, ZHAO G L, et al. Advances on preparation, modification and application of nanosized titanium dioxide [J]. Journal of Northeast Dianli University, 2014, 34(2): 52-56(in Chinese). doi: 10.3969/j.issn.1005-2992.2014.02.011
[9] 刁润丽. 纳米二氧化钛的光催化性及其应用研究进展 [J]. 佛山陶瓷, 2021, 31(6): 5-7. DIAO R L. Research progress on photocatalytic activity and application of nano-titanium dioxide [J]. Foshan Ceramics, 2021, 31(6): 5-7(in Chinese).
[10] 李润虎. 纳米技术的风险问题及对策研究 [J]. 工程研究-跨学科视野中的工程, 2020, 12(4): 380-387. LI R H. Research on risks and countermeasures of nanotechnology [J]. Journal of Engineering Studies, 2020, 12(4): 380-387(in Chinese).
[11] 刘诚, 孙强, 张刚, 等. 磷酸盐对纳米TiO2在土壤中迁移的影响 [J]. 苏州科技大学学报(自然科学版), 2018, 35(3): 44-50. LIU C, SUN Q, ZHANG G, et al. Influence of phosphate on the transport of nTiO2 in soil column [J]. Journal of Suzhou University of Science and Technology (Natural Science Edition), 2018, 35(3): 44-50(in Chinese).
[12] 韩鹏, 王雪婷. 不同环境因子对纳米二氧化钛在饱和多孔介质中迁移持留行为的影响 [J]. 应用基础与工程科学学报, 2013, 21(3): 512-521. doi: 10.3969/j.issn.1005-0930.2013.03.013 HAN P, WANG X T. Transport and retention behaviors of TiO2 nanoparticles in saturated porous media: Effects of ionic strength, iron(hydr)oxide and humic acid [J]. Journal of Basic Science and Engineering, 2013, 21(3): 512-521(in Chinese). doi: 10.3969/j.issn.1005-0930.2013.03.013
[13] 彭爱夏, 占敬敬, 吴明火. 纳米粒子在多孔介质中迁移模型的优化 [J]. 化工学报, 2021, 72(10): 5114-5122. PENG A X, ZHAN J J, WU M H. Optimization of nanoparticles transport model in porous media [J]. CIESC Journal, 2021, 72(10): 5114-5122(in Chinese).
[14] 方婧, 余博阳. 3种金属氧化物纳米材料在不同土壤中运移行为研究 [J]. 环境科学, 2013, 34(10): 4050-4057. FANG J, YU B Y. Transport behaviors of metal oxide nanomaterials in various soils [J]. Environmental Science, 2013, 34(10): 4050-4057(in Chinese).
[15] 张柯柯. 纳米TiO2在土壤中迁移及其与Pb共迁移机制研究[D]. 杭州: 浙江工商大学, 2015: 72. ZHANG K K. Transport of nano-TiO2 in soil and the CO-transport mechanisms of Pb with nano-TiO2[D]. Hangzhou: Zhejiang Gongshang University, 2015: 72(in Chinese).
[16] FANG J, SHAN X Q, WEN B, et al. Stability of titania nanoparticles in soil suspensions and transport in saturated homogeneous soil columns [J]. Environmental Pollution, 2009, 157(4): 1101-1109. doi: 10.1016/j.envpol.2008.11.006 [17] DUNPHY GUZMAN K A, FINNEGAN M P, BANFIELD J F. Influence of surface potential on aggregation and transport of titania nanoparticles [J]. Environmental Science & Technology, 2006, 40(24): 7688-7693. [18] CAI L, PENG S N, WU D, et al. Effect of different-sized colloids on the transport and deposition of titanium dioxide nanoparticles in quartz sand [J]. Environmental Pollution, 2016, 208: 637-644. doi: 10.1016/j.envpol.2015.10.040 [19] 罗小廷, 吴丹, 梁嘉良, 等. 典型阴离子对纳米二氧化钛在载铁石英砂多孔介质中迁移行为的影响 [J]. 北京大学学报(自然科学版), 2017, 53(4): 749-757. LUO X T, WU D, LIANG J L, et al. Influence of typical anions on the transport of titanium dioxide nanoparticles in iron oxide-coated porous media [J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 2017, 53(4): 749-757(in Chinese).
[20] 吴其圣, 杨琛, 胡秀敏, 等. 环境因素对纳米二氧化钛颗粒在水体中沉降性能的影响 [J]. 环境科学学报, 2012, 32(7): 1596-1603. WU Q S, YANG C, HU X M, et al. Influences of environmental factors on aggregation of titanium dioxide nanoparticles [J]. Acta Scientiae Circumstantiae, 2012, 32(7): 1596-1603(in Chinese).
[21] 张瑞昌, 章海波, 涂晨, 等. 腐殖酸作用下酸性多孔介质中纳米TiO2的迁移与滞留机制 [J]. 中国环境科学, 2018, 38(9): 3542-3551. doi: 10.3969/j.issn.1000-6923.2018.09.041 ZHANG R C, ZHANG H B, TU C, et al. Mechanisms of mobility and retention of nano-TiO2 in acidic porous media in the presence of humus acids [J]. China Environmental Science, 2018, 38(9): 3542-3551(in Chinese). doi: 10.3969/j.issn.1000-6923.2018.09.041
[22] SUN P D, ZHANG K K, FANG J, et al. Transport of TiO2 nanoparticles in soil in the presence of surfactants [J]. Science of the Total Environment, 2015, 527/528: 420-428. doi: 10.1016/j.scitotenv.2015.05.031 [23] GODINEZ I G, DARNAULT C J G. Aggregation and transport of nano-TiO2 in saturated porous media: Effects of pH, surfactants and flow velocity [J]. Water Research, 2011, 45(2): 839-851. doi: 10.1016/j.watres.2010.09.013 [24] 张柯柯, 汪敏浩. 离子强度和富里酸对纳米TiO2在土壤中迁移的影响 [J]. 广州化工, 2014, 42(23): 110-111. ZHANG K K, WANG M H. Effect of ionic strength and fulvic acid on the transport of TiO2 nanoparticles in soil [J]. Guangzhou Chemical Industry, 2014, 42(23): 110-111(in Chinese).
[25] 杜宜春, 梁志卿. 多孔介质中纳米颗粒及其携带重金属污染物迁移研究进展 [J]. 河南科技, 2019(32): 133-136. doi: 10.3969/j.issn.1003-5168.2019.32.047 DU Y C, LIANG Z Q. Review of migration of nanoparticle and its facilitated heavy metal contaminants in porous media [J]. Henan Science and Technology, 2019(32): 133-136(in Chinese). doi: 10.3969/j.issn.1003-5168.2019.32.047
[26] PIRSON S J. Performance of fractured oil reservoirs [J]. AAPG Bulletin, 1953, 377(2): 232-244. doi: 10.1306/5ceadc4a-16bb-11d7-8645000102c1865d [27] PHILIP J R. The theory of absorption in aggregated media [J]. Soil Research, 1968, 6(1): 1. doi: 10.1071/SR9680001 [28] FANG J, XU M J, WANG D J, et al. Modeling the transport of TiO2 nanoparticle aggregates in saturated and unsaturated granular media: Effects of ionic strength and pH [J]. Water Research, 2013, 47(3): 1399-1408. doi: 10.1016/j.watres.2012.12.005 [29] PACKMAN A I, BROOKS N H, MORGAN J J. A physicochemical model for colloid exchange between a stream and a sand streambed with bed forms [J]. Water Resources Research, 2000, 36(8): 2351-2361. doi: 10.1029/2000WR900059 [30] WANG D J, CHU L Y, PARADELO M, et al. Transport behavior of humic acid-modified nano-hydroxyapatite in saturated packed column: Effects of Cu, ionic strength, and ionic composition [J]. Journal of Colloid and Interface Science, 2011, 360(2): 398-407. doi: 10.1016/j.jcis.2011.04.064 [31] 方华祥. 溶液性质对金属纳米颗粒在多孔介质中迁移持留行为的影响机制研究[D]. 杭州: 浙江大学, 2017: 87. FANG H X. The influence mechanisms of solutions properties on the transport and retention behaviors of metal-based nanoparticles in porous media[D]. Hangzhou: Zhejiang University, 2017: 87(in Chinese).
-
点击查看大图
计量
- 文章访问数: 2449
- HTML全文浏览数: 2449
- PDF下载数: 33
- 施引文献: 0
