濕地松活立木彈性模量無(wú)損評(píng)估技術(shù)研究

張帥楠，姜景民，徐永勤，欒啟福*

(1.中國(guó)林業(yè)科學(xué)研究院亞熱帶林業(yè)研究所，浙江省林木育種技術(shù)研究重點(diǎn)實(shí)驗(yàn)室，浙江 杭州 311400；2.浙江省杭州市余杭區(qū)長(zhǎng)樂(lè)林場(chǎng)，浙江 杭州 311123)

濕地松活立木彈性模量無(wú)損評(píng)估技術(shù)研究

張帥楠1，姜景民1，徐永勤2，欒啟福1*

(1.中國(guó)林業(yè)科學(xué)研究院亞熱帶林業(yè)研究所，浙江省林木育種技術(shù)研究重點(diǎn)實(shí)驗(yàn)室，浙江 杭州 311400；2.浙江省杭州市余杭區(qū)長(zhǎng)樂(lè)林場(chǎng)，浙江 杭州 311123)

[目的]開(kāi)發(fā)操作簡(jiǎn)便、破壞性小且結(jié)果可靠的彈性模量無(wú)損評(píng)估技術(shù)，為實(shí)現(xiàn)以大群體試驗(yàn)材料為基礎(chǔ)的遺傳選擇和評(píng)估，培育速生優(yōu)質(zhì)林木良種奠定基礎(chǔ)。[方法]以22年生濕地松自由授粉家系測(cè)定林為研究對(duì)象，采用力學(xué)阻力儀(Pilodyn)和應(yīng)力波速(V)記錄儀(AV)等活立木無(wú)損評(píng)估儀器和方法以及木材物理力學(xué)常規(guī)檢測(cè)方法，建立活立木彈性模量評(píng)估技術(shù)體系。[結(jié)果]Pilodyn測(cè)定值Pr與實(shí)測(cè)基本密度(ρ)間呈顯著負(fù)相關(guān)(P<0.01)；彈性模量測(cè)定值MOEB與彈性模量(計(jì)算值)MOE(MOE=ρV2)之間呈極顯著線性相關(guān)關(guān)系，其線性回歸模型為MOEB=2.013 3+0.828 6MOE，決定系數(shù)R2=0.854 1；利用ρ和V計(jì)算的彈性模量值MOE與利用Pr和V計(jì)算的彈性模量相對(duì)值MOEP(MOEP=PrV2)間呈極顯著正相關(guān)關(guān)系(P<0.001)，回歸模型為MOE= 0.400 7+0.231 1MOEP, 決定系數(shù)R2=0.734 0；彈性模量MOE與V值相關(guān)極顯著(P<0.001)，V與MOE的回歸方程為MOE= -5.121 8+3.239 4V, 決定系數(shù)R2=0.921 9。[結(jié)論]在遺傳選擇時(shí)彈性模量(計(jì)算值)完全可以代表彈性模量測(cè)定值進(jìn)行評(píng)估與選育；可直接利用MOEP作為各樣木的彈性模量值參與遺傳分析；同時(shí)建立了MOE與V值之間的回歸方程，在比較同一測(cè)定群體彈性模量值相對(duì)大小時(shí)也可直接比較V值，大大減少了工作量。建立的模型對(duì)其它樹(shù)種測(cè)定具有參考意義，但尚需檢驗(yàn)。

濕地松；活立木；彈性模量；無(wú)損評(píng)估

為培育速生優(yōu)質(zhì)林木良種，需要大規(guī)模對(duì)遺傳測(cè)定群體活立木材性等重要品質(zhì)性狀開(kāi)展無(wú)損評(píng)估。其中，彈性模量(MOE)是林木建筑用材最重要的力學(xué)性質(zhì)，能夠反映木材的韌性和抵抗變形的能力，木材的彈性模量越大，在承受荷載時(shí)其變形越小[1-2]。目前，已成功應(yīng)用于活立木彈性模量無(wú)損評(píng)估的是基于應(yīng)力波速法的AV技術(shù)，它通過(guò)測(cè)量應(yīng)力波在開(kāi)始探頭和結(jié)束探頭之間的傳播時(shí)間來(lái)測(cè)得AV值(波速V)，該數(shù)值與立木強(qiáng)度具有顯著的相關(guān)性，可以作為活立木強(qiáng)度的檢測(cè)工具[3-5]。Ross等[6]采用應(yīng)力波檢測(cè)方法得出木材的彈性模量(MOE)、木材密度(ρ)和應(yīng)力波速度的關(guān)系式：MOE=ρV2，這為應(yīng)力波速技術(shù)在木材彈性模量無(wú)損檢測(cè)提供了研究基礎(chǔ)[7-9]?；盍⒛灸静幕久芏?ρ)測(cè)定的常規(guī)方法是用生長(zhǎng)錐在其胸徑處取木芯，然后利用飽和含水量法進(jìn)行測(cè)定[10]，雖然該方法獲得的數(shù)據(jù)直觀、可靠，但是成本高、耗時(shí)長(zhǎng)、具破壞性。Pilodyn是應(yīng)用最為普遍的一種基于阻力的活立木密度無(wú)損快速檢測(cè)儀器[11-12]，有效地克服了傳統(tǒng)檢測(cè)方法的缺點(diǎn)，在活立木木材密度評(píng)估中已有一定的進(jìn)展[13-17]。在實(shí)現(xiàn)了應(yīng)力波速(V)和基本密度(ρ)無(wú)損評(píng)估的基礎(chǔ)上，Chen等[5, 18]綜合利用了Pilodyn和AV技術(shù)，以Pilodyn探測(cè)值估測(cè)的木材外部密度作為木材基本密度，對(duì)挪威云杉(Piceaabies(L.) H. Karst.)大群體MOE和基本密度(ρ)進(jìn)行了估測(cè)，首次實(shí)現(xiàn)了數(shù)千株大群體活立木材性性狀的遺傳變異分析。

濕地松(PinuselliottiiEngelman var.elliottii)原產(chǎn)于美國(guó)東南部，現(xiàn)已成為我國(guó)南方重要的用材樹(shù)種[19-21]。經(jīng)過(guò)數(shù)十年的育種工作，濕地松生長(zhǎng)性狀獲得較大的改良[22-24]，以彈性模量為代表的力學(xué)性狀由于受測(cè)定技術(shù)影響改良進(jìn)程緩慢[25-28]。因此，開(kāi)發(fā)操作簡(jiǎn)單方便、破壞性小且結(jié)果可靠的彈性模量無(wú)損評(píng)估技術(shù)，實(shí)現(xiàn)以大群體試驗(yàn)材料為基礎(chǔ)的遺傳選擇和評(píng)估，對(duì)提高遺傳測(cè)定的精度與效率、促進(jìn)遺傳改良進(jìn)程具有重大現(xiàn)實(shí)意義。本研究采用力學(xué)阻力儀(Pilodyn)和應(yīng)力波速記錄儀(AV)等活立木無(wú)損評(píng)估儀器和方法以及木材物理力學(xué)常規(guī)檢測(cè)方法，初步建立濕地松活立木彈性模量評(píng)估技術(shù)體系，為濕地松彈性模量等材性性狀遺傳改良奠定基礎(chǔ)，同時(shí)對(duì)其它樹(shù)種材性測(cè)定具有參考意義。

1 試驗(yàn)林概況

試驗(yàn)地位于浙江省杭州市余杭區(qū)長(zhǎng)樂(lè)林場(chǎng)(30°27′ N，119°48′E)，亞熱帶濕潤(rùn)季風(fēng)性氣候，年均降水量1 399 mm，年均氣溫16.1℃，屬低山緩坡，土壤以紅壤為主。

試驗(yàn)林為濕地松子代測(cè)定林，營(yíng)建于1994年春季，含33個(gè)半同胞家系，采取單行6株小區(qū)、6個(gè)重復(fù)的隨機(jī)完全區(qū)組設(shè)計(jì)。

2 研究方法

2.1 材性性狀指標(biāo)測(cè)定方法

2.1.1 Pilodyn測(cè)定 2015年12月對(duì)上述33個(gè)濕地松半同胞家系進(jìn)行取樣測(cè)定，隨機(jī)選取6個(gè)重復(fù)中的3個(gè)重復(fù)，其中，每個(gè)小區(qū)選取1株平均木作為樣木，總計(jì)99個(gè)單株，進(jìn)行Pilodyn(6J，PROCEQ，Switzerland)單株測(cè)定[5, 29]。在樣木胸高南、北向分別測(cè)定2次(若2次探測(cè)值差值大于2 mm，則進(jìn)行第3次探測(cè)，選擇3次探測(cè)差異最小的2個(gè)值作為有效探測(cè)值)，取2次測(cè)定的平均值作為南向的Pilodyn值(Ps)和北向的Pilodyn值(Pn)，并取南、北向的平均值作為Pilodyn測(cè)定值(Pr)。

2.1.2 基本密度生長(zhǎng)錐測(cè)定 Pilodyn測(cè)定后，用內(nèi)徑為5 mm的生長(zhǎng)錐在樣木胸高處南北向鉆取木芯(木芯穿透樣木)?；久芏?ρ)采用飽和含水量法測(cè)定[10]，其計(jì)算公式如下：

(1)

式(1)中:ρ為基本密度(g·cm-3);Gmw為飽和含水率時(shí)試樣的質(zhì)量(g);Gh為絕干時(shí)試樣的質(zhì)量(g)。

2.1.4 彈性模量評(píng)估

2.1.4.1 無(wú)損評(píng)估 根據(jù)生長(zhǎng)錐法和AV技術(shù)得出的測(cè)定值ρ和V，應(yīng)用公式MOE=ρV2計(jì)算出各樣木的彈性模量(即彈性模量計(jì)算值，本文中以MOE表示)；根據(jù)Pilodyn法和AV技術(shù)得出的測(cè)定值Pr和V，以Pr代替ρ，應(yīng)用上述公式計(jì)算出各樣木的彈性模量相對(duì)值MOEP(MOEP=PrV2)。

2.1.4.2 常規(guī)測(cè)定 在99株樣木中隨機(jī)選取30株樣木依據(jù)國(guó)家標(biāo)準(zhǔn)《GB/T 1929-2009 木材物理力學(xué)試材鋸解及試樣截取方法》進(jìn)行破壞性取樣，并委托浙江省林科院質(zhì)檢站按照國(guó)家標(biāo)準(zhǔn)《GB/T 1936.

2-2009 木材抗彎彈性模量測(cè)定方法》測(cè)定其抗彎彈性模量(即彈性模量測(cè)定值，本文中以MOEB表示)。

2.2 統(tǒng)計(jì)分析方法

數(shù)據(jù)分析在R軟件(v 3.1.3，CRAN，http://cran.r-project.org)上進(jìn)行。

(1)各指標(biāo)間的相關(guān)性分析利用R軟件中的corr. test()函數(shù)計(jì)算Pearson相關(guān)系數(shù)r并進(jìn)行顯著性檢驗(yàn)。

(2)MOEB和MOE、MOE和MOEP以及MOE和V之間的回歸分析均基于普通最小二乘法(OLS)，利用R軟件中的lm()函數(shù)擬合其簡(jiǎn)單線性回歸模型：

MOEB= b1+k1MOE

(2)

MOE = b2+k2MOEP

(3)

(4)

式(2)～(4)中:b1、b2、b3代表回歸常數(shù)，k1、k2、k3代表回歸系數(shù)。

3 結(jié)果與分析

3.1 各測(cè)定指標(biāo)的基本情況

99個(gè)濕地松單株材性性狀各測(cè)定指標(biāo)的基本情況見(jiàn)表1。濕地松木材基本密度和彈性模量分別為0.502 0 g·cm-3、4.495 6 GPa，其中，MOEP極差與方差均較大，說(shuō)明彈性模量值變異豐富，這可能是由于MOEP計(jì)算時(shí)因子Pr與V2的乘積引起的MOEP值波動(dòng)較大，可能會(huì)有較好的遺傳選擇成效。基本密度性狀測(cè)定指標(biāo)(ρ、Pr)變異系數(shù)較小(CV<10%)，彈性模量性狀測(cè)定指標(biāo)(MOE、MOEP)變異系數(shù)較大(CV>20%)，也說(shuō)明濕地松彈性模量性狀相對(duì)于基本密度性狀具有更豐富的遺傳變異性。

表1 濕地松基本密度和彈性模量性狀各測(cè)定指標(biāo)的基本情況Table 1 Basic information of determination of wood density and MOE of slash pine

3.2 各測(cè)定指標(biāo)間相關(guān)性分析

對(duì)濕地松木材基本密度性狀測(cè)定指標(biāo)(ρ、Pr)和彈性模量性狀測(cè)定指標(biāo)(V、MOE、MOEP)進(jìn)行相關(guān)分析，結(jié)果見(jiàn)表2。從表2中可以看出:ρ與Pr呈顯著負(fù)相關(guān)關(guān)系(P<0.01)，而這二者基本密度性狀測(cè)定指標(biāo)與3個(gè)彈性模量性狀測(cè)定指標(biāo)之間相關(guān)關(guān)系不顯著或顯著性低(P<0.1)；V與MOE及MOEP達(dá)到極顯著的高度正相關(guān)(r>0.90，P<0.001)。這說(shuō)明波速是影響濕地松彈性模量估測(cè)的關(guān)鍵因子，木材密度對(duì)彈性模量影響相對(duì)較小,而MOE與MOEP之間也呈極顯著高度正相關(guān)(r=0.85，P<0.001)，因此，可以認(rèn)為MOEP也可以對(duì)MOE進(jìn)行預(yù)測(cè)和評(píng)估。

表2 濕地松基本密度和彈性模量性狀各測(cè)定 指標(biāo)間相關(guān)關(guān)系Table 2 The correlations of determination of wood density and MOE indexes of slash pine

注：括號(hào)中為該相關(guān)系數(shù)顯著性水平；*代表顯著性水平為0.1，**代表顯著性水平為0.01，***代表顯著性水平為0.001。

Note: Significance probability of the correlations were given in parenthesis; *, **, *** representP< 0.05,P< 0.01 andP< 0.001, respectively.

3.3 彈性模量測(cè)定值和彈性模量(計(jì)算值)的回歸方程

濕地松彈性模量測(cè)定值與彈性模量(計(jì)算值)的線性回歸模型為MOEB=2.013 3+0.828 6MOE，其決定系數(shù)R2=0.854 1(圖1)。模型中，回歸常數(shù)和回歸系數(shù)都達(dá)極顯著水平(P<0.001)，即MOEB與MOE呈極顯著線性相關(guān)關(guān)系。因此，利用關(guān)系式MOE=ρV2得出的濕地松彈性模量計(jì)算值MOE雖然不能作為其標(biāo)準(zhǔn)值，卻完全可以代表彈性模量測(cè)定值進(jìn)行評(píng)估與選育。

圖1 濕地松彈性模量測(cè)定值與彈性模量(計(jì)算值)的
簡(jiǎn)單線性回歸模型
Fig.1 Simple linear regression model between MOEBand MOE

3.4 彈性模量(計(jì)算值)和彈性模量相對(duì)值的回歸方程

濕地松彈性模量與彈性模量相對(duì)值的線性回歸模型為MOE=0.400 7+0.231 1MOEP，其決定系數(shù)R2為0.734 0(圖2)。模型中，回歸常數(shù)(P<0.05)和回歸系數(shù)(P<0.001)都達(dá)顯著水平，這說(shuō)明該線性模型適合于MOE與MOEP之間的回歸分析，也就是說(shuō)利用Pilodyn法和AV技術(shù)得出的測(cè)定值Pr和V，計(jì)算出MOEP并應(yīng)用以上模型即能快速評(píng)估彈性模量。

圖2 濕地松彈性模量值與彈性模量相對(duì)值的簡(jiǎn)單
線性回歸模型
Fig.2 Simple linear regression model between MOE and MOEP

3.5 彈性模量(計(jì)算值)和波速的回歸方程

濕地松彈性模量與波速的線性回歸模型為MOE=-5.121 8+3.239 4V，其決定系數(shù)R2為0.921 9(圖3)。模型中，回歸常數(shù)和回歸系數(shù)也都達(dá)極顯著水平(P<0.001)，該模型的方差解釋率為91.19%，這說(shuō)明該線性模型能很好地?cái)M合MOE與V之間的線性回歸關(guān)系，即僅利用AV技術(shù)也能實(shí)現(xiàn)彈性模量評(píng)估。

圖3 濕地松彈性模量值與波速V的簡(jiǎn)單線性回歸模型
Fig.3 Simple linear regression model between MOE and V

4 討論

彈性模量與樹(shù)木本身抵抗倒伏能力密切相關(guān)[31]。濕地松作為當(dāng)前我國(guó)南方最重要的用材、采脂及荒山綠化優(yōu)良樹(shù)種，種植面積超過(guò)200萬(wàn)hm2，然而，由于前期研究工作對(duì)其抗風(fēng)雪等力學(xué)性質(zhì)的選擇力度不夠，部分地區(qū)濕地松遇到異常天氣情況出現(xiàn)雪壓、風(fēng)倒等災(zāi)害，成為制約該樹(shù)種發(fā)展的一個(gè)重要方面[32]，因此，從樹(shù)種健康角度考慮，需要繼續(xù)加強(qiáng)其彈性模量等物理力學(xué)性質(zhì)的評(píng)估與選育研究。隨著無(wú)損評(píng)估技術(shù)的發(fā)展，全林選擇測(cè)定濕地松彈性模量成為可能。本研究以力學(xué)阻力儀(Pilodyn)和應(yīng)力波速記錄儀(AV)等活立木無(wú)損評(píng)估儀器和方法，初步建立濕地松活立木彈性模量評(píng)估技術(shù)體系，為濕地松彈性模量等材性性狀遺傳改良奠定基礎(chǔ)，同時(shí)對(duì)其它樹(shù)種測(cè)定具有參考意義。

自Ross等[6]研究得出關(guān)系式MOE=ρV2以來(lái)，基于應(yīng)力波速法的AV技術(shù)在雪松(Cedrusdeodara(Roxb.) G. Don.)[7]、歐洲赤松(PinussylvestrisL.)[8]、火炬松(PinustaedaL.)[33]、杉木(Cunninghamialanceolata(Lamb.) Hook.)[9]等樹(shù)種活立木MOE評(píng)估中的研究應(yīng)用逐漸成熟，為濕地松大群體活立木彈性模量的快速評(píng)估奠定了基礎(chǔ)。本研究結(jié)果也顯示出濕地松彈性模量測(cè)定值與彈性模量(計(jì)算值)之間存在著極顯著的線性相關(guān)關(guān)系，其線性回歸模型為MOEB=2.013 3+0.828 6MOE，決定系數(shù)R2=0.854 1?？梢钥闯?，在濕地松彈性模量評(píng)估中，關(guān)系式MOE=ρV2不能得出其標(biāo)準(zhǔn)彈性模量值，但這并不能否定該關(guān)系式的應(yīng)用價(jià)值，在濕地松無(wú)損、快速評(píng)估及其遺傳改良研究進(jìn)程中，其仍具有相當(dāng)大的優(yōu)勢(shì)，完全可以代表彈性模量測(cè)定值進(jìn)行評(píng)估與選育。

濕地松木材基本密度性狀測(cè)定指標(biāo)和彈性模量性狀測(cè)定指標(biāo)相關(guān)分析結(jié)果顯示：基本密度ρ與Pilodyn測(cè)定值Pr呈顯著負(fù)相關(guān)關(guān)系。直接以Pr值代替ρ計(jì)算的彈性模量相對(duì)值MOEP與彈性模量MOE之間呈極顯著的高度正相關(guān)關(guān)系，回歸模型MOE=0.400 7+0.231 1MOEP，因此，在不要求精確獲取各樣木彈性模量絕對(duì)值時(shí)，可直接利用MOEP作為各樣木的彈性模量值參與遺傳分析。以Pr值代替ρ計(jì)算的彈性模量相對(duì)值MOEP較之先前彈性模量無(wú)損評(píng)估研究而言，節(jié)省了測(cè)定相應(yīng)木材基本密度的工作量[8-9]，具有一定的便利優(yōu)勢(shì)，可在現(xiàn)實(shí)工作中應(yīng)用。同時(shí)，相關(guān)性研究中還發(fā)現(xiàn),基本密度性狀測(cè)定指標(biāo)(ρ、Pr)與彈性模量性狀測(cè)定指標(biāo)(MOE、MOEP)之間相關(guān)關(guān)系不顯著或顯著性較低，波速V與彈性模量MOE以及彈性模量相對(duì)值MOEP則均達(dá)到極顯著高度正相關(guān)，說(shuō)明波速V是彈性模量間接估測(cè)中的關(guān)鍵因子，而木材密度對(duì)其影響相對(duì)較小。彈性模量MOE與波速V之間的線性回歸模型MOE= -5.121 8+3.239 4V，該模型很好地?cái)M合了本研究中MOE與V之間的線性關(guān)系?？赡苁茉囼?yàn)材料影響，Wu等[34]以22個(gè)4年生桉樹(shù)(Eucalyptusspp.)無(wú)性系為研究對(duì)象，結(jié)果認(rèn)為Pilodyn與彈性模量呈顯著負(fù)相關(guān)，而本研究中結(jié)果顯示，盡管這二者之間存在負(fù)相關(guān)關(guān)系，但相關(guān)是不顯著的，與Chen等[5]對(duì)5 618株挪威云杉材性無(wú)損測(cè)定的結(jié)果相同。因此，木材基本密度與彈性模量之間的相關(guān)性受試驗(yàn)對(duì)象或環(huán)境條件影響較大，不同研究對(duì)象評(píng)估時(shí)需要建立對(duì)應(yīng)的評(píng)估模型。

一般認(rèn)為彈性模量由物體基本密度(ρ)和應(yīng)力波速(V)共同決定，由經(jīng)驗(yàn)方程MOE=ρV2可以看出V的影響應(yīng)該更大，本研究中僅應(yīng)用V就能很好地預(yù)測(cè)彈性模量，可以將公式中ρ看作是一個(gè)常數(shù)，這和Kennedy等[9]的觀點(diǎn)一致，為活立木彈性模量的簡(jiǎn)單快速預(yù)測(cè)奠定了基礎(chǔ)，同時(shí)研究中直接用Pilodyn測(cè)定值Pr代替ρ計(jì)算的彈性模量估計(jì)值MOEP也可以作為活立木彈性模量的參考值。

5 結(jié)論

(1)濕地松彈性模量測(cè)定值MOEB與彈性模量(計(jì)算值)MOE之間存在極顯著線性相關(guān)關(guān)系，在遺傳選擇時(shí)MOE完全可以代表彈性模量測(cè)定值進(jìn)行評(píng)估與選育。

(2)彈性模量相對(duì)值MOEP與彈性模量MOE之間呈極顯著高度正相關(guān)關(guān)系，回歸模型MOE=0.400 7+0.231 1MOEP，可直接利用MOEP作為各樣木的彈性模量值參與遺傳分析。

(3)彈性模量MOE與波速V之間的線性回歸模型MOE=-5.121 8+3.239 4V，在比較同一測(cè)定群體彈性模量值相對(duì)大小時(shí)也可直接比較V值，大大減少了工作量。

以上建立的模型對(duì)其它樹(shù)種測(cè)定具有參考意義，但尚需檢驗(yàn)。

[1] Raymond C A, Kube P D, Pinkard L,etal. Evaluation of non-destructive methods of measuring growth stress inEucalyptusglobulus: relationships between strain, wood properties and stress[J]. Forest Ecology and Management, 2004, 190(2): 187-200.

[2] Lasserre J P, Mason E G, Watt M S. Assessing corewood acoustic velocity and modulus of elasticity with two impact based instruments in 11-year-old trees from a clonal-spacing experiment ofPinusradiataD. Don.[J]. Forest Ecology and Management, 2007, 239(1): 217-221.

[3] Grabianowski M, Manley B, Walker J C F. Acoustic measurements on standing trees, logs and green lumber[J]. Wood Sci Technol, 2006, 40(3): 205-216.

[4] Wessels C B, Malan F S, Rypstra T. A review of measurement methods used on standing trees for the prediction of some mechanical properties of timber[J]. European Journal of Forest Research, 2011, 130(6): 881-893.

[5] Chen Z Q, Karlsson B, Lundqvist S O,etal. Estimating solid wood properties using Pilodyn and acoustic velocity on standing trees of Norway spruce[J]. Annals of Forest Science, 2015, 72(4): 499-508.

[6] Ross R J, Pellerin R F. NDE of wood-based composites with longitudinal stress wave[J]. Forest Products Journal, 1988, 38(5):39-45.

[7] Nanami N, Nakamura N, Arima T,etal. Measuring the properties of standing trees with stress waves, 3: Evaluating the properties of standing trees for some forest stands[J]. Journal of the Japan Wood Research Society, 1993, 39(8):903-909.

[8] Auty D, Achim A. The relationship between standing tree acoustic assessment and timber quality in Scots pine and the practical implications for assessing timber quality from naturally regenerated stands[J]. Forestry, 2008, 81(4): 475-487.

[9] Kennedy S G, Cameron A D, Lee S J. Genetic relationships between wood quality traits and diameter growth of juvenile core wood in Sitka spruce[J]. Can J For Res, 2013, 43(1): 1-6.

[10] 成俊卿. 木材學(xué)[M]. 北京: 中國(guó)林業(yè)出版社，1985: 282-464.

[11] Yanchuk A D, and Kiss G K. Genetic variation in growth and wood specific gravity and its utility in the improvement of interior spruce in British Columbia[J]. Silvae Genet, 1993,42(2-3): 141-148.

[12] Koch L, Fins L. Genetic variation in wood specific gravity from progeny tests of ponderosa pine (PinusponderosaLaws.) in northern Idaho and Western Montana[J]. Silvae Genet, 2000, 49(4-5): 174-181.

[13] Cown D J. Comparison of the Pilodyn and torsiometer methods for the rapid assessment of wood density in living trees[J]. N Z J For Sci, 1978, 8(3): 384-391.

[14] Taylor F W. Rapid determination of southern pine specific gravity with a Pilodyn tester[J]. For Sci,1981, 27(1): 59-61.

[15] 朱景樂(lè)，王軍輝，張守攻，等. 毛白楊材性指標(biāo)預(yù)測(cè)及選擇[J]. 林業(yè)科學(xué)，2008，44(7): 23-28.

[16] 欒啟福，盧 萍，井振華，等. Pilodyn評(píng)估雜交松活立木的基本密度及其性狀相關(guān)分析[J]. 江西農(nóng)業(yè)大學(xué)學(xué)報(bào)，2011，33(3): 548-552.

[17] 易 敏，賴 猛，孫曉梅，等. 無(wú)損檢測(cè)在日本落葉松材性育種中的應(yīng)用前景探討[J]. 林業(yè)科學(xué)，2014，50(11): 96-103.

[18] Chen Z Q, Gil M R G, Karlsson B,etal. Inheritance of growth and solid wood quality traits in a large Norway spruce population tested at two locations in southern Sweden[J]. Tree Genetics & Genomes, 2014, 10(5): 1291-1303.

[19] 曲桂林，吳志民. 國(guó)外林業(yè)技術(shù)考察報(bào)告選編[M]. 北京: 中國(guó)林業(yè)出版社，1996: 63-72.

[20] Wen X, Kuang Y, Shi M,etal. Biology of Hylobitelus xiaoi (Coleoptera: Curculionidae), a new pest of slash pine, Pinus elliottii[J]. Journal of Economic Entomology, 2004, 97(6): 1958-1964.

[21] Barnett J P, Sheffield R M. Slash pine: characteristics, history, status and trends[C]// Dickens E D, Barnett J P, Hubbard W G,etal. Slash pine: still growing and growing! Proceedings of the slash pine symposium, Jekyll Island, Georgia, USA, 23-25 April 2002. Gen. Tech. Rep. SRS-76. Asheville, NC: U.S. Department of Agriculture, Forest Service, Southern Research Station, 2005:1-6.

[22] 潘志剛，鄭勇奇. 濕地松種源試驗(yàn)研究[J]. 林業(yè)科學(xué)研究，1990，3(4): 314-322.

[23] White T L, Hodge G R, Powell G L. An advanced-generation tree improvement plan for slash pine in the southeastern United States[J]. Silvae Genetica, 1993(42):359-359.

[24] 姜景民，孫海菁，劉昭息. 濕地松自由授粉家系的早期評(píng)定[J]. 林業(yè)科學(xué)研究，1995，8(5):574-581.

[25] Zobel, B J. Inheritance of wood properties in conifers[J]. Silvae Genet, 1961,10(3): 65-96.

[26] Wei X, Borralho N M G. Genetic control of basic density and bark thickness and their relationships with growth traits ofEucalyptusurophyllain south east China[J]. Silvae Genetica, 1997, 46(4):245-250.

[27] Hansen C P. Application of the Pilodyn in forest tree improvement. DFSC Series of Technical Notes. TN55[R]//Humlebaek, Denmark: DANIDA Forest Seed Centre, 2000.

[28] 姜笑梅，駱秀琴，殷亞方，等. 不同濕地松種源木材材性遺傳變異的研究[J]. 林業(yè)科學(xué)，2002，38(3):130-135.

[29] 朱景樂(lè)，王軍輝，張守攻，等. Pilodyn在林木遺傳改良應(yīng)用中的研究進(jìn)展[J]. 浙江林學(xué)院學(xué)報(bào)，2008，25(5): 661-665.

[30] Carter P, Briggs D, Ross R J,etal. Acoustic testing to enhance western forest values and meet customer wood quality needs[R] //Harrington C A, Schoenholtz S H. Productivity of western forests: a forest products focus. General technical report PNW-GTR-642. Portland: U.S. Department of Agriculture, Forest Service, Pacific Northwest Research, 2005: 121-129.

[31] 許秀玉，王明懷，仲崇祿，等. 不同樹(shù)種木材性質(zhì)及其抗臺(tái)風(fēng)性能[J]. 浙江農(nóng)林大學(xué)學(xué)報(bào)，2014，31(5): 751-757.

[32] 欒啟福，盧 萍，肖復(fù)明，等. 雨雪冰凍重災(zāi)區(qū)濕地松受害情況調(diào)查及其原因初步分析[J]. 林業(yè)科學(xué)，2008，44(11): 50-54.

[33] Mora C R, Schimleck L R, Isik F,etal. Relationships between acoustic variables and different measures of stiffness in standing Pinus taeda trees[J]. Canadian Journal of Forest Research, 2009, 39(8): 1421-1429.

[34] Wu S, Xu J, Li G,etal. Use of the Pilodyn for assessing wood properties in standing trees ofEucalyptusclones[J]. Journal of Forestry Research, 2010, 21(1): 68-72.

(責(zé)任編輯：金立新)

Study on the Modulus of Elasticity Non-destructive Evaluation Technique of Slash Pine Standing Tree

ZHANGShuai-nan1,JIANGJing-min1,XUYong-qin2,LUANQi-fu1

(1.Research Institute of Subtropical Forestry, Chinese Academy of Forestry, Key Laboratory of Tree Breeding of Zhejiang Province,Hangzhou 311400, Zhejiang, China; 2.Changle Forestry Farm of Yuhang District, Hangzhou, Hangzhou 311123, Zhejiang, China)

[Objective]To develop the non-destructive evaluation technique of modulus of elasticity (MOE) for realizing the genetic selection and evaluation of forest tree with the traits of fast growing and high quality based on a large population of test materials. [Method]The evaluation systems of MOE of standing tree were established by non-destructive evaluation technique, Pilodyn and acoustic velocity (AV), and conventional testing methods based on the open pollinated family of 22-year-old slash pine (Pinuselliottii).[Result]The Pilodyn value (Pr) and the basic density (ρ) showed a significant negative correlation (P<0.01). The measurement value ofMOE(MOEB), determined by conventional methods, had a very significant linear correlation with theMOEvalue calculated by the formulaMOE=ρV2. The regression model wasMOEB=2.013 3+0.828 6MOE(Coefficient of determinationR2=0.854 1). A significant positive correlation (P<0.001) was observed betweenMOEandMOEP, the relative value ofMOE, (MOEP=PrV2) (Vmeans the value ofAV), and the linear regression model wasMOE= 0.400 7+0.231 1MOEP(coefficient of determinationR2=0.734 0). The positive correlation betweenMOEandVwas extremely significant (P<0.001) and the linear regression model betweenMOEandVwasMOE= -5.121 8+3.239 4V(coefficient of determinationR2=0.921 9).[Conclusion]The modulus of elasticity calculated byMOE=ρV2can represent the determination values ofMOE. TheMOEPcan be used directly as theMOEvalue for each sample tree to participate in genetic analysis. It is feasible to compare only theVwhen evaluating the relative value ofMOEof the sample trees in the same population. Therefore, these linear models are also useful for other tree species.

slash pine; standing tree; modulus of elasticity; non-destructive evaluation

10.13275/j.cnki.lykxyj.2017.01.011

2016-05-12

國(guó)家自然科學(xué)基金項(xiàng)目(31570668)；“十二·五”國(guó)家林業(yè)科技支撐計(jì)劃專題(2012BAD01B0203)；國(guó)家林業(yè)局“948”引進(jìn)項(xiàng)目(2012-4-46)

張帥楠(1993—)，男，河南安陽(yáng)人，碩士研究生，主要從事林木遺傳育種研究. * 通訊作者：欒啟福，副研究員，博士，從事林木遺傳育種研究. E-mail: qifu.luan@caf.ac.cn

S791.246

A

1001-1498(2017)01-0075-06