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

    基于K-means聚類和熵值法的二次供水水質(zhì)評估

    2023-11-17 07:59:36朱延平吳霖璟于佳麗顧趙福舒詩湖
    關(guān)鍵詞:水質(zhì)設(shè)備分析

    朱延平,吳霖璟,于佳麗,顧趙福,舒詩湖

    (1.東華大學(xué)環(huán)境科學(xué)與工程學(xué)院, 上海 201620;2.上海市供水管理事務(wù)中心, 上海 200092;3.上海污染控制與生態(tài)安全研究院, 上海 200092)

    在城市化進(jìn)程中,我國部分城市通過構(gòu)建水質(zhì)在線監(jiān)測系統(tǒng)以應(yīng)對日益凸顯的二次供水環(huán)節(jié)的問題。目前,上海市已在部分住宅小區(qū)安裝超過200個(gè)二次供水水質(zhì)在線監(jiān)測儀表,遠(yuǎn)期將實(shí)現(xiàn)每小區(qū)至少安裝1套水質(zhì)在線儀表的目標(biāo)。在線水質(zhì)監(jiān)測儀可實(shí)現(xiàn)對二次供水水質(zhì)連續(xù)采樣、實(shí)時(shí)檢測、數(shù)據(jù)反饋和信息處理等任務(wù)[1]。然而由于設(shè)備異常、外界干擾等原因,采集的數(shù)據(jù)有時(shí)不能滿足分析的需求,因此,需要梳理二次供水在線監(jiān)測數(shù)據(jù)中常見的錯(cuò)誤,并進(jìn)行數(shù)據(jù)質(zhì)量評估,從而支撐二次供水的智慧化管理。

    經(jīng)過多年水質(zhì)提升工程的實(shí)施,二次供水的水質(zhì)可滿足微生物指標(biāo)、毒理指標(biāo)、放射性指標(biāo)等國家標(biāo)準(zhǔn)[2]。二次供水水質(zhì)在線監(jiān)測系統(tǒng)需要滿足監(jiān)測數(shù)據(jù)容量大、實(shí)時(shí)性強(qiáng)等特點(diǎn),同時(shí)具有監(jiān)測對象有限、信息密度低、穩(wěn)定性不足等問題。常見的飲用水水質(zhì)分析方法有神經(jīng)網(wǎng)絡(luò)法[3]、投影尋蹤方法[4]、灰色指數(shù)法[5]、物元分析法[6-7]等,不適用于分析二次供水在線監(jiān)測水質(zhì)。因此,二次供水水質(zhì)在線監(jiān)測大數(shù)據(jù)應(yīng)選擇快速、有效的處理方法,以便后續(xù)分析和水質(zhì)預(yù)警工作的展開。K-means聚類分析被廣泛應(yīng)用于水質(zhì)分析[8],其顯著優(yōu)點(diǎn)是可以定量和客觀地應(yīng)用聚類有效性指標(biāo)評估聚類結(jié)果[9]。熵值法在確定權(quán)重上可以從某種程度上避免傳統(tǒng)方法在權(quán)重確定中主觀因素的影響[10-11]。本文研究基于K-means聚類法和熵值法的二次供水水質(zhì)在線監(jiān)測評價(jià)方法及其應(yīng)用。

    1 二次供水在線監(jiān)測數(shù)據(jù)質(zhì)量分析

    二次供水監(jiān)測點(diǎn)通常設(shè)置在泵房泵后、屋頂水箱、一體化箱式泵站泵后的出水口,方便在配電箱或電梯機(jī)房就近取電。因此,致錯(cuò)因素大概可歸納為3類:1)環(huán)境因素,主要影響因子為通訊信號(hào)和工作環(huán)境溫度;2)設(shè)備因素,主要影響因子為設(shè)備時(shí)鐘故障、遠(yuǎn)程傳輸單元、電池狀態(tài)和設(shè)備穩(wěn)定性能;3)其他因素,主要?dú)w納為自互聯(lián)網(wǎng)攻擊風(fēng)險(xiǎn)和設(shè)備運(yùn)行管理兩方面的影響[12]。

    常見的二次供水水質(zhì)在線監(jiān)測數(shù)據(jù)根據(jù)以上3類致錯(cuò)因素又可分為4種數(shù)據(jù)錯(cuò)誤類型,分別為時(shí)間間斷型、時(shí)間錯(cuò)位型、記錄缺失型、短暫數(shù)值尖峰型(又稱“毛刺”)。4類數(shù)據(jù)錯(cuò)誤主要可歸納為2大數(shù)據(jù)特征:數(shù)據(jù)特征對數(shù)據(jù)質(zhì)量的影響較大,但易于識(shí)別和發(fā)現(xiàn);數(shù)據(jù)特征對數(shù)據(jù)質(zhì)量的影響較小,但排查困難、致錯(cuò)環(huán)節(jié)復(fù)雜。具體的數(shù)據(jù)錯(cuò)誤類型、致錯(cuò)因素和數(shù)據(jù)錯(cuò)誤特征三者的關(guān)系如圖1所示。

    圖1 致錯(cuò)因素、數(shù)據(jù)錯(cuò)誤類型、數(shù)據(jù)錯(cuò)誤特征關(guān)系圖Fig.1 Factors of errors, types of wrong data and data error characteristics diagrams

    1.1 時(shí)間間斷型數(shù)據(jù)錯(cuò)誤

    運(yùn)行正常的智能設(shè)備根據(jù)設(shè)定時(shí)間記錄數(shù)據(jù),當(dāng)上傳數(shù)據(jù)的時(shí)間間隔超過用戶設(shè)定時(shí),可視為時(shí)間間斷型數(shù)據(jù)錯(cuò)誤,根據(jù)間斷時(shí)間的長短可以分為數(shù)分鐘、數(shù)小時(shí)和數(shù)天。

    時(shí)間間斷型數(shù)據(jù)錯(cuò)誤主要由2種因素引起數(shù)據(jù)故障:1)較弱的通訊信號(hào)可能造成監(jiān)測設(shè)備未能在傳輸時(shí)間內(nèi)上傳全部數(shù)據(jù)[13];2)盡管在線監(jiān)測設(shè)備的就近取電方式可以保障穩(wěn)定的運(yùn)行,但電池狀態(tài)可能會(huì)導(dǎo)致數(shù)據(jù)中斷,一般水質(zhì)在線監(jiān)測設(shè)備的電池壽命為5年左右。

    1.2 時(shí)間錯(cuò)位型數(shù)據(jù)錯(cuò)誤

    相鄰數(shù)據(jù)時(shí)間間隔正確,但數(shù)據(jù)記錄的時(shí)刻發(fā)生了改變,這種數(shù)據(jù)錯(cuò)誤被認(rèn)為是時(shí)間錯(cuò)位型錯(cuò)誤,主要?dú)w因于設(shè)備的遠(yuǎn)程傳輸單元的時(shí)鐘故障。

    1.3 記錄缺失型數(shù)據(jù)錯(cuò)誤

    當(dāng)對應(yīng)記錄時(shí)間的記錄值為空白或“0”認(rèn)定為缺失型數(shù)據(jù)錯(cuò)誤,此類錯(cuò)誤通常是由于設(shè)備故障或環(huán)境惡劣引起的數(shù)據(jù)“丟包”導(dǎo)致。該類數(shù)據(jù)錯(cuò)誤一般包含2種情況:1)水質(zhì)在線監(jiān)測儀器對溫度敏感,通常工作環(huán)境溫度為5~40 ℃,當(dāng)工作環(huán)境溫度超出適用范圍則可能導(dǎo)致設(shè)備靈敏度下降,甚至失靈(這類問題可以通過設(shè)備的維護(hù)管理工作得到改善);2)由于設(shè)備自身的遠(yuǎn)程傳輸單元故障,導(dǎo)致設(shè)備上傳錯(cuò)誤的檢測值至工作站。

    1.4 短暫數(shù)值尖峰型數(shù)據(jù)錯(cuò)誤

    短暫數(shù)值尖峰類型的錯(cuò)誤表現(xiàn)為部分記錄數(shù)值在時(shí)間序列上呈“毛刺”狀,這通常與設(shè)備穩(wěn)定性和適應(yīng)監(jiān)測環(huán)境能力有關(guān):1)剛運(yùn)行的智能監(jiān)測設(shè)備,監(jiān)測探頭或分析單元的穩(wěn)定性不足導(dǎo)致出現(xiàn)許多極端數(shù)據(jù);2)剛運(yùn)行遠(yuǎn)傳水表的監(jiān)測靈敏性較強(qiáng),對監(jiān)測環(huán)境敏感。

    1.5 數(shù)據(jù)長期無波動(dòng)型數(shù)據(jù)錯(cuò)誤

    長期無波動(dòng)型數(shù)據(jù)錯(cuò)誤表現(xiàn)為同一監(jiān)測點(diǎn)的時(shí)間相鄰檢測差值超過1 h的情況下始終為0,這與工作環(huán)境(環(huán)境溫度)和運(yùn)維管理有關(guān):1)監(jiān)測設(shè)備未及時(shí)進(jìn)行設(shè)備管理導(dǎo)致靈敏度下降;2)工作環(huán)境的惡劣造成智能設(shè)備故障(如檢測探頭失靈)導(dǎo)致分析結(jié)果始終不變。

    1.6 數(shù)據(jù)質(zhì)量結(jié)果與分析

    本研究對L市的36個(gè)監(jiān)測點(diǎn)(檢測時(shí)間為2020年1—10月)的數(shù)據(jù)監(jiān)測結(jié)果進(jìn)行質(zhì)量評估分析,具體數(shù)據(jù)質(zhì)量說明和數(shù)據(jù)錯(cuò)誤的判定標(biāo)準(zhǔn)如表1所示。

    表1 數(shù)據(jù)質(zhì)量分析總體結(jié)果匯總Table 1 Summary of overall results of data quality analysis

    由表1的數(shù)據(jù)質(zhì)量分析結(jié)果可知:1)間斷數(shù)分鐘的時(shí)間間斷型錯(cuò)誤是普遍存在的數(shù)據(jù)錯(cuò)誤類型;2)記錄缺失型數(shù)據(jù)是數(shù)據(jù)錯(cuò)誤的首要類型,使用GPRS會(huì)發(fā)生包丟失現(xiàn)象;3)不同在線監(jiān)測設(shè)備所導(dǎo)致的數(shù)據(jù)錯(cuò)誤程度也不同,但此類錯(cuò)誤一般可以通過加強(qiáng)設(shè)備運(yùn)行維護(hù)管理和網(wǎng)絡(luò)安全管理的方式優(yōu)化。

    2 二次供水水質(zhì)在線監(jiān)測數(shù)據(jù)評價(jià)

    L市由4個(gè)不同的水源地供水,按二次供水方式分為水箱水池聯(lián)合供水、水池加變頻泵供水、無負(fù)壓供水等方式。供水企業(yè)全面接管2000年以前建設(shè)并在近年完成改造的二次供水設(shè)施,同步建設(shè)了共36個(gè)二次供水水質(zhì)在線監(jiān)測系統(tǒng)。政府出于對二次供水監(jiān)管的目的,自主安裝了若干水質(zhì)在線監(jiān)管點(diǎn),并接入供水企業(yè)的在線監(jiān)測點(diǎn)數(shù)據(jù)。本研究采用基于K-means聚類法和熵值法的二次供水水質(zhì)在線監(jiān)測數(shù)據(jù)評價(jià)方法,技術(shù)路線如圖2所示。

    圖2 二次供水水質(zhì)評價(jià)方法及相關(guān)分析技術(shù)路線圖Fig.2 Water quality assessment methods and related analysis technology roadmap of secondary water supply

    2.1 K-means聚類分析法在水質(zhì)評價(jià)的應(yīng)用

    大數(shù)據(jù)分析的理論核心是數(shù)據(jù)挖掘算法,包括聚類分析、回歸分析、因子分析等。聚類算法是常用的大數(shù)據(jù)分析方法,最大應(yīng)用優(yōu)勢在于能夠簡單快速地處理大數(shù)據(jù),并可自主設(shè)定初始聚類中心K,通常選用常規(guī)的歐幾里得距離作為相似度指標(biāo)[14],如式(1)所示。

    (1)

    式中:xi、yi分別為不同樣本x、y的N維數(shù)值;N為二次供水監(jiān)測指標(biāo)數(shù),此處為2維。

    現(xiàn)以L市二次供水在線監(jiān)測結(jié)果為數(shù)據(jù)樣本(監(jiān)測指標(biāo)為余氯和濁度),對比設(shè)定聚類中心分別為2、3時(shí)的兩種聚類結(jié)果(k=2,k=3):

    1)當(dāng)設(shè)定k=2類時(shí),余氯對水質(zhì)分類的影響較大,分界線約為0.35 mg/L,工程經(jīng)驗(yàn)表明,針對小區(qū)內(nèi)而言,余氯值為0.05~0.30 mg/L,較為理想。

    2)當(dāng)設(shè)定k=3類時(shí),濁度與余氯共同作用水質(zhì)分類,二次供水水質(zhì)聚類分析結(jié)果如圖3所示。由圖3可知,目前存在的3類二次供水水質(zhì)特點(diǎn):Ⅰ類水的余氯值適中且濁度控制最好,Ⅰ類水樣本占比為79.1%;Ⅱ類水的余氯值相對較低(處于合適的范圍),同時(shí)濁度適中,Ⅱ類水樣本占比為14.5%;Ⅲ類水的余氯值相對較高(處于合適的范圍),同時(shí)濁度普遍較高,且水質(zhì)季節(jié)性變化最明顯,Ⅲ類水樣本占比為6.5%。濁度和余氯會(huì)共同影響水質(zhì)分類,Ⅰ類、Ⅱ類和Ⅱ類、Ⅲ類水濁度的界限分別為0.2 NTU左右、0.4 NTU左右,Ⅰ類、Ⅱ類和Ⅱ類、Ⅲ類水余氯的界限分別為0.3 mg/L和0.8 mg/L。

    圖3 二次供水水質(zhì)聚類分析結(jié)果(k=3)Fig.3 Cluster analysis results of water supply of secondary water supply (k=3)

    2.2 水質(zhì)評分劃分標(biāo)準(zhǔn)

    根據(jù)上述的兩種聚類中心(k=2,3)的分析結(jié)果,并結(jié)合實(shí)際經(jīng)驗(yàn),人為確定了水質(zhì)評價(jià)的劃分標(biāo)準(zhǔn),具體如表2所示。并依據(jù)水質(zhì)評價(jià)的劃分標(biāo)準(zhǔn)和所有監(jiān)測點(diǎn)的檢測平均值進(jìn)行評分,初步得出計(jì)算結(jié)果。

    表2 水質(zhì)評價(jià)劃分標(biāo)準(zhǔn)Table 2 Criteria for water quality assessment

    2.3 熵值法細(xì)化同類型水質(zhì)評價(jià)

    熵值法通過各項(xiàng)評價(jià)指標(biāo)觀測值所提供信息量的大小來確定評價(jià)指標(biāo)權(quán)重系數(shù),也可用作指標(biāo)離散程度的評判方式,離散值越大代表該指標(biāo)對綜合評價(jià)的影響越大[15]。

    通過熵值法細(xì)化上述劃分的4類水質(zhì)可較好地對數(shù)量較多的監(jiān)測點(diǎn)進(jìn)行后續(xù)分析,規(guī)定熵值法得分(S波動(dòng))不超過1從而保證表2中劃分的水質(zhì)區(qū)別,將表2計(jì)算所得的初步得分結(jié)合熵值法,可得出各點(diǎn)的水質(zhì)波動(dòng)得分S余氯和S濁度,代入式(2)和式(3)得出最終水質(zhì)得分結(jié)果S得分。

    S波動(dòng)=S余氯+S濁度

    (2)

    S得分=S標(biāo)準(zhǔn)-S波動(dòng)

    (3)

    2.4 單因素/多因素相關(guān)性分析

    分別采用單因素方差分析(one-way analysis of variance, ANOVA)和交叉列聯(lián)表法,利用LSD(least significant difference)、皮爾遜數(shù)等顯著性計(jì)算方法,對水源地、水處理工藝、供水方式和管網(wǎng)年限等4個(gè)因素與二次供水水質(zhì)得分單獨(dú)和綜合的相關(guān)性分析,當(dāng)p<0.05時(shí),單因素與水質(zhì)得分存在相關(guān)性,具體結(jié)果如表3所示。表4為主體間效應(yīng)檢驗(yàn)。

    表3 單因素與水質(zhì)得分的相關(guān)性分析

    表4 主體間效應(yīng)檢驗(yàn)Table 4 Intersubjective effect test

    由表3和表4可知:

    1)供水方式將直接影響監(jiān)測點(diǎn)的濁度和余氯指標(biāo)(p<0.05),水源地也存在直接影響監(jiān)測點(diǎn)二次供水水質(zhì)的可能(p=0.07)。盡管水處理工藝和管網(wǎng)使用年限并未直接影響二次供水水質(zhì),但二次供水方式和管網(wǎng)使用年限與水質(zhì)存在協(xié)同效應(yīng),共同影響二次供水水質(zhì)。

    2)結(jié)合4種因素的綜合影響上看,按照影響比重從大到小排序分別為供水方式、管網(wǎng)使用年限和水源地共同影響二次供水水質(zhì)情況,水處理工藝并不能對二次供水的濁度和余氯造成直接或間接影響。

    2.5 協(xié)方差分析

    將2020年1月、2月二次供水水質(zhì)得分協(xié)方差分析法與2019年1月、2月二次供水水質(zhì)得分進(jìn)行對比(即分析相隔一年的同月份二次供水水質(zhì)得分,當(dāng)p<0.05時(shí),存在相關(guān)性)。分析結(jié)果表明4種因素導(dǎo)致水質(zhì)在2年內(nèi)水質(zhì)變化的影響十分有限,因此,水質(zhì)的長期監(jiān)測顯得更有意義。

    結(jié)合第2.4小節(jié)相關(guān)性分析和協(xié)方差的結(jié)果可知,2年內(nèi)水質(zhì)的變化與以上4類因素不具備相關(guān)性,其中管網(wǎng)使用年限的顯著性最小(p=0.19,但是遠(yuǎn)大于0.05)。從協(xié)同影響角度上看,管網(wǎng)使用年限和水源地的共同作用最有可能對逐年二次供水水質(zhì)產(chǎn)生影響(p=0.12,大于0.05),但并不能認(rèn)為兩者間存在相關(guān)性;其次為供水方式和水源地的協(xié)同影響(p=0.17,大于0.05),因此,水源地對二次供水水質(zhì)的間接影響也不可忽視。

    2.6 指標(biāo)散度分析

    箱形圖是一種顯示數(shù)據(jù)分散情況資料的統(tǒng)計(jì)圖,主要用于反映原始數(shù)據(jù)分布的特征,并進(jìn)行多組數(shù)據(jù)分布特征的比較,其數(shù)據(jù)上邊緣和下邊緣分別代表上、下兩個(gè)四分位數(shù),箱體中間為中位數(shù)。將余氯、濁度兩個(gè)指標(biāo)以月份為組別進(jìn)行對比,以觀察供水環(huán)節(jié)對水質(zhì)分散度的影響。L市12個(gè)月份總體余氯水平如圖4所示,L市12個(gè)月份總體濁度水平如圖5所示。

    圖5 L市12個(gè)月總體濁度水平Fig.5 Overall turbidity level of L city in 12 months

    由圖4的余氯散度水平可知:L市12個(gè)月余氯波動(dòng)為0.25~0.75 mg/L,異常值較少,總體余氯水平良好;其中6月、7月的余氯合格率最低,但也達(dá)到95%以上;普遍溫度高的月份余氯水平低于其他月份,但從散度分布上看差別不大,這表明二次供水環(huán)節(jié)中對余氯的管控較為到位,基本可以滿足余氯標(biāo)準(zhǔn)。

    由圖5的濁度散度水平可知:L市12個(gè)月總體濁度水平良好,各月份濁度普遍低于0.2 NTU,合格率幾乎達(dá)到96%以上,各月濁度水平和濁度散度水平相差不大,但異常值較多,這表明部分點(diǎn)的供水設(shè)施需要管理,某些檢測點(diǎn)需要重點(diǎn)關(guān)注。

    3 結(jié)論與建議

    1)不同水質(zhì)在線監(jiān)測儀器對應(yīng)的數(shù)據(jù)質(zhì)量有所不同,本研究發(fā)現(xiàn)濁度儀的數(shù)據(jù)質(zhì)量明顯高于余氯儀。以間斷數(shù)分鐘型為主的時(shí)間間斷型數(shù)據(jù)錯(cuò)誤是普遍存在的(占比約1.5%),區(qū)別于余氯儀和濁度儀的主要數(shù)據(jù)錯(cuò)誤類型為記錄缺失型。

    2)供水方式和水源地改造工程應(yīng)當(dāng)作為二次供水水質(zhì)提升工作的方向,其中供水方式的改造工程更為重要。目前并未發(fā)現(xiàn)城市采用的水處理工藝(主要為“常規(guī)處理+深度處理或深度處理)與二次供水水質(zhì)存在明顯的直接或間接聯(lián)系,說明二次供水水質(zhì)目前不受水處理工藝影響。

    3)供水方式、水源地、水處理工藝和管網(wǎng)使用年限均單獨(dú)對逐年二次供水水質(zhì)變化的影響有限,但水源地的影響不可忽略,因此未來可以從供水方式、水處理工藝和管網(wǎng)使用年限3個(gè)方面保證二次供水水質(zhì)穩(wěn)定性,其中管網(wǎng)運(yùn)維管理是保障逐年水質(zhì)穩(wěn)定性的關(guān)鍵。從水質(zhì)散度分析上看,目前供水企業(yè)對余氯指標(biāo)的管控較為到位,而濁度指標(biāo)應(yīng)是未來重點(diǎn)優(yōu)化的對象,尤其應(yīng)關(guān)注少數(shù)濁度異常的檢測點(diǎn)。

    此外,建議除了余氯和濁度,可適當(dāng)增加高錳酸鹽指數(shù)(CODMn)、pH等少數(shù)檢測簡單、代表性強(qiáng)的水質(zhì)指標(biāo)作為二次供水在線監(jiān)測對象。

    猜你喜歡
    水質(zhì)設(shè)備分析
    水質(zhì)抽檢豈容造假
    環(huán)境(2023年5期)2023-06-30 01:20:01
    諧響應(yīng)分析在設(shè)備減振中的應(yīng)用
    隱蔽失效適航要求符合性驗(yàn)證分析
    一月冬棚養(yǎng)蝦常見水質(zhì)渾濁,要如何解決?這9大原因及處理方法你要知曉
    電力系統(tǒng)不平衡分析
    電子制作(2018年18期)2018-11-14 01:48:24
    基于MPU6050簡單控制設(shè)備
    電子制作(2018年11期)2018-08-04 03:26:08
    電力系統(tǒng)及其自動(dòng)化發(fā)展趨勢分析
    500kV輸變電設(shè)備運(yùn)行維護(hù)探討
    原來他們都是可穿戴設(shè)備
    水質(zhì)總磷測定存在的問題初探
    河南科技(2014年23期)2014-02-27 14:19:07
    bbb黄色大片| 人成视频在线观看免费观看| 啪啪无遮挡十八禁网站| 国产精品美女特级片免费视频播放器 | 亚洲国产精品合色在线| 久久久国产欧美日韩av| 天天添夜夜摸| 欧美色欧美亚洲另类二区| x7x7x7水蜜桃| 99国产精品99久久久久| 神马国产精品三级电影在线观看 | 俺也久久电影网| 国产精品久久电影中文字幕| 又黄又粗又硬又大视频| 美女大奶头视频| av在线天堂中文字幕| 国产成人一区二区三区免费视频网站| 91在线观看av| 亚洲av成人不卡在线观看播放网| 变态另类丝袜制服| 亚洲一卡2卡3卡4卡5卡精品中文| 真人做人爱边吃奶动态| 神马国产精品三级电影在线观看 | 欧美人与性动交α欧美精品济南到| 正在播放国产对白刺激| 久久国产精品影院| 88av欧美| 操出白浆在线播放| 麻豆一二三区av精品| 欧美亚洲日本最大视频资源| 国产一区二区在线av高清观看| 香蕉国产在线看| 久久人妻av系列| 亚洲五月色婷婷综合| 男女下面进入的视频免费午夜 | 人人妻人人看人人澡| 国产精品综合久久久久久久免费| 90打野战视频偷拍视频| 在线av久久热| 欧美不卡视频在线免费观看 | 99久久无色码亚洲精品果冻| 午夜成年电影在线免费观看| 白带黄色成豆腐渣| 久久精品aⅴ一区二区三区四区| 满18在线观看网站| 法律面前人人平等表现在哪些方面| 色老头精品视频在线观看| 97人妻精品一区二区三区麻豆 | 波多野结衣巨乳人妻| 午夜福利高清视频| 一个人免费在线观看的高清视频| aaaaa片日本免费| 国产成人精品久久二区二区免费| 韩国精品一区二区三区| 免费女性裸体啪啪无遮挡网站| 两人在一起打扑克的视频| 黑人操中国人逼视频| 神马国产精品三级电影在线观看 | 国产精品爽爽va在线观看网站 | 国产精品永久免费网站| 男男h啪啪无遮挡| 一本大道久久a久久精品| 美女国产高潮福利片在线看| 999久久久国产精品视频| 人人妻人人澡欧美一区二区| 亚洲第一青青草原| 亚洲av第一区精品v没综合| 母亲3免费完整高清在线观看| 久久久久九九精品影院| 免费在线观看成人毛片| 成人18禁高潮啪啪吃奶动态图| 久久青草综合色| 99久久99久久久精品蜜桃| 欧美绝顶高潮抽搐喷水| 亚洲精品中文字幕在线视频| 久久国产精品男人的天堂亚洲| 在线永久观看黄色视频| 俄罗斯特黄特色一大片| 青草久久国产| 99久久综合精品五月天人人| 窝窝影院91人妻| 无限看片的www在线观看| 成人国产综合亚洲| 久久久久亚洲av毛片大全| 美女免费视频网站| 一个人免费在线观看的高清视频| 性欧美人与动物交配| 国产伦在线观看视频一区| 波多野结衣av一区二区av| 国产亚洲精品久久久久5区| 亚洲一区高清亚洲精品| 欧美精品亚洲一区二区| 久久久久久久久久黄片| 非洲黑人性xxxx精品又粗又长| 色综合婷婷激情| 日本五十路高清| 淫秽高清视频在线观看| 校园春色视频在线观看| 夜夜看夜夜爽夜夜摸| 99精品欧美一区二区三区四区| 1024手机看黄色片| 欧美乱码精品一区二区三区| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品野战在线观看| 免费看a级黄色片| 日韩欧美免费精品| 亚洲欧洲精品一区二区精品久久久| 国产亚洲精品久久久久久毛片| 欧美激情 高清一区二区三区| 亚洲五月婷婷丁香| 日日摸夜夜添夜夜添小说| 韩国精品一区二区三区| 男女视频在线观看网站免费 | 亚洲 欧美 日韩 在线 免费| 久久久精品欧美日韩精品| 亚洲人成伊人成综合网2020| 免费无遮挡裸体视频| 精品国产美女av久久久久小说| 桃红色精品国产亚洲av| 夜夜躁狠狠躁天天躁| 午夜视频精品福利| or卡值多少钱| 在线播放国产精品三级| 欧美亚洲日本最大视频资源| 久久人妻福利社区极品人妻图片| 国产97色在线日韩免费| 国产精品亚洲一级av第二区| a级毛片a级免费在线| 亚洲一区高清亚洲精品| 亚洲成人免费电影在线观看| 琪琪午夜伦伦电影理论片6080| av片东京热男人的天堂| 曰老女人黄片| 欧美激情极品国产一区二区三区| 精品日产1卡2卡| 婷婷精品国产亚洲av| 88av欧美| www国产在线视频色| 欧美日韩精品网址| 国产精品九九99| 国产成人影院久久av| 亚洲人成伊人成综合网2020| 人妻丰满熟妇av一区二区三区| 免费在线观看视频国产中文字幕亚洲| 中文亚洲av片在线观看爽| 精品久久久久久,| 国产亚洲欧美在线一区二区| 高清毛片免费观看视频网站| 午夜久久久在线观看| 国产精品亚洲美女久久久| 成人特级黄色片久久久久久久| 美女午夜性视频免费| a级毛片a级免费在线| 麻豆成人午夜福利视频| 久久久久久久久免费视频了| 曰老女人黄片| av中文乱码字幕在线| 国产三级黄色录像| 亚洲成国产人片在线观看| 亚洲一区中文字幕在线| 美女高潮喷水抽搐中文字幕| 亚洲一区高清亚洲精品| 麻豆成人av在线观看| 男女床上黄色一级片免费看| 国产精品亚洲av一区麻豆| av电影中文网址| 午夜福利18| 亚洲av第一区精品v没综合| 欧美中文综合在线视频| 91字幕亚洲| а√天堂www在线а√下载| 欧洲精品卡2卡3卡4卡5卡区| 此物有八面人人有两片| 日本 欧美在线| 两个人看的免费小视频| 国产野战对白在线观看| 中亚洲国语对白在线视频| 久久青草综合色| 无遮挡黄片免费观看| 国产精品久久视频播放| 黑人操中国人逼视频| 久久精品aⅴ一区二区三区四区| 男人操女人黄网站| 国产免费男女视频| 亚洲精品美女久久久久99蜜臀| 非洲黑人性xxxx精品又粗又长| 亚洲一区中文字幕在线| 色综合亚洲欧美另类图片| 99国产精品一区二区蜜桃av| 亚洲成a人片在线一区二区| 黑人操中国人逼视频| 久热这里只有精品99| 麻豆国产av国片精品| 国产伦在线观看视频一区| 欧美午夜高清在线| www国产在线视频色| 熟女电影av网| 精品免费久久久久久久清纯| 午夜免费鲁丝| 日韩精品免费视频一区二区三区| 1024手机看黄色片| 亚洲精品一卡2卡三卡4卡5卡| 国产91精品成人一区二区三区| 俺也久久电影网| 最近最新中文字幕大全免费视频| 制服丝袜大香蕉在线| 久久精品aⅴ一区二区三区四区| 国内少妇人妻偷人精品xxx网站 | 欧美最黄视频在线播放免费| 亚洲 国产 在线| 欧美激情高清一区二区三区| 亚洲国产中文字幕在线视频| 欧美日韩瑟瑟在线播放| www日本黄色视频网| 欧美zozozo另类| 国产久久久一区二区三区| 搞女人的毛片| 在线视频色国产色| 免费搜索国产男女视频| 日本三级黄在线观看| 国产成+人综合+亚洲专区| 在线视频色国产色| а√天堂www在线а√下载| 欧美成人免费av一区二区三区| 欧美激情久久久久久爽电影| 国产精品久久久久久精品电影 | 国语自产精品视频在线第100页| 精品久久久久久,| 丝袜美腿诱惑在线| 在线观看日韩欧美| 亚洲在线自拍视频| 淫妇啪啪啪对白视频| 日日摸夜夜添夜夜添小说| 色哟哟哟哟哟哟| 欧美日韩亚洲综合一区二区三区_| 男人操女人黄网站| 欧洲精品卡2卡3卡4卡5卡区| 欧美乱色亚洲激情| 波多野结衣高清无吗| 哪里可以看免费的av片| 亚洲一码二码三码区别大吗| 久久久久精品国产欧美久久久| 免费在线观看影片大全网站| 欧美黑人巨大hd| 国产1区2区3区精品| 男人舔女人下体高潮全视频| xxxwww97欧美| 亚洲国产欧美日韩在线播放| 精品久久久久久久末码| 国产成人精品久久二区二区免费| 亚洲av美国av| 在线av久久热| 妹子高潮喷水视频| 99久久无色码亚洲精品果冻| 黑人操中国人逼视频| 欧美黄色片欧美黄色片| 午夜久久久在线观看| 日本免费一区二区三区高清不卡| 丝袜人妻中文字幕| 久久精品亚洲精品国产色婷小说| 老汉色av国产亚洲站长工具| 桃色一区二区三区在线观看| 精品一区二区三区视频在线观看免费| 久久精品91无色码中文字幕| 757午夜福利合集在线观看| 美女高潮到喷水免费观看| 天堂影院成人在线观看| 日韩欧美一区视频在线观看| 亚洲第一欧美日韩一区二区三区| 国产高清视频在线播放一区| 在线观看午夜福利视频| av免费在线观看网站| 久久人妻福利社区极品人妻图片| 9191精品国产免费久久| 琪琪午夜伦伦电影理论片6080| 日韩欧美国产在线观看| 在线观看一区二区三区| 精品欧美一区二区三区在线| 亚洲欧美日韩无卡精品| 老司机午夜十八禁免费视频| 国产午夜福利久久久久久| 亚洲一码二码三码区别大吗| 成人永久免费在线观看视频| 午夜福利免费观看在线| 国产激情偷乱视频一区二区| 老汉色av国产亚洲站长工具| 国内精品久久久久精免费| 男女做爰动态图高潮gif福利片| 十分钟在线观看高清视频www| 少妇粗大呻吟视频| 一夜夜www| 久久久久久久久免费视频了| 国产av在哪里看| 欧美zozozo另类| 一级毛片精品| 亚洲精品美女久久av网站| 美女扒开内裤让男人捅视频| 美女午夜性视频免费| 两人在一起打扑克的视频| 国产精品综合久久久久久久免费| 麻豆国产av国片精品| 一边摸一边做爽爽视频免费| 亚洲精品久久国产高清桃花| 美女国产高潮福利片在线看| 黄频高清免费视频| 熟妇人妻久久中文字幕3abv| 久久中文字幕一级| 欧美日本视频| 国产精品自产拍在线观看55亚洲| 一级毛片精品| 国产精品日韩av在线免费观看| 久久热在线av| 少妇被粗大的猛进出69影院| 中文字幕人妻熟女乱码| 国产精品亚洲av一区麻豆| 最近最新中文字幕大全电影3 | 亚洲美女黄片视频| 久久久久久久精品吃奶| 色尼玛亚洲综合影院| 国内少妇人妻偷人精品xxx网站 | 国产高清视频在线播放一区| 亚洲九九香蕉| 午夜成年电影在线免费观看| 亚洲成人久久性| 性色av乱码一区二区三区2| 亚洲精品国产一区二区精华液| 欧美精品啪啪一区二区三区| 两个人免费观看高清视频| 青草久久国产| 在线观看www视频免费| 1024视频免费在线观看| 又紧又爽又黄一区二区| 在线观看午夜福利视频| 久热爱精品视频在线9| 免费搜索国产男女视频| 亚洲专区国产一区二区| 俺也久久电影网| 亚洲专区国产一区二区| 这个男人来自地球电影免费观看| 一边摸一边做爽爽视频免费| 在线观看免费日韩欧美大片| 亚洲一区二区三区色噜噜| 嫁个100分男人电影在线观看| 精品不卡国产一区二区三区| 一本久久中文字幕| 国产一级毛片七仙女欲春2 | 国产男靠女视频免费网站| 99热6这里只有精品| 久久婷婷成人综合色麻豆| 国产精品精品国产色婷婷| 黄片播放在线免费| 日韩欧美三级三区| 精品国产乱子伦一区二区三区| 国产免费av片在线观看野外av| 非洲黑人性xxxx精品又粗又长| 精品福利观看| 久久久久久久精品吃奶| 又黄又爽又免费观看的视频| 久久久久久久午夜电影| 两个人看的免费小视频| 国产av不卡久久| 亚洲国产精品久久男人天堂| 亚洲全国av大片| 日韩精品中文字幕看吧| 久久天堂一区二区三区四区| 人人妻人人澡人人看| 日韩精品中文字幕看吧| 欧美黑人巨大hd| 国产黄片美女视频| 精品久久久久久久末码| 亚洲国产高清在线一区二区三 | 在线播放国产精品三级| 少妇被粗大的猛进出69影院| 成人18禁高潮啪啪吃奶动态图| 亚洲成人国产一区在线观看| 人妻丰满熟妇av一区二区三区| 亚洲熟妇中文字幕五十中出| 亚洲 欧美一区二区三区| АⅤ资源中文在线天堂| 麻豆一二三区av精品| 国产精品九九99| 日日摸夜夜添夜夜添小说| 亚洲精品国产精品久久久不卡| 免费一级毛片在线播放高清视频| 国产精品自产拍在线观看55亚洲| 老司机深夜福利视频在线观看| 51午夜福利影视在线观看| 欧美性长视频在线观看| 老鸭窝网址在线观看| 午夜激情av网站| 男女视频在线观看网站免费 | 国产91精品成人一区二区三区| 亚洲av电影在线进入| 午夜久久久久精精品| 成人三级黄色视频| 国产精品久久视频播放| 国产成+人综合+亚洲专区| www.999成人在线观看| 老司机福利观看| 精品电影一区二区在线| 高清毛片免费观看视频网站| 精品一区二区三区视频在线观看免费| 亚洲人成电影免费在线| 老熟妇仑乱视频hdxx| 18禁黄网站禁片午夜丰满| 18禁观看日本| 夜夜夜夜夜久久久久| 午夜影院日韩av| 久久草成人影院| 一a级毛片在线观看| √禁漫天堂资源中文www| 男女做爰动态图高潮gif福利片| 日本在线视频免费播放| 国产精品,欧美在线| 最新美女视频免费是黄的| 亚洲熟妇中文字幕五十中出| 一级毛片精品| 91av网站免费观看| 国产精品98久久久久久宅男小说| 日韩高清综合在线| 老汉色∧v一级毛片| 免费在线观看黄色视频的| 国产精品影院久久| 精品日产1卡2卡| 国产蜜桃级精品一区二区三区| a级毛片在线看网站| 最近最新中文字幕大全电影3 | 国产精品自产拍在线观看55亚洲| 黄片播放在线免费| 草草在线视频免费看| 久久精品国产99精品国产亚洲性色| 国内少妇人妻偷人精品xxx网站 | 国产又爽黄色视频| 99久久国产精品久久久| 日本黄色视频三级网站网址| 国产成人精品无人区| 香蕉久久夜色| 哪里可以看免费的av片| 国产精品永久免费网站| 久久国产乱子伦精品免费另类| 99在线视频只有这里精品首页| 久久婷婷成人综合色麻豆| 国产精品美女特级片免费视频播放器 | 1024手机看黄色片| 午夜免费成人在线视频| 亚洲av电影不卡..在线观看| 精品国产美女av久久久久小说| 精品一区二区三区av网在线观看| 女性生殖器流出的白浆| 在线十欧美十亚洲十日本专区| 男女视频在线观看网站免费 | 99热只有精品国产| 久久 成人 亚洲| 国产精品av久久久久免费| 婷婷亚洲欧美| 国产色视频综合| 成人国产一区最新在线观看| 婷婷丁香在线五月| 精品乱码久久久久久99久播| 国产1区2区3区精品| 欧美日韩乱码在线| 国产亚洲av高清不卡| 国产熟女午夜一区二区三区| 亚洲成人精品中文字幕电影| 久久伊人香网站| 亚洲久久久国产精品| 亚洲av片天天在线观看| 亚洲黑人精品在线| 999久久久精品免费观看国产| 热99re8久久精品国产| 亚洲国产毛片av蜜桃av| 伦理电影免费视频| 色综合站精品国产| www日本在线高清视频| 视频在线观看一区二区三区| 午夜视频精品福利| 亚洲欧美精品综合久久99| 国产欧美日韩一区二区三| 精品久久久久久久久久免费视频| 免费高清在线观看日韩| 777久久人妻少妇嫩草av网站| 免费搜索国产男女视频| 国产精品一区二区精品视频观看| www.www免费av| 国产精品综合久久久久久久免费| 欧美一级毛片孕妇| 中文字幕人妻丝袜一区二区| 国产伦一二天堂av在线观看| 在线观看午夜福利视频| 亚洲av成人不卡在线观看播放网| 法律面前人人平等表现在哪些方面| 亚洲精品美女久久久久99蜜臀| 91在线观看av| 亚洲精品久久国产高清桃花| 亚洲全国av大片| 色av中文字幕| 亚洲精品美女久久av网站| 熟妇人妻久久中文字幕3abv| 一区二区三区精品91| 日韩欧美国产一区二区入口| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品美女久久av网站| 无限看片的www在线观看| 日韩大码丰满熟妇| 亚洲av成人不卡在线观看播放网| 日本在线视频免费播放| 久久久久久久午夜电影| 激情在线观看视频在线高清| 麻豆国产av国片精品| 天天躁狠狠躁夜夜躁狠狠躁| 人人澡人人妻人| 99国产综合亚洲精品| 在线观看免费视频日本深夜| 老司机午夜十八禁免费视频| 一进一出好大好爽视频| 91在线观看av| 国产成人影院久久av| 午夜福利在线在线| 日韩欧美国产在线观看| 国产精品精品国产色婷婷| 日本成人三级电影网站| 非洲黑人性xxxx精品又粗又长| 欧美黑人精品巨大| 国产1区2区3区精品| 国产精品久久久久久人妻精品电影| 国产午夜福利久久久久久| xxxwww97欧美| 国产私拍福利视频在线观看| 中文资源天堂在线| a在线观看视频网站| 此物有八面人人有两片| 精品久久久久久久久久久久久 | 欧洲精品卡2卡3卡4卡5卡区| 亚洲熟妇熟女久久| 欧美激情久久久久久爽电影| 桃红色精品国产亚洲av| 欧美av亚洲av综合av国产av| 亚洲人成网站高清观看| 亚洲av第一区精品v没综合| 精品一区二区三区视频在线观看免费| 美女 人体艺术 gogo| 国产主播在线观看一区二区| 久久国产亚洲av麻豆专区| 69av精品久久久久久| 妹子高潮喷水视频| 麻豆国产av国片精品| 身体一侧抽搐| 免费在线观看黄色视频的| 精品国产亚洲在线| 男人操女人黄网站| 亚洲 国产 在线| 欧美日韩一级在线毛片| 成人国产一区最新在线观看| 99久久精品国产亚洲精品| 国产午夜福利久久久久久| 黄片大片在线免费观看| 在线观看www视频免费| 成年人黄色毛片网站| 中文亚洲av片在线观看爽| 久久中文字幕一级| 久久精品亚洲精品国产色婷小说| 少妇裸体淫交视频免费看高清 | 老司机午夜福利在线观看视频| 巨乳人妻的诱惑在线观看| 久9热在线精品视频| 精品国产超薄肉色丝袜足j| 亚洲色图 男人天堂 中文字幕| 每晚都被弄得嗷嗷叫到高潮| 久9热在线精品视频| 黄网站色视频无遮挡免费观看| 精品欧美国产一区二区三| 韩国精品一区二区三区| 精品一区二区三区四区五区乱码| 男人舔奶头视频| 国产精品亚洲av一区麻豆| 日本黄色视频三级网站网址| 最近在线观看免费完整版| 国产高清视频在线播放一区| 久久九九热精品免费| 天天躁夜夜躁狠狠躁躁| 伦理电影免费视频| 国产伦人伦偷精品视频| 日韩欧美 国产精品| 日本三级黄在线观看| 成在线人永久免费视频| 国产99白浆流出| 午夜福利成人在线免费观看| 女性生殖器流出的白浆| 熟女电影av网| www国产在线视频色| 国产一区二区在线av高清观看| 黄色成人免费大全| 国内精品久久久久久久电影| 亚洲精品色激情综合| 成人国产综合亚洲| 免费在线观看成人毛片| 亚洲成av片中文字幕在线观看| 成人国产综合亚洲| 在线看三级毛片| 国产精品免费一区二区三区在线| 99国产精品一区二区三区| 欧美日韩福利视频一区二区| 久久精品国产亚洲av香蕉五月| 国产伦在线观看视频一区| 亚洲一卡2卡3卡4卡5卡精品中文| 久久99热这里只有精品18| 成人永久免费在线观看视频| 99精品久久久久人妻精品| 亚洲激情在线av| 欧美精品啪啪一区二区三区| 国产成人一区二区三区免费视频网站| 99热只有精品国产| 白带黄色成豆腐渣| 久久久久久亚洲精品国产蜜桃av|