油菜低温冻害天气指数保险研究
孔维财１，高 苹２，徐 敏２
（１．江苏省南京市高淳区气象局，江苏南京２１１３００；２．江苏省气候中心，江苏南京２１０００８）
摘要：利用１９６１—２０１９年油菜生育期内逐日气象资料、油菜产量和低温冻害数据，确定由最低气温构成的低温冻害天气指数指标，分析减产率和天气指数的关系，建立天气指数模型，使用燃烧定价法厘定了天气指数保险的纯保费率，并根据相关政策和实际情况，确定不同触发条件下的差异保费率、赔付标准，设计油菜低温冻害天气指数保险合同。低温冻害天气指数保险产品依托于相关国标、文献和实际情况，指标客观、操作方便，设计思路可为相似区域油菜低温冻害天气指数研究提供借鉴。
关键词：油菜；低温冻害；天气指数；保险；纯保费率；气象灾害
中图分类号：Ｆ８４０．６６ 文献标志码：Ａ 文章编号：１００２－１３０２（２０２１）０７－
０２４４－０５
收稿日期：２０２０－０７－２５
基金项目：国家重点研发计划（编号：２０１９ＹＦＤ１００２２０１）；南京市气象局自立科研项目（编号：
ＮＪ２０１８０２）。作者简介：孔维财（１９８６—），男，江苏南京人，硕士，工程师，主要从事农业气象研究。Ｅ－ｍａｉｌ：９４６９５９２７７＠ｑｑ．ｃｏｍ。
通信作者：徐 敏，硕士，高级工程师，主要从事农业气象研究。Ｅ－ｍａｉｌ：ａｍｉｎ０５０６＠１６３．ｃｏｍ。
我国油菜可分为春油菜和冬油菜，其中冬油菜
种植面积占绝大多数，且主要集中在长江流域［１］。
江苏省气象条件优越，适宜油菜生长，是我国重要的油菜种植区域。２０１９年江苏省油菜种植面积为
１４．４万ｈｍ２
，总产量为５０．５万ｔ，约占全国总产量的１３％［２］。近年来，受粮油价格上涨影响，农民种植
油菜的积极性提高，外加油菜花观光旅游经济的发展，种植面积不降反增。
天气指数保险是一种区别于传统的、基于损失赔付的新型保险产品，其赔付标准为保险合同中约定的天气指数数值或级别。天气指数保险能排除传统统计方法计算中存在的不确定性。灾害评估是研究天气指数保险的前提条件。我国学者已经做了大量相关研究，栾庆祖等研究了水果冰雹灾损
评估方法［３］
，金志凤等研究了浙江茶叶天气指数保险［４］，李德等研究了砀山酥梨花期低温冷害风险评估方法［５］，任义方等研究了江苏水稻高温热害天气指数保险［６］，金志凤等研究了浙江茶叶冻害评估方法［７］，这些方法的广泛应用，使天气指数的研究取
得了阶段性的成果。如何把农业灾害风险评估方法应用于天气指数保险是防灾减灾的重点之一。
杨太明等研究了安徽夏玉米干旱、水稻高温热害、
冬小麦种植天气指数保险［
８－１０］
，曹雯等研究了宁夏枸杞病害天气指数保险，这些农业灾害风险评估方法与天气指数保险相结合的产品得到了良好的推
广和应用，并取得了较好的服务效果［
１１］
。受全球气候变暖影响，江苏低温冻害频次增多、强度增强，低温冻害是油菜的重要自然灾害之一，油菜苗期低温、蕾薹期倒春寒、开花期低温均会造成油菜不同程度的减产。江苏省曾在１９６１年、１９６９年、１９７６年、１９８０年、１９９４年、２００８年、２０１８年出现过较严重的低温冻害，对油菜产量产生较大的影响。因此，开展提高油菜低温冻害的防灾减灾能力、灾后恢复生产的研究十分必要。农业保险是常用的风险转移的手段，常规农业保险在查险、定损、理赔、估价等方面存在较大分歧，基于天气指数保险能较好地解决这些问题。因不同区域气候的差异性，针对江苏省油菜天气指数保险的研究较少，本研究将为江苏省油菜低温冻害天气指数保险设计提供依据，对于拓展天气指数保险领域，规避油菜种植低温冻害风险很有必要。１ 资料与方法１．１ 资料来源
选取高淳国家气象站１９６１—２０１９年油菜播种到成熟期（根据生育期观测资料，确定多年平均
起止日期为９月２１日至翌年５月２２日）的气象数据；收集整理了南京市高淳区油菜种植面积、总产、单产资料和灾情信息，油菜数据来源于南京市统计年
鉴，灾情信息来源于《中国气象灾害大典（江苏卷）》和南京市气象局的历年灾情记录。
１．２ 减产率计算
产量主要受农业生产水平和气象条件的影响［１２］，可将油菜单产分解为随生产力变化的趋势产量、随气象条件变化的气象产量、随机变量３个部分，计算公式如下：
Ｙ＝Ｙ
ｗ＋Ｙ
ｔ
＋ε。（１）
式中，Ｙ为实际单产；Ｙ
ｗ为气象产量；Ｙ
ｔ
为趋势产
量，趋势产量由５年滑动平均值求得；ε为随机变量，由政策等因素导致，占比小，可忽略。
别那么骄傲2低温冻害对油菜产量的影响因地区、年份而异，为便于对比分析，本研究通过气象产量计算相对气象产量，把相对气象产量中负值部分定义为减产率，计算公式如下：
Ｓ＝
Ｙ
ｉ
－Ｙ
ｔ
Ｙ陈道明的第二任妻子
ｔ
×１００％ Ｙ
ｉ
＜Ｙ
ｔ
０Ｙ
ｉ
≥Ｙｔ
{。（２）式中：Ｓ为减产率；Ｙ
ｉ
为第ｉ年实际单产。
舞林大会走光１．３ 天气指数选取
油菜极易受苗期低温、蕾薹期倒春寒、开花期低温影响，且各生育期同一低温造成的灾害程度不一致［１３－１４］。油菜苗期生长最适温度为１０～２０℃，当最低气温降至－５～－３℃时，叶片轻度受冻；降至－９～－６℃时，叶片严重受冻，外围大叶基本冻死，低于－１０℃时心叶冻死，植株死亡。冬油菜一般在开春后，气温稳定在５℃以上时现蕾，气温在１０℃以上时可迅速抽薹，若蕾薹期气温低于－２℃将严重受冻，造成裂薹和死蕾；开花期遇到５℃以下低温，开花会明显减少［１５］。
天气保险指数的选取既要准确反映实际损失、受人为影响小，也要客观稳定、有效控制基差风险，在实际操作中要易于保险公司定损和理赔、易于投保人理解和推广［１６］。相关研究表明，油菜对低温敏感，减产与低温关系密切，且低温不受人为影响，便于投保人理解和操作，适合构建天气指数。本研究根据油菜对低温的敏感性、不同生育期低温致灾性，参考ＧＢ／Ｔ２７９５９—２０１１《南方水稻、油菜和柑桔低温灾害》《实用农业气象指标》《高淳县农业气候手册》，把低温冻害临界指标分为３级：１级冻害影响较轻，２级冻害影响中等，３级冻害严重、减产明显，每个等级有对应的最低气温范围，见表１。其中苗期分为３个等级，蕾薹期和开花期各设置１个等级。
表１ 油菜低温冻害临界指标
羽生结弦结婚了吗生育期低温冻害等级
温度范围
（℃）
平均起止日期
苗期 １（轻）－５～０９月２１日至翌年１月２９日
２（中）－９～－６
３（重）＜－１０
蕾薹期３（重）＜－２１月３０日至３月１７日
开花期３（重）＜５３月１８日至５月８日
油菜低温冻害不仅与极端最低气温相关，还与持续时间有关，在油菜生育期出现的冻害，造成的灾害一般呈现叠加的效应［１７］。为反映低温冻害的累计效应，综合考虑低温强度和持续时间对油菜的影响，设计低温冻害天气指数，公式如下：
Ｆ＝∑ｎ
ｊ＝１
Ｈ
ｊ
。（３）
式中：Ｆ为低温冻害天气指数；Ｈ
ｊ
为某年的第ｊ次低温冻害的灾害等级；ｎ为１个油菜生育期内低温冻害灾害出现的次数。
１．４ 纯保费率计算
纯保费率的计算是保险合同设计的重要组成部分，也是低温冻害风险的主要体现。天气指数纯保费率的计算一般采用燃烧定价法来计算，该方法是通过假设未来期望损失率与历史损失分布相同，通过历史数据来估算期望值的一种方法［１８－２０］。计算公式如下：
Ｒ＝
Ｅ［ｌｏｓｓ］
λμ
＝Ｅ［ｌｏｓｓ］＝∑ｎ
ｋ＝１
Ｓ
ｋ
Ｐ
ｋ
。（４）式中，Ｒ为纯保费率；Ｅ［ｌｏｓｓ］为期望损失率；λ为保障比例；μ为预期油菜单产比例；λ、μ均取值为
１００％［２１］；Ｓ
ｋ
为第ｋ级低温冻害的减产率；Ｐ
ｋ
为第ｋ级低温冻害的发生概率；ｋ为低温冻害等级。１．５ 天气指数保险合同设计
保险赔偿的触发条件及相应的赔付标准是油菜低温冻害天气指数保险合同的核心部分，其赔付公式如下：
Ｉ＝
Ｓ－Ｓ
ｍｉｎ
Ｓ
ｍａｘ
－Ｓ
()
ｍｉｎ
Ｑ Ｆ
ｍ
≥Ｆ０
０Ｆ
ｍ
＜Ｆ
{
０
。（５）
式中：Ｉ为单位面积赔付金额；Ｑ为保额；Ｓ
ｍａｘ
为最高
张嘉译资料减产率（历史最高减产率为４０．３％）；Ｓ
ｍｉｎ
为赔付触悼词范文
发值对应的减产率；Ｆ
ｍ
为第ｍ年天气指数；Ｆ
０
为保险赔付触发值。当天气指数大于或等于触发值时，投保人可获得保险赔偿。投保人可以根据历史灾情及实际情况，选择适合的触发值条件进行投保，
以获得更好的经济效益。２ 结果与分析
２．１ 油菜产量及低温变化规律分析
随着科学技术的发展、农业机械化的推进、优良种子的推广和政策的支持，油菜产量呈现波动上升的
趋势，线性趋势倾向率为４４５．４ｋｇ／（ｈｍ２·１０年）。
近５
９年油菜产量经历了３个明显的上升阶段（图１），分别为２０世纪６０年代到７０年代末、２０世纪８０年代、２０世纪９０年代到２１世纪１０年代末；近５９年油菜产量有３个明显的波谷，当实际产量低于趋势产量时为减产年份，１９７７年、１９８８年、２０１１年减产明显，分析其原因为１９７７年出现了－１２．２℃的极端最低气温，０℃以下低温有５２ｄ，低温冻害严重，减产明显；
１９８８年油菜生育期内极端最低气温为－６．３℃，０℃以下低温有３０ｄ，减产率最小；２０１１年油菜生育期内极端最低气温为－８．３℃，０℃以下低温有５１ｄ，减产率小于１９７７
年。
气温是影响油菜生长发育和产量形成的重要气象因子。南京市高淳区１
９６１—２０１９年年极端最低气温呈升高的趋势（图２），气候倾向率为
０．５７℃／１０年（ｒ２＝０．１８７，Ｐ＜０．００１），历史极端最
低气温为－１４．０℃，出现在１９６９年，年最低气温最高值为－
１．５℃，出现在２０１９年。油菜生育期低温持续时间对产量也有很大的影响。油菜生育期低温天数呈现波动下降的趋势（图３），其中＜５℃天数下降趋势最大，线性趋势倾向率为－４．７ｄ／１０年，下降趋势最小的为＜－５℃天数，线性趋势倾向率为－１．０ｄ／１０年。１９７７年＜－５℃天数最多，此年油菜出现明显减产，与图１中产量下降的趋势吻合；１９６８年＜０℃天数最多，为７２ｄ，但１９６８年极端最低气温相对较高，为－５．６℃，油菜产量未出现明
显下降，
２００１年＜０℃天数最少，为８ｄ；１９７０年＜５℃天数最多，为１２４ｄ，此年极端最低气温为－１０．１℃，油菜出现小幅减产。年极端最低气温的升高趋势及油菜生育期低温天数下降趋势可能与
全球的气候变暖背景有关。
２．２ 天气指数模型构建
由于天气指数保险存在基差风险，本研究根据《中国气象灾害大典（江苏卷）》和南京市气象局灾情记录对油菜样本进行筛选，去除非低温冻害导致
油菜减产的年份。根据公式（３）计算低温冻害天气指数，结果（图４）表明，近５９年以来，油菜低温冻害天气指数最大值为１
０３，出现在１９６８年；最小值为１０，出现在２０１７年，其中７５％的天气指数在２０～６０之间，Ｆ≤２０出现的概率为１０％，Ｆ≥８０出现的概率为５％，极值的出现概率较小。
为建立油菜低温冻害天气指数模型，采用油菜减产率与天气指数进行回归分析（图５），得到减产率与低温冻害天气指数的关系模型，如公式（６）所示，减产率与低温冻害天气指数相关系数０．３０，通过了０
．０５显著性水平检验。Ｓ＝０．１５９Ｆ＋３．６８５。
（６）
根据天气指数的影响程度，将天气指数分为５级（表２），Ⅱ、Ⅲ级发生概率为７
７．２％；Ⅰ级发生概
率较小，且对油菜生长影响较小，平均减产率为５ ４％；Ⅳ级影响较大，平均减产率为１４．９％，但发生概率较低，概率仅为９．１％；Ⅴ级对油菜影响极大，但发生概率最小，仅为４．６％。由此可见，低温冻害一般会造成油菜减产３．８％～１９．６％，这与实际生产较一致，因此，本研究构建的天气指数合理可行。
表２ 油菜低温冻害天气指数与减产率的关系
等级Ｆ取值范围发生概率
（％）
减产率
（％）
平均减产率
（％）
Ⅰ１≤Ｆ≤２０９．１３．８～６．９５．４
Ⅱ２１≤Ｆ≤４０３６．３７．０～１０．０８．５
Ⅲ４１≤Ｆ≤６０４０．９１０．１～１３．２１１．７Ⅳ６１≤Ｆ≤８０９．１１３．３～１６．４１４．９Ⅴ８１≤Ｆ≤１００４．６１６．５～１９．６１８．１
２．３ 天气指数合同设计
保额一般通过生产成本或产量确定，一般为最近连续５年平均产量的４０％～８０％。参照江苏农业保险网（ｈｔｔｐ：／／ｗｗｗ．ｊｓｎｘ．ｏｒｇ／）及相关研究，本研究将２０１５—２０１９年的平均产量、油菜最新平均单价和历史最大减产率相乘获得保额。油菜历史最高减产率为４０．３％，按照现在油菜的市场价计算，保额为１０５９１．８元／ｈｍ２，这与《江苏省政策性农业保险油菜种植保险条款》《省政府办公厅关于做好２０１４年全省农业保险工作的通知》［苏政发（２０１４）２１号］中规定的１ｈｍ２油菜保险金额３个档次（６０００、８２５０、１０５００元）中第３个档次十分吻合。
免赔额是指在保险合同中规定的损失在一定限度内保险人不负赔偿责任的额度，主要是为了减少频繁
发生的小额赔付，提高被保险人的责任心和注意力，降低保险机构的经营成本，按照相关规定，本研究定义损失率１０％为免赔额，代入公式（５），计算出赔付标准，见表３。
表３ 油菜低温冻害天气指数赔付方案
天气指数
取值范围
赔付比率
（％）
赔付金额
（元／ｈｍ２）Ｆ≤２００．００．０
２１≤Ｆ≤４０２４．４２５８４．４
４１≤Ｆ≤６０４８．２５１０５．２
６１≤Ｆ≤８０７２．０７６２６．１
８１≤Ｆ≤１００９５．８１０１４６．９
Ｆ＞１００１００．０１０５９１．８２．４ 天气指数保险费率的厘定
基于分析１９６１—２０１９年油菜生育期气象资料，得出低温冻害天气指数与灾害发生概率，代入公式（４）计算出不同触发条件下纯保费率和保费，见表
４。当触发值Ｆ
０
为８１时，保费最低为８４．７元／ｈｍ２，在历史数据中只触发过１次。相比较而言，触发值
Ｆ
０
为４１或６１时，灾害发生概率大、获得赔付可能性高，且保费少、经济效益高。投保人可根据实际情况，选择合适的触发条件，以便在低温冻害中获得更大收益。
３ 结论与讨论
本研究在国内外学者研究成果的基础上，通过分析１９６１—２０１９年油菜生育期气象数据、油菜产量
表４ 油菜低温冻害天气指数不同触发条件下纯保费率和保费
Ｆ０取值纯保费率
（％）
保费（元／ｈｍ２）
２１１０．４１１０１．５
４１６．８７２０．２
６１２．１２２２．４
８１０．８８４．７
和低温冻害灾害数据，针对油菜生育期灾害分析了减产率和低温冻害天气指数的关系，建立了低温冻害天气指数模型，线性模型相关系数０．３０，通过了０．０５水平显著性检验，厘定了天气指数保险的赔付标准，设计了天气指数保险产品，以期为江苏省保险公司开展油菜天气指数保险提供科学的依据。
本研究在确定油菜低温冻害指标时参考了农业气象专家建议和实际工作经验，也查阅了相关国标、气象指标书籍；天气指数选择是依据长时间序列的气象与产量数据，避免了因数据少而导致的不稳定现象；天气指数是依据最低气温建立的，也考虑了持续低温的叠加效应，能综合反映油菜低温冻害的情况，在实际操作中，方便投保人与保险公司进行投保和赔付，操作方便简单的同时也减少了争议的发生；本指数保险产品以１０％为免赔额，设计了不同出发条件下纯保费率和保费，纯保费率为０ ８％～１０．４％，保费为８４．７～１１０１．５元／ｈｍ２，方便投保人根据实际情况选择适合的保险触发条件，以获得更大收益。
２０１８年全国性特大雪灾持续时间长、降水量大、最低温度低，油菜低温冻害严重，据相关统计，南京市高淳区油菜赔付金额达２２１８．０万元，受此影响，２０１９年高淳油菜投保面积有所上升，为２．６６万ｈｍ２，仍不到总种植面积的一半。在实际工作中发现，高淳区一些油菜种植大户参与积极性较高，但大多数地区仍以农户为种植单位，且种植地域不集中、管理技术水平差异较大，造成了油菜保险推广困难，在以后的工作中，有待结合区域自动站数据开发精细化天气指数，寻保险的突破口，以更好地适应农业保险的实际需要。
油菜生育期低温冻害不仅与低温强度和持续时间有关，还与１ｄ中低温出现时段、降温幅度、降水量、日照时间等有关，为方便计算和应用，本研究只选择低温及持续时间作为天气指数指标，同时，文中天气指数保险产品只适用于单个灾害，而在油菜生育期中，多种气象灾害会起到叠加效果，在以
后的研究中，有待研发更全面且具符合生物学意义的天气指数。
参考文献：
［１］鞠英芹，杨霏云，马德栗，等．江淮地区油菜渍害的时空分析［Ｊ］．自然灾害学报，２０１７，２６（６）：１３６－１４６．
［２］张 佩，吴 田，姚 薇，等．江苏省油菜全生育期气候特征的变化分析［Ｊ］．江苏农业科学，２０１５，４３（９）：１１８－１２１．
［３］栾庆祖，董鹏捷，叶彩华．面向气象指数保险的水果冰雹灾害灾损评估方法［Ｊ］．中国农业气象，２０１９，４０（６）：４０２－４１０．
［４］金志凤，叶建刚，杨再强，等．浙江省茶叶生长的气候适宜性［Ｊ］．应用生态学报，２０１４，２５（４）：９６７－９７３．
［５］李 德，杨太明，戚尚恩，等．黄河故道砀山酥梨花期低温冷害风险的评估方法初探［Ｊ］．中国农业气象，２００９，３０（４）：６１１－６１５．［６］任义方，赵艳霞，张旭晖，等．江苏水稻高温热害气象指数保险风险综合区划［Ｊ］．中国农业气象，２０１９，４０（６）：３９１－４０１．
［７］金志凤，胡 波，严甲真，等．浙江省茶叶农业气象灾害风险评价［Ｊ］．生态学杂志，２０１４，３３（３）：７７１－７７７．
［８］杨太明，许 莹，孙喜波，等．安徽省夏玉米干旱天气指数保险产品设计及应用［Ｊ］．气象，２０１６，４２（４）：４５０－４５５．
［９］杨太明，孙喜波，刘布春，等．安徽省水稻高温热害保险天气指数模型设计［Ｊ］．中国农业气象，２０１５，３６（２）：２２０－２２６．
［１０］杨太明，刘布春，孙喜波，等．安徽省冬小麦种植保险天气指数设计与应用［Ｊ］．中国农业气象，２０１３，３４（２）：２２９－２３５．
［１１］曹 雯，武万里，杨太明，等．宁夏枸杞病害天气指数保险研究［Ｊ］．干旱气象，２０１９，３７（５）：８５７－８６５．
［１２］徐 敏，徐经纬，高 苹，等．不同统计模型在冬小麦产量预报中的预报能力评估———以江苏麦区为例［Ｊ］．中国生态农业学报（中英文），２０２０，２８（３）：４３８－４４７．
［１３］张学昆，张春雷，廖 星，等．２００８年长江流域油菜低温冻害调查分析［Ｊ］．中国油料作物学报，２００８，３０（１）：１２２－１２６．
［１４］ＸｕＳＪ，ＣｈｏｎｇＫ．Ｒｅｍｅｍｂｅｒｉｎｇｗｉｎｔｅｒｔｈｒｏｕｇｈｖｅｒｎａｌｉｓａｔｉｏｎ［Ｊ］．ＮａｔｕｒｅＰｌａｎｔｓ，２０１８，４（１２）：９９７－１００９．
［１５］杨霏云．实用农业气象指标［Ｍ］．北京：气象出版社，２０１５．［１６］李懿珈，叶 涛，德庆卓嘎，等．基于遥感植被指数的西藏那曲地区畜牧业旱灾指数保险产品设计研究［Ｊ］．农业现代化研究，２０１８，３９（４）：６８０－６８８．
［１７］周 嘉，王小巍，周安宁，等．平凉市苹果花期冻害政策性农业保险气象理赔指数设计［Ｊ］．湖北农业科学，２０１７，５６（１６）：３１４４－３１４８．
［１８］王新伟，杜明哲，王 丽，等．河南省花生连阴雨灾害气象指数保险设计［Ｊ］．生态学杂志，２０１８，３７（１１）：３３９０－３３９５．
［１９］ＺｈｏｕＬＨ，ＹａｎＴ，ＣｈｅｎＸ，ｅｔａｌ．Ｅｆｆｅｃｔｏｆｈｉｇｈｎｉｇｈｔｔｅｍｐｅｒａｔｕｒｅｏｎｓｔｏｒａｇｅｌｉｐｉｄｓａｎｄｔｒａｎｓｃｒｉｐｔｏｍｅｃｈａｎｇｅｓｉｎｄｅｖｅｌｏｐｉｎｇｓｅｅｄｓｏｆｏｉｌｓｅｅｄｒａｐｅ［Ｊ］．ＪｏｕｒｎａｌｏｆＥｘｐｅｒｉｍｅｎｔａｌＢｏｔａｎｙ，２０１８，６９（７）：１７２１－１７３３．
［２０］ＧｕＸＢ，ＬｉＹＮ，ＤｕＹＤ．Ｂｉｏｄｅｇｒａｄａｂｌｅｆｉｌｍｍｕｌｃｈｉｎｇｉｍｐｒｏｖｅｓｓｏｉｌｔｅｍｐｅｒａｔｕｒｅ，ｍｏｉｓｔｕｒｅａｎｄｓｅｅｄｙｉｅｌｄｏｆｗｉｎｔｅｒｏｉｌｓｅｅｄｒａｐｅ（ＢｒａｓｓｉｃａｎａｐｕｓＬ．）［Ｊ］．Ｓｏｉｌ＆ＴｉｌｌａｇｅＲｅｓｅａｒｃｈ，２０１７，１７１：４２－５０．［２１］郑小琴，赖焕雄，徐宗焕．台湾热带优良水果（寒）冻害气象保险指数设计［Ｊ］．西南农业学报，２０１１，２４（４）：１５９８－１６０３．
发布评论