• <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)移特征
    日日撸夜夜添| 欧美性猛交╳xxx乱大交人| 最近视频中文字幕2019在线8| or卡值多少钱| 99riav亚洲国产免费| 成人鲁丝片一二三区免费| 在线播放国产精品三级| 又紧又爽又黄一区二区| 91麻豆精品激情在线观看国产| 一区二区三区四区激情视频 | 婷婷色综合大香蕉| 免费大片18禁| 男人狂女人下面高潮的视频| 特级一级黄色大片| 国产男靠女视频免费网站| 九色成人免费人妻av| 性插视频无遮挡在线免费观看| 亚洲欧美激情综合另类| 色av中文字幕| 少妇人妻精品综合一区二区 | 在线天堂最新版资源| 成人高潮视频无遮挡免费网站| ponron亚洲| 久久精品国产自在天天线| 97超视频在线观看视频| 有码 亚洲区| 国产国拍精品亚洲av在线观看| 亚洲乱码一区二区免费版| 不卡一级毛片| 日韩欧美国产一区二区入口| 99视频精品全部免费 在线| 日韩人妻高清精品专区| 人人妻人人澡欧美一区二区| 国产91精品成人一区二区三区| 久久精品国产亚洲av香蕉五月| 97人妻精品一区二区三区麻豆| 久久久成人免费电影| 国产精品日韩av在线免费观看| 一个人免费在线观看电影| 一区二区三区四区激情视频 | 欧美黑人巨大hd| 国产v大片淫在线免费观看| 亚洲,欧美,日韩| 日本与韩国留学比较| 变态另类丝袜制服| 国产中年淑女户外野战色| 日日夜夜操网爽| 免费在线观看影片大全网站| 久久人人精品亚洲av| 国产高清视频在线播放一区| 国产淫片久久久久久久久| 日本一二三区视频观看| 禁无遮挡网站| 男人狂女人下面高潮的视频| 午夜久久久久精精品| 国内精品久久久久精免费| 窝窝影院91人妻| 99热这里只有是精品50| 夜夜夜夜夜久久久久| 成年免费大片在线观看| 亚洲一区高清亚洲精品| 国产精品久久电影中文字幕| 国产亚洲欧美98| 国产 一区 欧美 日韩| 非洲黑人性xxxx精品又粗又长| 国产精品伦人一区二区| 99久久精品一区二区三区| 91狼人影院| 中文亚洲av片在线观看爽| 亚洲 国产 在线| 婷婷六月久久综合丁香| 美女高潮喷水抽搐中文字幕| 亚洲内射少妇av| 日韩强制内射视频| or卡值多少钱| 久久久久九九精品影院| 欧美3d第一页| 成人美女网站在线观看视频| 国产免费男女视频| 悠悠久久av| 午夜福利欧美成人| 99在线人妻在线中文字幕| 欧美黑人巨大hd| 国产成人aa在线观看| 欧美高清性xxxxhd video| 天堂av国产一区二区熟女人妻| 成人国产综合亚洲| 成人国产一区最新在线观看| 久久婷婷人人爽人人干人人爱| 亚洲国产欧美人成| 校园人妻丝袜中文字幕| 久久午夜福利片| 99久久精品国产国产毛片| 乱系列少妇在线播放| 91精品国产九色| 日本一本二区三区精品| 天堂网av新在线| 我要搜黄色片| 午夜亚洲福利在线播放| 亚洲国产欧洲综合997久久,| 我要看日韩黄色一级片| 欧美在线一区亚洲| 热99re8久久精品国产| 天堂影院成人在线观看| 欧美+日韩+精品| 91久久精品电影网| av在线天堂中文字幕| 夜夜夜夜夜久久久久| 久久精品国产亚洲av涩爱 | 熟妇人妻久久中文字幕3abv| 国产高清激情床上av| 午夜老司机福利剧场| 99久久成人亚洲精品观看| 色综合色国产| 国产高清视频在线观看网站| 免费在线观看成人毛片| 国产黄片美女视频| 日本精品一区二区三区蜜桃| 在线免费十八禁| 少妇的逼水好多| 午夜福利在线在线| 国内毛片毛片毛片毛片毛片| 亚洲国产精品久久男人天堂| 成人av一区二区三区在线看| 最近中文字幕高清免费大全6 | 国产成年人精品一区二区| 国产精品亚洲美女久久久| 亚洲人成网站在线播放欧美日韩| 国产不卡一卡二| 国产一区二区亚洲精品在线观看| 国产精品爽爽va在线观看网站| 亚洲va在线va天堂va国产| 成人鲁丝片一二三区免费| 精品国内亚洲2022精品成人| 99久久九九国产精品国产免费| 91午夜精品亚洲一区二区三区 | 国产在视频线在精品| 国产老妇女一区| 国产欧美日韩精品亚洲av| 午夜激情欧美在线| 不卡视频在线观看欧美| 成人特级黄色片久久久久久久| 国内精品一区二区在线观看| 欧美黑人巨大hd| 岛国在线免费视频观看| 国产成人aa在线观看| 直男gayav资源| 天堂av国产一区二区熟女人妻| 亚洲精品一区av在线观看| 国产黄色小视频在线观看| 国产精品伦人一区二区| 国内精品美女久久久久久| 嫩草影院精品99| 精品日产1卡2卡| 婷婷亚洲欧美| 少妇的逼好多水| bbb黄色大片| 欧美人与善性xxx| 九九久久精品国产亚洲av麻豆| 春色校园在线视频观看| 日韩在线高清观看一区二区三区 | 啪啪无遮挡十八禁网站| av专区在线播放| 深夜a级毛片| 一区二区三区免费毛片| 成人一区二区视频在线观看| 永久网站在线| 国产一区二区三区在线臀色熟女| 一级av片app| 国产在线男女| 成人三级黄色视频| 给我免费播放毛片高清在线观看| 亚洲无线观看免费| 国产精品一及| 成人鲁丝片一二三区免费| 又紧又爽又黄一区二区| 久久久成人免费电影| 免费电影在线观看免费观看| 中文字幕高清在线视频| 国产黄片美女视频| 黄色配什么色好看| 又粗又爽又猛毛片免费看| 色综合亚洲欧美另类图片| av在线天堂中文字幕| 国产美女午夜福利| 亚洲在线自拍视频| 小说图片视频综合网站| 久久亚洲精品不卡| 国产精华一区二区三区| 亚洲国产精品合色在线| 亚洲精品国产成人久久av| 中出人妻视频一区二区| 成年女人看的毛片在线观看| 欧美性猛交黑人性爽| 国产欧美日韩一区二区精品| 黄色配什么色好看| 我要看日韩黄色一级片| 99热这里只有精品一区| 熟妇人妻久久中文字幕3abv| 国产在线男女| 国产精品不卡视频一区二区| 在线a可以看的网站| 国内少妇人妻偷人精品xxx网站| 高清日韩中文字幕在线| 一区二区三区免费毛片| 别揉我奶头 嗯啊视频| 婷婷亚洲欧美| 国内精品美女久久久久久| 国产不卡一卡二| 成年女人永久免费观看视频| 91久久精品国产一区二区成人| 最新在线观看一区二区三区| 日韩在线高清观看一区二区三区 | 91午夜精品亚洲一区二区三区 | 熟女电影av网| 听说在线观看完整版免费高清| 国产在视频线在精品| 午夜免费成人在线视频| 亚洲欧美清纯卡通| 欧美人与善性xxx| 国产91精品成人一区二区三区| 在线免费观看不下载黄p国产 | www.www免费av| 在线国产一区二区在线| 国产精品免费一区二区三区在线| 精品久久久久久成人av| av在线观看视频网站免费| 春色校园在线视频观看| 窝窝影院91人妻| 精品人妻偷拍中文字幕| 干丝袜人妻中文字幕| 国产视频内射| 国内精品久久久久精免费| 一个人看视频在线观看www免费| av中文乱码字幕在线| 国产高清视频在线观看网站| 日日干狠狠操夜夜爽| 婷婷丁香在线五月| 亚洲,欧美,日韩| 久久久午夜欧美精品| 波多野结衣高清作品| 欧美中文日本在线观看视频| 午夜福利高清视频| 给我免费播放毛片高清在线观看| 少妇熟女aⅴ在线视频| 十八禁网站免费在线| 特大巨黑吊av在线直播| 日本免费一区二区三区高清不卡| 99久久精品国产国产毛片| 少妇熟女aⅴ在线视频| 有码 亚洲区| 欧美国产日韩亚洲一区| 搡老岳熟女国产| 国产淫片久久久久久久久| 老熟妇乱子伦视频在线观看| 乱人视频在线观看| 欧美最黄视频在线播放免费| 日本 av在线| 精品无人区乱码1区二区| 国产视频一区二区在线看| av在线天堂中文字幕| 午夜精品一区二区三区免费看| 成人三级黄色视频| 熟女人妻精品中文字幕| 变态另类成人亚洲欧美熟女| 色播亚洲综合网| 国产在线男女| 变态另类丝袜制服| 成人美女网站在线观看视频| 精品国产三级普通话版| 国产精品伦人一区二区| 精品人妻偷拍中文字幕| 亚洲欧美日韩高清专用| 日韩中文字幕欧美一区二区| 欧美在线一区亚洲| 欧美日本视频| 国产一区二区激情短视频| 欧美国产日韩亚洲一区| 久久精品久久久久久噜噜老黄 | 最近视频中文字幕2019在线8| 在线免费观看的www视频| 日本 欧美在线| 午夜精品一区二区三区免费看| 午夜福利在线在线| 搡老妇女老女人老熟妇| 好男人在线观看高清免费视频| 亚洲天堂国产精品一区在线| 精品福利观看| 久久这里只有精品中国| 91久久精品国产一区二区成人| 国产综合懂色| 男女啪啪激烈高潮av片| 在线观看一区二区三区| 少妇猛男粗大的猛烈进出视频 | 久久亚洲真实| 91在线观看av| 亚洲av成人av| 久久久久久久久久久丰满 | 亚洲va在线va天堂va国产| 中国美白少妇内射xxxbb| 美女高潮喷水抽搐中文字幕| 国产色爽女视频免费观看| 午夜影院日韩av| 人妻少妇偷人精品九色| 亚洲第一区二区三区不卡| 欧美成人a在线观看| 人人妻人人看人人澡| 免费不卡的大黄色大毛片视频在线观看 | 国产伦人伦偷精品视频| 国产亚洲精品综合一区在线观看| 亚洲av美国av| av在线老鸭窝| 禁无遮挡网站| 成人国产综合亚洲| 女人十人毛片免费观看3o分钟| netflix在线观看网站| 哪里可以看免费的av片| 精品一区二区三区av网在线观看| 亚洲一级一片aⅴ在线观看| 亚洲精品乱码久久久v下载方式| 日本与韩国留学比较| 亚洲人与动物交配视频| 亚洲 国产 在线| 日韩亚洲欧美综合| 永久网站在线| 久久午夜亚洲精品久久| 三级国产精品欧美在线观看| 好男人在线观看高清免费视频| 国产老妇女一区| 天天躁日日操中文字幕| 国产精品一区二区三区四区免费观看 | 黄色女人牲交| 亚洲av日韩精品久久久久久密| 给我免费播放毛片高清在线观看| 久久人妻av系列| 国产高清视频在线播放一区| 国产精品精品国产色婷婷| 老司机福利观看| 色哟哟哟哟哟哟| 国产精品精品国产色婷婷| 亚洲精品乱码久久久v下载方式| 日本黄大片高清| 久久6这里有精品| 熟妇人妻久久中文字幕3abv| а√天堂www在线а√下载| 91在线观看av| 欧美一区二区亚洲| 国产精品亚洲美女久久久| or卡值多少钱| 欧美日韩国产亚洲二区| 嫩草影院新地址| 午夜免费男女啪啪视频观看| 久久ye,这里只有精品| 午夜视频国产福利| 国产男女内射视频| 国产伦精品一区二区三区四那| 在线观看三级黄色| 免费少妇av软件| 亚洲自偷自拍三级| 美女中出高潮动态图| 国产成人精品婷婷| 国产免费又黄又爽又色| 高清av免费在线| 国产淫语在线视频| 免费av中文字幕在线| 久久国产精品男人的天堂亚洲 | 青春草视频在线免费观看| 日韩大片免费观看网站| 欧美97在线视频| 一本色道久久久久久精品综合| 亚洲精品一区蜜桃| 亚洲一区二区三区欧美精品| 亚洲色图综合在线观看| 涩涩av久久男人的天堂| 久久久午夜欧美精品| 夜夜骑夜夜射夜夜干| 色网站视频免费| 国产精品久久久久久av不卡| 插逼视频在线观看| 热99国产精品久久久久久7| 久久人妻熟女aⅴ| 中文字幕精品免费在线观看视频 | 亚洲国产欧美人成| 亚洲美女搞黄在线观看| 中文字幕av成人在线电影| 日本欧美国产在线视频| av在线蜜桃| 毛片女人毛片| 国产精品一区二区性色av| 女性被躁到高潮视频| 午夜福利网站1000一区二区三区| 亚洲精品日本国产第一区| 成人亚洲欧美一区二区av| 国产精品免费大片| 妹子高潮喷水视频| 三级经典国产精品| 午夜激情久久久久久久| 极品少妇高潮喷水抽搐| 精品国产乱码久久久久久小说| 日本一二三区视频观看| 欧美精品国产亚洲| 亚洲成人手机| 精品亚洲成国产av| www.色视频.com| 亚洲av综合色区一区| 在线观看av片永久免费下载| 中文字幕精品免费在线观看视频 | 高清黄色对白视频在线免费看 | 一级毛片久久久久久久久女| 日韩视频在线欧美| 我的女老师完整版在线观看| 成人二区视频| 亚洲精品aⅴ在线观看| 高清不卡的av网站| 麻豆成人av视频| 有码 亚洲区| videos熟女内射| 亚洲天堂av无毛| 日本一二三区视频观看| 亚洲国产精品999| www.色视频.com| 男人狂女人下面高潮的视频| 午夜福利网站1000一区二区三区| 毛片一级片免费看久久久久| 另类亚洲欧美激情| 欧美精品国产亚洲| 国产深夜福利视频在线观看| 男人舔奶头视频| 99热全是精品| 国产真实伦视频高清在线观看| 久热久热在线精品观看| 久久久成人免费电影| 高清毛片免费看| 亚洲真实伦在线观看| 男男h啪啪无遮挡| 日韩国内少妇激情av| 欧美精品国产亚洲| 亚洲图色成人| 久久久午夜欧美精品| 日本欧美国产在线视频| 看十八女毛片水多多多| 18禁在线播放成人免费| 久久婷婷青草| 18禁裸乳无遮挡动漫免费视频| 成年av动漫网址| 街头女战士在线观看网站| 久久久久精品性色| 十八禁网站网址无遮挡 | 亚洲成人av在线免费| 亚州av有码| 国产69精品久久久久777片| 高清黄色对白视频在线免费看 | 在线亚洲精品国产二区图片欧美 | 精品少妇黑人巨大在线播放| 91精品伊人久久大香线蕉| 三级经典国产精品| 91久久精品国产一区二区成人| 色5月婷婷丁香| 免费观看无遮挡的男女| 欧美xxⅹ黑人| 人妻系列 视频| 美女主播在线视频| 亚洲国产毛片av蜜桃av| 91久久精品国产一区二区三区| 国产精品99久久久久久久久| 在线观看免费高清a一片| 欧美日韩一区二区视频在线观看视频在线| 国产 精品1| av播播在线观看一区| 久久久久久人妻| 久久久久视频综合| kizo精华| 精品久久久久久久末码| 丝袜脚勾引网站| 国产精品.久久久| 国产成人免费无遮挡视频| 熟女电影av网| 秋霞在线观看毛片| 国产一区二区在线观看日韩| 涩涩av久久男人的天堂| 亚洲精品国产色婷婷电影| 一个人看视频在线观看www免费| 欧美成人一区二区免费高清观看| 边亲边吃奶的免费视频| 51国产日韩欧美| 亚洲精品日韩av片在线观看| 精品人妻一区二区三区麻豆| 免费av中文字幕在线| 日韩,欧美,国产一区二区三区| 一区二区三区乱码不卡18| 国产乱人视频| 中文字幕亚洲精品专区| 国产精品人妻久久久影院| 亚洲欧美成人精品一区二区| 嫩草影院入口| 亚洲av欧美aⅴ国产| 国产精品久久久久久精品古装| www.色视频.com| 久久精品熟女亚洲av麻豆精品| 免费大片18禁| 你懂的网址亚洲精品在线观看| 51国产日韩欧美| 日韩伦理黄色片| 精品久久久噜噜| 99九九线精品视频在线观看视频| .国产精品久久| 国产人妻一区二区三区在| 国产精品无大码| 新久久久久国产一级毛片| 夫妻午夜视频| 少妇高潮的动态图| 亚洲国产色片| 亚洲精品自拍成人| 少妇人妻精品综合一区二区| 亚洲精品久久久久久婷婷小说| 日韩av不卡免费在线播放| 亚洲婷婷狠狠爱综合网| 免费av不卡在线播放| 国产黄片视频在线免费观看| 熟女av电影| 伦理电影免费视频| 日韩一区二区视频免费看| 国产成人精品婷婷| 夜夜爽夜夜爽视频| 国产精品人妻久久久影院| 中文字幕人妻熟人妻熟丝袜美| 久久精品国产自在天天线| 干丝袜人妻中文字幕| 欧美精品一区二区免费开放| 韩国高清视频一区二区三区| 久久久久久伊人网av| 18+在线观看网站| 美女福利国产在线 | 国产视频首页在线观看| 久久女婷五月综合色啪小说| 51国产日韩欧美| 国产精品av视频在线免费观看| av专区在线播放| 美女福利国产在线 | 少妇的逼好多水| 国产在线男女| 美女脱内裤让男人舔精品视频| 啦啦啦中文免费视频观看日本| av在线老鸭窝| 2018国产大陆天天弄谢| 80岁老熟妇乱子伦牲交| 91久久精品国产一区二区三区| 国产欧美日韩精品一区二区| 夫妻午夜视频| 亚洲久久久国产精品| 人人妻人人添人人爽欧美一区卜 | 91午夜精品亚洲一区二区三区| 日韩 亚洲 欧美在线| 日韩一本色道免费dvd| 在线免费观看不下载黄p国产| 汤姆久久久久久久影院中文字幕| 免费在线观看成人毛片| 国产精品成人在线| 免费看日本二区| 视频区图区小说| 直男gayav资源| 七月丁香在线播放| 国产伦理片在线播放av一区| 2021少妇久久久久久久久久久| 欧美日韩在线观看h| 国产色爽女视频免费观看| 日韩成人av中文字幕在线观看| 国产免费福利视频在线观看| 在线亚洲精品国产二区图片欧美 | 国产成人91sexporn| 97超碰精品成人国产| 亚洲精品国产av成人精品| 少妇人妻一区二区三区视频| 中文字幕久久专区| 午夜福利视频精品| 久久久午夜欧美精品| 国产色爽女视频免费观看| 日本av手机在线免费观看| 交换朋友夫妻互换小说| 午夜精品国产一区二区电影| a级毛片免费高清观看在线播放| 欧美最新免费一区二区三区| 性色av一级| 亚洲av福利一区| 岛国毛片在线播放| 狂野欧美白嫩少妇大欣赏| 久久久久久伊人网av| 精品少妇黑人巨大在线播放| 亚洲国产精品成人久久小说| 亚洲怡红院男人天堂| 18+在线观看网站| 大香蕉97超碰在线| 日韩伦理黄色片| 97超碰精品成人国产| 大香蕉97超碰在线| 日韩精品有码人妻一区| 熟女av电影| 在现免费观看毛片| 国产精品av视频在线免费观看| 日日摸夜夜添夜夜添av毛片| 亚洲四区av| 精品国产三级普通话版| 亚洲国产精品999| 女性被躁到高潮视频| 亚洲国产av新网站| 99re6热这里在线精品视频| 国产精品久久久久久精品电影小说 | 黄色视频在线播放观看不卡| 亚洲伊人久久精品综合| 国产探花极品一区二区| 欧美xxxx性猛交bbbb| 日韩大片免费观看网站| 国产精品一区二区三区四区免费观看| 一级黄片播放器| 精品午夜福利在线看|