張智察,倪長(zhǎng)健*,鄧 也,湯津贏,朱育雷,楊寅山,鄧佩云

?

氣溶膠等效復(fù)折射率反演的免疫進(jìn)化算法

張智察1,2,倪長(zhǎng)健1,2*,鄧 也3,湯津贏1,2,朱育雷1,2,楊寅山1,2,鄧佩云1,2

(1.成都信息工程大學(xué)大氣科學(xué)學(xué)院,四川 成都 610225；2.高原大氣與環(huán)境四川省重點(diǎn)實(shí)驗(yàn)室,四川 成都 610225；3.成都市環(huán)境保護(hù)科學(xué)研究院,四川 成都 610072)

通過(guò)Mie散射理論公式構(gòu)建目標(biāo)函數(shù),利用免疫進(jìn)化算法對(duì)氣溶膠等效復(fù)折射率的實(shí)部和虛部進(jìn)行協(xié)同優(yōu)化,據(jù)此創(chuàng)新性地提出了氣溶膠等效復(fù)折射率反演的新途徑.基于成都市2017年9~12月逐時(shí)的氣溶膠散射系數(shù)和吸收系數(shù)觀測(cè)數(shù)據(jù)以及該時(shí)段同時(shí)次GRIMM180大氣顆粒物監(jiān)測(cè)儀的連續(xù)監(jiān)測(cè)資料,研究結(jié)果表明,氣溶膠等效復(fù)折射率反演的免疫進(jìn)化算法不僅是普適的,而且還具有收斂速度快?計(jì)算穩(wěn)定和求解精度高等特點(diǎn).通過(guò)與其它氣溶膠等效復(fù)折射率反演方法的對(duì)比分析,進(jìn)一步論證了新方法的優(yōu)勢(shì),這為氣溶膠等效復(fù)折射率演變機(jī)理以及氣溶膠吸濕性增長(zhǎng)模型的后續(xù)研究提供了算法保障.

氣溶膠；復(fù)折射率；反演；Mie散射理論；免疫進(jìn)化算法

作為由液態(tài)或固態(tài)微粒在空氣中組成的懸浮體系,氣溶膠不僅在天氣氣候演變過(guò)程中扮演著重要角色,也是目前大氣復(fù)合污染研究的熱點(diǎn)問(wèn)題[1-3].諸多研究表明[4-7],氣溶膠消光是大氣消光的主體,也是能見(jiàn)度降低的主要誘因,并與顆粒物質(zhì)量濃度、化學(xué)組分以及濕度條件密切相關(guān).等效復(fù)折射率是表征氣溶膠光學(xué)特性的重要參數(shù),其實(shí)部主要與光的散射有關(guān),虛部主要與光的

吸收有關(guān),由此決定了氣溶膠環(huán)境氣候效應(yīng)[8-11].目前氣溶膠等效復(fù)折射率的測(cè)量方法主要包括直接測(cè)量法[12-13]和光聲法[14],但兩種方法在實(shí)際應(yīng)用中仍面臨精度和技術(shù)設(shè)備復(fù)雜性之間的矛盾.反演法是計(jì)算復(fù)折射率的一種有效手段,它基于粒子數(shù)濃度的譜分布以及散射系數(shù)和吸收系數(shù)的測(cè)量,通過(guò)Mie散射理論對(duì)氣溶膠等效復(fù)折射率進(jìn)行求解.進(jìn)一步分析表明,反演法雖然物理意義清晰,但由于其計(jì)算對(duì)象非常復(fù)雜,傳統(tǒng)方法無(wú)法精確獲得等效復(fù)折射率的實(shí)部和虛部,常用的枚舉法[15]和人工神經(jīng)網(wǎng)絡(luò)算法[16-17]分別存在計(jì)算量大以及計(jì)算復(fù)雜和過(guò)擬合等不足.

20世紀(jì)90年代起,遺傳算法的出現(xiàn)為諸多領(lǐng)域復(fù)雜優(yōu)化問(wèn)題的解決提供了新途徑[18].在深入理解現(xiàn)有進(jìn)化算法的基礎(chǔ)上,受生物免疫機(jī)制啟發(fā),倪長(zhǎng)健等[19]提出了免疫進(jìn)化算法,多峰函數(shù)的測(cè)試結(jié)果表明,該算法具有收斂速度快,計(jì)算穩(wěn)定以及普適性強(qiáng)等優(yōu)點(diǎn).基于馬爾可夫鏈的分析指出[20],免疫進(jìn)化算法不僅是全局收斂的,而且比遺傳算法具備更快的收斂速度.近年來(lái),免疫進(jìn)化算法也獲得了較為廣泛應(yīng)用,楊懷金等[21]利用該算法對(duì)鶴望蘭葉面積指數(shù)進(jìn)行模擬,計(jì)算結(jié)果的平均相對(duì)誤差只有3.44%,取得了滿意的模擬效果;李祚泳等[22]將免疫進(jìn)化算法用于優(yōu)化地下水水質(zhì)評(píng)價(jià)冪函數(shù)加權(quán)加和型綜合指數(shù)公式中的參數(shù),得到的普適指數(shù)公式具有可比性、通用性和適用性等特點(diǎn);投影尋蹤動(dòng)態(tài)聚類模型是一種數(shù)據(jù)自驅(qū)動(dòng)的統(tǒng)計(jì)新方法,基于免疫進(jìn)化算法確定其最佳投影方向可極大地簡(jiǎn)化模型的計(jì)算量,由此提升模型的實(shí)用性[23].

本文針對(duì)目前氣溶膠等效復(fù)折射率反演算法存在的問(wèn)題,通過(guò)Mie散射理論公式構(gòu)建目標(biāo)函數(shù),利用免疫進(jìn)化算法對(duì)其中氣溶膠等效復(fù)折射率的實(shí)部和虛部進(jìn)行協(xié)同優(yōu)化,據(jù)此創(chuàng)新性地提出了氣溶膠等效復(fù)折射率反演的新途徑,并對(duì)算法性能進(jìn)行了系統(tǒng)的測(cè)試分析.

1 資料介紹

本文使用的觀測(cè)資料包括成都市2017年9~12月期間AURORA 3000型濁度計(jì)觀測(cè)的氣溶膠散射系數(shù)?AE-31型黑碳檢測(cè)儀反演的氣溶膠吸收系數(shù)以及GRIMM180大氣顆粒物監(jiān)測(cè)儀監(jiān)測(cè)的氣溶膠粒徑譜資料,具體說(shuō)明如下:

(1) AURORA 3000型濁度計(jì)觀測(cè)波長(zhǎng)為520nm,數(shù)據(jù)采集頻率為5min/次,檢測(cè)范圍>0.25,每24h進(jìn)行零點(diǎn)檢查,24h進(jìn)行零點(diǎn)漂移<±1%,每周用R134a氣體進(jìn)行跨度標(biāo)定.此外,通過(guò)內(nèi)部溫濕度傳感器來(lái)控制內(nèi)部加熱系統(tǒng),使得儀器內(nèi)部腔室中氣溶膠相對(duì)濕度控制在40%左右.由于濁度計(jì)觀測(cè)的是520nm波段的散射系數(shù)sca520,Mm-,按文獻(xiàn)[24]提出的訂正公式,得到550nm波段氣溶膠散射系數(shù)sca550,Mm-,見(jiàn)式(1).

(2) AE-31黑碳檢測(cè)儀觀測(cè)黑碳質(zhì)量濃度BC,μg/m3,數(shù)據(jù)采集頻率為5min/次.在所有氣溶膠粒子中,黑碳對(duì)光的吸收效應(yīng)最強(qiáng)[25].基于聲光法測(cè)量吸收系數(shù)的代表性已得到目前國(guó)際的普遍認(rèn)可,吳兌等[26]利用該方法得到532nm波段的氣溶膠吸收系數(shù)abs532,Mm-,構(gòu)建了其與880nm波段所測(cè)量黑碳濃度,μg/m3,的統(tǒng)計(jì)學(xué)關(guān)系,見(jiàn)式(2).對(duì)比分析表明,該統(tǒng)計(jì)學(xué)關(guān)系與其他地區(qū)的類似研究結(jié)果總體一致[27-28],并在相關(guān)研究中得到了應(yīng)用[29].再將532nm波段的吸收系數(shù)abs532,Mm-,進(jìn)一步訂正到550nm波段的吸收系數(shù)abs550,Mm-,見(jiàn)式(3).

(3) GRIMM180大氣顆粒物監(jiān)測(cè)儀經(jīng)除濕后可以實(shí)時(shí)測(cè)量大氣中PM10、PM2.5和PM1的質(zhì)量濃度和31個(gè)粒徑段的“干”氣溶膠數(shù)濃度,其中各粒徑段粒子直徑的中值分別為0.25、0.28、0.3、0.35、0.4、0.45、0.5、0.58、0.65、0.7、0.8、1.0、1.3、1.6、2.0、2.5、3.0、3.5、4.0、5.0、6.5、7.5、8.0、10.0、12.5、15.0、17.5、20.0、25.0、30.0、32.0μm,儀器數(shù)據(jù)采集頻率為1min/次.

AURORA 3000型濁度計(jì)和AE-31型黑碳檢測(cè)儀位于成都市環(huán)境保護(hù)科學(xué)研究院(104°02′E, 30°39′N)頂樓,距地面25m;GRIMM180大氣顆粒物監(jiān)測(cè)儀位于成都市一環(huán)路聯(lián)益大廈(104°02′E, 30°39′N)頂樓,距地面81m.兩處監(jiān)測(cè)點(diǎn)位周圍無(wú)明顯大氣污染源,二者直線距離為410m,環(huán)境氣象條件基本一致.

2 算法的設(shè)計(jì)

2.1 目標(biāo)函數(shù)的構(gòu)建

基于Mie散射理論[30],大氣多粒子氣溶膠散射系數(shù)和吸收系數(shù)計(jì)算公式分別表示如下:

2.2 算法流程

式(6)是一個(gè)非常復(fù)雜的非線性優(yōu)化問(wèn)題,利用傳統(tǒng)手段很難在求解精度和計(jì)算效率之間尋求平衡.為此,將免疫進(jìn)化算法應(yīng)用于氣溶膠等效復(fù)折射率的反演,為方便起見(jiàn),將氣溶膠等效復(fù)折射率的實(shí)部(re)和虛部(i)分別表示為1和2.免疫進(jìn)化算法中的子代個(gè)體生殖方式為:

免疫進(jìn)化算法把進(jìn)化建立在最優(yōu)個(gè)體的基礎(chǔ)上,通過(guò)標(biāo)準(zhǔn)差的調(diào)整把局部搜索和全局搜索有機(jī)結(jié)合起來(lái),能較好地克服現(xiàn)有進(jìn)化算法不成熟收斂的缺點(diǎn),提高算法在中后期的搜索效率.設(shè)群體規(guī)模為,復(fù)折射率反演的免疫進(jìn)化算法計(jì)算步驟如下:

(1)確定復(fù)折射率反演問(wèn)題的表達(dá)方式為:

(3)根據(jù)式(7)進(jìn)行進(jìn)化操作,在復(fù)折射率解空間內(nèi)生成子代群體.

3 結(jié)果分析

3.1 新算法的性能測(cè)試

根據(jù)目前對(duì)大氣氣溶膠等效復(fù)折射率測(cè)定的相關(guān)研究成果,確定其實(shí)部和虛部的尋優(yōu)區(qū)間分別為[1.000, 2.000]和[0.001, 0.500],據(jù)此給出算法的相關(guān)計(jì)算參數(shù),見(jiàn)表1.

表1 算法的相關(guān)計(jì)算參數(shù)

基于成都市2017年9月至2017年12月期間AURORA 3000型濁度計(jì)觀測(cè)的氣溶膠散射系數(shù)?AE-31型黑碳檢測(cè)儀反演的氣溶膠吸收系數(shù)以及GRIMM180大氣顆粒物監(jiān)測(cè)儀監(jiān)測(cè)的氣溶膠粒徑譜數(shù)據(jù),通過(guò)對(duì)三類資料進(jìn)行質(zhì)量控制,共獲得研究樣本163個(gè).對(duì)目標(biāo)函數(shù)(1,2)進(jìn)行變換,得到相應(yīng)的適應(yīng)度函數(shù)(1,2),見(jiàn)式(9).

為進(jìn)一步展示算法的計(jì)算過(guò)程,基于空氣質(zhì)量指數(shù)(AQI)選擇2017年11月10日01:00、2017年12月2日04:00、2017年12月21日09:00、2017年12月26日10:00、2017年12月23日08:00、2017年12月29日04:00 6個(gè)時(shí)刻的觀測(cè)數(shù)據(jù)作為樣本,該6個(gè)時(shí)刻涵蓋了6個(gè)等級(jí)的空氣質(zhì)量狀況,對(duì)應(yīng)的空氣質(zhì)量分別為優(yōu)?良?輕度污染?中度污染?重度污染和嚴(yán)重污染,分別繪制了6個(gè)樣本每代最優(yōu)個(gè)體適應(yīng)度F隨進(jìn)化序列的演變曲線,見(jiàn)圖1.由圖可見(jiàn),計(jì)算六個(gè)樣本全局最優(yōu)解所需的總進(jìn)化序列分別為7、19、24、16、19和13,對(duì)應(yīng)的相對(duì)誤差(12)分別為0.9%、0.3%、0.7%、0.4%、0.8%和0.6%.上述測(cè)試結(jié)果表明,氣溶膠等效復(fù)折射率反演的免疫進(jìn)化算法能夠穩(wěn)定地得到全局最優(yōu)解,并具有收斂速度快和計(jì)算精度高等特點(diǎn).

圖1 每代最優(yōu)個(gè)體的適應(yīng)度F隨進(jìn)化序列的演變曲線

Fig.1 Variation of the best fitnessFof each generation with evolution sequence

a: 2017年11月10日01:00; b: 2017年12月2日04:00; c: 2017年12月21日09:00; d: 2017年12月26日10:00; e: 2017年12月23日08:00; f: 2017年12月29日04:00

3.2 不同反演算法的對(duì)比分析

為對(duì)比起見(jiàn),表2同時(shí)也給出了人工神經(jīng)網(wǎng)絡(luò)BP算法的相對(duì)誤差.結(jié)果發(fā)現(xiàn),(1)氣溶膠等效復(fù)折射率反演的免疫進(jìn)化算法不僅計(jì)算結(jié)果優(yōu)于人工神經(jīng)網(wǎng)絡(luò)BP算法,而且不同樣本的相對(duì)誤差也比較均勻,表明該算法具有更好的魯棒性;(2)基于免疫進(jìn)化算法求解3個(gè)樣本全局最優(yōu)解所需的總進(jìn)化序列平均為48代,相對(duì)3.1在計(jì)算精度提高10倍的前提下,計(jì)算量只增加了1.4倍左右.這表明該算法具有良好的尋優(yōu)能力,能在計(jì)算精度和計(jì)算效率之間保持很好的平衡;(3)作為一種黑箱模型,氣溶膠等效復(fù)折射率反演的人工神經(jīng)網(wǎng)絡(luò)BP算法受訓(xùn)練樣本完備性影響大,同時(shí)也未能真正表征Mie散射理論的非線性映射關(guān)系.因此,免疫進(jìn)化算法求解氣溶膠等效復(fù)折射率具有普適性.另外,假定實(shí)部和虛部小數(shù)點(diǎn)后均保留3位有效數(shù)字,在上述尋優(yōu)區(qū)間內(nèi)枚舉法的計(jì)算量及其龐大,每一復(fù)折射率的求解需要進(jìn)行1.0×106次循環(huán).基于表1給出的計(jì)算參數(shù),利用氣溶膠等效復(fù)折射率反演的免疫進(jìn)化算法對(duì)每一復(fù)折射率的求解平均只需進(jìn)行2.0×103次循環(huán),其計(jì)算量約為枚舉法的1/500.因此,氣溶膠等效復(fù)折射率反演的免疫進(jìn)化算法具有極優(yōu)的計(jì)算效率.

表2 兩種復(fù)折射率反演方法的對(duì)比

4 結(jié)論

4.1 通過(guò)Mie散射理論公式構(gòu)建目標(biāo)函數(shù),利用免疫進(jìn)化算法對(duì)其中氣溶膠等效復(fù)折射率的實(shí)部和虛部進(jìn)行協(xié)同優(yōu)化,創(chuàng)新性地提出了氣溶膠等效復(fù)折射率反演的新途徑.

4.2 基于成都市2017年9~12月逐時(shí)的氣溶膠散射系數(shù)和吸收系數(shù)數(shù)據(jù),結(jié)合該時(shí)段同時(shí)次GRIMM180大氣顆粒物監(jiān)測(cè)儀的連續(xù)監(jiān)測(cè)資料,計(jì)算結(jié)果表明,氣溶膠等效復(fù)折射率反演的免疫進(jìn)化算法不僅是普適的,而且還具有收斂速度快、計(jì)算穩(wěn)定和求解精度高等特點(diǎn).

4.3 氣溶膠等效復(fù)折射率反演的免疫進(jìn)化算法系統(tǒng)克服了目前枚舉法和人工神經(jīng)網(wǎng)絡(luò)方法反演氣溶膠等效復(fù)折射率存在的不足,并為氣溶膠等效復(fù)折射率演變機(jī)理以及氣溶膠吸濕性增長(zhǎng)模型的后續(xù)研究提供了算法保障.

[1] Chylek P, Coakley J A. Aerosols and Climate [J]. Science, 1974,183(4120):75-77.

[2] Toon O B, Pollack J B. Atmospheric Aerosols and Climate: Small particles in the Earth's atmosphere interact with visible and infrared light, altering the radiation balance and the climate [J]. American Scientist, 1980,68(3):268-278.

[3] Toon O B, Pollack J B. Atmospheric Aerosols and Climate: Small particles in the Earth's atmosphere interact with visible and infrared light, altering the radiation balance and the climate [J]. American Scientist, 1980,68(3):268-278.

[4] 王 英,李令軍,李成才.北京大氣能見(jiàn)度和消光特性變化規(guī)律及影響因素 [J]. 中國(guó)環(huán)境科學(xué), 2015,35(5):1310-1318.Wang Ying, Li Ling-jun, Li Cheng-cai. The variation characteristic and influence factors of atmospheric visibility and extinction effect in Beijing [J]. China Environmental Science, 2015,35(5):1310-1318.

[5] 賈 佳,韓力慧,程水源,等.京津冀區(qū)域PM2.5及二次無(wú)機(jī)組分污染特征研究[J]. 中國(guó)環(huán)境科學(xué), 2018,38(3):801-811.Jia Jia, Han Lihui, Cheng Shuiyuan, et al. Pollution characteristic of PM2.5and secondary inorganic ions in Beijing- Tianjin-Hebei region [J]. China Environmental Science, 2018,38(3): 801-811.

[6] 浦靜姣,徐宏輝,馬千里.長(zhǎng)江三角洲背景地區(qū)大氣污染對(duì)能見(jiàn)度的影響[J]. 中國(guó)環(huán)境科學(xué), 2017,37(12):4435-4441. Pu Jingjiao, Xu Honghui, Ma Qianli. Impacts of atmospheric pollution on visibility in the background area of Yangtze River Delta [J]. China Environmental Science, 2017,37(12):4435-4441.

[7] Ouimette J R, Flagan R C. The extinction coefficient of multicomponent aerosols [J]. Atmospheric Environment, 1982,16(10): 2405-2419.

[8] 張學(xué)海,戴聰明,武鵬飛,等.折射率和粒子尺度對(duì)大氣氣溶膠光散射特性的影響 [J]. 紅外與激光工程, 2017,001(12):104-110.Zhang Xuehai, Dai Congming, Wu Pengfei, et al. Effect of the refractive index and particle size parameter on light scattering properties of atmosphere aerosol [J]. Infrared and Laser Engineering, 2017,001(12):104-110.

[9] 孫 剛,翁寧泉,肖黎明,等.大氣溫度分布特性及對(duì)折射率結(jié)構(gòu)常數(shù)的影響 [J]. 光學(xué)學(xué)報(bào), 2004,24(5):592-596. Sun Gang, Wen Ningquan, Xiao Liming, et al., Profile and character of atmospheric temperature [J]. Acta Optical Sinica, 2004,24(5):592- 596.

[10] Colarco P R, Kim K, Nowottnick E P, et al. Importance of Particle Shape and Refractive Index in Climate Simulations of the Dust Aerosol Lifecycle [J]. Anticancer Research, 2002,22(5):2805-2809.

[11] Banakh V A, Mironov V L, Smalikho I N. Light reflection in a turbulent atmosphere in the presence of induced temperature inhomogeneity of the refractive index [J]. Journal of Pediatrics, 2009,154(4):588-591.

[12] Lin C I, Baker M, Charlson R J. Absorption coefficient of atmospheric aerosol: a method for measurement [J]. Applied Optics, 1973,12(6): 1356-63.

[13] Wolff G, Stroup C, Stroup D. The coefficient of haze as a measure of particulate elemental carbon [J]. AirRepair, 1983,33(8):746-750.

[14] Adams K M. Real-time in situ measurements of atmospheric optical absorption in the visible via photoacoustic spectroscopy. 1: Evaluation of photoacoustic cells [J]. Journal of Plant Nutrition & Soil Science, 2003,166(6):4052-6.

[15] 李學(xué)彬,黃印博,徐赤東,等.廈門地區(qū)氣溶膠折射率的測(cè)量 [J]. 光學(xué)精密工程, 2008,16(10):1831-1835. Li Xuebin, Huang Yinbo, Xu Chidong, et al.Measurement of refractive index for aerosol particle in Xiamen [J]. Optical and Precision Engineering, 2008,16(10):1831-1835.

[16] 亓秀永,黃印博,方 帥,等.氣溶膠粒子等效復(fù)折射率的BP神經(jīng)網(wǎng)絡(luò)建模 [J]. 大氣與環(huán)境光學(xué)學(xué)報(bào), 2009,4(1):30-37. Qi Xiuyong, Huang Yinbo, Fang Shuai. BP neural network model of the equivalent complex refractive index of aerosol [J]. Journal of Atmospheric and Environmental Optics, 2009,4(1):30-37.

[17] 韓道文,劉文清,陸亦懷,等.基于Madaline網(wǎng)絡(luò)的氣溶膠散射系數(shù)的反演算法 [J]. 光學(xué)學(xué)報(bào), 2007,27(3):384-390.Han Daowen, Liu Wenqing, Lu Yihuai, et al. A retrieve method for aerosol extinction coefficient based on madaline network [J]. Acta OpticaSinica, 2007,27(3):384-390.

[18] 葛繼科,邱玉輝,吳春明,等.遺傳算法研究綜述[J]. 計(jì)算機(jī)應(yīng)用研究, 2008,25(10):2911-2916. Ge Jike, Qiu Yuhui, Wu ChunMing, et al. Summary of genetic algorithms research [J]. Application Research of Computers, 2008, 25(10):2911-2916.

[19] 倪長(zhǎng)健,丁 晶,李祚泳.免疫進(jìn)化算法 [J]. 西南交通大學(xué)學(xué)報(bào), 2003,38(1):87-91.Ni Changjian, Ding Jing, Li Zuoyong. Immune evolutionary algorithm [J]. Journal of Southwest Jiaotong University, 2003,38(1): 87-91.

[20] 倪長(zhǎng)健,丁 晶,李祚泳.基于優(yōu)秀抗體的免疫算法及其收斂性問(wèn)題的研究 [J]. 系統(tǒng)工程, 2002,20(3):72-76.Ni Changjian, Ding Jing, Li Zuoyong. Study on immune algorithm based on superior antibodies and its convergence property [J]. Systems Engineering, 2002,20(3):72-76.

[21] 楊懷金,葉芝祥,朱克云,等.基于免疫進(jìn)化算法的鶴望蘭葉面積指數(shù)模擬[J]. 生態(tài)學(xué)報(bào), 2006,26(8):2744-2748. Yang Huaijin, Ye Zhixiang, Zhu Keyun, et al. The simulation of strelitziareginaelai (Leaf Area Index) based on the optimization of immune evolutionary algorithms (IEA) [J]. Acta. Ecologica. Sinica, 2006,26(8):2744-2748.

[22] 李祚泳,汪嘉楊,程會(huì)珍.基于免疫進(jìn)化算法優(yōu)化的地下水水質(zhì)評(píng)價(jià)普適公式[J]. 水科學(xué)進(jìn)展, 2008,19(5):707-713. Li Zuoyong, Wang Jiayang, Cheng Huizhen. Universal exponential formula for underground water assessment base on immune evolutionary algorithm [J]. Advances in Water Science, 2008,19(5): 707-713.

[23] 王順久,倪長(zhǎng)健.投影尋蹤動(dòng)態(tài)聚類模型及其應(yīng)用[J]. 哈爾濱工業(yè)大學(xué)學(xué)報(bào), 2009,(1):178-180. Wang Shunjiu, Ni Changjian. A projection pursuit dynamic cluster model and its application [J]. Journal of Harbin Insititute of Technology, 2009,(1):178-180.

[24] 劉新罡,張遠(yuǎn)航,曾立民,等.廣州市大氣能見(jiàn)度影響因子的貢獻(xiàn)研究 [J]. 氣候與環(huán)境研究, 2006,11(6):733-738. Liu Xingang, Zhang Yuanhang, Zeng Limin, et al. Research on the influential factor of visibility reduction in Guangzhou [J]. Climate and Environmental research, 2006,11(6):733-738.

[25] Japar S M, Brachaczek W W, Jr R A G, et al. The contribution of elemental carbon to the optical properties of rural atmospheric aerosols [J]. Atmospheric Environment, 1986,20(6):1281-1289.

[26] 吳 兌,毛節(jié)泰,鄧雪嬌,等.珠江三角洲黑碳?xì)馊苣z及其輻射特性的觀測(cè)研究 [J]. 中國(guó)科學(xué):地球科學(xué), 2009,(11):1542-1553.Wu Dui, Mao Jietai, Deng Xuejiao, et al. Black carbon aerosols and their radiative properties in the pearl river Delta region [J]. ScientiaSinica Terrae, 2009,(11):1542-1553.

[27] Arnott W P, Moosmüller H, Sheridan P J, et al. Photoacoustic and filter-based ambient aerosol light absorption measurements: Instrument comparisons and the role of relative humidity [J]. Journal of Geophysical Research Atmospheres, 2003,108(D1):AAC-1-AAC 15-11.

[28] Barnard J C, Kassianov E I, Ackerman T P, et al. Measurements of Black Carbon Specific Absorption in the Mexico City Metropolitan Area during the MCMA 2003 Field Campaign [J]. Atmospheric Chemistry & Physics Discussions, 2005,5(3):4083-4113.

[29] 姚 青,韓素芹,蔡子穎,等.天津城區(qū)春季大氣氣溶膠消光特性研究 [J]. 中國(guó)環(huán)境科學(xué), 2012,32(5):795-802. Yao Qing, Han Suqin, Cai Ziying, et al. Study on characteristic of aerosol extinction at Tianjin City in the spring [J]. China Environmental Science, 2012,32(5):795-802.

[30] 沈建琪,劉 蕾.經(jīng)典Mie散射的數(shù)值計(jì)算方法改進(jìn) [J]. 中國(guó)粉體技術(shù), 2005,11(s1):45-50.Shen Jianqi, Liu Lei. An improved algorithm of classical Mie scattering calculation [J]. China Powder Science and Technology, 2005,11(s1):45-50.

致謝：本論文的Mie散射理論算法受上海理工大學(xué)上海-漢堡國(guó)際工程學(xué)院沈建琪院長(zhǎng)提出的Mie散射理論改進(jìn)算法啟發(fā),并在沈建琪院長(zhǎng)鼓勵(lì)與支持下成功對(duì)其進(jìn)行了改進(jìn),在此對(duì)沈建琪院長(zhǎng)表示衷心感謝.

Retrieval of equivalent complex refractive index of aerosol particles based on immune evolution algorithm.

ZHANG Zhi-cha1,2, NI Chang-jian1,2*, DENG Ye3, TANG Jin-ying1,2, ZHU Yu-lei1,2, YANG Yin-shan1,2, DENG Pei-yun1,2

(1.College of Atmospheric Science, Chengdu University of Information Technology, Chengdu 610225, China；2.Plateau Atmosphere and Environment Key Laboratory of Sichuan Province, Chengdu 610225, China；3.Chengdu Academy of Environmental Sciences, Chengdu 610072, China)., 2019,39(2)：554~559

Based on immune evolutionary algorithm to optimize the objective function derived from Mie scattering theory, a new approach, which aimed at retrieving equivalent complex refractive index of aerosol particles (including both real part and imaginary part as a whole), was innovatively proposed. By utilizing the hourly aerosol scattering coefficient, aerosol absorption coefficient, and the coincidental data obtained from GRIMM180 in Chengdu from September to December 2017, the results showed thatthe new approach featured universal, fast convergent, robust, and precise. Compared with other available methods of retrieving equivalent complex refractive index of aerosol particles, the advantages of the new approach were further confirmed. The above achievement indicated that thenew algorithm should help to more clearly understand the evolution mechanism of equivalent complex refractive index of aerosol particles and aerosol hygroscopic growth model.

aerosol；complex refractive index；retrieval；Mie scattering theory；immune evolutionary algorithm

X513

A

1000-6923(2019)02-0554-06

張智察(1995-),男,浙江麗水人,成都信息工程大學(xué)碩士研究生,主要從事大氣物理學(xué)與大氣環(huán)境方面研究.

2018-07-13

國(guó)家重點(diǎn)研發(fā)計(jì)劃(2018YFC0214004);四川省科技廳重點(diǎn)研發(fā)項(xiàng)目(2018SZ0287)

* 責(zé)任作者, 教授, ncj1970@163.com