999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

渾河沈撫段水體污染季節性變化研究

2017-03-27 02:28:48禹雪迪
水資源保護 2017年2期
關鍵詞:水質污染

劉 強,禹雪迪

(沈陽建筑大學市政與環境工程學院,遼寧 沈陽 110168)

渾河沈撫段水體污染季節性變化研究

劉 強,禹雪迪

(沈陽建筑大學市政與環境工程學院,遼寧 沈陽 110168)

為明確渾河沈撫段水污染特征,提升水環境治理有效性,2015年3月—2016年2月對渾河沈撫段進行水環境調查。根據豐、平、枯3個水期的pH值、電導率、COD、NH3-N、BOD及CODMn等水質指標,運用改進的綜合水質標識指數法識別出主要污染因子,并分析渾河沈撫段水質污染時空變化特征。結果表明,研究區域的pH值在一個水文年變化范圍為6.42~8.06,較適宜水生生物生長;水體在枯水期的電導率最高,豐水期最低。渾河沈撫段水體受沿岸工業廢水和生活污水的點源污染影響,流域枯水期和平水期各監測斷面的綜合水質均未達標;從上游到下游4項水質指標的綜合水質指數增大,污染逐漸加重,主要污染指標為COD和NH3-N。

綜合水質標識指數法;水質評價;污染源;季節性變化;渾河沈撫段

隨著沈撫新城的快速發展,渾河沈撫段已經由自然原生態河流向城市內河轉變,面臨著激增的城市點源與原有農業面源對支流和干流的疊加污染、上游污水廠出水水質較差、河流污徑比較大、生態需水量難以保證等問題,河流自凈能力逐漸降低,生態功能不斷退化,河流兩岸自然生態帶急劇萎縮。

現階段,國內學者們已對渾河水質進行了一定的研究,王輝等[1]對渾河中11項污染物的健康風險的時空變化特征進行分析,得出NH3-N超標嚴重的結論。張鴻齡等[2]利用NH3-N、TP、重金屬3個指標分析渾河水體的污染變化特征,結果表明渾河流域的TP、重金屬污染較輕,渾河下游NH3-N污染較重,水質屬于Ⅴ類甚至劣Ⅴ類。陳艷麗等[3]選用COD、NH3-N、BOD等7項監測項目對渾河干流水質進行評價,得出COD及NH3-N超標為河流部分水質斷面超標的主要原因。

圖1 采樣點位置示意圖

東北地區的氣溫、徑流、降雨量季節性變化較明顯,影響水體污染的因素眾多,監測污染指標的年內平均值無法明確反映水體污染的季節性變化特征[4]。筆者選取2015年3月—2016年2月,對渾河沈撫段枯、豐、平3個水期的污染指標進行分析,利用改進的綜合水質標識指數法識別主要污染因子,以期獲得更具意義的渾河沈撫段水體污染時空變化特征,為鞏固遼河流域重污染“摘帽”成果、進一步提升渾河沈撫段干流水質、實現渾河沈撫段景觀用水水質功能要求提供技術支撐。

1 研究區域概況與方法

1.1 研究區概況

渾河源于遼寧省撫順市清原縣滾馬嶺,流經撫順、沈陽等市縣,本次研究區域為渾河沈撫段區域,指沈陽市東部和撫順市西部的沈撫連接帶,流經該區域的渾河干流水域約11 km2,區域內干流全長32.7 km,干流平均水面寬度約411 m。區域內年徑流量分布不均勻,冬季河面結冰,主要靠地下水補給徑流,4月份之后氣溫逐漸升高,積雪融化,同時上游大伙房水庫放水補給,流量顯著增大。渾河沈撫段接納了沈陽和撫順沿渾河岸邊排放的大量污染物,其環境狀況直接影響著區域居民的健康安全。

1.2 樣品采集和分析

采樣點沿渾河沈撫段上游撫順和平橋開始布控,共設置22個監測斷面,包括10個干流斷面、8個支流斷面、3個排污口斷面,見圖1。從2015年3月到2016年2月對各監測斷面每月進行1次樣品采集,根據《水和廢水監測分析方法》,采用事先凈化好的塑料瓶采集水面下0.5 m水樣,并現場采用YSI6600 V2型多參數水質參數儀測定水溫、濁度、pH值及電導率等指標。每個采集的樣品經過濾后立即加H2SO4(1 mg/L)酸化,在分析樣品前將其保存在4℃的冰箱中,用于測定COD、NH3-N、BOD、CODMn。

1.3 水質評價方法

綜合評價水環境質量的方法有很多,目前較常見的有污染指數法、模糊綜合指數法、人工神經網絡法等,但其共有的不足之處是不能顯示各個水質指標中超標指標的個數,更不能表明河流是否滿足水體功能區水質要求。因此,選用綜合水質標識指數法進行水質評價,其顯著優勢為既可以定性定量評價水體,又不會因個別水質指標較差就否定水體綜合水質,從而對河流綜合水質做出合理評價,計算公式為

IWQ=X1.X2X3X4

(1)

式中:IWQ為綜合水質標識指數;X1為綜合水質級別;X2為綜合水質在該級別水質變化區間中所處的位置;X3為參與綜合水質評價的單項水質指標中,劣于水環境功能目標的指標個數;X4為綜合水質類別與水體功能區類別的比較結果,視綜合水質的污染程度,X4為一位或兩位有效數字[5-6]。

綜合水質標識指數法的重點是計算綜合水質指數X1.X2,為使其更確切反映現狀,初步選定以下3種改進方案[7]。

方案一:

(2)

方案二:

(3)

方案三:

(4)

2 水質指標變化及分析

2.1 pH值在一個水文年內的變化

水體pH值是水體水化學特征的綜合反映。如圖2所示,研究區域在一個水文年內pH值為6.42~8.06,屬中性水平。枯水期,水體pH值變化范圍為7.02~7.56,屬中性水平;豐水期,pH>7,屬中性偏堿水平,最高值出現在和平橋點位;平水期,pH變化范圍為6.42~6.97,屬中性偏酸水平。總體來說,河流水體pH值較適宜水生生物的生長。

圖2 渾河沈撫段枯、豐、平水期pH值變化

從水期變化來看,研究區域水體枯水期和豐水期pH值略高于平水期,且大部分監測斷面豐水期pH值最高。主要由3個原因造成:①對于空氣直接接觸的表層,豐水期溫度較高,CO2在水中的溶解度降低,pH值則上升;②豐水期浮游生物的光合作用強烈消耗水中的CO2導致pH值升高[8];③豐水期上游大伙房水庫放水,大量水流對河底產生沖擊作用,加速水體的解離速度,加速pH值上升。從河道走向來看,豐水期和平水期撫順西段內水體pH值高于沈陽東段,而枯水期規律不明顯,主要因為枯水期水量小、污染重、水質差,水體主要以排污口及支流排放的流量為主,導致各監測區段pH值差異不顯著。

2.2 電導率在一個水文年內的變化

電導率(σ)主要取決于水中的總離子濃度,是水化學分析的重要指標之一。如圖3所示,從水期變化來看,研究區域水體枯水期的電導率最高,其次為平水期,豐水期最低。Anderson等[9]在研究阿拉斯加Kennicott河流流量和電導率時發現,一般水量增大時電導率則降低,水量下降時電導率則升高。本研究中,枯水期的河流電導率明顯高于平水期和豐水期,與前述規律一致[10]。從河道走向來看,枯水期和豐水期撫順西段內電導率略高于沈陽東段,特別是在枯水期高坎大橋監測點達到峰值6.15 mS/m,這主要是因為高坎大橋上游存在下伯官排污口,排放的污水中含有可溶性鹽類,從而使下游電導率升高[11-12]。

圖3 渾河沈撫段枯、豐、平水期電導率變化

3 水質指標評價

渾河水環境功能區目標執行GB3838—2002《地表水環境質量標準》Ⅳ類標準,采用COD、NH3-N、BOD和CODMn共4個項目作為主要評價指標[13],對其水環境狀況進行綜合評價,枯、豐、平3個水期各指標的單因子水質標識指數和水體綜合水質指數X1.X2的3種方案的評價結果見表1。

表1 渾河沈撫段枯、豐、平水期水質綜合評價結果

由表1可知,方案一的各監測斷面的綜合水質指數最小,主要是因其未充分考慮超標因子對水環境質量的影響,將所有因子賦予相同的權重;方案二賦予超標因子過多的權重,導致各監測斷面的綜合水質指數偏大;方案三合理考慮了超標因子對水環境的貢獻力,從而更加貼切地反映河流水體的綜合水質情況,因此本研究采用方案三對渾河沈撫段進行綜合水質標識評價。

3.1 單因子水質標識指數評價

渾河沈撫段枯、豐、平3個水期單因子水質標識指數評價結果見表1。可以看出,在枯水期,10個監測斷面的COD和NH3-N水質污染指標基本上達劣Ⅴ類水標準,絕大多數斷面BOD也處于Ⅴ類水標準,僅CODMn符合水環境功能區目標;豐水期各監測斷面的NH3-N、BOD及CODMn3個水質污染指標均為Ⅳ類水標準,水環境狀況良好,但COD仍均超標于水環境功能區目標,個別監測點劣于Ⅴ類水標準;平水期10個監測斷面的COD幾乎均劣于Ⅴ類水標準,BOD及CODMn水質污染指標均為Ⅳ類水標準,NH3-N變化差異性較大。

對研究區域各個監測斷面4個污染指標的單因子水質標識指數Pi取平均值(圖4)可看出,嚴重超標的污染指標是COD和NH3-N,3個水期中COD均未達到水環境功能區目標,超過水環境功能區2個等級的出現20次;NH3-N指標達到水環境功能區目標的僅有14次,有7次超過水環境功能區2個等級;BOD對河流水體的污染較輕,達到水環境功能區目標的占82.5%。

圖4 各監測斷面4項指標單因子水質標識指數平均值

3.2 綜合水質標識指數評價

采用方案三計算的綜合水質標識指數IWQ結果見表1。從水期變化來看,研究區域枯水期的IWQ最大,水質最差,X1.X2變幅為5.0~5.6,其中最大值出現在高坎大橋斷面,最小值出現在和平橋斷面,10個監測斷面均處于Ⅴ類水標準;平水期次之,X1.X2變幅為5.0~5.6,與枯水期差異較小,表明該區域污染源主要為點源污染;豐水期污染較輕,X1.X2變幅為4.0~5.5,僅高坎大橋和伯官大橋斷面未達到水環境功能區目標。根據表1中水質定性評價結果,枯水期和平水期整個研究流域均處于輕度污染狀態,豐水期情況較好。渾河受上游大伙房水庫影響較大,每年5月份左右水庫持續放水一個月供農業灌溉,之后斷續放水2~3個月,其他時間多為關閘蓄水[14],這是導致豐水期水質良好的主要原因。

研究流域各個監測斷面的IWQ在不同水期之間的變化率見圖5。可看出從枯水期到豐水期的過程中,水質情況改善幅度平均值為15.41%;而當轉為平水期時,水質情況惡化幅度平均值為16.08%,說明影響渾河水質的因素較復雜;從平水期到枯水期,水質情況惡化幅度最大為4.24%,水質基本不變。從圖5中可明顯看出,在高坎大橋斷面附近水質變化相對不明顯,說明此處水質情況主要受污染源影響。

圖5 各監測斷面IWQ隨水期變化率

從河道走向看,研究區域水體污染程度沿程變化較明顯。其中和平橋斷面位于研究區段上游,接納污水量相對較少,3個水期IWQ為4.010~5.021。下游污染情況逐漸加重,流經沈陽后在高坎大橋斷面達到峰值,3個水期IWQ為5.110~5.631,隨后污染情況呈緩慢好轉趨勢。

研究區域各監測斷面之間IWQ的空間變化率見表2。可看出豐水期高陽橡膠壩斷面至高坎大橋斷面IWQ惡化幅度為14.49%;豐水期渾河大橋斷面至高陽橡膠壩斷面IWQ惡化幅度為8.59%;豐水期新立堡立交橋斷面至王家灣攔河壩斷面的IWQ改善幅度為6.51%。

表2 各監測斷面之間IWQ的空間變化率 %

從渾河沈撫段水體IWQ變化規律來看,造成渾河水質變化的污染源較為復雜[15]。研究區域內主要為點源污染,從撫順西到沈陽東沿途有3處排污口、8條支流匯入渾河沈撫段區域[16]。各排污口及支流水質情況見表3。

表3 渾河沈撫段排污口及支流水質調查和評價結果

由表3可知,研究區段3個排污口均劣于Ⅴ類水標準,水質污染嚴重,排放的污水匯入干流,直接導致下游水質變差。8條支流中仁鏡河及舊站河堪比排污口,COD和NH3-N嚴重超標,由于周圍居民基本生活設施匱乏,生活污水和垃圾未經處理直接排放丟棄,河道逐漸成為天然垃圾場,導致水體呈黑色,基本不透明,是使下游高坎大橋斷面水質惡化的直接原因[17-18];蓮島河及滿堂河污染也較重,蓮島河沿岸存在多家化工工廠,其排放的污水導致蓮島河水體渾濁、有泡沫,劣于Ⅴ類水標準;只有東陵大橋斷面附近的白沙河符合水環境功能區目標。此外,渾河沿岸還布有多個生活污水直排口,偷排漏排現象也時有發生,這些都會加劇渾河水質的污染程度[19]。

4 結 論

為解決渾河水生態環境問題,摘掉流域重度污染帽子,本文利用綜合水質標識指數法對渾河沈撫段2015年3月—2016年2月的水質污染指標數據進行分析發現:

a. 研究區域在一個水文年內pH值為6.42~8.06,屬中性水平,較適宜水生生物的生長;水體枯水期的電導率最高,其次為平水期,豐水期最低,電導率沿河道走向變化不顯著。

b. 影響渾河沈撫段的主要污染指標為COD和NH3-N,BOD及CODMn基本上可達水功能區水質要求——Ⅳ類水標準。渾河沈撫段水體僅豐水期的80%監測斷面達到水環境功能區要求,枯水期及平水期各監測斷面的綜合水質均未達標;沿渾河沈撫段河道走向水污染程度表現出一定的規律性,撫順段的綜合水質標識指數總體低于沈陽段。

c. 造成渾河沈撫段水體污染的原因較復雜,從總體來看,渾河污染源主要來自兩岸工業和生活排污口、市政管網污水溢流口和渾河市區段主要支流排入的污水。因此,為防止水質進一步惡化,應加強對渾河的控制管理,以控源截污為重點,加強重點污染企業整治提升,減少污染物排放量;加快污水管網建設,提高污水納管率,從而實現科學治水。

[ 1 ] 王輝,孫麗娜,劉哲,等.渾河水環境健康風險特征研究[J].生態毒理學報,2015,10(2):394-402.(WANG Hui,SUN Lina,LIU Zhe, et al. Study on the character of health risk in the water environment of Hun River[J].Asian Journal of Ecotoxicology,2015, 10(2):394-402.(in Chinese))

[ 2 ] 張鴻齡,孫麗娜,羅慶,等.渾河流域水體污染的季節性變化及來源[J].生態學雜志,2011,30(1):119-125.(ZHANG Hongling, SUN Lina, LUO Qing, et al.Seasonal variations and sources of pollutants in surface water of Hunhe River[J]. Journal of Ecology, 2011,30(1):119-125.(in Chinese))

[ 3 ] 陳艷麗,范遠東,孫鵬.應用模糊綜合評判法評價渾河干流水質[J].水生態學雜志,2012,33(2):39-42.(CHEN Yanli, FAN Yuandong, SUN Peng. Fuzzy comprehensive evaluation on water quality of mainstream of Hun River[J].Journal of Hydroecology, 2012,33(2):39-42.(in Chinese))

[ 4 ] 趙倩,馬建,問青春,等.渾河上游大蘇河鄉農業非點源污染負荷及現狀評價[J].生態與農村環境學報,2010, 26(2):126-131.(ZHAO Qian, MA Jian, WEN Qingchun, et al. Load and status quo of agricultural non-point source pollution in Dasuhe Township on the upper streams of Hun River[J].Journal of Ecology and Rural Environment, 2010, 26(2):126-131.(in Chinese))

[ 5 ] 尹海龍,徐祖信.河流綜合水質評價方法比較研究[J].長江流域資源與環境,2008,17(5):729-733.(YIN Hailong, XU Zuxin. Comparative study on typical river comprehensive water quality assessment methods[J].Resources and Environment in the Yangtze Basin, 2008,17(5):729-733.(in Chinese))

[ 6 ] 徐祖信.我國河流綜合水質標識指數評價方法研究[J].同濟大學學報(自然科學版),2005,33(4):482-488.(XU Zuxin. Comprehensive water quality identification index for environmental quality assessment of surface water[J].Journal of Tongji University(Natural Science), 2005,33(4):482-488.(in Chinese))

[ 7 ] 郭勁松,龍騰銳,霍國友,等.四種水質綜合評價方法的比較[J].重慶建筑大學學報,2000,22(4):6-12.(GUO Jingsong, LONG Tengrui, HUO Guoyou, et al. A comparision of four methods of water quality assessment[J].Journal of Chongqing Jianzhu University, 2000,22(4):6-12.(in Chinese))

[ 8 ] ABSALON D, MATYSIK M. Changes in water quality and runoff in the Upper Oder River Basin. Gecmorphology, 2001,92:106-118.

[ 9 ] ANDERSON S P, LONGACRE S A, KRAAL E R. Patteins of water chemistry and discharge in the glacier fed Kennicott River, Alaska evidence for subglacial water storage cycles[J]. Chemical Geology, 2003, 202:297-312.

[10] SHRESTHA S, KAZAMA F. Assessment of surface water quality using multivariate statistical techniques: a case study of the Fuji river basin, Japan [J]. Environmental Modeling & Software,2007,22:464-476.

[11] VINCENZO N, TIZIANO Z, VINCENZO B. Optimization of sampling frequency for river water quality assessment according to Italian implementation of the EU Water Framework Directive[J].Environmental Science & Policy, 2007,10: 243-249.

[12] 范寧江,劉玉虹,安樹青,等.電導率作為流域水文變化指標初探[J].應用生態學報,2006,17(11):2127-2131.(FAN Ningjing, LIU Yuhong, AN Shuqing, et al. Electrical conductivity as an indicator of hydrological characteristics in catchment scale[J].Journal of Applied Ecology, 2006,17(11):2127-2131.(in Chinese))

[13] 陳利項,傅伯杰,張淑榮,等.于橋水庫流域地表水中水溶性氮季節變化特征[J].中國環境科學,2003,23(2):210-214.(CHEN Lixiang, FU Bojie, ZHANG Shurong, et al. Seasonal change of solvable nitrogen in surface water of Yuqiao Reservoir Basin[J]. China Environment Science, 2003,23(2):210-214.(in Chinese))

[14] 付意成,魏傳江,臧文斌,等.渾太河污染物入河控制量研究[J].水電能源科學,2010, 28(12):21-25.(FU Yicheng, WEI chuangjiang, ZANG Wenbin, et al. Research on contaminant control quantity of entering into Huntai River[J]. Water Resources and Power, 2010, 28(12):21-25.(in Chinese))

[15] 孟偉.中國流域水環境污染綜合防治戰略[J].中國環境科學,2007,27(5):712-716.(MENG Wei. The strategy of comprehensive pollution and cure of water environment in Chinese watershed[J]. China Environment Science, 2007,27(5):712-716.(in Chinese))

[16] 蘇丹,王彤,劉蘭嵐,等.遼河流域工業廢水污染物排放的時空變化規律研究[J]. 生態環境學報,2010,19(12):2953-2959.(SU Dan, WANG tong, LIU Lanlan, et al. Research on the spatio-temporal variation of pollutant discharged from industrial wastewater in the Liaohe River Basin[J].Ecology and Environmental Sciences, 2010,19(12):2953-2959.(in Chinese))

[17] 蘇丹,王志紅,王彤,等.遼河流域工業廢水主要污染物排放強度的單元差異分析[J].生態環境學報,2010, 19(2):275-280.(SU Dan, WANG Zhihong, WANG tong, et al. Analysis on the regional differences of industrial wastewater main pollutants emission intensity in Liao River Basin[J]. Ecology and Environmental Sciences, 2010, 19(2):275-280.(in Chinese))

[18] 陽平堅,吳為中,孟偉,等.基于生態管理的流域水環境功能區劃:以渾河流域為例[J]. 環境科學學報,2007,27(6): 944-952.(YANG Pingjian, WU Weizhong, MENG Wei, et al. Ecosystem management oriented water environmental function zoning at watershed scale: Hun River Basin as a case study[J]. Acta Scientiae Circumstantiae, 2007,27(6): 944-952.(in Chinese))

[19] 武玉峰.遼寧省典型河流健康問題分析及保護對策[J].水利規劃與設計,2015,(9):42-44.(WU Yufeng. Typical river health problem analysis and protection countermeasures in Liaoning Province[J].Water Resources Planning and Design, 2015,(9):42-44.(in Chinese))

Study of seasonal variation of water pollution of Shenfu section of Hunhe River

LIU Qiang, YU Xuedi

(SchoolofMunicipalandEnvironmentalEngineering,ShenyangJianzhuUniversity,Shenyang110168,China)

In order to determine water pollution characteristics of the Shenfu section of the Hunhe River and improve the effectiveness of water environmental treatment, a water environmental investigation was conducted from March 2015 to February 2016. The survey data included pH, conductivity, COD, NH3-N, BOD, and CODMnin the wet, normal, and dry seasons, and the improved comprehensive water quality identification index method was used to identify the main pollution factors and analyze the characteristics of temporal and spatial variations of water pollution in the Shenfu section of the Hunhe River. The results show that the pH value ranged from 6.42 to 7.92 within a hydrological year, which is suitable for the growth of aquatic organisms. The conductivity was highest in the dry season and lowest in the wet season. Affected by point source pollution from industrial wastewater and domestic sewage, the comprehensive water quality of the monitoring sections in the study area in the dry and normal seasons did not reach the standards. The overall water quality index of four water quality indices from upstream to downstream was increased, and pollution was gradually aggravated. The main pollution indices were COD and NH3-N.

comprehensive water quality identification index method; water quality assessment; pollution source; seasonal variation; Shenfu section of Hunhe River

10.3880/j.issn.1004-6933.2017.02.010

國家水體污染控制與治理科技重大專項(2014ZX07202-011)

劉強(1966—),男,教授,主要從事水污染控制理論與技術研究。E-mail:771780475@qq.com

X522

A

1004-6933(2017)02-0052-07

2016-06-07 編輯:王 芳)

猜你喜歡
水質污染
水質抽檢豈容造假
環境(2023年5期)2023-06-30 01:20:01
什么是污染?
什么是污染?
一月冬棚養蝦常見水質渾濁,要如何解決?這9大原因及處理方法你要知曉
當代水產(2019年1期)2019-05-16 02:42:04
這條魚供不應求!蝦蟹養殖戶、垂釣者的最愛,不用投喂,還能凈化水質
當代水產(2019年3期)2019-05-14 05:42:48
堅決打好污染防治攻堅戰
當代陜西(2019年7期)2019-04-25 00:22:18
堅決打好污染防治攻堅戰
圖像識別在水質檢測中的應用
電子制作(2018年14期)2018-08-21 01:38:16
濟下水庫徑流水質和垂向水質分析及評價
對抗塵污染,遠離“霾”伏
都市麗人(2015年5期)2015-03-20 13:33:49
主站蜘蛛池模板: 秘书高跟黑色丝袜国产91在线| 亚洲精品无码不卡在线播放| 色天堂无毒不卡| 成人综合在线观看| 欧美一级黄色影院| 欧美一级爱操视频| 欧美h在线观看| 青青久在线视频免费观看| 日本一区二区三区精品国产| 亚洲伦理一区二区| 在线看片国产| 在线免费亚洲无码视频| 波多野结衣在线se| 亚洲国产精品一区二区高清无码久久| 国产欧美视频在线观看| 亚洲国产91人成在线| 国产乱码精品一区二区三区中文 | 九九热精品在线视频| 国产欧美中文字幕| 九九香蕉视频| 日本免费高清一区| 老司机久久99久久精品播放| 国产精品吹潮在线观看中文| 国产一级一级毛片永久| 久久久波多野结衣av一区二区| 亚洲国产在一区二区三区| 欧美无专区| 国产精品网曝门免费视频| 亚洲美女高潮久久久久久久| 欧美狠狠干| 成人免费午间影院在线观看| 永久免费无码成人网站| 综合亚洲网| 伊在人亚洲香蕉精品播放| 亚洲午夜国产精品无卡| 婷婷99视频精品全部在线观看| 91九色国产porny| 最新国产成人剧情在线播放| 亚洲天堂网在线播放| 美女一级毛片无遮挡内谢| 色精品视频| 国产亚洲欧美另类一区二区| 国产综合精品一区二区| 国产成人精品一区二区三区| 男女男精品视频| 国产你懂得| 久久一本精品久久久ー99| 亚洲人妖在线| 久久国产精品波多野结衣| 亚洲av无码专区久久蜜芽| 无码乱人伦一区二区亚洲一| 亚洲va在线∨a天堂va欧美va| 国产情侣一区二区三区| 丁香婷婷激情网| 日韩在线2020专区| 国产日韩欧美精品区性色| 日韩精品成人在线| 激情综合网址| 一级一毛片a级毛片| 国产欧美日韩va另类在线播放| 国产成人精品一区二区免费看京| 人人爱天天做夜夜爽| 五月天丁香婷婷综合久久| 日本影院一区| 国产精品女熟高潮视频| 亚洲天堂网视频| 国产大全韩国亚洲一区二区三区| 国产色网站| 国产日本一线在线观看免费| 亚洲国产理论片在线播放| 毛片网站观看| 五月天天天色| 五月天在线网站| 欧美成人午夜在线全部免费| 欧美在线一二区| 久久精品丝袜| 在线观看视频99| 黄色网站在线观看无码| 久久久久人妻一区精品色奶水| 视频一区视频二区日韩专区 | 亚洲精品成人片在线观看| 国内丰满少妇猛烈精品播|