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

    基于離心機法獲取定體積質量下的土壤水分特征曲線

    2020-03-20 06:35:42任健張吳平王國芳黃明鏡喬磊
    灌溉排水學報 2020年1期
    關鍵詞:特征質量模型

    任健,張吳平,王國芳,黃明鏡,喬磊

    基于離心機法獲取定體積質量下的土壤水分特征曲線

    任健1,張吳平2*,王國芳1,黃明鏡3,喬磊1

    (1.山西農業(yè)大學 資源環(huán)境學院,山西 太谷 030800;2.山西農業(yè)大學 軟件學院,山西 太谷 030800;3.山西省農業(yè)科學研究院 旱地農業(yè)研究中心,太原 030000)

    【】克服離心機法測定土壤水分特征曲線時“變體積”的缺陷,實現用離心機法獲取擾動土壤準確的土壤水分特征曲線。選取了山西省16種粒徑組成、有機質量差別顯著的擾動土壤作為研究對象,采用離心機法測定了初始體積相同的各土壤樣品在不同轉速下的土壤體積含水率、動態(tài)變化的土壤體積;采用多元非線性回歸模型,描述了以土壤體積含水率、粒徑組成、體積質量和有機質質量分數作為自變量描述土壤水吸力的函數關系,并對模型進行了驗證與評價。土壤水吸力的自然對數函數可以表示為土壤砂粒質量分數、黏粒質量分數、土壤體積質量和含水率乘積以及有機質質量分數等土壤屬性的線性函數,模型的決定系數達到0.740;模型建模樣本平均相對誤差為14.460%,驗證樣本平均相對誤差為13.237%,模型預測精度較高;同時發(fā)現土壤屬性對于曲線的形狀與走勢的影響不能忽視,其中砂粒質量分數、土壤體積質量和含水率乘積與土壤水吸力負相關,黏粒質量分數、有機質質量分數與土壤水吸力正相關。因此,本方法構建的模型更符合實際情況,能夠提高土壤水力學參數的精度,同時可以克服離心機測量土壤水分特征曲線體積變化的缺點,得到更合理的土壤水吸力和土壤含水率的關系。

    土壤水分特征曲線;多元非線性預測模型;離心機法;仿真數值模擬

    0 引 言

    土壤水分特征曲線是土壤含水率與土壤水吸力的函數,反映了土壤數量與能量之間的關系,是獲取土壤非飽和導水率、定量模擬土壤水分溶質運移及科學灌溉的重要參數,在科學研究和生產實踐中有重要意義[1-5]。

    目前對于土壤水分特征曲線的測定包括間接推求法和直接測定法,一些研究通過土壤顆粒分布曲線、土壤體積質量、質地、田間持水率和萎蔫含水率等土壤屬性來間接地推求土壤水分特征曲線[6-9]。雖然間接推求法成本較低,但是其精度不高,在生產實踐當中難以推廣[10]。

    當前對于農田尺度的研究更多的是采用直接測定法,常用的方法有張力計法、離心機法、壓力膜儀法等。其中,離心機法操作較簡便,測得的吸力范圍較廣,省時省力,土壤中的水分更易排出[11]。但是土樣在離心的過程中,體積會越來越小,因此采用該方法獲取的土壤水分特征曲線是在“變體積”情況下測得的。呂殿青等[12]發(fā)現武功粉砂質黏土在用離心機測定土壤水分特征曲線時,體積質量在1.02~100.2 m水柱范圍內從1.398 g/cm3增加到了1.858 g/cm3,變幅達到了0.565 g/cm3;聶坤堃等[11]發(fā)現用離心機法獲取土壤收縮特性時,擾動土樣初始設計體積質量分別為1.30、1.35、1.40 g/cm3的情況下,體積質量變幅分別達到了37.69%、33.33%、28.57%。邵明安[13]認為土壤持水性質不是一條單一的曲線,而是由土壤含水率、體積質量、吸力共同組成的一個空間曲面,可見體積質量對于土壤水分特征曲線的影響不能夠被忽視。

    為了克服離心機“變體積”的缺陷,呂殿青[14]采用了“石蠟法”獲取了定體積質量下的土壤水分特征曲線,通過設定不同離心機轉速使土樣達到不同轉速下的平衡狀態(tài),并在平衡狀態(tài)的土樣中注入了熔化的石蠟,待石蠟冷凝后會使土壤體積固定,從而得到定體積質量下的土壤水分特征曲線,但此方法只能測定土樣平衡狀態(tài)下該體積質量下的定體積質量曲線,如果要測定不同體積質量梯度下的曲線還需多次設定離心機轉速來獲得不同平衡狀態(tài)、不同初始體積質量的土樣,并且體積質量不易控制,因此“石蠟法”有自身的局限性,操作起來較為不便。

    除此之外,離心機法獲取的土壤含水率范圍有限。由于離心機轉子的轉速有限,因此不能夠獲取在極低含水率下的土壤水吸力值。同時,相關研究指出,土壤水分特征曲線除了受到體積質量變化的影響外,還受到土壤粒徑組成、有機質質量分數的影響[15-16]。已有研究表明土壤質地、有機質質量分數對土壤水分特征曲線的影響極大,且與Van Genuchten模型當中的參數相關性極高[15,17]。

    基于此,本研究擬采用離心機法直接獲取不同粒徑組成、有機質質量分數土壤的“變體積”下的土壤水吸力和含水率。根據實測結果,采用統(tǒng)計回歸分析模型,構建以土壤水吸力對數為因變量,土壤體積質量、含水率、粒徑組成、有機質及各因素變量的兩兩乘積作為自變量,建立土壤水分特征曲線多元非線性預測模型;采用仿真數值分析方法克服“石蠟法”所帶來的缺陷,排除體積質量對其曲線的影響,克服離心機法的缺陷,并同時考慮粒徑組成和有機質對土壤水分特征曲線的定量影響,其建模結果可以自由設定土壤體積質量的值,從而最終獲取“定體積質量”下的土壤水分特征曲線,并對所建預測模型進行驗證與評價分析。以作為黃土高原區(qū)開展的旱作農業(yè)中土壤水分運動的定量模擬和預測,調節(jié)土壤水分有效性及儲水量,對指導旱作農業(yè)的灌溉與排水具有重要意義。

    1 材料與方法

    1.1 供試土壤

    此次研究選擇山西省16種不同質地、有機質質量分數的土樣。在每個采樣區(qū)選取代表性的地塊將部分擾動土用塑封袋帶回。將采集到的擾動土樣帶入室內自然風干,研磨過2 mm篩,作為此次離心機當中的供試土樣,并按照初始體積質量為1.3 g/cm3均勻裝入直徑為5 cm、高度為3 cm的環(huán)刀中。采用烘干法測定土樣殘余含水率。根據殘余含水率的數值可得到土樣中干土的質量。將裝有擾動土的環(huán)刀浸泡于蒸餾水中24 h,水面高度低于環(huán)刀高度1 mm,之后將環(huán)刀滯空放置一段時間,直到不再有重力水排出,此時認為土樣達到飽和,并用高精度電子天平稱質量計算獲取土樣飽和含水率。土壤質地的測定采用比重計法,并按照國際制土壤質地分類為分級標準。土壤有機質質量分數的測定采用濃硫酸-重鉻酸鉀外加熱法。

    對山西省土壤16種基本性質的測定,如表1所示,所選的土壤質地與有機質質量分數差異變化較大,砂粒的變化范圍是5%~60%,黏粒變化范圍是10%~60%,有機質質量分數變化范圍是0.248%~2.475%,此變化范圍可以代表山西省土壤質地、有機質質量分數的變化,從而大大減少工作量,并且仍能得到試驗對象主要特征的試驗結果。

    1.2 土壤水分特征曲線的獲取

    土壤水分特征曲線的獲取采用離心機法。離心機采用日立CR22N型高速恒溫冷凍離心機,并附有土壤水分特征曲線定制轉子、土樣盒、環(huán)刀與排水裝置。定制轉子最高轉速可達11 000 rpm,最大能夠同時裝填4個離心盒,并具有轉頭自動識別與轉頭自鎖功能,保證在離心過程中土樣盒保持相對靜止。試驗通過設置離心機不同的轉速梯度從而計算得到不同轉速下對應的土壤水吸力,為了保證試驗的精度,各研究區(qū)每層土壤均設置3次重復,每次設置9個不同轉速,分別為300、600、1 000、1 600、2 500、3 000、3 500、6 500、11 000 rpm,對應9個不同土壤水吸力。待每個轉速下土樣平衡之后,取出離心盒,靜置一段時間后用高精度電子天平稱量土樣質量[18],用游標卡尺測量土樣高度,從而計算土樣平衡時的質量與體積。根據計算出的土壤體積含水率和所對應的土壤水吸力,即可繪制出土壤水分特征曲線[19]。

    1.3 土壤水吸力的計算

    邢旭光等[20]在研究土壤水分特征曲線測定過程中土壤收縮特性時,土壤水吸力計算式為:

    式中:為土壤水吸力(cm);1為離心機轉子軸心至土樣中心之間的距離(cm);2為離心機轉子軸心至土樣底部的距離(cm);w為水的密度(g/cm3);g為重力加速度(cm/s2);為角速度(rad/s)。

    1.4 土壤體積質量的計算

    離心機中的土樣達到平衡狀態(tài)之后,根據游標卡尺測得的土樣高度可計算體積,根據殘余含水率可以計算土樣的干土質量,計算式為:

    =/[(r+1)π2],(2)

    式中:為土樣質量(g);為體積質量(g/cm3);r為殘余含水率(%);為環(huán)刀內半徑(cm)。

    1.5 模型精度評估的統(tǒng)計學方法

    本次研究結合16種土壤的建模樣本數據與另外23組驗證樣本數據,驗證樣本數據來自山西省太谷縣孟家莊村表層擾動土壤共23個不同轉速下的土壤水分特征曲線數據,同建模樣本數據的獲取方式相同。通過比較模型預測值與實測值的相對誤差值、驗證樣本數據的均方根誤差()、決定系數(2)、誤差比的幾何平均數()來分析建模精度[21]:

    表1 供試土壤基本理化性質

    2 結果與分析

    2.1 多元非線性預測模型的建立與驗證

    2.1.1 多元非線性預測模型的建立

    多元非線性模型是反映自變量與因變量之間的非線性關系的數學表達式。離心機在獲取土壤水分特征曲線時體積質量和含水率在不斷變化,將體積質量與含水率的乘積作為自變量,同時將16種供試土壤不同梯度的砂粒質量分數、黏粒質量分數、有機質質量分數作為自變量,通過回歸分析預測土壤水吸力,建立如下非線性關系的預測模型:

    ln=11.073-0.041+0.057+0.609-14.203,(7)

    式中:為土壤砂粒質量分數(%);為土壤黏粒質量分數(%);為土壤有機質質量分數(%);為土壤體積含水率(cm3/cm3);其余參數同前。該模型決定系數為0.740,相關性較高。

    2.1.2 多元非線性預測模型自變量的檢驗

    對建立該預測模型的16組數據進行檢驗,通過SPSS回歸結果進行分析,結果見表2。通過Excel軟件中的tinv函數計算得到(0.05,143)的值為1.977。由表2結果得知,4項參數的統(tǒng)計量值的絕對值全部大于1.977,值全部小于0.05,均通過檢驗,說明具有統(tǒng)計學意義,最終可確定式(7)所建預測模型。

    表2 非線性預測模型t值檢驗表

    2.1.3 多元非線性預測模型檢驗

    檢驗實際為對回歸關系的方差分析,對建立該預測模型的16組數據進行顯著性檢驗,通過SPSS軟件中的回歸統(tǒng)計模塊求出值,并比較顯著性水平為0.05下的值,結果見表3。

    表3 非線性預測模型F值檢驗表

    由表3可知,回歸分析建立的非線性預測模型通過了檢驗,檢驗統(tǒng)計量值為99.928,大于顯著性水平為0.05下的值,表明該模型的土壤水吸力與其余4個自變量具有統(tǒng)計學上極顯著的意義,是有效的。

    2.1.4 預測模型的驗證

    上述模型的建立共16條土壤水分特征曲線數據,每條曲線含有9個不同梯度的含水率—吸力數據,總共144個樣本數據,做出此144個樣本吸力對數的實測值、預測值、相對誤差值。由圖1可知,建模樣本相對誤差平均值為14.460%,精度較高,說明用離心機法結合土壤基本屬性推求吸力預測模型是可行的。

    另外做出23個驗證樣本吸力對數的實測值、預測值、相對誤差值如圖2所示。由圖2可知,驗證樣本的平均相對誤差平均值為13.237%,其精度與建模樣本的精度相當。無論在模型的建立階段和驗證階段,其相對誤差平均值都較低,說明該預報模型精度較高,具有較強的適用性,可以用于土壤水分特征曲線的預測。

    將23組驗證數據代入預測模型中,對比分析預測吸力對數與實測吸力對數,并根據式(3)—式(6)計算各統(tǒng)計量,結果繪于圖3中。由圖3可知,預測值與實測值比較接近,基本均勻分布在對稱線左右,說明模型預測精度良好;預測值與實測值的決定系數為0.958,說明預測值與實測值的變化趨勢基本一致;為1.10,該值大于1,說明預測值較實測值偏大;的值越小,表明預測值與真值的偏差越小,精度更高,其值為0.67,表明預測模型的預測值與實測值吻合度較高。

    2.2 根據預測模型模擬各種參數變化下的土壤水分特征曲線

    2.2.1 定體積質量下的土壤水分特征曲線的模擬與分析

    土壤體積質量的變化影響土壤孔隙度[22],且還會影響土壤孔隙分布,從而改變了土壤的持水特性,勢必會影響土壤水分特征曲線。采用仿真數值模擬的方法,將砂粒黏粒配比、有機質質量分數設定為定值,體積質量設定為1.1~1.7 g/cm3,即可獲取不同體積質量下土壤水分特征曲線(圖4)。由圖4可知,隨著土壤體積質量的增大,土壤水分特征曲線呈現出向左移動的趨勢,且坡度更小;在土壤含水率相同的情況下,高體積質量的土壤表現為更低的土壤水吸力,也就意味著在離心過程中受到較小的離心力的作用即可排出土壤中殘留的水分;在相同的吸力下,體積質量較大的土壤所持有的土壤含水率較低。

    2.2.2 粒徑變化對土壤水分特征曲線的影響分析

    土壤質地的差異會改變土壤水分特征曲線的特性,砂粒質量分數高的土壤,其通透性良好,飽和導水率較高,導致土壤對水分的保蓄能力差,在相同的土壤含水率下,砂粒質量分數越高,土壤水吸力越小。在所建預測模型當中,砂粒質量分數的系數為負,即砂粒質量分數與土壤水吸力為負相關。黏粒質量分數高的土壤,土壤的通透性差,飽和導水率較低,土壤對水分的保蓄能力強。由圖4可知,在相同的含水率下,土壤水分特征曲線右移,黏粒質量分數越高,土壤水吸力越大,曲線更為陡峭。這驗證了預測模型黏粒質量分數與土壤水吸力的正相關關系。

    2.2.3 有機質變化對土壤水分特征曲線的影響分析

    土壤有機質是土壤固相部分的重要組成,疏松多孔的性質使得能夠改善土壤的物理性質,促進微生物和土壤生物的活動,大大提高了土壤儲水性質,從而改變土壤水分特征曲線的走勢。通過此次仿真數值模擬得知,土壤有機質質量分數越高,土壤水分特征曲線的走勢越陡,說明在有機質質量分數高的土壤中,土壤吸附能力更強,這與舒凱民[15]的相關研究結果是一致的。

    圖1 建模樣本相對誤差圖

    圖2 驗證樣本相對誤差圖

    圖3 吸力對數預測值與實測值對比

    3 討 論

    通過離心機法直接獲取的“變體積質量”下土壤水分特征曲線實測結果,并利用回歸分析方法建立其多元非線性預報模型,最后采用仿真數值模擬方法得到“定體積質量”下的土壤水分特征曲線,能夠完全實現“石蠟法”的效果[14],克服了離心過程中“變體積質量”的缺陷。土壤水分特征曲線的多元非線性預測模型相比于Van Genuchten模型而言更加直接。Van Genuchten模型需要眾多系數,包括土壤形狀系數及進氣值系數,這些參數較難獲得;而本文提出的模型只需要土壤基本的屬性即可推求高精度的土壤水分特征曲線,更方便用于土壤水分運動與農田科學合理灌溉的相關研究。

    然而,用離心機法獲取土壤水分特征曲線時,由于轉子轉速有限、極低的土壤含水率對應的土壤水吸力非常大,結果導致缺失土壤含水率極低情況下對應的土壤水吸力數據。在這16種不同的實測土壤水分特征曲線當中,其最大吸力對應的體積含水率最小值為0.083 cm3/cm3,最大值為0.274 cm3/cm3,平均值為0.160 cm3/cm3。對于大田耕作土壤而言,大部分農作物適墑的土壤質量含水率范圍在12%~30%之間[2],換算為體積含水率之后其二者數值相當,可見離心機法所獲取的土壤含水率范圍與大田耕作土壤含水率基本吻合,缺失的那一段曲線基本可以忽略,因此用離心機法獲取土壤水分特征曲線的方法是可行的。

    該預報模型建立了土壤多種屬性影響下的水分特征曲線,土壤體積質量、質地、有機質質量分數都會顯著影響土壤水分特征曲線的形狀與走勢。多元非線性預報模型的建立表明了通過土壤質地、有機質質量分數、體積質量、含水率即可反映出土壤水吸力,其相關系數為0.740,說明該模型能夠預測出74%的土壤水分特征曲線。相比于直接方法獲取土壤水分特征曲線,該方法更為便捷。該預測模型可看作土壤水吸力受砂粒質量分數、黏粒質量分數、有機質質量分數、體積質量與含水率的乘積4個因素的決定與影響,各項因素的系數可表達出對模型結果的波動狀況,其中體積質量與含水率乘積的系數明顯大于其余系數,表明體積質量與含水率的乘積對土壤水吸力的波動之大,更堅定的說明了在用離心機法獲取土壤水分特征曲線不能夠忽視“變體積質量”對其的影響。

    4 結 論

    本文所建預測模型相關系數為0.740,通過了檢驗與檢驗,建模樣本、驗證樣本的相對誤差平均值分別為14.460%、13.237%,值為0.67,所建立的預測模型精度較高;通過仿真數值模擬的方法,得出基于離心機法下獲取的定體積質量土壤水分特征曲線,并隨著體積質量的減小,曲線有向上移動的趨勢;同時也定量表示了預測模型中土壤質地、有機質質量分數變化下的土壤水分特征曲線。

    [1] 邵明安, 王全九, 黃明斌. 土壤物理學[M]. 北京: 高等教育出版社, 2006: 67-69.

    SHAO Mingan, WANG Quanjiu, HUANG Mingbin. Soil Physics[M]. Beijing: Higher Education Press, 2006: 67-69.

    [2] 龔振平, 邵孝侯, 張富倉, 等. 土壤學與農作學[M]. 北京: 中國水利水電出版社, 2009.

    GONG Zhenping, SHAO Xiaohou, ZHANG Fucang, et al. Soil science and agronomy[M]. Beijing: China Water Power Press, 2009.

    [3] 車政, 王仰仁, 王永紅, 等. 農田土壤水分特征曲線參數擬合及其剖面變異特性研究[J]. 灌溉排水學報, 2016, 35(7): 22-27.

    CHE Zheng, WANG Yangren, WANG Yonghong, et al. Research on Parameter Fitting of Soil Water Characteristic Curve and Its Variation Property with Depth of Farmland Soil Profile[J]. Journal of Irrigation and Drainage, 2016, 35(7): 22-27.

    [4] 于博, 于曉芳, 高聚林, 等. 秸稈全量深翻還田和施加生物炭對不同土壤持水性的影響[J]. 灌溉排水學報, 2018, 37(5): 25-32.

    YU Bo, YU Xiaofang, GAO Julin, et al. Change in Hydraulic Properties of Soils Amended with Biochar Following Plough of Straw Stalk into Deep Soil Horizon[J]. Journal of Irrigation and Drainage, 2018, 37(5): 25-32.

    [5] 陳安強, 雷寶坤, 胡萬里, 等. 洱海近岸菜地包氣帶土壤水分特征曲線參數變化及其影響因素[J]. 灌溉排水學報, 2018, 37(10): 48-54.

    CHEN Anqiang, LEI Baokun, HU Wanli, et al. Water Retention Curve of Vegetative Soils in Erhai Lake Catchment and the Factors Affecting It[J]. Journal of Irrigation and Drainage, 2018,37(10):48-54.

    [6] AHUJA L R, NANEY J W, WILLIANMS R D. Estimating soil water retention form simpler properties or limited data[J]. Soil Science Socciety of America Journal, 1985, 49: 1 100-1 105.

    [7] MINASNY B, MCBRATNEY A B, BRISTOW K L. Comparison of different approaches to the development of pedotransfer functions for water-retention curves[J]. Geoderma, 1999, 93(3/4): 225-253.

    [8] NAVEED M, MOLDRUP P, TULLER M, et al. Prediction of the soil water characteristic from soil particle volume fractions[J]. Soil Science Society of America Journal, 2012, 76(6): 1 946-1 956.

    [9] RAJKAI K, KABOS S, VAN Genuchten M T. Estimating the water retention curve from soil properties: Comparison of linear, nonlinear and concomitant variable methods[J]. Soil and Tillage Research, 2004, 79(2): 145-152.

    [10] 呂殿青, 邵明安. 非飽和土壤水力參數的模型及確定方法[J]. 應用生態(tài)學報, 2004, 15(1): 163-166.

    LV Dianqing, SHAO Mingan. Models and determining methods of unsaturated soil hydraulic parameters[J]. Chinese Journal of Applied Ecology, 2004, 15(1): 163-166.

    [11] 聶坤堃, 聶衛(wèi)波, 馬孝義. 離心機法測定土壤水分特征曲線中的收縮特性[J]. 排灌機械工程學報, 2018(3): 1-2.

    NIE Kunkun, NIE Weibo, MA Xiaoyi. Characteristic of soil shrinkage during soil water characteristic curve measurement by the centrifuge method[J]. Journal of Drainage and Irrigation Mechanical Engineering, 2018(3): 1-2.

    [12] 呂殿青, 邵明安, 王全九. 土壤持水特征測定中的容重變化及其確定方法[J]. 水利學報, 2003(3): 110-114.

    LV Dianqing, SHAO Mingan, WANG Quanjiu. Application of fuzzy clustering method in designing fuzzy neural exciting controller of generator[J]. Journal of Hydraulic Engineering, 2003(3): 110-114.

    [13] 邵明安, 呂殿青, 付曉莉, 等. 土壤持水特征測定中質量含水量、吸力和容重三者間定量關系I. 填裝土壤[J]. 土壤學報, 2007, 44(6): 1 003-1 009.

    SHAO Mingan, LV Dianqing, FU Xiaoli, et al. Quantitative relationship between mass water content, pressure head and bulk density in determination of soil water retention characteristics i. packed soils[J]. Acta Pedologica Sinica, 2007, 44(6): 1 003-1 009.

    [14] 呂殿青, 邵明安, 潘云. 容重變化與土壤水分特征的依賴關系研究[J]. 水土保持學報, 2009, 23(3): 209-212.

    LV Dianqing, SHAO Mingan, PAN Yun. Dependent Relationship Between Bulk Density Changes and Soil Water Characteristics[J]. Journal of Soil and Water Conservation, 2009, 23(3): 209-212.

    [15] 舒凱民. 基于土壤基本理化參數的土壤水分特征曲線Van-Genuechten模型預報研究[D]. 太原: 太原理工大學, 2017.

    SHU Kaiming. Research on Prediction of Soil Water Characteristic Curve Van-Genuechten Model Based on Soil Basic Physical and Chemical Parameters[D]. Taiyuan: Taiyuan University of Technology, 2017.

    [16] 楊緒, 趙廷超, 王豐, 等. 基于Gardner模型的黔西南土壤水分特征分析[J]. 灌溉排水學報, 2018, 37(7): 35-42.

    YANG Xu, ZHAO Tingchao, WANG Feng, et al. Water Retention Characteristics of Typical Soils in Southwestern Guizhou Province[J]. Journal of Irrigation and Drainage, 2018, 37(7): 35-42.

    [17] 舒凱民, 樊貴盛. 基于質地的土壤水分特征曲線參數非線性預測[J]. 人民黃河, 2016, 38(7): 138-141.

    SHU Kaiming, FAN Guisheng. The Nonlinear Prediction of Soil Water Characteristic Curve Parameters Based on the Soil Texture[J]. Yellow River, 2016, 38(7): 138-141.

    [18] 劉暢, 陳俊英, 蔡耀輝, 等. 斥水土壤水分特征曲線及擬合模型分析[J]. 中國農村水利水電, 2017(7): 45-48.

    LIU Chang, CHEN Junying, CAI Yaohui, et al. Soil Water Characteristic Curves of Different Kinds of Repellency Soil and an Analyses of the Fitting Model[J]. China Rural Water and Hydropower, 2017(7): 45-48.

    [19] 牛曉彤, 劉目興, 易軍, 等. 三峽山地土壤水分特征曲線及模型擬合[J]. 灌溉排水學報, 2017, 36(9): 75-80.

    NIU Xiaotong, LIU Muxing, YI Jun, et al. Water Release Curves and Their Fitting for Soils in the Hills of the Three Gorges Reservoir[J]. Journal of Irrigation and Drainage, 2017, 36(9): 75-80.

    [20] 邢旭光, 趙文剛, 馬孝義, 等. 土壤水分特征曲線測定過程中土壤收縮特性研究[J]. 水利學報, 2015, 46(10): 1 181-1 188.

    XING Xuguang, ZHAO Wengang, MA Xiaoyi, et al. Study on soil shrinkage characteristics during soil water characteristic curve measurement[J]. Journal of Hydraulic Engineering, 2015, 46(10):1 181-1 188.

    [21] 姚姣轉, 劉廷璽, 王天帥, 等. 科爾沁沙地土壤水分特征曲線傳遞函數的構建與評估[J]. 農業(yè)工程學報, 2014, 30(20): 98-108.

    YAO Jiaozhuan, LIU Tingxi, WANG Tianshuai, et al. Development and evaluation of pedo-transfer functions of soil water characteristic curves in Horqin sandy land[J]. Transactions of the Chinese Society of Agricultural Engineering, 2014, 30(20): 98-108.

    [22] 洪成, 尹殿勝, 陳俊英, 等. 容重對黏壤土土壤水分特征曲線的影響[J]. 節(jié)水灌溉, 2018(10): 5-8.

    HONG Cheng, YIN Diansheng, CHEN Junying, et al. Impact of Bulk Density on Soil Water Characteristic Curve of Clay Loam[J]. Water Saving Irrigation, 2018(10): 5-8.

    Method to Adjust the Soil Water Retention Curve Measured Using Centrifuge Apparatus

    REN Jian1, ZHANG Wuping2*, WANG Guofang1, HUANG Mingjing3, QIAO Lei1

    (1.Shanxi Agricultural University, College of Resources and Environment, Taigu 030800, China;2.Shanxi Agricultural University, Software Institute, Taigu 030800, China;3. Arid Farming Research Center, Shanxi Academy of Agricultural Sciences, Taiyuan 030000, China)

    【】Centrifuge is a common apparatus to measure soil water retention curve, but the resultant soil compression is an issue that might compromise its accuracy and reliability. The purpose of this paper it to present a method to correct such inferiorities.【】We measured the water retentions of 16 repacked soils with contrasting texture and organic matter (OM) taken from Shanxi province using a centrifuge apparatus. In the experiment, we measured the initial soil volume, dynamic soil volume, as well as water content under different spinning speeds. A multiple nonlinear regression model was derived to describe the dependence of the matric potential on volumetric moisture content, soil particle sizes, volumetric soil mass and OM; it was then verified against the experimental data. 【】The logarithm of the matric potential was linearly related to mass fractions of the sand, clay, OM and volume of the water. The accuracy of the model was 74%, with an average relative error of 14.460% for the modelled samples, and of 13.237% for the samples used for verification. It was also found that soil texture affected predictivity of the model. The matric potential was negatively correlated to the mass fraction of sand, volumetric soil mass and water content, but positively correlated to the mass fraction of clay and organic matter.【】The nonlinear regression model linking the matric potential to soil water content, with the change in soil volume modified, was able to improve the accuracy of the centrifuge method; it has implications in soil hydraulic parameter measurement.

    soil water retention; multiple nonlinear regression model; centrifuge method; virtual numerical simulation

    S152.7+1

    A

    10.13522/j.cnki.ggps.2019236

    1672 - 3317(2020)01 - 0084 - 07

    2019-09-08

    山西省重點研發(fā)計劃重點項目(201703D211002-2);山西省科技攻關項目(20130311008-5)

    任?。?995-),男。碩士研究生,主要從事土地信息技術研究。E-mail:549484055@qq.com

    張吳平(1973-),男。教授,主要從事土地信息技術研究。E-mail:zwping@126.com

    任健, 張吳平, 王國芳, 等. 基于離心機法獲取定體積質量下的土壤水分特征曲線[J]. 灌溉排水學報, 2020, 39(1):84-90.

    REN Jian, ZHANG Wuping, WANG Guofang, et al. Method to adjust the soil water retention curve measured using centrifuge apparatus [J]. Journal of Irrigation and Drainage, 2020, 39(1): 84-90.

    責任編輯:陸紅飛

    猜你喜歡
    特征質量模型
    一半模型
    “質量”知識鞏固
    質量守恒定律考什么
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權M-估計的漸近分布
    如何表達“特征”
    做夢導致睡眠質量差嗎
    不忠誠的四個特征
    當代陜西(2019年10期)2019-06-03 10:12:04
    抓住特征巧觀察
    3D打印中的模型分割與打包
    a级毛片在线看网站| 国产欧美日韩综合在线一区二区| 国产成人午夜福利电影在线观看| 国产亚洲av片在线观看秒播厂| av女优亚洲男人天堂| 欧美精品av麻豆av| 日本vs欧美在线观看视频| 国产精品无大码| 亚洲欧洲国产日韩| 国产精品欧美亚洲77777| 久久精品久久久久久噜噜老黄| 国产免费福利视频在线观看| 亚洲成av片中文字幕在线观看 | 免费在线观看视频国产中文字幕亚洲 | 亚洲欧美日韩另类电影网站| 综合色丁香网| 色哟哟·www| 欧美人与性动交α欧美精品济南到 | 最近中文字幕高清免费大全6| 欧美人与性动交α欧美精品济南到 | 午夜免费观看性视频| 侵犯人妻中文字幕一二三四区| 丝袜美腿诱惑在线| 黄色配什么色好看| tube8黄色片| 国产精品99久久99久久久不卡 | 中文字幕人妻丝袜制服| 国产欧美亚洲国产| 成年人免费黄色播放视频| 五月开心婷婷网| 亚洲av欧美aⅴ国产| 最黄视频免费看| 又粗又硬又长又爽又黄的视频| 国产av国产精品国产| 人妻少妇偷人精品九色| 在线观看人妻少妇| 久久精品国产亚洲av涩爱| 成人毛片a级毛片在线播放| 看非洲黑人一级黄片| 亚洲av.av天堂| 久久久国产一区二区| 人人妻人人添人人爽欧美一区卜| 亚洲国产欧美网| 亚洲精品第二区| 乱人伦中国视频| 国产又色又爽无遮挡免| 成人毛片60女人毛片免费| 中文字幕人妻熟女乱码| 免费播放大片免费观看视频在线观看| 日韩精品有码人妻一区| 国产亚洲精品第一综合不卡| 久久久精品94久久精品| 人人妻人人澡人人看| 日韩人妻精品一区2区三区| 国产精品国产av在线观看| 国产欧美日韩一区二区三区在线| 精品少妇一区二区三区视频日本电影 | 可以免费在线观看a视频的电影网站 | 黄色毛片三级朝国网站| 丝袜脚勾引网站| 欧美日韩亚洲高清精品| 免费黄网站久久成人精品| 男女国产视频网站| freevideosex欧美| 亚洲精品aⅴ在线观看| 精品国产露脸久久av麻豆| 永久网站在线| 日本-黄色视频高清免费观看| 成年美女黄网站色视频大全免费| 香蕉丝袜av| 国产av码专区亚洲av| 成年美女黄网站色视频大全免费| 欧美精品人与动牲交sv欧美| 2022亚洲国产成人精品| 亚洲美女黄色视频免费看| 免费大片黄手机在线观看| 免费黄频网站在线观看国产| 国产亚洲av片在线观看秒播厂| 伊人久久国产一区二区| 久久人人97超碰香蕉20202| 欧美日韩亚洲高清精品| 久久久久久久国产电影| 欧美激情 高清一区二区三区| 黄频高清免费视频| 七月丁香在线播放| 久久国产亚洲av麻豆专区| 国产黄色视频一区二区在线观看| 日本欧美视频一区| 久久女婷五月综合色啪小说| 一区在线观看完整版| √禁漫天堂资源中文www| 久久韩国三级中文字幕| 久久国内精品自在自线图片| 毛片一级片免费看久久久久| 午夜福利视频精品| 一区二区三区乱码不卡18| 亚洲婷婷狠狠爱综合网| 美国免费a级毛片| 咕卡用的链子| 亚洲美女视频黄频| 美女国产视频在线观看| 国产一区二区在线观看av| 亚洲综合色网址| 久久午夜福利片| 飞空精品影院首页| 亚洲国产日韩一区二区| 黄片无遮挡物在线观看| 9色porny在线观看| 久久热在线av| 在线观看美女被高潮喷水网站| 亚洲国产欧美在线一区| 各种免费的搞黄视频| 又大又黄又爽视频免费| 国产极品粉嫩免费观看在线| 国产一区二区三区av在线| 大香蕉久久成人网| 久久99精品国语久久久| 欧美日韩视频精品一区| 国产成人精品无人区| 丝袜人妻中文字幕| 婷婷色综合大香蕉| 免费在线观看黄色视频的| 久久久精品免费免费高清| 久久久久国产精品人妻一区二区| 午夜精品国产一区二区电影| h视频一区二区三区| 在线天堂中文资源库| 成人亚洲精品一区在线观看| 国产av码专区亚洲av| 日韩av在线免费看完整版不卡| 少妇人妻 视频| 又粗又硬又长又爽又黄的视频| 黑人猛操日本美女一级片| 国产一区二区 视频在线| 毛片一级片免费看久久久久| 美女福利国产在线| 好男人视频免费观看在线| 卡戴珊不雅视频在线播放| 亚洲精品日本国产第一区| 一区二区三区激情视频| 国产精品不卡视频一区二区| 欧美日韩一区二区视频在线观看视频在线| av福利片在线| www日本在线高清视频| 精品亚洲成a人片在线观看| 国产男女超爽视频在线观看| 少妇的丰满在线观看| 性色avwww在线观看| 亚洲图色成人| 久久女婷五月综合色啪小说| 日韩精品有码人妻一区| 国产亚洲欧美精品永久| 大片电影免费在线观看免费| 国产精品二区激情视频| 一区福利在线观看| 亚洲视频免费观看视频| 国产淫语在线视频| 久久亚洲国产成人精品v| 久久av网站| 成年美女黄网站色视频大全免费| 日韩中文字幕视频在线看片| 天堂8中文在线网| 一级a爱视频在线免费观看| 久久这里有精品视频免费| 寂寞人妻少妇视频99o| av免费在线看不卡| 99久久精品国产国产毛片| 99久久人妻综合| 国产在线免费精品| 少妇人妻 视频| 久久韩国三级中文字幕| 十八禁高潮呻吟视频| 国产免费现黄频在线看| 亚洲美女视频黄频| 国产1区2区3区精品| 成年av动漫网址| 少妇的逼水好多| 午夜福利一区二区在线看| 亚洲欧美成人精品一区二区| 国产xxxxx性猛交| 国产在线视频一区二区| 国产福利在线免费观看视频| 久久久久久久久久人人人人人人| 成人亚洲精品一区在线观看| 亚洲一码二码三码区别大吗| 亚洲一码二码三码区别大吗| 国产精品成人在线| 桃花免费在线播放| 国产精品欧美亚洲77777| av天堂久久9| 精品卡一卡二卡四卡免费| 久久久久久人人人人人| 欧美日韩一区二区视频在线观看视频在线| 久久久久视频综合| 丝袜脚勾引网站| 久久精品aⅴ一区二区三区四区 | 成人毛片60女人毛片免费| 99热网站在线观看| 国语对白做爰xxxⅹ性视频网站| 丁香六月天网| 寂寞人妻少妇视频99o| 国产男女内射视频| 97精品久久久久久久久久精品| 亚洲人成网站在线观看播放| 亚洲av国产av综合av卡| 亚洲国产精品一区三区| 九九爱精品视频在线观看| 国产精品国产av在线观看| 亚洲av国产av综合av卡| 我的亚洲天堂| 久久精品国产亚洲av天美| 久久久久久久久久久免费av| 日韩中文字幕视频在线看片| 国产精品国产三级国产专区5o| 国产伦理片在线播放av一区| 精品一区二区三区四区五区乱码 | 欧美+日韩+精品| 两性夫妻黄色片| 亚洲色图综合在线观看| 欧美成人午夜免费资源| 熟女电影av网| 91久久精品国产一区二区三区| 国产精品女同一区二区软件| 最近的中文字幕免费完整| 宅男免费午夜| 国产精品无大码| 久久99热这里只频精品6学生| 欧美 日韩 精品 国产| 久久国产精品大桥未久av| 亚洲精品乱久久久久久| 欧美日韩av久久| 考比视频在线观看| 一区在线观看完整版| 9色porny在线观看| av天堂久久9| 精品国产国语对白av| 国产精品久久久久久av不卡| 不卡av一区二区三区| 老汉色∧v一级毛片| av视频免费观看在线观看| 天天躁夜夜躁狠狠久久av| 男女边吃奶边做爰视频| 91精品伊人久久大香线蕉| 国产成人欧美| 欧美黄色片欧美黄色片| 欧美日韩成人在线一区二区| av天堂久久9| 久久人人97超碰香蕉20202| 亚洲精品美女久久久久99蜜臀 | 久久久久久久大尺度免费视频| 18+在线观看网站| 肉色欧美久久久久久久蜜桃| 美女大奶头黄色视频| 99香蕉大伊视频| 精品一区二区免费观看| 色吧在线观看| 国产成人a∨麻豆精品| 97精品久久久久久久久久精品| 国产野战对白在线观看| 在线天堂中文资源库| 免费高清在线观看视频在线观看| 国产视频首页在线观看| 伊人久久大香线蕉亚洲五| 黄频高清免费视频| 色视频在线一区二区三区| 亚洲av在线观看美女高潮| 纵有疾风起免费观看全集完整版| xxxhd国产人妻xxx| 午夜影院在线不卡| 国产精品一国产av| 色播在线永久视频| 久久久欧美国产精品| 观看av在线不卡| av电影中文网址| 免费看不卡的av| 欧美日韩一级在线毛片| 2021少妇久久久久久久久久久| 国产精品国产av在线观看| freevideosex欧美| 黑人猛操日本美女一级片| 国产精品免费大片| 国产欧美亚洲国产| 美女脱内裤让男人舔精品视频| 男女边摸边吃奶| 少妇被粗大猛烈的视频| 天堂8中文在线网| 亚洲精品国产av成人精品| 激情视频va一区二区三区| 涩涩av久久男人的天堂| 女人高潮潮喷娇喘18禁视频| 国产男女超爽视频在线观看| 日本-黄色视频高清免费观看| 国产福利在线免费观看视频| 纵有疾风起免费观看全集完整版| 亚洲av男天堂| 91久久精品国产一区二区三区| 一本大道久久a久久精品| 亚洲国产成人一精品久久久| 80岁老熟妇乱子伦牲交| 久久国产精品大桥未久av| 亚洲激情五月婷婷啪啪| 爱豆传媒免费全集在线观看| 免费黄色在线免费观看| 成人免费观看视频高清| 国产亚洲午夜精品一区二区久久| 亚洲国产色片| 久久精品熟女亚洲av麻豆精品| 国产成人一区二区在线| 男女下面插进去视频免费观看| 午夜福利网站1000一区二区三区| 精品午夜福利在线看| 又黄又粗又硬又大视频| 香蕉精品网在线| 久热这里只有精品99| 成人亚洲欧美一区二区av| 日韩中字成人| 国产精品免费大片| 观看美女的网站| 亚洲五月色婷婷综合| 亚洲欧美日韩另类电影网站| 亚洲国产精品国产精品| 成人毛片60女人毛片免费| 黄色怎么调成土黄色| av线在线观看网站| 国产黄色免费在线视频| 亚洲成人手机| 熟女av电影| 日日撸夜夜添| 人妻系列 视频| 波野结衣二区三区在线| 欧美另类一区| 亚洲av在线观看美女高潮| 大话2 男鬼变身卡| 男男h啪啪无遮挡| 午夜福利视频在线观看免费| 久久精品亚洲av国产电影网| 一本大道久久a久久精品| 午夜福利,免费看| 18禁裸乳无遮挡动漫免费视频| 国产亚洲一区二区精品| 熟女电影av网| 看免费av毛片| 亚洲人成网站在线观看播放| 国产精品亚洲av一区麻豆 | 狂野欧美激情性bbbbbb| 久久国产精品大桥未久av| 制服人妻中文乱码| 国产精品香港三级国产av潘金莲 | 18在线观看网站| 日韩不卡一区二区三区视频在线| 精品少妇内射三级| 午夜免费观看性视频| 亚洲国产av影院在线观看| 日韩中文字幕欧美一区二区 | 999精品在线视频| 肉色欧美久久久久久久蜜桃| 在线观看美女被高潮喷水网站| 国产精品久久久av美女十八| 免费观看a级毛片全部| 国语对白做爰xxxⅹ性视频网站| 国产精品二区激情视频| 国产精品免费大片| 欧美 日韩 精品 国产| 男女边摸边吃奶| 婷婷色综合www| 成年人免费黄色播放视频| 精品午夜福利在线看| 伊人久久大香线蕉亚洲五| 久久99精品国语久久久| 老司机影院成人| 亚洲精品美女久久久久99蜜臀 | 免费观看在线日韩| av在线播放精品| 赤兔流量卡办理| 亚洲成国产人片在线观看| 婷婷成人精品国产| 电影成人av| 亚洲欧美清纯卡通| 亚洲精品美女久久av网站| 欧美中文综合在线视频| 国产精品国产三级国产专区5o| 日韩熟女老妇一区二区性免费视频| 亚洲欧洲日产国产| 亚洲精品一区蜜桃| 久久精品亚洲av国产电影网| 日韩免费高清中文字幕av| 久久这里有精品视频免费| 日韩一卡2卡3卡4卡2021年| 性色avwww在线观看| 亚洲在久久综合| 在线观看免费日韩欧美大片| 亚洲婷婷狠狠爱综合网| 国产精品久久久av美女十八| 在线观看免费日韩欧美大片| 久久韩国三级中文字幕| 久热这里只有精品99| 成人黄色视频免费在线看| 中文字幕精品免费在线观看视频| 国产精品不卡视频一区二区| 最近最新中文字幕免费大全7| 久久久久久免费高清国产稀缺| 国产精品一国产av| 国产亚洲一区二区精品| 午夜福利网站1000一区二区三区| 91aial.com中文字幕在线观看| 国产精品人妻久久久影院| 女的被弄到高潮叫床怎么办| 天天操日日干夜夜撸| 青春草亚洲视频在线观看| 久久久久久久大尺度免费视频| 男女边摸边吃奶| 亚洲欧美清纯卡通| 精品酒店卫生间| 亚洲国产精品一区三区| 十八禁高潮呻吟视频| 亚洲精品自拍成人| 久久青草综合色| 亚洲国产精品国产精品| 亚洲av国产av综合av卡| 色播在线永久视频| 亚洲精品美女久久久久99蜜臀 | 久久久久国产网址| 性色avwww在线观看| 国产成人精品久久久久久| 9热在线视频观看99| 亚洲精品视频女| 久久精品久久久久久久性| 日韩中字成人| 一个人免费看片子| 色吧在线观看| 国产在线一区二区三区精| 亚洲色图 男人天堂 中文字幕| 欧美国产精品一级二级三级| 亚洲av福利一区| 蜜桃国产av成人99| 人人妻人人澡人人爽人人夜夜| 久久国产亚洲av麻豆专区| 久久久久精品人妻al黑| 日本猛色少妇xxxxx猛交久久| 久久精品国产亚洲av天美| 在线观看免费日韩欧美大片| 曰老女人黄片| 精品视频人人做人人爽| 日韩一卡2卡3卡4卡2021年| 一级爰片在线观看| 精品一区在线观看国产| 亚洲精品美女久久av网站| 国产免费一区二区三区四区乱码| 成人毛片60女人毛片免费| 日韩 亚洲 欧美在线| 国产精品偷伦视频观看了| 久久久久久久国产电影| 日韩三级伦理在线观看| 日韩在线高清观看一区二区三区| av国产精品久久久久影院| 成人亚洲欧美一区二区av| 在线 av 中文字幕| 中文欧美无线码| 黑人巨大精品欧美一区二区蜜桃| av又黄又爽大尺度在线免费看| 两个人免费观看高清视频| 国产一区二区 视频在线| 看十八女毛片水多多多| 在现免费观看毛片| 黄片无遮挡物在线观看| 亚洲少妇的诱惑av| 国产亚洲午夜精品一区二区久久| 亚洲精品美女久久av网站| 人妻人人澡人人爽人人| 亚洲av中文av极速乱| 如何舔出高潮| 精品午夜福利在线看| 18禁国产床啪视频网站| 国产成人a∨麻豆精品| 国产成人精品久久久久久| 国产乱来视频区| 国产国语露脸激情在线看| 少妇精品久久久久久久| 久久热在线av| 叶爱在线成人免费视频播放| 国产精品99久久99久久久不卡 | 欧美精品国产亚洲| 国产精品二区激情视频| 国产在视频线精品| 免费高清在线观看视频在线观看| 久久精品国产亚洲av天美| 又粗又硬又长又爽又黄的视频| 欧美激情 高清一区二区三区| 精品少妇内射三级| a级毛片在线看网站| 18禁动态无遮挡网站| 久久久久视频综合| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 精品一区二区三卡| 国产成人av激情在线播放| 国产成人精品在线电影| 国产精品免费视频内射| 亚洲欧洲国产日韩| tube8黄色片| xxx大片免费视频| 18禁观看日本| 男女无遮挡免费网站观看| 黄片小视频在线播放| 一级毛片电影观看| 热99久久久久精品小说推荐| 午夜久久久在线观看| 日韩精品免费视频一区二区三区| 大片电影免费在线观看免费| 美女午夜性视频免费| 亚洲美女视频黄频| 如何舔出高潮| 国产极品天堂在线| av女优亚洲男人天堂| 成人毛片60女人毛片免费| 欧美 亚洲 国产 日韩一| 亚洲欧美精品自产自拍| 热99国产精品久久久久久7| 永久网站在线| 国产老妇伦熟女老妇高清| 国产成人精品久久久久久| 丰满少妇做爰视频| 最近最新中文字幕大全免费视频 | 一二三四中文在线观看免费高清| 纵有疾风起免费观看全集完整版| 欧美 日韩 精品 国产| 777米奇影视久久| 免费看不卡的av| 亚洲国产欧美网| 亚洲精品视频女| 国产在线一区二区三区精| 丰满乱子伦码专区| 大陆偷拍与自拍| 国产乱来视频区| 一级毛片 在线播放| 18禁国产床啪视频网站| 午夜福利影视在线免费观看| 两性夫妻黄色片| 欧美精品高潮呻吟av久久| 久久久久久久久免费视频了| 国产熟女欧美一区二区| 婷婷色av中文字幕| 一级片'在线观看视频| 九草在线视频观看| 国产精品99久久99久久久不卡 | 久久韩国三级中文字幕| 日韩av在线免费看完整版不卡| 老熟女久久久| 制服人妻中文乱码| 香蕉丝袜av| 亚洲精品在线美女| 欧美精品国产亚洲| 伦理电影大哥的女人| 午夜激情av网站| 日韩精品有码人妻一区| 一区福利在线观看| 精品卡一卡二卡四卡免费| 亚洲av国产av综合av卡| 老女人水多毛片| 日韩成人av中文字幕在线观看| 在线看a的网站| 亚洲国产成人一精品久久久| 又粗又硬又长又爽又黄的视频| 亚洲国产日韩一区二区| a级毛片黄视频| 欧美日韩一级在线毛片| 欧美97在线视频| 美女福利国产在线| 免费高清在线观看视频在线观看| 欧美av亚洲av综合av国产av | 在线亚洲精品国产二区图片欧美| 欧美中文综合在线视频| 大陆偷拍与自拍| 国产无遮挡羞羞视频在线观看| 亚洲综合色网址| 久久久久久久久久久免费av| 国产乱来视频区| 国产精品人妻久久久影院| 欧美激情极品国产一区二区三区| 最近的中文字幕免费完整| 桃花免费在线播放| 亚洲精品美女久久av网站| 九色亚洲精品在线播放| 日本vs欧美在线观看视频| 深夜精品福利| a级毛片黄视频| 精品一区二区免费观看| 国产精品二区激情视频| 欧美日韩精品网址| 国产成人免费观看mmmm| 国产探花极品一区二区| 亚洲欧美一区二区三区国产| av免费在线看不卡| 久久久久国产网址| 亚洲精品成人av观看孕妇| 一本—道久久a久久精品蜜桃钙片| 亚洲精品美女久久av网站| 一级a爱视频在线免费观看| 最近最新中文字幕免费大全7| 日韩,欧美,国产一区二区三区| 亚洲欧美精品综合一区二区三区 | 亚洲美女搞黄在线观看| 精品第一国产精品| 交换朋友夫妻互换小说| 男女午夜视频在线观看| av国产精品久久久久影院| 亚洲av国产av综合av卡| 天天影视国产精品| 高清不卡的av网站| 只有这里有精品99| 日韩精品免费视频一区二区三区| 欧美日韩国产mv在线观看视频| 国产精品一国产av| 国产午夜精品一二区理论片| 亚洲综合精品二区| 久久久久精品久久久久真实原创|