• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    水稻土深松阻力與土壤擾動效果研究

    2017-02-08 01:26:44丁啟朔葛雙洋李毅念何瑞銀
    農(nóng)業(yè)機械學(xué)報 2017年1期
    關(guān)鍵詞:耕深土塊耕作

    丁啟朔 葛雙洋 任 駿 李毅念 何瑞銀

    (南京農(nóng)業(yè)大學(xué)工學(xué)院, 南京 210031)

    水稻土深松阻力與土壤擾動效果研究

    丁啟朔 葛雙洋 任 駿 李毅念 何瑞銀

    (南京農(nóng)業(yè)大學(xué)工學(xué)院, 南京 210031)

    為探討機械深松過程的水稻土擾動、土壤結(jié)構(gòu)及能量的效應(yīng),利用原位土壤耕作綜合測試平臺進行精準控制條件的深松試驗,采用鑿形鏟,設(shè)置5種耕作深度(10、15、20、25、30 cm),從耕作阻力及比阻、土壤宏觀擾動輪廓、土壤破碎體尺度分布指標綜合評價水稻土深松作業(yè)的效果。結(jié)果表明,深松耕作阻力與耕作深度符合二階函數(shù)擬合的遞增關(guān)系。耕深20 cm時,深松鏟對土壤擾動程度最大,地表隆起、開裂和縱向起伏程度最強,此時土壤隆起高度、隆起寬度、橫剖面擾動面積、平均土塊徑均達到最大值,分別為16.3 cm、42.2 cm、0.030 5 m2、28.77 cm,耕作比阻則相對較低,為62 kN/m2;然而在耕深超過20 cm之后,深松鏟的土壤擾動效果顯著降低,地表隆起、開裂和縱向起伏程度減弱,土壤隆起高度、隆起寬度、橫剖面擾動面積、平均土塊徑數(shù)值均明顯減小,耕深30 cm時分別降至10.3 cm、31.2 cm、0.026 8 m2、19.12 cm,相比耕深20 cm降幅分別為36.8%、26.1%、12.1%和33.54%,而此時耕作比阻急劇增大至195 kN/m2,相比耕深20 cm增幅高達214.5%。因此,在水稻土條件下,耕作深度20 cm時能夠獲得最佳的土壤擾動、土壤結(jié)構(gòu)及耕作能量的綜合效應(yīng)。

    水稻土; 深松; 土壤擾動; 耕作阻力

    引言

    合理的耕作系統(tǒng)能夠改善耕層土壤結(jié)構(gòu)、調(diào)節(jié)土壤三相比例、協(xié)調(diào)土壤的水肥氣熱關(guān)系并為作物生長發(fā)育創(chuàng)造良好的環(huán)境和條件[1-7]。機械耕作引起的土壤物理狀態(tài)變化及其效應(yīng)的定量是耕作系統(tǒng)機理研究的重要內(nèi)容[8],也是優(yōu)化耕作部件設(shè)計、選擇最佳耕作方式、指導(dǎo)田間耕作管理的重要理論依據(jù)[9-10]。

    機械化深松是目前廣泛應(yīng)用的保護性耕作方式,該耕作系統(tǒng)利用深松鏟疏松土壤、加深耕層而不翻轉(zhuǎn)土壤、降低土壤容重、提高蓄水能力,促進作物的根系生長并提高作物產(chǎn)量[11-12]。深松鏟與土壤的作用機理屬于窄齒類觸土部件的力學(xué)范疇,近年來圍繞深松鏟的幾何設(shè)計、仿生設(shè)計研究形成了大量的成果[13-18]。在深松作業(yè)的機理研究方面多使用2類評價指標:深松部件的阻力特性和耕后土壤結(jié)構(gòu)(包括宏觀松土結(jié)構(gòu)和土壤松碎狀況)[19-21],分別涉及深松作業(yè)的不同方面,如土壤的宏觀擾動輪廓和阻力的關(guān)系[22]、鏟寬及作業(yè)參數(shù)(速度、耕深等)對耕作阻力的影響[23]、基于土壤擾動質(zhì)量的深松機具評價指標(土壤硬度變化系數(shù)、土壤體積膨松系數(shù)、單位松土帶寬度耕作阻力系數(shù)和土壤相互擾動系數(shù)等)[24]、振動型深松機的擾動剖面特征[15]、以土壤隆起及開裂等指標進行的不同鏟型深松機碎土質(zhì)量評價[20]等。

    然而,不同的土壤條件引起的深松鏟受力及耕后土壤擾動狀態(tài)完全不同[21,25],進而決定了某一個土壤條件的研究成果并不能普遍應(yīng)用于其他土壤條件。另外,多數(shù)研究工作是在室內(nèi)土槽的重塑土條件下完成的[22-24],測試過程無法體現(xiàn)田間真實土壤的結(jié)構(gòu)性及特有的約束條件。因此,針對具體的土壤類型及輪作制度的深松研究是完善保護性農(nóng)業(yè)理論和實踐的一個重要方面。

    稻麥輪作制是廣泛分布于長江中下游的高產(chǎn)種植制度,覆蓋面積高達1 300萬hm2[26],周年理論產(chǎn)量高達21 t/hm2[27]。然而該輪作制的周年干濕交替、水旱輪替耕作、稻麥根系的低效穿插及不良水肥運籌等導(dǎo)致水稻土特殊的結(jié)構(gòu)狀態(tài)和耕作力學(xué)屬性。研究發(fā)現(xiàn)長期以旋耕代替深耕造成水稻土耕層變淺、土壤硬化嚴重、孔隙率低、機械阻力增大等問題[7,28-33],雖然水稻土的機械深松能夠有效提升土壤物理質(zhì)量,但濕粘型水稻土的深松作業(yè)十分困難[34],這也是南方水田制度難于開展深松推廣的根本原因之一。

    目前針對水稻土的機械化深松研究已開始探討水稻土的機械深松作業(yè)的適應(yīng)性等問題[35]。然而水稻土的機械化深松過程的土壤結(jié)構(gòu)變化、耕層構(gòu)造及擾動剖面響應(yīng)狀況、深松的能效、適宜深松深度、最佳深松作業(yè)參數(shù)等指標都缺乏相應(yīng)的理論依據(jù)。因此,本文擬利用最新研發(fā)的土壤耕作原位綜合測試平臺[36]對水稻土的機械化深松進行精準控制試驗,探討水稻土在機械深松過程的土壤擾動、土壤結(jié)構(gòu)及能效特征,為水稻土條件的合理深松機械設(shè)計及耕作技術(shù)提供參考依據(jù)。

    1 材料與方法

    1.1 試驗地概況

    田間試驗于2015年12月水稻收獲后進行,試驗地點位于南京市六合區(qū)八百橋試驗田(118°59′E,31°98′N),地勢平坦,土壤類型為壤質(zhì)粘土,該地區(qū)常年稻麥輪作。土壤pH值7.6,土壤砂粒、壤粒、黏粒、有機物質(zhì)量分數(shù)分別為24.06%、40.39%、35.55%和2.27%,液限41.94%,塑限25.18%。0~50 cm土層土壤容重、含水率、總孔隙度見表1。試驗前進行人工地表清茬。

    表1 0~50 cm土層基礎(chǔ)物理參數(shù)
    Tab.1 Basic physical properties in 0~50 cm soil layer

    土層深度/cm容重/(g·cm-3)含水率/%總孔隙度/%0~51.2535.6952.665~101.2634.8552.4110~151.3233.4050.1115~201.5325.4342.3420~251.5623.6841.1125~301.5624.8941.2430~351.5624.4541.0735~401.5624.1341.0040~451.5924.0440.1145~501.6423.3138.16

    1.2 試驗方法

    所用深松鏟類型為鑿形深松鏟,采用單因素試驗方法,設(shè)置5種耕作深度:10、15、20、25、30 cm,牽引速度均控制在0.1 m/s,使用單鏟進行耕作試驗。

    圖1 土壤耕作原位綜合測試平臺結(jié)構(gòu)圖Fig.1 Schematic of in-situ experimental platform1.牽引電動機 2.控制柜 3.鏈條 4.信號拖鏈 5.臺車機構(gòu)6.升降電動機 7.臺車導(dǎo)軌 8.電源拖鏈 9.信號終端計算機10.移動導(dǎo)軌

    為保證田間試驗的可靠性,在最新設(shè)計的土壤耕作原位綜合測試平臺[36](圖1)上進行深松耕作試驗。試驗臺架基本參數(shù)為:長度8 m,寬度1.8 m。耕作時由功率為4 kW牽引電動機提供動力,牽引測試臺車(圖2)整體向前、向后移動??刂乒窨烧{(diào)節(jié)牽引電動機轉(zhuǎn)速和轉(zhuǎn)向,從而調(diào)節(jié)臺車移動速度和移動方向。深松鏟懸掛于臺車下方的升降懸架上,升降懸架與升降絲桿連接,由升降電動機帶動絲桿上下移動,從而達到精準調(diào)節(jié)耕深的目的。臺車前、后各裝2個拉壓力傳感器,工作時,數(shù)據(jù)采集卡采集傳感器信號,通過信號線傳輸至終端計算機內(nèi),可實時動態(tài)監(jiān)測并記錄耕作阻力數(shù)據(jù)。臺架4個柱腳下方鋪設(shè)導(dǎo)軌,4個0.75 kW的升降電動機控制柱腳內(nèi)部萬向輪伸出或收起,從而實現(xiàn)臺架在田間橫向和縱向移動。

    圖2 臺車結(jié)構(gòu)圖Fig.2 Schematic of moving car1.升降絲桿 2.升降電動機 3.移動懸架 4.深松鏟 5.升降懸架 6.拉力傳感器 7.信號采集柜 8.導(dǎo)向柱

    每次測試的有效深松行程為5 m,去除前后各1 m的非穩(wěn)定工作區(qū),取中間3 m為試驗數(shù)據(jù)采集和土壤擾動質(zhì)量測試區(qū)。每次深松后試驗臺橫向移動2 m再進行下一次深松作業(yè),以防止擾動區(qū)域相交或試驗臺不良擾動對下一次測試的破壞。

    1.3 測定項目與方法

    1.3.1 耕作阻力及比阻

    試驗臺內(nèi)部臺車前、后各安裝2個LKL-101/1T型拉壓力傳感器,分別編號:前1(f1)、前2(f2)、后1(r1)、后2(r2),耕作阻力F=Ff1+Ff2-(Fr1+Fr2),利用Excel 2010軟件進行信號處理,得出原始耕作阻力波動曲線(圖3a),對原始數(shù)據(jù)進行濾波處理后(圖3b),截取波動穩(wěn)定段并求其平均值,作為相應(yīng)耕作深度下的耕作阻力,并由此計算耕作比阻(耕作阻力與擾動截面積比值)[37]。

    圖3 耕作阻力波動曲線Fig.3 Wave curves of tillage iraction

    1.3.2 土壤宏觀擾動輪廓

    根據(jù)田間試驗觀測到的具體情況,并參考前人研究方法[20,22-23,38],將深松后土壤宏觀擾動輪廓分為:橫剖面、俯視平面、縱向地表起伏,繪制簡易輪廓圖(圖4a、圖5),并規(guī)定適當?shù)膸缀螀?shù),以此來反映不同耕作深度下的土壤宏觀擾動輪廓變化規(guī)律。

    圖4 橫剖面擾動輪廓Fig.4 Disturbed profile of cross-section

    圖5 俯視平面擾動輪廓Fig.5 Disturbed profile from plan view

    圖6 土壤輪廓測量儀結(jié)構(gòu)圖Fig.6 Schematic of soil profile meter1.底架 2.水平移動架 3.讀數(shù)面板 4.測量鋼纖 5.抬針柵條

    圖7 縱向地表起伏測量Fig.7 Test of earth’s surface smooth in lengthways direction

    利用卷尺和自制的土壤輪廓測量儀(圖6)測量土壤宏觀擾動輪廓。土壤輪廓測量儀主要參數(shù)為:長3 m,寬1.2 m,鋼釬間隔10 mm,鋼釬總數(shù)100個,測量高度范圍為-450~400 mm。深松耕作后,土壤輪廓測量儀鋼釬沿深松溝橫向放置(圖4b),測量土壤隆起輪廓和對應(yīng)的擾動邊界輪廓,兩者合并組成土壤剖面擾動輪廓;利用卷尺手工測量俯視平面的地表裂紋以及擾動土垡取出后的俯視剖面輪廓參數(shù),每項參數(shù)取5次重復(fù);土壤輪廓測量儀鋼釬縱向放置(圖7),測量縱向地表起伏情況,為防止鋼釬全部落入隆起裂縫中,鋼釬排列向深松溝兩側(cè)偏離5~10 cm,總測量長度2 m。

    1.3.3 土壤破碎體尺度分布

    每次5 m深松行程完成后,取數(shù)據(jù)測量區(qū)域內(nèi)的2 m進行土壤取樣。沿著裂紋方向用雙手緩慢取出大尺度土垡,平放于海綿墊上,利用卷尺測量土垡直徑。然后收集深松溝內(nèi)的小尺度土壤破碎體,進行篩分,篩孔直徑分別為64、32、16、8 mm。對留在64 mm孔篩內(nèi)的土塊再次手工測量土塊直徑,其他孔篩內(nèi)土壤以孔篩分級區(qū)間上、下限的算術(shù)平均值作為其土塊直徑[39],依次為48、24、12 mm,0~8 mm土壤粒徑默認為4 mm。所有土塊自然風(fēng)干3 d后稱量,計算平均土塊徑[40]和質(zhì)量分形維數(shù)[41],公式為

    (1)

    式中D——平均土塊徑di——土塊粒徑mi——對應(yīng)粒徑的土塊質(zhì)量

    (2)

    式中dj——篩分孔徑dmax——最大土塊粒徑M(δ

    2 結(jié)果與分析

    2.1 耕作阻力及比阻

    耕作阻力在一定程度上反映了實際耕作中的拖拉機功率消耗情況[42]。圖8a展示了耕作阻力與耕作深度的擬合關(guān)系曲線,回歸函數(shù)與試驗數(shù)據(jù)擬合的決定系數(shù)R2=0.991 8,擬合可靠。耕作阻力y與耕作深度x滿足二階函數(shù)擬合關(guān)系:y=0.009 4x2-0.16x+1.503 9。從圖中可以看出,耕作阻力隨著耕作深度增大而增大,并且耕作深度越大,增長速度越快,這與IBRAHMI等[43]、齊關(guān)宇等[44]的研究結(jié)果相一致。以上結(jié)果與水稻土土層性質(zhì)變化有關(guān)(表1),即隨著耕作深度的增加,土壤容重逐漸增大,孔隙度逐漸減小,土壤愈加緊實,且土壤擾動量增加,從而工作阻力迅速增大[44-45]。

    對于耕作系統(tǒng)的優(yōu)化來說,僅知道耕作阻力是一方面,更重要的是耕作比阻信息(耕作阻力與擾動面積比值)[37],耕作比阻能更全面地反映機具耕作效率[9]。圖8b展示了耕作比阻隨著耕作深度的變化趨勢。從圖中可以看出,耕作比阻隨著耕作深度變化缺乏規(guī)律性,當耕深為10、15、20 cm時,耕作比阻呈平緩的線性增長趨勢,比阻在50 kN/m2水平。然而當耕深增大到25、30 cm時,耕作比阻急劇增大,分別達192、195 kN/m2,相比耕深20 cm增幅依次為209.7%和214.5%。水稻土條件下深松耕作比阻呈現(xiàn)這一變化特征的原因是,當耕深不超過20 cm時,土壤容重較低、孔隙度較大、土壤較為松散、易于松動,耕作阻力隨耕深增長的同時土壤擾動面積S也成比例增長(表2),因此耕作比阻穩(wěn)定在相對較低水平;但當耕深超過20 cm后,土壤容重增大、孔隙度減小、土壤變得緊實、不易于松動,耕作阻力隨耕深大幅增長,但土壤擾動面積S卻沒有相應(yīng)比例地增長,甚至?xí)霈F(xiàn)下降現(xiàn)象,因此造成耕作比阻急劇增大,耕作能效顯著降低,無法保證深松作業(yè)的經(jīng)濟性。因此,在只考慮能量投入和土壤擾動程度的情況下,本文水稻土條件下的深松耕作深度應(yīng)該控制在20 cm左右。

    圖8 不同耕深下耕作阻力和耕作比阻的變化曲線Fig.8 Changing curves of tillage traction and specific traction at different tillage depths

    耕作深度H/cm壟起寬度W/cm裂縫寬度W1/cm擾動寬度W2/cm溝槽寬度W3/cm壟起高度H1/cm臨界深度H2/cm擾動面積S/m21026.49.026.311.70.01681535.013.232.715.30.02722042.221.929.12.616.315.20.03052537.410.122.92.511.314.80.01753031.213.724.12.510.314.80.0268

    2.2 土壤宏觀擾動輪廓

    2.2.1 橫剖面擾動輪廓

    以往研究表明,深松鏟耕作時存在擾動臨界深度[46],當耕深小于臨界深度時,土壤向前、向上、向兩側(cè)移動,形成新月形破壞區(qū);當耕深超過臨界深度時,臨界深度下方土壤只向前、向兩側(cè)移動,由新月形破壞區(qū)轉(zhuǎn)變?yōu)閭?cè)向破壞區(qū),此深度下方發(fā)生更多的是鏟具對土壤的擠壓,而不是有效松土(圖9)[45, 47-49]。本文研究也同樣發(fā)現(xiàn)這一規(guī)律,由圖10所示的不同耕作深度下的土壤橫剖面擾動輪廓可知,耕深為20、25、30 cm時,橫剖面擾動輪廓存在擾動臨界深度,大致位于15 cm深度處。深松鏟耕作過程中,臨界深度以上土壤被擾動且向兩側(cè)抬起,中間形成裂縫,地表以上呈“山峰狀”的隆起形態(tài),擾動邊界呈“上寬下窄”的扇形;而臨界深度下方土壤則受深松鏟側(cè)向擠壓作用,形成寬度與鏟寬大致相等的矩形溝槽。因此當耕深小于或接近臨界深度時,擾動邊界輪廓整體呈“扇形”狀;當耕深大于臨界深度時,擾動邊界輪廓則整體呈“漏斗”狀。

    圖9 深松土壤擾動示意圖Fig.9 Schematic of soil failure produced by subsoiler tines

    圖10 不同耕深下土壤橫剖面擾動輪廓Fig.10 Profiles cross-sections of soil at different tillage depths

    深松鏟入土深度的變化造成了整體剖面輪廓形態(tài)的變化,同時也造成了隆起參數(shù)和擾動邊界參數(shù)的變化。表2展示了深松后橫剖面擾動輪廓各項參數(shù)情況,從表中數(shù)據(jù)可以看出,地表以上的壟起高度H1、壟起寬度W、裂縫寬度W1均隨著耕深的增大整體呈“先增大后減小”的趨勢,均在耕深20 cm時達最大值;地表處的擾動寬度W2于耕深15 cm時最大,耕深20 cm次之,但差別不大;地表下方的土壤擾動面積S隨著耕深增大同樣呈現(xiàn)“先增大后減小”的趨勢,于耕深20 cm時達最大值0.030 5 m2,耕深25 cm時擾動截面積相比耕深20 cm不升反降,只有0.017 5 m2。

    以上分析表明,水稻土條件下深松耕作深度對深松鏟的松土、抬土能力有著較大影響。當耕深超過臨界深度(15 cm左右)時,土壤破壞失效模式發(fā)生改變,由臨界深度上方的新月形破壞失效轉(zhuǎn)變成臨界深度下方的側(cè)向破壞失效。這一臨界深度位置可能與水稻土犁底層的位置有關(guān)。由表1數(shù)據(jù)可知,犁底層同樣處于15 cm左右深度處,此深度附近土壤物理性質(zhì)發(fā)生劇烈變化,15~20 cm土層相比10~15 cm土層容重增加15.9%、孔隙度降低15.5%、含水率降低23.9%,犁底層下方土壤相比上方耕層土壤更加緊實,難以被松動、抬升,只能受到深松鏟的側(cè)向擠壓破壞??傮w上,深松鏟的整體土壤擾動能力在耕深20 cm時達到最大,此時地表隆起程度和土壤擾動量均達到最大值。但耕深超過20 cm后,深松鏟更多部分處于臨界深度下方,形成對土壤側(cè)向擠壓而不是有效松土,整體土壤擾動能力減弱,地表隆起程度下降,土壤擾動量減少。因此在本文水稻土條件下,為獲得最大土壤擾動效果,深松鏟耕作深度應(yīng)在20 cm左右。

    2.2.2 俯視平面裂紋及剖面擾動輪廓

    深松耕作后,地表形成裂紋,沿前進方向左右交替分布(圖5a),這與GODWIN等[50]研究結(jié)果相一致。將擾動土垡取出后,俯視剖面兩側(cè)的擾動邊界呈平行的曲線波動形態(tài)(圖5b)。而在LISOWSKI等[51]的土槽試驗中,擾動土壤取出后深松鏟兩側(cè)的擾動邊界呈平行的直線形態(tài)。因此,在高含水率的水稻土條件下,土壤團粒間粘性較大,深松鏟對兩側(cè)的土壤難以進行同步的擾動破壞,而是沿前進方向左右交替依次進行,沿前進方向形成獨特的曲線波動的擾動邊界形態(tài)。

    表3展示了俯視平面裂紋及剖面擾動輪廓參數(shù)。裂紋長度a和裂紋寬度b反映了深松鏟前進過程中的表層土壤破裂程度,裂紋間隔c和波峰間隔d反映了深松鏟完成一次左右交替擾動的間隔長度。從表中數(shù)據(jù)可以看出,深松鏟入土深度對地表破裂程度和擾動間隔均存在影響。隨著耕作深度的增大,裂紋長度a和裂紋寬度b均呈現(xiàn)“先增大后減小”的趨勢,均于耕深20 cm時達到最大值,耕深15 cm次之,耕深超過20 cm后均大幅下降;裂紋間隔c和波峰間隔d均在耕深15 cm時達最大值,耕深20 cm次之,且耕深超過20 cm后同樣出現(xiàn)大幅下降現(xiàn)象。

    表3 俯視平面裂紋及剖面輪廓參數(shù)Tab.3 Properties of cracks and profiles from plan view cm

    以上分析表明,當耕深不超過20 cm時,深松鏟對土壤的擾動能力隨入土深度增大而逐漸增大,左右交替擾動間隔長,地表破裂程度大,擾動充分。然而,當耕深超過20 cm后,深松鏟對土壤的擾動能力大幅下降,左右交替擾動間隔縮短,地表破裂程度下降,擾動不充分。

    2.2.3 縱向地表起伏

    由前文得知,深松后地表土壤隆起,且在前進方向上深松鏟對土壤的擾動是左右交替依次進行,因此縱向的土垡壟起排列呈現(xiàn)一定的起伏波動狀態(tài)(圖11)。從圖中可以看出,耕深15、20 cm時,地表高低起伏劇烈,波動間隔大。但耕深達到25、30 cm時,地表高低起伏程度急劇減弱,波動間隔有所減小。總體上,隨著耕作深度的增加,地表起伏劇烈性也大致呈“先增大后減小”趨勢,耕深20 cm時達最大程度,耕深超過20 cm后則急劇下降。

    以上結(jié)果與橫剖面擾動輪廓、俯視平面裂紋及剖面擾動輪廓的研究結(jié)果相吻合。在土層20 cm深度內(nèi),土壤容重相對較低,孔隙度較大,含水率較高,土壤松散且易于松動、抬升,深松鏟入土深度的增加導(dǎo)致土壤擾動量的增加,地表破裂程度加劇,土壤隆起增高,深松鏟對土壤擾動能力增強,從而導(dǎo)致縱向地表起伏劇烈性增大。然而,當耕作深度超過20 cm后,深松鏟更多部分處于臨界深度下方緊實且不易于松動、抬升的土壤中,深松鏟對土壤擾動能力有所減弱,土壤擾動量減少,地表破裂程度減弱,土壤隆起降低,從而導(dǎo)致縱向地表起伏程度較小。

    圖11 不同耕深下縱向地表起伏情況Fig.11 Earth’s surface smooth of different tillage depth in lengthways direction

    2.3 土壤破碎體尺度分布

    對土壤破碎體(或土塊)進行尺度分析,是土壤宏觀力學(xué)結(jié)構(gòu)研究的基本內(nèi)容之一[52]。平均土塊徑從整體上反映了耕作后的土塊大小[40],質(zhì)量分形維數(shù)則是對土壤破碎體的定量描述[53]。機械耕作后,質(zhì)量分形維數(shù)越大,土壤破碎程度就越高,平均土塊徑越小;質(zhì)量分形維數(shù)越小,土壤破碎程度就越低,平均土塊徑越大[54]。本文水稻土深松耕作后,土壤破碎體粒徑主要分布在10~40 cm區(qū)間(質(zhì)量分數(shù)大于50%)(圖12)。隨著耕作深度的增加,粒徑大于20 cm的土壤破碎體質(zhì)量分數(shù)先增大后減小,耕深20 cm時達最大值78.37%,耕深30 cm時最小,為38.80%;平均土塊徑先增大后減小(表4),耕深20 cm時達最大值28.77 cm,耕深30 cm時達最小值19.12 cm,降低33.54%;質(zhì)量分形維數(shù)于耕深20 cm時達最小值1.75,耕深10 cm時達最大值1.97。

    圖12 土壤破碎體徑級分布Fig.12 Relative size distribution of soil clods by weight

    以上結(jié)果表明,在土層20 cm深度內(nèi),隨著耕作深度的增加,深松擾動后大尺度土壤破碎體比例上升,平均土塊徑增大,質(zhì)量分形維數(shù)降低,整體破碎程度下降;但耕深超過20 cm后,大尺度土壤破碎體比例下降,平均土塊徑減小,質(zhì)量分形維數(shù)增大,整體破碎程度升高。這是因為,在水稻土土層20 cm深度內(nèi),耕作深度的增加導(dǎo)致了土壤擾動量的增加,更多土壤被松動、抬升,但由于水稻土的濕粘性,被擾動土壤不易破碎,仍能保持較大的土垡結(jié)構(gòu),因此深松后的土壤破碎體尺度逐漸增大,整體破碎程度下降;當耕深超過20 cm后,深松鏟更多部分處于臨界深度(15 cm左右)下方,形成對土壤側(cè)向擠壓而不是有效松土,土壤擾動量減少,被松動、抬升起的土壤也相應(yīng)減少,土垡結(jié)構(gòu)尺度下降,因此深松后形成的土壤破碎體尺度較小,整體破碎程度較高。

    表4 不同耕深下的平均土塊徑和分形維數(shù)
    Tab.4 Mean diameters of soilDand mass fractal dimensionDm

    耕作深度H/cm1015202530平均土塊徑D/cm20.9528.2228.7719.3219.12質(zhì)量分形維數(shù)Dm1.971.841.751.911.79

    3 結(jié)論

    使用土壤耕作原位綜合測試平臺針對南方稻麥輪作制開展機械深松的定量化研究表明,水稻土條件的深松作業(yè)對土壤擾動、土壤結(jié)構(gòu)及耕作用能的影響特征如下:

    (1)深松耕作阻力隨耕作深度的增大而增大,且兩者呈二階函數(shù)擬合關(guān)系:y=0.009 4x2-0.16x+1.503 9,決定系數(shù)R2=0.991 8,擬合關(guān)系可靠。

    (2) 深松鏟對土壤的擾動程度隨著耕作深度的變化而變化,整體呈“先增大后減小”的趨勢。在土層深度20 cm內(nèi),隨著耕作深度的增加,土壤擾動量逐漸增加,破碎體尺度增大,地表隆起、開裂以及縱向起伏程度加劇,整體擾動程度增大;當耕深超過20 cm后,深松鏟更多部分處于臨界深度(15 cm左右深度處)下方,土壤擾動量變少,破碎體尺度減小,地表隆起、開裂以及縱向起伏程度減弱,整體擾動程度下降。

    (3) 耕深20 cm時,土壤擾動量最大,耕作比阻處于相對較低水平。但當耕深超過20 cm后,土壤擾動量減少,但耕作阻力仍大幅增長,耕作比阻急劇增大,能效大大降低。因此,在僅考慮耕作能效和土壤擾動程度的前提下,水稻土條件下深松最佳耕作深度在20 cm左右,但此時土壤結(jié)構(gòu)尺度大,破碎程度低,地表平整性較差,因此后期應(yīng)配合旋耕或耙地等整地措施,以便于后期作物播種。

    1 張銀平,杜瑞成,刁培松,等. 機械化生態(tài)沃土耕作模式提高土壤質(zhì)量及作物產(chǎn)量[J]. 農(nóng)業(yè)工程學(xué)報, 2015,31(7): 33-38. ZHANG Yinping, DU Ruicheng, DIAO Peisong, et al. Mechanical and ecological tillage pattern improving soil quality and crop yields[J]. Transactions of the CSAE, 2015, 31(7): 33-38. (in Chinese)

    2 聶良鵬,郭利偉,牛海燕,等. 輪耕對小麥-玉米兩熟農(nóng)田耕層構(gòu)造及作物產(chǎn)量與品質(zhì)的影響[J]. 作物學(xué)報, 2015,41(3): 468-478.http:∥zwxb.chinacrops.org/CN/article/searchArticleResult.DOI: 10.3724/SP.J.1006.2015.00468. NIE Liangpeng, GUO Liwei, NIU Haiyan, et al. Effects of rotational tillage on tilth soil structure and crop yield and quality in maize-wheat cropping system[J]. Acta Agronomica Sinica, 2015, 41(3): 468-478. (in Chinese)

    3 王玉玲,李軍,柏煒霞. 輪耕體系對黃土臺塬麥玉輪作土壤生產(chǎn)性能的影響[J/OL]. 農(nóng)業(yè)工程學(xué)報, 2015,31(1): 107-116.http:∥www.tcsae.org/nygcxb/ch/reader/view_abstract.aspx?flag=1&file_no=20150116&journal_id=nygcxb. DOI:10.3969/j.issn.1002-6819.2015.01.016. WANG Yuling, LI Jun, BAI Weixia. Effects of rotational tillage systems on soil production performance in wheat-maize rotation field in Loess Platform region of China[J/OL]. Transactions of the CSAE, 2015, 31(1): 107-116. (in Chinese)

    4 鄭洪兵,鄭金玉,羅洋,等. 農(nóng)田不同耕層構(gòu)造對玉米生長發(fā)育及產(chǎn)量的影響[J/OL]. 干旱地區(qū)農(nóng)業(yè)研究, 2015,33(5): 41-45.http:∥ghdqnyyj.ijournal.cn/ghdqnyyj/ch/reader/view_abstract.aspx?flag=1&file_no=20150508&journal_id=ghdqnyyj. DOI: 10.7606/j.issn.1000-7601.2015.05.08. ZHENG Hongbing, ZHENG Jinyu, LUO Yang, et al. Effects of different tillage layer structures on growth and yield of maize in cropland zone in Northeast of China[J/OL]. Agricultural Research in the Arid Areas, 2015, 33(5): 41-45. (in Chinese)

    5 張玉嬌,李軍,郭正,等. 渭北旱塬麥田保護性輪耕方式的產(chǎn)量和土壤水分效應(yīng)長周期模擬研究[J/OL]. 中國農(nóng)業(yè)科學(xué), 2015,48(14): 2730-2746.http:∥www.chinaagrisci.com/CN/article/searchArticleResultByKeyword.DOI:10.3864/j.issn.0578-1752.2015.14.005. ZHANG Yujiao, LI Jun, GUO Zheng, et al. Long-term simulation of winter wheat yield and soil water response to conservation tillage rotation in Weibei Highland[J/OL]. Scientia Agricultura Sinica, 2015, 48(14): 2730-2746. (in Chinese)

    6 鄭洪兵,鄭金玉,羅洋,等. 長期不同耕作方式下的土壤硬度變化特征[J/OL]. 農(nóng)業(yè)工程學(xué)報, 2015,31(9): 63-70. http:∥www.tcsae.org/nygcxb/ch/reader/view_abstract.aspx?flag=1&file_no=20150911&journal_id=nygcxb. DOI: 10.11975/j.issn.1002-6819.2015.09.011. ZHENG Hongbing, ZHENG Jinyu, LUO Yang, et al. Change characteristic of soil compaction of long-term different tillage methods in cropland[J/OL]. Transactions of the CSAE, 2015, 31(9): 63-70. (in Chinese)

    7 何進,李洪文,高煥文. 中國北方保護性耕作條件下深松效應(yīng)與經(jīng)濟效益研究[J]. 農(nóng)業(yè)工程學(xué)報, 2006,22(10): 62-67. HE Jin, LI Hongwen, GAO Huanwen. Subsoiling effect and economic benefit under conservation tillage mode in Northern China[J]. Transactions of the CSAE, 2006, 22(10): 62-67. (in Chinese)

    8 丁啟朔,丁為民. 現(xiàn)代土壤機械耕作研究的綜述[J]. 土壤通報, 2006,37(1): 149-153. DING Qishuo, DING Weimin. Modern soil tillage research: a review[J]. Chinese Journal of Soil Science, 2016, 37(1): 149-153. (in Chinese)

    9 GODWIN R J. A review of the effect of implement geometry on soil failure and implement forces[J]. Soil and Tillage Research, 2007, 97(2): 331-340.

    10 曾德超. 機械土壤動力學(xué)[M]. 北京: 北京科學(xué)技術(shù)出版社, 1995.

    11 李洪文,陳君達,李問盈. 保護性耕作條件下深松技術(shù)研究[J]. 農(nóng)業(yè)機械學(xué)報, 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 for Agricultural Machinery, 2000, 31(6): 42-45. (in Chinese)

    12 丁昆侖,HANN M J. 深松耕作對土壤水分物理特性及作物生長的影響[J]. 中國農(nóng)村水利水電, 1997(11): 13-16. DING Kunlun, HANN M J. Effect of deep loosen cultivation to soil moisture physical characteristics and crop growing[J]. China Rural Water and Hydropower, 1997(11): 13-16. (in Chinese)

    13 張金波,佟金,馬云海. 仿生減阻深松鏟設(shè)計與試驗[J/OL]. 農(nóng)業(yè)機械學(xué)報, 2014,45(4): 141-145. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20140422&flag=1.DOI: 10.6041/j.issn.1000-1298.2014.04.022. ZHANG Jinbo, TONG Jin, MA Yunhai. Design and experiment of bionic anti-drag subsoiler[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(4): 141-145. (in Chinese)

    14 張強,張璐,于海業(yè),等. 復(fù)合形態(tài)深松鏟耕作阻力有限元分析與試驗[J/OL]. 農(nóng)業(yè)機械學(xué)報, 2012,43(8): 61-65. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20120812&flag=1. DOI: 10.6041/j.issn.1000-1298.2012.08.012. ZHANG Qiang, ZHANG Lu, YU Haiye, et al. Finite element analysis and experiment of soil resistance of multiplex-modality subsoiler[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2012, 43(8): 61-65. (in Chinese)

    15 李霞,付俊峰,張東興,等. 基于振動減阻原理的深松機牽引阻力試驗[J]. 農(nóng)業(yè)工程學(xué)報, 2012,28(1): 32-36. http:∥www.tcsae.org/nygcxb/ch/reader/view_abstract.aspx?flag=1&file_no=20120107&journal_id=nygcxb. DOI: 10.3969/j.issn.1002-6819.2012.01.007. LI Xia, FU Junfeng, ZHANG Dongxing, et al. Study on traction resistance based on the experiment analysis of a vibration ripper[J]. Transactions of the CSAE, 2012, 28(1): 32-36. (in Chinese)

    16 辛麗麗,李傳琦,梁繼輝,等. 考慮分段土壤作用力的振動減阻分析[J/OL]. 農(nóng)業(yè)機械學(xué)報, 2014,45(2): 136-140. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20140223&flag=1. DOI: 10.6041/j.issn.1000-1298.2014.02.023. XIN Lili, LI Chuanqi, LIANG Jihui, et al. Vibrating drag reduction considering acting force of piecewise soil[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(2): 136-140. (in Chinese)

    17 李霞,張東興,王維新,等. 受迫振動深松機性能參數(shù)優(yōu)化與試驗[J/OL]. 農(nóng)業(yè)工程學(xué)報, 2015,31(21): 17-24. http:∥www.tcsae.org/nygcxb/ch/reader/view_abstract.aspx?flag=1&file_no=20152103&journal_id=nygcxb. DOI:10.11975/j.issn.1002-6819.2015.21.003. LI Xia, ZHANG Dongxing, WANG Weixin, et al. Performance parameter optimization and experiment of forced-vibration subsoiler[J/OL]. Transactions of the CSAE, 2015, 31(21): 17-24. (in Chinese)

    18 郭志軍,周志立,張毅,等. 土壤耕作部件仿生優(yōu)化設(shè)計研究[J]. 中國科學(xué)E輯:技術(shù)科學(xué), 2009,39(4): 720-728. GUO Zhijun, ZHOU Zhili, ZHANG Yi, et al. Bionic optimization research of soil cultivating component design[J]. Science in China Series E: Technological Sciences, 2009, 39(4): 720-728. (in Chinese)

    19 CONTE O, LEVIEN R, DEBIASI H, et al. Soil disturbance index as an indicator of seed drill efficiency in no-tillage agrosystems[J]. Soil and Tillage Research, 2011, 114(1): 37-42.

    20 丁永凱,丁健偉. 深松部件的工作性能分析及其應(yīng)用[J]. 農(nóng)機化研究, 1990(2): 41-47.

    21 郭志軍,佟金,周志立,等. 深松技術(shù)研究現(xiàn)狀與展望[J]. 農(nóng)業(yè)工程學(xué)報, 2001,17(6): 169-174. GUO Zhijun, TONG Jin, ZHOU Zhili, et al. Review of subsoiling techniques and their applications[J].Transactions of the CSAE, 2001, 17(6): 169-174. (in Chinese)

    22 郭志軍,周志立,張毅,等. 深松耕作土壤宏觀擾動輪廓分析[J]. 拖拉機與農(nóng)用運輸車, 2003(3): 27-30. GUO Zhijun, ZHOU Zhili, ZHANG Yi, et al. The analysis of the macroscopical disturbing profile of soil after subsoiling[J]. Tractor & Farm Transporter, 2003(3): 27-30. (in Chinese)

    23 MANUWA S I. Performance evaluation of tillage tines operating under different depths in a sandy clay loam soil[J]. Soil and Tillage Research, 2009, 103(2): 399-405.

    24 黃玉祥,杭程光,李偉,等. 深松作業(yè)效果試驗及評價方法研究[J]. 西北農(nóng)林科技大學(xué)學(xué)報:自然科學(xué)版, 2015,43(11): 228-234. HUANG Yuxiang, HANG Chengguang, LI Wei, et al. Subsoiling test and evaluation methodology study of tillage quality[J]. Journal of Northwest A&F University: Natural Science Edition, 2015, 43(11): 228-234. (in Chinese)

    25 王燕. 基于離散元法的深松鏟結(jié)構(gòu)與松土效果研究[D]. 長春:吉林農(nóng)業(yè)大學(xué), 2014. WANG Yan. Simulation analysis of structure and effect of the subsoiler based on DEM[D]. Changchun: Jilin Agricultural University, 2014. (in Chinese)

    26 吳成龍. 稻—麥輪作系統(tǒng)有機無機肥配施協(xié)同土壤氮素轉(zhuǎn)化的機制研究[D]. 南京:南京農(nóng)業(yè)大學(xué), 2010. WU Chenglong. Mechanisms on soil nitrogen transformation as affected by integrated use of inorganic and organic fertilizer in a rice-wheat cropping system[D]. Nanjing: Nanjing Agricultural University, 2010. (in Chinese)

    27 楊建昌,杜永,劉輝. 長江下游稻麥周年超高產(chǎn)栽培途徑與技術(shù)[J]. 中國農(nóng)業(yè)科學(xué), 2008,41(6): 1611-1621. YANG Jianchang, DU Yong, LIU Hui. Cultivation approaches and techniques for annual super-high-yielding of rice and wheat in the lower reaches of Yangtze River[J]. Scientia Agricultura Sinica, 2008, 41(6): 1611-1621. (in Chinese)

    28 TAYLOR J H. Reduction of traffic-induced soil compaction[J]. Soil & Tillage Research, 1992, 24(4): 301-302.

    29 金軻,蔡典雄,呂軍杰,等. 耕作對坡耕地水土流失和冬小麥產(chǎn)量的影響[J]. 水土保持學(xué)報, 2006,20(4): 1-5. JIN Ke, CAI Dianxiong, Lü Junjie, et al. Effects of tillage practices on erosion and winter wheat yield on sloping dryland[J]. Journal of Soil and Water Conservation, 2006, 20(4): 1-5. (in Chinese)

    30 雷友,曹國鑫,牛新勝,等. 土壤深耕對冬小麥根系在土壤剖面分布的影響[J]. 現(xiàn)代農(nóng)業(yè)科技, 2011(8): 272-273. LEI You, CAO Guoxin, NIU Xinsheng, et al. Influence of deep tillage on distribution of wheat root system in soil profile[J].Modern Agricultural Sciences and Technology, 2011(8): 272-273. (in Chinese)

    31 孫忠英,李寶筏. 農(nóng)業(yè)機器行走裝置對土壤壓實作用的研究[J]. 農(nóng)業(yè)機械學(xué)報, 1998,29(3): 173-174. SUN Zhongying, LI Baofa. Research on the effect of soil compaction caused by the device of running agricultural machinery[J]. Transactions of the Chinese Society for Agricultural Machinery, 1998, 29(3): 173-174. (in Chinese)

    32 王恩姮,柴亞凡,陳祥偉. 大機械作業(yè)對黑土區(qū)耕地土壤結(jié)構(gòu)性特征的影響[J]. 應(yīng)用生態(tài)學(xué)報, 2008,19(2): 351-356. WANG Enheng, CHAI Yafan, CHEN Xiangwei. Effects of heavy machinery operation the structural characters of cultivated soils in black soil region of Northeast China[J]. Chinese Journal of Applied Ecology, 2008, 19(2): 351-356. (in Chinese)

    33 王恩姮,陳祥偉. 大機械作業(yè)對黑土區(qū)耕地土壤三相比與速效養(yǎng)分的影響[J]. 水土保持學(xué)報, 2007,21(4): 98-102. WANG Enheng, CHEN Xiangwei. Effect of heavy machinery operation on soil three phases and available nutrient in phaeozem region[J]. Journal of Soil and Water Conservation, 2007, 21(4): 98-102. (in Chinese)

    34 張軍,丁啟朔,丁為民,等. 分層旋耕對水稻土耕層物理性狀的影響[J/OL]. 南京農(nóng)業(yè)大學(xué)學(xué)報, 2015,38(6): 1016-1022. http:∥nauxb.njau.edu.cn/oa/DArticle.aspx?type=view&id=201506022. DOI: 10.7685/j.issn.1000-2030.2015.06.022. ZHANG Jun, DING Qishuo, DING Weimin, et al.Effect of stratified rotary tillage on paddy soil physical properties[J/OL]. Journal of Nanjing Agricultural University, 2015, 38(6): 1016-1022. (in Chinese)

    35 孫少華,張文斌. 機械化深松技術(shù)及蘇州地區(qū)適應(yīng)性研究[J]. 安徽農(nóng)業(yè)科學(xué), 2012,40(20): 10701-10702. SUN Shaohua, ZHANG Wenbin. Research of subsoiling techniques and their adaptability in Suzhou[J]. Journal of Anhui Agricultural Sciences, 2012, 40(20): 10701-10702. (in Chinese)

    36 楊艷山,丁啟朔,丁為民,等. 田間原位綜合耕作試驗臺設(shè)計與應(yīng)用[J/OL]. 農(nóng)業(yè)機械學(xué)報, 2016,47(1): 68-74. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20160110&flag=1.DOI:10.6041/j.issn.1000-1298.2016.01.010. YANG Yanshan, DING Qishuo, DING Weimin, et al. Design and application of multi-purpose in-situ tillage tool testing platform[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(1): 68-74. (in Chinese)

    37 ARVIDSSON J, KELLER T, GUSTAFSSON K. Specific draught for mouldboard plough, chisel plough and disc harrow at different water contents[J]. Soil and Tillage Research, 2004, 79(2): 221-231.

    39 WATTS C W, DEXTER A R. Traffic and seasonal influences on the energy required for cultivation and on the subsequent tilth[J]. Soil & Tillage Research, 1994, 31(4): 303-322.

    40 辜松,劉立意,小松崎將一,等. 覆蓋作物對土壤物理特性及旋耕作業(yè)的影響分析[J]. 農(nóng)業(yè)機械學(xué)報, 2005,36(5): 41-45. GU Song, LIU Liyi, KOMATSUZAKI Masakazu, et al.Effect of cover crops on soil physical properties and rotary tillage[J]. Transactions of the Chinese Society for Agricultural Machinery, 2005, 36(5): 41-45. (in Chinese)

    41 楊金玲,李德成,張甘霖,等. 土壤顆粒粒徑分布質(zhì)量分形維數(shù)和體積分形維數(shù)的對比[J]. 土壤學(xué)報, 2008,45(3): 413-419. YANG Jinling, LI Decheng, ZHANG Ganlin, et al. Comparsion of mass and volume fractal dimensions of soil particle size distributions[J]. Acta Pedologica Sinica, 2008, 45(3): 413-419. (in Chinese)

    42 黃虎,王曉燕,李洪文,等. 固定道保護性耕作節(jié)能效果試驗研究[J]. 農(nóng)業(yè)工程學(xué)報, 2007,23(12): 140-143. HUANG Hu, WANG Xiaoyan, LI Hongwen, et al. Experimental investigation on energy saving of controlled traffic conservation tillage[J]. Transactions of the CSAE, 2007, 23(12): 140-143. (in Chinese)

    43 IBRAHMI A, BENTAHER H, HBAIEB M, et al. Study the effect of tool geometry and operational conditions on mouldboard plough forces and energy requirement: Part 1. Finite element simulation[J]. Computers and Electronics in Agriculture, 2015, 117: 258-267.

    44 齊關(guān)宇,劉林,趙艷忠,等. 深松鏟入土深度及鏟形對耕作阻力影響[J]. 農(nóng)機化研究, 2015(11): 161-165. QI Guanyu, LIU Lin, ZHAO Yanzhong, et al. Effects of subsoiler’s penetrating depth and spade shape on traction resistance[J]. Journal of Agricultural Mechanization Research, 2015(11): 161-165. (in Chinese)45 郭志軍,杜干,周志立,等. 土壤耕作部件宏觀觸土曲面減阻性能研究現(xiàn)狀分析[J]. 農(nóng)業(yè)機械學(xué)報, 2011,42(6): 47-52. GUO Zhijun, DU Gan, ZHOU Zhili, et al. Actuality analysis of resistance-reducing properties on soil cultivating components with different macroscopic soil-engaging surfaces[J]. Transactions of the Chinese Society for Agricultural Machinery, 2011, 42(6): 47-52. (in Chinese)

    46 SPOOR G, FRY R K. Soil disturbance generated by deep-working low rake angle narrow tines[J]. Journal of Agricultural Engineering Research, 1983, 28(3): 217-234.

    47 ALUKO O B, SEIG D A. An experimental investigation of the characteristics of and conditions for brittle fracture in two-dimensional soil cutting[J]. Soil and Tillage Research, 2000, 57(3): 143-157.

    48 GODWIN R J, SPOOR G. Soil failure with narrow tines[J]. Journal of Agricultural Engineering Research, 1977, 22(3): 213-228.

    49 MCKYES E, ALI O S. The cutting of soil by narrow blades[J]. Journal of Terramechanics, 1977, 14(2): 43-58.

    50 GODWIN R J, SPOOR G, LEEDS-HARRISON P. An experimental investigation into the force mechanics and resulting soil disturbance of mole ploughs[J]. Journal of Agricultural Engineering Research, 1981, 26(6): 477-497.

    51 LISOWSKI A, KLONOWSKI J, GREEN O, et al. Duckfoot tools connected with flexible and stiff tines: three components of resistances and soil disturbance[J]. Soil and Tillage Research, 2016, 158: 76-90.

    52 丁啟朔,丁為民,潘根興,等. 土壤宏觀力學(xué)結(jié)構(gòu)與精準耕作[J/OL]. 中國農(nóng)業(yè)科學(xué), 2012,45(1): 26-33. http:∥www.chinaagrisci.com/CN/article/searchArticleResultByKeyword.DOI:10.3864/j.issn.0578-1752.2012.01.004. DING Qishuo, DING Weimin, PAN Genxing, et al. Soil macro-structure and precision tillage[J/OL]. Scientia Agricultura Sinica, 2012, 45(1): 26-33. (in Chinese)

    53 孟鳳英,丁啟朔,鹿飛,等. 沖擊作用下粘性土壤破碎體的分形維數(shù)與影響因素[J]. 農(nóng)業(yè)機械學(xué)報, 2009,40(3): 108-111. MENG Fengying, DING Qishuo, LU Fei, et al. Fragmentation fractal dimensions of cohesive soil under impact and its influencing factors[J]. Transactions of the Chinese Society for Agricultural Machinery, 2009, 40(3): 108-111. (in Chinese)

    54 韓秋萍. 水稻土條件的旋耕機性能測試及小麥響應(yīng)研究[D]. 南京:南京農(nóng)業(yè)大學(xué), 2011. HAN Qiuping. Measurement of rotary tillor performance in paddy soil condition and wheat response research[D]. Nanjing: Nanjing Agricultural University, 2011. (in Chinese)

    Characteristics of Subsoiler Traction and Soil Disturbance in Paddy Soil

    DING Qishuo GE Shuangyang REN Jun LI Yinian HE Ruiyin

    (CollegeofEngineering,NanjingAgriculturalUniversity,Nanjing210031,China)

    Optimization of tillage system for paddy soil requires a precise quantification of soil disturbance, soil structure and tillage energy efficiency. An in-situ tillage test plateform was applied for precise control on the working parameters of a subsoiler, which was a chisel design and was evaluated in five depths (i.e., 10 cm, 15 cm, 20 cm, 25 cm and 30 cm). Traction, specific resistance, disturbed soil profile and soil fragment distribution were measured to evaluate the performance of the subsoiler. Results showed an increased traction with respect to the engaging depth, which satisfied a second-order equation. In the 20 cm tilling depth the degree of soil disturbance, soil surface humping, soil cracking and longitudinal humping were found the maximum. The bumping height, width, disturbed area of soil profile, mean weight diameter were 16.3 cm, 42.2 cm, 0.030 5 m2and 28.77 cm, respectively. Specific resistance in this depth was the minimum, i.e. 62 kN/m2. When surpassing 20 cm depth, the bumping height, width, disturbed soil profile area, mean weight diameter were significantly decreased, which were 10.3 cm, 31.2 cm, 0.026 8 m2and 19.12 cm, respectively, for the 30 cm depth. Compared with the 20 cm depth, these parameters were decreased by 36.8%,26.1%,12.1% and 33.54%, respectively. Meanwhile, the specific resistance was 195 kN/m2for the 30 cm depth which was increased by 214.5% compared with 20 cm depth. In considerating the overall effect of tillage-induced soil struction,soil tilth quality and energy use efficiency, tilling the paddy soil into 20 cm provided the best results.

    paddy soil; subsoiler; soil disturbance; tillage traction

    10.6041/j.issn.1000-1298.2017.01.007

    2016-06-01

    2016-06-25

    國家重點研發(fā)計劃項目(2016YFD0300900)和江蘇省農(nóng)機基金項目(201-051028)

    丁啟朔(1968—),男,教授,博士生導(dǎo)師,主要從事土壤耕作力學(xué)研究,E-mail: qsding@njau.edu.cn

    S222

    A

    1000-1298(2017)01-0047-10

    猜你喜歡
    耕深土塊耕作
    拖拉機多重模糊PID變論域耕深調(diào)節(jié)研究
    叩土籠
    賊的藏身術(shù)
    北方人(2021年24期)2021-02-26 08:52:06
    賊的藏身術(shù)
    賊的藏身術(shù)
    視野(2020年20期)2020-12-23 09:32:32
    基于卡爾曼濾波融合算法的深松耕深檢測裝置研究
    懸掛式深松機耕深監(jiān)測系統(tǒng)的設(shè)計與試驗
    價值工程(2019年31期)2019-12-04 04:11:05
    線性擬合與Kalman預(yù)測法修正耕深測量誤差
    耕作深度對紫色土坡地旋耕機耕作侵蝕的影響
    玉米保護性耕作的技術(shù)要領(lǐng)
    精品午夜福利在线看| 日日撸夜夜添| 我的女老师完整版在线观看| 欧美老熟妇乱子伦牲交| 国产亚洲av片在线观看秒播厂| 一级毛片电影观看| 国产爽快片一区二区三区| 成人二区视频| 色视频www国产| 成人高潮视频无遮挡免费网站| 亚洲久久久国产精品| 亚洲电影在线观看av| 中文天堂在线官网| 大话2 男鬼变身卡| 女性生殖器流出的白浆| 精品国产露脸久久av麻豆| 亚洲av欧美aⅴ国产| 国产亚洲5aaaaa淫片| 日本色播在线视频| av又黄又爽大尺度在线免费看| 一级毛片电影观看| 国产免费视频播放在线视频| 18禁在线无遮挡免费观看视频| 麻豆成人午夜福利视频| 午夜激情福利司机影院| 国产无遮挡羞羞视频在线观看| 精品99又大又爽又粗少妇毛片| 在线观看人妻少妇| 亚洲av中文字字幕乱码综合| 久久精品国产亚洲av涩爱| 人人妻人人澡人人爽人人夜夜| 精品熟女少妇av免费看| 亚洲精品日韩在线中文字幕| 欧美bdsm另类| 各种免费的搞黄视频| 国产高清国产精品国产三级 | 国产精品久久久久久精品电影小说 | 免费观看性生交大片5| 美女内射精品一级片tv| 免费观看的影片在线观看| av女优亚洲男人天堂| 天堂俺去俺来也www色官网| 久久精品国产鲁丝片午夜精品| 国产成人a区在线观看| 国产精品人妻久久久久久| 少妇高潮的动态图| 永久免费av网站大全| 黄片无遮挡物在线观看| 亚洲怡红院男人天堂| 日韩三级伦理在线观看| 偷拍熟女少妇极品色| 国产高清有码在线观看视频| 高清欧美精品videossex| 麻豆精品久久久久久蜜桃| 九九爱精品视频在线观看| 日日啪夜夜爽| 国产毛片在线视频| 国产精品久久久久久久电影| 午夜日本视频在线| 久久久久精品久久久久真实原创| 丰满人妻一区二区三区视频av| 国产精品人妻久久久影院| 久久精品熟女亚洲av麻豆精品| 纯流量卡能插随身wifi吗| 久久人人爽人人片av| 干丝袜人妻中文字幕| 性色av一级| 亚洲精品一区蜜桃| 少妇的逼水好多| 午夜精品国产一区二区电影| 尤物成人国产欧美一区二区三区| 久久婷婷青草| 高清视频免费观看一区二区| 国产白丝娇喘喷水9色精品| 人体艺术视频欧美日本| 亚洲欧美一区二区三区黑人 | 极品教师在线视频| 亚洲av男天堂| 国产成人精品久久久久久| 汤姆久久久久久久影院中文字幕| 日韩中字成人| 免费不卡的大黄色大毛片视频在线观看| 黑人猛操日本美女一级片| 亚洲自偷自拍三级| 久久青草综合色| 国内精品宾馆在线| 久久精品久久久久久久性| 精品亚洲成国产av| 亚洲精品一二三| 人体艺术视频欧美日本| 亚洲熟女精品中文字幕| 18禁裸乳无遮挡动漫免费视频| 99九九线精品视频在线观看视频| 亚洲欧美清纯卡通| 韩国高清视频一区二区三区| 观看免费一级毛片| 一个人看视频在线观看www免费| 日本-黄色视频高清免费观看| 欧美bdsm另类| 欧美一区二区亚洲| 成年人午夜在线观看视频| 国产91av在线免费观看| 一级毛片aaaaaa免费看小| 久久久久久久久久成人| 欧美日韩一区二区视频在线观看视频在线| 午夜福利网站1000一区二区三区| 亚洲欧美日韩无卡精品| 我的女老师完整版在线观看| 久久久久人妻精品一区果冻| 91久久精品国产一区二区成人| 麻豆成人午夜福利视频| 麻豆乱淫一区二区| 高清毛片免费看| 一级毛片aaaaaa免费看小| 亚洲中文av在线| 国产高潮美女av| 久久精品国产自在天天线| 久久久精品94久久精品| 女人十人毛片免费观看3o分钟| 国产精品久久久久久精品电影小说 | 亚洲人与动物交配视频| 久久精品国产鲁丝片午夜精品| 另类亚洲欧美激情| 欧美日韩一区二区视频在线观看视频在线| 在线观看免费日韩欧美大片 | 乱码一卡2卡4卡精品| 国产精品99久久久久久久久| 在线看a的网站| 免费高清在线观看视频在线观看| 国产精品国产三级国产av玫瑰| 中国国产av一级| 中文字幕精品免费在线观看视频 | 久久国产精品男人的天堂亚洲 | 国产一区二区三区综合在线观看 | 3wmmmm亚洲av在线观看| 黑人猛操日本美女一级片| 女人十人毛片免费观看3o分钟| 欧美精品亚洲一区二区| 九九爱精品视频在线观看| kizo精华| 欧美xxⅹ黑人| 欧美日韩精品成人综合77777| 亚洲欧洲日产国产| 六月丁香七月| 在线免费十八禁| 香蕉精品网在线| 99久久精品一区二区三区| 你懂的网址亚洲精品在线观看| 又大又黄又爽视频免费| 国产视频首页在线观看| 午夜老司机福利剧场| 建设人人有责人人尽责人人享有的 | 纯流量卡能插随身wifi吗| 久久国产精品男人的天堂亚洲 | av在线播放精品| 欧美成人a在线观看| 久久人人爽人人片av| 男女无遮挡免费网站观看| 九九在线视频观看精品| 一区二区三区精品91| 国产成人精品婷婷| 久久韩国三级中文字幕| 久久久色成人| 国产淫片久久久久久久久| 欧美xxxx性猛交bbbb| 亚洲av不卡在线观看| 亚洲图色成人| 熟女av电影| 亚洲一区二区三区欧美精品| 最近最新中文字幕免费大全7| 日本-黄色视频高清免费观看| 国精品久久久久久国模美| av播播在线观看一区| 亚洲经典国产精华液单| 国产欧美日韩精品一区二区| 国产成人免费无遮挡视频| 寂寞人妻少妇视频99o| 在线精品无人区一区二区三 | 精品人妻一区二区三区麻豆| 国产免费又黄又爽又色| 国产在线免费精品| av免费观看日本| 丰满迷人的少妇在线观看| 丝袜脚勾引网站| av线在线观看网站| 久久久久国产网址| 成人黄色视频免费在线看| 久久影院123| 亚洲综合精品二区| 熟妇人妻不卡中文字幕| 在线精品无人区一区二区三 | 久久久久网色| 久久精品国产鲁丝片午夜精品| 成人高潮视频无遮挡免费网站| 亚洲精品日韩在线中文字幕| 夫妻午夜视频| 免费观看的影片在线观看| 久久久久久久久久久丰满| 免费大片黄手机在线观看| 在线观看免费日韩欧美大片 | 国产一区有黄有色的免费视频| 女性生殖器流出的白浆| 美女主播在线视频| 99九九线精品视频在线观看视频| 最近2019中文字幕mv第一页| 国产日韩欧美亚洲二区| 大香蕉久久网| 在线观看一区二区三区| 一区二区三区四区激情视频| 一级黄片播放器| 一级爰片在线观看| 欧美最新免费一区二区三区| 一区在线观看完整版| 91在线精品国自产拍蜜月| 亚洲欧美日韩另类电影网站 | 插逼视频在线观看| 黄片无遮挡物在线观看| 激情 狠狠 欧美| 最近2019中文字幕mv第一页| 少妇的逼水好多| 午夜福利高清视频| 高清日韩中文字幕在线| 久久午夜福利片| 亚洲图色成人| 不卡视频在线观看欧美| 亚洲成人中文字幕在线播放| 观看美女的网站| 蜜桃亚洲精品一区二区三区| 日韩成人伦理影院| 日本与韩国留学比较| 我要看日韩黄色一级片| 久久精品久久久久久久性| 在线精品无人区一区二区三 | 在线免费观看不下载黄p国产| 日本欧美视频一区| 久久久久久久久久人人人人人人| 日本wwww免费看| 好男人视频免费观看在线| 精品人妻一区二区三区麻豆| 黑丝袜美女国产一区| 久久久久久伊人网av| 最后的刺客免费高清国语| 久久精品国产鲁丝片午夜精品| 三级经典国产精品| 日本色播在线视频| 国产精品爽爽va在线观看网站| 欧美xxxx性猛交bbbb| 欧美高清成人免费视频www| 日韩av不卡免费在线播放| 久久久国产一区二区| 国产精品嫩草影院av在线观看| 老女人水多毛片| 免费观看在线日韩| 五月玫瑰六月丁香| 婷婷色综合www| 边亲边吃奶的免费视频| 欧美极品一区二区三区四区| 小蜜桃在线观看免费完整版高清| 舔av片在线| 天堂俺去俺来也www色官网| 天堂中文最新版在线下载| 婷婷色麻豆天堂久久| 色综合色国产| 久久人人爽人人片av| 少妇人妻一区二区三区视频| 26uuu在线亚洲综合色| 搡老乐熟女国产| 最近的中文字幕免费完整| 亚洲精华国产精华液的使用体验| 一级毛片我不卡| av在线老鸭窝| 两个人的视频大全免费| 夫妻午夜视频| 青春草国产在线视频| 久久久久国产网址| 日韩不卡一区二区三区视频在线| 精品久久久久久久末码| 日日啪夜夜爽| 国产av一区二区精品久久 | 精品久久国产蜜桃| 国产精品99久久99久久久不卡 | 国产永久视频网站| 国产视频首页在线观看| 亚洲欧美日韩另类电影网站 | 亚洲av中文av极速乱| 九色成人免费人妻av| 国产精品99久久99久久久不卡 | 亚洲第一区二区三区不卡| 亚洲伊人久久精品综合| 性色avwww在线观看| 国产深夜福利视频在线观看| 久久久久久久久久久免费av| 成人国产av品久久久| 亚洲精品乱码久久久久久按摩| 欧美少妇被猛烈插入视频| 少妇人妻精品综合一区二区| 欧美日本视频| 欧美国产精品一级二级三级 | 好男人视频免费观看在线| 国产色爽女视频免费观看| 免费观看av网站的网址| 777米奇影视久久| 国产成人aa在线观看| 国产av一区二区精品久久 | 免费大片18禁| 中文字幕制服av| 国产在线免费精品| 国产精品秋霞免费鲁丝片| 午夜老司机福利剧场| av专区在线播放| 婷婷色麻豆天堂久久| 深爱激情五月婷婷| 久热久热在线精品观看| 亚洲人成网站在线观看播放| 欧美精品一区二区免费开放| 午夜福利视频精品| 丝袜喷水一区| 亚洲精品aⅴ在线观看| 国产精品女同一区二区软件| 亚洲国产成人一精品久久久| 大陆偷拍与自拍| 国产欧美另类精品又又久久亚洲欧美| 熟女人妻精品中文字幕| 久久久欧美国产精品| 久久国产乱子免费精品| 91久久精品国产一区二区成人| 乱码一卡2卡4卡精品| 亚洲精品aⅴ在线观看| 美女xxoo啪啪120秒动态图| 女的被弄到高潮叫床怎么办| 免费av不卡在线播放| 一级片'在线观看视频| 成年人午夜在线观看视频| 欧美人与善性xxx| 国内揄拍国产精品人妻在线| 国产亚洲精品久久久com| 精品一品国产午夜福利视频| 亚洲精品,欧美精品| 六月丁香七月| 成人综合一区亚洲| 亚洲一区二区三区欧美精品| 久久99热这里只频精品6学生| 老师上课跳d突然被开到最大视频| 91久久精品电影网| 国产欧美另类精品又又久久亚洲欧美| 秋霞在线观看毛片| 亚洲经典国产精华液单| xxx大片免费视频| 国产精品国产三级国产专区5o| 成人免费观看视频高清| 各种免费的搞黄视频| 国产黄色视频一区二区在线观看| 麻豆成人午夜福利视频| 亚洲人成网站在线播| 寂寞人妻少妇视频99o| 搡老乐熟女国产| 少妇精品久久久久久久| 国产乱人视频| 精品久久久久久久末码| 97在线人人人人妻| 国产乱人视频| 久久久久人妻精品一区果冻| 麻豆精品久久久久久蜜桃| 高清午夜精品一区二区三区| 亚洲精品,欧美精品| 亚洲成人手机| 青春草视频在线免费观看| 最黄视频免费看| 日本一二三区视频观看| 午夜免费男女啪啪视频观看| 精品国产三级普通话版| 麻豆乱淫一区二区| av福利片在线观看| 免费大片18禁| videos熟女内射| 91aial.com中文字幕在线观看| 国产免费一区二区三区四区乱码| 国产欧美日韩精品一区二区| 亚洲精品日韩av片在线观看| 久久国内精品自在自线图片| 日韩成人伦理影院| 97超视频在线观看视频| 青春草国产在线视频| 国产有黄有色有爽视频| 国产精品女同一区二区软件| 国产大屁股一区二区在线视频| 国产精品福利在线免费观看| 免费观看的影片在线观看| 中国美白少妇内射xxxbb| 亚洲av.av天堂| 嘟嘟电影网在线观看| 天堂中文最新版在线下载| 视频中文字幕在线观看| 国产日韩欧美在线精品| 少妇人妻一区二区三区视频| 热99国产精品久久久久久7| 高清视频免费观看一区二区| 色吧在线观看| 韩国高清视频一区二区三区| 欧美区成人在线视频| 99九九线精品视频在线观看视频| 国产 一区 欧美 日韩| 国产黄频视频在线观看| 久久午夜福利片| 国产精品99久久久久久久久| 国产精品国产三级国产av玫瑰| 国产精品国产三级国产专区5o| 波野结衣二区三区在线| 一边亲一边摸免费视频| 女的被弄到高潮叫床怎么办| 男女啪啪激烈高潮av片| av卡一久久| 国产精品人妻久久久久久| 五月玫瑰六月丁香| 免费黄网站久久成人精品| 精品人妻偷拍中文字幕| 国产免费视频播放在线视频| a级一级毛片免费在线观看| 3wmmmm亚洲av在线观看| 久久精品夜色国产| 女性生殖器流出的白浆| 国产高清国产精品国产三级 | 成人黄色视频免费在线看| av女优亚洲男人天堂| 亚洲色图av天堂| 国产精品精品国产色婷婷| 男人添女人高潮全过程视频| 搡女人真爽免费视频火全软件| 国产精品一二三区在线看| 国产免费福利视频在线观看| 久久国产亚洲av麻豆专区| 国产午夜精品一二区理论片| 久久久久久久久大av| 男人狂女人下面高潮的视频| 纵有疾风起免费观看全集完整版| 国产成人a∨麻豆精品| 亚洲精华国产精华液的使用体验| 国产成人免费观看mmmm| 如何舔出高潮| 成人毛片a级毛片在线播放| av卡一久久| 中文资源天堂在线| 日日撸夜夜添| 亚洲av不卡在线观看| 身体一侧抽搐| 王馨瑶露胸无遮挡在线观看| 国产在线视频一区二区| 一边亲一边摸免费视频| 狠狠精品人妻久久久久久综合| 日日啪夜夜撸| 日韩亚洲欧美综合| 三级经典国产精品| 女性被躁到高潮视频| 亚洲欧美日韩另类电影网站| 免费观看人在逋| 国产av国产精品国产| 国产精品麻豆人妻色哟哟久久| 男人操女人黄网站| 建设人人有责人人尽责人人享有的| 国产伦人伦偷精品视频| 少妇猛男粗大的猛烈进出视频| 欧美黄色片欧美黄色片| 精品久久蜜臀av无| 亚洲 国产 在线| 亚洲精品久久久久久婷婷小说| 欧美精品啪啪一区二区三区 | 亚洲专区中文字幕在线| 丝袜人妻中文字幕| 一本大道久久a久久精品| 波野结衣二区三区在线| 黄色a级毛片大全视频| 国产不卡av网站在线观看| 国产亚洲欧美精品永久| 女人精品久久久久毛片| 亚洲欧美日韩另类电影网站| 满18在线观看网站| 9191精品国产免费久久| 亚洲国产最新在线播放| 蜜桃在线观看..| 亚洲情色 制服丝袜| 国产淫语在线视频| 操美女的视频在线观看| 每晚都被弄得嗷嗷叫到高潮| 深夜精品福利| 精品一区在线观看国产| 波多野结衣av一区二区av| 午夜老司机福利片| 午夜视频精品福利| 青草久久国产| 一级黄片播放器| 亚洲国产欧美一区二区综合| 老司机深夜福利视频在线观看 | 国产男女超爽视频在线观看| 黄色a级毛片大全视频| 蜜桃在线观看..| 美女扒开内裤让男人捅视频| 久久久国产欧美日韩av| 国产极品粉嫩免费观看在线| 激情视频va一区二区三区| 日韩av在线免费看完整版不卡| 九色亚洲精品在线播放| 黄频高清免费视频| 午夜91福利影院| 大香蕉久久网| 日日夜夜操网爽| 91精品三级在线观看| av一本久久久久| 久热这里只有精品99| 欧美黄色淫秽网站| 成人国产一区最新在线观看 | 欧美日韩亚洲高清精品| 亚洲色图 男人天堂 中文字幕| 一边摸一边做爽爽视频免费| av网站免费在线观看视频| 国产精品熟女久久久久浪| 中文字幕人妻丝袜制服| 欧美乱码精品一区二区三区| 视频区欧美日本亚洲| 亚洲色图综合在线观看| 1024香蕉在线观看| 欧美少妇被猛烈插入视频| 成年动漫av网址| 免费久久久久久久精品成人欧美视频| 女人高潮潮喷娇喘18禁视频| 亚洲欧美成人综合另类久久久| 狠狠精品人妻久久久久久综合| 精品国产乱码久久久久久小说| 老汉色∧v一级毛片| 制服诱惑二区| 午夜两性在线视频| 麻豆av在线久日| 久久久久国产精品人妻一区二区| 人成视频在线观看免费观看| svipshipincom国产片| 久久亚洲精品不卡| 免费久久久久久久精品成人欧美视频| 久久国产精品男人的天堂亚洲| 爱豆传媒免费全集在线观看| 波多野结衣av一区二区av| 日韩大片免费观看网站| 无遮挡黄片免费观看| 精品欧美一区二区三区在线| 天天躁狠狠躁夜夜躁狠狠躁| 又粗又硬又长又爽又黄的视频| 免费少妇av软件| 欧美日韩亚洲国产一区二区在线观看 | 七月丁香在线播放| 亚洲图色成人| e午夜精品久久久久久久| 国产日韩欧美亚洲二区| 最新在线观看一区二区三区 | 大型av网站在线播放| 中文乱码字字幕精品一区二区三区| 99国产综合亚洲精品| 高清欧美精品videossex| 国产伦人伦偷精品视频| 日本欧美国产在线视频| 亚洲七黄色美女视频| 在线 av 中文字幕| 国产精品 国内视频| 午夜两性在线视频| 麻豆国产av国片精品| 久久人人97超碰香蕉20202| 操美女的视频在线观看| tube8黄色片| 国产精品免费大片| 91麻豆精品激情在线观看国产 | 纵有疾风起免费观看全集完整版| 色婷婷久久久亚洲欧美| 一级片免费观看大全| 黄网站色视频无遮挡免费观看| 国产精品 国内视频| 亚洲欧美日韩另类电影网站| 久久精品亚洲av国产电影网| 久久久久久人人人人人| 久久久精品国产亚洲av高清涩受| 亚洲伊人色综图| 国产精品亚洲av一区麻豆| 久久精品熟女亚洲av麻豆精品| 成人影院久久| 免费在线观看完整版高清| 亚洲精品av麻豆狂野| 色94色欧美一区二区| 99久久99久久久精品蜜桃| 日本a在线网址| 欧美日韩精品网址| 国产精品偷伦视频观看了| 婷婷色av中文字幕| 99热国产这里只有精品6| 精品亚洲成a人片在线观看| 色综合欧美亚洲国产小说| 男女无遮挡免费网站观看| av天堂久久9| 亚洲久久久国产精品| 日日爽夜夜爽网站| 老司机影院毛片| 我要看黄色一级片免费的| 亚洲av日韩在线播放| 97人妻天天添夜夜摸| 国产成人影院久久av| 丝袜在线中文字幕| 国产成人av教育| 91成人精品电影| 高清av免费在线| 另类亚洲欧美激情| 人人妻人人添人人爽欧美一区卜| videosex国产| 女性被躁到高潮视频| 只有这里有精品99| 亚洲专区中文字幕在线| 99久久精品国产亚洲精品| 在线观看一区二区三区激情| 美女扒开内裤让男人捅视频| 欧美激情 高清一区二区三区| av在线播放精品|