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

    金屬離子污染砂土復(fù)電阻率的時變特征及形成機制

    2019-03-29 07:40:50王澤亞能昌信潘永泰劉玉強
    中國環(huán)境科學(xué) 2019年3期
    關(guān)鍵詞:土柱極化電阻率

    王澤亞,徐 亞,董 路,能昌信,潘永泰,劉玉強

    ?

    金屬離子污染砂土復(fù)電阻率的時變特征及形成機制

    王澤亞1,徐 亞2*,董 路2,能昌信2,潘永泰1,劉玉強2

    (1.中國礦業(yè)大學(xué)(北京)化學(xué)與環(huán)境工程學(xué)院,北京 100083;2.中國環(huán)境科學(xué)研究院土壤與固體廢物研究所,北京 100012)

    為研究不同污染情況下土壤復(fù)電阻率的時間變異特征及控制機制,設(shè)計了人工砂土和水2種不同的介質(zhì)作對照,以NaCl、MnCl2和模擬海水作為污染物分別模擬海水入侵和金屬離子污染場景,測量100~103Hz范圍內(nèi)的相位變化并觀測其時變特性.結(jié)果顯示,無論在純液相介質(zhì)中還是在土壤介質(zhì)中,相位值都隨著污染物濃度的增大而減小,隨頻率的增大而增大.隨著污染物濃度的變化,土壤介質(zhì)中的相位比水相介質(zhì)中更加敏感.在土壤和水的背景環(huán)境下,通入1g/L的NaCl溶液,在1Hz條件下測量,土壤相位下降了約42%,而水相介質(zhì)幾乎看不出相位變化.針對固定頻率進行多次測量的結(jié)果表明,水相介質(zhì)中的相位值穩(wěn)定而土壤介質(zhì)中的相位值在中高頻段(>16Hz)表現(xiàn)出了明顯的時變特征,這種時變性也隨著污染物濃度的增大而減小.在102Hz時,測得沙柱背景值的相位浮動范圍在-4.05~-11.67mrad,方差達到4.8,而在水相溶液中方差僅為0.09.復(fù)電阻率法對金屬離子污染土壤的反應(yīng)很靈敏,有較好的指示意義,其參數(shù)變化主要與孔隙水中金屬離子含量有關(guān).但由于多相介質(zhì)的相界面特征以及孔隙結(jié)構(gòu)的復(fù)雜性,使得交變電流下電荷載流子在介質(zhì)中的遷移極化規(guī)律混亂,造成相位波動.

    相位;介電特性;物探;聯(lián)合反演

    環(huán)境地球物理探測具有無損、高效、低成本和高分辨率等特點,目前已被廣泛用于污染土壤和地下水中的重金屬和POPs類污染物的檢測和動態(tài)監(jiān)測[1-2].電法勘探作為地球物理探測的一個重要分支,根據(jù)受污染場地與周邊未污染的環(huán)境介質(zhì)存在導(dǎo)電性、電化學(xué)活性等特性的差異,在地質(zhì)勘探、污染探測等領(lǐng)域中發(fā)揮著重要作用[3-4].

    高密度電阻率法(ERT)能有效識別并劃定異常低阻區(qū)域,無論是室內(nèi)模擬還是現(xiàn)場調(diào)查,都已經(jīng)有大量研究和應(yīng)用.早期的研究關(guān)注土壤類型、含水率、溫度等非污染因素對電阻率特性的影響[5-7];近些年的研究主要討論土壤中的無機鹽[8-10]和有機物[11-12]組分與電阻率之間的關(guān)系問題,證實ERT能夠圈定污染范圍、追蹤污染羽流.然而電阻率法的弊端也在研究中逐漸凸顯,由于獲得的信息較為單一,電阻率與污染情況(污染物濃度和污染物類型)的對應(yīng)關(guān)系很難明確,因而在有些情況下無法判斷是否受到污染.究其原因是因為低阻體的異參同效效應(yīng),即不同的場地參數(shù)或條件,如含水率、污染物類型和濃度,可能導(dǎo)致相同的低阻特性.

    在傳統(tǒng)電阻率法基礎(chǔ)上發(fā)展起來的復(fù)電阻率法可以提供幅值和相位兩個參數(shù),因而很好地彌補了電阻率法單一參數(shù)的缺點[13].復(fù)電阻率法最初用于礦產(chǎn)資源的探測,但近年來也逐漸開始用于水文地質(zhì)和環(huán)境污染探測.研究表明,當(dāng)污染物入侵到土壤中時,會對土壤環(huán)境的理化特性造成影響,包括土壤孔隙流體特性和相界面上的顆粒表面特性[14-15],進而改變土壤的復(fù)電阻率特性.基于這種特性,國內(nèi)外學(xué)者開展了一系列研究.基礎(chǔ)性的研究包括重金屬污染土壤的復(fù)電阻率特性研究[16-17]以及有機污染物對土壤復(fù)電阻率參數(shù)的影響[18-20]等;國外學(xué)者在此基礎(chǔ)上進一步研究了復(fù)電阻率模型和反演算法[21-23],完善復(fù)電阻率法理論基礎(chǔ).

    現(xiàn)有研究主要集中于復(fù)電阻率參數(shù)與污染特征之間的關(guān)系問題,而很少有關(guān)于土壤復(fù)電阻率參數(shù)時變特性的研究.就導(dǎo)電機理而言,土壤和巖石的電學(xué)特性與金屬等有著較大差異[24].土壤主要依靠孔隙溶液的自由離子導(dǎo)電,導(dǎo)電離子會在持續(xù)電流的作用下發(fā)生動態(tài)流動,在相界面處發(fā)生累積,導(dǎo)致復(fù)電阻率指標(biāo)發(fā)生變化.而測量過程中的頻率變化可能進一步導(dǎo)致帶電離子的重分布,使得其時間變異進一步復(fù)雜化.土壤的這種特性嚴重影響了電法應(yīng)用于環(huán)境污染探測的重復(fù)性、精確性和動態(tài)可比性,因此對其時間變異性研究顯得尤其重要.本研究的主要目的就是討論環(huán)境污染探測過程中土壤復(fù)電阻率的時變特性,探索其變化規(guī)律和范圍,嘗試對其進行合理性的解釋.

    1 材料與方法

    1.1 實驗裝置

    實驗裝置如圖1所示,最左側(cè)為高200mm、內(nèi)徑38mm的有機玻璃柱.柱體最底部為入水口,頂部為出水口,水流從下往上流動,以保證均勻、穩(wěn)定的水飽和條件.土柱底部和頂部放置螺旋銅線圈作為供電電極,柱體右側(cè)布設(shè)86mm間距的采樣電極.采樣電極采用氯化銀不極化電極(直徑3.8mm)以減少極化效應(yīng);柱體右側(cè)為復(fù)電阻率采集裝置,與供電電極和采樣電極相連.裝置最右端為信號輸入輸出設(shè)備,在信號輸出端與玻璃柱間串聯(lián)精密電阻ref,采集電阻兩端的信號,作為參考電流信號.

    圖1 實驗裝置 Fig.1 Schematic of experimental setup

    1.2 測量原理

    1.3 模擬電路實驗

    圖2 模擬電路設(shè)計 Fig.2 Analogous circuit schematic

    為了驗證測量裝置的有效性,進而排除裝置本身的時變效應(yīng)干擾,設(shè)計了圖2所示的模擬電路代替圖1中最左端玻璃柱.因為模擬電路主要用電子元器件,其時間變異可以忽略.因此如果針對模擬電路的測試結(jié)果顯示出時變特性,那就可以判斷是裝置本身的問題;反之則說明測量裝置沒有問題.圖中,AC為供電信號,1為電流參考電阻(即ref),2和1兩端作為復(fù)電阻率采集端.

    設(shè)2和1兩端阻抗為2,則有

    復(fù)電阻率的實部()和虛部()分別為

    由式(3)和(4)可得復(fù)電阻率的幅值|Z|和相位角理論值分別為:

    實驗所選頻率為0~1000Hz,測量過程中保持除頻率以外其他條件不變,在一個頻率下連續(xù)采集50次數(shù)據(jù)后,換到下一個選定頻率繼續(xù)重復(fù)操作.下面的土柱實驗和水柱實驗均遵循這種采集方法.

    1.4 土柱實驗

    表1 模擬海水組成成分[27]

    由于自由離子在土壤中復(fù)雜的遷移規(guī)律,使得復(fù)電阻率的時變特性復(fù)雜化,這給數(shù)據(jù)分析帶來了很大的困難.針對這一現(xiàn)象,本研究設(shè)計向土柱中通入不同種類及濃度的污染物(氯化鈉、氯化錳及表1的模擬海水),分別模擬海水入侵以及金屬陽離子污染場地的土壤環(huán)境.人工土壤使用砂和黏土(由95%細砂顆粒和5%膨潤土組成)以模擬常見砂土[25-26],細砂顆粒經(jīng)過100目紗網(wǎng)濾出,加入質(zhì)量分數(shù)為5%的膨潤土采用機械攪拌使之混合均勻,填充至裝置中.采用統(tǒng)一灌裝的方法將土樣填充至3個裝置中,按表2中的步驟依次泵入溶液至飽和狀態(tài)并測量.

    表2 土柱實驗各階段污染物 Table 2 Pollutants in different stages of soil column experiment

    1.5 水柱實驗

    為了排除供電電極、測量電極以及裝置的影響,同時考慮到土壤與水之間形成的多相介質(zhì)的復(fù)雜性,作為對照,將上述土柱中的土壤掏空,用純?nèi)芤禾娲?參照土柱實驗,水柱實驗選擇NaCl和MnCl2作為污染物進行實驗.

    2 結(jié)果與討論

    2.1 模擬電路相位測量結(jié)果

    圖3 模擬電路相位測量值的絕對誤差(Ea)與相對誤差(Er) Fig.3 Absolute error (Ea) and relative error (Er) of phase of analogous circuit

    根據(jù)式(2)~(6)可計算出模擬電路每個頻率的相位理論值,與實驗所得相位期望值相比較可計算出實驗誤差.實驗相位期望值取50次的平均值,因此在一定程度上消除了隨機誤差.但是由于實驗室所用電阻和電容本身等存在一定誤差,在計算的時候會有系統(tǒng)誤差.由圖3可見測量值的相對誤差基本都小于1%,在可接受范圍之內(nèi).

    圖4 模擬電路相位期望及方差 Fig.4 The expectations and variances of phase of analogous circuit

    如圖4所示,相位期望在所選的頻率范圍內(nèi)都是負值,與頻率近似成一次函數(shù)關(guān)系,即采集的電壓相位要滯后于電流相位,且隨著頻率的增大,滯后的現(xiàn)象越明顯;方差都小于0.3,相比于期望值微乎其微.因而認為模擬電路的相位測量值很精確,相位穩(wěn)定,信號輸入輸出設(shè)備、實驗電路、參考電阻ref等因素對實驗的精確性沒有影響.

    2.2 土柱相位測量結(jié)果

    如圖5所示,在低頻段(<10Hz),數(shù)據(jù)的波動性很小,而在102Hz和103Hz下的相位值卻變得很不穩(wěn)定.圖中每條折線的變化趨勢不具有規(guī)律性,即每次的測量值會隨機落在一定范圍的區(qū)間內(nèi),因此可以認為每次測量的相位值之間相互獨立.為了更好地對結(jié)果進行分析,選擇了圖4的方法分別計算數(shù)學(xué)期望和方差.

    如圖6所示,各個曲線的變化不再像模擬電路一樣與直線的擬合度那么高,且變化的幅度小了很多.以背景值為例,從-6mard(16Hz)到-10.7mrad (1000Hz),而模擬電路則是從-5mrad(16Hz)到-307mrad (1000Hz).在土柱中只通水的情況下,其相位最大(這里說的大指的是絕對值,該值越大表明電壓相位落后電流相位越明顯,下同).在通入污染物之后,相位開始變小,且濃度越高值越小.對比不同污染物的相位曲線,在相同的質(zhì)量濃度下,由于海水以NaCl為主要組成部分,其表現(xiàn)出的電學(xué)特性與NaCl幾乎相同,而MnCl2的相位相對較小.由于電場的作用,材料中的各種電荷載流子能夠移動和極化,極化作用越強,測得的相位差越大.本文主要分析水分子的取向極化和介質(zhì)中離子在雙電層上的極化兩個方面.

    圖5 水飽和土柱中單一頻率相位重復(fù)測量值 Fig.5 Phase repeated measurements in water saturated sand column at single frequency

    在含水的介質(zhì)中,水分子的取向極化是普遍存在的.由于水分子是極性分子,在外加電場的作用下會旋轉(zhuǎn)并朝著一個方向排列,因而產(chǎn)生取向極化[28].在只有水通入柱中的情況下,介質(zhì)中的自由離子含量少,水分子在電場方向的取向幾乎不受限制,取向極化作用較強,因而有較大的相位差.而當(dāng)溶液中含有NaCl或MnCl2時,陰陽離子與水分子相互作用,陽離子(Na+、Mn2+)被水分子的氧吸引,Cl-被水分子的氫吸引,使得水偶極子的旋轉(zhuǎn)受到抑制,導(dǎo)致較低的取向極化[29].且溶液中電解質(zhì)濃度越高,對水偶極子旋轉(zhuǎn)的抑制作用越明顯,取向極化效果越低.在相同的質(zhì)量濃度下,MnCl2的離子數(shù)比NaCl多,因而對水分子取向極化的抑制作用也強于NaCl,所以相位差小于NaCl溶液.

    圖6 土柱中通入不同溶液的相位期望 Fig.6 The expectations of phase with different solutions injected into sand column

    對于多相介質(zhì)來說,相界面極化在整個極化作用中占據(jù)主導(dǎo)地位.在固液相界面上,當(dāng)污染物濃度變化時,自由離子含量也隨之變化,導(dǎo)致雙電層結(jié)構(gòu)改變[14, 30].隨著污染物濃度增大,溶液中與定位離子同電性的離子濃度雖大,但由于靜電斥力都被排斥到離相界面較遠處;而由于靜電引力的作用,反號離子濃度的增大則會大量進入滑動面內(nèi)(斯特恩層),以平衡土壤顆粒本身所帶的負電荷,使得滑動面以內(nèi)的凈電荷減小.隨著滑動面內(nèi)靜電荷的減小,擴散層中離子濃度降低,同時也使得擴散層的厚度減小,電動電位(z電位)減小,導(dǎo)致雙電層上的極化作用被抑制.

    關(guān)于Mn2+和Na+所造成的相位差異, Kaya等[29]研究表明相同濃度的CaCl2比NaCl對擴散層厚度的削弱作用更強,然而并沒有解釋原因;而Vaudelet等[31]則對CuSO4和NaCl溶液對土壤相位的影響給出了詳細的解釋,由于銅是2價,其比1價的鈉有更強的吸附能力,因而在雙電層上的流動性變差,極化效應(yīng)弱.對于本文中的錳離子,其價態(tài)總是高于鈉離子,因而MnCl2的極化效應(yīng)要弱于NaCl.

    此外還注意到,在100Hz附近,各曲線的走勢都有了一些變化,且污染物濃度越低,這種現(xiàn)象越明顯.小于100Hz的時候,相位的上升幅度(斜率)大;而在100Hz之后,上升幅度明顯減緩.這與此前的一些研究結(jié)果一致[26,30],在1000Hz以內(nèi)多孔介質(zhì)的極化作用主要由2種機制主導(dǎo).較低頻率(<100Hz)的極化作用可以用雙電層解釋,土壤顆粒的極化主要歸因于雙電層極化,其來源于顆粒表面吸附的配衡離子和電荷沿顆粒表面的切向運動;而在較高頻率(3100Hz)下,雙電層的極化作用已不隨頻率的變化而改變,此時則以麥克斯韋-瓦格納(MW)極化為主,其與多相介質(zhì)之間的孔隙電荷的積累有關(guān).

    由圖7可知,土柱的相位測量值出現(xiàn)了明顯的波動性,相比于模擬電路,方差要大很多.與期望類似,在土柱中只通水的情況下,方差達到最大,在通入鹽溶液之后,方差開始變小,且濃度越高該值越小.然而方差隨頻率的關(guān)系也比較特殊,是呈先增大后減小的趨勢.以背景值為例,在100Hz的時候方差達到最大值4.8,這對相位的測量精度造成了很大影響.

    有研究曾將土壤的介電特性與楊氏模量相對照,認為介電常數(shù)對外加電場的響應(yīng)就像楊氏模量對施加外力的響應(yīng)[29].隨著污染物濃度的提高,介質(zhì)中離子濃度增加,對環(huán)境因素變化的靈敏度降低,系統(tǒng)的穩(wěn)定性提高.但由于土壤的極化機理較為復(fù)雜,簡單的將其與楊氏模量相對應(yīng)很難解釋土壤復(fù)電阻率的變化規(guī)律,進而猜測這種不穩(wěn)定性與土壤的孔隙分布不規(guī)律有關(guān).由于土壤結(jié)構(gòu)復(fù)雜,其間孔隙的大小結(jié)構(gòu)都各不相同,因而水分充滿孔隙間的時候,各種電荷載流子在不同區(qū)域的水中的運動規(guī)律也不相同.再將其與雙電層聯(lián)合考慮,金屬陽離子在雙電層上的吸附與解吸在不同頻率下的表現(xiàn)可能有所差異.低頻段,離子有足夠的時間移動,吸附與解吸的速度保持平衡,因而相位穩(wěn)定;中頻段,由于孔隙分布不均勻,且電場方向變化較快,很可能出現(xiàn)陽離子吸附與解吸的速度不平衡的情況,因而相位浮動較大;高頻段,二者的速度都跟不上電場變化頻率的時候,產(chǎn)生了相對平衡.

    2.3 水柱相位測量結(jié)果

    由圖8可知,期望曲線相比于土柱而言,在較低頻率(<100Hz)下的相位表現(xiàn)有著比較明顯的區(qū)別,土壤中隨著濃度的變化各點位的區(qū)別很明顯,而在水中各濃度的相位基本都為0,即沒有極化現(xiàn)象.說明在低頻段,土壤中的極化主要由雙電層極化主導(dǎo).對比水相溶液相同濃度的NaCl和MnCl2曲線,沒有特別明顯的差異,甚至濃度相差10倍的差異也很小,不足以分辨.除此之外,水柱整體的相位隨頻率變化規(guī)律與在土壤中的類似.

    方差較土壤小了很多,基本與模擬電路保持在一個水平,表明水中相位測量值的穩(wěn)定性好,這種相位波動性是土壤顆粒與水組成的多相介質(zhì)的特性.因為在液相中,各種電荷載流子的運動不會受到土壤顆粒的阻礙,同時也沒有界面間雙電層的作用,只在電場力的作用下規(guī)律排列.

    由于水柱中相位較為穩(wěn)定,因而猜想造成土柱相位不穩(wěn)定的原因也可能是帶電粒子的電泳現(xiàn)象.在只通自來水的情況下,介質(zhì)中自由離子濃度低,z電位較高,電動現(xiàn)象明顯,因而相位出現(xiàn)較強的時變性;在通入鹽溶液后,隨著離子濃度的增大,斯特恩層中反號離子的濃度逐漸增大,導(dǎo)致z電位降低,電動現(xiàn)象變?nèi)?因而相位也更穩(wěn)定.

    3 結(jié)論

    3.1 復(fù)電阻率相位的滯后對金屬陽離子的存在非常敏感,復(fù)電阻率法在判別金屬離子污染程度和不同的污染物有很好的指示作用.

    3.2 污染物濃度對相位值有較大影響,隨著鹽含量的增加,土壤和水相溶液中測得的相位差都呈下降趨勢,且在土壤中更明顯.水相介質(zhì)的相位主要受水分子的取向極化作用影響,而土壤的相位主要受水與土壤顆粒之間的相界面極化影響.

    3.3 相位測量值的時變特性是土壤與水形成的混合介質(zhì)的特性,鹽濃度的增加會抑制這種性質(zhì),其可能與孔隙分布不規(guī)則導(dǎo)致的交變電流下孔隙離子遷移混亂有關(guān),也可能與土壤中膠體顆粒的電泳現(xiàn)象有關(guān).相位的波動性還有待進一步研究,其亦可能與土壤顆粒在Ag-AgCl電極上發(fā)生的物理、化學(xué)變化有關(guān).

    [1] Knight R, Pyrak-Nolte L J, Slater L, et al. Geophysics at the interface: Response of geophysical properties to solid-fluid, fluid-fluid, and solid-solid interfaces [J]. Reviews of Geophysics, 2010,48(4):324-336.

    [2] 徐 輝,詹良通,穆青翼,等.高有機質(zhì)含量垃圾的含水量監(jiān)測試驗研究——利用表面處理的TDR探頭[J]. 中國環(huán)境科學(xué), 2014,34(8): 2030-2039. Xu H, Zhan L T, Mu Q Y, et al. Experimental study on monitoring moisture content in municipal solid wastes with high organic content: using coated TDR probe [J].China Environmental Science, 2014,34(8): 2030-2039.

    [3] 李 忠,能昌信,寧書年,等.物探技術(shù)在固體廢棄物探測的應(yīng)用及前景展望[J]. 環(huán)境科學(xué)與技術(shù), 2006,29(12):93-95. LI Zhong, Nai C X, Ning S N, et al. Application and prospect of physical exploringtechnology for solid waste [J]. Environmental Science and Technology, 2006,29(12):93-95.

    [4] 李金銘.地電場與電法勘探[M]. 北京:地質(zhì)出版社, 2005:1-27. Li J M. Geoelectric field and electrical exploration [M]. Beijing: Geological Publishing House, 2005:1-27.

    [5] 孫亞坤,能昌信,劉玉強,等.鉻污染土壤電阻率特性及其影響因素研究[J]. 環(huán)境科學(xué)學(xué)報, 2011,31(9):1992-1998. Sun Y K, Nai C X, Liu Y Q, et al. Investigation on the electrical resistivity of chromium contaminated soil [J]. Acta Scientiae Circumstantiae, 2011,31(9):1992-1998.

    [6] 劉國華,王振宇,黃建平.土的電阻率特性及其工程應(yīng)用研究[J]. 巖土工程學(xué)報, 2004,26(1):83-87. Liu G H, Wang Z Y, Huang J P. Research on electrical resistivity feature of soil and it's application [J]. Chinese Journal of Geotechnical Engineering, 2004,26(1):83-87.

    [7] Delaney A, Peapples P, Arcone S. Electrical resistivity of frozen and petroleum-contaminated fine-grained soil [J]. Cold Regions Science & Technology, 2001,32(2/3):107-119.

    [8] 郭高山,李永濤,張 巖,等.垃圾填埋場污染土的電性與磁性響應(yīng)研究[J]. 中國環(huán)境科學(xué), 2015,39(9):2737-2744. Guo G S, Li Y T, Zhang Y, et al. The study of electric and magnetic response on polluted soil at waste landfill [J].China Environmental Science, 2015,39(9):2737-2744.

    [9] 管紹朋,王玉玲,能昌信,等.電學(xué)方法在雙襯層填埋場滲漏檢測中的應(yīng)用[J]. 中國環(huán)境科學(xué), 2011,31(12):2013-2017. Guan S P, Wang Y L, Nai C X, et al. Application of electrical leak detection method in double-lined landfills [J]. China Environmental Science, 2011,31(12):2013-2017.

    [10] 董 路,葉騰飛,能昌信,等.ERT技術(shù)在無機酸污染場地調(diào)查中的應(yīng)用[J]. 環(huán)境科學(xué)研究, 2008,21(6):67-71. Dong L, Ye T F, Nai C X, et al. Application of ERT technology to invest igate inorganic polluted sites [J].Research of Environmental Sciences, 2008,21(6):67-71.

    [11] Power C, Gerhard J I, Tsourlos P, et al. Improved time-lapse electrical resistivity tomography monitoring of dense non-aqueous phase liquids with surface-to-horizontal borehole arrays [J]. Journal of Applied Geophysics, 2015,112:1-13.

    [12] Liu Z B, Liu S Y, Cai Y, et al. Electrical resistivity characteristics of diesel oil-contaminated kaolin clay and a resistivity-based detection method [J]. Environmental Science & Pollution Research, 2015,22(11): 8216-8223.

    [13] 楊振威,許江濤,趙秋芳,等.復(fù)電阻率法(CR)發(fā)展現(xiàn)狀與評述[J]. 地球物理學(xué)進展, 2015,30(2):899-904. Yang Z W, Xu J T, Zhao Q F, et al. Current situation and review of complex resistivity [J]. Progress in Geophysics, 2015,30(2):899-904.

    [14] Revil A, Glover P W J. Theory of ionic-surface electrical conduction in porous media [J]. Physical Review B, 1997,55(3):1757-1773.

    [15] Kemna A, Binley A, Ramirez A, et al. Complex resistivity tomography for environmental applications [J]. Chemical Engineering Journal, 2000,77(1):11-18.

    [16] Revil A, Skold M. Salinity dependence of spectral induced polarization in sands and sandstones [J]. Geophysical Journal International, 2011,187(2):813-824.

    [17] 能昌信,劉玉強,劉豪睿,等.鉻污染土壤的導(dǎo)電性、頻譜激電性和介電特性的實驗結(jié)果[J]. 環(huán)境科學(xué), 2011,32(3):758-765. Nai C X, Liu Y Q, Liu H R, et al. Experiment results of conduction, spectral induced polarization and dielectric characteristics for chrome- contaminated soil [J].Environmental Science, 2011,32(3):758-765.

    [18] Schmutz M, Revil A, Vaudelet P, et al. Influence of oil saturation upon spectral induced polarization of oil-bearing sands [J]. Geophysical Journal International, 2010,183(1):211-224.

    [19] Grimm R E, Olhoeft G R, Mckinley K, et al. Nonlinear complex- resistivity survey for DNAPL at the Savannah River Site A-014outfall [J]. Journal of Environmental & Engineering Geophysics, 2005,10(4):351-364.

    [20] Orlando L, Renzi B. Electrical permittivity and resistivity time lapses of multiphase DNAPLs in a lab test [J]. Water Resources Research, 2015,51(1):377-389.

    [21] Revil A, Coperey A, Shao Z, et al. Complex conductivity of soils [J]. Water Resources Research, 2017,53(8):7121-7147.

    [22] Bérubé C, Chouteau M, Shamsipour P, et al. Bayesian inference of spectral induced polarization parameters for laboratory complex resistivity measurements of rocks and soils [J]. Computers & Geosciences, 2017,105:51-64.

    [23] Loke M, Chambers J, Ogilvy R. Inversion of 2D spectral induced polarization imaging data [J]. Geophysical Prospecting, 2010,54(3):287-301.

    [24] Revil A, Linde N. Chemico-electromechanical coupling in microporous media [J]. Journal of Colloid and Interface Science, 2006,302(2):682-694.

    [25] Yoon G L, Park J B. Sensitivity of leachate and fine contents on electrical resistivity variations of sandy soils [J]. Journal of Hazardous Materials, 2001,84(2):147-161.

    [26] Personna Y R, Slater L, Ntarlagiannis D, et al. Complex resistivity signatures of ethanol in sand-clay mixtures [J]. Journal of Contaminant Hydrology, 2013,149:76-87.

    [27] ASTM D1141-98(2003), Standard Practice for the Preparation of Substitute Ocean Water [S].

    [28] Lee J H, Oh M H, Park J, et al. Dielectric dispersion characteristics of sand contaminated by heavy metal, landfill leachate and BTEX (02- 104B) [J]. Journal of Hazardous Materials, 2003,105(1-3):83-102.

    [29] Kaya A, Fang Y H. Identification of contaminated soils by dielectric constant and electrical conductivity [J]. Journal of Environmentol Engineering, 1997,123(2):169-177.

    [30] Leroy P, Revil A, Kemna A, et al. Complex conductivity of water- saturated packs of glass beads [J]. Journal of Colloid and Interface Science, 2008,321(1):103-117.

    [31] Vaudelet P, Revil A, Schmutz M, et al. Induced polarization signatures of cations exhibiting differential sorption behaviors in saturated sands [J]. Water Resources Research, 2011,47(2):247-255.

    Complex resistivity of cationic metal contaminated sandy soils: Time-varying characteristics and formation mechanism.

    WANG Ze-ya1, XU Ya2*, DONG Lu2, NAI Chang-xin2, PAN Yong-tai1, LIU Yu-qiang2

    (1.School of Chemical & Environmental Engineering, China University of Mining & Technology, Beijing 100083, China;2.Research Institute of Soil and Solid Waste, Chinese Research Academy of Environment Sciences, Beijing 100012, China)., 2019,39(3):1147~1153

    In order to understand the time-varying characteristics and control mechanisms of complex resistivity in soils with different contamination degree, the experiment of sand column and water column was designed to measure the phase in 100~103Hz by taking NaCl, MnCl2and artificial seawater as contaminants. The results showed that the phase decreased with the pollutant concentration and increased with the frequency. The phase in sand column was more sensitive than that in water column with the change of pollutant concentration. Under the background of soil and water, when contaminated by 1g/L NaCl, the phase of soil decreased by about 42%, while the change in water column was almost invisible at 1Hz. The results of repeated measurements for fixed frequencies showed that the phase in water column was stable while it took on obvious time-varying property in sand column (>16Hz). This property also decreased with pollutant content in solution. At 102Hz, the background value of phase measurement in sand column was between -4.05~-11.67mrad, the variance reached 4.8 while it was only 0.09 in water column. Complex resistivity method has better performance in detection and monitoring of cations contaminated sites, the change of its parameters is mainly related to the content of cations inpore water.However, due to the interface properties and the complexity of the pore structure,migration and polarizationof charge carrier in the medium are chaoticby the alternating current, and finally cause phase fluctuation.

    phase;dielectric properties;joint inversion;geophysical prospecting

    X53

    A

    1000-6923(2019)03-1147-07

    王澤亞(1993-),男,江蘇泗洪人,中國礦業(yè)大學(xué)(北京)碩士研究生,主要從事固體廢物污染控制研究.

    2018-08-06

    國家重點研發(fā)計劃(2018YFC1800902);國家自然科學(xué)基金資助項目(51708529);中央級院所基本科研業(yè)務(wù)專項重點項目(2016YSKY14)

    * 責(zé)任作者, 副研究員, xuya@craes.org.cn

    猜你喜歡
    土柱極化電阻率
    降雨條件下植物修復(fù)分層尾礦土壤重金屬遷移的模擬分析
    認知能力、技術(shù)進步與就業(yè)極化
    分層土壤的持水性能研究
    雙頻帶隔板極化器
    電子測試(2017年15期)2017-12-18 07:18:51
    不同化學(xué)浸取劑對土壤鎘淋溶過程影響
    化工管理(2017年1期)2017-03-05 23:32:20
    不同灌水量對2種鹽堿土的洗鹽效果比較
    三維電阻率成像與高聚物注漿在水閘加固中的應(yīng)用
    基于PWM控制的新型極化電源設(shè)計與實現(xiàn)
    隨鉆電阻率測井的固定探測深度合成方法
    海洋可控源電磁場視電阻率計算方法
    欧美精品一区二区大全| 麻豆av在线久日| 久久久久久人人人人人| 热re99久久国产66热| 亚洲国产最新在线播放| 国产精品一区二区在线不卡| 制服人妻中文乱码| 中国国产av一级| 观看美女的网站| 精品卡一卡二卡四卡免费| 亚洲国产欧美网| 亚洲综合精品二区| 亚洲成人手机| 色婷婷久久久亚洲欧美| 久久久精品区二区三区| 妹子高潮喷水视频| 欧美亚洲日本最大视频资源| 99久久人妻综合| 精品亚洲乱码少妇综合久久| 国产在线一区二区三区精| 成年女人在线观看亚洲视频| 亚洲成av片中文字幕在线观看 | 久久精品久久久久久久性| 亚洲av综合色区一区| 女人高潮潮喷娇喘18禁视频| 久久久国产欧美日韩av| 两个人免费观看高清视频| 伊人久久大香线蕉亚洲五| 午夜影院在线不卡| 巨乳人妻的诱惑在线观看| 亚洲三级黄色毛片| 欧美日本中文国产一区发布| 你懂的网址亚洲精品在线观看| 亚洲经典国产精华液单| 欧美最新免费一区二区三区| 亚洲经典国产精华液单| 精品国产乱码久久久久久男人| 国产精品久久久久久av不卡| av天堂久久9| av在线播放精品| 我要看黄色一级片免费的| 亚洲一区中文字幕在线| 免费看av在线观看网站| 涩涩av久久男人的天堂| 少妇的逼水好多| 久久鲁丝午夜福利片| 亚洲av欧美aⅴ国产| 国产成人免费观看mmmm| 国产老妇伦熟女老妇高清| 老女人水多毛片| 欧美最新免费一区二区三区| 一本色道久久久久久精品综合| 午夜91福利影院| 国产深夜福利视频在线观看| 国产av码专区亚洲av| 国产精品秋霞免费鲁丝片| 高清不卡的av网站| 黄色毛片三级朝国网站| 黑人巨大精品欧美一区二区蜜桃| 成人免费观看视频高清| 18禁裸乳无遮挡动漫免费视频| 亚洲欧美一区二区三区久久| 亚洲欧美色中文字幕在线| 国产av码专区亚洲av| 男人操女人黄网站| 午夜激情久久久久久久| 欧美国产精品一级二级三级| 成人18禁高潮啪啪吃奶动态图| 两个人看的免费小视频| 久久久久久久精品精品| 久久国内精品自在自线图片| 亚洲国产色片| 飞空精品影院首页| tube8黄色片| 日韩一区二区视频免费看| av有码第一页| 在线看a的网站| 老鸭窝网址在线观看| 免费播放大片免费观看视频在线观看| 国产麻豆69| 多毛熟女@视频| 久久青草综合色| 国产老妇伦熟女老妇高清| 国产亚洲一区二区精品| 天天躁狠狠躁夜夜躁狠狠躁| 九九爱精品视频在线观看| 99久久人妻综合| 在线观看一区二区三区激情| 欧美人与善性xxx| 中文字幕另类日韩欧美亚洲嫩草| 日产精品乱码卡一卡2卡三| 狠狠婷婷综合久久久久久88av| 日韩一本色道免费dvd| 制服诱惑二区| 黄片无遮挡物在线观看| 欧美日韩视频精品一区| 欧美少妇被猛烈插入视频| 国产亚洲最大av| www.熟女人妻精品国产| av又黄又爽大尺度在线免费看| 亚洲综合色网址| 亚洲国产最新在线播放| 麻豆精品久久久久久蜜桃| 亚洲少妇的诱惑av| 美女午夜性视频免费| 国产毛片在线视频| 永久网站在线| 国产极品粉嫩免费观看在线| 成人漫画全彩无遮挡| 丰满乱子伦码专区| 成人漫画全彩无遮挡| 国产日韩一区二区三区精品不卡| 最近中文字幕高清免费大全6| 2021少妇久久久久久久久久久| 亚洲精品在线美女| 国产亚洲精品第一综合不卡| 不卡视频在线观看欧美| 亚洲欧美一区二区三区黑人 | 狂野欧美激情性bbbbbb| 女人被躁到高潮嗷嗷叫费观| 欧美在线黄色| 欧美激情极品国产一区二区三区| 男女啪啪激烈高潮av片| 亚洲成国产人片在线观看| 丰满乱子伦码专区| 中文字幕最新亚洲高清| 国产1区2区3区精品| 久久午夜福利片| 男人爽女人下面视频在线观看| 飞空精品影院首页| 精品视频人人做人人爽| 欧美人与善性xxx| 精品国产一区二区三区久久久樱花| 国产精品国产三级专区第一集| 久久久久精品久久久久真实原创| xxx大片免费视频| 久久精品国产综合久久久| 国精品久久久久久国模美| 亚洲国产精品成人久久小说| 亚洲伊人色综图| 久久久久精品人妻al黑| 电影成人av| 日韩大片免费观看网站| 国产人伦9x9x在线观看 | 黄片小视频在线播放| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 一边摸一边做爽爽视频免费| 亚洲欧美清纯卡通| 精品国产国语对白av| 亚洲精品av麻豆狂野| 国产成人91sexporn| 国产精品成人在线| 中文字幕另类日韩欧美亚洲嫩草| 成年动漫av网址| 免费黄网站久久成人精品| 在线 av 中文字幕| 色哟哟·www| 只有这里有精品99| 男女边摸边吃奶| 国产精品av久久久久免费| 国产精品久久久av美女十八| 人妻少妇偷人精品九色| 久久97久久精品| 成人二区视频| 黑人巨大精品欧美一区二区蜜桃| 老女人水多毛片| 亚洲经典国产精华液单| 你懂的网址亚洲精品在线观看| 亚洲美女视频黄频| 久久久久久久久免费视频了| 侵犯人妻中文字幕一二三四区| 美女午夜性视频免费| 成年美女黄网站色视频大全免费| 五月伊人婷婷丁香| 国产激情久久老熟女| 亚洲成人手机| 国产精品久久久久久久久免| 精品人妻在线不人妻| 久久精品人人爽人人爽视色| 边亲边吃奶的免费视频| 男人操女人黄网站| 久久久精品国产亚洲av高清涩受| 亚洲国产精品一区三区| 欧美日韩国产mv在线观看视频| 伊人久久大香线蕉亚洲五| 伦精品一区二区三区| 大话2 男鬼变身卡| 日韩av在线免费看完整版不卡| 日韩制服骚丝袜av| 免费观看无遮挡的男女| 国产精品国产三级国产专区5o| 久久久久久人人人人人| 又黄又粗又硬又大视频| 如何舔出高潮| 亚洲国产毛片av蜜桃av| 叶爱在线成人免费视频播放| 性少妇av在线| 在线观看免费高清a一片| 免费黄色在线免费观看| 午夜免费男女啪啪视频观看| 国产精品 国内视频| 国产亚洲午夜精品一区二区久久| 国产精品熟女久久久久浪| 丁香六月天网| 亚洲人成77777在线视频| 99热国产这里只有精品6| 中国国产av一级| 侵犯人妻中文字幕一二三四区| 欧美97在线视频| 久久人妻熟女aⅴ| 各种免费的搞黄视频| videossex国产| av线在线观看网站| 欧美最新免费一区二区三区| 久久久久精品久久久久真实原创| 久久国产亚洲av麻豆专区| 亚洲国产成人一精品久久久| 久久久久国产网址| 成人毛片a级毛片在线播放| 国产亚洲午夜精品一区二区久久| 成年人午夜在线观看视频| 伊人久久国产一区二区| 欧美日韩综合久久久久久| 一区二区av电影网| 国产精品久久久久久av不卡| 我要看黄色一级片免费的| av在线老鸭窝| 黄频高清免费视频| 精品久久蜜臀av无| 美女国产高潮福利片在线看| 国产色婷婷99| 日本-黄色视频高清免费观看| av免费在线看不卡| 卡戴珊不雅视频在线播放| 国产亚洲一区二区精品| 欧美老熟妇乱子伦牲交| 伊人亚洲综合成人网| 成人影院久久| 麻豆乱淫一区二区| 国产av一区二区精品久久| 熟女电影av网| 亚洲综合色惰| 一二三四中文在线观看免费高清| 97人妻天天添夜夜摸| av女优亚洲男人天堂| 观看美女的网站| 毛片一级片免费看久久久久| 国产精品免费视频内射| 新久久久久国产一级毛片| 久久久国产欧美日韩av| 欧美+日韩+精品| 你懂的网址亚洲精品在线观看| 日韩一本色道免费dvd| 国产黄频视频在线观看| 国产精品亚洲av一区麻豆 | 亚洲欧美精品自产自拍| 少妇的丰满在线观看| 亚洲一区中文字幕在线| 欧美+日韩+精品| 日本wwww免费看| 成人影院久久| 久久精品国产a三级三级三级| 夫妻性生交免费视频一级片| 精品午夜福利在线看| 亚洲五月色婷婷综合| 波多野结衣av一区二区av| 精品一区在线观看国产| 少妇 在线观看| 伦精品一区二区三区| 桃花免费在线播放| 少妇熟女欧美另类| 波野结衣二区三区在线| 可以免费在线观看a视频的电影网站 | 久久女婷五月综合色啪小说| 2018国产大陆天天弄谢| 麻豆精品久久久久久蜜桃| 国产白丝娇喘喷水9色精品| 日韩大片免费观看网站| 最近的中文字幕免费完整| 成年人免费黄色播放视频| av片东京热男人的天堂| 美国免费a级毛片| 日韩电影二区| 婷婷色综合www| 男女无遮挡免费网站观看| 久久99蜜桃精品久久| 久久精品久久精品一区二区三区| 欧美成人午夜免费资源| 国产福利在线免费观看视频| 国产精品 欧美亚洲| 免费在线观看视频国产中文字幕亚洲 | 国产亚洲午夜精品一区二区久久| 欧美日韩精品网址| 亚洲在久久综合| 亚洲一级一片aⅴ在线观看| 热99国产精品久久久久久7| 亚洲欧美成人精品一区二区| 搡女人真爽免费视频火全软件| 成年人午夜在线观看视频| 黄色一级大片看看| 一级毛片电影观看| 亚洲av在线观看美女高潮| a级毛片在线看网站| 色视频在线一区二区三区| 乱人伦中国视频| 91午夜精品亚洲一区二区三区| 亚洲精品av麻豆狂野| 午夜影院在线不卡| 香蕉国产在线看| 日韩欧美精品免费久久| 999久久久国产精品视频| 国产男女超爽视频在线观看| 久久久久国产一级毛片高清牌| 久久久精品94久久精品| 美国免费a级毛片| 日韩欧美一区视频在线观看| 亚洲激情五月婷婷啪啪| 久久韩国三级中文字幕| 在线观看三级黄色| 各种免费的搞黄视频| 黄色视频在线播放观看不卡| 99国产精品免费福利视频| 精品久久久久久电影网| 午夜激情av网站| 欧美激情高清一区二区三区 | 亚洲精品视频女| 美女xxoo啪啪120秒动态图| 亚洲av综合色区一区| 性高湖久久久久久久久免费观看| 亚洲一级一片aⅴ在线观看| 曰老女人黄片| 欧美日本中文国产一区发布| av片东京热男人的天堂| 青春草国产在线视频| 五月开心婷婷网| 亚洲成国产人片在线观看| 亚洲五月色婷婷综合| 亚洲av电影在线进入| www.精华液| 三上悠亚av全集在线观看| 欧美激情高清一区二区三区 | 91在线精品国自产拍蜜月| 丁香六月天网| av在线观看视频网站免费| 王馨瑶露胸无遮挡在线观看| 人妻 亚洲 视频| 最新中文字幕久久久久| 在现免费观看毛片| 捣出白浆h1v1| 亚洲成国产人片在线观看| 久久久精品区二区三区| 黄频高清免费视频| 综合色丁香网| 午夜影院在线不卡| 欧美日韩精品成人综合77777| 国产精品久久久久久精品电影小说| 欧美激情极品国产一区二区三区| 欧美精品国产亚洲| 啦啦啦啦在线视频资源| 亚洲欧美精品自产自拍| 久久精品久久精品一区二区三区| 久久精品熟女亚洲av麻豆精品| 最新中文字幕久久久久| 国产亚洲一区二区精品| 一级,二级,三级黄色视频| 咕卡用的链子| 国产伦理片在线播放av一区| 亚洲av男天堂| 热re99久久精品国产66热6| 制服丝袜香蕉在线| 久久人人爽av亚洲精品天堂| 一级爰片在线观看| 26uuu在线亚洲综合色| 天堂俺去俺来也www色官网| 亚洲欧美精品自产自拍| 啦啦啦在线免费观看视频4| 十八禁高潮呻吟视频| 亚洲人成网站在线观看播放| 午夜av观看不卡| 老司机亚洲免费影院| 捣出白浆h1v1| 97在线视频观看| 97在线人人人人妻| 久久精品亚洲av国产电影网| 亚洲国产最新在线播放| 在线精品无人区一区二区三| 中文精品一卡2卡3卡4更新| 亚洲国产成人一精品久久久| 男女边摸边吃奶| 一本大道久久a久久精品| 国产探花极品一区二区| 日本猛色少妇xxxxx猛交久久| 午夜精品国产一区二区电影| 日韩av在线免费看完整版不卡| 日韩免费高清中文字幕av| 久久久久国产网址| 亚洲欧洲精品一区二区精品久久久 | 日本午夜av视频| 韩国av在线不卡| 国产成人精品无人区| 大片免费播放器 马上看| 中文字幕人妻丝袜一区二区 | 在线观看国产h片| 天天躁狠狠躁夜夜躁狠狠躁| 热99久久久久精品小说推荐| 婷婷色综合大香蕉| 水蜜桃什么品种好| 男女免费视频国产| 免费在线观看视频国产中文字幕亚洲 | 91精品三级在线观看| 久久久国产一区二区| 女人被躁到高潮嗷嗷叫费观| 亚洲精品自拍成人| 夜夜骑夜夜射夜夜干| 日韩,欧美,国产一区二区三区| 国产亚洲最大av| 一个人免费看片子| 国产欧美日韩综合在线一区二区| 国产日韩欧美视频二区| 亚洲第一区二区三区不卡| 男男h啪啪无遮挡| 99国产精品免费福利视频| 亚洲欧洲精品一区二区精品久久久 | 亚洲精品久久午夜乱码| av视频免费观看在线观看| 亚洲精品乱久久久久久| 熟妇人妻不卡中文字幕| 欧美亚洲 丝袜 人妻 在线| 97人妻天天添夜夜摸| 男女免费视频国产| 夜夜骑夜夜射夜夜干| 亚洲人成77777在线视频| 高清不卡的av网站| 中文字幕人妻丝袜一区二区 | 99九九在线精品视频| 在线天堂中文资源库| 国产精品欧美亚洲77777| 天天影视国产精品| 国产精品秋霞免费鲁丝片| 久久精品久久久久久久性| 亚洲精品日韩在线中文字幕| 人妻少妇偷人精品九色| 最近最新中文字幕大全免费视频 | 久久久久久免费高清国产稀缺| 永久免费av网站大全| 精品卡一卡二卡四卡免费| 成年动漫av网址| 各种免费的搞黄视频| 亚洲国产欧美日韩在线播放| 精品一区二区免费观看| 黄片播放在线免费| 久热这里只有精品99| 桃花免费在线播放| 久久精品国产亚洲av天美| 亚洲熟女精品中文字幕| 亚洲欧美精品综合一区二区三区 | 国产一区亚洲一区在线观看| 成年av动漫网址| 国产日韩欧美在线精品| 日韩精品有码人妻一区| 丰满饥渴人妻一区二区三| 久久免费观看电影| www日本在线高清视频| 99国产综合亚洲精品| 丝袜美足系列| 欧美激情极品国产一区二区三区| 观看美女的网站| 亚洲,欧美精品.| 国产视频首页在线观看| 精品人妻在线不人妻| 精品午夜福利在线看| 亚洲欧美色中文字幕在线| 国产成人免费观看mmmm| 97在线人人人人妻| 精品久久久久久电影网| 亚洲伊人色综图| 大片电影免费在线观看免费| 欧美黄色片欧美黄色片| 大片免费播放器 马上看| 高清黄色对白视频在线免费看| 免费黄色在线免费观看| 视频区图区小说| 精品一区二区三卡| 在线观看美女被高潮喷水网站| 夫妻性生交免费视频一级片| 亚洲欧美一区二区三区国产| 五月伊人婷婷丁香| 在线观看人妻少妇| 国产精品 欧美亚洲| 99国产精品免费福利视频| 9191精品国产免费久久| 一级片免费观看大全| 性色av一级| 亚洲美女视频黄频| 日韩视频在线欧美| 国产综合精华液| 国产高清不卡午夜福利| 9热在线视频观看99| 青草久久国产| 乱人伦中国视频| 亚洲内射少妇av| 精品少妇黑人巨大在线播放| av国产久精品久网站免费入址| 亚洲av欧美aⅴ国产| 久热久热在线精品观看| xxx大片免费视频| a级片在线免费高清观看视频| 国产成人a∨麻豆精品| 国产成人精品无人区| 欧美亚洲日本最大视频资源| 国产一区亚洲一区在线观看| 欧美人与善性xxx| 最近手机中文字幕大全| 青草久久国产| 亚洲国产毛片av蜜桃av| 亚洲国产精品一区三区| 在线天堂最新版资源| 免费黄色在线免费观看| 精品久久久久久电影网| 亚洲视频免费观看视频| 99九九在线精品视频| 亚洲av中文av极速乱| 亚洲图色成人| 伦理电影免费视频| 大陆偷拍与自拍| 少妇人妻精品综合一区二区| 精品久久久久久电影网| 亚洲精品一二三| 精品国产国语对白av| 91午夜精品亚洲一区二区三区| 亚洲美女搞黄在线观看| 欧美另类一区| 老司机影院毛片| 国产精品免费视频内射| 日韩制服骚丝袜av| 免费看不卡的av| 9191精品国产免费久久| 2022亚洲国产成人精品| 亚洲伊人色综图| 久久ye,这里只有精品| 成人亚洲欧美一区二区av| 久久精品夜色国产| 亚洲国产av新网站| 18禁国产床啪视频网站| 最近手机中文字幕大全| 制服诱惑二区| 亚洲国产看品久久| 国产av国产精品国产| 国产福利在线免费观看视频| 狠狠精品人妻久久久久久综合| 寂寞人妻少妇视频99o| 在线观看www视频免费| 各种免费的搞黄视频| 久久精品国产亚洲av涩爱| 五月天丁香电影| 精品国产乱码久久久久久小说| 久久这里有精品视频免费| 一二三四中文在线观看免费高清| av网站免费在线观看视频| www.自偷自拍.com| 国产精品一区二区在线不卡| 肉色欧美久久久久久久蜜桃| 婷婷色av中文字幕| 亚洲精品国产一区二区精华液| 国产欧美亚洲国产| 亚洲一区中文字幕在线| 成人国产麻豆网| kizo精华| 2021少妇久久久久久久久久久| 欧美精品一区二区免费开放| 久久国产精品男人的天堂亚洲| av国产精品久久久久影院| 母亲3免费完整高清在线观看 | 成人黄色视频免费在线看| 免费少妇av软件| 日日摸夜夜添夜夜爱| 视频在线观看一区二区三区| 国产亚洲av片在线观看秒播厂| 9色porny在线观看| 久久99一区二区三区| 超碰97精品在线观看| 免费久久久久久久精品成人欧美视频| 色94色欧美一区二区| 日韩中文字幕欧美一区二区 | 久久这里只有精品19| 久久亚洲国产成人精品v| 99久国产av精品国产电影| av福利片在线| 国产又色又爽无遮挡免| 美国免费a级毛片| av免费观看日本| 80岁老熟妇乱子伦牲交| 国产在视频线精品| 国产 一区精品| 国产淫语在线视频| 久久午夜福利片| 黄片无遮挡物在线观看| 国产淫语在线视频| 欧美日韩国产mv在线观看视频| 如日韩欧美国产精品一区二区三区| 久久精品熟女亚洲av麻豆精品| 精品卡一卡二卡四卡免费| 香蕉国产在线看| 亚洲精品久久成人aⅴ小说| 丰满少妇做爰视频| 两性夫妻黄色片| 美女午夜性视频免费| 亚洲婷婷狠狠爱综合网| 欧美亚洲 丝袜 人妻 在线| 日韩视频在线欧美| 国产免费一区二区三区四区乱码| 国产精品无大码| 久久这里有精品视频免费| 日韩免费高清中文字幕av| 99热国产这里只有精品6|