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

    暴雨徑流對五里湖水質(zhì)的影響

    2015-06-06 12:40:43劉曉東胡功宇周媛媛涂琦樂鄭孝宇
    水資源保護(hù) 2015年4期
    關(guān)鍵詞:調(diào)水通量暴雨

    劉曉東,胡功宇,周媛媛,涂琦樂,鄭孝宇

    (1.淺水湖泊綜合治理與資源開發(fā)教育部重點實驗室,江蘇 南京 210098;2.河海大學(xué)環(huán)境學(xué)院,江蘇 南京 210098; 3.長江水利委員會荊江水文水資源勘測局,湖北 荊州 434002)

    ?

    暴雨徑流對五里湖水質(zhì)的影響

    劉曉東1,2,胡功宇3,周媛媛2,涂琦樂2,鄭孝宇2

    (1.淺水湖泊綜合治理與資源開發(fā)教育部重點實驗室,江蘇 南京 210098;2.河海大學(xué)環(huán)境學(xué)院,江蘇 南京 210098; 3.長江水利委員會荊江水文水資源勘測局,湖北 荊州 434002)

    針對五里湖曲折不規(guī)則邊界,采用四邊形網(wǎng)格擬合,構(gòu)建了基于四邊形網(wǎng)格的平面二維通量向量分裂(FVS)格式水量水質(zhì)耦合模型,并進(jìn)行了模型驗證。利用該模型模擬了設(shè)計暴雨徑流對五里湖水質(zhì)影響的動態(tài)過程,分析東西五里湖寶界橋斷面CODMn、NH3-N、TP和TN的交換通量。結(jié)果表明:在設(shè)計典型暴雨條件下,罵蠡港大流量的進(jìn)水對西五里湖的水質(zhì)有較大影響,設(shè)計暴雨發(fā)生3 d后,徑流污水開始影響到西五里湖,在第4天后水質(zhì)濃度達(dá)到峰值,使部分湖區(qū)水質(zhì)由Ⅲ類變?yōu)棰纛?但持續(xù)時間較短,此后濃度迅速減小。

    暴雨徑流;五里湖;二維水流水質(zhì)耦合模型;四邊形網(wǎng)格;CODMn

    五里湖位于江蘇省無錫市境內(nèi),是太湖北部的一個內(nèi)湖,為典型的平原淺水湖泊,湖泊東西長約6 km,南北寬0.3~1.5 km,平均水深約1.8 m,面積約為5.4 km2,沿線支流眾多,主要有罵蠡港、梁塘河、梁溪河、長廣溪等河道(圖1)。近年來,太湖富營養(yǎng)化現(xiàn)象日趨嚴(yán)重,尤其是五里湖,成為太湖流域水質(zhì)最差的水域,CODMn、TP、TN等水質(zhì)指標(biāo)均為劣V類,藻類頻繁暴發(fā)。加之五里湖水體流動緩慢、稀釋自凈能力差, 湖體面積小、流程短;同時清淤后湖底地形比較平坦,五里湖流態(tài)相對單一,很容易受到進(jìn)出湖吞吐流的影響,造成五里湖水環(huán)境嚴(yán)重惡化[1]。

    圖1 五里湖區(qū)域

    近年來,我國對五里湖和大太湖的水動力和水質(zhì)數(shù)值模擬陸續(xù)開展了不少研究工作。例如,中國科學(xué)院南京地理與湖泊研究所重點研究了梅梁湖泵站不同調(diào)水流量和不同出流方式下五里湖水環(huán)境的變化[2];田向榮等[3]以改善梅梁湖、五里湖水環(huán)境為目標(biāo),重點分析了調(diào)水的敏感影響因子對梅梁湖、五里湖水環(huán)境改善效果的影響;丁玲等[4]建立了五里湖二維水流水質(zhì)耦合模型,并對模型參數(shù)進(jìn)行率定,應(yīng)用該模型模擬五里湖在不同計算工況下的CODMn濃度,分析調(diào)貢湖水入五里湖時五里湖水質(zhì)的改善狀況。駱輝煌[5]研究了平面二維有限元水動力學(xué)模型RMA2和RMA4在五里湖調(diào)水?dāng)?shù)值模擬中主要參數(shù)的敏感性,提出通過集中調(diào)水和維持性調(diào)水相結(jié)合的方式改善五里湖水質(zhì)的優(yōu)化調(diào)水方案。

    上述研究基本上都是針對常規(guī)水文條件下五里湖的水質(zhì)狀況,未考慮暴雨等極端水文情況下的影響。本文采用基于四邊形網(wǎng)格的二維水流水質(zhì)耦合模型模擬暴雨徑流對五里湖水質(zhì)影響的動態(tài)過程,分析了寶界橋斷面的交換通量。

    1 平面二維通量向量分裂(FVS)格式水量水質(zhì)耦合模型的構(gòu)建

    1.1 基本方程

    二維水流水質(zhì)模型基本方程組[6-7]的守恒形式可表達(dá)為

    (1)

    其中,b1=0,b2=gh(S0x-Sfx),b3=bh(S0y-Sfy),

    b4=·(Di(hC))-K·hC

    式中:x、y、t分別為空間及時間坐標(biāo)系;h為水深;g為重力加速度;u和v分別為x和y向沿水深積分平均流速分量;C為污染物沿水深積分的垂線平均濃度;q為守恒物理量;f(q)為x向通量;g(q)為y向通量;b(q)為源(匯)項;S0x和Sfx分別為x向的水底底坡和摩阻坡度;S0y和Sfy分別為y向的水底底坡和摩阻坡度;Di為水平不同方向上的混合系數(shù);K為降解系數(shù);·為Laplace算子。

    1.2 邊界條件

    1.2.1 入流邊界

    對于入流邊界Γ0,須給定水位、流速及污染物濃度隨時間的變化值:

    1.2.2 出流邊界

    采用自由出流邊界,即:

    1.2.3 固壁邊界

    在湖岸等固壁邊界采用不可入條件,即vη=0(η為邊界法線方向)。

    1.3 基本方程的有限體積法離散

    對于任意單元Ω,其邊界為?Ω,通過對方程(1)積分并利用散度定理可得到有限體積法的基本方程:

    (2)

    式中:n為單元邊界?Ω的外法向單位向量;dω和dL分別為面積分和線積分微元;F(q)·n為法向數(shù)值通量,F(q)=(f(q),g(q))T。

    采用一階精度離散和f、q具有的旋轉(zhuǎn)不變性,得到式(2)的等價公式:

    (3)

    其中b*(q)=(Ab1,Ab2,Ab3,∑Di(hC)L-AKhC))

    式中:T(Φ)為坐標(biāo)軸旋轉(zhuǎn)角度Φ的變換矩陣;m為單元邊總數(shù);Lj為單元邊j的長度;A為單元面積。

    2 模型的檢驗

    2.1 網(wǎng)格劃分

    為了更好地擬合復(fù)雜的湖泊邊界,采用四邊形無結(jié)構(gòu)網(wǎng)格,網(wǎng)格尺寸變幅范圍約10 m×10 m~40 m×40 m,共有網(wǎng)格單元8 000多個,節(jié)點9 000余個。計算網(wǎng)格如圖2。

    圖2 計算網(wǎng)格

    2.2 模型驗證

    采用2003年4月資料在實際入流量和平均風(fēng)速條件下進(jìn)行模擬計算,糙率系數(shù)0.018~0.022。模擬流場如圖3,在各監(jiān)測點位上的流速在0.015~0.045 m/s之間,計算值與監(jiān)測結(jié)果對比見圖4,與監(jiān)測結(jié)果基本吻合,平均誤差為15.2%。

    圖3 驗證流場

    圖4 計算流速與實測流速對比

    3 暴雨徑流影響模擬分析

    3.1 設(shè)計暴雨徑流

    參照2004年入湖河流的實測水文資料,設(shè)計典型暴雨徑流過程如圖5。

    圖5 設(shè)計暴雨徑流過程線

    3.2 計算條件

    計算的水質(zhì)指標(biāo)為CODMn、NH3-N、TP和TN。五里湖水質(zhì)取汛期平均質(zhì)量濃度,入湖河流罵蠡港取汛期平均質(zhì)量濃度,各水質(zhì)指標(biāo)設(shè)計質(zhì)量濃度見表1。風(fēng)速取暴雨期間風(fēng)速12 m/s,風(fēng)向取主導(dǎo)風(fēng)向ESE。

    表1 各水質(zhì)指標(biāo)的設(shè)計質(zhì)量濃度 mg/L

    3.3 計算結(jié)果分析

    五里湖暴雨期模擬典型流場見圖6。

    圖6 五里湖暴雨期模擬典型流場

    3.3.1 CODMn擴散過程

    圖7 五里湖暴雨1~7 d后CODMn質(zhì)量濃度變化等值線(單位:mg/L)

    在暴雨發(fā)生后1~7 d內(nèi)五里湖CODMn濃度變化等值線見圖7。進(jìn)水1 d后,進(jìn)入的污水只局限在入水口附近的較小范圍內(nèi)(圖7(a))。進(jìn)水3 d后,污水的前鋒越過寶界橋,西五里湖開始受到影響(圖7(c))。進(jìn)水4 d后,污染物質(zhì)大規(guī)模穿過寶界橋,進(jìn)入西五里湖,部分湖區(qū)水質(zhì)由Ⅲ類變?yōu)棰纛?。進(jìn)水5 d后,五里湖水質(zhì)濃度分布趨于穩(wěn)定。西五里湖CODMn的質(zhì)量濃度有所升高,寶界橋附近的質(zhì)量濃度升高了約0.8 mg/L(14%)

    寶界橋斷面是東五里湖與西五里湖的交界斷面,CODMn動態(tài)交換過程情況如表2。由表2可見,隨著污染鋒面的通過,寶界橋斷面上的濃度分布趨向均勻,第7天的凈進(jìn)入通量僅0.018 t/d。暴雨期間的最大凈進(jìn)入通量達(dá)到1.539 t/d,大大高于西五里湖的自凈能力。

    表2 寶界橋斷面CODMn動態(tài)交換過程情況

    表3 寶界橋斷面TP、NH3-N、TN動態(tài)交換過程情況

    暴雨徑流對湖泊水質(zhì)的影響有明顯的區(qū)域性,由圖7可見,隨著時間的推移,污染物濃度分布的范圍開始發(fā)生變化。暴雨徑流對湖泊水質(zhì)產(chǎn)生一定影響,暴雨期間,罵蠡港進(jìn)水口附近出現(xiàn)明顯污染擴散帶,在風(fēng)生環(huán)流與擴散的共同作用下向湖心區(qū)遷移。 隨著污染物濃度的降低和湖泊水體的稀釋、污染物濃度的降解和擴散,河口附近水質(zhì)開始恢復(fù)。

    3.3.2 TP、NH3-N、TN擴散過程

    NH3-N、TP、TN的質(zhì)量濃度等值線圖變化特性與CODMn相似,這里不再給出。寶界橋斷面TP、NH3-N、TN的動態(tài)交換過程情況見表3。

    由表3可見,自罵蠡港等河流進(jìn)來的暴雨徑流對西五里湖水質(zhì)影響很大,污染物凈通量在第4天影響達(dá)到峰值,此后迅速減小。

    4 結(jié) 論

    a.構(gòu)建了基于四邊形網(wǎng)格的平面二維通量向量分裂(FVS)格式水量水質(zhì)耦合模型,較好地擬合五里湖曲折不規(guī)則邊界。

    b.暴雨徑流影響模擬結(jié)果表明:暴雨發(fā)生1 d后,影響只局限在入水口附近的較小范圍內(nèi),3 d后污水的前鋒越過寶界橋,西五里湖開始受到影響,4 d后污染物質(zhì)大規(guī)模穿過寶界橋,進(jìn)入西五里湖,部分湖區(qū)水質(zhì)由Ⅲ類變?yōu)棰纛?5 d后五里湖水質(zhì)濃度分布趨于穩(wěn)定。

    c.典型暴雨條件下,寶界橋斷面污染物通量在第四天達(dá)到峰值,此后迅速減小。

    [1] 柏祥,陳開寧,黃蔚,等.五里湖水質(zhì)現(xiàn)狀與變化趨勢[J].水資源保護(hù),2010,26(5):6-10.(BAI Xiang, CHEN Kaining,HUANG Wei, et al.Current status and variation tendency of water quality in Wuli Lake, Jiangsu Province[J].Water Resources Protection,2010,26(5):6-10.(in Chinese))

    [2] 秦伯強,胡維平,陳偉民,等.太湖水環(huán)境演化過程與機理[M].北京:科學(xué)出版社,2004.

    [3] 田向榮,馬巍,廖文根,等.調(diào)水對梅梁湖、五里湖水環(huán)境影響研究[J].人民長江,2007,38(2): 69-72.(TIAN Xiangrong,MA Wei,LIAO Wengen, et al.Effect of water diversion on water environment of Meilianghu lake and Wulihu lake [J].Yangtze River,2007,38(2): 69-72.(in Chinese))

    [4] 丁玲,逄勇,趙棣華,等.調(diào)水工程對五里湖水環(huán)境影響分析[J].河海大學(xué)學(xué)報:自然科學(xué)版, 2003, 31(4): 366-369.(DING Ling,PANG Yong,ZHAO Dihua, et al.Effect of water diversion project on water environment of Wulihu Lake[J].Journal of Hohai University:Natural Sciences,2003, 31(4): 366-369.(in Chinese))

    [5] 駱輝煌.RMA2 和RMA4 在五里湖調(diào)水模擬中的開發(fā)應(yīng)用研究[D].北京: 中國水利水電科學(xué)研究院, 2005.

    [6] 趙棣華,姚琪,蔣艷,等.通量向量分裂格式的二維水流水質(zhì)模擬[J].水科學(xué)進(jìn)展,2002,13(6):701-706.(ZHAO Dihua,YAO Qi,JIANG Yan, et al.FVS scheme in2-D depth-averaged flow pollutants modeling[J].Advances in Water Science,2002,13(6):701-706.(in Chinese))

    [7] 趙棣華,戚晨,庾維德,等.平面二維水流-水質(zhì)有限體積法及黎曼近似解模型[J].水科學(xué)進(jìn)展,2000,11(4):368-373.(ZHAO Dihua,QI Chen,YU Weide, et al.Finite volume method and riemann solver for depth-averaged two-dimensional flow-pollutants coupled model[J].Advances in Water Science, 2000,11(4):368-373.(in Chinese))

    Impact of rainstorm runoff on water quality of Wuli Lake

    LIU Xiaodong1,2,HU Gongyu3,ZHOU Yuanyuan2,TU Qile2,ZHENG Xiaoyu2

    (1.KeyLaboratoryofIntegratedRegulationandResourceDevelopmentonShallowLakes,MinistryofEducation,Nanjing210098,China;2.CollegeofEnvironment,HohaiUniversity,Nanjing210098,China;3.JingjiangHydrologyandWaterResourcesSurveyBureau,ChangjiangWaterResourcesCommission,Jingzhou434002,China)

    To fit the irregular boundaries of the Wuli Lake, quadrilateral mesh fitting has been used and a two-dimensional flux vector splitting (FVS) scheme flow-pollutant coupled model has been established, which was verificated after then.By using this model, the dynamic process of the impact of designed storm runoff on the water quality of Wuli Lake was simulated, the CODMn, NH3-N, TP and TN exchange fluxes at Baojie bridge section were analyzed.The result showed that in the designed typical storm conditions, the Maligang River got large flow of water which had a great impact on the water quality of West Wuli Lake.After 3-day designed storm occurred, runoff sewage began to affect the West Wuli Lake, and the concentration reached the peak after the fourth day, which changed the water quality from class Ⅲ to class Ⅳin some lake district, but the duration is short, which then decreases rapidly.

    storm runoff; Wuli Lake; two-dimensional flow-pollutant coupled model; quadrilateral mesh; CODMn

    10.3880/j.issn.1004-6933.2015.04.014

    國家自然科學(xué)基金(51479064,51379060)

    劉曉東(1972—),男,副教授,博士,主要從事環(huán)境與生態(tài)水力學(xué)、環(huán)境模擬等研究。E-mail:lxd54321@163.com

    涂琦樂。E-mail:15295527223@163.com

    X52

    A

    1004-6933(2015)04-0077-05

    2014-12-10 編輯:徐 娟)

    猜你喜歡
    調(diào)水通量暴雨
    “80年未遇暴雨”襲首爾
    冬小麥田N2O通量研究
    暴雨
    當(dāng)暴雨突臨
    調(diào)水專家傾囊相授15年心得,解答底改重要10問
    暴雨襲擊
    支點(2017年8期)2017-08-22 17:18:27
    緩釋型固體二氧化氯的制備及其釋放通量的影響因素
    重大報道中的“微表達(dá)”——以湖北日報《瓊瑤話調(diào)水》專欄為例
    新聞前哨(2015年2期)2015-03-11 19:29:22
    膜下滴灌作物以肥調(diào)水作用研究
    春、夏季長江口及鄰近海域溶解甲烷的分布與釋放通量
    精品一区二区三区av网在线观看| 蜜桃亚洲精品一区二区三区| 亚洲人与动物交配视频| 久久久a久久爽久久v久久| 少妇熟女aⅴ在线视频| 一区二区三区高清视频在线| 男人舔女人下体高潮全视频| 观看美女的网站| 一级黄片播放器| 国产探花在线观看一区二区| 成年女人看的毛片在线观看| 亚洲一级一片aⅴ在线观看| 久久6这里有精品| 夜夜夜夜夜久久久久| 国产在线男女| 久久久午夜欧美精品| 天堂√8在线中文| 美女cb高潮喷水在线观看| 色av中文字幕| 白带黄色成豆腐渣| 草草在线视频免费看| 嫩草影院新地址| 久久久久性生活片| 亚洲,欧美,日韩| 亚洲精品日韩在线中文字幕 | 99riav亚洲国产免费| 亚洲av电影不卡..在线观看| 婷婷精品国产亚洲av| 少妇被粗大猛烈的视频| 一个人免费在线观看电影| 非洲黑人性xxxx精品又粗又长| 久久人人爽人人片av| 国产视频内射| 久久国内精品自在自线图片| 成年女人永久免费观看视频| 亚洲成人中文字幕在线播放| 99精品在免费线老司机午夜| 在线播放国产精品三级| 在线免费观看不下载黄p国产| 卡戴珊不雅视频在线播放| 91在线精品国自产拍蜜月| 嫩草影院入口| 一本久久中文字幕| 99久国产av精品| 日韩成人伦理影院| 日韩中字成人| 久久久久国内视频| 成人综合一区亚洲| 波多野结衣高清无吗| 久久久久久久午夜电影| 国产人妻一区二区三区在| 欧美另类亚洲清纯唯美| 亚洲精品456在线播放app| 熟女人妻精品中文字幕| 99在线视频只有这里精品首页| 一级毛片我不卡| 午夜福利在线观看吧| 亚洲精品国产av成人精品 | 亚洲色图av天堂| 寂寞人妻少妇视频99o| av在线观看视频网站免费| 狂野欧美激情性xxxx在线观看| 国产片特级美女逼逼视频| 国产成人a区在线观看| 亚洲人成网站在线播放欧美日韩| 亚洲av免费在线观看| 俺也久久电影网| 午夜久久久久精精品| 精品国产一区二区三区久久久樱花| 精品久久久久久电影网| 午夜福利影视在线免费观看| av黄色大香蕉| 五月开心婷婷网| 国产精品无大码| 亚洲情色 制服丝袜| 韩国av在线不卡| 欧美老熟妇乱子伦牲交| 在线观看免费日韩欧美大片 | 久久人人爽人人片av| 国产精品人妻久久久久久| 国产黄色免费在线视频| 日本爱情动作片www.在线观看| 丝袜在线中文字幕| 国产精品一区二区在线观看99| 久久久久精品久久久久真实原创| 大片免费播放器 马上看| 18禁动态无遮挡网站| 国产av一区二区精品久久| 乱人伦中国视频| 99久久精品国产国产毛片| 成人国产av品久久久| kizo精华| 一个人免费看片子| 亚洲国产色片| 久久99热这里只频精品6学生| 91精品伊人久久大香线蕉| 夜夜骑夜夜射夜夜干| 国产 精品1| 国产亚洲一区二区精品| 少妇精品久久久久久久| 国产探花极品一区二区| 国产成人精品久久久久久| 亚洲精品自拍成人| 少妇人妻久久综合中文| 亚洲av电影在线观看一区二区三区| 在现免费观看毛片| 国产午夜精品一二区理论片| 少妇被粗大猛烈的视频| 欧美 亚洲 国产 日韩一| 国产精品麻豆人妻色哟哟久久| 欧美三级亚洲精品| 亚洲av电影在线观看一区二区三区| 国产亚洲一区二区精品| av福利片在线观看| 国产成人精品婷婷| 国产日韩欧美视频二区| 色94色欧美一区二区| 极品人妻少妇av视频| 日韩欧美精品免费久久| av视频免费观看在线观看| 亚洲精品视频女| 搡老乐熟女国产| a级毛色黄片| 麻豆成人av视频| 国产一区二区三区av在线| 久久精品久久精品一区二区三区| 免费黄频网站在线观看国产| 亚洲av国产av综合av卡| 免费观看性生交大片5| 中文字幕亚洲精品专区| 亚洲欧美精品自产自拍| 一级毛片久久久久久久久女| 99久久人妻综合| 一级毛片黄色毛片免费观看视频| 久久ye,这里只有精品| 在线观看www视频免费| 精品人妻熟女毛片av久久网站| 精品久久久久久久久亚洲| 卡戴珊不雅视频在线播放| 国产乱人偷精品视频| 一级毛片黄色毛片免费观看视频| 自拍偷自拍亚洲精品老妇| 自拍偷自拍亚洲精品老妇| 五月伊人婷婷丁香| 99热网站在线观看| 国产男女内射视频| 亚洲一级一片aⅴ在线观看| 国产精品人妻久久久影院| 97在线视频观看| 久久午夜福利片| 下体分泌物呈黄色| 久久精品国产鲁丝片午夜精品| 一个人看视频在线观看www免费| 观看av在线不卡| 久久午夜福利片| 精品一区二区三卡| 国产精品国产三级专区第一集| 嫩草影院新地址| 高清在线视频一区二区三区| 男女边摸边吃奶| 日本爱情动作片www.在线观看| 久久精品国产亚洲av天美| 高清午夜精品一区二区三区| 久久午夜综合久久蜜桃| 丝瓜视频免费看黄片| 国产av码专区亚洲av| 日韩三级伦理在线观看| 国精品久久久久久国模美| 国产免费一区二区三区四区乱码| 大片免费播放器 马上看| 精品久久久精品久久久| 下体分泌物呈黄色| 久久精品国产亚洲av天美| 久久av网站| 在线看a的网站| www.色视频.com| 国产成人免费无遮挡视频| 国产精品99久久99久久久不卡 | 免费观看性生交大片5| 欧美日韩av久久| av在线app专区| 十八禁网站网址无遮挡 | 精品一区二区三区视频在线| 国产成人精品一,二区| 国产精品一区二区三区四区免费观看| 欧美日韩视频精品一区| 免费大片18禁| 国产69精品久久久久777片| 国产精品女同一区二区软件| 一本久久精品| 91精品一卡2卡3卡4卡| 最新中文字幕久久久久| 久久久久精品性色| 老司机影院毛片| 日本猛色少妇xxxxx猛交久久| 热re99久久精品国产66热6| 色网站视频免费| 看非洲黑人一级黄片| 夫妻午夜视频| 日韩亚洲欧美综合| 如何舔出高潮| 久热这里只有精品99| 免费看光身美女| .国产精品久久| 亚洲欧美精品专区久久| 亚洲美女视频黄频| 美女xxoo啪啪120秒动态图| 97精品久久久久久久久久精品| 精品久久久噜噜| 免费观看av网站的网址| 日韩大片免费观看网站| 久久99热6这里只有精品| 国产高清不卡午夜福利| 国产黄片视频在线免费观看| 日韩欧美一区视频在线观看 | av黄色大香蕉| 国产视频首页在线观看| 中文在线观看免费www的网站| 国产免费一级a男人的天堂| 免费观看性生交大片5| 国产色婷婷99| 国产伦精品一区二区三区视频9| 亚洲欧美中文字幕日韩二区| 久久久国产精品麻豆| 欧美国产精品一级二级三级 | 久久久国产一区二区| 如何舔出高潮| 22中文网久久字幕| 久久午夜福利片| 国产成人午夜福利电影在线观看| xxx大片免费视频| 男人狂女人下面高潮的视频| 日韩一本色道免费dvd| 亚洲婷婷狠狠爱综合网| 亚洲久久久国产精品| 26uuu在线亚洲综合色| 久久97久久精品| av黄色大香蕉| 男女免费视频国产| 一级片'在线观看视频| 一区二区三区乱码不卡18| 极品教师在线视频| 在线精品无人区一区二区三| 国产精品成人在线| 久久久精品94久久精品| 午夜激情福利司机影院| 国产毛片在线视频| av免费观看日本| 国产高清国产精品国产三级| 内射极品少妇av片p| 亚洲国产av新网站| 人人澡人人妻人| 欧美丝袜亚洲另类| 精品人妻熟女毛片av久久网站| 又大又黄又爽视频免费| 卡戴珊不雅视频在线播放| 免费黄频网站在线观看国产| 久久免费观看电影| 777米奇影视久久| 久久人妻熟女aⅴ| 国产午夜精品久久久久久一区二区三区| 日本猛色少妇xxxxx猛交久久| 亚洲欧美中文字幕日韩二区| 国产精品麻豆人妻色哟哟久久| 交换朋友夫妻互换小说| 日韩av免费高清视频| 国产淫片久久久久久久久| 欧美日韩av久久| 一级片'在线观看视频| videossex国产| 极品教师在线视频| 欧美日本中文国产一区发布| 这个男人来自地球电影免费观看 | 国产亚洲精品久久久com| 99热国产这里只有精品6| 免费不卡的大黄色大毛片视频在线观看| 在线观看免费高清a一片| 熟妇人妻不卡中文字幕| 不卡视频在线观看欧美| 又大又黄又爽视频免费| 最近手机中文字幕大全| 国产精品麻豆人妻色哟哟久久| 一区二区三区精品91| 男人爽女人下面视频在线观看| 国产欧美亚洲国产| 少妇裸体淫交视频免费看高清| h视频一区二区三区| 精品一区二区免费观看| 久久久久久久国产电影| 最新的欧美精品一区二区| av卡一久久| 久久精品国产亚洲网站| 免费av中文字幕在线| 国产成人免费无遮挡视频| 国产乱来视频区| 18+在线观看网站| 亚洲av国产av综合av卡| 黑人高潮一二区| 三上悠亚av全集在线观看 | 99视频精品全部免费 在线| 人人澡人人妻人| h视频一区二区三区| 少妇人妻久久综合中文| 香蕉精品网在线| 国产精品久久久久久久电影| 亚洲精品国产色婷婷电影| 国产成人a∨麻豆精品| 国产精品女同一区二区软件| 国内精品宾馆在线| 国产欧美日韩精品一区二区| 一级毛片aaaaaa免费看小| 免费黄频网站在线观看国产| 亚洲欧美一区二区三区黑人 | 男女边吃奶边做爰视频| 国产乱来视频区| 亚洲欧美日韩东京热| 亚洲综合色惰| 精品人妻熟女av久视频| 免费大片黄手机在线观看| 亚州av有码| 性色av一级| 国产男女超爽视频在线观看| 亚洲欧洲日产国产| 国产精品久久久久久久久免| 九九在线视频观看精品| 久久ye,这里只有精品| 熟妇人妻不卡中文字幕| 永久免费av网站大全| 国产综合精华液| 老女人水多毛片| 久久热精品热| 久久97久久精品| 中文字幕久久专区| 日本午夜av视频| 99久久精品国产国产毛片| 亚洲av男天堂| 97精品久久久久久久久久精品| 国产在线免费精品| 国产亚洲欧美精品永久| 美女视频免费永久观看网站| 免费观看在线日韩| 少妇裸体淫交视频免费看高清| 91成人精品电影| 只有这里有精品99| 午夜91福利影院| 91精品一卡2卡3卡4卡| 久久久久久久久大av| 51国产日韩欧美| 伦理电影大哥的女人| 国产成人精品婷婷| 特大巨黑吊av在线直播| 午夜福利在线观看免费完整高清在| av天堂久久9| av在线app专区| 只有这里有精品99| 中文字幕精品免费在线观看视频 | 乱系列少妇在线播放| 大香蕉久久网| 丰满乱子伦码专区| 国产亚洲av片在线观看秒播厂| 日韩人妻高清精品专区| 人人澡人人妻人| 91成人精品电影| 国产精品国产av在线观看| 两个人免费观看高清视频 | 一级片'在线观看视频| 最近2019中文字幕mv第一页| 麻豆成人av视频| 深夜a级毛片| 毛片一级片免费看久久久久| 少妇精品久久久久久久| 91久久精品电影网| 最近中文字幕高清免费大全6| 老女人水多毛片| 街头女战士在线观看网站| 一区二区三区免费毛片| 久热这里只有精品99| 草草在线视频免费看| 99热这里只有是精品在线观看| 国产黄色免费在线视频| 热re99久久精品国产66热6| 亚洲美女黄色视频免费看| 免费看av在线观看网站| 99热国产这里只有精品6| 国产欧美日韩一区二区三区在线 | 中文字幕av电影在线播放| 国产毛片在线视频| 国产精品人妻久久久影院| 亚洲精品国产av蜜桃| 少妇的逼好多水| 免费在线观看成人毛片| 寂寞人妻少妇视频99o| 美女福利国产在线| 免费少妇av软件| 两个人的视频大全免费| 一级毛片久久久久久久久女| 丝袜在线中文字幕| 欧美精品一区二区大全| 成人18禁高潮啪啪吃奶动态图 | 久久av网站| 成人国产麻豆网| 亚洲人与动物交配视频| 亚洲电影在线观看av| 亚洲精品国产成人久久av| 成年av动漫网址| 亚洲av日韩在线播放| 91在线精品国自产拍蜜月| 少妇的逼水好多| 一本色道久久久久久精品综合| 国产69精品久久久久777片| 永久免费av网站大全| 99九九在线精品视频 | 成人特级av手机在线观看| 亚洲av在线观看美女高潮| 22中文网久久字幕| 日韩欧美精品免费久久| 另类亚洲欧美激情| 国产一区二区三区av在线| 五月天丁香电影| 成人免费观看视频高清| 久久人人爽av亚洲精品天堂| 国产一区有黄有色的免费视频| 国产极品粉嫩免费观看在线 | 老司机影院成人| 永久免费av网站大全| 狂野欧美白嫩少妇大欣赏| 亚洲精品久久久久久婷婷小说| 男人狂女人下面高潮的视频| 精品久久久久久电影网| 男人添女人高潮全过程视频| 国产日韩一区二区三区精品不卡 | 最近中文字幕2019免费版| 99久久精品国产国产毛片| 亚洲av欧美aⅴ国产| 99热6这里只有精品| 只有这里有精品99| 午夜精品国产一区二区电影| 97超视频在线观看视频| 国产精品国产av在线观看| 午夜免费鲁丝| 波野结衣二区三区在线| 99久国产av精品国产电影| 22中文网久久字幕| 亚洲精品456在线播放app| 精品国产一区二区久久| 精品一区二区三卡| 偷拍熟女少妇极品色| 少妇高潮的动态图| 一本一本综合久久| 久久99精品国语久久久| 丝瓜视频免费看黄片| 日韩中文字幕视频在线看片| 亚洲精品一二三| 日本爱情动作片www.在线观看| .国产精品久久| 男人和女人高潮做爰伦理| 亚洲精品日韩在线中文字幕| 国产欧美日韩一区二区三区在线 | 一级毛片我不卡| 午夜福利网站1000一区二区三区| 国产亚洲一区二区精品| 国产伦精品一区二区三区视频9| 精品亚洲乱码少妇综合久久| 国产 一区精品| 日韩不卡一区二区三区视频在线| 一级毛片黄色毛片免费观看视频| 深夜a级毛片| 日韩成人av中文字幕在线观看| 精品99又大又爽又粗少妇毛片| 亚洲成人av在线免费| 最近的中文字幕免费完整| 一区二区三区免费毛片| 极品少妇高潮喷水抽搐| 汤姆久久久久久久影院中文字幕| 一区二区三区四区激情视频| 国产精品伦人一区二区| 成人亚洲欧美一区二区av| 一区二区三区四区激情视频| 亚洲性久久影院| 亚洲精品国产成人久久av| 涩涩av久久男人的天堂| 日韩制服骚丝袜av| 国产乱来视频区| 久久免费观看电影| 国产老妇伦熟女老妇高清| 永久网站在线| 一级毛片我不卡| 欧美区成人在线视频| 亚洲成人手机| 熟女av电影| 99热网站在线观看| 亚洲欧美一区二区三区国产| 黄色一级大片看看| 成人漫画全彩无遮挡| 国产伦在线观看视频一区| 噜噜噜噜噜久久久久久91| 日本91视频免费播放| 国产熟女午夜一区二区三区 | 卡戴珊不雅视频在线播放| 国产精品人妻久久久久久| 桃花免费在线播放| 黑人巨大精品欧美一区二区蜜桃 | 人妻系列 视频| 伊人亚洲综合成人网| 在线 av 中文字幕| 国产免费福利视频在线观看| 国产精品久久久久成人av| 黑丝袜美女国产一区| videossex国产| 成人亚洲欧美一区二区av| av国产久精品久网站免费入址| 成人毛片a级毛片在线播放| 欧美丝袜亚洲另类| 亚洲一级一片aⅴ在线观看| kizo精华| 亚洲国产最新在线播放| freevideosex欧美| 老熟女久久久| 激情五月婷婷亚洲| h日本视频在线播放| 日本欧美国产在线视频| 深夜a级毛片| tube8黄色片| 欧美性感艳星| 全区人妻精品视频| 性色av一级| 成人国产av品久久久| 成人亚洲精品一区在线观看| 国产片特级美女逼逼视频| 深夜a级毛片| 国产高清三级在线| 少妇裸体淫交视频免费看高清| 九九在线视频观看精品| 国产精品一二三区在线看| 极品人妻少妇av视频| 亚洲人成网站在线观看播放| 精品少妇内射三级| 亚洲精品自拍成人| 高清在线视频一区二区三区| 男女国产视频网站| 一个人看视频在线观看www免费| 久久精品久久精品一区二区三区| 国产探花极品一区二区| 男女国产视频网站| 免费看不卡的av| 9色porny在线观看| 99国产精品免费福利视频| 丰满乱子伦码专区| 久久婷婷青草| 爱豆传媒免费全集在线观看| 又爽又黄a免费视频| 午夜日本视频在线| 中文乱码字字幕精品一区二区三区| 高清黄色对白视频在线免费看 | 99视频精品全部免费 在线| 最近手机中文字幕大全| 国产成人精品久久久久久| 曰老女人黄片| 久久热精品热| 久久99热6这里只有精品| 午夜av观看不卡| 嘟嘟电影网在线观看| 91aial.com中文字幕在线观看| 亚洲精品第二区| 亚洲av国产av综合av卡| a 毛片基地| 91久久精品电影网| 亚洲国产精品成人久久小说| 国产免费福利视频在线观看| 久久国内精品自在自线图片| 亚洲精品aⅴ在线观看| 最近2019中文字幕mv第一页| 亚洲精品aⅴ在线观看| 色婷婷久久久亚洲欧美| 欧美精品人与动牲交sv欧美| 日本黄色片子视频| 国产欧美日韩综合在线一区二区 | 午夜福利在线观看免费完整高清在| 黄色视频在线播放观看不卡| 伊人亚洲综合成人网| 国产成人aa在线观看| 男人舔奶头视频| 少妇熟女欧美另类| 69精品国产乱码久久久| 久久久久久人妻| 亚洲精品乱久久久久久| 秋霞伦理黄片| 人人澡人人妻人| 久久久久久久久久久免费av| 美女cb高潮喷水在线观看| 在线观看av片永久免费下载| 极品少妇高潮喷水抽搐| 综合色丁香网| 欧美精品人与动牲交sv欧美| 97在线人人人人妻| 黄色怎么调成土黄色| 久久久久久久久久久久大奶| 精品人妻一区二区三区麻豆| 色视频在线一区二区三区| 交换朋友夫妻互换小说| 亚洲精品一区蜜桃| 成人漫画全彩无遮挡| 久久99一区二区三区| 成年av动漫网址| 丰满迷人的少妇在线观看| 国产成人精品无人区| 久久久久久人妻| 成人漫画全彩无遮挡| 我的老师免费观看完整版| 欧美成人精品欧美一级黄| 日本猛色少妇xxxxx猛交久久| 久久久久久久久久人人人人人人| 亚洲欧美日韩另类电影网站| 女人精品久久久久毛片| 午夜免费男女啪啪视频观看|