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

    水庫對海岸帶海水入侵影響的數(shù)值模擬研究

    2021-05-27 00:02:30季智靈周晨琦趙忠偉
    海洋科學(xué) 2021年4期
    關(guān)鍵詞:潮汐運(yùn)移鹽度

    季智靈, 周晨琦, 趙忠偉, 孔 俊

    水庫對海岸帶海水入侵影響的數(shù)值模擬研究

    季智靈, 周晨琦, 趙忠偉, 孔 俊

    (河海大學(xué) 水文水資源與水利工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室, 江蘇 南京 210098)

    修建海岸水庫是滿足濱海地區(qū)淡水資源的保障的一個(gè)重要方式, 為了進(jìn)一步分析海岸水庫的建立對濱海含水層的海水入侵的影響, 采用SUTRA(saturated-unsaturated transport)模型建立具有代表性的二維海岸模型開展定量和定性研究, 結(jié)果表明: 潮汐作用抑制底層鹽水楔的入侵, 加快海陸物質(zhì)交換; 實(shí)例中計(jì)算表明建壩后, 海岸水庫壩體的建立加劇了海水入侵; 水庫水位在平均海平面以下時(shí), 潮汐動力因素占主導(dǎo), 海水入侵嚴(yán)重且隨著水位抬高, 海水入侵程度基本保持不變, 應(yīng)注意海水入侵的防治; 庫區(qū)水位在平均海平面以上時(shí), 海水入侵程度隨著庫區(qū)水位的抬高而減少, 可以調(diào)節(jié)水庫水位進(jìn)行濱海含水層的海水入侵的防控; 海岸水庫在建設(shè)中存在一個(gè)控制水位, 當(dāng)水庫水位在控制水位以上, 庫區(qū)水體鹽度和含水層淡水儲量能得到保證; 庫區(qū)高水位時(shí)內(nèi)陸物質(zhì)運(yùn)移時(shí)間較長, 運(yùn)移區(qū)域較廣, 氮類污染物易在庫底含水層堆積, 需重視含水層的地下水污染防治。

    海水入侵; 鹽水楔; 海岸水庫; 物質(zhì)輸運(yùn)

    近年來, 在濱海地區(qū), 經(jīng)濟(jì)的快速發(fā)展加劇了水資源短缺的問題, 中國、韓國和新加坡等國家修建了大量海岸水庫工程來提供淡水資源[1]。而在濱海地區(qū), 海水入侵問題廣泛存在, 并且成為海內(nèi)外學(xué)者研究的重點(diǎn)。海水入侵會導(dǎo)致土壤鹽漬化, 地下水水質(zhì)惡化等問題, 這都是制約濱海地區(qū)經(jīng)濟(jì)發(fā)展的重要因素[2-3], 海水入侵的防治日趨關(guān)鍵, 常見的防治工程有淡水帷幕工程和地下壩工程[4]。海岸水庫的建立改變了海岸水文地質(zhì)環(huán)境, 對海水入侵的程度有著一定的影響[5]。為此深入研究海岸水庫建設(shè)下, 海水入侵的影響方式和程度, 并進(jìn)行海水入侵的防治尤為必要。

    海水入侵的研究聚焦于咸淡水界面的形狀、運(yùn)移機(jī)理和規(guī)律, 研究海水入侵的方法之一是建立數(shù)值模型來討論, 一種是突變界面模型[6-8], 另一種是過渡帶模型[9-11]。Robinson等[12]采用有限元法模擬了地下河口地下水流動和鹽分運(yùn)移過程, 武雅潔等[13]使用OpenGeoSys軟件模擬了潮汐波動對海水入侵規(guī)律的影響, Nishikawa等[14]使用過渡帶模型中的SUTRA模型研究了洛杉磯多明戈斯地區(qū)的海水入侵, Mahmoodzadeh等[15]利用SUTRA建立了模型, 對海島淡水透鏡體的海水入侵進(jìn)行了研究。本文采取的模型咸淡水過渡帶較寬且由于過渡帶模型能夠反映復(fù)雜水文地質(zhì)條件、人為活動條件等諸多因素下的水流流動和溶質(zhì)遷移, 故采用過渡帶模型中的SUTRA模型[16]。

    21世紀(jì)以來, 許多學(xué)者對海岸水庫與海水入侵的關(guān)系展開了大量的研究, 一方面聚焦于海水入侵對水庫水質(zhì)的影響, Liang等[17]通過同時(shí)求解地表水和地下水方程來模擬瀉湖的鹽度變化, Mao等[18]使用Delft3D在沿海水庫中模擬了海水淡化過程; 另一方面也有許多學(xué)者研究了海岸水庫對海水入侵的影響, Jin等[5]指出當(dāng)干旱時(shí)期水庫水位低時(shí), 海水更有可能侵入水庫, Mehdizadeh等[19]通過物理模型試驗(yàn)和數(shù)值模擬研究了淡水庫的建立對咸淡水界面的影響。對于海岸水庫對海水入侵的研究大部分沒有考慮潮汐作用, 而在潮汐作用下, 海水入侵的程度和方式發(fā)生改變, 會形成上混合區(qū)[20-23], 加快海水與淡水循環(huán), Xin等[24]指出潮汐減弱了含水層—海洋界面上的水通量相位滯后效應(yīng), 加強(qiáng)了鹽通量的相位滯后效應(yīng)。Robinson等[12]指出潮汐形成的上混合區(qū)具有環(huán)流特征, 增加了入海地下水水量。

    本文采用SUTRA模型建立具有代表性的二維海岸剖面模型, 考慮了潮汐作用, 并分不同庫區(qū)水位進(jìn)行數(shù)值模擬研究, 以找出庫區(qū)水位和海水入侵的關(guān)系, 為修建海岸水庫地區(qū)的海水入侵防治提供參考。

    1 數(shù)值模型

    1.1 控制方程

    由美國地質(zhì)調(diào)查局開發(fā)出來的飽和帶/非飽和帶運(yùn)移模型SUTRA(saturated-unsaturated transport)采用了有限元法和有限差分法, 來求解地下水的兩個(gè)物理過程, 一是飽和和非飽和地下水流動, 二是能量或者溶質(zhì)在含水層多孔介質(zhì)中的運(yùn)移, 來解決地下環(huán)境中水流溶質(zhì)和鹽度的飽和和不飽和運(yùn)輸。

    SUTRA模型的二維有限元地下水流動控制方程如式(1)所示:

    SUTRA的溶質(zhì)運(yùn)移方程如式(2)所示:

    式中:為水動力彌散常量,為流體鹽度,*為源流體鹽度。

    根據(jù)Van Genuchten[25]提出的經(jīng)驗(yàn)公式, 非飽和土壤的水力參數(shù)由下式確定:

    式中:和為形狀系數(shù);s為飽和系數(shù);wres為孔隙水殘余飽和度。

    1.2 概念模型

    海岸水庫指修建在海岸地區(qū)的水庫, 其主要作用是利用其堤壩抵擋海水, 并蓄積淡水, 為濱海地區(qū)提供淡水資源。為研究海岸水庫對海水入侵造成的影響, 本文對研究模型進(jìn)行了概化, 構(gòu)建了一個(gè)概化的海灘剖面進(jìn)行研究, 在低潮位處設(shè)立壩體, 壩后蓄淡, 形成海岸水庫。如圖1所示, 建模時(shí)將平均海平面與海岸的交點(diǎn)設(shè)為(150, 20), 在(155, 19)處設(shè)置壩體, 壩高4 m, 壩厚1 m,=30 m,=200 m,=100 m,=10 m。潮汐振幅=1 m, 周期=12 h。

    圖1 概念模型圖

    1.3 網(wǎng)格及參數(shù)設(shè)置

    在SUTRA模型中將計(jì)算模型分剖分為了24 000個(gè)單元, 24 461個(gè)節(jié)點(diǎn), 生成四邊形網(wǎng)格, 網(wǎng)格最小尺度為0.5 m, 在模型中假定為各向同性介質(zhì), 計(jì)算所采用的具體參數(shù)參照了Xin等[24]在模擬中使用的參數(shù), 見表1。

    表1 液體及傳輸介質(zhì)參數(shù)表

    1.4 邊界條件設(shè)置

    如圖1所示, 模型左側(cè)邊界為內(nèi)陸邊界, 設(shè)為定水位邊界, 淡水水位為21 m, 右邊界為海向邊界, 底部邊界設(shè)為不透水邊界, 模型灘面位置的表面邊界考慮三種工況:

    1) 無水庫、無潮汐情況, 我們把表面段設(shè)為定水頭定濃度邊界;

    2) 無水庫, 有潮汐情況,段設(shè)為潮汐動力邊界;

    3) 修建海岸水庫工況, 壩體為不透水邊界,為潮汐動力邊界, 其中段又考慮三種庫區(qū)水位情況: (a) 庫區(qū)不蓄水0 m, (b) 庫區(qū)水位為1 m, (c) 庫區(qū)水位為2 m, 都設(shè)置為定水頭零濃度邊界。

    通過使用SUTRA運(yùn)行計(jì)算模型, 得出每個(gè)節(jié)點(diǎn)的鹽度值, 流速, 繪制鹽度流場圖, 計(jì)算整個(gè)模型區(qū)域的鹽量, 分析計(jì)算結(jié)果, 總結(jié)規(guī)律, 分析不同水位的海岸水庫對海水入侵的影響。

    2 結(jié)果與分析

    2.1 潮汐對海水入侵的影響

    首先對比潮汐作用對濱海地區(qū)含水層的鹽淡水分布情況, 以無潮穩(wěn)定狀態(tài)為初始狀態(tài), 模擬出有潮汐作用下的地下含水層的時(shí)均穩(wěn)態(tài)鹽度分布。模型中內(nèi)陸水位邊界為21 m, 以周期平均的方法繪制在一個(gè)周期內(nèi)含水層的鹽度分布, 并采用拉格朗日質(zhì)點(diǎn)法研究含水層物質(zhì)輸移特征, 在灘面和內(nèi)陸邊界分別釋放一定數(shù)量的粒子, 繪制其運(yùn)動軌跡, 白線上的數(shù)字就是其從釋放點(diǎn)到排出點(diǎn)所經(jīng)歷的時(shí)間, 結(jié)果如圖2、圖3所示, 圖2是無潮汐條件下的穩(wěn)定狀態(tài), 圖3有潮汐條件下的穩(wěn)定狀態(tài)。圖中紅色代表的是海水, 黑線代表50%海水鹽度等值線(17.5)。將咸淡水界面(50%海水鹽度等值線即17.5)與底部邊界的交點(diǎn)稱為海水入侵點(diǎn), 其與海向邊界的水平距離稱為海水入侵距離。

    從圖3可以看出在潮汐條件下, 在潮間帶形成了上部混合區(qū)即USP(upper saline plume), 在潮汐作用下, 高潮期海水在灘面入滲, 并在低潮位線附近流出, 形成了較快的海水循環(huán), USP處的土壤和地下水的含鹽濃度較高, 在密度流作用下, 內(nèi)陸淡水出流通道被向下頂托, 從圖2可以看出無潮汐用下, 咸淡水過渡帶的寬度為2 m; 從圖3可以看出在潮汐作用下, 咸淡水過渡帶的寬度為10 m, 這說明了潮汐讓咸淡水混合作用加大。

    圖2 無潮汐鹽度分布和粒子路徑

    圖3 有潮汐鹽度分布和粒子路徑

    進(jìn)一步從圖3可以看出加上潮汐作用, 鹽水楔會向海側(cè)移動, 海水入侵距離從88.5 m減小到53.0 m, 即鹽水楔后退了35.5 m, 這是因?yàn)樗髟赨SP和鹽水楔之間形成了較窄的淡水排泄通道, 此處的流速較大, 使得鹽分從鹽水楔下部加快稀釋流出到海洋, 同時(shí)形成了向海的水位梯度, 抑制了底層鹽水的入侵, 導(dǎo)致底層鹽水楔后退。通過計(jì)算整個(gè)模型的鹽量變化, 可以發(fā)現(xiàn)在無潮汐作用下的單寬鹽量為6.448×104kg/m, 潮汐作用下含水層的單寬鹽量為4.901×104kg/m, 減少了24.0%, 這也進(jìn)一步說明了潮汐波動效應(yīng)降低了海水入侵的程度。

    觀察粒子運(yùn)移時(shí)間, 可以看出無潮汐作用下, 粒子先水平運(yùn)動到咸淡水交界面, 然后在垂直坡面方向運(yùn)動排出, 越靠近底部含水層的粒子所需時(shí)間越長。對比圖2和圖3可知, 由于潮汐作用形成了上混合區(qū), 并與鹽水楔之間形成了淡水排泄通道, 此處的流速較大, 加快了內(nèi)陸物質(zhì)輸運(yùn), 如從(0, 5)處釋放的粒子運(yùn)移時(shí)間從575.1 d減少到436.8 d, 減少了24.0%。而灘面釋放的粒子運(yùn)移時(shí)間由于鹽水楔的回退而減少。如從(170, 16)釋放的粒子運(yùn)移時(shí)間從6 007.4 d減小到116.7 d, 減小了98%, 明顯減小, 一方面是因?yàn)辂}水楔后退, 縮短了質(zhì)點(diǎn)運(yùn)移距離, 另一方面是因?yàn)槌毕饔么龠M(jìn)了底層咸淡水的混合循環(huán)過程, 導(dǎo)致海側(cè)物質(zhì)在潛水層內(nèi)的輸移時(shí)間變短。

    2.2 潮汐作用下海岸水庫工程對海水入侵的影響

    為了研究海岸水庫工程對海水入侵的影響, 在有潮穩(wěn)定工況的基礎(chǔ)上, 考慮建庫工程。在低潮位點(diǎn)設(shè)立壩體, 將有潮汐時(shí)的穩(wěn)定狀態(tài)設(shè)為初始狀態(tài), 模擬3種水文情況, 分別是壩后蓄水0 m、1 m、2 m, 計(jì)算出建庫后的鹽度分布, 模型模擬的鹽度分布和粒子路徑如圖4—圖6所示。

    首先對比圖3和圖4, 可以看出, 建壩以后, 鹽水楔向內(nèi)陸移動, 海水入侵距離增大, 這是因?yàn)閴误w阻隔了海水與灘面, 潮汐無法作用于潮間帶, 灘面上層鹽水羽逐漸消失, 從穩(wěn)定后的流場看, 之前在鹽水羽和鹽水楔之間形成的較窄的淡水排泄通道消失, 流速變小, 淡水對下層鹽水楔的擠壓作用減弱, 因此鹽水楔前進(jìn)。壩后蓄水0 m工況下的海水入侵距離為80.0 m, 增加了50.9%, 為了進(jìn)一步說明海水入侵程度的變化, 我們計(jì)算建壩后的模型的單寬鹽量, 壩后蓄水0 m時(shí), 潛水層內(nèi)的單寬鹽量為6.392×104kg/m, 增加了30.4%, 此時(shí)海水入侵程度較大。

    圖4 加壩蓄水0 m鹽度分布和粒子路徑

    圖5 加壩蓄水1 m鹽度分布和粒子路徑

    圖6 加壩蓄水2 m鹽度分布和粒子路徑

    進(jìn)一步抬高庫區(qū)水位至1 m, 從圖5看出海水入侵距離基本沒變, 單寬鹽量為5.871×104kg/m, 只減少了8.2%, 過渡帶變寬; 抬高水位至2 m, 從圖6可以看出海水入侵距離迅速減少到35.6 m, 比庫區(qū)水位0 m工況下減少了59.1%, 計(jì)算模型單寬鹽量為2.731×104kg/m, 比庫區(qū)水位0 m工況下減少了57.3%。

    為了進(jìn)一步研究不同水位工況下對海水入侵的影響, 設(shè)置6種庫區(qū)水位工況, 分別為0 m、0.5 m、1 m、1.5 m、2 m、2.5 m, 計(jì)算出對應(yīng)的海水入侵距離和模型單寬鹽量, 并擬合相關(guān)曲線, 如圖7, 圖8所示, 圖中×點(diǎn)為數(shù)據(jù)標(biāo)記??梢钥闯鰩靺^(qū)水位小于平均海平面時(shí), 即庫區(qū)水位小于1 m, 海水入侵距離和模型單寬鹽量變化較小, 這是因?yàn)樵摴r下, 含水層的水力梯度由潮汐波動控制, 潮汐動力因素占主導(dǎo), 而出現(xiàn)過渡帶變寬現(xiàn)象源于庫區(qū)存在水力梯度, 加強(qiáng)了過渡帶間的彌散作用; 當(dāng)庫區(qū)水位高于平均海平面, 海水入侵距離和含鹽量迅速減少, 水庫水位為2.5 m時(shí), 海水入侵距離減少了68.8%, 單寬鹽量減少了66.9%, 顯著減少, 這是因?yàn)榇藭r(shí)庫區(qū)水位高于平均海平面, 在含水層形成了新的水力梯度, 隨著水位的抬高導(dǎo)致上層淡水對下層鹽水楔的擠壓作用加強(qiáng), 使得鹽水楔向海一側(cè)推移, 抑制了海水入侵程度, 海水入侵距離降低, 含水層鹽量降低。

    圖7 庫區(qū)水位和海水入侵距離關(guān)系

    圖8 庫區(qū)水位和單寬鹽量關(guān)系

    觀察粒子運(yùn)移時(shí)間, 庫區(qū)水位在平均海平面以下時(shí), 隨著水位的抬高, 過渡帶的彌散作用加強(qiáng), 海側(cè)釋放粒子運(yùn)移路徑延長, 運(yùn)移時(shí)間增加, 如庫區(qū)水位為1.0 m時(shí), 在(190, 12)釋放的粒子運(yùn)移時(shí)間比水位為0 m時(shí)增加了66.9%。庫區(qū)水位在平均海平面以上時(shí), 隨著水庫水位的抬高, 內(nèi)陸粒子運(yùn)移時(shí)間增加, 海側(cè)邊界粒子運(yùn)移時(shí)間減少, 這是因?yàn)閹靺^(qū)高水位使得鹽水楔向海一側(cè)推移, 內(nèi)陸粒子運(yùn)移的路徑變長, 海側(cè)邊界粒子運(yùn)移的路徑相對變短。不難得出庫區(qū)高水位會導(dǎo)致陸源污染物運(yùn)移路徑變長, 運(yùn)移時(shí)間變長, 運(yùn)移區(qū)域變廣, 高水位抑制了鹽水楔的入侵, 過渡帶后退變窄, 過渡帶是硝化和反硝化的重要區(qū)域, 影響到陸源氮類污染物的硝化及隨后的反硝化過程, 庫區(qū)底部氮類污染物大量堆積, 容易對水庫水質(zhì)產(chǎn)生影響。

    從圖4可以看出, 在庫區(qū)水位較低時(shí), 庫區(qū)底部含水層鹽度較高, 此時(shí)庫區(qū)水體易受海水入侵的鹽分滲入, 為了進(jìn)一步說明水位和海水入侵的關(guān)系, 設(shè)置了庫區(qū)6種水位, 分別是0 m、0.5 m、1 m、1.5 m、2 m、2.5 m, 算出不同水位工況下模型庫底F點(diǎn)(155, 19)處鹽度值, 擬合出曲線, 繪制出水位與庫區(qū)底部鹽度關(guān)系曲線圖, 如圖9所示, 圖中×點(diǎn)為數(shù)據(jù)標(biāo)記。圖9說明, 水庫水位在平均海平面以上時(shí), 底鹽度變化較小, 都維持在相對較高的鹽度; 當(dāng)水庫水位在平均海平面以上時(shí), 隨著水位的抬高, 庫區(qū)底部邊界的鹽度迅速降低, 在庫區(qū)水位2 m時(shí)庫區(qū)底部鹽度接近0, 此后水位抬高, 庫區(qū)底部鹽度保持為0, 所以要保證庫區(qū)水體鹽度, 讓水庫不受鹽分滲入, 該模型的庫區(qū)控制水位應(yīng)在2 m以上。

    圖9 庫區(qū)水位和庫底鹽度關(guān)系

    從圖6可以觀察到海岸水庫在高水位時(shí)地下含水層淡水儲量較多, 而保證地下含水層淡水儲量是修建海岸水庫需要考慮的重要因素。一般來說將水體鹽度在0.5以下的當(dāng)作淡水, 為了分析庫區(qū)水位與地下淡水儲量關(guān)系, 設(shè)置了6個(gè)水位, 分別是0 m、0.5 m、1 m、1.5 m、2 m、2.5 m, 計(jì)算出各個(gè)水位工況下的地下含水層單寬淡水儲量, 擬合出曲線, 如圖10所示, 圖中×點(diǎn)為數(shù)據(jù)標(biāo)記??梢钥闯? 當(dāng)水庫水位在平均海平面以下時(shí), 隨著水位的抬高, 含水層單寬淡水儲量變化不大; 當(dāng)水庫水位在平均海平面以上時(shí), 隨著水庫水位的抬高, 單寬淡水儲量迅速增加, 水庫水位為2 m時(shí)含水層單寬淡水儲量為3 953 m2, 水庫水位增加到2.5 m, 淡水儲量只增加了13 m2。因此在考慮淡水儲量盡可能保持充足的情況下, 水庫控制水位應(yīng)在2 m以上。

    圖10 庫區(qū)水位和淡水儲量關(guān)系

    3 結(jié)論

    本文采用SUTRA(saturated-unsaturated transport)模型建立具有代表性的二維海岸模型, 研究潮汐作用下海岸水庫工程對海水入侵的影響, 結(jié)果表明:

    1) 潮汐作用抑制底層鹽水楔的入侵, 加快海陸物質(zhì)交換。

    2)海岸水庫壩體的建立加劇了海水入侵, 水庫水位在平均海平面以下時(shí), 含水層潮汐動力因素占主導(dǎo), 海水入侵程度基本保持不變, 海岸水庫在旱季時(shí)存在水位較低的情況, 應(yīng)注意海水入侵的防治, 以及對庫區(qū)水質(zhì)的監(jiān)測。

    3)庫區(qū)水位在平均海平面以上時(shí), 海水入侵程度隨著庫區(qū)水位的抬高而減少, 可以人為地調(diào)控水庫水位來進(jìn)行地下水的退咸, 從而進(jìn)行濱海地區(qū)海水入侵防治。

    4)通過研究水位與庫底鹽度、含水層淡水儲量關(guān)系, 得出海岸水庫在建設(shè)中存在一個(gè)控制水位, 水庫水位在控制水位以上, 庫區(qū)水體鹽度和含水層淡水儲量能得到保證。

    5)當(dāng)庫區(qū)水位較高時(shí), 內(nèi)陸物質(zhì)運(yùn)移時(shí)間較長, 運(yùn)移區(qū)域更廣, 需重視相關(guān)的濱海污染物的入海防治。庫區(qū)高水位削弱庫區(qū)底部陸源氮類污染物在鹽淡水混合區(qū)的硝化反硝化過程, 影響庫區(qū)底部含水層的水質(zhì)。

    [1] 袁文秀, 楊樹灘, 張明. 江蘇省如東縣沿海平原水庫水源論證[J]. 人民長江, 2007, 38(6): 35-37. YUAN Wenxiu, YANG Shutan, ZHANG Ming. Arguments of the coastal reservoir in Ru Dong JiangSu Province Yangtze River[J]. Yangtze River, 2007, 38(6): 35-37.

    [2] 李冬順, 楊勁松, 姚榮江. 生態(tài)風(fēng)險(xiǎn)分析用于蘇北灘涂土壤鹽漬化風(fēng)險(xiǎn)評估研究[J]. 土壤學(xué)報(bào), 2010, 47(5): 51-58.LI Dongshun, YANG Jinsong, YAO Rongjiang. Application of ecological risk analysis to soil salinization risk assessment of coastal tidal flat in north Jiangsu Provin-ce[J]. Acta Pedologica Sinica, 2010, 47(5): 51-58.

    [3] 姚榮江, 楊勁松, 陳小兵, 等. 蘇北海灘典型圍墾區(qū)土壤鹽漬化風(fēng)險(xiǎn)評估研究[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào), 2010, 18(5): 1000-1006. YAO Rongjiang, YANG Jinsong, CHEN Xiaobing, et al. Evaluating soil salinization risk in typical coastal reclaimed regions in North Jiangsu Province. Chinese Journal of Eco-Agriculture, 2010, 18(5): 1000-1006.

    [4] 王輝, 許學(xué)工. 淡水帷幕防治海水入侵的國際經(jīng)驗(yàn)[J]. 中國農(nóng)村水利水電, 2017(4): 62-66. WANG Hui, XU Xuegong. International experience on hydraulic barrier serving as a seawater intrusion coun-termeasure[J]. China Rural Water and Hydropower, 2017(4): 62-66.

    [5] JIN G Q, MO Y M, LI M D, et al. Desalinization and salinization: A review of major challenges for coastal reservoirs[J]. Journal of Coastal Research, 2019, 35(3): 664-672.

    [6] JAKOVOVIC D, WERNER A D, DE LOUW P G B, et al. Saltwater upconing zone of influence[J]. Advances in Water Resources, 2016, 94: 75-86.

    [7] KARATZAS G P, DOKOU Z. Optimal management of saltwater intrusion in the coastal aquifer of Malia, Crete (Greece), using particle swarm optimization[J]. Hydrogeology Journal, 2015, 23(6): 1181-1194.

    [8] PARK N, SHI L. A comprehensive sharp-interface simulation-optimization model for fresh and saline groundwater management in coastal areas[J]. Hydrogeology Journal, 2015, 23(6): 1195-1204.

    [9] COBANER M, YURTAL R, DOGAN A, et al. Three- di-mensional simulation of seawater intrusion in coastal aquifers: A case study in the Goksu Deltaic Plain[J]. Journal of Hydrology, 2012, 464-465: 262-280.

    [10] XIAO H, WANG D, HAGEN S C, et al. Assessing the impacts of sea-level rise and precipitation change on the surficial aquifer in the low-lying coastal alluvial plains and barrier islands, east-central Florida (USA)[J]. Hydrogeology Journal, 2016, 24(7): 1-16.

    [11] 薛顯武, 陳喜, 魏玲娜, 等. 基于GIS的海水入侵可視化數(shù)值模擬系統(tǒng)[J]. 地下水, 2006, 29(5): 35-39. XUE Xianwu, CHEN Xi, WEI Lingna, et al. Visualized numeric simulation system of seawater intrusion based on GIS[J]. Ground Water, 2006, 29(5): 35-39.

    [12] ROBINSON C, LI L, BARRY D A. Effect of tidal forcing on a subterranean estuary[J]. Advances in Water Resources, 2007, 30(4): 851-865.

    [13] 武雅潔, 楊自良, 程從敏, 等. 潮汐波動對潛水含水層海水入侵規(guī)律的影響研究[J]. 中國海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2020, 50(10): 91-98.WU Yajie, YANG Ziliang, CHENG Congmin. The influence of tidal oscillations on the behaviors of seawater intrusion in unconfined coastal aquifers[J]. Periodical of Ocean University of China, 2020, 50(10): 91-98.

    [14] NISHIKAWA T, SIADE A J, REICHARD E G, et al. Stratigraphic controls on seawater intrusion and implications for groundwater management, Dominguez Gap area of Los Angeles, California, USA[J]. Hydrogeology Journal, 2009, 17(7): 1699-1725.

    [15] MAHMOODZADEH D, KETABCHI H, ATAIE-ASHTIANI B, et al. Conceptualization of a fresh groundwater lens influenced by climate change: A modeling study of an arid- region island in the Persian Gulf, Iran[J]. Journal of Hydrology, 2014, 519: 399-413.

    [16] VOSS C I, PROVOST A M. SUTRA: A model for 2D or 3D saturated-unsaturated, variable-density ground-water flow with solute or energy transport[R]. U.S. Geological Survey Water-Resources Investigations Report 2002-4231, 2002. https://doi.org/10.3133/wri024231

    [17] LIANG D, FALCONER R A, LIN B. Coupling surface and subsurface flows in a depth averaged flood wave model[J]. Journal of Hydrology, 2007, 337(1-2): 147-158.

    [18] MAO X Z, CHEN F Y, YU Q W, et al. Numerical prediction of water desalination in polder reservoir[J]. Journal of Hydraulic Engineering, 2004.

    [19] MEHDIZADEH S S, KETABCHI H, GHOROQI M, et al. Experimental and numerical assessment of saltwater recession in coastal aquifers by constructing check dams[J]. Journal of Contaminant Hydrology, 2020, 231: 103637.

    [20] KUAN W K, JIN G, XIN P, et al. Tidal influence on seawater intrusion in unconfined coastal aquifers[J]. Water Resources Research, 2012, 48(2): 136-149.

    [21] MAO X, ENOT P, BARRY D A, et al. Tidal influence on behaviour of a coastal aquifer adjacent to a low-relief estuary[J]. Journal of Hydrology, 2006, 327(1): 110-127.

    [22] WERNER A D, LOCKINGTON D A. Tidal impacts on riparian salinities near estuaries[J]. Journal of Hydrology, 2006, 328(3/4): 511-522.

    [23] XIN P, ROBINSON C, LI L, et al. Effects of wave forcing on a subterranean estuary[J]. Water Resources Research, 2010, 46(12): W12505.

    [24] XIN P, WANG S S J, LU C, et al. Nonlinear interactions of waves and tides in a subterranean estuary[J]. Geophysical Research Letters, 2015, 42: 2277-2284.

    [25] VAN GENUCHTEN, M TH. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils[J]. Soil Science Society of America Journal, 1980, 44(5): 892-898.

    Numerical study of the influence of coastal reservoir on seawater intrusion

    JI Zhi-ling, ZHOU Chen-qi, ZHAO Zhong-wei, KONG Jun

    (State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering, Hohai University, Nanjing 210098, China)

    The construction of coastal reservoirs is one of alternative methods which plays an important role in guaranteeing freshwater resource supply in coastal areas. To analyze the impact of coastal reservoirs on the seawater intrusion in coastal aquifers, a SUTRA (saturated-unsaturated transport) model was adopted to establish a representative two-dimensional coastal case for the study. Quantitative and qualitative research has been carried out based on various scenarios. Results reveal that the tidal action can inhibit the intrusion of the bottom saltwater wedgeand accelerate the exchange between the subsurface and surface water before the construction of the coastal reservoir. After the establishment of the coastal reservoir, such features can be changed corresponding to different reservoir water levels. When the reservoir water level is below the mean sea level, tidal dynamic factors play a dominant role and lead to a more serious seawater intrusion. Such trends remain unchanged when the water level further rises. Muchattention should be given to the prevention and control of seawater intrusion under this condition. Moreover, when the water level of the reservoir area is above the mean sea level, the degree of seawater intrusion decreases with the elevation of the water level in the reservoir. A series of numerical results have confirmed that there is a critical water level above which the water level of the reservoir can effectively ensure the water quality and freshwater resource supply in the reservoir. Further study has also shown that when the water level in the reservoir area is high, the inland material transport time is lengthened, the transport area is widened, and the nitrogen pollutants are easy to accumulate in the reservoir’s bottom aquifer. Such potential risks deserve much attention in designing a coastal reservoir.

    seawater intrusion; saline wedge; coastal reservoir; mass transport

    Oct. 13, 2020

    P731.2

    A

    1000-3096(2021)04-0075-08

    10.11759/hykx20201013002

    2020-10-13;

    2021-01-15

    國家重點(diǎn)研發(fā)項(xiàng)目(2019YFC0409004)

    [National Key Research and Development Project, No. 2019YFC0409004]

    季智靈(1996—), 男, 江蘇南通人, 在讀研究生, 主要研究方向: 水動力研究等, E-mail: 1016074588@qq.com; 孔俊(1978—),通信作者, 男, 江蘇鎮(zhèn)江人, 教授, 主要研究方向: 河口海岸水動力學(xué)研究等, E-mail: kongjun999@126.com

    (本文編輯: 趙衛(wèi)紅)

    猜你喜歡
    潮汐運(yùn)移鹽度
    潮汐與戰(zhàn)爭(上)
    曲流河復(fù)合點(diǎn)壩砂體構(gòu)型表征及流體運(yùn)移機(jī)理
    東營凹陷北帶中淺層油氣運(yùn)移通道組合類型及成藏作用
    絕美海灘
    潮汐式灌溉控制系統(tǒng)的設(shè)計(jì)及應(yīng)用
    電子制作(2017年9期)2017-04-17 03:00:56
    干法紙的潮汐
    生活用紙(2016年6期)2017-01-19 07:36:25
    鹽度和pH對細(xì)角螺耗氧率和排氨率的影響
    開采過程中上覆急傾斜巖層運(yùn)移規(guī)律模擬與研究
    鹽度脅迫對入侵生物福壽螺的急性毒性效應(yīng)
    川西坳陷孝泉-新場地區(qū)陸相天然氣地球化學(xué)及運(yùn)移特征
    青青草视频在线视频观看| 亚洲怡红院男人天堂| 天天一区二区日本电影三级| 欧美成人精品欧美一级黄| 日本av手机在线免费观看| 精品久久久久久电影网 | 日日撸夜夜添| 久久欧美精品欧美久久欧美| 三级经典国产精品| 日本三级黄在线观看| 久99久视频精品免费| 国模一区二区三区四区视频| 九色成人免费人妻av| 美女大奶头视频| 国产熟女欧美一区二区| 亚洲精品色激情综合| 久久亚洲国产成人精品v| 亚洲人成网站在线播| 在线免费十八禁| 一本久久精品| 成人特级av手机在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 91av网一区二区| av天堂中文字幕网| 国产精品99久久久久久久久| 国产成人a区在线观看| 极品教师在线视频| 麻豆成人午夜福利视频| 精品免费久久久久久久清纯| 国产精品久久久久久精品电影| 大话2 男鬼变身卡| 免费一级毛片在线播放高清视频| 三级国产精品片| 1024手机看黄色片| 欧美一区二区亚洲| 久久久久久久久久黄片| 亚洲欧美精品专区久久| 亚洲国产高清在线一区二区三| 波多野结衣巨乳人妻| 美女高潮的动态| 最近的中文字幕免费完整| 日本av手机在线免费观看| 麻豆国产97在线/欧美| 亚洲精品成人久久久久久| 国产午夜福利久久久久久| 国产成人一区二区在线| 国产av一区在线观看免费| 午夜激情欧美在线| 好男人视频免费观看在线| 国产亚洲5aaaaa淫片| 欧美不卡视频在线免费观看| 欧美97在线视频| 国产成人aa在线观看| 久久6这里有精品| 亚洲欧洲日产国产| 97人妻精品一区二区三区麻豆| 一区二区三区高清视频在线| 小说图片视频综合网站| 欧美日韩国产亚洲二区| 免费黄网站久久成人精品| 欧美另类亚洲清纯唯美| 麻豆精品久久久久久蜜桃| 老司机影院成人| 成人漫画全彩无遮挡| 欧美精品国产亚洲| 成人国产麻豆网| 亚洲欧美精品专区久久| 少妇人妻一区二区三区视频| 两个人的视频大全免费| 亚洲第一区二区三区不卡| 美女脱内裤让男人舔精品视频| 欧美日韩国产亚洲二区| 久久久欧美国产精品| 日韩一区二区视频免费看| 亚洲av成人精品一区久久| 级片在线观看| 国产亚洲av片在线观看秒播厂 | 国产高清视频在线观看网站| 九九久久精品国产亚洲av麻豆| 亚洲激情五月婷婷啪啪| 精品熟女少妇av免费看| 国产精品一区二区三区四区免费观看| 女人久久www免费人成看片 | 国产一区二区在线av高清观看| 美女黄网站色视频| 免费播放大片免费观看视频在线观看 | 夜夜爽夜夜爽视频| 亚洲美女搞黄在线观看| 热99re8久久精品国产| 久久久久网色| 亚洲欧美日韩高清专用| 一区二区三区四区激情视频| 欧美激情在线99| 最后的刺客免费高清国语| 97热精品久久久久久| 97热精品久久久久久| 嫩草影院入口| 亚洲四区av| 成人av在线播放网站| 在线免费十八禁| 亚洲欧美日韩高清专用| 韩国高清视频一区二区三区| 我要看日韩黄色一级片| 97人妻精品一区二区三区麻豆| av免费观看日本| 成人三级黄色视频| 五月伊人婷婷丁香| 国产成人免费观看mmmm| 国产男人的电影天堂91| 精品久久国产蜜桃| 亚洲18禁久久av| 91久久精品国产一区二区三区| 亚洲最大成人av| 亚洲av熟女| 丝袜美腿在线中文| 男插女下体视频免费在线播放| 日韩亚洲欧美综合| 1024手机看黄色片| 午夜激情欧美在线| 久久人人爽人人片av| 国产精品久久久久久精品电影小说 | 夫妻性生交免费视频一级片| 国产精品99久久久久久久久| 内射极品少妇av片p| 黄色一级大片看看| 国产老妇女一区| 午夜福利高清视频| 亚洲av中文av极速乱| 日产精品乱码卡一卡2卡三| 国产精品无大码| 国产一区二区在线观看日韩| 午夜a级毛片| 美女国产视频在线观看| 国产国拍精品亚洲av在线观看| 一边亲一边摸免费视频| 视频中文字幕在线观看| 狠狠狠狠99中文字幕| 69av精品久久久久久| 久久午夜福利片| 麻豆成人午夜福利视频| 秋霞在线观看毛片| 韩国高清视频一区二区三区| 直男gayav资源| 国产成人a区在线观看| 亚洲精华国产精华液的使用体验| 国产毛片a区久久久久| h日本视频在线播放| 日本黄色视频三级网站网址| 成人毛片60女人毛片免费| 日韩一本色道免费dvd| 村上凉子中文字幕在线| 美女大奶头视频| 久久婷婷人人爽人人干人人爱| 国产精品人妻久久久影院| a级毛片免费高清观看在线播放| 美女xxoo啪啪120秒动态图| 在线免费十八禁| 免费看av在线观看网站| 人妻系列 视频| 晚上一个人看的免费电影| 如何舔出高潮| 亚洲三级黄色毛片| 一个人看视频在线观看www免费| 小说图片视频综合网站| 国产免费福利视频在线观看| 国产午夜精品久久久久久一区二区三区| 亚洲av中文av极速乱| 成人国产麻豆网| 国产精品三级大全| av在线播放精品| 亚洲欧美中文字幕日韩二区| 啦啦啦啦在线视频资源| 亚洲不卡免费看| 乱码一卡2卡4卡精品| 看黄色毛片网站| 亚洲成人av在线免费| 亚洲av日韩在线播放| 深夜a级毛片| 亚洲自偷自拍三级| 国产亚洲最大av| 一二三四中文在线观看免费高清| 性色avwww在线观看| 中文乱码字字幕精品一区二区三区 | 久久欧美精品欧美久久欧美| 哪个播放器可以免费观看大片| 99久久精品热视频| 亚洲怡红院男人天堂| 九九热线精品视视频播放| 卡戴珊不雅视频在线播放| 亚洲欧美清纯卡通| 国产精华一区二区三区| 国产真实乱freesex| 婷婷色av中文字幕| 天堂网av新在线| 伦理电影大哥的女人| 亚洲精品aⅴ在线观看| 2021天堂中文幕一二区在线观| 热99在线观看视频| 综合色av麻豆| 久久韩国三级中文字幕| 在线免费十八禁| 国产久久久一区二区三区| 日韩欧美在线乱码| 黄色日韩在线| 亚洲成人精品中文字幕电影| 麻豆久久精品国产亚洲av| 看非洲黑人一级黄片| 日本av手机在线免费观看| 麻豆乱淫一区二区| 国产免费又黄又爽又色| 亚洲国产高清在线一区二区三| 亚洲美女视频黄频| 一夜夜www| 极品教师在线视频| 亚洲最大成人中文| 少妇人妻一区二区三区视频| 国产精品女同一区二区软件| 麻豆乱淫一区二区| 亚洲国产日韩欧美精品在线观看| 又粗又爽又猛毛片免费看| 淫秽高清视频在线观看| 国产亚洲一区二区精品| 天天一区二区日本电影三级| 99在线人妻在线中文字幕| 99久国产av精品| 1000部很黄的大片| 午夜老司机福利剧场| 欧美日本视频| 精品久久久久久久久av| 丝袜喷水一区| 日本wwww免费看| 国产中年淑女户外野战色| 在线观看66精品国产| 青春草亚洲视频在线观看| 国产又黄又爽又无遮挡在线| 亚洲成人精品中文字幕电影| 99热6这里只有精品| 欧美区成人在线视频| 听说在线观看完整版免费高清| 久久久a久久爽久久v久久| 精品人妻视频免费看| 男人和女人高潮做爰伦理| 亚洲最大成人中文| 深爱激情五月婷婷| 久久久久久久久中文| 男人舔奶头视频| 少妇裸体淫交视频免费看高清| 国产 一区 欧美 日韩| 美女脱内裤让男人舔精品视频| 精品久久久久久久久av| 高清毛片免费看| 国产精品久久电影中文字幕| 国产精品99久久久久久久久| 亚洲最大成人av| av在线老鸭窝| 久久久久精品久久久久真实原创| 久久精品国产99精品国产亚洲性色| 身体一侧抽搐| 国语对白做爰xxxⅹ性视频网站| 搞女人的毛片| 国产一区二区亚洲精品在线观看| 成人美女网站在线观看视频| 国产精品av视频在线免费观看| 久久人人爽人人片av| 久久精品91蜜桃| av视频在线观看入口| 赤兔流量卡办理| av免费在线看不卡| 精品99又大又爽又粗少妇毛片| 汤姆久久久久久久影院中文字幕 | 亚洲人与动物交配视频| 欧美区成人在线视频| 69人妻影院| 日韩欧美在线乱码| 久久久精品94久久精品| 免费观看精品视频网站| 亚洲18禁久久av| 亚洲av电影不卡..在线观看| 亚洲中文字幕一区二区三区有码在线看| 日日干狠狠操夜夜爽| 成人午夜精彩视频在线观看| 国产亚洲91精品色在线| 亚洲人成网站高清观看| 97超视频在线观看视频| 国产大屁股一区二区在线视频| eeuss影院久久| 精品无人区乱码1区二区| 国内精品美女久久久久久| 国产精品久久视频播放| 一个人看的www免费观看视频| 中国国产av一级| 99热精品在线国产| 亚洲自偷自拍三级| 国产成人a区在线观看| 亚洲怡红院男人天堂| 女人被狂操c到高潮| 欧美一区二区精品小视频在线| 天天躁日日操中文字幕| 国产精品.久久久| 少妇裸体淫交视频免费看高清| 少妇的逼水好多| 国产高清视频在线观看网站| 日本一本二区三区精品| 爱豆传媒免费全集在线观看| 久久精品国产亚洲av天美| 中文欧美无线码| 成人av在线播放网站| 国语自产精品视频在线第100页| 99九九线精品视频在线观看视频| 国产精品不卡视频一区二区| 精品久久国产蜜桃| 国产淫语在线视频| 亚洲av成人av| 国产 一区 欧美 日韩| 久久这里只有精品中国| 别揉我奶头 嗯啊视频| 精品久久久久久久久亚洲| 久久精品久久久久久久性| 亚洲精品自拍成人| 午夜激情欧美在线| 欧美又色又爽又黄视频| 人妻系列 视频| 久久久亚洲精品成人影院| 国产老妇伦熟女老妇高清| 老司机影院毛片| 99国产精品一区二区蜜桃av| 欧美激情国产日韩精品一区| 听说在线观看完整版免费高清| 欧美潮喷喷水| 欧美一区二区国产精品久久精品| av专区在线播放| 天堂√8在线中文| 国产又黄又爽又无遮挡在线| 爱豆传媒免费全集在线观看| 十八禁国产超污无遮挡网站| av又黄又爽大尺度在线免费看 | 国产午夜精品久久久久久一区二区三区| 中文字幕精品亚洲无线码一区| 1000部很黄的大片| 能在线免费观看的黄片| 国产av码专区亚洲av| 欧美成人一区二区免费高清观看| 美女脱内裤让男人舔精品视频| 亚洲婷婷狠狠爱综合网| av福利片在线观看| 毛片一级片免费看久久久久| 一个人看的www免费观看视频| 69av精品久久久久久| 亚洲欧美日韩东京热| 一个人免费在线观看电影| av在线亚洲专区| 久久久久久大精品| 中国美白少妇内射xxxbb| 激情 狠狠 欧美| 国产又色又爽无遮挡免| 国产精品精品国产色婷婷| 在线播放无遮挡| 日韩强制内射视频| 亚洲无线观看免费| 免费观看人在逋| 国产午夜精品论理片| 美女xxoo啪啪120秒动态图| 亚洲国产日韩欧美精品在线观看| 亚洲av成人av| 国产免费又黄又爽又色| 精品久久国产蜜桃| 亚洲成av人片在线播放无| 九九热线精品视视频播放| 久久久欧美国产精品| 国产亚洲最大av| a级一级毛片免费在线观看| 亚洲国产欧美人成| 日韩欧美在线乱码| 美女xxoo啪啪120秒动态图| 国产午夜精品论理片| 欧美一区二区亚洲| 久久鲁丝午夜福利片| 69人妻影院| 看免费成人av毛片| 久久午夜福利片| 小说图片视频综合网站| 亚洲av免费高清在线观看| 国语自产精品视频在线第100页| 91精品伊人久久大香线蕉| 色综合色国产| 能在线免费看毛片的网站| 欧美日韩综合久久久久久| 一区二区三区乱码不卡18| 黄色日韩在线| 国产亚洲精品久久久com| 日本黄大片高清| 亚洲av免费高清在线观看| 我的老师免费观看完整版| 精品一区二区三区视频在线| 日韩精品青青久久久久久| 人妻系列 视频| 黄片无遮挡物在线观看| 成人亚洲精品av一区二区| 草草在线视频免费看| 国产在视频线在精品| 亚洲精品乱久久久久久| 免费一级毛片在线播放高清视频| 亚洲精品国产成人久久av| 久久亚洲国产成人精品v| 少妇裸体淫交视频免费看高清| 天天一区二区日本电影三级| 午夜福利在线观看吧| 九九在线视频观看精品| 国产精品久久久久久精品电影小说 | 日韩av在线大香蕉| 一个人看的www免费观看视频| 丰满少妇做爰视频| 国产一级毛片七仙女欲春2| 成人av在线播放网站| 婷婷色麻豆天堂久久 | 日日撸夜夜添| 亚洲三级黄色毛片| 乱人视频在线观看| 国产精品一区二区三区四区免费观看| 日韩高清综合在线| 欧美一区二区精品小视频在线| 午夜激情福利司机影院| av国产免费在线观看| 亚洲精品乱码久久久v下载方式| 亚洲av男天堂| 国产精品电影一区二区三区| 偷拍熟女少妇极品色| www.av在线官网国产| 国模一区二区三区四区视频| 我的老师免费观看完整版| 51国产日韩欧美| 国产三级在线视频| 蜜桃久久精品国产亚洲av| 色吧在线观看| 国产在线男女| 午夜亚洲福利在线播放| 国产精品乱码一区二三区的特点| 日本av手机在线免费观看| 国语自产精品视频在线第100页| 亚洲欧美清纯卡通| 91午夜精品亚洲一区二区三区| 日本黄色片子视频| 国产伦精品一区二区三区视频9| 国产 一区 欧美 日韩| 床上黄色一级片| 尾随美女入室| 亚洲婷婷狠狠爱综合网| 亚洲精品aⅴ在线观看| 久久精品国产亚洲网站| 日韩精品有码人妻一区| 国产亚洲av片在线观看秒播厂 | 少妇裸体淫交视频免费看高清| 白带黄色成豆腐渣| 亚洲av免费高清在线观看| 亚洲精品乱码久久久久久按摩| 成人国产麻豆网| 午夜福利成人在线免费观看| 亚洲国产精品专区欧美| 青春草视频在线免费观看| 久久久精品欧美日韩精品| 日韩,欧美,国产一区二区三区 | 久久热精品热| 国产单亲对白刺激| 国产成人91sexporn| 国产淫语在线视频| 一级毛片久久久久久久久女| 亚洲成色77777| 久久久色成人| 人人妻人人澡欧美一区二区| 亚洲美女搞黄在线观看| 亚洲av福利一区| 日韩成人伦理影院| 天天躁日日操中文字幕| 身体一侧抽搐| 韩国高清视频一区二区三区| 亚洲怡红院男人天堂| 亚洲三级黄色毛片| 青青草视频在线视频观看| av福利片在线观看| 精品人妻熟女av久视频| 欧美日韩在线观看h| 国产乱人偷精品视频| 麻豆国产97在线/欧美| 搞女人的毛片| 亚洲综合精品二区| 天天躁夜夜躁狠狠久久av| 国产v大片淫在线免费观看| 最近最新中文字幕免费大全7| 麻豆一二三区av精品| 日韩成人av中文字幕在线观看| 国产高清不卡午夜福利| 能在线免费观看的黄片| 日本午夜av视频| 老女人水多毛片| 国产精品爽爽va在线观看网站| 色综合色国产| 久久久国产成人免费| 亚洲欧美日韩无卡精品| 观看免费一级毛片| 听说在线观看完整版免费高清| 黄色日韩在线| 国产黄色视频一区二区在线观看 | 欧美一级a爱片免费观看看| 噜噜噜噜噜久久久久久91| 网址你懂的国产日韩在线| av女优亚洲男人天堂| 日本色播在线视频| 亚洲av成人精品一区久久| 18禁在线播放成人免费| 美女脱内裤让男人舔精品视频| 少妇人妻精品综合一区二区| 亚洲美女视频黄频| 亚洲av中文字字幕乱码综合| 国产免费视频播放在线视频 | 麻豆国产97在线/欧美| 久久鲁丝午夜福利片| 在现免费观看毛片| 成人三级黄色视频| 国产成人精品一,二区| 美女xxoo啪啪120秒动态图| 九九在线视频观看精品| 国产精品国产三级国产专区5o | 看黄色毛片网站| 丝袜美腿在线中文| 久久久久久九九精品二区国产| 中文字幕制服av| АⅤ资源中文在线天堂| 黄色配什么色好看| 简卡轻食公司| 午夜a级毛片| 联通29元200g的流量卡| 女人十人毛片免费观看3o分钟| 永久网站在线| 嘟嘟电影网在线观看| av又黄又爽大尺度在线免费看 | 综合色av麻豆| 色网站视频免费| 国产乱来视频区| 国产私拍福利视频在线观看| 亚洲无线观看免费| 久久国内精品自在自线图片| 深夜a级毛片| 午夜精品一区二区三区免费看| 嫩草影院新地址| 国内精品美女久久久久久| 欧美成人精品欧美一级黄| 高清av免费在线| 伦理电影大哥的女人| 熟女人妻精品中文字幕| 国产伦精品一区二区三区视频9| 小蜜桃在线观看免费完整版高清| 18禁在线无遮挡免费观看视频| 国产三级中文精品| 精品久久久久久久久久久久久| 亚洲av电影在线观看一区二区三区 | 九九爱精品视频在线观看| a级毛片免费高清观看在线播放| 国产高清不卡午夜福利| 国产亚洲一区二区精品| 亚洲国产欧洲综合997久久,| 大香蕉久久网| 免费黄色在线免费观看| 我的女老师完整版在线观看| 国产伦一二天堂av在线观看| 日韩亚洲欧美综合| 午夜a级毛片| 国产成人午夜福利电影在线观看| 18禁裸乳无遮挡免费网站照片| videos熟女内射| 视频中文字幕在线观看| 成人毛片60女人毛片免费| 免费看av在线观看网站| 久久精品影院6| 三级经典国产精品| 青春草国产在线视频| 国产精品人妻久久久影院| 黄色配什么色好看| 午夜激情福利司机影院| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 99热6这里只有精品| 丰满人妻一区二区三区视频av| 特级一级黄色大片| 午夜爱爱视频在线播放| 久久婷婷人人爽人人干人人爱| 日日摸夜夜添夜夜添av毛片| 好男人视频免费观看在线| 国产精品福利在线免费观看| 免费观看精品视频网站| 亚洲高清免费不卡视频| 91aial.com中文字幕在线观看| 欧美一级a爱片免费观看看| 卡戴珊不雅视频在线播放| 国产亚洲5aaaaa淫片| 午夜精品一区二区三区免费看| 99久久无色码亚洲精品果冻| 国产精品国产三级专区第一集| 亚洲在线观看片| av卡一久久| 黄片无遮挡物在线观看| av专区在线播放| 国语对白做爰xxxⅹ性视频网站| 亚洲成人中文字幕在线播放| 赤兔流量卡办理| 国产亚洲精品久久久com| 久久精品久久久久久噜噜老黄 | 毛片一级片免费看久久久久| 人人妻人人澡欧美一区二区| 亚洲三级黄色毛片| 成人无遮挡网站| 欧美日韩精品成人综合77777| 国产成人午夜福利电影在线观看| 99久国产av精品国产电影| 男女下面进入的视频免费午夜| 中文资源天堂在线| av女优亚洲男人天堂|