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

    北京冬季城市人為水汽排放的定量估算研究

    2021-08-11 10:07:28徐一平苗世光趙德龍丁德平
    氣候與環(huán)境研究 2021年4期
    關(guān)鍵詞:邊界層人為過量

    徐一平 苗世光 趙德龍 丁德平

    1 中國氣象科學(xué)研究院,北京 1000812 北京城市氣象研究院,北京 1000893 北京市人工影響天氣辦公室,北京 1000894 云降水物理研究和云水資源開發(fā)北京市重點(diǎn)實(shí)驗(yàn)室,北京 1000895 中國氣象局華北云降水野外科學(xué)試驗(yàn)基地,北京 101200

    1 引言

    城市化是人類社會現(xiàn)代化發(fā)展的必然趨勢,它使得土地利用和土地覆蓋都發(fā)生了很大的變化,人為熱排放、城市氣溶膠和人為溫室氣體排放也大大增加,對天氣、氣候和環(huán)境均產(chǎn)生了不同程度的影響,天氣、氣候和環(huán)境反過來也會對城市的運(yùn)行和居民的生產(chǎn)生活產(chǎn)生影響,城市化程度體現(xiàn)了一個國家的經(jīng)濟(jì)發(fā)達(dá)程度,是工業(yè)化水平高低的重要指標(biāo)(苗世光等, 2014b)。截至2019 年底,我國城市化率已達(dá)到60.6%,城鎮(zhèn)人口突破8.4×109(中華人民共和國國家統(tǒng)計(jì)局, 2020)。

    目前國內(nèi)大部分的研究主要集中在城市熱島、城市大氣環(huán)境上(汪光燾等, 2005; 壽亦萱和張大林, 2012; Deosthali, 2000),而對于城市人為水汽研究的相對較少的。國外在這方面已經(jīng)開展了一些相關(guān)研究。Moriwaki et al.(2008)使用地理信息系統(tǒng)和能源消耗數(shù)據(jù)庫的數(shù)據(jù),估算了東京人為水汽排放量和人為熱的時空變化,發(fā)現(xiàn)夏季東京市中心的人為水汽最大值超過500 W/m-2。冬季東京居民區(qū)的人為熱排放量大于夏季商業(yè)區(qū)的人為熱排放量,并且人為熱和水汽,使氣溫升高了約0.5°C,水汽混合比增加0.1 g/kg(Senoo et al., 2004)。在美國進(jìn)行的飛機(jī)觀測實(shí)驗(yàn)中,在城市下風(fēng)方向觀察到較高的水汽摩爾分?jǐn)?shù),水汽的排放量的估計(jì)值在1.6×104和1.7×105kg/s之間,占城市研究點(diǎn)大氣邊界層總的水汽平流的8.4%(Salmon et al., 2017)。

    人為水汽的排放也會對城市的天氣、氣候和環(huán)境等帶來多方面的影響。天然氣燃燒產(chǎn)生的水汽會使空氣飽和,在低溫的共同作用下,產(chǎn)生了埃德蒙頓的低溫霧(Rage, 1972)。城市的過量水汽與城市熱島之間也存在著正相關(guān),但是城市過量水汽對城市熱島的增強(qiáng)作用較弱,并且存在相當(dāng)大的負(fù)反饋(Holmer and Eliasson, 1999)。通常城市熱島最強(qiáng)時,經(jīng)常觀察到城市過量水汽(Kuttler et al.,2007)。另一方面,來自城市的過量水汽也會導(dǎo)致城區(qū)下風(fēng)向的降水增加(M?lders and Olson, 2004)。由此可以看出,城市水汽的排放對于城市的天氣、環(huán)境等也有不可低估的影響。

    作為我國的首都和政治文化中心,北京市人口規(guī)模大,經(jīng)濟(jì)活動發(fā)達(dá),屬于特大城市,近幾十年來,北京的城市建設(shè)十分迅速,經(jīng)歷了高速發(fā)展期,城市結(jié)構(gòu)和建筑高度等都發(fā)生巨變(徐陽陽等,2009)。隨著城市規(guī)模迅速擴(kuò)大和城市化步伐的加快,北京的氣象條件發(fā)生了一系列的變化。因城市效應(yīng),北京冬季表現(xiàn)為弱濕島,城區(qū)與郊區(qū)比濕日變化有明顯差異(竇晶晶等, 2014)。北京城區(qū)冬季寒冷干燥,波恩比較大,白天儲熱在凈輻射中所占比例最大,人為熱中潛熱和感熱表現(xiàn)為同等重要,多云時,冬季夜間的感熱略大于潛熱(苗世光等,2012)。在城市區(qū)域,人為潛熱排放與人類活動的規(guī)律基本吻合,比如空調(diào)系統(tǒng)、居民生活排放和汽車尾氣排放等均會導(dǎo)致人為潛熱的大量排放(苗世光等, 2014a)。進(jìn)一步研究發(fā)現(xiàn),城市對于北京地區(qū)的濕度影響顯著(Liu et al, 2009)??傮w上來說,城市在白天會比周圍地區(qū)更加干燥(Arnfield, 2003)。

    根據(jù)北京市統(tǒng)計(jì)局資料顯示,2019年北京市常住人口2153.6×104,常住人口密度達(dá)到1312 km-2(北京市統(tǒng)計(jì)局, 2020)。人口密度大,人類活動強(qiáng)烈,人為水汽排放量巨大,十分具有研究意義。本文利用飛機(jī)觀測數(shù)據(jù)量化了北京整個城市內(nèi)水汽摩爾分?jǐn)?shù)的增強(qiáng)情況,并討論了可能的來源,將為數(shù)值模式中人為潛熱的設(shè)置提供參考。

    2 研究方法

    本研究采用篩選出的2016年12月16日、17日、18日、19日符合質(zhì)量平衡飛行的飛機(jī)觀測數(shù)據(jù),計(jì)算了城市過量水汽通量、城市過量水汽排放率和城市水汽總排放率。

    2.1 航線位置

    北京市氣象局人工影響天氣辦公室的飛行觀測實(shí)驗(yàn),每次均從北京北部郊區(qū)出發(fā),環(huán)繞北京中心城區(qū)飛行觀測,并飛行多圈,最后返回北部郊區(qū),完成一次飛行觀測實(shí)驗(yàn),每次飛行的航線基本類似。

    2.2 機(jī)載平臺和機(jī)載儀器

    2.2.1 機(jī)載平臺和機(jī)載儀器的具體型號

    本文觀測儀器的機(jī)載平臺為北京市人工影響天氣辦公室空中國王3587飛機(jī)。本文使用的數(shù)據(jù)來自機(jī)載氣象探測系統(tǒng)AIMMS-20(Aircraft-Integrated Meteorological Measurement System-20)。該系統(tǒng)可用于飛行時的溫度、濕度、風(fēng)的三維空間狀況。該儀器包括空氣數(shù)據(jù)探頭,通過壓差和逆流溫度/濕度測量室測量氣流和空氣溫度濕度;慣性測量單元測量整個航跡中飛機(jī)高度數(shù)據(jù);GPS測量飛機(jī)高度和機(jī)身傾斜角度。

    2.2.2 機(jī)載儀器的指標(biāo)

    風(fēng)速精度水平分量為0.50 m/s,垂直分量為0.75 m/s;溫度精度為0.3°C,分辨率為0.01°C;相對濕度精度為2%,分辨率為0.1%;工作溫度為-20~50°C。

    2.3 數(shù)據(jù)

    本文使用數(shù)據(jù)是北京市氣象局人工影響天氣辦公室3830飛機(jī)觀測資料,篩選出溫度、相對濕度、氣壓、位置(經(jīng)緯度)、飛行高度、東西向和南北向飛行速度等數(shù)據(jù)。因組織一次專門飛機(jī)觀測實(shí)驗(yàn)的成本較高,本文研究所選用的數(shù)據(jù)是從北京市氣象局人工影響天氣辦公室的日常作業(yè)觀測的3830飛機(jī)觀測數(shù)據(jù)中篩選出來的,因此存在數(shù)據(jù)針對性不強(qiáng),很多數(shù)據(jù)不滿足計(jì)算條件而無法使用的情況,但是通過一定的篩選條件,依然能夠從廣泛的數(shù)據(jù)中可以找到可使用的部分,為本次研究提供數(shù)據(jù)基礎(chǔ)。

    2.4 質(zhì)量平衡飛行

    篩選出的飛機(jī)航線所用時間均為數(shù)分鐘,以盡量減少飛行期間大氣邊界層的增長。飛機(jī)在北京市的上風(fēng)和下風(fēng)向飛行,城市過量水汽排放率根據(jù)中心城區(qū)相對于背景的摩爾分?jǐn)?shù)的增加、相對于飛行軌跡的風(fēng)速垂直分量和當(dāng)時的邊界層高度計(jì)算。圖1、圖2 以2016年12月16日的一次飛行為例,顯示了飛機(jī)觀測航線和飛行高度變化的大致情況,所有日期飛行觀測的航線大致類似,并且在該次飛行中觀測到了過量水汽排放信號,即人為水汽排放。

    圖1 2016年12月16日的飛機(jī)觀測航線其中一圈?;疑珔^(qū)域表示城市區(qū)域,紅色曲線表示航線,黑色環(huán)線為北京二至五環(huán)路環(huán)線Fig. 1 One circle of the aircraft observation route on 16 December 2016; the gray area represents the urban area, the red curve represents the route,and the Beijing Second to Fifth Ring Road are drawn

    圖2 2016年12月16日的一次觀測飛行中飛行高度隨飛行時間變化Fig. 2 Change of flight altitude with flight time in a flight on 16 December 2016

    2.5 背景濃度的確定

    該研究的目的是確定城市過量水汽排放率。背景濃度可作為確定中心城區(qū)水汽摩爾分?jǐn)?shù)增量的參考。本次研究中采用中心城區(qū)周邊區(qū)域作為背景,周邊區(qū)域人口聚集度較小,相對中心城區(qū),人為潛熱排放量也較小。具體計(jì)算的過程中,即根據(jù)位于城區(qū)下風(fēng)向航線兩側(cè)邊緣的數(shù)據(jù)確定背景水汽摩爾分?jǐn)?shù),背景濃度與水平方向飛行距離線性相關(guān)。

    2.6 天氣形勢

    根據(jù)篩選出的航線的飛行時間,從北京市氣象臺獲取了2016年12月16~19日每天08:00(北京時間,下同)海平面、850 hPa、700 hPa、500 hPa的天氣圖。從海平面氣壓場等壓線分布圖可見,12月16日北京地區(qū)處于弱高壓前部的均壓場中,弱高壓的中心位于北京的西部;12月17日高壓系統(tǒng)東移入海,北京處于兩個高壓之間的低壓輻合區(qū)控制;12月18日北京持續(xù)處于兩個高壓之間;12月19日北京上游高壓逐漸東移南下,北京處于高壓前部的弱氣壓場中。大氣層結(jié)垂直方向上,16~19日,500、700、850 hPa高度上,北京均位于槽后,為下沉氣流??傮w看來,2016年12月16~19日,北京受高壓控制,無降水,水平方向風(fēng)速較小,無明顯暖濕氣流,天氣穩(wěn)定為主。

    2.7 城市過量水汽通量、城市過量水汽排放率、城市水汽總排放率的定義和計(jì)算方法

    城市過量水汽通量是指單位時間內(nèi),水平方向上流經(jīng)單位面積城市比周圍鄉(xiāng)村多的水汽量,單位用mol m-2s-1表示;城市過量水汽排放率指的是單位時間內(nèi),城市比周圍鄉(xiāng)村多排放的水汽質(zhì)量,單位用kg/s表示;城市水汽總排放率指的是單位時間內(nèi),整個城市排放的水汽質(zhì)量,單位用kg/s表示。對于城市過量水汽通量FUrban,i,計(jì)算公式如下:

    其中,i表示每一個下風(fēng)向數(shù)據(jù)點(diǎn)的標(biāo)號,U⊥,i表示垂直于航線的水平風(fēng)分量,cdw,i為下風(fēng)向水汽摩爾分?jǐn)?shù),cbg,i為背景水汽摩爾分?jǐn)?shù),使用理想氣體定律,可將水汽摩爾分?jǐn)?shù)轉(zhuǎn)化成水汽摩爾密度。相同氣壓和溫度下,下風(fēng)向水汽摩爾密度減去相應(yīng)的背景水汽摩爾密度,是城市過量水汽摩爾密度,城市過量水汽摩爾密度乘以風(fēng)速相對于航線的垂直分量就可以得到城市過量水汽通量FUrban,i,同理,城市水汽通量FTotal,i計(jì)算公式如下:

    為了計(jì)算城市范圍內(nèi)的水汽排放率,需要知道城市過量水汽通量矩陣,并做積分

    其中,H為大氣邊界層高度,x為水平飛行距離的一半,M是指城市過量水汽通量矩陣或者城市水汽總通量矩陣,進(jìn)一步比較可以確定城市過量水汽排放率(EUrban) 與城市水汽總排放率(ETotal)的比值。

    3 計(jì)算過程

    3.1 飽和水汽壓計(jì)算

    飛機(jī)觀測的氣壓范圍大致在640~1020 hPa,仍以0°C作為水的凝固點(diǎn),飽和水汽壓es的計(jì)算,溫度在0°C以上使用水面飽和水汽壓公式:

    溫度在0°C以下使用冰面飽和水汽壓公式:其中,t為攝氏溫度,T為熱力學(xué)溫度。

    3.2 水汽摩爾分?jǐn)?shù)的計(jì)算

    其中,χV為水汽摩爾分?jǐn)?shù),e為水汽壓,p為氣壓,UW為 相對濕度;es(T)為飽和水汽壓,由公式(4)、公式(5)可得;χVS為飽和水汽摩爾分?jǐn)?shù)。根據(jù)公式(6)和公式(7)則

    3.3 將水汽摩爾分?jǐn)?shù)轉(zhuǎn)化為摩爾密度

    利用水汽的理想氣體狀態(tài)方程:

    其中,n表示水汽的物質(zhì)的量,V表示氣體體積,可得水汽摩爾密度,R=8.314 J mol-1K-1。

    3.4 數(shù)據(jù)的篩選

    北京供暖時間為11月15日左右,飛機(jī)飛行日期為2016年12月16日、17日、18日、19日,北京地區(qū)已完全進(jìn)入冬季,之所以選擇冬季的數(shù)據(jù),是因?yàn)槎镜恼舭l(fā)速率和飽和水汽壓低于夏季,并且冬季人為水汽排放量大、背景水汽排放量小,易于估算城市人為水汽排放通量和排放率。

    根據(jù)公式(1)~(3)判斷,飛機(jī)觀測數(shù)據(jù)的篩選條件主要包括以下6點(diǎn):(1)飛機(jī)觀測過程中,風(fēng)向需指向城市外側(cè),即使得城市位于下風(fēng)向; (2)飛機(jī)觀測飛行高度需要在探測時大氣邊界層高度以下;(3)飛行路徑需要圍繞北京市區(qū),郊區(qū)的飛行路徑屬于無效路徑;(4)飛行路徑的側(cè)邊界盡量遠(yuǎn)離城區(qū);(5)飛機(jī)飛行路徑可根據(jù)經(jīng)緯度信息進(jìn)行判斷,當(dāng)飛行軌跡不為大致南北、東西向時,需要考慮飛行軌跡與正北向的夾角;(6)數(shù)據(jù)需無缺測。利用以上的這些篩選條件,本研究對2016年12月16日、12月17日、12月18日、12月19日的飛機(jī)探測數(shù)據(jù)進(jìn)行了篩選,最后得到了可用數(shù)據(jù)。

    3.5 大氣邊界層高度的計(jì)算

    在本研究中邊界層高度的確定采用位溫梯度法,位溫梯度法計(jì)算邊界層高度的依據(jù)是:近地層以上 (不穩(wěn)定邊界層時為混合層以上),位溫梯度最大的高度定義為邊界層高度。由于沒有連續(xù)的位溫資料,此處采用中央差分:

    其中,n=2,3,4,···,θn+1表示n+1點(diǎn)處的位溫,θn-1表示n-1處的位溫, Δz表示n+1點(diǎn)與n-1點(diǎn)之間的距離,根據(jù)飛機(jī)探空曲線可算出邊界層高度。之所以采用1次觀測代表大氣邊界層高度,是因?yàn)樵跊]有大尺度天氣系統(tǒng)影響下,混合層蓋頂逆溫位置變化相對比較緩慢,1天中1次觀測的蓋頂逆溫底高度可以代表這1天的情況,混合層頂覆蓋逆溫的位置就是大氣邊界層高度(袁仁民等, 2003)。

    4 結(jié)果分析

    2016年12月16日位溫廓線和可用航線選擇的結(jié)果如圖3所示,利用上升階段飛機(jī)觀測到的氣壓、溫度、高度數(shù)據(jù),可繪制位溫廓線,通過位溫梯度法,920 m處位溫梯度最大,為0.022 K/m,確定的大氣邊界層高度為920 m,并通過6點(diǎn)篩選條件確定了07:30:57至07:33:30自東向西、07:40:15至07:42:45自南向北的可用航線。

    圖3 (a)2016年12月16日北京的位溫廓線;(b)本文選取的2016年12月16日07:30:57至07:33:30(紅色)、07:40:15至07:42:45 (綠色)可用飛行路徑Fig. 3 (a) Potential temperature profile on 16 December 2016 in Beijing; (b) the red and green are the selected available flight path maps during 0730:57–0733:30 BT (Beijing time) (red) and 0740:15–0742:45 BT (green) on 16 December 2016, respectively

    12月17日位溫廓線和可用航線選擇的結(jié)果如圖4所示,同樣利用上升階段飛機(jī)觀測到的氣壓、溫度、高度數(shù)據(jù),通過計(jì)算可繪制位溫廓線,通過位溫梯度法,888 m處位溫梯度最大,為0.053 K/m,確定的大氣邊界層高度為888 m,并確定了06:39:33至06:43:00有一定傾角的自北向南、07:03:38至07:06:16和07:06:58至07:10:26兩條自西向東的可用航線。

    圖4 (a)2016年12月17日北京的位溫廓線;(b)本文選取的2016年12月17日06:39:33至06:43:00(紅色)、07:03:38至07:06:16( 綠色)、07:06:58至07:10:26(黃色)可用飛行路徑Fig. 4 (a) Potential temperature profile on 17 December 2016 in Beijing; (b) the red, green, and yellow are selected during 0639:33–0643:00 BT,0703:38–0706:16 BT, and 0706:58–0710:26 BT on 17 December 2016, respectively

    12月18日位溫廓線和可用航線選擇的結(jié)果如圖5所示,類似的,利用上升階段飛機(jī)觀測到的氣壓、溫度、高度數(shù)據(jù),繪制位溫廓線,570 m處位溫梯度最大,為0.031 K/m,確定的大氣邊界層高度為570 m,并確定了05:33:56至05:36:05較小傾角自北向南、05:37:04至05:39:57自西向東、05:40:32至05:43:24自南向北的可用航線。

    圖5 (a)2016年12月18日北京的位溫廓線;(b)本文選取的2016年12月18日05:33:56至05:36:05(紅色)、05:37:04至05:39:57 (綠色)、05:40:32至05:43:24(橙色)可用飛行路徑Fig. 5 (a) Potential temperature profile on 18 December 2016 in Beijing; (b) the red, green, and orange are selected available flight path diagrams during 0533:56–0536:05 BT, 0537:04–0539:57 BT, and 0540:32–0543:24 BT on 18 December 2016, respectively

    12月19日位溫廓線和可用航線選擇的結(jié)果如圖6所示,通過位溫廓線,617 m處位溫梯度最大,為0.023 K/m,確定的大氣邊界層高度為617 m,并確定了06:02:49至06:05:38和06:16:42至06:19:28自西向東,06:03:31至06:09:17自南向北的可用航線。

    圖6 (a)2016年12月19日北京的位溫廓線;(b)本文選取的2016年12月19日06:02:49至06:05:38(黃色)、06:03:31至06:09:17 (綠色)、06:16:42至06:19:28(紅色)可用飛行路徑Fig. 6 (a) Potential temperature profile on 19 December 2016 in Beijing; (b) the yellow, green, and red are the selected available flight path diagrams during 0602:49–0605:38 BT, 0603:31–0609:17 BT, and 0616:42–0619:28 BT on 19 December 2016, respectively

    在確定了可以使用的觀測數(shù)據(jù)之后,利用公式 (1)~(3),計(jì)算可得城市過量水汽通量(FUrban)、城市過量水汽排放率(EUrban)和城市水汽總排放率(ETotal)、城市過量水汽排放率(EUrban)占城市水汽總排放率(ETotal)的比例,如表1所示。

    表1 飛行時間、城市過量水汽通量(FUrban)、城市過量水汽排放率(EUrban)和城市水汽總排放率(ETotal)及其比例Table 1 Flight time, urban anthropogenic water vapor flux (FUrban), urban anthropogenic water vapor emission rate (EUrban),and total urban water vapor emission rate (ETotal) and their proportions

    在北京附近連續(xù)4天的飛機(jī)探測中選取的11次可用的飛行路徑中,有5次(2016年12月16日07:30:57至07:33:30、2016年12月18日05:33:56至05:36:05、2016年12月18日05:40:32至05:43:24、2016年12月19日06:06:31至06:09:17、2016年12月19日06:16:42至06:19:28)觀察到中心城區(qū)的城市水汽信號升高,可認(rèn)為這些時間內(nèi)中心城區(qū)相對于周邊區(qū)域存在水汽的過量排放。在篩選出的中心城區(qū)水汽過量排放的實(shí)例中,表現(xiàn)為中心城區(qū)的水汽摩爾分?jǐn)?shù)超過周邊地區(qū),過量水汽排放的占比在0.64%~13.12%。過量水汽排放的信號大小每時每刻都是在變化的,在有的飛行探測試驗(yàn)中,邊界層的垂直混合程度也會產(chǎn)生影響。

    從表格中可以看到在過量水汽排放的計(jì)算中,也會出現(xiàn)出一些負(fù)值,說明隨著北京周邊城市的快速發(fā)展,被作為背景濃度考慮的四環(huán)到五環(huán)之間的水汽排放也逐漸增多,甚至在某些時候已經(jīng)超過中心城區(qū)的水汽排放量??赡茉蚴请S著經(jīng)濟(jì)的發(fā)展,北京四環(huán)與五環(huán)之間成為新的經(jīng)濟(jì)中心,城市建設(shè)的快速發(fā)展與人口的大量集聚,導(dǎo)致人為潛熱的排放大量增加。對于2016年12月16日07:40:15至07:42:45飛行中計(jì)算出的城市過量水汽通量和城市過量水汽排放率的負(fù)值,主要原因是其作為背景濃度的一個取值點(diǎn)位于北京的望京區(qū)域,該區(qū)域商業(yè)和寫作樓密集,存在較大的人為水汽排放;另一個點(diǎn)位于北京經(jīng)濟(jì)技術(shù)開發(fā)區(qū)附近,產(chǎn)業(yè)密集,人口眾多,同樣存在較大的人為水汽排放,因此使得計(jì)算過程中確定的背景濃度相對較高,導(dǎo)致了城市過量水汽通量和排放率負(fù)值的出現(xiàn)。對于2016年12月17日06:39:33至06:43:00飛行中計(jì)算出的城市過量水汽通量和城市過量水汽排放率的負(fù)值,主要原因是作為背景濃度的一個取值點(diǎn)位于中關(guān)村科技園附近,經(jīng)濟(jì)和科技發(fā)達(dá),人員和商業(yè)密集,人為水汽排放量巨大;另一個背景濃度取值點(diǎn)位于永定河附近,有眾多旅游景點(diǎn)、企業(yè)和居民區(qū),同樣有較大的人為水汽排放量,同16日的飛行類似,計(jì)算中確定的背景濃度相對較大,導(dǎo)致了負(fù)值的出現(xiàn)。其他負(fù)值的分析類似。

    本文根據(jù)城區(qū)下風(fēng)向航線兩個側(cè)邊緣的數(shù)據(jù)確定背景水汽摩爾分?jǐn)?shù),若航線側(cè)邊緣處城市化水平較高,人為水汽排放量大于中心城區(qū),就會導(dǎo)致負(fù)值的出現(xiàn)。

    如果有專門的飛行試驗(yàn),使得飛機(jī)航線的兩個側(cè)邊緣遠(yuǎn)離城區(qū)和大型水體,更好地代表鄉(xiāng)村,則不會出現(xiàn)負(fù)值。

    將所有觀測到水汽過量信號的航線及所對應(yīng)的城市過量水汽排放率占城市水汽總排放率的比例繪制在一張圖上,如圖7所示,在觀測到水汽過量信號的航線中,北京城市北部地區(qū)的城市過量水汽排放率占城市水汽總排放率的比例顯著大于北京城市南部地區(qū),而東西兩側(cè)差異較小,在城市水汽總排放率相差不大的前提下,表明北京城市北部地區(qū)的人為潛熱排放顯著大于北京城市南部地區(qū)。

    圖7 2016年12月16日07:30:57至07:33:30(紅色)、18日05:33:56至05:36:05(黃色)、18日05:40:32至05:43:24(綠色)、19日06:06:31至06:09:17(藍(lán)色)、19日06:16:42至06:19:28(橙色)測得城市過量水汽信號的航線及由此計(jì)算出的城市過量水汽排放率占城市水汽總排放率的比例(EUrban/ETotal)Fig. 7 Route of the measured urban water vapor excess signal during 0730:57-0733:30 BT (red) on 16 December 2016, 0533:56-0536:05 BT(yellow), 0540:32-0543:24 BT (green) on 18 December 2016, 0606:31-0609:17 BT (blue), 0616:42-0619:28 BT (orange) on 19 December 2016 and the ratio of the calculated urban water vapor excess emission rate to the total urban water vapor emission rate (EUrban/ETotal)

    對于北京城市北部地區(qū)人為潛熱大于南部的現(xiàn)象,主要從天氣背景和城市內(nèi)部結(jié)構(gòu)兩個方面分析。天氣背景方面,如果飛行當(dāng)日北京相對濕度較大,人為水汽占比小,計(jì)算會出現(xiàn)一定的誤差。城市內(nèi)部結(jié)構(gòu)方面,這可能跟北京城市北部城市化程度高、人口密度較大有關(guān)。

    5 結(jié)論與討論

    本文利用飛機(jī)觀測資料對北京冬季人為水汽排放進(jìn)行了定量估算,主要計(jì)算了城市過量水汽通量、城市過量水汽排放率、城市水汽總排放率以及城市過量水汽排放率占城市水汽總排放率的比例,得到如下的主要結(jié)論:

    (1)在北京城區(qū)下風(fēng)方向進(jìn)行的飛機(jī)觀測中多次觀測到較高的水汽摩爾分?jǐn)?shù),當(dāng)觀測到城市水汽過量信號時,水汽排放率的估計(jì)值在57.64到9538.87 kg/s之間,占城市大氣邊界層水汽總平流 (背景+城市過量)的0.64%到13.12%。

    (2)城市下風(fēng)向觀察到的水汽摩爾分?jǐn)?shù)過量的原因,可能有以下3個方面,第一是人為的水汽來源,如燃燒和蒸發(fā)冷卻,近年來北京市推行的煤改氣政策可能起到了增加中心城區(qū)人為水汽排放的作用;第二是局地水體的貢獻(xiàn),如玉淵潭等城區(qū)水體;第三是中心城區(qū)和周邊地區(qū)融雪率和蒸發(fā)率的差異,在本研究中,這個差異較小。

    本研究觀測到的城市過量水汽排放發(fā)生在一年中蒸發(fā)速率比較低的冬季,而且,由于較低的氣溫,冬季的飽和水汽壓低于夏季,相對于噪聲和自然變率,城市水汽過量信號更加明顯。使用機(jī)載平臺在北京城市下風(fēng)方向?qū)^量的水汽摩爾分?jǐn)?shù)的研究,顯示的是全市范圍內(nèi)的水汽濃度,而不是局部的或者相對來說較小尺度的觀測結(jié)果。

    由于飛行航線不是專門為本研究設(shè)計(jì),因此存在可用數(shù)據(jù)量少、篩選困難等情況。如果要進(jìn)一步對該問題進(jìn)行研究,可設(shè)計(jì)專門的飛行試驗(yàn),進(jìn)一步擴(kuò)大航線范圍,使得背景濃度更能代表受城市影響較小的區(qū)域;在計(jì)算水汽排放率時也可以增加垂直方向上高度的變化等。另一方面,本文使用的飛機(jī)觀測數(shù)據(jù)是從北京市氣象局人工影響天氣辦公室日常飛機(jī)探測的數(shù)據(jù)中篩選出來的,沒有設(shè)計(jì)和進(jìn)行專門的飛行觀測試驗(yàn),成本較低,可操作性強(qiáng),為在國內(nèi)不同城市開展城市尺度人為潛熱的計(jì)算提供了新的思路和方法,最終的計(jì)算結(jié)果也為城市氣象環(huán)境數(shù)值模擬中人為潛熱的設(shè)置提供了一個參考。

    猜你喜歡
    邊界層人為過量
    過量食水果會加速衰老
    中老年保健(2021年7期)2021-08-22 07:43:44
    基于HIFiRE-2超燃發(fā)動機(jī)內(nèi)流道的激波邊界層干擾分析
    B3M4 Sandstorms in Asia Teaching Plan
    山高人為峰
    請勿過量飲酒
    吃糖過量也會“醉”?
    遵義(2018年15期)2018-08-28 12:20:14
    源正泉自清 山高人為峰
    中國篆刻(2017年5期)2017-07-18 11:09:30
    一類具有邊界層性質(zhì)的二次奇攝動邊值問題
    非特征邊界的MHD方程的邊界層
    山高人為峰
    詩歌月刊(2014年1期)2014-03-11 17:26:03
    大又大粗又爽又黄少妇毛片口| 噜噜噜噜噜久久久久久91| 人妻制服诱惑在线中文字幕| 乱人伦中国视频| 亚洲欧美中文字幕日韩二区| 免费大片18禁| 少妇被粗大猛烈的视频| 久久国产精品大桥未久av | 天美传媒精品一区二区| 纵有疾风起免费观看全集完整版| 国产日韩欧美亚洲二区| 国产精品女同一区二区软件| 国产精品福利在线免费观看| 少妇高潮的动态图| av在线app专区| 国产高清有码在线观看视频| 少妇丰满av| h日本视频在线播放| 在线免费观看不下载黄p国产| 免费av中文字幕在线| 久久99蜜桃精品久久| 插逼视频在线观看| 国产淫语在线视频| 蜜桃久久精品国产亚洲av| 久久综合国产亚洲精品| 亚洲综合色惰| 91久久精品电影网| 十八禁高潮呻吟视频 | 亚洲欧美日韩东京热| 中文字幕人妻熟人妻熟丝袜美| 国产黄色视频一区二区在线观看| 日日爽夜夜爽网站| 在线观看av片永久免费下载| 看免费成人av毛片| 精品人妻一区二区三区麻豆| 国产成人aa在线观看| 国产精品三级大全| 99久久中文字幕三级久久日本| 久久久午夜欧美精品| 人人妻人人澡人人爽人人夜夜| 国产在线男女| 国产精品一二三区在线看| 男人和女人高潮做爰伦理| 午夜福利视频精品| √禁漫天堂资源中文www| 精品人妻偷拍中文字幕| 国产白丝娇喘喷水9色精品| 欧美97在线视频| 成年av动漫网址| 亚洲精品国产av成人精品| 黄色怎么调成土黄色| 亚洲av电影在线观看一区二区三区| 少妇人妻 视频| 婷婷色麻豆天堂久久| 国产黄色免费在线视频| 人人妻人人澡人人爽人人夜夜| 熟女av电影| 99热全是精品| 最近手机中文字幕大全| 亚洲欧美日韩卡通动漫| 搡老乐熟女国产| 丰满人妻一区二区三区视频av| 久久99蜜桃精品久久| av天堂中文字幕网| 欧美人与善性xxx| 久久毛片免费看一区二区三区| 丝袜喷水一区| 狠狠精品人妻久久久久久综合| 一级黄片播放器| 久久女婷五月综合色啪小说| 国产精品麻豆人妻色哟哟久久| 国产69精品久久久久777片| 国产精品99久久99久久久不卡 | 在线观看免费视频网站a站| 新久久久久国产一级毛片| 纯流量卡能插随身wifi吗| 韩国高清视频一区二区三区| 丝袜脚勾引网站| 2021少妇久久久久久久久久久| 久久午夜福利片| 免费观看性生交大片5| 丰满人妻一区二区三区视频av| 精品国产乱码久久久久久小说| 中文天堂在线官网| 熟女电影av网| 久久久久国产网址| a级一级毛片免费在线观看| 99久久人妻综合| 亚洲精品日本国产第一区| 观看美女的网站| 日本午夜av视频| 日韩人妻高清精品专区| 人人妻人人澡人人爽人人夜夜| 国产一区二区三区综合在线观看 | 国产精品福利在线免费观看| 日本-黄色视频高清免费观看| 欧美激情国产日韩精品一区| 男女无遮挡免费网站观看| 欧美成人午夜免费资源| 亚洲av综合色区一区| 简卡轻食公司| 18禁在线播放成人免费| 日韩中字成人| 边亲边吃奶的免费视频| 人妻制服诱惑在线中文字幕| 在线精品无人区一区二区三| 欧美变态另类bdsm刘玥| 免费黄频网站在线观看国产| 国产成人精品一,二区| 晚上一个人看的免费电影| 亚洲精品第二区| 免费黄色在线免费观看| 亚洲四区av| 一本大道久久a久久精品| 99九九在线精品视频 | 久久久欧美国产精品| 久久久国产欧美日韩av| 久久精品国产亚洲av天美| 亚洲精品中文字幕在线视频 | 欧美 日韩 精品 国产| 亚洲av不卡在线观看| 人妻系列 视频| 国产 精品1| 日本欧美国产在线视频| 王馨瑶露胸无遮挡在线观看| 一级毛片我不卡| 卡戴珊不雅视频在线播放| 国产成人a∨麻豆精品| 中文字幕av电影在线播放| 久久青草综合色| 亚洲欧美日韩卡通动漫| 国产成人免费无遮挡视频| 国产精品久久久久久av不卡| 免费观看无遮挡的男女| 免费观看性生交大片5| 国产精品人妻久久久影院| 日韩不卡一区二区三区视频在线| 亚洲欧美精品专区久久| 国产黄频视频在线观看| 国产欧美日韩精品一区二区| 精品亚洲成a人片在线观看| 精品亚洲成a人片在线观看| av在线app专区| 精品亚洲乱码少妇综合久久| 精品亚洲乱码少妇综合久久| 99热这里只有是精品50| 在线观看一区二区三区激情| 成人美女网站在线观看视频| 亚洲丝袜综合中文字幕| 两个人免费观看高清视频 | av又黄又爽大尺度在线免费看| 日日撸夜夜添| 国产精品国产三级国产av玫瑰| 久久久久视频综合| 日韩三级伦理在线观看| 在线观看三级黄色| 免费久久久久久久精品成人欧美视频 | 激情五月婷婷亚洲| 日本爱情动作片www.在线观看| 亚洲av福利一区| 亚洲第一区二区三区不卡| 免费在线观看成人毛片| 欧美日韩国产mv在线观看视频| 一二三四中文在线观看免费高清| 黄色毛片三级朝国网站 | 国产精品国产av在线观看| 亚洲国产日韩一区二区| 大香蕉久久网| 女性生殖器流出的白浆| 老司机影院毛片| 免费观看在线日韩| 亚洲一区二区三区欧美精品| 如日韩欧美国产精品一区二区三区 | 亚洲精品国产成人久久av| 久久午夜综合久久蜜桃| 综合色丁香网| 国产精品久久久久久av不卡| 亚洲成人手机| 最近手机中文字幕大全| 下体分泌物呈黄色| 国产免费一区二区三区四区乱码| 插阴视频在线观看视频| 天堂8中文在线网| 视频中文字幕在线观看| 色视频在线一区二区三区| 国产高清有码在线观看视频| 一区二区三区乱码不卡18| 99九九在线精品视频 | tube8黄色片| 日韩欧美一区视频在线观看 | 高清av免费在线| 在线观看美女被高潮喷水网站| 亚洲婷婷狠狠爱综合网| 熟女人妻精品中文字幕| 久久久久久久大尺度免费视频| 日本91视频免费播放| 国产精品国产av在线观看| 在线天堂最新版资源| 国产精品成人在线| www.色视频.com| 国产亚洲午夜精品一区二区久久| 亚洲情色 制服丝袜| 久久国产亚洲av麻豆专区| 性色avwww在线观看| 国产中年淑女户外野战色| 嘟嘟电影网在线观看| 亚洲精品乱码久久久v下载方式| 精品视频人人做人人爽| 一个人看视频在线观看www免费| 亚洲av二区三区四区| 国产永久视频网站| 国产黄频视频在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | 蜜桃在线观看..| 国产在线视频一区二区| 三级国产精品欧美在线观看| av免费观看日本| 免费久久久久久久精品成人欧美视频 | 国产精品免费大片| 91精品国产国语对白视频| 简卡轻食公司| 高清毛片免费看| 国产精品人妻久久久久久| 成人免费观看视频高清| 男女边摸边吃奶| 日韩中字成人| 18禁裸乳无遮挡动漫免费视频| 久久久国产一区二区| 偷拍熟女少妇极品色| 国产精品久久久久久av不卡| 亚洲国产精品一区二区三区在线| av福利片在线| 黄色毛片三级朝国网站 | 亚洲,欧美,日韩| 欧美xxⅹ黑人| 中文欧美无线码| 日韩制服骚丝袜av| √禁漫天堂资源中文www| 精品少妇久久久久久888优播| 夜夜骑夜夜射夜夜干| 免费黄色在线免费观看| 久久久久久人妻| 精品一品国产午夜福利视频| 人人妻人人爽人人添夜夜欢视频 | 日本av免费视频播放| 亚洲精品色激情综合| 久久ye,这里只有精品| 人妻人人澡人人爽人人| 性色av一级| 一个人免费看片子| 久久久久久久久久成人| av福利片在线观看| 国产黄色免费在线视频| 26uuu在线亚洲综合色| 亚洲国产欧美在线一区| 一区二区三区四区激情视频| 久久人人爽人人片av| 午夜精品国产一区二区电影| 中文在线观看免费www的网站| 尾随美女入室| 我的老师免费观看完整版| 日韩中字成人| 国产视频首页在线观看| 春色校园在线视频观看| 在线亚洲精品国产二区图片欧美 | 中文字幕人妻熟人妻熟丝袜美| 亚洲欧美成人精品一区二区| 国产精品久久久久久久电影| 成人黄色视频免费在线看| 亚洲欧洲日产国产| 一级毛片我不卡| 国产爽快片一区二区三区| 久久毛片免费看一区二区三区| 亚洲经典国产精华液单| 欧美日韩视频精品一区| 一边亲一边摸免费视频| 美女内射精品一级片tv| 免费在线观看成人毛片| 日韩电影二区| 精品国产露脸久久av麻豆| 五月伊人婷婷丁香| 国产高清国产精品国产三级| 亚洲欧美清纯卡通| 久久久久久人妻| 麻豆成人av视频| 午夜日本视频在线| 爱豆传媒免费全集在线观看| 嫩草影院新地址| 亚洲欧美中文字幕日韩二区| 老女人水多毛片| 亚洲中文av在线| 久久99精品国语久久久| 国产国拍精品亚洲av在线观看| 18禁在线无遮挡免费观看视频| 亚洲精品成人av观看孕妇| tube8黄色片| 日韩熟女老妇一区二区性免费视频| 91久久精品电影网| 一级毛片 在线播放| 嫩草影院新地址| 一级黄片播放器| 两个人免费观看高清视频 | 成人漫画全彩无遮挡| 熟女人妻精品中文字幕| av黄色大香蕉| h日本视频在线播放| 国产极品粉嫩免费观看在线 | 欧美精品人与动牲交sv欧美| 美女国产视频在线观看| 日日摸夜夜添夜夜爱| 超碰97精品在线观看| 日韩欧美一区视频在线观看 | 观看美女的网站| 亚洲精品第二区| 国产淫语在线视频| 中文字幕人妻丝袜制服| 免费播放大片免费观看视频在线观看| 欧美性感艳星| 三级国产精品片| 在线观看免费视频网站a站| 高清午夜精品一区二区三区| 亚洲欧美成人综合另类久久久| 日韩成人av中文字幕在线观看| 免费看av在线观看网站| 久久亚洲国产成人精品v| 日韩av不卡免费在线播放| 中文精品一卡2卡3卡4更新| 亚洲欧美清纯卡通| av黄色大香蕉| 全区人妻精品视频| 在线观看三级黄色| 亚洲精品国产av成人精品| 极品少妇高潮喷水抽搐| av福利片在线观看| 少妇人妻 视频| 日韩欧美精品免费久久| 高清毛片免费看| 欧美日本中文国产一区发布| 免费观看a级毛片全部| 久久 成人 亚洲| 欧美成人午夜免费资源| av有码第一页| 精品国产一区二区久久| 久久国产亚洲av麻豆专区| 久久久久久久久久久丰满| 岛国毛片在线播放| 国产免费一级a男人的天堂| av一本久久久久| 水蜜桃什么品种好| 亚洲精华国产精华液的使用体验| 亚洲欧美一区二区三区黑人 | 国产精品偷伦视频观看了| 成人漫画全彩无遮挡| 国产亚洲一区二区精品| 99国产精品免费福利视频| 亚洲精品自拍成人| 中文精品一卡2卡3卡4更新| 日本黄色片子视频| 人人妻人人看人人澡| 欧美精品一区二区免费开放| 国产黄片视频在线免费观看| 亚洲精品456在线播放app| 国产成人freesex在线| 亚洲av成人精品一二三区| 国产成人a∨麻豆精品| 涩涩av久久男人的天堂| 国产免费福利视频在线观看| 大片电影免费在线观看免费| 在线 av 中文字幕| 午夜影院在线不卡| 亚洲电影在线观看av| 黑人高潮一二区| av不卡在线播放| 观看免费一级毛片| 国产免费一级a男人的天堂| av又黄又爽大尺度在线免费看| 人妻人人澡人人爽人人| 十八禁高潮呻吟视频 | 天美传媒精品一区二区| 成年av动漫网址| 国产成人freesex在线| 亚洲激情五月婷婷啪啪| 久久久精品94久久精品| 国产乱来视频区| 内射极品少妇av片p| 看非洲黑人一级黄片| 国产综合精华液| 亚洲av综合色区一区| 国产成人freesex在线| 久久精品久久久久久噜噜老黄| 一区二区三区四区激情视频| 美女视频免费永久观看网站| 五月开心婷婷网| 国产精品人妻久久久影院| 在线免费观看不下载黄p国产| 91精品一卡2卡3卡4卡| 日韩精品有码人妻一区| 亚洲综合精品二区| 国产精品麻豆人妻色哟哟久久| 一级,二级,三级黄色视频| 天美传媒精品一区二区| 高清在线视频一区二区三区| 国产成人精品无人区| 亚洲情色 制服丝袜| 久久久久久久久久久丰满| 国产精品久久久久久久电影| 少妇的逼好多水| 久久99热6这里只有精品| 国产探花极品一区二区| 国产一区二区在线观看av| 大香蕉97超碰在线| 日韩视频在线欧美| 亚洲美女黄色视频免费看| 中文字幕av电影在线播放| 国产成人精品一,二区| 国产 一区精品| 一区二区三区乱码不卡18| 一级爰片在线观看| 婷婷色麻豆天堂久久| 国产成人a∨麻豆精品| 男女无遮挡免费网站观看| 亚洲欧美精品专区久久| 亚洲精品国产av成人精品| 日韩一区二区视频免费看| 最黄视频免费看| 在线 av 中文字幕| 日日摸夜夜添夜夜添av毛片| 视频区图区小说| 亚洲精品第二区| 三级经典国产精品| 国产亚洲午夜精品一区二区久久| 国产亚洲一区二区精品| 国产黄片美女视频| 一边亲一边摸免费视频| 成人二区视频| 三级经典国产精品| 国国产精品蜜臀av免费| 熟妇人妻不卡中文字幕| 你懂的网址亚洲精品在线观看| 一边亲一边摸免费视频| 亚洲综合精品二区| 日本欧美国产在线视频| 99九九在线精品视频 | 天堂俺去俺来也www色官网| 国产精品国产三级国产专区5o| 少妇 在线观看| 国产淫语在线视频| 久久99一区二区三区| 久久99热这里只频精品6学生| 国产高清三级在线| 国产午夜精品久久久久久一区二区三区| 91精品伊人久久大香线蕉| 亚洲第一av免费看| 国产精品一区二区性色av| 久久人人爽av亚洲精品天堂| 免费黄频网站在线观看国产| 街头女战士在线观看网站| 插阴视频在线观看视频| 亚洲无线观看免费| 国产欧美另类精品又又久久亚洲欧美| 免费看日本二区| 国产女主播在线喷水免费视频网站| 在线观看美女被高潮喷水网站| 日本91视频免费播放| 国产精品无大码| 国产精品久久久久久精品电影小说| 亚洲中文av在线| 欧美少妇被猛烈插入视频| av黄色大香蕉| 自线自在国产av| 亚洲av国产av综合av卡| 丰满少妇做爰视频| 精品人妻一区二区三区麻豆| 男人舔奶头视频| 一本大道久久a久久精品| 人人妻人人爽人人添夜夜欢视频 | 国产精品免费大片| 内地一区二区视频在线| 亚洲av中文av极速乱| 国内精品宾馆在线| 亚洲成色77777| 一级毛片久久久久久久久女| 国产亚洲av片在线观看秒播厂| 久久国产乱子免费精品| 欧美日韩av久久| 少妇高潮的动态图| 噜噜噜噜噜久久久久久91| 高清在线视频一区二区三区| 国产精品成人在线| a 毛片基地| 亚洲无线观看免费| 高清av免费在线| 九九爱精品视频在线观看| 一级,二级,三级黄色视频| 插阴视频在线观看视频| 亚洲欧美日韩东京热| 国产亚洲午夜精品一区二区久久| 丝瓜视频免费看黄片| 亚洲欧美精品专区久久| 中文乱码字字幕精品一区二区三区| 成人免费观看视频高清| 欧美少妇被猛烈插入视频| 久久精品国产a三级三级三级| 国产永久视频网站| videossex国产| 欧美日本中文国产一区发布| 国产亚洲欧美精品永久| 97在线人人人人妻| 午夜av观看不卡| 最近的中文字幕免费完整| 免费观看av网站的网址| av在线播放精品| 亚洲国产精品999| 亚洲欧洲日产国产| 国产欧美日韩一区二区三区在线 | 人体艺术视频欧美日本| 人人澡人人妻人| 久久婷婷青草| 一本色道久久久久久精品综合| 久久久久久人妻| 国产午夜精品久久久久久一区二区三区| 人妻人人澡人人爽人人| 精品酒店卫生间| 亚洲欧洲日产国产| 久久久久久久亚洲中文字幕| 99热国产这里只有精品6| 免费在线观看成人毛片| 日韩av在线免费看完整版不卡| 亚洲精品国产av成人精品| 国产在线免费精品| 在线观看美女被高潮喷水网站| 少妇被粗大的猛进出69影院 | 美女内射精品一级片tv| 成年女人在线观看亚洲视频| 亚洲一级一片aⅴ在线观看| 晚上一个人看的免费电影| 亚洲精品,欧美精品| 亚洲欧美一区二区三区国产| 老司机影院毛片| 嫩草影院新地址| 尾随美女入室| 久久久国产一区二区| 观看免费一级毛片| 亚洲国产精品专区欧美| 纯流量卡能插随身wifi吗| 亚洲精品日韩av片在线观看| 秋霞在线观看毛片| 在线观看av片永久免费下载| 中文字幕精品免费在线观看视频 | 亚洲婷婷狠狠爱综合网| 女性生殖器流出的白浆| 日日爽夜夜爽网站| 亚洲在久久综合| 久久午夜综合久久蜜桃| 赤兔流量卡办理| 亚洲欧洲日产国产| 亚洲三级黄色毛片| 亚洲自偷自拍三级| 亚洲国产毛片av蜜桃av| 亚洲av男天堂| 国产精品伦人一区二区| av在线观看视频网站免费| 最近的中文字幕免费完整| 91午夜精品亚洲一区二区三区| 丝袜脚勾引网站| 交换朋友夫妻互换小说| 男人和女人高潮做爰伦理| 精品熟女少妇av免费看| 最近手机中文字幕大全| 在线天堂最新版资源| 性色avwww在线观看| 成人亚洲欧美一区二区av| 大片免费播放器 马上看| 国产精品一区二区在线观看99| 啦啦啦视频在线资源免费观看| 人人妻人人爽人人添夜夜欢视频 | 亚洲欧美一区二区三区国产| 日韩欧美 国产精品| 日韩亚洲欧美综合| 永久免费av网站大全| 男人爽女人下面视频在线观看| 欧美xxⅹ黑人| 在线 av 中文字幕| 成人毛片60女人毛片免费| 成年人午夜在线观看视频| 国产女主播在线喷水免费视频网站| 亚洲精品日本国产第一区| 伊人亚洲综合成人网| 久久精品夜色国产| 亚洲精品自拍成人| 免费在线观看成人毛片| 国产黄频视频在线观看| 婷婷色av中文字幕| 美女主播在线视频| 九九爱精品视频在线观看| 亚洲伊人久久精品综合| 观看免费一级毛片| 午夜激情久久久久久久| 99久久精品一区二区三区| 精品午夜福利在线看| 在线观看免费视频网站a站| 欧美少妇被猛烈插入视频| 亚洲欧美成人精品一区二区| 观看美女的网站| 国产一区二区三区综合在线观看 | 伦理电影大哥的女人| 久久午夜福利片| 夜夜看夜夜爽夜夜摸| 国产精品久久久久久精品电影小说| 两个人的视频大全免费| 精品一区二区三区视频在线| 91午夜精品亚洲一区二区三区| 免费在线观看成人毛片| 街头女战士在线观看网站| 一级av片app|