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

    貴州省糧食產(chǎn)量影響因素的偏最小二乘回歸分析

    2016-05-06 02:07:26李福奪楊興洪
    關(guān)鍵詞:糧食生產(chǎn)貴州省影響因素

    李福奪, 楊興洪

    (貴州大學(xué) 管理學(xué)院, 貴州 貴陽(yáng) 550025)

    ?

    貴州省糧食產(chǎn)量影響因素的偏最小二乘回歸分析

    李福奪,楊興洪*

    (貴州大學(xué) 管理學(xué)院, 貴州 貴陽(yáng) 550025)

    摘要:以貴州省為研究區(qū)域,運(yùn)用偏最小二乘回歸分析方法,構(gòu)建貴州省糧食生產(chǎn)影響因素的偏最小二乘回歸模型.通過(guò)對(duì)模型進(jìn)行檢驗(yàn),發(fā)現(xiàn)該回歸模型對(duì)自變量X和因變量Y的解釋能力較強(qiáng)自變量集與因變量集之間存在著顯著的線性關(guān)系,且擬合優(yōu)度高達(dá)R2=0.988 6,表明回歸模型精度較高,可靠性強(qiáng).研究結(jié)果表明,在所選取的指標(biāo)中,糧食播種面積、有效灌溉面積、農(nóng)業(yè)機(jī)械總動(dòng)力、化肥施用量和成災(zāi)面積等是影響貴州省糧食生產(chǎn)的關(guān)鍵因素.

    關(guān)鍵詞:貴州省; 糧食生產(chǎn); 影響因素; 偏最小二乘回歸法

    世界氣候變化、人口增長(zhǎng)及環(huán)境問(wèn)題的日益凸顯,糧食生產(chǎn)與安全問(wèn)題受到世界各國(guó)的關(guān)注,這方面的研究也越來(lái)越多,越來(lái)越深入.國(guó)內(nèi)圍繞糧食產(chǎn)量影響因素問(wèn)題,趙慧江[1]在1990—2005年統(tǒng)計(jì)數(shù)據(jù)的基礎(chǔ)上,采用回歸分析方法通過(guò)建立以5種可量化的影響因素為自變量,以糧食產(chǎn)量為因變量的多元線性回歸模型,對(duì)影響我國(guó)糧食產(chǎn)量主要因子進(jìn)行了分析,并以此為依據(jù)提出了關(guān)于我國(guó)穩(wěn)定發(fā)展糧食生產(chǎn)的參考意見(jiàn).肖海峰[2]通過(guò)建立糧食總產(chǎn)量和影響因素的柯布-道格拉斯生產(chǎn)函數(shù),對(duì)我國(guó)1978—2002年的數(shù)據(jù)進(jìn)行了分段研究,來(lái)分析我國(guó)糧食綜合生產(chǎn)能力.高倩倩等[3]運(yùn)用統(tǒng)計(jì)學(xué)的知識(shí)、采用Eviews 6.0統(tǒng)計(jì)軟件建立了糧食產(chǎn)量評(píng)估模型,對(duì)建國(guó)以來(lái)山東省糧食增產(chǎn)的主要因素進(jìn)行了分析,得出了提高糧食單產(chǎn)是糧食增產(chǎn)的最有效途徑的結(jié)論.杜國(guó)明等[4]以黑龍江省1981—2011年糧食總產(chǎn)、單產(chǎn)等方面的統(tǒng)計(jì)資料為數(shù)據(jù)源,運(yùn)用數(shù)理分析的方法對(duì)該省糧食產(chǎn)量的影響因素進(jìn)行了比較分析,認(rèn)為糧食增產(chǎn)主要依靠糧食播種面積等3個(gè)方面的提升.廖婧琳等[5]則更加關(guān)注自然災(zāi)害對(duì)糧食生產(chǎn)的影響,他們分析了水旱、冰雹、病蟲(chóng)害和低溫冷凍等農(nóng)業(yè)災(zāi)害對(duì)貴州省糧食產(chǎn)量的影響機(jī)理,并提出了相應(yīng)的應(yīng)對(duì)政策.盡管目前學(xué)術(shù)界對(duì)糧食產(chǎn)量影響因素的分析方法很多,但幾乎都存在著樣本量少,前提假設(shè)可靠性難以確認(rèn),模型難以充分包絡(luò)影響因子的有效信息、容易出現(xiàn)“偽回歸”現(xiàn)象等問(wèn)題.為了解決上述問(wèn)題,為糧食增產(chǎn)影響因素的分析提供一種可靠的計(jì)量模型,本文提出了一種新的統(tǒng)計(jì)回歸分析方法——偏最小二乘回歸分析( PLS),并以貴州省為研究區(qū)域,構(gòu)建偏最小二乘回歸分析模型,對(duì)影響貴州省糧食產(chǎn)量的主要因素進(jìn)行了定量分析,以期為貴州省糧食生產(chǎn)與安全戰(zhàn)略的實(shí)施提供科學(xué)的決策依據(jù).

    1研究區(qū)域與研究方法

    1.1研究區(qū)域概況貴州省地形以山地為主,山地、高原占全省土地面積的87%,丘陵地占10%,盆地、河流灘涂等其他地貌占3%,是我國(guó)唯一沒(méi)有平原地貌的省份,境內(nèi)地勢(shì)西高東低,自中部向南、北、東三面傾斜,平均海拔在1 100 m左右.這種地理特點(diǎn),使得貴州省可用于農(nóng)業(yè)開(kāi)發(fā)的土地資源不多,耕地面積少、質(zhì)量差,開(kāi)發(fā)利用比例大,后備耕地嚴(yán)重不足,有限的耕地資源給貴州省的糧食生產(chǎn)與安全造成了極大的壓力,也嚴(yán)重制約著貴州省經(jīng)濟(jì)的持續(xù)發(fā)展.2005年貴州省委省政府一號(hào)文件指出:“當(dāng)前和今后一個(gè)時(shí)期,必須堅(jiān)定不移地把加強(qiáng)農(nóng)業(yè)基礎(chǔ)設(shè)施建設(shè),加快農(nóng)業(yè)科技進(jìn)步,提高農(nóng)業(yè)綜合生產(chǎn)能力作為一項(xiàng)緊迫的戰(zhàn)略任務(wù),切實(shí)抓緊抓好,以確保我省糧食產(chǎn)量穩(wěn)定增長(zhǎng)和安全消費(fèi).”[6]實(shí)現(xiàn)糧食穩(wěn)定生產(chǎn)、增強(qiáng)糧食自給能力是貴州省三農(nóng)工作的重點(diǎn),也是必須堅(jiān)持的重大原則.

    1.2數(shù)據(jù)來(lái)源與指標(biāo)體系現(xiàn)實(shí)中影響糧食產(chǎn)量的因素有很多,對(duì)這些因素進(jìn)行量化是進(jìn)行分析研究的基礎(chǔ),但這同時(shí)也是一個(gè)比較復(fù)雜和困難的工作.一些學(xué)者利用灰色系統(tǒng)理論,對(duì)影響糧食產(chǎn)量的因素進(jìn)行了動(dòng)態(tài)灰色關(guān)聯(lián)量化處理;還有學(xué)者通過(guò)建立糧食生產(chǎn)的時(shí)空回歸模型,對(duì)糧食產(chǎn)量的影響因素在時(shí)間和空間二維坐標(biāo)上進(jìn)行了量化分析[7].本文從數(shù)據(jù)可獲取性、指標(biāo)量化的可行性、主導(dǎo)因素重點(diǎn)凸顯性等指標(biāo)選取的原則出發(fā),參考前人研究成果,并結(jié)合貴州省農(nóng)業(yè)系統(tǒng)的實(shí)際狀況,選取糧食總產(chǎn)量Y作為要預(yù)測(cè)的目標(biāo)數(shù)據(jù),選取糧食播種面積x1、有效灌溉面積x2、農(nóng)業(yè)生產(chǎn)資料價(jià)格指數(shù)x3、農(nóng)業(yè)機(jī)械總動(dòng)力x4、農(nóng)業(yè)從業(yè)人員x5、家庭農(nóng)業(yè)支出x6、農(nóng)村用電量x7、化肥施用量x8、政府財(cái)政支農(nóng)支出x9和成災(zāi)面積x1010個(gè)對(duì)貴州省糧食產(chǎn)量有較大影響的指標(biāo)來(lái)構(gòu)建農(nóng)業(yè)系統(tǒng)糧食生產(chǎn)的指標(biāo)體系.其中,x1~x4、x7、x8數(shù)據(jù)來(lái)源于1992—2012年的《貴州統(tǒng)計(jì)年鑒》和《新中國(guó)60年統(tǒng)計(jì)資料匯編》,x5、x6、x9和x104個(gè)指標(biāo)數(shù)據(jù)由《貴州統(tǒng)計(jì)年鑒》(1992—2012)、《中國(guó)統(tǒng)計(jì)年鑒》和《貴州年鑒》等國(guó)家正式發(fā)布的統(tǒng)計(jì)資料整理計(jì)算而來(lái).

    1.3研究方法偏最小二乘回歸分析是用來(lái)解決變量多重相關(guān)性對(duì)系統(tǒng)回歸建模干擾問(wèn)題的一種方法,在處理樣本容量小、解釋變量個(gè)數(shù)多、變量間存在嚴(yán)重多重相關(guān)性問(wèn)題方面具有獨(dú)特的優(yōu)勢(shì)[8]. 偏最小二乘回歸分析可以實(shí)現(xiàn)回歸建模、數(shù)據(jù)結(jié)構(gòu)簡(jiǎn)化以及2組變量間的相關(guān)性分析.在決定是否采用PLS方法建模前,首先要進(jìn)行預(yù)備分析,以判斷自變量(因變量)是否存在多重相關(guān)性,判斷自變量和因變量間是否存在相關(guān)關(guān)系,具體計(jì)算方法是:對(duì)于矩陣Z=(X,Y),記

    當(dāng)因變量Y的階數(shù)為1時(shí),為單變量偏最小二乘回歸分析模型(PLS1).偏最小二乘回歸分析模型的構(gòu)建步驟如下:

    1) 標(biāo)準(zhǔn)化原始數(shù)據(jù).標(biāo)準(zhǔn)化后的數(shù)據(jù)矩陣E0=(fij)n×q(即X) 和F0=(eij)n×p(即Y),其中,

    (1)

    (2)

    (3)

    其中W1是E0的第一個(gè)軸,稱(chēng)為模型的效應(yīng)權(quán)重,且等于1.

    (4)

    其中C1是F0的第一個(gè)軸,稱(chēng)為模型的因變量權(quán)重,且等于1. 這里要求t1和u1能分別較好的表達(dá)X和Y中的數(shù)據(jù)信息,且t1對(duì)u1具有較強(qiáng)的解釋能力. 之后求殘差矩陣E1和F1:

    其中

    從而得到F0、E0對(duì)t1的回歸方程:

    (5)

    (6)

    3) 第2成分t2的提取.令

    回到以上第3)步,對(duì)殘差矩陣進(jìn)行新一輪的成分提取和回歸分析,重復(fù)第1成分t1提取步驟,最終得到F1、E1對(duì)t2的回歸方程:

    (7)

    (8)

    第h成分th的提取同理.

    4) 構(gòu)建偏最小二乘回歸模型.根據(jù)以上分析,構(gòu)建的偏最小二乘回歸模型如下:

    (9)

    其中

    F2為殘差矩陣.

    (10)

    PRESS越小說(shuō)明模型的擬合效果越好,PRESS取最小值時(shí),模型的擬合效果達(dá)到最優(yōu),這時(shí)提取的成分個(gè)數(shù)h即為最佳成分個(gè)數(shù).

    2結(jié)果與分析

    2.1研究結(jié)果利用SIMCA-P11.5軟件構(gòu)建偏最小二乘回歸模型,對(duì)模型進(jìn)行自動(dòng)擬合,相關(guān)結(jié)果見(jiàn)表1.

    表 1 被偏最小二乘因子解釋的變差的百分比

    從表1可以看出,從自變量中提取的5個(gè)成分可解釋變差的百分比分別是75.331 2%、15.013 4%、4.550 3%、3.275 1%和2.447 9%;計(jì)算可得,t1和u1的PRESS=0.663 4,t2和u2的PRESS=0.528 6,t3和u3的PRESS=0.374 5,t4和u4的PRESS=0.155 8,t5和u5的PRESS=-0.201 5,根據(jù)PRESS最小原則,最終確定應(yīng)該提取4個(gè)PLS成分.組合模型為:

    t1=0.331 27x1+0.326 55x2+0.333 16x3+

    0.317 62x4+0.291 95x5+0.247 81x6+0.309 97x7+

    0.350 55x8+0.267 82x9-0.291 69x10, u1=0.303 68Y;

    t2=0.320 06x1-0.291 25x2+0.287 51x3+

    0.321 30x4+0.280 05x5-0.265 75x6+0.321 22x7+

    0.340 79x8+0.265 43x9-0.240 07x10, u2=0.310 03Y;

    t3=0.255 35x1+0.301 38x2+0.191 98x3-

    0.331 29x4+0.290 15x5+0.301 23x6-0.321 45x7+

    0.320 09x8+0.223 34x9-0.313 36x10, u3=0.307 65Y;

    t4=0.341 21x1+0.240 07x2-0.231 23x3+

    0.314 55x4-0.292 12x5+0.314 56x6+0.276 57x7+

    0.334 40x8+0.199 89x9-0.285 58x10,

    u4=0.291 24Y.

    根據(jù)以上組合模型,可得回歸模型的檢驗(yàn)統(tǒng)計(jì)量,其調(diào)整的R2值、F檢驗(yàn)值和t檢驗(yàn)值,見(jiàn)表2.

    表 2 回歸模型檢驗(yàn)統(tǒng)計(jì)量

    從表2可以看出:

    1) 擬合優(yōu)度R2=0.988 6,即t1、t2、t3和t4解釋了因變量98.86%的信息;

    2) 方程顯著性檢驗(yàn)(F檢驗(yàn)):F=70.287 33,通過(guò)查閱F分布表得F0.05(3.30)=2.92,又因?yàn)镕>F0.05(3.30),因此在99%的置信度下通過(guò)了F檢驗(yàn),所以因變量Y與t1、t2、t3、t4存在較好的線性關(guān)系如圖1所示;

    3) 解釋變量t檢驗(yàn):由于t0.05(30)=2.092,經(jīng)回歸分析所得t值大于這一數(shù)值,因此所有解釋變量都通過(guò)了t檢驗(yàn).

    4) 序列相關(guān)性杜賓·瓦森(Durbin-Watson)檢驗(yàn):由SIMCA-P11.5軟件構(gòu)建的偏最小二乘回歸模型的回歸結(jié)果可知,該模型的Durbin-Watson 檢驗(yàn)值為1.528 431,查表得

    dL=1.19,dU=1.47,

    故有

    dU

    所以x1、x2、…、x10序列不存在1階自相關(guān)性.通過(guò)以上分析得出,模型通過(guò)了可靠性檢驗(yàn),且對(duì)自變量X和因變量Y解釋能力是很強(qiáng)的,分別達(dá)到了97.45%和98.86%,因此回歸模型精度較高,可靠性也是比較強(qiáng)的.由此最終得到的標(biāo)準(zhǔn)化偏最小二乘回歸方程為:

    Y=0.193 25x1+0.183 76x2+0.020 17x3+

    0.174 43x4+0.052 01x5+0.031 52x6+0.019 95x7+

    0.160 44x8+0.017 879x9-0.184 03x10+1.776 25 ,

    2.2結(jié)果分析結(jié)果表明,在上述選取的10個(gè)對(duì)貴州省糧食生產(chǎn)有較大影響的因素中,糧食播種面積、有效灌溉面積、農(nóng)業(yè)機(jī)械總動(dòng)力、化肥施用量和成災(zāi)面積等是影響貴州省糧食生產(chǎn)的關(guān)鍵因素.其中,糧食播種面積每變動(dòng)1個(gè)單位,糧食總產(chǎn)量就會(huì)變動(dòng)0.193 25個(gè)單位;有效灌溉面積變動(dòng)1單位,糧食總產(chǎn)量則會(huì)變動(dòng)0.183 76個(gè)單位;農(nóng)業(yè)機(jī)械總動(dòng)力每提升1單位,就會(huì)引起糧食總產(chǎn)量以0.174 43的倍數(shù)增加;化肥施用量每增加1單位,會(huì)使糧食總產(chǎn)量增加0.160 44單位;而成災(zāi)面積每減少1單位,就會(huì)促使糧食產(chǎn)量增加0.184 03倍.圖2和表3分別為貴州省糧食總產(chǎn)實(shí)際值與模型擬合值的對(duì)比曲線和貴州省1992—2012年影響糧食生產(chǎn)的5個(gè)關(guān)鍵指標(biāo)的相關(guān)統(tǒng)計(jì)數(shù)據(jù),由圖2可以看出擬合值與實(shí)際值具有很高的擬合度,最大擬合誤差僅5.5%.

    表 3 貴州省1992—2012年影響糧食生產(chǎn)關(guān)鍵指標(biāo)統(tǒng)計(jì)數(shù)據(jù)

    2.2.1糧食播種面積由圖3可以看出,貴州省糧食總產(chǎn)量和糧食播種面積之間呈現(xiàn)出較強(qiáng)的正相關(guān)關(guān)系,由貴州省糧食產(chǎn)量影響因素的偏最小二乘回歸方程可知相關(guān)系數(shù)為0.193 25,是對(duì)該省糧食產(chǎn)量影響最大的一個(gè)因素.貴州省糧食播種面積從1992年以來(lái)總體上呈現(xiàn)出逐年遞增趨勢(shì),但增加幅度很小,增加速度也十分緩慢,這與貴州省多山多丘陵少平地的地理?xiàng)l件有關(guān).貴州省可用于農(nóng)業(yè)開(kāi)發(fā)的土地資源不多,耕地面積少、質(zhì)量差,開(kāi)發(fā)利用比例大,后備耕地嚴(yán)重不足,這極大地制約了糧食播種面積的擴(kuò)大和整個(gè)農(nóng)業(yè)經(jīng)濟(jì)的發(fā)展.同時(shí),由于貴州省農(nóng)村勞動(dòng)力大量向城市轉(zhuǎn)移,導(dǎo)致農(nóng)業(yè)勞動(dòng)資源外流,使得農(nóng)業(yè)生產(chǎn)出現(xiàn)勞動(dòng)力短缺的現(xiàn)象,這也是制約糧食播種面積進(jìn)一步增加的原因[9].2000年以后,隨著貴州省農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)戰(zhàn)略性調(diào)整,加之糧食比較利益的下降,貴州省糧食播種面積開(kāi)始減少,相應(yīng)的糧食總產(chǎn)量也出現(xiàn)了自1992年8連增以來(lái)的首次下滑,如圖3所示.特別是2006年,貴州省全省糧食播種面積僅為389 萬(wàn)hm2,比2000年減少了89萬(wàn)hm2,糧食產(chǎn)量減少36萬(wàn)t,6年間糧食播種面積減少比例達(dá)18.62%,年均減少3.1%,相應(yīng)的糧食產(chǎn)量同比下降4個(gè)百分點(diǎn).正由于貴州省糧食總產(chǎn)量和糧食播種面積之間存在的這種關(guān)系,在目前的農(nóng)業(yè)技術(shù)水平下,貴州省要增加糧食產(chǎn)量確保糧食安全,就必須保證一定的糧食播種面積來(lái)作為基本保障.

    2.2.2有效灌溉面積貴州省目前的農(nóng)田水利設(shè)施普遍存在著維護(hù)不足、老化嚴(yán)重,建設(shè)速度趕不上折舊速度等問(wèn)題,導(dǎo)致農(nóng)田有效灌溉面積難以滿(mǎn)足提高糧食產(chǎn)量的需要[10].由圖4可以看出,貴州省農(nóng)田有效灌溉面積在1992—2006這10多年間增長(zhǎng)速度經(jīng)歷了一個(gè)緩慢的波動(dòng),有效灌溉水平并沒(méi)有得到實(shí)質(zhì)性的提高,只是到了2006年以后這種情況才得到初步改善.

    據(jù)貴州省水利廳公布的資料顯示,2012年貴州全省農(nóng)田有效灌溉面積為154.84 萬(wàn)hm2,其中稻田92.15 萬(wàn)hm2,占貴州省稻田總面積的42.6%;旱地62.69 萬(wàn)hm2,占全省旱地總面積的27.83%,全省仍有57.4%的稻田和72.17%的旱地要依靠自然降雨和提灌等方式來(lái)保證農(nóng)業(yè)生產(chǎn)用水.但是,近些年來(lái),隨著貴州省開(kāi)始重視生態(tài)治理,作為配套工程修建了一些現(xiàn)代化的農(nóng)業(yè)灌溉基礎(chǔ)設(shè)施,再加之國(guó)家西部扶貧開(kāi)發(fā)資金源源不斷的引入,加強(qiáng)了貴州省農(nóng)業(yè)綜合開(kāi)發(fā)力度,一批節(jié)水灌溉水利設(shè)施也逐步完成建設(shè)并投入使用,這些都增加了農(nóng)田有效灌溉面積、提高了有效灌溉比例[11],由圖4和表3可看出,2006年到2012年貴州省的農(nóng)田有效灌溉面積從109.38萬(wàn)hm2提高到154.84萬(wàn)hm2,6年間增加了45.46萬(wàn)hm2,增加比例達(dá)41.56%,原因正在于此.雖然近年來(lái)貴州省農(nóng)田有效灌溉面積開(kāi)始有所增加,但要達(dá)到相當(dāng)規(guī)模以滿(mǎn)意糧食生產(chǎn)的需求,仍將需要相當(dāng)長(zhǎng)的一段時(shí)間.

    2.2.3農(nóng)業(yè)機(jī)械總動(dòng)力由貴州省糧食產(chǎn)量影響因素的偏最小二乘回歸方程可以看出,農(nóng)業(yè)機(jī)械總動(dòng)力對(duì)貴州省糧食生產(chǎn)的貢獻(xiàn)率為0.174 43,也就是說(shuō)農(nóng)業(yè)機(jī)械總動(dòng)力每增加1個(gè)單位,就會(huì)引起糧食產(chǎn)量以0.174 43的倍數(shù)增加,可見(jiàn)農(nóng)業(yè)機(jī)械總動(dòng)力也是貴州省糧食增產(chǎn)的關(guān)鍵因素之一.1992—2012年這20年間,貴州省農(nóng)業(yè)機(jī)械總動(dòng)力從0.8 萬(wàn)臺(tái)提升到11.3萬(wàn)臺(tái),同比增長(zhǎng)131.25%,其中50 馬力以上的大中型機(jī)械增加5.6 萬(wàn)臺(tái),占到全部增加量的44.8%,尤其是2008年以后,貴州省的農(nóng)業(yè)機(jī)械總動(dòng)力呈現(xiàn)出直線型增長(zhǎng)趨勢(shì),短短4年就增加了7.9 萬(wàn)臺(tái),增幅達(dá)232.35%,農(nóng)業(yè)機(jī)械總動(dòng)力的提升對(duì)促進(jìn)貴州省糧食產(chǎn)量的增加起到了重要的作用.然而,必須指出的是由于貴州省多山多坡多丘陵的地理狀況,耕地分散且種植條件普遍較差,導(dǎo)致農(nóng)業(yè)機(jī)械化進(jìn)程受到很大制約,農(nóng)業(yè)機(jī)械總動(dòng)力雖然在不斷增加,但發(fā)展瓶頸也越來(lái)越大,對(duì)糧食產(chǎn)量的促進(jìn)作用減少,相關(guān)問(wèn)題值得進(jìn)一步探討.

    2.2.4化肥施用量由偏最小二乘回歸擬合方程可知,化肥施用量對(duì)貴州省糧食生產(chǎn)的貢獻(xiàn)率為0.160 44,可見(jiàn)適當(dāng)增加化肥施用量也是促進(jìn)糧食增產(chǎn)的有效途徑之一.從下圖6可以看出,自1992年起,貴州省農(nóng)業(yè)化肥施用量總體上處在不斷增加的狀態(tài),除2001—2006年出現(xiàn)局部時(shí)間段小幅下降外,其他年份增幅都比較平穩(wěn),到2012年貴州全省農(nóng)業(yè)化肥施用量達(dá)到156.36 萬(wàn)t,是1992年的1.79倍.貴州省農(nóng)業(yè)化肥施用量與糧食產(chǎn)量的關(guān)聯(lián)度較高,對(duì)糧食產(chǎn)量影響也比較大,化肥在農(nóng)業(yè)生產(chǎn)中的廣泛應(yīng)用,對(duì)糧食穩(wěn)產(chǎn)和增產(chǎn)都起到了比較大的作用.然而另一方面要指出的是,由于受邊際效益遞減規(guī)律的影響,化肥投入在貴州省糧食生產(chǎn)方面的促進(jìn)作用也必將逐步下降,因此一味地靠增加化肥施用量來(lái)增加糧食產(chǎn)量的做法并不可行,關(guān)鍵是要提高化肥的使用效率.

    2.2.5成災(zāi)面積由圖7可以看出,貴州省糧食總產(chǎn)量和自然災(zāi)害成災(zāi)面積之間呈現(xiàn)出一定的負(fù)相關(guān)關(guān)系,相關(guān)系數(shù)為0.184 03,是對(duì)該省糧食產(chǎn)量影響第2大的因素.因受農(nóng)業(yè)自然災(zāi)害的影響, 貴州省糧食產(chǎn)量多年來(lái)波動(dòng)較大.貴州省常發(fā)性農(nóng)業(yè)自然災(zāi)害有冰雹、春旱、暴雨、倒春寒和作物病蟲(chóng)害等,其中寒災(zāi)和旱災(zāi)是對(duì)該省糧食生產(chǎn)影響最大的2種自然災(zāi)害[12].貴州省農(nóng)業(yè)自然災(zāi)害成災(zāi)面積在1992—2012年這20年里經(jīng)歷了一個(gè)“減-增-減-增-減”的波動(dòng)變化過(guò)程,1992年,貴州省農(nóng)業(yè)自然災(zāi)害成災(zāi)面積為56.36 hm2,約占當(dāng)年全省糧食播種面積的15%,之后的7年里,貴州省糧食成災(zāi)面積逐年減少,到1999年,成災(zāi)面積僅為20.54 hm2,減少量為35.82 hm2,降幅達(dá)63.56%,是8年里的最小值.

    農(nóng)業(yè)自然災(zāi)害成災(zāi)面積的減少是糧食穩(wěn)產(chǎn)和增產(chǎn)的有力保障,從表3可以看出,1992—1999年,貴州省的糧食產(chǎn)量是逐年遞增的,從772 萬(wàn)t增加到880 萬(wàn)t,增幅達(dá)14%.1999年以后,貴州省又經(jīng)歷了2個(gè)自然災(zāi)害多發(fā)期(分別為1999—2002年和2005—2006年)和2個(gè)自然災(zāi)害少發(fā)期(分別為2002—2005年和2006—2012年),在2個(gè)自然災(zāi)害多發(fā)期內(nèi)和2個(gè)自然災(zāi)害少發(fā)期內(nèi),貴州省糧食作物年均成災(zāi)面積分別為56.10 萬(wàn)hm2和21.34 萬(wàn)hm2,糧食產(chǎn)量年均增長(zhǎng)7.6%和15.8%,可見(jiàn)自然災(zāi)害成災(zāi)面積對(duì)貴州省糧食生產(chǎn)的影響程度是很大的.

    3結(jié)論與建議

    3.1結(jié)論1) 利用SIMCA-P11.5軟件構(gòu)建的偏最小二乘回歸模型對(duì)X和Y的解釋能力較強(qiáng),通過(guò)對(duì)自變量集和因變量的相關(guān)關(guān)系進(jìn)行分析,發(fā)現(xiàn)兩者之間存在著明顯的線性關(guān)系,通過(guò)對(duì)模型進(jìn)行自動(dòng)擬合,然后把擬合值與實(shí)際值進(jìn)行比較,發(fā)現(xiàn)最大誤差只有5.5%,因此認(rèn)為擬合效果很好,構(gòu)建的模型具有可靠性,在貴州省糧食生產(chǎn)的影響因素的相關(guān)性方面具有良好的分析能力.

    2) 貴州省糧食產(chǎn)量影響因素的偏最小二乘回歸模型表明,糧食播種面積、有效灌溉面積、農(nóng)業(yè)機(jī)械總動(dòng)力、化肥施用量和成災(zāi)面積5個(gè)指標(biāo)是影響貴州省糧食生產(chǎn)的關(guān)鍵因素.其中,前4個(gè)指標(biāo)與糧食產(chǎn)量正相關(guān),第5個(gè)指標(biāo)與糧食產(chǎn)量負(fù)相關(guān);糧食播種面積對(duì)糧食產(chǎn)量影響最大,化肥施用量對(duì)糧食產(chǎn)量影響相對(duì)在減弱.

    3) 由于貴州省特殊的地理地貌和氣候環(huán)境制約,在農(nóng)業(yè)現(xiàn)代化和機(jī)械化發(fā)展初期,提高農(nóng)業(yè)生產(chǎn)機(jī)械化水平可以提高生產(chǎn)力,但當(dāng)機(jī)械化水平達(dá)到一定程度后其發(fā)展便會(huì)遇到瓶頸;在一定的時(shí)期和一定的水平下,增加化肥施用量會(huì)提高糧食產(chǎn)量,但受邊際效益遞減規(guī)律的影響,在目前貴州省化肥施用量已很大的前提下,繼續(xù)增加化肥施用量不一定能夠繼續(xù)增加糧食產(chǎn)量.

    4) 成災(zāi)面積對(duì)貴州省糧食生產(chǎn)的影響較大,因此政府應(yīng)該加大投入,做好防災(zāi)減災(zāi)預(yù)案,提高糧食生產(chǎn)的可能性.

    3.2建議

    3.2.1加強(qiáng)耕地保護(hù),穩(wěn)定糧食播種面積2012年貴州省政府1號(hào)文件要求要保證全省糧食播種面積穩(wěn)定在423.3 萬(wàn)hm2、糧食總產(chǎn)量穩(wěn)定在1 150 萬(wàn)t左右,為此就必須進(jìn)一步挖掘耕地潛力.近5年來(lái),貴州省糧食播種面積一直保持在441~499 萬(wàn)hm2之間,平均為475.8 萬(wàn)hm2,變異系數(shù)為5.72%,總體糧食播種面積還算比較穩(wěn)定.糧食播種面積會(huì)受到諸如糧價(jià)、氣候環(huán)境、耕地?cái)?shù)量與質(zhì)量等多方面因素的影響[13],只有確保糧食價(jià)格穩(wěn)定增加農(nóng)民種糧收入、加強(qiáng)預(yù)測(cè)預(yù)警以減小氣候環(huán)境對(duì)播種的影響、實(shí)行嚴(yán)格的耕地保護(hù)政策提高耕地質(zhì)量,才能進(jìn)一步挖掘挖掘農(nóng)地潛力,糧食播種面積才會(huì)有所保障.

    3.2.2加強(qiáng)農(nóng)田基礎(chǔ)設(shè)施建設(shè),提升有效灌溉水平貴州省目前的農(nóng)田水利基礎(chǔ)設(shè)施建設(shè)比較薄弱,普遍存在著維護(hù)不足、老化嚴(yán)重、建設(shè)速度趕不上折舊速度的問(wèn)題,導(dǎo)致農(nóng)田有效灌溉面積難以滿(mǎn)足提高糧食產(chǎn)量的需要.但有效灌溉面積對(duì)貴州省糧食產(chǎn)量影響較大,因此必須要加強(qiáng)農(nóng)田水利設(shè)施建設(shè),改善農(nóng)業(yè)灌溉條件,提升農(nóng)田有效灌溉水平,增加有效灌溉面積,才能增加糧食產(chǎn)量.

    3.2.3加強(qiáng)農(nóng)業(yè)災(zāi)害預(yù)警預(yù)報(bào),降低災(zāi)害損失要加快完善自然災(zāi)害預(yù)警預(yù)報(bào)系統(tǒng),最大限度的攻克突發(fā)性農(nóng)業(yè)自然災(zāi)害在預(yù)報(bào)方面的難關(guān).可以考慮將地理信息系統(tǒng)中的遙感遙測(cè)技術(shù)引入到自然災(zāi)害的監(jiān)測(cè)中,提高預(yù)報(bào)精度.要加快貴州省防汛抗旱領(lǐng)導(dǎo)辦公室的職能建設(shè),對(duì)災(zāi)后農(nóng)業(yè)減災(zāi)工作開(kāi)展專(zhuān)業(yè)指導(dǎo)、對(duì)口援助,最大限度的降低糧食產(chǎn)量損失.另一方面,還要深化貴州省綜合防治體系建設(shè),農(nóng)業(yè)自然災(zāi)害是一個(gè)復(fù)雜的系統(tǒng)工程,只有有的放矢、全面規(guī)劃、綜合防治,才能有效減少自然災(zāi)害對(duì)糧食產(chǎn)量的影響.

    3.2.4加快農(nóng)業(yè)科技創(chuàng)新,提高單糧食單產(chǎn)在可用耕地面積有限、糧食播種面積也不可能大幅增加的情況下,貴州省糧食總產(chǎn)量的增加就必須要依靠提高單位面積產(chǎn)量.在2006—2012年這6年間,貴州省糧食單產(chǎn)平均值為268 kg,變異系數(shù)為4.37%,而同期全國(guó)糧食單產(chǎn)平均值平均為312 kg,是貴州省的1.16倍,由此可以看出,貴州省糧食單位面積產(chǎn)量依舊很低,農(nóng)業(yè)生產(chǎn)效率不高.要想提高糧食單產(chǎn),就必須加快農(nóng)業(yè)科技創(chuàng)新,提高糧食生產(chǎn)中的科技含量,具體可以考慮選育適合貴州省生長(zhǎng)的優(yōu)良糧食作物品種,提高生物技術(shù)應(yīng)用能力、擴(kuò)大生物技術(shù)應(yīng)用范圍,研制不會(huì)破壞土壤結(jié)構(gòu)的綠色新型肥料提高土壤肥力等[14].

    3.2.5改革現(xiàn)有耕地制度,進(jìn)一步挖掘糧食生產(chǎn)潛力目前貴州省政府和各下級(jí)地方政府在耕地保護(hù)上目標(biāo)不一致:省政府更注重耕地的可持續(xù)利用和全省的糧食生產(chǎn)與安全,地方政府更注重短期經(jīng)濟(jì)效益而不太關(guān)注對(duì)耕地的保護(hù)和永續(xù)利用;貴州全省耕地規(guī)模化經(jīng)營(yíng)水平很低,不利于提高糧食生產(chǎn)效率.針對(duì)這些個(gè)問(wèn)題,必須果斷改革現(xiàn)有耕地制度,可考慮通過(guò)把耕地有效保護(hù)和地方政府政績(jī)掛鉤來(lái)規(guī)范其非理性的供地行為;要提高土地規(guī)模經(jīng)營(yíng)程度,通過(guò)立法立規(guī)為其提供法律保障,以進(jìn)一步挖掘糧食生產(chǎn)潛力,提高糧食產(chǎn)出效益[15].

    致謝貴州大學(xué)人文社科重點(diǎn)特色學(xué)科重大項(xiàng)目(GDZT12007)對(duì)本文給予了資助,謹(jǐn)致謝意.

    參考文獻(xiàn)

    [1] 趙慧江. 基于回歸分析的糧食產(chǎn)量影響因素分析[J]. 懷化學(xué)院學(xué)報(bào),2009,28(2):31-35.

    [2] 肖海峰. 我國(guó)糧食綜合生產(chǎn)能力影響因素分析[J]. 農(nóng)業(yè)技術(shù)經(jīng)濟(jì),2004(6):45-49.

    [3] 高倩倩,邢秀峰,姚傳進(jìn). 基于逐步回歸分析的糧食產(chǎn)量影響因素研究[J]. 當(dāng)代經(jīng)濟(jì),2010(9):145-147.

    [4] 杜國(guó)明,劉彥隨,劉閣,等. 黑龍江省近30年來(lái)糧食生產(chǎn)變化及增產(chǎn)因素分析[J]. 農(nóng)業(yè)現(xiàn)代化研究,2014,35(5):519-524.

    [5] 廖婧琳,蘇月,馮澤蔚,等. 54年來(lái)農(nóng)業(yè)自然災(zāi)害對(duì)貴州省農(nóng)業(yè)經(jīng)濟(jì)的影響分析[J]. 安徽農(nóng)業(yè)科學(xué),2008,36(25):114-117.

    [6] 趙子瑩,謝頌磊. 1978—2008年貴州省糧食產(chǎn)量影響因素分析[J]. 農(nóng)業(yè)現(xiàn)代科技,2013(8):336-337.

    [7] 謝杰. 中國(guó)糧食生產(chǎn)影響因素研究[J]. 經(jīng)濟(jì)問(wèn)題探索,2007(9):36-40.

    [8] 馬明德,馬學(xué)娟. 農(nóng)業(yè)經(jīng)濟(jì)增長(zhǎng)影響因子的偏最小二乘回歸分析研究[J]. 中國(guó)農(nóng)業(yè)資源與區(qū)劃,2014,35(5):123-128.

    [9] 肖厚軍. 貴州糧食綜合生產(chǎn)能力的影響因素與對(duì)策分析[J]. 貴州農(nóng)業(yè)科學(xué),2005,33(4):101-103.

    [10] 王天生,王瑤. 貴州省糧食生產(chǎn)發(fā)展60年回顧與展望[J]. 貴州農(nóng)業(yè)科學(xué),2009,37(9):1-9.

    [11] 楊靜,張和喜,蔡長(zhǎng)舉,等. 貴州省節(jié)水灌溉發(fā)展現(xiàn)狀及對(duì)策[J]. 中國(guó)農(nóng)村水利水電,2012(7):29-30.

    [12] 蔣彥. 貴州省自然災(zāi)害對(duì)農(nóng)業(yè)經(jīng)濟(jì)的影響分析[J]. 北京農(nóng)業(yè),2012(18):223.

    [13] 郭柏林. 我國(guó)糧食播種面積變化的地理研究[J]. 經(jīng)濟(jì)地理,1995,15(1):55-60.

    [14] 劉春香,閆國(guó)慶. 我國(guó)農(nóng)業(yè)技術(shù)創(chuàng)新成效研究[J]. 農(nóng)業(yè)經(jīng)濟(jì)問(wèn)題,2012(2):32-37.

    [15] 許慶, 尹榮梁, 章輝. 規(guī)模經(jīng)濟(jì)、規(guī)模報(bào)酬與農(nóng)業(yè)適度規(guī)模經(jīng)營(yíng):基于我國(guó)糧食生產(chǎn)的實(shí)證研究[J]. 經(jīng)濟(jì)研究,2011(3):59-71.

    (編輯陶志寧)

    Analysis of the Partial Least-Square Regression (PLS) on Influencing Factors of Grain Output in Guizhou Province

    LI Fuduo,YANG Xinghong

    (Collegeofmanagement,GuizhouUniversity,Guiyang550025,Guizhou)

    Abstract:Taking Guizhou province as the study area, the influence factors of grain production partial least-square regression model were constructed. It shows that there is a significant linear relationship between independent variable set and dependent variable set. The goodness of fit R2=0.9886, which indicates the regression model with high precision and strong reliability.The research results also show that in the selection of indicators, the sown area to grain, effective irrigated area, total power of agricultural machinery, fertilizer application and disaster area are the key factors affecting grain production in Guizhou.

    Key words:Guizhou province; grain output; influencing factors; least-square regression method

    doi:10.3969/j.issn.1001-8395.2016.01.025

    中圖分類(lèi)號(hào):F321

    文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1001-8395(2016)01-0151-08

    *通信作者簡(jiǎn)介:楊興洪(1971—),女,教授,主要從事農(nóng)業(yè)經(jīng)濟(jì)理論與政策、農(nóng)業(yè)保險(xiǎn)等研究,E-mail:18366622500@163.com

    基金項(xiàng)目:貴州省教育廳高校人文社會(huì)科學(xué)基地研究項(xiàng)目(12JDO22)

    收稿日期:2014-11-11

    猜你喜歡
    糧食生產(chǎn)貴州省影響因素
    貴州省種公牛站
    不動(dòng)產(chǎn)登記地方立法的思考——以貴州省為例
    貴州省黨代會(huì)開(kāi)得最成功
    中原經(jīng)濟(jì)區(qū)農(nóng)業(yè)的可持續(xù)發(fā)展之路
    人民論壇(2016年22期)2016-12-13 11:07:41
    鄭州市糧食生產(chǎn)能力分析
    環(huán)衛(wèi)工人生存狀況的調(diào)查分析
    農(nóng)業(yè)生產(chǎn)性服務(wù)業(yè)需求影響因素分析
    商(2016年27期)2016-10-17 07:09:07
    村級(jí)發(fā)展互助資金組織的運(yùn)行效率研究
    商(2016年27期)2016-10-17 04:40:12
    基于系統(tǒng)論的煤層瓦斯壓力測(cè)定影響因素分析
    科技視界(2016年20期)2016-09-29 13:45:52
    淺論雷州市廣墾東西洋米業(yè)糧食產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展
    亚洲第一av免费看| 一级爰片在线观看| 午夜免费男女啪啪视频观看| 寂寞人妻少妇视频99o| 欧美人与性动交α欧美精品济南到 | 少妇 在线观看| 伊人久久国产一区二区| 成人免费观看视频高清| 免费在线观看黄色视频的| 免费大片18禁| 久热久热在线精品观看| 国产又爽黄色视频| 三上悠亚av全集在线观看| 国产熟女欧美一区二区| 插逼视频在线观看| 天天躁夜夜躁狠狠久久av| 欧美日韩视频精品一区| 黑人高潮一二区| 欧美国产精品一级二级三级| 免费观看无遮挡的男女| 99久国产av精品国产电影| 亚洲综合精品二区| 一本久久精品| 18禁在线无遮挡免费观看视频| 国产精品久久久久久av不卡| 热re99久久国产66热| 大香蕉97超碰在线| 最近手机中文字幕大全| 国内精品宾馆在线| 九色亚洲精品在线播放| 国产无遮挡羞羞视频在线观看| 中文字幕免费在线视频6| 伊人久久国产一区二区| 天天影视国产精品| 少妇被粗大的猛进出69影院 | 一级片免费观看大全| 人体艺术视频欧美日本| 久久国产亚洲av麻豆专区| 国产视频首页在线观看| 亚洲内射少妇av| 亚洲国产精品一区三区| 男男h啪啪无遮挡| 亚洲国产最新在线播放| 国产精品久久久久成人av| 成年人午夜在线观看视频| 少妇 在线观看| 国产精品蜜桃在线观看| 韩国精品一区二区三区 | 2022亚洲国产成人精品| 免费播放大片免费观看视频在线观看| 日韩欧美精品免费久久| 男人添女人高潮全过程视频| 蜜臀久久99精品久久宅男| 中国国产av一级| 在线免费观看不下载黄p国产| 一区二区日韩欧美中文字幕 | 亚洲成人手机| 久久久久久久亚洲中文字幕| 满18在线观看网站| 国产高清国产精品国产三级| 国产成人欧美| 免费av不卡在线播放| 日本wwww免费看| 男女边吃奶边做爰视频| 精品第一国产精品| 免费大片黄手机在线观看| 秋霞伦理黄片| 国产男女超爽视频在线观看| av福利片在线| 久久国产精品大桥未久av| 久久99蜜桃精品久久| 久久这里只有精品19| 亚洲国产欧美在线一区| 日本av手机在线免费观看| 久久久久网色| 97在线人人人人妻| 波野结衣二区三区在线| 99香蕉大伊视频| 国产激情久久老熟女| 好男人视频免费观看在线| 亚洲av综合色区一区| 另类亚洲欧美激情| 九九在线视频观看精品| 国产熟女午夜一区二区三区| 在线观看三级黄色| 亚洲 欧美一区二区三区| 久久久久网色| 美女视频免费永久观看网站| 十八禁高潮呻吟视频| 久久狼人影院| av有码第一页| 天天躁夜夜躁狠狠久久av| 久久久久视频综合| 丰满乱子伦码专区| av电影中文网址| 999精品在线视频| 夜夜爽夜夜爽视频| 午夜福利,免费看| 黄色一级大片看看| 成人亚洲欧美一区二区av| 全区人妻精品视频| 免费高清在线观看视频在线观看| 99香蕉大伊视频| 国产一区亚洲一区在线观看| 视频中文字幕在线观看| 一区二区三区乱码不卡18| 三级国产精品片| 丝瓜视频免费看黄片| 亚洲少妇的诱惑av| 亚洲精华国产精华液的使用体验| 不卡视频在线观看欧美| 国产精品久久久久久久电影| 欧美亚洲 丝袜 人妻 在线| 十分钟在线观看高清视频www| 在线观看美女被高潮喷水网站| 欧美精品亚洲一区二区| 啦啦啦啦在线视频资源| 99久久人妻综合| 精品国产露脸久久av麻豆| 日韩中文字幕视频在线看片| 99热这里只有是精品在线观看| 高清视频免费观看一区二区| 又大又黄又爽视频免费| 午夜福利影视在线免费观看| 久久久精品免费免费高清| av在线观看视频网站免费| 高清不卡的av网站| 欧美97在线视频| 亚洲av国产av综合av卡| 91精品国产国语对白视频| 亚洲伊人色综图| 肉色欧美久久久久久久蜜桃| 综合色丁香网| 久久精品国产a三级三级三级| 国产精品一区二区在线观看99| 日韩中文字幕视频在线看片| 蜜桃国产av成人99| 国产免费福利视频在线观看| av国产精品久久久久影院| 国产欧美日韩一区二区三区在线| 中文字幕制服av| 国产亚洲欧美精品永久| 男女高潮啪啪啪动态图| av卡一久久| 亚洲国产色片| 韩国精品一区二区三区 | 超色免费av| 国产欧美另类精品又又久久亚洲欧美| 91在线精品国自产拍蜜月| 女性被躁到高潮视频| 亚洲图色成人| 亚洲性久久影院| 亚洲精品456在线播放app| 国产色婷婷99| 又大又黄又爽视频免费| 一二三四在线观看免费中文在 | 国产成人精品婷婷| 欧美 日韩 精品 国产| 欧美日韩综合久久久久久| 国产欧美另类精品又又久久亚洲欧美| 国产一级毛片在线| 一级毛片黄色毛片免费观看视频| 亚洲国产精品专区欧美| 色哟哟·www| 免费人妻精品一区二区三区视频| 国产精品人妻久久久影院| 9色porny在线观看| 亚洲精品成人av观看孕妇| 国产免费又黄又爽又色| 亚洲精品aⅴ在线观看| 一本色道久久久久久精品综合| 美女福利国产在线| 黄色 视频免费看| 亚洲三级黄色毛片| 色网站视频免费| 如日韩欧美国产精品一区二区三区| √禁漫天堂资源中文www| 亚洲国产最新在线播放| 国产精品久久久久久精品古装| 国产av一区二区精品久久| 一级毛片黄色毛片免费观看视频| 亚洲精品乱码久久久久久按摩| 中国三级夫妇交换| 男男h啪啪无遮挡| 高清毛片免费看| 精品人妻熟女毛片av久久网站| 国产永久视频网站| 18禁裸乳无遮挡动漫免费视频| 国产成人a∨麻豆精品| 九九爱精品视频在线观看| 国产免费福利视频在线观看| 热re99久久国产66热| 午夜免费观看性视频| 免费在线观看黄色视频的| 国产老妇伦熟女老妇高清| 亚洲,欧美,日韩| 美女国产视频在线观看| 国产精品久久久久久久久免| 亚洲av日韩在线播放| 综合色丁香网| 国产一区有黄有色的免费视频| 蜜臀久久99精品久久宅男| 日日爽夜夜爽网站| 国产成人午夜福利电影在线观看| 一本色道久久久久久精品综合| 天堂俺去俺来也www色官网| 成人亚洲欧美一区二区av| 伦理电影免费视频| 新久久久久国产一级毛片| 亚洲精品乱码久久久久久按摩| 精品少妇内射三级| 国产毛片在线视频| 亚洲综合色网址| 999精品在线视频| 在线看a的网站| 久久精品久久精品一区二区三区| 少妇的逼水好多| 男人操女人黄网站| 亚洲欧美一区二区三区国产| 亚洲欧美精品自产自拍| 亚洲婷婷狠狠爱综合网| 亚洲丝袜综合中文字幕| 18禁裸乳无遮挡动漫免费视频| 超碰97精品在线观看| 香蕉丝袜av| 9色porny在线观看| 巨乳人妻的诱惑在线观看| 一区在线观看完整版| 午夜激情av网站| 91成年电影在线观看| 多毛熟女@视频| 欧美+亚洲+日韩+国产| 成在线人永久免费视频| 中出人妻视频一区二区| 黄色 视频免费看| 久久久久久久精品吃奶| 国产精品电影一区二区三区 | 狠狠狠狠99中文字幕| 久久人妻av系列| 身体一侧抽搐| 亚洲性夜色夜夜综合| 国产精品永久免费网站| 日韩成人在线观看一区二区三区| 一级作爱视频免费观看| 中文亚洲av片在线观看爽 | 国产乱人伦免费视频| 国产精品.久久久| 伦理电影免费视频| 精品一品国产午夜福利视频| 中文字幕色久视频| 日本五十路高清| 久久九九热精品免费| 无限看片的www在线观看| 美女高潮喷水抽搐中文字幕| 成人18禁在线播放| 精品一区二区三卡| 首页视频小说图片口味搜索| 老熟女久久久| 国产精品美女特级片免费视频播放器 | 免费观看a级毛片全部| 老司机在亚洲福利影院| 性色av乱码一区二区三区2| 午夜免费鲁丝| 亚洲视频免费观看视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产亚洲一区二区精品| 99精品欧美一区二区三区四区| 国产精品国产av在线观看| 12—13女人毛片做爰片一| 他把我摸到了高潮在线观看| 高清毛片免费观看视频网站 | 极品教师在线免费播放| 操美女的视频在线观看| 咕卡用的链子| 国产一区有黄有色的免费视频| 亚洲九九香蕉| 国产成人欧美在线观看 | videosex国产| 国产单亲对白刺激| 欧美日韩av久久| 精品国产美女av久久久久小说| 麻豆av在线久日| 久久影院123| 不卡一级毛片| 一个人免费在线观看的高清视频| 亚洲第一青青草原| 欧美日本中文国产一区发布| 亚洲色图 男人天堂 中文字幕| 国产成人av激情在线播放| 18禁观看日本| 国产精品av久久久久免费| 久久天堂一区二区三区四区| 欧美老熟妇乱子伦牲交| 大香蕉久久成人网| 久久热在线av| 中文字幕制服av| 中文字幕人妻丝袜制服| 国内久久婷婷六月综合欲色啪| 男人舔女人的私密视频| 中文字幕人妻丝袜制服| 多毛熟女@视频| 国产99白浆流出| 精品久久久久久久毛片微露脸| 欧美黄色淫秽网站| 久久久久国产精品人妻aⅴ院 | 嫁个100分男人电影在线观看| 欧美最黄视频在线播放免费 | 999久久久国产精品视频| 亚洲国产精品一区二区三区在线| 超碰成人久久| 性色av乱码一区二区三区2| 亚洲av片天天在线观看| 亚洲人成77777在线视频| 精品国内亚洲2022精品成人 | 亚洲avbb在线观看| 国产又色又爽无遮挡免费看| 欧美一级毛片孕妇| 亚洲欧美色中文字幕在线| 性少妇av在线| 国产野战对白在线观看| 免费不卡黄色视频| 精品亚洲成a人片在线观看| 在线观看日韩欧美| 久久婷婷成人综合色麻豆| 黄片大片在线免费观看| 久久香蕉精品热| 中亚洲国语对白在线视频| 亚洲一区中文字幕在线| av网站免费在线观看视频| 老司机深夜福利视频在线观看| 午夜福利免费观看在线| 午夜日韩欧美国产| 露出奶头的视频| 啦啦啦视频在线资源免费观看| 十八禁网站免费在线| 久久午夜综合久久蜜桃| 宅男免费午夜| 亚洲在线自拍视频| 天堂动漫精品| 啦啦啦在线免费观看视频4| 黄片小视频在线播放| 亚洲精品国产色婷婷电影| 午夜福利乱码中文字幕| 高清在线国产一区| 国产91精品成人一区二区三区| 亚洲av欧美aⅴ国产| 王馨瑶露胸无遮挡在线观看| 久久中文字幕人妻熟女| av电影中文网址| videosex国产| 精品卡一卡二卡四卡免费| 国产精品九九99| 人人妻人人添人人爽欧美一区卜| 中文字幕人妻丝袜一区二区| 高潮久久久久久久久久久不卡| 午夜福利,免费看| 最近最新中文字幕大全电影3 | 国产极品粉嫩免费观看在线| 日韩欧美国产一区二区入口| 亚洲成a人片在线一区二区| 精品人妻1区二区| 老司机午夜十八禁免费视频| 色在线成人网| 欧美成人免费av一区二区三区 | 18禁黄网站禁片午夜丰满| 国产欧美日韩一区二区三区在线| 久久久国产成人免费| 国产无遮挡羞羞视频在线观看| 欧美乱妇无乱码| 999久久久精品免费观看国产| 久久九九热精品免费| 久久久久久免费高清国产稀缺| 两性夫妻黄色片| 黄网站色视频无遮挡免费观看| 黄片大片在线免费观看| 久久婷婷成人综合色麻豆| 免费观看a级毛片全部| 每晚都被弄得嗷嗷叫到高潮| 女警被强在线播放| 日本五十路高清| e午夜精品久久久久久久| 涩涩av久久男人的天堂| 久久精品国产99精品国产亚洲性色 | 黄色a级毛片大全视频| 日韩成人在线观看一区二区三区| 免费在线观看亚洲国产| 搡老岳熟女国产| 狂野欧美激情性xxxx| 久久国产精品男人的天堂亚洲| 中出人妻视频一区二区| 国产蜜桃级精品一区二区三区 | 亚洲中文av在线| 亚洲成人手机| 男女之事视频高清在线观看| 国产精品99久久99久久久不卡| 1024视频免费在线观看| 国产精品久久久久久精品古装| 别揉我奶头~嗯~啊~动态视频| 悠悠久久av| 午夜精品久久久久久毛片777| 免费在线观看影片大全网站| 精品久久久久久电影网| 国产xxxxx性猛交| www.精华液| 黄色女人牲交| 黄色a级毛片大全视频| 啦啦啦 在线观看视频| 中文字幕av电影在线播放| 国产亚洲欧美在线一区二区| 欧美日韩国产mv在线观看视频| 欧美日韩乱码在线| 国产精品一区二区精品视频观看| 久久精品亚洲av国产电影网| 大型黄色视频在线免费观看| 国产区一区二久久| 久久精品aⅴ一区二区三区四区| 正在播放国产对白刺激| 欧美黑人欧美精品刺激| 国产成人精品久久二区二区免费| 亚洲av第一区精品v没综合| 成年版毛片免费区| 成人手机av| 成人av一区二区三区在线看| 亚洲欧美激情综合另类| 久久狼人影院| 99热网站在线观看| 一级黄色大片毛片| 无遮挡黄片免费观看| 村上凉子中文字幕在线| 亚洲成人国产一区在线观看| 熟女少妇亚洲综合色aaa.| 男女午夜视频在线观看| 如日韩欧美国产精品一区二区三区| 女性被躁到高潮视频| 久久久久久久午夜电影 | 欧美黄色淫秽网站| 亚洲精品中文字幕在线视频| 免费看十八禁软件| 色婷婷久久久亚洲欧美| 天天添夜夜摸| 中文字幕av电影在线播放| 欧美黑人欧美精品刺激| 欧美国产精品一级二级三级| 啦啦啦免费观看视频1| aaaaa片日本免费| 久久婷婷成人综合色麻豆| 夜夜爽天天搞| 免费在线观看视频国产中文字幕亚洲| 国产精品一区二区在线观看99| av国产精品久久久久影院| 99精国产麻豆久久婷婷| 精品福利观看| 欧美色视频一区免费| 精品人妻1区二区| 国产亚洲欧美精品永久| 高清视频免费观看一区二区| 9191精品国产免费久久| 亚洲精品中文字幕一二三四区| 侵犯人妻中文字幕一二三四区| 国产精品成人在线| 婷婷精品国产亚洲av在线 | 国产精品二区激情视频| 午夜日韩欧美国产| 夜夜爽天天搞| 亚洲精品中文字幕在线视频| 午夜福利一区二区在线看| av天堂久久9| 国产三级黄色录像| 女性被躁到高潮视频| 桃红色精品国产亚洲av| 人人妻,人人澡人人爽秒播| 亚洲综合色网址| 精品乱码久久久久久99久播| 王馨瑶露胸无遮挡在线观看| 啦啦啦 在线观看视频| 成人永久免费在线观看视频| 国产黄色免费在线视频| 久久人人97超碰香蕉20202| 亚洲精品乱久久久久久| 亚洲第一青青草原| 色播在线永久视频| 色精品久久人妻99蜜桃| 国产成人欧美| 大码成人一级视频| 男女床上黄色一级片免费看| 国产在视频线精品| 最新的欧美精品一区二区| 久久久久久久国产电影| 欧美日韩亚洲国产一区二区在线观看 | 男女免费视频国产| 精品国产超薄肉色丝袜足j| 日本撒尿小便嘘嘘汇集6| 日韩欧美一区视频在线观看| 国产99久久九九免费精品| 狠狠婷婷综合久久久久久88av| 国产精品免费视频内射| 成人黄色视频免费在线看| 国产深夜福利视频在线观看| 亚洲第一欧美日韩一区二区三区| 80岁老熟妇乱子伦牲交| 免费在线观看亚洲国产| 亚洲全国av大片| 男女之事视频高清在线观看| 香蕉丝袜av| 国产片内射在线| 日本五十路高清| 最新的欧美精品一区二区| 精品第一国产精品| 国产精品久久久av美女十八| 日本一区二区免费在线视频| 老司机亚洲免费影院| 黄色视频不卡| 欧美另类亚洲清纯唯美| 搡老岳熟女国产| 中文字幕人妻熟女乱码| 妹子高潮喷水视频| 黄色 视频免费看| 在线观看午夜福利视频| 久久国产精品影院| 一级,二级,三级黄色视频| 亚洲自偷自拍图片 自拍| 19禁男女啪啪无遮挡网站| 美女福利国产在线| 亚洲av成人av| 夫妻午夜视频| 日韩 欧美 亚洲 中文字幕| 久久久久视频综合| 日韩欧美国产一区二区入口| 日本一区二区免费在线视频| 久久久久久人人人人人| 成人手机av| 中文亚洲av片在线观看爽 | 婷婷丁香在线五月| 亚洲国产精品sss在线观看 | 村上凉子中文字幕在线| 在线av久久热| 天天躁狠狠躁夜夜躁狠狠躁| 桃红色精品国产亚洲av| 老熟妇乱子伦视频在线观看| 久久性视频一级片| 精品人妻1区二区| 久久久久国产一级毛片高清牌| 国产一区在线观看成人免费| 夫妻午夜视频| 99re6热这里在线精品视频| 免费在线观看日本一区| 人人妻人人澡人人爽人人夜夜| 欧美精品啪啪一区二区三区| 高清毛片免费观看视频网站 | 欧美性长视频在线观看| 久久精品成人免费网站| 精品国产亚洲在线| 久久性视频一级片| 久久婷婷成人综合色麻豆| 免费在线观看黄色视频的| 日本五十路高清| 欧洲精品卡2卡3卡4卡5卡区| 国产欧美日韩精品亚洲av| 美国免费a级毛片| cao死你这个sao货| 亚洲欧美精品综合一区二区三区| 天天操日日干夜夜撸| 国产亚洲av高清不卡| 亚洲国产欧美一区二区综合| 亚洲精品成人av观看孕妇| 两性夫妻黄色片| 天天躁夜夜躁狠狠躁躁| 女人精品久久久久毛片| 欧美精品高潮呻吟av久久| 亚洲第一欧美日韩一区二区三区| 精品视频人人做人人爽| 国产片内射在线| 亚洲成av片中文字幕在线观看| 久久 成人 亚洲| 日韩欧美免费精品| 水蜜桃什么品种好| 国产精品偷伦视频观看了| 12—13女人毛片做爰片一| aaaaa片日本免费| 精品一区二区三卡| 国产精品一区二区在线不卡| 久久久水蜜桃国产精品网| 欧美日韩国产mv在线观看视频| 国产成人精品久久二区二区91| 国产亚洲精品久久久久久毛片 | 老熟妇仑乱视频hdxx| av片东京热男人的天堂| 在线观看免费视频网站a站| 成人18禁高潮啪啪吃奶动态图| 日本黄色日本黄色录像| 欧美乱码精品一区二区三区| 亚洲一码二码三码区别大吗| 天堂俺去俺来也www色官网| 美女 人体艺术 gogo| 一a级毛片在线观看| 精品少妇久久久久久888优播| 免费高清在线观看日韩| 一级毛片高清免费大全| 极品教师在线免费播放| 老司机在亚洲福利影院| 中国美女看黄片| 曰老女人黄片| 99riav亚洲国产免费| 黄色 视频免费看| 在线av久久热| 人人妻人人爽人人添夜夜欢视频| 午夜福利影视在线免费观看| 大香蕉久久成人网| 黄色 视频免费看| 国产精品久久久久成人av| 国产三级黄色录像| 热99国产精品久久久久久7| 日韩大码丰满熟妇| 成人国产一区最新在线观看| 大片电影免费在线观看免费|