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

    天津大氣擴(kuò)散條件對污染物垂直分布的影響研究

    2022-04-24 11:09:06劉敬樂蔡子穎韓素芹
    中國環(huán)境科學(xué) 2022年4期
    關(guān)鍵詞:穩(wěn)定度湍流大氣

    劉敬樂,史 靜,姚 青,蔡子穎*,韓素芹,姜 明,崔 曄

    天津大氣擴(kuò)散條件對污染物垂直分布的影響研究

    劉敬樂1,2,4,史 靜1,姚 青3,4,蔡子穎3,4*,韓素芹3,姜 明1,崔 曄1

    (1.天津市氣象探測中心,天津 300074;2.天津市海洋氣象重點(diǎn)實(shí)驗(yàn)室,天津 300074;3.天津市環(huán)境氣象中心,天津 300074;4.天津市氣象科學(xué)研究所,天津 300074)

    利用2017~2019年夏、冬季天津市大氣污染物監(jiān)測和氣象觀測數(shù)據(jù),基于天津氣象鐵塔垂直觀測,針對大氣垂直擴(kuò)散條件對PM2.5和O3的影響進(jìn)行研究.結(jié)果顯示:近地面PM2.5濃度隨高度的升高而下降,O3濃度則隨高度的升高而上升,受大氣垂直擴(kuò)散條件的季節(jié)和日變化影響,冬季,地面與120m PM2.5質(zhì)量濃度相關(guān)明顯,與200m PM2.5質(zhì)量濃度無明顯相關(guān).夏季,120m和200m PM2.5質(zhì)量濃度相關(guān)系數(shù)為0.72,午后通常出現(xiàn)120m和200m PM2.5質(zhì)量濃度高于地面的情況.夏季,不同高度O3濃度差異小于冬季,地面與120m高度O3濃度接近.以大氣穩(wěn)定度、逆溫強(qiáng)度和氣溫遞減率作為大氣垂直擴(kuò)散指標(biāo),對地面PM2.5和O3垂直分布具有指示作用.冬季,TKE與PM2.5質(zhì)量濃度相關(guān)系數(shù)為到-0.65,夏季,TKE與ΔPM2.5相關(guān)系數(shù)為-0.39.夏、冬季TKE與地面O3濃度的相關(guān)系數(shù)分別為0.46和0.53,與ΔO3的相關(guān)系數(shù)分別為0.73和0.70.弱下沉運(yùn)動對地面O3濃度影響較強(qiáng),40m高度垂直運(yùn)動速度與地面O3濃度的相關(guān)系數(shù)在冬、夏季分別為-0.54和-0.61.對冬季典型PM2.5重污染過程的分析發(fā)現(xiàn),霧霾的生消維持和PM2.5濃度的變化與大氣穩(wěn)定度、氣溫垂直遞減率和TKE的變化有直接關(guān)系.對夏季典型O3污染過程的分析發(fā)現(xiàn),近地面的O3污染的形成與有利光化學(xué)反應(yīng)的氣象條件密切相關(guān),同時,垂直向下輸送和有利垂直擴(kuò)散條件對O3污染的形成和爆發(fā)影響明顯.

    PM2.5;O3;氣象塔;大氣垂直擴(kuò)散條件;天津

    PM2.5和O3污染是我國主要的大氣環(huán)境問題,近年來在PM2.5質(zhì)量濃度總體下降的趨勢背景下,對流層的O3濃度明顯上升,受到廣泛的關(guān)注[1-2].污染發(fā)生的內(nèi)因是排放量大,外因是氣象條件,氣象條件可通過影響通風(fēng)率、降水沉降、干沉積、化學(xué)轉(zhuǎn)化損失率、自然排放量和背景濃度等影響大氣氣溶膠濃度和臭氧濃度.與人類活動關(guān)系最為密切的大氣邊界層結(jié)構(gòu)不僅直接影響一次污染物的輸送、積累和擴(kuò)散,溫濕條件以及大氣氧化性等邊界層內(nèi)理化特性也會對二次氣溶膠的生成和轉(zhuǎn)化產(chǎn)生重要影響[3],城市邊界層特別是城市近地層的風(fēng)、溫、濕和湍流結(jié)構(gòu)存在明顯的時空變化[4],對城市大量污染物的輸送、擴(kuò)散、化學(xué)轉(zhuǎn)化和干濕清除起著重要作用[5-6].

    天津位于華北平原北部,東臨渤海,地形和超大城市下墊面造成了其邊界層結(jié)構(gòu)具有獨(dú)特的特點(diǎn).隨著天津城市化進(jìn)程的加快以及濱海新區(qū)石化、化工、冶金等產(chǎn)業(yè)規(guī)模的擴(kuò)大,不可避免地帶來高排放的人為污染物,并在適宜的氣象條件下形成重污染天氣[7].以往針對天津地區(qū)重污染天氣中邊界層的污染物濃度和化學(xué)組分[8]污染形成機(jī)制[9]以及利用地基遙感觀測、氣象塔和無人機(jī)等方式開展的地面污染物和氣象要素垂直分布特征[10-11]的研究已有諸多研究成果.其中,利用氣象塔進(jìn)行氣象要素觀測和大氣污染物的監(jiān)測具有較高的觀測分辨率和時間連續(xù)性,是最為可靠的一種方式,以往利用氣象塔開展針對典型霧霾天氣過程中溫濕和湍流特征等精細(xì)結(jié)構(gòu)的變化規(guī)律研究[12],大氣穩(wěn)定度、逆溫等變化特征對PM2.5污染的影響分析[13],O3、NOx、CO、氣溶膠等污染物的梯度觀測及變化規(guī)律研究[14-15],湍流強(qiáng)度等微氣象要素對氣溶膠和氣態(tài)污染物的影響研究[16-17],發(fā)現(xiàn)邊界層高度[18]、大氣穩(wěn)定度[19]、理查森數(shù)[20]、溫度層結(jié)[21]和湍流動能[22]等表征大氣垂直擴(kuò)散能力的指標(biāo)對重污染天氣的形成發(fā)展具有指示作用.例如,當(dāng)天津地區(qū)日均氣溫遞減率小于0.4℃/100m時,垂直擴(kuò)散條件不利于大氣污染物擴(kuò)散,中度以上污染出現(xiàn)概率明顯增加[23].天津地區(qū)的湍流擴(kuò)散系數(shù)與PM2.5質(zhì)量濃度相關(guān)系數(shù)可以達(dá)到-0.56[24].為更進(jìn)一步深入了解天津大氣垂直擴(kuò)散條件對重污染發(fā)生、維持和消散過程的作用機(jī)理,本研究依托天津氣象鐵塔(以下簡稱氣象塔),通過對2017~2019年天津主要污染物PM2.5和O3的立體監(jiān)測數(shù)據(jù)以及同期氣象觀測資料的分析,旨在揭示湍流強(qiáng)度、大氣穩(wěn)定度、溫度層結(jié)等近地面大氣垂直擴(kuò)散條件對污染物的影響規(guī)律,從微觀角度進(jìn)一步增進(jìn)對天津地區(qū)邊界層結(jié)構(gòu)和污染物垂直分布特征的深入了解,為建立重污染預(yù)報預(yù)警中大氣邊界層的識別條件因子提供支撐.

    1 數(shù)據(jù)和方法

    氣象和空氣質(zhì)量監(jiān)測數(shù)據(jù)來自2017~2019年冬(11月~次年2月)、夏季(6~8月)開展的綜合立體觀測,獲得的包括天津13個國家級自動氣象站和26個環(huán)保局監(jiān)測站的觀測數(shù)據(jù)以及氣象塔的相關(guān)觀測數(shù)據(jù),站點(diǎn)位置如圖1所示.天津氣象塔塔高255m,位于中國氣象局天津大氣邊界層觀測站內(nèi)(39°04′N, 117°12′E,海拔2.2m),屬于商業(yè)、居民混合點(diǎn)位,其周圍環(huán)境基本代表了大城市復(fù)雜下墊面,氣象塔5~250m安裝15層風(fēng)、溫、濕度傳感器,可實(shí)現(xiàn)風(fēng)速、風(fēng)向、溫度和相對濕度的實(shí)時在線觀測,氣象塔40m高度平臺安裝有超聲風(fēng)速儀和H2O/CO2觀測儀組成的渦動協(xié)方差觀測系統(tǒng),用以測量大氣下墊面界面上的物質(zhì)和能量交換,并能夠從宏觀和微觀氣象場上對污染過程進(jìn)行解析.氣象塔地面、120m和200m安裝有3套空氣質(zhì)量監(jiān)測站,包括PM10、PM2.5、SO2、NO2和O3濃度監(jiān)測儀器,可實(shí)現(xiàn)塔層內(nèi)大氣污染物的立體觀測.本文應(yīng)用的主要觀測儀器參數(shù)如表1所示.

    研究中的湍流動能(TKE)基于氣象塔超聲風(fēng)速儀,在獲取湍流數(shù)據(jù)基礎(chǔ)上,基于德國拜耳羅大學(xué)研發(fā)的TK3軟件計(jì)算獲得,該軟件包含了側(cè)向風(fēng)效應(yīng)校正、野點(diǎn)值去除、坐標(biāo)旋轉(zhuǎn)處理、WPL修正等步驟,以及湍流平穩(wěn)性檢驗(yàn)、近地層湍流譜特征分析、湍流發(fā)展充分性檢驗(yàn)和通量觀測的空間代表性分析等數(shù)據(jù)質(zhì)量評價方法,湍流計(jì)算時間長度選擇30min.平均湍流動能的計(jì)算公式如下:

    式中、和分別為三維超聲風(fēng)速脈動量,TKE的單位為m2/s2.

    研究中的大氣穩(wěn)定度的計(jì)算采用了溫差—風(fēng)速法[24],此方法能較好的反映大氣熱力湍流和機(jī)械湍流的影響.對大氣穩(wěn)定度分類標(biāo)準(zhǔn)如表2所示,大氣穩(wěn)定度A表示強(qiáng)不穩(wěn)定,B表示不穩(wěn)定,C表示弱不穩(wěn)定,D表示中性,E表示較穩(wěn)定,F表示穩(wěn)定.以往研究結(jié)果證明[11],對于天津地區(qū),此種大氣穩(wěn)定度分類標(biāo)準(zhǔn)較為合理.

    圖1 觀測站點(diǎn)位置示意

    表1 觀測儀器參數(shù)

    表2 △T和u組合法確定穩(wěn)定度的分類標(biāo)準(zhǔn)

    注:表中代表40m高度風(fēng)速,△代表100m和30m高度溫差.

    2 結(jié)果與討論

    2.1 觀測期間天津市空氣質(zhì)量特征

    根據(jù)觀測期間的天津生態(tài)環(huán)境局監(jiān)測空氣質(zhì)量數(shù)據(jù)顯示(圖2),天津地區(qū)冬季多以PM2.5為首要污染物.2017~2019年冬季PM2.5質(zhì)量濃度分別為60,75和70μg/m3,其中出現(xiàn)中度以上污染天氣日數(shù)分別為15,25和24d.2017~2019年冬季PM10質(zhì)量濃度分別為92,110和89μg/m3,其中污染物超標(biāo)日數(shù)分別為14,24和22d,均為輕度或中度污染水平.近3年冬季NO2質(zhì)量濃度分別為51,54和52μg/m3,污染物超標(biāo)日數(shù)分別為10,19和11d,超標(biāo)日均為輕度污染水平.近3年冬季中O3和SO2污染相對較低,其空氣質(zhì)量等級均達(dá)到優(yōu)良水平.

    天津地區(qū)夏季多以O(shè)3為首要污染物.2018年和2019年夏季中 O3日最大8h質(zhì)量濃度滑動平均值(MDA8)分別為165和155μg/m3,污染物超標(biāo)日數(shù)分別為45和41d,其中達(dá)到中度以上污染分別為13和12d.夏季的PM2.5和PM10質(zhì)量濃度明顯低于冬季,近2年夏季PM2.5質(zhì)量濃度分別為42和38μg/m3,超標(biāo)日數(shù)分別為2和4d;PM10質(zhì)量濃度分別為68和60μg/m3.夏季NO2和SO2污染相對較低,其空氣質(zhì)量等級均達(dá)到優(yōu)良水平.

    圖2 觀測期間空氣質(zhì)量等級日數(shù)

    a: 2017年冬季;b: 2018年夏季;c: 2018年冬季;d: 2019年夏季;e:2019年冬季

    2.2 觀測期間氣象塔PM2.5和O3質(zhì)量濃度垂直分布特征

    冬季觀測期間,地面、120m和200m高度的PM2.5質(zhì)量濃度分別為74,63和42 μg/m3.受邊界層結(jié)構(gòu)以及污染來源差異的影響,在不同空氣質(zhì)量條件下,近地面PM2.5質(zhì)量濃度的垂直分布存在差異(圖3),當(dāng)空氣質(zhì)量為優(yōu)良時,較好的垂直擴(kuò)散條件使各高度PM2.5質(zhì)量濃度接近.當(dāng)出現(xiàn)污染天氣時,通常伴隨大氣垂直擴(kuò)散條件減弱,因此各層PM2.5質(zhì)量濃度差異增大,輕、中度污染時,地面較120m和200m高度的PM2.5質(zhì)量濃度高約為17%~35%和40%~64%,當(dāng)出現(xiàn)重度污染時,地面較120m和200m高度的PM2.5質(zhì)量濃度高約為25%~42%和40%~62%.可見,120m高度與地面的PM2.5質(zhì)量濃度差異與地面污染程度明顯相關(guān),但200m高度由于超過常通量層,湍流輸送下降,且受逆溫、霧頂?shù)冗吔鐚哟怪苯Y(jié)構(gòu)的影響[10],PM2.5質(zhì)量濃度與地面污染程度無明顯相關(guān),這意味著冬季中120m到200m之間污染物無法充分的混合,導(dǎo)致大氣污染物更多的被滯留在近地層內(nèi)無法向上擴(kuò)散.夏季觀測期間,地面、120m和200m高度的PM2.5質(zhì)量濃度分別為41,34和33 μg/m3,由于大氣垂直擴(kuò)散條件普遍優(yōu)于冬季,因此PM2.5垂直混合均勻,120m和200m高度的PM2.5質(zhì)量濃度變化較為一致,兩者相關(guān)系數(shù)達(dá)到0.72**(**表示通過99%置信度檢驗(yàn)).

    受大氣垂直擴(kuò)散條件變化的影響,不同高度PM2.5質(zhì)量濃度日變化也具有明顯差異.冬季(圖4),受城市早晚交通高峰排放源影響,地面和120m高度的PM2.5質(zhì)量濃度均呈現(xiàn)明顯的早晚高峰,午后波谷的特征,200m高度的PM2.5質(zhì)量濃度日變化則較為平緩,主要是由于這一高度湍流輸送作用明顯減弱,因此局地污染對其影響較低.同時,由于熱力湍流作用對120m高度的影響較地面滯后,因此120m高度PM2.5質(zhì)量濃度峰值出現(xiàn)較地面滯后約2h.夏季,由于午后大氣湍流運(yùn)動較強(qiáng),各高度PM2.5質(zhì)量濃度差異明顯小于冬季,甚至可出現(xiàn)120m和200m的PM2.5濃度高于地面的情況,在交通晚高峰階段,熱力湍流作用將更多的PM2.5從地面帶向城市上空使地面PM2.5質(zhì)量濃度仍維持較低,但并不意味大氣中PM2.5含量減少[3].

    圖3 冬季不同空氣質(zhì)量條件下PM2.5質(zhì)量濃度垂直分布

    箱型圖表示10至90分位范圍和中位值,廓線表示平均值

    由于對流層內(nèi)的O3主要來自對流層的光化學(xué)反應(yīng)生成和平流層內(nèi)O3的向下輸送,因此,O3質(zhì)量濃度的垂直分布特征與PM2.5差異明顯,在大氣邊界層內(nèi)通常呈現(xiàn)濃度隨高度上升的特征?[20].冬季,地面、120m和200m高度的O3(MDA8)值分別為10,19和25μg/m3.天津地區(qū)O3污染日主要出現(xiàn)在5~9月[25],夏季,地面、120m和200m高度的O3(MDA8)值分別為153、151和178μg/m3.各高度O3質(zhì)量濃度日變化(圖4)顯示,由于近地面O3濃度主要受光化學(xué)反應(yīng)和NO、VOCs等前體物質(zhì)濃度影響[26],因此O3質(zhì)量濃度日變化有明顯的單峰型特征,日出后隨著光化學(xué)反應(yīng)強(qiáng)度增長和夜間殘留層內(nèi)O3向下垂直混合作用的影響,近地面O3質(zhì)量濃度迅速上升,并在中午到達(dá)一天中的峰值,午后隨著光化學(xué)反應(yīng)的減弱及地面排放的NO對O3的滴定作用增強(qiáng),O3質(zhì)量濃度逐步下降,夜間由于區(qū)域輸送作用很弱同時沒有光化學(xué)反應(yīng)生成,因此各高度O3濃度通常較低且保持相對穩(wěn)定.受地面排放源影響,地面的O3質(zhì)量濃度日變化最為顯著,隨著高度的升高,地面排放源影響逐漸減小,日變化趨于平緩.冬季,湍流運(yùn)動相對較強(qiáng)的午間時段,不同高度O3質(zhì)量濃度相對接近,地面與120m和200m高度O3濃度比分別為0.76和0.71.午后,高層O3質(zhì)量濃度下降速度明顯慢于低層,夜間也維持相對較高的濃度,造成了夜間O3濃度垂直差異增大,這是由于夜間湍流運(yùn)動減弱,地面排放的NO等還原物質(zhì)較難向上輸送,對高層的O3消耗量遠(yuǎn)小于低層[20].夏季,各高度O3質(zhì)量濃度差異低于冬季,與夏季大氣垂直擴(kuò)散條件普遍優(yōu)于冬季有關(guān),日間,地面與120m高度的O3質(zhì)量濃度接近,在9:00~14:00,地面O3質(zhì)量濃度甚至高于120m,夜間濃度差異雖增加,但仍明顯低于冬季同時段中的差異.由此可見,大氣垂直擴(kuò)散條件對O3質(zhì)量濃度的垂直差異具有顯著的影響.

    2.3 邊界層垂直擴(kuò)散條件對大氣污染物濃度的影響

    大氣污染物的垂直擴(kuò)散能力可由湍流強(qiáng)度表征,湍流動能(TKE)是湍流強(qiáng)度的量度[27]分析觀測期間PM2.5質(zhì)量濃度與TKE的相關(guān)性發(fā)現(xiàn)(為降低風(fēng)速和風(fēng)向?qū)M2.5水平擴(kuò)散的影響,僅分析水平風(fēng)速小于1m/s的條件),冬季,兩者相關(guān)系數(shù)達(dá)到-0.65**(圖5),夏季,由于PM2.5質(zhì)量濃度普遍較低,因此兩者相關(guān)系數(shù)僅為-0.17,但TKE與ΔPM2.5(地面與200m高度濃度差值)的相關(guān)系數(shù)仍可達(dá)到-0.39**.以往針對天津地區(qū)TKE對PM2.5影響的研究中指出[14],當(dāng)TKE為0~0.3m2/s2時,湍流受到抑制不利于污染物擴(kuò)散;當(dāng)TKE為0.3~0.5m2/s2時,湍流較弱,大氣污染物日變化不顯著;當(dāng)TKE為0.5~0.8m2/s2時,湍流偏強(qiáng),污染物濃度有明顯的日變化,不易持續(xù)累積;當(dāng)TKE為0.8~2.0m2/s2時,湍流較強(qiáng),垂直擴(kuò)散條件非常利于污染物擴(kuò)散;當(dāng)TKE大于2.0m2/s2時,通常有大風(fēng)或強(qiáng)對流天氣過程,空氣質(zhì)量為優(yōu)良等級.2018年冬季觀測期間出現(xiàn)的污染過程中,當(dāng)TKE在上述5個區(qū)間時PM2.5質(zhì)量濃度分別為(172±50) μg/m3,(169±58)μg/m3,(122±41)μg/m3,(70±31)μg/m3和(26±18)μg/m3.可見在水平風(fēng)場影響較小的情況下,TKE對PM2.5質(zhì)量濃度及其垂直分布影響顯著.

    圖5 天津冬季PM2.5和夏季O3質(zhì)量濃度與TKE分布關(guān)系

    表3 冬季不同大氣穩(wěn)定度等級下氣溫遞減率、TKE和PM2.5濃度

    表4 夏季不同大氣穩(wěn)定度等級下氣溫遞減率、TKE和O3濃度

    近地面大氣熱力層結(jié)對污染的影響作用同樣明顯[13],以往研究顯示,大氣層結(jié)穩(wěn)定條件下多出現(xiàn)貼地逆溫層,且逆溫強(qiáng)度和逆溫層厚度較高?[28],當(dāng)夜間大氣穩(wěn)定度為E以上以及白天大氣穩(wěn)定度為D及以上時,PM2.5質(zhì)量濃度會顯著高于平均值[24].冬季觀測期間,逆溫出現(xiàn)概率為40%,平均逆溫強(qiáng)度為0.95℃/100m.當(dāng)大氣穩(wěn)定度為A~C時未出現(xiàn)逆溫,隨著大氣穩(wěn)定度增長,逆溫出現(xiàn)概率和強(qiáng)度明顯增加(表3).同時,根據(jù)湍流促發(fā)條件,當(dāng)出現(xiàn)位溫逆溫時湍流的發(fā)生也將受到抑制[14],如表3所示,隨著大氣穩(wěn)定度的增加,氣溫遞減率和TKE均有明顯下降,同時地面PM2.5濃度和ΔPM2.5有明顯升高,當(dāng)大氣穩(wěn)定度為D時,出現(xiàn)輕度和中度污染的概率接近50%;當(dāng)大氣穩(wěn)定度為E時,出現(xiàn)中度以上污染的概率為67%;當(dāng)大氣穩(wěn)定度為F時,出現(xiàn)中度以上污染的概率為85%. TKE對O3質(zhì)量濃度同樣具有指示作用,在O3主要生成的白天時段內(nèi)(8:00~20:00)兩者相關(guān)系數(shù)在夏季和冬季分別為0.46**(圖5)和0.53**. TKE對ΔO3(200m與地面濃度差值)的影響更為明顯,其相關(guān)系數(shù)在夏季和冬季分別達(dá)到了0.73**和0.70**.通過對夏季TKE和O3質(zhì)量濃度的聚類分析發(fā)現(xiàn),當(dāng)TKE為0~0.5m2/s2時,湍流受到抑制,殘留層內(nèi)O3的向下傳輸和地面O3前體物的向上輸送能力弱;當(dāng)TKE為0.5~1.0m2/s2時,湍流較弱,O3質(zhì)量濃度的垂直差異較大;當(dāng)TKE為1.0~3.0m2/s2時,湍流較強(qiáng),利于殘留層內(nèi)O3的向下傳輸和地面O3前體物的向上輸送能力,O3質(zhì)量濃度垂直差異相對較小;當(dāng)TKE大于3.0m2/s2時,湍流強(qiáng)盛,通常伴有對流天氣過程出現(xiàn)且近地面風(fēng)速大,O3的水平和垂直擴(kuò)散較快.2018年夏季觀測期間,當(dāng)TKE在以上4個區(qū)間時,地面O3質(zhì)量濃度分別為(57±42)μg/m3,(72±70)μg/m3,(115±80)μg/ m3和(88±50)μg/m3.夏季觀測期間,大氣層結(jié)穩(wěn)定度對地面O3濃度和DO3具有一定的影響,如表4所示,隨著大氣穩(wěn)定度的增加,地面O3濃度和DO3明顯降低,當(dāng)大氣穩(wěn)定度等級為A和B時,出現(xiàn)O3污染天氣概率為37%;當(dāng)大氣穩(wěn)定度等級為C和D時,出現(xiàn)O3污染天氣概率為12%;當(dāng)大氣穩(wěn)定度等級為E和F時,未出現(xiàn)O3污染.可見,湍流混合作用對近地面O3濃度有著重要的影響,同時,這也是近年來大氣垂直擴(kuò)散能力改善后,華北平原O3濃度持續(xù)升高的原因之一[25].地面O3濃度的增加與近地面O3的垂直向下輸送運(yùn)動直接相關(guān)[9],通過分析40m高度的垂直運(yùn)動速度與地面O3質(zhì)量濃度的具有明顯相關(guān)性,冬、夏季的相關(guān)系數(shù)分別為-0.54**和-0.61**.夏季,當(dāng)垂直風(fēng)速大于0m/s,即為上升運(yùn)動時,O3平均濃度為(60±29)μg/m3;當(dāng)垂直速度為0~-0.1m/s時,O3平均濃度為(87±36)μg/m3;當(dāng)垂直速度為-0.1~-0.3m/s時,O3平均濃度為(121±45)μg/m3;當(dāng)垂直速度小于-0.3m/s時,O3平均濃度為(101±30)μg/m3. 2018年夏季,40m高度的平均垂直速度為-0.06? m/s,但是在O3中度污染日的平均垂直速度為-0.15m/s.可見,弱下沉運(yùn)動對近地面O3濃度增加的影響作用明顯.

    2.4 觀測期間典型PM2.5和O3污染過程分析

    2.4.1 典型PM2.5重污染過程分析 2019年1月10~14日,天津地區(qū)出現(xiàn)一次連續(xù)四日AQI超過200的重污染天氣過程.污染過程期間(10日7:00~14日23:00),首要污染物PM2.5質(zhì)量濃度均值達(dá)到213μg/ m3,其中重污染時間占污染總時長的84%.從天氣形勢場分析,1月10日,我國中東部地區(qū)高空受弱高壓控制,850hPa以偏西氣流為主,京津冀中南部處于低壓槽區(qū),天氣形勢穩(wěn)定,伴隨西南風(fēng)的輸送和本地不利的氣象條件,污染物快速積累,10日20:00空氣質(zhì)量達(dá)到重污染水平.11日夜間,我國中東部地區(qū)出現(xiàn)大范圍霧區(qū),12日03:00天津開始出現(xiàn)霧天氣,并維持至13日09:00,其間由于霧頂?shù)哪鏈匾约靶★L(fēng)、高濕度環(huán)境,使細(xì)顆粒物快速累積和吸濕增長,PM2.5質(zhì)量濃度從242μg/m3升至310μg/m3,空氣質(zhì)量達(dá)到嚴(yán)重污染水平.13日午間,天津受弱冷空氣影響,空氣質(zhì)量有所好轉(zhuǎn),但PM2.5質(zhì)量濃度仍維持在200μg/m3以上.12日下半夜至13日上午,一股弱冷空氣開始影響京津冀北部地區(qū),受鋒前輸送影響,空氣質(zhì)量再度達(dá)到嚴(yán)重污染水平,PM2.5質(zhì)量濃度升至330μg/m3.13日午間,受冷空氣影響,大霧瓦解,相對濕度迅速降低,能見度有所回升,PM2.5濃度略有下降.13日夜間,天津地區(qū)轉(zhuǎn)為弱高壓場控制,隨著風(fēng)力減弱,PM2.5質(zhì)量濃度持續(xù)升高,天津再次出現(xiàn)嚴(yán)重污染天氣.14日夜間隨著冷空氣侵入天津,空氣質(zhì)量迅速好轉(zhuǎn),污染過程隨之結(jié)束.

    從此次污染過程的氣象條件和垂直擴(kuò)散條件可以看到(圖6),10日白天,大氣由中性層結(jié)轉(zhuǎn)向穩(wěn)定,大氣穩(wěn)定度以C和D為主,TKE為0.5~1.0m2/s2, PM2.5濃度開始逐步升高,10日夜間,由于貼地逆溫層的形成,大氣垂直擴(kuò)散條件下降明顯,TKE大部分時間在0~0.3m2/s2,大氣穩(wěn)定度達(dá)到E,至11日03:00,近地面氣溫垂直遞減率已下降至-0.94℃/100m,大氣垂直交換能力下降到一個極弱的水平,PM2.5質(zhì)量濃度達(dá)到214μg/m3.11日白天,大氣垂直擴(kuò)散條件仍維持較低水平,TKE僅在12:00前后達(dá)到約0.3m2/s2,由于大氣垂直交換能力下降,地面PM2.5持續(xù)升高至250μg/m3,120m和200m高度的PM2.5質(zhì)量濃度達(dá)到了150和130μg/m3.11日夜間至13日上午,隨著霧天氣形成和維持,近地面氣溫垂直遞減率明顯下降,由于霧頂高度維持在120m左右,且霧頂存在4~6℃的強(qiáng)逆溫,因此,霧過程中120m和200m高度的PM2.5質(zhì)量濃度明顯低于地面,但由于污染物的持續(xù)積累,仍分別達(dá)到150和128μg/m3.13日午間,弱冷空氣的侵入使霧天氣過程結(jié)束,但由于近地面仍維持著強(qiáng)度約為1℃/100m的逆溫層結(jié),TKE沒有明顯躍升,因此空氣質(zhì)量雖有所好轉(zhuǎn)但仍維持在重度污染水平.13日夜間,隨著逆溫強(qiáng)度增至3℃/100m,TKE降至0.2m2/s2以下,PM2.5質(zhì)量濃度再次超過300μg/m3.14日午后,TKE升至0.5m2/s2以上,由于大氣垂直擴(kuò)散條件轉(zhuǎn)好,地面PM2.5質(zhì)量濃度有所下降,但120m和200m高度的PM2.5濃度有所升高.14日21:00前后,受冷空氣影響,逆溫層結(jié)快速瓦解,TKE升至3m2/s2以上,由于大氣水平和垂直擴(kuò)散能力的提升,空氣質(zhì)量迅速恢復(fù)至優(yōu)良水平.

    綜合而言,此次重污染過程中大氣穩(wěn)定度、氣溫垂直遞減率和TKE的變化與霧霾的生消維持和PM2.5質(zhì)量濃度的變化有較好的對應(yīng),大氣垂直擴(kuò)散條件的變化對此次霧霾生消維持中的影響作用明顯.

    2.4.2 典型O3污染過程分析 2018年6月6日和8月10日,天津出現(xiàn)了兩次O3重污染過程,以下通過對兩次過程中影響O3的主要?dú)庀笠蛩睾痛髿獯怪睌U(kuò)散條件的分析,可以清晰其在O3污染過程中的影響作用.

    6月6日當(dāng)日最高氣溫達(dá)到38.5℃,10:00~14:00平均紫外輻射強(qiáng)度為43.7W/m2,白天平均相對濕度為35%(圖7),氣象條件利于光化學(xué)反應(yīng)生成O3.同時,當(dāng)日天津地面受低壓場控制,天津地區(qū)有明顯輻合風(fēng)場,水平風(fēng)場的輻合對O3的累積有重要的影響.此次O3污染過程前期(5月31日~6月5日),天津一直處于O3中度污染,6月5日O3(MDA8)濃度已達(dá)到202?μg/m3.5日午間平均垂直速度為-0.49m/s,TKE為3.8m2/s2,大氣垂直擴(kuò)散條件對高層O3向下輸送非常有利,6日午間,垂直速度和TKE雖有所減弱,但近地面較強(qiáng)的水平輻合使O3濃度劇增,上述條件對O3的積累極其有利,最終導(dǎo)致6日O3(MDA8)濃度達(dá)到337μg/m3.7日的氣溫和濕度條件相對6日較不利于O3生成,同時,7日白天的平均垂直速度接近0m/s,TKE約為1.5m2/s2,大氣垂直擴(kuò)散條件不利于O3的向下輸送,同時隨著地面風(fēng)場輻合減弱,O3濃度較前日明顯下降.

    圖6 2019年1月10~14日天津市區(qū)地面能見度、相對濕度、風(fēng)向風(fēng)速,PM2.5質(zhì)量濃度和TKE、大氣穩(wěn)定度、氣溫垂直遞減率變化情況

    圖7 2018年6月5~7日天津市區(qū)氣溫、相對濕度、風(fēng)向風(fēng)速、紫外輻射強(qiáng)度,O3質(zhì)量濃度和TKE、垂直速度變化情況

    8月10日的O3污染過程呈現(xiàn)明顯的爆發(fā)特征,此次污染過程中地面O3(MDA8)濃度達(dá)到了299 μg/m3,為觀測期間的最高值.8月10日最高氣溫達(dá)到33.5℃,10:00~14:00平均紫外輻射強(qiáng)度為33.0W/m2,日間平均相對濕度為57%(圖8),相較于6月6日,當(dāng)日的光化學(xué)反應(yīng)氣象條件相對較差.同時,此次污染過程的O3前期積累同樣弱于6月6日,8月9日O3(MDA8)濃度為151 μg/m3,10日8時O3濃度僅為63 μg/m3.8月10日的O3濃度陡增與O3向下輸送密切相關(guān),由于8月10日天津受高壓脊控制,地面呈現(xiàn)北部弱高壓,存在明顯的下沉氣流,10日9時,隨著向下氣流的垂直速度由0.23m/s迅速升降至0.56m/s及TKE的快速增長,午后地面O3濃度出現(xiàn)了爆發(fā)性增長,17:00地面O3濃度達(dá)到374μg/m3,從天津風(fēng)廓線雷達(dá)和臭氧激光雷達(dá)上同樣可以清晰地看到此次O3重污染過程中存在明顯的下沉氣流及高空高濃度O3的向下輸送.

    綜合以上兩次O3污染過程的分析可以發(fā)現(xiàn),O3污染的形成與有利的光化學(xué)反應(yīng)氣象條件密切相關(guān),O3的向下輸送和利于地面O3增長的大氣垂直擴(kuò)散條件對O3污染的形成和爆發(fā)影響明顯.

    圖8 2018年8月9~11日天津市區(qū)氣溫、相對濕度、風(fēng)向風(fēng)速、紫外輻射強(qiáng)度,O3質(zhì)量濃度和TKE、垂直速度變化

    3 結(jié)論

    3.1 對2017~2019年冬夏季天津地區(qū)首要污染物PM2.5和O3質(zhì)量濃度和氣象要素的綜合分析發(fā)現(xiàn),PM2.5質(zhì)量濃度通常隨高度上升而下降,冬季,受大氣垂直擴(kuò)散條件影響,120m與地面PM2.5質(zhì)量濃度相關(guān)明顯,而200m則與地面PM2.5質(zhì)量濃度相關(guān)不明顯.夏季,在大氣湍流發(fā)展最為強(qiáng)盛的午后時段,通常出現(xiàn)120m和200m高于地面PM2.5質(zhì)量濃度的情況,120m和200m 的PM2.5質(zhì)量濃度變化較為一致,兩者相關(guān)系數(shù)達(dá)到0.72**.受地面排放源和大氣垂直擴(kuò)散條件影響,O3質(zhì)量濃度呈現(xiàn)隨高度上升的特征,地面O3質(zhì)量濃度日變化最為顯著.夏季,各高度O3質(zhì)量濃度差異低于冬季,地面與120m高度O3質(zhì)量濃度接近.

    3.2 對觀測期間氣象塔的TKE、垂直速度、大氣穩(wěn)定度和溫度層結(jié)的分析發(fā)現(xiàn),在水平風(fēng)場影響較小情況下,TKE與PM2.5和O3質(zhì)量濃度垂直分布有明顯影響.冬季,TKE與PM2.5質(zhì)量濃度的相關(guān)系數(shù)為-0.65**,夏季,TKE與ΔPM2.5的相關(guān)系數(shù)為-0.39**.夏季和冬季,TKE與地面O3濃度的相關(guān)系數(shù)分別為0.46**和0.53**,與ΔO3的相關(guān)系數(shù)分別為0.73**和0.70**.弱下沉運(yùn)動對地面O3濃度影響較強(qiáng),40m高度垂直運(yùn)動速度與地面O3濃度的相關(guān)系數(shù)在冬、夏季分別為-0.54**和-0.61**.以大氣穩(wěn)定度、TKE、逆溫強(qiáng)度及氣溫遞減率作為大氣垂直擴(kuò)散的指標(biāo),對PM2.5和O3濃度垂直分布具有一定指示作用.

    3.3 對冬季典型PM2.5重污染過程的分析發(fā)現(xiàn),霧霾的生消維持和PM2.5濃度的變化與大氣穩(wěn)定度、氣溫垂直遞減率和TKE的變化有直接關(guān)系.對夏季典型O3污染過程的分析發(fā)現(xiàn),O3污染形成與有利的光化學(xué)反應(yīng)氣象條件密切相關(guān),同時,垂直向下輸送和有利垂直擴(kuò)散條件對O3污染的形成和爆發(fā)影響明顯.

    [1] Wang T, Wei X, Ding A J, et al. Increasing surface ozone concentrations in the background atmosphere of Southern China, 1994~2007 [J]. Atmos. Chem. Phys., 2009,(9):6216-6226.

    [2] 程麟鈞,王 帥,宮正宇,等.中國臭氧濃度的時空變化特征及分區(qū) [J]. 中國環(huán)境科學(xué), 2017,37(11):5-14.

    Cheng L J, Wang S, Gong Z Y, et al. Spatial and seasonal variation and regionalization of ozone concentrations in China [J]. China Environ. Sci., 2017,37(11):5-14.

    [3] 蔣維楣.空氣污染氣象學(xué) [M]. 北京:氣象出版社, 2003:1-50.

    Jiang W M. Air pollution meteorology [M]. China Meteorological Press, 2003:1-50.

    [4] 李旭暉.邊界層氣象學(xué)基本原理 [M]. 北京:氣象出版社, 2018: 1-100.

    Li X H. Basic principles of boundary layer meteorology [M]. Beijing: China Meteorological Press, 2018:1-100.

    [5] 張人禾,李 強(qiáng),張若楠. 2013年1月中國東部持續(xù)性強(qiáng)霧霾天氣產(chǎn)生的氣象條件分析 [J]. 中國科學(xué):地球科學(xué), 2014,44(1):27-36.

    Zhang R H, Li Q, Zhang R N. Meteorological conditions for the persistent severe fog and haze event over Eastern China in January 2013 [J]. Sci. China Earth Sci., 2014,44(1):27-36.

    [6] 張小曳,孫俊英,王亞強(qiáng),等.我國霧-霾成因及其治理的思考 [J]. 科學(xué)通報, 2013,58(13):1178-1187.

    Zhang X Y, Sun J Y, Wang Y Q, et al. Factors contributing to haze and fog in China [J]. Chin. Sci. Bull. (Chin. Ver.), 2013,58(13):1178-1187.

    [7] 韓素芹.大氣邊界層對天津重污染天氣影響研究及應(yīng)用 [M]. 北京:氣象出版社, 2020:1-168.

    Han S Q. Influence of atmospheric boundary layer on heavy pollution weather in Tianjin and its application [M]. Beijing: China Meteorological Press, 2020:1-168.

    [8] Wang J, Zhou M, Liu B S, et al. Characterization and source apportionment of size-segregated atmospheric particulate matter collected at ground level and from the urban canopy in Tianjin [J]. Environ. Pollut., 2016,219:982-992.

    [9] Quan J, Zhang Q, He H, et al. Analysis of the formation of fog and haze in North China Plain (NCP) [J]. Atmos. Chem. Phys., 2011,(11): 8205-8214.

    [10] Han S Q, Liu J L, Hao T Y, et al. Boundary layer structure and scavenging effect during a typical winter haze-fog episode in a core city of BTH region, China [J]. Atmos. Environ., 2018,179:187-200.

    [11] Han S Q, Zhang Y F, Wu J H, et al. Evaluation of regional background particulate matter concentration based on vertical distribution characteristics [J]. Atmos. Chem. Phys., 2015,(15):11165-11177.

    [12] 姚 青,劉敬樂,蔡子穎,等.天津一次霧-霾天氣過程的近地層溫濕結(jié)構(gòu)和湍流特征分析[J]. 環(huán)境科學(xué)學(xué)報, 2018,38(10):3856-3867.

    Yao Q, Liu J L, Cai Z Y, et al. Analysis of temperature and moisture structure and turbulence characteristics of a smog and haze weather process in Tianjin [J]. Acta Sci. Circumst., 2018,38(10):3856-3867.

    [13] 姚 青,劉敬樂,蔡子穎,等.天津大氣穩(wěn)定度和逆溫特征對PM2.5污染的影響 [J]. 中國環(huán)境科學(xué), 2018,38(8):2865-2872.

    Yao Q, Liu J L, Cai Z Y, et al. Atmospheric stability and inverse temperature characteristics of Tianjin and its effect on PM2.5pollution [J]. China Environ. Sci., 2018,38(8):2865-2872.

    [14] 丁國安,陳尊裕,高志球,等.北京城區(qū)低層大氣PM10和PM2.5垂直結(jié)構(gòu)及其動力特征[J]. 中國科學(xué)(D輯:地球科學(xué)), 2005,35:31-44.

    Ding G A, Chen Z Y, Gao Z Q, et al. Vertical structure of PM10and PM2.5and their dynamical character in low atmosphere in Beijing urban areas [J]. Sci. China (Ser. D), 2005,35:31-44.

    [15] 李 昕,安俊琳,王躍思,等.北京氣象塔夏季大氣臭氧觀測研究[J]. 中國環(huán)境科學(xué), 2003,23(4):353-357.

    Li X, An J L, Wang Y S, et al. Studies on the measurement of atmospheric ozone in summer with Beijing meteorological tower [J]. China Environ. Sci., 2003,23(4):353-357.

    [16] Zhao H J, Wu B G, Liu J L, et al. Two air pollution events in the coastal city of Tianjin, north China plain [J]. Atmos. Pollut. Res., 2019,(10):1780-1794.

    [17] Li Q H, Wu B G, Liu J L, et al. Characteristics of the atmospheric boundary layer and its relation with PM2.5during haze episodes in winter in the North China Plain [J]. Atmos. Environ., 2020,223: 117265.

    [18] Liu Y S, Tang G Q, Zhou L B, et al. Mixing layer transport flux of particulate matter in Beijing, China [J]. Atmos. Chem. Phys., 2019, 19:9531-9540.

    [19] 孫兆彬,李梓銘,廖曉農(nóng),等.北京大氣熱力和動力結(jié)構(gòu)對污染物輸送和擴(kuò)散條件的影響 [J]. 中國環(huán)境科學(xué), 2017,37(5):1693-1705.

    Sun Z B, Li Z M, Liao X N, et al. Influence of Beijing atmospheric thermal and dynamics structure on pollutants transport and diffusion conditions [J]. China Environ. Sci., 2017,37(5):1693-1705.

    [20] Ren Y, Zhang H S, Wei W, et al. Comparison of the turbulence structure during light and heavy haze pollution episodes [J]. Atmos. Res., 2019,230:104645.

    [21] 李東宸,林慈哲,銀 燕.強(qiáng)對流天氣對O3和CO的垂直輸送作用 [J]. 應(yīng)用氣象學(xué)報, 2019,30(1):82-92.

    Li D C, Lin C Z, Yin Y. The vertical transport of the ozone and carbon monoxide by severe convective weather [J]. J. Appl. Meteor. Sci., 2019,30(1):82-92.

    [22] 蔡子穎,張 敏,韓素芹,等.天津重污染天氣混合層厚度閾值及應(yīng)用研究 [J]. 氣象, 2018,44(7):911-920.

    Cai Z Y, Zhang M, Han S Q, et al. Research on threshold and regularity of mixed layer thickness in heavy pollution weather in Tianjin [J]. Meteor. Mon., 2018,44(7):911-920.

    [23] Jiang Y Y, Xin J Y, Wang Y, et al. The thermodynamic structures of the planetary boundary layer dominated by synoptic circulations and the regular effect on air pollution in Beijing [J]. Atmos. Chem. Phys., 2021, (21):6111-6128.

    [24] 蔡子穎,韓素芹,張 敏,等.天津地區(qū)污染天氣分析中垂直擴(kuò)散指標(biāo)構(gòu)建及運(yùn)用 [J]. 環(huán)境科學(xué), 2018,39(6):2548-2556.

    Cai Z Y, Han S Q, Zhang M, et al. Construction and application of vertical diffusion index for analyzing weather during pollution events in Tianjin [J]. Chin. J. Envir. Sci., 2018,39(6):2548-2556.

    [25] 王 枚,鄭有飛,柳艷菊,等.京津冀臭氧變化特征及與氣象要素的關(guān)系 [J]. 中國環(huán)境科學(xué), 2019,37(7):2689-2698.

    Wang M, Zheng Y F, Liu Y J, et al. Characteristics of ozone and its relationship with meteorological factors in Beijing-Tianjin-Hebei Region [J]. China Environ. Sci., 2019,37(7):2689-2698.

    [26] 劉 健,吳 兌,范紹佳,等.前體物與氣象因子對珠江三角洲臭氧污染的影響 [J]. 中國環(huán)境科學(xué), 2017,37(3):813-820.

    Liu J, Wu D, Fan S J, et al. Impacts of precursors and meteorological factors on ozone pollution in Pearl River Delta [J]. China Environ. Sci., 2017,37(3):813-820.

    [27] 鄧雪嬌,李 菲,吳 兌,等.廣州地區(qū)典型清潔與污染過程的大氣湍流與物質(zhì)交換特征 [J]. 中國環(huán)境科學(xué), 2011,31(9):1424-1430.

    Deng X J, Li F, Wu D, et al. Turbulence and mass exchange characteristics of typical clean and pollution process over Guangzhou region [J]. China Environ. Sci., 2011,31(9):1424-1430.

    [28] 劉敬樂,姚 青,蔡子穎,等.基于天津255m氣象塔對近年天津近地面風(fēng)和氣溫變化特征的研究 [J]. 氣象, 2020,46(9):1235-1244.

    Liu J L, Yao Q, Cai Z Y, et al. Analysis of temperature and wind variation characteristics based on the 255m meteorological tower in Tianjin [J]. Meteor. Mon., 2020,46(9):1235-1244.

    Effects of atmospheric diffusion conditions on vertical distribution of pollutants in Tianjin.

    LIU Jing-le1,2,4, SHI Jing1, YAO Qing3,4, CAI Zi-ying3,4*, HAN Su-qin3, JIANG Ming1, CUI Ye1

    (1.Tianjin Meteorological Observation Center, Tianjin 300061;2.Tianjin Key Laboratory for Oceanic Meteorology, Tianjin 300074;3.Tianjin Environmental Meteorological Center, Tianjin 300074;4.Tianjin Institute of Meteorology, Tianjin 300074)., 2022,42(4):1575~1584

    Using the air pollutant monitoring data and meteorological observation data in summer and winter during 2017~2019 in Tianjin, the effects of atmospheric vertical diffusion conditions on PM2.5and O3were studied based on the vertical observation of Tianjin Meteorological Tower. The results showed that the mass concentration of near-surface PM2.5decreased with height, while the concentration of O3increased with height. Influenced by the seasonal and diurnal variation of atmospheric turbulence intensity, there was a significant correlation between the PM2.5mass concentration on the ground and at 120m height, but insignificant correlation between the PM2.5mass concentration on the ground and at 200m height in winter. In summer, the correlation coefficient of PM2.5mass concentration at 120m and 200m height was 0.72 and the PM2.5mass concentration at 120m and 200m height were even higher than the ground in the afternoon. In summer, the difference of O3mass concentration at different heights is lower than that in winter, and the O3mass concentration were close to the ground at 120m height. Atmospheric stability, temperature inversion intensity, and temperature decline rate were taken as atmospheric vertical diffusion indexes, which had certain indicator effects on the vertical gradient of PM2.5and O3. The correlation coefficient between TKE and PM2.5mass concentration was -0.65 in winter, and that between TKE and ΔPM2.5in summer was -0.39. The correlation coefficients between TKE and O3concentration were 0.46 and 0.53, in summer and winter respectively. The correlation coefficients between TKE and ΔO3concentration were 0.73 and 0.70, in summer and winter respectively. The weak downdraft has an obvious influence on the increase of O3concentration. The correlation coefficients between the vertical velocity at 40m height and O3concentration were -0.54 and -0.61 in winter and summer, respectively. Through the analysis of a typical PM2.5heavy pollution process and typical O3pollution processes, it is found that the changes of atmospheric stability, temperature decline rate and TKE were closely related to the generation, extinction, maintenance and variation of PM2.5pollutants. The formation of O3pollution near the surface was closely related to favorable photochemical reactions. In addition, the transport influence of downdraft should also be paid attention to during O3pollution processes.

    PM2.5;O3;meteorological tower;atmospheric vertical diffusion conditions;Tianjin

    X513

    A

    1000-6923(2022)04-1575-10

    劉敬樂(1986-),男,天津人,工程師,碩士,主要從事大氣邊界層和環(huán)境氣象研究.發(fā)表論文30篇.

    2021-09-02

    天津市自然科學(xué)基金(19JCQNJC08000);國家自然科學(xué)基金項(xiàng)目(41771242);天津市重大專項(xiàng)(18ZXAQSF00130,18ZXSZSF00160)

    *責(zé)任作者, 高級工程師, 120078030@163.com

    猜你喜歡
    穩(wěn)定度湍流大氣
    大氣的呵護(hù)
    軍事文摘(2023年10期)2023-06-09 09:15:06
    高穩(wěn)晶振短期頻率穩(wěn)定度的仿真分析
    重氣瞬時泄漏擴(kuò)散的湍流模型驗(yàn)證
    大氣古樸揮灑自如
    大氣、水之后,土十條來了
    多MOSFET并聯(lián)均流的高穩(wěn)定度恒流源研究
    工藝參數(shù)對橡膠球鉸徑向剛度穩(wěn)定度的影響
    旱澇不穩(wěn)定度的定量化研究及應(yīng)用
    “青春期”湍流中的智慧引渡(三)
    “青春期”湍流中的智慧引渡(二)
    国产午夜精品久久久久久一区二区三区| 亚洲精品,欧美精品| 中文字幕免费在线视频6| a级毛色黄片| 在现免费观看毛片| 2021少妇久久久久久久久久久| 色视频在线一区二区三区| 黑人高潮一二区| 观看av在线不卡| 日本av免费视频播放| 久久鲁丝午夜福利片| 久久青草综合色| 国产伦精品一区二区三区视频9| 天堂中文最新版在线下载| 久久狼人影院| 丝袜在线中文字幕| 日本爱情动作片www.在线观看| 日韩熟女老妇一区二区性免费视频| 国产极品粉嫩免费观看在线 | 久热久热在线精品观看| 一区二区三区四区激情视频| 亚洲av国产av综合av卡| 国产有黄有色有爽视频| 国产在线男女| 国产成人aa在线观看| 国产免费视频播放在线视频| 免费大片黄手机在线观看| 丝袜在线中文字幕| 久久青草综合色| 亚洲欧美精品自产自拍| 大又大粗又爽又黄少妇毛片口| 狂野欧美白嫩少妇大欣赏| 美女xxoo啪啪120秒动态图| 亚洲欧美成人综合另类久久久| 亚洲第一区二区三区不卡| 国产淫语在线视频| 麻豆成人午夜福利视频| 亚洲欧美日韩卡通动漫| 噜噜噜噜噜久久久久久91| 色婷婷av一区二区三区视频| 亚洲欧美一区二区三区黑人 | 秋霞在线观看毛片| www.色视频.com| 午夜福利网站1000一区二区三区| 久久婷婷青草| 国产深夜福利视频在线观看| 国产乱来视频区| 最后的刺客免费高清国语| 91久久精品国产一区二区三区| 亚洲天堂av无毛| 两个人的视频大全免费| 大香蕉久久网| 久久久久久久亚洲中文字幕| 久久久久久久国产电影| 亚洲精品456在线播放app| 国产精品久久久久成人av| 美女xxoo啪啪120秒动态图| 日本wwww免费看| 22中文网久久字幕| 精品熟女少妇av免费看| 汤姆久久久久久久影院中文字幕| 一区在线观看完整版| 成人国产av品久久久| 亚洲精品久久久久久婷婷小说| 国产精品不卡视频一区二区| 乱码一卡2卡4卡精品| 亚洲av不卡在线观看| 亚洲精品乱码久久久v下载方式| 国产成人免费无遮挡视频| 又粗又硬又长又爽又黄的视频| 国产精品秋霞免费鲁丝片| 一区二区三区精品91| 欧美精品一区二区大全| 五月开心婷婷网| 蜜桃在线观看..| 国产精品麻豆人妻色哟哟久久| 女的被弄到高潮叫床怎么办| 国精品久久久久久国模美| 精品一区在线观看国产| 欧美亚洲 丝袜 人妻 在线| 伦理电影大哥的女人| 国产亚洲av片在线观看秒播厂| av不卡在线播放| 国精品久久久久久国模美| 多毛熟女@视频| 菩萨蛮人人尽说江南好唐韦庄| 免费高清在线观看视频在线观看| 99视频精品全部免费 在线| 亚洲内射少妇av| tube8黄色片| 亚洲人成网站在线播| 免费av中文字幕在线| 久久久久国产网址| 亚洲精品中文字幕在线视频 | 国产精品国产三级国产av玫瑰| 日韩av不卡免费在线播放| 高清在线视频一区二区三区| 超碰97精品在线观看| 国产亚洲精品久久久com| av播播在线观看一区| 国产av精品麻豆| 日本黄色日本黄色录像| 老熟女久久久| 少妇人妻精品综合一区二区| 国模一区二区三区四区视频| 中文欧美无线码| 国产午夜精品久久久久久一区二区三区| 纵有疾风起免费观看全集完整版| 91精品一卡2卡3卡4卡| 在线观看免费视频网站a站| 欧美成人午夜免费资源| 欧美3d第一页| 人妻少妇偷人精品九色| 18禁在线无遮挡免费观看视频| 亚洲内射少妇av| 精品久久久久久电影网| 国产av精品麻豆| av线在线观看网站| av国产久精品久网站免费入址| 亚洲内射少妇av| 国产成人精品婷婷| 十分钟在线观看高清视频www | 91久久精品电影网| 边亲边吃奶的免费视频| 亚洲欧美日韩东京热| 久久鲁丝午夜福利片| 日本黄色日本黄色录像| 国产毛片在线视频| 高清视频免费观看一区二区| h视频一区二区三区| 亚洲av成人精品一区久久| 亚洲人成网站在线播| 国产无遮挡羞羞视频在线观看| 国产老妇伦熟女老妇高清| 国产免费一区二区三区四区乱码| 卡戴珊不雅视频在线播放| 国产高清国产精品国产三级| 草草在线视频免费看| 午夜激情久久久久久久| 丰满乱子伦码专区| 观看美女的网站| 菩萨蛮人人尽说江南好唐韦庄| 只有这里有精品99| 日韩一本色道免费dvd| 色视频www国产| 亚洲综合精品二区| 黄色毛片三级朝国网站 | 伦理电影免费视频| 国语对白做爰xxxⅹ性视频网站| 少妇的逼好多水| 岛国毛片在线播放| 精品一品国产午夜福利视频| 亚洲三级黄色毛片| 特大巨黑吊av在线直播| 国产乱来视频区| 女性被躁到高潮视频| 亚洲四区av| 一本大道久久a久久精品| av国产久精品久网站免费入址| 国产黄片视频在线免费观看| 免费av中文字幕在线| 成人影院久久| 狂野欧美激情性xxxx在线观看| 日本vs欧美在线观看视频 | 亚洲精品色激情综合| 少妇人妻 视频| 亚洲美女黄色视频免费看| 美女主播在线视频| 搡女人真爽免费视频火全软件| av免费观看日本| 国产 一区精品| 国产熟女欧美一区二区| 看免费成人av毛片| kizo精华| 久久久久久久久大av| 欧美老熟妇乱子伦牲交| 亚洲久久久国产精品| 欧美 日韩 精品 国产| 女的被弄到高潮叫床怎么办| 精品亚洲乱码少妇综合久久| 国产视频内射| 精品一区二区三区视频在线| 亚洲国产最新在线播放| 国国产精品蜜臀av免费| 国产亚洲5aaaaa淫片| 久久99蜜桃精品久久| 成人国产麻豆网| 中文精品一卡2卡3卡4更新| 极品人妻少妇av视频| 欧美日韩视频精品一区| 亚洲av二区三区四区| 免费观看av网站的网址| 免费看不卡的av| 又黄又爽又刺激的免费视频.| 亚洲国产毛片av蜜桃av| 婷婷色av中文字幕| 美女视频免费永久观看网站| 亚洲图色成人| 99热这里只有是精品在线观看| 欧美丝袜亚洲另类| 欧美日韩国产mv在线观看视频| 男男h啪啪无遮挡| 搡老乐熟女国产| 亚洲国产精品一区三区| 国产高清有码在线观看视频| 免费人妻精品一区二区三区视频| 国产黄频视频在线观看| av在线观看视频网站免费| 欧美日韩在线观看h| 人妻夜夜爽99麻豆av| 欧美精品人与动牲交sv欧美| 亚洲精品,欧美精品| 少妇被粗大猛烈的视频| av免费在线看不卡| 国产一区亚洲一区在线观看| 国产精品国产三级国产av玫瑰| 一边亲一边摸免费视频| 国产男女内射视频| 久久ye,这里只有精品| 亚洲情色 制服丝袜| 国产成人a∨麻豆精品| 免费观看的影片在线观看| 欧美丝袜亚洲另类| 国产日韩欧美在线精品| 日韩视频在线欧美| 亚洲精品一二三| 日韩,欧美,国产一区二区三区| 九草在线视频观看| 日本免费在线观看一区| 交换朋友夫妻互换小说| 国内揄拍国产精品人妻在线| 久久久久久伊人网av| 亚洲欧洲精品一区二区精品久久久 | 欧美一级a爱片免费观看看| 国产色婷婷99| kizo精华| .国产精品久久| 日韩欧美 国产精品| 日韩伦理黄色片| 日韩大片免费观看网站| 男女国产视频网站| 91精品国产国语对白视频| 99久久精品国产国产毛片| 另类精品久久| 高清在线视频一区二区三区| 在线播放无遮挡| 久久精品国产a三级三级三级| av黄色大香蕉| 日韩视频在线欧美| 天天操日日干夜夜撸| 国产亚洲91精品色在线| 乱系列少妇在线播放| 成人午夜精彩视频在线观看| 久久久久久久久久人人人人人人| 免费播放大片免费观看视频在线观看| 在线观看国产h片| 五月玫瑰六月丁香| a级毛片免费高清观看在线播放| 亚洲,一卡二卡三卡| 老司机影院成人| 超碰97精品在线观看| 九草在线视频观看| 免费播放大片免费观看视频在线观看| 日韩大片免费观看网站| 亚洲国产精品一区三区| 亚洲va在线va天堂va国产| 少妇人妻精品综合一区二区| 日本黄色片子视频| 亚洲成人手机| 下体分泌物呈黄色| 日本黄色片子视频| 久久久a久久爽久久v久久| 成年女人在线观看亚洲视频| 97超碰精品成人国产| 看免费成人av毛片| 国产av一区二区精品久久| 亚洲精品视频女| 久久99精品国语久久久| 日韩av在线免费看完整版不卡| 亚洲国产精品一区二区三区在线| 亚洲精品久久久久久婷婷小说| 亚洲综合色惰| 一本—道久久a久久精品蜜桃钙片| 丰满乱子伦码专区| videossex国产| 老司机影院成人| 久久久久久久大尺度免费视频| 久久人人爽人人片av| 中文字幕亚洲精品专区| 亚洲精品乱久久久久久| 日本色播在线视频| 观看av在线不卡| 国产精品人妻久久久影院| 777米奇影视久久| 久久精品久久精品一区二区三区| 在现免费观看毛片| 精品99又大又爽又粗少妇毛片| 久久久精品免费免费高清| 午夜视频国产福利| 久久精品国产亚洲av涩爱| 欧美日韩一区二区视频在线观看视频在线| 亚洲精品aⅴ在线观看| 汤姆久久久久久久影院中文字幕| 国产精品成人在线| 国产国拍精品亚洲av在线观看| 噜噜噜噜噜久久久久久91| 精品国产国语对白av| 人妻系列 视频| 久久久久久久久久久丰满| 夜夜爽夜夜爽视频| 少妇精品久久久久久久| 高清视频免费观看一区二区| 成人亚洲欧美一区二区av| 日韩人妻高清精品专区| 亚洲欧美精品专区久久| 国产精品国产三级国产av玫瑰| 女人精品久久久久毛片| xxx大片免费视频| 亚洲国产精品专区欧美| 狂野欧美激情性xxxx在线观看| 国产成人免费无遮挡视频| 成年人午夜在线观看视频| 韩国高清视频一区二区三区| 亚洲精品日本国产第一区| 久热久热在线精品观看| av视频免费观看在线观看| 亚洲欧美日韩东京热| 国产男女超爽视频在线观看| 国产亚洲5aaaaa淫片| 十分钟在线观看高清视频www | 亚洲av国产av综合av卡| 汤姆久久久久久久影院中文字幕| 爱豆传媒免费全集在线观看| 美女xxoo啪啪120秒动态图| 男人和女人高潮做爰伦理| 精品99又大又爽又粗少妇毛片| 18+在线观看网站| 国产精品久久久久久久电影| 午夜日本视频在线| 视频区图区小说| 亚洲av电影在线观看一区二区三区| 免费观看无遮挡的男女| 欧美日韩av久久| 国产91av在线免费观看| 久久毛片免费看一区二区三区| 久久av网站| 亚洲欧美一区二区三区黑人 | 亚洲av.av天堂| 中文字幕av电影在线播放| 亚洲欧美中文字幕日韩二区| 高清在线视频一区二区三区| 日韩欧美一区视频在线观看 | 欧美日韩亚洲高清精品| 日韩免费高清中文字幕av| 91精品伊人久久大香线蕉| 曰老女人黄片| 下体分泌物呈黄色| 丁香六月天网| 日韩成人伦理影院| 一个人看视频在线观看www免费| 男的添女的下面高潮视频| 久久久久久久久大av| 国产日韩一区二区三区精品不卡 | 91在线精品国自产拍蜜月| a级毛片在线看网站| 久久久久久伊人网av| a级毛片免费高清观看在线播放| 一本久久精品| 久久99热6这里只有精品| 国产高清三级在线| 亚洲欧美成人精品一区二区| 免费观看在线日韩| 久久久久久久久久人人人人人人| 国产真实伦视频高清在线观看| 青青草视频在线视频观看| 在线观看免费高清a一片| 卡戴珊不雅视频在线播放| 在线免费观看不下载黄p国产| 日本91视频免费播放| 精品99又大又爽又粗少妇毛片| 十八禁高潮呻吟视频 | 久久这里有精品视频免费| 看十八女毛片水多多多| 久久久久久久久久久久大奶| 天天操日日干夜夜撸| 五月伊人婷婷丁香| 狂野欧美激情性xxxx在线观看| 中文字幕精品免费在线观看视频 | 国产美女午夜福利| 国产一级毛片在线| 亚洲久久久国产精品| 高清av免费在线| 国产精品久久久久久久久免| 精品国产露脸久久av麻豆| 国产视频首页在线观看| 十八禁高潮呻吟视频 | 久久青草综合色| 亚洲三级黄色毛片| 五月伊人婷婷丁香| 在现免费观看毛片| 91精品一卡2卡3卡4卡| av网站免费在线观看视频| 一区二区av电影网| 免费观看性生交大片5| 日本爱情动作片www.在线观看| 午夜免费鲁丝| 色婷婷久久久亚洲欧美| 午夜激情福利司机影院| 精品亚洲成国产av| 大香蕉久久网| 高清黄色对白视频在线免费看 | 欧美精品国产亚洲| 又大又黄又爽视频免费| 噜噜噜噜噜久久久久久91| 日本欧美国产在线视频| 一个人免费看片子| 欧美xxⅹ黑人| 精华霜和精华液先用哪个| 欧美日韩视频高清一区二区三区二| 三级国产精品片| 夜夜爽夜夜爽视频| 国产精品一区二区在线观看99| 亚洲欧美成人综合另类久久久| 成人特级av手机在线观看| 国产亚洲精品久久久com| 国产日韩欧美在线精品| 成年美女黄网站色视频大全免费 | 91在线精品国自产拍蜜月| 亚洲欧美一区二区三区国产| a级一级毛片免费在线观看| 亚洲第一区二区三区不卡| 久久久久久伊人网av| 少妇猛男粗大的猛烈进出视频| 亚洲人与动物交配视频| 久久久久久久亚洲中文字幕| 国产熟女欧美一区二区| 视频区图区小说| 亚洲色图综合在线观看| 在线亚洲精品国产二区图片欧美 | 中文资源天堂在线| 夜夜看夜夜爽夜夜摸| 深夜a级毛片| 国产亚洲欧美精品永久| 亚洲国产精品成人久久小说| 97在线人人人人妻| 国产精品.久久久| 观看av在线不卡| 亚洲欧美日韩卡通动漫| 国产精品秋霞免费鲁丝片| 女人精品久久久久毛片| 国产精品久久久久久av不卡| 精品99又大又爽又粗少妇毛片| 亚洲精品一二三| 插逼视频在线观看| 99久久中文字幕三级久久日本| 精品少妇内射三级| 又爽又黄a免费视频| 熟女电影av网| 日日撸夜夜添| 美女福利国产在线| 国产亚洲av片在线观看秒播厂| 亚洲成人av在线免费| 亚洲自偷自拍三级| 亚洲欧美中文字幕日韩二区| av黄色大香蕉| 日韩欧美 国产精品| 噜噜噜噜噜久久久久久91| 简卡轻食公司| 亚洲精华国产精华液的使用体验| 中文天堂在线官网| 国产成人免费观看mmmm| 男人和女人高潮做爰伦理| 大又大粗又爽又黄少妇毛片口| 九九在线视频观看精品| 99九九线精品视频在线观看视频| 成人综合一区亚洲| 三级国产精品欧美在线观看| 久久国产亚洲av麻豆专区| 国产免费一区二区三区四区乱码| 国产毛片在线视频| 啦啦啦中文免费视频观看日本| 波野结衣二区三区在线| 日韩欧美精品免费久久| 永久网站在线| 国产黄色视频一区二区在线观看| 少妇人妻一区二区三区视频| 精品99又大又爽又粗少妇毛片| 精华霜和精华液先用哪个| 高清在线视频一区二区三区| 国产精品免费大片| a级片在线免费高清观看视频| 久久综合国产亚洲精品| 免费av不卡在线播放| 男人和女人高潮做爰伦理| 日韩av在线免费看完整版不卡| 国产无遮挡羞羞视频在线观看| 国产一区二区在线观看日韩| 国产精品麻豆人妻色哟哟久久| 免费av中文字幕在线| 精品人妻熟女毛片av久久网站| 精品国产国语对白av| 欧美日本中文国产一区发布| 我的女老师完整版在线观看| 国产亚洲一区二区精品| 男女免费视频国产| 亚洲在久久综合| 国产黄频视频在线观看| 久久久精品94久久精品| 三级国产精品片| 少妇被粗大的猛进出69影院 | 亚洲av在线观看美女高潮| 18禁裸乳无遮挡动漫免费视频| 80岁老熟妇乱子伦牲交| 亚洲av电影在线观看一区二区三区| 黄色一级大片看看| 色婷婷久久久亚洲欧美| 黄片无遮挡物在线观看| 精品亚洲乱码少妇综合久久| 一级片'在线观看视频| 亚洲精华国产精华液的使用体验| av卡一久久| 国产免费一级a男人的天堂| 丰满饥渴人妻一区二区三| 国产精品秋霞免费鲁丝片| 美女中出高潮动态图| 综合色丁香网| 一级毛片黄色毛片免费观看视频| a级毛色黄片| 精品国产一区二区久久| 99久久精品一区二区三区| 久久青草综合色| 精品国产露脸久久av麻豆| 最近的中文字幕免费完整| 观看av在线不卡| 久久久久久久亚洲中文字幕| 一级二级三级毛片免费看| 国产男女超爽视频在线观看| 最近的中文字幕免费完整| 人人妻人人澡人人看| 国产淫语在线视频| av在线播放精品| 亚洲欧美精品自产自拍| 青春草国产在线视频| 日本vs欧美在线观看视频 | 国产av国产精品国产| 欧美xxⅹ黑人| 特大巨黑吊av在线直播| 亚洲av福利一区| 内射极品少妇av片p| 国精品久久久久久国模美| 人妻一区二区av| 精品久久国产蜜桃| 亚洲国产精品成人久久小说| 成人综合一区亚洲| 亚洲国产av新网站| 三级国产精品片| 亚洲av日韩在线播放| 一级毛片我不卡| 在线精品无人区一区二区三| xxx大片免费视频| 国产爽快片一区二区三区| 久久久午夜欧美精品| 久久久久久久精品精品| av视频免费观看在线观看| 丰满人妻一区二区三区视频av| 久久久久国产精品人妻一区二区| 久久久国产精品麻豆| 亚洲内射少妇av| 精品亚洲乱码少妇综合久久| 国产黄频视频在线观看| 午夜福利视频精品| 国产成人一区二区在线| 激情五月婷婷亚洲| h视频一区二区三区| 色吧在线观看| 乱人伦中国视频| 少妇被粗大猛烈的视频| 亚洲精品乱码久久久v下载方式| √禁漫天堂资源中文www| 三级经典国产精品| 最近的中文字幕免费完整| 99热网站在线观看| 熟女av电影| 男人和女人高潮做爰伦理| 国产欧美日韩综合在线一区二区 | 在线看a的网站| 三上悠亚av全集在线观看 | 国产成人一区二区在线| 老司机亚洲免费影院| 最黄视频免费看| av专区在线播放| 亚洲成人手机| 蜜臀久久99精品久久宅男| 国产一区有黄有色的免费视频| 久久97久久精品| 天天躁夜夜躁狠狠久久av| 国产在线一区二区三区精| 久久久久久久久大av| 十八禁高潮呻吟视频 | 如何舔出高潮| 日本与韩国留学比较| 最近2019中文字幕mv第一页| 亚洲av中文av极速乱| 欧美xxxx性猛交bbbb| 国内少妇人妻偷人精品xxx网站| 曰老女人黄片| 久久久久久久久久久久大奶| 狂野欧美白嫩少妇大欣赏| 国产黄色免费在线视频| 国产 精品1| 汤姆久久久久久久影院中文字幕| 性色av一级| 国产精品无大码| 99热这里只有是精品在线观看| 国产精品国产三级专区第一集|