CN111460377B - 一种晴天条件下温室作物冠层饱和水汽压差计算方法 - Google Patents
一种晴天条件下温室作物冠层饱和水汽压差计算方法 Download PDFInfo
- Publication number
- CN111460377B CN111460377B CN202010183568.7A CN202010183568A CN111460377B CN 111460377 B CN111460377 B CN 111460377B CN 202010183568 A CN202010183568 A CN 202010183568A CN 111460377 B CN111460377 B CN 111460377B
- Authority
- CN
- China
- Prior art keywords
- greenhouse
- canopy
- vapor pressure
- water vapor
- saturated water
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Algebra (AREA)
- Evolutionary Biology (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Greenhouses (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
本发明属于冠层饱和水汽压差计算领域,具体涉及一种温室作物冠层饱和水汽压差的计算方法。该方法通过测量温室内作物冠层高(2.0m)处的空气温度、湿度、总辐射,计算得到温室内冠层水汽压亏缺(VPD),采用分段分析的方法对VPD进行了数值模拟,得到0:00~12:00和12:30~23:30两个时间段的饱和水汽压亏缺与室内冠层太阳辐射的拟合函数,为温室作物冠层水汽压差的计算提供了新的方法,对日光温室小气候变化特征及小气候模拟的研究具有指导意义。
Description
技术领域
本发明属于冠层饱和水汽压差技术领域,具体涉及一种温室作物冠层饱和水汽压差的计算方法。
背景技术
水汽压(Vapour Pressure)是大气压力中水汽的分压力,空气中的水分含量与水汽压直接相关。在一定温度下空气中水汽达到饱和时的分压力,称为饱和水汽压e0(T)。气温越高,空气储存水分的能力越大,饱和水汽压也越大。实际水汽压ea是实际大气压力中水汽的分压力。当空气没有达到饱和时,实际水汽压ea比饱和水汽压e0(T)小。饱和水汽压与实际水汽压之间的差称为水汽压亏缺量或饱和水汽压差(Vapour Pressure Deficit,即VPD),是表征空气实际蒸发能力的准确指示器。目前VPD的计算方法是根据定义计算的,需要对温室内冠层最高、最低气温、温度及干湿等数据进行测量计算后才能得到。该方法的其主要问题是需要测量的环境参数相对较多,特别是实际饱和水汽压的计算方法会随着环境参数的可靠性有不同的选择,如实际水汽压的计算可以通过露点温度计算,也可以根据需求通过干湿表资料计算,还可以通过相对湿度进行计算,不同的计算方法对测量时段、参数种类、可靠度都有不同的要求,这使得计算VPD的工作量有所增加,且温室内长期处于高温高湿环境下,需要的参数越多就越容易产生误差。
发明内容
针对目前冠层饱和水汽压差计算方法存在参数多、测量复杂且易产生误差的缺陷和问题,本发明提供一种晴天条件下温室作物冠层饱和水汽压差的计算方法。
本发明解决其技术问题所采用的方案是:一种晴天条件下温室作物冠层饱和水汽压差的计算方法,包括以下步骤:
步骤一、采用植物生长监测仪采集0:00~12:00和12:30~23:30两个时段中不同时间点的日光温室大棚内的温室环境参数和温室内作物的冠层太阳辐射并实时监测作物的生长状况,并计算两个时段中不同时间点的冠层饱和水汽压差;
步骤二、分别对0:00~12:00以及12:30~23:30两时段的室内冠层太阳辐射与饱和水汽压亏缺的关系进行模拟得到两时段的温室内饱和水汽压亏缺与室内冠层太阳辐射的回归模型:
0:00~12:00时段,VPD=a TIRi 2+b TIRi+c
12:30~23:30时段,VPD=a’TIRi 2+b’TIRi+c‘
式中:VPD为室内冠层水汽压亏缺,TIRi为室内冠层太阳辐射,a、b、c、a’、b’、c‘为常数;
步骤三、将步骤一得到的不同时间点的室内冠层太阳辐射分别代入步骤二中对应时段的回归模型中计算得到a、b、c、a’、b’、c‘常数数值,并得到两个时段具体的回归模型;
步骤四、将测得的不同时段的室内冠层太阳辐射数值代入步骤三对应时段的回归模型中即可得到温室内的饱和水汽压亏缺。
上述的一种晴天条件下温室作物冠层饱和水汽压差的计算方法,温室内温室环境参数和冠层太阳辐射的采集时间为每日的7:00、10:00、12:00、15:00和18:00。
上述的一种晴天条件下温室作物冠层饱和水汽压差的计算方法,所述温室环境参数包括温室内作物冠层高处的空气温度和湿度。
上述的一种晴天条件下温室作物冠层饱和水汽压差的计算方法,所述温室作物为日光温室番茄,0:00~12:00和12:30~23:30两个时间段的二元回归模型分别为:
0:00~12:00时段,
VPD=0.00002TIRi 2-0.00411TIRi-0.0056 R2=0.9901
12:30~23:30时段,
VPD=-0.00004TIRi 2+0.02TIRi+0.2012 R2=0.9583
式中:VPD为室内冠层水汽压亏缺,TIRi为室内冠层太阳辐射。
本发明的有益效果:本发明通过试验证明,在晴天VPD变化与总辐射直接有较强的拟合关系,采用分段分析的方法建立起晴天条件下0:00~12:00以及12:30~23:30两个时段的冠层饱和水汽压差与冠层太阳辐射的二元回归模型,从而可通过冠层太阳辐射实现对VPD的准确计算。该方法只涉及到总辐射一个环境参数,不涉及较多的环境参数,且冠层太阳辐射的测量成本较低,测量方法简单科学,测量结果误差较小,为温室作物冠层水气压差的计算提供了新的方法,对日光温室小气候变化特征及小气候模拟的研究有重要的指导意义。
附图说明
图1为晴天条件下室内冠层水汽压亏缺与太阳辐射对比图。
图2为晴天条件下室内冠层水汽压亏缺与太阳辐射关系图。
图3为0:00~12:00时段室内冠层水汽压亏缺与太阳辐射关系图。
图4为12:30~23:30时段室内冠层水汽压亏缺与太阳辐射关系图。
具体实施方式
水汽压亏缺VPD(Vapor Pressure Deficit)反映了空气中的实际水汽距离饱和的程度,值愈大,说明空气愈干燥;值愈小,空气愈接近于饱和。无论是研究植株的蒸腾还是水面蒸发,水汽压亏缺都是一个很重要的指标。但是目前水汽压亏缺的计算方法涉及较多参数,计算方法复杂的问题,本发明提供一种温室作物冠层饱和水汽压差的计算方法,下面结合附图和实施例对本发明作进一步阐述。
实施例1:本实施例提供一种晴天条件下温室作物冠层饱和水汽压差的计算方法,该方法包括以下步骤:
步骤一、采用植物生长监测仪采集0:00~12:00和12:30~23:30两个时段不同时间点的日光温室大棚内的温室环境参数和冠层太阳辐射并实时监测作物的生长状况,采集时间为每日的7:00、10:00、12:00、15:00和18:00,并计算两个时段不同时间点的冠层饱和水汽压差;
步骤二、分别对0:00~12:00以及12:30~23:30两时段的室内冠层太阳辐射与饱和水汽压亏缺的关系进行模拟得到两时段的温室内饱和水汽压亏缺与室内冠层太阳辐射的回归模型:
0:00~12:00时段,VPD=a TIRi 2+b TIRi+c
12:30~23:30时段,VPD=a’TIRi 2+b’TIRi+c‘
式中:VPD为室内冠层水汽压亏缺,TIRi为室内冠层太阳辐射,a、b、c、a’、b’、c‘为常数;
步骤三、将步骤一得到的两个时段不同时间点的室内冠层太阳辐射分别代入步骤二中对应时段的回归模型中计算得到a、b、c、a’、b’、c‘常数数值,并得到两个时段具体的回归模型;
步骤四、将测得的不同时段的室内冠层太阳辐射数值代入步骤三对应时段的回归模型中从而完成温室内的饱和水汽压亏缺的计算。
实施例2:本实施例提供一种晴天条件下温室作物冠层饱和水汽压差的计算方法,为了便于理解本发明,本实施例以2019年3月~4月在湖北省鄂州节水示范基地的温室内进行。该区地理位置东经114.52°,北纬30.23°,为季风气候区,冬季盛行偏北风,夏季盛行偏南风,属亚热带气候,无霜期约236天,年平均气温16.3℃;年降雨量831.8mm,年内分布不均,夏、秋降雨少。
试验温室属于高效节能性日光温室,剖面呈扇形,占地面积584㎡,北墙高2.5m(矢高3.2米)。温室采用钢结构支撑,薄膜为EVA 3层复合膜,顶部和侧面可人工开启用于通气,室外备有防寒被。试验期内温室内种植作物为番茄,灌溉方式为膜下滴灌,选取果实采摘期进行分析。
首先,在温室内采用LPS-05型植物生长监测仪自动采集温室环境参数并实时监测作物生长状况,自动采集数据时间步长为1h。采集项目有:温室内作物冠层高(2.0m)处的空气温度、湿度、总辐射;另外在温室内悬挂最高、最低气温计,温度计及干湿表以作对比测量,每天观测时间为7:00、10:00、12:00、15:00、18:00。
饱和水汽压与实际水汽压之间的差称为水汽压亏缺量或饱和水汽压差(VPD),VPD=e0(T)-ea,其中e0(T)为饱和水汽压,ea为实际水汽压。
e0(T):气温为T时的饱和水汽压,单位为千帕(kPa);T:空气温度,单位为摄氏度(℃)。
经监测发现在阴雨天,室内湿度接近饱和,VPD变化幅度不大,但外界太阳辐射受云层影响较大变化剧烈,两者间的相关性受到干扰。然而,在晴天,VPD与太阳辐射之间存在着很好的关联。如图1所示,12:00以前太阳辐射与VPD一直以几乎平行的坡度直线上升,太阳辐射在12:00左右达到最大值;12:00以后两者又以几乎平行的坡度直线下降,VPD在14:00左右达到最大值,室内冠层水汽压亏缺与太阳辐射关系如图2所示。
采用分段分析的方法对VPD进行了数值模拟,以12:00为分界点,分为0:00~12:00和12:30~23:30两个时间段分别对室内外辐射与饱和水汽压亏缺关系进行模拟,具体如下,两个时段的室内冠层水汽压亏缺与太阳辐射关系分别参见图3和图4。
在0:00~12:00时段,
VPD=0.00002TIRi 2-0.00411TIRi-0.0056 R2=0.9901 (1)
在12:30~23:30时段,
VPD=-0.00004TIRi 2+0.02TIRi+0.2012 R2=0.9583 (2)
式中:VPD为室内冠层水汽压亏缺,TIRi为室内冠层太阳辐射。
可以看出,0:00~12:00时段回归模型的相关系数为0.9901,12:30~23:30时段回归模型的相关系数为0.9583,两个时段的相关系数都非常高,可见晴天条件下可以根据室内辐射量和经验参数来确定室内冠层水汽压亏缺VPD。
后期针对该温室大棚,只需测出不同时间点的冠层太阳辐射,将测得的冠层太阳辐射值代入相对应时段的回归模型,即可计算得到该时段的室内冠层水汽压差,不需要按照目前现有的方法进行繁琐的计算,简便易行。
以上所述仅为本发明的较佳实施例,并不限制本发明,凡在本发明的精神和原则范围内所做的任何修改、等同替换和改进,均应包含在本发明的保护范围之内。
Claims (4)
1.一种晴天条件下温室作物冠层饱和水汽压差的计算方法,其特征在于:包括以下步骤,
步骤一、采用植物生长监测仪采集温室0:00~12:00和12:30~23:30两个时段中不同时间点的日光温室大棚内的温室环境参数和温室内作物的冠层太阳辐射并实时监测作物的生长状况,并计算两个时段中不同时间点的冠层饱和水汽压差;
步骤二、分别对0:00~12:00以及12:30~23:30两时段的室内冠层太阳辐射与饱和水汽压亏缺的关系进行模拟得到两时段的温室内饱和水汽压亏缺与室内冠层太阳辐射的回归模型:
0:00~12:00时段,VPD=a TIRi 2+b TIRi+c
12:30~23:30时段,VPD=a’TIRi 2+b’TIRi+c‘
式中:VPD为室内冠层水汽压亏缺,TIRi为室内冠层太阳辐射,a、b、c、a’、b’、c‘为常数;
步骤三、将步骤一得到的两个时段的室内冠层太阳辐射分别代入步骤二中对应时段的回归模型中计算得到a、b、c、a’、b’、c‘常数数值,得到不同时段具体的回归方程;
步骤四、将测得的不同时段的室内冠层太阳辐射数值代入步骤三对应时段的回归方程中计算得到温室内的饱和水汽压亏缺。
2.根据权利要求1所述的一种晴天条件下温室作物冠层饱和水汽压差的计算方法,其特征在于:温室内温室环境参数和冠层太阳辐射的采集时间为每日的7:00、10:00、12:00、15:00和18:00。
3.根据权利要求1或2所述的一种晴天条件下温室作物冠层饱和水汽压差的计算方法,其特征在于:所述温室环境参数包括温室内作物冠层高处的空气温度和湿度。
4.根据权利要求1所述的一种晴天条件下温室作物冠层饱和水汽压差的计算方法,其特征在于:所述温室作物为日光温室番茄,
0:00~12:00时段,
VPD=0.00002TIRi 2-0.00411TIRi-0.0056 R2=0.9901;
12:30~23:30时段,
VPD=-0.00004TIRi 2+0.02TIRi+0.2012 R2=0.9583。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010183568.7A CN111460377B (zh) | 2020-03-16 | 2020-03-16 | 一种晴天条件下温室作物冠层饱和水汽压差计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010183568.7A CN111460377B (zh) | 2020-03-16 | 2020-03-16 | 一种晴天条件下温室作物冠层饱和水汽压差计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111460377A CN111460377A (zh) | 2020-07-28 |
CN111460377B true CN111460377B (zh) | 2023-03-24 |
Family
ID=71684319
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010183568.7A Active CN111460377B (zh) | 2020-03-16 | 2020-03-16 | 一种晴天条件下温室作物冠层饱和水汽压差计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111460377B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112734244B (zh) * | 2021-01-14 | 2021-09-14 | 中国科学院地理科学与资源研究所 | 一种基于饱和水汽压差的干旱指数计算方法 |
CN118068897A (zh) * | 2024-04-17 | 2024-05-24 | 北京天创金农科技有限公司 | 一种日光温室环境调控系统及方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7155377B2 (en) * | 2001-08-23 | 2006-12-26 | Wisconsin Alumni Research Foundation | Method and system for calculating the spatial-temporal effects of climate and other environmental conditions on animals |
US20130308675A1 (en) * | 2012-05-18 | 2013-11-21 | Smartfield, Inc. | Optimum plant canopy temperature |
CN103234691B (zh) * | 2013-04-24 | 2015-02-18 | 北京农业智能装备技术研究中心 | 一种温室植物叶气饱和水汽压差测量装置及测量方法 |
CN110083954B (zh) * | 2019-05-05 | 2023-03-24 | 华北水利水电大学 | 基于边界层阻力测量技术的温室作物需水量计算方法 |
-
2020
- 2020-03-16 CN CN202010183568.7A patent/CN111460377B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111460377A (zh) | 2020-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Grattan et al. | New crop coefficients estimate water use of vegetables, row crops | |
Van Iersel et al. | Sensors for improved efficiency of irrigation in greenhouse and nursery production | |
Hollinger et al. | Carbon dioxide exchange between an undisturbed old‐growth temperate forest and the atmosphere | |
Parmele et al. | Micrometeorological measurement of pesticide vapor flux from bare soil and corn under field conditions | |
CN111460377B (zh) | 一种晴天条件下温室作物冠层饱和水汽压差计算方法 | |
CN113049750B (zh) | 一种基于高通量气孔导度诊断植物水分胁迫的方法及系统 | |
Gao et al. | Canopy and leaf gas exchange of Haloxylon ammodendron under different soil moisture regimes | |
Zhang et al. | Estimation of cucumber evapotranspiration in solar greenhouse in Northeast China | |
CN109447426B (zh) | 基于作物需水机理的灌溉需水对变化环境的响应分析方法 | |
CN103207258B (zh) | 利用指示植物的需水信息确定被考察植物需水量的方法 | |
CN103308094A (zh) | 一种温室作物耗水量估算方法 | |
Holwerda et al. | Surface energy exchange in a tropical montane cloud forest environment: Flux partitioning, and seasonal and land cover-related variations | |
CN110432046B (zh) | 一种温室内的智能灌溉系统 | |
CN110210142B (zh) | 一种南方大型灌区水稻实时需水量测算方法 | |
Buttar et al. | Height effect of air temperature measurement on sensible heat flux estimation using flux variance method. | |
CN114365682A (zh) | 一种设施栽培土壤水分预测方法、装置及电子设备 | |
CN114418235A (zh) | 一种农作物生长状况确定方法及系统 | |
Balbontin-Nesvara et al. | Vineyard evapotranspiration measuraments in a semiarid environment: Eddy covariance and bowen ratio comparison | |
Treder et al. | An hourly reference evapotranspiration model as a tool for estimating plant water requirements | |
CN114642159B (zh) | 一种设施作物灌溉定额计算方法及系统 | |
Maltese et al. | Comparing actual evapotranspiration and plant water potential on a vineyard | |
Mashonjowa et al. | Measurement and simulation of the ventilation rates in a naturally ventilated Azrom-type greenhouse in Zimbabwe | |
CN109239027A (zh) | 一种作物生理性缺水诊断方法及系统 | |
CN114186425A (zh) | 一种冠层导度与气象变量相关性评估方法及系统 | |
CN113080035A (zh) | 农田作物需水量的确定方法、装置、灌水方法及灌水系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |