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

高海拔高溫差深切峽谷橋址區日常大風成因

2015-01-07 07:59:28李永樂張明金徐昕宇陶齊宇朱樂東宋麗莉
西南交通大學學報 2015年4期
關鍵詞:風速

李永樂, 張明金, 徐昕宇, 陶齊宇, 朱樂東, 宋麗莉

(1.西南交通大學橋梁工程系,四川成都610031;2.四川省交通運輸廳公路規劃勘察設計研究院,四川成都610041;3.同濟大學橋梁工程系,上海200092;4.中國氣象局公共氣象服務中心,北京100081)

高海拔高溫差深切峽谷橋址區日常大風成因

李永樂1, 張明金1, 徐昕宇1, 陶齊宇2, 朱樂東3, 宋麗莉4

(1.西南交通大學橋梁工程系,四川成都610031;2.四川省交通運輸廳公路規劃勘察設計研究院,四川成都610041;3.同濟大學橋梁工程系,上海200092;4.中國氣象局公共氣象服務中心,北京100081)

為探討高海拔高溫差深切峽谷橋址區日常大風的成因,采用CAW600-RT型四要素自動氣象站、手持風速儀及便攜式溫度計,對大渡河大橋橋址區風特性進行實測,分析了橋位處平均風速與溫度、日照及地形地貌等的相關性.結果表明:大渡河大橋位于高海拔高溫差深切峽谷內,橋址區幾乎每天下午起風,平均風速常達10 m/s以上;根據成因,橋位處的大風可分為2類,一類受大尺度大氣環流影響,另一類受小尺度范圍內熱力驅動而產生日常大風,并受局部地形及隨時間變化的日照的影響;橋位處日常大風出現的頻率較高,雖不控制橋梁的設計基準風速,但影響橋梁的耐久性和行車舒適性.

深切峽谷;橋址區;日常大風;現場實測;成因分析

Key words:deep gorge;bridge site;daily strong wind;field measurement;cause analysis

通常,大風過程對國民經濟各個方面都會產生較大影響,因此,研究大風的成因具有重要的現實意義,也得到了許多氣象學者和相關專家的極大關注[1-4].目前,針對復雜地形地貌區風場的研究主要有4種手段,即理論推導、數值模擬、模型試驗和現場觀測.其中,理論推導是對復雜地形進行一定簡化,建立空氣運動方程和熱傳導方程,通過求解微分方程組,得到相關的風場特征[1].數值模擬通常是以理論推導為基礎,借助計算機求解的一種分析方法[2].

無論是理論推導還是數值模擬,均采用了一系列假定(邊界條件)和計算模型,而這樣的簡化會導致一定程度的失真[3-4].模型試驗是借助風洞實驗室對關心的地形地貌進行模擬試驗,但實驗室中風場的模擬也具有一定的局限性.

現場觀測是目前較有效也是采用較多的一種手段.馬玉堂、陳凱、張人文、宋麗莉、朱樂東、王凱等通過在山區建立風觀測站,對復雜山區地形的風場特性進行了觀測分析[5-10].李永樂、王文勇等采用計算機數值模擬方法,對山區風特性的分析表明,山區峽谷地形條件下的風環境較復雜,且不同橋址區的風特性差異較大[11-13].譚波、廖曉春等通過現場實測發現,在西部高海拔高溫差山區,每天下午均會出現規律性波動的大風,部分地區風速可達10 m/s以上[14-15],但相關文獻中均沒有對這種風場形成的原因進行討論.

山區峽谷地形在我國中西部地區較為常見,其地形復雜多變,現有相關規范、文獻中對山區風特性及成因的描述較少,對高海拔高溫差深切峽谷橋址區風特性的研究尚未見報道.

本文以位于四川西部的大渡河大橋為工程背景,通過在大橋橋址區進行現場實測,取得了一些關于高海拔高溫差深切峽谷地形風場特性的資料,通過對相關資料進行分析,探討了日常大風的特點及其成因,研究結果對類似地貌風特性的分析有一定參考意義.

1 觀測概況

1.1 地形地貌

大渡河大橋距離瀘定縣城約5 km.大橋地處高山峽谷之間,橋面距離大渡河溝底約300 m,連接橋梁兩端的均是陡峭的山脈.橋軸線向康定側延伸10 km后,地面海拔由橋位處的1 608 m升高到約4 500 m;橋軸線向雅安側延伸5 km后,地面海拔由橋位處的1 608 m升高到約3 700 m.峽谷兩側10 km范圍內均有終年不化的雪山,而峽谷內屬于典型的干熱河谷氣候,溫暖干燥,峽底和頂部溫差較大.

此外,橋位處晝夜溫差較大,加之橋位處海拔較高,地形呈現明顯的深切峽谷特性,影響風場的因素較多.

1.2 觀測站點設置和試驗儀器

根據大橋所處的位置走向以及橋位附近地區的地形特點,在橋位處安裝了1套四要素自動氣象站(CAW600-RT).自動氣象站位于橋址區風速較大的咱里村大風崗上,該觀測站點基本位于大橋縱向中心軸線上,在大橋上游約30 m處,大橋跨中偏離康定側橋塔約100 m.風速、風向傳感器的海拔高程為1 530 m,距離橋面設計高度78 m.儀器安裝和測點布置見圖1.

圖1 儀器安裝和測點布置Fig.1 Equipment installation and measuring points'layout

1.3 沿河谷溫度和風速測量

為了考察橋位處河谷溫度場和風場的分布情況,以橋位處為界,將河谷分為上游河谷和下游河谷.在上游河谷,沿國道(G318)走向布置9個固定

圖2 溫度測點布置Fig.2 Layout of temperature measuring points

2 日常大風成因分析

冬季和夏季典型大風天10 min平均風速與溫度的變化趨勢分別見圖3(a)和圖3(b).可見,無論是冬季還是夏季,風速都呈現出以天為周期的規律性波動,這種規律性波動在夏季更明顯.

圖3 10 min平均風速與溫度F20ig13.年3 8月M9日ean wind speed and temperature in 10 min

對比圖3(a)和圖3(b)可見,橋址區的大風過程可以分為兩大類:一類是冬季出現頻率較低的大風降溫過程,定義這類大風降溫過程為第Ⅰ類大風,持續時間一般可達2 d左右;另一類是夏季出現頻率較高、以天為周期的日常大風過程,定義這類大風為第Ⅱ類大風,其10 min平均風速可以達到10 m/s,對橋上車輛的正常行駛有一定影響.下面重點就第Ⅱ類日常大風的成因進行討論.

2.1 溫度影響

2.1.1 風速、溫度隨時間的變化

自動氣象站從2012年12月起開始采集數據,時間跨度包括冬季和夏季,具有較好的代表性.圖4為2013年2月典型日常大風天平均風速變化趨勢.

圖4 典型大風天10 min平均風速Fig.4 Typical mean wind speed in 10 min

可見,每天凌晨至中午時段內風速較小,在每天下午至上半夜風速較大.為分析橋位處每天的起風規律,分別對10 min平均風速大于5和10 m/s的天數進行統計分析.觀測站處10 min平均風速大于5 m/s的大風過程共143 d,10 min平均風速大于10 m/s的大風過程共54 d.

圖5(a)為每天10 min平均風速大于5 m/s的大風起止時間,大風起風時間基本上都在中午12點左右,止風時間約在晚上10點,10 min平均風速大于5 m/s的大風過程平均持續時間約9 h.

圖5(b)為每天10 min平均風速大于10 m/s的大風起止時間和持續時間,大風起風時間基本上都在下午2點半,晚8點半左右停止,平均持續時間約6 h.

從圖3(a)和圖3(b)可見,風速和溫度以天為周期的規律性波動明顯,并且二者的波動趨勢基本一致.每天風速較大的時段溫度也較高,由此說明熱力作用是大橋橋位處日常大風形成的一個重要因素.

圖5 大風起止時間Fig.5 Starting and ending times of strong wind

2.1.2 風速、溫度沿河谷的變化

為了考察橋位處河谷溫度場和風場的分布情況,沿河谷布置了18個測點,2013年4月19日下午對固定測量點相同時刻的溫度和風速進行了3次測量.測量儀器采用溫度計和手持風速儀,在開始測量前和測量結束后,對采用的溫度計進行統一標定和修正.測量過程中,為減小地面熱輻射引起的誤差,溫度計離地高度在2.0 m以上.

不同時刻測點的溫度和風速變化見圖6.可見,測點N9和S1的溫度在3次測量過程中均比其他測點低,以這2點為中心分別往上、下游延伸,溫度均出現上升的趨勢.

對比3次測量的溫度曲線可以看出,越接近晚上,上游側河谷(測點N5~N9)的溫度降低越慢,下游側河谷(測點N9~S9)的溫度下降越快.

分析同一時刻測點風速的變化規律,溫度較低的測點N8、N9和S1在3次測量中的風速均比兩側測點的大,而溫度較高兩側測點的風速相對較低.這表明,同一時刻沿河谷的溫差是形成日常大風的重要原因.

2.2 局部地形影響

大橋橋址區下午4點以前,整個河谷內風向通常是以由南向北的東南風為主;下午4點之后,橋位上游河谷仍然是以沿河道向上游方向吹的東南風為主,但橋位下游水壩以下河谷的風向則與上游河谷相反.圖7為橋位區域地形的三維視圖,可以看出,橋位處地形較復雜,河谷東西岸均有山脊,并且在大壩處河谷斷面收縮較大.

圖6 溫度和風速Fig.6 Temperature and wind speed at different measuring points

結合局部地形對河谷內大風風向的分布進行分析:由于東岸雪山上的來流同時受到東岸山脊和水電站大壩處西岸山脊的雙重阻擋,在水壩附近形成了一個較大的氣流漩渦,且由于山峰表面的不均勻性和主導來流風向的不確定性,該氣旋不一定是平面水平氣旋,應是一個立體氣旋;氣流漩渦導致橋位下游咱里村處風向較亂,而且在橋位附近的高空風速也不一定比低空風速大,甚至出現高空風速小于地面局部峽谷口和山峰處的風速;越接近晚上,瀘定縣城往上游的來流風增大程度不如東岸雪山上來流增大程度大,下午4點以后河谷內的風向以瀘定水電站大壩為分界,水壩以上河谷內以東南風為主,且風速較大,水壩以下的瀘定縣城以西北風為主,風速相對較小.這也是瀘定縣氣象站歷史記錄中北風偏多的原因.鑒于橋位處大風的風速和風向均在較大程度上受到橋位處局部地形的影響,因此,利用瀘定縣已有的基本氣象臺站的歷史數據時,應充分考慮到地形可能造成的影響.

圖7 橋位區域地形三維視圖Fig.7 3-D view of the terrain on the bridge site

2.3 日照影響

在上午晴朗少云的天氣下,午后常出現東南方向的大風,這是由于太陽照射對橋位處東南側、北側不規則山體表面輻射增溫隨時間變化不均勻,致使各個山坡面、山頂、山谷底部出現較大的溫度變化率和氣壓梯度.具體表現為上游河谷在下午5點以后陽光仍然可以直射,而此時在橋位處已經不能被陽光照射.實測表明,橋位處同一時刻的溫度比上游河谷低3~5℃.橋位下游5 km處為瀘定縣城,縣城內水泥路面、建筑物等較多,在整個白天的日照過程中吸收的熱能也較多,并且混凝土散熱較慢,使得瀘定縣城出現一定的熱島效應,導致橋位處同一時刻的溫度也比下游河谷低2~4℃.在該熱力溫差的驅動下,河谷這種非均勻下墊面的熱力、動力共同作用導致山谷內形成不穩定的擾動氣流,從而形成以天為周期的大風.

圖8為同一時刻上游河谷測點N1和下游河谷測點S1的照片,可見,上、下游河谷陽光照射情況差別明顯.

圖8 同一時刻(17:00)不同測點日照情況Fig.8 Sunlight at different points at 5:00 p.m.

2.4 綜合分析

通過觀察和分析發現,與常規的季風或強對流天氣出現的大風不同,橋位處每天下午出現的日常大風是由于局部溫差、局部地形和日照不均勻共同作用形成的局部小尺度大風,不平衡熱力作用對大風的產生及其強弱起控制作用.已有的現場觀測也表明,溫度越高、溫差越大,風速就越大.沒有較強的天氣系統影響時,橋位區域早晨、上午及夜間的風速較小.當河谷東岸雪山上來流冷空氣足夠強時,河谷內會出現以水電站大壩為分界線的反向風.

通過大量實測和現場考察,基本弄清了該橋位處風場流向,見圖9.可見,河谷內的大風主要是兩岸山坡上的山風和河谷中局部溫差形成的風.越接近傍晚,兩岸山坡上的來流越強,當東岸山坡上的來流強度超過河谷內本身向上游流動的強度時,就出現以水電站大壩為分界的風向截然相反的情況.此時,在大壩上游側河谷內的風是河谷風、山坡風和局部地形加速效應的三重疊加,所以在橋位附近區域出現了一個明顯的大風區.大壩下游沒有明顯的河谷加速效應,同時山坡風和河谷風的流動方向相反,河谷中的風速是山坡風與河谷風相互抵消后形成的,故大壩下游側河谷內的風速明顯小于大壩上游側河谷.

圖9 橋址區風場示意Fig.9 Sketch map of wind field on the bridge site

3 結 論

針對大渡河大橋橋址區幾乎每天下午出現的大風,通過現場實測,分析了高海拔高溫差深切峽谷區日常大風的成因,得到以下結論:

(1)橋址區大風可以分為2類,一類是受大尺度氣候環境影響,定義為第Ⅰ類大風;另一類是受小尺度范圍內熱力驅動而產生的日常大風,并受局部地形及隨時間變化的日照的影響,定義為第Ⅱ類大風.

(2)日常大風(第Ⅱ類大風)出現頻率較高,基本上每天下午都會出現,特別是在夏季,這類大風每天持續時間約6 h,其10 min平均風速一般在10 m/s左右,對日常行車有一定影響,但不控制大橋的設計基準風速.

(3)每天不同時段太陽對河谷同一位置的日照不同,形成以天為周期的溫度波動,此外,同一時刻對不同地點的日照不同,使得河谷內出現較大的局部溫差,因此,熱力不平衡作用是大橋橋位處日常大風形成的一個重要原因.

(4)由于橋位處地形的特殊性,導致以水電站大壩為分界的風向可能截然相反,此時在大壩上游河谷內的風是河谷風、山坡風和局部地形加速效應的疊加,所以在橋位附近區域出現了一個明顯的大風區.大壩下游沒有明顯的河谷加速效應,山坡風和河谷風的流動方向相反,河谷中的風因山坡風和河谷風相互抵消,所以大壩下游河谷內的風速明顯小于大壩上游河谷.

(5)橋址區日常大風的產生機理與常規地形不同,加上局部地形的影響,橋址區日常大風的風特性可能與抗風規范及通常研究的風特性存在顯著差異.

[1] 傅抱璞.山谷風[J].氣象科學,1980,1(2):1-14.FU Baopu.Mountain and valley breezes[J].Scientia Meteorological Sinica,1980,1(2):1-14.

[2] 余錦華,傅抱璞.山谷地形對盛行氣流影響的數值模擬[J].氣象學報,1995,53(1):50-61.YU Jinhua,FU Baopu.Numerical simulation of the airflow over and in valley terrain[J].Acta Meteorological Sinica,1995,53(1):50-61.

[3] 傅抱璞.山地氣候文集[M].北京:氣象出版社,1984:15-22.

[4] 徐大海.大氣邊界層內風的若干特性及其應用[J].空氣動力學學報,1984,2(3):75-87.XUDahai.Windcharacteristicsintheplanetary boundary layer and applications[J].Academy of Meteorological Science,1984,2(3):75-87.

[5] 馬玉堂,徐兆生.青藏高原大氣貼地層的基本物理特點[J].氣象學報,1987,45(2):210-217.MAYutang,XUZhaosheng.Thebasicfhysical characteristics of the atmosphere near the ground over the Qinghai-Xizang Plateau[J].Acta Meteorological Sinica,1987,45(2):210-217.

[6] 陳凱,余永生,賈叢賢.傍山地區的強風場特性實測研究[J].流體力學實驗與測量,2003,17(3):18-22.CHEN Kai,YU Yongsheng,JIA Congxian.Full scale measurement of a strong wind near mountains[J].Experiments and Measurements in Fluid Mechanics,2003,17(3):18-22.

[7] 張人文,范紹佳,李穎敏.2008年秋季從化山谷風觀測研究[J].熱帶氣象學報,2012,28(1):134-139.ZHANG Renwen,FAN Shaojia,LI Yingmin.Observational study of mountain-valley breeze over Conghua in 2008 fall[J].Journal of Tropical Meteorology,2012,28(1):134-139.

[8] 宋麗莉,吳戰平,秦鵬,等.復雜山地近地層強風特性分析[J].氣象學報,2009,67(3):452-460.SONG Lili,WU Zhanping,QIN Peng,et al.An analysis of the characteristics of strong winds in the surface layer overacomplex terrain[J].ActaMeteorologica Sinica,2009,67(3):452-460.

[9] 朱樂東,任鵬杰,陳偉,等.壩陵河大橋橋位深切峽谷風剖面實測研究[J].實驗流體力學,2011,25(4):15-21.ZHU Ledong,REN Pengjie,CHEN Wei,et al.Investigation on wind profiles in the deep gorge at the Balinghe bridge site via field measurement[J].Journal of Experiments in Fluid Mechanics,2011,25(4):15-21.

[10] 王凱,廖海黎,李明水,等.山區峽谷橋梁設計基準風速的確定方法[J].西南交通大學學報,2013,48(1):29-35.WANG Kai,LIAOHaili,LIMingshui,etal.Determination method for basic design wind speed of mountainous-valley bridge[J].Journal of Southwest Jiaotong University,2013,48(1):29-35.

[11] 李永樂,蔡憲棠,唐康,等.深切峽谷橋址區風場空間分布特性的數值模擬研究[J].土木工程學報,2011,44(2):116-122.LI Yongle,CAI Xiantang,TANG Kang,et al.Study of spatial distribution feature of wind fields over bridge sitewithadeep-cuttinggorgeusingnumerical simulation[J].ChinaCivilEngineeringJournal,2011,44(2):116-122.

[12] 李永樂,唐康,蔡憲棠,等.深切峽谷區大跨度橋梁的復合風速標準[J].西南交通大學學報,2010,45(2):167-173.LI Yongle,TANG Kang,CAI Xiantang,et al.Integrated wind speed standard for long-span bridges over deep-cuttinggorge[J].JournalofSouthwestJiaotong University,2010,45(2):167-173.

[13] 王文勇,陳楠.山區風場的WRF/CALMET模式模擬[J].西南交通大學學報,2010,45(6):990-996.WANG Wenyong,CHEN Nan.WRF/CALMET-based wind field simulation of mountain regions[J].Journal of Southwest Jiaotong University,2010,45(6):990-996.

[14] 譚波,趙志軍.林芝機場地面大風特征分析及其對飛行的影響[J].資源開發與市場,2010,26(12):1071-1073.

TAN Bo,ZHAO Zhijun.Effects of wind on flight in Linzhi airport[J].Resource Development&Market,2010,26(12):1071-1073.

[15] 廖曉春.九黃機場側風分布特征及其機理研究[D].廣漢:中國民用航空飛行學院空中交通管理學院,2012.

(中、英文編輯:付國彬)

Causes of Daily Strong Wind on Bridge Site in Deep Gorge with High Altitude and High Temperature Difference

LI Yongle1, ZHANG Mingjin1, XU Xinyu1, TAO Qiyu2, ZHU Ledong3, SONG Lili4
(1.Department of Bridge Engineering,Southwest Jiaotong University,Chengdu 610031,China;2.Sichuan Provincial Transport Department Highway Planning,Survey,Design and Research Institute,Chengdu 610041,China;
3.Department of Bridge Engineering,Tongji University,Shanghai 200092,China;4.Public Weather Service Center of China Meteorological Administration,Beijing 100081,China)

In order to investigate the causes of daily strong wind on a bridge site in a deep gorge with a high altitude and a high temperature difference,the field measurement of wind characteristics was carried out on the Dadu River bridge site with the CAW600-RT automatic-weather-station(AWS),hand anemometer and portable thermometer.The property of wind environment on the bridge site was explored through analyzing its influence factors,including mean wind speed,temperature,sunlight,topography and so on.The research results show that the Dadu River bridge is located in a deep gorge with a high altitude and a high temperature difference,where strong wind occurs nearly every afternoon and its mean wind speed is often more than 10 m/s.The causes of strong wind can be classified into two types.The first type of wind is induced by large-scale general atmospheric circulation,while the second type of wind,i.e.daily strong wind,is driven by local thermal disequilibrium,and the daily strong wind is also influenced by local terrain and time-varying sunshine.Though frequent daily strong wind has not effect on the design wind speed of the bridge,it can affect bridge durability and driving comfort to some extent.

U442

:A

0258-2724(2014)06-0935-07

10.3969/j.issn.0258-2724.2014.06.001

2013-12-12

國家自然科學基金資助項目(U1334201);交通運輸部建設科技計劃項目(2014318800240)

李永樂(1972-),男,教授,博士,博士生導師,研究方向為大跨度橋梁風致振動及車橋耦合振動,電話:028-87601119,E-mail:lele@swjtu.edu.cn

李永樂,張明金,徐昕宇,等.高海拔高溫差深切峽谷橋址區日常大風成因[J].西南交通大學學報,2014,49(6):935-941.

猜你喜歡
風速
邯鄲市近46年風向風速特征分析
基于Kmeans-VMD-LSTM的短期風速預測
基于最優TS評分和頻率匹配的江蘇近海風速訂正
海洋通報(2020年5期)2021-01-14 09:26:54
基于時間相關性的風速威布爾分布優化方法
陜西黃土高原地區日極大風速的統計推算方法
陜西氣象(2020年2期)2020-06-08 00:54:38
基于GARCH的短時風速預測方法
快速評估風電場50年一遇最大風速的算法
風能(2016年11期)2016-03-04 05:24:00
考慮風切和塔影效應的風力機風速模型
電測與儀表(2015年8期)2015-04-09 11:50:06
GE在中國發布2.3-116低風速智能風機
考慮風速分布與日非平穩性的風速數據預處理方法研究
主站蜘蛛池模板: 国产91精选在线观看| 91精品情国产情侣高潮对白蜜| 亚洲综合极品香蕉久久网| 香蕉久久国产精品免| 91福利国产成人精品导航| 91热爆在线| 国产自产视频一区二区三区| 国产永久无码观看在线| 国产一级片网址| 久久精品女人天堂aaa| 国产在线无码av完整版在线观看| 色综合国产| 精品伊人久久久久7777人| 国产成人高精品免费视频| 人人看人人鲁狠狠高清| 欧美午夜网站| 亚洲另类国产欧美一区二区| 亚洲精品无码成人片在线观看 | 精品国产香蕉伊思人在线| 欧美一级高清片久久99| 国产成人精品18| 国产精品视频第一专区| 欧美激情成人网| 日韩欧美视频第一区在线观看| 亚洲a级在线观看| 日韩成人在线网站| 色欲不卡无码一区二区| 伊人久久久大香线蕉综合直播| a欧美在线| 97人人做人人爽香蕉精品| 91av国产在线| 无码综合天天久久综合网| 免费一级成人毛片| 全部无卡免费的毛片在线看| 成人午夜精品一级毛片| 欲色天天综合网| 中国精品自拍| 五月婷婷激情四射| 中文字幕在线免费看| 高清无码手机在线观看| 国产视频久久久久| 亚洲va视频| 黄色网在线| 亚洲狼网站狼狼鲁亚洲下载| 精品无码一区二区三区在线视频| 日韩在线观看网站| 高清亚洲欧美在线看| 99re热精品视频国产免费| 国产黄色爱视频| 欧美无遮挡国产欧美另类| 亚洲三级色| 国产一区二区三区视频| 日本在线欧美在线| 亚洲男人天堂网址| 成人午夜视频在线| 国产日韩欧美在线播放| 伊在人亚洲香蕉精品播放| 91人妻日韩人妻无码专区精品| 国产精品免费露脸视频| 国产永久在线视频| 操国产美女| 91精品福利自产拍在线观看| 国产日韩精品欧美一区喷| 久久这里只有精品2| 亚洲天堂日韩在线| 久久精品日日躁夜夜躁欧美| 91成人免费观看在线观看| 国产精品黑色丝袜的老师| 片在线无码观看| 国产AV无码专区亚洲A∨毛片| 91麻豆精品国产91久久久久| 国产成人调教在线视频| 奇米影视狠狠精品7777| 456亚洲人成高清在线| 婷婷成人综合| 日本一区二区三区精品国产| 亚洲天堂视频在线观看| 成人福利在线免费观看| 国产一在线观看| 无码福利视频| 中文字幕日韩丝袜一区| 国产成人AV大片大片在线播放 |