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

        ?

        基于SPSS分析的土地生態(tài)安全評價

        2022-07-18 09:06:08劉泳昕
        農(nóng)業(yè)與技術(shù) 2022年13期
        關(guān)鍵詞:生態(tài)評價研究

        劉泳昕

        (長安大學(xué)土地工程學(xué)院,陜西 西安 710061)

        近年來,我國在經(jīng)濟(jì)社會各領(lǐng)域均成績斐然,新時期下的更多發(fā)展成就也不斷為民眾享用,隨著我國民眾的物質(zhì)生活水平日益改善,更多人的消費要求也逐漸由對傳統(tǒng)物質(zhì)商品和新文化產(chǎn)品的需要轉(zhuǎn)變?yōu)閷ΜF(xiàn)代生態(tài)商品的需要。自黨的十七大報告中明確了構(gòu)建生態(tài)文明的總體目標(biāo)之后,黨的十八大也在報告中指出要將社會主義生態(tài)文明建設(shè)擺在突出地位,黨的十九大報告進(jìn)一步明確提出,“我們所要建立的社會主義現(xiàn)代化是人與自然和諧共存的社會主義現(xiàn)代化”?;谶@種背景,對我國生態(tài)安全現(xiàn)狀進(jìn)行深入研究顯得尤為重要。

        “絕對數(shù)量大,相對數(shù)量小”已經(jīng)成為并將長期成為我國土地資源狀況的基本國情。改革開放以來,人口的持續(xù)增長和經(jīng)濟(jì)社會的飛速發(fā)展在對我國各行各業(yè)不斷增長的土地需求中體現(xiàn)出來,土地資源數(shù)量和質(zhì)量方面的供給短缺限制了我國經(jīng)濟(jì)社會高質(zhì)量發(fā)展和可持續(xù)發(fā)展的步伐。目前,由于人類活動對土地資源的不合理開發(fā),土地生態(tài)問題日益突出,水土流失、土地生態(tài)惡化、土壤污染等影響著人類生產(chǎn)生活的各個方面。

        為了能夠有序有效地管理這些問題,實現(xiàn)可持續(xù)發(fā)展,了解土地生態(tài)系統(tǒng)的基本狀況是十分必要的,因此開展土地生態(tài)安全評價已經(jīng)變得非常重要。

        1 國內(nèi)外研究現(xiàn)狀

        1.1 國內(nèi)研究現(xiàn)狀

        本研究以“土地生態(tài)安全評價”為主題詞,通過文獻(xiàn)管理軟件Noteexprss從CNKI中國知網(wǎng)中檢索篩選出與土地生態(tài)安全評價相關(guān)文獻(xiàn)共436篇。從論文作者、研究區(qū)域、研究方法、發(fā)表年限4個方面對樣本數(shù)據(jù)進(jìn)行統(tǒng)計分析,并對當(dāng)前土地生態(tài)安全評價的研究現(xiàn)狀做出總結(jié)性論述。

        通過Excel軟件對不同時期與土地生態(tài)安全評價相關(guān)文獻(xiàn)發(fā)表數(shù)量進(jìn)行統(tǒng)計,分析隨年代變化的文獻(xiàn)發(fā)表數(shù)量變化規(guī)律,總結(jié)土地生態(tài)系統(tǒng)研究的發(fā)展趨勢,如圖1所示。由圖1可知,土地生態(tài)系統(tǒng)研究的文獻(xiàn)數(shù)量總體呈增長態(tài)勢,其中,2000—2005年文獻(xiàn)數(shù)量<10篇,反映土地生態(tài)系統(tǒng)研究處于起步階段;2005—2015年文獻(xiàn)數(shù)量位于10~50篇,土地生態(tài)系統(tǒng)研究處于發(fā)展階段;2010—2020年文獻(xiàn)數(shù)量>50篇,土地生態(tài)系統(tǒng)研究處于高速發(fā)展階段,根據(jù)線性趨勢預(yù)測,未來這一研究領(lǐng)域的文獻(xiàn)數(shù)量將會持續(xù)增長。

        通過Noteexpress統(tǒng)計,目前國內(nèi)對土地生態(tài)系統(tǒng)評價的主要研究區(qū)集中在湖南省、安徽省、河南省、黑龍江省以及珠三角等東部、南部地區(qū),對西北地區(qū)的研究較少。發(fā)文量超過5篇的作者共6人,其中何如海、王鵬、雷國平、張小虎等被認(rèn)為是該領(lǐng)域的核心研究作者,其余作者均未對土地生態(tài)系統(tǒng)評價進(jìn)行深入的研究。研究方法以PSR模型、主成分分析法、熵權(quán)法為主,評價模型及方法上實現(xiàn)了多元化發(fā)展,如表1所示。

        圖1 2000—2021年土地生態(tài)安全評價相關(guān)文獻(xiàn)發(fā)表數(shù)量年度變化

        表1 土地生態(tài)安全評價相關(guān)文獻(xiàn)研究基本情況

        綜上所述,國內(nèi)相較國外對土地生態(tài)安全評價有關(guān)的研究起步晚。21世紀(jì)初,由于我國生態(tài)環(huán)境問題日益突出,土地生態(tài)安全水平對地區(qū)乃至國家發(fā)展和安全的影響才逐漸受到學(xué)者的關(guān)注,一直到21世紀(jì)10年代才逐漸成為國內(nèi)學(xué)者研究的熱點,研究步入高速發(fā)展階段。

        1.2 國外研究現(xiàn)狀

        1941年Aldo Leopold在土地功能狀況的評價中提出“土地健康”的概念,將國外對土地生態(tài)安全方面的研究最早追溯期限推至20世紀(jì)40年代[6]。土地生態(tài)安全評價指標(biāo)體系最具代表性的成果是1996年由世界銀行、聯(lián)合國開發(fā)署、聯(lián)合國糧農(nóng)組織、聯(lián)合國環(huán)境署共同提出的“壓力-狀態(tài)-響應(yīng)”(PSR)框架模型,該模型主要目的是利用土地質(zhì)量指標(biāo)檢測人類活動對土地資源環(huán)境的壓力和這些壓力對土地質(zhì)量狀態(tài)的影響及社會對這些壓力的響應(yīng),以此來反映土地質(zhì)量的好壞。在“壓力-狀態(tài)-響應(yīng)”(PSR)模型和“驅(qū)動力-狀態(tài)-響應(yīng)”(DSR)模型的基礎(chǔ)上,歐洲環(huán)境署提出了“驅(qū)動力-壓力-狀態(tài)-影響-響應(yīng)”(DPSIR)框架模型。Paracchini ML在2011年基于“經(jīng)濟(jì)-環(huán)境-社會”(EES)模型構(gòu)建評價指標(biāo)體系[5]。PSR、DSR、DPSIR、EES等諸多模型的廣泛應(yīng)用為國內(nèi)外土地生態(tài)安全評價提供了可以借鑒和參考的模型框架和指標(biāo)因子。

        2 土地生態(tài)安全評價指標(biāo)體系

        2.1 構(gòu)建指標(biāo)體系

        根據(jù)選取方法的合理性、評價指標(biāo)可比性、樣本數(shù)據(jù)易得性等原則,以生態(tài)文明建設(shè)和可持續(xù)發(fā)展的視角,對國內(nèi)外研究成果進(jìn)行綜合參考,采用德爾菲法篩選評價因子,在剔除無效和相關(guān)性較弱的因子后,構(gòu)建本文所采用的評價指標(biāo)體系,如表2所示。

        表2 土地生態(tài)安全評價指標(biāo)體系

        2.2 數(shù)據(jù)來源

        本文從《中國統(tǒng)計年鑒—2021》中搜集數(shù)據(jù),研究對象是全國通過行政區(qū)劃劃分的土地狀況,采用因子-聚類分析方法對數(shù)據(jù)進(jìn)行挖掘和運算。

        2.3 數(shù)據(jù)標(biāo)準(zhǔn)化處理

        本文利用Z-score標(biāo)準(zhǔn)化對指標(biāo)數(shù)據(jù)進(jìn)行分析處理,旨在減少誤差,先獲得標(biāo)準(zhǔn)化矩陣再進(jìn)行綜合對比分析。

        3 因子分析

        3.1 適用性檢驗

        本文對數(shù)據(jù)的適用性運用KMO和巴特利特的方法檢驗標(biāo)準(zhǔn)化后,結(jié)果如表3所示。KMO(Kaiser-Meyer-Olkin)檢驗結(jié)果為0.544>0.5;巴特利特球形度檢驗(Bartlett)的P值0.000<0.05,證明標(biāo)準(zhǔn)化后的指標(biāo)數(shù)據(jù)符合因子分析條件。

        表3 KMO和巴特利特檢驗

        3.2 公因子方差

        公因子方差表示原始變量被公因子解釋的比例和表達(dá)的程度,是成分矩陣中每一行中每個載荷值的平方和。公因子方差表中的“提取”值與變量指標(biāo)可以被公因子表達(dá)的程度呈正相關(guān),一般當(dāng)提取值>0.5時說明可以被表達(dá),當(dāng)提取值>0.7時可以說明變量能被公因子表達(dá)得很合理。依據(jù)表4可以看出,本研究對象的14個指標(biāo)的“提取”值最小為0.765,14個指標(biāo)的“提取”值都大于0.7。因此得出,選取的14個指標(biāo)都可以很好地被公因子表達(dá)。

        表4 公因子方差

        3.3 提取因子

        主成分的累計方差貢獻(xiàn)率要求在85%~95%的范圍內(nèi),如表5所示,主累計方差貢獻(xiàn)率86.084%,特征根值>1,可見提取的5個因子作為主成分是可行的。

        表5 總方差解釋

        3.4 因子命名

        為了讓主成分因子定義更直觀更具體,對主成分因子定義進(jìn)行說明,由于第1因子在人均地區(qū)生產(chǎn)總值(X4)、城市化水平(X3)系數(shù)較大,故命名為“經(jīng)濟(jì)響應(yīng)”因子;第2因子中系數(shù)較大的是總?cè)丝跀?shù)(X1)和生活垃圾清運量(X10),因而命名為“公共壓力”因子;第3因子在水土流失治理面積(X13)、一般工業(yè)固體廢物綜合利用量(X11)系數(shù)較大,故命名為“公共響應(yīng)”因子;第4因子在失業(yè)率(X2)、第二產(chǎn)業(yè)比重(X8)系數(shù)較大,從經(jīng)濟(jì)層面考慮命名為“經(jīng)濟(jì)狀態(tài)”;森林覆蓋率(X6)在第5因子中系數(shù)較大,故命名為“環(huán)境影響”因子。

        3.5 建立因子得分模型

        基于如表6所示的成分得分系數(shù),得出“經(jīng)濟(jì)響應(yīng)”、“公共壓力”、“公共響應(yīng)”、“經(jīng)濟(jì)狀態(tài)”、“環(huán)境影響”5大因子的得分模型。

        F1=-0.054X1+0.222X2+0.298X3+0.283X4-0.074X5+0.001X6-0.276X7+0.020X8+0.180X9+0.047X10-0.054X11+0.045X12-0.027X13-0.066X14

        F2=0.337X1-0.223X2+0.005X3+0.013X4+0.023X5+0.047X6+0.024X7-0.039X8+0.019X9+0.306X10+0.208X11+0.236X12-0.001X13+0.139X14

        F3=-0.035X1+0.404X2+0.152X3+0.025X4-0.264X5+0.138X6+0.048X7-0.048X8+0.010X9-0.044X10+0.077X11-0.020X12+0.335X13+0.268X14

        F4=-0.110X1+0.293X2-0.082X3+0.019X4+0.162X5-0.060X6-0.335X7+0.606X8-0.319X9-0.087X10-0.065X11+0.097X12-0.060X13-0.145X14

        F5=0.082X1+0.262X2+0.019X3+0.049X4+0.337X5+0.604X6+0.095X7-0.011X8-0.058X9+0.125X10-0.278X11-0.048X12+0.103X13-0.020X14

        通過因子分析法,追求損失最少的最優(yōu)解,壓縮多個指標(biāo)為5個綜合指標(biāo);用各個因子占全部因子的權(quán)重進(jìn)行賦權(quán)整合,算出土地生態(tài)安全綜合得分F;將F值按照從大到小順序進(jìn)行排序,得出各省F綜合排名,如表7所示。

        F=0.3438F1+0.2817F2+0.1641F3+0.1171F4+0.0933F5

        表6 成分得分系數(shù)矩陣

        表7 土地生態(tài)安全的綜合排名

        4 聚類分析

        4.1 K-均值聚類分析

        本文使用SPSS 25.0進(jìn)行K-均值聚類(K-Means clustering algorithm),將31個省市(自治區(qū))的土地生態(tài)安全水平分為低風(fēng)險區(qū)(一類)、較低風(fēng)險區(qū)(二類)、中風(fēng)險區(qū)(三類)、高風(fēng)險區(qū)(四類)共4類,結(jié)果如表8所示。

        4.2 系統(tǒng)聚類分析

        本文利用SPSS 25.0進(jìn)行系統(tǒng)聚類,將31個省市(自治區(qū))的土地生態(tài)安全水平分為低風(fēng)險區(qū)(一類)、較低風(fēng)險區(qū)(二類)、中風(fēng)險區(qū)(三類)、高風(fēng)險區(qū)(四類)共4類,結(jié)果如表9所示,譜系圖如圖2所示。

        表8 K-均值聚類分析結(jié)果

        表9 系統(tǒng)聚類分析結(jié)果

        圖2 基于系統(tǒng)聚類分析的土地生態(tài)安全譜系圖

        4.3 2種聚類分析結(jié)果比較

        4.3.1 關(guān)于聚類個數(shù)

        盡管K-均值聚類和系統(tǒng)聚類選擇聚類個數(shù)方面都需要人為進(jìn)行選擇,但是系統(tǒng)聚類優(yōu)于K-均值聚類的部分在于其可以根據(jù)生成的譜系圖自發(fā)選擇出分類效果相對較好的類別。

        4.3.2 關(guān)于聚類效果

        遼寧、上海、江蘇、浙江、福建、山東、廣東在2種方法中都自成一類,可見對于差異比較大的個案,2種方法都能區(qū)分出來,因而聚類效果基本一致。但是通過譜系圖,系統(tǒng)聚類能看出個案之間關(guān)系的遠(yuǎn)近。

        4.3.3 關(guān)于聚類結(jié)果

        2種方法的最終聚類結(jié)果基本相同,但也存在差異。黑龍江、廣西、甘肅在K-均值聚類中屬于三類,但在系統(tǒng)聚類中屬于四類。安徽、貴州在K-均值聚類中屬于二類,但在系統(tǒng)聚類中屬于三類。

        4.3.4 關(guān)于聚類方法不足

        系統(tǒng)聚類和K均值聚類分析都存在各自的局限性和不足之處。在進(jìn)行系統(tǒng)聚類分析時,每個個案一經(jīng)歸為某一類,之后該個案就不允許再被改變類別,但是個案最初的分類可能不是個案歸類的最優(yōu)選擇。K均值聚類的分類數(shù)需要提前確定固定的值,而且分類數(shù)確定后不能更改。因此,K均值聚類的分類數(shù)的確定需要根據(jù)行業(yè)特點、對研究對象的預(yù)先了解、對研究事物的判斷來確定預(yù)期的分類數(shù)。因此,在之后的實際研究中,為了更合理地確定分類數(shù),可以先對部分代表性樣本進(jìn)行系統(tǒng)聚類分析,從而輔助決定應(yīng)該在K均值聚類分析中分為幾類,使分類結(jié)果更加符合實際、更加合理。

        5 結(jié)論與建議

        5.1 結(jié)論

        運用因子分析法將選取的影響土地生態(tài)安全程度的14個評價指標(biāo)綜合劃分成經(jīng)濟(jì)響應(yīng)、公共壓力、公共響應(yīng)、經(jīng)濟(jì)狀態(tài)、環(huán)境影響共5個主要因子指標(biāo)。第一類因子主要由“人均地區(qū)生產(chǎn)總值”、“城市化水平”構(gòu)成,命名為“經(jīng)濟(jì)響應(yīng)”因子;第二類因子主要包括“生活垃圾清運量”、“總?cè)丝跀?shù)”,命名為“公共壓力”因子;第三類因子主要包括“一般工業(yè)固體廢物綜合利用量”、“水土流失治理面積”,命名為“公共響應(yīng)”因子;第四類因子主要包括“第二產(chǎn)業(yè)比重”、“失業(yè)率”,從經(jīng)濟(jì)層面考慮命名為“經(jīng)濟(jì)狀態(tài)”因子;第五類因子主要包含“森林覆蓋率”,命名為“環(huán)境影響”因子。

        運用SPSS 25.0軟件進(jìn)行運算得出所有省市(自治區(qū))的綜合得分,并對得分按照由高到低的原則進(jìn)行排序,得分與排名成正相關(guān)關(guān)系,分值越高排名越靠前,排名越前的地區(qū)土地生態(tài)安全越好,反之亦然。

        K-均值聚類和系統(tǒng)聚類2種方法分析得出的分類結(jié)果大致吻合。現(xiàn)將31個省市(自治區(qū))依據(jù)F值的高低劃分為高風(fēng)險區(qū)、中風(fēng)險區(qū)、較低風(fēng)險區(qū)、低風(fēng)險區(qū)4類,2種分類結(jié)果不一致的省市(自治區(qū))視為處于搖擺狀態(tài)。土地生態(tài)安全狀態(tài)處于低風(fēng)險的省市,包括廣東、山東、遼寧、上海、浙江、江蘇、福建;北京、天津、河北、山西、內(nèi)蒙古、江西、河南、湖北、重慶、四川、陜西、云南的土地生態(tài)安全水平基本上達(dá)到一般安全狀態(tài),處于較低風(fēng)險范圍;其他還未達(dá)到安全狀態(tài)的省市(自治區(qū)),需要積極推進(jìn)土地生態(tài)文明建設(shè),轉(zhuǎn)變發(fā)展思路和方式,在規(guī)劃和實施中更加注重土地生態(tài)環(huán)境的改善。安徽、貴州屬于在較低風(fēng)險和中風(fēng)險之間搖擺的省份,黑龍江、廣西、甘肅屬于在中風(fēng)險和高風(fēng)險之間搖擺的省份(自治區(qū)),需要維護(hù)好和管理好土地生態(tài)系統(tǒng),因地制宜地推進(jìn)城鄉(xiāng)綠化覆蓋和植被恢復(fù)工程,不斷治理水土流失和土地退化,不斷完善污水、環(huán)衛(wèi)等基礎(chǔ)設(shè)施建設(shè),加強自身水土保持能力與防治環(huán)境污染并重,優(yōu)化土地生態(tài),總體上提高本地區(qū)土地生態(tài)安全水平,使土地生態(tài)安全水平向更安全的方向發(fā)展、搖擺到更高標(biāo)準(zhǔn)的區(qū)間內(nèi)。

        5.2 建議

        5.2.1 將習(xí)近平生態(tài)文明思想融入到生態(tài)文化建設(shè)中

        需要通過顯性手段和隱性手段雙管齊下,不斷培養(yǎng)公眾對社會主義生態(tài)文明思想的心理認(rèn)同,不斷提升生態(tài)文明建設(shè)的實踐性,使公眾將心理認(rèn)同感逐步轉(zhuǎn)化成為保護(hù)生態(tài)環(huán)境身體力行的內(nèi)生動力。堅持以人為本,舉辦世界環(huán)境日、世界地球日、生物多樣性國際日、地球一小時等相關(guān)活動,鼓勵創(chuàng)作和傳播弘揚正能量的自然、災(zāi)難等類型的文學(xué)影視作品,不斷通過啟發(fā)自覺、培養(yǎng)自覺、督促自覺,讓生態(tài)文化在全社會扎根。

        5.2.2 科學(xué)編制保護(hù)規(guī)劃,大力治理生態(tài)環(huán)境

        貫徹“五大發(fā)展理念”,堅持戰(zhàn)略引領(lǐng)和剛性管控并重,在規(guī)劃層面積極推進(jìn)多規(guī)合一、協(xié)調(diào)發(fā)展的實現(xiàn),切實解決現(xiàn)有各類規(guī)劃存在沖突和缺乏體系的問題,健全有效協(xié)調(diào)機制,并結(jié)合各地區(qū)特點和發(fā)展需求,突破行政和地域限制,形成多方合力,推動保護(hù)和治理措施得到層層落實。必須要維護(hù)和管理好土地生態(tài)系統(tǒng)的質(zhì)量水平和數(shù)量水平,保證其質(zhì)量不降低、數(shù)量不減少,對生態(tài)污染地區(qū)要進(jìn)行修復(fù)整治,堅持以預(yù)防為主、防消結(jié)合的方針,做好土地生態(tài)良好地區(qū)、土地生態(tài)安全地區(qū)的長期保護(hù)工作,總體上提高土地生態(tài)安全的水平?!笆濉逼陂g實施的封山育林、水土流失綜合治理、退化林修復(fù)、水源涵養(yǎng)區(qū)建設(shè)等工程使我國國土整治修復(fù)工作取得一定進(jìn)展。在“十四五”期間需要對“十三五”期間取得的成果進(jìn)行鞏固,以保質(zhì)增量為原則,又快又穩(wěn)地推進(jìn)生態(tài)環(huán)境治理,積極推進(jìn)雨污水管網(wǎng)覆蓋率提高、廢氣廢物廢水回收再利用和凈化后排放以及垃圾無害化處理等類型項目,不斷提高土地生態(tài)安全綜合水平。

        5.2.3 建立健全土地生態(tài)安全價值實現(xiàn)機制,實現(xiàn)經(jīng)濟(jì)、社會和生態(tài)效益的統(tǒng)一

        以政府為主導(dǎo),按照行政區(qū)劃分級進(jìn)行生態(tài)產(chǎn)品價值評價,解碼各地區(qū)生態(tài)本底,為保護(hù)規(guī)劃編制、利用方式選取提供資源依據(jù)。充分借鑒和參考國內(nèi)外成功經(jīng)驗,結(jié)合本國實際情況,通過生態(tài)銀行、碳達(dá)峰、碳中和等方式來拓寬土地生態(tài)安全維護(hù)資金來源渠道,調(diào)動土地生態(tài)安全供應(yīng)方保護(hù)與提高土地生態(tài)安全水平的積極性和主動性,提升土地生態(tài)安全需求方(受益方)的保護(hù)意識與維權(quán)意識,實現(xiàn)生態(tài)效益、經(jīng)濟(jì)效益以及社會效益的有機統(tǒng)一,真正徹底地實現(xiàn)生態(tài)產(chǎn)品價值。

        猜你喜歡
        生態(tài)評價研究
        FMS與YBT相關(guān)性的實證研究
        SBR改性瀝青的穩(wěn)定性評價
        石油瀝青(2021年4期)2021-10-14 08:50:44
        “生態(tài)養(yǎng)生”娛晚年
        遼代千人邑研究述論
        住進(jìn)呆萌生態(tài)房
        生態(tài)之旅
        視錯覺在平面設(shè)計中的應(yīng)用與研究
        科技傳播(2019年22期)2020-01-14 03:06:54
        EMA伺服控制系統(tǒng)研究
        基于Moodle的學(xué)習(xí)評價
        生態(tài)
        精品亚洲国产成人av色哟哟| 囯产精品无码一区二区三区AV| 男子把美女裙子脱了摸她内裤| 中文字幕一区二区三区日日骚 | 青青草国产成人99久久| 免费大学生国产在线观看p| 中文字幕亚洲一区视频| 亚洲av一二三区成人影片| 大肉大捧一进一出好爽视频mba| 国产精品成人午夜久久| 亚洲国产精品成人一区| 青青草免费手机视频在线观看| 女人张开腿让男桶喷水高潮| 亚洲国产精品自拍一区| 日韩熟女一区二区三区 | 一区二区三区精彩视频在线观看| 就爱射视频在线视频在线| 亚洲av无码成人网站在线观看| 国内精品久久久久久久久久影院| 日本一区二区视频免费观看| 亚洲国产综合在线亚洲区亚洲av| 成人国产精品一区二区网站公司| 日韩无码无播放器视频| 亚洲人av毛片一区二区| 不卡一区二区视频日本| 麻豆成人精品国产免费| 少妇无码av无码去区钱| 一区二区三区免费自拍偷拍视频| 少妇高潮av久久久久久| 区久久aaa片69亚洲| 久久HEZYO色综合| 丝袜美腿视频一区二区 | 国产精品国产三级国产AvkTV| 亚洲av午夜福利精品一区不卡| 亚洲av综合av国产av中文| 国产熟妇高潮呻吟喷水| 久久国产成人免费网站| 国产一区二区不卡av| 一本色道久久88加勒比—综合| 国产乱人伦偷精品视频| 久久久99久久久国产自输拍|