CN114200986A - 一种考虑作物生产效益和节省能耗的温室环境优化设计方法 - Google Patents
一种考虑作物生产效益和节省能耗的温室环境优化设计方法 Download PDFInfo
- Publication number
- CN114200986A CN114200986A CN202111367396.XA CN202111367396A CN114200986A CN 114200986 A CN114200986 A CN 114200986A CN 202111367396 A CN202111367396 A CN 202111367396A CN 114200986 A CN114200986 A CN 114200986A
- Authority
- CN
- China
- Prior art keywords
- crop
- greenhouse
- environment
- environmental
- parameter
- 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.)
- Granted
Links
- 230000008901 benefit Effects 0.000 title claims abstract description 41
- 238000005457 optimization Methods 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000012272 crop production Methods 0.000 title claims abstract description 25
- 230000007613 environmental effect Effects 0.000 claims abstract description 46
- 230000004044 response Effects 0.000 claims abstract description 44
- 238000005070 sampling Methods 0.000 claims abstract description 14
- 230000008859 change Effects 0.000 claims abstract description 11
- 230000009467 reduction Effects 0.000 claims abstract description 10
- 238000004088 simulation Methods 0.000 claims abstract description 9
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 7
- 238000004134 energy conservation Methods 0.000 claims abstract description 7
- 238000005516 engineering process Methods 0.000 claims abstract description 7
- 239000012530 fluid Substances 0.000 claims abstract description 7
- 230000003993 interaction Effects 0.000 claims abstract description 7
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims description 58
- 229910002092 carbon dioxide Inorganic materials 0.000 claims description 35
- 239000001569 carbon dioxide Substances 0.000 claims description 29
- 239000013598 vector Substances 0.000 claims description 23
- 230000005855 radiation Effects 0.000 claims description 15
- 239000011159 matrix material Substances 0.000 claims description 14
- 238000004422 calculation algorithm Methods 0.000 claims description 13
- 238000005265 energy consumption Methods 0.000 claims description 12
- 230000006870 function Effects 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 230000029553 photosynthesis Effects 0.000 claims description 7
- 238000010672 photosynthesis Methods 0.000 claims description 7
- 230000000243 photosynthetic effect Effects 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 238000009826 distribution Methods 0.000 claims description 5
- 238000004519 manufacturing process Methods 0.000 claims description 5
- 241000208822 Lactuca Species 0.000 claims description 4
- 235000003228 Lactuca sativa Nutrition 0.000 claims description 4
- 230000005611 electricity Effects 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 230000004907 flux Effects 0.000 claims description 3
- 238000005286 illumination Methods 0.000 claims description 3
- 230000000241 respiratory effect Effects 0.000 claims description 3
- 230000036387 respiratory rate Effects 0.000 claims description 3
- 230000005068 transpiration Effects 0.000 claims description 3
- 230000009471 action Effects 0.000 claims description 2
- 238000002347 injection Methods 0.000 claims description 2
- 239000007924 injection Substances 0.000 claims description 2
- 238000006467 substitution reaction Methods 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 241000196324 Embryophyta Species 0.000 description 1
- 241000238631 Hexapoda Species 0.000 description 1
- 241000607479 Yersinia pestis Species 0.000 description 1
- 238000012271 agricultural production Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000005329 float glass Substances 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D27/00—Simultaneous control of variables covered by two or more of main groups G05D1/00 - G05D25/00
- G05D27/02—Simultaneous control of variables covered by two or more of main groups G05D1/00 - G05D25/00 characterised by the use of electric means
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Mining
-
- A—HUMAN NECESSITIES
- A01—AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
- A01G—HORTICULTURE; CULTIVATION OF VEGETABLES, FLOWERS, RICE, FRUIT, VINES, HOPS OR SEAWEED; FORESTRY; WATERING
- A01G2/00—Vegetative propagation
-
- 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
- Y02A40/00—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
- Y02A40/10—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in agriculture
- Y02A40/25—Greenhouse technology, e.g. cooling systems therefor
-
- 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
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P60/00—Technologies relating to agriculture, livestock or agroalimentary industries
- Y02P60/14—Measures for saving energy, e.g. in green houses
Abstract
本发明公开了一种考虑作物生产效益和节省能耗的温室环境优化设计方法，根据温室实际确定关键设计参数及可行范围；使用计算流体力学(CFD)工具建立温室环境的三维模型；构建多孔介质模型，模拟作物与周围环境的热交互；在可行范围内，对设计参数多维采样进行CFD仿真；提取环境参数稳态响应结果，利用奇异值分解技术重构环境替代模型；建立作物生长模型，设置作物生长效益指标和环境参数的均匀性指标；利用环境替代模型和多维插值快速解算环境响应，实施关键设计参数的多目标优化。本发明通过引入作物模型来测算候选环境参数引起的相关环境变化及经济效益，使得温室环境的优化设计更具针对性，能够最大化作物生产效益，同时促进设施农业节能减排。
Description
技术领域
本发明涉及设施农业生产领域，具体涉及一种考虑作物生产效益和节省能耗的温室环境优化设计方法。
背景技术
作为一种半封闭生态系统，温室能减轻传统农业受外界环境的依赖，但其自我调节能力有限，内部环境的优劣对于作物生长、病虫害防治、节能减排等影响显著。针对温室内环境进行合理布局与优化能有效改善作物生长条件、促进作物增产同时节约能耗，对于促进设施农业进一步发展具有重要的现实意义和实际应用价值。
当前农业温室优化设计研究大多假设环境参数均匀分布，用集总参数来评价作物的生长环境。例如中国专利公开号CN113079881A的文献公开的大跨度可变空间日光温室优化设计方法，通过计算能量平衡方程得到冬季温室温度及供热负荷，进而推算出集热系统及墙体内空气温度场等热性能参数。其未考虑温室内部环境参数的空间变化，以及调控设备的优化布局；中国专利授权公告号CN107545100B的文献公开的温室环境高分辨率降阶建模方法，通过本征正交分解(POD)技术构造温室环境特征向量，从而提取出高分辨率的温室环境时空信息，但其没有与作物生长模型结合，导致由此专利衍生出的温室优化策略不能达到作物生产效益与节能效果的最优。
发明内容
本发明的目的是解决目前农业温室环境设计方法过度依赖专家经验，导致优化精度低、设计结果并非经济效益最优。提出一种考虑作物生产效益和节省能耗的农业温室优化设计方法。特点在于：1.将基于奇异值分解(SVD)的模型降阶技术与计算流体力学(CFD)稳态仿真相结合，快速求取温室作物区域的高分辨率环境响应，提高优化设计的空间分辨率；2.考虑温室环境与作物之间的交互作用，通过建立作物生长模型，使温室环境设计能最大化作物生产效益，同时促进节能减排。
本发明采用的技术方案是：
一种考虑作物生产效益和节省能耗的温室环境优化设计方法，包括以下步骤：
(1)根据温室外部天气条件和围护结构特性，确定影响温室环境的关键设计参数及可行范围。
(2)在可行范围内，对设计参数进行多维采样，每个采样点组成一个参数向量。
(3)根据步骤(1)所述外部天气条件与温室围护结构，以及步骤(2)所述参数向量，使用计算流体力学工具建立温室环境的三维模型。
(4)温室内部的作物区域由多孔介质构建，用于模拟作物与周围环境的热交互作用。
(5)针对每个参数向量分别运行CFD稳态仿真。
(6)提取每次仿真结果中作物区域的环境参数稳态响应值，包括温湿度场、二氧化碳分布、风速流场等，组成对应的环境响应变化空间。
(7)利用奇异值分解技术对上述空间进行降维，重构出对应低维子空间，用于环境响应的快速解算。
(8)建立作物生长模型，用于求取各个环境响应下温室作物的生产效益。
(9)设置作物区域的环境均匀性指标，并结合作物生长效益指标建立优化设计的目标函数。
(10)根据步骤(9)所述目标函数，利用全局优化算法对步骤(1)所述设计参数进行多目标优化。每次迭代过程的温室环境响应由步骤(7)所得的低维子空间通过多维插值快速解算。
所述步骤(1)中，温室外部天气条件包括温湿度、光照强度和风速；影响温室作物生长环境的设计参数包括顶棚遮阳率范围、风机位置和二氧化碳加注位置等。
所述步骤(2)中，为保证环境参数模型降维的精度，参数采样方法采用多维等间距采样，见附图3。
所述步骤(3)中，CFD工具采用ANSYS软件；温室环境建模的具体步骤包括：
A、采用ANSYS里的DesignModeler软件对温室结构进行几何建模；根据步骤(2)所述参数向量对温室顶棚遮阳帘、负压风机和二氧化碳加注设备分别进行参数设置；
B、将几何模型导入mesh软件进行网格划分，并将划分后的网格文件导入fluent软件；
C、设置CFD模型边界条件；选择动量、湍流和DO辐射模型模拟温室内的热对流与热辐射现象。
所述步骤(4)中，设置多孔介质构建温室作物区域，运用显热、潜热模型模拟作物与环境之间的热交互作用。具体步骤包括：
A、由于温室中作物冠层与室内空气存在温差，作物-环境之间显热交换量Qsem由以下作物冠层空气动力学特性决定：
式中ρ为温室内空气密度，cp为空气的比热容，LAI为作物叶面积指数，Tc为冠层温度，Ti为室温，ra为作物叶片空气动力学阻力。
B、作物将太阳辐射产生的显热转换为潜热释放到温室，潜热量Qlat可用下式表示：
式中，λ为蒸腾通量，wleaf，wi为作物冠层和空气的绝对湿度，rsl为作物叶片的气孔阻力。
所述步骤(6)中，提取作物区域环境参数稳态响应值的具体步骤包括：
A、根据温室作物生长区域，选取作物冠层高度的平面，提取平面上所有网格点的环境响应值；
B、提取方法：通过fluent\export导出选定区域的物理量，根据导出的profile文件格式整理得到各环境参数的响应矩阵；
所述步骤(7)中，采用SVD技术对环境响应变化空间进行降维，具体步骤包括：
A、使用步骤(6)中的环境响应矩阵组成参数矩阵R。其中，R为一个n×m阶的矩阵，n为温室内部作物区域选定的网格点个数，m为参数向量的数目。
B、利用特征值方程(RTR)vi＝λivi解出n个特征值与对应的n个特征向量vi，所有特征向量张成一个n×n的矩阵V。利用(RRT)ui＝λiui解出n个特征值与对应的n个特征向量ui。
C、利用σi＝RVi/ui解出每个奇异值，对应的SVD基表示为：
D、选择合适的截断值k，使得前k个奇异值包含的系统特征值占总能的比值达到99％以上。由此，环境响应变化空间可近似表述为k个奇异向量
及其系数a的线性组合
其中P(n)表示n个网格点上的稳态环境响应结果。
所述步骤(8)中，采用作物生长模型中的干物质质量表征作物的生产效益，时间段选取作物日间光合作用时段。以生菜类作物模型为例，表征干物质质量的一阶差分方程如下：
其中，xd为作物干物质质量，Cαβ为作物产量因子，
是光合反应速率，Cresp,d是以呼吸干物质量表示的呼吸速率参数；zt为作物区域平均温度，由CFD环境响应得到；k为时间步。上式中的光合反应速率
由下式计算：
其中，Cpl,d是单位质量下的作物有效冠层面积参数(m2kg-1)，crad,phot是光利用效率参数(kgJ-1)；vrad为太阳辐射量(Wm-2)，由CFD室外太阳辐射和温室遮阳率变换得到；
和
分别表示温度对光合作用的影响参数；zc(k)是作物区域的平均二氧化碳浓度，由CFD环境响应得到；cr是二氧化碳补偿点参数(kgm-3)。
所述步骤(9)中，作物区域的生长环境均匀性指标Jeven由作物区域网格点处的温度、湿度和二氧化碳浓度值的方差和表征；作物生长效益指标由作物生产的经济收益减去风机电能耗和二氧化碳的加注能耗。为统一单位，收益和能耗都折算为市场价格，以指标Jecon表示。由此，优化设计的目标函数集为：
生长环境均匀性指标Jeven中，Ti，Hi，Ci表示作物区域网格点处的温度、湿度和CO2浓度，Zt，Zh，Zc为其平均值；Np为网格点数目；作物生长效益指标Jecon中，Prcxd为作物市场价格，Vair,k为风机送风速率(m3/s)，ηfan为风机效率，ΔP为温室进出口压强差(Pa)，Prce为当地电价；VCO2为二氧化碳加注速率(m3/s)，ωρ分别为CO2的质量分数和密度，Prcco2为CO2加注价格，Thour为作用时间。
所述步骤(10)中，优化算法采用多目标MOEA/D优化算法，优化目标为作物生长环境均匀性指标Jeven最小，作物生长效益指标Jecon最大。算法种群数设置为500，迭代次数设置为70。每次迭代过程的温室环境响应通过步骤(7)对奇异向量系数的多维插值快速解算得到。
本发明采用上述技术方案后的有益效果是：
1、本发明考虑作物生产效益和节省能耗的温室环境优化设计方法通过引入作物模型来测算候选环境参数引起的作物干物质质量及经济效益，同时通过折算温室运行能耗的经济代价，使得温室环境的优化设计更具针对性，优化结果有经济效益数据作支持，能够最大化作物生产效益，同时促进设施农业节能减排。
2、本发明考虑温室环境的空间分布影响，采用模型降阶的思想快速提取环境特征信息，使得本发明的温室环境设计方法能够保证最大化作物区域的环境适宜性，同时缩短优化设计的时长，提高优化效率。
附图说明
图1是本发明考虑作物生产效益和节省能耗的温室环境优化设计流程图；
图2是一个三维多跨温室模型示意图。图中绿色区域模拟为作物生长区域，二氧化碳加注的候选位置有三个，标注为圆圈1、2、3；负压风机的候选位置由相互距离h决定；
图3是多维设计参数的采样点分布示意图。
具体实施方式
为了更为具体的描述本发明，下面结合附图和具体实施案例对本发明进行详细说明。
图1描述了考虑作物生产效益和节省能耗的温室环境优化设计流程。
图2为一个多跨温室三维模型。温室长40m，宽19m，屋顶高5m。温室内部有一宽8m、长20m、占地160m2的作物生长区域。
下面详细描述本发明方法的实施步骤：
步骤1，根据温室外部温湿度、光照强度和风速等气候条件以及围护结构特性，确定影响温室环境的关键设计参数及可行范围。本例中，在温室几何模型已知的情况下，选取对温室作物生长影响较大的设计参数为：顶棚遮阳率、风机位置和二氧化碳加注位置。
步骤2，在可行范围内，对上述设计参数进行多维等间距采样，具体包括：遮阳率采样值为0.4、0.6、0.8，风机位置间隔(即附图2中风机相互距离h)采样值为1m、3m、5m，二氧化碳加注位置间隔(与南墙距离)采样值为3m、9m、15m，每组采样点组成一个参数向量参与CFD仿真，如附图3所示。
步骤3，根据步骤1所述外部天气条件与温室围护结构，使用计算流体力学工具建立温室环境的三维模型。具体步骤包括：
步骤3.1，确定温室、湿帘入口以及作物区域的位置和尺寸，利用DesignModeler软件搭建温室几何模型，如附图2所示；根据步骤2所述参数向量对温室顶棚遮阳帘、负压风机和二氧化碳加注设备分别进行参数设置；
步骤3.2，将几何模型导入mesh软件进行网格划分，对特殊边界(如湿帘入口、负压风机和二氧化碳加注口)局部加密。本例中共划分出182491个网格，并将网格文件导入fluent软件；
步骤3.3，设置CFD模型边界条件，包括：四周与屋顶的围护设置为温度边界(浮法玻璃材质)，风机设置为速度进口边界，湿帘设置为压力入口边界；温室内气体假设为低速流动的不可压缩粘性牛顿流体，湍流模型选用标准k-ε模型，近壁处理采用标准壁面函数，辐射模型采用DO辐射模型；使用组分运输方程以模拟二氧化碳等多组分气体的流动，采用压力速度耦合方式，采用SIMPLE算法求解。
步骤4，采用多孔介质构建温室作物区域，运用显热、潜热模型模拟作物与环境之间的热交互作用。具体步骤包括：
步骤4.1，根据作物特点设置孔隙率；由于温室中作物冠层与室内空气存在温差，作物-环境之间显热交换量Qsem由以下作物冠层空气动力学特性决定：
式中ρ为温室内空气密度，cp为空气的比热容，LAI为作物叶面积指数，Tc为冠层温度，Ti为室温，ra为作物叶片空气动力学阻力；
步骤4.2，作物将太阳辐射产生的显热转换为潜热释放到温室，潜热量Qlat用下式表示：
式中，λ为蒸腾通量，wleaf，wi为作物冠层和空气的绝对湿度，rsl为作物叶片的气孔阻力。
步骤5，根据步骤2所述，本例候选参数向量共3×3×3＝27个；根据步骤3和步骤4所述，分别针对每个参数向量进行CFD设置，并完成离线的CFD稳态仿真。
步骤6，提取每次仿真结果中作物区域的环境参数稳态响应值，包括温度场、二氧化碳分布、风速流场等，组成对应的环境响应变化空间。具体步骤包括：
步骤6.1，根据温室作物生长区域，选取作物冠层平面高度为1.3米，提取平面上所有网格点的环境响应值，本例网格点数为4500；
步骤6.2，提取方法：通过fluent\export导出选定区域的物理量，根据导出的profile文件格式整理得到各环境参数的响应矩阵，并组成变化空间。
步骤7，利用奇异值分解技术对上述空间进行降维，重构出对应低维子空间，具体步骤包括：
步骤7.1，使用步骤6中的环境响应变化空间组成参数矩阵R。其中，R为一个n×m阶的矩阵，n为温室内部作物区域选定的网格点个数，m为参数向量的数目。本例中，n＝4500，m＝27；
步骤7.2，利用(RTR)vi＝λivi解出n个特征值与对应的n个特征向量vi，所有特征向量张成一个n×n的矩阵V。利用(RRT)ui＝λiui解出n个特征值与对应的n个特征向量ui；
步骤7.3，利用σi＝RVi/ui解出每个奇异值，对应的SVD基表示为：
步骤7.4，选择合适的截断值k，使得前k个奇异值包含的系统特征值占总能的比值I(k)达到99％以上，本例中k＝5。由此，环境响应变化空间可近似表述为k个奇异向量
及其系数a的线性组合
其中P(n)表示n个网格点上的稳态环境响应结果。
步骤8，建立作物生长模型，采用作物生长模型中的干物质质量表征作物产量，以此求取各个环境响应下温室作物的生产效益。本例将一周内光合作用时段作物干物质产生的累计量表征作物产量。以生菜类作物模型为例，关于干物质质量产生的一阶差分方程如下：
其中，xd为作物干物质质量；Cαβ为作物产量因子，取值0.544；
是光合反应速率；Cresp,d是以呼吸干物质量表示的呼吸速率参数，取值2.65×10-7s-1；zt为作物区域平均温度，由CFD环境响应得到；k为时间步。进一步，上式中的光合反应速率
由下式计算：
其中，Cpl,d是单位质量下的作物有效冠层面积参数(m2kg-1)，取值53；crad,phot是光利用效率参数(kgJ-1)，取值3.55×10-9；vrad为太阳辐射量(Wm-2)，由CFD室外太阳辐射和温室遮阳率变换得到；温度对光合作用的影响参数：
取值5.11×10-6，
取值2.30×10-4，
取值6.29×10-4；zc(k)是作物区域的平均二氧化碳浓度，由CFD环境响应得到；cr是二氧化碳补偿点参数(kgm-3)，取值5.2×10-5。
步骤9，设置作物区域的环境均匀性指标，并结合作物生长效益指标建立优化设计的目标函数。其中，生长环境均匀性指标Jeven由作物区域网格点处的温度、湿度和二氧化碳浓度值的方差和表征；作物生长效益指标由作物生产的经济收益减去风机电能耗和二氧化碳的加注能耗。为统一单位，收益和能耗都折算为市场价格，以指标Jecon表示。由此，优化设计的目标函数为：
在生长环境均匀性指标Jeven中，Ti，Hi，Ci表示作物区域网格点处的温度、湿度和CO2浓度，Zt，Zh，Zc为其平均值，均可由CFD结果计算得到；Np为网格点数目，本例为4600；作物生长效益指标Jecon中，Prcxd为作物市场价格；Vair,k为风机送风速率(m/s)，本例设为3；ηfan为风机效率，本例设为0.75；ΔP为温室进出口压强差(Pa)，本例设为562.5；Prce为当地电价；VCO2为二氧化碳加注速率(m/s)，本例设为0.02；ωρ分别为CO2的质量分数和密度，Prcco2为CO2加注价格；本例考虑一周光合作用时段内作物干物质的产生量，作用时间Thour取值7×8＝56小时。
步骤10，根据步骤9所述目标函数，利用全局优化算法对步骤1所述设计参数进行多目标优化。为提高优化速度，每次迭代过程的温室环境响应通过步骤7对奇异向量系数a的多维插值解算得到。优化算法可采用多目标MOEA/D优化算法，优化目标为作物生长环境均匀性指标Jeven最小，作物生长效益指标Jecon最大。算法种群数设置为500，迭代次数设置为70。
为了验证本发明效果，以华东地区生菜种植的温室环境为例，优化对温室作物生长影响较大的三个设计参数，即顶棚遮阳率范围、风机位置和二氧化碳加注位置。运用此方法可得到最佳温室环境设计参数，使得作物区域内的作物生长效益最优且维持温室的节能运行。
Claims (10)
1.一种考虑作物生产效益和节省能耗的温室环境优化设计方法，其特征在于，包括以下步骤：
(1)根据温室外部天气条件和围护结构特性，确定影响温室环境的关键设计参数及可行范围；
(2)在可行范围内，对设计参数进行多维采样，每个采样点组成一个参数向量；
(3)根据步骤(1)所述外部天气条件与温室围护结构，以及步骤(2)所述参数向量，使用计算流体力学工具CFD建立温室环境的三维模型；
(4)温室内部的作物区域由多孔介质构建，用于模拟作物与周围环境的热交互作用；
(5)针对每个参数向量分别运行CFD稳态仿真；
(6)提取每次仿真结果中作物区域的环境参数稳态响应值，包括温湿度场、二氧化碳分布、风速流场等，组成对应的环境响应变化空间；
(7)利用奇异值分解技术对上述空间进行降维，重构出对应低维子空间，用于环境响应的快速解算；
(8)建立作物生长模型，用于求取各个环境响应下温室作物的生产效益；
(9)设置作物区域的环境均匀性指标，并结合作物生长效益指标建立优化设计的目标函数；
(10)根据步骤(9)所述目标函数，利用全局优化算法对步骤(1)所述设计参数进行多目标优化。每次迭代过程的温室环境响应由步骤(7)所得的低维子空间通过多维插值快速解算。
2.根据权利要求1所述的一种考虑作物生产效益和节省能耗的温室环境优化设计方法，其特征在于，所述步骤(1)中，温室外部天气条件包括温湿度、光照强度和风速；影响温室作物生长环境的设计参数包括顶棚遮阳率范围、风机位置和二氧化碳加注位置。
3.根据权利要求1所述的一种考虑作物生产效益和节省能耗的温室环境优化设计方法，其特征在于，所述步骤(2)中，为保证环境参数模型降维的精度，参数采样方法采用多维等间距采样。
4.根据权利要求1所述的一种考虑作物生产效益和节省能耗的温室环境优化设计方法，其特征在于，所述步骤(3)中，计算流体力学工具采用ANSYS软件；温室环境建模的具体步骤包括：
A、采用ANSYS里的DesignModeler软件对温室结构进行几何建模；根据步骤(2)所述参数向量对温室顶棚遮阳帘、负压风机和二氧化碳加注设备分别进行参数设置；
B、将几何模型导入mesh软件进行网格划分，并将划分后的网格文件导入fluent软件；
C、设置CFD模型边界条件；选择动量、湍流和DO辐射模型模拟温室内的热对流与热辐射现象。
5.根据权利要求1所述的一种考虑作物生产效益和节省能耗的温室环境优化设计方法，其特征在于，所述步骤(4)中，设置多孔介质构建温室作物区域，运用显热、潜热模型模拟作物与环境之间的热交互作用，具体步骤包括：
4.1、由于温室中作物冠层与室内空气存在温差，作物-环境之间显热交换量Qsem由以下作物冠层空气动力学特性决定：
式中ρ为温室内空气密度，cp为空气的比热容，LAI为作物叶面积指数，Tc为冠层温度，Ti为室温，ra为作物叶片空气动力学阻力；
4.2、作物将太阳辐射产生的显热转换为潜热释放到温室，潜热量Qlat可用下式表示：
式中，λ为蒸腾通量，wleaf，wi为作物冠层和空气的绝对湿度，rsl为作物叶片的气孔阻力。
6.根据权利要求1所述的一种考虑作物生产效益和节省能耗的温室环境优化设计方法，其特征在于，所述步骤(6)中，提取作物区域环境参数稳态响应值的具体步骤包括：
6.1、根据温室作物生长区域，选取作物冠层高度的平面，提取平面上所有网格点的环境响应值；
6.2、提取方法：通过fluent\export导出选定区域的物理量，根据导出的profile文件格式整理得到各环境参数的响应矩阵。
7.根据权利要求1所述的一种考虑作物生产效益和节省能耗的温室环境优化设计方法，其特征在于，所述步骤(7)中，采用奇异值分解SVD技术对环境响应变化空间进行降维，具体步骤包括：
7.1、使用步骤(6)中的环境响应矩阵组成参数矩阵R，其中，R为一个n×m阶的矩阵，n为温室内部作物区域选定的网格点个数，m为参数向量的数目；
7.2、利用特征值方程(RTR)vi＝λivi解出n个特征值与对应的n个特征向量vi，所有特征向量张成一个n×n的矩阵V；利用(RRT)ui＝λiui解出n个特征值与对应的n个特征向量ui；
7.3、利用σi＝RVi/ui解出每个奇异值，对应的SVD基表示为：
7.4、选择合适的截断值k，使得前k个奇异值包含的系统特征值占总能的比值达到99％以上；由此，环境响应变化空间可近似表述为k个奇异向量
及其系数a的线性组合；
其中P(n)表示n个网格点上的稳态环境响应结果。
8.根据权利要求1所述的一种考虑作物生产效益和节省能耗的温室环境优化设计方法，其特征在于，所述步骤(8)中，采用作物生长模型中的干物质质量表征作物的生产效益，时间段选取作物日间光合作用时段，以生菜类作物模型为例，表征干物质质量的一阶差分方程如下：
其中，xd为作物干物质质量，Cαβ为作物产量因子，
是光合反应速率，Cresp,d是以呼吸干物质量表示的呼吸速率参数；zt为作物区域平均温度，由CFD环境响应得到；k为时间步，上式中的光合反应速率
由下式计算：
其中，Cpl,d是单位质量下的作物有效冠层面积参数(m2kg-1)，crad,phot是光利用效率参数(kgJ-1)；vrad为太阳辐射量(Wm-2)，由CFD室外太阳辐射和温室遮阳率变换得到；
和
分别表示温度对光合作用的影响参数；zc(k)是作物区域的平均二氧化碳浓度，由CFD环境响应得到；cr是二氧化碳补偿点参数(kgm-3)。
9.根据权利要求1所述的一种考虑作物生产效益和节省能耗的温室环境优化设计方法，其特征在于，所述步骤(9)中，作物区域的生长环境均匀性指标Jeven由作物区域网格点处的温度、湿度和二氧化碳浓度值的方差和表征；作物生长效益指标由作物生产的经济收益减去风机电能耗和二氧化碳的加注能耗,为统一单位，收益和能耗都折算为市场价格，以指标Jecon表示,由此，优化设计的目标函数集为：
生长环境均匀性指标Jeven中，Ti，Hi，Ci表示作物区域网格点处的温度、湿度和CO2浓度，Zt，Zh，Zc为其平均值；Np为网格点数目；作物生长效益指标Jecon中，Prcxd为作物市场价格，Vair,k为风机送风速率(m3/s)，ηfan为风机效率，ΔP为温室进出口压强差(Pa)，Prce为当地电价；VCO2为二氧化碳加注速率(m3/s)，ωρ分别为CO2的质量分数和密度，Prcco2为CO2加注价格，Thour为作用时间。
10.根据权利要求1所述的一种考虑作物生产效益和节省能耗的温室环境优化设计方法，其特征在于，所述步骤(10)中，优化算法采用多目标MOEA/D优化算法，优化目标为作物生长环境均匀性指标Jeven最小，作物生长效益指标Jecon最大；算法种群数设置为500，迭代次数设置为70,每次迭代过程的温室环境响应通过步骤(7)对奇异向量系数的多维插值快速解算得到。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111367396.XA CN114200986B (zh) | 2021-11-18 | 2021-11-18 | 一种考虑作物生产效益和节省能耗的温室环境优化设计方法 |
| PCT/CN2021/140329 WO2023087483A1 (zh) | 2021-11-18 | 2021-12-22 | 一种考虑作物生产效益和节省能耗的温室环境优化设计方法 |
| GB2317885.8A GB2621942A (en) | 2021-11-18 | 2021-12-22 | Greenhouse environment optimization design method considering crop production benefits and energy consumption reduction |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111367396.XA CN114200986B (zh) | 2021-11-18 | 2021-11-18 | 一种考虑作物生产效益和节省能耗的温室环境优化设计方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114200986A true CN114200986A (zh) | 2022-03-18 |
| CN114200986B CN114200986B (zh) | 2023-03-21 |
Family
ID=80648035
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111367396.XA Active CN114200986B (zh) | 2021-11-18 | 2021-11-18 | 一种考虑作物生产效益和节省能耗的温室环境优化设计方法 |
Country Status (3)
| Country | Link |
|---|---|
| CN (1) | CN114200986B (zh) |
| GB (1) | GB2621942A (zh) |
| WO (1) | WO2023087483A1 (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115265633A (zh) * | 2022-09-05 | 2022-11-01 | 中车长春轨道客车股份有限公司 | 一种低温恒温模型工作状态的确定方法及装置 |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116449897B (zh) * | 2023-06-08 | 2023-10-20 | 北京市农林科学院智能装备技术研究中心 | 温室环境最优控制方法、服务器及系统 |
| CN116451597B (zh) * | 2023-06-19 | 2023-08-18 | 广东电网有限责任公司佛山供电局 | 一种变电站高压室内控温优化方法与系统 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2130426A2 (de) * | 2008-06-05 | 2009-12-09 | Institut für Gemüse- & Zierpflanzenbau e.V. | Verfahren und Anordnung zur Kohlendioxid-Versorgung von in einem Gewächshaus angebauten Pflanzen |
| CN103049612A (zh) * | 2012-12-24 | 2013-04-17 | 江苏大学 | 一种基于模型降阶技术的建筑室内环境优化方法 |
| CN103884075A (zh) * | 2014-01-06 | 2014-06-25 | 浙江工业大学 | 基于计算流体力学和能量预测混合的温室节能控制方法 |
| CN106842923A (zh) * | 2017-01-17 | 2017-06-13 | 同济大学 | 基于作物生理和能耗优化的温室环境多因子协调控制方法 |
| CN107545100A (zh) * | 2017-07-31 | 2018-01-05 | 江苏大学 | 一种venlo型温室环境的高分辨率降阶建模方法 |
| CN107562996A (zh) * | 2017-07-31 | 2018-01-09 | 江苏大学 | 一种Venlo型温室环境参数的多目标集成优化方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102937785A (zh) * | 2012-11-09 | 2013-02-20 | 江苏大学 | 一种温室多模型信息融合的温室环境调控方法及系统 |
| CN105159257B (zh) * | 2015-08-25 | 2017-09-19 | 浙江大学 | 一种植物工厂集成控制系统及方法 |
| US20170094920A1 (en) * | 2015-10-02 | 2017-04-06 | Craig Ellins | Integrated incubation, cultivation and curing system and controls for optimizing and enhancing plant growth, development and performance of plant-based medical therapies |
| CN107045367B (zh) * | 2017-04-20 | 2019-01-25 | 同济大学 | 一种温室环境多因子协调节能优化控制方法 |
-
2021
- 2021-11-18 CN CN202111367396.XA patent/CN114200986B/zh active Active
- 2021-12-22 GB GB2317885.8A patent/GB2621942A/en active Pending
- 2021-12-22 WO PCT/CN2021/140329 patent/WO2023087483A1/zh unknown
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2130426A2 (de) * | 2008-06-05 | 2009-12-09 | Institut für Gemüse- & Zierpflanzenbau e.V. | Verfahren und Anordnung zur Kohlendioxid-Versorgung von in einem Gewächshaus angebauten Pflanzen |
| CN103049612A (zh) * | 2012-12-24 | 2013-04-17 | 江苏大学 | 一种基于模型降阶技术的建筑室内环境优化方法 |
| CN103884075A (zh) * | 2014-01-06 | 2014-06-25 | 浙江工业大学 | 基于计算流体力学和能量预测混合的温室节能控制方法 |
| CN106842923A (zh) * | 2017-01-17 | 2017-06-13 | 同济大学 | 基于作物生理和能耗优化的温室环境多因子协调控制方法 |
| CN107545100A (zh) * | 2017-07-31 | 2018-01-05 | 江苏大学 | 一种venlo型温室环境的高分辨率降阶建模方法 |
| CN107562996A (zh) * | 2017-07-31 | 2018-01-09 | 江苏大学 | 一种Venlo型温室环境参数的多目标集成优化方法 |
Non-Patent Citations (3)
| Title |
|---|
| 李康吉: "建筑室内环境建模、控制与优化及能耗预测", 《浙江大学博士学位论文》 * |
| 王定成等: "基于作物响应的温室环境SVMR控制仿真", 《农业工程学报》 * |
| 胥芳等: "湿帘-风机降温下的温室热/流场模拟及降温系统参数优化", 《农业工程学报》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115265633A (zh) * | 2022-09-05 | 2022-11-01 | 中车长春轨道客车股份有限公司 | 一种低温恒温模型工作状态的确定方法及装置 |
| CN115265633B (zh) * | 2022-09-05 | 2023-10-27 | 中车长春轨道客车股份有限公司 | 一种低温恒温模型工作状态的确定方法及装置 |
| WO2024050985A1 (zh) * | 2022-09-05 | 2024-03-14 | 中车长春轨道客车股份有限公司 | 一种低温恒温模型工作状态的确定方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| GB2621942A (en) | 2024-02-28 |
| WO2023087483A1 (zh) | 2023-05-25 |
| CN114200986B (zh) | 2023-03-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114200986B (zh) | 一种考虑作物生产效益和节省能耗的温室环境优化设计方法 | |
| Kacira et al. | Optimization of vent configuration by evaluating greenhouse and plant canopy ventilation rates under wind-induced ventilation | |
| CN107545100B (zh) | 一种venlo型温室环境的高分辨率降阶建模方法 | |
| CN105117787B (zh) | 一种温室能耗的模型优化预测方法 | |
| Xu et al. | Optimal control of Chinese solar greenhouse cultivation | |
| CN106842914A (zh) | 一种温度控制节能处理方法、装置及系统 | |
| Lee et al. | Assessing the performance potential of climate adaptive greenhouse shells | |
| CN113203187A (zh) | 基于部分线性模型的建筑暖通空调负荷优化控制方法 | |
| Campen et al. | Greenhouse cooling and heat recovery using fine wire heat exchangers in a closed pot plant greenhouse: design of an energy producing greenhouse | |
| CN105373849A (zh) | 光伏电池组件温度的分步预测方法 | |
| Wu et al. | A systematic analysis of multiple structural parameters of Chinese solar greenhouse based on the thermal performance | |
| CN116992548B (zh) | 一种建筑及区域综合能源系统的设计优化方法及系统 | |
| CN103903068B (zh) | 一种混合优化算法的温室能量预测方法 | |
| Altes-Buch et al. | A modeling framework for the integration of electrical and thermal energy systems in greenhouses | |
| Cemek et al. | Evaluation of temperature distribution in different greenhouse models using computational fluid dynamics (CFD). | |
| CN103823976B (zh) | 日光温室光热环境计算方法 | |
| Lyu et al. | Effects of vent opening, wind speed, and crop height on microenvironment in three-span arched greenhouse under natural ventilation | |
| Fu et al. | A temperature and vent opening couple model in solar greenhouses for vegetable cultivation based on dynamic solar heat load using computational fluid dynamics simulations | |
| CN107562996A (zh) | 一种Venlo型温室环境参数的多目标集成优化方法 | |
| CN112931297B (zh) | 基于计算机仿真的密闭式蛋鸡育雏舍小气候模型搭建方法 | |
| CN114662423A (zh) | 一种风机排布优化方法 | |
| Van Henten et al. | Improving the efficiency of greenhouse climate control: an optimal control approach | |
| Zhou et al. | Multi-objective optimization of fan-pad system operation for venlo greenhouse using CFD model based data interactive mechanism | |
| Wei et al. | CFD simulation and optimization of ventilation for the layout of community architecture inspired by fishbone form. | |
| Labihi et al. | Transient Thermal Simulation of a Greenhouse for Tomato Plants: A TRNSYS Model with Experimental Validation |
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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20240123 Address after: 230000 Woye Garden Commercial Building B-1017, 81 Ganquan Road, Shushan District, Hefei City, Anhui Province Patentee after: HEFEI JINGLONG ENVIRONMENTAL PROTECTION TECHNOLOGY Co.,Ltd. Country or region after: China Address before: Zhenjiang City, Jiangsu Province, 212013 Jingkou District Road No. 301 Patentee before: JIANGSU University Country or region before: China |