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

    基于貝葉斯理論的水文多模型預(yù)報(bào)

    2017-02-27 23:18:17周正葉愛中馬鳳杜超
    南水北調(diào)與水利科技 2017年1期
    關(guān)鍵詞:不確定性

    周正+葉愛中+馬鳳+杜超

    摘要:水文預(yù)報(bào)對于防洪、抗旱以及水資源調(diào)度等具有重要意義。水文預(yù)報(bào)通常依靠水文模型來完成,由于受到不同流域特點(diǎn)、產(chǎn)匯流機(jī)制等的限制,每個(gè)水文模型都具有各自的特點(diǎn)及適用區(qū)域。單一模型具有非常大的水文預(yù)報(bào)不確定性,為了解決單一模型局限性的問題,多模型水文預(yù)報(bào)常作為降低水文預(yù)報(bào)不確定性有效方法之一。選用三種常見的水文模型:時(shí)變增益水文模型、新安江模型和薩克拉門托模型,在珠江飛來峽流域進(jìn)行分布式建模,采用相同的輸入與初始場,三個(gè)模型獨(dú)立進(jìn)行模擬,然后對比三個(gè)模型的結(jié)果,并進(jìn)行貝葉斯多模型加權(quán)平均和簡單平均得到多模型平均結(jié)果,研究結(jié)果表明,貝葉斯模型處理后的結(jié)果要比單個(gè)模型模擬結(jié)果和簡單平均處理后的結(jié)果準(zhǔn)確率高。

    關(guān)鍵詞:飛來峽;水文預(yù)報(bào);不確定性;貝葉斯多模型

    中圖分類號:TV124文獻(xiàn)標(biāo)志碼:A文章編號:1672-1683(2017)01-0043-06

    Abstract:Hydrological forecast is very important for flood control,drought resistance and water resources regulation.It is usually based on model simulation.Each hydrological model has its own characteristics and feasible basin.Multi-model hydrological forecast is one of the effective methods to reduce the forecasting uncertainty.This paper chose three commonly used hydrological models:time-variant gain hydrological model,Xin′anjiang model,and Sacramento model.The case study was on Feilaixia basin of Pearl River.The three hydrological models were used for independent simulation with the same input and initial value.Then,BMA and SMA were run on the three models′results.The research results show that the BMA results are better than the results from single model simulation.

    Key words:Feilaixia;hydrological forecast;uncertainty;Bayesian model averagin

    水文預(yù)報(bào)不確定性主要來源于水文模型的不確定性。為了降低水文模型的不確定性,越來越多的水文模型被提出并應(yīng)用到不同流域[1]。其中應(yīng)用比較廣泛的有薩克拉門托模型(Sacramento model,SAC)[2],新安江模型(Xin An Jiang,XAJ)[3],和時(shí)變增益水文模型(Time Variation Gain Model,DTVGM)[4]等。然而每個(gè)模型都具有各自的優(yōu)勢與局限性,而且由于水文過程的復(fù)雜性,導(dǎo)致水文模型不確定性日益凸顯。對此,廣大水文學(xué)者提出多模型信息融合的方法。該方法可以降低單個(gè)模型的偏差,使模擬更接近于觀測值。目前,水文預(yù)報(bào)綜合處理方法主要有簡單平均法(Simple Model Averaging,SMA)、加權(quán)平均法、神經(jīng)網(wǎng)絡(luò)和模糊推理法等[5-10]。其中,貝葉斯加權(quán)平均方法(Bayesian Model Averaging,BMA)最為廣泛使用。貝葉斯方法是基于貝葉斯理論而發(fā)展起來用于解決統(tǒng)計(jì)問題的方法。隨著科學(xué)技術(shù)的發(fā)展,貝葉斯方法在統(tǒng)計(jì)學(xué)[11]成功應(yīng)用之后,逐漸發(fā)展到計(jì)量經(jīng)濟(jì)學(xué)[12],社會學(xué)[13]等方面。在氣象領(lǐng)域,貝葉斯最大似然方法被應(yīng)用于預(yù)測溫度及降水的概率分布[14],Ma[15]等人通過改變BMA目標(biāo)函數(shù),將其應(yīng)用于NMME季節(jié)性降水集合預(yù)報(bào)。近年來,BMA方法被應(yīng)用于水文領(lǐng)域:Neuman[16]提出利用最大似然貝葉斯平均方法分析地下水模型不確定性來優(yōu)化水文預(yù)報(bào)的方法,Kavetski[17]等人通過使用貝葉斯誤差分析方法評估輸入降雨不確定性對VIC模型模擬的影響。Zhang[18]利用BMA方法對SWAT水文模型的擬定和不確定性分析。為了得到對水文模型更為優(yōu)化的處理結(jié)果,Raftery[19]等采用BMA對多個(gè)數(shù)字天氣預(yù)報(bào)模型權(quán)重進(jìn)行最優(yōu)化處理,Duan[20]通過貝葉斯方法進(jìn)行多模型水文集合預(yù)報(bào),并得到比較好的結(jié)果。

    基于BMA處理后的高精度預(yù)報(bào)[21],本文將BMA方法應(yīng)用于飛來峽流域進(jìn)行模擬預(yù)報(bào),并與DTVGM模型、XAJ模型、SAC模型及SMA方法處理結(jié)果進(jìn)行對比,分析多模型方法的預(yù)報(bào)結(jié)果。

    1 流域概況

    飛來峽流域位于珠江第二大水系北江的中上游,控制集水面積為34 097 km2,大約為北江流域面積的73%。該流域?qū)儆跓釒Ш蛠啛釒夂騾^(qū),且由于東南和西南季風(fēng)的影響,降水不均,汛期平均降水量為1 224 mm,約占全年降水量的70%[22-23],因此時(shí)常會遭受洪水災(zāi)害。流域概況見圖1。

    本文以飛來峽流域?yàn)檠芯繉ο螅捎?00 km2河網(wǎng)閾值將流域劃分為146個(gè)子流域。收集到氣象數(shù)據(jù)為1980年-2009年日降水、氣溫和蒸發(fā)數(shù)據(jù)。流量采用流域出口1980年-2009年日流量數(shù)據(jù),時(shí)間序列長度為30年,其中1980年-1999年作為擬定期,2000年-2009年作為驗(yàn)證期。

    2 模型方法介紹

    2.1 水文模型介紹

    2.1.1 分布式時(shí)變增益水文模型

    時(shí)變增益模型是夏軍[24-25]于1989年-1995年期間提出的一種研究方法,建立了降水-土壤濕度-產(chǎn)流之間的非線性關(guān)系。隨著GIS技術(shù)的發(fā)展,時(shí)變增益模型發(fā)展為分布式時(shí)變增益水文模型,并且仍在不斷地改進(jìn)。該模型可以根據(jù)時(shí)間尺度分為月尺度、日尺度和小時(shí)尺度。

    2.1.2 新安江模型

    新安江模型是于20世紀(jì)60、70年代由我國水文工作者研制開發(fā)的水文模型,主要包括蒸發(fā)計(jì)算模型、產(chǎn)流模型、三水源劃分模型以及匯流模型。其理論主要基于蓄滿產(chǎn)流理念,即土壤濕度達(dá)到田間持水量以后才開始進(jìn)行產(chǎn)流,在未達(dá)到田間持水量時(shí),則不進(jìn)行產(chǎn)流,所有的降水都會被土壤吸收,該模型被廣泛的應(yīng)用在濕潤與半濕潤的中型或者大型流域[26-28]。

    2.1.3 薩克拉門托模型

    薩克拉門托流域水文模型是集總參數(shù)型、連續(xù)運(yùn)算的確定性流域水文模型,由美國天氣局水文辦公室薩克拉門托預(yù)報(bào)中心于20世紀(jì)70年代初期,于第IV號斯坦福模型基礎(chǔ)上改進(jìn)和發(fā)展。薩克門托模型的核心主要為土壤含水量模型,結(jié)構(gòu)分為上、下兩層,通過計(jì)算輸入降水模擬蒸散發(fā)量、地表徑流以及下滲水量[29-31]。

    2.2 BMA和SMA介紹

    BMA (Bayesian Model Averaging,BMA)是一個(gè)通過加權(quán)平均不同模型的預(yù)報(bào)值得到更可靠的綜合預(yù)報(bào)值的數(shù)學(xué)方法[32]。貝葉斯所選定的模型空間中,最佳模型并不確定,根據(jù)貝葉斯模型平均法基本理論,可以給出模型變量的概率密度函數(shù),通過對實(shí)測和模擬流量數(shù)據(jù)的正態(tài)分位變化,進(jìn)行期望最大化(EM)算法,進(jìn)而通過正態(tài)反轉(zhuǎn)換得到空間對應(yīng)關(guān)系,即對不同模型所占權(quán)重進(jìn)行加權(quán)平均,從而計(jì)算出BMA模擬流量[33]。

    簡單平均法(Simple Model Averaging,SMA)即為模型模擬流量的簡單平均,每個(gè)模型所占的權(quán)重相同。

    2.3 模型驗(yàn)證指標(biāo)

    為了對模型模擬結(jié)果進(jìn)行合理的評價(jià),采用相關(guān)系數(shù)(R)、納什效率系數(shù)(NSE、NSEB[35]、NSEI[36])、偏差(rBias)和均方根誤差(RMSE)評價(jià)指標(biāo)(見表1)。其中rBias接近于0說明模型模擬總量和觀測總量一致,NSE越接近于1說明模型整體可靠性越高。

    由于NSE受模擬洪峰影響較大,即NSE的高低由模擬洪峰的準(zhǔn)確度主導(dǎo):洪峰模擬越好,則NSE越高,而枯水期的模擬對NSE的影響較小,因此本文增加NSEI指標(biāo),即NSEI越接近于1則枯水流量模擬精度越高。通過NSE和NSEI兩個(gè)指標(biāo)可以判斷整體模擬效果。此外,本文增加NSEB指標(biāo),NSEB越接近于1說明相對于參考系統(tǒng)(本文采用多年平均流量)更優(yōu)。

    3 結(jié)果分析與討論

    3.1 模型擬定與驗(yàn)證

    本文選定1980年-1999年為擬定期,2000年-2009年為驗(yàn)證期,分別采用分布式時(shí)變增益水文模型、新安江模型以及薩克拉門托模型進(jìn)行模擬[34],在得到模擬流量之后分別用BMA和簡單算術(shù)平均對模型結(jié)果進(jìn)行后處理,從而判斷模型的模擬性能,結(jié)果見表2。

    在擬定期,三個(gè)模型模擬的相關(guān)系數(shù)均達(dá)到了0.9以上,納什效率系數(shù)在0.85以上,偏差均在±5%左右,表明三個(gè)模型較高的可靠性和較小的模擬偏差。

    在驗(yàn)證期,DTVGM和新安江模型模擬結(jié)果相關(guān)系數(shù)在0.9以上,納什效率系數(shù)在0.75以上,略優(yōu)于SAC模型,三個(gè)模型均表現(xiàn)出比較好的模擬結(jié)果說明了擬定期參數(shù)的合理性。倒數(shù)納什效率系數(shù)(NSEI)在擬定和驗(yàn)證期均不太高,尤其是SAC模型出現(xiàn)負(fù)值,說明水文模型對枯季流量模擬精度較差。

    圖2和圖3分別給出了擬定期和驗(yàn)證期幾個(gè)模型的日徑流模擬過程圖,可以看出,分布式時(shí)變增益水文模型、新安江模型以及薩克拉門托模型整體模擬效果比較好,洪峰及峰現(xiàn)時(shí)間同觀測一致。

    3.2 多模型比較分析

    從表2可看出BMA和SMA的結(jié)果相對優(yōu)于三個(gè)模型獨(dú)立的模擬結(jié)果,而BMA的結(jié)果要優(yōu)于簡單平均的處理結(jié)果。通過BMA權(quán)重分析,充分考慮三個(gè)模型的精度,DTVGM、新安江模型和SAC模型所占權(quán)重分別為0.729、0.252、0.019,其中DTVGM模型所占權(quán)重最大,說明該模型模擬結(jié)果要優(yōu)于其余兩個(gè)模型,符合表2指標(biāo)。SAC模型由于在枯水期流量模擬精度較差,所以權(quán)重較低。簡單算術(shù)平均方法三個(gè)模型的權(quán)重相等,雖然在一定程度上對最終模擬效果進(jìn)行了優(yōu)化,但是缺乏對模型信息最大化的利用,并未考慮每個(gè)模型各自的精度,因此處理效果比BMA差。

    為了更清楚的展示BMA和單個(gè)模型模擬結(jié)果,圖4做出了多個(gè)模型在驗(yàn)證期10年數(shù)據(jù)的日平均流量,由圖4可以發(fā)現(xiàn)BMA處理后的結(jié)果更接近于實(shí)測流量,尤其在低水期:1月1日-3月1日及9月15日-12月31日時(shí)間段內(nèi)BMA優(yōu)化效果更為明顯。

    4 小結(jié)

    本研究探討了多模型水文預(yù)報(bào)應(yīng)用于實(shí)際流域的可行性,飛來峽流域?qū)儆跐駶櫟貐^(qū),三個(gè)模型對洪水都有比較好的模擬結(jié)果,其中DTVGM模擬結(jié)果最好,其次為XAJ模型和SAC模型,但在枯水期模擬存在一定的誤差,經(jīng)BMA和SMA處理后的結(jié)果要優(yōu)于單個(gè)模型的模擬結(jié)果。BMA處理方法充分考慮了每個(gè)模型的精度,最大程度利用了模型信息,給出了每個(gè)模型所占權(quán)重,并具有較高相關(guān)系數(shù)和最小的均方根誤差,較之于簡單的算數(shù)平均方法,具有更高的精度和穩(wěn)定性,并且在枯水期模擬效果更為明顯。

    參考文獻(xiàn)(References):

    [1] 李懷恩,沈晉.現(xiàn)行幾個(gè)主要產(chǎn)流模型的剖析[J].水文,1996(1):14-23.(LI Huai-en,SHEN Jin.Analysis of the main runoff yield models[J].Hydrology,1996(1):14-23.(in Chinese))

    [2] Burnash R J C.Chapter 10:the NWS river forecast system-catchment modeling[M],in Computer Models of Water Shed Hydrology,edited by Singh V P,Water Resources Publications,Littleton Colo,1995.

    [3] Zhao R J.The Xinanjiang model applied in China[J].Journal of Hydrology,1992,135(1):371-381.DOI:10.1016/0022-1694(92)90096-E.

    [4] Xia J.A system approach to real time hydrological forecasts in watersheds[J].Water International,27(1):87-97.DOI:10.1080/02508060208686981.

    [5] 路志英,趙智超,郝為,等.基于人工神經(jīng)網(wǎng)絡(luò)的多模型綜合預(yù)報(bào)方法[J].計(jì)算機(jī)應(yīng)用,2004,24(4):50-51,88.(LU Zhi-ying,ZHAO Zhi-chao,HAO Wei,et al.Multi-model ensemble forecast method based on ANN[J].Computer Applications,2004,24(4):50-51,88.(in Chinese))

    [6] 熊立華,郭生練,葉凌云.自適應(yīng)神經(jīng)模糊推理系統(tǒng)(ANFIS)在水文模型綜合中的應(yīng)用[J].水文,2006,26(1):38-41.(XIONG Li-hua,GUO Sheng-lian,YE Ling-yun.Application of ANFIS in the combination of hydrological models[J].Journal of China Hydrology,2006,26(1):38-41.(in Chinese)) DOI:10.3969/j.issn.1000-0852.2006.01.009.

    [7] Shamseldin A Y,O′Connor K M,Liang G C.Methods for combining the outputs of different rainfall-runoff models[J].Journal of Hydrology,1997(1):203-229.DOI:10.1016/S0022-1694(96)03259-3.

    [8] See L,Abrahart R J.Multi-model data fusion for hydrological forecasting[J].Computers & Geosciences,2001,27(8):987-994.DOI:10.1016/S0098-3004(00)00136-9.

    [9] Talei A,Chua L H C.Influence of lag time on event-based rainfall-runoff modeling using the data driven approach[J].Journal of Hydrology,438(1):223-233.DOI:10.1016/j.jhydrol.2012.03.027.

    [10] Astray G,Soto B,Lopez D,et al.Application of transit data analysis and artificial neural network in the prediction of discharge of Lor River,NW Spain[J].Water Science and Technology,73(7):1756-1767.DOI:10.2166/wst.2016.002.

    [11] Efron B.Bayesian,frequentists,and scientists[J].Journal of the American Statistical Association,2005,100(469):1-5.DOI:10.1198/016214505000000033.

    [12] Stock J H,Watson M W.Macroeconomic forecasting using diffusion indexes[J].Journal of Business & Economic Statistics,2002,20(2):147-162.DOI:10.1198/073500102317351921.

    [13] Zheng T,Salganik M J,Gelman A.How many people do you know in prison? Using Overdispersion in count data to estimate social structure in networks[J].Journal of the American Statistical Association,2006,101(474):409-423.DOI:10.1198/016214505000001168.

    [14] Duan Q Y,Phillips T J.Bayesian estimation of local signal and noise in multimodel simulations of climate change[J].Journal of Geophysical Research,2010,115(D18):D18123.DOI:10.1029/2009JD013654.

    [15] Ma F,Ye A Z,Deng X X,et al.Evaluating the skill of NMME seasonal precipitation ensemble predictions for 17 hydroclimatic regions in continental China[J].International Journal of Climatology,2016,36(1):132-144.DOI:10.1002/joc.4333.

    [16] Neuman S P.Maximum likelihood Bayesian averaging of uncertain model predictions[J].Stochastic Environmental Research and Risk Assessment,2003,17(5):291-305.DOI:10.1007/s00477-003-0151-7.

    [17] Kavetski D,Kuczera G,F(xiàn)ranks S W.Bayesian analysis of input uncertainty in hydrological modeling:2.Application[J].Water Resources Research,2006,42(3).DOI:10.1029/2005WR004376.

    [18] Zhang X S,Srinivasan R,Bosch D.Calibration and uncertainty analysis of the SWAT model using genetic algorithms and Bayesian model averaging[J].Journal of Hydrology,2009,374:307-317.DOI:10.1016/j.jhydrol.2009.06.023.

    [19] Raftery A E,Gneiting T,Balabdaooui F,et al.Using Bayesian model averaging to calibrate forecast ensembles[J].Monthly Weather Review,2005,133(5):1155-1174.DOI:10.1175/MWR2906.1.

    [20] Duan Q Y,Ajami N K,Gao X G,et al.Multi-model ensemble hydrologic prediction using Bayesian model averaging[J].Advances in Water Resources,2007,30(5):1317-1386.DOI:10.1016/j.advwatres.2006.11.014.

    [21] 梁忠民,戴榮,王軍,等.基于貝葉斯模型平均理論的水文模型合成預(yù)報(bào)研究[J].水力發(fā)電學(xué)報(bào),2010,29(2):114-118.(LIANG Zhong-min,DAI Rong,WANG Jun,et al.Study on forecast combination of different hydrological models by Bayesian model averaging[J].Journal of Hydroelectric Engineering,2010,29(2):114-118.(in Chinese))

    [22] Wu C H,Huang G R.Changes in heavy precipitation and floods in the upstream of the Beijiang River Basin,South China[J].International Journal of Climatology,2015,35(10):2978-2992.DOI:10.1002/joc.4187.

    [23] Wu C H,Huang G R,Yu H J,et al.Spatial and temporal distributions of trends in climate extremes of the Feilaixia catchment in the upstream area of the Beijiang River Basin,South China[J].International Journal of Climatology,2014,34(11):3161-3178.DOI:10.1002/joc.3900.

    [24] 葉愛中,段青云,徐靜.水文集合預(yù)報(bào)概述及模型案例[M].北京:中國水利水電出版社,2014.(YE Ai-zhong,DUAN Qing-yun,XU Jing.Overview of hydrological ensemble forecast and model cases[M].Beijing:China Water & Power Press,2014.(in Chinese))

    [25] Ye A Z,Duan Q Y,Zeng H J,et al.A distributed time-variant gain hydrological model based on remote sensing[J].Journal of Resources and Ecology,2010,1(3):222-229.DOI:10.3969/j.issn.1674-764x.2010.03.005.

    [26] 詹道江,葉守澤.工程水文學(xué)[M].北京:水利電力出版社,2004.(ZHAN Dao-jiang,YE Shou-ze.Engineering hydrology[M].Beijing:China Water & Power Press,2008.(in Chinese))

    [27] 李藝婷,張萬順,朱齊艷,等.新安江模型產(chǎn)流模式在山區(qū)型和平原型水庫入庫徑流量中的模擬與應(yīng)用[J].水資源與水工程學(xué)報(bào),2007,18(1):33-36.(LI Yi-ting,ZHANG Wan-shun,ZHU Qi-yan,et al.Application of runoff generation mode of Xinanjiang model in mountain and plain reservoir runoff simulation[J].Journal of Water Resources and Water Engineering,2007,18(1):33-36.(in Chinese)) DOI:10.3969/j.issn.1672-643X.2007.01.008.

    [28] Ju Q,Yu Z B,Hao Z C,et al.Division-based rainfall-runoff simulations with BP neural nnetworks and Xinanjiang model[J].Neurocomputing,2009,72(13):2873-2883.DOI:10.1016/j.neucom.2008.12.032.

    [29] 熊劍鋒,張樹全.薩克拉門托水文模型在天生橋洪水預(yù)報(bào)的應(yīng)用[J].云南水力發(fā)電,1999,15(1):4-7.(XIONG Jian-feng,ZHANG Shu-quan.Application of Sacramento hydrological model in flood forecast for Tianshengqiao-I hydropower project[J].Yunnan Water Power,1999,15(1):4-7.(in Chinese))

    [30] 張衛(wèi)華,李雨,魏朝富,等.不同水文模型在Broken流域的比較研究[J].西南師范大學(xué)學(xué)報(bào):自然科學(xué)版,2011,36(4):211-216.(ZHANG Wei-hua,LI Yu,WEI Chao-fu,et al.A comparative study of different hydrological model and their application in Broken Catchment[J].Journal of Southwest China Normal University (Natural Science Edition),2011,36(4):211-216.(in Chinese)) DOI:10.13718/j.cnki.xsxb.2011.04.007.

    [31] Carpenter T M,Georgakakos K P.Intercomparison of lumped versus distributed hydrologic model ensemble simulations on operational forecast scales[J].Journal of Hydrology,2006,329(1):174-185.DOI:10.1016/j.jhydrol.2006.02.013.

    [32] 董磊華,熊立華,萬民.基于貝葉斯模型加權(quán)平均方法的水文模型不確定性分析[J].水利學(xué)報(bào),2011,42(9):1065-1074.(DONG Lei-hua,XIONG Li-hua,WAN Min.Uncertainty analysis of hydrological modeling using the Bayesian Model Averaging Method[J].Journal of Hydraulic Engineering,2011,42(9):1065-1074.(in Chinese)) DOI:10.13243/j.cnki.slxb.2011.09.002.

    [33] Hoteling J A,Raftery M D,Volinsky A E.Bayesian Model Averaging:a tutorial[J].Statistical Science,1999,14(4):382-401.DOI:10.1214/ss/1009212814.

    [34] 夏軍,王綱勝,談戈,等.水文非線性系統(tǒng)與分布式時(shí)變增益模型[J].中國科學(xué),2004,34(11):1062-1071.(XIA Jun,WANG Gang-sheng,TAN Ge,et al.Nonlinear hydrological system and distributed time variant gain modeling[J].Science in China Series D:Earth Sciences,2004,34(11):1062-1071.(in Chinese)) DOI:10.3969/j.issn.1674-7240.2004.11.009.

    [35] Schaefli B,Gupta H V.Do Nash values have value?[J].Hydrological Processes,2007,21(15):2075-2080.DOI:10.1002/hyp.6825.

    [36] Pushpalatha R,Perrin C,Le Moine N,et al.A review of efficiency criteria suitable for evaluating low-flow simulations[J].Journal of Hydrology,2012,420:171-182.DOI:10.1016/j.jhydrol.2011.11.055.

    猜你喜歡
    不確定性
    法律的兩種不確定性
    法律方法(2022年2期)2022-10-20 06:41:56
    不確定性下的生態(tài)治理——以三江源草地修復(fù)為例
    科學(xué)(2020年6期)2020-02-06 08:59:50
    英鎊或繼續(xù)面臨不確定性風(fēng)險(xiǎn)
    中國外匯(2019年7期)2019-07-13 05:45:04
    英國“脫歐”不確定性增加 玩具店囤貨防漲價(jià)
    玩具世界(2019年6期)2019-05-21 03:30:02
    具有凸多面體不確定性的混雜隨機(jī)微分方程的鎮(zhèn)定分析
    考慮風(fēng)電功率與需求響應(yīng)不確定性的備用容量配置
    考慮系統(tǒng)不確定性的高超聲速飛行器容錯(cuò)控制
    具有不可測動態(tài)不確定性非線性系統(tǒng)的控制
    個(gè)人投資理財(cái)中不確定性問題的探討
    中國市場(2016年44期)2016-05-17 05:15:00
    不確定性與農(nóng)民專業(yè)合作社縱向一體化經(jīng)營
    久久久久国产一级毛片高清牌| 精品久久久久久电影网| 亚洲国产中文字幕在线视频| www日本在线高清视频| 国产高清视频在线播放一区| 黑人巨大精品欧美一区二区mp4| 黄色片一级片一级黄色片| 在线视频色国产色| 最好的美女福利视频网| 亚洲国产欧美一区二区综合| 久久精品91无色码中文字幕| 国产精品爽爽va在线观看网站 | 久久久久久久精品吃奶| 天堂动漫精品| 欧美乱码精品一区二区三区| 欧美最黄视频在线播放免费 | 天堂影院成人在线观看| 性欧美人与动物交配| 两人在一起打扑克的视频| 久久草成人影院| 女人被狂操c到高潮| 麻豆一二三区av精品| 在线观看免费视频日本深夜| 男女做爰动态图高潮gif福利片 | 欧美色视频一区免费| 熟女少妇亚洲综合色aaa.| 国产亚洲精品一区二区www| 女人被躁到高潮嗷嗷叫费观| 国产欧美日韩精品亚洲av| 亚洲精华国产精华精| 最新在线观看一区二区三区| 色播在线永久视频| 日韩欧美国产一区二区入口| 视频区欧美日本亚洲| 青草久久国产| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品在线美女| 99精品欧美一区二区三区四区| 免费搜索国产男女视频| 午夜久久久在线观看| 日韩国内少妇激情av| 日韩精品青青久久久久久| 亚洲熟妇中文字幕五十中出 | 999久久久国产精品视频| 99精品欧美一区二区三区四区| 亚洲五月天丁香| 中出人妻视频一区二区| 丝袜美足系列| 国产亚洲精品第一综合不卡| 男女之事视频高清在线观看| 日韩欧美在线二视频| 欧美性长视频在线观看| 午夜a级毛片| 亚洲欧美一区二区三区黑人| 亚洲avbb在线观看| 久热爱精品视频在线9| 国产一区在线观看成人免费| 一级片免费观看大全| 久久久国产一区二区| 国产精品久久久久成人av| 国产亚洲精品第一综合不卡| 国产成人精品无人区| 国产97色在线日韩免费| www日本在线高清视频| 免费少妇av软件| 亚洲狠狠婷婷综合久久图片| 国产色视频综合| 日本vs欧美在线观看视频| 欧美日韩乱码在线| 亚洲av成人一区二区三| 精品久久久久久久久久免费视频 | 99精品久久久久人妻精品| 久久精品成人免费网站| 亚洲精华国产精华精| 亚洲色图 男人天堂 中文字幕| 久久九九热精品免费| 极品人妻少妇av视频| 欧美中文综合在线视频| av国产精品久久久久影院| 一进一出抽搐动态| 欧美成狂野欧美在线观看| 亚洲av熟女| 高清黄色对白视频在线免费看| 91国产中文字幕| aaaaa片日本免费| 免费av中文字幕在线| av国产精品久久久久影院| 真人做人爱边吃奶动态| 亚洲精品国产色婷婷电影| 国产成人精品久久二区二区免费| 水蜜桃什么品种好| 精品人妻1区二区| 久久精品影院6| 99国产综合亚洲精品| 免费在线观看完整版高清| 成年人免费黄色播放视频| 男女之事视频高清在线观看| 亚洲国产欧美一区二区综合| 18禁裸乳无遮挡免费网站照片 | 久久精品亚洲熟妇少妇任你| 一区二区日韩欧美中文字幕| 黄片大片在线免费观看| 国产精品电影一区二区三区| 国产野战对白在线观看| 免费搜索国产男女视频| 欧美久久黑人一区二区| 欧美+亚洲+日韩+国产| 视频区欧美日本亚洲| x7x7x7水蜜桃| 国产亚洲精品综合一区在线观看 | 久久久久久人人人人人| 九色亚洲精品在线播放| 欧美日韩精品网址| 精品一区二区三区视频在线观看免费 | 国产精品野战在线观看 | 国产熟女午夜一区二区三区| 黄色a级毛片大全视频| 亚洲三区欧美一区| 国产精品九九99| 一进一出抽搐动态| 国产精品久久久久成人av| 搡老岳熟女国产| 夜夜爽天天搞| 精品久久久久久久久久免费视频 | 电影成人av| 高清在线国产一区| 激情视频va一区二区三区| 亚洲国产精品sss在线观看 | 日韩一卡2卡3卡4卡2021年| 午夜久久久在线观看| 老司机亚洲免费影院| 一级a爱片免费观看的视频| 国产免费现黄频在线看| 久久久久九九精品影院| 日本三级黄在线观看| 欧美人与性动交α欧美软件| 涩涩av久久男人的天堂| 国产精品爽爽va在线观看网站 | 亚洲 欧美一区二区三区| 美女福利国产在线| 中文字幕av电影在线播放| 亚洲一区中文字幕在线| 超碰97精品在线观看| 日本欧美视频一区| 精品一品国产午夜福利视频| 国产精品久久久人人做人人爽| 91麻豆精品激情在线观看国产 | 法律面前人人平等表现在哪些方面| 正在播放国产对白刺激| 黄片播放在线免费| 国产精品野战在线观看 | 午夜福利免费观看在线| 俄罗斯特黄特色一大片| 亚洲七黄色美女视频| 久久天躁狠狠躁夜夜2o2o| 久久精品亚洲av国产电影网| 国产国语露脸激情在线看| 欧美丝袜亚洲另类 | 涩涩av久久男人的天堂| 久久天躁狠狠躁夜夜2o2o| 国产aⅴ精品一区二区三区波| 不卡一级毛片| 黄频高清免费视频| 久久影院123| 国产精品综合久久久久久久免费 | 满18在线观看网站| 精品第一国产精品| 国产精品久久久av美女十八| 在线观看免费日韩欧美大片| 日韩免费高清中文字幕av| 老熟妇仑乱视频hdxx| 亚洲精品一二三| 大型av网站在线播放| 日本wwww免费看| 精品国产乱子伦一区二区三区| 悠悠久久av| 亚洲狠狠婷婷综合久久图片| 搡老乐熟女国产| 12—13女人毛片做爰片一| 午夜福利免费观看在线| 国产高清国产精品国产三级| 国产三级在线视频| 久久久国产成人免费| 亚洲av片天天在线观看| 黄色女人牲交| 亚洲人成电影免费在线| 欧美激情久久久久久爽电影 | 久久久久国内视频| 亚洲av片天天在线观看| 男女做爰动态图高潮gif福利片 | 国产xxxxx性猛交| 国产免费av片在线观看野外av| 最近最新中文字幕大全免费视频| 国产精品亚洲av一区麻豆| 50天的宝宝边吃奶边哭怎么回事| 美女国产高潮福利片在线看| 男人舔女人的私密视频| 在线免费观看的www视频| 亚洲aⅴ乱码一区二区在线播放 | 国产视频一区二区在线看| 99精品欧美一区二区三区四区| 亚洲欧美激情综合另类| 亚洲人成网站在线播放欧美日韩| 亚洲精品国产精品久久久不卡| 久久草成人影院| 亚洲熟女毛片儿| 久久久久久大精品| 在线播放国产精品三级| 国产男靠女视频免费网站| 男女下面插进去视频免费观看| 国产aⅴ精品一区二区三区波| 亚洲av成人不卡在线观看播放网| 国产成人影院久久av| 在线国产一区二区在线| 大型黄色视频在线免费观看| av欧美777| 精品第一国产精品| 国产成人精品久久二区二区免费| 精品少妇一区二区三区视频日本电影| 丁香欧美五月| 99国产综合亚洲精品| 极品人妻少妇av视频| 中文字幕av电影在线播放| 三上悠亚av全集在线观看| 免费女性裸体啪啪无遮挡网站| 免费久久久久久久精品成人欧美视频| av天堂在线播放| 久久精品成人免费网站| 不卡av一区二区三区| www.999成人在线观看| 亚洲一区高清亚洲精品| 亚洲一区二区三区色噜噜 | 免费在线观看亚洲国产| 欧美乱色亚洲激情| 宅男免费午夜| 久久久精品国产亚洲av高清涩受| 亚洲精品av麻豆狂野| 国产精品久久久av美女十八| 99久久国产精品久久久| 一级黄色大片毛片| 少妇裸体淫交视频免费看高清 | 男女做爰动态图高潮gif福利片 | 亚洲在线自拍视频| 午夜福利欧美成人| 视频区图区小说| 97人妻天天添夜夜摸| videosex国产| 黄色a级毛片大全视频| 波多野结衣av一区二区av| 91国产中文字幕| 超碰成人久久| videosex国产| 看免费av毛片| 久久中文看片网| 久久久久久亚洲精品国产蜜桃av| 纯流量卡能插随身wifi吗| 大型av网站在线播放| 一级毛片女人18水好多| 一进一出好大好爽视频| 成人亚洲精品av一区二区 | 欧美另类亚洲清纯唯美| 色在线成人网| 亚洲欧美激情综合另类| 欧美日韩亚洲高清精品| 丝袜人妻中文字幕| 欧美精品啪啪一区二区三区| 久久中文字幕人妻熟女| 麻豆一二三区av精品| 亚洲精品一二三| 国产单亲对白刺激| 99国产极品粉嫩在线观看| 国产色视频综合| 欧美激情久久久久久爽电影 | 国产成人一区二区三区免费视频网站| 身体一侧抽搐| 超色免费av| 国产伦一二天堂av在线观看| 19禁男女啪啪无遮挡网站| 亚洲avbb在线观看| 亚洲第一av免费看| 无人区码免费观看不卡| 亚洲成人久久性| 亚洲精华国产精华精| 美国免费a级毛片| 啦啦啦 在线观看视频| www.精华液| 最好的美女福利视频网| 午夜老司机福利片| 午夜视频精品福利| 香蕉国产在线看| 777久久人妻少妇嫩草av网站| av在线天堂中文字幕 | 国产蜜桃级精品一区二区三区| 桃色一区二区三区在线观看| 99久久久亚洲精品蜜臀av| 18美女黄网站色大片免费观看| 黑人操中国人逼视频| 欧美激情 高清一区二区三区| 在线av久久热| 97碰自拍视频| 国产成人av教育| 色哟哟哟哟哟哟| 免费搜索国产男女视频| a级毛片黄视频| 露出奶头的视频| xxx96com| 80岁老熟妇乱子伦牲交| 亚洲男人的天堂狠狠| 97碰自拍视频| 三级毛片av免费| 夜夜看夜夜爽夜夜摸 | 免费一级毛片在线播放高清视频 | 在线观看免费日韩欧美大片| 中文亚洲av片在线观看爽| 免费女性裸体啪啪无遮挡网站| 在线观看免费高清a一片| 岛国视频午夜一区免费看| 欧美日韩亚洲国产一区二区在线观看| 一级a爱片免费观看的视频| x7x7x7水蜜桃| 人人妻,人人澡人人爽秒播| 一级作爱视频免费观看| 免费在线观看影片大全网站| 后天国语完整版免费观看| 欧美精品亚洲一区二区| 亚洲aⅴ乱码一区二区在线播放 | 成人亚洲精品av一区二区 | 国产人伦9x9x在线观看| 亚洲少妇的诱惑av| 后天国语完整版免费观看| 久热爱精品视频在线9| 亚洲国产精品一区二区三区在线| 亚洲欧美一区二区三区久久| 黑丝袜美女国产一区| 精品国产乱子伦一区二区三区| 免费av毛片视频| 成人三级做爰电影| 亚洲伊人色综图| 久9热在线精品视频| 日本黄色视频三级网站网址| 亚洲五月婷婷丁香| 久久香蕉精品热| 交换朋友夫妻互换小说| 757午夜福利合集在线观看| www.999成人在线观看| 91大片在线观看| 两个人免费观看高清视频| 可以免费在线观看a视频的电影网站| 亚洲 国产 在线| tocl精华| 不卡一级毛片| 麻豆国产av国片精品| 丰满饥渴人妻一区二区三| 交换朋友夫妻互换小说| 热99国产精品久久久久久7| 国产av精品麻豆| 色播在线永久视频| 成在线人永久免费视频| 色播在线永久视频| videosex国产| netflix在线观看网站| 在线天堂中文资源库| 国产在线精品亚洲第一网站| 日本精品一区二区三区蜜桃| 亚洲精品久久午夜乱码| 国产欧美日韩精品亚洲av| 黑人欧美特级aaaaaa片| 天堂√8在线中文| 50天的宝宝边吃奶边哭怎么回事| 欧美午夜高清在线| 伦理电影免费视频| 国产一卡二卡三卡精品| 欧美丝袜亚洲另类 | 老鸭窝网址在线观看| 久久久久久久久中文| 亚洲精品粉嫩美女一区| 国产精品电影一区二区三区| 国产一区二区激情短视频| 可以在线观看毛片的网站| 欧美日本中文国产一区发布| 久久影院123| 国产91精品成人一区二区三区| 成人影院久久| 99热只有精品国产| 人妻丰满熟妇av一区二区三区| 88av欧美| 精品一区二区三区av网在线观看| 长腿黑丝高跟| 美女福利国产在线| 国产精品自产拍在线观看55亚洲| 国产免费男女视频| 中文字幕另类日韩欧美亚洲嫩草| 欧美乱色亚洲激情| 日韩av在线大香蕉| 午夜激情av网站| 日韩人妻精品一区2区三区| 亚洲一区中文字幕在线| av超薄肉色丝袜交足视频| 国产精品99久久99久久久不卡| 一个人观看的视频www高清免费观看 | 久久久久久久久免费视频了| 中文字幕人妻丝袜一区二区| 丁香六月欧美| 成人亚洲精品一区在线观看| 精品久久久精品久久久| 日韩欧美免费精品| 免费在线观看视频国产中文字幕亚洲| 亚洲精品中文字幕一二三四区| 在线十欧美十亚洲十日本专区| 亚洲av第一区精品v没综合| 午夜福利影视在线免费观看| 最近最新中文字幕大全免费视频| 国产精品一区二区精品视频观看| av网站免费在线观看视频| 国产精品国产高清国产av| 亚洲黑人精品在线| 女同久久另类99精品国产91| 大香蕉久久成人网| 欧美亚洲日本最大视频资源| 国产国语露脸激情在线看| 久久欧美精品欧美久久欧美| e午夜精品久久久久久久| 亚洲国产精品999在线| 久久草成人影院| 美女午夜性视频免费| 婷婷丁香在线五月| 婷婷六月久久综合丁香| 国产精品一区二区在线不卡| 男女床上黄色一级片免费看| 看黄色毛片网站| 亚洲av电影在线进入| 一本综合久久免费| 亚洲七黄色美女视频| 老熟妇乱子伦视频在线观看| 久久精品成人免费网站| 精品一区二区三区视频在线观看免费 | 制服人妻中文乱码| 欧美黄色片欧美黄色片| 怎么达到女性高潮| 黄色视频,在线免费观看| 国产免费男女视频| av网站免费在线观看视频| 麻豆成人av在线观看| 91国产中文字幕| 纯流量卡能插随身wifi吗| 精品久久久久久电影网| 久久国产精品男人的天堂亚洲| 精品久久久久久成人av| 久久久久久久精品吃奶| 精品福利观看| 亚洲人成电影观看| 麻豆久久精品国产亚洲av | 美女大奶头视频| 久久久久国产一级毛片高清牌| 80岁老熟妇乱子伦牲交| 国产99久久九九免费精品| 国产一区二区在线av高清观看| 久久热在线av| 亚洲第一青青草原| 可以在线观看毛片的网站| 90打野战视频偷拍视频| 精品高清国产在线一区| 亚洲av日韩精品久久久久久密| av电影中文网址| 高清黄色对白视频在线免费看| 午夜两性在线视频| 久久久精品国产亚洲av高清涩受| 国产又爽黄色视频| 搡老乐熟女国产| 美女高潮到喷水免费观看| av电影中文网址| 久久天躁狠狠躁夜夜2o2o| 亚洲国产欧美日韩在线播放| 老司机午夜十八禁免费视频| 国产在线观看jvid| 黑人欧美特级aaaaaa片| 99精国产麻豆久久婷婷| 久热爱精品视频在线9| 欧美人与性动交α欧美精品济南到| 亚洲精品一区av在线观看| 法律面前人人平等表现在哪些方面| 99精品欧美一区二区三区四区| 丝袜美足系列| 电影成人av| 每晚都被弄得嗷嗷叫到高潮| 99久久精品国产亚洲精品| 国产不卡一卡二| 我的亚洲天堂| 黄频高清免费视频| 亚洲 欧美 日韩 在线 免费| 亚洲人成伊人成综合网2020| 亚洲精品久久成人aⅴ小说| 国产精品一区二区免费欧美| 久久久精品欧美日韩精品| 女性被躁到高潮视频| 色精品久久人妻99蜜桃| 亚洲精品中文字幕在线视频| 十分钟在线观看高清视频www| 日韩精品免费视频一区二区三区| 怎么达到女性高潮| 亚洲av第一区精品v没综合| 成人特级黄色片久久久久久久| 极品人妻少妇av视频| 久久精品aⅴ一区二区三区四区| 精品国内亚洲2022精品成人| 少妇裸体淫交视频免费看高清 | 久久精品影院6| 国产片内射在线| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲va日本ⅴa欧美va伊人久久| 久久久久久久午夜电影 | 亚洲熟妇熟女久久| 欧美 亚洲 国产 日韩一| 男女之事视频高清在线观看| 精品国产亚洲在线| 十八禁人妻一区二区| 久久精品亚洲精品国产色婷小说| 十八禁人妻一区二区| 国产99白浆流出| 97超级碰碰碰精品色视频在线观看| 日韩欧美国产一区二区入口| 日本三级黄在线观看| 亚洲第一青青草原| 国产伦一二天堂av在线观看| cao死你这个sao货| 午夜精品国产一区二区电影| 操美女的视频在线观看| 他把我摸到了高潮在线观看| 高潮久久久久久久久久久不卡| 午夜亚洲福利在线播放| www日本在线高清视频| 国产欧美日韩一区二区精品| 久久欧美精品欧美久久欧美| 伦理电影免费视频| 99精品欧美一区二区三区四区| 欧美乱色亚洲激情| 国产精品成人在线| ponron亚洲| 亚洲一卡2卡3卡4卡5卡精品中文| 黑人巨大精品欧美一区二区蜜桃| 男女高潮啪啪啪动态图| 亚洲av美国av| 黑人操中国人逼视频| 69av精品久久久久久| ponron亚洲| 两人在一起打扑克的视频| 亚洲第一青青草原| 男男h啪啪无遮挡| 麻豆av在线久日| 欧美亚洲日本最大视频资源| 精品日产1卡2卡| 亚洲成人国产一区在线观看| 亚洲午夜理论影院| 日韩欧美一区二区三区在线观看| 高清av免费在线| 天堂俺去俺来也www色官网| 亚洲精品美女久久av网站| 精品国产乱子伦一区二区三区| 婷婷六月久久综合丁香| 99久久综合精品五月天人人| 亚洲国产毛片av蜜桃av| 亚洲成a人片在线一区二区| 精品国产超薄肉色丝袜足j| 电影成人av| 老司机在亚洲福利影院| 99riav亚洲国产免费| 欧美精品亚洲一区二区| 18美女黄网站色大片免费观看| 老司机深夜福利视频在线观看| 久久中文字幕一级| www日本在线高清视频| 国产有黄有色有爽视频| 香蕉久久夜色| 99国产精品一区二区蜜桃av| 无限看片的www在线观看| 精品久久久久久久毛片微露脸| 久久精品亚洲精品国产色婷小说| 俄罗斯特黄特色一大片| 欧美国产精品va在线观看不卡| 叶爱在线成人免费视频播放| 热re99久久国产66热| 国内毛片毛片毛片毛片毛片| 亚洲av成人一区二区三| 性色av乱码一区二区三区2| 亚洲欧美一区二区三区久久| 日韩免费高清中文字幕av| www.999成人在线观看| 久久久久久久久免费视频了| 亚洲精品美女久久av网站| 精品国产乱子伦一区二区三区| 国产欧美日韩精品亚洲av| av国产精品久久久久影院| 精品乱码久久久久久99久播| 亚洲国产中文字幕在线视频| 午夜免费观看网址| 级片在线观看| 一级a爱片免费观看的视频| 日本vs欧美在线观看视频| 99久久人妻综合| 国产单亲对白刺激| 亚洲国产精品999在线| 久久这里只有精品19| 999久久久国产精品视频| 人人妻人人澡人人看| 久久国产乱子伦精品免费另类| 十分钟在线观看高清视频www| 日韩成人在线观看一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 一a级毛片在线观看| 淫妇啪啪啪对白视频| 俄罗斯特黄特色一大片| 九色亚洲精品在线播放| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日本亚洲视频在线播放| 亚洲avbb在线观看|