999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

激光在春季江淮氣旋降水中的衰減特性研究

2017-03-14 08:00:28陳凱諾邵利民魏海亮戴仁威
海洋技術(shù)學(xué)報(bào) 2017年1期

陳凱諾,邵利民,魏海亮,戴仁威,艾 陽(yáng),穆 楊

(1.海軍航空兵學(xué)院,遼寧 葫蘆島 125100;2.海軍大連艦艇學(xué)院 軍事海洋系,遼寧 大連 116018;3.中國(guó)人民解放軍92956部隊(duì) 旅順裝備技術(shù)保障大隊(duì)三中隊(duì) ,遼寧 大連 116041;4.海軍舟山水警區(qū)司令部氣象臺(tái),浙江 舟山 316000;5.東海艦隊(duì)水文氣象中心 參謀部37分隊(duì),浙江 寧波 315000)

激光在春季江淮氣旋降水中的衰減特性研究

陳凱諾1,2,邵利民2,魏海亮3,戴仁威4,艾 陽(yáng)5,穆 楊1

(1.海軍航空兵學(xué)院,遼寧 葫蘆島 125100;2.海軍大連艦艇學(xué)院 軍事海洋系,遼寧 大連 116018;3.中國(guó)人民解放軍92956部隊(duì) 旅順裝備技術(shù)保障大隊(duì)三中隊(duì) ,遼寧 大連 116041;4.海軍舟山水警區(qū)司令部氣象臺(tái),浙江 舟山 316000;5.東海艦隊(duì)水文氣象中心 參謀部37分隊(duì),浙江 寧波 315000)

在雨中艦載激光武器和激光雷達(dá)的工作效能受到制約。因此,研究激光在降水中的衰減特性對(duì)軍事行動(dòng)具有重要的科學(xué)意義。基于WRF中尺度氣象研究模式,對(duì)2015年3月30日-4月2日的降水過程進(jìn)行模擬分析,江淮氣旋由于倒槽鋒生產(chǎn)生,江蘇地區(qū)和朝鮮半島、日本等的降水過程就是高空西風(fēng)急流與低空急流耦合的結(jié)果,與此同時(shí)還有700 hPa切變線的配合。在氣旋降水條件下,傳輸性能較好的10.6μm遠(yuǎn)紅外波仍然受到衰減,在小雨的情況下 (降雨率為0.25 mm/h)探測(cè)距離可損失10%;在中雨的情況下(降雨率為2.7mm/h)探測(cè)距離為正常情況下的70%;而在大雨情況下(降雨率為6mm/h)探測(cè)距離僅為正常情況下的50%。

江淮氣旋;低空急流;倒槽鋒生;激光衰減系數(shù);探測(cè)距離

隨著激光技術(shù)的不斷發(fā)展,激光的應(yīng)用也越來越廣泛,特別是在軍事通信領(lǐng)域有著廣闊的發(fā)展空間。各種特殊天氣的影響是制約無線光通信發(fā)展的重要因素,雨霧等天氣對(duì)激光信號(hào)傳輸?shù)乃p影響很大。李娣等[1]認(rèn)為雨滴對(duì)激光信號(hào)的衰減主要取決于雨滴的尺寸分布、散射和吸收截面。研究表明,在可見光到遠(yuǎn)紅外光區(qū)約10μm的波長(zhǎng)范圍內(nèi),雨滴引起的激光衰減與光波長(zhǎng)的依賴關(guān)系不明顯。劉西川等[2]將雨滴的非球形效應(yīng)考慮在內(nèi),在雨滴的近似橢球模型、歸一化譜分布的基礎(chǔ)上,利用射線追蹤法計(jì)算群雨滴在可見光和近紅外波段的散射和衰減特性,討論分析雨滴譜分布和降雨強(qiáng)度對(duì)可見光和近紅外波段激光傳輸衰減的影響。Marshall和Palmer[3]提出了一種負(fù)指數(shù)分布模型即M-P分布。M-P譜分布適用于穩(wěn)定的層狀云降水,但是對(duì)強(qiáng)度較大降水?dāng)M合效果并不太好。郭婧等[4-5]以幾何光學(xué)散射理論為基礎(chǔ),研究了降雨引起激光的后向散射特性并給出不同降雨條件下的后向散射系數(shù),結(jié)合某型號(hào)激光定距系統(tǒng)計(jì)算了不同降雨條件和傳輸距離下的系統(tǒng)信噪比,分析了雨滴后向散射對(duì)系統(tǒng)信噪比的影響。汪亭玉[6]結(jié)合雨滴尺寸分布模型,建立了船舶大氣激光通信中雨滴遮擋作用的數(shù)理模型,對(duì)光在雨中的最大傳輸距離及其光強(qiáng)的衰減做了初步計(jì)算與實(shí)驗(yàn)。結(jié)果證明,該衰減模型是合理的,理論模型基本符合實(shí)驗(yàn)結(jié)果。該項(xiàng)研究為進(jìn)一步探索船舶大氣激光通信系統(tǒng)的性能奠定了基礎(chǔ)。

Weibull分布是一種從雨滴形成過程中嚴(yán)格推導(dǎo)出來的雨滴譜,它能夠充分反映各種雨滴尺寸的分布情況[7]。數(shù)值模擬結(jié)果表明,不同譜分布的群雨滴散射能力各不相同。激光衰減不僅受降雨強(qiáng)度的影響,也與雨滴的數(shù)密度有關(guān),數(shù)密度越大,激光衰減越大;而形狀對(duì)群雨滴的影響較小。但是,前人都是在脫離大氣環(huán)流背景來研究激光在雨中的衰減,而降水往往都是由于鋒面氣旋、臺(tái)風(fēng)等天氣系統(tǒng)造成的;因此本文利用1°×1°FNL氣象再分析資料和地面降水資料對(duì)2015年3月30日—4月2日的江淮氣旋的成因和環(huán)流形勢(shì)做簡(jiǎn)要分析,而后基于WRF中尺度氣象研究模式對(duì)該氣旋降水過程進(jìn)行模擬,通過Weibull雨滴尺寸分布模型得到激光在不同降水條件下的衰減特性(注:本文所說的時(shí)間均為北京時(shí))。

1 天氣過程概況

2015年3 月31日,受江淮氣旋的影響,江蘇全省平均降水量較常年同期異常偏多8.7%,降水大值區(qū)域位于蘇南地區(qū),其中徐州西部、長(zhǎng)江中部和蘇南南部地區(qū)偏多1~2.6倍。總體來說,該次江淮氣旋的形成過程是倒槽鋒生氣旋(圖略),開始時(shí),地面變性高壓東移入海后,由于高空南支鋒區(qū)上西南氣流將暖空氣向北輸送,地面減壓形成倒槽并東伸。這時(shí)在北支鋒區(qū)上有一小槽從西北移來,在地面上配合有一條冷鋒和鋒后冷高壓。而后,由于高空暖平流不斷增強(qiáng),地面倒槽進(jìn)一步發(fā)展并在槽中江淮地區(qū)有暖鋒生成,并形成了暖鋒。此時(shí),西北小槽繼續(xù)東移,南北兩支鋒區(qū)在江淮流域逐漸接近。冷鋒及其后部高壓也向東南移動(dòng),向倒槽靠近。最后,高空南北鋒區(qū)疊加,小槽發(fā)展,地面上冷鋒進(jìn)入倒槽與暖鋒接合,在高空槽前的正渦度平流下方,形成江淮氣旋。

2 江淮氣旋的發(fā)生發(fā)展機(jī)制

2.1 中高層環(huán)流形勢(shì)分析

2.1.1 200 hPa西風(fēng)急流分析 2015年3月30日08時(shí)(圖1(a)),200 hPa西風(fēng)急流中心位于45°N,急流橫跨15個(gè)緯度,南伸至35°N,氣旋的活動(dòng)與行星鋒區(qū)密切相關(guān):地面鋒線常常是極鋒行星鋒區(qū)在地面上的反映。氣旋的發(fā)生、發(fā)展一般都是在鋒區(qū)上進(jìn)行的,其出現(xiàn)的最大頻數(shù)以及主要路徑和鋒區(qū)的平均位置基本一致。在高空急流的南側(cè),強(qiáng)反氣旋式切變渦度造成的氣流輻散有利于地面氣旋的發(fā)展,在高空急流的北側(cè),強(qiáng)氣旋式渦度有利于地面反氣旋發(fā)展[8]。剛好該急流軸的南側(cè)為江蘇省上空,配合著江淮氣旋的發(fā)生與加強(qiáng)。

2.1.2 500 hPa環(huán)流分析 本次天氣過程屬于兩脊一槽型形成的江淮氣旋(圖1(b))。500 hPa圖上烏拉爾山附近為暖性高壓脊,中國(guó)沿海大陸也有一明顯的高壓脊,貝加爾湖、蒙古和中國(guó)北部地區(qū)為兩脊之間的大槽控制。南支小槽在偏西氣流的引導(dǎo)下向東移動(dòng),槽后為冷平流,槽前為暖平流,利于槽的發(fā)展,低空對(duì)應(yīng)氣旋的加強(qiáng)。

2.2 低空環(huán)流形勢(shì)分析

氣旋剛開始發(fā)展時(shí),700 hPa(圖1(c))同500 hPa高空一樣,也是兩脊一槽型,南支槽前輸送正渦度平流,使得槽得以發(fā)展,對(duì)應(yīng)地面氣旋的發(fā)展。與此同時(shí),850 hPa(圖1(d))上的南支小槽槽前為暖平流,槽后冷平流,等高線與等溫線有交角,說明槽線附近有斜壓不穩(wěn)定的趨勢(shì)。

圖1 2015年3月30日08時(shí)高低空環(huán)流形勢(shì)及物理量場(chǎng)分析

3 降水的發(fā)生機(jī)制分析

2015年4 月2日08時(shí)500 hPa上(圖2(a))該南支槽與東亞大槽合并,對(duì)應(yīng)地面氣旋入海,由于非絕熱加熱作用,海洋釋放凝結(jié)潛熱,提供能量,使得氣旋得以發(fā)展,在江蘇南部、朝鮮半島和日本形成大范圍降水。700 hPa上以上地區(qū)位于槽前濕區(qū)(圖2(b)),而且有江淮切變線的配合,水汽在該切變線附近輻合,降水地區(qū)位于地面鋒線的北部、700 hPa切變線以南的地區(qū)。這是因?yàn)?00 hPa切變線以南的偏南氣流一方面可以將南方的水汽不斷輸送過來,另一方面這股氣流沿著鋒面向上滑升,使水汽冷卻凝結(jié)成雨。而且該地區(qū)的強(qiáng)烈的垂直上升運(yùn)動(dòng)導(dǎo)致了江蘇一帶的強(qiáng)降水[9]。低空850 hPa有一支很強(qiáng)的西南急流(圖2(c)),我國(guó)與暴雨相聯(lián)系的西南風(fēng)低空急流存在于副熱帶高壓西側(cè)或北側(cè),低空急流多位于高空西風(fēng)急流入口區(qū)的右側(cè)或南亞高壓東部脊線附近。低空急流的形成與維持跟高低空環(huán)流的耦合發(fā)展有關(guān),是大氣環(huán)流演變的產(chǎn)物,同時(shí)與降水有密切的聯(lián)系。低空急流左側(cè)對(duì)應(yīng)降水區(qū),利于氣旋的發(fā)展;而且水汽通量很大,超過20×10-8g·cm-2hPa-1s-1,水汽通量散度為負(fù)值,低空水汽輻合,有垂直上升運(yùn)動(dòng),使得降水過程得以持續(xù)。該次降水過程就是高空西風(fēng)急流與低空急流的耦合的結(jié)果。

圖2 2015年4月2日08時(shí)高低空環(huán)流形勢(shì)及物理量場(chǎng)分析

綜上所述,擾動(dòng)小槽、高低空急流為氣旋發(fā)展加強(qiáng)提供動(dòng)力機(jī)制。西風(fēng)槽和高空急流隨氣旋同步東移,使得氣旋在東移中不斷加強(qiáng)。氣旋南側(cè)的低空西南風(fēng)急流以脈動(dòng)形式向其南部輸送水汽和不穩(wěn)定能量[10]。

4 數(shù)值模擬分析激光在氣旋降水中的衰減特性

網(wǎng)格設(shè)置為子母網(wǎng)格比率為3,在水平方向上:粗網(wǎng)格的分辨率30 km,細(xì)網(wǎng)格的分辨率10 km,垂直方向分為44層,采用η坐標(biāo)系,中心位置為38.63°N,122.90°E。微物理過程方案采用WSM3類簡(jiǎn)單冰方案,輻射方案采用rrtm長(zhǎng)波輻射方案,Dudhia短波輻射方案。近地面層方案采用Monin-Obukhov方案,陸面過程方案采用熱量擴(kuò)散方案,邊界層方案采用YSU方案,積云參數(shù)化方案采用淺對(duì)流Kain-Fritsch(new Eta)方案,初始邊界條件使用的是2015年3月30日08時(shí)—4月02日08時(shí)FNL再分析資料,積分時(shí)長(zhǎng)為72 h。

與此同時(shí),在選用擾動(dòng)邊界層和陸面物理過程時(shí)(sf_sfclay_physics=1)考慮了地面熱量和水汽通量,輻射光學(xué)厚度中考慮了云的影響(僅當(dāng)ra_sw_physics=1和ra_lw_physics=1時(shí)有效)。通過模擬和實(shí)測(cè)降水資料的對(duì)比相差不多,得出該參數(shù)化方案較為準(zhǔn)確,模擬效果可信。

4.1 雨滴衰減系數(shù)的計(jì)算

雨滴對(duì)激光信號(hào)的衰減主要取決于雨滴的尺寸分布、散射和吸收截面。雨滴的尺寸分布是指不同降雨率,不同尺寸的雨滴在空間的分布狀況,也稱為雨滴譜。雨滴的大小及在空間中的分布都是隨機(jī)的,一般用雨滴尺寸分布函數(shù)N(D)來描述雨滴的尺寸分布情況,單位為m-3·mm-1,D表示雨滴直徑,單位為mm。

降雨引起的光波信號(hào)的衰減系數(shù)可由式(1)計(jì)算得到:

式中:Dmin,Dmax分別為雨滴的最小和最大直徑;N(D)d D為單位體積空氣中,雨滴直徑在D和D+ d D之間的個(gè)數(shù),單位為m-3。Qe(x,m)為球形雨滴的總的衰減截面,它們都是粒子折射率和尺度參數(shù)的函數(shù)。因?yàn)橛甑蔚陌霃皆?0μm~5mm之間,遠(yuǎn)大于激光波長(zhǎng),根據(jù)Van de Hulst近似公式[11],可認(rèn)為衰減效率因子Qe(x,m)≈2。

雨滴尺寸分布函數(shù)N(D)利用Weibull分布,它的表達(dá)式為

式中:N0=1 000,η=0.95R0.14,σ=0.26R0.42,a=1。4.2 10.6μm激光在春季江淮氣旋降水中探測(cè)距離研究

對(duì)于不同的紅外武器系統(tǒng),加上系統(tǒng)的修正因子K,同時(shí)考慮各種因素的影響,其作用距離一般可寫成式(3)的形式[12]:

式中:Jλ1-λ2為波段范圍λ1-λ2內(nèi)的輻射強(qiáng)度值;為波段范圍內(nèi)的大氣透過率;τ0為光學(xué)系統(tǒng)的光譜透過率;D*為敏感元件在波段區(qū)間內(nèi)的平均光譜探測(cè)度;ω為系統(tǒng)的瞬時(shí)視場(chǎng);Δf為電路等效噪聲帶寬為系統(tǒng)的信噪比;K為系統(tǒng)的修正因子[13];D0(NA)為光學(xué)系統(tǒng)的有效入射孔徑的面積;NEFD為噪聲等效通量密度。

由式(3)可知,當(dāng)系統(tǒng)穩(wěn)定、目標(biāo)參數(shù)一定時(shí),紅外系統(tǒng)作用距離與大氣透過率有如下關(guān)系:

式中:K0為式(3)中各已知參數(shù)的計(jì)算值。

對(duì)于衰減率分布不均勻的降水區(qū)域,若已知降水衰減率的分布,便可利用式(4)求出探測(cè)距離。設(shè)a(d B/km)為降水衰減率,A(d B)為至距離d處的衰減度,則:

令:

由式(5)和式(6)可得:

則當(dāng)Φ(d)=K0時(shí),d=Dmax,為最大探測(cè)距離,即:

Φ(d)是一個(gè)單調(diào)增函數(shù),在探測(cè)距離以內(nèi),其值小于K0,即dDmax時(shí),Φ(d)>K0。設(shè)某型紅外探測(cè)裝備晴天下對(duì)某特定目標(biāo)的最大探測(cè)距離為14 km,測(cè)得此時(shí)大氣透過率為0.92,可估算出K0=14.596;K0的值也可以根據(jù)系統(tǒng)參數(shù)由式(3)計(jì)算得到。本文由式(1)計(jì)算了2015年4月01日12時(shí)10.6μm紅外波在降水區(qū)域的衰減率分布,設(shè)定紅外探測(cè)器經(jīng)緯度,由式(8)和式(9)求得每個(gè)方向(間隔1度)的探測(cè)距離,得到紅外系統(tǒng)探測(cè)距離評(píng)估圖象,如圖3所示。

圖3中間的十字交叉點(diǎn)表示紅外探測(cè)器的位置,外圈表示在正常情況下的探測(cè)距離,內(nèi)圈表示在不同降水強(qiáng)度情況下各個(gè)方向的探測(cè)距離。由圖3可知,在氣旋降水條件下,傳輸性能較好的10.6μm遠(yuǎn)紅外波仍然受到衰減,在小雨的情況下(圖3(a),降雨率為0.25mm/h)探測(cè)距離可損失10%;在中雨的情況下(圖3(b),降雨率為2.7mm/h)探測(cè)距離為正常情況下的70%;而在大雨情況下(圖3(c),降雨率為6mm/h)探測(cè)距離僅為正常情況下的50%。圖3的探測(cè)距離評(píng)估圖象中也可以將中間的點(diǎn)設(shè)為目標(biāo),則圖3中的內(nèi)圈表示為對(duì)目標(biāo)在各個(gè)方位的探測(cè)距離。

圖3 10.6μm紅外系統(tǒng)探測(cè)距離評(píng)估圖

5 結(jié)論

(1)該次江淮氣旋由于倒槽鋒生產(chǎn)生,江蘇地區(qū)和朝鮮半島、日本等的降水過程就是高空西風(fēng)急流與低空急流的耦合的結(jié)果,與此同時(shí)還有700 hPa切變線的配合。700 hPa切變線以南的偏南氣流一方面可以將南方的水汽不斷輸送過來,另一方面這股氣流沿著鋒面向上滑升,使水汽冷卻凝結(jié)成雨。

(2)在氣旋降水條件下,傳輸性能較好的10.6 μm遠(yuǎn)紅外波仍然受到衰減,在小雨的情況下(降雨率為0.25mm/h)探測(cè)距離可損失10%;在中雨的情況下(降雨率為2.7 mm/h)探測(cè)距離為正常情況下的70%;而在大雨情況下(降雨率為6mm/h)探測(cè)距離僅為正常情況下的50%。

[1]李娣,陳輝.激光大氣傳輸?shù)挠觎F衰減特性研究[J].電子設(shè)計(jì)工程,2011,09:1-5.

[2]劉西川,高太長(zhǎng),劉磊.基于非球形雨滴的降雨對(duì)激光傳輸衰減的影響[J].紅外與激光工程,2013,01:167-173.

[3]KimⅡ,McArthur B,Korevaar E.Comparison of Laser Beam Propagation at785 and 1550 nm in Fog and Haze for OpticalWireless Communications[J].OpticalWireless CommunicationsⅢ,2001(4214):26-37.

[4]郭婧,張合,王曉鋒.雨滴對(duì)波長(zhǎng)0.532μm激光光束的衰減特性[J].兵工學(xué)報(bào),2011,07:878-883.

[5]郭婧,張合,趙學(xué)健.降雨對(duì)激光定距系統(tǒng)后向散射特性的影響[J].光子學(xué)報(bào),2014,11:90-94.

[6]汪亭玉,陳家璧.雨滴在船舶大氣激光通信中的遮擋作用[J].光學(xué)精密工程,2014,10:2652-2658.

[7]馬春林,楊瑞科.激光衰減隨降雨率的變化關(guān)系[J].淮陰師范學(xué)院學(xué)報(bào):自然科學(xué)版,2009,04:302-305.

[8]朱乾根,林錦瑞,壽紹文等.天氣學(xué)原理和方法(第四版)[M].北京:氣象出版社,2007.

[9]徐娟,周春雨,高天赤.梅雨鋒江淮氣旋發(fā)展機(jī)制及其與暴雨的關(guān)系[J].科技通報(bào),2013,05:24-29+86.

[10]秦曾灝,李永平,黃立文.中國(guó)近海和西太平洋溫帶氣旋的氣候?qū)W研究[J].海洋學(xué)報(bào),2002(S1):106-112.

[11]Van de HulstH C.LightScatteringby SmallParticles[M].New York:Dover Publications,1981.

[12]王威,李進(jìn)杰,顧德均,高偉.紅外輻射的氣候環(huán)境傳輸特性分析[J].激光與紅外.2000,30(5):290-292.

[13]JAMilbrandt,M K Yau.A Multimoment Bulk Microphysics Parameterization.Part I:Analysis of the Role of the Spectral Parameter [J].Journalof the Atmosphere Sciences,2005,62(9):3051-3064.

Analysis on Laser Attenuation Characteristics in Spring Cyclonic Precipitation in the Yangtze and Huai River Basin

CHEN Kai-nuo1,2,SHAO Li-min2,WEIHai-liang3,DAIRen-wei4,AIYang5,MU Yang1
1.Department of Graduate Management,Dalian Naval Academy,Dalian 116018,Liaoning Province,China; 2.Department ofMilitary Oceanography,Dalian Naval Academy,Dalian 116018,Liaoning Province,China; 3.The 3rd Squadron of the Lvshun Equipment Technical Support Team,Troop 92956 of the PLA,Dalian 116041,Liaoning Province, China; 4.The Navy Zhoushan Maritime Garrison Command Headquarters Observatory,Zhoushan 316000,Zhejiang Province,China; 5.The 37th Team of the Staff,East China Sea Fleet Hydrologic Meteorological Center,Ningbo 315000,Zhejiang Province,China

In the rain,the performance of shipboard laser weapons and laser radars is constrained.Therefore,it is of great scientific and military significance to study the characteristics of laser attenuation in the rain.Based on the WRF mesoscale meteorological model,this paper analyzes and simulates the process of precipitation from March 30 to April 2,2015.It is found that the Jiang-huai cyclone was generated from inverted trough.The process of precipitation across Jiangsu Province as well as the Korean Peninsula and Japan resulted from the coupling of high altitude westerly jet and low-altitude jet.At the same time,the shear line on 700 hPa also played a role in this event.Under cyclone precipitation conditions,10.6μm far-infrared laser with good transmission properties also weaken.Under the drizzle (rain rate is 0.25 mm/h),the loss of detection range can reach 10%.In the case ofmoderate rain (rain rate is 2.7 mm/h),the detection range is only 70%of the normal, and in the case of the heavy rain(rainfall rate of 6mm/h)the detection range is only 50%of the normal.

cyclone in the Yangtze and Huai River Basin;low altitude jet stream;trough frontogenesis;laser attenuation coefficient;detection range

P732

A

1003-2029(2017)01-0092-05

10.3969/j.issn.1003-2029.2017.01.017

2016-04-06

陳凱諾(1991-),男,碩士研究生,主要從事海霧數(shù)值模擬和激光衰減特性研究。E-mail:837677852@qq.com

主站蜘蛛池模板: 欧美精品三级在线| 一级片免费网站| 日韩在线成年视频人网站观看| 日韩大乳视频中文字幕| 就去吻亚洲精品国产欧美| 美女无遮挡免费视频网站| 丁香婷婷激情综合激情| 国产亚洲欧美日韩在线一区| 天堂中文在线资源| 国产成人精品一区二区三区| 国产在线观看91精品亚瑟| 国产va在线观看免费| 国产自在线拍| 老司国产精品视频91| 国产在线八区| 热re99久久精品国99热| 国产精品太粉嫩高中在线观看| 欧美一区二区三区不卡免费| 国产精品视频第一专区| 国产成人精品男人的天堂下载| 久夜色精品国产噜噜| 色天堂无毒不卡| 国产产在线精品亚洲aavv| 亚洲色精品国产一区二区三区| 超碰91免费人妻| 一级毛片基地| 99国产在线视频| 青草娱乐极品免费视频| aa级毛片毛片免费观看久| 丁香婷婷在线视频| 亚洲第七页| 亚洲一本大道在线| 亚洲av日韩综合一区尤物| 国产一区二区免费播放| 天堂网亚洲系列亚洲系列| 国模极品一区二区三区| 欧美怡红院视频一区二区三区| 91久久偷偷做嫩草影院电| 欧美在线视频不卡| 日韩A∨精品日韩精品无码| 亚洲最大情网站在线观看 | 国产制服丝袜91在线| 毛片三级在线观看| 精品人妻一区二区三区蜜桃AⅤ| 国产精品伦视频观看免费| 欧美精品成人| 91丝袜在线观看| 狠狠干综合| 亚洲制服丝袜第一页| 成人亚洲天堂| 全午夜免费一级毛片| 日韩精品成人在线| 2020最新国产精品视频| 自拍亚洲欧美精品| 亚洲视频二| 最新午夜男女福利片视频| 国产成人精品一区二区三在线观看| 91视频日本| 中文字幕丝袜一区二区| 国产偷倩视频| 免费毛片a| 中文成人在线| 亚洲国产欧美自拍| 亚洲激情区| 香蕉久久国产超碰青草| 国产精品不卡片视频免费观看| 第一区免费在线观看| 四虎影视国产精品| 国产极品嫩模在线观看91| 国产一级在线观看www色 | 欧类av怡春院| 91日本在线观看亚洲精品| 操国产美女| 国产精品免费福利久久播放| 国产成人成人一区二区| 无套av在线| 国产精品亚欧美一区二区三区| 毛片免费网址| 黄片一区二区三区| 国产一区二区三区免费观看| 97在线免费视频| 亚洲国产精品无码AV|