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

    北京地區(qū)秋季日光溫室黃瓜白粉病預(yù)測(cè)的貝葉斯網(wǎng)絡(luò)模型研究

    2022-03-12 19:10:12魏少偉任愛(ài)新趙靖暄楊信廷李明劉慧英
    中國(guó)瓜菜 2022年2期
    關(guān)鍵詞:貝葉斯網(wǎng)絡(luò)預(yù)測(cè)模型白粉病

    魏少偉 任愛(ài)新 趙靖暄 楊信廷 李明 劉慧英

    摘 ? ?要:為了準(zhǔn)確預(yù)測(cè)日光溫室黃瓜白粉病的發(fā)生,以日光溫室水果黃瓜為試驗(yàn)材料,于2020年9—11月在北京市不同方位的4個(gè)日光溫室內(nèi),運(yùn)用無(wú)線網(wǎng)絡(luò)環(huán)境監(jiān)測(cè)系統(tǒng)對(duì)日光溫室黃瓜的生長(zhǎng)環(huán)境(空氣溫度、相對(duì)濕度、光照強(qiáng)度)進(jìn)行了實(shí)時(shí)動(dòng)態(tài)監(jiān)測(cè),并同步進(jìn)行白粉病流行調(diào)查。采用貝葉斯網(wǎng)絡(luò)模型建立日光溫室黃瓜白粉病預(yù)測(cè)模型,預(yù)測(cè)黃瓜白粉病是否發(fā)病,并與田間觀測(cè)的實(shí)際發(fā)病情況進(jìn)行比較。結(jié)果表明,模型在4個(gè)溫室和整體預(yù)測(cè)發(fā)生準(zhǔn)確度 ACC(Accuracy)分別為0.95、0.92、0.91、0.87、0.87,約登指數(shù)J(Youden Index)分別為0.90、0.86、0.84、0.70、0.74。表明模型在預(yù)測(cè)整體和各個(gè)溫室效果較好,具有良好的普適性,可為實(shí)際生產(chǎn)中黃瓜白粉病預(yù)測(cè)提供參考。

    關(guān)鍵詞: 黃瓜;白粉病;貝葉斯網(wǎng)絡(luò);預(yù)測(cè)模型

    中圖分類號(hào):S642.2 文獻(xiàn)標(biāo)志碼:A 文章編號(hào):1673-2871(2022)02-020-08

    Prediction of cucumber powdery mildew in autumn solar greenhouse in Beijing based on Bayesian network model

    WEI Shaowei1,2, REN Aixin2, ZHAO Jingxuan2, YANG Xinting2, LI Ming1,2, LIU Huiying1

    (1. Key Laboratory of Cultivation Physiology and Germplasm Resources Utilization of Featured Fruits and Vegetables of Xinjiang Production and Construction Corps/College of Agriculture, Shihezi University, Shihezi 832003, Xinjiang, China; 2. National Engineering Laboratory for Quality and Safety Traceability Technology and Application of Agricultural Products/Information Technology Research Center, Beijing Academy of Agriculture and Forestry Sciences /National Engineering Research Center of Information Technology in Agriculture/Meteorological Service Center of Urban Agriculture, China Meteorological Administration-Ministry of Agriculture and Rural Affairs, Beijing 100097, China)

    Abstract: To accurately predict the occurrence of cucumber powdery mildew in solar greenhouse, four solar greenhouses were chosen in Beijing from September to November 2020, using the wireless network environmental monitoring system and real-time dynamic monitoring to collect cucumber growth environment (air temperature, humidity, light intensity) and powdery mildew data. A forecast model of cucumber powdery mildew in solar greenhouse was developed using the Bayesian network model, the incidence of cucumber powdery mildew predicted was and compared with field observation. The model performed well in four greenhouses with the accuracy ACC= 0.95, 0.92, 0.91, 0.87, 0.87, and the Youden Index J=0.90, 0.86, 0.84, 0.70, 0.74; The results showed that the model had a good prediction at all greenhouse and each greenhouse level, and could provide a guide for powdery mildew management in the actual production.

    Key words: Cucumber; Powdery mildew; Bayesian network; Prediction model

    白粉病(Erysiphe cucurbitacearum Zheng & Chen)是瓜類生產(chǎn)上的重要病害,在我國(guó)各地保護(hù)地和陸地的多種瓜菜上普遍發(fā)生[1]。一般在黃瓜生長(zhǎng)中、后期病情發(fā)展傳播迅速,會(huì)導(dǎo)致大量減產(chǎn),造成較為嚴(yán)重的經(jīng)濟(jì)損失[2]。若不及時(shí)防治會(huì)造成整株被白粉病層覆蓋,不僅會(huì)降低植株的光合效率,還會(huì)導(dǎo)致整株枯萎和植株過(guò)早衰老,嚴(yán)重影響黃瓜的品質(zhì)和產(chǎn)量[3]。病害的早期預(yù)防及防治對(duì)病害的發(fā)生及控制極其關(guān)鍵,因此在病害發(fā)生之前對(duì)黃瓜白粉病進(jìn)行準(zhǔn)確快速的預(yù)測(cè)具有重要意義。

    黃瓜白粉病流行趨勢(shì)主要受寄主、氣候條件、病原三方面的因素影響[4]。其中,溫度[5]、光照[6]和濕度[7]是影響該病害流行的重要因素?;谶@些影響因素,植保專家和學(xué)者結(jié)合多年的田間生產(chǎn)經(jīng)驗(yàn)總結(jié)和建立了多種黃瓜白粉病的田間預(yù)測(cè)預(yù)報(bào)方法,如卷積神經(jīng)網(wǎng)絡(luò)[8]、逐步回歸法[9]、圖像識(shí)別[10]、物聯(lián)網(wǎng)技術(shù)[11]等。徐寧[9]通過(guò)人工接種和模擬高溫高濕環(huán)境的方法研究了黃瓜白粉病的發(fā)生和流行趨勢(shì),采用回歸分析,建立黃瓜白粉病始病期與流行程度預(yù)測(cè),結(jié)果表明效果模型較好。李磊福等[12]采用逐步回歸分析法構(gòu)建溫室甜瓜白粉病的流行預(yù)測(cè)模型。結(jié)果表明,該模型具有很強(qiáng)的實(shí)用價(jià)值。尹哲等[13]采用回歸分析技術(shù),建立了溫室大棚中黃瓜白粉病流行趨勢(shì)(病情指數(shù))與濕度、溫度間關(guān)系的預(yù)測(cè)模型,結(jié)果表明模型預(yù)測(cè)效果比較好。

    針對(duì)黃瓜白粉病的預(yù)測(cè)大多采用多元回歸等方法,但是這些方法無(wú)法揭示變量之間的因果關(guān)系,甚至不能進(jìn)行推理驗(yàn)證。而貝葉斯網(wǎng)絡(luò)卻可以反映變量之間的因果關(guān)系,還能進(jìn)行推理[14],并且在解決問(wèn)題方面具有明確、直觀等優(yōu)點(diǎn),已成為各界學(xué)者的研究熱點(diǎn),在農(nóng)牧業(yè)領(lǐng)域研究比較多,包括農(nóng)產(chǎn)品質(zhì)量安全預(yù)測(cè)[15]、農(nóng)業(yè)環(huán)境分析[16]、大田作物病害預(yù)測(cè)[17-18]以及農(nóng)業(yè)工程中的故障診斷[19]等;迄今,在包括小麥[20]、番茄[21]、葡萄[22]等多種植物上,已經(jīng)成功建立了白粉病的貝葉斯網(wǎng)絡(luò)預(yù)測(cè)模型,均在白粉病的預(yù)測(cè)預(yù)報(bào)方面取得了超越傳統(tǒng)方法的良好效果。然而,通過(guò)貝葉斯網(wǎng)絡(luò)建模來(lái)預(yù)測(cè)溫室黃瓜白粉病的發(fā)生情況尚未見(jiàn)報(bào)道。因此,筆者以溫室黃瓜種植期間的環(huán)境和病害發(fā)生情況數(shù)據(jù)作為樣本,采用貝葉斯網(wǎng)絡(luò)建模技術(shù),建立了溫室黃瓜白粉病預(yù)測(cè)貝葉斯網(wǎng)絡(luò)模型,為病害預(yù)防奠定基礎(chǔ)。

    1 材料與方法

    1.1 材料

    供試黃瓜品種:京研迷你2號(hào)水果黃瓜,屬于中抗白粉病,由國(guó)家蔬菜工程技術(shù)研究中心提供。儀器:溫室小氣候監(jiān)測(cè)點(diǎn)(Enviro monitor,Davis instruments,Hayward,California,USA),Enviro Monitor 無(wú)線網(wǎng)絡(luò)環(huán)境監(jiān)測(cè)系統(tǒng)(美國(guó) Davis Instruments)。

    1.2 溫室黃瓜定植

    試驗(yàn)于2020年9—11月在北京市昌平區(qū)小湯山國(guó)家精準(zhǔn)農(nóng)業(yè)試驗(yàn)基地(5號(hào))溫室、海淀區(qū)首農(nóng)莊園(9號(hào))溫室、房山區(qū)弘科農(nóng)場(chǎng)日光(6號(hào))溫室、大興區(qū)裕農(nóng)公司(7號(hào))溫室內(nèi)中進(jìn)行,溫室由聚乙烯薄膜覆蓋,東西延長(zhǎng)坐北朝南。于9月定植,并建立生產(chǎn)檔案,如實(shí)記錄灌溉、施肥、夜間通風(fēng)等栽培管理措施。

    1.3 溫室黃瓜白粉病調(diào)查

    調(diào)查采用五點(diǎn)取樣法,沿對(duì)角線選取5個(gè)調(diào)查樣點(diǎn),每個(gè)樣點(diǎn)標(biāo)記12株,發(fā)病前每隔2 d觀察1次,待發(fā)現(xiàn)白粉病粉層后每隔4~5 d觀察1次,統(tǒng)計(jì)發(fā)病情況并根據(jù) GB/T 17980.26—2000 對(duì)病害嚴(yán)重度進(jìn)行分級(jí)和記錄。所有溫室黃瓜植株的白粉病病情指數(shù)≥5%時(shí)結(jié)束調(diào)查。

    發(fā)病程度以病情指數(shù)表示,病情指數(shù)計(jì)算公式為:

    病情指數(shù)=[Σ(各級(jí)病葉數(shù)×相對(duì)代表值)調(diào)查總?cè)~數(shù)×9]×

    100%。 ? ? ? ? ? ? ? ?(1)

    黃瓜白粉病分級(jí)標(biāo)準(zhǔn)如下[2]:

    0 級(jí):無(wú)病斑;1 級(jí):病斑面積占整個(gè)葉面積 5%及以下;3 級(jí):病斑面積占整個(gè)葉面積 6%~10%;5 級(jí):病斑面積占整個(gè)葉面積 l1%~25%;7 級(jí):病斑面積占整個(gè)葉面積 26%~50%;9 級(jí):病斑面積占整個(gè)葉面積 51%及以上。

    筆者調(diào)查小湯山基地、弘科農(nóng)場(chǎng)、首農(nóng)莊園、裕農(nóng)公司白粉病發(fā)病情況,待發(fā)現(xiàn)白粉病粉層到病害調(diào)查結(jié)束共調(diào)查4次,根據(jù)調(diào)查情況進(jìn)行病情指數(shù)計(jì)算,所得病情指數(shù)為 9.90%、8.37%、9.07%、7.12%,結(jié)束調(diào)查,并及時(shí)進(jìn)行防治。

    1.4 溫室環(huán)境數(shù)據(jù)采集

    為了采集溫室氣溫、相對(duì)濕度、太陽(yáng)輻射強(qiáng)度和地溫等環(huán)境數(shù)據(jù),在溫室病害調(diào)查點(diǎn)安裝溫室小氣候監(jiān)測(cè)節(jié)點(diǎn)(Enviromonitor,Davis instruments,Hayward,California,USA)。儀器采用 Enviro Monitor 無(wú)線網(wǎng)絡(luò)環(huán)境監(jiān)測(cè)系統(tǒng)(美國(guó) Davis Instruments),溫室內(nèi)儀器設(shè)置記錄數(shù)據(jù)時(shí)間間隔均為 15 min,每 15 min對(duì)溫室內(nèi)的氣溫、相對(duì)濕度、太陽(yáng)輻射強(qiáng)度和地溫等數(shù)據(jù)實(shí)時(shí)觀測(cè)并儲(chǔ)存,數(shù)據(jù)通過(guò)設(shè)備上的無(wú)線傳輸模塊和手動(dòng)下載相結(jié)合來(lái)收集。

    1.5 模型驅(qū)動(dòng)因子的篩選

    不同黃瓜品種在抗病原菌侵入方面差異不顯著,而在病害擴(kuò)展速度上有顯著差異,所以品種抗性不作為預(yù)測(cè)模型的驅(qū)動(dòng)因子,選用溫度、濕度作為預(yù)測(cè)模型的驅(qū)動(dòng)因子[23]。此外,太陽(yáng)輻射強(qiáng)度和地溫也是黃瓜白粉病發(fā)生的影響因素。因此,選擇氣溫、相對(duì)濕度、太陽(yáng)輻射強(qiáng)度和土壤溫度作為模型的驅(qū)動(dòng)因子。

    1.6 模型評(píng)價(jià)

    筆者在本試驗(yàn)中采用國(guó)際上廣泛使用的模型檢驗(yàn)方法:R2(決定系數(shù))計(jì)算模型的擬合度,均方誤差MSE(Mean squared error)評(píng)估模型的精度來(lái)描述試驗(yàn)樣本。均方根誤差RMSE(Root mean squared error)對(duì)預(yù)測(cè)值與實(shí)際值之間的符合度進(jìn)行統(tǒng)計(jì)分析,用準(zhǔn)確度ACC(Accuracy,公式(2))和約登指數(shù) J(Youden Index, 公式(3))來(lái)評(píng)價(jià),其中TP(True positive)是真陽(yáng)性比例,F(xiàn)N(False negative)是假陰性比例,F(xiàn)P(False positive)是假陽(yáng)性比例,TN(True negative)是真陰性比例[24]。

    [ACC=TP±TNTP±FN±FP±TN]; ? ? ?(2)

    [J=TP-FP。] ? ? ?(3)

    1.7 數(shù)據(jù)分析

    貝葉斯網(wǎng)絡(luò)建模采用Pyton 3.7和Pycham 2019.3軟件,Micro-soft Excel 2019 軟件進(jìn)行數(shù)據(jù)整理分析。

    2 結(jié)果與分析

    2.1 黃瓜白粉病預(yù)測(cè)模型構(gòu)建

    2.1.1? ? 模型選取和變量定義? ? 將黃瓜白粉病是否發(fā)生作為因變量,空氣溫度、空氣相對(duì)濕度、土壤溫度、太陽(yáng)輻射強(qiáng)度作為自變量,進(jìn)行模型分析。黃瓜白粉病是否發(fā)生的分析屬于二分類因變量的分析,本試驗(yàn)中選擇貝葉斯網(wǎng)絡(luò)模型進(jìn)行分析。由于變量參數(shù)的記錄值數(shù)量范圍不一致,因此首先對(duì)原始數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理(表1)。

    2.1.2? ? 數(shù)據(jù)的選取及相關(guān)分析 ? ?選取2020年9—11月小湯山基地、首農(nóng)莊園、弘科農(nóng)場(chǎng)、裕農(nóng)公司各溫室內(nèi)定植1周后至發(fā)病前5個(gè)監(jiān)測(cè)點(diǎn)的環(huán)境數(shù)據(jù),將傳感器每15 min采集一次數(shù)據(jù)的環(huán)境數(shù)據(jù)通過(guò)計(jì)算處理成每1 h試驗(yàn)數(shù)據(jù)共7560次數(shù)據(jù),將80%的數(shù)據(jù)作為訓(xùn)練樣本,20%的數(shù)據(jù)作為測(cè)試樣本。采用Pearson相關(guān)分析方法對(duì)數(shù)據(jù)進(jìn)行分析,分析各變量之間的相關(guān)程度并確定任意兩個(gè)變量之間是否存在邊,在貝葉斯網(wǎng)絡(luò)圖中,邊代表了兩個(gè)變量之間具有關(guān)聯(lián)性。根據(jù)表2的結(jié)果并結(jié)合病害發(fā)生規(guī)律,按相關(guān)系數(shù)從大到小的順序給邊定向,如果在網(wǎng)絡(luò)中出現(xiàn)環(huán),舍棄相關(guān)性比較小的兩個(gè)變量之間的邊[28]。

    2.1.3? ? 貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)學(xué)習(xí)? ? 本試驗(yàn)中利用節(jié)點(diǎn)間的因果關(guān)系,采用爬山搜索算法計(jì)算每個(gè)節(jié)點(diǎn)元素的條件獨(dú)立性,確定初始貝葉斯網(wǎng)絡(luò)結(jié)構(gòu),然后利用MDL(minimum definition length)評(píng)分爬山搜索算法對(duì)貝葉斯網(wǎng)絡(luò)進(jìn)行全局優(yōu)化。為了消除模型學(xué)習(xí)中的“過(guò)擬合”問(wèn)題,采用十折交叉驗(yàn)證的方法來(lái)學(xué)習(xí)貝葉斯網(wǎng)絡(luò)結(jié)構(gòu),然后得到MDL評(píng)分估計(jì)值最高的網(wǎng)絡(luò)結(jié)構(gòu)作為最終的貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)[29](如圖1所示)。

    2.1.4? ? 貝葉斯網(wǎng)絡(luò)參數(shù)學(xué)習(xí)? ? ?貝葉斯網(wǎng)絡(luò)學(xué)習(xí)參數(shù)是在已知貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)的基礎(chǔ)上,根據(jù)數(shù)據(jù)對(duì)網(wǎng)絡(luò)中各個(gè)節(jié)點(diǎn)變量進(jìn)行概率計(jì)算。貝葉斯估計(jì)和極大似然估計(jì)是目前參數(shù)學(xué)習(xí)的主要方法 [29]。筆者利用頻數(shù)法[30]的方法建立 BN(Bayesian network)各節(jié)點(diǎn)變量的概率表(CPT)。如圖1所示,本文中用SR 表示太陽(yáng)輻射強(qiáng)度,Temp 表示溫度,Hum表示相對(duì)濕度;Disease 表示黃瓜白粉病發(fā)病狀態(tài)。例如,圖1中的太陽(yáng)輻射強(qiáng)度沒(méi)有父節(jié)點(diǎn),那么它的邊緣概率就等于太陽(yáng)輻射強(qiáng)度大于115 W·m-2的次數(shù)除以總數(shù)據(jù)數(shù)。

    [P(SR=1)=(SR=1)(SR)=0.234]。 (4)

    這樣,[PSR=0=1-PSR=1=1-0.234=0.766]。

    而對(duì)于節(jié)點(diǎn)Disease黃瓜白粉病發(fā)病狀態(tài),其父節(jié)點(diǎn)為空氣溫度和相對(duì)濕度情況,那么它的條件概率:

    [PDisease=1|Temp=1,Hum=1=(Disease=1|Temp=1,Hum=1)(Temp=1,Hum=1)=0.886]。(5)

    同理可計(jì)算出其他節(jié)點(diǎn)的概率表(圖1)。

    2.1.5? ? 貝葉斯網(wǎng)絡(luò)推理? ? 貝葉斯網(wǎng)絡(luò)推理是利用已知節(jié)點(diǎn)的貝葉斯公式來(lái)求解其他節(jié)點(diǎn)的條件似然性。推理方法可分為兩類:精確推理和近似推理[30]。圖1是1個(gè)黃瓜白粉病病害貝葉斯網(wǎng)絡(luò)模型,包含4個(gè)節(jié)點(diǎn)變量,每個(gè)節(jié)點(diǎn)有2個(gè)狀態(tài)(0或1)。表3是節(jié)點(diǎn)SR的先驗(yàn)概率表,表4、表5、表6是節(jié)點(diǎn) Temp、Hum、Disease的條件概率表。利用這個(gè)模型可以推算出給定證據(jù)下任何節(jié)點(diǎn)的概率,其基本原理是貝葉斯理論。例如,如果知道節(jié)點(diǎn)Temp發(fā)生,即Temp=1,要計(jì)算節(jié)點(diǎn)SR 的概率,因?yàn)楣?jié)點(diǎn)Temp 和SR 相關(guān),根據(jù)貝葉斯公式計(jì)算條件概率

    [PSR|Temp=PTemp|SR PSRPTemp]。 (6)

    從表3和表4知道 P(SR)=0.234和P(Temp |SR)= 0.920,所以公式(6)的分子等于0.215,其分母是邊界概率 [PTemp],由于節(jié)點(diǎn)Temp只和節(jié)點(diǎn)SR 有關(guān)系,所以有

    P(Temp=1)= P(Temp=1|SR=1)P(SR=1)+P(Temp=1|SR=0)P(SR=0)=0.276,則 P(Temp=0)=0.724。

    將值P(Temp=1)=0.276代入式(6)得到 [PSR|Temp=0.779]。

    根據(jù)全概率公式(7)計(jì)算得出P(Hum=1)=0.338,則 P(Hum=0)=0.662。

    P(Hum=1)=P(Hum=1|Temp=1,SR=1)P(Temp=1)P(SR=1)

    +P(Hum=1|Temp=1,SR=0)P(Temp=1)P(SR=0)

    +P(Hum=1|Temp=0,SR=1)P(Temp=0)P(SR=1)

    +P(Hum=1|Temp=0,SR=0)P(Temp=0)P(SR=0)? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?(7)

    計(jì)算節(jié)點(diǎn) Disease 的概率,由于節(jié)點(diǎn) Disease 和節(jié)點(diǎn)Temp,Hum 相關(guān),根據(jù)全概率公式(8)得出P(Disease=1)= 0.694。

    P(Disease=1)=P(Disease=1|Temp=1, Hum=1)P(Temp=1)P(Hum=1)

    +P(Disease=1|Temp=1,Hum=0)P(Temp=1)P(Hum=0)

    +P(Disease=1|Temp=0,Hum=1)P(Temp=0)P(Hum=1)

    +P(Disease=1|Temp=0,Hum=0)P(Temp=0)P(Hum=0)? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?(8)

    [PTemp|Disease=PDisease|Temp PTempPDisease]。

    (9)

    根據(jù)公式(9),當(dāng)病害發(fā)生,則滿足病害發(fā)生的溫度概率為0.449。

    根據(jù)公式(10),病害發(fā)生,則滿足病害發(fā)生的相對(duì)濕度概率為0.730。

    [PHumD|isease=PDisease|Hum PHumPDisease]。

    (10)

    貝葉斯網(wǎng)絡(luò)可以進(jìn)行雙向推理,既可以從結(jié)果推理原因也可以從原因推理結(jié)果。例如可以從模型結(jié)構(gòu)概率表進(jìn)行推理,當(dāng)溫度、濕度在滿足什么條件下,病害發(fā)生的概率是多少,或者當(dāng)病害發(fā)生時(shí),溫度和濕度發(fā)生概率是多少;已知一個(gè)節(jié)點(diǎn)概率去推理另外一個(gè)節(jié)點(diǎn)發(fā)生的概率;當(dāng)新的變量發(fā)生改變、加入或減少時(shí),將改變貝葉斯網(wǎng)絡(luò)中的其他變量值。假設(shè)節(jié)點(diǎn)Temp發(fā)生新的變化,即Temp=1,知道P(SR=1)從0.234變成0.435,P(SR=0)也就從0.766變成0.565,將這些值重新代入公式中,得到修改后的邊界概率即為條件概率,而其他節(jié)點(diǎn)及病害發(fā)生的概率隨之發(fā)生變化。將新的發(fā)生數(shù)據(jù)輸入網(wǎng)絡(luò)中更新各個(gè)節(jié)點(diǎn)的概率這個(gè)過(guò)程稱為概率繁殖[30]。這樣的雙向推理能力是其他模型方法無(wú)法實(shí)現(xiàn)的[30]。

    2.2 模型預(yù)測(cè)結(jié)果

    為了檢驗(yàn)建立的貝葉斯網(wǎng)絡(luò)模型的預(yù)測(cè)的準(zhǔn)確性,本試驗(yàn)選取小湯山基地260次數(shù)據(jù)、首農(nóng)莊園274次數(shù)據(jù)、弘科農(nóng)場(chǎng)806次數(shù)據(jù)、裕農(nóng)公司273次數(shù)據(jù),共計(jì)1613次數(shù)據(jù)進(jìn)行混淆矩陣(Confusion Matrix)對(duì)比分析對(duì)預(yù)測(cè)病害發(fā)生和實(shí)際發(fā)生的結(jié)果進(jìn)行評(píng)價(jià),還建立了logistic回歸模型并進(jìn)行比較,兩種預(yù)測(cè)結(jié)果的比較見(jiàn)表7,以進(jìn)一步驗(yàn)證貝葉斯網(wǎng)絡(luò)模型的預(yù)測(cè)效果。從表7可以看出平均BN模型的真陽(yáng)性比例(TP=0.88)高于LR模型(TP=0.57);BN模型的假陰性比例(FN=0.12)低于LR模型(FN=0.43),BN模型比LR模型預(yù)測(cè)病害發(fā)生的準(zhǔn)確率更高,漏報(bào)率更低。BN模型的準(zhǔn)確度(ACC=0.87)優(yōu)于LR模型(ACC=0.66),BN模型的約登指數(shù)(J=0.74)高于LR模型(J=0.25),表明BN模型預(yù)測(cè)結(jié)果真實(shí)性高于LR模型。這說(shuō)明模型在預(yù)測(cè)整體和各個(gè)溫室表現(xiàn)較好,具有良好的普適性,為實(shí)際生產(chǎn)中黃瓜白粉病預(yù)測(cè)提供參考。

    同時(shí)繪制各溫室病害發(fā)生的ROC曲線( Receiver Operating Characteristic Curve,即受試者工作特征曲線)綜合評(píng)價(jià)黃瓜白粉病病害預(yù)測(cè)效果(圖 2)。以真陽(yáng)性率(靈敏度)為縱坐標(biāo),假陽(yáng)性率(1-特異性)為橫坐標(biāo),繪制ROC曲線。曲線下區(qū)域越大(AUC范圍為0~1),則表示該模型的綜合預(yù)測(cè)效果越好[30]。由圖2可以看出,本試驗(yàn)中建立的BN模型預(yù)測(cè)平均命中率在85%以上,反應(yīng)該模型具有較好的預(yù)測(cè)精度,通常80%以上命中率即可被認(rèn)為是理想結(jié)果[31]。同時(shí)根據(jù)表3的結(jié)果和各地區(qū)ROC曲線面積均大于85%,該模型在整體病害預(yù)測(cè)和各地區(qū)預(yù)測(cè)病害的效果比較好,具有一定的普適性,可以進(jìn)一步利用貝葉斯網(wǎng)絡(luò)構(gòu)建病害預(yù)警模型。

    2.3 貝葉斯網(wǎng)絡(luò)模型效果評(píng)估

    BN和LR的統(tǒng)計(jì)指標(biāo)和誤差如表8所示,RMSE衡量觀測(cè)值與真實(shí)值之間的偏差,作為機(jī)器學(xué)習(xí)模型預(yù)測(cè)結(jié)果衡量的標(biāo)準(zhǔn)之一,其值越小,表明模型的偏差就越小,預(yù)測(cè)的準(zhǔn)確性越好, MSE評(píng)價(jià)模型的準(zhǔn)確性,MSE值越小,反映模型預(yù)測(cè)效果越好,準(zhǔn)確性越高。從表8中可以得出BN模型的總體,小湯山基地、首農(nóng)莊園、裕農(nóng)公司、弘科農(nóng)場(chǎng)的MSE、RMSE值均低于LR模型值,其值均小于1。小湯山基地BN決定系數(shù)(R2=0.78)>LR模型(R2=0.38)、首農(nóng)莊園BN決定系數(shù)(R2=0.61)>LR模型(R2=0.49)、弘科農(nóng)場(chǎng)BN決定系數(shù)(R2=0.61)>LR模型(R2=0.06),BN模型總體決定系數(shù)大于LR模型。結(jié)果表明BN模型的預(yù)測(cè)效果優(yōu)于LR模型。說(shuō)明本試驗(yàn)構(gòu)建的貝葉斯網(wǎng)絡(luò)模型可信度較高。

    3 討論與結(jié)論

    保護(hù)地栽培條件下瓜類白粉病的發(fā)生受到多種因素的影響,筆者通過(guò)相關(guān)分析發(fā)現(xiàn),溫度濕度、光照度與病害發(fā)生相關(guān)。通過(guò)查閱相關(guān)文獻(xiàn)得知,分生孢子在10~30 ℃都能萌發(fā),而20~25 ℃最為適宜,超過(guò)30 ℃或者低于-1 ℃,則難以萌發(fā)并失去活力[32],病菌分生孢子萌發(fā)所要求的濕度范圍比較大,即發(fā)病適宜的相對(duì)濕度為45%~95%,但是濕度不能過(guò)高,因分生孢子吸水后膨壓過(guò)大,會(huì)引起破裂,反而對(duì)萌發(fā)不利[26]。設(shè)施環(huán)境屬于密閉環(huán)境,空氣流動(dòng)比較小,溫度和濕度很容易達(dá)到病害所需的條件,筆者分析得出,溫度符合病害發(fā)生概率為0.449,濕度符合病害發(fā)生的概率為0.730。病害發(fā)生的概率為0.695。這與高峰[25]和徐寧等[26]的研究病害發(fā)生規(guī)律一致。筆者發(fā)現(xiàn)光照度影響發(fā)生概率比較小,但是在設(shè)施環(huán)境下,對(duì)溫度的變化影響比較大。而對(duì)露地作物來(lái)說(shuō),光照同樣影響病害的發(fā)生[6]。為此,筆者將光照條件作為影響溫室黃瓜白粉病發(fā)生的因素之一。通過(guò)對(duì)日光溫室黃瓜白粉病的病害調(diào)查結(jié)果分析,白粉病發(fā)生前期病害傳播較慢,中期、后期傳播速度快、危害性較大,利用模型預(yù)測(cè)病害發(fā)生的準(zhǔn)確率為95%、92%、91%、87%,整體預(yù)測(cè)的準(zhǔn)確度為87%。表明模型在整體病害預(yù)測(cè)和各地區(qū)預(yù)測(cè)病害的效果比較好。

    筆者將貝葉斯網(wǎng)絡(luò)與黃瓜白粉病發(fā)生預(yù)測(cè)模型相結(jié)合,建立了基于貝葉斯網(wǎng)絡(luò)的日光溫室黃瓜白粉病預(yù)警系統(tǒng)。通過(guò)構(gòu)建黃瓜白粉病發(fā)生預(yù)測(cè)模型,探明了黃瓜白粉病發(fā)生概率與溫室內(nèi)主要環(huán)境因子之間的關(guān)系,對(duì)建立的模型進(jìn)行檢驗(yàn),進(jìn)一步證明模型檢驗(yàn)結(jié)果較優(yōu),表明系統(tǒng)模型在預(yù)測(cè)溫室黃瓜白粉病發(fā)生方面具有較強(qiáng)的普適性。此外,黃瓜白粉病的發(fā)生除與環(huán)境氣象因子有密切關(guān)系外,還與品種抗性、水肥管理、田間管理等因素密切相關(guān),今后應(yīng)在綜合考慮其他因子的情況下,進(jìn)一步完善預(yù)測(cè)模型,保持模型的滾動(dòng)發(fā)展,以達(dá)到預(yù)警效果最佳,為溫室大棚蔬菜病害防治提供參考。

    參考文獻(xiàn)

    [1] 張怡,徐克東,陳佩佩,等.河南省周口地區(qū)瓜類白粉病的分子鑒定與分析[J].華中師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,47(6):830-835.

    [2] 呂雄杰,王曉蓉,賈寶紅.基于物聯(lián)網(wǎng)技術(shù)的日光溫室黃瓜白粉病預(yù)警系統(tǒng)研究(英文)[J].Agricultural Science & Technology,2016,17(12):2873-2876.

    [3] 沈麗平,徐強(qiáng),MOUAMMAR A,等.黃瓜白粉病抗性遺傳模型分析[J].江蘇農(nóng)業(yè)學(xué)報(bào),2011,27(2):361-365.

    [4] 蔣龍燕.氣候變化與海南橡膠樹(shù)白粉病流行關(guān)系的分析[D].海口:海南大學(xué),2015.

    [5] TRECATE L,SEDLALOVA B,MIESLEROVA B,et al. Effect of temperature on infection and development of powdery mildew on cucumber[J]. Plant Pathology,2019,68(6):1165-1178.

    [6] WANG H,JIANG Y P,YU H J,et al.Light quality affects incidence of powdery mildew,expression of defence-related genes and associated metabolism in cucumber plants[J].European journal of Plant Pathology,2010,127(1):125-135.

    [7] WHIPPS J M,BUDGE S P.Effect of humidity on development of tomato powdery mildew(Oidium lycopersici)in the glasshouse[J].European Journal of Plant Pathology,2000,106(4):395-397.

    [8] 王國(guó)偉,劉嘉欣.基于卷積神經(jīng)網(wǎng)絡(luò)的玉米病害識(shí)別方法研究[J].中國(guó)農(nóng)機(jī)化學(xué)報(bào),2021,42(2):139-145.

    [9] 徐寧.塑料大棚黃瓜白粉病和霜霉病流行預(yù)測(cè)和管理系統(tǒng)的研究[D].南京:南京農(nóng)業(yè)大學(xué),2003.

    [10] KAMPMANN H H,HANSEN O B.Using colour image analysis for quantitative assessment of powdery mildew on cucumber[J].Euphytica,1994,79(12):19-27.

    [11] 葉小林.探析物聯(lián)網(wǎng)技術(shù)在桃樹(shù)栽培與病蟲(chóng)害防治中的應(yīng)用[J]. 農(nóng)業(yè)開(kāi)發(fā)與裝備,2021(3):208-209.

    [12] 李磊福,孫秋玉,史娟,等.寧夏溫室瓜菜白粉病菌鑒定及病害流行預(yù)測(cè)模型構(gòu)建[J].植物保護(hù)學(xué)報(bào), 2017,44(5):788-795.

    [13] 尹哲,原鍇,谷培云,等.京郊溫室大棚黃瓜白粉病流行預(yù)測(cè)模型構(gòu)建[J].中國(guó)植保導(dǎo)刊,2014,34(3):46-48.

    [14] 趙森棟,劉挺.因果關(guān)系及其在社會(huì)媒體上的應(yīng)用研究綜述[J].軟件學(xué)報(bào),2014,25(12):2733-2752.

    [15] 謝樂(lè).基于貝葉斯網(wǎng)絡(luò)的農(nóng)超對(duì)接供應(yīng)鏈風(fēng)險(xiǎn)預(yù)警模型研究[D].南昌:江西財(cái)經(jīng)大學(xué),2020.

    [16] 張妍,張秋英,李發(fā)東,等.基于穩(wěn)定同位素和貝葉斯模型的引黃灌區(qū)地下水硝酸鹽污染源解析[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)(中英文),2019,27(3):484-493.

    [17] 馮思宇.河北省冬小麥白粉病氣候預(yù)報(bào)的貝葉斯網(wǎng)絡(luò)模型研究[D].石家莊:河北師范大學(xué),2018.

    [18] 聶臣巍.基于貝葉斯網(wǎng)絡(luò)的小麥條銹病預(yù)測(cè)研究[J].安徽農(nóng)業(yè)科學(xué),2014,42(16):5027-5030.

    [19] JMP A,JAA B,JAG A,et al.Efficient and accurate structural fusion of Bayesian networks[J].Information Fusion,2021,66:155-169.

    [20] ZHU Y,LIU D,CHEN G,et al.Mathematical modeling for active and dynamic diagnosis of crop diseases based on Bayesian networks and incremental learning[J].Mathematical and Computer Modelling,2013,58(3/4):514-523.

    [21] SINGH N,GUPTA N.Bayesian network for decisionsupport on pest management of tomato fruit borer,H. armigera[J].International Journal of Engineering & Technology,2017,6(4):168-170.

    [22] LU W,NEWLANDS N K,CARISSE O,et al. Disease risk forecasting with Bayesian learning networks: Application to grape powdery mildew(Erysiphe necator)in vineyards[J].Agronomy,2020,10(5):622.

    [23] 杜玉寧,陳杭,白小軍,等.寧夏黃瓜新品種及主栽品種對(duì)白粉病的抗性評(píng)價(jià)[J].中國(guó)瓜菜,2019,32(11):27-30.

    [24] 劉鑒,任愛(ài)新,劉冉,等.考慮日光溫室空間異質(zhì)性的黃瓜葉片濕潤(rùn)時(shí)間估算模型研究[J].智慧農(nóng)業(yè)(中英文),2020,2(2):135-144.

    [25] 高鋒.吉林省主要瓜類白粉病發(fā)生規(guī)律及綜合防治措施的研究[D].長(zhǎng)春:吉林農(nóng)業(yè)大學(xué),2007.

    [26] 徐寧,劉鳳權(quán),胡白石,等.室內(nèi)模擬塑料大棚黃瓜白粉病流行預(yù)測(cè)模型的建立[J].植物保護(hù), 2004(6):52-55.

    [27] 李卓,金恭璽,郎寧,等.溫濕度與光照對(duì)葡萄霜霉病菌孢子囊萌發(fā)及存活的影響[J].新疆農(nóng)業(yè)科學(xué),2016,53(11):2090-2097.

    [28] 劉艷杰,李霞. 基于貝葉斯網(wǎng)絡(luò)的學(xué)生成績(jī)預(yù)測(cè)[J].山東理工大學(xué)學(xué)報(bào)(自然科學(xué)版),2019,33(5): 75-78.

    [29] 熊曉夏,陳龍,梁軍,等.基于貝葉斯網(wǎng)絡(luò)模型的道路交通事故鏈生成與演化研究[J].公路交通科技, 2018,35(5):99-107.

    [30] 黃影平.貝葉斯網(wǎng)絡(luò)發(fā)展及其應(yīng)用綜述[J].北京理工大學(xué)學(xué)報(bào),2013,33(12):1211-1219.

    [31] 黃建明.貝葉斯網(wǎng)絡(luò)在學(xué)生成績(jī)預(yù)測(cè)中的應(yīng)用[J].計(jì)算機(jī)科學(xué),2012,39(S3):280-282.

    [32] 范紹強(qiáng),謝咸升,李峰,等.山西省小麥條銹病流行趨勢(shì)預(yù)測(cè)模型研究[J].中國(guó)農(nóng)業(yè)生態(tài)學(xué)報(bào),2007,15(4):113-155.

    猜你喜歡
    貝葉斯網(wǎng)絡(luò)預(yù)測(cè)模型白粉病
    一到春季就流行 蔬菜白粉病該咋防
    無(wú)人機(jī)數(shù)據(jù)鏈測(cè)試與評(píng)估研究
    科技資訊(2016年25期)2016-12-27 16:22:32
    基于貝葉斯網(wǎng)絡(luò)的流域內(nèi)水文事件豐枯遭遇研究
    基于神經(jīng)網(wǎng)絡(luò)的北京市房?jī)r(jià)預(yù)測(cè)研究
    商情(2016年43期)2016-12-23 14:23:13
    中國(guó)石化J分公司油氣開(kāi)發(fā)投資分析與預(yù)測(cè)模型研究
    基于興趣預(yù)測(cè)和熱點(diǎn)分析的聯(lián)合推薦算法研究 
    拉薩設(shè)施月季白粉病的發(fā)生與防治
    西藏科技(2016年8期)2016-09-26 09:00:21
    基于IOWHA法的物流需求組合改善與預(yù)測(cè)模型構(gòu)建
    基于預(yù)測(cè)模型加擾動(dòng)控制的最大功率點(diǎn)跟蹤研究
    科技視界(2016年1期)2016-03-30 13:37:45
    基于貝葉斯網(wǎng)絡(luò)的城市居民出行方式研究
    最好的美女福利视频网| 亚洲男人天堂网一区| 一级a爱视频在线免费观看| 热re99久久国产66热| 免费在线观看日本一区| xxx96com| 欧美色欧美亚洲另类二区 | 99精品在免费线老司机午夜| 婷婷丁香在线五月| 两人在一起打扑克的视频| 国产亚洲av嫩草精品影院| 1024视频免费在线观看| 久久久久久国产a免费观看| 国产精品,欧美在线| 村上凉子中文字幕在线| 亚洲国产精品999在线| 久久人人精品亚洲av| 伊人久久大香线蕉亚洲五| 精品电影一区二区在线| 日本免费一区二区三区高清不卡 | 久久精品人人爽人人爽视色| 深夜精品福利| 欧美午夜高清在线| 国产高清激情床上av| 亚洲国产精品sss在线观看| 国产不卡一卡二| 巨乳人妻的诱惑在线观看| 精品国产一区二区三区四区第35| 嫩草影视91久久| 日韩中文字幕欧美一区二区| 激情在线观看视频在线高清| 成人国产综合亚洲| 国产精品99久久99久久久不卡| 女性被躁到高潮视频| 欧美绝顶高潮抽搐喷水| 丝袜美足系列| 亚洲午夜理论影院| 韩国av一区二区三区四区| 99久久99久久久精品蜜桃| 一级毛片精品| 午夜久久久久精精品| 精品第一国产精品| 久久久久久亚洲精品国产蜜桃av| 日本一区二区免费在线视频| 午夜精品久久久久久毛片777| 国产精品影院久久| 黄色 视频免费看| 久热爱精品视频在线9| 97超级碰碰碰精品色视频在线观看| 国产91精品成人一区二区三区| 中文字幕最新亚洲高清| 国产精品久久电影中文字幕| 国产亚洲欧美精品永久| 色综合欧美亚洲国产小说| 欧洲精品卡2卡3卡4卡5卡区| 国产精品亚洲av一区麻豆| 成年女人毛片免费观看观看9| 精品久久蜜臀av无| 国产精品久久电影中文字幕| 久久久久国内视频| 啪啪无遮挡十八禁网站| 久久久久久久久久久久大奶| 夜夜爽天天搞| 亚洲一码二码三码区别大吗| 亚洲黑人精品在线| 禁无遮挡网站| 午夜精品久久久久久毛片777| 操出白浆在线播放| 日本vs欧美在线观看视频| 女人精品久久久久毛片| 国产麻豆69| 久久久国产成人精品二区| 十分钟在线观看高清视频www| 欧美丝袜亚洲另类 | 在线天堂中文资源库| 深夜精品福利| 成人免费观看视频高清| 美女高潮到喷水免费观看| 欧美在线黄色| 日本 av在线| 亚洲,欧美精品.| 日本精品一区二区三区蜜桃| 国产精品九九99| 国产成人精品无人区| 欧美人与性动交α欧美精品济南到| 欧美黄色淫秽网站| 亚洲avbb在线观看| 老司机午夜福利在线观看视频| 欧美成人免费av一区二区三区| 欧美大码av| 精品久久久久久久人妻蜜臀av | 国产精品99久久99久久久不卡| 老司机在亚洲福利影院| 久久精品91蜜桃| 欧美一区二区精品小视频在线| 神马国产精品三级电影在线观看 | av片东京热男人的天堂| 精品国产一区二区久久| 欧美色欧美亚洲另类二区 | 免费在线观看影片大全网站| 日韩精品青青久久久久久| 一二三四在线观看免费中文在| 黄色成人免费大全| 国产成人免费无遮挡视频| 成人永久免费在线观看视频| 国产精品一区二区精品视频观看| 欧美亚洲日本最大视频资源| 久久人人97超碰香蕉20202| 啦啦啦免费观看视频1| 亚洲久久久国产精品| 亚洲人成77777在线视频| 1024视频免费在线观看| 国产色视频综合| 首页视频小说图片口味搜索| 国产av一区在线观看免费| 人人澡人人妻人| 韩国av一区二区三区四区| 亚洲av五月六月丁香网| 女性被躁到高潮视频| 黄网站色视频无遮挡免费观看| 久久久久久人人人人人| 一级a爱片免费观看的视频| 人人妻人人爽人人添夜夜欢视频| 极品教师在线免费播放| 中国美女看黄片| 国产极品粉嫩免费观看在线| √禁漫天堂资源中文www| 精品国产乱子伦一区二区三区| 搡老熟女国产l中国老女人| 日韩中文字幕欧美一区二区| 亚洲国产毛片av蜜桃av| 亚洲久久久国产精品| 亚洲欧美一区二区三区黑人| 99在线视频只有这里精品首页| 免费不卡黄色视频| 精品午夜福利视频在线观看一区| 久久人妻av系列| 美国免费a级毛片| 久久精品91蜜桃| 日韩精品免费视频一区二区三区| 亚洲无线在线观看| a级毛片在线看网站| 免费在线观看亚洲国产| 黄片播放在线免费| 欧美激情久久久久久爽电影 | www.www免费av| 十分钟在线观看高清视频www| 怎么达到女性高潮| 亚洲欧美精品综合一区二区三区| 婷婷精品国产亚洲av在线| 宅男免费午夜| 大陆偷拍与自拍| 大香蕉久久成人网| 成熟少妇高潮喷水视频| 成人欧美大片| 亚洲欧美精品综合一区二区三区| 男人舔女人的私密视频| 中文字幕人成人乱码亚洲影| 国产精品日韩av在线免费观看 | 亚洲精品国产一区二区精华液| 久久国产精品人妻蜜桃| 看片在线看免费视频| 在线十欧美十亚洲十日本专区| 日本免费一区二区三区高清不卡 | 色综合欧美亚洲国产小说| 中文字幕高清在线视频| 中文亚洲av片在线观看爽| 国产精品影院久久| 亚洲激情在线av| 一级片免费观看大全| 久久久精品国产亚洲av高清涩受| 亚洲天堂国产精品一区在线| 超碰成人久久| 精品一区二区三区视频在线观看免费| 亚洲国产精品999在线| 久久久水蜜桃国产精品网| 色播亚洲综合网| 亚洲人成伊人成综合网2020| 正在播放国产对白刺激| 欧美日本视频| av免费在线观看网站| 99国产极品粉嫩在线观看| 国产精品久久久久久精品电影 | 亚洲熟妇熟女久久| 亚洲五月天丁香| 非洲黑人性xxxx精品又粗又长| 久久伊人香网站| 中文字幕高清在线视频| 国产午夜福利久久久久久| 国产又爽黄色视频| 成人免费观看视频高清| 国产精品 国内视频| 国内精品久久久久久久电影| 亚洲情色 制服丝袜| 国产精品一区二区精品视频观看| 男女午夜视频在线观看| 亚洲精品国产精品久久久不卡| 母亲3免费完整高清在线观看| 性欧美人与动物交配| 日韩大尺度精品在线看网址 | 嫩草影院精品99| 可以在线观看的亚洲视频| 国产av一区在线观看免费| 少妇熟女aⅴ在线视频| 涩涩av久久男人的天堂| 国产午夜精品久久久久久| 色精品久久人妻99蜜桃| 久久久国产成人精品二区| 丁香六月欧美| 亚洲三区欧美一区| 一a级毛片在线观看| 欧美日韩亚洲综合一区二区三区_| 无人区码免费观看不卡| 精品一区二区三区av网在线观看| 欧美国产精品va在线观看不卡| 老熟妇仑乱视频hdxx| 国产欧美日韩精品亚洲av| 大香蕉久久成人网| 色av中文字幕| 欧美绝顶高潮抽搐喷水| www.自偷自拍.com| 天堂√8在线中文| 国产亚洲av嫩草精品影院| 熟女少妇亚洲综合色aaa.| 天天一区二区日本电影三级 | 香蕉国产在线看| 多毛熟女@视频| 熟妇人妻久久中文字幕3abv| 此物有八面人人有两片| 天堂√8在线中文| 欧美人与性动交α欧美精品济南到| 国产在线精品亚洲第一网站| 无限看片的www在线观看| 日韩欧美国产一区二区入口| 深夜精品福利| 欧美成人一区二区免费高清观看 | 91麻豆av在线| 99国产精品一区二区三区| 久久欧美精品欧美久久欧美| 日韩三级视频一区二区三区| 后天国语完整版免费观看| xxx96com| 国产视频一区二区在线看| 亚洲国产看品久久| 男女之事视频高清在线观看| 国产区一区二久久| 亚洲中文日韩欧美视频| 97碰自拍视频| 国产午夜精品久久久久久| 日韩有码中文字幕| 看免费av毛片| 午夜视频精品福利| 老汉色av国产亚洲站长工具| 一本久久中文字幕| 亚洲在线自拍视频| 免费在线观看影片大全网站| 国产欧美日韩精品亚洲av| 如日韩欧美国产精品一区二区三区| 午夜成年电影在线免费观看| 又紧又爽又黄一区二区| 亚洲精品久久成人aⅴ小说| 69av精品久久久久久| 在线av久久热| 熟妇人妻久久中文字幕3abv| www.自偷自拍.com| 亚洲,欧美精品.| av超薄肉色丝袜交足视频| 色尼玛亚洲综合影院| 精品不卡国产一区二区三区| 一区二区日韩欧美中文字幕| 母亲3免费完整高清在线观看| 午夜福利成人在线免费观看| 757午夜福利合集在线观看| 国产成人系列免费观看| 国内毛片毛片毛片毛片毛片| 亚洲成a人片在线一区二区| 欧美日本视频| 精品人妻1区二区| 久热爱精品视频在线9| 天天躁狠狠躁夜夜躁狠狠躁| 国内毛片毛片毛片毛片毛片| 亚洲一区高清亚洲精品| 欧美精品啪啪一区二区三区| 国产又爽黄色视频| 亚洲黑人精品在线| 成年女人毛片免费观看观看9| 成人永久免费在线观看视频| 亚洲全国av大片| 国产黄a三级三级三级人| 日本一区二区免费在线视频| 国产又色又爽无遮挡免费看| 深夜精品福利| 高清毛片免费观看视频网站| 国产日韩一区二区三区精品不卡| 精品一区二区三区视频在线观看免费| 99国产精品免费福利视频| 一区在线观看完整版| 国产成人免费无遮挡视频| 嫩草影视91久久| 国产av精品麻豆| 精品日产1卡2卡| 丝袜美足系列| 久久婷婷成人综合色麻豆| 亚洲中文av在线| 欧美黄色淫秽网站| 熟妇人妻久久中文字幕3abv| 免费少妇av软件| 欧美在线一区亚洲| 日韩欧美国产在线观看| 在线观看www视频免费| 亚洲精品美女久久久久99蜜臀| 黄色片一级片一级黄色片| 精品卡一卡二卡四卡免费| 波多野结衣高清无吗| 淫秽高清视频在线观看| 99久久精品国产亚洲精品| 免费在线观看视频国产中文字幕亚洲| 久久国产乱子伦精品免费另类| 国产一区二区激情短视频| 很黄的视频免费| 欧美精品啪啪一区二区三区| 男女做爰动态图高潮gif福利片 | 久久久久国产一级毛片高清牌| 欧美激情久久久久久爽电影 | 国产激情久久老熟女| 亚洲中文字幕日韩| 国产三级黄色录像| 欧美黄色淫秽网站| 成人国语在线视频| 女人爽到高潮嗷嗷叫在线视频| 久久久国产欧美日韩av| 高清黄色对白视频在线免费看| 伦理电影免费视频| 伊人久久大香线蕉亚洲五| 免费看美女性在线毛片视频| 少妇被粗大的猛进出69影院| 中国美女看黄片| 国产精品野战在线观看| 国产成人系列免费观看| 又大又爽又粗| 精品久久久久久久人妻蜜臀av | 欧美激情高清一区二区三区| 日韩欧美在线二视频| 欧美一区二区精品小视频在线| 黄色视频,在线免费观看| 色尼玛亚洲综合影院| 一区二区三区高清视频在线| 欧美激情高清一区二区三区| 亚洲男人的天堂狠狠| 精品不卡国产一区二区三区| 啦啦啦 在线观看视频| 午夜福利成人在线免费观看| 妹子高潮喷水视频| 大陆偷拍与自拍| 狂野欧美激情性xxxx| 精品少妇一区二区三区视频日本电影| 男女床上黄色一级片免费看| e午夜精品久久久久久久| 手机成人av网站| 中文字幕人妻丝袜一区二区| 很黄的视频免费| 不卡av一区二区三区| 国产91精品成人一区二区三区| 成人手机av| 国产精品一区二区免费欧美| 此物有八面人人有两片| 免费看美女性在线毛片视频| 一卡2卡三卡四卡精品乱码亚洲| 亚洲熟女毛片儿| 天天添夜夜摸| 精品午夜福利视频在线观看一区| 一二三四在线观看免费中文在| 国产一区二区三区综合在线观看| 97超级碰碰碰精品色视频在线观看| 激情在线观看视频在线高清| 精品乱码久久久久久99久播| 国产成人av教育| 熟妇人妻久久中文字幕3abv| 久久精品国产清高在天天线| 青草久久国产| 亚洲 欧美一区二区三区| 免费看a级黄色片| 高清黄色对白视频在线免费看| 9色porny在线观看| 亚洲国产毛片av蜜桃av| 亚洲五月婷婷丁香| 亚洲午夜理论影院| 国产欧美日韩一区二区三| 国产麻豆成人av免费视频| 日韩精品中文字幕看吧| 日本精品一区二区三区蜜桃| 亚洲成国产人片在线观看| 亚洲国产欧美日韩在线播放| 51午夜福利影视在线观看| 人妻丰满熟妇av一区二区三区| 搡老岳熟女国产| 国产精品美女特级片免费视频播放器 | 一区二区三区高清视频在线| 精品免费久久久久久久清纯| 最近最新中文字幕大全电影3 | 超碰成人久久| 18禁国产床啪视频网站| 亚洲在线自拍视频| 日韩中文字幕欧美一区二区| 性色av乱码一区二区三区2| 丁香六月欧美| 国产高清激情床上av| 成人国语在线视频| 99久久久亚洲精品蜜臀av| 亚洲精品中文字幕一二三四区| 9色porny在线观看| 成人三级黄色视频| 国产精华一区二区三区| 久久精品国产综合久久久| av欧美777| 免费在线观看黄色视频的| 亚洲一区中文字幕在线| 天堂动漫精品| 亚洲av电影不卡..在线观看| 丝袜在线中文字幕| 日本a在线网址| 麻豆一二三区av精品| 国产欧美日韩一区二区三| 97超级碰碰碰精品色视频在线观看| 亚洲中文字幕日韩| 91av网站免费观看| 午夜久久久久精精品| 在线视频色国产色| 亚洲色图 男人天堂 中文字幕| 欧美 亚洲 国产 日韩一| 长腿黑丝高跟| 日韩欧美在线二视频| av中文乱码字幕在线| 亚洲人成电影观看| 国产伦一二天堂av在线观看| 两个人视频免费观看高清| 欧美色欧美亚洲另类二区 | 满18在线观看网站| 亚洲色图综合在线观看| 一进一出好大好爽视频| 国产精品亚洲av一区麻豆| 亚洲成人久久性| 色播亚洲综合网| 在线观看舔阴道视频| 长腿黑丝高跟| 丰满人妻熟妇乱又伦精品不卡| 91成人精品电影| 一级作爱视频免费观看| 亚洲欧美精品综合一区二区三区| 亚洲九九香蕉| 国产精品电影一区二区三区| 一区福利在线观看| e午夜精品久久久久久久| 老司机靠b影院| 在线观看午夜福利视频| 欧美精品啪啪一区二区三区| 国产精品日韩av在线免费观看 | 欧美在线黄色| av福利片在线| 国产精品亚洲美女久久久| 久久中文看片网| 国产97色在线日韩免费| 不卡av一区二区三区| 亚洲熟妇中文字幕五十中出| 亚洲精品国产一区二区精华液| av天堂在线播放| 欧美人与性动交α欧美精品济南到| 一边摸一边做爽爽视频免费| а√天堂www在线а√下载| 免费一级毛片在线播放高清视频 | 丰满人妻熟妇乱又伦精品不卡| 国产成人系列免费观看| 69精品国产乱码久久久| 女同久久另类99精品国产91| а√天堂www在线а√下载| 麻豆一二三区av精品| 中文字幕人妻丝袜一区二区| 色尼玛亚洲综合影院| 欧美在线黄色| 桃色一区二区三区在线观看| 久久人人97超碰香蕉20202| 亚洲成av人片免费观看| 亚洲av五月六月丁香网| 高清黄色对白视频在线免费看| 精品卡一卡二卡四卡免费| 99精品欧美一区二区三区四区| 免费久久久久久久精品成人欧美视频| 无遮挡黄片免费观看| 国产极品粉嫩免费观看在线| 999精品在线视频| 一本久久中文字幕| 97人妻天天添夜夜摸| 两个人视频免费观看高清| 欧美乱码精品一区二区三区| 国产亚洲欧美在线一区二区| 国产在线精品亚洲第一网站| 国产私拍福利视频在线观看| 国产精品久久电影中文字幕| 国产在线观看jvid| 老司机午夜十八禁免费视频| 嫩草影视91久久| 无人区码免费观看不卡| 少妇粗大呻吟视频| 青草久久国产| 亚洲三区欧美一区| 99riav亚洲国产免费| 在线观看66精品国产| 亚洲一码二码三码区别大吗| 国产伦一二天堂av在线观看| 国产一区二区三区在线臀色熟女| 麻豆一二三区av精品| 啦啦啦 在线观看视频| 啦啦啦观看免费观看视频高清 | 校园春色视频在线观看| 美女午夜性视频免费| 亚洲第一欧美日韩一区二区三区| 午夜福利18| 操美女的视频在线观看| 自线自在国产av| 可以在线观看的亚洲视频| 亚洲一区二区三区色噜噜| 一区二区三区精品91| 欧美成人午夜精品| 丝袜人妻中文字幕| 麻豆成人av在线观看| 亚洲精品久久国产高清桃花| 精品少妇一区二区三区视频日本电影| 亚洲国产欧美一区二区综合| 亚洲色图av天堂| 欧美成人一区二区免费高清观看 | 日本五十路高清| 97超级碰碰碰精品色视频在线观看| 午夜免费成人在线视频| x7x7x7水蜜桃| 嫩草影视91久久| 免费不卡黄色视频| 亚洲美女黄片视频| 少妇的丰满在线观看| 巨乳人妻的诱惑在线观看| 可以免费在线观看a视频的电影网站| 午夜福利18| 性欧美人与动物交配| 天天躁狠狠躁夜夜躁狠狠躁| 久久精品国产99精品国产亚洲性色 | 久久精品91蜜桃| 亚洲色图 男人天堂 中文字幕| 悠悠久久av| 国产精品秋霞免费鲁丝片| 亚洲精品在线观看二区| 非洲黑人性xxxx精品又粗又长| 精品国内亚洲2022精品成人| 亚洲久久久国产精品| 黄色a级毛片大全视频| 长腿黑丝高跟| 免费一级毛片在线播放高清视频 | 麻豆一二三区av精品| 91成年电影在线观看| 欧美乱码精品一区二区三区| 91av网站免费观看| 又黄又粗又硬又大视频| 18禁美女被吸乳视频| 少妇熟女aⅴ在线视频| 一个人观看的视频www高清免费观看 | 日韩有码中文字幕| 亚洲精品在线观看二区| 两人在一起打扑克的视频| 亚洲 国产 在线| 国产成人一区二区三区免费视频网站| 97碰自拍视频| 非洲黑人性xxxx精品又粗又长| 精品人妻1区二区| 久久久久国产一级毛片高清牌| 搡老岳熟女国产| 亚洲精品国产一区二区精华液| 亚洲视频免费观看视频| 91大片在线观看| 深夜精品福利| 精品国内亚洲2022精品成人| 国产精品日韩av在线免费观看 | 亚洲av日韩精品久久久久久密| 国产亚洲精品av在线| 久久久精品欧美日韩精品| 久久精品国产综合久久久| 搡老妇女老女人老熟妇| 免费无遮挡裸体视频| 一级,二级,三级黄色视频| 亚洲熟女毛片儿| 日本一区二区免费在线视频| 神马国产精品三级电影在线观看 | 久久人人精品亚洲av| 日本撒尿小便嘘嘘汇集6| 午夜精品国产一区二区电影| 久久久国产成人免费| 欧美亚洲日本最大视频资源| 久久人妻熟女aⅴ| 91在线观看av| 成人亚洲精品一区在线观看| 欧美大码av| 日韩精品免费视频一区二区三区| 如日韩欧美国产精品一区二区三区| 久久精品亚洲精品国产色婷小说| 国产欧美日韩一区二区三| 97超级碰碰碰精品色视频在线观看| 国产真人三级小视频在线观看| 黄色丝袜av网址大全| 午夜福利,免费看| 国产1区2区3区精品| 叶爱在线成人免费视频播放| 搞女人的毛片| 亚洲aⅴ乱码一区二区在线播放 | 亚洲av电影在线进入| 不卡av一区二区三区| 99国产精品一区二区蜜桃av| 国产精品免费视频内射| 1024视频免费在线观看|