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

    土壤水分-基質(zhì)勢-溫度復(fù)合傳感器研制

    2023-05-15 05:38:36林昱檳向梓薇顏小飛
    農(nóng)業(yè)工程學(xué)報 2023年6期
    關(guān)鍵詞:傳輸線土壤水分標(biāo)定

    林昱檳,程 強(qiáng),向梓薇,顏小飛

    土壤水分-基質(zhì)勢-溫度復(fù)合傳感器研制

    林昱檳1,程 強(qiáng)2,向梓薇1,顏小飛1※

    (1. 北京林業(yè)大學(xué)工學(xué)院,北京 100083;2. 中國農(nóng)業(yè)大學(xué)信息與電氣工程學(xué)院,北京 100083)

    土壤水分、基質(zhì)勢是描述土壤持水特性的重要參數(shù),不同類型土壤的水分特征曲線存在明顯差異。目前在野外環(huán)境下通常采用傳感器分離測量法測定土壤水分特征曲線,但因土壤的空間異質(zhì)性導(dǎo)致較大的測量誤差。為此,該研究設(shè)計了一種可同步測量土壤水分、基質(zhì)勢和溫度的復(fù)合傳感器,主要組成部分包括32位MCU主控芯片、繼電器切換模塊、介電水分/基質(zhì)勢測量模塊和溫度測量模塊,主控芯片通過控制繼電器切換實現(xiàn)土壤水分、基質(zhì)勢和溫度的同步測定,同時進(jìn)行溫度校正。采用森林土壤、農(nóng)田壤土和砂土對復(fù)合傳感器進(jìn)行性能測試、標(biāo)定和觀測試驗,并與商業(yè)化儀器的測量結(jié)果進(jìn)行對比。結(jié)果表明:在干燥至近飽和范圍內(nèi),復(fù)合傳感器測定的土壤水分和基質(zhì)勢均具有良好的單調(diào)性,決定系數(shù)2均大于0.98,土壤體積含水率和基質(zhì)勢的測量范圍分別為0%~40%和-1 500~-15 kPa,水分和基質(zhì)勢的響應(yīng)時間分別為450 ms和150 s。采用商業(yè)化儀器分離測量方法與本文復(fù)合傳感器分別同步測定森林土壤、農(nóng)田壤土和砂土的水分特征曲線,2種方法測得的土壤水分特征曲線的皮爾遜相關(guān)系數(shù)均大于0.96。本文研制的復(fù)合傳感器可以實現(xiàn)不同類型土壤水分特征曲線的準(zhǔn)確測定,為農(nóng)業(yè)生產(chǎn)和環(huán)境監(jiān)測等領(lǐng)域提供可靠的技術(shù)支持。

    傳感器;溫度;土壤;基質(zhì)勢;特征曲線

    0 引 言

    土壤水分特征曲線(soil water characteristic curve,SWCC)描述了土壤水分與吸力(基質(zhì)勢)之間的關(guān)系,是土壤水力特性的重要指標(biāo),反映了土壤水的數(shù)量與能量之間的關(guān)系[1]。通過SWCC可探究不同質(zhì)地土壤的持水特性、植物可利用的有效土壤水分、土壤孔隙結(jié)構(gòu)的分布特征等,同時在利用數(shù)學(xué)物理方法分析土壤中水的運動規(guī)律時具有重要作用[2]。

    目前SWCC的測定主要采用水分和基質(zhì)勢分離測量的方法,土壤水分測量方法主要有烘干法、遙感法、放射法、介電法等[3-6];土壤基質(zhì)勢測量方法主要有張力計法、壓力室法、介電法、離心機(jī)法等[7]。在土壤水分測量方法中,烘干法作為最準(zhǔn)確的測量方法,一般用于實驗室測量,缺點在于對土壤取樣是破壞性的,且耗時較長[8];遙感法通常用于測量大片區(qū)域的土壤表層水分,但是容易受到地表植被覆蓋的影響[9];放射法的測定依據(jù)是土壤中的水分會影響放射線的衰減,測量準(zhǔn)確,其缺點在于放射線對人體有傷害[10];介電法具有操作方便、響應(yīng)速度快、可用于野外原位測量等優(yōu)點,成為目前主要的土壤水分測量方式[11-12]。在土壤基質(zhì)勢測量方法中,張力計可用于田間原位測量,但測量范圍較?。?80~0 kPa),且需要較長的平衡時間[13];壓力室測量最為準(zhǔn)確但平衡時間最長,需要對土壤進(jìn)行破壞性采樣且只能在實驗室中測量[14-15];離心機(jī)法通過設(shè)定轉(zhuǎn)速對土壤施加不同壓力,會導(dǎo)致土壤吸力發(fā)生變化,影響SWCC測量的準(zhǔn)確性,同時也需要破壞性采樣,只適用于實驗室測量[16];介電法可實現(xiàn)原位測定,通過測量與土壤達(dá)成水力平衡的介質(zhì)的基質(zhì)勢來間接獲得土壤基質(zhì)勢,如商用的TEROS 21傳感器[17-18],但其價格比較昂貴。本課題組研制的基于介電法的同軸多孔傳感器[19]可用于土壤基質(zhì)勢的低成本原位測量。

    采用分離測量法測定SWCC會不可避免地因土壤空間異質(zhì)性存在測量誤差,而且現(xiàn)有測量儀器存在價格昂貴、測量范圍小或無法進(jìn)行原位測定等缺點。因此,本文基于介電法研制了一種可同步測量土壤水分、基質(zhì)勢和溫度的復(fù)合傳感器,以降低分別放置傳感器時因土壤空間異質(zhì)性造成的測量誤差,實現(xiàn)SWCC的原位測定。

    1 材料與方法

    1.1 測量原理

    1.1.1 水分測量原理

    水的相對介電常數(shù)(常溫下為81)與干土的相對介電常數(shù)(約為3)差異較大,可通過測定土壤介電常數(shù)變化來反映土壤水分變化[20],因此,本文基于駐波比原理來實現(xiàn)土壤水分和基質(zhì)勢的測定。根據(jù)傳輸線理論,傳輸線阻抗與負(fù)載阻抗相等時,傳輸線上不會產(chǎn)生反射波,當(dāng)傳輸線與負(fù)載阻抗不匹配時,傳輸線上則會產(chǎn)生反射波,入射波與反射波疊加形成駐波。駐波比定義為駐波的波峰與波谷的比值,傳輸線阻抗與負(fù)載阻抗不匹配的程度影響反射波的大小,因此也會影響駐波比的大小。

    本文所設(shè)計的測量電路主要由有源晶體振蕩器、檢波器和放大器組成,如圖1所示。研究表明,土壤介電常數(shù)的測定受測量頻率影響,在低于100 MHz的測量頻率下,土壤中的鹽分會對測量結(jié)果造成較大干擾,而在100~500 MHz的測量頻率范圍內(nèi),測量結(jié)果受到的干擾較小[11,21]。因此,本文采用100 MHz作為測量頻率。在設(shè)計電路時,首先需要進(jìn)行阻抗匹配,使傳輸線的特征阻抗0與置于空氣中的探頭阻抗相等。當(dāng)探頭埋入土壤中,探頭及其周圍土壤的阻抗Z發(fā)生變化,此時0與Z不匹配,傳輸線上產(chǎn)生駐波,通過檢波器檢測電壓差并最終通過放大器輸出為電壓值U

    式中A為放大器的增益。通過放大器放大檢波輸出UU的電壓差值,得到傳感器的輸出電壓。由圖1可知,電壓與阻抗的關(guān)系為

    式中Z為探頭及其周圍土壤的混合阻抗,Ω,為土壤的介電常數(shù)。

    將式(2)代入式(1)可得:

    由式(3)可知,當(dāng)探頭周圍的土壤水分含量發(fā)生變化時,傳感器的輸出電壓U也隨之改變。

    注:U、U為檢波器A、B的電壓,V;U為傳感器的輸出電壓,V;0為傳輸線的特征阻抗,Ω;Z為探頭與周圍土壤的混合阻抗,Ω。

    Note: U、Uis the voltage of detector A and B,V. Uis the output of sensor, V;0is the transmission line characteristic impedance, Ω.Zis the mixed impedance of probe and surrounding soil, Ω.

    圖1 測量電路原理圖

    Fig.1 Schematic of measurement circuit

    1.1.2 基質(zhì)勢測量原理

    土壤基質(zhì)勢的測定同樣基于介電原理,在所設(shè)計的基質(zhì)勢測量探頭中填充水分特征曲線穩(wěn)定的多孔材料?;|(zhì)勢探頭埋置于土壤后,根據(jù)水力平衡原理,土壤與探頭中的多孔材料有勢能差,致使水分在兩者之間運動并最終達(dá)到平衡,此時土壤與探頭中多孔材料的基質(zhì)勢相等,即

    1.1.3 溫度測量原理

    常用的溫度傳感器有數(shù)字式溫度傳感器、集成溫度傳感器、熱電偶以及熱電阻等。熱電阻以其成本低廉、使用簡單、適用范圍廣的特點得到了廣泛應(yīng)用。在測量范圍(-200~850 ℃)內(nèi),PT100的溫度計算式為

    PT100有多種接線方式,其中兩線制最簡單,四線制精度最高,本文綜合接線和精度要求,采用三線接法,同時使用恒流源連接鉑電阻[26]。溫度變化會引起鉑電阻的阻值變化,從而轉(zhuǎn)變?yōu)槲⑷醯碾妷鹤兓偻ㄟ^放大電路以及A/D轉(zhuǎn)換電路和公式換算獲得溫度值。

    1.2 復(fù)合傳感器探頭

    土壤基質(zhì)勢探頭如圖2所示。同軸多孔探頭包含內(nèi)外2個電極,電極材料為不銹鋼,內(nèi)部電極作為發(fā)射極,外部環(huán)形多孔電極作為接收極,內(nèi)外電極之間填充的多孔材料為石膏,土壤中的水分可通過多孔電極與內(nèi)部填充的石膏進(jìn)行交換,從而達(dá)到內(nèi)外的基質(zhì)勢平衡,即電極外部土壤與內(nèi)部石膏的基質(zhì)勢相同。

    所設(shè)計的復(fù)合傳感器如圖2c所示。溫度測量采用PT100鉑電阻。采用探頭復(fù)用的方式進(jìn)行水分測量,將溫度探頭和基質(zhì)勢探頭作為土壤水分測量的2個電極,其中溫度探頭外殼作為發(fā)射極,與基質(zhì)勢探頭外部電極組成類雙針式測量結(jié)構(gòu),可測量溫度探頭與基質(zhì)勢探頭之間的土壤水分。傳感器工作電壓范圍為6~32 V,水分測量電路的輸出信號為模擬信號,輸出電壓為0~2 V,溫度測量范圍-25~45 ℃。

    復(fù)合傳感器的測量過程如下:首先進(jìn)行基質(zhì)勢和溫度測量。基質(zhì)勢探頭帶有中心發(fā)射極和接收極,可直接測量探頭中柱狀石膏的水分,在探頭內(nèi)的石膏與土壤達(dá)到水力平衡后,傳感器的輸出值反映了石膏的水含量,通過標(biāo)定建立傳感器輸出與土壤基質(zhì)勢的對應(yīng)關(guān)系,從而得到獲得土壤基質(zhì)勢。同時,溫度探頭測量此時的土壤溫度。然后,控制切換電路進(jìn)行土壤水分測量,此時溫度探頭作為發(fā)射極,基質(zhì)勢探頭的外殼作為接收極,測量兩個探頭之間的土壤水分,從而完成一次土壤水分、基質(zhì)勢和溫度的測量,并將測量數(shù)據(jù)存儲于TF卡中。

    1.供電模塊 2.控制模塊 3.存儲模塊 4.土壤水分測量電路 5.切換電路 6.土壤基質(zhì)勢探頭 7.溫度探頭 8.土壤水分探頭

    1.3 性能測試試驗

    1.3.1 復(fù)合傳感器的單調(diào)性

    復(fù)合傳感器使用繼電器切換測量通路,切換電路以及探頭分布造成的參數(shù)變化會使傳感器與探頭之間的傳輸線阻抗出現(xiàn)不匹配的現(xiàn)象,從而導(dǎo)致傳感器在測量水分以及基質(zhì)勢的過程中輸出值不單調(diào),影響傳感器測量的準(zhǔn)確性,因此需要設(shè)計不同長度的傳輸線來保證阻抗匹配。本文設(shè)計包含不同長度傳輸線的PCB板,調(diào)整傳感器至探頭間傳輸線的長度,然后依次連接探頭,測試探頭在干燥至濕潤過程中傳感器輸出電壓的變化。若傳感器輸出電壓不單調(diào),則調(diào)整傳輸線長度,直至復(fù)合傳感器的水分和基質(zhì)勢輸出電壓與土壤水分之間呈現(xiàn)單調(diào)性。

    1.3.2 復(fù)合傳感器的標(biāo)定

    為保證測量的準(zhǔn)確性,在實驗室條件下對傳感器進(jìn)行標(biāo)定。本文分別對3種不同類型的土壤(森林土壤、農(nóng)田壤土和砂土)進(jìn)行水分測量標(biāo)定。其中,森林土壤采集自北京市鷲峰國家森林公園(39°54′N,116°28′E)中鵝耳櫟、油松和栓皮櫟混交林,農(nóng)田壤土和砂土采集自中國農(nóng)業(yè)大學(xué)甘肅石羊河試驗站(45°93′N,114°83′E)。首先使用粗篩(孔徑6 mm)將土壤中的石塊、樹枝等雜物濾除,然后使用細(xì)篩(孔徑2 mm)將土壤過篩,在105 ℃下烘干24 h。

    取640 g烘干后的土壤,加入32 g水,攪拌均勻,分層放入測試容器中(長16.2 cm,寬7 cm,高5 cm)并逐層壓實,同時將復(fù)合傳感器埋置于土壤中,測量土壤體積,計算質(zhì)量含水率和土壤容重,進(jìn)而得到待測土壤的體積含水率。用保鮮膜覆蓋測試容器并靜置4 h,記錄水分探頭的輸出電壓。然后取出土壤并再次加入32 g水拌勻,重新分層放入測量容器中進(jìn)行水分測量,如此往復(fù)直至土壤水分接近飽和狀態(tài),完成各土壤樣本的水分標(biāo)定。

    采用商品化的基質(zhì)勢傳感器(TEROS 21,美國METER公司)對本文研制的復(fù)合傳感器的基質(zhì)勢測量模塊進(jìn)行標(biāo)定。TEROS 21的測量范圍為-2 000~-9 kPa,其中-100~-9 kPa范圍內(nèi)測量精度為±(讀數(shù)的10%+2 kPa);在-1 500~-100 kPa范圍內(nèi),TEROS 21也具有良好的測量準(zhǔn)確性[27-29]。稱取640 g烘干后的土壤,分層放入測試容器中(長16.2 cm,寬7 cm,高5 cm)并逐層壓實,將2個TEROS 21傳感器與本文研制的復(fù)合傳感器埋置于同一土樣平面。采用滴灌方式濕潤土壤,滴灌點分別位于TEROS21與復(fù)合傳感器的上方,滴灌流速為0.32 L/h,直至土壤水分達(dá)到近飽和狀態(tài),使用保鮮膜覆蓋測試容器并靜置4 h,等待土壤中水分運移趨于穩(wěn)定。之后使用風(fēng)扇風(fēng)干土壤,分別記錄土壤在濕潤至干燥過程中復(fù)合傳感器與TEROS 21的輸出值,重復(fù)此過程2次,完成基質(zhì)勢標(biāo)定。

    溫度標(biāo)定使用高精度溫度計(型號DT321S,深圳華盛昌)。將復(fù)合傳感器和DT321S同時埋置于土箱中,并將土箱放置于恒溫箱中,設(shè)置恒溫箱初始溫度為-15 ℃,每90 min溫度提高5 ℃,考慮到復(fù)合傳感器的使用環(huán)境以及野外的環(huán)境溫度,本文在-15~45 ℃的溫度范圍內(nèi)對復(fù)合傳感器的溫度測量模塊進(jìn)行校準(zhǔn)。

    1.3.3 復(fù)合傳感器的有效測量半徑

    將復(fù)合傳感器探頭埋置于土壤中,從表面逐層移除土壤,每次移除土壤厚度約3 mm,記錄并對比每次移除土壤后傳感器的輸出電壓,得到傳感器在水分和基質(zhì)勢檢測中的有效測量半徑。當(dāng)移除土壤后傳感器輸出電壓發(fā)生較大變化時,則認(rèn)為此時覆蓋在傳感器上方的土層厚度即為傳感器的有效測量半徑。

    1.3.4 傳感器的響應(yīng)時間

    在林區(qū)或農(nóng)田環(huán)境中,往往因為降水或灌溉導(dǎo)致土壤含水量快速變化。當(dāng)土壤含水量發(fā)生變化時,基質(zhì)勢探頭中的石膏需要一定時間與土壤建立水勢平衡。本文通過試驗測定干燥的復(fù)合傳感器探頭埋入飽和土壤中的響應(yīng)時間。首先將探頭置于烘干箱中烘干,配置壤土飽和土樣,隨后快速將復(fù)合傳感器探頭插入濕潤土樣中,并記錄傳感器的輸出。由于水分測量的響應(yīng)時間較短,本文將傳感器的輸出接入示波器,水勢探頭響應(yīng)時間通過測量復(fù)合傳感器的輸出電壓變化特征,采樣時間設(shè)置為1 s。

    1.3.5 傳感器的穩(wěn)定性

    將傳感器探頭放置于容器中,填滿濕潤的土壤,并使用保鮮膜將容器封好,防止水分揮發(fā)。靜置4 h后,每隔1 min采集72 h內(nèi)復(fù)合傳感器的水分、基質(zhì)勢、溫度輸出值。

    1.3.6 傳感器的溫度效應(yīng)

    在環(huán)境溫度發(fā)生變化時,測量輸出值也會隨之發(fā)生變化。該影響主要體現(xiàn)在電路板的溫度系數(shù)和介電常數(shù)的溫度效應(yīng)[19]。為了降低介電常數(shù)溫度效應(yīng)的影響,根據(jù)溫度對多孔材料影響的模型[19],得到水分和基質(zhì)勢的介電常數(shù)的溫度效應(yīng)校正系數(shù)分別為0.55和0.71。為了消除電路板的溫度系數(shù)對測量準(zhǔn)確性的影響,利用復(fù)合傳感器的溫度探頭對電路板的溫度進(jìn)行檢測,同時校正復(fù)合傳感器的輸出值。將復(fù)合傳感器放置于恒溫箱中,在-25~45 ℃的溫度范圍內(nèi),檢測并記錄傳感器的空載輸出,根據(jù)測試數(shù)據(jù)擬合結(jié)果校正傳感器的輸出值。

    1.4 不同類型土樣的SWCC測定

    為探究復(fù)合傳感器測量SWCC的效果,分別使用復(fù)合傳感器與傳統(tǒng)分離測量方法在實驗室環(huán)境中同時繪制森林土壤、農(nóng)田壤土和砂土的SWCC。分離測量方法使用稱重法測量土壤水分含量;在基質(zhì)勢測量中,使用TEROS 21同步測量土壤基質(zhì)勢。

    具體步驟如下:

    1)在測量容器中填入待測試的干燥土壤,并分別埋入復(fù)合傳感器探頭以及TEROS21基質(zhì)勢傳感器,二者保持相同高度,在測試容器下放置高精度電子秤,并連接采集盒定時采集數(shù)據(jù)。

    2)將干燥土壤滴灌(滴灌速度0.32 L/h)至近飽和狀態(tài),然后風(fēng)干至干燥狀態(tài)。分別獲得復(fù)合傳感器與分離測量法測得的數(shù)據(jù),并繪制SWCC。完成該類型土樣的SWCC繪制后,清理測量容器并填入另一種類型的干燥土壤。

    3)重復(fù)以上過程,完成3種土樣的SWCC測定并進(jìn)行對比。

    2 結(jié)果與分析

    2.1 單調(diào)性測試結(jié)果

    傳輸線長會直接影響反射波的相位,從而影響駐波的大小和變化趨勢,最終影響傳感器的輸出值。為保證整個測量過程中傳感器的輸出值保持單調(diào),本研究在探頭由干燥至濕潤的過程中,對不同傳輸線長的傳感器輸出進(jìn)行測試(圖3)。當(dāng)傳輸線長分別為4、6、12、15和24 cm時,傳感器在探頭干燥時的輸出值分別為0.013、0.026、0.246、0.697和1.842 V,可見,探頭干燥狀態(tài)下,傳感器輸出隨著傳輸線長度的增加而增大。在探頭濕潤過程中,傳輸線長為4、6、12和15 cm的傳感器輸出均呈現(xiàn)先增大后減小的趨勢,不滿足單調(diào)性需求;而傳輸線長為24 cm的傳感器輸出電壓由探頭干燥時的1.82 V逐漸降低至探頭完全濕潤時的0.013 V,具有良好的單調(diào)性。因此,本文傳感器的傳輸線長設(shè)定為24 cm。

    圖3 傳感器的單調(diào)性測試結(jié)果

    2.2 傳感器標(biāo)定結(jié)果

    分別使用森林土壤、農(nóng)田壤土和砂土對復(fù)合傳感器進(jìn)行標(biāo)定,控制容重和加水量,得到不同體積含水率的土壤樣本,記錄傳感器的輸出電壓,并對其進(jìn)行擬合。

    土壤樣本的水分標(biāo)定結(jié)果如圖4所示,傳感器輸出電壓與3種土壤的體積含水率之間呈現(xiàn)良好的線性關(guān)系,決定系數(shù)2均大于0.98,土壤體積含水率測量范圍為0%~40%。

    圖4 土壤水分含量標(biāo)定結(jié)果

    根據(jù)測量原理可知,本文基質(zhì)勢的測量僅與探頭內(nèi)部柱形石膏的結(jié)構(gòu)、孔隙度等因素有關(guān),與電極外部待測土壤的類型無關(guān),因此本文只采用森林土壤進(jìn)行基質(zhì)勢的標(biāo)定。土壤中放置復(fù)合傳感器與TEROS 21,在土壤由飽和到干燥的過程中,復(fù)合傳感器的基質(zhì)勢輸出電壓值與TEROS 21基質(zhì)勢輸出值的關(guān)系如圖5所示。根據(jù)使用手冊[30]可知,TEROS 21在-100~-9 kPa范圍內(nèi)的測量精度為±(讀數(shù)的 10%+2 kPa),在-1 500~-100 kPa范圍內(nèi)準(zhǔn)確性良好。因此,本研究使用分段擬合的方法在-100~-9 kPa、-260~-100 kPa以及小于-260 kPa三個范圍對復(fù)合傳感器的基質(zhì)勢測量范圍進(jìn)行分別標(biāo)定,由圖5可知,復(fù)合傳感器的輸出電壓與TEROS21測得的基質(zhì)勢具有良好的線性關(guān)系,3個測量范圍內(nèi)的擬合2分別為0.988 7、0.979 9和0.972 6,復(fù)合傳感器的基質(zhì)勢測量范圍為-1 500~-15 kPa。

    圖5 土壤基質(zhì)勢標(biāo)定結(jié)果

    將高精度溫度計與復(fù)合傳感器埋置于土箱中,并將土箱置于恒溫箱中,在-15~45 ℃范圍內(nèi)調(diào)節(jié)溫度。恒溫箱初始溫度設(shè)置為-15 ℃,每30 min將溫度調(diào)高5 ℃,記錄溫度穩(wěn)定后的高精度溫度計(標(biāo)準(zhǔn)溫度)與復(fù)合傳感器溫度模塊的測試結(jié)果,如圖6所示。可以看出,復(fù)合傳感器溫度模塊具有良好的測量精度,決定系數(shù)2為0.99,滿足測量要求。

    圖6 溫度標(biāo)定結(jié)果

    2.3 傳感器的測量半徑

    傳感器測量半徑測試結(jié)果如圖7所示,當(dāng)土層表面距離探頭17 mm時,土壤水分測量模塊的輸出電壓發(fā)生較大變化,因此,復(fù)合傳感器水分測量模塊的有效測量半徑為17 mm。由于基質(zhì)勢探頭使用了同軸多孔的設(shè)計結(jié)構(gòu),其測量半徑只局限在探頭內(nèi)部,不受外部土層移除的影響。

    圖7 傳感器有效測量半徑測試結(jié)果

    2.4 傳感器的響應(yīng)時間

    復(fù)合傳感器的響應(yīng)時間如圖8~圖9所示。由圖8可知,從傳感器開始響應(yīng)直至輸出電壓穩(wěn)定的時間約為450 ms。由圖9可知,在70 s左右,傳感器輸出電壓開始出現(xiàn)大幅快速變化,在100 s時,傳感器輸出電壓變化速度放緩,并在150 s時到達(dá)平衡。

    圖8 土壤水分探頭響應(yīng)時間

    圖9 土壤基質(zhì)勢探頭響應(yīng)時間

    試驗開始時基質(zhì)勢探頭與飽和土壤的基質(zhì)勢差距很大,在60~100 s時,基質(zhì)勢探頭處于快速吸水的階段,隨著基質(zhì)勢探頭與飽和土壤的水勢差減小,在100~150 s時基質(zhì)勢探頭的吸水速度減慢;在150 s后,基質(zhì)勢探頭與飽和土壤達(dá)到水勢動態(tài)平衡,傳感器輸出變化近似為0。全過程基質(zhì)勢探頭的水勢變化趨勢是先快后慢的,基質(zhì)勢探頭的響應(yīng)時間為150 s。

    2.5 傳感器的穩(wěn)定性

    將復(fù)合傳感器埋入濕潤土壤中靜置4 h,待土壤中水分運移完成后開始進(jìn)行復(fù)合傳感器的穩(wěn)定性測試,此后的72 h內(nèi),復(fù)合傳感器的基質(zhì)勢、水分輸出值保持穩(wěn)定,變化幅度均小于10 mV(圖10)。環(huán)境溫度由于空調(diào)的緣故有所起伏波動。結(jié)果可體現(xiàn)復(fù)合傳感器具有良好的穩(wěn)定性。

    2.6 溫度效應(yīng)校正

    由于溫度變化會導(dǎo)致電路相關(guān)參數(shù)的改變,進(jìn)而影響傳感器輸出,復(fù)合傳感器設(shè)計的溫度模塊可以準(zhǔn)確測量環(huán)境溫度,為土壤水分和基質(zhì)勢測量進(jìn)行溫度校正。為探究溫度與傳感器輸出電壓的關(guān)系,將傳感器空載放入恒溫箱中,在-25~45 ℃的范圍內(nèi),每隔15 min向上調(diào)整一次溫度,記錄傳感器空載輸出電壓發(fā)生的變化,結(jié)果如圖11,在空載條件下,溫度由-25 ℃增加至45 ℃時,傳感器的輸出電壓下降160 mV。根據(jù)經(jīng)驗,溫度上升時,輸出電壓會隨之升高,但復(fù)合傳感器接入的傳輸線使得其輸出呈負(fù)相關(guān),傳感器輸出電壓隨溫度升高而下降。根據(jù)不同溫度下傳感器實際輸出值相比于25 ℃標(biāo)定溫度下的偏離量,在最終輸出值加上當(dāng)前溫度產(chǎn)生的偏移量,可用于低溫或高溫環(huán)境下的溫度效應(yīng)校正。

    圖10 傳感器的穩(wěn)定性測試結(jié)果

    2.7 土壤水分特征線的測試結(jié)果

    SWCC具有溫度效應(yīng),在不同溫度下有所差異[31-32],因此本文在25 ℃下對土壤水分特征線進(jìn)行測定。分別使用商業(yè)化傳感器和本文研制的復(fù)合傳感器測定森林土壤、農(nóng)田壤土和砂土的SWCC,如圖12所示。由圖12可知,2種方法測得的SWCC具有較好的一致性。試驗過程的土壤變化是干燥-濕潤-干燥,由于SWCC具有滯后性[32-33],3種土壤的SWCC在相同水吸力情況下,干燥過程的土壤含水量大于濕潤過程的土壤含水量。對于森林土壤,2種方法測定的基質(zhì)勢與水分的皮爾遜相關(guān)系數(shù)分別為0.962與0.983;對于農(nóng)田壤土,2種方法得到的基質(zhì)勢與水分測量結(jié)果的皮爾遜相關(guān)系數(shù)分別為0.993與0.995;對于砂土,2種方法測得的基質(zhì)勢與水分的皮爾遜相關(guān)系數(shù)分別為0.979和0.998。結(jié)果表明,采用復(fù)合傳感器測量土壤水分特征曲線獲得的結(jié)果與傳統(tǒng)分離測量方法的結(jié)果相近,可代替分離測量方法實現(xiàn)不同類型土樣SWCC的測定。

    圖11 溫度效應(yīng)測試結(jié)果

    a. 森林土壤 a. Forest soilb. 農(nóng)田壤土 b. Farmland loamc. 砂土 c. Sandy soil

    3 結(jié) 論

    本文設(shè)計了一種同步觀測土壤水分、基質(zhì)勢和溫度的復(fù)合傳感器,主要研究結(jié)果如下:

    1)復(fù)合傳感器水分和基質(zhì)勢標(biāo)定后的決定系數(shù)2分別大于0.98和0.97,在測量范圍內(nèi)具有良好的線性輸出。

    2)復(fù)合傳感器對于土壤水分、基質(zhì)勢、溫度的測量范圍分別為0~0.4 cm3/cm3與-1 500~-15 kPa和-15~45 ℃,同軸結(jié)構(gòu)設(shè)計的基質(zhì)勢探頭測量內(nèi)外電極之間包裹的多孔石膏,水分探頭的測量半徑為17 mm。

    3)將干燥的復(fù)合傳感器探頭埋置于飽和土壤中,測得水分探頭的響應(yīng)時間為450 ms,基質(zhì)勢探頭的響應(yīng)時間為150 s。

    4)采用復(fù)合傳感器分別測定森林土壤、農(nóng)田壤土和砂土的土壤水分特征曲線(soil water characteristic curv,SWCC),與使用商業(yè)化傳感器的分離測量方法獲得的結(jié)果具有較好的一致性(皮爾遜相關(guān)系數(shù)大于0.96),本文研制的復(fù)合傳感器可代替分離測量的方法,實現(xiàn)SWCC的原位測量。

    [1] 韓祥偉,邵明安,王全九. 簡單入滲法在確定Brooks-Corey水分特征曲線模型參數(shù)中的應(yīng)用研究[J]. 土壤學(xué)報,2006, 43(3):506-508. HAN Xiangwei, SHAO Ming′an, WANG Quanjiu. Application study of simple infiltration method for determining the parameters of the Brooks-Corey model of soil water retention curve[J]. Acta Pedologica Sinica, 2006, 43(3): 506-508. (in Chinese with English abstract)

    [2] 王愿斌,王佳銘,樊媛媛,等. 土壤水分特征曲線模型模擬性能評價[J]. 冰川凍土,2019,41(6):1448-1455. WANG Yuanbin, WANG Jiaming, FAN Yuanyuan, et al. Performance evaluation of 12 models describing the soil water retention characteristics[J]. Journal of Glaciology and Geocryology, 2019, 41(6): 1448-1455. (in Chinese with English abstract)

    [3] 謝鵬宇,劉澤鑫. 土壤水分測量原理與技術(shù)方法研究[J]. 現(xiàn)代農(nóng)業(yè)科技,2020(23):166-168. XIE Pengyu, LIU Zexin. Study on principles and technical methods of soil moisture measurement[J]. Modern Agricultural Science and Technology, 2020(23): 166-168. (in Chinese with English abstract)

    [4] 邵長亮,吳東麗. 土壤水分測量方法適用性綜述[J]. 氣象科技,2019,47(1):1-9. SHAO Changliang, WU Dongli. Review of applicability of soil moisture measurement techniques[J]. Meteorological Science and Technology, 2019, 47(1): 1-9. (in Chinese with English abstract)

    [5] DOBRIYAL P, QURESHI A, BADOLA R, et al. A review of the methods available for estimating soil moisture and its implications for water resource management[J]. Journal of Hydrology, 2012, 458-459: 110-117.

    [6] ROBINSON D A, CAMPBELL C S, HOPMANS J W, et al. Soil moisture measurement for ecological and hydrological watershed-scale observatories: A review[J]. Vadose Zone Journal, 2008, 7(1): 358-389.

    [7] 萬書勤,李曉彬,康躍虎. 基于土壤基質(zhì)勢的溫室土壤墑情監(jiān)測技術(shù)[J]. 農(nóng)業(yè)工程技術(shù),2019,39(13):10-17.

    [8] 張學(xué)禮,胡振琪,初士立. 土壤含水量測定方法研究進(jìn)展[J]. 土壤通報,2005, 36(1):118-123. ZHANG Xueli, HU Zhenqi, CHU Shili. Methods for measuring soil water content: A review[J]. Chinese Journal of Soil Science, 2005, 36(1): 118-123. (in Chinese with English abstract)

    [9] 馬春芽,王景雷,黃修橋. 遙感監(jiān)測土壤水分研究進(jìn)展[J]. 節(jié)水灌溉,2018(5):70-74, 78.

    [10] JAYAWARDANE N S, MEYER W S, BARRS H D. Moisture measurement in a swelling clay soil using neutron moisture meters[J]. Soil Research, 1984, 22(2): 109-117.

    [11] 孫宇瑞,汪懋華,趙燕東. 一種基于駐波比原理測量土壤介電常數(shù)的方法[J]. 農(nóng)業(yè)工程學(xué)報,1999, 15(2):43-47. SUN Yurui, WANG Maohua, ZHAO Yandong. A kind of determinations of soil dielectric constant using the principle of standing-wave ratio[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 1999, 15(2): 43-47. (in Chinese with English abstract)

    [12] 程強(qiáng). 凍土冰/水含量同腔共射原位測定方法創(chuàng)新與凍融模型檢驗[D]. 北京:中國農(nóng)業(yè)大學(xué),2014. CHENG Qiang. A Co-Tube-Radiation Method for In-Situ Determination of Ice and Liquid Water Contents in Frozen Soil and Freeze-Thaw Model Validation[D]. Beijing: China Agricultural University, 2014. (in Chinese with English abstract)

    [13] TSAI Y Z, CHENG M L, HUANG Q Z, et al. Effect of effective saturation and ceramic cup properties on the response time of tensiometers[J]. Journal of Hydrology, 2020, 582: 124445.

    [14] BULUT R, LEONG E C. Indirect Measurement of Suction[M]. Dordrecht: Springer Netherlands, 2009: 21-32.

    [15] CRESSWELL H P, GREEN T W, MCKENZIE N J. The adequacy of pressure plate apparatus for determining soil water Retention[J]. Soil Science Society of America Journal, 2008, 72(1): 41-49.

    [16] 聶坤堃,聶衛(wèi)波,馬孝義. 離心機(jī)法測定土壤水分特征曲線中的收縮特性[J]. 排灌機(jī)械工程學(xué)報,2019,37(11):978-985. NIE Kunkun, NIE Weibo, MA Xiaoyi. Soil shrinkage property during soil water characteristic curves being measured by centrifugal method[J]. Journal of Drainage and Irrigation Machinery Engineering, 2019, 37(11): 978-985. (in Chinese with English abstract)

    [17] PERSSON M, WRAITH J M, DAHLIN T. A small-scale matric potential sensor based on time domain reflectometry[J]. Soil Science Society of America Journal, 2006, 70(2): 533-536.

    [18] CRESSWELL H P, Lilly A. A two-point method for determining the soil water characteristic of typical northern temperate boreal soils[J]. Geoderma Regional, 2015, 5: 71-76.

    [19] CHEN C, YAN X, MA Y, et al. Monitoring near-surface soil water content using an innovative perforated cylinder coaxial dielectric sensor[J]. Journal of Hydrology, 2019, 573: 746-754.

    [20] 向陽,于淞,徐嬙,等. 太陽能供電的土壤剖面水分動態(tài)原位自動監(jiān)測系統(tǒng)的研制[J]. 農(nóng)業(yè)工程學(xué)報,2021,37(7):150-157. XIANG Yang, YU Song, XU Qiang, et al. Development of a novel solar-powered system for in-situ automatic monitoring of soil water dynamics in a soil profile[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(7): 150-157. (in Chinese with English abstract)

    [21] 趙燕東,王一鳴. 基于駐波率原理的土壤含水率測量方法[J]. 農(nóng)業(yè)機(jī)械學(xué)報,2002, 33(4):109-111, 121. ZHAO Yandong, WANG Yiming. Study on the measurement of soil water content based on the principle of standing-wave ratio[J]. Transactions of the Chinese Society for Agricultural Machinery , 2002, 33(4): 109-111, 121. (in Chinese with English abstract)

    [22] OR D A, WRAITH J M. A new soil matric potential sensor based on time domain reflectometry[J]. Water Resources Research, 1999, 35(11): 3399-3407.

    [23] WHALLEY W R, CLARK L J, TAKE W A, et al. A porous-matrix sensor to measure the matric potential of soil water in the field[J]. European Journal of Soil Science, 2007, 58(1): 18-25.

    [24] CHEN C, YAN X, XU Q, et al. A novel frequency domain impedance sensor with a perforated cylinder coaxial design for in-situ measuring soil matric potential[J]. Sensors, 2019, 19(11): 2626.

    [25] 張修太,胡雪惠,翟亞芳,等. 基于PT100的高精度溫度采集系統(tǒng)設(shè)計與實驗研究[J]. 傳感技術(shù)學(xué)報,2010,23(6):812-815. ZHANG Xiutai, HU Xuehui, ZHAI Yafang, et al. The design and experimental study of high precise temperature acquisition system based on pt100[J]. Chinese Journal of Sensors and Actuators, 2010, 23(6): 812-815. (in Chinese with English abstract)

    [26] 郭少朋,吳嵐軍,劉鴻. 鉑電阻多點測溫控溫設(shè)計[J]. 傳感器技術(shù),2004, 23(1):67-69. GUO Shaopeng, WU Lanjun, LIU Hong. Design of multipoint temperature measurement and control system by platinum resistor[J]. Journal of Transducer Technology, 2004, 23(1): 67-69. (in Chinese with English abstract)

    [27] WALTHERT L, SCHLEPPI P. Equations to compensate for the temperature effect on readings from dielectricMPS-2 and MPS-6 water potential sensors in soils[J]. Journal of Plant Nutrition and Soil Science, 2018, 181(5): 749-759.

    [28] ZEITOUN R, VANDERGEEST M, VASAVA H B, et al. In-situ estimation of soil water retention curve in silt loam and loamy sand soils at different soil depths[J]. Sensors, 2021, 21(2): 447.

    [29] SHARMA K, IRMAK S, KUKAL M S, et al. Evaluation of soil moisture sensing technologies in silt loam and loamy sand soils: Assessment of performance, temperature sensitivity, and site- and sensor-specific calibration functions[J]. Transactions of the ASABE, 2020, 64(4): 1123-1139.

    [30] Meter Inc. TEROS21_Gen1_Manual[EB/OL]. (2020-06-30) [2023-02-25]. http://publications.metergroup.com/Manuals/ 20428_TEROS21_Gen1_Manual_Web.pdf

    [31] 陳宸,張灃,張華,等. 溫度效應(yīng)對粉砂土壤水分特征曲線的影響研究[J]. 中國農(nóng)村水利水電,2016(6):173-176,181. CHEN Chen, ZHANG Feng, ZHANG Hua, et al. Research on temperature effect on soil-water characteristic curve of sandy soils[J]. China Rural Water and Hydropower, 2016(6): 173-176, 181. (in Chinese with English abstract)

    [32] 秦冰,陳正漢,孫發(fā)鑫,等. 高吸力下持水曲線的溫度效應(yīng)及其吸附熱力學(xué)模型[J]. 巖土工程學(xué)報,2012,34(10):1877-1886. QIN Bing, CHEN Zhenghan, SUN Faxin, et al. Temperature effect on water retention curve under high suction and its modeling based on thermodynamics of sorption[J]. Chinese Journal of Geotechnical Engineering, 2012, 34(10): 1877-1886. (in Chinese with English abstract)

    [33] 尹玉鳳,張妙仙,沈海波,等. 土壤水分特征曲線滯后性與比水密度擬合研究[J]. 灌溉排水學(xué)報,2013,32(6):50-53. YIN Yufeng, ZHANG Miaoxian, SHEN Haibo, et al. Fitting research on hysteresis phenomenon of soil water characteristic curves and soil water capacity[J]. Journal of Irrigation and Drainage, 2013, 32(6): 50-53. (in Chinese with English abstract)

    Development of a composite sensor for the soil moisture, matric potential and temperature synchronous measurement

    LIN Yubin1, CHENG Qiang2, XIANG Ziwei1, YAN Xiaofei1※

    (1.,,100083,;2.,,100083,)

    Soil moisture content and matric potential are two important parameters for determining soil water characteristic curve (SWCC), which has significant differences for different types of soils. The common method to determine SWCC under field condition is using two separated sensors to measure moisture content and matric potential of the tested soils. However, measurement error may occur due to soil spatial heterogeneity around two sensors. In this study, a novel composite sensor was designed to simultaneously measure soil moisture content, matric potential, and temperature. The composite sensor consists of a 32-bit MCU main control chip, relay switching module, dielectric measurement module for determining soil moisture content and matric potential, and temperature measurement module. The main MCU control chip measures soil moisture content, matric potential, and temperature through controlling the relay switching, and performs temperature correction as well. The performances of the composite sensor were tested for soil moisture content and soil matric potential measurements including monotonicity, calibration, volume of sensitivity, response time, etc. To test the monotonicity of the dielectric mode, the length of high frequency transmission line was gradually increased and the sensor output was recorded when the probes of soil moisture content and soil water potential were moistened from dry to saturation. The three modes, soil moisture content, soil matric potential and temperature were calibrated respectively using forest soil, farmland loam and sandy soil. To determine the volume of sensitivity of the soil moisture content mode, the probe of soil moisture content was installed in a soil sample. The soil sample was gradually removed layer by layer and meanwhile the output of the soil moisture content was recorded. When the output started to decrease, the radius of volume of sensitivity was determined by measuring the thickness of residual soil layer from surface to the probe. To determine the response time of soil moisture content and soil matric potential modes, the composite sensor probe was firstly dried and then put into a saturated soil sample. Meanwhile, the output of the soil moisture content and soil matric potential was recorded as time elapse. When the recorded output of each mode did not remarkedly change, the elapsed time was the response time. Experiments for sensor observation were conducted using forest soil, farmland loam and sandy soil under laboratory condition at 25 ℃. The experimental results were compared with those measured using commercial instruments. The results showed that the composite sensor had a good monotonicity when measuring soil moisture content and matric potential in the range from dry to near saturation. After calibration, the composite sensor can determine soil moisture content, matric potential and temperature with high accuracies, achieving the determination coefficients larger than 0.98. The measurement ranges of volumetric soil water content and soil matric potential were 0%-40% and -1 500--15 kPa, respectively, which is feasible for most soil textures. The response time of the soil moisture content and soil matric potential was 450 ms and 150 s, respectively. The temperature coefficient of the circuit board of the composite sensor was yielded to correct the temperature induced output variation. The SWCCs of forest soil, farmland loam and sandy soil were determined using both the commercial instruments and the composite sensor during wetting and drying processes. The hysteresis between soil drying and wetting processes was also observed as expected. The Pearson correlation coefficients between the SWCCs measured by the two methods were larger than 0.96 during wetting and drying processes, i.e., 0.962 and 0.983 for the forest soil, 0.993 and 0.995 for the farmland loam, and 0.979 and 0.998 for the sandy soil. The composite sensor developed in this study can accurately determine soil water characteristic curves of different types of soils.

    sensor; temperature; soil; matric potential; characteristic curve

    10.11975/j.issn.1002-6819.202211187

    S237

    A

    1002-6819(2023)-06-0035-09

    林昱檳,程強(qiáng),向梓薇,等. 土壤水分-基質(zhì)勢-溫度復(fù)合傳感器研制[J]. 農(nóng)業(yè)工程學(xué)報,2023,39(6):35-43.doi:10.11975/j.issn.1002-6819.202211187 http://www.tcsae.org

    LIN Yubin, CHENG Qiang, XIANG Ziwei, et al. Development of a composite sensor for the soil moisture, matric potential and temperature synchronous measurement[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2023, 39(6): 35-43. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.202211187 http://www.tcsae.org

    2022-11-22

    2023-03-02

    國家自然科學(xué)基金資助項目(31971576)

    林昱檳,研究方向為林業(yè)信息檢測。Email:lybfjha@163.com

    顏小飛,博士,副教授,研究方向為林業(yè)生態(tài)信息智能監(jiān)測。Email:yanxf@bjfu.edu.cn

    猜你喜歡
    傳輸線土壤水分標(biāo)定
    使用朗仁H6 Pro標(biāo)定北汽紳寶轉(zhuǎn)向角傳感器
    基于勻速率26位置法的iIMU-FSAS光纖陀螺儀標(biāo)定
    西藏高原土壤水分遙感監(jiān)測方法研究
    船載高精度星敏感器安裝角的標(biāo)定
    略談λ/4傳輸線的實際應(yīng)用
    基于Harris-張正友平面標(biāo)定法的攝像機(jī)標(biāo)定算法
    不同覆蓋措施對棗園土壤水分和溫度的影響
    植被覆蓋區(qū)土壤水分反演研究——以北京市為例
    土壤水分的遙感監(jiān)測方法概述
    不等長非均勻有損耗傳輸線FDTD瞬態(tài)分析
    久热爱精品视频在线9| 国产极品粉嫩免费观看在线| 母亲3免费完整高清在线观看| 亚洲成人国产一区在线观看| 国产在线一区二区三区精| 搡老熟女国产l中国老女人| 久久青草综合色| 免费av中文字幕在线| 中文精品一卡2卡3卡4更新| 男女边摸边吃奶| 天堂中文最新版在线下载| 国产伦人伦偷精品视频| 性色av乱码一区二区三区2| 国内毛片毛片毛片毛片毛片| 肉色欧美久久久久久久蜜桃| 人人妻人人澡人人爽人人夜夜| 精品人妻熟女毛片av久久网站| 深夜精品福利| 免费少妇av软件| 久久狼人影院| 考比视频在线观看| 啦啦啦免费观看视频1| 午夜影院在线不卡| 国产高清视频在线播放一区 | 1024香蕉在线观看| 久久毛片免费看一区二区三区| 久久精品国产亚洲av香蕉五月 | 久久久精品国产亚洲av高清涩受| 一级,二级,三级黄色视频| 伊人亚洲综合成人网| 永久免费av网站大全| cao死你这个sao货| 99国产精品一区二区蜜桃av | 99国产精品一区二区三区| 丰满少妇做爰视频| 国产成人精品久久二区二区免费| 真人做人爱边吃奶动态| 日本av手机在线免费观看| 2018国产大陆天天弄谢| 亚洲第一欧美日韩一区二区三区 | 精品一区二区三卡| 国产淫语在线视频| 日韩电影二区| 精品一区二区三卡| 天天影视国产精品| 女人高潮潮喷娇喘18禁视频| 欧美日韩亚洲高清精品| 在线亚洲精品国产二区图片欧美| 自线自在国产av| 国产在线观看jvid| 人妻 亚洲 视频| 国产精品二区激情视频| 国产一区二区三区综合在线观看| 国产精品九九99| 19禁男女啪啪无遮挡网站| 国产男女超爽视频在线观看| 国产精品av久久久久免费| a级毛片在线看网站| a在线观看视频网站| 香蕉国产在线看| 国产欧美日韩一区二区精品| 国产97色在线日韩免费| 亚洲国产av新网站| 老司机福利观看| 男女国产视频网站| 久久国产精品人妻蜜桃| 国产精品一区二区在线不卡| 我的亚洲天堂| 又黄又粗又硬又大视频| 黑人操中国人逼视频| xxxhd国产人妻xxx| 午夜日韩欧美国产| 五月天丁香电影| 久久精品成人免费网站| 久久中文字幕一级| 久久精品国产亚洲av高清一级| 久久精品成人免费网站| tocl精华| 一本—道久久a久久精品蜜桃钙片| 国产激情久久老熟女| xxxhd国产人妻xxx| 高清欧美精品videossex| 一级a爱视频在线免费观看| 性色av一级| 午夜免费成人在线视频| 精品亚洲成国产av| 一区福利在线观看| 99精品欧美一区二区三区四区| 久久精品人人爽人人爽视色| 视频在线观看一区二区三区| 制服诱惑二区| 亚洲免费av在线视频| 9色porny在线观看| 欧美乱码精品一区二区三区| 国产精品一区二区精品视频观看| 久久国产精品影院| 国产一级毛片在线| 乱人伦中国视频| 日韩 欧美 亚洲 中文字幕| 亚洲成av片中文字幕在线观看| 成年人免费黄色播放视频| 丝袜美腿诱惑在线| 色精品久久人妻99蜜桃| 亚洲熟女精品中文字幕| 色婷婷av一区二区三区视频| 亚洲国产中文字幕在线视频| 高清在线国产一区| 午夜福利视频精品| 欧美日韩av久久| 国产精品久久久久久精品古装| av线在线观看网站| 久久久精品免费免费高清| 国产一区二区三区综合在线观看| 99国产精品一区二区蜜桃av | 日韩视频在线欧美| 五月开心婷婷网| a级毛片在线看网站| 国产一区二区在线观看av| 国产片内射在线| 窝窝影院91人妻| 一级a爱视频在线免费观看| 韩国精品一区二区三区| 国产精品久久久人人做人人爽| 岛国毛片在线播放| 天堂中文最新版在线下载| 91国产中文字幕| 精品国产一区二区久久| 亚洲综合色网址| 亚洲精品久久午夜乱码| 又黄又粗又硬又大视频| 国产精品久久久久久精品古装| 丝袜喷水一区| 亚洲av美国av| 国产精品一区二区在线观看99| 美女大奶头黄色视频| 日本vs欧美在线观看视频| 人人澡人人妻人| 亚洲专区字幕在线| 激情视频va一区二区三区| 波多野结衣av一区二区av| 99精品欧美一区二区三区四区| 99国产综合亚洲精品| 日韩大码丰满熟妇| 久久久久久亚洲精品国产蜜桃av| 久久影院123| 欧美国产精品va在线观看不卡| 国产高清视频在线播放一区 | 国产男人的电影天堂91| 国产免费视频播放在线视频| 亚洲三区欧美一区| 青春草视频在线免费观看| 国产精品九九99| cao死你这个sao货| 国产成人一区二区三区免费视频网站| 免费久久久久久久精品成人欧美视频| 看免费av毛片| 在线 av 中文字幕| 人人妻,人人澡人人爽秒播| 1024香蕉在线观看| 精品免费久久久久久久清纯 | 人妻 亚洲 视频| 91成年电影在线观看| 亚洲av欧美aⅴ国产| 日本vs欧美在线观看视频| 啦啦啦中文免费视频观看日本| 午夜激情av网站| 午夜福利视频在线观看免费| 国产精品二区激情视频| 精品免费久久久久久久清纯 | 久久亚洲精品不卡| 久久中文字幕一级| 国产国语露脸激情在线看| 搡老乐熟女国产| 国产精品久久久久成人av| 99久久人妻综合| 男女床上黄色一级片免费看| 人妻一区二区av| 男女边摸边吃奶| 欧美一级毛片孕妇| 精品少妇内射三级| 动漫黄色视频在线观看| 亚洲国产欧美网| 999精品在线视频| 成年动漫av网址| 久久人妻福利社区极品人妻图片| 美女主播在线视频| 欧美大码av| 秋霞在线观看毛片| 精品一品国产午夜福利视频| 欧美日韩亚洲综合一区二区三区_| 别揉我奶头~嗯~啊~动态视频 | 亚洲国产看品久久| 黄色视频,在线免费观看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产高清videossex| 亚洲午夜精品一区,二区,三区| 午夜福利视频在线观看免费| 大香蕉久久网| 视频在线观看一区二区三区| 亚洲国产欧美在线一区| 亚洲成人手机| 蜜桃国产av成人99| 国产精品久久久人人做人人爽| 黄频高清免费视频| 亚洲国产av影院在线观看| 久久综合国产亚洲精品| 欧美 亚洲 国产 日韩一| 日韩制服丝袜自拍偷拍| 国产一区二区三区综合在线观看| 亚洲精品久久久久久婷婷小说| 777米奇影视久久| 国产免费一区二区三区四区乱码| 亚洲精品国产av蜜桃| 美女高潮喷水抽搐中文字幕| 电影成人av| 可以免费在线观看a视频的电影网站| 中文字幕av电影在线播放| 蜜桃国产av成人99| 久久毛片免费看一区二区三区| 人成视频在线观看免费观看| 国产成人欧美在线观看 | 人人妻人人澡人人爽人人夜夜| 欧美亚洲 丝袜 人妻 在线| 精品欧美一区二区三区在线| 欧美日韩国产mv在线观看视频| 欧美亚洲日本最大视频资源| 少妇被粗大的猛进出69影院| 香蕉丝袜av| 亚洲国产av影院在线观看| 久久青草综合色| 麻豆乱淫一区二区| 青春草视频在线免费观看| 免费观看人在逋| 少妇人妻久久综合中文| 两个人看的免费小视频| 狂野欧美激情性bbbbbb| 精品人妻1区二区| 亚洲国产毛片av蜜桃av| 国产在线免费精品| 国产在线视频一区二区| a级毛片在线看网站| 熟女少妇亚洲综合色aaa.| 热99国产精品久久久久久7| 97在线人人人人妻| 亚洲精品国产精品久久久不卡| 久久精品亚洲av国产电影网| 大片电影免费在线观看免费| 黑人操中国人逼视频| 成年av动漫网址| 国产在线视频一区二区| 欧美精品一区二区免费开放| 91大片在线观看| 秋霞在线观看毛片| 亚洲精品久久成人aⅴ小说| 国产免费福利视频在线观看| 日本vs欧美在线观看视频| 中文字幕av电影在线播放| 另类精品久久| 精品少妇黑人巨大在线播放| 制服诱惑二区| 亚洲 欧美一区二区三区| 婷婷成人精品国产| 99re6热这里在线精品视频| 成人手机av| 亚洲情色 制服丝袜| 不卡一级毛片| 免费在线观看完整版高清| 成人手机av| 淫妇啪啪啪对白视频 | 一区二区日韩欧美中文字幕| 亚洲国产欧美一区二区综合| 欧美另类一区| 欧美日韩av久久| 美女午夜性视频免费| 99国产精品一区二区蜜桃av | 欧美黑人精品巨大| 欧美精品人与动牲交sv欧美| √禁漫天堂资源中文www| 久久国产精品影院| 男女床上黄色一级片免费看| 中文欧美无线码| 免费在线观看日本一区| 国产亚洲一区二区精品| 精品乱码久久久久久99久播| 免费看十八禁软件| 一级毛片电影观看| 97人妻天天添夜夜摸| 久久亚洲精品不卡| 最新的欧美精品一区二区| 久久ye,这里只有精品| 天堂俺去俺来也www色官网| 国产91精品成人一区二区三区 | 中国美女看黄片| 人人澡人人妻人| 国产高清视频在线播放一区 | 青青草视频在线视频观看| 黄色视频在线播放观看不卡| 一区二区日韩欧美中文字幕| 久久久久久久久久久久大奶| 国产人伦9x9x在线观看| 9191精品国产免费久久| 午夜福利视频在线观看免费| 欧美午夜高清在线| 日韩熟女老妇一区二区性免费视频| 亚洲黑人精品在线| 国产麻豆69| 超碰成人久久| 成年动漫av网址| 日韩大码丰满熟妇| 色视频在线一区二区三区| 999久久久精品免费观看国产| 十八禁网站免费在线| 操出白浆在线播放| 精品国产国语对白av| 久久精品国产a三级三级三级| 人人妻人人澡人人爽人人夜夜| 久久久精品国产亚洲av高清涩受| 国产在视频线精品| 新久久久久国产一级毛片| 天堂中文最新版在线下载| 大码成人一级视频| 曰老女人黄片| 国产又色又爽无遮挡免| 嫁个100分男人电影在线观看| 久久久久久久久免费视频了| 久久久久国产精品人妻一区二区| 精品人妻熟女毛片av久久网站| 国产99久久九九免费精品| 国产日韩欧美视频二区| 热99久久久久精品小说推荐| 免费人妻精品一区二区三区视频| 精品一区在线观看国产| 蜜桃国产av成人99| 久久久水蜜桃国产精品网| 日本一区二区免费在线视频| 亚洲av成人一区二区三| 新久久久久国产一级毛片| 精品免费久久久久久久清纯 | 久久精品国产亚洲av香蕉五月 | 下体分泌物呈黄色| 久久久精品免费免费高清| 欧美日韩视频精品一区| 久久国产亚洲av麻豆专区| 国产xxxxx性猛交| 丝袜人妻中文字幕| 中文字幕高清在线视频| 国产在线一区二区三区精| 天天操日日干夜夜撸| 精品第一国产精品| 三上悠亚av全集在线观看| 亚洲性夜色夜夜综合| 免费久久久久久久精品成人欧美视频| 黄网站色视频无遮挡免费观看| 免费观看a级毛片全部| 91大片在线观看| 久久精品国产a三级三级三级| 99国产精品一区二区蜜桃av | www.999成人在线观看| 十分钟在线观看高清视频www| 老司机在亚洲福利影院| 亚洲男人天堂网一区| 性少妇av在线| 水蜜桃什么品种好| 男女床上黄色一级片免费看| 精品国产一区二区久久| 美女国产高潮福利片在线看| 久久毛片免费看一区二区三区| 美女主播在线视频| 国产精品久久久久成人av| 国产一卡二卡三卡精品| 热99re8久久精品国产| 亚洲午夜精品一区,二区,三区| 久久久久精品人妻al黑| 亚洲成国产人片在线观看| 久久久久国产精品人妻一区二区| 欧美激情极品国产一区二区三区| av片东京热男人的天堂| 少妇 在线观看| 在线天堂中文资源库| 各种免费的搞黄视频| 欧美成人午夜精品| 女人久久www免费人成看片| 国产亚洲欧美精品永久| 12—13女人毛片做爰片一| 久久中文字幕一级| www.熟女人妻精品国产| 午夜精品国产一区二区电影| 午夜福利在线观看吧| 亚洲av片天天在线观看| 亚洲av美国av| 亚洲综合色网址| 黄色 视频免费看| 亚洲七黄色美女视频| 精品少妇久久久久久888优播| 免费黄频网站在线观看国产| 在线观看免费日韩欧美大片| 一区二区三区乱码不卡18| 国产91精品成人一区二区三区 | 女人精品久久久久毛片| 国产日韩欧美视频二区| 黄色视频,在线免费观看| 男女国产视频网站| 乱人伦中国视频| 国产黄色免费在线视频| 大片免费播放器 马上看| 国产又色又爽无遮挡免| bbb黄色大片| 777久久人妻少妇嫩草av网站| 五月开心婷婷网| 国产片内射在线| 乱人伦中国视频| 伊人亚洲综合成人网| 热re99久久国产66热| 三上悠亚av全集在线观看| 岛国毛片在线播放| svipshipincom国产片| 亚洲av成人不卡在线观看播放网 | a级毛片黄视频| 久久久水蜜桃国产精品网| 国产亚洲欧美精品永久| 欧美日韩亚洲高清精品| 亚洲一码二码三码区别大吗| 久久国产亚洲av麻豆专区| 欧美老熟妇乱子伦牲交| 99久久精品国产亚洲精品| 久久人人爽人人片av| 最近中文字幕2019免费版| 精品少妇久久久久久888优播| 一区福利在线观看| 精品国产一区二区久久| 亚洲熟女毛片儿| 天堂中文最新版在线下载| 亚洲熟女精品中文字幕| 免费日韩欧美在线观看| 一二三四社区在线视频社区8| 午夜激情久久久久久久| 嫩草影视91久久| 久久性视频一级片| 色老头精品视频在线观看| 色94色欧美一区二区| 国产精品一区二区免费欧美 | 极品人妻少妇av视频| 国产亚洲一区二区精品| 男女下面插进去视频免费观看| 亚洲av片天天在线观看| 欧美 日韩 精品 国产| 男女床上黄色一级片免费看| www.熟女人妻精品国产| 午夜免费成人在线视频| 别揉我奶头~嗯~啊~动态视频 | 天堂中文最新版在线下载| 日日摸夜夜添夜夜添小说| 久久天堂一区二区三区四区| 在线av久久热| 91字幕亚洲| 91国产中文字幕| 色老头精品视频在线观看| 亚洲性夜色夜夜综合| 日日夜夜操网爽| a级片在线免费高清观看视频| 在线av久久热| 纯流量卡能插随身wifi吗| 免费人妻精品一区二区三区视频| 亚洲午夜精品一区,二区,三区| 免费观看av网站的网址| 男男h啪啪无遮挡| 99精品欧美一区二区三区四区| 免费av中文字幕在线| 天天影视国产精品| 国产国语露脸激情在线看| 妹子高潮喷水视频| 十八禁人妻一区二区| 色婷婷久久久亚洲欧美| 久久精品人人爽人人爽视色| 亚洲欧美日韩高清在线视频 | 1024香蕉在线观看| 中文精品一卡2卡3卡4更新| 精品亚洲成a人片在线观看| 别揉我奶头~嗯~啊~动态视频 | 1024视频免费在线观看| 高清欧美精品videossex| 色94色欧美一区二区| 9色porny在线观看| 青春草视频在线免费观看| 在线av久久热| 精品亚洲成国产av| 亚洲国产av影院在线观看| 亚洲国产欧美日韩在线播放| 超色免费av| 亚洲欧美日韩另类电影网站| 欧美精品一区二区大全| 亚洲av成人一区二区三| 国产精品九九99| 黄色a级毛片大全视频| 国产不卡av网站在线观看| 国产男女内射视频| 日本91视频免费播放| 亚洲国产av新网站| 日韩制服丝袜自拍偷拍| 男人操女人黄网站| 国产视频一区二区在线看| 丝袜美腿诱惑在线| 国产亚洲欧美精品永久| 激情视频va一区二区三区| 亚洲第一av免费看| 欧美日本中文国产一区发布| a级毛片在线看网站| 精品国产一区二区久久| 国产精品欧美亚洲77777| av电影中文网址| 十八禁高潮呻吟视频| 欧美精品人与动牲交sv欧美| 精品卡一卡二卡四卡免费| 亚洲男人天堂网一区| 丰满迷人的少妇在线观看| 国产91精品成人一区二区三区 | 国产成人av激情在线播放| 人人妻人人爽人人添夜夜欢视频| 制服人妻中文乱码| 国产欧美日韩一区二区三 | 亚洲专区国产一区二区| 最近最新免费中文字幕在线| a级毛片黄视频| 亚洲av日韩在线播放| 老司机影院成人| 亚洲一码二码三码区别大吗| 欧美精品啪啪一区二区三区 | 亚洲国产中文字幕在线视频| 美女中出高潮动态图| 午夜激情久久久久久久| 国产一区有黄有色的免费视频| 夜夜夜夜夜久久久久| 亚洲av成人不卡在线观看播放网 | 欧美 日韩 精品 国产| 99精品久久久久人妻精品| 欧美 日韩 精品 国产| 欧美日韩福利视频一区二区| 少妇粗大呻吟视频| 亚洲人成电影免费在线| 动漫黄色视频在线观看| a在线观看视频网站| 老司机影院成人| 永久免费av网站大全| 色婷婷久久久亚洲欧美| 亚洲一码二码三码区别大吗| 岛国在线观看网站| 69av精品久久久久久 | 97精品久久久久久久久久精品| 汤姆久久久久久久影院中文字幕| 亚洲五月色婷婷综合| 大香蕉久久成人网| 妹子高潮喷水视频| 久9热在线精品视频| 在线亚洲精品国产二区图片欧美| 操出白浆在线播放| 国产精品影院久久| 国产精品亚洲av一区麻豆| 少妇精品久久久久久久| av网站在线播放免费| 精品人妻在线不人妻| 成人影院久久| 亚洲av美国av| 亚洲色图综合在线观看| 色94色欧美一区二区| 精品久久久精品久久久| 亚洲欧美精品自产自拍| 9热在线视频观看99| 最新的欧美精品一区二区| 另类精品久久| 最黄视频免费看| 精品人妻1区二区| 中文字幕色久视频| 欧美黑人精品巨大| 欧美精品av麻豆av| 日本黄色日本黄色录像| 久久国产亚洲av麻豆专区| 91精品三级在线观看| 看免费av毛片| 国产一区二区三区综合在线观看| 黄色毛片三级朝国网站| 最近中文字幕2019免费版| 亚洲专区字幕在线| 母亲3免费完整高清在线观看| 亚洲专区国产一区二区| 成年人黄色毛片网站| 欧美精品亚洲一区二区| 久久久久久久久久久久大奶| 日韩一区二区三区影片| 欧美人与性动交α欧美精品济南到| 男人舔女人的私密视频| 最新在线观看一区二区三区| 精品一区二区三卡| 成人影院久久| 亚洲成av片中文字幕在线观看| 少妇猛男粗大的猛烈进出视频| 久久人妻福利社区极品人妻图片| 美女国产高潮福利片在线看| 亚洲精品日韩在线中文字幕| 热re99久久精品国产66热6| 18禁黄网站禁片午夜丰满| 精品人妻一区二区三区麻豆| 爱豆传媒免费全集在线观看| 一区二区日韩欧美中文字幕| 9色porny在线观看| 91成年电影在线观看| www日本在线高清视频| 国产黄色免费在线视频| 欧美 日韩 精品 国产| 十分钟在线观看高清视频www| 欧美黑人精品巨大| 我的亚洲天堂| 大陆偷拍与自拍| 精品福利观看| 国产一区二区三区av在线| 国产精品影院久久|