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

    電法監(jiān)測城市淺層地下水污染可行性研究

    2020-06-10 01:41:14崔少北李曉斌楊雙安
    物探化探計算技術 2020年2期
    關鍵詞:曲線圖測線電導率

    崔少北, 李曉斌, 楊雙安

    (1.河南理工大學 資源環(huán)境學院,焦作 454000;2.中原經(jīng)濟煤層(頁巖)氣河南省協(xié)同創(chuàng)新中心,焦作 454000)

    0 引言

    2000 年以來,我國工業(yè)化快速發(fā)展,與此同時造成的地下水污染和破壞日益加劇,對生態(tài)環(huán)境也帶來了不利影響。據(jù)水利部2016年月發(fā)布的“地下水動態(tài)月報”顯示,在抽查的18個省市區(qū)地下水抽樣中,地下水水質(zhì)不達標的約占60%[1]。在全國范圍內(nèi),約有90%的城市地下水遭受到不同程度的污染[2]。城市淺層地下水水質(zhì)受各種污染源綜合作用的影響,包括污水處理廠排污、垃圾填埋場的瀝出液、工廠排污、排水管道滲漏、大氣沉降、建筑垃圾及游樂場排污等[3]。因此,能夠準確及時地對城市地下水污染情況進行有效地監(jiān)測,改善城市生態(tài)系統(tǒng)水環(huán)境質(zhì)量,對生態(tài)建設、工業(yè)發(fā)展及居民生活有重要意義。

    現(xiàn)有地下水監(jiān)測的方法主要采用定期在地下水監(jiān)測井中取樣,用水化學方法分析水樣中污染物成分和濃度,進而確定地下水的污染狀況。顯然,該方法費時、費力,采樣分析結果具有滯后性,也不能給出地下污染水體的邊界范圍。因此,尋找一種快速有效的方法,能對城市淺層地下水進行長期監(jiān)測已經(jīng)迫在眉睫。地球物理方法主要用于探測地球本體及近地空間的介質(zhì)結構、物質(zhì)組成、形成演化,在工程、環(huán)保、礦產(chǎn)等方面有非常廣泛的應用,其中電法勘探在地下水勘測方面具有顯著優(yōu)勢。在利用電阻率差異探測淺層地下水方面,劉宏偉等[4]基于水化學與地球物理法的萊州灣南岸海(咸)水入侵勘查,結合水化學法和地球物理探測方法,探索出EH4法和準三維高密度電法,較好地指示咸淡水體在空間上的分布,框定海(咸)入侵水界線的大致范圍[4];陳超啟等[5]高采用密度電法對北京市區(qū)周邊某廢棄垃圾填埋場垃圾滲濾液下滲情況進行勘查分析,取得較好效果。

    城市的各種污染主要造成地下水中的硝酸根離子濃度的增加[6],在一些工業(yè)發(fā)展水平較高的城市和地區(qū),地下水中致癌性重金屬嚴重超標[7-8]。硝酸根離子濃度的增加和重金屬滲透會使淺層地下水的電導率顯著提高,地下受污染水體的電學性能與未污染水體及周圍環(huán)境在電阻率等物理特性方面有顯著差異。地球物理電法的供電電源在地下空間中形成的穩(wěn)定電場,可以圈定出地下水低阻異常區(qū)的范圍,推測地下水體受污染情況,追蹤地下水污染源,為淺層地下水治理提供參考。根據(jù)這一特性,筆者針對城市淺層地下水水質(zhì)預警,提出通過在潛水層甚至更深的地層中布設足夠長的供電電極,長期監(jiān)測城市淺層水電阻率變化。由于長供電電極可以深埋地下,提供更多的電流進入深部地層,其地電場反映有更多深部信息[9]。因此,針對城市淺層地下水賦存條件,設計不同地質(zhì)物理模型,分別用點電極與長電極供電電源對模型激勵相同電壓,數(shù)值模擬地下水污染前后的電場分布特征,比較兩種電極在監(jiān)測城市淺層地下水污染狀況方面的差異,論證點電極和長電極電源在監(jiān)測預警城區(qū)淺層地下水體污染方面的可行性。

    1 理論與方法

    1.1 理論基礎

    常規(guī)直流電法采用不銹鋼或銅電極向地下提供直流電場,利用相同的電極測量地面兩點之間電位差,進一步反演獲得地下地質(zhì)體的電學特性。由于電極長度較短,地面電法勘探中近似為點電流源。然而,采用長電極供電時,其長度不能忽略。假設大地是一個無窮大的均質(zhì)半空間,地下水體為空間內(nèi)低阻異常區(qū),分別放置點電極和長電極兩種供電電源,分析其在均質(zhì)空間中電位的分布特征。

    在地表A(0,0,0)點處放置電流為I的點電極(如圖1),則電流密度J、介質(zhì)電導率σ與電位u滿足以下關系:

    J=-σ▽u

    (1)

    在地下的半空間中做任意閉合曲面Г,曲面包圍的區(qū)域為Ω。根據(jù)通量定律,電流I與電位u滿足微分方程:

    圖1 點電極空間電位示意圖Fig.1 Point electrode space potential diagram

    圖2 長電極空間電位示意圖Fig.2 The long electrode space potential diagram

    (2)

    可以采用混合邊界條件[10-14],即:

    (3)

    式中:a=cos(r,n)/r,r是點電源與邊界?!奚系狞c之間的距離,cos(r,n)是r與外法向n的夾角余弦。在地表處,沿地表法線方向沒有電流,即a=0,則式(3)轉化為黎曼(Neumann)邊界條件[13]。令▽·(σ▽u)=-2Iδ(A),構造泛函:

    (4)

    對式(4)作變分處理,點電極無限半空間區(qū)域的邊值問題與下列變分問題等價:

    (5)

    對長電極的泛函問題,因為長電極電阻率一般遠小于周圍地層電阻率,可視長電極為等勢體。將其分解為一系列的電流元dl,每個電流元可以近似為一個獨立的點電極(如圖2)。在柱坐標中,電流元dl在B(r,φ,d)點處的電位可由鏡像法得出。同理,長電極L在半空間中B點的電位與長電極2L在全空間中B點的電位相同。所以沿2L對電流元dl積分,可得長電極在半空間中的電位分布。

    由于電位關于Z軸軸對稱,不考慮角度φ,電流元dl在地下均質(zhì)半空間任意B點處的電位可由下式得出:

    (6)

    其中:J為長電極L的線電流密度。積分得出表達式如下:

    (7)

    圖2所示的半空間中,長電極L是一個線性系統(tǒng)[15],L的泛函可由式(5)積分給出:

    (8)

    當長電極L穿過不同介質(zhì)時,電極表面電流隨周圍介質(zhì)電導率變化而改變。電流密度分布的解析式可以由長電極表面電位近似導出,不同介質(zhì)的長電極近似為等勢體,由式(1)可得:

    (9)

    式中:Ji為電極電流密度;σ為長電極所在介質(zhì)電導率;i為長電極穿過介質(zhì)的序數(shù)。當長電極穿過的N個不同介質(zhì)時,有以下關系:

    (10)

    (11)

    在式(10)中:Li為長電極在第i區(qū)域的長度;σi為該區(qū)域的電導率。由以上兩式可得第k段長電極上的電導率:

    (12)

    1.2 有限元分析

    上世紀50年代,有限元方法在彈性力學領域大量使用,隨后發(fā)展為地球物理學中重要的正演數(shù)值模擬方法[16-18]。由于非結構網(wǎng)格沒有規(guī)則的拓撲關系,組成網(wǎng)格的單元之間相對獨立,可以構建任何復雜的幾何模型。

    采用多重網(wǎng)格方法將空間域剖分為正四面體空間離散單元,全空間中任一單元K的電導率設為常數(shù)σ,正四面體網(wǎng)格具有四個節(jié)點,使用線性插值[16],單元K中電位u由如下公式確定:

    (13)

    式中:ui為單元K上第i節(jié)點處電位的有限元近似解;φi為對應的基函數(shù)或形函數(shù)。將式(13)代入泛函式(8)可得其離散格式:

    (14)

    式中:

    u= (u1,u2,u3,u4)T

    φ= (φ1,φ2,φ3,φ4)T

    其中:NE為空間剖分的單元總數(shù);P為全局電源分布項,相對于點電極的電阻率有限元方法;單元矩陣A、B多了沿長電極的線積分。對式(14)取一階變分,并令其等于零,即可得到基于節(jié)點有限元方法的線性方程組:

    QU=P

    (15)

    式中:U代表所有的待求未知量,即所有單元節(jié)點處的電位;Q為全局系數(shù)矩陣,一般為對稱正定稀疏矩陣。各單元內(nèi)部電位均可利用單元網(wǎng)絡節(jié)點處的電位公式(13)插值得到。

    2 模型設計

    根據(jù)常見城市淺層地下水的賦存形態(tài),設計圓餅狀、條帶狀、球狀三種模型。圓餅狀模型模擬賦存于在地表有廣泛分布的第四系松散沉積物空隙中的潛水層,其形態(tài)受基巖表面特征與其上部自由水面控制。潛水層在近地表沿地面平行展布,層厚度遠小于延展半徑。條帶狀是常見的地下水賦存形態(tài),在層間錯動形成的破碎裂隙帶中或者巖溶形成的地下暗河中,地下水常呈條帶狀分布。在巖溶地區(qū)或一些構造復雜的地區(qū),地下賦水空間發(fā)育成洞室狀,長寬高數(shù)值相近,以球狀模型模擬此地下水形態(tài)。

    圖3 圓餅狀地下水模型IFig.3 Round cake groundwater model I

    圖4 圓餅狀地下水XZ軸切面圖Fig.4 XZ axis cutaway view of round cake groundwater

    圖5 條帶狀(線狀)地下水模型IIFig.5 Strip Groundwater Model II

    圖6 球狀(點狀)地下水模型IIIFig.6 Spherical groundwater model

    如圖3所示,圓餅狀模型I的基本形態(tài)。中心藍色區(qū)域為地下水污染區(qū),其電導率σ1取值位0.5 S/m,灰色區(qū)域為未污染潛水層,電導率σ2取值0.1 S/m,上層覆土壤和基底隔水層電導率σ3取值0.01 S/m。布置五條測線:SI-1、SI-2、SI-3、SI-4、SI-5。地表水平測線SI-1的起點坐標為E(20,0,0),終點為F(140,0,0),4條垂直測線距原點O的距離為:40 m、70 m、100 m、130 m,具體參數(shù)參見表1。

    條帶狀是常見的地下水賦存空間,圖5模擬地下暗河的地下水分布特性,中心藍色區(qū)域為地下水污染區(qū)。在模型II中,地下暗河中心坐標為(0,0,-30),點電極與長電極位于坐標原點O(0,0,0),長電極L穿過地下暗河。在模型上布置四條測線:SII-1、SII-2、SII-3、SII-4。測線SII-1平行地下暗河走向,端點坐標為(50,0,0)、(550,0,0)。測線SII-2垂直暗河走向,端點坐標為(250,200,0)、(250,-200,0)。兩條垂向測線SII-3、SII-4,地表出漏點坐標為(250,0,0)(250,50,0),具體參數(shù)參見表2。

    模型III模擬賦水空間發(fā)育成洞室狀的分布特性。在A(30,0,0)點放置兩種電極,同樣布置四條測線:SIII-1、SIII-2、SIII-3、SIII-4。水平測線SIII-1、SIII-2互相垂直,坐標為SIII-1:(10,0,0)、(-80,0,0),SIII-2:(0,80,0)、(0,-80,0)。垂向測線SIII-3、SIII-4地表出漏點坐標為(0,0,0)(-80,0,0),具體參數(shù)參見表3。

    3 數(shù)值計算分析

    3.1 圓餅狀地下水模型

    對圓餅狀地下水模型中點O點處分別放置點電極和長電極,施加200 V的直流電場,無窮遠處電位為零。經(jīng)過計算,可得地下水污染前后地表測線SI-1上電位變化特征(圖7),分析圖7可知,點電極附近地表電位數(shù)值較大,隨距離r的增加電位值急劇衰減,在距供電電極較遠處,電位曲線趨于平緩。

    表1 圓餅狀地下水模型參數(shù)值

    表2 條帶狀地下水模型參數(shù)值

    表3 球狀地下水模型參數(shù)值

    圖7 點電極地表測線SI-1電位曲線圖Fig.7 Point electrode surface survey line SI-1 potential curve

    圖8 點電極XZ軸面污染中心區(qū)電位等值線圖(有污染)Fig.8 Point electrode XZ axial pollution center area potential contour map(with pollution)

    被污染水體對應的地表電位數(shù)值較污染前變小、電位下降加快。這是由于地下水體被污染后,電導率增大,吸收更多電流,地表電位隨之降低。

    在模型I中XZ軸面上,作地下水體被污染后的等值線圖(圖8)。圖上由外及內(nèi)分別作:0.2 V、0.25 V、0.3 V、1 V,五個電位等值線。以0.2 V等值線所在等值面旋轉3600可得模型IV(圖9)。

    同理可得XZ軸切面上兩種供電模式在地下水被污染前后的電位分布(圖10、圖11、圖12)。在點電極供電條件下,電位0.23 V等值線分布于污染水體周圍空間,受水體電導率變化影響較小,空間電位的擾動主要出現(xiàn)在污染區(qū)內(nèi)部。地下水被污染后,0.3 V電位等值線分布區(qū)域減小,污染區(qū)內(nèi)電位值處在0.23 V~0.3 V之間,污染水體接更近等勢體。在長電極供電條件下,XZ軸切面上的電位值遠大于點電極。地下水被污染后,電場擾動范圍擴大,電導率增高使得高電位分布區(qū)明顯擴大。沿圖中五條測線繪制電位曲線,作地下水被污染前后各測線上電位變化曲線(圖13~圖16)。

    圖9 點電極XZ軸切面電位分布圖(有污染)Fig.9 Point electrode XZ axis cut surface potential distribution map (Pollution)

    圖10 長電極XZ軸切面電位分布圖(有污染)Fig.10 Long electrode XZ axis cut surface potential distribution map (Pollution)

    圖11 點電極XZ軸切面電位分布圖(無污染)Fig.11 Point electrode XZ axis cut surface potential distribution map (non-pollution)

    圖12 長電極XZ軸切面電位分布圖(無污染)Fig.12 Long electrode XZ axis cut surface potential distribution map (non-pollution)

    圖13 點電極地表測線SI-1電位曲線圖Fig.13 Point electrode surface measurement line SI-1 potential curve

    圖14 長電極地表測線SI-1電位曲線圖Fig.14 Long electrode surface measurement line SI-1 potential curve

    圖15 點電極垂直測線電位曲線圖Fig.15 Point electrode vertical line potential curve

    圖16 長電極垂直測線電位曲Fig.16 Long electrode vertical line potential curve

    圖17 點電極地表測線電位曲線圖Fig.17 Point electrode surface line potential curve

    圖18 長電極地表測線電位曲線圖Fig.18 Long electrode surface line potential curve

    圖19 點電極垂直測線電位曲線圖Fig.19 Point electrode vertical line potential curve

    圖20 長電極垂直測線電位曲線圖Fig.20 Long electrode vertical line potential curve

    在圖13中,點電極供電條件下,污染區(qū)界面以內(nèi)的電位減小,電位曲線下降速率減緩。將80 m~120 m段電位曲線圖放大,兩條電位曲線在污染區(qū)邊界以外基本重合,地下水被污染后地表電位曲線在污染區(qū)邊界處下降速率變大。這是由于在污染區(qū)邊界處,地下水電導率由σ1變?yōu)棣?,電位下降速率變大,電位曲線出現(xiàn)輕微彎折。在垂直方向上,污染區(qū)邊界(測線SI-4)以內(nèi)電位數(shù)值明顯減小,隨深度的增加電位下降速率減緩,潛水層中電位曲線與X軸近似平行。在污染區(qū)邊界以外,SI-4、SI-5兩條測線上的電位分布變化不大,這說明污染水體對遠處地下空間電位的影響較小,電導率變化引起的空間電位異常主要分布于被污染水體上部地表以及污染水體所在空間層位。

    相對而言,長電極供電條件下,水體污染后地表電位數(shù)值變大,電位下降速率減緩,如圖14所示,在污染區(qū)邊界處電位下降速率有明顯改變。在圖16中,潛水層下界面以上空間電位幾乎不隨深度發(fā)生變化,地下水體被污染后的電位遠大于污染前。污染水體對電位分布的擾動范圍擴大,增大的電位區(qū)域延伸至污染區(qū)邊界處(SI-4)以外。

    圖21 點電極地表測線電位曲線圖Fig.21 Point electrode surface line potential curve

    圖22 長電極地表測線電位曲線圖Fig.22 Long electrode surface line potential curve

    圖23 點電極垂直測線電位曲線圖Fig.23 Point electrode vertical line potential curve

    圖24 長電極垂直測線電位曲線圖Fig.24 Long electrode vertical line potential curve

    綜合以上分析,對于圓餅狀地下含水層,點電極電場受地下水污染狀況影響較小,電位變化不明顯。長電極穿過地下含水層,地下水體被污染后,電導率變大,電位下降速率減緩,對應的地表電位值大幅增加。所以,長電極在尋找地下水以及地下水體污染監(jiān)測方面具有顯著優(yōu)勢。

    3.2 條帶狀地下水模型

    對模型II模擬計算,繪制電位曲線圖(圖17~圖20)。對比圖17~圖20分析,長電極供電條件下電位數(shù)值遠大于點電極。兩種電極地表測線SII-1上的電位隨r(測點與原點的距離)的增大逐步降低。地下水體被污染后,點電極附近地表電位變小,電位下降速率減緩,隨r的增大電位又逐漸高于污染前電位,在賦水區(qū)邊界處,電位曲線下降速率加快,邊界以外兩條電位曲線趨向于重合。測線SII-2關于SII-1對稱,地下電位關于含水空間走向對稱分布。在垂直測線上,穿過含水層的測線SII-3上的電位曲線,在含水層下界面處曲線出現(xiàn)拐點,下界面以上空間中電位變化較小。在右下角放大圖中可以發(fā)現(xiàn),當測線深度超過下界面時,電位數(shù)值快速下降。

    3.3 球狀

    在模型III中,點電極和長電極的位置偏離球心,長電極插入含水層。對模型III進行模擬計算,可以得到球狀儲水空間地下水電位分布特征??梢园l(fā)現(xiàn),長電極電位值大于點電極。點電極供電條件下,測線SIII-1上的兩條電位曲線,差異微小。SIII-2測線上的電位在球心上部地表位置出現(xiàn)極值,污染前的電位大于污染后,這是由于地下水體對電場的擾動在球心處最大。在長電極供電條件下,地表兩條測線上電位曲線接近重合。測線SIII-1與長電極共面,隨著測點與電極距離的增大,電位曲線逐漸降低。SIII-2測線垂直SIII-1,電位曲線呈對稱狀,在原點處取得極值。觀察圖23,穿過模型III的測線SIII-3上,電位曲線有兩個拐點,球型模型所在空間的電位曲線下降平緩,上下兩端測線上電位變化速率較大。這是由于球形儲水空間中水體電導率較大,電位變化速率較小。測線SIII-4在球心等深度線上的電位曲線出現(xiàn)拐點,這是由于球形儲水空間體積較小,其對空間電場的擾動出現(xiàn)在球心以上空間。

    4 結論

    筆者結合常見的地下水形態(tài),設計了三種地下水模型。分別用長電極和點電極提供穩(wěn)定電壓,分析地表電位的差異,比較兩種電極在監(jiān)測地下水污染情況方面的優(yōu)劣。通過分析三個模型電位分布特征,得出以下結論:

    1)長電極和點電極通過提供穩(wěn)定的電場,都可以監(jiān)測地下水污染情況。由于長電極可以將更多電流傳入深部地層,其電場帶有更多信息,對地下水電導率變化更加敏感。長電極形成的電場強度更大,電位差異更加明顯,具有較高的信噪比。點電極供電電流受地表覆土阻隔,對地下水電導率變化不敏感。

    2)對于不同類型的地下水模型,其電場分布特征不同。地下水體被污染后,水體對電場的擾動受賦存形態(tài)影響。在對地下水污染情況進行監(jiān)測時,合理布置供電電極與測線的位置,可以有效獲取地下水污染情況。

    3)根據(jù)水文地質(zhì)特征具體分析,考慮到地下水深度和滲流速度的影響,結合長電極和點電極電源的電法設計城市淺層地下水污染狀況監(jiān)測預計系統(tǒng),可以長期有效的監(jiān)測地下水污染狀況。

    猜你喜歡
    曲線圖測線電導率
    極地海洋多波束測量測線布設系統(tǒng)設計及實現(xiàn)
    基于動態(tài)規(guī)劃的多波束測線布設模型
    秦皇島煤價周曲線圖
    秦皇島煤價周曲線圖
    秦皇島煤價周曲線圖
    秦皇島煤價周曲線圖
    基于比較測量法的冷卻循環(huán)水系統(tǒng)電導率檢測儀研究
    低溫脅迫葡萄新梢電導率和LT50值的研究
    高電導率改性聚苯胺的合成新工藝
    技術與教育(2014年2期)2014-04-18 09:21:33
    隧洞中雷達探測地質(zhì)構造的測線布置與三維地質(zhì)解譯
    在线观看一区二区三区激情| 少妇猛男粗大的猛烈进出视频| 一区二区三区国产精品乱码| 亚洲专区中文字幕在线| 日韩欧美三级三区| 久久久久国内视频| 国产男女内射视频| 成年版毛片免费区| 亚洲欧美色中文字幕在线| 成年人免费黄色播放视频| 亚洲国产欧美一区二区综合| 国产精品欧美亚洲77777| 我的亚洲天堂| 蜜桃国产av成人99| av视频免费观看在线观看| 久久人妻福利社区极品人妻图片| 亚洲人成伊人成综合网2020| 99久久99久久久精品蜜桃| 欧美日韩黄片免| 91麻豆精品激情在线观看国产 | a级毛片黄视频| 精品亚洲乱码少妇综合久久| 国产精品av久久久久免费| 国产国语露脸激情在线看| 制服人妻中文乱码| av不卡在线播放| 午夜视频精品福利| 十八禁网站免费在线| 99国产精品一区二区蜜桃av | 免费av中文字幕在线| 老司机亚洲免费影院| 免费黄频网站在线观看国产| 黄频高清免费视频| 久久久久久亚洲精品国产蜜桃av| 免费日韩欧美在线观看| 久久精品熟女亚洲av麻豆精品| 纯流量卡能插随身wifi吗| 欧美午夜高清在线| 嫁个100分男人电影在线观看| 多毛熟女@视频| 欧美国产精品va在线观看不卡| 亚洲成人免费av在线播放| 欧美激情极品国产一区二区三区| 久久午夜亚洲精品久久| 人人妻人人添人人爽欧美一区卜| 欧美亚洲日本最大视频资源| 国产亚洲欧美在线一区二区| 国产国语露脸激情在线看| 91精品国产国语对白视频| 动漫黄色视频在线观看| 黄色视频不卡| aaaaa片日本免费| 9191精品国产免费久久| 99久久人妻综合| 国产男女超爽视频在线观看| 一本色道久久久久久精品综合| 精品卡一卡二卡四卡免费| 三上悠亚av全集在线观看| 捣出白浆h1v1| 亚洲人成电影免费在线| netflix在线观看网站| 日韩制服丝袜自拍偷拍| 极品教师在线免费播放| 色综合欧美亚洲国产小说| 日韩免费高清中文字幕av| 十八禁高潮呻吟视频| 视频区欧美日本亚洲| 人人妻人人澡人人看| 免费在线观看日本一区| 亚洲av成人不卡在线观看播放网| 国产精品一区二区在线观看99| 国产精品成人在线| 热99re8久久精品国产| 精品久久蜜臀av无| 国产精品.久久久| 欧美精品一区二区大全| 18禁黄网站禁片午夜丰满| 制服人妻中文乱码| av不卡在线播放| 97在线人人人人妻| 91精品三级在线观看| 狠狠婷婷综合久久久久久88av| 国产黄色免费在线视频| 成年女人毛片免费观看观看9 | 亚洲欧美一区二区三区黑人| 国产一区二区三区综合在线观看| 亚洲av美国av| 午夜福利在线观看吧| 国产一区二区在线观看av| 亚洲av国产av综合av卡| 精品欧美一区二区三区在线| 久久热在线av| 午夜福利影视在线免费观看| 五月天丁香电影| 激情视频va一区二区三区| 午夜免费成人在线视频| 1024香蕉在线观看| 中文字幕人妻丝袜制服| 亚洲天堂av无毛| 一级片'在线观看视频| 成人精品一区二区免费| 国产无遮挡羞羞视频在线观看| 丝袜喷水一区| 性高湖久久久久久久久免费观看| 大型av网站在线播放| 99久久国产精品久久久| 夫妻午夜视频| 久久精品国产综合久久久| 日韩一区二区三区影片| 在线观看一区二区三区激情| 亚洲熟妇熟女久久| 欧美精品一区二区大全| 亚洲人成伊人成综合网2020| 国产av国产精品国产| 亚洲精品一卡2卡三卡4卡5卡| 他把我摸到了高潮在线观看 | 亚洲自偷自拍图片 自拍| 免费少妇av软件| 啪啪无遮挡十八禁网站| 精品国内亚洲2022精品成人 | 亚洲少妇的诱惑av| 成人免费观看视频高清| 99国产精品免费福利视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品一区二区在线观看99| 美女主播在线视频| 伦理电影免费视频| 在线 av 中文字幕| 女人高潮潮喷娇喘18禁视频| 久久久水蜜桃国产精品网| tocl精华| 看免费av毛片| 自线自在国产av| 亚洲精品粉嫩美女一区| 久久午夜亚洲精品久久| 高潮久久久久久久久久久不卡| 日本黄色视频三级网站网址 | 两个人看的免费小视频| 亚洲欧美日韩另类电影网站| 男女午夜视频在线观看| 精品国产一区二区三区久久久樱花| 视频区图区小说| 久久青草综合色| 国产高清视频在线播放一区| 12—13女人毛片做爰片一| 极品少妇高潮喷水抽搐| 飞空精品影院首页| 国产精品秋霞免费鲁丝片| 国产伦人伦偷精品视频| 久久久久国产一级毛片高清牌| 我要看黄色一级片免费的| 国产深夜福利视频在线观看| 免费在线观看完整版高清| 欧美精品一区二区免费开放| 人妻久久中文字幕网| 天天躁日日躁夜夜躁夜夜| 久久久久久久国产电影| 久久av网站| 夜夜爽天天搞| 欧美另类亚洲清纯唯美| 欧美人与性动交α欧美软件| 天天躁夜夜躁狠狠躁躁| 久久久国产精品麻豆| 欧美人与性动交α欧美精品济南到| 一级毛片电影观看| 水蜜桃什么品种好| 国产一区有黄有色的免费视频| 成人亚洲精品一区在线观看| 在线观看免费午夜福利视频| 搡老乐熟女国产| 老司机亚洲免费影院| 一个人免费看片子| 国产成人精品在线电影| 久久久国产成人免费| 色视频在线一区二区三区| 变态另类成人亚洲欧美熟女 | 欧美中文综合在线视频| 亚洲专区中文字幕在线| 亚洲精品自拍成人| 啦啦啦视频在线资源免费观看| 日韩成人在线观看一区二区三区| 一级片免费观看大全| 亚洲 欧美一区二区三区| 十八禁高潮呻吟视频| 男女边摸边吃奶| 亚洲欧美日韩高清在线视频 | 丰满迷人的少妇在线观看| 老熟女久久久| 男女下面插进去视频免费观看| 黄片播放在线免费| 啦啦啦免费观看视频1| 91国产中文字幕| av电影中文网址| 两个人免费观看高清视频| 叶爱在线成人免费视频播放| 男女高潮啪啪啪动态图| 亚洲国产毛片av蜜桃av| 最黄视频免费看| 超碰成人久久| 性高湖久久久久久久久免费观看| 丰满迷人的少妇在线观看| 亚洲伊人色综图| 麻豆av在线久日| 真人做人爱边吃奶动态| 欧美变态另类bdsm刘玥| 99久久精品国产亚洲精品| 变态另类成人亚洲欧美熟女 | 激情在线观看视频在线高清 | 女警被强在线播放| 成年人午夜在线观看视频| 国产欧美日韩精品亚洲av| 久久精品国产99精品国产亚洲性色 | 亚洲情色 制服丝袜| 极品少妇高潮喷水抽搐| 久久久久国内视频| 国产老妇伦熟女老妇高清| 一二三四社区在线视频社区8| 亚洲成国产人片在线观看| 午夜两性在线视频| 天天操日日干夜夜撸| 精品第一国产精品| 一个人免费看片子| 免费观看a级毛片全部| 久久久久网色| 一个人免费在线观看的高清视频| 久久人人爽av亚洲精品天堂| 午夜福利欧美成人| 国产区一区二久久| 12—13女人毛片做爰片一| 亚洲黑人精品在线| 在线观看免费高清a一片| 亚洲精品久久成人aⅴ小说| 久久香蕉激情| 又紧又爽又黄一区二区| 国产精品av久久久久免费| 欧美大码av| 波多野结衣av一区二区av| 黄片大片在线免费观看| 超碰成人久久| 精品久久久久久久毛片微露脸| 高清在线国产一区| 国产在线观看jvid| 亚洲国产av新网站| 国产无遮挡羞羞视频在线观看| 精品免费久久久久久久清纯 | 午夜精品国产一区二区电影| 交换朋友夫妻互换小说| 亚洲第一欧美日韩一区二区三区 | 久久久久国内视频| 欧美激情久久久久久爽电影 | 一边摸一边抽搐一进一出视频| 岛国毛片在线播放| 黄片播放在线免费| av欧美777| 日韩一卡2卡3卡4卡2021年| 亚洲熟女毛片儿| 免费观看av网站的网址| 亚洲av成人一区二区三| 国产在线视频一区二区| 欧美精品av麻豆av| 考比视频在线观看| 久久久国产一区二区| 久久人人97超碰香蕉20202| 母亲3免费完整高清在线观看| 亚洲国产欧美网| 91麻豆精品激情在线观看国产 | 一区二区av电影网| 激情视频va一区二区三区| 69av精品久久久久久 | 亚洲国产欧美一区二区综合| 亚洲av片天天在线观看| 啦啦啦 在线观看视频| 一级毛片电影观看| 久久精品aⅴ一区二区三区四区| 精品一区二区三区四区五区乱码| 午夜视频精品福利| 欧美成狂野欧美在线观看| 国产精品99久久99久久久不卡| 久久久国产精品麻豆| 人人澡人人妻人| 99re在线观看精品视频| av又黄又爽大尺度在线免费看| 成人av一区二区三区在线看| 丝袜人妻中文字幕| 亚洲美女黄片视频| 午夜老司机福利片| 男人舔女人的私密视频| 男人操女人黄网站| 亚洲成人免费电影在线观看| 99riav亚洲国产免费| www日本在线高清视频| videosex国产| 国精品久久久久久国模美| a级片在线免费高清观看视频| 亚洲精品乱久久久久久| 乱人伦中国视频| 久久九九热精品免费| 美女国产高潮福利片在线看| 午夜福利免费观看在线| 久久国产亚洲av麻豆专区| 蜜桃在线观看..| 色综合欧美亚洲国产小说| 精品亚洲成a人片在线观看| 麻豆国产av国片精品| 午夜激情av网站| 在线天堂中文资源库| 国产人伦9x9x在线观看| 国产主播在线观看一区二区| 美女扒开内裤让男人捅视频| 午夜视频精品福利| 免费在线观看黄色视频的| 一边摸一边做爽爽视频免费| 国产黄色免费在线视频| 欧美日韩精品网址| 老鸭窝网址在线观看| 一本综合久久免费| 叶爱在线成人免费视频播放| 欧美日韩国产mv在线观看视频| 亚洲精品美女久久av网站| 成人18禁在线播放| 国产欧美日韩一区二区三| 欧美日韩亚洲国产一区二区在线观看 | 丰满饥渴人妻一区二区三| 啪啪无遮挡十八禁网站| 亚洲欧美精品综合一区二区三区| 一级片免费观看大全| 岛国在线观看网站| 天堂中文最新版在线下载| 亚洲av成人一区二区三| 国产男女超爽视频在线观看| 欧美日韩国产mv在线观看视频| 国产精品.久久久| 日本av手机在线免费观看| 国产一卡二卡三卡精品| 9色porny在线观看| 久久久久精品人妻al黑| xxxhd国产人妻xxx| 国产一区二区激情短视频| 国产99久久九九免费精品| 一本一本久久a久久精品综合妖精| 午夜成年电影在线免费观看| 他把我摸到了高潮在线观看 | 日本黄色日本黄色录像| 日本av手机在线免费观看| 久久久精品免费免费高清| 欧美另类亚洲清纯唯美| 人妻 亚洲 视频| 欧美中文综合在线视频| 久久亚洲精品不卡| 在线播放国产精品三级| 欧美老熟妇乱子伦牲交| 久久影院123| 欧美变态另类bdsm刘玥| 久久国产精品人妻蜜桃| 欧美黑人欧美精品刺激| 精品国产乱子伦一区二区三区| 久热这里只有精品99| 国产成人av教育| 大型黄色视频在线免费观看| 日韩有码中文字幕| 国产深夜福利视频在线观看| 国产97色在线日韩免费| 国产一区有黄有色的免费视频| 国产精品亚洲av一区麻豆| 99国产极品粉嫩在线观看| 国产欧美日韩综合在线一区二区| 99国产极品粉嫩在线观看| 三级毛片av免费| 日韩欧美国产一区二区入口| 精品国产一区二区久久| 黑人巨大精品欧美一区二区蜜桃| 极品人妻少妇av视频| 日本五十路高清| 三级毛片av免费| av线在线观看网站| 夜夜夜夜夜久久久久| 在线观看人妻少妇| 自拍欧美九色日韩亚洲蝌蚪91| 涩涩av久久男人的天堂| 国产成人欧美在线观看 | 性色av乱码一区二区三区2| 丰满少妇做爰视频| 最近最新中文字幕大全免费视频| 五月天丁香电影| 欧美精品人与动牲交sv欧美| 欧美精品亚洲一区二区| 这个男人来自地球电影免费观看| 丰满饥渴人妻一区二区三| 一区二区三区精品91| 日韩中文字幕欧美一区二区| 久久性视频一级片| 悠悠久久av| 国产极品粉嫩免费观看在线| 欧美+亚洲+日韩+国产| 午夜老司机福利片| 在线观看人妻少妇| 国产精品久久电影中文字幕 | 亚洲精品粉嫩美女一区| 色播在线永久视频| 国产精品一区二区在线观看99| 丰满少妇做爰视频| 久久精品国产亚洲av高清一级| 国产精品国产高清国产av | 法律面前人人平等表现在哪些方面| 日韩精品免费视频一区二区三区| √禁漫天堂资源中文www| 久久精品国产亚洲av香蕉五月 | 成年人午夜在线观看视频| 大型av网站在线播放| 欧美日韩av久久| 日韩人妻精品一区2区三区| 黄色视频在线播放观看不卡| 国产精品成人在线| 国产欧美日韩一区二区三区在线| 亚洲免费av在线视频| 少妇精品久久久久久久| 黑人猛操日本美女一级片| 国产欧美日韩一区二区三| 天堂8中文在线网| 国产欧美日韩精品亚洲av| 在线观看免费高清a一片| 午夜福利免费观看在线| 欧美国产精品va在线观看不卡| av网站免费在线观看视频| 在线观看66精品国产| 成人手机av| 一边摸一边做爽爽视频免费| 亚洲人成电影免费在线| 新久久久久国产一级毛片| 欧美一级毛片孕妇| 少妇被粗大的猛进出69影院| 精品国内亚洲2022精品成人 | 亚洲黑人精品在线| 亚洲三区欧美一区| 黄色片一级片一级黄色片| 成人国产一区最新在线观看| 91麻豆av在线| 国产一卡二卡三卡精品| 真人做人爱边吃奶动态| 国产亚洲精品久久久久5区| 99精国产麻豆久久婷婷| 中文字幕色久视频| 欧美黑人精品巨大| 国产精品电影一区二区三区 | 免费不卡黄色视频| 日韩 欧美 亚洲 中文字幕| 色尼玛亚洲综合影院| 黄片小视频在线播放| 黄片播放在线免费| 国产精品自产拍在线观看55亚洲 | 久久亚洲精品不卡| 天堂8中文在线网| 成人亚洲精品一区在线观看| 热99国产精品久久久久久7| 中国美女看黄片| 成人特级黄色片久久久久久久 | 日韩一区二区三区影片| 精品国产乱码久久久久久小说| 亚洲少妇的诱惑av| 免费av中文字幕在线| 蜜桃在线观看..| 丰满迷人的少妇在线观看| 国产午夜精品久久久久久| 久久精品国产综合久久久| 免费女性裸体啪啪无遮挡网站| www日本在线高清视频| 1024视频免费在线观看| 999久久久国产精品视频| 欧美黄色淫秽网站| 考比视频在线观看| 99精品在免费线老司机午夜| 黄频高清免费视频| 一区二区三区激情视频| 免费在线观看黄色视频的| 日韩成人在线观看一区二区三区| 黄色丝袜av网址大全| 亚洲精品粉嫩美女一区| 一本大道久久a久久精品| 男女边摸边吃奶| 黄色 视频免费看| 日韩大片免费观看网站| 日韩欧美一区二区三区在线观看 | 午夜福利影视在线免费观看| 色精品久久人妻99蜜桃| 国产成人欧美| av又黄又爽大尺度在线免费看| 少妇被粗大的猛进出69影院| 黄色毛片三级朝国网站| 国产免费福利视频在线观看| 一边摸一边做爽爽视频免费| 99久久国产精品久久久| 午夜福利乱码中文字幕| xxxhd国产人妻xxx| 免费观看人在逋| 一级a爱视频在线免费观看| 亚洲色图 男人天堂 中文字幕| 新久久久久国产一级毛片| 欧美中文综合在线视频| 岛国毛片在线播放| 色婷婷久久久亚洲欧美| 日本av免费视频播放| 97在线人人人人妻| 黑人猛操日本美女一级片| 国产欧美日韩精品亚洲av| 高清欧美精品videossex| 这个男人来自地球电影免费观看| 久久人人97超碰香蕉20202| 高清av免费在线| 人人妻人人爽人人添夜夜欢视频| 亚洲欧美一区二区三区久久| 午夜福利在线观看吧| 美女视频免费永久观看网站| 一本—道久久a久久精品蜜桃钙片| 国产免费av片在线观看野外av| 精品国内亚洲2022精品成人 | 一级,二级,三级黄色视频| 亚洲精品在线美女| 国产区一区二久久| 亚洲av美国av| 国产在视频线精品| 亚洲av国产av综合av卡| 日韩制服丝袜自拍偷拍| 黄色丝袜av网址大全| 露出奶头的视频| 午夜激情av网站| 女警被强在线播放| 精品一品国产午夜福利视频| 少妇猛男粗大的猛烈进出视频| 亚洲国产精品一区二区三区在线| 免费在线观看影片大全网站| 777米奇影视久久| 亚洲男人天堂网一区| 9191精品国产免费久久| 久久精品成人免费网站| 新久久久久国产一级毛片| 国产片内射在线| videosex国产| 国产精品一区二区精品视频观看| 国产成人一区二区三区免费视频网站| 一区福利在线观看| 黄色毛片三级朝国网站| 色精品久久人妻99蜜桃| 国产深夜福利视频在线观看| 亚洲欧洲精品一区二区精品久久久| 亚洲精品国产精品久久久不卡| 下体分泌物呈黄色| 日本一区二区免费在线视频| av视频免费观看在线观看| 日日爽夜夜爽网站| 亚洲国产欧美一区二区综合| 亚洲伊人色综图| 伦理电影免费视频| 精品国产乱码久久久久久男人| 亚洲三区欧美一区| 老司机午夜福利在线观看视频 | 久久午夜综合久久蜜桃| 黑丝袜美女国产一区| 亚洲一区中文字幕在线| 亚洲人成电影免费在线| 久久亚洲精品不卡| 王馨瑶露胸无遮挡在线观看| 中亚洲国语对白在线视频| www.自偷自拍.com| 最新美女视频免费是黄的| 国产有黄有色有爽视频| 亚洲精品国产区一区二| 婷婷丁香在线五月| 成人免费观看视频高清| 在线观看免费视频日本深夜| 夜夜夜夜夜久久久久| 天堂动漫精品| 久久久久网色| 一区二区av电影网| 1024香蕉在线观看| 99国产精品一区二区三区| 亚洲天堂av无毛| 午夜福利乱码中文字幕| 午夜日韩欧美国产| a在线观看视频网站| 久久久国产欧美日韩av| 国产一区二区 视频在线| 亚洲国产精品一区二区三区在线| 久久久久久人人人人人| 亚洲国产中文字幕在线视频| 成人亚洲精品一区在线观看| 建设人人有责人人尽责人人享有的| 首页视频小说图片口味搜索| 国产欧美亚洲国产| 激情在线观看视频在线高清 | 久久精品亚洲精品国产色婷小说| 两性夫妻黄色片| av欧美777| 王馨瑶露胸无遮挡在线观看| 欧美成狂野欧美在线观看| 日韩 欧美 亚洲 中文字幕| 精品少妇内射三级| 丝袜人妻中文字幕| 啦啦啦免费观看视频1| 国产免费视频播放在线视频| 嫁个100分男人电影在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 日日摸夜夜添夜夜添小说| 久久久久久久久久久久大奶| 精品人妻1区二区| 青青草视频在线视频观看| 欧美精品人与动牲交sv欧美| 丁香六月天网| 久久青草综合色| 一边摸一边抽搐一进一小说 | 男女下面插进去视频免费观看| 好男人电影高清在线观看| 国产aⅴ精品一区二区三区波| 亚洲一码二码三码区别大吗| 日韩免费av在线播放| 精品第一国产精品|