王乙 高忠燕 田國(guó)雙

(東北林業(yè)大學(xué)，哈爾濱，150040) (扎龍國(guó)家級(jí)自然保護(hù)區(qū)管理局) (東北林業(yè)大學(xué))

野生動(dòng)物生態(tài)游憩價(jià)值，是在不干預(yù)野生動(dòng)物自然生活狀態(tài)、不消耗野生動(dòng)物資源、不破壞自然生態(tài)環(huán)境的前提下，人類參與野生動(dòng)物觀賞、旅游、休閑活動(dòng)產(chǎn)生的價(jià)值[1]。野生動(dòng)物生態(tài)游憩，是積極的野生動(dòng)物保護(hù)手段，也是發(fā)展區(qū)域經(jīng)濟(jì)的有效途徑。據(jù)統(tǒng)計(jì)，2016年國(guó)內(nèi)旅游44.4億人次，比上年增長(zhǎng)11.%[2]。研究野生動(dòng)物生態(tài)游憩價(jià)值，找出提高生態(tài)游憩價(jià)值的方法，對(duì)科學(xué)保護(hù)野生動(dòng)物、促進(jìn)旅游產(chǎn)業(yè)健康發(fā)展、提高人民福祉有重要意義。

目前對(duì)野生動(dòng)物生態(tài)游憩價(jià)值的評(píng)價(jià)，多采用旅行費(fèi)用支出法[3]、條件價(jià)值評(píng)估法[4]等方法評(píng)價(jià)。這些研究多以貨幣化結(jié)果為主，較少關(guān)注哪些因素對(duì)價(jià)值構(gòu)成影響；同時(shí)，這些方法來(lái)自于環(huán)境價(jià)值評(píng)價(jià)，用于研究野生動(dòng)物生態(tài)游憩價(jià)值，評(píng)價(jià)結(jié)果容易受到外在環(huán)境因素影響。選擇實(shí)驗(yàn)法(CE)是一種用于自然資源價(jià)值評(píng)價(jià)較為先進(jìn)的方法[5-10]。雖然國(guó)內(nèi)外學(xué)者利用CE對(duì)濕地、森林、海岸、耕地、景區(qū)等自然資源開展大量研究工作[11-15]，但應(yīng)用于野生動(dòng)物生態(tài)游憩價(jià)值的研究案例相對(duì)少見。筆者以扎龍自然保護(hù)區(qū)丹頂鶴(Grusjaponensis)為例，利用CE對(duì)野生動(dòng)物生態(tài)游憩價(jià)值進(jìn)行評(píng)估，旨在為保護(hù)區(qū)生態(tài)游憩功能建設(shè)、CE在野生動(dòng)物生態(tài)游憩價(jià)值評(píng)價(jià)中的應(yīng)用提供參考。

1 研究區(qū)自然概況

研究區(qū)域扎龍國(guó)家級(jí)自然保護(hù)區(qū)，位于黑龍江省齊齊哈爾市東南部，松嫩平原的西部。地處齊齊哈爾市鐵鋒區(qū)、昂昂溪區(qū)、泰來(lái)縣、富?？h和大慶市杜爾伯特蒙古族自治縣、林甸縣的交界處。地理坐標(biāo)為東經(jīng)123°47′～124°37′，北緯46°52′～47°32′；呈東北至西南不規(guī)則的長(zhǎng)方形，南北長(zhǎng)80.6 km，東西寬58.0 km，總面積21萬(wàn)hm2。扎龍自然保護(hù)區(qū)屬大陸性半干旱季風(fēng)氣候，年均降水量402.7 mm，年均氣溫3.9 ℃。保護(hù)區(qū)內(nèi)土壤主要為黑鈣土、草甸土和沼澤土。

該區(qū)是丹頂鶴等珍稀水禽棲息地。現(xiàn)有鳥類16目48科269種，其中國(guó)家重點(diǎn)保護(hù)的鳥類有35種，特別是有丹頂鶴、白枕鶴、白頭鶴、白鶴、灰鶴、蓑羽鶴6種。另外，有白鸛、黑鸛、天鵝以及眾多的雁鴨、鷗等鳥類棲息在這里。還有2種爬行類、6種兩棲類、21種獸類、46種魚類動(dòng)物，277種昆蟲和468種高等植物。該區(qū)具有較好的生態(tài)游憩功能。

2 研究方法

2.1 選擇實(shí)驗(yàn)法

選擇實(shí)驗(yàn)法(CE)的理論基礎(chǔ)，是消費(fèi)者理論和隨機(jī)效用理論。消費(fèi)者對(duì)野生動(dòng)物生態(tài)游憩的支付意愿受到多個(gè)狀態(tài)屬性影響，任何一個(gè)屬性的變化會(huì)導(dǎo)致消費(fèi)者的選擇發(fā)生變化。由于受訪者每個(gè)選擇總是以效用最大化作為目的，因此，可以通過(guò)觀察消費(fèi)者支付意愿分析野生動(dòng)物生態(tài)游憩價(jià)值；也可以通過(guò)觀察屬性變化對(duì)受訪者選擇的影響，分析提高哪種屬性水平，能夠最大程度的提高野生動(dòng)物生態(tài)游憩價(jià)值。

CE的調(diào)查問(wèn)卷設(shè)計(jì)核心是設(shè)計(jì)選擇集。選擇集，是由與某一野生動(dòng)物物種相關(guān)的不同狀態(tài)屬性組合而成，每個(gè)選擇集有一個(gè)參考方案和多個(gè)備選方案；通常以野生動(dòng)物的現(xiàn)狀作為參考方案；每個(gè)方案中都包含一個(gè)用貨幣量表達(dá)的屬性，同時(shí)含有若干個(gè)野生動(dòng)物狀態(tài)屬性；被調(diào)查者以不同選擇集中的貨幣量屬性，衡量狀態(tài)屬性進(jìn)行選擇。

CE中的效用函數(shù)，可由產(chǎn)品某些可識(shí)別屬性的固定項(xiàng)與那些無(wú)法觀測(cè)的影響因素共同決定。假設(shè)，Uij表示第i個(gè)人選擇j方案的效用；Vij表示第i個(gè)人選擇j方案的固定項(xiàng)；εij表示隨機(jī)擾動(dòng)項(xiàng)，即表示無(wú)法觀測(cè)的影響因素。而固定項(xiàng)又可設(shè)定為一個(gè)線性組合Vij=βiXij，式中Xij表示第i個(gè)人選擇j方案的屬性特征、βi表示備選屬性特征集Xij的系數(shù)向量。

對(duì)面臨J種選擇的第i個(gè)調(diào)查者，假定該調(diào)查者選擇j方案的效用為Uij=βiXij+εij；i=1、2、…、N，j=1、2、…、J。因每個(gè)被調(diào)查者都會(huì)選擇對(duì)自己效用最大的方案，如果被調(diào)查者i選擇了k方案，則有理由相信對(duì)調(diào)查者i而言，一定是方案選擇集C的J個(gè)備選方案中，第k個(gè)方案的效用最大。因此，用統(tǒng)計(jì)模型表示調(diào)查者i選擇k方案的概率為：P(Uik>Uij)=P[(Uik-Uij)>(εij-εik)]，?j≠k。

當(dāng)J個(gè)隨機(jī)干擾項(xiàng)之間相互獨(dú)立同分布，且服從極限分布時(shí)的分布函數(shù)為P(εij≤t)=F(t)=exp{-exp(t)}。于是得選擇實(shí)驗(yàn)中最常用的條件Logit模型：

式中：μ為標(biāo)量參數(shù)，通常取值1。

模型參數(shù)估計(jì)采用極大似然估計(jì)方法，其對(duì)數(shù)似然函數(shù)為：

式中：yij為指示變量。表示當(dāng)?shù)趇個(gè)人選擇第j個(gè)備選方案時(shí)，為1，否則為0；對(duì)每個(gè)i，有且僅有1個(gè)yij=1。

當(dāng)?shù)玫絽?shù)向量的估計(jì)值β，即可以測(cè)算因?qū)傩宰兓鴮?dǎo)致的福利變化值。即由初始狀態(tài)X0到替代狀態(tài)X1的福利變化為：

式中：α為收入福利值，通常用成本屬性的系數(shù)表示。

由此可以進(jìn)一步通過(guò)屬性系數(shù)與成本屬性系數(shù)之比，計(jì)算由野生動(dòng)物狀態(tài)屬性水平邊際變化所導(dǎo)致的價(jià)值變化(vi)vi=-βi/βM，βi為野生動(dòng)物各狀態(tài)屬性項(xiàng)的估計(jì)系數(shù)；βM為收入的福值，通常用成本項(xiàng)的估計(jì)系數(shù)表示。

2.2 調(diào)查問(wèn)卷設(shè)計(jì)

于2017年4月設(shè)計(jì)調(diào)查問(wèn)卷，并通過(guò)保護(hù)區(qū)工作人員進(jìn)行預(yù)調(diào)研，取得數(shù)據(jù)并對(duì)問(wèn)卷進(jìn)行調(diào)整，形成最終問(wèn)卷。5月開始正式問(wèn)卷調(diào)研。

調(diào)查問(wèn)卷分為兩個(gè)部分。第一部分是游客基礎(chǔ)信息部分，包括游客的性別、年齡、學(xué)歷、工作性質(zhì)、出游目的、總體花費(fèi)以及出游滿意度。第二部分為CE調(diào)查主體問(wèn)卷部分，圍繞“丹頂鶴數(shù)量”、“丹頂鶴年齡結(jié)構(gòu)”、“丹頂鶴展示方式”、“游客游覽方式”、“游客支出增加意愿”5個(gè)狀態(tài)屬性設(shè)計(jì)調(diào)查問(wèn)卷(見表1)。問(wèn)卷針對(duì)丹頂鶴生態(tài)游憩屬性設(shè)置，觀察旅行者對(duì)丹頂鶴生態(tài)游憩不同屬性提升方案帶來(lái)的支出意愿變化情況。

根據(jù)CE屬性設(shè)置水平，本調(diào)查問(wèn)卷共有405種(3×3×3×3×5)不同水平的組合方式，使用SPSS軟件進(jìn)行正交試驗(yàn)，選取具有代表性的組合方式，最終得到25個(gè)具有代表性的水平組合。將當(dāng)前的各項(xiàng)屬性組合作為基準(zhǔn)方案，剩余24個(gè)組合進(jìn)行隨機(jī)排序，然后按順序選取2個(gè)組合與基準(zhǔn)方案一同作為一個(gè)選擇集，得到12個(gè)選擇集，每個(gè)選擇集3個(gè)觀測(cè)值。選擇集1見表2。

注：*表示當(dāng)前水平。

表2 選擇集1的設(shè)置

2.3 數(shù)據(jù)收集

調(diào)查分為兩個(gè)階段。第一階段為預(yù)調(diào)研，針對(duì)游客來(lái)源、總花費(fèi)、游覽方式等屬性進(jìn)行調(diào)查并調(diào)整問(wèn)卷設(shè)計(jì)。第二階段為實(shí)地調(diào)查，發(fā)放問(wèn)卷300組，回收率100%，通過(guò)問(wèn)卷篩查，有效問(wèn)卷91%。

本研究共取得的273份有效問(wèn)卷，其中男性154人、女性119人；25～55歲游客占90%，其中25～35歲游客比例最高，學(xué)生群體占游客總數(shù)的49.45%。月收入選項(xiàng)，1 000～1 500元占據(jù)游客總數(shù)的41.03%、>1 500～3 000元占總數(shù)的35.16%，說(shuō)明中等收入游客是目前扎龍自然保護(hù)區(qū)的主要游客。出游目的方面，出于觀賞丹頂鶴為目的的游客占總游客的93.77%，此外1.47%的游客出于科學(xué)研究的目的來(lái)到扎龍自然保護(hù)區(qū)(見表3)。

扎龍自然保護(hù)區(qū)距離中心城市較遠(yuǎn)，因此其游憩價(jià)值應(yīng)包括旅行費(fèi)用和消費(fèi)者剩余。調(diào)研中發(fā)現(xiàn)旅行總花費(fèi)方面，200～400元區(qū)間占據(jù)消費(fèi)者支出的絕大比例，這與游客的經(jīng)濟(jì)屬性和出游城市相關(guān)。本次調(diào)查取平均值312元作為人均消費(fèi)金額(見表4)。

3 結(jié)果與分析

3.1 條件Logit模型分析

調(diào)研回收有效問(wèn)卷273份，每份問(wèn)卷有6個(gè)觀測(cè)值，共有1 638個(gè)觀測(cè)值。采用Stata軟件對(duì)樣本進(jìn)行條件Logit模型參數(shù)估計(jì)(見表5)。由表5可見：除A1、T1變量外，其余變量均在95%置信水平顯著，模型的擬合度指標(biāo)R2=0.194。

表3 扎龍自然保護(hù)區(qū)游客人口統(tǒng)計(jì)學(xué)特征

表4 扎龍自然保護(hù)區(qū)人均消費(fèi)描述統(tǒng)計(jì)

表5 樣本條件Logit模型估計(jì)結(jié)果

Logit估計(jì)結(jié)果中系數(shù)的大小，反應(yīng)各要素水平下對(duì)扎龍自然保護(hù)區(qū)丹頂鶴生態(tài)游憩價(jià)值的相對(duì)重要性。模型分析中，無(wú)法測(cè)算基礎(chǔ)場(chǎng)景的系數(shù)，但可以假設(shè)同種屬性在不同水平下的效用之和為零，因此各屬性系數(shù)之和的負(fù)數(shù)即為基準(zhǔn)屬性的系數(shù)。在所有屬性中，丹頂鶴數(shù)量(N)屬性是對(duì)游客效用影響最大的屬性。N1(1.001)是本研究屬性中游客最支持的狀態(tài)，N3(-1.232)是游客最不支持的狀態(tài)。這說(shuō)明，丹頂鶴數(shù)量的屬性變化，最能影響游客效用。通過(guò)系數(shù)可發(fā)現(xiàn)，丹頂鶴展示方式(M)，對(duì)游客效用影響處于第二位。在展示方式中，M3(-0.846)是游客最不滿意的狀態(tài)，M1(0.532)是游客在該屬性當(dāng)中最滿意的狀態(tài)。說(shuō)明游客對(duì)保護(hù)區(qū)內(nèi)丹頂鶴展示的自然程度具有很高的要求，丹頂鶴生活狀態(tài)越自然，游客越滿意，能給游客帶來(lái)的效用值也越高。游客游覽方式(T)屬性系數(shù)中，T1(0.317)能夠?yàn)橛慰蛶?lái)相對(duì)最好的效用，說(shuō)明游客對(duì)游覽方式多樣化的需求，特別是對(duì)開展個(gè)性化游覽方式的偏好。在所有屬性當(dāng)中，丹頂鶴年齡(A)的相對(duì)重要性最低，尤其是從當(dāng)前狀態(tài)A2(0.145)提高到最好水平A1(0.274)，給游客帶來(lái)的效用增幅最小。說(shuō)明隨著游客可見丹頂鶴年齡層級(jí)的豐富，游客效用會(huì)有提升，但提升幅度最小。

由于非線性模型的系數(shù)，不能直接解釋各屬性對(duì)游客決策的邊際影響，因此需要利用機(jī)會(huì)比進(jìn)行解釋和分析。機(jī)會(huì)比，可以表示各屬性相互作用下，由某一屬性水平而引起的被選擇相對(duì)概率變化值。從選擇樣本看，N、A、M、T屬性的改善，都會(huì)引起被選概率的增加，其中N升高到N1水平(100只以上)，被選擇的概率增加了2.713倍；而如果降為N3(10只以內(nèi))的水平，被選概率僅為原有概率的0.287倍；表明，游客對(duì)丹頂鶴數(shù)量水平非常敏感。同時(shí)，丹頂鶴展示方式(M)屬性，如果能夠提升到M1水平，游客的選擇概率增加1.852倍；T屬性，如果能夠提升到T1水平，游客的選擇概率增加1.346倍。

3.2 游客對(duì)丹頂鶴生態(tài)游憩價(jià)值的支付意愿

丹頂鶴屬性水平的變化會(huì)帶來(lái)額外支付意愿的變化。當(dāng)前屬性水平，游客的額外支付意愿為107.26元/人；當(dāng)所有屬性水平增加到最佳時(shí)，游客額外支付意愿為265.51元/人；當(dāng)所有屬性降低到最差時(shí)，游客的支付意愿為負(fù)數(shù)，即在現(xiàn)有支付水平下減少支付372元。在所有屬性水平中，從N2狀態(tài)提升為N1狀態(tài)，可額外增加95.26元的支付意愿，是所有屬性水平提高中支付意愿提高最高的屬性(見表6)。

表6 各屬性狀態(tài)下游客支付意愿

丹頂鶴當(dāng)前狀態(tài)水平，游客額外支付意愿107.26元/人，游客的額外支付意愿又稱為隱含價(jià)格，是現(xiàn)有支付價(jià)格基礎(chǔ)上的隱含價(jià)值。扎龍自然保護(hù)區(qū)人均消費(fèi)312元/人，在計(jì)算總游憩價(jià)值時(shí)需同時(shí)納入計(jì)算。因此可以得出，當(dāng)前丹頂鶴游客總支付意愿約為419元/人。

4 結(jié)論與討論

本文以扎龍自然保護(hù)區(qū)丹頂鶴為例，運(yùn)用選擇實(shí)驗(yàn)法(CE)評(píng)價(jià)其生態(tài)游憩價(jià)值。研究表明，丹頂鶴作為扎龍自然保護(hù)區(qū)保護(hù)旗艦物種，具有巨大的生態(tài)游憩價(jià)值。當(dāng)前狀態(tài)下，游人額外支付意愿為107.26元/人，如果各項(xiàng)屬性水平達(dá)到最佳狀態(tài)，游客的額外支付意愿為265.51元/人。當(dāng)前，扎龍自然保護(hù)區(qū)年接待游客約17萬(wàn)人，以人均實(shí)際消費(fèi)平均值312.56元/人計(jì)算，當(dāng)前扎龍丹頂鶴生態(tài)游憩價(jià)值約7 136萬(wàn)元/a；如各屬性水平達(dá)到最佳，扎龍丹頂鶴游憩價(jià)值約9 674萬(wàn)元/a。

當(dāng)前狀態(tài)下，游客對(duì)丹頂鶴游憩屬性的選擇與額外支付意愿是一致的。選擇偏好由高至低依次為丹頂鶴展示方式、丹頂鶴數(shù)量、游客游覽方式、丹頂鶴年齡結(jié)構(gòu)，額外支付意愿為39.25、28.88、21.00、18.13元/人。游客對(duì)丹頂鶴展示方式的偏好程度最高，且愿意為此多支付39.25元/人；滿意程度最低的是丹頂鶴年齡結(jié)構(gòu)，支付意愿也相對(duì)較低，為8.13元/人；這符合扎龍自然保護(hù)區(qū)丹頂鶴游憩活動(dòng)實(shí)際情況。扎龍自然保護(hù)區(qū)以保護(hù)丹頂鶴種群及其生境為主要目標(biāo)，游人可以近距離觀賞野生丹頂鶴、觀看丹頂鶴放飛，人們能夠感受到的樂(lè)趣處于滿意水平之上。

游客對(duì)丹頂鶴數(shù)量屬性的提高愿望最強(qiáng)烈。一方面，為了使游憩過(guò)程中丹頂鶴數(shù)量屬性達(dá)到研究的最佳狀態(tài)，游客愿意額外支付125.13元/人，相當(dāng)于在當(dāng)前人均實(shí)際支付水平下提高40.1%；另一方面，一旦游覽過(guò)程中丹頂鶴數(shù)量下降到最差水平，游客支付意愿將降低154元/人，相當(dāng)于在當(dāng)前人均實(shí)際支付水平下降低49.35%。丹頂鶴數(shù)量受保護(hù)區(qū)生態(tài)環(huán)境影響最明顯，因此，必須加強(qiáng)保護(hù)區(qū)生態(tài)管理，努力提高生態(tài)撫育水平，這與當(dāng)前保護(hù)區(qū)設(shè)置的基本目的和建設(shè)目標(biāo)是一致的。

選擇實(shí)驗(yàn)法(CE)在野生動(dòng)物生態(tài)游憩價(jià)值評(píng)估中具有良好的適用性，評(píng)價(jià)結(jié)果具有較高的可信度。本研究選取“丹頂鶴數(shù)量”、“丹頂鶴年齡結(jié)構(gòu)”、“丹頂鶴展示方式”、“游客游覽方式”、“游客支出增加意愿”5個(gè)屬性進(jìn)行討論，這些屬性代表生態(tài)旅游過(guò)程中野生動(dòng)物的數(shù)量、野生動(dòng)物種群結(jié)構(gòu)、游客感受和游覽路徑選擇。這些屬性符合選擇實(shí)驗(yàn)法對(duì)野生動(dòng)物游憩價(jià)值的評(píng)價(jià)要求，可以推廣至其他野生動(dòng)物物種的生態(tài)游憩價(jià)值評(píng)價(jià)。

在問(wèn)卷調(diào)查的基礎(chǔ)上獲取數(shù)據(jù)，仍然難以避免傳統(tǒng)問(wèn)卷帶來(lái)的諸多偏差，如何找到更加客觀的數(shù)據(jù)來(lái)源，避免受訪者主觀因素對(duì)研究結(jié)果的影響，仍值得深入探討。

[1] 張飛.基于SWOT分析的納木錯(cuò)國(guó)家公園生態(tài)旅游開發(fā)對(duì)策研究[J].中南林業(yè)科技大學(xué)學(xué)報(bào)(社會(huì)科學(xué)版)，2017，11(4)：71-74.

[2] 中國(guó)旅游研究院.2016年全年旅游統(tǒng)計(jì)數(shù)據(jù)報(bào)告及2017年旅游經(jīng)濟(jì)形勢(shì)預(yù)測(cè).[EB/OL].(2017-02-24)[2017-10-15]http://www.ctaweb.org/html/2017-2/2017-2-24-9-10-97907.html.

[3] 施德群.基于旅行費(fèi)用法的觀鳥游憩價(jià)值評(píng)估[D].北京：北京林業(yè)大學(xué)，2005.

[4] 周學(xué)紅，馬建章，張偉.我國(guó)東北虎保護(hù)的經(jīng)濟(jì)價(jià)值評(píng)估：以哈爾濱市居民的支付意愿研究為例[J].東北林業(yè)大學(xué)學(xué)報(bào)，2007，35(5)：81-83，86.

[5] 樊輝，趙敏娟.自然資源非市場(chǎng)價(jià)值評(píng)估的選擇實(shí)驗(yàn)法：原理及應(yīng)用分析[J].資源科學(xué)，2013，35(7)：1347-1354.

[6] 王文智，武拉平.選擇實(shí)驗(yàn)理論及其在食品需求研究中的應(yīng)用：文獻(xiàn)綜述[J].技術(shù)經(jīng)濟(jì)，2014，33(1)：110-117.

[7] JETTE BREDAHL J， JOHN HALFDAN B， BO JELLESMARK T， et al. What’s in a name? The use of quantitative measures versus ‘Iconised’ species when valuing biodiversity[J]. Environmental and Resource Economics，2008，39(3)：247-263.

[8] HORNE P， BOXALL P， ADAMOWICZ W. Multiple-use management of forest recreation sites： a spatially explicit choice experiment[J]. Forest Ecology and Management，2005，207(1/2)：189-199.

[9] KATARIA M， BATEMAN I， CHRISTENSEN T， et al. Scenario realism and welfare estimates in choice experiments-A non-market valuation study on the European water framework directive[J]. Journal of Environmental Management，2012，94(1)：25-33.

[10] DIAS N， BELCHER K. Value and provision of ecosystem services from prairie wetlands： A choice experiment approach[J]. Ecosystem Services，2015，15：35-44.

[11] 史恒通，趙敏娟.基于選擇試驗(yàn)?zāi)Ｐ偷纳鷳B(tài)系統(tǒng)服務(wù)支付意愿差異及全價(jià)值評(píng)估：以渭河流域?yàn)槔齕J].資源科學(xué)，2015，37(2)：351-359.

[12] 呂歡歡.基于選擇實(shí)驗(yàn)法的國(guó)家森林公園游憩資源價(jià)值評(píng)價(jià)研究[D].大連：大連理工大學(xué)，2013.

[13] 蘇紅巖，李京梅.基于改進(jìn)選擇實(shí)驗(yàn)法的廣西紅樹林濕地修復(fù)意愿評(píng)估[J].資源科學(xué)，2016，38(9)：1810-1819.

[14] 馬愛慧，蔡銀鶯，張安錄.基于選擇實(shí)驗(yàn)法的耕地生態(tài)補(bǔ)償額度測(cè)算[J].自然資源學(xué)報(bào)，2012，27(7)：1154-1163.

[15] 王晉楠.基于選擇實(shí)驗(yàn)法的村鎮(zhèn)景觀游憩價(jià)值評(píng)估[D].北京：北京林業(yè)大學(xué)，2014.