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

    山東南四湖上級(jí)湖水質(zhì)變化評(píng)價(jià)(2008-2014年)及成因分析*

    2016-09-01 07:46:38解文靜婁山崇曹升樂孫秀玲王月敏
    湖泊科學(xué) 2016年3期
    關(guān)鍵詞:南四湖入湖監(jiān)測點(diǎn)

    解文靜,王 松,婁山崇,曹升樂,孫秀玲,王月敏,林 潔

    (1:山東大學(xué)土建與水利學(xué)院,濟(jì)南 250061)(2:山東省水文局,濟(jì)南 250002)

    ?

    山東南四湖上級(jí)湖水質(zhì)變化評(píng)價(jià)(2008-2014年)及成因分析*

    解文靜1,王松2,婁山崇2,曹升樂1,孫秀玲1,王月敏1,林潔1

    (1:山東大學(xué)土建與水利學(xué)院,濟(jì)南 250061)(2:山東省水文局,濟(jì)南 250002)

    基于南四湖上級(jí)湖6個(gè)監(jiān)測點(diǎn)2008-2014年42個(gè)奇數(shù)月份的水質(zhì)資料,選取溶解氧、化學(xué)需氧量(COD)、高錳酸鹽指數(shù)(CODMn)、五日生化需氧量(BOD5)、氨氮、總磷和氟化物作為水質(zhì)評(píng)價(jià)的指標(biāo). 對(duì)水質(zhì)指標(biāo)的歸一化方法進(jìn)行了改進(jìn),采用主成分分析法對(duì)6個(gè)監(jiān)測點(diǎn)的水質(zhì)綜合狀況進(jìn)行評(píng)價(jià),根據(jù)得出的各月綜合得分序列分析水質(zhì)變化趨勢. 與2000-2007年南四湖上級(jí)湖的水質(zhì)綜合狀況進(jìn)行比較, 2008-2014年水質(zhì)明顯變好. 評(píng)價(jià)結(jié)果表明,同一監(jiān)測點(diǎn)的水質(zhì)在不同時(shí)間點(diǎn)上有明顯變化,同一時(shí)間序列不同監(jiān)測點(diǎn)的水質(zhì)狀況也有不同變化. 此外,分析得出南四湖上級(jí)湖的主要水質(zhì)影響指標(biāo)為COD、CODMn和BOD5. 結(jié)合評(píng)價(jià)結(jié)果,分析水質(zhì)變化原因,入湖河流水質(zhì)狀況和湖區(qū)底泥是其重要影響因素.

    主成分分析法;歸一化;南四湖上級(jí)湖;水質(zhì)時(shí)間序列;水質(zhì)評(píng)價(jià)

    湖泊水質(zhì)主要受外源污染和內(nèi)源污染的雙重影響,入湖河流和湖區(qū)底泥是兩個(gè)重要的污染來源. 南四湖上級(jí)湖的入湖河流眾多,河流沿途多經(jīng)過城鎮(zhèn),可能攜帶生活、工業(yè)污水;由于南水北調(diào)東線工程完成之前南四湖污染較為嚴(yán)重,含氮、磷等污染物的底泥蓄積深度達(dá)10~15 cm,而這層沉積物與水體間存在頻繁的物質(zhì)交換,也將成為水體的污染源[1]. 對(duì)南四湖的水質(zhì)變化趨勢進(jìn)行深入研究,可以指導(dǎo)人們在社會(huì)與經(jīng)濟(jì)發(fā)展中能夠更加有效合理地開發(fā)、利用湖泊水資源,保護(hù)湖泊水環(huán)境,維持生態(tài)平衡[2]. 為進(jìn)一步有效地改善與提升湖泊水質(zhì),綜合全面分析南四湖的水質(zhì),特別是分析南水北調(diào)工程實(shí)施以來水質(zhì)的變化成因及變化規(guī)律是非常必要和迫切的.

    主成分分析法(PCA)是一種將多維因子納入統(tǒng)一系統(tǒng)中進(jìn)行定量化研究且理論較完善的多元統(tǒng)計(jì)方法[3]. 在此前的研究中,主成分分析法一般用于對(duì)同時(shí)間、不同空間監(jiān)測點(diǎn)的水質(zhì)評(píng)價(jià)進(jìn)行比較,本文嘗試采用該方法對(duì)同一空間監(jiān)測點(diǎn)進(jìn)行時(shí)間上的分析評(píng)價(jià). 通過對(duì)評(píng)價(jià)結(jié)果的分析,尋找南四湖上級(jí)湖水質(zhì)變化趨勢.

    1 評(píng)價(jià)資料與分析方法

    1.1研究區(qū)概況

    南四湖(34°27′~35°20′N,116°34′ ~117°21′E)位于山東省西南部,上級(jí)湖主要包括南陽湖、獨(dú)山湖和昭陽湖,南四湖為我國北方最大的淡水湖,二級(jí)湖壩將南四湖分為上級(jí)湖與下級(jí)湖,上級(jí)湖面積約602 km2,下級(jí)湖644 km2,平均水深1.5m,汛期最大水深2.5~3.0 m;湖東入湖河流為山溪性河流,湖西入湖為平原坡水性河流. 上級(jí)湖來水量占南四湖全湖來水量的88.4%. 南四湖是南水北調(diào)東線工程的調(diào)蓄湖泊,其水質(zhì)的好壞直接關(guān)系到濟(jì)寧、棗莊等地經(jīng)濟(jì)社會(huì)可持續(xù)發(fā)展和水生態(tài)環(huán)境的改善. 國務(wù)院明確指出,南水北調(diào)東線工程順利實(shí)施的前提是節(jié)水,治污是關(guān)鍵,而南四湖就是治污的重點(diǎn)[4]. 了解南四湖的水質(zhì)變化規(guī)律,對(duì)南四湖的整治管理及其周邊地區(qū)經(jīng)濟(jì)社會(huì)的可持續(xù)發(fā)展具有重要意義.

    研究采用南四湖上級(jí)湖的6個(gè)水質(zhì)監(jiān)測點(diǎn),標(biāo)記為S1~S6,分別為王廟、前白口、南陽、獨(dú)山、沙堤和二級(jí)壩閘上(圖1),其中,王廟和沙堤監(jiān)測點(diǎn)分別鄰近洙水河和城郭河的入湖口;洙水河主要接納嘉祥縣和任城區(qū)唐口鎮(zhèn)的工業(yè)和城市污水,是濟(jì)寧市的主要納污河流之一;城郭河是流經(jīng)滕州市最主要的河流.

    圖1 南四湖流域和監(jiān)測點(diǎn)分布Fig.1 The distribution of Lake Nansi and the 6 stations in upper-reach lakes of Lake Nansi

    1.2資料選取

    因已有數(shù)據(jù)資料中,所有年份均有奇數(shù)月監(jiān)測數(shù)據(jù),部分年份有偶數(shù)月的監(jiān)測數(shù)據(jù),為保證數(shù)據(jù)序列時(shí)間間隔的一致性,本研究所采用的數(shù)據(jù)為南陽等6個(gè)監(jiān)測點(diǎn)2008年1月-2014年11月期間奇數(shù)月份共42個(gè)月的監(jiān)測數(shù)據(jù).

    《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB 3838-2002)[5]中共有24項(xiàng)水質(zhì)指標(biāo). 但根據(jù)已有監(jiān)測資料,考慮到監(jiān)測數(shù)據(jù)的合理性以及某些指標(biāo)監(jiān)測值小于檢出限,且pH值都在6~9之間,符合水質(zhì)標(biāo)準(zhǔn)限值,篩選出溶解氧(DO)、化學(xué)需氧量(COD)、高錳酸鹽指數(shù)(CODMn)、五日生化需氧量(BOD5)、氨氮(NH3-N)、總磷(TP)和氟化物7項(xiàng)指標(biāo)組成本次水質(zhì)評(píng)價(jià)的指標(biāo)體系.

    1.3改進(jìn)的歸一化方法

    常用的水質(zhì)指標(biāo)歸一化方法都不同程度地放大或縮小了水質(zhì)指標(biāo)所屬的水質(zhì)類別,特別是當(dāng)越小越優(yōu)型指標(biāo)的水質(zhì)類別屬于劣Ⅴ類,且超出Ⅴ類標(biāo)準(zhǔn)很多時(shí),常用的歸一化方法便不再適用. 為了更好地將水質(zhì)指標(biāo)的監(jiān)測值和水質(zhì)類別標(biāo)準(zhǔn)保持一致,本文采取如下的歸一化方法:

    根據(jù)水質(zhì)評(píng)價(jià)標(biāo)準(zhǔn)將數(shù)值等間距劃分(表1),介于兩類標(biāo)準(zhǔn)之間的采用線性插值計(jì)算. 效益型(越大越好)和成本型(越小越好)的指標(biāo)歸一化計(jì)算公式分別見式(1)和式(2). 本節(jié)中的歸一化方法主要側(cè)重于水質(zhì)指標(biāo)監(jiān)測值無量綱化的處理,因此,歸一化后的數(shù)值并不僅限于[0,1].

    表1 水質(zhì)監(jiān)測值歸一化標(biāo)準(zhǔn)

    (1)

    (2)

    式中,xij為水質(zhì)指標(biāo)的原始監(jiān)測數(shù)據(jù);b1j、b2j、b3j、b4j和b5j分別為水質(zhì)指標(biāo)的Ⅰ、Ⅱ、Ⅲ、Ⅳ和Ⅴ類相對(duì)應(yīng)的水質(zhì)標(biāo)準(zhǔn).

    經(jīng)過上述歸一化后,可以直觀地通過數(shù)值的大小來反映水質(zhì)狀況,數(shù)值越大代表水質(zhì)越好;數(shù)值越小代表水質(zhì)越差. 以王廟監(jiān)測點(diǎn)為例,2008年1月7項(xiàng)指標(biāo)的監(jiān)測值和歸一化值見表2. 由表2可見,各指標(biāo)水質(zhì)好壞程度依次為:DO>氟化物>NH3-N> CODMn>COD>TP>BOD5.

    1.4主成分分析法

    本文在已有研究的基礎(chǔ)上,建立水質(zhì)綜合評(píng)價(jià)指標(biāo)體系,評(píng)價(jià)某一測站隨時(shí)間的水質(zhì)污染程度變化. 分析步驟如下:

    表2 王廟監(jiān)測點(diǎn)2008年1月7項(xiàng)指標(biāo)的監(jiān)測值和歸一化值

    ①將監(jiān)測值矩陣X中的各變量xij進(jìn)行無量綱處理,具體方法見1.3節(jié);②在標(biāo)準(zhǔn)化矩陣的基礎(chǔ)上計(jì)算原始指標(biāo)相關(guān)矩陣R;③計(jì)算相關(guān)矩陣R的特征值、方差貢獻(xiàn)率、累計(jì)方差貢獻(xiàn)率、因子荷載矩陣,按累計(jì)貢獻(xiàn)率達(dá)70%,確定主成分個(gè)數(shù)n;④選取前n個(gè)主成分,計(jì)算相應(yīng)的單位特征向量,即可得主成分計(jì)算公式;⑤將各待評(píng)月份的標(biāo)準(zhǔn)化數(shù)據(jù)分別代入各主成分的表達(dá)式中,計(jì)算得出各月的各主成分得分Fi,以方差貢獻(xiàn)率(di)為權(quán)重求和計(jì)算各月綜合得分[6],水質(zhì)綜合得分序列以下簡稱為得分序列.

    2 評(píng)價(jià)結(jié)果及分析

    2.1綜合得分時(shí)間序列

    對(duì)南四湖上級(jí)湖王廟等6個(gè)監(jiān)測點(diǎn)的水質(zhì)綜合狀況進(jìn)行評(píng)價(jià),各項(xiàng)綜合得分值即是某一監(jiān)測點(diǎn)各月份水質(zhì)污染程度的定量化描述,得分大小可以反映其水質(zhì)好壞,根據(jù)評(píng)價(jià)結(jié)果以及監(jiān)測點(diǎn)地理位置的分布,6個(gè)監(jiān)測點(diǎn)的水質(zhì)綜合得分序列見圖2.

    圖2 各監(jiān)測點(diǎn)綜合得分序列 Fig.2 The score sequence of each station

    2.2水質(zhì)變化趨勢分析

    2.2.16個(gè)監(jiān)測點(diǎn)水質(zhì)變化趨勢王廟和沙堤監(jiān)測點(diǎn)的水質(zhì)在2010年末之后得分變大,出現(xiàn)穩(wěn)定變好的趨勢(圖2). 王廟和沙堤監(jiān)測點(diǎn)分別靠近洙水河和城郭河的入湖口,這2個(gè)監(jiān)測點(diǎn)的水質(zhì)狀況直接受該兩條河流水質(zhì)的影響(圖1). 自2008年以來,洙水河水質(zhì)呈現(xiàn)逐步變好趨勢,以NH3-N和TP為例,其年均濃度分別由9.44和0.64 mg/L減小為1.80和0.20 mg/L;城郭河水質(zhì)指標(biāo)變化幅度小于洙水河,但也明顯變好,其NH3-N和TP年均濃度分別由2.55和0.49 mg/L減小為1.16和0.23 mg/L;DO濃度明顯增大變好(圖3). 因此,王廟和沙堤監(jiān)測點(diǎn)水質(zhì)受入湖河流水質(zhì)影響,穩(wěn)定變好.

    前白口、南陽和獨(dú)山3個(gè)監(jiān)測點(diǎn)水質(zhì)變好趨勢不明顯,且呈現(xiàn)出明顯的季節(jié)差異,水質(zhì)冬季好、夏季差(圖2). 分析單項(xiàng)水質(zhì)指標(biāo),3個(gè)監(jiān)測點(diǎn)的DO和NH3-N濃度均呈現(xiàn)明顯的季節(jié)性,DO濃度冬季高、夏季低,NH3-N濃度冬季低、夏季高,其余指標(biāo)年際變化無明顯規(guī)律,這與整體水質(zhì)狀況是對(duì)應(yīng)的. 因數(shù)據(jù)僅選取奇數(shù)月份,故分析3個(gè)監(jiān)測點(diǎn)2011-2014年5、7、9月各項(xiàng)指標(biāo)變化,尤以DO、NH3-N、BOD5、TP濃度的變差趨勢最為明顯. 分析入湖河流水質(zhì)狀況,以距離3個(gè)監(jiān)測點(diǎn)較近的東魚河、泗河、萬福河為例,3條河流單項(xiàng)水質(zhì)指標(biāo)年均濃度大部分有變差趨勢(圖3),因此,受入湖河流水質(zhì)的影響,前白口、南陽和獨(dú)山3個(gè)監(jiān)測點(diǎn)的水質(zhì)并無明顯變好趨勢. 南四湖流域四季分明,流域內(nèi)平均降水量為620~740 mm左右,年內(nèi)降水分配不均勻,有60%~70%集中于6-9月[7]. 在汛期,入湖河流的流量和流速均明顯變大,沖刷河流沿途城鎮(zhèn)的生活垃圾進(jìn)入湖區(qū)內(nèi),且極易引起湖區(qū)底泥擾動(dòng),加大水體與湖底沉積物的物質(zhì)交換,從而使水體質(zhì)量變差.

    二級(jí)壩閘上監(jiān)測點(diǎn)的水質(zhì)變化趨勢不明顯,2010年底水質(zhì)變好之后,在2013年11月水質(zhì)突然變差(圖2). 分析水質(zhì)單項(xiàng)指標(biāo),2013年9月NH3-N濃度為0.42 mg/L,在2013年11月變?yōu)?.56 mg/L,TP濃度由小于檢出限到0.1mg/L,其余水質(zhì)指標(biāo)濃度變化不大. 當(dāng)湖水流速不大時(shí),南四湖上覆湖水的有機(jī)質(zhì)、氮磷濃度與湖底表層底泥有機(jī)質(zhì)、氮磷濃度的相關(guān)性并不是很大[8],南水北調(diào)東線一期工程在2013年10月試通水,二級(jí)壩閘門打開,使得二級(jí)壩附近的水體流量達(dá)到125 m3/s,基于南四湖南北狹長的地理特性,流速變大,水流的沖擊力導(dǎo)致該處底泥擾動(dòng),底泥中豐富的有機(jī)質(zhì)和氮、磷等污染物進(jìn)入水體中,從而導(dǎo)致二級(jí)壩閘上水質(zhì)明顯變差;且二級(jí)壩閘下極易蓄積污染物,調(diào)水使得這些污染物進(jìn)入上級(jí)湖.

    圖3 南四湖上級(jí)湖入湖河流水質(zhì)指標(biāo)年均濃度變化Fig.3 The indexes’ annual variation of the rivers into upper-reach lakes of Lake Nansi

    表3 南四湖上級(jí)湖2000-2014年水質(zhì)類別

    2.2.2南四湖上級(jí)湖整體水質(zhì)變化趨勢采用全參數(shù)評(píng)價(jià)法[9]計(jì)算南四湖上級(jí)湖2008-2014年的水質(zhì)類別. 2000-2007年水質(zhì)類別查閱山東省水資源公報(bào)(2000-2007年)[10]. 由表3可見,南四湖上級(jí)湖2008-2014年的水質(zhì)較2000-2007年有明顯改善,特別是2011-2014年,水質(zhì)達(dá)到Ⅲ類水標(biāo)準(zhǔn).

    2010年,由于山東省南水北調(diào)截蓄導(dǎo)用工程的實(shí)施,該工程對(duì)南四湖的水質(zhì)狀況產(chǎn)生了明顯的影響. 中水截蓄導(dǎo)用工程的建設(shè),使中水變成了工業(yè)、農(nóng)業(yè)灌溉、城市景觀等用水,從而要求污水處理廠、工業(yè)點(diǎn)源做到達(dá)標(biāo)排放. 以滕州市為例,滕州市城郭河中水截蓄在干線輸水期回用水量為1201.2×104t,通過中水回用可直接減少的COD和NH3-N入河量分別為1133.1和98.7t/a[11],從而使得城郭河入南四湖的COD和NH3-N量減少;在濟(jì)寧市環(huán)南四湖大生態(tài)帶工程建設(shè)規(guī)劃中,2008-2010年為規(guī)劃的第1階段,湖區(qū)修復(fù)和水質(zhì)凈化工程的建設(shè)實(shí)施[12],保證了南四湖上級(jí)湖水質(zhì)的穩(wěn)定.

    2.3水體主要污染指標(biāo)分析

    在利用主成分分析法進(jìn)行水質(zhì)評(píng)價(jià)的過程中,由主成分荷載矩陣可以得出王廟等6個(gè)監(jiān)測點(diǎn)共有的主要水質(zhì)污染指標(biāo)為COD、CODMn和BOD5,且在各月份最差指標(biāo)所占比例的統(tǒng)計(jì)中,該3項(xiàng)指標(biāo)在42個(gè)奇數(shù)月中所占比重均在10%以上,特別是BOD5比重均在30%以上. 因此,南四湖上級(jí)湖主要是有機(jī)物污染.

    南四湖整體水質(zhì)變差,受近年來沿湖地區(qū)無序開發(fā),以及強(qiáng)度過大影響較為顯著,文獻(xiàn)[13]中詳細(xì)分析了南四湖入湖河流的水質(zhì)及主要污染指標(biāo),即南陽湖上級(jí)湖周邊入湖河流的有機(jī)物污染嚴(yán)重,造成了接納入湖河流的湖區(qū)有機(jī)物污染嚴(yán)重,本文分析得出的主要污染指標(biāo)與其一致. 所以河流是南四湖上級(jí)湖的主要污染源,湖區(qū)水環(huán)境的恢復(fù)治理首先要解決入湖河流的水污染問題.

    3 結(jié)論及建議

    1)為了更好地使水質(zhì)指標(biāo)的監(jiān)測值和水質(zhì)類別標(biāo)準(zhǔn)保持一致,本文對(duì)常用的歸一化方法進(jìn)行了改進(jìn),該方法客觀地反映各類指標(biāo)所屬的水質(zhì)類別. 采用主成分分析法進(jìn)行同一監(jiān)測點(diǎn)水質(zhì)在時(shí)間上的評(píng)價(jià),得出的評(píng)價(jià)結(jié)果符合物理機(jī)制,能夠較好地考慮各項(xiàng)指標(biāo)之間的相關(guān)性,可用于水質(zhì)時(shí)間變化分析.

    2)6個(gè)監(jiān)測點(diǎn)在湖區(qū)分布位置不同,水質(zhì)變化趨勢也存在較大差異,水質(zhì)主要影響因素也不同. 王廟和沙堤監(jiān)測點(diǎn)水質(zhì)在2010年末呈現(xiàn)穩(wěn)定變好趨勢;前白口、南陽和獨(dú)山監(jiān)測點(diǎn)水質(zhì)并未變好,但呈現(xiàn)出明顯的季節(jié)性;二級(jí)壩閘上水質(zhì)變化趨勢不明顯,但在2013年11月水質(zhì)明顯變差.

    3)南四湖上級(jí)湖水質(zhì)較2010年之前有明顯改善,南水北調(diào)山東省截蓄導(dǎo)用工程的實(shí)施對(duì)水質(zhì)改善起到了重要作用;但在南四湖調(diào)水期間,底泥擾動(dòng)等因素會(huì)影響水體質(zhì)量,應(yīng)在此期間嚴(yán)格控制水體流速,改善水力特性,減少底泥及外源對(duì)水體的雙重污染,從而保證輸水通道的水質(zhì)安全.

    4)南四湖上級(jí)湖湖區(qū)水質(zhì)的主要影響指標(biāo)為COD、CODMn和BOD5,其有機(jī)污染的特點(diǎn)較明顯,由于缺少入湖河流水質(zhì)監(jiān)測數(shù)據(jù),本研究僅根據(jù)沿湖地區(qū)經(jīng)濟(jì)社會(huì)發(fā)展現(xiàn)狀,推斷入湖河流的水質(zhì)狀況對(duì)湖區(qū)及出湖水質(zhì)的影響應(yīng)該是起決定性作用的因素. 南四湖的治理需要針對(duì)不同水域的水質(zhì)差別,查明污染源,實(shí)施入湖河流流域污染的源頭控制與湖泊內(nèi)源污染防治相結(jié)合,以及工程治污、結(jié)構(gòu)減排與監(jiān)管控污相結(jié)合.

    [1]辛穎, 武周虎, 慕金波等. 南四湖入湖、入干線河流與輸水干線的水質(zhì)動(dòng)態(tài)分析. 安全與環(huán)境工程, 2013, 20(3): 36-41.

    [2]宮瑩. 城市中小型淺水湖泊水環(huán)境預(yù)測方法應(yīng)用研究[學(xué)位論文]. 南京: 河海大學(xué), 2004: 10-11.

    [3]劉小楠, 崔巍. 主成分分析法在汾河水質(zhì)評(píng)價(jià)中的應(yīng)用. 中國給水排水, 2009, 25(18): 105-108.

    [4]倪曉. 南四湖流域水污染總量控制及水質(zhì)綜合改善方案研究[學(xué)位論文]. 濟(jì)南: 山東師范大學(xué), 2013: 19-22.

    [5]地表水環(huán)境質(zhì)量標(biāo)準(zhǔn). GB 3838-2002.

    [6]楊競, 童禎恭, 劉玉哲. SPSS軟件對(duì)飲用水水質(zhì)進(jìn)行主成分分析評(píng)價(jià)的運(yùn)用. 環(huán)境科學(xué)與技術(shù), 2011, 34(7): 171-174.

    [7]趙群群, 楊凱. 南四湖污染物排放對(duì)南水北調(diào)東線水質(zhì)的影響及治理措施. 價(jià)值工程, 2010, (8): 105-106.

    [8]張雙雙, 武周虎, 張潔等. 南四湖表層底泥有機(jī)質(zhì)及氮磷分布特征. 青島理工大學(xué)學(xué)報(bào), 2014, 35(5): 64-69.

    [9]山東省水功能區(qū)水質(zhì)通報(bào), 2013.

    [10]山東省水資源公報(bào), 2000-2007.

    [11]鄭從奇. 南水北調(diào)山東段中水截蓄導(dǎo)用工程綜合效益評(píng)價(jià)研究[學(xué)位論文]. 濟(jì)南: 山東大學(xué), 2013: 22.

    [12]南水北調(diào)東線濟(jì)寧市環(huán)南四湖大生態(tài)帶工程建設(shè)規(guī)劃. 濟(jì)寧市人民政府, 2009: 15.

    [13]李文猛. 南四湖及其入湖河流的水質(zhì)評(píng)價(jià)分析[學(xué)位論文]. 天津: 天津大學(xué), 2012: 50-56.

    [14]趙世新, 張晨, 高學(xué)平等. 南水北調(diào)東線調(diào)度對(duì)南四湖水質(zhì)的影響. 湖泊科學(xué), 2012, 24(6): 923-931. DOI 10. 18307/2012.0616.

    [15]李靖華, 郭耀煌. 主成分分析用于多指標(biāo)評(píng)價(jià)的方法研究——主成分評(píng)價(jià). 管理工程學(xué)報(bào), 2002, 16(1): 39-43.

    [16]鄒海明, 蔣良富, 李粉茹. 基于主成分分析的水質(zhì)評(píng)價(jià)方法. 數(shù)學(xué)的實(shí)踐與認(rèn)識(shí), 2008, 38(8): 85-90.

    [17]劉臣輝, 呂信紅, 范海燕. 主成分分析法用于環(huán)境質(zhì)量評(píng)價(jià)的探討. 環(huán)境科學(xué)與管理, 2011, 26(3): 183-186.

    [18]劉遠(yuǎn)書, 杜婷, 黃躍飛等. 南四湖截污導(dǎo)流工程水量水質(zhì)聯(lián)合調(diào)度運(yùn)行研究. 河海大學(xué)學(xué)報(bào): 自然工學(xué)版, 2013, 41(5): 395-399.

    [19]高學(xué)平, 李文猛, 張晨等. 入湖河流對(duì)南四湖水質(zhì)的影響. 水資源保護(hù), 2013, 29(2): 1-5.

    Evaluation of temporal water quality change (2008-2014) and the cause analysis in the upper-reach lakes of Lake Nansi, Shandong Province

    XIE Wenjing1, WANG Song2, LOU Shanchong2, CAO Shengle1, SUN Xiuling1, WANG Yuemin1& LIN Jie1

    (1:SchoolofCivilEngineering,ShandongUniversity,Jinan250061,P.R.China)(2:HydrologyandWaterResourcesBureauofShandongProvince,Jinan250002,P.R.China)

    Based on water quality information from 6 stations in the upper-reach lakes of Lake Nansi which was sampled from a total of 42 odd months during 2008 and 2014, choosing dissolved oxygen, COD, CODMn, BOD5, ammonia nitrogen, total phosphorus and fluoride as the indexes of water quality evaluation. The water quality index normalization was used to improve the method, and the principal component analysis was used to evaluate water quality comprehensively for the data from the six monitoring stations. According to the obtained comprehensive score of the monthly sequence, the trend of water quality status was analyzed. Compared with the general condition of water quality during 2000 and 2007, it is found that water quality was significantly better during 2008 and 2014. The evaluation results showed that the water quality at any monitoring station had obvious changes at different time, and had different changes at any time sequence on different water quality monitoring stations. Besides, the main influenced indicators of water quality in the upper-reach lakes of Lake Nansi were COD, CODMnand BOD5. According to the evaluation results and the analyzed changes in the water quality, we found that the river quality into the lake and lake sediments were the two important factors influencing on the water quality.

    Principal component analysis; normalized; upper-reach lakes of Lake Nansi; time series of water quality; water quality evaluation

    J.LakeSci.(湖泊科學(xué)), 2016, 28(3): 513-519

    10.18307/2016.0306

    ?2016 byJournalofLakeSciences

    *山東省省級(jí)水利科研及技術(shù)推廣項(xiàng)目(SDSLSK201302)資助.2015-06-09收稿;2015-08-19收修改稿.解文靜(1990~),女,碩士研究生;E-mail: 994306933@qq.com.

    猜你喜歡
    南四湖入湖監(jiān)測點(diǎn)
    降雨事件下洱海北部主要入湖河流污染物特征分析*
    天津南港LNG接收站沉降監(jiān)測點(diǎn)位布設(shè)
    煤氣與熱力(2022年4期)2022-05-23 12:44:56
    環(huán)洱海主要入湖河流水質(zhì)特征及入湖污染負(fù)荷估算
    人民長江(2022年1期)2022-04-01 02:44:35
    撫河流域綜合治理監(jiān)測布局優(yōu)化
    不忘初心創(chuàng)沂沭泗基建新篇 除險(xiǎn)加固保南四湖治水安瀾
    治淮(2021年12期)2021-12-31 05:46:30
    全站儀極坐標(biāo)法監(jiān)測點(diǎn)穩(wěn)定性分析方法研究
    入湖河口濕地恢復(fù)與重建規(guī)劃設(shè)計(jì)初探——以資興市興寧河入湖河口濕地為例
    滇池入湖河流磷負(fù)荷時(shí)空變化及形態(tài)組成貢獻(xiàn)*
    我省舉辦家畜血吸蟲病監(jiān)測點(diǎn)培訓(xùn)班
    沂沭泗河洪水東調(diào)南下續(xù)建工程南四湖湖西大堤加固工程(山東實(shí)施段)暨南四湖二級(jí)壩閘下西段股河上段剩余工程通過竣工驗(yàn)收
    中國水利(2015年13期)2015-02-28 15:14:14
    欧美性感艳星| 成人漫画全彩无遮挡| 午夜亚洲福利在线播放| 日本爱情动作片www.在线观看 | 欧美成人免费av一区二区三区| 尤物成人国产欧美一区二区三区| 男女那种视频在线观看| 亚洲三级黄色毛片| 国产伦在线观看视频一区| 欧美不卡视频在线免费观看| 久久久久性生活片| 久久午夜亚洲精品久久| 午夜福利成人在线免费观看| 国产日本99.免费观看| 99热这里只有是精品50| 精品久久久久久久久av| 亚洲高清免费不卡视频| 亚洲人与动物交配视频| 给我免费播放毛片高清在线观看| 麻豆成人午夜福利视频| 久久久久性生活片| 晚上一个人看的免费电影| 精品99又大又爽又粗少妇毛片| 亚洲在线观看片| 国产色爽女视频免费观看| 久久久久久久久中文| 在线天堂最新版资源| 亚洲精品日韩在线中文字幕 | 一级a爱片免费观看的视频| 国产精品国产三级国产av玫瑰| 国产成人福利小说| 国模一区二区三区四区视频| 99热6这里只有精品| 久久精品国产亚洲网站| 亚洲av电影不卡..在线观看| 久久人妻av系列| 看片在线看免费视频| 久久久久久久久中文| 99久久九九国产精品国产免费| 不卡一级毛片| 亚洲成人精品中文字幕电影| www.色视频.com| 一级毛片电影观看 | 亚洲国产精品合色在线| av中文乱码字幕在线| 99热只有精品国产| 97在线视频观看| 夜夜爽天天搞| 啦啦啦啦在线视频资源| 国产探花极品一区二区| 国产伦精品一区二区三区视频9| 欧美在线一区亚洲| 国产又黄又爽又无遮挡在线| 12—13女人毛片做爰片一| 国产成年人精品一区二区| 午夜免费激情av| 一本精品99久久精品77| 人妻夜夜爽99麻豆av| 久久精品影院6| 亚洲国产色片| 在线播放无遮挡| av天堂在线播放| 综合色av麻豆| 亚洲av不卡在线观看| 亚洲欧美日韩东京热| 91av网一区二区| 狠狠狠狠99中文字幕| 中出人妻视频一区二区| 成人高潮视频无遮挡免费网站| 精华霜和精华液先用哪个| 99九九线精品视频在线观看视频| 成人亚洲精品av一区二区| 99久久九九国产精品国产免费| 精品国产三级普通话版| 老司机影院成人| 精品国内亚洲2022精品成人| 国产av在哪里看| 国产精品女同一区二区软件| 人妻久久中文字幕网| 看免费成人av毛片| 久久精品国产清高在天天线| 春色校园在线视频观看| 成人综合一区亚洲| 蜜臀久久99精品久久宅男| 听说在线观看完整版免费高清| 国产视频内射| 婷婷精品国产亚洲av| av免费在线看不卡| 午夜爱爱视频在线播放| 特级一级黄色大片| 香蕉av资源在线| 插阴视频在线观看视频| 99久久精品一区二区三区| 亚洲精品亚洲一区二区| 午夜亚洲福利在线播放| 丝袜喷水一区| 亚洲精品日韩在线中文字幕 | 97超级碰碰碰精品色视频在线观看| 日日摸夜夜添夜夜添av毛片| 免费看a级黄色片| 国产精品一区二区三区四区久久| 亚洲人成网站在线播放欧美日韩| 两性午夜刺激爽爽歪歪视频在线观看| 国产黄片美女视频| 永久网站在线| 久久精品久久久久久噜噜老黄 | 国产探花在线观看一区二区| 伊人久久精品亚洲午夜| 国产精品一区二区性色av| 欧美日韩精品成人综合77777| 亚洲熟妇熟女久久| 美女内射精品一级片tv| 亚洲人与动物交配视频| 成人美女网站在线观看视频| 婷婷精品国产亚洲av| 亚洲国产精品sss在线观看| 黄色视频,在线免费观看| 精品一区二区三区视频在线| 亚洲无线观看免费| 久久久成人免费电影| 男插女下体视频免费在线播放| 青春草视频在线免费观看| 精品人妻偷拍中文字幕| 卡戴珊不雅视频在线播放| 亚洲天堂国产精品一区在线| 此物有八面人人有两片| 免费搜索国产男女视频| 啦啦啦韩国在线观看视频| 97碰自拍视频| 综合色丁香网| 十八禁国产超污无遮挡网站| 人人妻人人澡欧美一区二区| 99riav亚洲国产免费| 久久久精品大字幕| 欧美日韩在线观看h| 午夜亚洲福利在线播放| 亚洲第一电影网av| 日日啪夜夜撸| 久久午夜福利片| 精品一区二区三区av网在线观看| 欧美日本亚洲视频在线播放| 俺也久久电影网| 国产精品电影一区二区三区| www日本黄色视频网| 欧美成人一区二区免费高清观看| 中文字幕av在线有码专区| 国产 一区精品| 亚洲精品影视一区二区三区av| 99热这里只有是精品50| 亚洲精品成人久久久久久| 18禁在线无遮挡免费观看视频 | 91久久精品国产一区二区三区| 欧美丝袜亚洲另类| 在线观看免费视频日本深夜| 狂野欧美激情性xxxx在线观看| 国产单亲对白刺激| 国产成人精品久久久久久| 欧美性猛交黑人性爽| 秋霞在线观看毛片| 亚洲色图av天堂| 亚洲国产欧洲综合997久久,| 中出人妻视频一区二区| 久久久久久九九精品二区国产| 国产黄色小视频在线观看| 日韩欧美国产在线观看| 色在线成人网| 国产淫片久久久久久久久| 在线国产一区二区在线| 免费电影在线观看免费观看| 色视频www国产| 国语自产精品视频在线第100页| 亚洲精品色激情综合| 美女 人体艺术 gogo| 国内精品一区二区在线观看| 亚洲成人av在线免费| 女生性感内裤真人,穿戴方法视频| 欧美色欧美亚洲另类二区| 看免费成人av毛片| 国产一区二区亚洲精品在线观看| 国产精品伦人一区二区| 国产女主播在线喷水免费视频网站 | 久久国内精品自在自线图片| 国产精品乱码一区二三区的特点| 有码 亚洲区| 又爽又黄无遮挡网站| 国产精品亚洲一级av第二区| 免费看美女性在线毛片视频| 大又大粗又爽又黄少妇毛片口| 自拍偷自拍亚洲精品老妇| 美女内射精品一级片tv| 国产精品亚洲美女久久久| 丰满的人妻完整版| 久久久国产成人精品二区| 小说图片视频综合网站| 最近中文字幕高清免费大全6| 国产高清三级在线| 久久精品久久久久久噜噜老黄 | 99riav亚洲国产免费| 午夜日韩欧美国产| 国产男人的电影天堂91| 亚洲av电影不卡..在线观看| av天堂在线播放| 一级毛片我不卡| 免费av观看视频| 国产伦精品一区二区三区视频9| 国产黄色视频一区二区在线观看 | 91在线观看av| 在线免费观看的www视频| 啦啦啦啦在线视频资源| 午夜福利高清视频| 美女cb高潮喷水在线观看| 久久久久免费精品人妻一区二区| 国产男人的电影天堂91| 校园春色视频在线观看| 人人妻人人澡欧美一区二区| 成人性生交大片免费视频hd| 九九在线视频观看精品| 亚洲av熟女| 日韩一区二区视频免费看| 3wmmmm亚洲av在线观看| 最近2019中文字幕mv第一页| 成年女人看的毛片在线观看| 日本爱情动作片www.在线观看 | 男女之事视频高清在线观看| a级一级毛片免费在线观看| avwww免费| 成人性生交大片免费视频hd| 亚洲专区国产一区二区| 一个人免费在线观看电影| 国产在线精品亚洲第一网站| 日本 av在线| 热99re8久久精品国产| 亚洲专区国产一区二区| 亚洲av熟女| 麻豆一二三区av精品| 亚洲美女视频黄频| 一级av片app| 不卡一级毛片| 一边摸一边抽搐一进一小说| 99久久九九国产精品国产免费| 色综合站精品国产| 99热这里只有是精品在线观看| 精品少妇黑人巨大在线播放 | 国产日本99.免费观看| av专区在线播放| 91午夜精品亚洲一区二区三区| 久久久久久国产a免费观看| 久99久视频精品免费| 啦啦啦观看免费观看视频高清| 国产爱豆传媒在线观看| 亚洲18禁久久av| 国产片特级美女逼逼视频| 日韩大尺度精品在线看网址| 色哟哟哟哟哟哟| 国国产精品蜜臀av免费| 国产成人aa在线观看| 中国美白少妇内射xxxbb| 自拍偷自拍亚洲精品老妇| 午夜亚洲福利在线播放| 天堂影院成人在线观看| 国产女主播在线喷水免费视频网站 | 成人欧美大片| 亚洲自偷自拍三级| 中文在线观看免费www的网站| 国产真实伦视频高清在线观看| 亚洲欧美中文字幕日韩二区| 寂寞人妻少妇视频99o| 99久久无色码亚洲精品果冻| 人妻制服诱惑在线中文字幕| 偷拍熟女少妇极品色| 免费观看精品视频网站| 丰满的人妻完整版| 99久国产av精品| 一进一出抽搐gif免费好疼| 18禁在线无遮挡免费观看视频 | 亚洲欧美清纯卡通| 老熟妇乱子伦视频在线观看| 日日干狠狠操夜夜爽| 在线国产一区二区在线| 一级毛片电影观看 | 男人和女人高潮做爰伦理| 网址你懂的国产日韩在线| 听说在线观看完整版免费高清| 成人综合一区亚洲| 亚洲18禁久久av| 亚洲第一电影网av| 草草在线视频免费看| 亚洲一区高清亚洲精品| av中文乱码字幕在线| 亚洲欧美中文字幕日韩二区| 全区人妻精品视频| 久久久久国内视频| 亚洲专区国产一区二区| 亚洲婷婷狠狠爱综合网| 精品福利观看| 99国产精品一区二区蜜桃av| 成人高潮视频无遮挡免费网站| 噜噜噜噜噜久久久久久91| 亚洲三级黄色毛片| 久久人人爽人人爽人人片va| 99热这里只有是精品50| 亚洲av美国av| 最近最新中文字幕大全电影3| 国内揄拍国产精品人妻在线| 免费黄网站久久成人精品| 2021天堂中文幕一二区在线观| 亚洲国产精品成人久久小说 | 亚洲综合色惰| 国内久久婷婷六月综合欲色啪| 精华霜和精华液先用哪个| 欧美色欧美亚洲另类二区| 欧美一区二区精品小视频在线| 免费观看的影片在线观看| 国产黄片美女视频| 亚洲精品影视一区二区三区av| 亚洲欧美日韩东京热| 欧美性感艳星| 欧美人与善性xxx| 久久久精品欧美日韩精品| 看片在线看免费视频| 久久人人爽人人爽人人片va| 99热精品在线国产| 欧美精品国产亚洲| 成人综合一区亚洲| 精品久久久久久久久av| 日韩大尺度精品在线看网址| 在线观看av片永久免费下载| 熟女人妻精品中文字幕| 午夜福利高清视频| 最近的中文字幕免费完整| 91久久精品电影网| 久久久国产成人免费| 中文亚洲av片在线观看爽| 久久韩国三级中文字幕| 波野结衣二区三区在线| videossex国产| 三级经典国产精品| 97热精品久久久久久| 网址你懂的国产日韩在线| 岛国在线免费视频观看| 久久久久国产精品人妻aⅴ院| 一个人免费在线观看电影| 久久久久免费精品人妻一区二区| 亚洲av一区综合| 欧美高清成人免费视频www| 有码 亚洲区| 99精品在免费线老司机午夜| 国国产精品蜜臀av免费| 午夜日韩欧美国产| 18禁在线播放成人免费| 不卡视频在线观看欧美| 91在线精品国自产拍蜜月| 成人永久免费在线观看视频| 天美传媒精品一区二区| 国产成人freesex在线 | 精品一区二区三区人妻视频| 国产真实乱freesex| 国产成人精品久久久久久| 国产亚洲精品久久久com| 熟女电影av网| 在现免费观看毛片| 欧美一区二区精品小视频在线| 精品一区二区三区视频在线| 大型黄色视频在线免费观看| 欧美色欧美亚洲另类二区| av天堂在线播放| aaaaa片日本免费| 久久欧美精品欧美久久欧美| 国产亚洲欧美98| 蜜臀久久99精品久久宅男| 麻豆乱淫一区二区| 成人三级黄色视频| 免费黄网站久久成人精品| 搡女人真爽免费视频火全软件 | 亚洲精品亚洲一区二区| 久久99热6这里只有精品| 久久久午夜欧美精品| 午夜激情欧美在线| 日本与韩国留学比较| 成人鲁丝片一二三区免费| 日本与韩国留学比较| 欧美3d第一页| 国产一区二区在线观看日韩| 直男gayav资源| 久久久久久久久大av| 最近最新中文字幕大全电影3| 日韩欧美三级三区| 国产av麻豆久久久久久久| 床上黄色一级片| 麻豆一二三区av精品| 国产精品爽爽va在线观看网站| 日韩欧美在线乱码| 搡老妇女老女人老熟妇| 国产白丝娇喘喷水9色精品| 亚洲五月天丁香| av天堂中文字幕网| 国产亚洲精品久久久com| 两个人的视频大全免费| 少妇人妻精品综合一区二区 | 国产伦精品一区二区三区视频9| 欧美国产日韩亚洲一区| 久久久久久伊人网av| 真实男女啪啪啪动态图| 波多野结衣高清作品| 精品国内亚洲2022精品成人| 最近视频中文字幕2019在线8| 色5月婷婷丁香| 我要搜黄色片| 成年版毛片免费区| 免费观看的影片在线观看| 我的老师免费观看完整版| 国产精品一区www在线观看| 日韩精品青青久久久久久| 国产亚洲欧美98| 97人妻精品一区二区三区麻豆| 中出人妻视频一区二区| 国产成人91sexporn| 真实男女啪啪啪动态图| 国产一区二区亚洲精品在线观看| 国产伦在线观看视频一区| 成人午夜高清在线视频| 日本成人三级电影网站| 午夜精品在线福利| 亚洲精品在线观看二区| 亚洲av电影不卡..在线观看| 国产精品,欧美在线| 淫妇啪啪啪对白视频| 久久人人爽人人爽人人片va| 国产亚洲精品av在线| 69人妻影院| 国产精品1区2区在线观看.| 男插女下体视频免费在线播放| 两个人的视频大全免费| 成年免费大片在线观看| 少妇丰满av| 国产综合懂色| 久久久精品94久久精品| 国产成人freesex在线 | 亚洲精品在线观看二区| 一卡2卡三卡四卡精品乱码亚洲| 99久久精品热视频| 啦啦啦韩国在线观看视频| 久久久精品欧美日韩精品| 不卡视频在线观看欧美| 中文亚洲av片在线观看爽| 欧美性感艳星| 人人妻人人澡人人爽人人夜夜 | 亚洲精品一卡2卡三卡4卡5卡| 免费观看的影片在线观看| 18禁黄网站禁片免费观看直播| 美女cb高潮喷水在线观看| 成人av一区二区三区在线看| 色视频www国产| 一区福利在线观看| 日本a在线网址| 免费看av在线观看网站| 色噜噜av男人的天堂激情| 午夜激情欧美在线| 国产不卡一卡二| 亚洲不卡免费看| 97超视频在线观看视频| 欧美高清成人免费视频www| 日本a在线网址| 国产老妇女一区| 最近在线观看免费完整版| av天堂在线播放| 99热这里只有精品一区| 亚洲综合色惰| 赤兔流量卡办理| 亚洲av第一区精品v没综合| 波野结衣二区三区在线| 国产人妻一区二区三区在| 亚洲精品一区av在线观看| 亚洲人成网站在线播| 亚洲一级一片aⅴ在线观看| 麻豆一二三区av精品| а√天堂www在线а√下载| 亚洲专区国产一区二区| 亚洲人与动物交配视频| 少妇丰满av| 成人美女网站在线观看视频| 欧美一区二区亚洲| 日本成人三级电影网站| 最近的中文字幕免费完整| av国产免费在线观看| 精品人妻一区二区三区麻豆 | 精品日产1卡2卡| 色吧在线观看| 欧美日本视频| 嫩草影院入口| 久久久久久九九精品二区国产| 免费看日本二区| 美女黄网站色视频| 日日摸夜夜添夜夜添av毛片| 中国美女看黄片| 免费在线观看影片大全网站| 性欧美人与动物交配| 久久久成人免费电影| 中文字幕久久专区| 久久久精品94久久精品| 国产黄片美女视频| av在线天堂中文字幕| av免费在线看不卡| 最近手机中文字幕大全| 夜夜看夜夜爽夜夜摸| 久久午夜福利片| 亚洲人成网站在线播放欧美日韩| 狂野欧美激情性xxxx在线观看| 尾随美女入室| 亚洲国产精品成人综合色| 又黄又爽又免费观看的视频| 久久九九热精品免费| 久久欧美精品欧美久久欧美| 成人二区视频| 色吧在线观看| 麻豆一二三区av精品| 免费观看人在逋| 国产高清不卡午夜福利| 久久韩国三级中文字幕| 欧美日韩国产亚洲二区| 精品人妻偷拍中文字幕| 人人妻,人人澡人人爽秒播| 91在线精品国自产拍蜜月| 国产伦精品一区二区三区四那| 99久久精品热视频| 精品午夜福利视频在线观看一区| 午夜福利在线在线| 高清毛片免费观看视频网站| 亚洲欧美成人精品一区二区| 亚洲成人中文字幕在线播放| 国产精品无大码| 我要搜黄色片| 99精品在免费线老司机午夜| 黄片wwwwww| 夜夜夜夜夜久久久久| 免费观看在线日韩| 日韩强制内射视频| 日韩欧美国产在线观看| 亚洲精品日韩在线中文字幕 | 欧美日韩在线观看h| 18禁在线无遮挡免费观看视频 | 精品熟女少妇av免费看| 亚洲av第一区精品v没综合| 亚洲人成网站高清观看| 一a级毛片在线观看| 熟女人妻精品中文字幕| 日日摸夜夜添夜夜添av毛片| 变态另类丝袜制服| 国产在线精品亚洲第一网站| 少妇丰满av| 国产精品久久久久久久久免| 国产三级中文精品| 好男人在线观看高清免费视频| 波野结衣二区三区在线| av在线天堂中文字幕| 亚洲中文日韩欧美视频| 亚洲av成人av| 久久精品国产亚洲av天美| 神马国产精品三级电影在线观看| 国产成人a∨麻豆精品| 在线播放无遮挡| 日本精品一区二区三区蜜桃| 极品教师在线视频| 1000部很黄的大片| 日韩成人av中文字幕在线观看 | 欧美+日韩+精品| 有码 亚洲区| 日韩大尺度精品在线看网址| 免费不卡的大黄色大毛片视频在线观看 | 村上凉子中文字幕在线| 久久久久性生活片| 久久精品人妻少妇| 在线免费十八禁| 亚洲中文字幕一区二区三区有码在线看| 欧美成人a在线观看| 欧美不卡视频在线免费观看| 舔av片在线| 中文字幕人妻熟人妻熟丝袜美| 国产三级中文精品| 午夜福利视频1000在线观看| av在线观看视频网站免费| 观看免费一级毛片| 综合色丁香网| 国产精品一区二区性色av| 晚上一个人看的免费电影| 国产黄a三级三级三级人| 三级经典国产精品| 美女 人体艺术 gogo| 久久精品国产清高在天天线| av在线亚洲专区| 亚洲欧美日韩高清专用| 少妇丰满av| 97超碰精品成人国产| 插阴视频在线观看视频| 无遮挡黄片免费观看| 国产色婷婷99| 嫩草影院入口| 最好的美女福利视频网| 中文字幕久久专区| 免费人成视频x8x8入口观看| 国产三级中文精品| 十八禁国产超污无遮挡网站| 狠狠狠狠99中文字幕| 欧美一区二区精品小视频在线| 亚洲国产精品久久男人天堂| or卡值多少钱| 国产午夜福利久久久久久| 午夜福利成人在线免费观看| 国产精品日韩av在线免费观看| 亚洲成人av在线免费| 性欧美人与动物交配| 国产国拍精品亚洲av在线观看| 亚洲经典国产精华液单| aaaaa片日本免费| 亚洲18禁久久av|