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

    無機肥施用對農田干濕循環(huán)土壤收縮特性的影響

    2021-12-16 11:03:40劉暢齊偉王策張展羽
    灌溉排水學報 2021年11期
    關鍵詞:級數(shù)無機含水率

    劉暢,齊偉,王策,張展羽*

    無機肥施用對農田干濕循環(huán)土壤收縮特性的影響

    劉暢1,齊偉2,王策1,張展羽1*

    (1.河海大學 農業(yè)科學與工程學院,南京 211100;2.河海大學 水利水電學院,南京 210098)

    【】研究無機肥施用對農田干濕循環(huán)土壤收縮特性的影響。選取南方地區(qū)粉壤土為試驗土樣,設置尿素、硫酸鉀2種肥料配置浸潤液(質量濃度設置為1、2、3 g/L),開展多級循環(huán)脫濕試驗。①土壤收縮應變隨著循環(huán)級數(shù)的推進而減小,且逐漸趨于穩(wěn)定。無機肥的添加減少了土壤收縮應變,且隨著循環(huán)級數(shù)的推進愈加明顯。②VG-Peng模型擬合土壤收縮特性曲線效果較好,各處理均具有結構收縮、比例收縮、殘余收縮、零收縮4個階段,干濕循環(huán)使土壤結構收縮段增加,比例收縮段和殘余收縮段減小。無機肥的添加使土壤比例收縮段孔隙比范圍減少。③土壤收縮經歷了從完全的縱向收縮-以縱向收縮為主-以橫向收縮為主3個階段,多次循環(huán)后土壤收縮各向異性減弱。無機肥添加和干濕循環(huán)條件不同程度減弱了土壤的收縮特性,多次循環(huán)后無機肥累積是土壤收縮能力下降的主要原因。

    干濕循環(huán);無機肥;土壤;收縮特性曲線;幾何因子

    0 引言

    【研究意義】自然界中,土壤體積收縮主要因為土壤含水率降低,土壤基質吸力增大,土壤顆粒間產生吸力導致土壤顆粒和團聚體相互吸引并發(fā)生重排,最終使土壤孔隙體積減小[1]。由于土壤失水收縮總伴隨著土壤基質吸力的增加,當土壤薄弱區(qū)土體拉應力超過土壤顆粒間有效應力時,土壤開裂并逐漸發(fā)育成裂隙[2],嚴重影響土壤結構的穩(wěn)定性。在工程地質方面,土壤裂隙的產生降低了土體抗剪強度,易導致地基沉陷、邊坡失穩(wěn),引發(fā)一系列工程災害[3]。在農業(yè)水土工程方面,土壤干縮裂隙引發(fā)的優(yōu)先流效應增大了指流、深層滲漏的風險,導致灌溉水分、肥料利用率下降,阻礙農業(yè)生產發(fā)展[4-5]。因此,土壤收縮特性作為裂縫發(fā)育的關鍵因素,受到地質科學、土木工程、農業(yè)水土工程等多學科學者的重點關注。

    【研究進展】土壤收縮特性常采用土壤收縮特性曲線描述,國內外學者提出多種模型用以擬合土壤收縮特性曲線,包括三直線模型、多項式模型、Sigmoid模型、Logistic模型等[6]。邢旭光等[7]通過Logistic模型擬合土壤收縮曲線,發(fā)現(xiàn)土壤的收縮應變有隨體積質量增加而減小的趨勢,且土壤中砂粒量與土壤收縮特征指標呈負相關,粉粒和黏粒則相反。Wei等[8]通過分形維數(shù)法發(fā)現(xiàn)花崗巖土壤的收縮能力與土壤粒徑分布和壓實度有關。Li等[9]采用壓汞法和電子顯微鏡掃描研究了土壤收縮過程中微觀結構的變化,發(fā)現(xiàn)以離心機法測得的土壤收縮特性曲線并未呈現(xiàn)出典型的“S”形,土壤在干燥過程中大孔隙增加而團聚體間小孔隙無明顯變化。

    【切入點】日常農業(yè)生產活動中,為改善土壤肥力、提高作物產量,肥料的應用越來越廣泛。隨著農業(yè)技術的發(fā)展,傳統(tǒng)撒施技術也在向灌溉施肥轉變。灌溉施肥在一定程度上提高了化肥利用率[10],降低了養(yǎng)分流失引起的農田土壤生態(tài)環(huán)境污染風險[11]。同時也對土壤含鹽量[12]、有機質量、水力特性[13]、結構特性產生影響,從而影響土壤收縮特性。反之,土壤收縮特性又對土壤內水肥、溶質遷移有重要影響,因此開展無機肥添加對土壤收縮特性研究尤為重要。此外,考慮到實際農業(yè)生產中,作物生育期內通常會經歷數(shù)次干旱和灌排過程,土壤也因此在干濕循環(huán)中反復干濕縮脹,土壤結構穩(wěn)定性以及孔隙頻率分布也在循環(huán)中不斷改變,這也將導致土壤收縮特性的差異性[14]。

    【擬解決的關鍵問題】為研究干濕循環(huán)中農田土壤反復縮脹導致孔隙結構的變化以及無機肥添加對土壤收縮特性的共同影響,采用室內試驗,結合數(shù)字圖像處理技術,研究多級干濕循環(huán)下無機肥添加對農田土壤收縮特性影響,并引入幾何因子對土壤收縮過程中體變各向異性進行評價,以期為土壤水肥運移機理研究、水肥一體化灌溉制度制定以及防治面源污染提供理論基礎。

    1 試驗設計與方法

    1.1 試驗材料

    供試土壤選自河海大學江寧節(jié)水園區(qū)試驗田旱地表層土壤,采集深度為20 cm。原狀土碾碎、烘干、過2 mm過篩后,用吸管法測定土壤機械組成,用比重瓶法測土壤比重,依據(jù)國際土壤質地分類方法進行土壤分類。土壤物理特性如表1所示。試驗土壤全鉀量、全氮量分別為11.32、1.15 g/kg。

    表1 供試土壤物理特性

    1.2 試驗方案設計

    試驗在室內進行,將尿素、硫酸鉀溶于水配制成肥液并設置3種質量濃度梯度:1、2、3 g/L(記作N-1、N-2、N-3、K-1、K-2、K-3),以純水(CK)為對照組,共7組處理。供試土壤按1.35 g/cm3預設土壤體積質量填入土箱(20 cm×20 cm×20 cm),浸潤液飽和處理24 h。試驗采用環(huán)刀取樣,取200 cm3環(huán)刀(?70 mm×52 mm),環(huán)刀底部平鋪無紡布防止土粒從透水孔隙漏出,內壁涂抹凡士林以減小邊壁摩擦,外壁貼布基膠帶防止水分側滲,環(huán)刀取充分飽和后土樣,每組處理取3次重復。試樣制備完成后測量初始質量、截面積、高度,置于40 ℃烘箱中干燥至恒質量,后置于105 ℃烘箱完全烘干。干燥過程中每隔2 h測量試驗土樣質量和體積變化,電子天平測質量,數(shù)碼相機拍攝試樣截面(圖1(a)),游標卡尺測定試驗土樣高度(圖1(b))。干燥結束后測各組數(shù)據(jù)再將環(huán)刀置于各處理溶液中,控制溶液液面高度略低于環(huán)刀頂部,試樣通過底部透水孔吸濕飽和,飽和處理持續(xù)24 h。試樣從飽和干燥至恒定質量,又增濕至飽和的過程稱為1次干濕循環(huán),試驗共進行4級循環(huán)。試驗所用電子天平精度為0.1 g,電子游標卡尺精度為0.01 mm,相機為佳能EOS 850D,像素2 400萬。拍攝過程中相機和試樣位置均固定不變并阻隔室外光源,采用LED燈均勻照射試樣表面。

    1.3 分析方法

    1.3.1 數(shù)字圖像處理

    試樣橫截面積通過數(shù)字圖像處理技術提取計算,圖像處理過程如圖1(c)所示。首先對初始圖像進行批量裁剪預處理,取70 mm×70 mm作為研究區(qū)域,像素設置為1 600×1 600,然后將裁剪后圖像轉為灰度圖。再運用OTSU算法[14]獲取該灰度圖的閾值,由于干燥過程中試樣失水收縮截面積減小,試樣邊緣與環(huán)刀分離,其收縮區(qū)域與土壤基質區(qū)域像素灰度值對比明顯,二者類間方差最大時對應最佳閾值。將灰度值大于的像素設為1(白),代表土壤收縮區(qū)域,灰度值小于的像素設為0(黑),代表環(huán)刀外區(qū)域(計算中去除)和土壤基質區(qū)域,從而得到二值圖像,由此計算試樣橫截面積。以上處理過程均在Matlab R2018中進行,通過自編程序實現(xiàn)。

    圖1 試驗方法及圖像處理

    1.3.2 試樣收縮應變

    試樣干燥過程中隨著水分蒸發(fā)體積收縮,橫向表現(xiàn)為截面積減小、縱向表現(xiàn)為高度降低。采用橫向收縮應變r、縱向收縮應變s、體積收縮應變v描述土壤收縮程度,表達式為:

    r=(0-)/0×100%, (1)

    s=(0-)/0×100%, (2)

    v=(0-)/0×100%, (3)

    式中:、、分別為各循環(huán)脫水過程中試樣直徑(mm)、高度(mm)、體積(mm3);0、0、0分別為各循環(huán)前試樣初始直徑(mm)、高度(mm)、體積(mm3)。試樣收縮應變取3次重復平均值采用SPSS軟件進行統(tǒng)計學分析。

    1.3.3 土壤收縮特性曲線模型

    土壤收縮特性曲線采用2005年Peng等[16]提出的VG-Peng模型擬合,表達式為:

    1.3.4 收縮幾何因子

    土壤在收縮過程中具有各向異性,表現(xiàn)為橫向收縮與縱向沉降的差異性。幾何因子s將土壤體積收縮縱向收縮相結合,用以描述土壤體積變化過程中的各向異性[18]。圓柱形式樣s計算式為:

    式中:參數(shù)意義同上,s值在不同區(qū)域有不同物理意義:

    ①s=1,土體僅發(fā)生縱向沉降;

    ② 1≤s≤3,土體收縮以縱向沉降為主導;

    ③s=3,土體收縮呈各項同性;

    ④s>3,土體收縮軸向收縮為主導;

    2 結果與分析

    2.1 土壤收縮應變

    干燥過程中,隨著土壤孔隙水的蒸發(fā),土壤含水率降低,體積質量增加。圖2、圖3分別反映了試樣橫向收縮應變、縱向收縮應變隨質量含水率變化。表2反映了土壤各向極限收縮應變。首次循環(huán)中,土體橫向收縮隨含水率的降低而增大,變化速率呈慢-快-慢趨勢,表現(xiàn)為土體干燥過程橫向收縮應變先隨著含水率降低緩慢增大,在含水率降至0.31 g/g左右時變化速率加快,最后在含水率降至0.07 g/g左右時保持不變,完全干燥時各處理橫向收縮應變在9.93%~10.91%之間。首次循環(huán)中,縱向收縮也隨含水率降低而變大,與橫向收縮不同,干燥開始時縱向收縮率即隨含水率的降低而快速增加,在含水率降至0.10 g/g左右時基本不變,干燥結束時各組處理極限縱向收縮應變在7.97%~8.94%之間。隨著循環(huán)級數(shù)的推進,土壤橫、縱向收縮應變變化趨勢變緩,完全干燥時土樣橫向收縮應變、縱向收縮應變、體積收縮應變分別從初次循環(huán)的9.93%~10.91%、7.97%~8.94%、25.41%~27.68%降至末次循環(huán)的4.85%~6.53%、4.12%~5.16%、14.56%~17.99%,表明土壤收縮能力隨循環(huán)級數(shù)推進而減弱。

    圖2 土壤失水過程中橫向收縮應變

    圖3 土壤失水過程中縱向應變

    對比各處理在3、4次循環(huán)中土壤體積收縮應變發(fā)現(xiàn):CK在3、4次循環(huán)中體積收縮應變差值為0.17%,而無機肥處理組同比差值普遍在1%以上,這說明在3次循環(huán)后無機肥仍對土壤收縮起抑制作用。為研究肥液質量濃度對土壤各向收縮程度影響,采用SPSS軟件對肥液質量濃度和土體各向收縮應變進行相關性分析(表3)。從整體上來看,2種肥液質量濃度與土壤各向收縮應變均具有較高負相關性,肥液質量濃度與橫向收縮應變和體積收縮應變多呈現(xiàn)出極顯著負相關(<0.01),與縱向收縮應變多表現(xiàn)出顯著負相關(<0.05)。隨著循環(huán)級數(shù)的推進,肥液質量濃度與土體各向收縮應變顯著性增強,在最后1次循環(huán)中,2種肥料處理質量濃度與各向收縮應變均極顯著負相關(<0.01)。

    表2 各循環(huán)級數(shù)下不同肥料處理試樣極限收縮應變

    表3 肥液質量濃度與土壤各向極限收縮應變相關性分析結果

    注 *和**分別表示在<0.05和<0.01水平(雙側)顯著和極顯著相關。

    2.2 土壤收縮階段特性

    采用VG-Peng模型擬合土壤收縮特性曲線(圖4),擬合效果較好,相關系數(shù)2均大于0.99。由圖4可知,土壤收縮特性曲線整體呈“S”形,各組處理均有結構收縮、比例收縮、殘余收縮、0收縮4個階段,各收縮階段起止點見表4。首次循環(huán)中各組處理土樣s在1.12~1.22之間,r在0.54~0.66之間。各組處理比例收縮區(qū)水分、孔隙比分別占總損失的41.4%~51.5%、70.3%~76.9%,表明土壤收縮以比例收縮段為主。隨著循環(huán)級數(shù)的推進,各組處理s減小,r變化不大。

    s隨循環(huán)級數(shù)減小是由于土壤的不可逆收縮現(xiàn)象,即經歷了從飽和至完全干燥的土壤,即使再次吸水飽和,其體積膨脹也不能達到初次浸潤體積。由于s與土壤收縮應變均隨循環(huán)級數(shù)的推進而減小,導致r變化并不明顯。分析各處理土壤收縮區(qū)域變化發(fā)現(xiàn),隨著循環(huán)級數(shù)的推進,土壤結構收縮區(qū)范圍變大,結構收縮區(qū)孔隙比變化幅度分別從初次循環(huán)0.020~0.060增加到末次循環(huán)0.180~0.232,表明土壤團聚體間大孔隙隨干濕循環(huán)增多。比例收縮區(qū)在各級循環(huán)中均為收縮主體區(qū)域,但水分、孔隙比占比隨循環(huán)級數(shù)的推進不斷減少,表明土壤收縮程度逐漸減弱。殘余收縮區(qū)水分、孔隙比占比隨循環(huán)級數(shù)的推進也有所減少,表明土壤團聚體內孔隙減少。

    圖4 土壤收縮特性曲線擬合

    對比同級循環(huán)各組處理,發(fā)現(xiàn)無機肥處理土壤收縮特性曲線始終位于CK組上方,而曲線形態(tài)并未發(fā)生明顯改變。無機肥的添加使土壤飽和孔隙比s和殘余孔隙比r有不同程度的增加,在相同含水率下無機肥處理往往有更大的孔隙比。對比各處理收縮區(qū)域發(fā)現(xiàn),各處理土壤結構收縮區(qū)和殘余收縮區(qū)并無顯著差異。而比例收縮區(qū)在無機肥處理下孔隙比范圍減小,且在高質量濃度肥液處理下和多次循環(huán)中表現(xiàn)更加明顯。比例收縮區(qū)作為收縮的主體區(qū)域,其孔隙比范圍降低是無機肥導致土壤收縮應變減小的主要原因。

    表4 土壤各收縮階段起止點

    圖5 土壤失水過程中幾何因子變化

    2.3 土壤收縮各向異性

    收縮幾何因子反映了土壤橫、縱向收縮應變的差異性,應用圓柱形s公式計算得到幾何因子隨含水率變化如圖5所示。首次循環(huán)中,各處理s值以1為起點隨含水率的減小而增加,在含水率降至0.24~0.27 g/g左右時s增至3,繼續(xù)增加最終穩(wěn)定在3.22~3.66。結合實驗現(xiàn)象來看,干燥開始時s=1,土體表現(xiàn)為完全的縱向沉降。隨著土壤水分流失,基質吸力增加,開始出現(xiàn)橫向收縮,此時s介于1~3之間,說明此時土體仍以縱向收縮為主。隨著含水率的降低,土壤基質吸力進一步增加,由基質吸力引起的張拉應力的增加使土壤體積變化以縱向收縮為主導轉變?yōu)橐詸M向收縮為主導(s>3)。在含水率降至0.1 g/g左右時土壤體積幾乎不變,s在此時增至最大后保持不變。對比各組處理,發(fā)現(xiàn)CK與無機肥添加處理在開始失水時均表現(xiàn)出了完全的縱向收縮,且s均隨著含水率的減小而增加最后趨于穩(wěn)定,各處理s變化斜率穩(wěn)定時s值并無明顯差異。對比各級循環(huán)s變化發(fā)現(xiàn),隨著循環(huán)級數(shù)的推進,土壤體變以縱向收縮為主導的區(qū)域增多,s的變化更加緩和,完全干燥時s減小。

    3 討論

    3.1 干濕循環(huán)對土壤收縮特性影響

    土壤結構影響水、空氣、溶質的存儲和運移,是植物根系生長和微生物活動的媒介,是土壤最基本的物理性質[20]。試驗在反復干濕過程中觀察到土壤的不可逆收縮現(xiàn)象,說明土壤干濕循環(huán)過程中黏土顆粒和團聚體反復縮脹導致土壤結構發(fā)生改變,土壤收縮特性也因此改變。表現(xiàn)為隨著循環(huán)級數(shù)的推進,土壤橫、縱向收縮程度減弱,s減小,r增大,孔隙比變化幅度減小,土壤收縮能力降低。通過分析土壤收縮特性曲線發(fā)現(xiàn):干濕循環(huán)使土壤結構收縮段增加,比例收縮段、殘余收縮段減少。土壤各收縮區(qū)域代表土壤不同尺度孔隙在失水過程中作用范圍,從側面反映了土壤孔隙結構的變化規(guī)律,根據(jù)土壤各收縮區(qū)域變化趨勢和土壤結構在干濕循環(huán)中演變規(guī)律來看:干濕循環(huán)過程中土壤團聚體破碎重排形成了復雜的板狀、塊狀、團粒結構[21],這些結構較均質結構相對穩(wěn)定,但也導致團聚體間大孔隙增多。然而土壤結構收縮區(qū)的增大并不代表土壤收縮能力加強,因為團聚體間大孔隙在水分損失過程中并不會土壤體積大幅變化。作為收縮的主體部分,比例收縮區(qū)占據(jù)了土壤失水過程中大部分體積變化和水分損失,其范圍在干濕循環(huán)過程中不斷減小,表明土壤收縮能力不斷減弱。殘余收縮區(qū)的減少表明團聚體內孔隙占比減少,這一方面可能是由于干濕循環(huán)使團聚體破碎[22],導致團聚體數(shù)量減少,另一方面可能是由于原團聚體破碎重排后形成的新的團聚體孔徑較小。此外,試驗中還發(fā)現(xiàn):土壤橫、縱向收縮應變、體積收縮應變前3次循環(huán)間變化幅度較大,而在3、4次循環(huán)間變化不大,表明干濕循環(huán)對土壤結構和收縮能力影響隨循環(huán)級數(shù)的推進逐漸減弱。這是因為每次循環(huán)中土壤團聚體結構破碎重排都會形成更穩(wěn)定的土結構,干濕循環(huán)引起的體變不可逆程度隨級數(shù)推進逐漸降低。在經歷多次循環(huán)后,土壤收縮量趨于穩(wěn)定值因為此時土壤內團聚體多發(fā)生脹縮而少有團聚體的再造與滑移[23]。

    3.2 無機肥添加對土壤收縮特性影響

    無機肥添加在一定程度上改變了土壤理化性質[24]、結構特性,從而使土壤收縮特性發(fā)生改變,導致土壤增濕脫濕過程中體積變化行為的差異性。研究中發(fā)現(xiàn)無機肥的添加在各級循環(huán)中均不同程度地減小了土壤橫、縱向收縮程度,且隨著肥液質量濃度的增加作用效果更加明顯。通過觀察土壤收縮特性曲線發(fā)現(xiàn),無機肥的添加并未明顯改變結構收縮區(qū)和殘余收縮區(qū)的范圍,說明無機肥的添加并未導致土壤團聚體間大孔隙和團聚體內孔隙結構發(fā)生改變。無機肥添加導致土壤收縮應變的減少主要是由于比例收縮區(qū)孔隙比范圍的減小,即土壤黏粒收縮能力的減弱。這一方面可能是因為隨著土壤孔隙內流體質量濃度的增加,土壤雙電子層中擴散層厚度減小,土顆粒間距減小,土壤收縮能力減弱[25]。另一方面可能是因為隨著孔隙內流體質量濃度的增加,氣-液界面張力(表面張力)減小,導致溶液與土顆粒接觸角減小,曲率半徑增大,導致基質吸力減小,土壤收縮程度降低[26]。此外試驗中還發(fā)現(xiàn),干濕循環(huán)中土壤團聚體結構變化和肥料的累積都有降低土壤收縮應變的趨勢,二者共同作用使土壤的收縮能力明顯減弱。隨著循環(huán)級數(shù)的推進,由干濕循環(huán)引起的土壤收縮能力的變化逐漸減弱,多級循環(huán)后土壤收縮能力下降主要是由無機肥累積引起的,這意味著無機肥的添加在一定程度上加劇了土壤板結,裂隙固化,因此在水肥一體化灌溉過程中,應注意水肥統(tǒng)籌,避免過量施肥。不過,本研究為精確控制試驗變量,未能充分考慮作物根系對無機肥的吸收利用,這與農業(yè)生產大田實際情況有所不同,今后的研究需進一步分析在作物耕種條件下農田土壤收縮特性、干縮裂隙、水肥遷移規(guī)律三者的聯(lián)系。

    4 結論

    1)土壤各向收縮應變隨著循環(huán)級數(shù)的推進而減小,并逐漸趨近穩(wěn)定。無機肥的添加減少了橫、縱向收縮應變,且隨著循環(huán)級數(shù)的推進表現(xiàn)愈加明顯。

    2)VG-Peng模型擬合土壤收縮特性曲線效果較好。干濕循環(huán)中各處理均存在結構收縮、比例收縮、殘余收縮、零收縮4個階段。無機肥的添加使土壤殘余孔隙比、飽和孔隙比增加,同時也導致土壤收縮能力減弱,原因在于比例收縮段孔隙比變化幅度的減少。

    3)土壤收縮經歷了從完全的縱向收縮-以縱向收縮為主-以橫向收縮為主3個階段,多次干濕循環(huán)后土壤體變各向異性減弱。

    [1] TANG C S, SHI B, LIU C, et al. Experimental characterization of shrinkage and desiccation cracking in thin clay layer[J]. Applied Clay Science, 2011, 52(1/2): 69-77.

    [2] WANG C, ZHANG Z Y, LIU Y, et al. Geometric and fractal analysis of dynamic cracking patterns subjected to wetting-drying cycles[J]. Soil and Tillage Research, 2017, 170: 1-13.

    [3] KRISNANTO S, RAHARDJO H, FREDLUND D G, et al. Mapping of cracked soils and lateral water flow characteristics through a network of cracks[J]. Engineering Geology, 2014, 172: 12-25.

    [4] ALLAIRE S E, ROULIER S, CESSNA A J. Quantifying preferential flow in soils: A review of different techniques[J]. Journal of Hydrology, 2009, 378(1/2): 179-204.

    [5] ZHANG Z B, PENG X, ZHOU H, et al. Characterizing preferential flow in cracked paddy soils using computed tomography and breakthrough curve[J]. Soil and Tillage Research, 2015, 146: 53-65.

    [6] CHEN B, MEINERTZHAGEN I A, SHAW S R. Circadian rhythms in light-evoked responses of the fly's compound eye, and the effects of neuromodulators 5-HT and the peptide PDF[J]. Journal of Comparative Physiology A-neuroethology Sensory Neural & Behavioral Physiology, 1999, 185(5): 393-404.

    [7] 邢旭光, 馬孝義. 土壤失水過程對土體收縮特性的影響[J]. 農業(yè)機械學報, 2016, 47(6): 148-153.

    XING Xuguang, MA Xiaoyi. Effects of soil drying process on soil shrinkage characteristic[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(6): 148-153.

    [8] WEI Y J, WU X L, XIA J W, et al. Relationship between granitic soil particle-size distribution and shrinkage properties based on multifractal method[J]. Pedosphere, 2020, 30(6): 853-862.

    [9] LI L C, LI X, WANG L, et al. The effects of soil shrinkage during centrifuge tests on SWCC and soil microstructure measurements[J]. Bulletin of Engineering Geology and the Environment, 2020, 79(7): 3 879-3 895.

    [10] 劉德利, 吳文勇, 肖娟, 等. 不同灌溉施肥時機對稻田肥料分布和水稻生長的影響[J]. 灌溉排水學報, 2021, 40(4): 29-36.

    LIU Deli, WU Wenyong, XIAO Juan, et al. The timing of irrigation and fertilization affects fertilizer distribution and rice growth[J]. Journal of Irrigation and Drainage, 2021, 40(4): 29-36.

    [11] BRUNETTI G, ?IM?NEK J, BAUTISTA E. A hybrid finite volume-finite element model for the numerical analysis of furrow irrigation and fertigation[J]. Computers and Electronics in Agriculture, 2018, 150: 312-327.

    [12] 張文平, 聶衛(wèi)波, 馬敬東, 等. 肥液濃度對土壤水力特性參數(shù)的影響[J]. 水土保持學報, 2020, 34(4): 63-71.

    ZHANG Wenping, NIE Weibo, MA Jingdong, et al. Effect of fertilizer concentration on soil hydraulic parameters[J]. Journal of Soil and Water Conservation, 2020, 34(4): 63-71.

    [13] 孫樹臣, 趙鑫, 翟勝, 等. 生物質炭施用對土壤水分影響的研究進展[J]. 灌溉排水學報, 2020, 39(11): 112-119.

    SUN Shuchen, ZHAO Xin, ZHAI Sheng, et al. The effect of biochar amendment on soil moisture: A review[J]. Journal of Irrigation and Drainage, 2020, 39(11): 112-119.

    [14] 鞏學鵬, 唐朝生, 施斌, 等. 黏性土干/濕過程中土結構演化特征研究進展[J]. 工程地質學報, 2019, 27(4): 775-793.

    GONG Xuepeng, TANG Chaosheng, SHI Bin, et al. Evolution of soil microstructure during drying and wetting[J]. Journal of Engineering Geology, 2019, 27(4): 775-793.

    [15] 劉健莊, 栗文青. 灰度圖象的二維Otsu自動閾值分割法[J]. 自動化學報, 1993, 19(1): 101-105.

    LIU Jianzhuang, LI Wenqing. The automatic thresholding of gray-level pictures via two-dimensional otsu method[J]. Acta Automatica Sinica, 1993, 19(1): 101-105.

    [16] PENG X, HORN R. Modeling soil shrinkage curve across a wide range of soil types[J]. Soil Science Society of America Journal, 2005, 69(3): 584-592.

    [17] PENG X, HORN R, DEERY D, et al. Influence of soil structure on the shrinkage behaviour of a soil irrigated with saline - sodic water[J]. Soil Research, 2005, 43(4): 555.

    [18] MISHRA P N, ZHANG Y, BHUYAN M H, et al. Anisotropy in volume change behaviour of soils during shrinkage[J]. Acta Geotechnica, 2020, 15(12): 3399-3414.

    [19] BRAUDEAU E, COSTANTINI J M, BELLIER G, et al. New device and method for soil shrinkage curve measurement and characterization[J]. Soil Science Society of America Journal, 1999, 63(3): 525-535.

    [20] DIEL J, VOGEL H J, SCHLüTER S. Impact of wetting and drying cycles on soil structure dynamics[J]. Geoderma, 2019, 345: 63-71.

    [21] HUSSEIN J, ADEY M A. Changes of structure and tilth mellowing in a Vertisol due to wet/dry cycles in the liquid and vapour phases[J]. European Journal of Soil Science, 1995, 46(3): 357-368.

    [22] MA R M, CAI C F, LI Z X, et al. Evaluation of soil aggregate microstructure and stability under wetting and drying cycles in two Ultisols using synchrotron-based X-ray micro-computed tomography[J]. Soil and Tillage Research, 2015, 149: 1-11.

    [23] AIRò FARULLA C, FERRARI A, ROMERO E. Volume change behaviour of a compacted scaly clay during cyclic suction changes[J]. Canadian Geotechnical Journal, 2010, 47(6): 688-703.

    [24] 丁嘉寧, 梁利寶, 馮鵬艷, 等. 施肥對采煤塌陷復墾土壤理化性狀和細菌群落結構影響的研究[J]. 灌溉排水學報, 2020, 39(12): 64-70.

    DING Jianing, LIANG Libao, FENG Pengyan, et al. Effects of fertilization on physical and chemical properties and bacterial community of reclaimed soil from coal mining subsidence[J]. Journal of Irrigation and Drainage, 2020, 39(12): 64-70.

    [25] 孫凱強. 鹽分對土體干縮開裂的影響研究[D]. 南京: 南京大學, 2017.

    SUN Kaiqiang. The research on the influence of salt on soil desiccation cracking[D]. Nanjing: Nanjing University, 2017.

    [26] 羅茂泉, 楊松, 馬澤慧, 等. 干濕循環(huán)下氣-液界面張力對黏土收縮開裂的影響[J]. 土壤學報, 2018, 55(2): 369-379.

    LUO Maoquan, YANG Song, MA Zehui, et al. Effect of gas-liquid interfacial tension on shrinkage cracking of clay as affected by wetting-drying alternation[J]. Acta Pedologica Sinica, 2018, 55(2): 369-379.

    The Effects of Inorganic Fertilization on Soil Shrinkage in Wetting-drying Cycles

    LIU Chang1,QI Wei2, WANG Ce1, ZHANG Zhanyu1*

    (1. College of Agricultural Science and Engineering, Hohai University, Nanjing 211100, China;2. College of Water Conservancy and Hydropower Engineering, Hohai University, Nanjing 210098, China)

    【】Swelling and shrinkage are important process reshaping soil structure. It is impacted by many factors, and the objective of this paper is to investigate how fertilization impacts this process.【】We used a silty loam in the experiment. The soil was fertilized using different inorganic N and K fertilizers at concentrations ranging from 1 g/L to 3 g/L. The emergence and development of the shrinkage during and following multiple wetting-drying cycles were measured in each treatment.【】① The shrinkage in the transverse and longitudinal directions both decreased asymptotically with the increase in the number of drying-wetting cycles. Adding inorganic fertilizer reduced the shrinkage in both transverse and longitudinal directions, and the reduction became more significant as the wetting-dying cycles increased. ② The VG-Peng model can fit the effects of the fertilization on soil shrinkage, and the shrinkages in each treatment can be classified into four types: structural shrinkage, proportional shrinkage, residual shrinkage and zero shrinkage. Increasing drying-wetting cycles increased the structural shrinkage at the expense of proportional shrinkage and residual shrinkage. The addition of the inorganic fertilizers reduced the void ratio of the soil in the shrunk region. ③Soil shrinkage experienced three stages: complete longitudinal shrinkage, longitudinal-dominant shrinkage, and transverse-dominant shrinkage, with the anisotropy of the shrinkages decreasing after several drying-wetting cycles.【】Applying inorganic fertilizers reduced soil shrinkage due to wetting-drying cycles, and the inorganic fertilizers accumulated after multiple wetting-drying cycles due to the decline in soil shrinkage.

    wetting-drying cycle; inorganic fertilization; soil; contraction curve; geometric factors

    S152.4; S143

    A

    10.13522/j.cnki.ggps.2021263

    1672 - 3317(2021)11 - 0044 - 07

    2021-06-21

    國家自然科學基金面上項目(51879071,51579069)

    劉暢(1998-),男。碩士研究生,主要從事農田土壤干縮裂隙方向研究。E-mail:Liuchang_@hhu.edu.cn

    張展羽(1957-),男。教授,主要從事節(jié)水灌溉理論與技術研究。E-mail: zhanyu @hhu.edu.cn

    劉暢, 齊偉, 王策, 等. 無機肥施用對農田干濕循環(huán)土壤收縮特性的影響[J]. 灌溉排水學報, 2021, 40(11): 44-50.

    LIU Chang, QI Wei, WANG Ce, et al. The Effects of Inorganic Fertilization on Soil Shrinkage in Wetting-drying Cycles[J]. Journal of Irrigation and Drainage, 2021, 40(11): 44-50.

    責任編輯:白芳芳

    猜你喜歡
    級數(shù)無機含水率
    昆明森林可燃物燃燒機理研究
    無機滲透和促凝劑在石材防水中的應用
    石材(2020年9期)2021-01-07 09:30:04
    基于表面能的濕煤顆粒含水率表征方法
    加快無機原料藥產品開發(fā)的必要性和途徑
    Dirichlet級數(shù)及其Dirichlet-Hadamard乘積的增長性
    弱膨脹土增濕變形量試驗及路堤填筑分析
    幾個常數(shù)項級數(shù)的和
    原油含水率在線測量技術研究
    電子制作(2016年1期)2016-11-07 08:42:56
    有機心不如無機心
    山東青年(2016年2期)2016-02-28 14:25:31
    p級數(shù)求和的兩種方法
    国产伦一二天堂av在线观看| 色哟哟哟哟哟哟| 日本黄色片子视频| 噜噜噜噜噜久久久久久91| 嫩草影视91久久| 9191精品国产免费久久| 国内精品久久久久久久电影| 国产黄色小视频在线观看| 色噜噜av男人的天堂激情| 国产探花极品一区二区| 老司机午夜福利在线观看视频| 草草在线视频免费看| 免费av不卡在线播放| 色尼玛亚洲综合影院| 日本免费a在线| 亚洲欧美日韩无卡精品| 成人国产综合亚洲| av黄色大香蕉| 日韩高清综合在线| 91麻豆精品激情在线观看国产| 亚洲色图av天堂| 天堂影院成人在线观看| 男人和女人高潮做爰伦理| 国产成人福利小说| 此物有八面人人有两片| 成年版毛片免费区| 亚洲avbb在线观看| 亚洲欧美日韩卡通动漫| 国产主播在线观看一区二区| 国产 一区 欧美 日韩| 亚洲av.av天堂| 国产一区二区三区视频了| 又爽又黄a免费视频| 免费观看的影片在线观看| 精品免费久久久久久久清纯| 婷婷色综合大香蕉| 亚洲经典国产精华液单 | 日本一本二区三区精品| 国产真实乱freesex| 免费观看人在逋| 日韩欧美一区二区三区在线观看| 永久网站在线| 欧美绝顶高潮抽搐喷水| 老女人水多毛片| 中文字幕精品亚洲无线码一区| 亚洲最大成人av| 天堂√8在线中文| 亚洲av.av天堂| 中文亚洲av片在线观看爽| 国产精品一区二区免费欧美| 国产成人aa在线观看| 午夜亚洲福利在线播放| 丰满人妻一区二区三区视频av| 中文字幕久久专区| 麻豆成人午夜福利视频| 一个人免费在线观看电影| 色综合亚洲欧美另类图片| 91字幕亚洲| 国产黄片美女视频| 亚洲人成网站高清观看| 制服丝袜大香蕉在线| 精品人妻1区二区| 看十八女毛片水多多多| 日韩欧美国产在线观看| 国产精品一区二区免费欧美| 欧美成人免费av一区二区三区| 国产av一区在线观看免费| 亚洲人与动物交配视频| 亚洲无线在线观看| 久久精品国产自在天天线| 国产精品久久久久久久电影| 国产淫片久久久久久久久 | 嫩草影院入口| 久久久精品欧美日韩精品| 最新中文字幕久久久久| 99国产极品粉嫩在线观看| 亚洲精品在线美女| 日本精品一区二区三区蜜桃| 狠狠狠狠99中文字幕| 国产视频一区二区在线看| 久久久久免费精品人妻一区二区| 日韩欧美国产在线观看| 成年免费大片在线观看| 亚洲成人久久爱视频| 国产精品女同一区二区软件 | 欧美bdsm另类| 91久久精品电影网| 成人av一区二区三区在线看| 中文在线观看免费www的网站| 别揉我奶头~嗯~啊~动态视频| 免费看美女性在线毛片视频| 免费电影在线观看免费观看| av在线观看视频网站免费| 国产精品98久久久久久宅男小说| 欧美最黄视频在线播放免费| 亚洲人成伊人成综合网2020| 麻豆成人午夜福利视频| 一夜夜www| 欧美性猛交黑人性爽| www.999成人在线观看| 婷婷色综合大香蕉| 好男人电影高清在线观看| 小说图片视频综合网站| 国产精品99久久久久久久久| 美女被艹到高潮喷水动态| av黄色大香蕉| 国产精品永久免费网站| 欧美国产日韩亚洲一区| 亚洲av不卡在线观看| 久久精品夜夜夜夜夜久久蜜豆| 亚洲成人久久性| .国产精品久久| 小蜜桃在线观看免费完整版高清| 男女那种视频在线观看| 欧美日韩亚洲国产一区二区在线观看| 国产精品亚洲av一区麻豆| 午夜免费激情av| 舔av片在线| 久久亚洲精品不卡| 国产三级中文精品| 久久国产精品影院| 国产成+人综合+亚洲专区| 人妻久久中文字幕网| 一进一出抽搐动态| 午夜精品久久久久久毛片777| 国产成人av教育| 精品久久久久久久久久久久久| 美女被艹到高潮喷水动态| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产亚洲精品久久久com| 哪里可以看免费的av片| 国产精品久久视频播放| 天堂动漫精品| www.色视频.com| 亚洲成人精品中文字幕电影| 久久精品人妻少妇| netflix在线观看网站| 丰满人妻熟妇乱又伦精品不卡| 97超级碰碰碰精品色视频在线观看| 1024手机看黄色片| 色播亚洲综合网| 中文字幕熟女人妻在线| 波多野结衣高清作品| 国产欧美日韩一区二区三| 热99在线观看视频| www日本黄色视频网| 亚洲美女黄片视频| 午夜免费成人在线视频| 亚洲成人免费电影在线观看| 欧美日韩瑟瑟在线播放| 91九色精品人成在线观看| 亚洲av成人精品一区久久| 国产精品伦人一区二区| 午夜激情欧美在线| 精品久久久久久久末码| 99国产综合亚洲精品| АⅤ资源中文在线天堂| 一a级毛片在线观看| 日韩欧美国产一区二区入口| 午夜福利在线观看免费完整高清在 | 国产一级毛片七仙女欲春2| 国产精品不卡视频一区二区 | 精品国产亚洲在线| 免费在线观看亚洲国产| 久久精品综合一区二区三区| 久久久精品欧美日韩精品| 国产精品一区二区性色av| 丝袜美腿在线中文| 最新中文字幕久久久久| 国产高清激情床上av| 乱人视频在线观看| 好男人在线观看高清免费视频| 男人和女人高潮做爰伦理| 中国美女看黄片| 18美女黄网站色大片免费观看| 欧美zozozo另类| 日日干狠狠操夜夜爽| 欧美日韩乱码在线| 免费观看人在逋| 国产高清有码在线观看视频| 欧美zozozo另类| 男女视频在线观看网站免费| 一边摸一边抽搐一进一小说| 久久亚洲真实| 日本免费a在线| 欧美日韩乱码在线| 国产成人aa在线观看| 欧美极品一区二区三区四区| 怎么达到女性高潮| 国产成+人综合+亚洲专区| 国产又黄又爽又无遮挡在线| 欧美另类亚洲清纯唯美| 欧美成人免费av一区二区三区| 国产精品久久久久久久久免 | 国产麻豆成人av免费视频| 日韩有码中文字幕| 国产麻豆成人av免费视频| 欧美色欧美亚洲另类二区| 精品熟女少妇八av免费久了| 欧美成狂野欧美在线观看| 久久久久国内视频| 成年人黄色毛片网站| 1024手机看黄色片| 97超级碰碰碰精品色视频在线观看| 免费在线观看成人毛片| 91在线精品国自产拍蜜月| 在线观看免费视频日本深夜| 怎么达到女性高潮| 久久久久性生活片| 亚洲人成电影免费在线| 狂野欧美白嫩少妇大欣赏| 一进一出抽搐动态| 国产中年淑女户外野战色| 欧美黄色淫秽网站| 亚洲久久久久久中文字幕| 中文字幕人妻熟人妻熟丝袜美| 看片在线看免费视频| 精品不卡国产一区二区三区| 波多野结衣高清无吗| 美女大奶头视频| 亚洲国产精品sss在线观看| 久久婷婷人人爽人人干人人爱| 亚洲,欧美精品.| 一二三四社区在线视频社区8| 亚洲av日韩精品久久久久久密| a级毛片a级免费在线| 中文字幕精品亚洲无线码一区| 久久久国产成人免费| 国产欧美日韩一区二区精品| 美女 人体艺术 gogo| 桃色一区二区三区在线观看| 久久草成人影院| 午夜久久久久精精品| 国产在线精品亚洲第一网站| 真人做人爱边吃奶动态| 久久久久久久久中文| 国产精品乱码一区二三区的特点| 亚洲精品一卡2卡三卡4卡5卡| 国产一区二区三区视频了| 亚洲 国产 在线| 村上凉子中文字幕在线| 国产精品av视频在线免费观看| 久久精品国产清高在天天线| 亚洲18禁久久av| 久久精品综合一区二区三区| av在线天堂中文字幕| 黄色视频,在线免费观看| 欧美日韩中文字幕国产精品一区二区三区| 俄罗斯特黄特色一大片| 免费人成视频x8x8入口观看| 蜜桃亚洲精品一区二区三区| 久久久久国内视频| 色播亚洲综合网| 亚洲欧美日韩高清专用| 精品无人区乱码1区二区| 91麻豆av在线| 国产探花在线观看一区二区| 欧美+日韩+精品| 午夜两性在线视频| 日韩大尺度精品在线看网址| 欧美最黄视频在线播放免费| 757午夜福利合集在线观看| 不卡一级毛片| 日韩中字成人| 久久人人爽人人爽人人片va | 69人妻影院| 日本一二三区视频观看| 首页视频小说图片口味搜索| 乱码一卡2卡4卡精品| 9191精品国产免费久久| 午夜两性在线视频| 1000部很黄的大片| 人人妻人人澡欧美一区二区| 国产私拍福利视频在线观看| 美女xxoo啪啪120秒动态图 | 国产成+人综合+亚洲专区| 国产亚洲精品久久久久久毛片| 人人妻,人人澡人人爽秒播| 国产精品久久久久久精品电影| 淫秽高清视频在线观看| 亚洲精华国产精华精| 国内精品久久久久久久电影| 国产成人欧美在线观看| 国产精品电影一区二区三区| 亚洲人成网站在线播放欧美日韩| 波多野结衣高清无吗| 午夜福利18| 神马国产精品三级电影在线观看| 色噜噜av男人的天堂激情| 国产精品,欧美在线| a级毛片免费高清观看在线播放| 成人午夜高清在线视频| 久久久国产成人免费| 亚洲欧美精品综合久久99| 精品一区二区三区av网在线观看| 精品不卡国产一区二区三区| 国产一区二区三区视频了| 身体一侧抽搐| 亚洲国产精品999在线| 国产精品免费一区二区三区在线| 九九热线精品视视频播放| 精品久久久久久久末码| 国产麻豆成人av免费视频| 午夜福利成人在线免费观看| 精品一区二区三区av网在线观看| 亚洲最大成人中文| 亚洲 国产 在线| 婷婷六月久久综合丁香| 亚洲无线在线观看| 岛国在线免费视频观看| 欧美乱妇无乱码| 一级av片app| 窝窝影院91人妻| 白带黄色成豆腐渣| 丁香欧美五月| 成年女人永久免费观看视频| 亚洲av成人不卡在线观看播放网| 日韩国内少妇激情av| 成人高潮视频无遮挡免费网站| 午夜两性在线视频| www日本黄色视频网| 免费高清视频大片| 欧美激情在线99| bbb黄色大片| 国产综合懂色| 丁香六月欧美| 免费无遮挡裸体视频| 成年版毛片免费区| 欧美又色又爽又黄视频| 欧美午夜高清在线| 午夜亚洲福利在线播放| 国产人妻一区二区三区在| www.熟女人妻精品国产| 亚洲午夜理论影院| av天堂中文字幕网| 18禁黄网站禁片免费观看直播| 久久久久久久久久黄片| 国产黄色小视频在线观看| 久久午夜福利片| 精品99又大又爽又粗少妇毛片 | 国产成人欧美在线观看| 丰满人妻熟妇乱又伦精品不卡| 波多野结衣巨乳人妻| 国产乱人视频| 久久精品国产亚洲av香蕉五月| 我要搜黄色片| 亚洲成人久久爱视频| 国产av一区在线观看免费| 国产精品99久久久久久久久| ponron亚洲| 亚洲精品日韩av片在线观看| 国产一区二区在线av高清观看| 性色av乱码一区二区三区2| 国产真实乱freesex| 欧洲精品卡2卡3卡4卡5卡区| 99国产精品一区二区三区| 亚洲av成人不卡在线观看播放网| 免费电影在线观看免费观看| 欧美一区二区亚洲| 午夜福利视频1000在线观看| 久久国产精品影院| 在线观看舔阴道视频| 国内少妇人妻偷人精品xxx网站| 欧美午夜高清在线| 色精品久久人妻99蜜桃| 特大巨黑吊av在线直播| 国产伦一二天堂av在线观看| 国产爱豆传媒在线观看| 国产国拍精品亚洲av在线观看| 一边摸一边抽搐一进一小说| 丝袜美腿在线中文| 久久草成人影院| 欧美+日韩+精品| x7x7x7水蜜桃| 十八禁国产超污无遮挡网站| av福利片在线观看| 色综合亚洲欧美另类图片| 丝袜美腿在线中文| 91麻豆av在线| 欧美+日韩+精品| 久久欧美精品欧美久久欧美| 一二三四社区在线视频社区8| 男人的好看免费观看在线视频| 色av中文字幕| 国产成+人综合+亚洲专区| 日本五十路高清| 亚洲精品粉嫩美女一区| 久99久视频精品免费| 校园春色视频在线观看| 国产熟女xx| 两个人的视频大全免费| 亚洲精品色激情综合| 国产精品嫩草影院av在线观看 | 亚洲美女黄片视频| 免费大片18禁| 能在线免费观看的黄片| 亚洲在线自拍视频| 亚洲一区二区三区色噜噜| 十八禁网站免费在线| 一本一本综合久久| 亚洲三级黄色毛片| 人妻夜夜爽99麻豆av| 哪里可以看免费的av片| 亚洲成av人片免费观看| 天堂网av新在线| 一进一出抽搐gif免费好疼| 午夜影院日韩av| 国产黄色小视频在线观看| 亚洲美女黄片视频| 欧洲精品卡2卡3卡4卡5卡区| 黄色视频,在线免费观看| 夜夜爽天天搞| 一本一本综合久久| av在线老鸭窝| 麻豆成人av在线观看| 亚洲激情在线av| 亚洲五月婷婷丁香| 免费在线观看日本一区| 国产精品一区二区性色av| 黄片小视频在线播放| 久久天躁狠狠躁夜夜2o2o| 国内久久婷婷六月综合欲色啪| 亚洲av成人精品一区久久| 99久久九九国产精品国产免费| 久久精品久久久久久噜噜老黄 | 无遮挡黄片免费观看| 亚洲无线在线观看| 国产成人av教育| netflix在线观看网站| 国产三级在线视频| 免费av毛片视频| 99riav亚洲国产免费| 午夜福利成人在线免费观看| 婷婷亚洲欧美| 久久久成人免费电影| 亚洲成av人片免费观看| 国产探花在线观看一区二区| 大型黄色视频在线免费观看| 有码 亚洲区| 一级黄片播放器| 久久久久免费精品人妻一区二区| 亚洲欧美日韩无卡精品| 一个人免费在线观看电影| 久久久国产成人精品二区| 一区二区三区激情视频| АⅤ资源中文在线天堂| 中文字幕精品亚洲无线码一区| 亚洲欧美日韩东京热| 亚洲综合色惰| 国产真实乱freesex| 伊人久久精品亚洲午夜| 宅男免费午夜| 精品熟女少妇八av免费久了| 长腿黑丝高跟| 免费黄网站久久成人精品 | 日本精品一区二区三区蜜桃| 成年人黄色毛片网站| 中文资源天堂在线| 色噜噜av男人的天堂激情| 国产精华一区二区三区| 国产精品人妻久久久久久| 精华霜和精华液先用哪个| 毛片女人毛片| 91狼人影院| 免费电影在线观看免费观看| 在线看三级毛片| 变态另类成人亚洲欧美熟女| 国产av不卡久久| 日本黄色片子视频| 国产白丝娇喘喷水9色精品| 国产欧美日韩精品亚洲av| 久久久久性生活片| 精品不卡国产一区二区三区| 国产精品三级大全| 国产精品久久久久久久久免 | 色噜噜av男人的天堂激情| 亚洲精品久久国产高清桃花| 午夜免费男女啪啪视频观看 | bbb黄色大片| 国产毛片a区久久久久| 欧美黑人欧美精品刺激| 欧美色欧美亚洲另类二区| 欧美不卡视频在线免费观看| 女生性感内裤真人,穿戴方法视频| 欧美中文日本在线观看视频| 女同久久另类99精品国产91| 欧美性猛交黑人性爽| 特大巨黑吊av在线直播| 午夜免费成人在线视频| 日本免费a在线| netflix在线观看网站| 一二三四社区在线视频社区8| xxxwww97欧美| 欧美不卡视频在线免费观看| 色哟哟·www| 天堂影院成人在线观看| 精品人妻1区二区| 亚洲精华国产精华精| 国产精品人妻久久久久久| 国产蜜桃级精品一区二区三区| 久9热在线精品视频| 国产午夜精品论理片| 亚洲自偷自拍三级| 禁无遮挡网站| 国产激情偷乱视频一区二区| 日韩欧美一区二区三区在线观看| 精品日产1卡2卡| 亚洲专区中文字幕在线| 色精品久久人妻99蜜桃| 午夜免费男女啪啪视频观看 | 久久精品国产清高在天天线| 男女做爰动态图高潮gif福利片| 日本黄色视频三级网站网址| 91午夜精品亚洲一区二区三区 | 国产精品伦人一区二区| 一本综合久久免费| 国产主播在线观看一区二区| 日本撒尿小便嘘嘘汇集6| 国产精品99久久久久久久久| 亚洲自偷自拍三级| 熟女人妻精品中文字幕| 少妇高潮的动态图| 天堂动漫精品| 亚洲av不卡在线观看| 黄色日韩在线| 精品午夜福利在线看| 首页视频小说图片口味搜索| 我的老师免费观看完整版| 噜噜噜噜噜久久久久久91| 久久精品国产亚洲av天美| 免费在线观看影片大全网站| 亚洲精品影视一区二区三区av| 一级黄片播放器| 日本黄色视频三级网站网址| 亚洲精品亚洲一区二区| 美女cb高潮喷水在线观看| 99久久成人亚洲精品观看| 少妇人妻一区二区三区视频| 最近中文字幕高清免费大全6 | 色噜噜av男人的天堂激情| 欧美成人a在线观看| 欧美日韩国产亚洲二区| 婷婷六月久久综合丁香| 亚洲精品456在线播放app | 搡老熟女国产l中国老女人| 国内揄拍国产精品人妻在线| 久久99热6这里只有精品| 国产亚洲精品久久久久久毛片| 搡老妇女老女人老熟妇| 国内久久婷婷六月综合欲色啪| 亚洲 国产 在线| 桃红色精品国产亚洲av| 久久久久免费精品人妻一区二区| 亚洲内射少妇av| 亚洲欧美精品综合久久99| 欧美又色又爽又黄视频| 色播亚洲综合网| 99久久精品一区二区三区| 国产精品影院久久| 亚洲精品久久国产高清桃花| 亚洲va日本ⅴa欧美va伊人久久| www.999成人在线观看| 久久伊人香网站| 国产免费男女视频| 乱人视频在线观看| 波多野结衣巨乳人妻| 黄色女人牲交| 亚洲美女搞黄在线观看 | 少妇裸体淫交视频免费看高清| 99久久成人亚洲精品观看| 亚洲片人在线观看| 欧美中文日本在线观看视频| 精品一区二区三区视频在线观看免费| 91久久精品国产一区二区成人| 天堂av国产一区二区熟女人妻| 久久久久久久久久黄片| 亚洲第一欧美日韩一区二区三区| 中文在线观看免费www的网站| 国产三级在线视频| 国内精品久久久久久久电影| 国产亚洲精品久久久久久毛片| 久久亚洲真实| 怎么达到女性高潮| 精品乱码久久久久久99久播| 亚洲第一电影网av| 亚洲精品亚洲一区二区| av女优亚洲男人天堂| 免费无遮挡裸体视频| 亚洲性夜色夜夜综合| 99热这里只有是精品50| 国产伦精品一区二区三区视频9| 国产精品美女特级片免费视频播放器| 亚洲欧美日韩东京热| 一二三四社区在线视频社区8| 久久九九热精品免费| 国产亚洲精品综合一区在线观看| 色综合站精品国产| 亚洲片人在线观看| 久久天躁狠狠躁夜夜2o2o| 久久久久免费精品人妻一区二区| 丰满乱子伦码专区| 久久伊人香网站| 久久精品国产亚洲av香蕉五月| 免费观看的影片在线观看| h日本视频在线播放| 色综合站精品国产| 人人妻人人看人人澡| 久久久精品欧美日韩精品| 九色国产91popny在线| 国产精品免费一区二区三区在线| 国产精品伦人一区二区| 成人三级黄色视频| 久久久久性生活片| 国产男靠女视频免费网站| 久久久国产成人免费|