高级检索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于阻抗仪和ERT技术的红松内部腐朽定量检测

时小龙 王立海 徐华东 徐庆波 曹时凯 许明贤

时小龙, 王立海, 徐华东, 徐庆波, 曹时凯, 许明贤. 基于阻抗仪和ERT技术的红松内部腐朽定量检测[J]. 北京林业大学学报, 2017, 39(9): 102-111. doi: 10.13332/j.1000-1522.20170196
引用本文: 时小龙, 王立海, 徐华东, 徐庆波, 曹时凯, 许明贤. 基于阻抗仪和ERT技术的红松内部腐朽定量检测[J]. 北京林业大学学报, 2017, 39(9): 102-111. doi: 10.13332/j.1000-1522.20170196
SHI Xiao-long, WANG Li-hai, XU Hua-dong, XU Qing-bo, CAO Shi-kai, XU Ming-xian. Quantitative characterization of decay detection in standing trees of Korean pine based on Resistograph and ERT[J]. Journal of Beijing Forestry University, 2017, 39(9): 102-111. doi: 10.13332/j.1000-1522.20170196
Citation: SHI Xiao-long, WANG Li-hai, XU Hua-dong, XU Qing-bo, CAO Shi-kai, XU Ming-xian. Quantitative characterization of decay detection in standing trees of Korean pine based on Resistograph and ERT[J]. Journal of Beijing Forestry University, 2017, 39(9): 102-111. doi: 10.13332/j.1000-1522.20170196

基于阻抗仪和ERT技术的红松内部腐朽定量检测

doi: 10.13332/j.1000-1522.20170196
基金项目: 

国家自然科学基金项目 31570547

详细信息
    作者简介:

    时小龙,博士生。主要研究方向:木材无损检测。Email:1105318962@qq.com  地址:150040 黑龙江省哈尔滨市和兴路26号东北林业大学工程技术学院

    通讯作者:

    王立海,教授,博士生导师。主要研究方向:无损检测。Email:lihaiwang@yahoo.com  地址:同上

  • 中图分类号: S781.5

Quantitative characterization of decay detection in standing trees of Korean pine based on Resistograph and ERT

  • 摘要: 为了定量表征红松活立木内部腐朽情况,在黑龙江省五营国家森林公园选取了15棵有不同腐朽程度的红松,应用Resistograph树木微钻阻抗仪和PiCUS Treetronic树木电阻断层成像仪得到南北向与东西向的微钻阻力变化和同样高度上的截面电阻率分布图,对比分析阻抗仪测量方向上的电阻率变化和微钻阻力变化。用树木生长锥在样木的测量位置附近各取南北向与东西向的木芯,将其带回实验室测量质量损失率。将微钻阻力的下降幅度定义为阻力损失,对阻力损失和木芯质量损失率进行回归分析。还用腐朽和健康部分的电阻率均值定义了电阻率的异常程度, 并对其和阻力损失进行了回归分析。以阻抗仪为基准在电阻率分布图上对腐朽部位定位,应用监督分类的方法在图上将腐朽与健康部分区分出来,并计算腐朽面积所占比例。然后用最小二乘法将算得的比例与木芯质量损失率作回归分析。结果表明:平均阻力损失与木芯质量损失率之间存在显著的线性相关关系,阻抗仪可以较准确地表征一维直线上的腐朽情况;微钻阻力变化曲线上发生阻力值明显下降的位置对应的电阻率也会发生明显变化,阻力损失值和电阻率异常程度存在显著指数关系,还能推断出电阻断层成像检测对于轻度腐朽更加敏感,而阻抗仪检测对较严重的腐朽更加敏感;联合阻抗仪和电阻断层成像技术能够确定腐朽的位置和面积,得到的腐朽面积所占比例与木芯质量损失率有显著线性相关关系。
  • 图  1  木芯照片

    Figure  1.  Photo of wood core

    图  2  电阻率和微钻阻力随探头钻入深度的变化

    H1H2分别为波谷的左右两边在右侧微钻阻力坐标轴上的投影长度(resi);D1为波谷在横坐标上的投影长度(mm)。

    Figure  2.  Electrical resistivity and micro-drilling resistance varying with drilling depth

    H1 and H2 are projection lengths (resi) of left and right sides of the trough on the vertical axis of mocro-drilling resistance curve, respectively; D1 is projection length (mm) of the trough on the horizontal axis.

    图  3  质量损失率与平均阻力损失散点图

    Figure  3.  Scatter-plot of mass loss rate and average of resistance loss

    图  4  ERT图像

    Figure  4.  ERT image

    图  5  电阻率异常程度和阻力损失散点图

    Figure  5.  Scatter-plot of abnormal extent of electrical resistivity and resistance loss

    图  6  ERT图像分析过程

    Figure  6.  Process of analyzing ERT image

    图  7  质量损失率和ERT图像中腐朽面积比例散点图

    Figure  7.  Scatter-plot of mass loss rate and percentage of decay area in ERT image

    表  1  样木的基本情况

    Table  1.   Profile of 15 sample trees to be measured

    样木编号
    No. of sample tree
    胸径
    DBH/mm
    树高
    Height of tree/m
    含水率
    Moisture content/%
    141122.840.88
    279631.438.44
    361829.750.46
    460027.969.92
    571330.846.97
    648727.741.50
    746125.653.67
    835921.754.42
    943924.934.33
    1064030.349.17
    1156327.532.01
    1258328.434.13
    1353726.856.97
    1466529.235.39
    1546826.153.90
    下载: 导出CSV

    表  2  质量损失率与平均阻力损失

    Table  2.   Mass loss rate and average of resistance loss

    样木编号
    No. of sample tree
    平均阻力损失
    Average of resistance loss/resi
    质量损失率
    Mass loss rate/%
    119.0032.48
    221.5146.24
    37.2419.80
    45.0436.12
    52.1510.46
    622.5535.37
    74.4524.63
    84.6419.33
    926.0040.04
    106.377.81
    1132.1851.15
    1219.6443.47
    134.2223.95
    1415.4429.45
    1526.7468.92
    下载: 导出CSV

    表  3  阻力损失及对应的电阻率异常程度

    Table  3.   Resistance loss and corresponding abnormal extent of electrical resistivity

    样木编号
    No. of sample tree
    方向
    Direction
    波谷个数
    Number of troughs
    左端圆心距
    Center distance of left end (d1)/mm
    右端圆心距
    Center distance of right end (d2)/mm
    阻力损失
    Resistance loss/resi
    电阻率异常程度
    Abnormal extent of electrical resistivity/%
    1东西向East-west378.031.59.14103.83
    31.5-44.08.5043.73
    -62.5-77.01.3739.63
    2南北向North-south2266.5-43.527.36116.56
    257.0185.56.6857.53
    东西向East-west2185.5167.01.0242.73
    142.0-23.07.9477.21
    3东西向East-west2299.0250.04.769.72
    167.0132.02.4852.31
    4东西向East-west1160.0119.05.0496.23
    5南北向North-south2-10.5-30.51.4929.56
    -30.5-55.51.6623.60
    东西向East-west1314.0284.51.145.17
    6南北向North-south4107.558.55.0833.98
    51.013.55.5814.15
    13.50.51.2217.42
    -21.5-81.56.4821.20
    东西向East-west2156.571.07.1742.33
    39.5-83.519.56121.21
    7南北向North-south2-52.0-72.52.758.12
    -77.5-92.51.3719.72
    东西向East-west28.0-19.53.0042.00
    -108.5-125.51.7761.47
    8南北向North-south1127.0106.54.640.75
    9南北向North-south229.5-28.010.81111.24
    -28.0-155.520.69127.96
    东西向East-west4113.554.57.7390.79
    44.5-8.03.8063.61
    -8.0-74.57.8338.58
    -116.0-128.51.144.37
    10南北向North-south225.0-23.03.1562.00
    -34.0-76.03.9475.07
    东西向East-west1194.0106.05.6464.32
    11南北向North-south164.5-154.534.62129.35
    东西向East-west3160.5141.51.5121.56
    141.598.53.4010.99
    98.5-134.524.83111.51
    12南北向North-south2216.5151.56.836.83
    96.5-33.54.464.46
    东西向East-west1136.5-103.527.99112.42
    13南北向North-south1-121.5-152.03.4143.99
    东西向East-west1-121.5-166.55.0391.07
    14南北向North-south4238.5170.56.0123.65
    170.5155.50.7953.95
    140.5106.52.7536.76
    106.522.54.8926.78
    东西向East-west2161.55.013.06109.34
    -15.5-71.53.373.49
    15南北向North-south1188.0-167.032.03114.87
    东西向East-west2185.0140.05.1240.39
    140.0-56.016.33120.13
    下载: 导出CSV

    表  4  ERT图像中腐朽面积占截面比例

    Table  4.   Percentage of decay area in ERT image

    样木编号No. of sample tree123456789101112131415
    腐朽面积比例Percentage of decay area (Pd)/%33.6040.8828.9034.0926.8937.6523.2421.4242.5319.5955.2165.2522.3843.8175.37
    下载: 导出CSV
  • [1] 魏旸艳, 陆文明, 郎书平.中国木材进口影响因素分析[J].世界林业研究, 2009, 22(2): 78-80. http://www.sjlyyj.com/ch/reader/view_abstract.aspx?file_no=20090216&flag=1

    WEI Y Y, LU W M, LANG S P. An analysis of main influencing factors for timber import of China[J]. World Forestry Research, 2009, 22(2): 78-80. http://www.sjlyyj.com/ch/reader/view_abstract.aspx?file_no=20090216&flag=1
    [2] 王立海, 杨学春, 徐凯宏.木材无损检测技术的研究现状与进展[J].森林工程, 2001, 17(6): 1-3. doi:  10.3969/j.issn.1001-005X.2001.06.001

    WANG L H, YANG X C, XU K H. The research situations and development of non-destructive testing for wood properties[J]. Forest Engineering, 2001, 17(6): 1-3. doi:  10.3969/j.issn.1001-005X.2001.06.001
    [3] BRIANK B, VOICHITA B, FERENC D, et al. Nondestructive testing and evaluation of wood: a worldwide research update[J]. Forest Products Journal, 2009, 59(3): 7-14.
    [4] DIVOS F. Stress wave based tomography for trees evaluation[C]//Proceedings of the 12th International Symposium on Nondestructive Testing of Wood. 2000: 13-15.
    [5] WANG X, ALLISON R B. Decay detection in red oak trees using a combination of visual inspection, acoustic testing, and resistance microdrilling[J]. Arboriculture & Urban Forestry, 2008, 34(1): 1-4. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=b23f4a37293271e9180ada36d3e65048
    [6] 孙天用, 王立海.基于应力波与X射线二维CT图像原木内部腐朽无损检测[J].森林工程, 2011, 27(6): 26-29. doi:  10.3969/j.issn.1001-005X.2011.06.006

    SUN T Y, WANG L H. Non-destructive testing of log internal decay based on two-dimensional CT images of stress wave and X-ray testing[J]. Forest Engineering, 2011, 27(6): 26-29. doi:  10.3969/j.issn.1001-005X.2011.06.006
    [7] 葛晓雯, 王立海, 孙天用, 等.基于应力波和阻抗仪技术的旱柳内部腐朽定量检测[J].林业科技开发, 2014, 28(5): 87-91. http://d.old.wanfangdata.com.cn/Periodical/lykjkf201405023

    GE X W, WANG L H, SUN T Y, et al. Quantitative detection of Salix matsudana inner decay based on stress wave and resistograph techniques[J]. China Forestry Science and Technology, 2014, 28(5): 87-91. http://d.old.wanfangdata.com.cn/Periodical/lykjkf201405023
    [8] 安源, 殷亚方, 姜笑梅, 等.应力波和阻抗仪技术勘查木结构立柱腐朽分布[J].建筑材料学报, 2008, 11(4): 457-463. doi:  10.3969/j.issn.1007-9629.2008.04.015

    AN Y, YIN Y F, JIANG X M, et al. Inspection of decay distribution in wood column by stress wave and resistograph techniques[J]. Journal of Building Materials, 2008, 11(4): 457-463. doi:  10.3969/j.issn.1007-9629.2008.04.015
    [9] 鲍震宇, 王立海.电阻测试法在立木腐朽检测中的应用研究进展[J].森林工程, 2013, 29(6): 47-51. doi:  10.3969/j.issn.1001-005X.2013.06.012

    BAO Z Y, WANG L H. Application of electrical resistance testing method in detection of decay in standing trees[J]. Forest Engineering, 2013, 29(6): 47-51. doi:  10.3969/j.issn.1001-005X.2013.06.012
    [10] 孙天用, 王立海, 孙墨珑.小兴安岭红松活立木树干腐朽与立地土壤理化特性的关系[J].应用生态学报, 2013, 24(7): 1837-1842. http://d.old.wanfangdata.com.cn/Periodical/yystxb201307008

    SUN T Y, WANG L H, SUN M L. Correlations between standing trees trunk decay degree and soil physical-chemical properties in Korean pine-broadleaved mixed forest in Xiao Xing'an Mountains of northeast China[J]. Chinese Journal of Applied Ecology, 2013, 24(7): 1837-1842. http://d.old.wanfangdata.com.cn/Periodical/yystxb201307008
    [11] 李坚.木材科学研究[M].北京:科学出版社, 2009.

    LI J. Wood science research[M]. Beijing: Science Press, 2009.
    [12] 杨忠, 江泽慧, 任海青, 等.木材加速腐朽试验方法的研究[J].木材工业, 2007, 21(4): 12-14. doi:  10.3969/j.issn.1001-8654.2007.04.004

    YANG Z, JIANG Z H, REN H Q, et al. An accelerated laboratory test method to assess wood decay[J]. China Wood Industry, 2007, 21(4): 12-14. doi:  10.3969/j.issn.1001-8654.2007.04.004
    [13] 陶新民, 曲艺卓, 温少阳.基于电阻断层成像技术的立木探伤检测系统[J].应用科技, 2017, 44: 1-7. http://d.old.wanfangdata.com.cn/Periodical/yykj201704008

    TAO X M, QU Y Z, WEN S Y. Tree flaw detection system based on electrical resistivity tomography[J]. Applied Science and Technology, 2017, 44: 1-7. http://d.old.wanfangdata.com.cn/Periodical/yykj201704008
    [14] 鲍震宇.腐朽对木材电阻的影响研究[D].哈尔滨: 东北林业大学, 2015. http://cdmd.cnki.com.cn/Article/CDMD-10225-1015665582.htm

    BAO Z Y. The effect of decay on wood electrical resistance[D]. Harbin: Northeast Forestry University, 2015. http://cdmd.cnki.com.cn/Article/CDMD-10225-1015665582.htm
    [15] GUYOT A, OSTERGAARD K T, LENKOPANE M, et al. Using electrical resistivity tomography to differentiate sapwood from heartwood: application to conifers[J]. Tree Physiology, 2013, 33(2): 187-194. doi:  10.1093/treephys/tps128
    [16] 岳小泉, 王立海, 王兴龙, 等.电阻断层成像、应力波及阻抗仪3种无损检测方法对活立木腐朽程度的定量检测[J].林业科学, 2017, 53(3): 138-146. http://d.old.wanfangdata.com.cn/Periodical/lykx201703015

    YUE X Q, WANG L H, WANG X L, et al. Quantitative detection of internal decay degree for standing trees based on three NDT methods—electric resistance tomography, stress wave imaging and resistograph techniques[J]. Scientia Silvae Sinicae, 2017, 53(3): 138-146. http://d.old.wanfangdata.com.cn/Periodical/lykx201703015
    [17] SCHWARZE F W M R. Wood decay under the microscope[J]. Fungal Biology Reviews, 2007, 21(4): 133-170. doi:  10.1016/j.fbr.2007.09.001
    [18] 张新玲, 王立海, 高珊.小兴安岭4树种健康材和腐朽材结晶度的试验[J].东北林业大学学报, 2010, 38(7): 83-84. doi:  10.3969/j.issn.1000-5382.2010.07.026

    ZHANG X L, WANG L H, GAO X. Crystallinity of healthy and rotten wood of several tree species in Xiao Xing'an Mountain[J]. Journal of Northeast Forestry University, 2010, 38(7): 83-84. doi:  10.3969/j.issn.1000-5382.2010.07.026
    [19] WILCOX W W. Review of literature on the effects of early stages of decay on wood strength[J]. Wood & Fiber, 1978, 9(4): 252-257.
    [20] 尹思慈.木材学[M].北京:中国林业出版社, 1996.

    YIN S C. Wood science[M]. Beijing: China Forestry Publishing House, 1996.
    [21] 张富文, 许清风, 张治宇, 等.钻入阻抗法检测木材缺陷[J].无损检测, 2016, 38(1): 6-9. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=wsjc201601003

    ZHANG F W, XU Q F, ZHANG Z Y, et al. Wood defect inspection by drilling resistance method[J]. Nondestructive Testing, 2016, 38(1): 6-9. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=wsjc201601003
    [22] TATTAR T A, SHIGO A L, CHASE T. Relationship between the degree of resistance to a pulsed electric current and wood in progressive stages of discoloration and decay in living trees[J]. Canadian Journal of Forest Research, 1972, 2(3): 236-243. doi:  10.1139/x72-039
    [23] 蒋艳凰.遥感图像高精度并行监督分类技术研究[D].北京: 国防科学技术大学, 2004. http://cdmd.cnki.com.cn/Article/CDMD-90002-2005014429.htm

    JIANG Y H. High-accuracy and parallel algorithms for supervised remote-sensing image classification[D]. Beijing: National University of Defense Technology, 2004. http://cdmd.cnki.com.cn/Article/CDMD-90002-2005014429.htm
  • [1] 李美玉, 张守红, 王云琦, 谢朝帅, 李睿贤.  不同清理方式对北京市透水砖铺装渗透率衰减过程影响 . 北京林业大学学报, 2020, 42(3): 143-150. doi: 10.12171/j.1000-1522.20190062
    [2] 靳晓东, 姜立春.  基于树干不同形率的樟子松立木材积方程研建 . 北京林业大学学报, 2020, 42(3): 78-86. doi: 10.12171/j.1000-1522.20190047
    [3] 张静, 漆楚生, 母军.  热处理温度和时间对杉木质量损失和静曲强度的影响 . 北京林业大学学报, 2020, 42(10): 137-144. doi: 10.12171/j.1000-1522.20200257
    [4] 杨坤, 赵艳玲, 张建勇, 陈超, 赵鹏鹏.  利用无人机高分辨率影像进行树木高度提取 . 北京林业大学学报, 2017, 39(8): 17-23. doi: 10.13332/j.1000-1522.20160428
    [5] 时小龙, 王立海, 徐华东, 徐庆波, 曹时凯, 许明贤.  基于阻抗仪和ERT技术的红松内部腐朽定量检测 . 北京林业大学学报, 2017, 39(9): 102-111. doi: 10.13332/j.1000-1522.20170123
    [6] 费学海, 张一平, 宋清海, 李培广, 刘运通, 金艳强, 李敬, 
    林友兴, 张树斌.  元江干热河谷太阳辐射各分量及反照率变化特征 . 北京林业大学学报, 2016, 38(3): 1-10. doi: 10.13332/j.1000-1522.20150360
    [7] 赵荣玮, 张建军, 陈宝强, 王雅琼, 孙若修, 冯换成, 茹豪.  不同覆盖措施对土壤水分和当年造林成活率的影响 . 北京林业大学学报, 2016, 38(6): 87-93. doi: 10.13332/j.1000-1522.20150511
    [8] 解潍嘉, 黄侃, 李瑞平, 孙浩, 扈晶晶, 黄华国.  应用高分辨率卫星数据估算阔叶红松林乔木多样性 . 北京林业大学学报, 2015, 37(3): 20-26. doi: 10.13332/j.1000-1522.20140306
    [9] 王玉婷, 徐华东, 王立海, 李凤日, 孙虎.  小兴安岭天然林红松活立木腐朽率的调查研究 . 北京林业大学学报, 2015, 37(8): 97-104. doi: 10.13332/j.1000-1522.20150015
    [10] 刘昊, 高建民.  含水率和密度对木材应力波传播速度的影响 . 北京林业大学学报, 2014, 36(6): 154-158. doi: 10.13332/j.cnki.jbfu.2014.06.002
    [11] 孙锋, 周永东, 李晓玲, 吕建雄, 韩晨静, 詹满军.  桉木单板出材率的研究 . 北京林业大学学报, 2013, 35(4): 128-133.
    [12] 李超, 张明辉, 于建芳.  利用核磁共振自由感应衰减曲线测定木材含水率 . 北京林业大学学报, 2012, 34(4): 142-145.
    [13] 张厚江, 朱磊, 孙燕良, 王喜平.  古建筑木构件材料主要力学性能检测方法研究 . 北京林业大学学报, 2011, 33(5): 126-129.
    [14] 李贤军, 蔡智勇, 傅峰, 周永东.  间歇微波干燥过程中木材内含水率动态分布规律 . 北京林业大学学报, 2010, 32(4): 227-231.
    [15] 宋桂龙, 韩烈保, 张训忠, 辜再元.  土壤机械阻力对草地早熟禾根系生长的影响 . 北京林业大学学报, 2009, 31(1): 60-65.
    [16] 池玉杰, 志娟.  10种平伏形木材腐朽菌的培养特性 . 北京林业大学学报, 2008, 30(1): 85-92.
    [17] 周繇, 张运春, 闫德千, 谭健晖, 赵铁珍, 邢韶华, 张秀新, 王超, 尹增芳, 任云卯, 王春梅, 陈圆, 林勇明, 张玉兰, 李云开, 孙阁, 李春义, 张仁军, 王戈, 李昌晓, 梁善庆, 王蕾, 崔丽娟, 吴淑芳, 金莹杉, 张颖, 林娅, 刘杏娥, 周海宾, 余养伦, 马钦彦, 周荣伍, 王莲英, 张曼胤, 马履一, 王以红, 徐秋芳, 杨远芬, 洪滔, 吴普特, 高岚, 张桥英, 刘艳红, 黄华国, 翟明普, 罗建举, 于俊林, 樊汝汶, 张志强, 张明, 刘国经, 刘青林, 温亚利, 江泽慧, 赵勃, 钟章成, 江泽慧, 杨培岭, 安玉涛, 冯浩, 杨海军, 崔国发, 殷际松, 周国模, 罗鹏, 于文吉, 王小青, 田英杰, 邵彬, 陈学政, 刘俊昌, 汪晓峰, 柯水发, 张晓丽, 何春光, 吴承祯, 王希群, 费本华, 张本刚, 王玉涛, 周国逸, 蔡玲, 费本华, 刘爱青, 徐克学, 赵景刚, 徐昕, 骆有庆, 洪伟, 何松云, 魏晓华, 李敏, 王九中, 高贤明, 任树梅, 康峰峰, 马润国, 邬奇峰, 温亚利, 胡喜生, 林斌, 赵焕勋, 朱高浦, 安树杰, 赵弟行, 郑万建, 吴家森, 田平, 吴宁, 任海青, 宋萍, 卢俊峰, 李永祥, 范海兰.  美国白蛾入侵部分非经济损失评估案例分析 . 北京林业大学学报, 2007, 29(3): 171-177.
    [18] 韦艳葵, 党文杰, 汪杭军1, 方升佐, 段文霞, 李义良, 刘剑锋, 耿玉清, 李生宇, HUALi_zhong, 李雪华, 任强, 周传艳, 李国雷, 王立海, 吴丽娟, 薛康, 王兰珍, 刘鹏举, 朱小龙, 王旭, 赵铁珍, 张冬梅, 周宇飞, 朱波, 周国逸, 杨慧敏, 尹光彩, 杨娅, 黎明, 雷加强, 李建章, JIANGXi_dian, 宋永明, 韩士杰, 何茜, 余新晓, 李振基, 刘勇, 苏晓华, 高岚, 周亮, 刘勇, 方陆明, 阎秀峰, 崔同林, 鹿振友, 虞木奎, 周国逸, 唐小明, 徐扬, 王春林, HEXiu_bin, 王清文, 沈熙环, 玲, 程云清, 赖志华, 张冰玉, 柯水发, 李吉跃, 孙向阳, 宗文君, 喻理飞, 徐新文, 王新杰, 刘锐, 宋爱琴, 3, 李俊清, 齐涛, 孙阁, 张志毅, 王伟宏, 温亚利, 郭蓓, 李晓兰, 陈峻崎, , 茹广欣, 陈实, 国庆, 陈培金, 李丙文, 周晓梅, 张可栋, 刘志明, 长山, 周玉平, 蒋德明, 王晓静, 姚永刚, 王旭, 唐旭利, 王建林, 宋湛谦, 赵双荣, 关少华, 王春林, 陈放, 杨伟伟, 闫俊华, 郑凌峰.  美国白蛾入侵损失评估指标体系的构建 . 北京林业大学学报, 2007, 29(2): 156-160.
    [19] 张志山, 张春晓, 焦雯珺, 高克昌, 雷妮娅, 许景伟, 李俊, 金则新, 刘足根, 吴家兵, 邵杰, 马玲, 张建军, 武林, 周艳萍, 于海霞, 奚如春, 郎璞玫, 吕文华, 陆平, 周睿, 索安宁, 宋先亮, 郑景明, 李黎, 孙志蓉, 于文吉, 陈勇, 赵秀海, 葛剑平, 于志明, 张小由, 李钧敏, 朱清科, 纳磊, 翟明普, 马履一, 关德新, 赵广杰, 李传荣, 郑红娟, 戴伟, 赵文喆, 毕华兴, 余养伦, 习宝田, 蔡锡安, 韦方强, 朱教君, Kwei-NamLaw, 陈少良, 饶兴权, 盖颖, 杨永福, 朱艳燕, 李笑吟, 贾桂霞, 夏良放, 李增鸿, 王文全, 崔鹏, 方家强, 于波, 马履一, 张弥, 樊敏, 曾小平, 袁小兰, 谭会娟, 王天明, 李俊清, ClaudeDaneault, 王瑞刚, 江泽慧, 张宇清, 张春雨, 赵平, 李丽萍, 李庆卫, 韩士杰, 刘丽娟, 陈雪梅, 袁飞, 何明珠, 殷宁, 郭孟霞, 王卫东, 邓宗付, 吴秀芹, 王贺新, 贺润平, 张欣荣, 唐晓军, 刘鑫, 吴记贵, 王月海, 江杰, 孔俊杰, 郑敬刚, 王娜, 熊颖, 王旭琴, 于贵瑞, 蒋湘宁, 毛志宏, 葛剑平, 王贵霞, 王瑞辉, 林靓靓, 聂立水, 孙晓敏, 李新荣, 郭超颖, 董治良.  晋西黄土区坡面糙率的研究 . 北京林业大学学报, 2007, 29(1): 108-113.
    [20] 颜绍馗, 邓小文, 郑小贤, 杨平, 黄荣凤, 毛俊娟, 李瑞, 魏潇潇, 吴彩燕, 张洪江, 张莉俊, 刘杏娥, 王费新, 殷亚方, 胡万良, 何亚平, 秦爱光, 胡胜华, 高黎, 王芳, 白岗栓, 周永学, 袁怀文, 张璧光, 赵天忠, 杜社妮, 李猛, NagaoHirofumi, 樊军锋, 乔建平, 汪思龙, 王兆印, 戴思兰, 王晓欢, 王小青, 崔赛华, 罗晓芳, 张克斌, 谭学仁, 王胜华, 刘燕, 王正, 张岩, 孙向阳, 常旭, 费世民, 徐嘉, 王海燕, , 龚月桦, 李昀, 李华, 江玉林, KatoHideo, 张旭, 张占雄, 刘云芳, 陈放, 范冰, 孔祥文, 江泽慧, 韩士杰, 张双保, 高荣孚, , 郭树花, 陈秀明, 杨培华, 丁磊, 李媛良, 侯喜录, 陈宗伟, 常亮, 李晓峰, 任海青, IdoHirofumi, 刘秀英, , 张代贵, 费本华, 蒋俊明, 李考学, 张桂兰, 高建社, 薛岩, , 陈学平, 徐庆祥, 涂代伦, , 王晓东, 金鑫, 李雪峰, 续九如, 刘永红, , 张红丽, , 丁国权.  古建筑木材内部腐朽状况阻力仪检测结果的定量分析 . 北京林业大学学报, 2007, 29(6): 167-171.
  • 加载中
图(7) / 表 (4)
计量
  • 文章访问数:  2300
  • HTML全文浏览量:  114
  • PDF下载量:  5
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-06-02
  • 修回日期:  2017-08-29
  • 刊出日期:  2017-09-01

基于阻抗仪和ERT技术的红松内部腐朽定量检测

doi: 10.13332/j.1000-1522.20170196
    基金项目:

    国家自然科学基金项目 31570547

    作者简介:

    时小龙,博士生。主要研究方向:木材无损检测。Email:1105318962@qq.com  地址:150040 黑龙江省哈尔滨市和兴路26号东北林业大学工程技术学院

    通讯作者: 王立海,教授,博士生导师。主要研究方向:无损检测。Email:lihaiwang@yahoo.com  地址:同上
  • 中图分类号: S781.5

摘要: 为了定量表征红松活立木内部腐朽情况,在黑龙江省五营国家森林公园选取了15棵有不同腐朽程度的红松,应用Resistograph树木微钻阻抗仪和PiCUS Treetronic树木电阻断层成像仪得到南北向与东西向的微钻阻力变化和同样高度上的截面电阻率分布图,对比分析阻抗仪测量方向上的电阻率变化和微钻阻力变化。用树木生长锥在样木的测量位置附近各取南北向与东西向的木芯,将其带回实验室测量质量损失率。将微钻阻力的下降幅度定义为阻力损失,对阻力损失和木芯质量损失率进行回归分析。还用腐朽和健康部分的电阻率均值定义了电阻率的异常程度, 并对其和阻力损失进行了回归分析。以阻抗仪为基准在电阻率分布图上对腐朽部位定位,应用监督分类的方法在图上将腐朽与健康部分区分出来,并计算腐朽面积所占比例。然后用最小二乘法将算得的比例与木芯质量损失率作回归分析。结果表明:平均阻力损失与木芯质量损失率之间存在显著的线性相关关系,阻抗仪可以较准确地表征一维直线上的腐朽情况;微钻阻力变化曲线上发生阻力值明显下降的位置对应的电阻率也会发生明显变化,阻力损失值和电阻率异常程度存在显著指数关系,还能推断出电阻断层成像检测对于轻度腐朽更加敏感,而阻抗仪检测对较严重的腐朽更加敏感;联合阻抗仪和电阻断层成像技术能够确定腐朽的位置和面积,得到的腐朽面积所占比例与木芯质量损失率有显著线性相关关系。

English Abstract

时小龙, 王立海, 徐华东, 徐庆波, 曹时凯, 许明贤. 基于阻抗仪和ERT技术的红松内部腐朽定量检测[J]. 北京林业大学学报, 2017, 39(9): 102-111. doi: 10.13332/j.1000-1522.20170196
引用本文: 时小龙, 王立海, 徐华东, 徐庆波, 曹时凯, 许明贤. 基于阻抗仪和ERT技术的红松内部腐朽定量检测[J]. 北京林业大学学报, 2017, 39(9): 102-111. doi: 10.13332/j.1000-1522.20170196
SHI Xiao-long, WANG Li-hai, XU Hua-dong, XU Qing-bo, CAO Shi-kai, XU Ming-xian. Quantitative characterization of decay detection in standing trees of Korean pine based on Resistograph and ERT[J]. Journal of Beijing Forestry University, 2017, 39(9): 102-111. doi: 10.13332/j.1000-1522.20170196
Citation: SHI Xiao-long, WANG Li-hai, XU Hua-dong, XU Qing-bo, CAO Shi-kai, XU Ming-xian. Quantitative characterization of decay detection in standing trees of Korean pine based on Resistograph and ERT[J]. Journal of Beijing Forestry University, 2017, 39(9): 102-111. doi: 10.13332/j.1000-1522.20170196
  • 我国木材的综合利用率低,平均70%以下,与发达国家的平均80%以上有较大差距[1]。准确地检测木材缺陷对于提高木材利用率十分有效。木材无损检测能够在不破坏木材特性前提下,通过各种技术手段在较短时间内获得其理化特征来反映缺陷程度,是现今木材检测的重点研究方向[2]

    木材无损检测技术出现的相对较晚,但是近几十年中发展迅速,先后有学者将超声波、X射线、微波和电阻测量等技术应用到木材无损检测[3]。经过不断地探索和研究发现,每种方法都有各自的优点和不足。因此联合不同的测量方法对木材进行综合性评估有利于提高检测的精度,是未来木材无损检测的趋势[2]。国内外有学者采用联合方法对木材进行无损检测。Divos[4]使用声学断层成像和应力波技术,获得了树木内部截面上的缺陷形状。Wang等[5]比较了目测法、声学法及阻抗仪测量的特点,认为综合多种方法能够提高木材检测的准确性。孙天用和王立海[6]使用应力波和X射线二维CT图像技术测量了原木内部腐朽,发现两者都能反映腐朽区域的形状,且X射线能进一步反映腐朽程度。葛晓雯等[7]运用应力波和阻抗仪技术检测了旱柳立木的内部腐朽,结果表明两种方法都能确定腐朽的位置,微钻阻力和应力波速有一定相关性。

    应用阻抗仪和电阻断层成像仪(ERT)在野外环境下对内部有不同腐朽程度的红松(Pinus koraiensis)活立木进行测量。阻抗仪能够记录探针钻入木材时一维直线方向上的微钻阻力变化,一旦木材内部发生腐朽、裂纹、虫蛀等缺陷,其力学性能会有所下降,在微钻阻力变化曲线上相应地表现为阻力值的下降[8]。电阻断层成像技术是电阻法应用于木材无损检测的发展,基于电磁理论和有限元理论,结合计算机成像技术,能够反映木材的二维或三维的电阻率分布情况[9]。关于联合这两种检测手段进行木材无损检测的研究很少。本研究比较和评价这两种检测方法的测量结果,并根据各自的优点联合两种方法,明确划定红松活立木二维截面上的腐朽区域,为定量表征红松活立木腐朽提供理论依据。

    • 在黑龙江省五营国家森林公园内选择红松种植密集的区域,划定一块样地,其面积约为30 hm2,坡度平缓,一般不超过10°。由经验丰富的林场工作人员目测挑选15棵有不同腐朽程度的红松作为待测样木,相邻两棵样木的间距都大于10 m。15棵树的基本情况列于表 1

      表 1  样木的基本情况

      Table 1.  Profile of 15 sample trees to be measured

      样木编号
      No. of sample tree
      胸径
      DBH/mm
      树高
      Height of tree/m
      含水率
      Moisture content/%
      141122.840.88
      279631.438.44
      361829.750.46
      460027.969.92
      571330.846.97
      648727.741.50
      746125.653.67
      835921.754.42
      943924.934.33
      1064030.349.17
      1156327.532.01
      1258328.434.13
      1353726.856.97
      1466529.235.39
      1546826.153.90
    • 首先在距离地面130 cm处用胸径尺测量样木的直径, 记为D2,并使用指南针记录方向,同时用树木测高仪(瑞典制造)测量树高。再用PiCUS Treetronic树木电阻断层成像仪(德国制造)在同样高度的截面上测量,ERT的1号电压传感器置于正北(或南)方向,逆时针等间距布置12个传感器进行测量,测量完成后保存样木截面的电阻率图像。然后在同一高度上分别沿南北向和东西向用Resistograph树木微钻阻抗仪(德国制造)测量,记录微钻阻力变化曲线,文中所有阻力皆为微钻阻力。阻抗仪测量完成后,用树木生长锥(瑞典制造)分别在测量位置附近取1~2根木芯,放入密封袋中保存。图 1为6号样木的木芯照片。在实验室环境下取出木芯,用电子分析天平称其质量m1(g),再用电热鼓风干燥箱(70 ℃)将其烘干至恒重并记录质量m2(g)和长度a(a取树木生长锥的长度25 cm)。再选择木芯的一段健康部分测量其质量m3(g)和长度b(cm),用m3/b表示健康材的线密度,用am3/b表示健康材的质量,记为m4(g)。含水率MC(%)根据下式计算:

      图  1  木芯照片

      Figure 1.  Photo of wood core

      $$ \text{MC}=\frac{{{m}_{1}}-{{m}_{2}}}{{{m}_{2}}}\times 100 \% $$

      质量损失率L(%)由下面的公式计算得到。对同一棵树取得的多个木芯,以其均值代表它的质量损失率。

      $$ \mathit{L}=\frac{{{m}_{4}}-{{m}_{2}}}{{{m}_{4}}}\times 100 \% $$
    • 图 2实线部分为沿东西方向测得的6号样木的微钻阻力变化曲线图。发生腐朽等缺陷区域的力学强度会减小,阻抗仪测得的微钻阻力值相对于健康部分有所下降。从图上能够看出两处明显的微钻阻力下降区域,在横坐标上分别用线段EFGH标记其位置。用微钻阻力的下降幅度定量表征腐朽程度,定义阻力损失Z (resi)如下[10]

      图  2  电阻率和微钻阻力随探头钻入深度的变化

      Figure 2.  Electrical resistivity and micro-drilling resistance varying with drilling depth

      $$ Z=\frac{({{H}_{1}}+{{H}_{2}}){{D}_{1}}}{2{{D}_{2}}} $$

      式中:H1H2分别为波谷的左右两边在右侧微钻阻力坐标轴上的投影长度,resi;D1为波谷在横坐标上的投影长度,mm;D2为所测截面的直径,mm。表 2中的平均阻力损失由该样木东西方向上的总阻力损失与南北方向上总阻力损失之和除以2得到,15棵样木的不同测量方向上的阻力损失列于表 3

      表 2  质量损失率与平均阻力损失

      Table 2.  Mass loss rate and average of resistance loss

      样木编号
      No. of sample tree
      平均阻力损失
      Average of resistance loss/resi
      质量损失率
      Mass loss rate/%
      119.0032.48
      221.5146.24
      37.2419.80
      45.0436.12
      52.1510.46
      622.5535.37
      74.4524.63
      84.6419.33
      926.0040.04
      106.377.81
      1132.1851.15
      1219.6443.47
      134.2223.95
      1415.4429.45
      1526.7468.92

      表 3  阻力损失及对应的电阻率异常程度

      Table 3.  Resistance loss and corresponding abnormal extent of electrical resistivity

      样木编号
      No. of sample tree
      方向
      Direction
      波谷个数
      Number of troughs
      左端圆心距
      Center distance of left end (d1)/mm
      右端圆心距
      Center distance of right end (d2)/mm
      阻力损失
      Resistance loss/resi
      电阻率异常程度
      Abnormal extent of electrical resistivity/%
      1东西向East-west378.031.59.14103.83
      31.5-44.08.5043.73
      -62.5-77.01.3739.63
      2南北向North-south2266.5-43.527.36116.56
      257.0185.56.6857.53
      东西向East-west2185.5167.01.0242.73
      142.0-23.07.9477.21
      3东西向East-west2299.0250.04.769.72
      167.0132.02.4852.31
      4东西向East-west1160.0119.05.0496.23
      5南北向North-south2-10.5-30.51.4929.56
      -30.5-55.51.6623.60
      东西向East-west1314.0284.51.145.17
      6南北向North-south4107.558.55.0833.98
      51.013.55.5814.15
      13.50.51.2217.42
      -21.5-81.56.4821.20
      东西向East-west2156.571.07.1742.33
      39.5-83.519.56121.21
      7南北向North-south2-52.0-72.52.758.12
      -77.5-92.51.3719.72
      东西向East-west28.0-19.53.0042.00
      -108.5-125.51.7761.47
      8南北向North-south1127.0106.54.640.75
      9南北向North-south229.5-28.010.81111.24
      -28.0-155.520.69127.96
      东西向East-west4113.554.57.7390.79
      44.5-8.03.8063.61
      -8.0-74.57.8338.58
      -116.0-128.51.144.37
      10南北向North-south225.0-23.03.1562.00
      -34.0-76.03.9475.07
      东西向East-west1194.0106.05.6464.32
      11南北向North-south164.5-154.534.62129.35
      东西向East-west3160.5141.51.5121.56
      141.598.53.4010.99
      98.5-134.524.83111.51
      12南北向North-south2216.5151.56.836.83
      96.5-33.54.464.46
      东西向East-west1136.5-103.527.99112.42
      13南北向North-south1-121.5-152.03.4143.99
      东西向East-west1-121.5-166.55.0391.07
      14南北向North-south4238.5170.56.0123.65
      170.5155.50.7953.95
      140.5106.52.7536.76
      106.522.54.8926.78
      东西向East-west2161.55.013.06109.34
      -15.5-71.53.373.49
      15南北向North-south1188.0-167.032.03114.87
      东西向East-west2185.0140.05.1240.39
      140.0-56.016.33120.13

      质量损失率是反映木材性质的一个重要指标[11],在本试验中用它作为腐朽程度的真值。15棵样木的平均阻力损失和质量损失率列于表 2。以质量损失率为自变量,平均阻力损失为因变量进行一元回归分析,结果如图 3所示。回归方程Z=0.512L-2.232(R2=0.646,P<0.01)的拟合程度较高,表明了木芯的质量损失率越大,对应的阻力损失(即微钻阻力的下降程度)越大,两者呈线性正相关关系。

      图  3  质量损失率与平均阻力损失散点图

      Figure 3.  Scatter-plot of mass loss rate and average of resistance loss

      出现这种现象的原因是随着腐朽的持续进行,木材降解程度增加,绝干密度减小,木芯质量损失率增加。相关研究表明:随着质量损失率增加,木材的弹性模量和抗弯强度等力学性能也会显著下降[12]。而阻抗仪记录的是探头受到的微钻阻力值,也是木材力学性能的表现。因此木材质量损失率和阻力损失有较强的线性相关关系,在本次试验中能够将微钻阻力下降的位置和相应的阻力损失作为该一维直线方向上腐朽的参考位置和相应的腐朽程度。

    • ERT的工作原理是通过激励电极在立木表面施加交变电流,在立木内部形成一定交变电场,内部电阻率的任何变化都会引起立木表面电位的变化,通过测量电极收集立木周围的电位变化信息,根据一定的算法重建立木截面上的电阻率分布,得到截面上的电阻率二维图像[13]。其中木材电阻率受含水率、温度、密度和纹理方向等因素的影响[11]图 4为6号样木的ERT图像,能够展现立木二维截面上的电阻率分布,蓝色越深对应电阻率越小,红色越深则相反。ERT图像与微钻阻力变化曲线均在距地130 cm处测得。

      图  4  ERT图像

      Figure 4.  ERT image

      前文已经证明了阻抗仪检测能够准确反映活立木一维直线方向上的腐朽情况,因此以微钻阻力变化曲线为参考,用同一方向上的电阻率变化与其对比。在ERT图像上根据颜色渐变手工提取东西和南北方向的电阻率值并制成曲线,图 2中的虚线部分为6号样木东西方向上的电阻率变化曲线。对比阻抗仪记录的微钻阻力变化曲线和相应的电阻率变化曲线,能够发现发生微钻阻力值明显下降的位置,相应的电阻率也有明显变化,如图 2EF段电阻率陡然增大而GH段电阻率明显减小。因此这种现象定性地说明腐朽部分的电阻率值会变得与健康材质的电阻率明显不同,而变化趋势相反是因为这两处的腐朽程度不同或存在不同类型的缺陷。一般来说,与健康材的电阻率相比,处于前中期腐朽的木材的电阻率会发生明显下降[14],而非常严重腐朽或空洞部分的电阻率则会显著增大[15]。值得一提的是,由于阻抗仪检测是将探头从树皮逐渐打入木质部,因此得到的微钻阻力变化曲线包含一定厚度的树皮部分;而ERT检测则是将探测器直接钉入木质部,得到的图像是不包括树皮部分的,所以图 2中的电阻率与微钻阻力的曲线会受树皮厚度的影响,图 2只能用来定性判断微钻阻力下降与电阻率变化的关系。

      因为前述的阻力损失Z与质量损失率L有很好的线性关系,能够通过微钻阻力下降的位置和相应的阻力损失大小判断腐朽的位置和腐朽程度,所以在ERT图像上的相同直线方向上标记出微钻阻力下降的位置作为腐朽发生的位置,如图 4上的EFGH段。在ERT图像上判断腐朽时,为了消除树皮厚度对位置判断的影响,先从样木的胸径得出其半径(r=D2/2),再减去微钻阻力变化曲线上腐朽位置的两个坐标以消去树皮厚度,得到腐朽位置的圆心距,列于表 3。以6号样木为例,E点的圆心距d1=r-OE, F点圆心距d2=r-OF, 负值表示在圆心另一侧。再根据圆心距标记ERT图像上一维直线方向上的腐朽位置。在微钻阻力变化曲线范围内将其他部分标记为健康部分,如图 4中除了EF段和GH段以外的线段。为了研究腐朽部分与健康部分的电阻率差异,定义电阻异常程度V(%)如下:

      $$ V=\frac{|{{R}_{1}}-{{R}_{2}}|}{{{R}_{2}}}\times 100\% $$

      式中:R1为微钻阻力变化曲线中波谷对应线段(视为腐朽部分,如图 4中的EF段或GH段)的电阻率均值;R2表示在阻抗仪的测量长度范围内,除了波谷以外的区间(视为健康部分,如图 4中的除了EF段和GH段以外的部分)对应的电阻率均值。考虑到不同的样木的生长阶段和立地环境有差别,为了尽可能消除这些差别导致的电阻率差异,用每棵样木自身的腐朽部分与健康部分的电阻率差来抵消这些因素的影响,以便将不同样木腐朽与健康部分的电阻率差异进行比较。

      表 3为15棵样木东西和南北方向上的阻力损失及对应的电阻率异常程度。阻力损失值由微钻阻力变化曲线上的波谷(微钻阻力的明显下降)计算而来, 部分样木只在一个方向的微钻阻力变化曲线上有明显的波谷,波谷的数量和位置(以圆心距表示)同样列于表 3。以电阻率异常程度为自变量,阻力损失为因变量做两者的散点图,非线性回归得到两者呈显著的指数关系,回归方程为Z=2.813+0.182e0.039V(R2=0.755,P<0.01),结果如图 4所示。从图上看出除了极少数几个点的电阻率异常程度小于5%,绝大多数的电阻率异常程度都大于10%。说明微钻阻力发生明显下降的腐朽部分,其对应的电阻率与健康部分相差普遍较大,能够根据ERT图像上电阻率的异常部分判断是否发生腐朽,即该部分的电阻率是否明显异于周围部分。这与我们从两者的曲线变化关系图中观察到的情况一致,也说明阻力损失和电阻率异常程度存在联系。

      图  5  电阻率异常程度和阻力损失散点图

      Figure 5.  Scatter-plot of abnormal extent of electrical resistivity and resistance loss

      回归曲线前半部分(V≤90%)呈现很平缓的上升。这段曲线的阻力损失都小于10 resi,对应着阻抗仪测试方向上有不连续而且腐朽程度较小的部分,其力学强度受腐朽的影响不大;而这些部分的电阻率均值与没有腐朽的部分相差从0到近100%,变化幅度非常明显。说明力学性能稍微下降的轻微腐朽部分的阻力损失都很小且彼此之间差别不大,此时阻力损失只能判断是否有轻度的腐朽发生;虽然轻微腐朽部分的力学性能相差很小,但影响其电阻率的某些微观因素(如离子浓度)却会发生很大变化[14],导致几乎相同的阻力损失值会对应相差很大的电阻率异常程度。所以此时电阻率异常程度和阻力损失间的正相关关系不明显,前者相比后者能更细致地反映轻微腐朽的程度,电阻率检测对轻微腐朽更加敏感。曲线的后半部分(V>90%)的斜率明显增大,呈现快速上升趋势。这段曲线的阻力损失明显变大,对应测试方向上有连续且腐朽程度较大的部分,其力学强度显著下降;这些部分的电阻率异常程度的变化不再明显,从100%到130%。说明当腐朽严重到力学性能发生大幅下降时,影响木材电阻率的因素的变化速度减慢,电阻率异常程度有小幅的增加,阻抗仪检测对较严重的腐朽更加敏感。这与相关研究所得到的电阻法对早期腐朽敏感,而阻抗仪法对严重腐朽的测量最为准确的结论相符[16]

      这一现象是因为两种检测方法的原理不同, 导致它们对腐朽的反应不同。活立木刚遭到腐朽菌感染后,其纤维素含量相对不变[17],腐朽菌会先消耗细胞内含物及细胞间隙物质,力学性能变化不大,所以曲线前部分的阻力损失值变化很小。随着腐朽持续进行,腐朽菌逐渐把细胞壁中的纤维素、半纤维素分解为糖类,使木材纤维结晶区的结晶度降低[18],这个过程由点到面逐渐进行,导致其力学性能发生较明显下降[19],这与曲线上斜率明显增大相对应。腐朽进行到后期时,细胞壁结构完全水解,木材组织遭到严重破坏,宏观结构发生变化, 如发生筛孔或轮裂等[20],所以曲线后部分的阻力损失很大。此外,裂纹和空洞缺陷也会使微钻阻力值陡然下降[21],在微钻阻力变化曲线上形成波谷, 其阻力损失大小由缺陷的尺寸决定。

      相关研究表明木材电阻会在早期腐朽阶段大幅下降,因为在早期腐朽阶段,含水率虽然会上升但其并不是影响电阻率的主要因素,而作为主要因素的离子浓度会在木材腐朽菌作用下显著增加,导致其电阻率大幅下降[22]。所以曲线前部分的电阻率异常程度变化很大,表明此时腐朽使木材电阻率变化很快。而曲线的走势发生明显变化后,电阻率异常程度的变化程度比较小。此时腐朽部分的电阻率虽然比健康部分低很多,但是随着腐朽程度的增加其变化不再明显,说明木材电阻率在中后期腐朽时的变化则趋于平稳。这种现象与相关的腐朽对木材电阻影响的研究结果相符,其原因可能是木腐菌的分解作用会随着腐朽菌生长趋于平稳[14]。至于微钻阻力下降段的电阻率相对于健康部分增大的原因有两种:一是在腐朽非常严重的位置,树木组织失去功能,得不到水分供给,电阻率因水分减少而显著增大[15],此时电阻率异常程度和阻力损失值都很大;二是发生裂纹、空洞等缺陷部位的电阻率也会明显增大,这种情况下测得的电阻率其实近似为绝缘的空气电阻率,对应的电阻率异常程度都很大,而阻力损失的大小由缺陷的尺寸决定。

      以微钻阻力变化曲线上有明显阻力下降的位置为腐朽的参考位置,用这些位置与健康部分的电阻率均值的差值定义了腐朽部分的电阻率异常程度。通过非线性回归发现阻力损失和电阻率异常程度呈正相关指数关系,且电阻率异常程度基本都大于10%,说明ERT图像上电阻率高于或低于周围区域的位置存在不同程度腐朽, ERT能够反映红松活立木二维截面上的腐朽情况。两者的回归曲线可以分为变化趋势相差很大的两部分,由此推断出ERT检测对于轻微腐朽更加敏感而阻抗仪检测对较严重的腐朽更加敏感。

    • 由于微钻阻力变化曲线上的波谷位置与ERT图像上电阻率极大值或极小值的位置相对应,电阻率异常程度与阻力损失存在较强的相关关系,因此可以认为ERT图像上电阻率的极值位置存在不同程度的腐朽或其他缺陷。然而ERT图像上电阻率是逐渐变化的,仅根据图像上电阻率大小所对应的颜色深浅只能对腐朽区域进行粗略定位。以图 4所示的6号样木的ERT图像为例,EF段及其周围的深红色高电阻率区域为电阻率呈极大值的腐朽区域,该区域的边界位置的黄色和浅黄色区域是否为腐朽却无法仅通过电阻率判断。此外,不同样木的生长阶段和立地条件总有差异,导致不同样木截面上影响电阻率的含水率、温度和密度等众多因素也有所差异。即使是在距地面同样高度的截面上,不同样木的腐朽位置测得的电阻率也不同。因此很难仅根据ERT图像界定腐朽电阻率的统一范围,无法根据腐朽区域的电阻率范围在ERT图像上明确确定腐朽的位置。

      为了解决这个问题,考虑两种测量方法的优势,结合阻抗仪测得的微钻阻力变化曲线和ERT图像进行分析。在ERT图像上找出阻抗仪的测量方向及范围,用细的直条表示阻抗仪的测量方向和范围,微钻阻力值发生比较明显下降的位置的腐朽部分标记为白色,其他位置的健康区域标记为黑色,依据是阻抗仪能够在一维直线方向上准确地反应活立木的腐朽情况。然后运用监督分类的方法对ERT图像进行分析,以选出的腐朽部分为训练样区,区分腐朽部分与健康部分。监督分类法是以统计识别函数为理论基础,根据已知训练区提供的样本通过选择特征参数,求出特征参数作为决策规则,建立判别函数以对待分类影像进行的图像分类[23]。这样处理的依据是:ERT检测对轻微腐朽很敏感,即使是轻微的腐朽,其电阻率与健康部分相差较大,这些部分在ERT图像上与其他部分的颜色有较明显的差异,满足监督分类的要求。图 6为6号样木的ERT图像分析过程,左边为原ERT图像,其中白色直条为根据微钻阻力变化曲线标记的腐朽位置,黑色直条为健康位置,直条总长为阻抗仪测量长度;右图为经监督分类后图像,其中红色区域为处理后得到的腐朽区域,绿色为健康区域,蓝色为背景区域。用处理后的腐朽区域面积占总截面积的比例表示由这种方法得到的腐朽程度,具体公式为:

      图  6  ERT图像分析过程

      Figure 6.  Process of analyzing ERT image

      $$ {{P}_{\text{d}}}=\frac{{{\mathit{P}}_{\text{r}}}}{1-{{P}_{\text{b}}}} $$

      式中: Pd为腐朽区域占截面的比例,%;PrPb分别为红色所代表的腐朽区域和蓝色所代表的背景区域占整个图像的面积比,%。用这种方法对15棵样木的ERT图像进行分析,所得的腐朽面积占截面的比例见表 4

      表 4  ERT图像中腐朽面积占截面比例

      Table 4.  Percentage of decay area in ERT image

      样木编号No. of sample tree123456789101112131415
      腐朽面积比例Percentage of decay area (Pd)/%33.6040.8828.9034.0926.8937.6523.2421.4242.5319.5955.2165.2522.3843.8175.37

      以质量损失率为自变量,上述方法得到的腐朽面积占截面的比例Pd为因变量,对两者进行线性回归分析,结果如图 7所示。回归方程Pd=0.917L-8.616(R2=0.787,P<0.01)的拟合程度很高,说明木芯的质量损失率越大,结合两种检测方法得到的腐朽面积越大,两者呈显著正相关关系。这一现象是由于ERT图像上监督分类的结果是用阻抗仪的微钻阻力变化曲线定位得到的,红色区域实质上是一维直线上力学性能下降的区域在二维截面上的扩展,扩展的依据是前文讨论的腐朽部分的阻力损失与电阻率异常程度之间存在相关关系,微钻阻力下降段的腐朽位置电阻率与其他健康位置的电阻率相差较大。而质量损失率作为本文中腐朽程度的参考值,已经证明它与阻力损失有线性正相关关系,会随力学性能下降而减小。因为计算得到的腐朽面积比例Pd和质量损失率L都与力学强度下降有密切关系,所以两者表现出较好的线性相关关系。这种现象也印证了腐朽会使木材的力学和电学性能都发生变化,这种联合阻抗仪和ERT图像表征腐朽程度的方法具有可行性。

      图  7  质量损失率和ERT图像中腐朽面积比例散点图

      Figure 7.  Scatter-plot of mass loss rate and percentage of decay area in ERT image

      对比阻力损失与质量损失率的回归方程,此方法得到的腐朽程度与质量损失率的相关系数更大,拟合程度提高。说明与一维方向上的微钻阻力变化曲线相比,结合微钻阻力变化曲线对ERT图像分析后得到的二维截面图像,能更加全面地反映活立木内部的腐朽情况。而且此方法能明确地表示出红松活立木截面上的腐朽区域,充分利用了阻抗仪检测在一维方向上的准确性以及电阻率对轻微腐朽敏感的特点,克服了仅根据ERT图像难以明确确定腐朽区域的困难。

    • 对ERT图像和阻抗仪两种无损检测方法进行了比较和评价,结合两种方法的优势,并利用监督分类的图像分析法对活立木树干内部腐朽进行了检测和定量表征,结果表明:

      1) 根据微钻阻力下降的程度定义阻力损失,木芯的质量损失率与平均的阻力损失呈较显著的正相关关系,阻抗仪能够准确测量红松活立木一维直线方向上的腐朽情况。

      2) 微钻阻力下降区域的电阻率与其他区域相差很大,以这些区域为腐朽区域,用腐朽区域和健康区域的电阻率值定义电阻率异常程度。阻力损失值和对应区域的电阻率异常程度存在显著指数关系,说明了ERT能够检测木材二维截面上的腐朽情况。但是仅从ERT图像不能明确得到腐朽的区域,只能根据电阻率过高或过低大致判断腐朽部位。由上述指数函数关系还得到ERT检测对于轻度腐朽更加敏感,而阻抗仪检测对较严重的腐朽更加敏感。

      3) 结合阻抗仪和ERT图像检测各自的优点,根据微钻阻力变化曲线上阻力下降的位置,在ERT图像对应的一维方向上确定腐朽位置,运用监督分类的图像分析手段在红松活立木二维截面上明确确定腐朽区域,腐朽面积占截面比例与质量损失率呈显著正相关关系,证明了这种联合阻抗仪和ERT定量表征红松活立木腐朽的方法有相当的可行性。

参考文献 (23)

目录

    /

    返回文章
    返回