双颗粒凝聚气溶胶的吸湿增长模型及散射特性 下载: 844次
1 引言
大气环境中气溶胶颗粒物的日益增多,对人类及其生存的环境产生了诸多影响[1-3],如:能见度下降[4-6],引起呼吸道和肺部疾病[7-9],造成酸雨[10-12]和气候变化[13-14]等。气溶胶颗粒的吸湿性是影响其物理性质的重要因素,随着环境相对湿度(RH)的增加,空气中的水汽会被颗粒吸收或者包覆在颗粒表面,不仅使颗粒的形貌、粒径、折射率及粒度分布等物理性质发生改变[15-16],而且吸湿颗粒表面黏性的增加使得粒子间的碰并率显著提高,形成团聚粒子,进而影响粒子的散射特性[17]。因此,研究相对湿度对凝聚气溶胶散射特性的影响显得尤为重要,该项研究能够为大气环境监测和预测环境变化提供重要的理论依据。
按照气溶胶的吸湿性强弱,可将其分为亲水性和疏水性两类。亲水性气溶胶包括无机盐气溶胶(包括硫酸盐、硝酸盐、铵盐和海盐等)和吸湿性有机气溶胶;疏水性气溶胶包括黑碳、非吸湿性有机气溶胶和沙尘等[17-19]。由于无机盐气溶胶含量较多[20],无机盐气溶胶颗粒之间或者无机盐颗粒与其他气溶胶颗粒之间易形成凝聚粒子。众所周知,在吸湿过程中,由两种无机盐组成的亲水性凝聚颗粒通常会发生逐步相变,部分颗粒在相互潮解相对湿度(MDRH)时完全潮解,溶解部分的颗粒会随着RH的进一步增加而继续吸收水分,并且当RH达到对应的二次潮解相对湿度(SDRH)时,残余的固体颗粒会完全潮解,潮解相对湿度值仅取决于颗粒的组成[21]。Gupta等[22]利用电子扫描显微镜(SEM)研究了摩尔分数在0.1~0.9时氯化钠(NaCl)-硝酸钠(NaNO3)混合颗粒在3%~93% RH范围内的吸湿性能,发现在NaCl的摩尔分数大于0.38的情况下,NaCl摩尔分数的增加使得二次潮解相对湿度值趋向纯NaCl的潮解相对湿度,而富含NaNO3的情况则对应相反。上述实验规律在Li等[23]对NaCl-氯化钾(KCl)混合物吸湿性的研究中以及Gupta等[24]对NaCl-氯化镁(MgCl2)混合物吸湿性的研究中均有体现。Shi等[25]利用吸湿串联差动分析仪研究了在10%~90% RH范围内苯甲酸和无机盐(氯化钠和硫酸铵)组成的内混合颗粒吸湿行为,发现无机盐颗粒的潮解抑制了苯甲酸在较高RH下的进一步吸湿性生长。
综上所述,目前人们对于亲水性气溶胶之间凝聚粒子的研究处于对其潮解相对湿度的测量以及吸湿现象的观察阶段,尚未对这一体系中RH所引起的散射特性变化规律进行研究;对于亲疏水性凝聚的气溶胶粒子研究则处于基础的观测阶段,并未对RH所引起的散射特性规律进行详细研究。因此,本文将分别建立亲水性和亲疏水性双颗粒凝聚气溶胶的吸湿模型,并对它们的散射特性进行研究,这为进一步研究多颗粒凝聚粒子的吸湿散射特性变化提供一定的理论基础,将有助于研究人员更好地理解气溶胶吸湿所造成的物理现象,为基于光散射技术的大气颗粒物探测研究提供一定的参考。
2 理论模型
2.1 双颗粒凝聚的亲水性气溶胶吸湿模型
考虑到高湿度的大气环境中凝聚粒子存在的普遍性[26],且亲水性颗粒含量较多,下面将在单颗粒亲水性气溶胶吸湿模型的基础上,建立双颗粒凝聚的亲水性气溶胶吸湿模型,并讨论其散射特性随RH的变化规律。由Tang等[27]使用光学粒子计数器对不同湿度下凝聚气溶胶粒径分布的研究可知:两种不同成分的亲水性颗粒凝聚成一个粒子,随着RH的增加,凝聚粒子会产生两次潮解。因此,本文将亲水性凝聚粒子的吸湿过程分为一次潮解前、一次潮解后和二次潮解后三个阶段,并对这三个阶段分别进行建模,双颗粒凝聚的亲水性气溶胶吸湿模型如
式中:RRH为RH;k为吸湿性参数[29];Vdry为亲水性干粒子的体积。
图 1. 双颗粒凝聚的亲水性气溶胶的吸湿模型示意图。(a)干凝聚粒子;(b)一次潮解前的湿凝聚粒子;(c)一次潮解时的湿凝聚粒子;(d)二次潮解时的湿凝聚粒子
Fig. 1. Diagram of moisture absorption model of two-particle agglomerated hydrophilic aerosol. (a) Dry agglomerated aerosol; (b) wet agglomerated aerosol before first deliquescence; (c) wet agglomerated aerosol for first deliquescence; (d) wet agglomerated aerosol for second deliquescence
当环境湿度达到亲水性凝聚粒子的MDRH时,凝聚的亲水性粒子发生第一次潮解,一部分固态粒子与水溶解形成液滴,剩余的未溶解固相粒子仍与水层组成核壳模型,等待二次潮解,如
式中:Ddry为干粒子的直径;A为常量,A=
当RH达到亲水性凝聚粒子的二次潮解相对湿度时,剩余粒子发生潮解,并与已潮解部分凝聚为一颗更大的液滴,如
2.2 双颗粒凝聚的亲疏水性气溶胶吸湿模型
由于大气环境中的亲水性颗粒含量多,除了亲水性颗粒之间会凝聚成一个大颗粒外,亲水性颗粒还会与疏水性粒子发生碰并,形成亲疏水性的凝聚粒子。下面将建立双颗粒的亲疏水性凝聚粒子吸湿模型。考虑到疏水性粒子的吸湿性弱,不会发生潮解,而亲水性颗粒具有较强的吸湿性,在高湿度条件下有潮解特性,于是将亲疏水性凝聚粒子的吸湿过程分为潮解前和潮解后两个过程,如
图 2. 双颗粒凝聚的亲疏水性气溶胶吸湿模型示意图。(a)干凝聚粒子;(b)潮解前的湿凝聚粒子;(c)潮解后的湿凝聚粒子
Fig. 2. Diagrams of moisture absorption models of agglomerated hydrophilic particle and hydrophobic particle. (a) Dry aggregated aerosol; (b) wet aggregated aerosol before deliquescence; (c) wet aggregated aerosol after deliquescence
在潮解前,亲疏水性凝聚粒子的吸湿过程如
式中:m表示疏水性粒子表面的水层厚度;RB是疏水性干粒子的半径。
当大气环境RH达到亲疏水性凝聚粒子中亲水性粒子的DRH时,亲水性粒子潮解形成饱和液滴,疏水性粒子与已潮解的液态粒子形成
需要说明的是,对于亲疏水性的双颗粒凝聚粒子,因疏水性粒子不溶于液态粒子,潮解后两者的位置关系并不确定,可能存在部分吞没状态,也可能形成内混合状态[26-28,32]。因此,本文将通过改变疏水性粒子在内混合体中的位置,研究内部结构对亲疏水性凝聚粒子散射特性的影响,潮解后亲疏水性凝聚粒子的位置关系示意图如
式中:hB为q1的球冠高度,hB=
而Vx又可以表示为
比较(4)式和(5)式,得到了关于Rx的方程为
图 3. 潮解后的亲疏凝聚粒子的相对位置示意图。(a)外切;(b)部分吞没;(c)内切;(d)同心
Fig. 3. Diagrams of relative positions of agglomerated hydrophilic particle and hydrophobic particle after deliquescence. (a) External circumscribing; (b) partial engulfing; (c) internal circumscribing; (d) concentricity
3 散射特性的计算与讨论
3.1 计算方法
离散偶极子近似(DDA)方法允许任意形状的模型运算,被广泛运用到了电磁散射问题的求解中。该算法的核心思想是将散射体用一群离散的偶极子替代,通过求解偶极子的极化强度,最终由散射体上所有点在远场的辐射总和构成散射场,假设rj(j表示第j个电偶极子,j=1,2,3,…,N,N表示电偶极子的总数)处的电场为Ej,入射电场为
式中:Pk为第k个电偶极子的偶极矩;Ajk为一个3×3的矩阵,若定义Ajj=
由(8)式可以求出偶极矩Pk,再求解出散射截面为
从而获得散射系数为
式中:Csca为散射体的散射截面;Cext为消光截面;Cabs为吸收截面;Ω为散射立体角;n为散射方向的单位矢量;rj为第j个电偶极子的位置矢量大小;aeff为等效半径。
DDA方法最大的优势在于允许用户自定义任何几何形状的颗粒模型,以计算出颗粒的散射、吸收和消光系数、相位滞后效率因子及Mueller 散射强度矩阵等物理量。计算时,将颗粒用一群离散偶极子表示,设置的离散偶极子数目越多,计算结果越精确。本文在计算双颗粒凝聚粒子的散射特性时,均采用了50000以上的偶极子数,误差值设定为1.0×10-5。下面将采用DDA方法详细计算与讨论凝聚粒子散射系数随RH的变化规律。
3.2 双颗粒凝聚的亲水性粒子吸湿散射特性
根据双颗粒凝聚的亲水性气溶胶吸湿模型,采用DDA算法,以NaCl和NaNO3气溶胶不同体积比形成的凝聚粒子为例,详细计算了凝聚粒子散射系数Qsca随RH的变化关系。在计算过程中,干燥状态时,NaCl和NaNO3气溶胶的折射率分别为1.5480和1.4497,它们的总体积始终保持不变,为5.236×10-3 μm3,当体积比为1∶1时,两者的直径均为1 μm。潮解前,亲水性粒子的表面包裹一层水膜,此时分别采用干粒子折射率和水的折射率计算干粒子的体积与水层区域的体积;根据文献[
33],潮解后粒子的折射率是利用与湿度相关的溶解度和粒子折射率计算得到,两种亲水性粒子的折射率nNaCl(RRH)、
计算过程中,首先通过气溶胶无机-有机混合物官能团活性系数模型[34]计算出NaCl和NaNO3气溶胶凝聚粒子的一次和二次潮解相对湿度,计算结果如
当入射光的波长为0.55 μm时,利用DDA算法计算了不同体积比时,NaCl-NaNO3双颗粒凝聚粒子散射系数Qsca与RH的关系,如
上述结果表明,双颗粒凝聚的亲水性气溶胶散射系数Qsca随RH的变化规律与两颗单粒子的体积比密切相关,体积比不同,Qsca产生跃变时的相对湿度值及其增长倍数存在较大差异,特别是在亲水性凝聚粒子两次潮解的区间内时,该现象更明显。
图 5. NaCl的体积分数不同时,RH对NaCl-NaNO3亲水性凝聚粒子散射系数的影响
Fig. 5. Effects of RH on scattering coefficient of NaCl-NaNO3 hydrophilic aggregated particles under different NaCl volume fraction
3.3 双颗粒凝聚的亲疏水性粒子吸湿散射特性
利用上文建立的亲疏水性凝聚粒子的吸湿模型,根据DDA算法,详细计算亲疏水性双颗粒凝聚气溶胶散射系数Qsca与RH的变化规律。以氯化钠和烟尘凝聚粒子体系(NaCl-Soot)为例,在入射光波长为0.55 μm的干燥状态下,两粒子的总体积始终不变,为5.236×10-3 μm3,两者的折射率取1.548和1.75-0.04i。
图 6. 在不同的相对位置条件下,NaCl的体积分数和RH对NaCl-Soot亲疏水性凝聚粒子散射系数的影响。(a)外切;(b)部分吞没;(c)内切;(d)同心
Fig. 6. Effects of NaCl volume fraction and RH on scattering coefficient of NaCl-Soot hydrophobic agglomerated aerosols under different relative position conditions. (a) External circumscribing; (b) partial engulfing; (c) internal circumscribing; (d) concentricity
综上所述,在入射光波长和凝聚粒子体积不变的条件下,不同成分颗粒的体积比是影响亲水性和亲疏水性双颗粒凝聚气溶胶散射系数与RH的变化规律的主要因素。
4 结论
运用热力学原理与表面吸附理论,建立了两种双颗粒凝聚气溶胶的吸湿模型。采用离散偶极子近似方法,计算了不同相对湿度时两种气溶胶凝聚粒子的散射特性。结果表明:在40%~90%相对湿度范围内,当入射光波长为0.55 μm,两颗干粒子的总体积不变时,不同体积比的亲水性双颗粒凝聚气溶胶 (以NaCl-NaNO3颗粒为例)的散射系数存在两次跃变,分别发生在两次潮解相对湿度的位置处,而散射系数跃变的位置和增幅与NaCl-NaNO3颗粒的体积比密切相关;二次潮解后,不同体积比亲水性凝聚粒子的散射系数随相对湿度的增加均呈现出指数增长趋势,且NaNO3的体积分数越大,散射系数的增长幅度越大。对于亲疏水性双颗粒凝聚气溶胶 (以NaCl-烟尘颗粒为例),不同体积比凝聚粒子的散射系数随相对湿度的增加均呈现出指数增长趋势,NaCl的体积分数越大,散射系数的增长速度越快,即NaCl和Soot颗粒体积比也是影响该凝聚粒子散射特性的主要因素;此外,结果还表明潮解后亲疏水性两粒子的相对位置关系对凝聚粒子散射系数的影响较小。上述结果能够为气溶胶浓度和大气能见度的精确测量、正确认识大气能见度衰减规律及全球气候变化提供重要的理论基础。
致谢 特别感谢美国普林斯顿大学的Dranie教授提供DDA计算程序。
[1] Jia R, Luo M, Liu Y Z, et al. Anthropogenic aerosol pollution over the eastern slope of the Tibetan Plateau[J]. Advances in Atmospheric Sciences, 2019, 36(8): 847-862.
[2] Feng Y, Chen D M, Zhang X H. Atmospheric aerosol pollution across China: a spatiotemporal analysis of satellite-based aerosol optical depth during 2000—2016[J]. International Journal of Digital Earth, 2019, 12(7): 843-857.
[3] Yu M M, Yuan X L, He Q Q, et al. Temporal-spatial analysis of crop residue burning in China and its impact on aerosol pollution[J]. Environmental Pollution, 2019, 245: 616-626.
[5] Fei Y, Fu D S, Song Z J, et al. Spatiotemporal variability of surface extinction coefficient based on two-year hourly visibility data in mainland China[J]. Atmospheric Pollution Research, 2019, 10(6): 1944-1952.
[11] Craig R L, Nandy L, Axson J L, et al. Spectroscopic determination of aerosol pH from acid-base equilibria in inorganic, organic, and mixed systems[J]. The Journal of Physical Chemistry A, 2017, 121(30): 5690-5699.
[12] Shi G, Xu J, Peng X, et al. pH of aerosols in a polluted atmosphere: source contributions to highly acidic aerosol[J]. Environmental Science & Technology, 2017, 51(8): 4289-4296.
[13] Gustafsson Ö, Ramanathan V. Convergence on climate warming by black carbon aerosols[J]. Proceedings of the National Academy of Sciences of the United States of America, 2016, 113(16): 4243-4245.
[14] Najafi M R, Zwiers F W, Gillett N P. Attribution of arctic temperature change to greenhouse-gas and aerosol influences[J]. Nature Climate Change, 2015, 5(3): 246-249.
[15] 张辉, 陶宗明, 麻晓敏, 等. 合肥地区PM2.5质量浓度与气溶胶后向散射系数之间吸湿因子的拟合[J]. 中国激光, 2018, 45(7): 0704006.
[16] 张学海, 戴聪明, 张鑫, 等. 相对湿度和粒子形态对海盐气溶胶粒子散射特性的影响[J]. 红外与激光工程, 2019, 48(8): 253-260.
[17] 殷金英, 郑宇浓, 杨洪艳, 等. 大气气溶胶相对湿度影响碳黑团聚物辐射特性分析[J]. 光学学报, 2017, 37(5): 0529001.
[18] Textor C, Schulz M, Guibert S, et al. Analysis and quantification of the diversities of aerosol life cycles within AeroCom[J]. Atmospheric Chemistry and Physics, 2006, 6(7): 1777-1813.
[19] Zhu Q, He L Y, Huang X F, et al. Atmospheric aerosol compositions and sources at two national background sites in northern and Southern China[J]. Atmospheric Chemistry and Physics, 2016, 16(15): 10283-10297.
[20] 张智察, 倪长健, 汤津赢, 等. “干”气溶胶等效复折射率与其质量浓度指标的相关性研究[J]. 光学学报, 2019, 39(5): 0501002.
[25] Shi Y J, Ge M F, Wang W G. Hygroscopicity of internally mixed aerosol particles containing benzoic acid and inorganic salts[J]. Atmospheric Environment, 2012, 60: 9-17.
[28] Petters M D, Kreidenweis S M. A single parameter representation of hygroscopic growth and cloud condensation nucleus activity[J]. Atmospheric Chemistry and Physics, 2007, 7(8): 1961-1971.
[33] William MH. CRC handbook of chemistry and physics[M]. 95th ed. Boca Raton: Taylor & Francis Group, 2014: 123- 133.
[35] 赵佳佳, 顾芳, 张加宏, 等. 单颗粒气溶胶的吸湿增长模型及散射特性研究[J]. 光学学报, 2020, 40(5): 0501001.
Article Outline
顾芳, 张加宏, 陈云云, 赵佳佳, 崔芬萍, 李敏, 赵静. 双颗粒凝聚气溶胶的吸湿增长模型及散射特性[J]. 光学学报, 2021, 41(3): 0301001. Fang Gu, Jiahong Zhang, Yunyun Chen, Jiajia Zhao, Fenping Cui, Min Li, Jing Zhao. Hygroscopic Growth Model and Scattering Characteristics of Two-Particle Agglomerated Aerosol[J]. Acta Optica Sinica, 2021, 41(3): 0301001.