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

    基于卡爾曼濾波算法的稻縱卷葉螟短期預(yù)測模型*

    2016-10-13 06:00:06包云軒陳心怡謝曉金陸明紅
    中國農(nóng)業(yè)氣象 2016年5期
    關(guān)鍵詞:卡爾曼發(fā)生量卷葉螟

    包云軒,陳心怡,謝曉金,王 琳,陸明紅

    ?

    基于卡爾曼濾波算法的稻縱卷葉螟短期預(yù)測模型*

    包云軒1,2,陳心怡1,2,謝曉金1,2,王 琳1,2,陸明紅3

    (1.南京信息工程大學(xué)氣象災(zāi)害預(yù)報和評估協(xié)同創(chuàng)新中心,南京 210044;2.江蘇省農(nóng)業(yè)氣象重點實驗室/南京信息工程大學(xué),南京 210044;3.農(nóng)業(yè)部全國農(nóng)業(yè)技術(shù)推廣與服務(wù)中心,北京 100125)

    利用1994-2014年中國南方四大稻區(qū)(華南、西南、江嶺和江淮稻區(qū))代表性病蟲測報站的稻縱卷葉螟逐候田間趕蛾量資料,篩選出影響各站稻縱卷葉螟發(fā)生量的關(guān)鍵氣象因子,應(yīng)用卡爾曼濾波方法分別對各站建立稻縱卷葉螟遷入期候發(fā)生量的卡爾曼短期預(yù)測模型,并計算模型的準(zhǔn)確率、誤差大小和穩(wěn)定性。結(jié)果表明:(1)稻縱卷葉螟發(fā)生量與前一候和前兩候的田間蛾量呈極顯著正相關(guān)(P<0.01),與前一候的近地面最低氣溫、平均氣溫和最高氣溫呈極顯著正相關(guān)(P<0.01),與前一候的地面氣壓呈極顯著負(fù)相關(guān)(P<0.01)。(2)經(jīng)1994-2011年的回檢擬合和2012-2014年試報檢驗,卡爾曼模型的發(fā)生量預(yù)測綜合平均誤差為-88.63,平均絕對誤差為217.72,均方根誤差為605.04。發(fā)生量預(yù)測綜合準(zhǔn)確率為84.33%,平均歷史擬合率為83.33%,各站卡爾曼模型的預(yù)報結(jié)果與實測值基本吻合,表明模型可以應(yīng)用于稻縱卷葉螟候發(fā)生量的預(yù)測。

    稻縱卷葉螟;氣象因子;卡爾曼濾波算法;候發(fā)生量預(yù)報模型;準(zhǔn)確率

    稻縱卷葉螟(Guenee)是為害水稻的一種遷飛性害蟲,每年3-10月在中國有5次大規(guī)模的南北往返遷飛[1]。20世紀(jì)60年代中期以前在中國為局部偶發(fā)性害蟲,70-80年代開始其發(fā)生頻次和危害程度大幅升高,逐漸成為中國主要稻區(qū)的普遍性和常發(fā)性害蟲[2]。21世紀(jì)以來,稻縱卷葉螟在中國各稻區(qū)的發(fā)生量和危害面積明顯上升,2003年以來在全國連年大發(fā)生,給水稻生產(chǎn)造成了嚴(yán)重?fù)p失[3-4]。由于近10a來,稻縱卷葉螟的遷入與發(fā)生趨于復(fù)雜化、持續(xù)化,而其發(fā)生又具有間歇性和突發(fā)性的特點,其年際發(fā)生量、發(fā)生程度和發(fā)生區(qū)域差異均較大,分布極不均勻[5],更為稻縱卷葉螟的預(yù)測預(yù)報增加了難度。因此,研究稻縱卷葉螟的預(yù)測預(yù)報方法,對減少糧食損失,維護(hù)社會穩(wěn)定具有重要意義。

    稻縱卷葉螟的發(fā)生除與蟲源地、水稻種植制度、地形分布、天敵情況等因素相關(guān)外,還與氣象條件有著極為密切關(guān)系,研究表明[6-10],稻縱卷葉螟的遷入發(fā)生量受高空氣流、溫度、濕度、降水和光照等氣象條件的綜合影響。

    關(guān)于稻縱卷葉螟發(fā)生的預(yù)報方法已有很多,最早使用的是經(jīng)驗預(yù)測法[11],主要根據(jù)害蟲的生長發(fā)育狀況和天氣形勢對稻縱卷葉螟的主害代發(fā)生峰期和發(fā)生程度進(jìn)行預(yù)測。經(jīng)驗預(yù)測法簡單直觀,但是主觀因素帶來的誤差明顯,難以達(dá)到較高的準(zhǔn)確率。之后數(shù)理統(tǒng)計方法得到廣泛應(yīng)用,所用方法包括回歸分析、多維時間序列分析、模糊貼近度、灰色系統(tǒng)、周期分析等[12]。這類方法主要通過因子之間的簡單關(guān)系分析來建模,機(jī)理不明,隨機(jī)性“噪聲”導(dǎo)致的誤差也較明顯。近年來隨著信息技術(shù)的不斷發(fā)展,人工神經(jīng)網(wǎng)絡(luò)[13]等方法也被廣泛應(yīng)用,但神經(jīng)網(wǎng)絡(luò)建模的缺點是其建模存在“黑箱”性,網(wǎng)絡(luò)的學(xué)習(xí)和決策過程難以理解,不利于對模型的進(jìn)一步優(yōu)化。

    目前對病蟲害的預(yù)測多用回歸方法預(yù)測世代或年發(fā)生程度和發(fā)生高峰期[14-17],對實際發(fā)生量的預(yù)測不多,以候為單位的短期預(yù)測很少見報道。目前對病蟲害的預(yù)測仍以回歸方法為主,其預(yù)測參數(shù)不能隨實際情況進(jìn)行更新,而本文應(yīng)用的卡爾曼濾波算法是一個遞歸算法,能夠從一系列不完全和包含噪聲的量測中,估計動態(tài)系統(tǒng)的狀態(tài)[18]。該方法于1960年被提出,之后被廣泛應(yīng)用于航空、航天、工業(yè)過程、社會經(jīng)濟(jì)和天氣預(yù)報等不同領(lǐng)域[19-22],但在病蟲害測報領(lǐng)域的應(yīng)用報道極為罕見。2010年江勝國等[23]首次應(yīng)用卡爾曼濾波方法建立了稻飛虱發(fā)生適宜氣象等級預(yù)測模型,而利用卡爾曼濾波方法對稻縱卷葉螟發(fā)生量進(jìn)行預(yù)測的研究或應(yīng)用還未見報道。因此,本文嘗試?yán)每柭鼮V波算法與前期對稻縱卷葉螟發(fā)生量有顯著影響的氣象因子相結(jié)合,建立其候發(fā)生量的短期預(yù)報模型,旨在為中國稻縱卷葉螟的預(yù)測預(yù)報提供一種思路。

    1 資料與方法

    1.1 資料來源

    1.1.1 蟲情資料

    考慮到資料的完整性和代表性,從華南稻區(qū)、西南稻區(qū)、江嶺稻區(qū)和江淮稻區(qū)分別選取廣西全州、重慶秀山、湖南湘陰和江蘇張家港4個植物保護(hù)站作為代表站,利用其1994-2014年逐日田間趕蛾量數(shù)據(jù)。以候為時間單位,利用MatLab編程將各站逐日田間趕蛾量資料處理成候平均數(shù)據(jù),并進(jìn)行對數(shù)化處理。

    1.1.2 氣象資料

    選取1994-2014年全州、秀山、湘陰和張家港4個氣象站逐日常規(guī)氣象觀測資料,氣象要素為地面氣壓、風(fēng)速、最低氣溫、平均氣溫、最高氣溫、相對濕度、降水量、日照時數(shù)等,所有數(shù)據(jù)均進(jìn)行候平均處理。

    1.2 方法

    1.2.1 逐步回歸方法

    分別對各站稻縱卷葉螟逐候發(fā)生量與前期發(fā)生量及氣象因子之間的相關(guān)性進(jìn)行分析,選取通過P<0.01水平顯著性檢驗的因子,再利用逐步回歸分析方法篩選4個代表站稻縱卷葉螟發(fā)生量預(yù)測的預(yù)報因子。

    1.2.2 卡爾曼濾波模型

    卡爾曼濾波算法的基本思想是,以最小均方誤差為最佳估計準(zhǔn)則,采用信號與噪聲的狀態(tài)空間模型,利用前一時刻預(yù)測誤差來更新對狀態(tài)變量的估計,求出當(dāng)前時刻的估計值,根據(jù)建立的狀態(tài)方程和量測方程對需要處理的信號作出滿足最小均方誤差的估計[24]。

    本研究在氣象上常用的卡爾曼濾波系統(tǒng)方程基礎(chǔ)上加以改進(jìn),得到適用于病蟲害預(yù)測的卡爾曼濾波預(yù)測方程,即

    預(yù)測遞推方程為

    式(3)中的a是變動的,預(yù)測t+1候的at由式(8)得到,而式(8)又由式(5)-式(7)得到,需要用到Ct-1、St、At和Rt,預(yù)測t+2候的at時需要先計算Ct,可通過式(9)實現(xiàn),依次類推,以實現(xiàn)模型參數(shù)更新。

    1.2.3 預(yù)測效果評價指標(biāo)

    選取平均誤差(Mean Error,ME)、平均絕對誤差(Mean Abosolute Error,MAE)、均方根誤差(Root Mean Square Error,RMSE)和預(yù)測準(zhǔn)確率(Prediction Accurracy,PA)等指標(biāo)來評價預(yù)測模型的誤差大小、穩(wěn)定性和預(yù)測準(zhǔn)確率。

    式中,yi為預(yù)測候蟲量值,單位為頭·667m-2,xi為實測候蟲量值,單位為頭·667m-2。i=1,2,3,…,n,i代表遷入期內(nèi)候序,n為各站年遷入期總候數(shù)。ME反映綜合預(yù)測誤差,ME越大預(yù)測誤差越大。與此相比,MAE由于離差被絕對值化,不會出現(xiàn)正負(fù)相抵消的情況,因而,能更好地反映預(yù)測值誤差的實際情況。MAE越大表示預(yù)測誤差越大。RMSE越大,表示預(yù)測值相對于實測值的離散程度越大,RMSE越小,模型的穩(wěn)定性越好。

    用唐啟義等[25]提出的病蟲測報應(yīng)驗程度判定模式計算預(yù)測準(zhǔn)確性評價的分值,即

    式中,SD為發(fā)生量判定模式的分值;SD<40表明預(yù)報不準(zhǔn)確,40≤SD<60表明預(yù)報較準(zhǔn)確,SD≥60表明預(yù)報準(zhǔn)確;a為預(yù)報對象常年平均值,a1為實測值,a2為預(yù)測值;δ為預(yù)報對象常年標(biāo)準(zhǔn)差,t為自預(yù)報發(fā)出至實際發(fā)生時的期距(d),因為本研究是以候為時間單位預(yù)報,故d=5。

    根據(jù)分值計算預(yù)測準(zhǔn)確率(PA)和歷史擬合率(HA)

    式中,Ap為當(dāng)年預(yù)測評分準(zhǔn)確的候數(shù);Sp為當(dāng)年進(jìn)行預(yù)測的總候數(shù)。AY為歷史擬合準(zhǔn)確的年數(shù),SY為歷史擬合的總年數(shù)。PA越大,表明預(yù)測準(zhǔn)確率越高。HA越大,表明歷史擬合率越高。

    1.3 遷入期和預(yù)報時間

    根據(jù)1994-2014年各站歷年稻縱卷葉螟逐日蟲量資料,最早始見期為該站21a中稻縱卷葉螟每年第一次遷入最早的候序,最晚終見期為該站21a中稻縱卷葉螟每年最后一次遷入最晚的候序。遷入期即為最早始見期至最晚終見期之間的候數(shù)。4個代表站的遷入期見表1。本文4個站每年預(yù)報的時間分別為該站的遷入期起始候到終見候之間的候數(shù),提前1候預(yù)報下一候的發(fā)生量。

    表1 南方四大稻區(qū)代表站遷入期的起止候

    2 結(jié)果與分析

    2.1 稻縱卷葉螟發(fā)生量預(yù)報模型中因子的確定

    2.1.1 發(fā)生量影響因子

    全州、秀山、湘陰和張家港4個站選用建模的氣象因子有氣壓、風(fēng)速、最低氣溫、平均氣溫、最高氣溫、相對濕度、降水量和日照時數(shù)。由于前期蟲源基數(shù)對后期稻縱卷葉螟的發(fā)生有直接影響,所以加上前一候發(fā)生量、前兩候發(fā)生量這兩個影響因子,分別與各站稻縱卷葉螟發(fā)生量進(jìn)行相關(guān)性分析。結(jié)果見表2。

    表2 四站當(dāng)前候稻縱卷葉螟發(fā)生量與前期發(fā)生量、前期氣象因子的相關(guān)系數(shù)

    注:A1為前一候的蟲口數(shù)量(頭·667m-2),A2為前兩候的蟲口數(shù)量(頭·667m-2),P為前一候氣壓(hPa),TA為前一候地面平均溫度(℃),TX為前一候地面最高氣溫(℃),TN為前一候地面最低氣溫(℃),V為前一候風(fēng)速(m·s-1),RH為前一候相對濕度(%),R為前一候降水量(mm),SH為前一候日照時數(shù)(h)。**表示P<0.01,*表示P<0.05,?表示相關(guān)性不顯著。下同。

    Note:A1is thes moth amounts at the preceding pentad(Nr·667m-2), A2is the moth amounts at the two preceding pentad(Nr·667m-2),P is pressure(hPa) ,TAis mean temperature (℃),TXis average maximum temperature(℃), TNis average minimum temperature(℃), V is wind speed(m·s-1), RH is relative humidity(%), R is precipitation(mm), SH is sunshine duration(h), all these factors are in the preceding pentad .**is P<0.01,*is P<0.05, ? is P>0.05. The same as below.

    由表2可見,(1)4站稻縱卷葉螟發(fā)生量與前一候、前兩候的發(fā)生量的相關(guān)性均表現(xiàn)為極顯著,表明前期田間蟲源基數(shù)越大,越有利于后期稻縱卷葉螟的發(fā)生。(2)4站稻縱卷葉螟發(fā)生量與前一候平均最低氣溫、平均氣溫和平均最高氣溫呈極顯著正相關(guān),表明溫度越高,越有利于稻縱卷葉螟的發(fā)生,適溫范圍內(nèi)的較高溫的環(huán)境有利于稻縱卷葉螟的遷飛和發(fā)育,這與其喜溫習(xí)性相吻合。(3)4站稻縱卷葉螟發(fā)生量與前一候氣壓呈極顯著負(fù)相關(guān),實際上在稻縱卷葉螟發(fā)生季內(nèi),氣壓越低,天氣越悶熱潮濕,越利于喜溫喜濕的稻縱卷葉螟的發(fā)生和危害。(4)全州站稻縱卷葉螟發(fā)生量與前一候風(fēng)速呈顯著負(fù)相關(guān),表明風(fēng)速越大,越不利于該站稻縱卷葉螟的發(fā)生。因為該站地處西南偏南至東北偏北走向的湘桂走廊內(nèi),東西兩側(cè)均為較高的地形,風(fēng)速過大,近地層擾動過大,不利于稻縱卷葉螟在水稻田間駐留。(5)全州站稻縱卷葉螟發(fā)生量與前一候日照時數(shù)呈顯著正相關(guān),表明該站前一候日照時數(shù)越多,稻縱卷葉螟發(fā)生量越大。全州站常年春寒時間長,若日照時數(shù)多,則太陽輻射強(qiáng),有利于氣溫升高,能促進(jìn)水稻生長,對喜暖的稻縱卷葉螟的遷入和為害有利。(6)湘陰站稻縱卷葉螟發(fā)生量與前一候的相對濕度呈顯著正相關(guān),表明相對濕度越大,越有利于湘陰站稻縱卷葉螟的發(fā)生。湘陰站地處湖南省東北部的濱湖平原,境內(nèi)水系發(fā)達(dá),水陸環(huán)流旺盛,高濕環(huán)境既有利于降蟲,又利于田間取食、為害和繁殖,這與其喜濕習(xí)性十分吻合。(7)秀山站稻縱卷葉螟發(fā)生量與前一候降水呈顯著正相關(guān)。事實上,降水對空中遷飛的稻縱卷葉螟種群的降落有明顯的動力脅迫作用,一定強(qiáng)度以下的降水越多,越有利于稻縱卷葉螟的遷入,也越利于其田間的發(fā)生。秀山系武陵山二級隆起帶的南段,境內(nèi)地表起伏大,西南高,東北低,地形阻礙不利于西南氣流攜帶蟲群遷入,種群遷入主要靠降水動力脅迫來實現(xiàn),而秀山站充沛的降水,對稻縱卷葉螟的遷入為害十分有利。從表2還可看出,由于各站地理位置和地形不同,氣候特點不同,水稻生長條件和狀況也不一樣,不同的氣象因子組合呈現(xiàn)出對稻縱卷葉螟發(fā)生量影響的不同結(jié)果。

    2.1.2 預(yù)報因子選取

    選取上述與發(fā)生量相關(guān)性通過P<0.01顯著性檢驗的因子進(jìn)行逐步回歸,因子進(jìn)入方程條件為F=0.15,刪除因子條件為F=0.2,最后剔除顯著性檢驗P≥0.05的因子。各站入選的預(yù)報因子見表3。

    表3 逐步回歸法入選的各站稻縱卷葉螟發(fā)生量預(yù)測因子

    4站地理位置由南往北排序為全州(25.56°N,111.04°E)、秀山(28.27°N,108.59°E)、湘陰(28.41°N,112.53°E)、張家港(31.88°N,120.64°E),秀山站南鄰貴州寡照區(qū),因日照少和降水多的影響,全年氣溫不如全州和湘陰高,作為稻縱卷葉螟發(fā)生熱量條件的平均氣溫因子顯得較為重要;張家港位置偏北,年平均氣溫最低,所以最高氣溫和最低氣溫對其遷入田間為害較關(guān)鍵;全州和湘陰站年平均氣溫相對較高,但春、秋季的低溫會明顯影響稻縱卷葉螟的遷入和為害,因而最低氣溫是關(guān)鍵因子;全州站春寒時間長,如果日照時數(shù)多,太陽輻射強(qiáng),將有利于氣溫升高,能促進(jìn)水稻生長,也對喜暖的稻縱卷葉螟的遷入和為害有利,因此,日照時數(shù)是關(guān)鍵因子之一;秀山站和湘陰站都處在西南高、東北低的地形過渡區(qū),稻縱卷葉螟的遷入和為害受地形重力波引起的氣壓變化影響較大;湘陰站境內(nèi)河湖交匯,水汽條件較好,對喜濕的稻縱卷葉螟的遷入和為害十分有利。由此可見,由逐步回歸篩選的預(yù)報因子能較好地解釋限制性氣象條件對中國不同稻區(qū)稻縱卷葉螟發(fā)生量影響的機(jī)制。

    2.2 稻縱卷葉螟發(fā)生量預(yù)報模型

    將表3中確定的預(yù)報因子作為卡爾曼模型的預(yù)報因子,以逐步回歸方法確定的回歸系數(shù)作為各站卡爾曼模型的初始回歸系數(shù),分別以4個站的1994-2010年和1995-2011年候數(shù)據(jù)為樣本建立兩個量測方程,從而確定模型其它初始參數(shù)。用建立的各站模型分別預(yù)測4站2012-2014年的候發(fā)生量,同時以各站2012-2014年的實測數(shù)據(jù)檢驗預(yù)測效果,并對各站1994-2011年的數(shù)據(jù)進(jìn)行歷史回檢。本文僅列出秀山站的卡爾曼預(yù)測模型初始參數(shù)。

    初始回歸系數(shù)為通過逐步回歸方法對1995-2011年資料建立的回歸方程的回歸系數(shù),即

    用逐步回歸方法建立兩個量測方程中的回歸系數(shù)的變化對系統(tǒng)噪聲方差陣W進(jìn)行估算。

    量測噪聲的方差矩陣V以1994-2010年建立的回歸方程中得到的殘差確定,V=195.321。

    C0是a0的誤差方差陣,默認(rèn)

    轉(zhuǎn)移系數(shù)根據(jù)歷史擬合值與實際值的差距調(diào)整,本文d=0.1。

    確定了模型的初始參數(shù),按照式(3)-式(9),可得到每一候的預(yù)報量和at、Ct。at即為式(3)中預(yù)測每一候發(fā)生量的回歸系數(shù),Ct是更新回歸系數(shù)at時需要用到的誤差方差陣。

    2.3 稻縱卷葉螟發(fā)生量預(yù)報模型的檢驗

    2.3.1 歷史回代檢驗

    應(yīng)用卡爾曼模型對4個站1994-2011年18a遷入期的候數(shù)據(jù)進(jìn)行歷史回檢,結(jié)果見表4。總體看來,卡爾曼模型的4個站發(fā)生量歷史擬合率均在72%以上,平均歷史擬合率為83.33%。

    表4 四站1994-2011年發(fā)生量的歷史擬合率(HA,%)

    2.3.2 預(yù)測效果評價

    利用所建模型對全州站、秀山站、湘陰站和張家港站2012-2014年遷入期的候發(fā)生量進(jìn)行預(yù)測,與實際發(fā)生值的對比結(jié)果見圖1,相應(yīng)的評價指標(biāo)見表5。

    由圖1可見,卡爾曼模型預(yù)測的各站稻縱卷葉螟候發(fā)生量與實測值趨勢基本吻合,且能較好地反映低值區(qū)和高值區(qū)的變化,預(yù)測效果較好。表5顯示,全州站的ME、MAE和RMSE均最小,表明其預(yù)測誤差最小,模型穩(wěn)定性最好。事實上,全州站歷年發(fā)生情況也較平穩(wěn)。張家港站前兩年預(yù)測誤差偏大,且RMSE偏高,可能是因為2012年張家港站是異常大發(fā)生年,而2011年和2013年均為輕發(fā)生年,變化差異過大。2012-2014年卡爾曼模型4個站3a綜合平均ME為-88.63,綜合平均MAE為217.72,綜合平均RMSE為605.04。

    表5 卡爾曼模型預(yù)測的各站2012-2014年稻縱卷葉螟發(fā)生量評價指標(biāo)

    注:ME為平均誤差,MAE為平均絕對誤差,RMSE為均方根誤差,PA為預(yù)測準(zhǔn)確率。

    Note: ME is mean error, MAE is mean absolute error, RMSE is root mean square error, PA is predicting accuracy rate.

    由表5還可見,利用卡爾曼模型預(yù)測的2012-2014年4個站的稻縱卷葉螟發(fā)生量,以2012年張家港站年平均預(yù)測準(zhǔn)確率最低,為61.90%,而該站2013年的預(yù)測準(zhǔn)確率達(dá)到100%的極高準(zhǔn)確度;全州、秀山、湘陰和張家港站卡爾曼模型的預(yù)測準(zhǔn)確率3a平均值分別為82.93%、88.46%、88.17%、77.78%,4站3a平均預(yù)測準(zhǔn)確率為84.33%,表明卡爾曼模型的預(yù)測準(zhǔn)確率較理想。

    3 結(jié)論與討論

    (1)稻縱卷葉螟的發(fā)生量與前一候、前兩候發(fā)生量呈極顯著正相關(guān);與前一候平均氣溫、平均最低氣溫和平均最高氣溫呈極顯著正相關(guān);與前一候地面氣壓均呈極顯著負(fù)相關(guān)。

    (2)本文首次將卡爾曼濾波算法應(yīng)用于稻縱卷葉螟的候發(fā)生量預(yù)測,該模型可提前5d(1候)進(jìn)行預(yù)測,結(jié)果表明卡爾曼模型的綜合預(yù)測準(zhǔn)確率為84.33%,優(yōu)于汪泗水等[13]建立的準(zhǔn)確率為75%的稻縱卷葉螟發(fā)生程度神經(jīng)網(wǎng)絡(luò)預(yù)警模型。

    (3)用卡爾曼濾波算法預(yù)測最重要的是測量方程和狀態(tài)方程的建立,而預(yù)報因子的選取對卡爾曼模型的預(yù)測效果至關(guān)重要。本文雖然結(jié)合地面氣象資料和前期蟲情建立了預(yù)測模型,但是影響稻縱卷葉螟發(fā)生的因素還有很多,如水稻生長狀況、栽培措施、高空氣流、海溫等。下一步研究將綜合考慮其它顯著影響稻縱卷葉螟發(fā)生量的預(yù)報因子加入模型,以提高模型的預(yù)測準(zhǔn)確率。

    (4)卡爾曼模型對于非平穩(wěn)時間序列的預(yù)測效果不如平穩(wěn)時間序列,特別是異常大發(fā)生年的預(yù)測效果不夠理想,如何提高卡爾曼模型對突變點的預(yù)測準(zhǔn)確率,結(jié)合長期預(yù)測預(yù)先設(shè)置下年的預(yù)測權(quán)重是今后將要進(jìn)一步研究的問題。

    參考文獻(xiàn):References

    [1]王鳳英,胡高,陳曉,等.近年來廣西南寧稻縱卷葉螟大發(fā)生原因分析[J].中國水稻科學(xué),2009,23(5):537-545.

    Wang F Y,Hu G,Chen X,et al.Analysis on the causes of Recent outbreaks ofin Nanning, China[J].Chinese Journal of Rice Science, 2009, 23(5): 537-545.(in Chinese)

    [2]全國稻縱卷葉螟聯(lián)合測報網(wǎng).我國稻縱卷葉螟的發(fā)生動態(tài)及其預(yù)測預(yù)報[[J].中國農(nóng)業(yè)科學(xué),1981,24(5):17- 24.

    National Coordinated Network of Forecast for Rice Leaf Roller.The occurrence dynamic and forcast of national rice leaf roller[J].Agricultural Sciences in China,1981, 24(5): 17-24.(in Chinese)

    [3]齊國君,梁載林,辛德育,等.廣西永??h近年來早稻田稻縱卷葉螟連年大發(fā)生的原因分析[J].環(huán)境昆蟲學(xué)報,2009, 31(2):95-101.

    Qi G J,Liang Z L,Xin D Y,et al.Analysis of the successive outbreaks of(Guenee) in early rice paddies in Yongfu county,Guangxi[J].Journal of Environmental Entomology,2009,31(2):95-101.(in Chinese)

    [4]劉宇,王建強(qiáng),馮曉東,等.2007年全國稻縱卷葉螟發(fā)生實況分析與2008年發(fā)生趨勢預(yù)測[J].中國植保導(dǎo)刊,2008, 28(7):33-35.

    Liu Y,Wang J Q,Feng X D,et al.Analysis of national rice leaf roller’s occurrence in 2007 and prediction of occurrence trend in 2008[J].China Plant Protection,2008,28(7):33-35.(in Chinese)

    [5]石尚柏,李紹石,胡伯海.稻縱卷葉螟[M].北京:中國農(nóng)業(yè)出版社,1997:13-17.

    Shi S B,Li S S,Hu B H.[M]. Beijing:China Agricultural Press,1997:13-17.(in Chinese)

    [6]包云軒,曹云,謝曉金,等.中國稻縱卷葉螟發(fā)生特點及北遷的大氣背景[J].生態(tài)學(xué)報,2015,35(11):3519-3533.

    Bao Y X,Cao Y,Xie X J,et al.Migration pattern of rice leaf roller and impact of atmospheric conditions on a heavy migration event in China[J].Acta Ecologica Sinica,2015, 35(11):3519-3533.(in Chinese)

    [7]方源松,廖懷建,錢秋,等.溫濕度對稻縱卷葉螟卵的聯(lián)合作用[J].昆蟲學(xué)報,2013,56(7):786-791.

    Fang Y S,Liao H J,Qian Q,et al.Combined effects of temperature and relative humidity on eggs of the rice leaf roller,(Lepidoptera: Pyralidae)[J]. Acta Entomologica Sinica,2013,56(7):786-791.(in Chinese)

    [8]翟保平,張孝羲.遷飛過程中昆蟲的行為:對風(fēng)溫場的適應(yīng)與選擇[J].生態(tài)學(xué)報,1993,13(4):356-363.

    Zhai B P,Zhang X X.Insects behavior during migration: adaptation and selection of air temperature field[J].Acta Ecologica Sinica,1993,13(4):356-363.(in Chinese)

    [9]吳進(jìn)才.光照、溫度及食物的變化對稻縱卷葉螟遷飛的效應(yīng)[J].昆蟲學(xué)報,1985,28(4):398-405.

    Wu J C.Influence of light,temperature and food change on rice leaf roller migration[J].Acta Entomologica Sinica,1985, 28(4):398-405.(in Chinese)

    [10]王翠花,包云軒,王建強(qiáng),等.2003年稻縱卷葉螟重大遷入過程的大氣動力機(jī)制分析[J].昆蟲學(xué)報,2006,49(4):604-612.

    Wang C H,Bao Y X,Wang J Q,et al.Analysis on atmospheric dynamical mechanisms for the great immigration events of the rice leaf rollerGuenee in 2003[J].Acta Entomologica Sinica,2006,49(4):604-612.(in Chinese)

    [11]朱金城.稻縱卷葉螟低齡幼蟲防治適期的簡化測報法[J].廣東農(nóng)業(yè)科學(xué),1990,(1):35-37.

    Zhu J C.Simplified forecasting method of rice leaf roller during young larvae control period[J].Guangdong Agricultural Sciences,1990,(1):35-37.(in Chinese)

    [12]程極益.作物病蟲害數(shù)理統(tǒng)計預(yù)報[M].北京:農(nóng)業(yè)出版社,1992.

    Cheng J Y.Mathematical statistics forecasts of plant diseases and insect pests of crop[M].Beijing:Agriculture Press,1992. (in Chinese)

    [13]汪泗水,張孝羲,張夕林.稻縱卷葉螟發(fā)生程度的神經(jīng)網(wǎng)絡(luò)預(yù)警[J].生物數(shù)學(xué)學(xué)報,2003,18(1):93-97.

    Wang S S,Zhang X X,Zhang X L.Neural network warning of the occurrence grade ofGuenee[J]. Journal of Biomathematics,2003,18(1):93-97.(in Chinese)

    [14]劉靜,張宗山,馬力文,等.寧夏枸杞蚜蟲發(fā)生規(guī)律及其氣象等級預(yù)報[J].中國農(nóng)業(yè)氣象,2015,36(3):356-363.

    Liu J,Zhang Z S,Ma L W,et al.Occurrence ofsp. onL.and Its meteorological grades forecast in Ningxia[J].Chinese Journal of Agrometeorology, 2015, 36(3): 356-363.(in Chinese)

    [15]成兆金,薛曉萍,王建源,等.山東棉鈴蟲發(fā)生發(fā)展氣象預(yù)測模型[J].中國農(nóng)業(yè)氣象,2008,29(2):224-229.

    Cheng Z J,Xue X P,Wang J Y,et al.Meteorological prediction model for occurrence and development of bollworm () in Shandong province[J].Chinese Journal of Agrometeorology,2008,29(2):224-229.(in Chinese)

    [16]張淑杰,張玉書,吳微微,等.遼寧省玉米螟發(fā)生程度的氣象預(yù)報模型[J].中國農(nóng)業(yè)氣象,2009,30(2):262-266.

    Zhang S J,Zhang Y S,Wu W W,et al. Meteorological forecast model for emergence grade of corn borer in Liaoning province[J].Chinese Journal of Agrometeorology, 2009,30(2): 262-266.(in Chinese)

    [17]陳惠,關(guān)瑞鋒,楊凱,等.福建省稻飛虱氣象條件適宜程度等級預(yù)報[J].中國農(nóng)業(yè)氣象,2008,29(4):496-498.

    Chen H,Guan R F,Yang K,et al. Forecast of suitable grade of meteorological conditions for occurrences of rice plant hoppers in Fujian province[J].Chinese Journal of Agrometeorology, 2008,29(4):496-498.(in Chinese)

    [18]程進(jìn)偉,李建勛.卡爾曼濾波算法評估平臺的設(shè)計與實現(xiàn)[J].系統(tǒng)仿真學(xué)報,2013,(11):2567-2574.

    Cheng J W,Li J X.Design and implementation of software platform to evaluate Kalman filter algorithm[J].Journal of System Simulation,2013,(11):2567-2574.(in Chinese)

    [19]Battin R H.Astronautical guidance[M].New York:McGraw- Hill,1964:303-340.

    [20]Dove M J,Miller K M.Kalman filters in navigation systems[J].Journal of Navigation,1989,42(2):260-262.

    [21]Gol'tsov A S,Simbirskii D F,Kudryashov S V.Dynamic method of measuring heat fluxes by battery heat flowmeters using a Kalman filter[J].Journal of Engineering Physics,1978, 33(6):1469-1474.

    [22]Kilpinen J.The application of Kalman filter in statistical interpretation of numerical weather forecasts[A].Proceedings of the 12th conference on probability and statistics in atmospheric sciences[C].American Meteorological Society, 1992:11-16.

    [23]江勝國,楊太明,程林,等.卡爾曼濾波方法在稻飛虱發(fā)生等級預(yù)測中的應(yīng)用研究[J].氣象,2010,36(10):106-109.

    Jiang S G,Yang T M,Cheng L,et al.Application of the Kalman filtering method to rice planthopper occurrence grade prediction[J].Meteorological Monthly,2010,36(10): 106-109. (in Chinese)

    [24]彭丁聰.卡爾曼濾波的基本原理及應(yīng)用[J].軟件導(dǎo)刊,2009,8(11):32-34.

    Peng D C.Basic principle and application of Kalman filte[J]. Software Guide,2009,8(11):32-34.(in Chinese)

    [25]唐啟義,李紹石.病蟲測報應(yīng)驗程度判定模式[J].中國植保導(dǎo)刊,1988,8(1):1-5.

    Tang Q Y,Li S S.Evaluation model of degree of fulfillment for forecasting crop disease and pest[J].China Plant Protection, 1988,8(1):1-5.(in Chinese)

    Short-term Forecasting Models on Occurrence of Rice Leaf Roller Based on Kalman Filter Algorithm

    BAO Yun-xuan1, 2, CHEN Xin-yi1, 2, XIE Xiao-jin1, 2, WANG Lin1, 2, LU Ming-hong3

    (1.Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Nanjing University of Information Science & Technology, Nanjing 210044, China; 2.Jiangsu Key Laboratory of Agricultural Meteorology/Nanjing University of Information Science & Technology, Nanjing 210044; 3.National Agricultural Technology Extension and Service Center, Ministry of Agricultural, Beijing 100125)

    In this paper, the pentad systematic investigation data ofat the four representative plant protection stations of four main rice-growing regions (including the rice-growing region of the south China, the rice-growing region of the southwestern China, the rice-growing region between the Nanling mountains and the Yantze River valley and the rice-growing region between the Yantze River valley and the Huaihe River valley) in China was collected from 1994 to 2014, the key meteorological factors influencing onoccurrence amount were screened out and Kalman filter algorithm was used to establish the short-term forecasting models ofpentad occurrence amount at the four plant protection stations, including Quanzhou in the Guangxi Zhuang Autonomous Region, Xiushan in Chongqing city, Xiangyin in Hunan province and Zhangjiagang in Jiangsu province in the immigration and damage period ofrespectively. Based on the back substitution fittings and forecasting tests of the model, the errors and stability and accuracy rates of the Kalman model were calculated. The results showed as follows: (1) for four stations, the occurrence amount ofin the present pentad was significantly and positively correlated (P<0.01) with thes moth amounts of the preceding pentad and the preceding two pentads in the field respectively. There were significantly positive correlations (P<0.01) between the occurrence amounts ofin the present pentad and the minimum air temperature, mean air temperature and maximum air temperature in the preceding pentad. But the pentad occurrence amount was significantly and negatively correlated with the surface pressure in the preceding pentad. (2) The back substitution fitting calculations from the Kalman model on the occurrence amount offrom 1994 to 2011 and the trial forecast tests from 2012 to 2014 showed that the comprehensive mean error (ME) of the occurrence amounts by the Kalman model was -88.63, the mean absolute error (MAE) was 217.72, the comprehensive root mean square error (RMSE) was 605.04, the comprehensive mean accuracies (MA) was 84.33%, and the fitting rate was 83.33%. The Kalman model’s forecasting results were basically consistent with measured values, which indicated that the model could be applied to the prediction of occurrence amount of.

    Guenee; Meteorologic elements; Kalman Filter Algorithm; Pentad forecasting model of occurrence amount; Accuracy rate

    10.3969/j.issn.1000-6362.2016.05.010

    2016-02-23

    國家自然科學(xué)基金面上項目(41475106;41075086);國家公益性行業(yè)(氣象)科研專項(GYHY201306053);江蘇省農(nóng)業(yè)科技自主創(chuàng)新項目[SCX(12)3058];江蘇省高校優(yōu)勢學(xué)科建設(shè)工程

    包云軒(1963-),博士,教授,主要研究方向為氣候變化與防災(zāi)減災(zāi)、應(yīng)用氣象、病蟲害測報學(xué)。E-mail:baoyx@nuist.edu.cn; baoyunxuan@163.com

    猜你喜歡
    卡爾曼發(fā)生量卷葉螟
    太倉市近年稻縱卷葉螟發(fā)生情況及防治進(jìn)展
    棉大卷葉螟識別與為害癥狀
    長江蔬菜(2022年11期)2022-07-29 02:15:42
    基于眼高刻槽數(shù)預(yù)測松褐天牛發(fā)生量分析
    狀態(tài)變換擴(kuò)展卡爾曼平滑算法在AUV水下航跡修正中的應(yīng)用
    六盤山地區(qū)落葉松紅腹葉蜂發(fā)生量對溫度與降水的響應(yīng)研究
    縱卷葉螟絨繭蜂對稻縱卷葉螟的控制作用
    基于卡爾曼算法的動力鋰電池SOC估算
    基于卡爾曼預(yù)測的逆變焊接電源信號處理方法
    焊接(2016年4期)2016-02-27 13:04:27
    基于擴(kuò)展卡爾曼估計的飛機(jī)防滑剎車系統(tǒng)模糊控制
    5%環(huán)蟲酰肼防治稻縱卷葉螟田間藥效評價
    一级毛片久久久久久久久女| 男人和女人高潮做爰伦理| 亚洲av.av天堂| 日本午夜av视频| 熟女电影av网| 国产av精品麻豆| 国产探花极品一区二区| 久久久久久久亚洲中文字幕| 精品人妻一区二区三区麻豆| 丝袜在线中文字幕| 少妇 在线观看| 最后的刺客免费高清国语| 日韩成人av中文字幕在线观看| 日韩中文字幕视频在线看片| 九九久久精品国产亚洲av麻豆| 午夜av观看不卡| a级毛片在线看网站| 日韩欧美 国产精品| 精品熟女少妇av免费看| 日本免费在线观看一区| 精品一区在线观看国产| 午夜免费男女啪啪视频观看| 精品亚洲成a人片在线观看| 伊人久久国产一区二区| 久久精品国产a三级三级三级| 丝袜在线中文字幕| 女的被弄到高潮叫床怎么办| 一级av片app| 免费黄网站久久成人精品| 成年av动漫网址| 成年av动漫网址| 国产欧美日韩一区二区三区在线 | 一级av片app| 建设人人有责人人尽责人人享有的| 色婷婷久久久亚洲欧美| 精品卡一卡二卡四卡免费| 亚洲高清免费不卡视频| 一区二区三区免费毛片| 精品少妇内射三级| 成人无遮挡网站| 日韩欧美精品免费久久| 大香蕉97超碰在线| 国内精品宾馆在线| 最后的刺客免费高清国语| 精品久久久久久电影网| 亚洲电影在线观看av| 亚洲三级黄色毛片| 日韩一区二区视频免费看| 久久狼人影院| 人妻 亚洲 视频| 两个人免费观看高清视频 | 最新中文字幕久久久久| 少妇丰满av| 一区在线观看完整版| 欧美日韩精品成人综合77777| 精品久久久久久久久亚洲| 日韩人妻高清精品专区| 欧美xxⅹ黑人| 亚洲电影在线观看av| 在线天堂最新版资源| 亚洲第一区二区三区不卡| 国产精品一区二区性色av| 亚洲一区二区三区欧美精品| 国产亚洲5aaaaa淫片| 少妇被粗大猛烈的视频| 日韩中字成人| 日本色播在线视频| 少妇被粗大猛烈的视频| 国产精品国产三级专区第一集| 日日啪夜夜撸| 如日韩欧美国产精品一区二区三区 | 日日啪夜夜撸| 日本av免费视频播放| 最新中文字幕久久久久| av福利片在线| 国产精品秋霞免费鲁丝片| 美女主播在线视频| 又大又黄又爽视频免费| 菩萨蛮人人尽说江南好唐韦庄| av免费在线看不卡| 精品久久久久久电影网| av福利片在线| 亚洲av成人精品一二三区| 日韩不卡一区二区三区视频在线| av免费观看日本| 精品国产一区二区久久| 极品人妻少妇av视频| 爱豆传媒免费全集在线观看| 一边亲一边摸免费视频| 精品午夜福利在线看| 夫妻午夜视频| 中国国产av一级| 午夜免费鲁丝| 亚洲欧美日韩卡通动漫| 在线观看免费高清a一片| av国产久精品久网站免费入址| 高清欧美精品videossex| 一区二区三区乱码不卡18| 久久久久网色| 99热6这里只有精品| 国产一级毛片在线| 国产爽快片一区二区三区| 亚洲av在线观看美女高潮| 一级爰片在线观看| 大又大粗又爽又黄少妇毛片口| 日韩视频在线欧美| a级一级毛片免费在线观看| 两个人免费观看高清视频 | 亚洲精品乱久久久久久| 一级毛片电影观看| 国产黄片视频在线免费观看| 美女内射精品一级片tv| 欧美 亚洲 国产 日韩一| 五月天丁香电影| 免费黄色在线免费观看| 成人二区视频| 国产欧美日韩一区二区三区在线 | 春色校园在线视频观看| 亚洲va在线va天堂va国产| 曰老女人黄片| 99热全是精品| 欧美日韩亚洲高清精品| 自线自在国产av| 中文在线观看免费www的网站| 老女人水多毛片| 国产在线一区二区三区精| 日韩电影二区| 亚洲欧美成人精品一区二区| 一个人看视频在线观看www免费| 色网站视频免费| 我的老师免费观看完整版| 一区二区三区精品91| 日韩精品有码人妻一区| 日韩大片免费观看网站| av一本久久久久| 国产精品女同一区二区软件| 91久久精品国产一区二区成人| 亚洲国产精品成人久久小说| a级毛色黄片| 自线自在国产av| 午夜视频国产福利| 色94色欧美一区二区| 最新的欧美精品一区二区| 亚洲怡红院男人天堂| 一边亲一边摸免费视频| 六月丁香七月| 韩国av在线不卡| av线在线观看网站| 久久久精品免费免费高清| 亚洲在久久综合| 日本黄大片高清| 亚洲av国产av综合av卡| 少妇人妻久久综合中文| 国产亚洲91精品色在线| 97在线视频观看| 午夜老司机福利剧场| 91久久精品国产一区二区成人| 精品一区二区三区视频在线| 欧美成人精品欧美一级黄| tube8黄色片| 99久久精品一区二区三区| 日韩,欧美,国产一区二区三区| 久久人妻熟女aⅴ| 欧美日韩在线观看h| 色吧在线观看| 一级二级三级毛片免费看| 人妻人人澡人人爽人人| 国产精品久久久久久精品古装| 欧美一级a爱片免费观看看| 国产伦在线观看视频一区| 精品熟女少妇av免费看| 日本av免费视频播放| 国产成人精品婷婷| 在线观看美女被高潮喷水网站| 亚洲四区av| 久久久久久人妻| 18+在线观看网站| 欧美区成人在线视频| 婷婷色综合大香蕉| 免费av中文字幕在线| 插阴视频在线观看视频| 丝袜在线中文字幕| 国产欧美日韩精品一区二区| 国产av精品麻豆| 3wmmmm亚洲av在线观看| 午夜久久久在线观看| 六月丁香七月| 啦啦啦啦在线视频资源| 欧美精品高潮呻吟av久久| 成人午夜精彩视频在线观看| .国产精品久久| 亚洲一区二区三区欧美精品| 一本色道久久久久久精品综合| 国产综合精华液| 黄色视频在线播放观看不卡| 精品亚洲成a人片在线观看| 在线精品无人区一区二区三| 久久99一区二区三区| 免费播放大片免费观看视频在线观看| 国产男人的电影天堂91| 在线观看免费高清a一片| 色视频www国产| 久久这里有精品视频免费| av天堂中文字幕网| 久久99蜜桃精品久久| 日韩中文字幕视频在线看片| 99九九在线精品视频 | 大陆偷拍与自拍| 亚洲精品,欧美精品| 亚洲内射少妇av| 我要看黄色一级片免费的| 国产精品蜜桃在线观看| 又黄又爽又刺激的免费视频.| 国产一区二区三区综合在线观看 | 国产乱来视频区| 欧美97在线视频| 亚洲成色77777| 最黄视频免费看| 久久人人爽av亚洲精品天堂| 日韩亚洲欧美综合| 久久久久久久国产电影| 国产伦精品一区二区三区四那| av天堂中文字幕网| 观看av在线不卡| 国产视频首页在线观看| 2021少妇久久久久久久久久久| 午夜福利视频精品| 精品久久久久久电影网| 国产黄片美女视频| 三级国产精品片| 美女xxoo啪啪120秒动态图| 国产高清不卡午夜福利| 成年人午夜在线观看视频| videossex国产| 中文字幕av电影在线播放| 51国产日韩欧美| 日韩中文字幕视频在线看片| av一本久久久久| 午夜91福利影院| 成人二区视频| 秋霞伦理黄片| 久久久午夜欧美精品| 精品人妻偷拍中文字幕| 成人国产av品久久久| 黑人猛操日本美女一级片| 亚洲欧美日韩另类电影网站| 日本91视频免费播放| 国产精品一区二区在线观看99| 五月伊人婷婷丁香| 少妇的逼好多水| 国产黄片视频在线免费观看| 免费黄色在线免费观看| 大又大粗又爽又黄少妇毛片口| 高清视频免费观看一区二区| 色网站视频免费| 日韩制服骚丝袜av| 亚洲精品久久午夜乱码| 免费在线观看成人毛片| 韩国av在线不卡| 国产亚洲av片在线观看秒播厂| 国产色婷婷99| 国产毛片在线视频| 午夜日本视频在线| 国产男人的电影天堂91| 夜夜骑夜夜射夜夜干| 日日啪夜夜爽| 看免费成人av毛片| 中文欧美无线码| 免费av中文字幕在线| 男人和女人高潮做爰伦理| 极品教师在线视频| 欧美高清成人免费视频www| av网站免费在线观看视频| 一级片'在线观看视频| 插阴视频在线观看视频| 亚洲av在线观看美女高潮| 精品少妇黑人巨大在线播放| 亚洲精品自拍成人| 国国产精品蜜臀av免费| 下体分泌物呈黄色| 中文字幕久久专区| 大话2 男鬼变身卡| 丰满人妻一区二区三区视频av| 99九九在线精品视频 | 美女主播在线视频| 午夜免费鲁丝| 成人二区视频| 国产一区二区三区综合在线观看 | 极品教师在线视频| www.av在线官网国产| 99热国产这里只有精品6| 免费观看性生交大片5| 香蕉精品网在线| 国产成人精品福利久久| 久久99蜜桃精品久久| 国产综合精华液| 国产一级毛片在线| 国产伦理片在线播放av一区| 欧美亚洲 丝袜 人妻 在线| 精品一区在线观看国产| 亚洲精品成人av观看孕妇| 日韩欧美一区视频在线观看 | 18禁在线无遮挡免费观看视频| 亚洲av不卡在线观看| 人妻一区二区av| 亚洲av国产av综合av卡| 国产精品一区二区在线不卡| 麻豆精品久久久久久蜜桃| 国产成人午夜福利电影在线观看| 午夜日本视频在线| 搡老乐熟女国产| 在线观看www视频免费| 在线精品无人区一区二区三| 免费观看a级毛片全部| 国产高清有码在线观看视频| 国产精品久久久久久久久免| 欧美+日韩+精品| 日日啪夜夜撸| 欧美性感艳星| 亚洲国产欧美日韩在线播放 | 18禁在线播放成人免费| 亚洲成人一二三区av| 久久精品国产亚洲av涩爱| 成人黄色视频免费在线看| 国产成人免费无遮挡视频| 精品人妻一区二区三区麻豆| 亚洲av二区三区四区| 少妇人妻 视频| 色婷婷久久久亚洲欧美| 男人添女人高潮全过程视频| 日韩av在线免费看完整版不卡| 亚洲精品国产色婷婷电影| 国产精品国产三级国产专区5o| 国产精品99久久久久久久久| 亚洲天堂av无毛| 91久久精品国产一区二区成人| videos熟女内射| 国产成人精品无人区| 一级毛片电影观看| 国产无遮挡羞羞视频在线观看| 亚洲中文av在线| 亚洲综合色惰| videossex国产| 亚洲图色成人| 少妇精品久久久久久久| 国产 精品1| 久久精品久久久久久久性| 丝瓜视频免费看黄片| 国产在线男女| av天堂中文字幕网| 国产精品福利在线免费观看| 又粗又硬又长又爽又黄的视频| 久久99热这里只频精品6学生| 热99国产精品久久久久久7| 日本黄色片子视频| 精品亚洲成国产av| 日韩成人伦理影院| 精品国产露脸久久av麻豆| 中文欧美无线码| 亚洲内射少妇av| 王馨瑶露胸无遮挡在线观看| 韩国av在线不卡| 精品亚洲成国产av| 在线观看av片永久免费下载| 亚洲欧洲日产国产| 简卡轻食公司| 亚洲无线观看免费| 99热国产这里只有精品6| 日韩欧美精品免费久久| 中文字幕亚洲精品专区| 黄色一级大片看看| 亚洲一级一片aⅴ在线观看| 国产黄片视频在线免费观看| 欧美日韩亚洲高清精品| 欧美日韩综合久久久久久| 国产国拍精品亚洲av在线观看| 在线 av 中文字幕| 另类亚洲欧美激情| 欧美bdsm另类| 日韩亚洲欧美综合| 内射极品少妇av片p| 国产成人免费观看mmmm| 免费在线观看成人毛片| 岛国毛片在线播放| 又粗又硬又长又爽又黄的视频| 三级国产精品片| 亚洲av二区三区四区| 精品久久久久久久久av| 男人狂女人下面高潮的视频| 伦理电影大哥的女人| 亚洲av福利一区| a级毛色黄片| 欧美区成人在线视频| 一级毛片久久久久久久久女| 丝袜喷水一区| 亚洲美女黄色视频免费看| 视频中文字幕在线观看| 精品久久国产蜜桃| 欧美激情国产日韩精品一区| 五月开心婷婷网| 久久久久视频综合| 建设人人有责人人尽责人人享有的| 青春草视频在线免费观看| 久久久久久久久久久丰满| 22中文网久久字幕| 精品卡一卡二卡四卡免费| 在线观看免费视频网站a站| 搡女人真爽免费视频火全软件| 夜夜骑夜夜射夜夜干| 久久亚洲国产成人精品v| 亚洲一级一片aⅴ在线观看| 熟女人妻精品中文字幕| av.在线天堂| 久久人人爽av亚洲精品天堂| 一级爰片在线观看| 亚洲美女黄色视频免费看| 能在线免费看毛片的网站| 亚洲欧美日韩东京热| 日韩三级伦理在线观看| 三级国产精品片| 最近中文字幕2019免费版| 夫妻午夜视频| 有码 亚洲区| 夜夜骑夜夜射夜夜干| 精品少妇黑人巨大在线播放| 黄色日韩在线| 成人黄色视频免费在线看| 成人漫画全彩无遮挡| 老熟女久久久| 国产亚洲5aaaaa淫片| 久久国内精品自在自线图片| 人人妻人人看人人澡| 我的女老师完整版在线观看| 青春草国产在线视频| 亚洲经典国产精华液单| 久久人人爽人人爽人人片va| 日日啪夜夜爽| 欧美bdsm另类| 久久精品久久久久久久性| 天天躁夜夜躁狠狠久久av| 能在线免费看毛片的网站| 免费大片黄手机在线观看| 成人18禁高潮啪啪吃奶动态图 | 啦啦啦视频在线资源免费观看| 在线播放无遮挡| 色哟哟·www| 日本91视频免费播放| 搡女人真爽免费视频火全软件| 亚洲精品中文字幕在线视频 | 青青草视频在线视频观看| 成人午夜精彩视频在线观看| 亚洲精品色激情综合| 又爽又黄a免费视频| 美女主播在线视频| 在线观看www视频免费| 亚洲人成网站在线观看播放| 欧美精品国产亚洲| 精品亚洲成国产av| 亚洲国产av新网站| 国产91av在线免费观看| 亚洲熟女精品中文字幕| 人妻系列 视频| 欧美精品一区二区大全| 性色avwww在线观看| 久久久久久久久久久免费av| 精品少妇久久久久久888优播| 久久亚洲国产成人精品v| 97在线视频观看| 精品视频人人做人人爽| 国产一区亚洲一区在线观看| 亚洲国产精品国产精品| 最新中文字幕久久久久| 免费看日本二区| 免费高清在线观看视频在线观看| 人妻少妇偷人精品九色| 我的老师免费观看完整版| 一本大道久久a久久精品| 欧美区成人在线视频| 少妇猛男粗大的猛烈进出视频| 精品国产露脸久久av麻豆| 精品99又大又爽又粗少妇毛片| 在线看a的网站| 最后的刺客免费高清国语| 亚洲欧美日韩东京热| 中国国产av一级| h视频一区二区三区| 七月丁香在线播放| 国产欧美日韩一区二区三区在线 | 成年人免费黄色播放视频 | 国产 精品1| 亚洲欧美精品自产自拍| 成人无遮挡网站| 免费看光身美女| 欧美性感艳星| 日韩欧美一区视频在线观看 | 极品少妇高潮喷水抽搐| 观看av在线不卡| 秋霞在线观看毛片| 国内少妇人妻偷人精品xxx网站| 一本久久精品| 伦理电影大哥的女人| 观看美女的网站| 久久av网站| 亚洲欧洲日产国产| 久久国产精品男人的天堂亚洲 | 午夜影院在线不卡| 日本wwww免费看| 亚洲精品456在线播放app| 中文天堂在线官网| 日本欧美国产在线视频| 免费av不卡在线播放| 亚洲色图综合在线观看| 亚洲性久久影院| 一级二级三级毛片免费看| 男女无遮挡免费网站观看| 桃花免费在线播放| 97在线人人人人妻| 精品久久久久久电影网| 成年av动漫网址| 视频区图区小说| 国产一级毛片在线| 亚洲精品,欧美精品| 免费在线观看成人毛片| 亚洲精品aⅴ在线观看| 成年人午夜在线观看视频| 免费观看在线日韩| 一本大道久久a久久精品| 大陆偷拍与自拍| 热re99久久国产66热| 久久久国产精品麻豆| 卡戴珊不雅视频在线播放| 欧美激情极品国产一区二区三区 | 看免费成人av毛片| 青青草视频在线视频观看| 日日摸夜夜添夜夜添av毛片| 在线观看国产h片| 91精品一卡2卡3卡4卡| 久久午夜综合久久蜜桃| 亚洲精华国产精华液的使用体验| 欧美另类一区| 国产淫语在线视频| av又黄又爽大尺度在线免费看| 少妇精品久久久久久久| 一本色道久久久久久精品综合| tube8黄色片| 久久久久视频综合| 青春草视频在线免费观看| 九色成人免费人妻av| 成年女人在线观看亚洲视频| 少妇的逼好多水| 国产欧美日韩精品一区二区| 18禁裸乳无遮挡动漫免费视频| 又黄又爽又刺激的免费视频.| 爱豆传媒免费全集在线观看| 久久国产精品男人的天堂亚洲 | 人妻 亚洲 视频| 国产亚洲欧美精品永久| 国产淫语在线视频| 午夜久久久在线观看| 在线观看免费视频网站a站| 国产精品欧美亚洲77777| 国产欧美日韩精品一区二区| 91久久精品国产一区二区三区| 另类亚洲欧美激情| 亚洲av电影在线观看一区二区三区| 人人妻人人看人人澡| 亚洲av电影在线观看一区二区三区| 国产精品不卡视频一区二区| 一级毛片久久久久久久久女| 国产精品99久久久久久久久| av国产精品久久久久影院| h日本视频在线播放| 好男人视频免费观看在线| 最近中文字幕2019免费版| 夜夜骑夜夜射夜夜干| 男女免费视频国产| 国产成人aa在线观看| 日韩av不卡免费在线播放| 国产美女午夜福利| 汤姆久久久久久久影院中文字幕| 青春草国产在线视频| 观看av在线不卡| 久久久久久人妻| 亚洲av男天堂| 欧美精品高潮呻吟av久久| 日韩av免费高清视频| 亚洲图色成人| 亚洲精品乱久久久久久| 男人添女人高潮全过程视频| 91久久精品国产一区二区三区| 中国美白少妇内射xxxbb| 日本爱情动作片www.在线观看| 午夜福利,免费看| 久久97久久精品| 美女脱内裤让男人舔精品视频| av不卡在线播放| 国产免费一级a男人的天堂| 国产精品一区www在线观看| 国产精品国产三级专区第一集| 超碰97精品在线观看| 国产黄片美女视频| 国模一区二区三区四区视频| 精品熟女少妇av免费看| 亚洲人与动物交配视频| 国产精品福利在线免费观看| 亚洲av不卡在线观看| 99热国产这里只有精品6| 丁香六月天网| 大片电影免费在线观看免费| 80岁老熟妇乱子伦牲交| 黄色怎么调成土黄色| 国产白丝娇喘喷水9色精品| 中文字幕精品免费在线观看视频 | 少妇 在线观看| 免费观看a级毛片全部|