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

    電阻斷層成像、應(yīng)力波及阻抗儀3種無損檢測(cè)方法對(duì)活立木腐朽程度的定量檢測(cè)*

    2017-04-27 06:44:11岳小泉王立海王興龍榮賓賓葛曉雯劉澤旭陳清耀
    林業(yè)科學(xué) 2017年3期

    岳小泉 王立海 王興龍 榮賓賓 葛曉雯 劉澤旭 陳清耀

    (1.東北林業(yè)大學(xué)工程技術(shù)學(xué)院 哈爾濱 150040; 2.福建農(nóng)林大學(xué)交通與土木工程學(xué)院 福州 350002)

    ?

    電阻斷層成像、應(yīng)力波及阻抗儀3種無損檢測(cè)方法對(duì)活立木腐朽程度的定量檢測(cè)*

    岳小泉1,2王立海1王興龍1榮賓賓1葛曉雯1劉澤旭1陳清耀2

    (1.東北林業(yè)大學(xué)工程技術(shù)學(xué)院 哈爾濱 150040; 2.福建農(nóng)林大學(xué)交通與土木工程學(xué)院 福州 350002)

    【目的】應(yīng)用電阻斷層成像法、應(yīng)力波斷層成像法和阻抗儀法對(duì)活立木樹干腐朽程度進(jìn)行檢測(cè)和定量表征,判斷3種方法的可靠性,為野外立木健康情況檢測(cè)提供基礎(chǔ)數(shù)據(jù)和理論依據(jù)?!痉椒ā吭诤邶埥」枮I市東北林業(yè)大學(xué)試驗(yàn)林場(chǎng)內(nèi),選取水曲柳和北京楊各25株(15株腐朽,10株健康)共50株100個(gè)截面,采用PICUS Tree Tronic型樹木電阻斷層成像儀、Arbotom應(yīng)力波斷層成像系統(tǒng)和Resistograph針式阻抗儀檢測(cè)及木芯質(zhì)量損失率估測(cè)4種方法對(duì)斷面腐朽程度進(jìn)行檢測(cè)和定量表征。以木芯質(zhì)量損失率計(jì)算的樣木腐朽程度(Es)為真值,利用最小二乘法分別建立電阻測(cè)定的腐朽程度(Ed)、應(yīng)力波測(cè)定的腐朽程度(Ey)及阻抗儀測(cè)定的腐朽程度(Ez)與真值(Es)之間的線性關(guān)系?!窘Y(jié)果】電阻斷層成像、應(yīng)力波和阻抗儀3種無損檢測(cè)方法能在不同腐朽程度上表征活立木木材質(zhì)量損失率; 從整體上來看,電阻檢測(cè)和應(yīng)力波檢測(cè)結(jié)果與腐朽程度真值的擬合程度低于阻抗儀檢測(cè)結(jié)果與腐朽程度真值的擬合程度; 從不同腐朽階段來看,在Es<30%時(shí),Ed和Es的相關(guān)系數(shù)最高(R=0.823,P<0.01); 在30%≤Es<50%時(shí),Ey和Es的相關(guān)系數(shù)最高(R=0.658,P<0.01); 在Es≥50%時(shí),Ez和Es的相關(guān)系數(shù)最高(R=0.914,P<0.01)?!窘Y(jié)論】電阻斷層成像法對(duì)木材早期腐朽的檢測(cè)比較敏感,而應(yīng)力波斷層成像法則相對(duì)在木材腐朽稍嚴(yán)重時(shí)更準(zhǔn)確,阻抗儀法在腐朽各個(gè)階段的檢測(cè)結(jié)果都較為準(zhǔn)確,且在腐朽嚴(yán)重時(shí)最準(zhǔn)確,但會(huì)對(duì)木材造成微損傷。3種無損檢測(cè)方法均能有效檢測(cè)活立木腐朽,并有各自的特點(diǎn),在實(shí)際腐朽檢測(cè)中應(yīng)根據(jù)實(shí)際檢測(cè)需要來選擇。關(guān)鍵詞: 質(zhì)量損失率; 電阻斷層成像; 應(yīng)力波傳播速度; 阻抗儀阻力損失值; 腐朽程度

    木材無損檢測(cè)技術(shù)在避免對(duì)木材造成創(chuàng)傷和破壞的前提下,可以對(duì)木材腐朽缺陷進(jìn)行快速、準(zhǔn)確檢測(cè),便捷地獲得木材的特性指標(biāo)和內(nèi)部狀況,不僅可為立木評(píng)價(jià)和林區(qū)經(jīng)營(yíng)提供科學(xué)依據(jù),而且也能為優(yōu)化造材和合理加工提供重要參考。當(dāng)前,木材無損檢測(cè)方法很多,其中,應(yīng)用于活立木腐朽檢測(cè)的主要有應(yīng)力波、阻抗儀和電阻技術(shù)(Wangetal., 2004; 2008; 高珊等, 2013; 王立海等, 2001)。實(shí)際應(yīng)用表明,各種檢測(cè)方法都有一定的優(yōu)缺點(diǎn)。

    電阻斷層成像法在木材檢測(cè)中的應(yīng)用是近幾年才發(fā)展起來的,能便捷地對(duì)活立木缺陷進(jìn)行檢測(cè)。電阻斷層成像法在國(guó)外發(fā)展較早,Just等(1998)將其應(yīng)用于活立木腐朽檢測(cè)中,證明了電阻斷層成像檢測(cè)法的可行性。Brazee等(2011)將聲波層析成像和電阻抗成像EIT(electric impedance tomography)相結(jié)合,檢測(cè)并量化了活立木內(nèi)部的腐朽程度。徐速等(2006)、吳華橋等(2008)、周啟友等(2009)將ERT(electric resistance tomography)技術(shù)應(yīng)用于立木樹干水分分布和移動(dòng)研究,通過測(cè)量電阻值分析立木中的樹液流動(dòng)情況,并通過反演計(jì)算獲得了立木截面的二維電阻圖像,從而分析降雨、晝夜變化等因素對(duì)立木中水分變化的影響。

    應(yīng)力波可檢測(cè)出木質(zhì)材料的力學(xué)性能和內(nèi)部缺陷,是目前較常用的木材無損檢測(cè)技術(shù)之一(徐華東等, 2014a; 楊學(xué)春等, 2005; 2007)。針對(duì)應(yīng)力波斷層成像,徐華東等(2010)對(duì)40株旱柳(Salixmatsudana)進(jìn)行檢測(cè),研究應(yīng)力波在立木檢測(cè)中的傳播規(guī)律,并應(yīng)用應(yīng)力波二維成像技術(shù)對(duì)旱柳的安全狀況進(jìn)行了評(píng)估; Wang等(2007)分析傳感器數(shù)量對(duì)應(yīng)力波成像的擬合度和誤差率的影響,并指出圖像擬合度接近90%和誤差率在0.1左右時(shí),需12個(gè)傳感器才能滿足要求; 葛曉雯等(2014)聯(lián)合應(yīng)用應(yīng)力波和阻抗儀2種儀器對(duì)旱柳行道樹內(nèi)部腐朽的位置、程度等狀況進(jìn)行檢測(cè),結(jié)果發(fā)現(xiàn)2種檢測(cè)方法的對(duì)應(yīng)關(guān)系較好。

    阻抗儀法主要用于立木、木橋及木結(jié)構(gòu)等的探測(cè)。阻抗儀能夠以較高的效率獲得被檢測(cè)建筑木構(gòu)件、立木或木材內(nèi)部的缺陷情況(黃榮鳳等, 2007; 孫天用等, 2013; 2014),因此在美國(guó)、歐洲和日本,對(duì)古建筑木結(jié)構(gòu)進(jìn)行維修前的木結(jié)構(gòu)安全評(píng)價(jià)常采用阻力曲線圖。

    本文將野外測(cè)試和室內(nèi)試驗(yàn)相結(jié)合,對(duì)活立木樹干腐朽程度進(jìn)行檢測(cè)和定量表征。以木芯質(zhì)量損失率計(jì)算值作為樣木腐朽程度的真值,分析野外應(yīng)用電阻斷層成像法(下文簡(jiǎn)稱電阻法)、應(yīng)力波斷層成像法(下文簡(jiǎn)稱應(yīng)力波法)及阻抗儀法的檢測(cè)結(jié)果,判斷3種方法的可靠性,為野外立木健康情況檢測(cè)提供基礎(chǔ)數(shù)據(jù)和理論依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    研究區(qū)位于哈爾濱市東北林業(yè)大學(xué)試驗(yàn)林場(chǎng),地理坐標(biāo)為126°37′E,45°43′N,海拔136~140 m,地形略有起伏,南高北低,西高東低,坡度<5°,土地總面積43.95 hm2。研究區(qū)屬溫帶半濕潤(rùn)季風(fēng)氣候區(qū),年均溫3.6 ℃,7月最高溫36.4 ℃,1月最低溫-38.1 ℃,無霜期136天,≥10 ℃年積溫2 757 ℃,年降水量600 mm左右。原生植被為溝谷榆樹疏林草原,并于20世紀(jì)50年代末期和60年代初期進(jìn)行了人工林造林試驗(yàn),林場(chǎng)現(xiàn)有18種人工林分布在46個(gè)樣地內(nèi),每個(gè)樣地面積為0.5 hm2,每個(gè)樣地內(nèi)都栽植1種樹種。

    1.2 試驗(yàn)材料及設(shè)備

    2015年7月,在研究區(qū)內(nèi)選擇水曲柳(Fraxinusmandshurica)和北京楊(Populus×beijingensis)為測(cè)試對(duì)象。在樣地內(nèi)先目測(cè)可能存在內(nèi)部腐朽的活立木,主要觀察樹干上是否有空洞、枝葉是否枯落、樹皮是否有腐爛或外傷、樹干是否有臃腫或是否存在痂皮現(xiàn)象等。目測(cè)后,水曲柳和北京楊各選取15株可能存在內(nèi)部腐朽的活立木和10株健康活立木作為樣木。兩樹種樹齡均為50~60年,水曲柳胸徑為20~38 cm,北京楊胸徑為30~50 cm。

    試驗(yàn)主要采用3種儀器對(duì)立木進(jìn)行測(cè)試: PICUS Tree Tronic型樹木電阻斷層成像儀(德國(guó)Argus公司生產(chǎn))、Arbotom應(yīng)力波斷層成像系統(tǒng)(德國(guó)RINNTECH公司生產(chǎn))及型號(hào)為3450的Resistograph針式阻抗儀(德國(guó)FrankRinnIML公司生產(chǎn)),瑞典樹木生長(zhǎng)錐鉆,天津市泰斯特儀器有限公司生產(chǎn)的101-3A型鼓風(fēng)干燥機(jī)。

    1.3 測(cè)試方法

    (1)

    1.4 數(shù)據(jù)處理

    1.4.1 電阻斷層成像數(shù)據(jù)處理 圖1為電阻法測(cè)得的立木斷面電阻分布。該圖像表示樹木斷面的電阻率變化情況,紅色部分表示高電阻率區(qū),藍(lán)色部分表示低電阻率區(qū)。為準(zhǔn)確計(jì)算電阻斷層成像圖中的電阻值,先將彩色圖像進(jìn)行灰度化(灰階0-255),然后通過Matlab軟件將圖中的像素值轉(zhuǎn)化為電阻值,使電阻的數(shù)值顯示于圖像上,如圖2所示?,F(xiàn)定義電阻檢測(cè)的腐朽程度為Ed,公式如下:

    (2)

    式中:R0為同種樹所有健康立木斷面阻抗儀檢測(cè)方向上電阻率的平均值(Ω);Rd為某個(gè)腐朽斷面被檢測(cè)方向上的平均電阻率值(Ω)。

    圖1 立木斷面電阻分布Fig.1 Resistance distribution of cross-section

    1.4.2 應(yīng)力波數(shù)據(jù)處理 圖3為應(yīng)力波法測(cè)得的立木斷面波速分布圖。如果存在腐朽缺陷,則應(yīng)力波傳播速度就變小,相應(yīng)地在圖像中就呈現(xiàn)出顏色較深的區(qū)域(徐華東等, 2014a)。按照處理電阻斷層成像相同的圖像處理方式來計(jì)算應(yīng)力波傳播速度?,F(xiàn)定義應(yīng)力波檢測(cè)的腐朽程度為Ey,公式如下:

    (3)

    式中:Vj為同種樹所有健康立木斷面應(yīng)力波測(cè)試儀檢測(cè)方向上應(yīng)力波傳播速度的平均值(m·s-1);Vf為某腐朽斷面被檢測(cè)方向上的平均傳播速度(m·s-1)。

    1.4.3 阻抗儀數(shù)據(jù)處理 阻抗儀的檢測(cè)結(jié)果為阻力曲線圖,通過對(duì)阻力曲線圖進(jìn)行分析可實(shí)現(xiàn)腐朽程度的定量。在阻力曲線圖(圖4)上,橫坐標(biāo)表示

    圖2 電阻分布灰度圖及對(duì)應(yīng)的電阻值分布Fig.2 Grey-scale map and corresponding value of resistance distribution

    圖4 阻抗儀Resistograph測(cè)得的阻力曲線Fig.4 Resistance curve of cross-section by Resistograph

    探針鉆入木材的深度,縱坐標(biāo)表示探針?biāo)茏枇Φ南鄬?duì)大小(單位:resi)。木材的硬度、密度等物理力學(xué)特性影響著探針受到的阻力大小(黃榮鳳等, 2007; 安源等, 2008),當(dāng)活立木某處發(fā)生腐朽時(shí),該處木材的力學(xué)強(qiáng)度和密度就會(huì)顯著降低,探針鉆到該處時(shí)所受周圍木材的阻力就會(huì)降低,在阻力曲線圖上呈現(xiàn)為一個(gè)波谷,圖4中ABCD區(qū)域所包含的那段曲線表明活立木有腐朽。根據(jù)曲線圖上橫縱坐標(biāo)軸的實(shí)際意義,現(xiàn)定義腐朽程度E(單位resi),公式如下:

    (4)

    式中:L為下降段曲線在橫軸上投影的長(zhǎng)度(mm);H,h分別為下降段曲線左右兩邊AB和CD在縱軸上投影的長(zhǎng)度(resi);D為所測(cè)活立木橫截面的直徑(mm)。

    圖3 立木斷面應(yīng)力波波速分布Fig.3 The velocity distribution of stress wave about cross-section

    應(yīng)用此公式將每幅阻力曲線圖的腐朽程度(E)計(jì)算出來,找出所有斷面E的最大值(Emax),則所有斷面的腐朽程度(Ez)按下式計(jì)算:

    (5)

    腐朽程度(Ez)是一個(gè)百分?jǐn)?shù)(%),介于0~100之間。

    2 結(jié)果與分析

    根據(jù)上述算法得到木芯質(zhì)量損失率(作為腐朽程度真值Es)、電阻法腐朽程度(Ed)、應(yīng)力波法腐朽程度(Ey)及阻抗儀法腐朽程度(Ez)數(shù)據(jù),統(tǒng)計(jì)結(jié)果如表1所示。

    2.1 電阻檢測(cè)結(jié)果與腐朽程度真值之間的關(guān)系

    用最小二乘法對(duì)Ed和Es進(jìn)行一元線性回歸分析,得到回歸方程Ed=0.665 9Es+11.852(R2=0.516 6,R=0.718 8,P<0.01),表明Ed和Es之間呈顯著正相關(guān)性。若以Es=30%和Es=50%為界線,分成Es<30%、30%≤Es<50%及Es≥50% 3個(gè)區(qū)域進(jìn)行一元線性回歸分析(圖5),在Es<30%時(shí)得到相應(yīng)回歸方程Ed=1.303 3Es+4.285 5(R2=0.677 3,R=0.823 0,P<0.01); 在30%≤Es<50%時(shí)得到相應(yīng)回歸方程Ed=0.757 69Es+13.427(R2=0.067 7,R=0.260 2,P<0.01);在Es≥50%時(shí)得到相應(yīng)回歸方程Ed=1.362 2Es-30.292(R2=0.304 9,R=0.552 2,P<0.01)。以上表明,在Es的各個(gè)區(qū)域,Ed和Es之間均呈顯著正相關(guān)性,但Es<30%區(qū)域的相關(guān)系數(shù)高于30%≤Es<50%和Es≥50% 2個(gè)區(qū)域的相關(guān)系數(shù),因此,在Es<30%時(shí),電阻檢測(cè)法檢測(cè)結(jié)果更可靠。

    表1 各種檢測(cè)法測(cè)定結(jié)果統(tǒng)計(jì)

    圖5 電阻檢測(cè)結(jié)果Ed和腐朽程度真值Es之間的關(guān)系Fig.5 Relationship between result of measuring resistance Ed and the truth value of degree of decay Es

    2.2 應(yīng)力波檢測(cè)結(jié)果與腐朽程度真值之間的關(guān)系

    用最小二乘法對(duì)Ey和Es進(jìn)行一元線性回歸分析,得到回歸方程Ey=0.999 3Es+7.536 9(R2=0.637,R=0.799,P<0.01),此時(shí)Ey和Es之間呈顯著正相關(guān)性。同樣以Es=30%和Es=50%為界線,分成Es<30%、30%≤Es<50%及Es≥50% 3個(gè)區(qū)域進(jìn)行一元線性回歸分析(圖6),在Es<30%時(shí)得到相應(yīng)回歸方程Ey=1.250 1Es+5.949 7(R2=0.399,R=0.632,P<0.01); 在30%≤Es<50%時(shí)得到相應(yīng)回歸方程為Ey=1 085 7Es+1.268(R2=0.432,R=0.658,P<0.01); 在Es≥50%時(shí)得到相應(yīng)回歸方程Ey=2.011Es-47.8(R2=0.626,R=0.791,P<0.01)。以上表明,在Es的各個(gè)區(qū)域,Ey和Es之間均呈顯著正相關(guān)性,但在Es≥30%時(shí)的擬合度高于Es<30%時(shí)的擬合度,說明在腐朽程度大于等于30%時(shí),應(yīng)力波檢測(cè)結(jié)果能較清晰地表征木材腐朽程度。

    圖6 應(yīng)力波檢測(cè)結(jié)果Ey和腐朽程度真值Es之間的關(guān)系Fig.6 Relationship between result of stress wave Ey and the truth value of degree of decay Es

    2.3 阻抗儀檢測(cè)結(jié)果與腐朽程度真值之間的關(guān)系

    用最小二乘法對(duì)Ez和Es進(jìn)行一元線性回歸分析,得到相應(yīng)回歸方程Ez=1.08Es+5.317(R2=0.688,R=0.829,P<0.01),此時(shí)Ez和Es之間呈顯著正相關(guān)性,且擬合度較高。同樣以Es=30%和Es=50%為界線,分成Es<30%、30%≤Es<50%及Es≥50% 3個(gè)區(qū)域進(jìn)行一元線性回歸分析(圖7),在Es<30%時(shí)得到相應(yīng)回歸方程Ez=1.311Es+2.538(R2=0.506,R=0.711,P<0.01); 在30%≤Es<50%時(shí)得到相應(yīng)回歸方程Ez=1.184Es-2.398(R2=0.337,R=0.58,P<0.01); 在Es≥50%時(shí)得到相應(yīng)回歸方程Ez=4.177Es-182.63(R2=0.843,R=0.914,P<0.01)。以上表明,Ez和Es之間在Es≥50%時(shí)擬合度最高,此時(shí)阻抗儀檢測(cè)結(jié)果能更清晰地表征木材質(zhì)量損失率的情況。

    圖7 阻抗儀檢測(cè)結(jié)果Ez和腐朽程度真值Es之間的關(guān)系Fig.6 Relationship between result of resistance Ez and the truth value of degree of decay Es

    3 討論

    3.1 電阻檢測(cè)結(jié)果與腐朽程度真值之間的關(guān)系分析

    木材腐朽是指木材受到木腐菌感染,其細(xì)胞壁被分解引起的木材腐爛和解體。當(dāng)木材出現(xiàn)腐朽變色時(shí),菌絲生長(zhǎng)需要大量水分,將使腐朽變色區(qū)含水率升高;同時(shí),木腐菌分解木材的細(xì)胞壁,釋放出木材細(xì)胞中的離子。研究顯示,隨著立木變色腐朽的發(fā)生,其病變組織中的鉀、鈣、錳、鎂等金屬離子含量增加;而隨著陽(yáng)離子濃度增加,立木中腐朽變色的組織與健康組織相比,其電阻會(huì)明顯減小(Houston, 1971)。電阻法測(cè)得的腐朽程度Ed主要反映了腐朽區(qū)域在活立木體內(nèi)的含水率升高比例和其結(jié)構(gòu)中金屬離子含量增加情況,腐朽程度真值Es主要反映了木材質(zhì)量損失情況,質(zhì)量損失率與腐朽分布范圍、木材結(jié)構(gòu)破壞程度和力學(xué)強(qiáng)度密切相關(guān)(楊忠, 2005),所以2種方法都能反映立木腐朽程度。在本研究中,當(dāng)Es<30%時(shí),Ed和Es之間呈顯著相關(guān)性,且相關(guān)系數(shù)最高,說明電阻法在檢測(cè)木材初期腐朽方面具有較好的靈敏性,若能盡早地發(fā)現(xiàn)木材腐朽,就能盡早防治,與其他方法相比,這是電阻檢測(cè)法的優(yōu)勢(shì)。

    3.2 應(yīng)力波檢測(cè)結(jié)果與腐朽程度真值之間的關(guān)系分析

    木材中的大量纖維素、半纖維素和木質(zhì)素若被木腐菌腐蝕,則會(huì)發(fā)生腐朽,此時(shí)其密度相應(yīng)減小,木材內(nèi)部將形成空洞。而當(dāng)應(yīng)力波在有缺陷的木材中傳播時(shí),會(huì)繞過缺陷部位沿其邊緣傳播,傳播路徑由直線變?yōu)榍€,傳播時(shí)間增加,速度降低(徐華東等, 2014b)。應(yīng)力波法測(cè)得的腐朽程度Ey主要反映了立木內(nèi)部缺陷的大小情況,而木材的質(zhì)量損失率也與腐朽分布范圍、木材結(jié)構(gòu)破壞程度和力學(xué)強(qiáng)度密切相關(guān)(楊忠, 2005),所以2種方法都能反映立木腐朽程度,其間存在相關(guān)關(guān)系。本研究中,當(dāng)30%≤Es<50%時(shí),Ey和Es之間在3種方法中相關(guān)系數(shù)最大,說明腐朽程度在此區(qū)域時(shí),應(yīng)力波檢測(cè)結(jié)果更準(zhǔn)確。

    3.3 阻抗儀測(cè)定結(jié)果與腐朽程度真值之間的關(guān)系分析

    由于腐朽導(dǎo)致木材密度和強(qiáng)度下降(黃榮鳳等, 2007),因此阻抗儀法測(cè)得的腐朽程度Ez主要反映了腐朽區(qū)域在活立木體內(nèi)的分布比例和木材力學(xué)強(qiáng)度下降比例。而木芯質(zhì)量損失率計(jì)算法利用了腐朽導(dǎo)致木材顏色和密度發(fā)生變化的原理,根據(jù)木芯質(zhì)量損失率測(cè)得的腐朽程度Es主要反映了木材質(zhì)量損失情況,質(zhì)量損失率也與腐朽分布范圍、木材結(jié)構(gòu)破壞程度和力學(xué)強(qiáng)度密切相關(guān)。所以2種方法都從多個(gè)角度體現(xiàn)了腐朽程度,并且在理論上具有很強(qiáng)的相關(guān)性。在本研究中,在Es≥50%區(qū)域,Ez和Es之間的相關(guān)系數(shù)最大,說明在腐朽嚴(yán)重時(shí),阻抗儀檢測(cè)結(jié)果更準(zhǔn)確。

    3.4 3種無損檢測(cè)方法檢測(cè)結(jié)果的比較

    從整體上來看,電阻檢測(cè)和應(yīng)力波檢測(cè)結(jié)果與腐朽程度真值之間的擬合程度低于阻抗儀檢測(cè)結(jié)果與腐朽程度真值之間的擬合程度?;盍⒛倦娮柚凳芏喾N因素影響,如環(huán)境濕度、溫度、含水率、腐朽程度、生長(zhǎng)季節(jié)及測(cè)定部位等,且電阻檢測(cè)比較敏感,容易發(fā)生誤判(王立海等, 2001; Justetal., 1998),所以其擬合程度較低。同樣的,應(yīng)力波檢測(cè)結(jié)果受橫截面形狀、缺陷類型、木材腐朽程度及傳感器數(shù)量等因素的影響(徐華東等, 2011; 鮑震宇等, 2013a; 劉澤旭等, 2014),所以其擬合程度較低。對(duì)于阻抗儀,其探頭阻力與木材力學(xué)強(qiáng)度密切相關(guān)(安源等, 2008),不易受其他因素影響。而質(zhì)量損失率大小與木材絕干密度和木材力學(xué)強(qiáng)度直接相關(guān),故阻抗儀檢測(cè)結(jié)果和質(zhì)量損失率檢測(cè)真值之間的擬合程度較高。

    從腐朽不同階段來看,電阻檢測(cè)法在腐朽程度Es<30%時(shí),其檢測(cè)結(jié)果Ed和腐朽程度真值Es之間的擬合程度相對(duì)較高; 在30%≤Es<50%區(qū)域,應(yīng)力波檢測(cè)結(jié)果Ey和腐朽程度真值Es之間的擬合程度在3種方法中相對(duì)較高; 而在Es≥50%區(qū)域,則是阻抗儀檢測(cè)結(jié)果Ez與腐朽程度真值Es之間的擬合程度在3種方法中相對(duì)較高。這可能與木材的腐朽過程有關(guān),在木材腐朽初期,其質(zhì)量、外觀等特征變化微小,但化學(xué)成分發(fā)生了變化。池玉杰(2002)指出,木材腐朽菌可以通過菌絲或者根狀菌索的延伸和蔓延來繁殖,當(dāng)木材腐朽菌侵染進(jìn)入木材細(xì)胞并在木材細(xì)胞間定居后,便分泌多種酶,將木材細(xì)胞壁中的纖維素、半纖維素和木質(zhì)素分解為糖類,并進(jìn)一步將其作為養(yǎng)料消化分解。木材電阻主要與腐朽區(qū)域的含水率升高比例和其結(jié)構(gòu)中金屬離子含量的增加相關(guān),因此在腐朽初期,電阻檢測(cè)方法比較準(zhǔn)確。隨著木腐菌的分解過程開始穩(wěn)定,其電阻變化也就趨于平緩。 對(duì)于應(yīng)力波檢測(cè)方法,在腐朽初期,木材內(nèi)部并沒有形成空洞,因此應(yīng)力波檢測(cè)不那么準(zhǔn)確,而隨著腐朽程度的增加,導(dǎo)致木材內(nèi)部慢慢形成空洞(李堅(jiān), 2002)。有研究(徐華東等, 2014a; 2014b)指出,當(dāng)應(yīng)力波遇到空洞時(shí),不會(huì)沿直線傳播,而是繞著空洞缺陷沿空洞周圍進(jìn)行傳播,此時(shí)應(yīng)力波傳播路程增長(zhǎng),則相應(yīng)的應(yīng)力波傳播時(shí)間也增加。而在應(yīng)力波儀器計(jì)算過程中,應(yīng)力波傳播路程被設(shè)定為兩點(diǎn)之間最短的距離,應(yīng)力波傳播時(shí)間增加,使傳感器接受信號(hào)時(shí)間變長(zhǎng),故檢測(cè)得到的應(yīng)力波傳播速度變化能反映腐朽程度。對(duì)于阻抗儀檢測(cè)方法,無論腐朽程度如何,一旦發(fā)生腐朽均會(huì)使得木材密度發(fā)生變化,并導(dǎo)致力學(xué)性能下降,而阻抗儀探頭阻力與木材力學(xué)強(qiáng)度密切相關(guān),在腐朽程度更嚴(yán)重時(shí),阻力檢測(cè)值下降得更為明顯(黃榮鳳等, 2007; 張厚江等, 2011),本試驗(yàn)也表明在腐朽程度更大時(shí),阻抗儀檢測(cè)結(jié)果較為準(zhǔn)確。

    4 結(jié)論

    采用電阻斷層成像、應(yīng)力波及阻抗儀3種無損檢測(cè)方法對(duì)立木樹干內(nèi)部腐朽進(jìn)行了檢測(cè)和定量表征,結(jié)果表明:

    1) 電阻檢測(cè)結(jié)果(Ed)和腐朽程度真值(Es)之間呈顯著正相關(guān)性,且在腐朽程度Es<30%時(shí),相對(duì)其他2種方法,電阻測(cè)定的腐朽程度(Ed)和腐朽程度真值(Es)之間的相關(guān)系數(shù)最大(R=0.823,P<0.01)。

    2) 應(yīng)力波檢測(cè)結(jié)果(Ey)和腐朽程度真值(Es)之間呈顯著正相關(guān)性,在腐朽程度Es≥30%時(shí)的相關(guān)系數(shù)大于腐朽程度小于30%時(shí)的相關(guān)系數(shù)。3種方法對(duì)比來看,在30%≤Es<50%時(shí)應(yīng)力波檢測(cè)結(jié)果(Ey)和腐朽程度真值(Es)之間的相關(guān)系數(shù)最大(R=0.658,P<0.01),說明在此腐朽程度時(shí)應(yīng)力波檢測(cè)結(jié)果能更清晰地表征木材腐朽程度情況。

    3) 阻抗儀檢測(cè)結(jié)果(Ez)與腐朽程度真值(Es)之間呈顯著正相關(guān)性,且Es≥50%時(shí),在3種方法中其檢測(cè)結(jié)果(Ez)與腐朽程度真值(Es)之間的相關(guān)系數(shù)最大(R=0.914,P<0.01)。

    4) 電阻斷層成像、應(yīng)力波和阻抗儀3種無損檢測(cè)方法均能在不同腐朽程度上表征活立木木材質(zhì)量損失率,而質(zhì)量損失率是木材腐朽程度的主要表征指標(biāo)(池玉杰, 2002),因此3種無損檢測(cè)方法都能在一定范圍內(nèi)有效檢測(cè)活立木的腐朽程度。

    5) 電阻斷層成像、應(yīng)力波和阻抗儀3種無損檢測(cè)方法各有特點(diǎn),電阻法對(duì)早期腐朽(Es<30%)比較敏感, 而應(yīng)力波法在腐朽程度稍微嚴(yán)重(30%≤Es<50%)時(shí)檢測(cè)結(jié)果相對(duì)較為準(zhǔn)確, 阻抗儀法對(duì)不同腐朽程度的檢測(cè)結(jié)果都較為準(zhǔn)確,且在腐朽程度嚴(yán)重(Es≥50%)時(shí)最為準(zhǔn)確,但阻抗儀會(huì)對(duì)木材造成微損傷。因此,在實(shí)際腐朽檢測(cè)中應(yīng)當(dāng)根據(jù)實(shí)際檢測(cè)需要來選擇。

    安 源, 殷亞方, 姜笑梅, 等. 2008. 應(yīng)力波和阻抗儀技術(shù)勘查木結(jié)構(gòu)立柱腐朽分布. 建筑材料學(xué)報(bào), 11(4): 457-463.

    (An Y, Yin Y F, Jiang X M,etal. 2008. Inspection of decay distribution in wood column by stress wave and resistograph techniques. Journal of Building Materials, 11(4): 457-463. [in Chinese])

    鮑震宇, 邸向輝, 王立海, 等. 2013a. 小興安嶺天然林中冷杉立木含水率對(duì)應(yīng)力波傳播速度的影響. 林業(yè)科技, 38(1): 22-24.

    (Bao Z Y, Di X H, Wang L H,etal. 2013a. Effects of moisture content on propagation velocity of stress wave inAbiesfabricstanding trees in Xiaoxing’anling nature forest. Forestry Science & Technology, 38(1): 22-24. [in Chinese])

    鮑震宇, 王立海. 2013b. 電阻測(cè)試法在立木腐朽檢測(cè)中的應(yīng)用研究進(jìn)展.森林工程, 29(6): 47-51.

    (Bao Z Y, Wang L Y. 2013b. Application of electrical resistance testing method in detection of decay in standing trees.Forest Engineering, 29(6):47-51. [in Chinese])

    池玉杰. 2002. 木材腐朽與木材腐朽菌. 北京: 科學(xué)出版社.

    (Chi Y J. 2002. Wood decay and wood decay fungi. Beijing: Science Press. [in Chinese])

    高 珊, 王立海, 王 洋. 2013. 東北林區(qū)十個(gè)樹種凍結(jié)活立木中應(yīng)力波傳播速度試驗(yàn)研究. 森林工程, 29(4): 47-52.

    (Gao S, Wang L H, Wang Y. 2013. Experimental study on stress wave propagation velocity in frozen standing trees of 10 species in northeast China. Forest Engineering, 29(4): 47-52. [in Chinese])

    葛曉雯, 王立海, 孫天用, 等. 2014. 基于應(yīng)力波和阻抗儀技術(shù)的旱柳內(nèi)部腐朽定量檢測(cè).林業(yè)科技開發(fā),28(5): 87-91.

    (Ge X W, W L H, Sun T Y,etal. 2014. Quantitative detection ofSalixmatsudanainner decay based on stress wave and resistograph techniques. China Forestry Science and Technology, 28(5): 87-91. [in Chinese])

    黃榮鳳, 王曉歡, 李 華, 等. 2007. 古建筑木材內(nèi)部腐朽狀況阻力儀檢測(cè)結(jié)果的定量分析. 北京林業(yè)大學(xué)學(xué)報(bào), 29(6): 167-171.

    (Huang R F,Wang X H, Li H,etal. 2007. Quantitative analysis on the detected results by resisto-graph on inside wood decay of ancient architecture. Journal of Beijing Forestry University, 29(6): 167-171. [in Chinese])

    李 堅(jiān). 2002. 木材科學(xué). 北京: 高等教育出版社.

    (Li J. 2002. Wood science. Beijing: Higher Education Press.[in Chinese])

    李文彬. 2008. 第15屆國(guó)際木材無損檢測(cè)會(huì)議. 國(guó)際學(xué)術(shù)動(dòng)態(tài),(4): 41-42.

    (Li W B. 2008. 15thInternational symposium on nondestructive testing of wood. International Academic Developments, (4): 41-42. [in Chinese])

    林文樹,楊慧敏,王立海. 2005. 超聲波與應(yīng)力波在木材內(nèi)部缺陷檢測(cè)中的對(duì)比研究. 林業(yè)科技,30(2):39-41.

    (Lin W S, Yang H M, Wang L H. 2005. Comparative study on ultrasonic and stress wave for nondestructive test of wood defects.Forestry Science & Technology, 30(2):39-41. [in Chinese])

    劉澤旭, 邸向輝, 王立海, 等. 2014. 檢測(cè)角對(duì)健康立木中應(yīng)力波傳播速度的影響. 東北林業(yè)大學(xué)學(xué)報(bào), 42(4): 105-108.

    (Liu Z X, Di X H, Wang L H,etal. 2014. Effect of different detection angle on propagation velocity of stress wave in health standing trees.Journal of Northeast Forestry University, 42(4): 105-108. [in Chinese])

    孫天用, 王立海, 徐華東, 等. 2014. 小興安嶺土壤化學(xué)性質(zhì)對(duì)紅松活立木樹干腐朽的影響. 北京林業(yè)大學(xué)學(xué)報(bào), 36(2):30-37.

    (Sun T Y, Wang L H, Xu H D,etal. 2014. Effects of soil chemical properties on trunk decay of Korean pine standing trees in Xiaoxing’an Mountains of northeastern China. Journal of Beijing Forestry University, 36(2):30-37. [in Chinese])

    孫天用, 王立海, 孫墨瓏. 2013. 小興安嶺紅松活立木樹干腐朽與立地土壤理化特性的關(guān)系. 應(yīng)用生態(tài)學(xué)報(bào), 24(7): 1837-1842.

    (Sun T Y, Wang L H, Sun M L. 2013. Correlations between standing trees trunk decay degree and soil physical-chemical properties in Korean pine-broadleaved mixed forest in Xiaoxing’ an Mountains of northeast China. Chinese Journal of Applied Ecology, 24(7): 1837-1842. [in Chinese])

    王朝志,張厚江. 2006. 應(yīng)力波用于木材和活立木無損檢測(cè)的研究進(jìn)展. 林業(yè)機(jī)械與木工設(shè)備, 34(3): 9-13.

    (Wang C Z, Zhang H J. 2006. The research progress in wood and standing tree based on stress wave method. Forestry Machinery & Woodworking Equipment, 34(3): 9-13. [in Chinese])

    王立海, 楊學(xué)春, 徐凱宏. 2001. 木材無損檢測(cè)技術(shù)的研究現(xiàn)狀與進(jìn)展. 森林工程, 17(6): 1-3.

    (Wang L H, Yang X C, Xu K H. 2001. The present situation and advance in development of wood non-destructive testing techniques. Forest Engineering, 17(6): 1-3. [in Chinese])

    王立海,徐華東,閆在興, 等. 2008. 傳感器的數(shù)量與分布對(duì)應(yīng)力波檢測(cè)原木缺陷效果的影響. 林業(yè)科學(xué),44(55): 115-121.

    (Wang L H, Xu H D, Yan Z X,etal. 2008. Effects of sensor quantity and planar distribution on testing results of log defects based on stress wave.Scientia Silvae Sinicae, 44(5):115-121. [in Chinese])

    吳華橋, 周啟友, 劉漢樂. 2008. 樹干水分分布及運(yùn)移的高密度電阻率成像法時(shí)空監(jiān)測(cè). 地球物理學(xué)進(jìn)展, 23(4): 1310-1316.

    (Wu H Q, Zhou Q Y, Liu H L. 2008. Temporal and spatial monitoring of water distribution and migration in a tree trunk using electrical resistivity tomography. Progress in Geophysics, 23(4): 1310-1316. [in Chinese])

    徐華東, 王立海, 游祥飛, 等. 2010. 應(yīng)力波在旱柳立木內(nèi)的傳播規(guī)律分析及其安全評(píng)價(jià). 林業(yè)科學(xué), 46(8):145-150.

    (Xu H D, Wang L H, You X F. 2010. Analysis of stress wave propagation in hankow willow standing trees and stability assessment.Scientia Silvae Sinicae, 46(8):145-150. [in Chinese])

    徐華東, 王立海. 2014a. 空洞對(duì)木材應(yīng)力波傳播路徑的影響. 東北林業(yè)大學(xué)學(xué)報(bào), 42(4): 82-88.

    (Xu H D, Wang L H. 2014a. Effects of cavity on propagation path of stress wave in wood.Journal of Northeast Forestry University, 42(4): 82-88. [in Chinese])

    徐華東, 王立海. 2011.溫度和含水率對(duì)紅松木材中應(yīng)力波傳播速度的影響. 林業(yè)科學(xué), 47(9): 123-128.

    (Xu H D, Wang L H. 2011. Effects of moisture content and temperature on propagation velocity of stress waves in Korean pine wood.Scientia Silvae Sinicae, 47(9): 123-128. [in Chinese])

    徐華東, 徐國(guó)祺, 王立海, 等.2014b. 原木橫截面應(yīng)力波傳播時(shí)間等值線繪制及影響因素分析. 林業(yè)科學(xué),50(4): 95-100.

    (Xu H D, Xu G Q, Wang L H,etal. 2014b. Construction of stress wave time isolines on log cross section and analysis of its effect factors.Scientia Silvae Sinicae, 50(4):95-100. [in Chinese])

    徐凱宏. 2004. 壓阻式木材內(nèi)部缺陷類無損檢測(cè)理論與技術(shù)的研究. 哈爾濱: 東北林業(yè)大學(xué)博士學(xué)位論文

    (Xu K H. 2004. Theory and technology of detecting inner-defects in wood with a probing drill device. Harbin: PhD thesis of Northeast Forestry University. [in Chinese])

    徐 速, 周啟友, 劉漢樂. 2006. ERT 法在樹干水分運(yùn)移監(jiān)測(cè)中的運(yùn)用. 桂林工學(xué)院學(xué)報(bào), 26(3): 347-352.

    (Xu S, Zhou Q Y, Liu H L. 2006. Application of ERT method in monitoring tree water movement. Journal of Guilin University of Technology, 26(3): 347-352. [in Chinese])

    楊學(xué)春, 王立海. 2005. 應(yīng)力波在原木中傳播理論的研究. 林業(yè)科學(xué), 41(5): 132-138.

    (Yang X C, Wang L H. 2005. Study on the propagation theories of stress wave in log. Scientia Silvae Sinicae, 41(5): 132-138. [in Chinese])

    楊學(xué)春, 王立海. 2007. 原木內(nèi)部腐朽應(yīng)力波二維圖像的獲取與分析. 林業(yè)科學(xué), 43(11):93-97.

    (Yang X C, Wang L H. 2007. Gain and analysis of two-dimensional images of interior decay of logs with stress wave method.Scientia Silvae Sinicae, 43(11):93-97. [in Chinese])

    楊 忠. 2005. 近紅外光譜預(yù)測(cè)人工林濕地松木材性質(zhì)與腐朽特性的研究. 北京: 中國(guó)林業(yè)科學(xué)研究院博士學(xué)位論文.

    (Yang Z. 2005. Modeling slash pine plantation wood properties and decay using near infrared spectroscopy. Beijing: PhD thesis of Chinese Academy of Forestry. [in Chinese])

    張厚江, 朱 磊, 孫燕良, 等. 2011. 古建筑木構(gòu)件材料主要力學(xué)性能檢測(cè)方法研究. 北京林業(yè)大學(xué)學(xué)報(bào), 33(5): 126-129.

    (Zhang H J, Zhu L, Sun Y L,etal. 2011. Determining main mechanical properties of ancient architectural timber. Journal of Beijing Forestry University, 33(5): 126-129. [in Chinese])

    周啟友, 劉漢樂, 唐明杰. 2009. 電阻率成像法在樹干水分吸收過程研究中的應(yīng)用. 生態(tài)學(xué)雜志, 28(2): 350-356.

    (Zhou Q Y, Liu H L, Tang M J. 2009. Application of electrical resistivity tomography in studying water uptake process in tree trunk. Chinese Journal of Ecology, 28(2): 350-356. [in Chinese])

    Brazee N J, Marra R E, G?cke L,etal. 2011. Non-destructive assessment of internal decay in three hardwood species of northeastern North America using sonic and electrical impedance tomography. Forestry, 84(1): 33-39.

    Houston D R. 1971. Discoloration and decay in red maple and yellow birch. Forest Sci, 17(4):402-406.

    Just A, Jacbbs F.1998. Elektrische widerstandstomographie zur untersuchung des gesundheitszustandes von b?umen. Seminar, Hochaufl?sende Geoelektrik, Universit?t Leipzig, Institut für Geophysik und Geologie.

    Wang L H, Xu H D, Zhou C L,etal. 2007. Effect of sensor quantity on measurement accuracy of log inner defects by using stress wave. Journal of Forestry Research, 18(3):221-225, 251.

    Wang X P, Allison R B. 2008. Decay detection in red oak trees using a combination of visual inspection, acoustic testing, and resistance microdrilling. Arboriculture & Urban Forestry, 34(1): 1-4.

    Wang X P, Divos F, Pilon C,etal. 2004. Assessment of decay in standing timber using stress wave timing nondestructive evaluation tools. General Technical Report FPL-GTR-147. Madison, WI: Forest Products Laboratory, Forest Service, United States Department of Agriculture.

    (責(zé)任編輯 石紅青)

    Quantitative Detection of Internal Decay Degree for Standing Trees Based on Three NDT Methods—Electric Resistance Tomography,Stress Wave Imaging and Resistograph Techniques

    Yue Xiaoquan1, 2Wang Lihai1Wang Xinglong1Rong Binbin1Ge Xiaowen1Liu Zexu1Chen Qingyao2

    (1.CollegeofEngineeringandTechnology,NortheastForestryUniversityHarbin150040;2.CollegeofTransportationandCivilEngineering,FujianAgricultureandForestryUniversityFuzhou350002)

    【Objective】 Electric resistance tomography, stress wave imaging and resistograph were employed to detect and quantitatively characterize the internal decay for standing trees,then these three NDT(nondestructive testing)techniques were made comparisons in order to find out the appropriate technique matched with specific condition of standing trees in the forest field.【Method】An investigation about the accuracy of different NDT methods was carried out in the experiment forest of Northeast Forestry University in Harbin, Heilongjiang Province. 100 cross-sections ofFraxinusmandshuricaandPopulus×beijingensisstanding trees were tested by four methods: electric resistance tomography, stress wave tomography, resistograph and estimation weight loss ratios of wooden increment cores. Taking wood core samples as research objects, testing results of three other kinds of NDT methods were compared.Esdetermined by estimating weight loss ratios of wooden cores was regarded as the truth value of decay extent. Using ordinary least square regression to analyze the relationship betweenEsandEd(degree of decay determined by electric resistance tomography),EsandEy(degree of decay determined by stress wave),EsandEz(degree of decay determined by resistograph). 【Result】Results showed that three NDT methods were able to estimate the different degree of decay forFraxinusmandshuricaandPopulus×beijingensisstanding trees; As a whole, the fitting degree ofEdandEs,EyandEswere both lower than that ofEzandEs. In different decay degree, whenEs<30%,Edhad a strong positive correlation withEs(R=0.823,P<0.01), while when 30%≤Es<50%,Eyhave a significant positive correlation relationship withEs(R=0.658,P<0.01), electric resistance tomography had more sensitively reflect in incipient decay of standing trees, and stress wave had more accurate diagnosis on detecting middle decay of standing tress.Ezhad a strong positive correlation relationship withEs, and whenEs≥50%,Ezhave a more positive correlation relationship withEs(R=0.914,P<0.01).【Conclusion】 Electric resistance tomography showed better diagnosis than the other two methods for incipient decay of standing trees, while stress wave imaging method used in the middle stage of decay testing was best, and resistograph can be used in the different stages of decay. It was suggested that each technique could be employed in the practical internal decay testing for standing trees according to the decay stage and operational conditions.

    mass loss rate; electric resistance tomography; velocity of stress wave; wood core resistance loss; decay extent

    10.11707/j.1001-7488.20170315

    2015-12-01;

    2016-02-23。

    948項(xiàng)目“便攜式立木腐朽電阻斷層成像關(guān)鍵技術(shù)引進(jìn)”(2014-4-78); 福建省教育廳科技項(xiàng)目(JA15155); 福建省高水平大學(xué)建設(shè)專項(xiàng)基金(61201401801); 國(guó)家自然科學(xué)基金項(xiàng)目(31300474);黑龍江省自然科學(xué)基金面上項(xiàng)目(C201410)。

    S715.3;S791.247

    A

    1001-7488(2017)03-0138-09

    *王立海為通訊作者。

    国产精品三级大全| 视频在线观看一区二区三区| 熟妇人妻不卡中文字幕| 综合色丁香网| 香蕉国产在线看| 亚洲在久久综合| 国产无遮挡羞羞视频在线观看| 成人影院久久| 国产精品偷伦视频观看了| 久久久久精品人妻al黑| 亚洲av在线观看美女高潮| 最近手机中文字幕大全| 免费久久久久久久精品成人欧美视频| 精品国产一区二区三区久久久樱花| 一区二区三区激情视频| 一二三四中文在线观看免费高清| 日本wwww免费看| 黄色一级大片看看| 丝瓜视频免费看黄片| 最近最新中文字幕免费大全7| 欧美日韩综合久久久久久| 天天躁夜夜躁狠狠躁躁| 亚洲欧美清纯卡通| 久久热在线av| 国产精品秋霞免费鲁丝片| 国产精品一区二区在线观看99| 国产又爽黄色视频| 国产日韩欧美视频二区| 搡老乐熟女国产| 免费观看a级毛片全部| 99精国产麻豆久久婷婷| 久久久久国产一级毛片高清牌| 超碰成人久久| 国产精品一国产av| 亚洲男人天堂网一区| 老司机在亚洲福利影院| 最新在线观看一区二区三区 | 高清黄色对白视频在线免费看| 久久久久久人人人人人| 欧美国产精品一级二级三级| 黑丝袜美女国产一区| 亚洲精品国产av蜜桃| 人人妻人人添人人爽欧美一区卜| 夫妻性生交免费视频一级片| av不卡在线播放| 欧美精品一区二区免费开放| 美女午夜性视频免费| 婷婷色综合www| 99久久精品国产亚洲精品| 人人妻人人澡人人爽人人夜夜| 亚洲成人一二三区av| 国产亚洲最大av| 欧美精品人与动牲交sv欧美| 一区二区日韩欧美中文字幕| 黄色视频在线播放观看不卡| 久久久久精品国产欧美久久久 | 宅男免费午夜| 99九九在线精品视频| 波多野结衣一区麻豆| 另类亚洲欧美激情| 国产精品熟女久久久久浪| 99久久精品国产亚洲精品| 天堂中文最新版在线下载| 久久狼人影院| 美女国产高潮福利片在线看| 一级毛片 在线播放| 日日撸夜夜添| 天天添夜夜摸| av线在线观看网站| 欧美乱码精品一区二区三区| 国产成人91sexporn| 久久国产精品大桥未久av| 国产日韩欧美亚洲二区| 热re99久久国产66热| 一本—道久久a久久精品蜜桃钙片| 久久人妻熟女aⅴ| 91成人精品电影| 女性被躁到高潮视频| 高清av免费在线| 亚洲av欧美aⅴ国产| 一个人免费看片子| 热99国产精品久久久久久7| 777米奇影视久久| 1024视频免费在线观看| 十八禁网站网址无遮挡| 日韩视频在线欧美| 捣出白浆h1v1| 捣出白浆h1v1| 午夜福利在线免费观看网站| www.av在线官网国产| 久久国产精品大桥未久av| 人体艺术视频欧美日本| 99国产精品免费福利视频| 亚洲第一av免费看| www.自偷自拍.com| 免费黄色在线免费观看| 亚洲自偷自拍图片 自拍| 久久99热这里只频精品6学生| 久久久久视频综合| 操出白浆在线播放| 日韩av在线免费看完整版不卡| 日韩,欧美,国产一区二区三区| 搡老岳熟女国产| 午夜精品国产一区二区电影| 在线免费观看不下载黄p国产| 美女福利国产在线| 大码成人一级视频| 日本色播在线视频| 亚洲精品视频女| 成人亚洲欧美一区二区av| 国产精品麻豆人妻色哟哟久久| 丝袜脚勾引网站| 日本av手机在线免费观看| 中文精品一卡2卡3卡4更新| av福利片在线| 18在线观看网站| 黑丝袜美女国产一区| 欧美激情 高清一区二区三区| tube8黄色片| 亚洲专区中文字幕在线 | 亚洲久久久国产精品| 国产精品麻豆人妻色哟哟久久| 丝瓜视频免费看黄片| 亚洲成人一二三区av| 精品少妇一区二区三区视频日本电影 | 丝袜喷水一区| 亚洲熟女毛片儿| 韩国高清视频一区二区三区| 亚洲精品日本国产第一区| 国产在视频线精品| 人妻人人澡人人爽人人| a 毛片基地| 麻豆精品久久久久久蜜桃| 亚洲专区中文字幕在线 | 51午夜福利影视在线观看| 亚洲一码二码三码区别大吗| 免费高清在线观看视频在线观看| 视频区图区小说| 99香蕉大伊视频| 99热国产这里只有精品6| 亚洲伊人久久精品综合| 操出白浆在线播放| 操出白浆在线播放| 久久久久视频综合| 国产精品欧美亚洲77777| 久久久久精品国产欧美久久久 | 国语对白做爰xxxⅹ性视频网站| 婷婷色综合www| 老司机影院成人| 日韩中文字幕欧美一区二区 | 男人爽女人下面视频在线观看| 午夜福利视频精品| 大片免费播放器 马上看| 美女福利国产在线| 国产高清国产精品国产三级| 成人国产av品久久久| 成年动漫av网址| 在线观看人妻少妇| 欧美成人精品欧美一级黄| 国产精品国产三级国产专区5o| 国产精品99久久99久久久不卡 | 九草在线视频观看| 久久av网站| 熟女av电影| 蜜桃在线观看..| 国产伦理片在线播放av一区| 在线观看人妻少妇| 国产精品蜜桃在线观看| 国产精品久久久久久精品古装| 精品国产超薄肉色丝袜足j| 建设人人有责人人尽责人人享有的| 亚洲欧美一区二区三区国产| 下体分泌物呈黄色| 国产老妇伦熟女老妇高清| 男人爽女人下面视频在线观看| 无限看片的www在线观看| 美国免费a级毛片| 亚洲国产最新在线播放| 毛片一级片免费看久久久久| 亚洲国产欧美网| 久久精品国产综合久久久| 热re99久久国产66热| 精品国产国语对白av| 色播在线永久视频| 午夜日韩欧美国产| 国产免费视频播放在线视频| 中文字幕av电影在线播放| 嫩草影院入口| 人妻人人澡人人爽人人| 久久韩国三级中文字幕| 久久97久久精品| 午夜福利视频精品| 99九九在线精品视频| 亚洲自偷自拍图片 自拍| 亚洲精品国产av成人精品| 久久韩国三级中文字幕| av国产久精品久网站免费入址| 99久久99久久久精品蜜桃| 美国免费a级毛片| 欧美人与性动交α欧美软件| 久久久久网色| 欧美日韩精品网址| 97人妻天天添夜夜摸| 国产精品一区二区精品视频观看| 韩国av在线不卡| 丰满饥渴人妻一区二区三| 日韩中文字幕欧美一区二区 | 美国免费a级毛片| 侵犯人妻中文字幕一二三四区| 高清不卡的av网站| 美女扒开内裤让男人捅视频| 欧美精品av麻豆av| 久久青草综合色| 亚洲国产欧美网| 免费在线观看视频国产中文字幕亚洲 | 少妇 在线观看| 最近2019中文字幕mv第一页| 最近2019中文字幕mv第一页| 国产精品99久久99久久久不卡 | 亚洲欧美精品自产自拍| 久久久久久久久久久久大奶| 1024视频免费在线观看| 日本色播在线视频| 日韩一卡2卡3卡4卡2021年| 精品午夜福利在线看| 国产高清国产精品国产三级| 日本爱情动作片www.在线观看| 午夜福利视频精品| av电影中文网址| 老司机影院成人| 韩国av在线不卡| 午夜精品国产一区二区电影| 国产亚洲av片在线观看秒播厂| 在线观看免费视频网站a站| 人妻一区二区av| 波野结衣二区三区在线| 永久免费av网站大全| 亚洲av国产av综合av卡| 1024视频免费在线观看| 青草久久国产| 一区二区av电影网| 国产男人的电影天堂91| 国产日韩一区二区三区精品不卡| 日本爱情动作片www.在线观看| 丝袜美足系列| 9热在线视频观看99| 丰满饥渴人妻一区二区三| 国产深夜福利视频在线观看| 精品少妇黑人巨大在线播放| 亚洲国产欧美一区二区综合| 欧美在线一区亚洲| 欧美日韩视频精品一区| 亚洲欧美成人综合另类久久久| 色综合欧美亚洲国产小说| 国产人伦9x9x在线观看| 男男h啪啪无遮挡| 黄色视频在线播放观看不卡| 午夜激情av网站| 亚洲国产成人一精品久久久| 免费少妇av软件| 99re6热这里在线精品视频| 国产高清不卡午夜福利| 国产精品 国内视频| 观看av在线不卡| 一本—道久久a久久精品蜜桃钙片| 国产精品 欧美亚洲| 久久97久久精品| 曰老女人黄片| 纯流量卡能插随身wifi吗| 欧美97在线视频| 在线精品无人区一区二区三| 精品国产乱码久久久久久男人| 久久99精品国语久久久| 操美女的视频在线观看| 久久影院123| av不卡在线播放| 激情视频va一区二区三区| 波多野结衣av一区二区av| 人妻一区二区av| 香蕉丝袜av| 亚洲精品久久久久久婷婷小说| 菩萨蛮人人尽说江南好唐韦庄| 国产高清国产精品国产三级| 中文字幕精品免费在线观看视频| 精品国产超薄肉色丝袜足j| 国产成人系列免费观看| 99国产精品免费福利视频| 制服丝袜香蕉在线| 国产 一区精品| 日韩av免费高清视频| 久久国产亚洲av麻豆专区| 国产精品免费视频内射| 新久久久久国产一级毛片| 久久天堂一区二区三区四区| 天堂中文最新版在线下载| 观看av在线不卡| 欧美另类一区| 欧美日韩国产mv在线观看视频| 亚洲三区欧美一区| 在线观看一区二区三区激情| 久久鲁丝午夜福利片| 母亲3免费完整高清在线观看| 亚洲美女视频黄频| 伦理电影免费视频| 高清欧美精品videossex| 女人被躁到高潮嗷嗷叫费观| 成人午夜精彩视频在线观看| 最近中文字幕2019免费版| 免费人妻精品一区二区三区视频| 中文字幕色久视频| 少妇猛男粗大的猛烈进出视频| 一级a爱视频在线免费观看| videos熟女内射| 深夜精品福利| 满18在线观看网站| 国产成人91sexporn| 综合色丁香网| 九色亚洲精品在线播放| 一个人免费看片子| 国产成人精品在线电影| bbb黄色大片| 久久人人爽av亚洲精品天堂| 国产老妇伦熟女老妇高清| 乱人伦中国视频| 91精品三级在线观看| 亚洲精品自拍成人| 可以免费在线观看a视频的电影网站 | 国产精品.久久久| 午夜日本视频在线| 大码成人一级视频| 好男人视频免费观看在线| 精品国产超薄肉色丝袜足j| 久久国产精品男人的天堂亚洲| 午夜福利在线免费观看网站| 亚洲色图 男人天堂 中文字幕| 多毛熟女@视频| 超碰成人久久| 在线亚洲精品国产二区图片欧美| 尾随美女入室| 久久天躁狠狠躁夜夜2o2o | 日韩精品有码人妻一区| 欧美日韩视频高清一区二区三区二| 亚洲欧美日韩另类电影网站| 美女中出高潮动态图| 肉色欧美久久久久久久蜜桃| 亚洲国产av新网站| 国产成人免费无遮挡视频| 亚洲av综合色区一区| 午夜免费男女啪啪视频观看| 亚洲精品aⅴ在线观看| 日韩大片免费观看网站| 国产1区2区3区精品| 亚洲精品在线美女| 精品国产一区二区三区四区第35| 亚洲av电影在线进入| 日韩精品有码人妻一区| 观看美女的网站| 亚洲欧美激情在线| 赤兔流量卡办理| 一级,二级,三级黄色视频| 2018国产大陆天天弄谢| 五月天丁香电影| 少妇精品久久久久久久| 香蕉丝袜av| 七月丁香在线播放| 少妇人妻精品综合一区二区| a级片在线免费高清观看视频| 欧美日韩一区二区视频在线观看视频在线| 免费在线观看黄色视频的| 人人妻人人添人人爽欧美一区卜| 在线观看www视频免费| 国产一卡二卡三卡精品 | 日本色播在线视频| 国产日韩一区二区三区精品不卡| 国产乱来视频区| 国产午夜精品一二区理论片| av网站免费在线观看视频| 国产男人的电影天堂91| 99久国产av精品国产电影| 无遮挡黄片免费观看| 国产精品av久久久久免费| 国产一区二区在线观看av| 嫩草影院入口| 乱人伦中国视频| 老司机影院毛片| 极品人妻少妇av视频| 国产精品熟女久久久久浪| 黄片小视频在线播放| 亚洲精品aⅴ在线观看| 极品人妻少妇av视频| 一区二区三区四区激情视频| 国产乱来视频区| avwww免费| 国产一区二区激情短视频 | 黄色怎么调成土黄色| 亚洲国产最新在线播放| 18禁观看日本| 亚洲国产日韩一区二区| 亚洲精品,欧美精品| av一本久久久久| 免费在线观看视频国产中文字幕亚洲 | 精品第一国产精品| 女性被躁到高潮视频| 欧美国产精品va在线观看不卡| 亚洲av成人不卡在线观看播放网 | 亚洲一卡2卡3卡4卡5卡精品中文| 国产片内射在线| 在线免费观看不下载黄p国产| 人妻 亚洲 视频| 男女下面插进去视频免费观看| 精品久久久久久电影网| 欧美日韩视频高清一区二区三区二| 国产精品一区二区精品视频观看| 国产xxxxx性猛交| 一级毛片黄色毛片免费观看视频| 欧美日韩福利视频一区二区| 亚洲av电影在线观看一区二区三区| 久久久久国产精品人妻一区二区| 欧美97在线视频| 男女之事视频高清在线观看 | 热99久久久久精品小说推荐| 久久久久久久久久久久大奶| 国产亚洲午夜精品一区二区久久| 五月天丁香电影| 美女中出高潮动态图| 黄频高清免费视频| xxx大片免费视频| 男人爽女人下面视频在线观看| 久久国产精品男人的天堂亚洲| 亚洲七黄色美女视频| 大码成人一级视频| 欧美 亚洲 国产 日韩一| 国产一区有黄有色的免费视频| 精品国产乱码久久久久久男人| 少妇的丰满在线观看| 中文字幕精品免费在线观看视频| 亚洲少妇的诱惑av| 免费观看人在逋| 精品午夜福利在线看| 成年人午夜在线观看视频| 亚洲国产毛片av蜜桃av| 在线观看免费午夜福利视频| 精品少妇久久久久久888优播| 美女视频免费永久观看网站| 欧美激情极品国产一区二区三区| 亚洲天堂av无毛| 亚洲精品av麻豆狂野| 久久人人97超碰香蕉20202| 久久鲁丝午夜福利片| 国产国语露脸激情在线看| 美女福利国产在线| 十八禁高潮呻吟视频| 欧美日韩综合久久久久久| 亚洲伊人久久精品综合| 日韩欧美精品免费久久| 国产亚洲av片在线观看秒播厂| 欧美老熟妇乱子伦牲交| 日韩,欧美,国产一区二区三区| 日本猛色少妇xxxxx猛交久久| 国产成人欧美| 高清黄色对白视频在线免费看| 成人漫画全彩无遮挡| 久久精品国产a三级三级三级| 成人毛片60女人毛片免费| 天堂8中文在线网| 精品亚洲乱码少妇综合久久| 交换朋友夫妻互换小说| 免费观看性生交大片5| 欧美av亚洲av综合av国产av | 亚洲婷婷狠狠爱综合网| 成人午夜精彩视频在线观看| 国产成人免费观看mmmm| 97人妻天天添夜夜摸| 性高湖久久久久久久久免费观看| 一本色道久久久久久精品综合| 一级黄片播放器| videosex国产| 欧美中文综合在线视频| 女的被弄到高潮叫床怎么办| 成年美女黄网站色视频大全免费| 欧美久久黑人一区二区| 极品少妇高潮喷水抽搐| 欧美日韩成人在线一区二区| 校园人妻丝袜中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| av在线观看视频网站免费| 美国免费a级毛片| 久久av网站| 18禁裸乳无遮挡动漫免费视频| 天堂俺去俺来也www色官网| 你懂的网址亚洲精品在线观看| 制服人妻中文乱码| 伊人久久国产一区二区| 婷婷色麻豆天堂久久| 丝袜美足系列| 亚洲成人国产一区在线观看 | 久久天躁狠狠躁夜夜2o2o | 黑人猛操日本美女一级片| 久久久久久人妻| 又大又爽又粗| 黄色视频不卡| 成年女人毛片免费观看观看9 | videos熟女内射| 18在线观看网站| 午夜福利影视在线免费观看| 欧美激情 高清一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品成人av观看孕妇| 国产熟女欧美一区二区| 啦啦啦视频在线资源免费观看| 亚洲成人国产一区在线观看 | 日本wwww免费看| 亚洲欧美日韩另类电影网站| 欧美最新免费一区二区三区| 9色porny在线观看| 午夜福利一区二区在线看| 中文字幕色久视频| 亚洲av在线观看美女高潮| 别揉我奶头~嗯~啊~动态视频 | 国产精品久久久久久人妻精品电影 | 欧美激情高清一区二区三区 | 欧美激情 高清一区二区三区| 男女边吃奶边做爰视频| 搡老岳熟女国产| 精品一区二区免费观看| 国产精品久久久久久精品古装| 亚洲av中文av极速乱| a级毛片黄视频| 黄频高清免费视频| 亚洲综合精品二区| 亚洲av日韩精品久久久久久密 | 国产一区二区三区av在线| 亚洲国产精品999| 一本—道久久a久久精品蜜桃钙片| 九九爱精品视频在线观看| 亚洲色图 男人天堂 中文字幕| 97精品久久久久久久久久精品| 免费观看av网站的网址| 美女福利国产在线| 18禁国产床啪视频网站| 狠狠精品人妻久久久久久综合| 老汉色av国产亚洲站长工具| 欧美久久黑人一区二区| 国产精品国产av在线观看| 亚洲精品国产一区二区精华液| 精品人妻一区二区三区麻豆| 少妇被粗大猛烈的视频| 欧美亚洲日本最大视频资源| 大香蕉久久成人网| 亚洲精品成人av观看孕妇| av.在线天堂| 亚洲国产欧美网| 国产在线免费精品| 91老司机精品| 亚洲国产精品一区三区| 一区福利在线观看| 免费不卡黄色视频| 999久久久国产精品视频| 国产精品成人在线| 免费黄色在线免费观看| 国产成人精品久久二区二区91 | 亚洲成人一二三区av| 亚洲激情五月婷婷啪啪| www日本在线高清视频| 叶爱在线成人免费视频播放| 在线观看免费高清a一片| 亚洲久久久国产精品| 精品一区在线观看国产| h视频一区二区三区| 成人国产av品久久久| 9色porny在线观看| 少妇被粗大的猛进出69影院| 国产精品免费大片| 一个人免费看片子| 日日撸夜夜添| av在线app专区| 久久这里只有精品19| 老汉色∧v一级毛片| 欧美精品一区二区免费开放| 亚洲精品第二区| 亚洲一区二区三区欧美精品| 午夜激情av网站| 欧美精品一区二区大全| 亚洲精品日本国产第一区| 免费观看av网站的网址| 国语对白做爰xxxⅹ性视频网站| 亚洲av日韩精品久久久久久密 | 国产女主播在线喷水免费视频网站| 中文字幕精品免费在线观看视频| 人体艺术视频欧美日本| 好男人视频免费观看在线| 国产成人精品在线电影| 国产成人a∨麻豆精品| 国产在线免费精品| 亚洲,一卡二卡三卡| 欧美激情极品国产一区二区三区| 看免费av毛片| www.av在线官网国产| 一二三四中文在线观看免费高清| 九草在线视频观看| 国产成人免费无遮挡视频| 亚洲中文av在线| 久久国产亚洲av麻豆专区| 大香蕉久久网| 99久久99久久久精品蜜桃| 黄色一级大片看看| 黄色怎么调成土黄色| 午夜福利,免费看| 午夜免费鲁丝| 国产在线视频一区二区| av国产精品久久久久影院| 九草在线视频观看| 国产1区2区3区精品| 午夜福利免费观看在线| 亚洲国产最新在线播放|