• <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à)論述
    伊人久久精品亚洲午夜| 在线观看av片永久免费下载| 亚洲色图av天堂| 91狼人影院| 国产伦一二天堂av在线观看| 尾随美女入室| 国产又黄又爽又无遮挡在线| 99热这里只有精品一区| 国产单亲对白刺激| 亚洲国产成人一精品久久久| 高清毛片免费看| 特大巨黑吊av在线直播| 午夜激情福利司机影院| 亚洲欧美中文字幕日韩二区| 插阴视频在线观看视频| 老司机影院成人| 久久精品国产99精品国产亚洲性色| 秋霞伦理黄片| 精品熟女少妇av免费看| 日韩大片免费观看网站 | 国产黄色视频一区二区在线观看 | 亚洲精品乱久久久久久| 国产在视频线精品| 最近最新中文字幕大全电影3| 中国美白少妇内射xxxbb| 黄色配什么色好看| 免费av观看视频| 内射极品少妇av片p| 伦理电影大哥的女人| 小蜜桃在线观看免费完整版高清| a级毛色黄片| 国产国拍精品亚洲av在线观看| 小说图片视频综合网站| 国产在线一区二区三区精 | 精品欧美国产一区二区三| 国产视频首页在线观看| 亚洲国产精品合色在线| 中国国产av一级| 美女国产视频在线观看| 亚洲激情五月婷婷啪啪| 国语自产精品视频在线第100页| 天堂网av新在线| 免费看av在线观看网站| 亚洲最大成人手机在线| 听说在线观看完整版免费高清| 国产精品国产三级专区第一集| 成人综合一区亚洲| 成人欧美大片| 久久久久久久久久成人| av天堂中文字幕网| 日本免费一区二区三区高清不卡| 黄片无遮挡物在线观看| 久久婷婷人人爽人人干人人爱| 国产亚洲91精品色在线| 一级二级三级毛片免费看| 91在线精品国自产拍蜜月| av国产久精品久网站免费入址| videossex国产| 乱人视频在线观看| av卡一久久| 久久久久网色| 成人二区视频| 亚洲人成网站高清观看| 欧美日本视频| 国产精品1区2区在线观看.| 1000部很黄的大片| 亚洲精品自拍成人| av国产免费在线观看| 亚洲一级一片aⅴ在线观看| 亚洲国产欧洲综合997久久,| 亚洲精品国产成人久久av| 丰满人妻一区二区三区视频av| 天天躁夜夜躁狠狠久久av| 老司机影院成人| 春色校园在线视频观看| 国产精品久久久久久精品电影小说 | 亚洲av二区三区四区| 精品熟女少妇av免费看| 国产三级在线视频| 亚洲国产精品专区欧美| 亚洲国产欧美人成| 99国产精品一区二区蜜桃av| 免费观看的影片在线观看| 韩国av在线不卡| 在线播放国产精品三级| 在线播放无遮挡| h日本视频在线播放| 国产精品av视频在线免费观看| 亚洲欧美日韩卡通动漫| 最后的刺客免费高清国语| 日韩欧美国产在线观看| 国内少妇人妻偷人精品xxx网站| 亚洲国产欧美在线一区| 国产精品美女特级片免费视频播放器| 国产午夜精品论理片| 国产一区二区三区av在线| 观看美女的网站| 国产精品久久久久久久久免| 人人妻人人澡人人爽人人夜夜 | 精品午夜福利在线看| 国产亚洲av嫩草精品影院| 久久久a久久爽久久v久久| 亚洲第一区二区三区不卡| av又黄又爽大尺度在线免费看 | 一个人看的www免费观看视频| 亚洲av成人精品一二三区| 日本av手机在线免费观看| 成人二区视频| 网址你懂的国产日韩在线| 久久99热这里只有精品18| 日韩欧美精品免费久久| 国产欧美日韩精品一区二区| 蜜臀久久99精品久久宅男| 一级毛片久久久久久久久女| 欧美高清性xxxxhd video| 国产精品1区2区在线观看.| 久久久久久伊人网av| 最近手机中文字幕大全| 国产精品电影一区二区三区| 美女大奶头视频| 最近中文字幕高清免费大全6| 午夜精品国产一区二区电影 | 国产av一区在线观看免费| 全区人妻精品视频| 久久久久久大精品| 久久久久久久国产电影| 成人毛片60女人毛片免费| 成年女人永久免费观看视频| 免费看日本二区| 免费无遮挡裸体视频| www.色视频.com| 18禁裸乳无遮挡免费网站照片| 亚洲在线观看片| 国产精品综合久久久久久久免费| 国产黄片视频在线免费观看| 国产一区有黄有色的免费视频 | 国产精品无大码| 99国产精品一区二区蜜桃av| 99热精品在线国产| av免费观看日本| 欧美xxxx黑人xx丫x性爽| 老女人水多毛片| 精品国产一区二区三区久久久樱花 | 久久韩国三级中文字幕| videos熟女内射| av在线蜜桃| 激情 狠狠 欧美| 国产v大片淫在线免费观看| 夫妻性生交免费视频一级片| 男女国产视频网站| 亚洲人成网站在线播| 亚洲成人av在线免费| 国产成年人精品一区二区| 麻豆成人午夜福利视频| 人妻制服诱惑在线中文字幕| 国产精品一区二区三区四区久久| 成人鲁丝片一二三区免费| 国产毛片a区久久久久| 99热网站在线观看| 国产成人aa在线观看| 国产成人一区二区在线| 九九久久精品国产亚洲av麻豆| 国产亚洲av嫩草精品影院| 免费一级毛片在线播放高清视频| 国产不卡一卡二| 男人的好看免费观看在线视频| 亚洲天堂国产精品一区在线| 女人久久www免费人成看片 | 亚洲美女视频黄频| 三级国产精品欧美在线观看| 亚洲国产欧美在线一区| 精品欧美国产一区二区三| 精品99又大又爽又粗少妇毛片| 一个人看视频在线观看www免费| 亚洲精品久久久久久婷婷小说 | 国产国拍精品亚洲av在线观看| 亚洲av熟女| 色视频www国产| 久久久久久久久中文| 又粗又爽又猛毛片免费看| 国产成人91sexporn| 久久韩国三级中文字幕| 春色校园在线视频观看| 国产亚洲av片在线观看秒播厂 | 欧美不卡视频在线免费观看| 久久久午夜欧美精品| 久久久久久大精品| 久久久久性生活片| 国产乱人偷精品视频| 白带黄色成豆腐渣| 日韩一本色道免费dvd| 久久久久久大精品| 婷婷色av中文字幕| 久久鲁丝午夜福利片| 亚洲精品aⅴ在线观看| 99久久中文字幕三级久久日本| 欧美一区二区国产精品久久精品| 久99久视频精品免费| 啦啦啦韩国在线观看视频| ponron亚洲| 男女边吃奶边做爰视频| 热99re8久久精品国产| 99久久精品一区二区三区| 最近中文字幕2019免费版| 亚洲精品成人久久久久久| 搡老妇女老女人老熟妇| 国产中年淑女户外野战色| 日韩视频在线欧美| 三级国产精品欧美在线观看| 亚洲精品乱久久久久久| 一区二区三区四区激情视频| 国产精品美女特级片免费视频播放器| 免费看av在线观看网站| 亚洲欧美清纯卡通| 亚洲婷婷狠狠爱综合网| 毛片女人毛片| 国产v大片淫在线免费观看| 国产高清三级在线| 少妇的逼水好多| 伊人久久精品亚洲午夜| 免费观看在线日韩| 人人妻人人澡欧美一区二区| 99久久精品热视频| 观看美女的网站| 一区二区三区乱码不卡18| 三级毛片av免费| 老司机福利观看| 青春草视频在线免费观看| 日韩,欧美,国产一区二区三区 | 欧美xxxx黑人xx丫x性爽| 色尼玛亚洲综合影院| av福利片在线观看| 亚洲av.av天堂| 人妻少妇偷人精品九色| 国产淫语在线视频| 精品一区二区三区人妻视频| 精品久久久久久电影网 | 国产伦在线观看视频一区| 色网站视频免费| 婷婷色综合大香蕉| 国产视频首页在线观看| 久久久久网色| 亚洲国产精品久久男人天堂| 亚洲国产日韩欧美精品在线观看| 精品不卡国产一区二区三区| 亚洲精品乱码久久久久久按摩| 国产精品人妻久久久影院| 嫩草影院入口| 日本黄色片子视频| 99久久精品一区二区三区| 国产免费福利视频在线观看| 热99在线观看视频| 少妇裸体淫交视频免费看高清| 好男人在线观看高清免费视频| 亚洲图色成人| 日韩一区二区视频免费看| 日本-黄色视频高清免费观看| 日韩中字成人| 免费大片18禁| 麻豆成人午夜福利视频| 99热这里只有是精品50| 免费av不卡在线播放| 日日撸夜夜添| 一区二区三区免费毛片| 精品久久久噜噜| 日韩av不卡免费在线播放| 18禁在线播放成人免费| 中文字幕免费在线视频6| www.av在线官网国产| 亚洲精品日韩在线中文字幕| 久久精品久久精品一区二区三区| 一级av片app| 久久这里有精品视频免费| 一个人看的www免费观看视频| 在线观看av片永久免费下载| 久久这里只有精品中国| 噜噜噜噜噜久久久久久91| 亚洲婷婷狠狠爱综合网| 国产中年淑女户外野战色| 别揉我奶头 嗯啊视频| 国产成人午夜福利电影在线观看| 成年版毛片免费区| 边亲边吃奶的免费视频| 亚洲成人中文字幕在线播放| 亚洲熟妇中文字幕五十中出| 日本av手机在线免费观看| 久久久色成人| 夜夜爽夜夜爽视频| av.在线天堂| 国产一区二区在线观看日韩| 欧美一区二区国产精品久久精品| 国产精品美女特级片免费视频播放器| 欧美97在线视频| 国产极品精品免费视频能看的| АⅤ资源中文在线天堂| 老司机影院成人| 级片在线观看| 人体艺术视频欧美日本| 国产精品久久久久久精品电影| 午夜免费激情av| 99热这里只有精品一区| 成人高潮视频无遮挡免费网站| 哪个播放器可以免费观看大片| 亚洲av熟女| 日韩欧美国产在线观看| 黄色日韩在线| eeuss影院久久| 我要看日韩黄色一级片| 女的被弄到高潮叫床怎么办| 日本熟妇午夜| 亚洲自拍偷在线| 久久久精品94久久精品| 国产精品国产三级国产专区5o | 国产爱豆传媒在线观看| 国产精品,欧美在线| 日本欧美国产在线视频| 日韩欧美国产在线观看| 三级经典国产精品| 久久99精品国语久久久| 精品午夜福利在线看| 男女视频在线观看网站免费| 中文字幕精品亚洲无线码一区| 亚洲人成网站在线播| 国产精品一区二区三区四区免费观看| 国产一区二区三区av在线| 99国产精品一区二区蜜桃av| 免费黄色在线免费观看| 久久久精品大字幕| 日本免费在线观看一区| 亚洲av成人精品一二三区| 欧美xxxx黑人xx丫x性爽| 青春草国产在线视频| 尾随美女入室| 国产午夜精品久久久久久一区二区三区| 亚洲最大成人av| 18+在线观看网站| 色综合色国产| 精品一区二区免费观看| 性插视频无遮挡在线免费观看| av在线老鸭窝| 成人毛片a级毛片在线播放| 国产视频首页在线观看| 人妻系列 视频| av在线亚洲专区| 亚洲美女视频黄频| 国产真实乱freesex| 中文字幕久久专区| 精品国产露脸久久av麻豆 | 国产成人精品一,二区| 国产白丝娇喘喷水9色精品| 亚洲精品日韩av片在线观看| 国产亚洲最大av| 秋霞在线观看毛片| 色5月婷婷丁香| 日本与韩国留学比较| av在线天堂中文字幕| 日本黄色片子视频| 国产午夜精品论理片| 日产精品乱码卡一卡2卡三| 国产精品女同一区二区软件| 免费在线观看成人毛片| 直男gayav资源| 精品少妇黑人巨大在线播放 | 日本av手机在线免费观看| 国产麻豆成人av免费视频| 日韩一区二区视频免费看| 国产av在哪里看| 丝袜美腿在线中文| 真实男女啪啪啪动态图| 亚洲精品久久久久久婷婷小说 | 婷婷色av中文字幕| 波野结衣二区三区在线| 亚洲国产欧洲综合997久久,| 能在线免费看毛片的网站| 在线播放国产精品三级| av福利片在线观看| 国产精品久久久久久精品电影| 亚洲综合精品二区| 最近最新中文字幕大全电影3| 一级黄色大片毛片| av又黄又爽大尺度在线免费看 | 男人狂女人下面高潮的视频| 国产亚洲av片在线观看秒播厂 | 国产高清有码在线观看视频| 日本一本二区三区精品| 国产精品99久久久久久久久| 久久热精品热| 久久久久久久久久黄片| 久久综合国产亚洲精品| 亚洲精品456在线播放app| 内射极品少妇av片p| 欧美日韩精品成人综合77777| 成人午夜高清在线视频| 又爽又黄a免费视频| 日本免费在线观看一区| av又黄又爽大尺度在线免费看 | av线在线观看网站| 欧美三级亚洲精品| 亚洲国产欧洲综合997久久,| 热99re8久久精品国产| 大香蕉久久网| 日韩三级伦理在线观看| 中文字幕精品亚洲无线码一区| 人人妻人人看人人澡| 在线观看一区二区三区| 最后的刺客免费高清国语| 天堂av国产一区二区熟女人妻| 深夜a级毛片| 久久午夜福利片| 免费观看精品视频网站| 我的老师免费观看完整版| 亚洲欧洲日产国产| 久久精品国产鲁丝片午夜精品| 国产爱豆传媒在线观看| 69av精品久久久久久| 精品久久久久久久久久久久久| 国产精品久久久久久久电影| 亚洲精品乱久久久久久| 99在线视频只有这里精品首页| 一级毛片aaaaaa免费看小| 免费看a级黄色片| 春色校园在线视频观看| 十八禁国产超污无遮挡网站| 熟女电影av网| 黑人高潮一二区| 亚洲人与动物交配视频| 国产精品不卡视频一区二区| 久久这里只有精品中国| 亚洲国产精品成人久久小说| 欧美日韩国产亚洲二区| 免费观看a级毛片全部| 国产精华一区二区三区| 人人妻人人看人人澡| 91av网一区二区| 欧美成人a在线观看| 欧美变态另类bdsm刘玥| 国产一区二区三区av在线| videossex国产| 91狼人影院| 国产亚洲91精品色在线| 午夜福利成人在线免费观看| 成人毛片60女人毛片免费| 搡老妇女老女人老熟妇| 久久久精品大字幕| 亚洲精品成人久久久久久| 成人av在线播放网站| 九草在线视频观看| 亚洲欧洲日产国产| 色尼玛亚洲综合影院| 自拍偷自拍亚洲精品老妇| 联通29元200g的流量卡| 亚洲欧美日韩高清专用| av在线亚洲专区| 97超视频在线观看视频| 午夜激情福利司机影院| 欧美激情久久久久久爽电影| 亚洲一级一片aⅴ在线观看| 午夜激情欧美在线| 午夜福利在线观看免费完整高清在| 激情 狠狠 欧美| 亚洲精品乱码久久久v下载方式| 国产爱豆传媒在线观看| 男插女下体视频免费在线播放| 成年版毛片免费区| 狂野欧美激情性xxxx在线观看| 色尼玛亚洲综合影院| 午夜日本视频在线| 成人午夜高清在线视频| 99热这里只有是精品50| 欧美日本亚洲视频在线播放| 一区二区三区免费毛片| 国产午夜精品一二区理论片| 男人舔女人下体高潮全视频| 1024手机看黄色片| 你懂的网址亚洲精品在线观看 | 在线播放无遮挡| 日韩大片免费观看网站 | 日韩成人伦理影院| 国产成人免费观看mmmm| 国产乱来视频区| 夜夜看夜夜爽夜夜摸| 搡老妇女老女人老熟妇| 国产精品电影一区二区三区| 久久婷婷人人爽人人干人人爱| 久99久视频精品免费| 亚洲成色77777| 国产亚洲午夜精品一区二区久久 | 国产爱豆传媒在线观看| 亚洲av熟女| 日韩一区二区视频免费看| 免费播放大片免费观看视频在线观看 | 久久精品国产自在天天线| 亚洲欧美中文字幕日韩二区| 国产探花极品一区二区| 国产男人的电影天堂91| 久热久热在线精品观看| 三级国产精品欧美在线观看| 日本爱情动作片www.在线观看| 级片在线观看| 久久欧美精品欧美久久欧美| 午夜福利在线观看吧| 人人妻人人看人人澡| 亚洲综合色惰| 色视频www国产| 亚洲精品aⅴ在线观看| 男女视频在线观看网站免费| 亚洲av成人精品一区久久| 久久精品国产自在天天线| 久久精品夜色国产| 99视频精品全部免费 在线| 综合色丁香网| 久久这里有精品视频免费| 久久婷婷人人爽人人干人人爱| 久久99热6这里只有精品| 国产午夜福利久久久久久| 日本免费a在线| 99久久精品一区二区三区| 97超视频在线观看视频| 免费黄网站久久成人精品| 亚洲欧美精品专区久久| 简卡轻食公司| 日本黄大片高清| 亚洲va在线va天堂va国产| 白带黄色成豆腐渣| 免费大片18禁| 国产在线男女| 国产一区亚洲一区在线观看| 色尼玛亚洲综合影院| 91av网一区二区| 午夜激情福利司机影院| 久久精品91蜜桃| 99久久无色码亚洲精品果冻| 内射极品少妇av片p| a级毛片免费高清观看在线播放| 国产精品麻豆人妻色哟哟久久 | 免费看a级黄色片| 一本一本综合久久| 日本黄色视频三级网站网址| 嫩草影院入口| 亚洲国产精品合色在线| 免费大片18禁| 熟妇人妻久久中文字幕3abv| 国产精品精品国产色婷婷| 国产精品野战在线观看| 中文精品一卡2卡3卡4更新| 精品久久久久久成人av| 国产午夜精品久久久久久一区二区三区| 最近的中文字幕免费完整| 国内少妇人妻偷人精品xxx网站| 免费观看的影片在线观看| 夜夜看夜夜爽夜夜摸| 亚洲四区av| 国产精品久久久久久久电影| 六月丁香七月| 国产精品.久久久| 久久6这里有精品| 国产成人freesex在线| 亚洲欧美日韩卡通动漫| 联通29元200g的流量卡| 国产亚洲91精品色在线| 国产激情偷乱视频一区二区| 亚洲内射少妇av| 成人综合一区亚洲| 黄片wwwwww| 久久99热6这里只有精品| 又粗又爽又猛毛片免费看| 午夜精品在线福利| 九九在线视频观看精品| 成人午夜高清在线视频| 国产成人一区二区在线| 91久久精品国产一区二区三区| 国产成人精品一,二区| 色综合亚洲欧美另类图片| 18禁动态无遮挡网站| 久久99蜜桃精品久久| 99热6这里只有精品| 精品人妻偷拍中文字幕| 一本一本综合久久| 日本黄色视频三级网站网址| 好男人视频免费观看在线| 人妻制服诱惑在线中文字幕| 美女xxoo啪啪120秒动态图| 美女被艹到高潮喷水动态| 亚洲av成人av| 天堂影院成人在线观看| 搡女人真爽免费视频火全软件| 中文字幕亚洲精品专区| 亚洲精品日韩av片在线观看| 桃色一区二区三区在线观看| 久久久午夜欧美精品| 国产美女午夜福利| 中文天堂在线官网| 亚洲国产精品成人久久小说| 三级国产精品欧美在线观看| 禁无遮挡网站| 欧美日韩一区二区视频在线观看视频在线 | 男人舔女人下体高潮全视频| 极品教师在线视频| 91av网一区二区| 国产91av在线免费观看| 亚洲av.av天堂| 一级爰片在线观看| 国产国拍精品亚洲av在线观看| 日韩一区二区视频免费看| 国产一区二区在线av高清观看| 小蜜桃在线观看免费完整版高清| 精品午夜福利在线看| 午夜精品国产一区二区电影 | 亚洲精品国产成人久久av| 午夜免费男女啪啪视频观看| 亚洲最大成人av| 欧美性感艳星| 国产免费男女视频| 久久精品久久精品一区二区三区| 九草在线视频观看|