劉中新 宋云

摘要：根據(jù)收集到的獼猴桃果實(shí)生長后期采集的獼猴桃可溶性固形物含量檢測數(shù)據(jù)資料，使用前期氣象資料設(shè)計(jì)并計(jì)算獼猴桃成熟熵。選取獼猴桃當(dāng)年末花期次日至檢測日前一日≥0℃積溫和累計(jì)日照時(shí)數(shù)作為計(jì)算獼猴桃成熟熵的氣象因子，按其與獼猴桃可固含量的回歸系數(shù)權(quán)重加權(quán)求和，得到成熟熵指數(shù)（H）的計(jì)算公式為：H=0.28ΣT+0.72ΣS，對獼猴桃成熟熵與可溶性固形物檢測含量進(jìn)行回歸分析，建立可溶性固形物檢測含量與獼猴桃成熟熵回歸方程：G=-3.9875+0.0060H。根據(jù)回歸方程計(jì)算得到的可溶性固形物含量值（G）與檢測值最小誤差0.03%，最大誤差-0.96%，平均誤差0.083%。參照獼猴桃成熟要求可溶性固形物含量≥8.0%標(biāo)準(zhǔn)，計(jì)算得到獼猴桃所需成熟熵為1949.06，預(yù)報(bào)2017—2019年獼猴桃成熟期與實(shí)況誤差-1～2天，平均誤差0.7天（偏遲）。運(yùn)用本方法可根據(jù)氣象因子計(jì)算當(dāng)期獼猴桃可溶性固形物含量，結(jié)合后期氣象要素預(yù)報(bào)，可提前計(jì)算預(yù)報(bào)獼猴桃成熟期，籍此可以擺脫對折光分析檢測儀器的依賴，又能提前預(yù)報(bào)獼猴桃成熟期和上市期，為避免早采而造成對獼猴桃品牌形象的損害提供科學(xué)依據(jù)。

關(guān)鍵詞：獼猴桃；成熟期；預(yù)報(bào)；可溶性固形物；成熟熵；積溫；日照

中圖分類號：S165+.26文獻(xiàn)標(biāo)志碼：A論文編號：cjas2020-0237

Forecasting Method of Kiwifruit Ripening Period Under Facility Irrigation

Liu Zhongxin1， Song Yun2

（1Weather Bureau of Xishui County， Xishui 438200， Hubei， China；

2Meteorological Bureau of Huanggang， Hubei Province， Huanggang 438000， Hubei， China）

Abstract： Based on the data of soluble solid content of kiwifruit collected in the later period of kiwifruit growth， the maturity entropy of kiwifruit was designed and calculated by using the meteorological data in the earlier period. The accumulated temperature≥0℃and the accumulated sunshine time from the second day of the end flowering period to the day before the test date were selected as the meteorological factors to calculate the entropy of kiwifruit ripening， the formula for calculating the maturity entropy index （H） was obtained as： H= 0.28ΣT+0.72ΣS， to establish the regression equation between the content of soluble solid and the entropy of kiwifruit ripening： G=-3.9875+0.0060H. According to the regression equation， the minimum error was 0.03%， the maximum error was -0.96%， and the average error was 0.083%. According to the requirement of kiwifruit ripeness for soluble solid content≥8.0%， the maturity entropy of kiwifruit was calculated to be 1949.06， and the error of kiwifruit ripeness prediction was - 1～2 days in 2017-2019， with an average error of 0.7 days（delayed） . This method can be used to calculate the soluble solid content of kiwifruit in the current period according to the meteorological factors. Combined with the forecast of the later meteorological elements， the ripening period of kiwifruit could be calculated and predicted in advance， it could also forecast the ripening period and the marketing period of kiwifruit in advance， and avoid the damage to the brand image of kiwifruit caused by early harvest.

Keywords： Kiwifruit； Maturity Period； Forecast； Soluble Solids； Maturity Entropy； Accumulated Temperature； Sunshine

0引言

獼猴桃（Actinidia chinensis Planch），也稱狐貍桃、羊桃、毛木果、奇異果等，因獼猴喜食，故名獼猴桃，是一種品質(zhì)鮮嫩，營養(yǎng)豐富，風(fēng)味鮮美的水果。獼猴桃是20世紀(jì)人工馴化栽培野生果樹最有成就的四大果種之一[1-3]，其漿果具有特殊的風(fēng)味，含多種氨基酸、多種礦物質(zhì)、維生素C，藤、根、葉及果均可入藥，具有很高的經(jīng)濟(jì)價(jià)值和使用價(jià)值。

獼猴桃現(xiàn)代種植技術(shù)已不斷發(fā)展，人工規(guī)?；N植面積不斷擴(kuò)大。有關(guān)獼猴桃種植技術(shù)及與氣象條件的關(guān)系，許多學(xué)者作了大量的研究，這些研究主要集中在獼猴桃生長季氣象條件、獼猴桃產(chǎn)量與氣象條件的關(guān)系[4-6]；獼猴桃種植氣候生態(tài)適宜性區(qū)劃[7]；獼猴桃高溫干旱災(zāi)害風(fēng)險(xiǎn)評估[8]；獼猴桃異常落果的原因及防御對策研究[9-12]等。近年來，較多研究還注重獼猴桃果品氣候品質(zhì)認(rèn)證模型和方法技術(shù)[13-15]。

適時(shí)采收是獼猴桃生產(chǎn)的關(guān)鍵技術(shù)之一。獼猴桃品質(zhì)與成熟度密切相關(guān)，如過早采收，風(fēng)味不佳。早采的獼猴桃流入市場，將大大損害獼猴桃在消費(fèi)者心目中的形象，造成品牌形象的損失。獼猴桃成熟期主要由當(dāng)?shù)貧夂驐l件決定，其次與品種有關(guān)。獼猴桃的產(chǎn)地分布較廣，品種各異，因此，成熟時(shí)間也相差較大。目前，獼猴桃的采收標(biāo)準(zhǔn)主要是依據(jù)獼猴桃生理成熟期，生產(chǎn)上最直接的判斷方法是借助折光分析儀測定其可溶性固形物含量來確定。一般中華獼猴桃可溶性固形物含量（G）達(dá)6.2%～6.5%，美味獼猴桃可溶性固形物含量（G）達(dá)6.6%～7.0%，就可以認(rèn)為達(dá)到采收成熟期。2013年以來，有些新培育的優(yōu)良品種可達(dá)8.0%以上[16]。

使用氣象資料對農(nóng)作物成熟期進(jìn)行預(yù)報(bào)，已較早在其他農(nóng)產(chǎn)品應(yīng)用[16-22]。而對獼猴桃成熟期預(yù)報(bào)的方法研究，目前還是空白。對于獼猴桃成熟期的判斷，目前仍然是主要依賴折光分析檢測儀，在缺乏檢測儀器的情況下，一般采用生育期天數(shù)判斷成熟期[23]，但同一品種，每年氣象條件差異造成每年生育期相差較大，少則相差10天左右，多則20天左右。因此，生育期天數(shù)判斷法，在實(shí)際生產(chǎn)中基本無法掌握，而廣大專業(yè)獼猴桃種植戶又缺少儀器檢測設(shè)備和檢測技術(shù)，開展獼猴桃成熟期預(yù)報(bào)方法研究，具有重要意義。

本研究采用湖北全苑生態(tài)農(nóng)業(yè)科技公司浠水獼猴桃種植基地2017—2019年不同時(shí)期采樣檢測的獼猴桃可溶性性固形物含量（G）資料，與氣象資料進(jìn)行相關(guān)分析，并求取成熟熵，以成熟熵作為當(dāng)年獼猴桃成熟期的預(yù)報(bào)因子，探討獼猴桃成熟期預(yù)報(bào)方法，為獼猴桃規(guī)模化人工種植基地，根據(jù)當(dāng)年氣象條件，準(zhǔn)確掌握成熟期，從而擺脫對折光分析檢測儀的依賴性，提供有效方法。

1材料和方法

1.1獼猴桃可溶性固形物檢測

獼猴桃可溶性性固形物含量資料來自湖北全苑生態(tài)農(nóng)業(yè)科技公司，資料時(shí)間2017—2019年。公司獼猴桃種植基地位于湖北省浠水縣散花鎮(zhèn)石牛山村，115.09°E，30.31°N，海拔高度45 m?；赜?013年開始規(guī)模化人工種植獼猴桃，種植品種主要為中科院武漢植物園培育的‘金玉’、‘金桃’和‘滿天紅’，按照生產(chǎn)有機(jī)產(chǎn)品的標(biāo)準(zhǔn)，通過土壤測定、抽槽整地、苗木定植、嫁接、棚架架設(shè)、配套建設(shè)滴灌、噴灌系統(tǒng)、防風(fēng)系統(tǒng)，種植面積100 hm2，2017年掛果面積約45 hm2?！鹩瘛?、‘金桃’和‘滿天紅’成熟采收標(biāo)準(zhǔn)根據(jù)品種屬性確定為可溶性固形物含量（G）≥8.0%。

2017年起，在獼猴桃生長果實(shí)營養(yǎng)積累期后期，即8月底或9月上中旬開始，每隔5～7天，對獼猴桃果實(shí)進(jìn)行成熟度測定，每次采樣數(shù)40～60果，求取平均。檢測項(xiàng)目含單果重、可溶性固形物含量（G）、干物質(zhì)含量、色彩角等，其中，可溶性固形物含量（G）檢測分別對花萼端及花柱端進(jìn)行檢測，再求取平均。

1.2氣象資料

氣象資料采用浠水縣國家氣象站觀測資料。浠水縣國家氣象站位于115.22°E，30.47°N，海拔高度61.8m，與獼猴桃種植基地直線距離18 km。資料時(shí)間與2017—2019年獼猴桃種植基地物候觀測相對應(yīng)。

2結(jié)果和分析

2.1獼猴桃成熟熵的設(shè)計(jì)與計(jì)算

先對獼猴桃不同時(shí)期測定的可溶性固形物含量與前階段不同時(shí)段發(fā)育期溫度、光照因子進(jìn)行相關(guān)分析，選取相關(guān)性最好的2個(gè)因子：獼猴桃末花期次日至檢測日前一天的≥0℃有效積溫（ΣT）和累計(jì)日照總時(shí)數(shù)（ΣS）。

對ΣT、ΣS進(jìn)行W檢驗(yàn)，W計(jì)算值分別為：2.0493、2.129，均>Wα（20，α=0.01） 0.868，表明數(shù)據(jù)符合正態(tài)分布。

對ΣT、ΣS進(jìn)行共線性檢驗(yàn)，R2=0.7672，方差膨脹因子VIFi=4.2955，表明ΣT、ΣS存在一定的共線性關(guān)系。使用逐步回歸法，擬合優(yōu)度變化顯著，計(jì)算F值也變化顯著，因此，忽略ΣT、ΣS共線性問題。

ΣT、ΣS與與可溶性固形物含量（G）相關(guān)系數(shù)分別為0.8577、0.8266，復(fù)相關(guān)系數(shù)為0.8759，F(xiàn)計(jì)算值26.9265>F查表值6.11 （α=0.001），通過F檢驗(yàn)。

為了更加直接和更簡便的使用氣象因子計(jì)算獼猴桃的成熟度，且又能統(tǒng)籌考慮所選氣象因子對作物（獼猴桃）成熟度的貢獻(xiàn)值大小，這里先設(shè)計(jì)獼猴桃成熟熵的計(jì)算方法并計(jì)算不同檢測時(shí)期成熟熵。

農(nóng)作物生長發(fā)育及正常成熟與光、溫、水等氣象條件密切相關(guān)，但，獼猴桃標(biāo)準(zhǔn)化種植基地在配套建設(shè)滴灌、噴灌設(shè)施的條件下，獼猴桃生長發(fā)育及果實(shí)成熟所需的水分條件能夠得到充分滿足，致使前期降水量與獼猴桃成熟度的相關(guān)性不明顯（計(jì)算失真），而溫度條件和光照條件仍然是依靠自然氣象條件，因此，這里不考慮降水量對成熟期的影響，只考慮溫度和光照2個(gè)條件。

以獼猴桃末花期次日至統(tǒng)計(jì)日間的≥0℃有效積溫（ΣT）和累計(jì)日照時(shí)數(shù)（ΣS）為獼猴桃成熟熵（H）的計(jì)算變量X1、X2。因積溫與日照時(shí)數(shù)單位和量綱均不同，需對2個(gè)變量進(jìn)行加權(quán)處理。用二元回歸分析計(jì)算得到的回歸系數(shù)b1=0.0016、b2=0.0041，分別計(jì)算b1、b2所占它們的總和的比重[24]，得到2個(gè)變量的權(quán)重系數(shù)分別為：0.28、0.72，用檢測日前一日ΣT和ΣS統(tǒng)計(jì)值分別與它們的權(quán)重系數(shù)相乘再求和，得到檢測日期的前一日的成熟熵（H），如式（1）所示。

H=0.28ΣT+0.72ΣS……………………………（1）

2.2獼猴桃成熟熵與可溶性固形物含量回歸分析

對各檢測日計(jì)算得到的獼猴桃成熟熵再次進(jìn)行W檢驗(yàn)，W計(jì)算值為2.230，通過W檢驗(yàn)，符合正態(tài)分布。

對成熟熵與可溶性固形物進(jìn)行直線回歸分析（表1），相關(guān)系數(shù)r=0.8718，并進(jìn)行T檢驗(yàn)，T計(jì)算值7.4747>查表值3.9920 （α=0.001），決定系數(shù)為0.76，相關(guān)顯著。

按表1成熟熵計(jì)算數(shù)據(jù)系列建立獼猴桃成熟熵直線回歸方程如（2）所示。

G=-3.9875+0.0060H…………………………（2）

圖1為獼猴桃可固含量（G）與成熟熵散點(diǎn)分布及回歸分析效果比較圖。

根據(jù)回歸方程計(jì)算得到的2017—2019年可溶性固形物含量值（G）與檢測實(shí)況值比較見表2。從結(jié)果看，兩者最小誤差0.03%，最大誤差-0.96%，平均誤差0.083%，均方差為0.31。

根據(jù)誤差分析正態(tài)分布原理，可溶性固形物含量（G）回歸計(jì)算值Gi置信度為0.05時(shí)的誤差范圍為：Gi±0.6076。

從表2分析誤差值情況，當(dāng)G≤8.0時(shí)，誤差值較小，誤差較大值主要是當(dāng)可溶性固形物含量（G）值>8.0%以后，即出現(xiàn)在達(dá)到采收期標(biāo)準(zhǔn)后。在成熟期預(yù)報(bào)應(yīng)用中，當(dāng)計(jì)算值G≥8.0時(shí)，已達(dá)到預(yù)報(bào)目的，停止計(jì)算，因此，這種較大誤差對回歸方程使用不產(chǎn)生實(shí)際影響。

按公式（1）計(jì)算得到3年成熟期的成熟熵分別為：1957.03、1961.43、1928.73，平均1949.06。按最大熵原則，取3個(gè)值中的最大值，即獼猴桃可溶性固形物含量（G）≥8.0%所需成熟熵為1961.43，其中，‘金玉’品種成熟熵為1961.43，‘金桃’品種為1928.73。

將G=8.0代入公式（2），計(jì)算得到的獼猴桃成熟熵為1948.38，計(jì)算值與3年實(shí)況平均值1949.06誤差較小（0.68），誤差率0.035%。

2.3獼猴桃成熟期預(yù)報(bào)與檢驗(yàn)

根據(jù)實(shí)測獼猴桃可溶性固形物含量（G）達(dá)到≥8.0%的前后相鄰兩個(gè)檢測日期內(nèi)插求得2017—2019年獼猴桃實(shí)際成熟期（實(shí)況）分別為：10月14日（‘金玉’）、9月28日（‘金玉’）、9月29日（‘金桃’）。

根據(jù)預(yù)報(bào)方程計(jì)算得到的2017—2019年獼猴桃成熟期日期分別為：10月15日（‘金玉’）、9月27日（‘金玉’）、10月1日（‘金桃’），與實(shí)況誤差-1～2天，平均誤差0.7（偏遲）天（見表3）。

3結(jié)論與討論

本研究基于收集到的獼猴桃可溶性固形物含量（G）檢測資料與氣象資料進(jìn)行相關(guān)分析，設(shè)計(jì)獼猴桃成熟熵計(jì)算方法，所選取的氣象因子：≥0℃有效積溫和累計(jì)日照時(shí)數(shù)，是作物生長發(fā)育及成熟的最重要?dú)庀笠蜃樱哂忻黠@的生物學(xué)意義。值得注意的是，在不具備灌溉設(shè)施條件下的獼猴桃人工種植基地，還需考慮將自然降水量納入成熟熵計(jì)算因子，對降水變量的加權(quán)處理方法相同。

本研究對積溫與日照數(shù)據(jù)的加權(quán)處理采用回歸系數(shù)權(quán)重法，充分考慮了2個(gè)氣象因子與獼猴桃成熟度的相關(guān)系數(shù)，消除了2個(gè)因子因單位和量綱的不同的影響，避免了人為主觀確定權(quán)重帶來的的分析誤差，計(jì)算結(jié)果與實(shí)況高度近似。

在獼猴桃果實(shí)生長接近成熟時(shí)，即在獼猴桃果實(shí)營養(yǎng)積累期后期，可隨時(shí)統(tǒng)計(jì)當(dāng)年獼猴桃末花期次日至當(dāng)前階段氣象資料，計(jì)算獼猴桃成熟熵，從而得到獼猴桃可溶性固形物含量（G），判斷當(dāng)前獼猴桃成熟度；也可按照獼猴桃不同品種可溶性固形物含量（G）采收標(biāo)準(zhǔn)，根據(jù)后期氣溫、日照數(shù)值預(yù)報(bào)，提前進(jìn)行獼猴桃成熟期預(yù)報(bào)，從而準(zhǔn)確預(yù)測獼猴桃商品上市時(shí)間。

參考文獻(xiàn)

[1]松會能.氣候條件對我國獼猴桃主要生態(tài)類型代表種地理分布的影響I—?dú)夂驐l件對我國獼猴桃水平分布的影響[J].貴州農(nóng)業(yè)科學(xué)，1987（2）：43-51，54.

[2]松會能.氣候條件對我國獼猴桃主要生態(tài)類型代表種地理分布的影響II—?dú)夂驐l件對我國獼猴桃垂直分布的影響[J].貴州農(nóng)業(yè)科學(xué)，1987（3）：40-48，30.

[3]姚春潮，張林森，劉旭峰.世界獼猴桃產(chǎn)業(yè)生產(chǎn)研究現(xiàn)狀[J].西北園藝（果樹?？?，2003（1）：54-55.

[4]宋云，劉中新，中國獼猴桃種植與氣象條件研究綜述[J].江蘇農(nóng)業(yè)科學(xué)，2020，48（8）：41-46.

[5]張舒姝，王艷妮，顏萍.都江堰市獼猴桃產(chǎn)量與氣象條件的關(guān)系[J].鄉(xiāng)村科技，2018（8）：57-58.

[6]池再香，張錦，李貴瓊，等.貴州紅心獼猴桃膨大期對氣溫與降水變化的響應(yīng)[J].氣象科技，2016，44（4）：669-674.

[7]賀文麗，李星敏，朱琳等.基于GIS的關(guān)中獼猴桃氣候生態(tài)適宜性區(qū)劃[J].中國農(nóng)學(xué)通報(bào)，2011，27（22）：202-207.

[8]王景紅，梁軼，柏秦鳳，等.陜西獼猴桃高溫干旱災(zāi)害風(fēng)險(xiǎn)區(qū)劃研究[J].中國農(nóng)學(xué)通報(bào)，2013，29（7）：105-110.

[9]葉茵，王天鏡.人工種植獼猴桃的氣候條件分析[J].貴州氣象， 2013，37（4）：41-44.

[10]屈振江，周廣勝.中國主栽獼猴桃品種的氣候適宜性區(qū)劃[J].中國農(nóng)業(yè)氣象，2017，38（4）：257-266.

[11]彭永宏，章文才.獼猴桃異常落果的原因及防御對策研究[J].中國農(nóng)業(yè)氣象，1994，15（2）：5-7.

[12]徐丹丹，黃桔梅，谷曉平，等.貴州獼猴桃人工栽培生態(tài)氣候條件調(diào)查[J].貴州農(nóng)業(yè)科學(xué)，2007，35（3）：30-31.

[13]劉璐，屈振江，張勇，等.陜西獼猴桃果品氣候品質(zhì)認(rèn)證模型構(gòu)建[J].陜西氣象，2017（4）：21-25.

[14]張向榮，何可杰，雷雯，等，獼猴桃果品氣候品質(zhì)認(rèn)證技術(shù)研究[J].陜西農(nóng)業(yè)科學(xué)，2015，61（10）：65-68.

[15]付芳婧，谷曉平，于飛.水城縣紅心獼猴桃氣候品質(zhì)認(rèn)證[J].耕作與栽培，2017（2）：47-50.

[16]陳素華，姚麗桃.內(nèi)蒙古自治區(qū)玉米成熟期預(yù)報(bào)方法的初步探討[J].內(nèi)蒙古氣象，1996（4）：24-27.

[17]張利華，孟錦，張永強(qiáng)，等.徐州地區(qū)冬小麥成熟期預(yù)報(bào)[J].江蘇農(nóng)業(yè)科學(xué)，2011，39（6）：157-159.

[18]楊棟，丁燁毅，金志鳳，等.浙江水蜜桃成熟期集合預(yù)報(bào)模型[J].中國農(nóng)業(yè)氣象，2018，39（6）：421-430.

[19]黃健熙，牛文豪，馬洪元，等，衛(wèi)星遙感和積溫-輻射模型預(yù)測區(qū)域冬小麥成熟期[J].農(nóng)業(yè)工程學(xué)報(bào)，2016，32（7）：152-157.

[20]段金省.隴東塬區(qū)冬小麥?zhǔn)斋@期預(yù)報(bào)[J].中國農(nóng)業(yè)氣象，2007，28（1）：17-20.

[21]李友勇，杜黎君，楊靖，等，基于特定時(shí)期溫度因子的小麥成熟期預(yù)測[J].中國農(nóng)業(yè)氣象，2012，33（1）：104-108.

[22]顧品強(qiáng)，姚瑤.黃桃始花期和成熟期統(tǒng)計(jì)預(yù)報(bào)模型研究[J].上海農(nóng)業(yè)學(xué)報(bào)，2013，29（3）：54-58.

[23]陜西省地方標(biāo)準(zhǔn)DB61《獼猴桃采收質(zhì)量與技術(shù)規(guī)范》.

[24]劉中新，羅和清，周汝寶，等.高溫脅迫對雜交中秈稻生育特性和產(chǎn)量的影響[J].湖北農(nóng)業(yè)科學(xué)，2018，57（12）：12-15.