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

    有限供水條件下玉米灌水時間確定方法的研究

    2022-11-09 12:15:04時晴晴王仰仁李炎楊麗霞申孝軍李金玉劉宏武
    灌溉排水學(xué)報 2022年10期
    關(guān)鍵詞:降雨量灌水定額

    時晴晴,王仰仁,李炎,楊麗霞,申孝軍,李金玉,劉宏武

    有限供水條件下玉米灌水時間確定方法的研究

    時晴晴1,王仰仁1*,李炎1*,楊麗霞2,申孝軍1,李金玉3,劉宏武3

    (1.天津農(nóng)學(xué)院 水利工程學(xué)院,天津 300392;2.太原理工大學(xué) 水利科學(xué)與工程學(xué)院,太原 030024;3.山西省中心灌溉試驗站,太原 030012)

    【目的】充分發(fā)揮自然降水的生產(chǎn)力,最大限度提高有限灌溉水資源利用效率,研究灌溉預(yù)報中確定灌水時間的方法?!痉椒ā炕谧魑锷L模型和土壤水熱動態(tài)模擬模型,建立了玉米產(chǎn)量與水分關(guān)系模型,利用山西省文峪河試驗站2017—2020年玉米生長季試驗資料,通過最小二乘法確定了模型參數(shù)。并利用該站1957—2018年氣象資料以灌溉需水量從小到大的順序做頻率分析,確定了5個典型年(分別對應(yīng)5%、25%、50%、75%和95%這5個頻率);針對不同的灌水定額(45、75 mm和105 mm)和灌水次數(shù)(1、2、3、4次和5次)組合,優(yōu)化確定了相應(yīng)的灌水時間,以累計有效降水量為自變量,以考慮可供水量修正的累計蒸散量為因變量,分析構(gòu)建了相鄰2次灌水期間累計蒸散量和累計有效降水量的關(guān)系,基于該關(guān)系提出了確定灌水時間的K-M法;在此基礎(chǔ)上,利用2009—2018年的氣象資料,通過模擬分析比較了3種方法(K-M法、動態(tài)灌水下限法和經(jīng)驗灌水法)的增產(chǎn)效益。【結(jié)果】結(jié)果表明,不同灌水定額下,累計蒸散量和累計有效降水量關(guān)系的相關(guān)系數(shù)均能達(dá)到0.92以上,表明該模型在年際間和不同灌水定額之間均有較好的適用性;隨著灌水定額的增大,K-M法中參數(shù)呈增大的趨勢,表明參數(shù)一定程度上具有灌水定額的物理特性?!窘Y(jié)論】基于產(chǎn)量、效益和用水效率3方面對3種確定灌水時間方法的評價表明,與經(jīng)驗灌水法相比,以K-M法最好,動態(tài)灌水下限法次之。

    作物生長模型;灌水下限;蒸散量;灌水時間;限量灌溉

    0 引言

    【研究意義】灌水時間是影響作物產(chǎn)量的重要因素,灌水量一定時確定合理的灌水時間可以顯著提高作物產(chǎn)量,如冬小麥提前灌水,產(chǎn)量可增加475.3 kg/hm2,產(chǎn)量增幅10.1%;玉米推后灌水產(chǎn)量可增加734.6 kg/hm2,產(chǎn)量增幅9.5%[1]。但是由于年際間降雨量及其分布的不確定性,不同年得出的試驗結(jié)果是不相同的,為此,人們對作物合理灌水時間的確定進(jìn)行了大量的理論研究?!狙芯窟M(jìn)展】首先,利用作物水分生產(chǎn)函數(shù)(產(chǎn)量與供水量的關(guān)系),基于產(chǎn)量最大,對于給定的灌水量,通過優(yōu)化方法確定作物某一典型年的灌水時間,即通過優(yōu)化灌溉制度確定灌水時間[2-3]??紤]到年際間降雨量及其分布的隨機(jī)特性,崔遠(yuǎn)來等[4]基于水分敏感指數(shù)累積函數(shù)方法將作物生長期劃分為更短的時段,以降雨量為隨機(jī)變量,灌水定額和灌水時間為決策變量,以可供灌溉水量和田面蓄水深度為狀態(tài)變量,采用隨機(jī)動態(tài)規(guī)劃的方法確定了水稻的優(yōu)化灌溉制度。Raju等[5]將作物的干物質(zhì)量及土壤含水率作為狀態(tài)變量,灌水量作為決策變量,并用降水量作為隨機(jī)變量來修正動態(tài)規(guī)劃模型。采用隨機(jī)動態(tài)規(guī)劃方法確定的灌水時間,實(shí)際上給出了當(dāng)?shù)亟涤昙捌浞植紬l件下多年平均的優(yōu)化灌水時間,是現(xiàn)在廣泛采用的經(jīng)驗灌水方法中灌水時間確定的理論化方法。由于年際間降水量和蒸發(fā)量的隨機(jī)變化,該灌水時間不同于某一具體年份的最優(yōu)灌水時間。為了使得灌水時間盡可能地接近于具體年份的最優(yōu)灌水時間,對墑情預(yù)測和灌溉預(yù)報進(jìn)行了廣泛研究,提出了實(shí)時灌溉預(yù)報的方法。其核心是充分利用降雨、氣溫、光照等氣象、作物生長和土壤墑情等實(shí)時信息,對作物根系層土壤含水率做出精準(zhǔn)預(yù)測,當(dāng)土壤含水率小于給定的灌水下限值時做出灌水預(yù)報。實(shí)時灌溉預(yù)報有助于節(jié)水灌溉、增加作物產(chǎn)量、提高灌溉經(jīng)濟(jì)效益,但是,現(xiàn)狀的灌溉預(yù)報僅適用水資源能夠滿足作物需水要求的情況[6]。在限量供水條件下,孫景生等[7]通過確定小麥葉片光合速率的高值區(qū)域?qū)?yīng)的土壤水分,來確定適宜土壤水分上下限指標(biāo)??紤]到有限供水條件下,土壤含水率調(diào)控下限值應(yīng)該是變化的,王仰仁等[8]提出了動態(tài)灌水下限方法,按照該方法確定的灌水時間進(jìn)行灌水,較經(jīng)驗灌水有明顯的增產(chǎn)效果[9]?!厩腥朦c(diǎn)】但是基于土壤水分監(jiān)測的方法往往需要大量傳感器,且由于年際間降雨量強(qiáng)烈隨機(jī)變化特性,有限灌溉供水量條件下,干旱年份難以維持既定的灌水下限值。鑒于此,人們在分析利用多種方法求得實(shí)際蒸散量[10-13]的基礎(chǔ)上,通過蒸散量來確定灌水時間[14-21]。如:顧哲等[22]基于蒸散量和水量平衡,設(shè)計當(dāng)田間蒸發(fā)蒸騰總量大于土壤中可供作物利用水分時便啟動灌溉的灌水決策方法。杜江濤等[23]基于DSSAT模型,對新疆棉花膜下滴灌進(jìn)行了分析,提出了當(dāng)階段實(shí)際蒸散量與同期降雨量之差達(dá)到25 mm時灌溉的灌水決策方法。由于降雨量較大時會產(chǎn)生深層滲漏,甚至徑流等原因,簡單利用蒸散量與降雨量的日差值確定灌水時間會帶來較大誤差;另外,現(xiàn)在利用蒸散量確定灌水時間的方法仍然只適用于充分供水灌溉?!緮M解決的關(guān)鍵問題】針對以上問題,本研究依據(jù)作物蒸散量隨灌溉供水量變化規(guī)律,通過典型年灌溉制度優(yōu)化,給出不同灌水定額不同灌水次數(shù)條件下的優(yōu)化灌水時間;以優(yōu)化灌水時間為依據(jù),分析確定相鄰2次灌水之間(對于第1次灌水是指播種日到第1次灌水日之間)的累計蒸散量和累計有效降雨量的關(guān)系?;谠撽P(guān)系構(gòu)建提出確定灌水時間的方法,稱為蒸散量法。并利用2009—2018年的氣象資料分析比較了蒸散量法、動態(tài)灌水下限法和經(jīng)驗灌水方法的增產(chǎn)增收效益。

    1 材料與方法

    1.1 試驗區(qū)概況

    文峪河試驗區(qū)(E112°01′,N37°28′,海拔788.2 m),地處山西省中部呂梁山東側(cè)文峪河灌區(qū),研究區(qū)屬暖溫帶大陸性半干旱季風(fēng)氣候區(qū)。多年平均降雨量453.5 mm,平均無霜期160~180 d,年平均日照時間為3 004.0 h,多年平均氣溫10.1 ℃。中等肥力,土壤質(zhì)地為中壤土,0~100 cm 平均土壤體積質(zhì)量1.38 g/cm3,地下水埋深大于10 m。

    1.2 研究方法

    1.2.1 有限供水條件下灌水時間的優(yōu)化

    1)典型年的選取

    由于降雨量年際間和年內(nèi)分布不均,每年的灌溉需水量差異較大,所以在進(jìn)行灌溉制度設(shè)計的時候,需確定一個特定的水文年作為規(guī)劃設(shè)計的依據(jù),即:典型年。根據(jù)歷年灌溉需水量(作物潛在最大蒸發(fā)蒸騰量與降雨量的差值)從小到大進(jìn)行排序,做頻率分析,選擇與頻率5%、25%、50%、75%、95%相同或相近的年份作為典型年。

    2)灌水時間的優(yōu)化

    在選出典型年的基礎(chǔ)上,用規(guī)劃求解法進(jìn)行典型年的灌水時間(以播種日開始算起的天數(shù)表示)優(yōu)化計算,具體方法是針對某一個典型年,假定可供水量分別為1次水、2次水、…、等,利用Excel軟件進(jìn)行規(guī)劃求解,灌水按照全生育期均勻分布、播種后10~20 d、收獲前10~20 d、生育期中間前后7~8 d、生育期前末、生育期后期末和經(jīng)驗灌水7種不同的灌水時間初始值,例如灌1次水時,灌水時間按全生育期均勻分布初始值為18、70、128;灌2次水時,初始值為11/15、27/59、61/110、67/120、68/72、85/102、127/129(上述數(shù)據(jù)分別表示第幾日,下同);灌3次水時,初始值為11/15/18、27/45/57、58/71/90、61/85/110、65/85/105、70/74/78、123/125/127。對每種灌水次數(shù)進(jìn)行優(yōu)化計算,直到灌水次數(shù)對應(yīng)的凈效益開始減小時停止計算。灌溉制度優(yōu)化模型(式(1)),灌水定額設(shè)置45、75和105 mm,給定灌水次數(shù)=1、2、3、…,m,其中,m為最大灌水次數(shù),根據(jù)凈效益最大確定。

    目標(biāo)函數(shù)為:

    1.2.2 蒸散量法

    蒸散量法是基于蒸散量確定灌水時間的一種方法。本研究利用5個典型年灌溉制度優(yōu)化結(jié)果確定相鄰2次灌水期間的累計蒸散量和累計有效降雨量,考慮可供灌水次數(shù)對累計蒸散量的修正,建立了有效降雨量和蒸散量之間的定量關(guān)系,見式(2),并用修正的累計蒸散量與累計有效降雨量的相關(guān)系數(shù)[25]評價模型的擬合精度。

    1.2.3 其他確定灌水時間方法

    1)動態(tài)灌水下限法

    在有限供水條件下灌水下限值隨灌水時間和可供灌溉水量變化,式(3)給出了灌水下限值與可供水量和灌水時間的關(guān)系[26]。以式(3)計算的值為模擬值,以灌溉制度優(yōu)化確定的灌水前的含水率作為灌水下限實(shí)測值,以二者之間誤差平方和最小為目標(biāo)函數(shù)進(jìn)行規(guī)劃求解,可求得參數(shù)0、1、2、1、2。為了應(yīng)用方便,以相對含水率tr(占田間持水率的比例)表示含水率下限值,并用灌水下限實(shí)測值與式(3)計算的模擬值的相關(guān)系數(shù)[25]進(jìn)行擬合結(jié)果的評價。該方法以下簡稱為D-T法。

    式中:r為相對生長時間(d),r=/m,其中,為從播種日期算起的生長時間(d);m為最大生長天數(shù);0、1、2、1、2為待定參數(shù);為每次灌水前的可供灌溉水量(mm)。

    2)經(jīng)驗灌溉法

    經(jīng)驗灌溉法是指按照大量的試驗研究和生產(chǎn)實(shí)踐確定的灌水時間進(jìn)行灌溉的一種方法。本研究根據(jù)文獻(xiàn)[1]確定的灌水時間段的中間日作為經(jīng)驗灌水時間。以下將該方法簡稱為J-Y法。

    2 結(jié)果與分析

    2.1 產(chǎn)量與水分關(guān)系模型

    產(chǎn)量與水分關(guān)系模型是灌溉制度優(yōu)化的依據(jù),其模型參數(shù)包括土壤水分動態(tài)模擬模型參數(shù)和作物生長模型參數(shù),利用文峪河試驗站2017—2020年玉米生長季試驗資料,確定了相關(guān)參數(shù)(表1、表2)。其中利用分層土壤含水率進(jìn)行了土壤水分動態(tài)模型參數(shù)的率定,其0~60 cm土壤含水率平均值的模擬值與實(shí)測值的相關(guān)系數(shù)能達(dá)到0.876 4(樣本數(shù)為128);利用莖、葉、籽粒干物質(zhì)質(zhì)量的實(shí)測值對作物生長模型參數(shù)進(jìn)行了率定,其籽粒產(chǎn)量實(shí)測值與模擬值的相關(guān)系數(shù)能達(dá)到0.854 1(樣本數(shù)46個),說明2種模型的模擬精度均較高。由表2可見,g為0.472 0,與文獻(xiàn)[27]計算的轉(zhuǎn)化效率值(0.341 8~0.498 6)較為接近。

    表1 土壤水分動態(tài)模擬模型參數(shù)優(yōu)化結(jié)果

    注s為飽和含水率;r為殘余含水率;、為土壤水分特征曲線形狀參數(shù);s為土壤飽和導(dǎo)水率;為經(jīng)驗指數(shù);50為作物蒸騰量減小到最大可能蒸騰量50%時所對應(yīng)的土壤水勢[28]。

    表2 作物生長模型關(guān)鍵參數(shù)優(yōu)化結(jié)果

    注1234均為經(jīng)驗系數(shù);g為光合產(chǎn)物轉(zhuǎn)化效率[20]。

    2.2 典型年灌溉制度優(yōu)化

    根據(jù)玉米灌溉需水量進(jìn)行頻率分析確定了5個典型年(1971、1983、2017、1994年和2015年),不同典型年的水文信息見表3。從表3可以看出,隨著頻率的增大,玉米的灌溉需水量由-150.0 mm增加到301.3 mm;相應(yīng)地,生長期降雨量呈減少的變化趨勢,由5%典型年的571.3 mm減小到95%典型年的173.0 mm,變化幅度較大;參考作物蒸散量和潛在作物蒸散量隨頻率增大,沒有明顯的變化趨勢。

    表3 典型年水文信息

    當(dāng)灌水定額為45、75 mm和105 mm時,對于5%典型年灌1次水與不灌水比較,產(chǎn)量和效益均有所降低,即不需要灌溉(本研究所有典型年均未考慮澇漬災(zāi)害),因而表4只給出了其他4個典型年的優(yōu)化灌溉制度,表中灌水時間以播種后的時間表示。

    從表4可以看出,隨著灌水次數(shù)的增加,產(chǎn)量均增加,總體效益也在增加,但增加的幅度逐漸減小。如灌水定額為45 mm時,75%典型年,灌2次水較灌1次水產(chǎn)量增加了1.34 t/hm2,效益增加了1 783元/hm2;灌3次水較灌2次水產(chǎn)量增加了0.89 t/hm2,效益增加了978元/hm2;灌4次水較灌3次水,產(chǎn)量增加了0.53 t/hm2,效益增加了208元/hm2。同一典型年,隨著灌水定額的增大,最大灌水次數(shù)對應(yīng)的產(chǎn)量和效益都在減小。就75%典型年而言,灌水定額為45、75 mm和105 mm時,最大灌水次數(shù)分別為4、2次和2次,對應(yīng)的產(chǎn)量分別為11.53、10.69 t/hm2和10.57 t/hm2,對應(yīng)的效益分別為6 878、5 798元/hm2和4 863元/hm2。不同灌水次數(shù)條件下,不同典型年適宜的灌水定額不同。如灌1次水時,25%和50%典型年均以灌水定額為45 mm時,產(chǎn)量效益最高;75%和95%典型年均以灌水定額為105 mm時,產(chǎn)量效益最高。25%和50%典型年中隨著灌水定額增大,產(chǎn)量卻在降低。

    表4 不同典型年不同灌水定額的灌溉制度優(yōu)化結(jié)果

    2.3 K-M法和D-T法參數(shù)的確定結(jié)果及分析

    表5給出45、75、105 mm這3種灌水定額條件下2種決策方法的模型參數(shù)。由表5可知,對于K-M法,3種灌水定額的值都在0.9以上,該公式在不同灌水定額下擬合效果均較好,表明該方法在不同年份和不同灌水定額二方面均具有較好的適用性;且參數(shù)隨灌水定額的增加呈增大的趨勢,表明參數(shù)一定程度上具有灌水定額的物理特性。在D-T法中,當(dāng)灌水定額為75 mm時,值最大,為0.810 4,擬合效果最好,表明當(dāng)灌水定額為75 mm時,該法在不同年份上有較好的適用性。

    表5 不同灌水定額條件下K-M法和D-T法的參數(shù)

    2.4 增產(chǎn)效益對比分析

    利用2009—2018年的氣象資料,在給定3種灌水次數(shù)的條件下,按照3種確定灌水時間的方法,求出了不同灌水定額下相應(yīng)的灌水時間(表6)和產(chǎn)量效益(表7),并給出了2009—2018年的產(chǎn)量和效益的平均值。

    由于J-Y法確定的灌水時間每年都是相同的(灌1次水時間為播種后第61天,灌2次水為播種后第61、110天,灌3次水為播種后第61、85、110天),所以表6僅給出了K-M法和D-T法確定的灌水時間。由表6可以看出,不同灌水定額不同灌水次數(shù)條件下,K-M法確定的灌水時間均比D-T法確定的灌水時間晚。同一灌水次數(shù)情況下,灌水定額增大,出現(xiàn)不灌水的年份也增多,如灌3次水時,45 mm和75 mm的灌水定額,不灌水年份有2 a;105 mm的灌水定額,不灌水年份有3 a。主要原因是濕潤年份玉米生長期降雨量較多,灌水定額越大,深層滲漏越多,同樣施肥條件下作物吸收的養(yǎng)分越少,因而會造成不同程度的減產(chǎn),所以灌水定額越大,不灌水的年份越多。

    由表7可知,不同灌水定額和灌水次數(shù)條件下,K-M法的用水效率(灌溉水的增產(chǎn)量與灌溉水量的比值)均是3種方法中最高的。當(dāng)灌水定額為45 mm和105 mm時,無論灌幾次水,從產(chǎn)量和效益來看,K-M法均最高;當(dāng)灌水定額為75 mm時,從產(chǎn)量來看,不同灌水次數(shù)時,D-T法產(chǎn)量均高于K-M法和J-Y法;從效益來看,1次水和2次水時,D-T法最高;3次水時,K-M法最高。綜合分析,K-M法效果最好,其次是D-T法,J-Y法最差。

    表6 不同灌水定額在不同方法下的灌水時間

    表7 不同灌水定額在不同方法下的產(chǎn)量效益

    3 討論

    1)隨著灌水定額的增加,玉米千粒質(zhì)量、籽粒產(chǎn)量均呈先增大后減小的趨勢[29],與此不同,本研究發(fā)現(xiàn)隨著灌水定額的增大,產(chǎn)量有降低的趨勢(表4),主要原因是在產(chǎn)量與水分關(guān)系模型中,考慮了灌溉和降雨對土壤氮素的淋失使根系層氮素的濃度降低,由此導(dǎo)致作物產(chǎn)量的降低。

    2)蒸散量法與經(jīng)驗法相比較,優(yōu)勢在于考慮了年際間氣象要素變化對作物生長及灌水時間的影響,因而增產(chǎn)增收效果顯著;與灌水下限法比較,蒸散量法考慮了蒸散量隨供水量的減小而減小的特性,且克服了灌水下限的確定依賴于田間試驗的弊端,因而描述蒸散量法的模型有更高的擬合精度,更適合于限量供水條件下灌水時間的確定。經(jīng)驗法較差的主要原因是該法確定的灌水時間就多年平均而言,是最優(yōu)的,但對于某一個具體年份而言,一般不是最優(yōu)的。

    3)以蒸散量為基礎(chǔ)確定灌水時間的研究中,有些只考慮蒸散量對灌溉的影響,忽略了土壤和作物本身因素[30],有些只適合于溫室及膜下滴灌[22-23],有些以蒸發(fā)蒸騰量和深層滲漏量來預(yù)測灌溉時間[31]。以上研究均沒有考慮有限供水對灌溉時間的影響。為了能更好地適用于限量灌溉,本研究引入可供灌溉水次數(shù),更好地提高蒸發(fā)蒸騰量與有效降雨量之間的關(guān)系精度,從而更精準(zhǔn)地確定灌水時間。

    4)本研究在確定K-M法和D-T法的參數(shù)時發(fā)現(xiàn),不同灌水定額條件下各個參數(shù)均有差異,所以,對于不同的灌水定額,應(yīng)當(dāng)運(yùn)用與之相對應(yīng)的模型參數(shù),以便由此確定的灌水時間更準(zhǔn)確。且在確定K-M法參數(shù)時,值均能達(dá)到0.92以上,表明模型的擬合精度較高,結(jié)果可信度也較高,這為大田確定灌水時間的實(shí)際應(yīng)用提供了依據(jù)?;诖_定K-M法參數(shù)的方法,提供了一種確定灌水時間的理論化方法,避免了傳統(tǒng)灌水下限法單純依賴田間試驗的弊端,豐富了非充分灌溉理論。用K-M法代替以土壤含水率為指標(biāo)的灌水下限法,不僅適用于有限供水灌溉,也適用于充分供水灌溉,同時也可以將K-M法編程,植入灌溉控制器[32],實(shí)施自動化灌水。

    4 結(jié)論

    1)基于灌溉制度優(yōu)化結(jié)果確定了研究區(qū)不同灌水定額條件下,玉米灌水控制下限與可供灌溉水量和灌水時間的關(guān)系(D-T法);提出了基于蒸散量和有效降雨量確定灌水時間的方法(K-M法)。為有限供水條件下灌區(qū)農(nóng)田水分精準(zhǔn)調(diào)控提供依據(jù)。

    2)基于產(chǎn)量、效益和用水效率等對3種確定灌水時間方法綜合評價結(jié)果表明,K-M法確定灌水時間更為精確,動態(tài)灌水下限法次之,經(jīng)驗法最差。

    [1] 王仰仁, 孫小平. 山西農(nóng)業(yè)節(jié)水理論與作物高效用水模式[M]. 北京: 中國科學(xué)技術(shù)出版社, 2003.

    WANG Yangren, SUN Xiaoping. Shanxi agricultural water-saving theory and crop high-efficiency water use model[M]. Beijing: China Science and Technology Press, 2003.

    [2] FLINN J C, MUSGRAVE W F Development and analysis of input-output relations for irrigation water[J]. Australian Journal of Agricultural Economics. 1967, 11(1): 1-19.

    [3] 榮豐濤.節(jié)水型農(nóng)田灌溉制度的初步研究[J].水利水電技術(shù),1986, 17(7): 17-21.

    RONG Fengtao.Preliminary study on irrigation degree of water-saving farmland[J].Water Resources and Hydropower Engineering, 1986, 17(7): 17-21.

    [4] 崔遠(yuǎn)來, 袁宏源, 李遠(yuǎn)華. 考慮隨機(jī)降雨時稻田高效節(jié)水灌溉制度[J]. 水利學(xué)報,1999, 30(7): 40-45.

    CUI Yuanlai, YUAN Hongyuan, LI Yuanhua. The optimum irrigation scheduling for rice with stochastic rainfall under limited water supply[J]. Journal of Hydraulic Engineering, 1999, 30(7): 40-45.

    [5] RAJU K S, BIERE A W, ANEMASU E T, et al. Irrigation scheduling based on a dynamic crop response model[J]. In D.Hi11e1, ed. Advances in irrigation, Orlando, Florida, USA, Academic press, 1983, 54(2): 134-141.

    [6] 茆智, 李遠(yuǎn)華, 李會昌. 實(shí)時灌溉預(yù)報[J]. 中國工程科學(xué),2002, 4(5): 24-33.

    MAO Zhi, LI Yuanhua, LI Huichang. Real-time forecast of irrigation[J]. Engineering Science, 2002, 4(5): 24-33.

    [7] 孫景生, 劉祖貴, 肖俊夫, 等. 冬小麥節(jié)水灌溉的適宜土壤水分上、下限指標(biāo)研究[J]. 中國農(nóng)村水利水電, 1998(9): 10-12.

    SUN Jingsheng, LIU Zugui, XIAO Junfu, et al. On the research of upper and lower limit indicators of suitable soil moisture for water-saving irrigation of winter wheat[J]. China Rural Nater and Hydropower, 1998 (9): 10-12.

    [8] 王仰仁, 杜秀文, 張紹強(qiáng). 限量供水條件下灌溉預(yù)報增產(chǎn)效益分析[J]. 中國農(nóng)村水利水電, 2016(12): 1-7.

    WANG Yangren, DU Xiuwen, ZHANG Shaoqiang. Analytical investigation on the irrigation forecast and its benefit of increasing yield under limited supply water[J]. China Rural Water and Hydropower, 2016(12): 1-7.

    [9] 杜娟娟, 王仰仁, 李粉嬋. 基于動態(tài)灌水下限值的冬小麥非充分灌溉預(yù)報研究[J]. 節(jié)水灌溉, 2017(5): 93-97.

    DU Juanjuan, WANG Yangren, LI Fenchan. Research on insufficient irrigation forecasting for winter wheat based on lower limit value of dynamic irrigation amount[J]. Water Saving Irrigation, 2017(5): 93-97.

    [10] MCCAUELY D M, NACKLEY L L. Development of mini-lysimeter system for use in irrigation automation of container-grown crops[J] HardwareX, 2022(11): 1-10.

    [11] SINGH R P, PARAMANIK S, BHATTACHARYA B K, et al. Modelling of evapotranspiration using land surface energy balance and thermal infrared remote sensing[J]. Tropical Ecology, 2020,61(6): 42-50.

    [12] KHAN M S, JEONG J, CHOI M. An improved remote sensing based approach for predicting actual Evapotranspiration by integrating LiDAR[J]. Advances in Space Research, 2021, 68(4): 1 732-1 753.

    [13] REZAEI M, GHASEMIEH H, ABDOLLAHI K. Simplified version of the METRIC model for estimation of actual evapotranspiration[J]. International Journal of Remote Sensing, 2021, 42(14): 5 568-5 599.

    [14] 趙穎, 紀(jì)建偉, 崔會坤, 等. 基于作物蒸散量模型的新型滑蓋溫室智能灌溉系統(tǒng)設(shè)計[J]. 節(jié)水灌溉, 2017(8): 83-87, 91.

    ZHAO Ying, JI Jianwei, CUI Huikun, et al. Design of A new type intelligent irrigation system based on crop evapotranspiration model[J]. Water Saving Irrigation, 2017(8): 83-87, 91.

    [15] 劉小飛, 張寄陽, 孫景生, 等.自動補(bǔ)水蒸發(fā)皿裝置的原理及應(yīng)用[J]. 排灌機(jī)械工程學(xué)報, 2012, 30(1): 80-84.

    LIU Xiaofei, ZHANG Jiyang, SUN Jingsheng, et al. Principle and application of evaporation pan with device of water automatic replenishing[J]. Journal of Drainage and Irrigation Machinery Engineering, 2012, 30(1): 80-84.

    [16] 任夢之. 基于蒸發(fā)互補(bǔ)原理的新疆干旱灌區(qū)節(jié)水灌溉效果評估[D]. 西安: 西安理工大學(xué), 2021.

    REN Mengzhi. Evaluation of water-saving irrigation effect in Xinjiang arid irrigation area based on evaporation complementary principle[D]. Xi’an: Xi’an University of Technology, 2021.

    [17] 張伶鳦, 王潤濤, 張長利, 等. 基于調(diào)虧理論和模糊控制的寒地水稻智能灌溉策略[J]. 農(nóng)業(yè)工程學(xué)報, 2016, 32(13): 52-58.

    ZHANG Lingyi, WANG Runtao, ZHANG Changli, et al. Intelligent irrigation strategy based on regulated deficit theory and fuzzy control for rice in cold region[J]. Transactions of the Chinese Society of Agricultural Engineering, 2016, 32(13): 52-58.

    [18] PEREA R G, GARCíA I F, ARROYO M M, et al. Multiplatform application for precision irrigation scheduling in strawberries[J]. Agricultural Water Management, 2017, 183: 194-201.

    [19] 田強(qiáng)明, 溫宗周, 李麗敏, 等. 基于APSO-ELM和模糊邏輯的灌溉時間決策[J]. 中國農(nóng)村水利水電, 2020(4): 124-128.

    TIAN Qiangming, WEN Zongzhou, LI Limin, et al. Irrigation time decision based on APSO-ELM and fuzzy logic system[J]. China Rural Water and Hydropower, 2020(4): 124-128.

    [20] THYSEN I, DETLEFSEN N K. Online decision support for irrigation for farmers[J]. Agricultural Water Management, 2006, 86(3): 269-276.

    [21] 張明岳, 李麗敏, 溫宗周, 等.基于改進(jìn)Elman神經(jīng)網(wǎng)絡(luò)和模糊控制的智能灌溉算法設(shè)計[J]. 國外電子測量技術(shù), 2021, 40(11): 155-160.

    ZHANG Mingyue, LI Limin, WEN Zongzhou, et al. Design of intelligent irrigation algorithm based on improved Elman neural network and fuzzy control[J]. Foreign Electronic Measurement Technology, 2021, 40(11): 155-160.

    [22] 顧哲, 袁壽其, 齊志明, 等. 基于ET和水量平衡的日光溫室實(shí)時精準(zhǔn)灌溉決策及控制系統(tǒng)[J]. 農(nóng)業(yè)工程學(xué)報, 2018, 34(23): 101-108.

    GU Zhe, YUAN Shouqi, QI Zhiming, et al. Real-time precise irrigation scheduling and control system in solar greenhouse based on ET and water balance[J]. Transactions of the Chinese Society of Agricultural Engineering, 2018, 34(23): 101-108.

    [23] 杜江濤, 張楠, 龔珂寧, 等. 基于DSSAT模型的南疆膜下滴灌棉花灌溉制度優(yōu)化[J]. 生態(tài)學(xué)雜志, 2021, 40(11): 3 760-3 768.

    DU Jiangtao, ZHANG Nan, GONG Kening, et al. Optimization of cotton irrigation schedule under mulch drip irrigation in southern Xinjiang based on DSSAT model[J]. Chinese Journal of Ecology, 2021, 40(11): 3 760-3 768.

    [24] 王鐵英, 王仰仁, 柴俊芳, 等. 根系水分脅迫響應(yīng)函數(shù)對土壤水及作物生長動態(tài)和產(chǎn)量模擬影響的研究[J]. 干旱區(qū)地理, 2022, 45(2): 566-577.

    WANG Tieying, WANG Yangren, CHAI Junfang, et al. Effect of root water stress response function on soil water, crop growth dynamics and yield simulation[J]. Arid Land Geography, 2022, 45(2): 566-577.

    [25] 明道緒. 田間試驗與統(tǒng)計分析[M]. 3版. 北京: 科學(xué)出版社, 2013.

    MING Daoxu. Field experiments and statistical analysis[M]. Beijing: Science Press, 2013.

    [26] 姚麗. 限量供水條件下精準(zhǔn)灌溉技術(shù)集成效益分析[D]. 天津: 天津農(nóng)學(xué)院, 2020.

    YAO Li. Benefit analysis of precision irrigation technology integration under limited water supply[D].Tianjin: Tianjin Agricultural University, 2020.

    [27] 王佳旭, 王宏偉, 姜文野, 等. 不同種植方式對玉米干物質(zhì)積累、分配和產(chǎn)量的影響[J]. 玉米科學(xué), 2021, 29(5): 128-136.

    WANG Jiaxu, WANG Hongwei, JIANG Wenye, et al. Effects of different planting patterns on dry matter accumulation, distribution and yield of maize[J]. Journal of Maize Sciences, 2021, 29(5): 128-136.

    [28] 尚松浩, 毛曉敏, 雷志棟. 土壤水分動態(tài)模擬模型及其應(yīng)用[M]. 北京: 科學(xué)出版社, 2009.

    SHANG Songhao, MAO Xiaomin, LEI Zhidong. Dynamic simulation model of soil moisture and its application[M]. Beijing: Science Press, 2009.

    [29] 曹雪松, 鄭和祥, 苗平, 等. 引黃滴灌條件下水氮互作對玉米耗水特性及產(chǎn)量的影響[J]. 灌溉排水學(xué)報, 2022, 41(3): 33-39.

    CAO Xuesong, ZHENG Hexiang, MIAO Ping, et al.The effect of drip fertigation with yellow river water on water consumption and yield of summer maize[J].Journal of Irrigation and Drainage, 2022, 41(3): 33-39.

    [30] BARTLETT A C, ANDALES A A, ARABI M, et al. A smartphone app to extend use of a cloud - based irrigation scheduling tool[J].Computers and Electronics in Agriculture,2015, 111: 127-130.

    [31] 武開福. 基于調(diào)虧理論的小麥智能灌溉設(shè)計[J]. 水利技術(shù)監(jiān)督, 2017, 25(3): 34-37.

    WU Kaifu. Intelligent irrigation design of wheat based on deficit adjustment theory[J]. Technical Supervision in Water Resources. 2017, 25(3): 34-37.

    [32] 曲頌, 劉玉敏, 宋博, 等. 自動灌溉控制器設(shè)計與實(shí)現(xiàn)[J]. 黑龍江大學(xué)工程學(xué)報, 2021, 12(4): 72-77.

    QU Song, LIU Yumin, SONG Bo, et al. Design and implementation of automatic irrigation controller[J]. Journal of Engineering of Heilongjiang University, 2021, 12(4): 72-77.

    Methods to Determine the Timing of Maize Irrigation When Irrigation Amount is Restricted

    SHI Qingqing1, WANG Yangren1*, LI Yan1*, YANG Lixia2, SHEN Xiaojun1, LI Jinyu3, LIU Hongwu3

    (1. College of Hydraulic Engineering, Tianjin Agricultural University, Tianjin 300392, China;2. School of Hydraulic Science and Engineering, Taiyuan University of Technology, Taiyuan 030024, China;3. Shanxi Provincial Central Irrigation Experimental Station, Taiyuan 030012, China)

    【Objective】To reverse the continued decline of dwindling water resources in northern China, most regions have proposed to restrict groundwater extraction for irrigation. Optimizing irrigation timing and making most of natural precipitation and irrigation water is essential to minimizing crop yield loss due to the reduced irrigation amount. The purpose of this paper is to study the optimization of timing of maize irrigation. 【Method】Models relating maize yield and water usage were established based on crop growth model and soil water-heat coupling model. The model parameters were determined based on data measured from 2017—2020 from a maize field at Wenyuhe Experiment Station in Shanxi province. Meteorological data measured from 1957—2018 at the station was used in the frequency analysis for determining the irrigation amount by dividing the frequency from 5% to 95%. We optimized the timing of the irrigation assuming that the irrigation amount varied from 45 to 105 mm, and the times of irrigation varied from 1 to 5. The K-M model was used to determine the relationship between cumulative evapotranspiration and cumulative precipitation between two adjacent irrigations. Using the meteorological data measured from 2009 to 2018, the model was applied to optimize the irrigation timing, from which we compared the yield-increasing benefits obtained from the K-M method, the dynamic irrigation water limit method, and the empirical irrigation method. 【Result】The correlation coefficient between the cumulative evapotranspiration and the cumulative effective precipitation under different irrigation amount was greater than 0.92. With the increase in irrigation amount, the parameter M in the K-M method, which characterizes the irrigation amount, increased. 【Conclusion】Evaluation of the three methods for determining the optimal irrigation timing based on yield, benefit and water use efficiency showed that, compared with the empirical irrigation method, the K-M method was the best, followed by the dynamic irrigation water limit method.

    crop growth model; low limit of irrigation; evapotranspiration method; irrigation time; limited irrigation

    時晴晴, 王仰仁, 李炎, 等. 有限供水條件下玉米灌水時間確定方法的研究[J]. 灌溉排水學(xué)報, 2022, 41(10): 51-57.

    SHI Qingqing, WANG Yangren, LI Yan, et al. Methods to Determine the Timing of Maize Irrigation When Irrigation Amount is Restricted[J]. Journal of Irrigation and Drainage, 2022, 41(10): 51-57.

    1672 - 3317(2022)10 - 0051 - 07

    S275

    A

    10.13522/j.cnki.ggps.2021598

    2021-12-03

    國家自然科學(xué)基金項目(51779174);天津市科技計劃項目(20YDTPJC01450);天津市“項目+團(tuán)隊”重點(diǎn)培養(yǎng)專項(XB202016)

    時晴晴(1996-),女。碩士研究生,主要研究方向為灌溉排水理論與新技術(shù)。E-mail: 1240777143@qq.com

    王仰仁(1962-),男。教授,博士,主要研究方向為灌溉排水理論與新技術(shù)。E-mail: wyrf@163.com

    李炎(1973-),男。副教授,博士,主要研究方向為灌溉排水理論與新技術(shù)。E-mail: bfsqz0922@163.com

    責(zé)任編輯:趙宇龍

    猜你喜歡
    降雨量灌水定額
    公立醫(yī)院零修工程預(yù)算與定額庫編制問題探討
    建筑工程預(yù)結(jié)算中的定額與清單計價應(yīng)用
    降雨量與面積的關(guān)系
    灌水取球
    番茄灌水掌握技巧
    冬季棚菜灌水四關(guān)鍵
    “三新三化”在LNG接收站概算定額標(biāo)準(zhǔn)中應(yīng)用的探討
    損害賠償舉證責(zé)任的雙刃作用——對提高定額賠償立法趨勢的質(zhì)疑
    洞庭湖區(qū)降雨特性分析
    灌水秘笈
    美女脱内裤让男人舔精品视频| 一本大道久久a久久精品| 亚洲精品成人av观看孕妇| 精品国产一区二区三区四区第35| 成在线人永久免费视频| 十八禁人妻一区二区| 在线av久久热| 久久精品亚洲熟妇少妇任你| 男女午夜视频在线观看| 男女之事视频高清在线观看 | 国产精品一区二区精品视频观看| 日韩视频在线欧美| 99精国产麻豆久久婷婷| 一级毛片 在线播放| 高清不卡的av网站| 亚洲国产欧美日韩在线播放| 亚洲中文av在线| 午夜av观看不卡| 亚洲欧洲日产国产| 亚洲精品一卡2卡三卡4卡5卡 | 一本一本久久a久久精品综合妖精| 国产精品久久久久久人妻精品电影 | 老汉色av国产亚洲站长工具| 操出白浆在线播放| 丁香六月欧美| 午夜视频精品福利| 99国产精品免费福利视频| 老司机靠b影院| 你懂的网址亚洲精品在线观看| 视频在线观看一区二区三区| 精品一区在线观看国产| 亚洲欧美日韩另类电影网站| 亚洲男人天堂网一区| 婷婷色麻豆天堂久久| 国产亚洲午夜精品一区二区久久| 亚洲国产精品成人久久小说| 国产欧美日韩一区二区三 | 只有这里有精品99| 校园人妻丝袜中文字幕| www.自偷自拍.com| 中文字幕人妻丝袜制服| 18禁裸乳无遮挡动漫免费视频| 午夜91福利影院| 99久久人妻综合| 女性生殖器流出的白浆| 国产精品秋霞免费鲁丝片| 国产av国产精品国产| 黄色 视频免费看| 中文字幕人妻丝袜一区二区| 女警被强在线播放| 99九九在线精品视频| 久久久久久久大尺度免费视频| 国产精品二区激情视频| 欧美日韩视频精品一区| 黑人欧美特级aaaaaa片| 欧美 日韩 精品 国产| 国产欧美日韩精品亚洲av| 欧美成狂野欧美在线观看| 国产一区二区 视频在线| 亚洲中文字幕日韩| 亚洲中文av在线| 日本一区二区免费在线视频| 一本综合久久免费| 日韩视频在线欧美| 久久99精品国语久久久| 久久国产亚洲av麻豆专区| 悠悠久久av| 九草在线视频观看| 最新在线观看一区二区三区 | 熟女av电影| 人妻一区二区av| cao死你这个sao货| 色综合欧美亚洲国产小说| 国产精品久久久久成人av| 日本91视频免费播放| 久久中文字幕一级| 一边亲一边摸免费视频| 亚洲熟女毛片儿| 亚洲男人天堂网一区| 日本av免费视频播放| 校园人妻丝袜中文字幕| 91字幕亚洲| 免费高清在线观看日韩| 国产亚洲欧美精品永久| 亚洲五月色婷婷综合| 国产又色又爽无遮挡免| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日韩国产mv在线观看视频| 男女午夜视频在线观看| 国产免费福利视频在线观看| svipshipincom国产片| 亚洲欧美日韩另类电影网站| 亚洲,欧美,日韩| 日本色播在线视频| av电影中文网址| 精品少妇内射三级| 亚洲精品一区蜜桃| 母亲3免费完整高清在线观看| 午夜激情久久久久久久| 亚洲专区国产一区二区| 黄色视频不卡| 久久久精品国产亚洲av高清涩受| 在现免费观看毛片| 免费看不卡的av| 久久九九热精品免费| 50天的宝宝边吃奶边哭怎么回事| 老熟女久久久| 国产有黄有色有爽视频| 亚洲精品国产色婷婷电影| 日韩中文字幕欧美一区二区 | 中文字幕精品免费在线观看视频| 热99久久久久精品小说推荐| 女警被强在线播放| 一级毛片我不卡| 欧美在线黄色| 超色免费av| 久久青草综合色| 丝瓜视频免费看黄片| 国产有黄有色有爽视频| 青春草亚洲视频在线观看| 超碰97精品在线观看| 欧美激情极品国产一区二区三区| 亚洲国产欧美在线一区| 久久久久久免费高清国产稀缺| 五月开心婷婷网| 七月丁香在线播放| 久久国产精品影院| 亚洲视频免费观看视频| www.av在线官网国产| 99久久综合免费| 午夜91福利影院| 欧美日韩精品网址| 人人妻,人人澡人人爽秒播 | 80岁老熟妇乱子伦牲交| 在线天堂中文资源库| 国产有黄有色有爽视频| 国产色视频综合| 国产男女内射视频| 99久久99久久久精品蜜桃| 建设人人有责人人尽责人人享有的| 国产xxxxx性猛交| 国产精品一区二区在线观看99| 久久性视频一级片| 国产视频首页在线观看| 老司机影院成人| 亚洲国产精品国产精品| 欧美日本中文国产一区发布| 色精品久久人妻99蜜桃| 2018国产大陆天天弄谢| 在线观看免费高清a一片| 亚洲少妇的诱惑av| 国产成人影院久久av| 亚洲色图综合在线观看| 丝袜喷水一区| 亚洲色图 男人天堂 中文字幕| 91精品伊人久久大香线蕉| 国产精品麻豆人妻色哟哟久久| 日韩人妻精品一区2区三区| 国产成人av激情在线播放| 伊人久久大香线蕉亚洲五| 性色av乱码一区二区三区2| 伦理电影免费视频| 国产野战对白在线观看| 日本av免费视频播放| 久久精品国产综合久久久| 男女无遮挡免费网站观看| 午夜av观看不卡| 精品久久久精品久久久| 国产成人一区二区在线| 中文精品一卡2卡3卡4更新| 99热全是精品| 狠狠婷婷综合久久久久久88av| 欧美另类一区| 久久人人97超碰香蕉20202| 中文字幕人妻丝袜制服| 啦啦啦在线免费观看视频4| 少妇精品久久久久久久| 91国产中文字幕| 一级毛片我不卡| 考比视频在线观看| 两人在一起打扑克的视频| 女性被躁到高潮视频| 国产成人一区二区在线| 国产亚洲欧美精品永久| av天堂在线播放| 日韩人妻精品一区2区三区| 啦啦啦在线免费观看视频4| 亚洲第一av免费看| 免费在线观看黄色视频的| 成年人午夜在线观看视频| 男女床上黄色一级片免费看| 国产亚洲欧美精品永久| 十八禁高潮呻吟视频| 大话2 男鬼变身卡| 丁香六月欧美| av欧美777| 精品国产乱码久久久久久小说| www日本在线高清视频| 欧美黄色片欧美黄色片| 色婷婷av一区二区三区视频| 飞空精品影院首页| 亚洲精品自拍成人| 老司机深夜福利视频在线观看 | 嫁个100分男人电影在线观看 | 丝袜人妻中文字幕| 热99国产精品久久久久久7| 中文字幕人妻丝袜一区二区| 欧美国产精品一级二级三级| 午夜福利视频在线观看免费| 国产欧美日韩综合在线一区二区| 99热全是精品| 纵有疾风起免费观看全集完整版| 精品卡一卡二卡四卡免费| 午夜福利乱码中文字幕| 女人精品久久久久毛片| 天天躁夜夜躁狠狠躁躁| 欧美精品一区二区免费开放| 国产黄色视频一区二区在线观看| 成年美女黄网站色视频大全免费| 黄片播放在线免费| 美女午夜性视频免费| 欧美人与性动交α欧美精品济南到| 亚洲自偷自拍图片 自拍| 中文字幕人妻丝袜一区二区| 激情五月婷婷亚洲| 欧美97在线视频| 观看av在线不卡| 国产片特级美女逼逼视频| 国产一区二区三区综合在线观看| 黑人猛操日本美女一级片| 精品亚洲乱码少妇综合久久| 只有这里有精品99| 国产精品99久久99久久久不卡| 久久精品人人爽人人爽视色| 欧美久久黑人一区二区| avwww免费| 国产精品秋霞免费鲁丝片| 欧美成人午夜精品| av天堂在线播放| 母亲3免费完整高清在线观看| 少妇猛男粗大的猛烈进出视频| 亚洲人成网站在线观看播放| 大码成人一级视频| 丁香六月天网| 久久久精品国产亚洲av高清涩受| 亚洲国产欧美日韩在线播放| 久久性视频一级片| 亚洲精品日本国产第一区| 国产精品九九99| 男女边吃奶边做爰视频| 天天影视国产精品| 99久久精品国产亚洲精品| 性色av一级| 国产免费福利视频在线观看| 亚洲国产欧美网| 男女床上黄色一级片免费看| 中文字幕最新亚洲高清| 精品视频人人做人人爽| 精品欧美一区二区三区在线| 国产野战对白在线观看| 精品一区二区三区四区五区乱码 | 国产精品久久久久久人妻精品电影 | 亚洲美女黄色视频免费看| 热99国产精品久久久久久7| 成人免费观看视频高清| 脱女人内裤的视频| av视频免费观看在线观看| 桃花免费在线播放| 欧美+亚洲+日韩+国产| 18在线观看网站| 国产老妇伦熟女老妇高清| 免费高清在线观看日韩| 亚洲中文日韩欧美视频| 成年人免费黄色播放视频| 国产精品三级大全| 精品一品国产午夜福利视频| 777米奇影视久久| 亚洲图色成人| 亚洲av成人精品一二三区| 国产亚洲一区二区精品| 欧美亚洲 丝袜 人妻 在线| 久久久久网色| 国产福利在线免费观看视频| 高清欧美精品videossex| 五月开心婷婷网| 黑人欧美特级aaaaaa片| 国产伦理片在线播放av一区| av天堂久久9| 国产免费又黄又爽又色| 成人国语在线视频| 黑人巨大精品欧美一区二区蜜桃| 日本一区二区免费在线视频| 午夜精品国产一区二区电影| 成人黄色视频免费在线看| 久久人人爽人人片av| 久久99一区二区三区| 亚洲国产最新在线播放| 夫妻性生交免费视频一级片| 中文字幕高清在线视频| 九草在线视频观看| 久久 成人 亚洲| 深夜精品福利| 久久热在线av| 尾随美女入室| 久久热在线av| 亚洲,欧美,日韩| 中文字幕最新亚洲高清| 1024香蕉在线观看| 丰满人妻熟妇乱又伦精品不卡| 午夜影院在线不卡| xxx大片免费视频| 女人高潮潮喷娇喘18禁视频| 操美女的视频在线观看| 女人高潮潮喷娇喘18禁视频| 亚洲av男天堂| 久久久国产欧美日韩av| 久久久久视频综合| 在线观看一区二区三区激情| 狂野欧美激情性xxxx| 精品人妻1区二区| 免费在线观看黄色视频的| 国产无遮挡羞羞视频在线观看| 1024视频免费在线观看| 侵犯人妻中文字幕一二三四区| 肉色欧美久久久久久久蜜桃| 婷婷色综合www| 成人亚洲欧美一区二区av| av国产精品久久久久影院| 欧美人与善性xxx| 国产高清国产精品国产三级| a级片在线免费高清观看视频| 国产精品欧美亚洲77777| 巨乳人妻的诱惑在线观看| 久久久久久亚洲精品国产蜜桃av| 丁香六月天网| 精品久久蜜臀av无| 久久国产精品男人的天堂亚洲| av片东京热男人的天堂| 天天躁狠狠躁夜夜躁狠狠躁| 91麻豆av在线| 精品视频人人做人人爽| 欧美国产精品一级二级三级| 国产不卡av网站在线观看| www.熟女人妻精品国产| 久久天堂一区二区三区四区| 两个人看的免费小视频| 欧美黑人欧美精品刺激| av国产精品久久久久影院| 精品人妻熟女毛片av久久网站| 黄片播放在线免费| 国产伦理片在线播放av一区| 色视频在线一区二区三区| 一本一本久久a久久精品综合妖精| 制服人妻中文乱码| 亚洲第一青青草原| 久久影院123| 欧美变态另类bdsm刘玥| 久久国产精品人妻蜜桃| 亚洲久久久国产精品| 午夜免费鲁丝| 在线精品无人区一区二区三| 啦啦啦视频在线资源免费观看| 亚洲国产日韩一区二区| 亚洲成人手机| 精品国产乱码久久久久久男人| 女人爽到高潮嗷嗷叫在线视频| 免费高清在线观看日韩| 婷婷色综合www| 免费女性裸体啪啪无遮挡网站| 亚洲欧美清纯卡通| 国产女主播在线喷水免费视频网站| 在线 av 中文字幕| cao死你这个sao货| 欧美日韩av久久| 久久免费观看电影| 久久久精品区二区三区| 最黄视频免费看| 久久久久久人人人人人| 亚洲伊人久久精品综合| 91成人精品电影| 久久久久视频综合| 欧美黑人精品巨大| 成人18禁高潮啪啪吃奶动态图| 欧美黄色片欧美黄色片| 亚洲精品av麻豆狂野| 大香蕉久久网| 国语对白做爰xxxⅹ性视频网站| 欧美日韩亚洲综合一区二区三区_| 咕卡用的链子| 男女之事视频高清在线观看 | 欧美日韩综合久久久久久| 天天躁日日躁夜夜躁夜夜| 中文字幕人妻丝袜制服| 欧美日韩精品网址| 免费人妻精品一区二区三区视频| 国产一级毛片在线| 日韩制服骚丝袜av| 免费少妇av软件| 午夜激情久久久久久久| 亚洲视频免费观看视频| 看免费av毛片| 亚洲人成77777在线视频| 一本综合久久免费| 亚洲av欧美aⅴ国产| 国产精品99久久99久久久不卡| 国产黄色免费在线视频| 亚洲国产欧美日韩在线播放| 日韩制服骚丝袜av| 国产xxxxx性猛交| √禁漫天堂资源中文www| 另类亚洲欧美激情| 黑人欧美特级aaaaaa片| 国产精品.久久久| 99久久99久久久精品蜜桃| 精品国产乱码久久久久久男人| 国产精品熟女久久久久浪| 欧美精品人与动牲交sv欧美| 一区二区三区激情视频| 久久精品久久久久久噜噜老黄| av天堂在线播放| 日韩中文字幕视频在线看片| 少妇被粗大的猛进出69影院| 精品少妇久久久久久888优播| 性高湖久久久久久久久免费观看| 日韩免费高清中文字幕av| 国产福利在线免费观看视频| 日韩中文字幕欧美一区二区 | 国产日韩欧美在线精品| 多毛熟女@视频| 国产伦理片在线播放av一区| 美女福利国产在线| 亚洲国产精品一区二区三区在线| a 毛片基地| 超色免费av| 中文精品一卡2卡3卡4更新| 久热这里只有精品99| 丰满迷人的少妇在线观看| 国产欧美亚洲国产| 亚洲国产最新在线播放| 一级片免费观看大全| 久久人人爽人人片av| 免费观看a级毛片全部| 久久精品国产a三级三级三级| 欧美黄色淫秽网站| 成人午夜精彩视频在线观看| 青青草视频在线视频观看| 国产一区二区激情短视频 | 午夜影院在线不卡| 中文字幕最新亚洲高清| 日韩熟女老妇一区二区性免费视频| 日韩免费高清中文字幕av| av网站免费在线观看视频| 亚洲三区欧美一区| 日本wwww免费看| 日本91视频免费播放| 亚洲国产中文字幕在线视频| 免费观看人在逋| 欧美精品啪啪一区二区三区 | 别揉我奶头~嗯~啊~动态视频 | 狠狠婷婷综合久久久久久88av| 欧美日韩av久久| 青草久久国产| 女人高潮潮喷娇喘18禁视频| 亚洲国产av新网站| 国产爽快片一区二区三区| 成年美女黄网站色视频大全免费| 肉色欧美久久久久久久蜜桃| 久久99热这里只频精品6学生| 久久精品aⅴ一区二区三区四区| 亚洲免费av在线视频| 国产一级毛片在线| 大香蕉久久成人网| 欧美成人午夜精品| 久久天躁狠狠躁夜夜2o2o | 黄色一级大片看看| 午夜免费鲁丝| 高清视频免费观看一区二区| 国产日韩一区二区三区精品不卡| 考比视频在线观看| www.精华液| av欧美777| 国产男女内射视频| 男人添女人高潮全过程视频| 可以免费在线观看a视频的电影网站| 精品亚洲乱码少妇综合久久| 亚洲国产av新网站| 成人国产一区最新在线观看 | 天天躁夜夜躁狠狠躁躁| 在线观看免费午夜福利视频| 在线亚洲精品国产二区图片欧美| 精品熟女少妇八av免费久了| 欧美亚洲日本最大视频资源| 最黄视频免费看| tube8黄色片| 亚洲av电影在线观看一区二区三区| 成人亚洲精品一区在线观看| 欧美日韩av久久| 亚洲一卡2卡3卡4卡5卡精品中文| 免费少妇av软件| 国产成人av激情在线播放| 欧美日韩亚洲综合一区二区三区_| 少妇 在线观看| 国产极品粉嫩免费观看在线| 国产又爽黄色视频| xxxhd国产人妻xxx| 亚洲欧美一区二区三区黑人| 妹子高潮喷水视频| 久久久久视频综合| 欧美黑人欧美精品刺激| 我要看黄色一级片免费的| 国产精品国产三级专区第一集| 亚洲av在线观看美女高潮| 欧美国产精品一级二级三级| 成人亚洲精品一区在线观看| 亚洲欧美精品综合一区二区三区| 女警被强在线播放| 久热这里只有精品99| 久久国产亚洲av麻豆专区| 欧美在线一区亚洲| 三上悠亚av全集在线观看| 高清av免费在线| 国产欧美亚洲国产| a 毛片基地| 日本一区二区免费在线视频| 欧美黑人精品巨大| 男男h啪啪无遮挡| 久久久亚洲精品成人影院| 成年av动漫网址| 国产亚洲午夜精品一区二区久久| 看十八女毛片水多多多| 少妇 在线观看| 精品久久久久久电影网| 国产91精品成人一区二区三区 | 久久 成人 亚洲| 国产精品久久久久久人妻精品电影 | 色视频在线一区二区三区| 黄频高清免费视频| 黄色 视频免费看| 又粗又硬又长又爽又黄的视频| 日本黄色日本黄色录像| 成年人午夜在线观看视频| 精品熟女少妇八av免费久了| 欧美黑人精品巨大| 色播在线永久视频| 国产精品国产三级专区第一集| 精品人妻一区二区三区麻豆| 国产精品一区二区免费欧美 | 国产激情久久老熟女| 19禁男女啪啪无遮挡网站| 曰老女人黄片| 一区二区av电影网| 午夜福利一区二区在线看| 亚洲成人手机| 久久99精品国语久久久| 午夜影院在线不卡| 中文字幕av电影在线播放| 天天躁夜夜躁狠狠久久av| 欧美+亚洲+日韩+国产| 亚洲欧洲精品一区二区精品久久久| 国产av精品麻豆| 韩国高清视频一区二区三区| 亚洲国产欧美在线一区| 美女高潮到喷水免费观看| 狂野欧美激情性bbbbbb| 少妇人妻久久综合中文| 欧美黄色片欧美黄色片| 女性被躁到高潮视频| 18禁黄网站禁片午夜丰满| 国产精品久久久人人做人人爽| 99久久综合免费| 久久99精品国语久久久| 久久国产精品男人的天堂亚洲| 国精品久久久久久国模美| 亚洲激情五月婷婷啪啪| 国产野战对白在线观看| xxxhd国产人妻xxx| 久久久亚洲精品成人影院| 涩涩av久久男人的天堂| 日本av手机在线免费观看| av网站免费在线观看视频| av线在线观看网站| 两个人看的免费小视频| 高清黄色对白视频在线免费看| 久久人人97超碰香蕉20202| 欧美大码av| 亚洲国产精品一区二区三区在线| 国产爽快片一区二区三区| 国产人伦9x9x在线观看| 国产亚洲精品久久久久5区| 亚洲国产成人一精品久久久| 国产高清videossex| 久久人妻熟女aⅴ| 久久久久久人人人人人| 亚洲欧洲精品一区二区精品久久久| 中文字幕另类日韩欧美亚洲嫩草| 国产精品久久久久久精品古装| 日本猛色少妇xxxxx猛交久久| 大陆偷拍与自拍| 丝袜在线中文字幕| 男女高潮啪啪啪动态图| 久久人妻熟女aⅴ| 好男人电影高清在线观看| 久久热在线av| 交换朋友夫妻互换小说| 精品人妻一区二区三区麻豆| 脱女人内裤的视频| 亚洲一区二区三区欧美精品| 搡老岳熟女国产| 91精品国产国语对白视频| 久久精品成人免费网站| 亚洲国产毛片av蜜桃av| 中文字幕高清在线视频| 成年人黄色毛片网站| 欧美黑人精品巨大| 一本大道久久a久久精品|