• <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| 老司机午夜十八禁免费视频| 母亲3免费完整高清在线观看| 久久精品成人免费网站| 美女高潮到喷水免费观看| 久久久国产精品麻豆| 一级毛片黄色毛片免费观看视频| 欧美亚洲 丝袜 人妻 在线| 性色av乱码一区二区三区2| 亚洲国产欧美一区二区综合| 欧美日韩成人在线一区二区| 国产成人影院久久av| 精品久久久久久久毛片微露脸 | 两人在一起打扑克的视频| 国产亚洲欧美在线一区二区| 成人国产一区最新在线观看 | 美女国产高潮福利片在线看| 黑人欧美特级aaaaaa片| 久久久久久久久免费视频了| 男女高潮啪啪啪动态图| 麻豆国产av国片精品| 亚洲人成电影观看| 国产福利在线免费观看视频| 国产免费福利视频在线观看| 波野结衣二区三区在线| 亚洲国产日韩一区二区| 精品少妇黑人巨大在线播放| 久久久欧美国产精品| 免费看不卡的av| h视频一区二区三区| 丁香六月天网| 亚洲精品国产av蜜桃| 亚洲 国产 在线| 国产精品久久久久久精品电影小说| 午夜福利影视在线免费观看| 在线 av 中文字幕| 久久午夜综合久久蜜桃| 久久久久久久久久久久大奶| 好男人电影高清在线观看| 无限看片的www在线观看| av一本久久久久| 国产极品粉嫩免费观看在线| 高清视频免费观看一区二区| 国产在线观看jvid| 国产爽快片一区二区三区| 亚洲国产欧美日韩在线播放| 一本综合久久免费| 午夜福利影视在线免费观看| 黄色毛片三级朝国网站| 日韩av免费高清视频| 日本猛色少妇xxxxx猛交久久| 国产精品免费大片| 狂野欧美激情性xxxx| 十分钟在线观看高清视频www| 久久久久久亚洲精品国产蜜桃av| 一区二区av电影网| 国产福利在线免费观看视频| 亚洲欧美精品自产自拍| 成人免费观看视频高清| 国产午夜精品一二区理论片| 男女无遮挡免费网站观看| 日本欧美视频一区| 亚洲成人免费av在线播放| 狂野欧美激情性bbbbbb| 少妇猛男粗大的猛烈进出视频| 久久影院123| 搡老乐熟女国产| 久久影院123| 亚洲av电影在线观看一区二区三区| 成人三级做爰电影| 免费女性裸体啪啪无遮挡网站| 精品少妇黑人巨大在线播放| 欧美激情极品国产一区二区三区| 精品亚洲成国产av| 亚洲欧美精品综合一区二区三区| 亚洲自偷自拍图片 自拍| 亚洲国产中文字幕在线视频| 国产一区亚洲一区在线观看| 91麻豆av在线| 国产一区亚洲一区在线观看| 国产黄色免费在线视频| 久久久精品免费免费高清| 亚洲精品久久成人aⅴ小说| 侵犯人妻中文字幕一二三四区| 国语对白做爰xxxⅹ性视频网站| 精品第一国产精品| 美女扒开内裤让男人捅视频| 久久这里只有精品19| 99热国产这里只有精品6| 日韩av不卡免费在线播放| 男女无遮挡免费网站观看| avwww免费| 一区二区三区激情视频| 亚洲欧美中文字幕日韩二区| av国产久精品久网站免费入址| av国产精品久久久久影院| 黄色片一级片一级黄色片| 一本一本久久a久久精品综合妖精| 纯流量卡能插随身wifi吗| 麻豆乱淫一区二区| 女性被躁到高潮视频| 午夜91福利影院| 久久午夜综合久久蜜桃| 国产欧美日韩综合在线一区二区| 91麻豆精品激情在线观看国产 | 久久国产亚洲av麻豆专区| 中文字幕高清在线视频| 激情视频va一区二区三区| 亚洲成人免费电影在线观看 | 国产精品久久久久久精品电影小说| 中文字幕最新亚洲高清| 极品人妻少妇av视频| 日韩av在线免费看完整版不卡| 国产av精品麻豆| 深夜精品福利| 成年动漫av网址| 久久人人97超碰香蕉20202| 国产亚洲av高清不卡| 最近手机中文字幕大全| av国产精品久久久久影院| 亚洲国产中文字幕在线视频| 老司机午夜十八禁免费视频| 国产精品国产三级专区第一集| 国产成人精品久久久久久| 欧美97在线视频| 黄色 视频免费看| 欧美国产精品va在线观看不卡| 国产日韩一区二区三区精品不卡| 最新在线观看一区二区三区 | 天天添夜夜摸| 国产淫语在线视频| 999精品在线视频| 一边摸一边抽搐一进一出视频| 黄色毛片三级朝国网站| 男女边吃奶边做爰视频| 久久久久久亚洲精品国产蜜桃av| 亚洲五月婷婷丁香| cao死你这个sao货| 波野结衣二区三区在线| 日韩精品免费视频一区二区三区| 好男人电影高清在线观看| 另类精品久久| 午夜免费鲁丝| 啦啦啦在线观看免费高清www| 视频区图区小说| 日本一区二区免费在线视频| 两性夫妻黄色片| cao死你这个sao货| 99国产精品一区二区蜜桃av | 免费人妻精品一区二区三区视频| 精品免费久久久久久久清纯 | 日韩熟女老妇一区二区性免费视频| 国产亚洲一区二区精品| 国产精品国产三级国产专区5o| 中文字幕精品免费在线观看视频| 亚洲成人手机| 不卡av一区二区三区| 狂野欧美激情性bbbbbb| 飞空精品影院首页| 精品亚洲成国产av| 久久久精品免费免费高清| 亚洲久久久国产精品| 天堂中文最新版在线下载| 国产精品av久久久久免费| 大香蕉久久网| 91成人精品电影| 激情视频va一区二区三区| 女人被躁到高潮嗷嗷叫费观| 色婷婷久久久亚洲欧美| 精品高清国产在线一区| 一区二区三区四区激情视频| 操美女的视频在线观看| 搡老乐熟女国产| 久久精品久久精品一区二区三区| 人妻人人澡人人爽人人| 在线观看免费高清a一片| 校园人妻丝袜中文字幕| 国产免费一区二区三区四区乱码| 不卡av一区二区三区| 中文字幕精品免费在线观看视频| 在线 av 中文字幕| 久热爱精品视频在线9| 国产视频一区二区在线看| 美女视频免费永久观看网站| 久久精品国产a三级三级三级| 人人妻人人添人人爽欧美一区卜| 欧美性长视频在线观看| videos熟女内射| 成年动漫av网址| 国产欧美亚洲国产| 国产一区二区激情短视频 | 国产在线一区二区三区精| 三上悠亚av全集在线观看| 国产精品九九99| 国产欧美日韩综合在线一区二区| 91九色精品人成在线观看| 在线亚洲精品国产二区图片欧美| 亚洲av美国av| 欧美亚洲 丝袜 人妻 在线| 成人影院久久| 欧美成人精品欧美一级黄| 国产在线视频一区二区| 国产黄频视频在线观看| av福利片在线| 亚洲情色 制服丝袜| 亚洲成人免费电影在线观看 | 成人手机av| 在线观看国产h片| a 毛片基地| 赤兔流量卡办理| 国产成人影院久久av| 青春草亚洲视频在线观看| 国产视频一区二区在线看| 日本五十路高清| 在线观看免费视频网站a站| 国产高清不卡午夜福利| 精品欧美一区二区三区在线| 曰老女人黄片| 蜜桃在线观看..| 亚洲av国产av综合av卡| 校园人妻丝袜中文字幕| 久久久精品区二区三区| 大陆偷拍与自拍| 首页视频小说图片口味搜索 | 亚洲,欧美,日韩| 欧美成狂野欧美在线观看| 久久久久视频综合| 国产在视频线精品| 在线观看一区二区三区激情| 黄频高清免费视频| 国产精品 欧美亚洲| 国产精品免费视频内射| tube8黄色片| 成年人黄色毛片网站| 脱女人内裤的视频| av线在线观看网站| 伊人亚洲综合成人网| 又大又黄又爽视频免费| 亚洲 欧美一区二区三区| 18禁观看日本| 交换朋友夫妻互换小说| 黄色怎么调成土黄色| 飞空精品影院首页| 美女高潮到喷水免费观看| 久久久精品94久久精品| 飞空精品影院首页| 狠狠婷婷综合久久久久久88av| 91国产中文字幕| 在线 av 中文字幕| 国产成人一区二区在线| 人人妻人人澡人人看| 一区二区三区乱码不卡18| 巨乳人妻的诱惑在线观看| 色精品久久人妻99蜜桃| 久久ye,这里只有精品| 亚洲精品在线美女| 亚洲成色77777| 亚洲av男天堂| 亚洲欧美一区二区三区国产| 黄色一级大片看看| 亚洲伊人久久精品综合| 亚洲欧美中文字幕日韩二区| 日韩,欧美,国产一区二区三区| 18禁黄网站禁片午夜丰满| 9191精品国产免费久久| 日韩制服丝袜自拍偷拍| 大型av网站在线播放| 欧美国产精品一级二级三级| 国产无遮挡羞羞视频在线观看| 看免费成人av毛片| 亚洲免费av在线视频| 久久精品亚洲熟妇少妇任你| 赤兔流量卡办理| 国产精品成人在线| 中文字幕av电影在线播放| 在线亚洲精品国产二区图片欧美| 国产精品av久久久久免费| 精品欧美一区二区三区在线| av不卡在线播放| 精品视频人人做人人爽| 欧美日韩成人在线一区二区| 亚洲欧美一区二区三区久久| 91麻豆精品激情在线观看国产 | 50天的宝宝边吃奶边哭怎么回事| 国产真人三级小视频在线观看| 欧美黑人欧美精品刺激| 人成视频在线观看免费观看| 亚洲伊人久久精品综合| 777米奇影视久久| 青青草视频在线视频观看| 一区二区三区乱码不卡18| 黄网站色视频无遮挡免费观看| 久久久久久久国产电影| 精品亚洲乱码少妇综合久久| 黄片小视频在线播放| videosex国产| 亚洲熟女精品中文字幕| 午夜久久久在线观看| 91精品伊人久久大香线蕉| 亚洲欧美日韩另类电影网站| 亚洲欧美色中文字幕在线| 欧美激情 高清一区二区三区| 亚洲中文av在线| 丝袜人妻中文字幕| 亚洲精品自拍成人| 99热网站在线观看| 国产99久久九九免费精品| 亚洲av成人不卡在线观看播放网 | 国产精品人妻久久久影院| 欧美日韩国产mv在线观看视频| 亚洲国产毛片av蜜桃av| 精品久久久久久久毛片微露脸 | 国产成人一区二区在线| 少妇人妻 视频| 热99国产精品久久久久久7| avwww免费| 婷婷色麻豆天堂久久| 亚洲欧美日韩高清在线视频 | 在线av久久热| 嫩草影视91久久| 欧美日韩亚洲综合一区二区三区_| 18禁国产床啪视频网站| 免费看十八禁软件| 啦啦啦中文免费视频观看日本| 在线观看人妻少妇| 99精品久久久久人妻精品| 亚洲自偷自拍图片 自拍| 国产午夜精品一二区理论片| 亚洲欧美成人综合另类久久久| 免费一级毛片在线播放高清视频 | 亚洲国产欧美日韩在线播放| 男女之事视频高清在线观看 | 老司机深夜福利视频在线观看 | 在线精品无人区一区二区三| 亚洲欧洲国产日韩| 中文乱码字字幕精品一区二区三区| 肉色欧美久久久久久久蜜桃| 无限看片的www在线观看| 成人亚洲精品一区在线观看| 精品卡一卡二卡四卡免费| 捣出白浆h1v1| 丰满迷人的少妇在线观看| 国产不卡av网站在线观看| 免费久久久久久久精品成人欧美视频| 欧美黑人精品巨大| av一本久久久久| 考比视频在线观看| 男女边吃奶边做爰视频| 国产又色又爽无遮挡免| 国产视频首页在线观看| av视频免费观看在线观看| 成人三级做爰电影| 日韩 欧美 亚洲 中文字幕| 黄色片一级片一级黄色片| 国产黄色视频一区二区在线观看| 尾随美女入室| 超色免费av| 亚洲情色 制服丝袜| 欧美成人午夜精品| 一区二区日韩欧美中文字幕| 高清黄色对白视频在线免费看| 久久久久视频综合| 王馨瑶露胸无遮挡在线观看| 三上悠亚av全集在线观看| 日本黄色日本黄色录像| 99精品久久久久人妻精品| 大片免费播放器 马上看| 亚洲精品一卡2卡三卡4卡5卡 | 国产熟女欧美一区二区| 丰满人妻熟妇乱又伦精品不卡| 欧美 日韩 精品 国产| 精品免费久久久久久久清纯 | 女人精品久久久久毛片| 9热在线视频观看99| www.999成人在线观看| 中文字幕人妻熟女乱码| 黄色视频在线播放观看不卡| 啦啦啦 在线观看视频| 国产国语露脸激情在线看| 午夜激情av网站| 精品久久久久久久毛片微露脸 | 久久人人97超碰香蕉20202| 亚洲av欧美aⅴ国产| av国产精品久久久久影院| 精品一区二区三卡| 国产成人一区二区三区免费视频网站 | 久久久久国产精品人妻一区二区| 日韩欧美一区视频在线观看| 狠狠婷婷综合久久久久久88av| 日本wwww免费看| 久久 成人 亚洲| 91字幕亚洲| 老汉色av国产亚洲站长工具| 国产一区亚洲一区在线观看| 看免费av毛片| 国产av国产精品国产| 叶爱在线成人免费视频播放| 国产成人影院久久av| 99久久99久久久精品蜜桃| 亚洲,欧美精品.| 99热网站在线观看| 一边摸一边做爽爽视频免费| 亚洲人成网站在线观看播放| 18禁裸乳无遮挡动漫免费视频| 在线 av 中文字幕| 国产男女超爽视频在线观看| 国产一区二区 视频在线| videos熟女内射| 波多野结衣一区麻豆| 亚洲国产精品成人久久小说| www.av在线官网国产| 日韩av免费高清视频| 999精品在线视频| 久久久亚洲精品成人影院| 人妻人人澡人人爽人人| 老汉色av国产亚洲站长工具| www.熟女人妻精品国产| 国产真人三级小视频在线观看| 热re99久久精品国产66热6| 久久国产精品大桥未久av| 中文字幕最新亚洲高清| 精品一区二区三卡| 日本vs欧美在线观看视频| 一级a爱视频在线免费观看| 最近中文字幕2019免费版| 国产成人av激情在线播放| 日韩免费高清中文字幕av| kizo精华| 欧美黄色淫秽网站| 精品国产超薄肉色丝袜足j| 国产伦人伦偷精品视频| 亚洲中文日韩欧美视频| 男女午夜视频在线观看| 欧美激情高清一区二区三区| 在线观看免费日韩欧美大片| 色视频在线一区二区三区| 亚洲精品国产区一区二| 亚洲精品一二三| 午夜老司机福利片| 天堂俺去俺来也www色官网| 18禁黄网站禁片午夜丰满| 真人做人爱边吃奶动态| 少妇被粗大的猛进出69影院| 两人在一起打扑克的视频| 性色av一级| 这个男人来自地球电影免费观看| 亚洲自偷自拍图片 自拍| cao死你这个sao货| 在线观看www视频免费| 99热国产这里只有精品6| 精品亚洲成国产av| www.熟女人妻精品国产| 国产欧美亚洲国产| 91九色精品人成在线观看| 天天操日日干夜夜撸| 久久国产精品影院| 日韩伦理黄色片| 午夜福利视频在线观看免费| 午夜久久久在线观看| 人妻一区二区av| 日日摸夜夜添夜夜爱| 香蕉丝袜av| 九色亚洲精品在线播放| 99香蕉大伊视频| 精品亚洲成国产av| 久久久久精品国产欧美久久久 | 亚洲精品久久成人aⅴ小说| 久久热在线av| 久久久国产精品麻豆| 国产精品熟女久久久久浪| 日韩一本色道免费dvd| 黄色 视频免费看| 亚洲国产av新网站| h视频一区二区三区| 亚洲国产精品一区三区| 成人免费观看视频高清| 亚洲午夜精品一区,二区,三区| 国产精品成人在线| av在线app专区| 国产亚洲午夜精品一区二区久久| 亚洲成av片中文字幕在线观看| 少妇猛男粗大的猛烈进出视频| 亚洲九九香蕉| 国产黄色视频一区二区在线观看| 在线观看免费午夜福利视频| 少妇精品久久久久久久| 伦理电影免费视频| videosex国产| 777久久人妻少妇嫩草av网站| 男人操女人黄网站| av在线app专区| 欧美另类一区| videosex国产| 日韩 亚洲 欧美在线| 如日韩欧美国产精品一区二区三区| 亚洲精品一区蜜桃| www.熟女人妻精品国产| 成人免费观看视频高清| 下体分泌物呈黄色| 99香蕉大伊视频| 天天影视国产精品| 欧美日韩视频高清一区二区三区二| 久久精品人人爽人人爽视色| 水蜜桃什么品种好| 视频区图区小说| 日韩电影二区| 亚洲精品av麻豆狂野| 亚洲av日韩在线播放| 亚洲中文字幕日韩| 国产成人免费观看mmmm| 宅男免费午夜| 美国免费a级毛片| 久久免费观看电影| 久久人人97超碰香蕉20202| 大片电影免费在线观看免费| 99re6热这里在线精品视频| 亚洲精品国产av成人精品| 国产精品国产三级专区第一集| 99久久人妻综合| 欧美另类一区| 免费高清在线观看日韩| 一级黄片播放器| 视频区图区小说| 日韩欧美一区视频在线观看| 午夜福利视频精品| 久久免费观看电影| 99久久人妻综合| 男女高潮啪啪啪动态图| 午夜福利在线免费观看网站| 久热爱精品视频在线9| 一区二区av电影网| 久久久欧美国产精品| 色播在线永久视频| 亚洲av国产av综合av卡| 97在线人人人人妻| 国产日韩一区二区三区精品不卡| 久久中文字幕一级| 久久久久久免费高清国产稀缺| 国产精品亚洲av一区麻豆| kizo精华| 无限看片的www在线观看| 国产精品成人在线| 国产一区有黄有色的免费视频| 久久影院123| 国产有黄有色有爽视频| 欧美精品高潮呻吟av久久| 美女扒开内裤让男人捅视频| 看十八女毛片水多多多| 国产精品99久久99久久久不卡| 别揉我奶头~嗯~啊~动态视频 | 国产精品久久久久久人妻精品电影 | 国产精品一区二区在线不卡| 精品国产一区二区三区久久久樱花| 亚洲成人免费电影在线观看 | 宅男免费午夜| 夜夜骑夜夜射夜夜干| 免费人妻精品一区二区三区视频| av在线老鸭窝| 久久精品久久精品一区二区三区| 交换朋友夫妻互换小说| www.自偷自拍.com| 成人影院久久| 狠狠婷婷综合久久久久久88av| 国产精品久久久久久精品古装| 91麻豆av在线| 女性生殖器流出的白浆| 亚洲欧美日韩另类电影网站| 在线av久久热| 国产在线视频一区二区| 水蜜桃什么品种好| 极品人妻少妇av视频| 黄片小视频在线播放| 超碰成人久久| 久9热在线精品视频| 欧美黑人精品巨大| 午夜91福利影院| 成年人黄色毛片网站| 日韩电影二区| 王馨瑶露胸无遮挡在线观看| 久久精品亚洲av国产电影网| 欧美97在线视频| 午夜福利视频精品| 嫁个100分男人电影在线观看 | 精品一区二区三区av网在线观看 | av国产久精品久网站免费入址| 男女下面插进去视频免费观看| 热re99久久精品国产66热6| 一级毛片 在线播放| 免费看不卡的av| 中文字幕高清在线视频| 精品福利永久在线观看| 国产精品一二三区在线看| 一区二区av电影网| 国产亚洲av片在线观看秒播厂| 首页视频小说图片口味搜索 | 欧美中文综合在线视频| 人妻一区二区av| 91精品伊人久久大香线蕉| 免费观看人在逋| 九色亚洲精品在线播放| 制服诱惑二区| 久热爱精品视频在线9| 国产熟女午夜一区二区三区| 亚洲伊人色综图| 精品一区二区三区四区五区乱码 | 婷婷色麻豆天堂久久| 久久久久久久大尺度免费视频| 少妇的丰满在线观看| 亚洲精品久久成人aⅴ小说| 午夜免费鲁丝| 美女脱内裤让男人舔精品视频| 欧美精品亚洲一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 丁香六月天网| √禁漫天堂资源中文www|