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

    基于PSO-DE算法的突發(fā)水域污染溯源研究

    2017-11-07 04:47:51曹宏桂贠衛(wèi)國
    中國環(huán)境科學(xué) 2017年10期
    關(guān)鍵詞:優(yōu)化

    曹宏桂,贠衛(wèi)國

    ?

    基于PSO-DE算法的突發(fā)水域污染溯源研究

    曹宏桂,贠衛(wèi)國*

    (西安建筑科技大學(xué)信息與控制工程學(xué)院,陜西西安710055)

    利用PSO-DE混合優(yōu)化算法結(jié)合移動(dòng)監(jiān)測平臺(tái)研究了污染物源項(xiàng)識別問題,包括單點(diǎn)固定源和多點(diǎn)固定源位置的反演.該方法把源項(xiàng)識別反問題轉(zhuǎn)化為非線性優(yōu)化問題,用個(gè)移動(dòng)平臺(tái)檢測并記錄所在水域的污染物濃度,將各自位置的坐標(biāo)值記為此移動(dòng)平臺(tái)的best,每一個(gè)移動(dòng)平臺(tái)均對應(yīng)一個(gè)best,即共有個(gè)best,將個(gè)移動(dòng)平臺(tái)獲取的污染物濃度值進(jìn)行對比,選擇最大污染物濃度值對應(yīng)的水域坐標(biāo),記為best,以此作為初始種群先進(jìn)行PSO優(yōu)化獲得的種群,再進(jìn)行DE優(yōu)化,取兩者濃度高的作為best,直到獲得濃度值最高的點(diǎn),即污染物初始投放點(diǎn).多個(gè)算例的計(jì)算結(jié)果表明,采用該算法對含點(diǎn)源的二維水域污染源溯源問題能夠得到精度較高的反演結(jié)果.

    PSO-DE;污染物溯源;移動(dòng)監(jiān)測平臺(tái);二維水域

    污染物源項(xiàng)識別反問題是指根據(jù)環(huán)境水力學(xué)測量數(shù)據(jù),來反推污染源信息的問題,如獲取點(diǎn)源的個(gè)數(shù)、強(qiáng)度和位置等[1],是環(huán)境水力學(xué)反問題中一類重要的反問題,具有重要的研究價(jià)值和應(yīng)用背景.目前相關(guān)方面的研究很多但尚不成熟,根據(jù)研究方法的特點(diǎn),可以分為三類[2-3],第一類方法將污染源項(xiàng)識別問題轉(zhuǎn)化為概率計(jì)算問題,利用概率和數(shù)理統(tǒng)計(jì)理論估算污染源的源項(xiàng)信息[4-5],如陳媛華等[4]使用相關(guān)系數(shù)法進(jìn)行污染物源項(xiàng)信息識別,王家彪等[5]使用耦合概率密度溯源河渠突發(fā)水污染,朱嵩等[6]基于貝葉斯推理實(shí)現(xiàn)標(biāo)準(zhǔn)k-e湍流模型的參數(shù)識別,楊海東等[7]在微分進(jìn)化和蒙特卡羅的基礎(chǔ)上利用貝葉斯算法改進(jìn),提高其模型的穩(wěn)定性和可靠性.但這類方法對于模型的準(zhǔn)確性有非常高的要求,而由于環(huán)境水力模型在建模時(shí)有很多條件都是理想化的,所以最終的計(jì)算精度誤差較大.鑒于這種情況,第二類方法采用對模型準(zhǔn)確性不高的智能尋優(yōu)算法,如辛小康等[8]采用遺傳算法對水污染事故的源項(xiàng)信息進(jìn)行識別,但是這種方法需要有足夠的訓(xùn)練樣本對模型進(jìn)行訓(xùn)練,對于突發(fā)水域污染的情況無法獲得足夠數(shù)量的樣本,所以可行性不高.第三類方法借鑒地下水污染的方法,將污染源溯源問題轉(zhuǎn)化為尋找有限元的最優(yōu)解問題[9-11],使用DE差分進(jìn)化算法和PSO粒子群算法進(jìn)行全局最優(yōu)解的尋找.如牟行洋等[10], Jong等[12], Ednah等[13]都有類似的研究,但由于算法本身的特點(diǎn),粒子群算法擁有較好的全局搜索能力,差分進(jìn)化算法具有較強(qiáng)的局部搜索能力,如果可以將兩個(gè)算法的優(yōu)點(diǎn)結(jié)合將會(huì)獲得更精確的尋優(yōu)結(jié)果[12-15].

    本研究采用有限元法作為含點(diǎn)源的二維環(huán)境水力學(xué)模型正演問題的數(shù)值計(jì)算方法,采用PSO-DE粒子群與差進(jìn)化混合優(yōu)化算法結(jié)合小型移動(dòng)監(jiān)測平臺(tái)進(jìn)行相應(yīng)的參數(shù)反演研究,包括單點(diǎn)固定源反演和多點(diǎn)固定源反演.

    1 數(shù)學(xué)模型

    1.1 二維水域的單點(diǎn)穩(wěn)態(tài)源污染物水質(zhì)模型

    模型為無限邊際水體假設(shè)下的二維(水域縱向和橫向)穩(wěn)態(tài)點(diǎn)源排放的水質(zhì)模型,假設(shè)當(dāng)污染物投入水中時(shí),在深度方向迅速擴(kuò)散均衡.由于污染物是穩(wěn)態(tài)點(diǎn)源污染物,既以恒定速率往水中注入污染物,當(dāng)某水域出現(xiàn)污染物后,該水域中污染物的質(zhì)量總量并不隨著時(shí)間而改變,只和與污染源的相對坐標(biāo)有關(guān),與時(shí)間無關(guān).

    二維水域的單點(diǎn)穩(wěn)態(tài)源(恒定速率排放污染物)污染物(惰性污染物)濃度的解析解如式1所示[16].

    1.2 二維無限水域多點(diǎn)穩(wěn)態(tài)源污染物水質(zhì)模型

    在二維的單點(diǎn)源污染模型的基礎(chǔ)上,為了能夠更加貼近實(shí)際的水域情境,將單點(diǎn)源污染物增至多點(diǎn)源污染物,在式1的基礎(chǔ)上,推導(dǎo)出此時(shí)二維多點(diǎn)穩(wěn)態(tài)源(恒定速率排放污染物)污染物水域的污染物濃度解析解為:

    同樣的,由于此時(shí)仍舊為穩(wěn)態(tài)源排放污染物,因此水域中污染物的質(zhì)量總量為恒定值,與時(shí)間無關(guān).式2可以用來模擬大型二維水域同時(shí)存在多個(gè)同類點(diǎn)源污染物穩(wěn)定排污情景下的水域各個(gè)斷面的污染物濃度.

    2 PSO-DE粒子群與差進(jìn)化混合優(yōu)化算法

    粒子群算法是由Keimedyl和Eberhart 1995年提出,源于對鳥群捕食行為的研究.用于解決優(yōu)化問題.粒子群進(jìn)化算法雖然收斂速度較快、容易實(shí)現(xiàn),在各種問題的求解和應(yīng)用已經(jīng)展現(xiàn)其獨(dú)特的優(yōu)點(diǎn),但在收斂過程中容易停滯,收斂精度較低.微分進(jìn)化算法作為一種新的算法,近年來才引起廣大學(xué)者的關(guān)注,算法具有易于操作,優(yōu)化性能好等特點(diǎn),但容易陷入局部最優(yōu),特別是到了進(jìn)化后期收斂速度明顯降低[17-18].

    為了改善粒子群進(jìn)化算法和微分進(jìn)化算法獨(dú)自求解優(yōu)化問題的不足.本文引用文獻(xiàn)中的粒子群算法和微分進(jìn)化算法混合的方法.該算法是對種群先后做PSO和DE優(yōu)化的策略,對于一個(gè)種群,先通過PSO優(yōu)化找到最優(yōu)解得大概位置,再運(yùn)用DE優(yōu)化找到最優(yōu)解,為了協(xié)調(diào)全局和局部搜索能力,新算法采取慣性權(quán)重和學(xué)習(xí)因子非線性變化策略;而為了提高算法跳出局部最優(yōu)的能力,在DE優(yōu)化時(shí)采取放縮因子非線性變化策略.

    PSO-DE算法的具體實(shí)現(xiàn)框架如下:

    步驟1:粒子群初始化.初始速度和位置,將每個(gè)小型移動(dòng)平臺(tái)的初始坐標(biāo)賦值為種群的初始值,種群規(guī)模NP即為平臺(tái)個(gè)數(shù),最大迭代次數(shù)itermax,慣性權(quán)重最大值max,慣性權(quán)重最小值min,學(xué)習(xí)因子1、2,變異概率CR,每個(gè)粒子的適應(yīng)度最好的點(diǎn)best被設(shè)置為其初始位置,best中的適應(yīng)度最好的點(diǎn)被設(shè)置成為全局最佳粒子best.

    (1)慣性權(quán)重因子

    慣性權(quán)重是影響算法性能的關(guān)鍵參數(shù),根據(jù)其特點(diǎn),在粒子群迭代的初期因具有較好的全局搜索能力,到達(dá)后期,則更希望粒子具有較好的局部搜索能力[19].提出了以下動(dòng)態(tài)變化的慣性權(quán)重:

    式中:iter和itermax分別是當(dāng)前迭代次數(shù)和最大迭代次數(shù),max和min分別為慣性權(quán)重的初始最大值和最小值,一般取max=0.9,min=0.4.

    (2)時(shí)變學(xué)習(xí)因子

    在粒子群進(jìn)化算法中,學(xué)習(xí)因子1和2分別影響著每個(gè)粒子的個(gè)體最優(yōu)和全局最優(yōu).較小的學(xué)習(xí)因子會(huì)限制粒子的運(yùn)動(dòng),而較大的學(xué)習(xí)因子會(huì)使粒子過于發(fā)散.基于此提出了時(shí)變學(xué)習(xí)因子策略:

    式中:11、12、21和22是常數(shù),iter和itermax分別是當(dāng)前及最大迭代次數(shù).當(dāng)11=21=1.5時(shí),12=22=0.5算法的優(yōu)化性能最好.

    步驟 2: PSO優(yōu)化階段.

    (1):通過下列公式修正速度與位置:

    步驟5:如果滿足結(jié)束條件,即達(dá)到最大迭代次數(shù),則循環(huán)結(jié)束輸出結(jié)果,否則跳轉(zhuǎn)回步驟2.

    3 數(shù)值算例

    3.1 單點(diǎn)源算例分析

    3.1.1 算例來源 假設(shè)某個(gè)大型的湖泊,在位置A排放2000kg惰性污染物,該湖泊的平均水深為l0m,平均水流速度0.5m/s,橫向彌散系數(shù)50m2/ min[20].該水域相關(guān)斷面的污染物濃度可以通過公式(1)來計(jì)算.令位置A為坐標(biāo)原點(diǎn)(0,0),選取水域縱向100m(步長為0.5)橫向200m(步長為1)的范圍繪制污染物濃度分布如圖1所示.

    圖1 二維水域單點(diǎn)穩(wěn)態(tài)源污染物濃度分布

    圖2 單點(diǎn)源最佳粒子運(yùn)動(dòng)軌跡

    表1 單點(diǎn)源算例的最佳粒子坐標(biāo)值

    3.2 多點(diǎn)源算例分析

    3.2.1 算例來源 假設(shè)污染源在單點(diǎn)固定源算例的基礎(chǔ)上增加到3個(gè)點(diǎn)源,即=3.3個(gè)點(diǎn)源的位置分別為(0m,0m)、(50m,-50m)和(60m,60m),3個(gè)點(diǎn)源的污染物質(zhì)量分別為2000,1000和300kg.其他條件保持不變.其濃度分布圖如圖3所示.

    3.2.2 結(jié)果分析 對于多點(diǎn)源的情況,依然搜索全局范圍內(nèi)污染物濃度最大的點(diǎn),即坐標(biāo)值為(0,0)的點(diǎn).最佳粒子的運(yùn)動(dòng)軌跡如圖4所示.其坐標(biāo)值如表2所示.

    圖3 二維水域多點(diǎn)穩(wěn)態(tài)源污染物濃度分布

    圖4 多點(diǎn)源最佳粒子運(yùn)動(dòng)軌跡

    表2 多點(diǎn)源算例的最佳粒子坐標(biāo)值

    通過圖4可以看出,多點(diǎn)源情況下在開始搜索的時(shí)候出現(xiàn)了往復(fù),但最終收斂在污染物濃度最大點(diǎn).通過表2數(shù)據(jù)顯示,在迭代到第40次時(shí)達(dá)到污染物最大的點(diǎn)附近并保持穩(wěn)定.所以PSO- DE優(yōu)化算法在多點(diǎn)源情況下依然可以達(dá)到要求.對于另外的2個(gè)污染源暫時(shí)不能準(zhǔn)確獲取其坐標(biāo),但是其大致范圍基本可以確定在粒子出現(xiàn)多次往復(fù)的區(qū)域.

    3.3 反演結(jié)果影響因子分析

    3.3.1 粒子個(gè)數(shù)的影響 實(shí)際應(yīng)用過程中初始種群的粒子數(shù)關(guān)系到成本問題,因此考察單點(diǎn)源反演過程中不同的粒子數(shù)對反演結(jié)果的影響.

    由表3可知,在粒子數(shù)不同的情況下,本算例中隨著粒子數(shù)的增加,其反演結(jié)果的精度也隨之提高,但是考慮到實(shí)際應(yīng)用中成本問題,結(jié)合反演結(jié)果的誤差,建議實(shí)際應(yīng)用時(shí)粒子數(shù)設(shè)置為5.

    表3 粒子個(gè)數(shù)對反演結(jié)果的影響

    表4 水流速對反演結(jié)果的影響 Table 4 Influence of inversion results for water velocity

    由表4可知,在其他情況相同的情況下,水流的流速越快對反演計(jì)算的影響越大,溯源精度越低.在流速較快的水域可以通過增加粒子數(shù)量來提高反演精度,或者增加迭代次數(shù)來解決.

    4 結(jié)論

    本研究考慮了污染物水質(zhì)特征,給出二維水域下適用于移動(dòng)平臺(tái)的PSO-DE優(yōu)化溯源算法.對PSO-DE算法中的慣性權(quán)重因子和時(shí)變學(xué)習(xí)因子進(jìn)行改進(jìn),通過MATLAB進(jìn)行算例仿真,在單點(diǎn)源污染的情況下,污染源定位坐標(biāo)值為(0.001818,0.3972),與實(shí)際位置誤差僅為0.3972m;多點(diǎn)源污染的情況下,確定了排放量最大的點(diǎn)坐標(biāo),同時(shí)確定了其他污染源的范圍.結(jié)果表明本方法的計(jì)算結(jié)果接近真實(shí)情況.探討和分析論證了所提出的算法在實(shí)際運(yùn)用時(shí)可能遇到的問題,具有一定的應(yīng)用價(jià)值.

    [1] 朱 劍,楊 江,鄧尚超.適用于移動(dòng)平的河流突發(fā)污染算法研究 [J]. 安全與環(huán)境學(xué)報(bào), 2015,15(6):295-299.

    [2] Gurhan G, Halil K. Solving inverse problems of groundwater- pollution-source identification using a differential evolution algorithm [J]. Hydrogeology Journal, 2015,23:1109–1119.

    [3] Jyoti C, Deepak K. Groundwater pollution source identification through inverse modeling [J]. Water Resources and River Engineering, 2015.

    [4] 陳媛華,王 鵬,姜繼平,等.基于相關(guān)系數(shù)優(yōu)化法的河流突發(fā)污染源項(xiàng)識別 [J]. 中國環(huán)境科學(xué), 2011,31(11):1802-1807.

    [5] 王家彪,雷曉輝,廖衛(wèi)紅,等.基于耦合概率密度方法的河渠突發(fā)水污染溯源 [J]. 水利學(xué)報(bào), 2015,46(11):1280-1289.

    [6] 朱 嵩,劉國華,毛欣煒,等.基于貝葉斯推理的標(biāo)準(zhǔn)k-ε湍流模型參數(shù)識別 [J]. 四川大學(xué)學(xué)報(bào)(工程科學(xué)版), 2010,42(4):78- 82.

    [7] 楊海東,肖 宜,王卓民,等.突發(fā)性水污染事件溯源方法[J]. 水科學(xué)進(jìn)展, 2014,25(1):122-129.

    [8] 辛小康,韓小波,李 建,等.基于遺傳算法的水污染事故污染源識別模型 [J]. 水電能源科學(xué), 2014,(7):52-55.

    [9] 葉永偉,葛沈浩,任設(shè)東,等.基于PSO-DE算法的污水處理優(yōu)化控制研究[J]. 計(jì)算機(jī)測量與控制, 2016,24(2):68-76.

    [10] 牟行洋.基于微分進(jìn)化算法的污染物源項(xiàng)識別反問題研究術(shù) [J]. 水動(dòng)力學(xué)研究與進(jìn)展, 2011,26(1):24-30.

    [11] 彭 亮,于 暉.基于濃度梯度-逆流搜索的污染源定位策略 [J]. 智能計(jì)算機(jī)與應(yīng)用, 2012,2(6):64-66.

    [12] Jong C, Meng G W, et al. Structural reliability-based optimization design using PSO-DE hybrid algorithm [J]. Journal of South China University of Technology, 2014,42(9): 41-45.

    [13] Ednah O, Akpofure T, John N. Groundwater pollution source identification by optimization and the green element method [J]. World Environmental and Water Resources Congress, 2016, 309-318.

    [14] MT A. A hybrid simulation–optimization approach for solving the areal groundwater pollution source identification problems [J]. Journal of Hydrology, 2016,538:161–176.

    [15] 肖傳寧,盧文喜,趙 瑩,等.基于徑向基函數(shù)模型的優(yōu)化方法在地下水污染源識別中的應(yīng)用 [J]. 中國環(huán)境科學(xué), 2016,36(7): 2067-2072.

    [16] 傅國偉.河流水質(zhì)數(shù)學(xué)模型及其模擬計(jì)算[M]. 北京:中國環(huán)境科學(xué)出版社, 1987.

    [17] Yu F Y, Chen X, Guo L. Hybrid algorithm based on particle swarm optimization and differential evolution [M]. Journal of computational information system, 2012.

    [18] Singh R M, Datta B. Identification of groundwater pollution sources using GA-based linked simulation optimization model [J]. Journal of Hydrologic Engineering, 2006,11(2):1216-1227.

    [19] 高岳林,劉軍民.差分進(jìn)化算法的參數(shù)研究 [J]. 黑龍江大學(xué)自然科學(xué)學(xué)報(bào), 2009,26(1):81-85.

    [20] 鄭 彤.環(huán)境系統(tǒng)數(shù)學(xué)模型[M]. 北京:化學(xué)工業(yè)出版社, 2003.

    Research of the abrupt waters pollution source based on optimization algorithm of PSO-DE.

    CAO Hong-gui, YUN Wei-guo*

    (School of Information and Control Engineering, Xi'an University of Architecture and Technology, Xi’an 710055, China)., 2017,37(10):3807~3812

    The optimization algorithm of PSO-DE cooperated with mobile monitoring platform was studied to solve the inverse problem of pollution, which included inversion of the position of the single point and multiple-point stationary sources. The inverse problem of pollution source was transformed into nonlinear optimization problem. The pollutant concentration of waters were detected and recorded bymobile platforms; the coordinate of mobile platform was denoted bybest, and they were corresponded one by one, there would be Nbestaltogether. The pollutant concentration of waters which attained by the mobile platform would be compared with each other, and the coordinate of maximum pollutant concentration would be chosen and marked asbest. First, thebestwould be the initial population for the PSO optimization. Second, the population would be used for DE optimization. Finally, thebestwould be chosen from the high concentration of both until the highest point of pollutant concentration was obtained, which was the initial point of pollutant. The calculation results of examples showed that the algorithm could attained a high precision inversion results for pollutant source traceability problem of two-dimensional waters.

    PSO-DE;pollutant source traceability;mobile monitoring platform;two-dimensional waters

    X52

    A

    1000-6923(2017)10-3807-06

    曹宏桂(1991-),男,陜西榆林人,西安建筑科技大學(xué)碩士研究生,主要研究領(lǐng)域?yàn)榍度胧介_發(fā),智能控制.發(fā)表論文2篇.

    2017-03-14

    住房和城鄉(xiāng)建設(shè)部科學(xué)項(xiàng)目計(jì)劃(2016-R2-045)

    * 責(zé)任作者, 教授, yun_weiguo@aliyun.com

    猜你喜歡
    優(yōu)化
    超限高層建筑結(jié)構(gòu)設(shè)計(jì)與優(yōu)化思考
    PEMFC流道的多目標(biāo)優(yōu)化
    能源工程(2022年1期)2022-03-29 01:06:28
    民用建筑防煙排煙設(shè)計(jì)優(yōu)化探討
    關(guān)于優(yōu)化消防安全告知承諾的一些思考
    一道優(yōu)化題的幾何解法
    由“形”啟“數(shù)”優(yōu)化運(yùn)算——以2021年解析幾何高考題為例
    圍繞“地、業(yè)、人”優(yōu)化產(chǎn)業(yè)扶貧
    事業(yè)單位中固定資產(chǎn)會(huì)計(jì)處理的優(yōu)化
    4K HDR性能大幅度優(yōu)化 JVC DLA-X8 18 BC
    幾種常見的負(fù)載均衡算法的優(yōu)化
    電子制作(2017年20期)2017-04-26 06:57:45
    老司机午夜十八禁免费视频| 午夜精品久久久久久毛片777| 精品国产三级普通话版| 18+在线观看网站| 国产麻豆成人av免费视频| 亚洲无线在线观看| 亚洲精品亚洲一区二区| 波多野结衣巨乳人妻| 亚洲成av人片免费观看| 国产一区二区激情短视频| 无遮挡黄片免费观看| 亚洲国产精品合色在线| 欧美高清成人免费视频www| 亚洲国产色片| 亚洲av电影在线进入| 午夜福利18| 两性午夜刺激爽爽歪歪视频在线观看| 丰满人妻熟妇乱又伦精品不卡| 久久九九热精品免费| 国产欧美日韩精品亚洲av| 免费在线观看影片大全网站| 国产亚洲精品av在线| 99久国产av精品| 一个人观看的视频www高清免费观看| 一进一出好大好爽视频| 精品一区二区三区人妻视频| 免费在线观看日本一区| 日韩有码中文字幕| 女生性感内裤真人,穿戴方法视频| 国产乱人视频| 国产探花极品一区二区| 精品久久久久久成人av| 热99re8久久精品国产| 俺也久久电影网| 国产高清视频在线播放一区| 看片在线看免费视频| 午夜福利成人在线免费观看| 欧美成狂野欧美在线观看| 身体一侧抽搐| 欧美xxxx性猛交bbbb| 深爱激情五月婷婷| 亚洲片人在线观看| 免费在线观看影片大全网站| 国产亚洲精品av在线| 天堂网av新在线| 美女高潮的动态| 性色avwww在线观看| 国产亚洲欧美在线一区二区| 国产午夜精品久久久久久一区二区三区 | 亚洲自偷自拍三级| 亚洲自拍偷在线| 午夜福利在线观看免费完整高清在 | 一区二区三区免费毛片| 亚洲性夜色夜夜综合| 久久中文看片网| 久久中文看片网| 欧美一区二区国产精品久久精品| 亚洲国产欧洲综合997久久,| 天堂动漫精品| 国产一区二区在线av高清观看| 欧美一区二区国产精品久久精品| 久久久久久久久大av| 午夜影院日韩av| 波多野结衣高清无吗| 日韩人妻高清精品专区| 2021天堂中文幕一二区在线观| 国产午夜福利久久久久久| 久久精品91蜜桃| 伦理电影大哥的女人| 国产精品不卡视频一区二区 | 91午夜精品亚洲一区二区三区 | 淫妇啪啪啪对白视频| 亚洲国产精品成人综合色| 一级a爱片免费观看的视频| 国产黄片美女视频| 国产白丝娇喘喷水9色精品| 久久久久免费精品人妻一区二区| 国产在视频线在精品| 真实男女啪啪啪动态图| 国产成人欧美在线观看| 老司机午夜十八禁免费视频| 一级毛片久久久久久久久女| 亚洲av.av天堂| 国产主播在线观看一区二区| 美女大奶头视频| 神马国产精品三级电影在线观看| 午夜福利在线观看免费完整高清在 | 精华霜和精华液先用哪个| 欧美+亚洲+日韩+国产| 1000部很黄的大片| 免费无遮挡裸体视频| 在线观看舔阴道视频| 国产色爽女视频免费观看| 老司机午夜福利在线观看视频| 久久久久国内视频| 色播亚洲综合网| 久久久久久国产a免费观看| avwww免费| 国产乱人视频| 久久精品91蜜桃| 亚洲va日本ⅴa欧美va伊人久久| 人妻丰满熟妇av一区二区三区| 国产日本99.免费观看| 欧美成人一区二区免费高清观看| 90打野战视频偷拍视频| 丰满人妻一区二区三区视频av| 搡老岳熟女国产| 波多野结衣高清无吗| 国内少妇人妻偷人精品xxx网站| 永久网站在线| 久久香蕉精品热| 亚洲最大成人中文| 亚洲中文日韩欧美视频| 国内精品久久久久精免费| 久久久久精品国产欧美久久久| 亚洲男人的天堂狠狠| 成人鲁丝片一二三区免费| 三级毛片av免费| 国产美女午夜福利| 欧美色欧美亚洲另类二区| 中文字幕精品亚洲无线码一区| 非洲黑人性xxxx精品又粗又长| 国产真实乱freesex| 久久午夜亚洲精品久久| 男人和女人高潮做爰伦理| 成人无遮挡网站| 一级毛片久久久久久久久女| 极品教师在线免费播放| 毛片一级片免费看久久久久 | 国产成+人综合+亚洲专区| 老熟妇乱子伦视频在线观看| 九九在线视频观看精品| 亚洲国产高清在线一区二区三| 亚洲欧美日韩东京热| netflix在线观看网站| 色综合婷婷激情| 亚洲成av人片在线播放无| 国产熟女xx| 免费av观看视频| 欧美激情国产日韩精品一区| 国产精品国产高清国产av| 丁香六月欧美| netflix在线观看网站| 一本综合久久免费| 国产精品98久久久久久宅男小说| 在线a可以看的网站| av天堂中文字幕网| 中文亚洲av片在线观看爽| 激情在线观看视频在线高清| 久久精品91蜜桃| 国产高潮美女av| 1000部很黄的大片| 97碰自拍视频| 亚洲真实伦在线观看| 午夜免费成人在线视频| 国产乱人视频| 最好的美女福利视频网| 日本与韩国留学比较| 亚洲美女黄片视频| 一级a爱片免费观看的视频| 亚洲美女视频黄频| 久久99热这里只有精品18| 国产91精品成人一区二区三区| 一区福利在线观看| 99国产综合亚洲精品| 欧美成人性av电影在线观看| 搡老妇女老女人老熟妇| 国产精品精品国产色婷婷| 黄色视频,在线免费观看| 欧美精品国产亚洲| 国产不卡一卡二| 国产一级毛片七仙女欲春2| 高潮久久久久久久久久久不卡| 男人狂女人下面高潮的视频| 亚洲人成电影免费在线| 国模一区二区三区四区视频| 欧美在线一区亚洲| 久久九九热精品免费| 亚洲国产色片| 成人鲁丝片一二三区免费| 欧美国产日韩亚洲一区| 夜夜躁狠狠躁天天躁| 黄色丝袜av网址大全| 欧美最黄视频在线播放免费| 亚洲天堂国产精品一区在线| 别揉我奶头 嗯啊视频| 久久久久精品国产欧美久久久| 一卡2卡三卡四卡精品乱码亚洲| 中出人妻视频一区二区| 午夜激情欧美在线| x7x7x7水蜜桃| 久久99热6这里只有精品| 亚洲人成伊人成综合网2020| 日本撒尿小便嘘嘘汇集6| 亚洲自偷自拍三级| 国产高清视频在线观看网站| 天堂√8在线中文| 亚洲精品成人久久久久久| 精品国产亚洲在线| .国产精品久久| 国产黄片美女视频| 天天一区二区日本电影三级| 91九色精品人成在线观看| 国产aⅴ精品一区二区三区波| 日韩 亚洲 欧美在线| 夜夜爽天天搞| 99国产精品一区二区蜜桃av| 免费av毛片视频| 免费高清视频大片| 亚洲成人精品中文字幕电影| 亚洲av成人不卡在线观看播放网| h日本视频在线播放| 搡老熟女国产l中国老女人| 99热精品在线国产| 天堂网av新在线| 男女那种视频在线观看| 欧美国产日韩亚洲一区| 精品午夜福利视频在线观看一区| 成年免费大片在线观看| 成年免费大片在线观看| 日韩有码中文字幕| 久久九九热精品免费| 天天躁日日操中文字幕| 久久6这里有精品| 极品教师在线免费播放| 亚洲av电影在线进入| 一个人免费在线观看的高清视频| 听说在线观看完整版免费高清| 国产精品日韩av在线免费观看| 两个人视频免费观看高清| av视频在线观看入口| 欧美乱妇无乱码| 亚洲五月天丁香| 亚洲七黄色美女视频| 久久久国产成人免费| 成年女人毛片免费观看观看9| 黄色丝袜av网址大全| 久久久久免费精品人妻一区二区| 欧美绝顶高潮抽搐喷水| 亚洲最大成人中文| 一个人看视频在线观看www免费| 中文字幕高清在线视频| 国产精品98久久久久久宅男小说| 美女cb高潮喷水在线观看| 他把我摸到了高潮在线观看| 欧美午夜高清在线| 国内揄拍国产精品人妻在线| 久久精品国产清高在天天线| 深夜精品福利| 久久久精品欧美日韩精品| 亚洲七黄色美女视频| 国产精品久久久久久精品电影| 丁香六月欧美| 一本精品99久久精品77| 在线国产一区二区在线| 搞女人的毛片| 999久久久精品免费观看国产| 日本黄色视频三级网站网址| 免费搜索国产男女视频| 性色av乱码一区二区三区2| 欧美三级亚洲精品| .国产精品久久| 午夜视频国产福利| 国产av麻豆久久久久久久| 99久久成人亚洲精品观看| 久久久国产成人精品二区| 欧美在线一区亚洲| 少妇裸体淫交视频免费看高清| 亚洲综合色惰| 日韩欧美精品v在线| www.色视频.com| 国产高清视频在线观看网站| 黄色视频,在线免费观看| 女人十人毛片免费观看3o分钟| 成熟少妇高潮喷水视频| 欧美中文日本在线观看视频| 精品久久久久久久久久久久久| 久久人妻av系列| 亚洲av不卡在线观看| 高潮久久久久久久久久久不卡| 九色国产91popny在线| 亚洲三级黄色毛片| 国产精品亚洲av一区麻豆| 亚洲一区二区三区不卡视频| 国产亚洲精品综合一区在线观看| 在线看三级毛片| 色视频www国产| aaaaa片日本免费| 俺也久久电影网| 观看美女的网站| 丝袜美腿在线中文| 久久久久九九精品影院| 一本久久中文字幕| 99精品久久久久人妻精品| 精品日产1卡2卡| 国产精品亚洲av一区麻豆| 亚洲男人的天堂狠狠| 亚洲中文日韩欧美视频| 特级一级黄色大片| 看黄色毛片网站| 国产精品亚洲av一区麻豆| 十八禁国产超污无遮挡网站| 国产av麻豆久久久久久久| 淫妇啪啪啪对白视频| 亚洲一区二区三区不卡视频| 伦理电影大哥的女人| 亚洲第一区二区三区不卡| 成人av一区二区三区在线看| 神马国产精品三级电影在线观看| 午夜福利欧美成人| 黄片小视频在线播放| 亚洲av日韩精品久久久久久密| 国产伦人伦偷精品视频| 午夜福利视频1000在线观看| 亚洲,欧美精品.| 搡老岳熟女国产| 亚洲专区中文字幕在线| 精品久久久久久久久亚洲 | 国产不卡一卡二| 成人三级黄色视频| 国产蜜桃级精品一区二区三区| 色噜噜av男人的天堂激情| 欧美黄色片欧美黄色片| 色吧在线观看| 欧美性猛交╳xxx乱大交人| 我要搜黄色片| 国产黄色小视频在线观看| 亚洲av一区综合| 在线观看美女被高潮喷水网站 | 一进一出抽搐gif免费好疼| 精品国产亚洲在线| 波野结衣二区三区在线| 91久久精品电影网| 国产亚洲精品久久久久久毛片| 久久热精品热| 婷婷六月久久综合丁香| 三级男女做爰猛烈吃奶摸视频| 欧美zozozo另类| 免费观看人在逋| 久久久久国内视频| 成年女人毛片免费观看观看9| 国产av一区在线观看免费| 精品久久久久久久久久久久久| 亚洲国产精品合色在线| 动漫黄色视频在线观看| 国产精品国产高清国产av| 中文字幕久久专区| 欧美一级a爱片免费观看看| 亚洲无线观看免费| 成人鲁丝片一二三区免费| 国产精品永久免费网站| 国产亚洲精品综合一区在线观看| 国产成人欧美在线观看| 精品久久久久久成人av| 亚洲avbb在线观看| av专区在线播放| 内射极品少妇av片p| 99久久久亚洲精品蜜臀av| 国产免费一级a男人的天堂| 亚洲 国产 在线| 精品日产1卡2卡| 国产单亲对白刺激| 亚洲avbb在线观看| 乱码一卡2卡4卡精品| 国产午夜精品论理片| 国产毛片a区久久久久| 日本黄大片高清| av在线蜜桃| 亚洲真实伦在线观看| 欧美bdsm另类| 亚洲自偷自拍三级| 国产国拍精品亚洲av在线观看| 好男人电影高清在线观看| 免费看a级黄色片| 国产精品久久久久久精品电影| 变态另类丝袜制服| 久久久精品欧美日韩精品| a在线观看视频网站| 97超级碰碰碰精品色视频在线观看| 日韩高清综合在线| 美女高潮的动态| 国产探花在线观看一区二区| 国产白丝娇喘喷水9色精品| 可以在线观看毛片的网站| 欧美+日韩+精品| 2021天堂中文幕一二区在线观| 日本精品一区二区三区蜜桃| 成年人黄色毛片网站| 校园春色视频在线观看| 成人美女网站在线观看视频| 小蜜桃在线观看免费完整版高清| 成人永久免费在线观看视频| av中文乱码字幕在线| 香蕉av资源在线| 88av欧美| 丁香欧美五月| 在线免费观看的www视频| 亚洲熟妇熟女久久| 午夜老司机福利剧场| а√天堂www在线а√下载| 国产人妻一区二区三区在| 又黄又爽又刺激的免费视频.| 麻豆成人av在线观看| 深夜a级毛片| 琪琪午夜伦伦电影理论片6080| 欧美性猛交╳xxx乱大交人| 好男人在线观看高清免费视频| 午夜激情欧美在线| 日韩大尺度精品在线看网址| 日本熟妇午夜| 亚洲国产高清在线一区二区三| 国产免费一级a男人的天堂| 大型黄色视频在线免费观看| 国产单亲对白刺激| 啦啦啦观看免费观看视频高清| 看免费av毛片| 日韩中文字幕欧美一区二区| 国产精品亚洲美女久久久| 久久热精品热| 天堂动漫精品| 舔av片在线| bbb黄色大片| 亚洲av成人精品一区久久| 婷婷精品国产亚洲av| 国产不卡一卡二| 免费观看人在逋| 桃色一区二区三区在线观看| 欧美日本视频| 99久久99久久久精品蜜桃| 久久久久久九九精品二区国产| 天堂√8在线中文| 好男人在线观看高清免费视频| 一级av片app| 天天一区二区日本电影三级| 欧美中文日本在线观看视频| 99久久精品热视频| 久久久成人免费电影| 国产精品免费一区二区三区在线| 一进一出好大好爽视频| 岛国在线免费视频观看| 日本 欧美在线| av在线老鸭窝| 淫妇啪啪啪对白视频| 黄色一级大片看看| 国产精品综合久久久久久久免费| 嫩草影视91久久| 国产精品一区二区三区四区久久| 桃红色精品国产亚洲av| 久久精品影院6| 国产淫片久久久久久久久 | 成人国产一区最新在线观看| 国产精品久久久久久精品电影| 美女xxoo啪啪120秒动态图 | 成年人黄色毛片网站| 久久久久精品国产欧美久久久| 直男gayav资源| 又紧又爽又黄一区二区| 中文字幕人妻熟人妻熟丝袜美| 国产亚洲精品久久久com| 国产成人aa在线观看| 久久久久久久亚洲中文字幕 | 久久精品国产亚洲av天美| 一个人看的www免费观看视频| 成人国产一区最新在线观看| 亚洲精品粉嫩美女一区| 免费大片18禁| 久9热在线精品视频| 成人亚洲精品av一区二区| 在线十欧美十亚洲十日本专区| 久久久久亚洲av毛片大全| 亚洲av电影在线进入| x7x7x7水蜜桃| 日本三级黄在线观看| 天堂影院成人在线观看| 日日摸夜夜添夜夜添av毛片 | 亚洲狠狠婷婷综合久久图片| xxxwww97欧美| 2021天堂中文幕一二区在线观| 此物有八面人人有两片| 久久这里只有精品中国| 97热精品久久久久久| 级片在线观看| 中文亚洲av片在线观看爽| 欧美黄色片欧美黄色片| 亚洲欧美日韩高清专用| 精品久久久久久久人妻蜜臀av| 国产精品综合久久久久久久免费| 国产私拍福利视频在线观看| 美女高潮的动态| 91九色精品人成在线观看| 婷婷亚洲欧美| 精品久久久久久,| 色播亚洲综合网| 老司机福利观看| 国产视频一区二区在线看| 黄片小视频在线播放| 色精品久久人妻99蜜桃| 国产 一区 欧美 日韩| 久久精品91蜜桃| 亚洲激情在线av| 一级作爱视频免费观看| 看黄色毛片网站| 国产三级在线视频| 一区二区三区四区激情视频 | 又黄又爽又免费观看的视频| 午夜福利在线在线| 午夜亚洲福利在线播放| 非洲黑人性xxxx精品又粗又长| 99热这里只有是精品在线观看 | 老熟妇乱子伦视频在线观看| 国产一区二区在线观看日韩| 毛片一级片免费看久久久久 | 听说在线观看完整版免费高清| 国产精品爽爽va在线观看网站| 在线观看午夜福利视频| 国产精品一区二区免费欧美| 免费av不卡在线播放| 日本五十路高清| 免费高清视频大片| 可以在线观看的亚洲视频| 成人永久免费在线观看视频| 人人妻人人澡欧美一区二区| 一夜夜www| 日本一本二区三区精品| 大型黄色视频在线免费观看| 美女免费视频网站| 精品久久久久久久末码| 精品一区二区三区av网在线观看| 网址你懂的国产日韩在线| 久久精品影院6| 一进一出抽搐动态| 天堂av国产一区二区熟女人妻| 日本五十路高清| 91久久精品国产一区二区成人| 一边摸一边抽搐一进一小说| 欧美一级a爱片免费观看看| 国产精品久久久久久亚洲av鲁大| 在线免费观看的www视频| 成人av在线播放网站| 国产精品99久久久久久久久| 国产免费一级a男人的天堂| 欧美精品国产亚洲| 亚洲欧美激情综合另类| 成人国产综合亚洲| 在线观看舔阴道视频| 身体一侧抽搐| 色综合亚洲欧美另类图片| 久久久国产成人精品二区| 亚洲人成网站在线播放欧美日韩| 久久久精品大字幕| 欧美潮喷喷水| 久久精品人妻少妇| 深夜精品福利| 欧美色欧美亚洲另类二区| 欧美极品一区二区三区四区| 国产精品一区二区三区四区免费观看 | 成人特级av手机在线观看| 亚洲精品456在线播放app | 国产精品自产拍在线观看55亚洲| 级片在线观看| 国产三级在线视频| 日本精品一区二区三区蜜桃| 一个人免费在线观看的高清视频| 亚洲七黄色美女视频| 免费黄网站久久成人精品 | 窝窝影院91人妻| 五月伊人婷婷丁香| 免费观看的影片在线观看| 精品午夜福利视频在线观看一区| 波野结衣二区三区在线| 男人舔奶头视频| 色精品久久人妻99蜜桃| 国产aⅴ精品一区二区三区波| 黄色丝袜av网址大全| 一级a爱片免费观看的视频| 亚洲欧美清纯卡通| 午夜免费男女啪啪视频观看 | av在线观看视频网站免费| 久久亚洲真实| 91九色精品人成在线观看| 久久久久久久久久黄片| 国产亚洲欧美98| 在线观看66精品国产| 日韩欧美免费精品| 成年人黄色毛片网站| 内射极品少妇av片p| 十八禁网站免费在线| 亚洲性夜色夜夜综合| 最近最新中文字幕大全电影3| 国产乱人视频| 免费人成视频x8x8入口观看| 精品久久久久久久人妻蜜臀av| 亚洲,欧美,日韩| 欧美日本亚洲视频在线播放| 男插女下体视频免费在线播放| 亚洲va日本ⅴa欧美va伊人久久| 亚洲熟妇熟女久久| 一卡2卡三卡四卡精品乱码亚洲| 悠悠久久av| 国产精品三级大全| 欧美另类亚洲清纯唯美| 人妻夜夜爽99麻豆av| 嫩草影院新地址| 成人一区二区视频在线观看| 亚洲国产精品999在线| 精品人妻熟女av久视频| 日韩精品中文字幕看吧| 免费看光身美女| 女生性感内裤真人,穿戴方法视频| 日本在线视频免费播放| 村上凉子中文字幕在线| 男人和女人高潮做爰伦理| 757午夜福利合集在线观看| 国产三级黄色录像| 18禁裸乳无遮挡免费网站照片|