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

    黃河調(diào)水調(diào)沙期間入海泥沙沉積動(dòng)力過程的數(shù)值模擬*

    2017-01-12 10:59:52孫效功范德江
    關(guān)鍵詞:調(diào)沙落潮漲潮

    曾 明 , 孫效功, 范德江

    (1.中國海洋大學(xué)海底科學(xué)與探測(cè)技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,山東 青島 266100; 2.中國氣象科學(xué)研究院災(zāi)害天氣國家重點(diǎn)實(shí)驗(yàn)室, 北京 100081)

    黃河調(diào)水調(diào)沙期間入海泥沙沉積動(dòng)力過程的數(shù)值模擬*

    曾 明1, 孫效功2**, 范德江1

    (1.中國海洋大學(xué)海底科學(xué)與探測(cè)技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,山東 青島 266100; 2.中國氣象科學(xué)研究院災(zāi)害天氣國家重點(diǎn)實(shí)驗(yàn)室, 北京 100081)

    本文根據(jù)2007年黃河調(diào)水調(diào)沙期間黃河三角洲及其毗鄰海域的水文、泥沙調(diào)查數(shù)據(jù)資料,使用Mike 3三維水沙數(shù)值模型模擬了黃河調(diào)水調(diào)沙全過程(低流量、高流量、高含沙量)入海水沙輸運(yùn)過程,揭示了人造洪峰下水沙擴(kuò)散范圍和特征,分析了泥沙入海后的沉積動(dòng)力學(xué)機(jī)制。研究發(fā)現(xiàn):入海淡水和表層細(xì)顆粒泥沙主要以異輕羽狀流形式向外海和河口兩側(cè)沿岸擴(kuò)散,中、底層粗顆粒泥沙向外海擴(kuò)散非常有限,有少量向河口兩側(cè)沿岸擴(kuò)散,擴(kuò)散范圍小于表層;河口區(qū)在漲-落、落-漲潮流轉(zhuǎn)換時(shí)刻分別存在著潮流切變鋒,切變鋒均形成在淺水區(qū),而后向深水區(qū)移動(dòng),切變鋒的發(fā)生對(duì)泥沙向外海擴(kuò)散有著明顯的阻隔作用。

    黃河口;調(diào)水調(diào)沙;三維數(shù)值模擬;水沙輸運(yùn)規(guī)律;沉積動(dòng)力過程

    黃河是中國第二大河,以水少沙多聞名于世,多年年平均入海的泥沙約為10.8億t[1]。黃河水利委員會(huì)于2002年7月開始了第一次“調(diào)水調(diào)沙”試驗(yàn),該試驗(yàn)利用黃河干流的大型水庫,以制造人造洪峰的方式將占全年20%以上的徑流和30%以上的泥沙在20天左右輸送入海[2]。調(diào)水調(diào)沙實(shí)施以來,現(xiàn)行河口的葉瓣由侵蝕狀態(tài)轉(zhuǎn)變?yōu)橛俜e狀態(tài)[3],河口泥沙異重流變的不常見,異輕流成為泥沙向海擴(kuò)散的主要方式[4]。調(diào)水調(diào)沙的實(shí)施改變了黃河入海水沙的自然形態(tài),改變了泥沙向海擴(kuò)散的方式,對(duì)黃河三角洲近岸帶地形、地貌演化以及近海海洋生態(tài)環(huán)境產(chǎn)生了深遠(yuǎn)影響[5-10]。

    近年來,有的學(xué)者基于現(xiàn)場(chǎng)調(diào)查數(shù)據(jù),研究了調(diào)水調(diào)沙期間黃河入海水沙的擴(kuò)散規(guī)律,對(duì)調(diào)水調(diào)沙期間黃河入海水沙的基本擴(kuò)散態(tài)勢(shì)有了一定的認(rèn)識(shí)[11-13];有的也基于衛(wèi)星遙感資料進(jìn)一步研究了黃河河口區(qū)的蝕-淤狀態(tài),并得到了新的認(rèn)識(shí)[3];個(gè)別學(xué)者還利用數(shù)值模擬方法對(duì)黃河洪、枯季的水沙輸運(yùn)[14]和調(diào)水調(diào)沙期間入海泥沙輸運(yùn)[15]進(jìn)行了模擬,也取得了較好的結(jié)果。迄今為止,雖然許多學(xué)者對(duì)調(diào)水調(diào)沙期間黃河入海水沙的擴(kuò)散規(guī)律以及河口沉積動(dòng)力過程進(jìn)行了系列研究,但尚缺乏對(duì)調(diào)水調(diào)沙全過程河口沉積動(dòng)力過程的完整認(rèn)識(shí)。

    本文立足于調(diào)水調(diào)沙期間黃河三角洲及其毗鄰海域的水文泥沙調(diào)查資料,結(jié)合實(shí)測(cè)資料驗(yàn)證的三維水沙數(shù)值模型,通過對(duì)黃河調(diào)水調(diào)沙輸送入海全過程的數(shù)值模擬,力圖詳細(xì)了解調(diào)水調(diào)沙期間入海水沙的輸運(yùn)規(guī)律,探討其沉積動(dòng)力過程的控制機(jī)制。

    1 現(xiàn)場(chǎng)調(diào)查及數(shù)據(jù)處理

    1.1 調(diào)水調(diào)沙期間入海水沙形態(tài)

    2007年調(diào)水調(diào)沙開始于6月19日,一直延續(xù)到7月7日,共18 d。這期間,入海水沙不同步,2007年黃河調(diào)水調(diào)沙期間利津站入海逐日徑流量、輸沙率如圖1所示,據(jù)此可以分為3個(gè)階段: 6月19~6月25日的調(diào)水調(diào)沙之前階段(低流量、低沙量),其入海徑流量不超過300m3/s,輸沙率為0.1t/s左右;6月26日~7月6日的調(diào)水調(diào)沙第一階段(高流量、中沙量),7月3日徑流量達(dá)到最高值,為3800m3/s,其輸沙率為46.4t/s;7月7日的調(diào)水調(diào)沙第二階段(中流量、高沙量),其徑流量為2310m3/s,輸沙率達(dá)到調(diào)水調(diào)沙期間最高值,為73.1t/s。其中6月26—7月7日入海水量均超過2000m3/s,7月3日入海水量最高,達(dá)到3800m3/s。

    1.2 調(diào)查站位布設(shè)和數(shù)據(jù)分析

    2007年6月29—7月6日黃河調(diào)水調(diào)沙期間,在黃河三角洲海域進(jìn)行了3斷面3船同步25h(1個(gè)潮周期)的連續(xù)觀測(cè),另布設(shè)了24個(gè)大面站(見圖2)。采用YSI公司1000kHz ADCP測(cè)流速、流向,采用ALEC水質(zhì)儀及RBR自容式CTD測(cè)溫度、鹽度及濁度,獲得了研究海域的流速、流向、溫度、鹽度、濁度等實(shí)測(cè)數(shù)據(jù)。同時(shí),在6個(gè)預(yù)設(shè)層位(表層,0.2,0.4,0.6,0.8h,底層)采集水樣各500mL左右,每小時(shí)1次,經(jīng)抽濾、洗鹽后,并在50℃條件下烘干24h,然后使用1/100000天平稱重,獲得水體懸浮泥沙含量。

    2 黃河口三維水沙數(shù)值模型的率定和驗(yàn)證

    采用MIKE 3三維水沙模型來模擬黃河調(diào)水調(diào)沙期間入海水沙的輸運(yùn)。MIKE 3是丹麥水力研究所(DHI)開發(fā)的MIKE系列水沙模型的三維版本,有著幾十年的開發(fā)和運(yùn)用經(jīng)驗(yàn),能較好地模擬近海水動(dòng)力和泥沙輸運(yùn),在全世界河口海岸有諸多成功的應(yīng)用[16-19]。該模型包括水動(dòng)力模塊(HD)、對(duì)流模塊(AD)、波浪模塊(SW)、粒子追蹤模塊(PT)、泥輸運(yùn)模塊(MT)、沙輸運(yùn)模塊(ST)等幾大部分。

    2.1 三維水動(dòng)力模塊

    笛卡爾坐標(biāo)系統(tǒng)中的水動(dòng)力的連續(xù)方程如下:

    x、y向水平動(dòng)量方程:

    其中:t為時(shí)間;x、y、z為笛卡爾坐標(biāo)軸;η為表層潮位;d為靜水深;h=η+d為表示總水深;u、v、w為x、y、z方向的速度分量;f=2Ωsinφ為科氏力系數(shù);g為重力加速度;ρ為水的密度;Sxx、Sxy、Syx、Syy為輻射應(yīng)力張量要素;Vt為垂向擾動(dòng)黏滯系數(shù);Pa為大氣壓力;ρ0為水的相對(duì)密度;S為點(diǎn)源的流量;us、vs為水體間相對(duì)運(yùn)動(dòng)的的速度分量。

    2.2 泥沙輸運(yùn)模塊

    Mike3泥沙模塊中,主要研究的是粘性和非粘性泥沙的運(yùn)移過程,該模塊對(duì)泥質(zhì)和砂質(zhì)等不同類型的沉積物分別進(jìn)行定義,包括絮凝、沉降、再懸浮、層間運(yùn)移等過程。其中,黏性泥沙是指平均粒徑小于60μm的泥沙顆粒,而非黏性泥沙則指平均粒徑在60~500μm的泥沙顆粒。由于渤海沉積物及黃河入海泥沙的平均粒徑大部分都是小于60μm,因此本文采用Mike3粘性泥沙模塊(MT)來進(jìn)行數(shù)值研究。

    泥沙運(yùn)移數(shù)值模擬是建立在水動(dòng)力模塊中對(duì)流-擴(kuò)散計(jì)算基礎(chǔ)之上的,其中泥沙運(yùn)移的控制方程如下:

    其中:c為水體含沙量(kg/m3);t為時(shí)間(s);x、y、z為原點(diǎn)o置于某一水平基面,z軸垂直向上的直角坐標(biāo)系坐標(biāo);u、v、w為空間流速矢量分別沿x、y、z方向的速度分量(m/s);Dx、Dy、Dz為泥沙沿x、y、z方向上的紊動(dòng)擴(kuò)散系數(shù)(m2/s);S為沉積/侵蝕量(kg·m-3·s-1); QL為泥沙源強(qiáng)(m3·s-1·m-2); CL為泥沙源強(qiáng)濃度(kg·m-3)。

    2.3 水沙數(shù)值模型的設(shè)置及驗(yàn)證

    2.3.1 模型設(shè)置

    (1)水深和岸界

    模型使用的是ETOP5的1/12度水深資料,渤海灣及黃河口附近海域水深由2000年渤海灣實(shí)測(cè)水深和2004年黃河口36斷面實(shí)測(cè)水深替換,黃河三角洲岸線采用2004年資料(見圖3)。

    (2) 開邊界輸入

    開邊界:利用遼寧登沙河(A點(diǎn))、山東雞鳴島(B點(diǎn))潮位觀測(cè)資料調(diào)和求得的M2、S2、K1和O14個(gè)主要分潮調(diào)和常數(shù)值輸入計(jì)算。

    (3)徑流邊界條件

    本論文模型里只考慮黃河注入渤海,黃河徑流邊界條件采用的是黃河利津站2007年6和7月逐日徑流量、輸沙率數(shù)據(jù)。

    (4)渤海底質(zhì)類型資料

    渤海底質(zhì)類型資料來自2005年出版的《中國東部海域海底沉積物成因環(huán)境圖》[20]數(shù)字化結(jié)果。渤海底質(zhì)類型以粉砂質(zhì)粘土和粘土質(zhì)粉砂為主,細(xì)砂主要分布在水動(dòng)力較強(qiáng)的渤海海峽北部和遼東淺灘,黃河水下三角洲的表層沉積物以粉砂質(zhì)粘土為主。

    (5)溫鹽條件

    模型初始鹽度值設(shè)為32 psu,徑流鹽度值為0;研究區(qū)初始溫度設(shè)為20℃,模型計(jì)算期間溫度不變。

    2.3.2 模型驗(yàn)證 模擬結(jié)果經(jīng)2007年調(diào)水調(diào)沙期間三角洲近岸海域連續(xù)站實(shí)測(cè)資料驗(yàn)證,證明其潮位、流速、流向、懸沙濃度基本能反映調(diào)水調(diào)沙期間黃河三角洲及其鄰近海域的水動(dòng)力泥沙輸運(yùn)特征,可用于入海水沙的擴(kuò)散規(guī)律分析。

    (1)潮位驗(yàn)證

    模擬得到的潮位變化和三角洲近岸海域6個(gè)連續(xù)站位的潮位實(shí)測(cè)數(shù)據(jù)吻合良好,基本趨勢(shì)一致,但是實(shí)測(cè)數(shù)據(jù)波動(dòng)較大,表明水文模型模擬潮位效果較好。潮位驗(yàn)證如圖4。

    (2)流速、流向驗(yàn)證

    模型得到的潮流的流速、流向與實(shí)測(cè)結(jié)果符合度較高,兩者的趨勢(shì)完全一致(見圖5)。

    (3)懸沙濃度驗(yàn)證

    模型得到的懸沙濃度與測(cè)試懸沙濃度在不同層位符合程度不一,在底層,模擬結(jié)果和實(shí)測(cè)結(jié)果相關(guān)性較好,兩者的變化趨勢(shì)基本一致。在表層,兩者符合度差,實(shí)測(cè)數(shù)據(jù)明顯高于模型的數(shù)據(jù)。如圖6。

    3 調(diào)水調(diào)沙全過程入海水沙擴(kuò)散的數(shù)值模擬

    3.1 調(diào)水調(diào)沙入海淡水的輸運(yùn)

    3.1.1 低流量狀態(tài)下黃河入海淡水的擴(kuò)散(調(diào)水調(diào)沙之前) 低流量漲潮時(shí),表層淡水主要向河口南側(cè)擴(kuò)散,底層水體基本不受淡水?dāng)U散影響,始終為海水控制(見圖7A、B);落潮時(shí)淡水隨落潮流向西北方向擴(kuò)散,擴(kuò)散范圍大于漲潮時(shí)(見圖7C、D)。從漲、落潮時(shí)河口剖面鹽度分布可以看出(見圖7E、F),淡水浮于海水表層向海有少量擴(kuò)散,落潮時(shí)擴(kuò)散量高于漲潮時(shí),口門外水體垂向分層較為明顯,河道內(nèi)海水入侵程度較漲潮時(shí)有所減弱。

    3.1.2 高流量時(shí)入海淡水的擴(kuò)散(調(diào)水調(diào)沙第一階段) 6月26—7月7日黃河逐日入海水量均在2000m3/s以上,其中7月3日入海水量最高,達(dá)到3 800 m3/s。通過漲潮時(shí)鹽度等值線可以看出,高流量時(shí)表、底層淡水?dāng)U散范圍均大于低流量時(shí),表層擴(kuò)散范圍遠(yuǎn)大于底層,低于4 psu的低鹽度水以舌狀偏向河口南側(cè),10 psu的等鹽度線能到達(dá)原清水溝廢棄河口,底層淡水主要集中在河口附近,少量向河口兩側(cè)沿岸擴(kuò)散,向海擴(kuò)散有限(見圖8A、B);落潮時(shí)表層擴(kuò)散范圍遠(yuǎn)大于底層,受落潮流影響,淡水?dāng)U散范圍較漲潮時(shí)整體向北偏移,低于4 psu的低鹽度水以舌狀偏向河口北側(cè),10 psu的等鹽度線能到達(dá)神仙溝外海(見圖8C、D)。從河口斷面的鹽度分布可以看出,因流量較大,漲潮時(shí)河道內(nèi)基本為淡水控制,口門外淡水?dāng)U散范圍較大,淡水浮于表層向外擴(kuò)散,鹽度垂向分層明顯, 落潮時(shí)河道內(nèi)基本為淡水控制,落潮時(shí)淡水?dāng)U散量比漲潮時(shí)大,口門處垂向混合較為均勻(見圖8E、F)。

    高含沙量時(shí)(調(diào)水調(diào)沙第二階段),入海淡水量較高流量時(shí)(調(diào)水調(diào)沙第一階段)變小,淡水隨漲、落潮流的擴(kuò)散方式與高流量時(shí)相似,其詳細(xì)的擴(kuò)散過程在此略過。

    3.2 調(diào)水調(diào)沙期間泥沙的輸運(yùn)

    3.2.1 低流量時(shí)黃河入海泥沙的擴(kuò)散(調(diào)水調(diào)沙之前) 低流量時(shí)入海泥沙量較低,受漲、落潮流影響,漲潮時(shí)表層泥沙主要向河口南側(cè)擴(kuò)散,落潮時(shí)則主要向河口北側(cè)擴(kuò)散,落潮時(shí)入海泥沙量和擴(kuò)散范圍大于漲潮時(shí)(見圖9A、B、C、D、E、F)。

    (A-漲潮時(shí)表層鹽度分布A-salinity distribution of the surface layer at flood tide phase;B-漲潮時(shí)底層鹽度分布B-salinity distribution of the bottom layer at flood tide phase;C-落潮時(shí)表層鹽度分布C-salinity distribution of the surface layer at ebb tide phase;D-落潮時(shí)底層鹽度分布D-salinity distribution of the bottom layer at ebb tide phase;E-漲潮時(shí)河口剖面鹽度E-salinity distribution of cross-section in the estuary at flood tide phase;F-落潮時(shí)河口剖面鹽度F-salinity distribution of cross-section in the estuary at ebb tide phase.)

    圖7 低流量狀態(tài)下入海淡水的擴(kuò)散

    Fig.7 Disperse patterns of river-delivered fresh water in low-flow period

    (A-漲潮時(shí)表層鹽度分布A-salinity distribution of the surface layer at flood tide phase;B-漲潮時(shí)底層鹽度分布B-salinity distribution of the bottom layer at flood tide phase;C-落潮時(shí)表層鹽度分布C-salinity distribution of the surface layer at ebb tide phase;D-落潮時(shí)底層鹽度分布D-salinity distribution of the bottom layer at ebb tide phase;E-漲潮時(shí)河口剖面鹽度E-salinity distribution of cross-section in the estuary at flood tide phase;F-落潮時(shí)河口剖面鹽度F-salinity distribution of cross-section in the estuary at ebb tide phase.)

    圖8 高流量時(shí)入海淡水的擴(kuò)散

    Fig.8 Disperse patterns of river-delivered fresh water in high-flow period

    (A-漲潮時(shí)表層泥沙濃度A-distribution of the Suspended Sediment Concentration(SSC) in the surface layer at flood tide phase;B-漲潮時(shí)底層泥沙濃度B-distribution of the Suspended Sediment Concentration(SSC) in the bottom layer at flood tide phase;C-落潮時(shí)表層泥沙濃度C-distribution of the Suspended Sediment Concentration(SSC) in the surface layer at ebb tide phase;D-落潮時(shí)底層泥沙濃度D-distribution of the Suspended Sediment Concentration(SSC) in the bottom layer at ebb tide phase;E-漲潮時(shí)河口斷面泥沙濃度E-the Suspended Sediment Concentration(SSC) distribution of cross-section in the estuary at flood tide phase;F-落潮時(shí)河口斷面泥沙濃度F-the Suspended Sediment Concentration(SSC) distribution of cross-section in the estuary at ebb tide phase.)

    圖9 低流量時(shí)入海泥沙擴(kuò)散

    Fig.9 Disperse patterns of river-delivered suspended sediment in low-flow period

    3.2.2 高流量時(shí)黃河入海泥沙的擴(kuò)散(調(diào)水調(diào)沙第一階段) 調(diào)水調(diào)沙第一階段時(shí),入海水沙量較大,漲、落潮時(shí)表、底層泥沙擴(kuò)散范圍遠(yuǎn)大于較低流量時(shí),底層泥沙濃度高于表層(見圖10A、B、C、D)。漲潮時(shí)表層泥沙主要向河口南側(cè)擴(kuò)散,能擴(kuò)散至原清水溝廢棄河口近岸海域,而落潮時(shí)表層泥沙隨落潮流可擴(kuò)散至神仙溝近岸海域;底層泥沙漲潮時(shí)主要向河口南側(cè)擴(kuò)散,落潮時(shí)向河口北側(cè)擴(kuò)散,底層泥沙向外海擴(kuò)散有限,因河口區(qū)漲潮流大于落潮流,從而導(dǎo)致漲潮時(shí)泥沙擴(kuò)散距離大于落潮時(shí)。通過高流量漲、落潮時(shí)河口斷面泥沙濃度分布也可看出(見圖10E、F),入海泥沙向外海擴(kuò)散有限,主要沉積在河口近岸海域,底層泥沙濃度明顯高于表、中層。

    3.2.3 高含沙量狀態(tài)下黃河入海泥沙的擴(kuò)散(調(diào)水調(diào)沙第二階段) 本次調(diào)水調(diào)沙接近尾聲時(shí),入海泥沙量在7月7日達(dá)到了最高值,含沙量為31.645kg/m3。這一階段的泥沙擴(kuò)散特征與高流量時(shí)類似,主要沿岸線向河口兩側(cè)擴(kuò)散,受徑流量減小及河口潮流動(dòng)力減弱的影響,高含沙量狀態(tài)下入海泥沙擴(kuò)散距離小于高流量時(shí)(圖11A、B、C、D、E、F)。

    (A-漲潮時(shí)表層泥沙濃度A-distribution of the Suspended Sediment Concentration(SSC) in the surface layer at flood tide phase;B-漲潮時(shí)底層泥沙濃度B-distribution of the Suspended Sediment Concentration(SSC) in the bottom layer at flood tide phase;C-落潮時(shí)表層泥沙濃度C-distribution of the Suspended Sediment Concentration(SSC) in the surface layer at ebb tide phase;D-落潮時(shí)底層泥沙濃度D-distribution of the Suspended Sediment Concentration(SSC) in the bottom layer at ebb tide phase;E-漲潮時(shí)河口斷面泥沙濃度E-the Suspended Sediment Concentration(SSC) distribution of cross-section in the estuary at flood tide phase;F-落潮時(shí)河口斷面泥沙濃度F-the Suspended Sediment Concentration(SSC) distribution of cross-section in the estuary at ebb tide phase.)

    圖10 高流量時(shí)入海泥沙擴(kuò)散

    Fig.10 Disperse patterns of river-delivered suspended sediment in high-flow period

    4 調(diào)水調(diào)沙期間黃河入海水沙擴(kuò)散的動(dòng)力學(xué)分析

    4.1 河口羽狀流

    黃河入海水沙以異輕、異重羽狀流兩種形式擴(kuò)散和沉積。黃河口區(qū)洪水期水體包含潮汐鹽水層(主要分布在三角洲前緣中下部深水區(qū),當(dāng)位于水體底部時(shí),含沙量較高 )、異輕羽狀流(主要分布在水體中上層,含沙量較低)和異重羽狀流層(主要分布在水體中下層、底層,含沙量高)[ 21-22]。

    本文通過數(shù)值模擬發(fā)現(xiàn)異輕羽狀流是黃河調(diào)水調(diào)沙期間入海泥沙向海輸送的重要形式。在黃河低流量狀態(tài)下,入海泥沙主要是以異輕羽狀流形式向外輸送(見圖9E、F),表層含沙量明顯高于底層;在高流量和高含沙量狀態(tài)下,河口區(qū)底層泥沙濃度高于表層,這可能是由于入海粗顆粒泥沙含量升高或發(fā)生了泥沙絮凝作用所導(dǎo)致的,而異輕羽狀流亦是細(xì)顆粒泥沙向外輸送的重要形式(見圖10E、11E)。

    (A-漲潮時(shí)表層泥沙濃度 A-distribution of the Suspended Sediment Concentration(SSC) in the surface layer at flood tide phase;B-漲潮時(shí)底層泥沙濃度B-distribution of the Suspended Sediment Concentration(SSC) in the bottom layer at flood tide phase;C-落潮時(shí)表層泥沙濃度C-distribution of the Suspended Sediment Concentration(SSC) in the surface layer at ebb tide phase;D-落潮時(shí)底層泥沙濃度D-distribution of the Suspended Sediment Concentration(SSC) in the bottom layer at ebb tide phase;E-漲潮時(shí)河口斷面泥沙濃度 E-the Suspended Sediment Concentration(SSC) distribution of cross-section in the estuary at flood tide phase;F-落潮時(shí)河口斷面泥沙濃度 F-the Suspended Sediment Concentration(SSC) distribution of cross-section in the estuary at ebb tide phase)

    圖11 高含沙量下入海泥沙擴(kuò)散

    Fig.11 Disperse patterns of river-delivered suspended sediment in high sediment concentration period

    4.2 河口切變鋒

    在河口區(qū)漲、落潮轉(zhuǎn)換階段,受岸線和海底地形的影響,出現(xiàn)了內(nèi)漲-外落型和外漲-內(nèi)落型切變鋒,Wang等[23]根據(jù)實(shí)測(cè)資料分析,認(rèn)為黃河口切變鋒發(fā)生在潮流轉(zhuǎn)流時(shí),形成于淺水區(qū),并向深水區(qū)移動(dòng),在一個(gè)潮周期內(nèi)歷時(shí)約4h。本文數(shù)值模型顯示在調(diào)水調(diào)沙期間黃河口近岸海域漲-落、落-漲潮流轉(zhuǎn)換時(shí)刻皆存在著河口切變鋒,其切變鋒在漲、落潮時(shí)皆持續(xù)約2h。

    內(nèi)漲-外落型和外漲-內(nèi)落型切變鋒均形成在淺水區(qū),而后向深水區(qū)移動(dòng)(見圖12A、B、C、D)。此時(shí)研究海區(qū)潮流速度較小,三角洲近岸區(qū)流向一致,與外海區(qū)流向基本相反,兩者間存在著流速低值區(qū),入海泥沙主要是沿岸輸運(yùn)。

    河口斷面鹽度和懸沙濃度分布(見圖12E、F、G、H)所示,切變鋒兩側(cè)水沙擴(kuò)散形態(tài)不同,內(nèi)漲-外落型切變鋒向岸側(cè)處于漲潮狀態(tài),口門外淺水區(qū)鹽水入侵,淡水浮于表層,向海側(cè)處于落潮狀態(tài);外漲-內(nèi)落型切變鋒向岸側(cè)處于落潮狀態(tài),入海淡水?dāng)U散量增加,而向海側(cè)處于漲潮狀態(tài),河道內(nèi)泥沙向海擴(kuò)散量高于內(nèi)漲-外落型切變鋒時(shí)。切變鋒的存在阻隔了入海泥沙的向海擴(kuò)散,使得入海泥沙主要沉積在近岸淺水區(qū)。

    (A-內(nèi)漲-外落切變鋒1時(shí)刻表層流場(chǎng)A-velocity distribution of the surface layer flow for the shear front of inner-flood-outer-ebb at the 1st hour;B-內(nèi)漲-外落切變鋒2時(shí)刻表層流場(chǎng)B-velocity distribution of the surface layer flow for the shear front of inner-flood-outer-ebb at the 2nd hour;C-外漲-內(nèi)落切變鋒1時(shí)刻表層流場(chǎng)C-velocity distribution of the surface layer flow for the shear front of inner-ebb-outer-flood at the 1st hour;D-外漲-內(nèi)落切變鋒2時(shí)刻表層流場(chǎng)D-velocity distribution of the surface layer flow for the shear front of inner-ebb-outer-flood at the 2nd hour;E-內(nèi)漲-外落切變鋒第1時(shí)刻河口斷面鹽度E-salinity distribution of cross-section in the estuary for the shear front of inner-flood-outer-ebb at the 1st hour;F-外漲-內(nèi)落切變鋒第1時(shí)刻河口斷面鹽度F-salinity distribution of cross-section in the estuary for the shear front of inner-ebb-outer-flood at the 1st hour;G-內(nèi)漲-落第1時(shí)刻河口斷面泥沙濃度G-Suspended Sediment Concentration(SSC) distribution of cross-section in the estuary for the shear front of inner-flood-outer-ebb at the 1st hour;H-外漲-內(nèi)落第1時(shí)刻河口斷面泥沙濃度H-Suspended Sediment Concentration(SSC) distribution of cross-section in the estuary for the shear front of inner-ebb-outer-flood at the 1st hour)

    圖12 河口切變鋒發(fā)育狀態(tài)

    Fig.12 The generation of different patterns of the shear front

    5 結(jié)論

    本文基于2007年黃河調(diào)水調(diào)沙期間黃河三角洲及其鄰近海域的實(shí)際調(diào)查資料,利用Mike 3三維水沙數(shù)值模型模擬了黃河調(diào)水調(diào)沙全過程(低流量、高流量、高含沙量)的水沙輸運(yùn)特征,并探討了黃河入海水沙輸運(yùn)的動(dòng)力學(xué)機(jī)制。主要結(jié)論如下:

    (1)黃河調(diào)水調(diào)沙期間入海淡水主要浮在海水表層向海擴(kuò)散,低流量(調(diào)水調(diào)沙前)漲潮時(shí)淡水隨漲潮流向河口南側(cè)擴(kuò)散,落潮時(shí)則隨落潮流向河口北側(cè)擴(kuò)散,落潮時(shí)淡水?dāng)U散范圍大于漲潮時(shí),口門外水體垂向分層較為明顯;高流量時(shí)(調(diào)水調(diào)沙第一階段)表、底層淡水?dāng)U散范圍均大于低流量時(shí),淡水浮于表層向外擴(kuò)散,擴(kuò)散范圍大于底層,口門外鹽度分層明顯;高含沙量時(shí)(調(diào)水調(diào)沙第二階段)入海淡水的擴(kuò)散方式與高流量時(shí)相似,擴(kuò)散幅度較高流量時(shí)變小。

    (2)黃河低流量時(shí)(調(diào)水調(diào)沙前)入海泥沙以細(xì)顆粒泥沙為主,主要以異輕羽狀流形式向海擴(kuò)散,漲潮時(shí)泥沙主要向河口南側(cè)擴(kuò)散,落潮時(shí)則主要向河口北側(cè)擴(kuò)散,泥沙擴(kuò)散范圍和擴(kuò)散量大于漲潮時(shí);高流量時(shí)入海泥沙量較大,漲、落潮時(shí)表、底層泥沙擴(kuò)散范圍遠(yuǎn)大于低流量時(shí),底層泥沙受河口切變鋒阻隔影響,向外海擴(kuò)散有限。高含沙量時(shí)泥沙擴(kuò)散特征與高流量時(shí)類似,主要沿岸線向河口兩側(cè)擴(kuò)散,高含沙量狀態(tài)下入海泥沙擴(kuò)散距離小于高流量時(shí)。

    (3)在漲-落、落-漲潮流轉(zhuǎn)換時(shí)刻分別存在著河口潮流切變鋒,這兩型切變鋒均形成在淺水區(qū),而后向深水區(qū)移動(dòng),在一個(gè)潮周期內(nèi)歷時(shí)約4 h,此時(shí)研究海區(qū)潮流速度較小,三角洲近岸區(qū)流向一致,與外海區(qū)流向基本相反,兩者間存在著流速低值區(qū),切變鋒的發(fā)生對(duì)泥沙向外海擴(kuò)散有著明顯的阻隔作用。

    [1] Milliman J D, Meade R H. World-wide delivery of sediment to the oceans[J]. Journal of Geology, 1983, 91(1):1-21.

    [2] 畢乃雙. 黃河三角洲毗鄰海域懸浮泥沙擴(kuò)散和季節(jié)性變化及沖淤效應(yīng)[D]. 青島:中國海洋大學(xué), 2009. Bi N S. Suspended Sediment Dispersal off the Huanghe(Yellow River) Delta and in Its Adjacent Bohai Sea, Its Seasonal Variation and Effect on the Delta Erosion-Accumulation[D]. Qingdao: Ocean University of China, 2009.

    [3] Bi N, Wang H, Yang Z. Recent changes in the erosion-accretion patterns of the active Huanghe(Yellow River) delta lobe caused by human activities[J]. Continental Shelf Research, 2014, 90: 70-78.

    [4] Wang H, Bi N, Saito Y, et al. Recent changes in sediment delivery by the Huanghe(Yellow River) to the sea: causes and environmental implications in its estuary[J]. Journal of Hydrology, 2010, 391(3): 302-313.

    [5] 李少文,張瑩,李凡,等. 調(diào)水調(diào)沙對(duì)黃河口海域大型底棲動(dòng)物群落的影響[J]. 環(huán)境科學(xué)研究, 2015, 28(2): 259-266. Li S W, Zhang Y, Li F, et al. Effects of water and sediment discharge regulation on macrobenthic community in the Yellow River Estuary [J]. Research of Environmental Sciences, 2015, 28(2): 259-266.

    [6] 王婷. 2002~2004年及調(diào)水調(diào)沙期間黃河下游營養(yǎng)鹽的變化特征[D]. 青島:中國海洋大學(xué),2007. Wang T. The Variations of Nutrients in the Lower Main Channel of the Yellow River from 2002 to 2004 and Water-Sediment Regulation[D]. Qingdao: Ocean University of China, 2007.

    [7] 于帥. 黃河調(diào)水調(diào)沙影響下河口入海泥沙擴(kuò)散及地貌效應(yīng)[D]. 青島:中國海洋大學(xué),2014. Yu S. Suspended Sediment Dispersal off the Huanghe(Yellow River) River Mouth and Its Morphological Effects Under Impact of the Water-Sediment Regulation Scheme [D]. Qingdao: Ocean University of China, 2014.

    [8] 鄒榮婕,鄧旭修,王斌,等. 黃河調(diào)水調(diào)沙對(duì)黃河口海域雙酚A的影響[J]. 生態(tài)學(xué)報(bào), 2015, 35(2): 263-269. Zou R J, Deng X X, Wang B, et al. Effect of water-sediment flushing events on bisphenol A contamination in the Yellow River estuary [J]. Acta Ecologica Sinica, 2015, 35(2): 263-269.

    [9] 顧培培,張桂玲,李佩佩,等. 調(diào)水調(diào)沙對(duì)黃河口及鄰近海域溶存甲烷的影響[J]. 中國環(huán)境科學(xué), 2011, 31(11):1821-1828. Gu P P, Zhang G L, Li P P, et al. Effect of the water-sediment regulation on dissolved methane in the lower Yellow River estuary and it's adjacent marine area [J]. Environmental Science, 2011,31(11): 1821-1828.

    [10] 姚前前. 小浪底調(diào)水調(diào)沙對(duì)黃河口來水來沙條件和海岸線演變規(guī)律的影響研究[D]. 泰安:山東農(nóng)業(yè)大學(xué),2014. Yao Q Q. Influence of Flow/Sediment Regulation in Xiaolangdi on Incoming Flow/Sediment Condition and Coastline Evolution of the Yellow River Estuary [D]. Taian: Shandong Agricultural University, 2014.

    [11] 史文靜. 黃河口懸浮泥沙擴(kuò)散規(guī)律及其數(shù)值模擬研究[D]. 青島:中國海洋大學(xué),2008. Shi W J. 3D Numerical Simulation and Analysis of the Transportation of Suspended Sediment in the Huanghe Estuary[D]. Qingdao: Ocean University of China, 2008.

    [12] Wang H J, Yang Z S, Bi N S, et al. Rapid shifts of the river plume pathway off the Huanghe(Yellow) River mouth in response to water-sediment regulation scheme in 2005[J]. Chinese Science Bulletin, 2005, 50(24): 2878-2884.

    [13] 劉鋒,陳沈良,周永東,等. 黃河2009年調(diào)水調(diào)沙期間河口水動(dòng)力及懸沙輸移變化特征[J]. 泥沙研究,2010, 6: 1-8. Liu F, Chen S L, Zhou Y D, et al. Effect of water-sediment regulation in Yellow River on hydrodynamics and suspended sediment transport in its estuary [J]. Journal of Sediment Research, 2010, 6: 1-8.

    [14] 胡小雷,陳沈良,劉小喜,等. 2012年調(diào)水調(diào)沙期間黃河口水沙擴(kuò)散途徑及范圍[J]. 泥沙研究, 2014, 3: 49-56. Hu X L, Chen S L, Liu X X, et al. Diffusion path and range of water flow and sediment in Yellow River estuary during water-sediment regulation in 2012 [J]. Journal of Sediment Research, 2014, 3: 49-56.

    [15] 李平,豐愛平,陳義中,等. 2005年黃河調(diào)水調(diào)沙期間入海泥沙擴(kuò)散研究[J]. 海洋湖沼通報(bào), 2010, 4: 72-78. Li P, Feng A P, Chen Y Z, et al. The diffusion studies of sediment into sea during the period of Yellow River water and sediment regulation, 2005 [J]. Transactions of Oceanology and Limnology, 2010, 4: 72-78.

    [16] 許婷. MIKE 21 HD計(jì)算原理及應(yīng)用實(shí)例[J]. 港工技術(shù),2010,47(5): 1-5. Xu T. Calculation principle and application example of a two-dimensional flow model-MIKE21 HD [J]. Port Engineering Technology, 2010, 47(5): 1-5.

    [17] 陳雪峰,王桂萱. MIKE 21計(jì)算軟件及其在長興島海域改造工程上的應(yīng)用[J]. 大連大學(xué)學(xué)報(bào),2007,28(6):93-98. Chen X F, Wang G X. MIKE21 software and its application on the offshore reconstruction engineering of Changxing Islands [J]. Journal of Dalian University, 2007, 28(6): 93-98.

    [18] 朱龍海. 遼東淺灘潮流沉積動(dòng)力地貌學(xué)研究[D].青島:中國海洋大學(xué),2010. Zhu L H. Research on Dynamic Geomorphology of Tidal Deposition in Liaodong Shoal[D]. Qingdao: Ocean University of China, 2010.

    [19] Lumborg U, Pejrup M. Modelling of cohesive sediment transport in a tidal lagoon—an annual budget Original[J]. Marine Geology, 2005, 218(1-4):1-16.

    [20] 李廣雪,楊子庚,劉勇,等. 中國東部海域海底沉積物成因環(huán)境圖[M]. 北京:科學(xué)出版社,2005. Li G X, Yang Z G, Liu Y, et al. Genetic Environment of Sea Floor Sediments in China Eastern Sea Areas[M]. Beijing: Science Press, 2005.

    [21] Wiseman W J, Fan Y B, Bornhold B D. Suspended sediment advection by tidal current of the Huanghe Delta[J]. Geo-Marine Letters, 1986, (6):107-113.

    [22] 李廣雪. 黃河入海泥沙擴(kuò)散與河海相互作用[J]. 海洋地質(zhì)與第四紀(jì)地質(zhì), 1999, 19(3): 1-10. Li G X. Suspended sediment dispersal and interaction of river-sea off the Yellow River mouth[J]. Marine Geology & Quaternary Geology, 1999, 19(3): 1-10.

    [23] Wang H J, Yang Z S, Li Y H, et al. Dispersal pattern of suspended sediment in the shear frontal zone off the Huanghe(Yellow River) mouth[J]. Continental Shelf Research, 2007, 27:854-871.

    責(zé)任編輯 徐 環(huán)

    Numerical Simulating on Sediment Transportation in the Estuary and its Adjacent Areas During the Water and Sediment Regulation Period of the Huanghe River

    ZENG Ming1, SUN Xiao-Gong2, FAN De-Jiang1

    (1.Key Laboratory of Submarine Geosciences and Technology, Ministry of Education, Ocean University of China, Qingdao 266100, China; 2.State Key Laboratory of Severe Weather,Chinese Academy of Meteorological Sciences, Beijing 100081, China)

    Based on a synthetic survey which was carried out in the Huanghe estuary and its adjacent area during the water-sediment regulation period in the summer season of 2007, the MIKE 3 numerical model was applied to the Huanghe estuary and its adjacent areas. With the model, the transportation processes of the fresh water and suspended sediment were successfully simulated in the Huanghe estuary, and their dispersal patterns in different discharge phase (Low-flow, high-flow, high sediment concentration) were analyzed carefully. The results of numerical model showed that whether the river discharge was low or high, the fresh water could disperse in the surface layer as the hypopycnal flow. When the river discharge was low, river-delivered suspended sediment dispersed towards the offshore was mainly through the hypopycnal flow; When the river discharge was high, the hypopycnal flow was also an important controlling factor for the dispersal of fine sediment, most of the coarse sediment deposited in the nearshore area of the Huanghe estuary, a small amount of coarse sediment was transported along the coastline. Two types of shear front were identified in the offshore region near the river mouth at the conversion time of tidal flow, which was formed in the shallow water, and then moved to the deep water. Sediment was transported along the coast, the offshore transportation was significant constrained.

    Huanghe Estuary; water-sediment regulation; 3D numerical simulation; dispersal patterns of fresh water and suspended sediment; sedimentary dynamic processes

    國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃項(xiàng)目(2010CB951202)資助

    2015-10-27;

    2015-12-01

    曾 明(1985-),男,博士生。研究方向:海洋沉積動(dòng)力學(xué)。

    ** 通訊作者:E-mail:sunxg@camscma.cn

    P641.2; X523

    A

    1672-5174(2017)04-081-12

    10.16441/j.cnki.hdxb.20150367

    曾明, 孫效功, 范德江. 黃河調(diào)水調(diào)沙期間入海泥沙沉積動(dòng)力過程的數(shù)值模擬[J]. 中國海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2017, 47(4): 81-92.

    ZENG Ming, SUN Xiao-Gong, FAN De-Jiang. Numerical simulating on sediment transportation in the estuary and its adjacent areas during the water and sediment regulation period of the Huanghe River[J]. Periodical of Ocean University of China, 2017, 47(4): 81-92.

    Supported by the National Basic Research Program(2010CB951202)

    猜你喜歡
    調(diào)沙落潮漲潮
    迎戰(zhàn)調(diào)水調(diào)沙
    漲潮
    散文詩世界(2023年3期)2023-04-15 17:12:36
    2011—2020年調(diào)水調(diào)沙前后黃河口海域浮游動(dòng)物年間變化及影響因子
    潮 汐
    “水城”被淹
    中資出海并購“落潮”
    為什么會(huì)有漲潮和落潮
    神奇的潮汐
    漲潮啦
    長江口北槽落潮分流比變化原因分析
    亚洲国产日韩欧美精品在线观看 | 欧美精品亚洲一区二区| or卡值多少钱| 国产1区2区3区精品| 亚洲人成网站在线播放欧美日韩| 一个人观看的视频www高清免费观看 | 搡老妇女老女人老熟妇| 亚洲av成人av| 在线观看免费视频日本深夜| 国产成人欧美在线观看| 国产精品av久久久久免费| 国产亚洲精品第一综合不卡| 亚洲av片天天在线观看| 午夜免费激情av| 国产私拍福利视频在线观看| 国产精品综合久久久久久久免费| 看片在线看免费视频| 国产视频内射| 少妇裸体淫交视频免费看高清 | 久久婷婷人人爽人人干人人爱| 精品少妇一区二区三区视频日本电影| 在线看三级毛片| 午夜影院日韩av| netflix在线观看网站| 欧美极品一区二区三区四区| 丝袜人妻中文字幕| 午夜福利免费观看在线| 日韩免费av在线播放| 午夜久久久久精精品| 欧美在线一区亚洲| 亚洲成人中文字幕在线播放| 婷婷亚洲欧美| 天堂av国产一区二区熟女人妻 | 欧美日韩福利视频一区二区| 免费在线观看完整版高清| 国产视频一区二区在线看| 午夜精品一区二区三区免费看| 成年人黄色毛片网站| 18禁国产床啪视频网站| 国产精品综合久久久久久久免费| 91在线观看av| 亚洲第一欧美日韩一区二区三区| 两个人视频免费观看高清| 亚洲精品久久成人aⅴ小说| 女同久久另类99精品国产91| 国产aⅴ精品一区二区三区波| 青草久久国产| 丝袜美腿诱惑在线| 国产高清激情床上av| 一二三四在线观看免费中文在| 伦理电影免费视频| 黄色片一级片一级黄色片| 色播亚洲综合网| 岛国视频午夜一区免费看| 神马国产精品三级电影在线观看 | 亚洲成人精品中文字幕电影| 亚洲中文字幕日韩| 99精品在免费线老司机午夜| 欧美另类亚洲清纯唯美| 久久久久九九精品影院| 免费电影在线观看免费观看| 999久久久国产精品视频| 亚洲一区中文字幕在线| 国内久久婷婷六月综合欲色啪| 1024视频免费在线观看| 黄色视频不卡| 久久久国产欧美日韩av| 欧美日韩精品网址| 一级毛片女人18水好多| 亚洲人成伊人成综合网2020| 99riav亚洲国产免费| 亚洲性夜色夜夜综合| av欧美777| 国产精品99久久99久久久不卡| 久久久久久久久免费视频了| 日本撒尿小便嘘嘘汇集6| 国产黄色小视频在线观看| 国产欧美日韩一区二区三| 日韩免费av在线播放| 日本熟妇午夜| 国产麻豆成人av免费视频| 岛国在线免费视频观看| 亚洲熟妇中文字幕五十中出| 男人舔奶头视频| 亚洲国产精品合色在线| 亚洲欧美日韩高清在线视频| 国产视频内射| 俄罗斯特黄特色一大片| 国产三级中文精品| 精品国产超薄肉色丝袜足j| 婷婷亚洲欧美| 91在线观看av| 国产又黄又爽又无遮挡在线| 日韩有码中文字幕| 18美女黄网站色大片免费观看| 国产成人精品无人区| 色综合亚洲欧美另类图片| 久久国产乱子伦精品免费另类| 中文亚洲av片在线观看爽| 国产黄片美女视频| 熟妇人妻久久中文字幕3abv| 91在线观看av| 99精品在免费线老司机午夜| 蜜桃久久精品国产亚洲av| 国产精品 国内视频| 欧美黑人欧美精品刺激| 国产真人三级小视频在线观看| 一边摸一边做爽爽视频免费| 中文资源天堂在线| 制服人妻中文乱码| 99热6这里只有精品| 国产精华一区二区三区| 亚洲国产精品成人综合色| 亚洲欧美一区二区三区黑人| 日韩三级视频一区二区三区| 一进一出好大好爽视频| 日韩成人在线观看一区二区三区| 欧美一区二区精品小视频在线| 欧美日韩精品网址| 精品国产乱子伦一区二区三区| 蜜桃久久精品国产亚洲av| 最近在线观看免费完整版| 久久精品91无色码中文字幕| 午夜福利免费观看在线| 人人妻人人澡欧美一区二区| 成人特级黄色片久久久久久久| 一本大道久久a久久精品| 伦理电影免费视频| 成年女人毛片免费观看观看9| av福利片在线| 日韩高清综合在线| 亚洲aⅴ乱码一区二区在线播放 | 久久精品aⅴ一区二区三区四区| 成人国产一区最新在线观看| 女同久久另类99精品国产91| 国产精品亚洲美女久久久| 亚洲欧美日韩高清在线视频| 少妇的丰满在线观看| 最近在线观看免费完整版| av视频在线观看入口| 岛国视频午夜一区免费看| 人人妻人人澡欧美一区二区| 色老头精品视频在线观看| 人人妻人人看人人澡| 欧美日韩亚洲综合一区二区三区_| 97人妻精品一区二区三区麻豆| 国产一区二区三区视频了| 一级毛片高清免费大全| 黄色a级毛片大全视频| 熟女少妇亚洲综合色aaa.| 欧美色视频一区免费| 麻豆av在线久日| 亚洲一码二码三码区别大吗| 国产精品免费一区二区三区在线| 亚洲熟妇中文字幕五十中出| 国产高清激情床上av| 人人妻人人澡欧美一区二区| 天堂√8在线中文| 日韩三级视频一区二区三区| 亚洲国产高清在线一区二区三| 岛国在线观看网站| 欧美绝顶高潮抽搐喷水| 久久久久免费精品人妻一区二区| 一级毛片精品| www.熟女人妻精品国产| 亚洲天堂国产精品一区在线| 久久草成人影院| 中文字幕精品亚洲无线码一区| 美女 人体艺术 gogo| 俺也久久电影网| 叶爱在线成人免费视频播放| 在线永久观看黄色视频| 村上凉子中文字幕在线| 999精品在线视频| 麻豆成人av在线观看| 亚洲av五月六月丁香网| 一本综合久久免费| 一级毛片高清免费大全| 国产99久久九九免费精品| 色尼玛亚洲综合影院| 午夜福利成人在线免费观看| 国产精品亚洲av一区麻豆| 最近在线观看免费完整版| 国产精品久久视频播放| 国产区一区二久久| 香蕉丝袜av| 亚洲国产欧洲综合997久久,| 两个人免费观看高清视频| 国产成人系列免费观看| 国产av在哪里看| 国产成人系列免费观看| 久久久久亚洲av毛片大全| 欧美日韩黄片免| 熟女电影av网| 日本成人三级电影网站| 国产精品99久久99久久久不卡| 在线国产一区二区在线| 91老司机精品| 成人欧美大片| 亚洲av成人一区二区三| 亚洲午夜精品一区,二区,三区| 日韩精品中文字幕看吧| 亚洲国产精品合色在线| 成人特级黄色片久久久久久久| 精品国产超薄肉色丝袜足j| 国产精品亚洲一级av第二区| 亚洲午夜精品一区,二区,三区| 成年版毛片免费区| 老熟妇仑乱视频hdxx| 男人舔女人的私密视频| 伊人久久大香线蕉亚洲五| 亚洲 欧美 日韩 在线 免费| 在线观看免费视频日本深夜| 中文亚洲av片在线观看爽| 在线观看www视频免费| 一二三四在线观看免费中文在| 又黄又粗又硬又大视频| 国产亚洲精品av在线| 亚洲国产日韩欧美精品在线观看 | 又粗又爽又猛毛片免费看| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲性夜色夜夜综合| 欧美成人一区二区免费高清观看 | 午夜激情av网站| 亚洲国产精品合色在线| 全区人妻精品视频| www.自偷自拍.com| 香蕉av资源在线| 露出奶头的视频| 我要搜黄色片| 一级作爱视频免费观看| 日韩高清综合在线| 最近在线观看免费完整版| 婷婷丁香在线五月| 欧美中文日本在线观看视频| 首页视频小说图片口味搜索| 久久久久久久久免费视频了| 精品高清国产在线一区| 精品国产乱子伦一区二区三区| 亚洲专区字幕在线| 日韩欧美 国产精品| 色在线成人网| 老汉色av国产亚洲站长工具| 免费人成视频x8x8入口观看| 狂野欧美白嫩少妇大欣赏| 日本在线视频免费播放| 精品久久久久久久人妻蜜臀av| 国产亚洲av高清不卡| 国产在线观看jvid| 国产高清视频在线观看网站| 亚洲性夜色夜夜综合| 中国美女看黄片| 亚洲第一欧美日韩一区二区三区| 久久 成人 亚洲| 在线观看www视频免费| 国产精品,欧美在线| АⅤ资源中文在线天堂| 精品日产1卡2卡| av福利片在线观看| 欧美黄色片欧美黄色片| 女警被强在线播放| 又黄又爽又免费观看的视频| 国产精品免费一区二区三区在线| bbb黄色大片| 一本综合久久免费| 欧美午夜高清在线| 国产成人av教育| 不卡av一区二区三区| 在线永久观看黄色视频| 欧美人与性动交α欧美精品济南到| 真人做人爱边吃奶动态| 99re在线观看精品视频| 香蕉国产在线看| 国产亚洲精品av在线| 国产av一区二区精品久久| 午夜精品久久久久久毛片777| 国产伦在线观看视频一区| av欧美777| 99在线视频只有这里精品首页| 国产亚洲精品一区二区www| 在线视频色国产色| 岛国在线免费视频观看| 桃色一区二区三区在线观看| 日日摸夜夜添夜夜添小说| 免费一级毛片在线播放高清视频| 99久久99久久久精品蜜桃| 国产乱人伦免费视频| 身体一侧抽搐| 999精品在线视频| 在线观看免费视频日本深夜| 无限看片的www在线观看| 午夜精品一区二区三区免费看| 可以在线观看的亚洲视频| 国产亚洲精品一区二区www| 久久精品91蜜桃| 免费电影在线观看免费观看| 亚洲av成人av| 久久精品夜夜夜夜夜久久蜜豆 | 久久久久久久久免费视频了| 亚洲精品av麻豆狂野| 在线观看舔阴道视频| 老汉色∧v一级毛片| 少妇人妻一区二区三区视频| 亚洲片人在线观看| 久久久久久国产a免费观看| 女生性感内裤真人,穿戴方法视频| 亚洲狠狠婷婷综合久久图片| 久久天躁狠狠躁夜夜2o2o| 国内精品久久久久久久电影| 制服人妻中文乱码| 天堂av国产一区二区熟女人妻 | 两人在一起打扑克的视频| 久久热在线av| 久久人妻av系列| 久久久久久久久中文| 欧美精品啪啪一区二区三区| 亚洲精品一区av在线观看| 亚洲av成人不卡在线观看播放网| 亚洲美女黄片视频| 88av欧美| 高清毛片免费观看视频网站| 美女扒开内裤让男人捅视频| 男男h啪啪无遮挡| 欧美成人午夜精品| 国产主播在线观看一区二区| 精品国产乱码久久久久久男人| 午夜精品久久久久久毛片777| 最近最新中文字幕大全电影3| 国产精品99久久99久久久不卡| 中文字幕人成人乱码亚洲影| 51午夜福利影视在线观看| 国产亚洲欧美98| 午夜日韩欧美国产| 又粗又爽又猛毛片免费看| 亚洲精品美女久久久久99蜜臀| 日韩高清综合在线| 国产成人啪精品午夜网站| 欧美精品啪啪一区二区三区| 亚洲av美国av| 在线观看日韩欧美| 91老司机精品| 在线观看免费视频日本深夜| 国产精品香港三级国产av潘金莲| 日韩欧美免费精品| 国产精品美女特级片免费视频播放器 | 欧美+亚洲+日韩+国产| 亚洲中文av在线| 别揉我奶头~嗯~啊~动态视频| 久久久久国产一级毛片高清牌| 岛国在线观看网站| 18禁黄网站禁片免费观看直播| 一夜夜www| 亚洲成av人片免费观看| 国产日本99.免费观看| 啦啦啦观看免费观看视频高清| 女同久久另类99精品国产91| 亚洲av美国av| 亚洲精品久久国产高清桃花| 成人18禁在线播放| 国产精品一区二区精品视频观看| 1024香蕉在线观看| 91国产中文字幕| 91大片在线观看| 欧美三级亚洲精品| 欧美黄色片欧美黄色片| 欧美黑人巨大hd| 性欧美人与动物交配| 一区二区三区国产精品乱码| 国产一级毛片七仙女欲春2| 波多野结衣高清作品| 国产视频内射| 欧美精品啪啪一区二区三区| 久久亚洲精品不卡| 午夜福利视频1000在线观看| 成人欧美大片| 亚洲中文字幕日韩| 午夜福利高清视频| 国产精品免费一区二区三区在线| 亚洲av第一区精品v没综合| 欧美3d第一页| 欧美激情久久久久久爽电影| 女同久久另类99精品国产91| 美女高潮喷水抽搐中文字幕| 一级作爱视频免费观看| 久久精品国产综合久久久| 欧美黑人精品巨大| 久久性视频一级片| 老司机深夜福利视频在线观看| 91麻豆精品激情在线观看国产| 成年人黄色毛片网站| 啦啦啦韩国在线观看视频| 99久久久亚洲精品蜜臀av| 亚洲人成77777在线视频| 97超级碰碰碰精品色视频在线观看| av免费在线观看网站| www.999成人在线观看| 久久久久性生活片| av福利片在线观看| 黄色视频不卡| 在线观看日韩欧美| 色老头精品视频在线观看| 欧美黄色片欧美黄色片| 国内精品久久久久久久电影| 欧美激情久久久久久爽电影| 男女床上黄色一级片免费看| 国产麻豆成人av免费视频| 超碰成人久久| 久久99热这里只有精品18| 精品欧美国产一区二区三| 又黄又爽又免费观看的视频| 人人妻,人人澡人人爽秒播| 国产av不卡久久| 国内精品久久久久久久电影| 日本熟妇午夜| 日韩欧美一区二区三区在线观看| 99久久无色码亚洲精品果冻| xxx96com| 好男人电影高清在线观看| 在线观看美女被高潮喷水网站 | 十八禁网站免费在线| 中文字幕精品亚洲无线码一区| 欧美绝顶高潮抽搐喷水| 国产免费av片在线观看野外av| 1024手机看黄色片| 亚洲国产日韩欧美精品在线观看 | 99久久精品国产亚洲精品| 中文亚洲av片在线观看爽| 香蕉丝袜av| 国产精品电影一区二区三区| 99久久精品热视频| 国产97色在线日韩免费| 91字幕亚洲| 午夜福利视频1000在线观看| 亚洲国产精品合色在线| 亚洲狠狠婷婷综合久久图片| 亚洲国产欧美网| 19禁男女啪啪无遮挡网站| 久久久久久久久免费视频了| 国产av一区二区精品久久| 首页视频小说图片口味搜索| 老熟妇仑乱视频hdxx| 18禁美女被吸乳视频| 国产人伦9x9x在线观看| 亚洲乱码一区二区免费版| 中国美女看黄片| 男男h啪啪无遮挡| 国产主播在线观看一区二区| 香蕉国产在线看| 成人高潮视频无遮挡免费网站| 欧美一区二区国产精品久久精品 | 黄色毛片三级朝国网站| 18禁黄网站禁片免费观看直播| 黄色丝袜av网址大全| 蜜桃久久精品国产亚洲av| 天天躁夜夜躁狠狠躁躁| 久久天躁狠狠躁夜夜2o2o| 搡老熟女国产l中国老女人| 日韩欧美一区二区三区在线观看| 午夜精品久久久久久毛片777| 韩国av一区二区三区四区| 久久这里只有精品中国| 狂野欧美白嫩少妇大欣赏| 亚洲 欧美 日韩 在线 免费| 两个人视频免费观看高清| 中文字幕久久专区| 一a级毛片在线观看| 国产欧美日韩一区二区精品| 757午夜福利合集在线观看| 国产精品久久久人人做人人爽| 最近在线观看免费完整版| 欧美日韩精品网址| 91av网站免费观看| 丁香六月欧美| 国产高清激情床上av| 美女 人体艺术 gogo| 小说图片视频综合网站| 叶爱在线成人免费视频播放| 久久久久久久久免费视频了| 成人三级黄色视频| 国产av又大| 在线十欧美十亚洲十日本专区| 亚洲av成人av| 欧美日韩亚洲国产一区二区在线观看| 麻豆久久精品国产亚洲av| 精品久久久久久久末码| 国产熟女xx| 精品国产美女av久久久久小说| 久久久久久亚洲精品国产蜜桃av| 不卡一级毛片| 久久精品成人免费网站| 久久久久免费精品人妻一区二区| 国产高清激情床上av| 在线看三级毛片| 亚洲一区二区三区色噜噜| 久久亚洲真实| 国产精品电影一区二区三区| 久久婷婷人人爽人人干人人爱| 少妇的丰满在线观看| 又粗又爽又猛毛片免费看| 精品久久久久久久末码| 欧美日韩瑟瑟在线播放| 亚洲一码二码三码区别大吗| 欧美日韩精品网址| av超薄肉色丝袜交足视频| av免费在线观看网站| 国产高清有码在线观看视频 | 亚洲成人免费电影在线观看| 中亚洲国语对白在线视频| 99国产精品一区二区三区| 国产一区二区三区在线臀色熟女| 真人做人爱边吃奶动态| 国产午夜精品久久久久久| 中文字幕久久专区| 亚洲av美国av| 久久人妻av系列| 18禁国产床啪视频网站| 日本在线视频免费播放| av欧美777| 熟女少妇亚洲综合色aaa.| 国产真人三级小视频在线观看| 欧美中文日本在线观看视频| 国产男靠女视频免费网站| 黄色片一级片一级黄色片| 欧美一区二区精品小视频在线| 麻豆av在线久日| 精品久久久久久久人妻蜜臀av| 亚洲精品色激情综合| 成熟少妇高潮喷水视频| 亚洲熟女毛片儿| 狠狠狠狠99中文字幕| 国产亚洲精品综合一区在线观看 | 亚洲第一欧美日韩一区二区三区| 亚洲成人久久性| 男人的好看免费观看在线视频 | 91成年电影在线观看| 久久久久精品国产欧美久久久| 亚洲成人中文字幕在线播放| 亚洲免费av在线视频| 身体一侧抽搐| a级毛片在线看网站| 99国产精品99久久久久| 国产精品,欧美在线| 精品久久蜜臀av无| 日韩精品青青久久久久久| 午夜精品在线福利| 久久国产精品影院| 51午夜福利影视在线观看| 国产成年人精品一区二区| 三级毛片av免费| 每晚都被弄得嗷嗷叫到高潮| 又大又爽又粗| 国内精品一区二区在线观看| 最近视频中文字幕2019在线8| 国产在线精品亚洲第一网站| 日韩精品免费视频一区二区三区| 午夜激情福利司机影院| 一边摸一边做爽爽视频免费| 亚洲av熟女| 搞女人的毛片| 国产精品av久久久久免费| 国产aⅴ精品一区二区三区波| 久久这里只有精品中国| 特大巨黑吊av在线直播| 欧美又色又爽又黄视频| 国产真实乱freesex| av福利片在线观看| 国产成人精品无人区| 女人爽到高潮嗷嗷叫在线视频| 变态另类丝袜制服| 国产成人av激情在线播放| 久久精品aⅴ一区二区三区四区| 99久久综合精品五月天人人| 亚洲自拍偷在线| 两人在一起打扑克的视频| 欧美又色又爽又黄视频| 国产成人av教育| 天天躁夜夜躁狠狠躁躁| 亚洲人成电影免费在线| 成人av一区二区三区在线看| tocl精华| 黄色 视频免费看| 久久久国产成人免费| 久久婷婷成人综合色麻豆| 欧美日韩一级在线毛片| 欧美一区二区精品小视频在线| 99在线视频只有这里精品首页| 色哟哟哟哟哟哟| 最近在线观看免费完整版| 91字幕亚洲| 欧美3d第一页| 国产精品香港三级国产av潘金莲| www.精华液| 欧美黑人欧美精品刺激| 成人欧美大片| 中文字幕熟女人妻在线| 不卡av一区二区三区| av超薄肉色丝袜交足视频| 亚洲精品一卡2卡三卡4卡5卡| tocl精华| 午夜福利在线在线| 亚洲av电影不卡..在线观看| 日韩欧美三级三区| 精品久久久久久久人妻蜜臀av| 久久久国产精品麻豆| 欧美日韩国产亚洲二区| 国产亚洲精品久久久久5区| 啦啦啦韩国在线观看视频| 琪琪午夜伦伦电影理论片6080| 亚洲av片天天在线观看| 国产成人精品无人区| 亚洲专区国产一区二区| 俺也久久电影网| 国产一区二区三区视频了| 日本一区二区免费在线视频| 神马国产精品三级电影在线观看 |