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

    單軸壓縮鹽巖聲發(fā)射特征及損傷演化探討

    2016-05-30 02:59:02周志威劉建鋒徐楊夢迪
    長江科學(xué)院院報(bào) 2016年5期
    關(guān)鍵詞:聲發(fā)射鹽巖

    周志威,劉建鋒,鄒 航,卓 越,徐楊夢迪

    (四川大學(xué) 水力學(xué)與山區(qū)河流開發(fā)保護(hù)國家重點(diǎn)實(shí)驗(yàn)室,成都 610065)

    ?

    單軸壓縮鹽巖聲發(fā)射特征及損傷演化探討

    周志威,劉建鋒,鄒航,卓越,徐楊夢迪

    (四川大學(xué) 水力學(xué)與山區(qū)河流開發(fā)保護(hù)國家重點(diǎn)實(shí)驗(yàn)室,成都610065)

    摘要:為了探究鹽巖壓縮變形破壞過程中,聲發(fā)射現(xiàn)象與內(nèi)部損傷的相關(guān)關(guān)系,利用四川大學(xué)MTS815巖石力學(xué)試驗(yàn)系統(tǒng),對(duì)云南安寧鹽巖進(jìn)行單軸壓縮試驗(yàn),分析了鹽巖聲發(fā)射特征。依據(jù)損傷理論,得到了基于聲發(fā)射參數(shù)的鹽巖的損傷演化模型,并引入修正系數(shù)對(duì)其進(jìn)行理論修正,分析損傷終值對(duì)損傷本構(gòu)模型的影響。研究結(jié)果表明:鹽巖的聲發(fā)射(acoustic emission,簡稱AE)振鈴計(jì)數(shù)率和能量率在彈性極限附近達(dá)到最大值,之后呈現(xiàn)高頻低幅值的聲發(fā)射現(xiàn)象;基于聲發(fā)射的鹽巖損傷變量值在壓縮變形前期增加較快;基于AE累計(jì)振鈴計(jì)數(shù)得到的損傷值<基于AE累計(jì)能量計(jì)算得到的鹽巖損傷值;修正后的聲發(fā)射損傷理論模型曲線與試驗(yàn)曲線全過程吻合較好,能夠較好地反映鹽巖內(nèi)部損傷演化過程;振鈴計(jì)數(shù)參數(shù)模型比能量參數(shù)模型能更好地模擬鹽巖的損傷演化特征。

    關(guān)鍵詞:鹽巖;聲發(fā)射;損傷模型;損傷變量;修正系數(shù)

    1研究背景

    鹽巖是國際上公認(rèn)的石油、天然氣等能源地下儲(chǔ)備的理想介質(zhì),也是中國能源地下儲(chǔ)備的主要建庫場所。鹽巖損傷會(huì)直接危及儲(chǔ)庫的安全運(yùn)營,因此鹽巖損傷研究是鹽穴能源儲(chǔ)庫研究的重要內(nèi)容之一。謝和平等[1-4]基于循環(huán)荷載試驗(yàn)對(duì)鹽巖損傷演化特征進(jìn)行了探討,揭示了鹽巖應(yīng)力狀態(tài)對(duì)損傷發(fā)展的影響。Schulze等[5-6]根據(jù)鹽巖變形過程中滲透性的變化規(guī)律,對(duì)鹽巖損傷發(fā)展特征進(jìn)行了研究。

    聲發(fā)射(acoustic emission,簡稱AE)作為一種無損監(jiān)測技術(shù),被廣泛應(yīng)用于材料變形破壞全過程的損傷研究。國內(nèi)外利用聲發(fā)射對(duì)巖石和混凝土等材料的損傷進(jìn)行了大量的研究。Ohtsu[7-8]運(yùn)用速率過程理論將聲發(fā)射用于混凝土的損傷評(píng)估當(dāng)中,認(rèn)為聲發(fā)射活性與材料所含裂紋的多少成正相關(guān)。Heiple等[9-10]通過大量實(shí)驗(yàn)分析,認(rèn)為聲發(fā)射振鈴計(jì)數(shù)更能反映材料內(nèi)部損傷演化的過程。Tang和Xu[11]采用連續(xù)損傷力學(xué)方法,從細(xì)觀角度出發(fā),得到了單軸壓縮條件下巖石材料基于聲發(fā)射的損傷模型。國內(nèi)一些學(xué)者利用聲發(fā)射對(duì)煤等巖石的損傷進(jìn)行了研究,揭示了研究巖石的聲發(fā)射特征及基于聲發(fā)射的應(yīng)力-應(yīng)變損傷演化全過程[12-14]。

    損傷理論將橫截面破壞面積與無損狀態(tài)下橫截面積之比定義為損傷變量D。在對(duì)壓縮應(yīng)力狀態(tài)下的巖石損傷演化研究中,為獲得損傷演化模型中的材料參數(shù),通常假定加載前巖石無損傷,即損傷變量初值D=0,當(dāng)巖石達(dá)到抗壓強(qiáng)度時(shí)完全損傷,即損傷變量D=1[15]。事實(shí)上,巖石在壓縮應(yīng)力狀態(tài)下,當(dāng)壓縮應(yīng)力達(dá)到抗壓強(qiáng)度時(shí),橫截面的損傷呈單裂紋或多裂紋破壞特征,仍具有有效承載面積,橫截面并非完全失效,因此,根據(jù)損傷變量定義得到的破壞后的損傷變量D(即損傷終值)應(yīng)<1。假定D=1必然不能反映巖石的真實(shí)損傷特征,為盡可能真實(shí)地反映巖石的損傷演化特征,文獻(xiàn)[12-14]通過聲發(fā)射試驗(yàn),提出了單軸壓縮應(yīng)力下殘余應(yīng)力與峰值應(yīng)力之比的損傷變量修正方法,并對(duì)巖石破壞過程中的損傷本構(gòu)模型進(jìn)行了研究。

    大量研究表明,鹽巖變形破壞過程中的聲發(fā)射特征更為顯著[16-18],因此,基于聲發(fā)射對(duì)鹽巖的損傷研究具有較強(qiáng)的現(xiàn)實(shí)意義,而目前針對(duì)聲發(fā)射特征對(duì)鹽巖損傷演化的研究還相對(duì)較少。

    本文通過鹽巖單軸壓縮聲發(fā)射試驗(yàn),利用損傷修正系數(shù)η對(duì)損傷本構(gòu)模型進(jìn)行修正,分析損傷終值對(duì)損傷本構(gòu)模型的影響,探討利用聲發(fā)射作為描述損傷演化的參量進(jìn)行損傷研究時(shí),聲發(fā)射能量和振鈴計(jì)數(shù)哪一個(gè)更能反映鹽巖的損傷演化特征。

    2基于聲發(fā)射特征參數(shù)的損傷模型

    20世紀(jì)60年代,Kachanov[19]首次提出了連續(xù)度φ這一概念,并用連續(xù)度來描述材料內(nèi)部的逐漸損傷。在連續(xù)度φ的基礎(chǔ)上,對(duì)于單軸應(yīng)力下的各向同性損傷,Kachanov隨后又提出損傷因子(即損傷變量)D,即

    (1)

    式中:A′為損傷后截面有效面積;A為無損狀態(tài)時(shí)的橫截面面積;Ad為破壞截面面積。

    對(duì)于上述表述材料損傷的公式,由于材料的有效承載面積需通過復(fù)雜的統(tǒng)計(jì)分析才能確定,故為了簡化考慮損傷因素后的結(jié)構(gòu)計(jì)算物性方程,1972年,Lemaitre[20]提出了應(yīng)變等效性假設(shè),通過簡單的等效變換可在原有無損材料的本構(gòu)方程基礎(chǔ)上建立損傷材料的本構(gòu)方程,從而在一定程度上促進(jìn)了損傷力學(xué)的快速發(fā)展。該假設(shè)認(rèn)為,對(duì)于彈性材料,在真實(shí)應(yīng)力σ作用下,受損狀態(tài)的應(yīng)變等效于在有效應(yīng)力σ′作用下虛構(gòu)的無損狀態(tài)的應(yīng)變,那么則有

    (2)

    式中:E為變形模量;ε為應(yīng)變。

    假設(shè)斷面微缺陷面積達(dá)到A時(shí)的各累積聲發(fā)射參數(shù)值為N;斷面面積A完全破壞時(shí),即巖石發(fā)生破壞時(shí)的聲發(fā)射參數(shù)累積值為Nm,那么當(dāng)截面損傷破壞面積達(dá)到Ad時(shí),這一時(shí)刻聲發(fā)射參數(shù)累計(jì)值即為

    (3)

    聯(lián)立式(1)和式(3),可以得到損傷變量D與聲發(fā)射參數(shù)值之間的歸一化方程,即

    (4)

    將式(4)代入式(2),即可得出應(yīng)力-應(yīng)變與聲發(fā)射參數(shù)之間的相關(guān)關(guān)系方程,即

    (5)

    通常假定巖石破壞后,試件承載面積全部破壞,事實(shí)上并非如此。

    圖1 鹽巖試驗(yàn)前、后的CT掃描結(jié)果Fig.1 CT scanning images of rock saltbefore and after test

    圖1給出了鹽巖在單軸應(yīng)力下破壞后CT掃描得到的試件內(nèi)部裂紋面分布特征,雖然試件破壞后,產(chǎn)生了大量的裂紋,然而從統(tǒng)計(jì)分析學(xué)來講,試驗(yàn)結(jié)束后試件承載面并未發(fā)生完全破壞,有效面積A′確實(shí)是存在的。該圖表明破壞后的橫截面并非完全失效,故此時(shí)的損傷變量D(損傷終值)應(yīng)<1,同時(shí)也進(jìn)一步說明了單軸壓縮應(yīng)力下,巖石峰后仍具有承載能力的原因。因此,針對(duì)損傷終值<1這一實(shí)際情況,引入損傷修正系數(shù)η(η≤1,本文假定初始時(shí)鹽巖處于無損狀態(tài))對(duì)鹽巖的損傷作出修正,那么可得

    (6)

    (7)

    式(7)為修正后基于聲發(fā)射參數(shù)值的損傷本構(gòu)方程。

    圖2 標(biāo)準(zhǔn)試驗(yàn)試件Fig.2 Standard testspecimens

    3試件制備及試驗(yàn)設(shè)備

    3.1試驗(yàn)試件及制備

    本文鹽巖試件取自云

    南安寧擬建天然儲(chǔ)氣庫鹽礦,平均埋深約500 m。試驗(yàn)參照相關(guān)規(guī)定,采用車床干車法將巖樣制成尺寸為φ90 mm×180 mm的標(biāo)準(zhǔn)圓柱形試件,編號(hào)分別為1#,2#(如圖2)。

    3.2試驗(yàn)設(shè)備

    加載設(shè)備為四川大學(xué)MTS815 Flex Test GT巖石力學(xué)測試系統(tǒng)和PCI-Ⅱ聲發(fā)射(AE)三維定位系統(tǒng)(圖3)。

    圖3 MTS815巖石力學(xué)測試系統(tǒng)Fig.3 Rock mechanics testing system MTS815

    試驗(yàn)軸向荷載傳感器1 000 kN,軸向引伸計(jì)量程-4~4 mm,環(huán)向引伸計(jì)量程為-2.5~8 mm,各測量傳感器精度均達(dá)到標(biāo)定量程點(diǎn)的0.5%。聲發(fā)射監(jiān)測采用8個(gè)傳感器對(duì)稱分布于試件柱面上下兩端,每端4個(gè),傳感器頻率為100 kHz。為保證傳感器與鹽巖試件的耦合效果,在二者接觸部位涂凡士林。

    4試驗(yàn)結(jié)果及分析

    4.1聲發(fā)射特征分析

    鹽巖的應(yīng)力-應(yīng)變?nèi)^程示例曲線和鹽巖試件破壞前、后照片的對(duì)比見圖4和圖5。

    圖4 應(yīng)力-應(yīng)變?nèi)^程曲線Fig.4 Complete stress-strain curves

    圖5 鹽巖破壞后示例Fig.5 Photos of salt rock specimens after failure

    圖6 AE各特征參數(shù)曲線Fig.6 AE characteristic curves

    從圖4可看出,鹽巖的應(yīng)力-應(yīng)變?nèi)^程大致分為3個(gè)階段:①彈性變形階段(0~10 MPa),應(yīng)力增加較快而應(yīng)變?cè)鲩L較慢;②塑性變形階段,隨著應(yīng)力增加,應(yīng)變開始快速增大,加載應(yīng)力上升速度減緩,此過程持續(xù)時(shí)間較長,應(yīng)變值變化較大,塑性應(yīng)變量占整個(gè)應(yīng)變值的80%~90%;③峰后破壞階段,應(yīng)力到達(dá)峰值應(yīng)力后,試件發(fā)生破壞,應(yīng)力值下降一定幅度之后即終止試驗(yàn)。1#試件的峰值應(yīng)力(26.3 MPa)比2#試件的峰值應(yīng)力(23.73 MPa)大10%左右,而峰值應(yīng)力對(duì)應(yīng)的應(yīng)變值大概為2#試件的75%。對(duì)試驗(yàn)后的鹽巖試件進(jìn)行化學(xué)成分分析,得到1#試件的雜質(zhì)含量(約20%)高于2#試件的雜質(zhì)含量(約10%),通常情況下,雜質(zhì)含量越高,鹽巖的強(qiáng)度越大,塑性變形能力相應(yīng)越低[21]。

    由圖5可知,在單軸壓縮條件下,鹽巖試件主要表現(xiàn)為破裂面平行于軸向力方向的貫通張性破壞特征。

    圖6為鹽巖的AE特征曲線。由圖6可看出,對(duì)于鹽巖這種塑性較強(qiáng)的巖石,其聲發(fā)射總體表現(xiàn)為頻度高、強(qiáng)度低。整個(gè)加載試驗(yàn)過程,聲發(fā)射振鈴計(jì)數(shù)率均處于一個(gè)比較高的計(jì)數(shù)水平,最大值都出現(xiàn)在600~700次/s之間,整個(gè)過程并沒有太大的起伏變化,1#試件平均維持在300~400次/s,2#試件整個(gè)過程則維持在400~500次/s,因而累計(jì)振鈴計(jì)數(shù)接近于線性增長。測試鹽巖的聲發(fā)射能量率變化特征基本一致,都是先增大,后逐步減小并最終趨于低能穩(wěn)定,最大值都在28左右,因而累計(jì)能量值在加載前期(彈性階段)快速增長,占整個(gè)試驗(yàn)總值的40%以上,后期則逐漸趨于低速率緩慢增加。

    對(duì)圖6(a)和圖6(b)進(jìn)行細(xì)致分析,發(fā)現(xiàn)鹽巖的聲發(fā)射現(xiàn)象在各個(gè)階段又不大相同。加載初期AE活動(dòng)較為活躍,聲發(fā)射振鈴計(jì)數(shù)和能量隨著加載應(yīng)力增大而增多,試件處在彈性壓密階段,原生缺陷出現(xiàn)閉合、鹽晶晶粒間產(chǎn)生摩擦和滑移而產(chǎn)生聲發(fā)射信號(hào),同時(shí)鹽巖的各向異性和非均勻性也會(huì)引起聲發(fā)射現(xiàn)象。當(dāng)應(yīng)力加載到彈性極限承載力左右時(shí),能量率達(dá)到最大值,此時(shí)鹽巖內(nèi)部初始缺陷基本閉合,開始產(chǎn)生新裂紋以及初始裂紋由于應(yīng)力集中開始擴(kuò)展,試件出現(xiàn)擴(kuò)容現(xiàn)象,因而聲發(fā)射比較活躍,并出現(xiàn)高幅值的能量率[16];隨著應(yīng)力增加,應(yīng)變開始快速增大,試件進(jìn)入塑性變形階段,加載應(yīng)力上升速度減緩,微裂紋的衍生和擴(kuò)展逐漸趨于穩(wěn)定,振鈴計(jì)數(shù)率相對(duì)塑性變形前期稍微減小并最終趨于一個(gè)穩(wěn)定值,能量率顯著降低。王祖蔭[22]認(rèn)為:材料處于穩(wěn)定的位錯(cuò)運(yùn)動(dòng)時(shí),內(nèi)部處于最低能狀態(tài),因而產(chǎn)生高頻度而低(強(qiáng)度)幅值的聲發(fā)射現(xiàn)象。到了加載后期,應(yīng)力逐漸達(dá)到峰值強(qiáng)度,微裂紋逐漸擴(kuò)展并貫通,AE信號(hào)主要由微裂紋貫通及破壞面晶?;飘a(chǎn)生并有所增強(qiáng),但相對(duì)加載前期的聲發(fā)射計(jì)數(shù)頻度和能量率,特別是能量率,已大為降低。

    4.2損傷演化分析

    圖7 煤巖試驗(yàn)擬合結(jié)果[13]Fig.7 Fitting result of coal testing[13]

    文獻(xiàn)[12-13]利用聲發(fā)射研究單軸壓縮下煤損傷本構(gòu)模型時(shí),用峰后殘余應(yīng)力與峰值應(yīng)力之比,并將1-σc/σm(σc為殘余應(yīng)力,σm為峰值應(yīng)力)作為損傷修正因子,對(duì)損傷邊界作出了修正,圖7為修正后理論曲線與試驗(yàn)曲線的對(duì)比。需要指出的是,巖石壓縮試驗(yàn)時(shí),殘余強(qiáng)度的取值本身帶有一定的主觀性,導(dǎo)致試驗(yàn)曲線與理論曲線差異較大。同時(shí),由于巖石峰后繼續(xù)破壞過程中仍有聲發(fā)射信號(hào),峰后取值點(diǎn)不同,不僅導(dǎo)致?lián)p傷理論公式得到不同結(jié)果,而且也不能反映鹽巖的真實(shí)損傷特征,故本文理論公式中Nm均取值至峰值應(yīng)力點(diǎn)。

    根據(jù)文獻(xiàn)[2-3]對(duì)鹽巖損傷的研究,單軸循環(huán)荷載下,鹽巖的損傷終值為0.85~0.95,由于循環(huán)荷載下試件損傷一般較常規(guī)荷載下試件損傷更為嚴(yán)重,并結(jié)合已有鹽巖破壞后CT掃描試驗(yàn)結(jié)果,假定損傷修正系數(shù)η取值為0.80~0.85,得出基于振鈴計(jì)數(shù)和能量參數(shù)修正前、后的損傷變量D與應(yīng)變的關(guān)系如圖8。

    圖8 修正前、后損傷因子D與應(yīng)變的關(guān)系Fig.8 Relationship between damage variable D andstrain before and after modification

    由于聲發(fā)射能量主要集中于加載前期,振鈴計(jì)數(shù)在整個(gè)加載過程比較均勻分布,因而計(jì)算所得的基于振鈴計(jì)數(shù)的理論損傷值要<基于能量的損傷值。式(8)和式(9)分別是基于圖8數(shù)據(jù)擬合得到的1#和2#鹽巖單軸壓縮荷載下修正后的基于聲發(fā)射振鈴計(jì)數(shù)和能量的損傷演化方程。

    1#鹽巖:

    (8)

    2#鹽巖:

    (9)

    將其分別代入式(7)即可得單軸壓縮荷載下鹽巖一維損傷的本構(gòu)方程。

    由圖8可知,基于聲發(fā)射的鹽巖損傷演化大致分為2個(gè)階段:

    (1) 初始損傷階段,主要對(duì)應(yīng)鹽巖的彈性變形階段。損傷理論認(rèn)為,此時(shí)微裂紋和裂隙的發(fā)展較少,變形屬于可恢復(fù)變形,因此,損傷變量應(yīng)該是趨于0的。然而,從圖8可以很清楚地看出,鹽巖在這一階段的損傷變量變化幅度較大,占整個(gè)過程的20%~40%,損傷值在壓縮變形早期增加較快。究其原因:一是受內(nèi)部晶粒之間的膠結(jié)情況以及其他物質(zhì)成分影響,鹽巖的強(qiáng)度、致密程度以及內(nèi)部結(jié)構(gòu)特征與其他巖石具有一定的差異性;二是由于鹽巖在這一階段并不是完全單純的彈性體,原始缺陷的壓密閉合和新裂紋的擴(kuò)展會(huì)產(chǎn)生大量的聲發(fā)射信號(hào)并產(chǎn)生不可逆變形,因而導(dǎo)致其初始損傷變量迅速增加。

    (2) 損傷逐步穩(wěn)定演化發(fā)展階段。從塑性變形前期一直延續(xù)到峰值應(yīng)力出現(xiàn),損傷變量增長幅度逐漸減小并最終趨于穩(wěn)定增加,鹽巖的微裂紋開始擴(kuò)展,同時(shí)產(chǎn)生新的微裂紋,并緩慢發(fā)展貫通,此階段損傷是連續(xù)穩(wěn)定的。由于鹽巖的延展性強(qiáng),密集及高幅值的聲發(fā)射信號(hào)發(fā)生在屈服階段,在破壞斷裂時(shí)產(chǎn)生可察覺但幅值不高的聲發(fā)射信號(hào)[22],因而在破壞階段,損傷變量變化不顯著,即沒有明顯的損傷加速階段出現(xiàn)。

    根據(jù)基于聲發(fā)射參數(shù)的損傷模型理論公式(5)和修正公式(7),得到基于聲發(fā)射參數(shù)值的應(yīng)力-應(yīng)變修正前后的曲線,并將其與試驗(yàn)結(jié)果作對(duì)比(見圖9)。

    圖9 修正前、后理論曲線與試驗(yàn)結(jié)果對(duì)比Fig.9 Comparison between theoretical curves andtest results before and after modification

    從圖9可以看出,未經(jīng)修正的聲發(fā)射模型與試驗(yàn)情況相差較大,得到的峰值應(yīng)力值也明顯小于試驗(yàn)所測得的峰值應(yīng)力;而經(jīng)過修正后的理論曲線在整個(gè)階段都能較好地與試驗(yàn)曲線吻合,修正后的理論峰值應(yīng)力值也與試驗(yàn)值非常接近。對(duì)比圖7和圖9理論曲線與試驗(yàn)曲線的吻合度可知,本文針對(duì)損傷終值提出的修正系數(shù)η及取值是有效的??梢?,聲發(fā)射作為鹽巖變形破壞過程中的伴生現(xiàn)象,如果獲得破壞時(shí)的真實(shí)損傷,便能夠較好地用來描述鹽巖壓縮狀態(tài)下的應(yīng)力-應(yīng)變損傷演化發(fā)展特征。

    由圖9對(duì)修正后基于振鈴計(jì)數(shù)和能量的曲線作比較,發(fā)現(xiàn)相對(duì)于聲發(fā)射能量參數(shù)的理論模型,基于聲發(fā)射振鈴計(jì)數(shù)的理論模型與試驗(yàn)曲線更為吻合,表明基于振鈴計(jì)數(shù)參數(shù)的模型比基于能量參數(shù)的模型更適合用于描述鹽巖應(yīng)力-應(yīng)變損傷全過程。Heiple和Carpenter[9]認(rèn)為:振鈴計(jì)數(shù)與材料中位錯(cuò)運(yùn)動(dòng)、夾雜物、第二相粒子剝離斷裂及裂紋擴(kuò)展所釋放的應(yīng)變能呈一定的相關(guān)關(guān)系,因而在一定程度上能夠更好地反映材料的損傷演化過程。

    4.3討論

    巖石損傷量的準(zhǔn)確試驗(yàn)測試具有很大困難。為準(zhǔn)確測定損傷,任建喜、葛修潤等[23-24]利用CT對(duì)巖石損傷分布特征、變形破壞過程中的損傷閾值等進(jìn)行了分析,取得了豐富研究成果。本文僅根據(jù)現(xiàn)有的鹽巖損傷研究文獻(xiàn)及已有CT掃描結(jié)果,對(duì)損傷修正系數(shù)取值進(jìn)行假定,未能實(shí)現(xiàn)聲發(fā)射與CT掃描的同步試驗(yàn)測試。鑒于CT掃描可為準(zhǔn)確獲得巖石破壞后各橫截面損傷及整體損傷量值提供有效測試手段,若能利用CT掃描獲得真實(shí)的損傷修正系數(shù),得到的本構(gòu)模型將更能反映實(shí)際損傷演化特征。筆者所在實(shí)驗(yàn)室已引進(jìn)了加載式工業(yè)CT,對(duì)巖石變形破壞全過程的損傷演化特征進(jìn)行CT掃描與聲發(fā)射三維定位的同步實(shí)時(shí)監(jiān)測將在今后進(jìn)行系統(tǒng)深入研究。

    5結(jié)論

    (1) 單軸壓縮下,鹽巖的聲發(fā)射振鈴計(jì)數(shù)率在整個(gè)試驗(yàn)過程維持在一個(gè)較高的水平,整體變化幅度不大,能量率在彈性極限附近達(dá)到最大值,塑性變形后期至破壞階段均呈現(xiàn)高頻低幅值的聲發(fā)射現(xiàn)象。

    (2) 變形初期階段,損傷變量值增長速率較快;基于能量計(jì)算得到的損傷變量值要大于基于振鈴計(jì)數(shù)的損傷變量值;引入修正系數(shù)前的理論曲線與試驗(yàn)曲線相差較大,且峰值應(yīng)力明顯小于試驗(yàn)值,針對(duì)損傷終值,引入修正系數(shù)后的理論曲線與試驗(yàn)曲線吻合較好,能更好地模擬鹽巖應(yīng)力-應(yīng)變損傷全過程。

    (3) 基于振鈴計(jì)數(shù)參數(shù)的模型比基于能量參數(shù)的理論模型計(jì)算得到的理論曲線更接近于試驗(yàn)結(jié)果,說明采用振鈴計(jì)數(shù)這一聲發(fā)射參數(shù)能夠較好地模擬鹽巖的應(yīng)力-應(yīng)變損傷過程。

    參考文獻(xiàn):

    [1]XIE He-ping,LIU Jian-feng,JU Yang,etal. Fractal Property of Spatial Distribution of Acoustic Emissions during the Failure Process of Bedded Rock Salt[J]. International Journal of Rock Mechanics and Mining Sciences, 2011, 48(8): 1344-1351.

    [2]LIU Jian-feng,XIE He-ping,HOU Zheng-meng,etal. Damage Evolution of Rock Salt under Cyclic Loading in Uniaxial Tests[J]. Acta Geotechnica,2014, 9(1):153-160.

    [3]劉建鋒,裴建良,張茹,等. 基于多級(jí)荷載試驗(yàn)的巖石損傷模量探討[J]. 巖石力學(xué)與工程學(xué)報(bào),2012,31(增1):3145-3151. [4]郭印同,趙克烈,孫冠華,等.周期荷載下鹽巖的疲勞變形及損傷特性研究[J]. 巖土力學(xué),2011,32(5):1353-1359.

    [5]SCHULZE O, POPP T, KERN H. Development of Damage and Permeability in Deforming Rock Salt[J]. Engineering Geology, 2001, 61(2/3): 163-180.

    [6]STORMONT J C. In situ Gas Permeability Measurements to Delineate Damage in Rock Salt[J]. International Journal of Rock Mechanics and Mining Sciences, 1997, 34(7): 1055-1064.[7]OHTSU M. Acoustic Emission Characteristics in Concrete and Diagnostic Application [J]. Acoustic Emission, 1987, 6(2): 99-108.[8]OHTSU M. Rate Process Analysis of Acoustic Emission Activity in Core Test of Concrete[J]. Concrete, 1992 (20):143-153 .

    [9]HEIPLE C R,CARPENTER S H.Acoustic Emission from Dislocation Motion[M]. In: Acoustic Emission, New York: Gordon and Breach Science Publishers, 1983.

    [10]WADLEY H N G, SCRUBY C B, SPEAKE J H. Acoustic Emission for Physical Examination of Metals[J]. International Metals Reviews, 1980, 25(1): 41-64.

    [11]TANG C A,XU X H. Evolution and Propagation of Material Defects and Kaiser Effect Function[J]. Journal of Seismological Research,1990,13(2):203-213.

    [12]劉保縣,黃敬林,王澤云,等.單軸壓縮煤巖損傷演化及聲發(fā)射特征研究[J].巖石力學(xué)與工程學(xué)報(bào),2009,28(增1):3234-3238.

    [13]高保彬,李回貴. 基于聲發(fā)射參量的煤樣損傷模型研究[J].防災(zāi)減災(zāi)工程學(xué)報(bào),2014,34(1):101-106.

    [14]楊永杰,王德超,郭明福,等. 基于三軸壓縮聲發(fā)射試驗(yàn)的巖石損傷特征研究[J].巖石力學(xué)與工程學(xué)報(bào),2014,33(1):98-104.

    [15]謝和平. 巖石混凝土損傷力學(xué)[M]. 徐州:中國礦業(yè)大學(xué)出版社,1990.

    [16]姜德義,陳結(jié),任松,等. 鹽巖單軸應(yīng)變率效應(yīng)與聲發(fā)射特征試驗(yàn)研究[J].巖石力學(xué)與工程學(xué)報(bào),2012,31(2):326-336.

    [17]任松,白月明,姜德義,等. 周期荷載作用下鹽巖聲發(fā)射特征試驗(yàn)研究[J].巖土力學(xué),2012,33(6):1613-1618.

    [18]李浩然,楊春和,劉玉剛,等.單軸荷載作用下鹽巖聲波與聲發(fā)射特征試驗(yàn)研究[J].巖石力學(xué)與工程學(xué)報(bào),2014,33(10):2107-2116.

    [19]KACHANOV L M. On the Time to Failure under Creep Conditions [J]. Izv AN SSSR Otd Tekhn Nauk,1958, 8(2): 26-31.

    [20]LEMAITURE J. Evaluation of Dissipation and Damage in Metals Submitted to Dynamic Loading[C]∥Proceedings of the ICM-1. Kyoto, Japan, 1972:540-549.

    [21]鄭雅麗,張華賓,王芝銀,等. 含雜質(zhì)鹽巖力學(xué)特性對(duì)比試驗(yàn)研究[J].煤炭學(xué)報(bào),2012,37(1):17-20.

    [22]王祖蔭.聲發(fā)射技術(shù)基礎(chǔ)[M].濟(jì)南:山東科學(xué)技術(shù)出版社,1990.

    [23]任建喜,葛修潤.單軸壓縮巖石損傷演化細(xì)觀機(jī)理及其本構(gòu)模型研究[J].巖石力學(xué)與工程學(xué)報(bào),2001,20(4):425-431.

    [24]葛修潤,任建喜,蒲毅彬,等. 巖石疲勞損傷擴(kuò)展規(guī)律CT細(xì)觀分析初探[J].巖土工程學(xué)報(bào),2001,23(2):191-195.

    (編輯:趙衛(wèi)兵)

    Acoustic Emission Characteristics and Damage Evolution ofRock Salt under Uniaxial Compression

    ZHOU Zhi-wei , LIU Jian-feng, ZOU Hang, ZHUO Yue, XU Yang-meng-di

    (State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu610065, China)

    Abstract:In order to explore the relation between acoustic emission phenomenon and internal damage in the failure process of compression deformation of rock salt, we conducted uniaxial compression test on salt rock from Anning, Yunnan province using MTS815 rock mechanics test system of Sichuan University. The damage evolution model of rock salt was established with acoustic emission parameters based on damage theory and the damage variable D was modified by coefficient of correction. The research results show that the maximum values of ring count rate and energy rate appeared near the elastic limit and the phenomenon of high frequency and low amplitude appeared later; the damage value based on AE increased rapidly in the early stage of compressive deformation, and the damage value based on AE accumulated ring counts is smaller than that based on cumulative energy; the revised damage model curve is in agreement with the whole process of test curve and could reflect the damage evolution process of rock salt well.Simulation using ring count can reflect the damage evolution characteristics of rock salt better than using energy.

    Key words:rock salt;acoustic emission;damage model ;damage variable;coefficient of correction

    中圖分類號(hào):TU45

    文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1001-5485(2016)05-0063-06

    doi:10.11988/ckyyb.201501682016,33(05):63-68,75

    作者簡介:周志威(1990-),男,四川樂山人,碩士研究生,主要從事巖石力學(xué)與工程方面的研究工作,(電話)15882485585(電子信箱)154762469@qq.com。通訊作者:劉建鋒(1979-),男,河南新鄉(xiāng)人,副教授,博士,研究方向?yàn)閹r石力學(xué)與工程,(電話)13550106392(電子信箱)ljfscu@163.com。

    基金項(xiàng)目:國家自然科學(xué)基金項(xiàng)目(51120145001,51104101,51374148)

    收稿日期:2015-03-10;修回日期:2015-04-14

    猜你喜歡
    聲發(fā)射鹽巖
    水熱綜合作用下鈣芒硝鹽巖強(qiáng)度等參數(shù)的衰減規(guī)律研究*
    鹽巖路基工程特性研究進(jìn)展
    公路工程(2021年6期)2021-02-14 12:34:06
    鹽巖巴西劈裂損傷愈合特性實(shí)驗(yàn)研究
    鹽巖儲(chǔ)庫腔底堆積物空隙體積試驗(yàn)與計(jì)算
    聲發(fā)射技術(shù)在壓力容器中的應(yīng)用與發(fā)展
    濾波器在聲發(fā)射信號(hào)噪聲預(yù)處理中的應(yīng)用
    典型鹽巖力學(xué)特性分析
    修正的鹽巖擴(kuò)容模型及擴(kuò)容界限研究
    基于LPC1768的聲發(fā)射數(shù)據(jù)采集系統(tǒng)設(shè)計(jì)
    聲發(fā)射信號(hào)處理關(guān)鍵技術(shù)研究
    国产麻豆成人av免费视频| 伊人久久精品亚洲午夜| 亚洲中文字幕一区二区三区有码在线看| 久久久久久久久久久免费av| 男人舔奶头视频| 神马国产精品三级电影在线观看| 色综合亚洲欧美另类图片| 国产精品女同一区二区软件| 久久精品国产亚洲网站| 成人午夜精彩视频在线观看| 日韩人妻高清精品专区| 最近最新中文字幕免费大全7| 少妇猛男粗大的猛烈进出视频 | 91精品伊人久久大香线蕉| 岛国毛片在线播放| 久久精品人妻少妇| 亚洲18禁久久av| 日日干狠狠操夜夜爽| 精品一区二区三区视频在线| 免费人成在线观看视频色| 在线a可以看的网站| 26uuu在线亚洲综合色| 综合色丁香网| 99久国产av精品| 色5月婷婷丁香| 亚洲精品乱码久久久久久按摩| 亚洲av男天堂| av免费观看日本| 亚洲成色77777| 久久精品影院6| 久久久久久久午夜电影| 亚洲av电影在线观看一区二区三区 | 少妇被粗大猛烈的视频| 欧美一区二区精品小视频在线| 色网站视频免费| 在线a可以看的网站| 亚洲成色77777| 国产精品综合久久久久久久免费| 国产极品天堂在线| 中文欧美无线码| 日韩精品青青久久久久久| 久久久久久久久久成人| 高清在线视频一区二区三区 | 久久久欧美国产精品| 欧美日本视频| 成人一区二区视频在线观看| 免费一级毛片在线播放高清视频| 亚洲最大成人中文| 欧美三级亚洲精品| 久久欧美精品欧美久久欧美| 色噜噜av男人的天堂激情| 爱豆传媒免费全集在线观看| 久久久久久大精品| 一二三四中文在线观看免费高清| 一区二区三区四区激情视频| 99热全是精品| 国产探花在线观看一区二区| 亚洲精品久久久久久婷婷小说 | 国产av一区在线观看免费| 一个人看的www免费观看视频| 26uuu在线亚洲综合色| 欧美激情在线99| 特大巨黑吊av在线直播| 国产老妇女一区| 亚洲真实伦在线观看| 国产黄a三级三级三级人| 毛片一级片免费看久久久久| 欧美极品一区二区三区四区| 国产熟女欧美一区二区| 国产一区二区在线av高清观看| 大香蕉97超碰在线| 天堂网av新在线| 日本黄大片高清| 又粗又硬又长又爽又黄的视频| 麻豆成人午夜福利视频| 天堂av国产一区二区熟女人妻| 欧美一区二区精品小视频在线| 九色成人免费人妻av| 国产精品嫩草影院av在线观看| 日本一本二区三区精品| 亚洲av中文av极速乱| 欧美性猛交黑人性爽| 精品无人区乱码1区二区| 欧美高清成人免费视频www| ponron亚洲| 亚洲国产欧洲综合997久久,| 久久久久性生活片| 日韩,欧美,国产一区二区三区 | 精品人妻一区二区三区麻豆| av天堂中文字幕网| 高清av免费在线| 亚洲成人久久爱视频| 国产午夜精品一二区理论片| 日日啪夜夜撸| 日本av手机在线免费观看| 如何舔出高潮| 国产在线男女| 国产免费一级a男人的天堂| 变态另类丝袜制服| 热99re8久久精品国产| 亚洲国产日韩欧美精品在线观看| 搡女人真爽免费视频火全软件| 国产精品久久视频播放| av专区在线播放| 蜜臀久久99精品久久宅男| 偷拍熟女少妇极品色| 亚洲欧美精品综合久久99| 国产精品.久久久| 婷婷六月久久综合丁香| 日本欧美国产在线视频| 成人一区二区视频在线观看| av国产免费在线观看| 六月丁香七月| 成人午夜精彩视频在线观看| 美女大奶头视频| 少妇猛男粗大的猛烈进出视频 | 午夜福利成人在线免费观看| 人妻夜夜爽99麻豆av| 久久久久久久午夜电影| 亚洲国产成人一精品久久久| 午夜福利在线观看免费完整高清在| 99在线人妻在线中文字幕| 中文字幕精品亚洲无线码一区| 中文资源天堂在线| 美女国产视频在线观看| 最近的中文字幕免费完整| 91狼人影院| 日韩中字成人| 亚洲精品国产av成人精品| 青春草视频在线免费观看| 亚洲av中文av极速乱| 日本爱情动作片www.在线观看| 色尼玛亚洲综合影院| 大话2 男鬼变身卡| 国产精品久久久久久精品电影小说 | 国产一级毛片七仙女欲春2| 国产精品三级大全| 啦啦啦啦在线视频资源| 搞女人的毛片| 国产私拍福利视频在线观看| 热99re8久久精品国产| 不卡视频在线观看欧美| 久久精品影院6| 建设人人有责人人尽责人人享有的 | 国产亚洲5aaaaa淫片| 深爱激情五月婷婷| 直男gayav资源| 国产午夜福利久久久久久| a级毛片免费高清观看在线播放| 国产精品伦人一区二区| 大话2 男鬼变身卡| 丝袜美腿在线中文| a级一级毛片免费在线观看| 国产真实伦视频高清在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产大屁股一区二区在线视频| 纵有疾风起免费观看全集完整版 | 一边摸一边抽搐一进一小说| 日韩欧美在线乱码| eeuss影院久久| 亚洲欧美清纯卡通| 波多野结衣巨乳人妻| 哪个播放器可以免费观看大片| 18禁在线无遮挡免费观看视频| 青春草视频在线免费观看| 自拍偷自拍亚洲精品老妇| 久久韩国三级中文字幕| 大又大粗又爽又黄少妇毛片口| 亚洲真实伦在线观看| 国产成人a∨麻豆精品| 亚洲在线观看片| 亚洲精品,欧美精品| 欧美成人a在线观看| 中文字幕久久专区| 欧美日韩在线观看h| 欧美xxxx黑人xx丫x性爽| 联通29元200g的流量卡| 亚洲怡红院男人天堂| 能在线免费观看的黄片| 小说图片视频综合网站| 青青草视频在线视频观看| 国产av一区在线观看免费| 久久久久久九九精品二区国产| 免费观看性生交大片5| 日韩欧美国产在线观看| 天天躁夜夜躁狠狠久久av| 国产黄片视频在线免费观看| 国产伦精品一区二区三区四那| 日韩中字成人| 秋霞在线观看毛片| 国产熟女欧美一区二区| 岛国在线免费视频观看| 精品久久久久久久人妻蜜臀av| 3wmmmm亚洲av在线观看| 亚洲精品国产av成人精品| 在线播放无遮挡| 亚洲av成人精品一二三区| 成人性生交大片免费视频hd| 国产精品一区二区在线观看99 | av视频在线观看入口| 亚洲在线自拍视频| 26uuu在线亚洲综合色| 国产老妇伦熟女老妇高清| 久久精品国产亚洲av涩爱| 亚洲精品乱久久久久久| 18禁动态无遮挡网站| 寂寞人妻少妇视频99o| 午夜福利在线观看免费完整高清在| 18禁在线播放成人免费| 国产淫语在线视频| 亚洲av.av天堂| 一级黄片播放器| 国产白丝娇喘喷水9色精品| 亚洲内射少妇av| 午夜福利成人在线免费观看| 自拍偷自拍亚洲精品老妇| 插阴视频在线观看视频| 国产精品一区二区三区四区免费观看| 九九热线精品视视频播放| 99久久成人亚洲精品观看| 亚洲av中文av极速乱| av国产久精品久网站免费入址| av专区在线播放| 热99在线观看视频| 国产中年淑女户外野战色| 色网站视频免费| 深爱激情五月婷婷| 麻豆一二三区av精品| 国产精品久久久久久精品电影| 少妇人妻精品综合一区二区| 国内精品美女久久久久久| 一区二区三区四区激情视频| 日韩成人伦理影院| 精品久久国产蜜桃| 又爽又黄a免费视频| 免费看a级黄色片| 99久久精品国产国产毛片| 免费看光身美女| 国产国拍精品亚洲av在线观看| 欧美极品一区二区三区四区| 91aial.com中文字幕在线观看| 国产精品熟女久久久久浪| 日韩,欧美,国产一区二区三区 | 亚洲最大成人中文| 成人漫画全彩无遮挡| 99热这里只有是精品50| 日韩制服骚丝袜av| 男人狂女人下面高潮的视频| 日韩av在线免费看完整版不卡| 欧美zozozo另类| 十八禁国产超污无遮挡网站| 久久婷婷人人爽人人干人人爱| 国产视频内射| 国产三级中文精品| 1024手机看黄色片| 熟女电影av网| 国产精品久久视频播放| 国产精品.久久久| 国产高清国产精品国产三级 | 国产亚洲最大av| 精品久久久久久久久久久久久| av在线观看视频网站免费| 久久久色成人| 亚洲成色77777| 18禁在线播放成人免费| 联通29元200g的流量卡| 国产精品人妻久久久影院| 久久久久久久国产电影| 美女大奶头视频| 国产麻豆成人av免费视频| 高清午夜精品一区二区三区| 99久久无色码亚洲精品果冻| 九草在线视频观看| 搡老妇女老女人老熟妇| 91狼人影院| 精品久久国产蜜桃| 免费人成在线观看视频色| 久久久久久久久大av| 欧美成人精品欧美一级黄| 91午夜精品亚洲一区二区三区| 日韩欧美精品免费久久| .国产精品久久| 一夜夜www| 天堂av国产一区二区熟女人妻| av免费在线看不卡| 精品人妻一区二区三区麻豆| 又粗又硬又长又爽又黄的视频| 国产精品永久免费网站| 免费播放大片免费观看视频在线观看 | av专区在线播放| 丰满人妻一区二区三区视频av| 午夜福利成人在线免费观看| www日本黄色视频网| 99热6这里只有精品| 亚洲av日韩在线播放| 亚洲av熟女| 日日摸夜夜添夜夜添av毛片| 国产亚洲最大av| 青春草视频在线免费观看| 七月丁香在线播放| 欧美丝袜亚洲另类| 亚洲精品亚洲一区二区| 午夜爱爱视频在线播放| 全区人妻精品视频| 亚洲精品国产av成人精品| 国产白丝娇喘喷水9色精品| 视频中文字幕在线观看| 国产伦理片在线播放av一区| 中文字幕熟女人妻在线| 国产又色又爽无遮挡免| 亚洲精品国产成人久久av| 国产探花在线观看一区二区| 长腿黑丝高跟| 国产淫语在线视频| 男女下面进入的视频免费午夜| 国产成人91sexporn| 寂寞人妻少妇视频99o| 亚洲自拍偷在线| 天美传媒精品一区二区| 乱码一卡2卡4卡精品| 亚洲精品久久久久久婷婷小说 | 亚洲国产最新在线播放| 天堂网av新在线| 日韩中字成人| 亚洲人与动物交配视频| 免费av观看视频| 99热网站在线观看| 国内精品宾馆在线| 亚洲精品一区蜜桃| 小蜜桃在线观看免费完整版高清| 国产乱人视频| 欧美激情在线99| 欧美一级a爱片免费观看看| 自拍偷自拍亚洲精品老妇| videossex国产| 国产伦精品一区二区三区四那| 日本午夜av视频| 国产在视频线精品| 亚洲国产日韩欧美精品在线观看| 国产私拍福利视频在线观看| 建设人人有责人人尽责人人享有的 | 中文字幕免费在线视频6| 国内精品美女久久久久久| 成人特级av手机在线观看| 直男gayav资源| 久久久久精品久久久久真实原创| 亚洲精品乱码久久久v下载方式| 91av网一区二区| 美女高潮的动态| 久久久午夜欧美精品| 人妻少妇偷人精品九色| 91午夜精品亚洲一区二区三区| 能在线免费观看的黄片| 亚洲精品自拍成人| 少妇熟女欧美另类| av播播在线观看一区| 麻豆精品久久久久久蜜桃| 久久久午夜欧美精品| 中文字幕精品亚洲无线码一区| 91午夜精品亚洲一区二区三区| 精品人妻熟女av久视频| 2022亚洲国产成人精品| 国产亚洲一区二区精品| 亚洲成人精品中文字幕电影| 成人二区视频| 最新中文字幕久久久久| 又爽又黄a免费视频| 美女国产视频在线观看| 国产成人精品一,二区| 美女黄网站色视频| 精品国内亚洲2022精品成人| 欧美性猛交黑人性爽| 午夜视频国产福利| 99九九线精品视频在线观看视频| 久久精品国产亚洲av涩爱| 久久综合国产亚洲精品| 九草在线视频观看| av国产久精品久网站免费入址| 成人亚洲欧美一区二区av| 丰满少妇做爰视频| 精品一区二区三区人妻视频| 最近2019中文字幕mv第一页| 男人舔奶头视频| 男人和女人高潮做爰伦理| www.av在线官网国产| 狂野欧美激情性xxxx在线观看| 久久这里有精品视频免费| 久久精品综合一区二区三区| av在线亚洲专区| 日本免费a在线| 熟妇人妻久久中文字幕3abv| 欧美激情在线99| 国产白丝娇喘喷水9色精品| 嫩草影院入口| 高清午夜精品一区二区三区| 亚洲国产色片| 亚洲精品日韩在线中文字幕| 啦啦啦啦在线视频资源| 午夜a级毛片| 嘟嘟电影网在线观看| 国产乱人视频| 免费观看的影片在线观看| 亚洲三级黄色毛片| 国产亚洲av嫩草精品影院| 色5月婷婷丁香| 综合色av麻豆| 丝袜喷水一区| 国产一区二区在线观看日韩| 国产精品久久久久久久久免| 日韩人妻高清精品专区| 亚洲va在线va天堂va国产| 亚洲国产精品合色在线| 国国产精品蜜臀av免费| 欧美另类亚洲清纯唯美| 深夜a级毛片| 国产成人aa在线观看| 日本黄大片高清| 男插女下体视频免费在线播放| 亚洲最大成人中文| 美女xxoo啪啪120秒动态图| 看黄色毛片网站| 亚洲av福利一区| 国产91av在线免费观看| 永久免费av网站大全| 少妇猛男粗大的猛烈进出视频 | 18禁在线播放成人免费| 国产激情偷乱视频一区二区| 国产精品久久久久久久电影| 91久久精品国产一区二区成人| 91精品一卡2卡3卡4卡| 精品午夜福利在线看| 高清毛片免费看| 少妇人妻一区二区三区视频| 最近2019中文字幕mv第一页| 天堂中文最新版在线下载 | 大香蕉久久网| 一个人观看的视频www高清免费观看| 最近中文字幕2019免费版| 一本一本综合久久| 久久久午夜欧美精品| 亚洲精品456在线播放app| 精品国产三级普通话版| 日本一本二区三区精品| 狂野欧美激情性xxxx在线观看| 搞女人的毛片| 能在线免费看毛片的网站| 午夜福利在线观看免费完整高清在| 水蜜桃什么品种好| 亚洲精品成人久久久久久| 99久久成人亚洲精品观看| 亚洲欧美日韩东京热| 亚洲av一区综合| 三级经典国产精品| 日韩亚洲欧美综合| 中文字幕av成人在线电影| 成年免费大片在线观看| 久久99热这里只有精品18| 黄色欧美视频在线观看| 欧美一区二区亚洲| 人妻制服诱惑在线中文字幕| 一级毛片我不卡| 国产高潮美女av| 麻豆精品久久久久久蜜桃| 少妇丰满av| 国产综合懂色| 日韩强制内射视频| 女人久久www免费人成看片 | 蜜臀久久99精品久久宅男| 欧美极品一区二区三区四区| 国产精品人妻久久久影院| 97超碰精品成人国产| 日日干狠狠操夜夜爽| 在线天堂最新版资源| 日韩 亚洲 欧美在线| 久久国产乱子免费精品| 国产又黄又爽又无遮挡在线| 精品久久国产蜜桃| 看免费成人av毛片| 国产精品嫩草影院av在线观看| 汤姆久久久久久久影院中文字幕 | 夜夜爽夜夜爽视频| 听说在线观看完整版免费高清| 国产白丝娇喘喷水9色精品| 一本一本综合久久| 欧美日韩在线观看h| 高清在线视频一区二区三区 | av黄色大香蕉| 久久久精品94久久精品| 国产精品.久久久| 国产精品嫩草影院av在线观看| av在线观看视频网站免费| 亚洲不卡免费看| 国产男人的电影天堂91| 特级一级黄色大片| 亚洲综合精品二区| 久久99精品国语久久久| 亚洲精品一区蜜桃| 亚洲精品国产成人久久av| 少妇高潮的动态图| 男人和女人高潮做爰伦理| 国产精品福利在线免费观看| 老师上课跳d突然被开到最大视频| 日日撸夜夜添| a级一级毛片免费在线观看| 女人十人毛片免费观看3o分钟| 插逼视频在线观看| 建设人人有责人人尽责人人享有的 | 嫩草影院新地址| 日本黄色视频三级网站网址| 日韩一区二区三区影片| 亚洲综合色惰| kizo精华| 午夜老司机福利剧场| 久久欧美精品欧美久久欧美| 久久久久久久久久久丰满| 精品久久国产蜜桃| 国产白丝娇喘喷水9色精品| 国产精品嫩草影院av在线观看| 亚洲成人精品中文字幕电影| 久久久久久久亚洲中文字幕| 一个人免费在线观看电影| 插逼视频在线观看| 国产69精品久久久久777片| 免费电影在线观看免费观看| 91狼人影院| 日本一二三区视频观看| 日韩欧美精品v在线| 国产精品蜜桃在线观看| 九色成人免费人妻av| 国内精品一区二区在线观看| 精品一区二区三区视频在线| 看非洲黑人一级黄片| 精品国内亚洲2022精品成人| 丝袜美腿在线中文| 久久久亚洲精品成人影院| 99热这里只有精品一区| 久久久久久久午夜电影| 欧美最新免费一区二区三区| 秋霞在线观看毛片| 我的女老师完整版在线观看| 欧美日本视频| 春色校园在线视频观看| 深夜a级毛片| 欧美一级a爱片免费观看看| 精品免费久久久久久久清纯| 午夜福利网站1000一区二区三区| 欧美另类亚洲清纯唯美| 欧美一区二区国产精品久久精品| av天堂中文字幕网| 麻豆精品久久久久久蜜桃| 日韩大片免费观看网站 | 男女边吃奶边做爰视频| 亚洲av二区三区四区| 国产中年淑女户外野战色| 亚洲综合色惰| 欧美人与善性xxx| 亚洲欧美日韩高清专用| 久久久久精品久久久久真实原创| 欧美变态另类bdsm刘玥| 97热精品久久久久久| 婷婷色麻豆天堂久久 | 国产高清国产精品国产三级 | 国产乱来视频区| 亚洲一级一片aⅴ在线观看| 特级一级黄色大片| 麻豆成人午夜福利视频| 亚洲av成人精品一二三区| 午夜免费男女啪啪视频观看| 一本久久精品| av免费观看日本| 久久这里有精品视频免费| 国产精品久久视频播放| 国产精品福利在线免费观看| 午夜精品一区二区三区免费看| 最近中文字幕高清免费大全6| 国产女主播在线喷水免费视频网站 | 久久久a久久爽久久v久久| 精品免费久久久久久久清纯| 久久久久久九九精品二区国产| 国产久久久一区二区三区| 国产视频首页在线观看| 国产精品久久久久久久电影| 国产v大片淫在线免费观看| 麻豆国产97在线/欧美| 大香蕉久久网| 韩国av在线不卡| 亚洲精品乱久久久久久| 国产毛片a区久久久久| 精品人妻熟女av久视频| 日日撸夜夜添| 国产精品乱码一区二三区的特点| 成人性生交大片免费视频hd| 久久久久久伊人网av| 两性午夜刺激爽爽歪歪视频在线观看| 国产片特级美女逼逼视频| 亚洲怡红院男人天堂| 国产伦一二天堂av在线观看| 亚洲欧美中文字幕日韩二区| 国产成人91sexporn| 不卡视频在线观看欧美| 51国产日韩欧美| 国产精品,欧美在线| 国内精品一区二区在线观看| 国产中年淑女户外野战色| 国产精品人妻久久久影院| 国产白丝娇喘喷水9色精品| 成年免费大片在线观看| 变态另类丝袜制服| 小说图片视频综合网站| 蜜桃亚洲精品一区二区三区| 国产一区二区亚洲精品在线观看| 色综合色国产| 久久精品国产亚洲网站| АⅤ资源中文在线天堂| 一级黄色大片毛片|