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

    微咸水微潤(rùn)灌溉下土壤水鹽運(yùn)移特性研究

    2017-02-08 01:26:52張珂萌牛文全汪有科薛萬(wàn)來(lái)張子卓
    關(guān)鍵詞:微潤(rùn)咸水礦化度

    張珂萌 牛文全,2 汪有科,2 薛萬(wàn)來(lái) 張子卓

    (1.西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院, 陜西楊凌 712100; 2.西北農(nóng)林科技大學(xué)水土保持研究所, 陜西楊凌 712100;3.北京市水科學(xué)技術(shù)研究院, 北京 100048)

    微咸水微潤(rùn)灌溉下土壤水鹽運(yùn)移特性研究

    張珂萌1牛文全1,2汪有科1,2薛萬(wàn)來(lái)2,3張子卓1

    (1.西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院, 陜西楊凌 712100; 2.西北農(nóng)林科技大學(xué)水土保持研究所, 陜西楊凌 712100;3.北京市水科學(xué)技術(shù)研究院, 北京 100048)

    為探明土壤水分和鹽分在利用微咸水進(jìn)行微潤(rùn)灌溉條件下的運(yùn)移情況,采用室內(nèi)土箱模擬試驗(yàn)方式,設(shè)置2.0、2.5、3.0、3.5、5.0 g/L 5種不同礦化度處理,以蒸餾水處理作對(duì)照,共入滲72 h。結(jié)果表明:入滲結(jié)束時(shí)在不同方向上的最大運(yùn)移距離隨礦化度增大呈先增大后減小趨勢(shì),在3.0 g/L處理下達(dá)到最大值,且微咸水處理的濕潤(rùn)鋒運(yùn)移距離均大于蒸餾水處理;將累積入滲量代入Kostiakov入滲公式,入滲系數(shù)隨礦化度的增大呈先增大后減小趨勢(shì),入滲指數(shù)不斷減??;土壤電導(dǎo)率以微潤(rùn)帶為軸心向四周不斷增大,在濕潤(rùn)鋒處達(dá)到最大值,脫鹽區(qū)與濕潤(rùn)體形狀相關(guān),呈圓環(huán)狀分布;入滲結(jié)束后土壤剖面平均含鹽量與蒸餾水處理之間無(wú)顯著性差異,脫鹽半徑隨礦化度的增大呈線性遞減趨勢(shì);利用微潤(rùn)灌進(jìn)行灌溉,土壤鹽分存在表聚和底聚現(xiàn)象,且表層積鹽更為嚴(yán)重。

    礦化度; 入滲; 水鹽分布; 微潤(rùn)灌溉

    引言

    我國(guó)人均淡水資源不足且污染嚴(yán)重,農(nóng)業(yè)用水利用效率低下[1]致使灌溉水資源日益短缺,而微咸水分布廣、儲(chǔ)量大[2],已成為重要的灌溉水源[3],在淡水十分緊缺的華北和西北地區(qū)有著廣闊的開(kāi)發(fā)利用前景,如何合理利用微咸水進(jìn)行農(nóng)業(yè)灌溉對(duì)于緩解我國(guó)淡水資源緊缺具有重要意義。近年來(lái),部分學(xué)者利用微咸水針對(duì)不同灌溉方式[4-8]、土壤質(zhì)地以及作物[9-13]開(kāi)展了大量研究,普遍認(rèn)為:①微咸水灌溉對(duì)土壤結(jié)構(gòu)的影響主要表現(xiàn)在土壤交換性鈉和土壤溶液電導(dǎo)率2個(gè)方面,提高土壤鹽分濃度有利于促進(jìn)土壤顆粒的紊凝作用和團(tuán)聚性、加強(qiáng)土壤穩(wěn)定性、增加土壤中大孔隙、增強(qiáng)滲透性能;同時(shí),Na+過(guò)量會(huì)造成土壤顆粒松散,降低土壤滲透能力[14-18],利用微咸水灌溉存在礦化度臨界值。②利用微咸水灌溉產(chǎn)生輕微的鹽分脅迫可以在一定程度上提高小麥[9]、棉花[10-11]、甜瓜[12]等作物的產(chǎn)量和品質(zhì),但對(duì)番茄[13]沒(méi)有顯著影響,不同作物對(duì)鹽分的耐受性呈現(xiàn)差異,應(yīng)將植物根區(qū)的鹽分控制在植物對(duì)鹽分的耐受范圍內(nèi),才不至造成減產(chǎn)。③滴灌、涌泉灌溉等間歇灌溉對(duì)土壤鹽分淋洗效果較漫灌、溝灌等一次灌溉淋洗效果明顯,鹽分在濕潤(rùn)鋒處累積[19],脫鹽區(qū)的形狀根據(jù)不同灌溉方式下濕潤(rùn)區(qū)形狀呈現(xiàn)差異,如漫灌呈淺寬式、滴灌為寬窄式[20];另有研究者對(duì)微咸水進(jìn)行磁化處理,發(fā)現(xiàn)磁化處理后的微咸水可以有效增益鹽敏植物對(duì)鹽漬化土壤環(huán)境的適應(yīng)性[21]。微潤(rùn)灌溉屬于連續(xù)灌溉,通過(guò)地埋微潤(rùn)管,減小地表蒸發(fā)并將水緩慢輸送到作物根區(qū)土壤[18],形成近似圓柱體形[22-23]濕潤(rùn)體以滿足植物需水。然而目前有關(guān)利用微咸水進(jìn)行微潤(rùn)灌溉時(shí)土壤中水-鹽運(yùn)移過(guò)程、脫鹽區(qū)域的大小、鹽分累積位置等研究尚不足。本文通過(guò)室內(nèi)土箱模擬試驗(yàn),分析不同礦化度微咸水微潤(rùn)灌溉下土壤水-鹽運(yùn)移特性、濕潤(rùn)土壤脫鹽區(qū)域等,以期為微咸水微潤(rùn)灌溉提供一定的理論依據(jù)。

    1 材料與方法

    1.1 試驗(yàn)土壤

    試驗(yàn)用土取自陜西楊凌渭河三級(jí)階地農(nóng)田塿土。取土深度為0~40 cm,將取得的土壤經(jīng)風(fēng)干、碾壓、過(guò)篩(孔徑2 mm)、均勻混合制成室內(nèi)試驗(yàn)土樣。試驗(yàn)土壤的顆粒組成采用MS2000型激光分析粒度儀(英國(guó)馬爾文公司)測(cè)定,按國(guó)際制土壤質(zhì)地三角形分類,粉粒(粒徑0.002~0.02 mm)質(zhì)量分?jǐn)?shù)在50%以上,粘粒(粒徑0~0.002 mm)質(zhì)量分?jǐn)?shù)低于15%,試驗(yàn)所用土壤類型為粉質(zhì)粘壤土。供試土壤的具體理化性質(zhì)見(jiàn)表1。

    表1 供試土樣的理化性質(zhì)
    Tab.1 Physical and chemical properties of soil

    1.2 試驗(yàn)水質(zhì)

    研究灌溉水水質(zhì)對(duì)土壤入滲的影響主要有灌溉水礦化度和鈉吸附比(Sodium adsorption ratio,SAR)2個(gè)指標(biāo),NIKOS等[24]發(fā)現(xiàn)礦化度小于等于3 g/L條件下,可以認(rèn)為入滲水的SAR對(duì)土壤鹽分分布無(wú)顯著影響,肖振華等[25]通過(guò)試驗(yàn)初步提出了引起堿害的SAR臨界值為14(mmol/L)1/2,為避免堿害發(fā)生,本次試驗(yàn)用水的SAR采用14(mmol/L)1/2。礦化水根據(jù)蒲城鹵陽(yáng)湖地區(qū)地下微咸水鹽分組成特點(diǎn),由NaHCO3、Na2SO4、MgSO4、MgCl2、CaCl2配制而成,各礦化度下基本離子濃度如表2所示。

    1.3 試驗(yàn)裝置

    試驗(yàn)裝置由土箱和供水裝置2部分組成。試驗(yàn)土箱尺寸(長(zhǎng)×寬×高)為120 cm×50 cm×50 cm,由厚度為10 mm有機(jī)玻璃制成。土箱兩側(cè)打有直徑為25 mm的對(duì)稱小孔(距土箱上沿250 mm側(cè)沿15 mm)為微潤(rùn)帶穿過(guò)該土箱供水的進(jìn)出口。供水裝置為馬氏瓶,提供恒定水頭,其橫截面積25 cm2,高70 cm。微潤(rùn)帶為雙層結(jié)構(gòu),直徑25 mm,內(nèi)層厚度0.06 mm,微孔直徑10~900 nm不等,外層為無(wú)紡布保護(hù)層以增加微潤(rùn)帶強(qiáng)度。微潤(rùn)帶長(zhǎng)度與土箱長(zhǎng)度均為120 cm,進(jìn)口端連接供水系統(tǒng),出口端封閉,試驗(yàn)土箱與馬氏瓶利用橡皮軟管連接,土箱表面塑封,如圖1所示。

    表2 試驗(yàn)配置不同礦化度微咸水的基本離子濃度Tab.2 Ion content of water with different mineralization degrees mol/L

    圖1 試驗(yàn)裝置示意圖Fig.1 Schematic of experimental equipment1.馬氏瓶 2.輸水軟管 3.可移動(dòng)支架 4.微潤(rùn)帶 5.試驗(yàn)土箱

    1.4 試驗(yàn)方法及測(cè)定內(nèi)容

    試驗(yàn)共設(shè)定5種礦化度,分別為2.0、2.5、3.0、3.5、5.0 g/L作微咸水處理,以蒸餾水(0 g/L)處理作為對(duì)照。本次試驗(yàn)設(shè)置微潤(rùn)帶埋深為15 cm[22],壓力水頭設(shè)置為180 cm。將試驗(yàn)土樣(容重1.35 g/cm3,風(fēng)干含水率2.9%)裝入土箱,層間打毛,以便于充分接觸,表面用塑料膜覆蓋,防止土壤水分蒸發(fā)。試驗(yàn)過(guò)程中記錄不同時(shí)刻馬氏瓶的水位線。微潤(rùn)帶充水埋入土后,立即開(kāi)始用秒表記錄灌水時(shí)間,按照先密后疏的原則,入滲72 h后停止供水,用1 cm土鉆取樣,取樣點(diǎn)位置沿管帶的水平和垂直方向同時(shí)進(jìn)行,間距為5 cm。

    用干燥法測(cè)定土壤含水率,將測(cè)定完含水率的土壤按照土水質(zhì)量比1∶5進(jìn)行浸提,使用哈納多參數(shù)分析儀Hi-4522型(意大利)測(cè)定土壤浸提液電導(dǎo)率。試驗(yàn)結(jié)束后,累積入滲量根據(jù)馬氏瓶刻度和橫截面積計(jì)算,平均入滲率為單位時(shí)間內(nèi)累積入滲量。

    各處理重復(fù)3次,分析數(shù)據(jù)時(shí)取平均值,采用SPSS 17.0對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行雙尾-t檢驗(yàn)和Duncan多重比較分析,用Excel進(jìn)行圖形繪制。

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

    2.1 微咸水對(duì)土壤入滲特征的影響

    2.1.1 微咸水對(duì)土壤濕潤(rùn)鋒運(yùn)移距離的影響

    圖2是連續(xù)入滲0.5、30、72 h時(shí),土壤濕潤(rùn)鋒在水平、垂直向上和垂直向下3個(gè)方向的運(yùn)移距離隨入滲水礦化度增加的變化趨勢(shì)。圖2a表明灌水初期(0.5 h)不同礦化度下的微咸水在水平方向的運(yùn)移距離大于其他2個(gè)方向,濕潤(rùn)鋒運(yùn)移距離在3個(gè)方向上基本隨入滲水礦化度的增加而呈先增大后減小的趨勢(shì);圖2b、2c表明灌水中期(30 h)和灌水后期(72 h)礦化度為3.0 g/L處理的濕潤(rùn)鋒運(yùn)移距離在不同運(yùn)移方向上均大于2.5 g/L和3.5 g/L,達(dá)到最大,說(shuō)明在相同的土壤質(zhì)地、容重、初始含水率和入滲時(shí)間下利用微咸水進(jìn)行微潤(rùn)灌溉時(shí)濕潤(rùn)鋒運(yùn)移距離的最大值處于3.0 g/L附近,本次試驗(yàn)中3.0 g/L處理的濕潤(rùn)鋒運(yùn)移距離在試驗(yàn)結(jié)束時(shí)為最大值;入滲72 h末3.0 g/L處理在垂直向下方向的濕潤(rùn)鋒運(yùn)移距離較蒸餾水處理隨灌水增加的百分比分別為50%、30%和25%,微咸水在連續(xù)灌溉下的濕潤(rùn)鋒運(yùn)移距離的增加程度逐漸降低。

    圖2 入滲0.5、30、72 h時(shí)濕潤(rùn)鋒運(yùn)移距離與礦化度之間的變化曲線Fig.2 Changing curves of maximum wetting front with different mineralization degrees of saline water in different time

    表3為入滲72 h后不同礦化度濕潤(rùn)鋒運(yùn)移距離。垂直向下、向上和水平方向上不同礦化度濕潤(rùn)鋒運(yùn)移距離的標(biāo)準(zhǔn)差分別為1.50、1.51和1.07,對(duì)樣本均值雙尾-t檢驗(yàn),P≤0.05,垂直向下和向上方向的決定系數(shù)R2為0.825,垂直向上和水平方向的決定系數(shù)R2為0.946,垂直向下和水平方向的決定系數(shù)R2為0.875,表明當(dāng)?shù)V化度相同時(shí),濕潤(rùn)鋒在不同方向上的運(yùn)移距離差異較小,礦化度對(duì)不同方向運(yùn)移距離的影響不顯著,利用微潤(rùn)灌進(jìn)行微咸水灌溉的濕潤(rùn)體近似圓形;不同礦化度對(duì)濕潤(rùn)鋒運(yùn)移距離影響的顯著性分析結(jié)果表明0、3.5、5.0 g/L之間的差異不顯著,2.0、2.5、3.5 g/L之間的差異不顯著,但3.0 g/L礦化度條件下濕潤(rùn)鋒運(yùn)移距離顯著高于其他處理,并在表層出現(xiàn)一定的積水,積水深小于1 cm。

    表3 入滲72 h濕潤(rùn)鋒運(yùn)移距離Tab.3 Transport distance of wetting front under different salinities in different directions after 72 h cm

    注:同列數(shù)值后不同字母表示差異顯著,P<0.05,下同。

    2.1.2 微咸水對(duì)土壤入滲參數(shù)的影響

    土壤入滲參數(shù)是反映入滲條件下入滲速率的重要參數(shù),受土壤質(zhì)地、結(jié)構(gòu)、容重等因素的影響。對(duì)于同一土壤,入滲水礦化度也影響土壤的入滲性能,進(jìn)而影響其入滲參數(shù)。

    Kostiakov經(jīng)驗(yàn)入滲公式能很好地描述微咸水的土壤入滲特征[26-27],即土壤累積入滲量與入滲時(shí)間的關(guān)系,擬合決定系數(shù)R2在0.99 以上。Kostiakov經(jīng)驗(yàn)入滲公式為

    I=kta

    (1)

    式中I——累積入滲量,cmt——入滲時(shí)間,hk——入滲系數(shù)a——入滲指數(shù)

    由表4可以得到,入滲系數(shù)k隨礦化度的增加呈先增加后減小趨勢(shì),并在3.0 g/L處達(dá)到最大值,且微咸水處理的k值均大于蒸餾水處理,灌溉水的礦化度對(duì)入滲系數(shù)k具有顯著影響。入滲指數(shù)a隨著礦化度的增大呈先減小后在5.0 g/L時(shí)與蒸餾水處理相近。利用微咸水灌溉時(shí),存在一個(gè)礦化度的臨界值,隨礦化度的增大,k值逐漸增大,a值逐漸減小,但超過(guò)該臨界值時(shí),隨著礦化度的繼續(xù)增大k值呈遞減趨勢(shì),本次試驗(yàn)驗(yàn)證該值在微潤(rùn)灌條件下位于3.0 g/L附近,表中a值在5.0 g/L時(shí)增大,說(shuō)明當(dāng)微咸水礦化度增加到5.0 g/L時(shí)不再促進(jìn)微潤(rùn)灌溉條件下微咸水的入滲。

    表4 礦化度對(duì)Kostiakov 入滲參數(shù)的影響
    Tab.4 Effect of mineralization of saline water on infiltration parameters of Kostiakov model

    礦化度/(g·L-1)kaR202.086a0.850a0.9992.03.162b0.801b0.9972.54.162c0.791c0.9993.05.546d0.775d0.9983.53.391e0.773d0.9945.02.748f0.834e0.999

    2.2 微咸水入滲對(duì)土壤脫鹽區(qū)的影響

    入滲72 h后,沿微潤(rùn)管垂直方向上取樣測(cè)土壤電導(dǎo)率,當(dāng)不同深度實(shí)測(cè)電導(dǎo)率低于初始電導(dǎo)率時(shí)稱為脫鹽,脫鹽區(qū)域邊沿與微潤(rùn)管的距離為脫鹽半徑。圖3為灌溉水礦化度與脫鹽半徑的關(guān)系。

    圖3 礦化度與脫鹽半徑的關(guān)系Fig.3 Relationship between salinity and radius of desalination in soil

    從圖3可以看出,微潤(rùn)灌脫鹽半徑隨礦化度的增加呈線性遞減趨勢(shì),蒸餾水處理的脫鹽半徑最大,但礦化度2.5 g/L與3.0 g/L處理的差異較小,而礦化度3.0 g/L與3.5 g/L的差異較大,分析其原因是處理3.0 g/L和處理3.5 g/L之間含鹽量差異并不大,但3.0 g/L的入滲量遠(yuǎn)大于3.5 g/L,鹽隨水運(yùn)移而造成3.0 g/L處理的脫鹽半徑大于3.5 g/L。從脫鹽半徑的變化趨勢(shì)分析,礦化度為3.0 g/L是一個(gè)拐點(diǎn),當(dāng)?shù)V化度小于3.0 g/L時(shí),濕潤(rùn)體脫鹽占優(yōu),反之則積鹽占優(yōu)。

    2.3 微咸水入滲對(duì)土壤剖面鹽分積累的影響

    土壤電導(dǎo)率與土壤含鹽量之間存在正相關(guān)關(guān)系,可用電導(dǎo)率反映土壤含鹽量的變化情況[26]。圖4是距離微潤(rùn)管帶水平方向5 cm處,垂直方向上不同礦化度條件下入滲72 h時(shí)土壤電導(dǎo)率的變化。

    圖4 不同礦化度下土壤電導(dǎo)率隨土壤深度的變化Fig.4 Changes of soil EC value with soil depth for different mineralization degrees of saline water

    由圖4可知,微潤(rùn)灌土壤電導(dǎo)率分布以微潤(rùn)帶(埋深15 cm處)為軸心向四周逐漸增大,在濕潤(rùn)鋒處達(dá)到最大值。微潤(rùn)帶以上部分(0~15 cm),相同深度的濕潤(rùn)土壤電導(dǎo)率隨著灌溉水礦化度的增大而增大;微潤(rùn)帶以下部分(15~30 cm),總體與上部相似,相同深度濕潤(rùn)土壤電導(dǎo)率隨著灌溉水礦化度的增大而增大,但3.0 g/L處理下的電導(dǎo)率變化曲線與2.5 g/L處理出現(xiàn)交叉,3.0 g/L處理的脫鹽區(qū)域略大于2.5 g/L處理,在2.5~3.0 g/L之間即使增加灌溉水的礦化度也不會(huì)造成進(jìn)一步積鹽。不同礦化度(0~5.0 g/L)濕潤(rùn)鋒垂直向下與表層運(yùn)移距離分別為26.8、28.82、29.93、31.94、28.47、27.65 cm,對(duì)照?qǐng)D4垂直向下土壤電導(dǎo)率隨礦化度變化情況發(fā)現(xiàn)濕潤(rùn)土壤中的電導(dǎo)率隨微咸水的入滲向濕潤(rùn)鋒附近增大,即鹽分隨入滲向濕潤(rùn)鋒處累積,垂直向上的變化情況與之一致,微咸水微潤(rùn)灌鹽分分布與濕潤(rùn)體形狀。圖中分界線267 μS/cm為土壤初始電導(dǎo)率,不同深度濕潤(rùn)土壤的電導(dǎo)率大于土壤初始電導(dǎo)率視為積鹽,隨灌溉水礦化度的增加表層積鹽深度也隨之增加并分布在距表層11.9 cm內(nèi),底部同時(shí)存在積鹽且蒸餾水處理較其他微咸水處理輕微,微潤(rùn)灌溉鹽分隨水上移和底層遷移同時(shí)存在,但淺層土壤更易積鹽。

    由圖5試驗(yàn)結(jié)束時(shí)不同礦化度下水平距離微潤(rùn)帶5 cm處土壤剖面平均電導(dǎo)率分布,可進(jìn)一步分析礦化度對(duì)土壤剖面電導(dǎo)率的影響。由圖5可知,蒸餾水處理的平均電導(dǎo)率低于土壤初始電導(dǎo)率,具有一定的洗鹽作用,而利用微咸水灌溉則不同程度地增加了土壤的含鹽量。單因素方差分析得到不同礦化度條件下土壤電導(dǎo)率變化不顯著(P=0.149),在礦化度小于等于5.0 g/L且距表層25 cm內(nèi)土壤中的平均電導(dǎo)率受礦化度影響有升高的趨勢(shì),但之間的差異不顯著,說(shuō)明利用微咸水進(jìn)行灌溉礦化度在2.0~5.0 g/L以內(nèi)對(duì)土壤總含鹽量的影響差異不大,但在微潤(rùn)灌溉入滲的方式下不同深度土層鹽分分布存在一定的差異(表5)。

    表5為不同深度土層5個(gè)礦化度的平均電導(dǎo)率??煽闯霾煌疃韧翆又g的電導(dǎo)率有顯著性差異(P≤0.05),5 cm處平均電導(dǎo)率最大,與其他各土層差異顯著,其次是25 cm處,而15 cm和20 cm差異不顯著,20 cm和10 cm差異不顯著。隨微咸水的入滲,鹽分隨水分以埋管15 cm處為中心呈對(duì)稱分布,在5 cm和25 cm處積鹽較為嚴(yán)重,10~20 cm內(nèi)的電導(dǎo)率均小于初始電導(dǎo)率(表1,267 μS/cm),說(shuō)明位于濕潤(rùn)體內(nèi)鹽分降低,得到淋洗;鹽分在濕潤(rùn)鋒處積累最為嚴(yán)重,產(chǎn)生鹽分表層和深層雙聚的現(xiàn)象,且鹽分的表聚程度大于底層聚集情況。

    圖5 礦化度對(duì)土壤剖面平均電導(dǎo)率的影響Fig.5 Effects of different mineralization degrees on averaged EC value in soil profile

    土層深度/cm510152025平均電導(dǎo)率/(μS·cm-1)335.00a259.00c222.08d244.44cd298.90b

    3 討論

    (1)微潤(rùn)灌條件下,入滲水礦化度對(duì)入滲及濕潤(rùn)鋒運(yùn)移距離的影響

    濕潤(rùn)鋒運(yùn)移距離、入滲系數(shù)k通常受土壤質(zhì)地[28]、灌溉水水質(zhì)[29]、初始含水率[14]、供水壓力[26]等的影響;而入滲指數(shù)a主要與土壤的機(jī)械組成、結(jié)構(gòu)、供水強(qiáng)度和初始含水率等因素有關(guān),當(dāng)土壤結(jié)構(gòu)改變,孔隙度增大,導(dǎo)水性增強(qiáng)時(shí)濕潤(rùn)峰運(yùn)移距離和k值就會(huì)增大[30],土壤骨架發(fā)生一定的變形,使得入滲界面的土壤結(jié)構(gòu)密實(shí),入滲指數(shù)則減小[25]。微咸水的水鹽運(yùn)移特征主要是由入滲水的礦化度和土壤中Na+濃度共同決定的[25,31],在礦化度小于等于3.0 g/L時(shí),礦化度占主要作用,能夠增強(qiáng)土壤的凝絮作用[14],增加土壤的有效孔隙,改善導(dǎo)水性能,故k值和濕潤(rùn)鋒運(yùn)移距離隨礦化度的增大而增大;當(dāng)?shù)V化度繼續(xù)增加到大于3.0 g/L時(shí),Na+濃度由表2中的0.644 7 mol/L增加到1.074 5 mol/L,土壤結(jié)構(gòu)主要由Na+作用,使得土壤團(tuán)粒結(jié)構(gòu)分散膨脹,連通性孔隙變小和阻塞[26],降低導(dǎo)水性能,土壤更加密實(shí),此時(shí)k值和濕潤(rùn)鋒運(yùn)移距離會(huì)逐漸減小,a值減小。表4中蒸餾水處理和5.0 g/L處理a值相近,說(shuō)明當(dāng)?shù)V化度升高到5.0 g/L時(shí),土壤結(jié)構(gòu)發(fā)生較大改變[32],發(fā)生崩解等現(xiàn)象造成土壤松動(dòng)致使a值升高。

    (2)微潤(rùn)灌溉條件下,入滲水礦化度對(duì)濕潤(rùn)體含鹽量及脫鹽半徑的影響

    微咸水本身含有一定的鹽分,利用微咸水灌溉必然會(huì)給土壤帶來(lái)額外的鹽分積累,鹽分在土壤中的分布狀況對(duì)植物生長(zhǎng)有極重要的影響[33]。土壤鹽分的表聚現(xiàn)象在地表灌溉方式下主要受土壤表面蒸發(fā)和植物根系作用[10],而在微潤(rùn)灌溉方式下濕潤(rùn)土壤中鹽分隨水分入滲以微潤(rùn)帶為圓心向四周擴(kuò)散,與濕潤(rùn)體形狀相關(guān),呈圓環(huán)狀分布[34],故加劇了鹽分的表聚。

    4 結(jié)論

    (1)微潤(rùn)灌溉條件下,微咸水入滲結(jié)束時(shí)的土壤濕潤(rùn)鋒運(yùn)移距離、入滲系數(shù)k值隨礦化度的增大先增大再減小,入滲指數(shù)a值逐漸減小,灌溉水礦化度存在臨界值3.0 g/L使?jié)駶?rùn)鋒運(yùn)移距離和k值達(dá)到最大。

    (2)利用微潤(rùn)灌進(jìn)行灌溉時(shí),土壤鹽分分布存在向表層和土壤深層聚積現(xiàn)象,且表層更易積鹽。

    (3)微咸水微潤(rùn)灌溉土壤鹽分以微潤(rùn)帶為中心對(duì)稱分布且逐漸增大,并在濕潤(rùn)鋒處達(dá)到最大值;微咸水灌溉增加土壤的含鹽量,且同一土層深度土壤含鹽量隨礦化度的增大而增大,土壤平均含鹽量在不同礦化度處理之間無(wú)顯著性差異,微咸水入滲對(duì)不同土層鹽分分布影響顯著;脫鹽半徑與灌溉水礦化度呈線性負(fù)相關(guān)關(guān)系;脫鹽區(qū)圍繞微潤(rùn)帶分布,積鹽區(qū)在濕潤(rùn)鋒附近,濕潤(rùn)體外側(cè),呈圓環(huán)狀分布。

    1 韓洪云,趙連閣.中國(guó)灌溉農(nóng)業(yè)發(fā)展——問(wèn)題與挑戰(zhàn)[J].水利經(jīng)濟(jì),2004,22(1):54-58. HAN Hongyun, ZHAO Lian’ge. Development of science and technology on conservancy under the market economy[J]. Journal of Economics of Water Resources,2004,22(1):54-58. (in Chinese)

    2 張新民.干旱區(qū)水資源量與質(zhì)統(tǒng)一管理研究[D]. 西安:西安理工大學(xué),2000. ZHANG Xinmin. Water quantity and quality unified management research in arid areas[D]. Xi’an: Xi’an University of Technology, 2000. (in Chinese)

    3 王艷娜,侯振安,龔江,等.咸水資源農(nóng)業(yè)灌溉應(yīng)用研究進(jìn)展與展望[J].中國(guó)農(nóng)學(xué)通報(bào),2007,23(2):393-397. WANG Yanna, HOU Zhen’an, GONG Jiang, et al. Development and expectation of utilization of saline water resources in agriculture irrigation[J]. Chinese Agricultural Science Bulletin,2007,23(2):393-397.(in Chinese)

    4 吳昌娟,張玉龍,虞娜,等.保護(hù)地節(jié)點(diǎn)滲灌下水肥耦合對(duì)土壤全鹽量的影響[J].土壤通報(bào),2013,44(4):897-904. WU Changjuan, ZHANG Yulong, YU Na,et al. Effects of the coupling of water and fertilization on total salt of the soil with node subsurface irrigation in protected field[J]. Chinese Journal of Soil Science,2013,44(4):897-904.(in Chinese)

    5 HU Hongchang, TIAN Fuqiang, HU Heping. Soil particle size distribution and its relationship with soil water and salt under mulched drip irrigation in Xinjiang of China [J]. Science China Technological Sciences, 2011,54(6): 1568-1574.

    6 陳樹(shù)飛,何新林,汪宗飛,等.微咸水滴灌研究進(jìn)展[J].節(jié)水灌溉,2010(2):6-9. CHEN Shufei, HE Xinlin, WANG Zongfei, et al. Development of research on drip irrigation of saline water[J]. Water Saving Irrigation,2010(2):6-9.(in Chinese)

    7 雪靜,王全九,畢遠(yuǎn)杰.微咸水間歇供水土壤入滲特征[J].農(nóng)業(yè)工程學(xué)報(bào),2009,25(5):14-19. XUE Jing, WANG Quanjiu, BI Yuanjie. Soil infiltration properties with saline water intermittent application[J]. Transactions of the CSAE,2009, 25(5): 14-19.(in Chinese)

    8 王丹,康躍虎.微咸水滴灌條件下不同鹽分離子在土壤中的分布特征[J].農(nóng)業(yè)工程學(xué)報(bào),2007,23(2):83-87. WANG Dan, KANG Yuehu. Distribution characteristics of different salt ions in soil under drip irrigation with saline water[J]. Transactions of the CSAE,2007, 23(2): 83-87.(in Chinese)

    9 巨龍,王全九,王琳芳,等.灌水量對(duì)半干旱區(qū)土壤水鹽分布特征及冬小麥產(chǎn)量的影響[J].農(nóng)業(yè)工程學(xué)報(bào),2007,23(1):86-90. JU Long, WANG Quanjiu, WANG Linfang, et al. Effects of irrigation mounts on yield of winter wheat and distribution characteristics of soil water-salt in semi-arid region[J]. Transactions of the CSAE,2007, 23(1): 86-90.(in Chinese)

    10 張偉,呂新,李魯華,等.新疆棉田膜下滴灌鹽分運(yùn)移規(guī)律[J].農(nóng)業(yè)工程學(xué)報(bào),2008,24(8):15-19. ZHANG Wei, Lü Xin, LI Luhua, et al. Salt transfer law for cotton field with drip irrigation under the plastic mulch in Xinjiang region[J]. Transactions of the CSAE,2008, 24(8): 15-19.(in Chinese)

    11 宰松梅,仵峰,溫季,等.不同滴灌方式對(duì)棉田土壤鹽分的影響[J].水利學(xué)報(bào),2011,42(12):1496-1503. ZAI Songmei, WU Feng, WEN Ji, et al. Effect of drip fustigation on soil salinity of cotton field in Northwest China[J]. Journal of Hydraulic Engineering,2011,42(12):1496-1503.(in Chinese)

    12 AMNON B, SHABTAI C, YOEL D M, et al. Effects of timing and duration of brackish irrigation water on fruit yield and quality of late summer melons[J]. Agricultural Water Management, 2005,74(2):123-134.

    13 萬(wàn)書勤,康躍虎,王丹,等.微咸水滴灌對(duì)黃瓜產(chǎn)量及灌溉水利用效率的影響[J].農(nóng)業(yè)工程學(xué)報(bào),2007,23(3):30-35. WAN Shuqin, KANG Yuehu, WANG Dan, et al. Effects of saline water on cucumber yields and irrigation water use efficiency under drip irrigation[J]. Transactions of the CSAE, 2007, 23(3):30-35.(in Chinese)

    14 王全九,葉海燕,史曉楠,等.土壤初始含水量對(duì)微咸水入滲特征影響[J].水土保持學(xué)報(bào),2004,18(1):51-53. WANG Quanjiu, YE Haiyan, SHI Xiaonan, et al. Influence of initial water contention slight saline water infiltration[J]. Journal of Soil and Water Conservation,2004,18(1): 51-53.(in Chinese)

    15 王全九,單魚(yú)洋.微咸水灌溉與土壤水鹽調(diào)控研究進(jìn)展[J/OL].農(nóng)業(yè)機(jī)械學(xué)報(bào),2015,46(12):117-127. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20151217&flag=1. DOI:10.6041/j.issn.1000-1298.2015.12.017. WANG Quanjiu, SHAN Yuyang. Review of research development on water and soil regulation with brackish water irrigation[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2015,46(12):117-127. (in Chinese)

    16 郭太龍, 遲道才, 王全九,等.入滲水礦化度對(duì)土壤水鹽運(yùn)移影響的試驗(yàn)研究[J].農(nóng)業(yè)工程學(xué)報(bào),2005,21(增刊):84-87. GUO Tailong, CHI Daocai, WANG Quanjiu, et al. Experimental study on salt and water movement affected by mineralization degree of infiltration water[J]. Transactions of the CSAE,2005, 21(Supp.): 84-87.(in Chinese)

    17 肖娟,雷廷武,李光永.水質(zhì)及流量對(duì)鹽堿土滴灌濕潤(rùn)鋒運(yùn)移影響的室內(nèi)試驗(yàn)研究[J].農(nóng)業(yè)工程學(xué)報(bào),2007,23(2):88-91. XIAO Juan, LEI Tingwu, LI Guangyong. Laboratory experiment study on the effects of water quality and emitter discharge on wetted soil geometry and volume[J]. Transactions of the CSAE,2007, 23(2): 88-91.(in Chinese)

    18 史曉楠,王全九,蘇瑩.微咸水水質(zhì)對(duì)土壤水鹽運(yùn)移特征的影響[J].干旱區(qū)地理,2005,28(4):516-520. SHI Xiaonan, WANG Quanjiu, SU Ying. Effects of slight salt water quality on the characteristics of soil water and salt transference[J]. Arid Land Geography,2005,28(4):516-520.(in Chinese)

    19 彭振陽(yáng),伍靖?jìng)ィS介生.采用間歇灌溉進(jìn)行土壤鹽分淋洗的適用性[J].水科學(xué)進(jìn)展,2016,27(1):31-39. PENG Zhenyang, WU Jingwei, HUANG Jiesheng. Feasibility of applying intermittent irrigation on solute leaching[J].Advances in Water Science,2016,27(1):31-39. (in Chinese)

    20 魯為華,任愛(ài)天,楊潔晶,等.滴灌苜蓿田間土壤水鹽及苜蓿細(xì)根的空間分布[J].農(nóng)業(yè)工程學(xué)報(bào),2014,30(23):128-137. LU Weihua, REN Aitian, YANG Jiejing, et al. Soil water and salt movement and spatial distribution of fine alfalfa roots under drip irrigation[J]. Transaction of the CSAE, 2014,30(23):128-137.(in Chinese)

    21 劉秀梅,王華田,王延平,等.磁化微咸水灌溉促進(jìn)歐美楊I(lǐng)-107生長(zhǎng)及其光合特性分析[J].農(nóng)業(yè)工程學(xué)報(bào),2016,32(增刊1):1-7. LIU Xiumei, WANG Huatian, WANG Yanping, et al. Analysis of magnetic salinity water irrigation promoting growth and photosynthetic characteristics ofPopulus×euramericanna‘Neva’[J]. Transactions of the CSAE, 2016, 32(Supp.1):1-7.(in Chinese)

    22 張俊,牛文全,張琳琳,等.微潤(rùn)灌溉線源入滲濕潤(rùn)體特性試驗(yàn)研究[J].中國(guó)水土保持科學(xué),2012,10(6):32-38. ZHANG Jun, NIU Wenquan, ZHANG Linlin, et al. Experimental study on characters of wetted soil in moistube irrigation[J]. Science of Soil and Water Conservation,2012,10(6):32-38.(in Chinese)

    23 牛文全,張俊,張琳琳,等.埋深與壓力對(duì)微潤(rùn)灌濕潤(rùn)體水分運(yùn)移的影響[J/OL].農(nóng)業(yè)機(jī)械學(xué)報(bào),2013,44(12):128-134. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20131221&flag=1. DOI:10.6041/j.issn.1000-1298.2013.12.021. NIU Wenquan, ZHANG Jun, ZHANG Linlin,et al. Effects of buried depth and pressure head on water movement of wetted soil during moistube-irrigation[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2013,44(12):128-134.(in Chinese)

    24 NIKOS J W, JAMES W B, KRISTA E P, et al. Basics of salinity and solidity effects on soil physical properties[R]. Bozeman: Montana State University, 2003.

    25 肖振華,萬(wàn)洪福,鄭連芬.灌溉水質(zhì)對(duì)土壤化學(xué)特征和作物生長(zhǎng)的影響[J].土壤學(xué)報(bào),1997,34(3):272-285. XIAO Zhenhua, WANG Hongfu, ZHENG Lianfen. Effect of irrigation water quality on soil chemical charateritics and crop growth[J]. Acta Pedologica Sinica,1997,34(3):272-285.(in Chinese)

    26 吳忠東,王全九.入滲水礦化度對(duì)土壤入滲特征和離子遷移特性的影響[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2010,41(7):64-75. WU Zhongdong, WANG Quanjiu. Effect on both infiltration characteristics and ion mobility features by mineralization degree of infiltration water[J]. Transactions of the Chinese Society for Agricultural Machinery, 2010,41(7):64-75.(in Chinese)

    27 仵峰,吳普特,范永申,等. 地下滴灌條件下土壤水能態(tài)研究[J].農(nóng)業(yè)工程學(xué)報(bào),2008,24(12):31-35. WU Feng, WU Pute, FAN Yongshen, et al. Distribution of soil water potential energy under subsurface drip irrigation[J]. Transactions of the CSAE,2008, 24(12): 31-35.(in Chinese)

    28 陳麗娟,馮超,王昱,等.微咸水灌溉條件下含黏土夾層土壤的水鹽運(yùn)移規(guī)律[J].農(nóng)業(yè)工程學(xué)報(bào),2012,28(8):44-51. CHEN Lijuan, FENG Chao, WANG Yu, et al. Water and salt movement under saline water irrigation in soil with clay interlayer[J]. Transactions of the CSAE,2012,28(8):44-51.(in Chinese)

    29 MOHAMMAD R, EMDAD S R, RAINE R J, et al. Effect of water quality on soil structure and infiltration under furrow irrigation[J]. Irrigation Science,2004,23(2):55-60.

    30 肖振華,萬(wàn)洪富.灌溉水質(zhì)對(duì)土壤水力性質(zhì)和物理性質(zhì)的影響[J].土壤學(xué)報(bào),1998,35(3):359-366. XIAO Zhenhua, WAN Hongfu. Effect of irrigation water quality on soil hydraulic and physical properties[J]. Acta Pedologica Sinica,1998,35(3):359-366.(in Chinese)

    31 吳忠東,王全九.微咸水鈉吸附比對(duì)土壤理化性質(zhì)和入滲特性的影響研究[J].干旱地區(qū)農(nóng)業(yè)研究,2008,26(1):231-236. WU Zhongdong, WANG Quanjiu. Study on impact of sodium adsorption ratio of saline water on soil physical and chemical properties and infiltration characteristics[J]. Agricultural Research in the Arid Areas,2008,26(1):231-236.(in Chinese)

    32 HEYDARI N, DAS A, GUPTA R L. Salinity and solidity influences on infiltration during surge flow irrigation[J]. Irrigation Science,2001,20(4):165-173.

    33 ABDEL G G, ARSLAN A, GAIHBE A, et al. The effects of saline irrigation water management and salt tomato varieties on sustainable production of tomato in Syria [J].Agricultural Water Management,2005,78(12):39-53.

    34 牛文全,薛萬(wàn)來(lái).礦化度對(duì)微潤(rùn)灌溉土壤入滲特性的影響[J/OL].農(nóng)業(yè)機(jī)械學(xué)報(bào),2014,45(4):163-172.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20140426&flag=1. DOI:10.6041/j.issn.1000-1298.2014.04.026. NIU Wenquan, XUE Wanlai. Effects of mineralization degrees on soil infiltration under moistube-irrigation[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2014,45(4):163-172.(in Chinese)

    35 劉勝堯,范鳳翠,李志宏,等.咸水負(fù)壓滲灌對(duì)番茄生長(zhǎng)和土壤鹽分的影響[J].農(nóng)業(yè)工程學(xué)報(bào),2013,29(22):108-117. LIU Shengyao, FAN Fengcui, LI Zhihong, et al. Effects of negative pressure irrigation with saline water on tomato yield and soil salt[J]. Transactions of the CSAE,2013,29(22): 108-117.(in Chinese)

    36 馬東豪,王全九,來(lái)劍斌.膜下滴灌條件下灌水水質(zhì)和流量對(duì)鹽分分布影響的田間試驗(yàn)研究[J].農(nóng)業(yè)工程學(xué)報(bào),2005,21(3):42-46. MA Donghao, WANG Quanjiu, LAI Jianbin. Field experimental studies on the effects of water quality and drip rate on soil salt distribution in drip irrigation under film[J]. Transactions of the CSAE,2005, 21(3): 42-46.(in Chinese)

    Characteristics of Water and Salt Movement in Soil under Moistube-irrigation with Brackish Water

    ZHANG Kemeng1NIU Wenquan1,2WANG Youke1,2XUE Wanlai2,3ZHANG Zizhuo1

    (1.CollegeofWaterResourcesandArchitecturalEngineering,NorthwestA&FUniversity,Yangling,Shaanxi712100,China2.InstituteofSoilandWaterConservation,NorthwestA&FUniversity,Yangling,Shaanxi712100,China
    3.BeijingWaterScienceandTechnologyInstitute,Beijing100048,China)

    In order to explore the movement characteristics of water and salt in soil under moistube-irrgation with brackish water, an indoor soil box simulation experiment was conducted under five different mineralization degree treatments with distilled water as control in 72 h. The results showed that with the rise of mineralization degree, the maximum migration distance presented a trend of first increase and then decrease in different directions and reached the highest when the degree of mineralization was 3.0 g/L. And the distances of wetted front under treatments of brackish water were all greater than that of distilled water. Then the cumulative infiltration was substituted into the Kostiakov infiltration equation and the infiltration coefficient showed a trend of rise first and then fall with the increase of mineralization degree, while infiltration index continued to decline. The soil conductivity was increased constantly to all around with the moistube-irrgation zone as axis and reached the peak at the wet front. The desalting area was related to the shape of the wetted body and distributed in circular shape. There was no significant difference in the average salt content of soil profile between brackish and distilled water after infiltration. The desalination radius was decreased linearly with the increase of mineralization degree. Both surface and bottom convergence phenomenon existed and the salification in soil surface was more serious under moistube-irrgation.

    mineralization degree; infiltration; distributions of water and salt; moistube-irrigation

    10.6041/j.issn.1000-1298.2017.01.023

    2016-07-05

    2016-10-28

    “十二五”國(guó)家科技支撐計(jì)劃項(xiàng)目(2015BAC01B03)和陜北水蝕區(qū)植被功能調(diào)控技術(shù)與示范項(xiàng)目(2015BAC01B03)

    張珂萌(1989—),女,博士生,主要從事水土資源高效利用研究,E-mail: zkmsummerfool@163.com

    汪有科(1956—),男,研究員,博士生導(dǎo)師,主要從事作物高效用水和水土保持研究,E-mail: gjzwyk@vip.sina.com

    S152.7

    A

    1000-1298(2017)01-0175-08

    猜你喜歡
    微潤(rùn)咸水礦化度
    基于水勢(shì)差驅(qū)動(dòng)出流機(jī)制的微潤(rùn)灌水分運(yùn)動(dòng)模型構(gòu)建及模擬
    主題微視頻“微潤(rùn)”學(xué)生成長(zhǎng)
    離子組成及礦化度對(duì)低礦化度水驅(qū)采收率的影響
    巖性油氣藏(2019年5期)2019-09-11 03:54:44
    聊城市地下咸水地質(zhì)特征與綜合開(kāi)發(fā)利用分析
    微潤(rùn)灌水器深埋條件下土壤水分入滲濕潤(rùn)特性試驗(yàn)研究
    微潤(rùn)灌溉系統(tǒng)堵塞評(píng)價(jià)及處理方法研究
    微咸水滴灌能提高紅棗果實(shí)品質(zhì)
    利用綜合測(cè)井資料計(jì)算地下水等效NaCl溶液礦化度方法
    低礦化度水驅(qū)技術(shù)增產(chǎn)機(jī)理與適用條件
    特高礦化度Cr3+交聯(lián)聚合物溶液滲流特性及其機(jī)制
    亚洲av一区综合| 久久久久久久久久久丰满 | 精品不卡国产一区二区三区| 成人二区视频| 欧美另类亚洲清纯唯美| 美女高潮的动态| 亚洲男人的天堂狠狠| 亚洲中文日韩欧美视频| 又粗又爽又猛毛片免费看| 国产高清激情床上av| 久久精品综合一区二区三区| 欧美色视频一区免费| 日本在线视频免费播放| 一级黄色大片毛片| 深夜a级毛片| 日日撸夜夜添| 国产精品一及| 亚洲精品一区av在线观看| 小说图片视频综合网站| 久久久久久久久久成人| 国内精品一区二区在线观看| 免费av不卡在线播放| 人妻夜夜爽99麻豆av| 神马国产精品三级电影在线观看| 女人十人毛片免费观看3o分钟| 精品一区二区免费观看| 他把我摸到了高潮在线观看| 久久这里只有精品中国| 香蕉av资源在线| 午夜福利在线在线| 变态另类成人亚洲欧美熟女| 国产毛片a区久久久久| 成人二区视频| 人妻制服诱惑在线中文字幕| 国产色婷婷99| 中文字幕av在线有码专区| 日韩大尺度精品在线看网址| 成人特级av手机在线观看| 精品午夜福利在线看| 色在线成人网| 变态另类成人亚洲欧美熟女| 国国产精品蜜臀av免费| 日日撸夜夜添| 男女边吃奶边做爰视频| 日本免费a在线| 亚洲不卡免费看| 成人亚洲精品av一区二区| 成人国产麻豆网| 国产伦在线观看视频一区| 欧洲精品卡2卡3卡4卡5卡区| 久久久久免费精品人妻一区二区| 观看免费一级毛片| 又黄又爽又刺激的免费视频.| 最近最新中文字幕大全电影3| 免费看a级黄色片| 日日啪夜夜撸| av.在线天堂| 99九九线精品视频在线观看视频| 草草在线视频免费看| 精品人妻视频免费看| 免费观看在线日韩| 日本黄色片子视频| 欧美日韩黄片免| 国产视频一区二区在线看| 精品久久久久久久久av| 精品久久久久久久久av| 亚洲av免费在线观看| 18禁黄网站禁片免费观看直播| 国产av在哪里看| 婷婷精品国产亚洲av| 欧美人与善性xxx| netflix在线观看网站| 黄色配什么色好看| 成年女人看的毛片在线观看| 波多野结衣高清无吗| 久久99热这里只有精品18| 97热精品久久久久久| 日韩欧美在线二视频| 国产白丝娇喘喷水9色精品| eeuss影院久久| 制服丝袜大香蕉在线| 夜夜夜夜夜久久久久| 91狼人影院| 亚洲乱码一区二区免费版| 国产成年人精品一区二区| 精品久久久久久久末码| 淫妇啪啪啪对白视频| 亚洲第一区二区三区不卡| 免费不卡的大黄色大毛片视频在线观看 | 国产亚洲91精品色在线| 十八禁网站免费在线| 国产精品av视频在线免费观看| 99久久中文字幕三级久久日本| 婷婷精品国产亚洲av在线| 欧美+亚洲+日韩+国产| 悠悠久久av| 欧美又色又爽又黄视频| 亚洲精品一区av在线观看| 免费看a级黄色片| 亚洲精品色激情综合| 乱人视频在线观看| 国产精品久久视频播放| 又爽又黄无遮挡网站| 国产精品99久久久久久久久| 国内揄拍国产精品人妻在线| 国产欧美日韩一区二区精品| 亚洲成av人片在线播放无| 午夜福利18| 亚洲av一区综合| 免费看av在线观看网站| 亚洲经典国产精华液单| 日日撸夜夜添| 国产成人av教育| 日韩欧美在线乱码| 国内揄拍国产精品人妻在线| 观看美女的网站| 真人一进一出gif抽搐免费| 亚洲精品一区av在线观看| 色噜噜av男人的天堂激情| 一区福利在线观看| 丰满的人妻完整版| 国产精品av视频在线免费观看| 国产精品无大码| 亚洲国产欧美人成| 狂野欧美激情性xxxx在线观看| av在线天堂中文字幕| 男女之事视频高清在线观看| а√天堂www在线а√下载| 日本黄大片高清| 91久久精品电影网| 亚洲av中文av极速乱 | 久久久久国产精品人妻aⅴ院| 欧美日韩黄片免| 国产精品久久视频播放| 亚洲精品影视一区二区三区av| 淫秽高清视频在线观看| 亚洲综合色惰| 18+在线观看网站| 日日摸夜夜添夜夜添av毛片 | 波多野结衣高清无吗| 又紧又爽又黄一区二区| 亚洲电影在线观看av| 色哟哟哟哟哟哟| 国内少妇人妻偷人精品xxx网站| 亚洲成人精品中文字幕电影| 天堂动漫精品| 日本a在线网址| 人妻少妇偷人精品九色| 可以在线观看毛片的网站| 一个人看的www免费观看视频| 亚洲 国产 在线| 国产精品av视频在线免费观看| 亚洲自偷自拍三级| 国产一区二区三区在线臀色熟女| 日本一二三区视频观看| 午夜福利成人在线免费观看| 欧美高清成人免费视频www| 国产伦在线观看视频一区| 国产高清不卡午夜福利| 91麻豆精品激情在线观看国产| 狠狠狠狠99中文字幕| 国产亚洲精品久久久com| 色综合色国产| 亚洲精品乱码久久久v下载方式| 亚洲欧美精品综合久久99| 少妇裸体淫交视频免费看高清| 啪啪无遮挡十八禁网站| 国产精品久久久久久久久免| 国产伦精品一区二区三区四那| 熟妇人妻久久中文字幕3abv| 国产亚洲av嫩草精品影院| 日韩亚洲欧美综合| 99久久成人亚洲精品观看| 亚洲无线观看免费| 国产精品综合久久久久久久免费| 日本与韩国留学比较| 精品国产三级普通话版| 亚洲欧美日韩高清在线视频| 特级一级黄色大片| 一区二区三区四区激情视频 | 欧美日韩瑟瑟在线播放| 亚洲一区二区三区色噜噜| 国产一区二区三区av在线 | 日韩强制内射视频| 又紧又爽又黄一区二区| 99国产极品粉嫩在线观看| av中文乱码字幕在线| 99在线视频只有这里精品首页| 很黄的视频免费| 少妇人妻精品综合一区二区 | 深夜精品福利| 91久久精品国产一区二区三区| 国产一区二区亚洲精品在线观看| 国产精品电影一区二区三区| 看免费成人av毛片| 国语自产精品视频在线第100页| 中文字幕免费在线视频6| 国产aⅴ精品一区二区三区波| 别揉我奶头 嗯啊视频| 如何舔出高潮| 在线国产一区二区在线| 国产毛片a区久久久久| 成人毛片a级毛片在线播放| 97人妻精品一区二区三区麻豆| 美女高潮喷水抽搐中文字幕| 久久精品夜夜夜夜夜久久蜜豆| 色综合亚洲欧美另类图片| 国产精品女同一区二区软件 | 欧美日韩国产亚洲二区| 国产精品免费一区二区三区在线| 啦啦啦啦在线视频资源| 少妇熟女aⅴ在线视频| 99国产极品粉嫩在线观看| 一级毛片久久久久久久久女| 精品人妻熟女av久视频| 国产又黄又爽又无遮挡在线| 国产精品久久久久久久久免| 日韩一区二区视频免费看| 国产黄片美女视频| av国产免费在线观看| 亚洲精品粉嫩美女一区| 女的被弄到高潮叫床怎么办 | 亚洲最大成人手机在线| 国产成人a区在线观看| 三级男女做爰猛烈吃奶摸视频| 91久久精品电影网| www.色视频.com| 色综合站精品国产| 国产淫片久久久久久久久| 日本三级黄在线观看| 我要搜黄色片| 99在线视频只有这里精品首页| 免费av观看视频| 国产不卡一卡二| 男插女下体视频免费在线播放| 真人做人爱边吃奶动态| 精品久久久久久久久亚洲 | 一本精品99久久精品77| 在线观看66精品国产| 亚洲av免费在线观看| 日本一二三区视频观看| 一本精品99久久精品77| 神马国产精品三级电影在线观看| 亚洲av一区综合| 99热这里只有精品一区| 夜夜看夜夜爽夜夜摸| 长腿黑丝高跟| 亚洲va在线va天堂va国产| 直男gayav资源| 国产女主播在线喷水免费视频网站 | 久9热在线精品视频| 内地一区二区视频在线| 国产乱人伦免费视频| 国产真实伦视频高清在线观看 | 男女那种视频在线观看| 动漫黄色视频在线观看| 国产伦精品一区二区三区视频9| 色播亚洲综合网| 老女人水多毛片| 男人舔奶头视频| 观看美女的网站| 精品99又大又爽又粗少妇毛片 | 亚洲精品国产成人久久av| 久久精品国产亚洲网站| 久久精品国产亚洲av涩爱 | 国产在线精品亚洲第一网站| 亚洲专区国产一区二区| 亚洲成a人片在线一区二区| 精品国内亚洲2022精品成人| 日韩欧美 国产精品| 精品一区二区免费观看| 午夜老司机福利剧场| 91在线观看av| 国产成年人精品一区二区| 特大巨黑吊av在线直播| 亚洲美女黄片视频| 国产麻豆成人av免费视频| 在线天堂最新版资源| 精品免费久久久久久久清纯| 少妇的逼水好多| 久久中文看片网| 国产亚洲欧美98| 免费在线观看影片大全网站| 亚洲黑人精品在线| 色综合色国产| 日本a在线网址| 亚洲av成人精品一区久久| 国产一区二区三区av在线 | 成年人黄色毛片网站| 欧美激情久久久久久爽电影| 欧美一区二区国产精品久久精品| 国产在线精品亚洲第一网站| 久久久久久伊人网av| 夜夜夜夜夜久久久久| 不卡一级毛片| 久久国内精品自在自线图片| 婷婷精品国产亚洲av| 超碰av人人做人人爽久久| 黄色欧美视频在线观看| 久久精品国产清高在天天线| 一卡2卡三卡四卡精品乱码亚洲| 欧美绝顶高潮抽搐喷水| 国内毛片毛片毛片毛片毛片| 国产亚洲精品av在线| 亚洲av成人av| 最近最新中文字幕大全电影3| 国产一区二区在线观看日韩| 很黄的视频免费| av在线蜜桃| 欧美性猛交黑人性爽| 中文在线观看免费www的网站| av中文乱码字幕在线| 午夜福利视频1000在线观看| 欧美最新免费一区二区三区| 免费搜索国产男女视频| 欧美人与善性xxx| 国产一区二区亚洲精品在线观看| 免费人成在线观看视频色| 久久精品国产亚洲av香蕉五月| 欧美中文日本在线观看视频| 亚洲美女黄片视频| 欧美色欧美亚洲另类二区| 动漫黄色视频在线观看| 国产成人影院久久av| 黄色配什么色好看| 日本免费a在线| 午夜视频国产福利| 欧美精品国产亚洲| 淫妇啪啪啪对白视频| 草草在线视频免费看| 国产视频内射| 22中文网久久字幕| 在现免费观看毛片| 国产亚洲精品久久久久久毛片| 嫩草影院新地址| 亚洲av中文字字幕乱码综合| 久久人人爽人人爽人人片va| 99久久成人亚洲精品观看| 18禁黄网站禁片免费观看直播| 美女 人体艺术 gogo| 成年女人看的毛片在线观看| 久久久久九九精品影院| 国产视频一区二区在线看| 丝袜美腿在线中文| 黄色一级大片看看| 日本爱情动作片www.在线观看 | 日本免费一区二区三区高清不卡| 麻豆国产av国片精品| 麻豆成人av在线观看| 一级毛片久久久久久久久女| 中国美白少妇内射xxxbb| 女生性感内裤真人,穿戴方法视频| 天堂影院成人在线观看| 99热网站在线观看| 亚洲不卡免费看| www日本黄色视频网| 亚洲欧美激情综合另类| 又粗又爽又猛毛片免费看| 亚洲av免费高清在线观看| 伊人久久精品亚洲午夜| 一个人免费在线观看电影| 精品久久久久久久人妻蜜臀av| 88av欧美| 国产色婷婷99| 亚洲精品粉嫩美女一区| 真人做人爱边吃奶动态| 狂野欧美激情性xxxx在线观看| 人人妻人人看人人澡| 日日撸夜夜添| 又爽又黄a免费视频| 草草在线视频免费看| 亚洲人成网站在线播| 最好的美女福利视频网| 久久久国产成人精品二区| 2021天堂中文幕一二区在线观| 日本欧美国产在线视频| 国产女主播在线喷水免费视频网站 | 一个人免费在线观看电影| a级一级毛片免费在线观看| 久久久成人免费电影| 九九爱精品视频在线观看| 黄片wwwwww| 免费看日本二区| 少妇人妻一区二区三区视频| 欧美中文日本在线观看视频| 精品久久久久久久久亚洲 | 久久精品夜夜夜夜夜久久蜜豆| 国产单亲对白刺激| 五月伊人婷婷丁香| 久99久视频精品免费| 成人国产一区最新在线观看| xxxwww97欧美| 免费看a级黄色片| 久久久久性生活片| 日本与韩国留学比较| 国产亚洲欧美98| 自拍偷自拍亚洲精品老妇| 亚洲成a人片在线一区二区| 人人妻人人澡欧美一区二区| 桃色一区二区三区在线观看| 一a级毛片在线观看| 一卡2卡三卡四卡精品乱码亚洲| 色噜噜av男人的天堂激情| 色av中文字幕| 国产精品久久久久久精品电影| 99热这里只有精品一区| 一进一出好大好爽视频| 亚洲人成网站高清观看| 亚洲18禁久久av| 国产成人a区在线观看| 五月玫瑰六月丁香| av黄色大香蕉| 99视频精品全部免费 在线| 观看美女的网站| 亚洲va日本ⅴa欧美va伊人久久| 国产精品爽爽va在线观看网站| 亚洲图色成人| 男人舔女人下体高潮全视频| 亚洲国产色片| 嫩草影视91久久| 国产精品嫩草影院av在线观看 | 嫁个100分男人电影在线观看| 欧美日韩国产亚洲二区| 男女视频在线观看网站免费| 一个人看的www免费观看视频| 在线观看美女被高潮喷水网站| 伦精品一区二区三区| 69人妻影院| 看十八女毛片水多多多| 久久久久九九精品影院| 国产老妇女一区| 欧美黑人巨大hd| 九九在线视频观看精品| 亚洲欧美日韩高清专用| 日本免费一区二区三区高清不卡| xxxwww97欧美| 国产高清激情床上av| 国内久久婷婷六月综合欲色啪| 亚洲av免费在线观看| 99热这里只有是精品50| 亚洲av日韩精品久久久久久密| 国产亚洲精品av在线| 乱码一卡2卡4卡精品| 自拍偷自拍亚洲精品老妇| 国产精品嫩草影院av在线观看 | 国产蜜桃级精品一区二区三区| 99精品久久久久人妻精品| 成人高潮视频无遮挡免费网站| 亚洲av.av天堂| 91麻豆av在线| 3wmmmm亚洲av在线观看| 国产亚洲欧美98| 欧美日韩精品成人综合77777| 欧美另类亚洲清纯唯美| 色吧在线观看| 男女那种视频在线观看| 亚洲无线在线观看| 欧美精品啪啪一区二区三区| 99riav亚洲国产免费| 国内精品美女久久久久久| 在线免费十八禁| 村上凉子中文字幕在线| 中文字幕免费在线视频6| 香蕉av资源在线| 精品一区二区三区av网在线观看| 丰满乱子伦码专区| 久久久久久久久久成人| 国产精品久久视频播放| 国产高清视频在线播放一区| 亚洲精品色激情综合| 午夜久久久久精精品| 国产精品三级大全| 搞女人的毛片| 一边摸一边抽搐一进一小说| 日本黄大片高清| 97超级碰碰碰精品色视频在线观看| 美女免费视频网站| 国产大屁股一区二区在线视频| 麻豆国产97在线/欧美| 一个人看视频在线观看www免费| 天天躁日日操中文字幕| 成人毛片a级毛片在线播放| 国产69精品久久久久777片| 精品国产三级普通话版| 窝窝影院91人妻| 校园人妻丝袜中文字幕| 亚洲电影在线观看av| 国内精品久久久久久久电影| 在线观看一区二区三区| 国产成人a区在线观看| 亚洲性夜色夜夜综合| 婷婷亚洲欧美| 成人性生交大片免费视频hd| 亚洲美女黄片视频| 国产单亲对白刺激| av国产免费在线观看| 亚洲av第一区精品v没综合| 韩国av在线不卡| 日韩欧美国产一区二区入口| 三级毛片av免费| 国产三级中文精品| 午夜激情福利司机影院| 午夜福利高清视频| 亚洲av免费高清在线观看| 亚洲欧美日韩东京热| 嫩草影院新地址| 熟妇人妻久久中文字幕3abv| 国内毛片毛片毛片毛片毛片| 欧美日本视频| 无遮挡黄片免费观看| 99热只有精品国产| 真实男女啪啪啪动态图| 国产私拍福利视频在线观看| 久久精品国产鲁丝片午夜精品 | 男人和女人高潮做爰伦理| 国产老妇女一区| 午夜老司机福利剧场| 亚洲精品一卡2卡三卡4卡5卡| 欧美激情在线99| 亚洲色图av天堂| 啪啪无遮挡十八禁网站| 欧美激情久久久久久爽电影| 日韩av在线大香蕉| 亚洲av五月六月丁香网| 97热精品久久久久久| 亚洲av二区三区四区| 夜夜爽天天搞| 黄色一级大片看看| 亚洲国产精品成人综合色| 天天一区二区日本电影三级| 校园人妻丝袜中文字幕| 内射极品少妇av片p| 一级毛片久久久久久久久女| 久久久久久久久久成人| 一区二区三区免费毛片| av在线观看视频网站免费| 国产综合懂色| 少妇裸体淫交视频免费看高清| 欧美三级亚洲精品| 老司机午夜福利在线观看视频| 国产精品一区二区三区四区久久| 国产伦一二天堂av在线观看| 亚洲四区av| 自拍偷自拍亚洲精品老妇| 国产 一区 欧美 日韩| 亚洲狠狠婷婷综合久久图片| 熟女电影av网| 国产精品伦人一区二区| 天堂√8在线中文| 老司机深夜福利视频在线观看| 精品欧美国产一区二区三| 国产视频一区二区在线看| bbb黄色大片| 丰满的人妻完整版| 国产伦一二天堂av在线观看| 久久久色成人| avwww免费| 国产伦精品一区二区三区四那| 欧美黑人巨大hd| 日韩欧美 国产精品| 亚洲精品456在线播放app | 老司机深夜福利视频在线观看| 成人午夜高清在线视频| 波野结衣二区三区在线| x7x7x7水蜜桃| 在线观看一区二区三区| 看黄色毛片网站| 夜夜爽天天搞| 身体一侧抽搐| 精品无人区乱码1区二区| 国产主播在线观看一区二区| 真人做人爱边吃奶动态| www.www免费av| 日本免费一区二区三区高清不卡| 99久国产av精品| 99国产精品一区二区蜜桃av| 欧美区成人在线视频| 性欧美人与动物交配| 天堂动漫精品| 免费无遮挡裸体视频| 久久精品综合一区二区三区| 国产伦在线观看视频一区| 1024手机看黄色片| 欧美日韩精品成人综合77777| 国产av麻豆久久久久久久| 99在线人妻在线中文字幕| 联通29元200g的流量卡| 精品久久久久久久末码| 亚洲美女黄片视频| 国产精品电影一区二区三区| 搡老岳熟女国产| 国产中年淑女户外野战色| 婷婷精品国产亚洲av在线| 禁无遮挡网站| 国产单亲对白刺激| 亚洲色图av天堂| 97人妻精品一区二区三区麻豆| 亚洲av不卡在线观看| 国产免费av片在线观看野外av| 中文字幕精品亚洲无线码一区| 1000部很黄的大片| 超碰av人人做人人爽久久| 少妇高潮的动态图| 欧美xxxx黑人xx丫x性爽| 欧美人与善性xxx| 成人午夜高清在线视频| 国产单亲对白刺激| 欧美日韩瑟瑟在线播放| 一个人免费在线观看电影| 精品人妻偷拍中文字幕| 精品久久久久久久久久久久久| 亚洲国产欧美人成| 免费观看在线日韩| 九九爱精品视频在线观看| 波多野结衣高清作品|