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

    洋河水質(zhì)現(xiàn)狀及其成因分析

    2019-09-21 02:28:32張敏李令軍趙文慧郭逍宇許金浩王鑫龍趙文吉
    生態(tài)科學(xué) 2019年4期
    關(guān)鍵詞:洋河土地利用流域

    張敏, 李令軍, 趙文慧, 郭逍宇, 許金浩, 王鑫龍, 趙文吉,*

    洋河水質(zhì)現(xiàn)狀及其成因分析

    張敏1, 李令軍2, *, 趙文慧2, 郭逍宇1, 許金浩1, 王鑫龍1, 趙文吉1,*

    1. 首都師范大學(xué)資源環(huán)境與旅游學(xué)院, 北京 100048 2. 北京市環(huán)境保護(hù)監(jiān)測(cè)中心, 北京 100048

    為探究洋河水質(zhì)現(xiàn)狀, 利用2017年11月至12月水質(zhì)實(shí)測(cè)數(shù)據(jù), 分析了洋河10項(xiàng)水質(zhì)指標(biāo)的異質(zhì)性及其成因。結(jié)果表明: (1)單因子指數(shù)評(píng)價(jià)結(jié)果: 洋河總體水質(zhì)較差, 不同地段差異明顯。基于各斷面水質(zhì)現(xiàn)狀與評(píng)價(jià)結(jié)果來看, 總體達(dá)標(biāo)率為25%, 下游污染較嚴(yán)重達(dá)標(biāo)率為0%, 主要超標(biāo)因子為COD、TP、DO、NH4+-N和濁度。(2)分析上下游水質(zhì)變化特征, 發(fā)現(xiàn)上游水體中濁度變異系數(shù)最大, 其次是TP和BOD5; 下游水體中TP變異系數(shù)最大, 其次是BOD5, 表明水質(zhì)呈現(xiàn)一定變異性。(3)基于遙感影像, 提取河岸兩側(cè)一定范圍內(nèi)的土地利用類型并與水質(zhì)要素進(jìn)行相關(guān)性分析, 洋河上游主要以面污染源為主, 下游主要以點(diǎn)源污染為主。耕地、城鎮(zhèn)、建設(shè)用地、水面和草地對(duì)水質(zhì)的影響較顯著。(4)統(tǒng)計(jì)分析流域區(qū)縣社會(huì)經(jīng)濟(jì)結(jié)構(gòu)并建立與水質(zhì)要素間相關(guān)關(guān)系, 結(jié)果顯示洋河流域水生生態(tài)系統(tǒng)較脆弱, 下游受工業(yè)企業(yè)的污染貢獻(xiàn)較大。人口、農(nóng)產(chǎn)品、水產(chǎn)品產(chǎn)量和綜合能源消耗量對(duì)水質(zhì)影響顯著。

    水質(zhì)現(xiàn)狀; 土地利用; 能源產(chǎn)業(yè)結(jié)構(gòu); 相關(guān)性

    0 前言

    洋河作為官?gòu)d水庫(kù)的主要入庫(kù)河流之一, 多年來倍受關(guān)注。由于洋河及其他河流流域管理措施不當(dāng)、生態(tài)環(huán)境脆弱以及大量污水排入河道等原因?qū)е氯霂?kù)河流水質(zhì)嚴(yán)重超標(biāo)[1], 1997年水庫(kù)被迫退出城市生活飲用水體系[2]。之后經(jīng)過各個(gè)政府相關(guān)部門的努力工作, 水質(zhì)得到改善。但是, 受人類活動(dòng)的影響, 水庫(kù)水質(zhì)現(xiàn)狀依然不容樂觀, 洋河流域非點(diǎn)源污染現(xiàn)象嚴(yán)重, 改善水庫(kù)水質(zhì)加強(qiáng)流域綜合管理的工作依然艱巨[3]。已有許多學(xué)者證明土地利用類型與水環(huán)境質(zhì)量之間存在明顯的相關(guān)關(guān)系[4-5], 土地利用的方式、強(qiáng)度、空間布局和覆蓋類型所導(dǎo)致水環(huán)境能量傳輸發(fā)生變化而顯著影響水環(huán)境質(zhì)量[6]。此外, 還有研究表明流域社會(huì)經(jīng)濟(jì)發(fā)展與水資源間存在緊密的聯(lián)系[7], 區(qū)域產(chǎn)業(yè)結(jié)構(gòu)、能源結(jié)構(gòu)的不合理會(huì)對(duì)水資源造成一定危機(jī)[8]。因此, 為應(yīng)對(duì)北京市水資源匱乏和洋河流水質(zhì)下降的現(xiàn)狀, 研究洋河水質(zhì)現(xiàn)狀與流域土地利用類型的響應(yīng)關(guān)系, 并分析洋河流域能源產(chǎn)業(yè)結(jié)構(gòu)以及水資源現(xiàn)狀有著重要意義。

    洋河流域早期研究集中于水體水質(zhì)與河岸不同土地利用類型的土壤中重金屬間相應(yīng)關(guān)系的研究。2014年譚冰等學(xué)者研究洋河流域?qū)偃f全縣一帶受重金屬污染存在潛在風(fēng)險(xiǎn), 主要以Cr和Hg污染嚴(yán)重, 導(dǎo)致水質(zhì)已不滿足漁業(yè)養(yǎng)殖用水[9]。商靖敏等學(xué)者研究了洋河流域不同土地利用類型的表層土壤中Se含量存在顯著的空間異質(zhì)性[10]。龐博等學(xué)者研究了洋河流域張家口段河流水質(zhì)演化及其驅(qū)動(dòng)因子, 表明該河段主要污染物有氨氮、高錳酸鹽指數(shù)、BOD5,重金屬Cr6+和Cd等, 主要的驅(qū)動(dòng)因子是土地利用類型的變化[11]。此外, 另有學(xué)者致力于研究洋河流域生態(tài)服務(wù)的價(jià)值, 羅維等學(xué)者研究了從1990—2013年洋河流域土地利用的時(shí)空變異對(duì)生態(tài)服務(wù)功能價(jià)值的影響, 發(fā)現(xiàn)洋河流域土地利用除了林地和未利用地, 其他類型用地的年總生態(tài)服務(wù)價(jià)值先增加后減少, 總體呈現(xiàn)減少的趨勢(shì), 林地、草地和耕地是洋河流域生態(tài)服務(wù)的主要貢獻(xiàn)者和生態(tài)服務(wù)價(jià)值的主要敏感因子[12]。綜上所所述, 針對(duì)洋河流域水環(huán)境保護(hù)現(xiàn)有的研究集中于水體中水質(zhì)要素與近岸土地利用類型間的相關(guān)性分析, 本文將在此基礎(chǔ)上進(jìn)一步分析流域能源產(chǎn)業(yè)結(jié)構(gòu)對(duì)洋河的污染貢獻(xiàn)。

    根據(jù)調(diào)查結(jié)果顯示, 洋河上游與下游水環(huán)境質(zhì)量存在較大差異, 上下游流域河岸帶土地利用和流域區(qū)縣能源產(chǎn)業(yè)結(jié)構(gòu)出現(xiàn)不同比例的變化。為有針對(duì)性的改善洋河水環(huán)境質(zhì)量, 本項(xiàng)研究將分上游和下游兩大區(qū)域探究上下游水質(zhì)空間異質(zhì)性并分析其成因。運(yùn)用單因子評(píng)價(jià)法評(píng)價(jià)洋河水質(zhì)現(xiàn)狀, 分析上下游水質(zhì)存在的差別, 以及上下游水質(zhì)與河岸土地利用類型的相關(guān)關(guān)系, 探究水質(zhì)要素對(duì)土地利用的響應(yīng)關(guān)系, 結(jié)合年鑒統(tǒng)計(jì)數(shù)據(jù)分析流域能源產(chǎn)業(yè)結(jié)構(gòu)對(duì)河流造成的污染, 最終確定上游和下游存在的主要污染源。結(jié)果顯示, 洋河上游主要以面源污染為主, 下游主要以點(diǎn)源污染為主。

    1 材料與方法

    1.1 研究區(qū)概況

    洋河是永定河的支流, 位于113°29′32″E—115° 43′40″E, 39°59′42″N—41°15′3″N之間。地處我國(guó)北方農(nóng)牧交錯(cuò)帶, 流域近河岸帶1500 m范圍內(nèi)土地覆蓋類型主要以耕地為主, 流經(jīng)村鎮(zhèn)河岸帶以居民地、草地和林地為主, 偶爾會(huì)有工廠和建設(shè)用地。洋河流域主要包括三大支流: 南洋河、東、西洋河。西洋河發(fā)源于烏蘭察布的興和縣, 南洋河發(fā)源于山西省的陽(yáng)高縣, 在懷安縣附近匯合, 稱為洋河。沿途流經(jīng)張家口市、宣化、下花園、懷來等地后, 在懷來縣朱官屯與桑干河一并匯入永定河, 注入北京官?gòu)d水庫(kù)。流域面積15078 km2, 總長(zhǎng)度262 km, 平均海拔1000 m以上[13]。洋河位于東亞大陸季風(fēng)氣候區(qū), 夏季涼爽, 冬季寒冷漫長(zhǎng), 年降水量(375±25) mm, 年均氣溫7.7℃, 屬于典型的干旱、半干旱氣候, 河流主要以降水為補(bǔ)給來源, 河岸帶植被覆蓋類型單一, 生態(tài)環(huán)境脆弱[14]。

    1.2 數(shù)據(jù)獲取與處理

    1.2.1 采樣及樣品測(cè)定

    根據(jù)對(duì)洋河流域的實(shí)地考察與遙感影像結(jié)合分析, 在洋河干流以及匯入干流的支流布設(shè)監(jiān)測(cè)斷面。對(duì)于經(jīng)過村莊、城鎮(zhèn)和工廠、農(nóng)田等面污染源的河段加設(shè)采樣點(diǎn)。 2017年11月至12月在洋河干流布設(shè)8個(gè)監(jiān)斷面。如圖1所示為洋河監(jiān)測(cè)點(diǎn)布設(shè)情況及河岸緩沖區(qū)內(nèi)土地利用類型信息。

    由于條件限制, 現(xiàn)場(chǎng)監(jiān)測(cè)所用的水質(zhì)儀只能測(cè)定所需指標(biāo)中7項(xiàng), 其余的樣品采集后進(jìn)行實(shí)驗(yàn)室分析?,F(xiàn)場(chǎng)監(jiān)測(cè)與采樣時(shí), 每個(gè)監(jiān)測(cè)斷面等距布設(shè)3個(gè)采樣點(diǎn), 水質(zhì)儀測(cè)定三次后取平均值作為該斷面水質(zhì)參數(shù)值; 實(shí)驗(yàn)室分析則用采樣器采集液面下0.5 m處的水樣, 將3個(gè)采樣點(diǎn)的水樣均勻混合后保存至聚乙烯瓶中, 帶回實(shí)驗(yàn)室后測(cè)定基礎(chǔ)指標(biāo)。水樣的采集、保存、運(yùn)輸依據(jù)《水和廢水監(jiān)測(cè)分析方法》( 第4版)?,F(xiàn)場(chǎng)測(cè)定使用OTT Hydrolab DS5X多參數(shù)水質(zhì)監(jiān)測(cè)儀和手持GPS定位儀, 對(duì)采樣點(diǎn)的水體進(jìn)行實(shí)時(shí)在線測(cè)定, 測(cè)定包括溶解氧(DO)、pH、氨氮(NH4+-N)、總?cè)芙夤腆w(TDS)、藍(lán)藻素(SPC)、氧化還原電位(ORP)、濁度等水質(zhì)參數(shù), 并同時(shí)記錄采樣點(diǎn)的經(jīng)緯度與高程。所要測(cè)定的化學(xué)需氧量(COD)、五日生化需氧量(BOD5)、總磷(TP)等指標(biāo)嚴(yán)格按照《地表水和污水監(jiān)測(cè)技術(shù)規(guī)范》(HJ/T91-2002)中規(guī)定的方法采集樣品后帶回實(shí)驗(yàn)室測(cè)定。

    圖1 研究區(qū)水質(zhì)監(jiān)測(cè)采樣點(diǎn)及土地利用類型信息

    Figure 1 Water quality monitoring sampling point and land use type information in the study area

    1.2.2 遙感數(shù)據(jù)處理及土地利用信息提取

    選用研究區(qū)2017年2月和8月兩期高分影像數(shù)據(jù), 在幾何糾正和大氣糾正基礎(chǔ)上, 以主干河流兩側(cè)1500 m為緩沖區(qū), 采用面向?qū)ο蠓诸愑?jì)數(shù), 提取區(qū)內(nèi)土地利用信息。根據(jù)《土地利用現(xiàn)狀分類標(biāo)準(zhǔn)》(GB/T21010-2007)一、二級(jí)分類標(biāo)準(zhǔn)進(jìn)行分類, 選用包括耕地、園地、林地、城鎮(zhèn)、草地、裸地、水面(本文為水產(chǎn)養(yǎng)殖池)和建設(shè)用地(包括工廠)八類土地覆蓋類型為研究對(duì)象, 并獲取緩沖區(qū)內(nèi)各類土地利用的矢量數(shù)據(jù)。

    1.2.3 洋河流域各區(qū)縣經(jīng)濟(jì)統(tǒng)計(jì)數(shù)據(jù)

    選取2016年《張家口經(jīng)濟(jì)年鑒》中洋河流域各區(qū)縣生產(chǎn)總值、規(guī)模以上工業(yè)企業(yè)能源消費(fèi)情況等作為原始數(shù)據(jù)進(jìn)行分析。

    1.2.4 研究方法

    單因子指數(shù)評(píng)價(jià)法單因子評(píng)價(jià)法是指在參與水質(zhì)評(píng)價(jià)的監(jiān)測(cè)點(diǎn)中, 各項(xiàng)水質(zhì)指標(biāo)的實(shí)測(cè)值與目標(biāo)水質(zhì)指標(biāo)濃度限值的比值[15]。

    描述統(tǒng)計(jì)分析利用spss19.0軟件的描述統(tǒng)計(jì)功能可分析數(shù)據(jù)樣本的結(jié)構(gòu)和總體情況, 反映數(shù)據(jù)背后的實(shí)際情況。對(duì)水質(zhì)要素的數(shù)據(jù)進(jìn)行總體分析, 所得平均值、中位數(shù)、最大最小值和變異系數(shù)可直觀的表達(dá)洋河水質(zhì)的集中趨勢(shì)和現(xiàn)狀。

    Pearson相關(guān)性[16]利用spss19.0軟件的相關(guān)性分析功能分析洋河監(jiān)測(cè)點(diǎn)周邊土地利用面積與水質(zhì)要素的相關(guān)性, 可反映土地利用類型對(duì)水質(zhì)產(chǎn)生的影響[17]。

    3 結(jié)果與分析

    3.1 洋河流域水質(zhì)現(xiàn)狀特征分析

    采用《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB3838-2002)中III類標(biāo)準(zhǔn)為水質(zhì)目標(biāo)利用單因子污染指數(shù)法評(píng)價(jià)水質(zhì)現(xiàn)狀。見表1, 從洋河整體流域來看, 水質(zhì)較差, 達(dá)標(biāo)率僅為25%; 主要超標(biāo)因子為COD、TP、DO。pH值范圍在8.1—11.23, 呈堿性; 遠(yuǎn)超出地表水環(huán)境III類標(biāo)準(zhǔn)。

    分上下游來看, 洋河上游達(dá)標(biāo)率為50%, 而下游達(dá)標(biāo)率僅為0 ; 說明洋河下游污染狀況比上游要嚴(yán)重。洋河上游IV類占河段50%; 洋河下游V類水質(zhì)占河段的50%, 洋河下游COD在各斷面均超標(biāo)。在87號(hào)斷面水質(zhì)指標(biāo)TP超編現(xiàn)象較為嚴(yán)重, 導(dǎo)致該河段歸為劣V類水質(zhì), 現(xiàn)場(chǎng)觀察到水體呈現(xiàn)淺褐色。結(jié)果顯示, 洋河流域自上游向下游各水質(zhì)要素呈現(xiàn)增加的趨勢(shì), COD值顯著增加, 表現(xiàn)出下游河段水質(zhì)污染較上游嚴(yán)重。為探究這一現(xiàn)象, 將洋河流域分上下游兩個(gè)單元進(jìn)行水質(zhì)污染特征分析。

    將洋河流域的各指標(biāo)均值與對(duì)應(yīng)的地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)對(duì)比后可知整個(gè)流域DO、COD、TP均超過了地表水Ⅲ類標(biāo)準(zhǔn), pH值呈堿性。是洋河流域應(yīng)重點(diǎn)加以控制的主要因子。從上下游河段統(tǒng)計(jì)特征情況看: (1)洋河上游表現(xiàn)為DO、COD、TP和濁度均超過了地表水Ⅲ類標(biāo)準(zhǔn), 除濁度超標(biāo)率達(dá)100%, 其余水質(zhì)因子超標(biāo)率均達(dá)25%。pH均值為9.895, 最小值為8.1, 上游呈弱堿性并超出正常值范圍, 是上游應(yīng)重點(diǎn)加以控制的主要因子。從變異系數(shù)來看, 表現(xiàn)為濁度變異系數(shù)最大(160%), 其次TP是(85.8%)和BOD5是(70%); 離散程度較高。(2)洋河下游表現(xiàn)為DO、COD、TP、NH4+-N和濁度均超出地表水Ⅲ類標(biāo)準(zhǔn)。其中TP、濁度和pH值超標(biāo)率達(dá)100%, 其余水質(zhì)因子超標(biāo)率為25%。pH均值為11.21, 最小值為10.83, 下游呈現(xiàn)強(qiáng)堿性嚴(yán)重超出正常范圍, 是洋河下游應(yīng)重點(diǎn)加以控制的水質(zhì)因子。從變異系數(shù)來看, 表現(xiàn)為TP變異系數(shù)最大(130.1%), 其次是BOD5(98.1%); 離散程度較高。

    3.2 水質(zhì)異質(zhì)性成因分析

    3.2.1 河岸土地利用類型與水質(zhì)變化的關(guān)系

    對(duì)近河岸帶500 m、1000 m、1500 m區(qū)內(nèi)土地利用類型信息樣本統(tǒng)計(jì)分析。由圖2可知, 洋河河岸帶主要以耕地為主, 其次是城鎮(zhèn), 建設(shè)用地及工廠占總統(tǒng)計(jì)面積最少; 上游緩沖區(qū)內(nèi)三個(gè)尺度下耕地面積均占分類總面積的88%以上, 林地、城鎮(zhèn)用地都占3%—7%, 而建設(shè)用地及工廠占總面積不足1%; 下游緩沖區(qū)內(nèi)三個(gè)尺度下耕地均占分類總面積76%—77%, 城鎮(zhèn)用地占4%—10%, 草地占2%—9%, 建設(shè)用地及工廠用地占4%—5%, 水面(水產(chǎn)養(yǎng)殖池)占2%。說明洋河上游流域主要以面源污染為主, 河岸帶土地利用類型較為單一。洋河下游河岸耕地面積減少, 草地、園地、林地和建設(shè)用地增加導(dǎo)致耕地面積所占比例顯著減少。土地利用豐富度增加, 說明受人類活動(dòng)影響增加。城鎮(zhèn)化率提高、人口增加是導(dǎo)致土地利用類型發(fā)生空間變異的主導(dǎo)因素[18]。

    注: “- ”為未超標(biāo)

    圖2 洋河河岸帶緩沖區(qū)內(nèi)不同尺度下土地覆蓋類型信息

    Figure 2 Information on land cover types at different scales in the buffer zone of the Yanghe riparian zone

    為探究水質(zhì)要素變異性與土地覆蓋類型間的響應(yīng)關(guān)系, 選取代表性的5項(xiàng)水質(zhì)指標(biāo)作為數(shù)據(jù)樣本。如圖3 所示, 發(fā)現(xiàn)5項(xiàng)指標(biāo)均呈現(xiàn)出向下游污染加重的現(xiàn)象, 并且最大值均出現(xiàn)在下游。從各項(xiàng)指標(biāo)變化情況來看, COD變化最顯著, 下游水體中該項(xiàng)指標(biāo)明顯高于上游, 并且最大值出現(xiàn)在下游; 其次是BOD5, 上游水體中該項(xiàng)指標(biāo)值變化緩慢且接近于1(mg/L),下游水體中該指標(biāo)值變化顯著且出現(xiàn)最大值。TP、DO和NH4+-N變化較緩, 其中由水質(zhì)現(xiàn)狀評(píng)價(jià)結(jié)果來看, 下游水體中TP的濃度值要高于上游水體中的。從水質(zhì)要素累積量來看, 下游總累計(jì)量明顯大于上游, 其中87號(hào)點(diǎn)累計(jì)量達(dá)最高值其次是86號(hào)點(diǎn)和85號(hào)點(diǎn), 表明下游水體中污染物的量大于上游水體中的。結(jié)合圖2結(jié)果, 發(fā)現(xiàn)上游河岸兩側(cè)以農(nóng)田耕地為主要土地覆蓋類型, 偶有工廠出現(xiàn)。河岸兩側(cè)以河灘地為主, 面源所產(chǎn)生的污染物可直接進(jìn)入水體污染水質(zhì), 由此可認(rèn)為上游主要受到面源污染為主。由于下游河岸兩側(cè)的工廠、城鎮(zhèn)等用地增多, 耕地面積所占比例減少。由于建設(shè)用地及工廠、城鎮(zhèn)等用地的增加, 下游水體受到點(diǎn)源污染加重??傮w來看, 洋河上游主要以面污染源為主; 洋河下游則傾向于點(diǎn)污染源。

    3.2.2 水質(zhì)要素與河岸帶土地利用類型間相關(guān)性分析

    以洋河兩岸不同緩沖區(qū)內(nèi)各土地利用類型面積和9項(xiàng)水質(zhì)指標(biāo)作為基礎(chǔ)分析數(shù)據(jù), 計(jì)算土地利用類型與水質(zhì)之間的Pearson相關(guān)性, 結(jié)果見表2。

    洋河河岸帶500 m緩沖區(qū)內(nèi), COD、ORP、DO和濁度與土地利用類型的面積未出現(xiàn)顯著相關(guān)性, BOD5與耕地呈現(xiàn)高度相關(guān)性; TP和NH4+-N與建設(shè)用地、裸地和城鎮(zhèn)的面積呈高度相關(guān)性, NH4+-N還與草地面積呈較高的相關(guān)性; TDS和SPC與水面(水產(chǎn)養(yǎng)殖池)面積呈現(xiàn)很高的相關(guān)性。

    洋河河岸帶1000 m緩沖區(qū)內(nèi), DO和濁度與土地利用類型的面積未出現(xiàn)顯著相關(guān)性; COD與草地呈現(xiàn)很高的相關(guān)性; BOD5與水面、城鎮(zhèn)和耕地呈現(xiàn)很高的相關(guān)性; TP和NH4+-N與建設(shè)用地、裸地、園地的面積呈高度相關(guān)性, 此外TP與城鎮(zhèn), NH4+-N與草地也呈現(xiàn)較高的相關(guān)性; TDS和ORP與林地和水面面積呈現(xiàn)較高的相關(guān)性; ORP與水面面積呈現(xiàn)較高的相關(guān)性。

    圖3 采樣點(diǎn)水質(zhì)要素濃度累積柱狀圖

    Figure 3 Concentration accumulation histogram of water Quality Element at Sampling Point

    洋河河岸帶1500 m緩沖區(qū)內(nèi), DO和NH4+-N與各土地利用類型未出現(xiàn)顯著相關(guān)性, 其余水質(zhì)要素均與土地利用類型面積均呈現(xiàn)一定正相關(guān)。其中, BOD5與水面、城鎮(zhèn)和耕地顯著相關(guān), 與耕地相關(guān)性最顯著; COD與建設(shè)用地及工廠、草地呈現(xiàn)出顯著的正相關(guān)關(guān)系; TP與建設(shè)用地及工廠、裸地和城鎮(zhèn)呈現(xiàn)出顯著正相關(guān)關(guān)系。TDS和SPC均與林地呈現(xiàn)出顯著正相關(guān); ORP與耕地呈現(xiàn)出現(xiàn)在正相關(guān)關(guān)系; 濁度只與城鎮(zhèn)表現(xiàn)正相關(guān)關(guān)系。說明河岸帶緩沖區(qū)內(nèi)土地利用類型對(duì)水質(zhì)有不同程度影響, 流域土地利用類型的變化可以引起水質(zhì)的變化。總體來看, 河岸1000m緩沖區(qū)內(nèi)土地利用面積與水質(zhì)呈現(xiàn)的相關(guān)性程度最高, 說明該范圍內(nèi)污染源對(duì)河流水質(zhì)產(chǎn)生最為復(fù)雜的影響。應(yīng)對(duì)該范圍內(nèi)景觀格局制定合理的生態(tài)治理與生態(tài)保護(hù)政策。

    3.2.3 流域區(qū)縣能源和產(chǎn)業(yè)結(jié)構(gòu)污染貢獻(xiàn)分析

    由2016年張家口市經(jīng)濟(jì)年鑒中統(tǒng)計(jì)數(shù)據(jù)進(jìn)行分析, 結(jié)果如圖4所示。可以看出洋河上、下游流域產(chǎn)業(yè)結(jié)構(gòu)的變化對(duì)比: 由上游至下游, 第二產(chǎn)業(yè)所占比重增加且增量最高, 其次是第一產(chǎn)業(yè)稍有增加; 第三產(chǎn)業(yè)所占比重顯著降低。上游流域區(qū)縣第一、二、三產(chǎn)業(yè)比重為14.92:34.28:50.79, 到下游變?yōu)?6.48:39.14:44.38其中第二產(chǎn)業(yè)增加了近5個(gè)百分點(diǎn), 表明工業(yè)成為經(jīng)濟(jì)增長(zhǎng)的主體。隨著經(jīng)濟(jì)不斷增長(zhǎng), 產(chǎn)業(yè)結(jié)構(gòu)不能得到合理調(diào)配, 再加上傳統(tǒng)的污水處理設(shè)備及技術(shù)不能滿足工業(yè)規(guī)模的擴(kuò)大速度, 區(qū)域污水處理效率低下, 可導(dǎo)致水環(huán)境質(zhì)量下降[19]。

    表2 不同緩沖區(qū)內(nèi)土地利用面積與水質(zhì)要素含量間Pearson相關(guān)性

    注: “**”為在 .01 水平(雙側(cè))上顯著相關(guān)?!?”為在 0.05 水平(雙側(cè))上顯著相關(guān)。“-”為負(fù)相關(guān)

    圖4 洋河流域各區(qū)縣產(chǎn)業(yè)結(jié)構(gòu)信息

    Figure 4 Industrial structure information of urban areas and county towns in the Yanghe River Basin

    能源消耗與用水間存在復(fù)雜關(guān)系, 能源的開采與使用過程中需消耗大量水資源, 并產(chǎn)生大量工業(yè)廢水。同時(shí), 在取水、用水和回收處理過程中, 需要消耗大量能源。面對(duì)區(qū)域水環(huán)境質(zhì)量下降的問題, 基于水資源與能源間的相互關(guān)系, 本文將以洋河流域區(qū)縣工業(yè)企業(yè)能源消費(fèi)量為原始數(shù)據(jù)進(jìn)行分析, 結(jié)果如圖5 所示。

    發(fā)現(xiàn)洋河上游與下游流域區(qū)縣能源結(jié)構(gòu)變化顯著。其中焦炭增加量最大, 其次是能源消費(fèi)總量、原煤量和電力消費(fèi)量, 最后是汽油和柴油用量。煤炭與煤制油、氣兩大企業(yè)是用水較大的企業(yè), 企業(yè)運(yùn)營(yíng)過程中產(chǎn)生大量廢水, 對(duì)水環(huán)境污染起到較大貢獻(xiàn)[20]。綜上所述, 可以認(rèn)為下游污染主要由于工業(yè)企業(yè)的活動(dòng)影響, 流域區(qū)縣能源、產(chǎn)業(yè)結(jié)構(gòu)的合理調(diào)控是水環(huán)境保護(hù)的重點(diǎn)工作。

    3.2.4 流域區(qū)縣經(jīng)濟(jì)因素與水質(zhì)間相關(guān)關(guān)系

    統(tǒng)計(jì)洋河流域各區(qū)縣經(jīng)濟(jì)年鑒中8項(xiàng)數(shù)據(jù)作為影響水質(zhì)的環(huán)境因子, 計(jì)算這些環(huán)境因子與水質(zhì)要素之間的相關(guān)關(guān)系, 結(jié)果如表3所示。

    圖5 洋河流域各區(qū)縣能源消費(fèi)現(xiàn)狀

    Figure 5 Current status of energy consumption in urban areas and counties of Yanghe River Basin

    表3 洋河流域社會(huì)經(jīng)濟(jì)與水質(zhì)要素間Pearson相關(guān)性

    注: “**”為在 .01 水平(雙側(cè))上顯著相關(guān)?!?”為在 0.05 水平(雙側(cè))上顯著相關(guān)。“-”為負(fù)相關(guān)

    流域各區(qū)縣的經(jīng)濟(jì)因素作為影響水質(zhì)要素的因子來看: BOD5受各類因子影響最復(fù)雜, 其中人口與該要素呈顯著正相關(guān)關(guān)系, 而且主要是與農(nóng)村人口呈現(xiàn)顯著正相關(guān)關(guān)系; 農(nóng)產(chǎn)品、水產(chǎn)品的產(chǎn)量以及區(qū)域GDP與該水質(zhì)要素呈現(xiàn)較顯著的正相關(guān)關(guān)系, 由此可以表明BOD5受農(nóng)村居民的生活活動(dòng)影響最大。NH4+-N和COD受綜合能源消耗量影響較大, 相互之間呈現(xiàn)出顯著正相關(guān)關(guān)系, 說明NH4+-N和COD受到工業(yè)活動(dòng)的影響最大。結(jié)合水質(zhì)評(píng)價(jià)結(jié)果來看, COD是洋河主要超標(biāo)的水質(zhì)要素之一, 保護(hù)水環(huán)境質(zhì)量重點(diǎn)應(yīng)該落實(shí)到對(duì)區(qū)域能源結(jié)構(gòu)調(diào)整的工作上。

    4 結(jié)論

    綜上所述, 洋河流域水質(zhì)污染較嚴(yán)重, 下游比上游污染更嚴(yán)重。上游主要以DO、COD污染物為主, 表明上游水體中有機(jī)污染物含量較高, 經(jīng)土地利用類型的分析結(jié)果顯示, 洋河流域整體以耕地面積為主, 尤其集中于上游, 多以面源污染為主, 是水環(huán)境保護(hù)的重點(diǎn)。下游主要以COD、TP污染物為主, 建設(shè)用地和城鎮(zhèn)面積增加, 主要表現(xiàn)為點(diǎn)污染源。通過水質(zhì)要素與土地利用類型間相關(guān)性分析, 表明河岸帶緩沖區(qū)內(nèi)各水質(zhì)要素與土地利用類型面積呈現(xiàn)出不同程度的正相關(guān)關(guān)系。導(dǎo)致水體污染的因素主要為耕地、城鎮(zhèn)、建設(shè)用地、水面和草地。對(duì)洋河流域各區(qū)縣能源、產(chǎn)業(yè)結(jié)構(gòu)的分析。發(fā)現(xiàn)洋河流域自上游向下游, 工業(yè)越來越成為流域經(jīng)濟(jì)增長(zhǎng)的主體, 能源消費(fèi)水平顯著增加。表明, 洋河下游流域受工業(yè)企業(yè)的污染較大。為保護(hù)水生態(tài)的穩(wěn)定需要大力調(diào)節(jié)當(dāng)前社會(huì)經(jīng)濟(jì)的發(fā)展方向, 將資源能源型社會(huì)轉(zhuǎn)向環(huán)境友好型。

    [1] 于淼, 魏源送, 劉俊國(guó), 等. 永定河(北京段)水資源、水環(huán)境的變遷及流域社會(huì)經(jīng)濟(jì)發(fā)展對(duì)其影響[J]. 環(huán)境科學(xué)學(xué)報(bào), 2011, 31(09): 1817–1825.

    [2] 孫寓姣, 陳程, 丁愛中, 等. 官?gòu)d水庫(kù)水質(zhì)特征及水體微生物多樣性的響應(yīng)[J]. 中國(guó)環(huán)境科學(xué), 2015, 35(05): 1547–1553.

    [3] YI Huiwen, GUI Huanliu, HAI Guanghao, et al. Analysis on the Temporal and Spatial Pattern of Ecosystem Services Based on Land-Use Changes: A Case Study of Guanting Reservoir Watershed in China[J]. Advanced Materials Research, 2013, 2695(807).

    [4] 王娟, 張飛, 張?jiān)? 等. 艾比湖區(qū)域水質(zhì)空間分布特征及其與土地利用/覆被類型的關(guān)系[J]. 生態(tài)學(xué)報(bào), 2016, 36(24): 7971–7980.

    [5] 項(xiàng)頌, 龐燕, 竇嘉順, 等. 不同時(shí)空尺度下土地利用對(duì)洱海入湖河流水質(zhì)的影響[J/OL]. 生態(tài)學(xué)報(bào), 2018(03): 1–10[2018-01-29].

    [6] 于術(shù)桐, 黃賢金, 程緒水. 土地利用變化的水環(huán)境效應(yīng)研究進(jìn)展[J]. 長(zhǎng)江流域資源與環(huán)境, 2010, 19(S2): 177– 181.

    [7] 孫艷芝, 魯春霞, 謝高地, 等. 北京城市發(fā)展與水資源利用關(guān)系分析[J]. 資源科學(xué), 2015, 37(06): 1124–1132.

    [8] 鐘銳, 趙鵬, 張宏偉, 等. 國(guó)際能源貿(mào)易中的能源隱含水流動(dòng)及其對(duì)各國(guó)水資源影響[J]. 資源科學(xué), 2014, 36(01): 103–109.

    [9] 譚冰, 王鐵宇, 朱朝云, 等. 洋河流域萬全段重金屬污染風(fēng)險(xiǎn)及控制對(duì)策[J]. 環(huán)境科學(xué), 2014, 35(02): 719–726.

    [10] 商靖敏, 羅維, 吳光紅, 等. 洋河流域不同土地利用類型土壤硒(Se)分布及影響因素[J]. 環(huán)境科學(xué), 2015, 36(01): 301–308.

    [11] 龐博, 王鐵宇, 呂永龍, 等. 洋河流域張家口段河流水質(zhì)演化及驅(qū)動(dòng)因子分析[J]. 環(huán)境科學(xué), 2013, 34(01): 379–384.

    [12] 羅維, 易海杰, 李紅舉, 等. 洋河流域土地利用時(shí)空變異及其對(duì)生態(tài)服務(wù)功能價(jià)值的影響[J]. 生態(tài)學(xué)報(bào), 2017, 37(16): 5342–5351.

    [13] 高佳佳, 羅維, 奚曉霞. 首都水源地——洋河流域人為源多環(huán)芳烴(PAHs)排放清單估算及其影響分析[J]. 環(huán)境科學(xué), 2014, 35(12): 4573–4581.

    [14] 徐蘭, 羅維, 周寶同. 基于土地利用變化的農(nóng)牧交錯(cuò)帶典型流域生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)——以洋河為例[J]. 自然資源學(xué)報(bào), 2015, 30(04): 580–590.

    [15] 曾永, 樊引琴, 王麗偉, 等. 水質(zhì)模糊綜合評(píng)價(jià)法與單因子指數(shù)評(píng)價(jià)法比較[J]. 人民黃河, 2007(02): 45+65.

    [16] 高永恒, 曾曉陽(yáng), 周國(guó)英, 王根緒. 長(zhǎng)江源區(qū)高寒濕地植物群落主要種群種間關(guān)系分析[J]. 濕地科學(xué), 2011, 9(01): 1–7.

    [17] XIA Jun, WANG Longfeng, YU Jingjie, et al. Impact of environmental factors on water quality at multiple spatial scales and its spatial variation in Huai River Basin, China[J]. Science China(Earth Sciences), 2018, 61(01): 82–92.

    [18] 韓會(huì)然, 楊成鳳, 宋金平. 北京市土地利用變化特征及驅(qū)動(dòng)機(jī)制[J]. 經(jīng)濟(jì)地理, 2015, 35(05): 148–154+197.

    [19] 李志濤, 黃河清, 張明慶, 等. 鄱陽(yáng)湖流域經(jīng)濟(jì)增長(zhǎng)與水環(huán)境污染關(guān)系研究[J]. 資源科學(xué), 2010, 32(02): 267– 273.

    [20] 顧阿倫, 姜冬梅, 張?jiān)? 能源–水關(guān)系研究現(xiàn)狀及對(duì)我國(guó)的啟示[J]. 生態(tài)經(jīng)濟(jì), 2016, 32(07): 20–23+28.

    Present situation of Yanghe water quality and analysis of its causes

    ZHANG Min1,LI Lingjun2,*, ZHAO Wenhui2, GUO Xiaoyu1, XU Jinhao1, WANG Xinlong1, ZHAO Wenji1,*

    1. College of Resource Environment and Tourism, Capital Normal University, Beijing 100048, China 2. Beijing Municipal Environment Monitoring Center, Beijing 100048, China

    The water quality of Yanghe River between November and December 2017 was monitored, with spatial heterogeneities and their dependencies of ten indices analyzed. The results show that: (1) Based on the results of single-factor pollution index, the overall water quality of Yanghe River was poor, and presented distinct spatial differences. The water quality evaluation results of only 25% of the sampling sites along the river met the standard. Water quality in the downstream became even worse, with none of the water samples collected meeting the standard. The main indices which exceeded the standard included COD, TP, DO, NH4+-N and turbidity. (2) Spatial variation of the water quality was observed. In the upstream, the coefficient of variation of turbidity was the greatest, followed by TP and BOD5; while in the downstream, the coefficient of variation of TP was the greatest, followed by BOD5. (3) The correlation between water quality and land use types within a certain range of river banks on both sides extracted by remote sensing data was conducted. In the upstream, non-pointsources were the dominant pollution sources; while in the downstream, point sources were dominant. The water quality was significantly affected by cropland, town and construction land, water surface and grassland. (4) The correlation between water quality and socio-economic factors of counties in the river basin was analyzed. The aquatic ecosystems in the Yanghe River Basin were relatively vulnerable, and the downstream of river was heavily polluted by the industry. Socio-economic factors such as population, agricultural products, aquatic production and integrated energy consumption had significant impacts on water quality.

    water quality status; land use; energy industry structure; correlation

    10.14108/j.cnki.1008-8873.2019.04.012

    X522

    A

    1008-8873(2019)04-077-08

    2018-09-02;

    2018-10-23

    國(guó)家青年科學(xué)基金項(xiàng)目(41101404); 水源保護(hù)區(qū)水質(zhì)情況評(píng)價(jià)方法研究(01117220010037)

    張敏(1994—), 女, 內(nèi)蒙古鄂爾多斯市人, 碩士生,主要從水土污染研究, E-mail:455944412@qq.com

    趙文吉E-mail: zhwenji1215@163.com; 李令軍lilj2000@126.com

    張敏, 李令軍, 趙文慧, 等. 洋河水質(zhì)現(xiàn)狀及其成因分析[J]. 生態(tài)科學(xué), 2019, 38(4): 77-84.

    ZHANG Min, LI Lingjun, ZHAO Wenhui, et al. Present situation of Yanghe water quality and analysis of its causes[J]. Ecological Science, 2019, 38(4): 77-84.

    猜你喜歡
    洋河土地利用流域
    Generation of hyperentangled photon pairs based on lithium niobate waveguide
    壓油溝小流域
    天鵝洋河舞翩躚
    考試與招生(2022年3期)2022-03-15 07:24:24
    千年酒鎮(zhèn) 醉美洋河
    堡子溝流域綜合治理
    羅堰小流域
    水利規(guī)劃與設(shè)計(jì)(2018年1期)2018-01-31 01:53:37
    土地利用生態(tài)系統(tǒng)服務(wù)研究進(jìn)展及啟示
    濱??h土地利用挖潛方向在哪里
    秦皇島市洋河水電站增效擴(kuò)容改造工程經(jīng)濟(jì)評(píng)價(jià)論述
    日日撸夜夜添| 真实男女啪啪啪动态图| 国产黄a三级三级三级人| 91午夜精品亚洲一区二区三区| 人体艺术视频欧美日本| 看黄色毛片网站| 国产探花在线观看一区二区| av黄色大香蕉| 精品99又大又爽又粗少妇毛片| 99久国产av精品| 麻豆乱淫一区二区| 麻豆乱淫一区二区| 干丝袜人妻中文字幕| 欧美极品一区二区三区四区| 伦理电影大哥的女人| 2022亚洲国产成人精品| 欧美性感艳星| 麻豆av噜噜一区二区三区| 小蜜桃在线观看免费完整版高清| 2022亚洲国产成人精品| 校园人妻丝袜中文字幕| 欧美另类亚洲清纯唯美| 亚洲色图av天堂| 久久国产乱子免费精品| 色播亚洲综合网| 亚洲成a人片在线一区二区| 免费一级毛片在线播放高清视频| 男人舔奶头视频| 亚洲精品粉嫩美女一区| 麻豆乱淫一区二区| 免费看光身美女| 最好的美女福利视频网| 亚洲av一区综合| 欧美色欧美亚洲另类二区| 中文字幕久久专区| 大型黄色视频在线免费观看| 一级二级三级毛片免费看| 三级毛片av免费| 可以在线观看毛片的网站| 永久网站在线| 国产一区二区亚洲精品在线观看| 99热这里只有是精品在线观看| 91精品一卡2卡3卡4卡| 一进一出抽搐gif免费好疼| 国产亚洲91精品色在线| 国产中年淑女户外野战色| 亚洲av中文字字幕乱码综合| 国产精品久久久久久亚洲av鲁大| 在线观看美女被高潮喷水网站| 免费在线观看成人毛片| 日韩,欧美,国产一区二区三区 | 哪个播放器可以免费观看大片| 乱人视频在线观看| 久久久国产成人精品二区| 中国美白少妇内射xxxbb| 日产精品乱码卡一卡2卡三| 天堂√8在线中文| 99久国产av精品国产电影| av在线老鸭窝| 22中文网久久字幕| 久久综合国产亚洲精品| 特大巨黑吊av在线直播| av福利片在线观看| 中文精品一卡2卡3卡4更新| 亚洲av免费高清在线观看| 丝袜喷水一区| 身体一侧抽搐| 亚洲国产欧洲综合997久久,| 久久精品夜色国产| 国产精品野战在线观看| 欧美日韩在线观看h| 日韩成人av中文字幕在线观看| 性欧美人与动物交配| 91狼人影院| 国产三级在线视频| 青青草视频在线视频观看| 女人被狂操c到高潮| 成人永久免费在线观看视频| 搡女人真爽免费视频火全软件| 国产成人a∨麻豆精品| 欧美高清性xxxxhd video| 国产精品.久久久| 男人舔女人下体高潮全视频| 国产亚洲91精品色在线| 给我免费播放毛片高清在线观看| 日本一本二区三区精品| 不卡视频在线观看欧美| 美女国产视频在线观看| 国产精品一区www在线观看| 国产色婷婷99| 久久热精品热| 69av精品久久久久久| 蜜桃久久精品国产亚洲av| 国语自产精品视频在线第100页| 亚洲欧美成人综合另类久久久 | 久久久久久久久久久丰满| 三级男女做爰猛烈吃奶摸视频| 久久亚洲精品不卡| 欧美性感艳星| 麻豆av噜噜一区二区三区| 国内久久婷婷六月综合欲色啪| 国产成人福利小说| 中文精品一卡2卡3卡4更新| 亚洲av熟女| 久久精品国产亚洲网站| 国产在线精品亚洲第一网站| 人妻制服诱惑在线中文字幕| 亚洲成人精品中文字幕电影| 国产黄a三级三级三级人| 亚洲精品国产av成人精品| 亚洲久久久久久中文字幕| 日韩在线高清观看一区二区三区| 夜夜看夜夜爽夜夜摸| 欧美bdsm另类| 国产一区二区三区在线臀色熟女| 日本一二三区视频观看| 午夜福利成人在线免费观看| 中文字幕av在线有码专区| 亚洲第一区二区三区不卡| 亚洲欧美成人综合另类久久久 | 精品久久久久久久久亚洲| 国产成年人精品一区二区| 中文字幕熟女人妻在线| 插阴视频在线观看视频| 春色校园在线视频观看| av女优亚洲男人天堂| 午夜精品国产一区二区电影 | 免费人成在线观看视频色| 91精品一卡2卡3卡4卡| 国产精品一区二区在线观看99 | 一区二区三区免费毛片| 成人鲁丝片一二三区免费| 国产日韩欧美在线精品| 成人永久免费在线观看视频| 成人av在线播放网站| 精品免费久久久久久久清纯| 我要搜黄色片| 成人特级av手机在线观看| 国产一区亚洲一区在线观看| 在线观看一区二区三区| 天堂中文最新版在线下载 | 国产午夜精品一二区理论片| 婷婷亚洲欧美| 三级国产精品欧美在线观看| 国产极品天堂在线| 国产 一区精品| 哪里可以看免费的av片| 中国美女看黄片| 国产单亲对白刺激| 国产午夜精品一二区理论片| 久久婷婷人人爽人人干人人爱| 国产亚洲精品久久久com| 国产激情偷乱视频一区二区| 国产一区二区在线观看日韩| 乱系列少妇在线播放| 亚洲无线观看免费| 亚洲欧美成人综合另类久久久 | 国产精品电影一区二区三区| 国产不卡一卡二| 久久九九热精品免费| 人人妻人人看人人澡| 日本色播在线视频| 99在线视频只有这里精品首页| 在线国产一区二区在线| 国产成人91sexporn| 亚洲欧美中文字幕日韩二区| 大香蕉久久网| 综合色av麻豆| 尤物成人国产欧美一区二区三区| 日产精品乱码卡一卡2卡三| 免费人成视频x8x8入口观看| 欧美变态另类bdsm刘玥| 热99在线观看视频| 久久久久网色| 欧美激情久久久久久爽电影| 色哟哟·www| 男人狂女人下面高潮的视频| 男的添女的下面高潮视频| 国产精品女同一区二区软件| 成人美女网站在线观看视频| 黄片无遮挡物在线观看| 久久久久网色| 夫妻性生交免费视频一级片| 欧美性猛交╳xxx乱大交人| 好男人视频免费观看在线| av在线老鸭窝| 男人舔女人下体高潮全视频| 在线免费十八禁| 精品少妇黑人巨大在线播放 | 午夜精品在线福利| 亚洲欧美日韩东京热| 久久久久久久久久黄片| 欧美另类亚洲清纯唯美| 国产极品天堂在线| 国产精品嫩草影院av在线观看| 国产午夜精品论理片| 久久精品国产清高在天天线| 热99re8久久精品国产| 黄色一级大片看看| 亚洲欧洲日产国产| 观看免费一级毛片| www日本黄色视频网| 精品欧美国产一区二区三| 国产免费男女视频| 黑人高潮一二区| 欧美色欧美亚洲另类二区| 精品不卡国产一区二区三区| 免费看美女性在线毛片视频| 亚洲av成人精品一区久久| 中文欧美无线码| 国产亚洲av嫩草精品影院| 欧美日韩国产亚洲二区| 日日撸夜夜添| 韩国av在线不卡| 一个人看视频在线观看www免费| 婷婷精品国产亚洲av| 又爽又黄a免费视频| 国内精品久久久久精免费| 小说图片视频综合网站| 少妇被粗大猛烈的视频| 我要看日韩黄色一级片| 国产真实乱freesex| 亚洲丝袜综合中文字幕| 日本成人三级电影网站| 在线观看一区二区三区| 丝袜美腿在线中文| 91精品国产九色| 欧美日本视频| 日韩欧美 国产精品| 亚州av有码| 亚洲精品久久国产高清桃花| 色尼玛亚洲综合影院| 黄色配什么色好看| av国产免费在线观看| 国产精品女同一区二区软件| 青春草国产在线视频 | 网址你懂的国产日韩在线| 久久这里只有精品中国| 韩国av在线不卡| 日本黄色视频三级网站网址| 国产黄色小视频在线观看| 大又大粗又爽又黄少妇毛片口| 亚洲成人中文字幕在线播放| 在线免费十八禁| 免费av不卡在线播放| 嫩草影院精品99| 最近最新中文字幕大全电影3| 只有这里有精品99| 成人性生交大片免费视频hd| 国产熟女欧美一区二区| 亚洲高清免费不卡视频| www.av在线官网国产| 国产精品三级大全| 久久久久久久久久久丰满| 免费看a级黄色片| 久久国内精品自在自线图片| 小说图片视频综合网站| 夜夜爽天天搞| 九草在线视频观看| 我要搜黄色片| 久久亚洲国产成人精品v| 丰满人妻一区二区三区视频av| 成人二区视频| 国产三级在线视频| 精品无人区乱码1区二区| 亚洲第一区二区三区不卡| 丝袜美腿在线中文| 亚洲最大成人中文| 午夜亚洲福利在线播放| 中文资源天堂在线| 亚洲精品国产成人久久av| 老师上课跳d突然被开到最大视频| 亚洲欧美精品专区久久| 免费在线观看成人毛片| 亚洲精品自拍成人| 久久国内精品自在自线图片| 国产熟女欧美一区二区| 成人无遮挡网站| 亚洲最大成人中文| 欧美日韩乱码在线| 99热这里只有是精品50| 永久网站在线| 国产精品99久久久久久久久| 亚洲精品自拍成人| 97超视频在线观看视频| 亚洲国产精品成人综合色| 特大巨黑吊av在线直播| 在线天堂最新版资源| 日韩精品青青久久久久久| 中文字幕熟女人妻在线| 人妻少妇偷人精品九色| 联通29元200g的流量卡| 我的老师免费观看完整版| 亚洲一级一片aⅴ在线观看| 99久久精品国产国产毛片| 亚洲在线自拍视频| 亚洲一区二区三区色噜噜| 国产成人精品婷婷| 久久久色成人| 97超碰精品成人国产| 日本五十路高清| 我要搜黄色片| 久久精品国产亚洲av天美| 寂寞人妻少妇视频99o| 久久精品综合一区二区三区| 欧美另类亚洲清纯唯美| ponron亚洲| 高清毛片免费观看视频网站| 在线a可以看的网站| 99国产精品一区二区蜜桃av| 国产毛片a区久久久久| 青青草视频在线视频观看| 波多野结衣巨乳人妻| 又粗又硬又长又爽又黄的视频 | 日本黄色视频三级网站网址| 男插女下体视频免费在线播放| 国产精品一区二区三区四区免费观看| а√天堂www在线а√下载| 国产精品人妻久久久久久| 国产成人精品婷婷| 青春草视频在线免费观看| 久久久久久久久大av| 99热这里只有精品一区| 精品国产三级普通话版| 亚洲不卡免费看| 五月伊人婷婷丁香| 成熟少妇高潮喷水视频| 五月伊人婷婷丁香| 老司机福利观看| 成人三级黄色视频| 国产精品久久久久久久久免| 国产 一区精品| 久久精品久久久久久久性| 又粗又爽又猛毛片免费看| 日韩亚洲欧美综合| 日韩欧美 国产精品| 精品欧美国产一区二区三| 成人特级黄色片久久久久久久| 51国产日韩欧美| 亚洲成人av在线免费| 夜夜爽天天搞| 亚洲成人av在线免费| 亚洲高清免费不卡视频| 久久久久久久久久久免费av| 国产精品久久久久久av不卡| 国产高清不卡午夜福利| 亚洲精品久久国产高清桃花| 国产精品久久久久久久电影| 久久九九热精品免费| 国产精品久久久久久久电影| 一本一本综合久久| 天堂中文最新版在线下载 | 久久久久久久久久久免费av| 好男人视频免费观看在线| 亚洲经典国产精华液单| 国产精品久久视频播放| 夜夜看夜夜爽夜夜摸| 亚洲精品国产av成人精品| 九九在线视频观看精品| 国产亚洲av嫩草精品影院| 国产熟女欧美一区二区| 在现免费观看毛片| 国内少妇人妻偷人精品xxx网站| 亚洲国产色片| 亚洲激情五月婷婷啪啪| 久久精品国产亚洲网站| 亚洲精品乱码久久久久久按摩| 一本久久中文字幕| 欧美+日韩+精品| 欧美精品一区二区大全| 欧美激情久久久久久爽电影| 亚洲欧洲国产日韩| 精品人妻熟女av久视频| 免费看光身美女| 久久99热这里只有精品18| 老司机影院成人| 91精品一卡2卡3卡4卡| 热99在线观看视频| 久久久久网色| 一个人看视频在线观看www免费| 午夜a级毛片| 久久久久九九精品影院| 99精品在免费线老司机午夜| 久久人人精品亚洲av| av卡一久久| 国产精品野战在线观看| 免费av毛片视频| 又爽又黄无遮挡网站| 成人毛片60女人毛片免费| 欧美在线一区亚洲| 国产在线男女| 久久中文看片网| 午夜福利成人在线免费观看| 亚洲精品日韩在线中文字幕 | av在线亚洲专区| 激情 狠狠 欧美| 日韩在线高清观看一区二区三区| 一个人免费在线观看电影| 一个人看视频在线观看www免费| 国产高清有码在线观看视频| 国产高清不卡午夜福利| 亚洲欧美日韩高清专用| 神马国产精品三级电影在线观看| 免费看日本二区| 久久精品国产自在天天线| 国产久久久一区二区三区| 亚洲综合色惰| 国产探花在线观看一区二区| 舔av片在线| 久久精品国产亚洲av涩爱 | 国产免费一级a男人的天堂| 国产老妇伦熟女老妇高清| 大型黄色视频在线免费观看| 免费看av在线观看网站| 色哟哟哟哟哟哟| 亚洲高清免费不卡视频| 伊人久久精品亚洲午夜| 亚洲av中文av极速乱| 3wmmmm亚洲av在线观看| a级毛色黄片| 国产精品福利在线免费观看| 春色校园在线视频观看| 白带黄色成豆腐渣| 亚洲av一区综合| 免费av观看视频| 99热只有精品国产| 九九爱精品视频在线观看| 免费观看a级毛片全部| 99热这里只有是精品在线观看| 22中文网久久字幕| 99久久无色码亚洲精品果冻| 欧美日本视频| 欧美性感艳星| 亚洲国产日韩欧美精品在线观看| 变态另类丝袜制服| 午夜久久久久精精品| 欧美bdsm另类| 秋霞在线观看毛片| 天堂中文最新版在线下载 | 中文字幕熟女人妻在线| 午夜精品一区二区三区免费看| 国产蜜桃级精品一区二区三区| 国产一级毛片在线| 国产av在哪里看| 国产成人aa在线观看| 国产一级毛片七仙女欲春2| videossex国产| 看十八女毛片水多多多| 日本av手机在线免费观看| 小蜜桃在线观看免费完整版高清| 欧美精品国产亚洲| 一本精品99久久精品77| 男女下面进入的视频免费午夜| 久久草成人影院| 国产 一区精品| 成人鲁丝片一二三区免费| av天堂在线播放| 欧美xxxx性猛交bbbb| ponron亚洲| 在线免费十八禁| 秋霞在线观看毛片| 国产精品一区二区三区四区免费观看| 日本av手机在线免费观看| 成人美女网站在线观看视频| 欧美最黄视频在线播放免费| 国产乱人视频| 午夜福利视频1000在线观看| 国产白丝娇喘喷水9色精品| 卡戴珊不雅视频在线播放| 国产精品不卡视频一区二区| 久久久精品欧美日韩精品| 精品久久久久久久人妻蜜臀av| 麻豆国产97在线/欧美| 天堂网av新在线| 99久久中文字幕三级久久日本| 国产av一区在线观看免费| 菩萨蛮人人尽说江南好唐韦庄 | 午夜福利视频1000在线观看| 国产精品乱码一区二三区的特点| 久久精品综合一区二区三区| av在线老鸭窝| www.色视频.com| 国产亚洲av片在线观看秒播厂 | 精品欧美国产一区二区三| 精品熟女少妇av免费看| 亚洲欧美日韩卡通动漫| 免费人成在线观看视频色| 少妇被粗大猛烈的视频| 乱系列少妇在线播放| 国产精品1区2区在线观看.| 免费搜索国产男女视频| 色综合色国产| 欧美日韩乱码在线| 国产精品野战在线观看| av又黄又爽大尺度在线免费看 | 久久6这里有精品| 日韩欧美精品免费久久| 色尼玛亚洲综合影院| 国产一区二区三区av在线 | 1000部很黄的大片| 午夜福利视频1000在线观看| 欧美区成人在线视频| 69av精品久久久久久| av专区在线播放| 有码 亚洲区| 国产在线精品亚洲第一网站| 男人舔女人下体高潮全视频| 国产精品蜜桃在线观看 | 毛片一级片免费看久久久久| 免费av不卡在线播放| 久久韩国三级中文字幕| 国产精品1区2区在线观看.| 亚洲高清免费不卡视频| 变态另类成人亚洲欧美熟女| 国产精品嫩草影院av在线观看| 中文欧美无线码| 最近最新中文字幕大全电影3| 中文字幕熟女人妻在线| 日本成人三级电影网站| 国产日本99.免费观看| 久久久久久伊人网av| 免费观看精品视频网站| 床上黄色一级片| 中文字幕av在线有码专区| 一级二级三级毛片免费看| 熟妇人妻久久中文字幕3abv| 亚洲综合色惰| 99久久中文字幕三级久久日本| 亚洲精品456在线播放app| 欧美+亚洲+日韩+国产| 可以在线观看的亚洲视频| 欧美一区二区精品小视频在线| 99久国产av精品国产电影| 色综合色国产| 天天躁日日操中文字幕| 亚洲在线自拍视频| 国产精品一区二区三区四区久久| 又黄又爽又刺激的免费视频.| 亚洲av免费在线观看| 国产大屁股一区二区在线视频| 久久精品国产亚洲av香蕉五月| 内地一区二区视频在线| 26uuu在线亚洲综合色| 晚上一个人看的免费电影| 日韩欧美精品v在线| 婷婷色av中文字幕| 亚洲国产精品sss在线观看| 国产单亲对白刺激| 亚洲精品自拍成人| 麻豆乱淫一区二区| 最好的美女福利视频网| av在线播放精品| 亚洲色图av天堂| 亚洲七黄色美女视频| 国产精品1区2区在线观看.| 91在线精品国自产拍蜜月| 在线观看免费视频日本深夜| av女优亚洲男人天堂| 精品人妻视频免费看| 国产大屁股一区二区在线视频| 又粗又硬又长又爽又黄的视频 | 国产高清不卡午夜福利| 春色校园在线视频观看| 亚洲av二区三区四区| 日日摸夜夜添夜夜爱| 色综合色国产| 能在线免费看毛片的网站| 少妇高潮的动态图| 好男人视频免费观看在线| 国产熟女欧美一区二区| 男人和女人高潮做爰伦理| 一进一出抽搐gif免费好疼| 国产人妻一区二区三区在| 亚洲精品乱码久久久久久按摩| 国产黄色视频一区二区在线观看 | avwww免费| 欧美另类亚洲清纯唯美| 免费av毛片视频| 亚洲在线观看片| 男女边吃奶边做爰视频| 国产精品.久久久| 久久精品国产亚洲网站| 国产精品人妻久久久影院| 少妇熟女aⅴ在线视频| 欧美xxxx性猛交bbbb| 国产精品一及| 日韩,欧美,国产一区二区三区 | 99久国产av精品| 2022亚洲国产成人精品| 国产亚洲91精品色在线| 看片在线看免费视频| 日韩强制内射视频| 国产成人freesex在线| 国产高清有码在线观看视频| 久久久久久久久久久免费av| 在线观看午夜福利视频| 九九在线视频观看精品| 成人无遮挡网站| 日本-黄色视频高清免费观看| 最近2019中文字幕mv第一页| 91aial.com中文字幕在线观看| 免费黄网站久久成人精品| 我的女老师完整版在线观看| 久久这里有精品视频免费| 综合色丁香网| 亚洲成人av在线免费| 在线免费十八禁| 精品人妻视频免费看| 一本精品99久久精品77| 一级毛片电影观看 | 亚洲成a人片在线一区二区| 搡女人真爽免费视频火全软件| 啦啦啦观看免费观看视频高清| 18禁在线播放成人免费| 波野结衣二区三区在线|