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

    基于布谷鳥搜索算法的水質(zhì)污染擴散仿真

    2022-04-19 00:47:26唐林旺李雪蓉麻欣宇王俊峰
    計算機仿真 2022年3期
    關(guān)鍵詞:河渠布谷鳥擴散系數(shù)

    唐林旺,李雪蓉,麻欣宇,王俊峰,2*

    (1. 江西理工大學(xué)資源與環(huán)境工程學(xué)院,江西 贛州 341000;2.江西省礦冶環(huán)境污染控制重點實驗室,江西 贛州 341000)

    1 引言

    水資源作為人類以及自然界中所有生物的生存根本,其用途廣泛,無法代替,同時也最容易受到外界污染。水體污染會導(dǎo)致圍繞水體生存的生態(tài)群體受到影響。且會導(dǎo)致污染所處區(qū)域的水文條件和污染物的種類、強度、時間、地點等各方面出現(xiàn)較強的不穩(wěn)定性,且污染物隨時處于變化的狀態(tài),隨著水體的流動發(fā)生遷移難以掌控,會對環(huán)境、經(jīng)濟以及社會造成嚴(yán)重的破壞性,甚至甚至?xí)ㄟ^生態(tài)鏈危害人類的身體健康[1-2]。

    為此本文提出一種基于CS(Cuckoo Search)算法的水質(zhì)污染垂向擴散跟蹤建模方法,是基于布谷鳥的繁衍方式以及萊維飛行演變所得到的算法,因布谷鳥會把鳥蛋生在宿主鳥的巢內(nèi),若被宿主鳥發(fā)現(xiàn),就會發(fā)生沖突,而宿主鳥事后發(fā)現(xiàn)布谷鳥蛋,那么會在其它地方建立一個鳥巢或者將其鳥蛋推出巢穴之外;而布谷鳥的蛋一旦被宿主鳥孵化,會將宿主鳥的蛋推出巢穴外,本文根據(jù)該特性,隨機地輸入數(shù)據(jù),尋找最優(yōu)二維水質(zhì)水流模型,隨后憑借計算橫向系數(shù),利用垂向系數(shù)觀察污染物對于水體垂向污染的影響,完成最終的跟蹤建模。

    2 CS算法的應(yīng)用計算原理

    將CS算法做出三種理想假設(shè),具體如下所示:

    1)布谷鳥隨機生產(chǎn)在某宿主鳥巢內(nèi);

    2)布谷鳥蛋會被孵化。

    3)布谷鳥能夠使用到的宿主鳥巢穴個數(shù)n為固定的,而布谷鳥鳥蛋在宿主鳥巢穴內(nèi)被發(fā)現(xiàn)概率為pa∈[0,1]。

    以上述三種假設(shè)作為前提,可以得到布谷鳥尋找鳥巢的位置以及路徑,具體公式為:

    (1)

    萊維飛行連續(xù)的跳躍路徑和時間關(guān)系是服從萊維分布情況[3],通過對分布函數(shù)進行簡化以及傅里葉變換之后,獲得冪次的行為概率密度函數(shù),具體可得公式為:

    Le′vy~u=t-λ;1<λ<3

    (2)

    上式中:λ代表冪次系數(shù)。

    由于通過模擬萊維飛行的跳躍路徑,可以獲得具體公式為:

    (3)

    上式中:s代表萊維飛行跳躍的路徑Le′vy;參數(shù)Le′vy(λ)與式(2)內(nèi)λ關(guān)系是λ=1+β,β的取值范圍是0<β<2,能夠在CS算法內(nèi)取β=1.5,參數(shù)o、j代表正態(tài)的分布隨機數(shù),可以服從式(4)正態(tài)的分布情況,在式(4)相應(yīng)正態(tài)的分布標(biāo)準(zhǔn)差為σo、σj取值,具體如式(5)所示:

    (4)

    式中,

    (5)

    布谷鳥按照萊維飛行方式,其搜索路徑是隨意改變的,此時會從一個區(qū)域轉(zhuǎn)換至另外一個區(qū)域,使其整體尋優(yōu)的能力增強,此外,CS方法憑借布谷鳥的生存方式,能夠?qū)ζ溥M行定義,在布谷鳥被宿主鳥所發(fā)現(xiàn)的概率為pa=0.25,以此可以適應(yīng)環(huán)境的布谷鳥則被孵化,反之,則會被淘汰掉,以此剩下最佳的布谷鳥。本文憑借CS方法的強收斂性對整體水域水流流動方式進行搜索,尋找最優(yōu)的流動路徑以后,通過構(gòu)建污染跟蹤模型,尋找出污染物的擴散情況[4]。

    3 水質(zhì)污染問題的垂向擴散跟蹤建模方法。

    3.1 水質(zhì)污染垂向擴散系數(shù)

    水質(zhì)污染的垂向擴散系數(shù)是衡量污染物在河渠內(nèi)運輸能力的重要參數(shù)[5],采用對流速分布的函數(shù)形式,以此能夠在各向同性的紊動條件之下,得出污染物垂向擴散的系數(shù)。

    以混合區(qū)域外邊界的標(biāo)準(zhǔn)曲線特征尺度基礎(chǔ),獲取順直寬闊的水體垂向擴散系數(shù),具體公式為:

    (6)

    上式中:Ez代表垂向的擴散系數(shù),單位為:m2/s,Ls代表河渠最大的寬度,單位為m,Sz代表河渠的垂直方向面積,單位為m2,e代表岸邊污染混合區(qū)域的最大長度以及對應(yīng)縱坐標(biāo)比值。

    在式(6)內(nèi)表明,河渠的垂直方向平均流速u恒定條件下,垂直方向上面積Sz越大,說明垂向擴散的系數(shù)Ez就越大,而最大的長度Ls越大,那么垂向擴散的系數(shù)Ez就越小[6]。

    3.2 二維污染物垂向跟蹤模型

    河渠的一維水流水質(zhì)模型,一般情況下要對其進行理想化以及大量簡化的操作,而實際過程內(nèi),要模擬水體種類有很大幾率不滿足一維模型的理想狀態(tài),例如:水平寬度較大或者水平尺度較短、橫向濃度的梯度較大等問題,這時不用對水力水質(zhì)參數(shù)沿著垂向變化進行考慮,可采用沿著水深均勻流動量來代表垂直方向流場,因此根據(jù)CS算法的應(yīng)用原理,本文采用CS算法獲取二維水流水質(zhì)的模型,該模型能夠跟蹤水體內(nèi)水流垂向與污染物濃度變化,從而完成跟蹤建模[7]。

    3.2.1 水動力模型

    在河渠內(nèi)獲得從水底至水面微小的水體。h0、z0分別代表某一個基準(zhǔn)面下河渠水位以及底部的高程,Ux、Uy、Uz分別代表x、y、z方向上沿著水深的平均流速,u、v、w分別代表x、y、z方向沿著水深時的平均流速,具體公式為:

    (7)

    自由表面和底部的運動學(xué)條件公式為:

    (8)

    (9)

    1)沿著水深均勻連續(xù)性的方程

    水流的連續(xù)方程公式為:

    (10)

    上式中:q代表單位面積上污染物進出河渠流量,其中,流入即為正,而流出為負(fù)[8]。

    2)沿著水深的平均運動方程

    (11)

    設(shè)定u、v分別代表河渠內(nèi)x、y方向流速分量,那么河渠x方向動量的輸運方程公式為:

    (12)

    與河渠x方向動量的輸運方程相似,河渠y方向上動量的輸運方程公式為:

    (13)

    河渠湍流動能以及湍流耗散輸運的方程公式為:

    (14)

    湍流動能的產(chǎn)生項公式為:

    (15)

    上式中:g代表重力的加速度,n代表河渠的底部糙率,k代表河渠的深度湍流動能,ε代表耗散率,μ、代表分子動力粘性的系數(shù),μt代表漩渦運動粘性的系數(shù),γeff代表有效的粘性系數(shù),C代表河渠輸送水質(zhì)變量的濃度,單位是g/m3,SC代表污染物強度,Sk代表河渠內(nèi)污染物的衰弱減項[9]。

    3.2.2 污染物跟蹤模型

    污染物在進入水體內(nèi),會發(fā)生擴散,在任意的dt時間段中,所有因素一起作用,會引發(fā)微分單元體內(nèi)的污染物質(zhì)量增量,包括水平面(x,y)方向的變化量。因此要構(gòu)建水體內(nèi)污染物的輸運微分方程,就需要許多因素一起作用下,所引發(fā)單元體積中,污染物的變化量與該水體中污染物處于dt時間段中的變化量[10]。

    1)在水流的作用下,處于x、y方向所引發(fā)污染物的增量通量,具體公式為:

    (16)

    2)在分子的分散作用下,處于x、y方向所引發(fā)污染物的增量通量,具體公式為:

    (17)

    3)在紊動擴散作用下,處于x、y方向所引發(fā)污染物的增量通量[11],具體公式為:

    (18)

    4)在彌散作用下,處于x、y方向所引發(fā)污染物的增量通量,具體公式為:

    (19)

    5)在微元體源項Sv和水體衰減作用下所引發(fā)污染物的增量通量,具體公式為:

    Sv|(x,y)dxdydt+Sk|(x,y)dxdydt

    (20)

    依據(jù)質(zhì)量守恒定律,獲得河渠二維污染垂向跟蹤的轉(zhuǎn)換基本模型,具體公式為:

    (21)

    上式中:t代表污染物的運移時間,單位是h代表河渠跟蹤高度,Ex、Ey分別代表河渠x,y方向上離散系數(shù)[12],其單位為m2/s。

    4 仿真實驗設(shè)計

    4.1 實驗一

    為驗證本文方法的有效性,選擇實驗環(huán)境,如表1所示:

    表1 實驗環(huán)境

    通過模擬測試,設(shè)定水面污染源的強度為:m=50g/m2,橫向以及垂向系數(shù)分別為:0.1m2/s、0.02m2/s。擴散時間為t=1000s,水深20m,在該次實驗分成三次進行,不考慮橫向污染的影響,僅觀察垂向污染擴散情況,具體結(jié)果如圖1所示:

    圖1 模擬污染物垂向擴散系數(shù)

    通過觀察圖1能夠看出,分別進行三次不同的污染濃度測試過程中,第一次污染物濃度較低,當(dāng)擴散至大約2米的位置,濃度擴散系數(shù)逐漸降低,而達到5米的位置,基本停止向下擴散,開始隨著時間的流失,橫向擴散,由于橫向擴散跟蹤系數(shù),不在本文實驗考慮范圍之內(nèi),因此第一次終止;第二次提升污染物濃度為第一次的一倍,在幾秒鐘內(nèi),擴散得深度就超越了第一次,然后隨著時間的推移,在約9米處擴散系數(shù)降低,由于該次濃度較大,直到擴散約12米處才基本停止向下擴散;第三次,再次提升一倍的污染物濃度,但污染物并沒有像第二次以很快的速度超越第一次一樣,超越第二次測試,而是在達到10米的時候,出現(xiàn)勻速下降,直到14米的時間,擴散系數(shù)降低;第四次濃度為12mg·L-1到15米的時候擴散系數(shù)降低,16米就已經(jīng)停止擴散。總結(jié)來說,污染物濃度越大,垂直擴散越深,隨著深度提升擴散系數(shù)會逐步降低,對污染物的影響力就越小。

    4.2 實驗二

    為進一步驗證本文方法的建模效果,通過實際污染擴散情況進行對比,在相同條件,觀察二者之間垂向擴散的變化情況,具體如圖2所示:

    圖2 本文跟蹤結(jié)果與實際跟蹤結(jié)果對比

    通過觀察圖2能夠看出,本文方法跟蹤建模的模擬情況與實際污染擴散情況,它們的曲線波動,基本相同、差距較小,都是隨著時間的增加,擴散的越深。這是由于本文通過CS算法可以找出水體流動的方式,再以污染物隨著水體流動的方式,找出污染物的擴散系數(shù),完成有效跟蹤。

    5 結(jié)束語

    本文采用CS算法獲取二維水流水質(zhì)的模型,跟蹤水體內(nèi)水流垂向與污染物濃度變化,

    根據(jù)混合區(qū)域外邊界的標(biāo)準(zhǔn)曲線特征尺度獲取流速分布函數(shù)計算水體垂向擴散系數(shù)。結(jié)合水體污染垂直方向上的流場,從而完成污染擴散跟蹤。最后設(shè)計仿真實驗,實驗分為兩部分。實驗首先獲取了模擬污染物垂向擴散系數(shù),基于此結(jié)果,將模擬的實驗結(jié)果與實際結(jié)果對比,二者之間的擴散系數(shù)基本相同,以此能夠證明本文方法能夠符合現(xiàn)實的精準(zhǔn)跟蹤需求。

    不過本文并未考慮特殊的環(huán)境場所,例如,水流速度極快的河渠、攜帶沙土的河渠等,因此,未來本文還需要進一步研究,爭取跟蹤的范圍更廣,構(gòu)建更精準(zhǔn)的模型。

    猜你喜歡
    河渠布谷鳥擴散系數(shù)
    布谷鳥讀信
    布谷鳥讀信
    西夏河渠技術(shù)述論
    西夏學(xué)(2020年2期)2020-01-24 07:42:52
    噓!布谷鳥來了
    大灰狼(2019年4期)2019-05-14 16:38:38
    河北省寧晉縣河渠鎮(zhèn)中心小學(xué)學(xué)生作品
    植被覆蓋與生態(tài)發(fā)展的研究
    布谷鳥叫醒的清晨
    基于河渠和坑塘聯(lián)通的雨洪資源綜合利用研究
    基于Sauer-Freise 方法的Co- Mn 體系fcc 相互擴散系數(shù)的研究
    上海金屬(2015年5期)2015-11-29 01:13:59
    FCC Ni-Cu 及Ni-Mn 合金互擴散系數(shù)測定
    上海金屬(2015年6期)2015-11-29 01:09:09
    久久精品国产亚洲av香蕉五月| 精品午夜福利视频在线观看一区| 久久精品亚洲精品国产色婷小说| 久久久久亚洲av毛片大全| 国产精品免费一区二区三区在线| 色播亚洲综合网| 在线av久久热| 中文字幕精品免费在线观看视频| 99久久久亚洲精品蜜臀av| 午夜免费观看网址| 亚洲国产欧美一区二区综合| 国产区一区二久久| 啦啦啦韩国在线观看视频| 国产精品一区二区精品视频观看| 美女免费视频网站| 久久久精品欧美日韩精品| 一区在线观看完整版| 91大片在线观看| 欧美精品啪啪一区二区三区| 亚洲自拍偷在线| 岛国视频午夜一区免费看| 久久天躁狠狠躁夜夜2o2o| 成熟少妇高潮喷水视频| 亚洲色图av天堂| 搡老岳熟女国产| 丝袜美腿诱惑在线| 久久精品国产99精品国产亚洲性色 | 精品日产1卡2卡| 亚洲午夜理论影院| 国产成人系列免费观看| 久久久久久久精品吃奶| 多毛熟女@视频| 又黄又爽又免费观看的视频| 久久精品aⅴ一区二区三区四区| 亚洲九九香蕉| 中文字幕色久视频| 成人国产综合亚洲| 丰满的人妻完整版| 亚洲欧美日韩另类电影网站| 国产在线观看jvid| 啪啪无遮挡十八禁网站| 伦理电影免费视频| 身体一侧抽搐| 国产一区二区三区综合在线观看| 18禁黄网站禁片午夜丰满| 脱女人内裤的视频| 久久婷婷成人综合色麻豆| 亚洲第一青青草原| 国产熟女午夜一区二区三区| 亚洲情色 制服丝袜| 亚洲欧美日韩另类电影网站| 中文字幕人妻丝袜一区二区| 91国产中文字幕| 欧美大码av| 精品不卡国产一区二区三区| 一二三四社区在线视频社区8| 禁无遮挡网站| 最近最新中文字幕大全电影3 | 在线观看www视频免费| 亚洲激情在线av| www.www免费av| 国产人伦9x9x在线观看| 999久久久精品免费观看国产| 国产精品免费一区二区三区在线| 欧美日韩福利视频一区二区| 国产亚洲精品一区二区www| 精品人妻1区二区| 天天添夜夜摸| 精品日产1卡2卡| 美女国产高潮福利片在线看| 国产又色又爽无遮挡免费看| a级毛片在线看网站| √禁漫天堂资源中文www| 免费无遮挡裸体视频| 日本欧美视频一区| 男女下面插进去视频免费观看| 亚洲欧美精品综合久久99| 不卡av一区二区三区| 国产精品久久久av美女十八| 国产精品久久视频播放| 免费在线观看日本一区| 日本五十路高清| 日韩视频一区二区在线观看| 性少妇av在线| 一本大道久久a久久精品| 亚洲成人精品中文字幕电影| 亚洲精华国产精华精| 在线播放国产精品三级| 欧美亚洲日本最大视频资源| av有码第一页| e午夜精品久久久久久久| 制服人妻中文乱码| 美女国产高潮福利片在线看| 精品一品国产午夜福利视频| 国产精品久久视频播放| 天堂影院成人在线观看| 99久久久亚洲精品蜜臀av| 国产一区二区在线av高清观看| 免费在线观看黄色视频的| 麻豆一二三区av精品| 波多野结衣高清无吗| 成人亚洲精品一区在线观看| 最新在线观看一区二区三区| 又黄又粗又硬又大视频| 女人被躁到高潮嗷嗷叫费观| 亚洲国产精品999在线| 黄频高清免费视频| 日本三级黄在线观看| 久久国产精品影院| 十八禁网站免费在线| 欧美人与性动交α欧美精品济南到| 日日爽夜夜爽网站| 亚洲国产精品sss在线观看| 欧美日韩乱码在线| 国产人伦9x9x在线观看| 午夜视频精品福利| 亚洲av熟女| 欧洲精品卡2卡3卡4卡5卡区| 亚洲成人久久性| 国产午夜精品久久久久久| 久久久精品国产亚洲av高清涩受| 搞女人的毛片| 日韩欧美一区视频在线观看| 国产av又大| 日日夜夜操网爽| 人人妻,人人澡人人爽秒播| 制服丝袜大香蕉在线| 亚洲国产中文字幕在线视频| 亚洲成人久久性| 久久久精品欧美日韩精品| 狂野欧美激情性xxxx| 欧美不卡视频在线免费观看 | 国产三级在线视频| 中文字幕另类日韩欧美亚洲嫩草| 亚洲欧洲精品一区二区精品久久久| 成人特级黄色片久久久久久久| www.999成人在线观看| 在线观看日韩欧美| 搡老熟女国产l中国老女人| 巨乳人妻的诱惑在线观看| 久久人妻熟女aⅴ| www.精华液| 999久久久精品免费观看国产| 久久久国产精品麻豆| www国产在线视频色| 国产精品1区2区在线观看.| 亚洲av美国av| 最新在线观看一区二区三区| 欧美成人午夜精品| 久久久精品欧美日韩精品| 99久久综合精品五月天人人| 午夜两性在线视频| 一个人免费在线观看的高清视频| 色尼玛亚洲综合影院| 不卡一级毛片| 高清黄色对白视频在线免费看| 色尼玛亚洲综合影院| 精品一区二区三区视频在线观看免费| 黑丝袜美女国产一区| 亚洲av电影不卡..在线观看| 亚洲av成人一区二区三| 亚洲电影在线观看av| www.999成人在线观看| 很黄的视频免费| 亚洲国产精品合色在线| 女生性感内裤真人,穿戴方法视频| 亚洲一区高清亚洲精品| 动漫黄色视频在线观看| а√天堂www在线а√下载| 国产高清videossex| 久久精品91蜜桃| 午夜福利免费观看在线| 欧洲精品卡2卡3卡4卡5卡区| 免费在线观看影片大全网站| 亚洲一码二码三码区别大吗| 一级毛片女人18水好多| 人人妻人人澡欧美一区二区 | 精品一区二区三区视频在线观看免费| 午夜福利免费观看在线| 国产成人啪精品午夜网站| 欧美日韩黄片免| 俄罗斯特黄特色一大片| 91大片在线观看| 大型黄色视频在线免费观看| 久久精品91无色码中文字幕| 男男h啪啪无遮挡| 国产在线观看jvid| 欧美日韩一级在线毛片| 色婷婷久久久亚洲欧美| 19禁男女啪啪无遮挡网站| 久久久久久亚洲精品国产蜜桃av| av免费在线观看网站| 在线天堂中文资源库| 两人在一起打扑克的视频| 欧美日韩精品网址| 亚洲人成77777在线视频| 色av中文字幕| 亚洲一区二区三区色噜噜| 欧美日韩亚洲综合一区二区三区_| 免费久久久久久久精品成人欧美视频| 男人舔女人的私密视频| а√天堂www在线а√下载| 久久久久久久久免费视频了| 精品一品国产午夜福利视频| 精品欧美国产一区二区三| 午夜福利影视在线免费观看| 国产欧美日韩精品亚洲av| 欧美丝袜亚洲另类 | 在线观看免费日韩欧美大片| 亚洲va日本ⅴa欧美va伊人久久| 宅男免费午夜| 久久精品国产99精品国产亚洲性色 | 亚洲精品久久成人aⅴ小说| 制服诱惑二区| 超碰成人久久| 精品久久久久久久毛片微露脸| 美女扒开内裤让男人捅视频| 欧美+亚洲+日韩+国产| 69av精品久久久久久| 欧美一级a爱片免费观看看 | 波多野结衣av一区二区av| 亚洲精品久久国产高清桃花| 亚洲中文日韩欧美视频| 午夜福利欧美成人| 老司机午夜十八禁免费视频| 久久精品91蜜桃| 日本免费a在线| 国产精品一区二区免费欧美| 巨乳人妻的诱惑在线观看| 别揉我奶头~嗯~啊~动态视频| 欧美性长视频在线观看| 久久久国产精品麻豆| 欧美丝袜亚洲另类 | 久99久视频精品免费| 天天一区二区日本电影三级 | 国产伦一二天堂av在线观看| 午夜免费观看网址| 精品午夜福利视频在线观看一区| 国产av在哪里看| 视频在线观看一区二区三区| 免费在线观看视频国产中文字幕亚洲| 亚洲成av片中文字幕在线观看| 久久这里只有精品19| 午夜激情av网站| 亚洲av日韩精品久久久久久密| 高清在线国产一区| 99re在线观看精品视频| 亚洲人成伊人成综合网2020| 不卡一级毛片| 99riav亚洲国产免费| 老司机深夜福利视频在线观看| 日本精品一区二区三区蜜桃| 老司机午夜福利在线观看视频| 亚洲一区二区三区色噜噜| 亚洲av熟女| 一区二区三区国产精品乱码| 免费高清视频大片| 欧美国产日韩亚洲一区| 亚洲欧美日韩无卡精品| 国产麻豆69| 国产午夜福利久久久久久| 亚洲av熟女| 国产亚洲精品综合一区在线观看 | 久久亚洲精品不卡| 操美女的视频在线观看| 大型黄色视频在线免费观看| 国产熟女午夜一区二区三区| 激情视频va一区二区三区| 亚洲欧美一区二区三区黑人| 亚洲国产日韩欧美精品在线观看 | 两性夫妻黄色片| 亚洲午夜精品一区,二区,三区| 99国产精品免费福利视频| 午夜久久久在线观看| av中文乱码字幕在线| 欧美乱妇无乱码| 国产精品99久久99久久久不卡| 俄罗斯特黄特色一大片| 国产成人精品久久二区二区免费| 亚洲精品粉嫩美女一区| 高潮久久久久久久久久久不卡| 一a级毛片在线观看| 非洲黑人性xxxx精品又粗又长| 久久香蕉精品热| 国产高清有码在线观看视频 | 免费在线观看亚洲国产| 久久久久久人人人人人| 午夜精品在线福利| 在线观看www视频免费| av网站免费在线观看视频| 国产熟女xx| 日韩三级视频一区二区三区| 夜夜躁狠狠躁天天躁| 日本 av在线| 两人在一起打扑克的视频| 久久久久久人人人人人| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲国产高清在线一区二区三 | 一区在线观看完整版| 精品少妇一区二区三区视频日本电影| 久久天躁狠狠躁夜夜2o2o| 99国产精品一区二区三区| 亚洲狠狠婷婷综合久久图片| 国产成人免费无遮挡视频| 91麻豆精品激情在线观看国产| 精品午夜福利视频在线观看一区| 国产在线观看jvid| www国产在线视频色| 日韩精品青青久久久久久| 国产色视频综合| 国产精品一区二区精品视频观看| 夜夜爽天天搞| 欧美国产精品va在线观看不卡| 亚洲av日韩精品久久久久久密| 精品一区二区三区av网在线观看| 国产欧美日韩一区二区三| 精品午夜福利视频在线观看一区| 免费少妇av软件| 免费在线观看亚洲国产| 亚洲精品在线观看二区| x7x7x7水蜜桃| 999久久久国产精品视频| 国产av又大| 亚洲成人免费电影在线观看| 久久久久九九精品影院| 日本黄色视频三级网站网址| 两性夫妻黄色片| 久久香蕉精品热| 色精品久久人妻99蜜桃| 亚洲一区中文字幕在线| 亚洲黑人精品在线| 国产精品98久久久久久宅男小说| 亚洲国产欧美一区二区综合| 变态另类成人亚洲欧美熟女 | 欧美日韩黄片免| 十八禁人妻一区二区| 国产精品一区二区精品视频观看| 久久天躁狠狠躁夜夜2o2o| 欧美不卡视频在线免费观看 | 90打野战视频偷拍视频| 黑人操中国人逼视频| 美女高潮到喷水免费观看| 看免费av毛片| 国产午夜精品久久久久久| 人妻丰满熟妇av一区二区三区| 欧美午夜高清在线| 久久热在线av| 88av欧美| 国产99久久九九免费精品| 日本精品一区二区三区蜜桃| 中文字幕久久专区| 老汉色av国产亚洲站长工具| 久久 成人 亚洲| 19禁男女啪啪无遮挡网站| 一级毛片女人18水好多| 亚洲中文av在线| 日日爽夜夜爽网站| 露出奶头的视频| 久久国产亚洲av麻豆专区| 老汉色av国产亚洲站长工具| 国产1区2区3区精品| 免费高清视频大片| 免费一级毛片在线播放高清视频 | 欧美国产精品va在线观看不卡| 国产伦人伦偷精品视频| 天天一区二区日本电影三级 | 亚洲精品在线美女| 69精品国产乱码久久久| 成人国产一区最新在线观看| 操美女的视频在线观看| 97人妻精品一区二区三区麻豆 | 欧美大码av| 国产91精品成人一区二区三区| 级片在线观看| 岛国视频午夜一区免费看| 黄片播放在线免费| 啦啦啦 在线观看视频| 久久青草综合色| 国产一区二区三区综合在线观看| 99国产精品99久久久久| 俄罗斯特黄特色一大片| www日本在线高清视频| 国产成人精品久久二区二区91| 亚洲自偷自拍图片 自拍| 搡老妇女老女人老熟妇| 欧美绝顶高潮抽搐喷水| 亚洲一码二码三码区别大吗| 变态另类成人亚洲欧美熟女 | 波多野结衣高清无吗| 亚洲人成伊人成综合网2020| 亚洲国产毛片av蜜桃av| 国产激情久久老熟女| 久久精品国产亚洲av香蕉五月| 丰满人妻熟妇乱又伦精品不卡| 日韩欧美一区二区三区在线观看| 久久午夜综合久久蜜桃| 天堂√8在线中文| 啦啦啦韩国在线观看视频| 国产三级黄色录像| 亚洲午夜理论影院| 亚洲精品中文字幕一二三四区| 岛国在线观看网站| 国产欧美日韩综合在线一区二区| 久久欧美精品欧美久久欧美| 亚洲成国产人片在线观看| 国产又爽黄色视频| 国产av一区在线观看免费| 国产精品久久久人人做人人爽| 高清在线国产一区| 亚洲精品美女久久久久99蜜臀| 日韩大尺度精品在线看网址 | 超碰成人久久| 国产精品 国内视频| 国语自产精品视频在线第100页| a级毛片在线看网站| 国产精品亚洲一级av第二区| 精品欧美国产一区二区三| ponron亚洲| 非洲黑人性xxxx精品又粗又长| 亚洲精品一区av在线观看| www国产在线视频色| av在线播放免费不卡| 丰满人妻熟妇乱又伦精品不卡| 18禁国产床啪视频网站| 精品国内亚洲2022精品成人| 久久久精品欧美日韩精品| 欧美绝顶高潮抽搐喷水| 1024视频免费在线观看| 国产一级毛片七仙女欲春2 | netflix在线观看网站| 日本免费a在线| 亚洲午夜精品一区,二区,三区| 久久精品91无色码中文字幕| 久久婷婷成人综合色麻豆| 妹子高潮喷水视频| 桃红色精品国产亚洲av| 亚洲 国产 在线| 十分钟在线观看高清视频www| 国产在线观看jvid| 精品久久蜜臀av无| 亚洲激情在线av| 免费在线观看日本一区| 国产精品av久久久久免费| 欧美性长视频在线观看| 99riav亚洲国产免费| 免费在线观看影片大全网站| 国产区一区二久久| 两个人免费观看高清视频| av天堂久久9| 亚洲久久久国产精品| 看黄色毛片网站| 国产伦人伦偷精品视频| 在线观看免费视频日本深夜| 在线观看日韩欧美| 亚洲国产精品sss在线观看| 国产精品 欧美亚洲| 欧美激情 高清一区二区三区| 啦啦啦观看免费观看视频高清 | 曰老女人黄片| 在线观看免费日韩欧美大片| 波多野结衣高清无吗| 精品欧美一区二区三区在线| 国产91精品成人一区二区三区| 国产精品电影一区二区三区| 国内毛片毛片毛片毛片毛片| 国语自产精品视频在线第100页| 两性夫妻黄色片| 波多野结衣高清无吗| 色综合欧美亚洲国产小说| 国产蜜桃级精品一区二区三区| 精品久久久久久久久久免费视频| 免费高清视频大片| 人人澡人人妻人| 欧美老熟妇乱子伦牲交| 9色porny在线观看| 黑丝袜美女国产一区| 色老头精品视频在线观看| 国产成人影院久久av| 久久人人97超碰香蕉20202| 欧美中文日本在线观看视频| 午夜免费观看网址| 色综合婷婷激情| 男女下面插进去视频免费观看| 亚洲 欧美 日韩 在线 免费| 精品无人区乱码1区二区| 一a级毛片在线观看| 人人澡人人妻人| 天天添夜夜摸| 久久国产精品影院| 黄色片一级片一级黄色片| 亚洲国产毛片av蜜桃av| 国产日韩一区二区三区精品不卡| 亚洲狠狠婷婷综合久久图片| 国产精品免费视频内射| 女警被强在线播放| 精品国产一区二区三区四区第35| 夜夜躁狠狠躁天天躁| 久久中文字幕人妻熟女| 久久国产精品影院| 精品国产一区二区久久| 男女做爰动态图高潮gif福利片 | 欧美黄色片欧美黄色片| 国产熟女xx| 在线观看午夜福利视频| 波多野结衣高清无吗| 亚洲国产中文字幕在线视频| 亚洲精品国产精品久久久不卡| 在线十欧美十亚洲十日本专区| 日韩大码丰满熟妇| 国产精品影院久久| 国产精品永久免费网站| 琪琪午夜伦伦电影理论片6080| 久久久久国产精品人妻aⅴ院| 成人手机av| 日本五十路高清| 黄色丝袜av网址大全| 国产精品日韩av在线免费观看 | 久久九九热精品免费| 国产精品久久久久久亚洲av鲁大| 色尼玛亚洲综合影院| 久久久久久国产a免费观看| 美女国产高潮福利片在线看| av福利片在线| 国产一区在线观看成人免费| 看免费av毛片| 亚洲av熟女| 99久久久亚洲精品蜜臀av| 自拍欧美九色日韩亚洲蝌蚪91| 中文字幕人妻丝袜一区二区| 欧美性长视频在线观看| 好男人在线观看高清免费视频 | 男女下面插进去视频免费观看| 欧美亚洲日本最大视频资源| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av成人av| 超碰成人久久| 久9热在线精品视频| 啦啦啦免费观看视频1| 亚洲在线自拍视频| 一夜夜www| 999精品在线视频| 精品欧美一区二区三区在线| 欧美日本中文国产一区发布| 欧美日韩瑟瑟在线播放| 久久香蕉激情| 欧美黑人精品巨大| 91精品国产国语对白视频| 电影成人av| 此物有八面人人有两片| 欧美国产日韩亚洲一区| 日韩视频一区二区在线观看| 黄网站色视频无遮挡免费观看| 亚洲免费av在线视频| 88av欧美| 女性被躁到高潮视频| 亚洲人成电影免费在线| 久久精品aⅴ一区二区三区四区| 99国产精品99久久久久| 欧美成人一区二区免费高清观看 | 国产精品永久免费网站| 国产一区二区三区在线臀色熟女| 在线观看66精品国产| 日本撒尿小便嘘嘘汇集6| 久久精品国产清高在天天线| 成人永久免费在线观看视频| 午夜免费激情av| 日韩精品青青久久久久久| 天堂影院成人在线观看| 国产亚洲精品久久久久5区| 欧美亚洲日本最大视频资源| 欧美另类亚洲清纯唯美| 久久久国产精品麻豆| 国产亚洲av高清不卡| 女人被狂操c到高潮| www.999成人在线观看| 国产精品国产高清国产av| bbb黄色大片| 黄色女人牲交| 美国免费a级毛片| 午夜精品在线福利| 在线观看一区二区三区| 午夜日韩欧美国产| 久久性视频一级片| 精品国产一区二区三区四区第35| 精品久久久久久久人妻蜜臀av | 亚洲最大成人中文| 国产精品美女特级片免费视频播放器 | 无人区码免费观看不卡| 国产蜜桃级精品一区二区三区| 好看av亚洲va欧美ⅴa在| 欧美色视频一区免费| 日韩欧美一区视频在线观看| 91九色精品人成在线观看| 国产精华一区二区三区| 成人三级黄色视频| 成人欧美大片| 久久天躁狠狠躁夜夜2o2o| 亚洲美女黄片视频| 欧美最黄视频在线播放免费| av超薄肉色丝袜交足视频| 欧美日本亚洲视频在线播放| 久久青草综合色| 午夜福利在线观看吧| 亚洲欧美日韩无卡精品| 免费一级毛片在线播放高清视频 | 久久久久亚洲av毛片大全| 两个人看的免费小视频| 人成视频在线观看免费观看| 亚洲少妇的诱惑av| 免费观看精品视频网站| 午夜福利在线观看吧| 国产在线观看jvid| 亚洲中文日韩欧美视频| 成年人黄色毛片网站| 久久精品国产亚洲av香蕉五月|