基于雷达检测的木材内部孔洞面积测算与修正

陈勇平，郭文静，唐启恒

（中国林业科学研究院 木材工业研究所，北京 100091）

雷达探测木材内部异常结构主要通过分析发射的电磁波的反射信息来实现。当木材内部存在孔洞时，因孔洞与木材本身具有不同的介电性能，进入木材的电磁波在穿过内部孔洞时会在两者的界面上发生反射。分析反射回波波形资料，可以推断木材内部孔洞情况。利用雷达对木材进行扫描可以快速获知木材内部是否存在孔洞。雷达波在不同含水率不同树种等条件下具有不同的传播速度，其检测图像若采取统一的方法进行分析会存在一定误差。大量研究表明：木材孔洞的探测精度与被测对象的介电常数密切相关，而木材的含水率、纤维密度和纤维方向等均会对介电常数产生一定影响[1]，其中木材含水率与介电常数密切相关[2]，因此分析木材含水率与介电常数的关系对开展雷达检测木材内部缺陷研究非常重要[3]。由此也可逆向通过雷达回波信号进行木材含水率的预估[4]。在实际应用过程中，针对雷达探测木材内部缺陷存在解析困难、探测精度不高、可视化程度低等不足，国内外学者展开了广泛研究。可行性及可靠性方面，LORENZO等[5]研究了雷达波对树干以及树根检测的可行性；邸向辉等[6]指出：电磁波信号解释分析的困难性是树木无损探测的主要制约因素；JEZOVA等[7]开展了树木内部残损的雷达图像特征研究；WU等[8]评估了探地雷达技术在探测不同树种内部损伤方面的可靠性。成像精度及算法方面，DEVARU等[9]阐述了探地雷达在树木无损检测方面的优点，并设计试验研究了探地雷达对树干内部缺陷的自动检测算法；BUTNOR等[10]采用阈值法对针叶树木内部缺陷进行了探测，并获得了较好的效果；李伟林等[11]对颐和园界河桥旁6株古柳Salixspp.进行雷达探测，并采用匹配滤波器算法计算回波时延以提高检测精度；肖夏阳等[12]通过比较阈值法、匹配滤波器法和希尔伯特积算法在树木内部层位探测的准确性，得出希尔伯特积算法对树木内缺陷区域的检测更加准确。陈勇平等[13]研究也发现：木材含水率等对雷达检测孔洞面积的影响较大，并开展树种、含水率等对雷达预估孔洞面积的初步分析。本研究旨在探究基于雷达检测的木材内部孔洞面积预估修正方法，提高雷达检测精度，从而为木材内部残损检测提供理论和数据支撑。

1 材料与方法

1.1 材料

马尾松Pinusmassoniana木材采集于湖南省攸县，造林年度为1966年，树龄50 a，胸径约51 cm，树高约22 m。杉木Cunninghamialanceolata木材采集于湖南省攸县，造林年度为1982年，树龄34 a，胸径约46 cm，树高约25 m。

1.2 试材加工

马尾松和杉木原木伐倒后，自树根位置80 cm起按每段长约80 cm逐个截取木段，以模拟短柱；截断后的部分木段通风干燥3 a，使其达到北京地区的平衡含水率。对干燥后的木段进行挑选并分别编号为S1、S2、S3、S4、S5、S6和C1用于实验研究，每个木段标记其起始检测点，其中S1～S6为马尾松试件，C1为杉木试件。取S2、S3、S4木段在其内部加工大小不一致的圆形孔洞；取S5木段在其内部加工三角形孔洞；取S6木段在其内部加工四边形孔洞；取C1木段在其内部加工圆形孔洞；S1为备用木段。以上木段的孔洞加工均位于中心位置，孔洞深度约30 cm(图1)。

1.3 测试方法

采用树木雷达对上述试材进行雷达扫描[13]，其中S3、S4、S5、S6为试验用材，S2和C1为数据分析验证试材。本次测试研究采用美国树木雷达检测系统[Tree Radar Unit(TRU)]，主要参数设置为：频率900 MHz，最大扫描直径4.2 m，深度分辨率1.9 mm，数据采集间隔5 mm。

2 结果与分析

2.1 扫描结果初步校正

采用雷达扫描技术对带孔试材进行检测。以S4试件为例，通过分析雷达波反射波形和检测图像可以直观判定木材内部孔洞存在(图2A)，但雷达预测孔洞面积(图2B)与实际孔洞面积(图2C)存在较大的偏差。

基于直观成像的误差，且希尔伯特积算法相较于阈值法和匹配滤波器法对木质体内部缺陷识别效果更好[12,14]，本研究首先对上述检测结果进行希尔伯特积算法修正，其部分检测结果如图3所示。本研究采用层剥反演法获得木材介电常数并对其内部孔洞进行定位，结合检测对象外轮廓数据，可实现内部孔洞较准确定位。

图 1 孔洞加工示意图Figure 1 Schematic of internal cavity in wood columns

图 2 雷达成像直观结果Figure 2 Intuitionistic results of the cavity areas in wood columns by TRU

图 3 希尔伯特积算法演算结果Figure 3 Analysis results of the cavity areas in wood columns by Hilbert

结合图2和图3可以看出：利用希尔伯特积算法对雷达检测图像进行修正，提高了检测精度，但还是存在一定误差，修正面积仍低于实际孔洞面积，其中试件S4修正图像的相对误差为20.57%，试件S6修正图像的相对误差为41.49%。这主要是因为希尔伯特积算法是根据层剥反演获得相对介电常数，从而对木材内部的传输速度进行在线估测，而实际过程该介电常数可能与估算的介电常数存在一定得误差。为此，本研究根据希尔伯特积算法原理，采取提取单路径检测数据，并使用曲线偏移、介电常数修正、辅助系数校正等方法尝试进一步的修正研究。

2.2 单路径数据提取分析

根据雷达检测原理，树木雷达工作时雷达天线绕该木柱一周相当于虚拟钻孔机对该木柱进行了360°的“钻孔”操作。本研究对检测对象进行圆周10等分，并提取等分点上每一次的“钻孔”结果。每一次“钻孔”结果即对应雷达检测中每一个位置的回波图形。通过回波图形分析可获得对应位置木材边缘与孔洞边缘距离，其中孔洞边缘即高介电常数物质(木材)与低介电常数物质(孔洞)交界面，交界面位置在回波图形中的特征为谷—峰—谷，叠加后在雷达检测图像中表现的形态为黑—白—黑。根据“钻孔”检测结果，将得到10个孔洞边界点，将边界点进行连线，即形成了内部孔洞轮廓。以试件S3为例，结果如图4。

图 4 雷达检测初步分析结果Figure 4 Preliminary results of the cavity areas in wood columns by TRU

图4A为试件S3；图4B为图4A的断面轮廓提取，其目的是便于检测对象横断面及孔洞面积的计算；图4C为雷达检测初步分析图像。为了便于数据比对，本研究对试件S3、S4、S5、S6的雷达初步检测结果进行了测算和统计。结果见表1。

从表1和图4可以看出：通过雷达检测可以快速判定木材内部孔洞的存在，但其预估的孔洞面积与实际面积相差较大。这主要是因为雷达检测内部孔洞的轮廓线产生了一定的偏移，且实际检测对象的介电常数与预设常数相差较大，导致预估孔洞面积与实际孔洞面积相差甚远。从表1和图4结果也可以看出：雷达检测预估内部孔洞面积若不加以修正，其数值远小于实际孔洞面积，无法直接应用于古建筑木构件内部残损的定量评估。

表 1 雷达检测初步测算与统计Table 1 Preliminary calculation and statistics of the cavity areas in wood columns by TRU

2.3 数据处理与成像精度改进

图5A为试件S3；图5B为根据偏移量进行第1次修正后的孔洞图像；图5C为根据介电常数进行第2次修正后的孔洞图像；图5D为3次修正后的孔洞图像。为了便于分析，同样对试件S3、S4、S5、S6的雷达检测数据修正结果进行了计算和统计。从表2和图5可以看出：根据偏移量进行第1次修正后，识别精度得到一定得提升，但其与实际孔洞相比均表现为偏大，这主要是因为修正距离s1c是根据雷达检测系统预设介电常数ε1计算得来，但实际检测对象的介电常数为ε2，在气干木材的检测中预设介电常数ε1大于检测对象介电常数ε2，故而s1c距离偏小，导致孔洞识别偏大。从表2和图5也可以看出：根据介电常数进行第2次修正后，其检测修正面积与实际孔洞面积仍存在一定偏差。从图5也可以看出：无论是否修正，雷达识别孔洞形状的能力均较差，这可能是因为外部形状等因素引起了内部图像的扭曲。

图 5 雷达检测数据修正结果Figure 5 Corrected results of the cavity areas in wood columns by TRU

表 2 雷达检测数据修正与统计Table 2 Corrected calculation and statistics of the cavity areas in wood columns by TRU

上述研究结果表明：单独依靠介电常数比值进行修正，其孔洞识别面积与实际孔洞仍存在一定偏差。为此，本研究设置了另一辅助修正系数k2，即雷达检测中木材边缘与孔洞边缘距离按照s3c=(sr−a)k1k2进行修正，其中辅助修正系数k2为常数。经过前期的研究分析和数据拟合，可初步确定常数k2取值为0.80～0.90，本研究k2取值为0.85。从表1和表2可以看出：雷达检测结果进行3次修正后，其识别精度得到明显提升，总体相对误差小于30%。

2.4 验证分析

为验证上述结论，本研究选取了具有不同大小孔洞马尾松试材S2和杉木试材C1进行补充试验，并根据上述研究步骤获取其断面轮廓、雷达检测初步结果、第1次修正结果、第2次修正结果和第3次修正结果，经数据整理见图6。

从验证试验可以得出同样结论，未修正或按照一定含水率下特定树种的介电常数修正的雷达图像其预估孔洞面积与孔洞实际面积均存在一定误差，按照s3c=(sr−a)k1k2对图像进行3次修正后精度较高。3次修正后，试材S2雷达检测预估孔洞面积相对误差为9.47%，试材C1预估相对误差为28.72%，检测结果可能会因为个体差异产生一定的偏差，但大致趋势相同。

3 结论

图 6 雷达检测数据修正过程Figure 6 Corrected process of the cavity areas in wood columns by TRU

通过分析雷达回波图形和检测图像可以快速判定松木和杉木木材内部孔洞的存在，但其直观成像技术在定量评估方面误差较大，需进行一定的数据处理和修正。①对带孔木材进行雷达扫描，未经修正的雷达图像其预估孔洞面积小于实际孔洞面积，按介电常数进行修正的雷达图像其预估孔洞面积与实际孔洞面积仍存在一定的误差。②无论是否进行数据处理和修正，利用雷达扫描进行内部孔洞的形状预估难度较大，难以准确识别三角形、四边形、圆形等各种形状的孔洞，应开展进一步研究。③实际应用中，根据树种和含水率所得到的介电常数应进行再次修正，即实际应用中根据树种和含水率所得到的介电常数仅为多参数校正中的一个因素。④利用该雷达检测与修正方法进行木材内部孔洞探测，综合树种和含水率等参数，可确保孔洞面积的识别误差小于30%。此项技术可应用于古建筑木构件内部残损的检测评估。⑤在开展古建筑木构件内部残损的雷达检测与定量评估前，应首先调查分析木结构建筑常用的树种类型，并建立常用树种、含水率与介电常数之间的关系模型。