-
木材无损检测技术可以对木材的力学性能、内部腐朽进行准确检测,对充分利用木材资源有着重要的意义。应力波技术是一种非常有效的木材无损检测方法[1-2],其原理是利用脉冲锤敲击木材,使其内部产生应力波的传播,通过测量应力波的传播时间和传播速度的变化确定木质材料的性能(如弹性模量、腐朽等)。国外诸多科研院所以及企业均投入大量精力,新产品层出不穷。孟瑞华等[3]、冯国红等[4]、肖江等[5]、Hayes等[6]采用不同方法,实现了应力波传播时间、传播速度的测量。在这些研究中,大多通过应力波信号的阈值比较或过零比较等操作,计算应力波在木材中的传播时间。由于应力波在木材中传播过程复杂,这种通过简单的比较运算获取的传播时间分散性较大,抗干扰能力弱。本研究提出一种高阶统计量[7]和小波变换[8]相结合的广义相关算法,通过抑制应力波信号中的噪声,提高应力波传播时间的测量精度。
-
互相关法利用空间上2个独立传感器接收到的信号获取时间差。根据应力波木材无损检测的原理,建立信号模型如下:
$$\begin{align} & {{y}_{1}}\left( n \right)=x\left( n \right)+{{s}_{1}}\left( n \right) \\ & {{y}_{2}}\left( n \right)=Ax\left( n-D \right)+{{s}_{2}}\left( n \right) \\ \end{align}$$ (1) 式(1)中:y1(n),y2(n),为2个传感器采集到的信号,x(n)表示纯净应力波信号,A为应力波信号在木材中传播过程中出现的衰减系数,D表示应力波从1号传感器传播至2号传感器所需时间,s1(n)和s1(n)表示y1(n)和y2(n)中混杂的噪声信号。
y1(n)和y2(n
)的互相关函数为: $${{R}_{12}}\left( \tau \right)=E\left[ {{y}_{1}}\left( n \right){{y}_{2}}\left( n-\tau \right) \right]={{R}_{xx}}\left( \tau -D \right)+{{R}_{xs1}}\left( \tau -D \right)+{{R}_{xs2}}\left( \tau \right)+{{R}_{s1s2}}\left( \tau \right)$$ (2) 式(2)中:Rxx为x(n)的自相关函数,Rxs1为x(n)和s1(n
)的互相关函数,Rxs2为x(n)和s2(n)的互相关函数,Rs1s2为s1(n)和s2(n)的互相关函数。 由于x(n)与s1(n),s2(n)之间互不相关,所以存在:
$$\begin{align} & {{R}_{xs1}}\left( \tau -D \right)=0; \\ & {{R}_{xs2}}\left( \tau \right)=0; \\ & {{R}_{s1s2}}\left( \tau \right)=0 \\ \end{align}$$ (3) 式(3)定义的互相关函数等于纯净应力波信号的自相关函数,即:
$${{R}_{12}}\left( \tau \right)={{R}_{xx}}\left( \tau -D \right)$$ (4) 由自相关函数性质可以得出$\left| {{R}_{xx}}\left( \tau -D \right) \right|\le {{R}_{xx}}\left( 0 \right)$,那么当τ-D=0时R
12(τ)取得最大值,通过寻找互相关R12(τ)函数峰值对应的时间可估计出应力波信号y1(n)传播至第2个传感器所需的时间D。 广义相关法是Knappy等[9]和Hero等[10]提出的,它是在基本相关法基础上对信号提前进行处理,提高信噪比,从而有效提高时延估计精度。图 1描述了广义相关法的原理,它的步骤包括对信号去噪、求两路信号的互相关函数、计算互相关函数的峰值时刻等,最后得到应力波信号的传播时间。
-
在广义相关算法中,需利用滤波算法对信号去噪以便提高时延估计的精度。由于在高阶统计量中,高斯噪声的双谱为0,能够彻底消除信号中的高斯噪声。小波变换又能够很好地保留信号的峰值特征和信号信息。因此将高阶统计量和小波变换结合的滤波方法具有更好的降噪滤波功能。
算法具体步骤如下:①对含噪信号进行小波变换,将原始波形分解为高频部分和低频部分;②使用小波分层阀值降噪对低频小波系数进行处理;③处理高斯噪声后的小波系数仍是高斯分布的,对高频部分进行高阶统计量的双谱估计法进行滤波;④利用小波反变换对去噪后的波形进行重构。
-
本节提出基于小波变换与高阶统计量的广义互相关应力波传播时延估计方法。时延估计的输入信号为初始应力波信号在经过高阶统计量和小波变换相结的滤波算法滤波去噪后的信号,然后对输入信号求互相关函数,找到互相关函数峰值对应的时刻。算法原理如图 2所示。
图 2 基于广义互相关的应力波时延估计算法原理
Figure 2. Scheme of wood nondestructive detection of stress wave time delay estimation
算法实现步骤如下:第1步,数据采集。利用搭建好的实验平台采集所需的2列应力波信号,通过Openchoice软件以Excel数据保存。第2步,数据处理。剔除有用信号之前的大量的0数据,保留长度为2 048的信号,分别记为y1(n),y2(n)。第3步,滤波。利用Matlab软件编写程序对波形进行3层尺度分解,分解后得到高频部分cD1,cD2,cD3和低频部分cA1,cA2,cA3的小波系数。对cA1,cA2,cA3进行阀值处理,阀值选取采用小波默认阀值。对cD1,cD2,cD3则采用高阶统计量中的双谱滤波法进行滤波,因为高阶统计量中的高斯噪声双谱为0。将处理完后的小波系数进行重构,反变换得到去噪波形$y_{1}^{'}\left( n \right),y_{2}^{'}\left( n \right)$。第4步,时延估计。求互相关函数R12(τ)。利用Matlab软件的工具箱函数max( )求出R12(τ)的最大值及最大值所在位置D。应力波信号的传播时间t传播就可以根据下式求出:
$${{t}_{传播}}=D/{{f}_{s}}$$ (5) 式(5)中fs为信号采集时所用的采样频率。
-
把原始信号y1(n)(抽样频率100 kHz,数据长度为2 048)向右平移10个点,得到延迟信号y2(n),用于测试提出的算法对应力波传播时延估计的效果。
首先进行滤波处理,去除噪声信号。本研究利用Matlab工具箱设计了实验。先用小波工具箱中的wavedec( )函数对小波信号进行3层分解,提取小波系数。输出部分包括分解结构中小波分解矢量C和相应的记录向量L。随后用小波变换和高阶累计量相结合的滤波方法对含噪信号进行滤波处理,得到小波系数。最后在Matlab软件平台上用互相关法求出应力波传播时延的估计值为100 μs。
为了比较本研究提出的广义相关法与传统相关法的效果,还利用相关法直接估计y1(n)和y2(n)的时延,估计结果为80 μs。因为y2(n)是y2(n)右移10个采样点得到的,滤波后的互相关算法比不滤波直接求互相关函数的算法更加准确,时延估计精度提高了81.5%。
-
时延估计算法引入应力波木材无损检测中,提高时间测量的准确性。利用项目组自行研发的Wopeck木材无损检测仪完成了本组实验,实验场景如图 3和图 4所示,Wopeck木材无损检测仪采用传统计时器计时方法。根据Wopeck仪器使用方法,将12个传感器均匀地安装在腐朽木材周围并根据编号记录下来,分别敲击传感器,利用Wopeck仪器测量出各对传感器之间的传播时间,并产生木材内部缺陷的断层图像,如图 5所示。
图 5 Wopeck检测仪(传统计时法)的检测结果
Figure 5. Experimental results of Wopeck tester (traditional timing method)
为了验证本研究提出的时延估计方法,拆除木材上的传感器,根据事先记录下来的编号,利用数据采集平台的2个传感器分别采集2个测量点的信号波形。然后在Matlab软件上用基于小波变换和高阶统计量滤波的互相关时延估计算法分别求取时延值,导入Wopeck仪器的软件上,生成木材内部缺陷断层图像。如图 6所示。
从图 5和图 6可以看出:本研究提出的时延估计法和传统的计时器计时法都能判断出木材的腐朽位置和程度,但采用本研究提出的时延估计方法后,生成的图像中腐朽位置更加突出,颜色对比更加明显,测量结果的准确性更高。
-
在实际的木材无损检测过程中,通常会遇到震动和声音对测量信号等干扰,为了验证本算法的抗干扰能力,采用高分贝喇叭在传感器附近模拟外界声音,并用震动机模拟来自地面的震动。按照高、中、低三级信噪比完成了3组不同实验。3种情况下采集到的信号波形如图 7所示。
采用本研究提出的方法估计不同情况下的采样点延迟以及传播时间如表 1所示:从表 1可以看出:3种不同信噪比下的延时估计方差仅为2/3×10-8,说明本研究提出的时延估计算法性能稳定,能够适应木材无损检测环境下不同等级的信噪比。该算法法提高了应力波传播时间的测量精度,使木材无损检测的准确性得到进一步的提高。
表 1 不同干扰情况下传播时间估计结果
Table 1. Results under there levels interference
干扰强度 采样点延迟(D) 传播时间(t传播)/s 低 1 059 0.021 2 中 1 052 0.021 0 高 1 053 0.021 0 -
为了克服木材无损检测技术中传统计时方法的不足,本研究提出了一种小波变换和高阶统计量相结合的滤波算法以及广义互相关时延估计算法,用于木材无损检测中的应力波传播时间估计。大量的实验结果表明:该算法能够准确测量木材内部应力波传播时间,断层成像效果好,并具有较强的抗干扰能力。
Stress wave propagation time estimation algorithm for nondestructive testing of wood
-
摘要: 为了准确测量应力波在木材内部的传播时间,提高木材内部缺陷识别的精度,提出了一种基于小波变换和高阶统计量的广义相关时延估计算法。该方法利用滤波算法,将采集到的应力波信号进行滤波去噪,然后输入到广义互相关时延模型,求出传播时延估计值。仿真和实际的木材检测实验结果表明:提出的算法比传统的定时器计时法时延估计精度提高81.5%,并具有较强的抗干扰能力。Abstract: The purpose of this study is to measure the stress wave propagation time and improve the accuracy of identifying the internal decay of wood. An algorithm of stress wave propagation time estimation combining high order cumulants with wavelet transformation was presented. The proposed algorithm removed the noise of stress wave signal with filtering algorithm, then input the filtered signal to the generalized cross-correlation time delay model, and calculated the estimation of the stress wave propagation time. The results of simulation and experiments showed that the accuracy of the measured stress wave propagation time was improved by 81.5 per cent. Further more, under low, medium and high signal to noise ratios (SNR), the experimental results demonstrated the anti-interference ability of the new algorithm.
-
植物的蒸腾作用是其水分利用的主要方式,拉动水分在“土壤—植物—大气”连续体体系中不断循环迁移。蒸腾耗水与植物生命表征直接联系,决定着植物的水分盈缺和灌溉与否[1]。对活立木蒸腾耗水的准确测定,可以为低耗水树种选择、合理密度配置以及城市园林绿化等工作提供理论依据和参考[2]。树木蒸腾耗水产生的水势差会拉动水分通过木质部向上运输进而形成液流,因此树干液流可作为评估树木蒸腾耗水能力的一项重要指标[3-4]。目前,有多种方法可以评估树木蒸腾耗水能力,多数是通过测定树木液流速率来估算蒸腾耗水量和耗水能力。不同树干液流测定方法测量精度不同,在选择树干液流速率测量方法时需要考虑实验研究目的、活立木树种的生理条件和实验研究所处的自然环境等因素。目前液流速率的测量主要有同位素示踪法和热技术法[5] 2类。其中同位素示踪法通过将化学同位素作为示踪剂注射到树木木质部,从而检测树木液流速率;但该方法在野外应用不便,且测定精度较低,有待改进[6]。利用热技术法测定树干液流不受外界环境和树木自身结构影响,安装布置操作相对简易,并且对树木组织结构损伤较小,具有一定的应用优势[7],因此被广泛应用于树干液流测定、液流速率与环境因子的关系研究中[8-10]。如王檬檬等[11]应用热技术法研究了晋西黄土区苹果Malus pumila树液流速率与太阳辐射、大气水份等的关系;温淑红等[12]应用热技术法分析了宁南黄土陵区山桃Amygdalus davidiana树树干液流速率与太阳辐射、温度、风速的关系;杨洁等[13]应用热技术法研究了树干液流时滞效应,并精确估算树木的蒸腾耗水;还有学者[14-15]应用热技术法探究树干木质部径向不同深度的液流速率和不同时间尺度下的液流速率特征等。目前常用的测量树干液流的热技术法主要有热脉冲法(Heat Pluse Velocity Method,HPVM)、热平衡法、热扩散法、热场变形法以及外热比法等,前3种方法在国内应用较多,后2种方法在国外有较为详细的应用描述,但在国内的研究有限。鉴于此,本研究综述了现有的树干液流无损检测方法,阐述这些方法的基本原理、装置布置、应用领域和最新研究案例,对不同热技术方法的测量精度、适用范围、潜在优势以及今后改进方向等方面进行讨论和比较,并对不同研究目标和实验条件下适用的测定方法给出建议,展望各自在液流研究方面的应用前景。
1. 基于热技术的液流检测方法
1.1 热脉冲法
热脉冲法最早由HUBER等[16]提出,首次用热作为液体流动的示踪剂,利用热脉冲测量植物液流速率。该装置由加热元件和2个热电偶组成探针块,通过测量加热器发出热脉冲随着液流上升到达热电偶处所需的时间,计算液流速率(图1A)。由于热传导和对流会使得测量结果偏高,MARSHALL[17]利用热流方程建立了热脉冲法的模型框架,为热脉冲法的进一步发展提供了理论基础。
基于脉冲加热的方法包括补偿热脉冲法、最大温差法(T-max法)以及热比率法。其中补偿热脉冲法(Compensation Heat Pulse Method,CHPM)通过测量2个对称放置在线性加热器两侧的温度传感器达到相同温度时的时间来计算液流密度,该装置安装探针时会对周围木材组织造成损伤而导致液流速率失真,因此需要根据不同探针间距设置校正参数[18],从而使液流速率的测量值更接近真实值。T-max法[19]的装置由加热器和1个温度探针组成,通过记录从发出热脉冲至温度探针到达最大温度时的时间,再根据MARSHALL的基础理论计算液流速率。该方法设备简单,仅需确定被测树干边材的导热系数,即可计算树干液流密度。热比率法(Heat Ratio Method,HRM)基于补偿热脉冲法提出[20],测量范围可以达到零值甚至延伸至负值,其装置由1个加热器和2个安装在加热器上下游的温度探针组成,通过测定2个探针的增温比即可计算液流速率。
1.2 热平衡法
热平衡法(Heat Balance Method,HBM)的原理是当树干内通过一定量液流时,加热元件作为热源会向树干提供已知的热量,直至树干温度趋于稳定。若不考虑热传导以及隔热层损失热量,热源提供的热量应与被液流带走的热量相等,可根据这种热平衡关系计算液流速率。热平衡法可分为茎热平衡法和树干热平衡法[21]。
1.2.1 茎热平衡法
茎热平衡法(Stem Heat Balance,SHB)[22]以环形加热元件作为热源,提供稳定的热量,热量散失途径包括树干液流带走、热传导向树干上下方散失和对流散失。装置(图1B)设计为包裹式,利用包裹式隔热层(通常为聚苯乙烯泡沫)防止热辐射造成的热量散失(树干周围的热辐射忽略不计)。加热元件上下方安装2对热电偶,用来测定液流通过后的温差,依据热量平衡关系计算液流。SHB法适用于测定胸径较小的树干,其优点是检测时不需要标定,不需要将热电偶插入树干中,对树木无直接损伤。
1.2.2 树干热平衡法
树干热平衡法(Trunk Heat Balance,THB)[23]的原理与茎热平衡法类似,均通过热量平衡关系计算液流,不同点是THB法测量装置(图1C)由插入树干的加热片和1对热电偶组成,2个热电偶分别安装在紧挨加热片上端(温度场最大值)和加热片下端(不受温度场影响)位置,通过记录液流通过前后树干温度差来计算液流。THB法同样不需要标定,并且可测定胸径较大的树干。但THB法设备较多,安装相对复杂,易对树干造成微损伤。目前热平衡法的应用较为广泛,通常用来研究环境因子与液流速率的关系以及耗水特性。
1.3 热扩散法
热扩散法(Thermal Dissipation Method,TDM)[24]又称Granier法,是目前应用最广泛的液流测定方法。该装置(图1D)包含2个传感器探针,沿液流方向插入树干中。下游(上部)探针包括加热元件(长约20 mm),并缠绕在装有热电偶的钢针上,热电偶尖端正对加热元件的中间;下游(上部)探针不加热,用作参考探头,以测量木质部的环境温度。工作时下游探针以恒定功率(0.2 W)连续加热,受液流的散热影响,2个探头间存在温度差异,因此可通过温差与液流速率间的关系计算液流速率。
TDM法常用来研究树干液流与环境因子的关系。万艳芳等[25]应用热扩散技术测定并分析了青海云杉Picea crassifolia树干液流密度与环境因子的关系,确定液流密度的主要环境影响因子是太阳辐射。朱敏捷等[26]利用热扩散法测定了尾叶桉Eucalyptus urophylla树干液流,研究了树干液流的方位差异以及与环境因子的关系。姚增旺等[27]应用热扩散探针测定梭梭Haloxylon ammodendron树干液流,研究了树干液流与环境因子之间的时滞效应。另外,通过测定单株树干液流还可以推算林分蒸腾量。王志超等[28] 研究了林分蒸腾耗水规律后发现:忽略夜间林分蒸腾耗水量会导致对林分蒸腾耗水量的估计不准确。
1.4 热场变形法
基于热脉冲法测量树木液流密度时,需要测量热脉冲前后温度,获得温度差,这就要求木材具有较高的热稳定性;热脉冲法两侧测量需要时间间隔,可以测得液流密度的最大值为45 cm3·cm−2·h−1,说明该方法具有一定的局限性。为解决以上问题,NADEZHDINA等[29]提出了热场变形法(Heat Field Deformation,HFD),通过记录线性加热器周围的木质部中不同径向位置的热场变化,将热场变形与树干木质部的液流密度联系起来。热场变形法液流检测系统(图1E)包括3个探针和1个加热器,其中2个探针沿轴向对称安装在加热器的上游和下游,另1个探针沿切向平行于加热器水平安装在加热器侧边,轴向探针测量对称温差,切向探针测量不对称温差。通过测定加热器周围轴向和切向的温度差来表征由树液流动而产生的热场变化,进而确定液流密度。液流密度q (cm3·cm−2·h−1)的一般计算公式为:
$$ q=3\; 600D_{{\rm{st}}}(K+T_{{\rm{s}}-{\rm{a}}})/T_{{\rm{as}}}Z_{{\rm{ax}}}Z_{{\rm{tg}}}L_{{\rm{sw}}}。 $$ 其中:Dst表示树干边材热扩散率(m2·s−1);(K+Ts-a)/Tas表示温差比率;ZaxZtg表示传感器探针间距离的校正因子;Lsw表示边材深度。K表示零液流下Ts-a的绝对值,其中Ts-a为Tsym与Tas的差。Tsym表示对称探针间的温差;Tas表示非对称探针间的温差;Zax表示轴向上游探针与加热器的距离;Ztg表示切向探针与加热器的距离。
HFD传感器也可以记录反向流量,即将Tsym改为负值。因此,计算公式转换为:
$$ q=-3 \;600D_{{\rm{st}}}(-K+T_{{\rm{s}}-{\rm{a}}})/T_{{\rm{as}}}Z_{{\rm{ax}}}Z_{{\rm{tg}}}L_{{\rm{sw}}}。 $$ 用HFD法测量液流密度,非零液流下,利用线性外推法能准确测得零液流密度,相比其他热技术方法优势显著。同时HFD法结合了对称与非对称温差测量,利用对称温差测量低液流密度较为有效,而测得的高液流密度与实际蒸腾量线性关系不显著,因此高液流密度准确性不够[30]。而利用非对称温差测量是中高液流密度准确性较高。因此,HFD法对于低液流量和高液流量都可以准确测定。
HFD法广泛应用于液流指数(the sap flow index,SFI) 的测定。SFI是植物水分状况的敏感指标,用来决定植物是否需要灌溉。SFI值通过在加热器周围轴向等距安装差动热电偶测得,是液流速率测定的原始数据之一[31]。此外,HFD法可以直接监测木质部的水分运动[32],通过沿着木质部半径的不同深度,用围绕普通线性加热器的传感器进行液流测定,具有快速响应和高度敏感的特性。NADEZHDINA[33]在对枫树Acer spp. 水运输路径的研究中,利用HFD法测定枫树木质部液流,证实了枫树的维管结构具有完整拓扑结构。
1.5 外热比法
外热比法(External Heat-Ratio,EHR)是在热比率法的基础上提出的,用外部加热元件代替插入式加热元件,其基本原理与激光脉冲法(laser heat-pulse gauge,LHPG)类似。不同点是后者用近红外激光源代替插入式加热元件,并通过红外摄像机从外部监控热量传播,热脉冲速度由温度数据确定,并与液流速率相关。HELFTER等[34]利用激光脉冲法对小茎木本植物的液流速率进行了测定,发现小茎木本植物韧皮部与木质部液流速率几乎一致。CLEARWATER等[35]首次提出了外热比法(图1F),将1个微型外部加热器(电子芯片电阻)和温度传感器(精密热电偶)粘在软木块上,并压在茎干表面。释放热脉冲后,根据2个热电偶的增温比来计算液流密度。利用外热比法可以测定灌木液流速率[36],研究植物水动力学,对直径较小的茎干具有良好的适用性。外热比法最小可测直径为5 mm,可测液流密度为0.36~50.00 cm3·cm−2·h−1,较少应用于直径较大的茎干。因此,下一步可改进EHR技术,用于测定较大茎段植物的液流密度。
2. 问题与建议
2.1 热技术方法的不足与改进
探针的使用会对树干边材造成一定的破坏,使得探针处树干边材的热均匀性改变,从而降低测量结果的准确性。GREEN等[37]用二维的“热-液流”模型确定不同伤口大小的校正因子,给出了补偿脉冲法和T-max法的校正因子表,并通过比较美洲黑杨Populus deltoides与白柳Salix alba的液流通量值与实际蒸腾速率值的关系证明了校正因子的有效性。TESTI等[38]在补偿热脉冲法的基础上提出了校准平均梯度法(calibrated average gradient,CAG),有效测定了低速液流,使用也较为简易。
LANGENSIEPEN等[39] 发现:为更好地适应小麦Triticum aestivum茎的解剖结构和热物理特性,在应用茎热平衡法测量小麦液流速率时,通过引入降噪方程可有效提高液流计的测量精度。TRCALA等[40]利用热场变形法的温度场理论,通过改变传感器的几何形状(从垂直到水平)来改善热平衡法的传感性能,实现了零液流和反向液流的测定。这种方法也被称为线性热平衡法[41],是从基础传导—对流传热方程解析得出的精确方程,不仅提高了液流测定的精度,而且基于热导率信息实现了水含量的估算。NAKANO等[42]发现:对金柑Fortunella crassifocia进行环剥处理后,可利用热平衡法测定其韧皮部和木质部的液流速率。
TDM法测定液流速率需要估算线性回归关系,确定零液流状态下的温差,而这个过程会产生一定误差[43],许多情况下准确性受到质疑[44]。因此用热扩散法确定树木的蒸腾量时,有必要对测量树种液流量估计方程进行校准[45-47]。
外热比法也存在一些不足。首先,大多数加热传感器从加热芯片的中心到两侧感温元件有一定的窄间距。随着热量沿横截面向内传播和沿茎轴上下传播,热量到达木质部导管时变得非常分散,来自液流的热比率信号会减弱。其次,加热器和温度传感器被安置在1个矩形的不导电硅酮/软木块中,无法有效隔绝环境温度对检测温度的影响,增加了液流检测结果的误差。再次,矩形加热芯片横压在圆柱形树干上,载荷不均匀,加热器元件使用窄的矩形芯片电阻,比圆形芯片更容易断开。为此,王胜[48]开发了1种新设计的EHR加热传感器,增加了加热元件至温度传感器的间距,使之更适应直径较大的茎干。改良后的装置茎干直径检测范围扩大,可用于胸径较小的树木测量。
2.2 基于热技术液流速度测量的常用方法比较
目前,基于热技术的树干液流测定方法日趋完善,不同方法具有相对应的优势和劣势。由表1可知:不同热技术方法液流测定装置均包括为加热器提供能量的能量供应单元和用于收集检测数据的数据记录仪。具体来看,热脉冲法不受环境条件以及树冠结构及根系特性的影响,装置简洁,但存在一定的灵敏度和精度问题。热平衡法无需标定,测量精度有所提高,但仅适于测定高液流密度。热扩散法是目前研究蒸腾耗水特性应用最广泛的方法,测定结果较准确,仪器成本较低,安装简单,有较成熟的商业化产品,但测定结果容易被低估。热场变形法操作复杂,应用较少,但该方法能够准确测定零液流以及逆向液流,测定精度与范围也有很大的提升。外热比法与激光脉冲法均可实现精确的零破坏检测,但仅适用于胸径较小的树干,另外,激光脉冲法装置成本昂贵,未能普及。
表 1 树干液流的测定方法对比Table 1 Comparison of methods for sap flow measurement方法 装置 优点 缺点 热脉冲法 加热器,2个热电偶 不受环境条件,树冠结构及根系特性的影响,简洁准确, 经济可行[49] 存在测定精度问题[37−38] 热平衡法 探针,加热元件 无需标定,进一步提高了测量精度[22−23] 不适用于液流速率较高的 植物[50] 热扩散法 加热探针,参照探针 测定结果较准确,仪器成本较低,安装简单,有较成熟的 商品化产品[51] 液流可能被低估[43] 激光热脉冲法 近红外激光源,红外摄像机 无须将热源插入植物茎干内,避免对茎干内组织破坏而造 成误差[34] 成本较高 热场变形法 加热器,3根探针 能够准确地测定零液流量以及逆向液流[52] 测定过程较复杂[29] 外热比法 微型外部加热传感器 精确、无损地测定胸径较小的树干中的双向液流[53] 微型外部量规的配置尚存 在问题[35−36] 2.3 活立木液流测定方法选择建议
利用热技术方法测定液流速率,可以精确估算树木蒸腾耗水量[54],但不同热技术方法的测量精度、测定范围以及适用性不尽相同,实际应用时应根据不同实验条件选择不同的热技术方法。选择热技术方法测定树干液流通常需要考虑树木胸径大小、热技术方法的误差范围、热技术方法的测定精度、热技术方法的可行性等因素。
首先,不同胸径活立木应选用不同的热技术方法,外热比法和茎热平衡法适合测定胸径较小的树木,树干热平衡法适合测定胸径较大的树木。其次,不同热技术方法可测得的液流速率范围不同,补偿热脉冲法、T-max法以及热扩散法测定低速液流的误差较大,茎热平衡法测定高液流速率时误差较大,热比率法和热场变形法测定液流范围较广。热场变形法和外热比法可以测定逆向液流,热场变形法还可以准确测定零液流。热脉冲法、热平衡法和热扩散法较成熟[55],应用较广泛,可行性较高。另外,将不同测定范围的热技术方法组合使用,可以有效提高测定精度。不同植物的液流速率不同。向日葵Helianthus annuus和玉米Zea mays等植物的液流速率相对较低,在利用T-max方法测定时,测量值总是略高于实际值[56];换成热比率法测定也不够精确,而采用T-max法与热比率法组合测量则较为准确。
3. 热技术方法应用展望
利用热技术方法测定液流速率约80多年的研究,方法不断得到改进与创新,在校准度、测定范围和测定精度上均有所提高,同时,实验操作不断简化,数据实现自动化采集和存储,并逐步实现连续时间以及多层空间的同步测定[57]。其中,热脉冲法、热平衡法、热扩散法经一系列的发展与完善,极大程度上减小了测量误差[58]。热扩散法还形成了成熟的商业化产品,并得到了广泛的应用。虽然利用外热比法和热场变形法测定活立木液流速率的研究有限,但外热比法实现了精确的零破坏检测,热场变形法液流速率测定范围广,并可准确的测定零液流和逆向液流。在利用外热比法测定液流速率时,需要针对不同样本以及实验条件设计不同的量规,这是外热比法的不足之处。因此,外热比法和热场变形法亟待更为深入研究。
目前,热技术方法成为液流测量的首要选择。在未来,应用热技术测定树干液流仍需关注以下热点:在完善研究活立木蒸腾耗水特点的同时,结合土壤生物因子和气象因子与树干液流的关系,进一步深入研究活立木生理作用;从微观和宏观方面监控水分运动,研究水分利用与树木生长的关系;在生产实际方面,进一步完善活立木单株和森林林区的数据监控,为实现高效的林区治理提供有力依据。
-
表 1 不同干扰情况下传播时间估计结果
Table 1. Results under there levels interference
干扰强度 采样点延迟(D) 传播时间(t传播)/s 低 1 059 0.021 2 中 1 052 0.021 0 高 1 053 0.021 0 -
[1] WANG Xiping,ROSS R J,MCCLELLAN M,et al. Nondestructive evaluation of standing trees with a stress wave method[J]. Wood Fiber Sci,2001,33(4):522-533. [2] 李光辉,尹建新. 基于应力波的木材无损检测技术研究进展[C]//中国自动化学会制造技术专门委员会. 第6届制造业自动化与信息化技术研讨会论文集. 桂林:[出版者不祥],2007:189-193. [3] 孟瑞华,胡勤龙,刘毅. 应力波时间仪的研制和应用[J]. 南京化工大学学报:自然科学版,2000,22(6):55-57. MENG Ruihua,HU Qinlong,LIU Yi. Research and design of stress wave time instrument and its applications[J]. J Nanjing Univ Chem Technol Nat Sci Ed,2000,22(6):55-57. [4] 冯国红,王立海,杨慧敏,等. 应力波木材无损检测信号采集系统[J]. 森林工程,2008,24(2):22-24. FENG Guohong,WANG Lihai,YANG Huimin,et al. Signal collection system of stress wave nondestructive testing of wood[J]. For Eng,2008,24(2):22-24. [5] 肖江,杨建华,李黎. 基于单片机的单板层积材抗弯强度应力波波速检测系统[J]. 木材加工机械,2008,19(1):1-5. XIAO Jiang,YANG Jianhua,LI Li. The detecting system of stress velocity for laminated veneer lumber bending strength based on single chip[J]. Wood Proc Mach,2008,19(1):1-5. [6] HAYES M P,CHEN J. A portable stress wave measurement system for timber inspection[C]//ENZC. Proceedings of the Electronics New Zealand Conference. Wellington:[s. n.],2003:1-6. [7] 何涛. 基于DSP的小波降噪技术研究及应用[D]. 郑州:河南工业大学,2011. HE Tao. Research and Application on Technology of Wavelet Denoising Based on DSP[D]. Zhengzhou:Henan University of Technology,2011. [8] 栾风虎,李玉峰,于学明,等. 基于高阶累计量的时延估计研究[J]. 黑龙江大学自然科学学报,2010,27(2):260-263. LUAN Fenghu,LI Yufeng,YU Xueming,et al. Research of time delay estimation based on higher-order cumulant[J]. J Nat Sci Heilongjiang Univ,2010,27(2):260-263. [9] KNAPPY C H,CARTER G C. Generalized correlation method or estimation of time delay[J]. IEEE Trans Acoust Speech Sign Proc,1976,24(4):320-327. [10] HERO A O,SCHWARTZ S C. New generalized cross correlator[J]. IEEE Trans Acoust Speech Sign Proc,1985,33(1):38-45. -
-
链接本文:
https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.2014.03.010