• <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不卡免费在线播放| 欧美国产精品va在线观看不卡| 欧美日韩黄片免| 青春草视频在线免费观看| 国产高清不卡午夜福利| 国产欧美日韩综合在线一区二区| 9色porny在线观看| 亚洲欧美激情在线| 国产精品国产av在线观看| 亚洲国产精品一区二区三区在线| 一本色道久久久久久精品综合| 免费一级毛片在线播放高清视频 | 一本—道久久a久久精品蜜桃钙片| 久久精品久久精品一区二区三区| 另类精品久久| av欧美777| 国产成人精品无人区| 免费高清在线观看视频在线观看| 搡老乐熟女国产| 国产精品一区二区精品视频观看| 日韩欧美一区视频在线观看| 高清黄色对白视频在线免费看| 成人国语在线视频| 午夜福利影视在线免费观看| 女人高潮潮喷娇喘18禁视频| 最近中文字幕2019免费版| 亚洲天堂av无毛| 精品第一国产精品| 国产亚洲一区二区精品| 亚洲第一青青草原| 国产成人影院久久av| 自拍欧美九色日韩亚洲蝌蚪91| 午夜91福利影院| 国产高清视频在线播放一区 | 性色av乱码一区二区三区2| 高清黄色对白视频在线免费看| 欧美黄色片欧美黄色片| 国产精品国产av在线观看| 国产不卡av网站在线观看| 亚洲视频免费观看视频| 免费在线观看完整版高清| 国产日韩欧美在线精品| 精品国产一区二区三区四区第35| 免费在线观看完整版高清| 丝袜美腿诱惑在线| 亚洲精品自拍成人| 婷婷色综合大香蕉| 一本—道久久a久久精品蜜桃钙片| 777米奇影视久久| av又黄又爽大尺度在线免费看| 首页视频小说图片口味搜索 | 国产99久久九九免费精品| 日日夜夜操网爽| 久久精品亚洲熟妇少妇任你| 97人妻天天添夜夜摸| 亚洲一区中文字幕在线| 2021少妇久久久久久久久久久| 2021少妇久久久久久久久久久| 在现免费观看毛片| 亚洲国产精品成人久久小说| 久久国产精品影院| 丝袜在线中文字幕| 日韩中文字幕欧美一区二区 | 国产极品粉嫩免费观看在线| 女人精品久久久久毛片| 精品一区在线观看国产| netflix在线观看网站| 精品少妇一区二区三区视频日本电影| 日韩制服丝袜自拍偷拍| 9热在线视频观看99| 9191精品国产免费久久| 亚洲精品久久午夜乱码| 国精品久久久久久国模美| 欧美精品亚洲一区二区| 亚洲精品一区蜜桃| 少妇被粗大的猛进出69影院| 欧美日韩综合久久久久久| 超色免费av| 午夜激情久久久久久久| 久久热在线av| 高清不卡的av网站| 免费在线观看视频国产中文字幕亚洲 | 人人妻,人人澡人人爽秒播 | av网站免费在线观看视频| 亚洲av成人精品一二三区| 777米奇影视久久| 高潮久久久久久久久久久不卡| 婷婷色综合www| 国产亚洲一区二区精品| 一区二区日韩欧美中文字幕| 午夜福利免费观看在线| 亚洲av国产av综合av卡| 免费人妻精品一区二区三区视频| av又黄又爽大尺度在线免费看| av天堂久久9| 国产精品 欧美亚洲| 国产在视频线精品| 欧美日韩黄片免| 国产深夜福利视频在线观看| 老汉色av国产亚洲站长工具| 国产成人精品在线电影| 精品亚洲成国产av| 精品少妇黑人巨大在线播放| 男女免费视频国产| 国产熟女欧美一区二区| 亚洲av欧美aⅴ国产| 国产成人免费观看mmmm| 啦啦啦在线免费观看视频4| 五月天丁香电影| 久久久久国产精品人妻一区二区| 国产精品一区二区精品视频观看| www.自偷自拍.com| 欧美黑人精品巨大| 男女床上黄色一级片免费看| 国产欧美日韩一区二区三 | 亚洲伊人色综图| 精品国产一区二区三区四区第35| 国产三级黄色录像| 日韩av在线免费看完整版不卡| 成人国语在线视频| 欧美精品人与动牲交sv欧美| 国产国语露脸激情在线看| 一本色道久久久久久精品综合| 丰满迷人的少妇在线观看| 亚洲欧美精品综合一区二区三区| 国产在线一区二区三区精| 亚洲伊人色综图| 成年动漫av网址| 国产午夜精品一二区理论片| av有码第一页| 免费观看a级毛片全部| 亚洲精品国产色婷婷电影| 国语对白做爰xxxⅹ性视频网站| 99久久人妻综合| av有码第一页| 亚洲欧美成人综合另类久久久| 人人妻人人添人人爽欧美一区卜| 久久人人爽人人片av| 国产精品香港三级国产av潘金莲 | 色播在线永久视频| 一边摸一边做爽爽视频免费| 考比视频在线观看| 免费看十八禁软件| 爱豆传媒免费全集在线观看| 国产成人免费无遮挡视频| 麻豆乱淫一区二区| 久久久久久人人人人人| 国产精品香港三级国产av潘金莲 | 亚洲 欧美一区二区三区| 亚洲欧美清纯卡通| 自线自在国产av| 19禁男女啪啪无遮挡网站| 亚洲黑人精品在线| 纵有疾风起免费观看全集完整版| 精品久久久久久电影网| 黄片播放在线免费| 免费日韩欧美在线观看| 久久国产精品男人的天堂亚洲| 中文字幕另类日韩欧美亚洲嫩草| 青草久久国产| 亚洲国产欧美在线一区| 国产在视频线精品| 老司机靠b影院| 国产xxxxx性猛交| 老汉色av国产亚洲站长工具| 狠狠婷婷综合久久久久久88av| 色精品久久人妻99蜜桃| 亚洲人成电影观看| 国产成人av教育| e午夜精品久久久久久久| 男女国产视频网站| 热99国产精品久久久久久7| 日韩制服丝袜自拍偷拍| 少妇被粗大的猛进出69影院| 亚洲欧美一区二区三区黑人| 后天国语完整版免费观看| 日本猛色少妇xxxxx猛交久久| 成年人午夜在线观看视频| 少妇 在线观看| 十八禁高潮呻吟视频| 日本a在线网址| 亚洲国产欧美在线一区| 欧美国产精品va在线观看不卡| h视频一区二区三区| 亚洲激情五月婷婷啪啪| 色综合欧美亚洲国产小说| 看免费av毛片| 侵犯人妻中文字幕一二三四区| 美女午夜性视频免费| 亚洲精品美女久久久久99蜜臀 | 亚洲欧美精品综合一区二区三区| 亚洲国产精品999| 视频区图区小说| 热99国产精品久久久久久7| 日韩欧美一区视频在线观看| 久久精品亚洲av国产电影网| 男男h啪啪无遮挡| 少妇精品久久久久久久| 国产黄色免费在线视频| 人人妻人人添人人爽欧美一区卜| 亚洲av成人不卡在线观看播放网 | 人人妻人人澡人人看| 老司机深夜福利视频在线观看 | 久久鲁丝午夜福利片| 麻豆国产av国片精品| 多毛熟女@视频| 一级毛片 在线播放| 午夜久久久在线观看| 一本综合久久免费| 黑人巨大精品欧美一区二区蜜桃| 人人妻人人爽人人添夜夜欢视频| 亚洲一码二码三码区别大吗| 两人在一起打扑克的视频| 丁香六月欧美| 精品国产超薄肉色丝袜足j| 赤兔流量卡办理| 你懂的网址亚洲精品在线观看| 女警被强在线播放| 黄色毛片三级朝国网站| 999久久久国产精品视频| 久久精品aⅴ一区二区三区四区| 久久精品国产亚洲av高清一级| 桃花免费在线播放| 欧美精品啪啪一区二区三区 | 最近手机中文字幕大全| 另类精品久久| www.av在线官网国产| 午夜av观看不卡| 搡老乐熟女国产| 每晚都被弄得嗷嗷叫到高潮| 叶爱在线成人免费视频播放| 操美女的视频在线观看| 一级毛片我不卡| 亚洲精品中文字幕在线视频| 91精品三级在线观看| 日韩av在线免费看完整版不卡| 欧美另类一区| 国产精品免费视频内射| 新久久久久国产一级毛片| 精品久久久久久久毛片微露脸 | 美女脱内裤让男人舔精品视频| 老司机午夜十八禁免费视频| 啦啦啦在线免费观看视频4| 悠悠久久av| 19禁男女啪啪无遮挡网站| 99热网站在线观看| 美国免费a级毛片| 美女午夜性视频免费| 亚洲激情五月婷婷啪啪| 一级毛片电影观看| 欧美成人精品欧美一级黄| 国产高清不卡午夜福利| 如日韩欧美国产精品一区二区三区| 国产人伦9x9x在线观看| 少妇 在线观看| 少妇被粗大的猛进出69影院| 日本午夜av视频| 黄色毛片三级朝国网站| 韩国精品一区二区三区| 成人亚洲欧美一区二区av| 亚洲欧洲国产日韩| 欧美国产精品va在线观看不卡| 国产一区有黄有色的免费视频| 1024视频免费在线观看| 色视频在线一区二区三区| 亚洲精品国产区一区二| 女性生殖器流出的白浆| 久久久久久亚洲精品国产蜜桃av| 不卡av一区二区三区| 亚洲少妇的诱惑av| 黑丝袜美女国产一区| av又黄又爽大尺度在线免费看| 人妻人人澡人人爽人人| 国产日韩欧美在线精品| 老司机影院毛片| 午夜av观看不卡| 久久99热这里只频精品6学生| 国产亚洲av片在线观看秒播厂| 黄频高清免费视频| 99久久精品国产亚洲精品| 午夜福利影视在线免费观看| 嫩草影视91久久| 黄频高清免费视频| 国产精品三级大全| 亚洲欧美精品综合一区二区三区| 午夜福利一区二区在线看| 男女边摸边吃奶| 亚洲精品久久午夜乱码| 日韩制服丝袜自拍偷拍| 色94色欧美一区二区| 午夜老司机福利片| 午夜91福利影院| 婷婷成人精品国产| 久久天堂一区二区三区四区| 久久免费观看电影| 亚洲综合色网址| 一区二区三区乱码不卡18| 国产欧美日韩一区二区三区在线| 欧美人与善性xxx| 制服诱惑二区| 久久青草综合色| 亚洲国产欧美日韩在线播放| 黄色一级大片看看| 亚洲欧洲国产日韩| 亚洲精品国产av成人精品| 精品一品国产午夜福利视频| 婷婷色麻豆天堂久久| 9色porny在线观看| 国产主播在线观看一区二区 | 国产主播在线观看一区二区 | 国产成人一区二区在线| 久久影院123| 日韩精品免费视频一区二区三区| 亚洲综合色网址| 精品亚洲成国产av| 亚洲精品乱久久久久久| 在线观看人妻少妇| 一本—道久久a久久精品蜜桃钙片| 自线自在国产av| 69精品国产乱码久久久| 国产成人欧美| 久久精品亚洲av国产电影网| 久久久久久免费高清国产稀缺| 一级,二级,三级黄色视频| 美女国产高潮福利片在线看| cao死你这个sao货| 手机成人av网站| 亚洲av日韩在线播放| 一本大道久久a久久精品| 老司机影院毛片| 丁香六月欧美| 日韩,欧美,国产一区二区三区| 久久久久久久久免费视频了| 十八禁网站网址无遮挡| 日日夜夜操网爽| 精品少妇黑人巨大在线播放| 亚洲色图综合在线观看| 国产欧美亚洲国产| 少妇猛男粗大的猛烈进出视频| 国产精品二区激情视频| 日本av手机在线免费观看| 国产欧美日韩一区二区三 | 亚洲精品av麻豆狂野| 欧美黄色片欧美黄色片| 天堂8中文在线网| 欧美变态另类bdsm刘玥| 欧美国产精品一级二级三级| 80岁老熟妇乱子伦牲交| 国产精品一区二区在线不卡| 国产伦理片在线播放av一区| svipshipincom国产片| 精品一区二区三卡| 2018国产大陆天天弄谢| 一区二区日韩欧美中文字幕| 夫妻午夜视频| av欧美777| 免费久久久久久久精品成人欧美视频| av视频免费观看在线观看| av网站在线播放免费| 日韩视频在线欧美| 欧美日韩亚洲国产一区二区在线观看 | 19禁男女啪啪无遮挡网站| 看免费av毛片| 日韩,欧美,国产一区二区三区| 日本av手机在线免费观看| 久久精品国产亚洲av涩爱| 色视频在线一区二区三区| 99国产精品免费福利视频| 午夜免费鲁丝| 一二三四社区在线视频社区8| 我的亚洲天堂| 国产在线视频一区二区| 午夜福利视频精品| 最近中文字幕2019免费版| 少妇猛男粗大的猛烈进出视频| 女性生殖器流出的白浆| 一区二区三区乱码不卡18| 欧美日韩亚洲综合一区二区三区_| 国产免费一区二区三区四区乱码| 免费女性裸体啪啪无遮挡网站| 欧美 亚洲 国产 日韩一| 亚洲欧美精品自产自拍| 亚洲精品日本国产第一区| 交换朋友夫妻互换小说| 九色亚洲精品在线播放| 免费高清在线观看视频在线观看| 黄色 视频免费看| 免费观看a级毛片全部| 亚洲国产成人一精品久久久| 国产成人精品久久久久久| 国产成人精品久久二区二区免费| 老司机影院毛片| 一本久久精品| 性色av乱码一区二区三区2| www日本在线高清视频| 51午夜福利影视在线观看| 欧美人与性动交α欧美软件| 欧美精品人与动牲交sv欧美| 亚洲熟女毛片儿| 在线观看免费日韩欧美大片| 悠悠久久av| 18禁黄网站禁片午夜丰满| 欧美+亚洲+日韩+国产| 日韩制服骚丝袜av| 日本欧美视频一区| 伦理电影免费视频| 日本黄色日本黄色录像| 亚洲欧美一区二区三区久久| 51午夜福利影视在线观看| 亚洲欧美中文字幕日韩二区| 一区二区三区精品91| 美女视频免费永久观看网站| 亚洲伊人久久精品综合| 国产日韩欧美亚洲二区| 成人三级做爰电影| 亚洲成人国产一区在线观看 | 久久人妻福利社区极品人妻图片 | 制服诱惑二区| 男人舔女人的私密视频| 可以免费在线观看a视频的电影网站| 精品国产一区二区三区四区第35| 搡老乐熟女国产| 久久精品国产亚洲av涩爱| 国产成人免费无遮挡视频| 欧美日韩综合久久久久久| 婷婷色综合www| cao死你这个sao货| 丰满饥渴人妻一区二区三| 成人黄色视频免费在线看| 在线亚洲精品国产二区图片欧美| 最新在线观看一区二区三区 | 久久这里只有精品19| 国产在视频线精品| 国产欧美亚洲国产| 欧美乱码精品一区二区三区| 欧美精品一区二区大全| 秋霞在线观看毛片| 69精品国产乱码久久久| 久久久国产欧美日韩av| 日本猛色少妇xxxxx猛交久久| 热99久久久久精品小说推荐| 国产精品欧美亚洲77777| 久久毛片免费看一区二区三区| 欧美大码av| 成年人免费黄色播放视频| 中文字幕高清在线视频| 成年美女黄网站色视频大全免费| 老司机在亚洲福利影院| 久久九九热精品免费| 久久久久久人人人人人| 亚洲成国产人片在线观看| 精品视频人人做人人爽| 我的亚洲天堂| 亚洲九九香蕉| 又紧又爽又黄一区二区| 国产极品粉嫩免费观看在线| 久久国产精品男人的天堂亚洲| 伦理电影免费视频| 美女高潮到喷水免费观看| 欧美乱码精品一区二区三区| 日韩人妻精品一区2区三区| 国产精品国产三级国产专区5o| 日本色播在线视频| 亚洲精品国产av蜜桃| 亚洲av男天堂| 国产精品久久久人人做人人爽| 久久综合国产亚洲精品| 99精国产麻豆久久婷婷| 精品福利观看| 日本av免费视频播放| 亚洲国产精品国产精品| svipshipincom国产片| 国产欧美日韩精品亚洲av| 午夜日韩欧美国产| 亚洲一码二码三码区别大吗| 婷婷色av中文字幕| 极品人妻少妇av视频| 成年动漫av网址| 黄色a级毛片大全视频| 人人妻人人澡人人爽人人夜夜| 老司机影院成人| 这个男人来自地球电影免费观看| 在线亚洲精品国产二区图片欧美| 777米奇影视久久| 免费在线观看完整版高清| 丝袜美足系列| 久久免费观看电影| 亚洲专区中文字幕在线| 91麻豆av在线| 精品少妇黑人巨大在线播放| 高清欧美精品videossex| 又大又爽又粗| 亚洲国产精品国产精品| 激情五月婷婷亚洲| 看十八女毛片水多多多| 欧美日韩精品网址| 一区在线观看完整版| tube8黄色片| 中国国产av一级| 亚洲av电影在线进入| 老司机在亚洲福利影院| 99热国产这里只有精品6| 欧美在线黄色| 国产99久久九九免费精品| 国产熟女欧美一区二区| 黄网站色视频无遮挡免费观看| 精品一区二区三区四区五区乱码 | 日韩大片免费观看网站| 亚洲人成网站在线观看播放| 99久久综合免费| cao死你这个sao货| 纯流量卡能插随身wifi吗| 少妇的丰满在线观看| 亚洲国产精品999| 91九色精品人成在线观看| 亚洲av国产av综合av卡| 午夜免费成人在线视频| 中文字幕最新亚洲高清| 国产av精品麻豆| 99热全是精品| 七月丁香在线播放| 尾随美女入室| 国产精品久久久久成人av| 久久精品成人免费网站| 精品亚洲乱码少妇综合久久| 久久ye,这里只有精品| 免费高清在线观看日韩| 99热全是精品| 亚洲免费av在线视频| 久久天堂一区二区三区四区| 9191精品国产免费久久| 亚洲精品日韩在线中文字幕| 黑人巨大精品欧美一区二区蜜桃| 国产精品 欧美亚洲| videosex国产| 日韩 亚洲 欧美在线| 国产成人一区二区三区免费视频网站 | 欧美在线一区亚洲| 久久久久久久国产电影| 多毛熟女@视频| 成人亚洲精品一区在线观看| 亚洲欧美日韩另类电影网站| 亚洲av成人精品一二三区| 国产精品久久久人人做人人爽| 丝袜美足系列| 十八禁人妻一区二区| 在线观看人妻少妇| 黄色a级毛片大全视频| 久久精品国产亚洲av高清一级| 每晚都被弄得嗷嗷叫到高潮| 国产激情久久老熟女| 视频区欧美日本亚洲| 久久国产精品影院| 国产成人免费观看mmmm| 男女床上黄色一级片免费看| 亚洲精品日本国产第一区| 中文字幕人妻熟女乱码| 久久狼人影院| 日日爽夜夜爽网站| 日韩一本色道免费dvd| 亚洲av在线观看美女高潮| 国产精品一区二区在线不卡| 免费观看av网站的网址| 晚上一个人看的免费电影| 亚洲国产中文字幕在线视频| 久久女婷五月综合色啪小说| 秋霞在线观看毛片| 一本色道久久久久久精品综合| 亚洲精品美女久久久久99蜜臀 | 国产熟女午夜一区二区三区| 久久狼人影院| 好男人视频免费观看在线| 超色免费av| 国产一区有黄有色的免费视频| 久久99热这里只频精品6学生| 欧美精品一区二区大全| 亚洲av在线观看美女高潮| 欧美久久黑人一区二区| 久久99一区二区三区| 久久久久久久国产电影| 成人影院久久| 亚洲欧美一区二区三区黑人| 男人舔女人的私密视频| 在线观看国产h片| 欧美xxⅹ黑人| 91字幕亚洲| 青春草亚洲视频在线观看| 欧美xxⅹ黑人| 久久人人97超碰香蕉20202| 欧美精品av麻豆av| 少妇裸体淫交视频免费看高清 | 精品国产一区二区三区四区第35| 真人做人爱边吃奶动态| 国产精品国产三级专区第一集| 在线观看一区二区三区激情| 国产免费视频播放在线视频| 亚洲人成电影观看| 国产又爽黄色视频| 精品少妇内射三级| 免费高清在线观看视频在线观看| 国产又爽黄色视频| 精品少妇内射三级| 国产一区二区三区av在线| 欧美在线黄色| 亚洲国产欧美网| 午夜av观看不卡| av视频免费观看在线观看| 欧美人与性动交α欧美软件| 免费在线观看黄色视频的| 另类亚洲欧美激情| 亚洲精品国产av成人精品| 操出白浆在线播放| 黄片小视频在线播放| 亚洲午夜精品一区,二区,三区| 中国美女看黄片|