基于CFD的发酵床网养鸭舍冬季环境模拟与优化

吴照学 梁伟 鲍恩财 陈菁 柏宗春 应诗家

摘要：针对冬季规模化养殖鸭群集聚问题，对发酵床网养鸭舍温度、相对湿度和风速测试，构建鸭舍CFD模型，开展温度场、湿度场和气流场的环境模拟，提出鸭舍结构与装备的优化方案。结果表明：测试期间，舍内温度为9.80℃～17.68℃，舍内外平均温差8.29℃；舍内平均相对湿度为81.90%，略低于舍外；舍内风速为0.04～0.57m/s。对原模型温度场、湿度场和气流场进行模拟，与实测结果进行对比，其平均相对误差为2.96%～7.67%。进一步将通风小窗开启角度从45°增大至50°，在位于湿帘侧的发酵床初始端和发酵床中心位置增设2台暖风机对原有鸭舍环境优化，模拟发现优化后模型的舍内温度从15.32℃提高到16.85℃，达到育成期蛋鸭生长最适温度；舍内平均湿度降至66.52%，较原来下降11.07%，低于高湿环境阈值；舍内平均风速从0.29m/s提高到0.38m/s，提高排除有害气体的能力。为同类型鸭舍改善舍内环境质量提供理论依据。

关键词：鸭舍；发酵床网养；计算流体力学；环境优化

中图分类号：S834

文献标识码：A

文章编号：20955553 （2023） 12003307

Simulation and optimization of winter environment for duck house with

fermentation bed net based on CFD

Wu Zhaoxue1， Liang Wei1， 2， Bao Encai2， Chen Jing3， Bai Zongchun2， Ying Shijia2

（1. School of Engineering， Anhui Agricultural University， Hefei， 230036， China; 2. Jiangsu Academy of

Agricultural Sciences， Nanjing， 210014， China; 3. Jiangsu Gaoyou Duck Development Group Co.， Ltd.，

Yangzhou， 211200， China）

Abstract：

In response to the problem of large-scale duck farming in winter， temperature， relative humidity， and wind speed tests were conducted on the fermentation bed net duck breeding house. A CFD model of the duck house was constructed， and environmental simulations of temperature， humidity， and airflow fields were carried out. Optimization plans for the structure and equipment of the duck house were proposed. The results showed that during the testing period， the temperature inside the room was 9.80℃-17.68℃， and the average temperature difference between inside and outside the room was 8.29℃. The average relative humidity inside the house was 81.90%， slightly lower than outside. The wind speed inside the building was 0.04-0.57m/s. The temperature， humidity， and airflow fields of the original model were simulated and compared with the measured results， with an average relative error of 2.96% to 7.67%. Further， the opening angle of the ventilation was increased from 45° to 50°， and two warm air fans were added to the initial end and center of the fermentation bed located on the wet curtain side to optimize the environment of the original duck house. Simulation showed that the temperature inside the optimized model was increased from 15.32℃ to 16.85℃， reaching the optimal temperature for the growth of egg ducks during the breeding period. The average humidity inside the house was decreased to 66.52%， a decrease of 11.07% compared to the original， which was lower than the threshold of high humidity environment. The average wind speed inside the building was increased from 0.29m/s to 0.38m/s， improving the ability to remove harmful gases. This study can provide a theoretical basis for improving the environmental quality of similar duck houses.

Keywords：

duck coop; fermentation bed net; CFD; environmental optimization

0 引言

设施环境显著影响鸭的生长性能，近年来，鸭养殖逐渐向标准化、集约化、规模化和福利化方向发展，传统的临水养殖和旱地平养由于设施简陋、饲养环境差逐渐被淘汰［1］。发酵床网养鸭舍利用发酵床实现鸭粪原位发酵解决鸭粪处理难的问题，该模式下鸭在网床上自由活动，摆脱鸭笼的束缚，有效地将鸭群与鸭粪分开，减少鸭群病菌的感染率，此外粪便通过网孔下落会进行原位发酵，将粪污转化成有机肥料，可用于农业生产［2］。但由于发酵床网养鸭舍集约化程度高，会造成舍内气流不通畅，加之冬季舍外温度低，为了达到保温的效果鸭舍门和通风窗多为关闭状态，会导致舍内气流蓄积有害气体浓度变高［3］，进而影响鸭的健康［4］，降低鸭的生产效益。因此，对冬季发酵床网养鸭舍环境评估与优化等问题亟待解决。

目前，畜禽舍内环境测试优化的重要手段是计算流体力学（Computational Fluid Dynamic，CFD）技术，CFD技术是通过数据分析模拟为试验人员提供详尽的气动数据和图像［5］。CFD技术在国外起步比较早，并在近年来有了较大的突破。1990年，Choi等［6］最先使用k-ε湍流模型对畜禽舍内环境进行了初步模拟，模拟的过程相对较长，由于自然条件的欠缺和技术的不成熟，所得的气流模拟结果与实际所测误差较大。1992年，Hoff等［7］将畜禽舍模型假定成一个能够发热的板状结构，模拟了畜禽舍内的热交换现象，结果表明模拟结果具有很好的准确性。此后，大部分研究都会将畜禽舍气流的模拟与温度的模拟结合起来研究。Kic等［8］利用CFD技术对夏冬两季肉鸡舍内的通风状况进行二维和三维模拟，试验证明相比夏冬两季肉鸡舍二维通风模型，三维模拟数据精确度更加精确。Sun等［9］又利用CFD技术模拟了猪舍的气流和有害气体，构建猪舍的气流场三维模型，结果发现模拟值与实测值无明显差异，从而表明CFD模型是评估畜禽舍气流场重要手段。国内针对CFD技术对畜禽内环境评估起步比较晚，多用于鸡舍、牛舍等温度场和气流场的研究。李文良等［10］利用CFD技术模拟了冬季密闭式鸡舍纵向通风过程，对鸡舍内进风口大门的安装位置、高度进行模拟，探究其对鸡舍气流的影响。邓书辉等［11］运用CFD对低屋面横向通风牛舍气流场进行模拟，旨在为低屋面横向通风牛舍的结构提供合理的优化。近年来，国内畜禽调控技术发展快速，畜禽生长性能等一系列问题得到了有效的解决，显著提高养殖水平。

本文拟通过现场实测结合CFD模拟的方式，对发酵床网养鸭舍内温湿度、风速进行测试与分析，并对鸭舍内环境进行CFD模拟，针对模拟的温度场、湿度场和气流场存在的问题进行设施结构与装备优化。

1 材料与方法

1.1 供试鸭舍

1.1.1 鸭舍结构

试验鸭舍位于江苏省扬州市高邮鸭集团养殖示范基地（119°46′E、32°78′N）。鸭舍平面结构与尺寸如图1（a）所示，鸭舍长100m，宽15m，檐高3.4m，屋脊高5.3m，呈南北走向。南北两侧纵墙上各安装有24个通风小窗和12扇窗户，通风小窗尺寸为0.6m×0.3m，通风小窗导流板开启角度为0°～100°，窗户尺寸为长1.95m、宽0.97m、距地面1.5m，冬季正常情况下窗户全关闭，湿帘安装于西侧墙头，尺寸为12.4m×2m，厚度为0.2m，距离地面1m。污道一侧山墙并列安装6台风机，从左至右依次标号为1、2、3、4、5、6，每台间距1.03m，风机长宽都是1.3m，距离地面1.5m，理论风量为31 800m3/h［12］，如图1（b）所示。

网床长90m，宽15m，厚0.1m，距离地面1.9m。网床分为活动区和休息区，活动区11.5m，休息区3.5m。鸭舍内配套自动饮水和喂料系统以及粪便自动翻耙系统。自动饮水系统采用乳头式位于发酵床两侧，自动喂料系统位于发酵床中间，白天每隔3h自动喂料一次，鸭群在网床上层自由活动，粪便通过漏缝板下落。

1.1.2 通风模式

舍内饲养苏邮1号蛋鸭约4000羽。采用负压纵向通风模式，测试期间舍内蛋鸭处于250日龄，属于育成期蛋鸭，理想温度范围是16℃～18℃［13］，在满足最小通风量的前提下结合当地气候与多年养殖经验进行调整。根据江苏地区冬季气候特点，设定试验鸭舍冬季理想温度为16.0℃。（1）当舍内温度≤16.0℃时，执行最小通风量的风机运行状态，开启3号风机，执行顺序为：开40s，停50s，再开40s，再停50s，如此循环往复。（2）当16.0℃＜舍内温度＜18.0℃时，开启3号和5号两台风机，执行顺序为：开50s，停60s，再开50s，再停60s，如此循环往复。舍内温度每升高0.5℃，该2台风机多运行25s，如舍内气温为16.5℃，则两台风机开75s，停60s，如此循环往复。（3）当舍内温度=18.0℃时，则开启1、3、5号三台风机，执行顺序为：开40s，停50s，如此循环往复。（4）当舍内温度＞18.0℃，1、3、5号三台风机连续运行，并增开4号风机，直到降至18.0℃，关闭4号风机。

1.2 试验设计

本试验主要测量舍内不同区域的温度、相对湿度、风速以及舍外温湿度。舍内和舍外共布置8个测点，舍外放置在鸭舍门前，测量舍外的温度和相对湿度，舍内网床四等分点设置5个测点，中间两侧各放置1个。

温度和相对湿度采用HOBO温湿度记录仪UX100-011实时监测记录（温度量程-20℃～70℃，精度为±0.2℃；相对湿度量程1%～100%，精度为±2.5%），记录时间为30min一次；风速采用Testo425热敏风速仪（测量范围0～20m/s，分辨率为0.01m/s）测定，测定方式为手持读数，测定时间为每天8：00、10：00、12：00、14：00、16：00、18：00、20：00。

2 结果与分析

2.1 舍内外温湿度分析

图2和图3显示连续8d（2022-01-18 0：00—2022-01-26 0：00）舍内外温度和相对湿度变化情况，测试期间舍内纵向温湿度采集点共有5处，选取5处取平均值作为舍内整体温湿度。

由图2可知，连续8d舍外温度范围为2.22℃～11.93℃，平均温度为5.77℃；舍内温度范围为9.80℃～17.68℃，平均温度为14.05℃。舍内外平均温差8.29℃，最大温差12.39℃，说明供试鸭舍冬季具有保温效果，但相比育成期蛋鸭最适生长温度16℃～18℃，舍内温度偏低。

由图3可知，连续8d舍外相对湿度最高为95.55%，最低为59.65%，平均相对湿度82.57%；舍内相对湿度最高为90.53%，最低为67.57%，平均相对湿度81.90%。供试鸭舍内相对湿度曲线变化趋势与舍外基本一致，日变化幅度均小于舍外，变化幅度小减少鸭的应激反应，更有利于鸭的生长。畜禽养殖行业将相对湿度>75%定义为高湿环境，舍内大部分时间内处于高湿环境。

2.2 舍内风速分析

图4显示不同时刻舍内不同区域风速变化。由图4可知，舍内风速变化范围为0.04～0.57m/s。西侧（湿帘侧）显著高于中部和东侧（风机侧）（P＜0.05），东、西、中三个区域风速都是先向上升后下降再上升，在正午时刻达到顶峰，并且由湿帘至风机呈上升趋势；整体来看，舍内风速远低于适合家禽生长的最佳风速范围1.5～2.0m/s［14］，风速低会造成舍内气流不通畅，有害气体聚集。

3 CFD模拟

3.1 三维模型建立

假设鸭舍墙壁和屋顶为绝热壁面，同一侧墙体的不同区域温度看作相同；地面定义为恒温壁面；进风口气流匀速等温，流体的流动遵循四大定律包括；舍内各维护结构定义为热流量壁面，忽略物线水线等结构对舍内气流的影响；由于蛋鸭在活动区较分散，模型建立忽略蛋鸭对舍内气流的影响。

根据实地测量和模拟工况的数据，在上述内部结构与装备简化的基础上利用SolidWorks对几何进行建模后，导入Ansys Design model进行模型布尔操作，同时对各边界条件进行命名以用于后期数值求解，Z轴正向为东向，X轴正向为北向，建立如图5所示的三维鸭舍模型。

3.2 网格划分

将建立的鸭舍三维模型导入ANSYS fluent中，采用fluent meshing对三维几何模型进行离散。划分网格类型为多面体网格，对45°导向通风小窗位置进行网格细化，定义该区域局部面网格尺寸为30mm，风机和湿帘面网格尺寸定义为150mm，几何模型整体最大面网格尺寸定义为800mm，网格增长率为1.2，面网格生成采用曲率与近似的方法以进一步增强网格质量，最终几何模型生成网格数量为1381204，网格质量评判指标最大畸变度和最小正交质量分别为0.86和0.14，建立如图6所示的鸭舍网格划分图。

3.3 边界条件及求解器设置

划分网格之后，将网格文件导入求解器中，进行边界条件的设置，具体模拟工况边界条件设置为：（1）假定舍内热空气为连续、不可压定的理想气流。（2）假定鸭舍屋顶为绝热壁面，地面定义为恒温壁面，其数值为模拟时间内室内平均温度（16.02℃）；墙体定义为热流量壁面，以实测数据为准（东墙19.9W/m2；西墙10.1W/m2；南墙21.5W/m2；北墙9.9W/m2）。（3）忽略舍内供水供料等管道对舍内气流的影响；舍内封闭大窗小窗、门区域定义为与外界气温对流的壁面，对流换热系数7.5W/m2，具有45°通风导向小窗定义为自由进出口。（4）风机开1台（中间第3号风机），其他风机定义为绝热壁面；湿帘为均匀进风入口，冬季假定湿帘水泵未开启不需要降温，其温湿度及气流流速定义为测试值。（5）发酵床定义为多孔介质以减少仿真计算成本。

根据Forchheimer描述空气流经多孔介质定律［15］

式中：

1/α——黏性阻力系数，8341605.404m-2；

β——惯性阻力系数，120.5191837m-1；

ΔP——流体在介质中的压力损失，N；

d——多孔介质长度，m；

v——介质中空气流速，m/s；

ρ——介质的渗透率；

μ—孔黏度。

代入式（1）得假定多孔介质后压降

式中：

μ0——空气黏度，μ0=1.7894×10-5为15℃时的黏度；

T——温度，℃；

B——与气体种类有关的常数，空气的B=110.4。

发酵床估计孔隙率为20×25×13×10/300×300=0.72，由此可得发酵床网孔的孔隙率为0.72。

模拟过程中舍内处于稳定状态，风机正常运行，舍内入风口与出风口气流速率恒定不变，所以使用基于压力稳态（状态不随时间的变化而改变）求解数学模型；数值求解采用标准k-ε湍流模型，标准k-ε湍流模型具有计算速度快、对硬件条件低的优点，是近年针对畜禽舍数值模拟的最广泛完整是湍流模型，因此，本文所有的模拟计算均采用了标准的k-ε模型作为湍流模型。

3.4 模拟结果验证

因为2022年1月18日12：00这一时刻舍内外温差最大，所以选择这一时刻进行模拟。对应模拟的温度场、湿度场和气流场如图7所示。

由模拟结果可以看出供试鸭舍冬季内环境存在以下问题：（1）舍内整体温度偏低（15.32℃），低于育成期蛋鸭最适生长温度（16℃～18℃），由于养殖人员进出鸭舍，出口门经常打开，导致靠近湿帘侧温度较低（12.53℃）。（2）舍内整体湿度较高（77.32%），高于高湿环境阈值75%。（3）舍内整体风速较低（0.29m/s），特别是靠近风机一侧风速低于0.2m/s。

为验证CFD仿真模拟的准确性，详细对比分析实测值与模拟值，引入相对误差计算公式［16］，如式（3）所示。

式中：

Ev——实测值与模拟值之间的相对误差；

Cs——模拟值；

Cm——实测值。

评估模拟结果的准确性就是比较CFD仿真模拟中与测点同一位置、同一时刻的数据，本文选择鸭舍同一时刻5个测点（距离湿帘1m、距离湿帘25m、距离湿帘50m、距离湿帘75m、距离湿帘98m靠近风机）的温度、相对湿度和气流对比分析得出结果如图8所示：温度相对误差范围1.40%～6.09%，平均误差3.78%；湿度的相对误差范围1.71%～4.54%，平均误差2.96%；气流相对误差范围4.48%～10.00%。平均误差7.67%，说明模拟结果与实测数据吻合性较高，模拟值较准确反映试验鸭舍内温度、相对湿度和气流速度变化情况，并且以该模拟为基础可以对该鸭舍环境进行有效的评估与分析。

3.5 舍内温度、湿度和气流优化

为解决模拟结果出现的温度低、湿度大、风速低等问题，利用CFD数值仿真技术对鸭舍结构与装备进行优化设计，提出在改进通风小窗开启角度的基础上增加暖风机的优化方案，具体方案与优化结果如下。

暖风机边界条件设置：鸭舍总面积1500m2，预期设定温度从14℃提升至16℃所需热量按照式（4）计算。

Q=C×m×Δt（4）

式中：

Q——需要总热量，J；

C——比热容，J/（kg·℃）；

m——质量，kg；

Δt——需要提升温度的差值，℃。

鸭舍总面积1500m2，空气比热容1.004kJ/（kg·℃），空气密度为1.23kg/m3，求得鸭舍温度从14℃提升至16℃所需热量为9261.9kJ，由于鸭舍存在一些保温结构，假定鸭舍与外界进行热量交换损失50%［17］，所以需要总热量18523.8kJ，根据电力热量换算系数860［18］，所以需要20kW功率的工业柴油暖风机。出风口定义为速度出口，出风温度为45℃，进风口定义为速度入口，风速与出风口相同。为避免暖风机热风直接吹向鸭对鸭产生应激反应，仿真中将暖风机风口仰角调控至30°，鸭舍内布置两台，位于靠近湿帘侧的发酵床初始端和发酵床中心位置。

在其他原有模拟边界条件和求解设置不变的情况下，对2022年1月18日12：00这一时刻改进后的鸭舍几何模型经过仿真模拟得出优化后的仿真结果如图9所示。

由优化后温度场、湿度场和气流场可知：（1）温度场：增加通风小窗开启角度会加大舍内的能量损失，但通过在靠近湿帘侧的发酵床初始端和发酵床中心位置添加两台暖风机使得舍内各位置温度增高，并且达到育成期蛋鸭生长的最适温度。（2）湿度场：增添暖风机使得舍内温度增高，湿度降低，低于高湿环境阈值75%，特别是靠近风机侧温度明显上升，湿度明显降低。（3）气流场：优化后提高了舍内风速，并且由于暖风机在发酵床上方因此对下方气流分布影响忽略不计，优化后明显提高了发酵床上方鸭活动区域气流流速并且暖风机将原有的部分盛行气流导引至发酵床上方，有利于鸭活动范围内CO2和NH3的稀释。

4 讨论

Yahav等［14］研究发现，畜禽舍风速范围为1.5～2.0m/s时最适合禽类生长，本研究测试期间舍内风速范围0.04～0.57m/s，远低于禽类生长的最适风速范围，这也导致舍内CO2浓度远高于定期清粪鸭舍CO2浓度阈值1500mg/m3，因此还需要优化舍内通风模式，以减少因CO2浓度过高对鸭生长的影响，此外，由于前期试验准备不充分，没有测量舍内颗粒物浓度，颗粒物含量过高抑制鸭的生长，保证鸭舍颗粒物含量在适当范围也是改善鸭舍环境的主要措施之一，在后续的研究中，会积极与企业沟通交流，对舍内颗粒物浓度进行测量改善。

对鸭舍环境模拟发现，优化后模型增加了舍内风速特别是网床上端鸭活动区域，暖风机设置向上仰角30°，一方面降低暖风直吹造成鸭热应激，另一方面可以将鸭活动区域CO2和NH3排至网床上方从通风小窗排出，但具体排出量还未计算，这将是下一步的研究重点。另外，模拟值与测试值的个别误差比较大可能是由于建模过程中未考虑蛋鸭对舍内环境因子的影响所造成的。同时，本文对装备与结构优化的方案未进行实测验证，后续将开展相关研究。

5 结论

本文通过实地测量发酵床网养鸭舍内冬季温湿度、风速，构建了鸭舍的CFD模型，验证了模型的准确性，对舍内存在的温度低、湿度大和风速低等问题进行优化，提出改变通风小窗开启角度（从45°增大到50°）和增加暖风机（在位于靠近湿帘侧的发酵床初始端和发酵床中心位置增设两台暖风机）的方案。

1）测试期间，舍内温度9.80℃～17.68℃，舍外温度2.22℃～11.93℃，舍内外平均温差8.29℃，最大温差12.39℃，说明供试鸭舍冬季具有保温效果，但相比育成期蛋鸭最适宜的生长温度16℃～18℃，舍内温度偏低；舍内相对湿度67.57%～90.53%，舍外相对湿度59.65%～95.55%，舍内相对湿度日变化幅度均小于舍外，但大部分时间处于高湿环境。

2）根据现场实测环境因子和维护结构尺寸，对冬季发酵床网养模式鸭舍进行CFD仿真模拟，并将现场实测温度、相对湿度和风速与CFD模拟数值进行验证结果发现具有较好的吻合性，说明该模型可有效评估冬季发酵床网养鸭舍内环境，并可对鸭舍结构与装备进行优化模拟。

3）优化后结果表明：舍内温度从15.32℃提高到16.85℃，达到育成期蛋鸭生长最适温度，靠近湿帘侧温度达到14.87℃；舍内平均湿度降至66.52%，较原来下降11.07%；舍内平均风速从0.29m/s提高到0.38m/s，并且在优化过程中我们将暖风机出风角度风口仰角调控至30°，一方面降低因吹风过大对鸭的应激，另一方面，暖风机明显提高发酵床上层鸭活动区域气流速率，利于CO2和NH3通过通风小窗排出舍外，在有效提高发酵床上层鸭活动区域气流流速的同时，加速舍内CO2和NH3等有害气体排出舍外。

参 考 文 献

［1］ 侯水生， 刘灵芝. 2021年水禽产业现状、未来发展趋势与建议［J］. 中国畜牧杂志， 2022， 58（3）： 227-231， 238.

［2］ 李明阳， 应诗家， 戴子淳， 等. 新型肉鸭养殖模式生产性能及经济效益对比分析［J］. 中国家禽， 2020， 42（4）： 80-85.

Li Mingyang， Ying Shijia， Dai Zichun， et al. Comparative analysis of production performance and economic benefits of new meat duck production systems in China ［J］. China Poultry， 2020， 42（4）： 80-85.

［3］ Zhao Y， Zhao D， Ma H， et al. Environmental assessment of three egg production systems-Part Ⅲ： Airborne bacteria concentrations and emissions ［J］. Poultry Science， 2016， 95（7）： 1473-1481.

［4］ 黄炎坤， 刘健， 范佳英， 等. 鸡舍内环境季节性变化与鸡群死淘率的关联分析［J］. 中国畜牧杂志， 2012， 48（6）： 57-60， 64.

［5］ 王福军. 计算流体动力学分析： CFD软件原理与应用［M］. 北京： 清华大学出版社， 2004.

［6］ Choi H L， Albright L D， Timmons M B， et al. An application of the k-epsilon turbulence model to predict air distribution in a slot-ventilated enclosure ［J］. Transactions of the ASAE， 1988， 31（6）： 1804-1814.

［7］ Hoff S J， Janni K A ， Jacobson L D .Three-dimensional buoyant turbulent flows in a scaled model， slot-ventilated， livestock confinement facility ［J］. Transactions of the ASAE， 1992， 35（2）： 671-686.

［8］ Kic P， Zajicek M. A numerical CFD method for the broiler house ventilation analysis ［J］. Journal of Information Technology in Agriculture， 2011， 4（1）： 7.

［9］ Sun H， Keener H M， Deng W， et al. Development and validation of 3-D CFD models to simulate airflow and ammonia distribution in a high-riseTM hog building during summer and winter conditions ［J］. Agricultural Engineering International： CIGR Journal， 2004.

［10］ 李文良， 施正香， 王朝元. 密闭式平养鸡舍纵向通风的数值模拟［J］. 中国农业大学学报， 2007， 12（6）： 80-84.

Li Wenliang， Shi Zhengxiang， Wang Chaoyuan. Numerical simulation of tunnel ventilation system of plane-raising enclosed henhouse ［J］. Journal of China Agricultural University， 2007， 12（6）： 80-84.

［11］ 邓书辉， 施正香， 李保明， 等. 低屋面横向通风牛舍空气流场CFD模拟［J］. 农业工程学报， 2014， 30（6）： 139-146.

Deng Shuhui， Shi Zhengxiang， Li Baoming， et al. CFD simulation of airflow distribution in low profile cross ventilated dairy cattle barn ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2014， 30（6）： 139-146.

［12］ 丁露雨， 鄂雷， 李奇峰， 等. 畜舍自然通风理论分析与通风量估算［J］. 农业工程学报， 2020， 36（15）： 189-201.

Ding Luyu， E Lei， Li Qifeng， et al. Mechanism analysis and airflow rate estimation of natural ventilation in livestock buildings ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2020， 36（15）： 189-201.

［13］ Hutchings N J， Sommer S G， Andersen J M， et al. A detailed ammonia emission inventory for Denmark ［J］. Atmospheric Environment， 2001， 35（11）： 1959-1968.

［14］ Yahav S， Straschnow A， Vax E， et al. Air velocity alters broiler performance under harsh environmental conditions ［J］. Poultry Science， 2001， 80（6）： 724-726.

［15］ Skrivan M， Marounek M， Bubancova I， et al. Influence of limestone particle size on performance and egg quality in laying hens aged 24-36 weeks and 56-68 weeks ［J］. Animal Feed Science and Technology， 2010， 158（1-2）： 110-114.

［16］ 马瑞钰， 詹凯， 司绍宏， 等. 规模化蛋鸡场不同养殖模式冬季能耗差异性分析［C］. 中国畜牧业协会. 第七届（2015）中国蛋鸡行业发展大会会刊［A］. 中禽传媒（江苏）有限公司， 2015： 166-170.

［17］ 程琼仪. 叠层笼养蛋鸡舍夏季通风气流CFD模拟与优化［D］. 北京： 中国农业大学， 2018.

Cheng Qiongyi. The CFD simulation and improvements of ventilation in stacked-cage laying hen housing ［D］. Beijing： China Agricultural University， 2018.

［18］ 蔡文庆， 周仁才， 王亚男， 等. 风机盘管机组冷热量换算方程［J］. 集美大学学报（自然科学版）， 1990（2）： 35-41.

Cai Wenqing， Zhou Rencai， Wang Yanan， et al. Equation for calculating cooling and heating capacities of fan coil unit ［J］. Journal of Jimei University （Natural Science）， 1990（2）： 35-41.