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

    大連海上機場施工通道開口寬度對水環(huán)境影響的數(shù)值研究

    2014-12-05 02:17:18張寧川
    水道港口 2014年1期
    關(guān)鍵詞:模型施工

    猶 爽,張寧川

    (大連理工大學(xué)海岸和近海工程國家重點實驗室,大連 116024)

    大連海上機場位于渤海遼東灣南部金州灣海域,是我國內(nèi)地即將修建的第一座離岸式海上機場,竣工后將成為目前為止世界上最大的海上機場。機場人工島在軟土基礎(chǔ)上填筑形成[1],規(guī)劃人工島總面積為21 km2,平面布置為矩形,寬3.4 km,長6.4 km。機場物流區(qū)布置在機場人工島西側(cè)海岸,機場施工過程中擬修建從空港物流中心到機場人工島西南角的全封閉式施工通道。機場人工島填海完成后,為保障水質(zhì)環(huán)境,施工通道擬全部或部分拆除。

    金州灣距離大連市金州區(qū)約2 km,海灣呈橢圓形,灣口朝向西北,水域面積約為342 km2,海岸線長65.7 km,灣內(nèi)水深自東向西遞增至約10 m。金州灣是金州區(qū)工業(yè)和生活廢水以及北大河、紅旗河等流域的最終歸處。因為機場人工島位于近灣底水域,加之施工通道,在灣底形成半封閉水域。該水域水質(zhì)保障成為機場建設(shè)環(huán)境保護面臨的挑戰(zhàn)性課題。

    為保障大連海上機場建成后金州灣海域水流形態(tài)和污染物輸移順暢,進行海上機場施工通道布置的水動力優(yōu)化研究,期望通過優(yōu)化海上機場施工通道開口寬度,盡可能減小工程建設(shè)對金州灣海域的水動力及環(huán)境影響?;诖四康?,以工程水域物理自凈能力、水體半更換期、局部污染物擴散軌跡、局部污染物示蹤粒子第一次到達灣外的時長等參數(shù)為評價指標,建立非結(jié)構(gòu)網(wǎng)格下的水動力及污染物輸移擴散數(shù)值模型,對大連海上機場施工通道布置方案進行水動力學(xué)優(yōu)化研究。

    1 水動力模型的建立及驗證

    1.1 模型計算范圍及網(wǎng)格劃分

    金州灣位于遼南半島前端黃渤海分界線附近的渤海一側(cè)水域,海灣呈橢圓形,灣口朝向西北開放。漲潮水動力主要來自北黃海漲潮流通過遼南半島老鐵山水道的水流動力;落潮水動力主要來自渤海遼東灣落潮流動力。灣內(nèi)潮流受進出灣口的水動力控制。灣口2個岬角地理坐標分別為39°03′30"N、121°22′25"E和39°11′42"N、121°34′50"E(圖 1)。

    人工島工程位于金州灣灣底,距南岸線最短距離僅2 km,距東岸線最短距離約4 km,長軸走向為131°~311°。人工島西南角與西側(cè)空港物流中心由施工通道連接(圖 2)。

    針對上述地理形勢和工程布置方案,考慮到開邊界對工程所在海域潮流計算的影響,將模型計算區(qū)域開邊界設(shè)置在距離工程區(qū)域西側(cè)約30 km附近的2個岬角南北垂直連線處。計算范圍東西跨度約58 km,南北跨度約55 km,計算范圍的經(jīng)緯度坐標為121.266°~121.927°E;38.950°~39.446°N。

    為更好地擬合岸線,網(wǎng)格劃分采用三角形網(wǎng)格單元[2-4],并對重點關(guān)注區(qū)域進行加密,網(wǎng)格單元總數(shù)為19 491,節(jié)點總數(shù)為10 187,最小網(wǎng)格步長為100 m。

    圖1 金州灣地理位置及測點布置Fig.1 Geographic position of Jinzhou Bay and layout of sites

    圖2 機場人工島及物流中心平面規(guī)劃布置Fig.2 Plane layout of artificial island and logistics center of airport

    1.2 計算模型選取及參數(shù)設(shè)置

    計算模型選取MIKE21中的flow model FM模型。該模型采用淺水連續(xù)性方程和動量方程作為控制方程,適用于近岸海域下的水動力計算。本研究不考慮溫度和鹽度影響,模型為正壓模型。

    計算時,模型渦黏性系數(shù)通過Smagorinsky方程確定,取Smagorinsky系數(shù)Cs為0.28;通常情況下阻力系數(shù)Cf可由曼寧系數(shù)M或謝才系數(shù)C確定,本模型采用曼寧系數(shù)來決定底摩阻,鑒于計算區(qū)域范圍內(nèi)水深和地質(zhì)類型變化不大,故將整個計算范圍內(nèi)的曼寧系數(shù)值設(shè)置為常數(shù),取值為50;MIKE21采用“凍結(jié)法”處理邊界,通過定義臨界水深來確定干濕單元[5],本研究中干濕單元的臨界水深取為0.05 m。

    模型計算歷時長度取為34 d,啟動時間設(shè)為1 d,模型時間步長為10 s。計算得到的流場如圖3所示。

    圖3 金州灣海域流場Fig.3 Flow field of Jinzhou Bay

    1.3 模型邊界條件及初始條件

    模型計算區(qū)域驗潮站有限,難以提供開邊界資料,本文采用大小區(qū)域嵌套法獲取開邊界潮位。大區(qū)域計算范圍選取潮位資料豐富的整個渤海海域,驗潮站位包括煙臺、蓬萊港、龍口港、老虎灘、長興島等10個站點,開邊界為威?!¢L山島連線及其延長線與岸線相交部分,大范圍模型計算結(jié)果與實測潮位吻合較好,可以為小范圍模型提供可靠的開邊界條件。初始水位及流速均設(shè)置為0。

    1.4 模型驗證

    選取靠近岸邊的棋盤磨潮位驗證點(實測資料為2009年11月7日~2009年12月11日)和1#~6#潮流驗證點分別對模型的潮位和潮流計算結(jié)果進行驗證,各驗證點布置如圖1所示。潮流驗證分大潮時段(2009年12月3日~4日)和小潮時段(2009年11月27日~28日)。

    驗證結(jié)果如圖4所示,這里潮流驗證以2#點和5#點作為示例,從圖4中可以看出,計算值和實測值吻合較好,具有良好的一致性,由于本模型沒有考慮溫度、鹽度及風(fēng)的影響,在個別時刻實測值和計算值具有一定差別,但總體符合較好,且能正確地刻畫漲落潮的情況,通過潮位和流場的驗證結(jié)果表明,模型是合理的,模擬結(jié)果可以反映計算海域的水動力場特征。

    圖4 各站點潮位和潮流驗證Fig.4 Tidal level and flow verification of each site

    2 施工通道開口寬度對水環(huán)境的影響

    2.1 施工通道開口寬度設(shè)定

    工程水域人工島—施工通道—海岸之間形成的半封閉水體與外海的交換,主要取決于施工通道的開口寬度和人工島—海岸之間形成的流口寬度B(圖2)。在人工島寬度和走向確定的條件下,人工島—海岸之間形成的流口寬度B為一確定值。為了使海灣自凈能力、局部水體交換率等指標達到最優(yōu)值,擬將連接人工島與陸地的施工通道從靠近人工島一側(cè)進行開挖,并確定最佳開挖寬度。

    施工通道總長度約3 000 m(記為L,圖2)。為了確定最佳開口寬度,設(shè)定5種開口方案(表1),開口寬度分別為施工通道總長的1/10、1/3、1/2、2/3以及將施工通道全部開挖,即開口寬度分別為300 m、1 000 m、1 500 m、2 000 m及3 000 m。

    表1 施工通道開口方案Tab.1 Schemes of opening-size of construction access

    2.2 施工通道開口寬度對水體交換情況的影響

    采用濃度對流擴散模型研究工程水域的水體交換情況。在潮流模型的基礎(chǔ)上,添加對流擴散模塊,求解對流擴散方程得到水體中指標物質(zhì)的濃度分布,計算時不考慮污染物降解。

    在工程水域人工島—施工通道—海岸之間形成的半封閉水域網(wǎng)格點上均布濃度為1 mg/L的保守污染物,邊界污染物初始濃度設(shè)定為0,計算污染物經(jīng)過30 d輸運擴散后的濃度分布,并采用物理自凈率和水體半更換期作為水體交換情況的量化指標。

    (1)物理自凈率。設(shè)經(jīng)過一段時間i后灣內(nèi)的污染物平均濃度為n%,則物理自凈率定義為

    圖5給出了施工通道開口方案四條件下,污染物經(jīng)歷不同時長的輸移擴散后在全灣水域的數(shù)值計算結(jié)果示例。其中圖5-a~圖5-f分別為污染物輸移擴散5 d、10 d、15 d、20 d、25 d以及30 d后全灣的濃度分布。

    圖6匯總給出了不同施工通道開口方案、不同時間點上的物理自凈率之比較。

    圖5 某時刻污染物濃度分布示例Fig.5 Pollutant distribution in the bay at a certain time

    從圖6中可以看出,總體而言,5種施工通道開口寬度的物理自凈率差別不大。具體比較可以看出:施工通道開口寬度與灣內(nèi)污染物擴散率并不呈現(xiàn)簡單的單增或單減關(guān)系,方案一、方案四和方案五的物理自凈率較方案二和方案三略大。這主要是因為當(dāng)施工通道開口較小時,開口處的水流湍急,流速較大,在灣內(nèi)形成環(huán)流,從而有利于灣內(nèi)污染物的輸移;當(dāng)施工通道開口較大時,盡管此時開口處的流速相對較小,但由于開口寬度大,流量也隨之增大,增強了灣內(nèi)水體與灣外水體的交換。

    (2)水體半更換期。類比于化學(xué)元素的半衰期,把灣內(nèi)污染物總量降低到起始污染物總量一半所需的時間定義為水體的半更換期,表明此時灣內(nèi)水體已被更換一半[6]。根據(jù)此定義,分別統(tǒng)計各方案的半更換期,分別為 313 h、350 h、361 h、299 h和288 h。其中,方案五的半更換期最短,其次為方案四。

    圖6 不同方案之間物理自凈率之比較Fig.6 Comparison of physical self-purification rate among different schemes

    圖7 追蹤粒子釋放點Fig.7 Release position of tracer particles

    2.3 對局部污染物擴散規(guī)律的影響

    灣內(nèi)局部范圍的污染物擴散規(guī)律直接影響到部分海域的環(huán)境質(zhì)量,灣內(nèi)水體的交換情況僅反映海灣整體水交換情況,并不能反應(yīng)局部區(qū)域的水體質(zhì)量,因而有必要對局部污染物的擴散規(guī)律進行討論。

    施工通道掩護的近岸海域流速較小,若在該海域排放污染物,將對海域水環(huán)境造成較大不良影響,此處著重關(guān)注污染物近岸排放的情況。為追蹤污染物擴散輸移規(guī)律,在施工通道掩護的水域內(nèi)釋放5個追蹤粒子,考慮污染物近岸排放,追蹤粒子均在靠近海岸500 m范圍內(nèi)釋放,粒子釋放地點如圖7所示。根據(jù)前人研究結(jié)果[7],在漲急時刻投放粒子更為合理,故本模型釋放追蹤粒子的時刻為2009年11月10日3:00的漲急時刻。

    限于篇幅,圖8中以各方案中3#粒子運動軌跡作為示例。比較圖8中各圖可以看出:對應(yīng)開口寬度較小的方案一和方案二,施工通道掩護的近岸水域流速極小,出現(xiàn)近乎死水區(qū),該局部污染物輸移非常緩慢,長時間不能排出;對應(yīng)方案三,落潮時可以將污染物較順利地排出灣外,但漲潮時將已經(jīng)排出灣外的污染物又重新帶回灣內(nèi),易造成局部潮流通道水域反復(fù)污染;對應(yīng)方案四,落潮時排出到灣外的污染物漲潮時不會再帶回灣內(nèi);對應(yīng)方案五,由于其開口最大,流速最小,落潮流速幾乎不足以將污染物帶出灣外。

    綜上所述,開口寬度為2/3 L時,更有利于局部污染物的排放,即從施工通道掩護的近岸局部水域污染物輸移情況看,方案四最佳。

    圖8 施工通道掩護的近岸水域污染物擴散軌跡Fig.8 Diffusion trajectory of pollutant in offshore area shielded by construction access

    表2匯總給出了各方案中示蹤粒子第一次到達灣外所經(jīng)歷的時長。由表2中可見,除方案三中的4#粒子第一次到達灣外時間較長外,從總體上而言,方案三和方案四中各追蹤粒子第一次到達灣外的時間明顯短于其他各方案,說明方案三和方案四施工通道掩護的近岸水域內(nèi)污染物擴散輸移速度較其他方案快。對于單個方案來講,這一統(tǒng)計結(jié)果和追蹤粒子釋放的位置之間存在一定的關(guān)系,粒子釋放的位置不同,統(tǒng)計結(jié)果也存在一定差異,但從整體規(guī)律上來講,方案一和方案二施工通道較長,其掩護的近岸海域中靠近海岸的區(qū)域流速較其他方案小,再加上施工通道開口寬度較小,阻礙了追蹤粒子向灣外的運動;方案五中由于施工通道全部開挖,開口寬度較大,其中一些粒子借助落潮流能很快到達灣外,但工程所在海域漲潮流流速大于落潮流流速,因此方案五中追蹤粒子很快返回灣內(nèi),且在計算時間內(nèi)追蹤粒子長時間停留在灣內(nèi)。

    2.4 方案推薦

    在不同的評價指標下得到的最優(yōu)方案并不完全一致,若以物理自凈能力為評價指標,方案一、四、五的計算結(jié)果較為突出;若以半更換期為評價指標,方案四和方案五較其他方案好;從局部區(qū)域污染物擴散軌跡來看,方案四較其他方案理想;根據(jù)局部區(qū)域追蹤粒子第一次到達灣外的歷時來看,方案三和方案四較優(yōu)。綜合以上各評價指標,方案四是較為理想的方案,因而推薦方案四(即開挖2/3 L)為最終方案。

    表2 各方案追蹤粒子第一次到達灣外經(jīng)歷時長Tab.2 The time when tracer particles in local area first get out of the bay of each scheme

    3 結(jié)語

    本文以大連海上機場工程為背景,基于MIKE21 flow model FM模塊建立了金州灣二維潮流模型和二維污染物輸移擴散模型,計算結(jié)果與實測值的對比表明模型是可靠的。在此基礎(chǔ)上計算機場施工通道不同開口寬度下的物理自凈能力、水體半更換期、局部污染物擴散軌跡、局部污染物第一次到達灣外所經(jīng)歷的時長,并以此為指標確定施工通道開挖全長的2/3(即2 000 m)時為最佳方案。本文計算結(jié)果可為大連海上機場施工通道布置的進一步改進提供可靠的依據(jù)。

    [1]王諾,顏華錕,左書華,等.大連海上機場人工島建設(shè)對區(qū)域水動力及海床沖淤影響分析[J].水運工程,2012(4):5-11.WANG N,YAN H K,ZUO S H,et al.Impact of Exploitation of Dalian Offshore Airport Artificial Island on Regional Hydrodynamic Condition and Erosion-deposition Variation of Jinzhou bay[J].Port&Waterway Engineering,2012(4):5-11.

    [2]王志力,陸永軍,耿艷芬.基于非結(jié)構(gòu)網(wǎng)格有限體積法的二維高精度物質(zhì)輸運模擬[J].水科學(xué)進展,2008,19(4):531-536.WANG Z L,LU Y J,GENG Y F.High Resolution Two-dimensional Numerical Model of Scalar Transport Based on Unstructured Grid Finite Volume Method[J].Advances in Water Science,2008,19(4):531-536.

    [3]賴錫軍,曲卓杰,周杰,等.非結(jié)構(gòu)網(wǎng)格上的三維淺水流動數(shù)值模型[J].水科學(xué)進展,2006,17(5):693-699.LAI X J,QU Z J,ZHOU J,et al.3-D Hydrodynamic Model for Shallow Water on Unstructured Grids[J].Advances in Water Science,2006,17(5):693-699.

    [4]WANG J H,SHEN Y M,GUO Y K.Seasonal circulation and influence factors of the Bohai Sea:a numerical study based on Lagrangian particle tracking method[J].Ocean Dynamics,2010(60):1 581-1 596.

    [5]許婷.丹麥 MIKE21 模型概述及應(yīng)用實例[J].水利科技與經(jīng)濟,2010,16(8):867-869.XU T.Calculation Principle and Application Example of a Two-dimensional Flow Model-MIKE 21 HD[J].Water Conservancy Science and Technology and Economy,2010,16(8):867-869.

    [6]胡建宇,傅子瑯.羅源灣潮流及海水半更換期的數(shù)值研究[J].廈門大學(xué)學(xué)報:自然科學(xué)版,1989,28(S1):24-39.HU J Y,FU Z L.Numerical Studies on Half-exchange Period and Tidal Current in Luoyuan Bay[J].Journal of Xiamen University:Natural science,1989,28(S1):24-39.

    [7]趙亮,魏皓,趙建中.膠州灣水交換數(shù)值研究[J].海洋與湖沼,2002,33(3):23-29.ZHAO L,WEI H,ZHAO J Z.Numerical Study on Water Exchange in Jiaozhou Bay [J].Oceanlogia et Limnologia Sinica,2002,33(3):23-29.

    猜你喜歡
    模型施工
    一半模型
    后澆帶施工技術(shù)在房建施工中的踐行探索
    后澆帶施工技術(shù)在房建施工中的應(yīng)用
    土木工程施工技術(shù)創(chuàng)新探討
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    3D打印中的模型分割與打包
    土木工程施工實習(xí)的探討與實踐
    扶貧村里施工忙
    河南電力(2016年5期)2016-02-06 02:11:34
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    精品亚洲成a人片在线观看 | 一区在线观看完整版| 婷婷色综合大香蕉| 搡老乐熟女国产| 色婷婷久久久亚洲欧美| 丝袜喷水一区| 精品国产露脸久久av麻豆| 免费观看a级毛片全部| 国产视频首页在线观看| 中文乱码字字幕精品一区二区三区| 国产深夜福利视频在线观看| .国产精品久久| 国产真实伦视频高清在线观看| 亚洲精品乱久久久久久| 欧美zozozo另类| 韩国高清视频一区二区三区| 日日摸夜夜添夜夜添av毛片| 精品人妻偷拍中文字幕| 欧美另类一区| 免费黄色在线免费观看| 中文字幕久久专区| 纵有疾风起免费观看全集完整版| 丰满人妻一区二区三区视频av| 91精品国产国语对白视频| 久久久久性生活片| 成人午夜精彩视频在线观看| 亚洲精品视频女| 久久久国产一区二区| 18禁裸乳无遮挡免费网站照片| 久久久久久久大尺度免费视频| 在线播放无遮挡| 欧美精品人与动牲交sv欧美| 一级a做视频免费观看| 国产亚洲精品久久久com| 九九爱精品视频在线观看| 精品少妇黑人巨大在线播放| 免费人妻精品一区二区三区视频| 成人特级av手机在线观看| 男人爽女人下面视频在线观看| 在线观看人妻少妇| 性高湖久久久久久久久免费观看| 亚洲怡红院男人天堂| 精品少妇久久久久久888优播| 久久精品夜色国产| 欧美老熟妇乱子伦牲交| 亚洲一区二区三区欧美精品| 丰满迷人的少妇在线观看| 少妇 在线观看| 成年女人在线观看亚洲视频| 身体一侧抽搐| 亚州av有码| 午夜福利影视在线免费观看| 精品久久久噜噜| 高清午夜精品一区二区三区| 精品久久久久久电影网| 一二三四中文在线观看免费高清| 国产深夜福利视频在线观看| 久久国产精品大桥未久av | 最近中文字幕高清免费大全6| 韩国av在线不卡| 少妇的逼水好多| 纵有疾风起免费观看全集完整版| 男女国产视频网站| 热re99久久精品国产66热6| kizo精华| 亚洲一区二区三区欧美精品| 午夜免费鲁丝| 亚洲美女搞黄在线观看| 99视频精品全部免费 在线| 国产伦精品一区二区三区视频9| 男女国产视频网站| 少妇人妻精品综合一区二区| 在线观看免费视频网站a站| 人妻夜夜爽99麻豆av| 中文精品一卡2卡3卡4更新| 美女cb高潮喷水在线观看| 涩涩av久久男人的天堂| 永久免费av网站大全| 亚州av有码| 热99国产精品久久久久久7| 亚洲av成人精品一区久久| 国产精品久久久久久久电影| 日韩av免费高清视频| 老司机影院成人| 精品久久久噜噜| 国产免费一级a男人的天堂| 在线观看三级黄色| 亚洲国产精品成人久久小说| 久久亚洲国产成人精品v| 久久久午夜欧美精品| 搡老乐熟女国产| 韩国高清视频一区二区三区| 亚洲中文av在线| 日韩精品有码人妻一区| 亚洲精品国产av成人精品| 国产大屁股一区二区在线视频| 久久毛片免费看一区二区三区| 国产毛片在线视频| 欧美人与善性xxx| 亚洲在久久综合| 国产在线男女| 在线观看免费视频网站a站| 久久久久久久久久人人人人人人| 日韩三级伦理在线观看| 在线免费观看不下载黄p国产| 赤兔流量卡办理| 这个男人来自地球电影免费观看 | 久久精品国产自在天天线| 亚洲欧美日韩卡通动漫| 亚洲精品国产av蜜桃| 看免费成人av毛片| 国产成人免费观看mmmm| 精品人妻熟女av久视频| 小蜜桃在线观看免费完整版高清| 午夜福利高清视频| 在线观看美女被高潮喷水网站| 欧美3d第一页| 国产成人freesex在线| 精品久久久久久久末码| 九九爱精品视频在线观看| av福利片在线观看| 又粗又硬又长又爽又黄的视频| 免费观看无遮挡的男女| 久久99热6这里只有精品| 亚洲人成网站在线观看播放| 欧美区成人在线视频| 看免费成人av毛片| 亚洲电影在线观看av| 国产精品国产三级专区第一集| 少妇高潮的动态图| 欧美三级亚洲精品| 日本猛色少妇xxxxx猛交久久| 男人爽女人下面视频在线观看| 99久久精品一区二区三区| 亚洲激情五月婷婷啪啪| 欧美zozozo另类| 极品教师在线视频| 国产高清不卡午夜福利| 精品久久国产蜜桃| 热re99久久精品国产66热6| 免费看光身美女| 黄片wwwwww| 国产成人a∨麻豆精品| 中文字幕人妻熟人妻熟丝袜美| 全区人妻精品视频| 亚洲欧洲日产国产| 特大巨黑吊av在线直播| 国产一区二区三区综合在线观看 | 国产欧美日韩一区二区三区在线 | 久久精品国产亚洲av涩爱| 日韩中字成人| 一区二区三区免费毛片| 蜜桃在线观看..| 3wmmmm亚洲av在线观看| 下体分泌物呈黄色| 看十八女毛片水多多多| a级毛片免费高清观看在线播放| 免费看日本二区| 成人亚洲欧美一区二区av| 在线观看三级黄色| 亚洲av不卡在线观看| 美女脱内裤让男人舔精品视频| 男女边吃奶边做爰视频| 国产精品蜜桃在线观看| 日本免费在线观看一区| 色视频在线一区二区三区| 极品少妇高潮喷水抽搐| 91久久精品电影网| 午夜福利视频精品| 18+在线观看网站| 精品国产乱码久久久久久小说| 婷婷色综合www| 蜜桃在线观看..| 亚洲怡红院男人天堂| 麻豆成人午夜福利视频| h日本视频在线播放| 午夜激情久久久久久久| 亚洲精品国产av蜜桃| 亚洲欧美一区二区三区黑人 | a级毛色黄片| 黄色日韩在线| 欧美性感艳星| 精品视频人人做人人爽| 一区二区av电影网| 国产精品99久久99久久久不卡 | 成人午夜精彩视频在线观看| 精品国产三级普通话版| 国产淫片久久久久久久久| 蜜桃在线观看..| 精品人妻视频免费看| 欧美国产精品一级二级三级 | 纵有疾风起免费观看全集完整版| 亚洲国产精品一区三区| 国产淫片久久久久久久久| 成人综合一区亚洲| 国产精品秋霞免费鲁丝片| 国产精品久久久久久精品电影小说 | 亚洲成人一二三区av| 国产精品不卡视频一区二区| www.色视频.com| 国产美女午夜福利| 亚洲在久久综合| 久久99蜜桃精品久久| 天美传媒精品一区二区| 一级av片app| 欧美日韩一区二区视频在线观看视频在线| 黑人高潮一二区| 欧美日韩国产mv在线观看视频 | 精品久久久久久电影网| 欧美精品亚洲一区二区| 欧美激情国产日韩精品一区| 26uuu在线亚洲综合色| 久久国内精品自在自线图片| 国产黄色视频一区二区在线观看| 在线 av 中文字幕| 97精品久久久久久久久久精品| 我要看黄色一级片免费的| 国产精品久久久久成人av| av国产久精品久网站免费入址| 婷婷色av中文字幕| 男女无遮挡免费网站观看| 网址你懂的国产日韩在线| 少妇人妻一区二区三区视频| 激情 狠狠 欧美| 18禁裸乳无遮挡免费网站照片| 一级黄片播放器| 纯流量卡能插随身wifi吗| h视频一区二区三区| 伦理电影大哥的女人| 自拍偷自拍亚洲精品老妇| 26uuu在线亚洲综合色| 欧美精品国产亚洲| 夫妻午夜视频| 卡戴珊不雅视频在线播放| 国产淫片久久久久久久久| 在线观看一区二区三区激情| h日本视频在线播放| 日本-黄色视频高清免费观看| 国产精品三级大全| 国产中年淑女户外野战色| .国产精品久久| 国产一区亚洲一区在线观看| 韩国高清视频一区二区三区| 久久女婷五月综合色啪小说| 欧美性感艳星| 蜜桃久久精品国产亚洲av| 最近中文字幕2019免费版| 精品午夜福利在线看| 下体分泌物呈黄色| 国产男女超爽视频在线观看| 伦理电影免费视频| 免费黄网站久久成人精品| 亚洲无线观看免费| 黑人猛操日本美女一级片| 成人二区视频| 波野结衣二区三区在线| 国产成人精品婷婷| av在线老鸭窝| 国产极品天堂在线| 国产免费一区二区三区四区乱码| 日韩电影二区| 99热网站在线观看| 少妇熟女欧美另类| 五月开心婷婷网| 十分钟在线观看高清视频www | 亚洲av免费高清在线观看| 欧美一级a爱片免费观看看| 男女下面进入的视频免费午夜| 一级毛片我不卡| 美女中出高潮动态图| 国产亚洲5aaaaa淫片| 高清不卡的av网站| 在线精品无人区一区二区三 | 国语对白做爰xxxⅹ性视频网站| 97超碰精品成人国产| 婷婷色综合大香蕉| 国产片特级美女逼逼视频| 中文资源天堂在线| 日韩大片免费观看网站| 狠狠精品人妻久久久久久综合| 国产有黄有色有爽视频| 国产视频首页在线观看| 亚洲不卡免费看| 国产精品成人在线| 高清欧美精品videossex| 秋霞伦理黄片| 亚洲av中文av极速乱| 毛片女人毛片| 精品少妇黑人巨大在线播放| 国产男人的电影天堂91| 欧美成人一区二区免费高清观看| 久久鲁丝午夜福利片| 色吧在线观看| 干丝袜人妻中文字幕| 成人亚洲精品一区在线观看 | 特大巨黑吊av在线直播| 肉色欧美久久久久久久蜜桃| 超碰av人人做人人爽久久| 国产亚洲一区二区精品| 日韩av不卡免费在线播放| 久久99热这里只频精品6学生| 高清在线视频一区二区三区| 免费高清在线观看视频在线观看| 久久精品夜色国产| 色网站视频免费| 大片电影免费在线观看免费| 久久精品国产鲁丝片午夜精品| 一区在线观看完整版| 中文字幕精品免费在线观看视频 | 最近中文字幕2019免费版| 国产精品秋霞免费鲁丝片| 国产成人免费无遮挡视频| 国产免费又黄又爽又色| 亚洲精品自拍成人| 18禁在线无遮挡免费观看视频| 中文字幕久久专区| 人人妻人人澡人人爽人人夜夜| 久久久久视频综合| 日本一二三区视频观看| av专区在线播放| 哪个播放器可以免费观看大片| 啦啦啦中文免费视频观看日本| 色视频www国产| 欧美高清成人免费视频www| 亚洲av成人精品一二三区| 日韩欧美精品免费久久| 久久精品国产亚洲av天美| 免费看光身美女| 亚洲天堂av无毛| 噜噜噜噜噜久久久久久91| 亚洲欧美日韩另类电影网站 | 国产在线视频一区二区| 日韩一本色道免费dvd| 亚洲精品国产色婷婷电影| 黄色视频在线播放观看不卡| 欧美一级a爱片免费观看看| 国产黄频视频在线观看| 免费大片18禁| 亚洲精品456在线播放app| kizo精华| 日韩人妻高清精品专区| 久久久国产一区二区| 一二三四中文在线观看免费高清| 韩国高清视频一区二区三区| 久久亚洲国产成人精品v| 亚洲真实伦在线观看| 丰满迷人的少妇在线观看| 少妇熟女欧美另类| 午夜视频国产福利| 国产有黄有色有爽视频| 日韩不卡一区二区三区视频在线| 国产成人免费无遮挡视频| 香蕉精品网在线| 国产精品人妻久久久久久| 午夜福利高清视频| 久久久久久伊人网av| 美女国产视频在线观看| 精品久久久久久久久av| 毛片一级片免费看久久久久| 97超碰精品成人国产| 一级毛片 在线播放| 国产中年淑女户外野战色| 国产熟女欧美一区二区| 美女cb高潮喷水在线观看| 亚洲久久久国产精品| 777米奇影视久久| 国产一级毛片在线| 高清黄色对白视频在线免费看 | 一级毛片久久久久久久久女| 国产精品不卡视频一区二区| 多毛熟女@视频| 亚洲精品一区蜜桃| 国产精品国产三级专区第一集| 免费在线观看成人毛片| 肉色欧美久久久久久久蜜桃| 美女内射精品一级片tv| 三级经典国产精品| 日韩强制内射视频| 国产 一区 欧美 日韩| 国产精品福利在线免费观看| 简卡轻食公司| 中文天堂在线官网| 久久国内精品自在自线图片| 看免费成人av毛片| 少妇裸体淫交视频免费看高清| 麻豆成人av视频| 亚洲精品色激情综合| 亚洲人成网站高清观看| 一区二区av电影网| tube8黄色片| 观看av在线不卡| 久久久久久久国产电影| 在线观看免费日韩欧美大片 | 久久99精品国语久久久| 色婷婷久久久亚洲欧美| 午夜福利高清视频| 老司机影院成人| 成人毛片a级毛片在线播放| 欧美xxxx黑人xx丫x性爽| 18+在线观看网站| 亚洲内射少妇av| 久久这里有精品视频免费| av又黄又爽大尺度在线免费看| 婷婷色av中文字幕| 一区二区av电影网| 在线亚洲精品国产二区图片欧美 | av专区在线播放| 国产亚洲午夜精品一区二区久久| 久久国产精品男人的天堂亚洲 | 久久久久久久大尺度免费视频| 国产爱豆传媒在线观看| 一级黄片播放器| 久久人妻熟女aⅴ| 免费大片18禁| 黄色配什么色好看| 欧美一区二区亚洲| 国产伦在线观看视频一区| 精品久久国产蜜桃| 欧美97在线视频| 欧美bdsm另类| 国产久久久一区二区三区| 欧美精品一区二区免费开放| 在线免费观看不下载黄p国产| 在线观看免费日韩欧美大片 | 亚洲国产欧美在线一区| 欧美极品一区二区三区四区| 亚洲国产精品专区欧美| 欧美成人精品欧美一级黄| 国产精品一区二区性色av| 国产亚洲午夜精品一区二区久久| 国产无遮挡羞羞视频在线观看| 18+在线观看网站| 高清黄色对白视频在线免费看 | 一本色道久久久久久精品综合| av不卡在线播放| 久久这里有精品视频免费| 久久99蜜桃精品久久| 18禁在线播放成人免费| 国产v大片淫在线免费观看| 三级经典国产精品| 九九久久精品国产亚洲av麻豆| 日本一二三区视频观看| 亚洲,欧美,日韩| 欧美精品国产亚洲| 建设人人有责人人尽责人人享有的 | 日韩强制内射视频| 国产爽快片一区二区三区| 亚洲精品自拍成人| 搡老乐熟女国产| 久久久久性生活片| 两个人的视频大全免费| a级毛片免费高清观看在线播放| 欧美精品亚洲一区二区| 日本黄色日本黄色录像| 精品少妇久久久久久888优播| 一级毛片aaaaaa免费看小| 国产高清国产精品国产三级 | 伦理电影免费视频| 久久这里有精品视频免费| 视频区图区小说| 一个人免费看片子| 一级毛片 在线播放| 精品亚洲乱码少妇综合久久| 久久热精品热| 视频中文字幕在线观看| 青春草视频在线免费观看| 午夜福利在线在线| 视频中文字幕在线观看| 国产日韩欧美亚洲二区| 国产爱豆传媒在线观看| 日韩成人av中文字幕在线观看| 99热这里只有精品一区| 午夜福利在线观看免费完整高清在| 国产亚洲欧美精品永久| 国产亚洲一区二区精品| 大香蕉久久网| av免费在线看不卡| 国产人妻一区二区三区在| 国产精品人妻久久久久久| 观看av在线不卡| 午夜福利视频精品| 亚洲国产毛片av蜜桃av| 日韩av在线免费看完整版不卡| 日本av手机在线免费观看| 一级爰片在线观看| 99久久中文字幕三级久久日本| 久久久久久久久久久丰满| 男女下面进入的视频免费午夜| 欧美成人午夜免费资源| 精品久久久久久电影网| 伊人久久国产一区二区| 免费播放大片免费观看视频在线观看| 黄色欧美视频在线观看| 国产亚洲欧美精品永久| 国产在线视频一区二区| 美女脱内裤让男人舔精品视频| 国产美女午夜福利| 国产精品无大码| 中文字幕人妻熟人妻熟丝袜美| 欧美少妇被猛烈插入视频| 日韩电影二区| 精品少妇久久久久久888优播| 在线精品无人区一区二区三 | 五月开心婷婷网| 女性被躁到高潮视频| 国产黄频视频在线观看| 精品亚洲成a人片在线观看 | 精品酒店卫生间| 国产成人a∨麻豆精品| 久久女婷五月综合色啪小说| 狂野欧美激情性xxxx在线观看| 国产91av在线免费观看| 男人舔奶头视频| 狂野欧美白嫩少妇大欣赏| 国产欧美亚洲国产| 嫩草影院新地址| 亚洲国产毛片av蜜桃av| 日本午夜av视频| 欧美 日韩 精品 国产| 人妻夜夜爽99麻豆av| 欧美性感艳星| 九九爱精品视频在线观看| 街头女战士在线观看网站| 久久韩国三级中文字幕| 另类亚洲欧美激情| 日本av免费视频播放| 亚洲va在线va天堂va国产| 人人妻人人爽人人添夜夜欢视频 | 国产淫语在线视频| 这个男人来自地球电影免费观看 | 久久影院123| 在线播放无遮挡| 成人无遮挡网站| 熟女电影av网| 麻豆精品久久久久久蜜桃| 国产精品偷伦视频观看了| 国产亚洲av片在线观看秒播厂| 18禁裸乳无遮挡动漫免费视频| 亚洲人成网站在线观看播放| 久久ye,这里只有精品| 国产毛片在线视频| 最近中文字幕高清免费大全6| 亚洲中文av在线| 国产美女午夜福利| 亚洲精品,欧美精品| 一个人免费看片子| 国产男女内射视频| 九草在线视频观看| 亚洲国产最新在线播放| 亚洲av中文av极速乱| 久热久热在线精品观看| 啦啦啦在线观看免费高清www| 国产在线男女| 欧美xxxx黑人xx丫x性爽| 精品国产一区二区三区久久久樱花 | 国产精品一二三区在线看| 黄片无遮挡物在线观看| 精品久久久噜噜| 激情五月婷婷亚洲| av福利片在线观看| 国产淫片久久久久久久久| 男女下面进入的视频免费午夜| 涩涩av久久男人的天堂| 青青草视频在线视频观看| 伊人久久精品亚洲午夜| 丝袜脚勾引网站| 亚洲美女视频黄频| av在线老鸭窝| 久久久久久久精品精品| 大陆偷拍与自拍| 欧美丝袜亚洲另类| 欧美日韩视频精品一区| 直男gayav资源| 自拍欧美九色日韩亚洲蝌蚪91 | av专区在线播放| 国产亚洲欧美精品永久| 国产成人一区二区在线| 99久国产av精品国产电影| 成人毛片a级毛片在线播放| 久久精品国产亚洲av天美| 黑人高潮一二区| www.av在线官网国产| 国语对白做爰xxxⅹ性视频网站| 亚洲美女黄色视频免费看| 丝瓜视频免费看黄片| 丰满迷人的少妇在线观看| 女人久久www免费人成看片| 免费黄色在线免费观看| 亚洲成人中文字幕在线播放| 久久久精品免费免费高清| 18禁在线无遮挡免费观看视频| 亚洲精品aⅴ在线观看| 水蜜桃什么品种好| 高清日韩中文字幕在线| 国产大屁股一区二区在线视频| 观看美女的网站| 人人妻人人爽人人添夜夜欢视频 | 国模一区二区三区四区视频| 午夜福利在线观看免费完整高清在| 亚洲国产成人一精品久久久| 日韩大片免费观看网站| 看免费成人av毛片| 欧美日韩国产mv在线观看视频 | 狂野欧美激情性bbbbbb| 久久热精品热| 国产精品久久久久久精品古装| 2018国产大陆天天弄谢| 日日摸夜夜添夜夜添av毛片| 成人漫画全彩无遮挡| 性高湖久久久久久久久免费观看| 日本欧美国产在线视频| 国产精品成人在线| 精品亚洲成a人片在线观看 | 国产成人一区二区在线| 一区二区三区四区激情视频|