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

    星載激光測風技術體制及其評估方法綜述

    2022-07-06 07:33:14孫學金張傳亮方樂陸文趙世軍葉松
    遙感學報 2022年6期
    關鍵詞:風場體制大氣

    孫學金,張傳亮,,方樂,陸文,趙世軍,葉松

    1.國防科技大學氣象海洋學院,南京 211101;

    2.61741部隊,北京 100081

    1 引言

    風是描述大氣狀態(tài)的基本參量之一,全球范圍內(nèi)大氣風場的探測對于研究中尺度動力學過程、大氣輸送和交換作用以及極端天氣的預報預防等具有重要意義(Baker 等,1995)。從20世紀70年代起,人們就開始對風場衛(wèi)星探測技術進行研究。1971年Endlich 等(1971)通過提取云的特征量實現(xiàn)了衛(wèi)星云導風的探測。Leese 等(1971)通過引入快速傅里葉變化技術,提高了云導風的計算速度。許健民等(1997)對云導風計算的相關問題開展了研究。王振會和曾維麟(1996)、朱平等(2007)、龍智勇等(2011)、馬俠霖等(2014)進一步發(fā)展了云導風的處理算法,提高了云導風的計算效率和業(yè)務化應用。目前云導風的探測技術已經(jīng)比較成熟,然而該技術只能實現(xiàn)有云大氣層的單層風信息的提取。同時,人們也發(fā)展了利用微波輻射計、微波散射計獲取海表風場的技術。1978年,美國發(fā)射的SeaSAT 裝載了微波輻射計載荷,初步實現(xiàn)了海面風速的遙感(Njoku等,1980)。后來人們在傳統(tǒng)微波輻射計的基礎上,引入了對于電磁波第三、第四Stokes分量的測量,使得探測海面風速風向成為可能。2003年,第一臺星載全極化微波輻射計WindSat 發(fā)射升空,成功實現(xiàn)了海表風速和風向的遙感(Gaiser 等,2004)。歐洲空間局研制的SMOS(Soil Moisture and Ocean Salinity)也搭載了極化微波輻射計,實現(xiàn)了海面風場的探測(Barre 等,2008)。在國內(nèi),中國航天科技集團公司第504 研究所(王曉海和李浩,2005)、解放軍理工大學(施健康等,2009;陸文等,2010;嚴衛(wèi)等,2010)等單位也對全極化微波輻射計的設計和應用展開了深入研究,并取得了豐碩成果。然而,星載微波輻射計只能實現(xiàn)全球海表風場的探測,無法實現(xiàn)空中風場的廓線探測。在星載激光測風雷達發(fā)射之前,探測空中風場廓線的主要手段包括無線電探空儀和風廓線雷達,然而無線電探空儀站點主要分布在北半球的陸地地區(qū)(Atlas等,1985),南半球、赤道地區(qū)以及海洋分布較少,無法實現(xiàn)全球風場廓線的有效和持續(xù)探測。世界氣象組織將全球三維風場的探測視為21世紀最具挑戰(zhàn)性和最有意義的氣象觀測任務之一(WMO,2001)。

    通過對不同衛(wèi)星測風技術的研究,星載激光測風雷達被認為是探測全球風場的最佳手段。Abreu(1979)首次提出在衛(wèi)星或者空間站上搭載激光測風雷達進行全球風場探測的概念,并研究了雷達參數(shù)對測風精度的影響。1985年,美國啟動了大氣風場探測計劃LAWS(Laser Atmospheric Wind Sounder),開展了2 μm 相干激光多普勒雷達的相關技術研究,該項目由于受制于當時技術條件的限制而被迫擱淺(Beranek 等,1989)。經(jīng)過20 多年的探索,美國于2007年再次啟動星載激光雷達測風計劃:全球風場探測系統(tǒng)GWOS(Global Wind Observing System),采用直接探測和相干探測的混合技術體制的多普勒激光測風雷達HDWL(Hybrid Doppler Wind Lidar)實現(xiàn)全球風場的探測,目前該項目已推遲發(fā)射(Baker,2008)。歐空局對于星載激光測風雷達的研究也起步較早,可追溯至20世紀80年代,并于1999年正式啟動大氣動力學任務ADM-Aeolus ADM (Atmospheric Dynamics Mission),并計劃搭載基于直接探測體制的星載激光測風雷達ALADIN(Atmospheric LAser Doppler INstrument)來實現(xiàn)全球風場的探測。最早計劃于2007年發(fā)射,后因激光器和望遠鏡等硬件方面的技術難題幾經(jīng)推遲(Wernham 等,2016)。2018年8月22日(世界時),Aeolus 衛(wèi)星順利升空,對于人類全球風場廓線的探測具有里程碑式的意義。1998年日本開始研制運行于國際空間站上的相干激光多普勒測風雷達JEM/CDL(Japanese Experiment Module/Coherent Doppler Lidar)(Itabe等,2001)。2011年日本又開展了運行于獨立衛(wèi)星軌道的相干多普勒激光測風雷達CDWL(Coherent Doppler Wind Lidar)的研究,該技術體制可實現(xiàn)邊界層和云場內(nèi)部風場的探測(Baron 等,2017;Ishii等,2017)。

    歐洲空間局、美國和日本都投入了大量人力、物力、財力積極探索和研制星載激光測風雷達,提出了3種不同的技術體制。本文通過大量調(diào)研相關文獻簡要概述了歐洲空間局、美國和日本的星載激光測風雷達技術體制,并從數(shù)據(jù)獲取率和獲取量、風場探測精度以及對數(shù)值預報系統(tǒng)的改善效果3個方面對這3種技術體制進行了評估和論證。

    2 星載激光測風技術體制

    本節(jié)分別對歐洲空間局Aeolus、美國HDWL和日本CDWL 研制的星載激光測風技術體制進行介紹。

    2.1 歐洲空間局Aeolus技術體制

    歐洲空間局研制的Aeolus 星載激光測風技術體制如圖1所示(Marseille 和Stoffelen,2003;Reitebuch 等,2006)。衛(wèi)星擬運行在400 km 高度的太陽同步晨昏軌道上,以減小太陽背景輻射對測風精度的影響(Zhang 等,2019)。后來為了增加回波信號的信噪比,將軌道高度調(diào)整至320 km。衛(wèi)星觀測天頂角為35°。Aeolus 激光雷達工作在355 nm 波段,采用直接探測技術分別通過Rayleigh通道和Mie 通道進行大氣風場的探測,其中Rayleigh 通道主要以大氣分子為示蹤物探測對流層和平流層下層的大氣風場,Mie 通道主要以氣溶膠和云滴粒子為示蹤物探測邊界層和云內(nèi)風場。歐洲空間局對Aeolus 風場探測的垂直分辨率和精度要求如表1所示。為了提高測風精度,Aeolus 探測單元需累積若干個激光束后再進行徑向風速的反演,其中一次測量(measurement)過程累積的激光脈沖為20個,一次觀測(observation)是由30次測量組成的,造成的一次觀測過程的水平分辨率約為90 km。

    表1 歐洲空間局對星載激光雷達測風不確定度的要求Table 1 The observation accuracy expectation of ESA for spaceborne DWL

    圖1 Aeolus激光測風技術體制(Reitebuch等,2006)Fig.1 The technical system of spaceborne DWL Aeolus(Reitebuch et al.,2006)

    Aeolus 的系統(tǒng)結構如圖2所示,系統(tǒng)主要由激光器、發(fā)射接收系統(tǒng)和鑒頻系統(tǒng)組成,其詳細工作流程參見張日偉等(2014),本文進行簡要描述。首先,激光器產(chǎn)生355 nm 的激光,并發(fā)射到大氣中,大氣中分子、氣溶膠或云滴粒子通過后向散射作用將部分激光信號散射回激光雷達。氣溶膠和云滴粒子的后向散射信號進入到Mie通道鑒頻器中,即Fizeau干涉儀。不同頻率的回波信號會在Fizeau干涉儀不同位置產(chǎn)生干涉條紋。鑒頻系統(tǒng)通過條紋形成的位置反演氣溶膠和云滴粒子回波信號的多普勒頻移。大氣分子的后向散射信號會進入到Rayleigh 通道鑒頻器中,Rayleigh 通道鑒頻器由2 個F-P 標準具組成。回波信號一部分進入F-P 標準具1,并在相應位置處的CCD 探測器上形成亮斑,另外一部分信號反射進入到F-P 標準具2中,在CCD探測器上形成另外一個亮斑。不同頻率的信號產(chǎn)生的兩個亮斑的相對亮度會產(chǎn)生變化,鑒頻系統(tǒng)通過兩個亮斑之間的相對亮度差異反演大氣分子回波信號的多普勒頻移。歐洲空間局完成了星載激光測風雷達關鍵部件,包括激光二極管、發(fā)射系統(tǒng)以及光子探測器等的主要技術攻關工作,并為星載激光雷達部件設定了相應的技術標準(Zahir和Durand,2011)。

    圖2 星載激光多普勒測風雷達系統(tǒng)結構圖(張日偉等,2014)Fig.2 The system structure diagram of spaceborne DWL(Zhang et al.,2014)

    2.2 美國HDWL技術體制

    在美國HDWL 技術體制中,衛(wèi)星擬運行于400 km 高度的太陽同步晨昏軌道。該技術體制采用相干探測技術和非相干探測技術相結合的混合多普勒激光測風雷達方案(Emmitt,2004),其相干探測技術采用的波長為2 μm,非相干探測技術采用的波長為355 nm。與Aeolus不同的是,HDWL可探測多個方向的徑向風速值。在不采用望遠鏡掃描技術的情況下,一個望遠鏡只能探測一個方向的徑向風速值。而望遠鏡的體積一般較大,以Aeolus 為例,其望遠鏡直徑達1.5 m,如果采用掃描技術的話,會給衛(wèi)星的姿態(tài)控制帶來新的挑戰(zhàn)。為了探測多個方向的徑向風速值,HDWL 擬安裝4個不同觀測方向的望遠鏡,其方位角分別為45°、135°、225°和315°,其觀測天頂角為45°,如圖3所示。Marx等(2013)已經(jīng)在實驗室研制出來由4個望遠鏡組成的簡潔、輕量化的星載激光測風雷達發(fā)射接收系統(tǒng),并且每個望遠鏡都可以接收355 nm和2 μm 的后向散射回波信號。Masutani 等(2010)基于衛(wèi)星軌道仿真軟件對HDWL 軌道進行了模擬,得到了激光束與地表的交叉點軌跡。模擬結果表示這種觀測方式保證了每個探測點都會以不同的觀測角度被先后探測兩次,從而實現(xiàn)了不同方向徑向風速的探測。在進行風速合成時,假設垂直方向的風速分量為0 m·s-1,即可實現(xiàn)二維水平風矢量的探測。

    圖3 HDWL觀測體制示意圖(Masutani等,2010)Fig.3 The observation geometry of HDWL(Masutani et al.,2010)

    2.3 日本CDWL技術體制

    日本CDWL計劃是從ISS/JE計劃發(fā)展而來,二者技術體制比較相似,CDWL 計劃在ISS/JE 計劃的基礎上進行了部分改動。ISS/JE星載激光測風雷達觀測方式示意圖如圖4所示(Ishii等,2012)。ISS/JE 擬在國際空間站上運行,采用相干探測體制,該激光雷達包含兩個望遠鏡,其方位角分別為45°和135°,觀測天頂角為30°。CDWL 運行于獨立軌道,軌道高度為220 km,并將兩個望遠鏡的觀測天頂角調(diào)整為35°。其觀測方式與美國HDWL 觀測方式類似,其前向望遠鏡和后向望遠鏡將分別探測同一探測點的徑向風場,如圖5所示。假設垂直風速為0 m·s-1,可反演大氣風場的二維水平風矢量。

    圖4 ISS/JE觀測體制示意圖(Ishii等,2012)Fig.4 The observation geometry of ISS/JE(Ishii et al.,2012)

    圖5 CDWL觀測方式示意圖(Ishii等,2017)Fig.5 The observation geometry of CDWL(Ishii et al.,2017)

    3 星載激光測風雷達技術體制的評估

    星載激光測風雷達研發(fā)費用比較高昂,所以在發(fā)射之前必須對采用的技術體制進行充分論證。針對不同的目的,需要從不同方面進行論證。目前對于星載激光測風雷達主要的評估方法有基于數(shù)據(jù)獲取率和獲取量的評估、基于測風精度的評估以及基于對數(shù)值預報系統(tǒng)改善效果的評估,本節(jié)將對以上3種評估方法進行分別論述。

    3.1 基于數(shù)據(jù)獲取率和獲取量的評估

    全球云的分布是影響星載激光測風雷達數(shù)據(jù)獲取率的關鍵因素,云滴粒子是星載激光測風雷達探測云內(nèi)風場的主要示蹤物,但由于激光在云中的穿透能力有限,云也是限制星載激光測風雷達對云層以下大氣探測能力的重要因子。

    目前的星載激光測風雷達可利用分子的Rayleigh 散射信號進行大氣風場的探測,也可利用云和氣溶膠粒子的Mie散射信號對大氣風場進行探測。在晴空大氣下,在大部分高度層,星載激光測風雷達主要以分子為示蹤物進行大氣風場的探測,在邊界層,氣溶膠數(shù)密度較大,星載激光測風雷達主要以氣溶膠為示蹤物進行大氣風場的探測。在有云大氣下,星載激光測風雷達可利用云粒子的Mie散射信號探測云內(nèi)風場,云粒子數(shù)密度越大,Mie 散射信號越強,探測的風速越準確,但云粒子同時對激光具有較強的衰減作用,激光在云中的穿透能力取決于單束激光能量和云的光學厚度,如果云的光學厚度較大,星載激光測風雷達則無法探測云層底部以及云層以下的大氣風場,所以星載激光測風雷達的數(shù)據(jù)獲取率與云場的分布情況有關。在評估過程中,全球云場模型是否能夠代表真實的大氣狀況至關重要。另外,由于云的分布具有明顯的日變化特征(Wood等,2002;Eastman 和Warren,2014;Min 和Zhang,2014),所以星載激光測風雷達的數(shù)據(jù)獲取率還與衛(wèi)星軌道的升交點有關。同時,數(shù)據(jù)獲取率還與激光測風雷達的探測技術有關,以氣溶膠和云粒子為示蹤物的相干探測技術可獲取約30%的全球風場數(shù)據(jù),以分子為示蹤物的直接探測技術可獲取約70%的全球風場數(shù)據(jù)(Masutani等,2010)。

    與數(shù)據(jù)獲取率相關的量為數(shù)據(jù)獲取量,本文中的數(shù)據(jù)獲取量是指是指徑向風速的探測量。假設HDWL 和CDWL 每個望遠鏡都對應一套獨立的激光發(fā)射與接收單元,則在激光發(fā)射頻率一致的條件下,HDWL 和CDWL 測得的徑向風速廓線的數(shù)量應該分別是Aeolus 的4 倍和2 倍。Masutani 等(2010)的評估實驗也驗證了這一點,他們基于ISCCP(International Satellite Cloud Climatology Project)全球云場的分布數(shù)據(jù)庫,在相同的軌道參數(shù)下,即高度為400 km 的太陽同步晨昏軌道,分別模擬了基于Aeolus、HDWL以及CDWL星載激光測風雷達技術體制的衛(wèi)星運行情況,并統(tǒng)計了一個月內(nèi)探測到的measurement 數(shù)據(jù)量情況,如圖6所示。由圖可得,3 種技術體制的星載激光測風雷達數(shù)據(jù)獲取量存在較大的差異。其中以分子和氣溶膠/云滴粒子為示蹤物進行大氣風場探測的Aeolus 和HDWL 能夠探測邊界層至對流層的大氣風場。美國HDWL 體制探測的數(shù)據(jù)量最多,其數(shù)據(jù)獲取量約為Aeolus的4倍,在對流層中下層大氣約為日本CDWL 體制的2 倍;由于使用相干探測技術,CDWL體制無法探測對流層上層大氣風場情況,只能獲取邊界層至對流層下層的大氣風場,在該高度層內(nèi),其數(shù)據(jù)獲取量約為Aeolus的2倍。

    圖6 1個月內(nèi)3種星載激光測風雷達體制的有效探測數(shù)量格點分辨率為2.5°×2.5°(Masutani等,2010)Fig.6 The number of measurements obtained from the three technical systems of spaceborne DWLs in a 2.5 grid box(Masutani et al.,2010)

    3.2 基于風場探測精度的評估

    影響星載激光測風雷達測風精度的主要因素包括泊松噪聲和大氣不均勻性。泊松噪聲,也被稱為光子噪聲,主要由于光學探測器的量子特性導致的,在星載激光測風雷達探測器將回波信號的能量轉化為電子的過程中,探測器上產(chǎn)生的電子數(shù)并不等于回波信號的能量所能激發(fā)的電子數(shù),而是服從泊松分布(Liu 等,2006;Hasinoff 等,2010;Zhang等,2018)。按照泊松分布的特性,方差等于平均值。如果回波信號的能量等于10000個電子對應的能量,則探測器上產(chǎn)生的電子數(shù)的標準差為100 個,其引起的相對不確定度為1%;如果回波信號的能量等于100個電子對應的能量,則探測器上產(chǎn)生電子數(shù)的標準差為10 個,其引起的相對不確定度為10%。所以回波信號能量越大,CCD 上產(chǎn)生電子數(shù)的相對不確定度就越小,反演的測風不確定度也就越小。由于Mie后向散射回波的能量一般高于Rayleigh 后向散射回波,所以利用Mie 后向散射回波的測風精度一般高于Rayleigh 后向散射回波的測風精度。一般來說,降低泊松噪聲對測風精度的影響的方法包括增加激光的發(fā)射能量或在探測器上將回波信號能量轉化為電子之前累積更多的信號能量,這即是Aeolus的一個“測量值”累積了20個激光束原因。激光發(fā)射能量的降低會增加測風的不確定度,Aeolus激光器單束激光能量原計劃為80 mJ,發(fā)射升空之后由于散熱不理想的原因,其最大能量為65 mJ,逐漸降低并最終穩(wěn)定至50 mJ,導致其Rayleigh通道測風不確定度約為4.75 m·s-1,相比歐洲空間局預期(如表1所示)存在一定差距(Kanitz等,2019;Reitebuch等,2019),但數(shù)值預報系統(tǒng)實驗表明,該風場探測結果仍然能夠較好地改善預報精度(Straume等,2019)。

    影響Aeolus 測風不確定度最直接的因素是信噪比,Marseille 和Stoffelen(2003)在計算Aeolus的信噪比時,考慮了泊松噪聲的特點,計算公式為

    式中,SNRMie(z)和SNRRay(z)分別代表在高度z處Mie 通道和Rayleigh 通道的信噪比。對于Mie 通道,NMie為Mie 通道探測器探測到的總的電子數(shù),NB_Ray為由于串擾作用,Mie 通道接收到的Rayleigh 回波信號激發(fā)的電子數(shù),NB_BKG為由太陽背景輻射能量所激發(fā)的電子數(shù),N為探測器像元數(shù),對于Aeolus來說,N= 16,Nnoise為探測器暗電流激發(fā)的電子數(shù)。對于Rayleigh通道,NFP(z,1)和NFP(z,2)分別為F-P標準具1 和2 所接收到的Rayleigh 光譜能量激發(fā)的電子數(shù),和分別為F-P 標準具1和2接收到的總信號的能量激發(fā)的電子數(shù)。

    大氣不均勻性會也會對測風不確定度產(chǎn)生重要影響(Baars 等,2019;Sun 等,2014)。星載激光測風雷達探測到的風速的垂直分辨率較低,介于250 m—2 km。在數(shù)據(jù)處理時,默認探測到的風速數(shù)據(jù)為該高度層幾何中心位置的風速,而實際上探測的風速為該高度層的光學重心位置的風速。當大氣垂直均勻時,幾何中心位置與光學重心位置重合,此時探測到風速的位置是準確的。但是由于大氣的垂直不均勻性,該高度層的幾何中心與光學重心往往不重合。而且由于垂直風切變是普遍存在的,從而導致了測風誤差的產(chǎn)生,統(tǒng)計表明,大氣垂直風切變的平均值約為0.01 s-1。Sun等(2014)的研究表明,對于Rayleigh 通道而言,垂直不均勻性導致的Rayleigh 通道測風誤差約為0.5 m·s-1,對于Mie 通道而言,垂直不均勻性導致的測風誤差約為0.4 m·s-1。一般來說,垂直不均勻性導致的測風誤差在Rayleigh 通道可進行修正,但是在Mie通道難以修正。

    基于Aeolus、HDWL和CDWL的系統(tǒng)結構,并利用大氣光學特性數(shù)據(jù)庫,考慮泊松噪聲和大氣不均勻性對測風精度的影響,歐洲空間局、美國和日本分別模擬了不同高度處不同質量的風速探測結果所占百分比情況,其中歐洲空間局Aeolus體制測風質量模擬仿真結果如圖7所示(Stoffelen等,2005),美國HDWL 體制測風質量的模擬結果如圖8所示(Masutani等,2010),日本CDWL測風質量模擬結果如圖9所示(Baron等,2017)。

    圖7 Aeolus測風質量的模擬仿真結果(Stoffelen等,2005)Fig.7 The simulated wind observation accuracy of Aeolus(Stoffelen et al.,2005)

    圖8 HDWL測風質量的模擬仿真結果(Masutani等,2010)Fig.8 The simulated wind observation accuracy of HDWL(Masutani et al.,2010)

    圖9 CDWL測風質量模擬結果(Baron等,2017)Fig.9 The simulated wind observation accuracy of CDWL(Baron et al.,2017)

    由圖7可得,在歐洲空間局Aeolus測風體制下大部分的風速測量值的精度介于1.4—2.8 m·s-1,其高度范圍多位于對流層的中高層,主要是Rayleigh通道的探測結果。質量較好的風速探測結果多集中在邊界層和對流層底部,主要為Mie通道的探測結果。質量較差的風速測量結果也多集中在在邊界層和對流層底部。這主要是由于光學厚度較厚的云的存在,導致激光穿透云層后能量較弱,風場探測結果具有較大的不確定度。對比圖7(a)(b)可知,當考慮風切變時,質量較差的風速探測值所占比例明顯增加,說明風切變對星載激光雷達的測風質量具有較大影響,這部分測風不確定度主要是由于大氣不均勻性導致的。

    由圖8(a)可得,美國HDWL技術體制下,相干探測測得的風場數(shù)據(jù)約占觀測到總數(shù)據(jù)的30%,測風精度較高,測風不確定度多介于0—2 m·s-1,少部分測風不確定度介于2—3 m·s-1。其中15 km 以上高度層的探測示蹤物主要是氣溶膠;3—15 km高度層探測示蹤物多為卷云和不透明云;0—3 km高度層的探測示蹤物主要是氣溶膠。由圖8(b)可得,直接探測測得的風場數(shù)據(jù)約占總數(shù)據(jù)的70%,其中14 km 以上高度層的測風不確定度多介于2—3 m·s-1,3—14 km 高度層測風不確定度多介于1—2 m·s-1,0—3 km 測風不確定度較大,且不可用數(shù)據(jù)所占比例較大,這主要是由于激光穿透能力有限所致。結合圖8(a)(b)可得,相干探測測風精度較高,但探測范圍小,直接探測測風精度低,但覆蓋范圍廣,二者結合可以較好地實現(xiàn)全球風場的探測。

    由圖9可得,衛(wèi)星升軌軌道和降軌軌道對應的平均風速大小相近,但是方向相反,存在著互為鏡像的關系;在日本CDWL 測風體制下,升軌軌道和降軌軌道測風不確定度的全球分布也比較相似;對于相干測風來說,由于低層氣溶膠數(shù)密度相對較大,測風精度較高,在中高層大氣中,云滴粒子或氣溶膠粒子數(shù)密度相對較小,測風精度較低,該仿真結果與圖8(a)具有一定的相似之處。Baron等(2017)和Okamoto等(2018)的研究表明,對于星載相干激光測風雷達而言,測風不確定度主要與大氣不均勻性有關。

    3.3 基于對數(shù)值預報系統(tǒng)改進作用的評估

    OSSE(Observing System Simulation Experiments)系統(tǒng)是評估某個新引入的觀測量能否有助于改善數(shù)值預報系統(tǒng)精度的有力工具。Rohaly和Krishnamurti(1993)最早使用了OSSE 系統(tǒng)評估了全球風場探測數(shù)據(jù)是否能夠改善數(shù)值預報系統(tǒng)的預報結果。OSSE 運行時主要分兩步進行,首先利用現(xiàn)有的數(shù)值預報系統(tǒng)進行預報,然后將新的觀測量同化到預報系統(tǒng)的觀測場中,再次進行預報,最后將兩次預報結果進行對比,即可評估該觀測量對數(shù)值預報系統(tǒng)的改善作用(Lahoz 等,2005)。在運行OSSE 之前,為了同化新的觀測量,需要對數(shù)據(jù)同化系統(tǒng)進行重構。另外在同化數(shù)據(jù)時,觀測量的不確定度是一個不可或缺的量,還需建立新觀測量的模擬仿真系統(tǒng)以計算其不確定度。所以OSSE實驗需要消耗較大的人力和計算力,但相比建立新的觀測系統(tǒng),這些消耗只是其中一小部分。

    Stoffelen 等(2006)利用OSSE 系統(tǒng)對星載激光測風雷達能否改善數(shù)值預報產(chǎn)品精度進行了評估。在評估時分別使用了兩種資料作為觀測場,一是無星載激光雷達風場的模式預報資料,二是使用了模式預報資料加模擬的星載激光雷達觀測風場資料(Marseille 和Stoffelen,2003;Reitebuch等,2017)。將上述兩種資料分別利用4D-Var 技術同化到觀測場中進行預報。結果表明,增加了星載激光雷達風場觀測資料后,預報模式的結果得到了明顯的改善,但是當星載激光雷達測風結果精度較低時,對預報結果產(chǎn)生了副作用。Tan 等(2007)的研究也驗證了這一點,而且其結果表明,現(xiàn)在的數(shù)值同化系統(tǒng)可從Aeolus 探測的徑向風場中提取有用信息,能夠減小預報的誤差和再分析風場的不確定度。

    Masutani 等(2010)使用了OSSE 評估了HDWL獲取的風場資料對預報系統(tǒng)的改善作用,在實驗過程中使用的風場、溫度、氣壓和濕度等數(shù)據(jù)取自NCEP 的再分析資料,其中模擬風場資料的系統(tǒng)誤差和不確定度是由DLSM(Defense Lidar Simulation Model)基于泊松噪聲的原理產(chǎn)生的。實驗結果表明,HDWL 對預報系統(tǒng)改善作用的置信度達到95%以上,尤其是在熱帶地區(qū)的高層大氣下,改善作用更加明顯。

    Okamoto 等(2018)評估了CDWL 對預報系統(tǒng)的改善作用,在將風場資料加入到觀測場之前,首先利用仿真系統(tǒng)基于氣溶膠和云場的觀測資料模擬了風場觀測量的信噪比、不確定度以及誤差,并且評估了極地軌道和熱帶軌道對預報系統(tǒng)的改善作用,結果表明兩種軌道都能較好地改善天氣預報的結果。在熱帶地區(qū),加入星載激光測風雷達全球風場觀測資料后,預報相對誤差能夠減小2%左右。另外該研究表明,低估CDWL 對預報系統(tǒng)改善作用的很大一部分原因是因為不能較好地模擬觀測風場的不確定度。

    美國和日本的星載激光雷達技術體制通過在一顆衛(wèi)星上安裝多個望遠鏡的方法實現(xiàn)水平風場矢量的探測,這種技術體制也可通過多顆衛(wèi)星聯(lián)合探測來實現(xiàn)。Marseille 等(2008)評估了5 種聯(lián)合探測體制對預報系統(tǒng)的改善作用,這5種聯(lián)合探測體制如圖10所示,分別為:圖10(a)Aeolus 單顆衛(wèi)星探測;圖10(b)Dual-perspective:其軌道與Aeolus 軌道一致,前后兩顆星載激光測風雷達的探測方位角相差90°,可實現(xiàn)水平二維風場的探測;圖10(c)Tandem-Aeolus:一顆衛(wèi)星運行在Aeolus的軌道,另一顆衛(wèi)星軌道的升交點與Aeolus軌道升交點相差50 min 左右,其他參數(shù)與Aeolus軌道一致,兩顆星載激光測風雷達的技術體制與Aeolus相同,探測時,兩顆衛(wèi)星基本上位于相同緯度面上;圖10 (d)Triple-Aeolus:與Tandem-Aeolus 類似,3 顆采用Aeolus 技術體制的星載激光測風雷達進行聯(lián)合探測,其軌道升交點相差30 min左右;圖10(e)Dual-inclination:兩顆衛(wèi)星運行于不同傾角的軌道,其傾角分別為97.2°和70°,第二顆衛(wèi)星主要用于探測中低緯風暴多發(fā)地區(qū)的大氣風場。

    圖10 5種星載激光雷達聯(lián)合測風技術方案探測數(shù)據(jù)覆蓋情況示意圖(Marseille等,2008)Fig.10 The data coverage of five observation scenarios of spaceborne DWLs(Marseille et al.,2008)

    Marseille等(2008)評估結果表明:(1)聯(lián)合探測技術體制對數(shù)值預報系統(tǒng)的改善作用要明顯優(yōu)于單顆星載激光測風雷達;(2)提高探測風場的覆蓋面更有助于提高預報系統(tǒng)的精度,如Tandem-Aeolus 體制、 Dual-perspective 體制和Dual-inclination 體制探測到的徑向風速數(shù)量都是Aeolus 的兩倍,但從對預報系統(tǒng)的改善作用來說,Tandem-Aeolus 體制的改善效果在三者中最好,后兩種技術體制在特定情況下對數(shù)值預報系統(tǒng)的改善效果好,如Dual-inclination 體制的探測結果有助于改善中低緯地區(qū),尤其是大西洋地區(qū)風暴的預測結果,在Dual-perspective 體制下,衛(wèi)星星下點附近的數(shù)值預報結果改善作用明顯,對比其他4 種探測體制,該體制其改善效果最佳,但是由于探測覆蓋面較低,從全球整體來看,其對數(shù)值預報系統(tǒng)的改善作用比Tandem-Aeolus 體制低;(3)Triple-Aeolus 體制相比Tandem-Aeolus 體制,對數(shù)值預報系統(tǒng)有一定的改善作用,Tandem-Aeolus 體制相比Aeolus,其改善作用大約增加了66%,增加第3 顆衛(wèi)星之后,相比兩顆衛(wèi)星,其增加的改善作用約為20%,該結果表明,3 顆以上星載激光測風雷達對預報系統(tǒng)的改善作用有限,太多的衛(wèi)星可能帶來數(shù)據(jù)冗余。

    OSSE 結果表明,3 種星載激光雷達測風技術體制都能有效地改善數(shù)值預報系統(tǒng)的預報結果。水平二維風場的探測與徑向風場的探測相比,改善星下點附近數(shù)值預報結果的作用更加明顯。從對全球預報的整體改善效果來看,更廣泛的探測覆蓋面相比水平二維風場的探測,其效果更好。另外在運行OSSE 實驗時,需要通過星載激光測風雷達仿真系統(tǒng)對觀測風場的不確定度進行模擬,否則可能無法準確評估其對數(shù)值預報系統(tǒng)的改善效果。

    通過Aeolus、HDWL 和CDWL 這3 種星載激光測風雷達技術體制的數(shù)據(jù)獲取率和獲取量、風場探測精度以及對數(shù)值預報系統(tǒng)的改進作用方面的對比可得,美國HDWL 技術體制可以探測平行于星下點軌跡上兩條路徑的風場數(shù)據(jù),而Aeolus 和CDWL 只能探測一條路徑的風場數(shù)據(jù),故HDWL具有更廣的數(shù)據(jù)覆蓋面,其測風精度也相對較高,而且HDWL 可以實現(xiàn)水平二維風場的探測,日本CDWL雖然能實現(xiàn)水平二維風場探測,但其數(shù)據(jù)獲取率較低,結合Marseille 等(2008)對5 種多星聯(lián)合探測技術體制的評估結論可以預估,HDWL 技術體制對數(shù)值預報系統(tǒng)的改善作用相比其他兩種體制較明顯。但上述結論主要是通過前人所做實驗的結論進行推理的,目前缺乏將3種技術體制置于完全相同的大氣環(huán)境中進行模擬評估的實驗,而且美國和日本的星載激光測風雷達目前仍處于論證階段,其參數(shù)還在調(diào)整中,參數(shù)的調(diào)整會對測風精度產(chǎn)生較大的影響,并且雷達在太空運行后,由于環(huán)境的改變,其發(fā)揮的性能還有待評估。所以這兩種技術體制的激光雷達測風精度以及對預報系統(tǒng)的改善效果還有待進一步研究。

    4 結論

    2018年全球首顆星載激光測風雷達衛(wèi)星Aeolus的成功發(fā)射標志著人類在全球風場的探測方面邁入了新的階段。目前,美國和日本也在積極論證和研制HDWL 和CDWL 星載激光測風雷達技術體制,中國也在開展相關研究工作。不同星載激光測風雷達技術體制決定了研制成本和發(fā)揮的作用。目前的評估結果表明星載激光測風雷達對數(shù)值預報系統(tǒng)具有較好地改善作用,使用單顆衛(wèi)星進行全球風場探測時,美國HDWL 技術體制相比歐洲空間局Aeolus 和日本CDWL 技術體制能夠獲取更多的風場探測數(shù)據(jù),其測風精度較高,將更有助于數(shù)值預報系統(tǒng)結果的改善。使用多顆衛(wèi)星進行聯(lián)合探測時,更廣泛的探測覆蓋面相比水平二維風場的探測,更有助于改善全球數(shù)值預報系統(tǒng)的預報效果。水平二維風場的探測,對改善衛(wèi)星軌道附近的數(shù)值預報系統(tǒng)結果效果最佳。

    星載激光測風資料對全球數(shù)值天氣預報效果改善作用的評估是一項復雜的工作,其中牽涉到對測風不確定度的評估,而星載激光測風資料的不確定度不僅與激光雷達技術參數(shù)有關,還與衛(wèi)星軌道參數(shù)和大氣不均勻性有關。因此,開展星載激光測風技術體制的論證還需要做深入的研究,以便獲得最優(yōu)的技術方案。

    猜你喜歡
    風場體制大氣
    大氣的呵護
    軍事文摘(2023年10期)2023-06-09 09:15:06
    基于FLUENT的下?lián)舯┝魅S風場建模
    試論烏俄案對多邊貿(mào)易體制的維護
    “最美風場”的贏利法則
    能源(2017年8期)2017-10-18 00:47:39
    建立“大健康”體制是當務之急
    為“三醫(yī)聯(lián)動”提供體制保障
    大氣古樸揮灑自如
    大氣、水之后,土十條來了
    側向風場中無人機的飛行研究
    建立高效的政府辦醫(yī)體制
    午夜免费激情av| 久久精品国产亚洲网站| 91久久精品国产一区二区三区| 97超视频在线观看视频| 国产午夜精品久久久久久一区二区三区 | 九九热线精品视视频播放| av免费在线看不卡| 日韩欧美精品v在线| 日韩人妻高清精品专区| 国产老妇女一区| 天天一区二区日本电影三级| 国产精品人妻久久久久久| 国产久久久一区二区三区| 欧美色欧美亚洲另类二区| 美女cb高潮喷水在线观看| 久久综合国产亚洲精品| 国产精品一区二区三区四区久久| 91午夜精品亚洲一区二区三区| 日本-黄色视频高清免费观看| 久久鲁丝午夜福利片| 亚洲在线自拍视频| 国产精品久久久久久久久免| 日本一二三区视频观看| 九色成人免费人妻av| 搡女人真爽免费视频火全软件 | 免费高清视频大片| 一级a爱片免费观看的视频| 免费看a级黄色片| 午夜福利高清视频| 亚洲av中文av极速乱| 男女啪啪激烈高潮av片| 亚洲专区国产一区二区| 久久久久国产精品人妻aⅴ院| 99热精品在线国产| 久久精品国产鲁丝片午夜精品| aaaaa片日本免费| 1024手机看黄色片| 91久久精品国产一区二区成人| 亚洲精品国产成人久久av| 国产精品国产三级国产av玫瑰| 国产精品av视频在线免费观看| 国产三级在线视频| 国产高潮美女av| 熟女电影av网| 黄色一级大片看看| 中国美白少妇内射xxxbb| 日韩大尺度精品在线看网址| 久久久久久久久久久丰满| 成年免费大片在线观看| 国产欧美日韩一区二区精品| 禁无遮挡网站| 亚洲高清免费不卡视频| 免费看美女性在线毛片视频| 色综合站精品国产| 在线免费观看的www视频| 97人妻精品一区二区三区麻豆| 国产三级中文精品| 亚洲人成网站在线播放欧美日韩| 99久久精品一区二区三区| 欧美国产日韩亚洲一区| 热99re8久久精品国产| 久久久国产成人免费| 18禁裸乳无遮挡免费网站照片| 嫩草影院新地址| 淫秽高清视频在线观看| 五月玫瑰六月丁香| 免费电影在线观看免费观看| 搡女人真爽免费视频火全软件 | 久久精品综合一区二区三区| 免费高清视频大片| 免费不卡的大黄色大毛片视频在线观看 | 中文字幕av在线有码专区| 成人综合一区亚洲| 天堂影院成人在线观看| 国产精品久久久久久久久免| 亚洲成a人片在线一区二区| 国产成人aa在线观看| 亚洲最大成人手机在线| 一边摸一边抽搐一进一小说| 欧美zozozo另类| 美女xxoo啪啪120秒动态图| 久久久久久国产a免费观看| 亚洲自偷自拍三级| 精品一区二区免费观看| 露出奶头的视频| 日本一二三区视频观看| 国产亚洲av嫩草精品影院| 九九爱精品视频在线观看| 午夜日韩欧美国产| 女同久久另类99精品国产91| 欧美日韩在线观看h| 国产av不卡久久| 最近手机中文字幕大全| 一级黄片播放器| 国产av一区在线观看免费| 看片在线看免费视频| 国产精品电影一区二区三区| 可以在线观看毛片的网站| 亚洲,欧美,日韩| 一进一出抽搐gif免费好疼| 草草在线视频免费看| 久久久久久久午夜电影| 六月丁香七月| 最近2019中文字幕mv第一页| 国产 一区 欧美 日韩| 亚洲无线观看免费| 色哟哟哟哟哟哟| 国产一级毛片七仙女欲春2| 久久人人精品亚洲av| 午夜a级毛片| 亚洲18禁久久av| 99九九线精品视频在线观看视频| 高清毛片免费观看视频网站| aaaaa片日本免费| 国产综合懂色| 91精品国产九色| 亚洲av不卡在线观看| 亚洲精品456在线播放app| 女的被弄到高潮叫床怎么办| 特大巨黑吊av在线直播| 欧美成人精品欧美一级黄| 成人亚洲精品av一区二区| 99久久成人亚洲精品观看| 亚洲三级黄色毛片| 亚洲av成人精品一区久久| 精品一区二区三区av网在线观看| 国内精品一区二区在线观看| 欧美日韩一区二区视频在线观看视频在线 | 综合色丁香网| 一进一出好大好爽视频| 国产一区二区三区av在线 | 久久久久久九九精品二区国产| 亚洲激情五月婷婷啪啪| 自拍偷自拍亚洲精品老妇| 久久人人精品亚洲av| 国产三级中文精品| 色哟哟哟哟哟哟| 精品熟女少妇av免费看| 九九在线视频观看精品| 免费观看在线日韩| 国产白丝娇喘喷水9色精品| 日本熟妇午夜| 熟妇人妻久久中文字幕3abv| 精品久久久久久久人妻蜜臀av| 欧美国产日韩亚洲一区| 国产91av在线免费观看| 桃色一区二区三区在线观看| 麻豆精品久久久久久蜜桃| 在现免费观看毛片| 久久九九热精品免费| 毛片女人毛片| 免费观看的影片在线观看| 国产精品久久久久久av不卡| 一区福利在线观看| 午夜福利18| av在线观看视频网站免费| 日韩欧美精品v在线| www日本黄色视频网| 日韩精品有码人妻一区| 深夜a级毛片| 久久国产乱子免费精品| 欧美又色又爽又黄视频| 日韩精品青青久久久久久| 一进一出抽搐gif免费好疼| 国产高清视频在线播放一区| 精品一区二区三区人妻视频| 日本成人三级电影网站| 网址你懂的国产日韩在线| 日韩精品中文字幕看吧| 黄色配什么色好看| 人妻制服诱惑在线中文字幕| 色5月婷婷丁香| 日日摸夜夜添夜夜添av毛片| 国产精品亚洲美女久久久| 欧美日韩国产亚洲二区| 97超碰精品成人国产| 日韩欧美在线乱码| 久久久欧美国产精品| 国产精品人妻久久久久久| 午夜日韩欧美国产| 国产精品亚洲美女久久久| 国产伦在线观看视频一区| 日韩在线高清观看一区二区三区| 欧美日本视频| 国产成人一区二区在线| 两性午夜刺激爽爽歪歪视频在线观看| 特级一级黄色大片| 一个人看视频在线观看www免费| 伊人久久精品亚洲午夜| 亚洲国产精品合色在线| 三级经典国产精品| 亚洲专区国产一区二区| av福利片在线观看| 日本黄色视频三级网站网址| 成年女人看的毛片在线观看| 精品一区二区三区视频在线观看免费| 日本-黄色视频高清免费观看| 成年女人看的毛片在线观看| 悠悠久久av| 男人舔奶头视频| 亚洲乱码一区二区免费版| 日本与韩国留学比较| 亚洲精品成人久久久久久| 久久精品国产自在天天线| 菩萨蛮人人尽说江南好唐韦庄 | 99热网站在线观看| 亚洲av免费在线观看| 五月玫瑰六月丁香| 欧美一区二区国产精品久久精品| 亚洲无线在线观看| 黄色视频,在线免费观看| 美女xxoo啪啪120秒动态图| 在线观看美女被高潮喷水网站| 一个人免费在线观看电影| 69av精品久久久久久| 六月丁香七月| 久久九九热精品免费| 可以在线观看的亚洲视频| 国产69精品久久久久777片| 中出人妻视频一区二区| 蜜桃久久精品国产亚洲av| 精品久久久久久久人妻蜜臀av| 久久久久国产网址| 久久鲁丝午夜福利片| 人人妻人人澡人人爽人人夜夜 | 麻豆精品久久久久久蜜桃| 欧美色欧美亚洲另类二区| 91精品国产九色| 在线播放国产精品三级| 中文字幕久久专区| 高清午夜精品一区二区三区 | 嫩草影视91久久| avwww免费| 搡老妇女老女人老熟妇| 亚洲最大成人手机在线| 国产精品久久久久久久久免| 日韩在线高清观看一区二区三区| 高清午夜精品一区二区三区 | 久久精品91蜜桃| 亚洲美女黄片视频| 深夜精品福利| 一级毛片久久久久久久久女| 亚洲国产色片| 99热6这里只有精品| 亚洲精品乱码久久久v下载方式| 精品久久久噜噜| 又黄又爽又免费观看的视频| 又粗又爽又猛毛片免费看| 久久亚洲精品不卡| 露出奶头的视频| 久久天躁狠狠躁夜夜2o2o| 国产精品1区2区在线观看.| 麻豆久久精品国产亚洲av| 精品久久久久久久久亚洲| 我的女老师完整版在线观看| 超碰av人人做人人爽久久| 毛片女人毛片| 极品教师在线视频| 大又大粗又爽又黄少妇毛片口| 久久久久免费精品人妻一区二区| 女同久久另类99精品国产91| 久久久午夜欧美精品| 国产精品免费一区二区三区在线| or卡值多少钱| 97超级碰碰碰精品色视频在线观看| 赤兔流量卡办理| 国产白丝娇喘喷水9色精品| 97碰自拍视频| 美女大奶头视频| 亚洲精品色激情综合| 亚洲av一区综合| 国产片特级美女逼逼视频| 99久久成人亚洲精品观看| 91在线观看av| 国产熟女欧美一区二区| 日韩精品有码人妻一区| 国产精品爽爽va在线观看网站| 给我免费播放毛片高清在线观看| 亚洲一区高清亚洲精品| 少妇的逼水好多| 亚洲精品一区av在线观看| 91久久精品国产一区二区三区| 别揉我奶头~嗯~啊~动态视频| 成人亚洲精品av一区二区| 小蜜桃在线观看免费完整版高清| 亚洲国产精品合色在线| 国产女主播在线喷水免费视频网站 | 18禁在线无遮挡免费观看视频 | 老熟妇乱子伦视频在线观看| 三级国产精品欧美在线观看| 欧美bdsm另类| 国产精品福利在线免费观看| 久久99热这里只有精品18| 亚洲熟妇熟女久久| 狠狠狠狠99中文字幕| 精品人妻熟女av久视频| 尤物成人国产欧美一区二区三区| 久久亚洲国产成人精品v| 欧美日韩一区二区视频在线观看视频在线 | 日本熟妇午夜| 欧美一区二区精品小视频在线| www日本黄色视频网| 成年女人永久免费观看视频| 亚洲五月天丁香| 麻豆精品久久久久久蜜桃| 老师上课跳d突然被开到最大视频| 久久久久久久久中文| 变态另类成人亚洲欧美熟女| 国产真实伦视频高清在线观看| 床上黄色一级片| 中国美女看黄片| 国产单亲对白刺激| 免费看光身美女| eeuss影院久久| 两个人视频免费观看高清| 久久久久精品国产欧美久久久| 熟女电影av网| eeuss影院久久| 日本精品一区二区三区蜜桃| 天美传媒精品一区二区| 欧美色视频一区免费| 一卡2卡三卡四卡精品乱码亚洲| 不卡视频在线观看欧美| 12—13女人毛片做爰片一| 久久国产乱子免费精品| videossex国产| 亚洲国产精品久久男人天堂| 国语自产精品视频在线第100页| av免费在线看不卡| 免费人成视频x8x8入口观看| 亚洲精品乱码久久久v下载方式| 两性午夜刺激爽爽歪歪视频在线观看| 久久人人精品亚洲av| 偷拍熟女少妇极品色| 韩国av在线不卡| 18+在线观看网站| 亚洲精品亚洲一区二区| 欧美一区二区国产精品久久精品| 国产成人aa在线观看| 午夜视频国产福利| 国产女主播在线喷水免费视频网站 | 哪里可以看免费的av片| 少妇的逼好多水| 男女边吃奶边做爰视频| 久久精品综合一区二区三区| or卡值多少钱| 国产精品日韩av在线免费观看| 天堂√8在线中文| 听说在线观看完整版免费高清| 一夜夜www| 性色avwww在线观看| 欧美色视频一区免费| 国产中年淑女户外野战色| 国产精品久久久久久久久免| 国内久久婷婷六月综合欲色啪| 精品一区二区三区av网在线观看| 波多野结衣巨乳人妻| 蜜臀久久99精品久久宅男| av天堂中文字幕网| 亚洲精品日韩av片在线观看| 日韩欧美一区二区三区在线观看| 成年女人永久免费观看视频| 国产精品精品国产色婷婷| 久久午夜福利片| 老司机福利观看| 亚洲成人av在线免费| 国产高清激情床上av| 色播亚洲综合网| 亚洲久久久久久中文字幕| 亚洲av五月六月丁香网| 精品人妻视频免费看| 我要看日韩黄色一级片| 99热全是精品| 亚洲成av人片在线播放无| 亚洲精品一区av在线观看| 最近视频中文字幕2019在线8| 人妻少妇偷人精品九色| 看片在线看免费视频| 国产亚洲精品久久久久久毛片| 国产蜜桃级精品一区二区三区| 99久久无色码亚洲精品果冻| 99热这里只有精品一区| 成年版毛片免费区| 日日啪夜夜撸| 国产色爽女视频免费观看| 欧美色欧美亚洲另类二区| 搡老岳熟女国产| 我的女老师完整版在线观看| 亚洲激情五月婷婷啪啪| 寂寞人妻少妇视频99o| 成人av一区二区三区在线看| 国产av一区在线观看免费| 亚洲人成网站高清观看| 国产成人freesex在线 | 国内揄拍国产精品人妻在线| 黑人高潮一二区| 精品免费久久久久久久清纯| 成人无遮挡网站| 成人鲁丝片一二三区免费| 久久精品国产亚洲av天美| 一个人看的www免费观看视频| 老熟妇仑乱视频hdxx| 国产精品无大码| 午夜影院日韩av| 淫秽高清视频在线观看| 夜夜夜夜夜久久久久| 国产精品乱码一区二三区的特点| 午夜福利在线在线| 亚洲av熟女| 精品久久久久久成人av| 亚洲第一区二区三区不卡| 亚洲欧美成人综合另类久久久 | 国产精品一区二区免费欧美| 亚洲中文字幕日韩| 有码 亚洲区| 女人被狂操c到高潮| 国产成人freesex在线 | 国产黄色视频一区二区在线观看 | 色哟哟·www| 色播亚洲综合网| 特大巨黑吊av在线直播| 亚洲三级黄色毛片| 身体一侧抽搐| 久久人人爽人人爽人人片va| 欧美在线一区亚洲| 中文字幕人妻熟人妻熟丝袜美| 免费看日本二区| 看片在线看免费视频| 大型黄色视频在线免费观看| 久久精品影院6| 99久久精品国产国产毛片| 国产aⅴ精品一区二区三区波| h日本视频在线播放| 91av网一区二区| 亚洲av一区综合| 极品教师在线视频| 亚洲中文日韩欧美视频| 夜夜看夜夜爽夜夜摸| 亚洲自拍偷在线| 少妇猛男粗大的猛烈进出视频 | 国产 一区 欧美 日韩| 真实男女啪啪啪动态图| 国产高清不卡午夜福利| 国产精品无大码| 哪里可以看免费的av片| 我要看日韩黄色一级片| 别揉我奶头~嗯~啊~动态视频| 人妻丰满熟妇av一区二区三区| 成年女人看的毛片在线观看| 波多野结衣高清作品| 日本免费a在线| 亚洲国产精品成人综合色| 亚洲成a人片在线一区二区| 久久精品国产亚洲av涩爱 | 日韩精品有码人妻一区| 一本一本综合久久| 亚洲,欧美,日韩| 国产精品av视频在线免费观看| 欧美日韩乱码在线| 免费人成在线观看视频色| 精品久久久久久久久亚洲| 最近在线观看免费完整版| 老熟妇乱子伦视频在线观看| 免费不卡的大黄色大毛片视频在线观看 | 免费电影在线观看免费观看| 成年av动漫网址| 黄片wwwwww| 特大巨黑吊av在线直播| 亚洲国产精品sss在线观看| 99久久精品国产国产毛片| 久久鲁丝午夜福利片| 欧美一区二区国产精品久久精品| 男女那种视频在线观看| 无遮挡黄片免费观看| 99热6这里只有精品| 永久网站在线| 亚洲aⅴ乱码一区二区在线播放| 美女被艹到高潮喷水动态| 亚洲电影在线观看av| 少妇裸体淫交视频免费看高清| av黄色大香蕉| 日韩欧美精品v在线| 性色avwww在线观看| 亚洲欧美日韩东京热| 在线观看美女被高潮喷水网站| 日日摸夜夜添夜夜添小说| 亚州av有码| 国产乱人偷精品视频| 久久精品久久久久久噜噜老黄 | 国产高清激情床上av| 中文字幕免费在线视频6| 少妇人妻一区二区三区视频| 久久久久性生活片| 国产三级中文精品| 麻豆成人午夜福利视频| 寂寞人妻少妇视频99o| 国内精品宾馆在线| 午夜福利在线观看免费完整高清在 | 亚州av有码| 日本一本二区三区精品| 精品乱码久久久久久99久播| 老师上课跳d突然被开到最大视频| www.色视频.com| 在线播放国产精品三级| 精华霜和精华液先用哪个| 乱人视频在线观看| 在线免费十八禁| 亚洲第一区二区三区不卡| 岛国在线免费视频观看| 日韩欧美在线乱码| 国产国拍精品亚洲av在线观看| 国产精品永久免费网站| 精品一区二区免费观看| av在线老鸭窝| 欧美性猛交黑人性爽| 午夜日韩欧美国产| 亚洲精品国产成人久久av| 精品久久久久久久久久免费视频| 亚洲综合色惰| 精品99又大又爽又粗少妇毛片| 国产乱人视频| 久99久视频精品免费| 日韩欧美 国产精品| 国产欧美日韩一区二区精品| 不卡一级毛片| 久久久色成人| 身体一侧抽搐| 成人性生交大片免费视频hd| 日韩欧美 国产精品| 51国产日韩欧美| 亚洲真实伦在线观看| 成人二区视频| 一a级毛片在线观看| 亚洲久久久久久中文字幕| 成人国产麻豆网| 日本在线视频免费播放| 菩萨蛮人人尽说江南好唐韦庄 | 女同久久另类99精品国产91| 乱系列少妇在线播放| 91av网一区二区| 国内精品久久久久精免费| 在线观看美女被高潮喷水网站| 精品一区二区三区视频在线观看免费| 最新在线观看一区二区三区| 丰满人妻一区二区三区视频av| 日韩在线高清观看一区二区三区| 美女内射精品一级片tv| 久久精品人妻少妇| 中文字幕免费在线视频6| 国产一区二区亚洲精品在线观看| 色综合色国产| 亚洲成人久久爱视频| 国产视频内射| 看非洲黑人一级黄片| 国产在线精品亚洲第一网站| 一级毛片久久久久久久久女| 全区人妻精品视频| 老司机午夜福利在线观看视频| 亚洲自偷自拍三级| 久久久久久伊人网av| 晚上一个人看的免费电影| 免费高清视频大片| 国产探花极品一区二区| 国内精品宾馆在线| 精品久久久久久久久亚洲| 国产精品电影一区二区三区| 日日摸夜夜添夜夜爱| 久久久久久九九精品二区国产| 亚洲欧美精品综合久久99| 亚洲国产欧洲综合997久久,| 日韩精品中文字幕看吧| 亚洲国产精品成人久久小说 | 丰满的人妻完整版| 日韩在线高清观看一区二区三区| 春色校园在线视频观看| 精品久久久久久久久亚洲| 真人做人爱边吃奶动态| 一区二区三区免费毛片| 国产午夜精品久久久久久一区二区三区 | 看免费成人av毛片| 亚洲国产色片| 久久亚洲精品不卡| 国产成人a区在线观看| 亚洲内射少妇av| 国产成人一区二区在线| 给我免费播放毛片高清在线观看| 欧美日韩精品成人综合77777| 国产大屁股一区二区在线视频| 国产一区二区亚洲精品在线观看| 日韩国内少妇激情av| 最近2019中文字幕mv第一页| 国产蜜桃级精品一区二区三区| 97超碰精品成人国产| 日本欧美国产在线视频| 久久亚洲精品不卡| 亚洲性久久影院| 国产精品嫩草影院av在线观看| 久久亚洲精品不卡| 国产精品亚洲美女久久久| 麻豆成人午夜福利视频| 国产精品久久久久久精品电影| 日日摸夜夜添夜夜添av毛片| 免费观看的影片在线观看| av国产免费在线观看| 日韩欧美在线乱码| 日韩亚洲欧美综合| 国内揄拍国产精品人妻在线| 中国美白少妇内射xxxbb| 精品一区二区三区视频在线| 亚洲欧美成人综合另类久久久 | 亚洲人成网站在线观看播放| 国产亚洲精品久久久com| 日本 av在线| 欧美激情国产日韩精品一区| 国产真实乱freesex|