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

    伶仃洋洪季懸沙平面分布特征及成因探討

    2017-09-05 11:28:53何嘉偉賈良文韋獻(xiàn)革賈妍紅程聰
    海洋學(xué)報(bào) 2017年9期
    關(guān)鍵詞:懸沙伶仃洋沙場(chǎng)

    何嘉偉,賈良文,韋獻(xiàn)革,賈妍紅,程聰

    (1.中山大學(xué) 海洋科學(xué)學(xué)院 近岸海洋科學(xué)與技術(shù)研究中心,廣東 廣州 510275; 2.中山大學(xué) 環(huán)境科學(xué)與工程學(xué)院,廣東 廣州 510275)

    伶仃洋洪季懸沙平面分布特征及成因探討

    何嘉偉1,賈良文1,韋獻(xiàn)革2*,賈妍紅1,程聰1

    (1.中山大學(xué) 海洋科學(xué)學(xué)院 近岸海洋科學(xué)與技術(shù)研究中心,廣東 廣州 510275; 2.中山大學(xué) 環(huán)境科學(xué)與工程學(xué)院,廣東 廣州 510275)

    基于伶仃洋河口2007年洪季(8月)潮流和泥沙現(xiàn)場(chǎng)觀測(cè)資料,通過(guò)實(shí)測(cè)資料分析結(jié)合二維泥沙數(shù)值模型的方法,分析其懸沙濃度平面分布特征。結(jié)果表明,懸沙濃度從伶仃洋灣頂虎門至灣口沿西槽向海在縱向上呈現(xiàn)出先減小再增大后減小的變化趨勢(shì);橫向上,西槽最大,東槽和中灘次之;內(nèi)伶仃島西北側(cè)海域?yàn)楹沉扛咧抵行?。懸沙縱向輸移機(jī)制分析表明,伶仃洋洪季輸沙主要貢獻(xiàn)項(xiàng)為平流向海輸沙、潮汐捕集和垂向環(huán)流的向陸輸沙,其中后兩項(xiàng)在內(nèi)伶仃島附近的貢獻(xiàn)較大。潮流、徑流來(lái)沙、徑潮強(qiáng)度對(duì)比以及地形邊界作用為懸沙場(chǎng)平面分布差異的主要成因。

    伶仃洋;懸沙濃度;平面分布特征;成因探討

    1 引言

    河口的發(fā)育主要依賴于徑流、潮流、波浪和泥沙等自然條件的影響,但近幾十年來(lái),我國(guó)河口地區(qū)的人類活動(dòng)愈發(fā)劇烈,甚至成為了河口發(fā)育演變的主導(dǎo)因素。關(guān)于在人類活動(dòng)影響下我國(guó)河口地區(qū)泥沙新時(shí)期特征的研究也逐漸增多,Dai等[1]認(rèn)為長(zhǎng)江三峽大壩工程建成后,雖然上游輸沙量減少,但河口攔門沙附近區(qū)域的年平均懸沙濃度變化不大,攔門沙內(nèi)側(cè)懸沙濃度的季節(jié)變化為洪高枯低,攔門沙外側(cè)則洪低枯高。Kong等[2]指出黃河下游在小浪底水庫(kù)水沙調(diào)控影響下,輸沙量總體相對(duì)穩(wěn)定,其中洪季有所下降,枯季有所增加,泥沙中值粒徑表現(xiàn)為上游細(xì)下游粗。胡德禮等[3]以1978年和1999年洪季作為對(duì)比,研究了人類活動(dòng)影響下珠江口門水沙分配調(diào)整后伶仃洋的水沙場(chǎng)的響應(yīng),結(jié)果表明,水動(dòng)力和懸沙濃度整體上均有增大。

    伶仃洋位于珠江三角洲的東南部,為呈南北走向的典型喇叭狀河口灣(圖1),東、西、北江的來(lái)水來(lái)沙通過(guò)東四口門即虎門、蕉門、洪奇瀝以及橫門注入伶仃洋灣內(nèi)。珠江三角洲注入伶仃洋的徑流和泥沙分別從20世紀(jì)80年代的占比53%和47%,到90年代增大至64%和57%[4]。伶仃洋灣內(nèi)和灣口存在天然的島嶼屏障,其水下地形的整體格局呈現(xiàn)為“三灘夾兩槽”,即東灘、中灘、西灘和東槽(礬石水道)、西槽(伶仃航道)。

    自20世紀(jì)80年代起,珠江三角洲的社會(huì)經(jīng)濟(jì)發(fā)展帶動(dòng)了河道采砂業(yè)的繁榮,而大規(guī)模且無(wú)序的采砂活動(dòng)使部分河床發(fā)生非自然深切,不同口門河段的水沙對(duì)比情況則發(fā)生了相應(yīng)調(diào)整。伶仃洋的分流比和分沙比逐步增大,同時(shí)出現(xiàn)大規(guī)模的圍墾活動(dòng),近30年來(lái)伶仃洋海域?yàn)┩棵娣e減少約132 km2,而且海岸線呈現(xiàn)出明顯的海進(jìn)趨勢(shì),人類活動(dòng)的影響幾乎掩蓋了其自然演變的趨勢(shì)[5]。西槽作為伶仃洋的出海主航道,通過(guò)航道疏浚工程,在2000年時(shí)已經(jīng)從平均水深6.8 m被挖深至11.5 m。故近30年來(lái),伶仃洋河口灣在河床采砂、灘涂圍墾、航道疏浚等強(qiáng)烈人類活動(dòng)影響下,其來(lái)水來(lái)沙以及邊界條件快速變化。懸沙分布及成因是伶仃洋河口灣開(kāi)發(fā)保護(hù)的基礎(chǔ)性研究工作,因此,需要在新的來(lái)水來(lái)沙、地形邊界條件下對(duì)其進(jìn)行研究,可為伶仃洋的綜合整治規(guī)劃、出海航道的維護(hù)等提供一定的科學(xué)依據(jù)。

    圖1 伶仃洋地形示意圖Fig.1 Topography of the Lingdingyang Estuary

    2 資料和方法

    2.1 資料來(lái)源

    水文資料觀測(cè)時(shí)間為2007年8月13-14日(大潮期間)以及16-17日(中潮期間),共布設(shè)了21條垂線進(jìn)行同步水文信息(潮流、含沙量和鹽度)測(cè)量以及13個(gè)水位站的同步潮位觀測(cè)。本文主要應(yīng)用其中14條垂線的潮流、含沙量以及8個(gè)水位站的同步潮位實(shí)測(cè)數(shù)據(jù)進(jìn)行分析以及模型驗(yàn)證(站點(diǎn)分布見(jiàn)圖2)。

    模型地形資料為2008年伶仃洋河口灣以及三角洲河網(wǎng)的水深數(shù)據(jù),采用北京54坐標(biāo)系,水深數(shù)據(jù)為珠江基面。

    圖2 伶仃洋洪季實(shí)測(cè)站點(diǎn)位置圖Fig.2 Observation stations in the Lingdingyang Estuary

    2.2 研究方法

    本文主要采用實(shí)測(cè)資料分析與TELEMAC數(shù)值模擬相結(jié)合的研究方法,分析河口灣內(nèi)洪季懸沙場(chǎng)的平面分布特征,并探討懸沙場(chǎng)的平面分布差異的成因。

    實(shí)測(cè)資料分析包括對(duì)實(shí)測(cè)流速、含沙量、鹽度等信息進(jìn)行數(shù)據(jù)處理以及作圖,如計(jì)算余流、滯流點(diǎn)和泥沙縱向輸移通量機(jī)制分解等。數(shù)值模擬方面,應(yīng)用TELEMAC水動(dòng)力模塊以及泥沙模塊SISYPHE模擬伶仃洋河口灣的水動(dòng)力場(chǎng)以及懸沙場(chǎng)。

    TELEMAC-MASCARET是一套用于模擬自由表面流的水動(dòng)力、水質(zhì)、沉積等的模型系統(tǒng),其由法國(guó)國(guó)家水力學(xué)和環(huán)境實(shí)驗(yàn)室開(kāi)發(fā)。模型采用非結(jié)構(gòu)三角形網(wǎng)格,數(shù)值模擬計(jì)算基于有限元法,穩(wěn)定性相對(duì)較好。

    伶仃洋區(qū)域泥沙大多數(shù)屬于黏性砂,關(guān)于水體與底床泥沙的交換,本文在泥沙模塊SISYPHE中采用Einstein和Krone[6]、Partheniades和Paaswell[7]的黏性砂經(jīng)典公式進(jìn)行侵蝕、沉積通量計(jì)算:

    (1)

    式中,E為河床泥沙侵蝕通量;D為河床泥沙沉積通量;M為系數(shù);u*為水流的摩阻流速;u*e和u*d分別為泥沙侵蝕和沉積的臨界剪切流速;Ws為泥沙沉降速度;C為懸沙濃度。當(dāng)水流的摩阻流速大于泥沙侵蝕的臨界剪切流速時(shí),河床發(fā)生侵蝕,河床泥沙進(jìn)入水體;當(dāng)水流的摩阻流速小于泥沙沉積的臨界剪切流速時(shí),懸移質(zhì)沉降至河床。

    模型中采用的二維懸移質(zhì)輸沙方程如下:

    (2)

    式中,h為水深;U和V分別為x、y方向的垂向平均流速;C為垂向平均懸沙濃度;εs為泥沙擴(kuò)散系數(shù);E-D為侵蝕量減去沉積量的凈泥沙通量。

    3 實(shí)測(cè)資料分析

    3.1 潮流分析

    3.1.1 縱向流速變化

    由灣頂至灣口沿西槽向海選取L1、L5、L7、L8、L11及L13站的大潮期間漲、落潮平均流速進(jìn)行比較,如圖3,伶仃洋海域潮流動(dòng)力從外海到灣中部增強(qiáng),灣中部過(guò)后開(kāi)始減弱,接近灣頂時(shí)潮流動(dòng)力再度增強(qiáng)。漲潮時(shí),潮流從外海進(jìn)入河口灣,受到地形橫斷面收縮的作用而增大,至L8站附近時(shí)受到的內(nèi)伶仃島的阻隔作用,潮流分支進(jìn)入西槽和東槽,故流速減小,在灣上部L5站,西、東槽漲潮流匯合進(jìn)入虎門,同時(shí)橫斷面收窄,故流速再次增強(qiáng)。落潮時(shí),潮流出虎門后橫斷面放寬,并分成兩支沿西、東槽下泄,故流速減小,在西槽L7站受到灣上部西部口門下泄流匯入的影響而流速增大,隨著河口橫斷面放寬和島嶼的阻流效應(yīng),潮流開(kāi)始不斷減小,最后流入外海。

    圖3 大潮期間流速縱向變化圖Fig.3 Longitudinal variation of current speed in the spring tide

    3.1.2 余流

    潮流余流是指從實(shí)測(cè)的潮流矢量疊加后剩下的部分,其對(duì)泥沙凈輸移以及污染物擴(kuò)散方向具有一定的指示意義。

    根據(jù)本次大、中潮的水文實(shí)測(cè)數(shù)據(jù)計(jì)算分析,各站點(diǎn)的表、底層余流情況如圖4所示。河口灣上部?jī)?nèi)伶仃島以北區(qū)域的余流基本指向海,偏南向,其中洪奇瀝出口L3站表層余流流速能達(dá)到0.48 m/s;L8站出現(xiàn)表層指向海,底層指向陸的典型密度環(huán)流結(jié)構(gòu);灣口大濠島西側(cè)海域的L11和L12站,底層余流較表層大,指向外海;中潮期間灣口外海的L13和L14站存在相對(duì)較強(qiáng)的向西的表層余流,主要因?yàn)成喜靠陂T淡水徑流沿著河口灣西側(cè)而下形成。

    3.1.3 滯流點(diǎn)

    滯流點(diǎn)概念最早見(jiàn)于Schultz和Simmons[8]以及Simmons和Broun[9]對(duì)水流優(yōu)勢(shì)程度理論的論述中。在一個(gè)潮周期內(nèi),由于河口存在密度梯度,沿程向陸的底部水流會(huì)從上溯流逐漸變成下泄流,中途會(huì)出現(xiàn)某處凈泄量為0的地點(diǎn),稱為“滯流點(diǎn)”。而滯流點(diǎn)所在的位置并不是固定的,它隨著徑流量與潮流的對(duì)比變化而在一定的范圍內(nèi)有規(guī)律地來(lái)回移動(dòng),而滯流帶通常為河口內(nèi)泥沙顯著淤積的區(qū)域。

    由灣頂至灣口沿西槽向海選取L1、L5、L7、L8及L11站作為縱向代表站位來(lái)繪制優(yōu)勢(shì)流曲線圖5。大潮期間,底層潮流滯流點(diǎn)位于L8站附近,而中潮期間底層潮流滯流點(diǎn)位置則往上游移動(dòng)至L7站附近。表明在L7站和L8站之間,即內(nèi)伶仃島附近為泥沙聚集沉降的區(qū)域。

    3.2 懸沙分析

    3.2.1 縱向懸沙濃度變化

    由灣頂至灣口沿西槽向海選取L1、L5、L7、L8、L11以及L13站作為代表站比較伶仃洋內(nèi)縱向平均懸沙濃度的變化。如圖6所示,河口灣中部L7站和L8站之間的懸沙濃度最大,恰好是潮流分析中滯流點(diǎn)的位置,懸沙濃度往外海迅速減小,往灣頂方向先減小,到灣上部口門后增大。

    圖4 各站實(shí)測(cè)表、底層余流圖Fig.4 Surface and bottom residual currents at the observation stations

    圖5 優(yōu)勢(shì)流曲線圖Fig.5 Dominant current curves

    圖6 懸沙濃度縱向變化圖Fig.6 Longitudinal variation of suspended sediment concentration

    3.2.2 懸沙縱向輸移分解

    根據(jù)單寬懸沙縱向輸運(yùn)模式[10-11]將懸沙輸移量分解為8項(xiàng):

    (3)

    上述具體算法以及公式中各項(xiàng)參數(shù)意義見(jiàn)文獻(xiàn)[12]。其中,T1和T2共為平流項(xiàng)輸移;T3和T4分別表示潮汐、潮流與懸沙濃度的潮變化項(xiàng);T5表示垂向流速與懸沙濃度變化相關(guān)項(xiàng);T6、T7和T8項(xiàng)的貢獻(xiàn)一般較小不予考慮。計(jì)算時(shí),灣上部口門河段選取的流速值投影到河道方向,其他站的流速值則投影到南北方向,負(fù)號(hào)表示向陸輸移,正號(hào)表示向海輸移。

    計(jì)算結(jié)果如表1所示,河口灣中上部口門站點(diǎn)基本為向海輸沙為主,以向海平流輸沙項(xiàng)T1和T2的貢獻(xiàn)為主,但在內(nèi)伶仃島東西側(cè)附近海域的L7、L9、L6和L10站的T4項(xiàng)和T5項(xiàng)的貢獻(xiàn)有所增大。位于內(nèi)伶仃島南側(cè)附近海域的L8站,T4和T5項(xiàng)的向陸輸沙貢獻(xiàn)接近向海平流輸沙項(xiàng)的貢獻(xiàn)。其中T4項(xiàng)體現(xiàn)河口近底過(guò)程對(duì)最大渾濁帶高濃度懸沙的貢獻(xiàn),如床沙的再懸浮過(guò)程和潮流不對(duì)稱引起的輸沙[13],而T5項(xiàng)則表示垂向環(huán)流對(duì)輸沙的貢獻(xiàn)。L8站對(duì)應(yīng)潮流分析中滯流點(diǎn)和垂向環(huán)流發(fā)育的位置。

    表1 大潮期間懸沙輸移分解項(xiàng)統(tǒng)計(jì)表

    從各輸沙因子的重要性來(lái)說(shuō),對(duì)伶仃洋河口灣洪季輸沙影響最大的為T1、T2、T4和T5項(xiàng),主要為平流向海輸沙,潮汐捕集和垂向環(huán)流的向陸輸沙,其他項(xiàng)的貢獻(xiàn)相對(duì)十分小。

    4 模型建立與驗(yàn)證

    4.1 模型建立

    本研究基于TELEMAC系統(tǒng)內(nèi)的水動(dòng)力模塊TELEMAC-2D和泥沙模塊SISYPHE建立了兩個(gè)模型,一為包含珠江網(wǎng)河、河口灣及相鄰近海的大范圍水動(dòng)力模型 Ⅰ,二為伶仃洋河口灣水沙耦合模型 Ⅱ。

    模型Ⅰ計(jì)算區(qū)域如圖7a所示,網(wǎng)格單元數(shù)共有202 321個(gè),節(jié)點(diǎn)個(gè)數(shù)為114 665,外海單元空間步長(zhǎng)最大為7 000 m,部分上游內(nèi)河段網(wǎng)格空間步長(zhǎng)達(dá)50 m。模型上游邊界從東向西分別到達(dá)了東江博羅、流溪河老鴉崗、北江石角、西江高要以及潭江石咀水文(位)站,輸入實(shí)測(cè)流量或潮位,外海邊界采用OSU預(yù)報(bào)的潮位數(shù)據(jù),并應(yīng)用大萬(wàn)山T9站實(shí)測(cè)水位進(jìn)行校正。

    模型Ⅱ計(jì)算區(qū)域?yàn)榱尕暄蠛涌跒?圖7b),網(wǎng)格單元數(shù)共有44 507個(gè),節(jié)點(diǎn)個(gè)數(shù)為23 114,網(wǎng)格空間步長(zhǎng)200~500 m。泥沙參數(shù)中值粒徑給定為0.002 cm,沉降速度給定為0.05 cm/s,各開(kāi)邊界懸沙輸入取本次實(shí)測(cè)懸沙濃度的平均值。

    4.2 模型驗(yàn)證

    采用實(shí)測(cè)數(shù)據(jù)分別對(duì)大范圍水動(dòng)力模型Ⅰ結(jié)果的潮位(9個(gè)站點(diǎn))、潮流(14個(gè)站點(diǎn))以及伶仃洋水沙耦合模型Ⅱ結(jié)果的懸沙濃度(8個(gè)站點(diǎn))進(jìn)行驗(yàn)證。其中驗(yàn)證點(diǎn)位置見(jiàn)圖2,大潮驗(yàn)證時(shí)間為2007年8月13-14日,中潮潮驗(yàn)證時(shí)間為16-17日。

    潮位驗(yàn)證結(jié)果如圖8所示,RMSE平均為0.137 m,相關(guān)系數(shù)r為0.988;流速驗(yàn)證結(jié)果如圖9所示,RMSE平均為0.101 m/s,相關(guān)系數(shù)r為0.967;懸沙濃度驗(yàn)證結(jié)果如圖10所示,RMSE平均為0.051 kg/m3。總體來(lái)說(shuō),模擬結(jié)果較好,模型基本能反映伶仃洋洪季水動(dòng)力特點(diǎn)和懸沙變化過(guò)程。

    圖7 模型計(jì)算范圍Fig.7 TELEMAC model computational domains

    圖8 潮位驗(yàn)證圖Fig.8 Comparison between the model calculated and observed water level

    圖9 潮流流速驗(yàn)證圖Fig.9 Comparison between the model calculated and observed tidal current velocity

    圖10 懸沙濃度驗(yàn)證圖Fig.10 Comparison between the model calculated and observed suspended sediment concentration

    圖11 大、小潮平均懸沙場(chǎng)平面分布Fig.11 Disribution of suspended sediment concentration in spring and neap tide

    圖12 大、小潮各斷面懸沙場(chǎng)變化圖Fig.12 Longitudinal variation of suspended sediment concentration in the spring and neap tide

    圖13 伶仃洋洪季各特征時(shí)刻懸沙場(chǎng)Fig.13 Distribution of suspended sediment concentration at characteristic time in the Lingdingyang Estuary

    5 懸沙場(chǎng)平面分布

    5.1 平均懸沙場(chǎng)

    伶仃洋河口灣洪季大潮(13-14日)以及小潮期間(20-21日)的平均懸沙場(chǎng)平面分布模擬結(jié)果如圖11所示。大潮期間,懸沙濃度最高的區(qū)域位于內(nèi)伶仃島西北側(cè)海域,約為0.20 kg/m3,其次在東槽以及中灘西北部區(qū)域的也相對(duì)較大,達(dá)到0.12 kg/m3以上,而虎門河道、川鼻水道以及燈籠水道為0.08 kg/m3,其余地方的懸沙濃度則相對(duì)較小。小潮期間整體平均懸沙濃度水平明顯降低,內(nèi)伶仃島西北側(cè)的海域含沙量?jī)H約0.05 kg/m3,其中虎門段含沙量較高可能是受到模型邊界懸沙輸入的影響。

    在大潮平均懸沙場(chǎng)中取2個(gè)縱斷面和4個(gè)橫斷面進(jìn)行分析(圖11a),其懸沙濃度的沿程變化如圖12所示。沿D1縱斷面,灣頂虎門處含沙量達(dá)0.10 kg/m3,隨后沿西槽向海迅速增大,至內(nèi)伶仃島以北接近0.20 kg/m3,在內(nèi)伶仃島以南開(kāi)始逐漸減?。谎貣|槽的縱斷面D2,從龍穴島附近向海含沙量變化不大,而過(guò)了內(nèi)伶仃島之后含沙量則逐漸下降。灣上部橫斷面D3則表現(xiàn)為中部礬石水道較高;灣中部的橫斷面D4以及D5,從西往東,西灘水體含沙量較小,在西槽達(dá)最大值約0.20 kg/m3,至中灘懸沙濃度略減,而后至東槽又略增,然后向東灘減小;位于內(nèi)伶仃島南側(cè)和大濠島北側(cè)的斷面D6,從西灘往東含沙量逐漸增大,西槽處達(dá)最大,隨后在內(nèi)伶仃洋島南側(cè)海域出現(xiàn)急劇下降,再往東緩慢下降。

    從各斷面的懸沙濃度變化情況可知,縱向上出虎門口后懸沙濃度先減小后增大,在內(nèi)伶仃島附近懸沙濃度最大,隨后向海逐漸減小。橫向上看,西槽最大,東槽和中灘次之。

    5.2 各特征時(shí)刻懸沙場(chǎng)

    伶仃洋河口灣洪季大潮以及小潮的漲急時(shí)刻、漲憩時(shí)刻、落急時(shí)刻和落憩時(shí)刻懸沙場(chǎng)分布模擬結(jié)果如圖13所示。

    大潮落憩時(shí)刻伶仃洋整體的懸沙濃度最大,內(nèi)伶仃島西側(cè)海域能達(dá)到0.40 kg/m3,平面上表現(xiàn)為西、東槽以及中灘的含沙量相對(duì)較大。落急時(shí)刻整體區(qū)域的懸沙濃度也比較大,但是比落憩時(shí)刻略小,高值區(qū)域可達(dá)0.30 kg/m3,亦出現(xiàn)在內(nèi)伶仃島西北側(cè)的海域。漲急和漲憩時(shí)刻整體含沙量明顯比落急和落憩時(shí)刻小,含沙量高值區(qū)域出現(xiàn)在內(nèi)伶仃島西北側(cè)的海域,但不到0.20 kg/m3。小潮期間懸沙平面分布相對(duì)情況與大潮期間類似,但量值明顯減小,最大值出現(xiàn)在落急時(shí)刻,約0.12 kg/m3。

    6 懸沙場(chǎng)分布成因探討

    6.1 潮流作用

    一般來(lái)說(shuō),較大的流速條件下會(huì)產(chǎn)生相對(duì)大的底切應(yīng)力,從而能在底床起動(dòng)更多的泥沙,增大水體的懸沙濃度。灣內(nèi)漲、落潮流不對(duì)稱現(xiàn)象明顯,大、小潮潮流強(qiáng)度變化亦明顯,這必然影響懸沙的潮內(nèi)變化情況。選取西槽上、中、下游3個(gè)站點(diǎn)作為代表(位置見(jiàn)圖11a),其流速與懸沙濃度過(guò)程線如圖14所示,由圖可見(jiàn),該站懸沙濃度和流速大小變化趨勢(shì)較一致并具有周期性變化,特別是當(dāng)流速在落急時(shí)候達(dá)到極大值時(shí),懸沙濃度隨后迅速增大達(dá)到極大值,當(dāng)流速下降時(shí),懸沙濃度亦有一個(gè)相對(duì)緩慢的下降過(guò)程,懸沙濃度對(duì)流速的響應(yīng)存在滯后現(xiàn)象,延遲約2 h。此外,大潮期間流速極值小于小潮期間,大潮期間的懸沙濃度亦相對(duì)大于小潮期間。故潮流動(dòng)力是伶仃洋泥沙再懸浮的主要?jiǎng)恿Α?/p>

    圖14 流速與懸沙濃度過(guò)程線Fig.14 Time series variation of velocity and suspended sediment concentration

    圖15 有無(wú)徑流來(lái)沙的大、小潮平均懸沙場(chǎng)Fig.15 Comparison of the suspended sediment concentration distribution without and with upstream sediment transport

    6.2 徑流來(lái)沙作用

    伶仃洋河口灣的來(lái)水來(lái)沙集中在洪季,其上游的西江高要站、北江石角站、東江博羅站8月份的月徑流量分別約為240億m3、30億m3和40億m3,各口門來(lái)水來(lái)水沙存在差異,故洪季的徑流來(lái)沙,以及灣上部東西來(lái)沙不平衡必將影響懸沙場(chǎng)的空間分布。

    為了探討徑流來(lái)沙對(duì)懸沙場(chǎng)的影響,本文在模型Ⅱ中關(guān)閉河口灣上游口門的懸沙輸入,計(jì)算了一個(gè)對(duì)比案例,并選取了大潮、小潮平均的狀態(tài)下,對(duì)有上游來(lái)沙和無(wú)上游來(lái)沙的懸沙分布進(jìn)行了比較。如圖15所示,大潮平均懸沙場(chǎng),在無(wú)上游來(lái)沙的情況下,僅內(nèi)伶仃島西北側(cè)以及東槽的懸沙濃度比較大,最大不超過(guò)0.15 kg/m3,當(dāng)有上游來(lái)沙時(shí),內(nèi)伶仃島西北側(cè)的高值范圍明顯往四周擴(kuò)散,懸沙濃度最大能達(dá)到0.20 kg/m3,同時(shí)東槽的含沙量高值也不再局限于窄的深槽。其次,受到虎門河口徑流來(lái)沙的影響,中灘北部含沙量明顯增大。另外,各口門河段深槽的含沙量亦有明顯增大,反映了泥沙隨著徑流沿河道下泄入注灣內(nèi)的路徑。小潮期間的情況與大潮類似。

    6.3 徑、潮流對(duì)比作用

    伶仃洋河口灣存在四大口門,其中東側(cè)的虎門屬于潮流優(yōu)勢(shì)型,山潮比值僅為0.25,而西側(cè)的蕉門、洪奇門以及橫門則屬?gòu)搅鲀?yōu)勢(shì)型,山潮比分別為1.67、2.06以及2.63。

    根據(jù)3.1.3的實(shí)測(cè)資料分析得出,滯流點(diǎn)在內(nèi)伶仃島南北附近移動(dòng),而滯流點(diǎn)的位置是下泄流與上溯流的對(duì)比情況決定的,徑流會(huì)加大下泄流,滯流點(diǎn)位置是徑流、潮流作用的平衡點(diǎn),泥沙易于在此淤積,為懸沙活躍的中心。滯流點(diǎn)上游以下泄流作用為主,下游以上流溯流作用為主。根據(jù)3.2.2泥沙縱向輸移機(jī)制分析的結(jié)果,灣上部各站總體以徑流的平流凈向海輸沙為主,而滯流點(diǎn)以南的內(nèi)伶仃島南側(cè)L8站,潮流作用引起的T4項(xiàng)以及垂向環(huán)流作用引起的T5項(xiàng)的向陸輸沙貢獻(xiàn)較大?;㈤T為灣內(nèi)流速高值區(qū),其以平流向海凈輸沙貢獻(xiàn)為主,泥沙不在此淤積集中,故并未在此形成懸沙高值中心。

    6.4 地形邊界作用

    6.4.1 口門

    “門”是珠江三角洲獨(dú)特而典型的地貌單元,指的是通過(guò)狹窄的基巖峽口流入河口灣的口門?;㈤T兩側(cè)基巖島丘使其成為的峽口,而其上下游的獅子洋和伶仃洋河口灣的橫斷面均有放寬,并具有較大的納潮量,故在此發(fā)育雙向射流系統(tǒng),這特別有利于徑流攜帶大量泥沙往下游輸移[14]。而蕉門、洪奇瀝和橫門河道出口后橫斷面亦大幅放寬,導(dǎo)致流速逐漸減小,水流挾沙力下降,從而在出口門處形成局部懸沙濃度高值區(qū)。

    6.4.2 島嶼

    伶仃洋河口灣灣口分布大大小小的島嶼,如大濠島、萬(wàn)山群島等,這些灣口的島嶼阻擋外海波浪的傳入,使得灣內(nèi)的波浪作用減弱,懸沙場(chǎng)的分布受到波浪的改造作用也比較弱。灣口東側(cè)由于受到大濠島的影響,通過(guò)香港暗士頓水道出口的路徑幾乎與灣內(nèi)漲落潮主流向垂直,故懸沙更加集中分布在伶仃洋西側(cè)即靠近珠海-澳門一側(cè),但灣南部東側(cè)海域含沙量處于較低的狀態(tài)。位于灣中間的內(nèi)伶仃島,具有阻隔水、沙的作用,由于泥沙主要來(lái)源于東四口門,受到內(nèi)伶仃島的阻隔,東四口門下泄的泥沙更易在其北側(cè)淤積,這也是內(nèi)伶仃島北側(cè)形成中灘的原因之一。

    6.4.3 深槽

    伶仃洋河口灣內(nèi)水深較大的水道有:虎門口的川鼻水道、伶仃水道(西槽)、礬石水道(東槽)以及洪奇瀝出口的燈籠水道。特別是虎門河道,其上游獅子洋具備較大納潮量,當(dāng)處于大潮落潮時(shí),水流沿著川鼻水道深槽迅速下泄,隨后分成西、東兩路分別順著西槽和東槽攜帶大量泥沙快速下泄。深槽是漲落潮流的主要通道,具有較大的流速,同時(shí)也能更多地?cái)y帶泥沙下泄,導(dǎo)致其懸沙濃度增大。

    6.5 分、匯流作用

    鳧洲水道為蕉門的分汊,其出口幾乎以垂直方向匯入虎門水道,在交匯處懸沙在此形成強(qiáng)烈的混合聚集,所以在鳧洲水道東側(cè)出口處小范圍內(nèi)形成懸沙濃度高值中心(圖13);由于橫門出口分東西二汊,東汊呈東西向,因此,洪奇瀝、橫門出口處存在較強(qiáng)的橫向向東的輸沙流,其與來(lái)自蕉門的南向落潮流以幾乎以垂直方向交匯,懸沙在此形成強(qiáng)烈的混合聚集,在萬(wàn)頃沙南面海域形成一含沙量高值中心(圖13);淇澳島東側(cè),燈籠水道攜帶泥沙向東南匯入西槽,故在西槽中段高值中心的西側(cè),其含沙量亦同樣保持在較高值。挾沙水流的分、匯作用是灣內(nèi)局部區(qū)域懸沙濃度增強(qiáng)的原因之一。

    7 結(jié)論

    通過(guò)實(shí)測(cè)資料分析與數(shù)值模擬相結(jié)合,本文探討了伶仃洋河口灣洪季懸沙場(chǎng)平面分布特征及成因,主要結(jié)論如下:

    (1)實(shí)測(cè)水沙資料分析表明:洪季滯流點(diǎn)位于內(nèi)伶仃島附近;縱向輸移機(jī)制分析表明,伶仃洋洪季輸沙主要貢獻(xiàn)項(xiàng)為平流向海輸沙,潮汐捕集和垂向環(huán)流的向陸輸沙,其中后兩者在內(nèi)伶仃島附近的貢獻(xiàn)較大。

    (2)伶仃洋洪季懸沙場(chǎng)平面分布:縱向上,虎門出口后懸沙濃度先減小后增大,在內(nèi)伶仃島附近的懸沙濃度最大,隨后向海逐漸減??;橫向上看,西槽最大,東槽和中灘次之;高值中心出現(xiàn)在內(nèi)伶仃島西北側(cè)的海域。

    (3)從潮流、徑流來(lái)沙、徑潮強(qiáng)度對(duì)比、地形邊界作用分析了伶仃洋洪季懸沙場(chǎng)的平面分布差異成因:大潮落急時(shí)刻的高流速是泥沙再懸浮的主要?jiǎng)恿Γ粡搅鱽?lái)沙明顯會(huì)增大整個(gè)灣內(nèi)懸沙含量,而灣上部東西輸沙量不平衡則導(dǎo)致西部懸沙濃度高值和范圍比東部大;徑、潮流強(qiáng)度的對(duì)比變化在內(nèi)伶仃島附近形成滯流點(diǎn),利于形成懸沙濃度高值區(qū);地形邊界、分匯流作用是灣內(nèi)局部區(qū)域懸沙濃度增大或減小的原因之一。

    [1] Dai Zhijun, Chu Ao, Li Weihua, et al. Has suspended sediment concentration near the mouth bar of the Yangtze (Changjiang) Estuary been declining in recent years?[J]. Journal of Coastal Research, 2013, 29(4): 809-818.

    [2] Kong Dongxian, Miao Chiyuan, Wu Jingwen, et al. The hydro-environmental response on the lower Yellow River to the water-sediment regulation scheme[J]. Ecological Engineering, 2015, 79: 69-79.

    [3] 胡德禮, 楊清書(shū), 吳超羽, 等. 珠江網(wǎng)河水沙分配變化及其對(duì)伶仃洋水沙場(chǎng)的影響[J]. 水科學(xué)進(jìn)展, 2010, 21(1): 69-76.

    Hu Deli, Yang Qingshu, Wu Chaoyu, et al. Changing water and sediment dynamics in the Pearl River network and consequences on water and sediment regimes in the Lingdingyang estury[J]. Advances in Water Science, 2010, 21(1): 69-76.

    [4] 廣東省水利電力勘測(cè)設(shè)計(jì)研究院. 珠江河口近期防洪實(shí)施工程可行性研究報(bào)告[R]. 廣州: 廣東省水利電力勘測(cè)設(shè)計(jì)研究院, 2002: 20-35.

    Guangdong Hydropower Planning & Design Institute. Feasibility study report of recent flood control project implementation in Pearl River Estuary[R]. Guangzhou: Guangdong Hydropower Planning & Design Institute, 2002: 20-35.

    [5] 姚才華, 吳自銀. 30a來(lái)伶仃洋海岸線變遷及海底沖淤變化[J]. 海洋學(xué)研究, 2012, 30(3): 44-55.

    Yao Caihua,Wu Ziyin. Shoreline migration and scour-and-silting alternation in the Lingdingyang estuary over past 30 years[J]. Journal of Marine Sciences, 2012, 30(3): 44-55.

    [6] Einstein H A, Krone R B. Experiments to determine modes of cohesive sediment transport in salt water[J]. Journal of Geophysical Research, 1962, 67(4): 1451-1461.

    [7] Partheniades E, Paaswell R E. Erodibility of channels with cohesive boundary[J]. Journal of the Hydraulics Division, 1970, 96(3): 755-771.

    [8] Schultz E A, Simmons H B. Fresh water-salt water density currents: a major cause of siltation in estuaries[M]. Waterways Experiment Station, 1957.

    [9] Simmons H B, Broun F R. Salinity effect on hydraulics and shoaling in estuary[J]. IAHR 13th congress 1969,1969:311-326.

    [10] Ingram R G. Characteristics of the Great Whale River plume[J]. Journal of Geophysical Research:Oceans, 1981, 86(C3): 2017-2023.

    [11] Uncles R J, Elliott R C A, Weston S A. Observed fluxes of water, salt and suspended sediment in a partly mixed estuary[J]. Estuarine, Coastal and Shelf Science, 1985, 20(2): 147-167.

    [12] 任杰, 包蕓, 林衛(wèi)強(qiáng). 珠江口伶仃洋水沙縱向輸移特征分析[J]. 熱帶海洋學(xué)報(bào), 2001, 20(3): 35-40.

    Ren Jie, Bao Yun, Lin Weiqiang. Analyses on water and suspended sediment fluxes in Lingdingyang estuary of Zhujiang river mouth[J]. Journal of Tropical Oceanography, 2001, 20(3): 35-40.

    [13] 沈健, 沈煥庭, 潘定安, 等. 長(zhǎng)江河口最大渾濁帶水沙輸運(yùn)機(jī)制分析[J]. 地理學(xué)報(bào), 1995, 50(5): 411-420.

    Shen Jian, Shen Huanting, Pan Dingan, et al. Analysis of transport mechanism of water and suspended sediment in the turbidity maximum of the Changjiang estuary[J]. Acta Geographica Sinica, 1995, 50(5): 411-420.

    [14] 吳超羽, 任杰, 包蕓, 等. 珠江河口 “門” 的地貌動(dòng)力學(xué)初探[J]. 地理學(xué)報(bào), 2006, 61(5): 537-548.

    Wu Chaoyu, Ren Jie, Bao Yun, et al. A preliminary study on the morphodynamic evolution of the ‘Gate’ of the Pearl River Delta, China[J]. Acta Geographica Sinica, 2006, 61(5): 537-548.

    The 2D distribution characteristics and formation mechanism of suspended sediment concentration during the flood season in the Lingdingyang Estuary

    He Jiawei1,Jia Liangwen1,Wei Xiange2,Jia Yanhong1,Cheng Cong1

    (1.CoastalOceanResearchCenter,SchoolofMarineScience,SunYen-satUniversity,Guangzhou510275,China;2.SchoolofEnvironmentScienceandEngineering,SunYen-satUniversity,Guangzhou510275,China)

    Based on the observed data of current and suspended sediment concentration (SSC) during the flood season (13-14 and 16-17, August, 2007), by the means of measured data analysis and the TELEMAC-2D numerical model, this article reveals the spatial distribution characteristics of SSC in the Lingdingyang Esturay: in the longitudinal direction along west channel, SSC decrease first, then increase, and then decrease again from bay head to open sea; laterally, west channel, east channel and the north of the middle flat are higher SSC areas; the maximal center of SSC is located on the northwest of Neilingding Island. The model of longitudinal material transport reveals that the main reasons of net suspended sediment transport are advection, tidal asymmetry and vertical circulation. The effects of tidal current, runoff with sediment transport, solid boundary and topography condition are the main reasons for the 2D distribution characteristics of SSC in the Lingdingyang Estuary.

    Lingdingyang Estuary; suspended sediment concentration; 2D distribution characteristics; mechanism

    10.3969/j.issn.0253-4193.2017.09.003

    2016-11-07;

    2016-12-21。

    國(guó)家自然科學(xué)基金“復(fù)雜動(dòng)力條件下磨刀門河口長(zhǎng)周期地貌演變研究”(41376101);國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目“珠江河口與網(wǎng)河演變機(jī)制及治理研究”(2016YFC0402603);廣東省海洋漁業(yè)科技與產(chǎn)業(yè)發(fā)展專項(xiàng)“珠江口伶仃洋圍填海工程的累積環(huán)境影響及生態(tài)修復(fù)措施研究”(A201401D01);水利公益性行業(yè)科研專項(xiàng)經(jīng)費(fèi)項(xiàng)目“珠江河口咸情變化及抑咸對(duì)策研究”(201501010)。

    何嘉偉(1990—),男,廣東省深圳市人,研究助理,從事河口海岸動(dòng)力、泥沙研究。E-mail:jiawei_he@foxmail.com

    *通信作者:韋獻(xiàn)革(1966-),男,廣西省柳州市人,博士,從事環(huán)境科學(xué)研究。E-mail:eeswxg@mail.sysu.edu.cn

    TV142

    A

    0253-4193(2017)09-0026-14

    何嘉偉, 賈良文, 韋獻(xiàn)革, 等. 伶仃洋洪季懸沙平面分布特征及成因探討[J]. 海洋學(xué)報(bào), 2017, 39(9): 26-39,

    He Jiawei, Jia Liangwen, Wei Xiange, et al. The 2D distribution characteristics and formation mechanism of suspended sediment concentration during the flood season in the Lingdingyang Estuary[J]. Haiyang Xuebao, 2017, 39(9):26-39, doi:10.3969/j.issn.0253-4193.2017.09.003

    猜你喜歡
    懸沙伶仃洋沙場(chǎng)
    世界最高全飄浮伶仃洋大橋成功合龍
    軍事文摘(2023年16期)2023-02-19 04:13:14
    深中通道伶仃洋大橋貓道貫通
    近岸懸沙垂線分布多元線性回歸分析
    圍填海對(duì)伶仃洋水流動(dòng)力的短期影響模擬研究
    2017,朱日和沙場(chǎng)閱兵
    臺(tái)風(fēng)對(duì)長(zhǎng)江口表層懸沙濃度的影響
    港珠澳大橋:伶仃洋上架起的巨龍
    金橋(2018年12期)2019-01-29 02:47:28
    強(qiáng)軍夢(mèng)——朱日和沙場(chǎng)大閱兵
    沙場(chǎng)點(diǎn)兵
    沙場(chǎng)大點(diǎn)兵
    青年歌聲(2017年12期)2017-03-15 05:20:26
    国产欧美另类精品又又久久亚洲欧美| 在线精品无人区一区二区三| 永久免费av网站大全| 一级毛片 在线播放| 啦啦啦中文免费视频观看日本| 国产精品99久久99久久久不卡 | 日本欧美视频一区| 人体艺术视频欧美日本| 亚洲怡红院男人天堂| 久久综合国产亚洲精品| 色婷婷久久久亚洲欧美| 男女无遮挡免费网站观看| 夫妻性生交免费视频一级片| 日本vs欧美在线观看视频 | 日本-黄色视频高清免费观看| 久久久午夜欧美精品| 一级a做视频免费观看| 大香蕉97超碰在线| tube8黄色片| 黑丝袜美女国产一区| 看非洲黑人一级黄片| 校园人妻丝袜中文字幕| 欧美成人精品欧美一级黄| 高清视频免费观看一区二区| 美女主播在线视频| 99久久综合免费| 热re99久久精品国产66热6| 赤兔流量卡办理| 九九在线视频观看精品| 伊人亚洲综合成人网| 男女啪啪激烈高潮av片| 国产一区二区在线观看日韩| 国产亚洲91精品色在线| 色视频在线一区二区三区| 久久精品久久久久久噜噜老黄| 寂寞人妻少妇视频99o| 久久久午夜欧美精品| 精品久久久久久久久亚洲| 精品国产一区二区久久| 亚洲av电影在线观看一区二区三区| 午夜免费鲁丝| 国模一区二区三区四区视频| 韩国av在线不卡| 亚洲av在线观看美女高潮| 黄色欧美视频在线观看| 日韩欧美一区视频在线观看 | 国产精品成人在线| 又大又黄又爽视频免费| 免费不卡的大黄色大毛片视频在线观看| 免费久久久久久久精品成人欧美视频 | 美女xxoo啪啪120秒动态图| 一二三四中文在线观看免费高清| 大香蕉97超碰在线| 建设人人有责人人尽责人人享有的| 熟妇人妻不卡中文字幕| 国产乱人偷精品视频| 秋霞伦理黄片| 成人亚洲精品一区在线观看| 一本色道久久久久久精品综合| 国产在视频线精品| 日韩欧美一区视频在线观看 | 晚上一个人看的免费电影| 人妻制服诱惑在线中文字幕| 亚洲综合色惰| 简卡轻食公司| 亚洲无线观看免费| 少妇人妻 视频| 高清不卡的av网站| 黑人巨大精品欧美一区二区蜜桃 | 欧美激情国产日韩精品一区| 欧美日韩精品成人综合77777| 国产高清有码在线观看视频| 中文在线观看免费www的网站| 高清午夜精品一区二区三区| 国国产精品蜜臀av免费| 国产精品伦人一区二区| 黄色一级大片看看| 一本—道久久a久久精品蜜桃钙片| 亚洲婷婷狠狠爱综合网| 在线观看人妻少妇| 人妻夜夜爽99麻豆av| 精品久久久久久久久av| 亚洲精品中文字幕在线视频 | 我要看日韩黄色一级片| 女的被弄到高潮叫床怎么办| 九九在线视频观看精品| 丰满乱子伦码专区| 最近2019中文字幕mv第一页| 亚洲不卡免费看| 久久精品熟女亚洲av麻豆精品| 国产又色又爽无遮挡免| 日本黄大片高清| 日本vs欧美在线观看视频 | 日韩电影二区| 免费观看av网站的网址| 午夜日本视频在线| 一区在线观看完整版| 亚洲一区二区三区欧美精品| 蜜臀久久99精品久久宅男| av福利片在线| 免费久久久久久久精品成人欧美视频 | 国产成人aa在线观看| 久久久久久久久久久丰满| 色5月婷婷丁香| 亚洲精品乱码久久久久久按摩| 只有这里有精品99| 国产亚洲最大av| 日本黄色片子视频| 国产在线免费精品| 亚洲怡红院男人天堂| 高清欧美精品videossex| 久久婷婷青草| 3wmmmm亚洲av在线观看| 美女福利国产在线| 国产淫语在线视频| 69精品国产乱码久久久| 亚洲国产最新在线播放| 国产日韩欧美在线精品| 精品人妻偷拍中文字幕| 亚洲欧美日韩卡通动漫| 美女中出高潮动态图| 日韩av在线免费看完整版不卡| 国产精品一区二区在线不卡| 乱人伦中国视频| 中文字幕人妻丝袜制服| 欧美性感艳星| videos熟女内射| 永久免费av网站大全| 日韩强制内射视频| 高清在线视频一区二区三区| 久久久午夜欧美精品| 最近中文字幕高清免费大全6| 黑人巨大精品欧美一区二区蜜桃 | 久久久国产精品麻豆| 国产精品国产三级国产av玫瑰| 下体分泌物呈黄色| 亚洲性久久影院| 人人妻人人澡人人看| 精品久久久精品久久久| 中国美白少妇内射xxxbb| 久久国产乱子免费精品| 五月天丁香电影| 熟女av电影| 久久精品国产自在天天线| 国语对白做爰xxxⅹ性视频网站| 美女中出高潮动态图| 国产精品福利在线免费观看| 王馨瑶露胸无遮挡在线观看| 国产精品国产av在线观看| 精华霜和精华液先用哪个| 久久午夜福利片| 久久国产精品男人的天堂亚洲 | 五月开心婷婷网| 在线观看免费日韩欧美大片 | 一区二区三区乱码不卡18| 国产淫语在线视频| 草草在线视频免费看| 插阴视频在线观看视频| 精品少妇黑人巨大在线播放| 久久久久久久久大av| 人妻少妇偷人精品九色| 老司机影院毛片| 免费久久久久久久精品成人欧美视频 | 久久久久精品久久久久真实原创| 乱系列少妇在线播放| 国产一区二区三区综合在线观看 | 欧美bdsm另类| 国产成人aa在线观看| 久久久久国产网址| 国产亚洲午夜精品一区二区久久| 亚洲av综合色区一区| 99精国产麻豆久久婷婷| 中国国产av一级| 国产爽快片一区二区三区| 人妻夜夜爽99麻豆av| 少妇裸体淫交视频免费看高清| 乱人伦中国视频| a级毛片在线看网站| 国产高清有码在线观看视频| 赤兔流量卡办理| 国内揄拍国产精品人妻在线| 欧美日韩av久久| 成年人免费黄色播放视频 | 久久久久久久国产电影| 91午夜精品亚洲一区二区三区| 啦啦啦中文免费视频观看日本| 少妇精品久久久久久久| 蜜臀久久99精品久久宅男| 18禁在线无遮挡免费观看视频| 日韩电影二区| 久久久久久久久大av| 亚洲欧美日韩另类电影网站| 欧美区成人在线视频| 国国产精品蜜臀av免费| 97在线人人人人妻| 亚洲欧美精品专区久久| 亚洲四区av| 国产伦在线观看视频一区| 亚洲av.av天堂| 国产日韩欧美视频二区| 国产成人精品久久久久久| 国产精品女同一区二区软件| 美女中出高潮动态图| av福利片在线| 青春草国产在线视频| 色网站视频免费| 在线观看国产h片| av线在线观看网站| 哪个播放器可以免费观看大片| 亚洲欧美清纯卡通| 五月开心婷婷网| 欧美国产精品一级二级三级 | 久久久国产一区二区| 女的被弄到高潮叫床怎么办| 日产精品乱码卡一卡2卡三| 国产免费一区二区三区四区乱码| 日本黄色日本黄色录像| 国产精品三级大全| 在线观看人妻少妇| 亚洲欧美成人综合另类久久久| 国产精品久久久久成人av| 观看免费一级毛片| 黄色一级大片看看| 亚洲国产精品专区欧美| 欧美日韩国产mv在线观看视频| 久久久国产欧美日韩av| 国产精品久久久久久久电影| 亚洲美女视频黄频| 91精品国产国语对白视频| 亚洲av不卡在线观看| 亚洲性久久影院| 欧美亚洲 丝袜 人妻 在线| 伦理电影免费视频| 丝袜脚勾引网站| 蜜桃在线观看..| 久久99精品国语久久久| 99九九线精品视频在线观看视频| 一级av片app| 少妇高潮的动态图| 美女主播在线视频| 亚洲婷婷狠狠爱综合网| 国产免费又黄又爽又色| 我的老师免费观看完整版| 成人影院久久| a级毛色黄片| 一个人看视频在线观看www免费| 日韩av在线免费看完整版不卡| av一本久久久久| 秋霞在线观看毛片| 免费看不卡的av| 午夜福利网站1000一区二区三区| 成人国产av品久久久| av视频免费观看在线观看| 亚洲精品aⅴ在线观看| 国产精品三级大全| 涩涩av久久男人的天堂| 亚洲欧美成人综合另类久久久| 99热这里只有是精品在线观看| 永久网站在线| 午夜老司机福利剧场| 午夜福利,免费看| 一级片'在线观看视频| 国产爽快片一区二区三区| 日产精品乱码卡一卡2卡三| 日本欧美国产在线视频| 免费av不卡在线播放| 偷拍熟女少妇极品色| 国产色婷婷99| 色吧在线观看| 99久久人妻综合| 在线观看免费日韩欧美大片 | 99久久精品一区二区三区| 精品久久久久久电影网| 狠狠精品人妻久久久久久综合| 99精国产麻豆久久婷婷| 欧美另类一区| 另类精品久久| 国产黄片美女视频| 国产一区亚洲一区在线观看| 最近2019中文字幕mv第一页| kizo精华| 少妇人妻一区二区三区视频| 成人漫画全彩无遮挡| 曰老女人黄片| 欧美高清成人免费视频www| 国产精品成人在线| 高清不卡的av网站| av在线app专区| 永久免费av网站大全| 一级毛片电影观看| 亚洲精品中文字幕在线视频 | 中国三级夫妇交换| 熟女人妻精品中文字幕| .国产精品久久| 嫩草影院新地址| 亚洲欧美成人精品一区二区| 国产亚洲最大av| 中文字幕久久专区| 丰满人妻一区二区三区视频av| 熟妇人妻不卡中文字幕| 少妇丰满av| 亚洲av综合色区一区| 男女边摸边吃奶| 2021少妇久久久久久久久久久| 性高湖久久久久久久久免费观看| 91精品国产国语对白视频| 久久精品国产亚洲av涩爱| 夫妻午夜视频| a级一级毛片免费在线观看| 大码成人一级视频| 成人美女网站在线观看视频| 色婷婷av一区二区三区视频| 久久久久久久国产电影| 日本午夜av视频| 久久精品国产亚洲网站| 美女大奶头黄色视频| 日韩亚洲欧美综合| 寂寞人妻少妇视频99o| 一本久久精品| 欧美三级亚洲精品| 国产精品熟女久久久久浪| 大陆偷拍与自拍| 国精品久久久久久国模美| 麻豆成人av视频| 日韩强制内射视频| 91精品一卡2卡3卡4卡| 黑人猛操日本美女一级片| 美女中出高潮动态图| 各种免费的搞黄视频| 亚洲国产精品国产精品| 桃花免费在线播放| 91成人精品电影| 色婷婷av一区二区三区视频| 乱码一卡2卡4卡精品| 精品久久久久久电影网| 涩涩av久久男人的天堂| 一本色道久久久久久精品综合| 十八禁高潮呻吟视频 | av专区在线播放| 啦啦啦在线观看免费高清www| 新久久久久国产一级毛片| 亚洲美女视频黄频| 99久久人妻综合| 成人国产麻豆网| 在线亚洲精品国产二区图片欧美 | 一二三四中文在线观看免费高清| 国产精品.久久久| 啦啦啦中文免费视频观看日本| 亚洲久久久国产精品| 美女xxoo啪啪120秒动态图| kizo精华| 丰满迷人的少妇在线观看| 80岁老熟妇乱子伦牲交| √禁漫天堂资源中文www| 男女无遮挡免费网站观看| 亚洲av福利一区| 夜夜爽夜夜爽视频| 久久久久视频综合| 欧美成人精品欧美一级黄| a级毛片在线看网站| 91在线精品国自产拍蜜月| 99热这里只有是精品50| 91aial.com中文字幕在线观看| 丰满少妇做爰视频| 少妇的逼水好多| 国内揄拍国产精品人妻在线| 亚洲精品日本国产第一区| 欧美人与善性xxx| 国产精品国产三级国产av玫瑰| 免费在线观看成人毛片| 美女脱内裤让男人舔精品视频| 国国产精品蜜臀av免费| 成年av动漫网址| 日韩一本色道免费dvd| 嫩草影院新地址| 一个人免费看片子| 午夜av观看不卡| 男人和女人高潮做爰伦理| 亚洲久久久国产精品| 有码 亚洲区| 午夜精品国产一区二区电影| 天天操日日干夜夜撸| 天堂8中文在线网| 亚洲精品一区蜜桃| freevideosex欧美| 最后的刺客免费高清国语| 久久av网站| 国产国拍精品亚洲av在线观看| 特大巨黑吊av在线直播| 成人毛片60女人毛片免费| 日日爽夜夜爽网站| 男人和女人高潮做爰伦理| 国产成人freesex在线| 亚洲av福利一区| 久久精品久久久久久久性| 亚洲第一区二区三区不卡| 亚洲国产欧美在线一区| 热re99久久精品国产66热6| 大话2 男鬼变身卡| 啦啦啦视频在线资源免费观看| 国产精品人妻久久久久久| 九色成人免费人妻av| h视频一区二区三区| 黑人猛操日本美女一级片| 亚洲精品一二三| 天堂8中文在线网| 免费大片黄手机在线观看| 久久韩国三级中文字幕| 成人亚洲欧美一区二区av| 最近手机中文字幕大全| 国产熟女欧美一区二区| 日韩成人av中文字幕在线观看| 日韩一区二区视频免费看| 久久久久久伊人网av| 日日啪夜夜爽| 一边亲一边摸免费视频| 日日摸夜夜添夜夜添av毛片| 国产成人免费无遮挡视频| 国产高清不卡午夜福利| 老司机亚洲免费影院| 久久久亚洲精品成人影院| 日韩一区二区三区影片| 日本午夜av视频| av线在线观看网站| 又粗又硬又长又爽又黄的视频| 一级毛片aaaaaa免费看小| 国产亚洲最大av| 在线精品无人区一区二区三| 大片免费播放器 马上看| 亚洲电影在线观看av| 亚洲精品日本国产第一区| 亚洲精品自拍成人| 女的被弄到高潮叫床怎么办| av免费在线看不卡| 两个人免费观看高清视频 | 国产在线免费精品| 欧美 亚洲 国产 日韩一| 天天躁夜夜躁狠狠久久av| 人妻系列 视频| 日本午夜av视频| 伊人久久国产一区二区| 亚洲内射少妇av| 97超碰精品成人国产| 一区二区三区四区激情视频| 美女脱内裤让男人舔精品视频| 久久久久网色| 伊人亚洲综合成人网| 国产片特级美女逼逼视频| 免费观看在线日韩| 日韩欧美精品免费久久| 人妻一区二区av| 高清av免费在线| 国产一区有黄有色的免费视频| 日本vs欧美在线观看视频 | 观看av在线不卡| 欧美3d第一页| 蜜臀久久99精品久久宅男| 乱系列少妇在线播放| 在线看a的网站| 夜夜骑夜夜射夜夜干| av福利片在线观看| 亚洲欧美日韩卡通动漫| 伦精品一区二区三区| 亚洲国产精品一区三区| 在线观看免费视频网站a站| 国产男人的电影天堂91| 久久久a久久爽久久v久久| 亚洲伊人久久精品综合| h视频一区二区三区| 国产综合精华液| 黄色日韩在线| 水蜜桃什么品种好| 欧美bdsm另类| 国产色婷婷99| 亚洲国产色片| 亚洲真实伦在线观看| 蜜桃久久精品国产亚洲av| 中国三级夫妇交换| a级片在线免费高清观看视频| 国产av国产精品国产| 免费观看av网站的网址| 麻豆乱淫一区二区| 国产精品秋霞免费鲁丝片| 男人添女人高潮全过程视频| 在线观看美女被高潮喷水网站| 国产白丝娇喘喷水9色精品| 日本欧美国产在线视频| 欧美精品一区二区免费开放| 精品亚洲成国产av| 多毛熟女@视频| 黄色配什么色好看| 成年人午夜在线观看视频| 亚洲精品,欧美精品| 丝瓜视频免费看黄片| 少妇的逼好多水| 十八禁网站网址无遮挡 | 国产精品国产三级专区第一集| 国产爽快片一区二区三区| 国产亚洲91精品色在线| 免费看光身美女| 国语对白做爰xxxⅹ性视频网站| 日韩av在线免费看完整版不卡| 啦啦啦啦在线视频资源| 日本av免费视频播放| 亚洲精品中文字幕在线视频 | 中文字幕av电影在线播放| 日本黄色日本黄色录像| 一二三四中文在线观看免费高清| 黄色日韩在线| 最后的刺客免费高清国语| 国产av码专区亚洲av| 偷拍熟女少妇极品色| 水蜜桃什么品种好| 国内少妇人妻偷人精品xxx网站| 中文乱码字字幕精品一区二区三区| 99视频精品全部免费 在线| 精品卡一卡二卡四卡免费| av福利片在线观看| 高清欧美精品videossex| 两个人的视频大全免费| 国产美女午夜福利| 亚洲av中文av极速乱| 亚洲欧美日韩另类电影网站| 亚洲欧美精品自产自拍| 国产男人的电影天堂91| 午夜视频国产福利| 久久青草综合色| 夫妻性生交免费视频一级片| 日韩欧美精品免费久久| 国产男女超爽视频在线观看| 婷婷色综合www| 天天躁夜夜躁狠狠久久av| 国产日韩欧美视频二区| av福利片在线| 精品亚洲乱码少妇综合久久| av网站免费在线观看视频| 简卡轻食公司| 国产精品偷伦视频观看了| 国产精品一区www在线观看| 国产在线一区二区三区精| 中国三级夫妇交换| 在线观看一区二区三区激情| 亚洲国产欧美在线一区| 欧美变态另类bdsm刘玥| 精品一区在线观看国产| 国产成人午夜福利电影在线观看| 男女边吃奶边做爰视频| 人妻 亚洲 视频| 黑人猛操日本美女一级片| 免费看不卡的av| 欧美性感艳星| 中国美白少妇内射xxxbb| 成人影院久久| 高清视频免费观看一区二区| 青青草视频在线视频观看| 国产精品人妻久久久久久| 激情五月婷婷亚洲| 人妻 亚洲 视频| 亚洲精品国产色婷婷电影| 中文资源天堂在线| 国产精品一区二区在线不卡| 伦精品一区二区三区| 日日摸夜夜添夜夜爱| 人体艺术视频欧美日本| 男男h啪啪无遮挡| 亚洲成色77777| 王馨瑶露胸无遮挡在线观看| 最近中文字幕2019免费版| 我要看黄色一级片免费的| 成人毛片60女人毛片免费| 一级黄片播放器| 如日韩欧美国产精品一区二区三区 | 日韩欧美 国产精品| 最近的中文字幕免费完整| 自拍偷自拍亚洲精品老妇| 女的被弄到高潮叫床怎么办| 99热网站在线观看| 日韩精品免费视频一区二区三区 | 精品人妻熟女毛片av久久网站| 久久久久人妻精品一区果冻| 欧美激情极品国产一区二区三区 | 亚洲激情五月婷婷啪啪| 色视频www国产| 又大又黄又爽视频免费| 91成人精品电影| 精品久久国产蜜桃| 国产精品人妻久久久久久| 校园人妻丝袜中文字幕| 99精国产麻豆久久婷婷| 久久久久精品性色| 亚洲精品自拍成人| 美女视频免费永久观看网站| 极品人妻少妇av视频| 亚洲三级黄色毛片| 一级二级三级毛片免费看| av福利片在线观看| 亚洲国产精品成人久久小说| 汤姆久久久久久久影院中文字幕| 亚洲国产欧美日韩在线播放 | 久久精品国产a三级三级三级| 爱豆传媒免费全集在线观看| 大又大粗又爽又黄少妇毛片口| 精品久久国产蜜桃| 亚洲久久久国产精品| 嫩草影院入口| 制服丝袜香蕉在线| 狂野欧美白嫩少妇大欣赏| 欧美精品一区二区免费开放| 久久久久久久国产电影| 丰满人妻一区二区三区视频av| 国产欧美日韩精品一区二区| 搡老乐熟女国产| 国产精品一区二区三区四区免费观看| 亚洲高清免费不卡视频| 亚洲电影在线观看av|