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

    4個竹種高生長模型的比較

    2018-10-25 02:05:38凡美玲方水元馮俊嬌張春霞梁大洪時培建
    竹子學報 2018年2期
    關鍵詞:竹種曲率速率

    凡美玲,方水元,馮俊嬌,張春霞,梁大洪,時培建

    (1.南京林業(yè)大學南方現(xiàn)代林業(yè)協(xié)同創(chuàng)新中心,竹類研究所,江蘇 南京 210037;2.貴州省桐梓縣林業(yè)局,貴州 桐梓 563200)

    竹子屬于禾本科竹亞科植物,竹林為我國南方重要的森林資源之一,其生長不同于木本植物,生長速度快、種群密度高使得竹子得以在我國南方廣泛分布[1]。用合理的模型對植物生長數(shù)據(jù)進行擬合,能夠模擬植物生長的動態(tài)變化,這將有益于探索和比較生物個體的生長模式[2]。上世紀以來,相繼許多用于描述溫度變化對外溫生物發(fā)育速率的影響非線性模型被構建,許多模型被運用于擬合動植物的動態(tài)生長。澳大利亞學者Ratkowsky等[3]提出了平方根模型,此模型最初被用于描述細菌的種群增長速率,后被擴展到對節(jié)肢動物溫發(fā)育速率的研究中。Shi等發(fā)現(xiàn)去平方根模型和除平方根的Ratkowsky模型對昆蟲的溫發(fā)育數(shù)據(jù)的擬合效果較好[4],該模型是否能夠用于植物生長發(fā)育有待研究。Yin等[5]首次提出beta模型用于描述溫度對農(nóng)作物的影響,Yin等[6]通過對beta模型的積分提出了一個新模型(BSE),該模型能預測植物達到最大生物量的時間,為農(nóng)作物的生長培育提供理論依據(jù)。Shi等[7]用BSE模型與其它4個模型對6種農(nóng)作物的干重進行擬合對比,發(fā)現(xiàn)BSE模型擬合效果優(yōu)于其它4個模型,該模型能適應于描述多種農(nóng)作物的生長。Shi等用BSE模型[2]和LRF模型[8]對6種作物的干重、4個竹種高生長、2種樹木原始年輪數(shù)據(jù)和16個動物的鮮重數(shù)據(jù)進行擬合比較,擬合隨著時間推移各個動植物生物量(或者體型大小)指標變化的數(shù)據(jù),發(fā)現(xiàn)2種增長方程與實際上動植物的觀測數(shù)據(jù)吻合,且可以通過這2個模型預測動植物達到最大生物量和最高高度的時間,揭示了這2種模型功效的普適性[2]。

    唐紅等[9]2014年對湖南省植物園的圣音竹(Phyllostachysedulisf.tubaeformis)高生長研究發(fā)現(xiàn):用Logistic方程能很好的擬合幼竹平均高生長數(shù)據(jù),晝夜生長量可達16.48 cm。周芳純[1]對毛竹高生長與出土時間關系研究發(fā)現(xiàn):兩者關系呈“S”形曲線,可用Logistic方程擬合。曾慶南等[10]對江西林業(yè)科學院竹種園花毛竹(Ph.edulis‘Tao Kiang’)的高生長節(jié)律研究表明Logistic方程同樣能夠較好地擬合竹高生長過程。竹類植物在幼竹高生長時期的生長速度迅速,不同竹種每天的高生長量不同,通過對4個幼竹高生長速度的比較,以期為4個竹種的栽培管理提供理論依據(jù)。以往用Logistic方程對竹子的高生長擬合較多,本研究旨在比較不同增長方程對竹子高生長數(shù)據(jù)的擬合效果的優(yōu)劣,用模型擬合生物個體間不同的生長模式,并對不同模型的適用性和擬合效果進行比較,探討適用于擬合竹類植物高生長的最優(yōu)方程。

    1 材料與方法

    1.1 研究區(qū)概況

    研究地點位于南京林業(yè)大學南京本部校園內(nèi)(32°04′N,118°48′E),海拔約25 m,屬于亞熱帶季風氣候,研究地年平均氣溫為(15.6±0.9)℃,年平均降雨量為(1 045±242.7)mm左右,該地季節(jié)分明,每年6-7月的梅雨季節(jié)約22 d(1986-2016年氣候數(shù)據(jù)來源于中國氣象數(shù)據(jù)網(wǎng))。

    1.2 研究材料

    研究試驗材料為4個竹種的竹筍-幼竹,從竹筍露出地面1~2 cm為基準,開始對竹筍高度進行跟蹤調(diào)查,于2016年春季選擇4個竹種(每個竹種4株)進行跟蹤調(diào)查16株幼竹高生長。竹種分別為紅殼竹(Ph.iridescens)、美竹(Ph.mannii)、唐竹(Sinobambusatootsik)、斑苦竹(Pleioblastusmaculatus),發(fā)筍期均為春季。

    1.3 研究方法

    研究調(diào)查日期從2016年4月4日開始,直至6月25日結束調(diào)查。調(diào)查中,選用卷尺、可伸縮魚竿測量幼竹高度(每次測定幼竹全高,即為每個幼竹從地面到幼竹的最高處),調(diào)查初期為每天調(diào)查1次,后期直至幼竹基本完成高生長后,調(diào)查周期為2-3 d,調(diào)查時間均為中午12點,直至幼竹高度值穩(wěn)定,結束調(diào)查。

    研究共選擇3個非線性模型對4個竹種高生長數(shù)據(jù)進行擬合,每個竹種均選擇4個個體,每個模型都對16個個體的高生長數(shù)據(jù)進行擬合,3個模型的參數(shù)均為3個,比較3個模型對16個個體高生長數(shù)據(jù)的擬合優(yōu)度及模型參數(shù)的線性近似表現(xiàn),3個模型如下所示。

    1.3.1von Bertalanffy模型 von Bertalanffy模型[11]首次由von Bertalanffy提出,2001年West等[12]提出了一個新的個體增長模型(OGM)。2014年Shi等將OGM模型用于擬合3種淡水魚身體長度的變化,指出不同的物種應該有不同的新陳代謝率指數(shù)取值[13]。研究使用的von Bertalanffy模型,其具體形式如下:

    h=H·(1-exp(-3/4·K(t-t0)))4/3

    此處,h表示在時間t處生物體一維幾何參數(shù)的大小(長度或者高度)。k為常數(shù);t0亦是一個常量,表示是當竹子高度為0時的起始時間;參數(shù)H表示漸近線趨于最大值的一維幾何參數(shù)(即最大長度或者最大高度)。

    1.3.2BSM模型 Yin等[5]提出一個beta模型,最初是為了描述溫度對玉米發(fā)育速率的影響。Yin等[6]在beta模型的基礎上建立了1個3個參數(shù)的個體生長模型BSE,該模型生物量為0時對應的時間被定義為0。該模型具體形式如下:

    h=copt·t·(2tmax-topt-t)/(2tmax-topt)(t/topt)topt/(tmax-topt)

    此處,h表示在時間t處生物體一維幾何參數(shù)的大小(長度或者高度)。模型共有3個待擬合參數(shù),tmax表示幼竹生長到最高高度的時間,topt表示最適幼竹生長的時間(生長速率最大時),copt表示幼竹最大生長速率的值。

    1.3.3Logistic模型 Logistic模型描述動物和植物的生物量和身體長度數(shù)據(jù)運用最廣泛的模型,曲線的對稱性為該模型最典型的特點。該模型具體形式如下:

    h=H/(1+exp(b-rt))

    此處,共有3個參數(shù),h表示隨時間t變化的幼竹高度,H代表待擬合的幼竹最高高度,r代表幼竹的瞬時增長率,b代表一個常數(shù)。

    1.4 參數(shù)擬合方法及模型評估的指標

    研究使用R語言自編程序?qū)τ字窀呱L進行擬合,使用進化遺傳算法[14]來執(zhí)行對模型參數(shù)的估計,該方法是針對于解決求一個函數(shù)最大值的優(yōu)化問題。研究使用決定系數(shù)(R2)和曲率評判模型的優(yōu)劣,R2越大擬合優(yōu)度越好[15]。在非線性回歸中曲率分為均方根內(nèi)稟曲率和均方根參數(shù)效果曲率(下文將省略“均方根”),其中內(nèi)稟曲率是受到模型本身制約的無法通過對模型參數(shù)的的再參數(shù)化來減小,而參數(shù)效果曲率則可以通過對模型參數(shù)的再參數(shù)化來實現(xiàn)變化,參數(shù)效果曲率相對于臨界曲率越小,表明模型參數(shù)近似線性的表現(xiàn)越好,模型越合理[16-17]。

    2 結果與分析

    2.1 3個模型的擬合參數(shù)比較

    研究共用3種模型分別對4個竹種不同個體的幼竹高生長數(shù)據(jù)進行擬合,計算結果(R2)如表1所示,從表1中的擬合結果可以看出:3個模型分別對16株幼竹個體高生長數(shù)據(jù)擬合效果都很好,R2最小也可以達到0.981 66,為von Bertalanffy模型對紅殼竹4號個體的擬合結果,而擬合效果最好的出現(xiàn)在BSM模型對斑苦竹1號個體的擬合效果,R2最大值可達0.999 28。對擬合結果R2分析發(fā)現(xiàn),Logstic模型擬合后的R2的計算值有6個數(shù)值達0.998 00,BSM模型有3個數(shù)值達0.998 00,但von Bertalanffy模型擬合后R2最大為0.996 52,沒有數(shù)值達0.998 00。不同模型對不同竹種的擬合結果也不同,Logstic模型和BSM模型對斑苦竹擬合效果最佳,但von Bertalanffy模型對美竹(Ph.mannii)擬合效果較差,該模型2014年被應用于描述淡魚體長變化,研究將其應用于植物,效果欠缺,有待改善模型結構及參數(shù)使其適應性更廣泛。綜上所述,Logistic模型對16株幼竹個體高生長數(shù)據(jù)擬合效果最優(yōu),BSM模型擬合效果較好,von Bertalanffy模型的擬合結果較差, Logistic模型對于諸多竹種幼竹高生長的擬合效果最好,可以被廣泛運用于動植物個體增長研究。

    表1 3個模型對4個竹種16株個體高生長數(shù)據(jù)集的R2擬合結果Tab.1 Coefficient of determination (R2) obtained from using the three models to fit the height growth data sets of sixteen individual shoots of four bamboo species

    說明:編號Ⅰ、Ⅱ、Ⅲ、Ⅳ分別代表紅殼竹、美竹、唐竹和斑苦竹,編號1、2、3、4代表每個竹種的不同個體。
    The number of Ⅰ、Ⅱ、Ⅲ、Ⅳ representPh.iridescens、Ph.Mannii、Sinobambusatootsik、Pleioblastusmaculatus,and the number of 1、2、3、4 represent different individuals of each bamboo species.

    2.2 3個模型參數(shù)曲率的比較

    表2為3個模型對16個幼竹高生長數(shù)據(jù)擬合后的曲率對比,內(nèi)稟曲率和參數(shù)效果曲率用于評估模型的非線性表現(xiàn)和參數(shù)結構的合理性。此處,曲率指標共有3個:內(nèi)稟曲率(RMS IC)、參數(shù)效果曲率(RMS PEC)和臨界曲率(CC),當模型擬合后的內(nèi)稟曲率和參數(shù)效果曲率都小于臨界曲率,表明該模型近似線性的表現(xiàn)越強,參數(shù)結構更為合理。由表2的擬合結果可以看出:3個模型Logstic模型擬合結果中,內(nèi)稟曲率和參數(shù)效果曲率均小于臨界曲率,表明該模型參數(shù)的線性近似表現(xiàn)較好,模型參數(shù)設置合理;對于BSM模型,內(nèi)稟曲率和參數(shù)效果曲率均小于臨界曲率,與Logstic模型擬合結果相比,BSE模型擬合得到的內(nèi)稟曲率和參數(shù)效果曲率均大于Logstic模型,表明Logstic模型的擬合結果優(yōu)于BSM模型;von Bertalanffy模型擬合結果中,內(nèi)稟曲率均小于臨界曲率,但有4個參數(shù)效果曲率大于臨界曲率,參數(shù)效果曲率均>BSM模型>Logstic模型。綜上所述,Logistic模型對16株幼竹個體高生長數(shù)據(jù)擬合優(yōu)度最高,參數(shù)設置更為合理,其次為BSM模型,von Bertalanffy模型的擬合優(yōu)度最低[17]。

    2.3 3個模型的擬合曲線比較

    圖1為BSM模型、von Bertalanffy模型和Logistic模型3種不同模型對斑苦竹中1號幼竹高生長數(shù)據(jù)的擬合結果(3個模型對該個體擬合效果最好)。研究中的斑苦竹種植于南京林業(yè)大學竹種園內(nèi),1號斑苦竹于2016年4月29日破土而出,直至2016年6月25日完成其高生長,最高達12 m。開始生長速度較慢,最快速度出現(xiàn)于6月6日,此時增長速率達到72 cm·d-1,隨后生長速度減慢,其生長曲線呈典型的“慢-快-慢”趨勢。

    表2 3個模型對16個個體高生長數(shù)據(jù)集擬合后得到的曲率比較Tab.2 Curvatures obtained from using three models to fit the height growth data sets of sixteen bamboo shoots

    研究共選擇了3個非線性模型對同一竹種不同個體高生長數(shù)據(jù)進行擬合,3個模型都設置了3個待擬合參數(shù)(BSM模型中tmin設定為0)。von Bertalanffy模型是一個通用型模型,可以用于諸多生物的生長,研究該模型新陳代謝率指數(shù)選用的為3/4,但von Bertalanffy模型對幼竹高度的擬合曲線與實際數(shù)據(jù)有很大出入(圖1c和圖1d),該擬合曲線在前10 d左右零生長,隨后生長速度驟增,但在曲線末端未見曲線平緩,有繼續(xù)增加的趨勢,這與實際觀測數(shù)據(jù)不符,幼竹在完成高生長后,高度會趨于一個穩(wěn)定值,其下方對幼竹生長速率的擬合曲線同樣可以看出,該幼竹前期生長速率為0,隨后在10 d左右達到最大生長速率,達到最大生長速率共用24 d,其后幼竹繼續(xù)生長,生長速率線始終沒有趨于0,且整條生長速率曲線不對稱,擬合后決定系數(shù)計算結果為0.988 45,該模型對其幼竹高生長擬合效果差。BSM模型最大的特點是能夠包容對稱和不對稱曲線,能準確估計植物到達最大生物量的時間。圖1a和圖1b展現(xiàn)了BSM模型對斑苦竹中1號幼竹高生長數(shù)據(jù)及生長速率的擬合結果,該擬合結果表現(xiàn)為前期生長緩慢,隨后進入高生長盛期,曲線最后趨于平緩,幼竹高度達到最大高度后稍有下降(幼竹開始伸枝展葉),最后穩(wěn)定于一個值,擬合曲線與幼竹實際高度重合率高。從該模型對其生長速率擬合曲線可以發(fā)現(xiàn):高生長階段達到最大生長速率的時間為距離出筍35 d,52 d后結束高生長,該生長速率曲線左側寬于右側表明,竹子前期到達最大生長速率的時間略多于后期生長的時間,擬合后決定系數(shù)計算結果可達0.999 28,該模型對其擬合效果較優(yōu)。Logistic模型被廣泛應用于描述動植物的生物量,從圖1e和圖1f可以看出,幼竹高生長呈“慢-快-慢”趨勢,擬合曲線與幼竹高生長數(shù)據(jù)高度重合,與BSM模型不同的是,該模型對生長速率的擬合曲線更平緩,且曲線兩側對稱,幼竹最大生長速率時間為距出筍33 d,與實際數(shù)據(jù)相符,擬合后決定系數(shù)計算結果能夠達到0.999 01,綜上所述,Logistic模型擬合結果優(yōu)于BSM模型優(yōu)于von Bertalanffy模型,在本研究中,3個模型對不同幼竹個體高生長數(shù)據(jù)的擬合結果表明,Logistic模型為最優(yōu)模型。

    圖1 用3種模型對斑苦竹(Pleioblastus maculatus)1號個體高生長數(shù)據(jù)的擬合結果Fig.1 Fitted results of the height growth for No.1 bamboo of Pleioblastus maculatus by using three models說明:a、b為BSM模型;c、d為von Bertalanffy模型;e、f為Logistic模型。a、b are the BSM model; c、d are the von Bertalanffy model; and e、f are Logistic model

    3 討論

    竹子每年的發(fā)筍時間主要受當年的溫度和降雨量等因素影響,按發(fā)筍期發(fā)筍數(shù)量可分為初期、中期和末期,末期筍常因營養(yǎng)不足或蟲害退筍[18-20]。實驗選擇是4個竹種發(fā)筍期差異主要因為竹種差異造成,選擇發(fā)筍前期到中期的筍為研究對象,調(diào)查發(fā)現(xiàn)紅殼竹與美竹發(fā)筍日期與高生長期相近,于4月初開始發(fā)筍,6月伊始結束高生長,唐竹和斑苦竹發(fā)筍期略晚,分別為4月中期和末期,高生長期約為40-55 d。

    研究使用的von Bertalanffy模型,參數(shù)=3/4,但總體擬合效果差,同時也嘗試參數(shù)=2/3,擬合優(yōu)度低于前者,2014年Shi等[2]用該模型和OGM模型同時對3種淡魚體長數(shù)據(jù)進行擬合,他們發(fā)現(xiàn)當OGM模型中的參數(shù)α越接近1,模型的擬合效果越好,與本研究的結果一致。本研究中,使用的BSM模型中限制δ=1、tmin為0,但并不會影響該模型描述數(shù)據(jù)的有效性,結果表明BSM模型對4種幼竹高生長數(shù)據(jù)擬合效果較好,對幼竹生長到最大生物量的時間估計與實際相符,Shi等[2]用2種生長模型分別對動植物生長進行擬合,作者指出:如果可以預先確定動植物增長的起始時間,那對于BSM模型的參數(shù)δ可以不限制為1,可以放寬δ作為一個參數(shù),這樣將幫助對比動植物種間增長模式的差異。對于常用的生長模型Logistic模型,擬合結果同前人研究一致,對本研究中的4種竹種個體高生長數(shù)據(jù)擬合效果好,Shi等13用4種模型同時對相同數(shù)據(jù)擬合比較發(fā)現(xiàn),Logistic模型對于描述動植物生物量和身體長度都是通用的,且擬合優(yōu)度都優(yōu)于其它模型,是最適合描述動植物的生長模型。Richards[21]提出一個廣義的Logstic模型,即Richards模型,該模型更加靈活,不再限制生長速率曲線的對稱性,且該模型對于在其最小值到最大值之間的任何值都能有其最大斜率,但該模型也存在一些爭議,Ratkowsky指出該模型對系數(shù)估計值不準確[21],Birch[23]也指出該模型不能反映一些蔬菜和牧草的生長,如黑麥草(LoliumperenneL.)在低密度時就能達到其最大凈產(chǎn)量,與實際不符。研究表明:許多非線性模型不僅適用于溫度對動物生長速度和內(nèi)稟增長率的影響,對于植物的生長也是適用的。

    最初一些非線性模型被提出用于描述溫度對昆蟲發(fā)育速率的影響,如表現(xiàn)型模型[24]、Taylor模型[25]和Ratkowsky模型[4]等等,這些模型可以幫助預測一些重要參數(shù),對農(nóng)業(yè)生產(chǎn)及科學研究有著重要意義,但因為模型結構設置及參數(shù)選擇不同,對不同數(shù)據(jù)的擬合優(yōu)度不同,選擇一個通用性的最優(yōu)模型具有重要意義。評價模型擬合優(yōu)度的好壞可以從幾個方面評估,模型結構的復雜程度、模型適用性以及模型預測曲線與實際是否相符等[2]。隨后,王如松等將數(shù)學模型用于描述溫度變化對昆蟲內(nèi)稟增長率的影響,擬合結果顯示相關系數(shù)高,擬合結果較優(yōu),指出可以利用數(shù)學模型研究一些與溫度相關的現(xiàn)象[26]。2016年Shi等提出BSM的通用表達式(即生長發(fā)育初始時間不為0時的BSM模型)能夠擬合植物的動態(tài)生長,并能幫助預測植物生長到最大生長量的時間,這為森林資源的合理利用提供了參考價值[7]。Shi等[13]將1種個體增長模型運用于作物的生長,結果表明Logistic模型能夠更好地擬合動植物生長的數(shù)據(jù),這與本研究結果一致。本研究中,BSM模型與Logistic模型對幼竹生長到最大生物量時間的預測相近,與實際數(shù)據(jù)吻合,而von Bertalanffy模型預測幼竹最大生物量的時間提前10 d,正如圖1d所展示的生長速率曲線左偏。同樣,Logistic模型擬合的生長曲線同前人研究一樣,曲線兩側關于對稱軸對稱,對數(shù)據(jù)缺乏靈活性[27],但總體Logistic模型適合描述植物生長過程中生物量的變化趨勢。

    研究使用3個非線性模型對不同竹種幼竹高生長數(shù)據(jù)進行擬合,將以往描述溫度對昆蟲發(fā)育速率影響的模型用于擬合植物的生長,將自變量溫度調(diào)換為時間,結果表明,生長模型也能很好地描述植物的生長[2]。調(diào)查發(fā)現(xiàn),不同竹種前期與盛期筍生長為成竹的概率較大,末期筍易發(fā)生退筍現(xiàn)象,所以在選擇留筍養(yǎng)竹時,應選擇初期和盛期萌發(fā)的筍為宜。利用3個模型對不同幼竹高生長數(shù)據(jù)進行擬合,4個竹種高生長均符合“S”型生長,符合“慢-快-慢”是生長規(guī)律,這為4個竹種的撫育管理提供了重要的理論依據(jù)。本研究發(fā)現(xiàn),Logistic模型對動植物生長的擬合效果最優(yōu),BSM模型與von Bertalanffy模型有待改進,這些生長模型的適用性有待進一步探索,生物有機體生長模式的比較也將成為一個研究熱點。

    猜你喜歡
    竹種曲率速率
    大曲率沉管安裝關鍵技術研究
    一類雙曲平均曲率流的對稱與整體解
    “化學反應的速率與限度”知識與能力提升
    半正迷向曲率的四維Shrinking Gradient Ricci Solitons
    速度和速率有什么不同
    河南淅川竹博園竹種引種試驗初報
    牡竹屬3個竹種地上部分生物量及特性研究
    不同冷卻速率下低壓轉子鋼30Cr2Ni4MoV的凝固組織
    上海金屬(2014年5期)2014-12-20 07:58:39
    蓮心超微粉碎提高有效成分的溶出速率
    Esn+1中具有至多兩個不同主曲率的2-調(diào)和超曲面
    午夜福利在线观看免费完整高清在 | 欧美高清成人免费视频www| 亚洲精品在线观看二区| 中文资源天堂在线| 国产精品野战在线观看| 国产aⅴ精品一区二区三区波| 精品久久久久久久久久免费视频| 一本精品99久久精品77| 国产视频一区二区在线看| 午夜免费成人在线视频| 亚洲av日韩精品久久久久久密| 日本欧美国产在线视频| 日本三级黄在线观看| 国产精品人妻久久久久久| 精品一区二区三区人妻视频| 国产av不卡久久| 精品久久久久久久末码| 别揉我奶头 嗯啊视频| 日本熟妇午夜| 国产真实伦视频高清在线观看 | 欧美国产日韩亚洲一区| 一级黄色大片毛片| 国产精品免费一区二区三区在线| 久久久久国内视频| 欧美一级a爱片免费观看看| 国产精品亚洲一级av第二区| 97超视频在线观看视频| 日韩亚洲欧美综合| 免费观看精品视频网站| 中文字幕av成人在线电影| 亚洲av成人av| 熟女人妻精品中文字幕| 中文字幕高清在线视频| 日韩欧美精品v在线| 中文字幕熟女人妻在线| 日本黄色片子视频| 国产亚洲精品av在线| 精品人妻偷拍中文字幕| 嫩草影院精品99| 51国产日韩欧美| 国产高清三级在线| 悠悠久久av| 88av欧美| 亚洲成人久久爱视频| 久久人妻av系列| 国产高清不卡午夜福利| 在线免费观看的www视频| 在线观看午夜福利视频| 国产 一区 欧美 日韩| 在线观看66精品国产| 精品久久久久久久久亚洲 | 九九在线视频观看精品| 国产欧美日韩精品亚洲av| 禁无遮挡网站| 免费观看精品视频网站| 最好的美女福利视频网| 深爱激情五月婷婷| 亚洲男人的天堂狠狠| 欧美人与善性xxx| 亚洲精品一区av在线观看| 九九热线精品视视频播放| 午夜爱爱视频在线播放| 中国美白少妇内射xxxbb| 久久久午夜欧美精品| 成人特级av手机在线观看| 国产精品一区www在线观看 | 国产黄色小视频在线观看| 欧美丝袜亚洲另类 | 亚洲欧美日韩卡通动漫| 国产av一区在线观看免费| 小说图片视频综合网站| 99久久中文字幕三级久久日本| 久久久久国内视频| 美女被艹到高潮喷水动态| 12—13女人毛片做爰片一| 国产av不卡久久| 熟女人妻精品中文字幕| 日韩国内少妇激情av| 美女大奶头视频| 亚洲狠狠婷婷综合久久图片| 乱人视频在线观看| 久久久久久久久久久丰满 | 一本久久中文字幕| 成年免费大片在线观看| 精品久久久久久久人妻蜜臀av| 国产精品久久久久久久电影| 免费看a级黄色片| 天堂√8在线中文| 91久久精品国产一区二区成人| 中文字幕精品亚洲无线码一区| 久久久久精品国产欧美久久久| 看片在线看免费视频| 亚洲午夜理论影院| 欧美最新免费一区二区三区| 村上凉子中文字幕在线| 欧美一区二区国产精品久久精品| 18禁裸乳无遮挡免费网站照片| 成年女人毛片免费观看观看9| 看黄色毛片网站| 12—13女人毛片做爰片一| 日韩欧美三级三区| 久久亚洲真实| xxxwww97欧美| 人妻丰满熟妇av一区二区三区| 黄色一级大片看看| 亚洲真实伦在线观看| 五月玫瑰六月丁香| 最后的刺客免费高清国语| 国产69精品久久久久777片| 国产精品野战在线观看| 两个人的视频大全免费| 免费av毛片视频| 成人午夜高清在线视频| 中文字幕精品亚洲无线码一区| 久久精品国产鲁丝片午夜精品 | 中文字幕人妻熟人妻熟丝袜美| 午夜久久久久精精品| 蜜桃亚洲精品一区二区三区| 国产大屁股一区二区在线视频| 搞女人的毛片| 12—13女人毛片做爰片一| 免费av毛片视频| 久久6这里有精品| 亚洲av免费高清在线观看| 免费一级毛片在线播放高清视频| 亚洲欧美日韩无卡精品| 少妇猛男粗大的猛烈进出视频 | 免费看av在线观看网站| 亚洲专区国产一区二区| 国产一区二区三区在线臀色熟女| 极品教师在线视频| 亚洲中文字幕日韩| 亚洲狠狠婷婷综合久久图片| 女人被狂操c到高潮| 亚洲av成人精品一区久久| 我的老师免费观看完整版| 色噜噜av男人的天堂激情| 两个人视频免费观看高清| 精品久久久久久久久久久久久| 成人特级黄色片久久久久久久| 亚洲国产欧洲综合997久久,| 久久草成人影院| 亚洲av二区三区四区| 精品乱码久久久久久99久播| 非洲黑人性xxxx精品又粗又长| 久久久久久久久久成人| 99久久中文字幕三级久久日本| 别揉我奶头~嗯~啊~动态视频| 色综合婷婷激情| 成人三级黄色视频| 1024手机看黄色片| 欧美日本视频| 国产蜜桃级精品一区二区三区| 国产亚洲精品久久久com| a级一级毛片免费在线观看| 无遮挡黄片免费观看| 狠狠狠狠99中文字幕| 久久精品国产亚洲av天美| 国产免费一级a男人的天堂| 露出奶头的视频| 亚洲精品456在线播放app | 熟女人妻精品中文字幕| 免费无遮挡裸体视频| 亚洲国产欧洲综合997久久,| 最新在线观看一区二区三区| 亚洲经典国产精华液单| 丰满乱子伦码专区| 久久久久精品国产欧美久久久| 欧美一区二区精品小视频在线| 偷拍熟女少妇极品色| 别揉我奶头 嗯啊视频| 91久久精品电影网| 免费大片18禁| 淫妇啪啪啪对白视频| 亚洲精品亚洲一区二区| 中国美女看黄片| 成人美女网站在线观看视频| 神马国产精品三级电影在线观看| 成人无遮挡网站| 窝窝影院91人妻| 亚洲精品成人久久久久久| 亚洲va在线va天堂va国产| 窝窝影院91人妻| www日本黄色视频网| 成人无遮挡网站| 真人做人爱边吃奶动态| 色综合亚洲欧美另类图片| 久99久视频精品免费| 免费看光身美女| 色吧在线观看| 蜜桃亚洲精品一区二区三区| 日韩一本色道免费dvd| 亚洲av第一区精品v没综合| 一本精品99久久精品77| www日本黄色视频网| 亚洲av免费高清在线观看| 国产成人一区二区在线| 深夜精品福利| 久久草成人影院| 午夜福利在线观看免费完整高清在 | 人妻夜夜爽99麻豆av| av在线天堂中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 一区二区三区高清视频在线| 男女之事视频高清在线观看| 国内精品久久久久精免费| 亚洲精品日韩av片在线观看| 无人区码免费观看不卡| 在线天堂最新版资源| av在线天堂中文字幕| 天堂影院成人在线观看| 少妇被粗大猛烈的视频| 精品无人区乱码1区二区| 国内精品一区二区在线观看| 舔av片在线| 乱码一卡2卡4卡精品| 免费看a级黄色片| 日本一二三区视频观看| 国产午夜精品论理片| 久久久久久伊人网av| 国产美女午夜福利| 又紧又爽又黄一区二区| 亚洲av熟女| 国产成人a区在线观看| 亚洲自拍偷在线| 亚洲va在线va天堂va国产| 亚洲无线观看免费| 美女大奶头视频| 深爱激情五月婷婷| 日韩欧美精品免费久久| 亚洲精品456在线播放app | 狂野欧美激情性xxxx在线观看| 乱系列少妇在线播放| 大又大粗又爽又黄少妇毛片口| 国产视频一区二区在线看| 亚州av有码| 男人和女人高潮做爰伦理| 日韩大尺度精品在线看网址| 国产高潮美女av| 男女边吃奶边做爰视频| 中文字幕精品亚洲无线码一区| 国产亚洲精品av在线| 午夜福利欧美成人| 久久亚洲真实| 少妇高潮的动态图| 亚洲熟妇中文字幕五十中出| 国产久久久一区二区三区| 嫩草影视91久久| 国产精品久久视频播放| 亚洲无线观看免费| 91久久精品国产一区二区三区| 中文在线观看免费www的网站| av中文乱码字幕在线| 亚洲最大成人手机在线| 午夜久久久久精精品| 久久人人精品亚洲av| 国产亚洲av嫩草精品影院| 老司机深夜福利视频在线观看| 床上黄色一级片| 九九爱精品视频在线观看| 亚洲专区中文字幕在线| 成年人黄色毛片网站| 伦理电影大哥的女人| 久久精品影院6| 女人被狂操c到高潮| 国产一区二区三区视频了| 亚洲av熟女| 免费av观看视频| 精品乱码久久久久久99久播| 99久久成人亚洲精品观看| 久久久午夜欧美精品| 午夜福利在线观看吧| 特级一级黄色大片| 露出奶头的视频| 亚洲在线观看片| 亚洲精品成人久久久久久| 热99re8久久精品国产| 欧美三级亚洲精品| 嫩草影院入口| 伦精品一区二区三区| 色吧在线观看| 日韩欧美国产在线观看| 国产蜜桃级精品一区二区三区| 国产单亲对白刺激| 日本与韩国留学比较| 别揉我奶头~嗯~啊~动态视频| .国产精品久久| 一夜夜www| 亚洲精华国产精华液的使用体验 | 国产一区二区在线观看日韩| 最新在线观看一区二区三区| 亚洲真实伦在线观看| 成人av在线播放网站| 国产高清不卡午夜福利| 中亚洲国语对白在线视频| 91狼人影院| 国产国拍精品亚洲av在线观看| 直男gayav资源| 精品人妻视频免费看| 无人区码免费观看不卡| 成人国产一区最新在线观看| 久久九九热精品免费| 亚洲av五月六月丁香网| 日本a在线网址| 午夜视频国产福利| 亚洲av成人av| 日韩大尺度精品在线看网址| 亚洲av成人av| 国内精品久久久久精免费| 他把我摸到了高潮在线观看| 亚洲精品乱码久久久v下载方式| 国产三级在线视频| 久久久久性生活片| 久久热精品热| 一进一出抽搐动态| 欧美一区二区精品小视频在线| 国产亚洲av嫩草精品影院| 亚洲精品一卡2卡三卡4卡5卡| 黄色视频,在线免费观看| 成人综合一区亚洲| 免费不卡的大黄色大毛片视频在线观看 | 尤物成人国产欧美一区二区三区| 男人舔女人下体高潮全视频| 老司机午夜福利在线观看视频| 国产av在哪里看| 一进一出抽搐动态| 男人舔奶头视频| 99在线视频只有这里精品首页| 欧美最新免费一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 色5月婷婷丁香| 狂野欧美激情性xxxx在线观看| 美女被艹到高潮喷水动态| 亚洲国产欧洲综合997久久,| 免费在线观看日本一区| 亚洲熟妇中文字幕五十中出| 成人综合一区亚洲| 亚洲精品粉嫩美女一区| 窝窝影院91人妻| 免费高清视频大片| 日韩大尺度精品在线看网址| 九色成人免费人妻av| 国产精品一区二区免费欧美| 亚洲18禁久久av| 国产在视频线在精品| 久久精品国产亚洲av香蕉五月| 亚洲久久久久久中文字幕| 黄片wwwwww| 欧美3d第一页| 欧美日韩精品成人综合77777| 欧美中文日本在线观看视频| 亚洲欧美日韩无卡精品| 精品午夜福利视频在线观看一区| 午夜a级毛片| 午夜福利成人在线免费观看| 一a级毛片在线观看| 老司机午夜福利在线观看视频| 无遮挡黄片免费观看| 最新在线观看一区二区三区| 美女免费视频网站| 欧美另类亚洲清纯唯美| 久久亚洲精品不卡| 日韩欧美国产在线观看| 亚洲精品色激情综合| 国产高清三级在线| 99在线人妻在线中文字幕| 亚洲欧美日韩高清在线视频| 99在线人妻在线中文字幕| 99精品在免费线老司机午夜| 国产黄a三级三级三级人| 午夜精品一区二区三区免费看| 国产精品国产三级国产av玫瑰| 国产又黄又爽又无遮挡在线| 亚洲熟妇中文字幕五十中出| 看片在线看免费视频| 欧美一区二区国产精品久久精品| 亚洲国产日韩欧美精品在线观看| 搡老妇女老女人老熟妇| 免费无遮挡裸体视频| 久久亚洲真实| 全区人妻精品视频| 久久午夜亚洲精品久久| 两个人的视频大全免费| 午夜激情福利司机影院| 久久久午夜欧美精品| 在线看三级毛片| 日本五十路高清| 日本撒尿小便嘘嘘汇集6| 亚洲精品在线观看二区| 亚洲精品456在线播放app | 久久久精品大字幕| 99在线人妻在线中文字幕| 日韩一区二区视频免费看| АⅤ资源中文在线天堂| av专区在线播放| 亚洲精品粉嫩美女一区| netflix在线观看网站| 18禁在线播放成人免费| 看免费成人av毛片| 国模一区二区三区四区视频| 国内精品美女久久久久久| 久久午夜福利片| 一个人看视频在线观看www免费| 美女高潮喷水抽搐中文字幕| 最近在线观看免费完整版| av在线亚洲专区| 一区二区三区激情视频| 嫩草影院精品99| 国产在线精品亚洲第一网站| 国产一区二区在线观看日韩| 久久久久性生活片| 成人毛片a级毛片在线播放| 一区二区三区免费毛片| 亚洲午夜理论影院| 欧美精品国产亚洲| 亚洲av电影不卡..在线观看| av天堂在线播放| 午夜福利视频1000在线观看| 婷婷亚洲欧美| 国内毛片毛片毛片毛片毛片| 日日夜夜操网爽| 色综合亚洲欧美另类图片| 久久精品国产99精品国产亚洲性色| 老熟妇仑乱视频hdxx| 午夜老司机福利剧场| a级一级毛片免费在线观看| 日韩强制内射视频| 亚洲国产精品合色在线| 久久久久九九精品影院| 免费观看在线日韩| 国产精品一区www在线观看 | 亚洲avbb在线观看| 免费不卡的大黄色大毛片视频在线观看 | 婷婷亚洲欧美| 中国美女看黄片| 日本黄大片高清| 午夜激情欧美在线| 少妇人妻一区二区三区视频| 别揉我奶头 嗯啊视频| 精品乱码久久久久久99久播| av专区在线播放| 欧美性感艳星| 少妇高潮的动态图| 亚洲最大成人手机在线| 神马国产精品三级电影在线观看| 亚洲av不卡在线观看| 一进一出好大好爽视频| 国产av麻豆久久久久久久| 三级男女做爰猛烈吃奶摸视频| 欧美zozozo另类| 日韩中文字幕欧美一区二区| 亚洲成人久久爱视频| 亚洲成人精品中文字幕电影| 九九爱精品视频在线观看| 亚洲av美国av| netflix在线观看网站| 久久6这里有精品| 九九在线视频观看精品| 麻豆国产av国片精品| av.在线天堂| 男女视频在线观看网站免费| 一个人看视频在线观看www免费| 久久精品夜夜夜夜夜久久蜜豆| 女人被狂操c到高潮| 波多野结衣高清无吗| bbb黄色大片| 小说图片视频综合网站| 午夜老司机福利剧场| 亚洲人成网站在线播放欧美日韩| 亚洲欧美清纯卡通| av.在线天堂| 高清毛片免费观看视频网站| 亚洲专区国产一区二区| 性插视频无遮挡在线免费观看| 欧美高清成人免费视频www| 精品午夜福利在线看| 久久精品国产鲁丝片午夜精品 | 成年女人毛片免费观看观看9| 国产一区二区三区av在线 | 久久久久久伊人网av| 老司机福利观看| 最好的美女福利视频网| 成人av在线播放网站| 成年女人看的毛片在线观看| 一边摸一边抽搐一进一小说| 亚洲经典国产精华液单| 欧美成人a在线观看| 内射极品少妇av片p| 欧美一区二区亚洲| 精品久久久久久久人妻蜜臀av| 在线免费观看的www视频| 人妻久久中文字幕网| 又粗又爽又猛毛片免费看| 欧美激情在线99| 成人永久免费在线观看视频| 深夜精品福利| 国产aⅴ精品一区二区三区波| eeuss影院久久| 美女高潮的动态| 欧美又色又爽又黄视频| 一进一出抽搐动态| 老熟妇仑乱视频hdxx| 亚洲国产色片| 国内精品美女久久久久久| 桃红色精品国产亚洲av| 在线天堂最新版资源| 国产成人影院久久av| 国产亚洲精品久久久久久毛片| 久久亚洲真实| 网址你懂的国产日韩在线| 毛片一级片免费看久久久久 | 国产高清激情床上av| 日韩欧美在线二视频| 成人国产麻豆网| 蜜桃亚洲精品一区二区三区| 婷婷色综合大香蕉| 欧美3d第一页| 精华霜和精华液先用哪个| 国产69精品久久久久777片| 成人三级黄色视频| 国产成人a区在线观看| 免费大片18禁| 在线天堂最新版资源| 国产视频内射| 国产私拍福利视频在线观看| 免费av观看视频| videossex国产| 波多野结衣高清无吗| 亚洲第一电影网av| or卡值多少钱| 国产熟女欧美一区二区| 最新在线观看一区二区三区| 久久久久久久久久黄片| 国产免费男女视频| 男人和女人高潮做爰伦理| 亚洲国产精品sss在线观看| 亚洲av成人精品一区久久| 少妇裸体淫交视频免费看高清| 日本黄色视频三级网站网址| 国产亚洲欧美98| 99热只有精品国产| 两性午夜刺激爽爽歪歪视频在线观看| 成年人黄色毛片网站| 一级a爱片免费观看的视频| 国产色婷婷99| 狠狠狠狠99中文字幕| 成人av一区二区三区在线看| 免费黄网站久久成人精品| 国产国拍精品亚洲av在线观看| 久久久久国产精品人妻aⅴ院| 国产真实乱freesex| 午夜a级毛片| 午夜视频国产福利| 97热精品久久久久久| 国产男靠女视频免费网站| 真实男女啪啪啪动态图| 美女xxoo啪啪120秒动态图| 亚洲国产高清在线一区二区三| 又爽又黄无遮挡网站| 久久人人爽人人爽人人片va| 两性午夜刺激爽爽歪歪视频在线观看| 久久久国产成人免费| 国内揄拍国产精品人妻在线| 午夜日韩欧美国产| 国模一区二区三区四区视频| 欧美日本亚洲视频在线播放| 亚洲在线观看片| 岛国在线免费视频观看| 日韩一本色道免费dvd| 久久精品国产鲁丝片午夜精品 | av天堂中文字幕网| 99热6这里只有精品| 色吧在线观看| 不卡视频在线观看欧美| 久99久视频精品免费| 色综合亚洲欧美另类图片| 亚洲性久久影院| 午夜福利视频1000在线观看| 亚洲精品成人久久久久久| 国产激情偷乱视频一区二区| 亚洲欧美日韩卡通动漫| 又粗又爽又猛毛片免费看| 一本精品99久久精品77| 超碰av人人做人人爽久久| 国产白丝娇喘喷水9色精品| 性色avwww在线观看| 久久久久久九九精品二区国产| 午夜精品久久久久久毛片777| 99久久久亚洲精品蜜臀av| 非洲黑人性xxxx精品又粗又长| 无遮挡黄片免费观看| 亚洲美女搞黄在线观看 | 特大巨黑吊av在线直播| 久久人人精品亚洲av| 成熟少妇高潮喷水视频| 久久精品国产清高在天天线| 国产伦精品一区二区三区四那| 舔av片在线| 欧美日韩国产亚洲二区| 国产精品久久久久久av不卡| 亚洲一区二区三区色噜噜| xxxwww97欧美| 国产精品av视频在线免费观看| 精品一区二区三区av网在线观看| 久久人人爽人人爽人人片va| 韩国av一区二区三区四区| 欧美黑人巨大hd| 国产一区二区三区视频了| 两人在一起打扑克的视频| 亚洲美女视频黄频| 免费人成在线观看视频色| 欧美日韩乱码在线| 亚洲人成网站在线播放欧美日韩| 日韩 亚洲 欧美在线| www日本黄色视频网|