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

        ?

        巢湖流域土地整治與面源污染時(shí)空過程及關(guān)系

        2021-01-12 03:28:38張淼淼徐建飛李素萃呂雪嬌趙彬程
        水土保持研究 2021年1期
        關(guān)鍵詞:巢湖面源整治

        張淼淼, 肖 武, 徐建飛, 李素萃, 呂雪嬌, 趙彬程

        (1.中國礦業(yè)大學(xué)(北京) 土地復(fù)墾與生態(tài)重建研究所, 北京 100083;2.浙江大學(xué) 公共管理學(xué)院, 杭州 310058; 3.中國資源衛(wèi)星應(yīng)用中心, 北京 100083)

        中國土地整治開始于20世紀(jì)90年代末,經(jīng)歷了“重?cái)?shù)量”、“數(shù)量與質(zhì)量并重”和“數(shù)量、質(zhì)量與生態(tài)并重”的3個(gè)階段[1]。土地整治的作用不僅僅是提高土地利用效率,促進(jìn)農(nóng)業(yè)生產(chǎn)活動(dòng),它更提升了一個(gè)維度的內(nèi)涵[2],不僅要通過“田水路林村鎮(zhèn)”改善人類的生存條件,還要通過“山水林田湖草”整治以實(shí)現(xiàn)生態(tài)化國土空間,即建設(shè)生態(tài)良田、修復(fù)受損土地和改善環(huán)境為中心的生態(tài)景觀型的土地整治已成為主線[3]。

        巢湖地處江淮,毗鄰長江,是安徽省長江流域最大的河流系統(tǒng),也是國家糧食安全戰(zhàn)略布局的重要區(qū)域。農(nóng)田安全是糧食安全的重要保障,而農(nóng)田水環(huán)境是保障農(nóng)田安全的關(guān)鍵環(huán)節(jié)。隨著農(nóng)業(yè)經(jīng)濟(jì)飛速發(fā)展,糧食大幅增產(chǎn),面源污染問題成為巢湖流域水環(huán)境治理的焦點(diǎn),從1978年首次發(fā)生富營養(yǎng)化,水質(zhì)不斷惡化,到2003年污染狀態(tài)峰值,逐漸的引起政府及專家學(xué)者的重視。目前,中國農(nóng)業(yè)面源污染主要來源于化肥施用、農(nóng)田固體廢棄物、畜禽養(yǎng)殖糞污、水產(chǎn)養(yǎng)殖污水和農(nóng)村生活污染等5個(gè)方面[4]。相關(guān)研究表明農(nóng)業(yè)布局調(diào)整措施[5],防護(hù)林種植措施[6],改善農(nóng)業(yè)種植結(jié)構(gòu)措施[7],生態(tài)排灌措施等[8]土地整治措施會(huì)對面源污染產(chǎn)生間接影響,以面源污染治理為中心的生態(tài)景觀型的土地整治是實(shí)現(xiàn)巢湖流域生態(tài)農(nóng)田建設(shè)的必經(jīng)之路。

        然而,從區(qū)域角度來看,土地整治多是點(diǎn)狀分布,且以增加耕地面積與提高耕地質(zhì)量為核心,對面源污染境與生態(tài)環(huán)境的改善無顯著作用,面源污染治理往往更側(cè)重于農(nóng)村污染防控[9]、農(nóng)用地化肥施用的控制[10],以及水體的治理等[11]方面。本研究以巢湖為例,基于SWAT模型模擬的2000—2010年面源污染負(fù)荷(總氮、總磷)和2000—2010年土地整治項(xiàng)目數(shù)據(jù),結(jié)合GWR模型,從宏觀角度探究區(qū)域尺度的土地整治項(xiàng)目實(shí)施與面源污染的空間關(guān)系,分析土地整治對面源污染的作用機(jī)制,為巢湖以及同類地區(qū)以面源污染治理為中心的生態(tài)景觀型土地綜合整治提供科學(xué)依據(jù)及建議。

        1 研究區(qū)概況與研究方法

        1.1 研究區(qū)概況

        巢湖是我國第五大淡水湖,位于安徽省中部,是全國重要的水源地和生態(tài)保護(hù)區(qū)。地理坐標(biāo)為116°25′30″—118°29′52″E,30°5′18″—32°32′30″N。巢湖北有柘皋河、店埠河入湖,西有南淝河、派河、豐樂河、杭埠河、白石天河流入,南與兆河、白湖相通,東經(jīng)裕溪河下泄長江,水資源豐富。巢湖地處江淮丘陵地帶,地勢西高東低、中間低洼平坦,且氣候溫和濕潤,年平均氣溫15~16℃,年平均降水量為1 100 mm,夏季降雨居多。土壤類型以黃棕壤為主,土地利用類型以耕地、林地、建設(shè)用地與水域?yàn)橹鳌?/p>

        考慮到巢湖流域土地整治項(xiàng)目以完整行政區(qū)為單位,本研究選取巢湖流域所在的14個(gè)縣(市、區(qū))為研究區(qū)域,即無為縣、和縣、含山縣、舒城縣、金安區(qū)、長豐縣、瑤海區(qū)、蜀山區(qū)、廬陽區(qū)、廬江區(qū)、肥西縣、肥東縣、巢湖市、包河區(qū),總面積約為1.97萬km2,見圖1。

        圖1 研究區(qū)位置

        1.2 數(shù)據(jù)來源及數(shù)據(jù)庫的建立

        本文參考以往研究[12-13]和多次試驗(yàn)結(jié)果,最終確定適合本研究單元格網(wǎng)大小為1 km×1 km,以此單元格網(wǎng)為單位進(jìn)行數(shù)據(jù)庫建立及結(jié)果與分析,并選擇3個(gè)時(shí)間點(diǎn)即2000年、2005年、2010年進(jìn)行節(jié)點(diǎn)性和階段性的研究。

        1.2.1 面源污染數(shù)據(jù) 巢湖流域中,面源污染的各方面研究已經(jīng)很成熟[14-15],SWAT模型也成為了面源污染問題研究的重要手段。SWAT(Soil and Water Assessment Tool)模型是20世紀(jì)90年代初美國(USDA)開發(fā)的分布式農(nóng)業(yè)面源污染負(fù)荷模型[16],可根據(jù)不同的土地利用、土壤以及水利措施對流域的水、泥沙、化學(xué)物質(zhì)等進(jìn)行模擬。

        SWAT模型所需要的數(shù)據(jù)分為兩類,空間數(shù)據(jù)和屬性數(shù)據(jù),采用阿爾伯斯雙標(biāo)準(zhǔn)緯線多圓錐投影(Asia_North_Albers_Equal_Area_Conic),具體情況見表1。

        表1 面源污染數(shù)據(jù)情況

        為了后續(xù)的定量分析和模型模擬,將SWAT模型模擬的巢湖流域總氮和總磷負(fù)荷,按照國家《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》修訂[19]分為6個(gè)級別進(jìn)行空間展示,劃分等級見表2。

        1.2.2 土地整治數(shù)據(jù) 巢湖流域土地整治指標(biāo)構(gòu)建的基礎(chǔ)數(shù)據(jù)來自于2000—2010年研究區(qū)范圍內(nèi)的土地整治項(xiàng)目。2000—2010年巢湖土地整治項(xiàng)目主要分布于巢湖市、肥東縣、肥西縣、廬江縣、長豐縣、無為縣和金安區(qū)等14個(gè)區(qū)縣,共2 609個(gè),2000—2010年的實(shí)際總投資達(dá)到9.48億,新增耕地1.43萬hm2,實(shí)際建設(shè)總規(guī)模達(dá)到3.12萬hm2??紤]數(shù)據(jù)本身特性以及土地整治項(xiàng)目驗(yàn)收規(guī)程的綜合性原則,借鑒以往土地整治指標(biāo)的構(gòu)建[20-21],選取了整治項(xiàng)目強(qiáng)度、建設(shè)強(qiáng)度、新增耕地、總投資以及投資強(qiáng)度5個(gè)基礎(chǔ)指標(biāo)。5個(gè)指標(biāo)均可作為土地整治力度的代表,故采用同等權(quán)重法綜合形成整治因子指標(biāo),代表2000—2010年的整治力度,詳情見表3,結(jié)果見圖2。

        表2 面源污染劃分等級

        表3 土地整治因子計(jì)算

        圖2 整治因子及整治項(xiàng)目分布

        1.3 統(tǒng)計(jì)方法

        本文首先利用全局Moran′sI和局部空間關(guān)聯(lián)指數(shù)(LISA)計(jì)算了2000年、2005年、2010年及2000—2005年、2005—2010年兩個(gè)時(shí)段面源污染全局空間自相關(guān)指數(shù)Moran′sI(顯著性水平p=0.001)。通過分析巢湖面源污染的空間自相關(guān)性及空間異質(zhì)性,探測其熱點(diǎn)區(qū)域,再結(jié)合地理加權(quán)回歸模型(GWR)進(jìn)行土地整治對2000年、2005年、2010年面源污染現(xiàn)狀及土地整治對2000—2005年,2005—2010年面源污染變化的影響研究。

        1.3.1 空間自相關(guān)分析 空間自相關(guān)是揭示空間數(shù)據(jù)分布的一個(gè)重要的方法,是測度某位置上的數(shù)據(jù)與其空間位置數(shù)據(jù)間的相互依賴程度,分為全域型和區(qū)域型兩種。該方法已經(jīng)廣泛應(yīng)用于空間數(shù)據(jù)挖掘,探索自然與社會(huì)現(xiàn)象的空間模式和分布。

        空間自相關(guān)的Moran′sI統(tǒng)計(jì)可表示為:

        (1)

        1.3.2 地理加權(quán)回歸模型(GWR) 地理加權(quán)回歸模型(Geographically weighted regression,GWR)[22],是普通線性回歸(OLS)的改進(jìn),將數(shù)據(jù)的地理位置帶入到回歸參數(shù)中,并考慮相鄰點(diǎn)的空間權(quán)重,允許局部的參數(shù)估計(jì)。已被廣泛地應(yīng)用到經(jīng)濟(jì)、自然資源管理、醫(yī)學(xué)多個(gè)等研究領(lǐng)域,近年來也逐漸成為水質(zhì)研究領(lǐng)域的新熱點(diǎn)。GWR模型計(jì)算公式為:

        (2)

        式中:(ui,vi)為第i點(diǎn)的空間位置;βk(ui,vi)為第i點(diǎn)的第k個(gè)回歸系數(shù);εi為第i點(diǎn)的隨機(jī)誤差。

        2 結(jié)果與分析

        2.1 面源污染分布及演變結(jié)果與分析

        面源污染變化趨勢分為3類:面源污染無變化區(qū)、面源污染惡化區(qū)和面源污染恢復(fù)區(qū)。圖3為2000—2010年面源污染(總氮、總磷)模擬結(jié)果,表4—7為面源污染風(fēng)險(xiǎn)轉(zhuǎn)移矩陣。由圖3可看出巢湖流域北部長豐縣、肥東縣、瑤海區(qū)、蜀山區(qū)、廬陽區(qū)、包河區(qū)在2000—2010年總氮等級下降較劇烈,而總磷等級在瑤海區(qū)、蜀山區(qū)、廬陽區(qū)、包河區(qū)下降趨勢較明顯。

        圖3 2000-2010年巢湖流域總氮和總磷負(fù)荷的空間分布

        從表4—7中可以看出,2000—2005年和2005—2010年兩個(gè)階段期間,總氮等級無變化區(qū)的比重從92.85%,下降至91.65%,總磷等級無變化區(qū)比重也從95.06%下降至94.91%,說明2005—2010年階段巢湖區(qū)域總氮和總磷變化更為劇烈。

        表4 2000-2005年面源污染(總氮)風(fēng)險(xiǎn)轉(zhuǎn)移矩陣

        表5 2005-2010年面源污染(總氮)風(fēng)險(xiǎn)轉(zhuǎn)移矩陣

        表6 2000-2005年面源污染(總磷)風(fēng)險(xiǎn)轉(zhuǎn)移矩陣

        表7 2005-2010年面源污染(總磷)風(fēng)險(xiǎn)轉(zhuǎn)移矩陣

        巢湖面源污染變化是惡化與恢復(fù)相互交織的過程,從變化區(qū)的面積來看,2005—2010年總氮恢復(fù)區(qū)所占面積比重為6.88%較2000—2005年上升了0.2%,總氮惡化區(qū)的比例也從2000—2005年的0.49%增加至2005—2010年的1.46%??偭谆謴?fù)區(qū)所占比重從2000—2005年的4.73%下降至2005—2010年的4.33%,惡化區(qū)從2000—2005年的0.22%增加至2005—2010年的0.76%。但2000—2010年期間,縱觀巢湖流域面源污染(總氮、總磷)整體情況看恢復(fù)面積始終大于惡化面積。

        2.2 空間自相關(guān)結(jié)果與分析

        圖4所示為面源污染局部空間關(guān)聯(lián)指數(shù)(LISA)分類圖(顯著性水平p=0.001),由圖4可知總氮的高風(fēng)險(xiǎn)區(qū)(HH指標(biāo))在2000,2005,2010分別占總面積的22.37%,22.71%,22.36%,呈穩(wěn)定趨勢。總氮的低風(fēng)險(xiǎn)區(qū)(LL指標(biāo))在2000年、2005年、2010年分別占總面積的29.87%,30.67%,31.12%,呈上升趨勢。表明2000—2010年,巢湖流域內(nèi)總氮的面源污染情況呈現(xiàn)逐步恢復(fù)趨勢。圖5所示的總氮變化的熱點(diǎn)區(qū)域的統(tǒng)計(jì)結(jié)果也證實(shí)了該點(diǎn),總氮變化的HH型的比重在2000—2005年、2005—2010年兩個(gè)時(shí)期均小于LL型的區(qū)域,2005—2010總氮變化的LL型較2000—2005上升了3.28%。研究期間內(nèi),總磷也呈相同的趨勢。即面源污染各指標(biāo)存在較為顯著的正自相關(guān),總氮和總磷的空間分布表現(xiàn)出較為明顯的空間聚集。

        圖4 面源污染局部空間關(guān)聯(lián)指數(shù)(LISA)分類

        圖5 面源污染變化局部空間關(guān)聯(lián)指數(shù)(LISA)分類

        2.3 巢湖流域土地整治與面源污染的關(guān)系與響應(yīng)

        面源污染與土地整治的響應(yīng)關(guān)系進(jìn)研究主要包括面源污染對土地整治現(xiàn)狀的響應(yīng)及面源污染對土地整治變化的響應(yīng)兩部分。

        2.3.1 巢湖面源污染對土地整治現(xiàn)狀的響應(yīng) 巢湖面源污染對土地整治現(xiàn)狀的響應(yīng)是分別以2000年、2005年、2010年面源污染指數(shù)(總氮、總磷)為因變量,相應(yīng)年份的土地整治因子為自變量進(jìn)行分析。GWR的R2、回歸系數(shù)的空間分布見圖6、圖7和表8。

        表8 GWR模型R2

        圖6 2000-2010年現(xiàn)狀R2空間分布

        圖7 2000-2010年現(xiàn)狀回歸系數(shù)空間分布

        2000年、2005年、2010年的面源污染總氮R2從0.720升至0.730,總磷R(shí)2從0.717升至0.731,整治因子的擬合大致在0.71~0.74之間,說明整治因子大致解釋面源污染當(dāng)年現(xiàn)狀的71%~74%。這表明了在研究期間,研究區(qū)整治因子對面源污染現(xiàn)狀的影響起重要性作用,且逐漸增大的趨勢,也證明土地整治對面源污染的影響不容忽視。

        從面源污染與土地整治關(guān)系的空間分布看,面源污染對土地整治的響應(yīng)在空間上具有較強(qiáng)的異質(zhì)性。其R2較大的主要分為3部分在研究區(qū)西北、東南、東北中部呈斑塊狀分布(肥西縣、肥東縣、無為縣),然后以這些區(qū)域?yàn)橹行南蛲膺f減。相關(guān)系數(shù)較大區(qū)與土地整治集中區(qū)基本重合。從回歸系數(shù)看,2000年和2005年回歸系數(shù)大于0的區(qū)域仍主要分布于西北和東南兩部分,其余大部分系數(shù)處于小于0的狀態(tài);在2010年的回歸系數(shù)中,在西南、東南以及巢湖中部是回歸系數(shù)大于0的集中區(qū),并以此3個(gè)斑塊為中心向周遞減,且絕大部分區(qū)域系數(shù)小于0。2000年、2005年、2010年的回歸系數(shù)分布雖然存在著差異,但大部分區(qū)域系數(shù)處于小于0的狀態(tài),說明在增加土地整治強(qiáng)度的情況下,面源污染的值反而降低,即巢湖流域內(nèi)的面源污染得到了一定的改善。

        2.3.2 巢湖面源污染對土地整治變化的響應(yīng) 巢湖面源污染對土地整治變化的響應(yīng)是分別以2000—2005年、2005—2010年面源污染指數(shù)(總氮、總磷)為因變量,相應(yīng)階段的土地整治因子為自變量進(jìn)行模擬。GWR分析的R2為0.38~0.50,即土地整治的變化大致解釋了面源污染變化的38%~50%,這說明研究期間土地整治的實(shí)施在一定程度上對面源污染的分布造成影響,其變化趨勢影響了該區(qū)水生態(tài)環(huán)境變化的走向。

        2000—2005年R2在全區(qū)的分布基本上以中北部以及中南的塊狀高值區(qū)向外逐漸遞減,其余大部分區(qū)域以片狀分布,R2處于0~0.1的狀態(tài),所占比重最大。2005—2010年R2于0.3~0.5的區(qū)域大多分布于北部、西北和東南部(包河區(qū)、肥西縣、無為縣),并以此為中心向中部巢湖方向遞減,說明邊緣地區(qū)面源污染變化受土地整治影響較大。從回歸系數(shù)可以看出,研究區(qū)系數(shù)以小于0為主。

        3 討 論

        本研究發(fā)現(xiàn)隨著土地整治的進(jìn)行,面源污染有恢復(fù)的傾向。面源污染的緩解與多方面因素有關(guān),本文僅就土地整治角度,討論其對面源污染的作用機(jī)制:

        (1) 耕地細(xì)碎化降低,集約化農(nóng)業(yè)活動(dòng)降低化肥使用量。過量及不恰當(dāng)?shù)幕适┯檬窃斐擅嬖次廴镜闹匾颍埠饔蛲恋卣雾?xiàng)目多,投資大,全域是安徽省土地整治的重點(diǎn)區(qū)域,通過土地整治,打造了耕地集中連片分布的糧食主產(chǎn)區(qū)和基本農(nóng)田保護(hù)示范區(qū),耕地的細(xì)碎化程度大大降低,研究表明耕地細(xì)碎化對農(nóng)戶化肥使用效率具有顯著負(fù)影響[23],降低了面源污染的輸入源頭。

        (2) 四大土地整治工程作用。防護(hù)林工程可以有效的改善水土流失,一定程度上緩解生態(tài)保護(hù)壓力。土地平整工程中不僅可以田塊合并規(guī)整,還可通過坡改梯方式改變原有的坡面過程,形成田坎,改善生產(chǎn)耕作條件和水土保持能力。農(nóng)田水利工程也可通過建立生態(tài)的農(nóng)田灌排網(wǎng)絡(luò)系統(tǒng)和水利工程配套設(shè)施,合理布局,減少洪澇災(zāi)害帶來的水土流失。

        (3) 村莊集約規(guī)劃管理,加強(qiáng)監(jiān)管農(nóng)村生活污水的統(tǒng)一處理排放。加強(qiáng)村莊規(guī)劃管理,穩(wěn)步推進(jìn)村莊整治,引導(dǎo)農(nóng)民集中居住,有序搬遷撤并空心村和過于分散、生態(tài)惡劣的村莊。改善農(nóng)村生活環(huán)境,集中處理農(nóng)村生活污水廢水[24],有效減少面源污染外源。

        (4) 土地整治帶動(dòng)面源污染治理項(xiàng)目落實(shí),綜合治理效果佳。根據(jù)自然資源部土地整治中心發(fā)布的土地整治藍(lán)皮書,土地整治投資的帶動(dòng)效應(yīng)達(dá)到將近4.0,土地整治工作的開展提供了相較更加有利的耕作條件,也增加了農(nóng)業(yè)面源污染治理相關(guān)項(xiàng)目在土地整治區(qū)落地的可能性,其復(fù)合作用下的面源污染治理效果得到增強(qiáng)。

        4 結(jié)論與建議

        4.1 結(jié) 論

        (1) 巢湖流域面源污染(總氮、總磷)的空間分布主要呈現(xiàn)東西為軸線條帶狀以及南部區(qū)域高值,北部區(qū)域低值的空間分布。該空間分布特征是自然地理環(huán)境(西北部低山丘陵圍繞)與人類活動(dòng)(農(nóng)業(yè)耕作)共同作用的結(jié)果。

        (2) 面源污染恢復(fù)與惡化的情況并存,且存在空間差異?,幒^(qū)、蜀山區(qū)、廬陽區(qū)、包河區(qū)為合肥市4個(gè)區(qū),在土地整治的作用下仍能保持耕地斑塊形狀的規(guī)則,耕地面積小,化肥使用量小,2000—2010年總氮和總磷在此區(qū)域改善趨勢較明顯。除此之外,居巢區(qū)位于巢湖東側(cè),土地整治種植了大量防護(hù)林,改善了當(dāng)?shù)厮亮魇Лh(huán)境,故氮磷有所改善。

        (3) 巢湖流域面源污染(總氮和總磷)及面源污染變化區(qū)域全局自相關(guān)指數(shù)隨著年份增加,持續(xù)上升,影響范圍和程度越來越大,成為影響該區(qū)域環(huán)境問題的重要因素。

        (4) 面源污染對土地整治現(xiàn)狀的響應(yīng)R2在0.71~0.74,面源污染對土地整治變化的響應(yīng)R2在0.38~0.50,呈負(fù)相關(guān)。

        4.2 建 議

        (1) 生態(tài)田埂技術(shù),適當(dāng)提高田埂高度,可減少大部分的農(nóng)田地表徑流。因地制宜的在田埂的兩側(cè)栽種植物,形成隔離帶,有效吸附攔截氮磷養(yǎng)分和殘留農(nóng)藥。

        (2) 建設(shè)生態(tài)溝渠,在溝底、溝板種植植物,可過濾吸附農(nóng)田徑流污染物,達(dá)到控制污染物向水體遷移和氮磷養(yǎng)分再利用目的。

        (3) 推進(jìn)村莊集并,加強(qiáng)生活污水集中處理,避免出現(xiàn)由不良生活習(xí)慣和落后的排污體系直接排放到溝渠、坑塘、農(nóng)田等水體。

        (4) 土地平整過程調(diào)整農(nóng)田的土壤質(zhì)量,對相對貧瘠的土壤進(jìn)行改良,因地制宜科學(xué)規(guī)劃農(nóng)業(yè)產(chǎn)業(yè)種植結(jié)構(gòu)。

        (5) 類似于巢湖同類大型湖泊或河流區(qū)域,規(guī)劃種植防護(hù)林,改善土壤侵蝕和面源污染情況。

        本文以2000—2010年作為研究時(shí)間段,盡管在一定程度上揭示了土地整治與面源污染的響應(yīng)關(guān)系,但仍需探索更長時(shí)間序列下兩者的空間關(guān)系。生態(tài)型土地整治是實(shí)現(xiàn)巢湖流域生態(tài)農(nóng)田建設(shè)的必經(jīng)之路,本文為巢湖及類似共同考慮面源污染治理和農(nóng)田安全保障的區(qū)域提供土地整治的建議和依據(jù)。

        猜你喜歡
        巢湖面源整治
        篤行不怠 奮力深化提升三年專項(xiàng)整治
        農(nóng)業(yè)面源污染的危害與治理
        澄江市農(nóng)業(yè)面源污染成因及對策
        專項(xiàng)整治
        專項(xiàng)整治
        巢湖頌歌
        基于SWAT模型的漳河流域面源污染模擬研究
        養(yǎng)殖業(yè)整治需多措并舉
        農(nóng)業(yè)面源污染對水質(zhì)的影響及防治對策
        春季和夏季巢湖浮游生物群落組成及其動(dòng)態(tài)分析
        欧美一级在线全免费| 大陆国产乱人伦| 亚洲啪av永久无码精品放毛片| 这里有精品可以观看| 婷婷激情五月综合在线观看| 日本免费看一区二区三区| 国产激情一区二区三区| 国产肉丝袜在线观看| 香港三级欧美国产精品| 亚洲av狠狠爱一区二区三区| 熟妇人妻无乱码中文字幕真矢织江| 99精品国产综合久久久久五月天| av资源在线看免费观看| 国产免费人成视频在线观看播放播| 天天躁日日躁aaaaxxxx| 亚洲欧美日韩综合久久| 狠狠色欧美亚洲综合色黑a | 久久精品av在线视频| 亚洲av日韩av女同同性| 日日噜噜夜夜狠狠久久无码区| 精品日韩欧美| 精品精品国产三级av在线 | 亚洲av熟女中文字幕| 东北老女人高潮大喊舒服死了| 亚洲成人福利在线观看| 中文字幕人妻少妇久久| av影院在线免费观看不卡| 日韩制服国产精品一区| 久久精品国产亚洲综合色| 亚洲av成熟国产精品一区二区| 少妇裸体性生交| 国产女精品视频网站免费| 宅男久久精品国产亚洲av麻豆 | 亚洲国产高清在线一区二区三区| 国内精品视频成人一区二区| 国产午夜免费啪视频观看| 和黑人邻居中文字幕在线| 亚洲一区综合精品狠狠爱| 手机在线观看成年人视频| 国产极品粉嫩福利姬萌白酱| 夜爽8888视频在线观看|