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

    汾河蘭村—柴村段潛流帶水動(dòng)力特征研究

    2021-07-08 23:27:26季晨龍李海明肖瀚李夢(mèng)雨王瑞芬楊麗霞張雪麗
    人民黃河 2021年5期

    季晨龍 李海明 肖瀚 李夢(mèng)雨 王瑞芬 楊麗霞 張雪麗

    摘?要:采用數(shù)值模擬方法研究汾河蘭村—柴村段的潛流帶水動(dòng)力特征,量化潛水含水層對(duì)河流水位瞬時(shí)抬升1.5 m的響應(yīng)幅度與影響距離,計(jì)算該條件下河流的側(cè)滲補(bǔ)給強(qiáng)度。結(jié)果表明:隨著與河岸線垂直距離的增大,潛水響應(yīng)幅度以及河流的側(cè)滲補(bǔ)給流量、側(cè)滲速度逐漸減小,相同時(shí)刻起始點(diǎn)x=0處的側(cè)滲速度為上游<中游<下游,側(cè)滲補(bǔ)給流量為上游>中游>下游;在5~90 d內(nèi),各斷面潛水響應(yīng)幅度、側(cè)滲速度、側(cè)滲補(bǔ)給流量在與河岸線距離1 000 m之后從峰值逐漸降為0;潛水響應(yīng)幅度隨時(shí)間延長(zhǎng)呈增大趨勢(shì),距河岸線500 m處90 d潛水響應(yīng)幅度為0.18~0.25 m;距河岸線20~1 000 m處的側(cè)滲速度與側(cè)滲補(bǔ)給流量在5~30 d達(dá)到峰值,之后逐漸降低;相同時(shí)刻的影響距離為1—1′上蘭村斷面最遠(yuǎn)、2—2′汾河大橋斷面與3—3′柴村橋斷面相近。

    關(guān)鍵詞:潛流帶;水動(dòng)力特征;側(cè)滲;汾河蘭村—柴村段

    中圖分類號(hào):P641.2;TV138?文獻(xiàn)標(biāo)志碼:A

    doi:10.3969/j.issn.1000-1379.2021.05.016

    Abstract: Numerical modeling method was used to quantify the hydrodynamic characteristics of the subsurface flow in the Lancun-Chaicun section of the Fenhe River such as its response to an instantaneous rise of 1.5 m of Fenhe River water stage and seepage velocity and flow rate from the Fenhe River to the adjacent phreatic aquifer. The results show that the response range of the phreatic aquifer and the lateral seepage rate decreases gradually with an increasing distance to the riparian line. Lateral seepage velocity in x=0 at the same time is upstreammidstream>downstream. Within 5 to 90 d, the response amplitude and lateral seepage flow rate gradually decrease from the peak value to zero after 1 000 m vertical distance from the river riparian line. The response range of phreatic aquifer water level increases with time and the response range of phreatic water reaches 0.18 -0.25 m at 90 days at 500 m away from the river riparian line. The lateral seepage velocity and the flow rate reach the peak value within 5-30 d at the distance of 20-1 000 m from each section to the riparian line of the river and then decrease with time. The influence at the same time is the farthest from section 1-1′ and 2-2′ and 3-3′ are roughly equal.

    Key words: subsurface flow; hydrodynamic characteristics; lateral seepage; Lancun-Chaicun section of Fenhe River

    隨著汾河流域經(jīng)濟(jì)社會(huì)發(fā)展速度加快,地下水位急劇下降,生態(tài)環(huán)境遭到嚴(yán)重破壞。為解決汾河生態(tài)問(wèn)題,山西省通過(guò)調(diào)引客水、污染防治、源頭及干流兩側(cè)污染企業(yè)綜合整治等方式進(jìn)行大規(guī)模治理,流域地下水位得到持續(xù)回升[1],但流域自身的產(chǎn)流功能仍未得到有效恢復(fù)[2]。汾河流域早期曾擁有豐富的濕地資源,在涵養(yǎng)水源、降解污染物、防洪蓄水等方面起到重要作用[3],但后來(lái)大量灘涂、濕地被開(kāi)發(fā)利用,濕地資源不斷減少[4]。為恢復(fù)汾河流域濕地資源,山西省實(shí)施了汾河流域生態(tài)修復(fù)中游核心區(qū)干流蓄水工程,該工程主要通過(guò)河道蓄水、向堤外濕地補(bǔ)水,改善水生態(tài)環(huán)境。該工程蓄水期間使流域濕地水位抬升1~2 m,在濕地補(bǔ)給地下水的過(guò)程中存在諸多問(wèn)題有待深入研究。

    潛流帶是指河流河床內(nèi)水分飽和的沉積物層,是河流地表水與地下水相互作用的交匯區(qū)域,也是河床與河水進(jìn)行物質(zhì)和能量交換的區(qū)域[5]。影響潛流交換的因素有很多種,其中河水流速與流量是主要因素[6-7]。HUNT Bruce[8]研究表明,對(duì)于假設(shè)飽和模式下的河流與含水層系統(tǒng),河水與地下水之間的水量交換強(qiáng)度同河水位與地下水位差成線性關(guān)系;PACKMAN Aaron I.等[9]通過(guò)水槽試驗(yàn)分析得出了潛流交換率與河水流速等因素之間的關(guān)系;HARVEY Judson W.等[10]得出了在河床附近潛流交換隨河流徑流量變化的規(guī)律;LIERMANN Catherine Reidy[11]通過(guò)對(duì)Fauntleroy溪調(diào)查表明,不同流速對(duì)潛流交換勢(shì)的影響不同。對(duì)河床滲透性能時(shí)空變異性及其影響因素的識(shí)別不僅是研究河水與地下水水量交換的關(guān)鍵與難點(diǎn)[12-15],也是河流與含水層相互作用研究的熱點(diǎn)[16-17]。

    目前潛流交換的研究多是在小尺度上進(jìn)行的,例如ZHANG Guotao等[18]研究了潛流帶水交換與河岸線曲率的關(guān)系。如何量化大尺度潛流交換已成為一個(gè)重要課題,侯宏冰等[19]對(duì)黃河溫縣段不同水位變化情況下的側(cè)滲量進(jìn)行了模擬計(jì)算,得出黃河水位的變化會(huì)引起地下水位的相應(yīng)變化;XI Haiyang等[20]在黑河流域下游進(jìn)行的田間試驗(yàn)表明,河流側(cè)滲主要受流量及持續(xù)時(shí)間的控制。對(duì)河流與含水層系統(tǒng)的數(shù)值模擬是區(qū)域尺度上河水與地下水交換研究的重要手段[19,21-24],其中對(duì)河流與含水層系統(tǒng)模型的概化以及對(duì)模型參數(shù)的估算是影響研究結(jié)果的關(guān)鍵[17,25]。

    鑒于目前對(duì)大尺度潛流交換的研究較少,筆者采用數(shù)值模擬的方法,從側(cè)滲速度和側(cè)滲補(bǔ)給流量?jī)煞矫?,研究位于汾河流域生態(tài)修復(fù)中游核心區(qū)干流蓄水工程受水區(qū)的上蘭村—柴村橋段水位瞬時(shí)抬升1.5 m時(shí)的潛流帶水動(dòng)力特征,以期為汾河流域濕地生態(tài)修復(fù)提供參考。

    1?研究區(qū)水文地質(zhì)條件

    研究區(qū)位于太原市西張盆地,最高海拔為850 m,地勢(shì)總體呈北高南低、東西兩側(cè)高中間低的特點(diǎn),地質(zhì)年代為第四系全新統(tǒng)及上更新統(tǒng),地貌類型為沖洪積傾斜平原地貌、沖積平原地貌以及沖洪積交接洼地地貌。研究區(qū)內(nèi)汾河長(zhǎng)度約10 km,屬平原性河流。研究區(qū)地下水主要接收汾河的滲漏補(bǔ)給以及大氣降水入滲補(bǔ)給,地下水位低于河水位,河水同時(shí)向河兩岸含水層側(cè)向滲透補(bǔ)給地下水。降水冬春少而夏秋多,降水量集中在7—9月,多年(1956—2013年)平均降水量為507 mm。研究區(qū)潛水徑流主要受地形、含水層巖性、補(bǔ)給來(lái)源等控制,總體上由北向南、由東西兩側(cè)向中部流動(dòng)。研究區(qū)潛水的排泄方式主要有蒸發(fā)、越流排泄及含水層的側(cè)向徑流3種形式。研究區(qū)地下水為松散巖類孔隙水,含水層多為全新統(tǒng)砂卵礫石層,在0~160 m埋深范圍內(nèi)分為潛水含水層和弱透水層,包氣帶厚度為10~22 m,巖性多為粉土、亞砂土。潛水含水層頂板高程620~824 m,厚度30~150 m,巖性多為砂卵石。弱透水層厚度10~80 m,底板埋深40 ~160 m,巖性多為黏土。

    2?地下水流數(shù)值模型構(gòu)建

    2.1?水文地質(zhì)概念模型的建立

    2.1.1?邊界條件

    根據(jù)地形地貌、地層結(jié)構(gòu)及區(qū)域流場(chǎng)特征,把潛水含水層南邊界與東邊界中下部概化為定流量邊界,東邊界上部、西邊界以及北邊界概化為隔水邊界,弱透水層均概化為隔水邊界(見(jiàn)圖1)。

    2.1.2?含水層空間概化

    本文只關(guān)注與河流水力聯(lián)系密切的潛水含水層,根據(jù)已有的地質(zhì)圖、剖面圖、地質(zhì)報(bào)告等資料,將含水層概化為非均質(zhì)各向同性潛水含水層,地下水流以水平運(yùn)動(dòng)為主。研究區(qū)含水層三維結(jié)構(gòu)與剖面見(jiàn)圖2。

    2.2?地下水流數(shù)值模型的建立

    2.2.1?網(wǎng)格剖分

    對(duì)模擬研究區(qū)進(jìn)行網(wǎng)格剖分,共剖分為20 000個(gè)網(wǎng)格(100行×100列×2層),其中活動(dòng)單元格8 368個(gè),垂向上劃分為兩層(見(jiàn)圖3)。

    2.2.2?模型識(shí)別應(yīng)力期選取

    根據(jù)現(xiàn)有資料,模型以2012年1月—2012年12月為校正期,2013年1月—2013年12月為驗(yàn)證期,以月為應(yīng)力期,每個(gè)應(yīng)力期為一個(gè)時(shí)間步長(zhǎng)。

    2.2.3?源匯項(xiàng)處理

    本文模型中潛水含水層地下水補(bǔ)給項(xiàng)主要包括降水入滲補(bǔ)給、地表水滲漏補(bǔ)給,排泄項(xiàng)主要包括蒸發(fā)、徑流排泄。

    2.3?模型的校正與驗(yàn)證

    將研究區(qū)校正期和驗(yàn)證期潛水流場(chǎng)水位模擬值與實(shí)測(cè)值進(jìn)行比較,發(fā)現(xiàn)兩者的擬合程度較高,說(shuō)明所建立的水文地質(zhì)概念模型和數(shù)值模型基本合理,能較好地模擬研究區(qū)地下水動(dòng)態(tài)變化過(guò)程。以2014年1月地下水位為初始條件,模擬預(yù)測(cè)當(dāng)汾河水位瞬時(shí)抬升1.5 m時(shí)河水側(cè)滲補(bǔ)給地下水的速度和補(bǔ)給量以及地下水位的時(shí)空變化規(guī)律。

    3?結(jié)果與討論

    3.1?不同斷面地下水位變化特征

    為了更好地認(rèn)識(shí)河流補(bǔ)給地下水的水動(dòng)力特征,沿汾河水流方向選取不同斷面,研究其地下水位變化特征、河水補(bǔ)給地下水側(cè)滲速度和補(bǔ)給量的變化特征。在汾河上游(上蘭村)、中游(汾河大橋)、下游(柴村橋)各設(shè)置一個(gè)斷面,分別記為1—1′、2—2′、3—3′,以河流斷面與河岸線的交叉點(diǎn)為原點(diǎn),以垂直河流中心線且遠(yuǎn)離河流的方向?yàn)閤軸。將各斷面的潛水水位與初始水位的差值記為Δh。

    (1)不同時(shí)刻Δh與河岸線距離的關(guān)系曲線見(jiàn)圖4,由圖4可知:Δh隨與河岸線距離的增大總體呈下降趨勢(shì),在x=0處,Δh=1.5 m,隨后Δh值隨x增大在0~20 m急速下降,降幅占河流水位瞬時(shí)抬升高度的78.7%~97.3%,在20 m以后緩慢下降直至在距河岸線足夠遠(yuǎn)處為0。①1—1′上蘭村斷面河水位抬升5 d時(shí)Δh從河岸線處的1.5 m降為20 m處的0.04 m,降幅為1.46 m,占河流水位瞬時(shí)抬升高度的97.3%,之后從20 m到3 000 m的距離內(nèi),從0.04 m降至0 m;90 d時(shí)Δh在距河岸線20 m處從1.5 m降至0.32 m,降幅為1.18 m,占河流水位瞬時(shí)抬升高度的78.7%,之后從20 m到3 000 m的距離內(nèi),從0.32 m降至0.03 m,降幅為0.29 m。②2—2′汾河大橋斷面河水位抬升5 d時(shí)Δh在距河岸線20 m處降幅占河流水位瞬時(shí)抬升高度的96.7%,之后從20 m到3 000 m的距離內(nèi),降至0 m;90 d時(shí)Δh在距河岸線20 m處降幅占河流水位瞬時(shí)抬升高度的81.3%,之后從20 m到3 000 m的距離內(nèi)降至0 m。③3—3′柴村橋斷面河水位抬升5 d時(shí)Δh在距河岸線20 m處降幅占河流水位瞬時(shí)抬升高度的96.7%,之后從20 m到3 000 m的距離內(nèi)降至0 m;90 d時(shí)Δh在距河岸線20 m處降幅占河流水位瞬時(shí)抬升高度的80%,之后從20 m到3 000 m的距離內(nèi)降至0 m。

    (2)為了更直觀地觀察不同斷面Δh隨距離的變化規(guī)律,選取各斷面x=500 m處繪制Δh隨時(shí)間的變化曲線(見(jiàn)圖5),發(fā)現(xiàn)各斷面Δh隨時(shí)間推移均呈增大趨勢(shì),并且各斷面在同一時(shí)刻的Δh值1—1′上蘭村斷面最大,2—2′汾河大橋斷面與3—3′柴村橋斷面大致接近,說(shuō)明受河流水位抬升影響后的地下水位在上游漲幅較大。1—1′上蘭村斷面距河岸線500 m處經(jīng)過(guò)90 d后Δh增幅最大,為0.25 m。2—2′汾河大橋斷面與3—3′柴村橋斷面Δh增幅大致接近,2—2′汾河大橋斷面Δh在90 d內(nèi)由0增加到0.18 m,3—3′柴村橋斷面在90 d內(nèi)由0增加到0.19 m。故距河岸線500 m處90 d潛水響應(yīng)幅度范圍為0.18~0.25 m。

    (3)定義Δh=0.01 m所對(duì)應(yīng)的垂直于河岸線的距離為河流水位瞬時(shí)抬升1.5 m時(shí)對(duì)潛水水位的影響距離。根據(jù)不同斷面的影響距離(見(jiàn)表1)可得:相同時(shí)間下,影響距離1—1′上蘭村斷面最遠(yuǎn),2—2′汾河大橋斷面與3—3′柴村橋斷面相近。1—1′上蘭村斷面5、10、30 d影響距離分別為458、851、2 150 m,60 d時(shí)超出研究區(qū)范圍;2—2′汾河大橋斷面5、10、30、60、90 d影響距離分別為349、657、1 438、2 041、2 538 m;3—3′柴村橋斷面5、10、30、60、90 d影響距離分別為364、642、1 360、2 030、2 517 m。

    3.2?側(cè)滲速度的變化規(guī)律

    根據(jù)Darcy定律,河流的側(cè)滲速度v可用如下公式來(lái)表示[26]:

    式中:K為滲透系數(shù);Δh0,t為在距河岸線0 m的位置經(jīng)過(guò)t時(shí)間水位上升的幅度;a為壓力傳導(dǎo)系數(shù),a=Khm/μ(hm為含水層厚度,μ為給水度);x為距離河岸線距離。

    滲透系數(shù)K與給水度μ根據(jù)數(shù)值模擬結(jié)果輸入,含水層厚度hm根據(jù)數(shù)值模型的潛水水位及潛水含水層底板高程插值獲得。

    (1)同一時(shí)刻不同位置的側(cè)滲速度不同,起始點(diǎn)x=0處側(cè)滲速度最大:

    隨著x增大,遠(yuǎn)離河流的側(cè)滲速度逐漸變小。

    (2)河流水位瞬時(shí)抬升1.5 m時(shí)側(cè)滲速度的時(shí)空變化規(guī)律見(jiàn)圖6,可以看出:相同時(shí)刻起始點(diǎn)x=0處的側(cè)滲速度上游(1—1′上蘭村斷面)<中游(2—2′汾河大橋斷面)<下游(3—3′柴村橋斷面)。①1—1′上蘭村斷面(見(jiàn)圖6(a)),起始點(diǎn)x=0處的側(cè)滲速度隨時(shí)間的推移逐漸減小并且下降速率逐漸減小。該斷面5 d時(shí)的側(cè)滲速度從x=0處的0.11 m/d開(kāi)始快速下降,隨著距離的增加,在1 750 m處降為0。90 d時(shí)的側(cè)滲速度從x=0處的0.03 m/d開(kāi)始緩慢下降,結(jié)合其他時(shí)間曲線變化趨勢(shì),可以推測(cè)該斷面90 d時(shí)的側(cè)滲速度在距河岸線足夠遠(yuǎn)處降為0。②2—2′汾河大橋斷面(見(jiàn)圖6(b)),5 d時(shí)的側(cè)滲速度從x=0處的0.13 m/d開(kāi)始快速下降,隨著距離的增加,在1 596 m處降為0。90 d時(shí)的側(cè)滲速度從x=0處的0.03 m/d開(kāi)始緩慢下降,隨著距離的增加,在3 000 m處降為0。③3—3′柴村橋斷面(見(jiàn)圖6(c)),5 d時(shí)的側(cè)滲速度從x=0處的0.15 m/d開(kāi)始快速下降,隨著距離的增加,在1 512 m處降為0。90 d時(shí)的側(cè)滲速度從x=0處的0.03 m/d開(kāi)始緩慢下降,隨著距離的增加,在3 000 m處降為0.003 m/d。

    (3)距河岸線不同距離側(cè)滲速度隨時(shí)間的變化規(guī)律見(jiàn)圖7。同一位置(x固定)側(cè)滲速度先快速增大,而后隨時(shí)間的延長(zhǎng)而逐漸變小,直至為0。①1—1′上蘭村斷面(見(jiàn)圖7(a)),距河岸線20 m與100 m處的側(cè)滲速度先從0開(kāi)始快速增大,至6 d達(dá)到峰值,分別為0.095、0.094 m/d,隨后再緩慢減小。距河岸線500 m處的側(cè)滲速度從0開(kāi)始快速增大,至8 d達(dá)到峰值0.06 m/d;距河岸線1 000 m處的側(cè)滲速度從0開(kāi)始快速增大,至26 d達(dá)到峰值0.03 m/d;距河岸線2 000、3 000 m處的側(cè)滲速度從0開(kāi)始緩慢增大,至90 d時(shí)尚未達(dá)到峰值。②2—2′汾河大橋斷面(見(jiàn)圖7(b)),與河岸線相距20、100、500、1 000 m處的側(cè)滲速度在90 d內(nèi)首先從0開(kāi)始增大,達(dá)到峰值后再緩慢減小,峰值分別為0.11、0.10、0.06、0.03 m/d;與河岸線相距2 000、3 000 m處的側(cè)滲速度在90 d內(nèi)未達(dá)到峰值。③3—3′柴村橋斷面(見(jiàn)圖7(c)),側(cè)滲速度隨時(shí)間的變化規(guī)律與1—1′上蘭村斷面、2—2′汾河大橋斷面大致相同,與河岸線相距20、100、500、1 000 m處的側(cè)滲速度峰值分別為0.13、0.12、0.06、0.03 m/d。

    3.3?側(cè)滲補(bǔ)給流量的變化規(guī)律

    水位變幅不大時(shí),潛水含水層飽和厚度可以用其平均厚度hm表示,通過(guò)任一斷面處的單寬流量q[26]為

    根據(jù)側(cè)滲補(bǔ)給流量的時(shí)空變化規(guī)律可得:

    (1)同一時(shí)刻不同位置的單寬流量不等。當(dāng)x=0時(shí),單寬流量最大,其值即河流向一側(cè)補(bǔ)給地下水的單寬流量:

    隨著遠(yuǎn)離河岸線(x增大)單寬流量逐漸變小。隨著時(shí)間的延長(zhǎng),單寬流量的減小速率逐漸變小。同一位置的單寬流量隨著時(shí)間的延長(zhǎng)先增大然后逐漸減小。

    (2)不同時(shí)刻側(cè)滲補(bǔ)給流量隨距離的變化規(guī)律見(jiàn)圖8,可以看出:相同時(shí)刻起始點(diǎn)x=0處的側(cè)滲補(bǔ)給流量為上游(1—1′上蘭村斷面)>中游(2—2′汾河大橋斷面)>下游(3—3′柴村橋斷面)。①1—1′上蘭村斷面(見(jiàn)圖8(a)),5 d時(shí)的側(cè)滲補(bǔ)給流量在x=0處從13 m3/d開(kāi)始快速下降,在1 750 m處降為0。90 d時(shí),x=0處的側(cè)滲補(bǔ)給流量為3.09 m3/d,在與河岸線相距3 000 m處降為0.46 m3/d,結(jié)合其他時(shí)間曲線變化趨勢(shì),可以推測(cè)該斷面90 d時(shí)的側(cè)滲補(bǔ)給流量在距離河岸線足夠遠(yuǎn)處降為0。②2—2′汾河大橋斷面(見(jiàn)圖8(b)),5 d時(shí)的側(cè)滲補(bǔ)給流量在x=0處從11.4 m3/d開(kāi)始快速下降,在1 692 m處降為0。90 d時(shí),x=0處側(cè)滲補(bǔ)給流量為2.7 m3/d,在與河岸線相距3 000 m處降為0。③3—3′柴村橋斷面(見(jiàn)圖8(c)),5 d時(shí)的側(cè)滲補(bǔ)給流量在x=0處從10.1 m3/d開(kāi)始快速下降,在1 470 m處降為0。90 d時(shí),x=0處側(cè)滲補(bǔ)給流量為2.4 m3/d,在與河岸線相距3 000 m處降為0.2 m3/d。

    (3)距河岸線不同距離側(cè)滲補(bǔ)給流量隨時(shí)間的變化規(guī)律見(jiàn)圖9。①1—1′上蘭村斷面(見(jiàn)圖9(a)),距河岸線20、100 m處的側(cè)滲補(bǔ)給流量首先從0開(kāi)始快速上升,至6 d達(dá)到峰值,分別為11、10.6 m3/d,隨后再緩慢下降;距河岸線500 m處的側(cè)滲補(bǔ)給流量從0開(kāi)始快速上升,至8 d達(dá)到峰值6.7 m3/d;距河岸線1 000 m處的側(cè)滲補(bǔ)給流量從0開(kāi)始快速上升,至28 d達(dá)到峰值3.2 m3/d;距河岸線2 000、3 000 m處的側(cè)滲補(bǔ)給流量從0開(kāi)始緩慢上升,至90 d時(shí)尚未達(dá)到峰值。②2—2′汾河大橋斷面(見(jiàn)圖9(b)),與河岸線相距20、100、500、1 000 m處的側(cè)滲補(bǔ)給流量在90 d內(nèi)首先從0開(kāi)始上升,達(dá)到峰值后再緩慢下降,峰值分別為9.5、9.2、4.8、2.1 m3/d;該斷面與河岸線相距2 000、3 000 m處的側(cè)滲補(bǔ)給流量在90 d內(nèi)未達(dá)到峰值。③3—3′柴村橋斷面(見(jiàn)圖9(c)),側(cè)滲補(bǔ)給流量隨時(shí)間的變化規(guī)律與1—1′上蘭村斷面、2—2′汾河大橋斷面大致相同,該斷面20、100、500、1 000 m處的峰值分別為8.48、8.07、3.96、1.82 m3/d。

    (4)河流水位抬高后在t時(shí)間內(nèi),起始點(diǎn)x=0處,總的單寬側(cè)滲流量Qt(單位長(zhǎng)度上一側(cè)河流補(bǔ)給地下水的總量)為

    總的單寬側(cè)滲流量隨時(shí)間逐漸增大,并且相同時(shí)間內(nèi),上、中、下游各斷面總的單寬側(cè)滲流量排序?yàn)?—1′上蘭村斷面>2—2′汾河大橋斷面>3—3′柴村橋斷面。0 d時(shí)各斷面總的單寬側(cè)滲流量均為0,90 d內(nèi)上、中、下游3個(gè)斷面總的單寬側(cè)滲流量分別為552、484、426 m3(見(jiàn)圖10)。

    4?結(jié)?論

    (1)Δh的變化規(guī)律反映河流水位瞬時(shí)抬升1.5 m后對(duì)潛水水位的補(bǔ)給過(guò)程。距離河岸線較近的位置先接收河流水體補(bǔ)給,Δh在0~20 m范圍內(nèi)減小明顯,在20 m之后緩慢減小直至為0,同一位置處Δh隨時(shí)間的延長(zhǎng)而增大。

    (2)根據(jù)不同時(shí)間Δh隨距離的變化規(guī)律得出相同時(shí)間的影響距離:1—1′上蘭村斷面最遠(yuǎn),2—2′汾河大橋斷面與3—3′柴村橋斷面相近。

    (3)同一位置處的側(cè)滲速度和側(cè)滲補(bǔ)給流量在接收河流水體補(bǔ)給后先增大達(dá)到最大值,隨后隨時(shí)間延長(zhǎng)緩慢減小,相同時(shí)刻起始點(diǎn)x=0處的側(cè)滲速度排序?yàn)樯嫌?中游<下游,側(cè)滲補(bǔ)給流量排序?yàn)樯嫌?中游>下游;同一時(shí)刻的側(cè)滲速度和側(cè)滲補(bǔ)給流量隨與河岸線距離逐漸增大而減小。河流水位抬高后起始點(diǎn)x=0處總的單寬側(cè)滲流量隨時(shí)間延長(zhǎng)逐漸增大,相同時(shí)刻上、中、下游各斷面總的單寬側(cè)滲流量排序?yàn)?—1′上蘭村斷面>2—2′汾河大橋斷面>3—3′柴村橋斷面。

    參考文獻(xiàn):

    [1]?杜向潤(rùn),張建龍.實(shí)施汾河流域生態(tài)修復(fù),落實(shí)中央生態(tài)文明建設(shè)宏大實(shí)踐[J].山西水利,2016(10):4,7.

    [2]?李海軍.汾河流域生態(tài)修復(fù)可行性研究[J].山西水利,2016(11):5-6.

    [3]?劉建榮.汾河流域濕地動(dòng)態(tài)分析與恢復(fù)對(duì)策研究[J].山西林業(yè)科技,2016,45(3):14-16.

    [4]?陳剛.汾河流域生態(tài)環(huán)境現(xiàn)狀及原因分析[J].山西水利,2016(10):10-11,19.

    [5]?金光球,李凌.河流中潛流交換研究進(jìn)展[J].水科學(xué)進(jìn)展,2008,19(2):285-293.

    [6]?ELLIOTT Alexander H,BROOKS Norman H. Transfer of Nonsorbing Solutes to a Streambed with Bed Forms: Laboratory Experiments[J]. Water Resources Research,1997,33(1): 137-151.

    [7]?SALEHIN Mashfiqus,PACKMAN Aaron I,PARADIS Matthew. Hyporheic Exchange with Heterogeneous Streambeds:Laboratory Experiments and Modeling[J]. Water Resources Research,2004,40(11):W11504.

    [8]?HUNT Bruce. Unsteady Stream Depletion from Ground Water Pumping[J].Ground Water,1999,37(1):98-102.

    [9]?PACKMAN Aaron I,SALEHIN Mashfiqus. Relative Roles of Stream Flow and Sedimentary Conditions in Controlling Hyporheic Exchange[J].Hydrobiologia,2003,494:291-297.

    [10]?HARVEY Judson W,BENCALA Kenneth E.The Effect of Streambed Topography on Surface-Subsurface Water Exchange in Mountain Catchments[J].Water Resources Research,1993,29(1):89-98.

    [11]?LIERMANN Catherine Reidy. Variability of Hyporheic Zones in Puget Sound Lowland Streams[D].Washington:University of Washington,2004:70.

    [12]?束龍倉(cāng),魯程鵬,李偉.考慮參數(shù)不確定性的地表水與地下水交換量的計(jì)算方法[J].水文地質(zhì)工程地質(zhì),2008,35(5):68-71.

    [13]?ROSENBERRY Donald O,PITLICK John.Local-Scale Variability of Seepage and Hydraulic Conductivity in a Shallow Gravel-Bed River[J]. Hydrological Processes,2009,23(23):3306-3318.

    [14]?POZDNIAKOV Sergey P,WANG Ping,LEKHOV Mikhail V. A Semi-Analytical Generalized Hvorslev Formula for Estimating Riverbed Hydraulic Conductivity with an Open-Ended Standpipe Permeameter[J]. Journal of Hydrology,2016,540:736-743.

    [15]?束龍倉(cāng), CHEN Xunhong.美國(guó)內(nèi)布拉斯加州普拉特河河床沉積物滲透系數(shù)的現(xiàn)場(chǎng)測(cè)定[J].水科學(xué)進(jìn)展,2002,13(5):629-633.

    [16]?CONSTANTZ Jim.Streambeds Merit Recognition as a Scientific Discipline[J].Wiley Interdisciplinary Reviews:Water,2016,3(1):13-18.

    [17]?BRUNNER Philip,THERRIEN Rene,RENARD Philippe,et al. Advances in Understanding River-Groundwater Interactions[J]. Reviews of Geophysics,2017,55(3):818-854.

    [18]?ZHANG Guotao,SONG Jinxi,WEN Ming,et al. Effect of Bank Curvatures on Hyporheic Water Exchange at Meter Scale[J]. Hydrology Research,2017,48(2):355-369.

    [19]?侯宏冰,劉長(zhǎng)禮,葉浩,等.黃河下游懸河溫縣段河水位變化對(duì)側(cè)滲的影響研究[J]. 南水北調(diào)與水利科技,2007,5(5):109-112.

    [20]?XI Haiyang,F(xiàn)ENG Qi,ZHANG Lu,et al. Groundwater Storage Changes and Estimation of Stream Lateral Seepage to Groundwater in Desert Riparian Forest Region[J]. Hydrology Research,2017,49(3):861-877.

    [21]?YAO Yingying,ZHENG Chunmiao,TIAN Yong,et al. Numerical Modeling of Regional Groundwater Flow in the Heihe River Basin,China:Advances and New Insights[J]. Science China Earth Sciences,2015,58(1):3-15.

    [22]?SHU Longcang,CHEN Xunhong. Simulation of Water Quantity Exchange Between Groundwater and the Platte River Water,Central Nebraska[J]. Journal of Central South University,2002,9(3):212-215.

    [23]?CHEN Xunhong,SHU Longcang. Stream-Aquifer Interactions: Evaluation of Depletion Volume and Residual Effectsfrom Groundwater Pumping[J]. Ground Water,2002,40(3):284-290.

    [24]?束龍倉(cāng),CHEN Xunhong.地下水開(kāi)采對(duì)河流流量衰減的影響分析[J].水利學(xué)報(bào),2003,34(2):112-116.

    [25]?YAO Yingying,ZHENG Chunmiao,LIU Jie,et al. Conceptual and Numerical Models for Groundwater Flow in an Arid Inland River Basin[J].Hydrological Processes,2015,29(6):1480-1492.

    [26]?吳吉春,薛禹群.地下水動(dòng)力學(xué)[M].北京:中國(guó)水利水電出版社,2012:74-79.

    【責(zé)任編輯?張華興】

    国产91精品成人一区二区三区 | 中文字幕精品免费在线观看视频| 悠悠久久av| 国产av一区二区精品久久| 亚洲精品粉嫩美女一区| 香蕉国产在线看| 一级毛片电影观看| 老司机在亚洲福利影院| 精品视频人人做人人爽| 国产免费福利视频在线观看| 亚洲精品乱久久久久久| 精品国产国语对白av| 中文字幕制服av| 亚洲视频免费观看视频| 欧美在线一区亚洲| 熟女少妇亚洲综合色aaa.| 国产免费视频播放在线视频| 国产精品久久久久久精品古装| 国产精品久久久人人做人人爽| 国产伦人伦偷精品视频| av又黄又爽大尺度在线免费看| 在线观看免费高清a一片| 亚洲精品第二区| 久久国产精品男人的天堂亚洲| 精品国产一区二区三区久久久樱花| 别揉我奶头~嗯~啊~动态视频 | 午夜福利在线免费观看网站| 日韩制服骚丝袜av| 亚洲国产欧美网| 丝袜在线中文字幕| 日本精品一区二区三区蜜桃| 亚洲熟女精品中文字幕| 亚洲精品av麻豆狂野| 亚洲精品在线美女| 亚洲av美国av| 欧美日韩一级在线毛片| 热re99久久国产66热| 国产成人欧美| av有码第一页| 国产麻豆69| 免费在线观看视频国产中文字幕亚洲 | 亚洲精品在线美女| 亚洲精品粉嫩美女一区| 国产av又大| 中文字幕最新亚洲高清| 国产日韩欧美视频二区| 欧美黑人欧美精品刺激| 欧美日韩一级在线毛片| 咕卡用的链子| av在线老鸭窝| a级毛片在线看网站| av网站在线播放免费| 老司机亚洲免费影院| 一个人免费看片子| 日韩人妻精品一区2区三区| 成人三级做爰电影| 国产精品1区2区在线观看. | 亚洲欧美激情在线| 国产精品二区激情视频| 亚洲人成电影免费在线| 国产av国产精品国产| 悠悠久久av| 日韩,欧美,国产一区二区三区| 久久九九热精品免费| 精品一区在线观看国产| 天堂俺去俺来也www色官网| 18禁国产床啪视频网站| 一区在线观看完整版| 淫妇啪啪啪对白视频 | 母亲3免费完整高清在线观看| 亚洲国产欧美在线一区| 老司机亚洲免费影院| 国产又爽黄色视频| 纵有疾风起免费观看全集完整版| av天堂久久9| 麻豆乱淫一区二区| 久久久久精品人妻al黑| e午夜精品久久久久久久| 99久久99久久久精品蜜桃| 丝瓜视频免费看黄片| 国产成人精品在线电影| av又黄又爽大尺度在线免费看| 老司机影院成人| 一区福利在线观看| 日本av免费视频播放| 黄色a级毛片大全视频| 国产无遮挡羞羞视频在线观看| 免费在线观看完整版高清| 汤姆久久久久久久影院中文字幕| 50天的宝宝边吃奶边哭怎么回事| 亚洲 国产 在线| 国产免费现黄频在线看| av片东京热男人的天堂| 777久久人妻少妇嫩草av网站| 国产真人三级小视频在线观看| 久久精品国产亚洲av高清一级| 美女扒开内裤让男人捅视频| 三级国产精品欧美在线观看 | 最近在线观看免费完整版| 国产精品电影一区二区三区| 亚洲乱码一区二区免费版| 午夜亚洲福利在线播放| 九色国产91popny在线| 悠悠久久av| 亚洲人成77777在线视频| 久久这里只有精品19| 久久精品国产综合久久久| 无遮挡黄片免费观看| 岛国在线免费视频观看| 国产成人av激情在线播放| 12—13女人毛片做爰片一| 亚洲狠狠婷婷综合久久图片| 日韩欧美在线二视频| 免费无遮挡裸体视频| 99久久精品热视频| 女警被强在线播放| 欧美日韩精品网址| 亚洲精品美女久久av网站| 狂野欧美白嫩少妇大欣赏| 搡老熟女国产l中国老女人| 日韩成人在线观看一区二区三区| 人成视频在线观看免费观看| 亚洲人成网站高清观看| 大型黄色视频在线免费观看| 人成视频在线观看免费观看| 中文字幕久久专区| av天堂在线播放| 亚洲七黄色美女视频| 日韩欧美在线二视频| 欧美三级亚洲精品| 国产精品精品国产色婷婷| 亚洲欧美精品综合久久99| 久热爱精品视频在线9| 美女高潮喷水抽搐中文字幕| 欧美3d第一页| 999久久久精品免费观看国产| 变态另类丝袜制服| 很黄的视频免费| 少妇熟女aⅴ在线视频| 精品国产美女av久久久久小说| 制服诱惑二区| 搞女人的毛片| 国产蜜桃级精品一区二区三区| 亚洲成av人片免费观看| 欧美 亚洲 国产 日韩一| 女人高潮潮喷娇喘18禁视频| 亚洲精品粉嫩美女一区| 亚洲精品国产一区二区精华液| 午夜日韩欧美国产| 久久人人精品亚洲av| 免费人成视频x8x8入口观看| 国产av麻豆久久久久久久| 国产男靠女视频免费网站| 亚洲aⅴ乱码一区二区在线播放 | 69av精品久久久久久| 91大片在线观看| 亚洲专区中文字幕在线| 久久亚洲真实| 免费在线观看日本一区| 国产男靠女视频免费网站| 啪啪无遮挡十八禁网站| 成人精品一区二区免费| 久久亚洲精品不卡| 国产成人一区二区三区免费视频网站| 变态另类成人亚洲欧美熟女| 最近视频中文字幕2019在线8| 亚洲男人的天堂狠狠| 精品一区二区三区视频在线观看免费| 18禁美女被吸乳视频| 级片在线观看| 无遮挡黄片免费观看| 久久久久国内视频| 亚洲一区二区三区不卡视频| 亚洲国产欧洲综合997久久,| 亚洲精品国产一区二区精华液| 亚洲中文字幕一区二区三区有码在线看 | 一个人免费在线观看的高清视频| 看免费av毛片| 婷婷精品国产亚洲av在线| 亚洲七黄色美女视频| 欧美久久黑人一区二区| 一级a爱片免费观看的视频| 国产成人影院久久av| 午夜日韩欧美国产| 88av欧美| 美女扒开内裤让男人捅视频| 午夜免费成人在线视频| 一本一本综合久久| 日韩国内少妇激情av| 国产精品免费一区二区三区在线| 老司机午夜福利在线观看视频| 日韩中文字幕欧美一区二区| 最近最新中文字幕大全免费视频| 亚洲av中文字字幕乱码综合| 99re在线观看精品视频| 91国产中文字幕| 欧美乱码精品一区二区三区| 黄片小视频在线播放| 99久久精品国产亚洲精品| 曰老女人黄片| 成在线人永久免费视频| 国产在线精品亚洲第一网站| 五月伊人婷婷丁香| 久9热在线精品视频| 亚洲最大成人中文| 久久久精品欧美日韩精品| svipshipincom国产片| 国产高清视频在线播放一区| 又大又爽又粗| 亚洲va日本ⅴa欧美va伊人久久| 久久久久久久久久黄片| 欧美另类亚洲清纯唯美| 国产精品久久久久久亚洲av鲁大| 精品欧美国产一区二区三| 欧美绝顶高潮抽搐喷水| 女生性感内裤真人,穿戴方法视频| 老司机福利观看| 身体一侧抽搐| 国模一区二区三区四区视频 | 国产av在哪里看| 久久久久久久精品吃奶| 中文资源天堂在线| 亚洲18禁久久av| 久久精品国产99精品国产亚洲性色| 久久香蕉精品热| 我要搜黄色片| 精品高清国产在线一区| 一二三四社区在线视频社区8| 中出人妻视频一区二区| 亚洲全国av大片| 一级片免费观看大全| 真人做人爱边吃奶动态| 欧美黄色片欧美黄色片| www.熟女人妻精品国产| 一级a爱片免费观看的视频| 午夜精品久久久久久毛片777| 日韩精品中文字幕看吧| 国产亚洲欧美在线一区二区| 国产精品一区二区免费欧美| 国产精品久久久久久精品电影| 精品一区二区三区四区五区乱码| 日韩中文字幕欧美一区二区| 色精品久久人妻99蜜桃| 岛国在线观看网站| 久久中文字幕一级| 亚洲av中文字字幕乱码综合| 亚洲熟妇中文字幕五十中出| 在线播放国产精品三级| 一a级毛片在线观看| 操出白浆在线播放| 三级国产精品欧美在线观看 | 亚洲第一电影网av| 国产精品久久久av美女十八| 国内精品久久久久久久电影| 一进一出抽搐gif免费好疼| 一个人观看的视频www高清免费观看 | 九色国产91popny在线| 成人三级黄色视频| 日本一区二区免费在线视频| 黄色 视频免费看| 99久久无色码亚洲精品果冻| 超碰成人久久| 国产一区二区在线av高清观看| 日韩欧美在线二视频| 国产三级在线视频| 久久久精品国产亚洲av高清涩受| 在线观看免费视频日本深夜| av超薄肉色丝袜交足视频| 91在线观看av| 久久久国产精品麻豆| 欧美久久黑人一区二区| 久久精品国产清高在天天线| 亚洲国产精品久久男人天堂| av在线播放免费不卡| 亚洲欧美激情综合另类| 欧美+亚洲+日韩+国产| a在线观看视频网站| 在线十欧美十亚洲十日本专区| videosex国产| 窝窝影院91人妻| 国产91精品成人一区二区三区| 色综合欧美亚洲国产小说| 97碰自拍视频| 国产精品美女特级片免费视频播放器 | 亚洲熟女毛片儿| av有码第一页| 俺也久久电影网| cao死你这个sao货| 99久久精品国产亚洲精品| 国产精品国产高清国产av| 亚洲男人天堂网一区| 在线看三级毛片| 欧美成人午夜精品| 舔av片在线| 最新在线观看一区二区三区| 熟妇人妻久久中文字幕3abv| 精品国产美女av久久久久小说| 少妇的丰满在线观看| 好男人电影高清在线观看| 午夜精品在线福利| 动漫黄色视频在线观看| 午夜免费激情av| 亚洲午夜精品一区,二区,三区| 俄罗斯特黄特色一大片| 在线国产一区二区在线| 午夜福利免费观看在线| 人妻丰满熟妇av一区二区三区| 国产免费av片在线观看野外av| 精品久久蜜臀av无| 97人妻精品一区二区三区麻豆| 国产成人一区二区三区免费视频网站| 久久精品人妻少妇| 18禁黄网站禁片午夜丰满| 老汉色∧v一级毛片| 丝袜人妻中文字幕| 在线观看66精品国产| 最近最新中文字幕大全免费视频| 欧美精品啪啪一区二区三区| 18禁美女被吸乳视频| 久久久久免费精品人妻一区二区| or卡值多少钱| 韩国av一区二区三区四区| 国产精品九九99| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品美女久久久久99蜜臀| 国产aⅴ精品一区二区三区波| 国产一区在线观看成人免费| av片东京热男人的天堂| 色在线成人网| 老司机午夜福利在线观看视频| 国产三级在线视频| 女人被狂操c到高潮| 午夜免费成人在线视频| 成年免费大片在线观看| 村上凉子中文字幕在线| 成人国语在线视频| 亚洲国产精品999在线| 亚洲人成伊人成综合网2020| 日韩免费av在线播放| 制服诱惑二区| 午夜视频精品福利| 狂野欧美白嫩少妇大欣赏| 一区二区三区高清视频在线| 午夜精品久久久久久毛片777| 国产一区二区三区视频了| 成人手机av| 亚洲 欧美 日韩 在线 免费| 在线国产一区二区在线| 午夜福利免费观看在线| 韩国av一区二区三区四区| 久久久精品欧美日韩精品| 中文字幕av在线有码专区| 九色成人免费人妻av| 免费在线观看日本一区| 久久久久久久午夜电影| 久久久精品大字幕| 久久久水蜜桃国产精品网| 中文字幕久久专区| 床上黄色一级片| 一级毛片女人18水好多| 国产三级黄色录像| 久久久久久久久中文| 亚洲无线在线观看| 女人爽到高潮嗷嗷叫在线视频| 一个人免费在线观看的高清视频| 国产亚洲精品一区二区www| 99国产综合亚洲精品| 天堂av国产一区二区熟女人妻 | 在线观看免费午夜福利视频| 国产午夜福利久久久久久| 男人的好看免费观看在线视频 | 黄色片一级片一级黄色片| 久久久久久九九精品二区国产 | 国产真人三级小视频在线观看| 精品欧美一区二区三区在线| 久久亚洲精品不卡| 国产aⅴ精品一区二区三区波| 亚洲成人国产一区在线观看| 亚洲欧美日韩高清在线视频| 琪琪午夜伦伦电影理论片6080| 真人做人爱边吃奶动态| 叶爱在线成人免费视频播放| 亚洲最大成人中文| 精品久久久久久久末码| 国产精品亚洲一级av第二区| 嫁个100分男人电影在线观看| 中文亚洲av片在线观看爽| 十八禁人妻一区二区| 夜夜看夜夜爽夜夜摸| 国产亚洲av嫩草精品影院| 国产一级毛片七仙女欲春2| 亚洲欧美日韩高清专用| 一区二区三区高清视频在线| 老司机午夜十八禁免费视频| 搡老岳熟女国产| 亚洲欧美精品综合久久99| 色噜噜av男人的天堂激情| 欧美日韩亚洲国产一区二区在线观看| 他把我摸到了高潮在线观看| 亚洲av成人精品一区久久| АⅤ资源中文在线天堂| 国产成人精品久久二区二区免费| 一个人免费在线观看电影 | avwww免费| 日本在线视频免费播放| 亚洲18禁久久av| 看免费av毛片| 国产精品精品国产色婷婷| 亚洲自偷自拍图片 自拍| 最近最新中文字幕大全免费视频| 在线a可以看的网站| 午夜福利视频1000在线观看| 999久久久国产精品视频| 久久久久性生活片| 中文字幕人成人乱码亚洲影| 黄色丝袜av网址大全| 99热只有精品国产| 成人欧美大片| 久久中文字幕人妻熟女| 久久久久国产精品人妻aⅴ院| 一个人观看的视频www高清免费观看 | 精品熟女少妇八av免费久了| 黄色成人免费大全| 亚洲aⅴ乱码一区二区在线播放 | 男女视频在线观看网站免费 | 亚洲黑人精品在线| www日本在线高清视频| 免费观看人在逋| 香蕉久久夜色| 国产一区二区三区在线臀色熟女| 成人一区二区视频在线观看| 欧美黄色片欧美黄色片| 国产精品 欧美亚洲| 日韩 欧美 亚洲 中文字幕| 久久这里只有精品19| 久久香蕉精品热| 国产精品久久久久久人妻精品电影| 国产熟女xx| 性色av乱码一区二区三区2| 亚洲av五月六月丁香网| а√天堂www在线а√下载| 国产1区2区3区精品| a级毛片a级免费在线| 一本大道久久a久久精品| 欧美在线一区亚洲| 亚洲中文字幕一区二区三区有码在线看 | 97超级碰碰碰精品色视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 在线播放国产精品三级| 麻豆国产97在线/欧美 | videosex国产| 国内精品一区二区在线观看| 一个人免费在线观看电影 | 亚洲专区字幕在线| 亚洲精品国产精品久久久不卡| 最新美女视频免费是黄的| av欧美777| 中出人妻视频一区二区| 欧美日本亚洲视频在线播放| 亚洲精品国产一区二区精华液| 美女黄网站色视频| 久久久久久久久中文| 亚洲七黄色美女视频| 国产视频一区二区在线看| 别揉我奶头~嗯~啊~动态视频| 国产男靠女视频免费网站| 一卡2卡三卡四卡精品乱码亚洲| 色综合亚洲欧美另类图片| 亚洲 欧美一区二区三区| 欧美日韩亚洲综合一区二区三区_| 国产乱人伦免费视频| tocl精华| 人人妻人人澡欧美一区二区| 日本三级黄在线观看| 亚洲熟妇中文字幕五十中出| 成人永久免费在线观看视频| 欧美成人一区二区免费高清观看 | 国产黄a三级三级三级人| 91九色精品人成在线观看| 日韩av在线大香蕉| 999久久久精品免费观看国产| 很黄的视频免费| 在线观看一区二区三区| 精品一区二区三区av网在线观看| 精品久久久久久久毛片微露脸| 岛国在线观看网站| 亚洲色图 男人天堂 中文字幕| 此物有八面人人有两片| 亚洲精品国产一区二区精华液| 中文字幕人妻丝袜一区二区| 久久精品亚洲精品国产色婷小说| 这个男人来自地球电影免费观看| 好男人电影高清在线观看| 老司机福利观看| 国产成年人精品一区二区| 亚洲国产日韩欧美精品在线观看 | 香蕉丝袜av| 在线观看66精品国产| 在线观看免费午夜福利视频| 国产不卡一卡二| 99热这里只有是精品50| 国产精品av久久久久免费| 国语自产精品视频在线第100页| 1024手机看黄色片| 成人三级黄色视频| 国产精品一区二区免费欧美| 免费看美女性在线毛片视频| 亚洲中文字幕日韩| 亚洲成a人片在线一区二区| 国产精品久久久人人做人人爽| 天堂av国产一区二区熟女人妻 | 欧美乱码精品一区二区三区| 久久久久国内视频| 午夜免费成人在线视频| 亚洲欧洲精品一区二区精品久久久| 视频区欧美日本亚洲| 亚洲国产精品合色在线| 午夜精品久久久久久毛片777| 美女免费视频网站| 啦啦啦观看免费观看视频高清| 免费看美女性在线毛片视频| 日本黄色视频三级网站网址| 夜夜夜夜夜久久久久| 他把我摸到了高潮在线观看| 亚洲午夜理论影院| 国产精品久久久久久人妻精品电影| 欧美一区二区国产精品久久精品 | 美女午夜性视频免费| 亚洲人成网站高清观看| 免费看美女性在线毛片视频| 精品国产亚洲在线| 狂野欧美激情性xxxx| 亚洲男人天堂网一区| 成人永久免费在线观看视频| 91国产中文字幕| 亚洲自拍偷在线| 国产精品一区二区免费欧美| 性色av乱码一区二区三区2| 不卡一级毛片| 成人三级做爰电影| 国产一区二区激情短视频| 亚洲午夜精品一区,二区,三区| 欧美日韩福利视频一区二区| 色噜噜av男人的天堂激情| 成人欧美大片| 日本 欧美在线| 日韩精品青青久久久久久| 国产av一区在线观看免费| 亚洲狠狠婷婷综合久久图片| 熟女电影av网| 国产69精品久久久久777片 | 黄片大片在线免费观看| 亚洲美女黄片视频| 国产成人av激情在线播放| а√天堂www在线а√下载| 免费在线观看影片大全网站| 欧美极品一区二区三区四区| 亚洲熟妇熟女久久| 国产黄片美女视频| 热99re8久久精品国产| 少妇被粗大的猛进出69影院| 亚洲av片天天在线观看| 嫁个100分男人电影在线观看| 亚洲自拍偷在线| 黄色a级毛片大全视频| 18禁美女被吸乳视频| 日韩欧美 国产精品| 国产单亲对白刺激| 国产欧美日韩精品亚洲av| 国产免费男女视频| 欧美乱妇无乱码| 国产成人系列免费观看| 50天的宝宝边吃奶边哭怎么回事| 亚洲av第一区精品v没综合| 午夜视频精品福利| www.熟女人妻精品国产| 久久人妻福利社区极品人妻图片| 欧美日韩黄片免| 欧美一级毛片孕妇| 国产午夜精品久久久久久| 一二三四在线观看免费中文在| 两人在一起打扑克的视频| 88av欧美| 曰老女人黄片| 在线看三级毛片| 国产高清视频在线观看网站| 白带黄色成豆腐渣| 免费在线观看亚洲国产| 欧美黑人欧美精品刺激| 亚洲午夜精品一区,二区,三区| 精品一区二区三区四区五区乱码| 在线免费观看的www视频| 老司机深夜福利视频在线观看| 欧美日韩一级在线毛片| 一进一出好大好爽视频| 欧美 亚洲 国产 日韩一| 可以在线观看的亚洲视频| 色综合亚洲欧美另类图片| 欧美乱色亚洲激情| 欧美日韩福利视频一区二区| 香蕉av资源在线| 亚洲欧美日韩东京热| 99久久国产精品久久久| 国产精品,欧美在线| 757午夜福利合集在线观看| 中文字幕精品亚洲无线码一区| 国产私拍福利视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲av中文字字幕乱码综合| 一边摸一边抽搐一进一小说| 日韩欧美精品v在线| 真人做人爱边吃奶动态| 国产精华一区二区三区| 国内揄拍国产精品人妻在线| 中文资源天堂在线|