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

    汾河河津段主要污染物時間序列分析

    2016-06-09 08:53:55郭新亞張興奇
    中國環(huán)境監(jiān)測 2016年1期

    郭新亞,張興奇

    南京大學(xué)地理與海洋科學(xué)學(xué)院,江蘇 南京 210023

    汾河河津段主要污染物時間序列分析

    郭新亞,張興奇

    南京大學(xué)地理與海洋科學(xué)學(xué)院,江蘇 南京 210023

    基于Mann-Kendall趨勢檢驗法和小波分析法,研究了汾河河津段2004—2013年高錳酸鹽指數(shù)、氨氮、溶解氧含量時間序列的變化趨勢和突變特性。結(jié)果表明:高錳酸鹽指數(shù)和氨氮含量有顯著下降趨勢,溶解氧含量有顯著上升趨勢。高錳酸鹽指數(shù)和氨氮含量在1、2月較大,8、9月較??;溶解氧含量在1、2月較大,6、7月較小。高錳酸鹽指數(shù)含量在2004—2009年差異較大;氨氮含量在2004—2013年差異均較大;溶解氧含量在2010—2013年差異較大。

    高錳酸鹽指數(shù);氨氮;溶解氧;小波分析

    隨著城市化、工業(yè)化的發(fā)展,越來越多的“三廢”排入河流,使得河流污染愈發(fā)嚴重。通過河流監(jiān)測系統(tǒng)可獲取河流水質(zhì)數(shù)據(jù),用于防治和控制河流污染[1]。通過分析河流監(jiān)測長期數(shù)據(jù),研究其時空分布特征,可為流域水環(huán)境管理者提供動態(tài)信息[2]。

    國內(nèi)分析水質(zhì)時間序列的方法主要利用統(tǒng)計學(xué)和數(shù)學(xué)模型,如王菊翠等[3]用統(tǒng)計方法分析了陜西段涇河水質(zhì)時空分布特征;孫國紅等[4]基于Box-Jenkins方法分析了黃河水質(zhì)時間序列并進行了預(yù)測;李冬等[5]基于時間序列分析了黃水庫水質(zhì)漸變性特點。小波分析在水文學(xué)中的應(yīng)用主要集中在降水、徑流等的時間序列變化特性及預(yù)測預(yù)報分析上,用于分析河流污染物含量時間序列變化趨勢及突變特性的較少,張微微等[6]利用小波分析研究了密云水庫上游白河磷濃度的周期變化特征。本文基于Mann-Kendall趨勢檢驗法和小波分析法,研究了汾河河津段高錳酸鹽指數(shù)、氨氮、溶解氧含量時間序列的變化趨勢和突變特性,以期為分析汾河水環(huán)境的變化趨勢提供一定的參考,并對小波分析在河流污染物時間序列分析中的應(yīng)用提供一定的借鑒作用。

    1 實驗部分

    1.1 研究區(qū)概況

    汾河是山西省境內(nèi)最大的河流,黃河的第二大支流,全長713 km,流域面積39 721 km2。汾河發(fā)源于忻州市寧武縣管涔山,由北向南縱貫山西省境內(nèi),在萬榮縣榮河鎮(zhèn)廟前村匯入黃河,如圖1所示。汾河運城段全長145.2 km,是運城市最大河流,河津大橋站是汾河在運城市的一個國控斷面[7],入黃河前的水沙控制站。

    1.2 數(shù)據(jù)資料

    本文數(shù)據(jù)來源于環(huán)境保護部數(shù)據(jù)中心所提供的全國主要流域重點斷面水質(zhì)自動監(jiān)測周報,選取山西運城河津大橋站點2004—2013年每周的高錳酸鹽指數(shù)、氨氮、溶解氧含量數(shù)據(jù),除去斷流及未測到的,高錳酸鹽指數(shù)數(shù)據(jù)累計443周,氨氮數(shù)據(jù)累計446周,溶解氧數(shù)據(jù)累計368周。

    1.3 Mann-Kendall趨勢檢驗法

    Mann-Kendall法是一種非參數(shù)統(tǒng)計檢驗法,因變量可以不具有正態(tài)分布特征,適用于檢驗水文變量的趨勢。假設(shè)有一長度為n的時間序列X1,X2,X3,…,Xn,定義統(tǒng)計量S:

    當Xi-Xj分別大于、等于、小于0時,sign(Xi-Xj)分別為1、0、-1;當n>10時,定義統(tǒng)計量Z:

    Z值的正負性表征變化趨勢的上升與下降,Z的絕對值大小表征趨勢變化的水平,當顯著性水平通過了置信度為90%、95%、99%的檢驗時對應(yīng)的Z值分別為1.28、1.64、2.32[8]。

    1.4 小波分析

    小波分析是在傅里葉變換基礎(chǔ)上引入窗口函數(shù),且時間窗和頻率窗可改變的一種時頻局域分析法,具有時頻多分辨功能。小波分析的基本方法是小波變換,包括連續(xù)小波變換和離散小波變換,連續(xù)小波變換一般需離散化。離散小波變換是逐層對信號進行分解,首先分解長度為M的數(shù)據(jù)序列,得到長度為M/2的低頻部分A1和高頻部分D1;然后對A1進行分解,得到長度為M/4的低頻部分A2和高頻部分D2;一直分解下去,直到一個數(shù)目較小的低頻部分被保留下來[9-10]。本文選用具有良好時頻分析性能的Daubechies小波,Daubechies系中的小波基記為dbN,N為序號,且N=1,2,…,10。

    2 結(jié)果與討論

    2.1 污染物含量時間序列趨勢分析

    近10年汾河河津段高錳酸鹽指數(shù)、氨氮和溶解氧含量變化趨勢如圖2所示。

    從圖2可見,高錳酸鹽指數(shù)最大值出現(xiàn)在2005年第14周(3月28日至4月3日),為320 mg/L,最小值出現(xiàn)在2013年第29周(7月15日至7月21日),為3.7 mg/L,含量有下降趨勢,7 d線性下降幅度為0.301 mg/L;氨氮含量最大值出現(xiàn)在2011年第8周(2月14日至2月20日),為60.7 mg/L,最小值出現(xiàn)在2013年第34周(8月19日至8月25日),為0.2 mg/L,含量有下降趨勢,7 d線性下降幅度為0.018 mg/L;溶解氧含量最大值出現(xiàn)在2013年第52周(2月23日至2月29日),為12.5 mg/L,最小值為0.01 mg/L,集中出現(xiàn)在2004、2008、2009年,含量有上升趨勢,7 d線性上升幅度為0.015 mg/L。

    采用Mann-Kendall趨勢檢驗法得到高錳酸鹽指數(shù)、氨氮、溶解氧的Z值分別為-18.44、-7.04、13.51,三者均達到顯著水平,表明近10 年汾河河津段中高錳酸鹽指數(shù)和氨氮含量呈顯著下降趨勢,溶解氧含量呈顯著上升趨勢。

    雖然高錳酸鹽指數(shù)和氨氮含量有一定幅度的下降,溶解氧含量有一定幅度的上升,水質(zhì)得到了一定的改善,但汾河河津段仍屬于重度污染,其中高錳酸鹽指數(shù)、氨氮、溶解氧含量達到Ⅴ類標準的只占6.5%、3.8%、36.7%,且集中分布在2012、2013年,其余的都屬于劣Ⅴ類。河津位于汾河下游,經(jīng)過中上游的過量開采和超標排放污水,導(dǎo)致汾河進入河津已屬于嚴重污染[11],且河津大力發(fā)展重工業(yè),忽視環(huán)境容量,河流已成為生活污水的排放區(qū),使得汾河河津段的污染進一步加重。2012—2013年高錳酸鹽指數(shù)均值較2004—2011年均值下降62.4 mg/L,達到Ⅳ、Ⅴ類標準的次數(shù)增多,說明對生活廢水和工業(yè)廢水的治理有一定的成效。2012—2013年氨氮含量達到Ⅱ、Ⅴ類的次數(shù)增多,但平均仍為Ⅴ類標準的5.1倍,應(yīng)從點源和面源控制化肥等的使用。2012—2013年溶解氧含量達到Ⅲ類的增加6.3%,溶解氧可反映水體的自凈能力,保持水體自凈不僅需要控制污染物總量,還需要保持一定的流量。面對嚴峻的污染形勢,已取得的治污成果要進一步鞏固,同時應(yīng)采取其他治污措施,加強治污效果。

    小波分析中,對2004—2013年汾河河津段污染物含量的時間序列來說,低頻部分代表著近10年汾河河津段主要污染物的變化規(guī)律,隨著分解尺度的增加,高頻部分信息越來越少,低頻部分代表污染物含量時間序列的變化規(guī)律更加明顯。因此,選擇小波分解和重構(gòu)后低頻系數(shù)最高層可反映污染物含量時間序列的變化規(guī)律。

    采用Daubechies小波,通過實驗確定小波序號N及分解次數(shù),以期能更好地反映污染物含量時間序列的變化規(guī)律。利用db6小波對2004—2013年高錳酸鹽指數(shù)、氨氮、溶解氧含量時間序列進行4次小波分解,再使用第四層的低頻系數(shù)重構(gòu)來確定高錳酸鹽指數(shù)、氨氮、溶解氧含量的變化規(guī)律,如圖3所示。

    圖3 高錳酸鹽指數(shù)、氨氮、溶解氧含量時間序列變化規(guī)律

    高錳酸鹽指數(shù)含量除在2006、2010、2011、2013年變幅較小外,其余年份變幅較大,年內(nèi)最大值與最小值平均相差134.78 mg/L。高錳酸鹽指數(shù)除2006、2012年外,每年均有先減小后增大的趨勢,1、2月較大(均值為72.90 mg/L),8、9月較小(均值為58.10 mg/L)。2006、2012年高錳酸鹽指數(shù)有先增大后減少的趨勢,6、7月含量較大。

    氨氮含量變化趨勢除2006、2012年與高錳酸鹽指數(shù)變化趨勢不同外,其余年份較一致,1、2月含量較大(均值為22.98 mg/L),8、9月含量較小(均值為10.67 mg/L)。2004—2009年氨氮含量最大、最小值變化無一定規(guī)律,2010—2013年氨氮含量最大值逐年減小。

    溶解氧含量在2004—2009年較小,6年均值為0.76 mg/L,在2010—2013年較大,4年均值為3.59 mg/L。溶解氧含量在1、2月較大(均值為2.90 mg/L),6、7月較小(均值為1.08 mg/L)。

    汾河流域?qū)俅箨懶约撅L(fēng)氣候,降雨主要集中在7—9月,枯水期河流水量減小、流速減慢、水環(huán)境容量下降,所以高錳酸鹽指數(shù)和氨氮含量在枯水期1、2月較大,豐水期8、9月較小。2006、2012年高錳酸鹽指數(shù)含量在6、7月較大,說明這段時間內(nèi)人類活動對河流水質(zhì)影響大。溶解氧量與水溫有關(guān),水溫高時溶入的氧量降低,水溫低時溶入的氧量增多,所以溶解氧含量在水溫低的1、2月較大,水溫高的6、7月較小。

    2.2 污染物含量時間序列突變特性分析

    污染物含量時間序列中的突變部分能反映污染比較嚴重的情況,因此分析時間序列的突變部分具有重要的意義。小波分析中,對2004—2013年汾河污染物含量的時間序列來說,高頻部分可以反映污染物含量時間序列的突變部分。本文選用Daubechies小波,其中db1小波具有很好的正則性,即小波很有規(guī)律,可用于檢測污染物含量時間序列的突變特性。小波分解層數(shù)越多,包含的高頻部分越少,一般情況下,第一層和第二層的高頻系數(shù)中包含高頻成分,可清楚地反映出污染物含量時間序列的突變特性。

    利用db1小波對2004—2013年高錳酸鹽指數(shù)、氨氮、溶解氧含量時間序列進行2次小波分解,再使用第一層和第二層的高頻系數(shù)重構(gòu)來確定高錳酸鹽指數(shù)、氨氮、溶解氧含量的突變特性,如圖4~圖6所示。

    圖4 高錳酸鹽指數(shù)含量時間序列突變特性圖

    圖5 氨氮含量時間序列突變特性圖

    圖6 溶解氧含量時間序列突變特性圖

    高錳酸鹽指數(shù)突變位置與數(shù)值較大時的位置相同,說明高錳酸鹽指數(shù)增大一般是突變的。高錳酸鹽指數(shù)在2004—2009年差異較大,2010—2013年差異較小。從第一層和第二層系數(shù)中可看出有2個明顯突變,分別為2005年第14周(3月28日至4月3日)和2009年第14周(3月30日至4月5日),高錳酸鹽指數(shù)分別為320 mg/L和241 mg/L,分別超過Ⅴ類標準的21、16倍,是近10年所有高錳酸鹽指數(shù)中最高的2次,均屬嚴重污染。

    氨氮含量近10年來每年差異均較大,從第二層系數(shù)中可看出有一個明顯突變,位于2011年第8周(2月14日至2月20日),質(zhì)量濃度為60.7 mg/L,超過Ⅴ類標準的30倍。

    溶解氧含量在2004—2009年差異較小,在2010—2013年差異較大。從第一層和第二層系數(shù)中可看出有4個突變點,分別為2009年第24周(6月8—14日)、2011年第7周(2月7—13日)、2013年第7周(2月11—17日)和2013年第32周(8月5—11日),溶解氧分別為4.67、9.01、9.56、9.53 mg/L,分別為V類標準的2.3、4.5、4.78、4.77倍。

    通過觀察小波分解的高頻部分,可了解污染物的突變特性,這些突變部分的污染一般均較嚴重,與突變點左右部分的水質(zhì)相比,由于自然條件一般相差不大,說明這部分水質(zhì)的變化受人類的影響較大。

    3 結(jié)論

    1)汾河河津段中高錳酸鹽指數(shù)和氨氮含量在2004—2013年有顯著下降趨勢,溶解氧含量在2004—2013年有顯著上升趨勢。雖然水質(zhì)得到了一定的改善,但汾河河津段仍屬于重度污染。

    2)高錳酸鹽指數(shù)和氨氮含量在1、2月較大,8、9月較小。溶解氧含量在1、2月較大,6、7月較小。2006、2012年高錳酸鹽指數(shù)在6、7月較大,說明這段時間內(nèi)人類活動對河流水質(zhì)影響較大。

    3)高錳酸鹽指數(shù)在2004—2009年差異較大,在2010—2013年差異較?。话钡吭?004—2013年差異均較大;溶解氧含量在2004—2009年差異較小,在2010—2013年差異較大。

    4)小波分析法可用于研究水質(zhì)時間演化規(guī)律,用于分析水質(zhì)的趨勢變化和突變特性。

    [1] KANNEL P R,LEE S,KANEL S R,et al.Chemometric application in classification and assessment of monitoring locations of an urban river system[J].Analytica Chimica Acta,2007,582:390-399.

    [2] 徐華山,徐宗學(xué),唐芳芳,等.漳衛(wèi)南運河流域水質(zhì) 時空變化特征及其污染源識別[J].環(huán)境科學(xué),2012,33(2):359-369.

    [3] 王菊翠,仵彥卿,黨碧玲,等.基于統(tǒng)計分析的陜西段涇河水質(zhì)時空分布特征[J].自然資源學(xué)報,2012,27(4):674-685.

    [4] 孫國紅,沈躍,徐應(yīng)明,等.基于Box-Jenkins方法的黃河水質(zhì)時間序列分析與預(yù)測[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報,2011,30(9):1 888-1 895.

    [5] 李冬,周川,袁朋飛,等.基于時間序列分析的漸變性水源水質(zhì)預(yù)測研究[J].環(huán)境科學(xué)與技術(shù),2012,35(6):184-188.

    [6] 張微微,李紅,孫丹峰,等.密云水庫上游白河磷濃度時間序列的傅立葉與小波分析[J].中國農(nóng)業(yè)科學(xué),2011,44(10):2 060-2 069.

    [7] 劉宇峰,孫虎,原志華.基于小波分析的汾河河津站徑流與輸沙的多時間尺度特征[J].地理科學(xué),2012,32(6):764-770.

    [8] 琚彤軍,石輝,胡慶.延安市近50年來降水特征及趨勢變化的小波分析研究[J].干旱地區(qū)農(nóng)業(yè)研究,2008,26(4):230-235.

    [9] 陳柳,馬廣大.小波變換在大氣污染物時間序列分析中的應(yīng)用[J].西安科技大學(xué)學(xué)報,2006,26(1):58-61.

    [10] 徐鳴,王斌,呂愛華,等.大氣污染物多時間分辨率的小波分析[J].環(huán)境科學(xué)學(xué)報,2008,28(4):786-790.

    [11] 張引栓.汾河運城段水質(zhì)評價及排污口排放標準分析[J].人民黃河,2010,32(7):54-55.

    Time Series Analysis of Main Pollutants in Hejin Section of Fenhe River Based On Mann-Kendall Method and Wavelet Analysis

    GUO Xinya, ZHANG Xingqi

    School of Geographic and Oceanographic Sciences at Nanjing University, Nanjing 210023, China

    The time series trends and mutation features of CODMn, NH3-N and DO in Hejin section of Fenhe River were analyzed via Mann-Kendall method and Wavelet Analysis. The results indicated that:The concentrations of CODMnand NH3-N had a significant downward trend and the concentration of DO had a significant upward trend. The concentrations of CODMnand NH3-N were large in January and February and were low in August and September; the concentration of DO was large in January and February and was low in June and July. The concentration of CODMnhad significant differences between 2004 and 2009; there were differences on the concentration of NH3-N between 2004 and 2013; the concentration of DO had significant differences between 2010 and 2013.

    CODMn; NH3-N; DO; Wavelet Analysis

    2015-01-05;

    2015-05-24

    郭新亞(1990-),女,山西臨汾人,碩士。

    張興奇

    X824

    A

    1002-6002(2016)01- 0070- 05

    精品久久久精品久久久| av国产精品久久久久影院| 卡戴珊不雅视频在线播放| √禁漫天堂资源中文www| 精品久久久久久久久av| 亚洲国产av新网站| 狂野欧美激情性xxxx在线观看| 日本爱情动作片www.在线观看| 亚洲欧美日韩卡通动漫| 热re99久久精品国产66热6| 久久久国产欧美日韩av| 波野结衣二区三区在线| 亚洲av中文av极速乱| 日本wwww免费看| 免费大片18禁| 少妇熟女欧美另类| 精品人妻在线不人妻| 成人亚洲精品一区在线观看| 少妇丰满av| 18禁在线播放成人免费| 国产老妇伦熟女老妇高清| 亚洲精品456在线播放app| 91成人精品电影| 国产淫语在线视频| 亚洲av欧美aⅴ国产| 桃花免费在线播放| 久久 成人 亚洲| 一级黄片播放器| 大又大粗又爽又黄少妇毛片口| 免费大片黄手机在线观看| 超色免费av| 多毛熟女@视频| 99久久精品国产国产毛片| 美女视频免费永久观看网站| 成人午夜精彩视频在线观看| 婷婷色麻豆天堂久久| 亚洲综合色网址| av电影中文网址| 久久久久久人妻| 高清视频免费观看一区二区| 国产一区二区三区av在线| 国产精品欧美亚洲77777| 免费观看的影片在线观看| 自线自在国产av| 久久精品国产亚洲av涩爱| 亚洲av中文av极速乱| 国产男人的电影天堂91| 亚洲国产色片| 成人午夜精彩视频在线观看| 久久99精品国语久久久| 青青草视频在线视频观看| 国产精品99久久99久久久不卡 | 国产精品国产三级专区第一集| 啦啦啦在线观看免费高清www| 亚洲国产最新在线播放| 老女人水多毛片| 国产精品一区二区在线观看99| 一区在线观看完整版| 亚洲av电影在线观看一区二区三区| 国产成人aa在线观看| 男女无遮挡免费网站观看| 一级毛片黄色毛片免费观看视频| 亚洲综合色网址| 免费日韩欧美在线观看| 青春草视频在线免费观看| 欧美xxⅹ黑人| 一级片'在线观看视频| 中文字幕精品免费在线观看视频 | 国产白丝娇喘喷水9色精品| 日韩av在线免费看完整版不卡| 国产国拍精品亚洲av在线观看| 欧美 亚洲 国产 日韩一| 亚洲精品乱码久久久v下载方式| 国产精品一区二区三区四区免费观看| 草草在线视频免费看| 91精品国产国语对白视频| 18禁动态无遮挡网站| 国产毛片在线视频| 国产极品粉嫩免费观看在线 | 国产一区二区三区av在线| 免费久久久久久久精品成人欧美视频 | av福利片在线| 成人毛片60女人毛片免费| 在线精品无人区一区二区三| 国产成人精品一,二区| 亚洲精品中文字幕在线视频| av女优亚洲男人天堂| 午夜免费鲁丝| 亚洲婷婷狠狠爱综合网| 91成人精品电影| 国产av精品麻豆| 久久久久久久久大av| 亚洲国产最新在线播放| 午夜免费观看性视频| 高清黄色对白视频在线免费看| 久久久国产一区二区| 女性被躁到高潮视频| 成年女人在线观看亚洲视频| 最新的欧美精品一区二区| 亚洲av电影在线观看一区二区三区| 三级国产精品片| 久久av网站| 如日韩欧美国产精品一区二区三区 | 午夜激情av网站| 人妻 亚洲 视频| 我的老师免费观看完整版| 中文精品一卡2卡3卡4更新| 中文乱码字字幕精品一区二区三区| 国产免费又黄又爽又色| 最新的欧美精品一区二区| 亚洲国产精品一区三区| 久久久久久久久久久丰满| 精品国产露脸久久av麻豆| 桃花免费在线播放| 最近2019中文字幕mv第一页| 国产成人freesex在线| 婷婷色综合大香蕉| 最近手机中文字幕大全| 亚洲人成网站在线观看播放| 日本色播在线视频| 一本大道久久a久久精品| 午夜影院在线不卡| 国产精品偷伦视频观看了| 久久99热这里只频精品6学生| 在线亚洲精品国产二区图片欧美 | 国产 精品1| 久久久欧美国产精品| 久久99蜜桃精品久久| 亚洲精品久久午夜乱码| 国产亚洲av片在线观看秒播厂| 制服丝袜香蕉在线| 国产精品免费大片| 如何舔出高潮| 免费观看在线日韩| 成人18禁高潮啪啪吃奶动态图 | 亚洲av国产av综合av卡| 黄色毛片三级朝国网站| 国产色爽女视频免费观看| 亚洲欧美精品自产自拍| 大码成人一级视频| 亚洲精品色激情综合| 亚洲精华国产精华液的使用体验| 成人二区视频| 插阴视频在线观看视频| 免费观看的影片在线观看| 免费大片18禁| 欧美另类一区| 国产亚洲欧美精品永久| 最近最新中文字幕免费大全7| 亚洲美女视频黄频| 最近最新中文字幕免费大全7| 日韩大片免费观看网站| 大码成人一级视频| 九色成人免费人妻av| 国产精品三级大全| 青春草亚洲视频在线观看| 一级a做视频免费观看| 国产av国产精品国产| 亚洲精品乱码久久久v下载方式| 欧美三级亚洲精品| 91在线精品国自产拍蜜月| 国产深夜福利视频在线观看| 亚洲精品亚洲一区二区| 卡戴珊不雅视频在线播放| 国产精品麻豆人妻色哟哟久久| 美女大奶头黄色视频| 久久久久精品久久久久真实原创| 亚洲av中文av极速乱| videos熟女内射| 亚洲久久久国产精品| 99视频精品全部免费 在线| 男女边摸边吃奶| 日韩中文字幕视频在线看片| 欧美日韩国产mv在线观看视频| 亚洲不卡免费看| 秋霞伦理黄片| 国产成人精品一,二区| 亚洲精品久久成人aⅴ小说 | 人人妻人人添人人爽欧美一区卜| av黄色大香蕉| 国产精品无大码| 丝袜美足系列| 91精品一卡2卡3卡4卡| 国产 一区精品| 丰满少妇做爰视频| 国产精品99久久久久久久久| 午夜激情久久久久久久| 欧美 日韩 精品 国产| 高清午夜精品一区二区三区| 黄色怎么调成土黄色| 99久久人妻综合| 国产成人精品无人区| 午夜免费鲁丝| 亚洲精品久久成人aⅴ小说 | 美女福利国产在线| 色视频在线一区二区三区| 久久午夜综合久久蜜桃| 亚洲国产精品专区欧美| 极品人妻少妇av视频| 女人精品久久久久毛片| 久久久久久久精品精品| 99热国产这里只有精品6| 成人午夜精彩视频在线观看| 免费黄网站久久成人精品| 日韩免费高清中文字幕av| 国产成人免费无遮挡视频| 一本久久精品| 久热久热在线精品观看| 老女人水多毛片| 97超碰精品成人国产| 性高湖久久久久久久久免费观看| 永久网站在线| 国产精品国产三级专区第一集| 亚洲怡红院男人天堂| 免费观看无遮挡的男女| 久久人人爽人人片av| 赤兔流量卡办理| 天天操日日干夜夜撸| av国产精品久久久久影院| 美女中出高潮动态图| 哪个播放器可以免费观看大片| 婷婷色综合www| 熟女av电影| 美女国产视频在线观看| 亚洲成人一二三区av| 亚洲欧美精品自产自拍| 亚洲,一卡二卡三卡| 美女中出高潮动态图| 免费播放大片免费观看视频在线观看| 国产亚洲一区二区精品| 成人毛片a级毛片在线播放| 亚洲精品av麻豆狂野| 男女无遮挡免费网站观看| 国产国语露脸激情在线看| 欧美日韩av久久| 国产精品.久久久| 国产亚洲午夜精品一区二区久久| 免费黄网站久久成人精品| 视频在线观看一区二区三区| 一本大道久久a久久精品| 免费看光身美女| 一级片'在线观看视频| 久久久久久久国产电影| 这个男人来自地球电影免费观看 | 亚洲中文av在线| 成人国语在线视频| 少妇熟女欧美另类| 91久久精品电影网| 丰满少妇做爰视频| 女性生殖器流出的白浆| 国产黄片视频在线免费观看| 久久这里有精品视频免费| 一区二区av电影网| 国产日韩欧美视频二区| 亚洲精品aⅴ在线观看| 精品卡一卡二卡四卡免费| 色5月婷婷丁香| 亚洲国产欧美在线一区| 久久久国产欧美日韩av| 色婷婷av一区二区三区视频| 亚洲av综合色区一区| 两个人的视频大全免费| 夜夜骑夜夜射夜夜干| 九草在线视频观看| 亚洲精品久久午夜乱码| 18+在线观看网站| 亚洲精品色激情综合| 三上悠亚av全集在线观看| 2022亚洲国产成人精品| 亚洲国产精品成人久久小说| 夜夜看夜夜爽夜夜摸| 久久99热这里只频精品6学生| 成人午夜精彩视频在线观看| 国产亚洲欧美精品永久| 久久国产精品男人的天堂亚洲 | 能在线免费看毛片的网站| 国产成人一区二区在线| 高清av免费在线| 高清在线视频一区二区三区| 99热国产这里只有精品6| 99热这里只有是精品在线观看| 免费观看av网站的网址| 卡戴珊不雅视频在线播放| 亚洲性久久影院| 国产精品欧美亚洲77777| 丝袜喷水一区| 日本vs欧美在线观看视频| 国产精品久久久久久精品古装| 久久久亚洲精品成人影院| 中文字幕精品免费在线观看视频 | 久久久久久伊人网av| 亚洲精品456在线播放app| 亚洲av国产av综合av卡| 狂野欧美白嫩少妇大欣赏| av.在线天堂| 丝袜美足系列| 波野结衣二区三区在线| 一区二区三区四区激情视频| 人妻制服诱惑在线中文字幕| av国产精品久久久久影院| 国产伦理片在线播放av一区| 久久久精品94久久精品| 成人午夜精彩视频在线观看| 欧美国产精品一级二级三级| 亚洲精品久久成人aⅴ小说 | 欧美精品高潮呻吟av久久| 搡女人真爽免费视频火全软件| 建设人人有责人人尽责人人享有的| 成人毛片60女人毛片免费| 亚洲精品国产av成人精品| 国产成人精品在线电影| 精品少妇内射三级| 亚洲经典国产精华液单| 毛片一级片免费看久久久久| 久久久久久久久久久免费av| 国产精品久久久久久av不卡| 成人亚洲精品一区在线观看| 精品卡一卡二卡四卡免费| 看免费成人av毛片| 亚洲,一卡二卡三卡| 伦精品一区二区三区| 狂野欧美激情性xxxx在线观看| av专区在线播放| 日韩中文字幕视频在线看片| 天天影视国产精品| 国产一区亚洲一区在线观看| tube8黄色片| 女人精品久久久久毛片| av在线老鸭窝| 亚洲精品乱久久久久久| 精品亚洲成a人片在线观看| 亚洲美女视频黄频| 国产乱人偷精品视频| 看非洲黑人一级黄片| 欧美3d第一页| 观看av在线不卡| 国产无遮挡羞羞视频在线观看| 最新的欧美精品一区二区| 亚洲人与动物交配视频| 最新的欧美精品一区二区| 欧美日韩av久久| .国产精品久久| 久久 成人 亚洲| 色5月婷婷丁香| 久久ye,这里只有精品| 在线观看免费视频网站a站| 国产亚洲欧美精品永久| 成人午夜精彩视频在线观看| 国产高清国产精品国产三级| 国产又色又爽无遮挡免| av一本久久久久| 欧美变态另类bdsm刘玥| 最近最新中文字幕免费大全7| 日日啪夜夜爽| 日韩制服骚丝袜av| 亚洲精品乱码久久久v下载方式| 3wmmmm亚洲av在线观看| 中文字幕免费在线视频6| 3wmmmm亚洲av在线观看| 狠狠精品人妻久久久久久综合| 亚洲人成77777在线视频| 一区二区三区四区激情视频| 国产在线免费精品| 国产成人a∨麻豆精品| 女性生殖器流出的白浆| 丰满少妇做爰视频| 亚洲欧美精品自产自拍| 水蜜桃什么品种好| 欧美亚洲 丝袜 人妻 在线| 秋霞在线观看毛片| 日本与韩国留学比较| 麻豆精品久久久久久蜜桃| 免费人妻精品一区二区三区视频| 妹子高潮喷水视频| 69精品国产乱码久久久| 18禁观看日本| www.色视频.com| 涩涩av久久男人的天堂| 国产男人的电影天堂91| av视频免费观看在线观看| 丰满乱子伦码专区| 欧美+日韩+精品| 青青草视频在线视频观看| 一个人看视频在线观看www免费| 欧美日韩视频高清一区二区三区二| 99re6热这里在线精品视频| 在线观看美女被高潮喷水网站| 精品少妇黑人巨大在线播放| 黄色毛片三级朝国网站| 插逼视频在线观看| 色婷婷久久久亚洲欧美| 日韩av不卡免费在线播放| 亚洲精品第二区| 汤姆久久久久久久影院中文字幕| 蜜桃国产av成人99| 成人国产麻豆网| 伊人久久国产一区二区| 免费不卡的大黄色大毛片视频在线观看| 久久久久久久亚洲中文字幕| 精品久久国产蜜桃| 熟女av电影| 国产成人av激情在线播放 | 久久ye,这里只有精品| 在线播放无遮挡| 亚洲av电影在线观看一区二区三区| 麻豆成人av视频| 日韩成人av中文字幕在线观看| 九草在线视频观看| 一区二区三区乱码不卡18| 免费观看在线日韩| 欧美xxxx性猛交bbbb| 只有这里有精品99| 国产精品蜜桃在线观看| 久久精品久久久久久噜噜老黄| 另类精品久久| 日韩制服骚丝袜av| 精品一区二区三卡| 国产亚洲最大av| 亚洲无线观看免费| 91精品国产国语对白视频| 久久久久久久久久人人人人人人| 国内精品宾馆在线| 交换朋友夫妻互换小说| 青春草国产在线视频| 亚洲精品日韩av片在线观看| 亚洲精品一区蜜桃| 中文精品一卡2卡3卡4更新| 考比视频在线观看| 久久青草综合色| 成人亚洲精品一区在线观看| 在现免费观看毛片| freevideosex欧美| 欧美日韩亚洲高清精品| 国产欧美日韩综合在线一区二区| 99久久人妻综合| 国产精品人妻久久久影院| 成人午夜精彩视频在线观看| 最近中文字幕高清免费大全6| 成人亚洲精品一区在线观看| 国产精品嫩草影院av在线观看| 日韩人妻高清精品专区| 国产精品久久久久久久久免| 亚洲伊人久久精品综合| 免费人成在线观看视频色| 中国国产av一级| 岛国毛片在线播放| av不卡在线播放| 亚洲精品自拍成人| 国产亚洲一区二区精品| 亚洲怡红院男人天堂| 欧美一级a爱片免费观看看| 国产亚洲最大av| 亚洲国产日韩一区二区| 精品久久久精品久久久| 欧美精品一区二区免费开放| 最近手机中文字幕大全| 亚洲av中文av极速乱| 啦啦啦啦在线视频资源| 亚洲国产精品专区欧美| 在线观看免费视频网站a站| 18禁在线播放成人免费| 老司机影院毛片| 婷婷色综合www| 日韩一本色道免费dvd| 欧美亚洲日本最大视频资源| 色吧在线观看| 人妻 亚洲 视频| xxxhd国产人妻xxx| 黑丝袜美女国产一区| 老司机影院成人| 只有这里有精品99| 亚洲精品一区蜜桃| 97在线视频观看| 亚洲国产日韩一区二区| 大香蕉97超碰在线| 国产精品成人在线| 十八禁高潮呻吟视频| 国产男女内射视频| 成年人免费黄色播放视频| 国语对白做爰xxxⅹ性视频网站| 国产精品 国内视频| 亚洲国产欧美日韩在线播放| 亚洲精品乱久久久久久| 熟女av电影| 丝袜美足系列| 久久免费观看电影| 国产爽快片一区二区三区| 在线观看三级黄色| 亚洲久久久国产精品| 一级毛片我不卡| 久久免费观看电影| 熟妇人妻不卡中文字幕| 久久久国产一区二区| 欧美激情国产日韩精品一区| 最后的刺客免费高清国语| 久久热精品热| 狠狠精品人妻久久久久久综合| xxx大片免费视频| 国产不卡av网站在线观看| 久久 成人 亚洲| 一级毛片黄色毛片免费观看视频| 韩国av在线不卡| 国产极品粉嫩免费观看在线 | 草草在线视频免费看| 免费人成在线观看视频色| 久久99热这里只频精品6学生| 中文字幕人妻丝袜制服| 在线免费观看不下载黄p国产| 人人妻人人澡人人看| 精品一品国产午夜福利视频| 免费大片18禁| 欧美三级亚洲精品| 欧美丝袜亚洲另类| 一级毛片 在线播放| 亚洲精品一二三| 国模一区二区三区四区视频| 久久久午夜欧美精品| 亚洲美女搞黄在线观看| 97在线视频观看| 久久久久久久久久久免费av| 国产黄色视频一区二区在线观看| 少妇人妻 视频| 男女无遮挡免费网站观看| 91在线精品国自产拍蜜月| 亚洲国产精品一区二区三区在线| 亚洲精品456在线播放app| 国产欧美另类精品又又久久亚洲欧美| 熟女av电影| 插阴视频在线观看视频| 国产黄频视频在线观看| 国产精品麻豆人妻色哟哟久久| 亚洲情色 制服丝袜| 在线天堂最新版资源| 天堂俺去俺来也www色官网| 99热网站在线观看| 国产伦理片在线播放av一区| 高清欧美精品videossex| 国产精品成人在线| 国产成人精品一,二区| 黄色怎么调成土黄色| 在线看a的网站| 久热久热在线精品观看| 高清av免费在线| 国产精品一区二区三区四区免费观看| 国产精品一区www在线观看| 最新中文字幕久久久久| 人妻 亚洲 视频| 精品久久国产蜜桃| 成人亚洲精品一区在线观看| 一级黄片播放器| 亚洲欧美成人综合另类久久久| 青春草视频在线免费观看| 日日爽夜夜爽网站| 久久久亚洲精品成人影院| 亚洲精品一区蜜桃| 亚洲国产精品一区三区| 免费观看在线日韩| 欧美 日韩 精品 国产| 亚洲内射少妇av| 色94色欧美一区二区| 亚洲欧美日韩另类电影网站| 99久久中文字幕三级久久日本| 97超碰精品成人国产| 最近手机中文字幕大全| 一区在线观看完整版| 免费人成在线观看视频色| 久久韩国三级中文字幕| a级毛片在线看网站| 国产欧美日韩一区二区三区在线 | 国产在线免费精品| 国产黄色视频一区二区在线观看| 日本vs欧美在线观看视频| 亚洲性久久影院| 国产老妇伦熟女老妇高清| 免费播放大片免费观看视频在线观看| 国产精品国产av在线观看| 性色avwww在线观看| 制服诱惑二区| 亚洲精品久久成人aⅴ小说 | 成年人免费黄色播放视频| 欧美精品亚洲一区二区| 久久久久久久久大av| 亚洲久久久国产精品| 亚洲国产色片| 亚洲伊人久久精品综合| 在线精品无人区一区二区三| 黄色视频在线播放观看不卡| 亚洲综合色惰| 国产探花极品一区二区| 日韩 亚洲 欧美在线| 一级毛片 在线播放| 国产一区二区三区综合在线观看 | 亚洲欧美成人精品一区二区| 欧美 亚洲 国产 日韩一| 看非洲黑人一级黄片| 国产成人午夜福利电影在线观看| 美女cb高潮喷水在线观看| av免费观看日本| 赤兔流量卡办理| 我要看黄色一级片免费的| 日韩伦理黄色片| 青春草视频在线免费观看| 午夜免费鲁丝| 国产有黄有色有爽视频| 久久精品夜色国产| 麻豆成人av视频| 国产成人精品婷婷| 国产精品麻豆人妻色哟哟久久| 中文精品一卡2卡3卡4更新| 尾随美女入室| 国产精品久久久久久av不卡| 亚洲国产精品一区二区三区在线| 九色亚洲精品在线播放| 久久精品国产自在天天线| 亚洲欧美色中文字幕在线|