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

    河口咸潮入侵研究進(jìn)展

    2014-12-04 14:39:51黃洪城匡翠萍冒小丹
    海洋科學(xué) 2014年9期
    關(guān)鍵詞:咸潮長(zhǎng)江口海平面

    黃洪城, 匡翠萍, 顧 杰, 陳 維, 冒小丹

    (1. 同濟(jì)大學(xué) 土木工程學(xué)院, 上海200092; 2. 上海海洋大學(xué) 海洋科學(xué)學(xué)院, 上海201306)

    河口地區(qū)是海洋與河流的交匯地帶, 擁有豐富的海洋資源, 交通便利, 是人類生活、生產(chǎn)、貿(mào)易活動(dòng)頻繁的場(chǎng)所。沿海及河口地區(qū)居住著全球三分之二的人口, 孕育了城市, 且其中很多發(fā)展成為重要的經(jīng)濟(jì)、貿(mào)易、文化中心[1]。然而, 長(zhǎng)期以來, 河口地區(qū)飽受咸潮災(zāi)害侵襲, 包括荷蘭、德國、意大利、美國、中國在內(nèi)的很多國家, 每年均因河口地區(qū)咸潮入侵遭受巨大損失[2]。近年來, 人類對(duì)河口地區(qū)的改造、全球氣候變化引起的河流徑流改變, 特別是全球海平面的不斷上升, 都對(duì)河口地區(qū)的咸潮災(zāi)害造成了一定的影響, 不少河口地區(qū)的咸潮災(zāi)害有頻率增大、災(zāi)害加重的趨勢(shì)。因此, 咸潮入侵作為河口水環(huán)境中普遍存在的問題, 受到越來越多的關(guān)注, 不少學(xué)者都在進(jìn)行咸潮入侵的研究。

    20世紀(jì)30年代, 美國的水道試驗(yàn)站(WES)和荷蘭Delft水工實(shí)驗(yàn)所就在咸潮研究上做了很多基礎(chǔ)工作, 而對(duì)河口或三角洲地區(qū)的咸潮研究則始于20世紀(jì)50年代左右, 美國潮汐力學(xué)委員會(huì)在1954年對(duì)哥倫比亞河口的咸潮及其現(xiàn)象首先進(jìn)行調(diào)查分析。國內(nèi)對(duì)河口咸潮的研究起步較晚, 始于20世紀(jì)60年代,到20世紀(jì)80年代才有較為系統(tǒng)的研究, 且研究區(qū)域集中在長(zhǎng)江口、珠江口等區(qū)域[3]。

    本文歸納總結(jié)了河口咸潮入侵的研究方法、影響因素及其危害的分析研究, 闡述了與咸潮入侵相關(guān)的河口基本現(xiàn)象和過程。河口咸潮入侵的數(shù)學(xué)模型經(jīng)過30多年的發(fā)展, 從最初的一維模型到后來的二維模型, 直至今天的三維模型, 不斷地豐富進(jìn)步。最后, 在總結(jié)目前研究的基礎(chǔ)上, 展望了河口咸潮入侵未來的研究方向。

    1 河口咸潮入侵的研究方法

    咸潮入侵研究的方法主要有 3種, 分別是現(xiàn)場(chǎng)觀測(cè)分析、物理模型試驗(yàn)和數(shù)學(xué)模型模擬。

    1.1 現(xiàn)場(chǎng)觀測(cè)分析

    現(xiàn)場(chǎng)觀測(cè)分析是研究咸潮入侵的基本方法, 該方法主要是通過實(shí)地觀測(cè)獲得大量鹽度數(shù)據(jù), 進(jìn)行統(tǒng)計(jì)分析, 得到一些咸潮入侵的時(shí)空變化規(guī)律, 并且可以分析徑流、潮汐等因素對(duì)咸潮入侵的影響?,F(xiàn)場(chǎng)觀測(cè)的數(shù)據(jù)同時(shí)可應(yīng)用于物理模型和數(shù)值模擬的驗(yàn)證。

    現(xiàn)場(chǎng)觀測(cè)雖然能夠獲得大量實(shí)測(cè)資料, 但其研究成果易受觀測(cè)資料的時(shí)間跨度、觀測(cè)站點(diǎn)布置、觀測(cè)指標(biāo)個(gè)數(shù)等影響, 而且現(xiàn)場(chǎng)觀測(cè)還受測(cè)量準(zhǔn)確性等人為因素和天氣等自然因素的影響。

    1.2 物理模型試驗(yàn)

    物理模型試驗(yàn)主要是通過建立研究區(qū)域的物理模型, 設(shè)計(jì)不同的河口水文情況、工況來研究河口咸潮入侵, 分析鹽度的擴(kuò)散機(jī)理及其與河口水環(huán)境、泥沙輸運(yùn)等的關(guān)系。

    然而, 該方法作為一個(gè)模型試驗(yàn)方法, 對(duì)于地形地貌復(fù)雜的大河口以及具體的工程, 往往需要進(jìn)行模型概化, 導(dǎo)致物理模型與實(shí)際河口地形地貌存在一定偏差。同時(shí), 物理模型研究咸潮入侵過程中的咸水處理工藝、咸水混合、咸界控制等也存在一定困難, 且鹽水是一種強(qiáng)電解質(zhì), 易對(duì)試驗(yàn)儀器造成腐蝕, 增加了咸潮入侵物理模型試驗(yàn)的成本。

    1.3 數(shù)學(xué)模型模擬

    數(shù)學(xué)模型模擬方法是目前研究咸潮入侵最為普遍的方法, 該方法從鹽度守恒方程出發(fā), 結(jié)合水動(dòng)力的變化, 分析不同水文條件下的咸潮入侵情況。咸潮入侵的數(shù)值模擬方法包括一維、二維和三維的數(shù)值模擬, 其中二維咸潮入侵?jǐn)?shù)值模擬又可以分為平面二維數(shù)值模擬和垂向二維數(shù)值模擬[1,4-5]。

    利用數(shù)值模擬方法研究河口咸潮入侵, 可以直觀地展現(xiàn)出河口鹽度的時(shí)空變化過程, 揭示水文條件變化后河口系統(tǒng)的響應(yīng); 采用單因子敏感性試驗(yàn),能夠清楚地了解不同動(dòng)力因子對(duì)河口咸潮入侵的影響; 而分析數(shù)學(xué)模型的計(jì)算結(jié)果, 可以進(jìn)一步研究河口咸潮入侵的機(jī)制[6]。

    數(shù)值模擬研究河口咸潮入侵能夠很好地彌補(bǔ)現(xiàn)場(chǎng)觀測(cè)易受觀測(cè)站點(diǎn)布置、測(cè)量精度、天氣影響及物理模型的高投入等不足, 這也是其成為目前主流研究方法的主要原因。

    2 河口咸潮入侵相關(guān)的物理現(xiàn)象和過程

    20 世紀(jì) 50 年代開始, 以 Pritchard[7-9]、Ippen[10]、Simmons[11]、Dyer[12]、Officer[13]等為代表, 對(duì)與咸潮入侵相關(guān)的基本河口現(xiàn)象和過程, 主要包括: 咸潮入侵長(zhǎng)度、河口的混合過程、河口咸淡水混合類型、河口密度環(huán)流、河口環(huán)流與最大渾濁帶的關(guān)系等, 進(jìn)行了富有開創(chuàng)性和奠基性的研究探索。Simmons[11]將一個(gè)潮周期內(nèi)流入河口的徑流量與潮流量之比定義為混合模數(shù)M, 來判別咸淡水混合的類型。Hansen等[14]根據(jù)實(shí)測(cè)鹽度和流速, 將河口分為 4大類: 充分混合型(余流單一向海)、部分分層型(底部余流出現(xiàn)向陸)、鹽水楔類型(具有明顯的重力環(huán)流結(jié)構(gòu))、高度分層但無明顯重力環(huán)流結(jié)構(gòu)型。Wolanski等[15]通過二維數(shù)學(xué)模型證明 Fly河口最大渾濁帶的形成與潮泵作用、床面侵蝕和沉積過程、密度流、沉降速度等因素有關(guān)。

    國內(nèi)對(duì)河口咸潮入侵的相關(guān)現(xiàn)象、過程研究主要集中在長(zhǎng)江口、珠江口等大河口。茅志昌等[16-17]依據(jù)實(shí)測(cè)資料, 對(duì)長(zhǎng)江口咸潮入侵鋒及潮汐分汊河口入侵類型等進(jìn)行了探討。張重樂等[18]利用實(shí)測(cè)資料討論了咸淡水混合及潮差、高潮位的關(guān)系。沈煥庭等[19]在對(duì)長(zhǎng)江口咸潮入侵進(jìn)行大量研究的基礎(chǔ)上,編著了《長(zhǎng)江河口鹽水入侵》, 對(duì)長(zhǎng)江口咸淡水混合類型、咸潮入侵的來源以及鹽度的時(shí)空變化規(guī)律等作了全面、深入的探討。周濟(jì)福等[20]建立準(zhǔn)二維鹽度數(shù)學(xué)模型來研究長(zhǎng)江口混合過程, 得到鹽度分布、鹽度鋒強(qiáng)度隨徑流和潮差定量變化的規(guī)律。朱建榮等[21]研究了河口最大渾濁帶形成的動(dòng)力機(jī)制。應(yīng)秩甫等[22]研究珠江口伶仃洋的咸淡水混合時(shí)發(fā)現(xiàn), 伶仃洋橫向存在東咸西淡的鹽度梯度, 鹽度為 5的等鹽度線可侵入珠江口各口門。Mao等[23]基于1998年實(shí)測(cè)資料發(fā)現(xiàn)伶仃洋東槽有明顯的密度環(huán)流。Dong等[24]研究分析了 1999~2000年實(shí)測(cè)資料, 發(fā)現(xiàn)伶仃洋東部存在明顯的重力環(huán)流, 枯季則一致向海。

    3 河口咸潮入侵的影響因素

    河口地區(qū)是河流與海洋的過渡地帶, 受到二者的雙重作用, 而入侵河口的咸潮又主要來源于海洋咸水, 因此, 徑流與潮汐是河口咸潮入侵的兩大主要影響因素。近年來, 全球海平面不斷上升, 海平面上升對(duì)河口咸潮入侵的影響逐漸變大。此外, 河口區(qū)域歷來是人類的聚居地, 經(jīng)濟(jì)發(fā)達(dá), 愈來愈頻繁的人類活動(dòng)對(duì)其產(chǎn)生的影響不斷加大。風(fēng)也是影響河口咸潮入侵的重要因素。

    3.1 徑流

    徑流是影響河口咸潮入侵的一個(gè)重要因素, 它主要通過徑流量的大小、徑流的季節(jié)變化和徑流量變化幅度大小等影響河口的咸潮入侵[19]。朱建榮等[25]、胡松等[26]利用ECOM模式研究了徑流對(duì)河口咸潮入侵的影響, 發(fā)現(xiàn)徑流量增大后, 口門內(nèi)咸潮入侵減弱, 口外鹽度減小, 沖淡水?dāng)U展范圍增大??梢? 徑流量直接影響河口的咸潮入侵, 但因?yàn)閺搅髋c海水的混合需要一定時(shí)間, 使得徑流量變化和鹽度變化通常不同步, 存在一定的時(shí)間滯后, 例如, 肖成猷等[27]在分析徑流對(duì)長(zhǎng)江口咸潮入侵影響時(shí)發(fā)現(xiàn), 長(zhǎng)江口咸潮入侵存在明顯的季節(jié)變化, 徑流量變化與鹽度變化的時(shí)間滯后在小徑流量情況下更加明顯。

    3.2 潮汐和潮流

    潮汐和潮流分別是天體引潮力引起的海面垂直方向的漲落和海水水平方向的流動(dòng), 是咸淡水混合的“動(dòng)力源”, 對(duì)咸潮入侵的影響是至關(guān)重要的。潮汐、潮流對(duì)咸潮入侵的影響包括: 潮流對(duì)咸潮的對(duì)流輸運(yùn)、潮汐引起的紊動(dòng)混合、潮汐與地形共同作用引起的“潮汐捕集”和“潮汐輸送”[28]。

    3.3 風(fēng)

    風(fēng)對(duì)咸潮入侵具有較大影響, 風(fēng)速、風(fēng)向不同,河口地區(qū)漲、落潮流的強(qiáng)度就不同, 對(duì)河口地區(qū)的咸潮入侵影響也就會(huì)有差異。不同的風(fēng)速和風(fēng)向作用下, 河口地區(qū)可以產(chǎn)生不同的水平環(huán)流[29], 所產(chǎn)生的水平環(huán)流, 可能對(duì)河口地區(qū)的咸潮入侵產(chǎn)生一定作用。例如, 朱建榮等[30]探討了風(fēng)應(yīng)力對(duì)長(zhǎng)江口沒冒沙淡水帶的影響, 得出枯季北風(fēng)產(chǎn)生向岸的??寺斶\(yùn), 生成了北港流進(jìn)南港和南槽流出的水平風(fēng)生環(huán)流, 阻礙了南槽外海的咸潮入侵。通過建立珠江口三維水動(dòng)力模型, 匡翠萍等[31]研究風(fēng)速、風(fēng)向?qū)ο募竞涌诃h(huán)流和咸淡水混合的影響, 發(fā)現(xiàn)常南風(fēng)下河口產(chǎn)生強(qiáng)分層現(xiàn)象, 強(qiáng)西南風(fēng)則加強(qiáng)了河口表層的混合, 但河口仍處于弱分層狀態(tài), 而強(qiáng)東北風(fēng)使得香港水域鹽度大大增加, 水體形成強(qiáng)混合現(xiàn)象。

    3.4 海平面上升

    近年來, 全球變暖以及日益密集的人類活動(dòng)使得海平面持續(xù)上升, 這可能使得咸潮上溯的距離增加, 加劇河口的咸潮入侵災(zāi)害。Hong 等[32]利用三維水動(dòng)力—富營養(yǎng)化數(shù)學(xué)模型(HEM-3D)對(duì)美國的Chesapeake Bay的鹽度、咸潮入侵深度、鹽度分層進(jìn)行研究, 發(fā)現(xiàn)海平面上升后, Chesapeake Bay的平均鹽度、咸潮入侵深度等呈增加趨勢(shì)。Mohsen 等[33]利用2D-FED模型研究海平面上升對(duì)尼羅河三角洲的影響, 發(fā)現(xiàn)海平面上升50 cm將使尼羅河三角洲咸潮入侵距離增加9 km。Bhuiyan等[34]利用數(shù)學(xué)模型研究了海平面上升對(duì) Gorai河咸潮入侵的影響, 發(fā)現(xiàn)海平面上升59 cm后, 河口上游80 km處的鹽度增加了0.9。

    據(jù)2011年中國海平面公報(bào)統(tǒng)計(jì), 1980年至2011年, 中國沿海海平面平均上升速率為 2.7 mm/a, 高于全球平均海平面上升速率 1.8 mm/a[35], 海平面上升對(duì)中國河口(珠江口、長(zhǎng)江口等)咸潮入侵的影響受到越來越多學(xué)者的關(guān)注??滋m等[36]建立一維動(dòng)態(tài)潮流、含氯度數(shù)學(xué)模型, 計(jì)算了海平面上升對(duì)咸潮上溯的影響, 并對(duì)珠江口代表口門在海平面上升 10, 30和60 cm后的咸潮上溯距離進(jìn)行計(jì)算。

    3.5 人類活動(dòng)

    河口地區(qū)人類活動(dòng)頻繁, 修筑擋潮閘、大量采砂、建造導(dǎo)堤、開挖深水航道等, 均對(duì)河口咸潮入侵有很大影響。其中, 不少人類活動(dòng)會(huì)加劇河口的咸潮入侵, 如南水北調(diào)東線工程若按預(yù)期調(diào)水方案實(shí)施將加劇長(zhǎng)江口的咸潮入侵, 加大長(zhǎng)江口水庫取水口附近的鹽度[37]; 長(zhǎng)江口深水航道工程增大了長(zhǎng)江口南槽和北槽上段的鹽度[29]; 珠江口人工挖沙可能破壞東江咸、淡水的回蕩平衡[38]。還有一些人類活動(dòng)則會(huì)弱化河口的咸潮入侵, 如長(zhǎng)江口南北支整治工程能減少北支鹽水倒灌, 降低陳行水庫附近鹽度[39]; 三峽水庫修建后, 枯水期水庫下泄流量增加, 有利于緩解長(zhǎng)江口水源地的咸潮入侵[37]。也就是說, 人類活動(dòng)對(duì)河口咸潮入侵的影響不一定是負(fù)面的, 也有正面的影響。

    河口地區(qū)因其地理位置優(yōu)越, 人類活動(dòng)在未來肯定會(huì)更加頻繁、劇烈, 對(duì)于河口咸潮入侵的影響將不斷擴(kuò)大。

    4 河口咸潮入侵的危害及其相關(guān)研究

    河口地區(qū)經(jīng)濟(jì)發(fā)達(dá), 人口眾多, 供應(yīng)充足的日常用水是維持河口地區(qū)正常生產(chǎn)生活的必要條件,同時(shí), 河口地區(qū)往往是大型港口所在地, 需要足夠的港口水深, 以滿足船只的正常通航。河口地區(qū)一旦遭遇咸潮入侵影響, 不僅日常用水得不到正常供應(yīng),咸潮入侵造成的泥沙絮凝淤積亦將妨礙港口的正常通航, 甚至造成航運(yùn)事故。

    4.1 河口地區(qū)供水

    中國《生活飲用水水源水質(zhì)標(biāo)準(zhǔn)》[40]規(guī)定, 飲用水氯化物含氯應(yīng)小于 250 mg/L, 當(dāng)河口地區(qū)發(fā)生咸潮入侵, 導(dǎo)致河口水源地附近鹽度增大, 將威脅河口地區(qū)居民的生活飲水用水, 如 2007~2008年珠江口強(qiáng)咸潮事件影響了廣州地區(qū)和珠海各水廠的供水[41]、2004年珠江口持續(xù)了近5個(gè)月的海水倒灌影響了1 000多萬人的飲用水[42]等。

    很多學(xué)者研究了咸潮入侵對(duì)河口地區(qū)供水的影響, 顧玉亮等[43]研究了北支咸潮入侵對(duì)長(zhǎng)江口水源地的影響; 余奕昌等[44]分析長(zhǎng)江口咸潮入侵途徑與規(guī)律, 進(jìn)而探討咸潮入侵對(duì)上海水資源的影響; 李勇等[45]分析了咸潮的成因與危害, 發(fā)現(xiàn)咸潮來臨時(shí)不僅帶來高濃度的氯化物, 同時(shí)也導(dǎo)致了下游污染物的上溯, 使得原水中其他水質(zhì)指標(biāo)上升; 鑒于咸潮入侵對(duì)河口地區(qū)供水的影響越來越大, 孔蘭等[42]對(duì)咸潮影響下磨刀門水道取淡時(shí)機(jī)進(jìn)行了研究, 提出了“最早取淡日”與“最后取淡日”等概念。

    河口地區(qū)人口密集, 經(jīng)濟(jì)發(fā)達(dá), 一旦遭遇咸潮入侵的影響, 必將對(duì)該地區(qū)居民生產(chǎn)生活用水產(chǎn)生強(qiáng)烈影響, 并帶來巨大的經(jīng)濟(jì)損失, 因而河口地區(qū)咸潮入侵對(duì)供水的影響將會(huì)受到更多學(xué)者的關(guān)注。

    4.2 泥沙絮凝淤積

    河口地區(qū)往往是大型港口所在地, 要求能夠滿足良好的通航條件, 而咸潮入侵會(huì)促使細(xì)顆粒泥沙發(fā)生絮凝沉降, 進(jìn)而造成河口泥沙淤積, 改變河口水動(dòng)力, 影響河口通航。Meade[46]、Sholkovitz[47]和Kate[48-49]等研究發(fā)現(xiàn), 鹽水內(nèi)含有大量強(qiáng)電解質(zhì),使得細(xì)顆粒泥沙在鹽水中具有物理化學(xué)粘結(jié)力, 在水流內(nèi)相互碰撞后粘附在一起, 形成較大絮凝體,加速沉降。關(guān)許為等[50-51]通過室內(nèi)試驗(yàn)證實(shí), 當(dāng)鹽度小于5時(shí), 泥沙沉積速率隨鹽度增加而加快, 而當(dāng)鹽度大于5后, 沉積速度基本與鹽度變化無關(guān)。金鷹等[52]發(fā)現(xiàn), 在徑流與咸水混合初期, 咸水中高價(jià)陽離子的吸附與離子強(qiáng)度的增大, 細(xì)顆粒泥沙的電位突降,從而大大降低了雙電層間的排斥作用, 使細(xì)顆粒泥沙趨于不穩(wěn)定。陳慶強(qiáng)等[53]認(rèn)為鹽度在影響泥沙絮凝淤積的各個(gè)因子(鹽度、含沙量、水溫、有機(jī)質(zhì)等)中占據(jù)主導(dǎo)地位。

    關(guān)于河口泥沙絮凝淤積問題已進(jìn)行了大量研究,但缺少將河口泥沙絮凝淤積與實(shí)際咸潮入侵時(shí)空分布聯(lián)系起來的系統(tǒng)研究。

    5 河口咸潮入侵?jǐn)?shù)學(xué)模型的研究進(jìn)展

    20世紀(jì)80年代, 咸潮入侵相關(guān)的數(shù)學(xué)模型開始發(fā)展。Savenije等[54-56]通過對(duì)大量河口咸潮入侵曲線的研究, 發(fā)展了一系列咸潮上溯模型, 用來模擬預(yù)測(cè)河口的咸潮。Gillibrand等[57]利用一維數(shù)學(xué)模型對(duì)Ythan河口的水位、鹽度及總氧氮進(jìn)行了模擬, 并將迎風(fēng)差分與中心差分相結(jié)合以模擬強(qiáng)潮流。Wolanski等[15]模擬研究 Fly河口的咸潮入侵情況, 探討了淺水效應(yīng)、潮波等對(duì)咸潮入侵的影響。Essink[58]采用MOCDENS3D模擬了北部Nether島嶼附近的淡水、鹽水的分布, 研究了其三維空間的咸潮入侵情況。Thain等[59]對(duì)部分混合型河口的咸潮上溯和潮汐入侵鋒進(jìn)行了詳細(xì)的研究。

    現(xiàn)今國際上使用較廣泛、發(fā)展較成熟的三維水動(dòng)力海洋數(shù)學(xué)模型主要有: POM模型[60]、ECOM模型[60]、ROMS模型[61]、TRIM 模型[62]等結(jié)構(gòu)網(wǎng)格模式, 以及 FVCOM模型[63]、ELCIRC模型[64]、UNTRIM模型[65]等無結(jié)構(gòu)網(wǎng)格模式, 它們均可應(yīng)用于河口咸潮入侵的數(shù)值模擬研究。其中ECOM模式和FVCOM模式已在許多河口海岸地區(qū)得到廣泛應(yīng)用和驗(yàn)證。還有一些可視化的商業(yè)軟件, 如Delft 3D和Mike等, 都能很好地模擬水流、咸潮入侵、泥沙輸運(yùn)等物理過程。

    國內(nèi)學(xué)者同樣建立了大量一維、二維和三維數(shù)學(xué)模型。

    一維咸潮入侵?jǐn)?shù)值模擬主要是對(duì)一維鹽度擴(kuò)散方程進(jìn)行離散和數(shù)值計(jì)算, 并利用計(jì)算結(jié)果研究鹽度縱向分布及入侵長(zhǎng)度等。黃昌筑[66]利用一維鹽度擴(kuò)散方程分析了長(zhǎng)江口咸潮入侵現(xiàn)象; 韓乃斌[67]基于一維鹽度擴(kuò)散方程研究了南水北調(diào)工程對(duì)長(zhǎng)江口咸潮入侵深度的影響; 易家豪[68]采用一維河口分汊水流數(shù)學(xué)模型研究河口縱向各斷面水流、鹽度均值的變化。

    二維咸潮入侵?jǐn)?shù)學(xué)模型包括平面和垂向二維模型, 周濟(jì)福等[20]、王義剛[69]、匡翠萍[70]等建立的長(zhǎng)江口垂向二維咸潮入侵模型, 以及肖成猷等[71]、陶學(xué)為[72]、羅小峰[73]等建立的長(zhǎng)江口平面二維咸潮入侵模型, 對(duì)長(zhǎng)江口的咸潮入侵進(jìn)行了很好的模擬研究。

    關(guān)于三維咸潮入侵模型, 匡翠萍[74]、宋元平等[75]分別建立了長(zhǎng)江口咸潮入侵的三維數(shù)學(xué)模型, 研究長(zhǎng)江口鹽度分層現(xiàn)象。嚴(yán)以新等[76]建立了河口三維非線性斜壓水流鹽度數(shù)學(xué)模型。羅小峰[73]通過三維數(shù)值模擬研究了長(zhǎng)江口深水航道工程實(shí)施后北槽咸潮入侵的變化。

    近年來, 國內(nèi)也有很多學(xué)者基于國外數(shù)學(xué)模型和商業(yè)軟件(如Delft 3D、Mike 21等)來研究河口的咸潮入侵。龔政等[77-78]基于POM模型得到適合長(zhǎng)江口的 σ坐標(biāo)下的三維非線性斜壓流場(chǎng)及鹽度數(shù)學(xué)模型。朱建榮等[79]、馬鋼鋒等[80]在ECOM模式基礎(chǔ)上,分別建立了長(zhǎng)江口三維鹽度數(shù)值模式。劉均衛(wèi)[81]利用 ELCIRC模型對(duì)長(zhǎng)江口的咸潮入侵作了研究。王彪[6]基于 FVCOM 模型建立了珠江口三維咸潮入侵?jǐn)?shù)學(xué)模型, 模型中考慮了徑流、潮汐、風(fēng)、斜壓以及陸架環(huán)流等各種動(dòng)力因子。匡翠萍等[31]利用Delft 3D模型建立珠江口水動(dòng)力模型, 研究風(fēng)速、風(fēng)向?qū)ο愀鬯颦h(huán)流和咸淡水混合的影響。

    6 總結(jié)

    (1) 河口咸潮入侵的研究方法有現(xiàn)場(chǎng)觀測(cè)分析、物理模型試驗(yàn)和數(shù)學(xué)模型模擬 3種。數(shù)值模擬研究河口咸潮入侵能夠很好地彌補(bǔ)現(xiàn)場(chǎng)觀測(cè)易受觀測(cè)站點(diǎn)布置、測(cè)量精度、天氣影響及物理模型的高投入等不足。

    (2) 與咸潮入侵相關(guān)的現(xiàn)象和過程的機(jī)理研究主要包括: 咸潮入侵長(zhǎng)度、河口的混合過程、河口咸淡水混合類型、河口密度環(huán)流、河口環(huán)流與最大渾濁帶的關(guān)系等, 對(duì)于這些基本河口現(xiàn)象、過程的研究,國外早在20世紀(jì)50年代便已經(jīng)開始, 國內(nèi)的研究則相對(duì)起步較晚, 且主要集中在長(zhǎng)江口、珠江口等大河口。

    (3) 影響河口咸潮入侵的因素包括徑流、潮汐、潮流、風(fēng)、海平面上升、人類活動(dòng)等, 近年來由于全球氣候的變化, 海平面不斷上升, 海平面上升以及人類活動(dòng)對(duì)河口咸潮入侵的影響越來越大, 當(dāng)然,人類活動(dòng)對(duì)咸潮入侵有正面和負(fù)面兩方面的影響。

    (4) 河口咸潮入侵會(huì)造成河口地區(qū)生產(chǎn)、生活供水困難、河口泥沙絮凝淤積等危害, 進(jìn)而影響河口地區(qū)的生產(chǎn)生活和經(jīng)濟(jì)發(fā)展, 學(xué)者們?cè)谙坛比肭謱?duì)河口地區(qū)供水危害方面作了很多研究, 但關(guān)于咸潮入侵對(duì)河口泥沙絮凝淤積影響方面的研究, 主要集中在鹽度對(duì)泥沙絮凝淤積影響的機(jī)理研究, 而缺少將咸潮入侵過程、時(shí)空分布與河口泥沙絮凝淤積聯(lián)系起來的系統(tǒng)研究。

    (5) 河口咸潮入侵的數(shù)學(xué)模型在國外的發(fā)展已經(jīng)相當(dāng)成熟, 國際上目前也已經(jīng)有一大批先進(jìn)的咸潮入侵?jǐn)?shù)學(xué)模型得到廣泛應(yīng)用; 國內(nèi)的學(xué)者們同樣建立了大量一維、二維和三維咸潮入侵?jǐn)?shù)學(xué)模型, 而近年來國內(nèi)不少學(xué)者利用國外先進(jìn)的咸潮入侵?jǐn)?shù)學(xué)模型和可視化的商業(yè)軟件研究河口咸潮入侵, 并取得不少成果。

    7 展望

    (1) 利用數(shù)學(xué)模型方法研究河口的咸潮入侵,伴隨著計(jì)算機(jī)功能的快速發(fā)展, 未來將受到更多學(xué)者的青睞, 物理模型和現(xiàn)場(chǎng)觀測(cè)分析是機(jī)理性方面研究和驗(yàn)證數(shù)學(xué)模型不可缺少的研究手段。

    (2) 對(duì)于咸潮入侵相關(guān)的基本河口現(xiàn)象和過程的研究, 可以從對(duì)長(zhǎng)江口、珠江口等大河口地區(qū)的重點(diǎn)關(guān)注延伸擴(kuò)展到國內(nèi)的其他中小型河口, 為解決中國整體的咸潮入侵問題服務(wù)。

    (3) 伴隨著全球氣候變化, 以及人類日益密集的生產(chǎn)生活活動(dòng), 海平面上升、人類活動(dòng)對(duì)于河口咸潮入侵的影響將越來越大, 尤其是海平面的上升,近些年全球海平面上升速度不斷加快, 海平面上升對(duì)河口咸潮入侵的影響應(yīng)受到更多關(guān)注。

    (4) 河口咸潮入侵引起的生產(chǎn)生活用水問題將繼續(xù)保持研究熱度, 對(duì)河口咸潮入侵引起的泥沙絮凝淤積問題, 應(yīng)該加強(qiáng)泥沙絮凝淤積機(jī)理與咸潮入侵過程、時(shí)空分布之間關(guān)系的系統(tǒng)研究。

    [1] 楊莉玲. 河口鹽水入侵的數(shù)值模擬研究[D]. 上海:上海交通大學(xué), 2007.

    [2] 張更生. 海水入侵機(jī)理及防治措施的三維數(shù)值模擬[D]. 大連: 大連海事大學(xué), 2007.

    [3] 彭靖. 磨刀門咸潮上溯的數(shù)值模擬及其機(jī)理分析[D].廣州: 中山大學(xué), 2007.

    [4] 胡溪. 珠江口磨刀門水道咸潮入侵?jǐn)?shù)值模擬研究[D].北京: 清華大學(xué), 2010.

    [5] 張宗坤. 河口水域鹽水入侵的試驗(yàn)研究和數(shù)值模擬[D]. 太原: 太原理工大學(xué), 2011.

    [6] 王彪. 珠江河口鹽水入侵[D]. 上海: 華東師范大學(xué), 2011.

    [7] Pritchard D W. A study of the salt balance of a coastal plain estuary[J]. Marine Sciences, 1954, 13: 133-144.

    [8] Pritchard D W. Estuarine Hydrography[J]. Advances in Geophysics, 1952, 1: 243-280.

    [9] Pritchard D W. Salinity distribution and circulation in the Chesapeake Bay estuarine system[J]. Marine Research, 1952, 11: 106-123.

    [10] Ippen A T, Harleman D R F. One-dimensional analysis of salinity intrusion in estuaries: Technical Bulletin number 5[R]. U.S.A: Committee on Tidal Hydraulics, 1961.

    [11] Simmons H B. Salinity effects on estuarine hydraulics and sedimentation[C]. Tanaka H. Process of 13th IAHR (3).Kyoto, Japan: International Association for Hydro-Environment Engineering and Research, 1969: 311-325.

    [12] Dyer K R. Estuaries and estuarine sedimentation[M].London: Cambridge University Press, 1979.

    [13] Officer C B. Physical oceanography of estuaries [M].New York: Academic Press, 1976.

    [14] Hansen D V, Maurice R J. New dimensions in estuary classification[J]. Limnology and Oceanography, 1966,1: 319-326.

    [15] Wolanski E, King B, Galloway D. Salinity intrusion in the Fly River Estuary, Papua New Guinea[J]. Journal of Coastal Research, 1997, 13(4): 983-994.

    [16] 茅志昌, 沈煥庭. 潮汐分汊河口鹽水入侵類型探討——以長(zhǎng)江口為例[J]. 華東師范大學(xué)學(xué)報(bào)(自然科學(xué)版),1995, 2: 77-85.

    [17] 茅志昌. 長(zhǎng)江河口鹽水入侵鋒研究[J]. 海洋與湖沼,1995, 26(6): 643-649.

    [18] 張重樂, 沈煥庭. 長(zhǎng)江口咸淡水混合及其對(duì)懸沙的影響[J]. 華東師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 1988, 4: 83-88.

    [19] 沈煥庭, 茅志昌, 朱建榮. 長(zhǎng)江河口鹽水入侵[M].北京: 海洋出版社, 2003.

    [20] 周濟(jì)福, 劉青泉, 李家春. 河口混合過程的研究[J].中國科學(xué)(A輯), 1999, 29(9): 835-843.

    [21] 朱建榮, 傅德健, 吳輝, 等. 河口最大渾濁帶形成的動(dòng)力模式和數(shù)值試驗(yàn)[J]. 海洋工程, 2004, 22(1): 66-73.

    [22] 應(yīng)秩甫, 陳世光. 珠江口伶仃洋咸淡水混合特征[J].海洋學(xué)報(bào), 1983, 5(1): 1-10.

    [23] Mao Q, Shi P, Yin K, et al. Tides and tidal currents in the Pearl River Estuary[J]. Continental Shelf Research,2004, 24(16): 1797-1808.

    [24] Dong L, Su J, Ah Wong L, et al. Seasonal variation and dynamics of the Pearl River plume[J]. Continental Shelf Research, 2004, 24(16): 1761-1777.

    [25] 朱建榮, 胡松, 傅得健, 等. 河口環(huán)流和鹽水入侵I——模式及控制數(shù)值試驗(yàn)[J]. 青島海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2003, 33(2): 180-184.

    [26] 胡松, 朱建榮, 傅得健, 等. 河口環(huán)流和鹽水入侵Ⅱ——徑流量和海平面上升的影響[J]. 青島海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2003, 33(3): 337-342.

    [27] 肖成猷, 沈煥庭. 長(zhǎng)江河口鹽水入侵影響因子分析[J]. 華東師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 1998, 3: 74-80.

    [28] Fischer H B. Mixing in Inland and Coastal Waters[M].New York: Academic Press, 1979.

    [29] 吳輝. 長(zhǎng)江河口鹽水入侵研究[D].上海: 華東師范大學(xué), 2006.

    [30] 朱建榮, 傅利輝, 吳輝. 風(fēng)應(yīng)力和科氏力對(duì)長(zhǎng)江河口沒冒沙淡水帶的影響[J]. 華東師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2008, 6: 1-8.

    [31] Kuang C P, Lee J H W, Harrison P J, et al. Effect of wind speed and direction on summer tidal circulation and vertical mixing in Hong Kong waters[J]. Journal of Coastal Research, 2011, 27(6A): 74-86.

    [32] Hong B, Shen J. Responses of estuarine salinity and transport processes to potential future sea level rise in the Chesapeake Bay[J]. Estuarine, Coastal and Shelf Science, 2012, 104: 33-45.

    [33] Mohsen M S, Vijay P S. Effect of climate change on sea water intrusion in coastal aquifers[J]. Hydrological Processes, 1999, 13: 1277-1287.

    [34] Bhuiyan M J A N, Dutta D. Assessing impacts of sea level rise on river salinity in the Gorai river network,Bangladesh[J]. Estuarine, Coastal and Shelf Science,2012, 96: 219-227.

    [35] 國家海洋局, 2012年中國海平面公報(bào)[EB/OL].[2012-12-21]. http: //www.coi.gov.cn/gongbao/haipingmian/201303/t20130308_26217.html .

    [36] 孔蘭, 陳曉宏, 杜建, 等. 基于數(shù)學(xué)模型的海平面上升對(duì)咸潮上溯的影響[J]. 自然資源學(xué)報(bào), 2010, 25(7):1097-1104.

    [37] 孫波, 劉曙光, 顧杰, 等. 三峽與南水北調(diào)工程對(duì)長(zhǎng)江口水源地的影響[J]. 人民長(zhǎng)江, 2008, 39(16): 4-7.

    [38] 韓龍喜, 李偉, 陸永軍, 等. 人工挖沙對(duì)東江水系水動(dòng)力及環(huán)境影響分析[J]. 河海大學(xué)學(xué)報(bào)(自然科學(xué)版),2005, 33(2): 123-126.

    [39] 李禔來, 李誼純, 高祥宇, 等. 長(zhǎng)江口整治工程對(duì)鹽水入侵影響研究[J]. 海洋工程, 2005, 23(3): 31-38.

    [40] CJ 3020-1993, 生活飲用水水源水質(zhì)標(biāo)準(zhǔn)[S].

    [41] 羅琳, 陳舉, 楊威, 等. 2007-2008年冬季珠江三角洲強(qiáng)咸潮事件[J]. 熱帶海洋學(xué)報(bào), 2010, 29(6): 22-28.

    [42] 孔蘭, 陳曉宏, 劉斌, 等. 咸潮影響下磨刀門水道取淡時(shí)機(jī)初探[J]. 水資源保護(hù), 2011, 27(6): 24-27.

    [43] 顧玉亮, 吳守培, 樂勤. 北支鹽水入侵對(duì)長(zhǎng)江口水源地影響研究[J]. 人民長(zhǎng)江, 2003, 34(4): 1-3.

    [44] 余奕昌, 陳景山. 咸潮入侵對(duì)上海市水質(zhì)影響研究[J]. 華東師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 1991, 3: 70-78.

    [45] 李勇, 陳超, 張曉健, 等. 咸潮對(duì)城鎮(zhèn)供水水質(zhì)的影響及對(duì)策[J]. 水利水電技術(shù), 2008, 39(10): 21-23.

    [46] Meade R H. Transport and deposition of sediments in estuaries[J]. Mem Geol Soc Am, 1972, 133: 91-119.

    [47] Sholkovitz E R. Flocculation of dissolved organic and inorganic matter during the mixing of river water and seawater[J]. Geochimica et Cosmochimica Acta, 1976,40(7): 831-845.

    [48] Kate K. Sediment deposition flocculated suspensions[J].Sedimentology, 1975, 22(1): 111-123.

    [49] Kate K. Particulate matter grain-size characteristics and flocculation in a partially mixed estuary[J]. Sedimentology, 1981, 28(1): 107-114.

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

    [51] 關(guān)許為, 陳英祖. 長(zhǎng)江口泥沙絮凝靜水沉降動(dòng)力學(xué)模式的試驗(yàn)研究[J]. 海洋工程, 1995, 13(1): 46-50.

    [52] 金鷹, 王義剛, 李宇. 長(zhǎng)江口粘性細(xì)顆粒泥沙絮凝試驗(yàn)研究[J]. 河海大學(xué)學(xué)報(bào)(自然科學(xué)版), 2002, 30(3): 61-63.

    [53] 陳慶強(qiáng), 孟翊, 周菊珍, 等. 長(zhǎng)江口細(xì)顆粒泥沙絮凝作用及其制約因素研究[J].海洋工程, 2005, 23(1): 74-82.

    [54] Savenije H H G. Predictive model for salt intrusion in estuaries[J]. Journal of Hydrology, 1993, 148(1-4): 203-218.

    [55] Savenije H H G. Rapid assessment technique for salt intrusion in alluvial estuaries[D]. Netherland: Delft University of Technology, 1992.

    [56] Savenije H H G. A one-dimensional model for salinity intrusion in alluvial estuaries[J]. Journal of Hydrology,1986, 85(1): 87-109.

    [57] Gillibrand P A, Balls P W. Modelling salt intrusion and nitrate concentrations in the Ythan Estuary[J]. Estuarine,Coastal and Shelf Science, 1998, 47(6): 695-706.

    [58] Essink O G H P. Saltwater intrusion in 3D large-scale aquifers: a dutch case[J]. Physics and Chemistry of the Earth, Part B: Hydrology, Oceans and Atmosphere,2001, 26(4): 337-344.

    [59] Thain R H, Priestley A D, Davidson M A. The formation of a tidal intrusion front at the mouth of a macrotidal, partially mixed estuary: a field study of the Dart estuary, UK[J]. Estuarine, Coastal and Shelf Science, 2004, 61(1): 161-172.

    [60] Blumberg A F, Mellor G L. A description of a three-dimensional coastal ocean circulation model[R].American Geophysical Union. 1987.

    [61] Shchepetkin A F, Mcwilliams J C. The regional oceanic modeling system (ROMS): a split-explicit, free-surface,topography-following-coordinate oceanic model[J].Ocean Modelling, 2005, 9(4): 347-404.

    [62] Casulli V, Cheng R T. Semi-implicit finite difference methods for three-dimensional shallow water flow[J].Int J Numer Methods Fluids , 1992, 15: 629-648.

    [63] Chen C S, Liu H, Beardsley R C. An unstructured,finite-volume, three-dimensional, primitive equation ocean model: application to coastal ocean and estuaries[J]. Journal of Atmospheric and Oceanic Technology,2003, 20: 159-186.

    [64] Zhang Y, Baptista A M, Myers E P. A cross-scale model for 3D baroclinic circulation in estuary–plume–shelf systems: I. Formulation and skill assessment[J].Continental Shelf Research, 2004, 24(18): 2187-2214.

    [65] Casulli V, Walters R A. An unstructured grid, threedimensional model based on the shallow water equations[J]. Int J Numer Methods Fluids, 2000, 32: 331-348.

    [66] 黃昌筑. 長(zhǎng)江口鹽水入侵及其對(duì)河口攔門沙的作用[D].南京: 河海大學(xué), 1982.

    [67] 韓乃斌. 南水北調(diào)對(duì)長(zhǎng)江口鹽水入侵影響的預(yù)測(cè)[J].地理研究, 1983, 2(2): 99-107.

    [68] 易家豪. 長(zhǎng)江口南水北調(diào)鹽水模型計(jì)算研究[R].南京:南京水利科學(xué)研究院, 1987.

    [69] 王義剛. 河口鹽水入侵垂向二維數(shù)值計(jì)算[D]. 南京:河海大學(xué), 1989.

    [70] 匡翠萍. 長(zhǎng)江口攔門沙沖淤及懸沙沉降規(guī)律研究和水流鹽度泥沙數(shù)學(xué)模型[D]. 南京: 南京水利科學(xué)研究院, 1993.

    [71] 肖成猷, 朱建榮, 沈煥庭. 長(zhǎng)江口北支鹽水倒灌的數(shù)值模型研究[J]. 海洋學(xué)報(bào), 2000, 22(5): 124-132.

    [72] 陶學(xué)為. 長(zhǎng)江口海水入侵研究[J]. 水利學(xué)報(bào), 1991, 9:36-41.

    [73] 羅小峰. 長(zhǎng)江口水流鹽度數(shù)值模擬[D]. 南京: 南京水利科學(xué)研究院, 2003.

    [74] 匡翠萍. 長(zhǎng)江口鹽水入侵三維數(shù)值模擬[J]. 河海大學(xué)學(xué)報(bào), 1997, 25(4): 56-62.

    [75] 宋元平, 胡方西, 谷國傳, 等. 長(zhǎng)江口口外海濱鹽度擴(kuò)散的分層數(shù)學(xué)模型[J]. 華東師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 1990, 4: 74-84.

    [76] 嚴(yán)以新, 宋志堯. 長(zhǎng)江口深水航道整治期工程二、三期三維數(shù)值研究[R]. 南京: 河海大學(xué), 1998.

    [77] 龔政. 長(zhǎng)江口三維斜壓流場(chǎng)及鹽度場(chǎng)數(shù)值模擬[D].南京: 河海大學(xué), 2002.

    [78] 龔政, 張長(zhǎng)寬, 張東生, 等. 長(zhǎng)江口正壓、斜壓診斷及斜壓預(yù)報(bào)模式——三維流場(chǎng)數(shù)值模擬[J]. 海洋工程, 2004, 22(2): 39-45.

    [79] 朱建榮, 朱首賢. ECOM模式的改進(jìn)及在長(zhǎng)江河口、杭州灣及鄰近海區(qū)的應(yīng)用[J]. 海洋與湖沼, 2003,34(4): 364-374.

    [80] 馬鋼峰, 劉曙光, 戚定滿. 長(zhǎng)江口鹽水入侵?jǐn)?shù)值模型研究[J]. 水動(dòng)力學(xué)研究與進(jìn)展(A輯), 2006, 21(1):53-61.

    [81] 劉均衛(wèi). ELCIRC模型在長(zhǎng)江口鹽水入侵研究中的應(yīng)用[D]. 南京: 河海大學(xué), 2008.

    猜你喜歡
    咸潮長(zhǎng)江口海平面
    天文大潮+寒潮大風(fēng)為何會(huì)加劇咸潮
    冰山熔化會(huì)使海平面上升嗎
    海平面變化背景下三大河口咸潮入侵特征及變化淺析
    海平面上升 我們?nèi)绾螒?yīng)對(duì)
    中國海平面比去年升高38毫米
    長(zhǎng)江口鳳鱭繁殖群體的年齡結(jié)構(gòu)和生長(zhǎng)特性
    長(zhǎng)江口及其鄰近水域仔稚魚種類組成及分布特征
    珠江口咸潮影響因素分析
    珠江口咸潮影響因素分析
    風(fēng)場(chǎng)對(duì)夏季長(zhǎng)江口外東北海域低鹽水團(tuán)的作用
    精品熟女少妇八av免费久了| 日韩精品免费视频一区二区三区| 一本精品99久久精品77| 欧美性猛交╳xxx乱大交人| 女同久久另类99精品国产91| 欧美日韩亚洲国产一区二区在线观看| 精品久久久久久,| 中文在线观看免费www的网站 | 国产人伦9x9x在线观看| 国产精品爽爽va在线观看网站| 久久精品人妻少妇| 国产熟女xx| 午夜福利在线观看吧| 亚洲激情在线av| 夜夜看夜夜爽夜夜摸| 真人做人爱边吃奶动态| 一级作爱视频免费观看| 久久久久久久午夜电影| 日韩大码丰满熟妇| 日韩 欧美 亚洲 中文字幕| 国产精品一区二区三区四区免费观看 | 曰老女人黄片| 国产久久久一区二区三区| 最近最新免费中文字幕在线| 国产精品av视频在线免费观看| 国产99白浆流出| 久久 成人 亚洲| 97超级碰碰碰精品色视频在线观看| 亚洲aⅴ乱码一区二区在线播放 | 99久久精品国产亚洲精品| 免费看美女性在线毛片视频| 亚洲美女黄片视频| 又黄又爽又免费观看的视频| 欧美一级a爱片免费观看看 | 国产精品乱码一区二三区的特点| 成人18禁在线播放| 在线观看免费日韩欧美大片| 日韩三级视频一区二区三区| 村上凉子中文字幕在线| 女同久久另类99精品国产91| 亚洲av成人一区二区三| 精品欧美一区二区三区在线| 午夜免费观看网址| 亚洲国产欧洲综合997久久,| 日本成人三级电影网站| 国产精品1区2区在线观看.| 国产高清视频在线播放一区| 亚洲18禁久久av| 中亚洲国语对白在线视频| 午夜免费激情av| 欧美另类亚洲清纯唯美| 亚洲欧美精品综合一区二区三区| 国产精品综合久久久久久久免费| 精品国产美女av久久久久小说| 国产成人精品久久二区二区91| 精品欧美一区二区三区在线| 国产精品99久久99久久久不卡| 久久人人精品亚洲av| 淫秽高清视频在线观看| 久久久精品国产亚洲av高清涩受| 搡老熟女国产l中国老女人| 精品少妇一区二区三区视频日本电影| 天堂动漫精品| 十八禁人妻一区二区| 欧美中文综合在线视频| 日韩中文字幕欧美一区二区| 一级a爱片免费观看的视频| 欧美一区二区国产精品久久精品 | 久久久久久人人人人人| 黄色视频不卡| 两性夫妻黄色片| 成人18禁高潮啪啪吃奶动态图| 全区人妻精品视频| 国内揄拍国产精品人妻在线| 亚洲av电影在线进入| 欧美三级亚洲精品| 99热只有精品国产| 免费在线观看黄色视频的| 99精品在免费线老司机午夜| 久久久久久人人人人人| 国产一区二区三区在线臀色熟女| 亚洲第一欧美日韩一区二区三区| av免费在线观看网站| 黄色女人牲交| 一边摸一边抽搐一进一小说| 中文字幕久久专区| 又粗又爽又猛毛片免费看| 88av欧美| 99热6这里只有精品| 搡老熟女国产l中国老女人| 欧美av亚洲av综合av国产av| 亚洲五月天丁香| 黄色丝袜av网址大全| 女同久久另类99精品国产91| 18禁国产床啪视频网站| 国产真实乱freesex| 一级毛片精品| 少妇人妻一区二区三区视频| 欧美成人一区二区免费高清观看 | 亚洲色图av天堂| 欧美黑人精品巨大| 久久午夜亚洲精品久久| 韩国av一区二区三区四区| 国产精品影院久久| 99久久无色码亚洲精品果冻| 欧美一级a爱片免费观看看 | 日本在线视频免费播放| 日本黄色视频三级网站网址| 露出奶头的视频| 久久久久性生活片| 手机成人av网站| 久久婷婷人人爽人人干人人爱| 日韩av在线大香蕉| √禁漫天堂资源中文www| 三级男女做爰猛烈吃奶摸视频| 他把我摸到了高潮在线观看| 久久香蕉激情| 国产成+人综合+亚洲专区| 在线视频色国产色| 国产精品久久视频播放| 真人做人爱边吃奶动态| 国产在线观看jvid| 俺也久久电影网| 99热这里只有是精品50| 伊人久久大香线蕉亚洲五| 一个人免费在线观看的高清视频| 精品欧美一区二区三区在线| 亚洲中文字幕一区二区三区有码在线看 | 久久久久久久久中文| 精品人妻1区二区| 国产伦人伦偷精品视频| 中文字幕最新亚洲高清| ponron亚洲| 亚洲av电影在线进入| 国产精品亚洲av一区麻豆| 欧美 亚洲 国产 日韩一| av欧美777| 1024视频免费在线观看| 人妻丰满熟妇av一区二区三区| 免费在线观看影片大全网站| 免费看美女性在线毛片视频| 伦理电影免费视频| 亚洲全国av大片| 露出奶头的视频| 国产v大片淫在线免费观看| 亚洲午夜精品一区,二区,三区| 国产精品av视频在线免费观看| 午夜免费激情av| 成人欧美大片| 亚洲国产精品999在线| 亚洲国产精品成人综合色| 国产aⅴ精品一区二区三区波| 久久精品91蜜桃| 国产精品一区二区三区四区久久| 国产97色在线日韩免费| 在线观看免费视频日本深夜| 九九热线精品视视频播放| 国产69精品久久久久777片 | 久久99热这里只有精品18| 久久久久九九精品影院| 在线观看舔阴道视频| 国产在线观看jvid| 午夜精品在线福利| 韩国av一区二区三区四区| 午夜福利欧美成人| 麻豆成人午夜福利视频| 成年版毛片免费区| 黄色视频不卡| 成人18禁在线播放| 成人国产一区最新在线观看| 欧美黑人巨大hd| 淫妇啪啪啪对白视频| 人人妻,人人澡人人爽秒播| 日日干狠狠操夜夜爽| 高清在线国产一区| 波多野结衣高清无吗| 在线观看美女被高潮喷水网站 | 午夜免费成人在线视频| 亚洲精品久久国产高清桃花| 日韩欧美在线乱码| 久久精品国产亚洲av香蕉五月| 一本一本综合久久| 亚洲va日本ⅴa欧美va伊人久久| 久久亚洲精品不卡| 日韩有码中文字幕| 国产精品乱码一区二三区的特点| 天堂动漫精品| 99热这里只有精品一区 | 日本五十路高清| 欧美zozozo另类| 成熟少妇高潮喷水视频| 一进一出抽搐gif免费好疼| 久久中文字幕人妻熟女| 制服人妻中文乱码| 亚洲精品一区av在线观看| 欧美zozozo另类| 50天的宝宝边吃奶边哭怎么回事| 成人三级黄色视频| 亚洲欧美日韩高清专用| 亚洲国产精品合色在线| 久久香蕉国产精品| 日本三级黄在线观看| 不卡一级毛片| 久久精品国产亚洲av高清一级| 一本精品99久久精品77| 亚洲乱码一区二区免费版| 亚洲五月天丁香| 精品久久久久久久人妻蜜臀av| 成人三级黄色视频| 女警被强在线播放| 国产精品免费视频内射| 亚洲欧美日韩高清在线视频| 嫩草影视91久久| 岛国在线免费视频观看| 黄色毛片三级朝国网站| 国产三级在线视频| 久久中文字幕一级| 国产av一区在线观看免费| 精品电影一区二区在线| 日本免费a在线| 一本精品99久久精品77| 亚洲真实伦在线观看| 久久精品成人免费网站| 黄色毛片三级朝国网站| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲电影在线观看av| 叶爱在线成人免费视频播放| 丰满人妻熟妇乱又伦精品不卡| 午夜福利高清视频| 日日干狠狠操夜夜爽| 免费在线观看日本一区| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久国内视频| 成人精品一区二区免费| avwww免费| 成人av在线播放网站| 亚洲欧美日韩东京热| 99精品欧美一区二区三区四区| 两人在一起打扑克的视频| 日本成人三级电影网站| 天天添夜夜摸| 国产aⅴ精品一区二区三区波| 亚洲精品一区av在线观看| 嫁个100分男人电影在线观看| 人人妻人人看人人澡| 男女做爰动态图高潮gif福利片| 亚洲人成电影免费在线| 国产精品国产高清国产av| 国产av不卡久久| 精品久久久久久成人av| 日韩中文字幕欧美一区二区| 制服丝袜大香蕉在线| xxx96com| 人人妻人人澡欧美一区二区| 亚洲精华国产精华精| 高潮久久久久久久久久久不卡| 国产精品久久久久久亚洲av鲁大| 色播亚洲综合网| 啦啦啦观看免费观看视频高清| 精品人妻1区二区| 久久久久国内视频| 黄色a级毛片大全视频| 国产伦在线观看视频一区| 亚洲国产精品成人综合色| 视频区欧美日本亚洲| 久久香蕉国产精品| av超薄肉色丝袜交足视频| 狠狠狠狠99中文字幕| 1024手机看黄色片| 国产精品一区二区免费欧美| 国产黄片美女视频| 日韩欧美国产一区二区入口| 国产欧美日韩一区二区三| 国产精品久久视频播放| 精品国产乱子伦一区二区三区| 国产午夜福利久久久久久| 国产精品久久久人人做人人爽| 亚洲aⅴ乱码一区二区在线播放 | 久久久久精品国产欧美久久久| 国产精品国产高清国产av| 久久欧美精品欧美久久欧美| 久久久水蜜桃国产精品网| 男人的好看免费观看在线视频 | 日本免费一区二区三区高清不卡| 亚洲熟妇中文字幕五十中出| 青草久久国产| 夜夜躁狠狠躁天天躁| 999久久久国产精品视频| 九色成人免费人妻av| 欧美日韩一级在线毛片| 怎么达到女性高潮| 黄色 视频免费看| 麻豆国产97在线/欧美 | 亚洲天堂国产精品一区在线| 久久国产精品影院| 国产91精品成人一区二区三区| 久9热在线精品视频| 99国产精品一区二区三区| 美女黄网站色视频| 久久久国产成人精品二区| 亚洲电影在线观看av| 亚洲欧美一区二区三区黑人| 成人一区二区视频在线观看| 男女下面进入的视频免费午夜| 一进一出抽搐动态| av福利片在线观看| 久久久精品欧美日韩精品| 级片在线观看| 97人妻精品一区二区三区麻豆| 精品欧美国产一区二区三| 亚洲一区中文字幕在线| 色综合站精品国产| 校园春色视频在线观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲aⅴ乱码一区二区在线播放 | 淫秽高清视频在线观看| 成人18禁高潮啪啪吃奶动态图| 欧美日本亚洲视频在线播放| 亚洲成av人片免费观看| 女人爽到高潮嗷嗷叫在线视频| 男女做爰动态图高潮gif福利片| 欧美丝袜亚洲另类 | 久久精品国产清高在天天线| 亚洲精品国产一区二区精华液| 美女大奶头视频| 中文字幕精品亚洲无线码一区| 国产99白浆流出| 色精品久久人妻99蜜桃| 亚洲最大成人中文| 国产午夜精品久久久久久| 久久热在线av| 国产av一区二区精品久久| 国产v大片淫在线免费观看| 国产三级在线视频| 久久久久性生活片| 日韩欧美国产一区二区入口| 日本一区二区免费在线视频| 国产午夜福利久久久久久| 亚洲精品久久国产高清桃花| 欧美精品亚洲一区二区| 精品免费久久久久久久清纯| 一本综合久久免费| 亚洲av电影在线进入| 亚洲av中文字字幕乱码综合| 亚洲18禁久久av| 亚洲人成网站高清观看| 精品乱码久久久久久99久播| 国产成人啪精品午夜网站| 国产精品香港三级国产av潘金莲| 99久久精品国产亚洲精品| 午夜福利高清视频| av福利片在线| 成人特级黄色片久久久久久久| 日本 欧美在线| 精品福利观看| 成人午夜高清在线视频| 亚洲 欧美一区二区三区| 久久久国产欧美日韩av| videosex国产| 最近最新免费中文字幕在线| 免费在线观看日本一区| 中国美女看黄片| 成人欧美大片| 国内揄拍国产精品人妻在线| 真人一进一出gif抽搐免费| 国产午夜精品论理片| 亚洲国产精品999在线| 国产精华一区二区三区| 窝窝影院91人妻| 亚洲片人在线观看| 麻豆成人午夜福利视频| 亚洲专区字幕在线| 88av欧美| 成人三级做爰电影| 亚洲精品在线美女| 国产高清有码在线观看视频 | 成年版毛片免费区| 午夜免费观看网址| 国产久久久一区二区三区| 午夜影院日韩av| 久久久久性生活片| 久久欧美精品欧美久久欧美| 亚洲欧美激情综合另类| av天堂在线播放| 成年免费大片在线观看| av天堂在线播放| 国产亚洲欧美98| 男女午夜视频在线观看| 首页视频小说图片口味搜索| 丁香欧美五月| 夜夜夜夜夜久久久久| 人人妻人人澡欧美一区二区| 首页视频小说图片口味搜索| 美女免费视频网站| 欧美另类亚洲清纯唯美| 国产av一区二区精品久久| 午夜久久久久精精品| 美女高潮喷水抽搐中文字幕| 精品久久久久久久久久免费视频| 国产欧美日韩精品亚洲av| 可以在线观看毛片的网站| 欧美日韩中文字幕国产精品一区二区三区| 看黄色毛片网站| av在线天堂中文字幕| 国产亚洲精品一区二区www| 搡老妇女老女人老熟妇| 亚洲欧美日韩东京热| 啪啪无遮挡十八禁网站| 国产精品99久久99久久久不卡| 在线观看午夜福利视频| 国产99白浆流出| 精品熟女少妇八av免费久了| 18禁观看日本| 99热6这里只有精品| 亚洲成人精品中文字幕电影| xxxwww97欧美| 搡老岳熟女国产| √禁漫天堂资源中文www| 18禁国产床啪视频网站| 精品国产乱子伦一区二区三区| 亚洲国产欧美人成| 国产熟女xx| 免费观看人在逋| 婷婷亚洲欧美| 免费搜索国产男女视频| 亚洲专区国产一区二区| 欧美乱妇无乱码| 国产一级毛片七仙女欲春2| 最近在线观看免费完整版| 好男人电影高清在线观看| 免费在线观看亚洲国产| 亚洲自拍偷在线| 99热这里只有精品一区 | 嫩草影院精品99| 在线观看舔阴道视频| 亚洲精品粉嫩美女一区| 观看免费一级毛片| 欧美黑人精品巨大| 亚洲男人的天堂狠狠| 99国产极品粉嫩在线观看| 成人av在线播放网站| 日本成人三级电影网站| 男女视频在线观看网站免费 | 久久人妻福利社区极品人妻图片| 亚洲av电影不卡..在线观看| 男人舔奶头视频| 一本一本综合久久| 国产精品1区2区在线观看.| 午夜福利免费观看在线| 在线观看免费午夜福利视频| 亚洲人成网站在线播放欧美日韩| 少妇裸体淫交视频免费看高清 | 午夜免费成人在线视频| 好男人在线观看高清免费视频| 精品不卡国产一区二区三区| 久久精品国产亚洲av高清一级| 欧美日韩国产亚洲二区| 亚洲国产欧美人成| 日韩大尺度精品在线看网址| 国产私拍福利视频在线观看| 成人亚洲精品av一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 在线看三级毛片| www.自偷自拍.com| 在线观看日韩欧美| 悠悠久久av| 最近在线观看免费完整版| 成人特级黄色片久久久久久久| a级毛片在线看网站| 白带黄色成豆腐渣| 亚洲av熟女| 无限看片的www在线观看| 亚洲人成77777在线视频| 人人妻人人看人人澡| 久久久国产精品麻豆| 久久精品亚洲精品国产色婷小说| 免费看十八禁软件| 悠悠久久av| 女人被狂操c到高潮| av在线播放免费不卡| 久久午夜综合久久蜜桃| 国产精品爽爽va在线观看网站| 欧美绝顶高潮抽搐喷水| 亚洲五月婷婷丁香| 女人高潮潮喷娇喘18禁视频| 亚洲第一电影网av| 亚洲精品国产精品久久久不卡| 久久精品亚洲精品国产色婷小说| 18禁裸乳无遮挡免费网站照片| 12—13女人毛片做爰片一| 国产成人欧美在线观看| 精品国产美女av久久久久小说| 国产av麻豆久久久久久久| 一夜夜www| 亚洲 欧美 日韩 在线 免费| 欧美色欧美亚洲另类二区| 国产精品香港三级国产av潘金莲| 看免费av毛片| 香蕉丝袜av| 国产视频一区二区在线看| 午夜激情福利司机影院| 久久婷婷人人爽人人干人人爱| 亚洲精品av麻豆狂野| 两个人视频免费观看高清| 精品不卡国产一区二区三区| 很黄的视频免费| 男人舔女人下体高潮全视频| 国内精品久久久久精免费| 国产精品自产拍在线观看55亚洲| 精品少妇一区二区三区视频日本电影| 久久久精品国产亚洲av高清涩受| 久久亚洲精品不卡| 天堂av国产一区二区熟女人妻 | 欧美日韩国产亚洲二区| 桃红色精品国产亚洲av| 国产亚洲精品第一综合不卡| www.熟女人妻精品国产| 亚洲av电影在线进入| 一二三四在线观看免费中文在| 麻豆国产av国片精品| 亚洲成人中文字幕在线播放| 亚洲精品色激情综合| 日韩欧美在线乱码| 欧美成人一区二区免费高清观看 | 露出奶头的视频| 亚洲欧美精品综合一区二区三区| 午夜老司机福利片| 他把我摸到了高潮在线观看| 91在线观看av| 男人的好看免费观看在线视频 | 亚洲欧洲精品一区二区精品久久久| www.999成人在线观看| 不卡av一区二区三区| 国产av又大| 18禁黄网站禁片免费观看直播| 精品第一国产精品| 国产精品自产拍在线观看55亚洲| 一边摸一边做爽爽视频免费| xxxwww97欧美| 丰满人妻一区二区三区视频av | 国产精品亚洲av一区麻豆| 亚洲,欧美精品.| 国产成+人综合+亚洲专区| 给我免费播放毛片高清在线观看| 2021天堂中文幕一二区在线观| 国产又黄又爽又无遮挡在线| 嫩草影院精品99| 熟女电影av网| 国产成人欧美在线观看| 国产欧美日韩精品亚洲av| 久久久久国产精品人妻aⅴ院| 亚洲人成电影免费在线| 婷婷六月久久综合丁香| 啦啦啦观看免费观看视频高清| 欧美成人性av电影在线观看| 悠悠久久av| 国产精品久久电影中文字幕| 亚洲av日韩精品久久久久久密| 老熟妇乱子伦视频在线观看| 在线观看日韩欧美| 久久婷婷人人爽人人干人人爱| 亚洲av片天天在线观看| 桃红色精品国产亚洲av| 丁香欧美五月| 国产av麻豆久久久久久久| a级毛片在线看网站| 日本五十路高清| 国产在线精品亚洲第一网站| 亚洲欧美一区二区三区黑人| 我的老师免费观看完整版| 校园春色视频在线观看| 欧美性猛交黑人性爽| 国产亚洲av高清不卡| 久久久国产成人精品二区| АⅤ资源中文在线天堂| 日本熟妇午夜| 久久性视频一级片| 亚洲国产精品999在线| 成人一区二区视频在线观看| 亚洲国产欧美一区二区综合| 女人高潮潮喷娇喘18禁视频| 50天的宝宝边吃奶边哭怎么回事| 欧美在线一区亚洲| 麻豆久久精品国产亚洲av| 人妻久久中文字幕网| 九九热线精品视视频播放| 国产精品久久久久久久电影 | 久久久久性生活片| 亚洲精品色激情综合| 成人国产综合亚洲| 免费在线观看黄色视频的| 一进一出好大好爽视频| 日韩国内少妇激情av| 亚洲成人久久爱视频| 免费电影在线观看免费观看| 国产成+人综合+亚洲专区| 成人av在线播放网站| 精品国产美女av久久久久小说| 久久精品亚洲精品国产色婷小说| 一级作爱视频免费观看| 欧美黑人精品巨大| 国产伦一二天堂av在线观看| 成人av在线播放网站| 午夜免费观看网址| 757午夜福利合集在线观看| 午夜激情福利司机影院| 久久久久久久久久黄片| 老司机在亚洲福利影院| 国产区一区二久久| 国产午夜精品久久久久久| 国产精品亚洲一级av第二区| 久久久国产成人精品二区| 国产欧美日韩精品亚洲av| av福利片在线观看| 老司机午夜福利在线观看视频|