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

    重塑黃土核磁共振T2譜的定量分析與應(yīng)用

    2022-06-17 00:53:21曹雪梅詹永康袁磊舒中潘張軍
    人民長江 2022年5期
    關(guān)鍵詞:重塑黃土滲透率

    曹雪梅 詹永康 袁磊 舒中潘 張軍

    摘要:核磁共振技術(shù)(NMR)憑借其便捷性、低成本、無損傷的優(yōu)勢被廣泛應(yīng)用于黃土孔隙空間研究領(lǐng)域中。但目前對于重塑黃土核磁共振T2譜的分析應(yīng)用成果過于單一,從而造成了一定程度上的數(shù)據(jù)浪費。有鑒于此,基于核磁共振技術(shù),利用數(shù)學(xué)分析手段對不同干密度重塑黃土的T2譜開展了定量分析;將分析得到的數(shù)學(xué)參數(shù)應(yīng)用于石油領(lǐng)域中的SDR模型,并構(gòu)建了重塑黃土的滲透率計算模型。結(jié)果表明:① 不同干密度重塑黃土的核磁共振T2譜在數(shù)學(xué)分析參數(shù)(均值、幾何均值、標準差以及峰度)中存在不同程度的差異,反映到物理意義上,則是重塑黃土干密度會對土體內(nèi)部孔隙大小、分布及分選特征產(chǎn)生影響。相比傳統(tǒng)的方法僅獲取土體孔隙體積,數(shù)學(xué)統(tǒng)計方法獲得的數(shù)學(xué)參數(shù)對孔隙的描述更全面也更精確。② 借助于構(gòu)建的重塑黃土的SDR模型,求解了試驗黃土的滲透率,計算的最小滲透率約為0.3 mD,計算結(jié)果與前人總結(jié)的規(guī)律接近。③ 干密度會影響多個土體微觀空間特征,從而引起土體滲透特性的改變。④ 從相關(guān)關(guān)系方面的分析來看,T2譜的均值和幾何均值對土體滲透率影響較大。研究成果在一定程度上可為核磁共振技術(shù)在黃土領(lǐng)域中的應(yīng)用提供新的思路,可為黃土滲透性的預(yù)測提供新的方法,也可為了解土體滲透機制提供新的思考。

    關(guān) 鍵 詞:重塑黃土; 孔隙空間; 核磁共振技術(shù)(NMR); 微觀特征; 定量參數(shù)

    中圖法分類號: P66

    文獻標志碼: A

    DOI:10.16232/j.cnki.1001-4179.2022.05.033

    0 引 言

    在中國西北地區(qū),黃土極為重要,分布面積達63萬km2,主要分布在山西、陜西、甘肅和寧夏地區(qū)[1]。典型黃土是一種具有開放和亞穩(wěn)定結(jié)構(gòu)的風(fēng)沙沉積[2]。一些黃土表現(xiàn)出水敏性,在濕潤[3]時,由于內(nèi)部結(jié)構(gòu)塌陷趨于不可控沉降,黃土地基塌陷等巖土工程問題與壓實黃土的水分滲流特性密切相關(guān),而作為流體在壓實黃土中的流動通道,孔隙空間起著至關(guān)重要的作用[4-6]。因此,研究重塑黃土微觀孔隙特征,是了解工程中黃土宏觀滲透表現(xiàn)的基礎(chǔ)。

    核磁共振技術(shù)(NMR)憑借其便捷性、低成本、無損傷的優(yōu)勢,被廣泛應(yīng)用于研究分析巖土體孔隙空間結(jié)構(gòu)的巖土工程領(lǐng)域,其精確性與穩(wěn)定性均得到了研究人員的驗證。田慧會等[7]利用自由水與結(jié)合水冰點不同的特性獲取了T2截止值,通過T2截止值得到了T2譜代表的結(jié)合水含量。Yuan等[8]通過核磁共振T2譜,研究了不同干密度和干濕循環(huán)條件下壓實黃土的滲透特性及結(jié)構(gòu)演化。還有大量學(xué)者通過核磁共振技術(shù)研究了不同巖土體的孔隙特征,進而分析了土體某些宏觀表現(xiàn)行為,獲得了許多科研成果,效果良好[9-11]。

    T2譜作為核磁共振技術(shù)的重要成果,可以反映土體內(nèi)孔隙分布、組成等特點,通過對試樣T2譜的特征進行定性與定量分析,可進一步了解重塑黃土孔隙的空間結(jié)構(gòu)。有學(xué)者將T2譜與壓汞試驗相結(jié)合,從而提出了非飽和土的滲透系數(shù)的預(yù)測方法[12];還有學(xué)者對T2譜進行了數(shù)學(xué)分析,將數(shù)學(xué)參數(shù)應(yīng)用于石油開采率上,從新的角度了解了砂巖儲油層中孔隙空間特征[11]。

    針對黃土的孔隙結(jié)構(gòu)空間,已有部分學(xué)者利用核磁共振技術(shù)對其開展了較多的研究[13-14],但對于T2譜的應(yīng)用往往僅局限于獲取土體孔隙直徑大小的分布上,并未對T2譜本身所含有的其他空間信息進行進一步提取,更未進行定量分析,成果應(yīng)用過于單一,這就使得核磁共振T2譜在黃土領(lǐng)域中的應(yīng)用不夠充分,而且在一定程度上也造成了數(shù)據(jù)的浪費。

    本文基于核磁共振技術(shù)獲取了不同干密度重塑黃土的T2譜,并利用數(shù)學(xué)分析手段對其開展了定量分析,并將定量分析得到的數(shù)學(xué)參數(shù)應(yīng)用到求取不同干密度重塑黃土的滲透率計算模型的構(gòu)建中。該方法在一定程度上可為核磁共振技術(shù)在黃土領(lǐng)域中的應(yīng)用提供一個新的思路,也可為黃土滲透性的計算提供新的方法;同時,也有助于了解土體的滲透機制。

    1 NMR技術(shù)原理

    核磁共振(NMR)技術(shù)是基于H原子的自旋磁性和外加磁場之間的相互作用對艙內(nèi)試樣中的水進行檢測的。信號大小與氫原子的個數(shù)成正比,與橫向磁化矢量的衰減成正比。核磁共振弛豫時間與試樣的孔隙半徑成正比,這是因為大孔中儲水量大,而微孔中儲水量小。核磁共振是指在外部磁場作用下能級之間的共振躍遷[15]。在這種磁場下,可以得到偏轉(zhuǎn)后的質(zhì)子數(shù)和質(zhì)子自轉(zhuǎn)軸恢復(fù)到初始平衡態(tài)所需要的時間(T2)[16]。

    多孔介質(zhì)的孔隙內(nèi)部結(jié)構(gòu)、孔隙中的分子運動及其與黃土的反應(yīng)過程是多孔介質(zhì)研究領(lǐng)域的重要課題。對于土壤孔隙中的流體,存在3種類型的弛豫機制,即表面流體和自由流體的弛豫、分子的自擴散。因此,孔隙流體流動的弛豫時間T2可表示如下[17]:

    1T2=1T2f+1T2s+1T2d(1)

    式中:T2為核磁共振橫向弛豫時間,ms;T2f為大容器中孔隙流體的弛豫時間,ms;T2s和T2d分別為表面流體弛豫引起的孔隙流體弛豫時間和磁場梯度(ms)擴散引起的孔隙流體弛豫時間。

    1T2s=ρ2(SV)(2)

    式中:ρ2是橫向弛豫系數(shù),取決于土壤性質(zhì)[17];S/V是孔隙表面積與流體體積的比值。

    當只考慮飽和流體時,自由弛豫可以忽略不計,其值比表面弛豫小得多。同時,當磁場均勻且回波時間間隔足夠短時,擴散弛豫也可以忽略不計。因此,流體的松弛主要來源于土壤顆粒的表面松弛,式(2) 可近似為[18]

    1T2=ρ2(SV)=Fsρ2R(3)

    當孔隙為球形時,T2與孔隙半徑有關(guān),F(xiàn)s為孔隙形狀因子(當孔隙為圓柱形時Fs=3)。因此,通過式(3)可知核磁共振T2譜線能夠反演出土體內(nèi)部的孔隙空間結(jié)構(gòu)與大小分布特征[13]。

    2 試驗過程

    2.1 試樣準備

    本次試驗所用黃土取自陜西省延安市某挖方邊坡,其基本物理性質(zhì)如表1所列。

    依據(jù)GB/T 50123-2019《土工試驗方法標準》,將土樣在108 ℃下烘干24 h,然后碾細后過2 mm篩備用。正如表1所列,試驗所用黃土最大干密度為1.8 g/cm3,最優(yōu)含水量為16.3%,據(jù)此本文設(shè)置了5種不同干密度(1.80,1.74,1.68,1.62,1.56 g/cm3)的重塑黃土樣,為了減少試驗偶然誤差,另設(shè)置1組平行試驗。重塑黃土試樣含水量均為最優(yōu)含水量:16.3%。

    試樣配制時,為了控制制樣均勻,首先,將定量散狀篩后黃土與相應(yīng)的蒸餾水混合攪拌均勻,控制含水量為16.3%,用保鮮膜密封靜置24 h,待土中水分分布均勻。然后,根據(jù)目標干密度計算得到單個試樣所需的濕土質(zhì)量,先后分3層擊實進入環(huán)刀內(nèi),每層土間刮毛防止分層現(xiàn)象。由于核磁共振中磁場較強,傳統(tǒng)鋼制環(huán)刀對其影響較大,本試驗中所用環(huán)刀為特制聚四氟乙烯材質(zhì)環(huán)刀,聚四氟乙烯可有效消除金屬對NMR信號的干擾。按照上述步驟,依次完成了2組5個干密度重塑黃土試樣。

    2.2 試驗過程

    采用核磁共振掃描法測量重塑黃土的T2圖譜,所用核磁共振儀器型號為MacroMR12-150H-I。該核磁共振儀永久磁鐵磁場強度為(0.3±0.05)T,射頻發(fā)射功率的峰值輸出大于300 W,最大線性失真小于0.5%。本次試驗所用到的線圈探頭的直徑為60 mm。

    在核磁共振中,要得到能完整反映土體空間結(jié)構(gòu)的T2譜線,試樣內(nèi)孔隙必須要先飽水。所以,首先對聚四氟乙烯內(nèi)的土樣進行飽和,將試樣兩端用透水石封住,防止土體孔隙在飽和過程中發(fā)生改變,也可使飽和過程更加均勻。然后將試樣置于真空飽和器中,抽取真空2 h后再進行真空飽和12 h,飽和后立刻進行NMR常溫掃描,以獲取2組不同干密度重塑黃土的T2圖譜。圖1為飽和中的試樣。

    3 核磁共振T2譜圖形分布特征

    在核磁共振對土體的測試中,首先對試樣施加磁場,讓土體孔隙中的水的H原子產(chǎn)生共振,然后撤銷磁場。由于H原子從共振的狀態(tài)恢復(fù)到激發(fā)前的狀態(tài)時會發(fā)出一定的電信號,孔隙越大,恢復(fù)到激發(fā)前狀態(tài)的橫向弛豫時間也就越久。同一個孔徑大小的孔隙越多,體積就越多,產(chǎn)生的電信號也就越強。利用傳感器收集這些電信號,然后經(jīng)過數(shù)學(xué)反演就得到了T2譜線,橫坐標為弛豫時間,其物理意義是孔徑大小,弛豫時間越大,代表的孔隙越大??v坐標是信號強度,某個大小的孔隙越多,其信號強度越強,其物理意義為某個大小孔隙的體積。

    接下來,再將縱坐標的信號強度進行歸一化,即用單個弛豫時間下的信號強度去除所有弛豫時間下的信號強度的總和,這樣就將信號強度轉(zhuǎn)化為了孔隙度分量。最大孔隙度分量約為3%。

    在核磁共振試驗中,不同干密度樣品的橫向弛豫時間T2分布曲線均為雙峰分布(見圖2),弛豫時間短則代表了小孔隙,弛豫時間長則代表了大孔隙。雙峰峰型分布表明重塑黃土體內(nèi)孔隙具有宏觀、中、微觀的特征。

    根據(jù)T2譜的特征可分為圖形分布特征與數(shù)學(xué)分布特征,圖形分布特征可用:最大最小弛豫時間、譜峰弛豫時間、區(qū)間孔隙分量、最大孔隙度分量和歪度這幾個特征參數(shù)進行描述,下面依次分析其圖形主要特征。

    由于2組試樣為平行試驗,以下參數(shù)分析均取這2組試驗的平均值。

    3.1 最大最小弛豫時間

    T2譜中的弛豫時間越大,其代表的孔隙也就越大,因此,最大最小弛豫時間代表了重塑黃土內(nèi)部最大與最小的孔隙。由圖2可知:不同干密度的重塑黃土最小弛豫時間基本一致,為0.01 ms,這是由于儀器本身的限制,所以最小弛豫時間沒有對比意義。最大弛豫時間隨干密度的增大而減?。ㄒ妶D3),呈現(xiàn)出非線性的趨勢。

    3.2 半弛豫時間

    半弛豫時間指的是累計孔隙度達到50%時的橫向弛豫時間,該參數(shù)可在一定程度上反映土體中孔隙比的大小。如表2所列,試樣的半弛豫時間隨干密度的減小而增大,兩者呈負相關(guān)的關(guān)系。這也表明,隨著干密度的減小,土體內(nèi)總孔隙體積與孔隙比在不斷的增加。

    3.3 譜峰弛豫時間

    譜峰弛豫時間指的是在核磁共振T2譜中信號強度最大的弛豫時間,代表了孔隙大小最為集中的弛豫時間分量。不同干密度試樣譜峰弛豫時間基本相差不大,說明重塑黃土在弛豫時間在0.73 ms左右所代表的孔隙孔徑最為集中,與試樣的密度無明顯關(guān)系(見表3)。

    3.4 區(qū)間孔隙分量

    核磁共振T2譜線的橫軸弛豫時間代表了土體中不同大小的孔隙,孔隙越大,其弛豫時間越大,而譜線縱軸的信號強度則可以反映該孔隙的體積大小。將不同大小的孔隙分成大中小3類,3類孔隙的體積分別在全孔隙體積中所占的比例被稱為區(qū)間孔隙分量。

    根據(jù)譜線的形態(tài)將孔隙分為3類(見圖4):最小弛豫時間至1.80 g/cm3試樣半弛豫時間(0.79 ms),是小孔隙S1;1.80 g/cm3試樣半弛豫時間(0.79 ms)至次譜峰弛豫時間(11.1 ms),為中孔隙S2;次譜峰弛豫時(11.1 ms)至最大弛豫時間(1 000 ms)為大孔隙S3。

    T2譜與橫軸圍成的面積(無量綱)可代表不同大小孔隙的相對體積,所以通過統(tǒng)計3個區(qū)間(S1、S2、S3)的譜峰面積,可以計算得到不同大小的孔隙的相對體積,再分別計算出3類孔隙體積所占比例,計算結(jié)果如圖5所示。

    圖5中,S2所代表的中孔隙比例隨干密度的變化影響不大,而S1代表的小孔隙,S3代表的大孔隙隨干密度的增大分別呈正線性相關(guān)和負線性相關(guān),這表明干密度對于土體內(nèi)部孔隙大小分布影響巨大,干密度越大的土體孔隙越少,孔隙也越集中。

    4 核磁共振T2譜數(shù)學(xué)分布特征

    在傳統(tǒng)地質(zhì)數(shù)學(xué)統(tǒng)計模型中,最常用的概率統(tǒng)計模型是正態(tài)分布統(tǒng)計。前人應(yīng)用正態(tài)概率分布里的數(shù)學(xué)特征參數(shù)比如均值、方差、中值、變異系數(shù)等,去描述各類巖石與土體的孔隙大小分布情況;而T2譜作為土體內(nèi)部孔隙大小及分布的直接物理體現(xiàn)。因此,T2曲線也能使用數(shù)學(xué)統(tǒng)計參數(shù)來對其進行數(shù)學(xué)分布特征的定量分析。

    4.1 均 值

    均值是數(shù)學(xué)統(tǒng)計中非常常用的數(shù)學(xué)位置特征參數(shù),可由核磁共振T2譜與孔隙度加權(quán)平均求得[11]:

    T2=Ni=1T2ii(4)

    式中:i為第i個弛豫時間對應(yīng)的核磁共振孔隙度分量;

    T2i為第i個弛豫時間,ms。

    由式(4)可求得各重塑黃土試樣的T2譜的均值,能反映全孔隙中的平均孔隙大小。由表4可以看出:試樣T2譜均值隨干密度的減小而增大,呈負線性相關(guān)關(guān)系。

    4.2 幾何均值

    幾何均值和均值都是求一組數(shù)據(jù)中平均值的方法[19],而幾何均值則更多用于數(shù)據(jù)的平均增長率的研究中。幾何均值可由下式定義[20]:

    T2gm=(∏Ni=1T2i)i(5)

    式中:T2gm為T2的幾何均值;

    i為第i個弛豫時間對應(yīng)的核磁共振孔隙度分量;

    T2i為第i個弛豫時間,ms。

    按式(5)求得每組試樣的幾何均值。各試樣T2譜幾何均值如表5所列。

    幾何均值隨干密度的增大而減小,呈負相關(guān)關(guān)系,表明重塑黃土的干密度越大,T2曲線的平均增長率越小。

    4.3 標準差

    標準差是描述離均差平方的算術(shù)平均數(shù)。在孔隙系統(tǒng)中,可用于描述孔隙體積大小的分選性,標準差越小,說明孔隙大小的分選性越好,孔隙大小越集中。可由下式計算:

    σ=(iT2i-T2)2(6)

    式中:σ為標準差;

    i為第i個弛豫時間對應(yīng)的核磁共振孔隙度分量;

    T2i為第i個弛豫時間,ms;

    T2為均值;

    為均值對應(yīng)的孔隙度分量。

    按式(6)求得每組試樣的標準差。各試樣T2譜標準差如表6所列。

    標準差隨干密度的增大而減小,呈負相關(guān)關(guān)系,表明干密度越大的重塑黃土,其內(nèi)部孔隙分選性越好,孔隙大小越集中在某個孔徑區(qū)間。

    4.4 變異系數(shù)

    標準差與均值的比值被稱為變異系數(shù),是對數(shù)據(jù)值相對變化的體現(xiàn)。變異系數(shù)越小,孔隙結(jié)構(gòu)、分選性也就越好。其計算式為

    C=σ/T2(7)

    按式(7)求得每組試樣的變異系數(shù)。各試樣T2譜變異系數(shù)如表7所列。

    從表7可以得知:不同干密度試樣核磁共振T2譜的變異系數(shù)基本一致,說明干密度的大小并不會影響重塑黃土內(nèi)部變異系數(shù)。

    4.5 峰 度

    峰度在核磁共振中是用來衡量T2譜曲線的陡峭程度,可以反映重塑黃土內(nèi)部不同大小孔隙的集中程度,是用來評價孔隙結(jié)構(gòu)很好的一個定量參數(shù),計算公式如下:

    KG=(2iT2i-T2)4σ4(8)

    按式(8)求得每組試樣的峰度。各試樣T2譜峰度如表8所列。

    由表8可以看出:峰度大致隨干密度增大而增大,表明干密度越大,T2曲線越陡,孔隙越集中。

    5 定量參數(shù)分析與應(yīng)用

    5.1 定量參數(shù)分析

    前文對不同干密度重塑黃土的核磁共振T2譜進行了定性與定量研究,得到了各類豐富的定量參數(shù)與圖形信息,各類參數(shù)在T2譜中所表達的物理意義如圖6所示。這些特征參數(shù)對于描述孔隙空間也更為全面,充分挖掘了T2譜中蘊含的物理及數(shù)學(xué)信息。通過對每組試樣的核磁共振T2譜的分解,研究了不同干密度所導(dǎo)致重塑黃土孔隙結(jié)構(gòu)空間的異同,更進一步地了解了重塑黃土微觀組構(gòu)。

    通過對T2譜的圖形分析與定量表征,核磁共振特征參數(shù)可根據(jù)其代表的物理意義的不同而劃分為以下3類。

    (1) 均值、幾何均值、半弛豫時間、譜峰弛豫時間以及最大弛豫時間,這類參數(shù)可表征重塑黃土的孔隙體積的大小;

    (2) 大中小孔隙的區(qū)間孔隙分量,可代表每種孔隙以及孔喉的分布范圍與占比,與流體運動直接相關(guān);

    (3) 變異系數(shù)、標準差,則表達了重塑黃土孔隙分布的分選特征與分布規(guī)律,間接影響了水在重塑黃土中的滲流過程。

    圖7展示了重塑黃土的數(shù)學(xué)特征參數(shù)與圖像特征參數(shù)隨干密度的變化而變化的分布特征,這些參數(shù)直觀地描述了試樣內(nèi)部孔隙空間特征的具體差異,并將這種差異大小進行了定量,有利于對孔隙空間特征進行下一步的分析與應(yīng)用。由圖7可以發(fā)現(xiàn):在圖形特征參數(shù)中,譜峰弛豫時間沒有太大差異,而半弛豫時間卻隨干密度的增大而明顯減小,這表明干密度越大,孔隙更加集中,大孔隙也越少,這種差異在區(qū)間孔隙分量這一參數(shù)中得以體現(xiàn)。在數(shù)學(xué)統(tǒng)計分析參數(shù)中,均值、幾何均值、峰度以及標準差均有顯著變化,而這些差異可以體現(xiàn)在下文對重塑黃土滲透率的影響中。

    5.2 定量參數(shù)在滲透率模型求取中的應(yīng)用

    作為流體在土體中的流動通道,孔隙的大小、形態(tài)、分布等空間特征對于土體滲透性的影響巨大。前文獲取到了可定量表征重塑黃土內(nèi)孔隙空間特征差異的數(shù)學(xué)參數(shù),這節(jié)將這些參數(shù)應(yīng)用到土體滲透率計算模型的構(gòu)建中,可為理解土體滲透機制提供新的思路。

    目前,在黃土領(lǐng)域利用核磁共振構(gòu)建滲透率模型的文獻偏少或幾乎沒有,而在石油領(lǐng)域,利用核磁共振構(gòu)建滲透率模型,有多種構(gòu)建方法,Coates經(jīng)典模型的建立需要求得土體中的束縛水體積,而另一種由Klenberg[21]提出的SDR模型則需要求取T2譜的幾何均值。綜合前文獲取的數(shù)學(xué)參數(shù),本研究在SDR模型基礎(chǔ)上,引入在石油儲層滲透率模型中應(yīng)用較好的SDR模型[11]來對重塑黃土的滲透率計算模型進行研究,該計算模型公式如下:

    K=0.01T2gmσT2hi,maxKG(9)

    式中:K為滲透率,mD;

    為孔隙度,%;

    σ為標準差;

    i,max為最大孔隙度分量,%;

    T2h為譜峰弛豫時間,ms;

    T2gm為幾何均值,ms;

    Kg為峰度。

    式(9)中,每組試樣孔隙度可根據(jù)土的三相公式換算得到,其他參數(shù)在前文均已計算得到,對不同干密度重塑黃土的滲透率模型進行構(gòu)建,根據(jù)計算公式得到每組試樣的滲透率,計算結(jié)果如表9所列。

    重塑黃土的干密度與滲透率呈現(xiàn)良好的負指數(shù)相關(guān)關(guān)系(見圖6),擬合相關(guān)性好,與前人對重塑黃土滲透性的研究[22-23]基本一致,一般認為干密度是對重塑黃土滲透性最大的因素,因為干密度會在很大程度上改變土體內(nèi)部孔隙空間結(jié)構(gòu),導(dǎo)致流體在流動通道內(nèi)的運動過程發(fā)生改變,從而引起滲透特性的變化。

    6 結(jié) 論

    基于核磁共振技術(shù),獲取了不同干密度重塑黃土的T2譜;并利用數(shù)學(xué)分析手段對其進行定量分析,將定量分析得到的數(shù)學(xué)參數(shù)應(yīng)用于構(gòu)建重塑黃土的滲透率計算模型中,探索了一套T2譜在黃土孔隙空間分析中應(yīng)用的新思路??梢缘贸鲆韵轮饕Y(jié)論。

    (1) 利用數(shù)學(xué)統(tǒng)計手段對核磁共振T2譜進行定量分析,可以更直觀與更精確地分析不同試樣間的孔隙空間差異。不同干密度重塑黃土的核磁共振T2譜在數(shù)學(xué)統(tǒng)計參數(shù)中存在不同程度的差異,這些差異反映到物理意義上則是:重塑黃土干密度不同會對土體內(nèi)部孔隙大小、分布及分選特征產(chǎn)生影響,這些影響的大小則均可通過參數(shù)的數(shù)值差異進行更為準確的定量表達。相比傳統(tǒng)僅獲取土體孔隙體積,數(shù)學(xué)統(tǒng)計方法獲得的數(shù)學(xué)參數(shù)對孔隙的描述更全面與精確。

    (2) 將石油領(lǐng)域中SDR模型引入到了重塑黃土的滲透率模型的構(gòu)建中,得到了重塑黃土的滲透率模型;并根據(jù)模型求取了不同干密度重塑黃土的滲透率。根據(jù)模型計算得到試樣黃土的滲透率為0.30~3.81 mD。其次,干密度會導(dǎo)致土體滲透特性的改變:干密度越大,滲透率越低,且呈負指數(shù)的相關(guān)關(guān)系。另外,模型還反映了不同微觀孔隙特征參數(shù)對土體滲透率的貢獻大小,進一步理解了干密度變化導(dǎo)致土體滲透率變化的微觀因素。模型計算結(jié)果與前人獲得規(guī)律相似度高,這表明本文研究方法可信度較高,可行性較好。

    本文研究成果在一定程度上可為核磁共振在黃土領(lǐng)域中的應(yīng)用提供新的思路,可為黃土滲透性的預(yù)測提供新的方法,更可為了解土體滲透機制提供新的思考,具有較高的理論創(chuàng)新性與應(yīng)用實用性。

    參考文獻:

    [1] 王數(shù).地質(zhì)學(xué)與地貌學(xué)[M].北京:中國農(nóng)業(yè)大學(xué)出版社,2013.

    [2] ROLLINS K M,ROGERS G M.Mitigation measures for small structures on collapsible alluvial soils[J].Journal of Geotechnical Engineering,1994,120(9):1533-1553.

    [3] ASSALLAY A M,ROGERS C D F,SMALLEY I J.Formation and collapse of metastable particle packings and open structures in loess deposits[J].Engineering Geology,1997,48(1-2):101-115.

    [4] SIRVITHAYAPAKORN S,KELLER A.Transport of colloids in saturated porous media:A pore-scaleobservation of the size exclusion effect and colloid acceleration[J].Water Resources Resarch,2003,39(4):1109-1120.

    [5] 李鵬舉,谷宇峰.核磁共振T2譜轉(zhuǎn)換偽毛管壓力曲線的矩陣方法[J].天然氣地球科學(xué),2015,26(4):700-705.

    [6] 劉勇健,李彰明,郭凌峰,等.基于核磁共振技術(shù)的軟土三軸剪切微觀孔隙特征研究[J].巖石力學(xué)與工程學(xué)報,2018,37(8):1924-1932.

    [7] 田慧會,韋昌富.基于核磁共振技術(shù)的土體吸附水含量測試與分析[J].中國科學(xué):技術(shù)科學(xué),2014,44(3):295-305.

    [8] YUAN K Z,LU X F,WANG X J.Permeability characteristics and structural evolution of compacted loess under different dry densities and wetting-drying cycles[J].PLoS ONE,2021,16(6):e0253508.

    [9] 運華云,趙文杰,劉兵開,等.利用T2分布進行巖石孔隙結(jié)構(gòu)研究[J].測井技術(shù),2002,26(1):18-21,89.

    [10] 李海波,朱巨義,郭和坤.核磁共振T2譜換算孔隙半徑分布方法研究[J].波譜學(xué)雜志,2008,25(2):273-280.

    [11] 白松濤,程道解,萬金彬,等.砂巖巖石核磁共振 T2譜定量表征[J].石油學(xué)報,2016,37(3):282-391.

    [12] 陳銀.基于核磁共振技術(shù)的非飽和土滲透系數(shù)預(yù)測方法[D].武漢:湖北工業(yè)大學(xué),2019.

    [13] 何攀,許強,劉佳良,等.基于核磁共振技術(shù)的結(jié)合水含量對重塑黃土抗剪強度影響試驗研究[J].山地學(xué)報,2020,38(4):571-580.

    [14] 楊更社尤,吳迪.凍融環(huán)境下原狀黃土孔徑分布與其力學(xué)特性關(guān)系的試驗研究[J].煤炭工程,2019,51(3):107-112.

    [15] SCHAUMANN G E.Proton nuclear magnetic resonance(NMR)relaxometry in soil science[J].Encyclopedia of Agrophysics,2011:543-669.

    [16] CONTE P,STEFANO C D,F(xiàn)ERRO V,et al.Assessing hydrological connectivityinside a soil by fast-fieldcycling nuclear magnetic resonance relaxometry and its link to sediment delivery processes[J].Environ.Earth Sci.,2017,76:526.

    [17] XU J,LI Y,REN C,et al.Damage of saline intact loess after dry-wet and its interpretation based on SEM and NMR[J].Soils and Foundations,2020,60(4):911-928.

    [18] 何雨丹,毛志強,肖立志.核磁共振T2分布評價巖石孔徑分布的改進方法[J].地球物理學(xué)報,2005,48(2):373-378.

    [19] 何更生,唐海.油層物理[M].2 版.北京:石油工業(yè)出版社,2011.

    [20] 楊正明,李治碩,王學(xué)武,等.特低滲透油田相對滲透率曲線測試新方法[J].石油學(xué)報,2010,31(4):629-632.

    [21] KLEINBERG R L.Utility of NMR T2 distributions,connection with capillary pressure,clay effect,and determination of the surface relaxivity parameter ρ2[J].Magnetic Resonance Imaging,1996,14(7-8):761-767.

    [22] 高燕燕.重塑黃土的滲透特性研究[D].西安:長安大學(xué),2016.

    [23] 高文靜.細粒土滲透特性的影響因素試驗研究[D].西安:西安理工大學(xué),2018.

    (編輯:趙秋云)

    Quantitative analysis on NMR T2 spectrum in remolded loess and application

    CAO Xuemei1,ZHAN Yongkang2,YUAN Lei3,SHU Zhongpan4,ZHANG Jun4

    (1.Sichuan Vocational and Technical College of Communications,Chengdu 611130,China; 2.Sichuan Aba Traffic Construction Investment Co.,Ltd.,Aba 624000,China; 3.Transportation Bureau of Aba Tibetan and Qiang Autonomous Prefecture,Aba 624000 ,China; 4.Sichuan Institute of Geological Engineering Investigation Group Co.,Ltd.,Chengdu 610072 ,China)

    Abstract:

    Nuclear Magnetic Resonance (NMR) is widely used in the analysis of pore space characteristics of loess due to its convenience,low cost,and no damage.However,the previous analysis and application of remolded loess T2 spectrum were limited to obtaining the distribution of pore size,and other spatial information contained in the T2 spectrum was not further extracted,which made the application results of the NMR T2 spectrum too single,resulting in a certain degree of data waste.Based on this,we used mathematical analysis method to quantitatively analyze the T2 spectrum of remolded loess with different dry densities,and applied the mathematical parameters obtained by quantitative analysis to the SDR model,so as to construct a permeability calculation model for remodeled loess.The main conclusions were as follows:(1) There are various differences in the statistical analysis parameters (mean,geometric mean,standard deviation and kurtosis) of the NMR T2 spectrum of loess remolded with different dry densities,which has physical meanings:different dry densities of remolded loess will affect the size,distribution and sorting characteristics of pores in the soil.Compared with the traditional method to obtain the pore volume of soil,the mathematical parameters obtained by numerical statistics method can describe the empty pores more comprehensively and accurately.(2) The permeability of each remolded loess was calculated by constructing the SDR model of the remolded loess,the minimum permeability was about 0.3 mD,and the calculated results were close to the previous experimental results.(3) Dry density affected micro-spatial characteristics of soil,and causing the change of soil permeability characteristics.(4) From the correlation analysis,the mean and geometric mean of T2 spectrum had great influence on soil permeability.To a certain extent,the research results can provide new ideas for the application of NMR in the loess study field,new methods for predicting the permeability of loess,and new thinking for understanding the mechanism of soil permeability.

    Key words:

    remolded loess;pore space;Nuclear Magnetic Resonance (NMR);microscopic characteristics;quantitative parameters

    猜你喜歡
    重塑黃土滲透率
    重塑未來
    中國慈善家(2021年5期)2021-11-19 18:38:58
    各路創(chuàng)新人才涌向“黃土高坡”
    當代陜西(2020年15期)2021-01-07 09:23:14
    黃土成金
    河北畫報(2020年21期)2020-12-14 03:16:04
    只要有信心 黃土變成金
    河北畫報(2020年21期)2020-12-14 03:16:04
    《劉文西:繪不盡是黃土情》
    自動化正悄然無聲地重塑服務(wù)業(yè)
    英語文摘(2019年6期)2019-09-18 01:49:08
    中煤階煤層氣井排采階段劃分及滲透率變化
    不同滲透率巖芯孔徑分布與可動流體研究
    李滄:再造與重塑
    商周刊(2018年11期)2018-06-13 03:41:54
    SAGD井微壓裂儲層滲透率變化規(guī)律研究
    91国产中文字幕| av免费在线看不卡| 免费av不卡在线播放| 国国产精品蜜臀av免费| www.色视频.com| 亚洲色图 男人天堂 中文字幕 | 亚洲精品aⅴ在线观看| 亚洲精华国产精华液的使用体验| 欧美成人精品欧美一级黄| 国产黄色视频一区二区在线观看| 亚洲精品久久午夜乱码| 中国三级夫妇交换| 亚洲激情五月婷婷啪啪| 精品少妇黑人巨大在线播放| 久久久a久久爽久久v久久| 性色av一级| 色婷婷久久久亚洲欧美| 国产成人免费无遮挡视频| 男女边吃奶边做爰视频| 日韩不卡一区二区三区视频在线| 中文天堂在线官网| 美女中出高潮动态图| 妹子高潮喷水视频| 在线免费观看不下载黄p国产| 美国免费a级毛片| 国产男人的电影天堂91| 香蕉丝袜av| 亚洲成av片中文字幕在线观看 | 最近最新中文字幕免费大全7| 欧美精品一区二区大全| 国产成人精品福利久久| 国产淫语在线视频| 高清在线视频一区二区三区| 超碰97精品在线观看| 国产精品一区www在线观看| 多毛熟女@视频| 国产亚洲最大av| 欧美最新免费一区二区三区| 久久久a久久爽久久v久久| 国产一级毛片在线| 最近的中文字幕免费完整| 夜夜骑夜夜射夜夜干| 久久久久久伊人网av| 亚洲欧美一区二区三区黑人 | 久热久热在线精品观看| 一级黄片播放器| 亚洲高清免费不卡视频| 亚洲高清免费不卡视频| 欧美日韩视频高清一区二区三区二| 亚洲内射少妇av| 少妇人妻 视频| 最新的欧美精品一区二区| 婷婷色麻豆天堂久久| 久久久久久久久久久久大奶| 久久久久久久久久人人人人人人| 亚洲精品乱码久久久久久按摩| 日本色播在线视频| 久久久久视频综合| xxxhd国产人妻xxx| 国产爽快片一区二区三区| 国产一区有黄有色的免费视频| 日韩成人伦理影院| 日韩人妻精品一区2区三区| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧洲日产国产| av电影中文网址| 韩国av在线不卡| 亚洲少妇的诱惑av| 国产成人一区二区在线| 国产永久视频网站| 亚洲伊人久久精品综合| 视频中文字幕在线观看| 亚洲精品国产av蜜桃| 美女中出高潮动态图| 精品少妇黑人巨大在线播放| 91国产中文字幕| 成人国产av品久久久| 欧美丝袜亚洲另类| 国产免费视频播放在线视频| 亚洲 欧美一区二区三区| 九色成人免费人妻av| 大码成人一级视频| xxxhd国产人妻xxx| 日韩一区二区三区影片| 国产精品麻豆人妻色哟哟久久| 香蕉丝袜av| 成人亚洲精品一区在线观看| 18禁国产床啪视频网站| 涩涩av久久男人的天堂| 性色avwww在线观看| 精品第一国产精品| 国产成人免费观看mmmm| 亚洲av电影在线观看一区二区三区| 99热6这里只有精品| 男女午夜视频在线观看 | 男女高潮啪啪啪动态图| 2021少妇久久久久久久久久久| 国产精品久久久久久精品古装| 在线亚洲精品国产二区图片欧美| 国产亚洲一区二区精品| 亚洲情色 制服丝袜| 伦理电影免费视频| 男女国产视频网站| 2018国产大陆天天弄谢| 一级a做视频免费观看| 国产综合精华液| 国产av国产精品国产| 在线观看免费高清a一片| 精品一区二区三区视频在线| 极品少妇高潮喷水抽搐| 国产欧美亚洲国产| 国产极品粉嫩免费观看在线| 中国三级夫妇交换| av又黄又爽大尺度在线免费看| 丝袜人妻中文字幕| 国产亚洲精品第一综合不卡 | 亚洲国产精品成人久久小说| 嫩草影院入口| 亚洲精品av麻豆狂野| 成人影院久久| 久久久国产一区二区| 亚洲精品乱久久久久久| 欧美精品一区二区免费开放| 亚洲av中文av极速乱| 国产成人精品无人区| 热re99久久精品国产66热6| 亚洲国产欧美在线一区| 日韩大片免费观看网站| 丰满乱子伦码专区| 亚洲熟女精品中文字幕| 一二三四中文在线观看免费高清| a级毛色黄片| 亚洲精品中文字幕在线视频| 女人被躁到高潮嗷嗷叫费观| www.熟女人妻精品国产 | 国产亚洲欧美精品永久| av免费在线看不卡| 成人毛片60女人毛片免费| 肉色欧美久久久久久久蜜桃| 七月丁香在线播放| 国产欧美另类精品又又久久亚洲欧美| 国产一区二区激情短视频 | 精品熟女少妇av免费看| 日本黄大片高清| 男的添女的下面高潮视频| 这个男人来自地球电影免费观看 | 国产有黄有色有爽视频| 免费女性裸体啪啪无遮挡网站| 久久人妻熟女aⅴ| 99精国产麻豆久久婷婷| 在线看a的网站| 插逼视频在线观看| 最近最新中文字幕大全免费视频 | 欧美人与性动交α欧美精品济南到 | 中文精品一卡2卡3卡4更新| 男女高潮啪啪啪动态图| 最近手机中文字幕大全| 国产欧美日韩一区二区三区在线| 亚洲精品美女久久久久99蜜臀 | 日韩大片免费观看网站| 国产av码专区亚洲av| 考比视频在线观看| 午夜免费男女啪啪视频观看| 大香蕉久久网| 免费观看av网站的网址| 精品午夜福利在线看| 国产欧美另类精品又又久久亚洲欧美| 51国产日韩欧美| 久久国产精品男人的天堂亚洲 | 国产永久视频网站| 欧美日韩成人在线一区二区| 哪个播放器可以免费观看大片| 日本wwww免费看| 精品国产露脸久久av麻豆| 精品酒店卫生间| 狠狠婷婷综合久久久久久88av| 成年动漫av网址| 高清视频免费观看一区二区| 九九在线视频观看精品| 亚洲一码二码三码区别大吗| 视频区图区小说| 夫妻午夜视频| 一区二区三区四区激情视频| 国内精品宾馆在线| 99久久精品国产国产毛片| 亚洲精品av麻豆狂野| 黑人猛操日本美女一级片| 国产亚洲最大av| 日本av手机在线免费观看| av在线播放精品| 一二三四中文在线观看免费高清| 咕卡用的链子| av女优亚洲男人天堂| 国产 精品1| 涩涩av久久男人的天堂| 天堂中文最新版在线下载| 香蕉国产在线看| 久久精品久久久久久噜噜老黄| 一级片免费观看大全| 免费人妻精品一区二区三区视频| 免费观看a级毛片全部| 亚洲国产精品一区二区三区在线| 丝袜人妻中文字幕| av福利片在线| 极品人妻少妇av视频| 在线观看三级黄色| 男女边吃奶边做爰视频| 亚洲欧洲精品一区二区精品久久久 | 波野结衣二区三区在线| 亚洲国产精品专区欧美| 一本久久精品| 日韩 亚洲 欧美在线| 久久久精品区二区三区| 不卡视频在线观看欧美| 2021少妇久久久久久久久久久| 成人无遮挡网站| 黄色视频在线播放观看不卡| 婷婷色av中文字幕| 十八禁高潮呻吟视频| 久久精品人人爽人人爽视色| 一区二区av电影网| 国产成人av激情在线播放| 国产在线一区二区三区精| 精品人妻在线不人妻| 捣出白浆h1v1| 国产欧美亚洲国产| 寂寞人妻少妇视频99o| 日韩一本色道免费dvd| 这个男人来自地球电影免费观看 | av片东京热男人的天堂| 午夜91福利影院| 男人爽女人下面视频在线观看| 亚洲欧美中文字幕日韩二区| 麻豆乱淫一区二区| 亚洲精品一二三| 国产男靠女视频免费网站| 香蕉丝袜av| 欧美日韩一级在线毛片| 视频在线观看一区二区三区| 亚洲人成电影免费在线| 少妇裸体淫交视频免费看高清 | 婷婷丁香在线五月| 国产高清videossex| 窝窝影院91人妻| 男人舔女人的私密视频| 欧美黑人欧美精品刺激| 欧美日韩乱码在线| 超色免费av| av天堂久久9| 亚洲中文字幕日韩| 国产精品.久久久| 国产亚洲欧美精品永久| 国产视频一区二区在线看| 国产主播在线观看一区二区| 香蕉国产在线看| 国产成人免费观看mmmm| 欧美激情久久久久久爽电影 | 啦啦啦 在线观看视频| 亚洲av日韩精品久久久久久密| 国产又色又爽无遮挡免费看| 色婷婷久久久亚洲欧美| 欧美不卡视频在线免费观看 | 757午夜福利合集在线观看| 中文亚洲av片在线观看爽 | 男女午夜视频在线观看| 久久中文看片网| 99re在线观看精品视频| av视频免费观看在线观看| 亚洲精品国产精品久久久不卡| 国产伦人伦偷精品视频| 国产精品美女特级片免费视频播放器 | 国产av精品麻豆| 色综合欧美亚洲国产小说| 九色亚洲精品在线播放| 日本精品一区二区三区蜜桃| 国产成人一区二区三区免费视频网站| 91在线观看av| 啦啦啦免费观看视频1| 精品国产一区二区三区四区第35| 大香蕉久久成人网| 中文字幕高清在线视频| 婷婷精品国产亚洲av在线 | 制服诱惑二区| 自拍欧美九色日韩亚洲蝌蚪91| 水蜜桃什么品种好| 国产91精品成人一区二区三区| 巨乳人妻的诱惑在线观看| 免费在线观看完整版高清| 成人三级做爰电影| 999精品在线视频| 欧美日本中文国产一区发布| 两人在一起打扑克的视频| 真人做人爱边吃奶动态| 色老头精品视频在线观看| 亚洲专区中文字幕在线| 午夜两性在线视频| 老汉色∧v一级毛片| 村上凉子中文字幕在线| 午夜福利欧美成人| 搡老岳熟女国产| 国产成人av教育| 久久这里只有精品19| 国产精品永久免费网站| 日日摸夜夜添夜夜添小说| 97人妻天天添夜夜摸| 极品少妇高潮喷水抽搐| 亚洲av成人av| 99热国产这里只有精品6| 搡老熟女国产l中国老女人| 欧美av亚洲av综合av国产av| 亚洲一码二码三码区别大吗| 国产男女内射视频| 黑人操中国人逼视频| 另类亚洲欧美激情| 精品熟女少妇八av免费久了| 亚洲va日本ⅴa欧美va伊人久久| 欧美日韩黄片免| 国产xxxxx性猛交| 国产欧美日韩精品亚洲av| 99riav亚洲国产免费| 成人三级做爰电影| 50天的宝宝边吃奶边哭怎么回事| 在线十欧美十亚洲十日本专区| 18禁观看日本| 精品人妻1区二区| 国产欧美日韩精品亚洲av| 90打野战视频偷拍视频| 91麻豆av在线| 亚洲熟女毛片儿| 亚洲精品国产精品久久久不卡| 欧美av亚洲av综合av国产av| 午夜福利在线观看吧| 99香蕉大伊视频| 操美女的视频在线观看| 国产激情欧美一区二区| 又紧又爽又黄一区二区| 妹子高潮喷水视频| 精品国产乱码久久久久久男人| 亚洲黑人精品在线| 91精品国产国语对白视频| 国产精品1区2区在线观看. | 国产1区2区3区精品| 建设人人有责人人尽责人人享有的| 午夜免费观看网址| 久久人妻熟女aⅴ| 无人区码免费观看不卡| 欧美激情 高清一区二区三区| 国产亚洲欧美在线一区二区| 国产精品自产拍在线观看55亚洲 | 少妇裸体淫交视频免费看高清 | 成人18禁高潮啪啪吃奶动态图| a级毛片黄视频| 99re6热这里在线精品视频| 国产精品秋霞免费鲁丝片| 精品国产国语对白av| 黄色女人牲交| 男女之事视频高清在线观看| 亚洲人成伊人成综合网2020| 一级毛片高清免费大全| 一a级毛片在线观看| 欧美激情高清一区二区三区| 一进一出抽搐动态| 免费在线观看日本一区| 欧美日韩亚洲高清精品| 欧美老熟妇乱子伦牲交| 在线看a的网站| 校园春色视频在线观看| 丝袜在线中文字幕| 亚洲精品美女久久久久99蜜臀| 在线观看免费日韩欧美大片| 欧美激情久久久久久爽电影 | 啦啦啦视频在线资源免费观看| 国产蜜桃级精品一区二区三区 | 亚洲午夜理论影院| 黑人巨大精品欧美一区二区蜜桃| 中亚洲国语对白在线视频| 亚洲情色 制服丝袜| av线在线观看网站| 欧美黄色淫秽网站| 成年版毛片免费区| 免费看十八禁软件| 亚洲精品一卡2卡三卡4卡5卡| 91麻豆精品激情在线观看国产 | 90打野战视频偷拍视频| 午夜精品久久久久久毛片777| 啦啦啦视频在线资源免费观看| 亚洲国产精品合色在线| 国产真人三级小视频在线观看| 国产又爽黄色视频| 久久午夜亚洲精品久久| 91国产中文字幕| 亚洲精品久久午夜乱码| 精品国产超薄肉色丝袜足j| 久久久久久人人人人人| 国产欧美日韩一区二区三区在线| 最近最新免费中文字幕在线| 中文亚洲av片在线观看爽 | 久久精品国产亚洲av香蕉五月 | 亚洲欧洲精品一区二区精品久久久| 精品少妇一区二区三区视频日本电影| 人妻 亚洲 视频| 麻豆国产av国片精品| 免费日韩欧美在线观看| 国产真人三级小视频在线观看| 国产精品 国内视频| 精品第一国产精品| 成人影院久久| 一本综合久久免费| 精品视频人人做人人爽| 久久天躁狠狠躁夜夜2o2o| 啦啦啦在线免费观看视频4| 嫩草影视91久久| 亚洲成人免费av在线播放| 亚洲av日韩在线播放| 久99久视频精品免费| 99久久综合精品五月天人人| 欧美国产精品一级二级三级| 欧美最黄视频在线播放免费 | 欧美日韩黄片免| 久久人妻福利社区极品人妻图片| 欧美 亚洲 国产 日韩一| 一区福利在线观看| 久久精品91无色码中文字幕| 国产97色在线日韩免费| a级毛片黄视频| 亚洲欧美激情在线| 国产亚洲欧美在线一区二区| 久久久久久久久久久久大奶| av网站在线播放免费| 亚洲中文字幕日韩| 99国产精品一区二区蜜桃av | 女同久久另类99精品国产91| 嫁个100分男人电影在线观看| 国产精品一区二区在线不卡| 久久精品熟女亚洲av麻豆精品| 国产黄色免费在线视频| 精品少妇一区二区三区视频日本电影| 久久国产精品大桥未久av| 精品第一国产精品| 999久久久精品免费观看国产| 国产精品一区二区在线观看99| 久久久久久久久免费视频了| 少妇 在线观看| av欧美777| 亚洲五月婷婷丁香| 欧美人与性动交α欧美软件| 90打野战视频偷拍视频| 色精品久久人妻99蜜桃| 亚洲一区二区三区不卡视频| 大香蕉久久网| 亚洲精品国产色婷婷电影| 国产又爽黄色视频| 日本wwww免费看| 丁香欧美五月| 老司机午夜十八禁免费视频| 国产精品久久久久久精品古装| 中亚洲国语对白在线视频| 一区二区三区国产精品乱码| 久久国产精品男人的天堂亚洲| 午夜福利在线观看吧| 亚洲情色 制服丝袜| av一本久久久久| av片东京热男人的天堂| 国产亚洲av高清不卡| 极品少妇高潮喷水抽搐| 国产不卡一卡二| av网站免费在线观看视频| 亚洲人成伊人成综合网2020| 国产成人欧美| 国产成人影院久久av| 久久婷婷成人综合色麻豆| 久久久久国内视频| 亚洲精品国产区一区二| 在线观看一区二区三区激情| 夜夜夜夜夜久久久久| 国产精品国产高清国产av | 免费高清在线观看日韩| 在线av久久热| 亚洲av成人不卡在线观看播放网| 久9热在线精品视频| 午夜两性在线视频| 久久天躁狠狠躁夜夜2o2o| 正在播放国产对白刺激| 黑人巨大精品欧美一区二区蜜桃| a在线观看视频网站| 欧美精品亚洲一区二区| 人人妻人人爽人人添夜夜欢视频| 91大片在线观看| 亚洲欧美色中文字幕在线| x7x7x7水蜜桃| 窝窝影院91人妻| 免费在线观看完整版高清| av福利片在线| 久久婷婷成人综合色麻豆| 人人妻人人澡人人爽人人夜夜| a在线观看视频网站| 美女福利国产在线| 国产91精品成人一区二区三区| 国产极品粉嫩免费观看在线| 一区二区三区激情视频| 久久午夜亚洲精品久久| 国产免费男女视频| 一区二区三区精品91| 欧美最黄视频在线播放免费 | 亚洲七黄色美女视频| 真人做人爱边吃奶动态| 亚洲性夜色夜夜综合| 在线国产一区二区在线| 一二三四社区在线视频社区8| 9色porny在线观看| 人妻一区二区av| 中出人妻视频一区二区| 日韩欧美三级三区| e午夜精品久久久久久久| 成年人午夜在线观看视频| 亚洲精品一二三| 久久ye,这里只有精品| 妹子高潮喷水视频| 日韩一卡2卡3卡4卡2021年| 亚洲av成人不卡在线观看播放网| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲伊人色综图| 亚洲欧美激情综合另类| 国产av一区二区精品久久| 大码成人一级视频| 日本欧美视频一区| 淫妇啪啪啪对白视频| 亚洲第一av免费看| 99久久综合精品五月天人人| 97人妻天天添夜夜摸| 欧美日韩中文字幕国产精品一区二区三区 | 久久人人爽av亚洲精品天堂| 老汉色av国产亚洲站长工具| 日日夜夜操网爽| 人妻 亚洲 视频| 久久精品国产99精品国产亚洲性色 | 国产精品1区2区在线观看. | 久久青草综合色| 亚洲欧美精品综合一区二区三区| 波多野结衣av一区二区av| a级毛片在线看网站| 亚洲情色 制服丝袜| 黑人巨大精品欧美一区二区蜜桃| а√天堂www在线а√下载 | 国产精品综合久久久久久久免费 | 中文字幕av电影在线播放| 另类亚洲欧美激情| 亚洲成国产人片在线观看| 一边摸一边做爽爽视频免费| 亚洲五月色婷婷综合| 久久中文字幕人妻熟女| 欧美日韩福利视频一区二区| 欧美日韩一级在线毛片| 国产在视频线精品| 精品国产美女av久久久久小说| 母亲3免费完整高清在线观看| 最新美女视频免费是黄的| videos熟女内射| 五月开心婷婷网| 午夜成年电影在线免费观看| 老熟妇仑乱视频hdxx| 亚洲精品中文字幕一二三四区| 人妻丰满熟妇av一区二区三区 | 国产黄色免费在线视频| 成年版毛片免费区| 久久精品成人免费网站| av中文乱码字幕在线| 99久久精品国产亚洲精品| 在线观看www视频免费| 国产一区二区三区在线臀色熟女 | 日本vs欧美在线观看视频| 身体一侧抽搐| 69av精品久久久久久| 国产高清视频在线播放一区| 精品亚洲成a人片在线观看| 久久国产亚洲av麻豆专区| 久久精品国产a三级三级三级| 在线av久久热| 亚洲男人天堂网一区| 两个人看的免费小视频| 欧美日韩乱码在线| 91九色精品人成在线观看| 视频在线观看一区二区三区| av不卡在线播放| 久99久视频精品免费| 日韩欧美一区视频在线观看| 欧美久久黑人一区二区| 高清黄色对白视频在线免费看| 757午夜福利合集在线观看| 两性夫妻黄色片| 变态另类成人亚洲欧美熟女 | 久久久久久免费高清国产稀缺| 波多野结衣一区麻豆| 欧美黑人精品巨大| 好看av亚洲va欧美ⅴa在| 欧美日韩乱码在线| 中文字幕制服av| 午夜成年电影在线免费观看| 欧美黑人欧美精品刺激| 欧美日韩av久久| 这个男人来自地球电影免费观看| 亚洲国产精品sss在线观看 | 夫妻午夜视频| 国产不卡av网站在线观看| 国产免费男女视频| 亚洲精品粉嫩美女一区| 亚洲精品一卡2卡三卡4卡5卡| 亚洲免费av在线视频| 国产淫语在线视频| 亚洲欧美一区二区三区久久| 久久精品国产亚洲av高清一级| 在线观看www视频免费| 亚洲国产中文字幕在线视频| 热99久久久久精品小说推荐| 亚洲va日本ⅴa欧美va伊人久久|