• <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à)論述
    精品一区二区三区人妻视频| 99国产综合亚洲精品| 一区二区三区免费毛片| 精品久久久久久成人av| 又黄又粗又硬又大视频| 欧美日本亚洲视频在线播放| 亚洲性夜色夜夜综合| 欧美成狂野欧美在线观看| 国产久久久一区二区三区| 国产亚洲精品av在线| 欧美一区二区精品小视频在线| 99热只有精品国产| 男插女下体视频免费在线播放| 色综合欧美亚洲国产小说| 久久精品国产综合久久久| 国产精品综合久久久久久久免费| 国产亚洲欧美在线一区二区| 99热6这里只有精品| 国内精品久久久久精免费| 中文亚洲av片在线观看爽| 国产主播在线观看一区二区| 午夜福利在线观看吧| 成熟少妇高潮喷水视频| 精品福利观看| 亚洲国产精品sss在线观看| 久久久久久国产a免费观看| 精品欧美国产一区二区三| 午夜久久久久精精品| 国产黄片美女视频| 亚洲激情在线av| 欧美黑人巨大hd| 成年女人毛片免费观看观看9| 亚洲国产欧美网| 亚洲专区中文字幕在线| 精品久久久久久久毛片微露脸| 欧美色视频一区免费| 九色国产91popny在线| 三级男女做爰猛烈吃奶摸视频| 日韩 欧美 亚洲 中文字幕| 亚洲av电影在线进入| 精品熟女少妇八av免费久了| 久久精品国产综合久久久| 人人妻,人人澡人人爽秒播| 久久久久久久午夜电影| e午夜精品久久久久久久| 亚洲国产欧美网| 久久欧美精品欧美久久欧美| 国内毛片毛片毛片毛片毛片| 特大巨黑吊av在线直播| 亚洲精品一卡2卡三卡4卡5卡| eeuss影院久久| 一夜夜www| 国产成人欧美在线观看| 老司机在亚洲福利影院| 久久精品夜夜夜夜夜久久蜜豆| 性色av乱码一区二区三区2| 国产69精品久久久久777片| 国产午夜精品论理片| av中文乱码字幕在线| 国产一区二区在线av高清观看| 男女视频在线观看网站免费| 黄色丝袜av网址大全| 国产精品98久久久久久宅男小说| 精品乱码久久久久久99久播| 国产真实乱freesex| 国产精品电影一区二区三区| 国产野战对白在线观看| 啪啪无遮挡十八禁网站| 久久国产精品人妻蜜桃| 在线观看av片永久免费下载| 亚洲av电影在线进入| bbb黄色大片| 看黄色毛片网站| 嫩草影视91久久| 亚洲av电影在线进入| xxxwww97欧美| 亚洲国产精品合色在线| 国产免费av片在线观看野外av| 他把我摸到了高潮在线观看| 久久精品人妻少妇| 久久精品人妻少妇| 午夜福利18| 国产蜜桃级精品一区二区三区| 蜜桃久久精品国产亚洲av| 麻豆一二三区av精品| 成人av一区二区三区在线看| 午夜a级毛片| 丰满乱子伦码专区| 99在线视频只有这里精品首页| 蜜桃久久精品国产亚洲av| 18+在线观看网站| 夜夜看夜夜爽夜夜摸| www.色视频.com| 久久亚洲真实| 亚洲精品乱码久久久v下载方式 | 国产精品亚洲一级av第二区| 成人三级黄色视频| 搡女人真爽免费视频火全软件 | 舔av片在线| www.www免费av| 成人18禁在线播放| 99热精品在线国产| 久久国产精品影院| 少妇丰满av| 亚洲国产欧洲综合997久久,| 18禁在线播放成人免费| 女警被强在线播放| 亚洲人成网站高清观看| 亚洲av五月六月丁香网| 成人三级黄色视频| 嫩草影院精品99| 俺也久久电影网| 中文在线观看免费www的网站| 国产激情偷乱视频一区二区| 国产精品一区二区三区四区免费观看 | 国产精品乱码一区二三区的特点| 国产一级毛片七仙女欲春2| 国产免费av片在线观看野外av| 18美女黄网站色大片免费观看| www.熟女人妻精品国产| 色老头精品视频在线观看| 国产午夜福利久久久久久| 18禁黄网站禁片免费观看直播| 香蕉av资源在线| 哪里可以看免费的av片| 欧美性猛交╳xxx乱大交人| 久久九九热精品免费| 熟女少妇亚洲综合色aaa.| 日韩大尺度精品在线看网址| 免费看光身美女| 欧美日本亚洲视频在线播放| 19禁男女啪啪无遮挡网站| 日本 欧美在线| 亚洲人成网站在线播| 手机成人av网站| 日韩欧美一区二区三区在线观看| 国产精品久久视频播放| 精品人妻一区二区三区麻豆 | 亚洲精品成人久久久久久| 18禁在线播放成人免费| 真人一进一出gif抽搐免费| 免费在线观看成人毛片| 国产又黄又爽又无遮挡在线| 日本一本二区三区精品| 欧美日韩国产亚洲二区| av在线蜜桃| 国产高清激情床上av| 观看美女的网站| 三级毛片av免费| 99热这里只有是精品50| 国内精品一区二区在线观看| 此物有八面人人有两片| 亚洲国产欧洲综合997久久,| 性色avwww在线观看| 欧美中文综合在线视频| 久久精品亚洲精品国产色婷小说| 91麻豆精品激情在线观看国产| 久久久国产成人免费| 韩国av一区二区三区四区| 五月伊人婷婷丁香| 欧美日韩福利视频一区二区| 一个人看视频在线观看www免费 | 久久精品亚洲精品国产色婷小说| 中文字幕人妻丝袜一区二区| 色av中文字幕| 日韩精品青青久久久久久| 欧美成人性av电影在线观看| 可以在线观看毛片的网站| 2021天堂中文幕一二区在线观| 欧美黑人巨大hd| 免费在线观看日本一区| 亚洲片人在线观看| 在线免费观看不下载黄p国产 | 在线播放无遮挡| avwww免费| 亚洲欧美日韩高清专用| 亚洲国产精品合色在线| 欧美性猛交黑人性爽| 熟女少妇亚洲综合色aaa.| 身体一侧抽搐| 男女视频在线观看网站免费| 欧美日本视频| 亚洲18禁久久av| 桃色一区二区三区在线观看| 亚洲中文字幕一区二区三区有码在线看| 成人特级av手机在线观看| 色综合欧美亚洲国产小说| 亚洲中文字幕日韩| 亚洲av第一区精品v没综合| 午夜久久久久精精品| 亚洲av成人不卡在线观看播放网| 亚洲精品影视一区二区三区av| 看免费av毛片| 在线观看一区二区三区| 国产成人aa在线观看| 久久久久国产精品人妻aⅴ院| 欧美黄色片欧美黄色片| 少妇的丰满在线观看| 动漫黄色视频在线观看| 免费看美女性在线毛片视频| 国产av麻豆久久久久久久| 日本与韩国留学比较| 在线观看午夜福利视频| 午夜福利高清视频| 免费搜索国产男女视频| 欧美乱码精品一区二区三区| 天堂动漫精品| 久久伊人香网站| 午夜福利高清视频| 免费观看的影片在线观看| 国产午夜精品论理片| 久久中文看片网| 免费在线观看亚洲国产| 国产午夜福利久久久久久| 变态另类丝袜制服| 黄色片一级片一级黄色片| 亚洲国产欧美网| 国产真实伦视频高清在线观看 | 亚洲 国产 在线| 久久久久久久精品吃奶| 啦啦啦韩国在线观看视频| 久久精品国产自在天天线| 成人av在线播放网站| 少妇熟女aⅴ在线视频| 国产精品一区二区三区四区久久| 成年女人毛片免费观看观看9| 日本黄大片高清| 欧美成人性av电影在线观看| 成人一区二区视频在线观看| 亚洲激情在线av| 最新在线观看一区二区三区| 国产乱人视频| 精品欧美国产一区二区三| 成人高潮视频无遮挡免费网站| 亚洲内射少妇av| 国产成人啪精品午夜网站| 大型黄色视频在线免费观看| 亚洲五月婷婷丁香| 国产精品亚洲美女久久久| 99热只有精品国产| www.熟女人妻精品国产| 中文字幕人成人乱码亚洲影| 亚洲欧美激情综合另类| 亚洲一区二区三区色噜噜| av视频在线观看入口| 看片在线看免费视频| 久久这里只有精品中国| 久久久久性生活片| 悠悠久久av| 91久久精品国产一区二区成人 | 亚洲av免费在线观看| 婷婷丁香在线五月| 亚洲成人精品中文字幕电影| 51国产日韩欧美| 变态另类成人亚洲欧美熟女| 69av精品久久久久久| 一区二区三区国产精品乱码| 国产91精品成人一区二区三区| 性欧美人与动物交配| 两个人的视频大全免费| 国产精华一区二区三区| av福利片在线观看| 我要搜黄色片| 久久久久久久久中文| 国产私拍福利视频在线观看| 内地一区二区视频在线| 97人妻精品一区二区三区麻豆| 国产淫片久久久久久久久 | 日韩欧美 国产精品| 中亚洲国语对白在线视频| 国产熟女xx| 亚洲av成人av| 九色国产91popny在线| 搡老岳熟女国产| 国产成人aa在线观看| 全区人妻精品视频| 中亚洲国语对白在线视频| 欧美乱色亚洲激情| 综合色av麻豆| 99久久精品一区二区三区| 18禁国产床啪视频网站| 色av中文字幕| 啦啦啦免费观看视频1| 热99re8久久精品国产| 精品一区二区三区视频在线观看免费| 一区二区三区免费毛片| 黄色日韩在线| 亚洲黑人精品在线| 欧美日本视频| 天天添夜夜摸| 在线国产一区二区在线| 女生性感内裤真人,穿戴方法视频| 国产精品永久免费网站| 日本 av在线| 又爽又黄无遮挡网站| 亚洲专区中文字幕在线| 国产亚洲精品久久久com| 激情在线观看视频在线高清| 女人十人毛片免费观看3o分钟| 中文字幕av在线有码专区| 97超视频在线观看视频| 久久国产精品人妻蜜桃| 久久中文看片网| 欧美中文综合在线视频| 精品国产三级普通话版| 欧美极品一区二区三区四区| 国产一区二区在线观看日韩 | 免费在线观看影片大全网站| 久久九九热精品免费| 久久国产精品影院| 中国美女看黄片| 婷婷亚洲欧美| 亚洲欧美日韩卡通动漫| 欧美乱码精品一区二区三区| 国产精品日韩av在线免费观看| 国产高清三级在线| 美女高潮喷水抽搐中文字幕| 国产 一区 欧美 日韩| 最新中文字幕久久久久| 最新在线观看一区二区三区| 亚洲美女黄片视频| 国产精品亚洲av一区麻豆| 嫩草影视91久久| 精品一区二区三区视频在线 | 很黄的视频免费| 久久久久久人人人人人| 午夜精品久久久久久毛片777| 精品国产亚洲在线| 欧美乱妇无乱码| 亚洲中文字幕日韩| 免费观看的影片在线观看| 在线观看66精品国产| 夜夜爽天天搞| 亚洲精品一卡2卡三卡4卡5卡| av视频在线观看入口| 国产探花极品一区二区| 午夜亚洲福利在线播放| 少妇熟女aⅴ在线视频| 日本免费a在线| 韩国av一区二区三区四区| 久久久久久大精品| 欧美日本视频| 日本免费a在线| 亚洲无线观看免费| 一个人免费在线观看电影| 久久久久免费精品人妻一区二区| 叶爱在线成人免费视频播放| 老司机深夜福利视频在线观看| 噜噜噜噜噜久久久久久91| 一区二区三区国产精品乱码| 舔av片在线| 在线国产一区二区在线| 欧美性感艳星| 国产精品乱码一区二三区的特点| 岛国在线免费视频观看| 成年版毛片免费区| 久久精品91无色码中文字幕| 熟女电影av网| 毛片女人毛片| ponron亚洲| 午夜精品久久久久久毛片777| 老汉色∧v一级毛片| 久久午夜亚洲精品久久| 国内少妇人妻偷人精品xxx网站| tocl精华| 91久久精品电影网| 亚洲成人久久爱视频| ponron亚洲| 亚洲成人久久爱视频| 在线观看av片永久免费下载| 国产成人av教育| 又粗又爽又猛毛片免费看| 一进一出抽搐动态| 欧美日本亚洲视频在线播放| 中出人妻视频一区二区| 国产不卡一卡二| 午夜福利成人在线免费观看| 久久精品国产亚洲av香蕉五月| 亚洲精品乱码久久久v下载方式 | 一夜夜www| 搞女人的毛片| 男女视频在线观看网站免费| 国产精品亚洲一级av第二区| 久久久色成人| 欧美+亚洲+日韩+国产| 精品无人区乱码1区二区| 日韩欧美国产一区二区入口| 女人被狂操c到高潮| 国产成人系列免费观看| 亚洲不卡免费看| 搡女人真爽免费视频火全软件 | 久久精品国产自在天天线| 一卡2卡三卡四卡精品乱码亚洲| 中文字幕人成人乱码亚洲影| 黄色丝袜av网址大全| 黄片小视频在线播放| 精品99又大又爽又粗少妇毛片 | 五月伊人婷婷丁香| 日本三级黄在线观看| 他把我摸到了高潮在线观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 中文字幕av成人在线电影| 女警被强在线播放| 亚洲精品在线观看二区| 亚洲精品456在线播放app | 色在线成人网| 亚洲久久久久久中文字幕| 婷婷六月久久综合丁香| 最新美女视频免费是黄的| 日本精品一区二区三区蜜桃| 中文字幕高清在线视频| 精品免费久久久久久久清纯| 精品福利观看| 久久午夜亚洲精品久久| 亚洲av成人不卡在线观看播放网| 欧美3d第一页| 久久婷婷人人爽人人干人人爱| 午夜久久久久精精品| 国产精品一区二区三区四区久久| 色综合婷婷激情| 国内毛片毛片毛片毛片毛片| 老司机在亚洲福利影院| 亚洲成av人片免费观看| 老汉色∧v一级毛片| bbb黄色大片| 国产野战对白在线观看| 亚洲不卡免费看| 精品乱码久久久久久99久播| 真人做人爱边吃奶动态| 国产黄色小视频在线观看| 亚洲色图av天堂| 啪啪无遮挡十八禁网站| 不卡一级毛片| 亚洲人与动物交配视频| 亚洲成人精品中文字幕电影| 国产精品永久免费网站| 国产国拍精品亚洲av在线观看 | 99精品欧美一区二区三区四区| 啦啦啦观看免费观看视频高清| 19禁男女啪啪无遮挡网站| 韩国av一区二区三区四区| 最近最新中文字幕大全免费视频| 在线天堂最新版资源| 久久久色成人| 午夜影院日韩av| 内地一区二区视频在线| 网址你懂的国产日韩在线| 最近最新中文字幕大全免费视频| 国产视频一区二区在线看| av女优亚洲男人天堂| 精品人妻1区二区| 色视频www国产| 一a级毛片在线观看| svipshipincom国产片| 深夜精品福利| 亚洲国产日韩欧美精品在线观看 | 精品一区二区三区视频在线观看免费| 午夜激情福利司机影院| 欧美在线黄色| 精品久久久久久久久久免费视频| 噜噜噜噜噜久久久久久91| 成人高潮视频无遮挡免费网站| 18美女黄网站色大片免费观看| 色哟哟哟哟哟哟| 老司机午夜十八禁免费视频| 国产一区二区在线观看日韩 | 热99在线观看视频| 脱女人内裤的视频| 免费观看人在逋| 亚洲精品一卡2卡三卡4卡5卡| 国产精品嫩草影院av在线观看 | 欧美乱色亚洲激情| 在线观看午夜福利视频| 好男人电影高清在线观看| 女人被狂操c到高潮| 国产精品av视频在线免费观看| 男女那种视频在线观看| 九色成人免费人妻av| 久久精品亚洲精品国产色婷小说| 久久久久国内视频| 国产乱人视频| 欧美另类亚洲清纯唯美| 国产精品一区二区三区四区久久| 成人永久免费在线观看视频| 色精品久久人妻99蜜桃| 伊人久久大香线蕉亚洲五| 欧美成人免费av一区二区三区| 丰满的人妻完整版| 中文字幕高清在线视频| 一个人观看的视频www高清免费观看| 日韩中文字幕欧美一区二区| 国产日本99.免费观看| 热99在线观看视频| 99久久精品一区二区三区| а√天堂www在线а√下载| 亚洲成人免费电影在线观看| 午夜激情福利司机影院| 国产精品av视频在线免费观看| 国产探花在线观看一区二区| 观看美女的网站| 亚洲人成网站在线播放欧美日韩| 狂野欧美白嫩少妇大欣赏| 亚洲国产精品999在线| 精品人妻一区二区三区麻豆 | 免费人成在线观看视频色| 少妇裸体淫交视频免费看高清| 99视频精品全部免费 在线| 在线免费观看不下载黄p国产 | 一级a爱片免费观看的视频| 真实男女啪啪啪动态图| 18禁在线播放成人免费| 午夜两性在线视频| 亚洲成a人片在线一区二区| 国产精品久久久久久精品电影| 真人做人爱边吃奶动态| 夜夜爽天天搞| 久久精品人妻少妇| 色尼玛亚洲综合影院| 床上黄色一级片| 亚洲国产精品合色在线| av女优亚洲男人天堂| 啦啦啦免费观看视频1| 天堂网av新在线| 一个人免费在线观看电影| 动漫黄色视频在线观看| 国产单亲对白刺激| 亚洲国产欧洲综合997久久,| 日本熟妇午夜| 一本一本综合久久| 亚洲国产精品sss在线观看| 亚洲人成网站高清观看| 天美传媒精品一区二区| av在线蜜桃| 特大巨黑吊av在线直播| 19禁男女啪啪无遮挡网站| 岛国在线观看网站| 可以在线观看的亚洲视频| 国产成人系列免费观看| 一级毛片女人18水好多| 91在线精品国自产拍蜜月 | 亚洲第一欧美日韩一区二区三区| 岛国视频午夜一区免费看| 国产免费男女视频| 1024手机看黄色片| 色噜噜av男人的天堂激情| 美女被艹到高潮喷水动态| 免费人成视频x8x8入口观看| 国产乱人视频| 高清在线国产一区| 99国产综合亚洲精品| 欧美三级亚洲精品| 亚洲第一欧美日韩一区二区三区| 国产国拍精品亚洲av在线观看 | 国内毛片毛片毛片毛片毛片| 国产极品精品免费视频能看的| 国产97色在线日韩免费| 天堂动漫精品| 婷婷亚洲欧美| 在线免费观看不下载黄p国产 | 中文字幕人成人乱码亚洲影| 日日摸夜夜添夜夜添小说| 嫩草影院入口| 免费看a级黄色片| 亚洲人成网站高清观看| 久久亚洲真实| 久久久久亚洲av毛片大全| 午夜久久久久精精品| 超碰av人人做人人爽久久 | 日韩欧美免费精品| 老汉色∧v一级毛片| 婷婷丁香在线五月| 一二三四社区在线视频社区8| 国产精品国产高清国产av| 亚洲人与动物交配视频| 国产欧美日韩一区二区三| 亚洲一区二区三区不卡视频| 午夜精品在线福利| 午夜福利欧美成人| 中文字幕高清在线视频| 国模一区二区三区四区视频| 久久欧美精品欧美久久欧美| 黑人欧美特级aaaaaa片| 亚洲欧美一区二区三区黑人| 中文字幕av成人在线电影| 精品人妻一区二区三区麻豆 | netflix在线观看网站| 国产亚洲精品久久久com| 亚洲精品成人久久久久久| 看片在线看免费视频| 99riav亚洲国产免费| 亚洲精品日韩av片在线观看 | 亚洲精品色激情综合| 国产精品国产高清国产av| 很黄的视频免费| 亚洲18禁久久av| 欧美日韩精品网址| 在线看三级毛片| 亚洲熟妇中文字幕五十中出| 婷婷精品国产亚洲av| 亚洲一区高清亚洲精品| 国产成人福利小说| 变态另类成人亚洲欧美熟女| 国产单亲对白刺激| 色综合站精品国产| 国产亚洲精品久久久久久毛片| 国产69精品久久久久777片| 99在线人妻在线中文字幕| 国产精品99久久99久久久不卡| 亚洲成a人片在线一区二区| 国产国拍精品亚洲av在线观看 | av视频在线观看入口| 欧美又色又爽又黄视频| 午夜精品一区二区三区免费看| 天天躁日日操中文字幕| 久久国产乱子伦精品免费另类|