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

    天津滯緩流型城市河網(wǎng)水質(zhì)時空分布特征

    2017-07-05 11:38:29董立新白昊陽
    水利水電科技進展 2017年4期
    關鍵詞:河網(wǎng)外環(huán)氮磷

    董立新,白昊陽

    (天津市水利科學研究院,天津 300061)

    ?

    天津滯緩流型城市河網(wǎng)水質(zhì)時空分布特征

    董立新,白昊陽

    (天津市水利科學研究院,天津 300061)

    針對天津城市河網(wǎng)滯緩流特點突出、水質(zhì)變化與天然河流不同的特征,以海河、津河、衛(wèi)津河、外環(huán)河4條重點河道為研究對象,應用聚類分析、因子分析、多元線性回歸分析等多元統(tǒng)計分析方法對11個重點河道斷面2008年和2009年39周的水質(zhì)監(jiān)測指標進行了分析。結果表明:河道斷面水質(zhì)在時間和空間上具有一定的分布特征,時間上,各斷面5—8月水質(zhì)較其他時間更差,氨氮、總磷超標嚴重;空間上,呈現(xiàn)出從海河到津河、衛(wèi)津河再到外環(huán)河水質(zhì)逐步下降的特征;海河、津河、衛(wèi)津河等河道氮磷類匯入是主要污染源,外環(huán)河則以有機物污染匯入為主要污染源。

    城市河網(wǎng);滯緩流;水質(zhì)變化;多元統(tǒng)計分析;天津市

    滯緩流型城市河網(wǎng)是指水體整體呈現(xiàn)出流速滯緩、流動性差等特征的城市河網(wǎng)區(qū)域。在我國北方地區(qū),包括天津在內(nèi)的多個城市河網(wǎng)滯緩流特點均比較突出。由于上游的水資源量較少,在滿足生活、生產(chǎn)用水后才能進行城市河網(wǎng)補給,生態(tài)補水無法得到保障;此外,北方地區(qū)的天然來水多集中在汛期,為保障城市排瀝安全,也無法給河網(wǎng)充分補水。由此,為實現(xiàn)城市河網(wǎng)的景觀環(huán)境功能、存蓄河網(wǎng)水量,常通過建立閘壩等水工設施將城市河網(wǎng)與天然河道分割開來,進行人工調(diào)控。因此,城市河網(wǎng)無法像天然河流一樣自然流動,呈現(xiàn)人工調(diào)控下的滯緩流特征[1],造成了水體流動性差,抑制了河道復氧過程,溶解氧含量長期處于較低水平,水體自凈能力下降,再加上污染匯入、蒸發(fā)滲漏的長期作用,河網(wǎng)水體水質(zhì)總體上較差,多數(shù)處于劣Ⅴ類狀態(tài)[2]。

    河道水質(zhì)變化規(guī)律研究已有多年歷史,目前主要集中在水質(zhì)時空分布特征分析與水質(zhì)綜合評價上[3-5],現(xiàn)階段主要采用的方法是多元統(tǒng)計分析方法,其中,系統(tǒng)聚類分析法、因子分析法、多元線性回歸分析法既可用于河流水質(zhì)的時空變化分析,還可兼顧河流水質(zhì)綜合評價,應用較為廣泛[6-8]。國外Berkant等[9]、André等[10]、Bahman等[11]采用多元統(tǒng)計分析方法分別對土耳其、加拿大、日本的河流開展了水質(zhì)時空變化特征分析,得到了河道主要水質(zhì)指標的變化規(guī)律及其影響因子;國內(nèi)卜紅梅等[7]、張棋等[12]、徐華山等[13]采用系統(tǒng)聚類分析、因子分析和主成分多元線性回歸分析等多元統(tǒng)計分析方法分別對金水河流域、湘江流域、漳衛(wèi)南運河流域水質(zhì)時空變化特征進行了研究,分析了主要污染因子和可能的污染來源。

    綜合來看,目前,多元統(tǒng)計分析方法在天然河流上應用較多,還未見在城市滯緩流河道中的應用。關于北方城市滯緩流河道各方面的研究較少,且偏重于水資源管理、水生態(tài)恢復、污染治理等方面[14],而針對其水質(zhì)時空演變及污染源定性識別研究鮮見報道。本文以天津城市河網(wǎng)中重點河道為研究對象,采用系統(tǒng)聚類分析法分析河道水質(zhì)時空分布特征,采用因子分析法對河流主要污染源進行定性識別,采用多元線性回歸分析法確定主要公因子對各水質(zhì)指標的貢獻率,為北方滯緩流型城市河網(wǎng)水質(zhì)改善和修復提供科學依據(jù)。

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

    1.1 研究區(qū)概況

    天津城市河網(wǎng)目前共有行洪排瀝河道20條,再生水河道4條,河道總長度231.4 km,滯緩流特點突出。選取水環(huán)境影響較大、關注度較高的海河、津河、衛(wèi)津河、外環(huán)河4條重點河道作為代表,在這4條研究河道中選擇了11個重點斷面進行分析,包括海河的三岔口(S1)和四新橋(S2),津河的三元村(S3)、八里臺(S4)、廣東路橋(S5)和津河廣場(S6),衛(wèi)津河的紀莊子橋(S7)和天塔(S8),外環(huán)河的津淶公路橋(S9)、七號橋(S10)和津漢橋(S11)。河道及其重點斷面位置如圖1所示。

    水質(zhì)數(shù)據(jù)由現(xiàn)場采樣監(jiān)測獲得,采樣時間為39周(2008年5月下旬至2008年11月中旬以及2009年3月初至2009年5月下旬),采樣頻次為每周1次,監(jiān)測的水質(zhì)指標包括氨氮、總氮、總磷、高錳酸鹽指數(shù)(CODMn)、生化需氧量(BOD5)和溶解氧(DO)共6項,采樣監(jiān)測標準依據(jù)《水和廢水監(jiān)測分析方法》[15]。表1為采樣時間內(nèi)所有11個重點河道斷面6項水質(zhì)指標的監(jiān)測數(shù)據(jù)。

    表1 天津城市河網(wǎng)重點河道斷面主要水質(zhì)指標監(jiān)測數(shù)據(jù) mg/L

    1.2 研究方法

    a. 系統(tǒng)聚類分析法。根據(jù)變量或樣品的親疏程度,逐次聚合,將性質(zhì)最接近的對象結合在一起,直至聚成一類。通常,根據(jù)采樣點和采樣時間進行聚類以分析水質(zhì)的時空變化特征,根據(jù)評價指標進行聚類以分析指標相似性[16-17]。聚類分析采用的是組間連接法,利用兩組數(shù)據(jù)間各變量值平均平方距離作為衡量兩組數(shù)據(jù)差異的基準。

    b. 因子分析法。從多個實測變量中提取出較少的、互不相關的抽象綜合指標(因子),每個變量可用因子的線性組合表示;同時,根據(jù)各因子對變量的影響大小,也可將變量分為等同于因子數(shù)目的類數(shù)。此方法用于提取污染因子和識別污染源[18-19]。因子分析前首先要對數(shù)據(jù)進行KMO檢驗和Barlett球形檢驗,結果合格后才可進行因子分析。

    c. 多元線性回歸分析法。將變量值與因子得分進行多元線性回歸,根據(jù)回歸參數(shù)得到針對各因子的估計值,從而確定因子對各個變量的貢獻。該方法用于計算污染因子對各評價指標的貢獻率[20]。多元統(tǒng)計分析用SPSS完成。

    2 結果與分析

    依據(jù)河道監(jiān)測數(shù)據(jù),利用系統(tǒng)聚類分析、因子分析和多元線性回歸分析等方法,研究分析河流水質(zhì)指標,探討滯緩流型城市河網(wǎng)水質(zhì)分布特征。根據(jù)GB 3838—2002《地表水環(huán)境質(zhì)量標準》,監(jiān)測斷面的水質(zhì)類別以劣Ⅴ類和Ⅴ類為主,劣Ⅴ類斷面占監(jiān)測斷面的比例達到了63%,Ⅴ類斷面占比為25%,表明天津城市河網(wǎng)水質(zhì)污染問題依然十分嚴峻,其中,總氮、氨氮、總磷超標較為突出,這與水體滯緩流特性、呈現(xiàn)出一定的湖泊水質(zhì)特征有直接關系。

    2.1 河道水質(zhì)時空分布特征

    2.1.1 河道水質(zhì)時間分布特征

    從時間聚類結果(圖2)看,雖然各月水質(zhì)以劣Ⅴ類為主,但根據(jù)聚類中心距不同,可以分為4個不同的時段:3、9、10、11月為第Ⅰ階段,5、6、7月為第Ⅱ階段,4月為第Ⅲ階段,8月為第Ⅳ階段。結合各階段主要污染物質(zhì)量濃度指標(圖3),第Ⅰ階段為非汛期,水質(zhì)較好,氨氮、CODMn、BOD5質(zhì)量濃度值顯著低于其他時段;第Ⅲ階段為汛前,BOD5、氨氮質(zhì)量濃度均有一定上升,水質(zhì)有所下降;第Ⅱ階段和第Ⅳ階段為汛期,污染物質(zhì)量濃度均明顯高于前2個階段,非汛期到汛期水質(zhì)變差的趨勢較顯著,其中第Ⅳ階段污染物質(zhì)量濃度更高,氨氮、總氮超標最嚴重。

    圖2 時間聚類結果

    圖3 時間聚類各階段主要污染物質(zhì)量濃度平均值

    2.1.2 河道水質(zhì)空間分布特征

    a. 汛期。按空間聚類結果(圖4(a)),斷面可以分為3類:第Ⅰ類斷面是S1和S2,第Ⅱ類斷面是S3~S8和S10,第Ⅲ類斷面是S9和S11。如圖5(a)所示,第Ⅰ類斷面的CODMn、BOD5、氨氮基本滿足地表水Ⅴ類標準限值,只有氮磷指標超標;第Ⅱ類斷面屬重要景觀河道斷面,汛期水質(zhì)污染較重,只有BOD5、CODMn、氨氮基本滿足地表水Ⅴ類標準限值;第Ⅲ類斷面受外環(huán)沿線污染匯入影響,汛期水質(zhì)最差,所有指標均超過地表水Ⅴ類標準限值,有機污染物和氮磷均超標數(shù)倍。

    圖4 空間聚類結果

    圖5 空間聚類各類斷面主要污染物質(zhì)量濃度平均值

    b. 非汛期。按空間聚類結果(圖4(b)),斷面可以分為3類:第Ⅰ類斷面是S1~S3,第Ⅱ類斷面是S4~S8和S10,第Ⅲ類斷面是S9和S11。如圖5(b)所示,3類斷面水質(zhì)差別顯著,污染程度逐步加重,第Ⅰ類斷面除了總氮外其他指標均滿足地表水Ⅴ類標準限值,第Ⅱ類斷面氮磷污染超過地表水Ⅴ類標準限值,第Ⅲ類斷面各污染物均超過地表水Ⅴ類標準限值,且超標嚴重。

    綜合空間聚類分析結果,天津城市河網(wǎng)水質(zhì)具有明顯的空間分布特征,作為城市補水水源的海河干流斷面水量多,水質(zhì)最好;二級景觀河道流動性比海河差,水量少,受瀝水、污水匯入影響,水質(zhì)較差;市區(qū)外圍的外環(huán)河流動性更低,沿線可能的污染源分布較多,水質(zhì)最差。

    2.2 河道主要污染源定性識別

    采用因子分析法分別對全部斷面以及聚類后的每一類斷面的主要水質(zhì)影響因子進行分析。

    2.2.1 全部斷面主要水質(zhì)影響因子

    對全部監(jiān)測河流斷面進行統(tǒng)計分析,KMO檢驗值為0.621,變量間的偏相關性強,Barlett球形檢驗表明各變量間具有相關性,因此可以使用因子分析法進行水質(zhì)影響因子分析。

    根據(jù)因子分析結果,提取了對水質(zhì)影響較大的2個公因子F1和F2,因子載荷、特征值及方差貢獻率結果見表2。F1的方差貢獻率49.9%,其中氨氮、總氮、總磷、BOD5的因子載荷大,均與F1呈正相關關系,因此,F1主要反映了水體氮磷的污染程度。F2的方差貢獻率為23.3%,其中CODMn、DO的因子載荷較大,且與F2呈正相關關系,反映了水體中有機污染的水平??傮w來說,天津城市河網(wǎng)水質(zhì)狀態(tài)是以氮磷污染為主導因素,其次是其他種類的有機污染。

    表2 全部斷面因子載荷、特征值及方差貢獻率

    因子得分(表3)反映各監(jiān)測斷面的污染狀況,其值越高,說明水質(zhì)越差,外環(huán)河的津漢橋斷面(S11)和津淶公路橋斷面(S9)得分最高,污染最為嚴重,海河干流的2個斷面得分最低,水質(zhì)最好。這一點與空間聚類的分析結果一致。

    表3 監(jiān)測斷面因子得分

    2.2.2 3類斷面的主要水質(zhì)影響因子

    對空間聚類得到的3類斷面進行統(tǒng)計分析,KMO檢驗值均大于0.6,表明變量間的偏相關性強,Barlett球形檢驗表明各變量間具有相關性,因此可使用因子分析法進行水質(zhì)影響因子分析。根據(jù)特征值大于1的原則,3類斷面均提取了F1、F22個公因子,公因子累計方差貢獻率均在80%左右(表4)。

    表4 3類斷面因子載荷、特征值及方差貢獻率

    a. 第Ⅰ類斷面。2個公因子中,F1的方差貢獻率為43.7%,氨氮、總氮、總磷的因子載荷較大,與F1呈正相關關系;F2的方差貢獻率為34.1%,CODMn、BOD5的因子載荷較大,均超過了0.94,與F2呈正相關關系。綜合來看,氮磷污染、有機物排入仍然是海河水質(zhì)變化的主要因素,這與上游來水、沿岸二級河道的水體匯入有關。

    b. 第Ⅱ類斷面。2個公因子中,F1的方差貢獻率為56.8%,氨氮、總氮、總磷因子載荷為主,BOD5因子載荷值也較高,與F1呈正相關關系;F2的方差貢獻率為20.5%,DO為主要因子載荷,與F2呈正相關關系;CODMn則同時受2個公因子影響,均呈中度正相關關系。綜合來看,氮磷污染、可降解有機物污染是水質(zhì)污染的主要因素,且與海河相比,氮磷污染影響更為突出,這與河道長年補水較少、汛期沿岸雨污水集中匯入有直接關系。

    c. 第Ⅲ類斷面。2個公因子中,F1的方差貢獻率為47.1%,CODMn、BOD5、DO的因子載荷較大,CODMn、BOD5的因子載荷與F1呈正相關關系,DO的因子載荷與F1呈負相關關系;F2的方差貢獻率為30.4%,因子載荷較大的是總氮、總磷、氨氮,均與F2呈正相關關系。綜合來看,對外環(huán)河而言,生活、工業(yè)污染帶來的有機污染影響是河道斷面水質(zhì)變化的關鍵因子,氮磷的影響相對較弱。

    2.3 主要因子對水質(zhì)指標的貢獻率

    根據(jù)因子分析結果,采用多元線性回歸分析法,計算得到3類斷面不同的公因子對各水質(zhì)指標的貢獻率如表5所示。

    回歸過程中,相關系數(shù)R在0.886~0.988之間,說明回歸分析有統(tǒng)計學意義;各水質(zhì)指標估計值與實測值的比值在0.991~1.009之間,也證明了回歸結果的準確性。結果顯示,第Ⅰ類斷面,F1主要影響氨氮、總氮、總磷,F2主要影響CODMn、BOD5、DO;第Ⅱ類斷面,F1主要影響氨氮、總氮、總磷和BOD5,F2主要影響CODMn、DO;第Ⅲ類斷面,F1主要影響CODMn、BOD5、DO,F2主要影響氨氮、總氮、總磷。

    表5 公因子對水質(zhì)指標的貢獻率 %

    2.4 污染源分析

    綜合上述分析,氮磷污染和有機污染匯入及其他不確定因子共同影響著天津城市河網(wǎng)水質(zhì),其中,氮磷污染物是河網(wǎng)水質(zhì)的主要污染因子,有機污染物的影響也較為突出。海河、津河、衛(wèi)津河等主要河道的水域污染因子分布均呈類似的特征,這與城市河網(wǎng)天然補水較少、水體滯緩流、汛期雨污水匯入有密切關系,而在城市外圍河流區(qū)域,有機污染則超過了氮磷污染,成為首要污染因素,這與城市周邊污水收集設施不完善、各類污染源排入較多有直接關系。

    結合天津市實際情況,氮磷污染的主要污染源來自于汛期的雨污水。天津市城區(qū)70%以上的降水量集中在汛期,雨水通過各種途徑進入城區(qū)內(nèi)的景觀河道,同時大量污水也通過溢流井及其合流制泵站、管網(wǎng)進入了景觀河道,導致河道水質(zhì)在汛期氮磷污染嚴重。有機污染的主要污染源來自于工業(yè)點源,且多數(shù)未經(jīng)處理在城市外圍直接排入,雖然總量小于汛期污染匯入,但其有機污染物濃度較高,導致河道水質(zhì)在非汛期有機污染嚴重。

    3 水質(zhì)改善措施

    根據(jù)水質(zhì)時空分布特征及主要污染源的分析,結合天津市的具體特點,可采取以下水質(zhì)改善措施:

    a. 加快推進城市河網(wǎng)水系連通與循環(huán)。從水質(zhì)對比可以看出,海河水質(zhì)顯著優(yōu)于津河、衛(wèi)津河和外環(huán)河。而海河作為城市河網(wǎng)水系的核心,與所有河道都可以實現(xiàn)連通。目前,在局部區(qū)域還存在阻斷的情況,水系沒有實現(xiàn)完全的連通,因此,加快推進城市河網(wǎng)水系連通與循環(huán),對于景觀河道的水質(zhì)改善會起到積極的作用。特別是在非汛期,通過海河給津河、衛(wèi)津河等重點景觀河道補水,較好水質(zhì)水體的進入對其他景觀河道的水質(zhì)和生態(tài)將會帶來顯著的提升和改善。

    b. 推動雨污合流片區(qū)的改造。天津城區(qū)中仍存在部分雨污合流片區(qū),雨污水的混合流入使得這些區(qū)域的景觀河道水質(zhì)會迅速惡化,特別在汛期,水體黑臭,富營養(yǎng)化頻發(fā),景觀效果極差。為此,應積極推動雨污合流片區(qū)的改造,減少污水的直接匯入,將污水引入污水處理廠統(tǒng)一處理,從而從根本上解決津河、衛(wèi)津河等部分景觀河道汛期氮磷負荷過高的問題。

    c. 采取生物生態(tài)等多種手段改善城市景觀河道水質(zhì)。由于長期水質(zhì)較差,單靠外來補水和雨污合流制改造不能在短期內(nèi)改善津河、衛(wèi)津河等河道水質(zhì),可采取底泥清淤、生物投菌、生態(tài)浮島等技術來提升河道水體自凈能力,加速水質(zhì)的改善。其中,底泥清淤屬物理手段,操作簡單,可以直接減少污染層底泥厚度,降低底泥中污染物釋放速度。生物投菌方法屬化學手段,將生物制劑投入富營養(yǎng)化突出的河段,可迅速沉降水體中的磷,對水質(zhì)短時間的改善有顯著效果。生態(tài)浮島技術屬生態(tài)治理技術,利用植物對水體中污染物的吸收降解作用達到減污凈水效果。這些技術都經(jīng)過了大量的實踐檢驗,效果突出。

    d. 利用濕地技術改善外環(huán)河水質(zhì)。不同于城市內(nèi)部,外環(huán)河沿線有豐富的土地資源,建議開展人工濕地建設,通過在外環(huán)河沿線構建多個大型人工濕地,對外環(huán)河水體進行異位處理,來改善水質(zhì)。目前,外環(huán)河津淶公路橋附近已建成一個大型人工濕地,從初期運行情況看,對局部河段水質(zhì)特別是氮磷有顯著改善效果。

    e. 加強對外環(huán)河排污口監(jiān)管。外環(huán)河水質(zhì)較差,主要污染源是有機污染匯入,這與外環(huán)河沿線工業(yè)點源有密切關系。結合天津市開展的排污口治理,應加強對保留的排污口門的監(jiān)督和管理,控制排污量,調(diào)控排污對河道水質(zhì)的影響。

    4 結 論

    a. 根據(jù)河道水質(zhì)時間分布特征,年內(nèi)可劃分為4個階段:3、9、10、11月為第Ⅰ階段,4月為第Ⅱ階段,5、6、7月為第Ⅲ階段,8月為第Ⅳ階段,前2個階段水質(zhì)總體較好,后2個階段水質(zhì)較差,氨氮、總磷質(zhì)量濃度超標嚴重。

    b. 根據(jù)河道水質(zhì)空間分布特征,將11個監(jiān)測斷面分為3類:汛期海河2個斷面為第Ⅰ類斷面,津河、衛(wèi)津河6個斷面及外環(huán)河七號橋斷面為第Ⅱ類斷面,外環(huán)河其他2個斷面為第Ⅲ類斷面;非汛期與汛期相似,只是津河三元村斷面由汛期第Ⅱ類斷面變?yōu)榱说冖耦悢嗝?。從第Ⅰ類斷面到第Ⅲ類斷?水質(zhì)逐步下降。

    c. 天津滯緩流城市河網(wǎng)水系水質(zhì)以氮磷污染為主,其他有機污染次之,呈現(xiàn)出典型的湖泊水質(zhì)特征。其中,第Ⅰ類斷面以氮磷污染為主,第Ⅱ類斷面以氮磷污染、生化降解有機污染為主,第Ⅲ類斷面則主要受有機污染影響。3類斷面6個水質(zhì)指標分別不同程度地受到水體氮磷污染水平、有機污染水平以及其他不確定因素的影響。

    [ 1 ] 陳興茹.城市河流生態(tài)修復淺議[J].中國水利水電科學研究院學報,2006,4(3):226-231.(CHEN Xinru.Brief discussion on urban river course eco-restoration[J].Journal of China Institute of Water Resources and Hydropower Research,2006,4(3):226-231.(in Chinese))

    [ 2 ] 王秀朵.北方缺水城市景觀水體污染控制[J].給水排水,2011,37(7):1-3.(WANG Xiuduo.Scenic water pollution prevention measures of city water shortage in North China[J].Water & Wastewater Engineering,2011,37(7):1-3.(in Chinese))

    [ 3 ] 楊學福,王蕾,關建玲,等.基于多元統(tǒng)計分析的渭河西咸段水質(zhì)評價[J].環(huán)境工程學報,2016,10(3):1560-1565.(YANG Xuefu,WANG Lei,GUAN Jianling,et al.Comprehensive assessment of water quality in Xi’an-Xianyang section of Weihe River based on multivariate analysis method[J].Chinese Journal of Environmental Engineering,2016,10(3):1560-1565.(in Chinese))

    [ 4 ] 程琳琳,楊開宇,杜鵑,等.河北省河流水質(zhì)時空變化特征及受人類活動影響的分析[J].水資源與水工程學報,2015,26(1):1-7.(CHENG Linlin,YANG Kaiyu,DU Juan,et al.Feature of spatial land and temporal variation of water quality and influence of human activities on it in Hebei Province[J].Journal of Water Resources and Water Engineering,2015,26(1):1-7.(in Chinese))

    [ 5 ] 葛懷鳳,秦大庸,周祖昊,等.基于污染遷移轉化過程的海河干流天津段污染關鍵源區(qū)及污染類別分析[J].水利學報,2011,42(1):61-67.(GE Huaifeng,QIN Dayong,ZHOU Zuhao,et al.Analysis of key source areas and pollution type in the lower Haihe River based on pollution loading movement and transformation[J].Journal of Hydraulie Engineering,2011,42(1):61-67.(in Chinese))

    [ 6 ] 王剛,李兆富,萬榮榮,等.基于多元統(tǒng)計分析方法的西苕溪流域水質(zhì)時空變化研究[J].農(nóng)業(yè)環(huán)境科學學報,2015,34(9):1797-1803.(WANG Gang,LI Zhaofu,WAN Rongrong,et al.Analysis of temporal and spatial variations in water quality of Xitiaoxi watershed using multivariate statistical techniques[J].Journal of Agro-Environment Science,2015,34(9):1797-1803.(in Chinese))

    [ 7 ] BU Hongmei,TAN Xiang,LI Siyue,et al.Water quality assessment of the Jinshui River (China) using multivariate statistical techniques[J].Environmental Earth Sciences,2010,60(8):1631-1639.

    [ 8 ] 程學寧,湯云,盧毅敏.基于多元統(tǒng)計分析的閩江水質(zhì)時空變化特征[J].水資源與水工程學報,2016,27(6):89-94. (CHENG Xuening,TANG Yun,LU Yimin.Spatial and temporal charactersitics of water quality using multivariate statistical analysis in Min River[J].Journal of Water Resources and Water Engineering,2016,27(6):89-94.(in Chinese))

    [ 9 ] BERKANT O,FATIH E.Monitoring water quality and quantity of national watersheds in Turkey[J].Environmental Monitoring and Assessment,2007,133(1/2/3):215-229.

    [10] ANDRé S H,GUY B,SIMON C,et al.Multivariate analysis of water quality in the richibucto drainage basin (New Brunswick, Canada)[J].Journal of American Water Resources Association,2007,40(3):691-703.

    [11] BAHMAN J A,KANEYUKI N.Modeling the linkage between river water quality and landscape metrics in the Chugoku District of Japan[J].Water Resources Management,2008,23(5):931-956.

    [12] ZHANG Qi,LI Zhongwu,ZENG Guangming,et al.Assessment of surface waterquality using multivariate statistical techniques in red soil hillyregion:a case study of Xiangjiang Watershed,China[J].Environmental Monitoring and Assessment,2009,152 (1/2/3/4):123-131.

    [13] 徐華山,徐宗學,唐芳芳,等.漳衛(wèi)南運河流域水質(zhì)時空變化特征及其污染源識別[J].環(huán)境科學,2012,33(2):359-369.(XU Huashan,XU Zongxue,TANG Fangfang,et al.Spatiotemporal variation analysis and identification of water pollution sources in the Zhangweinan River Basin[J].Environmental Science,2012,33(2):359-369.(in Chinese))

    [14] 李婉,張娜,吳芳芳.北京轉河河岸帶生態(tài)修復對河流水質(zhì)的影響[J].環(huán)境科學,2011,32(1):80-87.(LI Wan,ZHANG Na,WU Fangfan.Influence of ecological restoration of riparian zone on water quality of Zhuanhe River in Beijing[J].Environmental Science,2011,32(1):80-87.(in Chinese))

    [15] 國家環(huán)境保護總局《水和廢水監(jiān)測分析方法》編委會.水和廢水監(jiān)測分析方法[M].4版.北京:中國環(huán)境科學出版社,2002.

    [16] SHRESTHA S,KAZAMA F.Assessment of surface water quality using multivariate statistical techniques:a case study of the Fuji River Basin,Japan[J].Environmental Modelling & Software,2007,22(4):464-475.

    [17] SUNDARAY S K,PANDA U C,NAYAK B B,et al.Multivariate statistical techniques for the evaluation of spatial and temporal variations in water quality of the Mahanadi river-estuarine system(India):a case study[J].Environmental Geochemistry and Health,2006,28(4):317-330.

    [18] 余輝,燕姝雯,徐軍.湯浦水庫及入庫支流水質(zhì)時空變化特征與影響因素分析[J].水生態(tài)學雜志,2013,34(5):9-15.(YU Hui,YAN Shuwen,XU Jun.Characteristics and influencing factors analysis of spatio-temporal variations of water quality in Tangpu Reservoir and its tributaries[J].Journal of Hydroecology,2013,34(5):9-15.(in Chinese))

    [19] JOHNSON H O,GUPTA S C,VECCHIA A V,et al.Assessment of water quality trends in the Minnesota River using non-parametricand parametric methods[J].Journal of Environmental Quality,2009,38(3):1018-1030.

    [20] PALMA P,ALVARENGA P,PALMA V L,et al.Assessment of anthropogenic sources of water pollution using multivariatestatistical techniques:a case study of the Alqueva’s Reservoir,Portugal[J].Environmental Monitoring and Assessment,2010,165(1/2/3/4):539-552.

    Temporal and spatial distribution characteristics of water quality of stagnant river network in Tianjin City, China

    DONG Lixin, BAI Haoyang

    (Tianjin Hydraulic Research Institute, Tianjin 300061, China)

    The urban river network in Tianjin City has stagnant or slow-flow water, which leads to water quality variations different from natural rivers. Multivariate statistical analysis methods, including the cluster analysis method, factor analysis method, and multiple linear regression method, were used to evaluate the water quality data from 11 sites in four main rivers of Tianjin, the Haihe River,Jinhe River, Weijin River, and Waihuai River, for 39 weeks during the period from 2008 to 2009. The results show that the water quality variation exhibited specific temporal and spatial characteristics. Water quality data from May to August were worse than during other months, especially for the ammonia nitrogen and total phosphorous concentrations, which severely exceeded the water quality standards. Water quality deteriorated from the Haihe River to Jinhe River, Weijin River, and, finally, Waihuan River gradually. Nitrogen and phosphorous pollutants were the main pollution sources of the Haihe River, Jinhe River, and Weijin River, while organic pollutants were the main pollution sources of the Waihuan River.

    urban river network; stagnant and slow-flow water; water quality variation; multivariate statistical analysis; Tianjin City

    國家水體污染控制與治理科技重大專項(2012ZX07203-002)

    董立新(1977—),男,高級工程師,碩士,主要從事水資源與水環(huán)境研究。E-mail:dlx@tjhri.com

    10.3880/j.issn.1006-7647.2017.04.002

    X522;X824

    A

    1006-7647(2017)04-0008-06

    2017-01-05 編輯:熊水斌)

    猜你喜歡
    河網(wǎng)外環(huán)氮磷
    基于小世界網(wǎng)絡的海河流域河網(wǎng)結構及功能響應
    水科學進展(2023年4期)2023-10-07 11:23:44
    玻利維亞擬建新的氮磷鉀肥料工廠
    深圳外環(huán)高速公路通車
    石油瀝青(2021年1期)2021-01-11 05:41:31
    基于PSR模型的上海地區(qū)河網(wǎng)脆弱性探討
    鞏膜外環(huán)扎術治療復發(fā)性視網(wǎng)膜脫離的臨床觀察
    不同引水水源對平原河網(wǎng)影響分析
    MIKE11模型在城市河網(wǎng)生態(tài)調(diào)水工程中的應用
    某型機低壓渦輪外環(huán)釬焊技術
    焊接(2015年5期)2015-07-18 11:03:41
    植物需要氮磷鉀
    不同氮磷水平下AM真菌對棉花的侵染率及氮磷含量的影響
    免费观看在线日韩| 国产乱来视频区| 亚洲国产精品成人久久小说| 午夜精品国产一区二区电影| 亚洲欧美成人精品一区二区| 视频在线观看一区二区三区| 丰满乱子伦码专区| 午夜福利视频精品| 特大巨黑吊av在线直播| a级片在线免费高清观看视频| 精品少妇黑人巨大在线播放| 亚洲精品自拍成人| 一级毛片电影观看| 久久人人爽人人片av| 久久婷婷青草| 免费av不卡在线播放| 最黄视频免费看| 伦理电影免费视频| 老司机亚洲免费影院| 亚洲少妇的诱惑av| 高清午夜精品一区二区三区| 少妇人妻 视频| 日韩强制内射视频| 少妇 在线观看| 伊人久久国产一区二区| 国产免费现黄频在线看| 七月丁香在线播放| 色视频在线一区二区三区| 少妇熟女欧美另类| 日韩成人伦理影院| 亚洲精品,欧美精品| 男女国产视频网站| 成人二区视频| 久久精品久久精品一区二区三区| 国产精品熟女久久久久浪| 91成人精品电影| 国产精品一区二区在线观看99| 日韩中字成人| 精品国产国语对白av| 精品99又大又爽又粗少妇毛片| 亚洲,欧美,日韩| 日韩熟女老妇一区二区性免费视频| av福利片在线| 一区二区三区乱码不卡18| 亚洲精品一区蜜桃| 久久久久久久久久久免费av| 97超视频在线观看视频| 国产精品久久久久久精品电影小说| 国产色爽女视频免费观看| 色网站视频免费| 哪个播放器可以免费观看大片| 久久ye,这里只有精品| 国产片特级美女逼逼视频| av在线app专区| 亚洲成色77777| 菩萨蛮人人尽说江南好唐韦庄| 在线观看免费高清a一片| 中文字幕免费在线视频6| 国产一区二区三区综合在线观看 | 国产一区二区三区综合在线观看 | 黄色怎么调成土黄色| 精品少妇内射三级| 99九九在线精品视频| 全区人妻精品视频| 男女边摸边吃奶| 亚洲精品日韩av片在线观看| 国产极品粉嫩免费观看在线 | 男女国产视频网站| 亚洲伊人久久精品综合| 亚洲av欧美aⅴ国产| 久久人妻熟女aⅴ| 精品国产一区二区三区久久久樱花| 91午夜精品亚洲一区二区三区| 欧美日韩视频精品一区| 美女国产视频在线观看| 国产亚洲午夜精品一区二区久久| 午夜免费鲁丝| 最近手机中文字幕大全| 丝袜脚勾引网站| 国产精品久久久久久久久免| 99久久人妻综合| 乱码一卡2卡4卡精品| 伦精品一区二区三区| 边亲边吃奶的免费视频| 少妇丰满av| 久久影院123| 成人18禁高潮啪啪吃奶动态图 | 亚洲第一av免费看| 黄色欧美视频在线观看| 亚洲精品美女久久av网站| 大话2 男鬼变身卡| 26uuu在线亚洲综合色| 成人影院久久| 精品一区二区三区视频在线| 免费黄频网站在线观看国产| 欧美+日韩+精品| 国产精品 国内视频| 亚洲国产欧美在线一区| 春色校园在线视频观看| 亚洲国产精品国产精品| 欧美激情国产日韩精品一区| 九草在线视频观看| 亚洲一级一片aⅴ在线观看| 免费日韩欧美在线观看| 人妻一区二区av| 如日韩欧美国产精品一区二区三区 | 少妇被粗大猛烈的视频| 美女内射精品一级片tv| 91成人精品电影| 国产成人精品久久久久久| av又黄又爽大尺度在线免费看| 在线 av 中文字幕| 午夜视频国产福利| 国产 一区精品| 18禁在线播放成人免费| 亚洲国产成人一精品久久久| 国产女主播在线喷水免费视频网站| 免费久久久久久久精品成人欧美视频 | 热re99久久国产66热| 精品亚洲乱码少妇综合久久| 久久久久久久久久人人人人人人| 免费高清在线观看日韩| 久久精品久久久久久久性| 熟妇人妻不卡中文字幕| 纯流量卡能插随身wifi吗| 老女人水多毛片| 日韩,欧美,国产一区二区三区| 中文字幕人妻丝袜制服| 欧美精品一区二区免费开放| av免费在线看不卡| 国产精品国产三级国产av玫瑰| 国产精品三级大全| 国产成人免费无遮挡视频| 一区二区三区免费毛片| 51国产日韩欧美| 高清毛片免费看| 亚洲欧美成人精品一区二区| 欧美最新免费一区二区三区| 亚洲精品视频女| 少妇被粗大猛烈的视频| 激情五月婷婷亚洲| 王馨瑶露胸无遮挡在线观看| 久久精品夜色国产| 精品视频人人做人人爽| 波野结衣二区三区在线| 极品人妻少妇av视频| 婷婷色综合www| 国产av一区二区精品久久| 欧美日韩在线观看h| 王馨瑶露胸无遮挡在线观看| 国产精品不卡视频一区二区| 一本—道久久a久久精品蜜桃钙片| av女优亚洲男人天堂| 伦精品一区二区三区| 日韩制服骚丝袜av| 三级国产精品片| 国产成人精品福利久久| 日本猛色少妇xxxxx猛交久久| 欧美丝袜亚洲另类| 亚洲欧美一区二区三区黑人 | 日本av免费视频播放| 丰满少妇做爰视频| 尾随美女入室| 啦啦啦在线观看免费高清www| 亚洲精品第二区| 国产一区二区三区综合在线观看 | 最后的刺客免费高清国语| 亚洲丝袜综合中文字幕| 成人国产麻豆网| 99久国产av精品国产电影| 蜜桃在线观看..| 亚洲欧洲国产日韩| 国产精品国产三级国产专区5o| 大香蕉久久网| 男女高潮啪啪啪动态图| 午夜福利视频在线观看免费| 18在线观看网站| 国产免费视频播放在线视频| 国产免费又黄又爽又色| 亚洲精品自拍成人| 日日摸夜夜添夜夜爱| 亚洲欧美色中文字幕在线| 免费高清在线观看视频在线观看| 日韩免费高清中文字幕av| 人妻 亚洲 视频| av免费观看日本| 99久国产av精品国产电影| 国产成人一区二区在线| 啦啦啦在线观看免费高清www| 久久精品久久久久久久性| 色网站视频免费| 狂野欧美激情性xxxx在线观看| 三上悠亚av全集在线观看| 精品熟女少妇av免费看| 午夜福利视频精品| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲无线观看免费| 成人国产av品久久久| 免费人妻精品一区二区三区视频| 人妻人人澡人人爽人人| 伊人久久精品亚洲午夜| 极品人妻少妇av视频| 日日爽夜夜爽网站| 18+在线观看网站| 亚洲国产日韩一区二区| 男女啪啪激烈高潮av片| 人人妻人人澡人人看| 久久人人爽av亚洲精品天堂| 97精品久久久久久久久久精品| 免费少妇av软件| 蜜桃在线观看..| 自线自在国产av| 80岁老熟妇乱子伦牲交| 精品熟女少妇av免费看| 如日韩欧美国产精品一区二区三区 | 99国产综合亚洲精品| 亚洲精品亚洲一区二区| 2018国产大陆天天弄谢| 午夜91福利影院| 丝袜在线中文字幕| 18+在线观看网站| 久久99一区二区三区| 久久久久久久久大av| 男女边吃奶边做爰视频| 欧美 亚洲 国产 日韩一| 国产欧美日韩一区二区三区在线 | 日日爽夜夜爽网站| 一个人免费看片子| 日韩精品有码人妻一区| 青春草国产在线视频| 精品国产露脸久久av麻豆| 久久久久人妻精品一区果冻| 黄色欧美视频在线观看| 美女国产视频在线观看| 成人二区视频| 天天躁夜夜躁狠狠久久av| 国产精品久久久久久av不卡| 日本黄色日本黄色录像| 久久ye,这里只有精品| 久久久国产欧美日韩av| 九九久久精品国产亚洲av麻豆| 久久精品国产亚洲网站| 母亲3免费完整高清在线观看 | 日韩制服骚丝袜av| 特大巨黑吊av在线直播| 日韩一本色道免费dvd| 日韩欧美精品免费久久| 亚洲中文av在线| 欧美另类一区| 青春草视频在线免费观看| 又大又黄又爽视频免费| 99热全是精品| 国产永久视频网站| 精品国产国语对白av| 麻豆成人av视频| 久久久国产一区二区| 国产深夜福利视频在线观看| 日本色播在线视频| 亚洲国产精品成人久久小说| 精品人妻一区二区三区麻豆| 国产精品嫩草影院av在线观看| 老女人水多毛片| 国产精品偷伦视频观看了| 国产成人a∨麻豆精品| 狠狠精品人妻久久久久久综合| av视频免费观看在线观看| 男女啪啪激烈高潮av片| 亚洲情色 制服丝袜| 999精品在线视频| 尾随美女入室| 免费播放大片免费观看视频在线观看| 在线观看美女被高潮喷水网站| 午夜精品国产一区二区电影| 国产精品一国产av| 国产日韩欧美视频二区| 22中文网久久字幕| 国产精品久久久久久精品电影小说| 亚洲成人av在线免费| 亚洲精品久久午夜乱码| 伦精品一区二区三区| 精品久久国产蜜桃| 免费人成在线观看视频色| 国产精品秋霞免费鲁丝片| 免费看av在线观看网站| av播播在线观看一区| 美女国产高潮福利片在线看| 熟女电影av网| videossex国产| 国产欧美亚洲国产| 一区二区三区四区激情视频| 人人妻人人爽人人添夜夜欢视频| av一本久久久久| 欧美一级a爱片免费观看看| 日韩精品有码人妻一区| 中文字幕免费在线视频6| 91国产中文字幕| 亚洲欧美一区二区三区国产| 三级国产精品欧美在线观看| 美女cb高潮喷水在线观看| 人人妻人人澡人人爽人人夜夜| 18禁在线播放成人免费| 尾随美女入室| 美女xxoo啪啪120秒动态图| 日韩欧美一区视频在线观看| 97超视频在线观看视频| 国产一区有黄有色的免费视频| 亚洲国产av影院在线观看| 丝袜在线中文字幕| 欧美人与性动交α欧美精品济南到 | 最近中文字幕2019免费版| 久久精品久久精品一区二区三区| 国产在线视频一区二区| 国产高清有码在线观看视频| 一个人看视频在线观看www免费| 黑人猛操日本美女一级片| 精品一区二区三卡| 少妇的逼好多水| 午夜影院在线不卡| 777米奇影视久久| 国产在线一区二区三区精| av福利片在线| 狠狠精品人妻久久久久久综合| 亚洲美女黄色视频免费看| 日韩一区二区三区影片| 久久99蜜桃精品久久| 99热6这里只有精品| 国产成人精品无人区| 中文字幕精品免费在线观看视频 | 少妇人妻久久综合中文| 极品少妇高潮喷水抽搐| 精品一区二区免费观看| 99久久人妻综合| 最新中文字幕久久久久| 久久久久久久国产电影| 精品国产一区二区久久| 九九爱精品视频在线观看| 男人操女人黄网站| 晚上一个人看的免费电影| 欧美精品一区二区免费开放| 久热久热在线精品观看| 少妇丰满av| 欧美日韩亚洲高清精品| 观看av在线不卡| 久久99热6这里只有精品| 亚洲一区二区三区欧美精品| 国产乱来视频区| 久久久午夜欧美精品| 九色亚洲精品在线播放| 婷婷色麻豆天堂久久| 精品人妻熟女av久视频| 这个男人来自地球电影免费观看 | 精品国产一区二区久久| 日日啪夜夜爽| 日本午夜av视频| 丝袜脚勾引网站| 在现免费观看毛片| 国模一区二区三区四区视频| 乱码一卡2卡4卡精品| 高清黄色对白视频在线免费看| 丰满迷人的少妇在线观看| 18禁动态无遮挡网站| 高清午夜精品一区二区三区| 欧美另类一区| 青春草国产在线视频| 全区人妻精品视频| 亚洲欧洲日产国产| 熟妇人妻不卡中文字幕| 你懂的网址亚洲精品在线观看| 中文字幕最新亚洲高清| 国产一级毛片在线| 免费观看av网站的网址| 丰满乱子伦码专区| 中文字幕制服av| 久久人妻熟女aⅴ| 在线观看国产h片| 91成人精品电影| 街头女战士在线观看网站| 久久久久精品性色| 国产高清不卡午夜福利| 精品久久蜜臀av无| 999精品在线视频| 亚洲国产精品一区二区三区在线| 啦啦啦视频在线资源免费观看| 欧美日韩av久久| 欧美丝袜亚洲另类| 永久网站在线| 亚洲av电影在线观看一区二区三区| av专区在线播放| 2018国产大陆天天弄谢| 夜夜爽夜夜爽视频| 久久久久视频综合| 妹子高潮喷水视频| 久久久亚洲精品成人影院| 性色avwww在线观看| 十八禁网站网址无遮挡| 18禁裸乳无遮挡动漫免费视频| 另类精品久久| 成人毛片60女人毛片免费| 欧美日韩国产mv在线观看视频| 日韩,欧美,国产一区二区三区| 精品一区在线观看国产| 亚洲国产欧美日韩在线播放| 精品国产一区二区三区久久久樱花| 国精品久久久久久国模美| 乱码一卡2卡4卡精品| 超色免费av| 色婷婷久久久亚洲欧美| 日本免费在线观看一区| 精品国产国语对白av| 女性生殖器流出的白浆| 一级二级三级毛片免费看| 国产精品人妻久久久久久| 亚洲av福利一区| tube8黄色片| 色婷婷av一区二区三区视频| 热re99久久精品国产66热6| 欧美xxⅹ黑人| 一本大道久久a久久精品| .国产精品久久| av在线播放精品| 国产不卡av网站在线观看| 国产精品免费大片| 精品人妻在线不人妻| 亚洲综合色网址| 夜夜爽夜夜爽视频| 亚洲精品色激情综合| 国产日韩欧美亚洲二区| 精品一区二区免费观看| 亚洲精品第二区| 三级国产精品欧美在线观看| videossex国产| 91精品三级在线观看| 五月玫瑰六月丁香| 97在线人人人人妻| 美女主播在线视频| 国产精品国产三级国产专区5o| 国产精品 国内视频| 久久久久网色| 欧美三级亚洲精品| 视频中文字幕在线观看| 在线播放无遮挡| 国产亚洲午夜精品一区二区久久| av在线老鸭窝| 草草在线视频免费看| 特大巨黑吊av在线直播| 在线观看免费日韩欧美大片 | 午夜免费男女啪啪视频观看| 国国产精品蜜臀av免费| 午夜福利视频在线观看免费| 十分钟在线观看高清视频www| 香蕉精品网在线| 欧美精品一区二区免费开放| 韩国高清视频一区二区三区| 国产成人精品无人区| 久久国内精品自在自线图片| 国产在线免费精品| 免费大片18禁| 少妇被粗大的猛进出69影院 | 大香蕉久久成人网| av电影中文网址| 国产探花极品一区二区| 一级毛片 在线播放| 日韩一区二区三区影片| 国产高清不卡午夜福利| 国产精品女同一区二区软件| 久久久国产欧美日韩av| 日韩中字成人| 下体分泌物呈黄色| 九草在线视频观看| av天堂久久9| 又黄又爽又刺激的免费视频.| 亚洲欧美中文字幕日韩二区| 中国美白少妇内射xxxbb| 自拍欧美九色日韩亚洲蝌蚪91| 99久久人妻综合| 国产精品秋霞免费鲁丝片| 麻豆乱淫一区二区| 天天躁夜夜躁狠狠久久av| 久久国内精品自在自线图片| 精品国产露脸久久av麻豆| 成人无遮挡网站| 亚洲经典国产精华液单| 中国国产av一级| 午夜福利,免费看| 亚洲欧洲精品一区二区精品久久久 | 人人澡人人妻人| 久久婷婷青草| 中文精品一卡2卡3卡4更新| 自线自在国产av| 国产午夜精品久久久久久一区二区三区| 97精品久久久久久久久久精品| 天天操日日干夜夜撸| 能在线免费看毛片的网站| 免费大片黄手机在线观看| 美女xxoo啪啪120秒动态图| 日本91视频免费播放| 少妇丰满av| 婷婷色av中文字幕| 欧美亚洲 丝袜 人妻 在线| av不卡在线播放| 色网站视频免费| 美女内射精品一级片tv| 一个人看视频在线观看www免费| 男女国产视频网站| 国产日韩欧美在线精品| 亚洲欧洲国产日韩| 亚洲国产精品国产精品| 国产成人免费观看mmmm| videos熟女内射| 亚洲伊人久久精品综合| 日日啪夜夜爽| 精品一区二区免费观看| 中文字幕免费在线视频6| 精品国产一区二区久久| 蜜桃久久精品国产亚洲av| 欧美日韩精品成人综合77777| 另类精品久久| 成人手机av| 久久99热6这里只有精品| 欧美日韩成人在线一区二区| 午夜激情福利司机影院| 亚洲精品久久久久久婷婷小说| 免费黄色在线免费观看| 大话2 男鬼变身卡| 亚洲四区av| 大话2 男鬼变身卡| 中文乱码字字幕精品一区二区三区| 午夜福利在线观看免费完整高清在| 婷婷色综合www| 大香蕉久久网| 在线观看免费高清a一片| 欧美日本中文国产一区发布| 国产黄片视频在线免费观看| 午夜福利,免费看| 欧美性感艳星| 欧美老熟妇乱子伦牲交| 久久99热这里只频精品6学生| 国产在线视频一区二区| 美女内射精品一级片tv| 大话2 男鬼变身卡| 久久人人爽人人片av| 婷婷色综合大香蕉| 99久久人妻综合| av又黄又爽大尺度在线免费看| 黄色欧美视频在线观看| 久久ye,这里只有精品| 亚洲久久久国产精品| 欧美日韩视频高清一区二区三区二| 亚洲久久久国产精品| 黑人欧美特级aaaaaa片| 国产伦理片在线播放av一区| 男女国产视频网站| 两个人免费观看高清视频| 国产精品国产三级专区第一集| 午夜激情av网站| 国产有黄有色有爽视频| 欧美 亚洲 国产 日韩一| 日本-黄色视频高清免费观看| 中国国产av一级| 日本黄色日本黄色录像| 亚洲国产av影院在线观看| 免费黄网站久久成人精品| 22中文网久久字幕| 丰满少妇做爰视频| 国产在线免费精品| 免费观看无遮挡的男女| 黄色欧美视频在线观看| 欧美性感艳星| 伦理电影大哥的女人| 我要看黄色一级片免费的| 美女大奶头黄色视频| 婷婷色av中文字幕| 精品国产一区二区久久| 精品亚洲成a人片在线观看| 色吧在线观看| 久久国产精品男人的天堂亚洲 | 美女xxoo啪啪120秒动态图| 午夜福利,免费看| 黑人欧美特级aaaaaa片| 午夜福利,免费看| 18禁在线播放成人免费| 久久亚洲国产成人精品v| 久久午夜综合久久蜜桃| 人人妻人人爽人人添夜夜欢视频| 久久久久人妻精品一区果冻| 久久精品国产鲁丝片午夜精品| 日韩人妻高清精品专区| 亚洲国产精品成人久久小说| 精品酒店卫生间| 精品一区在线观看国产| 欧美变态另类bdsm刘玥| 精品一区在线观看国产| 国产又色又爽无遮挡免| 制服人妻中文乱码| 日韩在线高清观看一区二区三区| av天堂久久9| 久久av网站| 日本av手机在线免费观看| 免费看光身美女| 日日爽夜夜爽网站| av.在线天堂| 51国产日韩欧美| 久久精品国产自在天天线| 亚洲图色成人| 黑人高潮一二区| 免费看不卡的av| 精品少妇久久久久久888优播| 成人免费观看视频高清| 国产免费视频播放在线视频| 你懂的网址亚洲精品在线观看| 最近2019中文字幕mv第一页| 久久99蜜桃精品久久| 精品亚洲乱码少妇综合久久| 国产亚洲欧美精品永久| 97超碰精品成人国产|