论文部分内容阅读
摘要：利用2012年4月至2013年5月辽宁省盘锦市大洼县日光温室内外小气候观测资料，采用逐步回归方法建立辽宁地区不同天气类型和不同季节日光温室内最高气温及最低气温预报模型，并对拟合效果进行检验。结果表明，辽宁省大洼地区不同季节和不同天气类型日光温室内最高及最低气温预报模型差异较大，均通过了显著性检验；不同季节晴天和冬季阴天拟合效果稍差，其余季节不同天气类型拟合效果较好。春季、秋季、冬季中晴天、多云天、阴天、降水天日光温室内最低气温预报的平均绝对误差分别为0420～1130、0005～0064、0043～0150、0190～0270 ℃，均方根误差分别为0540～1530、0005～0076、0050～0180、0200～0320 ℃，平均相对误差分别为880%～1990%、007%～070%、040%～260%、260%～410%；晴天、多云天、阴天、降水天日光温室内最高气温预报的平均绝对误差分别为195～223、026～095、005～155、001～057 ℃，均方根误差分别为261～278、029～103、006～188、0012～0710 ℃，平均相对误差分别为680%～770%、022%～100%、030%～1180%、005%～340%。可见模型预报精度较高，可用于辽宁地区日光温室内最高和最低气温的预报。
关键词：日光温室；最高气温；最低气温；逐步回归；拟合检验；预报模型；天气类型；季节
中图分类号： S6255 1 文献标志码：
文章编号：1002-1302（2017）16-0276-06
收稿日期：2016-04-06
基金项目：中央级公益性科研院所基本科研业务费（编号：2014IAE-CMA04）；公益性行业（气象）科研专项（编号：GYHY201206024）；辽宁省“十二五”科学技术重大项目（编号：2011210002）。
作者简介：赵先丽（1977—），女，山东聊城人，硕士，副研究员，主要从事应用气象和生态气象研究。Tel：（024）83890246；E-mail：zhaoxianli2001@163com。
通信作者：李丽光，博士，副研究员，主要从事城市环境和全球气候变化研究。Tel：（024）83893253；E-mail：liliguangyjs@foxmailcom。
辽宁省属温带季风气候，夏季高温多雨，冬季寒冷干燥，是中国日光温室冬季生产不须加温的最北部地区，光热条件适合进行设施农业生产。辽宁省是中国日光温室蔬菜生产的发源地，日光温室蔬菜的面积和产量均居全国首位，设施农业水平为全国前列，以日光温室蔬菜生产为主的设施农业已成为辽宁省现代化农业的重要特征。截至2014年，中国设施农业总面积为411×106 hm2，产值超过8 000亿元，辽宁省设施农业总面积约为746×105 hm2，其中蔬菜播种面积为588×105 hm2，总产量达323×107 t，总产值达721亿元，居全国第一。目前，设施农业是辽宁省抗灾避灾、农业增效、农民增收的优势主导产业，已成为辽宁省现代农业的重要标志之一。
日光温室结构相对简单且面积较大，防御自然灾害及适应逆境能力较弱，主要通过吸收太阳能提高温室土壤温度和气温，当外界环境温度较低时，日光温室作物易遭受低温胁迫，造成减产或品质降低，而中午时段易遭受高温胁迫[6-7]。因此，根据日光温室外天气条件合理调控日光温室小气候是设施农业生产的关键，开展日光温室最低和最高气温预报以进行生产管理、调节气象条件十分重要，是防御设施农业低温冻害和高温热害的重要措施之一[8]。日光温室生产与日光温室内气象条件密切相关，尤其温度条件十分重要。国内目前对日光温室内气象要素的特征和预报研究较多，20世纪90年代开始对日光温室小气候（光照、温度、湿度、土壤温度）特征进行了系统研究并建立了统计模型。李军等基于能量和质量平衡原理，建立了以塑料大棚外气象要素为变量预报塑料大棚内温湿度的模型[9]。魏瑞江等分析了河北省日光温室内各气候因子之间的关系，并建立了日光温室内外的温度、湿度预报模型[10]。李德等选取常规气象预报因子和实测值，利用多元回归统计方法建立了宿州地区日光温室内最高及最低气温预报模型[11]。高丽娜等基于BP神经网络统计方法，以大棚外气象要素为自变量，建立了冬春季塑料大棚内极端温度和湿度的预报模型[12-14]。刘可群等通过分析大棚内小气候特征及其与大气候的关系，建立了武漢地区大棚内温度与太阳高度角及大气温度的相关模型[15]。张索铁等分别建立了黑龙江省、吉林省温室内不同季节不同天气型的最高和最低温度相关预报模型[8，16]。范辽生等对杭州地区冬季晴天、多云、寡照3种天气型单双层棚内的最低气温建立了预报模型，并进行了拟合和试报检验[17]。
由于日光温室结构、管理和气候地域性的差异，各地建立的日光温室气温预报模型不能进行推广应用；且以往建立的模型大部分为温室内气温与外界气温之间的单因子模型，未考虑前一天的温室内小气候要素[18]。目前关于辽宁省日光温室内气温预报的研究较少[19-20]，本研究利用2012年4月至2013年5月辽宁省大洼县日光温室内外的观测资料，根据日光温室内小气候条件与温室外气温、相对湿度等气象要素的关系，采用逐步回归方法，建立可动态预测温室内最高气温和最低气温的小气候模型，以期为辽宁省日光温室蔬菜生产提供气象预报服务，为农民增产和增收提供技术支撑。
1材料与方法
11研究区概况
辽宁省大洼地区属温带大陆性半湿润季风气候，四季分明、气候温和、雨量集中、光照充足。年平均气温为93 ℃，年平均最高气温为327 ℃，年平均最低气温为238 ℃，年平均降水量为 6473 mm，年平均无霜期为 208 d，年日照时数为 2 816 h，气候特点适合开展设施农业生产[20]。 12试验设计
2012年4月至2013年5月在辽宁省盘锦市大洼县王家镇华侨村日光温室（122°03′45″E、40°56′08″N）内进行试验，温室结构为单面坡面，坐北朝南，东西走向，长度为100 m，跨度为10 m，脊高为35 m。东、北及西面墙为双层砖墙，中间填充草帘、纸被、棉被及薄膜，北墙厚050 m，东西墙厚 024 m，温室覆盖008 mm厚的聚乙烯塑料薄膜，膜上覆盖草帘和保温被，日出后揭开，日落前后盖上。温室设有上通风口，采用半自动放风方式，基本管理方式为每日08：00后揭开草帘，16：00覆盖草帘[20]。
温室内无供暖设施，无照明设备，采用人工方式灌溉和施肥。温室内种植作物为番茄，品种为靓粉1号，主要观测项目为番茄出苗期、移栽期、开花期、结果期和采摘期，记录各发育期开始期、普遍期、末期出现日期；田间管理主要记录整地、移栽、灌水、施肥、喷药、铲趟、除草、采摘的时间和措施。番茄观测时间为2012年4月16日至5月31日、2012年10月1日至2013年5月31日。番茄行距为103 m，株距为024 m，2011年9月移栽，2012年4月21日进入采摘期，2012年5月29日采摘结束；2012年10月15日番茄移栽，2013年1月1日进入结果初期，2013年3月6日进入采摘期，2013年6月18日采摘结束[20]。
13资料来源
温室内外小气候监测采用锦州阳光气象科技有限公司生产的TRM-ZS3型设施农业小气候观测仪，观测精度分别为±01 ℃、±2%，每10 min采集1次数据，温室内传感器距离地面10 m，位于温室中间，温室外传感器距离地面15 m，观测项目为温室内气温、温室外气温、温室内相对湿度、温室外相对湿度、太阳总辐射、光合辐射。日平均值（气温和相对湿度）取1 d内144个数据的平均值，日最高值和最低值（气温和相对湿度）取1 d内144个数据中的最高和最低值，剔除缺测和异常的观测数据。日光温室外日照时数和日降水量数据来源于大洼县气象站（122°04′E、41°01′N），距日光溫室约105 km[20]。
14研究方法
141预报时段划分
将日光温室生产季划分为秋季（10至11月）、冬季（12月至翌年2月）、春季（3至5月）3个时段，建立不同季节日光温室内最高气温和最低气温预报模型。
142天气型划分
天气型划分：日照时数≥6 h为晴天，3 h<日照时数<6 h为多云天，0 h<日照时数≤3 h为阴天，在阴天情况下如有降水，即为降水天[21-22]。
143模型建立
本研究采用2012年4月17日至2013年5月31日大洼日光温室的观测数据进行建模，剔除异常值，有效数据为237 d，晴天为181 d，多云天为15 d，阴天为25 d，降水天为16 d；春季晴天为88 d，多云天为7 d，阴天为9 d，降水天为7 d；秋季晴天为32 d，多云天为2 d，阴天为4 d，降水天为7 d；冬季晴天为61 d，多云天为6 d，阴天为12 d，降水天为2 d。采用前一日降水量x1、前一日日照时数x2、前一日温室内（最高气温x3、最低气温x4、平均气温x5）和温室外气温（最高气温x6、最低气温x7、平均气温x8）、前一日温室内（最大相对湿度x9、最小相对湿度x10、平均相对湿度x11）和温室外相对湿度（最大相对湿度x12、最小相对湿度x13、平均相对湿度x14）及当日预报的降水量x15、当日预报的温室外气温（最高气温x16、最低气温x17、平均气温x18）、当日预报的温室外相对湿度（最大相对湿度x19、最小相对湿度x20、平均相对湿度x21）共21个因素作为预报因子，采用逐步回归法建立预报方程。
144拟合检验
采用均方根误差RMSE、相对误差RE和绝对误差AE对建立的模型进行评价分析。
3）
式中：Pi为预测值，Ai为实测值，n为样本数。RMSE、RE、AE越小，说明预报误差越小，预报越准确。
2结果与分析
21春季气温预报模型
211春季最低气温预报模型
由引进的预报因子可见，春季晴天时，决定日光温室内最低气温的主要影响因子为当日和前一日日光温室外最低气温及前一日日光温室内最低气温，前一日温室内最低气温、温室外最低气温、当日温室外最低气温与温室内最低气温相关系数均通过了0001水平显著性检验，一般外界气温越低，日光温室内气温也越低。李德等研究认为，春季晴天前1～2 d日光温室外温度条件对日光温室最低气温具有1～2 d的滞后效应[11]，本研究结果与之一致。而春季非晴天时，日光温室内最低气温的主要影响因子为前一日日光温室外最高气温与最低气温、相对湿度，前一日日光温室内最低气温与平均气温，当日日光温室外最低气温和平均气温等，其中阴天时日光温室内最低气温与前一日温室内最低气温、前一日温室内平均气温、当日温室外最低气温、当日温室外平均气温、当日温室外最大相对湿度显著相关（表1）。
由春季日光温室内最低气温反演可知，晴天温室内最低气温预报平均绝对误差为1130 ℃，均方根误差为1530 ℃，平均相对误差为1990%；多云天温室内最低气温预报平均绝对误差为0064 ℃，均方根误差为0076 ℃，平均相对误差为070%；阴天温室内最低气温预报平均绝对误差为 0043 ℃，均方根误差为0050 ℃，平均相对误差为040%；降水天温室内最低气温预报平均绝对误差为 0190 ℃，均方根误差为0200 ℃，平均相对误差为260%（表1、图1）。由
此可见，大洼春季晴天时日光温室内最低气温预报误差较大，其次为降水天，多云天和阴天预报效果较好。
212春季最高气温预报模型由表2可见，日光温室晴天最高气温引入的预报因子主要为反映空气水汽含量的当日和前一日降水量，决定日光温室基础气温的前一日温室内平均气温、当日温室外最高气温和当日温室外平均气温。而春季非晴天时，日光温室内最高气温的主要影响因子为前一日日光温室内最高气温与平均气温及当日温室外的湿度条件等。 晴天温室内最高气温预报平均绝对误差为2170 ℃，均方根误差为2780 ℃，平均相对误差为760%；多云天温室内最高气温预报平均绝对误差为0260 ℃，均方根误差为 0290 ℃，平均相对误差为100%；阴天温室内最高气温预报平均绝对误差为0910 ℃，均方根误差为0960 ℃，平均相对误差为410%；降水天温室内最高气温预报平均绝对误差为001 ℃，均方根误差为0012 ℃，平均相对误差为005%（图2、表2）。由此可见，大洼春季晴天时日光温室内最高气温预报误差稍大，其次为阴天，多云天和降水天预报效果较好。
22秋季气温预报模型
221秋季最低气温预报模型
由表3可知，秋季晴天温室内最低气温主要由前一日温室内平均气温、前一日温室外最大相对湿度等决定，而前一日温室外最低气温为负贡献。分析表明，前一日温室外最低气温偏高时，温室内最低气温增高，同时使温室内最高气温偏高，通风换气时间延长，一定程度上导致温室内最低气温降低。而秋季非晴天时，日光温室内最低气温的主要影响因子为前一日温室内平均气温、前一日[CM（25]温室外平均气温及前一日温室外平均相对湿度等。可见，
注：多云天样本仅为2 d，本研究未进行建模。表4同。
秋季非晴天时前一日外界气象条件对温室内最低气温影响显著。
由表3、图3可知，秋季日光温室内最低气温预报平均绝对误差为0920 ℃，均方根误差为1180 ℃，平均相对误差为1880%；阴天温室内最低气温预报平均绝对误差为 0070 ℃，均方根误差为 0080 ℃，平均相对误差为090%；降水天温室内最低气温预报平均绝对误差为0270 ℃，均方根误差为0320 ℃，平均相对误差为410%。由此可见，大洼秋季晴天温室内最低气温预报误差较大，其次为降水天，阴天预报效果较好。
222秋季最高气温预报模型由表4可知，晴天温室内最高气温主要影响因子为前一日温室内最高气温、前一日温室外平均相对湿度、当日温室外最高气温，而当日温室外平均气温为负贡献，主要是由于外界温度高，导致温室内最高气温偏高，通风时间长，降低了日光温室内的基础温度。秋季非晴天时，日光温室内最高气温主要受前一日温室内最低气温、前一日温室外最高气温、前一日温室外最大相对湿度及当日温室外最高气温影响。
秋季晴天日光温室内最高气温预报平均绝对误差为2230 ℃，均方根误差为2780 ℃，平均相对误差为770%；阴天温室内最高气温预报平均绝对误差为0050 ℃，均方根误差为0060 ℃，平均相对误差为030%；降水天温室内最高气温预报平均绝对误差为0570 ℃，均方根误差为 0710 ℃，平均相对误差为340%（表4、图4）。由此可见，大洼秋季晴天温室内最高气温预报误差稍大，其次为降水天，阴天预报效果较好。
23冬季气温预报模型
231冬季最低气温预报模型
由入选预报因子可见，冬季晴天时，决定日光温室内最低气温的主要气象要素为前一日温室内最低气温与平均气温、前一日温室外最低气温、当日降水量及当日温室外最高气温与平均气温，相关系数均通过了0001水平的显著性检验。而冬季非晴天时，日光温室内最低气温主要受前一日外界相对湿度、当日温室外气温等的影响（表5）。分析认为，前一日外界环境的相对湿度高，空气温度下降受到抑制，使夜间外界温度不会降得过低，从而影响温室内外的热交换速度和强度，最终影响次日温室内最低气温。
注：冬季降水天样本仅为2 d，本研究未建模。表6同。
由表5、图5可见，秋季晴天温室内最低气温预报平均绝对误差为0420 ℃，均方根误差为0540 ℃，平均相对误差为880%；多云天温室内最低气温预报平均绝对误差为 0005 ℃，均方根误差为0005 ℃，平均相对误差为007%；阴天温室内最低气温预报平均绝对误差为0150 ℃，均方根误差为0180 ℃，平均相对误差为260%。
[FK（W11][TPZXL5tif][FK）]
232冬季最高气温预报模型由表6可见，冬季晴天温室内最高气温受前一日温室外最低气温与平均气温、当日降水量及当日温室外平均气温的影响，由此进一步表明，日光温室外气温对日光温室内气温影响显著。而冬季非晴天时，决定日光温室内最高气温的为前一日温室外最小相对湿度、当日温室外最大相对湿度、当日温室外最低气温等，还引入了日照时数预报因子，即外界的日照对日光温室内的气温也有一定的影响，进一步表明非晴天时，增加太阳散射光可提高日光温室内气温。
由表6、图6可见，冬季晴天日光温室内最高气温预报平均绝对误差为1950 ℃，均方根误差为2610 ℃，平均相对误差为680%；多云天温室内最高气温预报平均绝对误差为 0950 ℃，均方根误差为1030 ℃，平均相对误差为022%；阴天温室内最高气温预报平均绝对误差为1550 ℃，均方根误差为1880 ℃，平均相对误差为1180%。由此可见，大洼冬季阴天温室内最高气温预报误差稍大，其次为晴天，多云天预报效果较好。
3结论与讨论
本研究采用逐步回归方法，建立适合大洼地区不同季节不同天气类型的温室内最低气温、最高气温的预报模型，模型均通过了显著性检验。大洼不同季节晴天温室内最低气温和最高气温预测模型的准确率稍差，春季最低、最高气温平均相对误差为1990%、760%，秋季最低、最高气温平均相对误差为1880%、770%，冬季最低、最高气温平均相对误差为880%、680%；阴天温室内最高气温预报准确率稍差，平均相对误差为118%；其余温室内最高气温和最低气温预报准确率较高，平均相对误差为005%～410%，均低于50%，预测准确率在950%以上。部分温室内气温预测模型预报精度高于张索铁等模型的預报精度[8，11，17]，可为大洼地区预防温室蔬菜冻害和热害提供参考，以提高蔬菜生产效益。 在全球气候变暖的背景下，极端气象事件发生频率增大，北方地区日光温室生产的高温、冻害、大风灾害和雪灾频繁发生，因此须进行温室内最低气温、最高气温的预测预报技术研究，为防灾减灾和温室生产管理提供依据。
本研究中大洼地区日光温室内最低气温和最高气温预测模型是基于日光温室内的观测资料及常规天气预报资料建立的，模型精度取决于常规天气预报水平，因此提高日光温室内最低气温和最高气温预报模型的精度，须进一步提高天气预报水平和积累试验资料进行深入研究。由于观测作物和观测资料有限，本研究未进行模型试预报，须进一步积累观测资料进行模型检验。依据当日天气现象分为4种天气类型，但在天气类型转化时日光温室最高气温和最低气温预测须考虑前一日的天气情况进行天气类型划分，须进一步深入研究。
参考文献：
陈艳秋，陈宇，吴曼丽 辽宁省设施农业暴雪及大风灾情特征分析[J] 安徽农业科学，2013，41（4）：1660-1662
陈妮娜，蒋大凯，王瀛，等 辽宁省设施农业大风和暴雪致灾指标[J] 江苏农业科学，2013，41（11）：386-387
[3]郎立新，史书强，张鹏，等 辽宁省设施农业发展分析[J] 园艺与种苗，2011 （1）：54-57
[4]孙立德 温室气象与作物保护研究[M] 沈阳：辽宁科学技术出版社，2012：1-151
[5]吕霞 辽宁省设施农业生产存在问题及对策[J] 农业科技与装备，2013 （1）：74-75
[6]阎琦，陈妮娜，田莉，等 辽宁设施农业致灾暴雪时空分布及天气学模型[J] 江苏农业科学，2016，44 （1）：373-376
[7]马成芝，孙立德，穆春华 喀左日光温室内气温变化规律及其应用[J] 气象与环境学报，2007，23（5）：49-52
[8]张索铁，马树庆，庞义，等 东北地区大型日光温室最高、最低温度预报模式探讨[J] 中国农学通报，2014，30（23）：249-253
[9]李军，姚益平，徐蕊，等 长江下游防虫网覆盖塑料大棚内温湿度模拟[J] 农业工程学报，2010，26（6）：238-244
[10]魏瑞江，王春乙，范增禄石家庄地区日光温室冬季小气候特征及其与大气候的关系[J] 气象，2010，36（1）：97-103
[11]李德，张学贤，祁宦，等 宿州日光温室内部最高和最低气温的预报模型[J] 中国农业气象，2013，34（2）：170-178
[12]高丽娜，孙擎，郭翠荣，等 山西日光温室逐日极端气温预测模型研究[J] 中国农学通报，2015，31（15）：240-246
[13]刘淑梅，薛庆禹，黎贞发，等 基于BP神经网络的日光温室气温预报模型[J] 中国农业大学学报，2015，20（1）：176-184
[14]金志凤，符国槐，黄海静，等 基于BP神经网络的杨梅大棚内气[HJ]
温预测模型研究[J] 中国农业气象，2011，32（3）：362-367
[15]刘可群，魏明锋，杨文刚 大棚小气候特征及其与大气候的关系[J] 气象，2008，34（7）：101-107
[16]王萍，刘春雪，王秋京，等 黑龙江省温室小气候变化特征及预报模型的初步研究[J] 黑龙江农业科学，2014（5）：75-79
[17]范辽生，朱兰娟，柴伟国，等 杭州冬季塑料大棚内气温变化特征及日最低气温预报模型[J] 中国农业气象，2014，35（3）：268-275
[18]袁静，李树军，崔建云，等 山东寿光冬季日光温室内温度变化特征及低温预报[J] 中国农学通报，2012，28（3）：300-304
[19]张淑杰，杨再强，陈艳秋，等 低温、弱光、高湿胁迫对日光温室番茄花期生理生化指标的影响[J] 生态学杂志，2014，33（11）：2995-3001
[20]赵先丽，张淑杰，陈艳秋，等 辽宁大洼地区日光温室小气候特征[J] 生态学杂志，2014，33（10）：2656-2663
[21]魏瑞江 日光温室低温寡照灾害指标[J] 气象科技，2003，31（1）：50-53
[22]魏瑞江，李春强，康西言河北省日光温室低溫寡照灾害风险分析[J] 自然灾害学报，2008，17（3）：56-62
关键词：日光温室；最高气温；最低气温；逐步回归；拟合检验；预报模型；天气类型；季节
中图分类号： S6255 1 文献标志码：
文章编号：1002-1302（2017）16-0276-06
收稿日期：2016-04-06
基金项目：中央级公益性科研院所基本科研业务费（编号：2014IAE-CMA04）；公益性行业（气象）科研专项（编号：GYHY201206024）；辽宁省“十二五”科学技术重大项目（编号：2011210002）。
作者简介：赵先丽（1977—），女，山东聊城人，硕士，副研究员，主要从事应用气象和生态气象研究。Tel：（024）83890246；E-mail：zhaoxianli2001@163com。
通信作者：李丽光，博士，副研究员，主要从事城市环境和全球气候变化研究。Tel：（024）83893253；E-mail：liliguangyjs@foxmailcom。
辽宁省属温带季风气候，夏季高温多雨，冬季寒冷干燥，是中国日光温室冬季生产不须加温的最北部地区，光热条件适合进行设施农业生产。辽宁省是中国日光温室蔬菜生产的发源地，日光温室蔬菜的面积和产量均居全国首位，设施农业水平为全国前列，以日光温室蔬菜生产为主的设施农业已成为辽宁省现代化农业的重要特征。截至2014年，中国设施农业总面积为411×106 hm2，产值超过8 000亿元，辽宁省设施农业总面积约为746×105 hm2，其中蔬菜播种面积为588×105 hm2，总产量达323×107 t，总产值达721亿元，居全国第一。目前，设施农业是辽宁省抗灾避灾、农业增效、农民增收的优势主导产业，已成为辽宁省现代农业的重要标志之一。
日光温室结构相对简单且面积较大，防御自然灾害及适应逆境能力较弱，主要通过吸收太阳能提高温室土壤温度和气温，当外界环境温度较低时，日光温室作物易遭受低温胁迫，造成减产或品质降低，而中午时段易遭受高温胁迫[6-7]。因此，根据日光温室外天气条件合理调控日光温室小气候是设施农业生产的关键，开展日光温室最低和最高气温预报以进行生产管理、调节气象条件十分重要，是防御设施农业低温冻害和高温热害的重要措施之一[8]。日光温室生产与日光温室内气象条件密切相关，尤其温度条件十分重要。国内目前对日光温室内气象要素的特征和预报研究较多，20世纪90年代开始对日光温室小气候（光照、温度、湿度、土壤温度）特征进行了系统研究并建立了统计模型。李军等基于能量和质量平衡原理，建立了以塑料大棚外气象要素为变量预报塑料大棚内温湿度的模型[9]。魏瑞江等分析了河北省日光温室内各气候因子之间的关系，并建立了日光温室内外的温度、湿度预报模型[10]。李德等选取常规气象预报因子和实测值，利用多元回归统计方法建立了宿州地区日光温室内最高及最低气温预报模型[11]。高丽娜等基于BP神经网络统计方法，以大棚外气象要素为自变量，建立了冬春季塑料大棚内极端温度和湿度的预报模型[12-14]。刘可群等通过分析大棚内小气候特征及其与大气候的关系，建立了武漢地区大棚内温度与太阳高度角及大气温度的相关模型[15]。张索铁等分别建立了黑龙江省、吉林省温室内不同季节不同天气型的最高和最低温度相关预报模型[8，16]。范辽生等对杭州地区冬季晴天、多云、寡照3种天气型单双层棚内的最低气温建立了预报模型，并进行了拟合和试报检验[17]。
由于日光温室结构、管理和气候地域性的差异，各地建立的日光温室气温预报模型不能进行推广应用；且以往建立的模型大部分为温室内气温与外界气温之间的单因子模型，未考虑前一天的温室内小气候要素[18]。目前关于辽宁省日光温室内气温预报的研究较少[19-20]，本研究利用2012年4月至2013年5月辽宁省大洼县日光温室内外的观测资料，根据日光温室内小气候条件与温室外气温、相对湿度等气象要素的关系，采用逐步回归方法，建立可动态预测温室内最高气温和最低气温的小气候模型，以期为辽宁省日光温室蔬菜生产提供气象预报服务，为农民增产和增收提供技术支撑。
1材料与方法
11研究区概况
辽宁省大洼地区属温带大陆性半湿润季风气候，四季分明、气候温和、雨量集中、光照充足。年平均气温为93 ℃，年平均最高气温为327 ℃，年平均最低气温为238 ℃，年平均降水量为 6473 mm，年平均无霜期为 208 d，年日照时数为 2 816 h，气候特点适合开展设施农业生产[20]。 12试验设计
2012年4月至2013年5月在辽宁省盘锦市大洼县王家镇华侨村日光温室（122°03′45″E、40°56′08″N）内进行试验，温室结构为单面坡面，坐北朝南，东西走向，长度为100 m，跨度为10 m，脊高为35 m。东、北及西面墙为双层砖墙，中间填充草帘、纸被、棉被及薄膜，北墙厚050 m，东西墙厚 024 m，温室覆盖008 mm厚的聚乙烯塑料薄膜，膜上覆盖草帘和保温被，日出后揭开，日落前后盖上。温室设有上通风口，采用半自动放风方式，基本管理方式为每日08：00后揭开草帘，16：00覆盖草帘[20]。
温室内无供暖设施，无照明设备，采用人工方式灌溉和施肥。温室内种植作物为番茄，品种为靓粉1号，主要观测项目为番茄出苗期、移栽期、开花期、结果期和采摘期，记录各发育期开始期、普遍期、末期出现日期；田间管理主要记录整地、移栽、灌水、施肥、喷药、铲趟、除草、采摘的时间和措施。番茄观测时间为2012年4月16日至5月31日、2012年10月1日至2013年5月31日。番茄行距为103 m，株距为024 m，2011年9月移栽，2012年4月21日进入采摘期，2012年5月29日采摘结束；2012年10月15日番茄移栽，2013年1月1日进入结果初期，2013年3月6日进入采摘期，2013年6月18日采摘结束[20]。
13资料来源
温室内外小气候监测采用锦州阳光气象科技有限公司生产的TRM-ZS3型设施农业小气候观测仪，观测精度分别为±01 ℃、±2%，每10 min采集1次数据，温室内传感器距离地面10 m，位于温室中间，温室外传感器距离地面15 m，观测项目为温室内气温、温室外气温、温室内相对湿度、温室外相对湿度、太阳总辐射、光合辐射。日平均值（气温和相对湿度）取1 d内144个数据的平均值，日最高值和最低值（气温和相对湿度）取1 d内144个数据中的最高和最低值，剔除缺测和异常的观测数据。日光温室外日照时数和日降水量数据来源于大洼县气象站（122°04′E、41°01′N），距日光溫室约105 km[20]。
14研究方法
141预报时段划分
将日光温室生产季划分为秋季（10至11月）、冬季（12月至翌年2月）、春季（3至5月）3个时段，建立不同季节日光温室内最高气温和最低气温预报模型。
142天气型划分
天气型划分：日照时数≥6 h为晴天，3 h<日照时数<6 h为多云天，0 h<日照时数≤3 h为阴天，在阴天情况下如有降水，即为降水天[21-22]。
143模型建立
本研究采用2012年4月17日至2013年5月31日大洼日光温室的观测数据进行建模，剔除异常值，有效数据为237 d，晴天为181 d，多云天为15 d，阴天为25 d，降水天为16 d；春季晴天为88 d，多云天为7 d，阴天为9 d，降水天为7 d；秋季晴天为32 d，多云天为2 d，阴天为4 d，降水天为7 d；冬季晴天为61 d，多云天为6 d，阴天为12 d，降水天为2 d。采用前一日降水量x1、前一日日照时数x2、前一日温室内（最高气温x3、最低气温x4、平均气温x5）和温室外气温（最高气温x6、最低气温x7、平均气温x8）、前一日温室内（最大相对湿度x9、最小相对湿度x10、平均相对湿度x11）和温室外相对湿度（最大相对湿度x12、最小相对湿度x13、平均相对湿度x14）及当日预报的降水量x15、当日预报的温室外气温（最高气温x16、最低气温x17、平均气温x18）、当日预报的温室外相对湿度（最大相对湿度x19、最小相对湿度x20、平均相对湿度x21）共21个因素作为预报因子，采用逐步回归法建立预报方程。
144拟合检验
采用均方根误差RMSE、相对误差RE和绝对误差AE对建立的模型进行评价分析。
3）
式中：Pi为预测值，Ai为实测值，n为样本数。RMSE、RE、AE越小，说明预报误差越小，预报越准确。
2结果与分析
21春季气温预报模型
211春季最低气温预报模型
由引进的预报因子可见，春季晴天时，决定日光温室内最低气温的主要影响因子为当日和前一日日光温室外最低气温及前一日日光温室内最低气温，前一日温室内最低气温、温室外最低气温、当日温室外最低气温与温室内最低气温相关系数均通过了0001水平显著性检验，一般外界气温越低，日光温室内气温也越低。李德等研究认为，春季晴天前1～2 d日光温室外温度条件对日光温室最低气温具有1～2 d的滞后效应[11]，本研究结果与之一致。而春季非晴天时，日光温室内最低气温的主要影响因子为前一日日光温室外最高气温与最低气温、相对湿度，前一日日光温室内最低气温与平均气温，当日日光温室外最低气温和平均气温等，其中阴天时日光温室内最低气温与前一日温室内最低气温、前一日温室内平均气温、当日温室外最低气温、当日温室外平均气温、当日温室外最大相对湿度显著相关（表1）。
由春季日光温室内最低气温反演可知，晴天温室内最低气温预报平均绝对误差为1130 ℃，均方根误差为1530 ℃，平均相对误差为1990%；多云天温室内最低气温预报平均绝对误差为0064 ℃，均方根误差为0076 ℃，平均相对误差为070%；阴天温室内最低气温预报平均绝对误差为 0043 ℃，均方根误差为0050 ℃，平均相对误差为040%；降水天温室内最低气温预报平均绝对误差为 0190 ℃，均方根误差为0200 ℃，平均相对误差为260%（表1、图1）。由
此可见，大洼春季晴天时日光温室内最低气温预报误差较大，其次为降水天，多云天和阴天预报效果较好。
212春季最高气温预报模型由表2可见，日光温室晴天最高气温引入的预报因子主要为反映空气水汽含量的当日和前一日降水量，决定日光温室基础气温的前一日温室内平均气温、当日温室外最高气温和当日温室外平均气温。而春季非晴天时，日光温室内最高气温的主要影响因子为前一日日光温室内最高气温与平均气温及当日温室外的湿度条件等。 晴天温室内最高气温预报平均绝对误差为2170 ℃，均方根误差为2780 ℃，平均相对误差为760%；多云天温室内最高气温预报平均绝对误差为0260 ℃，均方根误差为 0290 ℃，平均相对误差为100%；阴天温室内最高气温预报平均绝对误差为0910 ℃，均方根误差为0960 ℃，平均相对误差为410%；降水天温室内最高气温预报平均绝对误差为001 ℃，均方根误差为0012 ℃，平均相对误差为005%（图2、表2）。由此可见，大洼春季晴天时日光温室内最高气温预报误差稍大，其次为阴天，多云天和降水天预报效果较好。
22秋季气温预报模型
221秋季最低气温预报模型
由表3可知，秋季晴天温室内最低气温主要由前一日温室内平均气温、前一日温室外最大相对湿度等决定，而前一日温室外最低气温为负贡献。分析表明，前一日温室外最低气温偏高时，温室内最低气温增高，同时使温室内最高气温偏高，通风换气时间延长，一定程度上导致温室内最低气温降低。而秋季非晴天时，日光温室内最低气温的主要影响因子为前一日温室内平均气温、前一日[CM（25]温室外平均气温及前一日温室外平均相对湿度等。可见，
注：多云天样本仅为2 d，本研究未进行建模。表4同。
秋季非晴天时前一日外界气象条件对温室内最低气温影响显著。
由表3、图3可知，秋季日光温室内最低气温预报平均绝对误差为0920 ℃，均方根误差为1180 ℃，平均相对误差为1880%；阴天温室内最低气温预报平均绝对误差为 0070 ℃，均方根误差为 0080 ℃，平均相对误差为090%；降水天温室内最低气温预报平均绝对误差为0270 ℃，均方根误差为0320 ℃，平均相对误差为410%。由此可见，大洼秋季晴天温室内最低气温预报误差较大，其次为降水天，阴天预报效果较好。
222秋季最高气温预报模型由表4可知，晴天温室内最高气温主要影响因子为前一日温室内最高气温、前一日温室外平均相对湿度、当日温室外最高气温，而当日温室外平均气温为负贡献，主要是由于外界温度高，导致温室内最高气温偏高，通风时间长，降低了日光温室内的基础温度。秋季非晴天时，日光温室内最高气温主要受前一日温室内最低气温、前一日温室外最高气温、前一日温室外最大相对湿度及当日温室外最高气温影响。
秋季晴天日光温室内最高气温预报平均绝对误差为2230 ℃，均方根误差为2780 ℃，平均相对误差为770%；阴天温室内最高气温预报平均绝对误差为0050 ℃，均方根误差为0060 ℃，平均相对误差为030%；降水天温室内最高气温预报平均绝对误差为0570 ℃，均方根误差为 0710 ℃，平均相对误差为340%（表4、图4）。由此可见，大洼秋季晴天温室内最高气温预报误差稍大，其次为降水天，阴天预报效果较好。
23冬季气温预报模型
231冬季最低气温预报模型
由入选预报因子可见，冬季晴天时，决定日光温室内最低气温的主要气象要素为前一日温室内最低气温与平均气温、前一日温室外最低气温、当日降水量及当日温室外最高气温与平均气温，相关系数均通过了0001水平的显著性检验。而冬季非晴天时，日光温室内最低气温主要受前一日外界相对湿度、当日温室外气温等的影响（表5）。分析认为，前一日外界环境的相对湿度高，空气温度下降受到抑制，使夜间外界温度不会降得过低，从而影响温室内外的热交换速度和强度，最终影响次日温室内最低气温。
注：冬季降水天样本仅为2 d，本研究未建模。表6同。
由表5、图5可见，秋季晴天温室内最低气温预报平均绝对误差为0420 ℃，均方根误差为0540 ℃，平均相对误差为880%；多云天温室内最低气温预报平均绝对误差为 0005 ℃，均方根误差为0005 ℃，平均相对误差为007%；阴天温室内最低气温预报平均绝对误差为0150 ℃，均方根误差为0180 ℃，平均相对误差为260%。
[FK（W11][TPZXL5tif][FK）]
232冬季最高气温预报模型由表6可见，冬季晴天温室内最高气温受前一日温室外最低气温与平均气温、当日降水量及当日温室外平均气温的影响，由此进一步表明，日光温室外气温对日光温室内气温影响显著。而冬季非晴天时，决定日光温室内最高气温的为前一日温室外最小相对湿度、当日温室外最大相对湿度、当日温室外最低气温等，还引入了日照时数预报因子，即外界的日照对日光温室内的气温也有一定的影响，进一步表明非晴天时，增加太阳散射光可提高日光温室内气温。
由表6、图6可见，冬季晴天日光温室内最高气温预报平均绝对误差为1950 ℃，均方根误差为2610 ℃，平均相对误差为680%；多云天温室内最高气温预报平均绝对误差为 0950 ℃，均方根误差为1030 ℃，平均相对误差为022%；阴天温室内最高气温预报平均绝对误差为1550 ℃，均方根误差为1880 ℃，平均相对误差为1180%。由此可见，大洼冬季阴天温室内最高气温预报误差稍大，其次为晴天，多云天预报效果较好。
3结论与讨论
本研究采用逐步回归方法，建立适合大洼地区不同季节不同天气类型的温室内最低气温、最高气温的预报模型，模型均通过了显著性检验。大洼不同季节晴天温室内最低气温和最高气温预测模型的准确率稍差，春季最低、最高气温平均相对误差为1990%、760%，秋季最低、最高气温平均相对误差为1880%、770%，冬季最低、最高气温平均相对误差为880%、680%；阴天温室内最高气温预报准确率稍差，平均相对误差为118%；其余温室内最高气温和最低气温预报准确率较高，平均相对误差为005%～410%，均低于50%，预测准确率在950%以上。部分温室内气温预测模型预报精度高于张索铁等模型的預报精度[8，11，17]，可为大洼地区预防温室蔬菜冻害和热害提供参考，以提高蔬菜生产效益。 在全球气候变暖的背景下，极端气象事件发生频率增大，北方地区日光温室生产的高温、冻害、大风灾害和雪灾频繁发生，因此须进行温室内最低气温、最高气温的预测预报技术研究，为防灾减灾和温室生产管理提供依据。
本研究中大洼地区日光温室内最低气温和最高气温预测模型是基于日光温室内的观测资料及常规天气预报资料建立的，模型精度取决于常规天气预报水平，因此提高日光温室内最低气温和最高气温预报模型的精度，须进一步提高天气预报水平和积累试验资料进行深入研究。由于观测作物和观测资料有限，本研究未进行模型试预报，须进一步积累观测资料进行模型检验。依据当日天气现象分为4种天气类型，但在天气类型转化时日光温室最高气温和最低气温预测须考虑前一日的天气情况进行天气类型划分，须进一步深入研究。
参考文献：
陈艳秋，陈宇，吴曼丽 辽宁省设施农业暴雪及大风灾情特征分析[J] 安徽农业科学，2013，41（4）：1660-1662
陈妮娜，蒋大凯，王瀛，等 辽宁省设施农业大风和暴雪致灾指标[J] 江苏农业科学，2013，41（11）：386-387
[3]郎立新，史书强，张鹏，等 辽宁省设施农业发展分析[J] 园艺与种苗，2011 （1）：54-57
[4]孙立德 温室气象与作物保护研究[M] 沈阳：辽宁科学技术出版社，2012：1-151
[5]吕霞 辽宁省设施农业生产存在问题及对策[J] 农业科技与装备，2013 （1）：74-75
[6]阎琦，陈妮娜，田莉，等 辽宁设施农业致灾暴雪时空分布及天气学模型[J] 江苏农业科学，2016，44 （1）：373-376
[7]马成芝，孙立德，穆春华 喀左日光温室内气温变化规律及其应用[J] 气象与环境学报，2007，23（5）：49-52
[8]张索铁，马树庆，庞义，等 东北地区大型日光温室最高、最低温度预报模式探讨[J] 中国农学通报，2014，30（23）：249-253
[9]李军，姚益平，徐蕊，等 长江下游防虫网覆盖塑料大棚内温湿度模拟[J] 农业工程学报，2010，26（6）：238-244
[10]魏瑞江，王春乙，范增禄石家庄地区日光温室冬季小气候特征及其与大气候的关系[J] 气象，2010，36（1）：97-103
[11]李德，张学贤，祁宦，等 宿州日光温室内部最高和最低气温的预报模型[J] 中国农业气象，2013，34（2）：170-178
[12]高丽娜，孙擎，郭翠荣，等 山西日光温室逐日极端气温预测模型研究[J] 中国农学通报，2015，31（15）：240-246
[13]刘淑梅，薛庆禹，黎贞发，等 基于BP神经网络的日光温室气温预报模型[J] 中国农业大学学报，2015，20（1）：176-184
[14]金志凤，符国槐，黄海静，等 基于BP神经网络的杨梅大棚内气[HJ]
温预测模型研究[J] 中国农业气象，2011，32（3）：362-367
[15]刘可群，魏明锋，杨文刚 大棚小气候特征及其与大气候的关系[J] 气象，2008，34（7）：101-107
[16]王萍，刘春雪，王秋京，等 黑龙江省温室小气候变化特征及预报模型的初步研究[J] 黑龙江农业科学，2014（5）：75-79
[17]范辽生，朱兰娟，柴伟国，等 杭州冬季塑料大棚内气温变化特征及日最低气温预报模型[J] 中国农业气象，2014，35（3）：268-275
[18]袁静，李树军，崔建云，等 山东寿光冬季日光温室内温度变化特征及低温预报[J] 中国农学通报，2012，28（3）：300-304
[19]张淑杰，杨再强，陈艳秋，等 低温、弱光、高湿胁迫对日光温室番茄花期生理生化指标的影响[J] 生态学杂志，2014，33（11）：2995-3001
[20]赵先丽，张淑杰，陈艳秋，等 辽宁大洼地区日光温室小气候特征[J] 生态学杂志，2014，33（10）：2656-2663
[21]魏瑞江 日光温室低温寡照灾害指标[J] 气象科技，2003，31（1）：50-53
[22]魏瑞江，李春强，康西言河北省日光温室低溫寡照灾害风险分析[J] 自然灾害学报，2008，17（3）：56-62