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

    運用地理加權泊松模型估測天然次生林枯損量分布1)

    2018-01-15 06:02:35張凌宇劉兆剛董靈波
    東北林業(yè)大學學報 2018年1期
    關鍵詞:莫蘭泊松次生林

    張凌宇 劉兆剛 董靈波

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

    天然次生林具有樹種繁多、分布范圍廣等特點,在我國東北、西南和東南地區(qū)均有分布,天然次生林不僅作為用材林以及林副產(chǎn)品的重要生產(chǎn)基地,在涵養(yǎng)水源以及維持生態(tài)平衡方面等都起到了重要的作用。我國次生林面積約占全國森林面積的46.2%,對天然次生林進行全面規(guī)劃、改善林分質量,使其盡快成為重要的森林資源。林分枯損影響著森林資源的動態(tài)變化,在森林演替過程中也起著至關重要的作用,因此,對天然次生林枯損株數(shù)分布的研究具有重要的理論意義。

    在20世紀60年代,我國學者就已經(jīng)開始了對林分枯損的研究,李克志[1]利用帶嶺地區(qū)4種天然林的實測數(shù)據(jù),對紅松等4種天然林枯損率進行了計算,對各齡組的枯損原因進行了分析;曾偉生[2]利用2006年西藏地區(qū)云杉林的實測固定樣地數(shù)據(jù),在非線性混合估算模型的基礎上,構建了云杉林消耗率、枯損率、收獲率三者一致兼容的估算模型,這一方法為西藏自治區(qū)天然云杉林的枯損量預估提供了參考依據(jù)。隨著統(tǒng)計模型在林業(yè)上的應用,一些學者開始利用全局模型對林分枯損株數(shù)進行預測和分析,張雄清等[3]利用落葉松(Larixolgensis)林分連續(xù)觀測數(shù)據(jù),分別利用泊松回歸模型、負二項模型、零膨脹模型擬合林分枯損株數(shù)。由于在林業(yè)調(diào)查中所收集的數(shù)據(jù)往往處于不同的地理位置,無論是研究枯損概率的logistic模型,還是研究枯損株數(shù)的泊松模型,實際上都沒有考慮到空間異質性,因此,空間異質性在有關林業(yè)的各項研究中普遍存在,也往往是研究者最容易忽略的問題。

    20世紀90年代,英國學者提出了地理加權回歸模型(GWR),解決空間非平穩(wěn)性以及空間異質性問題[4-5]。在國內(nèi)林業(yè)領域,有關GWR模型的研究一般集中在森林碳儲量以及生物量預估等方面[6-8],將地理加權回歸模型應用到林分枯損株數(shù)空間分布中的研究還未見報道。地理加權泊松模型(GWPR)是以地理加權回歸模型(GWR)為基礎擴展建立的局域泊松模型,可以用來構建以林木枯損株數(shù)為因變量的局域模型。本研究采用帽兒山實驗林場2期固定樣地數(shù)據(jù),分別建立了模擬天然次生林枯損株數(shù)的全局泊松模型和局域模型,將為大范圍內(nèi)天然次生林枯損株數(shù)預測以及解決計數(shù)型數(shù)據(jù)的空間非平穩(wěn)性問題提供理論基礎。

    1 研究區(qū)概況

    帽兒山實驗林場位于黑龍江省尚志市西部,地理坐標127°30′~127°34′E、45°20′~45°25′N,總面積26 496 hm2,共劃分10個森林經(jīng)營施業(yè)區(qū)、151個林班。帽兒山實驗林場的植被類型屬于長白山植物區(qū)系,是地帶性頂級植被闊葉紅松林被人為干擾破壞后形成的較典型天然次生林。主要林分類型有楊樺林、珍貴硬闊混交林、色木(Acermono)、柞樹(Quercusmongolica)為主的硬雜木林,其它主要樹種還有水曲柳(Fraxinusmandshurica)、胡桃楸(Junglusmandshurica)、黃菠蘿(Phellodendronamurense)、楓樺(Betualcostata)、白樺(Betualplatyphlla)、山楊(Populusdavidiana)等[9]。

    2 研究方法

    2.1 數(shù)據(jù)來源

    本研究數(shù)據(jù)來源于帽兒山實驗林場森林資源二類調(diào)查的2期固定樣地數(shù)據(jù)(2004年和2016年),在所有的固定樣地中剔除宜林沼澤、非林地、農(nóng)田等樣地,同時檢查數(shù)據(jù)中的異常值,對于異常值予以剔除,最終選擇101塊天然次生林樣地作為研究數(shù)據(jù),樣地面積均為0.06 hm2,數(shù)據(jù)涵蓋研究區(qū)內(nèi)的所有10個施業(yè)區(qū)。

    2.2 模型變量的選擇

    采用混合逐步選擇法對模型變量進行選擇,混合逐步選擇法是在常規(guī)模型的基礎上,根據(jù)一定的顯著性標準,逐步將所選變量剔除或添加進模型的方法[10],本研究將顯著性標準設為α=0.05,最終得到海拔、坡度、林分平均胸徑、林分密度和蓄積量共5個影響次生林枯損株數(shù)分布的變量。以往國內(nèi)外研究結果表明,林分枯損株數(shù)主要與林分密度、林分平均胸徑以及單位面積蓄積量有關[11-12],同時,林分單位面積枯損株數(shù)與土壤因子、立地條件和林分密度有著密切的關系[13],而坡度以及海拔作為重要的地形因子影響著土壤水分以及養(yǎng)分的含量[14],間接影響著天然次生林枯損株數(shù)分布。在所有的模型中對獨立變量均采取了標準化處理,模型各變量的統(tǒng)計量見表1。

    表1 模型各變量的基本統(tǒng)計量

    2.3 全局泊松模型

    全局泊松模型(GP)作為分析計數(shù)型數(shù)據(jù)的一種方法,來模擬天然次生林枯損株數(shù)的分布情況,其概率密度函數(shù)如下:

    (1)

    式中:P(Y=y)表示在一段時間內(nèi)時間Y發(fā)生次數(shù)的概率;λ為隨機變量Y的期望和方差,且期望和方差相等,即E(Y)=λ,Var(Y)=λ;當y=0時,表示天然次生林枯損株數(shù)為0的概率;當y=1時,表示天然次生林枯損株數(shù)為1的概率;當y=k時,即P(y=k)=e-λλk/k!,表示枯損株數(shù)為k的概率。全局泊松回歸模型的形式如下:

    log(E(Y))=logμ=β0+β1x1+β2x2+β3x3+β4x4+

    β5x5。

    (2)

    式中:log(.)為一個鏈接函數(shù);μ表示枯損株數(shù);β1~β5為模型的回歸系數(shù);x1~x5分別為海拔、坡度、平均胸徑、林分密度以及單位面積蓄積量。

    2.4 地理加權泊松模型

    地理加權泊松模型屬于全局泊松回歸模型在局域形式上的一種表達,在計算過程中可以將每塊樣地的地理坐標都納入到模型中,具體形式如下:

    log(e(Y))=logμ=β0(ui,vi)+β1(ui,vi)x1+β2(ui,vi)x2+β3(ui,vi)x3+β4(ui,vi)x4+

    β5(ui,vi)x5。

    (3)

    式中:log(.)為一個鏈接函數(shù);β1(ui,vi)~β5(ui,vi)為樣地點(ui,vi)上的回歸系數(shù);x1~x5分別為海拔、坡度、平均胸徑、林分密度以及單位面積蓄積量。

    2.5 模型殘差的空間自相關性分析

    真實值與預測值之間的差為模型殘差,通過引入莫蘭指數(shù)(MoranI)這一概念來檢驗模型殘差中的空間相關性,近年來,莫蘭指數(shù)被普遍運用到在林業(yè)研究中[6,15]。全局莫蘭指數(shù)計算公式如下:

    (4)

    (5)

    在研究中,首先設定一個零假設,即當莫蘭指數(shù)的值為零時,假設模型殘差服從隨機分布,對于本研究而言,即天然次生林枯損株數(shù)的空間隨機分布。通過計算莫蘭指數(shù),給出一個可以判斷出是否拒絕零假設的閾值,稱之為Z值,如果Z值在-1.96~1.96的區(qū)間內(nèi),則P>0.05(α=0.05),表示不能拒絕模型殘差服從隨機分布這一零假設,認為模型殘差表現(xiàn)出的模式在很大程度上是隨機空間過程產(chǎn)生的結果;反之,如果不在這一區(qū)間內(nèi),則認為模型殘差的空間模式不是隨機過程產(chǎn)生的,在這種情況下可以拒絕零假設,所表現(xiàn)出來的空間分布模式為具有統(tǒng)計性顯著的聚類或分散模式,這種方法可以判斷n個樣點之間的觀測值是否存在空間自相關性[6]。

    2.6 模型評價

    本研究使用均方誤差(MSE)以及赤池信息準則(AIC)值來比較2種模型的擬合效果,當MSE值越小時,說明模型的精度越高,模型的擬合效果越好。具體公式如下:

    (6)

    赤池信息準則是在熵的概念基礎上建立的一種衡量統(tǒng)計模型擬合效果優(yōu)良的指標[16],在模型比較時,AIC值越小的模型,說明該模型擬合效果越好。具體公式如下:

    AIC=2p-2ln(L)。

    (7)

    式中:L為最大似然函數(shù)值;p為變量個數(shù)。

    本文采用SAS9.4通過混合逐步選擇法對影響天然次生林枯損株數(shù)的變量進行篩選,并建立全局泊松模型;采用GWR4.0建立局域GWPR模型,采用ROOKCASE軟件來計算全局莫蘭指數(shù)以及局域莫蘭指數(shù),利用反距離權重法(IDW)對局域模型各變量對應參數(shù)的空間分布圖、局域莫蘭指數(shù)的空間分布圖、局域Z值的空間分布圖以及全局和局域模型枯損株數(shù)的空間分布圖分別進行繪制。

    3 結果與分析

    3.1 模型擬合

    3.1.1 全局泊松模型和局域模型的擬合結果

    表2分別給出了全局泊松模型參數(shù)估計值、標準誤差、顯著性檢驗以及局域模型系數(shù)的擬合結果。對于全局泊松模型來說,所有的模型系數(shù)均通過了α=0.05水平下的顯著性檢驗,海拔、坡度、林分密度以及單位面積蓄積量的參數(shù)估計值均為正值,說明這4個變量與天然次生林枯損株數(shù)呈正相關關系,相反地,只有林分平均胸徑的參數(shù)估計值為負值,說明林分平均胸徑與天然次生林枯損株數(shù)呈負相關關系。在所有影響天然次生林枯損株數(shù)的因子中,林分因子對枯損株數(shù)的影響最大,其中林分平均胸徑對次生林枯損株數(shù)的影響最大,其次是林分密度,最后是單位面積蓄積量,海拔和坡度等地形因子對于次生林枯損株數(shù)的影響要小于林分因子。

    表2 全局泊松模型和局域模型的參數(shù)估計值

    全局泊松模型只有一組固定的模型參數(shù),局域模型的參數(shù)則是連續(xù)變化的,表2給出了5組模型參數(shù)統(tǒng)計量,分別為最小值、Q1(四分位數(shù)下限值)、中值、Q3(四分位數(shù)上限值)以及最大值。

    由表2可知,局域模型參數(shù)隨著地理位置的不同產(chǎn)生了一定的變化范圍,在全局泊松模型中β0=1.668 3,β0的參數(shù)范圍為-0.127 2~3.043 2,同理,在全局泊松模型中β1=0.452 7,β2=0.404 7,β3=-1.427 2,β4=1.128 3,β5=0.966 4,對應局域模型中β1的范圍為-1.435 9~1.013 4,β2的范圍為-0.238 9~2.044 1,β3的范圍為-3.362 9~2.004 8,β4的范圍為-1.575 7~5.192 5,β5的范圍為-2.967 6~2.829 3。

    3.1.2 局域模型各參數(shù)的空間分布

    從圖1A可以看出,在研究區(qū)的東南部地區(qū)以及東北部小范圍內(nèi),海拔參數(shù)的估計值為負值,說明在這些區(qū)域內(nèi),隨著海拔的逐漸增高,天然次生林枯損株數(shù)逐漸減少,而在研究區(qū)的中部地區(qū)以及西部地區(qū),海拔參數(shù)的估計值為正值,說明在這些區(qū)域內(nèi)枯損株數(shù)隨著海拔高度的增加而增加。

    圖1 GWPR模型各變量對應參數(shù)的空間分布

    圖1B中,在研究區(qū)的北部以及南部地區(qū)很小范圍內(nèi)的坡度參數(shù)值為負值,所占區(qū)域的比例很小,在大部分范圍內(nèi),隨著坡度的逐漸增加,林分枯損株數(shù)逐漸增加。

    圖1C中,西部地區(qū)及東南部分地區(qū)林分平均胸徑的參數(shù)估計值為正,說明在這兩個區(qū)域內(nèi),林分平均胸徑與枯損株數(shù)呈正相關關系,在其他區(qū)域內(nèi)林分平均胸徑與枯損株數(shù)呈負相關關系,從總體上看,在整個研究區(qū)范圍內(nèi),自西南向東北方向林分平均胸徑參數(shù)估計值的變化范圍從正到負。

    在圖1D中,林分密度的參數(shù)估計值除了在東北部地區(qū)極小范圍內(nèi)為負值外,在其他大部分地區(qū)林分密度的參數(shù)估計值均為正值,說明在這些區(qū)域內(nèi),隨著林分密度的增加,天然次生林枯損株數(shù)也隨之增加。

    在圖1E中,林分單位面積蓄積量的參數(shù)估計值自東南至西北方向均為負值,參數(shù)估計值為負值的區(qū)域約占整個研究區(qū)域面積的50%,說明在此研究區(qū)域內(nèi),隨著蓄積量的增加,林分枯損株數(shù)也隨之增加,在其他區(qū)域內(nèi)林分單位面積蓄積量與枯損株數(shù)呈正相關關系??偟膩碚f,在局域模型下各參數(shù)分級的細化程度較為清晰,同時形成了較好范圍內(nèi)的局域化空間分布效果。

    3.2 模型擬合優(yōu)度比較

    根據(jù)兩種模型的均方誤差和赤池信息準則值的計算結果可知,全局泊松模型赤池信息準則值為702.14,局域泊松模型的赤池信息準則值為192.28,說明局域模型在擬合數(shù)據(jù)的優(yōu)良性上產(chǎn)生了較好的效果;同時局域模型的均方誤差為14.64,全局泊松模型的均方誤差值為37.28,說明局域模型的擬合效果要好于全局模型。

    3.3 空間自相關性分析

    3.3.1 模型殘差的全局空間自相關性

    圖2顯示了處于不同空間滯后距離(Lag)狀態(tài)下,全局泊松模型和局域模型殘差的莫蘭指數(shù)(MoranI)相關圖。從總體上看,在空間滯后距離大致相同的情況下,全局泊松模型的莫蘭指數(shù)值要明顯大于局域模型,表明在大部分情況下,全局泊松模型的全局空間相關性要明顯高于局域模型;隨著空間滯后距離的增加,全局泊松模型的莫蘭指數(shù)值逐漸減小,而局域模型的莫蘭指數(shù)值基本在0點上下浮動,起伏沒有全局泊松模型明顯,說明在降低模型殘差全局空間自相關性的效果上,局域模型要明顯好于全局泊松模型。此外,全局泊松模型(GP)在空間滯后距離為14 km時,產(chǎn)生最小的莫蘭指數(shù)值,局域模型(GWPR)在空間滯后距離為5 km時產(chǎn)生最小的莫蘭指數(shù)值。

    圖2 全局泊松模型和局域模型殘差的莫蘭指數(shù)相關圖

    3.3.2 模型殘差的局域空間自相關性

    圖3給出了全局泊松模型和局域模型殘差的莫蘭指數(shù)值的空間分布圖。

    由圖3A可知,在本研究區(qū)域的中部地區(qū),全局泊松模型產(chǎn)生了較大范圍且取值為正值的局域莫蘭指數(shù)值,說明模型殘差在該區(qū)域內(nèi)傾向于出現(xiàn)相同觀測值的聚類模式;在圖3B中,局域模型則在相同區(qū)域內(nèi)產(chǎn)生了取值為負且范圍較小的局域莫蘭指數(shù)值,說明在局域模型下模型殘差出現(xiàn)了不同觀測值的聚類模式,不同觀測值的聚類模式可以很好的揭示空間異質性,產(chǎn)生了較為理想的效果。

    圖3 不同模型殘差的局域莫蘭指數(shù)的空間分布

    圖4給出了在α=0.05顯著性水平下,全局泊松模型和局域模型殘差局域莫蘭指數(shù)的Z值空間分布。GWPR模型(圖4B)的局域Z值在大部分情況下均在-1.96~1.96,在P>0.05的情況下,各樣地點間無顯著性差異;全局泊松模型(圖4A)的Z值呈現(xiàn)顯著性的區(qū)域,說明局域模型在降低局域空間自相關性的效果上要明顯好于全局泊松模型。

    圖4 不同模型殘差的局域莫蘭指數(shù)的Z值空間分布

    3.4 枯損株數(shù)的空間分布

    由圖5可知,從總體上看,兩種模型在對天然次生林枯損株數(shù)空間分布的模擬上與樣地內(nèi)實際枯損株數(shù)的空間分布大致相同,即研究區(qū)的東北部地區(qū)林木枯損株數(shù)較多,在該地區(qū)內(nèi)兩種模型的林木枯損株數(shù)均在15株以上;在研究區(qū)的東南部地區(qū),全局泊松模型的林木枯損株數(shù)在10~15株,而對應局域模型的枯損株數(shù)為8~10株;此外,兩種模型的擬合偏差存在差異,全局泊松模型對樣地內(nèi)枯損株數(shù)模擬的擬合偏差在6株左右,局域模型的擬合偏差在3~4株,要明顯低于全局泊松模型,同時在分級的細化程度上,局域模型也要好于全局泊松模型。

    圖5 不同模型林木枯損量的空間分布

    4 結論與討論

    采用全局泊松模型和局域模型(GWPR)對帽兒山實驗林場101塊天然次生林樣地內(nèi)的林木枯損株數(shù)進行了模擬,并對2種模型的擬合效果以及空間自相關性進行了分析。對于全局泊松模型,林分平均胸徑與次生林枯損株數(shù)呈負相關關系,且影響最大,對于局域模型來說,隨著地理位置的不同,模型參數(shù)產(chǎn)生一定的變化范圍,各個參數(shù)對于枯損株數(shù)的影響大小也隨之變化;在模型擬合方面,GWPR模型擁有很高的擬合精度,AIC值明顯減小,得到了較好的模型參數(shù)的局域化空間分布效果;在解決空間自相關性上,GWPR模型殘差無論是在全局空間自相關性還是局域空間自相關性上都要小于全局模型,同時產(chǎn)生了不同觀測值少量聚類這一理想的空間分布模式;在對枯損株數(shù)空間分布的模擬上,局域模型的擬合偏差也要小于全局模型。

    在林業(yè)數(shù)據(jù)研究中,一般地勢地形不同都會導致研究的數(shù)據(jù)之間存在空間差異,因此空間自相關性在有關林業(yè)的各項研究中普遍存在,同時,泊松模型屬于計數(shù)模型的范疇,在以往的研究中將GWR模型應用到林分枯損株數(shù)空間分布中的研究還未見報道。通過2個模型的空間自相關性分析發(fā)現(xiàn),全局泊松模型的空間自相關性隨著空間滯后距離的增加逐漸下降,局域模型的莫蘭指數(shù)值在0點上下浮動,這是由于局域模型在每塊樣地內(nèi)都可以得到一組模型參數(shù),每組模型參數(shù)都會形成一個GWPR模型,最后綜合成一個模型,局域模型在運行過程中充分考慮了地形因子變化對林木枯損株數(shù)的影響,從根本上解決了空間異質性的問題,而全局泊松模型的所有樣地只有一組模型參數(shù),在可視化分布圖的平滑性以及精度上要小于局域模型(GWPR)。

    有關研究結果表明,選取不同的空間尺度可能對空間自相關性的大小產(chǎn)生一定的影響[17-18],因此,在以后的研究中,應適當結合研究區(qū)的面積,選取不同空間尺度下的局域模型,對天然次生林的枯損株數(shù)進行探索。此外,本文是以天然次生林為一個整體進行研究的,所選取的天然次生林樣地大多為軟闊林、硬闊混交林以及軟闊混交林。由于天然次生林樹種繁多,因此,不同的森林類型可能會導致林分枯損株數(shù)的不同,在今后的研究中,可以按照不同林型對天然次生林枯損株數(shù)的空間分布進行探討。

    [1] 李克志.小興安嶺帶嶺林區(qū)天然林材積枯損率的初步研究[J].林業(yè)科學,1958(1):97-111.

    [2] 曾偉生.西藏天然云杉林枯損率與采伐率模型系統(tǒng)研究[J].林業(yè)科學研究,2008,21(3):353-356.

    [3] 張雄清,雷淵才,雷相東,等.基于計數(shù)模型方法的林分枯損研究[J].林業(yè)科學,2012,48(8):54-61.

    [4] FOTHERINGHAM A S, CHARLTON M, BRUNSDON C. The geography of parameter space: an investigation of spatial non-stationarity[J]. Geographical Information Systems,1996,10(5):605-627.

    [5] BRUNSDON C, FOTHERINGHAM A S, CHARLTON M E. Geographically Weighted Regression: A method for exploring spatial nonstationarity[J]. Geographical Analysis,1996,28(4):281-298.

    [6] 劉暢.黑龍江省森林碳儲量空間分布研究[D].哈爾濱:東北林業(yè)大學,2014.

    [7] 劉暢,李鳳日,甄貞.空間誤差模型在黑龍江省森林碳儲量空間分布的應用[J].應用生態(tài)學報,2014,25(10):2779-2786.

    [8] 戚玉嬌.大興安嶺森林地上碳儲量遙感估算與分析[D].哈爾濱:東北林業(yè)大學,2014.

    [9] 盧軍.帽兒山天然次生林樹冠結構和空間優(yōu)化經(jīng)營[D].哈爾濱:東北林業(yè)大學,2008.

    [10] 王濟川.Logistic回歸模型:方法與應用[M].北京:高等教育出版社,2001.

    [12] STEPHENSON N L, FRANKLIN J F. Causes and implications of the correlation between forest productivity and tree mortality rates[J]. Ecological Monographs,2011,81(4):527-555.

    [13] 王岳,王海燕,楊小娟,等.基于土壤因子的近天然落葉松云冷杉林枯損模型研究[J].福建農(nóng)林大學學報(自然科學版),2015,44(4):379-383.

    [14] 劉鑫,畢華興,李笑吟,等.晉西黃土區(qū)基于地形因子的土壤水分分異規(guī)律研究[J].土壤學報,2007,44(3):411-417.

    [15] 劉暢,李鳳日,賈煒瑋,等.基于局域統(tǒng)計量的黑龍江省多尺度森林碳儲量空間分布變化[J].應用生態(tài)學報,2014,25(9):2493-2500.

    [16] 宋喜芳,李建平,胡希遠.模型選擇信息量準則AIC及其在方差分析中的應用[J].西北農(nóng)林科技大學學報(自然科學版),2009,37(2):88-92.

    [17] MA Z H, BENJAMI Z, WILLIAM F, et al. Use of localized descriptive statistics for exploring the spatial pattern changes of bird species richness at multiple scales[J]. Applied Geography,2012,32(2):185-194.

    [18] GUO L G, MA Z M, ZHANG L Z. Comparison of bandwidth selection in application of geographically weighted regression: A case study[J]. Canadian Journal of Forest Research,2008,38(9):2526-2534.

    猜你喜歡
    莫蘭泊松次生林
    張衛(wèi)《向莫蘭迪致敬之5》
    藝術學研究(2022年6期)2023-01-16 05:48:04
    基于泊松對相關的偽隨機數(shù)發(fā)生器的統(tǒng)計測試方法
    帶有雙臨界項的薛定諤-泊松系統(tǒng)非平凡解的存在性
    淺論北方次生林培育技術要點
    張衛(wèi)作品
    畫刊(2017年9期)2017-09-25 09:49:04
    次生林冠下紅松更新對光環(huán)境適應的研究進展
    阜新地區(qū)天然次生林食葉害蟲無公害防治應用
    張衛(wèi)作品
    畫刊(2017年2期)2017-04-12 06:43:58
    北方次生林培育技術措施探討
    泊松著色代數(shù)
    国内少妇人妻偷人精品xxx网站| 国产69精品久久久久777片| 国产不卡一卡二| av在线观看视频网站免费| 汤姆久久久久久久影院中文字幕 | 别揉我奶头 嗯啊视频| 国产91av在线免费观看| 亚洲av电影不卡..在线观看| 国产单亲对白刺激| 亚洲久久久久久中文字幕| 久久精品熟女亚洲av麻豆精品 | 又粗又爽又猛毛片免费看| 国产成人免费观看mmmm| 国产亚洲精品久久久com| 国产久久久一区二区三区| 国产探花极品一区二区| 欧美日韩在线观看h| 日本猛色少妇xxxxx猛交久久| 日本猛色少妇xxxxx猛交久久| 国产精品久久视频播放| 热99re8久久精品国产| 中文字幕免费在线视频6| 欧美激情国产日韩精品一区| 日本午夜av视频| 99九九线精品视频在线观看视频| 国产精品久久久久久久久免| 少妇高潮的动态图| 久久久久久伊人网av| 91精品伊人久久大香线蕉| 成人亚洲精品av一区二区| 尾随美女入室| 欧美又色又爽又黄视频| 我的女老师完整版在线观看| 3wmmmm亚洲av在线观看| 嫩草影院新地址| 成人午夜高清在线视频| 床上黄色一级片| av卡一久久| 91精品伊人久久大香线蕉| 国产精品嫩草影院av在线观看| 如何舔出高潮| 联通29元200g的流量卡| av在线亚洲专区| 国产av码专区亚洲av| 蜜臀久久99精品久久宅男| 男女下面进入的视频免费午夜| 国产国拍精品亚洲av在线观看| 一级二级三级毛片免费看| 国产精华一区二区三区| 能在线免费观看的黄片| 精品久久久久久久久久久久久| 中文字幕久久专区| 伦精品一区二区三区| 亚洲在久久综合| 日韩欧美在线乱码| 国产精品乱码一区二三区的特点| 白带黄色成豆腐渣| 精品国产露脸久久av麻豆 | 又爽又黄a免费视频| 亚洲三级黄色毛片| 中文字幕av在线有码专区| 国产成人freesex在线| 伦理电影大哥的女人| 在线天堂最新版资源| 日韩精品青青久久久久久| 91aial.com中文字幕在线观看| 国产白丝娇喘喷水9色精品| 成人鲁丝片一二三区免费| 简卡轻食公司| 国产69精品久久久久777片| 国产日韩欧美在线精品| 身体一侧抽搐| 一级毛片电影观看 | av福利片在线观看| 亚洲av免费在线观看| 久久午夜福利片| 久久精品国产鲁丝片午夜精品| av天堂中文字幕网| 99热网站在线观看| 能在线免费观看的黄片| 久久鲁丝午夜福利片| 色噜噜av男人的天堂激情| 亚洲五月天丁香| 精品久久久久久久久久久久久| 亚洲国产高清在线一区二区三| 97超视频在线观看视频| 最近的中文字幕免费完整| 少妇的逼水好多| 波多野结衣高清无吗| 久久久精品欧美日韩精品| 青春草国产在线视频| 欧美日韩国产亚洲二区| 亚洲欧美精品专区久久| 精品久久久久久久末码| 精品国内亚洲2022精品成人| 国产精品福利在线免费观看| 超碰av人人做人人爽久久| 如何舔出高潮| 国内精品美女久久久久久| 亚洲国产欧美在线一区| a级毛色黄片| 日韩国内少妇激情av| av国产免费在线观看| 欧美zozozo另类| 亚洲精品乱码久久久v下载方式| 最近中文字幕2019免费版| 九九久久精品国产亚洲av麻豆| 亚洲高清免费不卡视频| 亚洲在线观看片| 午夜福利高清视频| 两个人视频免费观看高清| 午夜福利成人在线免费观看| 不卡视频在线观看欧美| 男人舔女人下体高潮全视频| 2021少妇久久久久久久久久久| 成人一区二区视频在线观看| 国产三级中文精品| 日韩av在线大香蕉| 亚洲三级黄色毛片| 国产一区二区在线观看日韩| 18禁裸乳无遮挡免费网站照片| 亚洲在线自拍视频| 亚洲精品日韩在线中文字幕| 卡戴珊不雅视频在线播放| av黄色大香蕉| 日韩av在线免费看完整版不卡| 一级二级三级毛片免费看| 午夜精品国产一区二区电影 | 久久久a久久爽久久v久久| 午夜福利高清视频| av在线播放精品| 中文字幕人妻熟人妻熟丝袜美| 天堂中文最新版在线下载 | 亚洲精品亚洲一区二区| 精品免费久久久久久久清纯| 国产日韩欧美在线精品| 免费观看a级毛片全部| 国产成人a∨麻豆精品| 级片在线观看| 波多野结衣高清无吗| 91在线精品国自产拍蜜月| 在线播放国产精品三级| 男人舔奶头视频| 日本免费a在线| 精品99又大又爽又粗少妇毛片| 欧美区成人在线视频| 天堂√8在线中文| 桃色一区二区三区在线观看| 成人亚洲欧美一区二区av| 亚洲国产欧美在线一区| 亚洲欧美精品综合久久99| 老司机影院毛片| 久久人人爽人人爽人人片va| 久久久国产成人精品二区| 网址你懂的国产日韩在线| 深夜a级毛片| 99久久无色码亚洲精品果冻| 国产成人精品一,二区| 日韩欧美 国产精品| 欧美性感艳星| 亚洲图色成人| 午夜日本视频在线| 亚洲欧美清纯卡通| 亚洲成av人片在线播放无| 男人舔女人下体高潮全视频| 中文字幕精品亚洲无线码一区| 国产免费男女视频| 亚洲综合色惰| 国产单亲对白刺激| 22中文网久久字幕| 六月丁香七月| 欧美日本视频| 在线免费观看的www视频| 最近2019中文字幕mv第一页| 亚洲av二区三区四区| 国产成年人精品一区二区| 最近的中文字幕免费完整| 大香蕉久久网| 久久精品综合一区二区三区| av黄色大香蕉| 亚洲欧美日韩无卡精品| 日韩亚洲欧美综合| 精品酒店卫生间| 国产成年人精品一区二区| 韩国高清视频一区二区三区| .国产精品久久| 九色成人免费人妻av| 国产一级毛片在线| 三级国产精品欧美在线观看| 久久综合国产亚洲精品| 欧美一级a爱片免费观看看| 只有这里有精品99| 麻豆成人午夜福利视频| 欧美不卡视频在线免费观看| 亚洲国产欧美人成| 国产 一区 欧美 日韩| 又粗又爽又猛毛片免费看| 免费搜索国产男女视频| 大又大粗又爽又黄少妇毛片口| 亚洲欧美精品专区久久| 精品一区二区三区视频在线| 成年版毛片免费区| 一级毛片久久久久久久久女| 国产精品电影一区二区三区| 97人妻精品一区二区三区麻豆| 99久国产av精品国产电影| 久久久久久久久久成人| 99在线人妻在线中文字幕| 中文字幕制服av| 大话2 男鬼变身卡| videossex国产| 人人妻人人看人人澡| 久久婷婷人人爽人人干人人爱| 亚洲国产色片| 一级黄色大片毛片| 青青草视频在线视频观看| 久久久久久久亚洲中文字幕| 少妇人妻一区二区三区视频| 在线免费观看不下载黄p国产| 久久久久久久久久久免费av| 欧美性感艳星| 国产激情偷乱视频一区二区| av免费在线看不卡| 91精品国产九色| 日本黄大片高清| 视频中文字幕在线观看| 麻豆一二三区av精品| 热99在线观看视频| 国产亚洲一区二区精品| 日韩欧美 国产精品| 久久韩国三级中文字幕| 国产精品野战在线观看| 丰满人妻一区二区三区视频av| 又爽又黄无遮挡网站| 成年女人看的毛片在线观看| 毛片一级片免费看久久久久| 内射极品少妇av片p| 亚洲av福利一区| 老司机福利观看| 91久久精品电影网| 春色校园在线视频观看| 两性午夜刺激爽爽歪歪视频在线观看| 成人特级av手机在线观看| 国产精品永久免费网站| 日本午夜av视频| 国产美女午夜福利| 午夜福利在线观看免费完整高清在| 久久99热这里只频精品6学生 | 久久久久久伊人网av| 欧美3d第一页| 国产极品天堂在线| 久久久久网色| 精品国产露脸久久av麻豆 | 毛片女人毛片| 国产精品一区二区在线观看99 | 亚洲精品影视一区二区三区av| kizo精华| 成人特级av手机在线观看| 老师上课跳d突然被开到最大视频| 黄片无遮挡物在线观看| 久99久视频精品免费| .国产精品久久| 青青草视频在线视频观看| 久久精品熟女亚洲av麻豆精品 | 久久久久久大精品| 夜夜爽夜夜爽视频| 国产成人免费观看mmmm| 久久精品国产鲁丝片午夜精品| 村上凉子中文字幕在线| 最近的中文字幕免费完整| 日本黄色视频三级网站网址| 午夜老司机福利剧场| 久久久久久久久久成人| 偷拍熟女少妇极品色| 日韩视频在线欧美| 国产精品,欧美在线| 一卡2卡三卡四卡精品乱码亚洲| 男人舔奶头视频| 国模一区二区三区四区视频| 亚洲丝袜综合中文字幕| 白带黄色成豆腐渣| 亚洲人成网站在线观看播放| 一区二区三区高清视频在线| av播播在线观看一区| 嫩草影院新地址| 国产91av在线免费观看| 日韩在线高清观看一区二区三区| 中文欧美无线码| 久久6这里有精品| 国产又黄又爽又无遮挡在线| ponron亚洲| or卡值多少钱| 亚洲精品456在线播放app| 99热这里只有是精品50| 永久网站在线| 热99在线观看视频| av在线天堂中文字幕| 欧美成人免费av一区二区三区| 全区人妻精品视频| 久久久久久久午夜电影| 特大巨黑吊av在线直播| av.在线天堂| 我要看日韩黄色一级片| 内地一区二区视频在线| 亚洲电影在线观看av| 激情 狠狠 欧美| 国产在线一区二区三区精 | 精品免费久久久久久久清纯| 美女被艹到高潮喷水动态| 欧美日韩精品成人综合77777| 日韩 亚洲 欧美在线| .国产精品久久| 免费观看人在逋| 国产激情偷乱视频一区二区| 七月丁香在线播放| 晚上一个人看的免费电影| 亚洲精品亚洲一区二区| 国产亚洲av嫩草精品影院| 国产伦理片在线播放av一区| 日日撸夜夜添| 亚洲国产色片| 亚洲av电影不卡..在线观看| 干丝袜人妻中文字幕| 午夜免费男女啪啪视频观看| 精华霜和精华液先用哪个| 亚洲国产最新在线播放| 99热这里只有是精品50| 国产三级中文精品| 99久久无色码亚洲精品果冻| 直男gayav资源| 边亲边吃奶的免费视频| 日韩成人av中文字幕在线观看| 国产成人a区在线观看| 亚洲成人中文字幕在线播放| 国产淫片久久久久久久久| 国产又色又爽无遮挡免| 国产一区二区在线观看日韩| 日本欧美国产在线视频| 大香蕉久久网| 超碰av人人做人人爽久久| 天堂av国产一区二区熟女人妻| 日日撸夜夜添| 中文字幕免费在线视频6| 欧美日本亚洲视频在线播放| 国产老妇伦熟女老妇高清| 国产高清不卡午夜福利| 日日干狠狠操夜夜爽| 国产亚洲精品av在线| 亚洲欧洲日产国产| 99热这里只有是精品50| 99久久精品热视频| 男人狂女人下面高潮的视频| 波多野结衣高清无吗| av女优亚洲男人天堂| av.在线天堂| 18禁在线无遮挡免费观看视频| 熟妇人妻久久中文字幕3abv| 淫秽高清视频在线观看| 精品国内亚洲2022精品成人| 亚洲av成人精品一二三区| 在线免费观看不下载黄p国产| 日本黄大片高清| 69人妻影院| 久久这里有精品视频免费| 五月伊人婷婷丁香| 国产精华一区二区三区| 少妇高潮的动态图| 免费观看人在逋| 国产在线一区二区三区精 | 精品久久久久久电影网 | 国产成人精品久久久久久| 精品免费久久久久久久清纯| 国产成人a区在线观看| 精品人妻视频免费看| 97热精品久久久久久| 天堂av国产一区二区熟女人妻| 国产 一区精品| 能在线免费看毛片的网站| 国产在线男女| 丰满人妻一区二区三区视频av| 国产一级毛片在线| 寂寞人妻少妇视频99o| 高清av免费在线| 国产精品国产三级国产专区5o | av在线播放精品| 国产伦精品一区二区三区四那| 日本wwww免费看| 韩国av在线不卡| 日本黄色视频三级网站网址| 日日摸夜夜添夜夜爱| 小蜜桃在线观看免费完整版高清| 亚洲自偷自拍三级| 国产精品蜜桃在线观看| 波多野结衣高清无吗| 亚洲va在线va天堂va国产| 精品少妇黑人巨大在线播放 | 大又大粗又爽又黄少妇毛片口| 日本午夜av视频| 51国产日韩欧美| 国产精品一区二区三区四区免费观看| 色哟哟·www| 99久久成人亚洲精品观看| 日韩精品青青久久久久久| 一本久久精品| 九草在线视频观看| 我要看日韩黄色一级片| 联通29元200g的流量卡| 亚洲伊人久久精品综合 | 黄色一级大片看看| 91久久精品国产一区二区三区| 亚洲天堂国产精品一区在线| 中文欧美无线码| 欧美xxxx性猛交bbbb| 在线天堂最新版资源| 国产免费视频播放在线视频 | 三级经典国产精品| 一个人观看的视频www高清免费观看| 国产白丝娇喘喷水9色精品| 能在线免费看毛片的网站| 亚洲第一区二区三区不卡| 老司机影院成人| 天天躁夜夜躁狠狠久久av| 99热这里只有是精品在线观看| 97人妻精品一区二区三区麻豆| 亚洲久久久久久中文字幕| 七月丁香在线播放| 99久久精品一区二区三区| 欧美性猛交黑人性爽| 日韩成人伦理影院| 一个人观看的视频www高清免费观看| 国产老妇女一区| 18禁裸乳无遮挡免费网站照片| 我的老师免费观看完整版| 人妻少妇偷人精品九色| 欧美日韩在线观看h| 欧美高清性xxxxhd video| 美女黄网站色视频| 小蜜桃在线观看免费完整版高清| 久久精品久久精品一区二区三区| 欧美区成人在线视频| 亚洲国产色片| 黄片无遮挡物在线观看| 国产在线一区二区三区精 | 精品久久久久久久久久久久久| 亚洲自偷自拍三级| 91精品一卡2卡3卡4卡| 青青草视频在线视频观看| 久久久久久久久久久免费av| 日韩亚洲欧美综合| 亚洲乱码一区二区免费版| 精品99又大又爽又粗少妇毛片| 免费电影在线观看免费观看| 国产在线男女| 欧美激情在线99| av天堂中文字幕网| 国产精品一区www在线观看| 黄片无遮挡物在线观看| av黄色大香蕉| 一级毛片电影观看 | 嫩草影院精品99| 日韩国内少妇激情av| 一边摸一边抽搐一进一小说| 天天躁夜夜躁狠狠久久av| 国产单亲对白刺激| 色吧在线观看| kizo精华| 18禁在线无遮挡免费观看视频| 国产免费福利视频在线观看| 日本黄色片子视频| 99热网站在线观看| 小蜜桃在线观看免费完整版高清| 欧美激情在线99| 欧美三级亚洲精品| 国产精品国产三级国产专区5o | 国产精品一及| 精品久久久噜噜| 最近视频中文字幕2019在线8| 大香蕉久久网| 1024手机看黄色片| 欧美日韩在线观看h| 国产乱人偷精品视频| 国产黄a三级三级三级人| 熟妇人妻久久中文字幕3abv| 最近手机中文字幕大全| 91午夜精品亚洲一区二区三区| 国产精品伦人一区二区| 我的女老师完整版在线观看| 免费av观看视频| 毛片一级片免费看久久久久| ponron亚洲| 国产v大片淫在线免费观看| 美女国产视频在线观看| 久久久久久久久中文| av在线播放精品| 国产一级毛片七仙女欲春2| 99热6这里只有精品| 亚洲精品成人久久久久久| 日本一本二区三区精品| 18禁在线播放成人免费| 一级黄色大片毛片| 久久精品夜色国产| 婷婷六月久久综合丁香| 午夜视频国产福利| 性插视频无遮挡在线免费观看| 国产精品一区二区性色av| 又爽又黄无遮挡网站| 日韩大片免费观看网站 | 国产一区亚洲一区在线观看| 18+在线观看网站| 六月丁香七月| 又粗又爽又猛毛片免费看| 亚洲久久久久久中文字幕| 变态另类丝袜制服| 天天一区二区日本电影三级| 在现免费观看毛片| 哪个播放器可以免费观看大片| 国产亚洲av片在线观看秒播厂 | 卡戴珊不雅视频在线播放| 1000部很黄的大片| 国产一级毛片在线| 又粗又硬又长又爽又黄的视频| 午夜免费男女啪啪视频观看| 一区二区三区免费毛片| 国产精品国产三级专区第一集| 亚洲人成网站在线播| АⅤ资源中文在线天堂| www.av在线官网国产| 久久久国产成人免费| 久久国产乱子免费精品| 一级爰片在线观看| 国产黄片视频在线免费观看| 我的老师免费观看完整版| 亚洲欧美精品综合久久99| 亚洲av成人av| 黄色一级大片看看| 国产成人aa在线观看| 亚洲欧洲国产日韩| 午夜免费男女啪啪视频观看| 丰满少妇做爰视频| 欧美日本亚洲视频在线播放| 国产高清三级在线| 九色成人免费人妻av| kizo精华| av视频在线观看入口| 国产精品蜜桃在线观看| 人体艺术视频欧美日本| 两个人视频免费观看高清| 久久久久免费精品人妻一区二区| 搡老妇女老女人老熟妇| 欧美不卡视频在线免费观看| 久久午夜福利片| 九九在线视频观看精品| 毛片女人毛片| 亚洲精品日韩在线中文字幕| 麻豆乱淫一区二区| 日本免费一区二区三区高清不卡| 久久人人爽人人片av| av播播在线观看一区| 中文字幕制服av| 精品国产三级普通话版| 久久久久久久久久久丰满| 黑人高潮一二区| 国产精品野战在线观看| 中国国产av一级| 日韩,欧美,国产一区二区三区 | 国产激情偷乱视频一区二区| 中国国产av一级| 大香蕉久久网| 国产精华一区二区三区| 天天躁夜夜躁狠狠久久av| 成人鲁丝片一二三区免费| 热99在线观看视频| 99热这里只有是精品50| 九九热线精品视视频播放| 91精品一卡2卡3卡4卡| 久久6这里有精品| 全区人妻精品视频| 欧美xxxx性猛交bbbb| 国产 一区精品| 级片在线观看| 有码 亚洲区| 搞女人的毛片| 亚洲高清免费不卡视频| 三级国产精品欧美在线观看| 99热全是精品| 男人舔女人下体高潮全视频| 男女那种视频在线观看| 成人特级av手机在线观看| 日本熟妇午夜| 亚洲欧美清纯卡通| 久久欧美精品欧美久久欧美| 老司机福利观看| 国产不卡一卡二| 赤兔流量卡办理| 能在线免费看毛片的网站| 亚洲国产精品sss在线观看| 联通29元200g的流量卡| videossex国产| 日本色播在线视频| 少妇的逼水好多| 国内精品宾馆在线| 国产不卡一卡二| 久久亚洲国产成人精品v| 欧美色视频一区免费| www.色视频.com| 天天一区二区日本电影三级| 少妇人妻精品综合一区二区| 亚洲av免费在线观看| 久久久久久久久久久免费av| 色视频www国产| 亚洲熟妇中文字幕五十中出| 波野结衣二区三区在线| 国产精品一二三区在线看| 国产成人精品婷婷| 黄色日韩在线|