• <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)和超曲面
    狠狠狠狠99中文字幕| 色综合欧美亚洲国产小说| 精品午夜福利视频在线观看一区| 欧美绝顶高潮抽搐喷水| 无遮挡黄片免费观看| 一二三四社区在线视频社区8| 亚洲五月天丁香| 桃色一区二区三区在线观看| 日韩欧美在线乱码| 亚洲精品在线观看二区| 长腿黑丝高跟| 国产欧美日韩一区二区精品| 看免费av毛片| 一区二区三区高清视频在线| 国产午夜精品久久久久久| 一本久久中文字幕| 久久久久久久久免费视频了| 亚洲精品在线美女| 在线免费观看的www视频| 男女视频在线观看网站免费| 特大巨黑吊av在线直播| 欧美精品啪啪一区二区三区| avwww免费| 免费看美女性在线毛片视频| 又爽又黄无遮挡网站| 国产成人精品久久二区二区免费| 精品人妻1区二区| 亚洲精品国产精品久久久不卡| 午夜免费激情av| 老司机午夜福利在线观看视频| 亚洲国产欧美网| 一区二区三区高清视频在线| 欧美性猛交黑人性爽| 好男人在线观看高清免费视频| 伦理电影免费视频| 黄色日韩在线| 午夜a级毛片| 91在线精品国自产拍蜜月 | 老司机午夜福利在线观看视频| netflix在线观看网站| 12—13女人毛片做爰片一| 女人高潮潮喷娇喘18禁视频| 国内精品美女久久久久久| 国产精品亚洲一级av第二区| 久久久久久久久中文| 99热这里只有是精品50| 日本在线视频免费播放| 制服人妻中文乱码| 国产 一区 欧美 日韩| 级片在线观看| av黄色大香蕉| 国产欧美日韩精品一区二区| 一个人看视频在线观看www免费 | 久久久久精品国产欧美久久久| 亚洲avbb在线观看| 国产97色在线日韩免费| 亚洲第一欧美日韩一区二区三区| 欧美国产日韩亚洲一区| 法律面前人人平等表现在哪些方面| 成人18禁在线播放| 午夜精品一区二区三区免费看| 国产伦在线观看视频一区| 日韩中文字幕欧美一区二区| 不卡一级毛片| 日本熟妇午夜| 成人一区二区视频在线观看| 欧美日韩亚洲国产一区二区在线观看| 两个人看的免费小视频| 久久久久久久久中文| 亚洲精品一卡2卡三卡4卡5卡| av国产免费在线观看| 一夜夜www| 婷婷丁香在线五月| 国产视频一区二区在线看| 国产精品久久久人人做人人爽| 欧美xxxx黑人xx丫x性爽| 色综合欧美亚洲国产小说| 国产精品一区二区三区四区免费观看 | 动漫黄色视频在线观看| 亚洲中文av在线| 一进一出抽搐动态| 亚洲国产精品合色在线| 国产91精品成人一区二区三区| www.www免费av| 日本一本二区三区精品| 老熟妇仑乱视频hdxx| 九色成人免费人妻av| 精品熟女少妇八av免费久了| 无限看片的www在线观看| 久久国产乱子伦精品免费另类| 免费大片18禁| 亚洲男人的天堂狠狠| 欧美日韩精品网址| 日韩三级视频一区二区三区| 国产69精品久久久久777片 | 成人av在线播放网站| 免费一级毛片在线播放高清视频| 无遮挡黄片免费观看| 欧美成人性av电影在线观看| 制服人妻中文乱码| 日韩欧美国产一区二区入口| 91久久精品国产一区二区成人 | 免费观看精品视频网站| 亚洲真实伦在线观看| 97超视频在线观看视频| 久久精品国产清高在天天线| 久久中文看片网| 男女床上黄色一级片免费看| 国产成人精品久久二区二区91| 午夜福利免费观看在线| 亚洲七黄色美女视频| 午夜福利在线观看吧| 欧美日韩福利视频一区二区| 成人亚洲精品av一区二区| 精品久久久久久久人妻蜜臀av| 亚洲av成人精品一区久久| 免费观看精品视频网站| av天堂中文字幕网| 精品国产超薄肉色丝袜足j| 国产男靠女视频免费网站| 久久99热这里只有精品18| 国产精品爽爽va在线观看网站| 欧美在线一区亚洲| 美女被艹到高潮喷水动态| 久久欧美精品欧美久久欧美| 一个人免费在线观看电影 | 亚洲av成人av| 18禁黄网站禁片午夜丰满| 亚洲精品乱码久久久v下载方式 | 国产不卡一卡二| 亚洲五月天丁香| 亚洲第一欧美日韩一区二区三区| 国产成+人综合+亚洲专区| 人人妻人人澡欧美一区二区| 99在线人妻在线中文字幕| 精品欧美国产一区二区三| 中文字幕av在线有码专区| 天天躁狠狠躁夜夜躁狠狠躁| 国产精华一区二区三区| 少妇的逼水好多| 91av网一区二区| 国产精品 国内视频| 亚洲在线观看片| 午夜福利高清视频| 成人午夜高清在线视频| 国产精品电影一区二区三区| 少妇丰满av| 日韩高清综合在线| 免费一级毛片在线播放高清视频| 亚洲自拍偷在线| 69av精品久久久久久| 天堂影院成人在线观看| 人人妻,人人澡人人爽秒播| 狠狠狠狠99中文字幕| 久久午夜综合久久蜜桃| 精品电影一区二区在线| 母亲3免费完整高清在线观看| www.熟女人妻精品国产| 日本与韩国留学比较| 欧美激情在线99| 亚洲天堂国产精品一区在线| 日韩欧美 国产精品| 色综合站精品国产| 亚洲第一电影网av| 欧美又色又爽又黄视频| 国产成人影院久久av| 国产欧美日韩精品一区二区| 国产一区二区在线观看日韩 | 97超视频在线观看视频| 91av网站免费观看| 国产精品 国内视频| 中文字幕av在线有码专区| 免费在线观看影片大全网站| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品粉嫩美女一区| 色播亚洲综合网| 亚洲欧美激情综合另类| 国产高清三级在线| 熟女电影av网| 天堂动漫精品| 男人的好看免费观看在线视频| 午夜激情福利司机影院| 校园春色视频在线观看| 搡老熟女国产l中国老女人| 婷婷亚洲欧美| 老熟妇仑乱视频hdxx| 又粗又爽又猛毛片免费看| av国产免费在线观看| 欧美性猛交╳xxx乱大交人| 后天国语完整版免费观看| 成人三级做爰电影| 看黄色毛片网站| 麻豆久久精品国产亚洲av| 97超视频在线观看视频| 亚洲成av人片免费观看| 久久久久久久久中文| 亚洲片人在线观看| 很黄的视频免费| 成人国产综合亚洲| 国产成人福利小说| 亚洲天堂国产精品一区在线| 神马国产精品三级电影在线观看| 一夜夜www| 91老司机精品| 亚洲成人中文字幕在线播放| 国产亚洲av嫩草精品影院| 可以在线观看毛片的网站| 日韩欧美精品v在线| 欧美xxxx黑人xx丫x性爽| 一区二区三区国产精品乱码| 少妇人妻一区二区三区视频| 国产亚洲av嫩草精品影院| 国产伦一二天堂av在线观看| 国产精品久久久人人做人人爽| 99久久国产精品久久久| 老司机在亚洲福利影院| 国产三级黄色录像| 日韩欧美国产在线观看| 2021天堂中文幕一二区在线观| 日韩高清综合在线| 国内精品一区二区在线观看| 日本精品一区二区三区蜜桃| 99久久精品国产亚洲精品| 变态另类成人亚洲欧美熟女| 久久久久久久午夜电影| 99热这里只有是精品50| 久久精品国产亚洲av香蕉五月| 成人av在线播放网站| 婷婷亚洲欧美| 国产精品av久久久久免费| 制服人妻中文乱码| 国内揄拍国产精品人妻在线| 国产成人av教育| 国产乱人伦免费视频| 露出奶头的视频| 日韩欧美国产在线观看| 一a级毛片在线观看| 精品欧美国产一区二区三| 欧美不卡视频在线免费观看| 麻豆国产97在线/欧美| 免费看日本二区| 精品国产超薄肉色丝袜足j| 精品免费久久久久久久清纯| 欧美最黄视频在线播放免费| 日韩 欧美 亚洲 中文字幕| 熟女人妻精品中文字幕| 三级男女做爰猛烈吃奶摸视频| 午夜福利免费观看在线| 中文在线观看免费www的网站| 精品人妻1区二区| 中国美女看黄片| 国产亚洲精品av在线| 男女午夜视频在线观看| 久久中文看片网| 母亲3免费完整高清在线观看| 美女扒开内裤让男人捅视频| 欧美乱色亚洲激情| 一a级毛片在线观看| 久久中文字幕人妻熟女| 欧美成狂野欧美在线观看| 国产精品 欧美亚洲| 亚洲欧美精品综合一区二区三区| 99国产极品粉嫩在线观看| 成在线人永久免费视频| 少妇的丰满在线观看| 色老头精品视频在线观看| 国产高清激情床上av| 欧美黑人巨大hd| 久久午夜综合久久蜜桃| 12—13女人毛片做爰片一| 亚洲av电影在线进入| 欧美丝袜亚洲另类 | 免费在线观看视频国产中文字幕亚洲| 国产一区在线观看成人免费| 麻豆av在线久日| 夜夜躁狠狠躁天天躁| 人人妻,人人澡人人爽秒播| 国产精品女同一区二区软件 | 国产蜜桃级精品一区二区三区| 这个男人来自地球电影免费观看| 国内精品久久久久精免费| 亚洲真实伦在线观看| 99riav亚洲国产免费| 国产av麻豆久久久久久久| 好男人电影高清在线观看| 一个人免费在线观看电影 | 日韩精品青青久久久久久| 欧美激情在线99| 国产精品野战在线观看| 久久精品国产亚洲av香蕉五月| 国产成人aa在线观看| 国产精品亚洲美女久久久| 国产欧美日韩一区二区三| 亚洲国产色片| 亚洲va日本ⅴa欧美va伊人久久| 久久精品aⅴ一区二区三区四区| cao死你这个sao货| 午夜福利成人在线免费观看| 不卡av一区二区三区| 国产高清激情床上av| 天天添夜夜摸| 欧美3d第一页| 亚洲国产精品合色在线| 国产主播在线观看一区二区| 久久天躁狠狠躁夜夜2o2o| 欧洲精品卡2卡3卡4卡5卡区| 亚洲国产高清在线一区二区三| 久久久水蜜桃国产精品网| 美女 人体艺术 gogo| 午夜福利在线观看免费完整高清在 | 九九在线视频观看精品| 国模一区二区三区四区视频 | 一本久久中文字幕| 搡老熟女国产l中国老女人| 97人妻精品一区二区三区麻豆| 国产伦精品一区二区三区四那| 亚洲最大成人中文| 久久久久久人人人人人| 淫秽高清视频在线观看| 精品久久蜜臀av无| 女同久久另类99精品国产91| 琪琪午夜伦伦电影理论片6080| www国产在线视频色| 麻豆成人av在线观看| 香蕉丝袜av| 国产在线精品亚洲第一网站| 岛国在线观看网站| 午夜a级毛片| 午夜福利在线观看吧| 国产人伦9x9x在线观看| 又粗又爽又猛毛片免费看| 免费观看人在逋| 久久精品国产99精品国产亚洲性色| 国产精品1区2区在线观看.| 国产免费av片在线观看野外av| tocl精华| 一本综合久久免费| 久久精品国产清高在天天线| 无人区码免费观看不卡| 国语自产精品视频在线第100页| 香蕉av资源在线| 国产亚洲精品久久久com| 精品99又大又爽又粗少妇毛片 | 免费av不卡在线播放| 亚洲中文av在线| www国产在线视频色| 日本黄大片高清| 日本a在线网址| 看免费av毛片| 亚洲国产欧美人成| 人人妻人人看人人澡| 欧美黄色片欧美黄色片| 啦啦啦观看免费观看视频高清| 亚洲精品色激情综合| 国产高清videossex| 久久精品国产亚洲av香蕉五月| 亚洲男人的天堂狠狠| 男人舔女人下体高潮全视频| 色av中文字幕| 最新中文字幕久久久久 | 久久欧美精品欧美久久欧美| 99久久国产精品久久久| 亚洲精品一区av在线观看| 国产日本99.免费观看| 神马国产精品三级电影在线观看| 黄色片一级片一级黄色片| 男人舔女人的私密视频| 99久久精品国产亚洲精品| 亚洲精品一区av在线观看| 桃色一区二区三区在线观看| 亚洲精品一区av在线观看| 国产视频内射| 久久国产乱子伦精品免费另类| 女人被狂操c到高潮| 精品久久久久久,| e午夜精品久久久久久久| 国产三级中文精品| 色视频www国产| 国产男靠女视频免费网站| 两人在一起打扑克的视频| 欧美又色又爽又黄视频| 欧美3d第一页| avwww免费| 亚洲精品乱码久久久v下载方式 | 91久久精品国产一区二区成人 | 成人av一区二区三区在线看| 午夜精品久久久久久毛片777| 三级国产精品欧美在线观看 | 国产精品久久久av美女十八| 婷婷丁香在线五月| 久久久成人免费电影| 黄色丝袜av网址大全| 他把我摸到了高潮在线观看| 最新美女视频免费是黄的| 色吧在线观看| 日日夜夜操网爽| 久久久久久久午夜电影| 欧美日韩国产亚洲二区| 亚洲欧美一区二区三区黑人| 亚洲国产精品合色在线| 亚洲精品久久国产高清桃花| 中文亚洲av片在线观看爽| 久久久久九九精品影院| 一区二区三区高清视频在线| 国产真实乱freesex| 日韩免费av在线播放| 三级毛片av免费| 91字幕亚洲| 十八禁人妻一区二区| 好看av亚洲va欧美ⅴa在| 深夜精品福利| 桃色一区二区三区在线观看| 在线a可以看的网站| 久久久久久久午夜电影| 成人鲁丝片一二三区免费| 国产精品1区2区在线观看.| 19禁男女啪啪无遮挡网站| 老司机福利观看| 亚洲色图av天堂| 成在线人永久免费视频| 啦啦啦韩国在线观看视频| 国产精品av视频在线免费观看| 亚洲人成伊人成综合网2020| 国产精品久久久久久久电影 | 亚洲午夜理论影院| av天堂在线播放| 免费人成视频x8x8入口观看| 一级毛片高清免费大全| 巨乳人妻的诱惑在线观看| 97超级碰碰碰精品色视频在线观看| 国产一区二区三区视频了| 欧美性猛交黑人性爽| 最近最新中文字幕大全免费视频| 怎么达到女性高潮| 一进一出好大好爽视频| 操出白浆在线播放| 国产爱豆传媒在线观看| 精品久久久久久久末码| 国产精品一区二区免费欧美| 99久久精品热视频| 一卡2卡三卡四卡精品乱码亚洲| 午夜福利18| 青草久久国产| 亚洲欧美日韩东京热| 午夜亚洲福利在线播放| 99久久精品一区二区三区| 午夜a级毛片| 色尼玛亚洲综合影院| 欧美午夜高清在线| 精品久久久久久久久久久久久| 精品国产三级普通话版| 成人亚洲精品av一区二区| 欧美在线黄色| 可以在线观看的亚洲视频| 亚洲成人久久性| 免费在线观看日本一区| 国内精品一区二区在线观看| 国产一区二区三区在线臀色熟女| 亚洲av中文字字幕乱码综合| 一二三四社区在线视频社区8| 亚洲人成网站高清观看| 久久中文看片网| 国产高清videossex| 日本与韩国留学比较| 非洲黑人性xxxx精品又粗又长| 男人舔女人的私密视频| 狂野欧美白嫩少妇大欣赏| 国产高潮美女av| 男人和女人高潮做爰伦理| 不卡一级毛片| 成年人黄色毛片网站| av天堂中文字幕网| 色精品久久人妻99蜜桃| 欧美绝顶高潮抽搐喷水| 国产精品爽爽va在线观看网站| 日韩大尺度精品在线看网址| 国产精品久久久av美女十八| 亚洲精品456在线播放app | 亚洲欧美日韩高清在线视频| 欧美乱码精品一区二区三区| 国产精品久久久久久精品电影| 成人鲁丝片一二三区免费| 此物有八面人人有两片| 一本一本综合久久| 成人特级黄色片久久久久久久| 久久精品91无色码中文字幕| 无限看片的www在线观看| 日韩欧美三级三区| 在线播放国产精品三级| 亚洲国产看品久久| 久久久国产成人精品二区| 无遮挡黄片免费观看| 热99re8久久精品国产| av天堂中文字幕网| 精品免费久久久久久久清纯| 亚洲无线观看免费| 国产精品精品国产色婷婷| 国内精品久久久久精免费| 中文字幕人妻丝袜一区二区| 黑人欧美特级aaaaaa片| 一a级毛片在线观看| 熟女少妇亚洲综合色aaa.| 亚洲成人精品中文字幕电影| x7x7x7水蜜桃| 亚洲 国产 在线| 真人一进一出gif抽搐免费| 变态另类成人亚洲欧美熟女| 亚洲精华国产精华精| av天堂中文字幕网| 精品久久久久久久末码| 女同久久另类99精品国产91| 丝袜人妻中文字幕| 九九在线视频观看精品| 国产精华一区二区三区| 国产午夜福利久久久久久| 久99久视频精品免费| 美女午夜性视频免费| 午夜精品一区二区三区免费看| 免费大片18禁| 动漫黄色视频在线观看| 老司机午夜十八禁免费视频| 精品国产美女av久久久久小说| 亚洲中文日韩欧美视频| 午夜福利在线观看吧| 伊人久久大香线蕉亚洲五| 亚洲国产中文字幕在线视频| 国产亚洲欧美在线一区二区| 国产精品一区二区三区四区免费观看 | 啦啦啦免费观看视频1| 夜夜躁狠狠躁天天躁| 国产伦人伦偷精品视频| 亚洲 欧美 日韩 在线 免费| 国产视频内射| 美女黄网站色视频| 又黄又粗又硬又大视频| 亚洲色图 男人天堂 中文字幕| 在线视频色国产色| 天堂网av新在线| 最近视频中文字幕2019在线8| 久久国产精品影院| 在线免费观看不下载黄p国产 | 国产激情欧美一区二区| 一个人免费在线观看电影 | 美女高潮的动态| 美女黄网站色视频| 搞女人的毛片| 久久午夜亚洲精品久久| 村上凉子中文字幕在线| 欧美大码av| 黄色视频,在线免费观看| 一二三四在线观看免费中文在| 最近在线观看免费完整版| 99re在线观看精品视频| 一级毛片精品| 色精品久久人妻99蜜桃| 国产男靠女视频免费网站| 国产 一区 欧美 日韩| 亚洲狠狠婷婷综合久久图片| 国产精品一区二区精品视频观看| 很黄的视频免费| 久久中文看片网| 搞女人的毛片| 别揉我奶头~嗯~啊~动态视频| avwww免费| 人人妻人人看人人澡| 丝袜人妻中文字幕| 久久久国产欧美日韩av| 日本一本二区三区精品| 最近最新中文字幕大全免费视频| 日韩免费av在线播放| 成人18禁在线播放| av在线天堂中文字幕| 国产99白浆流出| 精品午夜福利视频在线观看一区| 老汉色av国产亚洲站长工具| 网址你懂的国产日韩在线| 99久久国产精品久久久| 18禁国产床啪视频网站| 国产伦一二天堂av在线观看| 国产欧美日韩一区二区精品| 国产又色又爽无遮挡免费看| 国产激情偷乱视频一区二区| 欧美黑人巨大hd| 免费看a级黄色片| 精品久久久久久久久久免费视频| av视频在线观看入口| 国产乱人视频| 制服人妻中文乱码| 在线观看美女被高潮喷水网站 | 一本久久中文字幕| 我的老师免费观看完整版| 久99久视频精品免费| 2021天堂中文幕一二区在线观| 欧美一级毛片孕妇| 在线观看一区二区三区| 99久久综合精品五月天人人| 少妇丰满av| 天堂影院成人在线观看| 99精品欧美一区二区三区四区| 国产av麻豆久久久久久久| 久久中文字幕人妻熟女| 国产精品一区二区三区四区免费观看 | 成人精品一区二区免费| 性色av乱码一区二区三区2| 日本成人三级电影网站| 久久久国产成人免费| 午夜福利在线在线| 狂野欧美激情性xxxx| 色综合亚洲欧美另类图片| 看免费av毛片| 亚洲性夜色夜夜综合| 欧美不卡视频在线免费观看| 桃红色精品国产亚洲av| 999久久久国产精品视频| 又爽又黄无遮挡网站|