• <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)換方法初步研究
    中国美白少妇内射xxxbb| 高清午夜精品一区二区三区| 欧美日韩在线观看h| 精品国产国语对白av| 赤兔流量卡办理| 三级国产精品片| 在线观看美女被高潮喷水网站| 久久久久久久亚洲中文字幕| 久久国内精品自在自线图片| 精品国产一区二区三区久久久樱花| 青春草视频在线免费观看| 青春草亚洲视频在线观看| 又大又黄又爽视频免费| 黄色欧美视频在线观看| 国产片特级美女逼逼视频| 久久综合国产亚洲精品| .国产精品久久| 久久精品国产亚洲网站| 国产成人免费无遮挡视频| 桃花免费在线播放| 欧美日韩在线观看h| 夜夜看夜夜爽夜夜摸| 国产淫语在线视频| 青春草国产在线视频| √禁漫天堂资源中文www| 国产免费一级a男人的天堂| 又黄又爽又刺激的免费视频.| 免费看日本二区| 亚洲国产最新在线播放| 国产伦精品一区二区三区视频9| 成年美女黄网站色视频大全免费 | 亚洲国产精品成人久久小说| 五月玫瑰六月丁香| 国产男女内射视频| 只有这里有精品99| 熟女电影av网| 欧美变态另类bdsm刘玥| 男人狂女人下面高潮的视频| 一区二区三区乱码不卡18| 少妇裸体淫交视频免费看高清| 亚洲国产精品专区欧美| 久久99蜜桃精品久久| 国产日韩欧美亚洲二区| 亚洲国产色片| 久久 成人 亚洲| 搡老乐熟女国产| 黑人高潮一二区| 日韩av免费高清视频| 天美传媒精品一区二区| 99九九在线精品视频 | 久久久久久久久久人人人人人人| av福利片在线观看| 一级毛片黄色毛片免费观看视频| 大香蕉久久网| 亚洲av二区三区四区| 国产熟女欧美一区二区| 国产欧美亚洲国产| 欧美高清成人免费视频www| 又爽又黄a免费视频| 麻豆成人午夜福利视频| 亚洲精品第二区| 日本欧美国产在线视频| 亚洲国产精品成人久久小说| 欧美3d第一页| 亚洲欧美一区二区三区国产| 一区二区三区免费毛片| 免费观看无遮挡的男女| 欧美 亚洲 国产 日韩一| 日日爽夜夜爽网站| 18禁在线播放成人免费| 国产亚洲91精品色在线| 中国三级夫妇交换| 国产精品人妻久久久影院| 久久久精品94久久精品| 成人国产av品久久久| 国产熟女欧美一区二区| 夜夜看夜夜爽夜夜摸| 欧美激情极品国产一区二区三区 | 久久久久视频综合| 少妇丰满av| 国产极品粉嫩免费观看在线 | 精品亚洲成a人片在线观看| 国产精品蜜桃在线观看| 国产精品久久久久久久久免| 欧美成人午夜免费资源| 少妇的逼水好多| 97超碰精品成人国产| 两个人免费观看高清视频 | 日韩三级伦理在线观看| 视频区图区小说| 婷婷色综合大香蕉| 3wmmmm亚洲av在线观看| 精华霜和精华液先用哪个| 深夜a级毛片| freevideosex欧美| 婷婷色麻豆天堂久久| 男的添女的下面高潮视频| 成年美女黄网站色视频大全免费 | 女人久久www免费人成看片| 亚洲在久久综合| 熟女av电影| 国产欧美日韩精品一区二区| 一级毛片aaaaaa免费看小| 成人综合一区亚洲| 国产精品麻豆人妻色哟哟久久| 久久久a久久爽久久v久久| 久久精品国产鲁丝片午夜精品| 内地一区二区视频在线| 18禁动态无遮挡网站| 亚洲国产成人一精品久久久| 欧美成人精品欧美一级黄| 国产精品久久久久久精品古装| 免费看av在线观看网站| 国产黄色视频一区二区在线观看| 日韩亚洲欧美综合| 免费观看的影片在线观看| 大片免费播放器 马上看| av国产精品久久久久影院| 久久人人爽av亚洲精品天堂| 国产视频首页在线观看| 久久影院123| 国精品久久久久久国模美| 男男h啪啪无遮挡| 午夜91福利影院| 春色校园在线视频观看| 美女内射精品一级片tv| 王馨瑶露胸无遮挡在线观看| 午夜老司机福利剧场| 国模一区二区三区四区视频| 一区二区三区乱码不卡18| 啦啦啦中文免费视频观看日本| 久久6这里有精品| 一区二区三区免费毛片| 久久久久网色| 午夜福利影视在线免费观看| 夜夜骑夜夜射夜夜干| 亚洲精品国产av成人精品| 一本久久精品| 日韩欧美精品免费久久| 我要看黄色一级片免费的| 亚洲国产精品成人久久小说| 大香蕉97超碰在线| 国产精品久久久久成人av| 毛片一级片免费看久久久久| 国产av国产精品国产| 国语对白做爰xxxⅹ性视频网站| 男男h啪啪无遮挡| av国产久精品久网站免费入址| 这个男人来自地球电影免费观看 | 亚洲高清免费不卡视频| 精品卡一卡二卡四卡免费| av国产久精品久网站免费入址| 王馨瑶露胸无遮挡在线观看| 高清欧美精品videossex| 国产精品嫩草影院av在线观看| 99热这里只有是精品在线观看| 香蕉精品网在线| 亚洲av成人精品一区久久| 99热这里只有精品一区| av女优亚洲男人天堂| 在线 av 中文字幕| 免费看av在线观看网站| 97精品久久久久久久久久精品| 人人澡人人妻人| 国产免费视频播放在线视频| 欧美日韩视频精品一区| 成年人免费黄色播放视频 | av网站免费在线观看视频| 伊人亚洲综合成人网| 国产午夜精品久久久久久一区二区三区| 一本—道久久a久久精品蜜桃钙片| 一本—道久久a久久精品蜜桃钙片| 啦啦啦在线观看免费高清www| 啦啦啦在线观看免费高清www| 夫妻午夜视频| 综合色丁香网| 美女内射精品一级片tv| 在线 av 中文字幕| 国产乱来视频区| 亚洲国产欧美日韩在线播放 | 日本与韩国留学比较| 亚洲av.av天堂| 高清不卡的av网站| 一级av片app| 国产精品99久久久久久久久| 高清午夜精品一区二区三区| 亚洲电影在线观看av| 如日韩欧美国产精品一区二区三区 | 男女免费视频国产| 久久久久久久久大av| 日日爽夜夜爽网站| 久久久久久久大尺度免费视频| 精品人妻熟女毛片av久久网站| 国产精品福利在线免费观看| 9色porny在线观看| 建设人人有责人人尽责人人享有的| 国产精品国产av在线观看| 美女大奶头黄色视频| 97在线视频观看| 成年人午夜在线观看视频| 青春草视频在线免费观看| 国产男女超爽视频在线观看| 亚洲美女视频黄频| 嫩草影院新地址| 女性被躁到高潮视频| 国产91av在线免费观看| 高清在线视频一区二区三区| 久久这里有精品视频免费| 国产成人午夜福利电影在线观看| 国产一区二区三区综合在线观看 | 一区二区三区免费毛片| 国产老妇伦熟女老妇高清| 午夜激情福利司机影院| 91精品国产九色| 久久久久久久亚洲中文字幕| 日日摸夜夜添夜夜添av毛片| 看十八女毛片水多多多| 边亲边吃奶的免费视频| 久久国产乱子免费精品| 女人精品久久久久毛片| 亚洲欧洲日产国产| 亚洲欧美日韩卡通动漫| 久久久久精品性色| 极品教师在线视频| 精品亚洲成国产av| 国产午夜精品一二区理论片| 国产伦在线观看视频一区| 久久女婷五月综合色啪小说| 国产免费福利视频在线观看| freevideosex欧美| 美女国产视频在线观看| 亚洲欧洲精品一区二区精品久久久 | 日韩av免费高清视频| 日本av手机在线免费观看| 亚洲国产最新在线播放| 免费黄色在线免费观看| 久久久久久人妻| 日韩一区二区视频免费看| 内射极品少妇av片p| 99视频精品全部免费 在线| 男人舔奶头视频| 精品国产一区二区久久| 男女国产视频网站| 一区二区三区免费毛片| 97精品久久久久久久久久精品| 美女中出高潮动态图| 国产成人免费无遮挡视频| 国产亚洲5aaaaa淫片| 99久久精品热视频| 国产亚洲最大av| 精品久久久精品久久久| 日本黄色片子视频| 青春草视频在线免费观看| 女人久久www免费人成看片| 亚洲精品亚洲一区二区| 亚洲一级一片aⅴ在线观看| 国产综合精华液| 少妇 在线观看| 免费人成在线观看视频色| 黑丝袜美女国产一区| 亚洲精品亚洲一区二区| 99久久中文字幕三级久久日本| 狂野欧美激情性xxxx在线观看| 国模一区二区三区四区视频| 天堂俺去俺来也www色官网| 极品少妇高潮喷水抽搐| kizo精华| 高清黄色对白视频在线免费看 | 国产中年淑女户外野战色| 少妇精品久久久久久久| 国产精品.久久久| 久久国产乱子免费精品| 丰满乱子伦码专区| 欧美日韩av久久| 国产色爽女视频免费观看| 欧美激情极品国产一区二区三区 | 欧美精品亚洲一区二区| 黑人巨大精品欧美一区二区蜜桃 | 丰满乱子伦码专区| 春色校园在线视频观看| 亚洲精品乱码久久久v下载方式| 在线播放无遮挡| 日本黄大片高清| www.色视频.com| 综合色丁香网| 亚洲精品成人av观看孕妇| 一个人免费看片子| 成人亚洲精品一区在线观看| 街头女战士在线观看网站| 国产亚洲午夜精品一区二区久久| 一级毛片我不卡| 亚洲精品乱码久久久v下载方式| 亚洲久久久国产精品| 亚洲精品视频女| 精品久久久久久久久av| 中文字幕精品免费在线观看视频 | 国产熟女午夜一区二区三区 | 久久久国产一区二区| 亚洲国产最新在线播放| 欧美日韩亚洲高清精品| 黄色配什么色好看| 国产精品偷伦视频观看了| a级毛片免费高清观看在线播放| 国产精品熟女久久久久浪| 免费黄网站久久成人精品| 99国产精品免费福利视频| 欧美亚洲 丝袜 人妻 在线| 又粗又硬又长又爽又黄的视频| 一区二区三区免费毛片| 如何舔出高潮| 久久亚洲国产成人精品v| 日韩一本色道免费dvd| 亚洲av福利一区| 最近手机中文字幕大全| 免费黄网站久久成人精品| 内射极品少妇av片p| 国产黄频视频在线观看| 免费av中文字幕在线| 国产精品久久久久久精品电影小说| av黄色大香蕉| 国产成人精品福利久久| 欧美日韩一区二区视频在线观看视频在线| 69精品国产乱码久久久| 青青草视频在线视频观看| 久久99热6这里只有精品| 午夜免费男女啪啪视频观看| 三级国产精品欧美在线观看| 久久久久久久久大av| 99久久综合免费| 亚洲精品国产av蜜桃| 国产成人一区二区在线| 亚洲av欧美aⅴ国产| 亚洲精品乱久久久久久| 99九九在线精品视频 | 一级,二级,三级黄色视频| 高清av免费在线| 女性被躁到高潮视频| 在线观看一区二区三区激情| 免费观看av网站的网址| 99久久中文字幕三级久久日本| 精品少妇内射三级| 纯流量卡能插随身wifi吗| 人妻一区二区av| 久久久午夜欧美精品| 亚洲欧洲日产国产| a 毛片基地| 亚洲丝袜综合中文字幕| 大话2 男鬼变身卡| 亚洲自偷自拍三级| 夜夜骑夜夜射夜夜干| 老女人水多毛片| 成年美女黄网站色视频大全免费 | 26uuu在线亚洲综合色| 久久99热这里只频精品6学生| 男人狂女人下面高潮的视频| 欧美老熟妇乱子伦牲交| freevideosex欧美| 黑人猛操日本美女一级片| 国产亚洲精品久久久com| 国产亚洲av片在线观看秒播厂| 成人亚洲精品一区在线观看| 99热国产这里只有精品6| 一区二区三区四区激情视频| 亚洲av综合色区一区| 国产91av在线免费观看| 国产欧美日韩精品一区二区| 伊人亚洲综合成人网| 男女国产视频网站| 久久久久国产精品人妻一区二区| 午夜福利,免费看| 色婷婷av一区二区三区视频| 久久韩国三级中文字幕| 免费黄网站久久成人精品| 伦理电影大哥的女人| 久久久久精品久久久久真实原创| 亚洲,一卡二卡三卡| 中文字幕亚洲精品专区| 五月开心婷婷网| 亚洲av成人精品一二三区| 成人无遮挡网站| 涩涩av久久男人的天堂| 亚洲欧美一区二区三区黑人 | 91精品一卡2卡3卡4卡| 国产精品久久久久久久电影| 一本—道久久a久久精品蜜桃钙片| 中文字幕精品免费在线观看视频 | 国内揄拍国产精品人妻在线| 国模一区二区三区四区视频| 日本91视频免费播放| 国产免费一级a男人的天堂| 午夜久久久在线观看| 日韩一区二区三区影片| 亚洲欧美日韩东京热| 国产欧美日韩综合在线一区二区 | 国产欧美日韩一区二区三区在线 | 国产在线免费精品| 日本wwww免费看| 最近手机中文字幕大全| 成人综合一区亚洲| 亚洲精品自拍成人| 一级毛片久久久久久久久女| 免费黄频网站在线观看国产| 大话2 男鬼变身卡| √禁漫天堂资源中文www| 亚洲欧美精品自产自拍| 一级毛片 在线播放| 五月天丁香电影| 亚洲人成网站在线观看播放| 午夜日本视频在线| 久久久久久久精品精品| 久久久久久久久久久丰满| 纵有疾风起免费观看全集完整版| 国产高清国产精品国产三级| 亚洲国产欧美日韩在线播放 | 777米奇影视久久| 日韩 亚洲 欧美在线| 亚洲无线观看免费| 日韩精品有码人妻一区| 女性被躁到高潮视频| 中文资源天堂在线| 九九久久精品国产亚洲av麻豆| 春色校园在线视频观看| 午夜免费男女啪啪视频观看| 亚洲精品日韩av片在线观看| 好男人视频免费观看在线| 免费人成在线观看视频色| 亚洲美女黄色视频免费看| 国产免费一区二区三区四区乱码| 国产日韩欧美在线精品| 一级,二级,三级黄色视频| 精品少妇久久久久久888优播| 亚洲国产精品国产精品| 成人毛片60女人毛片免费| 少妇人妻久久综合中文| 国模一区二区三区四区视频| 国产精品国产av在线观看| 69精品国产乱码久久久| 亚洲国产精品成人久久小说| 日韩,欧美,国产一区二区三区| 最近最新中文字幕免费大全7| 国产69精品久久久久777片| 国产熟女欧美一区二区| 狂野欧美激情性xxxx在线观看| 久久精品国产a三级三级三级| 亚洲激情五月婷婷啪啪| 91精品伊人久久大香线蕉| 国产成人午夜福利电影在线观看| 日韩不卡一区二区三区视频在线| 亚洲综合精品二区| 亚洲国产精品999| 七月丁香在线播放| 国产中年淑女户外野战色| 国产精品.久久久| 国产精品免费大片| 久久影院123| 老女人水多毛片| 亚洲天堂av无毛| 中文乱码字字幕精品一区二区三区| 伊人久久精品亚洲午夜| xxx大片免费视频| 晚上一个人看的免费电影| 制服丝袜香蕉在线| 一级二级三级毛片免费看| 新久久久久国产一级毛片| 只有这里有精品99| 国产在视频线精品| 99国产精品免费福利视频| 在线观看三级黄色| 国产精品蜜桃在线观看| 国内少妇人妻偷人精品xxx网站| 啦啦啦中文免费视频观看日本| 大又大粗又爽又黄少妇毛片口| av福利片在线| 曰老女人黄片| .国产精品久久| 国产综合精华液| 国产淫语在线视频| 成人国产麻豆网| 国产一级毛片在线| 婷婷色麻豆天堂久久| 下体分泌物呈黄色| 国产高清不卡午夜福利| 9色porny在线观看| 精品视频人人做人人爽| 人人妻人人澡人人看| 99热全是精品| 在线看a的网站| 亚洲婷婷狠狠爱综合网| 色5月婷婷丁香| 一边亲一边摸免费视频| 久久精品国产亚洲av天美| 亚洲一区二区三区欧美精品| 性色avwww在线观看| 国产伦在线观看视频一区| 国产亚洲91精品色在线| 高清毛片免费看| 亚洲欧美精品专区久久| 精品一区二区三区视频在线| 国产白丝娇喘喷水9色精品| 一级毛片电影观看| 久久久a久久爽久久v久久| 成人黄色视频免费在线看| 水蜜桃什么品种好| 大陆偷拍与自拍| 国产黄片视频在线免费观看| 日韩电影二区| 五月伊人婷婷丁香| 亚洲国产色片| 国产熟女欧美一区二区| 国产av国产精品国产| 国产伦精品一区二区三区四那| 精品一区二区免费观看| 丝袜在线中文字幕| 在线观看av片永久免费下载| 夜夜骑夜夜射夜夜干| 成人免费观看视频高清| 99九九线精品视频在线观看视频| 日韩视频在线欧美| 七月丁香在线播放| 精品99又大又爽又粗少妇毛片| 欧美日韩在线观看h| 免费黄色在线免费观看| 大话2 男鬼变身卡| 最近手机中文字幕大全| 亚洲经典国产精华液单| 91在线精品国自产拍蜜月| 亚洲欧美成人精品一区二区| 欧美国产精品一级二级三级 | 亚洲天堂av无毛| 国产 精品1| 亚洲精品aⅴ在线观看| 黄色怎么调成土黄色| 亚洲久久久国产精品| 亚洲精品国产色婷婷电影| 青春草视频在线免费观看| 亚洲一级一片aⅴ在线观看| 纯流量卡能插随身wifi吗| 只有这里有精品99| 熟女人妻精品中文字幕| 日本-黄色视频高清免费观看| 国产免费一区二区三区四区乱码| 99国产精品免费福利视频| 一级二级三级毛片免费看| 精品一品国产午夜福利视频| 久久久久视频综合| 丁香六月天网| 一级毛片 在线播放| 丝瓜视频免费看黄片| 69精品国产乱码久久久| 免费看av在线观看网站| 国产黄频视频在线观看| 97精品久久久久久久久久精品| 国产成人精品无人区| av不卡在线播放| 免费观看在线日韩| 一级毛片我不卡| 国产av码专区亚洲av| 99九九在线精品视频 | 亚洲情色 制服丝袜| 免费看不卡的av| 亚洲欧洲精品一区二区精品久久久 | 久久国产精品男人的天堂亚洲 | 99热这里只有是精品50| 一本—道久久a久久精品蜜桃钙片| 伊人久久国产一区二区| 国产伦精品一区二区三区视频9| 自线自在国产av| 久久久久精品久久久久真实原创| 极品教师在线视频| 国产伦在线观看视频一区| 国产亚洲最大av| 蜜桃在线观看..| 我要看日韩黄色一级片| 80岁老熟妇乱子伦牲交| 精品久久国产蜜桃| 在线观看免费视频网站a站| 九色成人免费人妻av| 青春草视频在线免费观看| 中文精品一卡2卡3卡4更新| 国产日韩欧美在线精品| xxx大片免费视频| 一级毛片aaaaaa免费看小| 91久久精品国产一区二区三区| 亚洲人成网站在线观看播放| 国产精品三级大全| 少妇被粗大的猛进出69影院 | 黄色视频在线播放观看不卡| 91久久精品电影网| 国产成人aa在线观看| 日韩av不卡免费在线播放| 亚洲精品成人av观看孕妇| av在线观看视频网站免费| 国产亚洲一区二区精品| 国产午夜精品久久久久久一区二区三区| 性色av一级| 欧美日韩视频精品一区| 亚洲国产欧美日韩在线播放 | 三级国产精品片| 国产精品无大码| 晚上一个人看的免费电影| 最新中文字幕久久久久| 国产深夜福利视频在线观看| 99久久精品热视频| 老女人水多毛片| 一边亲一边摸免费视频| h视频一区二区三区| 亚洲欧洲精品一区二区精品久久久 | 日韩电影二区| 精品一区在线观看国产| 国产精品国产av在线观看| av天堂久久9| 久久午夜福利片| h日本视频在线播放| 一本久久精品| 两个人的视频大全免费|