• <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
    最近最新中文字幕大全电影3| 国产三级在线视频| 舔av片在线| 亚洲欧美清纯卡通| 国产精品伦人一区二区| 色综合婷婷激情| 熟女电影av网| 国产精品av视频在线免费观看| 免费看美女性在线毛片视频| 成人性生交大片免费视频hd| 十八禁国产超污无遮挡网站| 夜夜看夜夜爽夜夜摸| 色在线成人网| 免费看美女性在线毛片视频| 99热这里只有是精品在线观看| 午夜精品一区二区三区免费看| 91狼人影院| 精品人妻偷拍中文字幕| 长腿黑丝高跟| 99久久精品国产国产毛片| 亚洲avbb在线观看| 婷婷精品国产亚洲av| 91av网一区二区| 国产精品自产拍在线观看55亚洲| 亚洲国产精品成人综合色| 最近视频中文字幕2019在线8| 国产一区二区三区视频了| 看黄色毛片网站| 99热精品在线国产| 两人在一起打扑克的视频| 中文字幕熟女人妻在线| 精品乱码久久久久久99久播| 搡女人真爽免费视频火全软件 | 乱系列少妇在线播放| 日韩一本色道免费dvd| 色在线成人网| 美女黄网站色视频| 久久久久久久午夜电影| 国内毛片毛片毛片毛片毛片| 丰满人妻一区二区三区视频av| 搡女人真爽免费视频火全软件 | 两性午夜刺激爽爽歪歪视频在线观看| 国产精品人妻久久久久久| 真实男女啪啪啪动态图| 精品一区二区三区人妻视频| 日日摸夜夜添夜夜添av毛片 | 午夜免费男女啪啪视频观看 | 精品一区二区三区av网在线观看| 人妻久久中文字幕网| 97热精品久久久久久| 变态另类丝袜制服| 美女高潮喷水抽搐中文字幕| АⅤ资源中文在线天堂| 亚洲av日韩精品久久久久久密| 舔av片在线| 日本免费a在线| 国产一区二区三区在线臀色熟女| 99国产精品一区二区蜜桃av| 久久久久国产精品人妻aⅴ院| 如何舔出高潮| 国产精品伦人一区二区| 我要搜黄色片| 免费大片18禁| 少妇裸体淫交视频免费看高清| 国产精品一区二区性色av| 国产主播在线观看一区二区| 深爱激情五月婷婷| 又紧又爽又黄一区二区| 国产精品综合久久久久久久免费| 少妇的逼水好多| 别揉我奶头 嗯啊视频| 国产久久久一区二区三区| 免费搜索国产男女视频| 自拍偷自拍亚洲精品老妇| 亚洲无线观看免费| 国产亚洲欧美98| 给我免费播放毛片高清在线观看| bbb黄色大片| 亚洲最大成人手机在线| АⅤ资源中文在线天堂| 免费观看人在逋| 18禁裸乳无遮挡免费网站照片| 99视频精品全部免费 在线| 亚洲人成伊人成综合网2020| 在线看三级毛片| 国产成人av教育| x7x7x7水蜜桃| 精品久久久久久久末码| 91精品国产九色| 国产高清三级在线| 亚洲最大成人中文| 亚洲性久久影院| 亚洲中文日韩欧美视频| 日韩国内少妇激情av| 日日啪夜夜撸| 波多野结衣高清作品| 久久热精品热| av中文乱码字幕在线| 99热只有精品国产| 亚洲av第一区精品v没综合| 无人区码免费观看不卡| 校园春色视频在线观看| 精品一区二区三区av网在线观看| 韩国av一区二区三区四区| 中国美白少妇内射xxxbb| 亚洲成人免费电影在线观看| eeuss影院久久| 欧美性猛交╳xxx乱大交人| 亚洲欧美日韩高清专用| 久久久久精品国产欧美久久久| 欧美绝顶高潮抽搐喷水| 日韩一本色道免费dvd| 高清日韩中文字幕在线| 亚洲av第一区精品v没综合| 麻豆国产av国片精品| 国产aⅴ精品一区二区三区波| 欧美人与善性xxx| 全区人妻精品视频| 精品人妻视频免费看| 国产成年人精品一区二区| 免费观看人在逋| 欧美xxxx性猛交bbbb| 免费在线观看成人毛片| 国产在视频线在精品| 欧美+日韩+精品| 国内精品一区二区在线观看| 欧美日韩国产亚洲二区| 在线国产一区二区在线| 日本一本二区三区精品| 久久久久久伊人网av| 丰满的人妻完整版| 亚洲av成人av| avwww免费| 少妇的逼好多水| 欧美精品国产亚洲| 免费高清视频大片| 少妇猛男粗大的猛烈进出视频 | 丝袜美腿在线中文| 一进一出抽搐动态| 国产欧美日韩精品一区二区| 在线a可以看的网站| 亚洲国产高清在线一区二区三| 免费无遮挡裸体视频| 日韩国内少妇激情av| 国产av一区在线观看免费| 最近最新中文字幕大全电影3| 婷婷精品国产亚洲av| 黄片wwwwww| 中文字幕熟女人妻在线| 最近在线观看免费完整版| 我要搜黄色片| 黄色视频,在线免费观看| 18禁黄网站禁片午夜丰满| 别揉我奶头 嗯啊视频| 内地一区二区视频在线| 欧美日韩国产亚洲二区| 婷婷精品国产亚洲av在线| 精品乱码久久久久久99久播| 国产高清有码在线观看视频| bbb黄色大片| 看片在线看免费视频| 国产成人一区二区在线| 精品日产1卡2卡| 久久久国产成人精品二区| 中文字幕人妻熟人妻熟丝袜美| 日韩一区二区视频免费看| 黄色视频,在线免费观看| 一区福利在线观看| 我要搜黄色片| 日本撒尿小便嘘嘘汇集6| 亚洲av中文字字幕乱码综合| 午夜免费男女啪啪视频观看 | 99热这里只有是精品在线观看| 亚洲av中文字字幕乱码综合| 亚洲人成网站在线播放欧美日韩| 久久精品夜夜夜夜夜久久蜜豆| 国产精品免费一区二区三区在线| 亚洲欧美日韩卡通动漫| 性欧美人与动物交配| 国产探花在线观看一区二区| 日韩中字成人| 99久久成人亚洲精品观看| 真实男女啪啪啪动态图| 在线观看66精品国产| 不卡视频在线观看欧美| 尾随美女入室| 国产亚洲精品av在线| 黄色一级大片看看| 在线天堂最新版资源| av在线蜜桃| 久久精品人妻少妇| 日韩欧美一区二区三区在线观看| 国产男人的电影天堂91| 69人妻影院| 免费在线观看影片大全网站| 大型黄色视频在线免费观看| а√天堂www在线а√下载| 国产 一区 欧美 日韩| 欧美又色又爽又黄视频| 国产精品久久视频播放| 中文资源天堂在线| 久久99热这里只有精品18| 真人一进一出gif抽搐免费| 欧美人与善性xxx| 草草在线视频免费看| 精品一区二区三区视频在线| 成人国产综合亚洲| a级毛片免费高清观看在线播放| 亚洲人成网站在线播| 国产极品精品免费视频能看的| 亚洲最大成人av| 精品99又大又爽又粗少妇毛片 | 亚洲国产精品成人综合色| 中文字幕久久专区| 在线观看一区二区三区| 精品无人区乱码1区二区| 国产精品1区2区在线观看.| 久久精品综合一区二区三区| 九九爱精品视频在线观看| 两个人的视频大全免费| 欧美成人a在线观看| 免费一级毛片在线播放高清视频| 草草在线视频免费看| 三级毛片av免费| 国产精品1区2区在线观看.| 国产一区二区亚洲精品在线观看| 久久久久久久久久久丰满 | 午夜老司机福利剧场| 真人做人爱边吃奶动态| 一a级毛片在线观看| 1024手机看黄色片| 婷婷精品国产亚洲av| a级毛片a级免费在线| 亚洲精品乱码久久久v下载方式| 成人毛片a级毛片在线播放| 丝袜美腿在线中文| 亚洲,欧美,日韩| 色吧在线观看| 干丝袜人妻中文字幕| 国产精品久久久久久久电影| 舔av片在线| av在线老鸭窝| 久久精品久久久久久噜噜老黄 | 村上凉子中文字幕在线| 午夜老司机福利剧场| 欧美高清成人免费视频www| 色吧在线观看| 日本三级黄在线观看| 成人国产一区最新在线观看| 午夜福利视频1000在线观看| 少妇高潮的动态图| 少妇裸体淫交视频免费看高清| 天堂网av新在线| 亚洲国产精品sss在线观看| 久9热在线精品视频| 欧美精品啪啪一区二区三区| 18禁在线播放成人免费| 小蜜桃在线观看免费完整版高清| 在线免费观看的www视频| 干丝袜人妻中文字幕| 我的老师免费观看完整版| av视频在线观看入口| 国产在线男女| 美女高潮喷水抽搐中文字幕| 国产亚洲欧美98| 身体一侧抽搐| 男插女下体视频免费在线播放| 不卡视频在线观看欧美| 18禁在线播放成人免费| 久久久久久久久中文| 一进一出好大好爽视频| 女的被弄到高潮叫床怎么办 | 成人鲁丝片一二三区免费| 12—13女人毛片做爰片一| 色播亚洲综合网| 韩国av一区二区三区四区| 精品午夜福利在线看| 午夜福利成人在线免费观看| 午夜影院日韩av| 午夜精品久久久久久毛片777| 国产白丝娇喘喷水9色精品| 我要看日韩黄色一级片| 窝窝影院91人妻| 91麻豆精品激情在线观看国产| 欧美丝袜亚洲另类 | 午夜激情福利司机影院| 老熟妇仑乱视频hdxx| 午夜福利在线观看吧| 亚洲无线观看免费| 噜噜噜噜噜久久久久久91| 亚洲人成网站高清观看| www.色视频.com| 久久99热6这里只有精品| 免费在线观看影片大全网站| 一本一本综合久久| 亚洲人成网站在线播放欧美日韩| 久久午夜福利片| 伦精品一区二区三区| АⅤ资源中文在线天堂| 黄色欧美视频在线观看| 我的女老师完整版在线观看| 国产精品永久免费网站| 日本色播在线视频| 两人在一起打扑克的视频| 日韩欧美 国产精品| 中国美女看黄片| 欧美精品啪啪一区二区三区| 观看美女的网站| 成人一区二区视频在线观看| 欧美精品国产亚洲| 人妻丰满熟妇av一区二区三区| 婷婷色综合大香蕉| 免费看av在线观看网站| 欧美又色又爽又黄视频| 在线免费观看不下载黄p国产 | 露出奶头的视频| 九九爱精品视频在线观看| 亚洲国产精品sss在线观看| 中文字幕精品亚洲无线码一区| 久久精品91蜜桃| 有码 亚洲区| 麻豆成人午夜福利视频| 亚洲人成网站在线播放欧美日韩| 制服丝袜大香蕉在线| 日本三级黄在线观看| 色在线成人网| 成人鲁丝片一二三区免费| 亚洲 国产 在线| 毛片女人毛片| 999久久久精品免费观看国产| 听说在线观看完整版免费高清| 日韩一区二区视频免费看| 欧美黑人欧美精品刺激| 亚洲在线自拍视频| 亚洲av免费在线观看| 亚洲精品乱码久久久v下载方式| 欧美中文日本在线观看视频| 亚洲欧美日韩东京热| 91麻豆精品激情在线观看国产| 国产淫片久久久久久久久| 久久国产精品人妻蜜桃| 日韩中文字幕欧美一区二区| 国产精品久久久久久精品电影| 国产欧美日韩精品亚洲av| 久久精品国产亚洲网站| 看免费成人av毛片| 亚洲精品色激情综合| 99热6这里只有精品| 99久久九九国产精品国产免费| 88av欧美| 欧美色欧美亚洲另类二区| 国产精品98久久久久久宅男小说| 国产亚洲91精品色在线| 一级av片app| 亚洲电影在线观看av| 噜噜噜噜噜久久久久久91| 免费黄网站久久成人精品| 精品一区二区三区人妻视频| 一夜夜www| 国产色爽女视频免费观看| 又粗又爽又猛毛片免费看| 在线观看午夜福利视频| 精品久久久久久久人妻蜜臀av| 久久精品夜夜夜夜夜久久蜜豆| 国产伦精品一区二区三区视频9| 日韩中文字幕欧美一区二区| 精品午夜福利在线看| 熟女人妻精品中文字幕| 国产极品精品免费视频能看的| 日本一二三区视频观看| 久久香蕉精品热| 国产精品久久视频播放| 国产精品久久久久久久电影| 日韩亚洲欧美综合| 亚洲中文字幕一区二区三区有码在线看| 成人无遮挡网站| 精品福利观看| 午夜福利在线观看吧| 免费av不卡在线播放| 日韩av在线大香蕉| 网址你懂的国产日韩在线| 一个人观看的视频www高清免费观看| 日韩,欧美,国产一区二区三区 | 99热精品在线国产| 丰满乱子伦码专区| 91狼人影院| 我要搜黄色片| 婷婷亚洲欧美| 国产精品久久视频播放| 亚洲电影在线观看av| 国产亚洲精品久久久com| 国产熟女欧美一区二区| 日韩高清综合在线| 中文亚洲av片在线观看爽| 男插女下体视频免费在线播放| 亚洲国产高清在线一区二区三| 99久久精品国产国产毛片| 国产极品精品免费视频能看的| 午夜精品在线福利| ponron亚洲| 国产美女午夜福利| 人人妻人人澡欧美一区二区| 99久国产av精品| 色噜噜av男人的天堂激情| 亚洲不卡免费看| 男女做爰动态图高潮gif福利片| 免费黄网站久久成人精品| 国产高清激情床上av| 欧美另类亚洲清纯唯美| 国产精品av视频在线免费观看| 亚洲精品456在线播放app | 中文在线观看免费www的网站| 久久久成人免费电影| 乱码一卡2卡4卡精品| 丰满人妻一区二区三区视频av| 国产综合懂色| 日本精品一区二区三区蜜桃| 国产精品电影一区二区三区| 欧美性猛交╳xxx乱大交人| 内射极品少妇av片p| 亚洲av熟女| 日韩精品中文字幕看吧| 亚洲精品久久国产高清桃花| 九九在线视频观看精品| 精品欧美国产一区二区三| 一区二区三区高清视频在线| 国产中年淑女户外野战色| 村上凉子中文字幕在线| 国产精品久久电影中文字幕| 国产91精品成人一区二区三区| 天堂影院成人在线观看| 99在线人妻在线中文字幕| 18+在线观看网站| 嫩草影院精品99| 欧美日韩国产亚洲二区| 久久亚洲真实| 久久香蕉精品热| 欧美色视频一区免费| 国产女主播在线喷水免费视频网站 | 国产精品日韩av在线免费观看| 精品久久久久久久人妻蜜臀av| 嫩草影院精品99| 亚洲欧美日韩东京热| 欧美黑人巨大hd| 久久香蕉精品热| 午夜精品久久久久久毛片777| 亚洲精品亚洲一区二区| 国产视频内射| 国内揄拍国产精品人妻在线| 天堂影院成人在线观看| 国产一区二区在线av高清观看| 在线播放无遮挡| 欧美成人性av电影在线观看| 欧美最新免费一区二区三区| 久久久国产成人免费| 成人国产一区最新在线观看| 最新在线观看一区二区三区| 搡女人真爽免费视频火全软件 | 99精品久久久久人妻精品| 窝窝影院91人妻| 俄罗斯特黄特色一大片| 黄色日韩在线| 日韩av在线大香蕉| av天堂中文字幕网| 亚洲av成人av| 精品久久久久久成人av| 亚洲国产日韩欧美精品在线观看| 不卡一级毛片| 99国产精品一区二区蜜桃av| 国产伦精品一区二区三区四那| 色综合亚洲欧美另类图片| 亚洲一区高清亚洲精品| 成人国产一区最新在线观看| 美女cb高潮喷水在线观看| 亚洲va日本ⅴa欧美va伊人久久| 国产国拍精品亚洲av在线观看| 别揉我奶头 嗯啊视频| 国产一区二区三区av在线 | 非洲黑人性xxxx精品又粗又长| 国产亚洲精品久久久久久毛片| 深爱激情五月婷婷| 欧美一级a爱片免费观看看| 女生性感内裤真人,穿戴方法视频| 国产在线精品亚洲第一网站| 久久99热这里只有精品18| 中文字幕熟女人妻在线| 久久99热6这里只有精品| 国内精品一区二区在线观看| 亚洲国产高清在线一区二区三| 中国美女看黄片| 又黄又爽又免费观看的视频| 日本a在线网址| 日韩国内少妇激情av| 嫩草影院新地址| 天天一区二区日本电影三级| 九九热线精品视视频播放| 窝窝影院91人妻| 国产av不卡久久| 久久草成人影院| 别揉我奶头 嗯啊视频| 日韩中字成人| 中出人妻视频一区二区| 亚洲第一区二区三区不卡| 精品午夜福利在线看| 久久久久久久久久久丰满 | 国内毛片毛片毛片毛片毛片| 国产精品乱码一区二三区的特点| 欧美性猛交╳xxx乱大交人| 日韩欧美在线乱码| 99九九线精品视频在线观看视频| 精品久久久久久久人妻蜜臀av| 久久久久久伊人网av| 中文字幕精品亚洲无线码一区| 国产极品精品免费视频能看的| 97热精品久久久久久| 又黄又爽又免费观看的视频| 嫩草影院精品99| 精品一区二区三区av网在线观看| 日韩欧美精品v在线| 欧美色视频一区免费| 国产高清不卡午夜福利| 校园人妻丝袜中文字幕| 国产av一区在线观看免费| 午夜福利在线在线| 午夜爱爱视频在线播放| 熟女电影av网| 亚洲av日韩精品久久久久久密| 黄色欧美视频在线观看| 日韩精品有码人妻一区| 韩国av在线不卡| 夜夜看夜夜爽夜夜摸| 欧美一区二区亚洲| 十八禁国产超污无遮挡网站| 色噜噜av男人的天堂激情| 成年女人毛片免费观看观看9| av在线蜜桃| 性插视频无遮挡在线免费观看| 久久久久久久久中文| 人人妻,人人澡人人爽秒播| 欧美一区二区精品小视频在线| 亚洲av一区综合| 免费人成在线观看视频色| 久久久久久久久大av| 亚洲精品456在线播放app | 草草在线视频免费看| 国产精品久久久久久久久免| 嫩草影视91久久| 美女cb高潮喷水在线观看| 又黄又爽又免费观看的视频| 九九热线精品视视频播放| netflix在线观看网站| 长腿黑丝高跟| 精品久久久久久久久久免费视频| 黄色女人牲交| 亚洲国产精品久久男人天堂| 精品人妻视频免费看| 久久99热6这里只有精品| 日韩欧美国产一区二区入口| 日本欧美国产在线视频| 亚洲av五月六月丁香网| 在线观看美女被高潮喷水网站| 亚洲天堂国产精品一区在线| 男女啪啪激烈高潮av片| 国产亚洲av嫩草精品影院| 国产免费男女视频| 婷婷精品国产亚洲av在线| 久久久久久久精品吃奶| 99热精品在线国产| 久久久久久伊人网av| 亚洲人成网站在线播| 99久久久亚洲精品蜜臀av| 久久天躁狠狠躁夜夜2o2o| 中文资源天堂在线| АⅤ资源中文在线天堂| 欧美高清性xxxxhd video| 午夜福利18| 国产亚洲精品久久久com| 黄色一级大片看看| 天天躁日日操中文字幕| 伦理电影大哥的女人| 一个人免费在线观看电影| 欧美绝顶高潮抽搐喷水| 超碰av人人做人人爽久久| 桃色一区二区三区在线观看| 日韩一区二区视频免费看| 天堂√8在线中文| 亚洲人成网站高清观看| 美女高潮喷水抽搐中文字幕| av在线老鸭窝| 精品乱码久久久久久99久播| av在线老鸭窝| 女生性感内裤真人,穿戴方法视频| 午夜视频国产福利| 嫩草影院精品99| 啦啦啦韩国在线观看视频| 精品午夜福利在线看| 免费在线观看日本一区| 精品久久国产蜜桃| 在线观看舔阴道视频| 欧美日韩乱码在线| 国产精品乱码一区二三区的特点| 亚洲欧美日韩卡通动漫| 精品福利观看| videossex国产| 直男gayav资源| 成人精品一区二区免费| 久久久成人免费电影| 亚洲av电影不卡..在线观看| 男人舔女人下体高潮全视频| 亚洲人成伊人成综合网2020| 村上凉子中文字幕在线| 精品人妻熟女av久视频| 无遮挡黄片免费观看| 内地一区二区视频在线| 久久中文看片网|