• <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真菌對棉花的侵染率及氮磷含量的影響
    老司机亚洲免费影院| 日本五十路高清| 岛国在线观看网站| 黄色 视频免费看| 亚洲成人免费电影在线观看| 日韩 欧美 亚洲 中文字幕| 日韩三级视频一区二区三区| 成人黄色视频免费在线看| 法律面前人人平等表现在哪些方面| 在线观看免费高清a一片| 韩国av一区二区三区四区| 欧美亚洲日本最大视频资源| 一二三四社区在线视频社区8| 精品福利永久在线观看| √禁漫天堂资源中文www| av天堂久久9| 在线观看免费高清a一片| 欧美 亚洲 国产 日韩一| 成人影院久久| 99热国产这里只有精品6| 深夜精品福利| а√天堂www在线а√下载 | 国产精品香港三级国产av潘金莲| 1024香蕉在线观看| av有码第一页| av有码第一页| 国产男女超爽视频在线观看| 中文字幕人妻丝袜制服| 久久九九热精品免费| 午夜成年电影在线免费观看| 色综合婷婷激情| 超碰97精品在线观看| 亚洲第一av免费看| 国产av精品麻豆| 久久精品亚洲av国产电影网| 狠狠婷婷综合久久久久久88av| 熟女少妇亚洲综合色aaa.| xxxhd国产人妻xxx| 99re在线观看精品视频| 久久国产亚洲av麻豆专区| 国产高清国产精品国产三级| 国产成人啪精品午夜网站| 在线播放国产精品三级| 欧美激情高清一区二区三区| 午夜影院日韩av| 精品电影一区二区在线| 久久中文看片网| 日韩人妻精品一区2区三区| 欧美人与性动交α欧美精品济南到| 欧美国产精品一级二级三级| 黑人巨大精品欧美一区二区mp4| 亚洲少妇的诱惑av| 一个人免费在线观看的高清视频| 美女 人体艺术 gogo| 国产成人av激情在线播放| 国产精品一区二区精品视频观看| 国产精品久久久av美女十八| 日本a在线网址| aaaaa片日本免费| 一级毛片高清免费大全| 久久国产精品影院| 黄色 视频免费看| 欧美黑人精品巨大| 麻豆成人av在线观看| 欧美性长视频在线观看| 亚洲av日韩精品久久久久久密| 91av网站免费观看| 91av网站免费观看| 免费高清在线观看日韩| av电影中文网址| 欧美色视频一区免费| 欧美在线黄色| 午夜两性在线视频| 精品一区二区三区四区五区乱码| e午夜精品久久久久久久| 久99久视频精品免费| 精品久久久久久久毛片微露脸| 欧美日韩黄片免| 中文字幕最新亚洲高清| 18在线观看网站| 午夜福利视频在线观看免费| 午夜两性在线视频| 一级黄色大片毛片| 黄色丝袜av网址大全| 99re6热这里在线精品视频| 亚洲av欧美aⅴ国产| 婷婷丁香在线五月| 亚洲欧美日韩另类电影网站| 美女高潮喷水抽搐中文字幕| 欧美日韩一级在线毛片| 看黄色毛片网站| 国产精品国产高清国产av | 亚洲美女黄片视频| 午夜福利欧美成人| 少妇被粗大的猛进出69影院| 成人手机av| 成人三级做爰电影| 制服人妻中文乱码| 美女福利国产在线| 俄罗斯特黄特色一大片| 美女高潮喷水抽搐中文字幕| 男人舔女人的私密视频| 午夜两性在线视频| 国产成+人综合+亚洲专区| 久久国产精品人妻蜜桃| 久久久久视频综合| 国产精品秋霞免费鲁丝片| 丝袜美腿诱惑在线| 99热网站在线观看| 老熟妇乱子伦视频在线观看| 男女床上黄色一级片免费看| 男人的好看免费观看在线视频 | 国产免费男女视频| 亚洲精品一二三| 欧美日韩成人在线一区二区| 国产1区2区3区精品| 黑人巨大精品欧美一区二区蜜桃| 欧美激情高清一区二区三区| 成年女人毛片免费观看观看9 | 国产免费现黄频在线看| 老司机福利观看| 激情在线观看视频在线高清 | 日日夜夜操网爽| 三上悠亚av全集在线观看| 丰满的人妻完整版| 久久午夜亚洲精品久久| 亚洲成人手机| 侵犯人妻中文字幕一二三四区| 国产91精品成人一区二区三区| 亚洲五月婷婷丁香| 丁香六月欧美| 国产成人系列免费观看| av超薄肉色丝袜交足视频| √禁漫天堂资源中文www| 国产精品一区二区精品视频观看| 在线天堂中文资源库| 午夜精品在线福利| 看免费av毛片| 天天影视国产精品| 成人手机av| 免费一级毛片在线播放高清视频 | 12—13女人毛片做爰片一| 黄色视频不卡| 大码成人一级视频| 精品一区二区三区视频在线观看免费 | 免费在线观看日本一区| 精品亚洲成国产av| 一边摸一边抽搐一进一小说 | 欧美日韩亚洲综合一区二区三区_| 一区在线观看完整版| 国产精品电影一区二区三区 | 午夜成年电影在线免费观看| 中文亚洲av片在线观看爽 | 高清欧美精品videossex| 在线观看一区二区三区激情| 亚洲人成电影观看| 王馨瑶露胸无遮挡在线观看| 久久精品国产亚洲av香蕉五月 | 乱人伦中国视频| 中文字幕人妻丝袜制服| 中文字幕最新亚洲高清| 黄色成人免费大全| 国产免费av片在线观看野外av| 国产av精品麻豆| 成年人黄色毛片网站| 热re99久久国产66热| 老熟女久久久| a在线观看视频网站| 脱女人内裤的视频| 免费黄频网站在线观看国产| 欧美成狂野欧美在线观看| 亚洲精品av麻豆狂野| 婷婷精品国产亚洲av在线 | 18禁观看日本| 久久精品成人免费网站| 亚洲成人手机| 亚洲男人天堂网一区| 国产麻豆69| 一边摸一边抽搐一进一出视频| 欧美黄色片欧美黄色片| 亚洲第一欧美日韩一区二区三区| 国产欧美日韩一区二区三区在线| 夜夜爽天天搞| 精品高清国产在线一区| 午夜福利一区二区在线看| 成人手机av| 午夜福利影视在线免费观看| 国产精品久久久人人做人人爽| 十八禁人妻一区二区| 亚洲国产欧美网| 一区二区日韩欧美中文字幕| 欧美日韩一级在线毛片| 亚洲精品成人av观看孕妇| 日韩三级视频一区二区三区| 777米奇影视久久| 999精品在线视频| 成人亚洲精品一区在线观看| svipshipincom国产片| 纯流量卡能插随身wifi吗| 一二三四社区在线视频社区8| 欧美黄色淫秽网站| 18禁观看日本| 青草久久国产| 国产乱人伦免费视频| 国产一区二区三区综合在线观看| 中出人妻视频一区二区| 黄色 视频免费看| 大型av网站在线播放| 50天的宝宝边吃奶边哭怎么回事| 精品人妻1区二区| 久久久久国内视频| 免费日韩欧美在线观看| 亚洲中文av在线| 美女扒开内裤让男人捅视频| 免费看a级黄色片| 大香蕉久久成人网| 欧美人与性动交α欧美软件| 国产免费现黄频在线看| 天堂√8在线中文| 99久久综合精品五月天人人| 老汉色∧v一级毛片| 免费观看人在逋| 9色porny在线观看| www.自偷自拍.com| 91国产中文字幕| 不卡av一区二区三区| 亚洲男人天堂网一区| 777米奇影视久久| 国产高清视频在线播放一区| 久久精品91无色码中文字幕| 国内久久婷婷六月综合欲色啪| 好男人电影高清在线观看| 精品福利永久在线观看| 午夜福利影视在线免费观看| 19禁男女啪啪无遮挡网站| 欧美日韩亚洲国产一区二区在线观看 | 午夜福利乱码中文字幕| 99久久国产精品久久久| 国产无遮挡羞羞视频在线观看| 51午夜福利影视在线观看| 亚洲色图综合在线观看| 成年动漫av网址| 视频区图区小说| 国产亚洲av高清不卡| 亚洲精品国产精品久久久不卡| 在线av久久热| 一级黄色大片毛片| 99在线人妻在线中文字幕 | 超碰成人久久| 女警被强在线播放| 中国美女看黄片| 欧美一级毛片孕妇| 丰满迷人的少妇在线观看| 亚洲一区二区三区不卡视频| 国产免费av片在线观看野外av| 欧美色视频一区免费| 亚洲一区二区三区不卡视频| 国产高清激情床上av| 久久久国产成人免费| 欧美精品啪啪一区二区三区| 亚洲精品久久成人aⅴ小说| 狠狠婷婷综合久久久久久88av| 久久久国产一区二区| 九色亚洲精品在线播放| www.自偷自拍.com| 欧美日韩亚洲综合一区二区三区_| 大码成人一级视频| 国产精品1区2区在线观看. | 成人特级黄色片久久久久久久| 国产深夜福利视频在线观看| 亚洲av成人av| 高清毛片免费观看视频网站 | 国产成人精品无人区| 亚洲精品国产一区二区精华液| 在线观看www视频免费| 90打野战视频偷拍视频| 18禁美女被吸乳视频| av一本久久久久| 妹子高潮喷水视频| 极品教师在线免费播放| 法律面前人人平等表现在哪些方面| 超碰97精品在线观看| 男女之事视频高清在线观看| 黄片播放在线免费| 久久久久久人人人人人| 精品国产一区二区三区久久久樱花| 亚洲人成伊人成综合网2020| 欧美 亚洲 国产 日韩一| 久久精品国产a三级三级三级| 欧美精品人与动牲交sv欧美| 亚洲精品久久午夜乱码| av线在线观看网站| 精品福利永久在线观看| 丰满迷人的少妇在线观看| ponron亚洲| 国产在线一区二区三区精| 欧美精品人与动牲交sv欧美| 最近最新中文字幕大全电影3 | 村上凉子中文字幕在线| 日韩视频一区二区在线观看| 成人国语在线视频| 变态另类成人亚洲欧美熟女 | 精品少妇一区二区三区视频日本电影| 国产三级黄色录像| 国产精品永久免费网站| 国产成人欧美| 男女之事视频高清在线观看| 久久午夜综合久久蜜桃| 日韩免费av在线播放| 大香蕉久久网| 多毛熟女@视频| 午夜免费成人在线视频| 欧美亚洲 丝袜 人妻 在线| 久久青草综合色| 一二三四在线观看免费中文在| 国产99白浆流出| 久久久国产精品麻豆| 国产激情久久老熟女| 亚洲熟女毛片儿| 99在线人妻在线中文字幕 | 国产野战对白在线观看| 久久久久久久久久久久大奶| 夜夜躁狠狠躁天天躁| 国产深夜福利视频在线观看| 亚洲av欧美aⅴ国产| av天堂久久9| 纯流量卡能插随身wifi吗| 日韩人妻精品一区2区三区| 亚洲色图综合在线观看| 久久亚洲精品不卡| 国产精品美女特级片免费视频播放器 | 国产亚洲欧美精品永久| 国产亚洲精品久久久久久毛片 | 午夜精品久久久久久毛片777| 亚洲情色 制服丝袜| 老司机深夜福利视频在线观看| 午夜福利,免费看| 久久精品国产99精品国产亚洲性色 | 国产成人欧美| 99国产精品一区二区蜜桃av | 黄片小视频在线播放| 久久久久久久久久久久大奶| 少妇被粗大的猛进出69影院| 叶爱在线成人免费视频播放| 村上凉子中文字幕在线| 露出奶头的视频| 老司机福利观看| 国产99白浆流出| cao死你这个sao货| 精品国产国语对白av| 人妻一区二区av| 亚洲国产中文字幕在线视频| 91字幕亚洲| 男女免费视频国产| 9色porny在线观看| av线在线观看网站| 在线观看午夜福利视频| 精品国内亚洲2022精品成人 | 亚洲av美国av| 日本wwww免费看| 不卡一级毛片| 最近最新免费中文字幕在线| 天天躁日日躁夜夜躁夜夜| 99在线人妻在线中文字幕 | 亚洲精品粉嫩美女一区| 99热网站在线观看| 亚洲精品国产区一区二| 最近最新中文字幕大全免费视频| 看免费av毛片| 伊人久久大香线蕉亚洲五| 极品教师在线免费播放| 可以免费在线观看a视频的电影网站| 很黄的视频免费| 色播在线永久视频| 1024香蕉在线观看| 99热国产这里只有精品6| 欧美乱码精品一区二区三区| 国产免费av片在线观看野外av| 操出白浆在线播放| 成年版毛片免费区| 91老司机精品| 老熟女久久久| av超薄肉色丝袜交足视频| 亚洲精品av麻豆狂野| 一级a爱片免费观看的视频| 国产成人影院久久av| 高潮久久久久久久久久久不卡| 一区二区三区国产精品乱码| 久热爱精品视频在线9| 日本欧美视频一区| 免费观看a级毛片全部| av免费在线观看网站| 午夜福利免费观看在线| 久久久久久亚洲精品国产蜜桃av| 亚洲专区字幕在线| 高清视频免费观看一区二区| 日本撒尿小便嘘嘘汇集6| 又黄又爽又免费观看的视频| 国产一区二区激情短视频| 夫妻午夜视频| 涩涩av久久男人的天堂| 久久99一区二区三区| 美女国产高潮福利片在线看| 久久精品91无色码中文字幕| 精品国产亚洲在线| 一进一出抽搐动态| 天天躁狠狠躁夜夜躁狠狠躁| 99久久国产精品久久久| 亚洲欧美色中文字幕在线| 十分钟在线观看高清视频www| 涩涩av久久男人的天堂| 美女扒开内裤让男人捅视频| 欧美激情 高清一区二区三区| 夜夜躁狠狠躁天天躁| 久久香蕉国产精品| 午夜福利影视在线免费观看| av网站免费在线观看视频| 亚洲精品成人av观看孕妇| 成人手机av| 国产欧美日韩一区二区三区在线| 国产成人免费观看mmmm| 涩涩av久久男人的天堂| 亚洲国产中文字幕在线视频| 美女高潮喷水抽搐中文字幕| 丝袜美足系列| 精品免费久久久久久久清纯 | 亚洲国产精品一区二区三区在线| 麻豆成人av在线观看| 中文欧美无线码| 久久久久久久国产电影| 国产成人av教育| 精品久久久久久,| av片东京热男人的天堂| 国产主播在线观看一区二区| 我的亚洲天堂| 国产成人精品在线电影| 美女国产高潮福利片在线看| 久久精品亚洲av国产电影网| 黑人巨大精品欧美一区二区蜜桃| 欧美成人免费av一区二区三区 | 国产单亲对白刺激| 日韩欧美一区二区三区在线观看 | 在线视频色国产色| 亚洲avbb在线观看| 91精品三级在线观看| 热re99久久国产66热| 久久人妻熟女aⅴ| 国产精品免费视频内射| 男人的好看免费观看在线视频 | 99精国产麻豆久久婷婷| 免费一级毛片在线播放高清视频 | 人人妻人人爽人人添夜夜欢视频| 我的亚洲天堂| 久久中文字幕一级| 女人爽到高潮嗷嗷叫在线视频| 精品国产一区二区三区久久久樱花| www日本在线高清视频| 男女下面插进去视频免费观看| 国产精品98久久久久久宅男小说| 麻豆国产av国片精品| 亚洲av成人不卡在线观看播放网| 国产亚洲欧美精品永久| 国产精品二区激情视频| 脱女人内裤的视频| 黄频高清免费视频| 亚洲欧美日韩另类电影网站| 18禁黄网站禁片午夜丰满| 搡老岳熟女国产| 久久久久精品国产欧美久久久| 亚洲精品美女久久av网站| av电影中文网址| 色精品久久人妻99蜜桃| 天天躁夜夜躁狠狠躁躁| 欧美日韩亚洲国产一区二区在线观看 | 中亚洲国语对白在线视频| 超碰97精品在线观看| 久久香蕉激情| 一进一出抽搐gif免费好疼 | 久久亚洲真实| 欧美 亚洲 国产 日韩一| 国产av精品麻豆| 在线国产一区二区在线| 精品人妻熟女毛片av久久网站| 久久影院123| 久久精品熟女亚洲av麻豆精品| 99精品久久久久人妻精品| 在线十欧美十亚洲十日本专区| 亚洲人成77777在线视频| 欧美精品av麻豆av| 一区福利在线观看| 老司机午夜福利在线观看视频| 久久久国产欧美日韩av| 久久久久久亚洲精品国产蜜桃av| 午夜免费成人在线视频| 精品一区二区三区四区五区乱码| www.自偷自拍.com| 欧美激情极品国产一区二区三区| 午夜免费鲁丝| 18禁国产床啪视频网站| 天堂动漫精品| 黄色丝袜av网址大全| 黑人巨大精品欧美一区二区蜜桃| 午夜老司机福利片| av天堂在线播放| 女人精品久久久久毛片| 亚洲少妇的诱惑av| 亚洲人成伊人成综合网2020| 丝袜在线中文字幕| 大香蕉久久成人网| 成人特级黄色片久久久久久久| 在线看a的网站| 精品久久久精品久久久| 亚洲人成77777在线视频| 国产成人精品久久二区二区免费| 午夜福利免费观看在线| 精品国产国语对白av| 极品教师在线免费播放| 午夜福利影视在线免费观看| 欧美日韩乱码在线| av电影中文网址| 久久精品国产亚洲av高清一级| www日本在线高清视频| 啪啪无遮挡十八禁网站| 757午夜福利合集在线观看| 国产精品99久久99久久久不卡| 欧美黑人欧美精品刺激| 怎么达到女性高潮| 欧美激情极品国产一区二区三区| 国产精品电影一区二区三区 | 亚洲色图 男人天堂 中文字幕| 精品国产超薄肉色丝袜足j| 中国美女看黄片| 国产精品乱码一区二三区的特点 | 久热爱精品视频在线9| 99riav亚洲国产免费| 午夜两性在线视频| 久久午夜亚洲精品久久| 久久久久久人人人人人| 老司机午夜十八禁免费视频| 涩涩av久久男人的天堂| 欧美成人免费av一区二区三区 | 亚洲人成电影观看| 老汉色∧v一级毛片| 无人区码免费观看不卡| 久久精品国产a三级三级三级| 夫妻午夜视频| 久久国产精品大桥未久av| 国产亚洲精品久久久久久毛片 | 欧美老熟妇乱子伦牲交| 9191精品国产免费久久| 18在线观看网站| 国产亚洲精品第一综合不卡| 91老司机精品| 99国产精品一区二区蜜桃av | 欧美另类亚洲清纯唯美| 欧美日韩中文字幕国产精品一区二区三区 | 嫁个100分男人电影在线观看| 久久久久久人人人人人| 成人18禁在线播放| 天天添夜夜摸| 国产无遮挡羞羞视频在线观看| 精品午夜福利视频在线观看一区| 99re6热这里在线精品视频| 亚洲国产精品合色在线| 黄色片一级片一级黄色片| 他把我摸到了高潮在线观看| 一本综合久久免费| 国产熟女午夜一区二区三区| 国产成人系列免费观看| 亚洲精品av麻豆狂野| 一级,二级,三级黄色视频| 亚洲伊人色综图| 制服人妻中文乱码| 日本a在线网址| 精品电影一区二区在线| 国产人伦9x9x在线观看| 天堂中文最新版在线下载| ponron亚洲| 伦理电影免费视频| 久久人妻熟女aⅴ| 国产欧美日韩一区二区三区在线| svipshipincom国产片| 丝袜美腿诱惑在线| 12—13女人毛片做爰片一| 免费av中文字幕在线| 亚洲九九香蕉| 女人被躁到高潮嗷嗷叫费观| 成熟少妇高潮喷水视频| 免费高清在线观看日韩| 超碰97精品在线观看| www.精华液| 丝袜美腿诱惑在线| 久久久久久人人人人人| 正在播放国产对白刺激| 18禁美女被吸乳视频| 伊人久久大香线蕉亚洲五| 最近最新中文字幕大全电影3 | 香蕉久久夜色| 69精品国产乱码久久久| 他把我摸到了高潮在线观看| 欧美人与性动交α欧美软件| 一进一出抽搐动态| 亚洲,欧美精品.| 美女国产高潮福利片在线看| aaaaa片日本免费| 国产欧美日韩精品亚洲av| 老司机在亚洲福利影院| 女人精品久久久久毛片| 国产午夜精品久久久久久| 亚洲欧美激情在线| 久久久国产一区二区| 极品教师在线免费播放| 欧美日韩亚洲高清精品| 黄色丝袜av网址大全| 亚洲精品中文字幕一二三四区|