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

    大興安嶺不同坡位地表可燃物含水率的動態(tài)變化與建模

    2017-09-27 09:16:52宋雨胡海清孫龍蘇漳文
    森林工程 2017年5期
    關(guān)鍵詞:陰坡陽坡白樺林

    宋雨,胡海清,孫龍,蘇漳文

    (東北林業(yè)大學(xué) 林學(xué)院,哈爾濱 150040)

    大興安嶺不同坡位地表可燃物含水率的動態(tài)變化與建模

    宋雨,胡海清*,孫龍,蘇漳文

    (東北林業(yè)大學(xué) 林學(xué)院,哈爾濱 150040)

    大興安嶺林區(qū)是我國森林資源最為豐富、森林火災(zāi)最為嚴(yán)重的地區(qū)。分析該地區(qū)氣象、地表可燃物含水率動態(tài)變化規(guī)律,準(zhǔn)確估測地表可燃物含水率可以降低森林火災(zāi)發(fā)生的幾率,對該地區(qū)的森林防火工作意義重大。本研究以大興安嶺地區(qū)三種典型林型(陰坡白樺林、陽坡樟子松林、陰坡落葉松林)的地表可燃物為研究對象,分析不同林型、不同坡位(上坡位、中坡位、下坡位)地表可燃物含水率和氣象參數(shù)的動態(tài)變化規(guī)律,計(jì)算氣象因子與地表可燃物含水率的Pearson相關(guān)系數(shù),分別采用Nelson法、Simard法和氣象要素回歸法建立地表可燃物含水率預(yù)測模型。研究結(jié)論如下:對于研究區(qū)不同林分,陽坡樟子松林空氣平均溫度高于陰坡白樺林,陰坡白樺林平均溫度高于陰坡落葉松林;陰坡落葉松空氣濕度高于陰坡白樺林空氣濕度,陰坡白樺林空氣濕度高于陽坡樟子松林空氣濕度;地表可燃物含水率的關(guān)系為陰坡落葉松林>陰坡白樺林>陽坡樟子松林;三種林分地表可燃物含水率均與空氣相對濕度、當(dāng)日最高空氣相對濕度呈正相關(guān),與空氣溫度、當(dāng)日最高空氣溫度呈負(fù)相關(guān),與風(fēng)速的相關(guān)性較低;Nelson法建立的模型決定系數(shù)R2為0.635~0.934,Simard法構(gòu)建的預(yù)測模型,決定系數(shù)R2為0.82~0.907,逐步回歸模型的決定系數(shù)R2為0.191~0.408,Nelson法和Simard法是更適于該地區(qū)的地表可燃物含水率預(yù)測方法。

    可燃物含水率;氣象因子;地形因子;回歸模型

    0 引言

    森林火災(zāi)是一種具有突發(fā)性和破壞性的森林燃燒現(xiàn)象,對森林結(jié)構(gòu)和森林環(huán)境的影響是深刻的。森林火災(zāi)的種類繁多,在各類森林火災(zāi)中,地表火發(fā)生的概率最高;而地表可燃物是地表火的主要載體,其含水率的大小是影響地表火發(fā)生概率大小及火行為指標(biāo)的直接因素[1-2]。對地表可燃物含水率動態(tài)變化規(guī)律進(jìn)行研究和預(yù)測是做好火險(xiǎn)預(yù)測工作的關(guān)鍵[3]。因此,通過建立適用于該地區(qū)的地表可燃物含水率模型并通過對其含水率進(jìn)行預(yù)測也日漸成為林火相關(guān)研究的重點(diǎn)。

    地表可燃物含水率受多種因素的綜合影響,包括可燃物自身的理化性質(zhì)、所處的地形以及周圍氣象條件、土壤中的水分等[3-4]。氣象因素是影響可燃物含水率動態(tài)變化的因素之一,其中空氣相對濕度和空氣溫度是主要的影響因子[5];降水量和降水時間作為重要的林火管理要素分別影響了森林枯落物層對降水的截留量和死可燃物含水率變化[6-7]。此外,地形、土壤中的水分也間接影響可燃物含水率的變化;地形差異會引起空氣濕度和空氣溫度等生態(tài)因子的重新分配,森林作為一種特殊的下墊面,與地形變化結(jié)合形成不同的局部氣候[8]。當(dāng)土壤中有足夠的水分,并且地表可燃物沒有達(dá)到飽和點(diǎn)時,地表可燃物就會從土壤上升的水汽中吸收水分[9]。

    國內(nèi)外學(xué)者對可燃物含水率與氣象因子關(guān)系及其預(yù)測進(jìn)行了大量的研究[3,10-15]。加拿大的Wright(1928年)提出了用相對濕度對林火發(fā)生進(jìn)行預(yù)報(bào)的方法,以相對濕度50%為界限,認(rèn)為小于50%時林火容易發(fā)生;美國的Gisborne提出了多因子預(yù)報(bào)的方法。Byram等人采用氣象因子對可燃物含水率進(jìn)行了估算。Gisborne和Jemison等大量學(xué)者基于降水量,以積累系數(shù)表示降水量,并用“火險(xiǎn)尺”模型來預(yù)測可燃物含水率。1963年,由Albina提出時滯的概念,他認(rèn)為可燃物干燥過程為一個指數(shù)衰減的過程并提供了干燥反應(yīng)方程。到了80年代,美國和加拿大基于氣象要素分別研制出了“火行為預(yù)報(bào)和可燃物模型系統(tǒng)(BEHAVE)”以及加拿大林火氣象指數(shù)系統(tǒng)(Canadian fire weather index,F(xiàn)WI),并依據(jù)可燃物平衡含水率法對細(xì)小可燃物含水率進(jìn)行了預(yù)測。

    目前林火管理的需求明顯提高,僅僅研究可燃物含水率及其影響因子的動態(tài)變化已經(jīng)無法適應(yīng)人們生產(chǎn)生活的要求。隨著數(shù)學(xué)模型的廣泛應(yīng)用,關(guān)于可燃物含水率模型的研究也越來越多。近年來我國在可燃物含水率預(yù)測模型領(lǐng)域發(fā)展較快,建立了一系列不同林型、不同地區(qū)的可燃物含水率預(yù)測模型[16-21]。

    依據(jù)當(dāng)前研究進(jìn)展,預(yù)測可燃物含水率的方法主要有以下4種:①平衡含水率法,②氣象要素回歸法,③遙感估測法,④過程模型法。其中平衡含水率法和氣象要素回歸法是當(dāng)前應(yīng)用比較廣泛的可燃物含水率預(yù)測方法[22]。物理上的可靠性使得平衡含水率法在小尺度的研究中較為精確,但對于大尺度則會大大增加工作量且精確度降低[23-24]。Nelson模型和Simard模型是當(dāng)前使用平衡含水率法對可燃物含水率進(jìn)行預(yù)測時,預(yù)測效果最好的兩種模型[25-26]。氣象要素回歸采用統(tǒng)計(jì)方法,相對簡單,使用廣泛[16,27-28]。

    黑龍江省大興安嶺地區(qū)是我國面積最大、火災(zāi)最嚴(yán)重的國有林區(qū)[29],因此建立適用于該地區(qū)的含水率預(yù)測模型,加強(qiáng)該地區(qū)可燃物含水率預(yù)測的準(zhǔn)確性,對降低森林火災(zāi)發(fā)生的可能性具有重大意義。本文通過實(shí)地監(jiān)測,以地表可燃物為研究對象,分析了研究區(qū)內(nèi)可燃物含水率及氣象因子的動態(tài)變化規(guī)律,并采用不同建模方法對研究區(qū)域不同林型地表可燃物含水率進(jìn)行建模,選出最優(yōu)建模方法,對提高該地區(qū)地表可燃物含水率預(yù)測的精度,進(jìn)而完善林火管理方法具有一定的意義。

    1 研究地區(qū)

    本文選取的研究地區(qū)為位于黑龍江省大興安嶺北部的西林吉林業(yè)局,東經(jīng)121°07′至124°20′,北緯52°10′至53°33′。該地區(qū)屬寒溫帶季風(fēng)性大陸氣候,年平均氣溫為-4.7℃,年均降水量為460.8 mm,植物生長期大約為100 d,無霜期85~115 d左右。冬季由于受到大陸季風(fēng)和蒙古高壓的控制,寒冷異常,氣溫最低可達(dá)-52℃,凍期漫長,可達(dá)7個月之久;夏季受到海洋氣團(tuán)的影響,濕潤且雨水充足,極端條件下最高氣溫達(dá)32℃。春季風(fēng)大干旱,是火災(zāi)頻發(fā)期。西林吉林業(yè)局森林總面積為16 281 km2,森林覆蓋率為90.6%。其中有林地面積覆蓋率為90.41%,灌木林面積覆蓋率為0.23%。林木蓄積量14 647萬m3。有林地的樹種組成結(jié)構(gòu)較為單一,主要為興安落葉松(Larixgmelinii)、樟子松(Pinussylvestrisvar.mongolica)和白樺(Betulaplatyphylla),其他樹種所占比重較少。

    圖1 研究區(qū)域位置圖Fig.1 Site map of research region

    2 研究方法

    2.1 數(shù)據(jù)獲取

    于2016年春季防火期(5~6月),在西林吉林業(yè)局選擇3種林分布置樣地,分別為陰坡白樺林、陽坡樟子松林和陰坡落葉松林(樣地信息見表1)。3種林分內(nèi)樣地大小均為40 m×80 m。在每個樣地的上坡位、中坡位和下坡位平行設(shè)置五個采樣點(diǎn),各個采樣點(diǎn)之間的間距為10 m。3塊樣地共設(shè)置45個采樣點(diǎn)。在每個采樣點(diǎn)挖取邊長為20 cm 的正方形樣方,將采集的地表可燃物裝入30 cm×20 cm的尼龍網(wǎng)中,尼龍網(wǎng)密度為40目。將尼龍網(wǎng)放回采樣點(diǎn),保持尼龍網(wǎng)與周圍環(huán)境充分接觸,即可以維持地表凋落物與周圍環(huán)境自然的水分交換狀態(tài),也避免了外界因素造成的樣品的增加或減少。

    表1 樣地信息Tab.1 Information of sample plots

    8點(diǎn)至16點(diǎn),每隔2 h 對采樣點(diǎn)的尼龍網(wǎng)進(jìn)行稱重,記錄樣品鮮重。監(jiān)測期內(nèi)的氣象數(shù)據(jù)獲取方式為手持氣象站手動測量,包括空氣溫度,空氣相對濕度、風(fēng)速和降雨量。監(jiān)測期結(jié)束后,將樣品帶回實(shí)驗(yàn)室烘干,并稱量樣品干重。

    3 數(shù)據(jù)處理

    本文分別采用時滯-平衡含水率法和氣象要素回歸法進(jìn)行建模。依據(jù)Catchpole等人提出的基于平衡含水率和氣象動態(tài)數(shù)據(jù)的直接估計(jì)法,首先分別采用Nelson模型法和Simard模型法[25]估算出可燃物的平衡含水率,然后將細(xì)小可燃物微分方程離散化,最后對Catchpole的直接估計(jì)法進(jìn)行參數(shù)估計(jì),得到可燃物含水率預(yù)測模型。

    Catchpole的直接估計(jì)法是基于細(xì)小可燃物微分方程(1)來對時滯平衡含水率模型進(jìn)行參數(shù)估計(jì)[13,30-32]。

    dm/dt=-(M-E)/τ。

    (1)

    式中:M為可燃物含水率,%;t為時間,h;E為平衡含水率,%;τ為時滯,h。

    Nelson法(1984)平衡含水率與氣象因子的函數(shù)關(guān)系:

    (2)

    式中:R為普適氣體常量(取8.314J/K·mol);T為環(huán)境空氣溫度,K;H為空氣相對濕度,%;m為H2O相對分子質(zhì)量(取18g/mol);α、β為待估計(jì)參數(shù)。

    Simard法平衡含水率與氣象因子的函數(shù)關(guān)系:

    (3)

    式中:T為環(huán)境空氣溫度,K;H為空氣相對濕度,%。

    將公式(1)離散化并分別把方程(2)、公式(3)帶入,即可得到直接估計(jì)方程。

    M(ti)=λ2Mi-1+λ(1-λ)Ei-1+(1-λ)Ei。

    (4)

    式中:Ei為ti時刻的平衡含水率值,%;λ=exp(-δt/(2τ)),時滯τ=-δt/(2lnλ)。

    氣象要素回歸法的形式為多元線性方程,并采用逐步回歸法進(jìn)行參數(shù)估計(jì),選取相關(guān)性最好的方程[33]。

    (5)

    式中:M為可燃物含水率,%;Xi為所選用的氣象,包括前n日和當(dāng)日14:00時的空氣相對濕度,%、空氣溫度,℃、風(fēng)速,m/s、降雨量,mm等;bi為待估計(jì)參數(shù)。

    4 結(jié)果與分析

    4.1 氣象因子動態(tài)分析

    根據(jù)手持氣象站測量的監(jiān)測點(diǎn)氣象數(shù)據(jù),得到研究區(qū)每日14時空氣溫度動態(tài)變化圖(圖2)。由結(jié)果可知,春季防火期5月20日至6月25日不同林型,坡位空氣溫度變化趨勢相近??諝鉁囟瓤傮w呈上升趨勢,隨著降雨的發(fā)生而波動。陰坡白樺林和陰坡落葉松林空氣溫度波動較大,陽坡樟子松林空氣溫度波動較小。降雨過后發(fā)生森林火災(zāi)的可能性較小,故沒有監(jiān)測降雨后的氣象參數(shù),在圖中以缺失數(shù)據(jù)斷線顯示。陰坡白樺林空氣溫度變化范圍為10.8~31.6℃,平均溫度約為21.43℃;陽坡樟子松林空氣溫度變化范圍為9.5~33.8℃,平均溫度約為22.81℃;陰坡落葉松林空氣溫度變化范圍為7.4~32.4℃,平均溫度約為19.83℃??傮w上,陽坡溫度高于陰坡,陰坡白樺林溫度高于陰坡落葉松林。大興安嶺地區(qū)春季防火期林內(nèi)積雪已完全融化,地表可燃物裸露于空氣中,地表可燃物水分蒸發(fā)快,含水率較低,易發(fā)生森林火災(zāi)。

    對空氣濕度數(shù)據(jù)進(jìn)行處理,得到研究區(qū)每日14時空氣濕度動態(tài)變化圖(圖3)。由結(jié)果可知,不同林型,不同坡位空氣濕度變化趨勢相近;陰坡白樺林空氣濕度的變化范圍為17.5%~87%,平均濕度約為38.46%;陽坡樟子松林空氣濕度變化范圍為12.2%~95.4%,平均濕度約為36.84%;陰坡落葉松林空氣濕度變化范圍為17.6%~98.2%,平均濕度約為39.3%。三種林型的空氣濕度陰坡落葉松林>陰坡白樺林>陽坡樟子松林。在春季防火期,研究區(qū)域內(nèi)降雨較為頻繁,可燃物含水率波動較大。一般而言,當(dāng)空氣相對濕度低于30%時,易發(fā)生重大森林火災(zāi)。

    (a)陰坡白樺林空氣溫度動態(tài)變化圖

    (b)陽坡樟子松林空氣溫度動態(tài)變化圖

    (c)陰坡落葉松林空氣溫度動態(tài)變化圖圖2 研究區(qū)每日14時空氣溫度動態(tài)變化圖Fig.2 The dynamic change of air temperature at 14:00 daily in the research region

    (a)陰坡白樺林空氣濕度動態(tài)變化圖

    (b)陽坡樟子松林空氣濕度動態(tài)變化圖

    (c)陰坡落葉松林空氣濕度動態(tài)變化圖圖3 研究區(qū)每日14時空氣濕度動態(tài)變化圖Fig.3 The dynamic change of air humidity at 14:00 daily in the research region

    4.2 不同坡位可燃物含水率動態(tài)變化

    圖4給出了監(jiān)測期內(nèi)每日14時可燃物含水率的動態(tài)變化情況。由結(jié)果可知,不同林型,不同坡位可燃物含水率變化趨勢相近。陰坡白樺林上坡位、中坡位和下坡位可燃物含水率平均值分別為76.13%、96.13%和128%,下坡位>中坡位>上坡位,陰坡白樺林可燃物含水率平均值為100.4%;陽坡樟子松林上坡位、中坡位和下坡位可燃物含水率平均值分別為48.16%、49.12%和57.07%,下坡位>中坡位>上坡位,陽坡樟子松林可燃物含水率平均值為51.45%;陰坡落葉松林上坡位、中坡位和下坡位可燃物含水率平均值分別為150.04%、141.01%和170.19%,下坡位>上坡位>中坡位,陰坡落葉松林可燃物含水率平均值為153.75%。三種林型可燃物含水率的大小關(guān)系為:陰坡落葉松林>陰坡白樺林>陽坡樟子松林。

    (a)陰坡白樺林可燃物含水率動態(tài)變化圖

    (b)陽坡樟子松林可燃物含水率動態(tài)變化圖

    (c)陰坡落葉松林可燃物含水率動態(tài)變化圖圖4 監(jiān)測期內(nèi)每日14時可燃物含水率的動態(tài)變化Fig.4 The dynamic change of the moisture content of fuel at 14:00 during monitoring period

    與陰坡白樺林相比,陰坡落葉松林郁閉度更高,水分蒸發(fā)較慢且落葉松林地表可燃物多由落葉松枯枝落葉組成,與陰坡白樺林相比,持水能力較強(qiáng),故陰坡落葉松林可燃物含水率高于陰坡白樺林可燃物含水率。陽坡樟子松日照時間長,太陽輻射量較大,水分蒸發(fā)較快,故陽坡樟子松林可燃物含水率低于陰坡白樺林和陰坡落葉松林可燃物含水率。

    4.3 可燃物含水率與氣象因子的關(guān)系

    表2給出了三種林分地表凋落物含水率與空氣相對濕度(H)、空氣溫度(T)、當(dāng)日最低空氣相對濕度(H1)、當(dāng)日最高空氣溫度(T1)和風(fēng)速(W)的Pearson相關(guān)系數(shù)。結(jié)果表明,三種林分地被凋落物含水率均與空氣相對濕度、當(dāng)日最高空氣相對濕度呈正相關(guān),與空氣溫度、當(dāng)日最高空氣溫度呈負(fù)相關(guān),與風(fēng)速的相關(guān)性較低。

    表2 Pearson相關(guān)系數(shù)Tab.2 Pearson correlation coefficient

    注:**表示在0.01水平上顯著相關(guān)。

    4.4 可燃物含水率模型的構(gòu)建

    表3給出了通過Nelson法和Simard法構(gòu)建的不同林型的可燃物含水率預(yù)測模型的參數(shù)。對于Nelson模型,參數(shù)中的β代表了平衡含水率對空氣溫度和空氣相對濕度的響應(yīng)速度,參數(shù)β的絕對值與可燃物的持水性大小成反比,陰坡落葉松林可燃物的持水性強(qiáng)于陽坡樟子松林,陽坡樟子松林可燃物持水性強(qiáng)于陰坡白樺林;時滯τ的最小值為15.3,最大值為28.5;通過Nelson法構(gòu)建的預(yù)測模型決定系數(shù)R2為0.635~0.934。對于Simard法構(gòu)建的預(yù)測模型,時滯最小值為38.5,時滯最大值為42.9,決定系數(shù)R2為0.82~0.907。使用Nelson法和Simard法算得的可燃物時滯差異較大。

    表4為采用逐步回歸法建立的可燃物含水率氣象回歸模型的參數(shù)。從表中可以得出,逐步回歸模型的決定系數(shù)R2為0.191~0.408。采用時滯平衡含水率法建立的模型相較于采用逐步回歸法建立的氣象回歸模型,更適于估計(jì)該研究區(qū)的可燃物含水率。

    表3 不同林型Nelson和Simard模型參數(shù)Tab.3 Parameters of Nelson and Simard models for different forest types

    表4 不同林型氣象要素回歸模型參數(shù)Tab.4 Parameters of meteorological element regression model of different forest types

    5 結(jié)論與討論

    (1)對于空氣溫度和空氣濕度的動態(tài)變化。陽坡樟子松林空氣平均溫度高于陰坡白樺林,陰坡白樺林平均溫度高于陰坡落葉松林??諝鉁囟扰c林分郁閉度成反比。陰坡落葉松林空氣濕度最高,陰坡白樺松林空氣濕度居中,陽坡樟子松林空氣濕度最低。

    (2)對于地表可燃物含水率的動態(tài)變化。2016年春季防火期,地表可燃物含水率的關(guān)系為陰坡落葉松林>陰坡白樺林>陽坡樟子松林。受太陽輻射影響,陽坡樟子松林內(nèi)空氣溫度較高,空氣濕度較低,林內(nèi)地表可燃物含水率最低;陰坡落葉松林地表可燃物主要有針葉構(gòu)成,結(jié)構(gòu)質(zhì)密,持水性較強(qiáng),且陰坡落葉松林郁閉度最高,發(fā)生降雨后水分揮發(fā)較慢,故陰坡落葉松林地表可燃物含水率最高。對于不同坡位地表可燃物含水率的關(guān)系,陰坡白樺林和陽坡樟子松林均為下坡位>中坡位>上坡位,陰坡落葉松林為下坡位>上坡位>中坡位。

    (3)對于地表可燃物含水率與氣象因子的關(guān)系與建模。研究區(qū)的三種林分地表凋落物含水率與空氣相對濕度、當(dāng)日最高空氣相對濕度呈正相關(guān),與空氣溫度、當(dāng)日最高空氣溫度和風(fēng)速呈負(fù)相關(guān),與風(fēng)速的相關(guān)性較低。采用Nelson法和Simard法建立的含水率預(yù)測模型適用性要強(qiáng)于采用氣象要素回歸法建立的可燃物含水率預(yù)測模型。

    [1] 舒立福,田曉瑞.國外森林防火工作現(xiàn)狀及展望[J].世界林業(yè)研究,1997,10(2):29-37.

    [2] 于宏洲,金森,邸雪穎.以小時為步長的大興安嶺興安落葉松林地表可燃物含水率預(yù)測模型[J].應(yīng)用生態(tài)學(xué)報(bào),2013,24(6):1565-1571.

    [3] 胡海清.林火生態(tài)與管理[M].北京:中國林業(yè)出版社,2005.

    [4] Viney N R.A Review of Fine Fuel Moisture Modelling[J].International Journal of Wildland Fire,1991,1(4):215-234.

    [5] 胡海清,陸昕,孫龍,等.大興安嶺典型林分地表死可燃物含水率動態(tài)變化及預(yù)測模型[J].應(yīng)用生態(tài)學(xué)報(bào),2016,27(7):2212-2224.

    [6] Walsh R P D,Voigt P J.Vegetation litter:an underestimated variable in hydrology and geomorphology[J].Journal of Biogeography,1977,4(3):253-274.

    [7] Fosberg M A.Theory of precipitation effects on dead cylindrical fuels[J].Forest Science,1972,18(2):98-108.

    [8] Holden Z A,Jolly W M.Modeling topographic influences on fuel moisture and fire danger in complex terrain to improve wildland fire management decision support[J].Forest Ecology&Management,2011,262(12):2133-2141.

    [9] Schaap M G,Bouten W.Forest floor evaporation in a dense Douglas fir stand[J].Journal of Hydrology,1997,193(1):97-113.

    [10] 白尚斌,張曉麗.林火預(yù)測預(yù)報(bào)研究綜述[J].森林防火,2008,26(2):22-25.

    [11] 舒立福,田曉瑞.國外森林防火工作現(xiàn)狀及展望[J].世界林業(yè)研究,1997,10(2):29-37.

    [12] Wotton B M,Beverly J L.Stand-specific litter moisture content calibrations for the canadian fine fuel moisture code[J].International Journal of Wildland Fire,2007,6(4):463-472.

    [13] Byram G,Nelson R M.An analysis of the drying process in forest fuel material[J].General Technical Report - Southern Research Station,USDA Forest Service,2015,19(45):1-56.

    [14] Forest Service United States Department of Agriculture.Modeling moisture content of fine dead wildland fuels:input to the behave fire prediction system[J].Usda Forest Service Intermountain Research Station Research Paper,1986,11(359):1-61.

    [15] 程藝龍,孫龍.東北三種典型針葉樹種可燃物烘干參數(shù)的確定及影響因素分析[J].森林工程,2016,32(4);22-27.

    [16] 何忠秋,李長勝,張成鋼,等.森林可燃物含水量模型的研究[J].森林防火,1995,13(2):15-16+20.

    [17] 張國防,林文革,花昆福,等.杉木人工林地表易燃物含水率變化規(guī)律[J].森林與環(huán)境學(xué)報(bào),2000,20(1):77-79.

    [18] 覃先林,張子輝,易浩若,等.一種預(yù)測森林可燃物含水率的方法[J].火災(zāi)科學(xué),2001,10(3):159-162.

    [19] 薛家翠,望勝玲,曾祥福,等.鄂西林地可燃物含水率及火險(xiǎn)等級的氣象預(yù)報(bào)研究[J].華中農(nóng)業(yè)大學(xué)學(xué)報(bào),2006,25(6):679-682.

    [20] 葉更新,葉希瑩.長白山區(qū)林下可燃物含水率的觀測特征和相關(guān)分析[J].中國農(nóng)學(xué)通報(bào),2011,27(22):22-27.

    [21] 李海洋,胡海清,孫龍.我國森林死可燃物含水率與氣象和土壤因子關(guān)系模型研究[J].森林工程,2016,32(3);1-6.

    [22] 劉昕,邸雪穎.三種方法對森林地表可燃物含水率的預(yù)測評價[J].森林工程,2013,29(2):8-13.

    [23] Jin S,Li X Y,Li Y X.Desorption processes of four fine fuels[J].Journal of Northeast Forestry University,2000,28(1):35-38.

    [24] Catchpole E A,Catchpole W R,Viney N R,et al.Estimating fuel response time and predicting fuel moisture content from field data[J].International Journal of Wildland Fire,2001,10(2):215-222.

    [25] Nelson M R.A method for describing equilibrium moisture content of forest fuels[J].Canadian Journal of Forest Research,2011,14(4):597-600.

    [26] 劉曦,金森.平衡含水率法預(yù)測死可燃物含水率的研究進(jìn)展[J].林業(yè)科學(xué),2007,43(12):126-133.

    [27] Saglam B,Bilgili E,Kucuk O,et al.Determination of surface fuels moisture contents based on weather conditions[J].Forest Ecology&Management,2006,234:75.

    [28] Viegas D X,Pinol J,Viegas M T,et al.Estimating live fine fuels moisture content using meteorologically-based indices[J].International Journal of Wildland Fire,2001,10(2):223-240.

    [29] 于宏洲,金森,邸雪穎.以小時為步長的大興安嶺興安落葉松林地表可燃物含水率預(yù)測模型[J].應(yīng)用生態(tài)學(xué)報(bào),2013,24(6):1565-1571.

    [30] Slijepcevic A,Anderson W.Hourly variation in fine fuel moisture in eucalypt forests in Tasmania[J].Forest Ecology&Management,2006,234:36.

    [31] Simard A J.The Moisture Content of Forest Fuels 1:a Review of the Basic Concepts[M].Canada:Forest Fire Research Institute,1968.

    [32] 姜宏志.楓樺紅松林地被死可燃物載量估測模型研究[J].林業(yè)科技,2016,41(6);21-23.

    [33] 居恩德,陳貴榮,王瑞君.可燃物含水率與氣象要素相關(guān)性的研究[J].森林防火,1993,11(1):17-21.

    DynamicChangeandModelingofMoistureContentofSurfaceFuelinDifferentSlopePositionsofDaxing’anling

    Song Yu,Hu Haiqing*,Sun Long,Su Zhangwen

    (College of Forestry,Northeast Forestry University,Harbin 150040)

    Daxing’anling is the forest which has the most abundant forest resources and has the most serious forest fires in China.It is of great significance for the forest fire prevention work in this area to analyze the dynamic change of meteorological elements and surface fuel moisture content,and to accurately estimate the moisture content of surface fuel moisture content to reduce the occurrence of forest fires.In this study,the surface fuel of three kinds of typical forest types(north slopeBetulaplatyphylla,south slopePinussylvestrisLinn,north slopeLarixgmelinii)in Daxing’anling was used as the research object.The dynamic changes of surface fuel moisture content and meteorological parameters in different forest types and different slope positions(high,mid and low slope-position)were analyzed.The Person correlation coefficient between meteorological factors and surface fuel moisture content was calculated.The Nelson method,Simard method and meteorological element regression method were selected to establish the fuel moisture content prediction models.The conclusion of the study were as follows:for different forest types in the study area,the average temperature of south slope Pinus sylvestris Linn was higher than that of north slopeBetulaplatyphyllaand the average temperature of north slope Betula platyphylla was higer than that of north slopeLarixgmelinii.The air humidity of north slope Betula platyphylla was higher than that of north slopeLarixgmeliniiand the air humidity of north slopeLarixgmeliniiwas higher than that of south slopePinussylvestrisLinn.The relationship between surface fuel moisture content was north slopeLarixgmelinii>north slopeBetulaplatyphylla>south slopePinussylvestrisLinn.The surface fuel moisture content in three types was positively correlated with air relative humidity and the highest daily air relative humidity,and negatively correlated with air temperature and the highest daily air temperature,but had low correlation with wind speed.The determination coefficient R2of stepwise regression model is 0.191~0.408.The determination coefficient R2of the model which is established with the Simard method is 0.82~0.907.The determination coefficient R2of stepwise regression model is 0.191~0.408.The Nelson and Simard methods were more suitable for the prediction of fuel moisture content in the study area.

    Fuel moisture content;meteorological factors;topographic factors;regression model

    S 762

    :A

    :1001-005X(2017)05-0001-07

    2017-04-16

    林業(yè)公益性行業(yè)科研專項(xiàng)(20144402)

    宋雨,碩士研究生。研究方向:森林防火。

    胡海清,教授。研究方向:森林防火。E-mail:517422138@qq.com

    宋雨,胡海清,孫龍,等.大興安嶺不同坡位地表可燃物含水率的動態(tài)變化與建模[J].森林工程,2017,33(5):1-7.

    猜你喜歡
    陰坡陽坡白樺林
    阿爾泰山天然林不同采伐坡向的植物群落特征分析
    山坡羊
    那一片白樺林
    草原歌聲(2021年4期)2021-06-10 02:25:48
    大地掠美,穿越白樺林
    文苑(2019年22期)2019-12-07 05:29:12
    接壩山區(qū)陡坡地段不同坡向油松造林生長狀況分析
    陽光灑滿白樺林
    北極光(2018年5期)2018-11-19 02:23:58
    不同坡向 坡位對祁連圓柏生長的影響
    冀西北陽坡硅藻土礦地質(zhì)特征及應(yīng)用前景
    有關(guān)林線高低問題
    憂傷的白樺林
    纯流量卡能插随身wifi吗| av不卡在线播放| 国产成人aa在线观看| 婷婷成人精品国产| 亚洲精品美女久久久久99蜜臀 | 亚洲av综合色区一区| 国产日韩欧美在线精品| 熟女电影av网| 久久精品国产自在天天线| 丰满迷人的少妇在线观看| 男人爽女人下面视频在线观看| 十八禁高潮呻吟视频| 九草在线视频观看| av免费观看日本| 制服人妻中文乱码| 欧美人与性动交α欧美软件| 黄色 视频免费看| 18禁动态无遮挡网站| 亚洲精品一二三| 一级毛片黄色毛片免费观看视频| 日韩制服骚丝袜av| 久久午夜福利片| 女人被躁到高潮嗷嗷叫费观| 亚洲在久久综合| 妹子高潮喷水视频| 男人舔女人的私密视频| 国产又爽黄色视频| 美女中出高潮动态图| 日韩在线高清观看一区二区三区| 国产成人精品福利久久| 精品国产乱码久久久久久男人| 肉色欧美久久久久久久蜜桃| 日韩熟女老妇一区二区性免费视频| 日韩视频在线欧美| 精品人妻在线不人妻| 热99久久久久精品小说推荐| 精品少妇黑人巨大在线播放| 精品国产一区二区三区久久久樱花| 肉色欧美久久久久久久蜜桃| 日韩一本色道免费dvd| 成人手机av| 欧美bdsm另类| 18禁裸乳无遮挡动漫免费视频| 男女边吃奶边做爰视频| 久久国产亚洲av麻豆专区| 在线 av 中文字幕| 国产精品久久久久久精品古装| 国产精品一二三区在线看| 亚洲精品中文字幕在线视频| 麻豆精品久久久久久蜜桃| av视频免费观看在线观看| 成人毛片a级毛片在线播放| 国产精品人妻久久久影院| videossex国产| 亚洲美女搞黄在线观看| 看十八女毛片水多多多| 中国三级夫妇交换| a 毛片基地| 日韩精品免费视频一区二区三区| 成人毛片a级毛片在线播放| 一边亲一边摸免费视频| a 毛片基地| 两个人免费观看高清视频| 成人毛片a级毛片在线播放| 国产精品免费大片| 丰满饥渴人妻一区二区三| 日韩在线高清观看一区二区三区| 最黄视频免费看| 一区二区三区四区激情视频| 国产爽快片一区二区三区| 精品酒店卫生间| 中文字幕人妻熟女乱码| 午夜日韩欧美国产| 国产色婷婷99| 免费高清在线观看视频在线观看| 日日摸夜夜添夜夜爱| 欧美 日韩 精品 国产| 久久av网站| 两性夫妻黄色片| 久久99热这里只频精品6学生| 如何舔出高潮| 夫妻性生交免费视频一级片| 爱豆传媒免费全集在线观看| 欧美成人午夜精品| 九草在线视频观看| 亚洲精品乱久久久久久| 久久久久国产精品人妻一区二区| 深夜精品福利| 少妇人妻 视频| 久久午夜福利片| 日本91视频免费播放| 欧美激情高清一区二区三区 | 色视频在线一区二区三区| 香蕉精品网在线| 国产欧美亚洲国产| 十分钟在线观看高清视频www| 午夜福利影视在线免费观看| 久久狼人影院| 91aial.com中文字幕在线观看| 人人妻人人添人人爽欧美一区卜| 久久久久久人人人人人| 欧美激情高清一区二区三区 | 在线 av 中文字幕| 尾随美女入室| 亚洲欧美成人综合另类久久久| 精品第一国产精品| 黑人巨大精品欧美一区二区蜜桃| 亚洲欧美精品自产自拍| 两个人看的免费小视频| 亚洲精品第二区| 国产伦理片在线播放av一区| 亚洲美女视频黄频| 亚洲精品日本国产第一区| 久久这里有精品视频免费| 2022亚洲国产成人精品| 国产精品久久久久久av不卡| 国产成人精品一,二区| 国产成人精品无人区| 大片电影免费在线观看免费| 九色亚洲精品在线播放| 成年人免费黄色播放视频| 制服诱惑二区| 天美传媒精品一区二区| 久久久a久久爽久久v久久| 久久精品国产亚洲av高清一级| 免费久久久久久久精品成人欧美视频| 999精品在线视频| 成人免费观看视频高清| 亚洲一级一片aⅴ在线观看| 国产片特级美女逼逼视频| 人妻 亚洲 视频| 免费人妻精品一区二区三区视频| av国产久精品久网站免费入址| 老女人水多毛片| 香蕉丝袜av| 国产精品免费大片| 精品99又大又爽又粗少妇毛片| 亚洲成人av在线免费| 久久久国产一区二区| 在线免费观看不下载黄p国产| 伦精品一区二区三区| 精品国产乱码久久久久久小说| 日韩欧美精品免费久久| 久久久欧美国产精品| 亚洲熟女精品中文字幕| 亚洲四区av| 国产av一区二区精品久久| 国产黄频视频在线观看| 日韩欧美精品免费久久| 少妇被粗大的猛进出69影院| 久久久久精品久久久久真实原创| 中文乱码字字幕精品一区二区三区| 国产精品香港三级国产av潘金莲 | 高清av免费在线| 国产 精品1| 人妻一区二区av| 国产极品粉嫩免费观看在线| 777久久人妻少妇嫩草av网站| 亚洲男人天堂网一区| 日韩精品有码人妻一区| 国产 精品1| 亚洲成av片中文字幕在线观看 | 国产精品国产三级国产专区5o| 看免费成人av毛片| 日本猛色少妇xxxxx猛交久久| 久久久久久免费高清国产稀缺| 亚洲av电影在线进入| av又黄又爽大尺度在线免费看| 亚洲欧美清纯卡通| 国产伦理片在线播放av一区| 777米奇影视久久| 少妇人妻精品综合一区二区| 波野结衣二区三区在线| 男女国产视频网站| 三上悠亚av全集在线观看| 老司机影院毛片| 国产成人精品婷婷| 另类精品久久| 丝袜美腿诱惑在线| 青青草视频在线视频观看| 十八禁高潮呻吟视频| 日韩一本色道免费dvd| 亚洲色图综合在线观看| 久久婷婷青草| av片东京热男人的天堂| 久久久久久免费高清国产稀缺| 老司机影院毛片| 国产精品欧美亚洲77777| 蜜桃国产av成人99| 狠狠婷婷综合久久久久久88av| 制服人妻中文乱码| 精品人妻一区二区三区麻豆| av在线观看视频网站免费| 国产1区2区3区精品| 国产成人精品一,二区| 只有这里有精品99| 亚洲一区中文字幕在线| 十八禁网站网址无遮挡| 国产亚洲欧美精品永久| www.自偷自拍.com| 超色免费av| 久久精品国产综合久久久| 久久婷婷青草| 亚洲欧洲国产日韩| 人妻系列 视频| 色94色欧美一区二区| 午夜久久久在线观看| 国产一级毛片在线| 麻豆av在线久日| 日日摸夜夜添夜夜爱| 男女无遮挡免费网站观看| 波多野结衣一区麻豆| 亚洲av免费高清在线观看| 国产毛片在线视频| 亚洲内射少妇av| 欧美 日韩 精品 国产| 大话2 男鬼变身卡| 天美传媒精品一区二区| 久热这里只有精品99| 五月伊人婷婷丁香| 少妇人妻精品综合一区二区| 久久国产精品大桥未久av| 亚洲欧洲日产国产| 丁香六月天网| 国产熟女午夜一区二区三区| 最近的中文字幕免费完整| 九草在线视频观看| 国产精品欧美亚洲77777| 韩国精品一区二区三区| 街头女战士在线观看网站| 亚洲av免费高清在线观看| 亚洲精品美女久久av网站| 精品午夜福利在线看| 亚洲国产欧美日韩在线播放| av.在线天堂| 黄片播放在线免费| 在线观看一区二区三区激情| 成人国产麻豆网| 午夜激情av网站| 亚洲经典国产精华液单| 午夜免费鲁丝| 赤兔流量卡办理| 极品人妻少妇av视频| 日韩中文字幕欧美一区二区 | 国产一区二区三区综合在线观看| 97在线视频观看| 亚洲av在线观看美女高潮| 18在线观看网站| av电影中文网址| 久久久久精品人妻al黑| 国产片内射在线| 国产野战对白在线观看| 亚洲色图 男人天堂 中文字幕| 一级毛片我不卡| 亚洲精品美女久久久久99蜜臀 | 少妇人妻久久综合中文| 国产精品一二三区在线看| 国产免费福利视频在线观看| 精品久久久精品久久久| 成人免费观看视频高清| av国产精品久久久久影院| 在线观看免费视频网站a站| 午夜av观看不卡| 五月天丁香电影| 成人免费观看视频高清| 亚洲国产精品一区三区| 日韩视频在线欧美| 久久久久视频综合| www.自偷自拍.com| 日韩伦理黄色片| 秋霞在线观看毛片| 天天操日日干夜夜撸| 午夜福利视频精品| 不卡视频在线观看欧美| 夫妻午夜视频| 国产色婷婷99| 美女国产视频在线观看| 女的被弄到高潮叫床怎么办| 大香蕉久久成人网| 最近2019中文字幕mv第一页| 亚洲欧美精品自产自拍| 亚洲欧美清纯卡通| 国产av精品麻豆| 午夜免费鲁丝| 欧美人与善性xxx| 国产精品 欧美亚洲| 男女午夜视频在线观看| 人人妻人人澡人人看| 麻豆乱淫一区二区| 一级片免费观看大全| 久久精品国产综合久久久| 国产一区二区三区综合在线观看| 国产人伦9x9x在线观看 | 黑人欧美特级aaaaaa片| 日韩在线高清观看一区二区三区| 午夜激情久久久久久久| 少妇精品久久久久久久| 亚洲美女视频黄频| 亚洲美女搞黄在线观看| 一区在线观看完整版| 在线免费观看不下载黄p国产| 大话2 男鬼变身卡| 日韩欧美精品免费久久| 国产高清不卡午夜福利| av国产久精品久网站免费入址| 十八禁网站网址无遮挡| 国产高清国产精品国产三级| 精品人妻在线不人妻| 精品久久蜜臀av无| 精品亚洲成a人片在线观看| 亚洲天堂av无毛| 一边摸一边做爽爽视频免费| 日韩一本色道免费dvd| 在现免费观看毛片| 精品国产露脸久久av麻豆| 久久 成人 亚洲| av片东京热男人的天堂| 天天躁夜夜躁狠狠躁躁| 又大又黄又爽视频免费| 天天躁夜夜躁狠狠躁躁| 成人国产麻豆网| 黄色 视频免费看| 1024视频免费在线观看| 9热在线视频观看99| 中文字幕人妻丝袜制服| 美女福利国产在线| 秋霞伦理黄片| 日韩制服骚丝袜av| 伊人久久国产一区二区| 欧美日韩成人在线一区二区| 热re99久久国产66热| 黄片小视频在线播放| 丰满乱子伦码专区| 考比视频在线观看| 中文字幕最新亚洲高清| 尾随美女入室| 在线观看美女被高潮喷水网站| kizo精华| 久久久久人妻精品一区果冻| 日韩视频在线欧美| 男人操女人黄网站| av女优亚洲男人天堂| 亚洲成色77777| 成人国语在线视频| 777久久人妻少妇嫩草av网站| 大陆偷拍与自拍| 国产视频首页在线观看| videos熟女内射| 如日韩欧美国产精品一区二区三区| 夜夜骑夜夜射夜夜干| 国产成人午夜福利电影在线观看| 国产一区二区在线观看av| 建设人人有责人人尽责人人享有的| 美女脱内裤让男人舔精品视频| av一本久久久久| 十八禁高潮呻吟视频| 如日韩欧美国产精品一区二区三区| 美女国产高潮福利片在线看| 精品少妇黑人巨大在线播放| 免费久久久久久久精品成人欧美视频| 久久久久久久大尺度免费视频| 亚洲欧美精品自产自拍| 国产免费视频播放在线视频| 成人午夜精彩视频在线观看| 日韩制服骚丝袜av| 热re99久久精品国产66热6| 天堂中文最新版在线下载| 啦啦啦中文免费视频观看日本| 亚洲国产毛片av蜜桃av| 纵有疾风起免费观看全集完整版| 久久99热这里只频精品6学生| 精品视频人人做人人爽| 日韩制服骚丝袜av| 曰老女人黄片| 国产欧美亚洲国产| 熟女少妇亚洲综合色aaa.| 亚洲精品在线美女| 高清黄色对白视频在线免费看| 国产成人av激情在线播放| 老汉色∧v一级毛片| 国产福利在线免费观看视频| 男女高潮啪啪啪动态图| 男女午夜视频在线观看| 日本av免费视频播放| 久久国内精品自在自线图片| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产精品欧美亚洲77777| 18禁观看日本| 这个男人来自地球电影免费观看 | 一级黄片播放器| 你懂的网址亚洲精品在线观看| 亚洲精品久久午夜乱码| 国产精品久久久av美女十八| 国产国语露脸激情在线看| 少妇熟女欧美另类| 国产av码专区亚洲av| 一区二区三区四区激情视频| 精品人妻在线不人妻| 女人精品久久久久毛片| 啦啦啦啦在线视频资源| 最近最新中文字幕大全免费视频 | 色播在线永久视频| 高清欧美精品videossex| 天堂俺去俺来也www色官网| 精品福利永久在线观看| 大片免费播放器 马上看| 日韩,欧美,国产一区二区三区| 国产成人免费观看mmmm| 飞空精品影院首页| 免费人妻精品一区二区三区视频| 国产精品熟女久久久久浪| 亚洲成色77777| 成人漫画全彩无遮挡| 日韩电影二区| 亚洲av男天堂| 久久久久久人妻| 国产精品一国产av| 三级国产精品片| 电影成人av| 人妻人人澡人人爽人人| 亚洲中文av在线| 观看美女的网站| 热99久久久久精品小说推荐| 一边亲一边摸免费视频| 波多野结衣av一区二区av| 国产国语露脸激情在线看| 另类精品久久| 久久99蜜桃精品久久| 久久精品夜色国产| 狠狠婷婷综合久久久久久88av| 男男h啪啪无遮挡| 久久精品国产亚洲av涩爱| √禁漫天堂资源中文www| 亚洲精品久久成人aⅴ小说| 大片电影免费在线观看免费| 黑人巨大精品欧美一区二区蜜桃| 亚洲色图 男人天堂 中文字幕| 99热全是精品| 男女啪啪激烈高潮av片| 国产 精品1| 国产福利在线免费观看视频| 免费在线观看视频国产中文字幕亚洲 | 亚洲av电影在线观看一区二区三区| 青春草视频在线免费观看| 亚洲精品乱久久久久久| 日日啪夜夜爽| 久久久久精品久久久久真实原创| 最近2019中文字幕mv第一页| 最近最新中文字幕大全免费视频 | 亚洲av中文av极速乱| 亚洲欧美中文字幕日韩二区| 久久久国产精品麻豆| 亚洲色图综合在线观看| 免费大片黄手机在线观看| 国产人伦9x9x在线观看 | 日韩不卡一区二区三区视频在线| 久久韩国三级中文字幕| 精品一品国产午夜福利视频| 在线观看三级黄色| 久久久久久久亚洲中文字幕| 人人妻人人澡人人看| 久久久久久人人人人人| av免费观看日本| 女的被弄到高潮叫床怎么办| 五月伊人婷婷丁香| 如何舔出高潮| 欧美日韩精品成人综合77777| 99热全是精品| 91午夜精品亚洲一区二区三区| 国产精品嫩草影院av在线观看| 老汉色av国产亚洲站长工具| 欧美日韩视频高清一区二区三区二| 宅男免费午夜| 国产精品一区二区在线不卡| 成人午夜精彩视频在线观看| 激情五月婷婷亚洲| 久久国产精品男人的天堂亚洲| 亚洲欧洲国产日韩| 三上悠亚av全集在线观看| 成人手机av| 日韩一区二区三区影片| 汤姆久久久久久久影院中文字幕| 成年人午夜在线观看视频| 免费观看在线日韩| videossex国产| 国产免费现黄频在线看| 一区二区三区四区激情视频| 亚洲国产精品999| 一级毛片电影观看| 午夜福利乱码中文字幕| 亚洲精品国产一区二区精华液| 国产精品久久久av美女十八| 日日撸夜夜添| www.自偷自拍.com| 国产极品粉嫩免费观看在线| 午夜精品国产一区二区电影| 黄片无遮挡物在线观看| 一二三四在线观看免费中文在| 男女下面插进去视频免费观看| 五月天丁香电影| 欧美xxⅹ黑人| 超碰成人久久| 在线天堂最新版资源| 人妻人人澡人人爽人人| 最近中文字幕2019免费版| 丰满少妇做爰视频| 成人亚洲欧美一区二区av| 99久久综合免费| 成人亚洲精品一区在线观看| 91精品国产国语对白视频| 久久久久久久久久久免费av| 91精品伊人久久大香线蕉| 免费大片黄手机在线观看| 国产1区2区3区精品| 日韩熟女老妇一区二区性免费视频| 久久国产亚洲av麻豆专区| 美女高潮到喷水免费观看| 最近中文字幕高清免费大全6| 精品酒店卫生间| 美女国产高潮福利片在线看| 美女主播在线视频| 亚洲三级黄色毛片| 亚洲色图综合在线观看| 永久免费av网站大全| 伦理电影免费视频| 亚洲国产日韩一区二区| 女性被躁到高潮视频| 日本av免费视频播放| av在线老鸭窝| 国产高清国产精品国产三级| 欧美日韩综合久久久久久| 精品少妇一区二区三区视频日本电影 | 色94色欧美一区二区| 欧美精品亚洲一区二区| 免费播放大片免费观看视频在线观看| 日韩一区二区视频免费看| 色视频在线一区二区三区| 黄频高清免费视频| 国产亚洲一区二区精品| 18禁裸乳无遮挡动漫免费视频| 免费在线观看完整版高清| 亚洲国产最新在线播放| 18在线观看网站| 男女午夜视频在线观看| 9色porny在线观看| 波野结衣二区三区在线| 99热全是精品| 最近最新中文字幕免费大全7| 18禁国产床啪视频网站| 咕卡用的链子| 黄片无遮挡物在线观看| 国产成人精品一,二区| 看免费av毛片| 亚洲欧美一区二区三区黑人 | 99久久中文字幕三级久久日本| 欧美激情 高清一区二区三区| 亚洲综合精品二区| 欧美精品一区二区免费开放| av在线老鸭窝| 久久99热这里只频精品6学生| 亚洲av日韩在线播放| 国产精品蜜桃在线观看| 亚洲欧洲日产国产| 亚洲av电影在线进入| 国产乱来视频区| 美女国产视频在线观看| 欧美人与善性xxx| 五月天丁香电影| 日韩伦理黄色片| 考比视频在线观看| 咕卡用的链子| 如日韩欧美国产精品一区二区三区| 国产又爽黄色视频| 91精品国产国语对白视频| 亚洲精品视频女| 老女人水多毛片| 日本爱情动作片www.在线观看| 亚洲一码二码三码区别大吗| 飞空精品影院首页| 黑人猛操日本美女一级片| 大陆偷拍与自拍| 另类精品久久| 免费久久久久久久精品成人欧美视频| 色94色欧美一区二区| 国产精品国产三级国产专区5o| 亚洲精品乱久久久久久| 国产欧美亚洲国产| 侵犯人妻中文字幕一二三四区| 色吧在线观看| 91精品三级在线观看| 在线 av 中文字幕| 搡女人真爽免费视频火全软件| 亚洲人成网站在线观看播放| 亚洲综合精品二区| 五月伊人婷婷丁香| 午夜福利一区二区在线看| 好男人视频免费观看在线| 中文欧美无线码| 成年女人在线观看亚洲视频| 国产免费一区二区三区四区乱码| 国产精品国产三级专区第一集| 在线观看一区二区三区激情| 国产1区2区3区精品| 汤姆久久久久久久影院中文字幕| 午夜影院在线不卡| 黄色毛片三级朝国网站| 18在线观看网站| 日本欧美视频一区| 丝袜美足系列| 如日韩欧美国产精品一区二区三区| 久久婷婷青草| av电影中文网址| 少妇被粗大猛烈的视频| 日韩在线高清观看一区二区三区| 热99国产精品久久久久久7|