——以河北省巨鹿縣為例"/>

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

        ?

        基于遙感數(shù)據(jù)的縣域土地利用模擬及經(jīng)濟(jì)效益評(píng)價(jià)
        ——以河北省巨鹿縣為例

        2021-04-22 03:22:40劉澤琴米夢(mèng)凡楊秀峰趙子輝解春苓
        衛(wèi)星應(yīng)用 2021年3期

        文 | 劉澤琴 米夢(mèng)凡 楊秀峰 趙子輝 解春苓

        1. 北華航天工業(yè)學(xué)院經(jīng)濟(jì)管理學(xué)院 2. 北華航天工業(yè)學(xué)院遙感信息工程學(xué)院 3. 河北省航天遙感信息處理與應(yīng)用協(xié)同創(chuàng)新中心

        一、前言

        在黨的組織結(jié)構(gòu)和國(guó)家政權(quán)結(jié)構(gòu)中,縣一級(jí)處在承上啟下的關(guān)鍵環(huán)節(jié),是發(fā)展經(jīng)濟(jì)、保障民生、維護(hù)穩(wěn)定的重要基礎(chǔ)。縣域經(jīng)濟(jì)帶有城鎮(zhèn)經(jīng)濟(jì)和農(nóng)村經(jīng)濟(jì)相結(jié)合的典型特征,在城鎮(zhèn)化發(fā)展進(jìn)程中的土地資源合理利用是具有重要研究?jī)r(jià)值的關(guān)鍵問(wèn)題之一。如何利用遙感、GIS等現(xiàn)代技術(shù)手段,對(duì)縣域土地資源進(jìn)行定性分析、精準(zhǔn)施策、科學(xué)管理,提升土地利用的效益,是縣、鄉(xiāng)政府等決策部門亟待思考的問(wèn)題。

        在區(qū)域經(jīng)濟(jì)研究中引入空間地理信息變量是近年來(lái)獲得迅速發(fā)展的經(jīng)濟(jì)地理學(xué)的一個(gè)熱門研究領(lǐng)域。作為經(jīng)濟(jì)學(xué)與遙感學(xué)、地理學(xué)相互融合的一門新興學(xué)科,經(jīng)濟(jì)地理學(xué)在產(chǎn)業(yè)經(jīng)濟(jì)和區(qū)域經(jīng)濟(jì)研究領(lǐng)域受到了越來(lái)越高的重視,其中土地利用與覆被變化(LUCC)屬于熱點(diǎn)領(lǐng)域。這類包括空間變量的模型主要特點(diǎn)是采用遙感信息作為數(shù)據(jù)來(lái)源,并且利用計(jì)算機(jī)技術(shù)和GIS,分析LUCC的空間顯性分布特征和格局構(gòu)成,同時(shí)重視對(duì)影響空間差異性的自然、人類活動(dòng)等驅(qū)動(dòng)因子的分析。目前廣泛采用的分析模型在原理構(gòu)造上大都源于地理元胞自動(dòng)機(jī)(CA)。CLUE-S是源于CA原理的一種比較成熟的模型形式。區(qū)域經(jīng)濟(jì)研究的另一個(gè)熱點(diǎn)領(lǐng)域是土地利用的效益評(píng)價(jià),涉及生態(tài)(環(huán)境)效益、經(jīng)濟(jì)效益和社會(huì)效益等方面。

        空間地理信息分布及格局演化方面的技術(shù)方法目前已取得較大進(jìn)展,空間地理模型分析早已有之,但是將其研究結(jié)果應(yīng)用于地方政府治理和民間應(yīng)用的案例不多,應(yīng)用于較小區(qū)域開展經(jīng)濟(jì)分析的研究數(shù)量更少,尤其是針對(duì)縣域這類國(guó)民經(jīng)濟(jì)基本單元的實(shí)證研究比較鮮見(jiàn)??h域經(jīng)濟(jì)所處的地位及其結(jié)構(gòu)具有重要研究?jī)r(jià)值,本文力圖在此方面做出探索。

        二、研究區(qū)、數(shù)據(jù)及研究方法

        1. 研究區(qū)概況

        巨鹿縣位于河北省中南部,邢臺(tái)市中部,地處太行山東麓,華北平原的邊緣,古黃河、漳河沖積平原上,地勢(shì)平坦。地理位置為37°N和115°E,四周與邢臺(tái)市7個(gè)縣交界,縣域面積631km。該縣總?cè)丝?2萬(wàn),下轄6個(gè)鎮(zhèn),4個(gè)鄉(xiāng),縣城駐地巨鹿鎮(zhèn)。巨鹿縣交通便捷,公路總里程逾1400km,邢衡高速貫穿全境,國(guó)家二級(jí)公路(省道)南郝線、邢德線橫穿東西,定魏線縱貫?zāi)媳?。?dāng)?shù)剞r(nóng)業(yè)發(fā)展具有特色,是中國(guó)金銀花道地藥材之鄉(xiāng)、中國(guó)枸杞產(chǎn)業(yè)基地、中國(guó)串枝紅杏之鄉(xiāng),中草藥材和園林水果種植兩項(xiàng)產(chǎn)值在農(nóng)業(yè)總產(chǎn)值中占比近1/4。面對(duì)土地資源以耕地為主、城鎮(zhèn)化帶來(lái)的建設(shè)用地?cái)U(kuò)張、新興的二、三產(chǎn)業(yè)發(fā)展對(duì)土地需求增加等土地利用的多重矛盾,如何實(shí)現(xiàn)土地節(jié)約、集約利用,提升土地利用的經(jīng)濟(jì)效益,是當(dāng)?shù)卣叫杵平獾碾y題。

        2. 數(shù)據(jù)

        本文基礎(chǔ)數(shù)據(jù)由衛(wèi)星遙感信息和經(jīng)濟(jì)社會(huì)發(fā)展統(tǒng)計(jì)資料兩部分組成。

        遙感數(shù)據(jù)用于提取土地利用現(xiàn)狀和交通數(shù)據(jù),根據(jù)研究目標(biāo)選擇了2010、2014和2018年三期遙感影像,其中2010年遙感數(shù)據(jù)來(lái)源于谷歌影像,數(shù)據(jù)精度7.36m;2014年和2018年遙感影像來(lái)源于高分一號(hào)衛(wèi)星,數(shù)據(jù)精度2m。影像時(shí)間均選取當(dāng)年7-8月份,數(shù)據(jù)質(zhì)量好。土地利用分類處理采用ArcGIS軟件通過(guò)人工解譯進(jìn)行矢量化,土地利用類型包括七種,分別為耕地、林地、城鎮(zhèn)建設(shè)用地、農(nóng)村居民點(diǎn)、水域、工礦企業(yè)用地和其他用地(主要是交通設(shè)施和未利用地)(表1)。圖1為150m×150m尺度的三期土地利用情況??h域內(nèi)交通設(shè)施包括省道、縣道、鄉(xiāng)村道路三類道路(線狀地理信息)(2014年10月邢衡高速公路邢臺(tái)段通車,鑒于本文所涉時(shí)間范圍自2010年開始,因此未將該高速公路和巨鹿縣高速站口作為交通因素納入研究范圍),以及縣客運(yùn)總站(點(diǎn)狀地理信息)。利用歐式距離法計(jì)算縣域內(nèi)各位置至各交通設(shè)施的距離。

        表1 土地利用類型及代碼

        圖1 巨鹿縣2010年、2014年、2018年土地利用狀況

        經(jīng)濟(jì)社會(huì)發(fā)展數(shù)據(jù)來(lái)源于2011-2018年《河北經(jīng)濟(jì)年鑒》《邢臺(tái)統(tǒng)計(jì)年鑒》和《巨鹿統(tǒng)計(jì)月報(bào)》,涉及國(guó)民經(jīng)濟(jì)核算、農(nóng)村經(jīng)濟(jì)、工業(yè)等部分,所有指標(biāo)均采用現(xiàn)價(jià)(當(dāng)年價(jià))口徑。在對(duì)2022年交通因素的變化推導(dǎo)中,參考了《巨鹿縣城鄉(xiāng)總體規(guī)劃(2013-2030)》,《巨鹿縣國(guó)民經(jīng)濟(jì)與社會(huì)發(fā)展第十三個(gè)五年規(guī)劃綱要》,以及巨鹿縣近年《政府工作報(bào)告》。

        3. 研究方法

        (1)土地利用變化模型CLUE-S

        CLUE-S模型是由CLUE模型演化而來(lái)的,更加適于開展較小區(qū)域的研究。對(duì)于CLUE-S模型的擬合效果檢驗(yàn)多采用Kappa系數(shù),一般認(rèn)為該系數(shù)值≥0.75時(shí),表示模型擬合效果較好。

        (2)經(jīng)濟(jì)效益評(píng)價(jià)模型

        土地利用經(jīng)濟(jì)效益評(píng)價(jià)主要利用研究區(qū)經(jīng)濟(jì)社會(huì)發(fā)展統(tǒng)計(jì)數(shù)據(jù),如地區(qū)生產(chǎn)總值及其構(gòu)成、主要產(chǎn)業(yè)產(chǎn)值等,結(jié)合各種土地利用類型的用地量數(shù)據(jù)進(jìn)行簡(jiǎn)單匯總。

        三、結(jié)果與分析

        1. 2010-2018年巨鹿縣土地利用模擬

        (1)Logistic回歸分析結(jié)果

        CLUE-S模型中,柵格大小一般應(yīng)在100~500m之間,結(jié)合研究區(qū)的規(guī)模和經(jīng)濟(jì)特征,本文采用150m×150m的柵格進(jìn)行分析。采用SPSS 22中的Binary Logistic回歸模型研究2010-2014年間交通因素對(duì)土地利用類型的驅(qū)動(dòng)關(guān)系,分析得到回歸系數(shù)。就結(jié)果來(lái)看,僅有農(nóng)村居民點(diǎn)的模擬效果欠佳,其他用地類型的ROC曲線下面積均超過(guò)0.6。因此Logistic回歸效果基本可以滿足當(dāng)前研究需要,可將其作為CLUE-S模型分析的位置布局參數(shù)。

        (2)地類間轉(zhuǎn)換規(guī)則設(shè)置

        針對(duì)巨鹿縣2022年的土地利用情況模擬預(yù)測(cè),本文擬作三種不同的假設(shè)情境設(shè)定,分別為:“自然發(fā)展”情境,即假設(shè)當(dāng)?shù)乩^續(xù)保持以往的發(fā)展方式,不對(duì)土地利用施加人為干預(yù);“耕地保護(hù)”情境,即假設(shè)當(dāng)?shù)氐母乜偭勘3植蛔?,并令已有耕地向其他用地類型轉(zhuǎn)換的難度加大,以保持和充分發(fā)揮巨鹿縣在特色農(nóng)產(chǎn)品種植業(yè)方面的優(yōu)勢(shì);“水域保護(hù)”情境,即假設(shè)當(dāng)?shù)厮蛎娣e保持不變,并令已有耕地向其他用地類型轉(zhuǎn)換的難度加大,以適應(yīng)巨鹿縣“千湖之縣”發(fā)展思路的需要。

        土地利用轉(zhuǎn)移彈性用于表示某種用地類型轉(zhuǎn)換為其他用地類型的難易程度,取值為0~1,其中0表示可轉(zhuǎn)換為任意其他地類,1表示不允許轉(zhuǎn)換為其他地類,該值越大表明轉(zhuǎn)換成其他用地類型的難度越高。針對(duì)這三種假設(shè)情境,結(jié)合各種用地類型自身的特征,分別設(shè)置七種地類的土地利用轉(zhuǎn)移彈性。土地利用轉(zhuǎn)移矩陣用于表示某種用地類型轉(zhuǎn)換為其他用地類型的可能性,取值為0或1,其中0表示不可轉(zhuǎn)換,1表示可以轉(zhuǎn)換。

        (3)基于CLUE-S的2018年土地利用空間模擬

        為了檢驗(yàn)上述參數(shù)設(shè)置的合理性,本文以2018年巨鹿縣遙感影像解譯得到的150m×150m尺度的土地利用情況為目標(biāo)時(shí)期的土地利用需求,見(jiàn)圖1(c),2014年的相應(yīng)數(shù)據(jù)作為基期,通過(guò)線性插值法確定中間年份各種地類的用地量,采用CLUE-S模擬“自然發(fā)展”情境下2018年的土地利用情況(鑒于巨鹿縣當(dāng)?shù)匚丛O(shè)有自然保護(hù)區(qū)或其他明確限制開發(fā)區(qū)域,因此本文在研究中未設(shè)置空間政策與約束條件),得到模擬正確的柵格數(shù)為25544個(gè),與2018年150m×150m尺度下總柵格數(shù)28019個(gè)相比,得到Kappa系數(shù)為91.17%,表明CLUE-S模型模擬效果好,以上參數(shù)設(shè)置合理,可以將其代入2018-2022年分析周期進(jìn)行2022年土地利用情況的預(yù)測(cè)。

        2. 不同假設(shè)情境下2022年巨鹿縣土地利用預(yù)測(cè)

        (1)2022年巨鹿縣土地利用需求預(yù)測(cè)

        對(duì)土地利用需求的預(yù)測(cè)可以采用簡(jiǎn)單趨勢(shì)外推法,也可以采用形式更為復(fù)雜的經(jīng)濟(jì)模型進(jìn)行分析。本文重點(diǎn)關(guān)注了在縣域經(jīng)濟(jì)發(fā)展中貢獻(xiàn)顯著、同時(shí)用地?cái)?shù)量變化較大的第二產(chǎn)業(yè)發(fā)展情況,以第二產(chǎn)業(yè)增加值作為其發(fā)展水平的測(cè)算指標(biāo)。接下來(lái)本文采用巨鹿縣第二產(chǎn)業(yè)增加值(X)和工礦企業(yè)用地面積(Y)進(jìn)行趨勢(shì)擬合的方式預(yù)測(cè)2022年該土地利用類型的土地需求量。此外,在“耕地保護(hù)”情境下假設(shè)耕地面積保持與2018年相同,在“水域保護(hù)”情境下假設(shè)水域面積保持與2018年相同。在保持土地總量不變的前提下,其他用地類型的土地需求量繼續(xù)延續(xù)原有趨勢(shì)并作適當(dāng)調(diào)整。

        (2)2022年巨鹿縣土地利用預(yù)測(cè)

        采用交通因素驅(qū)動(dòng)下土地利用類型的Logistic回歸系數(shù)和土地利用轉(zhuǎn)移矩陣,根據(jù)三種假設(shè)情境,分別采用不同的土地利用轉(zhuǎn)移彈性結(jié)合土地利用需求量,以2018年土地利用狀況為基礎(chǔ),運(yùn)行CLUE-S模型得到不同情境下的2022年巨鹿縣土地利用預(yù)測(cè),見(jiàn)圖2。

        (3)土地利用預(yù)測(cè)結(jié)果對(duì)比分析

        圖2 巨鹿縣2022年土地利用預(yù)測(cè)

        根據(jù)不同假設(shè)情境開展的2022年巨鹿縣土地利用預(yù)測(cè)結(jié)果顯示,在以下兩個(gè)方面存在比較顯著的差異。

        土地利用結(jié)構(gòu)的數(shù)量差異??v觀2010年以來(lái)的土地利用變化,可以發(fā)現(xiàn)耕地?cái)?shù)量變化的絕對(duì)值最大,對(duì)林地、工礦企業(yè)用地等用地量的增加貢獻(xiàn)最顯著。而在土地需求預(yù)測(cè)中,本文對(duì)工礦企業(yè)用地基于經(jīng)濟(jì)模型做出了數(shù)量規(guī)定,相當(dāng)于設(shè)置了較具剛性的土地需求條件。由于“耕地保護(hù)”情境假設(shè)耕地?cái)?shù)量維持在2018年現(xiàn)狀水平,直接導(dǎo)致的結(jié)果是林地?zé)o法有效增加,因此圖2(b)顯示的結(jié)果與圖1(c)差異最小。在“水域保護(hù)”情境下的情況類似,不過(guò)由于水域面積在土地總量中占比很小(2018年為5.85%),對(duì)其施加的數(shù)量限制條件影響較小,因此該情境與“自然發(fā)展”情境的土地利用需求相似度較高。

        土地利用格局的空間差異?!白匀话l(fā)展”和“水域保護(hù)”兩種情境與2018年土地利用現(xiàn)狀的差異最為顯著,共同的變化特征有兩個(gè):要求將目前縣城北部和東部周邊尤其是東北方向的一些耕地轉(zhuǎn)為林地,形成縣城以南和以北兩個(gè)林木富集區(qū)域;要求在目前縣城西北部和東北部部分地區(qū)安排工礦企業(yè)用地?cái)U(kuò)大,形成連片規(guī)模發(fā)展。此外,三種假設(shè)情境都要求將西部的西郭城鎮(zhèn)靠近與臨縣交界區(qū)(臨近省道S324邢德線)、蘇家營(yíng)鄉(xiāng)與臨縣交界區(qū)和官亭鎮(zhèn)臨近省道S327南郝線部分區(qū)域的農(nóng)村居民點(diǎn)轉(zhuǎn)為工礦企業(yè)用地,充分發(fā)揮交通設(shè)施助力企業(yè)發(fā)展的作用。

        3. 巨鹿縣土地利用經(jīng)濟(jì)效益評(píng)價(jià)

        考慮到各種土地利用類型與國(guó)民經(jīng)濟(jì)產(chǎn)業(yè)的對(duì)應(yīng)性,以及縣域統(tǒng)計(jì)資料的可得性,本文選取(狹義)農(nóng)業(yè)產(chǎn)值作為耕地的產(chǎn)出指標(biāo),林業(yè)產(chǎn)值、第三產(chǎn)業(yè)增加值、農(nóng)林牧漁服務(wù)業(yè)產(chǎn)值、漁業(yè)產(chǎn)值和規(guī)模以上工業(yè)總產(chǎn)值分別作為林地、城鎮(zhèn)建設(shè)用地、農(nóng)村居民點(diǎn)、水域和工礦建設(shè)用地的產(chǎn)出指標(biāo)。所有指標(biāo)均采用現(xiàn)價(jià)(當(dāng)年價(jià))口徑統(tǒng)計(jì)數(shù)據(jù)。根據(jù)2010-2018年統(tǒng)計(jì)資料和相應(yīng)地類的用地量,計(jì)算得到各地類的單位面積經(jīng)濟(jì)效益值,根據(jù)2010-2018年各地類單位面積經(jīng)濟(jì)效益值,通過(guò)趨勢(shì)外推法得到2022年各指標(biāo)預(yù)測(cè)值。以預(yù)測(cè)的2022年單位面積經(jīng)濟(jì)效益值分別結(jié)合三種假設(shè)情境下的各地類用地量,可以計(jì)算得到不同情境下的經(jīng)濟(jì)效益總量預(yù)測(cè)值(表2)。

        表2 土地利用經(jīng)濟(jì)效益總量 單位:萬(wàn)元

        由表2顯示的三種不同假設(shè)情境下巨鹿縣土地利用經(jīng)濟(jì)效益總量來(lái)看,“自然發(fā)展”情境的經(jīng)濟(jì)效益顯著高于其他兩種情境,這與該情境下單位面積經(jīng)濟(jì)效益值最高的工礦企業(yè)用地比其他二者多有密切關(guān)系。其他兩種假設(shè)情境之間的經(jīng)濟(jì)效益總量差異不顯著,雖然兩種情境下的各地類用地量存在差別,但是結(jié)合各種地類的單位面積經(jīng)濟(jì)效益值之后綜合來(lái)看,經(jīng)濟(jì)效果上是相當(dāng)?shù)?。如果僅從經(jīng)濟(jì)效益總量的角度出發(fā)考慮,當(dāng)?shù)厥紫葢?yīng)當(dāng)考慮選擇“自然發(fā)展”模式,但是如果補(bǔ)充更多的生態(tài)效益和社會(huì)效益資料,則綜合評(píng)價(jià)的結(jié)果可能會(huì)發(fā)生變化。

        四、巨鹿縣土地利用及產(chǎn)業(yè)發(fā)展策略啟示

        采用交通設(shè)施作為驅(qū)動(dòng)因素,運(yùn)用CLUE-S模型對(duì)2022年巨鹿縣土地利用情況進(jìn)行模擬預(yù)測(cè),對(duì)比2018年當(dāng)?shù)赝恋乩矛F(xiàn)狀發(fā)現(xiàn),工礦企業(yè)用地?cái)U(kuò)大和林地增加是預(yù)測(cè)中出現(xiàn)的主要特征,并且工礦企業(yè)的用地指向與當(dāng)?shù)爻青l(xiāng)總體規(guī)劃中指出的未來(lái)中長(zhǎng)期發(fā)展思路是一致的,但是用地類型的布點(diǎn)更加具體,體現(xiàn)了科學(xué)性和應(yīng)用性的結(jié)合。

        研究發(fā)現(xiàn),在“自然發(fā)展”和“水域保護(hù)”兩種假設(shè)情境下,應(yīng)在目前縣城東北方向成規(guī)模地建設(shè)林地,這與2010-2018年間巨鹿縣林地面積顯著增加有關(guān),也與模型的驅(qū)動(dòng)因素、限制條件等有關(guān)。這一發(fā)展方向與當(dāng)?shù)亟ㄔO(shè)“華北地區(qū)最綠最美平原縣”的發(fā)展目標(biāo)一致,也提示了林地建設(shè)的布局方向。但是就當(dāng)前林業(yè)產(chǎn)值較低的現(xiàn)狀來(lái)看,圍繞林業(yè)資源開發(fā)休閑旅游等第三產(chǎn)業(yè)是可以考慮的發(fā)展策略。

        “耕地保護(hù)”和“水域保護(hù)”兩種假設(shè)情境下,經(jīng)濟(jì)效益總量與“自然發(fā)展”情境存在一定差距,這固然與評(píng)價(jià)模型選擇的經(jīng)濟(jì)指標(biāo)有關(guān),但是也表明將種植業(yè)、漁業(yè)等囿于以往的第一產(chǎn)業(yè)部門是不夠的,應(yīng)該充分挖掘相應(yīng)自然資源、經(jīng)濟(jì)資源的潛力,將第一產(chǎn)業(yè)與工業(yè)、更重要地是與第三產(chǎn)業(yè)充分融合,將巨鹿縣道地中藥材產(chǎn)地的優(yōu)勢(shì)、河湖坑塘密集的優(yōu)勢(shì)充分發(fā)揮出來(lái),走出一條特色農(nóng)業(yè)發(fā)展之路。

        五、結(jié)論

        采用空間地理信息尤其是衛(wèi)星遙感信息與經(jīng)濟(jì)效益分析模型相結(jié)合的方法,對(duì)縣域土地利用及其效益開展模擬和評(píng)價(jià),這一思路具有可行性。這種研究方法突破了以往大量經(jīng)濟(jì)分析中僅僅依靠經(jīng)濟(jì)社會(huì)統(tǒng)計(jì)信息作為數(shù)據(jù)源的思路,將空間地理信息引入分析框架,可以充分發(fā)揮此類信息實(shí)時(shí)性、翔實(shí)性、動(dòng)態(tài)性等優(yōu)勢(shì),不僅實(shí)現(xiàn)了基礎(chǔ)數(shù)據(jù)的豐富性,更有利于提升經(jīng)濟(jì)分析的質(zhì)量。

        在本文所述新的研究框架下,模型分析所提供的結(jié)果更加豐富。根據(jù)不同的假設(shè)情境,可以模擬未來(lái)土地利用的多種格局,并且以圖形的形式呈現(xiàn)出來(lái),結(jié)果的直觀性更好,比以往常用的經(jīng)濟(jì)研究方法提供的結(jié)果更易于理解。結(jié)合經(jīng)濟(jì)效益分析之后,各種不同假設(shè)情境下的經(jīng)濟(jì)效益及其所對(duì)應(yīng)的產(chǎn)業(yè)發(fā)展策略更加明晰,政策指導(dǎo)性更強(qiáng)。

        根據(jù)本文所采用研究方法所提供的研究結(jié)果,對(duì)于地方政府調(diào)整土地利用格局、制定地方經(jīng)濟(jì)社會(huì)發(fā)展規(guī)劃、調(diào)整產(chǎn)業(yè)發(fā)展政策等均具有良好的參考價(jià)值。此類分析的結(jié)果明確、直觀、指向性強(qiáng),易于轉(zhuǎn)化為具體的政策措施,且研究所耗費(fèi)的時(shí)間、經(jīng)濟(jì)成本均較低,可以提高政府治理效率。

        欧美三级不卡在线观看| 久久精品国产亚洲av热九| 日本女u久久精品视频| 日本久久伊人特级黄色| s级爆乳玩具酱国产vip皮裤| 亚洲欧美日韩精品高清| 久久精品国产亚洲av成人擦边| 各类熟女熟妇激情自拍| 午夜爽爽爽男女污污污网站| 亚洲综合久久成人a片| 最新欧美一级视频| 美女被内射中出在线观看| 国产精品h片在线播放| 久久久亚洲色| 亚洲无线码一区在线观看| 高清中文字幕一区二区三区| 欧美日韩在线视频一区| 国产山东熟女48嗷嗷叫| 日本老年人精品久久中文字幕| 中文资源在线一区二区三区av| 永久免费人禽av在线观看| 国产一区二区三区四区五区vm| 日本护士一区二区三区高清热线 | 高潮内射主播自拍一区| 国产女人的高潮国语对白| 操B小视频国产| 亚洲熟女av一区少妇| 亚洲av永久无码精品网站| 久久婷婷成人综合色| 亚洲公开免费在线视频| aa日韩免费精品视频一| 狼狼综合久久久久综合网| 国产精品久久久久免费a∨| 日本精品国产1区2区3区| 精品国产一区二区三区三级| 欧美怡红院免费全部视频| 中文字幕在线观看国产双飞高清| 国产黄色三级三级三级看三级| 欧洲女人与公拘交酡视频| 久久免费看少妇高潮v片特黄| 插入中文字幕在线一区二区三区 |