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

    現(xiàn)代黃河三角洲沉積物粒度特征及其來源

    2024-01-01 00:00:00于霽懷介冬梅李平
    關(guān)鍵詞:端元黃河三角洲粒度

    摘要:現(xiàn)代黃河三角洲位于海陸交互地帶,為我國重要的生態(tài)保護(hù)區(qū),同時也是石油開采地。受河流、潮流、人類活動影響,現(xiàn)代黃河三角洲沉積物粒度組成復(fù)雜,物質(zhì)來源多樣。本文共采集了黃河三角洲表層沉積物樣品219個,首先采用粒度分析方法分析了現(xiàn)代黃河三角洲沉積物的粒度特征以及沉積物來源,然后采用參數(shù)端元模型分析方法對沉積物的物質(zhì)來源進(jìn)行了劃分,最后采用聚類分析方法分析了表層沉積物類型。結(jié)果表明:現(xiàn)代黃河三角洲沉積物粒度組成以粉砂(體積分?jǐn)?shù)71.10%)和砂(27.62%)為主,大部分沉積物粒度分選中等,粒度分布曲線呈近對稱分布。沉積物劃分為EM1、EM2、EM3三個沉積端元,分別代表三種不同物質(zhì)來源,其中EM1代表舊河流沉積作用,EM2代表海洋潮流沉積作用,EM3代表現(xiàn)黃河河流沉積作用。將表層沉積物分成三種沉積物類型,結(jié)合聚類分析結(jié)果將黃河三角洲分為三類沉積環(huán)境,其中:第一類主要分布在黃河三角洲等西部河流的中上游,包括馬新河、沾利河、草橋溝、挑河,主要為舊河流沉積,水動力較弱;第二類主要分布在刁口河、神仙溝、現(xiàn)黃河(清水溝)等東部河流,現(xiàn)黃河河流沉積作用較強;第三類主要分布在黃河三角洲北部沿岸及潮灘,受海洋潮流作用較強,海岸受到侵蝕。此外,人類活動對黃河三角洲岸線變化也產(chǎn)生了一定影響。

    關(guān)鍵詞:黃河三角洲;表層沉積物;粒度;端元;物源分析

    doi:10.13278/j.cnki.jjuese.2023020

    中圖分類號:P736.21

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

    Supported by the National Natural Science Foundation of China (42271162,41971100),the Natural Science Foundation of Shandong Province (ZR2020MD063) and the Joint Fund Key Support Project of NSFCShandong Province (U1706214)

    Grain Size Characteristic and Source of Sediments

    in Modern Yellow River Delta

    Yu Jihuai1,Jie Dongmei 3,4,5,Li Ping6,7

    1. School of Geographical Sciences, Northeast Normal University, Changchun 130024, China

    2. Institute for Peat and Mire Research, Northeast Normal University, Changchun 130024, China

    3. State Environmental Protection Key Laboratory of Wetland Ecology and Vegetation Restoration, Changchun 130024, China

    4. Key Laboratory of Geographical Processes and Ecological Security of Changbai Mountains (Northeast Normal University), "Ministry of Education, Changchun 130024, China

    5. Key Laboratory of Vegetation Ecology, Institute of Grassland (Northeast Normal University), Ministry of Education, "Changchun 130024, China

    6. First Institute of Oceanography, Ministry of Natural Resources, Qingdao 266061, Shandong,China

    7. Key Laboratory of Coastal Zone Science and Integrated Management, Ministry of Natural Resources, Qingdao 266061, "Shandong,China

    Abstract: There are some important ecological reserves and oil exploitations in modern Yellow River delta, which located at the junction of sea and land. Influenced by rivers, currents and human activities, the modern Yellow River delta has a complex composition of sediment grain size and sources. A total of 219 sediment samples from the Yellow River delta were collected in this paper. Firstly, the grain size distribution and material source of the modern Yellow River delta sediment were analyzed by particle size analysis method. Then the parametric endmember model analysis method was used to divide the material sources of sediments. Finally, the cluster analysis method was used to analyze the surface sediment types. The results show that the grain size composition of delta sediments is mainly silt (71.10%) and sand (27.62%), most of the sediments are of medium grain size sorted and the curve of the grain size distribution is nearly symmetrical. The delta sediments can be divided into three types of end-mumber such as EM1, EM2 and EM3, which represent three different material sources. Among them, EM1 is from fluvial sediments, EM2 is from tidal current sediments and EM3 is from the current Yellow River fluvial sediments. According to the results of cluster analysis, the Yellow River delta is divided into three types of sedimentary environments. The first type is mainly distributed in the middle and upper reaches of western part of the Yellow River delta rivers, including Maxin River, Zhanli River, Caoqiaogou River and Tiaohe River, which is mainly deposited by old river with weak hydrodynamic force. The second type is mainly distributed in the eastern part of rivers such as Diaokou River, Shenxiangou and the current Yellow River (Qingshuigou), and the current Yellow River has strong sedimentary dynamics. The third type is mainly distributed along the northern coast and tidal flat of the Yellow River delta, which is strongly affected by the marine tidal and the coastline is eroded. In addition, human activities also had a certain impact on the changes of the coastline in the Yellow River delta.

    Key words: Yellow River delta; surface sediments; grain size; end mumber; source analysis

    0 引言

    黃河作為我國第二大河,具有水少沙多、水沙異源、快速淤積的特點。黃河平均每年向渤海輸送8.7億t泥沙,其中70%沉積在河口附近并形成新的陸地[1]。1855年,黃河從銅瓦廂決口奪大清河匯入渤海,形成現(xiàn)代黃河三角洲[2]。近年來,受氣候變化以及人類活動的影響,黃河三角洲下游部分河流發(fā)生淤積,有的地段岸線受海洋侵蝕作用影響,發(fā)生蝕退[3]。為保護(hù)海岸生態(tài)環(huán)境,保障沿海居民生產(chǎn)生活安全,學(xué)者們對現(xiàn)代黃河三角洲潮間帶及水下三角洲地區(qū)展開了較多研究,研究內(nèi)容包括沉積物粒度特征及物質(zhì)來源[48]、沉積特征與沉積環(huán)境[913]、黃河泥沙輸送對地形地貌影響[1416]、人類活動影響[1718]、岸線變化[1921]、濕地健康評價[22]等。相比于黃河三角洲潮間帶及水下三角洲,黃河三角洲陸上部分沉積環(huán)境的研究工作較少,而且學(xué)界對于黃河三角洲地區(qū)沉積動力與物質(zhì)來源的關(guān)系認(rèn)識也不清楚,因此,有必要在黃河三角洲開展表層沉積物水動力條件及物質(zhì)來源分析。

    粒度作為沉積學(xué)研究中的重要指標(biāo),可以反映沉積物質(zhì)的搬運方式、物質(zhì)來源、沉積動力等多重信息,對研究沉積物輸運過程以及沉積環(huán)境演變具有重要意義[23]。由于河流沉積物粒度組成及分布受水動力環(huán)境和不同物質(zhì)來源的綜合影響,因此利用粒度數(shù)據(jù)分析沉積物特征及物質(zhì)來源時,常常需要對粒度組分進(jìn)行分離,其中粒度端元分析方法是區(qū)分粒度組分的有效手段。粒度端元分析通過運用數(shù)學(xué)算法將多種沉積動力作用下形成的沉積物粒度數(shù)據(jù)有效分離成多個獨立的具有某特定特征的粒度組分,從而判斷沉積物來源和識別沉積動力條件[24]。端元分析方法自提出以來,模型經(jīng)過了不斷地改進(jìn),在沉積學(xué)分析上應(yīng)用越來越廣泛[2526]。比如:張曉東等[27]對長江口臨近海域表層沉積物進(jìn)行了粒度端元分析,并對該地區(qū)沉積動力環(huán)境進(jìn)行了描述;林鎮(zhèn)坤等[28]對南流江河口水下三角洲進(jìn)行端元分析并判斷了當(dāng)?shù)氐某练e動力環(huán)境;薛成鳳等[29]運用端元分析推測了中小河流對長江水下三角洲遠(yuǎn)端泥沉積的貢獻(xiàn)。另外,在進(jìn)行沉積動力分區(qū)和探討沉積環(huán)境方面[3031],也常常運用聚類分析的方法,結(jié)合粒度特征參數(shù)(粒度參數(shù)、粒度組分)對沿海地區(qū)沉積物進(jìn)行沉積環(huán)境劃分[3233]。因此,粒度指標(biāo)可用于判斷黃河三角洲表層沉積物的物質(zhì)來源和動力環(huán)境。

    本文通過對黃河三角洲陸上沉積物的粒度特征及物質(zhì)來源進(jìn)行研究,探討不同沉積動力對黃河三角洲岸線及環(huán)境的影響,以期為黃河三角洲地區(qū)岸線防護(hù)及環(huán)境治理提供參考。

    1 區(qū)域概況與研究方法

    1.1 區(qū)域概況

    黃河三角洲位于渤海盆地東部的濟陽斷陷和埕寧隆起之間,基底構(gòu)造以強烈的褶皺和斷層為特征。黃河三角洲地勢總體西南高,東北低(圖1a)。1855年以來,黃河尾閭發(fā)生多次改流,并在不同時期形成多條分流河流。近幾十年來受人工改流控制,黃河三角洲存在多條河流并行現(xiàn)象。除現(xiàn)黃河外,河流多為南北走向,河流自西向東分別為馬新河、沾利河、草橋溝、挑河、刁口河、神仙溝、現(xiàn)黃河(圖1b),其中馬新河、沾利河、草橋溝為人工開挖的河流。黃河三角洲地區(qū)屬于暖溫帶半濕潤大陸性季風(fēng)氣候,四季差異分明,年均溫12 ℃,年降水量776 mm。海岸潮汐類型多樣,大部分為不正規(guī)半日潮[2]。植被組成以草本植物為主,木本植物所占比例較少。黃河三角洲植被的主要建群種和優(yōu)勢種為濕生植物和鹽生植物。植被形成時間較短,群落穩(wěn)定性較差[34]。

    1.2 研究方法

    研究區(qū)位于黃河下游、黃河口附近。采樣點分布于現(xiàn)代黃河三角洲的陸上部分及潮間帶區(qū)域(圖1c)。沿現(xiàn)黃河河流以及神仙溝、刁口河、挑河、草橋溝、沾利河、馬新河等黃河尾閭擺蕩形成較早的舊河流地區(qū)進(jìn)行等距采樣,采樣相鄰點間隔約為5 km,選取表層深度0~2 cm的沉積物樣品共219個。舊河流基本信息見表1。

    粒度測定方法:首先將樣品烘干,稱取約1 g樣品;然后量取10%的稀鹽酸溶液以除去樣品中的鈣質(zhì)膠結(jié),待反應(yīng)完全后加入蒸餾水并離心至中性;接著量取約5 mL雙氧水并進(jìn)行水浴加熱,以去除有機質(zhì),待反應(yīng)完全后加入蒸餾水進(jìn)行水洗并離心直至中性;再向離心管中加入0.05 mol/L的(NaPO3)6約5 mL,用超聲波清洗儀清洗樣品約10 min,使樣品充分分散;最后使用粒度分析儀Microtrac S3500進(jìn)行測試,測量范圍為0.025~2 000 μm,重復(fù)測量3次使相對誤差在±3%以內(nèi)。

    數(shù)據(jù)分析方法:首先,在Matlab中使用Paterson and Heslop提供的AnalySize程序包,采用Gen.Weibull參數(shù)模型對黃河三角洲表層沉積物粒度各粒級的體積分?jǐn)?shù)進(jìn)行分析,并對粒度數(shù)據(jù)進(jìn)行EMA(端元分析)矩陣計算。然后,用Folk和Ward公式計算各個樣品中的平均粒度、峰度、分選系數(shù)、偏度[36]。最后,運用IBM SPSS Statistics 中的K均值聚類方法,根據(jù)粒度參數(shù)Mz對樣品進(jìn)行聚類分析。

    2 結(jié)果與分析

    2.1 黃河三角洲表層沉積物粒度組成及分布

    按國際通用標(biāo)準(zhǔn)進(jìn)行分類,Φ≤4為砂,4<Φ≤8為粉砂,Φ>8為黏土(其中Φ=-log2D,D為顆粒直徑,mm)[37],得到沉積物組成按照體積分?jǐn)?shù)由高到低分別為粉砂(71.10%)、砂(27.62%)、黏土(1.28%)。從不同河流位置上表層沉積物的粒度組成中粉砂、砂和黏土的平均體積分?jǐn)?shù)來看,河流間粒度組成差異較明顯(表2)。從表2可知:整體上看,黃河三角洲表層沉積物粉砂體積分?jǐn)?shù)最高,范圍在59.93%~80.29%之間;刁口河、神仙溝、現(xiàn)黃河等東部河流地區(qū)砂體積分?jǐn)?shù)較高,分布于26.56%~38.84%之間;馬新河、沾利河、草橋溝、挑河等舊河流地區(qū)黏土體積分?jǐn)?shù)較高,在 1.82%~2.99%之間。

    黃河三角洲地區(qū)表層沉積物砂、粉砂、黏土粒度分布特征(圖2)顯示,現(xiàn)代黃河三角洲研究區(qū)范圍內(nèi)的樣品中,砂和粉砂的體積分?jǐn)?shù)在空間上的分布變化范圍較大,黏土的空間變化較小。不同河流表層沉積物粒度特征具體表現(xiàn)為:黃河三角洲表層沉積物粒度峰度平均值為1.09,近對稱;分選系數(shù)平均值為1.11,分選中等;偏度系數(shù)為0.10,近對稱。

    2.2 粒度分布曲線及概率累積曲線

    沉積物的粒度主要受搬運介質(zhì)、搬運方式、沉積環(huán)境等多種因素影響[23]。粒度概率累積曲線可以反映沉積物的搬運方式。根據(jù)累積曲線的截點以及分段可以判斷物質(zhì)的搬運是以懸移、躍移、底移中的哪種搬運方式為主,以及不同搬運方式所占的比例如何。原則上根據(jù)現(xiàn)代黃河三角洲研究區(qū)的大多數(shù)點的粒度分布曲線和概率累積曲線的形態(tài),將粒度分布曲線及概率累積曲線劃分為三組類型:類型一為馬新河、沾利河、草橋溝、挑河等舊河流;類型二為刁口河、神仙溝、現(xiàn)黃河等東部河流;類型三為沿海地區(qū)海岸及潮灘。將這三組不同類型的沉積物粒度組合取平均值,并做出相應(yīng)地區(qū)的概率累積曲線以及粒度分布曲線圖(圖3)。其中,粒度分布曲線可以反映粒度的分布,進(jìn)而可以推測沉積物的物質(zhì)來源。

    馬新河、沾利河、草橋溝、挑河等距離海岸有一定距離的舊河流平均粒度為5.26Φ(圖3a),粒級為細(xì)粉砂,分選系數(shù)為1.41,分選較差,偏度為0.09,粒度分布曲線呈近對稱分布,峰態(tài)為0.94,峰態(tài)中等。概率累積曲線為二段式(圖3a),可知物質(zhì)搬運方式主要為躍移和懸移。

    刁口河、神仙溝、現(xiàn)黃河等東部河流平均粒度為3.75Φ(圖3b),粒級為極細(xì)砂,分選系數(shù)為0.85,分選中等,偏度為0.08,粒度分布曲線呈近對稱分布,峰態(tài)為1.24,峰態(tài)窄。概率累積曲線為二段式(圖3b),物質(zhì)搬運方式主要為躍移和懸移?,F(xiàn)黃河在10Φ左右存在明顯截點,說明在沉積物的細(xì)端部分物質(zhì)的沉積方式發(fā)生了明顯變化。

    在黃河三角洲北部沿海地區(qū),平均粒度為4.62Φ(圖3c),粒級為粗粉砂,分選系數(shù)為1.03,分選較差,偏度為0.16,粒度分布曲線為正偏,峰態(tài)為1.19,峰態(tài)窄。概率累積曲線為四段式(圖3c),在6Φ和11Φ位置處存在拐點;其中躍移組分的斜率較高,反映了沉積物在能量相對較低的河口上快速堆積的特征,該處沉積物受海洋潮流的影響較為強烈。

    2.3 端元分析

    對黃河三角洲表層沉積物粒級數(shù)據(jù)進(jìn)行EMA(end member analysis)矩陣計算并得到端元數(shù)與線性相關(guān)性、角度偏差的關(guān)系(圖4),其中R2為各端元之間的判定系數(shù),角度分布曲線之間的偏離程度[24]。評估端元分解效果時,需要盡量保證 R2在0.9以上,角度偏差盡量小,同時選擇盡可能少的端元數(shù)目。隨著端元數(shù)量增加,線性相關(guān)增強,角度偏差下降。最終選擇3個端元(EM1、EM2和EM3)對該組粒度數(shù)據(jù)進(jìn)行反演并分析。

    沉積物粒度端元的性質(zhì)通常受沉積物物質(zhì)來源和沉積環(huán)境(括沉積物的運移特征以及水動力環(huán)境條件[27])的影響。3個不同端元所代表的粒度參數(shù)(表3)反映了黃河三角洲整體沉積環(huán)境受3種不同沉積動力來源的影響,結(jié)合每個樣品各端元的空間分布特征(圖5)可以分析沉積物粒度端元所指示的物質(zhì)來源及水動力環(huán)境。

    EM1代表沉積物的平均粒度和眾數(shù)粒度分別為6.13Φ和6.01Φ(表3),屬于細(xì)粉砂粒級,其中黏土體積分?jǐn)?shù)為4.07%,粉砂體積分?jǐn)?shù)為93.89%。分選系數(shù)為1.08,分選程度較中等;偏度為0.02,近對稱;峰度為0.97,峰態(tài)為中等。從空間位置分布上看,EM1物質(zhì)來源主要分布在現(xiàn)代黃河三角洲西北部沿岸潮灘以及馬新河、沾利河、草橋溝、挑河等舊河流地區(qū)。因此推測EM1為舊河流沉積動力,代表舊河流在較弱的水動力條件下逐漸沉積形成的河流沉積物。

    EM2代表沉積物的平均粒度和眾數(shù)粒度分別為4.74Φ和4.51Φ(表3),屬于粗粉砂粒級;分選系數(shù)為0.59,分選較好;偏度為0.11,偏度為正偏,顆粒在粗粒部分集中;峰度為1.00,峰態(tài)為中等。EM2主要分布在黃河三角洲北部沿海地區(qū)。北部廢棄河口岸線由于缺少河流泥沙補給,受波浪、潮流、風(fēng)暴潮等多種海洋因素作用,在刁口河口無潮點附近的海岸段岸線強烈蝕退[19]。因此推斷EM2受海洋影響強烈,為海洋潮流及波浪動力作用下形成的沉積物。

    EM3代表沉積物的平均粒度和眾數(shù)粒度分別為3.68Φ和3.76Φ(表3),屬于粗粉砂粒級;分選系數(shù)為0.56,分選較好;偏度為-0.13,為負(fù)偏,顆粒在細(xì)粒部分集中;峰度為0.96,峰態(tài)為中等。與EM1和EM2相比,EM3粒度更粗,且為負(fù)偏,在空間位置分布上看,EM3主要分布在黃河三角洲研究區(qū)的東北側(cè),以及現(xiàn)黃河。在神仙溝及刁口河流靠近陸地地區(qū)EM3含量較高,這些河流曾經(jīng)為黃河流經(jīng)河流,由于在這些地區(qū)黃河改流時間較短,且河流之間相互連通,因此仍受到黃河沉積作用影響。因此推斷EM3主要來源為現(xiàn)黃河沉積。

    從端元組成上看,EM2和EM3的占比較高,而EM1占比較低(表3)??傮w上是粉砂體積分?jǐn)?shù)高,黏土體積分?jǐn)?shù)低,沉積物粒度粗;說明黃河三角洲研究區(qū)沉積環(huán)境水動力強,沉積環(huán)境受河流和潮流等多種因素影響。近年來隨著入海泥沙減少,除了現(xiàn)黃河(清水溝)以外,其余地方受到海洋侵蝕作用加強,因而EM2占比較高(圖5)。

    3 討論

    3.1 沉積動力分區(qū)及沉積環(huán)境探討

    由于現(xiàn)代黃河三角洲粒度分布是與現(xiàn)代河口動力環(huán)境相適應(yīng)的,因此,分析沉積物的粒度分布特征可以對河口動力環(huán)境進(jìn)行推斷。

    結(jié)合現(xiàn)代黃河三角洲的粒度特征,運用K均值聚類方法對黃河三角洲表層沉積物的平均粒度進(jìn)行分類,共分為三種沉積類型,根據(jù)這三種類型反映的不同沉積特征將沉積物的空間分布劃分為三種沉積動力環(huán)境區(qū)(圖6)。對不同沉積動力環(huán)境區(qū)聚類占比及粒度特征進(jìn)行計算,結(jié)果見表4。

    沉積Ⅰ區(qū)主要分布在黃河三角洲西部河流,包括馬新河、沾利河、草橋溝、挑河等河流的上游及中游部分。從聚類組成上來看,沉積Ⅰ區(qū)中聚類1占比38.64%,聚類2占比13.64%,聚類3占比47.73%;聚類3占比最高,其平均粒度為5.46Φ(表4)。由于河流改流時間較早,河流廢棄時間較長,受人類活動用水增加以及人工堤壩的修建等影響,沉積Ⅰ區(qū)的沉積物粒度較細(xì),該地區(qū)物質(zhì)搬運方式主要以躍移和懸移為主。從岸線變化上看,結(jié)合前人[3]分析認(rèn)為該地區(qū)由于連接港口的公路修建,攔截了潮灘,岸線較為穩(wěn)定。沉積Ⅰ區(qū)的沉積環(huán)境以舊河流沉積為主。

    沉積Ⅱ區(qū)主要分布在刁口河、神仙溝、現(xiàn)黃河(清水溝)等東部河流。從聚類組成上來看,沉積Ⅱ區(qū)聚類1占比44.83%,聚類2占比48.28%,聚類3占比0.07%,主要以聚類2和聚類1為主。聚類2和聚類1平均粒度分別為4.03Φ和4.79Φ(表4)。該區(qū)域沉積環(huán)境主要為黃河河流沉積作用?,F(xiàn)行河口區(qū)主要受黃河水動力影響。一般來說,水動力環(huán)境越強,粒度越粗。由于黃河徑流量較大,流速較快,因而沉積物粒度較粗。其中,徑流區(qū)域影響下砂的含量較高。這可能是由于河流流速較快,細(xì)顆粒物質(zhì)隨水流被搬運至遠(yuǎn)處,粗顆粒物質(zhì)搬運距離較近,多在原地沉積。袁萍等[6]認(rèn)為現(xiàn)行河口三角洲葉瓣周圍的表層沉積物粒度較遠(yuǎn)離河口位置的沉積物粒度粗,且現(xiàn)黃河口沉積物粒度有粗化的趨勢。從岸線變化上看,結(jié)合前人[11]分析認(rèn)為黃河三角洲岸線在孤東以北和清水溝流路南側(cè)的區(qū)域岸線作用以侵蝕為主,而在孤東海域和清水溝流路河口處則以淤積作用為主。

    沉積Ⅲ區(qū)主要分布在黃河三角洲北部沿岸及潮灘地區(qū)。從聚類組成上來說,沉積Ⅲ區(qū)聚類1占比67.82%,聚類2占比12.64%,聚類3占比19.54%。主要以聚類1為主,其平均粒度為4.79Φ(表4)。該區(qū)域沉積環(huán)境主要為海洋沉積作用,岸線蝕退。廢棄刁口河口和神仙溝屬于波流控制型[16]。自1976年黃河改流清水溝流路以來,三角洲海岸線總的變化趨勢是現(xiàn)行河口海岸線明顯向海淤進(jìn),而北部廢棄河口岸線則明顯侵蝕后退[19]。由于潮灘地區(qū)地勢平坦,河流流速迅速減弱,水動力環(huán)境變?nèi)酰蚨淞6容^河流地區(qū)細(xì)。黃河改流后,黃河三角洲北部潮間帶部分岸灘波浪作用加強,侵蝕效應(yīng)亦增強[4]。在黃河三角洲北部沿岸沉積物主要受海洋潮流作用影響。其中在近岸,漲落潮的往復(fù)運動對細(xì)顆粒泥沙的輸移和分布有重要影響。當(dāng)來自海洋的潮流與受風(fēng)驅(qū)動下的水流、河流和海洋環(huán)流相互作用時,沉積物質(zhì)的凈運輸方向會發(fā)生相應(yīng)改變。隨著水深沿海岸變化,海底摩擦力增大,潮流速度向岸邊減小。細(xì)粒沉積物的搬運受侵蝕和沉積滯后效應(yīng)的控制,呈現(xiàn)出向岸邊搬運和沉積的趨勢。受潮流頂托作用,沉積物細(xì)顆粒物質(zhì)在海岸沉積,因而在黃河三角洲近岸粒度較細(xì)。近年來,黃河徑流減少導(dǎo)致海水潮流動力增強[15]。因而,在黃河三角洲北部海岸線有侵蝕的趨勢。

    3.2 物質(zhì)來源探討

    黃河三角洲沉積物物質(zhì)來源主要為現(xiàn)黃河河流沉積、舊河流沉積和海洋潮流沉積三個來源。黃河三角洲地區(qū)沉積物的粒度特征與物質(zhì)來源有密切關(guān)系。在黃河三角洲地區(qū),受現(xiàn)黃河河流沉積和受海洋潮流沉積影響的區(qū)域,沉積物粒度較粗,而受舊河流沉積影響的地區(qū)沉積物粒度較細(xì)。從黃河三角洲沉積物端元組成上看,現(xiàn)代黃河三角洲不同位置的端元構(gòu)成存在差異:在馬新河、沾利河、草橋溝、挑河等西部河流,端元組成以EM1為主,沉積物主要來源為舊河流沉積作用下的泥沙沉積;在刁口河、神仙溝、現(xiàn)黃河等東部河流,端元組成以EM3為主,沉積物主要物質(zhì)來源為現(xiàn)黃河河流沖刷下來的黃土高原的黃土[8];在黃河三角洲沿海地區(qū),端元組成以EM2為主,沉積物主要物質(zhì)來源可能為海洋和潮流挾沙堆積而成。這說明該區(qū)域主要受黃河河流沉積作用影響,海洋潮流僅影響黃河入??诒辈康貐^(qū)。

    本文所討論的陸上黃河三角洲沉積動力形式與前人關(guān)于黃河三角洲潮間帶及水下三角洲沉積動力形式較為一致,如陳小英等[13]認(rèn)為黃河三角洲物質(zhì)來源受到徑流、潮流以及海洋風(fēng)浪影響,喬淑卿等[7]認(rèn)為黃河三角洲沉積物物質(zhì)來源主要為河流沉積,同時受到海洋潮流動力影響;而且本研究發(fā)現(xiàn)黃河三角洲不同沉積區(qū)沉積動力形式具有明顯差別,其中本研究沉積Ⅲ區(qū)與陳小英等[13]劃分的廢棄三角洲沉積區(qū)距離較近,兩項研究均揭示了該地區(qū)主要受海洋沉積作用的影響。此外,本研究中黃河與陳小英等[13]

    研究的現(xiàn)行河口區(qū)位置一致,主要受河流作用以及沙嘴前方強潮流帶的作用??傊S河三角洲沉積動力同時受河流和海洋潮流影響,在黃河三角洲不同沉積區(qū)域沉積物主導(dǎo)的動力因素有所不同。

    黃河三角洲沉積物的粒度變化主要受到黃河影響,為河控型三角洲。1855年以來,黃河頻繁改流導(dǎo)致黃河三角洲不同河流上沉積物特征均不相同。黃學(xué)勇等[12]認(rèn)為黃河徑流是現(xiàn)代黃河口南岸沉積物的主要物質(zhì)來源,黃河改流對黃河三角洲沉積環(huán)境產(chǎn)生重要影響。對比黃河三角洲不同位置的黃河河流沉積物粒度特征并結(jié)合黃河改流歷史發(fā)現(xiàn),較早時期廢棄的河流粒度較細(xì),而現(xiàn)黃河以及廢棄時間較晚的河流由于仍主要受現(xiàn)河流流水搬運作用影響,河流流速較快,水動力環(huán)境較強,因而粒度更粗。在沾利河、草橋溝和挑河入海口處有少量代表現(xiàn)黃河的沉積物。這說明黃河改流后,仍有少量黃河沉積物通過西部舊河流運輸至入海口。

    近年來,人類工程活動、修建堤防、人工改流和開挖溝槽(尤其是河口疏浚)對黃河三角洲和其河口產(chǎn)生了巨大的影響。隨著河口向海延伸,人們在黃河三角洲陸地上種植樹木,用來防止海浪沖刷侵蝕,在油田附近建造防浪墻來防止海岸侵蝕[17]。在黃河三角洲北部,受修建堤壩、鞏固岸灘、石油開采和人工養(yǎng)殖業(yè)發(fā)展等人類活動的干擾,人工岸線向海推進(jìn)。許多島嶼并入陸地,海灣被攔截,原來曲折的天然韻律型海岸被平直的人工海岸取代[21]。因此,對黃河三角洲地區(qū)進(jìn)行開發(fā)時應(yīng)注意生態(tài)防護(hù)與經(jīng)濟建設(shè)協(xié)調(diào)發(fā)展,以減少黃河三角洲北部岸線侵蝕對生態(tài)環(huán)境的不利影響。

    4 結(jié)論與建議

    1)現(xiàn)代黃河三角洲沉積物組成以粉砂(體積分?jǐn)?shù)71.10%)、砂(27.62%)為主,黏土(1.28%)較少。在馬新河、沾利河、草橋溝、挑河等大部分舊河流地區(qū),黏土體積分?jǐn)?shù)相對其他河道較高;在刁口河、神仙溝、現(xiàn)黃河以及部分沿海地區(qū),砂體積分?jǐn)?shù)較高。在馬新河、沾利河、草橋溝、挑河等距離海岸有一定距離的舊河流,粒級為細(xì)粉砂,分選較差,粒度分布曲線呈近對稱分布,峰態(tài)中等,概率累積曲線為兩段式,物質(zhì)搬運方式主要為躍移和懸移;在刁口河、神仙溝、現(xiàn)黃河等東部河流,粒級為極細(xì)砂,分選中等,粒度分布曲線呈近對稱分布,峰態(tài)窄,概率累積曲線為二段式,物質(zhì)搬運方式主要為躍移和懸移;在黃河三角洲北部沿海地帶,粒級為粗粉砂,分選較差,粒度分布曲線為正偏,峰態(tài)窄,概率累積曲線為四段式,其中躍移組分的斜率較高。

    2)黃河三角洲表層沉積物由三個端元構(gòu)成,其中EM1為舊河流沉積動力,代表舊河流在較弱的水動力條件下逐漸沉積形成的河流沉積物,EM2為海洋潮流及波浪動力作用下形成的沉積物,EM3主要來源為現(xiàn)黃河沉積物。

    3)運用聚類分析方法將黃河三角洲分為三類沉積環(huán)境:沉積Ⅰ區(qū)主要分布在黃河三角洲等西部河流的中上游,包括馬新河、沾利河、草橋溝、挑河,主要為舊河流沉積,水動力較弱;沉積Ⅱ區(qū)主要分布在刁口河、神仙溝、現(xiàn)黃河(清水溝)等東部河流,現(xiàn)黃河河流搬運作用較強;沉積Ⅲ區(qū)主要分布在黃河三角洲北部沿岸及潮灘,受海洋潮流作用較強,海岸受到侵蝕。

    人類活動對黃河三角洲和其河口產(chǎn)生了一定的影響,在此背景下,黃河三角洲岸線發(fā)生動態(tài)改變。在黃河三角洲神仙溝以東部分岸線保存良好,但在刁口河以西岸線仍受到海洋侵蝕,應(yīng)加強岸線防護(hù)。

    參考文獻(xiàn)(References):

    [1] Zhang J, Huang W W, Shi M C. Huanghe (Yellow River) and Its Estuary: Sediment Origin, Transport and Deposition[J]. Journal of Hydrology, 1990, 120(1/2/3/4): 203223.

    [2] 成國棟. 黃河三角洲現(xiàn)代沉積作用及模式[M]. 北京:地質(zhì)出版社, 1991.

    Cheng Guodong. Modern Sedimentation and Model of Yellow River Delta [M]. Beijing: Geological Publishing House, 1991.

    [3] 尹延鴻, 周永青, 丁東. 現(xiàn)代黃河三角洲海岸演化研究[J]. 海洋通報, 2004, 23(2): 3240.

    Yin Yanhong, Zhou Yongqing, Ding Dong. Evolultion of Modern Yellow River Delt Coast[J]. Marine Science Bulletin, 2004, 23(2): 3240.

    [4] 薛允傳, 尹延鴻, 高抒. 黃河三角洲北部潮間帶沉積物的粒度特征[J]. 海洋科學(xué), 2001, 25(5): 5054.

    Xue Yunchuan, Yin Yanhong, Gao Shu. Grain Size Characteristics of Tidal Flat Sediments, Northern Yellow River Delta[J]. Marine Sciences, 2001, 25(5): 5054.

    [5] 李九發(fā), 李為華, 應(yīng)銘, 等. 黃河三角洲飛雁灘沉積物顆粒度分布和粒度參數(shù)特征及水動力解釋[J]. 海洋通報, 2006,25(3): 3844.

    Li Jiufa, Li Weihua, Ying Ming, et al. Characteristic and Hydrodynamic Explanation of Distribution and Parameters of Sediment Granularity in Feiyan Shoal of Yellow River Delta[J]. Marine Science Bulletin, 2006, 25(3): 3844.

    [6] 袁萍, 畢乃雙, 吳曉, 等. 現(xiàn)代黃河三角洲表層沉積物的空間分布特征[J]. 海洋地質(zhì)與第四紀(jì)地質(zhì), 2016, 36(2): 4957.

    Yuan Ping, Bi Naishuang, Wu Xiao, et al. Surface Sediments at the Subaqueous Yellow River Delta: Classification and Distribution[J]. Marine Geology amp; Quaternary Geology, 2016, 36(2): 4957.

    [7] 喬淑卿, 方習(xí)生, 石學(xué)法, 等. 黃河口及鄰近渤海海域表層沉積物中CaO和蒙皂石分布及其對黃河入海物質(zhì)運移的指示[J]. 海洋地質(zhì)與第四紀(jì)地質(zhì), 2010, 30(1): 1723.

    Qiao Shuqing, Fang Xisheng, Shi Xuefa, et al. Distribution of CaO and Smectite in Surface Sediments off the Yellow River Mouth and in the Nearby Bohai Sea and the Implications for Dispersion of the River Sediments to the Sea[J]. Marine Geology amp; Quaternary Geology, 2010, 30(1): 1723.

    [8] 王昆山, 石學(xué)法, 蔡善武, 等. 黃河口及萊州灣表層沉積物中重礦物分布與來源[J]. 海洋地質(zhì)與第四紀(jì)地質(zhì), 2010, 30(6): 18.

    Wang Kunshan, Shi Xuefa, Cai Shanwu, et al. Distribution and Provenance of the Surface Sediments of the Yellow River Mouth and Laizhou Bay Deduced from Heavy Minerals[J]. Marine Geology amp; Quaternary Geology, 2010, 30(6): 18.

    [9] 李栓科. 近代黃河三角洲的沉積特征[J]. 地理研究, 1989, 8(4): 4555.

    Li Shuanke. Sedimentary Characteristics of the Modern Yellow River Delta[J]. Geographical Research, 1989,8(4): 4555.

    [10] 任韌希子. 黃河三角洲沉積特征與環(huán)境演變研究[D]. 上海: 華東師范大學(xué), 2012.

    Ren Renxizi. Sedimentary Characteristics and Environment Evolution of the Yellow River Delta[D]. Shanghai: East China Normal University, 2012.

    [11] 趙廣明, 葉青, 薛春汀, 等. 現(xiàn)代黃河三角洲陸上表層沉積物類型與沉積環(huán)境分區(qū)及岸線演變[J]. 海洋地質(zhì)與第四紀(jì)地質(zhì), 2013, 33(5): 4752.

    Zhao Guangming, Ye Qing, Xue Chunting, et al. Surface Sediments, Sedimentary Subenvironments and Shoreline Evolution of Modern Yellow River Delta[J]. Marine Geology amp; Quaternary Geology, 2013, 33(5): 4752.

    [12] 黃學(xué)勇, 張戈, 高茂生, 等. 現(xiàn)代黃河口南岸海洋沉積環(huán)境分析[J]. 海洋地質(zhì)前沿, 2019, 35(1): 1221.

    Huang Xueyong, Zhang Ge, Gao Maosheng, et al. Analysis of Marine Sedimentary Environment on the South Bank of the Modern Yellow River Estuary[J]. Marine Geology Frontier, 2019, 35(1): 1221.

    [13] 陳小英, 陳沈良, 劉勇勝. 黃河三角洲濱海區(qū)沉積物的分異特征與規(guī)律[J]. 沉積學(xué)報, 2006, 24(5): 714721.

    Chen Xiaoying, Chen Shenliang, Liu Yongsheng. Sediment Differentiation Along Nearshore Zone of the Yellow River Delta[J]. Acta Sedimentologica Sinica, 2006, 24(5): 714721.

    [14] 王一鳴. 黃河水沙變化及其對黃河三角洲變化的影響分析[D]. 北京:中國科學(xué)院大學(xué), 2018.

    Wang Yiming. Analysis on the Water and Sediment Variation in Yellow River and Its Influence on the Change of Yellow River Delta[D]. Beijing:University of Chinese Academy of Sciences, 2018.

    [15] 樊輝. 黃河口泥沙輸移及三角洲的近期演變[D].青島:中國科學(xué)院研究生院, 2005.

    Fan Hui. Sediment Transport in the Huanghe River Mouth and the Recent Deltaic Evolution [D]. Qingdao: Graduate School of the Chinese Academy of Sciences, 2005.

    [16] 劉鋒. 黃河口及其鄰近海域泥沙輸運及其動力地貌過程[D]. 上海:華東師范大學(xué), 2012.

    Liu Feng. Sediment Transport and Dynamic Geomorphology Process in the Yellow River Estuary and Its Adjacent Sea[D]. Shanghai: East China Normal University, 2012.

    [17] 王兆印, 程東升, 劉成. 人類活動對典型三角洲演變的影響:Ⅱ:黃河和海河三角洲[J]. 泥沙研究,2006,25(1): 7680.

    Wang Zhaoyin, Cheng Dongsheng, Liu Cheng. Delta Processes and Management Strategies in China:Ⅱ :The Yellow and Haihe River Deltas[J]. Journal of Sediment Research, 2006,25(1): 7680.

    [18] 繆紅兵, 喬璐璐, 仲毅, 等. 人類活動和自然演變共同驅(qū)動下黃河三角洲海域潮波及物質(zhì)輸運變化[J]. 海洋學(xué)報, 2022, 44(9): 7386.

    Miao Hongbing, Qiao Lulu, Zhong Yi, et al. Evolution of Tidal System and Material Transport off the Huanghe River Delta Induced by Human Activities and Natural Evolution[J]. Haiyang Xuebao, 2022, 44(9): 7386.

    [19] 楊偉. 現(xiàn)代黃河三角洲海岸線變遷及灘涂演化[J].海洋地質(zhì)前沿, 2012, 28(7): 1723.

    Yang Wei. Shifting of Coastline and Evolution of Tidal Flat in Modern Yellow River Delta[J]. Marine Geology Frontiers, 2012, 28(7): 1723.

    [20] 黃波, 郭慧敏, 孫吉剛, 等. 黃河三角洲刁口河海岸侵蝕時空演變與防護(hù)對策[M]. 鄭州: 黃河水利出版社, 2020.

    Huang Bo,Guo Huimin, Sun Jigang, et al. Spatial-Temporal Evolution of Coastal Erosion and Protection Countermeasures of the Diaokou River in the Yellow River Delta[M]. Zhengzhou: Yellow River Water Conservancy Press, 2020.

    [21] 黃海軍, 樊輝. 1976年黃河改流以來三角洲近岸區(qū)變化遙感監(jiān)測[J]. 海洋與湖沼, 2004, 35(4): 306314.

    Huang Haijun, Fan Hui. Monitoring Changes of Nearshore Zones in the Huanghe (Yellow River) Delta Since 1976[J]. Oceanologia et Limnologia Sinica, 2004, 35(4): 306314.

    [22] 安樂生, 劉貫群, 葉思源, 等. 黃河三角洲濱海濕地健康條件評價[J]. 吉林大學(xué)學(xué)報(地球科學(xué)版), 2011, 41(4): 11571165.

    An Lesheng, Liu Guanqun, Ye Siyuan, et al. Health Assessment of Coastal Wetlands in the Yellow River Delta[J]. Journal of Jilin University(Earth Science Edition), 2011, 41(4): 11571165.

    [23] 成都地質(zhì)學(xué)院陜北隊. 沉積巖(物)粒度分析及其應(yīng)用[M]. 北京: 地質(zhì)出版社,1978.

    Shanbei Team,

    Chengdu College of Geology. Sedimentary Rock (Material) Particle-Size Analysis and Its Application[M]. Beijing: Geological Publishing House, 1978.

    [24] Weltje G J. End-Member Modeling of Compositional Data: Numerical-Statistical Algorithms for Solving the Explicit Mixing Problem[J]. Mathematical Geology, 1997, 29(4): 503549.

    [25] 白敏, 魯瑞潔, 丁之勇, 等.青海湖湖東沙地粒度端元分析及其指示意義[J]. 第四紀(jì)研究, 2020, 40(5): 12031215.

    Bai Min, Lu Ruijie, Ding Zhiyong, et al. End-Member Analysis of GrainSize in the East of Qinghai Lake and Its Environmental Implications[J]. Quaternary Sciences, 2020, 40(5): 12031215.

    [26] Paterson G A, Heslop D. New Methods for Unmixing Sediment Grain Size Data[J]. Technical Reports : Methods, 2015, 16(12): 44944506.

    [27] 張曉東, 翟世奎, 許淑梅. 端元分析模型在長江口鄰近海域沉積物粒度數(shù)據(jù)反演方面的應(yīng)用[J]. 海洋學(xué)報, 2006, 28(4): 159166.

    Zhang Xiaodong, Zhai Shikui, Xu Shumei. The Application of Grain-Size End-Member Modeling to the Shelf Near the Estuary of Changjiang River in China[J]. Haiyang Xuebao, 2006, 28(4): 159166.

    [28] 林鎮(zhèn)坤, 王愛軍, 葉翔. 南流江河口水下三角洲表層沉積物端元分析及其沉積動力環(huán)境意義[J]. 沉積學(xué)報, 2019, 37(1): 124134.

    Lin Zhenkun, Wang Aijun, Ye Xiang. End-Member Analysis for Surficial Sediment of Nanliujiang River Subaqueous Delta and Associated Sediment Dynamic Environmental Significance[J]. Acta Geologica Sinica, 2019, 37(1): 124134.

    [29] 薛成鳳, 賈建軍, 高抒, 等. 中小河流對長江水下三角洲遠(yuǎn)端泥沉積的貢獻(xiàn):以椒江和甌江為例[J]. 海洋學(xué)報, 2018, 40(5): 7589.

    Xue Chengfeng, Jia Jianjun, Gao Shu, et al. The Contribution of Middle and Small Rivers to the Distal Mud of Subaqueous Changjiang Delta: Results from Jiaojiang River and Oujiang River[J]. Haiyang Xuebao, 2018, 40(5): 7589.

    [30] 杜新強,王鈺升,冶雪艷,等. 大清河流域平原區(qū)地下水人工補給潛力與補給方式分析[J].吉林大學(xué)學(xué)報(地球科學(xué)版), 2022, 52(2):535549.

    Du Xinqiang, Wang Yusheng, Ye Xueyan,et al. Artificial Aquifer Recharge Potential and Methods in Plain Area of Daqing River Basin[J]. Journal of Jilin University (Earth Science Edition), 2022, 52(2):535549.

    [31] 洪梅,趙明明,魏濤,等.基于SWAT模型的翔陽河流域斷流河道復(fù)流研究[J].吉林大學(xué)學(xué)報(地球科學(xué)版),2024,54(1);253263.

    Hong Mei,Zhao Mingming,Wei Tao,et al. Study on Water Restoration in the Dried-up River of the Fuyang River Basin Based on the SWAT Model[J]. Journal of Jilin University (Earth Science Edition),2024,54( 1):253263.

    [32] 張盼, 吳建政, 胡日軍, 等. 萊州灣西南部表層沉積物粒度分布特征及其現(xiàn)代沉積環(huán)境分區(qū)[J]. 海洋地質(zhì)前沿, 2014, 30(9): 1117.

    Zhang Pan, Wu Jianzheng, Hu Rijun, et al. Distribution of Surface Sediments and Modern Sedimentary Environment in the Southwestern Laizhou Bay[J]. Marine Geology Frontiers, 2014, 30(9): 1117.

    [33] 楊陽, 高抒, 周亮, 等. 海南新村港潟湖表層沉積物粒度特征及其沉積環(huán)境[J]. 海洋學(xué)報, 2016, 38(1): 94105.

    Yang Yang, Gao Shu, Zhou Liang, et al. Grain Size Distribution of Surface Sediments and Sedimentary Environment in the Lagoon of Xincun, Hainan Island[J]. Haiyang Xuebao, 2016, 38(1): 94105.

    [34] 邢尚軍, 郗金標(biāo), 張建鋒, 等. 黃河三角洲植被基本特征及其主要類型[J]. 東北林業(yè)大學(xué)學(xué)報, 2003, 31(6): 8586.

    Xing Shangjun, Xi Jinbiao, Zhang Jianfeng, et al. Basic Characteristics and Main Types of Vegetation in the Yellow River Delta[J]. Journal of Northeast Forestry University, 2003, 31(6): 8586.

    [35] 東營市人民政府. 自然環(huán)境 [EB/OL]. (20230514).http://www.dongying.gov.cn/col/col307764/index.html.

    Dongying Municipal People’s Government. Natural Environment [EB/OL]. (20230514)." http://www.dongying.gov.cn/col/col307764/index.html.

    [36] Folk R L, Ward W C. Brazos River Bar (Texas): A Study in the Significance of Grain Size Parameters[J]. Journal of Sedimentary Research, 1957, 27(1): 326.

    [37] Wentworth C K. A Scale of Grade and Class Terms for Clastic Sediments[J]. Journal of Geology, 1922, 30(5): 377392.

    猜你喜歡
    端元黃河三角洲粒度
    基于優(yōu)化K-P-Means解混方法的高光譜圖像礦物識別
    粉末粒度對純Re坯顯微組織與力學(xué)性能的影響
    基于矩陣的多粒度粗糙集粒度約簡方法
    南昌地區(qū)不透水面遙感估算研究
    黃河三角洲保護(hù)區(qū)自然資源的開發(fā)與保護(hù)
    兩種基于異常權(quán)重的N-FINDR端元提取算法
    基于粒度矩陣的程度多粒度粗糙集粒度約簡
    基于Gram行列式的快速端元提取方法
    黃河三角洲濱海濕地維管植物多樣性現(xiàn)狀及保護(hù)策略
    城鎮(zhèn)化與生態(tài)環(huán)境交互協(xié)調(diào)行為研究——以黃河三角洲為例
    丁香六月天网| 亚洲av电影在线进入| 看十八女毛片水多多多| 在线 av 中文字幕| 18禁动态无遮挡网站| 永久免费av网站大全| 老女人水多毛片| 午夜福利影视在线免费观看| 青青草视频在线视频观看| 国产激情久久老熟女| 国产国语露脸激情在线看| 色吧在线观看| 国产男女内射视频| 两个人看的免费小视频| 日日爽夜夜爽网站| 老司机亚洲免费影院| 亚洲av综合色区一区| videos熟女内射| 秋霞伦理黄片| 大香蕉久久网| 国产精品久久久久久久电影| 久久国产精品大桥未久av| 赤兔流量卡办理| 51国产日韩欧美| 三级国产精品片| 亚洲国产成人一精品久久久| 亚洲第一av免费看| 亚洲久久久国产精品| 国产白丝娇喘喷水9色精品| 交换朋友夫妻互换小说| 哪个播放器可以免费观看大片| 十分钟在线观看高清视频www| 国产片内射在线| 啦啦啦中文免费视频观看日本| 免费观看无遮挡的男女| 中文乱码字字幕精品一区二区三区| 成人二区视频| 日本欧美视频一区| 亚洲,欧美,日韩| av免费观看日本| 波多野结衣一区麻豆| 日韩三级伦理在线观看| 丰满迷人的少妇在线观看| 亚洲婷婷狠狠爱综合网| 在线免费观看不下载黄p国产| 午夜精品国产一区二区电影| 男女无遮挡免费网站观看| 久久亚洲国产成人精品v| 亚洲性久久影院| 免费观看在线日韩| 日韩人妻精品一区2区三区| 波多野结衣一区麻豆| 七月丁香在线播放| 欧美性感艳星| 婷婷色综合www| 秋霞在线观看毛片| 咕卡用的链子| 三级国产精品片| 亚洲国产精品成人久久小说| 一级黄片播放器| 亚洲少妇的诱惑av| 欧美老熟妇乱子伦牲交| 免费观看在线日韩| 看非洲黑人一级黄片| 美女主播在线视频| 欧美另类一区| 美国免费a级毛片| 亚洲高清免费不卡视频| 久久人人爽人人片av| 亚洲av在线观看美女高潮| 亚洲人成网站在线观看播放| www.av在线官网国产| 亚洲久久久国产精品| 国产福利在线免费观看视频| 欧美日韩精品成人综合77777| 自拍欧美九色日韩亚洲蝌蚪91| 天美传媒精品一区二区| av女优亚洲男人天堂| 久久人人97超碰香蕉20202| 国产色爽女视频免费观看| 国产精品不卡视频一区二区| 成人毛片60女人毛片免费| av免费观看日本| 一级毛片电影观看| 亚洲国产最新在线播放| 美女视频免费永久观看网站| 免费人成在线观看视频色| 少妇高潮的动态图| 蜜臀久久99精品久久宅男| 亚洲一级一片aⅴ在线观看| 久久午夜综合久久蜜桃| 男女国产视频网站| 日本猛色少妇xxxxx猛交久久| 久久人人爽av亚洲精品天堂| 欧美人与性动交α欧美软件 | 国产精品99久久99久久久不卡 | 亚洲国产成人一精品久久久| 伦理电影大哥的女人| 26uuu在线亚洲综合色| 免费观看无遮挡的男女| 精品亚洲乱码少妇综合久久| 国产熟女欧美一区二区| 少妇的丰满在线观看| 又黄又粗又硬又大视频| 亚洲欧美色中文字幕在线| 久久人妻熟女aⅴ| 国产熟女午夜一区二区三区| 国产精品三级大全| 一级毛片我不卡| 欧美日韩亚洲高清精品| 午夜视频国产福利| 狠狠婷婷综合久久久久久88av| av免费观看日本| 久久久精品94久久精品| 99久久精品国产国产毛片| 久久久久网色| 热re99久久国产66热| 久久韩国三级中文字幕| 久久精品久久久久久噜噜老黄| 亚洲经典国产精华液单| 国产在线一区二区三区精| 欧美3d第一页| 欧美成人午夜免费资源| 亚洲成人一二三区av| 亚洲国产色片| 久久久久久人妻| 宅男免费午夜| 天天操日日干夜夜撸| 日本vs欧美在线观看视频| 黑人欧美特级aaaaaa片| 最近的中文字幕免费完整| 在线观看人妻少妇| 性高湖久久久久久久久免费观看| 欧美精品一区二区免费开放| 欧美国产精品一级二级三级| 欧美 亚洲 国产 日韩一| 日韩制服丝袜自拍偷拍| 欧美老熟妇乱子伦牲交| 成年美女黄网站色视频大全免费| 亚洲国产精品一区二区三区在线| 大香蕉久久网| 国产不卡av网站在线观看| tube8黄色片| 天天躁夜夜躁狠狠躁躁| 精品99又大又爽又粗少妇毛片| 久久精品国产亚洲av天美| av黄色大香蕉| 亚洲国产看品久久| 国产精品人妻久久久久久| 蜜桃在线观看..| 不卡视频在线观看欧美| 精品国产一区二区三区久久久樱花| 天天操日日干夜夜撸| 免费大片黄手机在线观看| 少妇精品久久久久久久| av视频免费观看在线观看| 99久久人妻综合| 午夜福利视频在线观看免费| 亚洲国产看品久久| 国产精品.久久久| 日本免费在线观看一区| 色网站视频免费| 久久久久久久久久成人| 色94色欧美一区二区| 国产精品 国内视频| av播播在线观看一区| 欧美bdsm另类| 精品久久国产蜜桃| 一区二区三区四区激情视频| 91久久精品国产一区二区三区| 久久99热这里只频精品6学生| 99九九在线精品视频| 美女内射精品一级片tv| 99久国产av精品国产电影| 亚洲精品国产av成人精品| 国精品久久久久久国模美| 黄色毛片三级朝国网站| 久久午夜福利片| 99久久综合免费| 在线亚洲精品国产二区图片欧美| 午夜视频国产福利| 黑人巨大精品欧美一区二区蜜桃 | 老司机影院成人| 99久久中文字幕三级久久日本| 国产精品免费大片| 久久影院123| 一级毛片我不卡| 香蕉精品网在线| a级片在线免费高清观看视频| 国产福利在线免费观看视频| av免费在线看不卡| 国产1区2区3区精品| 妹子高潮喷水视频| 赤兔流量卡办理| 成人毛片a级毛片在线播放| 国产成人免费观看mmmm| 亚洲欧美日韩卡通动漫| 国产无遮挡羞羞视频在线观看| 午夜免费观看性视频| 午夜福利在线观看免费完整高清在| 日韩制服骚丝袜av| 日本爱情动作片www.在线观看| 亚洲五月色婷婷综合| 久久久久久久亚洲中文字幕| 天堂中文最新版在线下载| 我的女老师完整版在线观看| a级毛片黄视频| 免费看光身美女| 各种免费的搞黄视频| 少妇人妻精品综合一区二区| 成人免费观看视频高清| 国产亚洲午夜精品一区二区久久| 在线亚洲精品国产二区图片欧美| 中文天堂在线官网| 亚洲欧美日韩卡通动漫| 在线观看一区二区三区激情| 亚洲精华国产精华液的使用体验| 人人澡人人妻人| 综合色丁香网| 欧美成人精品欧美一级黄| 久久久久精品人妻al黑| 哪个播放器可以免费观看大片| 自拍欧美九色日韩亚洲蝌蚪91| 婷婷色综合www| 国产爽快片一区二区三区| 国产成人午夜福利电影在线观看| 一级毛片黄色毛片免费观看视频| 大香蕉久久成人网| 久久人人爽人人爽人人片va| 亚洲四区av| 在线观看人妻少妇| 国产精品久久久久久久电影| 亚洲av男天堂| 99九九在线精品视频| 精品人妻一区二区三区麻豆| 成人亚洲欧美一区二区av| 国产国语露脸激情在线看| 国产1区2区3区精品| 欧美成人午夜免费资源| 日本黄色日本黄色录像| 午夜av观看不卡| 日本av手机在线免费观看| 国产精品久久久久久精品古装| 亚洲精品美女久久久久99蜜臀 | 久久久久精品久久久久真实原创| 日本av手机在线免费观看| 内地一区二区视频在线| 波野结衣二区三区在线| 欧美精品一区二区免费开放| 热99国产精品久久久久久7| 亚洲精品美女久久av网站| 国产精品无大码| 黄色毛片三级朝国网站| 另类亚洲欧美激情| 欧美日韩精品成人综合77777| 中文字幕精品免费在线观看视频 | 街头女战士在线观看网站| 亚洲内射少妇av| 国产片内射在线| 夜夜骑夜夜射夜夜干| 国产精品久久久久久久久免| 两个人免费观看高清视频| 女性生殖器流出的白浆| 91aial.com中文字幕在线观看| 久久久久久久精品精品| 黄色 视频免费看| 精品人妻熟女毛片av久久网站| 国产一区二区在线观看日韩| 日日爽夜夜爽网站| 国产深夜福利视频在线观看| 人体艺术视频欧美日本| 女人精品久久久久毛片| 免费日韩欧美在线观看| 国产又爽黄色视频| 精品人妻偷拍中文字幕| www.熟女人妻精品国产 | 91午夜精品亚洲一区二区三区| 久久影院123| 中文字幕亚洲精品专区| 国产av码专区亚洲av| 少妇人妻久久综合中文| 久久久久久久久久久免费av| 黄色一级大片看看| 亚洲av中文av极速乱| 多毛熟女@视频| 99热6这里只有精品| 国产成人一区二区在线| 美女大奶头黄色视频| 国产精品一区二区在线不卡| 丰满饥渴人妻一区二区三| 黄色一级大片看看| 美女国产视频在线观看| 最近手机中文字幕大全| 久久精品熟女亚洲av麻豆精品| 精品人妻偷拍中文字幕| 国产成人av激情在线播放| 熟妇人妻不卡中文字幕| av.在线天堂| 夫妻午夜视频| 超碰97精品在线观看| 色网站视频免费| 激情五月婷婷亚洲| 晚上一个人看的免费电影| 免费高清在线观看日韩| 中国三级夫妇交换| 女的被弄到高潮叫床怎么办| 欧美日韩成人在线一区二区| 亚洲国产精品成人久久小说| 久久人妻熟女aⅴ| 丝袜美足系列| 国产一区二区三区综合在线观看 | 一级爰片在线观看| 久久这里只有精品19| 另类精品久久| 亚洲一区二区三区欧美精品| 曰老女人黄片| 日本黄大片高清| 精品卡一卡二卡四卡免费| 91精品三级在线观看| www.熟女人妻精品国产 | 久久久久人妻精品一区果冻| 久热这里只有精品99| 十分钟在线观看高清视频www| 中文字幕制服av| 久久99精品国语久久久| 51国产日韩欧美| 肉色欧美久久久久久久蜜桃| 国产精品久久久久久av不卡| 国产又色又爽无遮挡免| 亚洲av男天堂| 99热全是精品| 欧美日韩一区二区视频在线观看视频在线| 王馨瑶露胸无遮挡在线观看| 80岁老熟妇乱子伦牲交| 99热国产这里只有精品6| 秋霞在线观看毛片| 一级,二级,三级黄色视频| 精品一区二区三区视频在线| 成人免费观看视频高清| 久久 成人 亚洲| 欧美+日韩+精品| 男女高潮啪啪啪动态图| 亚洲丝袜综合中文字幕| a级毛片黄视频| 免费av不卡在线播放| 久久99热这里只频精品6学生| 久久亚洲国产成人精品v| 多毛熟女@视频| 美女中出高潮动态图| 十分钟在线观看高清视频www| 亚洲图色成人| 搡女人真爽免费视频火全软件| 女的被弄到高潮叫床怎么办| 欧美日韩精品成人综合77777| 另类精品久久| av在线观看视频网站免费| 少妇 在线观看| 日韩人妻精品一区2区三区| a级毛片黄视频| 国产成人精品久久久久久| 性色avwww在线观看| 亚洲经典国产精华液单| 美女内射精品一级片tv| 亚洲人成77777在线视频| 亚洲欧美日韩卡通动漫| 久久鲁丝午夜福利片| 久久热在线av| 国产午夜精品一二区理论片| 黑人猛操日本美女一级片| 午夜免费观看性视频| 免费久久久久久久精品成人欧美视频 | 国产黄色视频一区二区在线观看| 精品人妻偷拍中文字幕| 少妇人妻精品综合一区二区| 最新的欧美精品一区二区| 欧美xxxx性猛交bbbb| 国产女主播在线喷水免费视频网站| 免费日韩欧美在线观看| 最近最新中文字幕大全免费视频 | 亚洲av国产av综合av卡| 国产精品熟女久久久久浪| 欧美+日韩+精品| 又粗又硬又长又爽又黄的视频| 日韩成人av中文字幕在线观看| 国产国语露脸激情在线看| 91久久精品国产一区二区三区| 国产日韩欧美视频二区| 成年人免费黄色播放视频| 国产 一区精品| 午夜老司机福利剧场| 欧美激情 高清一区二区三区| 99热网站在线观看| 国产69精品久久久久777片| 亚洲第一区二区三区不卡| 久久久久久久精品精品| 精品视频人人做人人爽| 日韩精品有码人妻一区| 午夜激情av网站| 亚洲av电影在线进入| 国产日韩一区二区三区精品不卡| 新久久久久国产一级毛片| 午夜日本视频在线| 2021少妇久久久久久久久久久| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 欧美老熟妇乱子伦牲交| 亚洲精品色激情综合| 女的被弄到高潮叫床怎么办| 永久免费av网站大全| 国语对白做爰xxxⅹ性视频网站| 亚洲成色77777| 精品久久久精品久久久| 久久99一区二区三区| a 毛片基地| 一级片免费观看大全| 亚洲精品第二区| 国产黄频视频在线观看| 午夜精品国产一区二区电影| 久久久久久伊人网av| av在线老鸭窝| 夫妻午夜视频| 亚洲成av片中文字幕在线观看 | 久久99热这里只频精品6学生| 丝袜喷水一区| 成人二区视频| 成人免费观看视频高清| 久久精品久久精品一区二区三区| 少妇被粗大的猛进出69影院 | 五月玫瑰六月丁香| 成年av动漫网址| 日本与韩国留学比较| 女性被躁到高潮视频| 亚洲国产日韩一区二区| 狂野欧美激情性xxxx在线观看| 免费少妇av软件| 一区二区av电影网| 国产一级毛片在线| freevideosex欧美| 午夜福利影视在线免费观看| 99久久人妻综合| 狠狠婷婷综合久久久久久88av| 国国产精品蜜臀av免费| 韩国精品一区二区三区 | 国产精品熟女久久久久浪| 成人亚洲精品一区在线观看| 制服人妻中文乱码| 国产淫语在线视频| 成人国产av品久久久| 亚洲成av片中文字幕在线观看 | av福利片在线| 日韩一区二区三区影片| 欧美精品人与动牲交sv欧美| 咕卡用的链子| 国产成人午夜福利电影在线观看| 欧美精品国产亚洲| tube8黄色片| 日本wwww免费看| 飞空精品影院首页| 国产精品 国内视频| 一边摸一边做爽爽视频免费| 精品视频人人做人人爽| 狂野欧美激情性xxxx在线观看| 久久久久久久久久成人| 一二三四中文在线观看免费高清| 美女脱内裤让男人舔精品视频| 在现免费观看毛片| 亚洲av中文av极速乱| 飞空精品影院首页| 18禁裸乳无遮挡动漫免费视频| 午夜福利乱码中文字幕| 天美传媒精品一区二区| 美女脱内裤让男人舔精品视频| av女优亚洲男人天堂| 国产黄频视频在线观看| 精品亚洲乱码少妇综合久久| 一级a做视频免费观看| 深夜精品福利| 丰满少妇做爰视频| 一级,二级,三级黄色视频| 国产一级毛片在线| 十分钟在线观看高清视频www| 在线观看免费视频网站a站| 黑人欧美特级aaaaaa片| 黑丝袜美女国产一区| 国产精品嫩草影院av在线观看| 丝袜人妻中文字幕| 日韩制服丝袜自拍偷拍| 亚洲内射少妇av| 亚洲成av片中文字幕在线观看 | 下体分泌物呈黄色| 一区在线观看完整版| 国产亚洲精品第一综合不卡 | 久久久久精品人妻al黑| 夜夜骑夜夜射夜夜干| 亚洲av国产av综合av卡| 欧美成人午夜免费资源| 在线天堂最新版资源| 精品一区二区三区视频在线| 亚洲熟女精品中文字幕| av在线老鸭窝| 美女xxoo啪啪120秒动态图| av有码第一页| 美女国产高潮福利片在线看| 制服诱惑二区| 国产一区有黄有色的免费视频| 成人影院久久| 女人久久www免费人成看片| 国产精品熟女久久久久浪| 免费人妻精品一区二区三区视频| 日日啪夜夜爽| 亚洲欧洲精品一区二区精品久久久 | 两个人看的免费小视频| 欧美日韩综合久久久久久| 日韩 亚洲 欧美在线| 最黄视频免费看| 久久久久久久亚洲中文字幕| 久久人人爽av亚洲精品天堂| 国产有黄有色有爽视频| 青春草国产在线视频| 极品人妻少妇av视频| 午夜视频国产福利| 国产精品久久久久久av不卡| 美女大奶头黄色视频| 18禁动态无遮挡网站| 女人精品久久久久毛片| 国产精品国产av在线观看| 美女福利国产在线| 日韩一本色道免费dvd| 18禁观看日本| 国产亚洲午夜精品一区二区久久| 少妇高潮的动态图| 欧美日韩av久久| 26uuu在线亚洲综合色| 美女中出高潮动态图| 久久青草综合色| 满18在线观看网站| 啦啦啦啦在线视频资源| 18禁裸乳无遮挡动漫免费视频| 亚洲,欧美精品.| 美国免费a级毛片| 久久99精品国语久久久| 精品亚洲成国产av| 亚洲情色 制服丝袜| 国产精品久久久久久久电影| 欧美3d第一页| xxx大片免费视频| 精品人妻熟女毛片av久久网站| 99久国产av精品国产电影| 午夜av观看不卡| 日本欧美视频一区| 亚洲色图综合在线观看| 成人毛片a级毛片在线播放| 中文字幕制服av| 免费久久久久久久精品成人欧美视频 | 成人手机av| 国产日韩欧美亚洲二区| 黑人欧美特级aaaaaa片| 午夜福利在线观看免费完整高清在| 免费看av在线观看网站| 777米奇影视久久| 人人妻人人爽人人添夜夜欢视频| 亚洲四区av| 五月伊人婷婷丁香| 秋霞伦理黄片| 精品卡一卡二卡四卡免费| 婷婷色麻豆天堂久久| 涩涩av久久男人的天堂| 国产日韩一区二区三区精品不卡| 中文字幕人妻丝袜制服| 国产69精品久久久久777片| 精品久久蜜臀av无| 岛国毛片在线播放| 欧美xxxx性猛交bbbb| 色婷婷av一区二区三区视频| 色5月婷婷丁香| av不卡在线播放| 最黄视频免费看| 考比视频在线观看| 一级毛片我不卡| 两性夫妻黄色片 | 极品人妻少妇av视频| 欧美老熟妇乱子伦牲交| 亚洲四区av| 日韩伦理黄色片| 80岁老熟妇乱子伦牲交| 午夜老司机福利剧场| av播播在线观看一区| 日韩视频在线欧美| 男男h啪啪无遮挡| 欧美日韩成人在线一区二区| 国产成人精品一,二区| 日日摸夜夜添夜夜爱| 一级片免费观看大全| av黄色大香蕉| 久久久久久久精品精品| 一级片免费观看大全| 涩涩av久久男人的天堂| 黄色怎么调成土黄色| 22中文网久久字幕| 亚洲精品一区蜜桃| 青青草视频在线视频观看| 男女下面插进去视频免费观看 | 国产激情久久老熟女| 老司机影院毛片| 三级国产精品片| 又粗又硬又长又爽又黄的视频| 看十八女毛片水多多多| 99热全是精品| 自拍欧美九色日韩亚洲蝌蚪91| 成人亚洲精品一区在线观看| 国产精品一二三区在线看| 久久人人爽人人爽人人片va| 欧美最新免费一区二区三区| 国产无遮挡羞羞视频在线观看| 高清av免费在线| 90打野战视频偷拍视频|