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

    基于遺傳算法的條件植被溫度指數(shù)的時(shí)間尺度轉(zhuǎn)換方法

    2018-05-05 08:44:49白雪嬌王鵬新張樹譽(yù)景毅剛劉峻明
    水土保持研究 2018年1期
    關(guān)鍵詞:關(guān)中平原時(shí)間尺度單產(chǎn)

    白雪嬌, 王鵬新, 張樹譽(yù), 李 俐, 景毅剛, 劉峻明

    (1.中國(guó)農(nóng)業(yè)大學(xué) 信息與電氣工程學(xué)院, 北京 100083; 2.陜西省氣象局, 西安 710014)

    干旱持續(xù)時(shí)間長(zhǎng)、強(qiáng)度大、影響范圍廣,是影響農(nóng)業(yè)生產(chǎn)能力的主要?dú)夂蛞蛩?,加?qiáng)干旱的監(jiān)測(cè)預(yù)警和影響評(píng)估,對(duì)農(nóng)業(yè)生產(chǎn)及其發(fā)展具有重要意義。隨著具有空間連續(xù)性的多源遙感技術(shù)的發(fā)展,為實(shí)時(shí)、精確、動(dòng)態(tài)和大范圍反演地面信息提供了技術(shù)手段,在旱情監(jiān)測(cè)預(yù)警和影響評(píng)估等方面得到廣泛的應(yīng)用[1-2]。關(guān)中平原是陜西省自然條件最好的地區(qū)和農(nóng)業(yè)、果業(yè)生產(chǎn)基地,但其地處內(nèi)陸、遠(yuǎn)離海洋,對(duì)干旱的響應(yīng)極為敏感,農(nóng)業(yè)干旱的季節(jié)性和區(qū)域性特征明顯,成為區(qū)域農(nóng)業(yè)生產(chǎn)的主要障礙因子[3-4]。國(guó)內(nèi)外研究者利用改進(jìn)型的垂直干旱指數(shù)(MPDI)、植被水分脅迫指數(shù)(VWSI)等指標(biāo)和CERES-Wheat模型等進(jìn)行了較多旱情的監(jiān)測(cè)預(yù)警和影響評(píng)估研究[5-7],發(fā)展了豐富的旱情遙感監(jiān)測(cè)預(yù)警和影響評(píng)估的理論和方法。條件植被溫度指數(shù)(VTCI)綜合了反射和輻射信息,可準(zhǔn)確地進(jìn)行區(qū)域級(jí)干旱監(jiān)測(cè)研究,已被國(guó)內(nèi)外學(xué)者廣泛應(yīng)用于旱情監(jiān)測(cè)、影響評(píng)估和作物產(chǎn)量估測(cè)[8-11]。

    冬小麥越冬后的生育時(shí)期包括返青、拔節(jié)、抽穗、灌漿和乳熟等時(shí)期。干旱發(fā)生時(shí)段和程度不同造成的減產(chǎn)率也不同,拔節(jié)期和抽穗期受旱會(huì)影響穗粒數(shù)及成熟期提前,灌漿期和成熟期受旱會(huì)影響千粒重和籽粒品質(zhì)[12]。研究冬小麥生長(zhǎng)規(guī)律及不同生育時(shí)期干旱對(duì)產(chǎn)量的影響程度能在一定程度上提供干旱對(duì)產(chǎn)量影響的有用信息,如能綜合利用不同生育時(shí)期多時(shí)相VTCI對(duì)產(chǎn)量的影響信息,對(duì)其賦予一定的權(quán)重進(jìn)行時(shí)間尺度轉(zhuǎn)換為一個(gè)綜合了多時(shí)相干旱信息的加權(quán)VTCI值,勢(shì)必可以提高冬小麥估產(chǎn)和干旱影響評(píng)估的精度。關(guān)鍵是如何確定不同生育時(shí)期旱情對(duì)冬小麥產(chǎn)量影響的權(quán)重系數(shù)。多指標(biāo)權(quán)重的確定方法大致分為主觀賦權(quán)法和客觀賦權(quán)法兩大類。主觀賦權(quán)法如專家調(diào)查法、層次分析法等,是基于決策者經(jīng)驗(yàn)判斷或主觀認(rèn)識(shí)決定指標(biāo)權(quán)重,主觀因素介入過(guò)多而使得評(píng)價(jià)結(jié)果主觀性和隨意性較大;客觀賦權(quán)法如變異系數(shù)法、主成分分析法,是基于指標(biāo)間關(guān)系和其包含的信息量通過(guò)建立一定的數(shù)學(xué)模型來(lái)確定指標(biāo)權(quán)重,評(píng)價(jià)結(jié)果數(shù)學(xué)基礎(chǔ)較強(qiáng)而忽略了實(shí)際情況[13-14]。為兼顧兩類方法的優(yōu)點(diǎn)同時(shí)規(guī)避其局限性,綜合集成主、客觀信息的組合賦權(quán)法和其他權(quán)重集成優(yōu)化模型得到了廣泛應(yīng)用。遲道才等[15]運(yùn)用博弈論的綜合賦權(quán)方法將主客觀賦權(quán)法相融合對(duì)東港灌區(qū)的運(yùn)行狀況進(jìn)行了評(píng)價(jià)研究。王書吉等[16]綜合主客觀賦權(quán)法的優(yōu)點(diǎn)提出基于單位化約束條件的綜合集成賦權(quán)法,確定了關(guān)中大型灌區(qū)的節(jié)水改造綜合評(píng)價(jià)指標(biāo)。曹永強(qiáng)等[17]采用一種基于離差最大化原理的組合賦權(quán)方法確定了農(nóng)業(yè)干旱風(fēng)險(xiǎn)綜合評(píng)價(jià)指標(biāo)的權(quán)重。

    本文采用窮舉法(EA)和遺傳算法(GA)求算各生育時(shí)期權(quán)重系數(shù)對(duì)冬小麥主要生育期內(nèi)多時(shí)相VTCI進(jìn)行時(shí)間尺度轉(zhuǎn)換以評(píng)估干旱對(duì)產(chǎn)量的影響。窮舉法的基本思想是依題目的條件確定問(wèn)題解的大致搜索范圍,對(duì)所有可能的情況逐個(gè)驗(yàn)證[18]。窮舉法應(yīng)用很多,如公元五世紀(jì)我國(guó)數(shù)學(xué)家張丘建在其《算經(jīng)》一書中提出的“百雞問(wèn)題”,一些密碼破譯、刪余卷積碼好碼的搜索等。窮舉法可得到問(wèn)題的最優(yōu)解且算法簡(jiǎn)單,但運(yùn)算量大導(dǎo)致計(jì)算效率低下。遺傳算法是一種基于自然選擇的模擬優(yōu)化算法,它模擬自然進(jìn)化的高效率搜索解決復(fù)雜的優(yōu)化問(wèn)題,因其通用、實(shí)用和高效的特點(diǎn)而被應(yīng)用于各個(gè)領(lǐng)域,尤其對(duì)解決一些機(jī)器學(xué)習(xí)、信號(hào)處理和組合優(yōu)化問(wèn)題如旅行商問(wèn)題[19]、作業(yè)車間調(diào)度問(wèn)題[20]等有獨(dú)特的優(yōu)越性。本文利用遺傳算法解決最優(yōu)權(quán)重搜索,可實(shí)現(xiàn)快速收斂至全局最優(yōu)解之目的。當(dāng)前基于窮舉法和遺傳算法原理的算法在電子、電力和人工生命等領(lǐng)域應(yīng)用較多,但應(yīng)用于多時(shí)相遙感數(shù)據(jù)的時(shí)間尺度轉(zhuǎn)換研究方面還不多見(jiàn)。

    以關(guān)中平原冬小麥為研究對(duì)象,基于關(guān)中平原2008—2014年Aqua MODIS VTCI時(shí)間序列數(shù)據(jù),以不同生育時(shí)期干旱對(duì)冬小麥單產(chǎn)的影響為基準(zhǔn),首先采用窮舉法、遺傳算法和課題組前期研究成果基于因子權(quán)重排序法與熵值法的歸一組合賦權(quán)法(CAFE)[10]求解不同生育時(shí)期干旱對(duì)產(chǎn)量的影響的權(quán)重值,然后根據(jù)3種方法獲得的權(quán)重值對(duì)冬小麥主要生育期多時(shí)相VTCI進(jìn)行時(shí)間尺度轉(zhuǎn)換獲取綜合加權(quán)VTCI值,對(duì)其進(jìn)行精度評(píng)價(jià)。以窮舉法所得結(jié)果為參考對(duì)比分析遺傳算法和基于因子權(quán)重排序法與熵值法的歸一組合賦權(quán)法,為高效求解不同生育時(shí)期干旱對(duì)產(chǎn)量的影響的最佳權(quán)重提供可行的方法,提高基于多時(shí)相VTCI遙感數(shù)據(jù)時(shí)間尺度轉(zhuǎn)換的干旱影響評(píng)估研究的精度。

    1 數(shù)據(jù)來(lái)源與研究方法

    1.1 研究區(qū)概況

    關(guān)中平原位于陜西省中部,海拔高度325~800 m,總面積約5.55萬(wàn)km2。關(guān)中平原質(zhì)疏松肥沃,又有涇、渭和洛等河流提供水源,十分適宜耕作,是陜西省最重要的農(nóng)耕區(qū)。關(guān)中平原年平均溫度6~13.6℃,年均降水量500~600 mm,屬暖溫帶半干旱氣候,四季分明,雨熱同期,冬春降水較少,是氣候變化的敏感區(qū),春、伏旱頻繁發(fā)生,受旱成災(zāi)面積總體在波動(dòng)中上升,造成糧食減產(chǎn)。將各市行政邊界的矢量圖與研究區(qū)域的遙感影像疊加,對(duì)研究區(qū)域按照各市的行政邊界進(jìn)行市屬的劃分,得到關(guān)中平原分市概況圖(圖1)。由于楊凌示范區(qū)面積較小,故將其劃分到咸陽(yáng)市境內(nèi)。

    圖1研究區(qū)概況圖

    1.2 數(shù)據(jù)來(lái)源及預(yù)處理

    1.2.1 數(shù)據(jù)來(lái)源 本文數(shù)據(jù)源主要是2008—2014年3—5月Aqua MODIS數(shù)據(jù)產(chǎn)品,包括日反射率產(chǎn)品(MYD09 GA)和日LST產(chǎn)品(MYD11 A1)。冬小麥單產(chǎn)數(shù)據(jù)來(lái)源于歷年的《陜西省統(tǒng)計(jì)年鑒》。

    1.2.2 條件植被溫度指數(shù)的生成 條件植被溫度指數(shù)(VTCI)是基于遙感反演的歸一化植被指數(shù)(NDVI)和地表溫度(LST)的散點(diǎn)圖呈三角形區(qū)域分布的基礎(chǔ)上提出的,主要用于監(jiān)測(cè)旱情。計(jì)算方法為[2,9]:

    (1)

    LSTmax(NDVIi)=a+bNDVIi

    (2)

    LSTmin(NDVIi)=a′+b′NDVIi

    (3)

    式中:NDVIi為第i個(gè)時(shí)期(旬)某一像素的NDVI值;LST(NDVIi)為某一像素的NDVI值為NDVIi時(shí)的LST值;LSTmax(NDVIi)為研究區(qū)域內(nèi)當(dāng)NDVIi值等于某一特定值時(shí)LST的最大值;LSTmin(NDVIi)為研究區(qū)域內(nèi)當(dāng)NDVIi值等于某一特定值時(shí)LST的最小值;a,b,a′和b′為待定系數(shù),由研究區(qū)域內(nèi)的NDVI和LST的散點(diǎn)圖近似獲得。

    應(yīng)用多年的旬NDVI和LST最大值合成產(chǎn)品確定每旬多年共同熱邊界LSTmax,應(yīng)用多年的旬NDVI最大值合成產(chǎn)品和LST最大、最小值合成產(chǎn)品確定每旬多年共同的冷邊界LSTmin。VTCI取值范圍為0,1,VTCI值越小,旱情越嚴(yán)重,VTCI值越大,則相反。

    計(jì)算關(guān)中平原從2008—2014年3—5月共63旬的VTCI數(shù)據(jù)。根據(jù)越冬后冬小麥的生長(zhǎng)情況,將越冬后的生育時(shí)期劃分為返青期(3月上旬—3月中旬)、拔節(jié)期(3月下旬—4月中旬)、抽穗—灌漿期(4月下旬—5月上旬)和乳熟期(5月中旬—5月下旬)4個(gè)主要生育時(shí)期,取某一生育時(shí)期內(nèi)多旬VTCI的平均值作為該生育時(shí)期的VTCI,再依據(jù)關(guān)中平原分市VTCI概況圖統(tǒng)計(jì)出關(guān)中平原各市主要生育時(shí)期的VTCI。銅川市由于占所選研究區(qū)域內(nèi)面積較小,統(tǒng)計(jì)結(jié)果代表性差,因此文中不做研究。

    1.3 研究方法

    不同生育時(shí)期發(fā)生的干旱對(duì)冬小麥生長(zhǎng)發(fā)育的影響是不同的,對(duì)最終產(chǎn)量的影響程度也是不同的。所以在實(shí)際應(yīng)用中需考慮不同生育時(shí)期的干旱對(duì)冬小麥生長(zhǎng)的影響程度,將不同生育時(shí)期多時(shí)相VTCI遙感數(shù)據(jù)進(jìn)行時(shí)間尺度轉(zhuǎn)換生成綜合VTCI值以合理評(píng)估干旱對(duì)冬小麥生產(chǎn)的影響。

    時(shí)間尺度轉(zhuǎn)換權(quán)重優(yōu)化的目標(biāo)是使基于各生育時(shí)期干旱對(duì)產(chǎn)量的影響權(quán)重獲得的加權(quán)VTCI與冬小麥單產(chǎn)的相關(guān)性最大,數(shù)學(xué)模型為:

    VTCIwi=wT×VTCIT+wE×VTCIE+wH×VTCIH+wM×VTCIM

    (4)

    (5)

    wT+wE+wH+wM=1

    (6)

    約束條件:干旱的發(fā)生時(shí)段和程度不同,對(duì)作物生長(zhǎng)發(fā)育和產(chǎn)量的影響也不同。返青期主要是生根、長(zhǎng)葉和分蘗,需水量較小,干旱對(duì)小麥產(chǎn)量的影響最小,甚至在某種情況下會(huì)促進(jìn)冬小麥生長(zhǎng)速率;拔節(jié)期為水分虧缺敏感期,這一時(shí)期干旱對(duì)小麥產(chǎn)量的影響最大;抽穗—灌漿期小麥以生殖生長(zhǎng)為主,也是需水的關(guān)鍵期;乳熟期由于穗粒結(jié)構(gòu)已經(jīng)形成,對(duì)水分虧缺表現(xiàn)出一定的耐旱性,干旱對(duì)產(chǎn)量影響相對(duì)較小。各生育時(shí)期干旱對(duì)產(chǎn)量影響的權(quán)重值范圍如下:

    -0.5

    0

    0

    -0.5

    (7)

    1.3.1 窮舉法 窮舉法(Exhaustive Attack Method, EA)是基于計(jì)算機(jī)特點(diǎn)而進(jìn)行解題的思維方法。一般是根據(jù)問(wèn)題中的約束條件,在研究對(duì)象是由有限個(gè)元素構(gòu)成的集合時(shí),將所有可能的情況一一列舉出來(lái),然后逐一加以分析、處理,并驗(yàn)證結(jié)果是否滿足給定的條件,列舉完所有對(duì)象,問(wèn)題將最終得以解決[18]。算法的工作流程:(1) 確定搜索空間,即列舉出所有符合情況的wT,wE,wHF和wM;(2) 確定符合窮舉規(guī)則的條件;(3) 確定能使搜索過(guò)程優(yōu)化的條件,以縮小搜索范圍,減少搜索時(shí)間;(4) 遍歷所有解,獲取符合條件的最優(yōu)解。

    1.3.2 遺傳算法 遺傳算法(Genetic Algorithm,GA)通過(guò)選擇、交叉和變異三個(gè)基本算子模擬自然界的選擇與遺傳機(jī)理尋求最優(yōu)解。以GA為代表的智能算法較傳統(tǒng)優(yōu)化算法如窮舉法和啟發(fā)式算法等,具有很好的收斂性,計(jì)算時(shí)間少且魯棒性高,在優(yōu)化計(jì)算中得到廣泛的應(yīng)用。算法的工作流程:(1) 隨機(jī)產(chǎn)生一組初始群體,其中的每個(gè)個(gè)體(染色體)代表問(wèn)題的一個(gè)潛在解,即一組wT,wE,wHF和wM,并設(shè)置迭代次數(shù)kmax=200;(2) 確定控制參數(shù)px和pm,計(jì)算群體中每個(gè)個(gè)體相應(yīng)的適應(yīng)度值以找出最優(yōu)個(gè)體;(3) 隨機(jī)從種群中選擇個(gè)體按一定概率交叉和變異;(4) 計(jì)算目標(biāo)函數(shù),判斷進(jìn)化是否結(jié)束,否則轉(zhuǎn)步驟(2)。

    1.4 精度評(píng)價(jià)方法

    本文應(yīng)用Pearson相關(guān)系數(shù)(Correlation coefficient)、估計(jì)值的標(biāo)準(zhǔn)誤差(Std.Error of the estimate)和F檢驗(yàn)分析比較回歸模型精度,并應(yīng)用相對(duì)誤差(relative error,RE)與均方根誤差(root mean square error,RMSE)比較冬小麥的估測(cè)單產(chǎn)與年鑒統(tǒng)計(jì)單產(chǎn)之間的差異,以此評(píng)價(jià)不同賦權(quán)方法確定的權(quán)重應(yīng)用于關(guān)中平原冬小麥主要生育期多時(shí)相VTCI時(shí)間尺度轉(zhuǎn)換的精度。

    (8)

    (9)

    式中:Gi是冬小麥估測(cè)單產(chǎn);n是樣本個(gè)數(shù)。

    2 結(jié)果與分析

    2.1 不同賦權(quán)方法獲取的權(quán)重結(jié)果對(duì)比

    分別采用CAFE,EA和GA獲取關(guān)中平原各市和4市整體的冬小麥返青期、拔節(jié)期、抽穗—灌漿期和乳熟期干旱對(duì)產(chǎn)量的影響權(quán)重值,見(jiàn)表1。3種方法獲得的拔節(jié)期干旱對(duì)產(chǎn)量的影響的權(quán)重值最大,其次是抽穗—灌漿期的權(quán)重值,乳熟期和返青期的權(quán)重值較小(除4市整體),這與冬小麥實(shí)際生長(zhǎng)規(guī)律和他人研究結(jié)果相一致[21]。在拔節(jié)期冬小麥的營(yíng)養(yǎng)器官生長(zhǎng)旺盛,對(duì)水分吸收利用越加明顯,是生長(zhǎng)的關(guān)鍵期,在該生育時(shí)期發(fā)生干旱會(huì)造成根、莖、葉生長(zhǎng)速率降低和影響冬小麥的有效小穗數(shù)。抽穗灌漿期,冬小麥由營(yíng)養(yǎng)生長(zhǎng)轉(zhuǎn)向生殖生長(zhǎng),在抽穗灌漿期發(fā)生干旱會(huì)使得穗粒數(shù)減少,籽粒癟小,粒重明顯降低,是影響小麥生長(zhǎng)發(fā)育的又一關(guān)鍵時(shí)期。在返青期冬小麥生長(zhǎng)緩慢,對(duì)水分的需求量較小,由于作物的補(bǔ)償生長(zhǎng)效應(yīng),在返青期發(fā)生干旱對(duì)產(chǎn)量影響不大,在某種情況下甚至可以促進(jìn)冬小麥生長(zhǎng)發(fā)育和提高產(chǎn)量。冬小麥在乳熟期,營(yíng)養(yǎng)物質(zhì)向籽粒運(yùn)輸?shù)倪^(guò)程結(jié)束,植株逐漸枯萎,根系開(kāi)始死亡,因此在乳熟期發(fā)生干旱對(duì)產(chǎn)量影響較小。各市同一生育時(shí)期干旱對(duì)產(chǎn)量影響的權(quán)重值不同,這是由于不同地理位置、地形地貌影響下的旱情的主要影響因子氣溫和降水的時(shí)間分布和空間分布不同,同時(shí)各市灌溉地分布和面積不同。寶雞和4市整體采用EA和GA獲得的返青期干旱對(duì)產(chǎn)量影響的權(quán)重值分別為-0.01和-0.05,但其絕對(duì)值較小,這是由于干旱后的復(fù)水補(bǔ)償效應(yīng),關(guān)中平原易發(fā)生春旱,嚴(yán)重影響小穗數(shù)和穗粒數(shù),而根據(jù)歷年實(shí)地調(diào)研發(fā)現(xiàn),關(guān)中平原大部分地區(qū)會(huì)在此期間進(jìn)行春灌,正好解除了返青期水分脅迫和緩解了土壤水分的不足,有利于冬小麥生長(zhǎng)發(fā)育。咸陽(yáng)采用EA和GA獲得的乳熟期干旱對(duì)產(chǎn)量影響的權(quán)重值均為0,西安采用EA和GA獲得的乳熟期干旱對(duì)產(chǎn)量影響的權(quán)重值分別為0,0.01,均小于其返青期的權(quán)重值,這是由于咸陽(yáng)和西安兩市地處關(guān)中平原中部,冬小麥生長(zhǎng)后期降水較多,同時(shí),地下水儲(chǔ)量豐富、灌溉渠系發(fā)達(dá),灌溉面積較大缺灌溉條件較好,冬小麥生長(zhǎng)后期水分供應(yīng)充分,使得其乳熟期干旱對(duì)產(chǎn)量影響的權(quán)重值小于返青期。渭南由于其境內(nèi)旱作地分布較多而使得拔節(jié)期干旱對(duì)產(chǎn)量影響的權(quán)重值遠(yuǎn)大于其他生育時(shí)期,采用EA和GA獲得權(quán)重值相等,為0.83。

    根據(jù)3種方法的計(jì)算原理,關(guān)中平原冬小麥不同生育時(shí)期干旱對(duì)產(chǎn)量影響的最優(yōu)權(quán)重值可由EA獲得。CAFE確定的干旱對(duì)產(chǎn)量影響的權(quán)重值與EA確定的權(quán)重值相差較大,同時(shí)其經(jīng)時(shí)間尺度轉(zhuǎn)換后獲得的加權(quán)VTCI與冬小麥單產(chǎn)的相關(guān)性整體低于EA和GA,這可能是受CAFE中因子權(quán)重排序法數(shù)據(jù)處理過(guò)程中權(quán)重范圍設(shè)置區(qū)間和熵值法忽略了不同生育時(shí)期發(fā)生干旱對(duì)冬小麥產(chǎn)量的不同影響程度的影響。而GA獲得的不同生育時(shí)期干旱對(duì)產(chǎn)量的影響的權(quán)重值與最優(yōu)權(quán)重值相等或相近,相關(guān)性亦均較高,可見(jiàn),賦權(quán)方法GA獲取的不同生育時(shí)期干旱對(duì)產(chǎn)量影響的權(quán)重值的精度和性能幾乎與窮舉法相同,且優(yōu)于CAFE。從數(shù)據(jù)處理過(guò)程及其耗時(shí)上看,GA搜索到的最優(yōu)解比EA和CAFE的方案簡(jiǎn)約易行,以咸陽(yáng)市為例,EA用了18h13m,而GA只用41 s,時(shí)間比為1 600倍。綜合實(shí)際算例的統(tǒng)計(jì)來(lái)看,GA獲得的干旱對(duì)產(chǎn)量影響的權(quán)重值與最優(yōu)權(quán)重值的最大相差不超過(guò)0.01,與EA相比,GA的時(shí)間代價(jià)比窮舉法大為降低,計(jì)算效率是EA的幾百倍不等。從實(shí)際問(wèn)題角度,GA求解冬小麥生育時(shí)期優(yōu)化賦權(quán)問(wèn)題遠(yuǎn)優(yōu)于EA和CAFE。

    表1 不同賦權(quán)方法獲得的權(quán)重值及其精度分析

    注:*表示在0.1水平顯著,**表示在0.05水平顯著,***表示在0.01水平顯著。渭南、咸陽(yáng)、寶雞和西安的樣本量n=6,4市綜合的樣本量n=24。F4,0.1=0.729,F(xiàn)4,0.05=0.811,F(xiàn)4,0.01=0.917,F(xiàn)22,0.1=0.344,F(xiàn)22,0.05=0.404,F(xiàn)22,0.01=0.515。

    除渭南市的權(quán)重值外,由3種方法確定的干旱對(duì)冬小麥各生育時(shí)期產(chǎn)量影響的權(quán)重值時(shí)間尺度轉(zhuǎn)換后獲得的VTCI與冬小麥單產(chǎn)的線性相關(guān)性均達(dá)到顯著相關(guān),均通過(guò)了F檢驗(yàn),其中EA和GA確定地權(quán)重值獲得的尺度轉(zhuǎn)換后的VTCI與冬小麥單產(chǎn)的相關(guān)系數(shù)均大于CAFE,估計(jì)值標(biāo)準(zhǔn)誤差均小于CAFE,因此EA和GA能更好的反映關(guān)中平原冬小麥各生育時(shí)期VTCI與產(chǎn)量的線性相關(guān)性,在分析冬小麥主要生育時(shí)期干旱對(duì)產(chǎn)量的影響中具有較大的適用性,而CAFE的計(jì)算結(jié)果略差。GA與EA僅估計(jì)值標(biāo)準(zhǔn)誤差和顯著系數(shù)值存在微小差異,以綜合實(shí)際算例的統(tǒng)計(jì)來(lái)看,GA與EA的估計(jì)值標(biāo)準(zhǔn)誤差最大相差為0.12 kg/hm2,顯著系數(shù)值均相等,結(jié)合上文時(shí)間復(fù)雜度和計(jì)算效率分析,說(shuō)明GA對(duì)關(guān)中平原各生育時(shí)期權(quán)重值的確定較為合理,能夠較準(zhǔn)確地反映關(guān)中平原冬小麥主要生育時(shí)期的旱情對(duì)產(chǎn)量的影響程度,可用于開(kāi)展基于多時(shí)相VTCI數(shù)據(jù)的時(shí)間尺度轉(zhuǎn)換的關(guān)中平原干旱影響評(píng)估研究。

    2.2 不同賦權(quán)方法獲取的冬小麥估測(cè)單產(chǎn)和精度評(píng)價(jià)

    采用CAFE,EA和GA確定的干旱對(duì)產(chǎn)量影響的權(quán)重值經(jīng)時(shí)間尺度轉(zhuǎn)換后獲得加權(quán)VTCI值,計(jì)算2008—2014年關(guān)中平原各市的估測(cè)單產(chǎn),結(jié)合歷年《陜西省統(tǒng)計(jì)年鑒》的單產(chǎn)統(tǒng)計(jì)數(shù)據(jù)對(duì)估產(chǎn)結(jié)果的精度進(jìn)行評(píng)價(jià)(表2)。結(jié)果表明,關(guān)中平原各市的單產(chǎn)估測(cè)結(jié)果與年鑒統(tǒng)計(jì)的單產(chǎn)較均為接近。GA獲取的單產(chǎn)的RMSE和RE均較低且皆小于CAFE,與EA的結(jié)果相近或相等。各市及4市綜合的估測(cè)單產(chǎn)的RMSE均小于300 kg/hm2,RE除4市綜合外均小于5.0%,其中咸陽(yáng)和西安的RMSE僅分別為102.98 kg/hm2和123.64 kg/hm2,RE僅分別為1.95%和1.93%。GA與EA相比RMSE的最大相差為0.18 kg/hm2,RE最大相差為0.02%。進(jìn)一步表明了采用GA獲取的關(guān)中平原冬小麥生育時(shí)期的干旱對(duì)產(chǎn)量的影響的權(quán)重較為合理,應(yīng)用于估產(chǎn)的精度較高,GA更適合于基于多時(shí)相VTCI數(shù)據(jù)的時(shí)間尺度轉(zhuǎn)換的關(guān)中平原干旱影響評(píng)估。

    表2 2008-2014年冬小麥估測(cè)單產(chǎn)與統(tǒng)計(jì)單產(chǎn)對(duì)比

    以上結(jié)果均表明采用GA求解冬小麥生育時(shí)期優(yōu)化賦權(quán)問(wèn)題優(yōu)于EA和CAFE,其結(jié)果應(yīng)用于估產(chǎn)的精度較高。估算分析干旱年2013年和2014年關(guān)中平原各市冬小麥的產(chǎn)量(據(jù)氣象部門報(bào)道,2013年陜西省遭遇了自1961年以來(lái)最嚴(yán)重的一次氣候干旱,所以采用2013年和2014年產(chǎn)量進(jìn)行精度評(píng)價(jià)),并與年鑒統(tǒng)計(jì)單產(chǎn)對(duì)比分析進(jìn)行精度評(píng)價(jià),見(jiàn)表3。結(jié)果表明,除2014年寶雞市估測(cè)結(jié)果外,各市的單產(chǎn)估測(cè)結(jié)果與年鑒統(tǒng)計(jì)的單產(chǎn)較為接近,精度較高。2013年各市估測(cè)單產(chǎn)與統(tǒng)計(jì)產(chǎn)量的絕對(duì)誤差的值域范圍為12.32~120.95 kg/hm2,相對(duì)誤差的值域范圍為0.32%~3.66%,精度均較高。2014年各市絕對(duì)誤差的值域范圍為7.70~440.93 kg/hm2,相對(duì)誤差的值域范圍為0.22%~10.46%,除寶雞市外,渭南市、咸陽(yáng)市和西安市的精度均較高。這是由于2014年3月中、下旬和4月上旬寶雞市地區(qū)發(fā)生了較嚴(yán)重的干旱,對(duì)處于返青期末期和拔節(jié)期的冬小麥的生長(zhǎng)產(chǎn)生了嚴(yán)重的影響,根、莖和葉生長(zhǎng)速率降低且有效小穗數(shù)減少,雖然抽穗—灌漿期和乳熟期水分條件較好,但仍對(duì)產(chǎn)量產(chǎn)生了嚴(yán)重的影響。綜上所述GA應(yīng)用于估產(chǎn)的精度較高且魯棒性較好,適合于基于多時(shí)相VTCI數(shù)據(jù)的時(shí)間尺度轉(zhuǎn)換的關(guān)中平原干旱影響評(píng)估。

    表3 基于遺傳算法的2013年、2014年冬小麥估測(cè)單產(chǎn)與統(tǒng)計(jì)單產(chǎn)對(duì)比

    3 討 論

    3種方法中EA可得到最優(yōu)解,GA和CAFE只能得到近似最優(yōu)解或滿意解。CAFE兼顧了主觀賦權(quán)法的主觀偏好和減小了客觀賦權(quán)法的客觀隨意性,其經(jīng)時(shí)間尺度轉(zhuǎn)換后獲得的加權(quán)VTCI除渭南市外均與冬小麥單產(chǎn)的相關(guān)性較高,但CAFE確定的權(quán)重值與EA確定的最優(yōu)權(quán)重值相差較大,且其經(jīng)時(shí)間尺度后獲得的加權(quán)VTCI與冬小麥單產(chǎn)的回歸分析結(jié)果精度較差,這可能是由于CAFE中因子權(quán)重排序法將返青期、拔節(jié)期、抽穗—灌漿期和乳熟期干旱對(duì)產(chǎn)量的影響權(quán)重范圍分別設(shè)置為0~0.1,0.4~0.5,0.3~0.4,0~0.2[10],區(qū)間范圍較小,賦權(quán)過(guò)程中不能遍歷到所有結(jié)果,同時(shí),熵值法不能體現(xiàn)不同生育時(shí)期的干旱對(duì)小麥的生長(zhǎng)發(fā)育的影響的不同程度,使得此方法只能獲得滿意解,結(jié)果精度較差。雖然EA可以得到最優(yōu)的權(quán)重,從而獲得最準(zhǔn)確的能夠反映關(guān)中平原冬小麥不同生育時(shí)期干旱對(duì)冬小麥生產(chǎn)及其單產(chǎn)的影響的權(quán)重值以進(jìn)行時(shí)間尺度轉(zhuǎn)換研究,且運(yùn)算簡(jiǎn)單,但是EA由于需要遍歷每一組可能的權(quán)重值尋找最優(yōu)結(jié)果,其運(yùn)算復(fù)雜度隨權(quán)重范圍呈指數(shù)增長(zhǎng),計(jì)算量大而使得有效性差,所以EA多用于評(píng)價(jià)其他方法的參考。而GA將問(wèn)題參數(shù)編碼為染色體后進(jìn)行優(yōu)化,從一個(gè)種群開(kāi)始,不斷地產(chǎn)生和測(cè)試新一代的種群,擴(kuò)大了搜索范圍,且可以利用染色體長(zhǎng)度控制數(shù)據(jù)精度,可較精確、快速地完成最優(yōu)權(quán)重值優(yōu)化問(wèn)題的求解,GA能較好地解決了EA和CAFE的限制問(wèn)題,在多時(shí)相VTCI數(shù)據(jù)時(shí)間尺度轉(zhuǎn)換的實(shí)際應(yīng)用中更簡(jiǎn)單、有效、適用。

    本文研究發(fā)現(xiàn)渭南市3種賦權(quán)方法的結(jié)果均未通過(guò)F檢驗(yàn),可能是渭南市旱作地較多,發(fā)生干旱后恢復(fù)能力差[4],同時(shí),近些年果業(yè)的發(fā)展,大量農(nóng)田變?yōu)閳@地或果樹與農(nóng)作物間種,使得3種賦權(quán)方法確定的權(quán)重值獲得尺度轉(zhuǎn)換后的VTCI與冬小麥產(chǎn)量的相關(guān)性較差。

    返青期適度旱后及時(shí)復(fù)水不會(huì)造成產(chǎn)量較低,反而由于復(fù)水補(bǔ)償效應(yīng)促進(jìn)冬小麥的生長(zhǎng)發(fā)育和收獲指數(shù)的提高,說(shuō)明返青期發(fā)生一定程度內(nèi)干旱后適度供水,有利于冬小麥產(chǎn)量和水分利用效率的提高,實(shí)現(xiàn)水資源的高效利用和指導(dǎo)農(nóng)業(yè)生產(chǎn)。不同生育期干旱對(duì)產(chǎn)量的影響權(quán)重值的大小能夠反映各生育時(shí)期VTCI數(shù)據(jù)所攜帶干旱對(duì)產(chǎn)量影響程度的信息,可依此權(quán)重值大小來(lái)確定能有效反映旱情對(duì)冬小麥產(chǎn)量影響的敏感VTCI時(shí)相,以指導(dǎo)田間管理和選擇少數(shù)關(guān)鍵遙感時(shí)相開(kāi)展干旱影響評(píng)估和冬小麥估產(chǎn)。同時(shí),考慮冬小麥生長(zhǎng)機(jī)理及非干旱因素對(duì)冬小麥產(chǎn)量的影響,更加準(zhǔn)確合理的確定冬小麥各主要生育時(shí)期的權(quán)重,結(jié)合氣候特征、人類活動(dòng)的影響等以得到更加科學(xué)、合理的干旱對(duì)農(nóng)業(yè)生產(chǎn)的影響評(píng)估。

    4 結(jié) 論

    作物長(zhǎng)勢(shì)和最終產(chǎn)量與其主要生育期發(fā)生的干旱程度密切相關(guān),VTCI是綜合地表主要參數(shù)歸一化植被指數(shù)和地表溫度的干旱指標(biāo),能夠準(zhǔn)確地監(jiān)測(cè)旱情狀況。本文針對(duì)關(guān)中平原冬小麥主要生育期VTCI的時(shí)間尺度轉(zhuǎn)化問(wèn)題進(jìn)行了建模分析,研究了求解該問(wèn)題的CAFE,EA和GA 3種方法獲取的權(quán)重值。3種方法獲取的冬小麥返青期、拔節(jié)期、抽穗—灌漿期和乳熟期干旱對(duì)產(chǎn)量的影響的權(quán)重值均較為合理。為提升求解效率,以EA計(jì)算得的全局最優(yōu)解為標(biāo)準(zhǔn),對(duì)比分析了3種賦權(quán)方法,CAFE確定的干旱對(duì)產(chǎn)量影響的權(quán)重值與最優(yōu)權(quán)重值相差較大,且其經(jīng)時(shí)間尺度轉(zhuǎn)換后獲得的VTCI與冬小麥單產(chǎn)的回歸分析精度整體低于EA和GA。GA獲得的權(quán)重與最優(yōu)權(quán)重值基本吻合,精度均較高,同時(shí),GA獲得的關(guān)中平原各市的單產(chǎn)估測(cè)結(jié)果與年鑒統(tǒng)計(jì)的單產(chǎn)較為接近,RMSE和RE均較低。GA接近于EA的精度和性能,而時(shí)間復(fù)雜度大大低于EA。GA能夠?yàn)榛诙鄷r(shí)相VTCI數(shù)據(jù)時(shí)間尺度轉(zhuǎn)換的干旱影響評(píng)估研究提供有效的賦權(quán)方法。

    參考文獻(xiàn):

    [1] Kongan F N. Application of vegetation index and brightness temperature for drought detection [J]. Advances in Space Research, 1995,15(11):91-100.

    [2] 王鵬新,龔健雅,李小文.條件植被溫度指數(shù)及其在干旱監(jiān)測(cè)中的應(yīng)用[J].武漢大學(xué)學(xué)報(bào):信息科學(xué)版,2001,26(5):412-418.

    [3] 白雪嬌,王鵬新,解毅,等.基于結(jié)構(gòu)相似度的關(guān)中平原旱情空間分布特征[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2015,46(11):348-354.

    [4] 李夢(mèng)娜,錢會(huì),喬亮.關(guān)中地區(qū)農(nóng)業(yè)干旱脆弱性評(píng)價(jià)[J].資源科學(xué),2016,38(1):155-165.

    [5] 馮海霞,秦其明,蔣洪波,等.基于HJ-1A/1BCCD數(shù)據(jù)的干旱監(jiān)測(cè)[J].農(nóng)業(yè)工程學(xué)報(bào),2011,27(S1):358-365.

    [6] Ghulam A, Li Z L, Qin Q M, et al. Estimating crop water stress with ETM+NIR and SWIR data [J]. Agricultural and Forest Meteorology, 2008,148(11):1679-1695.

    [7] 曹陽(yáng),楊婕,熊偉,等.1962—2010年潛在干旱對(duì)中國(guó)冬小麥產(chǎn)量影響的模擬分析[J].農(nóng)業(yè)工程學(xué)報(bào),2014(7):128-139.

    [8] Patel N R, Parida B R, Venus V, et al. Analysis of agricultural drought using vegetation temperature condition index(VTCI)from Terra/MODIS satellite data [J]. Environmental Monitoring & Assessment, 2011,184(12):7153-7163.

    [9] Sun W, Wang P X, Zhang S Y, et al. Using the vegetation temperature condition index for time series drought occurrence monitoring in the Guanzhong Plain, PR China[J]. International Journal of Remote Sensing, 2008,29(17/18):5133-5144.

    [10] 王蕾.基于條件植被溫度指數(shù)的冬小麥生育時(shí)期干旱影響評(píng)估研究[D].北京:中國(guó)農(nóng)業(yè)大學(xué)信息與電氣工程學(xué)院,2015.

    [11] 田苗,王鵬新,張樹譽(yù),等.基于條件植被溫度指數(shù)的冬小麥產(chǎn)量預(yù)測(cè)[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2014,45(2):239-245.

    [12] Foulkes M J, Sylvester-Bradley R, Weightman R, et al. Identifying physiological traits associated with improved resistance in winter wheat[J]. Field Crops Research, 2007,103(1):11-24.

    [13] 陳華友.多屬性決策中基于離差最大化的組合賦權(quán)方法[J].系統(tǒng)工程與電子技術(shù),2004,26(2):194-197.

    [14] 倪九派,李萍,魏朝富,等.基于AHP和熵權(quán)法賦權(quán)的區(qū)域土地開(kāi)發(fā)整理潛力評(píng)價(jià)[J].農(nóng)業(yè)工程學(xué)報(bào),2009,25(5):202-209.

    [15] 王書吉,費(fèi)良軍,雷雁斌,等.綜合集成賦權(quán)法在灌區(qū)節(jié)水改造效益評(píng)價(jià)中的應(yīng)用[J].農(nóng)業(yè)工程學(xué)報(bào),2008,24(12):48-51.

    [16] 遲道才,馬濤,李松.基于博弈論的可拓評(píng)價(jià)方法在灌區(qū)運(yùn)行狀況評(píng)價(jià)中的應(yīng)用[J].農(nóng)業(yè)工程學(xué)報(bào),2008,24(8):36-39.

    [17] 曹永強(qiáng),李香云,馬靜,等.基于可變模糊算法的大連市農(nóng)業(yè)干旱風(fēng)險(xiǎn)評(píng)價(jià)[J].資源科學(xué),2011,33(5):983-988.

    [18] 黃國(guó)強(qiáng),陸宇平,南英,等.多目標(biāo)連續(xù)小推力深空探測(cè)器軌道全局優(yōu)化[J].系統(tǒng)工程與電子技術(shù),2012,34(8):1652-1659.

    [19] Lin W, Delgadofirs Y G, Gause D C, et al. Hybrid newton-raphson genetic algorithm for the traveling salesman problem[J]. Cybernetics and Systems, 1995,26(4):387-412.

    [20] Park B J, Choi H R, Kim H S. A hybrid genetic algorithm for the job shop scheduling problems [J]. Computers & Industrial Engineering, 2003,45(4):597-613.

    [21] 祁有玲,張富倉(cāng),李開(kāi)峰,等.不同生育期水分虧缺及氮營(yíng)養(yǎng)對(duì)冬小麥生長(zhǎng)和產(chǎn)量的影響[J].灌溉排水學(xué)報(bào),2009,28(1):24-27.

    猜你喜歡
    關(guān)中平原時(shí)間尺度單產(chǎn)
    時(shí)間尺度上非完整系統(tǒng)的Noether準(zhǔn)對(duì)稱性與守恒量
    關(guān)中平原人為土形成的歷史探析
    時(shí)間尺度上Lagrange 系統(tǒng)的Hojman 守恒量1)
    農(nóng)大農(nóng)企聯(lián)手創(chuàng)山西小麥最高單產(chǎn)新紀(jì)錄
    油菜“不務(wù)正業(yè)”,單產(chǎn)3.4噸
    交直流混合微電網(wǎng)多時(shí)間尺度協(xié)同控制
    能源工程(2021年1期)2021-04-13 02:06:12
    單產(chǎn)948.48千克!“金種子”迸發(fā)大能量
    我國(guó)玉米單產(chǎn)紀(jì)錄第七次被刷新
    基于三角模型的關(guān)中5市土地生態(tài)安全區(qū)域差異分析
    《關(guān)中平原城市群發(fā)展規(guī)劃》獲批發(fā)布
    新西部(2018年3期)2018-03-21 10:09:08
    在线观看66精品国产| 精品一区二区三区四区五区乱码| 国产成人影院久久av| 国产区一区二久久| 美国免费a级毛片| 热re99久久精品国产66热6| 丁香欧美五月| 亚洲国产精品一区二区三区在线| 婷婷丁香在线五月| 在线观看舔阴道视频| 纯流量卡能插随身wifi吗| 深夜精品福利| 久久午夜亚洲精品久久| xxxhd国产人妻xxx| 蜜桃国产av成人99| 国产日韩欧美亚洲二区| 国产在线一区二区三区精| 在线亚洲精品国产二区图片欧美| 日韩欧美国产一区二区入口| www.自偷自拍.com| 一级a爱视频在线免费观看| 精品国产超薄肉色丝袜足j| 中文字幕制服av| 国产一区有黄有色的免费视频| 亚洲成人手机| 欧美+亚洲+日韩+国产| 一本一本久久a久久精品综合妖精| 老鸭窝网址在线观看| 婷婷丁香在线五月| www日本在线高清视频| 波多野结衣av一区二区av| 免费高清在线观看日韩| 精品一区二区三区四区五区乱码| 91精品国产国语对白视频| 国产日韩欧美视频二区| 亚洲精品美女久久久久99蜜臀| 欧美日韩成人在线一区二区| 伦理电影免费视频| 亚洲精品在线美女| 日韩欧美三级三区| 中文字幕人妻熟女乱码| 国产精品电影一区二区三区 | 亚洲欧美激情在线| 另类亚洲欧美激情| 成年人黄色毛片网站| 精品第一国产精品| 国产日韩欧美视频二区| 成年人免费黄色播放视频| 国产成人影院久久av| 欧美大码av| 在线观看66精品国产| av视频免费观看在线观看| 亚洲国产欧美一区二区综合| 黑人欧美特级aaaaaa片| 亚洲国产成人一精品久久久| 国产亚洲欧美在线一区二区| 操出白浆在线播放| 午夜福利影视在线免费观看| 一区二区日韩欧美中文字幕| 亚洲五月婷婷丁香| 黄片播放在线免费| 黄色 视频免费看| 一二三四社区在线视频社区8| 久久午夜亚洲精品久久| 亚洲av片天天在线观看| 国产xxxxx性猛交| 国产免费福利视频在线观看| 欧美日韩亚洲国产一区二区在线观看 | 视频在线观看一区二区三区| 久久99热这里只频精品6学生| 亚洲精品国产区一区二| 色94色欧美一区二区| 后天国语完整版免费观看| 国产精品免费大片| 青青草视频在线视频观看| 怎么达到女性高潮| 免费一级毛片在线播放高清视频 | 色94色欧美一区二区| 免费在线观看黄色视频的| 中文字幕av电影在线播放| 亚洲国产成人一精品久久久| 免费观看av网站的网址| 亚洲精品国产色婷婷电影| 在线观看免费视频日本深夜| 99精品欧美一区二区三区四区| 成人手机av| 看免费av毛片| 亚洲av成人不卡在线观看播放网| 大陆偷拍与自拍| 日韩大片免费观看网站| 在线观看免费视频日本深夜| av有码第一页| 女人精品久久久久毛片| 一二三四在线观看免费中文在| 日韩欧美一区二区三区在线观看 | 操出白浆在线播放| 欧美精品啪啪一区二区三区| 免费在线观看完整版高清| 国产av又大| 黄色毛片三级朝国网站| 国产一区有黄有色的免费视频| 高清黄色对白视频在线免费看| 亚洲avbb在线观看| 午夜两性在线视频| 国产男女超爽视频在线观看| 国产男女内射视频| 老司机靠b影院| 女人久久www免费人成看片| 精品一区二区三区四区五区乱码| 欧美精品亚洲一区二区| 久久久水蜜桃国产精品网| 女警被强在线播放| 岛国在线观看网站| 中文字幕人妻丝袜一区二区| 十分钟在线观看高清视频www| 色老头精品视频在线观看| 亚洲精品久久成人aⅴ小说| 国产免费现黄频在线看| 久久人妻熟女aⅴ| 人成视频在线观看免费观看| 亚洲天堂av无毛| 亚洲成a人片在线一区二区| 久久久久精品国产欧美久久久| 91精品国产国语对白视频| 激情视频va一区二区三区| 19禁男女啪啪无遮挡网站| 性少妇av在线| 麻豆乱淫一区二区| 99国产综合亚洲精品| 欧美日韩亚洲国产一区二区在线观看 | 色播在线永久视频| 精品人妻在线不人妻| 国产黄频视频在线观看| 热re99久久国产66热| 在线观看一区二区三区激情| tocl精华| 啦啦啦在线免费观看视频4| 欧美精品亚洲一区二区| 人人妻,人人澡人人爽秒播| 欧美变态另类bdsm刘玥| 国产欧美亚洲国产| 免费看十八禁软件| 亚洲情色 制服丝袜| 999久久久国产精品视频| 日韩欧美免费精品| 50天的宝宝边吃奶边哭怎么回事| 欧美精品高潮呻吟av久久| 久久久久久久久久久久大奶| 免费一级毛片在线播放高清视频 | 久久久久精品国产欧美久久久| 色婷婷久久久亚洲欧美| 人妻 亚洲 视频| 91国产中文字幕| 国产1区2区3区精品| 天天躁狠狠躁夜夜躁狠狠躁| 欧美日韩视频精品一区| 又黄又粗又硬又大视频| 天堂动漫精品| 91麻豆av在线| 女人被躁到高潮嗷嗷叫费观| 精品免费久久久久久久清纯 | 精品一区二区三区视频在线观看免费 | 一区二区三区国产精品乱码| 丁香欧美五月| tocl精华| 一本久久精品| 欧美黄色片欧美黄色片| 可以免费在线观看a视频的电影网站| 国产精品一区二区精品视频观看| 性高湖久久久久久久久免费观看| 伦理电影免费视频| 一个人免费看片子| 最黄视频免费看| 欧美 亚洲 国产 日韩一| 欧美老熟妇乱子伦牲交| 亚洲第一欧美日韩一区二区三区 | 母亲3免费完整高清在线观看| 欧美久久黑人一区二区| 99re在线观看精品视频| 精品一区二区三区四区五区乱码| 精品一区二区三区视频在线观看免费 | 欧美日韩精品网址| 午夜福利在线免费观看网站| www.熟女人妻精品国产| 欧美+亚洲+日韩+国产| 男女下面插进去视频免费观看| 淫妇啪啪啪对白视频| 亚洲全国av大片| 男女下面插进去视频免费观看| 午夜免费鲁丝| 两个人免费观看高清视频| 99热国产这里只有精品6| 亚洲va日本ⅴa欧美va伊人久久| 宅男免费午夜| 婷婷丁香在线五月| 乱人伦中国视频| 日韩人妻精品一区2区三区| 日日夜夜操网爽| 丁香欧美五月| 精品一区二区三区四区五区乱码| 999久久久精品免费观看国产| 国产欧美日韩精品亚洲av| 在线观看免费日韩欧美大片| 老司机福利观看| 久久国产亚洲av麻豆专区| 亚洲成人免费电影在线观看| 日韩视频在线欧美| 天天添夜夜摸| 美女主播在线视频| 午夜免费鲁丝| 亚洲精品国产区一区二| 亚洲精品一二三| 欧美日韩福利视频一区二区| 一边摸一边做爽爽视频免费| 香蕉丝袜av| 人人澡人人妻人| 天天躁狠狠躁夜夜躁狠狠躁| 99热国产这里只有精品6| 9热在线视频观看99| 大香蕉久久网| 91大片在线观看| 亚洲少妇的诱惑av| 欧美人与性动交α欧美精品济南到| 757午夜福利合集在线观看| 亚洲精品美女久久av网站| 熟女少妇亚洲综合色aaa.| 大型黄色视频在线免费观看| 欧美精品一区二区免费开放| 我要看黄色一级片免费的| www.999成人在线观看| 中文字幕色久视频| 12—13女人毛片做爰片一| 亚洲av日韩在线播放| 黄色怎么调成土黄色| 精品一区二区三卡| 女人精品久久久久毛片| 国产无遮挡羞羞视频在线观看| 性色av乱码一区二区三区2| 青青草视频在线视频观看| 久久国产精品男人的天堂亚洲| 一区在线观看完整版| 亚洲一区中文字幕在线| 精品国产乱码久久久久久男人| 人人妻人人添人人爽欧美一区卜| 国产精品国产av在线观看| 最新在线观看一区二区三区| 狠狠精品人妻久久久久久综合| 最近最新免费中文字幕在线| 免费日韩欧美在线观看| 国产av精品麻豆| 久久久久久久久免费视频了| 午夜福利一区二区在线看| 性少妇av在线| 好男人电影高清在线观看| 免费看十八禁软件| 欧美日韩av久久| 国产精品一区二区在线观看99| 国产人伦9x9x在线观看| 精品国产一区二区三区四区第35| 大香蕉久久成人网| 国产精品一区二区免费欧美| 999久久久精品免费观看国产| a级片在线免费高清观看视频| 精品国产一区二区三区久久久樱花| 少妇精品久久久久久久| 热99久久久久精品小说推荐| 精品免费久久久久久久清纯 | 1024视频免费在线观看| 丰满人妻熟妇乱又伦精品不卡| 在线观看66精品国产| 久久精品人人爽人人爽视色| 日本av手机在线免费观看| 人妻 亚洲 视频| 国产精品亚洲av一区麻豆| 久久久久久久大尺度免费视频| 大型av网站在线播放| 久久午夜综合久久蜜桃| 天堂俺去俺来也www色官网| 久久久精品免费免费高清| 亚洲七黄色美女视频| 精品久久蜜臀av无| 91九色精品人成在线观看| 乱人伦中国视频| 亚洲性夜色夜夜综合| av片东京热男人的天堂| 国产精品久久久久成人av| 日韩精品免费视频一区二区三区| 亚洲国产欧美一区二区综合| 欧美 亚洲 国产 日韩一| 视频区欧美日本亚洲| 国产在线精品亚洲第一网站| 18禁国产床啪视频网站| 一级毛片电影观看| 午夜福利乱码中文字幕| 久久久久国产一级毛片高清牌| av线在线观看网站| 视频区图区小说| xxxhd国产人妻xxx| 欧美黑人精品巨大| 国产成人av教育| 热re99久久精品国产66热6| 国产精品 国内视频| 欧美 亚洲 国产 日韩一| 国产精品一区二区在线观看99| 黄色怎么调成土黄色| 国产一卡二卡三卡精品| 国产高清国产精品国产三级| 少妇被粗大的猛进出69影院| 老司机福利观看| 午夜福利视频在线观看免费| 国产成人免费观看mmmm| 涩涩av久久男人的天堂| 成人国产av品久久久| 黄片大片在线免费观看| 亚洲欧美激情在线| 免费人妻精品一区二区三区视频| 老司机在亚洲福利影院| 国产精品久久久久久人妻精品电影 | kizo精华| 国产精品一区二区在线不卡| 热99久久久久精品小说推荐| 国产人伦9x9x在线观看| 亚洲成人国产一区在线观看| 丰满迷人的少妇在线观看| 国产av一区二区精品久久| 老司机福利观看| 国产成人系列免费观看| 精品久久蜜臀av无| 国产又色又爽无遮挡免费看| 人妻久久中文字幕网| 少妇精品久久久久久久| 777米奇影视久久| tube8黄色片| 日本vs欧美在线观看视频| 久久人人97超碰香蕉20202| 免费观看a级毛片全部| 一区二区三区激情视频| 成人手机av| 亚洲中文字幕日韩| 中文欧美无线码| 国产精品亚洲一级av第二区| 汤姆久久久久久久影院中文字幕| 欧美成人午夜精品| 青草久久国产| 美女视频免费永久观看网站| 日韩有码中文字幕| 久久人人97超碰香蕉20202| av福利片在线| 自线自在国产av| 女警被强在线播放| av天堂在线播放| 18禁黄网站禁片午夜丰满| 国产成人欧美| 最黄视频免费看| 欧美日韩亚洲高清精品| 一区在线观看完整版| 亚洲精品中文字幕在线视频| 两性夫妻黄色片| 妹子高潮喷水视频| 最黄视频免费看| 老熟妇仑乱视频hdxx| 亚洲av第一区精品v没综合| 久久青草综合色| 国产在线一区二区三区精| 久久久久久久精品吃奶| 亚洲综合色网址| 欧美精品亚洲一区二区| 久久av网站| 一区二区av电影网| 久久国产精品影院| 成人免费观看视频高清| 国产亚洲欧美精品永久| 日韩免费av在线播放| 亚洲精品自拍成人| 成人特级黄色片久久久久久久 | 亚洲专区字幕在线| 国产福利在线免费观看视频| 色老头精品视频在线观看| 纵有疾风起免费观看全集完整版| 一级a爱视频在线免费观看| 1024视频免费在线观看| 妹子高潮喷水视频| 99热国产这里只有精品6| 亚洲专区中文字幕在线| 汤姆久久久久久久影院中文字幕| 国产99久久九九免费精品| 国产单亲对白刺激| 热99国产精品久久久久久7| 91麻豆精品激情在线观看国产 | 国产精品国产高清国产av | 狠狠婷婷综合久久久久久88av| 久久99热这里只频精品6学生| 考比视频在线观看| 成年女人毛片免费观看观看9 | 亚洲精品在线美女| 精品福利永久在线观看| 欧美乱码精品一区二区三区| 国产精品二区激情视频| www.熟女人妻精品国产| 久久国产精品男人的天堂亚洲| 久久这里只有精品19| 老汉色∧v一级毛片| 亚洲av国产av综合av卡| 岛国在线观看网站| 久久中文字幕一级| 淫妇啪啪啪对白视频| 亚洲欧美色中文字幕在线| 国产极品粉嫩免费观看在线| 亚洲熟女毛片儿| 国产亚洲欧美在线一区二区| 精品久久久久久电影网| 少妇粗大呻吟视频| 国产有黄有色有爽视频| 国产精品久久电影中文字幕 | 久久久久久人人人人人| 一区二区三区精品91| 一边摸一边抽搐一进一出视频| 国产免费现黄频在线看| 中文字幕最新亚洲高清| 在线观看免费视频网站a站| 久久久精品国产亚洲av高清涩受| 97人妻天天添夜夜摸| 久久中文字幕人妻熟女| 国产免费视频播放在线视频| 国产一区二区三区综合在线观看| 国产精品一区二区在线不卡| 国产xxxxx性猛交| 亚洲精品乱久久久久久| 久久人人97超碰香蕉20202| 欧美精品一区二区免费开放| 黄色视频在线播放观看不卡| 19禁男女啪啪无遮挡网站| 国产日韩一区二区三区精品不卡| 色在线成人网| 欧美日韩福利视频一区二区| 亚洲成人免费电影在线观看| 大型av网站在线播放| 丰满少妇做爰视频| 国产欧美亚洲国产| 精品第一国产精品| 国产av一区二区精品久久| 成人三级做爰电影| 亚洲少妇的诱惑av| 男人操女人黄网站| 成年人午夜在线观看视频| av有码第一页| 两个人看的免费小视频| 女人高潮潮喷娇喘18禁视频| 精品一区二区三卡| 狂野欧美激情性xxxx| 久久精品熟女亚洲av麻豆精品| 中文字幕人妻丝袜制服| 人人妻人人添人人爽欧美一区卜| 精品一品国产午夜福利视频| 极品人妻少妇av视频| 老司机靠b影院| 丝袜人妻中文字幕| 一边摸一边抽搐一进一出视频| 国产精品av久久久久免费| 少妇猛男粗大的猛烈进出视频| 高清视频免费观看一区二区| 在线播放国产精品三级| h视频一区二区三区| 中文字幕色久视频| 午夜免费成人在线视频| 精品熟女少妇八av免费久了| 日本黄色视频三级网站网址 | 色婷婷av一区二区三区视频| 精品一品国产午夜福利视频| 看免费av毛片| 男女之事视频高清在线观看| 午夜两性在线视频| 菩萨蛮人人尽说江南好唐韦庄| 久久国产精品男人的天堂亚洲| 欧美大码av| 日韩熟女老妇一区二区性免费视频| 我要看黄色一级片免费的| 麻豆av在线久日| 久久久久久人人人人人| 大码成人一级视频| 欧美激情 高清一区二区三区| 男男h啪啪无遮挡| 高清黄色对白视频在线免费看| 建设人人有责人人尽责人人享有的| 人人妻人人澡人人看| 欧美在线黄色| 亚洲欧美日韩高清在线视频 | 麻豆成人av在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 一个人免费看片子| 精品第一国产精品| 俄罗斯特黄特色一大片| 黑人巨大精品欧美一区二区mp4| 欧美日韩一级在线毛片| 欧美在线黄色| 久久99热这里只频精品6学生| 中文字幕另类日韩欧美亚洲嫩草| 欧美成人午夜精品| 性少妇av在线| 久久久欧美国产精品| 丁香六月欧美| 精品福利永久在线观看| 久久久久久久大尺度免费视频| 成年女人毛片免费观看观看9 | 亚洲性夜色夜夜综合| 国产精品98久久久久久宅男小说| 中文字幕人妻熟女乱码| 日韩 欧美 亚洲 中文字幕| 50天的宝宝边吃奶边哭怎么回事| 狠狠婷婷综合久久久久久88av| 亚洲三区欧美一区| 99久久精品国产亚洲精品| 亚洲美女黄片视频| 首页视频小说图片口味搜索| 国产老妇伦熟女老妇高清| 黄色视频在线播放观看不卡| 久久精品国产亚洲av高清一级| 欧美性长视频在线观看| 精品免费久久久久久久清纯 | 真人做人爱边吃奶动态| 国产男女超爽视频在线观看| 99国产精品一区二区三区| 亚洲一码二码三码区别大吗| 在线观看66精品国产| 亚洲av国产av综合av卡| 国产亚洲av高清不卡| 久久天躁狠狠躁夜夜2o2o| 97人妻天天添夜夜摸| 亚洲av欧美aⅴ国产| 夜夜骑夜夜射夜夜干| 日本a在线网址| 黑人操中国人逼视频| 久久中文字幕人妻熟女| 亚洲专区字幕在线| 亚洲国产欧美一区二区综合| 欧美+亚洲+日韩+国产| svipshipincom国产片| 成人亚洲精品一区在线观看| 欧美精品亚洲一区二区| 国产麻豆69| 日韩欧美一区二区三区在线观看 | 精品福利观看| 黄色片一级片一级黄色片| 国产av一区二区精品久久| 男男h啪啪无遮挡| 色精品久久人妻99蜜桃| 老熟女久久久| 91大片在线观看| 欧美激情高清一区二区三区| 热re99久久精品国产66热6| 成年女人毛片免费观看观看9 | 亚洲av第一区精品v没综合| 婷婷成人精品国产| 2018国产大陆天天弄谢| 日韩熟女老妇一区二区性免费视频| 97人妻天天添夜夜摸| 久久这里只有精品19| 久久天堂一区二区三区四区| 午夜福利影视在线免费观看| 欧美日韩视频精品一区| av欧美777| 丝袜美足系列| 一区二区三区精品91| 国产在线观看jvid| 日日爽夜夜爽网站| 1024视频免费在线观看| 亚洲成人免费av在线播放| 50天的宝宝边吃奶边哭怎么回事| 亚洲一码二码三码区别大吗| 欧美+亚洲+日韩+国产| av天堂久久9| 丝袜美腿诱惑在线| 婷婷成人精品国产| 欧美中文综合在线视频| 国产男女超爽视频在线观看| 亚洲精品在线美女| 久久久国产一区二区| 色94色欧美一区二区| 亚洲三区欧美一区| 久久久久国产一级毛片高清牌| 欧美激情极品国产一区二区三区| 欧美精品亚洲一区二区| 国产成人免费观看mmmm| 老汉色av国产亚洲站长工具| 国产精品98久久久久久宅男小说| 两人在一起打扑克的视频| 99re6热这里在线精品视频| 亚洲午夜精品一区,二区,三区| 日本wwww免费看| 欧美激情 高清一区二区三区| 欧美国产精品va在线观看不卡| 日本一区二区免费在线视频| 欧美变态另类bdsm刘玥| 我要看黄色一级片免费的| 欧美日韩一级在线毛片| 欧美午夜高清在线| 我要看黄色一级片免费的| 在线观看免费午夜福利视频| 久久久欧美国产精品| 欧美 亚洲 国产 日韩一| 日本一区二区免费在线视频| 亚洲成人国产一区在线观看| 在线 av 中文字幕| 久久久精品国产亚洲av高清涩受| 人妻一区二区av| 国产老妇伦熟女老妇高清| 菩萨蛮人人尽说江南好唐韦庄| 亚洲久久久国产精品| av免费在线观看网站| 一本—道久久a久久精品蜜桃钙片| 国产午夜精品久久久久久| 欧美成人午夜精品| 一级片'在线观看视频| 亚洲精品美女久久av网站| 夜夜骑夜夜射夜夜干|