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

    單井抽出-回滲同步循環(huán)地下水水力控制模型研究

    2023-02-04 15:00:00王新港徐祥健姜永海
    環(huán)境科學(xué)研究 2023年1期
    關(guān)鍵詞:滲透系數(shù)水力含水層

    王新港,楊 昱,徐祥健,韓 旭,夏 甫,鄧 圣,肖 瀚*,姜永海

    1. 中國(guó)環(huán)境科學(xué)研究院,國(guó)家環(huán)境保護(hù)地下水污染模擬與控制重點(diǎn)實(shí)驗(yàn)室,北京 100012

    2. 河北工程大學(xué)能源與環(huán)境工程學(xué)院,河北 邯鄲 056038

    單井抽出-回滲同步循環(huán)地下水水力控制技術(shù)是從地下水污染羽范圍內(nèi)采用單井抽出污染地下水,經(jīng)過(guò)處理后,在抽水井周邊同步等量回滲至包氣帶和含水層,通過(guò)抽出-回滲同步循環(huán),人為改變污染羽及其周邊區(qū)域地下水流場(chǎng)和水力梯度,將受污染的地下水控制在一定范圍內(nèi),防止污染物進(jìn)一步擴(kuò)散蔓延,進(jìn)而將受污染的地下水體與未受污染的清潔地下水體分隔開(kāi),是地下水污染修復(fù)的重要手段[1-5]. 該技術(shù)的優(yōu)勢(shì)是操作簡(jiǎn)單、運(yùn)行成本較低、加速地下水循環(huán)流動(dòng)縮短地下水修復(fù)時(shí)間,影響該技術(shù)實(shí)施效果的重要因素是水動(dòng)力條件以及含水層水文地質(zhì)條件,水動(dòng)力條件主要由抽注水量控制,而含水層水文地質(zhì)條件主要指含水層滲透性[6-13].

    在地下水污染水力控制技術(shù)領(lǐng)域內(nèi),已有研究主要關(guān)注抽注技術(shù)參數(shù)的選取以及抽注水井布局和抽注水量的優(yōu)化對(duì)水力控制效果的影響,使抽注井的布局和抽注水量更為合理,從而縮短抽注時(shí)間、提高水力控制效率[14-19]. 如Bear等[20]定量地研究了注水井補(bǔ)給對(duì)含水層地下水運(yùn)動(dòng)規(guī)律的影響,為水力控制技術(shù)的發(fā)展奠定了理論基礎(chǔ);Shar等[21]、Christ等[22]采用解析、半解析法分別研究了相同流量下任意雙抽水井及不同流量下多抽水井水力控制區(qū)域的邊界形態(tài)特征,為數(shù)值模擬提供了理論基礎(chǔ);隨后,Wang等[23]用數(shù)值模擬和優(yōu)化耦合方法,采用拉普拉斯變換有限單元法對(duì)均勻等厚污染含水層的多個(gè)抽注井進(jìn)行位置布設(shè),提高了捕獲污染的有效性;Matott等[24]采用分析單元法對(duì)井群系統(tǒng)的抽水井?dāng)?shù)量、位置和抽水量進(jìn)行了優(yōu)化;劉苑[25]和萬(wàn)鵬[26]均采用MGO模塊結(jié)合遺傳算法對(duì)水力控制系統(tǒng)進(jìn)行優(yōu)化,通過(guò)約束抽水井?dāng)?shù)量、抽水流量、水頭和抽水總量等參數(shù),提高了水力控制能力. 目前地下水污染水力控制技術(shù)研究中數(shù)值模擬成為主要解決手段,常用的數(shù)值模擬軟件主要包括Groundwater Modeling Systems(GMS)、Visual MODFLOW、FEFLOW及Groundwater Vistas等[27-28]. 然而單井抽出-回滲技術(shù)特點(diǎn)是抽出-回滲循環(huán)同步進(jìn)行,通過(guò)人工干擾/控制實(shí)現(xiàn)水力控制下抽出-回滲同步進(jìn)行的研究相對(duì)較少[29].

    研究采取物理模型和數(shù)值模型相結(jié)合的方法,搭建實(shí)驗(yàn)室尺度砂柱孔隙含水層物理模型,基于該物理模型構(gòu)建MODFLOW地下水流數(shù)值模型并校準(zhǔn),使模擬地下水流可以較好地刻畫(huà)實(shí)際地下水流;利用校準(zhǔn)后的數(shù)值模型模擬不同抽注水量和含水層滲透性條件下實(shí)現(xiàn)水力控制時(shí)不同位置地下水水位的差異,基于潛水完整井穩(wěn)定流公式-裘布依公式的數(shù)據(jù)擬合,分析實(shí)現(xiàn)水力控制時(shí)含水層不同位置水位與抽注水量和含水層滲透系數(shù)的關(guān)系以及回滲半徑對(duì)水力控制的影響. 該研究論證了單井抽出-回滲同步循環(huán)地下水水力控制技術(shù)實(shí)現(xiàn)人工控制地下水流場(chǎng)和水力梯度的可行性,并量化了水動(dòng)力條件和水文地質(zhì)條件的影響,以期為該技術(shù)在實(shí)際場(chǎng)地的應(yīng)用提供重要參考.

    1 材料與方法

    1.1 物理模型

    1.1.1裝置搭建

    單井抽出-回滲同步循環(huán)地下水水力控制技術(shù)試驗(yàn)裝置如圖1所示. 試驗(yàn)裝置由砂柱、布設(shè)于砂柱中軸線的抽水管、砂柱表層的布水管以及砂柱含水層內(nèi)部不同位置、不同深度的測(cè)壓管組成.

    圖1 單井抽出-回滲同步循環(huán)地下水水力控制技術(shù)試驗(yàn)裝置Fig.1 Schematic diagram of single well pumpingrecharge synchronous cyclical groundwater hydraulic control technique device

    砂柱由1 cm厚有機(jī)玻璃制成,直徑為100 cm,高度為100 cm. 砂柱內(nèi)均勻填充石英砂〔采購(gòu)于天津某石料廠的80~120目(100目=0.150 mm)石英砂〕,通過(guò)下部進(jìn)水且每5 cm進(jìn)行飽和夯實(shí)的方式裝填,構(gòu)建厚度為70 cm的均質(zhì)潛水含水層. 抽水管(長(zhǎng)100 cm、直徑2 cm、開(kāi)篩70 cm)垂直固定于砂柱中軸線位置并包裹160目過(guò)濾篩網(wǎng)防止石英砂進(jìn)入,管底位于砂柱底部(用于模擬抽水完整井),管頂與蠕動(dòng)泵相連接. 布水管以抽水管為圓心,半徑32 cm范圍內(nèi)螺旋狀均勻鋪設(shè)于砂層表面(布水管內(nèi)徑0.6 cm,每隔5 cm開(kāi)孔),一端與蠕動(dòng)泵相連接,使抽出的水等量同步均勻回滲至砂柱含水層中. 蠕動(dòng)泵型號(hào)為保定蘭格BT100-1L,泵頭為YII25,可調(diào)整流量范圍為0.1~180 mL/min.

    1.1.2測(cè)試方法

    裝置監(jiān)測(cè)點(diǎn)位如圖2所示,采用上密下疏、內(nèi)密外疏的方式,通過(guò)砂柱側(cè)面48個(gè)監(jiān)測(cè)孔在砂柱內(nèi)不同位置、不同深度安裝48個(gè)測(cè)壓管,用于監(jiān)測(cè)含水層水位. 垂向共8層,測(cè)壓深度分別為1、3、5、7、14、28、48、70 cm,每層設(shè)置6個(gè)測(cè)壓點(diǎn),分別與中心抽水井井軸相距2、4、8、14、22、32 cm. 測(cè)壓管讀數(shù)實(shí)時(shí)反映砂柱內(nèi)不同位置、不同深度含水層水位的變化. 砂柱物理模型的滲透系數(shù)和孔隙度分別由抽水試驗(yàn)、補(bǔ)水試驗(yàn)測(cè)算[30]. 抽水試驗(yàn)測(cè)定滲透系數(shù)的原理是基于潛水完整井穩(wěn)定流公式-裘布依公式[31]:

    圖2 單井抽出-回滲同步循環(huán)地下水水力控制技術(shù)試驗(yàn)裝置概念Fig.2 Conceptual diagram of single well pumpingrecharge synchronous cyclical groundwater hydraulic control technique device

    式中:Q為抽水井流量,m3/d;K為含水層滲透系數(shù),m/d;h為含水層中任意位置地下水水位,m;hw為抽水井水位,m;r為含水層中某處與抽水井井軸距離,m;rw為抽水井半徑,m.

    實(shí)驗(yàn)室尺度砂柱物理模型搭建完成后,給定抽水量為2.5 cm3/s,運(yùn)行開(kāi)展抽水試驗(yàn),獲取抽水穩(wěn)定后不同位置、不同深度觀測(cè)水位,計(jì)算得到砂柱模型填土的滲透系數(shù)為0.01 cm/s,孔隙度為0.4.

    1.2 數(shù)值模型

    1.2.1模擬軟件選取

    GMS是一個(gè)界面友好、處理功能強(qiáng)大的地下水?dāng)?shù)值模擬軟件[28,32]. MODFLOW模塊離散方法簡(jiǎn)單、求解方式多樣、采用了模塊化結(jié)構(gòu),廣泛應(yīng)用于地下水滲流和溶質(zhì)運(yùn)移[33]. 該研究采用GMS中的MODFLOW模塊構(gòu)建實(shí)驗(yàn)室尺度砂柱孔隙含水層三維地下水穩(wěn)定流模型. MODFLOW模塊的水流控制方程[31]如下:

    式中:W為水流源匯項(xiàng),m3/d;t為模擬預(yù)測(cè)的時(shí)間,d;h為地下水水位值,m;μs為含水層儲(chǔ)水率,m-1;Kx、Ky、Kz為不同軸向上的滲透系數(shù)值,m/d.

    1.2.2模型構(gòu)建

    數(shù)值模型參照上述物理模型等比例構(gòu)建,模擬單井抽出-回滲循環(huán)過(guò)程含水層流場(chǎng),4種模型構(gòu)建參數(shù)條件如表1所示. 模擬區(qū)域?yàn)樵囼?yàn)砂柱的石英砂填裝區(qū)(直徑100 cm、高度70 cm),水平方向剖分為1 cm ×1 cm的正方形網(wǎng)格(100行、100列),垂向剖分為4層,從上至下各層厚度分別為20、16、24、10 cm.

    表1 模型構(gòu)建參數(shù)Table 1 Parameters for model construction

    模擬區(qū)域下部及側(cè)向邊界設(shè)定為零流量邊界;模擬區(qū)域中心有抽水井排泄,設(shè)定為定流量邊界(抽水量分別設(shè)定為1、2.5、5、10 cm3/s);模擬區(qū)域上部邊界有回滲補(bǔ)給(補(bǔ)給范圍為以抽水井井軸為圓心,半徑32 cm的圓),設(shè)定為定流量邊界(回滲量分別設(shè)定為1、2.5、5、10 cm3/s,對(duì)應(yīng)的補(bǔ)給速率分別為3.11×10-4、7.78×10-4、1.56×10-3、3.11×10-3cm/s).

    1.3 數(shù)值模型校準(zhǔn)

    為保證模型模擬結(jié)果可靠,對(duì)模型進(jìn)行校準(zhǔn)驗(yàn)證. 物理模型的48個(gè)觀測(cè)點(diǎn)分別在抽出-回滲量為1、2.5、5、10 cm3/s達(dá)到抽出-回滲平衡時(shí)記錄實(shí)測(cè)水位,對(duì)比抽出-回滲量分別為1、2.5、5、10 cm3/s時(shí)數(shù)值模型的模擬水位與物理模型的實(shí)測(cè)水位,微調(diào)滲透系數(shù),重復(fù)多次試算,直至大部分觀測(cè)點(diǎn)模擬水位與實(shí)測(cè)水位的誤差達(dá)到模擬精度要求.

    數(shù)值模擬結(jié)果的準(zhǔn)確度通過(guò)納什效率系數(shù)(Nash-Sutcliffe efficiency coefficient, NSE)驗(yàn)證[34].

    式中:Qo為實(shí)測(cè)水位值,cm;Qm為模擬水位值,cm;為t時(shí)的實(shí)測(cè)水位值,為實(shí)測(cè)水位值的總平均值,cm.

    NSE取值為負(fù)無(wú)窮至1,NSE接近1,表示模型質(zhì)量好,模型可信度高;NSE接近0,表示模擬結(jié)果接近觀測(cè)值的平均值水平,即總體結(jié)果可信,但誤差大;NSE遠(yuǎn)遠(yuǎn)小于0,則模型不可信.

    數(shù)值模型校準(zhǔn)后,滲透系數(shù)調(diào)整為0.005 cm/s,模擬水位和實(shí)測(cè)水位(見(jiàn)表2)分散比較結(jié)果如圖3所示,納什效率系數(shù)達(dá)到0.88,說(shuō)明模型精度較為理想,可以較好地刻畫(huà)物理模型中的實(shí)際地下水流場(chǎng).

    表2 不同抽出-回滲量條件下模擬水位和實(shí)測(cè)水位對(duì)比Table 2 Comparison between simulated and observed groundwater levels under different pumping-return seepage conditions

    圖3 不同抽出-回滲量條件下模擬水位和實(shí)測(cè)水位的分散比較Fig.3 Comparison diagram between simulated and observed groundwater levels under different pumping-return seepage conditions

    2 結(jié)果與分析

    2.1 實(shí)現(xiàn)水力控制時(shí)的模擬流場(chǎng)

    當(dāng)抽出-回滲量分別為1、2.5、5、10 cm3/s時(shí),達(dá)到抽出-回滲平衡后,模擬流場(chǎng)三維空間分布如圖4所示;距離抽水井軸0、2、4、8、14、22、32 cm處的模擬水位以及抽水影響半徑見(jiàn)表3. 抽出-回滲量對(duì)水位降落漏斗的深度影響較大,抽出-回滲量越大,降落漏斗的深度越大. 當(dāng)抽出-回滲量分別為1、2.5、5、10 cm3/s時(shí),最大降深分別達(dá)到1.76、4.55、9.75、18.65 cm,分別為含水層厚度的2.5%、6.5%、13.9%、26.6%. 然而,抽出-回滲量對(duì)水位降落漏斗的范圍幾乎無(wú)影響,抽水影響半徑幾乎不隨抽出-回滲量的變化而變化.

    圖4 不同抽出-回滲量條件下模擬流場(chǎng)三維空間分布特征Fig.4 Characteristics of simulated three-dimensional spatial distribution of the flow field under different pumping-return seepage conditions

    表3 不同抽出-回滲量條件下模擬水位及抽水影響半徑Table 3 Simulated water level and pumping influence radius under different pumping-return seepage conditions

    2.2 實(shí)現(xiàn)水力控制時(shí)水位與抽注水量和含水層滲透性的關(guān)系

    已有研究[6]表明,實(shí)現(xiàn)水力控制時(shí)含水層不同位置處的水位h主要取決于抽出-回滲量Q(水動(dòng)力條件)和含水層滲透系數(shù)K(水文地質(zhì)條件). 基于抽出-回滲量和含水層滲透性設(shè)置不同情景,利用校準(zhǔn)后的數(shù)值模型模擬不同情景實(shí)現(xiàn)水力控制時(shí)地下水水位的分布特征,采用潛水完整井穩(wěn)定流公式-裘布依公式開(kāi)展數(shù)據(jù)分析擬合,分析實(shí)現(xiàn)水力控制時(shí)含水層不同位置水位與抽注水量和含水層滲透系數(shù)的關(guān)系,獲取單井抽出-回滲同步循環(huán)地下水水力控制技術(shù)關(guān)鍵參數(shù)之間的定量關(guān)系.

    選取《工程地質(zhì)手冊(cè)》提供的不同巖土類型(細(xì)砂、中砂、粗砂、礫砂)滲透系數(shù)的經(jīng)驗(yàn)值[35]設(shè)定4種水文地質(zhì)條件情景,選取不同抽出-回滲量(1、2.5、5、10 cm3/s)設(shè)定4種水動(dòng)力條件情景,共計(jì)16種組合情景(見(jiàn)表4).

    表4 不同情景水位模擬結(jié)果Table 4 Simulation results of water level under different scenarios

    基于裘布依公式,對(duì)16種組合情景模擬得到的抽水井處水位以及距抽水井井軸4、8、14、22、32 cm處的模擬水位進(jìn)行數(shù)據(jù)分析擬合,獲取Q/K與(h2-hw2)/(lgr-lgrw)之間的定量關(guān)系(見(jiàn)圖5)和擬合方程〔見(jiàn)式(5)〕. 由圖5可見(jiàn),Q/K與(h2-hw2)/(lgr-lgrw)之間呈線性關(guān)系,擬合系數(shù)(R2)達(dá)到0.99. 擬合后得到的線性方程斜率(α)為154.84~197.85,遠(yuǎn)大于式(2)中的斜率(1.366),主要是因?yàn)轸貌家拦降倪m用條件是補(bǔ)給項(xiàng)僅為側(cè)向滲流補(bǔ)給,不含上層回滲補(bǔ)給,而單井抽出-回滲同步循環(huán)地下水水力控制技術(shù)在抽水的同時(shí)從上層同步等量回滲補(bǔ)水,導(dǎo)致砂柱邊界地下水水位明顯高于中部地下水水位,水力梯度明顯增大.

    圖5 16種情景下模擬水位的數(shù)據(jù)分析擬合結(jié)果Fig.5 Simulated groundwater levels analyses and fitting results under 16 various scenarios

    擬合方程(5)的斜率區(qū)間為154.84~197.85,存在一定程度波動(dòng),為探究斜率波動(dòng)原因,基于表4設(shè)置的不同情景和現(xiàn)有模擬結(jié)果分析,利用校準(zhǔn)后的數(shù)值模型模擬不同回滲半徑下實(shí)現(xiàn)水力控制時(shí)擬合方程(5)的斜率變化. 目前回滲半徑為32 cm,進(jìn)一步考慮回滲半徑為15、20、25、30、35、40、45、50 cm 8種情況,在此基礎(chǔ)上對(duì)16種組合情景模擬得到的抽水井處水位以及距井軸4、8、14、22、32 cm處的模擬水位進(jìn)行數(shù)據(jù)分析擬合,對(duì)式(5)中的斜率進(jìn)行校正,結(jié)果如表5所示. 由表5可見(jiàn),回滲半徑與斜率之間存在一定關(guān)系,且擬合度比較高,R2均達(dá)到0.99. 隨著回滲半徑增加,擬合后得到的線性方程斜率在一定范圍內(nèi)逐漸增大,且擬合方程呈明顯的線性正相關(guān)關(guān)系;當(dāng)回滲半徑接近邊界半徑時(shí)斜率范圍逐漸減小,其中回滲半徑等于邊界半徑時(shí),斜率范圍為178.18~183.9,而回滲半徑為15 cm時(shí),斜率范圍為164.99~206.34,范圍較大;此外,當(dāng)回滲半徑固定時(shí),隨著監(jiān)測(cè)點(diǎn)位距離的增加,方程斜率逐漸增加. 主要原因是對(duì)于不變的水文地質(zhì)條件,當(dāng)流量確定后,回滲半徑大小不同,直接影響回滲范圍內(nèi)地下水水位的變化,因此合理地調(diào)整回滲范圍大小,有利于實(shí)現(xiàn)水力控制.

    表5 16種情景下的方程斜率Table 5 Equation of the slope of equations under 16 various scenarios

    3 討論

    通過(guò)單井-抽出回滲同步循環(huán)地下水水力控制技術(shù)的數(shù)值模擬結(jié)果,使用裘布依公式進(jìn)行了擬合,對(duì)該問(wèn)題中水動(dòng)力參數(shù)、水文地質(zhì)參數(shù)以及回滲范圍進(jìn)行了初步分析,得到了參數(shù)之間的對(duì)應(yīng)關(guān)系,當(dāng)回滲量取值為零時(shí),方程中各參數(shù)之間變?yōu)榻?jīng)典裘布依公式中的參數(shù)關(guān)系. 對(duì)不同巖土類型均質(zhì)含水層的模擬旨在概化不同類型滲透性場(chǎng)地中的參數(shù)關(guān)系.

    在實(shí)際應(yīng)用上,對(duì)場(chǎng)地回滲過(guò)程要求為均勻穩(wěn)定入滲,然而自然場(chǎng)地一般為非均質(zhì)性,入滲一般是不均勻、不穩(wěn)定的,對(duì)該問(wèn)題的較好解決方式仍是使用數(shù)值模擬,通過(guò)不同尺度下大量水動(dòng)力條件和非均質(zhì)的水文地質(zhì)條件的組合,同時(shí)調(diào)整抽出-回滲量和回滲范圍進(jìn)行模擬研究,以此得出參數(shù)之間的進(jìn)一步關(guān)系,使之更貼近于實(shí)際應(yīng)用條件.

    4 結(jié)論

    a) 校準(zhǔn)后的數(shù)值模型模擬結(jié)果與物理模型實(shí)際觀測(cè)結(jié)果相吻合,表明數(shù)值模型可以較好地刻畫(huà)實(shí)現(xiàn)水力控制時(shí)物理模型中的實(shí)際地下水流場(chǎng),可利用數(shù)值模型進(jìn)行單井抽出-回滲同步循環(huán)地下水水力控制技術(shù)研究.

    b) 抽出-回滲量的變化對(duì)水位降落漏斗的范圍幾乎無(wú)影響,而對(duì)水位降落漏斗的深度影響較大,抽出-回滲量越大,降落漏斗越深.

    c) 基于潛水完整井穩(wěn)定流公式-裘布依公式,對(duì)不同抽出-回滲量(水動(dòng)力條件)和含水層滲透性(水文地質(zhì)條件)情景模擬水位的擬合,獲取了實(shí)現(xiàn)水力控制時(shí)含水層不同位置水位與抽出-回滲量和含水層滲透系數(shù)之間的定量關(guān)系,可以為該技術(shù)在實(shí)際場(chǎng)地的應(yīng)用提供重要參考.

    d) 實(shí)現(xiàn)水力控制時(shí)含水層不同位置水位與抽出-回滲量和含水層滲透系數(shù)之間的定量關(guān)系受到回滲半徑的影響,對(duì)于特定水文地質(zhì)條件的含水層,合理調(diào)整抽出-回滲量以及回滲范圍大小,有利于實(shí)現(xiàn)水力控制.

    猜你喜歡
    滲透系數(shù)水力含水層
    水力全開(kāi)
    全球多個(gè)含水層里的水正快速流失
    基于Origin的滲透系數(shù)衰減方程在地?zé)崴毓嘀械膽?yīng)用
    多孔材料水滲透系數(shù)預(yù)測(cè)的隨機(jī)行走法
    輸水渠防滲墻及基巖滲透系數(shù)敏感性分析
    球墨鑄鐵管的水力計(jì)算
    美國(guó)西部奧加拉拉含水層水位下降原因初探
    河北平原新近系熱儲(chǔ)層滲透系數(shù)規(guī)律性分析
    水力噴射壓裂中環(huán)空水力封隔全尺寸實(shí)驗(yàn)
    全球地下含水層下降驚人:要被抽干了
    听说在线观看完整版免费高清| 欧美高清成人免费视频www| 亚洲经典国产精华液单 | 男女之事视频高清在线观看| 国产精品亚洲av一区麻豆| 国产欧美日韩精品一区二区| 两个人的视频大全免费| 十八禁网站免费在线| 欧美xxxx黑人xx丫x性爽| 国产单亲对白刺激| 国产精品不卡视频一区二区 | 免费在线观看成人毛片| 人人妻人人看人人澡| 精品国产亚洲在线| 国产欧美日韩精品亚洲av| 欧美成人a在线观看| 成年女人毛片免费观看观看9| 久久久久久久久中文| 日韩欧美 国产精品| 精品乱码久久久久久99久播| 日韩高清综合在线| 亚洲av成人精品一区久久| 国产精品1区2区在线观看.| 91av网一区二区| 亚洲性夜色夜夜综合| 亚洲成人久久爱视频| 高潮久久久久久久久久久不卡| 午夜两性在线视频| 精品午夜福利在线看| 18禁裸乳无遮挡免费网站照片| 久久精品国产亚洲av天美| 成人高潮视频无遮挡免费网站| 亚洲精品一区av在线观看| 国产精品日韩av在线免费观看| 国产真实乱freesex| 亚洲在线观看片| 久久精品91蜜桃| 香蕉av资源在线| 亚洲精华国产精华精| 精品午夜福利视频在线观看一区| 人妻制服诱惑在线中文字幕| 久久精品久久久久久噜噜老黄 | 蜜桃亚洲精品一区二区三区| 赤兔流量卡办理| 99热这里只有是精品50| 草草在线视频免费看| 久久久成人免费电影| 97人妻精品一区二区三区麻豆| 久久精品91蜜桃| 久久久久免费精品人妻一区二区| 日日夜夜操网爽| 国产精品国产高清国产av| eeuss影院久久| 69人妻影院| 性欧美人与动物交配| 91狼人影院| 免费看a级黄色片| 亚洲最大成人av| 51国产日韩欧美| 国产精品野战在线观看| 亚洲精品一区av在线观看| 国产精品精品国产色婷婷| 国产人妻一区二区三区在| 精品国内亚洲2022精品成人| 999久久久精品免费观看国产| 美女高潮喷水抽搐中文字幕| 欧美日韩福利视频一区二区| 真人做人爱边吃奶动态| 村上凉子中文字幕在线| 69av精品久久久久久| 成熟少妇高潮喷水视频| 黄色日韩在线| 国产精品一区二区三区四区久久| 午夜激情福利司机影院| 亚洲午夜理论影院| 亚洲国产精品999在线| 两个人视频免费观看高清| 极品教师在线视频| 成人国产综合亚洲| 亚洲,欧美精品.| 一本一本综合久久| 欧美日本视频| 久久精品夜夜夜夜夜久久蜜豆| 久久6这里有精品| 老司机福利观看| 国产毛片a区久久久久| 午夜a级毛片| 一边摸一边抽搐一进一小说| 国产综合懂色| 中文字幕av成人在线电影| 国产一区二区激情短视频| 狂野欧美白嫩少妇大欣赏| 午夜免费男女啪啪视频观看 | 99久久99久久久精品蜜桃| 日韩欧美国产一区二区入口| 丝袜美腿在线中文| 国产精品一区二区免费欧美| 成人国产一区最新在线观看| 岛国在线免费视频观看| 亚洲人成网站高清观看| 人人妻,人人澡人人爽秒播| 偷拍熟女少妇极品色| 老司机午夜福利在线观看视频| 欧美精品啪啪一区二区三区| 美女高潮喷水抽搐中文字幕| 欧美色欧美亚洲另类二区| 国产精品亚洲美女久久久| 国产野战对白在线观看| 久久天躁狠狠躁夜夜2o2o| 久久国产精品人妻蜜桃| 男人舔奶头视频| 一级毛片久久久久久久久女| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | av在线天堂中文字幕| 色视频www国产| 欧美激情久久久久久爽电影| 亚洲熟妇中文字幕五十中出| 国产亚洲欧美98| 欧美激情久久久久久爽电影| 色吧在线观看| 国产美女午夜福利| 又黄又爽又免费观看的视频| 亚洲熟妇中文字幕五十中出| xxxwww97欧美| 亚洲色图av天堂| 久久中文看片网| 国产淫片久久久久久久久 | 精品午夜福利在线看| 国产成人aa在线观看| 国产亚洲精品综合一区在线观看| 丁香欧美五月| 欧美成狂野欧美在线观看| 美女被艹到高潮喷水动态| 人人妻人人澡欧美一区二区| 日韩欧美精品免费久久 | 午夜亚洲福利在线播放| 国产在视频线在精品| 精品一区二区三区人妻视频| 一级毛片久久久久久久久女| 久久伊人香网站| 国产一级毛片七仙女欲春2| 久久久久久国产a免费观看| 夜夜夜夜夜久久久久| 99久久九九国产精品国产免费| 欧美午夜高清在线| 亚洲欧美激情综合另类| 欧美不卡视频在线免费观看| 哪里可以看免费的av片| 欧美另类亚洲清纯唯美| 久久久久久九九精品二区国产| 国产精品自产拍在线观看55亚洲| 国产v大片淫在线免费观看| 在线观看免费视频日本深夜| 90打野战视频偷拍视频| 国产精品永久免费网站| 国产精品影院久久| 18禁裸乳无遮挡免费网站照片| 成人欧美大片| 午夜a级毛片| 国产亚洲精品久久久com| 丁香六月欧美| 婷婷六月久久综合丁香| 桃色一区二区三区在线观看| 久久精品91蜜桃| 淫秽高清视频在线观看| 精品久久久久久久久久久久久| av在线蜜桃| 黄色视频,在线免费观看| 人妻久久中文字幕网| 午夜两性在线视频| 亚洲av第一区精品v没综合| 最近最新中文字幕大全电影3| 我的老师免费观看完整版| 美女xxoo啪啪120秒动态图 | 超碰av人人做人人爽久久| 国产精品久久电影中文字幕| 此物有八面人人有两片| 午夜免费激情av| 波多野结衣巨乳人妻| 亚洲成人中文字幕在线播放| 亚洲人与动物交配视频| 小说图片视频综合网站| 美女xxoo啪啪120秒动态图 | 亚洲七黄色美女视频| 国产又黄又爽又无遮挡在线| 国产伦在线观看视频一区| 精华霜和精华液先用哪个| 亚洲黑人精品在线| 久久99热6这里只有精品| 丰满人妻熟妇乱又伦精品不卡| 夜夜躁狠狠躁天天躁| 亚洲,欧美,日韩| 伦理电影大哥的女人| 亚洲 欧美 日韩 在线 免费| 久久精品久久久久久噜噜老黄 | 久久亚洲精品不卡| 亚洲精品影视一区二区三区av| 亚洲avbb在线观看| 亚洲三级黄色毛片| 偷拍熟女少妇极品色| 夜夜爽天天搞| 99热这里只有是精品在线观看 | 国产精品98久久久久久宅男小说| 日韩国内少妇激情av| 成人特级黄色片久久久久久久| 国产真实伦视频高清在线观看 | 少妇人妻一区二区三区视频| 变态另类丝袜制服| 两个人视频免费观看高清| 欧美黄色淫秽网站| 日本成人三级电影网站| 久久精品国产亚洲av香蕉五月| 国产一区二区在线av高清观看| 久久久色成人| 窝窝影院91人妻| 国产黄片美女视频| 日韩成人在线观看一区二区三区| 毛片一级片免费看久久久久 | 中出人妻视频一区二区| 又黄又爽又刺激的免费视频.| 波野结衣二区三区在线| 一区福利在线观看| 色播亚洲综合网| 日韩高清综合在线| 深夜a级毛片| 无遮挡黄片免费观看| 国产精品久久久久久精品电影| 国产伦人伦偷精品视频| 狂野欧美白嫩少妇大欣赏| 真人做人爱边吃奶动态| 亚洲av成人不卡在线观看播放网| 国产一区二区三区在线臀色熟女| 国产日本99.免费观看| 搡老熟女国产l中国老女人| 欧美日本亚洲视频在线播放| 国产精品女同一区二区软件 | 老司机午夜十八禁免费视频| 搡老熟女国产l中国老女人| 亚洲欧美日韩无卡精品| 91在线精品国自产拍蜜月| 午夜福利成人在线免费观看| 看十八女毛片水多多多| 久久久成人免费电影| 国产亚洲av嫩草精品影院| 久久久精品大字幕| 一区二区三区高清视频在线| 国产高清视频在线观看网站| bbb黄色大片| 尤物成人国产欧美一区二区三区| 亚洲 欧美 日韩 在线 免费| 国产国拍精品亚洲av在线观看| 能在线免费观看的黄片| 欧美日韩国产亚洲二区| 国产一区二区在线观看日韩| 久久性视频一级片| 波多野结衣巨乳人妻| 亚洲成人中文字幕在线播放| 亚洲五月天丁香| 九色国产91popny在线| 久久香蕉精品热| 我的女老师完整版在线观看| 亚洲av免费高清在线观看| 内射极品少妇av片p| 免费高清视频大片| 日韩欧美国产在线观看| 男女做爰动态图高潮gif福利片| 狂野欧美白嫩少妇大欣赏| 亚洲人成网站在线播放欧美日韩| 亚洲国产精品成人综合色| 在线播放无遮挡| 欧美精品国产亚洲| 午夜久久久久精精品| 日韩欧美精品v在线| 麻豆国产av国片精品| 在线观看66精品国产| 最近视频中文字幕2019在线8| 国产 一区 欧美 日韩| 亚洲国产精品sss在线观看| 亚洲人成电影免费在线| 国产精品一区二区性色av| 亚洲精华国产精华精| 内射极品少妇av片p| 国产精品日韩av在线免费观看| 国产综合懂色| 欧美午夜高清在线| 直男gayav资源| 成人三级黄色视频| 成年免费大片在线观看| 精品人妻偷拍中文字幕| 最近在线观看免费完整版| 亚洲专区国产一区二区| 欧美黑人巨大hd| 搞女人的毛片| 欧美一区二区亚洲| 亚洲国产日韩欧美精品在线观看| 90打野战视频偷拍视频| 亚洲中文日韩欧美视频| 国产一区二区亚洲精品在线观看| 日韩大尺度精品在线看网址| 天堂动漫精品| 亚洲美女搞黄在线观看 | 国产高潮美女av| 亚洲欧美清纯卡通| 全区人妻精品视频| 高清毛片免费观看视频网站| 熟女电影av网| a在线观看视频网站| 又黄又爽又免费观看的视频| 女人十人毛片免费观看3o分钟| 久久久久免费精品人妻一区二区| 最近视频中文字幕2019在线8| 国产欧美日韩精品一区二区| 女生性感内裤真人,穿戴方法视频| 久久久久久久亚洲中文字幕 | 黄色丝袜av网址大全| 国产精品久久久久久久电影| 国产精品一及| a在线观看视频网站| 日本一本二区三区精品| 国产亚洲精品久久久久久毛片| av欧美777| 国产精品久久视频播放| 精品一区二区三区av网在线观看| av天堂在线播放| 欧美激情在线99| 日本一本二区三区精品| 国产精品影院久久| 欧美xxxx黑人xx丫x性爽| 成人一区二区视频在线观看| 深夜a级毛片| 精品欧美国产一区二区三| 最近最新免费中文字幕在线| 噜噜噜噜噜久久久久久91| 午夜久久久久精精品| 亚洲精品一区av在线观看| 我要搜黄色片| 日日摸夜夜添夜夜添小说| 午夜福利在线在线| 午夜免费激情av| 国产精品乱码一区二三区的特点| 国产精品日韩av在线免费观看| 精品欧美国产一区二区三| 中文字幕免费在线视频6| 亚洲,欧美,日韩| 国产欧美日韩精品亚洲av| 精品久久久久久成人av| 欧美又色又爽又黄视频| 男人和女人高潮做爰伦理| 老熟妇乱子伦视频在线观看| 亚洲最大成人av| 精品久久久久久久久久久久久| 日韩欧美精品v在线| 18禁黄网站禁片午夜丰满| 欧美区成人在线视频| 亚洲av电影在线进入| 级片在线观看| 亚洲美女黄片视频| 亚洲av免费高清在线观看| 国产探花极品一区二区| 黄片小视频在线播放| 欧美3d第一页| 亚洲精品在线美女| 国产黄色小视频在线观看| av在线老鸭窝| 18禁黄网站禁片免费观看直播| 真人一进一出gif抽搐免费| 日韩国内少妇激情av| 99在线人妻在线中文字幕| 人妻久久中文字幕网| av福利片在线观看| 亚洲熟妇熟女久久| 国产中年淑女户外野战色| 成人高潮视频无遮挡免费网站| 一进一出好大好爽视频| 欧美色欧美亚洲另类二区| 精品国产亚洲在线| а√天堂www在线а√下载| 国产探花在线观看一区二区| avwww免费| 日韩有码中文字幕| 亚洲av二区三区四区| 欧美一区二区亚洲| 淫妇啪啪啪对白视频| eeuss影院久久| 中文字幕精品亚洲无线码一区| 在线观看免费视频日本深夜| 国产黄色小视频在线观看| 欧美潮喷喷水| 亚洲av.av天堂| 99久久精品国产亚洲精品| 在线天堂最新版资源| 青草久久国产| 亚洲 欧美 日韩 在线 免费| 国产人妻一区二区三区在| 亚洲 国产 在线| 久久久久性生活片| 日韩中字成人| 久久久久精品国产欧美久久久| 国产私拍福利视频在线观看| 免费人成在线观看视频色| 久久久久久国产a免费观看| 午夜免费男女啪啪视频观看 | 免费看a级黄色片| 国产成人欧美在线观看| av欧美777| 尤物成人国产欧美一区二区三区| 精品午夜福利在线看| 成人一区二区视频在线观看| 听说在线观看完整版免费高清| 国产免费男女视频| 亚洲国产精品成人综合色| 亚洲专区国产一区二区| 国产精品嫩草影院av在线观看 | 久久精品国产自在天天线| 精品久久久久久成人av| 精品乱码久久久久久99久播| 国产激情偷乱视频一区二区| 欧美性猛交╳xxx乱大交人| 国产成+人综合+亚洲专区| 欧美zozozo另类| 变态另类成人亚洲欧美熟女| 亚洲最大成人手机在线| 波多野结衣高清无吗| 国产69精品久久久久777片| 精华霜和精华液先用哪个| 欧洲精品卡2卡3卡4卡5卡区| 亚洲av美国av| 久久久久久久久中文| 高清日韩中文字幕在线| 久久久精品大字幕| 国产精品爽爽va在线观看网站| 韩国av一区二区三区四区| 亚洲七黄色美女视频| 波野结衣二区三区在线| 最近在线观看免费完整版| 此物有八面人人有两片| 精品日产1卡2卡| 欧美午夜高清在线| 日日摸夜夜添夜夜添小说| av在线天堂中文字幕| 欧美激情久久久久久爽电影| 美女cb高潮喷水在线观看| 三级国产精品欧美在线观看| 99热这里只有是精品在线观看 | 精品久久久久久久久久免费视频| 国产三级黄色录像| 内射极品少妇av片p| 亚洲精品影视一区二区三区av| 动漫黄色视频在线观看| 午夜激情福利司机影院| 午夜福利欧美成人| 精品久久久久久久久久免费视频| or卡值多少钱| 国产毛片a区久久久久| 俺也久久电影网| 国产精品一区二区性色av| 久久久久性生活片| 国产精品av视频在线免费观看| 国产在线精品亚洲第一网站| 久久久精品欧美日韩精品| 黄色一级大片看看| 青草久久国产| 免费看a级黄色片| 内地一区二区视频在线| 国产欧美日韩精品亚洲av| 欧美日本视频| 精品久久久久久久末码| 久久久久久九九精品二区国产| 国产精品亚洲av一区麻豆| 黄色女人牲交| 国产探花在线观看一区二区| 欧美日韩乱码在线| 熟女电影av网| 免费在线观看亚洲国产| 久久精品国产亚洲av香蕉五月| 高潮久久久久久久久久久不卡| 国产精品自产拍在线观看55亚洲| 国语自产精品视频在线第100页| 久久性视频一级片| 最好的美女福利视频网| 久久伊人香网站| 变态另类丝袜制服| 伊人久久精品亚洲午夜| 十八禁人妻一区二区| 国产精品久久久久久久电影| 俺也久久电影网| 国产精品久久久久久久电影| 高清毛片免费观看视频网站| 国产精品一区二区三区四区免费观看 | 欧美xxxx性猛交bbbb| 亚洲精品在线美女| 精品欧美国产一区二区三| 成人av一区二区三区在线看| 99精品久久久久人妻精品| 老熟妇仑乱视频hdxx| av黄色大香蕉| 精品一区二区免费观看| 99久久成人亚洲精品观看| 婷婷色综合大香蕉| 夜夜夜夜夜久久久久| 亚洲精品一区av在线观看| xxxwww97欧美| 人人妻人人澡欧美一区二区| 亚洲欧美激情综合另类| 国产男靠女视频免费网站| 国产69精品久久久久777片| 国产三级在线视频| 99热这里只有是精品在线观看 | 国内精品久久久久久久电影| 久久精品人妻少妇| 国产麻豆成人av免费视频| 成人国产一区最新在线观看| 狂野欧美白嫩少妇大欣赏| 嫩草影视91久久| 一级黄片播放器| 亚洲精品亚洲一区二区| 在线看三级毛片| 18禁黄网站禁片免费观看直播| 久久国产精品人妻蜜桃| 日韩欧美免费精品| 精品国产亚洲在线| 国产成人欧美在线观看| 国产高清视频在线观看网站| 亚洲va日本ⅴa欧美va伊人久久| 国产亚洲av嫩草精品影院| 国产三级在线视频| 真实男女啪啪啪动态图| 深爱激情五月婷婷| 欧美三级亚洲精品| 午夜激情福利司机影院| 久久久久久久久大av| netflix在线观看网站| 中亚洲国语对白在线视频| 精品人妻偷拍中文字幕| 日本在线视频免费播放| 露出奶头的视频| 一个人看的www免费观看视频| 18禁黄网站禁片午夜丰满| 热99在线观看视频| a级一级毛片免费在线观看| 99久久无色码亚洲精品果冻| 久久久久久九九精品二区国产| 亚洲人成网站在线播| 国产欧美日韩精品亚洲av| 成人高潮视频无遮挡免费网站| 久9热在线精品视频| 精品久久久久久久久亚洲 | 好男人电影高清在线观看| 别揉我奶头 嗯啊视频| 亚洲av日韩精品久久久久久密| h日本视频在线播放| 亚洲人与动物交配视频| 精品一区二区免费观看| 国产熟女xx| 禁无遮挡网站| 亚洲精品成人久久久久久| 看免费av毛片| 色播亚洲综合网| 男人舔奶头视频| 国产大屁股一区二区在线视频| 我要看日韩黄色一级片| 精品国产亚洲在线| 精品日产1卡2卡| 国产精品日韩av在线免费观看| 婷婷精品国产亚洲av| 久久久久性生活片| 美女高潮喷水抽搐中文字幕| 成人高潮视频无遮挡免费网站| 欧美日韩福利视频一区二区| 九九久久精品国产亚洲av麻豆| 日日摸夜夜添夜夜添小说| 国产精品久久视频播放| 国产野战对白在线观看| 免费观看精品视频网站| 嫩草影院入口| 国产野战对白在线观看| 亚洲三级黄色毛片| 性插视频无遮挡在线免费观看| 精品人妻1区二区| 精品久久久久久久久久免费视频| 亚洲av熟女| 此物有八面人人有两片| 国产高清三级在线| 成人午夜高清在线视频| 97碰自拍视频| 久久人妻av系列| 性色av乱码一区二区三区2| 老熟妇仑乱视频hdxx| 亚洲第一电影网av| 国内毛片毛片毛片毛片毛片| 12—13女人毛片做爰片一| 99热这里只有精品一区| 乱人视频在线观看| 三级国产精品欧美在线观看| 欧洲精品卡2卡3卡4卡5卡区| 18禁在线播放成人免费| 每晚都被弄得嗷嗷叫到高潮| 国产av一区在线观看免费| 欧美日韩乱码在线| 精品人妻偷拍中文字幕| 老女人水多毛片| 女人十人毛片免费观看3o分钟| 久久久久亚洲av毛片大全| 国产精品一区二区免费欧美| 国产精品久久久久久人妻精品电影| 好男人电影高清在线观看| 最近最新中文字幕大全电影3| 一区二区三区四区激情视频 | 级片在线观看| 狠狠狠狠99中文字幕| 大型黄色视频在线免费观看| 级片在线观看| 欧美zozozo另类| 欧美日本亚洲视频在线播放| 成人欧美大片|