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

    港口工程影響下萊州灣西南側(cè)海域水動(dòng)力演化特征

    2016-07-11 08:54:44馮秀麗
    海洋科學(xué) 2016年3期
    關(guān)鍵詞:數(shù)值模擬

    劉 瀟,馮秀麗,劉 杰

    ?

    港口工程影響下萊州灣西南側(cè)海域水動(dòng)力演化特征

    劉瀟1,2,馮秀麗2,劉杰3

    (1.中國海洋大學(xué) 環(huán)境科學(xué)與工程學(xué)院,山東 青島 266100; 2.中國海洋大學(xué) 海地科學(xué)與探測技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,山東 青島 266100; 3.國家海洋局 第一海洋研究所,山東 青島 266061)

    摘要:河流輸沙、水深地形演變、大型海洋工程的建設(shè)均會(huì)對海洋水動(dòng)力產(chǎn)生明顯影響。本文基于MIKE21數(shù)值模型,對不同時(shí)期水深岸線條件下萊州灣西南側(cè)海域的水動(dòng)力演變特征及濰坊港建設(shè)對水動(dòng)力的影響進(jìn)行研究。結(jié)果表明,1984~2007年,黃河現(xiàn)行河口東南側(cè)高流速區(qū)位置向SE向移動(dòng),萊州灣灣口及西南側(cè)海域流速增大,增大幅度在10 cm/s左右。濰坊港10 km引堤工程的建設(shè)導(dǎo)致引堤端頭處流速增大10~20 cm/s,引堤兩側(cè)流速減小15~20 cm/s,流向偏轉(zhuǎn)量在10°~40°; 濰坊港擴(kuò)建工程的建設(shè)導(dǎo)致防波堤端頭處流速增大5~20 cm/s,兩側(cè)海域流速減小5~20 cm/s,流向偏轉(zhuǎn)5°~20°。黃河口沙嘴的形成對其南側(cè)海域形成有效掩護(hù),而濰坊港工程的建設(shè)對其兩側(cè)海域起到一定的掩護(hù)作用,導(dǎo)致該海域有效波高顯著減小,水動(dòng)力強(qiáng)度減弱。

    關(guān)鍵詞:海洋水動(dòng)力; 數(shù)值模擬; 港口建設(shè); 萊州灣南岸

    [Foundation: Marine Commonweal Scientific Research Foundation,No.201005009]

    河口海岸帶是人類活動(dòng)相對密集、經(jīng)濟(jì)相對發(fā)達(dá)的地區(qū)[1],據(jù)統(tǒng)計(jì),世界上67%的大中城市和60%的人口集中在沿海地區(qū)。河口海岸帶地區(qū)又是海陸交互作用相對集中的地帶[2],海岸地貌在波浪、潮汐、海流等相互作用下,各種過程相互耦合多變,海岸地形地貌及海洋動(dòng)力環(huán)境演變機(jī)制復(fù)雜。近年來日益增加的河口海岸帶開發(fā)利用活動(dòng)對海洋水動(dòng)力環(huán)境的改變具有重要作用[3],如圍填海工程、海洋堤壩、海上采砂等。因此,研究人類活動(dòng)對河口海岸沉積動(dòng)力環(huán)境的影響,尤其是分析海洋工程建設(shè)對海洋水動(dòng)力環(huán)境的影響成為沿海國家和地區(qū)的重要課題。

    本文基于收集的實(shí)測潮流數(shù)據(jù)和不同時(shí)期的水深岸線資料,利用目前已廣泛應(yīng)用[4-6]的 MIKE21數(shù)值模型,對研究區(qū)不同時(shí)期的潮流和波浪特征進(jìn)行模擬,研究港口工程影響下萊州灣西南側(cè)海域的海洋水動(dòng)力環(huán)境演變特征及濰坊港建設(shè)對其的影響。對保護(hù)周邊海洋環(huán)境、合理開發(fā)和利用海洋資源具有一定的理論和實(shí)際意義,為今后該類型海域海洋工程的選址和海洋環(huán)境保護(hù)等提供科學(xué)依據(jù)。

    1 研究區(qū)概況

    1.1自然環(huán)境概況

    研究區(qū)位于萊州灣西南側(cè)海域,受淄脈溝、小清河、白浪河、濰河等入海河流反復(fù)淤積、堆砌的影響[7-8],該海域形成了典型的粉砂質(zhì)海岸,海底地勢平坦,自西南向東北逐漸降低,起伏較小,水深等值線走向與岸線基本平行。

    研究區(qū)屬暖溫帶半濕潤季風(fēng)性氣候[9],多年平均風(fēng)速為3.1 m/s,主導(dǎo)風(fēng)向?yàn)镾SE向,出現(xiàn)頻率為10%,強(qiáng)風(fēng)向?yàn)镹E向。潮汐類型以不正規(guī)混合半日潮為主[10],漲潮流主流向?yàn)镾W向,落潮流主流向?yàn)镹E向。波浪以風(fēng)浪為主,常浪向?yàn)镹向,出現(xiàn)頻率為21.22%,次常浪向?yàn)镹NE向,出現(xiàn)頻率為16.14%,強(qiáng)浪向?yàn)镹NE向[11-12]。潮流屬于正規(guī)半日潮流,以往復(fù)流為主,落潮和漲潮平均流速相差不大,余流較?。ㄒ话阈∮?.5 cm/s)[13]。

    1.2濰坊港概況

    濰坊港位于萊州灣南側(cè)近岸海域,港區(qū)自 1996 年5月開工建設(shè),至1997年8月,10 km引堤基本建成,并在引堤端頭西側(cè)建成 5 000 t級碼頭。1998~2007年,濰坊港岸線形態(tài)基本沒有發(fā)生變化,2007 年 4月開始,濰坊港周邊海域人工構(gòu)筑物明顯增多,在10 km引堤的基礎(chǔ)上,建成了雙堤環(huán)抱、單一口門的港口布局(圖1)。目前濰坊港已順利完成了沿海防護(hù)堤港區(qū)工程段、兩道防波擋沙堤、萬噸級航道、西疏港公路路堤及萬噸級碼頭工程的建設(shè),成為萊州灣南岸最大的港口工程。

    圖1 研究區(qū)位置及濰坊港概況Fig.1 Location of study area and general situation of Weifang port

    1.3實(shí)測潮流特征

    本文收集了W1和W2站2007年10月26日~27日、B1~B4站2007年11月10日~11日的25 h連續(xù)實(shí)測潮流數(shù)據(jù),實(shí)測潮流觀測點(diǎn)位置及潮流玫瑰圖見圖2。

    從圖 2中可以看出,各實(shí)測潮流站位均表現(xiàn)出較強(qiáng)的往復(fù)流特性,潮流流向受岸線走向影響較大。其中,萊州灣西側(cè) B1和 B2站的主流向?yàn)?WNW~ESE向,濰坊港北側(cè)的B3、B4及W1、W2站的主流向?yàn)閃SW~ENE向。漲潮流主流向均為向岸方向,落潮流主流向?yàn)殡x岸方向。

    根據(jù)實(shí)測潮流資料,實(shí)測最大流速為73.4 cm/s,出現(xiàn)在B1站,對應(yīng)流向?yàn)?21.1°。漲潮流流速與落潮流相差不大,漲潮流略大于落潮流。萊州灣西側(cè)的 B1、B2站大潮期潮周期內(nèi)的平均流速在 35.3~42.3 cm/s,濰坊港北側(cè)的B3、B4及W1、W2站大潮期潮周期內(nèi)的平均流速在22.6~35.6 cm/s。

    潮流運(yùn)動(dòng)形式主要表現(xiàn)為往復(fù)流形式的半日潮流。根據(jù)實(shí)測潮流數(shù)據(jù),各站位潮流的旋轉(zhuǎn)方向有所差別,位于萊州灣西側(cè)開闊海域的B1、B2站潮流的旋轉(zhuǎn)方向?yàn)轫槙r(shí)針; 位于濰坊港北側(cè)的 B3、B4及W1、W2站潮流的旋轉(zhuǎn)方向?yàn)槟鏁r(shí)針,其中小清河口附近的B3站表現(xiàn)出較強(qiáng)的旋轉(zhuǎn)流特性。

    圖2 研究區(qū)各站位實(shí)測潮流玫瑰圖Fig.2 Observed tidal current at different stations

    2 潮流數(shù)值模擬

    2.1模型參數(shù)及驗(yàn)證

    MIKE21 Flow Model FM采用可隨意控制疏密的非結(jié)構(gòu)三角網(wǎng)格剖分計(jì)算域,采用標(biāo)準(zhǔn) Galerkin有限元法進(jìn)行水平空間離散[14],能夠較好的適用于岸線狀況復(fù)雜的近海區(qū)、河口區(qū)等。陳雪峰等[15]、王萬戰(zhàn)等[16]、陳昊等[17]對模型的建立、控制方程、輸入?yún)?shù)等均做了詳細(xì)闡述。

    本文基于MIKE21 Flow Model FM數(shù)值模型研究靜風(fēng)條件、港口工程影響下萊州灣西南側(cè)海域潮流演變特征。不同時(shí)期數(shù)值模型的水深岸線數(shù)據(jù)分別采用萊州灣海圖(11840,1984年測量,1986年 12月出版,比例尺1︰15萬; 11840,萊州灣膠萊河口以西海域的岸線及水深資料為2002年測量,2005年10月出版,比例尺1︰15萬)、萊州灣濰坊港及附近海圖(31501,岸線及水深資料為2007、2008年測量,2008 年11月出版,比例尺1︰3萬)及實(shí)際水深測量資料。

    利用研究區(qū)大潮期實(shí)測海流資料對模型進(jìn)行驗(yàn)證,驗(yàn)證點(diǎn)位置見圖 2,以距離濰坊港較遠(yuǎn)的 B1站和較近的W2站為例,流速和流向驗(yàn)證結(jié)果如圖3所示。驗(yàn)證結(jié)果表明,數(shù)值模擬結(jié)果與實(shí)測數(shù)據(jù)吻合良好,本文建立的水動(dòng)力數(shù)值模型能夠較好地反映研究區(qū)海域潮流特征。

    圖3 B1、W2站流速、流向驗(yàn)證曲線Fig.3 Verification of current speed and direction at the B1 and W2 stations

    2.2研究區(qū)流場特征(1984年)

    1984年研究區(qū)海域大潮期潮流特征如圖4所示,漲潮中間時(shí)潮流由灣外流向?yàn)硟?nèi),灣口處流向?yàn)镾SW向,向?yàn)硟?nèi)運(yùn)動(dòng)過程中,流向逐漸變?yōu)镾W向。黃河自1976年改道清水溝流路入海,至1984年已在河口處形成向東突出的沙嘴,受水下三角洲影響,附近海域潮流流向發(fā)生偏轉(zhuǎn)并形成高流速區(qū),最大流速為0.84 m/s,流速隨著水深的增加逐漸減小。落潮中間時(shí),潮流由灣內(nèi)向?yàn)惩饬鲃?dòng),流向?yàn)?NNE~NE向,黃河口附近海域潮流繞過沙嘴向?yàn)惩饬鲃?dòng),高流速區(qū)的最大流速為0.96 m/s。

    2.3研究區(qū)流場特征(2002年)

    2002年萊州灣西南側(cè)海域大潮期潮流特征如圖5所示,黃河河口處已形成向東南方向突出的巨大沙嘴,受其影響附近海域潮流流向發(fā)生明顯偏轉(zhuǎn)并在沙嘴東南側(cè)形成高流速區(qū),最大流速約為 1.0 m/s。落潮中間時(shí),高流速區(qū)的最大流速為 1.2 m/s,其他區(qū)域流速普遍介于0.2~0.8 m/s。受濰坊港10 km引堤工程的影響,萊州灣南側(cè)近岸海域潮流特征發(fā)生變化,引堤端頭處流速增大 10~20 cm/s,引堤兩側(cè)(特別是東側(cè))附近海域流速減小15~20 cm/s。

    2.4研究區(qū)流場特征(2007年)

    2007年萊州灣西南側(cè)海域大潮期潮流特征如圖6所示。與 2002年相比,研究區(qū)潮流變化主要出現(xiàn)在濰坊港附近海域。受濰坊港擴(kuò)建工程的影響,環(huán)抱式防波堤堤頭處流速增大5~20 cm/s,港口工程兩側(cè)(特別是東側(cè))附近海域流速減小 5~20 cm/s; 另外,濰坊港附近海域流向發(fā)生明顯偏轉(zhuǎn),環(huán)抱式防波堤堤頭處流向順時(shí)針偏轉(zhuǎn)約 5°~20°,港口工程兩側(cè)流向普遍發(fā)生逆時(shí)針方向偏轉(zhuǎn),偏轉(zhuǎn)量介于5°~20°。

    圖4 1984年萊州灣西南側(cè)海域大潮期潮流特征Fig.4 Current characteristics of the study area during spring tide in 1984

    圖5 2002年萊州灣西南側(cè)海域大潮期潮流特征Fig.5 Current characteristics of the study area during spring tide in 2002

    2.5潮流演變特征分析

    2002年與1984年大潮期漲潮中間時(shí)的流速、流向變化對比如圖 7所示。流速變化圖中正值代表流速增大,負(fù)值代表流速減??; 流向變化圖中,正值代表順時(shí)針方向旋轉(zhuǎn),負(fù)值代表逆時(shí)針方向旋轉(zhuǎn)。與1984年相比,2002年研究區(qū)潮流演變特征主要表現(xiàn)在以下幾個(gè)方面:

    1)與1984年相比,2002年黃河入??诟浇纬删薮笊匙?,岸線向東南推進(jìn)約18 km。受此影響,黃河三角洲附近海域潮流流速、流向均發(fā)生了較大改變,流向發(fā)生逆時(shí)針方向偏轉(zhuǎn),沙嘴北側(cè)海域偏轉(zhuǎn)量約為 20°~40°,南側(cè)海域偏轉(zhuǎn)量在 0°~200°。受黃河入??谏匙斓挠绊?,高流速區(qū)的位置向SE方向發(fā)生移動(dòng),流速較1984年增加約20~40 cm/s,黃河三角洲附近海域流速減小約0~30 cm/s。

    2)黃河清水溝流路沙嘴的形成,使得萊州灣灣口寬度變窄。與1984年相比,2002年黃河口至屺姆島之間的距離減少近20 km,萊州灣灣口束流現(xiàn)象增強(qiáng),灣口一線海域流速增大約 10 cm/s左右; 同時(shí),萊州灣西南側(cè)海域流速普遍比1984年稍大。

    3)濰坊港附近海域潮流特征發(fā)生改變,濰坊港10 km引堤堤頭處流速增大約10~20 cm/s,兩側(cè)流速減小約0~15 cm/s,工程周邊附近海域潮流流向發(fā)生順時(shí)針偏轉(zhuǎn),偏轉(zhuǎn)量約為10°~40°。

    圖6 2007年萊州灣西南側(cè)海域大潮期潮流特征Fig.6 Current characteristics of the study area during spring tide in 2007

    圖7 2002年與1984年漲潮中間時(shí)流速(a,單位: m/s)、流向(b,單位: °)變化量分布圖Fig.7 Contrast of current velocity(a)in meters per second and direction(b)in angle degrees at flood level between 1984 and 2002

    與2002年相比,2007年萊州灣西南側(cè)海域潮流特征主要在濰坊港附近海域發(fā)生了變化。2007年與2002年大潮期漲潮中間時(shí)流速、流向變化對比如圖8所示。濰坊港工程對附近海域潮流特征的影響主要集中在工程區(qū)附近海域,與2002年相比,環(huán)抱式防波堤堤頭處流速增大5~20 cm/s,港口工程兩側(cè)(特別是東側(cè))附近海域流速減小5~20 cm/s; 環(huán)抱式防波堤堤頭處流向發(fā)生順時(shí)針偏轉(zhuǎn),偏轉(zhuǎn)量約為5°~ 20°,港口工程兩側(cè)流向發(fā)生逆時(shí)針偏轉(zhuǎn),偏轉(zhuǎn)量在5°~20°。

    綜上所述,濰坊港位于萊州灣西南側(cè)近岸海域,潮流流速較小,水動(dòng)力條件較弱,又由于漲落潮流流向與港口的延伸方向基本一致,因此濰坊港工程對水動(dòng)力的影響主要集中在工程區(qū)周邊海域,對其他海域影響較小。

    3 波浪數(shù)值模擬

    波浪是近岸海域較為活躍的沉積動(dòng)力因素,在向近岸傳播過程中受水深與地形的影響,出現(xiàn)反射、繞射、破碎等現(xiàn)象,波長、波向、波高等會(huì)發(fā)生改變,進(jìn)而造成海底泥沙的起動(dòng)搬運(yùn)及岸灘的沖淤演變。

    萊州灣波浪主要受季風(fēng)控制,以風(fēng)浪為主[18],本文運(yùn)用 MIKE21 SW模型對研究區(qū)風(fēng)生波浪進(jìn)行模擬,計(jì)算了不同岸線條件、不同風(fēng)向作用下研究區(qū)的波浪特征,定性分析港口建設(shè)、岸線演變對附近海域波浪的影響程度和范圍。

    由于研究區(qū)常風(fēng)向?yàn)镾SE向、強(qiáng)風(fēng)向?yàn)槠玁向,本文分別模擬了濰坊港建設(shè)前(1984年)和建成后(2007年)岸線條件下,SSE和NW向大風(fēng)(8級風(fēng))作用24 h后的波浪特征。成導(dǎo)致其南側(cè)海域有效波高顯著減小,而濰坊港工程的建設(shè)對防波堤東側(cè)海域有較好的掩護(hù)作用,有效波高明顯減小。

    圖8 2007年與2002年漲潮中間時(shí)流速(a,單位: m/s)、流向(b,單位: °)變化量分布圖Fig.8 Contrast of current velocity(a)in meters per second and direction(b)in angle degrees at flood stage between 2002 and 2007

    3.2SSE向大風(fēng)作用下波浪特征

    3.1NW向大風(fēng)作用下波浪特征

    NW 向大風(fēng)作用下研究區(qū)有效波高分布特征如圖9所示,最大有效波高為2.3 m,主要分布在研究區(qū)東北側(cè)開闊海域,有效波高隨著水深的變淺而逐漸減小,近岸海域有效波高普遍小于0.3 m。與1984年相比,2007年黃河三角洲沙嘴及水下三角洲的形

    SSE向大風(fēng)作用下研究區(qū)有效波高分布特征如圖10所示,最大有效波高為1.95 m,主要分布在萊州灣灣口附近海域,自灣口向?yàn)硟?nèi)有效波高逐漸減小,近岸海域有效波高普遍介于 0~0.3 m。與 1984年相比,2007年黃河三角洲沙嘴南側(cè)海域以及濰坊港工程防波堤西側(cè)海域有效波高顯著減小。

    圖9 NW向大風(fēng)作用下研究區(qū)有效波高分布圖Fig.9 Distribution of significant wave height under the impact of northwesterly wind(20 m/s)

    綜上所述,在NW向和SSE向大風(fēng)作用下,萊州灣近岸淺水海域有效波高普遍較小。黃河水下三角洲及河口沙嘴的形成對其南側(cè)海域有一定的掩護(hù)作用;另外,由于濰坊港的主體走向與本文模擬的風(fēng)向近乎垂直,濰坊港工程建設(shè)對這兩種風(fēng)向下波浪的影響程度和范圍較大,導(dǎo)致港口兩側(cè)海域有效波高明顯減小。

    4 結(jié)論

    與1984年相比,2002年黃河口附近海域水深岸線、流速、流向發(fā)生了較大改變,高流速區(qū)位置向東南方向移動(dòng),萊州灣灣口及萊州灣西南側(cè)海域流速增大約10 cm/s; 濰坊港10 km引堤端頭處流速增大10~20 cm/s,引堤兩側(cè)附近海域流速減小15~20 cm/s,流向偏轉(zhuǎn)量在10°~40°之間。

    圖10 SSE向大風(fēng)作用下研究區(qū)有效波高分布圖Fig.10 Distribution of significant wave height under the impact of south-southeasterly wind(20 m/s)

    2007年濰坊港擴(kuò)建工程對附近海域潮流的影響主要集中在工程附近海域,環(huán)抱式防波堤堤頭處流速增大5~20 cm/s,工程兩側(cè)附近海域流速減小5~20 cm/s,流向偏轉(zhuǎn)量在5°~20°之間。

    在NW向和SSE向大風(fēng)作用下,萊州灣近岸海域有效波高普遍較小,黃河口沙嘴的形成對其南側(cè)海域形成有效掩護(hù),濰坊港工程則對其兩側(cè)海域起到一定的掩護(hù)作用,使有效波高減小,水動(dòng)力強(qiáng)度減弱。

    參考文獻(xiàn):

    [1]劉玉海.21世紀(jì)中國沿海城市面臨的主要地質(zhì)災(zāi)害

    問題[J].中國地質(zhì)災(zāi)害與防治學(xué)報(bào),1998,9(Supp):120-126.

    Liu Yuhai.Main geological hazards confronting the coastal cities of China in the twenty first century [J].The Chinese Journal of Geological Hazard and Control,1998,9(Supp): 120-126.

    [2]伊飛,張訓(xùn)華,胡克.海岸帶陸海相互作用研究綜

    述[J].海洋地質(zhì)前沿,2011,27(3): 28-34.

    Yi Fei,Zhang Xunhua,HU Ke.A review of researches on the land-ocean interaction in the coastal zone [J].Marine Geology Frontiers,2011,27(3): 28-34.

    [3]Fanos A M,Khafagy A A,Dean R G.Impacts of protective works on the Nile delta coast [J].Journal of Coastal Research,1995,11(2): 516-528.

    [4]馮秀麗,隋倩倩,林霖,等.威海靖海灣港區(qū)張家埠

    新港建設(shè)對泥沙沖淤影響預(yù)測分析[J].海洋科學(xué),2011,35(3): 72-76.

    Feng Xiuli,Sui Qianqian,LIN Lin,et al.Prediction of sediment souring and depositing driven by Zhangjiabu new port construction at Jinghai bay at Weihai [J].Marine Sciences,2011,35(3): 72-76.

    [5]劉瀟,馮秀麗,劉杰,等.威海張家埠新港建設(shè)對附近海域潮流影響預(yù)測分析[J].海岸工程,2013,32(1):10-18.Liu Xiao,F(xiàn)eng Xiuli,LIU Jie,et al.Prediction of impact of new port construction on tidal current in Zhangjiabu and its adjacent waters [J].Coastal Engineering,2013,32(1): 10-18.

    [6]周廣鎮(zhèn),馮秀麗,劉杰,等.萊州灣東岸近岸海域規(guī)劃圍填海后沖淤演變預(yù)測[J].海洋科學(xué),2014,38(1):15-19.Zhou Guangzhen,F(xiàn)eng Xiuli,LIU Jie,et al.Prediction of erosion evolution and deposition in the east coast of the Laizhou bay after the implemention of the coastal planning [J].Marine Sciences,2014,38(1):15-19.

    [7]豐愛平,夏東興,谷東起,等.萊州灣南岸海岸侵蝕過程與原因研究[J].海洋科學(xué)進(jìn)展,2006,24(1):83-90.Feng Aiping,Xia Dongxing,GU Dongqi,et al.Study on process and cause of the coastal erosion along the south coast of the Laizhou bay [J].Advances in Marine Science,2006,24(1): 83-90.

    [8]馬紹賽,辛福言,崔毅,等.黃河和小清河主要污染物入海量的估算[J].海洋水產(chǎn)研究,2004,25(5):47-51.Ma Shaosai,Xin Fuyan,CUI Yi,et al.Assessment of main pollution matter volume into the sea from Yellow River and Xiaoqing River [J].Marine Fisheries Research,2004,25(5): 47-51.

    [9]Luo X X,Zhang S S,Yang J Q,et al.Macrobenthic community in the Xiaoqing river estuary in Laizhou Bay,China [J].Journal of Ocean University of China,2013,12(3): 366-372.

    [10]程義吉,高菁.萊州灣海域水文特征及沖淤變化分析[J].海岸工程,2006,25(3): 1-6.Cheng Yiji,Gao Jing.Analysis of hydrographic charac-teristics and changes in scour and silting in the Laizhou bay area [J].Coastal Engineering,2006,25(3): 1-6.

    [11]齊惠林.濰坊港導(dǎo)堤工程實(shí)施后的地形變化分析[J].水道港口,2000,1: 36-40.Qi Huilin.Analysis on topographic change after setting up of Weifang harbor [J].Journal of Waterway and Harbor,2000,1: 36-40.

    [12]Chen G Q,Yi L,Chen S L,et al.Partitioning of grain-size components of estuarine sediments and implications for sediment transport in southwestern Laizhou Bay,China [J].Chinese Journal of Oceanology and Limnology,2013,31(4): 895-906.

    [13]劉國亭,閆新興,侯志強(qiáng).濰坊港沉積物特征及岸灘穩(wěn)定性研究[J].水道港口,2012,33(5): 380-386.Liu Guoting,Yan Xinxing,HOU Zhiqiang.Study on sedimentary characteristics and shore stability in Weifang port [J].Journal of Waterway and Harbor,2012,33(5): 380-386.

    [14]Janssen P,Lionello P,Zambresky L.On the interaction of wind and waves [J].Philosophical transactions of the royal society of London.Series A,Mathematical and Physical Sciences,1989,329: 289-301.

    [15]陳雪峰,王桂萱.MIKE21計(jì)算軟件及其在長興島海域改造工程上的應(yīng)用[J].大連大學(xué)學(xué)報(bào),2007,28(8):93-98.CHEN Xuefeng,WANG Guixuan.MIKE21 software and its application on the offshore reconstruction engineering of Changxing islands [J].Journal of Dalian university,2007,28(8): 93-98.

    [16]王萬戰(zhàn),董利瑾.渤海流場基本特征的MIKE21模擬研究[J].人民黃河,2007,29(10): 32-33.Wang Wanzhan,Dong Lijin.MIKE21 simulation of the basic characteristics of the flow field in Bohai [J].Yellow River,2007,29(10): 32-33.

    [17]陳昊,陳輔利,朱永英.基于MIKE 21 FM的普灣潮流場數(shù)值模擬[J].人民長江,2014,45(Supp): 132-134.Chen Hao,Chen Fuli,Zhu Yongying.Numerical simulation of tidal current based on MIKE 21 FM in Pu bay [J].Yangtze River,2014,45(Supp): 132-134.

    [18]盧曉東,劉艷霞,嚴(yán)立文.萊州灣西岸岸灘沖淤特征分析[J].海洋科學(xué),2008,32(10): 39-44.Lu Xiaodong,Liu Yanxia,Yan Liwen.Analysis of character of beach erosion and accumulation in the west coast of Laizhou bay [J].Marine Sciences,2008,32(10): 39-44.

    Received: Oct.9,2014
    Key words: hydrodynamic environment; numerical simulation; harbor construction; southern Laizhou Bay
    Abstract: Large-scale marine construction has a significant effect on the hydrodynamic environment in adjacent waters.In this paper,we report the results of our study of the hydrodynamic evolution of Laizhou Bay and the impact of Weifang port construction on the hydrodynamic environment.We used the MIKE 21 two-dimensional,depth-integrated hydrodynamic model.We found that,from 1984 to 2007,the position of a high-velocity zone located in the southeast side of the modern Yellow River estuary moved to the southeast,while the flow velocity of the mouth and the southwest side of Laizhou Bay increased by about 10 cm/s.Due to construction of the 10-km embankment of Weifang port,the flow velocity increased by 10?20 cm/s at the end of the causeway and decreased by 15?20 cm/s on both sides of the causeway.The flow velocity increased by 5?20 cm/s at the end of the breakwater and decreased by 5?20 cm/s on both sides of the port project due to the Weifang port expansion project.The formation of the Yellow River spit played a shielding role in the south-side water,and the construction of Weifang port played a shielding role in the water on both sides,leading to reduction of significant wave height.

    (本文編輯: 李曉燕)

    Hydrodynamic evolution characteristics of southwest Laizhou Bay under the effect of port construction

    LIU Xiao1,2,F(xiàn)ENG Xiu-li2,LIU Jie3
    (1.College of Environmental Science and Engineering,Ocean University of China,Qingdao 266100,China;2.Key Lab of Submarine Geosciences and Prospecting Techniques,Ocean University of China,Qingdao 266100,China; 3.First Institute of Oceanography,State Oceanic Administration,Qingdao 266061,China)

    中圖分類號:P736.2; P75

    文獻(xiàn)標(biāo)識碼:A

    文章編號:1000-3096(2016)03-0138-08

    doi:10.11759/hykx20141009001

    收稿日期:2014-10-09; 修回日期: 2014-12-10

    基金項(xiàng)目:海洋公益科研專項(xiàng)(201005009)

    作者簡介:劉瀟(1985-),女,內(nèi)蒙古包頭人,博士后,主要從事海洋沉積與工程環(huán)境方面研究,E-mail: liuxiao_0000@sina.com; 馮秀麗,通信作者,教授,主要從事海洋沉積與工程環(huán)境方面研究,E-mail:fengxiuli@ouc.edu.cn

    猜你喜歡
    數(shù)值模擬
    基于AMI的雙色注射成型模擬分析
    錐齒輪精密冷擺輾成形在“材料成型數(shù)值模擬”課程教學(xué)中的應(yīng)用
    基于氣象信息及風(fēng)場信息的風(fēng)機(jī)輪轂處風(fēng)速預(yù)測
    鉆孔灌注樁樁底沉渣對樁體承載特性影響的模擬分析
    西南地區(qū)氣象資料測試、預(yù)處理和加工研究報(bào)告
    科技資訊(2016年18期)2016-11-15 08:01:18
    張家灣煤礦巷道無支護(hù)條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開切眼錨桿支護(hù)參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機(jī)翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機(jī)的流場數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    精品国产三级普通话版| 97在线视频观看| 亚洲三级黄色毛片| 免费看日本二区| 免费大片18禁| 九九久久精品国产亚洲av麻豆| 人人妻人人看人人澡| 国产日本99.免费观看| 国产成人福利小说| 亚洲第一区二区三区不卡| 最近最新中文字幕大全电影3| 我要看日韩黄色一级片| 欧美日本亚洲视频在线播放| .国产精品久久| 成人一区二区视频在线观看| 中出人妻视频一区二区| 美女 人体艺术 gogo| 日本免费a在线| av免费观看日本| 日本黄色视频三级网站网址| 亚洲在线自拍视频| 亚洲精品国产av成人精品| 久久精品91蜜桃| 亚洲国产日韩欧美精品在线观看| 久久九九热精品免费| 日本爱情动作片www.在线观看| www日本黄色视频网| 欧美日韩一区二区视频在线观看视频在线 | 午夜精品在线福利| 国产成年人精品一区二区| 三级经典国产精品| 晚上一个人看的免费电影| videossex国产| 国产黄a三级三级三级人| 18禁黄网站禁片免费观看直播| 我的老师免费观看完整版| 成人午夜精彩视频在线观看| 熟女人妻精品中文字幕| 欧美日韩乱码在线| 波多野结衣高清作品| 亚洲精品久久久久久婷婷小说 | 成人无遮挡网站| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 一级av片app| 99久久精品国产国产毛片| 搡老妇女老女人老熟妇| 久久中文看片网| 久久99精品国语久久久| 婷婷精品国产亚洲av| 国产日韩欧美在线精品| 在线免费十八禁| 综合色丁香网| 美女国产视频在线观看| 22中文网久久字幕| 日韩,欧美,国产一区二区三区 | 午夜免费激情av| 天堂√8在线中文| 99九九线精品视频在线观看视频| 日韩欧美精品v在线| av在线蜜桃| 亚洲人成网站在线观看播放| 老司机福利观看| 成人鲁丝片一二三区免费| 美女内射精品一级片tv| 一个人免费在线观看电影| 亚洲精品日韩av片在线观看| 国产高清激情床上av| 97在线视频观看| 男女视频在线观看网站免费| 久久久精品94久久精品| 国产黄色小视频在线观看| 国产黄片视频在线免费观看| 成人特级av手机在线观看| 久久中文看片网| 久久亚洲精品不卡| 国产免费一级a男人的天堂| 午夜福利高清视频| 日韩欧美 国产精品| 日韩,欧美,国产一区二区三区 | 一本一本综合久久| 久久精品国产亚洲网站| 成年免费大片在线观看| 特大巨黑吊av在线直播| 国产精品永久免费网站| 老师上课跳d突然被开到最大视频| 乱码一卡2卡4卡精品| 亚洲图色成人| 日韩一区二区视频免费看| 99精品在免费线老司机午夜| 国产亚洲精品av在线| 国产黄片美女视频| 午夜福利高清视频| 91午夜精品亚洲一区二区三区| av在线蜜桃| 欧美日本亚洲视频在线播放| 在线免费观看的www视频| 高清日韩中文字幕在线| 国产精品国产三级国产av玫瑰| 国产一区二区在线av高清观看| 国产成人精品婷婷| 亚洲av第一区精品v没综合| 夜夜爽天天搞| 日韩 亚洲 欧美在线| 又爽又黄无遮挡网站| 亚洲图色成人| 日本三级黄在线观看| 女人被狂操c到高潮| а√天堂www在线а√下载| 国产爱豆传媒在线观看| а√天堂www在线а√下载| 成人三级黄色视频| 寂寞人妻少妇视频99o| 老女人水多毛片| 国产视频内射| 日本爱情动作片www.在线观看| 亚洲综合色惰| 久久精品人妻少妇| 天天一区二区日本电影三级| 亚洲天堂国产精品一区在线| 国产淫片久久久久久久久| 久久久精品94久久精品| 深夜a级毛片| 三级国产精品欧美在线观看| 精品不卡国产一区二区三区| 麻豆成人av视频| 黄色欧美视频在线观看| 男人狂女人下面高潮的视频| 高清毛片免费观看视频网站| 黄色视频,在线免费观看| 日韩制服骚丝袜av| 男的添女的下面高潮视频| 亚洲在久久综合| av在线老鸭窝| 91精品国产九色| 精品一区二区三区视频在线| 高清毛片免费看| 高清日韩中文字幕在线| 99国产极品粉嫩在线观看| 成人特级黄色片久久久久久久| 国产精品国产高清国产av| 在线国产一区二区在线| 国产高清三级在线| 色视频www国产| 国产高潮美女av| 久久久久久久午夜电影| 搡女人真爽免费视频火全软件| av国产免费在线观看| 亚洲精品自拍成人| 不卡一级毛片| 欧美不卡视频在线免费观看| 欧美日韩国产亚洲二区| 国产探花在线观看一区二区| 亚洲av成人精品一区久久| 老师上课跳d突然被开到最大视频| 久久精品国产鲁丝片午夜精品| 免费看av在线观看网站| 可以在线观看的亚洲视频| 久久久久九九精品影院| 最近中文字幕高清免费大全6| АⅤ资源中文在线天堂| 色吧在线观看| 草草在线视频免费看| av在线观看视频网站免费| 日本-黄色视频高清免费观看| 黄色配什么色好看| 成人综合一区亚洲| 精品久久久久久久久久久久久| 亚洲精品456在线播放app| 国产片特级美女逼逼视频| 日日摸夜夜添夜夜爱| 淫秽高清视频在线观看| 又粗又硬又长又爽又黄的视频 | 国产高潮美女av| 欧美精品一区二区大全| 久久久久久九九精品二区国产| .国产精品久久| 中文字幕av在线有码专区| 亚洲天堂国产精品一区在线| 欧美另类亚洲清纯唯美| 麻豆av噜噜一区二区三区| 麻豆国产av国片精品| 男插女下体视频免费在线播放| 女人十人毛片免费观看3o分钟| 久久久久久国产a免费观看| 又粗又硬又长又爽又黄的视频 | av在线天堂中文字幕| 男女边吃奶边做爰视频| 欧美精品国产亚洲| 黑人高潮一二区| 国产精品一区二区三区四区久久| 麻豆精品久久久久久蜜桃| 天天躁日日操中文字幕| 99国产精品一区二区蜜桃av| 少妇猛男粗大的猛烈进出视频 | 精品午夜福利在线看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 99热这里只有精品一区| 亚洲精品乱码久久久v下载方式| 精品久久久久久成人av| 嘟嘟电影网在线观看| 蜜桃亚洲精品一区二区三区| 亚洲成a人片在线一区二区| 99热6这里只有精品| av女优亚洲男人天堂| 久久久久九九精品影院| 18禁在线无遮挡免费观看视频| 99九九线精品视频在线观看视频| 亚洲无线在线观看| 麻豆国产97在线/欧美| 一级二级三级毛片免费看| 深爱激情五月婷婷| 99九九线精品视频在线观看视频| 美女xxoo啪啪120秒动态图| 校园人妻丝袜中文字幕| 亚洲av第一区精品v没综合| 插逼视频在线观看| 欧美激情久久久久久爽电影| 亚洲av二区三区四区| 一级黄色大片毛片| 在线观看66精品国产| 欧美区成人在线视频| 欧美日韩综合久久久久久| 亚洲欧美精品自产自拍| 三级国产精品欧美在线观看| 欧美精品一区二区大全| 不卡一级毛片| 精品一区二区三区人妻视频| 国产极品精品免费视频能看的| 亚洲国产高清在线一区二区三| av.在线天堂| 国产精品一区二区性色av| 亚洲成人久久性| 久久久久久九九精品二区国产| 12—13女人毛片做爰片一| 国内少妇人妻偷人精品xxx网站| 人妻夜夜爽99麻豆av| 亚洲经典国产精华液单| 伦理电影大哥的女人| 亚洲经典国产精华液单| 在线观看66精品国产| 成人特级黄色片久久久久久久| 国产成人精品久久久久久| 悠悠久久av| 此物有八面人人有两片| 深夜a级毛片| 亚洲av免费在线观看| 国产一区二区激情短视频| 在线观看美女被高潮喷水网站| 在线观看午夜福利视频| 午夜激情福利司机影院| 中文字幕免费在线视频6| 午夜久久久久精精品| 亚洲av中文字字幕乱码综合| 亚洲性久久影院| 免费无遮挡裸体视频| 三级男女做爰猛烈吃奶摸视频| 你懂的网址亚洲精品在线观看 | 色综合色国产| 国产一级毛片七仙女欲春2| 国产精品久久久久久亚洲av鲁大| 久久国产乱子免费精品| 久久久久久久久大av| 国产成人精品婷婷| 久久久精品94久久精品| 一级二级三级毛片免费看| 日日干狠狠操夜夜爽| 免费搜索国产男女视频| 日韩欧美精品免费久久| 国产片特级美女逼逼视频| 青春草国产在线视频 | 国产综合懂色| 成人欧美大片| 欧美激情国产日韩精品一区| h日本视频在线播放| 国产极品精品免费视频能看的| 久久99精品国语久久久| 国产精品av视频在线免费观看| 久久久久久久午夜电影| 欧美xxxx黑人xx丫x性爽| 久久精品夜色国产| 蜜桃久久精品国产亚洲av| 亚洲av中文字字幕乱码综合| 白带黄色成豆腐渣| 永久网站在线| 99国产精品一区二区蜜桃av| 精品午夜福利在线看| 日日干狠狠操夜夜爽| 成人av在线播放网站| 久久精品国产99精品国产亚洲性色| 亚洲国产精品国产精品| 国产精品人妻久久久久久| 欧美最新免费一区二区三区| 国产一区二区三区av在线 | 日本撒尿小便嘘嘘汇集6| 中出人妻视频一区二区| 久久久久久久亚洲中文字幕| 亚洲欧美成人精品一区二区| 婷婷色av中文字幕| 毛片女人毛片| 国产精品蜜桃在线观看 | 国产精品.久久久| 午夜福利在线在线| av又黄又爽大尺度在线免费看 | 一进一出抽搐gif免费好疼| 久久精品久久久久久噜噜老黄 | 久久久成人免费电影| 亚洲无线在线观看| videossex国产| 又爽又黄无遮挡网站| 免费一级毛片在线播放高清视频| 全区人妻精品视频| 26uuu在线亚洲综合色| or卡值多少钱| 国内精品宾馆在线| 欧美极品一区二区三区四区| 中文字幕久久专区| 亚洲在久久综合| 国产极品精品免费视频能看的| 97人妻精品一区二区三区麻豆| 五月玫瑰六月丁香| 欧美zozozo另类| 国产人妻一区二区三区在| 又爽又黄a免费视频| 久久久久免费精品人妻一区二区| 12—13女人毛片做爰片一| 日日撸夜夜添| 69人妻影院| 亚洲成a人片在线一区二区| 亚洲三级黄色毛片| 人人妻人人澡人人爽人人夜夜 | www.色视频.com| 亚洲七黄色美女视频| 久久精品夜色国产| 日本一本二区三区精品| 婷婷精品国产亚洲av| 日韩精品有码人妻一区| 美女xxoo啪啪120秒动态图| 51国产日韩欧美| 天堂网av新在线| 日本成人三级电影网站| 国产 一区精品| 久久久久久大精品| 99久国产av精品国产电影| 老司机福利观看| 91精品一卡2卡3卡4卡| 精品一区二区三区人妻视频| 久久亚洲精品不卡| 搡老妇女老女人老熟妇| a级毛色黄片| 免费av观看视频| 一区二区三区高清视频在线| 深夜精品福利| 亚洲精品影视一区二区三区av| 国产单亲对白刺激| 国产精品野战在线观看| 夜夜看夜夜爽夜夜摸| 男人舔女人下体高潮全视频| 最近中文字幕高清免费大全6| 桃色一区二区三区在线观看| 18+在线观看网站| 国内少妇人妻偷人精品xxx网站| 久久精品影院6| 精品欧美国产一区二区三| 久久久a久久爽久久v久久| 三级经典国产精品| 五月伊人婷婷丁香| 一级毛片久久久久久久久女| 可以在线观看的亚洲视频| 中文字幕人妻熟人妻熟丝袜美| 黄色一级大片看看| 亚洲欧美精品专区久久| 听说在线观看完整版免费高清| 久99久视频精品免费| 国产探花在线观看一区二区| 三级国产精品欧美在线观看| 亚洲国产欧洲综合997久久,| 国产老妇伦熟女老妇高清| 国产色婷婷99| 日本黄大片高清| 久久久久久久久久黄片| 欧美xxxx黑人xx丫x性爽| 亚洲美女搞黄在线观看| 欧美不卡视频在线免费观看| 欧美bdsm另类| 卡戴珊不雅视频在线播放| 欧洲精品卡2卡3卡4卡5卡区| 少妇高潮的动态图| 男女下面进入的视频免费午夜| 亚洲aⅴ乱码一区二区在线播放| 人妻久久中文字幕网| 国产精品1区2区在线观看.| avwww免费| 久久午夜福利片| 真实男女啪啪啪动态图| 日本黄大片高清| 国产精品三级大全| 乱系列少妇在线播放| 午夜精品一区二区三区免费看| 国产亚洲精品久久久久久毛片| 国产精品,欧美在线| 亚洲精品影视一区二区三区av| 又爽又黄无遮挡网站| 搞女人的毛片| 亚洲自拍偷在线| 免费观看人在逋| 校园春色视频在线观看| 亚洲成人av在线免费| 久久亚洲精品不卡| 夜夜爽天天搞| 日韩欧美国产在线观看| 久久精品91蜜桃| 中文精品一卡2卡3卡4更新| 国产精品一区www在线观看| 久久久欧美国产精品| 亚洲欧美精品专区久久| 一夜夜www| 久久久国产成人精品二区| 国产精品99久久久久久久久| 中文欧美无线码| 少妇猛男粗大的猛烈进出视频 | 能在线免费观看的黄片| 成人午夜精彩视频在线观看| 亚洲欧美精品专区久久| 免费人成视频x8x8入口观看| 老女人水多毛片| 日韩欧美 国产精品| 中文字幕熟女人妻在线| 成人鲁丝片一二三区免费| 日本熟妇午夜| 少妇丰满av| 午夜福利视频1000在线观看| 免费观看精品视频网站| 久久欧美精品欧美久久欧美| а√天堂www在线а√下载| 国产精品永久免费网站| 午夜福利在线观看免费完整高清在 | 我的老师免费观看完整版| 国产精品三级大全| 国产69精品久久久久777片| 久久久久久久久久久免费av| 99热只有精品国产| 非洲黑人性xxxx精品又粗又长| 国产免费男女视频| 成人国产麻豆网| 国产成人影院久久av| 久久久精品大字幕| eeuss影院久久| 久久久久网色| 啦啦啦韩国在线观看视频| 国产三级中文精品| 国产一区二区三区在线臀色熟女| 欧美3d第一页| 麻豆成人av视频| 丰满的人妻完整版| 久久久久网色| 91精品一卡2卡3卡4卡| 亚洲精品成人久久久久久| 国产精品久久视频播放| 精品国内亚洲2022精品成人| 国产三级在线视频| 日韩av在线大香蕉| 欧美一级a爱片免费观看看| 国产黄色小视频在线观看| 岛国在线免费视频观看| 久久人妻av系列| 亚洲av中文av极速乱| 99国产精品一区二区蜜桃av| 国产亚洲av片在线观看秒播厂 | 日韩av不卡免费在线播放| 免费大片18禁| av黄色大香蕉| 亚洲精品久久久久久婷婷小说 | 国产精品久久久久久亚洲av鲁大| 岛国在线免费视频观看| 久久久精品大字幕| 国产一区二区三区av在线 | av在线亚洲专区| 精品国产三级普通话版| 99在线视频只有这里精品首页| 国产成人一区二区在线| 国产成人91sexporn| 一级黄片播放器| 变态另类成人亚洲欧美熟女| 色综合站精品国产| 精品日产1卡2卡| 国产爱豆传媒在线观看| 女人十人毛片免费观看3o分钟| 日韩视频在线欧美| 91aial.com中文字幕在线观看| 男插女下体视频免费在线播放| 亚洲精华国产精华液的使用体验 | 非洲黑人性xxxx精品又粗又长| 婷婷精品国产亚洲av| 国国产精品蜜臀av免费| 亚洲久久久久久中文字幕| 国产乱人视频| 欧美成人a在线观看| 国产精品人妻久久久久久| 国产精品,欧美在线| 男人狂女人下面高潮的视频| ponron亚洲| 男插女下体视频免费在线播放| 亚洲欧美日韩东京热| 亚洲国产欧美人成| 老女人水多毛片| 最新中文字幕久久久久| 国产麻豆成人av免费视频| 亚洲第一区二区三区不卡| 久久午夜亚洲精品久久| 看免费成人av毛片| 久久久精品94久久精品| 在线播放国产精品三级| 国产精品久久久久久亚洲av鲁大| 欧美xxxx黑人xx丫x性爽| 免费av观看视频| av免费在线看不卡| 一级二级三级毛片免费看| 久久久久久久久久久丰满| 国产极品天堂在线| 99久久精品国产国产毛片| 亚洲无线观看免费| 深夜精品福利| 国产精品蜜桃在线观看 | 国产不卡一卡二| 狂野欧美白嫩少妇大欣赏| 欧美日韩综合久久久久久| 成人漫画全彩无遮挡| 男女那种视频在线观看| 哪个播放器可以免费观看大片| 午夜爱爱视频在线播放| 精品不卡国产一区二区三区| 国产黄片视频在线免费观看| 波多野结衣巨乳人妻| 国产私拍福利视频在线观看| 日本黄色片子视频| 边亲边吃奶的免费视频| 美女内射精品一级片tv| 日日摸夜夜添夜夜添av毛片| 国产av不卡久久| 亚洲国产精品国产精品| 日韩欧美精品v在线| 精品人妻熟女av久视频| 一级黄色大片毛片| 一级毛片久久久久久久久女| 成人午夜精彩视频在线观看| av免费在线看不卡| 午夜福利成人在线免费观看| 神马国产精品三级电影在线观看| 禁无遮挡网站| 国产精品一区二区三区四区免费观看| 久久久精品94久久精品| 联通29元200g的流量卡| 99热网站在线观看| 成熟少妇高潮喷水视频| 免费看a级黄色片| 国产人妻一区二区三区在| 在线观看免费视频日本深夜| 亚洲色图av天堂| 成人漫画全彩无遮挡| 成人一区二区视频在线观看| av免费在线看不卡| 蜜桃久久精品国产亚洲av| 少妇人妻一区二区三区视频| 午夜爱爱视频在线播放| 成年女人看的毛片在线观看| 免费搜索国产男女视频| 精品久久久久久久久av| 亚洲国产精品国产精品| 男人舔女人下体高潮全视频| 桃色一区二区三区在线观看| 久久精品人妻少妇| 精品一区二区三区视频在线| 熟女电影av网| 国产精品久久久久久av不卡| 国产探花极品一区二区| 亚洲最大成人av| 中文字幕av在线有码专区| 久久精品影院6| 国产高清三级在线| 久久精品国产清高在天天线| 欧美另类亚洲清纯唯美| 精品人妻视频免费看| 人妻系列 视频| 亚洲精品亚洲一区二区| av女优亚洲男人天堂| 1000部很黄的大片| 毛片女人毛片| 全区人妻精品视频| 99riav亚洲国产免费| 国产精品1区2区在线观看.| 欧美激情在线99| 日本熟妇午夜| 91久久精品国产一区二区成人| 国产成人freesex在线| 观看美女的网站| 91狼人影院| 18禁在线播放成人免费| 人人妻人人看人人澡| 九九热线精品视视频播放| 天堂av国产一区二区熟女人妻| 九色成人免费人妻av| 亚洲欧美日韩高清在线视频| 亚洲自偷自拍三级| 精品国内亚洲2022精品成人| 长腿黑丝高跟| 99精品在免费线老司机午夜| 亚洲熟妇中文字幕五十中出| 亚洲人成网站高清观看| 免费电影在线观看免费观看| 尤物成人国产欧美一区二区三区| 免费电影在线观看免费观看| 国产av一区在线观看免费| 国产成人91sexporn| 亚洲国产精品国产精品| www.色视频.com|