亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于灰色關(guān)聯(lián)度分析的畢節(jié)市引種朝天椒研究

        2023-07-17 02:18:57謝龍安張曼王欽安朝龍任佳佳吳培云李偉
        安徽農(nóng)業(yè)科學 2023年12期
        關(guān)鍵詞:灰色關(guān)聯(lián)度朝天椒引種

        謝龍安 張曼 王欽 安朝龍 任佳佳 吳培云 李偉

        摘要 [目的]優(yōu)化畢節(jié)市朝天椒栽培品種,篩選適宜畢節(jié)市朝天椒種植品種。[方法]利用灰色系統(tǒng)理論對引種畢節(jié)種植的13個朝天椒品種果實成熟率、生育期、分枝數(shù)、主莖高、主莖粗、果實橫徑、果長、單株果實數(shù)、單株鮮果重量、含水量、干辣總產(chǎn)量、紅熟干椒產(chǎn)量、植株抗倒伏率、植株抗病率、果實可溶性糖含量、果實可溶性蛋白含量、果實VC含量、辣椒素含量、二氫辣椒素含量等19個主要指標進行關(guān)聯(lián)度分析。[結(jié)果]13個參試品種的加權(quán)關(guān)聯(lián)度值比對照“畢節(jié)朝天椒”大,占參試品種的28.5%,其中“金嬌干6號”加權(quán)關(guān)聯(lián)度值和產(chǎn)量排名第一,綜合表現(xiàn)最好;“農(nóng)望指天1號”“艷椒425” 2個品種綜合表現(xiàn)次之。[結(jié)論]“金嬌干6號”“艷椒425”“農(nóng)望指天1號”3個朝天椒品種均適合在畢節(jié)種植,“金嬌干6號”最適合在畢節(jié)市引種種植。

        關(guān)鍵詞 灰色關(guān)聯(lián)度;畢節(jié);朝天椒;引種

        中圖分類號 S 641.3? 文獻標識碼 A

        文章編號 0517-6611(2023)12-0040-05

        doi:10.3969/j.issn.0517-6611.2023.12.009

        Study on Introduction of Chinese Prickly Ash in Bijie Based on Grey Relational Grade Analysis

        XIE Long-an1,2,3,ZHANG? Man1,WANG? Qin1 et al

        (1.College of Agriculture,Guizhou University,Guiyang,Guizhou 550025;2.Bijie Vocational and Technical College,Bijie,Guizhou 551700;3.Institute of Pepper Industry,Guizhou University,Guiyang,Guizhou 550025)

        Abstract [Objective]To optimize the cultivars of Capsicum Chinense in Bijie and select the suitable cultivars for Bijie.[Method] The 19 main indexes of fruit maturity,growth period,branch number,main stem height,main stem diameter,fruit diameter,fruit length,fruit number per plant,fresh fruit weight per plant,water content,total dry pepper yield,yield of red pepper,lodging resistance,disease resistance,soluble sugar content,soluble protein content,VC content,capsaicin content and dihydrocapsaicin content in 13 introduced cultivars were analyzed by using the grey system theory.[Result]The weighted correlation degree of the 13 tested varieties was larger than that of ‘Bijie Chaotianjiao,accounting for 28.5% of the tested varieties.Among them,‘Jinjiaogan 6 ranked first in weighted correlation degree and yield,and had the best comprehensive performance.‘Nongwangzhitian No.1‘ and ‘Yanjiao 425had the second comprehensive performance.[Conclusion]‘Jinjiaogan 6‘Yanjiao 425 and ‘Nongwangzhitian No. 1 were all suitable for cultivation in Bijie.‘Jinjiaogan 6was most suitable for introduction and cultivation in Bijie.

        Key words Grey correlation degree;Bijie;Capsicum annuum;Introduction

        基金項目 貴州省科技廳農(nóng)業(yè)攻關(guān)項目(黔科合NY〔2013〕3024號)。

        作者簡介 謝龍安(1985—),男,貴州畢節(jié)人,講師,碩士,從事蔬菜遺傳育種與栽培技術(shù)研究。*通信作者,副教授,博士,從事蔬菜遺傳育種與栽培技術(shù)研究。

        收稿日期 2022-07-19

        朝天椒(Capsicum annuum var.frutescens L.)為茄科(Solanaceae)辣椒屬(Capsicum L.),又名三櫻椒、天鷹椒,包括簇生椒、圓錐椒(小果型)、長辣椒(短指形)、櫻桃椒的4個變種[1-2],其椒果小、辣度高、易干制,既可鮮食、泡制,又可作為干辣椒食用,與羊角椒、線椒構(gòu)成我國三大干椒品種系列[3]。朝天椒種植面積逐年上升,已躍居干椒栽培面積首位[4-5],研究朝天椒引種栽培試驗,篩選適宜畢節(jié)市朝天椒種植品種意義重大。目前,朝天椒引種和品種比較試驗主要是應用統(tǒng)計分析,對產(chǎn)量進行豐產(chǎn)穩(wěn)定性分析[6-7]。自鄧聚龍[8]提出灰色系統(tǒng)理論以來,灰色關(guān)聯(lián)度分析廣泛運用于多種農(nóng)作物引種和育種篩選上[9]。該理論是一種用特定的方法描述信息,對不完全的系統(tǒng)進行預測、決策、控制的嶄新的系統(tǒng)理論[10-11],是一種新的多因素分析方法[12-13],為全面、合理地評價品種的優(yōu)劣提供了一條新途徑[14-19],其分析結(jié)果與定性分析結(jié)果比較吻合[8]。近年來,畢節(jié)市政府加大農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)調(diào)整力度,大力支持辣椒產(chǎn)業(yè)發(fā)展,助力鄉(xiāng)村產(chǎn)業(yè)振興,畢節(jié)市朝天椒市場交易量不斷增加,種植面積逐年擴大,朝天椒種植已經(jīng)成為部分鄉(xiāng)鎮(zhèn)農(nóng)民增收的主要渠道。但對畢節(jié)市朝天椒引種進行全面、客觀評價的報道未見。筆者利用畢節(jié)引種朝天椒品種為試驗資料,應用灰色關(guān)聯(lián)度分析法綜合評價14個朝天椒品種,旨在篩選適宜畢節(jié)市推廣種植的朝天椒品種,為畢節(jié)市朝天椒栽培提供理論依據(jù)。

        1 材料與方法

        1.1 試驗材料

        參試材料為“三都縣朝天椒1號”(X1)、“三都縣朝天椒2號”(X2)、“六枝特區(qū)朝天椒”(X3)、“沿河縣朝天椒1號”(X4)、“沿河縣朝天椒2號”(X5)、“艷椒425”(X6)、“農(nóng)望指天1號”(X7)、“泰國朝天椒”(X8)、“傲雪小米椒”(X9)、“川椒168”(X10)、“金嬌干6號”(X11)、“農(nóng)望玉嬌子”(X12)、“千斤子彈頭”(X13)、“畢節(jié)朝天椒”(X14,CK)。

        1.2 試驗方法

        試驗地選擇在貴州省畢節(jié)市七星關(guān)區(qū)生機鎮(zhèn)峨峰村,總面積約0.08 hm2,各試驗小區(qū)土地肥力均等。田間試驗采用隨機區(qū)組排列,每個品種3次重復,共42個小區(qū),每區(qū)行長15 m,寬1.2 cm,行距60 cm,株距40 cm,種植密度為34 125株/hm2,同期播種,同法栽培管理,對參試種在引種栽培地的生長發(fā)育期、植株農(nóng)藝性狀、果實外觀品質(zhì)、內(nèi)在品質(zhì)、產(chǎn)量、抗病性及抗逆性等進行記載、觀測及檢測。測定性狀包括果實成熟率(%)(K1)、生育期(d)(K2)、植株分枝數(shù)(個)(K3)、植株主莖高(cm)(K4)、植株主莖徑粗(cm)(K5)、果實橫徑(cm)(K6)、果長(cm)(K7)、單株果實數(shù)(個)(K8)、單株鮮果重量(g)(K9)、鮮果含水率(%)(K10)、干椒總產(chǎn)量(kg)(K11)、紅熟干椒產(chǎn)量(kg)(K12)、抗倒伏率(%)(K13)、抗病率(%)(K14)、可溶性糖含量(mg/kg)(K15)、可溶性蛋白含量(mg/kg)(K16)、VC含量(mg/kg)(K17)、辣椒素含量(mg/g)(K18)、二氫辣椒素含量(mg/g)(K19)。

        1.3 分析方法

        (1)確定理想品系。參考鄧聚龍[8]的灰色關(guān)聯(lián)分析法,將14個參試品系作為一個灰色系統(tǒng),每品系為系統(tǒng)中的一個因素。把引種朝天椒與參試品種的理想指標結(jié)合起來,構(gòu)建一個理想品種,理想品種的主要性狀指標為供試品種表現(xiàn)最好的指標,以理想品種的各性狀指標構(gòu)成參考數(shù)列,記為X0。其中,Xi(i =1,2,…,m),其中m為參試品種數(shù),各性狀用Ki表示(j=1,2,…,n),其中n 為選取分析的指標數(shù)。

        (2)數(shù)據(jù)無量綱化處理。因各分析因素量綱不一致,需要對原始數(shù)據(jù)進行無量綱化處理,利用初值化法將數(shù)據(jù)進行無量綱化處理,即用Xi數(shù)值分別除以X0,記作:

        xi(k)=xi′(k)xi′(1),i=0,1,…,n;k=1,2,…,m.(1)

        (3)計算絕對差值及最大、最小值。根據(jù)無量綱化值計算參考數(shù)列與比較數(shù)列的絕對差值,記作:

        △i(k)=︱X0(K) -Xi(k)︱(2)

        其中,X0(K)為理想品系初始化均值,Xi(k)為第i個參試品種第k個性狀的初始化均值。

        用公式(3)分別計算每個比較序列與參考序列對應元素的關(guān)聯(lián)系數(shù):

        ζi(k)=min|x0′(k)-xi′(k)|+ρ·maxi|x0′(k)-xi′(k)|

        |x0′(k)-xi′(k)|+ρ·maxi|x0′(k)-xi′(k)|(k=1,…,m)(3)

        其中,ζi(k) 為Xi對X0在K點的關(guān)聯(lián)系數(shù),ρ為分辨系數(shù),取值范圍在(0~1),該試驗數(shù)據(jù)分析取值為ρ=0.5。

        (4)計算關(guān)聯(lián)系數(shù)。由公式(4)計算試驗品系與理想品系之間的關(guān)聯(lián)系數(shù)ri:

        ri=1mmk=1ζi(k)(4)

        對所計算關(guān)聯(lián)度均值進行排序比較。

        2 結(jié)果與分析

        2.1 參試品種的性狀指標和理想品系

        14個品種的果實成熟率32.14%~96.65%,生育期175~220 d,植株分枝數(shù)3.2~10.0個,植株主莖高度4.58~15.95 cm,植株主莖徑粗0.595~1.320 cm,果實橫徑5.020~15.760 cm,果長40.375~80.300 cm,單株果實數(shù)34~153個,單株鮮果重122.50~472.00 g,鮮果含水量62.54%~72.93%,干椒總產(chǎn)量98.17~453.91 kg,紅熟干椒產(chǎn)量58.69~352.33 kg,抗倒伏率55.00%~97.00%,抗病率83.00%~99.30%,可溶性糖含量2 859.1~9 461.7 mg/kg,可溶性蛋白含量4 100.5~14 326.0 mg/kg,VC含量530.3~951.9 mg/kg,辣椒素含量0.440~2.591 mg/g,二氫辣椒素含量0.181~1.653 mg/g(表1)。

        根據(jù)灰色關(guān)聯(lián)分析方法,將理想品系X0和14個參試品種的指標數(shù)值列于表1。理想品系X0在該試驗中,果實成熟度取最大值96.65%,植株生育期取最小值175 d,植株分枝數(shù)取最大值10.0 cm,植株分枝點高度取10.00 cm(因植株分枝點過低,植株冠幅小,植株分枝點過高,植株分枝數(shù)較少,故取中間數(shù)字);植株主莖徑粗取最大值1.320 cm,果實橫徑取最大值15.760 cm,果長取最大值80.300 cm;鮮果可溶性糖含量取最大值9 461.7 mg/kg,鮮果可溶性蛋白含量取最大值14 326.0 mg/kg,鮮果VC含量取最大值951.9 mg/kg,干椒辣椒素含量取最大值2.951? mg/g,干椒二氫辣椒素含量取最大值1.653? mg/g,抗病率取最大值99.30%,抗倒伏率取最大值97.00%,鮮果含水量取最小值62.54%,單株果實數(shù)取最大值153個,單株鮮果重取最大值472.00 g,干椒總產(chǎn)量取最大值453.91 kg,紅熟干椒產(chǎn)量取最大值352.33 kg。

        2.2 數(shù)據(jù)無量綱化處理

        由于朝天椒不同性狀值差異較大,且不同性狀間單位不同,為了方便關(guān)聯(lián)分析,采用初值化法對14個供試品種的各性狀原始數(shù)據(jù)進行無量綱化處理,即用14個供試品種各個評價指標的原始數(shù)據(jù)分別去除理想品種對應的評價指標,得到表2的無量綱化數(shù)列。無量綱化后,K1為0.33~1.00,K2為1.00~1.26,K3為0.3~1.0,K4為0.46~1.60,K5為0.451~1.000,K6為0.319~1.000,K7為0.503~1.000,K8為0.2~1.0,K9為0.3~1.0,K10為1.00~1.17,K11為0.216~1.000,K12為0.167~1.000,K13為0.57~1.00,K14為0.811~1.000,K15為0.301~1.000,K16為0.286~1.000,K17為0.557~1.000,K18為0.323~1.000,K19為0.251~1.000。

        2.3 求絕對差值及最大、最小差值

        根據(jù)表2的數(shù)據(jù),計算參考數(shù)列與比較數(shù)列的絕對差值,即Δi(k)=| X|0(k)-Xi(k) (i=1,2,…,14;k=1,2,…,19),X0 (k)與Xi(k)各對應點的絕對差值形成數(shù)列,見表3。K1絕對差值0~0.67,K2絕對差值0~0.26,K3絕對差值0~0.7,K4絕對差值0~0.60,K5絕對差值0~0.549,K6絕對差值0~0.681,K7絕對差值0~0.497,K8絕對差值0~0.8,K9絕對差值0~0.7,K10絕對差值0~0.17,K11絕對差值0~0.784,K12絕對差值0~0.83,K13絕對差值0~0.43,K14絕對差值0~0.189,K15絕對差值0~0.699,K16絕對差值0~0.714,K17絕對差值0~0.401,K18絕對差值0~0.851,K19絕對差值0~0.891。

        2.4 關(guān)聯(lián)系數(shù)及加權(quán)關(guān)聯(lián)度計算

        將表3的數(shù)據(jù)帶入公式(3)中,計算出供試品種與理想品種的關(guān)聯(lián)系數(shù),其中,Δi(k)表示絕對差值,ρ表示分辨系數(shù)(0<ρ<1),該研究ρ=0.5,計算各參試品種與理想品種的關(guān)聯(lián)系數(shù)(表4)。K1關(guān)聯(lián)系數(shù)0.33~1.00,K2關(guān)聯(lián)系數(shù)0.34~1.00,K3關(guān)聯(lián)系數(shù)0.3~1.0,K4關(guān)聯(lián)系數(shù)0.34~0.85,K5關(guān)聯(lián)系數(shù)0.334~1.000,K6關(guān)聯(lián)系數(shù)0.333~1.000,K7關(guān)聯(lián)系數(shù)0.335~1.000,K8關(guān)聯(lián)系數(shù)0.3~1.0,K9關(guān)聯(lián)系數(shù)0.3~1.0,K10關(guān)聯(lián)系數(shù)0.33~1.00,K11關(guān)聯(lián)系數(shù)0.332~1.000,K12關(guān)聯(lián)系數(shù)0.33~1.00,K13關(guān)聯(lián)系數(shù)0.33~1.00,K14關(guān)聯(lián)系數(shù)0.322~1.000,K15關(guān)聯(lián)系數(shù)0.334~1.000,K16關(guān)聯(lián)系數(shù)0.333~1.000,K17關(guān)聯(lián)系數(shù)0.311~1.000,K18關(guān)聯(lián)系數(shù)0.333~1.000,K19關(guān)聯(lián)系數(shù)0.333~1.000。

        根據(jù)公式(4),計算供試品種與理想品種的關(guān)聯(lián)度,計算出供試品種與理想品種的關(guān)聯(lián)度,并排次位,排序順序見表5。根據(jù)灰色關(guān)聯(lián)度分析原則,關(guān)系度數(shù)值越大的品系與參考品系越接近,關(guān)聯(lián)度數(shù)值越大,品種綜合性狀越好,反之,綜合性狀越差。由表5可知,X11的加權(quán)關(guān)聯(lián)度最大,引種綜合表現(xiàn)最好,X6、X7加權(quán)關(guān)聯(lián)度大于0.70,引種綜合表現(xiàn)次之,X12加權(quán)關(guān)聯(lián)度大于0.60,X14、X13、X1、X2、X4、X3、X9、X10、X5、X8加權(quán)關(guān)聯(lián)度均小于0.60,加權(quán)關(guān)聯(lián)度較小,引種綜合表現(xiàn)差。

        3 討論

        通過對14個畢節(jié)市引種朝天椒物候期、植株農(nóng)藝性狀、果實外觀品質(zhì)、果實內(nèi)在品質(zhì)、產(chǎn)量、抗病蟲害、抗倒伏性等19指標進行灰色關(guān)聯(lián)分析,運用灰色關(guān)聯(lián)度對14個參試朝天椒品種進行綜合評價,加權(quán)關(guān)聯(lián)度位次排名與引種實際表現(xiàn)及產(chǎn)量基本保持一致,關(guān)聯(lián)系數(shù)大小為X11>X7>X6>X12>X14>X13>X1>X2>X4>X3>X9>X10>X5>X8。X11(r=0.776 7)加權(quán)關(guān)聯(lián)度最大,X7(r=0.716 2)、X6(r=0.713 4)加權(quán)關(guān)聯(lián)度次之,說明這3個品種綜合表現(xiàn)良好;X8(r=0.446 8)加權(quán)關(guān)聯(lián)度最小,說明該品種綜合表現(xiàn)最差。主要經(jīng)濟指標表現(xiàn)為X11產(chǎn)量排第一,X6次之,X14最差。該試驗結(jié)果與加權(quán)關(guān)聯(lián)度系數(shù)排序結(jié)果基本一致,灰色關(guān)聯(lián)度分析方法能全面、真實地反映參試品種的實際表現(xiàn),灰色關(guān)聯(lián)度分析法應用于朝天椒引種、品比試驗綜合評價中,有效地提高綜合評價的準確性,可更全面為引種試驗、優(yōu)質(zhì)品種推廣種植提供可靠依據(jù)。灰色關(guān)聯(lián)度分析能準確判斷參試品種優(yōu)劣,篩選出適宜推廣的參試品種;加權(quán)關(guān)聯(lián)度值能真實反映品種的綜合表現(xiàn),可應用于朝天椒新品種的選育和引種中。

        4 結(jié)論

        通過對畢節(jié)市引種朝天椒灰色關(guān)聯(lián)分析,若不考慮物候期、植株農(nóng)藝性狀、果實外觀品質(zhì)、果實內(nèi)在品質(zhì)、產(chǎn)量及抗病抗倒伏性各個指標之間的權(quán)重關(guān)系,各個觀測指標同等重要,“金嬌干6號”“艷椒425”“農(nóng)望指天1號”3種朝天椒品種均適合在畢節(jié)種植,其余朝天椒不適宜在畢節(jié)種植;單從朝天椒主要經(jīng)濟效益指標產(chǎn)量看,“金嬌干6號”最適宜畢節(jié)種植,建議在畢節(jié)地區(qū)推廣種植。因畢節(jié)市范圍廣,海拔、氣溫、日照和降雨等存在差異,該試驗存在一定的局限性,試驗地外的其他地點引種表現(xiàn)還需進一步研究。

        參考文獻

        [1] GUTSU E V,KINTYA P K,SHVETS S A,et al.Steroid glycosides of the roots of Capsicum annuum.Ⅰ.The structure of capsicosides A1,B1,and C1[J].Chemistry of natural compounds,1986,22(6):661-664.

        [2] 李璠.中國栽培植物發(fā)展史[M].北京:科學出版社,1984:130-131.

        [3] 戴雄澤.辣椒制種技術(shù)[M].北京:中國農(nóng)業(yè)出版社,2001:119-115.

        [4] 桂敏,龍洪進,鐘秋月,等.朝天椒種質(zhì)資源主要農(nóng)藝性狀的相關(guān)性和通徑分析[J].西南農(nóng)業(yè)學報,2017,30(1):176-182.

        [5] 高俊山,田曉遠,董高峰,等.河南省朝天椒產(chǎn)業(yè)發(fā)展前景分析與建議[J].安徽農(nóng)業(yè)科學,2016,44(15):225-227.

        [6] 趙加瑞,師曉丹,李敬勛,等.河南省春種辣椒新品種露地栽培比較試驗[J].安徽農(nóng)學通報,2019,25(15):58-60.

        [7] 趙伯蓮,孫松英,張忠武,等.瓦爾崗七星椒與13 個朝天椒品種的比較與評價[J].湖南農(nóng)業(yè)科學,2019(4):8-11.

        [8] 鄧聚龍.農(nóng)業(yè)系統(tǒng)灰色理論與方法[M].濟南:山東科學技術(shù)出版社,1988.

        [9] 陳雷,吳繼華,李可,等.基于灰色關(guān)聯(lián)度分析的河南省麥套花生研究[J].中國農(nóng)學通報,2019,35(8):5-11.

        [10] 鄧聚龍.灰色系統(tǒng)理論教程[M] .武漢:華中理工大學出版社,1990.

        [11] 鄧聚龍.灰色系統(tǒng)基本方法[M].武漢:華中理工大學出版社,1987.

        [12] 劉思峰,黨耀國,方志耕,等.灰色系統(tǒng)理論及其應用[M].北京:科學出版社,2004.

        [13] 王佳敏,劉敏,郭詠梅,等.AHP法和灰色關(guān)聯(lián)法在觀賞辣椒果實外觀品質(zhì)評價中的應用[J].山西農(nóng)業(yè)大學學報(自然科學版),2019,39(1):73-78.

        [14] 彭詩云,羅燚,謝挺,等.灰色關(guān)聯(lián)度分析法在遵椒4 號密度上的應用[J].辣椒雜志,2016,14(4):32-35.

        [15] 黃松,高路銀,楊森要,等.薄皮甜瓜糖酸含量及其與果實性狀的灰色關(guān)聯(lián)度分析[J].河南農(nóng)業(yè)科學,2018,47(7):112-115.

        [16] 呂霞,劉慧英,祁樂,等.基于灰色關(guān)聯(lián)度的加工番茄新品系綜合評價研究[J].安徽農(nóng)業(yè)科學,2013,41(11):4770-4772,4776.

        [17] 喬鋒,羅英,曾紹貴,等.基于灰色系統(tǒng)理論的福建辣椒地方品種產(chǎn)量性狀綜合評估[J].福建農(nóng)業(yè)學報,2009,24(4):318-322.

        [18] 張秀榮,張保才,劉崢.越冬甘藍主要農(nóng)藝性狀的灰色關(guān)聯(lián)度分析及綜合評價[J].江西農(nóng)業(yè)學報,2019,31(6):26-31.

        [19] 艾尼瓦爾,宋文勝,曹連莆,等.應用灰色關(guān)聯(lián)度分析評價制干辣椒品系區(qū)域試驗[J].北方園藝,2005(1):51-52.

        猜你喜歡
        灰色關(guān)聯(lián)度朝天椒引種
        朝天椒套種芝麻輕簡高效栽培技術(shù)
        節(jié)節(jié)紅板栗引種表現(xiàn)及其優(yōu)質(zhì)豐產(chǎn)栽培技術(shù)
        森林碳匯影響因素的灰色關(guān)聯(lián)分析
        基于灰色層次分析法的艦艇衛(wèi)星通信效能評估方法研究
        基于AHP灰色關(guān)聯(lián)分析的企業(yè)電子商務信用評價研究
        商(2016年25期)2016-07-29 21:07:14
        湖南省經(jīng)濟增長影響因素實證分析
        中國市場(2016年21期)2016-06-06 04:07:01
        韓國朝天椒新品種PE101的引種及栽培技術(shù)
        中國蔬菜(2015年9期)2015-12-21 13:04:44
        紅紅的朝天椒
        小說月刊(2014年1期)2014-04-23 08:59:51
        朝天椒包餃子
        小青蛙報(2009年2期)2009-04-17 09:33:12
        人人妻人人澡人人爽国产| AV中文字幕在线视| 男女后入式在线观看视频| 亚洲国产精品不卡av在线| 免费人妻精品一区二区三区| 亚洲国产成人久久综合一区77| 亚洲AV专区一专区二专区三 | 五月激情在线视频观看| 国产精品无码一区二区在线观一 | 亚洲av无码专区在线观看成人| 2019年92午夜视频福利| 一区二区三区国产高潮| 国产高清视频在线不卡一区| 人妻丰满熟妇岳av无码区hd| 国产精品丝袜黑色高跟鞋| 国产精品人成在线观看| 男女激情视频网站在线| 性无码专区无码| 亚洲AV电影天堂男人的天堂| 国产免费人成网站在线播放| 91自拍视频国产精品| 国产如狼似虎富婆找强壮黑人| 囯产精品无码va一区二区| 久久人妻精品中文字幕一区二区| 亚洲最新无码中文字幕久久 | 亚洲中文字幕诱惑第一页| 国产亚洲一区二区三区| 国产真实偷乱视频| 日本加勒比东京热日韩| 青青青爽在线视频免费播放| 男人和女人做爽爽免费视频| 国产精品国语对白露脸在线播放| 日本午夜理伦三级好看| 综合色免费在线精品视频| 一个人看的视频www免费| 9久9久女女热精品视频免费观看| 日韩激情视频一区在线观看| 亚洲中文字幕一区精品自拍| 国产精品国产三级农村妇女| 国内专区一区二区三区| 色与欲影视天天看综合网|