• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    耦合氣溶膠雙參數(shù)化方案的大氣能見(jiàn)度數(shù)值改進(jìn)算法

    2021-12-01 09:50:04張智察倪長(zhǎng)健趙軍平謝滟馨王楊菲
    中國(guó)環(huán)境科學(xué) 2021年11期
    關(guān)鍵詞:顆粒物大氣模型

    張智察,倪長(zhǎng)健,趙軍平,馮 淼,謝滟馨,王楊菲

    耦合氣溶膠雙參數(shù)化方案的大氣能見(jiàn)度數(shù)值改進(jìn)算法

    張智察1,2,倪長(zhǎng)健1*,趙軍平2,馮 淼3,謝滟馨1,王楊菲1

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

    將氣溶膠復(fù)折射率(Aerosol Complex Refractive Index, ACRI)和氣溶膠粒徑吸濕增長(zhǎng)因子(Growth Factor,(RH))參數(shù)化方案進(jìn)行耦合,提出了一種基于Mie散射模型的大氣能見(jiàn)度數(shù)值改進(jìn)算法.并利用成都市2017年10~12月WS600一體式氣象站、AURORA-3000積分濁度計(jì)、AE-31黑碳儀以及GRIMM180環(huán)境顆粒物監(jiān)測(cè)儀分別觀測(cè)獲得的相對(duì)濕度(RH),干氣溶膠散射系數(shù)(sp),干氣溶膠吸收系數(shù)(sp),氣溶膠質(zhì)量濃度(PM10, PM2.5, PM1)及其數(shù)濃度粒徑分布([(RH)])的地面逐時(shí)觀測(cè)資料,通過(guò)與兩種能見(jiàn)度計(jì)算模型(經(jīng)驗(yàn)參數(shù)的Mie散射模型和統(tǒng)計(jì)模型)在不同能見(jiàn)度區(qū)間(<2km, 2~5km, 5~10km, >10km)模擬結(jié)果的對(duì)比分析,評(píng)估了該改進(jìn)算法的適用性.結(jié)果表明:三種能見(jiàn)度計(jì)算方法均能較好地模擬出能見(jiàn)度的變化特征;改進(jìn)算法通過(guò)本地化參數(shù)化方案更準(zhǔn)確地估計(jì)出DACRI和(RH),從而可更準(zhǔn)確地模擬出四類能見(jiàn)度區(qū)間,對(duì)應(yīng)模擬值與實(shí)測(cè)值的相關(guān)系數(shù)()分別為0.62, 0.90, 0.89, 0.93,平均相對(duì)誤差(MRE)分別為9.86%, 10.39%, 9.94%, 14.06%.

    氣溶膠;復(fù)折射率;吸濕增長(zhǎng)因子;參數(shù)化方案;能見(jiàn)度

    隨著我國(guó)社會(huì)經(jīng)濟(jì)的迅猛發(fā)展以及城市化進(jìn)程城的不斷加快,以大氣能見(jiàn)度降低為主要特征的霾在秋冬季頻繁發(fā)生[1].能見(jiàn)度降低不僅表征環(huán)境空氣質(zhì)量的惡化,也對(duì)人體健康和交通運(yùn)輸安全產(chǎn)生嚴(yán)重影響[2].

    Koschmieder于1924年就提出了關(guān)于大氣能見(jiàn)度與消光系數(shù)的關(guān)系,即大氣能見(jiàn)度與消光系數(shù)之比為常數(shù)[3].因此,明晰大氣消光系數(shù)的演化機(jī)理就成為對(duì)能見(jiàn)度研究的關(guān)鍵技術(shù)環(huán)節(jié).諸多研究表明,氣溶膠消光系數(shù)是大氣消光系數(shù)的主體,其占比通常都在90%以上[4-5].目前,針對(duì)氣溶膠消光系數(shù)或能見(jiàn)度的參數(shù)化方案主要包括統(tǒng)計(jì)模型和利用Mie散射的計(jì)算模型兩種.

    統(tǒng)計(jì)模型的核心是建立大氣消光系數(shù)或能見(jiàn)度與其影響因子之間的統(tǒng)計(jì)關(guān)系[6].早在1985年,美國(guó)Improve計(jì)劃就率先建立了PM2.5中各化學(xué)組分質(zhì)量濃度與氣溶膠消光系數(shù)之間的定量關(guān)系,這為早期的能見(jiàn)度預(yù)報(bào)提供了有效的計(jì)算方案[7].進(jìn)一步研究指出[6,8-9],由于PM2.5中某些化學(xué)組分的散射消光能力存在顯著的時(shí)空差異性,利用Improve關(guān)系式計(jì)算的散射系數(shù)可能大大偏離實(shí)際的觀測(cè)數(shù)據(jù).影響能見(jiàn)度的因素非常復(fù)雜,針對(duì)多區(qū)域的研究一致表明,顆粒物濃度尤其是細(xì)顆粒物質(zhì)量濃度在區(qū)域能見(jiàn)度下降中占據(jù)主要地位[9-10].另外,Liu等[11]系統(tǒng)性研究了京津冀地區(qū)氣溶膠消光系數(shù)隨相對(duì)濕度的變化特征,發(fā)現(xiàn)當(dāng)相對(duì)濕度小于85%時(shí),氣溶膠散射吸濕增長(zhǎng)因子(RH)平緩增長(zhǎng),而當(dāng)相對(duì)濕度大于85%時(shí),(RH)急劇增大,表明氣溶膠的吸濕增長(zhǎng)特性是導(dǎo)致高濕條件下能見(jiàn)度急劇下降的重要誘因.白永清等[12]分析了武漢地區(qū)能見(jiàn)度與PM2.5質(zhì)量濃度和相對(duì)濕度之間的非線性關(guān)系,指出細(xì)顆粒物濃度和相對(duì)濕度共同制約了能見(jiàn)度的變化.Chen等[13]利用HaChi外場(chǎng)觀測(cè)數(shù)據(jù)的研究發(fā)現(xiàn),相比于僅利用氣溶膠體積濃度和相對(duì)濕度構(gòu)建的消光系數(shù)參數(shù)化方案,考慮以氣溶膠數(shù)濃度的消光系數(shù)參數(shù)化方案可以顯著提升低能見(jiàn)度的模擬效果.

    近年來(lái),隨著環(huán)境監(jiān)測(cè)技術(shù)的發(fā)展,基于Mie散射模型的氣溶膠消光或能見(jiàn)度模型也得到了深入的研究和應(yīng)用.Chen等[14]利用在楊村觀測(cè)的粒度譜、不同模態(tài)的氣溶膠吸濕增長(zhǎng)因子和環(huán)境消光系數(shù),基于Mie散射計(jì)算氣溶膠消光系數(shù),得到的計(jì)算值與觀測(cè)值吻合較好.雖然基于Mie散射模型的氣溶膠消光系數(shù)計(jì)算模型物理意義清晰,但模型關(guān)鍵參數(shù)的取值存在較大經(jīng)驗(yàn)性和人為任意性[11,13-14].氣溶膠光學(xué)參數(shù)的實(shí)際變化是非常復(fù)雜的,其經(jīng)驗(yàn)取值是大氣消光系數(shù)計(jì)算不確定性的重要來(lái)源[5].為此,張智察等[15-16]利用免疫進(jìn)化算法先后提出干燥條件下氣溶膠復(fù)折射率(DACRI)和氣溶膠粒徑吸濕增長(zhǎng)因子(RH)的反演方法,并基于成都地區(qū)DACRI與主要顆粒物質(zhì)量濃度指標(biāo)的統(tǒng)計(jì)分析,進(jìn)一步構(gòu)建了適用于該區(qū)域的DACRI和(RH)的參數(shù)化方案.

    因此,著眼于改善目前統(tǒng)計(jì)模型與Mie散射模型中大氣能見(jiàn)度模擬所需的關(guān)鍵參數(shù)取值的不確定性[17-18],以提升大氣能見(jiàn)度的模擬精度,將DACRI和(RH)參數(shù)化方案進(jìn)行耦合,提出了基于Mie散射模型的能見(jiàn)度改進(jìn)算法,并在成都地區(qū)進(jìn)行了實(shí)例驗(yàn)證,以期為區(qū)域能見(jiàn)度預(yù)報(bào)能力和大氣環(huán)境管理水平的提升提供技術(shù)支持.

    1 材料與方法

    1.1 觀測(cè)儀器

    使用的數(shù)據(jù)包括成都市2017年10~12月由WS600一體式氣象站、AURORA-3000積分濁度計(jì)、AE-31黑碳儀以及GRIMM180環(huán)境顆粒物監(jiān)測(cè)儀獲取的逐時(shí)觀資料.

    相關(guān)儀器介紹如下:(1)AURORA-3000型濁度計(jì)(Ecotech公司,澳大利亞)觀測(cè)波長(zhǎng)為525nm,采樣頻率為5min/次,TSP切割頭,檢測(cè)范圍>0.25Mm-1,每24h進(jìn)行零點(diǎn)檢查, 24h零點(diǎn)漂移<±1%,每周用R134a氣體進(jìn)行跨度標(biāo)定,通過(guò)內(nèi)部溫濕度傳感器來(lái)控制濁度計(jì)內(nèi)部加熱系統(tǒng),使得儀器內(nèi)部腔室中氣溶膠相對(duì)濕度控制在40%以下,將其作為氣溶膠的干燥狀態(tài).(2)AE-31型黑碳檢測(cè)儀(Magee Scientific公司,美國(guó))觀測(cè)黑碳(BC)質(zhì)量濃度,數(shù)據(jù)采集頻率為5min/次.黑碳儀采用TSP切割頭,采樣頭與儀器連接中間增設(shè)硅膠管減少水分對(duì)黑碳測(cè)量的影響.濁度計(jì)和黑碳儀的監(jiān)測(cè)資料經(jīng)過(guò)質(zhì)量控制后統(tǒng)一處理為小時(shí)均值數(shù)據(jù).(3)GRIMM180環(huán)境顆粒物監(jiān)測(cè)儀(GRIMM公司,德國(guó))可以實(shí)時(shí)測(cè)量大氣中PM10, PM2.5和PM1的顆粒物質(zhì)量濃度以及31個(gè)粒徑段的氣溶膠數(shù)濃度,據(jù)此獲得氣溶膠粒子譜分布n[(dry)],其數(shù)據(jù)頻率為5min/次,其中各粒徑段粒子直徑的起始值分別為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.0mm.(4)氣象要素(大氣能見(jiàn)度和相對(duì)濕度RH)由WS600(LUFFT公司,德國(guó))一體式氣象站進(jìn)行監(jiān)測(cè);氣態(tài)污染物NO2體積濃度由化學(xué)發(fā)光NO, NO2-NO分析儀(Thermo 42i,美國(guó))進(jìn)行監(jiān)測(cè).

    1.2 觀測(cè)地點(diǎn)

    GRIMM180大氣顆粒物監(jiān)測(cè)儀觀測(cè)點(diǎn)位于成都市一環(huán)路聯(lián)益大廈(104°02¢E,30°39¢N)頂樓,距地面81m.其余儀器觀測(cè)點(diǎn)位于成都市環(huán)境保護(hù)科學(xué)研究院綜合大樓樓頂(30°39¢N,104°02¢E),距離地面21m,四周2km內(nèi)無(wú)高大建筑物,視野開(kāi)闊,周圍是集中居住區(qū).兩處監(jiān)測(cè)點(diǎn)位周圍無(wú)明顯大氣污染源,二者直線距離為410m,環(huán)境氣象條件基本一致.

    1.3 資料處理方法

    將上述監(jiān)測(cè)數(shù)據(jù)統(tǒng)一處理成小時(shí)均值數(shù)據(jù).剔除了出現(xiàn)降水、沙塵以及大風(fēng)現(xiàn)象所在日的全部數(shù)據(jù);其次,剔除儀器烘干后相對(duì)濕度仍大于40%的異常數(shù)據(jù),以排除水汽影響;最后,剔除超出界限值數(shù)據(jù),連續(xù)無(wú)變化數(shù)據(jù),缺測(cè)數(shù)據(jù)以及氣溶膠質(zhì)量濃度存在倒掛等異常數(shù)據(jù),由此獲得匹配樣本1145個(gè).

    大氣消光系數(shù)代表光線在大氣中傳播單位距離時(shí)的相對(duì)衰減率,在550nm波長(zhǎng)處的環(huán)境大氣消光系數(shù)ext(RH) (km-1)與大氣能見(jiàn)度(km)的關(guān)系見(jiàn)式(1)和式(2)[19].

    ext(RH)=sp(RH)+ap+sg+ag(2)

    式中:sp(RH),ap,sg,ag,sp分別是550nm波長(zhǎng)處環(huán)境條件下氣溶膠散射系數(shù),氣溶膠吸收系數(shù),干潔大氣散射系數(shù),干潔大氣吸收系數(shù),單位為Mm-1.

    依據(jù)文獻(xiàn)[20-21],按式(3)將AURORA-3000積分濁度計(jì)觀測(cè)的525nm波長(zhǎng)處的干氣溶膠散射系數(shù)sp, 525nm訂正得到550nm波長(zhǎng)處的干氣溶膠散射系數(shù)(sp),式中= 1.36,代表了成都市Angstrom波長(zhǎng)指數(shù)[22].

    根據(jù)吳兌等[2]針對(duì)AE-31型黑碳檢測(cè)儀的對(duì)比觀測(cè)試驗(yàn)原理,利用AE - 31型黑碳檢測(cè)儀直接觀測(cè)得到未經(jīng)訂正的880nm波長(zhǎng)處黑碳(BC)質(zhì)量濃度(BC)/μg·m-3,按式(4)反演532nm波長(zhǎng)處的氣溶膠吸收系數(shù)ap, 532nm(Mm-1),并由式(5)訂正得到550nm波長(zhǎng)處的吸收系數(shù)ap(Mm-1),見(jiàn)式(4).

    參照Penndorf(1957)[23]的研究成果,對(duì)應(yīng)550nm波長(zhǎng)處的sg一般取值為13Mm-1.參照Sloane (1985)[24]的計(jì)算方法,對(duì)應(yīng)550nm波長(zhǎng)處的ag的計(jì)算見(jiàn)式(6).

    式中,NO2為NO2質(zhì)量濃度(10-9g/m3).

    2 能見(jiàn)度計(jì)算方案的設(shè)計(jì)

    作為大氣消光系數(shù)的主體,氣溶膠散射系數(shù)極易受顆粒物質(zhì)量濃度、顆粒物化學(xué)組分以及氣溶膠吸濕性的影響,一直以來(lái)是大氣消光系數(shù)和能見(jiàn)度預(yù)報(bào)不確定性的重要來(lái)源.因此,氣溶膠散射系數(shù)的計(jì)算方法是決定能見(jiàn)度計(jì)算精度的關(guān)鍵.

    本文將前2/3時(shí)間序列的樣本(763)作為下文三種大氣能見(jiàn)度計(jì)算方案的建模數(shù)據(jù)集.后1/3時(shí)間序列的樣本(382)作為驗(yàn)證數(shù)據(jù)集,用于比較各個(gè)方案在大氣能見(jiàn)度模擬中的適用性.首先將大氣能見(jiàn)度劃分為<2km, 2~5km, 5~10km以及>10km這4種范圍,并對(duì)應(yīng)給出了這4種范圍大氣能見(jiàn)度條件下常規(guī)氣象要素的統(tǒng)計(jì)結(jié)果(表1).從表1可見(jiàn),大氣能見(jiàn)度越低,對(duì)應(yīng)的PM2.5, BC, RH總體越高.

    表1 4種范圍大氣能見(jiàn)度下常規(guī)氣象要素的統(tǒng)計(jì)結(jié)果

    注:(, PM2.5),(, BC),(, RH)分別表示與PM2.5質(zhì)量濃度, BC, RH之間的相關(guān)系數(shù).

    2.1 能見(jiàn)度計(jì)算的統(tǒng)計(jì)模型(方案一)

    氣溶膠散射吸濕增長(zhǎng)因子(RH)是環(huán)境條件下氣溶膠散射系數(shù)與對(duì)應(yīng)干燥條件下氣溶膠散射系數(shù)的比值.根據(jù)已有的研究成果[25],成都地區(qū)秋冬季(RH)的計(jì)算公式見(jiàn)式(7).

    根據(jù)式(2),式(7)以及sp和PM2.5之間的統(tǒng)計(jì)關(guān)系[26],ext(RH)的計(jì)算公式見(jiàn)式(8).

    根據(jù)式(8),利用前述建模數(shù)據(jù)求解各待定參數(shù),得到成都地區(qū)秋冬季ext(RH)的統(tǒng)計(jì)模型,見(jiàn)式(9).

    至此,利用式(9)計(jì)算ext(RH),再根據(jù)式(1)計(jì)算獲得了大氣能見(jiàn)度,建模計(jì)算結(jié)果表明,大氣能見(jiàn)度的計(jì)算值總體較觀測(cè)值略低3.81%,兩者之間的和MRE分別為0.94和18.01%.其中,4種大氣能見(jiàn)度模擬值與其觀測(cè)值之間的分別為0.66、0.75、0.59和0.83, MRE則分別為12.42%、18.56%、17.21%和19.64%.

    2.2 能見(jiàn)度計(jì)算的經(jīng)驗(yàn)參數(shù)Mie散射模型(方案二)

    基于Mie散射理論[15-16,27-28],ext(RH)的計(jì)算公式見(jiàn)式(10).

    (RH) =re(RH) +i(RH) (12)

    式中:(RH)是環(huán)境條件下氣溶膠的尺度參數(shù);re(RH)和i(RH)是環(huán)境條件下氣溶膠復(fù)折射率(DACRI)(RH)的實(shí)部和虛部;sp[(RH),(RH)]是由Mie散射模型計(jì)算出的環(huán)境條件下氣溶膠散射效率因子,由(RH)和(RH)決定;(RH)是環(huán)境條件下的氣溶膠粒徑;[(RH)]是環(huán)境條件下氣溶膠數(shù)濃度粒徑分布;(RH)是環(huán)境條件下氣溶膠的粒徑.環(huán)境條件下氣溶膠尺度參數(shù)(RH)和復(fù)折射率(RH)可進(jìn)一步分解為式(13-17)[29-30].

    (dry) =re(dry) +i(dry) (16)

    (water) =re(water) +i(water)=1.33+0×i (17)

    式中:(RH)為氣溶膠粒徑吸濕增長(zhǎng)因子;(dry)為干燥條件下氣溶膠的粒徑;re(dry)和i(dry)為干燥條件下氣溶膠復(fù)折射率(dry)的實(shí)部和虛部;re(water)和i(water)為水復(fù)折射率(water)的實(shí)部和虛部.

    由于DACRI觀測(cè)與反演的復(fù)雜性,其實(shí)部re(dry)和虛部i(dry)在以往研究中通常分別被經(jīng)驗(yàn)假定為1.55與-0.005[13].另外,基于孫景群等提出的(RH)參數(shù)化方案,見(jiàn)式(18).

    DACRI和(RH)是基于Mie散射模型模擬大氣能見(jiàn)度所必須的關(guān)鍵光學(xué)參數(shù),但受制于DACRI和(RH)觀測(cè)與反演的復(fù)雜性,以往研究中re(dry)、i(dry)與(RH)參數(shù)化方案的參數(shù)()通常為經(jīng)驗(yàn)取值.為此,本文首先采用式(18)的參數(shù)化方案來(lái)計(jì)算(RH),由于成都地區(qū)的氣溶膠類型為污染型,再參考孫景群與Kasten等的研究成果[29-31],本文將式(18)中的參數(shù)假定為4.4.

    在上述re(dry)、i(dry)與(RH)經(jīng)驗(yàn)取值的情況下,基于建模數(shù)據(jù)集,根據(jù)式(1)可求得大氣能見(jiàn)度,建模計(jì)算結(jié)果表明,re(dry)和i(dry)的經(jīng)驗(yàn)取值與觀測(cè)值之間不存在相關(guān)性(由于DACRI的經(jīng)驗(yàn)定值),對(duì)應(yīng)的MRE則分別為3.37%和77.80%,根據(jù)式(18)(=4.4的經(jīng)驗(yàn)取值)所計(jì)算的(RH)與觀測(cè)值之間的和MRE則分別為0.92和18.19%,根據(jù)式(18)計(jì)算的re(RH)和i(RH)與其觀測(cè)值之間的分別為0.58和0.47,對(duì)應(yīng)的MRE則分別為6.23%和82.15%,大氣能見(jiàn)度的計(jì)算值較觀測(cè)值平均低23.80%,其對(duì)應(yīng)的和MRE分別為0.96和24.53%.其中,4種大氣能見(jiàn)度模擬值與其觀測(cè)值之間的分別為0.71,0.81,0.68和0.91, MRE則分別為29.63%,26.01%,26.81%和18.69%.從方案二的建模計(jì)算結(jié)果可見(jiàn),re(dry)經(jīng)驗(yàn)取值與其觀測(cè)值之間的誤差較小,但i(dry),(RH)以及大氣能見(jiàn)度與其觀測(cè)值之間的誤差均很大.

    2.3 能見(jiàn)度計(jì)算的雙參數(shù)化方案Mie散射模型(方案三)

    針對(duì)方案二中DACRI與(RH)經(jīng)驗(yàn)取值的不確定性,張智察(2020)在分析DACRI與顆粒物質(zhì)量濃度統(tǒng)計(jì)關(guān)系的基礎(chǔ)上,利用多元逐步線性回歸方法,提出了DACRI的參數(shù)化方案,見(jiàn)式(19)和式(20).

    式中:BC,PM1,PM2.5和PM10分別為BC, PM1, PM2.5和PM10的質(zhì)量濃度,單位為μg/m3.

    不同類型氣溶膠粒徑吸濕性差異很大,即便是同一類型的氣溶膠,受排放源以及氣溶膠物化過(guò)程的影響,氣溶膠粒徑吸濕性也有差異.為此,張智察等[16]通過(guò)結(jié)合免疫進(jìn)化算法與Mie散射理論反演(RH),分析了其與RH之間的關(guān)系,并基于式(18)提出了適用于成都地區(qū)秋冬季的(RH)參數(shù)化方案,見(jiàn)式(21).

    基于建模數(shù)據(jù)集,首先分別根據(jù)式(19)和式(20)分別計(jì)算re(dry)和i(dry),式(21)計(jì)算(RH).在此基礎(chǔ)上,大氣能見(jiàn)度的計(jì)算流程與方案二相同.最終計(jì)算結(jié)果表明,根據(jù)式(19)和式(20)計(jì)算的re(dry)和i(dry)與觀測(cè)值之間的分別為0.55和0.85,對(duì)應(yīng)的MRE分別為2.27%和14.72%,根據(jù)式(21)計(jì)算的(RH)與觀測(cè)值之間的和MRE則分別為5.03%和0.92,根據(jù)式(18)計(jì)算的re(RH)和i(RH)與其觀測(cè)值之間的分別為0.70和0.85,對(duì)應(yīng)的MRE則分別為1.46%和18.98%,大氣能見(jiàn)度的計(jì)算值較觀測(cè)值平均略低4.92%,其對(duì)應(yīng)的和MRE分別為0.93和14.26%.其中,4種大氣能見(jiàn)度模擬值與其觀測(cè)值之間的分別為0.66,0.79,0.65和0.91, MRE則分別為10.45%,13.72%,16.59%和13.21%.可見(jiàn)方案三針對(duì)復(fù)折射率和粒徑吸濕增長(zhǎng)因子的建模計(jì)算結(jié)果要顯著優(yōu)于方案二,且針對(duì)大氣能見(jiàn)度的建模計(jì)算效果則要顯著優(yōu)于方案一與方案二.

    3 三種能見(jiàn)度計(jì)算方案的適用性分析

    3.1 方案一

    圖1 不同范圍大氣能見(jiàn)度方案一模擬值與其觀測(cè)值之間的對(duì)比(a: <2km; b: 2~5km; c: 5~10km; d: >10km)

    基于驗(yàn)證數(shù)據(jù)集,根據(jù)方案一模擬了大氣能見(jiàn)度,并分別給出了上述不同范圍大氣能見(jiàn)度模擬值與其觀測(cè)值之間的散點(diǎn)圖(圖1). 模擬結(jié)果表明,4種范圍(<2km, 2~5km, 5~10km以及>10km)的大氣能見(jiàn)度模擬值與其觀測(cè)值之間的分別為0.72, 0.80, 0.64, 0.84, MRE分別為21.63%, 13.91%, 16.00%, 14.83%.并隨著大氣能見(jiàn)度范圍的增大,方案一的模擬效果逐漸降低,小于2km的低能見(jiàn)度(圖1a)模擬誤差最大(21.63%).結(jié)合表1與圖1分析可知,小于2km的低能見(jiàn)度與RH之間的(, RH)為-0.35,顯著大于此時(shí)大氣能見(jiàn)度與PM2.5和BC之間的(, PM2.5)和(, BC),因此推測(cè)高相對(duì)濕度條件下氣溶膠吸濕增長(zhǎng)的不確定性可能是引起統(tǒng)計(jì)模型模擬大氣能見(jiàn)度誤差的重要因素.

    3.2 方案二

    基于驗(yàn)證數(shù)據(jù)集,根據(jù)方案二中re(dry)、i(dry)與的經(jīng)驗(yàn)取值,首先根據(jù)式(18)及其=4.4的經(jīng)驗(yàn)取值模擬獲得了對(duì)應(yīng)的(RH),其模擬值與觀測(cè)值之間的和MRE分別為0.91和7.81%.雖然氣溶膠粒徑吸濕增長(zhǎng)模型的參數(shù)為經(jīng)驗(yàn)值4.4,但(RH)模擬值與其觀測(cè)值之間的散點(diǎn)分布(圖2)也能較好地體現(xiàn)不同RH條件下氣溶膠粒徑的吸濕增長(zhǎng)特征,且誤差尚在可接受范圍之內(nèi).

    圖2 Gf(RH)方案二模擬值與其觀測(cè)值的散點(diǎn)

    其次,re(dry)和i(dry)的經(jīng)驗(yàn)取值與觀測(cè)值之間的均為0, MRE則分別為2.79%和61.19%.在此基礎(chǔ)上,基于前述(RH)的模擬結(jié)果,根據(jù)式(13-18))分別模擬了ACRI(圖3)的re(RH)(圖3c)和i(RH)(圖3d),模擬結(jié)果表明,re(RH)和i(RH)的模擬值與其觀測(cè)值之間的分別為0.46和0.39, MRE分別為3.04%和78.78%.從上述模擬結(jié)果可知,即使re(dry)和i(dry)為經(jīng)驗(yàn)取值,其本身存在著很大的不確定性,但氣溶膠吸濕增長(zhǎng)之后,相應(yīng)光學(xué)參數(shù)變化的不確定性也會(huì)進(jìn)一步加大.另外,環(huán)境相對(duì)濕度越大,re(RH)和i(RH)與其觀測(cè)值之間的相關(guān)系數(shù)也越高,表明相對(duì)濕度的增大會(huì)加強(qiáng)相對(duì)濕度對(duì)于折射率的決定作用.

    根據(jù)前述模擬所得的(RH)、re(RH)和i(RH),依據(jù)方案二模擬獲得了大氣能見(jiàn)度,據(jù)此分別給出了上述不同范圍大氣能見(jiàn)度模擬值與其觀測(cè)值之間的散點(diǎn)圖(圖4).4種范圍(<2km, 2~5km, 5~10km以及>10km)的大氣能見(jiàn)度模擬值與其觀測(cè)值之間的分別為0.68, 0.91, 0.90, 0.92, MRE分別為22.79%, 25.06%, 22.85%, 12.01%.由此可見(jiàn),(RH)、re(RH)和i(RH)的經(jīng)驗(yàn)取值會(huì)導(dǎo)致大氣能見(jiàn)度模擬的重大誤差.另外,對(duì)比方案一和方案二的模擬結(jié)果可推測(cè),以相對(duì)濕度作為因變量的計(jì)算ACRI和(RH)的參數(shù)化方案,在進(jìn)一步結(jié)合Mie散射模型后,可更全面地表征大氣能見(jiàn)度的非線性演變特征,而且在不主要考慮顆粒物質(zhì)量濃度的本文前提下,相對(duì)濕度在大氣能見(jiàn)度的非線性演變中起著決定性作用,因此方案二中的要顯著高于方案一.當(dāng)然,由于(RH)、re(RH)和i(RH)來(lái)源于相關(guān)的經(jīng)驗(yàn)取值,這不可避免的會(huì)引起很大的模擬誤差.

    圖4 不同范圍大氣能見(jiàn)度方案二模擬值與其觀測(cè)值之間的對(duì)比

    3.3 方案三

    圖5 Gf(RH)方案三模擬值與其觀測(cè)值之間的對(duì)比

    基于驗(yàn)證數(shù)據(jù)集,首先根據(jù)方案三中(RH)的成都秋冬季本地化參數(shù)化方案(式(21))模擬獲得了驗(yàn)證集的(RH),其模擬值與觀測(cè)值之間的和MRE分別為0.81和4.24%.相較于方案二中的經(jīng)驗(yàn)參數(shù)化方案而言,采用本地化參數(shù)化方案模擬獲得的(RH)精度更高.

    基于上述(RH)的計(jì)算結(jié)果,再根據(jù)式(19)和式(20)分別模擬了干燥條件下氣溶膠復(fù)折射率DACRI的re(dry) (圖6a)和i(dry) (圖6b),各自的模擬值與其觀測(cè)值之間的分別為0.21和0.85, MRE分別為2.31%和13.36%.在DACRI和(RH)模擬結(jié)果的基礎(chǔ)上,根據(jù)式(13~18)分別模擬了環(huán)境條件下氣溶膠復(fù)折射率ACRI的re(RH) (圖6c)和i(RH) (圖6d),各自的模擬值與其觀測(cè)值之間的分別為0.51和0.83, MRE分別為3.46%和18.84%,可見(jiàn)采用了本地化參數(shù)化方案的復(fù)折射率模擬效果有了顯著的提升.

    張智察等研究中已指出,re(dry)和i(dry)對(duì)參數(shù)化方案中自變量因子敏感程度的不同是導(dǎo)致兩者模擬值與觀測(cè)值之間和MRE存在差異的重要因素.此外,本文的DACRI參數(shù)化方案雖然可以有效地模擬出re和i,但在氣溶膠化學(xué)組分信息的表征上可能仍不夠精細(xì),這很可能是導(dǎo)致ACRI模擬結(jié)果存在不確定性的主要原因.通過(guò)對(duì)比圖6a與圖6c以及圖6b和圖6d可知,re(dry)到re(RH)的模擬誤差從2.31%提高至3.46%,對(duì)應(yīng)i(dry)到i(RH)的計(jì)算誤差則從13.36%提高至18.84%,這表明氣溶膠的吸濕增長(zhǎng)過(guò)程會(huì)增大ACRI,(RH)等關(guān)鍵氣溶膠光學(xué)輻射參數(shù)的模擬誤差,這一推論于方案二中針對(duì)大氣能見(jiàn)度的模擬情況也有較好的體現(xiàn),諸多研究結(jié)論也較好地驗(yàn)證了這一點(diǎn)[24,32-35].

    圖6 ACRI方案三模擬值與觀測(cè)值之間的對(duì)比

    a:re(dry); b:i(dry); c:re(RH); d:i(RH)

    根據(jù)前述參數(shù)化方案模擬所得的(RH)和ACRI,再依據(jù)方案三模擬獲得了大氣能見(jiàn)度,并分別給出了上述不同范圍大氣能見(jiàn)度模擬值與其觀測(cè)值之間的散點(diǎn)圖(圖7). 4種范圍(<2km, 2~5km, 5~ 10km以及>10km)的大氣能見(jiàn)度模擬值與其觀測(cè)值之間的分別為0.62, 0.90, 0.89, 0.93, MRE分別為9.86%, 10.39%, 9.94%, 14.06%.由此可見(jiàn),較方案一和方案二而言,方案三的模擬效果在各方面均要顯著占優(yōu).另外,在不主要考慮顆粒物質(zhì)量濃度的本文前提下,方案三的與方案二基本一致,但模擬精度顯著更優(yōu),這也表明大氣能見(jiàn)度非線性演變特征受相對(duì)濕度的主導(dǎo)程度較大,光學(xué)參數(shù)的準(zhǔn)確性則主導(dǎo)了大氣能見(jiàn)度的模擬精度.

    通過(guò)耦合DACRI和(RH)的成都地區(qū)本地化參數(shù)化方案來(lái)估算氣溶膠光學(xué)輻射參數(shù),并基于Mie散射模型來(lái)模擬大氣能見(jiàn)度,最終取得了比傳統(tǒng)方法更優(yōu)的模擬效果,這反映了本文DACRI和(RH)的參數(shù)化方案在氣溶膠光學(xué)輻射強(qiáng)迫效應(yīng)模擬中的適用性.隨著大氣化學(xué)模式(GEOS-Chem/ WRF-Chem/WRF-CMAQ等)以及大數(shù)據(jù)機(jī)器學(xué)習(xí)算法的不斷發(fā)展,目前大氣污染物(BC, PM1, PM2.5, PM10and NO2),相對(duì)濕度RH以及氣溶膠數(shù)濃度粒徑分布()的可預(yù)報(bào)性也在不斷提高,這為大氣化學(xué)模式的改進(jìn)提供一定參考,例如在模式的輸出數(shù)據(jù)中直接運(yùn)用新參數(shù)化方案估算氣溶膠參數(shù),再進(jìn)行氣溶膠光學(xué)輻射效應(yīng)的模擬.

    圖7 不同范圍大氣能見(jiàn)度方案三模擬值與其觀測(cè)值之間的對(duì)比

    4 結(jié)論

    4.1 經(jīng)驗(yàn)參數(shù)的Mie散射模型、統(tǒng)計(jì)模型以及耦合氣溶膠復(fù)折射率(DACRI)和氣溶膠粒徑吸濕增長(zhǎng)因子((RH))參數(shù)化方案的大氣能見(jiàn)度數(shù)值改進(jìn)算法均能較好地模擬出能見(jiàn)度的變化特征.

    4.2 改進(jìn)算法通過(guò)本地化參數(shù)化方案更準(zhǔn)確地估計(jì)出DACRI和(RH),從而可更準(zhǔn)確地模擬出四類能見(jiàn)度區(qū)間(<2km, 2~5km, 5~10km, >10km),對(duì)應(yīng)模擬值與觀測(cè)值之間的相關(guān)系數(shù)分別為0.62,0.90, 0.89,0.93,平均相對(duì)誤差(MRE)分別為9.86%, 10.39%,9.94%,14.06%.

    [1] 張小曳.中國(guó)大氣氣溶膠及其氣候效應(yīng)的研究[J]. 地球科學(xué)進(jìn)展, 2007,22(1):12-16.

    Zhang X Y.Aerosol over China and their climate effect [J]. Advances in Earth Science, 2007,22(1):12-16.

    [2] 吳 兌,毛節(jié)泰,鄧雪嬌,等.珠江三角洲黑碳?xì)馊苣z及其輻射特性的觀測(cè)研究[J]. 中國(guó)科學(xué):地球科學(xué), 2009,(11):1542-1553.

    Wu D, Mao J T, Deng X J, et al. Black carbon aerosols and their radiative properties in the Pearl River Delta region [j]. Sci. China Ser. D-Earth Sci., 2009,52(8):1152-1163.

    [3] Koschmieder H. Theorie der horizontalen Sichtweite [J]. Beitrage zur Physik der freien Atmosphare, 1924:33-53.

    [4] 劉新罡,張遠(yuǎn)航.基于觀測(cè)的大氣氣溶膠散射吸濕增長(zhǎng)因子模型研究—以2006CAREBeijing加強(qiáng)觀測(cè)為例[J]. 中國(guó)環(huán)境科學(xué), 2009,29(12):1243-1248.

    Liu X G, Zhang Y H.Modelling research on the aerosol scattering hygroscopic growth factor based on measurement—Taking 2006 CAREBeijing campaign for example [J]. China Environment Science, 2009,29(12):1243-1248.

    [5] 楊寅山,倪長(zhǎng)健,鄧 也,等.成都市冬季大氣消光系數(shù)及其組成的特征研究[J]. 環(huán)境科學(xué)學(xué)報(bào), 2019,39(5):1425-1432.

    Yang Y S, Ni C J, Deng Y, et al. Characteristics of atmospheric extinction coefficient and its components in winter in Chengdu [J]. Acta Scientiae Circumstantiae, 2019,39(5):1425-1432.

    [6] Pitchford M, Maim W, Schichtel B, et al. Revised algorithm for estimating light extinction from IMPROVE particle speciation data [J]. Journal of the Air & Waste Management Association, 2007,57(11): 1326-1336.

    [7] Malm W C, et al. Spatial and seasonal patterns and temporal variability of haze and its constituents in the United States reports IV [R]. 2006:1-251.

    [8] 馬乾坤,成春雷,李 梅,等.鶴山氣溶膠光學(xué)性質(zhì)和單顆?;瘜W(xué)組分的研究[J]. 中國(guó)環(huán)境科學(xué), 2019,39(7):2710-2720.

    Ma Q K, Cheng C L, Li M, et al. The aerosol optical characteristics and chemical composition of single particles in Heshan [J]. China Environmental Science, 2019,39(7):2710-2720.

    [9] 江家豪,彭 杏,朱 波,等. 深圳大氣PM2.5化學(xué)組成的長(zhǎng)期變化特征[J]. 中國(guó)環(huán)境科學(xué), 2021,41(2):574-579.

    Jiang J H, Peng X, Zhu B, et al. Long-term variational characteristics of the chemical composition of PM2.5in Shenzhen. China Environmental Science, 2021,41(2):574-579.

    [10] 黃元龍.上海市大氣氣溶膠光學(xué)與化學(xué)性質(zhì)研究[D]. 上海:復(fù)旦大學(xué), 2014.

    Huang Y L. Field studies of aerosol optical and chemical properties in Shanghai [D]. Shanghai:Fudan University, 2014.

    [11] Liu X, Zhang Y, Jung J, et al. Research on the hygroscopic properties of aerosols by measurement and modeling during CARE Beijing-2006 [J]. Journal of Geophysical Research Atmospheres, 2009,114 (D00G16).

    [12] 白永清,祁海霞,劉 琳,等.武漢大氣能見(jiàn)度與PM2.5濃度及相對(duì)濕度關(guān)系的非線性分析及能見(jiàn)度預(yù)報(bào)[J]. 氣象學(xué)報(bào), 2016,74(2): 189-199.

    Bai Y Q, Qi H X, Liu L, et al. Study on the nonlinear relationship among the visibility, PM2.5concentration and relative humidity in Wuhan and the visibility prediction [J]. Journal of Meteorology, 2016, 74(2):189-199.

    [13] Chen J, Zhao C S, Ma N, et al. A parameterization of low visibilities for hazy days in the North China Plain [J]. Atmospheric Chemistry & Physics Discussions, 2012,12(11):4935-4950.

    [14] Chen J, Zhao C S, Ma N, et al. Aerosol hygroscopicity parameter derived from the light scattering enhancement factor measurements in the North China Plain [J]. Atmospheric Chemistry & Physics, 2014, 14(3):8105-8118.

    [15] 張智察,倪長(zhǎng)健,鄧 也,等.氣溶膠等效復(fù)折射率反演的免疫進(jìn)化算法[J]. 中國(guó)環(huán)境科學(xué), 2019,39(2):554-559.

    Zhang Z C, Ni C J, Deng Y, et al. Retrieval of equivalent complex refractive index of aerosol particles based on immune evolution algorithm [J]. China Environmental Science, 2019,39(2):554-559.

    [16] 張智察,倪長(zhǎng)健,鄧 也,等.免疫進(jìn)化算法反演均勻混合氣溶膠吸濕增長(zhǎng)因子[J]. 中國(guó)環(huán)境科學(xué), 2020,40(3):82-89.

    Zhang Z C, Ni C J, Deng Y, et al. Retrieval of hygroscopic growth factorof uniformly mixedaerosol particles based on immune evolution algorithm [J]. China Environmental Science, 2020,40(3):82-89.

    [17] 周 鑫,丁菊麗,田 偉,等.能見(jiàn)度參數(shù)化方案概述[J]. 氣象研究與應(yīng)用, 2009,30(3):20-23.

    Zhou X, Ding J L, Tian W, et al. introduction of visibility parameterization scheme [J]. Journal of Meteorological Research and Application, 2009,30(3):20-23.

    [18] 胡 俊,趙天良,張澤鋒,等.霾污染環(huán)境大氣能見(jiàn)度參數(shù)化方案的改進(jìn)[J]. 環(huán)境科學(xué)研究, 2017,30(11):1680-1688.

    Hu J, Zhao T L, Zhang Z F, et al. Upgrading atmospheric visibility parameterization scheme for haze pollution environment [J]. Research of Environmental Sciences, 2017,30(11):1680-1688.

    [19] Kotchenruther R A, Hobbs P V, Hegg D A. Humidification factors for atmospheric aerosols off the mid-Atlantic coast of the United States [J]. Journal of Geophysical Research Atmospheres, 1999,104(D2): 2239-2251.

    [20] Bodhaine B. Aerosol absorption measurements at Barrow, Mauna Loa and the south pole [J]. Journal of Geophysical Research Atmospheres, 1995,100(D5):8967-8975.

    [21] 伯廣宇,劉 東,吳德成,等.雙波長(zhǎng)激光雷達(dá)探測(cè)典型霧霾氣溶膠的光學(xué)和吸濕性質(zhì)[J]. 中國(guó)激光, 2014,41(1):0113001.

    Bo G Y, Liu D, Wu D C, et al. Two -wavelength lidar for observation of aerosol optical and hygroscopic properties in fog and haze days [J]. Chinese Journal of Lasers,2014,41(1):0113001.

    [22] 李梅芳,葉芝祥.基于太陽(yáng)光度計(jì)的成都雙流地區(qū)夏季氣溶膠光學(xué)特性研究[J]. 成都信息工程學(xué)院學(xué)報(bào), 2014,29(2):213-216.

    Li M F, Ye Z X. The studies of aerosol optical properties of Chengdu Shuangliu in summer based on the sun photometer [J]. Journal of Chengdu University of Information Technology, 2014,29(2):213-216.

    [23] Penndorf Rudolf. Tables of the refractive index for standard air and the rayleigh scattering coefficient for the spectral region between 0.2 and 20.0μ and their application to atmospheric optics [J]. Journal of the Optical Society of America, 1957,47(2):176-182.

    [24] Sloane C S, Wolff G T. Prediction of ambient light scattering using a physical model responsive to relative humidity: Validation with measurements from Detroit [J]. Atmospheric Environment, 1985, 19(4):669-680.

    [25] 張智察,倪長(zhǎng)健,尹單丹,等.兩種氣溶膠消光吸濕增長(zhǎng)因子的適用性分析[J]. 激光與光電子學(xué)進(jìn)展, 2020,57(9):090103.

    Zhang Z C, Ni C J, Yin D D, et al. Applicability of the two kinds of aerosol extinction hygroscopic growth factors [J]. Laser & Optoelectronics Progress, 2020,57(9):090103.

    [26] 陳一娜,趙普生,何 迪,等.北京地區(qū)大氣消光特征及參數(shù)化研究[J]. 環(huán)境科學(xué), 2015,36(10):3582-3589.

    Chen Y N, Zhao P S, He D, et al. Characteristics and parameterization for atmospheric extinction coefficient in Beijing [J]. Environmental Science, 2015,36(10):3582-3589.

    [27] Bohren C F, Huffman D R. Absorption and scattering of light by small particles [M]. John Wiley & Sons Inc. 1983.

    [28] 沈建琪,劉 蕾.經(jīng)典Mie散射的數(shù)值計(jì)算方法改進(jìn)[J]. 中國(guó)粉體技術(shù), 2005,11(4):45-50.

    Shen J Q, Liu L. An improved algorithm of classical Mie scattering calculation [J]. China Powder Science and Technology, 2005,11(4): 45-50.

    [29] 孫景群.濕氣溶膠的光散射特性[J]. 高原氣象, 1983,2(3):49-54.

    Sun J Q. Relationship between visibility and relative humidity [J]. Plateau Meteorology, 1983,2(3):49-54.

    [30] 孫景群.能見(jiàn)度與相對(duì)濕度的關(guān)系[J]. 氣象學(xué)報(bào), 1985,43(2):230- 234.

    Sun J Q. The relationship between visibility and relative humidity [J]. Journal of Meteorology, 1985,43(2):230-234.

    [31] Kasten F. Visibility forecast in the phase of pre-condensation [J]. Tellus, 1969,21(5):631-635.

    [32] Ycc A, Rws A, Ghm A, et al. Source apportionment of visibility degradation problems in Brisbane (Australia) using the multiple linear regression techniques [J]. Atmospheric Environment, 1999,33(19): 3237-3250.

    [33] Ebert M, Weinbruch S, Rausch A, et al. The complex refractive index of aerosols during LACE 98 as derived from the analysis of indivi-dual particles [J]. Journal of Geophysical Research: Atmospheres, 2002,107(D21):LAC 3.

    [34] Wex H, Neusü C, Koziar C, et al. Particle scattering, backscattering, and absorption coefficients: An in situ closure and sensitivity study [J]. Journal of Geophysical Research Atmospheres, 2002,107(21):LAC 4.

    [35] Cheng Y F, Eichler H, Wiedensohler A, et al. Mixing state of elemental carbon and non-light-absorbing aerosol components derived from in situ particle optical properties at Xinken in Pearl River Delta of China [J]. Journal of Geophysical Research Atmospheres, 2006,111(20): 4763-4773.

    An improved numerical algorithm for simulating atmosphere visibility by coupling two aerosol parameterization schemes.

    ZHANG Zhi-cha1,2, NI Chang-jian1*, ZHAO Jun-ping2, FENG Miao3, XIE Yan-xin1, WANG Yang-fei1

    (1.School of Atmospheric Sciences, Chengdu University of Information Technology, Chengdu 610225, China;2.Zhejiang Meteorological Observatory, Hangzhou 310017, China;3.Chengdu Academy of Environmental Sciences, Chengdu 610072, China)., 2021,41(11):5009~5018

    Based on the Mie theory, an improved numerical algorithm for simulating atmosphere visibility was proposed by coupling the parameterization schemes of aerosol complex refractive index (DACRI) and aerosol hygroscopic growth factor(RH). By utilizing the monitored data at an hourly time step recorded by WS600 integrated weather station, AURORA-3000 integrating nephelometer, AE-31aethalometer and GRIMM180 environment particle monitors from October to December in 2017 in Chengdu, the applicability of the improved algorithm was futher evaluated by compared with those of two visibility calculation methods, which are Mie theory model with empirical parameters and statistical model, in different visibility intervals (<2km, 2~5km, 5~10km, >10km). The results show that the three visibility calculation methods can capture the variation characteristics of visibility well; The improved algorithm better simulate visibility by adopting localization parameterization schemes to estimate DACRI and(RH). The corresponding correlation coefficients () between the simulation and the observation are 0.62, 0.90, 0.89, and 0.93, respectively. The mean relative error (MRE) are 9.86%, 10.39%, 9.94%, and 14.06%, respectively.

    aerosol;complex refractive index;hygroscopic growth factor;parameterization scheme;visibility

    X513

    A

    1000-6923(2021)11-5009-10

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

    2021-03-30

    國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2018YFC0214004;2018YFC1506006);四川省科技廳應(yīng)用基礎(chǔ)研究(2021YJ0314)

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

    猜你喜歡
    顆粒物大氣模型
    一半模型
    大氣的呵護(hù)
    軍事文摘(2023年10期)2023-06-09 09:15:06
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    南平市細(xì)顆粒物潛在來(lái)源分析
    3D打印中的模型分割與打包
    大氣古樸揮灑自如
    大氣、水之后,土十條來(lái)了
    錯(cuò)流旋轉(zhuǎn)填料床脫除細(xì)顆粒物研究
    多層介質(zhì)阻擋放電處理柴油機(jī)尾氣顆粒物
    国产精品偷伦视频观看了| 亚洲色图 男人天堂 中文字幕| tocl精华| 中文字幕人妻熟女乱码| 久久精品国产综合久久久| 人妻人人澡人人爽人人| 啦啦啦中文免费视频观看日本| 国产高清videossex| 50天的宝宝边吃奶边哭怎么回事| 制服人妻中文乱码| 中文精品一卡2卡3卡4更新| 日本一区二区免费在线视频| 最新在线观看一区二区三区| 国产精品免费大片| 国产精品一区二区免费欧美 | 丝袜美足系列| av不卡在线播放| 国产视频一区二区在线看| 日韩中文字幕欧美一区二区| 丝瓜视频免费看黄片| 国产日韩欧美亚洲二区| 欧美午夜高清在线| 人人妻人人澡人人爽人人夜夜| a级毛片在线看网站| 18在线观看网站| 国产亚洲欧美精品永久| 男人爽女人下面视频在线观看| 一个人免费看片子| 欧美精品人与动牲交sv欧美| 丰满人妻熟妇乱又伦精品不卡| 亚洲国产精品成人久久小说| 91麻豆精品激情在线观看国产 | 亚洲欧美激情在线| 国产xxxxx性猛交| 国产精品av久久久久免费| 中文欧美无线码| 久久人妻熟女aⅴ| 大香蕉久久成人网| 中文字幕色久视频| 久热爱精品视频在线9| 亚洲国产看品久久| 视频在线观看一区二区三区| 久久久久久久久免费视频了| 999久久久国产精品视频| 日本黄色日本黄色录像| 日本撒尿小便嘘嘘汇集6| 精品高清国产在线一区| 嫩草影视91久久| 美女高潮到喷水免费观看| 99香蕉大伊视频| 国产精品欧美亚洲77777| 女性被躁到高潮视频| 搡老岳熟女国产| 免费观看av网站的网址| 熟女少妇亚洲综合色aaa.| 91九色精品人成在线观看| 亚洲精品一区蜜桃| 久久99一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 美女福利国产在线| 亚洲av成人不卡在线观看播放网 | 国产老妇伦熟女老妇高清| 极品少妇高潮喷水抽搐| 精品福利永久在线观看| www.熟女人妻精品国产| 欧美亚洲日本最大视频资源| 国产av精品麻豆| e午夜精品久久久久久久| 搡老熟女国产l中国老女人| 极品少妇高潮喷水抽搐| 老鸭窝网址在线观看| 男女国产视频网站| 亚洲国产欧美日韩在线播放| 精品国产一区二区三区久久久樱花| 亚洲伊人久久精品综合| 热99国产精品久久久久久7| 欧美97在线视频| 纵有疾风起免费观看全集完整版| 多毛熟女@视频| 一级a爱视频在线免费观看| av视频免费观看在线观看| 最新的欧美精品一区二区| 搡老乐熟女国产| 午夜两性在线视频| 亚洲精品一二三| 黄片大片在线免费观看| 少妇猛男粗大的猛烈进出视频| 嫩草影视91久久| 国产精品自产拍在线观看55亚洲 | 国产极品粉嫩免费观看在线| 不卡av一区二区三区| 男人添女人高潮全过程视频| 美女大奶头黄色视频| 正在播放国产对白刺激| 亚洲中文字幕日韩| 久久精品人人爽人人爽视色| 日韩中文字幕视频在线看片| 日韩中文字幕欧美一区二区| 一个人免费在线观看的高清视频 | 欧美 日韩 精品 国产| 五月开心婷婷网| 99热国产这里只有精品6| 精品卡一卡二卡四卡免费| 精品人妻1区二区| 精品国产乱码久久久久久男人| 成年女人毛片免费观看观看9 | 精品一区在线观看国产| 亚洲专区中文字幕在线| 国产精品自产拍在线观看55亚洲 | 免费在线观看黄色视频的| 91精品伊人久久大香线蕉| 精品亚洲成a人片在线观看| 免费高清在线观看日韩| 99精品久久久久人妻精品| 国精品久久久久久国模美| 欧美人与性动交α欧美软件| 国产日韩欧美在线精品| 亚洲黑人精品在线| 亚洲欧美日韩高清在线视频 | 亚洲性夜色夜夜综合| 18禁裸乳无遮挡动漫免费视频| 欧美激情久久久久久爽电影 | a在线观看视频网站| 色视频在线一区二区三区| 狠狠精品人妻久久久久久综合| 亚洲人成77777在线视频| 精品久久蜜臀av无| 激情视频va一区二区三区| 黄色视频不卡| 高清av免费在线| 亚洲中文字幕日韩| 性高湖久久久久久久久免费观看| 日韩中文字幕视频在线看片| 欧美成人午夜精品| 成人亚洲精品一区在线观看| 夫妻午夜视频| 中文字幕制服av| 丝袜美腿诱惑在线| svipshipincom国产片| 妹子高潮喷水视频| 中文字幕人妻丝袜制服| 正在播放国产对白刺激| 欧美日韩亚洲高清精品| 国产日韩欧美在线精品| 成人亚洲精品一区在线观看| 国产亚洲欧美精品永久| 黄色 视频免费看| 亚洲av电影在线进入| 亚洲第一欧美日韩一区二区三区 | 婷婷成人精品国产| 国产欧美亚洲国产| 国产精品免费大片| 欧美+亚洲+日韩+国产| 亚洲欧美精品自产自拍| 看免费av毛片| 国产xxxxx性猛交| 嫁个100分男人电影在线观看| av国产精品久久久久影院| 亚洲少妇的诱惑av| 丝袜人妻中文字幕| 免费在线观看视频国产中文字幕亚洲 | 俄罗斯特黄特色一大片| 成人影院久久| 免费在线观看影片大全网站| 亚洲一码二码三码区别大吗| 成人18禁高潮啪啪吃奶动态图| 婷婷色av中文字幕| 岛国毛片在线播放| 国产亚洲午夜精品一区二区久久| 国产深夜福利视频在线观看| 国产成人影院久久av| 91精品国产国语对白视频| av在线老鸭窝| 99香蕉大伊视频| 不卡av一区二区三区| 精品乱码久久久久久99久播| 成年女人毛片免费观看观看9 | 国产精品一区二区免费欧美 | av一本久久久久| 欧美人与性动交α欧美软件| 欧美黄色片欧美黄色片| 亚洲国产欧美在线一区| 久久久久久久国产电影| 中亚洲国语对白在线视频| 啦啦啦视频在线资源免费观看| 国产成人系列免费观看| 精品熟女少妇八av免费久了| 少妇 在线观看| 国产精品.久久久| 亚洲精品国产一区二区精华液| 欧美亚洲日本最大视频资源| 一级片免费观看大全| 91老司机精品| avwww免费| 一区在线观看完整版| 亚洲欧美清纯卡通| 十八禁网站免费在线| 日韩中文字幕视频在线看片| 纯流量卡能插随身wifi吗| 国产男女超爽视频在线观看| 91老司机精品| 黄色视频,在线免费观看| 最新在线观看一区二区三区| 久久久精品区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品久久午夜乱码| 在线永久观看黄色视频| 亚洲成国产人片在线观看| 日韩免费高清中文字幕av| 99久久精品国产亚洲精品| 国产精品久久久久成人av| 18禁国产床啪视频网站| 亚洲国产欧美在线一区| 丝瓜视频免费看黄片| 最近最新免费中文字幕在线| 亚洲成人手机| 亚洲男人天堂网一区| 在线av久久热| 国产欧美日韩一区二区三区在线| 超色免费av| 亚洲欧美一区二区三区久久| 国产91精品成人一区二区三区 | 91九色精品人成在线观看| a级片在线免费高清观看视频| 99久久精品国产亚洲精品| 日韩一区二区三区影片| 男女免费视频国产| 精品卡一卡二卡四卡免费| 1024香蕉在线观看| 亚洲avbb在线观看| 涩涩av久久男人的天堂| 欧美日韩福利视频一区二区| 亚洲成人免费电影在线观看| 电影成人av| 王馨瑶露胸无遮挡在线观看| 成年人免费黄色播放视频| 亚洲国产欧美在线一区| 女人高潮潮喷娇喘18禁视频| 国产欧美日韩综合在线一区二区| 亚洲三区欧美一区| 亚洲精华国产精华精| 亚洲国产欧美网| 18禁国产床啪视频网站| 国产精品欧美亚洲77777| av在线app专区| 国产色视频综合| 欧美人与性动交α欧美软件| 国产男女超爽视频在线观看| 久久国产亚洲av麻豆专区| 久久久精品免费免费高清| 国产麻豆69| 天堂俺去俺来也www色官网| 男女免费视频国产| 亚洲欧美清纯卡通| 日韩视频一区二区在线观看| 日日摸夜夜添夜夜添小说| 久久精品国产亚洲av香蕉五月 | 最黄视频免费看| 国产精品99久久99久久久不卡| 午夜福利在线免费观看网站| 国产av国产精品国产| 午夜福利一区二区在线看| av片东京热男人的天堂| 老熟女久久久| 精品熟女少妇八av免费久了| 婷婷色av中文字幕| 国产黄频视频在线观看| 亚洲五月色婷婷综合| 日韩欧美免费精品| 9191精品国产免费久久| www.av在线官网国产| 精品一区在线观看国产| 日韩欧美一区视频在线观看| 十分钟在线观看高清视频www| 每晚都被弄得嗷嗷叫到高潮| 成年人黄色毛片网站| 丝袜脚勾引网站| 国产av精品麻豆| 亚洲色图综合在线观看| 国产91精品成人一区二区三区 | 丝袜脚勾引网站| 国产精品av久久久久免费| 一区二区三区四区激情视频| 久久久久精品人妻al黑| 50天的宝宝边吃奶边哭怎么回事| 1024香蕉在线观看| 久久久国产一区二区| 国产日韩欧美亚洲二区| 日本vs欧美在线观看视频| 每晚都被弄得嗷嗷叫到高潮| 男人操女人黄网站| 亚洲精品成人av观看孕妇| 欧美xxⅹ黑人| 一本大道久久a久久精品| 国产在线免费精品| 亚洲 欧美一区二区三区| 少妇被粗大的猛进出69影院| svipshipincom国产片| 国产又爽黄色视频| 久久久久视频综合| 国产精品麻豆人妻色哟哟久久| 国产极品粉嫩免费观看在线| 狠狠精品人妻久久久久久综合| 成人手机av| 亚洲专区国产一区二区| 亚洲精华国产精华精| 在线十欧美十亚洲十日本专区| av有码第一页| 99国产精品99久久久久| 视频区欧美日本亚洲| 亚洲精品久久午夜乱码| 国产一区二区三区综合在线观看| 久久99热这里只频精品6学生| 亚洲成人国产一区在线观看| 久久影院123| avwww免费| 亚洲第一欧美日韩一区二区三区 | 久久精品aⅴ一区二区三区四区| 精品第一国产精品| 国产91精品成人一区二区三区 | 又紧又爽又黄一区二区| 美女脱内裤让男人舔精品视频| 亚洲欧洲日产国产| 狠狠精品人妻久久久久久综合| 国产黄频视频在线观看| 国产亚洲一区二区精品| 国产黄频视频在线观看| 999久久久精品免费观看国产| 国产黄色免费在线视频| 欧美成人午夜精品| 亚洲五月婷婷丁香| 性色av一级| 国产成人影院久久av| 国产精品免费视频内射| 久久久久国内视频| 亚洲欧美激情在线| 欧美精品一区二区免费开放| 亚洲精品在线美女| 后天国语完整版免费观看| 欧美成人午夜精品| 青草久久国产| 超碰97精品在线观看| 嫁个100分男人电影在线观看| 侵犯人妻中文字幕一二三四区| 国产成人精品久久二区二区91| 亚洲va日本ⅴa欧美va伊人久久 | 99久久99久久久精品蜜桃| 人人妻人人爽人人添夜夜欢视频| 黄片播放在线免费| 老汉色av国产亚洲站长工具| 久久 成人 亚洲| 日韩,欧美,国产一区二区三区| 久久久久久久精品精品| www.av在线官网国产| 日日夜夜操网爽| 国产欧美日韩精品亚洲av| 黄色a级毛片大全视频| 深夜精品福利| 日日夜夜操网爽| 成年人免费黄色播放视频| 亚洲精品久久午夜乱码| 久久毛片免费看一区二区三区| 国产麻豆69| 99热全是精品| 免费看十八禁软件| 亚洲中文字幕日韩| tocl精华| 伊人久久大香线蕉亚洲五| 操出白浆在线播放| www.自偷自拍.com| 夜夜夜夜夜久久久久| 亚洲少妇的诱惑av| av又黄又爽大尺度在线免费看| tube8黄色片| 少妇的丰满在线观看| 日本91视频免费播放| 18禁裸乳无遮挡动漫免费视频| 欧美av亚洲av综合av国产av| 人妻久久中文字幕网| 日韩 欧美 亚洲 中文字幕| 亚洲综合色网址| 日韩视频在线欧美| www.熟女人妻精品国产| 亚洲男人天堂网一区| 色视频在线一区二区三区| 十八禁网站免费在线| 国产99久久九九免费精品| 欧美激情极品国产一区二区三区| 9191精品国产免费久久| 欧美xxⅹ黑人| 亚洲专区国产一区二区| 日本av免费视频播放| 黄色 视频免费看| 精品国产乱码久久久久久小说| 亚洲全国av大片| 精品一区二区三区四区五区乱码| 午夜福利视频在线观看免费| 菩萨蛮人人尽说江南好唐韦庄| 在线观看人妻少妇| www.精华液| 国产又色又爽无遮挡免| 久久99热这里只频精品6学生| 亚洲精品国产精品久久久不卡| 在线观看免费高清a一片| 国产亚洲欧美在线一区二区| 天天躁日日躁夜夜躁夜夜| 肉色欧美久久久久久久蜜桃| 下体分泌物呈黄色| 国产精品偷伦视频观看了| 啦啦啦啦在线视频资源| 99香蕉大伊视频| 成年av动漫网址| 人人妻人人澡人人看| 精品亚洲成a人片在线观看| 日韩中文字幕视频在线看片| 欧美变态另类bdsm刘玥| 黄色a级毛片大全视频| 岛国毛片在线播放| 无限看片的www在线观看| 国产成人精品久久二区二区91| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲一区中文字幕在线| 日韩大码丰满熟妇| 久久久久久人人人人人| 国产一区二区在线观看av| 欧美av亚洲av综合av国产av| 男女高潮啪啪啪动态图| 亚洲精品美女久久久久99蜜臀| 一区二区三区四区激情视频| 亚洲九九香蕉| 九色亚洲精品在线播放| 制服诱惑二区| 99久久人妻综合| 国产97色在线日韩免费| 久久精品aⅴ一区二区三区四区| 成人18禁高潮啪啪吃奶动态图| 国产在线视频一区二区| 咕卡用的链子| 久久久久久久久久久久大奶| 欧美日韩黄片免| 一进一出抽搐动态| 欧美av亚洲av综合av国产av| 亚洲精品国产色婷婷电影| 99国产精品免费福利视频| 91精品伊人久久大香线蕉| 中文字幕另类日韩欧美亚洲嫩草| 久热爱精品视频在线9| 香蕉国产在线看| 男人舔女人的私密视频| 精品久久蜜臀av无| 亚洲精品一二三| 国产av又大| 色婷婷av一区二区三区视频| 999精品在线视频| 天天操日日干夜夜撸| 一本大道久久a久久精品| 成人免费观看视频高清| 淫妇啪啪啪对白视频 | 欧美日韩成人在线一区二区| 精品福利观看| 亚洲七黄色美女视频| 国产成人免费无遮挡视频| 深夜精品福利| 久久久国产一区二区| 乱人伦中国视频| 最新的欧美精品一区二区| 亚洲成人免费电影在线观看| 中文欧美无线码| 悠悠久久av| 国产av又大| 丰满饥渴人妻一区二区三| 欧美日韩亚洲综合一区二区三区_| 亚洲精品国产色婷婷电影| 久久亚洲精品不卡| 成年女人毛片免费观看观看9 | 永久免费av网站大全| av不卡在线播放| 中文字幕另类日韩欧美亚洲嫩草| 99久久99久久久精品蜜桃| 久久久国产精品麻豆| 一级毛片精品| 老汉色av国产亚洲站长工具| 欧美久久黑人一区二区| 国产91精品成人一区二区三区 | 日本av手机在线免费观看| 欧美激情 高清一区二区三区| 男女国产视频网站| 两个人免费观看高清视频| 多毛熟女@视频| 一区二区日韩欧美中文字幕| 大片电影免费在线观看免费| 成年动漫av网址| 丰满人妻熟妇乱又伦精品不卡| 亚洲国产日韩一区二区| 欧美人与性动交α欧美软件| 国产精品偷伦视频观看了| 乱人伦中国视频| 法律面前人人平等表现在哪些方面 | 久久人人爽av亚洲精品天堂| 老司机影院毛片| a 毛片基地| 久9热在线精品视频| 午夜福利视频在线观看免费| 国产日韩一区二区三区精品不卡| 成在线人永久免费视频| 精品少妇一区二区三区视频日本电影| 久久久精品国产亚洲av高清涩受| 久久精品aⅴ一区二区三区四区| av在线播放精品| 亚洲精品中文字幕一二三四区 | 在线精品无人区一区二区三| 又黄又粗又硬又大视频| 亚洲精品自拍成人| 韩国精品一区二区三区| svipshipincom国产片| av有码第一页| 精品福利观看| 国产国语露脸激情在线看| 日韩电影二区| 丁香六月欧美| 人成视频在线观看免费观看| 国产欧美日韩一区二区三区在线| 欧美另类亚洲清纯唯美| 欧美激情极品国产一区二区三区| 亚洲国产精品成人久久小说| 欧美日韩一级在线毛片| 成年动漫av网址| 国产成人免费无遮挡视频| 欧美精品一区二区大全| 欧美 日韩 精品 国产| 亚洲综合色网址| 欧美日本中文国产一区发布| 国产精品熟女久久久久浪| 成人黄色视频免费在线看| 天天影视国产精品| 欧美少妇被猛烈插入视频| videosex国产| 亚洲五月色婷婷综合| 性少妇av在线| 免费在线观看完整版高清| 无限看片的www在线观看| 青草久久国产| 搡老岳熟女国产| 亚洲精品久久午夜乱码| 99国产精品一区二区蜜桃av | 欧美成人午夜精品| 女人精品久久久久毛片| 欧美黑人欧美精品刺激| 亚洲va日本ⅴa欧美va伊人久久 | 亚洲欧美激情在线| 在线av久久热| 国产欧美日韩精品亚洲av| 久久久久网色| 国产国语露脸激情在线看| 少妇被粗大的猛进出69影院| 亚洲欧美一区二区三区久久| 在线观看免费日韩欧美大片| 国产在线观看jvid| av国产精品久久久久影院| 国产精品久久久久久人妻精品电影 | 亚洲少妇的诱惑av| 两个人看的免费小视频| 51午夜福利影视在线观看| 悠悠久久av| 欧美中文综合在线视频| av有码第一页| 秋霞在线观看毛片| 国产精品亚洲av一区麻豆| 精品少妇黑人巨大在线播放| 手机成人av网站| 国产精品秋霞免费鲁丝片| 考比视频在线观看| 午夜91福利影院| 青草久久国产| 青春草视频在线免费观看| 国产免费视频播放在线视频| 久久久久国内视频| 欧美av亚洲av综合av国产av| 啦啦啦在线免费观看视频4| 老鸭窝网址在线观看| 女人高潮潮喷娇喘18禁视频| 淫妇啪啪啪对白视频 | 美国免费a级毛片| 男人操女人黄网站| 午夜免费观看性视频| 在线观看舔阴道视频| 两性夫妻黄色片| avwww免费| 99久久99久久久精品蜜桃| 久久性视频一级片| 日本欧美视频一区| 亚洲中文日韩欧美视频| 大香蕉久久成人网| 美女中出高潮动态图| 大片免费播放器 马上看| 久久性视频一级片| 国产精品香港三级国产av潘金莲| 精品卡一卡二卡四卡免费| 亚洲美女黄色视频免费看| 国产淫语在线视频| 午夜91福利影院| 免费观看a级毛片全部| 考比视频在线观看| 老司机深夜福利视频在线观看 | 精品亚洲成国产av| 国产男女内射视频| 深夜精品福利| 国产亚洲精品一区二区www | 欧美日本中文国产一区发布| av免费在线观看网站| 18禁裸乳无遮挡动漫免费视频| 欧美激情极品国产一区二区三区| 久久人人爽人人片av| 国产男人的电影天堂91| 老司机影院毛片|