基于模糊聚類的沂水雪茄煙種植生態(tài)適應(yīng)性評價

孫飛虎 柳平增 孟憲勇 陳秀齋 張超 劉莉

摘要：通過生態(tài)適應(yīng)性評價篩選適合雪茄煙種植生產(chǎn)的區(qū)域，是提升雪茄原料產(chǎn)量與質(zhì)量、促進(jìn)雪茄產(chǎn)業(yè)發(fā)展的前提。以山東沂水雪茄莊園為例，運(yùn)用生育期平均氣溫、晝夜氣溫差、生育期≥10℃活動積溫、生育期相對濕度、生育期降雨量以及生育期日均日照時數(shù)等雪茄種植生產(chǎn)的重要?dú)夂蛑笜?biāo)，對哈瓦那、比那爾得里奧、什邡等11個國內(nèi)外優(yōu)質(zhì)雪茄煙區(qū)進(jìn)行模糊聚類并通過隸屬度函數(shù)求得各地區(qū)相應(yīng)的CFI值。模糊聚類結(jié)果表明，可將11個區(qū)域劃分為4類，其中沂水氣候條件與哈瓦那、比那爾得里奧、什邡、多米尼加、巴伊亞州等國內(nèi)外優(yōu)質(zhì)雪茄種植地區(qū)相似，它們的氣候適宜性指數(shù)CFI值依次為沂水（0.931）、哈瓦那（0.963）、比那爾德里奧（0.956）、什邡（0.944）、多米尼加（0.932）、巴伊亞州（0.892），尤其沂水與多米尼加的CFI指數(shù)最為相近。由此得出，沂水雪茄莊園具有種植優(yōu)質(zhì)雪茄煙的生態(tài)氣候條件。

關(guān)鍵詞：模糊聚類；雪茄；氣候適宜性；生態(tài)適應(yīng)性；隸屬度函數(shù)

中圖分類號：S572

文獻(xiàn)標(biāo)識碼：A

文章編號：2095-5553 （2024） 06-0330-07

收稿日期：2022年11月21日

修回日期：2023年2月17日

*基金項(xiàng)目：山東省農(nóng)業(yè)重大應(yīng)用技術(shù)創(chuàng)新項(xiàng)目（SD2019ZZ019）；2019年度山東省重點(diǎn)研發(fā)計(jì)劃（公益類專項(xiàng)）項(xiàng)目（2019GNC106103）；山東省科技特派員項(xiàng)目（2020KJTPY078）；山東省重大科技創(chuàng)新工程項(xiàng)目（2019JZZY010713）

第一作者：孫飛虎，男，1999年生，安徽濉溪人，碩士研究生；研究方向?yàn)榇髷?shù)據(jù)與人工智能。E-mail： sunflytiger@163.com

通訊作者：柳平增，男，1968年生，山東濟(jì)南人，教授，博導(dǎo)；研究方向?yàn)檗r(nóng)業(yè)大數(shù)據(jù)、農(nóng)業(yè)物聯(lián)網(wǎng)與智慧農(nóng)業(yè)。E-mail： pzliu@sdau.edu.cn

Ecological suitability evaluation of cigar planting in Yishui based on fuzzy clustering

Sun Feihu^{1， 2， 3}， Liu Pingzeng^{1， 2， 3}， Meng Xianyong^{1， 2， 3}， Chen Xiuzhai⁴， Zhang Chao⁴， Liu Li⁴

（1. School of Information Science and Engineering， Shandong Agricultural University， Taian， 271018， China; 2. Key Laboratory of Intelligent Agriculture Technology of the Yellow and Huaihai Sea， Ministry of Agriculture and Rural Affairs， Taian， 271018， China; 3. Agricultural Big Data Research Center， Shandong Agricultural University， Taian， 271018， China; 4. Shandong Linyi Tobacco Co.， Ltd.， Linyi， 276003， China）

Abstract： Screening suitable areas for cigar planting and production through ecological adaptability evaluation is a prerequisite for improving the yield and quality of cigar raw materials and promoting the development of cigar industry. Taking Yishui Cigar Estate in Shandong Province as an example， using important climate indicators such as average temperature during the growth period， temperature difference between day and night， active accumulated temperature of ≥10℃ during the growth period， relative humidity during the growth period， rainfall during the growth period and average daily sunshine hours during the growth period， fuzzy clustering was carried out on 11 high-quality cigar tobacco areas at home and abroad， including Havana， Binard de Río and Shifang， and the corresponding CFI values of each region were obtained through the membership function. The fuzzy clustering results showed that 11 regions could be divided into 4 categories， among which the climatic conditions of Yishui were similar to those of Havana， Pinar de Río， Dominica， Bahia， Shifang and other high-quality cigar growing areas at home and abroad， and their climate suitability index CFI values were 0.931， 0.963， 0.956， 0.944， 0.932 and 0.892， respectively， especially Yishui and Dominicas CFI index was the most similar. It is concluded that Yishui Cigar Estate has the ecological and climatic conditions for growing high-quality cigars.

Keywords： fuzzy clustering; cigar; climatic suitability; ecological suitability; membership function

0 引言

近年來，我國雪茄產(chǎn)業(yè)發(fā)展迅速^[1]，已逐步成為我國煙草行業(yè)新的經(jīng)濟(jì)增長點(diǎn)^[2]。目前大部分國產(chǎn)雪茄煙原料依賴于國外進(jìn)口，雪茄煙葉特別是高質(zhì)量茄衣短缺一定程度上制約了我國雪茄產(chǎn)業(yè)的健康、可持續(xù)發(fā)展，已成為制約雪茄產(chǎn)業(yè)發(fā)展的最大瓶頸^[3]。雪茄煙作為晾煙的一種，對生態(tài)和氣候條件要求更高^[4]。生態(tài)因素對雪茄煙的產(chǎn)量與品質(zhì)皆有影響，尤其是對煙葉品質(zhì)影響最大^{[5， 6]}，雪茄煙生長的生態(tài)條件主要包括氣候條件與土壤條件^[7]，而氣候條件對雪茄煙的影響尤為顯著^[8]，因此篩選出適宜種植優(yōu)質(zhì)雪茄煙的產(chǎn)區(qū)對解決雪茄原料短缺、推動雪茄產(chǎn)業(yè)健康可持續(xù)發(fā)展具有重要意義。本研究擬在分析沂水氣候特點(diǎn)基礎(chǔ)上，研究煙草種植生態(tài)適宜性評價方法。

1 材料與方法

1.1 研究區(qū)概況

山東省沂水縣地處魯中南地區(qū)，沂蒙山區(qū)北部，沂山南麓，屬魯中南低山陵區(qū)，地勢西北高、東南低，西北部為低山區(qū)，東部和東北部為丘陵，中南部為平原。地理坐標(biāo)東經(jīng)118°11′～119°03′，北緯35°36′～36°12′。沂水縣屬暖溫帶季風(fēng)氣候區(qū)，年均氣溫為12.5℃～13.5℃，晝夜溫差大，無霜期長達(dá)240d以上，年均降雨量為785.1mm，7—9月降雨量充沛，約占全年降雨量的70%，年平均相對濕度為65%，年均日照率為 55%～59%，年日照總時長為2420.1h。雪茄莊園地處沂水縣西北部諸葛鎮(zhèn)跋山水庫北岸，境內(nèi)為低山區(qū)，海拔在159～524m之間，水資源充沛，黃煙為該地區(qū)的主要經(jīng)濟(jì)作物之一。

1.2 研究方法

模糊聚類分析是根據(jù)客觀事物間的特征、親疏程度、相似性，通過建立模糊相似關(guān)系對客觀事物進(jìn)行聚類分析，其在解決區(qū)域劃分、氣候分類等問題上都有一定的應(yīng)用^[9]。本文在文獻(xiàn)[10-13]等的基礎(chǔ)上，提出模糊聚類與隸屬函數(shù)相結(jié)合的氣候適宜性評價方法。首先利用模糊聚類對采集的數(shù)據(jù)樣本進(jìn)行聚類，然后采用隸屬度函數(shù)對和沂水聚到一類的雪茄煙區(qū)進(jìn)行評分。

1.2.1 模糊聚類分析

1） 數(shù)據(jù)標(biāo)準(zhǔn)化處理。由于各性狀指標(biāo)的量綱和數(shù)量級不同，不能直接計(jì)算分析。為消除單位的影響，使數(shù)據(jù)具有可比性，需采用標(biāo)準(zhǔn)差法對原始數(shù)據(jù)進(jìn)行處理。

F_ij=F′_ij-F′_jS_j（1）

式中：F_ij——初始數(shù)據(jù)；

F′_j——某一性狀的均值；

S_j——該性狀的標(biāo)準(zhǔn)差。

2） 建立模糊相似矩陣。設(shè)U為被分類對象的全體，r_ij為在U上建立的相似系數(shù)。采用歐式距離法由式（2）計(jì)算r_ij，并建立模糊相似矩陣R。

r_ij=1i=j

1-∑mk=1（F_ik-F_jk）²Di≠j（2）

其中D=max∑mk=1（F_ik-F_jk）²。

3） 求模糊等價矩陣R^*。模糊等價矩陣：即R自乘RR=R²，再自乘R²R²=R⁴，直到R^2k=R^k=R^*為止，則R^*便是一個模糊等價矩陣。

4） 模糊聚類分析。用R^*_λ=（λ_ij），取適宜實(shí)數(shù)λ∈[0，1]。當(dāng)r_ij≥λ時，λ_ij=1；當(dāng)r_ij≤λ時，λ_ij=0。對得到的結(jié)果矩陣進(jìn)行排序劃分，即得到在λ水平上的分類。

1.2.2 基于模糊數(shù)學(xué)求氣候適宜性指數(shù)

1） 評價指標(biāo)隸屬函數(shù)的建立。通過模糊統(tǒng)計(jì)法確定各評價指標(biāo)的隸屬函數(shù)類型、最優(yōu)值和最差值，然后計(jì)算出各評價因子的隸屬度。涉及的隸屬度函數(shù)為拋物線型（式（3））與升梯型（式（4）），其中拋物線型函數(shù)適用于雪茄煙生長既有上限又有下限的因子，升梯型函數(shù)適用于雪茄煙有下限無上限的因子。本研究涉及隸屬函數(shù)如式（3）所示，其中a為下限，b為最優(yōu)值下限，c為最優(yōu)值上限，d為上限，a＜b＜c＜d。

f（x）=1.0b≤x

x-ab-a×0.9+0.1a

d-xd-c×0.9+0.1c

0.1x≥d，x≤a（3）

f（x）=1.0x

x-ab-a×0.9+0.1a≤x

0.1x

本研究首先通過數(shù)據(jù)統(tǒng)計(jì)分析以及咨詢專家篩選出進(jìn)行適宜性評價的氣候因子，然后通過層次分析法構(gòu)建層次結(jié)構(gòu)模型，將與評價有關(guān)的要素分為目標(biāo)層（雪茄生態(tài)適應(yīng)性）、準(zhǔn)則層（氣候條件）和指標(biāo)層（生育期平均氣溫，生育期≥10℃活動積溫，生育期相對濕度，生育期平均降水量，生育期平均日照時數(shù)），結(jié)合專家賦值法來進(jìn)行判斷矩陣的構(gòu)造，最后進(jìn)行一致性檢驗(yàn)并計(jì)算各指標(biāo)權(quán)重值。所涉及的各氣候因子的權(quán)重值、隸屬函數(shù)類型及拐點(diǎn)如表1所示。

表1 各評價指標(biāo)的權(quán)重、隸屬函數(shù)類型及拐點(diǎn)

Tab. 1 Weight， membership function type and

inflection point of each evaluation index氣候因子權(quán)重隸屬函數(shù)

類型最優(yōu)值最差值生育期平均氣溫0.41拋物線型24℃～

30℃20℃，

35℃生育期≥10℃活動積溫0.29升梯型2 800℃2 600℃生育期相對濕度0.06拋物線型70％～

80％50％，

90％生育期平均降水量0.08拋物線型500～

600 mm300 mm，

800 mm生育期平均日照時數(shù)0.16拋物線型600～

800 h300 h，

1 000 h

2） 氣候適宜性指數(shù)計(jì)算。利用隸屬函數(shù)氣候因子指標(biāo)值與權(quán)重系數(shù)加成法則，分別計(jì)算通過聚類得到的與沂水同一分類的地區(qū)的氣候適宜性指數(shù)。氣候適宜性指數(shù)CFI計(jì)算如式（5）所示。

CFI=∑ni=1P_jX_ijj=1，2，3，…，m（5）

式中：n——評價指標(biāo)的總數(shù)；

P_j——第j個評價因子的權(quán)重值；

X_ij——第j個評價因子的隸屬度。

本研究將優(yōu)質(zhì)雪茄種植氣候適宜性按照CFI值分為4個等級，如表2所示。

CFI的取值在0.10～1.00之間，CFI值越高，說明該地區(qū)氣候條件越適宜種植雪茄煙；CFI值越低，說明該地區(qū)氣候條件越不適宜種植雪茄煙。通過計(jì)算得出的CFI值對沂水雪茄煙區(qū)生態(tài)適宜性進(jìn)行綜合評價。

1.3 數(shù)據(jù)來源

優(yōu)質(zhì)雪茄煙葉生產(chǎn)需要滿足特定的日照、降水以及溫濕度等條件^[14]，研究數(shù)據(jù)來源于weather-atlas（https：//www.weather-atlas.com/zh）。

2 結(jié)果與分析

將沂水與哈瓦那、比那爾得里奧、多米尼加、巴伊亞州、什邡、康涅狄格、津巴布韋、印度尼西亞、菲律賓、儋州根據(jù)生育期平均氣溫，生育期≥10℃活動積溫，生育期相對濕度，生育期平均降水量，生育期平均日照時數(shù)等指標(biāo)進(jìn)行模糊聚類，把與沂水聚為一類的地區(qū)通過隸屬度函數(shù)對各指標(biāo)進(jìn)行評價。

2.1 雪茄煙區(qū)模糊聚類分析

1） 數(shù)據(jù)標(biāo)準(zhǔn)化。根據(jù)式（1）進(jìn)行數(shù)據(jù)標(biāo)準(zhǔn)化處理的結(jié)果如表3所示。

2） 模糊相似矩陣和模糊等價矩陣。模糊相似矩陣見表4。

將標(biāo)準(zhǔn)化處理過的數(shù)據(jù)根據(jù)式（2）進(jìn)行模糊相似矩陣和模糊等價矩陣的求解。求得的矩陣為n階對稱方陣，其中n代表所進(jìn)行聚類的樣本總數(shù)。

當(dāng)k=8時，R⁸=R¹⁶，R⁸為模糊等價矩陣，模糊等價矩陣見表5。

3） 聚類分析。用R^*_λ=（λ_ij），取適宜實(shí)數(shù)λ∈[0，1]。當(dāng)r_ij≥λ時，λ_ij=1；當(dāng)r_ij≤λ時，λ_ij=0。對得到的結(jié)果矩陣進(jìn)行排序劃分，即得到在λ水平上的分類。

在模糊等價矩陣R^*中，取λ=0.68，當(dāng)r^*≥0.68時，r^*=1；當(dāng)r^*≤0.68時，r^*=0。即得到λ=0.68水平下的截距矩陣t（R），如表6所示。

從表6可得，當(dāng)λ=0.68時，11個區(qū)域可被劃分為4類，結(jié)果如表7所示。模糊聚類結(jié)果表明，沂水氣候條件與哈瓦那，比那爾得里奧，多米尼加，巴伊亞州，什邡等國內(nèi)外優(yōu)質(zhì)雪茄種植地區(qū)比較相似，具備生產(chǎn)優(yōu)質(zhì)雪茄的氣候條件。

2.2 氣候適宜性評價

將聚類分析得到的和沂水處在同一類的地區(qū)，即為Ⅰ類的地區(qū)進(jìn)行氣候適宜性評價。

2.2.1 氣候條件分析

雪茄煙葉是喜溫作物，在6℃～35℃均能生長，但優(yōu)質(zhì)的雪茄煙對溫度的要求非常嚴(yán)格，一般以20℃～30℃為宜；積溫反映了雪茄煙生長所需的熱量條件，積溫不足會導(dǎo)致煙葉內(nèi)物質(zhì)的形成和轉(zhuǎn)化效率降低，煙葉生育期氣溫≥10℃的活動積溫≥2800℃為煙葉適宜生長的熱量條件，最低積溫不低于2600℃。

圖1為沂水5—8月份平均氣溫變化，可以看出5—8月份沂水平均氣溫均處在20℃～30℃之間，總積溫達(dá)到了2967.5℃。說明沂水具有雪茄煙葉生長的適宜氣溫條件，且能為煙葉的生長發(fā)育提供充足的積溫。

煙葉生長大田期的水分供應(yīng)，對煙葉生長發(fā)育及煙葉質(zhì)量、產(chǎn)量影響很大。降水對烤煙的影響，主要決定于雨量的分布：移栽至旺長前，煙株的耗水量少，這段時間的月降雨量在70～100mm，土壤相對濕度為50%～60%較為理想；進(jìn)入旺長期后耗水增大，月降雨量在100～200mm，土壤相對濕度為70%～80%對烤煙生長、干物質(zhì)積累最為有利。

圖2為沂水5—8月份平均降雨量及平均相對濕度變化?？梢钥闯?，沂水5月份平均降雨量為76mm，相對濕度為57%；6—8月份，平均降雨量為158mm，平均相對濕度為77%。兩段時期內(nèi)的降雨量以及濕度均處在理想范圍內(nèi)，沂水地區(qū)滿足雪茄煙生長所需的降雨量以及濕度條件。

雪茄煙對日照的要求較高，優(yōu)質(zhì)煙葉的生產(chǎn)需要充足且不強(qiáng)烈的光照。沂水平均日照率約為57%，日均光照時長約為7.7h，滿足煙葉日照率40%以上、每天光照7～10h的需求。成熟期由于沂水降雨量較多，使其擁有充足而不強(qiáng)烈的光照，利于優(yōu)質(zhì)雪茄煙的生產(chǎn)。

圖3為沂水5—8月份的平均日照時長變化，可以看出沂水具有雪茄生長的良好光照條件。

2.2.2 氣候適宜性綜合評價

表8為氣候適宜性評價所需數(shù)據(jù)，結(jié)合表2氣候適宜性參評指標(biāo)計(jì)算得出各氣候指標(biāo)得分如表9所示。

由表9可知，CFI由高到低地區(qū)分別是哈瓦那、比那爾得里奧、什邡、多米尼加、沂水以及巴伊亞州，所選6個參照煙區(qū)按照適宜性等級劃分結(jié)果為：最適宜[哈瓦那，比那爾得里奧]，適宜[什邡，多米尼加，沂水，巴伊亞州]。從CFI數(shù)值層面看，沂水氣候與哈瓦那、比那爾得里奧、什邡等老牌雪茄煙生產(chǎn)區(qū)存在一定差距，和多米尼加的氣候得分非常相近，且強(qiáng)于巴伊亞州。可以考慮從以上地區(qū)進(jìn)行雪茄煙品種的選取以及學(xué)習(xí)雪茄煙的種植經(jīng)驗(yàn)。

3 討論

氣候是優(yōu)質(zhì)雪茄煙生長的基礎(chǔ)條件，不僅對雪茄的產(chǎn)量有很大影響，還和雪茄煙的品質(zhì)優(yōu)劣密切相關(guān)^{[15， 16]}。沂水與哈瓦那、比那爾得里奧、多米尼加、巴伊亞州、什邡具有較相似的氣候條件，其中沂水與多米尼加、巴伊亞州、什邡的適宜性等級為適宜。

3.1 氣溫與積溫

在氣溫與積溫方面，生育期的沂水平均溫度處在24℃～30℃之間，≥10℃活動積溫為2967.5℃，均滿足最優(yōu)的雪茄煙生長條件，這兩項(xiàng)指標(biāo)均為最優(yōu)，達(dá)到了最適宜的氣溫與積溫條件，能為雪茄煙的生長提供充足且溫和的熱能。雪茄煙生長期伴有極端高溫天氣發(fā)生的風(fēng)險，圖4為2011—2022年5—6月每月的極端氣溫（≥35℃）波動，可以看出沂水縣在5—8月均曾有極端氣溫出現(xiàn)，并且近些年所出現(xiàn)的極端氣溫溫度有上升的趨勢。其中6月的極端氣溫出現(xiàn)的頻率較高，極端氣溫溫度已連續(xù)四年處于一個較高水平，8月的極端氣溫也在2022年達(dá)到了歷史峰值。

3.2 降雨量與濕度

沂水的降雨量與濕度均符合雪茄煙生長所需的最優(yōu)條件。沂水屬魯中南低山陵區(qū)且屬近海地域，夏天容易發(fā)生強(qiáng)對流天氣和受到臺風(fēng)影響，出現(xiàn)暴雨以及大風(fēng)天氣。尤其地勢低洼的地區(qū)易受強(qiáng)降雨侵襲導(dǎo)致積水，從而影響煙葉中的還原糖、總糖、蛋白質(zhì)、氮以及煙堿等的含量，進(jìn)而導(dǎo)致煙葉的品質(zhì)與產(chǎn)量下降^[17]。

3.3 日照

沂水的日照時長得分相對其余雪茄煙產(chǎn)區(qū)較高，且超出了雪茄煙生長最優(yōu)值的上限。不過沂水縣在生長成熟期降雨較多，此時光照充足但不強(qiáng)烈，有利于優(yōu)質(zhì)雪茄煙的生長。

3.4 建議

沂水縣基本具備優(yōu)質(zhì)雪茄種植的氣候條件，但是由于沂水地處魯中南山區(qū)，夏季多強(qiáng)對流天氣，暴雨和冰雹時有發(fā)生，同時極端氣溫的出現(xiàn)會對煙葉的生長發(fā)育產(chǎn)生不利影響，因此在雪茄的品種選擇方面選擇能應(yīng)對氣象災(zāi)害的雪茄煙種、做好氣象災(zāi)害預(yù)測預(yù)警與防范措施對雪茄煙的種植具有重要意義；同時本研究只針對了沂水的氣候因素進(jìn)行雪茄種植適宜性的評價，最好還應(yīng)考慮地形、土壤以及社會經(jīng)濟(jì)等因素的影響，這樣可以更加準(zhǔn)確的對沂水雪茄煙的種植作適宜性評價，充分結(jié)合沂水氣候的特點(diǎn)，為做好合理種植以及提高經(jīng)濟(jì)效益提供合理意見。

4 結(jié)論

本研究創(chuàng)新性地提出一種基于模糊聚類的雪茄煙種植生產(chǎn)生態(tài)適宜性評價方法。首先通過聚類將氣候類別相近的雪茄產(chǎn)區(qū)劃分到一起，聚類結(jié)果將哈瓦那、比那爾得里奧、多米尼加、巴伊亞州、什邡以及沂水等產(chǎn)區(qū)劃分為一類。然后通過適應(yīng)性評價方法對同類的主產(chǎn)區(qū)進(jìn)行適宜性評價，評價結(jié)果氣候適宜性指數(shù)CFI值依次為沂水（0.931）、哈瓦那（0.963）、比那爾德里奧（0.956）、什邡（0.944）、多米尼加（0.932）、巴伊亞州（0.892）。

1） 沂水縣雪茄煙種植的光熱條件優(yōu)越，降雨量充沛，具有和國內(nèi)外等優(yōu)質(zhì)煙區(qū)較為接近的氣候條件。

2） 在進(jìn)行雪茄選種方面可嘗試與沂水CFI值較相近的地區(qū)的雪茄煙種以及學(xué)習(xí)其種植技術(shù)。

3） 本研究所提出的先聚類再評價的方法，為雪茄煙區(qū)的選取提供了一種新的思路與方法，該方法的推廣可提升雪茄煙產(chǎn)區(qū)選擇效率，并對推動雪茄產(chǎn)業(yè)的健康發(fā)展具有重大意義。

參 考 文 獻(xiàn)

[1]陳棟， 李猛， 王榮浩， 等. 國產(chǎn)雪茄茄芯煙葉研究進(jìn)展[J]. 揚(yáng)州大學(xué)學(xué)報(bào)（農(nóng)業(yè)與生命科學(xué)版）， 2019， 40（1）： 83-90.

Chen Dong， Li Meng， Wang Ronghao， et al. Progress of the domestic cigar filler tobacco [J]. Journal of Yangzhou University（Agricultural and Life Science Edition）， 2019， 40（1）： 83-90.

[2]閆新甫， 王欣， 孔勁松， 等. 全國近十年晾煙生產(chǎn)及市場變化分析[J]. 中國煙草科學(xué)， 2021， 42（2）： 98-104.

Yan Xinfu， Wang Xin， Kong Jinsong， et al. Changing trends of production and market of air-cured tobaccos of China in last 10 years [J]. Chinese Tobacco Science， 2021， 42（2）： 98-104.

[3]蔡斌， 耿召良， 高華軍， 等. 國產(chǎn)雪茄原料生產(chǎn)技術(shù)研究現(xiàn)狀[J]. 中國煙草學(xué)報(bào)， 2019， 25（6）： 110-119.

Cai Bin， Geng Zhaoliang， Gao Huajun， et al. Research progress of production technologies of cigar tobaccos in China [J]. Acta Tabacaria Sinica， 2019， 25（6）： 110-119.

[4]王琰琰， 劉國祥， 向小華， 等. 國內(nèi)外雪茄煙主產(chǎn)區(qū)及品種資源概況[J]. 中國煙草科學(xué)， 2020， 41（3）： 93-98.

Wang Yanyan， Liu Guoxiang， Xiang Xiaohua， et al. Overview of main cigar production areas and variety resources at domestic and overseas [J]. Chinese Tobacco Science， 2020， 41（3）： 93-98.

[5]秦艷青， 李愛軍， 范靜苑， 等. 優(yōu)質(zhì)雪茄茄衣生產(chǎn)技術(shù)探討[J]. 江西農(nóng)業(yè)學(xué)報(bào)， 2012， 24（7）： 101-103.

Qin Yanqing， Li Aijun， Fan Jingyuan， et al. Discussion on production technology of high-quality cigar wrapper [J]. Acta Agriculturae Jiangxi， 2012， 24（7）： 101-103.

[6]邵麗， 晉艷， 楊宇虹， 等. 生態(tài)條件對不同烤煙品種煙葉產(chǎn)質(zhì)量的影響[J]. 煙草科技， 2002（10）： 40-45.

Shao Li， Jin Yan， Yang Yuhong， et al. Influences of ecological conditions on the yield and quality of different flue-cured cultivars [J]. Tobacco Science & Technology， 2002（10）： 40-45.

[7]倪明， 濮永瑜， 何翔， 等. 保山市植煙土壤養(yǎng)分時空變異特征及土壤肥力適宜性評價[J]. 西南農(nóng)業(yè)學(xué)報(bào)， 2020， 33（10）： 2303-2309.

Ni Ming， Pu Yongyu， He Xiang， et al. Soil nutrients spatial-temporal variability and soil fertility suitability in baoshan tobacco-planting area [J]. Southwest China Journal of Agricultural Sciences， 2020， 33（10）： 2303-2309.

[8]薛劍波， 符云鵬， 尹永強(qiáng). 影響煙草中煙堿含量的因素及調(diào)控措施[J]. 安徽農(nóng)業(yè)科學(xué)， 2005（6）： 1053-1055.

Xue Jianbo， Fu Yunpeng， Yin Yongqiang. Factor and controlling measures influencing the nicotine content in tobacco leaf [J]. Journal of Anhui Agricultural Sciences， 2005（6）： 1053-1055.

[9]Moharana P C， Jena R K， Pradhan U K， et al. Geostatistical and fuzzy clustering approach for delineation of site-specific management zones and yield-limiting factors in irrigated hot arid environment of India [J]. Precision Agriculture， 2020， 21（2）： 426-448.

[10]張誼寒， 張晨東， 焦芳嬋， 等. 雪茄外包皮煙在云南種植的適宜氣候區(qū)初步篩選研究[J]. 西南農(nóng)業(yè)學(xué)報(bào)， 2012， 25（6）： 2005-2009.

Zhang Yihan， Zhang Chendong， Jiao Fangchan， et al. Preliminary study on suitable planting regions of cigar wrapper tobacco in Yunnan Province [J]. Southwest China Journal of Agricultural Sciences， 2012， 25（6）： 2005-2009.

[11]馮曉利， 何偉， 蔣貴國， 等. 基于模糊綜合評價法的雙流縣農(nóng)用地適宜性評價[J]. 西南農(nóng)業(yè)學(xué)報(bào)， 2012， 25（3）： 982-988.

Feng Xiaoli， He Wei， Jiang Guiguo， et al. Fuzzy comprehensive assessment analysis of agricultural land suitability evaluation in Shuangliu County [J]. Southwest China Journal of Agricultural Sciences， 2012， 25（3）： 982-988.

[12]李滟芳， 周桂園， 李輝， 等. 基于GIS技術(shù)和灰關(guān)聯(lián)度法的文山植煙生態(tài)適宜性評價[J]. 西南農(nóng)業(yè)學(xué)報(bào)， 2012， 25（2）： 408-411.

Li Yanfang， Zhou Guiyuan， Li Hui， et al. Ecological suitability evaluation of flue-cured tobacco planted in Wenshan prefecture based on GIS techniques and grey correlative analysis [J]. Southwest China Journal of Agricultural Sciences， 2012， 25（2）： 408-411.

[13]胡穎梅， 徐磊， 袁帥， 等. 陜西烤煙氣候適宜性時空變化特征及風(fēng)險分析[J]. 煙草科技， 2021， 54（7）： 41-50.

Hu Yingmei， Xu Lei， Yuan Shuai， et al. Spatiotemporal variation characteristics of climatic suitability for flue-cured tobacco in Shaanxi and corresponding risk analysis [J]. Tobacco Science & Technology， 2021， 54（7）： 41-50.

[14]許彥平， 姚曉紅， 王潤元， 等. 天水烤煙適生適產(chǎn)的氣候區(qū)劃[J]. 中國農(nóng)業(yè)氣象， 2008（2）： 188-190.

Xu Yanping， Yao Xiaohong， Wang Runyuan， et al. The climatic demarcation for flue-cured tobacco suitably lives in Tianshui [J]. Chinese Journal of Agrometeorology， 2008（2）： 188-190.

[15]張久權(quán)， 梁洪波， 董建新， 等. 基于GIS和模糊集理論的四川烤煙生態(tài)適應(yīng)性評價[J]. 中國煙草科學(xué)， 2016， 37（3）： 8-14.

[16]陳偉， 王三根， 唐遠(yuǎn)駒， 等. 不同煙區(qū)烤煙化學(xué)成分的主導(dǎo)氣候影響因子分析[J]. 植物營養(yǎng)與肥料學(xué)報(bào)， 2008（1）： 144-150.

Chen Wei， Wang Sangen， Tang Yuanju， et al. Analysis of the dominant climatic factors influencing the chemical compositions of flue-cured tobacco in different tobacco-growing areas [J]. Journal of Plant Nutrition and Fertilizers， 2008（1）： 144-150.

[17]孫延國， 馬興華， 黃擇祥， 等. 煙草溫光特性研究與利用： Ⅰ.氣象因素對山東煙區(qū)主栽品種生育期的影響[J]. 中國煙草科學(xué)， 2020， 41（1）： 30-37.

Sun Yanguo， Ma Xinghua， Huang Zexiang， et al. Investigation and utilization of temperature and light characteristics of tobacco： Ⅰ. Effects of meteorological factors on growth period of main varieties in Shandong tobacco-growing areas [J]. Chinese Tobacco Science， 2020， 41（1）： 30-37.