• <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)國(guó)
    中國(guó)環(huán)境科學(xué) 2017年10期
    關(guān)鍵詞:優(yōu)化

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

    ?

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

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

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

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

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

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

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

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

    模型為無(wú)限邊際水體假設(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í)間而改變,只和與污染源的相對(duì)坐標(biāo)有關(guān),與時(shí)間無(wú)關(guān).

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

    1.2 二維無(wú)限水域多點(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í)間無(wú)關(guān).式2可以用來(lái)模擬大型二維水域同時(shí)存在多個(gè)同類點(diǎn)源污染物穩(wěn)定排污情景下的水域各個(gè)斷面的污染物濃度.

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

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

    為了改善粒子群進(jìn)化算法和微分進(jìn)化算法獨(dú)自求解優(yōu)化問(wèn)題的不足.本文引用文獻(xiàn)中的粒子群算法和微分進(jìn)化算法混合的方法.該算法是對(duì)種群先后做PSO和DE優(yōu)化的策略,對(duì)于一個(gè)種群,先通過(guò)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ì)使粒子過(guò)于發(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):通過(guò)下列公式修正速度與位置:

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

    3 數(shù)值算例

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

    3.1.1 算例來(lái)源 假設(shè)某個(gè)大型的湖泊,在位置A排放2000kg惰性污染物,該湖泊的平均水深為l0m,平均水流速度0.5m/s,橫向彌散系數(shù)50m2/ min[20].該水域相關(guān)斷面的污染物濃度可以通過(guò)公式(1)來(lái)計(jì)算.令位置A為坐標(biāo)原點(diǎn)(0,0),選取水域縱向100m(步長(zhǎng)為0.5)橫向200m(步長(zhǎng)為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 算例來(lái)源 假設(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é)果分析 對(duì)于多點(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)值

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

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

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

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

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

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

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

    4 結(jié)論

    本研究考慮了污染物水質(zhì)特征,給出二維水域下適用于移動(dòng)平臺(tái)的PSO-DE優(yōu)化溯源算法.對(duì)PSO-DE算法中的慣性權(quán)重因子和時(shí)變學(xué)習(xí)因子進(jìn)行改進(jìn),通過(guò)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í)可能遇到的問(wèn)題,具有一定的應(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)識(shí)別 [J]. 中國(guó)環(huán)境科學(xué), 2011,31(11):1802-1807.

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

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

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

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

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

    [10] 牟行洋.基于微分進(jìn)化算法的污染物源項(xiàng)識(shí)別反問(wèn)題研究術(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)化方法在地下水污染源識(shí)別中的應(yīng)用 [J]. 中國(guó)環(huán)境科學(xué), 2016,36(7): 2067-2072.

    [16] 傅國(guó)偉.河流水質(zhì)數(shù)學(xué)模型及其模擬計(jì)算[M]. 北京:中國(guó)環(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)榍度胧介_(kāi)發(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
    幾種常見(jiàn)的負(fù)載均衡算法的優(yōu)化
    電子制作(2017年20期)2017-04-26 06:57:45
    天堂影院成人在线观看| 日本精品一区二区三区蜜桃| 91午夜精品亚洲一区二区三区 | 亚洲最大成人中文| 中文在线观看免费www的网站| 又黄又爽又免费观看的视频| 欧美人与善性xxx| 真实男女啪啪啪动态图| 啦啦啦啦在线视频资源| 亚洲va在线va天堂va国产| 99久国产av精品| 久久久久性生活片| 日本 欧美在线| 午夜亚洲福利在线播放| 亚洲18禁久久av| 香蕉av资源在线| 精品人妻视频免费看| 国产精品一区二区性色av| 亚洲av五月六月丁香网| 精品国内亚洲2022精品成人| 国产精品久久久久久久久免| 性插视频无遮挡在线免费观看| 97碰自拍视频| 国产精品亚洲美女久久久| 天堂av国产一区二区熟女人妻| 国国产精品蜜臀av免费| 国产老妇女一区| 亚洲经典国产精华液单| 国产免费男女视频| 欧美一级a爱片免费观看看| 欧美另类亚洲清纯唯美| 禁无遮挡网站| 久久午夜福利片| 天天躁日日操中文字幕| 亚洲aⅴ乱码一区二区在线播放| 嫩草影院精品99| 舔av片在线| 淫秽高清视频在线观看| 听说在线观看完整版免费高清| 亚洲最大成人av| 乱人视频在线观看| 久久久精品欧美日韩精品| 琪琪午夜伦伦电影理论片6080| 亚洲av成人av| 最新中文字幕久久久久| 国产免费av片在线观看野外av| 国产 一区精品| 一a级毛片在线观看| av在线观看视频网站免费| 亚洲精品亚洲一区二区| 精品久久久久久久久av| 精品人妻1区二区| 国产精品久久久久久亚洲av鲁大| 99国产精品一区二区蜜桃av| 国产成人福利小说| 九九热线精品视视频播放| 国产精品自产拍在线观看55亚洲| 麻豆国产97在线/欧美| 深爱激情五月婷婷| 国产主播在线观看一区二区| 亚洲成人精品中文字幕电影| 午夜影院日韩av| 男人舔奶头视频| 97碰自拍视频| 色在线成人网| 国产一区二区三区av在线 | 久久国内精品自在自线图片| 色播亚洲综合网| 亚洲黑人精品在线| 欧美zozozo另类| h日本视频在线播放| 亚洲欧美精品综合久久99| 别揉我奶头 嗯啊视频| 在线播放国产精品三级| 午夜久久久久精精品| 欧美一级a爱片免费观看看| 老师上课跳d突然被开到最大视频| 亚洲图色成人| 精品久久久久久成人av| 久久久久精品国产欧美久久久| 亚洲第一电影网av| 国产av一区在线观看免费| 国产综合懂色| av在线观看视频网站免费| 亚洲最大成人中文| 久久久久免费精品人妻一区二区| 欧美人与善性xxx| 欧美+日韩+精品| 亚洲av二区三区四区| 欧美日韩综合久久久久久 | 精品久久久久久,| 午夜福利欧美成人| 他把我摸到了高潮在线观看| 色综合亚洲欧美另类图片| 亚洲18禁久久av| 亚洲精品粉嫩美女一区| 国产黄色小视频在线观看| www.www免费av| 嫩草影院入口| 午夜福利在线观看吧| 免费观看的影片在线观看| 少妇人妻精品综合一区二区 | 亚洲国产日韩欧美精品在线观看| 老司机深夜福利视频在线观看| 黄色配什么色好看| 国产精品1区2区在线观看.| 在线观看免费视频日本深夜| 12—13女人毛片做爰片一| 亚洲精品亚洲一区二区| 国产精品无大码| 国产主播在线观看一区二区| 精品乱码久久久久久99久播| 久久人人爽人人爽人人片va| 乱人视频在线观看| 成人午夜高清在线视频| 一个人观看的视频www高清免费观看| 一本一本综合久久| 99热这里只有是精品在线观看| 特大巨黑吊av在线直播| 日韩欧美在线二视频| 人妻丰满熟妇av一区二区三区| 中文字幕av在线有码专区| 看十八女毛片水多多多| 身体一侧抽搐| 国产熟女欧美一区二区| 老熟妇乱子伦视频在线观看| 又紧又爽又黄一区二区| 成人av一区二区三区在线看| 亚洲美女黄片视频| 最好的美女福利视频网| 美女免费视频网站| 男女做爰动态图高潮gif福利片| 亚洲一级一片aⅴ在线观看| 成人特级av手机在线观看| 天美传媒精品一区二区| 日本一二三区视频观看| 日韩亚洲欧美综合| 精品人妻熟女av久视频| 亚洲欧美精品综合久久99| 日韩欧美国产在线观看| 日日摸夜夜添夜夜添小说| 免费看光身美女| 麻豆成人av在线观看| 成人av一区二区三区在线看| 久久精品久久久久久噜噜老黄 | 亚洲自拍偷在线| 国产av在哪里看| 国产蜜桃级精品一区二区三区| 国产精品免费一区二区三区在线| 色哟哟·www| av在线天堂中文字幕| 毛片女人毛片| 男插女下体视频免费在线播放| 麻豆国产av国片精品| 美女 人体艺术 gogo| 91麻豆av在线| 亚洲av成人精品一区久久| 亚洲性久久影院| av天堂中文字幕网| 99久久久亚洲精品蜜臀av| 国产 一区精品| 国产成年人精品一区二区| 亚洲无线在线观看| 中亚洲国语对白在线视频| 久久国内精品自在自线图片| 亚洲四区av| 97超级碰碰碰精品色视频在线观看| 女的被弄到高潮叫床怎么办 | 亚洲av不卡在线观看| 中文资源天堂在线| 成人综合一区亚洲| 国产色爽女视频免费观看| 美女高潮的动态| 日本黄色片子视频| 亚洲在线观看片| 国产毛片a区久久久久| 精品一区二区三区人妻视频| 亚洲性久久影院| 亚洲色图av天堂| 日韩欧美 国产精品| 18禁黄网站禁片免费观看直播| 97超视频在线观看视频| 欧美一区二区国产精品久久精品| 国产精品久久久久久av不卡| bbb黄色大片| 久久久久精品国产欧美久久久| 免费一级毛片在线播放高清视频| 亚洲专区国产一区二区| 无遮挡黄片免费观看| 久久久久久久久久黄片| 两个人视频免费观看高清| 精品人妻熟女av久视频| 色尼玛亚洲综合影院| 欧美精品啪啪一区二区三区| .国产精品久久| 久久热精品热| 男女做爰动态图高潮gif福利片| 久久精品影院6| 日本 欧美在线| 亚洲一区高清亚洲精品| 大又大粗又爽又黄少妇毛片口| 在线观看免费视频日本深夜| 99国产极品粉嫩在线观看| 国产亚洲91精品色在线| 性插视频无遮挡在线免费观看| 成年免费大片在线观看| 男女那种视频在线观看| 欧美不卡视频在线免费观看| 亚洲 国产 在线| 精品一区二区免费观看| 国产精品爽爽va在线观看网站| 日本一二三区视频观看| 一卡2卡三卡四卡精品乱码亚洲| 麻豆一二三区av精品| 精品不卡国产一区二区三区| 自拍偷自拍亚洲精品老妇| 亚洲av中文字字幕乱码综合| 婷婷亚洲欧美| 悠悠久久av| 国产精品一及| 一区二区三区四区激情视频 | 看十八女毛片水多多多| 亚洲自拍偷在线| 床上黄色一级片| 天天躁日日操中文字幕| 亚洲成a人片在线一区二区| 亚洲欧美日韩高清在线视频| 国产精品1区2区在线观看.| 成人三级黄色视频| 乱人视频在线观看| 非洲黑人性xxxx精品又粗又长| 十八禁国产超污无遮挡网站| 欧美不卡视频在线免费观看| 五月伊人婷婷丁香| 极品教师在线免费播放| 在线观看免费视频日本深夜| 一进一出抽搐gif免费好疼| 在线免费观看的www视频| 91在线精品国自产拍蜜月| 韩国av在线不卡| 真人做人爱边吃奶动态| 国产精品亚洲一级av第二区| 国产主播在线观看一区二区| 欧美成人性av电影在线观看| 国产精品久久久久久久久免| 亚洲欧美日韩高清在线视频| 熟女电影av网| 热99在线观看视频| 99国产极品粉嫩在线观看| 国内精品久久久久精免费| 一个人看的www免费观看视频| 国产一区二区在线观看日韩| 国产高清有码在线观看视频| 亚洲avbb在线观看| 看十八女毛片水多多多| 欧美激情在线99| 99久国产av精品| 欧美日韩瑟瑟在线播放| 婷婷色综合大香蕉| 观看美女的网站| 久久99热6这里只有精品| 中文字幕高清在线视频| 午夜激情福利司机影院| 精品久久久噜噜| 日本成人三级电影网站| 99热网站在线观看| 国模一区二区三区四区视频| 老司机午夜福利在线观看视频| 亚洲无线观看免费| 日韩一本色道免费dvd| 丰满乱子伦码专区| 久99久视频精品免费| 免费人成在线观看视频色| 看十八女毛片水多多多| 精品一区二区三区人妻视频| 国产精品,欧美在线| 最后的刺客免费高清国语| 能在线免费观看的黄片| 乱系列少妇在线播放| 老女人水多毛片| 全区人妻精品视频| 色吧在线观看| 亚洲av电影不卡..在线观看| 亚洲天堂国产精品一区在线| 欧美绝顶高潮抽搐喷水| 2021天堂中文幕一二区在线观| 亚洲性夜色夜夜综合| 成人二区视频| 国产精品一区二区三区四区免费观看 | 国产黄色小视频在线观看| 如何舔出高潮| 99久久九九国产精品国产免费| 久久这里只有精品中国| 成人无遮挡网站| 天堂影院成人在线观看| 男人舔奶头视频| 久久精品综合一区二区三区| 麻豆成人午夜福利视频| 亚洲三级黄色毛片| 久久热精品热| 露出奶头的视频| 国产v大片淫在线免费观看| 亚洲国产精品合色在线| 国产不卡一卡二| 久久天躁狠狠躁夜夜2o2o| 中文字幕人妻熟人妻熟丝袜美| 一区二区三区激情视频| 色哟哟·www| 无人区码免费观看不卡| 国产伦在线观看视频一区| 国产精品久久久久久亚洲av鲁大| 18禁黄网站禁片免费观看直播| 观看美女的网站| 国产精品久久电影中文字幕| 国产乱人伦免费视频| 国产精品久久久久久精品电影| 欧美黑人巨大hd| 在线播放国产精品三级| 国产精品98久久久久久宅男小说| 久久久久久久久大av| 久久久久国内视频| 亚洲人与动物交配视频| 99久久成人亚洲精品观看| 色哟哟哟哟哟哟| 免费高清视频大片| 久久久久久久精品吃奶| 午夜激情欧美在线| 国产精品野战在线观看| 欧美精品啪啪一区二区三区| 色综合亚洲欧美另类图片| 日本一二三区视频观看| 啦啦啦韩国在线观看视频| 国产乱人视频| 偷拍熟女少妇极品色| 老司机福利观看| 国产视频一区二区在线看| 免费不卡的大黄色大毛片视频在线观看 | 久久久久精品国产欧美久久久| 国产熟女欧美一区二区| 免费不卡的大黄色大毛片视频在线观看 | 午夜日韩欧美国产| 日韩中字成人| 久9热在线精品视频| 国产激情偷乱视频一区二区| 高清毛片免费观看视频网站| 午夜精品一区二区三区免费看| 又粗又爽又猛毛片免费看| 搡老妇女老女人老熟妇| 色尼玛亚洲综合影院| 国产女主播在线喷水免费视频网站 | 国产高清视频在线观看网站| 春色校园在线视频观看| 国产又黄又爽又无遮挡在线| 国产乱人视频| 亚洲欧美日韩卡通动漫| 真实男女啪啪啪动态图| 1000部很黄的大片| 在线观看美女被高潮喷水网站| 欧洲精品卡2卡3卡4卡5卡区| 三级男女做爰猛烈吃奶摸视频| 国产高清视频在线观看网站| 午夜激情福利司机影院| 日韩精品中文字幕看吧| 精品国产三级普通话版| 18禁裸乳无遮挡免费网站照片| 99热精品在线国产| 国产一区二区三区av在线 | 网址你懂的国产日韩在线| 日本 av在线| 日韩亚洲欧美综合| 国产成人一区二区在线| 国内揄拍国产精品人妻在线| 美女cb高潮喷水在线观看| 国产亚洲精品综合一区在线观看| 亚洲乱码一区二区免费版| 嫩草影院入口| 免费一级毛片在线播放高清视频| 精品久久久久久久久亚洲 | 在线a可以看的网站| 国产精品日韩av在线免费观看| 91精品国产九色| 成年女人看的毛片在线观看| 免费看美女性在线毛片视频| 国产av一区在线观看免费| 成人av一区二区三区在线看| 免费看美女性在线毛片视频| 久久国产精品人妻蜜桃| av天堂在线播放| 精品久久久久久久末码| av黄色大香蕉| 永久网站在线| 国内毛片毛片毛片毛片毛片| 最近最新中文字幕大全电影3| 午夜福利高清视频| 亚洲黑人精品在线| 国产色婷婷99| 美女大奶头视频| 日韩欧美精品免费久久| 午夜免费男女啪啪视频观看 | 97人妻精品一区二区三区麻豆| 99riav亚洲国产免费| 国产激情偷乱视频一区二区| 国产 一区精品| 美女cb高潮喷水在线观看| 国产午夜精品论理片| 午夜a级毛片| 国产人妻一区二区三区在| 黄色一级大片看看| 亚洲自偷自拍三级| 少妇人妻精品综合一区二区 | 久久亚洲精品不卡| 99精品久久久久人妻精品| 精品福利观看| 国国产精品蜜臀av免费| 亚洲精品乱码久久久v下载方式| 色哟哟哟哟哟哟| 日本免费a在线| 欧美最黄视频在线播放免费| 少妇人妻一区二区三区视频| 啪啪无遮挡十八禁网站| 国产伦精品一区二区三区四那| 黄片wwwwww| 精品人妻熟女av久视频| 精品福利观看| 亚洲最大成人av| av在线老鸭窝| 99国产极品粉嫩在线观看| 美女cb高潮喷水在线观看| 人妻制服诱惑在线中文字幕| 亚洲av美国av| 特大巨黑吊av在线直播| 又粗又爽又猛毛片免费看| 国产精品久久视频播放| 国产亚洲精品综合一区在线观看| 少妇丰满av| 精品无人区乱码1区二区| 小说图片视频综合网站| 精品乱码久久久久久99久播| 桃红色精品国产亚洲av| 国产精品亚洲美女久久久| 长腿黑丝高跟| 桃色一区二区三区在线观看| 色哟哟哟哟哟哟| 99国产精品一区二区蜜桃av| 天堂av国产一区二区熟女人妻| 国产蜜桃级精品一区二区三区| 成年免费大片在线观看| 国产女主播在线喷水免费视频网站 | 欧美成人a在线观看| 色5月婷婷丁香| 亚洲内射少妇av| 亚洲国产精品久久男人天堂| 99国产极品粉嫩在线观看| 久久国内精品自在自线图片| 色综合色国产| 欧美一区二区亚洲| 中文字幕久久专区| 国产伦精品一区二区三区视频9| 内射极品少妇av片p| 国产高清不卡午夜福利| 免费大片18禁| 日韩一区二区视频免费看| 久久久久久伊人网av| 久久精品国产99精品国产亚洲性色| 欧美最新免费一区二区三区| 国内毛片毛片毛片毛片毛片| 精品国内亚洲2022精品成人| 久久久久久伊人网av| 欧美成人免费av一区二区三区| 午夜免费成人在线视频| 国产人妻一区二区三区在| 亚洲中文字幕日韩| 丝袜美腿在线中文| 国产伦一二天堂av在线观看| 日本精品一区二区三区蜜桃| 一卡2卡三卡四卡精品乱码亚洲| 自拍偷自拍亚洲精品老妇| 国产人妻一区二区三区在| 九九在线视频观看精品| 最近中文字幕高清免费大全6 | 人妻制服诱惑在线中文字幕| 丰满人妻一区二区三区视频av| 桃色一区二区三区在线观看| 在线a可以看的网站| 女人十人毛片免费观看3o分钟| 赤兔流量卡办理| 日韩欧美三级三区| 亚洲国产精品久久男人天堂| 2021天堂中文幕一二区在线观| 可以在线观看的亚洲视频| 国内精品美女久久久久久| 亚洲精品影视一区二区三区av| 免费看光身美女| 一级毛片久久久久久久久女| 欧美潮喷喷水| 丰满人妻一区二区三区视频av| 色尼玛亚洲综合影院| 色噜噜av男人的天堂激情| 啦啦啦啦在线视频资源| 97碰自拍视频| 国产人妻一区二区三区在| 1000部很黄的大片| 欧美黑人巨大hd| 成人精品一区二区免费| 床上黄色一级片| 免费人成在线观看视频色| 国产欧美日韩精品亚洲av| 成人一区二区视频在线观看| 美女高潮喷水抽搐中文字幕| 日本三级黄在线观看| 村上凉子中文字幕在线| 亚洲欧美激情综合另类| 此物有八面人人有两片| 久久久国产成人免费| 日本黄大片高清| 国内精品宾馆在线| 日本-黄色视频高清免费观看| 最近中文字幕高清免费大全6 | 免费看a级黄色片| 日本撒尿小便嘘嘘汇集6| 欧美又色又爽又黄视频| 国产精品久久久久久av不卡| 免费看光身美女| 亚洲狠狠婷婷综合久久图片| 色综合婷婷激情| 琪琪午夜伦伦电影理论片6080| 成人一区二区视频在线观看| 亚洲国产欧洲综合997久久,| 国产精品一区二区免费欧美| 在线免费观看的www视频| 国产aⅴ精品一区二区三区波| 日本 av在线| 亚洲男人的天堂狠狠| 日本 av在线| 人人妻人人看人人澡| 波多野结衣高清作品| 国产色婷婷99| 久久精品影院6| 欧美一级a爱片免费观看看| 久久精品影院6| 成人一区二区视频在线观看| 国产av不卡久久| 欧美最黄视频在线播放免费| 国产亚洲91精品色在线| 成人美女网站在线观看视频| 国产成年人精品一区二区| 亚洲人成网站高清观看| 亚洲自偷自拍三级| 免费在线观看日本一区| 长腿黑丝高跟| 日日摸夜夜添夜夜添小说| 色综合站精品国产| 国产精品人妻久久久久久| 黄色视频,在线免费观看| 尾随美女入室| 嫩草影院新地址| 久久久久免费精品人妻一区二区| 黄色视频,在线免费观看| 国产精品亚洲美女久久久| 亚洲av第一区精品v没综合| 亚洲欧美日韩东京热| 国内毛片毛片毛片毛片毛片| 在线观看舔阴道视频| 一个人免费在线观看电影| 黄色配什么色好看| 一个人观看的视频www高清免费观看| 欧美精品啪啪一区二区三区| 日韩欧美免费精品| 99热精品在线国产| 国产伦精品一区二区三区视频9| 欧美另类亚洲清纯唯美| 国产伦精品一区二区三区视频9| 欧美高清性xxxxhd video| 亚洲美女视频黄频| 国产视频内射| 永久网站在线| 国产亚洲av嫩草精品影院| 日韩中文字幕欧美一区二区| 两个人视频免费观看高清| 99精品久久久久人妻精品| 在线观看av片永久免费下载| 人妻少妇偷人精品九色| 大型黄色视频在线免费观看| 韩国av在线不卡| 搡老岳熟女国产| 精品99又大又爽又粗少妇毛片 | 欧美丝袜亚洲另类 | 国产欧美日韩一区二区精品| ponron亚洲| 99热这里只有是精品在线观看| 精品久久久久久,| 成熟少妇高潮喷水视频| 免费av毛片视频| 国内精品久久久久精免费| 制服丝袜大香蕉在线| 免费看光身美女| 性插视频无遮挡在线免费观看| 欧美日韩综合久久久久久 | 日本一本二区三区精品| 亚洲精品粉嫩美女一区| 国产亚洲精品综合一区在线观看| 免费观看人在逋| 婷婷亚洲欧美| 午夜老司机福利剧场| 动漫黄色视频在线观看| 日韩 亚洲 欧美在线| 淫秽高清视频在线观看| 伦精品一区二区三区| 狠狠狠狠99中文字幕| 免费av不卡在线播放| 久久久久性生活片| 男人舔女人下体高潮全视频| 一本精品99久久精品77| 欧美精品国产亚洲| 在线播放无遮挡|