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

    地形環(huán)境因素對(duì)塔河林業(yè)局人為森林火災(zāi)發(fā)生的影響

    2015-12-22 12:33:46蘇漳文宋禹輝郭福濤楊婷婷王文輝劉愛琴福建農(nóng)林大學(xué)福州350002
    火災(zāi)科學(xué) 2015年1期
    關(guān)鍵詞:火點(diǎn)塔河林火

    蘇漳文,宋禹輝,郭福濤,楊婷婷,王文輝,劉愛琴(福建農(nóng)林大學(xué),福州,350002)

    ?

    地形環(huán)境因素對(duì)塔河林業(yè)局人為森林火災(zāi)發(fā)生的影響

    蘇漳文,宋禹輝,郭福濤*,楊婷婷,王文輝,劉愛琴*
    (福建農(nóng)林大學(xué),福州,350002)

    摘要:應(yīng)用1974年—2009年間人為火的空間地理坐標(biāo),結(jié)合研究地的基礎(chǔ)地理信息及矢量化林相圖,應(yīng)用Arc-GIS10.0中的空間分析工具和SPSS19.0的邏輯斯蒂回歸模型深入分析地形環(huán)境因素對(duì)人為火發(fā)生的影響。研究結(jié)果表明“到鐵路距離”對(duì)人為火的影響顯著,其次,在1∶1和1∶3選取隨機(jī)點(diǎn)條件下人為火發(fā)生分別與“優(yōu)勢(shì)樹種”和“林型”顯著相關(guān)。本文研究顯示在運(yùn)用“隨機(jī)點(diǎn)”方法進(jìn)行相關(guān)性分析時(shí),隨機(jī)點(diǎn)選取數(shù)量對(duì)模型擬合結(jié)果有一定影響。

    關(guān)鍵詞:北方針葉林;塔河地區(qū);人為火;地形環(huán)境;邏輯斯蒂回歸

    0 引言

    森林火災(zāi)(包括雷擊火和人為火)是北方針葉林中十分重要的生態(tài)干擾因子,影響森林更新和演替[1]。此外,林火還會(huì)直接導(dǎo)致森林資源和人民生命財(cái)產(chǎn)損失,也會(huì)大大增加森林經(jīng)營管理成本[2]。我國每年發(fā)生大約10,000次火災(zāi),年均過火面積為820,000公頃[3]。在森林火災(zāi)中,人為火(由于人類活動(dòng)或人為基礎(chǔ)設(shè)施而引起的火災(zāi))占很大比重[4,5]。位于北美阿拉斯加的北方針葉林帶,人為火比重也很大,年均發(fā)生次數(shù)約占到總森林火災(zāi)60%以上,過火面積由于人為撲救措施和技術(shù)手段的進(jìn)步,近年來呈遞減趨勢(shì)從38%(1950年—1959年)到5%(2000年—2005年),但火災(zāi)次數(shù)卻一直居高不下[6]。我國南方地區(qū)人為火比重也非常大,個(gè)別省份也可達(dá)到90%以上[7]。但由于南方全年濕度和降水量較大,人為火發(fā)生后如及時(shí)進(jìn)行撲救,不容易造成大面積森林損失。而塔河位于我國大興安嶺地區(qū),其植被和氣候特征與加拿大及美國北方針葉林帶十分相似,林分樹種組成多為針葉林,可燃性較強(qiáng)。秋季區(qū)域氣候較為干旱,一旦發(fā)生森林火災(zāi),撲救難度很大,容易造成大面積森林損失。統(tǒng)計(jì)表明塔河地區(qū)人為火的比重大約占到60%以上[8]。

    目前關(guān)于該地區(qū)人為火的研究,主要分為空間分布格局研究[9-11]和氣象因素影響研究[12-14]兩個(gè)方向。然而,國外很多研究表明,地形、林分特征等因素對(duì)人為火的發(fā)生同樣具有顯著的影響[15-19]。但由于研究地空間地理數(shù)據(jù)和林分調(diào)查數(shù)據(jù)的缺乏,國內(nèi)在探討該地區(qū)立地條件、地理因素對(duì)人為火影響的問題中還存在一些不確定性。本文在收集整理包含立地條件和地理因素(包括坡度、坡向等地形因素以及植被、腐質(zhì)層厚度、地被物蓋度等環(huán)境因素)在內(nèi)的塔河地區(qū)林相圖基礎(chǔ)上,應(yīng)用ArcGIS10.0中的空間分析工具和SPSS19.0的邏輯斯蒂回歸模型對(duì)影響人為火發(fā)生的主要立地條件和地理因素進(jìn)行分析,研究結(jié)論為進(jìn)一步解釋該區(qū)域人為火發(fā)生的主導(dǎo)因素提供依據(jù)。

    1 研究地區(qū)與研究方法

    1.1 研究區(qū)概況

    塔河地區(qū)位于我國黑龍江大興安嶺林區(qū),屬于典型的北方針葉林帶,具有植物多樣性低、喬木組成樹種非常簡(jiǎn)單等特點(diǎn)。植被以興安落葉松林和白樺林為主。屬寒溫帶季風(fēng)氣候,年均氣溫-2℃~4℃,年降水量350 mm~500 mm,相對(duì)濕度70%~75%。年積雪期達(dá)5個(gè)月,林內(nèi)積雪深可達(dá)30cm ~50cm。歷史上塔河地區(qū)是火災(zāi)多發(fā)區(qū),人為火源和雷擊火源共同存在,從次數(shù)來說,人為火約占火災(zāi)總數(shù)的2/3。從面積來說,雷擊火的過火面積高于人為火。該地區(qū)人為火在空間尺度上多呈現(xiàn)聚集分布,且多集中在鐵路和居民區(qū)附近,時(shí)間尺度上多集中在5月至9月,其中5月和6月是人為火高發(fā)期[8,20]。

    圖1 研究區(qū)域示意圖Fig.1 Sketch map of the study area

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

    林火發(fā)生數(shù)據(jù)來源于大興安嶺地區(qū)塔河縣森林防火辦公室。數(shù)據(jù)包括1974年—2009年塔河地區(qū)林火發(fā)生情況(起火地理坐標(biāo)、起火原因、發(fā)生時(shí)間、過火面積等)。本文提取出歷年人為火發(fā)生的經(jīng)緯度坐標(biāo)。研究中所需基礎(chǔ)地理地形圖和矢量化林相圖等數(shù)據(jù)來源于東北林業(yè)大學(xué)森林經(jīng)理教研室(http://forestry.nefu.edu.cn/dlxx/msg.php?c=13)。1∶100萬的基礎(chǔ)地理圖和矢量化林相圖數(shù)據(jù)庫分別建于2000年和2003年。其中林相圖數(shù)據(jù)庫精確到小班,并對(duì)林下植被類型和土壤條件進(jìn)行了很好的記錄,為本研究提供了有力的數(shù)據(jù)支持。本文參考國外相關(guān)研究[4,15,17,19],并結(jié)合自身數(shù)據(jù)內(nèi)容,通過ArcGIS軟件運(yùn)用空間分析技術(shù)從獲取的林火及林分?jǐn)?shù)據(jù)庫中提取火點(diǎn)和隨機(jī)點(diǎn)到“河流”,“街區(qū)”,“鐵路”和“機(jī)耕路”等人為基礎(chǔ)設(shè)施距離為第一類非氣象因子。選取“林型”,“坡度”,“坡向”,“地被物蓋度”和“優(yōu)勢(shì)樹種”等林分和地形特征為第二類非氣象因子。

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

    1.3.1 數(shù)據(jù)集屬性表轉(zhuǎn)化與合并

    本文應(yīng)用ArcGIS10.0對(duì)人為火數(shù)據(jù)與基礎(chǔ)地理和林相圖數(shù)據(jù)進(jìn)行空間分析,計(jì)算了火點(diǎn)到河流面、街區(qū)、鐵路、機(jī)耕路及等級(jí)公路的距離。并將林相圖中的字符型變量進(jìn)行屬性轉(zhuǎn)換以便進(jìn)行統(tǒng)計(jì)分析,包括坡向、林型、經(jīng)營措施和林下優(yōu)勢(shì)樹種(表1和表2),在給字符型因子附數(shù)值型屬性值時(shí),以整數(shù)形式并以1為計(jì)數(shù)間隔。這些被賦予的屬性數(shù)值雖是人為定值,但它們的差異真實(shí)反應(yīng)了相對(duì)應(yīng)的字符型屬性變量的差異。因此在實(shí)際邏輯斯蒂模型運(yùn)算中,變量的顯著性水平能有效地說明該因子是否對(duì)人為火有顯著影響。此外,坡度級(jí)、地被物蓋度、齡級(jí)和郁閉度也是本文考慮的因子,但由于其本身以數(shù)值型記錄因此不需要進(jìn)行屬性轉(zhuǎn)換。為了建模分析(邏輯斯蒂回歸模型)影響人為火發(fā)生的關(guān)鍵因子,需要隨機(jī)創(chuàng)建一定數(shù)量(≥人為火數(shù)量)的非著火點(diǎn)(即沒有林火發(fā)生的點(diǎn))與林火點(diǎn)一同進(jìn)行模型擬合。

    本文分別按照火點(diǎn)與非火點(diǎn)1∶1和1∶3的比例創(chuàng)建非火點(diǎn)。非火點(diǎn)(即隨機(jī)點(diǎn))的創(chuàng)建數(shù)量或比例,國際上并沒有統(tǒng)一的標(biāo)準(zhǔn)。有學(xué)者按照研究時(shí)間范圍內(nèi)平均火點(diǎn)數(shù)量選取非火點(diǎn)[21],也有學(xué)者以非火點(diǎn)之間的最小距離為標(biāo)準(zhǔn)在整個(gè)研究區(qū)域內(nèi)選取隨機(jī)點(diǎn)[22]。這些方法各有利弊,但非火點(diǎn)的“隨機(jī)性”是選擇的基本原則。此外,由于實(shí)際情況往往是研究區(qū)域的火點(diǎn)空間分布少于或大大少于非火點(diǎn),所以在隨機(jī)選取非火點(diǎn)時(shí),其數(shù)量通常不能低于火點(diǎn)數(shù)量,這也是反應(yīng)了林火發(fā)生的客觀規(guī)律。本文選取兩種不同比例的隨機(jī)點(diǎn),主要為了比較分析及探討隨機(jī)點(diǎn)數(shù)量的變化對(duì)研究結(jié)論的潛在影響。本文按照火點(diǎn)屬性表創(chuàng)建過程創(chuàng)建非火點(diǎn)屬性表,并將兩者合并。本文共選取15個(gè)因子進(jìn)行建模分

    表1 林相圖屬性轉(zhuǎn)化屬性表ITable 1 The convert of attribute value of stock map I

    表2 林相圖屬性轉(zhuǎn)化屬性表IITable 2 The convert of attribute value of stock map II

    析,采用點(diǎn)數(shù)據(jù)處理方式,通過數(shù)據(jù)預(yù)處理(剔除屬性殘缺點(diǎn))分別得到82和165個(gè)因變量數(shù)據(jù)點(diǎn)(表3<樣本數(shù)>)。

    1.3.2 二項(xiàng)邏輯斯蒂回歸模型(Logistic regression model)

    在Logistic回歸模型中,對(duì)兩類判別問題進(jìn)行取值。第一類為有林火發(fā)生,記為y=1,第二類為無林火發(fā)生,y=0。根據(jù)指標(biāo)可以建立LogitP關(guān)于自變量(相關(guān)林相屬性和火點(diǎn)到河流,鐵路,公路等設(shè)施距離)的Logistic回歸方程:

    Logit P=b0+b1x1+b2x2+…+bmxm(1)

    得到概率估計(jì)公式:

    根據(jù)估計(jì)概率進(jìn)行判別歸類,通常情況下,如果估計(jì)概率小于0.5,則認(rèn)為沒有林火發(fā)生,如果估計(jì)概率大于0.5,則認(rèn)為有林火發(fā)生;如果不進(jìn)行四舍五入對(duì)數(shù)值進(jìn)行取舍,一般不會(huì)出現(xiàn)恒等于0.5的情況。需要指出的是,關(guān)于這個(gè)林火的判別概率,國際上并沒有明確地規(guī)定,這主要根據(jù)研究者的實(shí)際工作和需求來劃定。更多情況是計(jì)算出區(qū)域的林火概率分布圖來指導(dǎo)實(shí)際林火管理工作。本文應(yīng)用SPSS19.0軟件進(jìn)行邏輯斯蒂回歸擬合。這種分析已廣泛應(yīng)用于林火發(fā)生的預(yù)測(cè)預(yù)報(bào)上[18,21,23-25]。

    2 結(jié)果與分析

    2.1 按不同比例創(chuàng)建隨機(jī)點(diǎn)

    本文擬應(yīng)用二元邏輯斯蒂回歸模型分析人為火發(fā)生與人為基礎(chǔ)設(shè)施和林分相關(guān)信息的關(guān)系。由于邏輯斯蒂模型的應(yīng)用條件,因變量的取值應(yīng)為1或0,即著火和不著火。我們的林火統(tǒng)計(jì)數(shù)據(jù)中只包含有林火發(fā)生的數(shù)據(jù),因此需要?jiǎng)?chuàng)建一些無林火發(fā)生數(shù)據(jù)來進(jìn)行模型擬合。創(chuàng)建過程按照隨機(jī)原則,且非火點(diǎn)的數(shù)量要呈遞增的趨勢(shì),以便更好地理解和探討各非氣象因子與人為火發(fā)生的關(guān)系,因此本文隨機(jī)點(diǎn)按照著火點(diǎn)1倍和3倍的數(shù)量創(chuàng)建了隨機(jī)點(diǎn)。

    2.2 屬性表整合及統(tǒng)計(jì)描述

    本文計(jì)算了火點(diǎn)和不同比例隨機(jī)點(diǎn)到河流面、街區(qū)、鐵路、機(jī)耕路、等外公路及等級(jí)公路的距離(單位:m)(圖2,圖3),并在ArcGIS中將火點(diǎn)和隨機(jī)點(diǎn)圖層與林相圖屬性圖層進(jìn)行疊加運(yùn)算,合并各圖層屬性。應(yīng)用SPSS19.0對(duì)整合后的屬性數(shù)據(jù)進(jìn)行統(tǒng)計(jì)描述(表3)。

    表3 模型變量的基本統(tǒng)計(jì)Table 3 The descriptive statistic of variables of model

    圖2 火點(diǎn)與1∶1隨機(jī)點(diǎn)到人為基礎(chǔ)設(shè)施及河流面的空間直線距離。圖中a,b,c,d,e分別代表“火點(diǎn)+隨機(jī)點(diǎn)”到鐵路、河流面、居住、機(jī)耕路、和等級(jí)公路的距離(單位:m),每個(gè)顏色條紋代表5000m緩沖帶Fig.2 The straight distance from fire points and 1∶1random no-fire points to infrastructures and rivers.a(chǎn),b,c,d,e represent the distance to railway,river,residential area,Jigeng road and Dengji Road,respectively(unit:m).Each color stripe represents 5000mbuffer

    圖3 火點(diǎn)與1∶3隨機(jī)點(diǎn)到人為基礎(chǔ)設(shè)施及河流面的空間直線距離。圖中a,b,c,d,e分別代表“火點(diǎn)+隨機(jī)點(diǎn)”到鐵路、河流面、居住、機(jī)耕路、和等級(jí)公路的距離(單位:m),每個(gè)顏色條紋代表5000m緩沖帶Fig.3 The straight distance from fire points and 1∶3random no-fire points to infrastructures and rivers.a(chǎn),b,c,d,e represent the distance to railway,river,residential area,Jigeng road and Dengji Road,respectively(unit:m).Each color stripe represents 5000mbuffer

    2.3 模型擬合結(jié)果

    2.3.1 火點(diǎn)與隨機(jī)點(diǎn)1∶1擬合結(jié)果

    本文對(duì)1∶1創(chuàng)建隨機(jī)點(diǎn)的數(shù)據(jù)進(jìn)行二元邏輯斯蒂回歸擬合。模型的全變量擬合結(jié)果顯示,模型整體擬合結(jié)果R2(Cox &Snell)為0.517,模型參數(shù)擬合結(jié)果顯示“到鐵路距離”與人為火發(fā)生之間有顯著相關(guān)性(p=0.009),而其它自變量因子在全模型中(含有所有自變量的邏輯斯蒂模型)與人為火發(fā)生之間并未表現(xiàn)出顯著相關(guān)(表4)。此外,本文應(yīng)用向后逐步回歸法,將不顯著變量逐一從全模型中剔除,最后得到最優(yōu)模型(剔除不顯著自變量的邏輯斯蒂模型)。模型整體擬合度指標(biāo)R2(Cox &Snell)為0.465,最優(yōu)模型中“到鐵路距離”,“坡向”和“優(yōu)勢(shì)樹種”這三個(gè)變量與人為火發(fā)生有顯著相關(guān)性(表5)。

    表4 邏輯斯蒂全模型的參數(shù)擬合(1∶1)Table 4 Parameter estimation of Logistic regression full model(1∶1)

    表5 邏輯斯蒂最優(yōu)模型的參數(shù)擬合(1∶1)Table 5 Parameter estimation of Logistic regression best model(1∶1)

    2.3.2 火點(diǎn)與隨機(jī)點(diǎn)1∶3擬合結(jié)果

    在驗(yàn)證隨機(jī)點(diǎn)與火點(diǎn)1∶1的情況后,為了深入了解隨機(jī)點(diǎn)數(shù)量對(duì)決策因子的影響,本文運(yùn)用相同的模型擬合方法對(duì)1∶3(火點(diǎn):隨機(jī)點(diǎn))的情況進(jìn)行了建模分析。模型的全變量擬合結(jié)果顯示,模型整體擬合結(jié)果R2(Cox &Snell)為0.408,模型參數(shù)擬合結(jié)果顯示“到鐵路距離”與人為火發(fā)生之間有顯著相關(guān)性(p=0.003),而其它自變量因子在全模型中(含有所有自變量的邏輯斯蒂模型)與人為火發(fā)生之間并未表現(xiàn)出顯著相關(guān)(表6)。此外,同樣運(yùn)用向后逐步回歸法,將不顯著變量逐一從全模型中剔除,最后得到最優(yōu)模型,模型整體擬合度指標(biāo)R2(Cox &Snell)為0.347,最優(yōu)模型中“到鐵路距離”和“林型”這兩個(gè)變量與人為火發(fā)生有顯著相關(guān)性(表7)。

    表6 邏輯斯蒂全模型的參數(shù)擬合(1∶3)Table 6 Parameter estimation of Logistic regression full model(1∶3)

    表7 邏輯斯蒂全模型的參數(shù)擬合(1∶3)Table 7 Parameter estimation of Logistic regression full model(1∶3)

    本文雖然隨機(jī)點(diǎn)的選取比例不同,模型擬合結(jié)果也有所差異,但兩個(gè)模型的擬合結(jié)果具有很強(qiáng)的規(guī)律性,模型顯示了鐵路對(duì)人為火的顯著影響,其次,林型和優(yōu)勢(shì)樹種代表了森林可燃物類型,研究結(jié)果表明人為火的發(fā)生受可燃物類型的影響顯著。此外,在1∶1選取隨機(jī)點(diǎn)條件下,坡向也與人為火發(fā)生表現(xiàn)了顯著相關(guān)性,然而在隨機(jī)點(diǎn)比例增加的情況下,坡向的顯著性水平有所下降,因此,在運(yùn)用“隨機(jī)點(diǎn)”方法進(jìn)行相關(guān)性分析時(shí),要充分注意和分析隨機(jī)點(diǎn)數(shù)量對(duì)模型擬合結(jié)果的影響。圖4為根據(jù)林火和隨機(jī)點(diǎn)創(chuàng)建的研究區(qū)域人為火發(fā)生概率分布圖,圖中直觀地顯示了火險(xiǎn)高發(fā)區(qū)。

    圖4 研究區(qū)域人為火發(fā)生概率克里格插值結(jié)果Fig.4 The Kriging interpolation of possibilityof human-caused fire in the study area

    3 結(jié)論與討論

    本文應(yīng)用塔河地區(qū)歷年人為火信息數(shù)據(jù),結(jié)合研究地的基礎(chǔ)地理信息及矢量化林相圖,應(yīng)用Arc-GIS10.0和SPSS19.0等工具,通過圖層疊加計(jì)算,屬性轉(zhuǎn)換,邏輯斯蒂回歸模型運(yùn)算等一系列手段,分析了對(duì)我國北方針葉林帶人為火發(fā)生有顯著影響的非氣象因子。研究結(jié)果顯示“到鐵路距離”“優(yōu)勢(shì)樹種”和“林型”對(duì)人為火的發(fā)生有顯著影響,這同Cardille等[4]在北美針葉林帶所進(jìn)行的研究結(jié)論相似,他們的研究指出“鐵路密度”以及“林型”對(duì)研究區(qū)域人為火有顯著影響。此外,Maingi和Henry[16]在美國肯塔西州的研究也表明“到公路距離”和“公路密度”與人為火發(fā)生顯著相關(guān),這也說明了人為基礎(chǔ)設(shè)施的空間分布對(duì)人為火的重要影響。塔河林業(yè)局的鐵路建設(shè),在過去70,80年代主要為木材和省內(nèi)旅客運(yùn)輸服務(wù),由于設(shè)備和人們防火認(rèn)識(shí)不夠完善,時(shí)常發(fā)生機(jī)車的閘瓦火,和由于旅客隨窗扔煙頭及在鐵路沿線活動(dòng)引發(fā)的森林火災(zāi),這在當(dāng)?shù)氐牧只鸢l(fā)生情況統(tǒng)計(jì)資料中可以發(fā)現(xiàn)。本研究也證實(shí)了塔河林業(yè)局人為火發(fā)生的特殊性。

    本研究主要選取河流及一些人為基礎(chǔ)設(shè)施為研究對(duì)象,與國外同類研究相比[4,16-18],本文沒有考慮當(dāng)?shù)厝丝诿芏?,GDP等經(jīng)濟(jì)發(fā)展因素,主要由于我國公開的GDP或人口密度等統(tǒng)計(jì)數(shù)據(jù),最小尺度單位為縣而本文的研究區(qū)域位于塔河縣內(nèi),因此無法體現(xiàn)出空間異質(zhì)性,進(jìn)而無法同其它因子一樣進(jìn)行分析。另外,塔河林區(qū)的經(jīng)濟(jì)發(fā)展較其它以工業(yè)為主的地區(qū)落后,發(fā)展速度較慢,加之我國人口政策的限制,多年來塔河地區(qū)的人口數(shù)量和居住區(qū)范圍并沒有較為明顯的改變。因此,本文用居住區(qū)的位置,而不是大小作為變量,來探討其對(duì)人為火發(fā)生的影響。另一方面,本研究時(shí)間尺度為1974—2009年,而選取的鐵路、公路、河流、居民點(diǎn)、機(jī)耕路等是基于2000年的基礎(chǔ)地理圖,由于客觀數(shù)據(jù)條件的限制,本研究假設(shè)各地理、人為基礎(chǔ)設(shè)施在研究時(shí)間尺度內(nèi)保持不變[25]。實(shí)際情況是該地區(qū)河流、居民點(diǎn)及鐵路位置的變化相對(duì)較小,而林區(qū)的公路和機(jī)耕路隨著時(shí)間的推移和當(dāng)?shù)亟?jīng)濟(jì)發(fā)展呈現(xiàn)上升趨勢(shì),這可能對(duì)研究結(jié)論產(chǎn)生一定的影響,因此,在未來的研究中,應(yīng)盡可能使用不同時(shí)間點(diǎn)的基礎(chǔ)地理圖,來減少研究過程中的不穩(wěn)定性。

    目前關(guān)于林火與影響因子的空間分析研究,主要有網(wǎng)格取點(diǎn)和隨機(jī)取點(diǎn)兩種方式。網(wǎng)格取點(diǎn)是指對(duì)研究區(qū)域進(jìn)行網(wǎng)格化,規(guī)定網(wǎng)格大小。每一個(gè)網(wǎng)格代表一個(gè)火點(diǎn)或?qū)?yīng)的非火點(diǎn)(隨機(jī)點(diǎn)),如果網(wǎng)格面積大,如5Km×5Km,那么坡向,林型等這類因子將不具備唯一性,沒有辦法考慮進(jìn)模型擬合中,若柵格面積小(1Km×1Km或更?。﹦t會(huì)出現(xiàn)“非火點(diǎn)”與“火點(diǎn)”比例嚴(yán)重失衡現(xiàn)象,會(huì)導(dǎo)致模型擬合的偏差,從而無法判斷影響因子。因此在火點(diǎn)和非火點(diǎn)比例失衡的情況下,隨機(jī)取點(diǎn)是一個(gè)比較好的選擇。因此本文采用了1∶1和1∶3兩種隨機(jī)選點(diǎn)方式。在模型運(yùn)算過程中,我們發(fā)現(xiàn)隨著隨機(jī)點(diǎn)比例的提高影響因子“坡向”的顯著性水平有所下降,從1∶1時(shí)的顯著相關(guān),到1∶3時(shí)的不顯著。這種情況說明,不應(yīng)該簡(jiǎn)單地根據(jù)某一種隨機(jī)點(diǎn)選取比例而直接下結(jié)論。而是應(yīng)該對(duì)比不同比例隨機(jī)點(diǎn)下同一自變量因子在相同模型下的顯著性水平變化趨勢(shì)。這樣可以較為客觀地反映出對(duì)人為火有顯著影響的決策因子。此外,本文根據(jù)火點(diǎn)和隨機(jī)點(diǎn)的空間分布規(guī)律,對(duì)塔河地區(qū)進(jìn)行了人為火空間分布插值,為直觀了解研究地人為火火險(xiǎn)分布情況提供依據(jù)。本文研究結(jié)論對(duì)深入揭示研究地人為火發(fā)生規(guī)律具有一定意義。

    參考文獻(xiàn)

    [1]Podur J,et al.Spatial patterns of lightning-caused forest fires in Ontario.1976-1998[J].Econ Model,2003,164 (1):1-20.

    [2]Chas-Amill ML,et al.Spatial distribution of humancaused forest fires in Galicia(NWSpain)[A].Transactions on ecology and the environment[C],2010,137.WIT Press,Southampton.ISSN 1743~3541(online)2010:247-258.

    [3]Zhong MH,et al.Statistical analysis on current status of China forest fire safety[J].Fire Safety Journal,2003,38(3):257-269.

    [4]Cardille JA,et al.Environmental and social factors influencing wildfires in the Upper Midwest.United States [J].Ecological Applications,2001,11(1):111-127.[5]Prestemon JP,Butry DT.Time to burn:modeling wildland arson as an autoregressive crime function[J].A-merican J.of Agricultural Economics,2005,87(3):756-770.

    [6]Susan K.Todd,Jewkes HA.Wildland fire in Alaska:a history of organized fire suppression and management in the last frontier[J].Agricultural and Forestry Experiment Station Bulletin No.114University of Alaska Fairbanks March,2006.

    [7]趙輝.福建三明2000~2009年森林火災(zāi)時(shí)空分布及影響因素研究[D].北京:北京林業(yè)大學(xué),2012.趙輝.福建三明2000~2009年森林火災(zāi)時(shí)空分布及影響因素研究[J].北京林業(yè)大學(xué)學(xué)報(bào),2012(1).

    [8]郭福濤.基于空間分析和模型理論的大興安嶺地區(qū)林火分布與預(yù)測(cè)模型研究[D].哈爾濱:東北林大學(xué),2010.[9]金森,胡海清.黑龍江省林火規(guī)律研究I林火時(shí)空動(dòng)態(tài)與分布[J].林業(yè)科學(xué),2002,38(1):88-94.

    [10]張正祥,等.呼倫貝爾草原人為火空間分布格局[J].生態(tài)學(xué)報(bào),2013,33(7):2023-2031.

    [11]王明玉,等.林火在空間上的波動(dòng)性及其對(duì)全球變化的響應(yīng)[J].火災(zāi)科學(xué),2003,12(3):165-176.

    [12]趙鳳君,等.氣候變化對(duì)林火動(dòng)態(tài)的影響研究進(jìn)展[J].氣候變化研究進(jìn)展,2009,5(1):50-55.

    [13]郭福濤,等.不同模型對(duì)擬合大興安嶺林火發(fā)生與氣象因素關(guān)系的適用性[J].應(yīng)用生態(tài)學(xué)報(bào),2010,21 (1):159-164.

    [14]魏書精,等.氣候變化對(duì)我國林火發(fā)生規(guī)律的影響[J].森林防火,2011,1:30-34.

    [15]Syphard AD,et al.Human influence on California fire regimes[J].Ecological Application,2007,17(5):1388-1402.

    [16]Maingi JK,Henry MC.Factors influencing wildfire occurrence and distribution in eastern Kentucky,USA [J].International Journal of Wildland Fire,2007,16 (1):23-33.

    [17]Romero-Calcerrada R,et al.GIS analysis of spatial patterns of human-caused wildfire ignition risk in the SW of Madrid(central Spain)[J].Landscape Ecology,2008,23(3):341-354.

    [18]Martínez J,et al.Human-caused wildfire risk rating for prevention planning in Spain[J].Journal of Environmental Management,2009,90(2):1241-1252.

    [19]Padilla M,Vega-Garcia C.On the comparative importance of fire danger rating indices and their interaction with spatial and temporal fire occurrences in Spain[J].International Journal ofWildland fire,2011,20:46-58.

    [20]于成龍,等.大興安嶺塔河林業(yè)局林火動(dòng)態(tài)氣象條件分析[J].東北林業(yè)大學(xué)學(xué)報(bào),2007,35(8):23-27.

    [21]Martell DL,et al.Modelling seasonal variation in daily people-caused forest fire occurrence[J].Canadian Journal of Forest Research,1989,19(12):1555-1563.

    [22]Vasconcelos MJP de,et al.Spatial Prediction of Fire Ignition Probabilities:Comparing Logistic Regression and Neural Networks[J].Photogrammetric Engineering &Remote Sensing,2001,67(1):73-81.

    [23]Garcia CV,et al.ALogit model for predicting the daily Occurrence of Human caused forest fires[J].International journal of wildland fire,1995,5(2):101-111.

    Keyword:Boreal forest;Tahe area;Human-caused fire;Geographical factors;Logistic regression

    [24]Andrews PL,et al.Evaluation of fire danger rating indexes using logistic regression and percentile analysis [J].International Journal of Wildland Fire,2003,12 (2):213-226.

    [25]Chang Y,et al.Predicting fire occurrence patterns with logistic regression in Heilongjiang Province,China[J].Landscape Ecology,2013,28(10):1989-2004.

    Influence of geographical and environmental factors on human-caused fires in Tahe area

    SU Zhangwen,SONG Yuhui,GUO Futao,YANG Tingting,WANG Wenhui,LIU Aiqin
    (Fujian Agriculture and Forestry University,F(xiàn)uzhou 350002,China)

    Abstract:In this paper,we used the spatial analysis tools in ArcGIS10.0and Logistic regression model imbedded in SPSS 19.0 to identify the impact of human-activity and geographical factors on forest fires,based on the coordinates of human-caused fires during 1974-2009and the basic geographical map.The results showed that“the distance to railway”has a significant impact on the human-caused fires.Moreover,“the dominant tree”and“the forest type”are significantly correlated with human-caused fires under 1∶1and 1∶3conditions,respectively.It was also found that the rate of random sample selection affects the parameter estimation of the model.

    通訊作者:郭福濤,E-mail:guofutao@126.com;劉愛琴,E-mail:lxylaq@126.com

    作者簡(jiǎn)介:蘇漳文(1987-),男,福建漳州。福建農(nóng)林大學(xué)碩士研究生,林業(yè)碩士專業(yè),研究方向?yàn)樯峙嘤夹g(shù)。

    基金項(xiàng)目:福建農(nóng)林大學(xué)重點(diǎn)大學(xué)建設(shè)基金(6112C035);國家自然科學(xué)基金(31400552)。

    收稿日期:2014-12-18;修改日期:2015-01-13

    DOI:10.3969/j.issn.1004-5309.2015.01.03

    文章編號(hào):1004-5309(2015)(-)0016-10

    文獻(xiàn)標(biāo)識(shí)碼:A

    中圖分類號(hào):S762.2;X954

    猜你喜歡
    火點(diǎn)塔河林火
    無錫惠山區(qū)堅(jiān)持“六抓六強(qiáng)” 構(gòu)建林火防治銅墻鐵壁
    林火監(jiān)測(cè)系統(tǒng)在森林防火中的應(yīng)用與發(fā)展
    亞像元火點(diǎn)對(duì)紅外預(yù)警衛(wèi)星的輻射干擾特性
    半邊天
    人身上有5個(gè)祛火點(diǎn)
    點(diǎn)煙頌
    鹿鳴(2018年1期)2018-01-30 12:05:42
    江蘇省海門市如何實(shí)現(xiàn)連續(xù)4年秸稈焚燒“零火點(diǎn)”?
    以石為絹塔河寫意
    中華奇石(2015年12期)2015-07-09 18:34:54
    以石為絹 塔河寫意
    中華奇石(2015年11期)2015-07-09 18:34:16
    井壁修整器在塔河油田TK1286井的應(yīng)用
    斷塊油氣田(2012年6期)2012-03-25 09:54:03
    成年女人永久免费观看视频| 国产精品不卡视频一区二区| 国产在线一区二区三区精 | 国产日韩欧美在线精品| 中文资源天堂在线| 你懂的网址亚洲精品在线观看 | 国产中年淑女户外野战色| 非洲黑人性xxxx精品又粗又长| 亚洲国产精品久久男人天堂| av在线亚洲专区| 网址你懂的国产日韩在线| 亚洲成人精品中文字幕电影| 高清av免费在线| 国产精品精品国产色婷婷| videossex国产| 99热全是精品| 中文字幕人妻熟人妻熟丝袜美| 在线播放国产精品三级| 尾随美女入室| 日本免费一区二区三区高清不卡| 我的女老师完整版在线观看| 丝袜美腿在线中文| 久久久久久久久久久丰满| 国产又色又爽无遮挡免| 精品一区二区三区人妻视频| 国产一区二区在线观看日韩| 欧美性猛交黑人性爽| 色播亚洲综合网| 亚洲成人久久爱视频| 亚洲人成网站在线播| 久久人妻av系列| 欧美高清成人免费视频www| 国产精品一区二区在线观看99 | 国产熟女欧美一区二区| 成人三级黄色视频| 成年女人永久免费观看视频| 狠狠狠狠99中文字幕| 97热精品久久久久久| 国内精品宾馆在线| 九九爱精品视频在线观看| 欧美潮喷喷水| 中文天堂在线官网| 免费不卡的大黄色大毛片视频在线观看 | 一级毛片久久久久久久久女| 国产免费视频播放在线视频 | 99久国产av精品国产电影| 又爽又黄无遮挡网站| 黄色欧美视频在线观看| 夜夜看夜夜爽夜夜摸| 亚洲成人精品中文字幕电影| 菩萨蛮人人尽说江南好唐韦庄 | 欧美另类亚洲清纯唯美| 亚洲性久久影院| 国产毛片a区久久久久| 中文天堂在线官网| 亚洲内射少妇av| 午夜爱爱视频在线播放| 天堂网av新在线| 国产精品国产高清国产av| 日本免费一区二区三区高清不卡| 国产男人的电影天堂91| 中国国产av一级| 亚洲熟妇中文字幕五十中出| 精品国产三级普通话版| 美女大奶头视频| 神马国产精品三级电影在线观看| 永久免费av网站大全| 夜夜看夜夜爽夜夜摸| 国产精品一区二区性色av| 精品人妻熟女av久视频| 日韩欧美三级三区| 婷婷色麻豆天堂久久 | 一级毛片我不卡| 能在线免费观看的黄片| 女人被狂操c到高潮| 精品99又大又爽又粗少妇毛片| 成人二区视频| 最近最新中文字幕大全电影3| 久久鲁丝午夜福利片| 精品人妻偷拍中文字幕| 一边亲一边摸免费视频| 水蜜桃什么品种好| 日本一本二区三区精品| 欧美又色又爽又黄视频| 精品久久国产蜜桃| 乱系列少妇在线播放| 国产精品永久免费网站| 在线a可以看的网站| 美女大奶头视频| 青春草亚洲视频在线观看| 亚洲欧美日韩高清专用| 女人久久www免费人成看片 | 欧美丝袜亚洲另类| 欧美97在线视频| h日本视频在线播放| 婷婷色麻豆天堂久久 | 国产亚洲精品av在线| 麻豆久久精品国产亚洲av| 中文欧美无线码| 天堂√8在线中文| 国产精品久久视频播放| 久久精品国产鲁丝片午夜精品| 少妇熟女aⅴ在线视频| 又粗又爽又猛毛片免费看| 亚洲五月天丁香| 狠狠狠狠99中文字幕| 色视频www国产| 午夜视频国产福利| 在线观看av片永久免费下载| 国产一区二区亚洲精品在线观看| 一区二区三区乱码不卡18| 久久草成人影院| 久久久久久久久久成人| 久久国产乱子免费精品| 深爱激情五月婷婷| 成人鲁丝片一二三区免费| 亚洲精华国产精华液的使用体验| 老司机影院毛片| 久久久成人免费电影| 日日干狠狠操夜夜爽| 久久人人爽人人片av| 免费一级毛片在线播放高清视频| ponron亚洲| 黑人高潮一二区| 一区二区三区免费毛片| 丝袜喷水一区| 精品人妻偷拍中文字幕| 一级毛片久久久久久久久女| 特大巨黑吊av在线直播| 国内精品美女久久久久久| 中文字幕免费在线视频6| 亚洲欧美中文字幕日韩二区| 亚洲国产成人一精品久久久| 一级毛片aaaaaa免费看小| 国产伦在线观看视频一区| 欧美区成人在线视频| av在线观看视频网站免费| 久久久久久久久久黄片| 久久久精品欧美日韩精品| 午夜福利在线观看吧| 在线免费十八禁| 免费看日本二区| 中文精品一卡2卡3卡4更新| 欧美精品国产亚洲| 国产视频首页在线观看| 欧美xxxx黑人xx丫x性爽| 成人性生交大片免费视频hd| 又爽又黄无遮挡网站| 成人美女网站在线观看视频| 最近最新中文字幕大全电影3| 在线天堂最新版资源| 又黄又爽又刺激的免费视频.| 最近最新中文字幕大全电影3| 联通29元200g的流量卡| 亚洲人成网站高清观看| 别揉我奶头 嗯啊视频| 干丝袜人妻中文字幕| 午夜激情欧美在线| 春色校园在线视频观看| 美女xxoo啪啪120秒动态图| 免费不卡的大黄色大毛片视频在线观看 | 久久这里只有精品中国| eeuss影院久久| 91av网一区二区| 精华霜和精华液先用哪个| 亚洲美女视频黄频| 国国产精品蜜臀av免费| 99国产精品一区二区蜜桃av| 水蜜桃什么品种好| 中文乱码字字幕精品一区二区三区 | 亚洲av成人精品一二三区| 日本与韩国留学比较| 国产午夜精品久久久久久一区二区三区| 男人舔女人下体高潮全视频| 国产高清视频在线观看网站| 免费看美女性在线毛片视频| 久久午夜福利片| 能在线免费观看的黄片| 晚上一个人看的免费电影| 亚洲天堂国产精品一区在线| 亚洲欧美成人精品一区二区| 亚洲经典国产精华液单| 中文天堂在线官网| 男人舔女人下体高潮全视频| 淫秽高清视频在线观看| 一区二区三区免费毛片| 日韩 亚洲 欧美在线| 午夜精品一区二区三区免费看| 亚洲中文字幕一区二区三区有码在线看| 国产精品人妻久久久影院| 久久久久网色| 日韩欧美在线乱码| 看非洲黑人一级黄片| 美女xxoo啪啪120秒动态图| 色综合色国产| 美女大奶头视频| av视频在线观看入口| 国产久久久一区二区三区| 小说图片视频综合网站| 少妇的逼好多水| 日韩亚洲欧美综合| 日韩在线高清观看一区二区三区| 亚洲国产精品sss在线观看| 亚洲国产最新在线播放| 成人美女网站在线观看视频| 国产三级在线视频| 国模一区二区三区四区视频| 免费人成在线观看视频色| 美女大奶头视频| 少妇猛男粗大的猛烈进出视频 | 免费看美女性在线毛片视频| 超碰97精品在线观看| 免费无遮挡裸体视频| 亚洲av免费在线观看| 国产老妇伦熟女老妇高清| 国产精品电影一区二区三区| 色综合亚洲欧美另类图片| 国产真实乱freesex| 波多野结衣巨乳人妻| 国产高潮美女av| 嫩草影院新地址| av线在线观看网站| 免费观看在线日韩| 狂野欧美激情性xxxx在线观看| 99久久人妻综合| 你懂的网址亚洲精品在线观看 | 国产色爽女视频免费观看| 亚洲av中文av极速乱| 观看美女的网站| 蜜桃久久精品国产亚洲av| 白带黄色成豆腐渣| 最近视频中文字幕2019在线8| 蜜桃亚洲精品一区二区三区| 中文亚洲av片在线观看爽| 一级二级三级毛片免费看| 亚洲美女搞黄在线观看| 国产女主播在线喷水免费视频网站 | 日韩av在线免费看完整版不卡| 偷拍熟女少妇极品色| 欧美高清性xxxxhd video| 精品熟女少妇av免费看| 欧美一区二区亚洲| 美女被艹到高潮喷水动态| 国产精品无大码| 一区二区三区免费毛片| 欧美日本视频| 国产精品av视频在线免费观看| 狂野欧美白嫩少妇大欣赏| 久久精品国产亚洲av涩爱| 久久这里只有精品中国| 国产精品久久久久久精品电影小说 | h日本视频在线播放| 亚洲精品乱久久久久久| 久久精品久久久久久久性| 亚洲精华国产精华液的使用体验| 七月丁香在线播放| 乱系列少妇在线播放| 一个人看视频在线观看www免费| 免费av毛片视频| 91久久精品国产一区二区成人| 久久人妻av系列| 欧美又色又爽又黄视频| 亚洲高清免费不卡视频| 国产免费又黄又爽又色| h日本视频在线播放| av免费在线看不卡| 三级毛片av免费| 欧美激情久久久久久爽电影| 成年女人永久免费观看视频| 99久久成人亚洲精品观看| 午夜福利高清视频| 欧美一区二区亚洲| 亚洲在线观看片| 国产三级中文精品| 99久久无色码亚洲精品果冻| 真实男女啪啪啪动态图| 欧美一区二区精品小视频在线| 嫩草影院新地址| 国产大屁股一区二区在线视频| 国产精品无大码| 国产美女午夜福利| 三级国产精品欧美在线观看| 一边摸一边抽搐一进一小说| av福利片在线观看| 亚洲精品影视一区二区三区av| 久久久久久久久久久丰满| 色网站视频免费| 国产精品一区二区三区四区免费观看| 不卡视频在线观看欧美| 国产极品天堂在线| 男女国产视频网站| 国产老妇伦熟女老妇高清| 国产亚洲av片在线观看秒播厂 | 波多野结衣高清无吗| 亚洲av福利一区| 欧美色视频一区免费| 亚洲成av人片在线播放无| 成人午夜精彩视频在线观看| 色5月婷婷丁香| 人体艺术视频欧美日本| 国产精品熟女久久久久浪| 国产精品国产三级国产专区5o | 国产高潮美女av| 久久久久久久午夜电影| 久久人人爽人人片av| 99久国产av精品| 亚洲精华国产精华液的使用体验| 精品一区二区三区人妻视频| 中文字幕亚洲精品专区| 男人狂女人下面高潮的视频| 一边摸一边抽搐一进一小说| 久久久久性生活片| 99久久精品一区二区三区| 两个人的视频大全免费| 五月伊人婷婷丁香| 午夜亚洲福利在线播放| 中文乱码字字幕精品一区二区三区 | 日韩在线高清观看一区二区三区| 搞女人的毛片| 神马国产精品三级电影在线观看| 国产黄a三级三级三级人| 丰满少妇做爰视频| 国产亚洲一区二区精品| 乱码一卡2卡4卡精品| 欧美激情久久久久久爽电影| 亚洲真实伦在线观看| 国产精品爽爽va在线观看网站| 久99久视频精品免费| 色尼玛亚洲综合影院| 白带黄色成豆腐渣| 国产av在哪里看| 欧美一区二区国产精品久久精品| 少妇熟女欧美另类| 亚洲人成网站在线观看播放| 国产乱人视频| 男女边吃奶边做爰视频| 久久国内精品自在自线图片| 国产一区有黄有色的免费视频 | 最近最新中文字幕免费大全7| 深爱激情五月婷婷| 在线免费十八禁| 中文字幕久久专区| 精品酒店卫生间| 国产伦精品一区二区三区视频9| 亚洲av熟女| 久久久久久久国产电影| 18+在线观看网站| 91久久精品电影网| 91久久精品国产一区二区三区| 老司机福利观看| 国产高潮美女av| 嫩草影院入口| 国产精品99久久久久久久久| 有码 亚洲区| 精品人妻熟女av久视频| 国产高清不卡午夜福利| 91久久精品国产一区二区三区| 国产亚洲av片在线观看秒播厂 | 国产免费又黄又爽又色| 久久99精品国语久久久| 欧美xxxx性猛交bbbb| 看免费成人av毛片| 听说在线观看完整版免费高清| 国产日韩欧美在线精品| 国产av在哪里看| 少妇人妻精品综合一区二区| 成年女人看的毛片在线观看| 亚洲国产精品成人综合色| 久久久精品大字幕| 自拍偷自拍亚洲精品老妇| 丰满少妇做爰视频| 亚洲国产色片| 美女黄网站色视频| 97在线视频观看| 久久国内精品自在自线图片| 99久国产av精品| 中文资源天堂在线| 看非洲黑人一级黄片| 国产精品久久视频播放| 久久精品国产亚洲av涩爱| 男人舔女人下体高潮全视频| 亚洲精品456在线播放app| 91aial.com中文字幕在线观看| 国产麻豆成人av免费视频| 美女脱内裤让男人舔精品视频| 亚洲av一区综合| 国产 一区 欧美 日韩| www日本黄色视频网| 一级黄片播放器| 国产探花极品一区二区| 淫秽高清视频在线观看| 99久久精品一区二区三区| 99久久中文字幕三级久久日本| 老司机福利观看| 精品人妻偷拍中文字幕| 少妇人妻精品综合一区二区| 精品一区二区三区人妻视频| 亚洲精品国产av成人精品| 岛国毛片在线播放| 舔av片在线| 91狼人影院| 久久久精品欧美日韩精品| 欧美97在线视频| 日日撸夜夜添| 成人亚洲精品av一区二区| 人人妻人人澡欧美一区二区| 亚洲人与动物交配视频| 国产成人精品久久久久久| videossex国产| 久久精品国产亚洲av涩爱| 天堂√8在线中文| 丝袜喷水一区| 99热网站在线观看| 国产精品久久电影中文字幕| 久久精品国产自在天天线| 久99久视频精品免费| 小说图片视频综合网站| 国产黄片视频在线免费观看| 一级爰片在线观看| 日本午夜av视频| 成年av动漫网址| 黄片无遮挡物在线观看| 在线免费观看不下载黄p国产| 麻豆久久精品国产亚洲av| 亚洲va在线va天堂va国产| 国产精品国产高清国产av| 超碰av人人做人人爽久久| 免费一级毛片在线播放高清视频| 亚洲av成人av| 久久人妻av系列| 99久久精品热视频| 亚洲久久久久久中文字幕| 婷婷色麻豆天堂久久 | 十八禁国产超污无遮挡网站| 色综合亚洲欧美另类图片| 啦啦啦啦在线视频资源| av播播在线观看一区| 麻豆成人午夜福利视频| 国产精品久久久久久久电影| 男人和女人高潮做爰伦理| 国产精品久久久久久久电影| 欧美xxxx黑人xx丫x性爽| 精品酒店卫生间| 免费观看性生交大片5| 伦精品一区二区三区| 中文字幕熟女人妻在线| 国产精品人妻久久久久久| 又粗又硬又长又爽又黄的视频| 高清在线视频一区二区三区 | 色综合站精品国产| 日韩高清综合在线| 欧美极品一区二区三区四区| 蜜臀久久99精品久久宅男| 国产亚洲av嫩草精品影院| 国产免费男女视频| 国产视频首页在线观看| 观看免费一级毛片| 日本wwww免费看| 久久精品国产99精品国产亚洲性色| 22中文网久久字幕| 亚洲av免费高清在线观看| 91午夜精品亚洲一区二区三区| 韩国高清视频一区二区三区| 69人妻影院| 色综合色国产| 1024手机看黄色片| 2021天堂中文幕一二区在线观| 中文字幕免费在线视频6| 国产精品国产三级国产av玫瑰| 日韩一区二区视频免费看| 亚洲av一区综合| 三级男女做爰猛烈吃奶摸视频| 国产片特级美女逼逼视频| 国产精品福利在线免费观看| 99久久成人亚洲精品观看| 国产成人福利小说| 国产精品乱码一区二三区的特点| 一级毛片我不卡| 国产乱人视频| 精品99又大又爽又粗少妇毛片| 黑人高潮一二区| 波多野结衣巨乳人妻| 国产黄片视频在线免费观看| av.在线天堂| 少妇熟女欧美另类| 国产伦精品一区二区三区视频9| 亚洲三级黄色毛片| 美女xxoo啪啪120秒动态图| 尾随美女入室| 听说在线观看完整版免费高清| 国产欧美日韩精品一区二区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 午夜免费男女啪啪视频观看| 国产成人91sexporn| 亚洲精品一区蜜桃| 久久婷婷人人爽人人干人人爱| 日韩一区二区视频免费看| 色综合亚洲欧美另类图片| 一卡2卡三卡四卡精品乱码亚洲| 成人三级黄色视频| 啦啦啦啦在线视频资源| 自拍偷自拍亚洲精品老妇| 韩国av在线不卡| 波多野结衣高清无吗| 三级毛片av免费| 日韩av不卡免费在线播放| av女优亚洲男人天堂| 三级国产精品欧美在线观看| 国产精品嫩草影院av在线观看| 免费观看在线日韩| 美女cb高潮喷水在线观看| 精品国产露脸久久av麻豆 | 男女国产视频网站| 午夜久久久久精精品| 国产精品熟女久久久久浪| 欧美激情国产日韩精品一区| 我要搜黄色片| 搞女人的毛片| av视频在线观看入口| www.av在线官网国产| 免费搜索国产男女视频| 免费av不卡在线播放| 18禁在线无遮挡免费观看视频| 国产老妇女一区| 99久国产av精品国产电影| 国产精品一区二区性色av| 简卡轻食公司| 亚洲激情五月婷婷啪啪| 少妇人妻精品综合一区二区| 欧美又色又爽又黄视频| 秋霞在线观看毛片| 久久精品人妻少妇| 热99re8久久精品国产| 国产爱豆传媒在线观看| 亚洲精品国产av成人精品| 日本-黄色视频高清免费观看| 一级毛片我不卡| 中文亚洲av片在线观看爽| 毛片女人毛片| 成年女人看的毛片在线观看| 精品一区二区三区视频在线| 日本免费一区二区三区高清不卡| 国产成人aa在线观看| 日韩欧美 国产精品| 久99久视频精品免费| 中文字幕av成人在线电影| 狂野欧美白嫩少妇大欣赏| 亚洲精品乱码久久久久久按摩| 在线免费观看不下载黄p国产| 日本五十路高清| 美女高潮的动态| av国产免费在线观看| 亚洲人成网站在线播| 国产私拍福利视频在线观看| 黄色欧美视频在线观看| 久久精品综合一区二区三区| 国产午夜精品论理片| 69av精品久久久久久| 色网站视频免费| 亚洲欧美日韩高清专用| 久久久亚洲精品成人影院| 韩国高清视频一区二区三区| www日本黄色视频网| 精品人妻一区二区三区麻豆| 国产精品一区二区三区四区久久| 可以在线观看毛片的网站| 综合色av麻豆| 麻豆成人av视频| 亚洲aⅴ乱码一区二区在线播放| 欧美最新免费一区二区三区| 午夜精品在线福利| 国内少妇人妻偷人精品xxx网站| 热99re8久久精品国产| 国产精品无大码| 亚洲国产精品久久男人天堂| 国产精品国产三级国产av玫瑰| 在现免费观看毛片| 亚洲无线观看免费| 精品免费久久久久久久清纯| 一个人看的www免费观看视频| 日韩欧美精品免费久久| 亚洲电影在线观看av| 热99re8久久精品国产| 26uuu在线亚洲综合色| 日韩视频在线欧美| 青青草视频在线视频观看| h日本视频在线播放| 99久久精品国产国产毛片| 亚洲精品亚洲一区二区| 亚洲国产成人一精品久久久| 国产精品人妻久久久影院| 国产伦理片在线播放av一区| 久久午夜福利片| 欧美精品一区二区大全| 国产精品久久久久久久久免| 国产私拍福利视频在线观看| 国产精品人妻久久久久久| 嫩草影院入口| 永久网站在线| 成人鲁丝片一二三区免费| 国产成人91sexporn| 午夜视频国产福利| 日韩av在线大香蕉| 我的老师免费观看完整版| 少妇人妻精品综合一区二区| 国产 一区精品| 精品久久久久久久久久久久久| 高清av免费在线| 日韩,欧美,国产一区二区三区 | 99久久精品一区二区三区| 国产高清国产精品国产三级 | 午夜福利高清视频| 男人狂女人下面高潮的视频| 国产精品久久久久久精品电影| 久久精品熟女亚洲av麻豆精品 |