基于simulink的猪粪好氧堆肥热量平衡仿真模型研究

基于SIMULINK的猪粪好氧堆肥热量平衡仿真研究王永江黄光群韩鲁佳（中国农业大学工学院，北京 100083）摘要：本研究在热量守恒定律的基础上，利用SIMULINK工具箱建立猪粪好氧堆肥过程热量平衡仿真模型，更加直观的对反应器好氧堆肥热量平衡的各个组成部分进行仿真分析。

热量平衡仿真模型考虑了进出口气体携带的热量、有机质生物降解产热量和堆肥反应器表面传导热量。

利用实验室堆肥反应器进行仿真模型验证试验，结果表明：仿真模型模拟值与堆肥试验实际测量值表现出较好的一致性，在通气量较小的情况下，反应器表面的热传递是主要的热量交换方式；变换试验参数可以进行热量动态仿真分析。

所建好氧堆肥过程热量平衡模型有一定的伦理价值和实际意义。

关键词：仿真模型，好氧堆肥，热量平衡，反应器堆肥中途分类号：S141.4文献标识码：A0 引言随着我国畜禽养殖业得集约化发展，畜禽粪便年产量约为30亿吨[1]。

高温好氧堆肥技术是畜禽粪便无害化，资源化利用的有效途径[2~4]。

好氧堆肥过程是一个有堆肥微生物参与的，复杂的物理和化学变化的过程[2,5]。

表现为有机质成分和含量的变化、含水率的变化、温度的变化、气体成分（如氧气、二氧化碳、甲烷、氨气和硫化氢等）的变化及微生物的生长和消亡等现象。

堆肥数学模型提供了一种描述堆肥系统中物理学、化学和生物学变量相互制约的关系的方法[3]。

过程模拟技术正越来越多的用来预测堆肥结果和指导工程实践[5~7]。

MATLAB经过长期不断的发展和完善，已经成为当今世界上最优秀的数值计算软件，其自带的SIMULINK工具箱由模块库、模型构造、指令分析和演示程序构成，是模块化和模型化的动态系统仿真环境。

本研究在建立数学模型的基础上，利用SIMULINK 工具箱实现猪粪好氧堆肥工程中的热量平衡仿真，使得对堆肥过程热量平衡的研究更收稿日期：修订日期：项目基金：“十一五”国家科技支撑计划资助项目（2006BAD10B05）；中国农业大学研究生科研创新专项资助项目。

作者简介：王永江（1985—），男（汉），山东寿光人，博士研究生，主要从事生物质资源与利用研究。

北京海淀区清华东路17号232信箱，100083。

Email:****************通讯作者：韩鲁佳（1964—），女（汉），籍贯，教授，博士，主要从事生物质资源与利用研究。

北京海淀区清华东路17号232信箱，100083。

Email:*************.cn 加便捷和直观，以期进一步探究堆肥过程动力学特征和推动堆肥过程实时控制的实现。

１好氧堆肥热量平衡模型根据能量守恒定律，堆肥系的热量平衡可以用方程描述为：热量积累=热量输入+热量产出-热量输出[8]，上述方程是堆肥过程热量平衡模拟和仿真的理论基础。

在能量守恒原理基础上建立的实验室反应器堆肥系统的热量平衡模型[8,9]如下：()i i o o CdT dVSmc G H G H UA T T H dt dt=---+(1)其中：G表示进气口或出气口气体质量流量，/gk d；T表示实验室环境温度；H 表示进气口或出气口气体的焓值；U表示堆肥反应器的总传热系数，W/㎡·℃；A表示反应器的总传热面积，2m，取UA=0.4077W/℃[9]；CH表示生化反应热，KJ/k g；VS表示堆肥过程中VS的含量，gk；t表示反应时间，d。

m表示堆体质量，gk；c表示堆体和反应器的比热，KJ/k g·℃。

假定堆肥过程中有机质的降解符合Monod一级降解动力学方程形式：()dVSk VSdt=-⋅(2) 其中：k表示反应速率，1d-。

反应速率受温度、含水率、氧浓度和自由空域等因素的影响[3,8]：2T M O Fask k k k k=⋅⋅⋅(3)其中：Ti k 表示温度对反应速率的影响；M k 表示含水率对反应速率的影响校正，取值范围：0~1；2O k 表示氧浓度对反应速率的影响校正，取值范围：0~1；Fas k 表示自由空域对反应速率的影响校正，取值范围：0~1。

温度对反应速率的影响可用Andrews-Kambhu 方程[10]表达如下：(30)(60)012[]T T T k k c c --=⋅-(4)其中：0k 表示不同反应阶段的反应速率常数；1c ，2c 表示温度系数； T 表示堆体温度，℃。

含水率、氧浓度和自由空域对反应速率的影响可用方程[2]分别表达如下： (17.6847.0622)11M M k e-⋅+=+ (5)2222O O O V k V =+(6)(23.675 3.4945)11Fas Fas k e -⋅+=+(7)其中：M 表示堆体含水率，%；2O V 表示堆体内部氧气体积浓度，%；Fas 表示堆体的自由空域，%。

2 SIMULINK 建模2.1 输入空气所携带的热量泵入堆肥反应器的空气所携带的热量为空气质量与输入空气所携带焓值的乘积。

空气的焓值是每千克空气及其所携带的水蒸气焓值之和[11]，可表达为： (1.01 1.88)2491H x t x =++(8) 其中：H 表示出口或入口气体的焓值，kJ/kg ；x 表示出口或入口气体的湿含量，又称绝对湿度，表示每千克干空气中含有水蒸气的质量，kg/kg ；t 表示出口或入口气体的温度，℃。

输入空气所携带热量的SIMULINK 模型如图1所示：图1输入气体热量的模型框图 Fig.1 Thermal model diagram of input gas2.2 出口气体所携带的热量根据堆肥过程气体流动过程分析，环境空气在通风泵作用下进入反应器参与生化反应，反应后可认为其总气体体积近似不变，所以进出口气体体积流量保持不变。

气体在流经高温高湿堆体之后，温度基本上与堆体温度相同，并且基本达到饱和温度，可近似认为出口气体的相对湿度为90%。

出口气体所携热量的SIMULINK 模型如图2所示。

图2 输出气体热量的模型框图Fig.2 Thermal model diagram of output gas2.3 有机质降解反应热好氧堆肥过程中能够降解的有机物主要有糖类、蛋白类、脂肪类和纤维类，不同的有机质在堆肥过程中的降解呈现梯次性[12]，由于不同的降解阶段主要参与降解的有机质的种类不同，有机质降解热系数亦应该是时间的函数，通过堆肥过程取样和文献数据[3,8,12]，得到堆肥过程中降解反应热系数的时间函数如下：7771.837100 3.5(0.0409 1.9802)10 3.591.61210914CtH t tt⎧⨯<<⎪=-+⨯≤≤⎨⎪⨯<<⎩(9)有机质降解反应热的SIMULINK模型如图3所示。

图3 有机质降解反应热的模型框图Fig.3 Thermal model diagram of organics degradation2.4 反应器表面热损量堆肥反应器表面的热损量用总传热系数、总传热面积和温差以乘积UA（T-T0）的形式表达。

取U值为0.4077W/℃[9]。

反应器表面热损量的SIMLINK模型如图4所示。

图4 反应器表面热损失模型框图Fig.4 Thermal model diagram of heat loss from reactor surface3 结果和讨论本研究利用实验室小型堆肥反应器进行猪粪、麦秸混合好氧堆肥试验研究。

堆肥反应器[13]有效容积15L，通风方式为连续通风（升温期/高温期/降温期的通风速率分别为0.4/0.3/0.2L/min），堆肥试验和数据采集周期为14d。

温度变化的实际测量值与模型仿真值的验证结果如图5所示。

温度变化的实际测量值与模型仿真值有一定的差距，这可能是因为在建立仿真模型时忽略了辐射散热的影响，或是因为堆体空间梯度性较大，与模型对堆体均匀性的假设不符，影响了模型的准确性。

图5 温度变化的实测值与模拟值对比Fig. 5 Predicted and observed values of Temperature 图6所示为进口输入空气所携带的热量，主要由通风量的大小决定，图7显示了出口气体所携带的热量，图8所示为通过反应器表面所散失的热量，图9显示了有机质降解的产热量，通过比对图6~9可知，有机质的降解是主要的热量来源形式，而热量的散失主要是通过反应器表面的热交换，因此，建造堆肥反应器是要考虑创造有利条件保障有机质的降解和注意保温设计。

有机质的有效降解与pH、含水率、调理剂、C/N 等诸多因素有关，不作进一步讨论，但热量的散失方式主要取决于堆肥反应器的型式和通风方式[14]，出口气体所携带的热量约占总热量的17~55%[15]，本仿真结果中出口气体所携带的热量占总热量的比例较小，可能是因为模型没有考虑到水蒸气潜热的影响，也可能与通气量小和气流平缓有关。

图6 进口气体的焓值Fig.6 Enthalpy value of inlet air图7 出口气体的焓值Fig. 7 Enthalpy value of exhaust air图8 反应器表面所散失的热量Fig. 8 Thermal loss through conduction duringcomposting图9 BVS降解的产热量Fig. 9 Thermal generation from BVS degradation4 结束语随着好氧堆肥技术研究的深入，对堆肥过程机理的研究提出了更高的要求。

采用SIMULINK工具箱可以简单快捷的建立好氧堆肥过程动力学仿真模型，直观的分析各个状态变量的变化，得到仿真结果，是进行好氧堆肥过程仿真分析的有力工具。

利用小型堆肥反应器对所建好氧堆肥仿真模型进行了验证试验，结果表明仿真结果与实际测量结果一致性良好，所建仿真模型对深入研究堆肥机理，指导工程实践有一定的理论价值和实际意义。

[参考文献][1]王晓燕, 黄光群, 韩鲁佳. 鸡粪工厂化堆肥过程中有机质含量预测模型[J].农业机械学报, 2010, 41(1): 101~105.Wang X Y, Hang G Q, Han L J. Prediction Model for Organic Matter Content in Chicken Manure during Plant-field Composting [J]. Transactions of the Chinese Society for Agricultural Machinery, 2010, 41(1): 101~105. (in Chinese with English abstract)[2]蔡建成. 堆肥工程与堆肥工厂[M].北京: 机械工业出版社, 1990.[3]Hang R T. The practical handbook of composting [M]. Boca Raton: Lewis Publishers, 1993.[4]李季, 彭生平. 堆肥工程实用手册[M]. 北京: 化学工业出版社, 2005.[5] Luong J H T, V olesky B, Heat evolution during microbial process – estimation, measurement,and applications [J]. Advances in Biochemical Engineering/Biotechnology, 1983, 28: 1~40.[6] Yeomans J S, Huang G H, Yoogalingam R. Combining simulation with evolutionaryalgorithms for optimal planning under uncertainty: an application to municipal solid waste management planning in the regional municipality of Hamilton-Wentworth [J]. Journal of Environmental Informatics, 2003, 2: 11~30.[7] Huang B, Xiong D, Li H. An integrated approach to real-time environmental simulation andvisualization [J]. Journal of Environmental Informatics, 2004, 3: 42~50.[8] Mason I G. Mathematical modeling of the composting process: A review [J].WasterManagement , 2006, 26: 3~21.[9] 吕黄珍. 猪粪麦秸好氧堆肥工艺参数优化及过程模拟[D].北京: 中国农业大学, 2007.LV H Z. Process Optimization and Simulation on Pig Slurry-Wheat Straw Composting [D].Beijing: China Agriculture University, 2007. (in Chinese with English abstract)[10] Andrews J F, Kambhu, K. Thermoplastic aerobic digestion of organic solid wastes [J]. U.S.EPA NTIS, 1973: 222~396.[11]潘永康. 现代干燥技术[M].北京: 化学工业出版社, 2002.[12]Kaiser J. Modeling composting as a microbial ecosystem: a simulation approach [J].Ecological Modeling, 1996, 91: 25~37.[13] 张锐. 小型试验室好氧堆肥反应器系统研究[D]. 北京: 中国农业大学, 2006.Zhang R. Study on the Bench-Scale Aerobic Composting Reactor System [D]. Beijing:China Agriculture University, 2006. (in Chinese with English abstract)[14] Petric I, Selimbasic V . Development and validation of mathematical model for aerobiccomposting process [J]. Chemical Engineering Journal. 2008, (139): 304~317.[15] Ahn H K, Richard T L, Choi H L. Mass and thermal balance during composting of poultrymanure-Wood shavings mixture at different aeration rates [J]. Process Biochemistry. 2007,(42): 215~223.Research on thermal balance models during swine manure composting based on SIMULINKWang Yongjiang Huang Guangqun Han Lujia(College of Engineering, China Agricultural University, Beijing 100083,China)Abstract: In this study, simulation models regarding composting heat were developed based on the law of thermal conservation during composting and built by SIMULINK tools. Forced convention and conduction were considered as the main heat exchanging form. An aerobic composting process of swine manure was conducted in order to verify the effectiveness of the developed models. Modeling solutions were consistent with experimental results with an acceptable accuracy level through verification-based composting experiment. Sensitivity analyses of the models showed that conduction was the main heat removal channel under low ventilation.Key words: Simulink, Aerobic composting, Thermal balance, Reactor composting。