• <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)和超曲面
    亚洲精品一二三| 国产成人精品福利久久| 春色校园在线视频观看| 免费看日本二区| 内地一区二区视频在线| 97超碰精品成人国产| 欧美区成人在线视频| 国产乱人视频| 观看av在线不卡| 日韩中文字幕视频在线看片 | 99久久精品国产国产毛片| 日韩欧美精品免费久久| 麻豆乱淫一区二区| 国产成人a区在线观看| 十八禁网站网址无遮挡 | 久久久久国产网址| 免费大片黄手机在线观看| 日韩在线高清观看一区二区三区| 国产大屁股一区二区在线视频| 欧美极品一区二区三区四区| 熟女av电影| 全区人妻精品视频| 搡女人真爽免费视频火全软件| 色哟哟·www| av在线播放精品| 人妻少妇偷人精品九色| 国产成人aa在线观看| 丝瓜视频免费看黄片| 亚洲av中文字字幕乱码综合| 久久人妻熟女aⅴ| 婷婷色综合大香蕉| 中国国产av一级| 毛片女人毛片| 国产大屁股一区二区在线视频| 国产成人一区二区在线| 在线观看av片永久免费下载| 五月伊人婷婷丁香| 在线观看国产h片| 国产精品嫩草影院av在线观看| 青春草国产在线视频| 国产色爽女视频免费观看| 国产色婷婷99| 在线天堂最新版资源| 亚洲精品乱码久久久久久按摩| 99热网站在线观看| 人人妻人人添人人爽欧美一区卜 | 欧美变态另类bdsm刘玥| 色网站视频免费| 少妇人妻一区二区三区视频| 精品久久国产蜜桃| 国内精品宾馆在线| 国产乱人偷精品视频| 在线播放无遮挡| 亚洲人成网站高清观看| 国产一区二区三区综合在线观看 | 少妇人妻 视频| 伦理电影免费视频| 天天躁日日操中文字幕| 欧美日韩亚洲高清精品| 人人妻人人爽人人添夜夜欢视频 | 91精品国产国语对白视频| 日韩一区二区三区影片| 伦精品一区二区三区| 自拍偷自拍亚洲精品老妇| 青春草亚洲视频在线观看| 久久精品国产亚洲网站| 不卡视频在线观看欧美| 国内揄拍国产精品人妻在线| 97热精品久久久久久| 直男gayav资源| 夫妻性生交免费视频一级片| 大码成人一级视频| 久久精品夜色国产| 内射极品少妇av片p| 亚洲欧美日韩无卡精品| 午夜福利视频精品| 国产 一区 欧美 日韩| 在线天堂最新版资源| 国产欧美另类精品又又久久亚洲欧美| 日韩 亚洲 欧美在线| 久热这里只有精品99| 日韩av在线免费看完整版不卡| 欧美另类一区| 晚上一个人看的免费电影| 黄色配什么色好看| 亚洲综合色惰| 精品一区二区三区视频在线| 日韩成人伦理影院| 最新中文字幕久久久久| 最近手机中文字幕大全| 亚洲av中文字字幕乱码综合| 精华霜和精华液先用哪个| 一区二区三区四区激情视频| 精品亚洲成国产av| 久久热精品热| 观看美女的网站| 晚上一个人看的免费电影| 性高湖久久久久久久久免费观看| 免费不卡的大黄色大毛片视频在线观看| 婷婷色综合www| 夜夜爽夜夜爽视频| 成人毛片60女人毛片免费| 有码 亚洲区| 免费黄色在线免费观看| 人妻制服诱惑在线中文字幕| 五月玫瑰六月丁香| 免费在线观看成人毛片| 亚洲精华国产精华液的使用体验| 99热国产这里只有精品6| 国产成人a∨麻豆精品| 高清毛片免费看| 全区人妻精品视频| 国产精品爽爽va在线观看网站| 久久精品久久久久久久性| 内射极品少妇av片p| 亚洲欧洲国产日韩| 在线 av 中文字幕| 一个人看视频在线观看www免费| 少妇被粗大猛烈的视频| 亚洲av免费高清在线观看| 成年女人在线观看亚洲视频| 精品久久久久久久久亚洲| 国产精品久久久久久av不卡| 国产伦精品一区二区三区四那| 久久99精品国语久久久| 久久午夜福利片| 国产精品麻豆人妻色哟哟久久| 边亲边吃奶的免费视频| 综合色丁香网| av线在线观看网站| 男的添女的下面高潮视频| 国产高清国产精品国产三级 | 日日撸夜夜添| av女优亚洲男人天堂| av视频免费观看在线观看| 久久久久久九九精品二区国产| 丰满人妻一区二区三区视频av| 九色成人免费人妻av| 国产男女内射视频| 久久 成人 亚洲| 国产精品99久久99久久久不卡 | 国产精品成人在线| 中文精品一卡2卡3卡4更新| 观看美女的网站| 久久精品国产亚洲网站| 80岁老熟妇乱子伦牲交| 久久精品国产自在天天线| videossex国产| 一级爰片在线观看| 国产黄频视频在线观看| videossex国产| 成人一区二区视频在线观看| 少妇被粗大猛烈的视频| 人妻少妇偷人精品九色| 亚洲国产日韩一区二区| 欧美精品国产亚洲| 亚洲av不卡在线观看| 国产av码专区亚洲av| 精品久久国产蜜桃| 蜜桃久久精品国产亚洲av| 免费观看av网站的网址| 成人影院久久| 免费黄频网站在线观看国产| 国产精品成人在线| 欧美成人a在线观看| 成人美女网站在线观看视频| 国产 一区精品| 亚洲国产最新在线播放| 国产精品人妻久久久久久| 一区在线观看完整版| 日韩欧美精品免费久久| 亚洲无线观看免费| 中国三级夫妇交换| 国产无遮挡羞羞视频在线观看| 久久久久久久久久久免费av| 高清av免费在线| 香蕉精品网在线| 国产欧美亚洲国产| 国产精品免费大片| 99热这里只有是精品在线观看| 亚洲中文av在线| 爱豆传媒免费全集在线观看| 干丝袜人妻中文字幕| 看十八女毛片水多多多| 黑人高潮一二区| 99热全是精品| 高清午夜精品一区二区三区| 在线观看三级黄色| 欧美极品一区二区三区四区| 日日啪夜夜撸| 高清欧美精品videossex| 中国美白少妇内射xxxbb| 日本vs欧美在线观看视频 | 亚洲人成网站在线观看播放| 亚洲天堂av无毛| 嫩草影院新地址| 久久国内精品自在自线图片| 国产亚洲最大av| 国产乱来视频区| 亚洲,欧美,日韩| 精华霜和精华液先用哪个| 好男人视频免费观看在线| av线在线观看网站| 久久女婷五月综合色啪小说| 高清日韩中文字幕在线| 国产免费一区二区三区四区乱码| 丝袜喷水一区| 91久久精品电影网| 久久ye,这里只有精品| 成年免费大片在线观看| 岛国毛片在线播放| 十分钟在线观看高清视频www | 精品久久久久久久久av| 久久99热6这里只有精品| 麻豆成人av视频| 亚洲欧美成人综合另类久久久| 偷拍熟女少妇极品色| 久久久久性生活片| 国产成人午夜福利电影在线观看| 午夜日本视频在线| 午夜免费鲁丝| 女的被弄到高潮叫床怎么办| 女人久久www免费人成看片| av国产免费在线观看| 国产成人精品婷婷| 成年人午夜在线观看视频| 成人美女网站在线观看视频| av在线蜜桃| 日韩av免费高清视频| 在线免费十八禁| 成年女人在线观看亚洲视频| 国产高潮美女av| 欧美另类一区| 99热这里只有是精品在线观看| 国产精品99久久久久久久久| 3wmmmm亚洲av在线观看| 欧美zozozo另类| 亚洲精品乱码久久久v下载方式| av福利片在线观看| 国产亚洲91精品色在线| 亚洲欧美一区二区三区国产| 色5月婷婷丁香| 日本色播在线视频| 国产视频首页在线观看| 亚洲一级一片aⅴ在线观看| 夜夜骑夜夜射夜夜干| 中文在线观看免费www的网站| av天堂中文字幕网| 七月丁香在线播放| 亚洲四区av| 国产成人精品一,二区| 制服丝袜香蕉在线| 国产男人的电影天堂91| 国产精品久久久久久精品电影小说 | 在线免费十八禁| 尤物成人国产欧美一区二区三区| .国产精品久久| 国产成人免费观看mmmm| 国产在视频线精品| 免费人成在线观看视频色| 三级国产精品欧美在线观看| 国产亚洲91精品色在线| 黄色视频在线播放观看不卡| 人妻 亚洲 视频| 亚洲一区二区三区欧美精品| 一级毛片电影观看| 久久99热6这里只有精品| 亚洲欧美成人精品一区二区| 街头女战士在线观看网站| 欧美zozozo另类| 丝袜喷水一区| 亚洲精品亚洲一区二区| 欧美日韩视频高清一区二区三区二| 久久国内精品自在自线图片| 国产精品秋霞免费鲁丝片| 最近中文字幕高清免费大全6| 国产69精品久久久久777片| 国产在线免费精品| 国产成人免费观看mmmm| 国产黄频视频在线观看| 91在线精品国自产拍蜜月| 丝瓜视频免费看黄片| 精品国产露脸久久av麻豆| 大片电影免费在线观看免费| 三级国产精品片| 伦理电影大哥的女人| 一区在线观看完整版| 18禁在线无遮挡免费观看视频| 免费看不卡的av| 久久久国产一区二区| 亚洲av在线观看美女高潮| 国产精品成人在线| 色视频www国产| 日韩人妻高清精品专区| 中文乱码字字幕精品一区二区三区| av天堂中文字幕网| 好男人视频免费观看在线| 日韩,欧美,国产一区二区三区| 黑人猛操日本美女一级片| 一级爰片在线观看| 美女国产视频在线观看| 精品亚洲乱码少妇综合久久| 视频区图区小说| 亚洲av综合色区一区| 日韩人妻高清精品专区| 国产黄色视频一区二区在线观看| av网站免费在线观看视频| 亚洲va在线va天堂va国产| 网址你懂的国产日韩在线| 偷拍熟女少妇极品色| 亚洲人成网站在线观看播放| 91久久精品国产一区二区三区| 成人美女网站在线观看视频| 亚洲av福利一区| 老女人水多毛片| 精品一区二区免费观看| 久久精品国产亚洲av天美| 欧美日韩视频精品一区| 啦啦啦啦在线视频资源| 街头女战士在线观看网站| 一级黄片播放器| 国产精品一区二区性色av| 亚洲欧美清纯卡通| 久久人妻熟女aⅴ| 日韩欧美精品免费久久| 亚洲成色77777| 免费观看无遮挡的男女| 国产白丝娇喘喷水9色精品| 国产高清有码在线观看视频| 黄色一级大片看看| 国产成人a∨麻豆精品| 成人高潮视频无遮挡免费网站| 一级毛片久久久久久久久女| 国产成人一区二区在线| 日本av手机在线免费观看| 欧美老熟妇乱子伦牲交| 超碰av人人做人人爽久久| 日韩成人伦理影院| 综合色丁香网| 黄色配什么色好看| av国产精品久久久久影院| 亚洲无线观看免费| 日本av手机在线免费观看| 亚洲,欧美,日韩| 偷拍熟女少妇极品色| 在线播放无遮挡| 狂野欧美激情性xxxx在线观看| av专区在线播放| 有码 亚洲区| 在线观看免费视频网站a站| 国产国拍精品亚洲av在线观看| 国产免费又黄又爽又色| 一级片'在线观看视频| 亚洲国产毛片av蜜桃av| 一级片'在线观看视频| 久久国产精品男人的天堂亚洲 | 涩涩av久久男人的天堂| 黑人高潮一二区| 亚洲国产毛片av蜜桃av| 又粗又硬又长又爽又黄的视频| 国产高清三级在线| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲精品aⅴ在线观看| 亚洲av免费高清在线观看| 免费观看在线日韩| 在现免费观看毛片| 天天躁夜夜躁狠狠久久av| 国产成人aa在线观看| 久久久亚洲精品成人影院| 简卡轻食公司| 成人毛片a级毛片在线播放| 人妻少妇偷人精品九色| 久久久午夜欧美精品| 久久久久国产精品人妻一区二区| 亚洲精品久久久久久婷婷小说| 51国产日韩欧美| 午夜日本视频在线| 国产无遮挡羞羞视频在线观看| 深夜a级毛片| 国产精品久久久久久精品古装| av播播在线观看一区| 成人二区视频| 免费播放大片免费观看视频在线观看| 精品亚洲乱码少妇综合久久| av又黄又爽大尺度在线免费看| 在线观看一区二区三区激情| 国产亚洲欧美精品永久| 18+在线观看网站| 国产视频首页在线观看| 亚洲久久久国产精品| 黄色怎么调成土黄色| 黄色配什么色好看| av卡一久久| 人体艺术视频欧美日本| 午夜老司机福利剧场| 熟女电影av网| 男人和女人高潮做爰伦理| 国产精品一区二区性色av| 国产 精品1| 亚洲av欧美aⅴ国产| 国产精品一区www在线观看| 男女啪啪激烈高潮av片| 亚洲精品中文字幕在线视频 | 各种免费的搞黄视频| av在线老鸭窝| 97在线视频观看| 欧美高清成人免费视频www| 久久人人爽av亚洲精品天堂 | 在现免费观看毛片| 午夜日本视频在线| 如何舔出高潮| 亚洲国产精品专区欧美| 精品一区二区三卡| 伦理电影免费视频| 久久鲁丝午夜福利片| 天天躁日日操中文字幕| 亚洲欧美精品自产自拍| 日韩免费高清中文字幕av| 日韩,欧美,国产一区二区三区| 午夜日本视频在线| 精华霜和精华液先用哪个| 欧美成人a在线观看| 精品一品国产午夜福利视频| 婷婷色麻豆天堂久久| 国产精品久久久久久av不卡| 午夜免费男女啪啪视频观看| 黑人高潮一二区| 久久精品国产亚洲av天美| 国产精品无大码| 国产黄色视频一区二区在线观看| 一级毛片黄色毛片免费观看视频| 国产免费视频播放在线视频| 国产成人一区二区在线| 一级a做视频免费观看| 搡老乐熟女国产| 国产美女午夜福利| 免费观看av网站的网址| 啦啦啦视频在线资源免费观看| h视频一区二区三区| 精品国产三级普通话版| 亚洲国产精品成人久久小说| 看十八女毛片水多多多| 18禁裸乳无遮挡动漫免费视频| 国产无遮挡羞羞视频在线观看| 日韩av不卡免费在线播放| 日韩免费高清中文字幕av| 亚洲色图av天堂| 成人18禁高潮啪啪吃奶动态图 | 男人舔奶头视频| 久久久久网色| 人体艺术视频欧美日本| 国产视频首页在线观看| 日韩欧美一区视频在线观看 | 丰满迷人的少妇在线观看| 久久精品国产亚洲av天美| 亚洲真实伦在线观看| 水蜜桃什么品种好| 97在线人人人人妻| 国产伦在线观看视频一区| 五月玫瑰六月丁香| 国产精品久久久久久久久免| 成人毛片a级毛片在线播放| 一本久久精品| 亚洲色图av天堂| 免费观看性生交大片5| 亚洲在久久综合| 美女xxoo啪啪120秒动态图| 18+在线观看网站| 夫妻午夜视频| 欧美一级a爱片免费观看看| 日韩大片免费观看网站| 久久精品国产鲁丝片午夜精品| 免费播放大片免费观看视频在线观看| 久热久热在线精品观看| 亚洲美女黄色视频免费看| a级毛片免费高清观看在线播放| 午夜免费观看性视频| 超碰av人人做人人爽久久| 波野结衣二区三区在线| 国产成人aa在线观看| av免费观看日本| 观看免费一级毛片| 伦理电影免费视频| 国产黄片美女视频| 精品一区二区三区视频在线| 久久国产乱子免费精品| 久久久久久久亚洲中文字幕| 美女xxoo啪啪120秒动态图| 免费av不卡在线播放| 麻豆成人午夜福利视频| 五月伊人婷婷丁香| 夜夜爽夜夜爽视频| 精品少妇久久久久久888优播| 国产av精品麻豆| 久久精品夜色国产| 亚洲欧美成人综合另类久久久| 日韩视频在线欧美| 嫩草影院新地址| 欧美三级亚洲精品| 妹子高潮喷水视频| 爱豆传媒免费全集在线观看| 国产成人一区二区在线| 婷婷色综合www| 亚洲真实伦在线观看| 亚洲国产精品国产精品| 亚洲欧洲日产国产| 免费大片黄手机在线观看| 亚洲av免费高清在线观看| av播播在线观看一区| 乱码一卡2卡4卡精品| 男女国产视频网站| 久久女婷五月综合色啪小说| 最近的中文字幕免费完整| 欧美bdsm另类| 精品一区二区三卡| 久久6这里有精品| 亚洲av日韩在线播放| 男人狂女人下面高潮的视频| 五月玫瑰六月丁香| 免费看光身美女| 在现免费观看毛片| 美女主播在线视频| 欧美一级a爱片免费观看看| 日韩成人av中文字幕在线观看| 亚洲电影在线观看av| 久久精品久久久久久久性| 中国三级夫妇交换| 人妻一区二区av| 最近手机中文字幕大全| 欧美精品一区二区免费开放| 亚洲av成人精品一区久久| 欧美日韩视频高清一区二区三区二| 全区人妻精品视频| 亚洲av综合色区一区| 国产欧美另类精品又又久久亚洲欧美| 国产精品99久久久久久久久| 99久国产av精品国产电影| 亚洲熟女精品中文字幕| 中文字幕久久专区| av国产精品久久久久影院| 国产高潮美女av| 国产爽快片一区二区三区| 夫妻午夜视频| 老司机影院毛片| 韩国av在线不卡| 99久久中文字幕三级久久日本| 久久av网站| 日韩欧美精品免费久久| 久久99蜜桃精品久久| 观看免费一级毛片| 日本av手机在线免费观看| 国产免费福利视频在线观看| 欧美日韩亚洲高清精品| 亚洲欧洲日产国产| 黑丝袜美女国产一区| 日韩一区二区三区影片| 各种免费的搞黄视频| 伊人久久精品亚洲午夜| 简卡轻食公司| 九九久久精品国产亚洲av麻豆| 午夜视频国产福利| 亚洲高清免费不卡视频| 蜜臀久久99精品久久宅男| 精品人妻一区二区三区麻豆| 一级毛片久久久久久久久女| av在线蜜桃| 成年av动漫网址| 国产视频首页在线观看| 全区人妻精品视频| 晚上一个人看的免费电影| 男女边吃奶边做爰视频| 成人亚洲精品一区在线观看 | 国产伦理片在线播放av一区| 亚洲精品国产av成人精品| 亚洲欧美精品专区久久| 亚洲国产毛片av蜜桃av| 国精品久久久久久国模美| 久久精品人妻少妇| 国产大屁股一区二区在线视频| 色视频在线一区二区三区| 欧美精品一区二区大全| 国产乱人偷精品视频| 欧美日韩视频高清一区二区三区二| 国产免费一级a男人的天堂| 亚洲人成网站在线播| 一级毛片电影观看| 夜夜看夜夜爽夜夜摸| 男人狂女人下面高潮的视频| 日本猛色少妇xxxxx猛交久久| 三级经典国产精品| 99热6这里只有精品| 中国三级夫妇交换| 九草在线视频观看| 99热全是精品| 亚洲国产欧美在线一区| 国产黄片美女视频| 啦啦啦啦在线视频资源| 九色成人免费人妻av| 黑丝袜美女国产一区| 亚洲精品456在线播放app| 欧美国产精品一级二级三级 | 亚洲精品,欧美精品| 亚洲精品乱码久久久v下载方式| 91久久精品国产一区二区三区| 亚洲国产高清在线一区二区三| 久久人人爽人人片av| 搡女人真爽免费视频火全软件| 联通29元200g的流量卡| 在线观看免费高清a一片| 日本与韩国留学比较| 联通29元200g的流量卡| 丰满少妇做爰视频| 精品久久久久久久久亚洲| kizo精华| 在线免费十八禁|