基于主成分分析的橡膠種植適宜性評(píng)價(jià)

蘇文地　張培松　羅微

摘 要 選取土壤及氣候資源，根據(jù)前人的研究成果及專家經(jīng)驗(yàn)選擇確定評(píng)價(jià)因子，運(yùn)用模糊數(shù)學(xué)方法來量化評(píng)價(jià)因子，采用主成分分析方法對(duì)海南省儋州市的橡膠種植適宜性進(jìn)行評(píng)價(jià)。結(jié)果表明：儋州市最適宜橡膠種植的面積為43 029.58 hm2，占農(nóng)用地面積的14.99%；適宜面積為57 770.90 hm2，占20.12%，次適宜面積為40 934.9 hm2，占14.26%。該評(píng)價(jià)結(jié)果可以從宏觀上實(shí)現(xiàn)儋州市的橡膠收益最大化以及為橡膠種植的區(qū)劃提供決策參考。

關(guān)鍵詞 橡膠樹 ；適宜性 ；種植 ；主成分分析

分類號(hào) S794.1

Abstract This paper selects the soil and climate resources， according to the previous research results and the expertise to determine the evaluation factor. Use the fuzzy mathematical method to quantify the evaluation factors， and then evaluated the planting suitability of rubber in Danzhou City by the method of principal component analysis. The results showed that the most suitable area for rubber planting is 43 029.58 ha， accounting for 14.99 % of the total agricultural land in Danzhou City； the suitable area for planting is 57 770.90 ha， accounting for 20.12 %； the less suitable area for planting is 40 934.9 ha， accounting for 14.26 %. The evaluation results not only bring the maximum benefit of the rubber at the macro， but also provide a reference for the zoning of the rubber plantations in Danzhou City.

Keywords rubber tree ； suitability ； planting ； principal component analysis

橡膠是關(guān)乎國計(jì)民生和國防安全的重要戰(zhàn)略物資，是四大工業(yè)原料中唯一可再生的資源[1]。海南省是我國最大的天然橡膠生產(chǎn)基地，橡膠產(chǎn)業(yè)不僅是海南農(nóng)民的重要經(jīng)濟(jì)來源，也是海南省的基礎(chǔ)產(chǎn)業(yè)、支柱產(chǎn)業(yè)和優(yōu)勢(shì)產(chǎn)業(yè)[2]。隨著我國人民的生活水平日益提高以及經(jīng)濟(jì)全球化和區(qū)域經(jīng)濟(jì)一體化，人們對(duì)天然橡膠的需求隨之增加。面臨國際市場(chǎng)的競(jìng)爭(zhēng)，如何在有限的土地上將橡膠的效益最大化，且最大限度地保護(hù)土壤資源、挖掘農(nóng)業(yè)資源的潛力，已成為橡膠生產(chǎn)管理部門迫切需要解決的問題[3]。雖然近年來海南橡膠種植的發(fā)展勢(shì)頭良好，但是在橡膠種植生產(chǎn)中還存在著一些問題：有些農(nóng)戶仍然過著靠天吃飯的日子或者憑著經(jīng)驗(yàn)從事橡膠生產(chǎn)，有些農(nóng)戶由于文化素質(zhì)較低而大量的施肥以片面地追求產(chǎn)量。這些生產(chǎn)模式并非是最好的，還忽略了生態(tài)效益；再有儋州市的橡膠實(shí)際種植面積及產(chǎn)量位居海南?。ㄞr(nóng)墾除外）第一，但是平均產(chǎn)量卻位居全省第六[4]。為了實(shí)現(xiàn)橡膠產(chǎn)量的提高，且在不增加種植面積的情況下，增加橡膠的總產(chǎn)量，這就要求逐步提高橡膠的單產(chǎn)。因此，因地制宜地種植橡膠已成為海南橡膠種植管理的熱點(diǎn)。

鑒于上述原因，本研究針對(duì)海南省儋州市的橡膠種植進(jìn)行了適宜性分析評(píng)價(jià)，以提高橡膠的單產(chǎn)，實(shí)現(xiàn)橡膠生產(chǎn)的經(jīng)濟(jì)效益最大化，增加農(nóng)民的收入。通過橡膠種植的適宜性評(píng)價(jià)為橡膠種植區(qū)劃作出理論依據(jù)，為農(nóng)業(yè)管理者提供科學(xué)準(zhǔn)確的區(qū)劃參考。

1 材料與方法

1.1 材料

儋州市位于海南島的西北部，是海南省土地面積最大、人口最多的市（縣級(jí)市），它瀕臨北部灣，地處東經(jīng)108°56′～109°46′，北緯19°11′～19°52′，是海南西部的經(jīng)濟(jì)、交通、通信和文化中心。儋州市陸地面積3 265 km2，人口106.86萬（2010年底）[4]。由于地處東亞大陸季風(fēng)氣候的南緣，屬熱帶季風(fēng)氣候。年平均氣溫23.3℃，年平均日照時(shí)數(shù)在2 000 h以上，西部沿海日照達(dá)2 500 h左右；降雨量充沛：全市各地降雨量為900～2 500 mm/a，平均降雨量為1 800 mm/a；全市全年雨量分布不均，冬春季節(jié)雨量稀少，夏秋季節(jié)雨量充沛。全市地表年徑流量達(dá)到189億m3。全市全年水熱資源豐富，適宜多種農(nóng)作物生長(zhǎng)繁殖。2010年全市熱帶作物面積46 646 hm2，其中橡膠種植面積45 605.07 hm2[4]。

空間數(shù)據(jù)：1∶50 000儋州市行政區(qū)界圖、土壤類型圖、土地利用現(xiàn)狀圖、數(shù)字高程模型等；

文本數(shù)據(jù)：儋州市氣象統(tǒng)計(jì)數(shù)據(jù)、儋州市“測(cè)土配方施肥”項(xiàng)目的采樣調(diào)查分析結(jié)果，包括采樣點(diǎn)坐標(biāo)、基本情況、土壤農(nóng)化分析數(shù)據(jù)，以及儋州市其他相關(guān)資料與數(shù)據(jù)。

1.2 數(shù)據(jù)處理方法

數(shù)據(jù)處理軟件：SPSS 統(tǒng)計(jì)分析軟件、GIS 地理信息系統(tǒng)軟件。

2 結(jié)果與分析

2.1 基于主成分方法的橡膠種植適宜性評(píng)價(jià)

2.1.1 評(píng)價(jià)單元的劃分

評(píng)價(jià)單元是橡膠種植適宜性評(píng)價(jià)的最基本單元，本研究的評(píng)價(jià)單元?jiǎng)澐植捎猛寥缊D、土地利用現(xiàn)狀圖、行政區(qū)劃圖疊置構(gòu)成。相同的土壤類型、土地利用現(xiàn)狀及行政區(qū)劃組成一個(gè)評(píng)價(jià)單元，以便于評(píng)價(jià)結(jié)果的應(yīng)用。其中土壤類型劃分到土屬，土地利用現(xiàn)狀提取農(nóng)用地，行政區(qū)劃劃分到行政村。通過疊加分析處理后得到儋州市的11 761個(gè)評(píng)價(jià)單元。endprint

2.1.2 評(píng)價(jià)指標(biāo)的選取及量化

參評(píng)的指標(biāo)主要分為兩類：一類是連續(xù)型指標(biāo)（即定量化的指標(biāo)），包括土壤養(yǎng)分值、土層厚度、容重、風(fēng)速、降雨量等；一類是離散型指標(biāo)（即定性、概念型指標(biāo)），包括質(zhì)地，成土母質(zhì)、剖面構(gòu)型等。由于參評(píng)指標(biāo)不能夠簡(jiǎn)潔地表明該因素對(duì)橡膠適宜性影響的大小，所以本研究采用“隸屬度”來界定[5]。

2.1.2.1 隸屬函數(shù)

由于參評(píng)的指標(biāo)是漸變的，固本文采用模糊數(shù)學(xué)的方法來建立隸屬函數(shù)。根據(jù)前人的研究經(jīng)驗(yàn)和數(shù)學(xué)原理來建立函數(shù)，為了避免計(jì)算時(shí)零值過多，隸屬函數(shù)定為0.1～1.0，最大值1.0表示該指標(biāo)完全適合作物生長(zhǎng)，最低值0.1表示該指標(biāo)不利于作物生長(zhǎng)。建立相應(yīng)的函數(shù)模型[6-8]：

2.1.2.2 指標(biāo)的選擇以及指標(biāo)臨界值的確定

根據(jù)橡膠正常生長(zhǎng)對(duì)環(huán)境條件的相關(guān)研究，結(jié)合專家經(jīng)驗(yàn)法和理論分析，確定各評(píng)價(jià)因子的隸屬函數(shù)的臨界值和閾值及散點(diǎn)型函數(shù)的隸屬度[9-13]。如表1。

2.1.3 評(píng)價(jià)因子屬性數(shù)據(jù)庫建立

利用儋州市各氣象站點(diǎn)的空間位置（緯度、經(jīng)度）及其觀測(cè)的數(shù)據(jù)，建立風(fēng)速、年降水量數(shù)據(jù)庫。采用反距離權(quán)重插值模型對(duì)氣候要素?cái)?shù)據(jù)進(jìn)行模擬空間處理，建立風(fēng)速、年降水量等柵格空間數(shù)據(jù)庫。利用儋州市DEM模型，借助GIS 相關(guān)模塊自動(dòng)生成坡度、坡向、海拔高度等地形因子?xùn)鸥窨臻g及其屬性數(shù)據(jù)庫。利用土壤普查、測(cè)土配方施肥調(diào)查樣點(diǎn)的資料和DEM模型，建立調(diào)查樣點(diǎn)及其相關(guān)屬性數(shù)據(jù)庫。采用克里格插值法生成儋州市土壤耕層有機(jī)質(zhì)、有效磷、速效鉀等柵格數(shù)據(jù)庫；利用ArcGIS進(jìn)行區(qū)域統(tǒng)計(jì)，作為評(píng)價(jià)單元相關(guān)因子的屬性數(shù)據(jù)。

2.1.4 主成分分析的計(jì)算步驟

（1）各量化指標(biāo)數(shù)據(jù)標(biāo)準(zhǔn)化

將采集p維隨機(jī)向量x=（x1，x2，...，xp）T）n個(gè)樣品

xi=（xi1，xi2，...，xipT ，i=1，2，…，n，n>p，構(gòu)造成樣本陣，接著對(duì)樣本陣元進(jìn)行標(biāo)準(zhǔn)化轉(zhuǎn)換；如下：zij=，i=1，2，…，n；j=1，2，…，p；其中=≥， 得到標(biāo)準(zhǔn)化陣Z。

（2）按照公式R=|Tij|pxp=；其中，rij=，i，j=1，2，…，p；求標(biāo)準(zhǔn)化陣Z的相關(guān)系數(shù)，并構(gòu)造相關(guān)系數(shù)矩陣R。

（3）計(jì)算出樣本的相關(guān)系數(shù)矩陣R的特征方程|R-λIp|=0得m個(gè)特征根，確定主成分的個(gè)數(shù)。按≥0.85確定所需的m值，使信息的利用率達(dá)85%以上，對(duì)每個(gè)λj，j=1，2，...，m，解方程組Rb=λjb得單位特征向量b0j。

（4）根據(jù)Uij=ZTiB0j，j=1，2，...，m；其中：U1為第一主成分；U2為第二主成分；…，Up為第p主成分。計(jì)算各個(gè)主成分的得分。

（5）對(duì)求得的m個(gè)主成分分別進(jìn)行評(píng)價(jià)，得到m個(gè)主成分的評(píng)價(jià)分值。再將每個(gè)主成分的評(píng)價(jià)分值與每個(gè)主成分的方差貢獻(xiàn)率（即權(quán)重）相乘累加，計(jì)算評(píng)價(jià)的綜合得分。

2.1.5 評(píng)價(jià)結(jié)果

根據(jù)指標(biāo)的臨界值和閾值及散點(diǎn)型函數(shù)的隸屬度，將指標(biāo)全部通過隸屬函數(shù)量化于同一水平上。然后利用數(shù)據(jù)分析軟件SPSS16.0和Excel完成主成分分析的相關(guān)計(jì)算。得到橡膠主成分分析對(duì)應(yīng)的特征根、貢獻(xiàn)率及其累計(jì)方差貢獻(xiàn)率（表2）。

從表2可以看出，當(dāng)m=8時(shí)，累計(jì)方差貢獻(xiàn)率85.207%>85%，因此可以利用前8個(gè)主成分代替原變量進(jìn)行分析[14]。從表2可知，橡膠的8個(gè)主成分與全部的14個(gè)指標(biāo)的線性關(guān)系公式為：

根據(jù)評(píng)價(jià)函數(shù)模型，將評(píng)價(jià)單元數(shù)據(jù)代入上述公式中，計(jì)算每個(gè)評(píng)價(jià)單元的綜合F值。將所有評(píng)價(jià)單元的綜合得分F作累積頻率曲線，并把頻率曲線的拐點(diǎn)作為相鄰等級(jí)的臨界值[15]。以最適宜、適宜、次適宜和不適宜對(duì)橡膠適宜性綜合指數(shù)進(jìn)行歸類，從而對(duì)適宜性進(jìn)行定性評(píng)價(jià)。橡膠的適宜性劃分標(biāo)準(zhǔn)見表3。

在橡膠種植適宜性綜合指數(shù)劃分基礎(chǔ)上，遵循種植區(qū)集中連片且顧及行政區(qū)劃界線等原則，結(jié)合海南限制橡膠生長(zhǎng)的自然環(huán)境條件：高程高于300 m，風(fēng)速高于3 m/s[10]等。將評(píng)價(jià)單元進(jìn)行適宜性定性評(píng)價(jià)，并在GIS中，將評(píng)價(jià)單元圖與定性評(píng)價(jià)結(jié)果相銜接，最后得到儋州市橡膠適宜性面積統(tǒng)計(jì)表（表4）及適應(yīng)性評(píng)價(jià)結(jié)果圖（圖1、2）。

2.1.6 結(jié)果分析

從圖1可以看出，儋州市中部區(qū)域橡膠種植適宜性較高，從中間往東南、西北兩邊呈適宜遞減趨勢(shì)。橡膠種植適宜區(qū)域主要集中在那大鎮(zhèn)附近，不適宜區(qū)主要集中在沿海及松濤水庫附近。

結(jié)合表4和圖1，可以清楚的看到儋州市橡膠種植的適宜區(qū)域：

（1）最適宜區(qū)域。橡膠種植最適宜區(qū)域的面積為43 029.58 hm2，占農(nóng)用地面積的14.99%。此類區(qū)域分布較為集中，主要分布在大成、和慶、那大、南豐四鎮(zhèn)以及八一、蘭洋、西聯(lián)、西流、西培、西慶等農(nóng)場(chǎng)。此區(qū)域農(nóng)用地的主要特點(diǎn)是地勢(shì)平坦、土層較厚、風(fēng)速較小、降雨量充足、土壤酸堿度適宜且養(yǎng)分充足。對(duì)于此類農(nóng)用地，應(yīng)當(dāng)切實(shí)加強(qiáng)保護(hù)，盡量避免建設(shè)用地等占用此類耕地，同時(shí)注意做好養(yǎng)地，合理調(diào)整化肥施用量和投入比，防止土壤酸化和鹽漬化。

（2）適宜區(qū)域。橡膠種植適宜區(qū)域的面積為57 770.9 hm2，占農(nóng)用地面積的20.13%。此類耕地主要分布在大成、和慶、那大、雅星等鄉(xiāng)鎮(zhèn)以及西華、西聯(lián)、西流、西培、西慶、新盈等農(nóng)場(chǎng)。此類農(nóng)用地環(huán)境與最適宜區(qū)域農(nóng)用地相差不大，影響橡膠種植的不利因素主要是水和養(yǎng)分，尤其是有機(jī)質(zhì)和全氮等養(yǎng)分的含量不高，處于中等或中等偏下水平。由于適宜種植面積較最適宜面積廣，因此加強(qiáng)此類耕地的管理和保護(hù)是儋州市實(shí)現(xiàn)橡膠高產(chǎn)、穩(wěn)產(chǎn)的保證。針對(duì)適宜耕地養(yǎng)分特別是有效磷含量低的特點(diǎn)，可因時(shí)適量增加磷肥用量，合理調(diào)配復(fù)合肥比例，以提高土壤有效磷含量，改善土壤養(yǎng)分結(jié)構(gòu)，達(dá)到以磷增氮的效果，消除農(nóng)業(yè)生產(chǎn)中的不利因素。endprint

（3）次適宜區(qū)域。次適宜種植橡膠的面積為40 934.9 hm2，占14.26%；該類農(nóng)用地集中分布在雅星、王五、東城、光村、峨蔓、木棠等鄉(xiāng)鎮(zhèn)以及紅嶺、西華、西聯(lián)、新盈等農(nóng)場(chǎng)。其影響橡膠種植的不利因素主要是降水量較小、風(fēng)速較大、土壤養(yǎng)分含量不高，雖然其地勢(shì)較低，但最重要的是此類土壤有機(jī)質(zhì)和全氮等養(yǎng)分含量低，難以維持橡膠正常生長(zhǎng)的生理需要。對(duì)于此類區(qū)域可以選擇抗風(fēng)性較強(qiáng)的橡膠品種種植，并且加強(qiáng)農(nóng)田水利設(shè)施的建設(shè)，調(diào)節(jié)水分狀況，增施有機(jī)肥等，使之滿足橡膠生長(zhǎng)需要。

（4）不適宜區(qū)域。不適宜種植橡膠的面積為145 317.6 hm2，占農(nóng)用地面積的50.62%。該類農(nóng)用地面積較大，主要分布在沿海以及松濤水庫一帶。沿海地區(qū)由于降雨量較小、風(fēng)速較大等氣候條件以及養(yǎng)分含量低等條件限制而無法種植；在松濤水庫雖然土壤養(yǎng)分充足，但由于其地勢(shì)較高、坡度較大等條件限制橡膠的種植，且缺少管理，這些都使得松濤水庫附近的大片區(qū)域不適宜種植橡膠。

（5）由于八一農(nóng)場(chǎng)在儋州市橡膠生產(chǎn)中占有較大比重，且農(nóng)場(chǎng)在規(guī)劃、管理上較為精細(xì)，因此以八一農(nóng)場(chǎng)來檢驗(yàn)橡膠適宜性評(píng)價(jià)結(jié)果。根據(jù)資料數(shù)據(jù)[4]計(jì)算得到儋州市的橡膠單產(chǎn)為2.5 kg/年·株。結(jié)合八一農(nóng)場(chǎng)的實(shí)際情況，將八一農(nóng)場(chǎng)的橡膠產(chǎn)量進(jìn)行分等定級(jí)，將干膠大于3.0 kg/年·株的地塊確定為最適宜種植，適宜種植為2.5～3.0 kg/年·株，次宜種植為2.0～2.5 kg/年·株，小于2.0 kg/年·株的為不適宜種植。將確定好的橡膠產(chǎn)量水平和評(píng)價(jià)結(jié)果進(jìn)行對(duì)應(yīng)性比較，得到圖2。

（6）從圖2中不難看出，對(duì)應(yīng)“相同”的區(qū)域（17 855.8 hm2）明顯多于“不相同”區(qū)域（2 454.9 hm2），“相同”區(qū)域占驗(yàn)證區(qū)域的87.9%。對(duì)應(yīng)性結(jié)果表明：圖1的評(píng)價(jià)結(jié)果對(duì)于八一農(nóng)場(chǎng)橡膠的種植具有較高的適用性。

3 討論與結(jié)論

（1）種植適宜性評(píng)價(jià)涉及土地眾多性狀，即使所選的因子相同，采用不同的評(píng)價(jià)方法，其評(píng)價(jià)結(jié)果也有所不同。本研究從儋州市出發(fā)，以村一級(jí)為研究尺度，采用模糊數(shù)學(xué)的方法，根據(jù)橡膠樹正常生長(zhǎng)對(duì)環(huán)境條件的相關(guān)研究來確定各評(píng)價(jià)因子的隸屬函數(shù)的臨界值和閾值及散點(diǎn)型函數(shù)的隸屬度。并運(yùn)用主成分分析法，對(duì)橡膠種植適宜性進(jìn)行綜合評(píng)價(jià)，并制作了適宜性種植圖。通過適宜性評(píng)價(jià)，因地制宜的發(fā)展橡膠樹種植產(chǎn)業(yè)，從而提高橡膠的產(chǎn)出率，進(jìn)而創(chuàng)造出最大的經(jīng)濟(jì)、生態(tài)和社會(huì)效益。

此方法在實(shí)際應(yīng)用中仍有一些問題需要注意。首先，數(shù)據(jù)資料的獲取是本方法應(yīng)用的制約因素，一些重要的評(píng)價(jià)因子數(shù)據(jù)將在很大程度上影響評(píng)價(jià)結(jié)果的科學(xué)性；其次，評(píng)價(jià)因子的確定具有一定的局限性，因?yàn)樵u(píng)價(jià)因子及隸屬函數(shù)是根據(jù)數(shù)據(jù)資料以及專家意見加以平衡得出的，受影響較大。另外， 在評(píng)價(jià)過程中獲得的基礎(chǔ)數(shù)據(jù)越多、真實(shí)性越高， 所得出的評(píng)價(jià)單元將越零碎。出于評(píng)價(jià)目的的考慮， 為了更好地為農(nóng)業(yè)管理部門提供決策參考，對(duì)于種植適宜性評(píng)價(jià)結(jié)果， 在評(píng)價(jià)單元（地塊）上，還有待進(jìn)一步的研究。

（2）在今后的研究中，建立、完善一套科學(xué)、可行的適宜性評(píng)價(jià)方法仍有很多工作需要去做。同時(shí)，還應(yīng)加強(qiáng)人員的管理以及肥料的合理使用，以確保橡膠穩(wěn)產(chǎn)、高產(chǎn)。

參考文獻(xiàn)

[1] 許海平，傅國華. 我國天然橡膠產(chǎn)業(yè)發(fā)展趨勢(shì)[J]. 中國熱帶農(nóng)業(yè)，2007，（2）：15-16.

[2] 許道順，許升鋒. 經(jīng)濟(jì)一體化背景下海南省橡膠產(chǎn)業(yè)的發(fā)展戰(zhàn)略[J]. 海南金融，2006，（9）：25-28.

[3] 黃昌勇. 土壤學(xué)[M]. 北京：中國農(nóng)業(yè)出版社，2000.

[4] 海南省統(tǒng)計(jì)局. 海南統(tǒng)計(jì)年鑒2011[M]. 中國統(tǒng)計(jì)出版社，2011.

[5] 朱國宇，黃川友，華國春. 層次分析法在水環(huán)境規(guī)劃中的應(yīng)用[J]. 東北水利水電，2003，21（4）：1-7.

[6] 曹志洪，周建民. 中國土壤質(zhì)量[M]. 北京：科學(xué)出版社，2008.

[7] 劉紹貴，張?zhí)伊郑跖d祥，等. 江西省余江縣土壤肥力變化和驅(qū)動(dòng)因素與對(duì)策研究[J]. 土壤通報(bào)，2006，（10）：869-874.

[8] 張 靜. 作物-地域多種組合中作物生態(tài)適宜性評(píng)價(jià)與權(quán)重配置方法的研究[D]. 南京：南京農(nóng)業(yè)大學(xué)，2005.

[9] 華南亞熱帶作物研究所. 中國橡膠栽培學(xué)[M]. 北京：中國科學(xué)出版社，1961.

[10] 何 康，黃宗道. 熱帶北緣橡膠所栽培[M]. 廣州：廣東科技出版社，1987.

[11] 王秉忠. 橡膠栽培學(xué)（第三版）[M]. 華南熱帶農(nóng)業(yè)大學(xué)出版社，2000.

[12] 韋娉婷. 耕地地力評(píng)價(jià)方法在橡膠園中的應(yīng)用研究——以海南農(nóng)墾八一片區(qū)為例[D]. ?？冢汉Ｄ洗髮W(xué)，2009.

[13] 歐 濱，羅 微，馬 利，等. 基于GIS的橡膠園地分等研究[J]. 熱帶農(nóng)業(yè)科學(xué)，2009，29（6）：9-14.

[14] 付 強(qiáng). 數(shù)據(jù)處理方法及其農(nóng)業(yè)應(yīng)用[M]. 北京：科學(xué)出版社，2006.

[15] 韓永學(xué). 特爾斐法與“拿來主義”[J]. 哈爾濱師專學(xué)報(bào)，2000（2）：67-68.endprint

（3）次適宜區(qū)域。次適宜種植橡膠的面積為40 934.9 hm2，占14.26%；該類農(nóng)用地集中分布在雅星、王五、東城、光村、峨蔓、木棠等鄉(xiāng)鎮(zhèn)以及紅嶺、西華、西聯(lián)、新盈等農(nóng)場(chǎng)。其影響橡膠種植的不利因素主要是降水量較小、風(fēng)速較大、土壤養(yǎng)分含量不高，雖然其地勢(shì)較低，但最重要的是此類土壤有機(jī)質(zhì)和全氮等養(yǎng)分含量低，難以維持橡膠正常生長(zhǎng)的生理需要。對(duì)于此類區(qū)域可以選擇抗風(fēng)性較強(qiáng)的橡膠品種種植，并且加強(qiáng)農(nóng)田水利設(shè)施的建設(shè)，調(diào)節(jié)水分狀況，增施有機(jī)肥等，使之滿足橡膠生長(zhǎng)需要。

（4）不適宜區(qū)域。不適宜種植橡膠的面積為145 317.6 hm2，占農(nóng)用地面積的50.62%。該類農(nóng)用地面積較大，主要分布在沿海以及松濤水庫一帶。沿海地區(qū)由于降雨量較小、風(fēng)速較大等氣候條件以及養(yǎng)分含量低等條件限制而無法種植；在松濤水庫雖然土壤養(yǎng)分充足，但由于其地勢(shì)較高、坡度較大等條件限制橡膠的種植，且缺少管理，這些都使得松濤水庫附近的大片區(qū)域不適宜種植橡膠。

（5）由于八一農(nóng)場(chǎng)在儋州市橡膠生產(chǎn)中占有較大比重，且農(nóng)場(chǎng)在規(guī)劃、管理上較為精細(xì)，因此以八一農(nóng)場(chǎng)來檢驗(yàn)橡膠適宜性評(píng)價(jià)結(jié)果。根據(jù)資料數(shù)據(jù)[4]計(jì)算得到儋州市的橡膠單產(chǎn)為2.5 kg/年·株。結(jié)合八一農(nóng)場(chǎng)的實(shí)際情況，將八一農(nóng)場(chǎng)的橡膠產(chǎn)量進(jìn)行分等定級(jí)，將干膠大于3.0 kg/年·株的地塊確定為最適宜種植，適宜種植為2.5～3.0 kg/年·株，次宜種植為2.0～2.5 kg/年·株，小于2.0 kg/年·株的為不適宜種植。將確定好的橡膠產(chǎn)量水平和評(píng)價(jià)結(jié)果進(jìn)行對(duì)應(yīng)性比較，得到圖2。

（6）從圖2中不難看出，對(duì)應(yīng)“相同”的區(qū)域（17 855.8 hm2）明顯多于“不相同”區(qū)域（2 454.9 hm2），“相同”區(qū)域占驗(yàn)證區(qū)域的87.9%。對(duì)應(yīng)性結(jié)果表明：圖1的評(píng)價(jià)結(jié)果對(duì)于八一農(nóng)場(chǎng)橡膠的種植具有較高的適用性。

3 討論與結(jié)論

（1）種植適宜性評(píng)價(jià)涉及土地眾多性狀，即使所選的因子相同，采用不同的評(píng)價(jià)方法，其評(píng)價(jià)結(jié)果也有所不同。本研究從儋州市出發(fā)，以村一級(jí)為研究尺度，采用模糊數(shù)學(xué)的方法，根據(jù)橡膠樹正常生長(zhǎng)對(duì)環(huán)境條件的相關(guān)研究來確定各評(píng)價(jià)因子的隸屬函數(shù)的臨界值和閾值及散點(diǎn)型函數(shù)的隸屬度。并運(yùn)用主成分分析法，對(duì)橡膠種植適宜性進(jìn)行綜合評(píng)價(jià)，并制作了適宜性種植圖。通過適宜性評(píng)價(jià)，因地制宜的發(fā)展橡膠樹種植產(chǎn)業(yè)，從而提高橡膠的產(chǎn)出率，進(jìn)而創(chuàng)造出最大的經(jīng)濟(jì)、生態(tài)和社會(huì)效益。

此方法在實(shí)際應(yīng)用中仍有一些問題需要注意。首先，數(shù)據(jù)資料的獲取是本方法應(yīng)用的制約因素，一些重要的評(píng)價(jià)因子數(shù)據(jù)將在很大程度上影響評(píng)價(jià)結(jié)果的科學(xué)性；其次，評(píng)價(jià)因子的確定具有一定的局限性，因?yàn)樵u(píng)價(jià)因子及隸屬函數(shù)是根據(jù)數(shù)據(jù)資料以及專家意見加以平衡得出的，受影響較大。另外， 在評(píng)價(jià)過程中獲得的基礎(chǔ)數(shù)據(jù)越多、真實(shí)性越高， 所得出的評(píng)價(jià)單元將越零碎。出于評(píng)價(jià)目的的考慮， 為了更好地為農(nóng)業(yè)管理部門提供決策參考，對(duì)于種植適宜性評(píng)價(jià)結(jié)果， 在評(píng)價(jià)單元（地塊）上，還有待進(jìn)一步的研究。

（2）在今后的研究中，建立、完善一套科學(xué)、可行的適宜性評(píng)價(jià)方法仍有很多工作需要去做。同時(shí)，還應(yīng)加強(qiáng)人員的管理以及肥料的合理使用，以確保橡膠穩(wěn)產(chǎn)、高產(chǎn)。

參考文獻(xiàn)

[1] 許海平，傅國華. 我國天然橡膠產(chǎn)業(yè)發(fā)展趨勢(shì)[J]. 中國熱帶農(nóng)業(yè)，2007，（2）：15-16.

[2] 許道順，許升鋒. 經(jīng)濟(jì)一體化背景下海南省橡膠產(chǎn)業(yè)的發(fā)展戰(zhàn)略[J]. 海南金融，2006，（9）：25-28.

[3] 黃昌勇. 土壤學(xué)[M]. 北京：中國農(nóng)業(yè)出版社，2000.

[4] 海南省統(tǒng)計(jì)局. 海南統(tǒng)計(jì)年鑒2011[M]. 中國統(tǒng)計(jì)出版社，2011.

[5] 朱國宇，黃川友，華國春. 層次分析法在水環(huán)境規(guī)劃中的應(yīng)用[J]. 東北水利水電，2003，21（4）：1-7.

[6] 曹志洪，周建民. 中國土壤質(zhì)量[M]. 北京：科學(xué)出版社，2008.

[7] 劉紹貴，張?zhí)伊?，王興祥，等. 江西省余江縣土壤肥力變化和驅(qū)動(dòng)因素與對(duì)策研究[J]. 土壤通報(bào)，2006，（10）：869-874.

[8] 張 靜. 作物-地域多種組合中作物生態(tài)適宜性評(píng)價(jià)與權(quán)重配置方法的研究[D]. 南京：南京農(nóng)業(yè)大學(xué)，2005.

[9] 華南亞熱帶作物研究所. 中國橡膠栽培學(xué)[M]. 北京：中國科學(xué)出版社，1961.

[10] 何 康，黃宗道. 熱帶北緣橡膠所栽培[M]. 廣州：廣東科技出版社，1987.

[11] 王秉忠. 橡膠栽培學(xué)（第三版）[M]. 華南熱帶農(nóng)業(yè)大學(xué)出版社，2000.

[12] 韋娉婷. 耕地地力評(píng)價(jià)方法在橡膠園中的應(yīng)用研究——以海南農(nóng)墾八一片區(qū)為例[D]. ?？冢汉Ｄ洗髮W(xué)，2009.

[13] 歐 濱，羅 微，馬 利，等. 基于GIS的橡膠園地分等研究[J]. 熱帶農(nóng)業(yè)科學(xué)，2009，29（6）：9-14.

[14] 付 強(qiáng). 數(shù)據(jù)處理方法及其農(nóng)業(yè)應(yīng)用[M]. 北京：科學(xué)出版社，2006.

[15] 韓永學(xué). 特爾斐法與“拿來主義”[J]. 哈爾濱師專學(xué)報(bào)，2000（2）：67-68.endprint

（3）次適宜區(qū)域。次適宜種植橡膠的面積為40 934.9 hm2，占14.26%；該類農(nóng)用地集中分布在雅星、王五、東城、光村、峨蔓、木棠等鄉(xiāng)鎮(zhèn)以及紅嶺、西華、西聯(lián)、新盈等農(nóng)場(chǎng)。其影響橡膠種植的不利因素主要是降水量較小、風(fēng)速較大、土壤養(yǎng)分含量不高，雖然其地勢(shì)較低，但最重要的是此類土壤有機(jī)質(zhì)和全氮等養(yǎng)分含量低，難以維持橡膠正常生長(zhǎng)的生理需要。對(duì)于此類區(qū)域可以選擇抗風(fēng)性較強(qiáng)的橡膠品種種植，并且加強(qiáng)農(nóng)田水利設(shè)施的建設(shè)，調(diào)節(jié)水分狀況，增施有機(jī)肥等，使之滿足橡膠生長(zhǎng)需要。

（4）不適宜區(qū)域。不適宜種植橡膠的面積為145 317.6 hm2，占農(nóng)用地面積的50.62%。該類農(nóng)用地面積較大，主要分布在沿海以及松濤水庫一帶。沿海地區(qū)由于降雨量較小、風(fēng)速較大等氣候條件以及養(yǎng)分含量低等條件限制而無法種植；在松濤水庫雖然土壤養(yǎng)分充足，但由于其地勢(shì)較高、坡度較大等條件限制橡膠的種植，且缺少管理，這些都使得松濤水庫附近的大片區(qū)域不適宜種植橡膠。

（5）由于八一農(nóng)場(chǎng)在儋州市橡膠生產(chǎn)中占有較大比重，且農(nóng)場(chǎng)在規(guī)劃、管理上較為精細(xì)，因此以八一農(nóng)場(chǎng)來檢驗(yàn)橡膠適宜性評(píng)價(jià)結(jié)果。根據(jù)資料數(shù)據(jù)[4]計(jì)算得到儋州市的橡膠單產(chǎn)為2.5 kg/年·株。結(jié)合八一農(nóng)場(chǎng)的實(shí)際情況，將八一農(nóng)場(chǎng)的橡膠產(chǎn)量進(jìn)行分等定級(jí)，將干膠大于3.0 kg/年·株的地塊確定為最適宜種植，適宜種植為2.5～3.0 kg/年·株，次宜種植為2.0～2.5 kg/年·株，小于2.0 kg/年·株的為不適宜種植。將確定好的橡膠產(chǎn)量水平和評(píng)價(jià)結(jié)果進(jìn)行對(duì)應(yīng)性比較，得到圖2。

（6）從圖2中不難看出，對(duì)應(yīng)“相同”的區(qū)域（17 855.8 hm2）明顯多于“不相同”區(qū)域（2 454.9 hm2），“相同”區(qū)域占驗(yàn)證區(qū)域的87.9%。對(duì)應(yīng)性結(jié)果表明：圖1的評(píng)價(jià)結(jié)果對(duì)于八一農(nóng)場(chǎng)橡膠的種植具有較高的適用性。

3 討論與結(jié)論

（1）種植適宜性評(píng)價(jià)涉及土地眾多性狀，即使所選的因子相同，采用不同的評(píng)價(jià)方法，其評(píng)價(jià)結(jié)果也有所不同。本研究從儋州市出發(fā)，以村一級(jí)為研究尺度，采用模糊數(shù)學(xué)的方法，根據(jù)橡膠樹正常生長(zhǎng)對(duì)環(huán)境條件的相關(guān)研究來確定各評(píng)價(jià)因子的隸屬函數(shù)的臨界值和閾值及散點(diǎn)型函數(shù)的隸屬度。并運(yùn)用主成分分析法，對(duì)橡膠種植適宜性進(jìn)行綜合評(píng)價(jià)，并制作了適宜性種植圖。通過適宜性評(píng)價(jià)，因地制宜的發(fā)展橡膠樹種植產(chǎn)業(yè)，從而提高橡膠的產(chǎn)出率，進(jìn)而創(chuàng)造出最大的經(jīng)濟(jì)、生態(tài)和社會(huì)效益。

此方法在實(shí)際應(yīng)用中仍有一些問題需要注意。首先，數(shù)據(jù)資料的獲取是本方法應(yīng)用的制約因素，一些重要的評(píng)價(jià)因子數(shù)據(jù)將在很大程度上影響評(píng)價(jià)結(jié)果的科學(xué)性；其次，評(píng)價(jià)因子的確定具有一定的局限性，因?yàn)樵u(píng)價(jià)因子及隸屬函數(shù)是根據(jù)數(shù)據(jù)資料以及專家意見加以平衡得出的，受影響較大。另外， 在評(píng)價(jià)過程中獲得的基礎(chǔ)數(shù)據(jù)越多、真實(shí)性越高， 所得出的評(píng)價(jià)單元將越零碎。出于評(píng)價(jià)目的的考慮， 為了更好地為農(nóng)業(yè)管理部門提供決策參考，對(duì)于種植適宜性評(píng)價(jià)結(jié)果， 在評(píng)價(jià)單元（地塊）上，還有待進(jìn)一步的研究。

（2）在今后的研究中，建立、完善一套科學(xué)、可行的適宜性評(píng)價(jià)方法仍有很多工作需要去做。同時(shí)，還應(yīng)加強(qiáng)人員的管理以及肥料的合理使用，以確保橡膠穩(wěn)產(chǎn)、高產(chǎn)。

參考文獻(xiàn)

[1] 許海平，傅國華. 我國天然橡膠產(chǎn)業(yè)發(fā)展趨勢(shì)[J]. 中國熱帶農(nóng)業(yè)，2007，（2）：15-16.

[2] 許道順，許升鋒. 經(jīng)濟(jì)一體化背景下海南省橡膠產(chǎn)業(yè)的發(fā)展戰(zhàn)略[J]. 海南金融，2006，（9）：25-28.

[3] 黃昌勇. 土壤學(xué)[M]. 北京：中國農(nóng)業(yè)出版社，2000.

[4] 海南省統(tǒng)計(jì)局. 海南統(tǒng)計(jì)年鑒2011[M]. 中國統(tǒng)計(jì)出版社，2011.

[5] 朱國宇，黃川友，華國春. 層次分析法在水環(huán)境規(guī)劃中的應(yīng)用[J]. 東北水利水電，2003，21（4）：1-7.

[6] 曹志洪，周建民. 中國土壤質(zhì)量[M]. 北京：科學(xué)出版社，2008.

[7] 劉紹貴，張?zhí)伊郑跖d祥，等. 江西省余江縣土壤肥力變化和驅(qū)動(dòng)因素與對(duì)策研究[J]. 土壤通報(bào)，2006，（10）：869-874.

[8] 張 靜. 作物-地域多種組合中作物生態(tài)適宜性評(píng)價(jià)與權(quán)重配置方法的研究[D]. 南京：南京農(nóng)業(yè)大學(xué)，2005.

[9] 華南亞熱帶作物研究所. 中國橡膠栽培學(xué)[M]. 北京：中國科學(xué)出版社，1961.

[10] 何 康，黃宗道. 熱帶北緣橡膠所栽培[M]. 廣州：廣東科技出版社，1987.

[11] 王秉忠. 橡膠栽培學(xué)（第三版）[M]. 華南熱帶農(nóng)業(yè)大學(xué)出版社，2000.

[12] 韋娉婷. 耕地地力評(píng)價(jià)方法在橡膠園中的應(yīng)用研究——以海南農(nóng)墾八一片區(qū)為例[D]. 海口：海南大學(xué)，2009.

[13] 歐 濱，羅 微，馬 利，等. 基于GIS的橡膠園地分等研究[J]. 熱帶農(nóng)業(yè)科學(xué)，2009，29（6）：9-14.

[14] 付 強(qiáng). 數(shù)據(jù)處理方法及其農(nóng)業(yè)應(yīng)用[M]. 北京：科學(xué)出版社，2006.

[15] 韓永學(xué). 特爾斐法與“拿來主義”[J]. 哈爾濱師專學(xué)報(bào)，2000（2）：67-68.endprint