设施葡萄低温灾害天气指数保险纯费率的确定方法
【技术领域】
[0001] 本发明涉及一种保险费率的计算，具体涉及设施葡萄低温灾害天气指数保险纯费 率的确定方法。
【背景技术】
[0002] 葡萄（K/iis 属落叶藤本植物，有"水果之神"的称号，是世界上最主要 水果之一，世界葡萄产量约占水果总产量的1/4。近年来，中国葡萄设施栽培发展十分迅速， 2014年全国葡萄栽培面积已达到300余万亩，江苏是葡萄栽培主要区域，栽培面积达45万 亩。其中设施栽培面积达到10万亩。由于江苏地处副热带高压季风气候区，夏季高温、高 湿和冬春季低温、寡照等不利气象条件严重影响设施葡萄生长及品质，因此，如何有效开展 设施葡萄气象服务、有效防御气象灾害成为气象部门和生产部门面临的技术问题。目前，关 于设施葡萄气象灾害监测预警和保障服务尚处于起步阶段，气象服务还不能满足设施葡萄 安全生产和快速发展的需要。
[0003] 国内外关于气象条件及气象灾害对葡萄生长和品质的影响有一定报道。欧洲和 南炜20°~40°之间是世界上葡萄栽培主要区域，前人研究表明高温、病害和缺乏足够的 温度诱发葡萄的休眠是限制葡萄生长发育的因素，高温、低温、寡照也是影响葡萄品质的主 要气候因子（Wen，2008 ;UlienChuche&DenisThi6ry，2009)。冬天低温潮湿条件能够诱发 葡萄病害发生，如低温寒潮和霜冻导致摩尔达维亚（罗马尼亚）葡萄顶芽、枝的严重伤害 (01ivier，2012)。Maurer等（2011)通过分析1523-2007年澳大利亚、法国、瑞士的历史气 象数据与产量数据，认为葡萄产量与春季和夏季平均气温成显著相关正相关。有学者研究 认为西班牙地中海气候区的温度和降水是影响葡萄收获期和产量的主要因素（Josepetal， 2012;Lorenzo, 2013)。Nyamdorj等（2014)研究澳大利亚西部气候变化对葡萄果实品质 (花青素和可滴定酸）的影响。刘明春等（2007)研究河西干旱区酿酒葡萄生长的气象条件， 果粒增长的适宜气温为20-21°C，超过21°C，增速明显变缓，40°C以上的高温将使叶片变黄 变褐，果实日灼。王瑞华（2013)研究也证实葡萄叶片脯氨酸的含量在高温胁迫下明显高于 对照，而叶绿素的含量和可溶性糖的含量则随胁迫时间的延长呈现低-高-低的波动，高温 高湿天气容易发生炭疽病、白腐病危害(汤照云等，2006)。
【发明内容】
[0004] 本发明的目的是为了解决现有技术中存在的缺陷，提供一种设施葡萄低温灾害天 气指数保险纯费率的计算，为设施葡萄的低温灾害风险管理提供了科学的指导作用。
[0005] 为了达到上述目的，本发明提供了一种设施葡萄低温灾害天气指数保险纯费率的 确定方法，包括以下步骤： (1) 对设施葡萄进行不同程度的低温处理； (2) 测定设施葡萄经步骤（1)不同程度低温处理下的各项生长指标及产量指标，计算设 施葡萄在不同程度低温处理下的产量损失率； (3) 利用历史气象数据，构建低温灾害天气发生的概率分布模型； (4) 测定设施葡萄相关的生长模型输入参数，结合步骤（2)得到不同程度低温处理下的 设施葡萄的产量损失率，建立产量损失率模型； (5) 根据步骤（4)得到的产量损失率模型和步骤（3)得到的低温灾害天气发生的概率 分布模型，计算保险纯费率。
[0006] 其中，步骤（5)中保险纯费率的计算公式为R=E[L0SS]= Σ LrXPi，其中R为保险 纯费率（％)，LOSS为作物损失，E[LOSS]为作物损失的期望值，Lr为不同低温灾害不同持续 天数下的产量损失率（％)，Pi为不同低温灾害天气的发生概率（％);Lr通过步骤（4)建立 的产量损失率模型计算，Pi通过步骤（3)构建的低温灾害天气发生的概率分布模型计算。
[0007] 步骤（3)中低温灾害天气发生的概率分布模型通过对各概率分布模型进行拟合， 通过拟合优度检验从中选择最优的概率分布模型；进行拟合的概率分布模型包括正态分 布、对数正态分布、指数分布、韦伯分布、伽马分布、Logistic分布和对数Logistic分布模 型；概率分布模型中各项参数的计算采用极大似然估计法；拟合优度检验采用AD检验、KS 检验或JB检验。
[0008] 步骤（4)中产量损失率模型的建立前，还需要进行设施小气候监测，利用测定的设 施小气候气象数据、所述设施葡萄相关的生长模型输入参数及不同程度低温处理下的设施 葡萄的产量损失率进行产量损失率模型的构建。步骤（4)中产量损失率模型通过DSSAT或 W0F0ST模型进行构建。
[0009] 当步骤（3)中历史气象数据选取为南京地区1957-2008年4月1日-4月30日期 间的日最低气温时，经本发明确定的保险纯费率为〇. 0327。
[0010] 本发明设施葡萄低温灾害天气指数保险纯费率的具体确定方法如下： 第一步：设施葡萄人工控制实验 将葡萄移栽至塑料盆中，选择盆栽葡萄10株，试验期间，设置低温5个处理：：日最低 温度设计〇°C、3°C、6°C、9°C、11°C，日最高温度设计为28°C;持续时间分别处理2、4、 6、8天；以日最高28°C/日最低18°C为对照。
[0011] 第二步：设施小气候观测方法 在选定的设施大棚内安装小气候观测站，收集室外气象数据和室内气象数据，设施环 境数据采集器为CR3000 (USA)，收集要素为1.5m高室内气温、总辐射、相对湿度20cm地 温、土壤水分、C02浓度等要素。气象要素采集频率为10s1次，30s储存1次取平均值。 如果当地没有监测的设施小气候气象资料，可以利用距离最近的气象台站资料，利用BP神 经网络模型进行将室外气象资料转换为设施小气候气象资料。
[0012] 第三步：生长发育指标测定及相关生长模型输入参数的测定 观测葡萄生育期和果实发育进程、测定叶片数、叶片面、枝蔓生长长度、测定果实生长 量及收获之后的产量，同时设施葡萄的作物品种、种植密度、行间距、种植深度、灌溉量、施 肥量、土壤种类、土壤有机碳、土壤中水的pH值、土壤深度等相关的作物模型输入参数。
[0013] 第四步：构建发生不同程度低温灾害的概率模型 利用气象数据对低温灾害不同程度发生的概率进行拟合，如利用不同地区的气象台的 资料（1950-2010年日最低温度的数据)，统计实际低温灾害发生的天数，并利用分布模型 (如正态、对数正态、βeta分布模型等）拟合，计算低温灾害持续天数的发生的概率。
[0014] 第五步：构建产量损失模型 根据人工控制实验获得的不同低温条件下的设施葡萄的产量损失率，结合设施小气候 气象数据及设施葡萄相关的作物模型输入参数，利用作物模型DSSAT与W0F0ST进行产量模 型的拟合，通过拟合优度的指标（RMSE、R2等）选择最优模型建立产量损失率模型，定量地确 定低温对设施葡萄产量的影响，然后利用在江苏省不同地区实时监测及历史的气象数据及 农业气象观测资料，输入到DSSAT或W0F0ST模型中进行减产率模型的调整，选择最符合当 地实际情况的产量损失率模型。
[0015] 第六步：计算保险纯费率 利用公式： R=E[L0SS]=ΣLrXPi (1) 计算保险纯费率，式中，为保险纯费率(％)，LOSS为作物损失，E[LOSS]为作物损失的期 望值即保险损失的期望值，Lr为不同灾害不同持续天数下的减产率（％)，Pi为不同灾害等 级的发生概率（％)。
[0016] 本发明相比现有具有以下优点：本发明综合采用人工控制实验方法与设施小气候 监测方法，确定了设施葡萄在低温灾害条件下的生长指标及产量指标，可以准确地得到低 温灾害条件下的天气指数保险纯费率，测定精度高，相对采用普通用历史产量数据进行估 算误差小，能够满足天气指数纯费率的计算要求。
【附图说明】
[0017] 图1为本发明设施葡萄低温灾害天气指数保险纯费率的确定方法流程图； 图2为本发明实施例1中BP神经网络模型结构； 图3为本发明实施例1设施小气候预报模型观测值与拟合值对比； 图4为本发明实施例1中历史气象数据所取四月日最低温度的经验累积函数分布图； 图5为本发明实施例1中历史气象数据所取四月日最低温度的正态分布概率图； 图6为本发明实施例1中历史气象数据所取四月日最低温度频数分布柱状图及正态分 布曲线拟合。
【具体实施方式】
[0018] 本发明以南京地区设施葡萄低温灾害天气指数保险纯费率的厘定为例。
[0019] 如图1所示，