氣壓深松土壤孔隙度測試與分析

左勝甲，孔德剛※，劉春生，李紫輝，張 超，陳 爽，吳艷鳳

（1. 東北農業大學工程學院，哈爾濱 150030；2. 黑龍江科技大學機械工程學院，哈爾濱 150022）

氣壓深松土壤孔隙度測試與分析

左勝甲1，孔德剛1※，劉春生2，李紫輝1，張 超1，陳 爽1，吳艷鳳1

（1. 東北農業大學工程學院，哈爾濱 150030；2. 黑龍江科技大學機械工程學院，哈爾濱 150022）

為了明晰氣壓深松土壤孔隙度的變化規律及其效果，利用溫納電阻率法測試氣壓深松前、后的土壤孔隙度，并以孔隙度增加率為指標進行分析，結果表明：深松氣壓和水平距離對孔隙度增加率有極顯著影響（P＜0.01），而犁底層容重沒有顯著影響；水平距離0.7 m內的土壤孔隙度增大，犁底層孔隙度增加率隨水平距離增加呈現先緩慢變化在0.4 m處出現峰值后，急速減小的變化趨勢，水平距離0.1～0.4 m的孔隙度增加率均值顯著大于0.5～0.6 m的孔隙度增加率均值（P＜0.05）；深松氣壓2.2和1.8 MPa作用下的犁底層孔隙度增加率差異不顯著，但均顯著大于1.4 MPa作用下的犁底層孔隙度增加率（P＜0.05），1.8 MPa為最佳深松氣壓。結論可為氣壓深松技術及其設備的研發提供依據。

孔隙度；壓力；土壤；氣壓深松；溫納測試法

0 引 言

實施深松作業可打破長期機械化作業形成的犁底層，使犁底層的土壤得到松動，能調節土壤三相比，提高土壤的孔隙度和透水、透氣性能，減少水土流失，有利于農業生態環境保護。傳統的深松方式主要是用深松鏟直接松動耕地犁底層土壤，這種方式存在深松深度淺，影響范圍小等問題[1]。氣壓深松利用氣壓劈裂原理，向耕地土壤內部注入高壓氣體，使土壤中空氣量增加，改變土壤三相比，同時也使土壤產生裂隙，導致土壤孔隙度增大，從而實現對耕地的深松，其具有深松深度深，影響范圍大等特點[2]。目前中國關于氣壓深松方面的研究很少，特別是氣壓深松土壤孔隙度隨深松參數的變化規律及其效果亟待研究。

深松效果可用深松后土壤孔隙度及其增加率來評價。目前針對土壤孔隙度的測量方法主要有傳感器測量、圖像分析及取樣分析 3 種[3]。其中傳感器測量法多為侵入式測量，探頭及傳輸線路會干擾被測土壤的物理特性[4]，圖像分析法測量效果及精度受土壤類型及觀測角度、距離等因素影響顯著[5-6]，取樣分析法需直接擾動土壤，無法實現實時準確測量及分析。近年來，溫納（Wenner）電阻率法在土壤測試中的應用增多，它是通過在被測量介質中布設電極，測定土體電阻借助反演計算方法獲得電阻率，并利用電阻率與土體物理狀態參數的定量關系來表征和測量某一物理特性，是一種無損、多維、多尺度的實時測量方法。如 C hambers等[7]、劉漢樂等[8]利用此法研究了非水相有機物在土層中的運移擴散規律；Greve等[9]、Hassan等[10]利用此法研究了土體裂縫展開規律；張志祥等[11]利用此法研究確定了土壤水分特征曲線。

因此，本文利用溫納電阻率法，測得氣壓深松前、后土壤中各測試點的土壤電阻率，并利用阿爾奇公式求其土壤孔隙度，進一步計算土壤孔隙度增加率，分析研究氣壓深松后土壤內的孔隙度增加率隨深松參數的變化規律，為氣壓深松技術及其設備的研發提供技術依據。

1 氣壓深松測試原理

土壤電阻率與含水量、孔隙度、結構及溫度等多種因素有關，在其他因素不變或相對穩定時，可通過測量土體的電阻率間接測量某單一因素的變化及分布特征[12-13]。

由拉普拉斯方程及歐姆定律可知一個點電流源電場中，均勻介質中任一點的電位[14-16]。

式中U為電位，V；ρ為電阻率，?·m；I為點電流源電流，A；r為任一點到點電流源的距離，m。在2個異性點電流源電場中（圖1），按疊加原理，任一點的電位為點電流源A（+I）和B（-I）在該點所產生電位的矢量和。

圖1 2個異性點電流源的電場Fig.1 Electric field of 2 opposite point current source

由式（1）可得，M點的電位為

式中UM為M點的電位，V；AM為M點到電源A的距離，m；BM為M點到電源B的距離，m。

利用電阻率法測定土壤電阻率時（圖 2），先通過電測儀器測定MN電極間的電位差和AB回路的供電電流，再計算其電阻率。

圖2 任意4電極裝置示意圖Fig.2 Schematic diagram of any quadrupole device

任意四極裝置中的AB為供電電源，MN之間的電位差為

式中ΔUMN為MN之間的電位差，V；AN為N點到電源A的距離，m；BN為N點到電源B的距離，m。

則電阻率的公式為

式中k為電極排列系數，m。本試驗采用的是溫納4電極測試法，其電極位于同一水平線上且間距相等，所以k=2π·a，a為相鄰電極間的間距，m。

由溫納法得到的電阻率是深度為電極間距、位置為電勢電極間連線中點處介質的綜合電阻率，也稱為視電阻率[17]。

2 氣壓深松試驗

2.1 供試土壤及試驗儀器

2016年1—4月進行試驗。試驗用土壤采自東北農業大學香坊農場（45°45′N、128°39′E），為使試驗用土壤的含水率與實際深松作業時的土壤含水率（15%～22%）[18-20]保持一致，試驗前將其含水率調制為18%±0.5%，為了模擬實際深松作業時的土層狀態，在土槽底部設置犁底層，試驗前用石磙碾壓使其容重為 1 .4～1.8 g/cm3[21-22]，厚度為0.2 m，碾壓過程中用環刀法監測土壤容重。然后在犁底層上面覆蓋約0.2 m厚的試驗用土壤作為耕作層，并使其表面平整光滑便于測試。

試驗所用設備為自制的氣壓可調式氣壓深松試驗裝置，主要由高壓氣泵、土槽、支撐架、調壓閥、氣槍開關和氣槍等組成（圖3）。高壓氣泵1產生的高壓氣體經過調壓閥2、高壓輸氣管4、氣槍開關5，由氣槍6噴氣孔注入土槽7的土壤之中，使土壤孔隙度增加實現深松。高壓氣泵提供的氣壓為0.8～3.0 MPa，通過調壓閥調節，氣槍可得到不同的深松氣壓。其中氣槍由無縫鋼管制成，頂端為圓錐體，前部均布4個直徑為5 mm的噴氣孔，后部通過高壓輸氣管、調壓閥與高壓氣泵連接。

圖3 氣壓深松試驗裝置示意圖Fig.3 Schematic of pneumatic subsoiling test equipment

試驗儀器有電測儀器、硬度計、環刀、天平和計算機等。

2.2 試驗設計及過程

為分析不同容重土壤中深松氣壓對土壤空隙度的影響，設置3個犁底層容重水平（1.4、1.6、1.8 g/cm3）和3個深松氣壓（1.4、1.8、2.2 kPa）水平，進行全組合氣壓深松試驗，并利用溫納法測試各測試點深松前、后的土壤電阻率。每項組合重復3次，取3次試驗結果的平均值作為分析計算值。

利用溫納法對各測試點的土壤電阻率進行測試，測試時在土面中心將氣槍插入犁底層中，使噴氣孔距土面0.35 m，并固定在支撐架上。測試點的位置由坐標（x，y）給出，其中x為水平距離（距土面中心的水平尺寸），分別取值0.1，0.2，0.3，0.4，0.5，0.6和0.7 m；y為深度（距土面的垂直尺寸），分別取值0.1，0.15，0.3和0.35 m。

2.3 指標測定及計算

2.3.1 土壤孔隙度計算

1）確定試驗用土壤電阻率與孔隙度的計算公式。利用阿爾奇建立的土壤電阻率模型[23]

式中Ф為土壤孔隙度；Sr為土壤飽和度，%；K為與土壤特性有關的常數系數；m和n分別為模型待定常數。

試驗前，采用供試土壤，設置不同飽和度和孔隙度試驗進行電阻率測試，對孔隙率與電阻率模型參數進行標定。通過回歸分析確定式（5）中m、n、K。得到其電阻率與孔隙度的計算公式

式中W為土壤含水率，%。決定系數R2為0.96（P＜0.05）。

本試驗W為18%，將W及各測試點深松前后的土壤電阻率ρ的數值代入到式（6）中，即可得到各測試點深松前后的土壤孔隙度。

2.3.2 土壤孔隙度增加率計算

氣壓深松使土壤中空氣的含量增加，引起土壤體積膨脹、土面抬升，同時使土壤產生裂隙，導致土壤孔隙度增大。為了便于分析比較不同條件下的深松效果，本文以土壤孔隙度增加率（試驗前、后土壤孔隙度之差與試驗前土壤孔隙度比值的百分比）為指標，分析研究土壤孔隙度增加率的變化規律并確定合理的技術參數。

3 結果與分析

3.1 土壤孔隙度增加率方差分析

通常深松作業時犁底層的容重為 1 .4～1.8 g/cm3[24-27]。對不同容重土壤犁底層在不同深松氣壓作用下各水平距離的孔隙度增加率進行方差分析，結果表明，深松氣壓極顯著影響土壤孔隙率增加率（F=6.751＞F(2,44)=4.98），水平距離也極顯著影響土壤孔隙率增加率（F=4.98＞F(5,44)=3.65），而犁底層容重沒有顯著影響。可見，深松氣壓和水平距離是影響土壤孔隙度的重要影響因素。 因此，下文主要針對犁底層容重為1.6 g/cm3（平均值）的土壤，分析不同深松氣壓和水平距離孔隙度增加率的分布規律。

3.2 不同水平距離孔隙度增加率分布規律

根據試驗所得數據，利用Matlab軟件，通過以三角形為基礎的三次方程內插法，得到不同深松氣壓下孔隙度增加率在土壤中不同深度不同水平距離的分布如圖4所示。

圖4 不同深松氣壓下土壤孔隙度增加率分布Fig.4 Porosity increase rate distribution under different pressures of pneumatic subsoiling

3.2.1 深松氣壓1.4 MPa下的孔隙度增加率分布

由圖4a可知，深松氣壓為1.4 MPa時，孔隙度增加率在耕作層（0.2 m以上部分，以下不再注明）內隨水平距離的增加呈減小趨勢，在水平距離0.7 m時減小為0；在犁底層（0.2～0.35 m，以下不再注明）內隨水平距離增加呈現先波動變化，在0.4 m處出現波峰后急速減小的變化趨勢。在相同水平距離時，孔隙度增加率隨深度的增加，在耕作層內變化很小，而在犁底層內變化較大；在深度約0.3 m，水平距離約0.4 m處，孔隙度增加率出現最大值，約為5.1 %。

3.2.2 深松氣壓1.8 MPa下的孔隙度增加率分布

由圖4b可知，深松氣壓為1.8 MPa時，在不同深度下，孔隙度增加率呈隨水平距離增加而減小的趨勢。在犁底層內呈現先是緩慢變化，在0.4 m處急速出現波峰后急速減小，之后再次緩慢變化的趨勢；而在耕作層內減小的速度相對較快，在水平距離0.7 m時減小為0。在相同水平距離下，孔隙度增加率隨深度增加呈增大趨勢，但在耕作層內增大的速度較快；而在犁底層內增大的速度相對較慢。在深度為0.3～0.35 m，水平距離約為0.4 m時，孔隙度增加率達到最大值，約為6.0 %。

3.2.3 深松氣壓2.2 MPa時的孔隙度增加率分布

由圖4c可知，深松氣壓為2.2 MPa時，孔隙度增加率在耕作層內隨水平距離的增加呈波動減小的趨勢，在0.7 m時減小為0；而在犁底層內隨水平距離增加呈現先緩慢變化，在0.4 m處急速出現波峰后急速減小的變化趨勢。在相同水平距離下，孔隙度增加率隨深度的增加，在耕作層內變化較大，而在犁底層內變化相對較小。在深度0.3～0.35 m，水平距離約0.4 m處，孔隙度增加率出現最大值，其值約為7.5%。

由上分析可知，在不同深松氣壓作用下，土壤孔隙度增加率在水平距離0.7 m范圍內均大于0，表明在此范圍得到了不同程度的深松；孔隙度增加率最大值均出現在水平距離約0.4 m，犁底層深度約0.3～0.35 m處；這個深度與氣槍噴氣孔深度位置相同，說明與氣槍噴氣孔深度相同的犁底層可得到最佳的深松效果。另外，孔隙度增加率在犁底層和耕作層無顯著差異，而在水平方向上則是0.1～0.4 m顯著大于0.5～0.6 m（P＜0.05），表明在0.4 m之后深松效果顯著減小。

3.3 深松氣壓對犁底層孔隙度增加率的影響

根據試驗結果，得到容重為1.6 g/cm3的犁底底層（深度0.3 m），不同深松氣壓作用下不同水平距離的孔隙度增加率，如圖5所示。由圖5可以看出，犁底層在深松氣壓1.4、1.8和2.2 MPa深松后，在半徑0.6 m范圍內任一水平距離下的孔隙度增加率，均以1.8 MPa的最大，1.4 MPa最小（P＜0.05），2.2 與1.8 MPa無顯著差異，但二者均顯著大于1.4 MPa的孔隙度增加率（P＜0.05）。表明在水平距離0.6 m范圍內，較大的深松氣壓產生較大孔隙度增加率，可獲得較好的深松效果。多項研究成果表明，氣壓劈裂時裂隙的擴展方向垂直于最小主應力方向，淺層土的氣壓劈裂其最小主應力為豎向，裂隙為水平裂隙[28-30]。本試驗的氣壓深松屬于淺層土的氣壓劈裂，其最小主應力為豎向，容重為1.6 g/cm3犁底層，因其土壤凝聚力相對較小，使其豎向主應力合力也相對較小，氣體在產生橫向水平擴散的同時，會有部分產生豎向擴散而泄漏出地表。 對于3種深松氣壓，1.4 MPa的深松氣壓相對較小，橫向擴散勢力較弱擴散到土壤中的氣體量較小，所以其孔隙度增加率較低；而1.8、2.2 MPa的深松氣壓相對較高，氣體橫向擴散勢力增強，在土壤中易于形成裂紋增多，擴散到土壤中的氣體量較多，孔隙度增加率較大。但同時其豎向擴散勢力相對較強，引起的氣體泄漏量較大，導致 2 種深松氣壓下的孔隙度增加率的差異不顯著。另外，由1.4、1.8和2.2 MPa引起的孔隙度增加率在半徑0.6 m范圍的平均值分別為2.80，4.93和4.41，主要考慮土壤疏松程度兼顧能耗，選用1.8 MPa的深松氣壓可獲得較好的深松效果。

圖5 不同深松氣壓下犁底層孔隙度增加率Fig.5 Porosity increase rate of plow pan under different pressures of pneumatic subsoiling

4 結 論

1）深松氣壓和水平距離對犁底層孔隙度增加率有極顯著影響，而犁底層的容重對其沒有顯著影響。

2）氣壓深松后，水平距離0.7 m內的土壤孔隙度增大，犁底層內孔隙度增加率隨水平距離增加呈現先緩慢變化，在0.4 m處出現峰值后急速減小的變化趨勢，且水平距離0.1～0.4 m的孔隙度增加率均值顯著大于0.5～0.6 m的孔隙度增加率均值。

3）深松氣壓2.2和1.8 MPa作用下的犁底層孔隙度增加率差異不顯著，但均顯著大于1.4 MPa作用下的孔隙度增加率。考慮土壤疏松程度兼顧能耗，選用1.8 MPa的深松氣壓最佳。

本文對氣壓深松技術的深入研究及其設備的研發有非常實際的意義，但試驗中深松噴氣點的數量較少，僅有一個噴氣點，在后續的研究中應該對多噴氣點下的氣壓深松效果及土壤內的孔隙度的變化規律作進一步研究。

[1] 左勝甲，孔德剛，劉立意，等. 基于氣壓劈裂原理的氣壓深松效果試驗研究[J]. 農業工程學報，2016，32(1)：54－60. Zuo Shengjia,Kong Denggang,Liu Liyi,et al. Experiment on effect of air-pressure subsoiling based on air-pressure cracking theory[J]. Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE),2016,32(1):54－60.(in Chinese with English abstract)

[2] 孔德剛. 一種機械式氣壓深松機：201420516606.6[P]. 2014－11－24.

[3] Lynch R J,Allersma H G B,Barker H,et al. Development of sensors,probes and imaging techniques for pollutant monitoring in geo-environmental model tests[J]. International Journal of Physical Modelling in Geotechnics,2001,1(4):17－27.

[4] Taylor R N. Geotechnical Centrifuge Technology[M]. London UK:Chapman &Hall,1995.

[5] Pasha A Y,Hu L M,Meegoda J N,et al. Centrifuge modeling of in situ surfactant enhanced flushing of diesel contaminated soil[J]. ASTM Geotechnical Testing Journal,2011,34(6):623－633.

[6] Friedman S P. Soil properties influencing apparent electrical conductivity:A review[J]. Computers and Electronics in Agriculture,2005,46(1):45－70.

[7] Chambers J E,Loke M H,Ogilvy R D,et al.Noninvasive monitoring of DNAPL migration through a saturated porous medium using electrical impedance tomography[J]. Journal of Contaminant Hydrology,2004,68(4):1－22.

[8] 劉漢樂，周啟友，吳華橋. 輕非水相液體污染過程的高密度電阻率成像法室內監測[J]. 地球物理學報，2008，51(4)：1246－1254. Liu Hanle,Zhou Qiyou,Wu Huaqiao. Laboratorial monitoring of the LNAPL contamination process using electrical resistivity tomography[J]. Chinese Journal of Geophysics,2008,51(4):1246－1254.(in Chinese with English abstract)

[9] Greve A K,Acworth R L,Kelly B F. Detection of subsurface soil cracks by vertical anisotropy profiles of apparent electrical resistivity[J]. Geophysics,2010,75(4):85－93.

[10] Hassan A,Toll D G. Electrical resistivity tomography for characterizing cracking of soils[C]//Stability and Performance of Slopes and Embankments III. San Diego,California:Geo-Congress,2013:818－827.

[11] 張志祥，徐紹輝，崔峻嶺，等. 電阻率法確定土壤水分特征曲線初探[J]. 土壤，2013，45(6)：1127－1132. Zhang Zhixiang,Xu Shaohui,Cui Junling. et al. Resistivity of soil moisture characteristic curve based on resistivity[J]. Soils,2013,45(6):1127－1132.(in Chinese with English abstract)

[12] Samouelian A,Cousin I,Tabbagh A,et al. Electrical resistivity survey in soil science:A review[J].Soil and Tillage Research,2005,83(2):173－193.

[13] 李金銘. 地電場與電法勘探[M]. 北京，地質出版社，2005.

[14] 劉國華，王振宇，黃建平. 土的電阻率特性及其工程應用研究[J]. 巖土工程學報，2004，26(1)：83－87. Liu Guohua,Wang Zhengyu,Huang Jianping. Research on electrical resistivity feature of soil and it’s application[J]. Chinese Journal of Geotechnical Engineering,2004,26(1):83－87.(in Chinese with English abstract)

[15] 劉盛東，吳榮新. 高密度電阻率法觀測煤層上覆巖層破壞[J]. 煤炭科學技術，2001，29(4)：18－21. Liu Shengdong,Wu Rongxin. High density electric resistance method applied to monitor and measure overburden failure above seam[J]. Coal Science and Technology,2001,29(4):18－21.(in Chinese with English abstract)

[16] 劉國興. 電法勘探原理與方法[M]. 北京：地質出版社，2005.

[17] 李瑛，龔曉南，郭彪，等. 電滲軟黏土電導率特性及其導電機制研究[J]. 巖石力學與工程學報，2009，29(2)：18－21. Li Ying,Gong Xiaonan,Guo Biao,et al. Research on conductivity characteristics of soft clay during electroosmosis and its conductive mechanism[J]. Chinese Journal ofRock Mechanics and Engineering,2009,29(2):18－21.(in Chinese with English abstract)

[18] 呂東山. 機械化深松技術的實施要點[J]. 農業科技與裝備，2012，220(10)：18－21. Lü Dongshan. The key points for the implementation of mechanized deep loosening Technology[J]. Agricultural Science &Technology and Equipment,2012,220(10):18－21.(in Chinese with English abstract)

[19] 李榮，侯賢清. 深松條件下不同地表覆蓋對馬鈴薯產量及水分利用效率的影響[J]. 農業工程學報，2015，31(20)：115－123. Li Rong,Hou Xianqing. Effects of different ground surface mulch under subsoiling on potato yield and water use efficiency[J]. Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE),2015,31(20):115－123.(in Chinese with English abstract)

[20] 李靜靜，李從鋒，李連祿，等. 苗帶深松條件下秸稈覆蓋對春玉米土壤水溫及產量的影響[J]. 作物學報，2014，40(10)：1787－1796. Li Jingjing,Li Congfeng,Li Lianlu,et al. Effect of straw mulching on soil temperature,soil moisture and spring maizeyield under seedling strip subsoiling[J]. Acta Agronomica Sinica,2014,40(10):1787－1796.(in Chinese with English abstract)

[21] 孫一源，高行方，余登苑. 農業土壤力學[M]. 北京：農業出版社，2000：8－15.

[22] 張瑞富，楊恒山，高聚林，等. 深松對春玉米根系形態特征和生理特性的影響[J]. 農業工程學報，2015，31(5)：78－84. Zhang Ruifu,Yang Hengshan,Gao Julin,et al. Effect of subsoiling on root morphological and physiological characteristics of spring maize[J]. Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE),2015,31(5):78－84.(in Chinese with English abstract)

[23] 李良福. 土壤電學[M]. 北京：氣象出版社，2008：100－103.

[24] 中國科學院南京土壤研究所土壤物理研室. 土壤物理性質測定法[M]. 北京：科學出版社，1978：10－11.

[25] 李洪文，陳君達，李問盈. 保護性耕作條件下深松技術研究[J]. 農業機械學報，2000，31(6)：42－45. Li Hongwen,Chen Junda,Li Wenying. Study on subsoiling technique for conservation tillage field[J]. Transactions of the Chinese Society of Agricultural Machinery,2000,31(6):42－45.(in Chinese with English abstract)

[26] 朱鳳武，王景利，潘世強. 土壤深松技術研究進展[J]. 吉林農業大學學報，2003，25(4)：57－61. Zhu Fengwu,Wang Jingli,Pan Shiqiang. Advances in researches of soil subsoiling technique[J]. Journal of Jilin Agricultural University,2003,25(4):57－61.(in Chinese with English abstract)

[27] 姚宇卿，呂軍杰，張潔，等. 深松覆蓋對旱地冬小麥產量和水分利用率的影響[J]. 河南農業科學，2012，41(4)：20－24. Yao Yuqing,Lü Junjie,Zhang Jie,et al. Effects on winter wheat yield and water using efficiency on dryland under long-term subsoiling and mulch tillage[J]. Journal of Henan Agricultural Sciences,2012,41(4):20－24.(in Chinese with English abstract)

[28] 韓文君，劉松玉，章定文，等. 壓力控制的圓孔擴張數值模擬分析[J]. 巖土力學，2010，31(1)：405－410. Han Wenjun,Liu Songyu,Zhang Dingwen,et al. Numerical simulation of pressure-controlled cavity expansion[J]. Rock and Soil Mechanics,2010,31(1):405－410.(in Chinese with English abstract)

[29] Suthersen S S. Remediation Engineering:Design Concept[M]. Boca Raton:CRC Press LLC,1999:238－239. [30] Murdoch L C. Mechanical analysis of idealized shallow hydraulic fracture[J]. Journal of Geotechnical and Geoenvironmental Engineering,2002,128(6):488－495.

Test and analysis of soil porosity with pneumatic subsoiling

Zuo Shengjia1,Kong Degang1※,Liu Chunsheng2,Li Zihui1,Zhang Chao1,Chen Shuang1,Wu Yanfeng1
(1. College of Engineering,Northeast Agricultural University,Harbin 150030,China;2. School of Mechanical Engineering,Heilongjiang University of Science and Technology,Harbin 150022,China)

Pneumatic subsoiling is a novel subsoiling method based on pneumatic split technology,which injects pressure-air into cultivated soil to form fractures. In order to study the change of soil porosity during pneumatic subsoiling based on Wenner electrical sounding method,a labroatory experiment was conducted in Agricultural Machinery Laboratory of the Engineering Research Center of Northeast Agricultural University from January to April in 2016. Different soil bulk density(1.4,1.6 and 1.8 g/cm3) and subsoiling pressure(1.4,1.8 and 2.2 kPa) were designed for the experiment. The initial soil moisture was kept at 18% for all the treatments. At the bottom of the soil bin,a plowpan in depth of 0.2 m was artificially designed by rolling. On the top,a plough layer was of the depth of 0.2 m. The resistivity was measured before and after subsoiling. Then the resistivity was converted into soil porosity by the model established based on a pre-experiment. Then,the porosity increase rate was used to evaluate the change of soil porosity. The results showed that the pressure had an extremely significant influence on soil porosity increase rate(P＜0.01) but the bulk density did not significantly affect the results. Thus,the change of soil porosity as affected by pressure was further investigated in soil with bulk density of 1.6 g/cm3. After subsoiling,the soil porosity increased within the horizontal distance of 0.7 m. In the plowpan(0.2-0.35 m),the porosity increase rate increased slowly with the horizontal distance before 0.4 m and reached its peak when the horizontal distance was 0.4 m. The maximum porosity increase rate was 5.1%,6.0% and 7.5% under the pressure of 1.4,1.8 and 2.2 Mpa,respectively. After the peak,the porosity increase rate rapidly decreased with the increase of the horizontal distance. The average of the porosity increase rate at the horizontal distance of 0.1-0.4 m was higher than that in the 0.5-0.6 m,indicating small subsoiling effect for the horizontal distance larger than 0.4 m. The porosity increase rate under the air pressure of 1.8 and 2.2 MPa were significantly higher than that with the air pressure of 1.4 MPa(P＜0.05). However,there was no significant difference between 1.8 and 2.2 MPa. Considering the effects of subsoiling and energy consumption for subsoiling,we suggested that the subsoiling air pressure was 1.8 MPa when the bulk density was 1.6 g/cm3. The results provides the technical support for development and application of pneumatic subsoiling equipment.

porosity;pressure;soils;pneumatic subsoiling;Wenner test method;

10.11975/j.issn.1002-6819.2017.01.022

S222.1+9

A

1002-6819(2017)-01-0162-05

左勝甲，孔德剛，劉春生，李紫輝，張 超，陳 爽，吳艷鳳. 氣壓深松土壤孔隙度測試與分析[J]. 農業工程學報，2017，33(1)：162－166.

10.11975/j.issn.1002-6819.2017.01.022 http://www.tcsae.org

Zuo Shengjia,Kong Degang,Liu Chunsheng,Li Zihui,Zhang Chao,Chen Shuang,Wu Yanfeng.Test and analysis of soil porosity with pneumatic subsoiling[J]. Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE),2017,33(1):162－166.(in Chinese with English abstract)doi：10.11975/j.issn.1002-6819.2017.01.022 http://www.tcsae.org

2016-05-18

2016-10-10

國家科技支撐計劃課題（2014BAD06B04）

左勝甲，男，吉林通化人，博士生，主要從事農業機械化工程研究。哈爾濱 東北農業大學工程學院，150030。

Email：zuoshengjia1984@163.com

※通信作者：孔德剛，男，吉林省白山市人，教授，博士生導師，主要從事農業機械化工程研究。哈爾濱 東北農業大學工程學院，150030。Email：kong-degang@hotmail.com