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

    基于第一著枝角度的羅田垂枝杉樹冠表面積預(yù)估模型研究

    2024-01-10 03:49:46邊更戰(zhàn)陳書杭
    湖南林業(yè)科技 2023年6期
    關(guān)鍵詞:羅田冠幅樹冠

    莫 嬃,易 烜,邊更戰(zhàn),陳書杭

    (1.中南林業(yè)科技大學(xué),湖南 長沙 410004;2.湖南省青羊湖國有林場,湖南 長沙 410600;3.中部林業(yè)產(chǎn)權(quán)交易服務(wù)中心,湖南 長沙 410007;4.咸寧市林業(yè)科學(xué)院,湖北 咸寧 437000)

    樹冠是樹木的重要組成部分,是樹木生長與生理過程的主要部位[1]。樹冠的大小能夠體現(xiàn)樹木的生長潛力以及競爭能力,也是預(yù)估樹木生長量的基本依據(jù)之一[2]。樹冠三維空間的大小通常用樹冠表面積和樹冠體積表示[3-7]。樹冠表面積、體積等復(fù)合指標(biāo)不僅與林木生長和健康有著密切的關(guān)系,而且能用其構(gòu)建生長空間競爭指數(shù)和生長空間指數(shù)[8-10],因此,研究樹冠表面積對進一步分析和優(yōu)化調(diào)整林分的空間結(jié)構(gòu)具有重要意義。目前,對于樹冠表面積量化方面的研究仍然存在許多挑戰(zhàn)[11]。傳統(tǒng)的樹冠表面積測量方法通常需要大量的人工測量和復(fù)雜的數(shù)據(jù)處理,需要消耗大量的人力、物力和財力[12]。幾何體法是常用的樹冠表面積計算方法,即將樹冠近似看作圓錐體進行計算[13]。早在1969年,Hamilton[14]對西特喀云杉(Stikaspruce)進行了研究,在計算時將樹冠看成圓錐體很容易得出樹冠的表面積。近年來,涂宏濤等[15]、吳明欽等[16]在構(gòu)建樹冠表面積模型時,為使得到的數(shù)據(jù)更加精確,常常將幾何體法與平均斷面求積法結(jié)合起來,也稱其為“分層切割法”。由于樹冠結(jié)構(gòu)的復(fù)雜性,很難對樹冠表面積進行準(zhǔn)確的測算[17],因此,在需要使用新方法測算樹冠表面積時,常常用“分層切割法”計算的值作為對照,以評判新方法的準(zhǔn)確性和可信度[18]。值得一提的是,不論是幾何體法還是分層切割法,其本質(zhì)都是將樹冠看作與其相近的幾何體,再進行拆解、估算。而在這兩種樹冠表面積測算法中,對于樹冠下層側(cè)枝的著枝狀態(tài)都未曾引起重視。我們以羅田垂枝杉(Cunninghamialanceolatavar.luotian)為研究對象[19-21],建立適用于該樹種的樹冠表面積預(yù)估模型,為其經(jīng)營管理提供參考。

    1 研究區(qū)概況

    羅田垂枝杉試驗林于2016年初造林,其造林密度為1.5m×1.5m,位于湖北省咸寧市東南部的賀勝橋鎮(zhèn)。該鎮(zhèn)地處114°16′—114°29′E,29°57′ — 30°02′N,其東西長31.5km,南北寬10.1km,最高山脈海拔264m,東部為丘陵崗地,西部為平原湖區(qū),地形東高西低,由東南向西北呈現(xiàn)為高丘、丘崗、平湖。研究區(qū)氣候?qū)俚湫偷膩啛釒Ъ撅L(fēng)氣候,冬冷夏涼,四季分明,光照充足,雨量充沛。冬季盛行偏北風(fēng),偏冷干燥;夏季盛行偏南風(fēng),高溫多雨。年平均氣溫16.8℃,極端最高氣溫41.4℃,極端最低氣溫-15.4℃。年平均降水量1577.4mm,年平均日照時數(shù)1754.5h,年平均無霜期245~258d[22]。

    2 研究方法

    2.1 數(shù)據(jù)采集

    2023年4—5月,在湖北省咸寧市賀勝橋鎮(zhèn)的垂枝杉試驗林隨機設(shè)置30塊15m×20m的標(biāo)準(zhǔn)地,并對標(biāo)準(zhǔn)地內(nèi)的林木進行每木檢尺,調(diào)查記錄胸徑(D)、樹高(H)、活枝下高(HCB)、冠幅(CW)、冠長(CL)等測樹因子,然后根據(jù)林分平均胸徑在每塊標(biāo)準(zhǔn)地內(nèi)選取2~3株平均木作為標(biāo)準(zhǔn)木,共選取71株標(biāo)準(zhǔn)木。將標(biāo)準(zhǔn)木伐倒后,采用分層切割法將樹冠按1 m的長度分層,每層內(nèi)選取3~5個標(biāo)準(zhǔn)枝作為樣枝測量枝長及其著枝角度。標(biāo)準(zhǔn)木主要測樹因子調(diào)查數(shù)據(jù)統(tǒng)計結(jié)果見表1。

    表1 羅田垂枝杉標(biāo)準(zhǔn)木主要測樹因子調(diào)查數(shù)據(jù)統(tǒng)計結(jié)果Tab.1Survey data statistics of main measuring tree factors of Cunninghamia lanceolata var.luotian model trees測樹因子模型擬合數(shù)據(jù)n=54模型檢驗數(shù)據(jù)n=17平均值最小值最大值標(biāo)準(zhǔn)差平均值最小值最大值標(biāo)準(zhǔn)差胸徑(D)/cm10.116.4013.201.739.745.3012.601.87樹高(H)/m7.816.609.000.527.545.608.800.73活枝下高(HCB)/m4.803.206.100.674.763.305.700.73冠長(CW)/m3.011.704.300.552.772.003.700.51冠幅(CL)/m0.780.401.170.160.850.451.100.16著枝角(θ)/°121.14100.30146.5811.85114.97101.09143.5915.13

    2.2 樹冠表面積計算

    采用分層切割法計算樹冠表面積[14],即將羅田垂枝杉的樹冠按1 m長度分成若干層,將其樹冠底部至中部每層視為圓臺體,將最頂層視為圓錐體;分別計算每個冠層的平均枝長,再根據(jù)枝長(L)與著枝角度(B)之間的三角函數(shù)關(guān)系計算每個冠層的頂部和底部半徑(C)。

    C=BLsinθ

    (1)

    用圓臺體、圓錐體表面積計算公式計算樹冠各層表面積,累加樹冠各層表面積得到全樹冠表面積。

    (2)

    (3)

    式中:h為冠層長度;CR為冠層底部半徑;Cr為冠層頂部半徑。

    2.3 樹冠表面積預(yù)估模型建立

    2.3.1 模型結(jié)構(gòu)設(shè)計

    根據(jù)實測所得的71株標(biāo)準(zhǔn)木調(diào)查數(shù)據(jù),建模數(shù)據(jù)與模型檢驗數(shù)據(jù)按約3∶1的比例隨機抽取,抽取其中54株標(biāo)準(zhǔn)木的解析數(shù)據(jù)作為建模數(shù)據(jù),剩余17株標(biāo)準(zhǔn)木的解析數(shù)據(jù)作為模型檢驗數(shù)據(jù)。相關(guān)研究表明,異速生長方程(式4)能夠很好地擬合樹冠表面積模型[23]。并運用R 4.3.1軟件進行回歸分析,建立預(yù)估方程。

    (4)

    式中:y為樹冠表面積;x1、x2、xi均為林木測樹因子變量;a、b、c、d均為待估參數(shù)。

    2.3.2 模型變量篩選

    根據(jù)實測數(shù)據(jù),將測樹因子變量與樹冠表面積進行 Pearson 相關(guān)性分析[24],結(jié)果見表2。由表2可知:胸徑、樹高、冠長、冠幅與樹冠表面積的相關(guān)性均達(dá)顯著或極顯著,其中樹冠表面積與冠幅的相關(guān)系數(shù)最大;活枝下高與樹冠表面積的相關(guān)系數(shù)最小。

    表2 樹冠表面積與標(biāo)準(zhǔn)木測樹因子變量的Pearson相關(guān)性分析結(jié)果Tab.2 Pearson relation analysis between main measuring tree dendrometry factor variables and crown surface area變量胸徑樹高冠長冠幅活枝下高樹冠表面積0.779**0.377**0.756**0.839**-2.81* 注:*表示相關(guān)顯著(P<0.05); **表示相關(guān)極顯著(P<0.01)。

    2.3.3 模型評價與檢驗

    對建模樣本進行檢驗,計算檢驗指標(biāo),評價模型優(yōu)度。使用確定系數(shù)(R2)、精度(EA)、剩余標(biāo)準(zhǔn)差(S)、平均相對誤差(E)和平均相對誤差絕對值(e)[15]對備選模型進行檢驗,選取最優(yōu)模型。其中,R2、EA越接近1越好;S、E、e越接近0越好。

    (5)

    (6)

    (7)

    (8)

    (9)

    3 結(jié)果與分析

    3.1 基于第一著枝角度的單木樹冠表面積計算

    3.1.1 傳統(tǒng)幾何體法

    幾何體法是計算樹冠表面積的傳統(tǒng)方法,即將樹冠近似看作圓錐體進行估算,其計算公式為式(10)。通常我們在使用幾何體法計算樹冠表面積時,對于冠長(圓錐體的高)的定義為樹冠頂部到第一活枝著枝點的垂直長度,即樹高與枝下高之差(CL=H-HCB)。其樹冠結(jié)構(gòu)示意圖見圖1。

    圖1 幾何體法樹冠結(jié)構(gòu)示意圖Fig.1 Schematicdiagram of tree crown structure using geometric method

    當(dāng)?=90°時,樹冠表面積傳統(tǒng)幾何體法的計算公式為:

    (10)

    (11)

    式中:CSAi為第i株木的樹冠表面積;CWi為第i株木的冠幅;CLi為第i株木的冠長;CWewi、CWsni分別為第i株木的東西、南北冠幅。冠長(CLi)為樹高與枝下高之差。

    3.1.2 基于第一著枝角度的幾何體法

    羅田垂枝杉與普通杉木相比有著十分顯著的形態(tài)差異,其側(cè)枝自然下垂,樹冠下層側(cè)枝著生角度達(dá)150°左右。因此,其冠長若仍以樹高與枝下高之差值來定義,在使用傳統(tǒng)幾何體法計算時得到的表面積將與實際表面積存在誤差,其誤差如圖2所示。

    圖2 幾何體法計算得到的羅田垂枝杉樹冠表面積與其實際表面積的誤差示意圖Fig.2 Schematic diagram of the error betweenthe surface area of Cunninghamia lanceolata var.luotian tree crown calculated by the geometric method and its actual surface area

    由圖2可知,在第一著枝角度?的影響下,垂枝杉冠底與著枝點在垂直方向上存在差值Δcl,與著枝角度?的關(guān)系為:

    (12)

    因此,垂枝杉實際冠長為傳統(tǒng)冠長與ΔCL之和,即CL+ΔCL。由此可推算出適用于在?>90°時羅田垂枝杉樹冠表面積的幾何體法計算公式為:

    (13)

    式中:CSAj為第j株木的樹冠表面積;CWj為第j株木的冠幅;CLj為第j株木的冠長;冠長(CLj)為樹高與枝下高之差。

    3.2 第一著枝角度對樹冠表面積預(yù)估精度的影響

    3.2.1 樹冠表面積不同計算方法的精度比較

    利用實測得到的71株垂枝杉標(biāo)準(zhǔn)木調(diào)查數(shù)據(jù),分別使用傳統(tǒng)幾何體法、基于第一著枝角度的幾何體法以及分層切割法計算垂枝杉樹冠的表面積,結(jié)果見表3。由表3可知,樹冠表面積的平均值、最大值、最小值和標(biāo)準(zhǔn)差都是分層切割法的>基于第一著枝角度的幾何體法的>傳統(tǒng)幾何體法的。

    表3 3種計算方法得到的羅田垂枝杉樹冠表面積Tab.3Tree crown surface area of Cunninghamia lanceo-lata var.luotian calculated by 3 methodsm2計算方法平均值最大值最小值標(biāo)準(zhǔn)差傳統(tǒng)幾何體法3.787.981.391.21基于第一著枝角度的幾何體法4.048.161.761.20分層切割法4.729.582.131.46

    由于垂枝杉樹冠形態(tài)的不規(guī)則性以及實際操作的困難性,目前無法獲取其樹冠的實際表面積。為比較傳統(tǒng)幾何體法與基于第一著枝角度的幾何體法的計算精度,以分層切割法得出的表面積作為參照,將分層切割法、傳統(tǒng)幾何體法與基于第一著枝角度的幾何體法計算得出的樹冠表面積繪制成圖3。圖3結(jié)果表明,分層切割法得出的表面積最大,基于第一著枝角度的幾何體法得出的次之,傳統(tǒng)幾何體法得出的最小??梢姡诘谝恢嵌鹊膸缀误w法得出的樹冠表面積較傳統(tǒng)幾何體法得出的更加接近分層切割法得出的,說明對于垂枝杉而言,基于第一著枝角度的幾何體法得出的樹冠表面積較傳統(tǒng)幾何體法得出的更加準(zhǔn)確。

    圖3 3種方法計算得出的樹冠表面積(小提琴圖與箱線圖)Fig.3 Crown surface area calculated by three methodes(violin chart and box chart)

    3.2.2 樹冠表面積與第一著枝角度的相關(guān)性

    樹冠第一著枝角度是樹冠表面積大小的決定性因素之一。將第一著枝角度大于90°的部分,即第一著枝角與垂直角的差值,和與其對應(yīng)帶來的表面積差值繪制成散點圖(圖4)。由圖4可知,第一著枝角度差與樹冠表面積差沒有明顯的線性關(guān)系。在樹冠表面積的相關(guān)測樹因子一致的前提下,繼續(xù)繪制角度差與表面積差的散點圖。

    圖4 角度差與表面積差的散點圖Fig.4 Scatter plot of angle difference versus surface area difference

    先考慮冠幅一致時第一著枝角度與樹冠表面積的關(guān)系。選取冠幅為0.7 m的9株標(biāo)準(zhǔn)木,繪制第一著枝角度差與樹冠表面積差的散點圖(圖5a)。由圖5a可知,在冠幅相同的情況下,角度差越大表面積差也越大,即第一著枝角度的大小與樹冠表面積的大小呈正相關(guān)。

    再考慮冠長相同時第一著枝角度與樹冠表面積的關(guān)系。選取冠長為3.1 m的9株標(biāo)準(zhǔn)木,繪制第一著枝角度差與樹冠表面積差的散點圖(圖5b)。由圖5b可知,在冠長相同的情況下,角度差越大表面積差也越大,說明第一著枝角度大小與樹冠表面積的大小也呈正相關(guān)。

    3.3 羅田垂枝杉樹冠表面積預(yù)估模型構(gòu)建

    以建模數(shù)據(jù)擬合樹冠表面積預(yù)估模型。先引入單變量胸徑(D)擬合為一元變量模型,再分別引入冠幅(CW)、冠長(CL)擬合為二元變量模型,然后用3個變量擬合成三元變量模型,最后分別選取3種擬合模型中R2(擬合優(yōu)度或判斷系數(shù))、EA精度(精度)最大的模型,并計算參數(shù)估計值、評價指標(biāo)值。擬合結(jié)果見表4。由表4可知:隨著模型變量的增加,模型方程的R2(擬合優(yōu)度或判斷系數(shù))不斷增大,EA也一直呈增加趨勢,說明模型的精度不斷提高。其中一元變量模型的擬合優(yōu)度和精度均最低,分別為0.62和81.99%;三元變量模型的均最大,分別為0.98、96.89%;在二元變量模型中,模型3的R2和EA均較大??紤]到實際調(diào)查中數(shù)據(jù)獲取的難易程度,選取模型3和模型4為備選模型。

    表4 羅田垂枝杉樹冠表面積模型擬合結(jié)果Tab.4 Fitting results of crown surface area model for Cunninghamia lanceolata var.luotian模型編號變量參數(shù)估計評價指標(biāo)abcdR2EA/%1D0.191.410.6281.992D、CW1.390.630.870.7786.093D、CL1.010.810.830.8187.144D、CW、CL1.970.070.940.870.9896.89

    模型3、4的表達(dá)式如下。

    模型3:CSA=1.01D0.81CL0.83

    (14)

    模型4:CSA=1.97D0.07CW0.94CL0.87

    (15)

    以模型檢驗數(shù)據(jù)對模型3、4進行檢驗,檢驗結(jié)果如表5所示。表5結(jié)果顯示:在2個樹冠表面積備選模型中,以模型4為優(yōu),其預(yù)估精度EA達(dá)96.95%,平均相對誤差E僅為0.01;模型3的預(yù)估精度EA也達(dá)到了80%以上,且平均相對誤差E為0.56,因此,2個備選模型都能較好的預(yù)估羅田垂枝杉的樹冠表面積。

    表5 羅田垂枝杉樹冠表面積模型檢驗結(jié)果Tab.5Test results of crown surface area model for Cun-ninghamia lanceolata var.luotian模型編號方程變量評價指標(biāo)R2EA/%SeE3D,CL0.7484.622.200.680.564D,CW,CL0.9896.950.620.130.01

    4 結(jié)論與討論

    (1)3種計算方法得到的羅田垂枝杉樹冠表面積的大小為分層切割法的>基于第一著枝角度的幾何體法的>傳統(tǒng)幾何體法的,表明以分層切割法作為參照時,基于第一著枝角度的幾何體法計算樹冠表面積的精度高于傳統(tǒng)幾何體法的,也進一步驗證了第一著枝角度是樹冠表面積的決定性因素之一。

    (2)羅田垂枝杉樹冠表面積與胸徑、冠長、冠幅等均呈顯著正相關(guān),其中與冠幅的相關(guān)系數(shù)最高;異速生長方程能很好地擬合樹冠表面積模型,測樹因子變量的增加能適當(dāng)?shù)靥岣吣P蛿M合的精度。其中以胸徑、冠長以及冠幅為變量的三元變量模型來擬合樹冠表面積的模型為最優(yōu)模型,模型的判斷系數(shù)R2達(dá)0.98,模型的精度EA也達(dá)到了96.95%,說明所建立的模型具有較高的估測精度,其模型為CSA=1.97D0.07CW0.94CL0.87。

    (3)在針對樹冠體積、表面積估測模型的研究中[25],林木的胸徑、冠幅等測樹因子都還需要通過實測手段來獲取數(shù)據(jù),并尋求與樹冠形態(tài)相近的表面積、體積公式進行計算,這種計算方式無疑過于粗放。隨著科學(xué)技術(shù)的發(fā)展,林業(yè)資源調(diào)查已不僅僅局限于傳統(tǒng)方式。如利用地面三維激光掃描技術(shù)對建模樣木進行不同方向的掃描來擬合獲得整株林木的完整點云數(shù)據(jù),再通過數(shù)據(jù)處理生成樹冠的三維點云模型[26];還有通過無人機傾斜攝影測量,對樹木的模型進行三維重建,從而得到樹冠表面積與體積[27]。顯然,因為新科技產(chǎn)生了新方法,使計算出的樹冠表面積、體積都更加精確,且大大節(jié)省了人力、物力。

    (4)本文在計算樹冠表面積時考慮了樹冠第一著枝角度,其表面積的計算更具靈活性。羅田垂枝杉的樹冠第一著枝角度大于90°,而在自然界中不乏枝條與地面反向生長的樹種,即樹冠第一著枝角度小于90°,對于此類型樹種樹冠表面積與第一著枝角度的關(guān)系,也有待進一步研究。

    猜你喜歡
    羅田冠幅樹冠
    無人機遙感影像提取的單木冠幅數(shù)據(jù)在桉樹林分蓄積量估測中的應(yīng)用1)
    徐建南
    城市綠地微環(huán)境對土壤動物群落多樣性的影響
    施肥對三江平原丘陵區(qū)長白落葉松人工林中齡林單木樹冠圓滿度影響
    樹冠羞避是什么原理?
    榕樹
    樹冠
    文學(xué)港(2019年5期)2019-05-24 14:19:42
    基于無人機高分影像的冠幅提取與樹高反演
    雪中浪漫 羅田薄刀峰
    幸福(2016年6期)2016-12-01 03:08:27
    羅田金秋
    黄色视频在线播放观看不卡| 国语对白做爰xxxⅹ性视频网站| 在线观看三级黄色| 欧美日韩在线观看h| 日韩人妻高清精品专区| 亚洲欧美一区二区三区黑人 | 欧美精品一区二区免费开放| 伦理电影免费视频| 日本黄色片子视频| 亚洲无线观看免费| 午夜av观看不卡| 国产欧美另类精品又又久久亚洲欧美| 中文精品一卡2卡3卡4更新| 深夜a级毛片| 久久久午夜欧美精品| 男女免费视频国产| 亚洲经典国产精华液单| av福利片在线观看| 在线看a的网站| 久久人妻熟女aⅴ| tube8黄色片| 一区二区三区乱码不卡18| 九色成人免费人妻av| av免费观看日本| 国产男女超爽视频在线观看| 尾随美女入室| 久热这里只有精品99| 国产日韩欧美亚洲二区| 国产精品.久久久| 国产精品福利在线免费观看| 校园人妻丝袜中文字幕| 欧美丝袜亚洲另类| 大香蕉久久网| 亚洲国产精品专区欧美| 美女cb高潮喷水在线观看| 久久久久久久久久久丰满| av有码第一页| www.色视频.com| 午夜福利视频精品| 国产在线一区二区三区精| 大香蕉久久网| 中国三级夫妇交换| 亚洲国产精品成人久久小说| 熟女av电影| 日日啪夜夜撸| 日韩成人av中文字幕在线观看| 麻豆精品久久久久久蜜桃| 欧美变态另类bdsm刘玥| 国产精品福利在线免费观看| 国产欧美另类精品又又久久亚洲欧美| 久久精品久久久久久噜噜老黄| 精品国产露脸久久av麻豆| 国产亚洲91精品色在线| 久久久久久久久久人人人人人人| 成年女人在线观看亚洲视频| 亚洲av成人精品一区久久| 久久久久久人妻| 久久精品国产鲁丝片午夜精品| 欧美97在线视频| 久久亚洲国产成人精品v| 岛国毛片在线播放| 一级毛片黄色毛片免费观看视频| 22中文网久久字幕| 26uuu在线亚洲综合色| 桃花免费在线播放| 中文精品一卡2卡3卡4更新| 成人18禁高潮啪啪吃奶动态图 | 免费观看av网站的网址| 亚洲精品乱码久久久久久按摩| 日本黄大片高清| 日韩欧美一区视频在线观看 | 精品一品国产午夜福利视频| 成年女人在线观看亚洲视频| 日韩成人伦理影院| 成人免费观看视频高清| 最近中文字幕2019免费版| 制服丝袜香蕉在线| 亚洲天堂av无毛| 一本—道久久a久久精品蜜桃钙片| 在线观看国产h片| 国产免费视频播放在线视频| 曰老女人黄片| 中文字幕亚洲精品专区| 亚洲精品一区蜜桃| 水蜜桃什么品种好| 亚洲高清免费不卡视频| 一本大道久久a久久精品| 久久韩国三级中文字幕| av国产精品久久久久影院| 人人妻人人添人人爽欧美一区卜| 国精品久久久久久国模美| 人人澡人人妻人| 国产成人a∨麻豆精品| 欧美性感艳星| 亚洲,一卡二卡三卡| 最近手机中文字幕大全| 韩国高清视频一区二区三区| 99热这里只有精品一区| 美女中出高潮动态图| 五月开心婷婷网| 精品国产一区二区久久| 精品一区二区三区视频在线| 久久精品国产亚洲av涩爱| 欧美精品一区二区免费开放| 美女福利国产在线| 精品久久久噜噜| 日韩欧美 国产精品| 国国产精品蜜臀av免费| 久热久热在线精品观看| 久久人妻熟女aⅴ| 有码 亚洲区| 少妇人妻精品综合一区二区| 五月开心婷婷网| 国产成人免费观看mmmm| 97精品久久久久久久久久精品| 麻豆成人av视频| 国产极品粉嫩免费观看在线 | 亚洲欧美清纯卡通| 久久亚洲国产成人精品v| 蜜臀久久99精品久久宅男| 国产国拍精品亚洲av在线观看| 日韩成人av中文字幕在线观看| 国产亚洲精品久久久com| 久久人人爽人人片av| 高清欧美精品videossex| 国产av国产精品国产| 91久久精品国产一区二区三区| 国产一区亚洲一区在线观看| 免费久久久久久久精品成人欧美视频 | 九九爱精品视频在线观看| 亚洲国产毛片av蜜桃av| 亚洲无线观看免费| 精品熟女少妇av免费看| 91aial.com中文字幕在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | 街头女战士在线观看网站| 在线观看美女被高潮喷水网站| 欧美精品亚洲一区二区| 亚洲图色成人| 成年美女黄网站色视频大全免费 | 偷拍熟女少妇极品色| 欧美成人精品欧美一级黄| 国产高清不卡午夜福利| 搡老乐熟女国产| 国产精品蜜桃在线观看| 在线看a的网站| 在线观看国产h片| 狠狠精品人妻久久久久久综合| 精品国产一区二区三区久久久樱花| 亚洲国产欧美在线一区| 女的被弄到高潮叫床怎么办| 全区人妻精品视频| 日韩成人av中文字幕在线观看| 中文字幕制服av| 国产毛片在线视频| 亚洲精品,欧美精品| 91久久精品国产一区二区成人| 新久久久久国产一级毛片| 国产女主播在线喷水免费视频网站| 欧美性感艳星| 精品久久久噜噜| 精品熟女少妇av免费看| 亚洲经典国产精华液单| 中文字幕人妻丝袜制服| 日韩av在线免费看完整版不卡| 亚洲精品,欧美精品| 一区在线观看完整版| 街头女战士在线观看网站| 欧美bdsm另类| √禁漫天堂资源中文www| 丰满迷人的少妇在线观看| 精品久久久精品久久久| 女的被弄到高潮叫床怎么办| 熟女av电影| 中文天堂在线官网| 亚洲综合色惰| 91精品国产国语对白视频| 99热6这里只有精品| 成人综合一区亚洲| 黑丝袜美女国产一区| 看非洲黑人一级黄片| 美女xxoo啪啪120秒动态图| 午夜激情福利司机影院| 国产黄片美女视频| 久久ye,这里只有精品| 三上悠亚av全集在线观看 | 一本一本综合久久| 大香蕉久久网| 激情五月婷婷亚洲| 国产欧美日韩一区二区三区在线 | 精品久久久久久久久av| 国产乱人偷精品视频| 免费观看在线日韩| 菩萨蛮人人尽说江南好唐韦庄| 99九九线精品视频在线观看视频| √禁漫天堂资源中文www| 一级毛片电影观看| 成人18禁高潮啪啪吃奶动态图 | 国产精品.久久久| 久久久久久久久久人人人人人人| 一级,二级,三级黄色视频| 国产一区亚洲一区在线观看| 高清欧美精品videossex| 国产高清三级在线| 久久人妻熟女aⅴ| 又粗又硬又长又爽又黄的视频| 在现免费观看毛片| 精品一品国产午夜福利视频| 人人妻人人澡人人看| 曰老女人黄片| 九九在线视频观看精品| 最近中文字幕高清免费大全6| 成人美女网站在线观看视频| 国产伦精品一区二区三区视频9| 国产精品一区www在线观看| 女人久久www免费人成看片| 亚洲国产av新网站| 亚洲精品日本国产第一区| a级一级毛片免费在线观看| 人人妻人人澡人人看| 日韩 亚洲 欧美在线| 日韩熟女老妇一区二区性免费视频| 亚洲内射少妇av| 中文字幕亚洲精品专区| 日本黄大片高清| 啦啦啦在线观看免费高清www| 最后的刺客免费高清国语| 韩国av在线不卡| 男人狂女人下面高潮的视频| 亚洲国产精品国产精品| 丁香六月天网| 久久精品熟女亚洲av麻豆精品| 国产白丝娇喘喷水9色精品| 亚洲av日韩在线播放| 少妇被粗大猛烈的视频| 中文字幕av电影在线播放| 女的被弄到高潮叫床怎么办| 欧美老熟妇乱子伦牲交| a级毛片在线看网站| 中国国产av一级| 国产淫片久久久久久久久| 国产精品国产av在线观看| 国产精品偷伦视频观看了| 黑人高潮一二区| 成人无遮挡网站| 久久精品国产亚洲av天美| 老司机影院毛片| 亚洲伊人久久精品综合| 国产精品熟女久久久久浪| 美女中出高潮动态图| 少妇人妻 视频| 免费在线观看成人毛片| 欧美精品一区二区大全| 亚洲四区av| 日韩,欧美,国产一区二区三区| 啦啦啦在线观看免费高清www| 久久久久久久国产电影| 中文欧美无线码| 曰老女人黄片| 午夜福利在线观看免费完整高清在| 亚洲色图综合在线观看| 国产精品99久久久久久久久| 我要看日韩黄色一级片| 国产色爽女视频免费观看| 日韩中字成人| 婷婷色av中文字幕| 国产黄频视频在线观看| 一级黄片播放器| 一级毛片黄色毛片免费观看视频| 国产高清有码在线观看视频| 在线播放无遮挡| 国产亚洲5aaaaa淫片| 久久韩国三级中文字幕| 麻豆成人av视频| xxx大片免费视频| 国产91av在线免费观看| 寂寞人妻少妇视频99o| 有码 亚洲区| 极品少妇高潮喷水抽搐| 久久久亚洲精品成人影院| 狂野欧美激情性xxxx在线观看| 一个人免费看片子| 日日撸夜夜添| 黑人巨大精品欧美一区二区蜜桃 | 岛国毛片在线播放| 综合色丁香网| 麻豆成人av视频| 欧美国产精品一级二级三级 | 久久午夜综合久久蜜桃| a级一级毛片免费在线观看| av免费在线看不卡| 乱人伦中国视频| 午夜av观看不卡| 久久精品夜色国产| 在线观看av片永久免费下载| 简卡轻食公司| 国产女主播在线喷水免费视频网站| 国模一区二区三区四区视频| 看十八女毛片水多多多| 涩涩av久久男人的天堂| 少妇人妻久久综合中文| 九九爱精品视频在线观看| 亚洲精品视频女| 啦啦啦中文免费视频观看日本| 男人狂女人下面高潮的视频| 日本欧美视频一区| 三级经典国产精品| 亚洲av国产av综合av卡| 午夜久久久在线观看| 丰满少妇做爰视频| 内地一区二区视频在线| 久久久久久久久久成人| 亚洲欧美精品专区久久| 免费大片18禁| 一级av片app| 国产在线免费精品| 午夜免费鲁丝| 亚洲无线观看免费| 久久人人爽av亚洲精品天堂| 熟女av电影| 多毛熟女@视频| 国产日韩欧美在线精品| 国产一区亚洲一区在线观看| 国产一区有黄有色的免费视频| www.色视频.com| 成年av动漫网址| 国产真实伦视频高清在线观看| 麻豆乱淫一区二区| 蜜臀久久99精品久久宅男| 狂野欧美白嫩少妇大欣赏| 两个人的视频大全免费| 丝袜脚勾引网站| 亚洲欧美一区二区三区黑人 | 久久女婷五月综合色啪小说| 精品卡一卡二卡四卡免费| 欧美激情极品国产一区二区三区 | 一级a做视频免费观看| 99热这里只有是精品50| 中文精品一卡2卡3卡4更新| 久久久a久久爽久久v久久| 亚洲国产精品专区欧美| 在线观看人妻少妇| 国产一区亚洲一区在线观看| 精品一区二区三卡| 国产亚洲欧美精品永久| 涩涩av久久男人的天堂| 久久热精品热| 中文精品一卡2卡3卡4更新| 精品一区二区三区视频在线| 国产精品麻豆人妻色哟哟久久| 久久久国产精品麻豆| 26uuu在线亚洲综合色| 免费黄色在线免费观看| 欧美亚洲 丝袜 人妻 在线| 少妇 在线观看| 欧美三级亚洲精品| 在线 av 中文字幕| 永久免费av网站大全| 免费观看的影片在线观看| 国产av精品麻豆| 六月丁香七月| 黑丝袜美女国产一区| 精品卡一卡二卡四卡免费| 国产黄片美女视频| 最后的刺客免费高清国语| 久久久久久久久久人人人人人人| 国产在视频线精品| 男女免费视频国产| 乱系列少妇在线播放| av网站免费在线观看视频| 亚洲国产色片| 日韩人妻高清精品专区| 国产乱来视频区| 亚洲真实伦在线观看| 麻豆精品久久久久久蜜桃| 免费高清在线观看视频在线观看| 超碰97精品在线观看| 欧美日韩在线观看h| 男人狂女人下面高潮的视频| 精品国产乱码久久久久久小说| 国产欧美日韩一区二区三区在线 | 久久97久久精品| 黄色配什么色好看| videossex国产| 欧美最新免费一区二区三区| 日韩免费高清中文字幕av| 夜夜看夜夜爽夜夜摸| 久久国内精品自在自线图片| 高清黄色对白视频在线免费看 | 一区在线观看完整版| 91精品国产九色| 久久国产乱子免费精品| 亚洲精品国产av蜜桃| 视频区图区小说| 国产在线一区二区三区精| 色婷婷av一区二区三区视频| 国产色爽女视频免费观看| 日韩中文字幕视频在线看片| 国产伦精品一区二区三区四那| 免费av中文字幕在线| 亚洲人成网站在线观看播放| 欧美日韩视频高清一区二区三区二| 欧美激情国产日韩精品一区| 国产成人a∨麻豆精品| av天堂中文字幕网| 男人和女人高潮做爰伦理| 纵有疾风起免费观看全集完整版| 少妇人妻精品综合一区二区| 成人黄色视频免费在线看| 麻豆乱淫一区二区| 亚洲av日韩在线播放| 9色porny在线观看| 99热这里只有是精品在线观看| 国产成人freesex在线| 免费不卡的大黄色大毛片视频在线观看| 啦啦啦啦在线视频资源| 高清黄色对白视频在线免费看 | 精品少妇黑人巨大在线播放| 亚洲内射少妇av| 亚洲av在线观看美女高潮| 欧美精品一区二区免费开放| 亚洲欧美中文字幕日韩二区| 久久国产精品大桥未久av | 一个人看视频在线观看www免费| 国产男女内射视频| 极品少妇高潮喷水抽搐| 一区在线观看完整版| 在线免费观看不下载黄p国产| 色网站视频免费| 国产av国产精品国产| 国产淫语在线视频| 中国国产av一级| 欧美日本中文国产一区发布| 午夜老司机福利剧场| 午夜福利在线观看免费完整高清在| 国产淫片久久久久久久久| freevideosex欧美| 纯流量卡能插随身wifi吗| 人人妻人人澡人人看| 在线观看人妻少妇| 97在线视频观看| 人人妻人人爽人人添夜夜欢视频 | 日韩成人伦理影院| 在线观看免费高清a一片| 久久久久久久久久人人人人人人| 国产淫片久久久久久久久| 婷婷色麻豆天堂久久| av国产精品久久久久影院| 欧美日本中文国产一区发布| 午夜福利影视在线免费观看| 91午夜精品亚洲一区二区三区| 青春草国产在线视频| 欧美bdsm另类| 午夜影院在线不卡| a级毛片免费高清观看在线播放| 国产淫语在线视频| 丰满人妻一区二区三区视频av| 久久精品国产a三级三级三级| 久久精品久久久久久噜噜老黄| 国产欧美日韩一区二区三区在线 | 亚洲精品国产av成人精品| 日韩欧美 国产精品| 狂野欧美激情性bbbbbb| 我的老师免费观看完整版| 国产老妇伦熟女老妇高清| 亚洲一区二区三区欧美精品| 女的被弄到高潮叫床怎么办| 欧美日韩av久久| 亚洲精品,欧美精品| 成年av动漫网址| av又黄又爽大尺度在线免费看| 少妇的逼水好多| 亚洲色图综合在线观看| 日本-黄色视频高清免费观看| 狂野欧美激情性bbbbbb| 看非洲黑人一级黄片| 亚洲国产毛片av蜜桃av| 青春草亚洲视频在线观看| 女性生殖器流出的白浆| 麻豆成人午夜福利视频| 少妇人妻一区二区三区视频| 国产亚洲av片在线观看秒播厂| √禁漫天堂资源中文www| 内射极品少妇av片p| 成人无遮挡网站| 在线观看av片永久免费下载| 纯流量卡能插随身wifi吗| 各种免费的搞黄视频| 久久毛片免费看一区二区三区| 国精品久久久久久国模美| 免费观看av网站的网址| 韩国高清视频一区二区三区| 99热网站在线观看| 国产成人午夜福利电影在线观看| 日本av手机在线免费观看| 欧美区成人在线视频| 狂野欧美激情性bbbbbb| 日韩不卡一区二区三区视频在线| 欧美97在线视频| 天堂8中文在线网| 免费av不卡在线播放| 国产淫片久久久久久久久| 深夜a级毛片| 免费黄频网站在线观看国产| 日韩强制内射视频| av福利片在线观看| 亚洲欧美日韩另类电影网站| 99久久精品热视频| 久久这里有精品视频免费| 啦啦啦中文免费视频观看日本| 免费观看在线日韩| 热re99久久精品国产66热6| 七月丁香在线播放| 中文资源天堂在线| 亚洲综合精品二区| 国产黄频视频在线观看| 中文字幕制服av| 亚洲av成人精品一二三区| av在线app专区| 国产亚洲午夜精品一区二区久久| 亚洲中文av在线| 简卡轻食公司| 又大又黄又爽视频免费| 亚洲丝袜综合中文字幕| 黄色毛片三级朝国网站 | 国产精品伦人一区二区| 亚洲精品色激情综合| 午夜福利,免费看| 午夜福利影视在线免费观看| 国产黄频视频在线观看| 又黄又爽又刺激的免费视频.| 国产精品女同一区二区软件| 九九在线视频观看精品| 久久99热6这里只有精品| 少妇人妻久久综合中文| 久久亚洲国产成人精品v| 人妻人人澡人人爽人人| 亚洲欧美精品专区久久| 午夜影院在线不卡| av国产精品久久久久影院| 永久免费av网站大全| 纯流量卡能插随身wifi吗| 九色成人免费人妻av| 国产成人精品福利久久| 日韩 亚洲 欧美在线| 亚洲怡红院男人天堂| a级一级毛片免费在线观看| av卡一久久| 国产永久视频网站| 国产在线一区二区三区精| 成人影院久久| 一级毛片黄色毛片免费观看视频| 另类精品久久| 又粗又硬又长又爽又黄的视频| 黄色配什么色好看| 日韩 亚洲 欧美在线| 国产一区二区在线观看av| kizo精华| 久久国内精品自在自线图片| 日本午夜av视频| 久久婷婷青草| 精品卡一卡二卡四卡免费| 国产极品天堂在线| 乱码一卡2卡4卡精品| 青青草视频在线视频观看| 亚洲av日韩在线播放| 久久久久久久国产电影| 亚洲精品国产成人久久av| 插逼视频在线观看| 一级av片app| 九色成人免费人妻av| 亚洲av二区三区四区| 日韩一区二区三区影片| 各种免费的搞黄视频| 国产成人freesex在线| 欧美日韩av久久| 极品人妻少妇av视频| 午夜福利,免费看| 热99国产精品久久久久久7| av福利片在线| 在线观看一区二区三区激情| 噜噜噜噜噜久久久久久91| 国产精品.久久久| 在线免费观看不下载黄p国产| av免费在线看不卡| 国产在线视频一区二区| 亚洲欧美一区二区三区国产| 亚洲av.av天堂| 亚洲,欧美,日韩| 丝袜脚勾引网站| 亚洲国产精品专区欧美| 亚洲在久久综合| 亚洲国产色片| www.色视频.com| 内射极品少妇av片p| 精品人妻熟女av久视频| 亚洲av.av天堂| 精品人妻熟女毛片av久久网站| 精品一区二区三区视频在线| 菩萨蛮人人尽说江南好唐韦庄| 麻豆成人av视频| 黄色视频在线播放观看不卡| 亚洲av.av天堂| 国产一区二区三区综合在线观看 | 男人和女人高潮做爰伦理| 成人国产麻豆网| 中文资源天堂在线| 国产淫片久久久久久久久| 日本vs欧美在线观看视频 | 高清欧美精品videossex| 精品国产乱码久久久久久小说| 欧美精品高潮呻吟av久久| 久久午夜综合久久蜜桃| 三级国产精品欧美在线观看| 久久精品国产a三级三级三级|