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

    農(nóng)村河道不同生態(tài)治理模式數(shù)值模擬

    2022-08-27 07:34:14謝夢(mèng)瑤劉曉璇高玉琴
    江蘇水利 2022年8期
    關(guān)鍵詞:河道水位斷面

    胡 樂,謝夢(mèng)瑤,劉曉璇,吳 迪,高玉琴

    (1.江蘇省農(nóng)村水利科技發(fā)展中心,江蘇南京 210029;2.河海大學(xué)水利水電學(xué)院,江蘇南京 210098)

    農(nóng)村河道承擔(dān)著防洪、排澇、灌溉、引水、景觀娛樂等功能,與人們的生產(chǎn)和生活息息相關(guān),但由于農(nóng)業(yè)生產(chǎn)、養(yǎng)殖業(yè)和鄉(xiāng)鎮(zhèn)制造業(yè)的快速發(fā)展以及農(nóng)村居民環(huán)保意識(shí)的薄弱,農(nóng)村河道普遍出現(xiàn)泥沙淤積、水體污染、河岸坍塌等問題。

    國內(nèi)學(xué)者運(yùn)用水質(zhì)模型模擬河湖及水庫的水質(zhì)變化情況,取得了大量研究成果。郭鵬程等[1]利用MIKE21 模擬了生態(tài)湖的流動(dòng)及水體交換情況,對(duì)人工生態(tài)湖進(jìn)行優(yōu)化設(shè)計(jì);程桂[2]運(yùn)用SWMM 模型結(jié)合海綿城市水文過程模擬不同組合措施下區(qū)域的水質(zhì)變化;楊坤等[3]采用EFDC 模型定量預(yù)測烏東德水庫蓄水后的水質(zhì)演變過程;程揚(yáng)等[4]基于二維對(duì)流擴(kuò)散模型模擬生態(tài)清淤工程下的水庫水質(zhì)變化情況;翟敏婷等[5]基于QUAL2k 模型模擬了夏季和冬季的六種工況下礫石濾床對(duì)水質(zhì)的凈化效果;陳軍等[6]運(yùn)用MIKE21 中HD 模塊耦合Ecolab模塊,構(gòu)建了灞河流域水環(huán)境模型,分析灞河口的水質(zhì)特征。

    國外學(xué)者對(duì)水質(zhì)模型也進(jìn)行多方面研究,Panda等[7]采用人工神經(jīng)網(wǎng)絡(luò)技術(shù)與MIKEI1相結(jié)合來預(yù)測水位和水質(zhì)的變化;Allam A[8]通過整合QUAL2K模型和遺傳算法模擬埃及尼羅河三角洲的Gharbia 河夏季和冬季的污染負(fù)荷減排量;Guangyi Deng等[9]采用MIKE11一維水動(dòng)力水質(zhì)模型和物理生境模型模擬了水庫運(yùn)行方案實(shí)施后的河流水質(zhì)和物理生境。

    目前水質(zhì)模擬模型和模擬方法具有多樣性,在河流、湖泊和海洋等水質(zhì)模擬研究中取得了大量研究成果,但總體研究偏向于大中型河流和湖泊的水質(zhì)水量研究,對(duì)中小型河道的水質(zhì)模擬研究較少,尤其在針對(duì)中小河流生態(tài)治理措施對(duì)河道水質(zhì)的改善效果方面的研究少之又少,不同河道生態(tài)治理措施對(duì)河道水環(huán)境的影響還有待深入研究。

    1 數(shù)據(jù)與方法

    1.1 研究區(qū)概況

    江蘇省無錫市錫山區(qū)位于長江三角洲腹地,地貌類型為長江三角洲沖洪積平原。境內(nèi)總體上構(gòu)成了“四縱四橫”骨干水系主框架,這些河道均是本區(qū)防洪、排澇、引水灌溉、通航的主要通道。本研究選取錫山區(qū)境內(nèi)一條未治理的農(nóng)村河道A 作為研究對(duì)象,其地理位置圖如圖1 所示。該河道河口寬19.6 m,橫斷面為概化為梯形,平均流速為0.1 m/s,河道為南北走向,規(guī)劃整治長度約為889 m,排澇標(biāo)準(zhǔn)為20年一遇。

    圖1 研究區(qū)地理位置

    1.2 研究數(shù)據(jù)

    由于研究區(qū)河道規(guī)模小,流速緩,水位變幅小,未設(shè)有水文與水質(zhì)監(jiān)測站點(diǎn),因此采用2021 年每10 d一次的實(shí)測水位、流量數(shù)據(jù)以及2021年每月一次的COD、NH3-N、TP實(shí)測質(zhì)量濃度數(shù)據(jù)進(jìn)行研究。

    1.3 模型原理

    水動(dòng)力模塊(HD)采用一維非恒定流方程組(圣維南方程組)對(duì)水體非連續(xù)性運(yùn)動(dòng)規(guī)律進(jìn)行模擬分析,該方程組反映了質(zhì)量守恒和動(dòng)量守恒定律。模型的控制方程為

    連續(xù)方程:

    動(dòng)量方程:

    式中:Q為流量;x為沿水流方向的距離;A為過水?dāng)嗝婷娣e;t為時(shí)間坐標(biāo);q為旁側(cè)入流量;g為重力加速度;h為水位;R為水力半徑;C為謝才系數(shù);α為動(dòng)量修正系數(shù)。

    水流中物質(zhì)在各處因質(zhì)量濃度差存在擴(kuò)散運(yùn)動(dòng),且隨河道平均流量運(yùn)移,對(duì)流擴(kuò)散方程模擬水體中溶解或懸浮的保守和非保守物質(zhì)的運(yùn)輸過程和一級(jí)線型衰減過程。

    對(duì)流擴(kuò)散方程:

    式中:x為距離坐標(biāo),m;t為時(shí)間坐標(biāo),s;C為物質(zhì)質(zhì)量濃度,mg/L;A為過水?dāng)嗝婷娣e,m2;D為縱向擴(kuò)散系數(shù),m2/s;Q為流量,m3/s;K為衰減系數(shù),1/d;C2為污染物源匯質(zhì)量濃度,mg/L。

    2 模型構(gòu)建

    2.1 模型構(gòu)建

    水動(dòng)力模型構(gòu)建:農(nóng)村河道缺少斷面數(shù)據(jù),在滿足模擬精度的前提下,根據(jù)實(shí)測數(shù)據(jù)及規(guī)劃資料,將河道斷面概化為梯形斷面。上游邊界為時(shí)間——流量序列,下游邊界為時(shí)間——水位序列。河長較短,河道形態(tài)變化小,因此河道糙率取全域值。

    水質(zhì)模型構(gòu)建:外部邊界為研究河段上下邊界實(shí)測污染物質(zhì)量濃度,內(nèi)部邊界條件為距河道上邊界200 m 處的排污口和河道概化的底泥面源污染,參考《全國水環(huán)境容量核定技術(shù)指南》中排污系數(shù)取值,結(jié)合現(xiàn)場對(duì)排污口的水質(zhì)檢測結(jié)果,確定污染物入河系數(shù)取為0.9。定期取樣測定并計(jì)算出底泥污染物中COD、NH3-N、TP 的入河貢獻(xiàn)量分別為14.4%、20.6%和18.3%。

    2.2 參數(shù)率定與驗(yàn)證

    河長較短,河道形態(tài)變化小,因此河道糙率取全域值0.032,水位模擬結(jié)果與實(shí)測數(shù)據(jù)如圖2所示。

    采用Nash-Sutcliffe 系數(shù)對(duì)模型精度進(jìn)行評(píng)估,公式為

    通過參考國內(nèi)部分河流擴(kuò)散系數(shù)取值,初步確定河道擴(kuò)散系數(shù),通過模擬值與實(shí)測值的比較不斷調(diào)整擴(kuò)散系數(shù),最終確定研究河道的擴(kuò)散系數(shù)為8 m2/s。利用2021年1—12月的水質(zhì)數(shù)據(jù)進(jìn)行率定,反復(fù)調(diào)試參數(shù)值優(yōu)化模型,最終確定所有的參數(shù)取值見表1。

    表1 衰減系數(shù)參數(shù)取值

    相應(yīng)的水質(zhì)指標(biāo)模擬結(jié)果與實(shí)測數(shù)據(jù)如圖3所示。

    根據(jù)圖3,河道斷面水質(zhì)模擬值與實(shí)測值變化趨勢基本吻合,相對(duì)誤差在0.9%~1.33%之間,在合理誤差范圍內(nèi),模型可用于河道的水質(zhì)模擬。

    圖3 各水質(zhì)指標(biāo)模擬值與實(shí)測值對(duì)比

    3 生態(tài)治理模式擬定

    根據(jù)研究區(qū)河道生態(tài)治理前數(shù)值模擬結(jié)果分析,將對(duì)A 河進(jìn)行截污控制方案模擬作為基礎(chǔ)模式,運(yùn)用河道生態(tài)治理技術(shù)針對(duì)性地?cái)M定3 種生態(tài)治理模式:模式1(生態(tài)浮床+曝氣增氧);模式2(仿木樁護(hù)坡+自然護(hù)坡);模式3(水下森林)。選取錫山區(qū)采用以上生態(tài)治理模式的3 條農(nóng)村河道B 河、C 河和D 河,分別建立這3 條河道的水動(dòng)力水質(zhì)模型,通過率定得到3 種生態(tài)治理模式的衰減系數(shù)如表2所示。

    表2 三種生態(tài)治理模式下的衰減系數(shù)率定值

    4 模擬結(jié)果分析

    4.1 污染源截污控制方案對(duì)河道水質(zhì)的影響

    A河的主要污染物來源為生活污水和底泥內(nèi)源釋放,本文選取水質(zhì)較差的10 月份水質(zhì)數(shù)據(jù),分別模擬實(shí)施排污口截流和清淤疏浚后的水質(zhì)情況,與河道治理前水質(zhì)狀況進(jìn)行比較,如圖4所示。

    由圖4可知,對(duì)生活污水進(jìn)行截流處理后,河道水質(zhì)有所提升,但效果不顯著,COD 平均降低0.05%,NH3-N 降低0.12%,TP 降低0.06%。進(jìn)行河道清淤后,河道沿程污染物質(zhì)量濃度降低,對(duì)污染物質(zhì)量濃度分布曲線進(jìn)行擬合,計(jì)算平均每百米

    圖4 污染物截污控制方案治理效果

    污染物質(zhì)量濃度降低率,COD 降低0.48%,NH3-N降低1.42%,TP 降低0.15%,河道水質(zhì)得到有效改善。清淤疏浚的治理效果優(yōu)于排污口截流,說明該河道水質(zhì)主要受底泥內(nèi)源污染的影響,河道沿程污染物逐漸累積,至河道下邊界達(dá)到最大值,清除底泥污染源后河道水質(zhì)沿程逐漸改善;而排污口由于流量較小,且入河污染負(fù)荷不大,因此排污口截流后河段的污染物質(zhì)量濃度較治理前變化較小。

    4.2 不同生態(tài)治理模式下的水位變化情況

    4.2.1 天然流量下的水位變化情況

    流量過程采用2021 年1—12 月的實(shí)測天然來水流量,通過水位模擬與現(xiàn)狀河道斷面對(duì)比分析,為河道的斷面設(shè)計(jì)提供參考。A-400斷面在3種河道生態(tài)治理模式下的水位變化情況如圖5所示。

    圖5 河道水位變化情況

    不同的河道生態(tài)治理模式對(duì)河道糙率的影響較大,在河道天然流量過程下,河道水位變化較大。根據(jù)水位計(jì)算結(jié)果,3 種模式下的河道水位相比于河道治理基礎(chǔ)模式分別升高0.3 m、0.78 m、0.92 m。因此在進(jìn)行河道生態(tài)治理時(shí),應(yīng)充分考慮河道的過流能力,避免因河道治理導(dǎo)致河道漫堤,無法承擔(dān)防洪和灌溉等河道功能。

    4.2.2 排澇標(biāo)準(zhǔn)下的水位變化情況

    A 河為錫北鎮(zhèn)圩內(nèi)小型農(nóng)村河道,河道功能為排澇及景觀作用,因此僅計(jì)算河道在排澇標(biāo)準(zhǔn)下的水位變化情況。該地區(qū)的內(nèi)河排澇標(biāo)準(zhǔn)為20 年一遇,A河集水面積為0.8 km2,來水流量按排澇模數(shù)M按每平方千米4.01 m3/s計(jì)算。

    設(shè)計(jì)流量的計(jì)算公式為

    計(jì)算得出A河的設(shè)計(jì)排澇流量為3.21 m3/s。

    A河橫斷面為梯形,底坡較緩,可作為明渠均勻流進(jìn)行計(jì)算。A河河底高程為2.1 m,底寬6 m,兩岸邊坡1∶2,岸頂高程為5.53 m,設(shè)計(jì)安全超高為0.5 m,則最大過流水深為2.93 m。據(jù)曼寧公式(式(6))計(jì)算河道在不同生態(tài)治理模式下的排澇設(shè)計(jì)水位值,并復(fù)核河道過流能力,如表3所示。

    表3 不同生態(tài)治理模式下河道過流能力計(jì)算

    根據(jù)計(jì)算結(jié)果,河道在基礎(chǔ)模式和模式2 兩種生態(tài)治理模式下,排澇設(shè)計(jì)水位為4.42 m和4.79 m,低于河道允許最高過流水位5.03 m,河道最大過流量為5.12 m3/s 和3.81 m3/s,大于河道設(shè)計(jì)排澇流量3.21 m3/s,河道斷面滿足排澇要求;河道在模式1 和模式3兩種生態(tài)治理模式下,排澇設(shè)計(jì)水位為5.4 m和5.56 m,高于河道允許最高過流水位5.03 m,河道最大過流量為2.52 m3/s 和2.27 m3/s,小于河道設(shè)計(jì)排澇流量3.21 m3/s,河道斷面不滿足排澇要求。因此在進(jìn)行河道整治時(shí),需要結(jié)合現(xiàn)場調(diào)查的河道具體情況,適當(dāng)采取增大河寬、降低河底高程、增大河道底坡等方法增大河道斷面面積,提高河道過流能力。

    4.3 不同生態(tài)治理模式對(duì)河道水質(zhì)的影響

    將B、C、D 三條河道的衰減系數(shù)率定值代入A河,模擬A 河A-400 斷面在4 種不同治理模式下的COD、NH3-N、TP 的變化情況,模擬時(shí)間為2021 年1—12月,模擬結(jié)果見圖6。

    由圖6分析可得,3種治理模式均對(duì)A河污染物有一定的降解作用。從3種模式對(duì)污染物的年平均降解效果來看,對(duì)于COD 指標(biāo),模式3(0.55%)>模式1(0.41%)>模式2(0.19%);對(duì)于NH3-N指標(biāo),模式3(0.80%)>模式1(0.56%)>模式2(0.39%);對(duì)于TP指標(biāo),模式3(0.68%)>模式1(0.57%)>模式2(0.26%)。從總體降解效果上看,模式3>模式1>模式2,這是由于模式3水下森林系統(tǒng)構(gòu)建出的復(fù)合生態(tài)系統(tǒng),通過物理吸附、植物吸收和微生物分解等物理、生物和化學(xué)反應(yīng)共同作用,各反應(yīng)之間相互影響,相互促進(jìn),形成了降解污染物的高效循環(huán)系統(tǒng)。模式1 生態(tài)浮床+曝氣增氧系統(tǒng)作用原理與模式3水下森林類似,模式1浮床面積較小,構(gòu)建的生態(tài)系統(tǒng)中動(dòng)植物和微生物數(shù)量遠(yuǎn)小于模式3,但增設(shè)的曝氣增氧設(shè)施增大了水體中的含氧量,提高了水生動(dòng)植物和微生物的活性,因此對(duì)河道污染物的降解效果僅次于模式1。

    圖6 不同治理模式下河道污染物濃度隨時(shí)間的相對(duì)變化率

    從時(shí)間上看,基礎(chǔ)模式下A 河河道污染物質(zhì)量濃度總體表現(xiàn)為夏季污染物質(zhì)量濃度較低,春秋冬三季污染物質(zhì)量濃度較高。A河水質(zhì)主要受上游來水影響,夏季水量相對(duì)充沛,動(dòng)植物活性較高,水質(zhì)相對(duì)較好;而春秋冬三季水量相對(duì)較小,春秋季節(jié)化肥農(nóng)藥使用量較大,殘留的污染物通過渠道排入河中,加上秋冬季節(jié)植物衰敗腐爛,混合在河道底泥中形成內(nèi)源污染,因此河道水質(zhì)較差。從3 種模式對(duì)污染物降解效果的時(shí)間變化規(guī)律上看,不同時(shí)期的污染物質(zhì)量濃度變化較大,而污染物降解速率變化較小,不同污染物質(zhì)量濃度下,降解速率相差在3%以內(nèi),說明當(dāng)污染物質(zhì)量濃度在一定范圍內(nèi)時(shí),污染物的降解速率不受污染物質(zhì)量濃度影響。

    5 結(jié) 論

    本文模擬了錫山區(qū)未治理的農(nóng)村河道A 河在排污口截流與清淤疏浚2 種措施下的水質(zhì)改善效果,并在此基礎(chǔ)上,分別模擬在3種不同生態(tài)治理模式下的水位、水質(zhì)變化情況。主要結(jié)論如下:

    (1)該河道水質(zhì)主要受底泥內(nèi)源污染的影響,入河污染負(fù)荷不大,清除底泥污染源后河道水質(zhì)沿程逐漸改善。

    (2)在河道排澇標(biāo)準(zhǔn)下,基礎(chǔ)模式和模式2兩種生態(tài)治理模式下的排澇設(shè)計(jì)水位低于河道允許最高過流水位,模式1 和模式3 兩種生態(tài)治理模式下的排澇設(shè)計(jì)水位高于河道允許最高過流水位,需要對(duì)河道斷面進(jìn)行擴(kuò)建。

    (3)3 種治理模式均對(duì)A 河污染物有一定的降解作用,從總體降解效果上看,模式3>模式1>模式2。對(duì)于同一污染指標(biāo),模式3 的降解速率最大,對(duì)于同一模式,NH3-N的降解速率最大。

    猜你喜歡
    河道水位斷面
    ??? ??? ?? ????? ? ??(斷面)
    河道里的垃圾
    山溪性河道植物選擇與應(yīng)用
    南京市鼓樓區(qū)黑臭河道的治理
    某河河道演變分析
    基于MFAC-PID的核電站蒸汽發(fā)生器水位控制
    基于Dijkstra算法的最優(yōu)解列斷面快速搜索方法
    國內(nèi)首臺(tái)全斷面煤巷高效掘進(jìn)機(jī)誕生并交付
    天鋼Ф210mm斷面圓坯的開發(fā)與生產(chǎn)實(shí)踐
    天津冶金(2014年4期)2014-02-28 16:52:32
    基于PLC的水位控制系統(tǒng)的設(shè)計(jì)與研究
    河南科技(2014年4期)2014-02-27 14:07:11
    宅男免费午夜| 一本久久精品| 国产精品久久久人人做人人爽| 午夜精品国产一区二区电影| 精品少妇一区二区三区视频日本电影 | 啦啦啦在线免费观看视频4| 午夜福利网站1000一区二区三区| 日韩熟女老妇一区二区性免费视频| 国产一卡二卡三卡精品 | 秋霞伦理黄片| 国产又爽黄色视频| 国产无遮挡羞羞视频在线观看| 欧美人与性动交α欧美软件| 99九九在线精品视频| 国产视频首页在线观看| 99热网站在线观看| 两性夫妻黄色片| 久久久久久久久久久免费av| 捣出白浆h1v1| 日韩欧美一区视频在线观看| 极品少妇高潮喷水抽搐| av线在线观看网站| 最近2019中文字幕mv第一页| 国产精品嫩草影院av在线观看| 亚洲欧美一区二区三区久久| 丁香六月欧美| 精品少妇久久久久久888优播| 久久久久久久精品精品| 黑人猛操日本美女一级片| 国产男女内射视频| e午夜精品久久久久久久| 亚洲国产av新网站| 中文字幕人妻熟女乱码| 久久免费观看电影| 永久免费av网站大全| 纵有疾风起免费观看全集完整版| 黑丝袜美女国产一区| 久久国产精品大桥未久av| 天天添夜夜摸| 女性生殖器流出的白浆| 日韩熟女老妇一区二区性免费视频| 成人免费观看视频高清| 成人亚洲精品一区在线观看| 精品国产乱码久久久久久男人| 久久久久精品人妻al黑| 19禁男女啪啪无遮挡网站| 久久鲁丝午夜福利片| 深夜精品福利| 成年人免费黄色播放视频| 国产伦理片在线播放av一区| 少妇被粗大的猛进出69影院| 国产淫语在线视频| 国产色婷婷99| 国产亚洲精品第一综合不卡| 涩涩av久久男人的天堂| 午夜精品国产一区二区电影| 国产爽快片一区二区三区| 三上悠亚av全集在线观看| 黄色 视频免费看| 久久久久久久大尺度免费视频| 青春草国产在线视频| 在线免费观看不下载黄p国产| 亚洲三区欧美一区| 亚洲三区欧美一区| 亚洲三区欧美一区| 婷婷色av中文字幕| 亚洲熟女毛片儿| 激情视频va一区二区三区| 久久久精品区二区三区| 国产精品免费视频内射| 亚洲三区欧美一区| 亚洲欧美激情在线| 如何舔出高潮| 晚上一个人看的免费电影| 成年av动漫网址| 高清在线视频一区二区三区| 国产在线一区二区三区精| 嫩草影视91久久| 日韩制服丝袜自拍偷拍| 18禁观看日本| 老熟女久久久| 欧美精品高潮呻吟av久久| 搡老岳熟女国产| 两个人免费观看高清视频| 亚洲欧美成人精品一区二区| 19禁男女啪啪无遮挡网站| 国产精品av久久久久免费| 一级毛片黄色毛片免费观看视频| 欧美精品人与动牲交sv欧美| 只有这里有精品99| 亚洲精品乱久久久久久| 亚洲国产精品一区三区| 丁香六月欧美| 黄色视频在线播放观看不卡| 久久久久国产精品人妻一区二区| 黄色 视频免费看| 国产亚洲av片在线观看秒播厂| 丁香六月天网| 青春草视频在线免费观看| 亚洲欧美一区二区三区国产| 国产男女内射视频| 亚洲情色 制服丝袜| 99久久99久久久精品蜜桃| 五月天丁香电影| 亚洲欧美中文字幕日韩二区| 国产成人欧美| 久久狼人影院| 亚洲成av片中文字幕在线观看| 精品酒店卫生间| 男女免费视频国产| 两性夫妻黄色片| 久久久精品国产亚洲av高清涩受| 夫妻午夜视频| av卡一久久| 亚洲精品国产av成人精品| 777久久人妻少妇嫩草av网站| 亚洲国产欧美日韩在线播放| videos熟女内射| 久久久欧美国产精品| 看十八女毛片水多多多| 一级片免费观看大全| 少妇人妻久久综合中文| 青春草亚洲视频在线观看| 亚洲,一卡二卡三卡| 五月天丁香电影| 另类精品久久| 日本91视频免费播放| 人妻 亚洲 视频| 亚洲av日韩精品久久久久久密 | 亚洲精品一区蜜桃| 日韩欧美一区视频在线观看| 一本一本久久a久久精品综合妖精| 老司机影院毛片| 免费高清在线观看视频在线观看| 操出白浆在线播放| 中国三级夫妇交换| 夫妻午夜视频| 在线观看免费高清a一片| 女性被躁到高潮视频| 欧美 日韩 精品 国产| 无遮挡黄片免费观看| 亚洲,欧美,日韩| 亚洲欧美一区二区三区黑人| 大香蕉久久成人网| 精品午夜福利在线看| 国产成人免费无遮挡视频| 一本大道久久a久久精品| 国产精品免费视频内射| 高清欧美精品videossex| 亚洲欧美精品综合一区二区三区| 国产精品蜜桃在线观看| 国产色婷婷99| 久久久国产精品麻豆| 国产精品偷伦视频观看了| 成人国产av品久久久| 大香蕉久久成人网| 视频区图区小说| 国产精品一区二区在线不卡| 亚洲欧美一区二区三区久久| 侵犯人妻中文字幕一二三四区| 亚洲精品国产av蜜桃| 男人添女人高潮全过程视频| 亚洲在久久综合| 欧美在线一区亚洲| 又大又黄又爽视频免费| 午夜久久久在线观看| 中文字幕av电影在线播放| 久久久久国产精品人妻一区二区| 三上悠亚av全集在线观看| 成年动漫av网址| 大香蕉久久网| 免费人妻精品一区二区三区视频| 少妇的丰满在线观看| 久久久久久久精品精品| 777久久人妻少妇嫩草av网站| 精品国产超薄肉色丝袜足j| 日本91视频免费播放| 19禁男女啪啪无遮挡网站| 美女大奶头黄色视频| 日韩伦理黄色片| 男女高潮啪啪啪动态图| 亚洲欧美清纯卡通| 亚洲国产欧美日韩在线播放| 少妇被粗大的猛进出69影院| 亚洲第一区二区三区不卡| 中文字幕人妻熟女乱码| av天堂久久9| 我要看黄色一级片免费的| 亚洲精品国产av蜜桃| 国产精品二区激情视频| 交换朋友夫妻互换小说| 99久久精品国产亚洲精品| 久久99热这里只频精品6学生| 老司机深夜福利视频在线观看 | 久久天躁狠狠躁夜夜2o2o | 国产麻豆69| 青青草视频在线视频观看| 亚洲国产精品成人久久小说| av线在线观看网站| 亚洲国产日韩一区二区| 国产极品天堂在线| 国产97色在线日韩免费| 天堂8中文在线网| 亚洲国产精品成人久久小说| 国产成人欧美| 精品一区二区三卡| 高清视频免费观看一区二区| 波多野结衣av一区二区av| 国产毛片在线视频| 一级片免费观看大全| 日韩一卡2卡3卡4卡2021年| 啦啦啦中文免费视频观看日本| 国产精品一区二区精品视频观看| 97在线人人人人妻| 热99久久久久精品小说推荐| 在线天堂中文资源库| 亚洲第一区二区三区不卡| 成人国产麻豆网| 青春草亚洲视频在线观看| 色播在线永久视频| 激情视频va一区二区三区| 日韩 欧美 亚洲 中文字幕| 成人亚洲精品一区在线观看| 黄片无遮挡物在线观看| 男女高潮啪啪啪动态图| www日本在线高清视频| 久久精品国产a三级三级三级| 亚洲视频免费观看视频| 国产熟女欧美一区二区| 69精品国产乱码久久久| 国产精品三级大全| 毛片一级片免费看久久久久| 老司机靠b影院| 国产精品久久久人人做人人爽| 欧美xxⅹ黑人| 性少妇av在线| 国产精品女同一区二区软件| 超碰成人久久| 少妇精品久久久久久久| 黄色一级大片看看| 18在线观看网站| 国产精品欧美亚洲77777| www日本在线高清视频| 精品亚洲乱码少妇综合久久| 19禁男女啪啪无遮挡网站| 国产精品成人在线| 久久久久精品国产欧美久久久 | 国产xxxxx性猛交| 国产免费现黄频在线看| 国产精品 国内视频| 欧美日韩亚洲国产一区二区在线观看 | 亚洲国产欧美网| 亚洲精品在线美女| 一区二区日韩欧美中文字幕| 欧美精品人与动牲交sv欧美| 久久久久精品国产欧美久久久 | 日韩电影二区| 在线天堂中文资源库| 国产精品一国产av| 夫妻午夜视频| 午夜激情av网站| 搡老岳熟女国产| 国产免费福利视频在线观看| 亚洲av福利一区| 老司机影院成人| 男女无遮挡免费网站观看| 交换朋友夫妻互换小说| 精品一区在线观看国产| 亚洲欧美日韩另类电影网站| 久久国产精品男人的天堂亚洲| 一级毛片黄色毛片免费观看视频| 国产一区亚洲一区在线观看| 亚洲精品日本国产第一区| 欧美97在线视频| 欧美精品一区二区大全| www.熟女人妻精品国产| 亚洲熟女毛片儿| 亚洲男人天堂网一区| 哪个播放器可以免费观看大片| 精品少妇久久久久久888优播| 亚洲专区中文字幕在线 | 午夜免费鲁丝| 蜜桃国产av成人99| 欧美日本中文国产一区发布| 国产精品一国产av| 久久国产亚洲av麻豆专区| 亚洲精品一区蜜桃| 精品国产乱码久久久久久男人| 一本色道久久久久久精品综合| 日本欧美国产在线视频| 国产午夜精品一二区理论片| 高清黄色对白视频在线免费看| 国产精品人妻久久久影院| 下体分泌物呈黄色| 亚洲视频免费观看视频| 成人手机av| 亚洲,一卡二卡三卡| 精品亚洲乱码少妇综合久久| 男女免费视频国产| 男女边吃奶边做爰视频| 国产乱来视频区| 国产伦理片在线播放av一区| av免费观看日本| 美女福利国产在线| 午夜福利视频在线观看免费| 免费观看av网站的网址| 久久久久精品国产欧美久久久 | 久久久久精品人妻al黑| 69精品国产乱码久久久| 嫩草影院入口| 亚洲成人免费av在线播放| av网站免费在线观看视频| 国产片内射在线| 久久久精品区二区三区| 18禁国产床啪视频网站| 美女主播在线视频| 夜夜骑夜夜射夜夜干| 日韩精品免费视频一区二区三区| 人成视频在线观看免费观看| 国产精品欧美亚洲77777| 亚洲精品国产区一区二| 精品人妻熟女毛片av久久网站| 亚洲伊人久久精品综合| 亚洲精品日本国产第一区| 亚洲人成网站在线观看播放| 婷婷色麻豆天堂久久| 亚洲精品自拍成人| 中文字幕色久视频| 午夜日本视频在线| 久久99热这里只频精品6学生| 久久久久久人人人人人| 免费高清在线观看日韩| 亚洲av福利一区| 激情视频va一区二区三区| videosex国产| 一级毛片 在线播放| 国产精品99久久99久久久不卡 | 一级片'在线观看视频| 亚洲三区欧美一区| 九草在线视频观看| av又黄又爽大尺度在线免费看| 最近2019中文字幕mv第一页| kizo精华| 久久精品久久久久久噜噜老黄| 国产免费一区二区三区四区乱码| tube8黄色片| 男的添女的下面高潮视频| av一本久久久久| 国产成人欧美| 亚洲av日韩精品久久久久久密 | 精品国产露脸久久av麻豆| 又粗又硬又长又爽又黄的视频| 亚洲国产精品成人久久小说| 1024香蕉在线观看| 91aial.com中文字幕在线观看| 国产成人免费无遮挡视频| 丝袜喷水一区| 国产极品天堂在线| 啦啦啦视频在线资源免费观看| 波多野结衣av一区二区av| 视频在线观看一区二区三区| 天美传媒精品一区二区| 国产精品蜜桃在线观看| 亚洲人成77777在线视频| 成年女人毛片免费观看观看9 | 成年人免费黄色播放视频| 99精品久久久久人妻精品| 校园人妻丝袜中文字幕| 精品亚洲成a人片在线观看| 久久天躁狠狠躁夜夜2o2o | 黄片播放在线免费| 国产深夜福利视频在线观看| 久久人妻熟女aⅴ| 男女边吃奶边做爰视频| 亚洲欧美日韩另类电影网站| 亚洲av日韩精品久久久久久密 | 天天操日日干夜夜撸| 黑人欧美特级aaaaaa片| av在线app专区| 国产黄频视频在线观看| 啦啦啦在线免费观看视频4| 欧美成人精品欧美一级黄| 国产男女超爽视频在线观看| 成人三级做爰电影| 水蜜桃什么品种好| 老汉色av国产亚洲站长工具| 女人精品久久久久毛片| 亚洲欧美中文字幕日韩二区| 亚洲国产日韩一区二区| 日韩中文字幕视频在线看片| 激情视频va一区二区三区| 90打野战视频偷拍视频| av网站在线播放免费| 热re99久久精品国产66热6| 最近中文字幕2019免费版| 中文精品一卡2卡3卡4更新| 亚洲精品第二区| 久久国产精品大桥未久av| 国产精品 欧美亚洲| 黄色 视频免费看| 色网站视频免费| 丰满乱子伦码专区| 久久久久久人妻| 亚洲激情五月婷婷啪啪| 久久97久久精品| 久久精品熟女亚洲av麻豆精品| 男人操女人黄网站| 精品久久久精品久久久| 午夜福利,免费看| 看非洲黑人一级黄片| 亚洲一区中文字幕在线| 中文字幕人妻丝袜一区二区 | 狠狠精品人妻久久久久久综合| 男女之事视频高清在线观看 | 街头女战士在线观看网站| 亚洲四区av| 尾随美女入室| 少妇人妻久久综合中文| 侵犯人妻中文字幕一二三四区| 国产av精品麻豆| av天堂久久9| 91成人精品电影| 中文字幕人妻熟女乱码| 成年人午夜在线观看视频| 90打野战视频偷拍视频| 欧美久久黑人一区二区| 搡老岳熟女国产| 又大又黄又爽视频免费| 九色亚洲精品在线播放| 国产成人精品在线电影| 黑人欧美特级aaaaaa片| 欧美av亚洲av综合av国产av | 国产精品蜜桃在线观看| av网站免费在线观看视频| 嫩草影院入口| 高清视频免费观看一区二区| 亚洲综合色网址| 午夜精品国产一区二区电影| 日韩av在线免费看完整版不卡| 久久精品国产综合久久久| 尾随美女入室| 免费久久久久久久精品成人欧美视频| 男人舔女人的私密视频| 国产黄频视频在线观看| 国产 精品1| 99国产综合亚洲精品| 欧美久久黑人一区二区| 男人添女人高潮全过程视频| 99久久精品国产亚洲精品| 国产精品久久久久久久久免| 伊人亚洲综合成人网| 亚洲精品国产av成人精品| 大香蕉久久网| 国产高清国产精品国产三级| 亚洲少妇的诱惑av| 观看av在线不卡| 亚洲精品成人av观看孕妇| 大话2 男鬼变身卡| 精品国产国语对白av| 日韩熟女老妇一区二区性免费视频| 久久精品人人爽人人爽视色| 精品一区在线观看国产| tube8黄色片| 亚洲国产av新网站| 国产精品蜜桃在线观看| 一级,二级,三级黄色视频| 无限看片的www在线观看| 91精品三级在线观看| 亚洲精品国产色婷婷电影| 精品亚洲成国产av| 美女主播在线视频| 建设人人有责人人尽责人人享有的| 午夜福利网站1000一区二区三区| 国产在视频线精品| 国产在线免费精品| av天堂久久9| 久久婷婷青草| 亚洲av福利一区| 伊人久久国产一区二区| 69精品国产乱码久久久| 亚洲欧美一区二区三区国产| 免费观看av网站的网址| 99久久99久久久精品蜜桃| 成年人免费黄色播放视频| 久久影院123| kizo精华| 久久久国产一区二区| 精品亚洲成a人片在线观看| 日韩伦理黄色片| 晚上一个人看的免费电影| 大香蕉久久成人网| 久热爱精品视频在线9| 亚洲一码二码三码区别大吗| 一区二区三区四区激情视频| 国产日韩欧美在线精品| 大香蕉久久成人网| 亚洲国产最新在线播放| 国产淫语在线视频| 中文欧美无线码| av卡一久久| 午夜影院在线不卡| 日韩一本色道免费dvd| 日韩av免费高清视频| 老司机在亚洲福利影院| 成年人免费黄色播放视频| 国产片内射在线| 日韩电影二区| 女人精品久久久久毛片| 国产深夜福利视频在线观看| 日韩伦理黄色片| 欧美亚洲日本最大视频资源| 黄片无遮挡物在线观看| 最近2019中文字幕mv第一页| 亚洲精品aⅴ在线观看| 久久精品亚洲熟妇少妇任你| 看十八女毛片水多多多| 日日爽夜夜爽网站| 亚洲,欧美,日韩| 9色porny在线观看| 蜜桃在线观看..| av天堂久久9| a级片在线免费高清观看视频| 免费av中文字幕在线| 亚洲精品国产av蜜桃| 国产一区二区三区av在线| 国产亚洲一区二区精品| 天天躁夜夜躁狠狠久久av| 中文字幕精品免费在线观看视频| 最近中文字幕2019免费版| av国产久精品久网站免费入址| 国产成人精品在线电影| 99久久综合免费| 90打野战视频偷拍视频| 91精品伊人久久大香线蕉| 五月开心婷婷网| av国产久精品久网站免费入址| 精品人妻熟女毛片av久久网站| 亚洲情色 制服丝袜| 在线精品无人区一区二区三| 国产日韩一区二区三区精品不卡| 亚洲精品一二三| 黄片播放在线免费| 久久av网站| 久久综合国产亚洲精品| 亚洲久久久国产精品| 国产成人av激情在线播放| 十八禁网站网址无遮挡| 高清不卡的av网站| 成年人免费黄色播放视频| 亚洲,一卡二卡三卡| 国产一区二区激情短视频 | 欧美激情高清一区二区三区 | 亚洲七黄色美女视频| 不卡av一区二区三区| 高清黄色对白视频在线免费看| 一边摸一边抽搐一进一出视频| 国产精品久久久久久人妻精品电影 | 青春草亚洲视频在线观看| 亚洲成人国产一区在线观看 | 久久精品久久精品一区二区三区| 国产一区二区三区综合在线观看| 日日爽夜夜爽网站| 性高湖久久久久久久久免费观看| 中国国产av一级| 国产探花极品一区二区| 午夜福利视频精品| 欧美激情高清一区二区三区 | 亚洲成人av在线免费| 亚洲精品国产av成人精品| 制服丝袜香蕉在线| 下体分泌物呈黄色| 这个男人来自地球电影免费观看 | 亚洲激情五月婷婷啪啪| 麻豆精品久久久久久蜜桃| 波多野结衣一区麻豆| 三上悠亚av全集在线观看| 精品一区二区三区四区五区乱码 | 下体分泌物呈黄色| 久久久亚洲精品成人影院| 少妇 在线观看| 只有这里有精品99| 老司机靠b影院| 又大又黄又爽视频免费| 亚洲少妇的诱惑av| 狂野欧美激情性xxxx| 在线观看一区二区三区激情| 日韩av免费高清视频| 亚洲欧美成人精品一区二区| 日韩成人av中文字幕在线观看| 黄色一级大片看看| 天天躁日日躁夜夜躁夜夜| 亚洲一码二码三码区别大吗| 欧美精品亚洲一区二区| 婷婷色麻豆天堂久久| 悠悠久久av| 99热网站在线观看| 国产又爽黄色视频| 激情五月婷婷亚洲| 国产成人一区二区在线| 男女高潮啪啪啪动态图| 亚洲,欧美,日韩| 色播在线永久视频| 国产av码专区亚洲av| 美国免费a级毛片| 男的添女的下面高潮视频| 国产精品国产三级国产专区5o| 日本av免费视频播放| av线在线观看网站| 18禁裸乳无遮挡动漫免费视频| 亚洲国产看品久久| 丰满乱子伦码专区| 老汉色av国产亚洲站长工具| 王馨瑶露胸无遮挡在线观看| 亚洲一码二码三码区别大吗| 秋霞在线观看毛片|