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

    遼寧省玉米地水分盈虧時(shí)空分布特征及灌溉模式分區(qū)研究

    2018-11-23 06:59:54魏新光王鐵良劉守陽(yáng)姚名澤鄭思宇景竹然
    關(guān)鍵詞:盈虧需水量生育期

    魏新光,王鐵良※,李 波,劉守陽(yáng),姚名澤,解 影,鄭思宇,景竹然

    ?

    遼寧省玉米地水分盈虧時(shí)空分布特征及灌溉模式分區(qū)研究

    魏新光1,王鐵良1※,李 波1,劉守陽(yáng)2,姚名澤1,解 影1,鄭思宇1,景竹然1

    (1. 沈陽(yáng)農(nóng)業(yè)大學(xué)水利學(xué)院,沈陽(yáng) 110866; 2. 法國(guó)國(guó)家農(nóng)業(yè)科學(xué)研究院地中海氣候與農(nóng)業(yè)水文過(guò)程模擬實(shí)驗(yàn)室,阿維尼翁 453003)

    根據(jù)不同區(qū)域農(nóng)作物需水規(guī)律和水資源供需狀況,對(duì)區(qū)域內(nèi)干旱區(qū)域類型、水分盈虧特性與干旱發(fā)生頻率等進(jìn)行綜合分析和科學(xué)區(qū)劃,對(duì)于提高整個(gè)地區(qū)(流域)農(nóng)作物水分利用效率和水資源高效利用具有重要意義。該文基于遼寧省27個(gè)氣象站1955—2014年逐日氣象數(shù)據(jù)和玉米生長(zhǎng)發(fā)育資料,對(duì)遼寧省不同玉米種植區(qū)域玉米生育期的需水過(guò)程、需水量、灌溉需水量、水分盈虧指數(shù)(crop water surplus deficit index,CWSDI)、干旱發(fā)生頻率等的時(shí)空分布特征進(jìn)行深入研究,得到以下結(jié)論:遼寧省不同區(qū)域玉米逐月需水量均呈現(xiàn)單峰變化趨勢(shì),7月需水量最大,全省生育期總需水量在335~391 mm之間;不同水文年玉米需水虧缺量在0~220 mm之間;CWSDI和干旱發(fā)生頻率在全省的空間分布規(guī)律類似,綜合兩指標(biāo)的數(shù)值分布特征,將遼寧省玉米灌區(qū)劃分為干旱區(qū)和易旱區(qū)2種類型,并結(jié)合不同水文年型玉米灌溉制度,將遼寧省玉米種植區(qū)劃分為7種灌溉模式。該研究成果可以為遼寧省區(qū)域農(nóng)業(yè)用水區(qū)劃與管理提供理論依據(jù)。

    需水量;灌溉;干旱;作物水分盈虧指數(shù);玉米

    0 引 言

    根據(jù)區(qū)域農(nóng)業(yè)布局、水資源分布狀況和氣象環(huán)境等特征,制定科學(xué)合理的灌溉制度和灌溉方式,對(duì)于提高區(qū)域水資源利用率、保障區(qū)域糧食安全具有重要意義[1]。玉米作為遼寧省第一大農(nóng)作物,在全省的13個(gè)地市均有分布。截止2015年,全省玉米種植面積達(dá)329.7×104hm2,產(chǎn)量2202.2×104t,約占全國(guó)總播種面積和總產(chǎn)量的8.6%和9.0%[2]。玉米作為遼寧省最主要的糧食作物和農(nóng)業(yè)用水大戶,對(duì)其水分供需特征、水量盈虧狀況進(jìn)行系統(tǒng)分析,并制定科學(xué)合理的灌溉制度與方式,進(jìn)而進(jìn)行分區(qū)灌溉勢(shì)在必行[3]。灌溉分區(qū)常用方法主要有:經(jīng)驗(yàn)法、指標(biāo)法、重疊發(fā)、聚類分析法等[4-6]。研究指標(biāo)主要有:作物參考蒸散量、作物需水量(evapotranspiration,ET)、灌溉需水量、作物水分盈虧指數(shù)(crop water surplus deficit index,CWSDI)、干旱發(fā)生頻率、濕潤(rùn)指數(shù)、干燥指數(shù)等[7-9]。在區(qū)域尺度上進(jìn)行灌溉分區(qū),前人有不少研究成果,康紹忠等[10]利用陜西省主要灌區(qū)灌溉資料和氣象數(shù)據(jù),確定了陜西省玉米需水量等值線圖,并在此基礎(chǔ)上,結(jié)合灌區(qū)灌溉制度,將陜西全省玉米種植區(qū)劃分為10個(gè)類型。高曉容等[11]對(duì)東北地區(qū)1961—2010年玉米生育階段的旱澇指數(shù)進(jìn)行分析,并基于旱澇發(fā)生頻率對(duì)東北地區(qū)玉米種植區(qū)域進(jìn)行分區(qū)。但這些區(qū)域特征和遼寧省差異較大。在全國(guó)尺度上,吳景社等[3]采用主成分和模糊聚類法,將全國(guó)節(jié)水灌區(qū)分為8大類型。仇寬彪等[12]利用遙感數(shù)據(jù)對(duì)中國(guó)蒸散比時(shí)空分布規(guī)律進(jìn)行研究,確定了不同時(shí)間段蒸散比較大的地區(qū),并確定影響該因子變化的主要驅(qū)動(dòng)因素。申雙和等[13]基1975—2004年長(zhǎng)時(shí)間氣象序列,對(duì)中國(guó)濕潤(rùn)指數(shù)時(shí)空分布狀況進(jìn)行分析,并利用該指數(shù)將中國(guó)濕潤(rùn)類型劃分為5個(gè)區(qū)域,但這些研究成果的空間尺度過(guò)大。遼寧省氣候類型復(fù)雜,水資源分布不均,物候條件差異較大。本文研究遼寧省玉米水分盈虧時(shí)空分布特征并對(duì)全省不同區(qū)域的灌溉模式進(jìn)行分區(qū),旨在為提高全省的玉米產(chǎn)量和水分利用效率、促進(jìn)區(qū)域水資源高效利用、保障遼寧省玉米產(chǎn)業(yè)的健康發(fā)展提供依據(jù)。

    1 材料與方法

    1.1 數(shù)據(jù)來(lái)源

    本研究選取遼寧省13個(gè)市的27個(gè)氣象站為代表站點(diǎn)(圖1),研究所用數(shù)據(jù)有:1)各站點(diǎn)經(jīng)緯度、高程信息;2)1955—2014年逐日氣象要素(最高氣溫、最低氣溫、平均氣溫、風(fēng)速、降水量、日照時(shí)數(shù)、水氣壓),基本氣象數(shù)據(jù)來(lái)源于中國(guó)氣象局國(guó)家氣象信息中心資料室;3)遼寧省玉米物候期、生長(zhǎng)發(fā)育狀況,根據(jù)文獻(xiàn)[14]確定。

    圖1 研究區(qū)域與氣象站點(diǎn)分布

    1.2 玉米需水量的確定

    玉米生育階段和全生育期需水量(evapotranspiration,ET)采用式(1)和式(2)確定:

    ET=K·ET0i(1)

    式中ETi為玉米第生育階段需水量,mm;ET為全生育期需水量,mm;K為第生育階段作物系數(shù);ET0i為第生育階段參考作物蒸發(fā)蒸騰量。玉米不同生育階段ET的確定采用FAO-56推薦的P-M公式[15]進(jìn)行確定,遼寧省不同區(qū)域分生育階段作物系數(shù)K采用李晶等[16]研究成果確定。

    1.3 玉米生育期灌溉需水量的確定

    玉米生育期灌溉需水量采用水量平衡方程[17]進(jìn)行計(jì)算。由于玉米種植區(qū)域?yàn)楹堤铮^大部分區(qū)域地下水埋深大于3 m,基本不考慮地下水補(bǔ)給。此外,東北地區(qū)玉米非生育期氣溫較低,大部分時(shí)期屬于土壤封凍時(shí)期,土壤水分狀況較為穩(wěn)定,春季解凍期積雪融化與土壤蒸發(fā)會(huì)對(duì)土壤水分狀況造成一定影響,但尚缺乏較為完整的全省尺度非生育期土壤水分盈虧數(shù)據(jù)。高歌等[18]的研究發(fā)現(xiàn),遼河流域冬春季月土壤水分虧缺量在 0~10 mm之間。張淑杰等[19]發(fā)現(xiàn)東北地區(qū)整個(gè)冬季水分虧缺量在10~40 mm之間,由此可見(jiàn)東北地區(qū)土壤水分虧缺量并不大,同時(shí)由于研究區(qū)域類型復(fù)雜,非生育期盈虧量空間差異明顯,且不易監(jiān)測(cè)。為研究方便,非生育期土壤水分盈虧量近似取0,這會(huì)使得玉米需灌水量在一定程度被低估,最多低估40 mm。生育期有效降水量0采用式(3)確定。

    0=·P(3)

    式中P為降水量,mm;為降水入滲系數(shù),受土壤性質(zhì)、地面覆蓋、地形、一次降水量、降水強(qiáng)度及降水延續(xù)時(shí)間等因素影響。一次降水量<5 mm時(shí),為0;當(dāng)一次降水時(shí)在5~50 mm時(shí),約為0.9;當(dāng)次降水量>50 mm時(shí),為0.8[20]。

    1.4 玉米水分盈虧指數(shù)與旱澇等級(jí)的確定

    作物水分虧缺指數(shù)[21]是以生育期需水量為需水指標(biāo),以有效降雨量為供水指標(biāo),來(lái)表征水分虧缺程度?;谠撝笖?shù)構(gòu)建玉米生育期水分盈虧指數(shù)(crop water surplus deficit index,CWSDI):

    CWSDI=(0-ET)/ET (4)

    當(dāng)CWSDI>0時(shí),表示玉米生育期水分盈余,當(dāng)CWSDI=0時(shí),表示玉米生育期水分供需平衡,當(dāng)CWSDI<0時(shí),表示玉米生育期水分虧缺。由于水分盈虧指數(shù)能夠真實(shí)反映作物供水量和需水量的平衡關(guān)系,可以較好地表征農(nóng)田的濕潤(rùn)程度與旱澇狀況[21]。參考前人基于水分虧缺指數(shù)的旱澇等級(jí)劃分[11,21],并結(jié)合遼寧玉米生育階段災(zāi)情統(tǒng)計(jì)資料[22],確定遼寧玉米生育階段玉米水分盈虧指數(shù)的旱澇等級(jí),如表1所示。

    表1 基于水分盈虧指數(shù)的遼寧省玉米旱澇等級(jí)劃分標(biāo)準(zhǔn)

    1.5 數(shù)據(jù)分析方法

    采用SPSS13.5進(jìn)行數(shù)據(jù)處理與統(tǒng)計(jì)分析,采用Excel2007、ArcGis10.2進(jìn)行相關(guān)圖形繪制。

    2 結(jié)果與分析

    2.1 玉米生育期需水量與水量盈虧

    2.1.1 玉米生育期需水規(guī)律

    為了較為全面地研究遼寧全省不同區(qū)域玉米生育期ET變化規(guī)律,在遼寧的東部、西部、南部、北部和中部分別選?。ǖ|、葉柏壽、大連、開(kāi)原和沈陽(yáng))5個(gè)代表站點(diǎn),對(duì)其逐月ET進(jìn)行分析。圖2為玉米生育期需水量逐月ET變化,由圖2可知,遼寧省不同區(qū)域ET變化規(guī)律基本一致,均呈單峰變化趨勢(shì)。在生育初期的5月份,不同區(qū)域ET值均較小,除遼南外,其他區(qū)域均無(wú)顯著差異(>0.05),而生育的中、后期(7—9月)遼西地區(qū)玉米ET均最大,特別是7月ET超過(guò)140 mm,其他地區(qū)則差異不顯著。

    2.1.2 玉米全生育期ET空間分布

    遼寧省不同地區(qū)全生育期ET變化如圖3所示,由圖3可知,全省ET變化范圍在335~391 mm, ET較大的地區(qū)主要集中在遼寧西部(ET≥380 mm),遼東的東部次之,遼寧北部、中部、南部和遼東大部需水量均較小,其中遼南南部最?。‥T<345 mm)。ET值波動(dòng)較大的地區(qū)集中在環(huán)渤海的遼寧中南部,其年際變化超15 mm,而遼寧東部和西部ET比較穩(wěn)定,年際變化一般在8 mm以下。

    注:同一月份不同小寫(xiě)字母表示區(qū)域間差異顯著(P<0.05)。

    圖3 遼寧省玉米全生育期需水量空間分布

    2.1.3 玉米全生育期水量盈虧

    為了準(zhǔn)確評(píng)估遼寧省玉米水分供需情況,按照生育期降水頻率將研究年份劃分為4種典型水文年:濕潤(rùn)年(=25%)、中水年(=50%)、干旱年(=75%)、特旱年(=90%)。圖4為遼寧省玉米全生育期水量盈虧空間分布,由圖4可知,不同水文年型下,玉米生育期水量盈虧狀況存在較大差異。在特旱年,僅遼寧東南部水量略有盈余,其他地區(qū)均出現(xiàn)水分虧缺,最大虧缺值達(dá)220 mm。隨著降水量的增加(水文年型由干旱向濕潤(rùn)變化),水分虧缺區(qū)域逐漸向遼寧西北部收縮。在濕潤(rùn)年,僅遼西西部出現(xiàn)水分虧缺,虧缺量在0~50 mm之間,而其他地方的降水量均能滿足玉米生育期耗水需求。

    2.2 玉米生育期灌溉制度與灌溉模式。

    遼寧省玉米一般4月底、5月初開(kāi)始播種、9月下旬開(kāi)始收獲,全生育期約為140~150 d,玉米生育期發(fā)育進(jìn)程較快,對(duì)水分需求緊迫,短期內(nèi)的干旱就可能造成難于挽回的損失,尤其是拔節(jié)、孕穗和開(kāi)花階段,短期的水分虧缺就會(huì)造成嚴(yán)重的減產(chǎn)甚至絕收[10,23]。玉米的播種期土壤水分狀況也對(duì)水玉米產(chǎn)量形成比較明顯的影響[24]。所以這幾個(gè)生育期是玉米灌溉需水的關(guān)鍵期??到B忠等[10]認(rèn)為,當(dāng)灌溉水補(bǔ)給不足時(shí),應(yīng)以保種(播種期)、?;ǎㄩ_(kāi)花期)和蓄墑灌溉為重點(diǎn)。根據(jù)遼寧省玉米灌溉需水規(guī)模在0~220 mm的現(xiàn)狀(圖4),考慮到遼寧省玉米水分虧缺嚴(yán)重的地區(qū)主要集中于遼西,水源取水緊張的特點(diǎn),參考遼寧省玉米水分敏感期相關(guān)研究成果[25],以及康紹忠等[10]春玉米非充分灌溉制度研究成果,補(bǔ)水灌水規(guī)??刂圃?0~135 mm之間,次灌水定額30~40 mm之間,保證玉米在極端干旱情況下(特旱年),灌水量135 mm(灌溉保證率為0.65)。在水源地比較豐富的地區(qū),可以采用充分灌溉,當(dāng)水源不足時(shí)可以采用非充分灌溉(灌溉保證率≥0.65)。據(jù)此,制定出遼寧省不同水文年玉米灌溉制度,如表2所示,4種水文年型下的灌溉定額分別為40、80、120和135 mm。

    注:Pp為水文頻率,下同。

    表2 不同水文年遼寧省推薦玉米灌溉制度

    根據(jù)玉米生育期水量盈虧研究成果(圖4),以及玉米灌溉制度(表2),確定了不同水文年遼寧省不同區(qū)域的灌溉制度空間分布如圖5所示。由圖5可知,在濕潤(rùn)年(圖5a),僅遼西的西部可采用A類灌溉制度。其他地方均無(wú)需灌溉;在中水年(圖5b),遼西大部地區(qū)均需進(jìn)行灌溉,其中遼西的西部需要采用B類灌溉制度,其余地區(qū)可采用A類灌溉制度;在干旱年(圖5c),遼寧中東部、遼東和遼南東部無(wú)需進(jìn)行灌溉,其余地區(qū)均應(yīng)進(jìn)行灌溉,其中遼西西部和北部應(yīng)采取類C灌溉制度;在特旱年(圖5b)僅遼東東部無(wú)需灌溉,其他地區(qū)均存在不同程度缺水,需要進(jìn)行灌溉補(bǔ)水,其中在遼西的西部和北部地區(qū)還應(yīng)采用D類灌溉制度。

    圖5 遼寧省不同水文年玉米灌溉制度分區(qū)

    2.3 基于作物CWSDI的玉米干旱類型區(qū)劃

    2.3.1 玉米CWSDI的時(shí)間變化

    根據(jù)式(4)得到遼寧省不同區(qū)域CWSDI年際變化如圖6所示,由圖6可知,不同地區(qū)CWSDI變化規(guī)律差異較大,遼東地區(qū)CWSDI普遍>0,僅在個(gè)別年份<0,且年際變化劇烈。其中,1985年其值超過(guò)2.5,水分供給遠(yuǎn)大于需求,說(shuō)明本區(qū)域水資源較為充沛,僅在個(gè)別年份出現(xiàn)供水量不足,基本不出現(xiàn)水量虧缺。遼西地區(qū)和遼東地區(qū)差異較大,不僅普遍CWSDI <0,而且大部分年份CWSDI<-0.30,個(gè)別年份CWSDI<-0.45。說(shuō)明遼西地區(qū)輕度和中度干旱均較為普遍,重度干旱也時(shí)有發(fā)生。

    注:CWSDI表示作物水分盈虧指數(shù),下同。

    2.3.2 玉米CWSDI的空間分布

    將全省各典型站點(diǎn)1955—2014年CWSDI數(shù)據(jù)進(jìn)行多年平均,并計(jì)算其標(biāo)準(zhǔn)偏差,得到其空間分布與變異狀況如圖7所示。由圖7可知,遼西大部分地區(qū)CWSDI<0,遼西西部甚至<–0.14,由此可見(jiàn)該區(qū)大部分普遍存在缺水,西部局地比較嚴(yán)重。但是由于該值年際波動(dòng)(圖7b)較為劇烈(全省普遍>0.48,遼東東部甚至>0.54),也就是說(shuō)該等值線年際漂移范圍過(guò)大,單純根據(jù)多年平均的CWSDI值進(jìn)行干旱區(qū)域劃分誤差較大。為此筆者對(duì)不同類型干旱發(fā)生頻率進(jìn)行進(jìn)一步分析。

    圖7 遼寧省玉米水分盈虧指數(shù)空間分布

    2.3.3 干旱發(fā)生頻率空間分布

    根據(jù)表1旱澇等級(jí)標(biāo)準(zhǔn)確定的旱澇等級(jí),選取輕旱、中旱和重旱3種干旱等級(jí)進(jìn)行重點(diǎn)分析。輕旱發(fā)生頻率[26],即為總研究序列(1955—2014年共60 a)中,輕旱發(fā)生次數(shù)與總研究年數(shù)的比值。將計(jì)算得到的頻率值用反距離加權(quán)插值法進(jìn)行空間插值,即可得到本區(qū)域輕旱、中旱、重旱發(fā)生頻率的空間分布,如圖8a所示。中旱和重旱發(fā)生頻率計(jì)算方法亦然,得到空間分布如圖8b和8c所示。從輕旱發(fā)生頻率的空間分布來(lái)看(圖8a)除遼東南一帶不發(fā)生輕度干旱,其他大部分地區(qū)(清源、岫巖、莊河一線以西地區(qū))均會(huì)發(fā)生輕度干旱,特別是遼西大部、遼北的大部分地區(qū),輕度干旱的頻率超過(guò)0.5,大約2a發(fā)生1次小旱;發(fā)生中度干旱的地區(qū)(圖8b)主要集中于遼寧西部,北部和遼中與遼南的西部,西豐、沈陽(yáng)、瓦房店一線以西區(qū)域,由東南向西北發(fā)生頻率依次增加。其中遼西和遼西北地區(qū)中旱發(fā)生頻率>0.3,該地大約3a 發(fā)生1次中旱;發(fā)生重度干旱的地區(qū)主要集中于遼寧西部和西北部,彰武,義縣、建昌一線以西的地區(qū)。該區(qū)域大部分地區(qū)中旱頻率超過(guò)0.1,總體而言,遼寧西部和西北部的大部分地區(qū)2a發(fā)生1次小旱,3a發(fā)生1次中旱,10a發(fā)生1次大旱。

    綜合圖7a、圖8可知,CWSDI值和干旱發(fā)生頻率空間分布規(guī)律基本類似。且輕度干旱發(fā)生頻率≥0.45,中度干旱發(fā)生頻率≥0.28,重度干旱發(fā)生頻率≥0.12和CWSDI≤0的地區(qū)基本重合,所以筆者將滿足以上條件之一的地區(qū)均定義為干旱區(qū)(圖9),該區(qū)主要集中于遼寧西部和西北部,水資源整體供小于求,而且經(jīng)常出現(xiàn)中度干旱。從圖7b和圖8可以看出,輕度干旱發(fā)生頻率≥0.1的區(qū)域和CWSDI≤0.45的地區(qū)也基本重合,將滿足2個(gè)條件之一、且干旱區(qū)以外地區(qū)定義為易旱區(qū)(圖9),該區(qū)域主要集中于遼北、遼中大部以及遼東西部,區(qū)域水資源供需基本平衡,容易出現(xiàn)輕旱,但不會(huì)出現(xiàn)澇田現(xiàn)象。綜上所述,遼寧省玉米種植區(qū)共可劃分為干旱區(qū)、易旱區(qū)和不旱區(qū)3塊區(qū)域,其中需要進(jìn)行補(bǔ)充灌溉的是干旱區(qū)和易旱區(qū)2個(gè)區(qū)域。

    圖8 遼寧省玉米生育期干旱發(fā)生頻率空間分布

    圖9 遼寧省玉米種植區(qū)干旱類型劃分

    2.4 玉米生育期灌溉模式分區(qū)

    根據(jù)圖9可知,遼寧省干旱類型可以劃分為(干旱區(qū)和易旱區(qū))2種類型區(qū)域,而各區(qū)域均存在4種灌溉制度,因而理論上可以形成8種研究區(qū)域與灌溉制度的組合,即形成8種灌溉模式。但是灌溉制度的實(shí)際執(zhí)行,需要依據(jù)水文年型進(jìn)行劃分,由于不同水文年型下,灌溉制度的執(zhí)行區(qū)域存在變化(圖5),因而不同水文年型下,灌溉模式出現(xiàn)頻率不一致。在特旱年出現(xiàn)的灌溉模式最多,為6種,在濕潤(rùn)年出現(xiàn)的模式最少,僅有1種,即干旱區(qū)A灌溉制度。4種水文年型實(shí)際上共出現(xiàn)了7種灌溉模式,不同水文年型出現(xiàn)的灌溉模式如表3所示。結(jié)合不同水文年灌溉模式(表3)和干旱區(qū)劃分(圖9),得到了遼寧省玉米種植區(qū)灌溉模式的空間分布如圖10所示。

    表3 遼寧省不同水文年玉米灌溉模式

    圖10 遼寧省玉米灌溉模式分區(qū)

    由圖10可知,濕潤(rùn)年僅有1種灌溉模式I,分布在遼西西部;中水年有3種灌溉模式,分別為I、II和V,分布于遼西大部和遼北局部;干旱年出現(xiàn)4種灌溉模式:II、III、V、VI,覆蓋遼西全部地區(qū),和遼北、遼中、遼南的西部地區(qū);特旱年出現(xiàn)6種灌溉模式II、III、IV、V、VI、VII,廣泛分布于除遼東南以外的大部分地區(qū)。

    灌溉模式的劃分,不僅僅要考慮區(qū)域農(nóng)田水量平衡,而且要考慮到其所處的水源與地形條件,在灌溉水源比較豐富的地區(qū),應(yīng)采用充分灌溉,以最大限度提高單產(chǎn)。本研究中玉米的灌溉保證率在0.65~1之間,通過(guò)張淑杰等[25]對(duì)遼寧省玉米需水量與產(chǎn)量關(guān)系的研究表明,當(dāng)玉米充分供水時(shí),玉米產(chǎn)量可達(dá)到9.0′106~1.2′106kg/hm2,當(dāng)灌溉保證率下降到0.65時(shí),產(chǎn)量下降到0.8′106kg/hm2左右,特旱年的減產(chǎn)率在11%~33.3%之間。在灌溉水源不足的地區(qū),則應(yīng)采用非充分灌溉制度,優(yōu)先保證需水關(guān)鍵期用水,提高水分利用效率,盡量做到減水不減產(chǎn),或者減水少減產(chǎn)[26-27]。此外,本研究確定的遼寧省不同區(qū)域的灌溉模式主要考慮不同水文年型的水量供需平衡與灌溉制度2大因素,所得到的分區(qū)灌溉模式中,灌溉用水量均為田間凈灌溉需水量,并未考慮需配水形式和灌溉用水方式,各區(qū)域在對(duì)玉米進(jìn)行灌溉時(shí)還應(yīng)在本研究成果的基礎(chǔ)上,結(jié)合本區(qū)域地形、供水水源等特點(diǎn),進(jìn)行進(jìn)一步補(bǔ)充與完善。

    結(jié)合遼寧省實(shí)際情況,在遼西北干旱區(qū),應(yīng)建立永久性灌溉措施,大力發(fā)展高效節(jié)水農(nóng)業(yè)。在地表水資源相對(duì)豐富的大凌河[28]、白石水庫(kù)一帶優(yōu)先通過(guò)水庫(kù)、河流取水,進(jìn)行渠道輸水灌溉或者渠井結(jié)合灌溉。其他沒(méi)有引水條件的地區(qū)則以井灌取水為主,但是要根據(jù)當(dāng)?shù)厮Y源狀況,嚴(yán)格控制機(jī)井?dāng)?shù)量、布設(shè)密度、灌溉取水量,合理利用地下水資源,防止超采[29]。田間灌溉用水還應(yīng)與膜下滴灌等高效節(jié)水技術(shù)相結(jié)合,盡可能提高玉米的灌溉水利用效率。在遼寧省易旱區(qū),也有可能發(fā)生輕度水分虧缺,但干旱發(fā)生的頻率較低、程度較輕,虧缺時(shí)間較短,若建立永久灌溉設(shè)施,投資成本較大,且設(shè)備使用率較低,因此不建議建立永久灌溉設(shè)施,主要以開(kāi)展地膜覆蓋、秸稈還田等旱作節(jié)水保墑措施為主。適當(dāng)提高土壤水分利用效率,即可以滿足玉米用水需要。與此同時(shí),有條件的地區(qū)也可采用簡(jiǎn)易或可移動(dòng)式灌溉設(shè)施,在個(gè)別虧水年進(jìn)行補(bǔ)充灌溉。

    3 結(jié) 論

    本文根據(jù)遼寧省玉米分布狀況、生育期需水規(guī)律與降水補(bǔ)給資料,對(duì)遼寧省不同區(qū)域玉米生育期需水量時(shí)空分布特征,生育期水量盈虧狀況、灌溉制度、干旱發(fā)生頻率等綜合分析,得到以下4條研究結(jié)論:

    1)遼寧全省玉米不同生育階段需水量呈單峰變化趨勢(shì),在生育初期的5月份,不同區(qū)域需水量均較小,且差異不大。7—9月遼西地區(qū)耗水顯著高于其他地區(qū),各區(qū)域需水峰值均出現(xiàn)在7月,其中遼西地區(qū)的峰值超過(guò)140 mm。全生育期ET較大的地區(qū)主要集中在遼寧西部,全生育期ET在380 mm以上。遼寧中部、南部ET較小,其中遼南南部小于345 mm。

    2)遼寧不同水文年玉米生育期灌溉需水量盈虧情況存在較大差異,不同水文年玉米需水的虧缺量在0~220 mm之間波動(dòng)。綜合考慮區(qū)域水資源供需狀況,制定了4種水文年的灌溉制度,灌溉定額分別為40、80、120和135 mm,灌溉保證率≥0.65。

    3)根據(jù)水分盈虧指數(shù)(crop water surplus deficit index, CWSDI)和干旱發(fā)生頻率將遼寧省劃分為2個(gè)灌溉類型區(qū):干旱區(qū)和易旱區(qū),干旱區(qū)主要集中在遼西和遼西北,該地發(fā)生輕度干旱發(fā)生頻率≥0.45,中度干旱發(fā)生頻率≥0.28,重度干旱發(fā)生頻率≥0.12且CWSDI≤0。易旱區(qū)集中在遼北,遼中大部和遼東西部,該地區(qū)輕度干旱發(fā)生頻率≥0.1且CWSDI≤0.45。

    4)綜合考慮灌溉制度和干旱區(qū)域類型,將遼寧省玉米種植區(qū)分為7種灌溉模式:其中特旱年、干旱年、中水年和濕潤(rùn)年分別出現(xiàn)6種、4種、3種和1種。

    [1] 水利部農(nóng)村水利司. 全國(guó)節(jié)水灌溉“十五”發(fā)展計(jì)劃及 2010年規(guī)劃[R]. 北京:水利部,2001.

    [2] 國(guó)家統(tǒng)計(jì)局. 國(guó)家統(tǒng)計(jì)局關(guān)于2015年糧食產(chǎn)量的公告[R/OL]. 2015-12-08[2017-10-10]. http://www.stats.gov.cn/tjsj/ zxfb/201512/t20151208_1286449.html

    [3] 吳景社,康紹忠,王景雷,等. 基于主成分分析和模糊聚類方法的全國(guó)節(jié)水灌溉分區(qū)研究[J]. 農(nóng)業(yè)工程學(xué)報(bào),2004,20(4):64-68. Wu Jingshe, Kang Shaozhong, Wang Jinglei, et al. Zoning of water-saving irrigation in China using principal components analysis and fuzzy-C-means[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE),2004, 20(4): 64-68. (in Chinese with English abstract)

    [4] 殷長(zhǎng)琛,齊廣平,康燕霞. 基于云模型的甘肅省參考作物蒸散量時(shí)空分布特征[J]. 農(nóng)業(yè)工程學(xué)報(bào),2015,31(8):152-158. Yin Changchen, Qi Guangping, Kang Yanxia. Analysis on characteristics of temporal-spatial potential evapotranspiration distribution in Gansu based on cloud mode[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(8): 152-158. (in Chinese with English abstract)

    [5] Tang Xuguang, Li Hengpeng, Desai A R, et al. How is water-use efficiency of terrestrial ecosystems distributed and changing on Earth[J]. Scientific Reports, 2014, 4: doi: 10.1038/srep07483

    [6] 穆少杰,周可新,齊楊,等. 內(nèi)蒙古植被降水利用效率的時(shí)空格局及其驅(qū)動(dòng)因素[J]. 植物生態(tài)學(xué)報(bào),2014,38(1):1-16. Mu Shaojie, Zhou Kexin, Qi Yang, et al. Spatio-temporal patterns of precipitation-use efficiency of vegetation and their controlling factors in Inner Mongolia[J]. Chinese Journal of Plant Ecology, 2014, 38(1): 1-16. (in Chinese with English abstract)

    [7] 張玉芳,王明田,劉娟,等. 基于水分盈虧指數(shù)的四川省玉米生育期干旱時(shí)空變化特征分析[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào),2013,21(2):236-242. Zhang Yufang, Wang Mingtian, Liu Juan, et al. Spatio-temporal characteristics of drought at different maize growth stages in Sichuan province as determined by water budget index[J]. Chinese Journal of Eco-Agriculture, 2013, 21(2): 236-242. (in Chinese with English abstract)

    [8] 王勁松,郭江勇,周躍武,等. 干旱指標(biāo)研究的進(jìn)展與展望[J]. 干旱區(qū)地理,2007,30(1):60-65. Wang Jinsong, Guo Jiangyong, Zhou Yuewu, et al. Progress and prospect on drought indices research[J]. Arid Land Geography, 2007, 30(1): 60-65. (in Chinese with English abstract)

    [9] Nutini F, Boschetti M, Candiani G, et al. Evaporative fraction as an indicator of moisture condition and water stress statusin semi-arid rangeland ecosystems[J]. Remote Sensing, 2014, 6(7): 6300-6323.

    [10] 康紹忠,賀正中,張學(xué). 陜西省作物需水量及分區(qū)灌溉模式[M ]. 北京:水利水電出版社,1992.

    [11] 高曉容,王春乙,張繼權(quán),等. 近50年?yáng)|北玉米生育階段需水量及旱澇時(shí)空變化[J]. 農(nóng)業(yè)工程學(xué)報(bào),2012,28(12):101-109. Gao Xiaorong, Wang Chunyi, Zhang Jiquan, et al. Crop water requirement and temporal-spatial variation of drought and flood disaster during growth stages for maize in Northeast during past 50 years[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(12): 101-109. (in Chinese with English abstract)

    [12] 仇寬彪,賈寶全,張志強(qiáng). 基于遙感的中國(guó)植被蒸散比時(shí)空分布及其對(duì)氣候因素的響應(yīng)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2015,31(19):151-158. Qiu Kuanbiao, Jia Baoquan, Zhang Zhiqiang. Spatio- temporal patterns of evaporative fraction in vegetation ecosystem based onremote sensing data and its responses to climatic factors in China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(19): 151-158. (in Chinese with English abstract)

    [13] 申雙和,張方敏,盛瓊. 1975-2004 年中國(guó)濕潤(rùn)指數(shù)時(shí)空變化特征[J]. 農(nóng)業(yè)工程學(xué)報(bào),2009,25(1):11-15. Shen Shuanghe, Zhang Fangmin, Sheng Qiong. Spatio- temporal changes of wetness index in China from 1975 to 2004[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2009, 25(1): 11-15. (in Chinese with English abstract)

    [14] 李正國(guó),楊鵬,唐華俊,等. 氣候變化背景下東北三省主要作物典型物候期變化趨勢(shì)分析[J]. 中國(guó)農(nóng)業(yè)科學(xué),2011,44(20):4180-4189. Li Zhengguo, Yang Peng, Tang Huajun, et al. Trend analysis of typical phenophases of major crops under climate change in the three provinces of Northeast China[J]. Scientia Agricultura Sinica, 2011, 44(20): 4180-4189. (in Chinese with English Abstract)

    [15] Allen R G, Pereira L S, Raes D, et al. Crop Evapotranspiration Guidelines for Computing Crop Water Requirements[M]. Rome: Food and Agriculture Organization of United Nation, 1998.

    [16] 李晶,張運(yùn)福,班顯秀. 遼寧玉米作物系數(shù)研究[J]. 遼寧氣象,2000(1):18-20.

    [17] 尚松浩. 作物非充分灌溉制度的模擬優(yōu)化方法[J]. 清華大學(xué)學(xué)報(bào):自然科學(xué)版,2005,45(9):1179-1183. Shang Songhao. Simulation-optimization method for crop irrigation scheduling with limited water supplies [J]. Journal of Tsinghua University: Science and Technology, 2005, 45(9): 1179-1183. (in Chinese with English Abstract)

    [18] 高歌,許崇育. 1961-2010年中國(guó)十大流域水分盈虧量時(shí)空變化特征[J]. 地理學(xué)報(bào),2015,70(3):380-391. Gao Ge, Xu Congyu. Characteristics of water surplus and deficit change in 10 major river basins in China during 1961-2010[J]. Acta Geographica Sinica, 2015, 70(3): 380-391. (in Chinese with English Abstract)

    [19] 張淑杰,張玉書(shū),孫龍彧,等. 東北地區(qū)玉米生育期干旱分布特征及其成因分析[J]. 中國(guó)農(nóng)業(yè)氣象,2013,34(3):350-357. Zhang Shujie, Zhang Yushu, Sun Longyu, et al. Analysis of distributional characteristics and primary causes of maize drought in northeast china[J]. Chinese Joural Agrometeorology, 2013, 34(3): 350-357. (in Chinese with English abstract)

    [20] 龐艷梅,陳超,潘學(xué)標(biāo). 1961-2010 年四川盆地玉米有效降水和需水量的變化特征[J]. 農(nóng)業(yè)工程學(xué)報(bào),2015,31(增刊1):133-141. Pang Yanmei, Chen Chao, Pan Xuebiao. Variation characteristics of maize effective precipitation and water requirement in Sichuan basin during 1961-2010[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(Supp. 1): 133-141. (in Chinese with English abstract)

    [21] 張艷紅,呂厚荃,李森. 作物水分虧缺指數(shù)在農(nóng)業(yè)干旱監(jiān)測(cè)中的適用性[J]. 氣象科技,2008,36(5):596-600. Zhang Yanhong, Lü Houquan, Li Sen. Applicability of crop water deficit index in agricultural drought monitoring[J]. Meteorological Science and Technology, 2008, 36(5): 596-600. (in Chinese with English abstract)

    [22] 張淑杰,張玉書(shū),李廣霞,等. 2005年氣象災(zāi)害對(duì)遼寧農(nóng)業(yè)生產(chǎn)影響評(píng)估[J]. 氣象與環(huán)境學(xué)報(bào),2006,22(6):46-49. Zhang Shujie, Zhang Yushu , Li Guangxia , et al. Evaluating of the effects of meteorological disaster on Liaoning agricultural production in 2005[J]. Journal of Meteorology and Environment 2006, 22(6): 46-49. (in Chinese with English abstract)

    [23] 劉鈺,汪林,倪廣恒,等. 中國(guó)主要作物灌溉需水量空間分布特征[J]. 農(nóng)業(yè)工程學(xué)報(bào),2009,25(12):6-12. Liu Yu, Wang Lin, Ni Guangheng, et al. Spatial distribution characteristics of irrigation water requirement for main crops inChina[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2009, 25(12): 6-12. (in Chinese with English abstract)

    [24] 王智威. 水分脅迫下玉米苗期的生長(zhǎng)及生理響應(yīng)[D]. 楊凌:西北農(nóng)林科技大學(xué),2013.

    Wang Zhiwei. Growth and Physiology of Mazie Seedling Response under Water Stress[D]. Yangling: Northwest A&F University, 2013.

    [25] 張淑杰,周廣勝,張玉書(shū),等. 遼寧地區(qū)玉米耗水量與產(chǎn)量的關(guān)系[J]. 干旱地區(qū)農(nóng)業(yè)研究,2018,36(1):244-250. Zhang Shujie, Zhou Guangsheng, Zhang Yushu, et al. The relationship between yield and water consumption of maize in Liaoning[J]. Agricultural Research in the Arid Areas, 2018, 36(1): 244-250 (in Chinese with English abstract)

    [26] 黃晚華,楊曉光,李茂松,等. 基于標(biāo)準(zhǔn)化降水指數(shù)的中國(guó)南方季節(jié)性干旱近58 a 演變特征[J]. 農(nóng)業(yè)工程學(xué)報(bào),2010,26(7):50-59. Huang Wanhua, Yang Xiaoguang, Li Maosong, et al. Evolution characteristics of seasonal drought in the south of China during the past 58 years based on standardized precipitation index[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(7): 50-59. (in Chinese with English abstract)

    [27] 康紹忠,史文娟,胡笑濤,等. 調(diào)虧灌溉對(duì)于玉米生理指標(biāo)及水分利用效率的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),1998,14(4):82-87. Kang Shaozhong, Shi Wenjuan, Hu Xiaotao, et al. Effects of regulated deficit irrigation on physiological indices and water use efficiency of maize[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 1998, 14(4): 82-87. (in Chinese with English abstract)

    [28] 熊偉,馮穎竹,高清竹,等. 氣候變化對(duì)石羊河、大凌河流域灌溉玉米生產(chǎn)的影響[J]. 干旱區(qū)地理,2011,34(1):150-159. Xiong Wei, Feng Yingzu, Gao Qingzhu, et al. Impacts of climate change on irrigated maize production in Shiyang and Daling River Basins[J] Arid Land Geography, 2011, 34(1): 150-159. (in Chinese with English abstract)

    [29] 齊悅. 大型水庫(kù)生態(tài)效應(yīng)研究[D]. 長(zhǎng)春:吉林大學(xué),2011. Qi Yue. Study on Ecological Effect of Large Reserior[D]. Changchun: Jilin University, 2011. (in Chinese with English abstract)

    Temporal and spatial distribution characteristics of maize water surplus deficit and irrigation mode partition in Liaoning province

    Wei Xinguang1, Wang Tieliang1※, Li Bo1, Liu Shouyang2, Yao Mingze1, Xie Ying1, Zheng Siyu1, Jing Zhuran1

    (1.110866,; 2.453003,)

    It is necessary to formulate a partition irrigation management mode based on the feature of water resource and meteorology of different sub areas in a large area. Maize is the main crop and major agricultural water consumer in Liaoning province. The objective of this study was to analyze the temporal and spatial distribution characteristics of maize water surplus deficit and lay out an irrigation mode partition scheme in Liaoning province. We collected the data of maize growth index, irrigation, and daily meteorological data from 1955 to 2014 in 27 agricultural meteorological stations in Liaoning province. The calculated index included evapotranspiration (ET), irrigation requirement, crop water surplus and deficit index (CWSDI)?the drought probability during the maize growth stage in Liaoning. The results showed that: the ET during the growth stage of maize varied in a single peak style in all the regions of Liaoning. The ET in western Liaoning was the highest in Liaoning, and the monthly ET of the western Liaoning was more than 140 mm in July. The high value of ET mainly appeared in the western region (ET≥380 mm). The ET was lower in the central and southern regions in Liaoning. The minimum value of ET was in the southern Liaoning, which was less than 345 mm. The standard deviation of ET were higher in the central and southern Liaoning around the Bohai sea (SD>15 mm), but it were stable in eastern and western Liaoning (SD<8 mm). A total of 4 typically deficient irrigation schedule were recommended and the irrigation norm in different hydrological years were 40, 80, 120 and 135 mm, respectively for the wet year, the normal year, the dry year and the drought year. In most years, CWSDI of maize were above 0 in the eastern Liaoning, and the peak was above 2.5 (in the year of 1985). It indicates that the crop water supply exceeded demand in most of years. It was significantly different between the western and eastern Liaoning. In the western Liaoning, the most of the CWSDI was below 0, and sometimes even below-0.14. There was a widespread mild dry in most of the years, sometimes was dry even drought in this region. The spatial distribution of the dryfrequency indicated that, in addition to the southeast Liaoning, the most regions may become mild dry, dry or drought, most especially in the western Liaoning. In the northwest Liaoning, the frequency of mild drought occurrence was high. The dry mainly appeared in the western Liaoning. The CWSDI and the dry frequency had the similar spatial distribution pattern. The area with mild dry frequency not less than 0.45, dry frequency not less than 0.28, the drought frequency not less than 0.12 and CWSDI higher than 0 was same. Thus, we defined the area meeting one of the conditions above as the arid area. Those areas were in the western and northern Liaoning where the water resource supply was less than demand and the dry condition appeared very often. We defined the area with mild dry frequency not less than 0.1 or CWSDI higher than 0.45 as the drouht-prone area. These areas were in the northern Liaoning, middle Liaoning and the west of eastern Liaoning. Based on the CWSDI and dry frequency, a total of 7 irrigation modes were set up. A total of 6, 4, 3 and 1 irrigation mode were for the drought, dry, normal and wet year, respectively. The study would improve not only the areas water resource utilization but also food security.

    evapotranspiration; irrigation; drought; CWSDI; maize

    魏新光,王鐵良,李 波,劉守陽(yáng),姚名澤,解 影,鄭思宇,景竹然. 遼寧省玉米地水分盈虧時(shí)空分布特征及灌溉模式分區(qū)研究[J]. 農(nóng)業(yè)工程學(xué)報(bào),2018,34(23):119-126. doi:10.11975/j.issn.1002-6819.2018.23.014 http://www.tcsae.org

    Wei Xinguang, Wang Tieliang, Li Bo, Liu Shouyang, Yao Mingze, Xie Ying, Zheng Siyu, Jing Zhuran. Temporal and spatial distribution characteristics of maize water surplus deficit and irrigation mode partition in Liaoning province[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(23): 119-126. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2018.23.014 http://www.tcsae.org

    2018-06-06

    2018-10-10

    國(guó)家自然科學(xué)基金項(xiàng)目(51709174);國(guó)家公益性行業(yè)(農(nóng)業(yè))科研專項(xiàng)(201303125);遼寧省博士科研啟動(dòng)基金(20170520169)

    魏新光,講師,博士,研究方向?yàn)樽魑锔咝в盟c蒸散發(fā)尺度效應(yīng)。Email:weixg_wi@163.com

    王鐵良,教授,博士,研究方向?yàn)楣?jié)水灌溉。 Email:tieliangwang@126.com

    10.11975/j.issn.1002-6819.2018.23.014

    S274.1

    A

    1002-6819(2018)-23-0119-08

    猜你喜歡
    盈虧需水量生育期
    南京城市中心區(qū)生態(tài)需水量探析
    巧用Excel盈虧圖分析差異一目了然
    兩種方法解盈虧問(wèn)題
    組合預(yù)測(cè)模型在區(qū)域需水量預(yù)測(cè)中的應(yīng)用
    帕滿灌區(qū)供需水量平衡分析
    大洋河生態(tài)需水量分析與評(píng)價(jià)
    基于作物生育期的潛在蒸散的時(shí)空演變特征及R/S 分析
    西藏科技(2015年5期)2015-09-26 11:55:25
    2013-2014年度二二二團(tuán)冬小麥各生育期氣象條件分析
    阿拉爾市2012年棉花生育期氣象條件分析
    棉花生育期調(diào)控關(guān)鍵技術(shù)
    日韩av免费高清视频| 男女高潮啪啪啪动态图| 夫妻性生交免费视频一级片| 久热这里只有精品99| 国产日韩欧美亚洲二区| 黄色 视频免费看| 黄网站色视频无遮挡免费观看| 咕卡用的链子| 一区二区三区四区激情视频| 丰满迷人的少妇在线观看| 777米奇影视久久| 久久久久国产一级毛片高清牌| 国产成人91sexporn| 一区二区av电影网| 在线观看免费日韩欧美大片| 国产福利在线免费观看视频| 国产高清不卡午夜福利| 黄色一级大片看看| 精品一品国产午夜福利视频| 国产野战对白在线观看| 男女午夜视频在线观看| 男女午夜视频在线观看| 成年av动漫网址| 久久久久久人妻| 狠狠精品人妻久久久久久综合| 久久久久久久国产电影| 久久久久久久大尺度免费视频| 欧美成人午夜精品| 一级a爱视频在线免费观看| 天天躁日日躁夜夜躁夜夜| 国产99久久九九免费精品| 人人妻人人添人人爽欧美一区卜| 精品一区在线观看国产| 国产亚洲最大av| 亚洲国产欧美在线一区| 欧美日韩亚洲国产一区二区在线观看 | 欧美日韩综合久久久久久| 美女午夜性视频免费| 日韩精品免费视频一区二区三区| 美女中出高潮动态图| 亚洲少妇的诱惑av| 欧美97在线视频| h视频一区二区三区| 在线观看免费视频网站a站| 考比视频在线观看| 亚洲精品一二三| 欧美激情 高清一区二区三区| 亚洲第一青青草原| 午夜福利,免费看| 午夜福利乱码中文字幕| 男人添女人高潮全过程视频| 精品国产一区二区三区四区第35| 国产精品久久久人人做人人爽| 男女床上黄色一级片免费看| 久久精品人人爽人人爽视色| 国产xxxxx性猛交| 欧美 日韩 精品 国产| 久久ye,这里只有精品| 精品一区在线观看国产| 久久久精品94久久精品| 激情五月婷婷亚洲| 国产片内射在线| 青春草亚洲视频在线观看| 黑人欧美特级aaaaaa片| 国产成人av激情在线播放| 啦啦啦视频在线资源免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 国产欧美日韩综合在线一区二区| 日本wwww免费看| 亚洲精品乱久久久久久| 高清av免费在线| 欧美精品av麻豆av| 国产精品一区二区在线不卡| 日本爱情动作片www.在线观看| 在线观看免费高清a一片| 精品国产国语对白av| 亚洲综合精品二区| 一本久久精品| 丝袜脚勾引网站| 秋霞在线观看毛片| 久久久久久免费高清国产稀缺| 自线自在国产av| 国产乱人偷精品视频| 日韩精品免费视频一区二区三区| 久久久欧美国产精品| 国产深夜福利视频在线观看| 精品亚洲成国产av| 日本91视频免费播放| 亚洲人成网站在线观看播放| 18禁国产床啪视频网站| 国产探花极品一区二区| 亚洲国产日韩一区二区| 国产成人精品在线电影| 亚洲七黄色美女视频| 丝袜人妻中文字幕| 欧美日韩福利视频一区二区| 日韩大片免费观看网站| 亚洲欧洲精品一区二区精品久久久 | kizo精华| 少妇被粗大的猛进出69影院| 久久精品国产亚洲av涩爱| 老司机深夜福利视频在线观看 | 男女无遮挡免费网站观看| 在线观看免费视频网站a站| 99re6热这里在线精品视频| 午夜福利在线免费观看网站| 亚洲国产精品一区二区三区在线| 丝袜人妻中文字幕| 在线观看免费午夜福利视频| 成人国语在线视频| 人人妻人人添人人爽欧美一区卜| 王馨瑶露胸无遮挡在线观看| 在线 av 中文字幕| 国产无遮挡羞羞视频在线观看| 精品亚洲成国产av| 熟女av电影| 国产av精品麻豆| 婷婷色麻豆天堂久久| 精品久久久久久电影网| 国产亚洲最大av| 国产在视频线精品| 国产精品久久久久久久久免| 国产色婷婷99| 90打野战视频偷拍视频| 国产精品久久久久久人妻精品电影 | 免费黄频网站在线观看国产| 久久这里只有精品19| 久久久久人妻精品一区果冻| 日日摸夜夜添夜夜爱| 丰满少妇做爰视频| 国产乱来视频区| 久久久国产欧美日韩av| 国产又色又爽无遮挡免| 岛国毛片在线播放| 男女免费视频国产| 别揉我奶头~嗯~啊~动态视频 | 在线观看免费视频网站a站| 欧美黑人欧美精品刺激| 夫妻性生交免费视频一级片| 亚洲激情五月婷婷啪啪| 欧美日韩视频精品一区| 国产熟女午夜一区二区三区| 国产成人欧美在线观看 | 大香蕉久久网| 亚洲第一区二区三区不卡| 汤姆久久久久久久影院中文字幕| 人人妻人人澡人人看| 久久久精品94久久精品| 欧美乱码精品一区二区三区| 黄色视频在线播放观看不卡| 青春草国产在线视频| 极品少妇高潮喷水抽搐| av视频免费观看在线观看| 男女床上黄色一级片免费看| 夫妻性生交免费视频一级片| 汤姆久久久久久久影院中文字幕| 视频区图区小说| 国产视频首页在线观看| 美女福利国产在线| 亚洲国产欧美网| 韩国av在线不卡| 久久精品国产a三级三级三级| 美女中出高潮动态图| 伦理电影大哥的女人| 丝袜人妻中文字幕| 免费不卡黄色视频| 欧美少妇被猛烈插入视频| 午夜福利视频精品| 亚洲一区中文字幕在线| 免费看av在线观看网站| 亚洲第一av免费看| 亚洲久久久国产精品| 欧美日韩亚洲综合一区二区三区_| 一级,二级,三级黄色视频| 欧美日韩视频精品一区| 各种免费的搞黄视频| 黑人欧美特级aaaaaa片| 免费日韩欧美在线观看| 如何舔出高潮| 日韩视频在线欧美| 大香蕉久久网| 在线观看免费日韩欧美大片| 亚洲综合色网址| 91精品伊人久久大香线蕉| 久久人人爽人人片av| av在线老鸭窝| 国产在线免费精品| 国产精品 国内视频| 亚洲伊人久久精品综合| 久久午夜综合久久蜜桃| 国产在视频线精品| 国产乱来视频区| 国产淫语在线视频| 国产在线免费精品| 欧美精品人与动牲交sv欧美| 国产成人免费观看mmmm| 91老司机精品| 丝袜脚勾引网站| 天天躁狠狠躁夜夜躁狠狠躁| 国产一区二区三区综合在线观看| 亚洲国产欧美一区二区综合| 精品国产一区二区久久| 国产视频首页在线观看| 热99国产精品久久久久久7| 纯流量卡能插随身wifi吗| 亚洲欧美成人精品一区二区| 性高湖久久久久久久久免费观看| 亚洲精品,欧美精品| 2018国产大陆天天弄谢| 国产免费一区二区三区四区乱码| 亚洲精品在线美女| 青春草视频在线免费观看| 久久女婷五月综合色啪小说| 两性夫妻黄色片| 国产精品人妻久久久影院| 国产成人免费无遮挡视频| 国产乱来视频区| 80岁老熟妇乱子伦牲交| 国产伦人伦偷精品视频| 亚洲av成人精品一二三区| 制服人妻中文乱码| 国产 精品1| videos熟女内射| 岛国毛片在线播放| 香蕉国产在线看| 亚洲国产中文字幕在线视频| 亚洲国产最新在线播放| 欧美日韩亚洲高清精品| 亚洲成人国产一区在线观看 | 国产精品.久久久| 国产精品成人在线| www.熟女人妻精品国产| 午夜日韩欧美国产| 亚洲成人一二三区av| 美女中出高潮动态图| 老司机影院毛片| 黄色怎么调成土黄色| 一本久久精品| 韩国精品一区二区三区| 欧美黄色片欧美黄色片| 久久久亚洲精品成人影院| 日韩大片免费观看网站| 丝袜美足系列| 晚上一个人看的免费电影| 熟女av电影| 日韩大码丰满熟妇| 国产亚洲午夜精品一区二区久久| 中文字幕高清在线视频| 老汉色∧v一级毛片| 午夜免费男女啪啪视频观看| 精品福利永久在线观看| 成人国产av品久久久| av网站在线播放免费| 亚洲精品国产av成人精品| 国产乱人偷精品视频| 捣出白浆h1v1| 久久性视频一级片| 久久热在线av| 超碰97精品在线观看| 日韩 欧美 亚洲 中文字幕| 欧美老熟妇乱子伦牲交| 亚洲一码二码三码区别大吗| 国产极品天堂在线| 日本猛色少妇xxxxx猛交久久| 亚洲四区av| 少妇精品久久久久久久| 在线天堂最新版资源| 精品一品国产午夜福利视频| 精品一品国产午夜福利视频| 人成视频在线观看免费观看| 在线观看三级黄色| av女优亚洲男人天堂| 日韩精品免费视频一区二区三区| 国产熟女欧美一区二区| 波多野结衣一区麻豆| 9色porny在线观看| 亚洲欧美一区二区三区国产| 亚洲视频免费观看视频| 天天操日日干夜夜撸| 亚洲欧洲精品一区二区精品久久久 | 欧美亚洲 丝袜 人妻 在线| 男女高潮啪啪啪动态图| 亚洲精品成人av观看孕妇| 亚洲,欧美精品.| 两性夫妻黄色片| 午夜免费鲁丝| 热re99久久精品国产66热6| 亚洲精品美女久久久久99蜜臀 | 一级,二级,三级黄色视频| 欧美另类一区| a 毛片基地| 1024香蕉在线观看| 婷婷色麻豆天堂久久| 久久综合国产亚洲精品| 日韩中文字幕欧美一区二区 | 美女中出高潮动态图| 免费观看av网站的网址| av又黄又爽大尺度在线免费看| 美女国产高潮福利片在线看| 两个人免费观看高清视频| 麻豆av在线久日| 亚洲国产精品一区三区| 亚洲三区欧美一区| 欧美日本中文国产一区发布| 午夜老司机福利片| 日韩成人av中文字幕在线观看| 婷婷成人精品国产| 日本黄色日本黄色录像| 国产一区二区在线观看av| 午夜日韩欧美国产| 亚洲国产毛片av蜜桃av| 日韩av不卡免费在线播放| 天堂8中文在线网| 三上悠亚av全集在线观看| 岛国毛片在线播放| 桃花免费在线播放| 国产精品久久久久成人av| 亚洲精品久久成人aⅴ小说| 纵有疾风起免费观看全集完整版| 亚洲欧洲日产国产| 国产精品99久久99久久久不卡 | 在线观看三级黄色| 一二三四在线观看免费中文在| 日韩,欧美,国产一区二区三区| 日本vs欧美在线观看视频| 中文乱码字字幕精品一区二区三区| 免费少妇av软件| 一级,二级,三级黄色视频| 久久99热这里只频精品6学生| 亚洲精品国产av蜜桃| 天堂中文最新版在线下载| 久久久精品94久久精品| 国产一级毛片在线| 午夜福利免费观看在线| 各种免费的搞黄视频| a 毛片基地| 欧美黑人欧美精品刺激| 国产国语露脸激情在线看| 高清欧美精品videossex| 国产精品香港三级国产av潘金莲 | 亚洲精品视频女| 人人澡人人妻人| 乱人伦中国视频| 大香蕉久久网| 老汉色av国产亚洲站长工具| 少妇猛男粗大的猛烈进出视频| 美女国产高潮福利片在线看| 久久久欧美国产精品| 国产精品无大码| 婷婷色综合www| 亚洲精品一二三| 视频区图区小说| 亚洲精品日本国产第一区| 捣出白浆h1v1| 亚洲伊人久久精品综合| 一本大道久久a久久精品| 日韩欧美一区视频在线观看| 亚洲熟女精品中文字幕| 精品久久久精品久久久| 日韩精品免费视频一区二区三区| 国产一区二区三区av在线| 王馨瑶露胸无遮挡在线观看| 精品一品国产午夜福利视频| 如日韩欧美国产精品一区二区三区| 亚洲国产欧美网| 巨乳人妻的诱惑在线观看| 免费观看a级毛片全部| 热99国产精品久久久久久7| 欧美老熟妇乱子伦牲交| 操美女的视频在线观看| 水蜜桃什么品种好| 不卡视频在线观看欧美| 又粗又硬又长又爽又黄的视频| 中文字幕亚洲精品专区| 中文天堂在线官网| 人人妻人人添人人爽欧美一区卜| 极品少妇高潮喷水抽搐| 十八禁人妻一区二区| 女人爽到高潮嗷嗷叫在线视频| 少妇被粗大的猛进出69影院| 日本色播在线视频| 超碰成人久久| 久久av网站| 一本色道久久久久久精品综合| 97精品久久久久久久久久精品| 国产男女超爽视频在线观看| 国产高清国产精品国产三级| 人人妻,人人澡人人爽秒播 | 日韩中文字幕欧美一区二区 | 亚洲精品aⅴ在线观看| 精品少妇久久久久久888优播| 亚洲av成人不卡在线观看播放网 | 亚洲国产成人一精品久久久| 香蕉丝袜av| 久久影院123| 天天躁狠狠躁夜夜躁狠狠躁| 国产极品天堂在线| 丁香六月欧美| 人妻 亚洲 视频| 国产一区二区在线观看av| 少妇 在线观看| 免费高清在线观看日韩| 国产精品三级大全| 亚洲精品一区蜜桃| 亚洲精品久久久久久婷婷小说| 日韩欧美精品免费久久| 国产不卡av网站在线观看| 亚洲av男天堂| 国产精品久久久久久精品古装| 欧美日韩亚洲综合一区二区三区_| 天堂8中文在线网| 日韩电影二区| 99香蕉大伊视频| 下体分泌物呈黄色| 久久精品人人爽人人爽视色| 欧美精品av麻豆av| 国产亚洲av片在线观看秒播厂| 国产有黄有色有爽视频| 狂野欧美激情性xxxx| 久久97久久精品| 少妇的丰满在线观看| xxx大片免费视频| 国产精品一区二区在线不卡| 国产高清国产精品国产三级| 欧美日韩视频高清一区二区三区二| 国语对白做爰xxxⅹ性视频网站| 国产精品久久久久久久久免| 波野结衣二区三区在线| 久久ye,这里只有精品| av在线app专区| 91aial.com中文字幕在线观看| 99久国产av精品国产电影| 看免费成人av毛片| 国产男女内射视频| 日韩制服丝袜自拍偷拍| 亚洲精品日韩在线中文字幕| 午夜福利影视在线免费观看| 久久av网站| 久久久国产一区二区| 久久精品国产a三级三级三级| 制服诱惑二区| 国产精品久久久久久精品电影小说| 亚洲成色77777| 国产高清国产精品国产三级| 国产亚洲av片在线观看秒播厂| 成人黄色视频免费在线看| 免费在线观看视频国产中文字幕亚洲 | 免费在线观看黄色视频的| 男人爽女人下面视频在线观看| 五月开心婷婷网| 一级毛片电影观看| 天天躁狠狠躁夜夜躁狠狠躁| 成人黄色视频免费在线看| 亚洲欧美一区二区三区国产| 欧美日韩亚洲综合一区二区三区_| 欧美 日韩 精品 国产| 女人久久www免费人成看片| 国产亚洲最大av| 国产精品久久久av美女十八| 欧美日韩一级在线毛片| 日韩一区二区三区影片| 满18在线观看网站| 少妇猛男粗大的猛烈进出视频| 免费在线观看完整版高清| 自拍欧美九色日韩亚洲蝌蚪91| 蜜桃国产av成人99| 色精品久久人妻99蜜桃| 2021少妇久久久久久久久久久| 黄色一级大片看看| 久久人人爽av亚洲精品天堂| avwww免费| 亚洲av男天堂| 男人舔女人的私密视频| 亚洲精品,欧美精品| 在线 av 中文字幕| 人体艺术视频欧美日本| 亚洲,一卡二卡三卡| 黄片小视频在线播放| 亚洲欧洲日产国产| 免费少妇av软件| 成人国语在线视频| 自线自在国产av| 如日韩欧美国产精品一区二区三区| 少妇 在线观看| 亚洲欧美一区二区三区黑人| 一本久久精品| 欧美日韩亚洲综合一区二区三区_| 亚洲色图综合在线观看| 午夜影院在线不卡| 一级毛片电影观看| 国产老妇伦熟女老妇高清| 国产高清不卡午夜福利| 久久 成人 亚洲| 久久精品国产亚洲av涩爱| 久久热在线av| 日韩电影二区| av又黄又爽大尺度在线免费看| 精品亚洲成国产av| 性少妇av在线| 99久久人妻综合| bbb黄色大片| 18禁观看日本| 操出白浆在线播放| 亚洲国产欧美日韩在线播放| 一区二区三区精品91| 久热爱精品视频在线9| 天天操日日干夜夜撸| 亚洲一区中文字幕在线| 亚洲欧美精品自产自拍| 一区二区三区乱码不卡18| 天天添夜夜摸| 丝袜美足系列| 看免费成人av毛片| 看免费av毛片| 十八禁网站网址无遮挡| 欧美日韩视频高清一区二区三区二| 1024香蕉在线观看| 亚洲精品国产色婷婷电影| 国产成人系列免费观看| 中国三级夫妇交换| 国产亚洲精品第一综合不卡| 色婷婷av一区二区三区视频| 亚洲少妇的诱惑av| 亚洲欧美成人综合另类久久久| 日韩制服丝袜自拍偷拍| 亚洲一卡2卡3卡4卡5卡精品中文| 建设人人有责人人尽责人人享有的| 一级,二级,三级黄色视频| 欧美日韩精品网址| 天天操日日干夜夜撸| 午夜激情av网站| 麻豆av在线久日| 亚洲精品成人av观看孕妇| 色精品久久人妻99蜜桃| 亚洲精品第二区| 只有这里有精品99| 999精品在线视频| 哪个播放器可以免费观看大片| 久久久久久人人人人人| 人人妻人人澡人人看| 涩涩av久久男人的天堂| 女人精品久久久久毛片| 亚洲精品一区蜜桃| 免费不卡黄色视频| 午夜福利视频在线观看免费| 最新在线观看一区二区三区 | av卡一久久| 久久精品人人爽人人爽视色| 亚洲国产欧美网| 欧美日韩一区二区视频在线观看视频在线| 亚洲欧美清纯卡通| 欧美日韩一区二区视频在线观看视频在线| 欧美人与性动交α欧美软件| 成人毛片60女人毛片免费| 国产精品99久久99久久久不卡 | 黄色视频在线播放观看不卡| 1024视频免费在线观看| 国产亚洲av高清不卡| 黄片无遮挡物在线观看| 999精品在线视频| 在线观看免费午夜福利视频| 女人被躁到高潮嗷嗷叫费观| 日韩精品有码人妻一区| 久久久久久人妻| 精品国产一区二区久久| 涩涩av久久男人的天堂| 777米奇影视久久| 人妻 亚洲 视频| 中国国产av一级| 熟女av电影| 久久综合国产亚洲精品| 看免费成人av毛片| 咕卡用的链子| 最近最新中文字幕大全免费视频 | 少妇 在线观看| 汤姆久久久久久久影院中文字幕| 丰满饥渴人妻一区二区三| 母亲3免费完整高清在线观看| 国产精品免费视频内射| 视频在线观看一区二区三区| 建设人人有责人人尽责人人享有的| 如日韩欧美国产精品一区二区三区| 女人被躁到高潮嗷嗷叫费观| 久久久精品国产亚洲av高清涩受| 精品人妻一区二区三区麻豆| 精品久久久久久电影网| 韩国精品一区二区三区| 精品一区二区免费观看| 天天躁狠狠躁夜夜躁狠狠躁| √禁漫天堂资源中文www| 两个人免费观看高清视频| 亚洲精品国产色婷婷电影| 欧美日本中文国产一区发布| 日韩精品免费视频一区二区三区| 精品一区在线观看国产| 黄片小视频在线播放| 亚洲,一卡二卡三卡| 国产av精品麻豆| 久久久久网色| 在线观看免费日韩欧美大片| 如何舔出高潮| 日韩一卡2卡3卡4卡2021年| 91成人精品电影| 欧美国产精品一级二级三级| 黄片小视频在线播放| 国产99久久九九免费精品| 亚洲av日韩精品久久久久久密 | 操美女的视频在线观看| 精品国产一区二区三区四区第35| 精品国产一区二区久久| 蜜桃国产av成人99| 又大又黄又爽视频免费| 国产亚洲av片在线观看秒播厂| 欧美日韩亚洲国产一区二区在线观看 | 三上悠亚av全集在线观看|