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

    珠江口磨刀門泥沙絮凝特征

    2015-01-05 02:33:06鄧智瑞何青楊清書林建良
    海洋學(xué)報(bào) 2015年9期
    關(guān)鍵詞:小潮絮團(tuán)大潮

    鄧智瑞,何青*,楊清書,林建良

    (1.華東師范大學(xué) 河口海岸學(xué)國家重點(diǎn)實(shí)驗(yàn)室,上海 200062;2.中山大學(xué) 河口海岸研究所,廣東 廣州 510006)

    珠江口磨刀門泥沙絮凝特征

    鄧智瑞1,何青1*,楊清書2,林建良1

    (1.華東師范大學(xué) 河口海岸學(xué)國家重點(diǎn)實(shí)驗(yàn)室,上海 200062;2.中山大學(xué) 河口海岸研究所,廣東 廣州 510006)

    本文利用激光粒度儀實(shí)測得到珠江口磨刀門河口2013年夏季懸浮泥沙現(xiàn)場絮凝及絮凝體特征,同時(shí)對比懸沙分散粒徑和含沙量,研究表明:懸沙分散粒徑平均值為27.9 μm,現(xiàn)場實(shí)測絮團(tuán)粒徑平均值為91.6 μm,表明磨刀門口外的懸浮泥沙絮凝現(xiàn)象顯著;實(shí)測絮團(tuán)平均粒徑變化范圍為13.0~273.8 μm,小潮期間絮團(tuán)粒徑平均值為131.5 μm,大于大潮平均值76.9 μm;絮凝體粒徑在垂向上的變化表現(xiàn)為由表及底先變大再變小。絮團(tuán)體積濃度、沉速與粒徑的關(guān)系在不同情況下有差異,體積濃度和絮團(tuán)粒徑在表層和中層有明顯正相關(guān)關(guān)系,絮團(tuán)沉速在大潮時(shí)刻隨著粒徑的增大而增大。綜合分析影響絮凝的因素,得知在珠江口鹽度對于絮團(tuán)大小影響不明顯;而流速大小的差異是影響大小潮之間絮團(tuán)大小不同的主要因素。研究結(jié)果有助于了解珠江口細(xì)顆粒泥沙輸移特性和相關(guān)生物化學(xué)過程。

    珠江口;磨刀門;懸浮泥沙;絮凝

    1 引言

    黏性細(xì)顆粒泥沙一般指小于62 μm的、具有顆粒間黏性的泥沙,其絮凝現(xiàn)象是河流和河口懸沙起動(dòng)、運(yùn)移和沉積的一個(gè)重要過程。懸浮泥沙顆粒在顆粒間電荷引力、水體膠質(zhì)等的作用下,多顆粒結(jié)合在一起,形成更大顆粒、低密度(相對于泥沙顆粒)的絮凝體,即絮團(tuán)。絮團(tuán)的沉降速度要遠(yuǎn)遠(yuǎn)大于分散的泥沙顆粒,因此研究懸浮泥沙絮凝一直是細(xì)顆粒泥沙研究中的一個(gè)重要內(nèi)容。同時(shí)絮凝體也促成物質(zhì)吸附作用,因此作為水體污染物和營養(yǎng)鹽等物質(zhì)輸送的載體,對河口水質(zhì)環(huán)境和生物化學(xué)物質(zhì)輸運(yùn)等具有重要作用。

    絮凝體的組成、動(dòng)力過程和形成機(jī)制在室內(nèi)實(shí)驗(yàn)、現(xiàn)場觀測和數(shù)值模擬等方面已經(jīng)有大量研究[1—3]。傳統(tǒng)室內(nèi)實(shí)驗(yàn)方法可獲取懸沙的級配和沉速數(shù)據(jù),但其工序相當(dāng)復(fù)雜,從水樣采集、保存到分析,脆弱的絮凝體均可能被打碎,因而分析得到的數(shù)據(jù)往往不能真實(shí)反映絮凝體的特性[4]。現(xiàn)場激光粒度儀LISST-100則可以現(xiàn)場實(shí)時(shí)獲取水中懸浮顆粒物現(xiàn)場絮團(tuán)粒徑譜和沉速譜[5],功能較全面,測驗(yàn)效率高,精度好。程江[6]使用LISST-100在長江口現(xiàn)場觀測到絮團(tuán),并使用譜分析方法研究了絮團(tuán)隨潮流的周期變化。唐建華[7]分析了水動(dòng)力、懸浮物濃度、鹽度等因素對絮團(tuán)的影響,并且提出了基于現(xiàn)場資料的絮團(tuán)沉速算法。陳錦山等[8]、Guo和He[9]從長江口絮凝現(xiàn)象的觀測延伸至整個(gè)長江流域。在珠江口,僅有夏小明等[10]使用LISST-ST實(shí)測分析了枯季磨刀門懸沙絮凝沉降特征,得到絮團(tuán)中值粒徑10~96 μm、現(xiàn)場沉速1~20 mm/s;同時(shí)認(rèn)為低懸沙濃度對于絮凝的影響不顯著,而水體污染物的存在可能增大絮凝體沉速。國內(nèi)對于河口區(qū)絮團(tuán)的研究主要集中在長江口,對于珠江口泥沙絮凝現(xiàn)象的研究較少,認(rèn)識還不夠深入,因此有待進(jìn)一步研究。

    珠江是由西江、北江、東江及珠江三角洲諸河匯聚而成的復(fù)合水系,年總徑流量為3.26×1012m3,西江作為珠江流域的最大支流,河流總長為2.21×103km,年徑流量和輸沙率分別為2.28×1011m3和6.45×109t,懸沙濃度為0.23 g/L;北江流域河流長468 km,年徑流量和輸沙率分別為4.49×1010m3和5.45×108t,懸沙濃度為0.13 g/L;東江河流長562 km,年徑流量和輸沙率分別為2.34×1010m3和2.36×108t,懸沙濃度為0.10 g/L。這3條河流總徑流量和輸沙率分別占整個(gè)珠江流域的80%和95%[11]。珠江河口由8個(gè)口門入海,其中磨刀門位于廣東省珠海市洪灣企人石,是西江徑流的主要出??陂T。磨刀門多年平均徑流量為9.23×1010m3,占馬口站徑流量的1/3,占西、北江徑流量的1/3,居珠江三角洲八大口門之首。磨刀門口門潮差小、潮流弱,年均潮差為0.86 m,多年最大潮差平均為2.04 m,在八大口門中是最小的。磨刀門屬強(qiáng)徑流弱潮流河口,多年平均山潮比為5.78,洪水期以徑流作用為主[12]。國內(nèi)學(xué)者已較全面地分析了珠江磨刀門河口懸浮泥沙的來源、路徑、通量以及河口沖淤演變規(guī)律[13—14],在此基礎(chǔ)上,本文利用LISST-100C在珠江磨刀門口外定點(diǎn)觀測洪季大潮和小潮懸浮絮凝及其絮團(tuán)的分布與變化情況。

    圖1 珠江口地形及磨刀門觀測點(diǎn)Fig.1 Map of the Pearl River Estuary topography and the sample site in Modaomen Estuary

    2 資料與方法

    2.1 數(shù)據(jù)來源

    在磨刀門攔門沙外側(cè)布設(shè)垂線進(jìn)行水文測驗(yàn),測點(diǎn)如圖1所示。觀測時(shí)間為2013年7月22日14:00至7月25日14:00(大潮)和7月28日8:00至7月29日11:00(小潮)。測量過程中,采用LISST-100現(xiàn)場激光粒度儀進(jìn)行水體懸浮物現(xiàn)場粒度的逐時(shí)垂線測量,垂線測量采樣間隔為2 s;同時(shí)用OBS-3A與LISST-100同步測量懸沙濃度(NTU)及鹽度,采樣間隔同樣為2 s。另外,每整點(diǎn)時(shí)刻按6層法采取各層水體水樣,用作室內(nèi)含沙量分析及鹽度標(biāo)定,同時(shí)部分特征時(shí)刻如漲急、漲憩、落急、落憩時(shí)刻,額外采集一份水樣,帶回實(shí)驗(yàn)室內(nèi)做懸浮泥沙分散粒度分析。

    2.2 樣品和數(shù)據(jù)處理

    水樣含沙量的分析方法是:用濾膜(孔徑為0.45 μm的醋酸纖維膜,先烘干稱重)進(jìn)行過濾,過濾后的膜在120℃的烘箱內(nèi)烘干8 h,然后稱重,由此重量差和水體體積計(jì)算得出懸浮含沙量。懸浮泥沙分散粒度的分析方法是:先加入0.16%的六偏磷酸鈉(NaPO3)6作為分散劑,浸泡24 h并用超聲波震蕩,用英國馬爾文公司生產(chǎn)的Master sizer 2000型激光粒度儀進(jìn)行粒度分析,獲得1/4Φ(Φ=-log2D)間隔的粒度分布。由粒度分布曲線可計(jì)算平均粒徑等代表參數(shù)。

    OBS濁度和含沙量之間的校正:根據(jù)現(xiàn)場獲取的水樣含沙量,與同一時(shí)刻、同一層位的OBS濁度相比較(見圖2),在濁度小于2 000 NTU的情況下,懸浮物濃度與濁度變化呈線性相關(guān)[15],因此可由OBS的濁度換算得到更多時(shí)間序列的含沙量。

    圖2 濁度和含沙量之間的關(guān)系Fig.2 Relationship between turbidity and suspended sediment concentration

    2.3 沉速計(jì)算方法

    Fettweis提出了現(xiàn)場絮團(tuán)有效密度的計(jì)算公式[16]:

    (1)

    式中,Δρ是有效密度(單位:kg/m3),ρW為水的密度,ρP為泥沙顆粒密度(取2 650 kg/m3),ρF為絮團(tuán)密度(單位:kg/m3),MP可用質(zhì)量濃度代替(單位:kg/m3);VF為體積濃度(單位:μL/L),由LISST-100直接測得。

    沉速計(jì)算根據(jù)Mikkelsen和Pejrup[17]的研究建議,采用Stokes公式:

    (2)

    式中,Dm為絮團(tuán)代表粒徑,此處取為平均粒徑(單位:μm),μ為水的動(dòng)力黏滯系數(shù)(單位:Pa·s),根據(jù)水溫選取相應(yīng)的值。

    3 結(jié)果

    3.1 基本水文泥沙特征

    觀測期間的水文數(shù)據(jù)項(xiàng)目有粒度、流速、含沙量、鹽度和溫度等。由于觀測期間為夏季,水體水溫較高,水溫變化范圍在27.9~29.0℃之間,同時(shí)觀測點(diǎn)水深較淺(約6~8 m),因而溫度在垂向上的變化不大。

    現(xiàn)場的含沙量在大潮和小潮的分布規(guī)律相似,底部的含沙量都大于表層和中層,表層和中層含沙量約在0.05 g/L以下,并且隨著潮周期的變化不顯著;但底層的含沙量平均為0.1 g/L,在大潮和小潮期間均出現(xiàn)多個(gè)峰值達(dá)到0.2 g/L,出現(xiàn)峰值的時(shí)刻接近急流特征時(shí)刻(漲急、落急)。

    分散粒徑隨時(shí)間的變化與特征時(shí)刻相對應(yīng)。泥沙的分散粒徑最大為60.7 μm,最小為8.3 μm;磨刀門口外的泥沙主要為粉砂和砂,黏土成分較少,泥沙顆粒較粗。在大潮時(shí)刻,分散泥沙的平均粒徑變化范圍約為10~140 μm,中層的分散粒徑比表層和底層的大,而表層和底層的分散粒徑較為接近,總體看來,分散粒徑隨著潮汐周期變化呈現(xiàn)出一定的規(guī)律性,在漲急、漲憩時(shí)刻分散粒徑較大,落憩時(shí)刻分散粒徑較小。

    3.2 絮團(tuán)大小及絮凝效率

    根據(jù)LISST-100所測數(shù)據(jù),現(xiàn)場絮團(tuán)平均粒徑基本特征如表1所示。大潮時(shí)最大絮團(tuán)出現(xiàn)在底層,平均粒徑為212.6 μm,最小絮團(tuán)也出現(xiàn)在底層,平均粒徑為13.0 μm,大潮時(shí)絮團(tuán)平均粒徑大小為76.9 μm;小潮時(shí)最大絮團(tuán)出現(xiàn)在表層,平均粒徑為273.8 μm,最小絮團(tuán)出現(xiàn)在底層,平均粒徑為25.9 μm,小潮時(shí)絮團(tuán)平均粒徑大小為131.5 μm。小潮絮團(tuán)平均粒徑大于大潮?,F(xiàn)場絮團(tuán)的粒徑在特征時(shí)刻變化范圍為14.0~212.6 μm;在大潮時(shí)刻表層絮團(tuán)大小多數(shù)大于底層,而在小潮時(shí)刻則相反,底層大多數(shù)絮團(tuán)粒徑大于表層絮團(tuán)。

    大潮時(shí)刻表層和底層絮團(tuán)大小較接近,平均大小約為60 μm,中層絮團(tuán)最大,平均為107 μm;表層和底層泥沙分散粒徑平均約為23 μm,中層泥沙顆粒最粗,平均粒徑為63 μm。絮團(tuán)粒徑與分散粒徑均呈現(xiàn)出中間大,表底小的規(guī)律,在小潮中也是如此。

    總體看來,垂向各水層懸浮物的絮團(tuán)粒徑基本大于其對應(yīng)時(shí)刻的分散粒徑,甚至達(dá)10倍以上,由此可推斷磨刀門河口的懸浮顆粒可能相互黏結(jié)形成絮團(tuán)而使懸浮物粒徑增大,即該區(qū)域發(fā)生明顯的絮凝作用。對應(yīng)特征時(shí)刻,發(fā)現(xiàn)絮團(tuán)粒徑與分散粒徑的分離在憩流時(shí)刻尤為明顯,表明潮周期動(dòng)力環(huán)境變化會(huì)影響磨刀門河口的泥沙絮凝效應(yīng)。此外,對觀測數(shù)據(jù)進(jìn)行垂向?qū)Ρ劝l(fā)現(xiàn),懸浮物絮團(tuán)粒徑與分散粒徑的分離在中層最為明顯,由此向表、底分離均減弱,表明絮凝作用隨水深的增加呈先增后減的變化趨勢。

    表1 磨刀門河口大小潮期間實(shí)測絮團(tuán)平均粒徑的比較Tab.1 Comparison about the mean floc diameter of two tidal cycle

    續(xù)表1

    圖3 磨刀門絮團(tuán)體積濃度、有效密度、沉速隨時(shí)間分布Fig.3 Variation of floc volume concentration,effective floc density and settling velocity along the Modaomen Estuary

    圖4 特征時(shí)刻現(xiàn)場和室內(nèi)粒度頻率曲線Fig.4 Representative accumulative size distributions of dispersed particles and flocs in the estuary

    3.3 絮團(tuán)特征隨時(shí)間變化

    圖3為絮凝體體積濃度、有效密度和沉速隨時(shí)間的變化。絮凝體體積包括絮凝體中水的體積和所有顆粒物的體積。絮凝體體積濃度是指單位體積水體中所含絮凝體的體積,單位為μL/L[6],即體積濃度越大,絮凝體之間的碰撞幾率越大,相互作用也越明顯。大潮時(shí)刻磨刀門水體底層的體積濃度要明顯大于表層和中層,而小潮時(shí)刻表、中、底層體積濃度較為接近。陳錦山等[8]研究了長江流域的絮凝現(xiàn)象指出,絮團(tuán)體積濃度與懸浮物濃度呈良好的線性關(guān)系,即懸浮物濃度高的區(qū)域,絮團(tuán)體積濃度也很高,而絮團(tuán)粒徑反而很小。

    絮凝體的有效密度是絮凝體在水中的真實(shí)密度,是絮團(tuán)的沉速計(jì)算的關(guān)鍵部分。圖中顯示有效密度變化范圍較大,在0~500 g/L之間,且底層的有效密度明顯大于表層和中層,大潮期間的表層和中層有效密度比小潮的略大。

    絮凝體沉速在大小潮時(shí)刻存在差異。大潮時(shí)刻,最大沉速為1.2 mm/s,大部分沉速小于0.8 mm/s,沉速垂向分布從大到小依次為底層、中層、表層,而小潮時(shí)刻底層沉速明顯大于表層和中層,最大可達(dá)1.3 mm/s,表層和中層沉速較接近,均小于0.6 mm/s。

    4 討論

    4.1 絮團(tuán)特性

    4.1.1 絮團(tuán)組成

    圖4a中展示了特征時(shí)刻的粒徑累計(jì)百分比。對比絮團(tuán)的現(xiàn)場觀測和泥沙室內(nèi)分散粒徑的累計(jì)頻率曲線可知,現(xiàn)場粒徑頻率和室內(nèi)粒徑頻率分布差別很大,形成絮團(tuán)后,細(xì)顆粒組分所占頻率明顯減小,而較大顆粒組分所占比例較高,可見在磨刀門口外,水體中的細(xì)顆粒參與了絮凝,并且更多的形成大絮團(tuán)。從垂向分布上看,表層的分散粒徑最細(xì),且現(xiàn)場形成的絮團(tuán)最小,而中層的分散粒徑最粗,絮團(tuán)也大于表層和底層,說明在珠江磨刀門河口泥沙分散粒徑的大小與形成的絮團(tuán)大小有密切關(guān)系。

    4.1.2 體積濃度與絮團(tuán)粒徑

    Van der Lee[18]認(rèn)為,在高體積濃度情況下,絮凝體顆粒增多,加大了絮凝體之間的碰撞幾率,有利于形成大的絮凝體,絮凝體粒徑隨著體積濃度的增大而增大。程江[6]觀測了長江口徐六涇的絮凝現(xiàn)象,認(rèn)為絮凝體平均粒徑和體積濃度的關(guān)系存在兩個(gè)階段,在臨界體積濃度(75 μL/L)之下時(shí)絮凝體平均粒徑隨體積濃度增加而增大,超過臨界值(75 μL/L)時(shí)絮凝體平均粒徑隨體積濃度的增加變化不大。本次的觀測結(jié)果與之略有不同,大、小潮中由于表層和中層的體積濃度基本處于75 μL/L以下,因而呈現(xiàn)出隨著體積濃度的增大而增大的關(guān)系,但是底層卻無明顯的相關(guān)性,這有可能是底部較大含沙量產(chǎn)生再懸浮作用,增加絮團(tuán)間碰撞,因而并未產(chǎn)生隨體積濃度增加而絮團(tuán)增大的現(xiàn)象。

    4.1.3 沉速與絮團(tuán)粒徑

    絮凝體沉速和其平均粒徑之間的關(guān)系可以用ωS=α·Dmβ來表示[17]。從圖5中可以看出,大小潮時(shí)刻的絮團(tuán)沉速均沒有明顯的分層現(xiàn)象,大潮時(shí)刻絮團(tuán)沉速隨著絮團(tuán)粒徑增大而增大,α值在0.000 005~0.000 7之間,β值在1.1~2.3之間,但小潮時(shí)刻則沒有這一明顯關(guān)系。根據(jù)Stokes公式,絮團(tuán)的沉速由有效密度和絮團(tuán)粒徑共同控制,同時(shí)也會(huì)受到其他物理、化學(xué)和生物因素的影響,因而小潮時(shí)刻的特殊性還需進(jìn)一步研究。

    圖5 絮團(tuán)粒徑和體積濃度、沉速的關(guān)系Fig.5 Ralationship between floc diameter,volume concentration and settling velocity

    4.1.4 珠江口絮凝特性

    表2展示了國內(nèi)外河口的絮團(tuán)研究成果,對比可知,珠江口磨刀門絮團(tuán)平均粒徑較大,與國外幾個(gè)河口相接近,變化范圍從10~600 μm不等,略大于長江河口的14~191 μm;同時(shí)絮團(tuán)的有效密度在5~500 kg/m3,比國外河口及長江口小;總體沉速較小與長江口絮團(tuán)沉速相接近。造成這些差異的原因主要與珠江口磨刀門河口的水動(dòng)力特性有關(guān),磨刀門河口含沙量為0.05~0.2 kg/m3,遠(yuǎn)低于長江口的0.1~0.6 kg/m3,并且流速小于長江口(最大0.8 m/s),易形成比長江口大的絮團(tuán),同時(shí)珠江口相比于國外的河口較大,徑流量大,水體交換劇烈,因而產(chǎn)生較低的有效密度,且沉速也相對較小。此外,珠江河口泥沙絮凝現(xiàn)象在不同季節(jié)也有不同的特征,對比夏小明等[10]在枯季的觀測結(jié)果,發(fā)現(xiàn)珠江河口洪季的絮團(tuán)粒徑,沉速等均大于枯季,這可能與徑流量、潮流及溫度等因素的季節(jié)性變化有關(guān)??傮w看來,珠江河口擁有特殊的動(dòng)力條件,其泥沙運(yùn)動(dòng)中絮凝過程是一個(gè)值得深入研究的課題。

    表2 國內(nèi)外河口絮凝研究結(jié)果比較Tab.2 Comparison about the floc from different researchers

    4.2 絮團(tuán)時(shí)空差異及影響因素

    4.2.1 垂向變化

    程江[6]觀測了長江口絮凝現(xiàn)象,發(fā)現(xiàn)絮凝大小在垂線上表現(xiàn)為從表層到底層逐漸增大的規(guī)律。但本次綜合磨刀門河口大、小潮絮凝體在不同特征潮位的垂線變化(圖6)發(fā)現(xiàn),泥沙絮凝體粒徑大小在垂線上基本表現(xiàn)為由表及底先增后減的趨勢,在中層出現(xiàn)最大值,表、底層粒徑相對較小且大小相接近。從鹽度垂向分布圖(圖6)可知,鹽度在垂向上表現(xiàn)為從表層到底層逐漸增大的趨勢,且水體分層較為明顯,小潮時(shí)刻水體層化現(xiàn)象大于大潮時(shí)刻(見圖7)。

    圖6 大小潮絮團(tuán)大小和鹽度的垂向分布Fig.6 Profile of floc diameter and salinity along the flow depth

    圖7 鹽度與絮團(tuán)大小的關(guān)系Fig.7 Relationship between mean floc diameter and salinity

    有對長江口的研究表明[24—26],鹽水會(huì)促進(jìn)絮團(tuán)的形成,鹽度對絮凝體大小的影響分兩個(gè)階段,在較低的鹽度下,鹽度的增加會(huì)促使大絮團(tuán)的產(chǎn)生,而當(dāng)超過一定鹽度之后,對于絮凝體大小的影響不明顯,有學(xué)者在長江口研究得到最佳絮凝鹽度范圍為3~16[6—7,26]。但夏小明等[10]在珠江口的研究中發(fā)現(xiàn)河口鹽度與現(xiàn)場平均粒徑和沉速之間似乎沒有較好的相關(guān)性,并且在不同的河口水域會(huì)呈現(xiàn)出不同的規(guī)律,鹽度的絮凝效應(yīng)具有明顯的區(qū)域變化特征。本次在珠江磨刀門河口的觀測中,研究區(qū)域鹽度較高,大潮鹽度在15~30之間,小潮鹽度在7~30之間,將絮團(tuán)粒徑與鹽度作相關(guān)性分析(見圖7),發(fā)現(xiàn)絮團(tuán)粒徑大小與鹽度沒有明顯相關(guān)關(guān)系,大絮團(tuán)多出現(xiàn)在20~30的鹽度范圍內(nèi),鹽度可能超過了影響絮凝形成的最佳范圍,因此鹽度對絮凝的形成影響較小,說明在珠江磨刀門河口,鹽度并非制約絮團(tuán)形成的主要因素。

    4.2.2 大小潮差異

    水動(dòng)力條件是制約水體中絮凝形成的關(guān)鍵因素。阮文杰[27]對絮凝形成機(jī)理的研究結(jié)果顯示,流速對于絮凝的影響應(yīng)分為兩個(gè)階段:在流速較小時(shí),水體的紊動(dòng)使得水體中泥沙顆粒碰撞增加,利于形成絮團(tuán);當(dāng)流速增加到一定值之后,絮團(tuán)體積也達(dá)到峰值,此時(shí)流速的增加使得絮團(tuán)受到的剪切力增加,絮團(tuán)自身會(huì)被分解,因而會(huì)阻礙大絮團(tuán)的產(chǎn)生?,F(xiàn)場流速和絮團(tuán)大小關(guān)系如圖8,由圖可知,小潮時(shí)刻總體的絮團(tuán)粒徑要大于大潮時(shí)刻,粒徑變化范圍也比大潮時(shí)刻大。流速與絮團(tuán)大小并無明顯相關(guān)性,但在大小潮時(shí)刻流速對絮團(tuán)的影響明顯不同,在小潮時(shí)刻,會(huì)有大于150 μm的絮團(tuán)產(chǎn)生,而在大潮流速較大時(shí),絮團(tuán)反而會(huì)變小,由此可見,在珠江磨刀門河口,大小潮之間的動(dòng)力差異是導(dǎo)致絮團(tuán)大小不同的主導(dǎo)因素,但由于其他因素的影響,流速與絮團(tuán)大小的關(guān)系較為復(fù)雜,絮團(tuán)大小的變化規(guī)律還需進(jìn)一步研究。

    圖8 流速與絮團(tuán)大小的關(guān)系Fig.8 Relationship between mean floc diameter and velocity along the Modaomen Estuary

    5 總結(jié)

    本文根據(jù)洪季珠江口磨刀門口外定點(diǎn)觀測數(shù)據(jù),分析討論了懸沙絮凝特征,得到以下主要結(jié)論:

    (1)珠江口磨刀門懸沙絮凝現(xiàn)象顯著,絮團(tuán)平均粒徑是分散粒徑的2.3倍。絮團(tuán)大小在13.0~273.8 μm之間。小潮絮團(tuán)大于大潮,大潮平均大小為76.9 μm,小潮平均大小為131.5 μm。絮團(tuán)的體積濃度受懸沙濃度影響,垂向分布規(guī)律為底層最大,最高可達(dá)200 μL/L,表層和中層則大多不超過50 μL/L。珠江口形成絮凝的泥沙以粒徑小的細(xì)顆粒為主,并且泥沙分散粒徑對于形成絮團(tuán)的大小有密切關(guān)系;

    (2)通過計(jì)算得到有效密度、沉速的分布規(guī)律,其在空間分布上表現(xiàn)為底層水體較大,中層和表層水體較小,有效密度變化范圍在5~500 g/L,沉速則在0~1.3 mm/s之間,在時(shí)間分布上,大小潮的差異不顯著。通過相關(guān)性分析可知,表層和中層的體積濃度和絮團(tuán)大小有較明顯的正相關(guān)關(guān)系,而底層則由于較高的含沙量使絮團(tuán)碰撞沉降,絮團(tuán)并未隨著體積濃度的增加而增大;沉速和絮團(tuán)粒徑大小的關(guān)系在大小潮不同,大潮時(shí)沉速隨粒徑增大而增大,但小潮時(shí)刻沒有明顯的相關(guān)性;

    (3)通過分析鹽度和流速對于絮團(tuán)大小的影響可知,在珠江口高鹽度(鹽度10~30)的環(huán)境下,鹽度對于絮團(tuán)大小的影響很小;絮團(tuán)粒徑在大小潮之間的差異主要由流速差異引起,大潮時(shí)刻流速較大,難以形成大于150 μm的絮團(tuán),而在流速較小的小潮時(shí)刻,可以形成大于150 μm的大絮團(tuán),由此可見流速對于絮團(tuán)大小的影響分兩個(gè)階段,當(dāng)流速較小時(shí)促進(jìn)絮團(tuán)的形成,當(dāng)流速較大時(shí)則會(huì)破壞絮團(tuán),抑制大絮團(tuán)的形成。

    致謝:本研究現(xiàn)場采樣得到中山大學(xué)河口海岸研究所的協(xié)助,特此感謝!

    [1] Manning A J, Langston W J, Jonas P J C. A review of sediment dynamics in the Severn Estuary: Influence of flocculation[J]. Marine Pollution Bulletin, 2010, 61(1/3): 37-51.

    [2] Manning A J, Dyer K R. The use of optics for the in situ determination of flocculated mud characteristics[J]. Journal of Optics a-Pure and Applied Optics, 2002, 4(4): S71-S81.

    [3] Winterwerp J C, Bale A J, Christie M C, et al. Flocculation and settling velocity of fine sediment[Z]. Elsevier, 2002,5: 25-40.

    [4] Gibbs R J. Estuarine flocs: their size, settling velocity and density[J]. Journal of Geophysical Research, 1985, 90(C2): 3249-3251.

    [5] Agrawal Y C, Pottsmith H C, Lynch J, et al. Laser instruments for particle size and settling velocity measurements in the coastal zone[G]// Prospects for the 21st Century Conference Proceedings. OCEANS 96 MTS/IEEE (Cat. No.96CH35967), 1996(3): 1135-1142.

    [6] 程江, 何青, 王元葉. 利用LISST觀測絮凝體粒徑、有效密度和沉速的垂線分布[J]. 泥沙研究, 2005(1):33-39.

    Cheng Jiang,He Qing,Wang Yuanye. Using LISST 100 for in situ estimates of floc size, density and settling velocity, Changjiang Estuary, China[J].Journal of Sediment Research, 2005(1):33-39.

    [7] 唐建華, 何青, 王元葉, 等. 長江口渾濁帶絮凝體特性[J]. 泥沙研究, 2008(2):27-33.

    Tang Jianhua, He Qing, Wang Yuanye, et al. Study on in-situ flocs size in turbidity maximum of the Changjiang Estuary[J].Journal of Sediment Research, 2008(2):27-33.

    [8] 陳錦山,何青,郭磊城. 長江懸浮物絮凝特征[J]. 泥沙研究, 2011(5): 11-18.

    Chen Jinshan, He Qing, Guo Leicheng. Flocculation characteristics of suspended particulate matter in Yangtze River[J].Journal of Sediment Research, 2011(5): 11-18.

    [9] Guo L C, He Q. Freshwater flocculation of suspended sediments in the Yangtze River, China[J]. Ocean Dynamics, 2011, 61(2/3): 371-386.

    [10] 夏小明,李炎,楊輝,等. 枯季珠江河口懸浮泥沙絮凝沉降特征的觀測與分析[J]. 海洋學(xué)研究, 2006, 24(1): 6-18.

    Xia Xiaoming, Li Yan, Yang Hui, et al. Observations of the size and settling velocity distributions of suspended sediment in the Zhujiang River Estuary during dry season[J]. Journal of Marine Sciences, 2006, 24(1): 6-18.

    [11] 吳創(chuàng)收. 華南流域人類活動(dòng)和氣候變化對入海水沙通量和三角洲演化的影響[D]. 上海:華東師范大學(xué),2013.

    Wu Chuangshou. Impacts of human activities and climate change on sediment flux and delta evolution in Southern China basin[D]. Shanghai:East China Normal University, 2013.

    [12] 譚超,黃本勝,龔文平,等. 珠江磨刀門河口排洪動(dòng)力特征及攔門沙演變響應(yīng)的初步研究[J]. 水利學(xué)報(bào), 2013(9): 1023-1029.

    Tan Chao, Huang Bensheng, Gong Wenping, et al. The dynamics of flood releasing and its effects on the morphological evolution of the mouth bar in the Modaomen Estuary of the Pearl River[J]. Shuili Xuebao, 2013(9): 1023-1029.

    [13] 賈良文,吳超羽,任杰,等. 珠江口磨刀門枯季水文特征及河口動(dòng)力過程[J]. 水科學(xué)進(jìn)展, 2006(1): 82-88.

    Jia Liangwen, Wu Chaoyu, Ren Jie, et al. Hydrologic characteristics and estuarine dynamic process during the dry season in Modaomen Estuary of the Pearl River[J]. Advances in Water Science, 2006(1): 82-88.

    [14] 方神光,王少波. 磨刀門水道枯季水動(dòng)力特性分析[J]. 水文, 2013(5): 70-74.

    Fang Shenguang, Wang Shaobo. Analysis of water dynamic characteristics of Modaomen waterway during dry season[J]. Journal of China Hydrology, 2013(5): 70-74.

    [15] 劉紅,何青,王元葉,等. 長江口渾濁帶海域OBS標(biāo)定的實(shí)驗(yàn)研究[J]. 泥沙研究, 2006(5): 52-58.

    Liu Hong, He Qing, Wang Yuanye, et al. OBS situ calibration research in the turbidity maximum of the Changjiang Estuary, China[J]. Journal of Sediment Research, 2006(5): 52-58.

    [16] Fettweis M. Uncertainty of excess density and settling velocity of mud flocs derived from in situ measurements[J]. Estuarine Coastal and Shelf Science, 2008, 78(2): 426-436.

    [17] Mikkelsen O A, Pejrup M. The use of a LISST-100 laser particle sizer for in-situ estimates of flee size, density and settling velocity[J]. Geo-Marine Letters, 2001, 20(4): 187-195.

    [18] Van der Lee W T B. Temporal variation of floc size and settling velocity in the Dollard estuary[J]. Continental Shelf Research, 2000, 20(12/13): 1495-1511.

    [19] Dyer K R, Bale A J, Christie M C, et al. The turbidity maximum in a Mesotidal Estuary, the Tamar Estuary, UK: Ⅱ. the floc properties[J]. Proceedings in Marine Science, 2002, 5: 219-232.

    [20] Dyer K R, Christie M C, Feates N, et al. An investigation into processes influencing the morphodynamics of an intertidal mudflat, the Dollard Estuary, the Netherlands: I. Hydrodynamics and suspended sediment[J]. Estuarine, Coastal and Shelf Science, 2000, 50(5): 607-625.

    [21] Manning A J, Dyer K R, Lafite R, et al. Flocculation measured by video based instruments in the Gironde Estuary during the European Commission SWAMIEE project[J]. Journal of Coastal Research, 2004(5): 58-69.

    [22] Manning A J, Martens C, De Mulder T, et al. Mud floc observations in the turbidity maximum zone of the Scheldt estuary during neap tides[J]. J Coast Res SI, 2007, 50: 832-836.

    [23] Manning A J, Schoellhamer D H. Factors controlling floc settling velocity along a longitudinal estuarine transect[J]. Marine Geology, 2013, 345: 266-280.

    [24] 關(guān)許為,陳英祖,杜心慧. 長江口絮凝機(jī)理的試驗(yàn)研究[J]. 水利學(xué)報(bào), 1996(6): 70-74.

    Guan Xuwei, Chen Yingzu, Du Xinhui. Experimental study on mechanism of flocculation in Yangtze Estuary[J]. Shuili Xuebao, 1996(6): 70-74.

    [25] 林以安,唐仁友,李炎. 長江口絮凝聚沉特征與顆粒表面理化因素作用——Ⅱ.顆粒表面性質(zhì)對聚沉的作用[J]. 泥沙研究, 1997(4): 78-86.

    Lin Yian, Tang Renyou, Li Yan. Coagulation, settling characteristics and action of surface physic-chemical factors of particles in the Changjiang Estuary——Ⅱ. Role of surface colloidal property variation on the coagulation and settling of particles[J]. Journal of Sediment Research, 1997(4): 78-86.

    [26] 張志忠,蔣國俊. 長江口陽離子濃度與細(xì)顆粒泥沙絮凝沉積[J]. 海洋學(xué)報(bào), 1995, 17(1): 76-82.

    Zhang Zhizhong, Jiang Guojun. Cation concentration and fine sediment flocculation deposition in Changjiang Estuary[J].Haiyang Xuebao, 1995, 17(1): 76-82.

    [27] 阮文杰. 細(xì)顆粒泥沙動(dòng)水絮凝的機(jī)理分析[J]. 海洋科學(xué), 1991(5): 46-49.

    Ruan Wenjie. Mechanistic analysis on hydrodynamic flocculation on fine sediments[J]. Marine Sciences, 1991(5): 46-49.

    Observations of in situ flocs characteristic in the Modaomen Estuary of the Pearl River

    Deng Zhirui1,He Qing1,Yang Qingshu2,Lin Jianliang1

    (1.StateKeyLaboratoryofEstuarineandCoastalResearch,EastChinaNormalUniversity,Shanghai200062,China;2.InstituteofEstuarineandCoastalResearch,SunYat-SenUniversity,Guangzhou510275,China)

    The Laser In Situ Scattering and Transmissometry-100(LISST-100)was used to the field observation in the Modaomen Estuary of the Pearl River during the July 2013(flood season). Based on the in situ data,we get some flocculation information: the dispersion of suspended sediment particle size is about 27.9 μm,the measured floccules mean diameter is about 91.6 μm,it indicated that suspended sediment flocculation phenomenon significantly. The in situ mean diameter of floc varied from 13.0 to 273.8 μm,the neap tidal floc was larger than spring tide. The variation of floc particle size was decreases at first,and then increases. The relationship between floc diameter and volume concentration,settling velocity were vary from case to case. The volume concentration and particle size of the floccules had positive correlation in surface and middle layer,floccules settling velocity increases with the increase of particle size in spring time. Comprehensive analysis of the factors affecting flocculation,the results show that the salinity have no significant relationship with floccules size,the influence of velocity to floccules had two kinds of circumstances,large flow motion will hinder the development of the floc size. These knowledge will contributes to a further understanding of fine sediment transport characteristics and biochemical process in the Pearl River Estuary.

    Pearl River Estuary; Modaomen Estuary; suspended sediment; flocculation

    2015-01-22;

    2015-04-09。

    國家自然科學(xué)基金項(xiàng)目(41276080,41476073,51320105005);科技部科技基礎(chǔ)性工作專項(xiàng)重點(diǎn)項(xiàng)目“中國典型河口動(dòng)力沉積地貌本底數(shù)據(jù)調(diào)查”(2013FY112000)。

    鄧智瑞(1990—),男,廣西省來賓市人,博士研究生,主要從事河口海岸泥沙研究。E-mail:52132601012@ecnu.cn

    *通信作者:何青,教授,主要從事河口海岸水沙運(yùn)動(dòng)研究工作。E-mail:qinghe@sklec.ecnu.edu.cn

    10.3969/j.issn.0253-4193.2015.09.015

    P737.14

    A

    0253-4193(2015)09-0152-10

    鄧智瑞,何青,楊清書,等. 珠江口磨刀門泥沙絮凝特征[J]. 海洋學(xué)報(bào),2015,37(9):152-161,

    Deng Zhirui,He Qing,Yang Qingshu,et al. Observations of in situ flocs characteristic in the Modaomen Estuary of the Pearl River[J]. Haiyang Xuebao,2015,37(9):152-161,doi:10.3969/j.issn.0253-4193.2015.09.015

    猜你喜歡
    小潮絮團(tuán)大潮
    第五課:中國青年投身改革開放的大潮
    基于FBRM和PVM技術(shù)的尾礦濃密過程絮團(tuán)演化規(guī)律
    絮團(tuán)稠化對全尾砂濃密性能的影響
    5G手機(jī)大潮即來
    海峽姐妹(2019年1期)2019-03-23 02:42:54
    希 望
    希望
    新一季流行色已發(fā)布?快來入手同色系數(shù)碼3C小潮物!
    一種改進(jìn)的聚合模型在污泥絮凝-沉降模擬中的應(yīng)用
    《神威啟示錄》系列報(bào)道四 大潮退去看“裸泳”
    測量絮團(tuán)分維數(shù)的新方法
    蜜桃在线观看..| 国产精品久久久久久av不卡| 男女啪啪激烈高潮av片| 看免费成人av毛片| 久久精品国产亚洲av天美| 亚洲国产av新网站| 精品少妇黑人巨大在线播放| av在线app专区| 插阴视频在线观看视频| 自拍欧美九色日韩亚洲蝌蚪91 | 高清视频免费观看一区二区| 国产淫片久久久久久久久| 成人毛片60女人毛片免费| 久久av网站| 亚洲av在线观看美女高潮| 51国产日韩欧美| 久久久精品94久久精品| 欧美日韩精品成人综合77777| 在线天堂最新版资源| 日韩视频在线欧美| 性色av一级| 亚洲av男天堂| 插逼视频在线观看| 中文字幕制服av| 天美传媒精品一区二区| 亚洲精品亚洲一区二区| 午夜免费观看性视频| 水蜜桃什么品种好| 一级毛片久久久久久久久女| 汤姆久久久久久久影院中文字幕| 欧美人与善性xxx| 少妇 在线观看| 我要看日韩黄色一级片| 国产爱豆传媒在线观看| 精品视频人人做人人爽| 亚洲精品日本国产第一区| 在线观看三级黄色| 欧美日韩综合久久久久久| 最近2019中文字幕mv第一页| 亚洲av国产av综合av卡| 全区人妻精品视频| 高清视频免费观看一区二区| 日韩电影二区| 丝瓜视频免费看黄片| 热re99久久精品国产66热6| 成人18禁高潮啪啪吃奶动态图 | 在线观看美女被高潮喷水网站| 99热网站在线观看| 国产精品成人在线| 1000部很黄的大片| 久久女婷五月综合色啪小说| 一个人看视频在线观看www免费| 久久国产精品男人的天堂亚洲 | 久久精品久久久久久久性| 最后的刺客免费高清国语| 亚洲高清免费不卡视频| 国产视频首页在线观看| 欧美日韩视频高清一区二区三区二| 国产精品人妻久久久影院| 自拍欧美九色日韩亚洲蝌蚪91 | 久久女婷五月综合色啪小说| 在线观看免费视频网站a站| 老熟女久久久| 韩国高清视频一区二区三区| 高清欧美精品videossex| 欧美zozozo另类| 国产精品欧美亚洲77777| 精品一区二区免费观看| 亚洲av国产av综合av卡| 永久网站在线| 亚洲经典国产精华液单| 亚洲国产精品999| 男人添女人高潮全过程视频| 亚洲av综合色区一区| 美女内射精品一级片tv| 在线观看美女被高潮喷水网站| 高清av免费在线| 在线观看一区二区三区| freevideosex欧美| av.在线天堂| 伊人久久精品亚洲午夜| 亚洲av成人精品一区久久| 午夜精品国产一区二区电影| 最近的中文字幕免费完整| 少妇猛男粗大的猛烈进出视频| 天堂8中文在线网| av又黄又爽大尺度在线免费看| 欧美97在线视频| 亚洲av不卡在线观看| 国产精品99久久99久久久不卡 | 亚洲一区二区三区欧美精品| 国产久久久一区二区三区| 97热精品久久久久久| 国产伦精品一区二区三区视频9| 国产精品成人在线| 身体一侧抽搐| 欧美变态另类bdsm刘玥| 亚洲精品一区蜜桃| 国产黄色视频一区二区在线观看| 六月丁香七月| 天美传媒精品一区二区| 欧美另类一区| 精品午夜福利在线看| 欧美bdsm另类| 国产极品天堂在线| 成人毛片60女人毛片免费| 成年美女黄网站色视频大全免费 | 各种免费的搞黄视频| 少妇高潮的动态图| 日本欧美视频一区| 超碰97精品在线观看| 大片免费播放器 马上看| 麻豆国产97在线/欧美| 日韩中字成人| 精品国产露脸久久av麻豆| 一级毛片aaaaaa免费看小| 国产精品蜜桃在线观看| 在线 av 中文字幕| 欧美一区二区亚洲| 欧美bdsm另类| 午夜福利在线观看免费完整高清在| 亚洲精品一二三| 国产在线一区二区三区精| 高清在线视频一区二区三区| 三级国产精品片| 啦啦啦啦在线视频资源| 中国国产av一级| 国产亚洲最大av| 久久久精品94久久精品| 国产亚洲午夜精品一区二区久久| 97在线视频观看| 激情 狠狠 欧美| 九色成人免费人妻av| 一级片'在线观看视频| 国产亚洲欧美精品永久| 麻豆国产97在线/欧美| 亚洲国产欧美人成| 色哟哟·www| 国产精品三级大全| 我要看日韩黄色一级片| 成人18禁高潮啪啪吃奶动态图 | 精品人妻熟女av久视频| 国产 精品1| 简卡轻食公司| 午夜激情福利司机影院| 久久人人爽人人片av| 免费在线观看成人毛片| 最近2019中文字幕mv第一页| 国产片特级美女逼逼视频| 一本一本综合久久| 乱系列少妇在线播放| 51国产日韩欧美| 国产乱人偷精品视频| 国产淫语在线视频| 亚洲成色77777| 男人舔奶头视频| 精品久久久久久电影网| 亚洲国产精品专区欧美| 国产亚洲91精品色在线| 有码 亚洲区| 啦啦啦中文免费视频观看日本| 最近手机中文字幕大全| 男人爽女人下面视频在线观看| 国产一级毛片在线| 国产精品久久久久久精品古装| 欧美成人精品欧美一级黄| 如何舔出高潮| 亚洲不卡免费看| 亚洲精品中文字幕在线视频 | 精品少妇久久久久久888优播| 噜噜噜噜噜久久久久久91| 青青草视频在线视频观看| 最近中文字幕高清免费大全6| 少妇裸体淫交视频免费看高清| 91久久精品国产一区二区成人| av卡一久久| 国产在视频线精品| av卡一久久| 在线观看三级黄色| 日本色播在线视频| 国产精品一区二区性色av| 蜜臀久久99精品久久宅男| 欧美精品一区二区免费开放| 亚洲精品日韩av片在线观看| 少妇人妻 视频| 中国三级夫妇交换| 国产精品三级大全| 校园人妻丝袜中文字幕| 婷婷色综合www| 国产美女午夜福利| 欧美xxxx性猛交bbbb| 国产v大片淫在线免费观看| 国产乱人偷精品视频| 免费人成在线观看视频色| 国产色婷婷99| 久久久欧美国产精品| 精品酒店卫生间| 制服丝袜香蕉在线| 久久热精品热| 精品久久久久久久久亚洲| 一本—道久久a久久精品蜜桃钙片| 麻豆国产97在线/欧美| 色视频在线一区二区三区| 成人无遮挡网站| 成年人午夜在线观看视频| 热re99久久精品国产66热6| 精品亚洲乱码少妇综合久久| 午夜免费观看性视频| 麻豆成人av视频| 亚洲成人av在线免费| 国产黄频视频在线观看| 国产精品av视频在线免费观看| 少妇的逼好多水| 深爱激情五月婷婷| videossex国产| 狂野欧美白嫩少妇大欣赏| 国产精品久久久久久精品古装| 乱码一卡2卡4卡精品| 精品久久久噜噜| 国产精品99久久久久久久久| 色网站视频免费| 午夜福利网站1000一区二区三区| 亚洲av国产av综合av卡| 精品国产三级普通话版| 精品久久国产蜜桃| 18禁裸乳无遮挡免费网站照片| 麻豆精品久久久久久蜜桃| 亚洲第一av免费看| 一个人看视频在线观看www免费| 寂寞人妻少妇视频99o| av免费观看日本| 一边亲一边摸免费视频| 天美传媒精品一区二区| 男人和女人高潮做爰伦理| 免费看av在线观看网站| 免费黄频网站在线观看国产| 人人妻人人看人人澡| 大香蕉久久网| 边亲边吃奶的免费视频| 国产一区二区在线观看日韩| 亚洲av福利一区| 欧美日韩视频高清一区二区三区二| 成人无遮挡网站| 韩国av在线不卡| 蜜桃久久精品国产亚洲av| 国产在线男女| 国产精品.久久久| 亚洲精品乱码久久久久久按摩| 免费播放大片免费观看视频在线观看| 国产国拍精品亚洲av在线观看| av国产精品久久久久影院| 麻豆成人午夜福利视频| 狠狠精品人妻久久久久久综合| 男人添女人高潮全过程视频| 国产在线免费精品| 中文字幕av成人在线电影| 国产伦精品一区二区三区视频9| 国产精品免费大片| 91久久精品国产一区二区成人| 国产美女午夜福利| 亚洲不卡免费看| 欧美老熟妇乱子伦牲交| 美女脱内裤让男人舔精品视频| 91aial.com中文字幕在线观看| 日韩三级伦理在线观看| 亚洲,一卡二卡三卡| 欧美xxxx性猛交bbbb| 丝瓜视频免费看黄片| 国产又色又爽无遮挡免| 在线天堂最新版资源| 五月伊人婷婷丁香| 丰满乱子伦码专区| 性色avwww在线观看| 91狼人影院| 久久精品国产a三级三级三级| 又爽又黄a免费视频| 99re6热这里在线精品视频| 高清黄色对白视频在线免费看 | 国产在线免费精品| 久久久久久九九精品二区国产| 欧美日韩精品成人综合77777| 一级毛片电影观看| 高清欧美精品videossex| 亚洲av不卡在线观看| 免费人成在线观看视频色| 久久国产精品大桥未久av | 身体一侧抽搐| 日本vs欧美在线观看视频 | av免费在线看不卡| 91在线精品国自产拍蜜月| 亚洲av综合色区一区| 亚洲在久久综合| 久久99蜜桃精品久久| 高清不卡的av网站| 中文字幕制服av| 黄色一级大片看看| 菩萨蛮人人尽说江南好唐韦庄| 精品久久久精品久久久| 国产69精品久久久久777片| 亚洲人成网站高清观看| 欧美日本视频| 精品国产露脸久久av麻豆| 免费播放大片免费观看视频在线观看| 欧美性感艳星| 日韩 亚洲 欧美在线| 免费黄频网站在线观看国产| 十八禁网站网址无遮挡 | 你懂的网址亚洲精品在线观看| 五月伊人婷婷丁香| 偷拍熟女少妇极品色| 丝袜脚勾引网站| 日韩视频在线欧美| 欧美日韩一区二区视频在线观看视频在线| 亚洲欧洲国产日韩| av在线播放精品| 热re99久久精品国产66热6| h视频一区二区三区| 在线看a的网站| 日韩制服骚丝袜av| 午夜免费鲁丝| 热re99久久精品国产66热6| 亚洲欧美成人精品一区二区| 天天躁夜夜躁狠狠久久av| 18禁裸乳无遮挡动漫免费视频| av卡一久久| 99热这里只有精品一区| 一个人看视频在线观看www免费| 久久久午夜欧美精品| 九色成人免费人妻av| 最近手机中文字幕大全| 精品人妻一区二区三区麻豆| h日本视频在线播放| 亚洲欧洲国产日韩| 啦啦啦中文免费视频观看日本| 妹子高潮喷水视频| a级一级毛片免费在线观看| 午夜福利网站1000一区二区三区| 国产精品一及| 一级片'在线观看视频| 日韩制服骚丝袜av| 美女内射精品一级片tv| 欧美成人午夜免费资源| 看十八女毛片水多多多| 视频区图区小说| 女人十人毛片免费观看3o分钟| 国产精品久久久久久久久免| 高清视频免费观看一区二区| 国产成人免费观看mmmm| 午夜福利影视在线免费观看| 午夜精品国产一区二区电影| 欧美人与善性xxx| 一区在线观看完整版| 一区二区三区四区激情视频| 国产一级毛片在线| 看免费成人av毛片| 少妇丰满av| 欧美成人a在线观看| 纯流量卡能插随身wifi吗| 伊人久久精品亚洲午夜| 91aial.com中文字幕在线观看| 校园人妻丝袜中文字幕| 中文字幕亚洲精品专区| 久久鲁丝午夜福利片| 国产精品偷伦视频观看了| 天美传媒精品一区二区| 美女高潮的动态| 亚洲欧美日韩东京热| 视频区图区小说| 能在线免费看毛片的网站| 亚洲av在线观看美女高潮| 欧美区成人在线视频| 七月丁香在线播放| 亚洲四区av| 大又大粗又爽又黄少妇毛片口| 天堂中文最新版在线下载| 久久人人爽人人爽人人片va| 国产 精品1| 日韩成人伦理影院| 日本av免费视频播放| 校园人妻丝袜中文字幕| 国产爽快片一区二区三区| av福利片在线观看| 色吧在线观看| av女优亚洲男人天堂| 亚洲精品一二三| 3wmmmm亚洲av在线观看| 国产成人免费无遮挡视频| 高清在线视频一区二区三区| 国产精品久久久久久久电影| 日韩中文字幕视频在线看片 | 大陆偷拍与自拍| 国产精品人妻久久久影院| 欧美3d第一页| 18+在线观看网站| kizo精华| 日韩欧美 国产精品| 国产一区亚洲一区在线观看| 午夜福利视频精品| 亚洲精品乱码久久久v下载方式| 看十八女毛片水多多多| 黄片wwwwww| 在线观看美女被高潮喷水网站| 国产乱来视频区| 久久这里有精品视频免费| 国产乱人视频| 国产精品人妻久久久影院| 久久久久网色| 久久这里有精品视频免费| 交换朋友夫妻互换小说| xxx大片免费视频| av女优亚洲男人天堂| 亚洲av中文字字幕乱码综合| 22中文网久久字幕| 国产男女内射视频| 少妇人妻精品综合一区二区| 亚洲精华国产精华液的使用体验| 亚洲av中文av极速乱| 搡老乐熟女国产| 国产精品久久久久久精品古装| 成人高潮视频无遮挡免费网站| av播播在线观看一区| av在线播放精品| 久久精品国产a三级三级三级| 插逼视频在线观看| 国产精品爽爽va在线观看网站| 超碰97精品在线观看| 99精国产麻豆久久婷婷| 多毛熟女@视频| 亚洲av综合色区一区| 成人亚洲精品一区在线观看 | 亚洲,一卡二卡三卡| 国产乱人偷精品视频| 中文天堂在线官网| 狂野欧美激情性xxxx在线观看| 超碰97精品在线观看| 日韩,欧美,国产一区二区三区| 亚洲av.av天堂| 简卡轻食公司| 91狼人影院| 精品国产三级普通话版| 王馨瑶露胸无遮挡在线观看| 少妇精品久久久久久久| 亚洲一区二区三区欧美精品| 国产av国产精品国产| 亚洲av在线观看美女高潮| 日韩欧美 国产精品| 综合色丁香网| 久久久久久久亚洲中文字幕| 亚洲精品久久久久久婷婷小说| av在线播放精品| 精品久久久噜噜| 免费人妻精品一区二区三区视频| 欧美成人a在线观看| 日韩欧美一区视频在线观看 | 国产黄色免费在线视频| 免费观看无遮挡的男女| 国产色婷婷99| 精品人妻视频免费看| 国内精品宾馆在线| 国产毛片在线视频| 18禁裸乳无遮挡动漫免费视频| 国产精品国产三级国产专区5o| 黄色日韩在线| a级毛片免费高清观看在线播放| 国产爽快片一区二区三区| 国内少妇人妻偷人精品xxx网站| 六月丁香七月| 国产在线视频一区二区| 国产亚洲欧美精品永久| 欧美3d第一页| 国产免费一级a男人的天堂| 人人妻人人添人人爽欧美一区卜 | 99视频精品全部免费 在线| 久久99蜜桃精品久久| 国产 精品1| 亚洲色图av天堂| 日韩av免费高清视频| 内地一区二区视频在线| 亚洲精品乱码久久久久久按摩| 精品久久久久久久末码| 国产精品人妻久久久久久| 草草在线视频免费看| 狂野欧美白嫩少妇大欣赏| 夫妻性生交免费视频一级片| 亚洲av电影在线观看一区二区三区| 国产 一区 欧美 日韩| 91aial.com中文字幕在线观看| 女性被躁到高潮视频| 日韩成人伦理影院| 高清日韩中文字幕在线| 国产亚洲一区二区精品| 男女免费视频国产| 美女福利国产在线 | 久久久久久久久久人人人人人人| 99九九线精品视频在线观看视频| 少妇精品久久久久久久| 少妇被粗大猛烈的视频| 精品人妻熟女av久视频| 网址你懂的国产日韩在线| 中国国产av一级| 久久国产精品男人的天堂亚洲 | 天天躁夜夜躁狠狠久久av| 一区二区三区免费毛片| 伦精品一区二区三区| 国产精品三级大全| videos熟女内射| 日韩中字成人| 秋霞在线观看毛片| 午夜福利在线观看免费完整高清在| 国精品久久久久久国模美| 久热这里只有精品99| 18禁裸乳无遮挡动漫免费视频| 国产乱人偷精品视频| 在线观看免费日韩欧美大片 | www.av在线官网国产| 亚洲欧洲日产国产| 久久精品国产a三级三级三级| 亚洲精品自拍成人| 欧美成人一区二区免费高清观看| 在线播放无遮挡| 美女国产视频在线观看| 国产亚洲一区二区精品| 国产av码专区亚洲av| 中国国产av一级| 免费少妇av软件| 久久久久精品久久久久真实原创| 最近的中文字幕免费完整| 成人亚洲欧美一区二区av| 噜噜噜噜噜久久久久久91| 精品视频人人做人人爽| 亚洲一级一片aⅴ在线观看| 国产日韩欧美亚洲二区| 99久久精品热视频| 美女视频免费永久观看网站| 日韩一本色道免费dvd| av不卡在线播放| 波野结衣二区三区在线| 国产精品女同一区二区软件| 午夜免费观看性视频| 国产精品国产三级国产av玫瑰| 九九爱精品视频在线观看| 深夜a级毛片| 青春草国产在线视频| 午夜福利高清视频| 亚洲av免费高清在线观看| 欧美日韩综合久久久久久| 91精品国产国语对白视频| 欧美人与善性xxx| 97在线视频观看| 联通29元200g的流量卡| 久久久久视频综合| 尤物成人国产欧美一区二区三区| 美女主播在线视频| 欧美3d第一页| 久久99精品国语久久久| 亚洲人与动物交配视频| 国产女主播在线喷水免费视频网站| 亚洲欧美清纯卡通| 成人免费观看视频高清| 成人18禁高潮啪啪吃奶动态图 | 欧美xxⅹ黑人| 能在线免费看毛片的网站| 免费大片黄手机在线观看| 欧美性感艳星| 赤兔流量卡办理| 亚洲欧洲日产国产| h视频一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91 | 久久精品国产亚洲av天美| 精华霜和精华液先用哪个| 妹子高潮喷水视频| 国产黄频视频在线观看| 一级黄片播放器| 日韩av在线免费看完整版不卡| 黄色视频在线播放观看不卡| 日韩人妻高清精品专区| 男人爽女人下面视频在线观看| 永久网站在线| 毛片女人毛片| 久久人妻熟女aⅴ| 中文乱码字字幕精品一区二区三区| 成人毛片a级毛片在线播放| 97精品久久久久久久久久精品| 久久99蜜桃精品久久| 香蕉精品网在线| 亚洲伊人久久精品综合| 欧美一级a爱片免费观看看| 国产男女内射视频| av专区在线播放| 成人特级av手机在线观看| 免费看不卡的av| 久久综合国产亚洲精品| 国产伦在线观看视频一区| 色视频在线一区二区三区| 不卡视频在线观看欧美| 亚洲欧美日韩卡通动漫| 三级经典国产精品| 成人无遮挡网站| 欧美三级亚洲精品| av女优亚洲男人天堂| 男女边摸边吃奶| 青春草国产在线视频| 欧美日韩在线观看h| 亚洲av不卡在线观看| 亚洲怡红院男人天堂| videos熟女内射| 精品午夜福利在线看| 秋霞伦理黄片| 亚洲国产精品一区三区| 免费看光身美女| 亚洲,一卡二卡三卡| 老司机影院毛片| 国内精品宾馆在线| 国产爱豆传媒在线观看| 日韩一区二区三区影片|