翟耀芳等

摘要：隨著養(yǎng)殖業(yè)向高度集約化方向發(fā)展，雞舍內(nèi)微氣候環(huán)境對雞的生產(chǎn)、生活性能有日漸重要的影響。本研究利用多孔介質(zhì)模型將單個雞籠及其雞體簡化為一個具有相同阻力特性的長方體，采用數(shù)值模擬方法研究了密閉式籠養(yǎng)蛋雞舍內(nèi)的氣流組織及熱濕環(huán)境。結(jié)果表明：在進口端氣流不穩(wěn)定，存在渦流，經(jīng)過一段距離后達到穩(wěn)定，氣流間相互影響變??；在雞籠范圍內(nèi)，溫度和濕度沿高度和寬度方向分布較為均勻，溫度沿縱向逐漸上升，相對濕度則逐漸下降。因此，為了使縱向通風能更好發(fā)揮作用，必須采取措施改善其氣流組織。

關(guān)鍵詞：縱向通風；籠養(yǎng)蛋雞舍；計算流體力學；氣流組織；熱濕環(huán)境

中圖分類號： TU264+.7 文獻標志碼： A 文章編號：1002-1302（2014）03-0337-06

雞舍內(nèi)微氣候環(huán)境包括氣流組織、溫度、相對濕度等，對雞的生產(chǎn)、生活性能有重要影響。由于養(yǎng)雞業(yè)向高密度、集約化發(fā)展，目前我國多數(shù)集約化蛋雞養(yǎng)殖場存在舍內(nèi)微環(huán)境不穩(wěn)定、空氣污濁、粉塵較大、缺乏現(xiàn)代化設(shè)施、環(huán)境控制能力差、存在疾病隱患等問題[1]。通風是密閉式雞舍內(nèi)環(huán)境調(diào)控的重要手段，良好的通風對改善舍內(nèi)熱環(huán)境條件和空氣質(zhì)量，提高雞的健康和生產(chǎn)水平有重要作用。與傳統(tǒng)橫向通風方式相比，縱向通風方式由于具有更高的通風效率、氣流速度，且噪音更低，加上污濁空氣便于集中消毒，相鄰雞舍間交叉污染少，從而在近年來的雞舍設(shè)計中得到了廣泛應(yīng)用。李保明對縱向通風系統(tǒng)的設(shè)計與應(yīng)用進行了研究，分別討論了雞舍縱向通風系統(tǒng)的通風換氣量與舍內(nèi)風速取值、進氣口面積與位置的確定，風機的選型及其不同季節(jié)的運行管理等問題[2]。但調(diào)查顯示，采用縱向通風的雞舍并不能達到滿意的效果，主要影響因素包括進風口設(shè)計、風機配置及安裝、過流斷面及漏風問題等，這些因素都會影響舍內(nèi)的氣流組織[3]，因此，研究舍內(nèi)氣流組織對縱向通風的設(shè)計和優(yōu)化具有重要意義。

目前關(guān)于籠養(yǎng)雞舍通風效果的研究主要集中于試驗方法[4-8]。由于安裝測試儀器及設(shè)備以及人員進出會影響雞群的生產(chǎn)和生活，導致試驗中監(jiān)測雞舍數(shù)量及舍內(nèi)監(jiān)測布點有限，另外監(jiān)測儀器也存在不可避免的干擾。隨著計算機技術(shù)的發(fā)展，計算流體力學（computer fluid dynamic，CFD）方法由于能得到計算區(qū)域內(nèi)任意點的相關(guān)參數(shù)，且后處理組件可以給出可視化結(jié)果，直觀反映流速或溫度的變化分布，其應(yīng)用越來越廣。目前在農(nóng)業(yè)領(lǐng)域，CFD方法在溫室[9]、豬舍[10-12]、牛舍[13]通風系統(tǒng)的研究中得到很好應(yīng)用；由于平養(yǎng)雞舍內(nèi)設(shè)備較少，計算模型較為簡單，CFD方法也在平養(yǎng)雞舍方面得到推廣[14-16]。而籠養(yǎng)雞舍由于其設(shè)備的復雜性，應(yīng)用CFD方法很少，且一般多采用二維模型[17-18]。本研究采用CFD方法研究中等規(guī)模疊層籠養(yǎng)蛋雞舍內(nèi)的縱向通風系統(tǒng)，分析雞舍內(nèi)的氣流組織及熱濕環(huán)境的三維分布，旨在為疊層籠養(yǎng)蛋雞舍的縱向通風系統(tǒng)設(shè)計提供參考。

1 標準雞舍的建立

根據(jù)調(diào)研及文獻查閱結(jié)果，目前蛋雞舍的飼養(yǎng)規(guī)模大致可分為3種：起步階段存欄2 000～5 000羽，中等規(guī)模飼養(yǎng)階段存欄1萬～5萬羽，大規(guī)模飼養(yǎng)階段存欄10萬羽以上，每棟雞舍最大飼養(yǎng)規(guī)模的推薦值為1.5萬～2萬羽。本研究將建立1個中等規(guī)模、采用縱向通風的標準雞舍，即存欄1萬羽蛋雞，雞舍總長67 m，跨度12 m，檐高3.2 m。雞舍內(nèi)采用2列4層的雞籠。

1.1 縱向通風系統(tǒng)的設(shè)計

通風換氣量根據(jù)雞舍內(nèi)斷面風速進行確定，取斷面風速為1.5 m/s[2]，則雞舍內(nèi)換氣量為：1.5×12×3.2=57.6 m3/s。采用6臺SFT-No10型號的風機，額定風量為32 100 m3/h，安裝在一側(cè)山墻上。進風口面積按風機面積的2倍確定，以減少進氣口阻力，風機和進風口分設(shè)在兩側(cè)山墻上[2]，具體布置見圖1。

1.2 雞籠和雞體模型的建立

為了計算雞籠和雞體對氣流組織的影響，應(yīng)盡量使模型與實際情況相符合，但若按雞籠和雞體的實際尺寸建立模型，須要劃分的網(wǎng)格數(shù)量很大，迭代過程很長，而且會影響迭代過程的收斂性。目前在模擬障礙物對氣流組織的影響時，多采用多孔介質(zhì)模型[19-21]。張?zhí)熘仍谀M二維雞舍的氣流組織時，也采用了多孔介質(zhì)模型簡化舍內(nèi)雞籠對氣流的阻力作用[17]。本研究將單個雞籠簡化為多孔介質(zhì)模型，設(shè)定其阻力系數(shù)來表示雞體對空氣流動的阻礙作用。由于雞籠在3個方向上對空氣的阻力特性均不同，所以應(yīng)分別計算3個方向的阻力系數(shù)。

多孔介質(zhì)模型采用經(jīng)驗公式定義多孔介質(zhì)上的流動阻力。從本質(zhì)上說，多孔介質(zhì)模型是在動量方程中增加了1個代表動量耗散的源項以模擬多孔介質(zhì)的作用。源項由2部分組成：黏性損失項，即方程（1）右端第1項；慣性損失項，方程（1）右端第2項。

3 模擬結(jié)果

通過計算得出，雞舍內(nèi)平均風速為1.4 m/s，與前面的理論計算值的相對誤差為6.7%，在可接受范圍內(nèi)，滿足設(shè)計要求。雞舍內(nèi)平均溫度為31.72 ℃，而從熱應(yīng)激的角度看，21～26 ℃是雞的舒適環(huán)境溫度范圍，26～32 ℃是不太舒適、但能維持正常生理功能的環(huán)境溫度范圍，高于32 ℃是雞生理功能區(qū)域紊亂的溫度范圍[23]。31.72 ℃的平均溫度可以維持蛋雞的基本生理功能，但為了保證其生產(chǎn)能力，還是應(yīng)該采取輔助降溫措施將舍內(nèi)溫度降至25 ℃以下。蛋雞的適宜相對濕度為60%，但相對濕度為45%～70%對蛋雞的生產(chǎn)性能影響不大。本研究計算得出的平均相對濕度為66%，能夠滿足雞生長所需的濕度環(huán)境。

3.1 1.5 m高度處的氣流組織及溫度、濕度分布

研究了第3層雞籠中心高度處即1.5 m處的氣流組織及熱濕環(huán)境分布。從圖4中可以清晰地看到，在雞舍進風端存在一段距離的渦流，該距離大約為8 m左右，經(jīng)過這段距離后，流線達到穩(wěn)定，平直向前；在雞舍末端，由于風機的抽吸作用，流線又匯到一起，最終經(jīng)由風機流出。通過分析風速矢量圖得出，雞舍內(nèi)流場的大部分區(qū)域氣流分布穩(wěn)定，相互間干擾較小。

從圖5可以看出，進口端的風速梯度較大，在達到y(tǒng)=18 m 的距離后，風速基本保持不變，中間走廊的平均風速可達1 m/s，兩側(cè)走廊的風速高于中間走廊，為1.5 m/s，而在雞所處的位置，風速僅為0.5 m/s。雖然雞舍的平均風速達到 1.4 m/s 左右，基本滿足設(shè)計要求，但在雞籠中間部位的風速遠不能滿足斷面風速為1.5 m/s的要求，這可能是由于雞舍前端存在渦流，使一部分新風停滯不前，有效利用的新風量減少，從而減少了斷面風速。因此在設(shè)計雞舍縱向通風系統(tǒng)時，

應(yīng)采取有效措施盡量減少雞舍前端渦流的存在，使通風盡可能被有效利用。

從圖6可以看出，沿長度方向雞舍內(nèi)溫度逐漸升高。在寬度方向上，雞所處位置由于有熱源存在，溫度相對走廊處的要高，且中間雞籠處的溫度要低于兩側(cè)雞籠，這可能是由于雞舍前端存在渦流，而中間位置處渦流范圍較大，氣流分布較為均勻，使得該處溫度梯度較小，即溫升小，導致該部分中間雞籠處溫度低于兩側(cè)雞籠，而由于雞舍中部氣流之間相互干擾較小，兩側(cè)氣流對中間的影響較小，導致沿雞籠的整個長度方向中間雞籠的溫度都低于兩側(cè)雞籠。

從圖7可見，雞舍內(nèi)相對濕度分布與溫度分布類似，與其不同的是，沿長度方向雞舍內(nèi)相對濕度逐漸減小。在寬度方向上，雞所在位置的相對濕度低于走廊處的相對濕度，中間雞籠處的相對濕度高于兩側(cè)雞籠。相對濕度的大小主要取決于舍內(nèi)空氣含濕量及溫度。對照溫度分布云圖可以得出，雞的散濕量對舍內(nèi)相對濕度的影響較小，相對濕度變化主要受溫度的影響，溫度越高，相對濕度越低。

不同高度處水平面的氣流組織及溫度、濕度分布云圖情況類似，因此不詳細列出。

3.2 不同高度處氣流組織分布差異

由于雞舍沿寬度方向是對稱的，所以只對其一半進行研究即可。取3個具有代表性的截面進行研究，分別為x=225 m（第1排雞籠中間位置），x=4.00 m（走道中間位置），x=5.75 m（第2排雞籠中間位置），具體位置如圖8所示。

由于雞舍內(nèi)通風主要是為了滿足雞的生產(chǎn)生活，因此在長度方向上主要研究雞籠存在的范圍，即y為3～64 m的范圍。從圖9至圖11可以看出，垂直平面x=2.25、5.75 m的充分發(fā)展段內(nèi)，風速沿長度方向即y方向的分布相似。在雞籠存在的高度范圍內(nèi)，包括z=0.5、1、1.5、2 m，其風速沿長度方向均保持不變，且4個高度處的風速差別不大，平均速度在0.7m/s左右。在x=2.25m的截面處，在其入口端，由于受渦流的影響，風速分布不穩(wěn)定，部分位置處風速較大，容易給雞造成較強的吹風感，應(yīng)對這部分氣流組織進行優(yōu)化，提高雞的舒適性。在x=5.75 m的截面處，其前部分正好處于渦流中心處，風速較小，在y=7 m左右的位置，風速接近于0，這對雞的散熱是不利的，故這部分氣流組織也應(yīng)進行優(yōu)化。在 x=4.00 m 處由于沒有安置雞籠，其風速分布與其他2個截面處有所不同。在充分發(fā)展段內(nèi)，風速在高度方向呈上升趨勢，且平均速度高于其他2個截面。2 m高度以下風速低于 2 m/s，不會引起雞的不舒適，故這部分氣流基本滿足要求。

3.3 不同高度處溫度、濕度分布差異

雞舍內(nèi)溫度分布受到熱源即雞的散熱量以及風速的影響。由圖12-a、圖13-a、圖14-a可見，沿高度方向，雞舍內(nèi)溫度呈下降趨勢，在雞籠高度范圍內(nèi)，溫度沿高度方向的梯度較小，這是由于位置較高處的風速較大。在散熱量一定的情況下，風速越大，其帶走的熱量越多，熱量在熱源處累積的量就小，故溫度較低。另外發(fā)現(xiàn)，在雞籠高度范圍內(nèi)的溫度差異不是很大，這與風速的分布是對應(yīng)的。在y=40 m之后，雞舍溫度達到了32 ℃以上，不能滿足雞的正常生理需求溫度，因此在采取輔助降溫措施時，應(yīng)將重點放在y=40 m之后的范圍。

從圖12-b、圖13-b、圖14-b中可以看出，沿高度方向，雞舍內(nèi)相對濕度逐漸升高，這與溫度分布剛好呈相反趨勢，且在雞籠高度范圍內(nèi)相對濕度分布差異不是很大，這與前述相對濕度受室內(nèi)溫度影響較大的結(jié)果一致。在整個雞舍內(nèi)相對濕度均能滿足雞舒適生長所需的濕度環(huán)境。

4 結(jié)論

采用CFD方法研究了縱向通風技術(shù)在疊層籠養(yǎng)蛋雞舍

的應(yīng)用效果。結(jié)果表明，雞舍內(nèi)平均溫度為31.72 ℃，平均相對濕度為66%，平均風速為1.4 m/s，為了滿足雞的正常生產(chǎn)需求，還應(yīng)采取一定的降溫措施。在雞籠位置范圍內(nèi)，氣流和熱濕環(huán)境在高度和寬度方向分布較為均勻；而沿長度方向溫度升高1.5 ℃，相對濕度降低5%；在氣流達到穩(wěn)定后走道平均風速為1 m/s，雞籠附近平均風速為0.5 m/s。本研究還發(fā)現(xiàn)一些問題：（1）在氣流入口端風速過大，有些地方達到 6 m/s，這對雞的生長極為不利；（2）在2個進風口之間由于存在渦流，使新風的有效利用率下降；（3）若只采用縱向通風技術(shù)對雞舍進行降溫，為了滿足雞的正常生理需求，在其他條件不變的情況下，雞舍長度應(yīng)小于40 m。綜上，在該雞舍內(nèi)單獨采用縱向通風技術(shù)不能滿足雞正常生產(chǎn)所需的溫度環(huán)境，還應(yīng)采取輔助降溫措施。由于氣流對雞舍內(nèi)熱濕環(huán)境的分布具有重要影響，對氣流分布較差的地方應(yīng)加以改善。

參考文獻：

[1]耿愛蓮，李保明，趙芙蓉，等. 集約化養(yǎng)殖生產(chǎn)系統(tǒng)下肉種雞健康與福利狀況的調(diào)查研究[J]. 中國家禽，2009，31（9）：10-15.

[2]李保明. 雞舍縱向通風系統(tǒng)設(shè)計與應(yīng)用中的幾個問題[J]. 農(nóng)業(yè)工程學報，1993，9（增刊）：82-85.

[3]黃炎坤. 雞舍縱向通風設(shè)計應(yīng)用情況的調(diào)查分析[J]. 山東家禽，2004，40（5）：15-17.

[4]楊 軍，雍芳梅，李 紅. 雞舍縱向通風系統(tǒng)在生產(chǎn)中的應(yīng)用效果[J]. 新疆畜牧業(yè)，1994（3）：49-51.

[5]張?zhí)熘?，黃之棟，潘 強，等. ?？诘貐^(qū)高密度疊層籠養(yǎng)蛋雞舍內(nèi)熱濕環(huán)境的實驗研究[J]. 農(nóng)業(yè)工程學報，1997，13（增刊）：80-84.

[6]姚惠文，周雁華.籠養(yǎng)蛋雞舍內(nèi)的溫度分布[J]. 國外畜牧學：豬與禽，2000（2）：32-34.

[7]趙淑梅，星典宏，黃仕偉，等. 混合通風方式下大型肉雞舍過渡期通風效果測試[J]. 農(nóng)業(yè)工程學報，2006，22（增刊）：96-99.

[8]中國農(nóng)業(yè)工程學會. 2007年中國農(nóng)業(yè)工程學會學術(shù)年會論文摘要集[C]. 北京：農(nóng)業(yè)工程學報雜志編輯部，2007.

[9]程秀花. 溫室環(huán)境因子時空分布CFD模型構(gòu)建及預測分析研究[D]. 鎮(zhèn)江：江蘇大學，2011.

[10]佟國紅，張國強，Morsing S，等. 豬舍內(nèi)氣流變化的模擬研究[J]. 沈陽農(nóng)業(yè)大學學報，2007，38（3）：379-382.

[11]Seo I H，Lee I B，Moon O K，et al. Modelling of internal environmental conditions in a full-scale commercial pig house containing animals[J]. Biosystems Engineering，2012，111（1）：91-106.

[12]賀 城，牛智有，齊德生. 豬舍溫度場和氣流場的CFD模擬比較分析[J]. 湖北農(nóng)業(yè)科學，2010，49（1）：134-136，155.

[13]Wu W，Zhai J，Zhang G，et al. Evaluation of methods for determining air exchange rate in a naturally ventilated dairy cattle building with large openings using computational fluid dynamics（CFD）[J]. Atmospheric Environment，2012，63：179-188.

[14]Seo I H，Lee I B，Moon O K，et al. Improvement of the ventilation system of a naturally ventilated broiler house in the cold season using computational simulations[J]. Biosystems Engineering，2009，104（1）：106-117.

[15]Blanes-Vidal V，Guijarro E，Balasch S，et al. Application of computational fluid dynamics to the prediction of airflow in a mechanically ventilated commercial poultry building[J]. Biosystems Engineering，2008，100（1）：105-116.

[16]李文良，施正香，王朝元. 密閉式平養(yǎng)雞舍縱向通風的數(shù)值模擬[J]. 中國農(nóng)業(yè)大學學報，2007，12（6）：80-84.

[17]張?zhí)熘?，?泳，黃之棟，等. 高密度疊層籠養(yǎng)蛋雞舍內(nèi)氣流場的數(shù)值模擬[J]. 農(nóng)業(yè)工程學報，1999，15（增刊）：38-43.

[18]潘 強，黃之棟，張?zhí)熘? 疊層籠養(yǎng)蛋雞舍夏季環(huán)境特性與熱濕環(huán)境模型[J]. 中國農(nóng)業(yè)大學學報，1998，3（4）：54-58.

[19]張 楠. 基于多孔介質(zhì)模型的城市濱江大道風環(huán)境數(shù)值模擬研究[D]. 長沙：中南大學，2005.

[20]劉元坤. 多孔介質(zhì)模型在大空間建筑氣流組織模擬中的應(yīng)用研究[D]. 長沙：湖南大學，2011.

[21]畢春偉. 養(yǎng)殖網(wǎng)箱周圍流場特性的三維數(shù)值模擬[D]. 大連：大連理工大學，2012.

[22]楊林亞，吳化方，張俊祥，等. 陜西關(guān)中農(nóng)戶籠養(yǎng)蛋雞舍夏季環(huán)境的測定[J]. 畜牧獸醫(yī)雜志，1993，12（4）：11-13.

[23]羅慶斌，何丹林. 雞熱應(yīng)激研究概述[J]. 養(yǎng)禽與禽病防治，2005（4）：2-5