刘中新 宋云

摘要：根據收集到的猕猴桃果实生长后期采集的猕猴桃可溶性固形物含量检测数据资料，使用前期气象资料设计并计算猕猴桃成熟熵。选取猕猴桃当年末花期次日至检测日前一日≥0℃积温和累计日照时数作为计算猕猴桃成熟熵的气象因子，按其与猕猴桃可固含量的回归系数权重加权求和，得到成熟熵指数（H）的计算公式为：H=0.28ΣT+0.72ΣS，对猕猴桃成熟熵与可溶性固形物检测含量进行回归分析，建立可溶性固形物检测含量与猕猴桃成熟熵回归方程：G=-3.9875+0.0060H。根据回归方程计算得到的可溶性固形物含量值（G）与检测值最小误差0.03%，最大误差-0.96%，平均误差0.083%。参照猕猴桃成熟要求可溶性固形物含量≥8.0%标准，计算得到猕猴桃所需成熟熵为1949.06，预报2017—2019年猕猴桃成熟期与实况误差-1～2天，平均误差0.7天（偏迟）。运用本方法可根据气象因子计算当期猕猴桃可溶性固形物含量，结合后期气象要素预报，可提前计算预报猕猴桃成熟期，籍此可以摆脱对折光分析检测仪器的依赖，又能提前预报猕猴桃成熟期和上市期，为避免早采而造成对猕猴桃品牌形象的损害提供科学依据。

关键词：猕猴桃；成熟期；预报；可溶性固形物；成熟熵；积温；日照

中图分类号：S165+.26文献标志码：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），也称狐狸桃、羊桃、毛木果、奇异果等，因猕猴喜食，故名猕猴桃，是一种品质鲜嫩，营养丰富，风味鲜美的水果。猕猴桃是20世纪人工驯化栽培野生果树最有成就的四大果种之一[1-3]，其浆果具有特殊的风味，含多种氨基酸、多种矿物质、维生素C，藤、根、叶及果均可入药，具有很高的经济价值和使用价值。

猕猴桃现代种植技术已不断发展，人工规模化种植面积不断扩大。有关猕猴桃种植技术及与气象条件的关系，许多学者作了大量的研究，这些研究主要集中在猕猴桃生长季气象条件、猕猴桃产量与气象条件的关系[4-6]；猕猴桃种植气候生态适宜性区划[7]；猕猴桃高温干旱灾害风险评估[8]；猕猴桃异常落果的原因及防御对策研究[9-12]等。近年来，较多研究还注重猕猴桃果品气候品质认证模型和方法技术[13-15]。

适时采收是猕猴桃生产的关键技术之一。猕猴桃品质与成熟度密切相关，如过早采收，风味不佳。早采的猕猴桃流入市场，将大大损害猕猴桃在消费者心目中的形象，造成品牌形象的损失。猕猴桃成熟期主要由当地气候条件决定，其次与品种有关。猕猴桃的产地分布较广，品种各异，因此，成熟时间也相差较大。目前，猕猴桃的采收标准主要是依据猕猴桃生理成熟期，生产上最直接的判断方法是借助折光分析仪测定其可溶性固形物含量来确定。一般中华猕猴桃可溶性固形物含量（G）达6.2%～6.5%，美味猕猴桃可溶性固形物含量（G）达6.6%～7.0%，就可以认为达到采收成熟期。2013年以来，有些新培育的优良品种可达8.0%以上[16]。

使用气象资料对农作物成熟期进行预报，已较早在其他农产品应用[16-22]。而对猕猴桃成熟期预报的方法研究，目前还是空白。对于猕猴桃成熟期的判断，目前仍然是主要依赖折光分析检测仪，在缺乏检测仪器的情况下，一般采用生育期天数判断成熟期[23]，但同一品种，每年气象条件差异造成每年生育期相差较大，少则相差10天左右，多则20天左右。因此，生育期天数判断法，在实际生产中基本无法掌握，而广大专业猕猴桃种植户又缺少仪器检测设备和检测技术，开展猕猴桃成熟期预报方法研究，具有重要意义。

本研究采用湖北全苑生态农业科技公司浠水猕猴桃种植基地2017—2019年不同时期采样检测的猕猴桃可溶性性固形物含量（G）资料，与气象资料进行相关分析，并求取成熟熵，以成熟熵作为当年猕猴桃成熟期的预报因子，探讨猕猴桃成熟期预报方法，为猕猴桃规模化人工种植基地，根据当年气象条件，准确掌握成熟期，从而摆脱对折光分析检测仪的依赖性，提供有效方法。

1材料和方法

1.1猕猴桃可溶性固形物检测

猕猴桃可溶性性固形物含量资料来自湖北全苑生态农业科技公司，资料时间2017—2019年。公司猕猴桃种植基地位于湖北省浠水县散花镇石牛山村，115.09°E，30.31°N，海拔高度45 m。基地于2013年开始规模化人工种植猕猴桃，种植品种主要为中科院武汉植物园培育的‘金玉’、‘金桃’和‘满天红’，按照生产有机产品的标准，通过土壤测定、抽槽整地、苗木定植、嫁接、棚架架设、配套建设滴灌、喷灌系统、防风系统，种植面积100 hm2，2017年挂果面积约45 hm2。‘金玉’、‘金桃’和‘满天红’成熟采收标准根据品种属性确定为可溶性固形物含量（G）≥8.0%。

2017年起，在猕猴桃生长果实营养积累期后期，即8月底或9月上中旬開始，每隔5～7天，对猕猴桃果实进行成熟度测定，每次采样数40～60果，求取平均。检测项目含单果重、可溶性固形物含量（G）、干物质含量、色彩角等，其中，可溶性固形物含量（G）检测分别对花萼端及花柱端进行检测，再求取平均。

1.2气象资料

气象资料采用浠水县国家气象站观测资料。浠水县国家气象站位于115.22°E，30.47°N，海拔高度61.8m，与猕猴桃种植基地直线距离18 km。资料时间与2017—2019年猕猴桃种植基地物候观测相对应。

2结果和分析

2.1猕猴桃成熟熵的设计与计算

先对猕猴桃不同时期测定的可溶性固形物含量与前阶段不同时段发育期温度、光照因子进行相关分析，选取相关性最好的2个因子：猕猴桃末花期次日至检测日前一天的≥0℃有效积温（ΣT）和累计日照总时数（ΣS）。

对ΣT、ΣS进行W检验，W计算值分别为：2.0493、2.129，均>Wα（20，α=0.01） 0.868，表明数据符合正态分布。

对ΣT、ΣS进行共线性检验，R2=0.7672，方差膨胀因子VIFi=4.2955，表明ΣT、ΣS存在一定的共线性关系。使用逐步回归法，拟合优度变化显著，计算F值也变化显著，因此，忽略ΣT、ΣS共线性问题。

ΣT、ΣS与与可溶性固形物含量（G）相关系数分别为0.8577、0.8266，复相关系数为0.8759，F计算值26.9265>F查表值6.11 （α=0.001），通过F检验。

为了更加直接和更简便的使用气象因子计算猕猴桃的成熟度，且又能统筹考虑所选气象因子对作物（猕猴桃）成熟度的贡献值大小，这里先设计猕猴桃成熟熵的计算方法并计算不同检测时期成熟熵。

农作物生长发育及正常成熟与光、温、水等气象条件密切相关，但，猕猴桃标准化种植基地在配套建设滴灌、喷灌设施的条件下，猕猴桃生长发育及果实成熟所需的水分条件能够得到充分满足，致使前期降水量与猕猴桃成熟度的相关性不明显（计算失真），而温度条件和光照条件仍然是依靠自然气象条件，因此，这里不考虑降水量对成熟期的影响，只考虑温度和光照2个条件。

以猕猴桃末花期次日至统计日间的≥0℃有效积温（ΣT）和累计日照时数（ΣS）为猕猴桃成熟熵（H）的计算变量X1、X2。因积温与日照时数单位和量纲均不同，需对2个变量进行加权处理。用二元回归分析计算得到的回归系数b1=0.0016、b2=0.0041，分别计算b1、b2所占它们的总和的比重[24]，得到2个变量的权重系数分别为：0.28、0.72，用检测日前一日ΣT和ΣS统计值分别与它们的权重系数相乘再求和，得到检测日期的前一日的成熟熵（H），如式（1）所示。

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

2.2猕猴桃成熟熵与可溶性固形物含量回归分析

对各检测日计算得到的猕猴桃成熟熵再次进行W检验，W计算值为2.230，通过W检验，符合正态分布。

对成熟熵与可溶性固形物进行直线回归分析（表1），相关系数r=0.8718，并进行T检验，T计算值7.4747>查表值3.9920 （α=0.001），决定系数为0.76，相关显著。

按表1成熟熵计算数据系列建立猕猴桃成熟熵直线回归方程如（2）所示。

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

图1为猕猴桃可固含量（G）与成熟熵散点分布及回归分析效果比较图。

根据回归方程计算得到的2017—2019年可溶性固形物含量值（G）与检测实况值比较见表2。从结果看，两者最小误差0.03%，最大误差-0.96%，平均误差0.083%，均方差为0.31。

根据误差分析正态分布原理，可溶性固形物含量（G）回归计算值Gi置信度为0.05时的误差范围为：Gi±0.6076。

从表2分析误差值情况，当G≤8.0时，误差值较小，误差较大值主要是当可溶性固形物含量（G）值>8.0%以后，即出现在达到采收期标准后。在成熟期预报应用中，当计算值G≥8.0时，已达到预报目的，停止计算，因此，这种较大误差对回归方程使用不产生实际影响。

按公式（1）计算得到3年成熟期的成熟熵分别为：1957.03、1961.43、1928.73，平均1949.06。按最大熵原则，取3个值中的最大值，即猕猴桃可溶性固形物含量（G）≥8.0%所需成熟熵为1961.43，其中，‘金玉’品种成熟熵为1961.43，‘金桃’品种为1928.73。

将G=8.0代入公式（2），计算得到的猕猴桃成熟熵为1948.38，计算值与3年实况平均值1949.06误差较小（0.68），误差率0.035%。

2.3猕猴桃成熟期预报与检验

根据实测猕猴桃可溶性固形物含量（G）达到≥8.0%的前后相邻两个检测日期内插求得2017—2019年猕猴桃实际成熟期（实况）分别为：10月14日（‘金玉’）、9月28日（‘金玉’）、9月29日（‘金桃’）。

根据预报方程计算得到的2017—2019年猕猴桃成熟期日期分别为：10月15日（‘金玉’）、9月27日（‘金玉’）、10月1日（‘金桃’），与实况误差-1～2天，平均误差0.7（偏迟）天（见表3）。

3结论与讨论

本研究基于收集到的猕猴桃可溶性固形物含量（G）檢测资料与气象资料进行相关分析，设计猕猴桃成熟熵计算方法，所选取的气象因子：≥0℃有效积温和累计日照时数，是作物生长发育及成熟的最重要气象因子，具有明显的生物学意义。值得注意的是，在不具备灌溉设施条件下的猕猴桃人工种植基地，还需考虑将自然降水量纳入成熟熵计算因子，对降水变量的加权处理方法相同。

本研究对积温与日照数据的加权处理采用回归系数权重法，充分考虑了2个气象因子与猕猴桃成熟度的相关系数，消除了2个因子因单位和量纲的不同的影响，避免了人为主观确定权重带来的的分析误差，计算结果与实况高度近似。

在猕猴桃果实生长接近成熟时，即在猕猴桃果实营养积累期后期，可随时统计当年猕猴桃末花期次日至当前阶段气象资料，计算猕猴桃成熟熵，从而得到猕猴桃可溶性固形物含量（G），判断当前猕猴桃成熟度；也可按照猕猴桃不同品种可溶性固形物含量（G）采收标准，根据后期气温、日照数值预报，提前进行猕猴桃成熟期预报，从而准确预测猕猴桃商品上市时间。

参考文献

[1]松会能.气候条件对我国猕猴桃主要生态类型代表种地理分布的影响I—气候条件对我国猕猴桃水平分布的影响[J].贵州农业科学，1987（2）：43-51，54.

[2]松会能.气候条件对我国猕猴桃主要生态类型代表种地理分布的影响II—气候条件对我国猕猴桃垂直分布的影响[J].贵州农业科学，1987（3）：40-48，30.

[3]姚春潮，张林森，刘旭峰.世界猕猴桃产业生产研究现状[J].西北园艺（果树专刊），2003（1）：54-55.

[4]宋云，刘中新，中国猕猴桃种植与气象条件研究综述[J].江苏农业科学，2020，48（8）：41-46.

[5]张舒姝，王艳妮，颜萍.都江堰市猕猴桃产量与气象条件的关系[J].乡村科技，2018（8）：57-58.

[6]池再香，张锦，李贵琼，等.贵州红心猕猴桃膨大期对气温与降水变化的响应[J].气象科技，2016，44（4）：669-674.

[7]贺文丽，李星敏，朱琳等.基于GIS的关中猕猴桃气候生态适宜性区划[J].中国农学通报，2011，27（22）：202-207.

[8]王景红，梁轶，柏秦凤，等.陕西猕猴桃高温干旱灾害风险区划研究[J].中国农学通报，2013，29（7）：105-110.

[9]叶茵，王天镜.人工种植猕猴桃的气候条件分析[J].贵州气象， 2013，37（4）：41-44.

[10]屈振江，周广胜.中国主栽猕猴桃品种的气候适宜性区划[J].中国农业气象，2017，38（4）：257-266.

[11]彭永宏，章文才.猕猴桃异常落果的原因及防御对策研究[J].中国农业气象，1994，15（2）：5-7.

[12]徐丹丹，黄桔梅，谷晓平，等.贵州猕猴桃人工栽培生态气候条件调查[J].贵州农业科学，2007，35（3）：30-31.

[13]刘璐，屈振江，张勇，等.陕西猕猴桃果品气候品质认证模型构建[J].陕西气象，2017（4）：21-25.

[14]张向荣，何可杰，雷雯，等，猕猴桃果品气候品质认证技术研究[J].陕西农业科学，2015，61（10）：65-68.

[15]付芳婧，谷晓平，于飛.水城县红心猕猴桃气候品质认证[J].耕作与栽培，2017（2）：47-50.

[16]陈素华，姚丽桃.内蒙古自治区玉米成熟期预报方法的初步探讨[J].内蒙古气象，1996（4）：24-27.

[17]张利华，孟锦，张永强，等.徐州地区冬小麦成熟期预报[J].江苏农业科学，2011，39（6）：157-159.

[18]杨栋，丁烨毅，金志凤，等.浙江水蜜桃成熟期集合预报模型[J].中国农业气象，2018，39（6）：421-430.

[19]黄健熙，牛文豪，马洪元，等，卫星遥感和积温-辐射模型预测区域冬小麦成熟期[J].农业工程学报，2016，32（7）：152-157.

[20]段金省.陇东塬区冬小麦收获期预报[J].中国农业气象，2007，28（1）：17-20.

[21]李友勇，杜黎君，杨靖，等，基于特定时期温度因子的小麦成熟期预测[J].中国农业气象，2012，33（1）：104-108.

[22]顾品强，姚瑶.黄桃始花期和成熟期统计预报模型研究[J].上海农业学报，2013，29（3）：54-58.

[23]陕西省地方标准DB61《猕猴桃采收质量与技术规范》.

[24]刘中新，罗和清，周汝宝，等.高温胁迫对杂交中籼稻生育特性和产量的影响[J].湖北农业科学，2018，57（12）：12-15.