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

    考慮水資源短缺風(fēng)險的不確定性水資源優(yōu)化配置

    2018-11-06 07:26:02王珊珊李天霄
    關(guān)鍵詞:降雨量水量生育

    付 強(qiáng),王珊珊,李天霄

    (東北農(nóng)業(yè)大學(xué)水利與土木工程學(xué)院,哈爾濱 150030)

    全球水資源需求以每年1%速度增長,未來20年還將繼續(xù)加快[1],關(guān)注水可獲得性、水質(zhì)和涉水風(fēng)險,改善水資源管理和利用可持續(xù)性。我國農(nóng)業(yè)灌溉用水短缺問題嚴(yán)峻,畝均水資源量僅占世界平均水平1/2[2-3]。水資源短缺突出,洪澇、干旱災(zāi)害與全球氣候變化變暖情況下,可利用水資源總量無法滿足農(nóng)業(yè)生產(chǎn)需求,農(nóng)業(yè)水資源供需矛盾加劇[4-5]。灌區(qū)作為農(nóng)業(yè)用水主體,優(yōu)化配置灌區(qū)農(nóng)業(yè)灌溉水資源,提高水資源綜合效益,保持水資源可持續(xù)利用具有重要意義。

    水資源優(yōu)化配置方法很多,灌區(qū)水資源管理系統(tǒng)受經(jīng)濟(jì)、資源等因素影響,為復(fù)雜系統(tǒng)動態(tài)配水過程。現(xiàn)實系統(tǒng)中水源降雨量時空分布不均勻性、可供水量、需水量、經(jīng)濟(jì)參數(shù)漲落不確定性,導(dǎo)致傳統(tǒng)確定性優(yōu)化方法無法準(zhǔn)確求解[6-7]。以往研究將水資源配置系統(tǒng)簡化為確定性系統(tǒng),忽略農(nóng)業(yè)水資源系統(tǒng)中不確定性,在水資源優(yōu)化配置上存在局限,決策方案使用價值有限[8-9]。學(xué)者針對系統(tǒng)中不確定性和復(fù)雜性,提出有效的系統(tǒng)分析方法。Chen等將隨機(jī)規(guī)劃框架中引入?yún)^(qū)間參數(shù),解決農(nóng)業(yè)用水系統(tǒng)中區(qū)間和概率分布不確定性[10];Li等將區(qū)間線性多目標(biāo)規(guī)劃應(yīng)用于農(nóng)業(yè)灌溉水資源配置[11];Xie等將魯棒性區(qū)間隨機(jī)規(guī)劃方法應(yīng)用于區(qū)域水資源管理[12];李晨洋等將區(qū)間兩階段模糊隨機(jī)規(guī)劃模型用于紅興隆灌區(qū)水資源管理[13];Yang等將區(qū)間參數(shù)線性規(guī)劃模型應(yīng)用于聯(lián)合使用地表水和地下水優(yōu)化調(diào)度[14],區(qū)間規(guī)劃不確定性理論方法可有效處理農(nóng)業(yè)系統(tǒng)復(fù)雜性問題。

    本文在前人研究基礎(chǔ)上,構(gòu)建區(qū)間線性模型,解決灌區(qū)水資源系統(tǒng)中存在的不確定性。降雨量和參考作物蒸騰量是灌溉系統(tǒng)中兩個具有相關(guān)關(guān)系的重要水文隨機(jī)變量,可表征灌區(qū)供需水變化關(guān)系[15-16]。由于二者變化規(guī)律存在相關(guān)性和隨機(jī)性,若僅分析單一變量變化規(guī)律和分布特征,不能反映灌區(qū)降雨量和參考作物蒸騰量的組合情況,無法科學(xué)指導(dǎo)灌溉規(guī)劃。以往研究很少將不同情景下水資源供需短缺風(fēng)險引入灌溉用水分配方案。采用Copula函數(shù)[17-18]對降雨量和參考作物蒸騰量二維聯(lián)合概率分布模型量化,分析灌區(qū)水資源供需短缺風(fēng)險組合。本文針對灌區(qū)水資源系統(tǒng)中不確定性,考慮供需關(guān)系,將區(qū)間線性規(guī)劃模型和水資源風(fēng)險短缺結(jié)合在各作物間優(yōu)化配水。模型中以灌區(qū)總收益最大為目標(biāo)函數(shù),采用Copula函數(shù)、離散區(qū)間解決不確定性問題,更適于灌區(qū)內(nèi)作物間優(yōu)化配置,合理優(yōu)化有限灌區(qū)農(nóng)業(yè)水資源,對提高黑龍江省灌區(qū)水資源利用效率具有重要意義。

    1 模型建立及求解

    1.1 Copula函數(shù)-聯(lián)合概率分布模擬

    假設(shè)(X1,X2,…,Xd)為d個隨機(jī)變量,其邊緣分布函數(shù)分別為FX1,F(xiàn)X2,…,F(xiàn)Xd?;赟klar定理Copula函數(shù)是多個隨機(jī)變量邊緣分布聯(lián)合分布函數(shù),定義域為[0,1],公式如下:

    式中,C()為Copula函數(shù);θ為Copula函數(shù)參數(shù)。

    將灌區(qū)降雨量和參考作物蒸騰量統(tǒng)計變化概率和概率分布特征結(jié)合,構(gòu)建水資源供需聯(lián)合概率分布模型?;诖?,引入Copula函數(shù)模擬其聯(lián)合分布。設(shè)X和Y分別為參考作物蒸騰量、降雨量隨機(jī)變量,其聯(lián)合分布函數(shù)為H(x,y),邊緣分布函數(shù)分別為Fx(x)和Fy(y),則存在聯(lián)合函數(shù)Copula C,滿足H(x,y)=C(Fx(x),F(xiàn)y(y))。邊緣分布模擬采用領(lǐng)域應(yīng)用較廣泛分布模型,即廣義極值(Generalized extreme value,GEV)和廣義正態(tài)(Generalized normal,GNO)。

    Copula函數(shù)各變量間相依結(jié)構(gòu)和邊緣分布可獨立考慮,應(yīng)用較優(yōu)優(yōu)勢。其中變量相關(guān)性函數(shù)有Frank、Clayton、Gumbel-Hougaard Copula 函 數(shù) 。通過采用非參數(shù)Kolmogorov-Simirnov(K-S)擬合度檢驗,離差平方和(OLS)最小準(zhǔn)則對Copula函數(shù)進(jìn)一步優(yōu)選,構(gòu)建和平灌區(qū)1958~2013年降雨量和參考作物蒸騰量聯(lián)合概率分布模型。

    采用頻率法劃分降雨量和參考作物蒸騰量豐枯。其中,

    1.2 區(qū)間線性規(guī)劃模型構(gòu)建

    灌區(qū)水資源配置目標(biāo)是將水資源合理配置到各種作物中,確保滿足作物生長所需最優(yōu)水量,獲得最大系統(tǒng)收益。但實際灌區(qū)研究中所獲得數(shù)據(jù)或信息具不確定性,將不確定相關(guān)參數(shù)值表達(dá)為一個區(qū)間范圍,可解決優(yōu)化模型中不確定性目標(biāo)函數(shù)和約束區(qū)間問題。通過區(qū)間規(guī)劃(Interval parameter programming,IPP)和線性模型(Linear programming,LP),構(gòu)建區(qū)間線性規(guī)劃模型(Interval parameter linear programming,IPLP)表達(dá)和求解。

    子模型下限:

    子模型上限:

    求解上限模型時考慮下限模型求解結(jié)果的約束,得到上限目標(biāo)函數(shù)決策目標(biāo)值和上限決策

    1.3 模型求解

    在優(yōu)化配置中,需考慮水源可供水量、灌溉需求、作物需水量等約束,考慮不同量級降雨量和參考作物蒸騰量遭遇組合,構(gòu)建區(qū)間線性優(yōu)化配置模型。

    模型目標(biāo)函數(shù)和約束表達(dá)如下:

    約束條件:

    ①水源可供水量約束

    ②灌溉需求約束

    ③作物需水量約束

    ④非負(fù)約束

    式中, f±為灌溉凈效益(元); j為不同作物,其中 j=1,2,3,分別代表水稻、玉米、大豆;為作物市場單價(元·kg-1);為灌區(qū)作物j種植面積(hm2);YWj為灌區(qū)作物 j水分生產(chǎn)率(kg·m-3);CW±為單方水價(元·m-3);為供水目標(biāo)(m3);為灌區(qū)初期蓄水量(m3);為灌區(qū)可利用來水量(m3); D±i為灌區(qū)蒸發(fā)、滲漏等損失水量(m3); EPj為灌區(qū)有效降雨量,m3·hm-2; Ri為灌區(qū)不同情景下灌溉水量(m3·hm-2);為灌區(qū)內(nèi)作物 j充分灌溉條件下最小需水量(m3);灌區(qū)內(nèi)作物 j充分灌溉條件下最大需水量(m3);η為灌溉水利用效率。

    1.4 作物灌溉水量優(yōu)化

    水資源有限情況下,一旦供水量未達(dá)到供水目標(biāo)值,需減少灌溉水量,但影響作物產(chǎn)量。因此,根據(jù)作物不同生育階段缺水對產(chǎn)量影響程度,判斷各生育階段需水程度,是否處于需水關(guān)鍵期,在生育階段應(yīng)用灌溉優(yōu)化模式[19-20]采用水分生產(chǎn)函數(shù)每個生育階段作物實際蒸騰量與作物最大蒸騰量比值,使產(chǎn)量最大化[21],保證灌溉水量在作物各生育階段合理配置,其中約束條件為每個生育階段

    2 案例研究

    2.1 研究區(qū)域概況

    黑龍江省和平灌區(qū)位于黑龍江省中部綏化市慶安縣境內(nèi),東經(jīng) 127°20'~127°49',北緯 46°41'~47°04'。地區(qū)氣候?qū)贉貛О敫珊蛋霛駶櫞箨懶约撅L(fēng)氣候,年平均氣溫1.7℃,多年平均蒸發(fā)量664.5 mm,多年平均降水量545.3 mm,主要降雨量集中在6~9月份,無霜期120~140 d。土壤以黑土和水稻土為主,灌區(qū)位居呼蘭河流域上游左岸一級階梯,灌區(qū)范圍由東向西呈帶狀分布,主要作物為水稻、玉米、大豆。和平灌區(qū)主要取水工程包含和平渠首和安邦河渠首。

    2.2 數(shù)據(jù)來源

    根據(jù)《呼蘭河灌區(qū)工程初期設(shè)計報告》及當(dāng)?shù)厮畡?wù)局提供調(diào)研數(shù)據(jù),得到不同作物(水稻、玉米、大豆)種植面積。通過黑龍江省用水定額標(biāo)準(zhǔn)及相關(guān)文獻(xiàn)田間試驗數(shù)據(jù)[22-24],得到3種作物單位面積最大和最小灌溉用水量(見表1)。

    基本氣象數(shù)據(jù)(降水量、平均風(fēng)速、最高氣溫、最低氣溫、平均相對濕度、日照時間等)來自中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)CMDC(http://data.cma.cn),利用FAO推薦Penman-Monteith公式[25]計算ET0,采用作物系數(shù)法[26]計算作物需水量。本文參照文獻(xiàn)[15]確定水稻、玉米、大豆水分生產(chǎn)率分別為1.61、0.82、0.72;灌溉水利用效率為0.4408。

    表1 和平灌區(qū)作物參數(shù)Table 1 Crops economic parameters in Heping irrigation district

    3 結(jié)果與分析

    3.1 邊緣分布函數(shù)模擬

    由圖1可知,年參考作物蒸騰量變化較小,而年降雨量與ET0間變化差異顯著。因此,對于灌溉水資源確定不同組合遭遇情景條件下ET0和降水有意義。應(yīng)用水文頻率分析方法,選用相應(yīng)隨機(jī)變量概率分布函數(shù),和平灌區(qū)降雨量服從正態(tài)分布,參考作物蒸騰量服從廣義極值分布。通過擬合結(jié)果判斷ET0和降水量理論分布函數(shù)與經(jīng)驗分布函數(shù)大致相同(見圖2),可假定和平灌區(qū)降雨量與ET0服從累積概率密度分布模型合理。

    借助變量相關(guān)性指標(biāo)分析水資源供需隨機(jī)變量間相關(guān)關(guān)系,采用Pearson相關(guān)系數(shù)γ、Kendall秩相關(guān)系數(shù)τ、Spearman秩相關(guān)系數(shù)ρ度量和平灌區(qū)1958~2013年降水量和參考作物蒸騰量間相關(guān)性,得出γ=-0.108、τ=-0.204、ρ=-0.146,模擬分析得到降雨量及參考作物蒸騰量邊緣分布函數(shù),其邊緣分布函數(shù)參數(shù)值見表3。可知和平灌區(qū)降雨量與ET0間存在負(fù)相關(guān)。因此,可采用Copula函數(shù)構(gòu)建水資源供需聯(lián)合分布模型描述和平灌區(qū)水資源供需相關(guān)特征。由Copula函數(shù)參數(shù)和Kendall秩相關(guān)系數(shù)間關(guān)系式,基于此運用Copula函數(shù)建立二維聯(lián)合分布模型,分別采用常見二維對稱Archimedean型(Frank、Clayton、Gumbel-Hougaard Copula函數(shù))構(gòu)建和平灌區(qū)1958~2013年降雨量和參考作物蒸騰量聯(lián)合概率分布模型。采用Kolmogorov-Simirnov(K-S)檢驗統(tǒng)計量D值均小于臨界值0.20056,因此3種Copula函數(shù)均可構(gòu)建聯(lián)合概率分布函數(shù)模型。通過離差平方和(OLS最小準(zhǔn)則)檢驗Copula函數(shù)擬合度,分別為0.1723、0.1931、0.1835。其中Frank Copula函數(shù)建立聯(lián)合分布模型下年降水量和參考作物蒸騰量聯(lián)合分布情況擬合較好,較好模擬降水量與ET0間相依關(guān)系??紤]灌區(qū)降雨量和參考作物蒸騰量不同組合條件下聯(lián)合分布概率:豐枯劃分為降雨量超過某一特定值663.5 mm、不超過某一特定值571.3 mm,其邊緣分布概率值分別為0.6209、0.3728;ET0超過某一特定值914 mm,不超過某一特定值852 mm,其邊緣分布概率值分別為0.6395、0.3615。采用頻率法分析和平灌區(qū)降雨量和ET0豐枯劃分情況頻率,豐枯遭遇情況分為9種。其中豐枯同步類型:P-ET0同豐,P-ET0同平,P-ET0同枯;分別為8.88%,7.54%,8.31%。豐枯異步類型:P-ET0豐平,PET0豐枯,P-ET0平枯,P-ET0平豐,P-ET0枯豐,P-ET0枯平;分別為9.13%,19.74%,8.33%,8.64%,19.92%,9.51%。降雨量和ET0豐枯異步頻率遠(yuǎn)大于其豐枯同步頻率,分別為75.27%,24.73%。P-ET0豐枯和P-ET0枯豐情況是9種遭遇組合中最不利于灌溉調(diào)度情況,其遭遇頻率約達(dá)20%。因此在農(nóng)業(yè)生產(chǎn)中,灌溉活動、科學(xué)調(diào)配水資源尤為必要。

    圖1 和平灌區(qū)1958~2013年降雨量、參考作物蒸騰量數(shù)據(jù)Fig.1 Long serial changes of precipitation and ET0for 1958-2013 years

    圖2 和平灌區(qū)降雨量、參考作物蒸騰量分布函數(shù)擬合Fig.2 Estimation of distribution functions of precipitation and ET0

    表2 水文變量基本統(tǒng)計參數(shù)Table 2 Basic statistical parameters of hydrological variables (mm)

    表3 邊緣分布函數(shù)參數(shù)值Table 3 Parameter values of marginal distribution function (mm)

    3.2 最優(yōu)配水方案

    由圖3可知,豐枯同步情境下,最優(yōu)總配水量變化較明顯,在P-ET0同枯8.31%,P-ET0同平7.54%,P-ET0同豐8.88%頻率下分別為[2 980,3 083]×104m3,[3 676,3 920]×104m3,[4 875,5 140]×104m3。其中,在P-ET0同豐可利用水量較大情景下,模型通過增大每種作物單位面積配水獲得目標(biāo)最優(yōu)化,總水量分配(凈灌溉水量和有效降雨量之和)達(dá)到作物最大需水量上限。而P-ET0同枯情景下為滿足作物最低需水量,確保作物基本生理需求,選擇加大經(jīng)濟(jì)效益獲得目標(biāo)最優(yōu)化,存在補(bǔ)給灌溉水量比例高達(dá)70.5%。

    豐枯異步情境下,頻率最高兩種遭遇組合P-ET0豐枯19.74%,P-ET0枯豐19.92%最優(yōu)總配水量分別為[4 470,4 462]×104m3,[3 420,3 785]×104m3。P-ET0平枯8.33%,P-ET0枯平9.51%頻率下最優(yōu)配水量分別為[4 602,4 754]×104m3,[3 874,3 990]×104m3。P-ET0平豐 8.64%,P-ET0豐平9.13%頻率下最優(yōu)配水量分別為[3 318,3 625]×104m3,[3 586,3 860]×104m3。表明不同情境下有效利用有限水資源分配方案重要性。

    整體看,P-ET0同枯,P-ET0枯豐,P-ET0枯平情境下,受最小需求量約束,模型將水量較為均勻分配。但P-ET0同豐,P-ET0豐平,P-ET0豐枯,總配水量相差有較大幅度差距,各類作物水量分配分化較明顯,實現(xiàn)模型目標(biāo)有利于作物獲得水量較多。因為降雨量具有明顯干濕特征,對作物水量配置貢獻(xiàn)影響較大,在降雨量較大濕潤時作物產(chǎn)量高于正常情況下,正常降雨下產(chǎn)量高于干旱條件。結(jié)果表明,作物對降水較敏感,通過可利用有效水資源應(yīng)對不同降雨量和參考作物蒸騰量聯(lián)合分布水資源短缺情況。

    圖3 不同情境下最優(yōu)水量分配及降雨量Fig.3 Optimal water allocation schemes and precipitation under different scenarios

    由圖4可知,由于市場價格和水分生產(chǎn)率相對較高,作物種植面積及配水量為水稻>玉米>大豆。分析所有遭遇組合頻率最大P-ET0豐枯和PET0枯豐情境對各作物水量分配,水量充足P-ET0豐枯19.74%情景下水稻、玉米、大豆配水量分別為 [3 702, 3 821]× 104m3, [359, 396] × 104m3,[409,443×104m3],對于模型結(jié)果較為敏感,水量區(qū)間分配變化幅度較大。

    各作物分配水量平均值比水量不足條件下PET0枯豐19.92%情景下高794×104m3,110×104m3,63×104m3。由于水稻種植面積過大,水稻配水必須滿足最低要求,而面積占比較小的玉米和大豆,少量水量即可滿足目標(biāo)優(yōu)化??蛇m當(dāng)調(diào)整農(nóng)作物種植結(jié)構(gòu),增加旱作物玉米、大豆種植面積,減少水稻種植面積,有效減少缺水量,提高水資源利用率。

    圖4 不同情境下作物最優(yōu)水量分配Fig.4 Optimal water allocation for each crop under different scenarios

    3.3 作物配水量優(yōu)化

    以耗水量最大水稻為典型作物,研究生育階段水量優(yōu)化情況,為便于定量分析水稻生育期降水量和需水量,將和平灌區(qū)水稻生育期確定為5月17日~9月20日,共127 d。其中返青期(5月17日~5月30日)、分蘗期(5月31日~7月7日)、拔節(jié)孕穗期(7月8日~7月25日)、抽穗開花期(7月26日~8月4日)、乳熟期(8月5日~8月24日)、黃熟期(8月25日~9月20日)。由于返青期稻苗較小、黃熟期自然落干,水稻蒸騰量較小,因此僅選取分蘗期、拔節(jié)期、抽穗期、乳熟期。采用相關(guān)文獻(xiàn)[21]分析各生育階段作物系數(shù),本文對所有情境中3種極端遭遇組合P-ET0同豐,P-ET0同平,P-ET0同枯計算水稻分蘗期、拔節(jié)期、抽穗期、乳熟期4個生育階段配水。在3種典型極端情境下滿足各生育期凈灌溉需水量時即達(dá)到ETa/ETm=1,作物從開始播種到成熟未遭受水分脅迫,產(chǎn)量達(dá)最高。如果作物生長發(fā)育階段,當(dāng)灌溉水量可用性達(dá)到限制時,產(chǎn)量下降。為避免水資源限制,考慮4個總體缺水目標(biāo)ETa/ETm(0.9,0.8,0.7,0.6),得到不同極端情境水分虧缺狀態(tài)時每個生育階段作物實際蒸騰量(ETa)和最大蒸騰量(ETm)間比值。

    由圖5可知,分蘗期、拔節(jié)期虧缺灌溉水平比例均低于抽穗期、乳熟期,因兩生育階段水分敏感指數(shù)高于其他,處于需水關(guān)鍵期。為獲得更高產(chǎn)量,分蘗期、拔節(jié)期分配水量顯著高于抽穗期、乳熟期。

    在P-ET0同豐情境下,ETa/ETm總體目標(biāo)為0.9、0.8時,分蘗期ETa/ETm為1.0,避免缺水灌溉。另外兩種目標(biāo)中,分蘗期、拔節(jié)期與其他生育階段ETa/ETm相比較大,抽穗期、乳熟期生育階段經(jīng)歷虧缺灌溉。在P-ET0同平和P-ET0同枯條件下,4個總體ETa/ETm目標(biāo)情況下,分蘗期無法達(dá)到充分灌溉。當(dāng)處于極度缺水條件ETa/ETm總體目標(biāo)為0.7、0.6時,由于3種情境下來水量年內(nèi)分配不均,在4個生育階段間水量分配間余水量供下一階段來水量較小生育階段,保證生育階段實際蒸騰量與最大蒸騰量比值(ETa/ETm)高于0.5。抽穗期、乳熟期ETa/ETm接近0.5,說明抽穗期、乳熟期是應(yīng)用虧缺灌溉主要生育階段,缺水水平較高。上述優(yōu)化結(jié)果表明,對于不同水分虧缺情景,分蘗期、拔節(jié)期屬于作物需水關(guān)鍵期,作物產(chǎn)量較大,因此為達(dá)高產(chǎn)目的,選擇較充分滿足作物水分基本需求,而在需水非關(guān)鍵期(抽穗期、乳熟期),為提高灌溉水利用效率,節(jié)約水資源,則選擇減少灌溉水量。P-ET0同枯,P-ET0同平,P-ET0同豐情景下配水量見圖6。

    圖5 3個典型情景4個缺水水平下水稻生育階段ETa/ETmFig.5 Optimized stage deficit targets(ETa/ETm)for four overall deficit targets for rice under different growth stages and scenarios

    圖6 3個典型情景4個缺水水平下水稻生育階段配水量變化情況Fig.6 Water amount change for four overall deficit targets for rice under different growth stages and scenarios

    3種情境下,P-ET0同豐配水量最大,P-ET0同枯反之;通過4個缺水水平對比,總體缺水目標(biāo)ETa/ETm為0.9時配水量最大。P-ET0同豐情境下供水目標(biāo)未達(dá)到初始目標(biāo)情況下缺水水平為0.9、0.8時,其中配水量分別為2 640.5×104m3、2 422.5×104m3,略下降,分蘗期、拔節(jié)期水量穩(wěn)定,抽穗期、乳熟期水量變化顯著。說明在保證產(chǎn)量同時,通過非工程措施應(yīng)用調(diào)虧灌溉優(yōu)化模式,管理者可根據(jù)不同缺水情境下判斷各生育階段實際蒸騰量與最大蒸騰量比值(ETa/ETm),使灌溉水資源在作物各生育階段合理配置,有效避免水資源浪費。而在其他情境缺水水平為0.7、0.6時配水量有顯著變化,產(chǎn)量均大幅下降。表明對于缺水水平較嚴(yán)重、模型穩(wěn)定性較差情況應(yīng)盡量避免。不同缺水條件下管理者應(yīng)優(yōu)先滿足作物分蘗期、拔節(jié)期用水需求,保證產(chǎn)量。

    4 討論

    本文通過建立區(qū)間不確定性模型解決水資源系統(tǒng)中的不確定性,考慮不同組合聯(lián)合概率,研究和平灌區(qū)作物水資源配置。以往水資源配置重點研究3種水平情境下配水方案,本文重點討論降雨量和ET09種量級組合下,應(yīng)用Copula函數(shù)定量分析灌區(qū)天然供水條件下水資源供需特征變量發(fā)生各種組合風(fēng)險,為灌區(qū)水資源短缺風(fēng)險分析提供技術(shù)指導(dǎo)。肖圓圓等僅研究高中低來水流量下水稻用水優(yōu)化[27],本文選擇3種作物不同條件展開研究,可更精確制定和調(diào)整灌溉規(guī)劃。通過將Stewart模型與水稻不同生育階段結(jié)合,減少灌溉水量達(dá)到較高產(chǎn)量。劉銀鳳研究認(rèn)為,水稻需水敏感期,應(yīng)充分滿足分蘗期、拔節(jié)期生育階段需水量,而需水非敏感期對產(chǎn)量影響較小,可選擇性減少灌溉水量,與本研究結(jié)果一致[28]。本研究假設(shè)水源同時向3種作物配水,與生產(chǎn)實際存在差異,對模型預(yù)測精度有一定影響。今后應(yīng)將環(huán)境效益、灌溉方式引入灌區(qū)作物水資源配置方案研究,使水資源優(yōu)化配置更符合實際,更具代表性。

    5 結(jié)論

    a.以黑龍江省和平灌區(qū)1958~2013年降水量P和ET0為研究對象,借助相關(guān)性指標(biāo)分析水資源隨機(jī)變量間相關(guān)性,兩者間存在負(fù)相關(guān)。根據(jù)邊緣分布函數(shù)極大似然法估計參數(shù),使用Copula函數(shù)構(gòu)建灌區(qū)自然供水條件下水資源聯(lián)合分布模型,應(yīng)用該模型分析相應(yīng)短缺風(fēng)險遭遇組合聯(lián)合概率。其中Frank Copula函數(shù)可較好擬合和平灌區(qū)降水量P和ET0理論聯(lián)合分布與經(jīng)驗聯(lián)合分布,運用聯(lián)合概率分布模型可較好描述不同遭遇組合事件,定量分析灌區(qū)水資源供需特征變量發(fā)生各種遭遇組合風(fēng)險情境。其中P-ET0豐枯和P-ET0枯豐情況是9種遭遇組合中最不利于灌溉調(diào)度情況,其遭遇頻率約達(dá)20%。

    b.用Copula函數(shù)及離散區(qū)間表示水資源系統(tǒng)不確定性,采用區(qū)間線性規(guī)劃模型優(yōu)化配置各種作物水量分配,確保滿足作物生長所需最優(yōu)水量,獲得最大系統(tǒng)收益。P-ET0同豐,P-ET0豐平,P-ET0豐枯總配水量較大,各類作物水量分配分化較明顯。從各種作物分配水量來看,耗水作物種植面積最大,水稻配水量最高,水量區(qū)間分配變化幅度較大。

    c.以水稻為典型作物優(yōu)化研究各生育階段用水量,考慮4個總體缺水目標(biāo)ETa/ETm(0.9,0.8,0.7,0.6),計算3種極端遭遇組合(P-ET0同豐,P-ET0同平,P-ET0同枯)下水稻分蘗期、拔節(jié)期、抽穗期、乳熟期階段配水,其中分蘗期、拔節(jié)期虧缺灌溉水平比例低于抽穗期、乳熟期,表明分蘗期、拔節(jié)期為需水關(guān)鍵期。在一定缺水條件,可據(jù)此模型在滿足產(chǎn)量同時對各生育階段灌溉水量做出調(diào)整,以高效利用水資源。

    猜你喜歡
    降雨量水量生育
    小水量超純水制備系統(tǒng)的最佳工藝選擇
    降雨量與面積的關(guān)系
    決不允許虐待不能生育的婦女
    紅土地(2018年12期)2018-04-29 09:16:40
    基于水力壓裂鉆孔的注水量及壓裂半徑的應(yīng)用研究
    應(yīng)對生育潮需早做準(zhǔn)備
    不能生育導(dǎo)致家庭破裂
    洞庭湖區(qū)降雨特性分析
    分散藥包千噸注水量的水壓爆破
    醫(yī)改和生育 兩大重點有看頭
    壓縮機(jī)級間冷卻冷凝析水量計算
    亚洲精品日韩av片在线观看| 嘟嘟电影网在线观看| av黄色大香蕉| 亚洲成人av在线免费| 蜜臀久久99精品久久宅男| 久久久久久伊人网av| 国精品久久久久久国模美| 下体分泌物呈黄色| 久久久a久久爽久久v久久| 国产精品99久久久久久久久| 国产精品 国内视频| 视频区图区小说| 青青草视频在线视频观看| 97精品久久久久久久久久精品| 国产精品麻豆人妻色哟哟久久| 国产精品.久久久| 街头女战士在线观看网站| 欧美国产精品一级二级三级| 最近的中文字幕免费完整| 丝袜喷水一区| 欧美日韩综合久久久久久| 3wmmmm亚洲av在线观看| 黄色毛片三级朝国网站| 丝袜在线中文字幕| 黄色配什么色好看| 美女主播在线视频| 在线观看美女被高潮喷水网站| 女性被躁到高潮视频| 久久久a久久爽久久v久久| 亚洲人成77777在线视频| 日本色播在线视频| 久久99热这里只频精品6学生| 少妇 在线观看| 男人添女人高潮全过程视频| 在线看a的网站| 亚洲av.av天堂| 热re99久久国产66热| 欧美成人精品欧美一级黄| 日韩成人伦理影院| 成人手机av| 日韩强制内射视频| 国产乱来视频区| 秋霞伦理黄片| 汤姆久久久久久久影院中文字幕| 18禁在线播放成人免费| 亚洲精品456在线播放app| a级片在线免费高清观看视频| 久久热精品热| 久久午夜综合久久蜜桃| 亚洲av综合色区一区| 国产伦理片在线播放av一区| 只有这里有精品99| 日本爱情动作片www.在线观看| 久久99热6这里只有精品| 赤兔流量卡办理| 91久久精品电影网| 国产视频首页在线观看| 热re99久久精品国产66热6| 亚洲欧美成人精品一区二区| 国产精品一国产av| 啦啦啦中文免费视频观看日本| 另类精品久久| 国产精品.久久久| 免费播放大片免费观看视频在线观看| 黄色毛片三级朝国网站| 精品亚洲成a人片在线观看| 国产精品一区二区在线观看99| 精品国产乱码久久久久久小说| 啦啦啦中文免费视频观看日本| 如何舔出高潮| 免费av不卡在线播放| 卡戴珊不雅视频在线播放| 国产欧美日韩一区二区三区在线 | 国产极品粉嫩免费观看在线 | 赤兔流量卡办理| 午夜日本视频在线| 午夜免费鲁丝| 高清毛片免费看| 久久久久久久久久久免费av| 黑人高潮一二区| 春色校园在线视频观看| 久久这里有精品视频免费| 边亲边吃奶的免费视频| av天堂久久9| 人妻系列 视频| 国产精品不卡视频一区二区| 欧美日韩精品成人综合77777| 观看美女的网站| 婷婷成人精品国产| 久久精品人人爽人人爽视色| 亚洲伊人久久精品综合| 国产片特级美女逼逼视频| 日韩视频在线欧美| 国产极品天堂在线| 啦啦啦中文免费视频观看日本| 18禁观看日本| 国产熟女欧美一区二区| 97在线视频观看| 午夜激情福利司机影院| 久久久久国产网址| 一级毛片 在线播放| 一本大道久久a久久精品| 成人综合一区亚洲| 亚洲精品久久久久久婷婷小说| 9色porny在线观看| 亚洲av在线观看美女高潮| 精品国产露脸久久av麻豆| 五月玫瑰六月丁香| 熟女人妻精品中文字幕| 国产探花极品一区二区| 亚洲av男天堂| 曰老女人黄片| 九草在线视频观看| 久久久久久久久大av| 免费观看的影片在线观看| 久久国产精品男人的天堂亚洲 | 精品久久久久久久久av| 另类精品久久| 伊人久久精品亚洲午夜| 美女国产视频在线观看| 在现免费观看毛片| 国产成人一区二区在线| 亚洲伊人久久精品综合| 丝袜脚勾引网站| 赤兔流量卡办理| 亚洲美女黄色视频免费看| 美女xxoo啪啪120秒动态图| 热99国产精品久久久久久7| 亚洲av不卡在线观看| 国产av一区二区精品久久| 国产av码专区亚洲av| 日韩中字成人| 天堂8中文在线网| 秋霞在线观看毛片| 十八禁高潮呻吟视频| 99re6热这里在线精品视频| 一级毛片黄色毛片免费观看视频| 亚洲av国产av综合av卡| 久久久久久久久久久久大奶| www.色视频.com| 99九九线精品视频在线观看视频| 少妇猛男粗大的猛烈进出视频| 国产精品99久久99久久久不卡 | 伦理电影大哥的女人| 乱码一卡2卡4卡精品| 国产一区有黄有色的免费视频| 97精品久久久久久久久久精品| 午夜免费男女啪啪视频观看| 国产av精品麻豆| 极品人妻少妇av视频| av国产久精品久网站免费入址| 内地一区二区视频在线| 国产精品一区二区三区四区免费观看| 亚洲av中文av极速乱| 精品久久蜜臀av无| 亚洲综合精品二区| 免费人成在线观看视频色| 亚洲精品日韩在线中文字幕| 美女xxoo啪啪120秒动态图| 51国产日韩欧美| 少妇的逼水好多| 汤姆久久久久久久影院中文字幕| 日本av手机在线免费观看| 国产69精品久久久久777片| 麻豆成人av视频| 亚洲色图 男人天堂 中文字幕 | 国产黄色免费在线视频| 大片电影免费在线观看免费| videossex国产| 乱人伦中国视频| 人妻制服诱惑在线中文字幕| 免费高清在线观看视频在线观看| 午夜免费鲁丝| 欧美亚洲日本最大视频资源| 尾随美女入室| 肉色欧美久久久久久久蜜桃| 18禁裸乳无遮挡动漫免费视频| 春色校园在线视频观看| 国产成人精品一,二区| 人人妻人人澡人人爽人人夜夜| 久久久久久久久久久久大奶| 色婷婷av一区二区三区视频| 国产片内射在线| 亚洲国产精品成人久久小说| 亚洲人成网站在线播| 精品人妻在线不人妻| 国产精品久久久久久久电影| 26uuu在线亚洲综合色| 久久人人爽人人爽人人片va| 精品卡一卡二卡四卡免费| 日韩一区二区三区影片| 亚洲三级黄色毛片| 国产男女超爽视频在线观看| 人人妻人人澡人人看| 色婷婷久久久亚洲欧美| 少妇高潮的动态图| 国产高清国产精品国产三级| 国产高清三级在线| 最黄视频免费看| 日本欧美国产在线视频| 国产成人精品福利久久| 久久ye,这里只有精品| 高清毛片免费看| 亚洲激情五月婷婷啪啪| 日韩欧美精品免费久久| 午夜福利,免费看| 十分钟在线观看高清视频www| 人妻系列 视频| 熟女电影av网| 国产免费视频播放在线视频| kizo精华| 色吧在线观看| 国产无遮挡羞羞视频在线观看| 99国产综合亚洲精品| 男女国产视频网站| 夜夜爽夜夜爽视频| 亚洲国产av新网站| 亚洲成人一二三区av| 在线天堂最新版资源| 成人漫画全彩无遮挡| 免费黄频网站在线观看国产| 视频区图区小说| 亚洲欧美日韩另类电影网站| 国产av精品麻豆| 99九九线精品视频在线观看视频| 国产淫语在线视频| 少妇人妻久久综合中文| 晚上一个人看的免费电影| 男女免费视频国产| 好男人视频免费观看在线| 在线播放无遮挡| 成人国产av品久久久| 欧美三级亚洲精品| 亚洲婷婷狠狠爱综合网| 亚洲国产日韩一区二区| 免费久久久久久久精品成人欧美视频 | 亚洲国产av新网站| 特大巨黑吊av在线直播| 视频在线观看一区二区三区| 在线 av 中文字幕| 大陆偷拍与自拍| 丰满迷人的少妇在线观看| 一个人免费看片子| 精品国产乱码久久久久久小说| 亚州av有码| 永久网站在线| 日本av免费视频播放| av国产精品久久久久影院| 亚洲第一av免费看| 国产成人精品久久久久久| 久久久久人妻精品一区果冻| 晚上一个人看的免费电影| 色网站视频免费| 国产精品国产三级国产av玫瑰| 我要看黄色一级片免费的| 日韩电影二区| 亚洲精品乱码久久久久久按摩| 日本av手机在线免费观看| 亚洲av在线观看美女高潮| 亚洲精品乱码久久久久久按摩| 国产成人a∨麻豆精品| 在线观看一区二区三区激情| 亚洲激情五月婷婷啪啪| 热re99久久国产66热| 秋霞伦理黄片| 欧美日韩在线观看h| 美女内射精品一级片tv| 久久国产精品男人的天堂亚洲 | 久久99一区二区三区| 国产精品三级大全| 99热全是精品| 女性生殖器流出的白浆| 九九久久精品国产亚洲av麻豆| 亚洲综合精品二区| 蜜桃在线观看..| 搡女人真爽免费视频火全软件| av电影中文网址| 两个人免费观看高清视频| 国产黄频视频在线观看| 国产男女内射视频| 丝袜脚勾引网站| 99久久精品一区二区三区| 精品人妻偷拍中文字幕| 草草在线视频免费看| 卡戴珊不雅视频在线播放| 精品一区二区三卡| 街头女战士在线观看网站| 久久人人爽av亚洲精品天堂| 国产亚洲一区二区精品| 亚洲,欧美,日韩| av福利片在线| 国产免费一区二区三区四区乱码| 亚洲少妇的诱惑av| 国产精品一区二区三区四区免费观看| 久久精品久久久久久噜噜老黄| 人妻人人澡人人爽人人| 在线观看一区二区三区激情| 一级毛片 在线播放| 狠狠精品人妻久久久久久综合| 国产极品天堂在线| 少妇熟女欧美另类| 国产极品天堂在线| 午夜激情福利司机影院| 一级,二级,三级黄色视频| av不卡在线播放| 看十八女毛片水多多多| 亚洲欧美日韩卡通动漫| 久久精品国产亚洲av涩爱| 国产精品欧美亚洲77777| 一级爰片在线观看| 色吧在线观看| 日本91视频免费播放| 国产极品粉嫩免费观看在线 | 永久免费av网站大全| av电影中文网址| 久久久国产精品麻豆| 校园人妻丝袜中文字幕| 人人妻人人爽人人添夜夜欢视频| 多毛熟女@视频| 精品一区二区三卡| 九色成人免费人妻av| 欧美精品高潮呻吟av久久| 日韩制服骚丝袜av| 国产在线一区二区三区精| 日产精品乱码卡一卡2卡三| 国语对白做爰xxxⅹ性视频网站| 亚洲精品一二三| 欧美一级a爱片免费观看看| 久久99一区二区三区| 观看美女的网站| 丝袜喷水一区| 制服诱惑二区| 国产在视频线精品| 久久亚洲国产成人精品v| 3wmmmm亚洲av在线观看| 国产在视频线精品| 中国美白少妇内射xxxbb| 日韩伦理黄色片| 欧美日韩av久久| 亚洲,欧美,日韩| 日日撸夜夜添| 999精品在线视频| 国产成人精品无人区| 日本爱情动作片www.在线观看| 国产精品麻豆人妻色哟哟久久| videosex国产| 婷婷色综合大香蕉| 毛片一级片免费看久久久久| 老女人水多毛片| 亚洲精品久久久久久婷婷小说| 91在线精品国自产拍蜜月| 在线观看人妻少妇| 亚洲激情五月婷婷啪啪| 一级毛片电影观看| 婷婷成人精品国产| 免费av不卡在线播放| 啦啦啦在线观看免费高清www| 最近2019中文字幕mv第一页| 99re6热这里在线精品视频| 少妇人妻久久综合中文| 一区二区三区免费毛片| 另类亚洲欧美激情| 婷婷色综合大香蕉| 欧美精品亚洲一区二区| 少妇高潮的动态图| 一区二区三区免费毛片| 日本免费在线观看一区| 色婷婷久久久亚洲欧美| 日韩一本色道免费dvd| 国产精品99久久99久久久不卡 | 亚洲精品国产av成人精品| 国产免费现黄频在线看| 国语对白做爰xxxⅹ性视频网站| 黑丝袜美女国产一区| 国模一区二区三区四区视频| 肉色欧美久久久久久久蜜桃| 久久久久网色| 22中文网久久字幕| 亚洲精品成人av观看孕妇| 久久免费观看电影| 亚洲五月色婷婷综合| 只有这里有精品99| 久久久久久久大尺度免费视频| 欧美bdsm另类| 亚洲不卡免费看| 最新中文字幕久久久久| 中文字幕亚洲精品专区| 草草在线视频免费看| 久久精品国产a三级三级三级| 免费播放大片免费观看视频在线观看| 久久久精品区二区三区| 日本-黄色视频高清免费观看| 国产精品三级大全| 99九九线精品视频在线观看视频| 在线观看一区二区三区激情| 国产午夜精品久久久久久一区二区三区| 我要看黄色一级片免费的| 亚洲经典国产精华液单| 汤姆久久久久久久影院中文字幕| 亚洲欧美清纯卡通| 亚洲久久久国产精品| 菩萨蛮人人尽说江南好唐韦庄| 嘟嘟电影网在线观看| 日本黄大片高清| 亚洲欧美日韩卡通动漫| 熟女av电影| 成人影院久久| 久久久久久人妻| 亚洲天堂av无毛| 水蜜桃什么品种好| 人人妻人人澡人人爽人人夜夜| 青青草视频在线视频观看| 在线亚洲精品国产二区图片欧美 | 91在线精品国自产拍蜜月| 亚洲精品乱码久久久v下载方式| 男女边摸边吃奶| 大片免费播放器 马上看| 日本欧美视频一区| 91成人精品电影| 9色porny在线观看| 91在线精品国自产拍蜜月| 久久99精品国语久久久| 狠狠婷婷综合久久久久久88av| 日本色播在线视频| 九色成人免费人妻av| 久久国产亚洲av麻豆专区| 最近的中文字幕免费完整| 97在线人人人人妻| 国产不卡av网站在线观看| 插逼视频在线观看| 精品久久久久久电影网| 成人漫画全彩无遮挡| 最近的中文字幕免费完整| 国产日韩欧美在线精品| 夫妻性生交免费视频一级片| 91aial.com中文字幕在线观看| 日韩中文字幕视频在线看片| 美女内射精品一级片tv| 999精品在线视频| h视频一区二区三区| 黑丝袜美女国产一区| 26uuu在线亚洲综合色| 久久久精品94久久精品| 精品久久久精品久久久| 成年av动漫网址| 欧美日本中文国产一区发布| 满18在线观看网站| 又大又黄又爽视频免费| 99九九线精品视频在线观看视频| 春色校园在线视频观看| 亚洲国产成人一精品久久久| 久久狼人影院| 全区人妻精品视频| 五月伊人婷婷丁香| 又粗又硬又长又爽又黄的视频| 国产69精品久久久久777片| 如何舔出高潮| av有码第一页| 久久久午夜欧美精品| 亚洲av福利一区| 99久久综合免费| 日韩 亚洲 欧美在线| 伊人久久国产一区二区| 80岁老熟妇乱子伦牲交| 免费av不卡在线播放| 蜜臀久久99精品久久宅男| 菩萨蛮人人尽说江南好唐韦庄| 成人影院久久| 有码 亚洲区| 搡老乐熟女国产| 免费观看的影片在线观看| 国产精品久久久久久久久免| 国产精品麻豆人妻色哟哟久久| 两个人免费观看高清视频| 在线精品无人区一区二区三| 免费观看在线日韩| 久久99蜜桃精品久久| 亚洲欧洲日产国产| av国产精品久久久久影院| a级毛片黄视频| 亚洲av在线观看美女高潮| 免费看不卡的av| 精品国产露脸久久av麻豆| 91aial.com中文字幕在线观看| 99热网站在线观看| 在线观看国产h片| 一区二区三区乱码不卡18| 青春草国产在线视频| 新久久久久国产一级毛片| 在线播放无遮挡| av.在线天堂| 汤姆久久久久久久影院中文字幕| 国产成人91sexporn| 国产高清国产精品国产三级| 人妻人人澡人人爽人人| 高清毛片免费看| 高清视频免费观看一区二区| 久久99一区二区三区| 日韩精品有码人妻一区| 啦啦啦视频在线资源免费观看| 99久久人妻综合| 免费高清在线观看日韩| 亚洲伊人久久精品综合| 亚洲成人手机| 国产日韩欧美亚洲二区| 满18在线观看网站| 在线观看国产h片| 欧美 日韩 精品 国产| 亚洲中文av在线| 久久久久久伊人网av| 国产片特级美女逼逼视频| 青春草视频在线免费观看| 欧美精品人与动牲交sv欧美| av天堂久久9| 嘟嘟电影网在线观看| www.av在线官网国产| 日韩大片免费观看网站| 99久国产av精品国产电影| 男女无遮挡免费网站观看| 欧美日韩成人在线一区二区| 美女福利国产在线| 免费黄频网站在线观看国产| 99国产综合亚洲精品| 男女高潮啪啪啪动态图| 99热这里只有是精品在线观看| 人妻制服诱惑在线中文字幕| 成人国产av品久久久| 少妇的逼好多水| 免费黄色在线免费观看| 午夜视频国产福利| 久久久久精品久久久久真实原创| 亚洲婷婷狠狠爱综合网| 美女中出高潮动态图| 制服诱惑二区| 在线观看免费高清a一片| 久久久国产欧美日韩av| 精品人妻偷拍中文字幕| 免费观看性生交大片5| 亚洲图色成人| 亚洲av日韩在线播放| 国产精品欧美亚洲77777| 国产无遮挡羞羞视频在线观看| 免费看不卡的av| 少妇丰满av| 欧美 亚洲 国产 日韩一| 久久精品熟女亚洲av麻豆精品| 免费高清在线观看视频在线观看| 日韩成人av中文字幕在线观看| 91精品国产九色| 精品国产一区二区三区久久久樱花| 精品久久蜜臀av无| 久久久久网色| 国产伦理片在线播放av一区| 26uuu在线亚洲综合色| 一区二区三区免费毛片| 一个人免费看片子| 91国产中文字幕| 一级毛片我不卡| 亚洲国产日韩一区二区| 丰满乱子伦码专区| 中文字幕亚洲精品专区| 亚洲国产精品一区二区三区在线| 国产亚洲午夜精品一区二区久久| 欧美日韩亚洲高清精品| 亚洲国产精品一区二区三区在线| 综合色丁香网| 精品一区二区三卡| 视频区图区小说| 国产日韩欧美视频二区| 你懂的网址亚洲精品在线观看| 亚洲精品久久午夜乱码| 久久婷婷青草| 亚洲综合色惰| 国产日韩欧美视频二区| 九九在线视频观看精品| 日本猛色少妇xxxxx猛交久久| 高清毛片免费看| 亚洲色图 男人天堂 中文字幕 | 一级,二级,三级黄色视频| 久久国产亚洲av麻豆专区| 热99久久久久精品小说推荐| 色婷婷av一区二区三区视频| 曰老女人黄片| videosex国产| 欧美性感艳星| 欧美一级a爱片免费观看看| av在线观看视频网站免费| 国产 一区精品| 中国三级夫妇交换| 啦啦啦啦在线视频资源| 永久免费av网站大全| 国产成人免费无遮挡视频| 国产免费又黄又爽又色| 免费播放大片免费观看视频在线观看| 日本欧美国产在线视频| 国产av精品麻豆| 热re99久久国产66热| 性色av一级| 黑丝袜美女国产一区| 久久影院123| 日韩 亚洲 欧美在线| 国产成人freesex在线| 精品久久国产蜜桃| 国产日韩欧美在线精品| 看非洲黑人一级黄片| 久久精品国产亚洲网站| 午夜日本视频在线| 成人二区视频| 成人国产麻豆网| 一级爰片在线观看| 免费av中文字幕在线| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 一区二区三区免费毛片| 日本免费在线观看一区| 成年人午夜在线观看视频|