999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

魯西北地區(qū)夏谷產(chǎn)量與主要農(nóng)藝性狀的灰色關(guān)聯(lián)分析

2025-01-17 00:00:00王春雨朱冠雄高祺曹鵬鵬李娜娜朱金英田藝心李春燕高鳳菊華方靜
安徽農(nóng)業(yè)科學(xué) 2025年1期
關(guān)鍵詞:灰色關(guān)聯(lián)分析產(chǎn)量

摘要為促進(jìn)魯西北地區(qū)優(yōu)良谷子資源的利用,采用灰色關(guān)聯(lián)分析方法,對(duì)魯西北地區(qū)17個(gè)夏谷新品種(系)產(chǎn)量與農(nóng)藝性狀進(jìn)行研究。調(diào)查和統(tǒng)計(jì)分析發(fā)現(xiàn),影響谷子產(chǎn)量的7個(gè)農(nóng)藝性狀的關(guān)聯(lián)度順序?yàn)樯趃t;株高gt;穗數(shù)gt;單穗重gt;穗粒重gt;穗長(zhǎng)gt;千粒重,其中生育期、株高、穗數(shù)、單穗重、穗粒重5個(gè)性狀對(duì)谷子產(chǎn)量的關(guān)聯(lián)度較高。因此,為促進(jìn)魯西北地區(qū)夏谷優(yōu)質(zhì)高產(chǎn)發(fā)展,應(yīng)在滿足正常發(fā)育成熟的前提下,選擇生育期較長(zhǎng)、株高較高、穗數(shù)較多的谷子品種。

關(guān)鍵詞產(chǎn)量;農(nóng)藝性狀;灰色關(guān)聯(lián)分析;夏谷

中圖分類號(hào)S515"文獻(xiàn)標(biāo)識(shí)碼A"文章編號(hào)0517-6611(2025)01-0022-04

doi:10.3969/j.issn.0517-6611.2025.01.006

開放科學(xué)(資源服務(wù))標(biāo)識(shí)碼(OSID):

GreyCorrelationAnalysisofSummerMilletYieldandMainAgronomicCharactersinNorthwestShandongProvince

WANGChun-yu,ZHUGuan-xiong,GAOQietal

(DezhouAcademyofAgriculturalSciences,Dezhou,Shandong253015)

AbstractTopromotetheutilizationofsuperiormilletresourcesinnorthwestShandong,theyieldandagronomiccharactersof17newsummermilletvarietiesinnorthwestShandongwerestudiedbygreycorrelationanalysis.Throughinvestigationandstatisticalanalysis,itwasfoundthatamongtheagronomictraitsaffectinggrainyield,thecorrelationsequenceofthesevenagronomictraitsrelatedtograinyieldwasintheorderofgrowthdategt;plantheightgt;earnumbergt;earweightperspikegt;earweightperspikegt;earlengthgt;1000-grainweight.Amongthem,growthdate,plantheight,earnumber,earweightperspikeandgrainweightperspikehadhighercorrelationwithgrainyield.Therefore,inordertopromotethedevelopmentofhighqualityandhighyieldofsummermilletinnorthwestShandong,weshouldchoosethemilletvarietieswithlongergrowthperiod,higherplantheightandmoreearunderthepremiseofmeetingthenormaldevelopmentandmaturity.

KeywordsYield;Agronomiccharacter;Greycorrelationanalysis;Summermillet

基金項(xiàng)目山東省現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系雜糧創(chuàng)新團(tuán)隊(duì)建設(shè)項(xiàng)目(SDAIT-5-01);山東省農(nóng)業(yè)良種工程項(xiàng)目(2021LZGC025);農(nóng)作物種質(zhì)資源挖掘與精準(zhǔn)鑒定項(xiàng)目(2021LZGC025)。

作者簡(jiǎn)介王春雨(1995—),女,山東德州人,助理農(nóng)藝師,碩士,從事大豆及雜糧栽培育種。*通信作者,高級(jí)農(nóng)藝師,碩士,從事食用豆及雜糧栽培育種研究。

谷子是起源于中國(guó)的禾本科狗尾草屬作物,又稱粟,剝殼后稱小米,在我國(guó)北方的干旱和半干旱地區(qū)均有種植[1-2。谷子具有耐瘠、耐旱、抗逆、節(jié)水的特點(diǎn),是環(huán)境友好型植物3-5。谷子成熟后含有豐富的蛋白質(zhì)、脂肪、維生素、氨基酸、膳食纖維等人體必需的營(yíng)養(yǎng)元素,其蛋白質(zhì)含量比玉米高30%左右,是糧飼兼用型作物[6-10。谷子的產(chǎn)量受遺傳特性、栽培條件、栽培方法、環(huán)境條件等多種因素的影響,這些因素對(duì)谷子產(chǎn)量的影響程度存在較大差別11-15

灰色關(guān)聯(lián)分析是由鄧聚龍教授創(chuàng)立的一種研究少數(shù)據(jù)、貧信息不確定性問(wèn)題的新方法[16。常用的方差分析是對(duì)產(chǎn)量進(jìn)行單因素評(píng)價(jià),而利用灰色關(guān)聯(lián)分析能夠更直觀地了解主要農(nóng)藝性狀對(duì)產(chǎn)量形成的相對(duì)重要性以及性狀與產(chǎn)量之間的相互作用,有利于育種者在實(shí)際育種過(guò)程中提高親本選配和后代選擇的效果,能夠一定程度減少因個(gè)人認(rèn)識(shí)程度和喜好不同等主觀因素對(duì)新品種選育的準(zhǔn)確性造成的負(fù)面影響。目前,灰色關(guān)聯(lián)分析方法在玉米17-18、大豆19-20、小麥21、水稻22等作物上應(yīng)用較多,但用于谷子育種上的報(bào)道較少。鑒于此,筆者對(duì)魯西北地區(qū)夏播谷子的主要農(nóng)藝性狀進(jìn)行灰色關(guān)聯(lián)分析,探討主要農(nóng)藝性狀對(duì)產(chǎn)量形成所起的作用和相對(duì)重要性,旨在為篩選適宜魯西北地區(qū)種植的優(yōu)良谷子品種以及谷子新品種選育提供理論依據(jù)。

1材料與方法

1.1試驗(yàn)地概況

試驗(yàn)于2022年在山東省德州市黃河涯鎮(zhèn)現(xiàn)代農(nóng)業(yè)科技園示范基地(37°35′N,116°30′E)進(jìn)行。該地地處黃淮海平原,年日照時(shí)數(shù)2451.21~2592.13h,年均氣溫12.3~12.9℃,年降水量540.5~562.3mm,試驗(yàn)地地勢(shì)平坦,排灌方便,土壤主要為砂質(zhì)壤土。基本理化性質(zhì)如下:土壤有機(jī)質(zhì)含量12.82g/kg、全氮0.91g/kg、速效氮57.18mg/kg、有效磷66.74mg/kg、速效鉀102.37mg/kg,pH7.02。前茬作物為小麥,平均產(chǎn)量9750kg/hm施復(fù)合基肥600kg/hm無(wú)追肥。

1.2試驗(yàn)材料

數(shù)據(jù)材料來(lái)源于山東谷子生態(tài)適應(yīng)性試驗(yàn),試驗(yàn)品種(系)17個(gè),來(lái)源于山東省農(nóng)業(yè)科學(xué)院和淄博市農(nóng)業(yè)科學(xué)研究院(表1)。

1.3試驗(yàn)方法

1.3.1試驗(yàn)設(shè)計(jì)。以17個(gè)谷子品種(系)為試驗(yàn)材料,采用隨機(jī)區(qū)組排列,3次重復(fù)。6行區(qū),行長(zhǎng)5.0m,寬3.6m,密度45萬(wàn)株/hm試驗(yàn)四周設(shè)保護(hù)行。

按照試驗(yàn)標(biāo)準(zhǔn)進(jìn)行各項(xiàng)農(nóng)藝性狀的調(diào)查記載,生育期調(diào)查記載播種期、出苗期、抽穗期、成熟期;收獲前,在每個(gè)小區(qū)選擇不缺苗、長(zhǎng)勢(shì)一致的行確定10株谷子取樣,調(diào)查株高并將其整株取回,以備考種;收獲后室內(nèi)考種,測(cè)定穗長(zhǎng)、單穗重、穗粒重、千粒重。整個(gè)小區(qū)全部收獲,測(cè)定產(chǎn)量和穗數(shù)。

1.3.2數(shù)據(jù)分析。采用Excel2019對(duì)數(shù)據(jù)進(jìn)行整理;采用DPS數(shù)據(jù)處理系統(tǒng)對(duì)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析。關(guān)聯(lián)系數(shù)、關(guān)聯(lián)度越大表示該性狀與產(chǎn)量關(guān)系越密切,反之,則該性狀與產(chǎn)量關(guān)系越弱。

2結(jié)果與分析

2.1數(shù)據(jù)無(wú)量綱化處理

表2為不同谷子品種(系)主要農(nóng)藝性狀比較。由表2可知,17個(gè)供試品種生育期為86.5~88.0d,產(chǎn)量為8025.75~9174.60kg/hm2。

將17個(gè)供試品種(系)視為1個(gè)灰色系統(tǒng),把產(chǎn)量作為參考數(shù)列,記為X0(k);把生育期、穗數(shù)、株高、穗長(zhǎng)、單穗重、穗粒重、千粒重為比較數(shù)列,記為Xi(k)(i="…,m;k="…,n),其中,m為與產(chǎn)量相關(guān)的性狀的個(gè)數(shù),n為參試品種(系)個(gè)數(shù)。由于各個(gè)性狀的量綱不同,無(wú)法直接進(jìn)行比較,因此需要對(duì)數(shù)據(jù)進(jìn)行無(wú)量綱化處理,轉(zhuǎn)化為可以直接比較的標(biāo)準(zhǔn)化數(shù)列。先對(duì)各個(gè)性狀進(jìn)行簡(jiǎn)單的統(tǒng)計(jì)分析,計(jì)算出平均值和標(biāo)準(zhǔn)差,然后將各個(gè)原始數(shù)據(jù)減去平均值后再除以標(biāo)準(zhǔn)差,得到標(biāo)準(zhǔn)化數(shù)列(表3)。

2.2比較數(shù)列與參考數(shù)列的絕對(duì)差值

根據(jù)標(biāo)準(zhǔn)化處理結(jié)果,按照公式Δ"0i(k)=∣X0(k)-Xi(k)∣(i="…,7;k="…,17),計(jì)算比較數(shù)列與參考數(shù)列的絕對(duì)差值(表4)。

2.3關(guān)聯(lián)系數(shù)

經(jīng)數(shù)據(jù)標(biāo)準(zhǔn)化處理后的參考數(shù)列記作X0(k),比較數(shù)列記作Xi(k),關(guān)聯(lián)系數(shù)ξ0i(k)按照如下公式計(jì)算:

ξ0i(k)=Δ(min)+ρΔ(max)ξ0i(k)+ρΔ(max)

式中:Δ表示比較序列與參考數(shù)列的絕對(duì)差值,關(guān)聯(lián)系數(shù)反映比較序列和參考數(shù)列的關(guān)聯(lián)程度;ρ取值0.5。

計(jì)算結(jié)果見表5。

2.4關(guān)聯(lián)度和排序計(jì)算公式如下:

ri=1Nkiξ0i(k)

7個(gè)農(nóng)藝性狀與產(chǎn)量的關(guān)聯(lián)度排序?yàn)樯趃t;株高gt;穗數(shù)gt;單穗重gt;穗粒重gt;穗長(zhǎng)gt;千粒重(表6),表明生育期、株高、穗數(shù)、單穗重、穗粒重對(duì)谷子產(chǎn)量影響較大;穗長(zhǎng)、千粒重對(duì)谷子產(chǎn)量影響較小。這與趙禹凱等[23研究結(jié)果類似。

2.5谷子產(chǎn)量相關(guān)性狀之間的相互關(guān)系

由表7可知,與生育期關(guān)系最密切的農(nóng)藝性狀是株高,其次是穗數(shù)和穗粒重;與穗數(shù)關(guān)系最密切的農(nóng)藝性狀是穗粒重,其次是千粒重和單穗重;與株高關(guān)系最密切的農(nóng)藝性狀是生育期,其次是穗數(shù)和單穗重;與穗長(zhǎng)關(guān)系最密切的農(nóng)藝性狀是穗粒重,其次是株高和穗數(shù);與單穗重關(guān)系最密切的農(nóng)藝性狀是穗粒重,其次是穗數(shù)和株高;與穗粒重關(guān)系最密切的農(nóng)藝性狀是單穗重,其次是穗數(shù)和生育期;與千粒重關(guān)系最密切的農(nóng)藝性狀是穗數(shù),其次是株高和穗粒重。

3結(jié)論與討論

通過(guò)灰色關(guān)聯(lián)分析方法,對(duì)供試的17個(gè)谷子品種(系)的7個(gè)主要農(nóng)藝性狀與產(chǎn)量關(guān)系的密切程度進(jìn)行分析。結(jié)果表明,與谷子產(chǎn)量相關(guān)的7個(gè)農(nóng)藝性狀的關(guān)聯(lián)度順序?yàn)樯趃t;株高gt;穗數(shù)gt;單穗重gt;穗粒重gt;穗長(zhǎng)gt;千粒重。其中生育期、株高、穗數(shù)、單穗重、穗粒重與產(chǎn)量的關(guān)聯(lián)系數(shù)較高,說(shuō)明生育期、株高、穗數(shù)、單穗重、穗粒重對(duì)產(chǎn)量形成產(chǎn)生的影響較大,而穗長(zhǎng)和千粒重對(duì)產(chǎn)量的影響較小。因此,在新品種選育過(guò)程中要綜合考慮生育期、株高、穗數(shù)、單穗重、穗粒重5個(gè)性狀,結(jié)合當(dāng)?shù)貧夂驐l件,在滿足谷子生長(zhǎng)發(fā)育、谷子能夠成熟的前提下,盡量選擇生育期較長(zhǎng)、株高較高、穗數(shù)較多的品種。

該研究結(jié)果表明,生育期與谷子產(chǎn)量的關(guān)系最為密切,其次是株高,趙禹凱等[23以東北地區(qū)的春播谷子為試驗(yàn)材料,研究了來(lái)源于黑龍江省、吉林省、遼寧省、河南省和內(nèi)蒙古自治區(qū)的22個(gè)谷子品種,結(jié)果表明生育期與谷子產(chǎn)量的關(guān)系較為密切,但株高與產(chǎn)量的密切關(guān)系較弱;而欒素榮等[24以西北地區(qū)的中晚熟谷子為試驗(yàn)材料,結(jié)果顯示生育期與谷子產(chǎn)量的關(guān)系最不密切。這種結(jié)果上的差異可能是試驗(yàn)材料不同、種植氣候條件不同等因素造成的。因此要針對(duì)育種者當(dāng)?shù)氐膶?shí)際情況進(jìn)行試驗(yàn)后分析得到的結(jié)果才更具指導(dǎo)意義。生育期是一個(gè)特殊的指標(biāo),受地域及種植當(dāng)年的氣象條件影響較大,因此生育期可以作為一個(gè)長(zhǎng)期選擇的重要指標(biāo),在滿足谷子成熟的條件下,盡可能選擇生育期較長(zhǎng)的品種是增加谷子產(chǎn)量的有效途徑。

灰色關(guān)聯(lián)分析具有計(jì)算簡(jiǎn)單、結(jié)果直觀的優(yōu)點(diǎn),是研究多元性狀相對(duì)重要性的良好方法[25。但由于其分析的系統(tǒng)為動(dòng)態(tài)系統(tǒng),分析結(jié)果在不同的時(shí)間地點(diǎn)有所變化,育種者在品種選育過(guò)程中應(yīng)結(jié)合實(shí)際情況綜合考慮。

參考文獻(xiàn)

[1]溫蕊,劉宇飛,王磊,等.內(nèi)蒙古西部地區(qū)谷子產(chǎn)業(yè)發(fā)展現(xiàn)狀及展望[J].現(xiàn)代農(nóng)業(yè),2021(5):77-78.

[2]趙志軍.中國(guó)農(nóng)業(yè)起源概述[J].遺產(chǎn)與保護(hù)研究,2019,4(1):1-7.

[3]張雪峰.中國(guó)谷子產(chǎn)業(yè)發(fā)展問(wèn)題研究[D].哈爾濱:東北農(nóng)業(yè)大學(xué),2013.

[4]董揚(yáng),李清泉,閆鋒,等.谷子在畜牧業(yè)中的應(yīng)用優(yōu)勢(shì)及栽培技術(shù)[J].現(xiàn)代畜牧科技,20240(6):56-57.

[5]XUBQ,GAOXL,DONGKJ,etal.Grainproteincontentcomparisonandproteomicanalysisoffoxtailmillet(Setariaitalica"L.)seedresponsetodifferentdroughtstresslevels[J].Actaphysiologiaeplantarum,2020,42(2):1-11.

[6]XIANGJL,ZHANGM,APEA-BAHFB,etal.Hydroxycinnamicacidamide(HCAA)derivatives,flavonoidC-glycosides,phenolicacidsandantioxidantpropertiesoffoxtailmillet[J].Foodchemistry,2019,295:214-223.

[7]SHARMAS,SAXENADC,RIARCS.Characteristicsofβ-glucanextractedfromrawandgerminatedfoxtail(Setariaitalica)andkodo(Paspalumscrobiculatum)millets[J].Internationaljournalofbiologicalmacromolecules,2018,118:141-148.

[8]劉建壘,常柳,段曉亮,等.谷子的生產(chǎn)概況及其保健功能與機(jī)理研究進(jìn)展[J].食品工業(yè)科技,20243(5):389-395.

[9]DEVIPB,VIJAYABHARATHIR,SATHYABAMAS,etal.Healthbenefitsoffingermillet(Eleusinecoracana"L.)polyphenolsanddietaryfiber:Areview[J].Journaloffoodscienceandtechnology,2014,51(6):1021-1040.

[10]HASSANZM,SEBOLANA,MABELEBELEM.Thenutritionaluseofmilletgrainforfoodandfeed:Areview[J].Agricultureamp;foodsecurity,20210(1):1-14.

[11]閆鋒,崔秀輝,李清泉,等.谷子主要農(nóng)藝性狀的遺傳參數(shù)分析[J].黑龍江農(nóng)業(yè)科學(xué),2010(3):28-30.

[12]劉正理,程汝宏,張鳳蓮,等.不同密度條件下3種類型谷子品種產(chǎn)量及其構(gòu)成要素變化特征研究[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào),2007,15(5):135-138.

[13]莊云,馬堯,牟金明.密度對(duì)谷子生長(zhǎng)及產(chǎn)量性狀的影響[J].安徽農(nóng)業(yè)科學(xué),2007,35(36):11795,11866.

[14]田伯紅,王建廣,李雅靜,等.空間誘變對(duì)谷子農(nóng)藝性狀效應(yīng)的研究[J].植物遺傳資源學(xué)報(bào),2008,9(3):340-345.

[15]田伯紅.鹽堿脅迫對(duì)谷子農(nóng)藝性狀及產(chǎn)量的影響[J].河北農(nóng)業(yè)科學(xué),2009,13(10):4-5.

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

[17]魏常敏,許衛(wèi)猛,邢永鋒,等.不同年份玉米雜交種主要農(nóng)藝性狀與產(chǎn)量的灰色關(guān)聯(lián)度分析[J].浙江農(nóng)業(yè)科學(xué),2019,60(10):1756-1759.

[18]吳榮華,莊克章,張春艷,等.魯南地區(qū)玉米產(chǎn)量與主要農(nóng)藝性狀的灰色關(guān)聯(lián)度分析[J].作物研究,2019,33(6):524-527.

[19]曹鵬鵬,王樂(lè)政.德州地區(qū)大豆優(yōu)異種質(zhì)資源的篩選[J].遼寧農(nóng)業(yè)科學(xué),2015(3):35-38.

[20]屈洋,王可珍,劉洋,等.大豆產(chǎn)量與主要農(nóng)藝性狀的灰色關(guān)聯(lián)度分析[J].陜西農(nóng)業(yè)科學(xué),2019,65(2):39-467.

[21]施立安,張艷軍,沈祥宏,等.玉溪市小麥新品系主要農(nóng)藝性狀與產(chǎn)量的灰色關(guān)聯(lián)度分析[J].農(nóng)業(yè)科技通訊,2020(8):77-79,293.

[22]代金英,張桂云,胡蕾,等.耐鹽水稻產(chǎn)量與主要農(nóng)藝性狀的灰色關(guān)聯(lián)度分析[J].大麥與谷類科學(xué),2020,37(6):9-13,20.

[23]趙禹凱,王顯瑞,張立媛,等.谷子產(chǎn)量與主要農(nóng)藝性狀的灰色關(guān)聯(lián)度分析[J].吉林農(nóng)業(yè)科學(xué),2014,39(2):9-12.

[24]欒素榮,王占廷,李青松.谷子產(chǎn)量與主要農(nóng)藝性狀的灰色關(guān)聯(lián)度分析[J].河北農(nóng)業(yè)科學(xué),2010,14(11):115-116,118.

[25]王立秋.玉米雜交種產(chǎn)量性狀與產(chǎn)量的灰色關(guān)聯(lián)度分析[J].玉米科學(xué),1997,5(4):23-25,29.

猜你喜歡
灰色關(guān)聯(lián)分析產(chǎn)量
2022年11月份我國(guó)鋅產(chǎn)量同比增長(zhǎng)2.9% 鉛產(chǎn)量同比增長(zhǎng)5.6%
提高玉米產(chǎn)量 膜下滴灌有效
世界致密油產(chǎn)量發(fā)展趨勢(shì)
海水稻產(chǎn)量測(cè)評(píng)平均產(chǎn)量逐年遞增
2018年我國(guó)主要水果產(chǎn)量按省(區(qū)、市)分布
2018上半年我國(guó)PVC產(chǎn)量數(shù)據(jù)
聚氯乙烯(2018年9期)2018-02-18 01:11:34
運(yùn)動(dòng)員組織承諾水平的評(píng)價(jià)與提升策略
新疆向西開放度與經(jīng)濟(jì)增長(zhǎng)灰色關(guān)聯(lián)分析
基于灰色關(guān)聯(lián)分析的制造企業(yè)跨國(guó)并購(gòu)財(cái)務(wù)決策
秸稈資源化對(duì)吉林省經(jīng)濟(jì)結(jié)構(gòu)的影響研究
主站蜘蛛池模板: 亚洲国产日韩在线观看| 又大又硬又爽免费视频| 色网站免费在线观看| 黄色网页在线观看| 尤物午夜福利视频| 丰满少妇αⅴ无码区| 亚洲AV无码乱码在线观看代蜜桃| 欧美在线网| 亚洲天堂首页| 亚洲一区色| 一级毛片免费播放视频| 亚洲精品日产精品乱码不卡| 岛国精品一区免费视频在线观看| 四虎AV麻豆| 国产三级视频网站| 国产精品露脸视频| 片在线无码观看| 国产乱子伦精品视频| 国产黄网站在线观看| 国产95在线 | 欧美一区国产| 青青草91视频| 亚洲男人天堂网址| 国产在线麻豆波多野结衣| 欧美日韩中文国产| 色偷偷av男人的天堂不卡| 超碰色了色| 中文字幕天无码久久精品视频免费| 一本大道香蕉高清久久| 国产成人综合日韩精品无码不卡| 国产精品视频系列专区| 亚洲天堂视频网站| 一级全黄毛片| 五月婷婷综合色| 亚洲欧美一区二区三区蜜芽| 在线观看无码av五月花| 青青草原国产| 国产三级视频网站| 欧美日韩资源| 国产丝袜啪啪| 婷婷激情亚洲| 国产一级毛片高清完整视频版| 91久草视频| 欧美精品伊人久久| 精品成人一区二区三区电影 | 欧美成人日韩| 亚洲国产91人成在线| 亚洲精品免费网站| 一本大道无码日韩精品影视 | 黄色a一级视频| 国产精品无码AV中文| 亚洲精选无码久久久| 97在线观看视频免费| 在线中文字幕日韩| 热re99久久精品国99热| 一本色道久久88| 欧美综合成人| 国产麻豆另类AV| 亚洲av无码人妻| 亚洲男人天堂2018| 国产三级精品三级在线观看| 99国产精品国产| 久久国产V一级毛多内射| 91破解版在线亚洲| www中文字幕在线观看| 国产亚洲精品va在线| 2021国产精品自产拍在线观看| 日韩中文无码av超清| 亚洲国产系列| 国产欧美网站| 男女性色大片免费网站| 日韩无码白| 91视频日本| 真人免费一级毛片一区二区| 亚洲妓女综合网995久久| 久久永久视频| 女同久久精品国产99国| av在线人妻熟妇| 成人午夜天| 免费在线播放毛片| 自慰网址在线观看| 亚洲色欲色欲www在线观看|