• <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ì)解譯
    成人18禁在线播放| 精品久久久久久成人av| 日韩大尺度精品在线看网址 | 久久久久久大精品| 黄色片一级片一级黄色片| 久久久久精品国产欧美久久久| 少妇熟女aⅴ在线视频| 亚洲成人久久性| 在线观看免费视频网站a站| 欧美日韩亚洲综合一区二区三区_| 久99久视频精品免费| 制服人妻中文乱码| 亚洲人成网站在线播放欧美日韩| 1024香蕉在线观看| 18禁裸乳无遮挡免费网站照片 | 怎么达到女性高潮| 欧美一区二区精品小视频在线| 欧美午夜高清在线| 大香蕉久久成人网| 高潮久久久久久久久久久不卡| 欧美精品亚洲一区二区| 久久婷婷人人爽人人干人人爱 | 日本黄色视频三级网站网址| 欧美老熟妇乱子伦牲交| 男人的好看免费观看在线视频 | 在线av久久热| 久久久久国产精品人妻aⅴ院| 久久欧美精品欧美久久欧美| 香蕉丝袜av| 国内精品久久久久久久电影| 又黄又粗又硬又大视频| 亚洲国产看品久久| 一边摸一边抽搐一进一小说| 高清在线国产一区| 亚洲 欧美 日韩 在线 免费| 天堂√8在线中文| 香蕉国产在线看| 精品第一国产精品| 老汉色∧v一级毛片| 岛国视频午夜一区免费看| 久久久久九九精品影院| 欧美最黄视频在线播放免费| 大型黄色视频在线免费观看| 国内毛片毛片毛片毛片毛片| 91麻豆av在线| 天天一区二区日本电影三级 | 两个人看的免费小视频| 99久久久亚洲精品蜜臀av| 日本精品一区二区三区蜜桃| 亚洲色图 男人天堂 中文字幕| 精品第一国产精品| 亚洲 国产 在线| 国产一区二区激情短视频| 亚洲一卡2卡3卡4卡5卡精品中文| 在线国产一区二区在线| 人人妻人人澡欧美一区二区 | 一二三四在线观看免费中文在| 欧美在线黄色| 亚洲黑人精品在线| 久久影院123| 国产亚洲精品久久久久久毛片| 一进一出好大好爽视频| 操出白浆在线播放| netflix在线观看网站| 欧美+亚洲+日韩+国产| 亚洲无线在线观看| 99久久精品国产亚洲精品| 99riav亚洲国产免费| 精品欧美一区二区三区在线| 不卡av一区二区三区| 久久国产精品男人的天堂亚洲| 久久人妻熟女aⅴ| 狂野欧美激情性xxxx| 久久精品亚洲熟妇少妇任你| 精品一品国产午夜福利视频| 中亚洲国语对白在线视频| 日本三级黄在线观看| 少妇 在线观看| 免费av毛片视频| 丝袜人妻中文字幕| 99精品在免费线老司机午夜| 久久久国产欧美日韩av| а√天堂www在线а√下载| 色av中文字幕| 欧美亚洲日本最大视频资源| 叶爱在线成人免费视频播放| 亚洲成国产人片在线观看| 757午夜福利合集在线观看| 午夜福利欧美成人| 国产色视频综合| 国内精品久久久久精免费| ponron亚洲| 国产成人精品在线电影| 校园春色视频在线观看| 亚洲欧美日韩另类电影网站| 巨乳人妻的诱惑在线观看| 久久久精品欧美日韩精品| 亚洲欧美激情综合另类| 黄色成人免费大全| 午夜免费观看网址| 热99re8久久精品国产| 久久精品国产综合久久久| 国产亚洲精品第一综合不卡| 免费搜索国产男女视频| 免费在线观看日本一区| 亚洲欧美日韩无卡精品| 黄片大片在线免费观看| 精品卡一卡二卡四卡免费| 精品电影一区二区在线| 日日干狠狠操夜夜爽| 亚洲国产精品久久男人天堂| 国产精品综合久久久久久久免费 | 免费观看精品视频网站| 波多野结衣av一区二区av| 国产成人精品无人区| 法律面前人人平等表现在哪些方面| 亚洲欧美精品综合一区二区三区| 国产一区二区在线av高清观看| 非洲黑人性xxxx精品又粗又长| 一本综合久久免费| 18禁国产床啪视频网站| 很黄的视频免费| 欧美 亚洲 国产 日韩一| 日韩精品免费视频一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 天天一区二区日本电影三级 | 国产1区2区3区精品| av天堂在线播放| 亚洲精品中文字幕一二三四区| 亚洲国产精品久久男人天堂| 嫩草影视91久久| 欧美日韩亚洲综合一区二区三区_| 国产亚洲欧美精品永久| 每晚都被弄得嗷嗷叫到高潮| 国产精品自产拍在线观看55亚洲| 欧美一级毛片孕妇| 国产成人欧美| 久久精品国产亚洲av高清一级| 日本三级黄在线观看| 国产午夜福利久久久久久| 欧美日韩中文字幕国产精品一区二区三区 | 99国产精品一区二区三区| av在线播放免费不卡| 十八禁人妻一区二区| 啦啦啦韩国在线观看视频| 国产三级在线视频| 99riav亚洲国产免费| 成人欧美大片| 亚洲精品久久成人aⅴ小说| 一边摸一边抽搐一进一小说| 久久精品国产综合久久久| 久久狼人影院| 黄色片一级片一级黄色片| 亚洲国产精品合色在线| 天堂√8在线中文| 十分钟在线观看高清视频www| 757午夜福利合集在线观看| 亚洲第一青青草原| 夜夜夜夜夜久久久久| 手机成人av网站| 人人妻人人爽人人添夜夜欢视频| 亚洲成国产人片在线观看| 麻豆成人av在线观看| 亚洲av成人av| 夜夜夜夜夜久久久久| 国产精品1区2区在线观看.| 亚洲成av片中文字幕在线观看| 国产成人一区二区三区免费视频网站| 777久久人妻少妇嫩草av网站| 99国产精品免费福利视频| 非洲黑人性xxxx精品又粗又长| 两性午夜刺激爽爽歪歪视频在线观看 | 热re99久久国产66热| 波多野结衣av一区二区av| 又大又爽又粗| 香蕉久久夜色| 88av欧美| 欧美日韩亚洲国产一区二区在线观看| 亚洲专区中文字幕在线| 色精品久久人妻99蜜桃| 韩国av一区二区三区四区| 欧美黑人精品巨大| av福利片在线| 麻豆成人av在线观看| 久久精品人人爽人人爽视色| 最近最新中文字幕大全电影3 | 午夜a级毛片| 好男人在线观看高清免费视频 | 亚洲国产毛片av蜜桃av| 麻豆国产av国片精品| 亚洲一区二区三区不卡视频| 长腿黑丝高跟| svipshipincom国产片| 午夜视频精品福利| 免费人成视频x8x8入口观看| 精品久久久久久久人妻蜜臀av | 桃色一区二区三区在线观看| 亚洲 欧美 日韩 在线 免费| 亚洲aⅴ乱码一区二区在线播放 | 国产高清videossex| 亚洲专区中文字幕在线| 黄片大片在线免费观看| 国产精品国产高清国产av| 日本精品一区二区三区蜜桃| 法律面前人人平等表现在哪些方面| 欧美日本视频| 99国产精品免费福利视频| www国产在线视频色| 久久久久久免费高清国产稀缺| 久久人妻福利社区极品人妻图片| 中文字幕久久专区| 在线永久观看黄色视频| 亚洲成国产人片在线观看| 国产一区二区三区视频了| 中文字幕另类日韩欧美亚洲嫩草| 99精品在免费线老司机午夜| avwww免费| 可以在线观看的亚洲视频| 手机成人av网站| 国产精品日韩av在线免费观看 | 日韩三级视频一区二区三区| 制服人妻中文乱码| 巨乳人妻的诱惑在线观看| 99精品久久久久人妻精品| 9热在线视频观看99| 亚洲无线在线观看| 一级毛片女人18水好多| 精品国产美女av久久久久小说| 精品国产一区二区久久| 999久久久精品免费观看国产| 午夜亚洲福利在线播放| 久久久久久久久中文| 少妇被粗大的猛进出69影院| 成人国语在线视频| 成人国产一区最新在线观看| 久久天躁狠狠躁夜夜2o2o| 久久国产亚洲av麻豆专区| 一区二区日韩欧美中文字幕| 日本vs欧美在线观看视频| 亚洲电影在线观看av| 久久人妻福利社区极品人妻图片| 精品高清国产在线一区| 欧美激情 高清一区二区三区| 欧美中文日本在线观看视频| 97碰自拍视频| 午夜免费激情av| 欧美乱色亚洲激情| 久久久精品欧美日韩精品| 91老司机精品| 丰满的人妻完整版| 成熟少妇高潮喷水视频| 久久国产乱子伦精品免费另类| tocl精华| 亚洲av五月六月丁香网| 国产精品av久久久久免费| 人妻丰满熟妇av一区二区三区| 在线观看一区二区三区| 精品久久蜜臀av无| av中文乱码字幕在线| 香蕉国产在线看| 老司机在亚洲福利影院| 最近最新中文字幕大全免费视频| 国产伦一二天堂av在线观看| 久久久精品欧美日韩精品| 亚洲在线自拍视频| 精品国产美女av久久久久小说| 久久久国产成人精品二区| 黄色 视频免费看| 一个人观看的视频www高清免费观看 | 黄网站色视频无遮挡免费观看| 极品教师在线免费播放| 国产激情欧美一区二区| 岛国视频午夜一区免费看| 免费久久久久久久精品成人欧美视频| 韩国精品一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲中文字幕日韩| 日韩av在线大香蕉| 久久精品亚洲熟妇少妇任你| 午夜亚洲福利在线播放| 午夜福利影视在线免费观看| 亚洲性夜色夜夜综合| 精品一区二区三区av网在线观看| 国产精品,欧美在线| 色av中文字幕| 99re在线观看精品视频| 自拍欧美九色日韩亚洲蝌蚪91| 欧美+亚洲+日韩+国产| 制服诱惑二区| 亚洲电影在线观看av| 老鸭窝网址在线观看| av电影中文网址| 在线天堂中文资源库| 黄色丝袜av网址大全| 国产一区在线观看成人免费| 亚洲欧美精品综合久久99| 亚洲精品在线观看二区| 欧美成狂野欧美在线观看| 在线国产一区二区在线| www.999成人在线观看| 精品熟女少妇八av免费久了| 国产熟女xx| 亚洲国产看品久久| 大香蕉久久成人网| 国产成+人综合+亚洲专区| 亚洲激情在线av| 日本一区二区免费在线视频| 黄色成人免费大全| 亚洲久久久国产精品| 国产单亲对白刺激| 一二三四在线观看免费中文在| 美女高潮到喷水免费观看| 亚洲欧美精品综合久久99| 天天躁夜夜躁狠狠躁躁| 中亚洲国语对白在线视频| 中文字幕久久专区| 悠悠久久av| 淫秽高清视频在线观看| 狂野欧美激情性xxxx| 精品福利观看| 男女床上黄色一级片免费看| 久久这里只有精品19| 欧美精品亚洲一区二区| 亚洲熟妇熟女久久| 精品欧美一区二区三区在线| 国产99白浆流出| 极品教师在线免费播放| 成人18禁在线播放| 久久人人97超碰香蕉20202| 亚洲七黄色美女视频| 黄片播放在线免费| 女人高潮潮喷娇喘18禁视频| 一级片免费观看大全| 一级毛片精品| 国产xxxxx性猛交| 给我免费播放毛片高清在线观看| 亚洲一区二区三区不卡视频| 在线免费观看的www视频| 高清在线国产一区| 国产一卡二卡三卡精品| 天堂影院成人在线观看| 日韩欧美三级三区| 淫妇啪啪啪对白视频| 午夜免费观看网址| 久久久久久久久免费视频了| 免费高清视频大片| 国产精品电影一区二区三区| 久久久国产精品麻豆| 91精品国产国语对白视频| 亚洲精品粉嫩美女一区| 热re99久久国产66热| 色综合站精品国产| 黑人操中国人逼视频| 久久久精品欧美日韩精品| 午夜日韩欧美国产| 亚洲无线在线观看| 国产精品电影一区二区三区| videosex国产| av电影中文网址| 国内毛片毛片毛片毛片毛片| 久久精品亚洲熟妇少妇任你| 免费搜索国产男女视频| 国产亚洲av高清不卡| 在线观看午夜福利视频| 日本黄色视频三级网站网址| 欧美色视频一区免费| 一二三四在线观看免费中文在| 日本 欧美在线| 9热在线视频观看99| 啪啪无遮挡十八禁网站| 侵犯人妻中文字幕一二三四区| 久久精品国产亚洲av高清一级| 久久精品亚洲熟妇少妇任你| 老鸭窝网址在线观看| 亚洲欧美精品综合一区二区三区| 十八禁人妻一区二区| 夜夜看夜夜爽夜夜摸| 亚洲人成电影观看| 一区二区三区激情视频| 欧美成狂野欧美在线观看| 亚洲精品中文字幕在线视频| 人人妻人人澡人人看| 黑人巨大精品欧美一区二区蜜桃| av片东京热男人的天堂| 亚洲欧美日韩无卡精品| 国产激情久久老熟女| 男女下面进入的视频免费午夜 | 性欧美人与动物交配| 国产成人系列免费观看| √禁漫天堂资源中文www| 国产精品久久久人人做人人爽| 黄色丝袜av网址大全| 1024香蕉在线观看| 脱女人内裤的视频| 亚洲性夜色夜夜综合| 欧美乱色亚洲激情| av网站免费在线观看视频| 国产精品一区二区精品视频观看| 亚洲五月婷婷丁香| 日本黄色视频三级网站网址| 亚洲人成电影免费在线| 国产野战对白在线观看| 精品日产1卡2卡| 日本 av在线| svipshipincom国产片| 1024香蕉在线观看| 日日夜夜操网爽| 91成年电影在线观看| 欧美乱码精品一区二区三区| 大型av网站在线播放| 一二三四在线观看免费中文在| a在线观看视频网站| 国产一区二区三区视频了| 在线观看免费午夜福利视频| 午夜福利18| 精品午夜福利视频在线观看一区| 久久亚洲精品不卡| 91大片在线观看| 深夜精品福利| 午夜日韩欧美国产| 亚洲av熟女| 变态另类成人亚洲欧美熟女 | 国产男靠女视频免费网站| 国产区一区二久久| 亚洲国产看品久久| 久久人妻av系列| 老汉色av国产亚洲站长工具| 亚洲午夜理论影院| 亚洲aⅴ乱码一区二区在线播放 | 黄色成人免费大全| 99在线视频只有这里精品首页| 欧美一区二区精品小视频在线| 极品人妻少妇av视频| 国产精品久久久av美女十八| svipshipincom国产片| 88av欧美| 国产亚洲精品综合一区在线观看 | 久久久国产精品麻豆| av电影中文网址| 大型av网站在线播放| 国产精品亚洲一级av第二区| 99香蕉大伊视频| 黑人巨大精品欧美一区二区蜜桃| 精品日产1卡2卡| 国产欧美日韩一区二区三| 精品国产乱子伦一区二区三区| 90打野战视频偷拍视频| 亚洲国产精品合色在线| 国产1区2区3区精品| 欧美激情极品国产一区二区三区| 中文字幕精品免费在线观看视频| 国产精品久久视频播放| 少妇粗大呻吟视频| 久久久国产成人精品二区| 久久精品aⅴ一区二区三区四区| 国产在线精品亚洲第一网站| 成人免费观看视频高清| 久久精品亚洲熟妇少妇任你| 男人的好看免费观看在线视频 | 身体一侧抽搐| 99久久精品国产亚洲精品| 国产又色又爽无遮挡免费看| 国产成人av激情在线播放| 亚洲精品国产精品久久久不卡| 欧美日本视频| 久久天躁狠狠躁夜夜2o2o| 国产成人欧美| 久久精品国产99精品国产亚洲性色 | 亚洲,欧美精品.| 欧美一级a爱片免费观看看 | 老汉色av国产亚洲站长工具| 国产精品亚洲美女久久久| 电影成人av| 精品国产超薄肉色丝袜足j| 午夜老司机福利片| 亚洲精品久久成人aⅴ小说| 亚洲免费av在线视频| 国产真人三级小视频在线观看| 久久人妻av系列| 90打野战视频偷拍视频| 午夜激情av网站| 麻豆国产av国片精品| 午夜影院日韩av| 国产精品久久电影中文字幕| 在线观看一区二区三区| 亚洲av电影不卡..在线观看| 精品卡一卡二卡四卡免费| 亚洲全国av大片| 高清毛片免费观看视频网站| 天天一区二区日本电影三级 | 精品久久久久久成人av| 午夜免费成人在线视频| 精品熟女少妇八av免费久了| 亚洲中文字幕一区二区三区有码在线看 | 亚洲精品美女久久av网站| 黄色丝袜av网址大全| 国产xxxxx性猛交| www国产在线视频色| 午夜影院日韩av| 国产区一区二久久| 国产男靠女视频免费网站| 三级毛片av免费| 97碰自拍视频| 精品无人区乱码1区二区| 神马国产精品三级电影在线观看 | 亚洲第一av免费看| 久久草成人影院| 嫩草影视91久久| 欧美成人性av电影在线观看| 精品久久久精品久久久| 国产伦人伦偷精品视频| 亚洲成人免费电影在线观看| 国产亚洲欧美98| 黄色女人牲交| 麻豆一二三区av精品| 国产精品日韩av在线免费观看 | 国产精华一区二区三区| 亚洲少妇的诱惑av| 制服人妻中文乱码| 精品国产乱码久久久久久男人| 中文亚洲av片在线观看爽| 国产成人av教育| 一a级毛片在线观看| 啦啦啦韩国在线观看视频| 国产高清有码在线观看视频 | 中文字幕av电影在线播放| 高清黄色对白视频在线免费看| 精品久久久久久成人av| 亚洲熟妇中文字幕五十中出| 无限看片的www在线观看| 久久性视频一级片| 91成人精品电影| 高清在线国产一区| 欧美激情极品国产一区二区三区| 一边摸一边抽搐一进一出视频| 别揉我奶头~嗯~啊~动态视频| 美女大奶头视频| 夜夜躁狠狠躁天天躁| 亚洲欧美一区二区三区黑人| 少妇粗大呻吟视频| 法律面前人人平等表现在哪些方面| 国产一区二区激情短视频| 亚洲色图av天堂| 人人妻人人澡人人看| 国产蜜桃级精品一区二区三区| 热re99久久国产66热| 国产成人精品久久二区二区91| 999精品在线视频| 精品日产1卡2卡| 国产午夜精品久久久久久| 欧美精品啪啪一区二区三区| 一级作爱视频免费观看| 国产主播在线观看一区二区| 久久国产亚洲av麻豆专区| 亚洲精华国产精华精| 亚洲伊人色综图| 欧美大码av| 午夜亚洲福利在线播放| 亚洲成人久久性| 国产精品久久电影中文字幕| 露出奶头的视频| 女同久久另类99精品国产91| 在线av久久热| 亚洲精品美女久久av网站| 手机成人av网站| 美女 人体艺术 gogo| 国产一卡二卡三卡精品| 无限看片的www在线观看| 真人一进一出gif抽搐免费| 亚洲第一av免费看| 国产1区2区3区精品| 视频在线观看一区二区三区| 久久久国产欧美日韩av| 脱女人内裤的视频| 亚洲五月天丁香| 神马国产精品三级电影在线观看 | 多毛熟女@视频| 天天添夜夜摸| 亚洲精品国产区一区二| 欧美中文综合在线视频| 一个人免费在线观看的高清视频| 亚洲国产高清在线一区二区三 | 日韩大尺度精品在线看网址 | 亚洲专区中文字幕在线| 级片在线观看| 亚洲av第一区精品v没综合| 日本免费一区二区三区高清不卡 | or卡值多少钱| 国产午夜福利久久久久久| 一卡2卡三卡四卡精品乱码亚洲| 美女高潮到喷水免费观看| 国产亚洲av高清不卡| 可以免费在线观看a视频的电影网站| 日韩三级视频一区二区三区| 侵犯人妻中文字幕一二三四区| 欧美中文日本在线观看视频| 欧美激情极品国产一区二区三区| 女同久久另类99精品国产91| 中文字幕色久视频| 色综合站精品国产| 一区二区三区激情视频| 国产亚洲欧美精品永久| 黑丝袜美女国产一区| 黄色 视频免费看| 黄色丝袜av网址大全| 18禁美女被吸乳视频| 自线自在国产av| 久久久久久久久免费视频了| av网站免费在线观看视频| 国产野战对白在线观看| 窝窝影院91人妻| 嫁个100分男人电影在线观看| 国产在线观看jvid| 在线国产一区二区在线| 精品欧美国产一区二区三|