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

    油松林分?jǐn)嗝娣e生長(zhǎng)預(yù)估模型研究

    2015-03-21 07:30:27高東啟鄧華鋒程志楚
    關(guān)鍵詞:株數(shù)間伐油松

    高東啟 鄧華鋒 蔣 益 程志楚

    (1.北京林業(yè)大學(xué)林學(xué)院,北京 100083;2.國(guó)家林業(yè)局調(diào)查規(guī)劃設(shè)計(jì)院,北京 100714)

    ?

    油松林分?jǐn)嗝娣e生長(zhǎng)預(yù)估模型研究

    高東啟1鄧華鋒1蔣 益1程志楚2

    (1.北京林業(yè)大學(xué)林學(xué)院,北京 100083;2.國(guó)家林業(yè)局調(diào)查規(guī)劃設(shè)計(jì)院,北京 100714)

    在森林經(jīng)營(yíng)中,間伐林分、未間伐林分分開單獨(dú)建模存在不相容的問(wèn)題,有必要對(duì)其進(jìn)行研究,并建立統(tǒng)一的模型。利用森林資源一類清查數(shù)據(jù),以Schumacher模型為基礎(chǔ),通過(guò)引入啞變量、間伐指標(biāo)分別建立北京市油松林分?jǐn)嗝娣e生長(zhǎng)預(yù)估模型,使得間伐林分和未間伐林分能夠整合在一起建立統(tǒng)一的模型。經(jīng)檢驗(yàn),引入啞變量模型對(duì)油松林分?jǐn)嗝娣e的預(yù)測(cè)精度最低,為87.99%;其次是引入間伐指標(biāo)模型,預(yù)估精度為88.36%;綜合引入啞變量和間伐指標(biāo)模型預(yù)估精度最高,為88.61%。所建模型對(duì)間伐林分和未間伐林分均適用,更方便于在實(shí)踐中應(yīng)用。

    林分?jǐn)嗝娣e;啞變量;間伐指標(biāo);油松

    油松(Pinustabulaeformis)是北京市主要的綠化樹種,多為人工林,天然林較少,在保護(hù)首都生態(tài)環(huán)境安全方面發(fā)揮著重要作用。據(jù)第七次全國(guó)森林資源清查統(tǒng)計(jì)[1],油松占北京市森林面積總量的18.7%,占蓄積總量的11.9%。為了提高林分的質(zhì)量和效益、科學(xué)合理的經(jīng)營(yíng)管理森林資源,有必要通過(guò)生長(zhǎng)模擬來(lái)掌握林分的生長(zhǎng)狀況?,F(xiàn)實(shí)林分由于受間伐措施的影響,使得間伐林分的生長(zhǎng)規(guī)律、空間結(jié)構(gòu)和生態(tài)環(huán)境相對(duì)未間伐林分都可能發(fā)生變化[2-4]。從理論上講,同一樹種的間伐林分與未間伐林分需要分開單獨(dú)建立模型,而間伐林分與未間伐林分之間是存在一定聯(lián)系的,因此如何將這兩種類型的林分整合在一起建立統(tǒng)一的模型是當(dāng)前急需解決的問(wèn)題,同時(shí)也使得間伐林分與未間伐林分生長(zhǎng)模型能夠緊密聯(lián)系形成一個(gè)有機(jī)的整體。關(guān)于間伐林分的研究,國(guó)外的Bailey等[5]在建立林分?jǐn)嗝娣e生長(zhǎng)預(yù)估模型時(shí)加入了一個(gè)表示間伐的自變量,Murray等[6]考慮用分布函數(shù)來(lái)表示間伐對(duì)林分?jǐn)嗝娣e生長(zhǎng)的影響,Pienaar等[7]提出了由未間伐林分預(yù)估間伐林分的新思路;國(guó)內(nèi)的杜紀(jì)山、李春明等研究了間伐對(duì)杉木(Cunninghamialarceolata)、落葉松(Larixgmelini)人工林生長(zhǎng)的影響[8-10]。

    本研究以油松人工林為例,從相容的角度出發(fā),以Schumacher模型為基礎(chǔ),分別引入啞變量[11-12]、間伐指標(biāo)[7],建立間伐林分和未間伐林分相容的、統(tǒng)一的油松林分?jǐn)嗝娣e生長(zhǎng)預(yù)估模型,并對(duì)模型的預(yù)測(cè)能力進(jìn)行分析比較。

    1 研究區(qū)概況

    北京市位于北緯39°28′~41°05′,東經(jīng)115°25′~117°30′,地處華北平原北端,北以燕山山地與內(nèi)蒙古高原接壤,西以太行山與河北省相連,東北與松遼平原相通,南與黃淮海平原連片;屬于暖溫帶半濕潤(rùn)大陸性季風(fēng)型氣候,四季分明,夏季炎熱多雨,冬季寒冷干燥,夏季降水量占全年降水量的74%。

    2 研究方法

    2.1 樣地調(diào)查及數(shù)據(jù)整理

    本研究所采用的數(shù)據(jù)為北京市一類清查數(shù)據(jù),每個(gè)樣地的面積為 0.066 7 hm2。樣地調(diào)查因子有:林木胸徑、林分年齡、林分平均高、林分蓄積、采伐蓄積等。所使用的數(shù)據(jù)分別調(diào)查于1996、2001、2006年,3期數(shù)據(jù)共計(jì)組成99個(gè)樣本,其中,間伐林分56個(gè),未間伐林分43個(gè),部分樣地在調(diào)查間隔期內(nèi)進(jìn)行了間伐,間伐強(qiáng)度不固定,間伐間隔期5 a。間伐林分與未間伐林分樣地的基本情況見表1。主要統(tǒng)計(jì)了林分的平均年齡、平均胸徑、平均樹高、單位面積株數(shù)、單位面積斷面積和單位面積蓄積量。

    表1 間伐林分和未間伐林分主要因子統(tǒng)計(jì)

    從表1中可以看出,油松間伐林分的平均年齡、平均胸徑、平均樹高、單位面積斷面積和單位面積蓄積量要明顯大于未間伐林分,而單位面積株數(shù)又小于未間伐林分,說(shuō)明間伐效果明顯。

    對(duì)99個(gè)樣本進(jìn)行隨機(jī)抽樣,其中50個(gè)樣本用于建模,包含間伐林分31個(gè)、未間伐林分19個(gè);49個(gè)樣本用于檢驗(yàn),包含間伐林分25個(gè)、未間伐林分24個(gè),樣地分布情況見表2。

    表2 隨機(jī)抽樣樣地分布情況

    2.2 林分優(yōu)勢(shì)木平均高計(jì)算

    林分優(yōu)勢(shì)木平均高與林分平均高呈線性關(guān)系,由于所用資料中沒(méi)有油松林分優(yōu)勢(shì)木平均高的相關(guān)記錄,參考賀姍姍[13]的研究結(jié)果,其通過(guò)解析優(yōu)勢(shì)木的方法,得到油松人工林優(yōu)勢(shì)木平均高與林分平均高的相關(guān)關(guān)系:

    H=1.115h+0.429 3

    (1)

    式中:H為林分優(yōu)勢(shì)木平均高;h為林分平均樹高;相關(guān)系數(shù)為0.952 1。

    通過(guò)(1)式計(jì)算出各樣地的林分優(yōu)勢(shì)木平均高,作為評(píng)價(jià)林分立地質(zhì)量的指標(biāo)。

    2.3 基于啞變量的方法

    通過(guò)引入啞變量[11]將油松間伐林分和未間伐林分2個(gè)類型的林分用定性代碼來(lái)表示,具體過(guò)程是將第i個(gè)類型的林分編號(hào)為Ki,將定性數(shù)據(jù)Ki轉(zhuǎn)化為(0,1):

    式中:i=1,2;K1、K2分別為間伐林分和未間伐林分的定性代碼。

    為了和下面的方法進(jìn)行對(duì)比,選用Schumacher模型建立林分?jǐn)嗝娣e生長(zhǎng)預(yù)估模型,其基本形式為[7]:

    G=exp(a1+a2/t)Na3+a5/tHa4+a6/t

    (2)

    為了減小建模時(shí)所產(chǎn)生的異方差,可采用對(duì)數(shù)回歸的方法[14],對(duì)Schumacher模型兩邊取對(duì)數(shù)變?yōu)槿缦滦问剑?/p>

    lnG=a1+a2/t+a3ln(N)a4ln(H)+a5ln(N)/t+a6ln(H)/t

    (3)

    參考李忠國(guó)等建立啞變量模型的方法[12],將啞變量引入到模型(3)式的參數(shù)中。經(jīng)過(guò)初步嘗試,發(fā)現(xiàn)在參數(shù)a0~a6中分別引入啞變量時(shí)模型的決定系數(shù)是不同的。其中,在a4中引入啞變量時(shí)模型的決定系數(shù)最高,說(shuō)明表示油松間伐、未間伐林分的啞變量對(duì)參數(shù)a4的影響最大,既將啞變量引入到模型(3)式的a4中最為合理,則確定啞變量模型的最終形式如下:

    lnG=a1+a2/t+a3ln(N)+(a4K1+a5K2) ln(H)+a6ln(N)/t+a7ln(H)/t

    (4)

    式(2)、(3)、(4)中:G為單位面積林分?jǐn)嗝娣e;lnG為林分?jǐn)嗝娣e的對(duì)數(shù)形式;H為林分優(yōu)勢(shì)木平均高;N為林分單位面積株數(shù);t為林分平均年齡;a1~a7為待定參數(shù);K1、K2分別為間伐林分和未間伐林分的定性代碼。

    采用含有啞變量的模型(4)式建立油松間伐林分和未間伐林分相容的林分?jǐn)嗝娣e生長(zhǎng)預(yù)估模型,并估計(jì)各個(gè)參數(shù)值。

    2.4 引入間伐指標(biāo)的方法

    有學(xué)者研究發(fā)現(xiàn)[7],間伐林分與未間伐林分?jǐn)嗝娣e差異的大小取決于間伐時(shí)的年齡和間伐強(qiáng)度,并且得到間伐林分和未間伐林分相容的斷面積預(yù)估方程:

    lnG=a1+a2/t+a3ln(N)+a4ln(N)+a5ln(H)/t+a6ln(H)/t+a7Nftf/(Nat)

    (5)

    式中:G為單位面積林分?jǐn)嗝娣e;lnG為斷面積的對(duì)數(shù)形式;N為林分單位面積株數(shù);H為林分優(yōu)勢(shì)木平均高;t為林分現(xiàn)在的平均年齡;Nf為最近一次間伐中伐去林木株數(shù);tf為林分最近一次間伐時(shí)的年齡;Na為最近一次間伐后保留林木株數(shù);a1~a7為待定參數(shù)。

    模型(5)式相對(duì)于對(duì)數(shù)形式的Schumacher模型(3)式而言,實(shí)際上是在(3)式的基礎(chǔ)上引入了一些表示間伐林分特征的相關(guān)指標(biāo),有表示間伐強(qiáng)度的Nf和Na,還有表示間伐年齡的tf。直接用(5)式去預(yù)估林分的斷面積,由于引入了間伐指標(biāo),也使得間伐林分與未間伐林分在建模時(shí)具有相容性。

    2.5 綜合引入啞變量和間伐指標(biāo)的方法

    在上述引入啞變量或間伐指標(biāo)的2種方法的基礎(chǔ)上,在Schumacher模型(3)式中同時(shí)引入啞變量和間伐指標(biāo),作為預(yù)估間伐林分與未間伐林分?jǐn)嗝娣e的第3種方法,其形式如下:

    lnG=a1+a2/t+a3ln(N)+(a4K1+a5K2)ln(H)+a6ln(N)/t+a7ln(H)/t+a8Nftf/(Nat)

    (6)

    式中:G為單位面積林分?jǐn)嗝娣e;lnG為林分?jǐn)嗝娣e的對(duì)數(shù)形式;N為林分單位面積株數(shù);H為林分優(yōu)勢(shì)木平均高;t為林分現(xiàn)在的平均年齡;Nf為最近一次間伐中伐去林木株數(shù);tf為林分最近一次間伐時(shí)的年齡;Na為最近一次間伐后保留林木株數(shù);a1~a8為待定參數(shù);K1、K2分別為間伐林分和未間伐林分的定性代碼。

    2.6 參數(shù)估計(jì)與模型檢驗(yàn)

    綜合應(yīng)用ForStat、Excel、SPSS進(jìn)行數(shù)據(jù)處理和參數(shù)估計(jì),計(jì)算均方根誤差(RMSE)、總相對(duì)誤差(TRE)、平均系統(tǒng)誤差(MSE)、預(yù)估精度(P)等指標(biāo)來(lái)檢驗(yàn)?zāi)P偷念A(yù)測(cè)能力,同時(shí)對(duì)所建模型進(jìn)行成對(duì)t檢驗(yàn)。

    3 結(jié)果與分析

    統(tǒng)計(jì)各模型的參數(shù)估計(jì)結(jié)果見表3。啞變量模型(4)式的擬合效果最差,決定系數(shù)R2為 0.845 0;其次是引入間伐指標(biāo)的模型(5)式,決定系數(shù)R2為 0.846 2;綜合引入了啞變量和間伐指標(biāo)的模型(6)式擬合效果最好,決定系數(shù)R2為 0.846 9,平均絕對(duì)偏差(MAD)和均方根誤差(RMSE)最小。

    表3 模型參數(shù)統(tǒng)計(jì)

    利用檢驗(yàn)數(shù)據(jù)對(duì)模型(4)、(5)、(6)式進(jìn)行檢驗(yàn),統(tǒng)計(jì)結(jié)果見表4。在表4中,啞變量模型(4)式的各項(xiàng)誤差最大,預(yù)估精度最低,為87.99%;引入間伐指標(biāo)的模型(5)式居中,預(yù)估精度為88.36%;綜合引入啞變量和間伐指標(biāo)的模型(6)式預(yù)估精度最高為88.61%,各項(xiàng)誤差均較小。進(jìn)行成對(duì)t檢驗(yàn),結(jié)果為Sig.<0.05,說(shuō)明在0.05顯著水平上模型的回歸效果顯著,實(shí)測(cè)值與預(yù)測(cè)值無(wú)顯著差異,所建的林分?jǐn)嗝娣e模型比較合理。

    進(jìn)一步檢驗(yàn)?zāi)P蛯?duì)間伐林分、未間伐林分?jǐn)嗝娣e的預(yù)測(cè)能力,統(tǒng)計(jì)結(jié)果見表5。從表5中可以看出:對(duì)于間伐林分而言,模型(4)、(5)、(6)式的預(yù)測(cè)精度分別為91.86%、91.90%、91.96%,預(yù)測(cè)效果較好,預(yù)測(cè)誤差以(4)式最大,(6)式最?。粚?duì)于未間伐林分而言,模型(4)、(5)、(6)式的預(yù)測(cè)精度分別為83.17%、84.47%、85.23%,以(6)式誤差最小,這說(shuō)明綜合引入啞變量和間伐指標(biāo)的模型(6)式要優(yōu)于模型(4)、(5)式,而引入間伐指標(biāo)的模型(5)式對(duì)間伐、未間伐林分的預(yù)測(cè)效果又要好于引入啞變量的模型(4)式。

    表4 模型檢驗(yàn)

    表5 模型對(duì)間伐、未間伐林分的預(yù)測(cè)效果檢驗(yàn)

    4 結(jié)論與討論

    以Schumacher模型為基礎(chǔ),采用3種方法將間伐林分與未間伐林分整合在一起分別建立北京市油松人工林的斷面積生長(zhǎng)預(yù)估模型,經(jīng)檢驗(yàn),模型對(duì)油松林分?jǐn)嗝娣e的預(yù)估效果都較好,預(yù)估精度在87%以上。無(wú)論是引入啞變量或間伐指標(biāo),還是綜合引入啞變量和間伐指標(biāo),都能使得間伐林分和未間伐林分整合在一起建立統(tǒng)一的模型,從而使2種不同類型的林分生長(zhǎng)模型形成一個(gè)有機(jī)聯(lián)系的整體,既適用于間伐林分,也適用于未間伐林分,便于在實(shí)踐中應(yīng)用。

    3種方法由于引入的因子不同,其預(yù)估效果也不相同。引入啞變量的模型對(duì)油松林分?jǐn)嗝娣e的預(yù)測(cè)精度最低,為87.99%;其次是引入間伐指標(biāo)的模型,預(yù)估精度為88.36%;綜合引入啞變量和間伐指標(biāo)的模型預(yù)估精度最高,為88.61%??赡艿脑蚴且腴g伐指標(biāo)增加了林分狀況的有效信息,使模型的預(yù)測(cè)效果更好;而引入啞變量并不能增加林分信息因子,只是調(diào)整了模型結(jié)構(gòu)建立了間伐林分和未間伐林分相統(tǒng)一的模型;綜合引入啞變量和間伐指標(biāo)不僅提高了模型對(duì)林分信息的利用程度,同時(shí)還優(yōu)化了模型結(jié)構(gòu),使其預(yù)測(cè)效果最好。

    由于本研究所用數(shù)據(jù)有限,且時(shí)間間隔較短,下一步對(duì)于油松間伐林分的研究可進(jìn)一步收據(jù)信息更加豐富、內(nèi)容更加全面的長(zhǎng)期觀測(cè)數(shù)據(jù),考慮區(qū)域效應(yīng)、樣地效應(yīng)、時(shí)間效應(yīng)等隨機(jī)因素和其他經(jīng)營(yíng)管理措施對(duì)模擬林分生長(zhǎng)的影響。

    [1] 賈忠奎,馬履一,徐程揚(yáng),等.北京市森林資源動(dòng)態(tài)及可持續(xù)經(jīng)營(yíng)對(duì)策[J].干旱區(qū)資源與環(huán)境,2006,20(3):30-36.

    [2] 李春明,杜紀(jì)山,張會(huì)儒.撫育間伐對(duì)森林生長(zhǎng)的影響及其模型研究[J].林業(yè)科學(xué)研究,2003,16(5):636-641.

    [3] 杜紀(jì)山,唐守正.撫育間伐對(duì)林分生長(zhǎng)的效應(yīng)及其模型研究[J].北京林業(yè)大學(xué)學(xué)報(bào),1996,16(1):80-84.

    [4] 岳永杰,余新曉,李鋼鐵,等.北京松山自然保護(hù)區(qū)蒙古櫟林的空間結(jié)構(gòu)特征[J].應(yīng)用生態(tài)學(xué)報(bào),2009,20(8):1811-1816.

    [5] Bailey R L, Ware K D. Compatible basal area growth and yield model for thinned and unthinned stands[J]. Canadian Journal Forest Research, 1983(13): 563-571.

    [6] Murray D M, Von G K. Relationships between the diameter distributions before and after thinning[J]. Forest Science, 1991, 37(2): 552-559.

    [7] Pienaar L V, Shiver B D. Basal area prediction and projection equations for pine plantations[J]. Forest Science, 1986, 32(3): 626-633.

    [8] 杜紀(jì)山,唐守正.杉木林分?jǐn)嗝娣e生長(zhǎng)預(yù)估模型及其應(yīng)用[J].北京林業(yè)大學(xué)學(xué)報(bào),1998,20(4):4-8.

    [9] 李春明,杜紀(jì)山,張會(huì)儒.間伐林分的斷面積生長(zhǎng)模型研究[J].林業(yè)資源管理,2004(3):52-55.

    [10] 李春明,杜紀(jì)山,張會(huì)儒.撫育間伐對(duì)人工落葉松斷面積和蓄積生長(zhǎng)的影響[J].林業(yè)資源管理,2007(3):90-93.

    [11] 唐守正,郎奎建,李海奎.統(tǒng)計(jì)與生物數(shù)學(xué)模型計(jì)算(ForStat教程) [M].北京:科學(xué)出版社,2009.

    [12] 李忠國(guó),孫曉梅,陳東升,等.基于啞變量的日本落葉松生長(zhǎng)模型研究[J].西北農(nóng)林科技大學(xué)學(xué)報(bào):自然科學(xué)版,2011,39(8):69-74.

    [13] 賀姍姍.北京山區(qū)油松人工林林分結(jié)構(gòu)與生長(zhǎng)模擬研究[D].北京:北京林業(yè)大學(xué),2009.

    [14] 曾偉生,唐守正.非線性模型對(duì)數(shù)回歸的偏差校正及與加權(quán)回歸的對(duì)比分析[J].林業(yè)科學(xué)研究,2011,24(2):137-143.

    (責(zé)任編輯 趙粉俠)

    Forecast Models Research of Stands Basal Area Growth forPinustabulaeformis

    GAO Dong-qi1,DENG Hua-feng1,JIANG Yi1,CHENG Zhi-chu2

    (1.College of Forestry, Beijing Forestry University, Beijing 100083, China; 2.Academy of Forest Inventory and Planning, State Forestry Administration, Beijing 100714, China)

    It was incompatible to model thinned stands and unthinned stands separately, so it was necessary to study and establish the unified models. Using forest inventory data, based on Schumacher model, through introducing the dummy variables and thinning indexes to establish the stands basal area growth forecast models ofPinustabulaeformisin Beijing respectively, and made the thinned stands and unthinned stands together to establish unified models.After inspection,on the forecast effect for stands basal area ofPinustabulaeformis, forecast accuracy of model with dummy variables was 87.99% as a minimum; followed by the model with thinning indexes, forecast accuracy was 88.36%;forecast effect of model with dummy variables and thinning indexes was the best, forecast accuracy was as high as 88.61%.The models was applicable for thinned stands and unthinned stands, it was more convenient for application in practice.

    stands basal area;dummy variables;thinning indexes;Pinustabulaeformis

    2014-04-29

    林業(yè)公益性行業(yè)科研專項(xiàng)(201204510)資助。

    鄧華鋒(1966—),男,博士,教授。研究方向:森林可持續(xù)經(jīng)營(yíng)。Email:denghuafeng@bjfu.edu.cn。

    10.11929/j.issn.2095-1914.2015.01.009

    S758.5

    A

    2095-1914(2015)01-0042-05

    第1作者:高東啟(1986—),男,碩士生。研究方向:森林可持續(xù)經(jīng)營(yíng)。Email:511602746@qq.com。

    猜你喜歡
    株數(shù)間伐油松
    梨園間伐改形技術(shù)探索與建議
    河北果樹(2021年4期)2021-12-02 01:14:56
    油松容器育苗技術(shù)
    不同間伐強(qiáng)度對(duì)香椿中齡林生長(zhǎng)和干形的影響
    森林撫育間伐技術(shù)的應(yīng)用
    ——以杉木撫育間伐為例
    綠色科技(2020年21期)2021-01-08 06:11:54
    間伐強(qiáng)度對(duì)杉木中齡林生長(zhǎng)和結(jié)構(gòu)的影響
    我院耐碳青霉烯類腸桿菌科細(xì)菌感染分布特點(diǎn)及耐藥性分析
    優(yōu)質(zhì)米水稻品種龍稻18配套栽培技術(shù)研究
    北方油松育苗技術(shù)
    巧解“植樹問(wèn)題”
    油松種植技術(shù)及主要病蟲害的防治對(duì)策
    丁香六月天网| 久久 成人 亚洲| 99九九在线精品视频| 亚洲熟女毛片儿| 精品亚洲成国产av| 岛国毛片在线播放| 国产精品偷伦视频观看了| 十分钟在线观看高清视频www| 国精品久久久久久国模美| av免费在线观看网站| 国产成人欧美| 19禁男女啪啪无遮挡网站| 久久久久国内视频| www.精华液| 欧美激情久久久久久爽电影 | 久久毛片免费看一区二区三区| 亚洲精品粉嫩美女一区| 一区在线观看完整版| 丰满人妻熟妇乱又伦精品不卡| 国产99久久九九免费精品| 九色亚洲精品在线播放| 一区福利在线观看| 亚洲国产欧美一区二区综合| 成人18禁在线播放| 九色亚洲精品在线播放| 国产精品免费一区二区三区在线 | 在线观看免费日韩欧美大片| 国产精品 欧美亚洲| 老司机亚洲免费影院| 人人妻,人人澡人人爽秒播| 91精品三级在线观看| 久久精品成人免费网站| 国产aⅴ精品一区二区三区波| 757午夜福利合集在线观看| 久久久久久人人人人人| 国产在线视频一区二区| 久久香蕉激情| 国精品久久久久久国模美| 国产成人欧美| 亚洲人成77777在线视频| 精品熟女少妇八av免费久了| 日日摸夜夜添夜夜添小说| 久久久精品区二区三区| 99国产精品免费福利视频| 黄片大片在线免费观看| 久久99热这里只频精品6学生| 亚洲成人免费电影在线观看| 免费在线观看完整版高清| 亚洲国产欧美一区二区综合| 亚洲avbb在线观看| 曰老女人黄片| 一进一出好大好爽视频| 久久精品亚洲av国产电影网| 国精品久久久久久国模美| 欧美 日韩 精品 国产| 另类亚洲欧美激情| 亚洲精品美女久久久久99蜜臀| 宅男免费午夜| 黄片播放在线免费| 国产成人av激情在线播放| 亚洲情色 制服丝袜| 午夜福利视频精品| 91精品国产国语对白视频| 1024视频免费在线观看| 久久精品成人免费网站| 日韩人妻精品一区2区三区| 这个男人来自地球电影免费观看| 黄片大片在线免费观看| 亚洲国产av新网站| 成年女人毛片免费观看观看9 | 亚洲黑人精品在线| 高清毛片免费观看视频网站 | 老司机深夜福利视频在线观看| 亚洲色图av天堂| 亚洲全国av大片| 乱人伦中国视频| 一个人免费在线观看的高清视频| 大码成人一级视频| 黄色视频在线播放观看不卡| 国产一卡二卡三卡精品| 国产黄频视频在线观看| 免费一级毛片在线播放高清视频 | 高清av免费在线| 曰老女人黄片| 99国产精品免费福利视频| 亚洲天堂av无毛| 日日爽夜夜爽网站| 亚洲国产毛片av蜜桃av| 欧美精品一区二区大全| 熟女少妇亚洲综合色aaa.| 美女国产高潮福利片在线看| 黄片小视频在线播放| 欧美日韩视频精品一区| 亚洲午夜理论影院| av线在线观看网站| 亚洲av第一区精品v没综合| av网站在线播放免费| 九色亚洲精品在线播放| 亚洲精品粉嫩美女一区| 91麻豆av在线| 中文亚洲av片在线观看爽 | 麻豆国产av国片精品| 啦啦啦中文免费视频观看日本| 飞空精品影院首页| 18禁美女被吸乳视频| 老司机福利观看| 免费在线观看视频国产中文字幕亚洲| 一区福利在线观看| 午夜精品久久久久久毛片777| 中文字幕高清在线视频| 国产男靠女视频免费网站| 国产伦人伦偷精品视频| 亚洲情色 制服丝袜| 国产高清视频在线播放一区| 性色av乱码一区二区三区2| 欧美成人免费av一区二区三区 | 每晚都被弄得嗷嗷叫到高潮| 国产av精品麻豆| 亚洲国产欧美在线一区| 亚洲三区欧美一区| 少妇精品久久久久久久| 久久av网站| 久热这里只有精品99| 一个人免费看片子| 欧美日韩亚洲综合一区二区三区_| 视频区欧美日本亚洲| 90打野战视频偷拍视频| 如日韩欧美国产精品一区二区三区| 一本色道久久久久久精品综合| 真人做人爱边吃奶动态| 婷婷成人精品国产| 亚洲成人国产一区在线观看| a在线观看视频网站| 欧美国产精品一级二级三级| 香蕉久久夜色| 精品一区二区三区视频在线观看免费 | 久久人妻av系列| 国精品久久久久久国模美| 在线 av 中文字幕| 纯流量卡能插随身wifi吗| 国产91精品成人一区二区三区 | 精品卡一卡二卡四卡免费| 亚洲va日本ⅴa欧美va伊人久久| 国产一区二区三区视频了| 99re在线观看精品视频| 91精品三级在线观看| 亚洲七黄色美女视频| 亚洲天堂av无毛| 正在播放国产对白刺激| 免费日韩欧美在线观看| 日韩视频在线欧美| 男人操女人黄网站| 国产精品久久久久久精品电影小说| 国产欧美日韩一区二区三| 如日韩欧美国产精品一区二区三区| 午夜日韩欧美国产| 久久毛片免费看一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美国产精品va在线观看不卡| 热99久久久久精品小说推荐| 十分钟在线观看高清视频www| 一二三四在线观看免费中文在| 久久久久精品人妻al黑| 91麻豆av在线| 18禁观看日本| 性高湖久久久久久久久免费观看| 亚洲伊人久久精品综合| 亚洲成人免费电影在线观看| 51午夜福利影视在线观看| 高清在线国产一区| 久久人人97超碰香蕉20202| 中文字幕另类日韩欧美亚洲嫩草| 久久这里只有精品19| 精品一区二区三区av网在线观看 | 国产精品久久电影中文字幕 | 99国产综合亚洲精品| 久久免费观看电影| 18在线观看网站| 国产欧美日韩精品亚洲av| 少妇裸体淫交视频免费看高清 | 黑人猛操日本美女一级片| 在线观看66精品国产| 欧美激情 高清一区二区三区| 久久精品亚洲av国产电影网| 怎么达到女性高潮| 天天躁夜夜躁狠狠躁躁| 亚洲九九香蕉| 中文字幕色久视频| 欧美 亚洲 国产 日韩一| 在线观看人妻少妇| 99精国产麻豆久久婷婷| 免费不卡黄色视频| 我要看黄色一级片免费的| 巨乳人妻的诱惑在线观看| 免费观看人在逋| av又黄又爽大尺度在线免费看| 91九色精品人成在线观看| 久久青草综合色| 色婷婷av一区二区三区视频| netflix在线观看网站| 亚洲伊人色综图| 天天添夜夜摸| 极品人妻少妇av视频| 国产成人精品久久二区二区免费| 美女视频免费永久观看网站| 成人亚洲精品一区在线观看| 18禁黄网站禁片午夜丰满| 黑人操中国人逼视频| 狠狠精品人妻久久久久久综合| 99re6热这里在线精品视频| 中国美女看黄片| 啦啦啦中文免费视频观看日本| 久久久久久久国产电影| 国产欧美亚洲国产| 色婷婷av一区二区三区视频| 不卡av一区二区三区| 免费少妇av软件| 老鸭窝网址在线观看| 欧美日韩亚洲国产一区二区在线观看 | www.999成人在线观看| 中文字幕另类日韩欧美亚洲嫩草| 无人区码免费观看不卡 | 国产欧美日韩精品亚洲av| 美女主播在线视频| 9191精品国产免费久久| 国产一区有黄有色的免费视频| 国产福利在线免费观看视频| 多毛熟女@视频| 国产精品久久久久久精品电影小说| 18禁观看日本| 国产91精品成人一区二区三区 | 日韩免费高清中文字幕av| 亚洲专区中文字幕在线| 日韩欧美一区二区三区在线观看 | 国产欧美亚洲国产| 如日韩欧美国产精品一区二区三区| 男女高潮啪啪啪动态图| 亚洲欧美一区二区三区久久| 精品亚洲成a人片在线观看| videosex国产| 视频在线观看一区二区三区| 人妻久久中文字幕网| 老司机深夜福利视频在线观看| 国产日韩欧美亚洲二区| 大陆偷拍与自拍| 日本vs欧美在线观看视频| 激情在线观看视频在线高清 | 亚洲伊人久久精品综合| 在线亚洲精品国产二区图片欧美| 久久精品亚洲av国产电影网| 少妇 在线观看| 国产一区二区三区在线臀色熟女 | svipshipincom国产片| 日韩制服丝袜自拍偷拍| 国产色视频综合| 在线 av 中文字幕| 狠狠婷婷综合久久久久久88av| 精品第一国产精品| 国产精品久久久人人做人人爽| aaaaa片日本免费| 建设人人有责人人尽责人人享有的| 黄频高清免费视频| 午夜福利视频精品| 日韩欧美三级三区| 正在播放国产对白刺激| 欧美黑人精品巨大| av又黄又爽大尺度在线免费看| 久久久久久久国产电影| 国产熟女午夜一区二区三区| 亚洲熟女毛片儿| 久久99一区二区三区| 中文字幕高清在线视频| 欧美激情极品国产一区二区三区| 色视频在线一区二区三区| 巨乳人妻的诱惑在线观看| 热99re8久久精品国产| 黑人巨大精品欧美一区二区mp4| 51午夜福利影视在线观看| 一级毛片高清免费大全| 麻豆国产97在线/欧美| 一本精品99久久精品77| 色综合亚洲欧美另类图片| www国产在线视频色| 又粗又爽又猛毛片免费看| 亚洲精华国产精华精| 亚洲自偷自拍图片 自拍| 波多野结衣高清作品| 欧美日韩亚洲国产一区二区在线观看| 97人妻精品一区二区三区麻豆| 亚洲国产日韩欧美精品在线观看 | 嫁个100分男人电影在线观看| 亚洲欧美日韩东京热| 嫩草影院精品99| 啦啦啦韩国在线观看视频| 久久婷婷人人爽人人干人人爱| 又爽又黄无遮挡网站| 黄色 视频免费看| 精品国产乱码久久久久久男人| 午夜成年电影在线免费观看| 日日摸夜夜添夜夜添小说| 成年人黄色毛片网站| av黄色大香蕉| 久久精品亚洲精品国产色婷小说| 国产欧美日韩一区二区三| 午夜精品久久久久久毛片777| 亚洲午夜理论影院| 亚洲一区高清亚洲精品| 久久久精品大字幕| 国产午夜福利久久久久久| 日本精品一区二区三区蜜桃| 日韩高清综合在线| 午夜精品一区二区三区免费看| 国产亚洲精品久久久久久毛片| 欧美性猛交黑人性爽| 日本撒尿小便嘘嘘汇集6| 国产精品一区二区精品视频观看| 亚洲欧美日韩无卡精品| 日本免费a在线| 国产aⅴ精品一区二区三区波| 国产亚洲av高清不卡| АⅤ资源中文在线天堂| 99久久精品国产亚洲精品| 免费在线观看视频国产中文字幕亚洲| 亚洲av电影在线进入| 亚洲在线观看片| 亚洲片人在线观看| 色av中文字幕| 精品国产三级普通话版| 法律面前人人平等表现在哪些方面| av天堂在线播放| 最近在线观看免费完整版| 亚洲 欧美 日韩 在线 免费| www.www免费av| 婷婷亚洲欧美| cao死你这个sao货| 国产精品一区二区精品视频观看| 午夜福利在线观看吧| 亚洲av五月六月丁香网| 日本 欧美在线| 老司机在亚洲福利影院| 99国产精品一区二区蜜桃av| h日本视频在线播放| 男女之事视频高清在线观看| 欧美精品啪啪一区二区三区| 综合色av麻豆| 亚洲 欧美一区二区三区| 宅男免费午夜| 九色国产91popny在线| 日韩欧美一区二区三区在线观看| 中文在线观看免费www的网站| 久久这里只有精品中国| 欧美极品一区二区三区四区| 搡老岳熟女国产| 久久精品国产99精品国产亚洲性色| 欧美又色又爽又黄视频| 欧美丝袜亚洲另类 | 国产成人影院久久av| 网址你懂的国产日韩在线| 性色av乱码一区二区三区2| 国产精品亚洲美女久久久| 偷拍熟女少妇极品色| 啦啦啦免费观看视频1| 欧美丝袜亚洲另类 | 亚洲人与动物交配视频| 99国产极品粉嫩在线观看| 别揉我奶头~嗯~啊~动态视频| 亚洲成人久久性| 香蕉国产在线看| 在线观看66精品国产| 18美女黄网站色大片免费观看| 黄色成人免费大全| 俺也久久电影网| 亚洲国产高清在线一区二区三| 男人的好看免费观看在线视频| 在线观看免费视频日本深夜| 久久性视频一级片| 国产精品综合久久久久久久免费| 给我免费播放毛片高清在线观看| 欧美日韩黄片免| 毛片女人毛片| 精品电影一区二区在线| 在线观看日韩欧美| 欧美+亚洲+日韩+国产| 99热6这里只有精品| 女同久久另类99精品国产91| 午夜福利视频1000在线观看| 国产一区二区三区视频了| 国产欧美日韩一区二区精品| 久久精品人妻少妇| 成人鲁丝片一二三区免费| 午夜影院日韩av| 精品久久蜜臀av无| 精品国产三级普通话版| xxx96com| 女生性感内裤真人,穿戴方法视频| 99国产精品99久久久久| 极品教师在线免费播放| 久久久久国产精品人妻aⅴ院| 免费搜索国产男女视频| 一区二区三区激情视频| 一本一本综合久久| 欧美+亚洲+日韩+国产| 在线看三级毛片| 日韩成人在线观看一区二区三区| 国产激情久久老熟女| 久久这里只有精品中国| 99热这里只有是精品50| 色综合站精品国产| 亚洲自偷自拍图片 自拍| 欧美zozozo另类| 曰老女人黄片| 免费观看的影片在线观看| 国产极品精品免费视频能看的| 好男人在线观看高清免费视频| 午夜福利高清视频| aaaaa片日本免费| 成人三级黄色视频| 亚洲第一电影网av| 长腿黑丝高跟| 国产亚洲精品av在线| 偷拍熟女少妇极品色| 757午夜福利合集在线观看| 日本 欧美在线| 国产精品久久久久久久电影 | 丁香六月欧美| 日韩欧美免费精品| 婷婷亚洲欧美| 两性夫妻黄色片| 午夜福利在线观看免费完整高清在 | 亚洲,欧美精品.| 精品一区二区三区av网在线观看| 90打野战视频偷拍视频| 国内精品一区二区在线观看| 俺也久久电影网| 国产欧美日韩精品亚洲av| 国产精品野战在线观看| 黄频高清免费视频| 日本成人三级电影网站| 成人性生交大片免费视频hd| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久精品国产欧美久久久| 天天躁狠狠躁夜夜躁狠狠躁| 久久午夜亚洲精品久久| 波多野结衣高清无吗| 亚洲第一欧美日韩一区二区三区| 日韩欧美一区二区三区在线观看| 99热精品在线国产| 91av网站免费观看| 又紧又爽又黄一区二区| 毛片女人毛片| 午夜免费成人在线视频| 欧美日韩福利视频一区二区| 国产精品99久久99久久久不卡| 90打野战视频偷拍视频| 亚洲人成电影免费在线| 欧洲精品卡2卡3卡4卡5卡区| 人人妻人人澡欧美一区二区| 日日干狠狠操夜夜爽| 少妇的逼水好多| 五月玫瑰六月丁香| 老司机福利观看| 国产精品久久久久久亚洲av鲁大| 久久欧美精品欧美久久欧美| 嫩草影院入口| 夜夜夜夜夜久久久久| 亚洲精品在线美女| 欧美乱妇无乱码| 亚洲专区国产一区二区| 男女视频在线观看网站免费| 两性午夜刺激爽爽歪歪视频在线观看| 90打野战视频偷拍视频| 高潮久久久久久久久久久不卡| 精品乱码久久久久久99久播| 国产高清三级在线| 国产精品一及| 最近最新中文字幕大全电影3| 99久久精品热视频| 色尼玛亚洲综合影院| 国产精品香港三级国产av潘金莲| 国产精品亚洲av一区麻豆| 每晚都被弄得嗷嗷叫到高潮| 久久九九热精品免费| 九九在线视频观看精品| 啪啪无遮挡十八禁网站| 免费看光身美女| av视频在线观看入口| 不卡av一区二区三区| 国产精品女同一区二区软件 | 亚洲欧美日韩东京热| 少妇裸体淫交视频免费看高清| avwww免费| 国产亚洲精品av在线| 久久久久久国产a免费观看| 国产激情欧美一区二区| 婷婷六月久久综合丁香| 美女扒开内裤让男人捅视频| 色噜噜av男人的天堂激情| 老司机福利观看| 美女午夜性视频免费| 老司机福利观看| 国产一区二区在线观看日韩 | 午夜福利高清视频| 俄罗斯特黄特色一大片| 亚洲无线观看免费| 琪琪午夜伦伦电影理论片6080| 看片在线看免费视频| 亚洲精品美女久久av网站| 窝窝影院91人妻| 九九热线精品视视频播放| 人人妻人人看人人澡| 亚洲av电影在线进入| 中文资源天堂在线| 极品教师在线免费播放| 无限看片的www在线观看| 日本免费a在线| 成人永久免费在线观看视频| 极品教师在线免费播放| 国产熟女xx| а√天堂www在线а√下载| 99久久99久久久精品蜜桃| 天天添夜夜摸| 韩国av一区二区三区四区| 狠狠狠狠99中文字幕| 日韩国内少妇激情av| 欧美在线一区亚洲| 国产黄片美女视频| 俺也久久电影网| 两性午夜刺激爽爽歪歪视频在线观看| 我要搜黄色片| 在线十欧美十亚洲十日本专区| 搡老妇女老女人老熟妇| 少妇熟女aⅴ在线视频| 亚洲七黄色美女视频| 曰老女人黄片| 又大又爽又粗| 女同久久另类99精品国产91| 十八禁人妻一区二区| 国产精品亚洲美女久久久| 国产成人精品久久二区二区91| 国产精品98久久久久久宅男小说| 两个人视频免费观看高清| 欧美日本视频| 亚洲av日韩精品久久久久久密| 国产av麻豆久久久久久久| 听说在线观看完整版免费高清| 久久久久免费精品人妻一区二区| 叶爱在线成人免费视频播放| 欧美精品啪啪一区二区三区| 免费人成视频x8x8入口观看| 黄色成人免费大全| 嫁个100分男人电影在线观看| 亚洲精品在线美女| 国产成人一区二区三区免费视频网站| 午夜免费成人在线视频| 美女高潮的动态| 丁香欧美五月| svipshipincom国产片| 美女高潮喷水抽搐中文字幕| 亚洲 欧美一区二区三区| 久久久久久大精品| 亚洲成av人片在线播放无| 18禁美女被吸乳视频| av福利片在线观看| 午夜免费激情av| av在线蜜桃| 免费在线观看成人毛片| 亚洲国产精品sss在线观看| 欧美另类亚洲清纯唯美| 宅男免费午夜| 脱女人内裤的视频| 天天添夜夜摸| 欧美激情久久久久久爽电影| 麻豆av在线久日| 一进一出抽搐gif免费好疼| 桃色一区二区三区在线观看| 丁香欧美五月| 级片在线观看| 哪里可以看免费的av片| 欧美高清成人免费视频www| 日本免费a在线| 黄色成人免费大全| av女优亚洲男人天堂 | 日本一二三区视频观看| 午夜福利在线观看免费完整高清在 | 精品熟女少妇八av免费久了| 午夜福利在线在线| 成人三级做爰电影| 亚洲,欧美精品.| 99久久99久久久精品蜜桃| 成人午夜高清在线视频| 午夜两性在线视频| 国产精品美女特级片免费视频播放器 | 色综合站精品国产| 亚洲专区字幕在线| a级毛片在线看网站| 两个人的视频大全免费| 美女高潮的动态| 色综合亚洲欧美另类图片| 欧美最黄视频在线播放免费| 极品教师在线免费播放| 婷婷精品国产亚洲av| 1024手机看黄色片| 亚洲自拍偷在线| 久久午夜综合久久蜜桃| 日韩三级视频一区二区三区| 欧美乱色亚洲激情| 高清毛片免费观看视频网站| 国产三级中文精品| 亚洲精品456在线播放app | 亚洲成av人片免费观看| 此物有八面人人有两片| 国产精品亚洲av一区麻豆| 国产蜜桃级精品一区二区三区| 真实男女啪啪啪动态图| 嫩草影院精品99|