小麥堆孔隙率測定方法及影響因素研究

李倩倩，陳 鑫，畢文雅，姜俊伊，石天玉?

（國家糧食和物資儲備局科學研究院，北京 100037）

糧食孔隙率是糧食流散特性之一，是糧食重要的物理特性[1]。糧堆孔隙分為種粒內孔隙和種粒間孔隙[2]。種粒內孔隙決定了籽粒硬度、易破損粒、研磨、干燥率和抗真菌發展[3]，種粒間孔隙影響糧食呼吸作用、熱量的釋放和蟲霉變化[4]。

糧堆孔隙率對糧堆內部濕熱轉移影響較大，影響通風干燥過程中氣體流動阻力，是糧食儲備管理作業中糧堆溫濕度的控制工藝選擇和優化的重要參數。孔隙率的測定研究對小麥的安全儲存有重要意義。

孔隙率的測定方法有比重瓶法，Chang[5]使用氣體比重瓶測定了不同含水率的玉米、小麥和高粱的孔隙率。目前測定糧食孔隙率的常用方法有氣體置換法，根據理想氣體狀態方程，直接測量糧食的孔隙率。張淑珍[6]依據理想氣體狀態方程研制了孔隙率測定裝置，指出此裝置對吸濕性強、易溶于水、比重測試較復雜的散粒物料更適用。田曉紅[1]通過對比儀器法和計算方法測定小麥、玉米、大豆、稻谷孔隙率，指出基于理想氣體狀態方程的儀器法樣品代表性強，檢測方法簡單、快捷、重復性好，且整個過程是無損操作，便于推廣。李長友[7]基于理想氣體狀態方程設計了定容型散體物料孔隙率的測定裝置，開發了自動檢測控制系統，實現了測定過程的自動控制和檢測參數的實時顯示與存儲。郝倩[8]基于理想氣體狀態方程研究玉米、小麥孔隙度就倉測量法，指出該測量系統測定時無需取樣，可以測出糧食在倉儲狀態下的孔隙度。唐福元[9]、陳雪[10]基于理想氣體狀態方程分別研究了筒倉中稻谷、小麥的空隙率分布規律，為糧倉孔隙率的分布研究提供了新的方法。許倩[3]利用孔隙率測定儀測定了油菜籽孔隙率，并分析含水率和品種對孔隙率的影響。

本文使用孔隙率測定儀和應變控制式三軸儀，根據理想氣體狀態方程，測量并計算不同品種、不同含水率、不同壓力下小麥的孔隙率值，并分析壓力、含水率、品種對孔隙率的影響。為小麥糧堆提高仿真精度，確定合理的通風方式提供理論依據。

1 材料與方法

1.1 實驗材料

見表1和表2。

表1 小麥樣品尺寸Table 1 The size of wheat mm

表2 小麥樣品的含水情況Table 2 The moisture content of wheat %

1.2 實驗儀器

DGG-9140A電熱恒溫鼓風干燥箱：上海森信實驗儀器有限公司；TSZ-6A型應變控制式三軸儀：南京土壤儀器廠；孔隙率測定儀：自制。

1.3 實驗原理

1.3.1 糧食堆密度測定原理

使用應變控制式三軸儀給予糧堆各方向應力，通過應力的改變，觀察糧堆體積的變化，從而得到不同加載壓力下糧堆的密度。糧堆受力情況如圖1所示：

圖1 糧堆應力分布情況Fig.1 The stress distribution of grain pile

根據糧堆質量m即可計算糧堆的初始密度和不同壓力條件下的壓縮密度：

ρ=1 000×m/v，ρ-不同壓力下糧堆壓縮密度（g/L）；v-不同壓力下糧堆的體積（cm3）。

1.3.2 糧食孔隙率測定原理

孔隙率測定儀采用兩個容積均為5 L，且相互可以通過閥門連接的壓力容器A和B，壓力容器A用于壓縮氣體填充，壓力容器B用于測試糧種填裝，記錄兩個壓力容器聯通前后的氣體壓力數據 P1、P2[1，3]，自制孔隙率測定儀裝置示意圖（圖 2）。

圖2 自制孔隙率測定儀裝置示意圖Fig.2 The schematic diagram of self-made porosity tester

1.4 實驗方法

（1）檢查儀器是否齊全，將待測試樣品裝入三軸儀橡皮膜，裝滿密封后安裝壓力室，對壓力室進行注水，壓力室注水完成后，關閉排水口和排氣閥。

（2）利用反壓調節手輪調節壓力室壓力，設置壓力梯度：0、50、100、150、200、250、300 kPa。并記錄加壓后三軸儀體量測試管的體積變化。

（3）加壓結束后將反壓調節手輪卸載壓力至0 kPa，打開排水口和排氣閥，卸載壓力室，進行樣品稱重記錄，用于堆密度計算。

（4）將不同水分樣品進行3組重復測試，計算小麥樣品不同含水率不同壓力下的堆密度。

（5）根據樣品0 kPa壓力下堆密度計算自制孔隙率測定儀樣品室所需樣品質量，將樣品裝滿樣品室，裝好桶蓋進行密封，打開閥門向氣壓室充氣，氣壓室內氣壓穩定后，關閉閥門，打開氣壓室與樣品室連接閥門，直至氣壓穩定，記錄初始及穩定時氣壓數值。

（6）根據自制孔隙率測定儀測定樣品在0 kPa下的孔隙率及樣品在不同壓力下的堆密度，進行樣品在不同壓力下的孔隙率計算。

2 結果與分析

2.1 小麥的孔隙率值

使用孔隙率測定儀和應變控制式三軸儀，測量并計算濟麥 22、寧麥 13、淮麥 20分別在不同含水率、不同壓力下的孔隙率，結果如表 3～5所示。

濟麥 22在含水率 11.40%～17.98%、壓力0～300 Pa時，小麥孔隙率變化范圍是 36.99%～49.73%；寧麥13在含水率11.23%～17.27%、壓力0～300 Pa時，小麥孔隙率的變化范圍是36.73%～48.49%；淮麥20在含水率10.27%～17.38%、壓力0～300 Pa時，小麥孔隙率的變化范圍是40.95%～52.83%。

表3 不同含水率、不同壓力下濟麥22小麥樣品孔隙率Table 3 The porosity of Jimai 22 wheat with different water content and pressure

表4 不同含水率、不同壓力下寧麥13小麥樣品孔隙率Table 4 The porosity of Ningmai 13 wheat with different water content and pressure

表5 不同水分、不同壓力下淮麥20小麥樣品孔隙率Table 5 The porosity of Huaimai 20 wheat with different water and pressure

2.2 壓力對孔隙率的影響

根據表3～5可做不同含水率小麥的孔隙率隨壓力變化圖，如圖3～5。

三種小麥的孔隙率隨壓力的變化均成冪函數關系，都可擬合成方程：y=kpm，y-樣品孔隙率（%），p-樣品所受壓力（Pa）。壓力與孔隙率的擬合方程系數如表6～8所示。

圖4 不同含水率寧麥13的孔隙率隨壓力的變化Fig.4 The porosity of different water content of Ningmai 13 under different pressure

圖5 不同含水率淮麥20的孔隙率隨壓力的變化Fig.5 The porosity of different water content of Huaimai 20 under different pressure

表6 濟麥22壓力與孔隙率的擬合方程系數Table 6 Coefficient of fitting equation between pressure and porosity of Jimai 22

表8 淮麥20壓力與孔隙率的擬合方程系數Table 8 Coefficient of fitting equation between pressure and porosity of Huaimai 20

從圖 3～5、表 6～8分析，可以看出在一定水分下，小麥的孔隙率隨壓力的增加而減小，小麥所受壓力與孔隙率成冪函數關系。濟麥22在壓力0～300 Pa下，含水率 11.40%的小麥孔隙率從45.14%下降到 36.99%，含水率 17.98%的小麥孔隙率從 49.73%下降到 39.81%；寧麥 13在壓力0～300 Pa下，含水率 11.23%的小麥孔隙率從44.85%下降到 36.73%，含水率 17.27%的小麥孔隙率從 48.49%下降到 40.42%；淮麥 20在壓力0～300 Pa下，含水率 10.27%的小麥孔隙率從47.98%下降到 40.95%，含水率 17.38%的小麥孔隙率從 52.83%下降到 42.45%。小麥的孔隙率隨壓力的變化幅度是逐漸減小，與程緒鐸[11]的研究結果類似。程緒鐸對玉米堆進行壓縮特性進行研究，發現壓縮密度隨著壓力的增大變化率變小，根據孔隙率與壓縮密度的關系，可以推導出孔隙率隨著壓力的增大變化率變小。導致這種現象的原因可能是小麥顆粒隨著壓力增加，前期種粒間孔隙空間重排占主導因素，后期種粒內孔隙空間重排占主導因素。

2.3 含水率對孔隙率的影響

根據表3～5可做不同壓力下小麥的孔隙率隨含水率變化圖，如圖6～8。

圖6 不同壓力下濟麥22的孔隙率隨含水率的變化Fig.6 The porosity of different pressure of Jimai 22 with different water content

圖7 不同壓力下寧麥13的孔隙率隨含水率的變化Fig.7 The porosity of different pressure of Nimai 13 with different water content

圖8 不同壓力下淮麥20的孔隙率隨含水率的變化Fig.8 The porosity of different pressure of Huaimai 20 with different water content

三種小麥的孔隙率隨含水率變化均成二次函數關系，都可擬合成方程：z=Ax2+Bx+C，z-樣品孔隙率（%），x-樣品含水率（%）。含水率與孔隙率的擬合方程系數如表9～11所示。

表9 濟麥22含水率與孔隙率的擬合方程系數Table 9 Coefficient of fitting equation between water content and porosity of Jimai 22

表10 寧麥13含水率與孔隙率的擬合方程系數Table 10 Coefficient of fitting equation between water content and porosity of Nimai 13

表11 淮麥20含水率與孔隙率的擬合方程系數Table 11 Coefficient of fitting equation between water content and porosity of Huaimai 20

從圖6～8、表9～11分析，可以看出在一定壓力下，小麥的孔隙率隨含水率增加有增加的趨勢。原因可能由于籽粒膨脹導致籽粒體積變大、籽粒間空隙變大，從而孔隙率增加，顆粒尺寸不均勻性也有一定影響。

2.4 品種對孔隙率的影響

實驗選取 3種品種的小麥進行孔隙率的測定，研究品種對孔隙率的影響。小麥尺寸在長軸6.275 6～6.514 0 mm、中軸 3.333 1～3.547 0 mm、短軸2.848 9 mm～3.262 3 mm時，孔隙率隨壓力的變化成冪函數關系，孔隙率隨水分的變化成二次函數關系。可以推測在一定尺寸范圍內，品種對孔隙率沒有明顯影響。這與徐倩[3]的研究結果相同，她認為油菜籽的品種對孔隙率無明顯影響。小麥樣品的幾何尺寸、體積與結構對孔隙率應該有一定的影響。

3 討論與結論

小麥的孔隙率隨壓力的增加而減小，成冪函數關系，變化幅度逐漸減小。濟麥 22在壓力 0～300 Pa下，含水率11.40%的小麥孔隙率從45.14%下降到 36.99%，含水率 17.98%的小麥孔隙率從49.73%下降到39.81%；寧麥13在壓力0～300 Pa下，含水率11.23%的小麥孔隙率從44.85%下降到36.73%，含水率 17.27%的小麥孔隙率從 48.49%下降到40.42%；淮麥20在壓力0～300 Pa下，含水率 10.27%的小麥孔隙率從 47.98%下降到40.95%，含水率 17.38%的小麥孔隙率從 52.83%下降到42.45%。

小麥的孔隙率隨水分的增加而增大，成二次函數關系；在一定尺寸范圍內，小麥的品種對孔隙率無明顯影響；筒倉中的小麥堆結構、應力和應變復雜，測量孔隙率的分布是個難題。本文測量的孔隙率在實驗后期，采用合適數值方法可模擬計算出筒倉中小麥堆孔隙率隨糧層深度的規律，從而得到筒倉中孔隙率值的分布。