• <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真菌對棉花的侵染率及氮磷含量的影響
    国产精品成人在线| 纯流量卡能插随身wifi吗| 国产亚洲最大av| 亚洲天堂av无毛| 精品国产露脸久久av麻豆| 午夜福利乱码中文字幕| 女的被弄到高潮叫床怎么办| 99久国产av精品国产电影| 卡戴珊不雅视频在线播放| 久久影院123| 自拍欧美九色日韩亚洲蝌蚪91| 精品少妇久久久久久888优播| av片东京热男人的天堂| 汤姆久久久久久久影院中文字幕| 久热这里只有精品99| 亚洲精品一区蜜桃| 久久久久久久国产电影| 男人爽女人下面视频在线观看| 高清不卡的av网站| 国产高清国产精品国产三级| 国产成人91sexporn| 街头女战士在线观看网站| 男女免费视频国产| 亚洲经典国产精华液单| 十八禁网站网址无遮挡| 欧美变态另类bdsm刘玥| 精品一区二区三卡| 国产亚洲精品第一综合不卡| 18禁观看日本| 国产欧美日韩一区二区三区在线| 激情视频va一区二区三区| 大片免费播放器 马上看| 综合色丁香网| 国产亚洲最大av| 黑人猛操日本美女一级片| 日韩三级伦理在线观看| 久久这里只有精品19| 如日韩欧美国产精品一区二区三区| 满18在线观看网站| 99精国产麻豆久久婷婷| 九草在线视频观看| 高清视频免费观看一区二区| 国产一区二区三区综合在线观看| 国产片内射在线| 成年女人毛片免费观看观看9 | 在线精品无人区一区二区三| 国产精品国产av在线观看| 天天影视国产精品| 久久久a久久爽久久v久久| 精品人妻一区二区三区麻豆| 婷婷色av中文字幕| 少妇人妻精品综合一区二区| 国产 精品1| 曰老女人黄片| 下体分泌物呈黄色| 国产精品99久久99久久久不卡 | 精品国产国语对白av| av在线播放精品| 日韩成人av中文字幕在线观看| 国产人伦9x9x在线观看 | 97在线人人人人妻| 久久久久久久大尺度免费视频| 蜜桃国产av成人99| 一级a爱视频在线免费观看| 在线免费观看不下载黄p国产| 下体分泌物呈黄色| 亚洲欧美一区二区三区国产| 亚洲欧美一区二区三区国产| 一级黄片播放器| 国产精品蜜桃在线观看| 纯流量卡能插随身wifi吗| 曰老女人黄片| a 毛片基地| 久久精品夜色国产| 在线免费观看不下载黄p国产| 国产乱来视频区| 老汉色∧v一级毛片| av线在线观看网站| 熟妇人妻不卡中文字幕| 日本欧美视频一区| 少妇的逼水好多| 考比视频在线观看| 欧美老熟妇乱子伦牲交| 久久av网站| 看免费成人av毛片| 国产成人一区二区在线| 国产一区二区三区综合在线观看| 国产一区二区三区综合在线观看| 女人久久www免费人成看片| 色网站视频免费| 久久人妻熟女aⅴ| 在线看a的网站| 交换朋友夫妻互换小说| 男的添女的下面高潮视频| 亚洲精品自拍成人| 免费高清在线观看视频在线观看| 欧美精品一区二区免费开放| 一区在线观看完整版| 美女主播在线视频| 国产成人91sexporn| 最近中文字幕2019免费版| 亚洲,欧美,日韩| 欧美日本中文国产一区发布| av有码第一页| 天堂中文最新版在线下载| 男女下面插进去视频免费观看| 交换朋友夫妻互换小说| freevideosex欧美| 一区二区三区精品91| 大香蕉久久成人网| 亚洲天堂av无毛| 亚洲精品,欧美精品| 18在线观看网站| 又大又黄又爽视频免费| 爱豆传媒免费全集在线观看| 各种免费的搞黄视频| 各种免费的搞黄视频| 欧美另类一区| 国产成人av激情在线播放| 亚洲国产欧美网| 一本大道久久a久久精品| 国产在视频线精品| av女优亚洲男人天堂| 久久久久国产精品人妻一区二区| 99久久精品国产国产毛片| 亚洲,一卡二卡三卡| 亚洲国产最新在线播放| 久久久久网色| 久久久久久免费高清国产稀缺| 超碰97精品在线观看| 黑人巨大精品欧美一区二区蜜桃| 国产在线一区二区三区精| 国产亚洲av片在线观看秒播厂| 久久这里有精品视频免费| 两个人免费观看高清视频| 久久精品熟女亚洲av麻豆精品| 精品一品国产午夜福利视频| 欧美+日韩+精品| 岛国毛片在线播放| 亚洲成人av在线免费| 美女大奶头黄色视频| 一级片免费观看大全| 91国产中文字幕| 国产欧美日韩综合在线一区二区| 日韩三级伦理在线观看| 久久久精品94久久精品| 久久人人爽av亚洲精品天堂| 看免费av毛片| 久久精品国产综合久久久| 欧美日韩国产mv在线观看视频| 亚洲av免费高清在线观看| 亚洲国产精品成人久久小说| 少妇 在线观看| 国产黄色视频一区二区在线观看| 一级a爱视频在线免费观看| 免费看av在线观看网站| 男人舔女人的私密视频| 欧美亚洲日本最大视频资源| 日韩精品免费视频一区二区三区| 日韩欧美一区视频在线观看| 国产精品不卡视频一区二区| 亚洲精品一二三| 亚洲精品美女久久久久99蜜臀 | 国产日韩一区二区三区精品不卡| 卡戴珊不雅视频在线播放| 少妇猛男粗大的猛烈进出视频| 亚洲欧美一区二区三区久久| 国产精品嫩草影院av在线观看| 免费黄色在线免费观看| 另类精品久久| 国产在线一区二区三区精| 久久精品熟女亚洲av麻豆精品| 国产精品.久久久| 波野结衣二区三区在线| 国产成人精品久久二区二区91 | 美女中出高潮动态图| 99精国产麻豆久久婷婷| 最新的欧美精品一区二区| 一区二区av电影网| 伦理电影免费视频| 亚洲 欧美一区二区三区| 2021少妇久久久久久久久久久| 999精品在线视频| 一级爰片在线观看| 伦精品一区二区三区| 亚洲第一区二区三区不卡| 一区在线观看完整版| 国产激情久久老熟女| 美女福利国产在线| 欧美日韩精品网址| 女人被躁到高潮嗷嗷叫费观| 午夜福利乱码中文字幕| 90打野战视频偷拍视频| 2022亚洲国产成人精品| 日韩一卡2卡3卡4卡2021年| 国产精品女同一区二区软件| 中文字幕精品免费在线观看视频| 亚洲精品第二区| 久久精品国产a三级三级三级| 亚洲精华国产精华液的使用体验| 丝袜美腿诱惑在线| 久久99蜜桃精品久久| 亚洲av福利一区| 好男人视频免费观看在线| 亚洲欧美日韩另类电影网站| 自线自在国产av| 少妇精品久久久久久久| 大香蕉久久网| 丝袜人妻中文字幕| 中文字幕色久视频| 欧美日韩视频高清一区二区三区二| 欧美 亚洲 国产 日韩一| 婷婷成人精品国产| 国产成人午夜福利电影在线观看| 国产av精品麻豆| 一区在线观看完整版| 男女边吃奶边做爰视频| 精品国产国语对白av| 一边摸一边做爽爽视频免费| 久久婷婷青草| 久久久久久久久久久免费av| 制服诱惑二区| 国产男女内射视频| 卡戴珊不雅视频在线播放| 麻豆精品久久久久久蜜桃| 亚洲欧美色中文字幕在线| 色94色欧美一区二区| 日韩不卡一区二区三区视频在线| 日日爽夜夜爽网站| 女人被躁到高潮嗷嗷叫费观| 久久毛片免费看一区二区三区| 少妇熟女欧美另类| a级片在线免费高清观看视频| 婷婷色综合大香蕉| av免费观看日本| 午夜久久久在线观看| 国产精品熟女久久久久浪| 亚洲国产看品久久| 捣出白浆h1v1| 国产成人91sexporn| 最近的中文字幕免费完整| 国产探花极品一区二区| 2022亚洲国产成人精品| 成人漫画全彩无遮挡| 国产成人欧美| 国产精品.久久久| 亚洲欧洲国产日韩| 亚洲第一区二区三区不卡| 中文天堂在线官网| av免费在线看不卡| 美女福利国产在线| 久久精品久久久久久久性| 精品一品国产午夜福利视频| 看免费成人av毛片| 久久精品亚洲av国产电影网| 国产成人精品婷婷| videossex国产| 天天躁夜夜躁狠狠躁躁| 美女高潮到喷水免费观看| 街头女战士在线观看网站| 男女边吃奶边做爰视频| 99久国产av精品国产电影| 国产在线一区二区三区精| 亚洲中文av在线| 免费在线观看黄色视频的| 最近中文字幕高清免费大全6| 香蕉精品网在线| 国产在线视频一区二区| 两性夫妻黄色片| 青春草亚洲视频在线观看| 亚洲国产色片| 伊人久久国产一区二区| 亚洲,欧美精品.| 国产成人精品福利久久| 制服诱惑二区| 中国国产av一级| 少妇人妻精品综合一区二区| 久热这里只有精品99| 亚洲av欧美aⅴ国产| 一边亲一边摸免费视频| 国产精品免费视频内射| 亚洲视频免费观看视频| 最新的欧美精品一区二区| 一区在线观看完整版| 制服诱惑二区| 天美传媒精品一区二区| 香蕉精品网在线| 五月开心婷婷网| 少妇人妻 视频| 国产片特级美女逼逼视频| 成年女人在线观看亚洲视频| 天天操日日干夜夜撸| 亚洲精品国产一区二区精华液| 国产精品久久久久成人av| 国产一区有黄有色的免费视频| 色94色欧美一区二区| av片东京热男人的天堂| 一区二区三区乱码不卡18| av在线app专区| 最近最新中文字幕大全免费视频 | 秋霞在线观看毛片| 黄片小视频在线播放| 国产精品成人在线| 国产综合精华液| 狠狠精品人妻久久久久久综合| 国产精品.久久久| 一区二区av电影网| 春色校园在线视频观看| 在线观看免费视频网站a站| 亚洲av中文av极速乱| 两个人看的免费小视频| 少妇熟女欧美另类| 日韩,欧美,国产一区二区三区| 看非洲黑人一级黄片| 亚洲欧美一区二区三区久久| 国产成人精品一,二区| 纵有疾风起免费观看全集完整版| √禁漫天堂资源中文www| 在线免费观看不下载黄p国产| 午夜福利,免费看| 免费在线观看黄色视频的| 国产日韩欧美视频二区| 日本黄色日本黄色录像| 成人国语在线视频| 免费日韩欧美在线观看| 久久这里有精品视频免费| 91午夜精品亚洲一区二区三区| 欧美中文综合在线视频| 9热在线视频观看99| 宅男免费午夜| 国产熟女午夜一区二区三区| 国产成人a∨麻豆精品| 国产精品久久久av美女十八| 在线天堂最新版资源| 国产国语露脸激情在线看| videossex国产| 午夜福利乱码中文字幕| 99香蕉大伊视频| 久久精品国产综合久久久| 国产免费又黄又爽又色| 久久99精品国语久久久| 黄色毛片三级朝国网站| 天天影视国产精品| 2022亚洲国产成人精品| 高清不卡的av网站| 在线观看一区二区三区激情| 男女边摸边吃奶| 欧美中文综合在线视频| 蜜桃国产av成人99| 波野结衣二区三区在线| 亚洲av福利一区| 少妇人妻精品综合一区二区| 国产伦理片在线播放av一区| 国产精品亚洲av一区麻豆 | 国产1区2区3区精品| 免费黄频网站在线观看国产| a级毛片在线看网站| 黑人欧美特级aaaaaa片| 建设人人有责人人尽责人人享有的| 久久av网站| 只有这里有精品99| 啦啦啦视频在线资源免费观看| 久久人人97超碰香蕉20202| 熟女电影av网| 国产精品三级大全| 老汉色∧v一级毛片| 丝袜美腿诱惑在线| 国产激情久久老熟女| 亚洲五月色婷婷综合| 日韩,欧美,国产一区二区三区| 久久久久久伊人网av| 国产精品99久久99久久久不卡 | 精品人妻在线不人妻| 如日韩欧美国产精品一区二区三区| 国精品久久久久久国模美| 精品亚洲成a人片在线观看| 天天躁夜夜躁狠狠久久av| 久久青草综合色| 国产色婷婷99| 日产精品乱码卡一卡2卡三| 亚洲色图综合在线观看| 免费少妇av软件| 日韩人妻精品一区2区三区| 免费观看性生交大片5| 午夜影院在线不卡| 欧美日韩一区二区视频在线观看视频在线| 国产精品偷伦视频观看了| 天天躁夜夜躁狠狠久久av| 亚洲美女黄色视频免费看| 在线观看免费日韩欧美大片| 日产精品乱码卡一卡2卡三| 亚洲精品aⅴ在线观看| a级毛片在线看网站| 夫妻午夜视频| 欧美精品亚洲一区二区| 欧美黄色片欧美黄色片| 久久精品久久久久久噜噜老黄| a级毛片黄视频| 国产精品.久久久| 大片电影免费在线观看免费| 日韩 亚洲 欧美在线| 丁香六月天网| 男的添女的下面高潮视频| 欧美日韩视频高清一区二区三区二| 熟妇人妻不卡中文字幕| 一级,二级,三级黄色视频| 男男h啪啪无遮挡| h视频一区二区三区| 五月伊人婷婷丁香| 水蜜桃什么品种好| 女性生殖器流出的白浆| 亚洲欧洲精品一区二区精品久久久 | 午夜影院在线不卡| 欧美日韩一级在线毛片| 99久久精品国产国产毛片| 欧美在线黄色| 国产精品一国产av| av不卡在线播放| 久久久欧美国产精品| 精品久久久久久电影网| 国产乱来视频区| 90打野战视频偷拍视频| 国产亚洲欧美精品永久| 999久久久国产精品视频| 国产成人精品久久久久久| 国产精品 欧美亚洲| 久久久久国产精品人妻一区二区| 亚洲中文av在线| 久久影院123| tube8黄色片| 日韩中字成人| 国产精品久久久久久精品电影小说| 国产在线视频一区二区| 大片电影免费在线观看免费| 欧美成人精品欧美一级黄| 永久免费av网站大全| 欧美中文综合在线视频| 亚洲国产成人一精品久久久| 交换朋友夫妻互换小说| 国产成人免费无遮挡视频| 1024香蕉在线观看| 少妇人妻久久综合中文| 午夜激情av网站| 久久国产精品大桥未久av| 国产精品一区二区在线不卡| 97在线人人人人妻| 国产精品无大码| 亚洲欧美一区二区三区国产| 欧美日韩一区二区视频在线观看视频在线| 黄色怎么调成土黄色| 成年动漫av网址| 日韩精品有码人妻一区| 少妇 在线观看| 亚洲第一区二区三区不卡| 日日爽夜夜爽网站| 国产午夜精品一二区理论片| 国产精品亚洲av一区麻豆 | 水蜜桃什么品种好| 老司机影院毛片| 最近中文字幕高清免费大全6| 精品一品国产午夜福利视频| 91aial.com中文字幕在线观看| 丝袜脚勾引网站| 欧美最新免费一区二区三区| 大码成人一级视频| 一级毛片我不卡| 不卡视频在线观看欧美| 啦啦啦在线观看免费高清www| 九九爱精品视频在线观看| 婷婷色麻豆天堂久久| 亚洲精品国产av成人精品| 捣出白浆h1v1| 欧美日韩亚洲国产一区二区在线观看 | 精品午夜福利在线看| 麻豆精品久久久久久蜜桃| 午夜福利在线观看免费完整高清在| 啦啦啦在线免费观看视频4| 赤兔流量卡办理| 精品久久蜜臀av无| 天天躁夜夜躁狠狠躁躁| 视频在线观看一区二区三区| 中文字幕精品免费在线观看视频| 亚洲精品美女久久久久99蜜臀 | 亚洲人成网站在线观看播放| 可以免费在线观看a视频的电影网站 | 爱豆传媒免费全集在线观看| 晚上一个人看的免费电影| 欧美精品一区二区大全| 咕卡用的链子| 欧美日韩综合久久久久久| 国产成人免费无遮挡视频| 18+在线观看网站| 中文精品一卡2卡3卡4更新| 高清在线视频一区二区三区| 免费人妻精品一区二区三区视频| 亚洲国产色片| 肉色欧美久久久久久久蜜桃| 中文字幕av电影在线播放| 18禁国产床啪视频网站| 大片电影免费在线观看免费| 久久人人97超碰香蕉20202| 波多野结衣一区麻豆| 麻豆乱淫一区二区| 国产女主播在线喷水免费视频网站| 美女福利国产在线| 免费在线观看完整版高清| 国产片内射在线| 国产日韩欧美在线精品| 久久久久视频综合| 天天躁夜夜躁狠狠躁躁| 欧美精品一区二区免费开放| 欧美成人午夜精品| 中文天堂在线官网| 春色校园在线视频观看| 你懂的网址亚洲精品在线观看| 卡戴珊不雅视频在线播放| 亚洲中文av在线| 色视频在线一区二区三区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久97久久精品| www.av在线官网国产| 亚洲精品av麻豆狂野| 国产毛片在线视频| 精品第一国产精品| 国产成人午夜福利电影在线观看| √禁漫天堂资源中文www| 成人午夜精彩视频在线观看| 大码成人一级视频| 欧美精品高潮呻吟av久久| 日韩中文字幕视频在线看片| 丝袜喷水一区| 久久久国产一区二区| 日韩欧美一区视频在线观看| 午夜久久久在线观看| 免费看不卡的av| 搡老乐熟女国产| 久久亚洲国产成人精品v| 免费av中文字幕在线| 另类亚洲欧美激情| 色哟哟·www| 欧美变态另类bdsm刘玥| 一级毛片 在线播放| 韩国高清视频一区二区三区| 欧美精品人与动牲交sv欧美| 在线 av 中文字幕| 免费女性裸体啪啪无遮挡网站| 亚洲欧美清纯卡通| 久久久久精品人妻al黑| 七月丁香在线播放| 精品亚洲成a人片在线观看| 又大又黄又爽视频免费| 亚洲精品日本国产第一区| 国产精品一区二区在线观看99| 91在线精品国自产拍蜜月| 看免费成人av毛片| xxxhd国产人妻xxx| 精品少妇黑人巨大在线播放| 国产成人av激情在线播放| 国产一区二区三区综合在线观看| 亚洲欧洲精品一区二区精品久久久 | 在线观看免费日韩欧美大片| 97精品久久久久久久久久精品| 欧美精品人与动牲交sv欧美| 成年人午夜在线观看视频| 纵有疾风起免费观看全集完整版| 美女国产高潮福利片在线看| 永久网站在线| 夜夜骑夜夜射夜夜干| 久久人妻熟女aⅴ| 色婷婷久久久亚洲欧美| 精品人妻偷拍中文字幕| 午夜激情av网站| 久久久久久久久久久久大奶| 欧美日韩精品网址| 超色免费av| 免费看av在线观看网站| 成人手机av| 欧美亚洲日本最大视频资源| 中文字幕人妻熟女乱码| 女人被躁到高潮嗷嗷叫费观| 一区二区av电影网| 日日啪夜夜爽| 又粗又硬又长又爽又黄的视频| 日韩一区二区视频免费看| 黑人猛操日本美女一级片| 欧美老熟妇乱子伦牲交| 日本av手机在线免费观看| 伦精品一区二区三区| 少妇猛男粗大的猛烈进出视频| 免费高清在线观看日韩| 99久久精品国产国产毛片| 成人毛片60女人毛片免费| 亚洲欧美一区二区三区国产| 亚洲精品日韩在线中文字幕| 建设人人有责人人尽责人人享有的| 成人黄色视频免费在线看| 妹子高潮喷水视频| 日韩伦理黄色片| 新久久久久国产一级毛片| 自拍欧美九色日韩亚洲蝌蚪91| 久久久精品免费免费高清| 丰满迷人的少妇在线观看| 精品酒店卫生间| 午夜福利影视在线免费观看| 天天操日日干夜夜撸| 欧美成人午夜精品| 欧美日韩一级在线毛片| 国产亚洲av片在线观看秒播厂| 自拍欧美九色日韩亚洲蝌蚪91| 1024香蕉在线观看| 中文字幕av电影在线播放| 大陆偷拍与自拍| 久久久久网色| 成人手机av| 日韩精品青青久久久久久|