• <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真菌對棉花的侵染率及氮磷含量的影響
    国产精华一区二区三区| 午夜福利在线观看免费完整高清在 | 长腿黑丝高跟| 中国美白少妇内射xxxbb| 最近的中文字幕免费完整| 久久久精品欧美日韩精品| 欧美性感艳星| 亚洲av美国av| 人妻丰满熟妇av一区二区三区| 久久久精品大字幕| 成人一区二区视频在线观看| 精品久久久噜噜| 国产激情偷乱视频一区二区| 午夜免费男女啪啪视频观看 | 男人舔奶头视频| 变态另类成人亚洲欧美熟女| 国产精品亚洲一级av第二区| 69av精品久久久久久| 人妻夜夜爽99麻豆av| 国产精品一二三区在线看| 菩萨蛮人人尽说江南好唐韦庄 | 此物有八面人人有两片| 亚洲精品日韩av片在线观看| 少妇人妻一区二区三区视频| 日韩人妻高清精品专区| 精品熟女少妇av免费看| 两个人的视频大全免费| 国产精品久久久久久久电影| 国产毛片a区久久久久| av在线亚洲专区| 国产精品不卡视频一区二区| 色5月婷婷丁香| 欧美激情国产日韩精品一区| 国产精品亚洲美女久久久| 欧美日本视频| 欧美+日韩+精品| 欧美+日韩+精品| 不卡视频在线观看欧美| 午夜免费男女啪啪视频观看 | 亚洲国产精品成人综合色| 日本黄色片子视频| 日韩欧美在线乱码| 日韩av不卡免费在线播放| 五月玫瑰六月丁香| 日本一二三区视频观看| 国内精品久久久久精免费| .国产精品久久| 日韩中字成人| 欧美日韩乱码在线| 欧美高清性xxxxhd video| 国产免费男女视频| 久久精品夜夜夜夜夜久久蜜豆| 免费大片18禁| 中文字幕精品亚洲无线码一区| 日韩大尺度精品在线看网址| 99久久久亚洲精品蜜臀av| 国产精品人妻久久久影院| 久久午夜福利片| 校园春色视频在线观看| 国产精品1区2区在线观看.| 3wmmmm亚洲av在线观看| 亚洲aⅴ乱码一区二区在线播放| 老司机午夜福利在线观看视频| 亚洲国产日韩欧美精品在线观看| 日日啪夜夜撸| 欧美性猛交黑人性爽| 日韩 亚洲 欧美在线| 大又大粗又爽又黄少妇毛片口| 亚洲第一区二区三区不卡| 我的女老师完整版在线观看| 麻豆成人午夜福利视频| 亚洲中文日韩欧美视频| 自拍偷自拍亚洲精品老妇| 亚洲性夜色夜夜综合| 熟女人妻精品中文字幕| 色5月婷婷丁香| 亚洲美女黄片视频| 免费不卡的大黄色大毛片视频在线观看 | 黄色欧美视频在线观看| 国产免费一级a男人的天堂| 亚洲av免费在线观看| 亚洲精品国产av成人精品 | 国产黄a三级三级三级人| 麻豆国产av国片精品| 国产高清三级在线| 亚洲性久久影院| 精品国内亚洲2022精品成人| 又黄又爽又免费观看的视频| 国产黄色视频一区二区在线观看 | 国产亚洲欧美98| 综合色av麻豆| 久久婷婷人人爽人人干人人爱| 1000部很黄的大片| 国产乱人偷精品视频| 日韩亚洲欧美综合| 国产精品日韩av在线免费观看| 欧美另类亚洲清纯唯美| 亚洲高清免费不卡视频| 给我免费播放毛片高清在线观看| 国产成人aa在线观看| 一级av片app| 国产在线精品亚洲第一网站| 日韩欧美在线乱码| av在线蜜桃| 美女免费视频网站| 国产精品一区二区三区四区久久| 亚洲av一区综合| 毛片女人毛片| 一级黄片播放器| 晚上一个人看的免费电影| 国内久久婷婷六月综合欲色啪| 赤兔流量卡办理| 女同久久另类99精品国产91| 国产精品久久久久久亚洲av鲁大| 久久99热这里只有精品18| 精华霜和精华液先用哪个| 午夜精品在线福利| 一个人免费在线观看电影| 有码 亚洲区| 中国美女看黄片| 免费在线观看成人毛片| 老师上课跳d突然被开到最大视频| 国产探花极品一区二区| 插逼视频在线观看| 自拍偷自拍亚洲精品老妇| 九色成人免费人妻av| а√天堂www在线а√下载| 熟女人妻精品中文字幕| 又爽又黄a免费视频| 国产探花极品一区二区| 91午夜精品亚洲一区二区三区| 亚洲av.av天堂| 91在线精品国自产拍蜜月| 桃色一区二区三区在线观看| 青春草视频在线免费观看| 亚洲中文字幕日韩| 亚洲精品成人久久久久久| 深爱激情五月婷婷| 成年av动漫网址| 国产精品亚洲一级av第二区| 狂野欧美激情性xxxx在线观看| 欧美中文日本在线观看视频| 你懂的网址亚洲精品在线观看 | 成人av一区二区三区在线看| 男女之事视频高清在线观看| 亚洲在线自拍视频| 色在线成人网| 我要看日韩黄色一级片| 国产日本99.免费观看| 伦精品一区二区三区| 精品久久久久久久久亚洲| 精品免费久久久久久久清纯| 国产精品永久免费网站| 日韩av不卡免费在线播放| 观看美女的网站| 精品日产1卡2卡| 亚洲欧美日韩东京热| 热99在线观看视频| 午夜爱爱视频在线播放| 国产免费男女视频| 欧美精品国产亚洲| 欧美xxxx性猛交bbbb| 中文资源天堂在线| 日日撸夜夜添| 久久精品夜色国产| av福利片在线观看| 乱人视频在线观看| 久久人人爽人人片av| 级片在线观看| 可以在线观看的亚洲视频| 久久欧美精品欧美久久欧美| 国产女主播在线喷水免费视频网站 | 色在线成人网| 精品一区二区三区人妻视频| 插逼视频在线观看| 国产三级中文精品| 少妇高潮的动态图| 国产av麻豆久久久久久久| 免费黄网站久久成人精品| 18禁裸乳无遮挡免费网站照片| 亚洲va在线va天堂va国产| 国产一区二区三区在线臀色熟女| 午夜福利视频1000在线观看| 亚洲不卡免费看| 亚洲最大成人中文| av天堂在线播放| 欧美性猛交╳xxx乱大交人| 人妻少妇偷人精品九色| 国内精品久久久久精免费| 国产熟女欧美一区二区| 日韩制服骚丝袜av| 亚洲av第一区精品v没综合| aaaaa片日本免费| 亚洲av中文字字幕乱码综合| 免费观看在线日韩| 成人三级黄色视频| 亚洲性夜色夜夜综合| 少妇熟女欧美另类| 亚洲成av人片在线播放无| 成人国产麻豆网| 亚洲av免费高清在线观看| 在线国产一区二区在线| 91麻豆精品激情在线观看国产| 精品人妻熟女av久视频| 插逼视频在线观看| 久久6这里有精品| 日韩大尺度精品在线看网址| 国产伦精品一区二区三区视频9| 给我免费播放毛片高清在线观看| 内射极品少妇av片p| 99热网站在线观看| 亚洲第一电影网av| 悠悠久久av| 精品免费久久久久久久清纯| 麻豆一二三区av精品| 少妇丰满av| 亚洲欧美中文字幕日韩二区| 国产精品亚洲一级av第二区| 国产精品一区二区三区四区久久| 一本一本综合久久| 欧美成人免费av一区二区三区| 舔av片在线| 中国美女看黄片| 国模一区二区三区四区视频| 婷婷六月久久综合丁香| 久久久久久久午夜电影| av在线观看视频网站免费| 国产伦在线观看视频一区| 精品久久久久久久人妻蜜臀av| 久久久国产成人精品二区| 亚洲高清免费不卡视频| 热99在线观看视频| 美女免费视频网站| 国产精品福利在线免费观看| 我的老师免费观看完整版| 亚洲最大成人中文| 亚洲国产精品合色在线| 久久久久免费精品人妻一区二区| 精品熟女少妇av免费看| 日韩三级伦理在线观看| 99国产精品一区二区蜜桃av| 亚洲最大成人手机在线| 九九爱精品视频在线观看| 欧美激情国产日韩精品一区| а√天堂www在线а√下载| 国产精品伦人一区二区| 精品少妇黑人巨大在线播放 | 久久久久性生活片| 免费av观看视频| 日产精品乱码卡一卡2卡三| 久久综合国产亚洲精品| 精品午夜福利视频在线观看一区| 国产精品亚洲美女久久久| 黄片wwwwww| 国产精品久久视频播放| 中文字幕人妻熟人妻熟丝袜美| 欧美日韩一区二区视频在线观看视频在线 | 亚洲在线自拍视频| 久久久成人免费电影| 蜜臀久久99精品久久宅男| 美女大奶头视频| 插阴视频在线观看视频| 国产黄片美女视频| 一夜夜www| 日韩大尺度精品在线看网址| 高清日韩中文字幕在线| 干丝袜人妻中文字幕| 欧美日韩一区二区视频在线观看视频在线 | 亚洲欧美日韩无卡精品| 久久久久久久午夜电影| 国产精品无大码| 九九久久精品国产亚洲av麻豆| 久久久久久久久久黄片| 少妇裸体淫交视频免费看高清| 亚洲四区av| 99久久精品一区二区三区| 亚洲精品一区av在线观看| 九九热线精品视视频播放| 亚洲成人精品中文字幕电影| 日韩精品有码人妻一区| 有码 亚洲区| 久久久国产成人精品二区| 伦精品一区二区三区| 国产精品无大码| 成人午夜高清在线视频| 美女高潮的动态| a级一级毛片免费在线观看| 12—13女人毛片做爰片一| 免费不卡的大黄色大毛片视频在线观看 | 国产私拍福利视频在线观看| 俺也久久电影网| 在线观看午夜福利视频| 可以在线观看的亚洲视频| 内地一区二区视频在线| 欧美区成人在线视频| 亚洲电影在线观看av| 日产精品乱码卡一卡2卡三| 国产不卡一卡二| 91麻豆精品激情在线观看国产| 久久久久久大精品| 嫩草影院入口| 亚洲欧美日韩高清专用| 99热全是精品| 亚洲美女搞黄在线观看 | 波多野结衣高清无吗| 国产蜜桃级精品一区二区三区| 99riav亚洲国产免费| 国产91av在线免费观看| 在线观看66精品国产| 国产精品一区二区三区四区久久| 亚洲人与动物交配视频| 又爽又黄a免费视频| 欧美中文日本在线观看视频| 人妻制服诱惑在线中文字幕| 国产伦精品一区二区三区视频9| 国产一区二区三区在线臀色熟女| 国产视频内射| 日韩精品中文字幕看吧| 99热只有精品国产| 日本黄色片子视频| 波野结衣二区三区在线| 亚洲国产欧洲综合997久久,| 男人狂女人下面高潮的视频| 亚洲专区国产一区二区| 少妇的逼好多水| 国产麻豆成人av免费视频| 亚洲av熟女| 亚洲欧美清纯卡通| 高清毛片免费看| 一级黄色大片毛片| 少妇丰满av| 男人和女人高潮做爰伦理| 午夜免费男女啪啪视频观看 | 久久久色成人| 久久久久免费精品人妻一区二区| 国产蜜桃级精品一区二区三区| 亚洲国产精品成人久久小说 | 国产高清不卡午夜福利| 国产精品三级大全| 欧美丝袜亚洲另类| 伦精品一区二区三区| 久久久久久久久大av| 乱系列少妇在线播放| 亚洲国产精品合色在线| 五月玫瑰六月丁香| 极品教师在线视频| 两个人的视频大全免费| 亚洲欧美日韩卡通动漫| 18+在线观看网站| 99热网站在线观看| 久久人人爽人人片av| 免费av观看视频| 日本黄色片子视频| 99久久精品一区二区三区| 2021天堂中文幕一二区在线观| 精品国产三级普通话版| 亚洲无线在线观看| 欧美xxxx性猛交bbbb| 校园春色视频在线观看| 18+在线观看网站| 永久网站在线| 日产精品乱码卡一卡2卡三| 精品午夜福利在线看| 欧美成人a在线观看| 狂野欧美白嫩少妇大欣赏| 天堂影院成人在线观看| 欧美日韩国产亚洲二区| 不卡一级毛片| 欧美+亚洲+日韩+国产| 露出奶头的视频| 中国美女看黄片| 日韩欧美在线乱码| 成人三级黄色视频| 男人和女人高潮做爰伦理| 亚洲最大成人中文| 长腿黑丝高跟| 色av中文字幕| 我的老师免费观看完整版| 欧美日韩乱码在线| 亚洲欧美日韩高清在线视频| 精品久久久久久久久亚洲| 毛片女人毛片| 搡老熟女国产l中国老女人| 高清日韩中文字幕在线| 桃色一区二区三区在线观看| 欧美日韩综合久久久久久| 九九热线精品视视频播放| 久久久久久久久久久丰满| 神马国产精品三级电影在线观看| 国产精品久久久久久av不卡| 国产av麻豆久久久久久久| 国内揄拍国产精品人妻在线| 夜夜夜夜夜久久久久| 午夜福利在线观看免费完整高清在 | 国产成人精品久久久久久| 欧美成人a在线观看| 婷婷精品国产亚洲av| 色播亚洲综合网| 国内精品一区二区在线观看| 亚洲五月天丁香| 亚洲专区国产一区二区| 亚洲最大成人手机在线| 国产高清视频在线播放一区| 成人性生交大片免费视频hd| 亚洲精品在线观看二区| 极品教师在线视频| 久久99热这里只有精品18| 别揉我奶头 嗯啊视频| 99热网站在线观看| 人妻丰满熟妇av一区二区三区| 亚洲熟妇熟女久久| 国产在线男女| 超碰av人人做人人爽久久| 国产精品亚洲美女久久久| 日本a在线网址| 大香蕉久久网| 亚洲高清免费不卡视频| 亚洲欧美日韩东京热| 久久久久久国产a免费观看| 欧洲精品卡2卡3卡4卡5卡区| 国产午夜福利久久久久久| 女人十人毛片免费观看3o分钟| 久久久久性生活片| 国产午夜精品论理片| 少妇丰满av| av女优亚洲男人天堂| 国产精品电影一区二区三区| 身体一侧抽搐| 国产精品三级大全| 51国产日韩欧美| 看片在线看免费视频| 国产三级在线视频| 久久久精品欧美日韩精品| 夜夜夜夜夜久久久久| 国产成人freesex在线 | 美女cb高潮喷水在线观看| 国产又黄又爽又无遮挡在线| 国产高潮美女av| 在线观看美女被高潮喷水网站| 中文字幕免费在线视频6| 亚洲成人久久性| 在线观看午夜福利视频| 成人高潮视频无遮挡免费网站| 亚洲精品粉嫩美女一区| 99久久成人亚洲精品观看| 亚洲精品日韩av片在线观看| a级毛色黄片| 日本免费一区二区三区高清不卡| 男插女下体视频免费在线播放| 黄色一级大片看看| 黄色一级大片看看| 亚洲精华国产精华液的使用体验 | 国产一区二区亚洲精品在线观看| 中文字幕久久专区| 一区福利在线观看| 在线看三级毛片| 12—13女人毛片做爰片一| 午夜爱爱视频在线播放| 亚洲成人久久性| 别揉我奶头~嗯~啊~动态视频| 女同久久另类99精品国产91| 中文字幕人妻熟人妻熟丝袜美| 两个人视频免费观看高清| 精品久久久久久久久久免费视频| 国产私拍福利视频在线观看| 久久久a久久爽久久v久久| 热99在线观看视频| 老司机午夜福利在线观看视频| 美女大奶头视频| 美女高潮的动态| 午夜福利18| 久久久久性生活片| 亚洲欧美精品自产自拍| 性色avwww在线观看| 亚洲经典国产精华液单| 在线观看美女被高潮喷水网站| 久久久久久久久中文| 男人舔奶头视频| 午夜激情福利司机影院| 久久精品人妻少妇| 可以在线观看的亚洲视频| 两个人的视频大全免费| 欧美色视频一区免费| 99久久久亚洲精品蜜臀av| 禁无遮挡网站| 欧美人与善性xxx| 亚洲最大成人av| 亚洲三级黄色毛片| 免费高清视频大片| 老女人水多毛片| 深夜精品福利| 啦啦啦观看免费观看视频高清| 波多野结衣高清作品| 女人被狂操c到高潮| 国产av在哪里看| 三级毛片av免费| 大型黄色视频在线免费观看| 美女 人体艺术 gogo| 国产av在哪里看| 国产精品一二三区在线看| 天天躁夜夜躁狠狠久久av| 国产成人aa在线观看| 日本黄色视频三级网站网址| 国产精品久久久久久av不卡| 午夜日韩欧美国产| 精品久久久久久久久av| 老师上课跳d突然被开到最大视频| 丰满人妻一区二区三区视频av| 亚洲久久久久久中文字幕| 一进一出好大好爽视频| 十八禁国产超污无遮挡网站| 日韩精品有码人妻一区| 内地一区二区视频在线| 国产蜜桃级精品一区二区三区| 国产高清三级在线| 日韩三级伦理在线观看| 国产 一区 欧美 日韩| 国产精品人妻久久久影院| av专区在线播放| 床上黄色一级片| 国产精品野战在线观看| 成人特级av手机在线观看| 中国美白少妇内射xxxbb| 嫩草影院新地址| 人妻制服诱惑在线中文字幕| 最近中文字幕高清免费大全6| 一个人看视频在线观看www免费| 日韩成人av中文字幕在线观看 | 最新中文字幕久久久久| 午夜福利成人在线免费观看| 国内揄拍国产精品人妻在线| 在线看三级毛片| 国产精品久久久久久久电影| 国产淫片久久久久久久久| 国产高清视频在线播放一区| aaaaa片日本免费| 午夜福利18| 校园春色视频在线观看| 一a级毛片在线观看| 亚洲第一电影网av| 久久久久久久久久黄片| 99热这里只有是精品50| 亚洲av美国av| 天堂网av新在线| 亚洲一级一片aⅴ在线观看| av在线老鸭窝| 久久久成人免费电影| 一本精品99久久精品77| 十八禁网站免费在线| 一级毛片aaaaaa免费看小| 亚洲国产欧美人成| 国产亚洲精品综合一区在线观看| 99在线人妻在线中文字幕| 人人妻人人看人人澡| 99九九线精品视频在线观看视频| 国产色爽女视频免费观看| 菩萨蛮人人尽说江南好唐韦庄 | 成人一区二区视频在线观看| 麻豆国产97在线/欧美| 99热精品在线国产| 黄色视频,在线免费观看| 国产免费男女视频| 久久6这里有精品| 亚洲最大成人中文| 久久久久九九精品影院| 高清午夜精品一区二区三区 | 中文字幕人妻熟人妻熟丝袜美| 亚洲一级一片aⅴ在线观看| 亚洲国产精品久久男人天堂| 无遮挡黄片免费观看| 国产亚洲av嫩草精品影院| 国产午夜精品论理片| 国产成人福利小说| 别揉我奶头 嗯啊视频| 亚洲自偷自拍三级| 精品国内亚洲2022精品成人| 亚州av有码| 男女那种视频在线观看| 免费不卡的大黄色大毛片视频在线观看 | 婷婷精品国产亚洲av在线| 国产人妻一区二区三区在| 国产精品久久久久久亚洲av鲁大| 精品久久久久久久末码| 国产精品久久久久久久电影| 国产亚洲av嫩草精品影院| 网址你懂的国产日韩在线| 国产一区二区激情短视频| 亚洲精品乱码久久久v下载方式| 99久久精品热视频| 女人被狂操c到高潮| 国产不卡一卡二| 一本精品99久久精品77| 国产精品日韩av在线免费观看| 亚洲av成人精品一区久久| 国产欧美日韩精品亚洲av| 国产欧美日韩精品一区二区| 欧美一区二区国产精品久久精品| 大又大粗又爽又黄少妇毛片口| 日本免费一区二区三区高清不卡| 国产视频内射| 亚洲av中文av极速乱| 亚洲高清免费不卡视频| 国产精品亚洲一级av第二区| 男女之事视频高清在线观看| 国产视频一区二区在线看| 一级毛片我不卡| 99久久久亚洲精品蜜臀av| 成年免费大片在线观看| 亚洲精品成人久久久久久| 免费黄网站久久成人精品| 亚州av有码| 国产亚洲精品综合一区在线观看| 久久综合国产亚洲精品| 在线观看免费视频日本深夜| 国产人妻一区二区三区在|