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

    自動(dòng)雙環(huán)入滲儀設(shè)計(jì)與試驗(yàn)

    2020-07-10 05:28:40張建豐胡克林
    關(guān)鍵詞:馬氏雙環(huán)恒定

    孫 權(quán),張建豐,張 杰,曾 輝,胡克林,劉 剛

    ·研究速報(bào)·

    自動(dòng)雙環(huán)入滲儀設(shè)計(jì)與試驗(yàn)

    孫 權(quán)1,2,張建豐3,張 杰1,曾 輝1,胡克林1,劉 剛1※

    (1. 中國農(nóng)業(yè)大學(xué)土地科學(xué)與技術(shù)學(xué)院,北京 100193;2. 南寧市農(nóng)業(yè)農(nóng)村局,南寧 530000;3. 西安理工大學(xué)水利水電學(xué)院,西安 710048)

    為克服傳統(tǒng)雙環(huán)入滲儀的供水不穩(wěn)定、讀取和記錄數(shù)據(jù)工作量大的問題,該文基于定水頭法測量導(dǎo)水率的原理,設(shè)計(jì)了一種高效的自動(dòng)雙環(huán)入滲儀。該入滲儀采用非接觸式電容感應(yīng)液位傳感器來監(jiān)測內(nèi)環(huán)和外環(huán)內(nèi)的水位(精度<1.5 mm);將激光測距傳感器(精度<1 mm)安裝在圓柱形蓄水器的頂部,以連續(xù)自動(dòng)監(jiān)測供水量;高精度液位控制系統(tǒng)與電磁閥組合以保持恒定水頭。通過理論分析,增大蓄水器中的液面與雙環(huán)中的液面之間的壓力差,可使自動(dòng)雙環(huán)入滲儀的初始供水速率遠(yuǎn)高于馬氏瓶的初始供水速率。自動(dòng)雙環(huán)入滲儀的精度和可靠性試驗(yàn)于2019年7月至2019年10月在中國農(nóng)業(yè)大學(xué)綠苑試驗(yàn)場進(jìn)行,并進(jìn)行了該自動(dòng)雙環(huán)入滲儀和傳統(tǒng)雙環(huán)入滲儀之間的比較。試驗(yàn)結(jié)果顯示,在選定土壤的穩(wěn)定入滲階段,用馬氏瓶測量的穩(wěn)定入滲速率為0.068 5 cm/min,該自動(dòng)雙環(huán)入滲儀的穩(wěn)定入滲速率為0.068 7 cm/min,自動(dòng)雙環(huán)入滲儀實(shí)測數(shù)據(jù)與Philip的理論入滲模型計(jì)算結(jié)果吻合(2>0.99)。結(jié)果表明,該自動(dòng)雙環(huán)入滲儀的自動(dòng)化測量結(jié)果可靠,對今后土壤水分入滲過程研究提供了更加高效的手段。

    土壤;水分;雙環(huán)入滲儀;激光傳感器模塊;非接觸式液位傳感器;電磁閥

    0 引 言

    土壤導(dǎo)水率是控制水分入滲和地表徑流及污染物向地下水遷移的土壤基本物理性質(zhì)之一[1],其大小有時(shí)在空間上差別很大[2-3]。確定田間尺度上的土壤導(dǎo)水率,需要大量的試驗(yàn)測量[4]。雙環(huán)入滲儀是田間原位測量土壤入滲情況的常用試驗(yàn)裝置。該裝置使用定水頭法測量土壤入滲,在測量期間要保持恒定水頭且必須記錄入滲水量隨時(shí)間的變化。獲得一次可靠的飽和導(dǎo)水率測量可能需要幾天的時(shí)間[5],特別是同時(shí)監(jiān)測多個(gè)原位測試點(diǎn),測量過程繁瑣且耗時(shí)。為了減小測量過程中的勞動(dòng)強(qiáng)度,研究雙環(huán)入滲儀測量過程的自動(dòng)化非常有價(jià)值。

    Casey等[6]首次使用差動(dòng)傳感器自動(dòng)測量馬氏瓶液位,單個(gè)差動(dòng)傳感器安裝在滲透儀儲(chǔ)液器的底部,另一端通過管道連接到馬氏瓶的頂部空間。使用馬氏瓶測量,存在2個(gè)缺陷:1)因田間地面不平,難以精確調(diào)整馬氏瓶發(fā)泡口高度與環(huán)內(nèi)所需恒定水頭到同一水平面,同時(shí)也難以精確統(tǒng)一內(nèi)環(huán)與外環(huán)之間的液面高度;2)馬氏瓶維持恒定水頭的原理是馬氏瓶發(fā)泡口與環(huán)內(nèi)液位構(gòu)成的連通器(連通器原理),供水速率由環(huán)內(nèi)液面下降與馬氏瓶發(fā)泡口之間的微小液位差決定,往往在土壤入滲的初始階段,馬氏瓶的供水速率會(huì)低于土壤的入滲速率,因此,不能保證環(huán)內(nèi)的恒定水頭。最早的改進(jìn)工作之一由Constantz等[7]完成,他們將壓力傳感器連接到馬氏瓶的頂部,通過氣體壓強(qiáng)的變化來測量流量,但該方法需要在內(nèi)環(huán)和外環(huán)中手動(dòng)控制液位,無法達(dá)到自動(dòng)化。為了提高供水流量的測量精度,Ankeny等[8]對該裝置改進(jìn),通過2個(gè)壓力傳感器之間的壓力差來計(jì)算馬氏瓶中液位隨時(shí)間的變化,一個(gè)壓力傳感器位于馬氏瓶頂部,另一個(gè)位于馬氏瓶底部,但同樣存在上述問題。Matula等[9]開發(fā)了一種用于恒定水頭測試的半自動(dòng)化測量系統(tǒng),該系統(tǒng)需要2個(gè)水容器,分別用于內(nèi)環(huán)和外環(huán)供水,內(nèi)環(huán)中的液位通過浮子和光敏晶體管來控制,但需要手動(dòng)記錄時(shí)間,也未實(shí)現(xiàn)自動(dòng)化。Prieksat等[10]基于前人的工作[7,9],使用壓力傳感器來確定馬氏瓶的供水流量,采用數(shù)據(jù)記錄儀記錄數(shù)據(jù),通過調(diào)節(jié)馬氏瓶發(fā)泡口高度可將恒定水頭精確到±1 mm,但仍存在上述的2個(gè)缺陷。Maheshwari[11]開發(fā)了一種自動(dòng)雙環(huán)入滲儀,該設(shè)計(jì)采用2個(gè)電極接觸式液位傳感器(用于保持內(nèi)環(huán)和外環(huán)恒定水頭)、電容式水位傳感器(用于計(jì)量蓄水罐對雙環(huán)的供水量)、電磁閥和12 V電池(用于為電磁閥供電),電腦控制閥門并從傳感器收集數(shù)據(jù)。但該裝置需要液位變化5 mm以上才會(huì)觸發(fā)該電極接觸式液位傳感器開關(guān),精度較低。Milla等[12]使用紅外液位測量傳感器系統(tǒng)(安裝在外環(huán)或馬氏瓶上)測量液位隨時(shí)間的變化,但未擺脫使用馬氏瓶。Lazarovitch等[13]開發(fā)的半自動(dòng)恒定水頭單環(huán)入滲儀,通過使用懸掛的S型稱重傳感器連續(xù)稱重蓄水池來測量流量,但單環(huán)不能保證水分垂直一維入滲。Arriaga等[14]利用連接到數(shù)據(jù)記錄器的小型壓力傳感器來自動(dòng)采集數(shù)據(jù),但該系統(tǒng)需操作員集中精力維持內(nèi)環(huán)和外環(huán)的液位一致,并根據(jù)需要重新補(bǔ)水。Ong等[15]通過使用壓力傳感器、微控制器、開源電子設(shè)備組合的自動(dòng)化設(shè)備和安裝在馬氏瓶上的壓力傳感器測量測試期間的水量,同樣使用了馬氏瓶。Fatehnia等[16]利用霍爾效用傳感器、蠕動(dòng)泵、液位傳感器和恒定水平浮閥開發(fā)的自動(dòng)雙環(huán)入滲儀,在入滲速率較大的土壤中需要另加蠕動(dòng)泵,同時(shí)采用的液位傳感器電阻易受溫度的影響。除此之外,諸如BEST之類的簡化方法的優(yōu)點(diǎn)是它們的簡單性和廉價(jià)性[17-18],但測量結(jié)果不可靠。

    雙環(huán)自動(dòng)化發(fā)展過程經(jīng)歷了使用接觸式、壓力式、電阻式和電容式液位傳感器的過程。由于水接觸式電容或電阻深度計(jì)對水質(zhì)敏感并易受腐蝕,導(dǎo)致過度維護(hù)和校準(zhǔn),從而導(dǎo)致結(jié)果不準(zhǔn)確。Casey等[6-8,10,12,15]開發(fā)的自動(dòng)入滲儀都使用馬氏瓶作為供水源,在田間使用馬氏瓶測量土壤入滲速率具有一定的局限性。基于上述問題,本研究設(shè)計(jì)的自動(dòng)雙環(huán)入滲儀采用液位控制系統(tǒng)和雙環(huán)一體化組合結(jié)構(gòu),從根本上解決傳統(tǒng)馬氏瓶供水速率低、馬氏瓶與雙環(huán)入滲儀配合使用時(shí)難以精確確定雙環(huán)所需的恒定水頭以及難以統(tǒng)一內(nèi)外環(huán)液位高度等問題,采用激光測距傳感器解決入滲水量的自動(dòng)記錄。該設(shè)計(jì)可為今后土壤水分入滲的過程研究提供更加高效的手段,非接觸式傳感器的使用為自動(dòng)雙環(huán)入滲儀的設(shè)計(jì)提供了新的思路。

    1 自動(dòng)雙環(huán)入滲儀整體結(jié)構(gòu)與工作原理

    1.1 整體結(jié)構(gòu)

    本文根據(jù)定水頭法測量土壤導(dǎo)水率的原理,設(shè)計(jì)了一種自動(dòng)記錄、自動(dòng)保持液位平衡、初始供水速率大的自動(dòng)雙環(huán)入滲儀。該自動(dòng)雙環(huán)入滲儀主要由非接觸式液體感應(yīng)電容液位傳感器、圓柱形水容器、帶串口調(diào)試軟件的激光測距傳感器模塊、電磁閥、12 V電源、筆記本電腦等組成。圓柱形水容器作為供水源;激光測距傳感器與筆記本電腦配套使用,能代替人工讀數(shù)完成數(shù)據(jù)采集;液位傳感器、液位控制器和電磁閥配套使用能保持雙環(huán)液位恒定。整體結(jié)構(gòu)如圖1所示,主要結(jié)構(gòu)參數(shù)如表1所示。

    1.2 工作原理

    傳統(tǒng)雙環(huán)入滲儀采用馬氏瓶保持恒定水頭測量飽和導(dǎo)水率,馬氏瓶利用連通器原理維持恒定水頭,如圖2a。隨著環(huán)內(nèi)水分的入滲,環(huán)內(nèi)液面與馬氏瓶發(fā)泡口產(chǎn)生微小液位差,導(dǎo)致馬氏瓶發(fā)泡口與環(huán)內(nèi)液面形成壓強(qiáng)差。該壓強(qiáng)差提供馬氏瓶的供水動(dòng)力。根據(jù)液體壓強(qiáng)公式有:

    ?1=(12)(1)

    式中?1是馬氏瓶發(fā)泡口與環(huán)內(nèi)液面形成的壓強(qiáng)差,Pa;是水的密度,g/cm3;是重力加速度,m/s2;1是發(fā)泡口離地表參考平面的距離,m;2是傳統(tǒng)雙環(huán)入滲儀的環(huán)內(nèi)液面離地表參考平面的高度,m;2圍繞恒定水頭動(dòng)態(tài)波動(dòng)[10]。

    1.激光測距傳感器 2.泡沫塑料浮 3.電磁閥 4.液位控制器 5.液位傳感器 6.電源 7.圓柱形水容器 8.內(nèi)環(huán) 9.外環(huán)1.Laser range sensor 2.Foam plastic float 3.Solenoid valve 4.Liquid level controller 5.Liquid level sensor 6.Power 7.Cylindrical water container 8.Inner ring 9.Outer ring

    表1 自動(dòng)雙環(huán)入滲儀的主要結(jié)構(gòu)參數(shù)

    本自動(dòng)雙環(huán)入滲儀采用非接觸式感應(yīng)電容液位控制器與電磁閥維持恒定水頭,利用圓柱形水容器蓄水,從而增大供水動(dòng)力,即增大圓柱形水容器的蓄水液面與環(huán)內(nèi)液面形成的壓強(qiáng)差,如圖2b。根據(jù)液體壓強(qiáng)公式有:

    ?2=(3) (2)

    式中?2是圓柱形水容器的蓄水液面與環(huán)內(nèi)液面形成的壓強(qiáng)差,Pa;是圓柱形水容器的液面相對于地表參考平面的高度,m;3是自動(dòng)雙環(huán)入滲儀的環(huán)內(nèi)液面離地表參考平面的高度,m;3圍繞恒定水頭動(dòng)態(tài)波動(dòng)。在利用傳統(tǒng)馬氏瓶與本自動(dòng)雙環(huán)入滲儀維持同一恒定水頭時(shí),因?yàn)?、3均圍繞恒定水頭動(dòng)態(tài)變化,即2=3,又(1-2)很微小,則1-3>>1-2,即?2>>?1。因此,該自動(dòng)雙環(huán)入滲儀的供水速率遠(yuǎn)遠(yuǎn)大于傳統(tǒng)馬氏瓶的供水速率。

    在使用傳統(tǒng)馬氏瓶保持雙環(huán)入滲儀恒定水頭時(shí),需要調(diào)節(jié)馬氏瓶發(fā)泡口與環(huán)內(nèi)液面到同一水平面(圖2a)。由于田間地面不平、馬氏瓶與雙環(huán)入滲儀的位置獨(dú)立,很難調(diào)節(jié)馬氏瓶發(fā)泡口與雙環(huán)入滲儀環(huán)內(nèi)液面到同一水平面,即很難精確調(diào)節(jié)環(huán)內(nèi)液面至所需恒定水頭。本雙環(huán)與非接觸式液位傳感器的一體化設(shè)計(jì),液位控制器結(jié)合電磁閥保持恒定水頭,可上下調(diào)節(jié)液位傳感器至所需恒定水頭(圖2b)。

    1.馬氏瓶 2.馬氏瓶發(fā)泡口 3.環(huán)內(nèi)液面 4.圓柱形水容器 5.泡沫塑料浮 6.激光測距傳感器 7.參考平面 8.內(nèi)環(huán) 9.外環(huán)1.Mariotte reservoir 2.Foam end of Mariotte reservoir 3.Liquid level of inner ring 4.Cylindrical water container 5.Foam plastic float 6.Laser range sensor 7.Reference plane 8.Inner ring 9.Outer ring

    注:為圓柱形水容器的液面相對于地表參考平面的高度,m;0為激光傳感器到圓柱形水容器的液面距離,m;1為發(fā)泡口離地表參考平面的距離,m;2為傳統(tǒng)雙環(huán)入滲儀的環(huán)內(nèi)液面離地表參考平面的高度,m;3為自動(dòng)雙環(huán)入滲儀的環(huán)內(nèi)液面離地表參考平面的高度,m。

    Note:is the height of the liquid surface of the cylindrical water reservoir relative to the ground reference plane, m;0is the distance from the laser sensor to the liquid surface of the cylindrical water reservoir, m;1is the distance of the foaming port from the surface reference plane, m ;2is the height of the liquid level in the ring of the traditional double-ring infiltrometer from the ground reference plane, m;3is the height of the liquid level in the ring of the traditional double-ring infiltrometer from the ground reference plane, m.

    圖2 雙環(huán)入滲儀工作原理示意圖

    Fig.2 Schematic diagram of working principle of double-ring infiltrometer

    2 關(guān)鍵部件設(shè)計(jì)

    2.1 圓柱形雙環(huán)

    2.1.1 雙環(huán)尺寸選擇依據(jù)

    ASTM(American Standard Test Method)雙環(huán)入滲儀裝置包括開放的內(nèi)外同心金屬環(huán)[19]。研究表明,內(nèi)外環(huán)的尺寸對飽和導(dǎo)水率測量有影響,尺寸較小的入滲儀通常會(huì)高估入滲速率,并且提出了幾個(gè)雙環(huán)的參考尺寸。Swartzendruber等[20]在沙土的研究發(fā)現(xiàn),外環(huán)和內(nèi)環(huán)半徑分別為60和50 cm的入滲測量結(jié)果最令人滿意。Ahuja等[21]研究表明外環(huán)直徑為90 cm、內(nèi)環(huán)直徑為30 cm時(shí)能夠消除側(cè)向流的影響,促進(jìn)水分向下垂直一維入滲。Bouwer等[22]建議,應(yīng)使用至少100 cm的直徑才能精確測量土壤的飽和導(dǎo)水率。Youngs[23]得出結(jié)論,當(dāng)內(nèi)環(huán)尺寸至少為15 cm時(shí),不同位置的測量結(jié)果差異性較小。Lai等[24]研究表明,內(nèi)環(huán)直徑大于80 cm才能獲得可靠原位測量的飽和導(dǎo)水率結(jié)果。Fatehnia等[25]還考慮了積水頭、雙環(huán)插入土壤深度、土壤初始含水量和土壤類型對穩(wěn)定入滲速率的影響。Milla[12]的研究表明,使用15~30 cm內(nèi)環(huán)直徑測量的結(jié)果較好且方法相對簡單,可以反映土壤入滲變化等有價(jià)值的信息。因此,為簡化操作,本研究參考Milla[12]的結(jié)論,采用內(nèi)環(huán)和外環(huán)直徑分別為17.6和30 cm(外環(huán)為直徑30 cm的圓桶,內(nèi)環(huán)由直徑17.6 cm的化工罐切割而成)。

    2.1.2 一體式雙環(huán)設(shè)計(jì)

    雙環(huán)入滲儀在入土的過程中均會(huì)破壞土壤結(jié)構(gòu),環(huán)壁和土壤之間會(huì)出現(xiàn)優(yōu)先流路徑[26]。因此,選擇由0.1 cm厚的鍍鋅鐵皮制成的內(nèi)外環(huán),將它們插入地面時(shí)阻力小,對土壤表面不會(huì)產(chǎn)生較大的擾動(dòng)。內(nèi)環(huán)在距離底部2 cm處打一個(gè)直徑為1 cm的圓孔,外環(huán)在距離底部2 cm處打2個(gè)直徑為1 cm的圓孔。內(nèi)環(huán)的圓孔通過孔徑為1 cm的直通管與外環(huán)的一個(gè)圓孔相連,外環(huán)上的2個(gè)孔分別裝上孔徑為1.2 cm的L型直角彎頭,直角彎頭接上外徑1.2 cm、壁厚0.1 cm的有機(jī)玻璃管(本試驗(yàn)采用的液位傳感器在管道外徑1~2 cm范圍內(nèi),能夠感應(yīng)到的最大管壁厚度的范圍是(1.5±0.5) mm[27])。開孔處與內(nèi)外環(huán)連接管相連時(shí),環(huán)的2面必須加裝橡膠墊,擰緊螺帽以防止漏水,如圖3所示。

    1.內(nèi)環(huán) 2.外環(huán) 3.內(nèi)外環(huán)連接管

    2.2 非接觸式感應(yīng)電容液位控制器

    非接觸式液位控制器由XKC-Y26-V液位傳感器探頭(以下簡稱液位傳感器)[28]與C352控制器組合而成,液位傳感器可間隔非金屬容器感應(yīng)液位的變化,從而控制電磁閥的開閉,安裝簡易且不影響液面。該液位傳感器利用水的感應(yīng)電容來檢測液體的存在,當(dāng)液位達(dá)到感應(yīng)點(diǎn)時(shí),液位控制器關(guān)閉電磁閥;當(dāng)液位離開感應(yīng)點(diǎn)時(shí),液位控制器打開電磁閥。

    2個(gè)液位傳感器(內(nèi)環(huán)和外環(huán)各一個(gè))用于保持內(nèi)環(huán)和外環(huán)中適當(dāng)?shù)乃^。液位傳感器裝在外徑1.2 cm、壁厚0.1 cm和長8 cm的有機(jī)玻璃管上,液位傳感器高度可根據(jù)水頭需要上下調(diào)節(jié)(如圖4所示)。將2個(gè)液位傳感器設(shè)定在同一恒定水頭,當(dāng)液位低于此高度時(shí),液位傳感器給控制器發(fā)送信號,打開電磁閥給環(huán)內(nèi)供水;當(dāng)液位高于此高度時(shí),液位傳感器給控制器發(fā)送信號,關(guān)閉電磁閥,停止向環(huán)內(nèi)供水。由于此傳感器靈敏度較高(液位誤差<1.5 mm),可將系統(tǒng)的水頭近似為恒定水頭。

    1.外環(huán) 2.液位傳感器

    2.3 激光測距傳感器

    入滲測量期間液位高度的下降由激光測距傳感器來測定。激光測距傳感器直接與筆記本電腦相連,在筆記本電腦上通過串口調(diào)試軟件實(shí)時(shí)監(jiān)測并記錄圓柱形水容器液位隨時(shí)間的變化。該激光測距傳感器包含一個(gè)激光發(fā)射器和一個(gè)大表面積光敏探測器,支持USB通信和WINDOWS操作系統(tǒng),采用 USB接口供電,測量范圍0.03~80 m,測量精度±1 mm。二級激光,功率小于1 mW;在10 m處的光斑直徑約為6 mm;可在-20~70 ℃的溫度范圍工作。整個(gè)傳感器體積小,質(zhì)量輕(約60 g),其安裝如圖5所示。

    2.4 電磁閥

    電磁閥為正泰常閉DC12 V,線圈功率14 W,接管口徑為9.5 mm。采用12 V鉛酸蓄電池供電。電磁閥根據(jù)水流方向安裝在供水端和出水端。因電磁閥需要水平安放,因此內(nèi)外環(huán)供水電磁閥安裝在同一平整木板上。

    3 田間試驗(yàn)

    3.1 試驗(yàn)設(shè)計(jì)

    野外測試試驗(yàn)在中國農(nóng)業(yè)大學(xué)綠苑試驗(yàn)田進(jìn)行,用100 cm3環(huán)刀在試驗(yàn)點(diǎn)垂直剖面0~30 cm深度每隔5 cm分層取土樣,采用烘干法分別測定土壤容重ρ(g/cm3)和初始體積含水率0(cm3/cm3)。用激光粒度儀(馬爾文激光粒度儀Mastersizer2000)測定土壤顆粒粒徑分布。土壤質(zhì)地為砂壤土,基本物理性質(zhì)見表2。土壤孔隙度(cm3/cm3)由ρ(g/cm3)計(jì)算,假定土壤顆粒密度為2.65 g/cm3。假設(shè)土壤飽和含水量θ(cm3/cm3)與一致[29-31]??紤]到在田間尺度上0和ρ的空間變異性[32],各參數(shù)值為測試結(jié)果的平均值。

    表2 不同深度土壤基本物理特性參數(shù)

    為檢驗(yàn)本自動(dòng)雙環(huán)入滲儀保持恒定水頭的可靠性,試驗(yàn)使用同樣大小的雙環(huán),原位測量土壤穩(wěn)定入滲階段的累積入滲量隨時(shí)間的變化。為充分保證土壤達(dá)到穩(wěn)定入滲階段,利用自動(dòng)雙環(huán)入滲儀保持恒定水頭入滲7 h后,使用馬氏瓶供水2 h并人工記錄數(shù)據(jù),再用自動(dòng)雙環(huán)入滲儀供水,保持同樣恒定水頭,自動(dòng)記錄數(shù)據(jù),保證2種方法原位測量條件相同。

    3.2 試驗(yàn)結(jié)果與分析

    3.2.1 土壤初始入滲階段的測量結(jié)果對比

    測量結(jié)果如圖6所示。由圖6可知,在入滲的初始階段,自動(dòng)雙環(huán)入滲儀的供水速率遠(yuǎn)遠(yuǎn)高于馬氏瓶的供水速率,能測量到土壤入滲初始階段較高的入滲速率。在田間使用馬氏瓶保持恒定水頭時(shí),由于土壤初始入滲速率很高(特別在土壤吸力很大或存在大孔隙時(shí)),馬氏瓶供水速率往往會(huì)低于土壤的初始入滲速率,即會(huì)出現(xiàn)環(huán)內(nèi)液面逐漸下降的過程。環(huán)內(nèi)液面逐漸下降,與馬氏瓶發(fā)泡口的液位差逐漸增大,導(dǎo)致馬氏瓶的供水速率逐漸增大,即測得的土壤入滲速率在初始階段逐漸增大,該階段測得的土壤入滲速率不能表征土壤的實(shí)際入滲速率。隨著入滲的進(jìn)行,入滲速率逐漸變小,當(dāng)?shù)陀隈R氏瓶的供水速率時(shí),環(huán)內(nèi)的液面逐漸上升到馬氏瓶發(fā)泡口附近,此刻開始馬氏瓶的水量變化才表征土壤的入滲情況。

    圖6 初始入滲階段馬氏瓶與自動(dòng)雙環(huán)入滲儀測量結(jié)果對比

    3.2.2 穩(wěn)定入滲階段的測量結(jié)果比較

    在土壤穩(wěn)定入滲狀態(tài)階段,土壤入滲速率相對初始入滲速率小很多,因此馬氏瓶的供水速率能滿足土壤入滲的需要。該入滲階段土壤累積入滲量隨時(shí)間的變化關(guān)系滿足一元一次方程,即

    =(3)

    式中是累積入滲量,cm;是入滲歷時(shí),min;是土壤穩(wěn)定入滲速率,cm/min。

    在穩(wěn)定入滲階段,使用常規(guī)馬氏瓶維持恒定水頭(人工記錄數(shù)據(jù))測量了125 min,然后再用自動(dòng)雙環(huán)入滲儀維持恒定水頭(自動(dòng)采集并記錄數(shù)據(jù))測量了125 min。結(jié)果如圖7所示。由圖7可知,馬氏瓶測量結(jié)果的曲線斜率為0.068 5,即穩(wěn)定入滲速率為0.068 5 cm/min,自動(dòng)雙環(huán)入滲儀的曲線斜率為0.068 7,即穩(wěn)定入滲速率為0.068 7 cm/min,二者2均在0.99以上。自動(dòng)雙環(huán)入滲儀在實(shí)現(xiàn)自動(dòng)化測量的同時(shí),與馬氏瓶在土壤穩(wěn)定入滲階段的測量結(jié)果一致,證明該自動(dòng)雙環(huán)入滲儀的精度可靠。

    圖7 穩(wěn)定入滲階段馬氏瓶與自動(dòng)雙環(huán)入滲儀測量結(jié)果對比

    3.2.3 Philip入滲模型校驗(yàn)

    Philip入滲模型在田間入滲試驗(yàn)中得到了很好的驗(yàn)證[33]。因此,本文引入Philip入滲模型校驗(yàn)該自動(dòng)雙環(huán)入滲儀的測量數(shù)據(jù)質(zhì)量,Philip入滲模型可表示為[34]

    ()0.5(4)

    式中為土壤吸滲率,cm/min0.5。圖8為雙環(huán)入滲儀田間測試的結(jié)果。Philip入滲模型擬合結(jié)果表明:該自動(dòng)雙環(huán)入滲儀測得的土壤累積入滲量隨時(shí)間的變化能很好地滿足Philip入滲模型,20.99。因此,該自動(dòng)雙環(huán)入滲儀測得的數(shù)據(jù)是可靠的,對傳統(tǒng)雙環(huán)入滲儀有改進(jìn)價(jià)值。

    圖8 累積入滲量

    4 討 論

    馬氏瓶保持恒定水頭是相對的。使用馬氏瓶供水時(shí),環(huán)內(nèi)液面隨著入滲的進(jìn)行與發(fā)泡口有一定液位差時(shí),才會(huì)激發(fā)馬氏瓶供水,供水達(dá)到平衡時(shí),馬氏瓶停止供水一段時(shí)間,當(dāng)環(huán)內(nèi)液面與發(fā)泡口再存在一定液位差時(shí),又激發(fā)馬氏瓶供水,即馬氏瓶保持恒定水頭是相對的。Prieksat等[10]在使用馬氏瓶時(shí)發(fā)現(xiàn)調(diào)節(jié)馬氏瓶發(fā)泡口高度可將恒定水頭精確到±1 mm。Maheshwari[11]采用電極傳感器控制電磁閥是否供水來保持恒定水頭,該傳感器上下電極之間的距離是5 mm。本文自動(dòng)雙環(huán)入滲儀的液位傳感器利用水的感應(yīng)電容來檢測液體的存在,反應(yīng)靈敏,液位誤差在1.5 mm以內(nèi),可近似一個(gè)恒定水頭,相對Maheshwari采用的電極傳感器誤差更小。

    在田間使用馬氏瓶保持恒定水頭時(shí),由于土壤初始入滲速率很高,馬氏瓶供水速率往往會(huì)低于土壤的初始入滲速率。該自動(dòng)雙環(huán)入滲儀利用液位控制器結(jié)合電磁閥保持液位恒定,充分利用水的重力勢能,提高了圓柱形水容器的蓄水液位與環(huán)內(nèi)液位之間的液位差。在初始入滲階段,圓柱形水容器的液位高,可提供足夠的供水速率來保持環(huán)內(nèi)的恒定水頭,滿足測量的需要。

    該自動(dòng)雙環(huán)入滲儀采用液位傳感器以及測量水量變化的激光測距傳感器,不必與液體直接接觸,不會(huì)受到腐蝕性液體的腐蝕,不受水垢等雜物影響,不受磁場、壓力等特殊環(huán)境影響,不需要接觸液體,不污染水源,穩(wěn)定性高,靈敏度高(液位傳感器誤差<1.5 mm,激光測距傳感器誤差<1 mm),抗干擾能力強(qiáng),使用壽命長,因此無需反復(fù)校準(zhǔn)傳感器。液位傳感器與雙環(huán)的一體化設(shè)計(jì),可方便快捷地調(diào)節(jié)至所需恒定水頭,避免了因田間地面不平恒定水頭難以精確調(diào)節(jié)的問題。通過原位測量土壤穩(wěn)定入滲階段的累積入滲量表明該自動(dòng)雙環(huán)儀與馬氏瓶的測量效果一致,測量精度高。

    5 結(jié) 論

    1)本研究利用高精度液位傳感器控制系統(tǒng)與電磁閥結(jié)合來保持水頭恒定,利用激光測距傳感器測量水量變化,該自動(dòng)雙環(huán)入滲儀實(shí)現(xiàn)了從土壤入滲初始階段到穩(wěn)定階段的連續(xù)自動(dòng)測量,率先將非接觸式感應(yīng)電容液位傳感器應(yīng)用于自動(dòng)雙環(huán)入滲儀。

    2)內(nèi)外環(huán)與非接觸式液體感應(yīng)電容液位傳感器(精度<1.5 mm)的一體式設(shè)計(jì)解決了難以精確確定恒定水頭的問題,解決了在入滲初始階段馬氏瓶供水速率低于土壤入滲速率的問題。

    3)該自動(dòng)雙環(huán)入滲儀測量的穩(wěn)定入滲速率與馬氏瓶的測量結(jié)果一致;測量的數(shù)據(jù)質(zhì)量能夠很好地滿足Philip入滲模型,R均大于0.99。

    在今后開發(fā)自動(dòng)雙環(huán)入滲儀時(shí),可考慮使用其他非接觸式的傳感器。圍繞供水系統(tǒng)、恒定水頭控制系統(tǒng)和自動(dòng)記錄系統(tǒng)3個(gè)方面開展研究,供水系統(tǒng)主要解決供水量與供水速率2方面問題,恒定水頭控制系統(tǒng)和自動(dòng)記錄系統(tǒng)可在提高精度方面做進(jìn)一步優(yōu)化,同時(shí)可探索在極低水頭時(shí)飽和導(dǎo)水率的測定。

    [1]Reynolds W D, Bowman B T, Brunke R R, et al. Comparison of tension infiltrometer, pressure infiltrometer, and soil core estimates of saturated hydraulic conductivity[J]. Soil Science Society of America Journal, 2000, 64(2): 478-484.

    [2]Ersahin S. Comparing ordinary kriging and cokriging to estimate infiltration rate[J]. Soil Science Society of America Journal, 2003, 67(6): 1848-1855.

    [3]Cichota R, van Lier Q de Jong, Leguizamón Rojas C A. Spatial variability of infiltration in an Alfisol[J]. Revista Brasileira de Ciencia do Solo, 2003, 27: 789-798.

    [4]Logsdon S D, Jaynes D B. Spatial variability of hydraulic conductivity in a cultivated field at different times[J]. Soil Science Society of America Journal, 1996, 60(3): 703-709.

    [5]Diamond J, Shanley T. Infiltration rate assessment of some major soils[J]. Irish Geography, 2003, 36(1): 32-46.

    [6]Casey F X M, Derby N E. Improved design for an automated tension infiltrometer[J]. Soil Science Society of America Journal, 2002, 66(1): 64-67.

    [7]Constantz J, Murphy F. An automated technique for flow measurements from Mariotte reservoirs[J]. Soil Science Society of America Journal, 1987, 51(1): 252-254.

    [8]Ankeny M D, Kaspar T C, Horton R. Design for an automated tension infiltrometer[J]. Soil Science Society of America Journal, 1988, 52: 893-896.

    [9]Matula S, Dirksen C. Automated regulating and recording system for cylinder infiltrometer[J]. Soil Science Society of America Journal, 1989, 53(1): 299-302.

    [10]Prieksat M A, Ankeny M D, Kaspar T C. Design for an automated, self-regulating, single-ring infiltrometer[J]. Soil Science Society of America Journal, 1992, 56(5): 1409-1411.

    [11]Maheshwari B. Development of an automated double-ring infiltrometer[J]. Australian Journal of Soil Research, 1996, 34(5): 709-714.

    [12]Milla K, Kish S. A low-cost microprocessor and infrared sensor system for automating water infiltration measurements[J]. Comput. Electron. Agric., 2006, 53(2): 122-129.

    [13]Lazarovitch N, Ben-Gal A, ?im?nek J, et al. Uniqueness of soil hydraulic parameters determined by a combined wooding inverse approach[J]. Soil Science Society of America Journal, 2007, 71(3): 860-865.

    [14]Arriaga F J, Kornecki T S, Balkcom K S, et al. A method for automating data collection from a double-ring infiltrometer under falling head conditions[J]. Soil Use and Management, 2010, 26(1): 61-67.

    [15]Ong J, Werkema D, Lane J. Revisiting the fully automated double-ring infiltrometer using open-source electronics[A]. American Geophysical Union Fall Meeting, San Francisco, 2012.

    [16]Fatehnia M, Paran S, Kish S, et al. Automating double ring infiltrometer with an Arduino microcontroller[J]. Geoderma, 2016, 262: 133-139.

    [17]Madsen, Matthew D, Chandler, et al. Automation and use of mini disk infiltrometers[J]. Soil Science Society of America Journal, 2007, 71(5): 1469-1472.

    [18]Dohnal M, Dusek J, Vogel T. Improving hydraulic conductivity estimates from minidisk infiltrometer measurements for soils with wide pore-size distributions[J]. Soil Science Society of America Journal, 2010, 74(3): 804-811.

    [19]ASTM, Standard Test Method for Infiltration Rate of Soils in Field Using Double-Ring Infiltrometer: D3385-09(2009) [S]. America: ASTM Standard Test Method, 2009.

    [20]Swartzendruber, Dale, Olsen, et al. Model study of the double ring infiltrometers as affected by depth of wetting and particle size[J]. Angewandte Chemie International Edition, 1961, 54(49/50): 505-507.

    [21]Ahuja L R, El-Swaify S A, Rahman A. Measuring hydrologic properties of soil with a double-ring infiltrometer and multiple-depth tensiometers[J]. Soil Science Society of America Journal, 1976, 40(4): 494-499.

    [22]Bouwer H. Intake rate: Cylinder infiltrometer[J]. Methods of Soil Amnalysis Part Physical & Mineralogical Methods, 1986, 1(1): 825-844.

    [23]Youngs E G. Estimating hydraulic conductivity values from ring infiltrometer measurements[J]. European Journal of Soil Science, 1987, 38(4): 623-632.

    [24]Lai Jianbin, Luo Yi, Ren Li. Buffer index effects on hydraulic conductivity measurements using numerical simulations of double-ring infiltration[J]. Soil Science Society of America Journal, 2010, 74(5): 1526-1536.

    [25]Fatehnia M, Tawfiq K, Abichou T. Comparison of the methods of hydraulic conductivity estimation from mini disk infiltrometer[J]. Electronic Journal of Geotechnical Engineering, 2014, 19: 1047-1063.

    [26]張婧. 土壤入滲與優(yōu)先流測量方法研究[D]. 北京: 中國農(nóng)業(yè)大學(xué),2017. Zhang Jing. Study on Measurement Methodology of Soil Infiltration and Preferential Flow[D]. Beijing: China Agricultural University, 2017. (in Chinese with English abstract)

    [27]深圳市星科創(chuàng)科技有限公司. 一種非接觸式液位傳感器:2016204270127[P]. 2016-11-30.

    [28]深圳市星科創(chuàng)科技有限公司. 一種專用于小管道的液位傳感器:2017213499271[P]. 2018-09-04.

    [29]Mubarak I, Mailhol J C, Angulo-Jaramillo R, et al. Temporal variability in soil hydraulic properties under drip irrigation[J]. Geoderma, 2009, 150(1/2): 158-165.

    [30]Bagarello V, Prima S D, Iovino M, et al. Testing different approaches to characterize Burundian soils by the BEST procedure[J]. Geoderma, 2011, 162(1/2): 141-150.

    [31]Bagarello V, Castellini M, Di Prima S, et al. Soil hydraulic properties determined by infiltration experiments and different heights of water pouring[J]. Geoderma, 2014, 213: 492-501.

    [32]Warrick A W. Environmental Soil Physics[M]. San Diego: Academic Press, 1998.

    [33]范嚴(yán)偉,趙文舉,王昱. 入滲水頭對垂直一維入滲 Philip 模型參數(shù)的影響[J]. 蘭州理工大學(xué)學(xué)報(bào),2015,41(1):65-70. Fan Yanwei, Zhao Wenju, Wang Yu. Influence of infiltration heads on parameters of Philip model of vertical one-dimensional infiltration[J]. Journal of Lanzhou University of Technology, 2015, 41(1): 65-70. (in Chinese with English abstract)

    [34]雷志棟,楊詩秀,謝森傳. 土壤水動(dòng)力學(xué)[M]. 北京:清華大學(xué)出版社,1988.

    Design and experiment of automatic double-ring infiltrometer

    Sun Quan1,2, Zhang Jianfeng3, Zhang Jie1, Zeng Hui1, Hu Kelin1, Liu Gang1※

    (1.,,100193,; 2.,530000,; 3.,,710048,)

    Double-ring infiltrometer is widely used to combine with the Mariotte reservoir for measuring soil hydraulic conductivity. However, there are various kinds of measurement error and shortcomings using the traditional double-ring infiltrometer, particularly on unstable water supply, high labor intense when reading and recording of data. In this study, a new automatic double-ring infiltrometer was designed to improve the measurement precision, while reduce various sorts of errors. A fixed water level method was develop a low-cost and efficient automatic double-ring infiltrometer according to the measuring theory of hydraulic conductivity. The new-developed double-ring infiltrometer was mainly composed of inner (diameter of 17.6 cm) and outer rings(diameter of 30 cm), non-contact liquid induction capacitor level sensor (accuracy <1.5 mm), cylindrical water container, laser range sensor module (accuracy <1 mm), solenoid valve and other components. A non-contact capacitor induction liquid level sensor was selected for monitoring the water level, in order to keep the water level constant inside the both inner and outer ring. A cylindrical water reservoir was served as water source for the container of infiltrometer, and a sensor module with laser range was installed on the top cover of the reservoir. A laptop was connected the laser range sensor module to automatically monitor the volume of water supplied continuously. The combination of the high-precision liquid level control system with the non-contact inductive capacitive liquid level sensor and solenoid valves can be used to maintain a fixed water level in the water reservoir. Therefore, it is possible to increase the pressure difference of liquid level between the water reservoir and the double-ring, and thereby to make the initial rate of water supply in the double-ring infiltrometer much higher than that in the Mariotte reservoir. The new infiltrometer can match the initial rate of water supply, thereby to improve the measuring accuracy of hydraulic conductivity. The accuracy and reliability of the automatic double-ring infiltrometer were verified from July 2019 to October 2019, in the Lüyuan experiment field of China Agricultural University. The experiment was also performed on comparison between the new infiltrometer and the traditional infiltrometer with fixed water level method (double-ring combined with the Mariotte reservoir, collecting data manually). The field experiment results demonstrated that the stable hydraulic conductivity measured by the Mariotte reservoir was 0.068 5 cm/min, and the stable hydraulic conductivity measured by our new automatic double-ring infiltrometer was 0.068 7 cm/min. The data measured by new designed double-ring infiltrometer agreed with the theoretical infiltration model of Philip, where the corresponding2was greater than 0.99. In addition, the new double-ring infiltrometer can continuously and automatically measure the infiltration process, which significantly reduced error caused by traditional double-ring infiltrometer. The new infiltration measurement system can be expected to extend for automatic data sampling with high precision, and low cost, efficient, and easy to use in the field.

    soil; moisture; double-ring infiltrometer; laser range sensor module; non-contact liquid level sensor; solenoid valve

    孫權(quán),張建豐,張杰,等. 自動(dòng)雙環(huán)入滲儀設(shè)計(jì)與試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2020,36(10):318-324.doi:10.11975/j.issn.1002-6819.2020.10.038 http://www.tcsae.org

    Sun Quan, Zhang Jianfeng, Zhang Jie, et al. Design and experiment of automatic double-ring infiltrometer[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(10): 318-324. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2020.10.038 http://www.tcsae.org

    2019-07-03

    2020-04-20

    國家自然科學(xué)基金項(xiàng)目(41771257)

    孫權(quán),主要研究方向?yàn)橥寥廊霛B與優(yōu)先流測量方法。Email:1724573073@qq.com

    劉剛,博士,博士生導(dǎo)師,主要從事多孔介質(zhì)中的能量與質(zhì)量傳遞輸運(yùn)過程、溶質(zhì)和水分的運(yùn)動(dòng)規(guī)律、熱脈沖探針方法的改進(jìn)研究。Email:liug@cau.edu.cn

    10.11975/j.issn.1002-6819.2020.10.038

    S152.7+2;S126

    A

    1002-6819(2020)-10-0318-07

    猜你喜歡
    馬氏雙環(huán)恒定
    一類時(shí)間變換的強(qiáng)馬氏過程
    有環(huán)的可逆馬氏鏈的統(tǒng)計(jì)確認(rèn)
    關(guān)于樹指標(biāo)非齊次馬氏鏈的廣義熵遍歷定理
    花花世界
    一致可數(shù)可加馬氏鏈不變測度的存在性
    漫畫十萬個(gè)為什么
    “單環(huán)學(xué)習(xí)”與“雙環(huán)學(xué)習(xí)”
    電流雙環(huán)控制的LCL單相并網(wǎng)逆變器逆變研究
    聚丙烯成核劑雙環(huán)[2.2.1]-庚烷-2,3-二羧酸鈉的合成
    雙環(huán)法結(jié)合雙“V”形乳腺切除法在乳房肥大整形術(shù)中的應(yīng)用
    国产三级在线视频| 久9热在线精品视频| 国产精品久久久久久久久免 | 亚洲精品456在线播放app | 欧美中文综合在线视频| 99久久九九国产精品国产免费| 国产三级中文精品| 天天添夜夜摸| 日本撒尿小便嘘嘘汇集6| 亚洲在线自拍视频| 麻豆久久精品国产亚洲av| 丰满的人妻完整版| 国产伦人伦偷精品视频| 成人鲁丝片一二三区免费| av天堂在线播放| 九色成人免费人妻av| 禁无遮挡网站| 一本精品99久久精品77| 婷婷精品国产亚洲av| 99国产极品粉嫩在线观看| 国产精品一区二区三区四区久久| 国产欧美日韩一区二区精品| 色综合欧美亚洲国产小说| 18+在线观看网站| 免费在线观看影片大全网站| 欧美黄色片欧美黄色片| 日韩高清综合在线| 日本三级黄在线观看| 少妇熟女aⅴ在线视频| 午夜影院日韩av| 国产三级黄色录像| 老司机午夜十八禁免费视频| 国产精品亚洲美女久久久| 此物有八面人人有两片| 国内毛片毛片毛片毛片毛片| 99久久综合精品五月天人人| 亚洲av五月六月丁香网| 国产精品香港三级国产av潘金莲| 丝袜美腿在线中文| 欧洲精品卡2卡3卡4卡5卡区| 99视频精品全部免费 在线| 国产一区在线观看成人免费| 欧美成人a在线观看| 欧美成人免费av一区二区三区| 国产伦在线观看视频一区| 欧美又色又爽又黄视频| 18禁在线播放成人免费| 精品日产1卡2卡| 久久精品国产亚洲av涩爱 | 国内毛片毛片毛片毛片毛片| 国产亚洲欧美98| 国产精品电影一区二区三区| 国产精品影院久久| 18禁美女被吸乳视频| 成人鲁丝片一二三区免费| 国产伦在线观看视频一区| 色在线成人网| 亚洲国产精品成人综合色| 人妻久久中文字幕网| 久久久久久久久中文| 亚洲 国产 在线| 噜噜噜噜噜久久久久久91| 热99re8久久精品国产| 精品国内亚洲2022精品成人| 久久久久久人人人人人| 亚洲av电影在线进入| 精品熟女少妇八av免费久了| 韩国av一区二区三区四区| 人人妻,人人澡人人爽秒播| 国产一区二区激情短视频| 成人性生交大片免费视频hd| 少妇的逼水好多| 精品电影一区二区在线| 偷拍熟女少妇极品色| 别揉我奶头~嗯~啊~动态视频| 日日摸夜夜添夜夜添小说| 人妻夜夜爽99麻豆av| 免费高清视频大片| 亚洲欧美激情综合另类| 亚洲无线观看免费| 日本 av在线| 久久亚洲真实| 天堂网av新在线| 三级国产精品欧美在线观看| 亚洲精品国产精品久久久不卡| 热99re8久久精品国产| 韩国av一区二区三区四区| 精品久久久久久,| 免费观看人在逋| www国产在线视频色| 老司机深夜福利视频在线观看| 久久久久亚洲av毛片大全| 亚洲国产高清在线一区二区三| 亚洲精品影视一区二区三区av| 日韩高清综合在线| 国产亚洲精品久久久com| 免费观看的影片在线观看| 亚洲五月天丁香| 青草久久国产| 久久草成人影院| 最近最新中文字幕大全电影3| 色综合婷婷激情| 久久香蕉精品热| 国产野战对白在线观看| 男女做爰动态图高潮gif福利片| 国产国拍精品亚洲av在线观看 | 欧美日韩一级在线毛片| 亚洲男人的天堂狠狠| 久久久久免费精品人妻一区二区| 亚洲成人中文字幕在线播放| 最近最新免费中文字幕在线| 午夜久久久久精精品| 精品国产亚洲在线| 欧美性猛交╳xxx乱大交人| 性欧美人与动物交配| 老司机在亚洲福利影院| 美女黄网站色视频| 99热6这里只有精品| 好看av亚洲va欧美ⅴa在| 97人妻精品一区二区三区麻豆| а√天堂www在线а√下载| 亚洲第一欧美日韩一区二区三区| 日本熟妇午夜| 国产主播在线观看一区二区| 最近在线观看免费完整版| 国产av麻豆久久久久久久| 久久久久久九九精品二区国产| 国产三级黄色录像| 久久久色成人| 亚洲最大成人中文| 久久中文看片网| 欧美色视频一区免费| 中文字幕久久专区| 国产午夜精品久久久久久一区二区三区 | 久久这里只有精品中国| 精品国产亚洲在线| 成人国产综合亚洲| 国产精品国产高清国产av| 夜夜看夜夜爽夜夜摸| 黄片小视频在线播放| 午夜免费激情av| 蜜桃久久精品国产亚洲av| 国产亚洲精品一区二区www| 欧美一区二区精品小视频在线| 少妇的逼好多水| 欧美色视频一区免费| 精品国内亚洲2022精品成人| 日韩免费av在线播放| 国产主播在线观看一区二区| 91久久精品国产一区二区成人 | 国产又黄又爽又无遮挡在线| 国产毛片a区久久久久| 国产野战对白在线观看| 最近最新中文字幕大全免费视频| 蜜桃亚洲精品一区二区三区| 在线国产一区二区在线| 日本成人三级电影网站| 97碰自拍视频| 久久中文看片网| 99久久无色码亚洲精品果冻| av片东京热男人的天堂| 国产在视频线在精品| 国产精品国产高清国产av| 亚洲国产精品sss在线观看| 美女高潮喷水抽搐中文字幕| 99久久精品一区二区三区| 不卡一级毛片| 国产真人三级小视频在线观看| 国产伦人伦偷精品视频| 最新在线观看一区二区三区| 啦啦啦观看免费观看视频高清| 动漫黄色视频在线观看| 亚洲,欧美精品.| 欧美日韩一级在线毛片| 国产伦一二天堂av在线观看| 免费看a级黄色片| 色播亚洲综合网| 欧美大码av| 国产精品98久久久久久宅男小说| 国产精品久久久久久亚洲av鲁大| 国产精品免费一区二区三区在线| 露出奶头的视频| 男人和女人高潮做爰伦理| 国产久久久一区二区三区| 亚洲国产欧美网| 亚洲专区国产一区二区| 亚洲成人中文字幕在线播放| 精品久久久久久久毛片微露脸| 午夜福利视频1000在线观看| 成人亚洲精品av一区二区| 成人永久免费在线观看视频| 久久久久久人人人人人| 中文字幕人妻丝袜一区二区| 99久久无色码亚洲精品果冻| 婷婷丁香在线五月| 免费人成在线观看视频色| 12—13女人毛片做爰片一| 99精品欧美一区二区三区四区| 国内毛片毛片毛片毛片毛片| x7x7x7水蜜桃| a级一级毛片免费在线观看| 亚洲av成人不卡在线观看播放网| 女人十人毛片免费观看3o分钟| aaaaa片日本免费| 欧美精品啪啪一区二区三区| 90打野战视频偷拍视频| 级片在线观看| 一级a爱片免费观看的视频| 黄片大片在线免费观看| 久久精品国产综合久久久| 免费看a级黄色片| 久久精品综合一区二区三区| 99在线人妻在线中文字幕| 日本免费一区二区三区高清不卡| 国产一区二区在线观看日韩 | 亚洲精品影视一区二区三区av| 国内久久婷婷六月综合欲色啪| 免费无遮挡裸体视频| 成年女人看的毛片在线观看| 麻豆一二三区av精品| 亚洲最大成人中文| 午夜福利在线观看免费完整高清在 | 日本a在线网址| 国产69精品久久久久777片| 中文字幕高清在线视频| 久久久久久人人人人人| 老汉色av国产亚洲站长工具| 1024手机看黄色片| 国产 一区 欧美 日韩| 国产黄片美女视频| 午夜a级毛片| 亚洲国产精品成人综合色| 国产成人啪精品午夜网站| 无限看片的www在线观看| 99国产极品粉嫩在线观看| 听说在线观看完整版免费高清| 欧美3d第一页| 人人妻人人看人人澡| 99久久精品热视频| 亚洲av不卡在线观看| 好男人电影高清在线观看| 天天躁日日操中文字幕| 在线观看舔阴道视频| 久久中文看片网| 免费看日本二区| 国产高清视频在线播放一区| 黄片大片在线免费观看| 国语自产精品视频在线第100页| 亚洲av免费在线观看| 又爽又黄无遮挡网站| 三级国产精品欧美在线观看| 久久久久久久精品吃奶| 国产精品99久久久久久久久| 国产一区二区在线av高清观看| 欧美极品一区二区三区四区| 国产精品久久久久久亚洲av鲁大| 男插女下体视频免费在线播放| 国产伦一二天堂av在线观看| 噜噜噜噜噜久久久久久91| 黑人欧美特级aaaaaa片| 国内揄拍国产精品人妻在线| 久久精品国产99精品国产亚洲性色| 国产极品精品免费视频能看的| 国内精品美女久久久久久| 一个人免费在线观看电影| 人妻丰满熟妇av一区二区三区| 少妇丰满av| 午夜激情福利司机影院| 悠悠久久av| 男女视频在线观看网站免费| 免费看光身美女| 丰满人妻一区二区三区视频av | 91久久精品电影网| 真人一进一出gif抽搐免费| 又黄又爽又免费观看的视频| 51午夜福利影视在线观看| 美女被艹到高潮喷水动态| 黄色女人牲交| 99热这里只有精品一区| 欧美bdsm另类| 色av中文字幕| 国产真人三级小视频在线观看| 非洲黑人性xxxx精品又粗又长| 国产99白浆流出| 嫩草影院入口| 又黄又爽又免费观看的视频| 69人妻影院| 亚洲午夜理论影院| 琪琪午夜伦伦电影理论片6080| 国产 一区 欧美 日韩| 亚洲片人在线观看| 亚洲欧美精品综合久久99| 在线观看av片永久免费下载| 日韩有码中文字幕| 国内精品美女久久久久久| 最近最新免费中文字幕在线| 黄色日韩在线| 99国产精品一区二区蜜桃av| 亚洲国产日韩欧美精品在线观看 | 香蕉久久夜色| 色综合站精品国产| 国产探花在线观看一区二区| 国产av不卡久久| 九色国产91popny在线| 国产一区在线观看成人免费| 老司机福利观看| 婷婷精品国产亚洲av在线| 日韩亚洲欧美综合| 动漫黄色视频在线观看| 给我免费播放毛片高清在线观看| 午夜两性在线视频| 亚洲精品久久国产高清桃花| 亚洲av二区三区四区| 欧美日韩乱码在线| 午夜免费男女啪啪视频观看 | 757午夜福利合集在线观看| 欧美性猛交╳xxx乱大交人| 麻豆一二三区av精品| 十八禁网站免费在线| 最近最新中文字幕大全电影3| 国产精品av视频在线免费观看| 亚洲一区高清亚洲精品| 亚洲国产精品999在线| bbb黄色大片| 精品一区二区三区视频在线观看免费| 老鸭窝网址在线观看| 在线十欧美十亚洲十日本专区| 天天添夜夜摸| 一本综合久久免费| 性色av乱码一区二区三区2| 午夜免费男女啪啪视频观看 | 免费无遮挡裸体视频| а√天堂www在线а√下载| 欧美日韩国产亚洲二区| 俄罗斯特黄特色一大片| 亚洲成av人片在线播放无| 亚洲熟妇熟女久久| 麻豆一二三区av精品| 国产黄a三级三级三级人| 亚洲欧美日韩高清在线视频| 亚洲国产欧洲综合997久久,| 国产精品电影一区二区三区| 丰满乱子伦码专区| 在线观看舔阴道视频| 精品久久久久久久毛片微露脸| 真人做人爱边吃奶动态| 久久久久久国产a免费观看| 窝窝影院91人妻| 日本精品一区二区三区蜜桃| 国产成人av激情在线播放| 午夜免费观看网址| 成人18禁在线播放| 精品乱码久久久久久99久播| 人妻丰满熟妇av一区二区三区| 欧美黑人欧美精品刺激| 亚洲欧美日韩高清在线视频| 亚洲精品国产精品久久久不卡| 99精品在免费线老司机午夜| 国产国拍精品亚洲av在线观看 | 午夜免费观看网址| 99热这里只有是精品50| 九色成人免费人妻av| 免费在线观看日本一区| 亚洲av日韩精品久久久久久密| 亚洲在线观看片| 亚洲美女视频黄频| aaaaa片日本免费| 制服人妻中文乱码| 天天一区二区日本电影三级| 成人三级黄色视频| 18禁裸乳无遮挡免费网站照片| 一级毛片高清免费大全| 国产主播在线观看一区二区| 男插女下体视频免费在线播放| 熟妇人妻久久中文字幕3abv| 久久精品夜夜夜夜夜久久蜜豆| 久久久色成人| 一个人看视频在线观看www免费 | 久久久久性生活片| 亚洲欧美日韩东京热| 天堂√8在线中文| 99久国产av精品| 我的老师免费观看完整版| 在线观看美女被高潮喷水网站 | 中文字幕精品亚洲无线码一区| 在线a可以看的网站| 少妇的逼好多水| 法律面前人人平等表现在哪些方面| 大型黄色视频在线免费观看| 亚洲人成伊人成综合网2020| 丁香欧美五月| 亚洲久久久久久中文字幕| 午夜a级毛片| 亚洲第一电影网av| 国产精品一区二区三区四区久久| 免费人成在线观看视频色| 国产一区二区亚洲精品在线观看| 91九色精品人成在线观看| 少妇裸体淫交视频免费看高清| 国产探花极品一区二区| 亚洲五月天丁香| 熟女电影av网| 毛片女人毛片| 亚洲精品亚洲一区二区| 国产亚洲欧美98| 18禁黄网站禁片免费观看直播| 亚洲欧美精品综合久久99| 日本 av在线| 一本精品99久久精品77| 少妇的逼好多水| 午夜日韩欧美国产| 日韩欧美免费精品| 香蕉久久夜色| 欧美成人a在线观看| 熟女电影av网| 18禁黄网站禁片午夜丰满| 亚洲 国产 在线| 最新在线观看一区二区三区| 亚洲片人在线观看| 国产成人av激情在线播放| 亚洲精品粉嫩美女一区| 在线播放无遮挡| 亚洲av一区综合| www日本在线高清视频| 欧美成人a在线观看| 国产亚洲精品久久久久久毛片| 午夜精品在线福利| 伊人久久大香线蕉亚洲五| 很黄的视频免费| 色在线成人网| 亚洲第一电影网av| 床上黄色一级片| 午夜免费男女啪啪视频观看 | 亚洲精品在线美女| 熟女人妻精品中文字幕| 国产野战对白在线观看| 亚洲不卡免费看| 88av欧美| 高清日韩中文字幕在线| 亚洲在线自拍视频| 中文字幕熟女人妻在线| 欧美精品啪啪一区二区三区| 高清毛片免费观看视频网站| 亚洲一区二区三区色噜噜| 久久久色成人| 人人妻人人澡欧美一区二区| 精品日产1卡2卡| 高清日韩中文字幕在线| 精品久久久久久久人妻蜜臀av| 乱人视频在线观看| 一个人看的www免费观看视频| 99riav亚洲国产免费| 91久久精品电影网| 狂野欧美激情性xxxx| 久久精品夜夜夜夜夜久久蜜豆| xxx96com| 天天添夜夜摸| 黄色日韩在线| 成人性生交大片免费视频hd| 97碰自拍视频| 亚洲avbb在线观看| 亚洲av成人av| 亚洲无线观看免费| 一边摸一边抽搐一进一小说| 国产免费一级a男人的天堂| 国产成人av教育| 亚洲成人中文字幕在线播放| 国产高清三级在线| 听说在线观看完整版免费高清| 白带黄色成豆腐渣| 不卡一级毛片| 51国产日韩欧美| 亚洲精品成人久久久久久| 国产精品 欧美亚洲| 亚洲av成人av| 又黄又粗又硬又大视频| 午夜免费激情av| 日韩精品青青久久久久久| 在线看三级毛片| 日本 av在线| 中文字幕久久专区| 男女之事视频高清在线观看| 亚洲中文字幕日韩| 黄色片一级片一级黄色片| 在线视频色国产色| 国产高清videossex| 久久精品影院6| 91在线观看av| 男女视频在线观看网站免费| 久久久久久九九精品二区国产| 亚洲一区高清亚洲精品| 美女黄网站色视频| 亚洲无线观看免费| 十八禁人妻一区二区| 男人舔奶头视频| 蜜桃久久精品国产亚洲av| 国产成年人精品一区二区| 搡老熟女国产l中国老女人| 日本 欧美在线| 国内少妇人妻偷人精品xxx网站| 蜜桃亚洲精品一区二区三区| 日本精品一区二区三区蜜桃| 欧美乱色亚洲激情| 青草久久国产| 亚洲欧美日韩无卡精品| 久久久久久久亚洲中文字幕 | 日本黄色片子视频| 丰满的人妻完整版| 国产亚洲精品综合一区在线观看| 国产精品久久视频播放| 此物有八面人人有两片| 高清在线国产一区| 亚洲无线在线观看| 久久国产精品影院| 三级毛片av免费| 精品国产三级普通话版| 国产探花在线观看一区二区| 两个人视频免费观看高清| 亚洲人成网站在线播放欧美日韩| 热99re8久久精品国产| 婷婷六月久久综合丁香| 久久久久久久久中文| 亚洲精品影视一区二区三区av| 日韩人妻高清精品专区| 一卡2卡三卡四卡精品乱码亚洲| 久久久久久久久久黄片| 老司机深夜福利视频在线观看| 给我免费播放毛片高清在线观看| 亚洲欧美日韩无卡精品| 午夜精品一区二区三区免费看| 国产黄色小视频在线观看| 中文亚洲av片在线观看爽| 国产视频内射| 亚洲av美国av| 日本熟妇午夜| 在线观看日韩欧美| 婷婷六月久久综合丁香| 日韩欧美在线乱码| 十八禁人妻一区二区| 久久人人精品亚洲av| 尤物成人国产欧美一区二区三区| 欧美性猛交黑人性爽| 免费av不卡在线播放| 亚洲欧美激情综合另类| 国产精品免费一区二区三区在线| 亚洲一区二区三区不卡视频| 午夜免费激情av| 精品久久久久久,| 夜夜躁狠狠躁天天躁| 少妇裸体淫交视频免费看高清| 国产一级毛片七仙女欲春2| 蜜桃亚洲精品一区二区三区| 少妇的逼好多水| 成年版毛片免费区| 很黄的视频免费| 美女黄网站色视频| 伊人久久精品亚洲午夜| 真人一进一出gif抽搐免费| 欧美bdsm另类| 中文字幕av在线有码专区| 国模一区二区三区四区视频| 免费看a级黄色片| 99热这里只有是精品50| 欧美成人性av电影在线观看| 成人av一区二区三区在线看| 欧美极品一区二区三区四区| 18禁国产床啪视频网站| 精品一区二区三区av网在线观看| 国产黄a三级三级三级人| 色综合婷婷激情| 色吧在线观看| 亚洲在线观看片| 女生性感内裤真人,穿戴方法视频| 天堂av国产一区二区熟女人妻| 舔av片在线| 伊人久久精品亚洲午夜| 母亲3免费完整高清在线观看| 老熟妇仑乱视频hdxx| 亚洲人成网站在线播| 欧美一区二区国产精品久久精品| 欧美区成人在线视频| 狂野欧美激情性xxxx| 成年女人永久免费观看视频| 国产精品98久久久久久宅男小说| 欧美色欧美亚洲另类二区| 五月玫瑰六月丁香| 亚洲色图av天堂| 国产视频内射| 天堂√8在线中文| 亚洲专区国产一区二区| 女人被狂操c到高潮| 九色成人免费人妻av| 又紧又爽又黄一区二区| 国产单亲对白刺激| 久久久久久国产a免费观看| 九九在线视频观看精品| 99久国产av精品| 久久久久久大精品| 国产午夜福利久久久久久| 亚洲人成网站在线播| 黄色视频,在线免费观看| 我的老师免费观看完整版| 国产精品永久免费网站| 日本免费a在线| 中文字幕av成人在线电影| 久久国产精品人妻蜜桃| 亚洲aⅴ乱码一区二区在线播放| 国产精品影院久久| 高潮久久久久久久久久久不卡| 亚洲精品影视一区二区三区av| 无遮挡黄片免费观看| 精品久久久久久久人妻蜜臀av| 免费观看的影片在线观看| 精品国产三级普通话版| 无限看片的www在线观看|