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

    城市人工湖動態(tài)換水水位對流速分布影響

    2018-09-13 01:14:06趙正文馮民權(quán)
    水利水運工程學(xué)報 2018年4期
    關(guān)鍵詞:雁鳴人工湖水深

    趙正文,馮民權(quán),程 剛, 2,唐 洋

    (1. 西安理工大學(xué)省部共建西北旱區(qū)生態(tài)水利工程國家重點實驗室,陜西西安 710048; 2. 西安市水利規(guī)劃勘測設(shè)計院,陜西西安 710054)

    城市人工湖是城市公眾休閑娛樂、親近自然的場所,具有調(diào)節(jié)小氣候、凈化水質(zhì)、維持生物多樣性等生態(tài)功能。城市人工湖具有水淺,流動性差、底泥污染嚴(yán)重,以及接納混雜污水的特點,導(dǎo)致水質(zhì)受到污染,水體富營養(yǎng)化嚴(yán)重。湖泊水動力過程對湖泊中物質(zhì)、能量的輸移擴散起著決定性作用,是研究泥沙、水質(zhì)和生態(tài)等的基礎(chǔ)[1-2],其水動力與湖體本身的形態(tài)、吞吐流、風(fēng)場和水位等物理參數(shù)有關(guān)[3],其中湖泊水位是影響湖泊水文過程和生態(tài)環(huán)境的主要因素[4]。

    人工調(diào)控性湖泊主要是對入湖流量的調(diào)控,自然要素只能對入湖流量產(chǎn)生影響。對于功能性湖泊,如灌溉、防洪、供水和凈化水質(zhì)等,入湖流量是最終影響水位及其變差范圍和時空變化特征的決定性因子[5]。目前,國內(nèi)外學(xué)者通過建立二維、三維淺水動力學(xué)模型[6-8],研究湖泊水動力過程的特征及形成機制,較多關(guān)注風(fēng)場、吞吐流和湖底地形等要素對湖泊水動力的影響[9-12],常常忽略水位的變化。

    水位是湖泊貯水量變化的度量,對于湖泊生態(tài)系統(tǒng)而言,存在一個適宜的水位[13]。研究者主要針對湖泊水位變動對水生植被影響方面取得了一定進(jìn)展,包含對天然大型湖泊水位變動的范圍[14-15]、發(fā)生時間等其他規(guī)律的研究[16-18],但未涉及對水動力影響的研究。天然湖泊水位一般在年內(nèi)隨季節(jié)變化,而人工湖水位則受人為調(diào)控。由于人工湖受人為調(diào)控的水位變化特別明顯,水位調(diào)控是湖泊生態(tài)恢復(fù)與管理的重要技術(shù)[7],因此研究水位對人工湖泊水動力的影響很有必要。

    1 研究方法和方案設(shè)計

    1.1 水動力模型原理

    由于人工湖水平尺度遠(yuǎn)大于垂直尺度,采用二維水動力模型。模型基于不可壓縮和雷諾值平均分布的N-S方程,并服從于Boussinesq假定和靜水壓力的假定[19]。二維淺水方程為:

    ?φ/?t+?p/?x+?q/?y=S

    (1)

    (2)

    (3)

    式中:φ為高程(m);p,q分別為沿x和y方向的單寬流量(m2/s);Vx,Vy分別為沿x,y方向的風(fēng)速(m/s);h為水深(m);C為Chezy阻力系數(shù)(m1/2/s);g為重力加速度(m/s2);τs為風(fēng)摩擦系數(shù);Ω為科氏力系數(shù),受維度影響(1/s);pa為大氣壓強(Pa);ρw為水密度(kg/m3);S,Sx,Sy分別為源質(zhì)量的動量和沿x,y方向的動量分量(kg·m/s);τxx,τxy,τyy為不同方向的有效剪切應(yīng)力;x,y為空間坐標(biāo)(m);t為時間(s)。

    針對人工湖數(shù)值計算面臨復(fù)雜幾何地形和不規(guī)則邊界的問題,建立非結(jié)構(gòu)化三角形網(wǎng)格下的二維水動力模型,將任意三角形作為計算單元,采用單元中心的有限體積法求解,有Ros’s近似黎曼法求解單元界面的對流流動,使用二階TVD格式的限制器函數(shù)避免數(shù)值震蕩。

    圖1 湖泊網(wǎng)格劃分和地形Fig.1 Model grids and bathymetry of lake

    圖2 流速測點分布Fig.2 Velocitymeasurement points

    1.2 研究區(qū)域

    雁鳴湖地處干旱半干旱地帶,湖泊供水由浐河提供,四季供水量變化明顯,其中雁鳴湖5號屬于小型淺水湖泊,面積不大,形狀不規(guī)則,東西寬300m(為x方向),南北長1 200m(為y方向,以西安80為基準(zhǔn)坐標(biāo)),位于浐河旁邊,屬于河道型湖泊,目前湖泊污染較嚴(yán)重。浐河多年平均最小流量0.286m3/s,水流瞬時流量0.235~0.450m3/s,滿足浐河補給雁鳴湖設(shè)計流量0.2m3/s要求。雁鳴湖正常蓄水位434.5m,最低水位433m,水深1.5~2.0m。采用非結(jié)構(gòu)三角形網(wǎng)格對湖泊區(qū)域進(jìn)行劃分,其中網(wǎng)格數(shù)8 084個,節(jié)點數(shù)3 421個,地形采用加權(quán)反離法插值到每個網(wǎng)格節(jié)點。非結(jié)構(gòu)網(wǎng)格能對復(fù)雜幾何地形提供最優(yōu)程度的擬合。為了更好地擬合邊界處的地形變化,對湖岸邊界附近進(jìn)行網(wǎng)格加密和光滑處理,如圖1所示。流速測點布置見圖2。

    1.3 時間步長

    時間步長的設(shè)定是模型穩(wěn)定運行的關(guān)鍵,必須保證CFL數(shù)小于1,對于笛卡爾坐標(biāo)系下的淺水方程式為:

    (4)

    式中:h為水深(m);u和v分別為x和y方向流速分量(m/s);g為重力加速度(m/s2);Δx和Δy分別為x和y方向特征長度,Δx和Δy近似于三角形網(wǎng)格的最小邊長,以三角形網(wǎng)格中心為水深和流速的取值點;Δt為時間間距。

    1.4 邊界條件

    湖泊的邊界條件分別為閉合邊界、開邊界和干濕邊界。閉合邊界指所有正交于邊界流動變量為0的邊界,對應(yīng)于湖泊的陸地邊界。開邊界指有隨時間變化的量的邊界,對于湖泊為進(jìn)出水口的流量過程和水位過程。干濕邊界處理動邊界問題,避免強淺水效應(yīng),對計算區(qū)域內(nèi)水位變化過程,當(dāng)退水水深小于0.005m時,令該網(wǎng)格點為“干點”,不參與水動力計算;當(dāng)水深大于0.005m但小于0.05m時,該網(wǎng)格點僅參與水流連續(xù)方程的計算;當(dāng)水深大于0.1m時,令該網(wǎng)格點為“濕點”,參與水流連續(xù)方程和動量方程的計算。對于雁鳴湖水位變化的動邊界,采用“干濕判別法”處理計算。

    1.5 模型驗證

    該模型水位、流量和流場等水動力指標(biāo)參數(shù)已在2007—2009年進(jìn)行了全面率定和驗證[20],表明模型具有關(guān)鍵過程模擬的優(yōu)勢和能力。此處對2016年10月至2017年9月的水動力過程進(jìn)行進(jìn)一步驗證,選擇湖區(qū)1#和3#及出水口(見圖2)實測水深變化過程進(jìn)行驗證(圖3),同時計算誤差指標(biāo):平均相對誤差(MAE)、均方根誤差(RMSE)、決定系數(shù)(r2)、納什效率系數(shù)(NSE),誤差指標(biāo)分別為:0.0344~0.8540,0.0684~0.8542,0.9132~0.9544,0.933 8~0.9759。對比曲線趨勢和誤差指標(biāo)統(tǒng)計可以看出,模型模擬精度較好,能真實反映雁鳴湖水動力變化過程。水位較高時,整體模擬效果好于低水位,這是因為較高水位時湖泊湖面較大,流場穩(wěn)定,低水位時受湖泊地形影響,流場受到一定程度阻滯作用,從而模擬誤差較大。

    圖3 模型驗證Fig.3 Model validation

    1.6 水位變化及模擬方案

    根據(jù)雁鳴湖調(diào)控?fù)Q水調(diào)研結(jié)果,每年3月為湖泊正常水位,4月湖泊放水,且無引水補給,水位下降明顯,5月在低水位下,開放進(jìn)出口閘門,引入浐河水源對湖底沉積物進(jìn)行沖刷,6月進(jìn)行引水補給,達(dá)到正常水位。湖泊的進(jìn)出流均為緩流,主要影響附近局部水域的流速變化[3],對湖泊整體流速的影響可忽略不計。設(shè)計模擬方案見表1。

    表1 模擬方案Tab.1 Model simulation schemes

    模擬方案主要考慮單因素水位對湖泊流速的影響,其中湖泊進(jìn)出流量為0.2m3/s,風(fēng)向為恒定常年主導(dǎo)風(fēng)、風(fēng)速為年平均風(fēng)速,不考慮降雨和蒸發(fā)。作為對比增加實測換水工況模擬,分析水位和流速時空分布的影響,實測工況中水文和氣象數(shù)據(jù)在2017-03-01—2017-06-30實測得到,水文數(shù)據(jù)包括實測湖泊進(jìn)出口流量(m3/s)和監(jiān)測點水位數(shù)據(jù)(m);氣象數(shù)據(jù)來自國家氣象中心西安站(ID:57131),氣象數(shù)據(jù)包括降雨(mm)、蒸發(fā)(mm)、溫度(℃)、風(fēng)速(m/s)和風(fēng)向。湖泊模擬的地形數(shù)據(jù)由西安市水利規(guī)劃勘測設(shè)計院提供。

    2 結(jié)果與討論

    水位變化是湖泊流場影響因素之一,特別是人工湖泊水位變化明顯,對于湖泊的影響更加顯著,根據(jù)模擬方案和實測換水工況探討水位對人工湖水動力的影響。

    2.1 模擬結(jié)果分析

    圖4為模擬結(jié)果,分別給出了高低水位流速穩(wěn)定時的流速分布與水位升高和降低過程中的湖泊平均流速分布。當(dāng)湖泊處于不同水位時,湖泊的流速差別較大,低水位湖泊流速大于高水位,但環(huán)流分布形態(tài)影響較少。正常蓄水位時,湖泊邊界凹區(qū)域和環(huán)流過渡區(qū)形成死水區(qū);低水位時,由于湖底地形作用在湖底低洼處形成高速流域,并且在進(jìn)出口區(qū)域由于水位頂托作用的減弱也形成高速流域,因而湖泊整體流速增大。在湖泊水位升高或降低過程中,湖泊的平均流速分布相近,高流速分布區(qū)域相近,相對水位升高或降低,流速在進(jìn)水口和出水口差別較大(因為湖泊水位變化主要由進(jìn)出水口流量不同引起,所以導(dǎo)致湖泊流速的差異)。

    圖4 不同方案的流速分布Fig.4 Distribution of flow velocities under different schemes

    圖5 湖泊不同位置的月均流速Fig.5 Monthly average velocities of different points in Yanming Lake

    為了說明水位對湖泊流速的影響,取了湖中不同位置的月平均流速值(圖5)來說明水位的影響程度,其中1#和3#是環(huán)流點位置,2#和4#是無環(huán)流點位置。由圖5可見,湖泊不同位置點的流速存在差異,同一位置點不同水位下的流速也存在差異。在正常水位下湖泊的流速最小,在低水位下流速最大,水位變化下平均流速介于兩者之間,說明低水位有利于湖泊流場,但是作為城市人工湖泊,水位太低,不能滿足生態(tài)環(huán)境的需求,達(dá)不到生態(tài)景觀的水位要求,所以湖泊不可以長時間低水位運行。

    2.2 換水期間水位和流速分析

    2.2.1水位時空差異 基于雁鳴湖水動力模型輸出的水位和水深時空分布如圖6所示。3月份整個湖泊保持較高水位,空間差異不明顯,湖面近似水平,湖泊高水位使水面相應(yīng)增大,湖泊水位對補水的頂托作用,致使湖水水面較平,水深空間分布直接受湖底地形的影響,自上游至下游呈逐級遞增趨勢,出水口的水深最深,局部區(qū)域水深接近2m,湖泊中心到邊界水深由深到淺。4月份湖泊水位出現(xiàn)明顯的水位梯度,水位自上游逐級降低,水域范圍有收縮態(tài)勢,水深變化明顯。5月湖泊水位達(dá)到最低,河灘裸露,水位變化梯度達(dá)到最大,水位差約1m。因水位不斷下降直至湖水歸槽,部分湖流近似河流特性,水位因湖底坡降作用而變化,上游大部分區(qū)域水深極小,呈現(xiàn)裸露狀態(tài);部分區(qū)域形成獨立的水域。6月份水位出現(xiàn)梯度變化,水深增大,水域范圍呈擴大態(tài)勢,水位達(dá)到正常水位??傮w來說,湖泊在3月和5月時水位分布相對穩(wěn)定,部分差異由湖底地形造成;在4月和6月湖泊水位分布從上游至下游呈梯度變化,且湖泊長寬比較大,屬狹長型湖泊,因而湖面比降較大。

    圖6 湖泊水位、水深時空分布Fig.6 Temporal and spatial variation of water level and depth in Yanming Lake

    2.2.2流速時空差異 湖泊流速受地形、吞吐流、風(fēng)場、水位等因素影響。根據(jù)實測換水水位的變化過程和實測風(fēng)場及吞吐流,分析模擬結(jié)果(如圖7)可得:3月湖泊水位較高,湖泊流速較小(0.002m/s左右),高水位使湖泊呈現(xiàn)典型湖相;4月湖泊水位下降,湖泊流速加快(0.005m/s左右),并且湖泊受降雨和風(fēng)場的影響較大;5月湖泊水位達(dá)到最低,湖泊流場為河相流場,由于湖泊地形的影響因素增加,其動力結(jié)構(gòu)更復(fù)雜,流速達(dá)到最大(0.007m/s左右);6月湖泊水位上升,流速減小(0.003m/s左右)。水位變化過程中湖泊形成的相對高流速區(qū)的位置基本不變,環(huán)流作用區(qū)域和進(jìn)水口區(qū)域平均流速相對較高,水位變化過程的平均流速相對大于正常水位的湖泊流速,水位的升高或降低過程實際是湖泊流場的湖相和河相相互轉(zhuǎn)化過程。

    圖7 雁鳴湖流場分布Fig.7 Spatial variation of flow field in Yanming Lake

    圖8 水位變化下日均流速變化Fig.8 Changes in daily average velocities during water level fluctuations

    各采樣點在連續(xù)水位變化下,流速變化見圖8。由圖可見,相對正常水位下(3月)各測點整體流速最??;低水位(5月)時流速最大;其次是水位變化下(4月和6月)的流速。分析實測資料可知,湖泊水動力是多因素綜合作用的結(jié)果,不同因素作用隨時間變化造成流速隨時間變化,不同位置多因素綜合作用導(dǎo)致湖泊流速的差異。在低水位情況下各點流速相近,變化不明顯;3#點流速較大,且流速變化明顯,究其原因是3#處于環(huán)流作用處;實際情況下,變化風(fēng)場的作用導(dǎo)致流速的瞬時增大,例如6月初湖泊流速突然增大,其他各點流速變化趨勢相近。

    2.3 討論

    針對雁鳴湖單因素水位和換水期間水位變化對流場進(jìn)行模擬,結(jié)果表明水位對湖泊水動力有一定影響。隨著湖泊水位降低,水體體積、質(zhì)量減少,需要的外界驅(qū)動力減弱,因而在相同的驅(qū)動力下湖區(qū)流速有所增加[21]。對補水水位升高而言,水面升高產(chǎn)生頂托作用隨接觸面積增大而增大;而在風(fēng)力保持不變的前提下,頂托作用使得補水驅(qū)動力無法滿足水體運動需要,因而整個湖區(qū)的流速便會降低[11]。另外,環(huán)流形成方式隨著水位變化而變化,在正常水位下主要由風(fēng)作用形成湖流,但水位較低時湖底地形和吞吐流作用增強共同形成與風(fēng)向相同的補償流,從而使湖流形式復(fù)雜,環(huán)流之間相互作用使得流速增加。水位不是單獨對水動力產(chǎn)生影響,而是多種湖泊驅(qū)動力共同作用的結(jié)果,不同水位通過外界驅(qū)動力對湖泊水動力產(chǎn)生影響,因此湖泊保持一定水位有利于湖泊水動力的改善。

    水位對人工湖生態(tài)有重要影響,若水位變化不適當(dāng),則會帶來負(fù)面影響,長期固定的高或低水位使湖泊多樣性降低[18]。人工湖在滿足景觀生態(tài)蓄水的要求下,增加水位變化既有利于水生態(tài),也有利于湖泊水動力改善。因此,可以通過調(diào)節(jié)人工湖進(jìn)出湖的流量使水位變化,從而改善生態(tài)和水動力環(huán)境。

    3 結(jié) 語

    基于二維水動力模型,模擬雁鳴湖2017年3—6月期間的水動力變化過程,探究單因素水位對湖泊流場影響,以及換水期間水位和流場的時空分布特征,結(jié)論如下:

    (1)水位是湖泊水動力影響因素之一。不同水位及其變化過程情況下,湖泊流速存在差異,低水位時湖泊平均流速大于正常水位時,水位降低和升高過程中湖泊平均流速介于兩者之間;隨著水位降低,環(huán)流的形成越復(fù)雜。

    (2)換水期間水位和流速分布存在差異。3月份湖泊處于較高水位期,湖區(qū)水位空間分布差異較小,水深自上游至下游呈逐級遞增趨勢,湖泊流速最??;4月和6月份水位出現(xiàn)明顯差異,水深變化明顯,湖泊平均流速大于3月份流速;5月份為低水位期,河灘裸露,水深空間分布差異較大,湖泊流速最大,驗證了單因素水位對湖泊流速的影響。

    (3)湖泊不同水位對流速大小及其分布影響有所不同,在保持生態(tài)景觀水位前提下,人工湖水位變化有利于改善湖泊水動力。

    猜你喜歡
    雁鳴人工湖水深
    書法靜水深流
    河北水利(2022年10期)2022-12-29 11:48:12
    基于水深分段選擇因子的多光譜影像反演水深
    海洋通報(2022年4期)2022-10-10 07:40:32
    Status and Development Trend of Evaluation of Post-Marketing Traditional Chinese Medicines in China
    Clinical Study on the Treatment Efficacy of Cerebral Hemorrhage with Xingnaojing Injection in Real World
    Enlightenment from WHO Pharmacovigilance in Construction of Chinese Materia Medica Pharmacovigilance System
    人工湖生態(tài)景觀建設(shè)和維護
    黑河河道氣動閘人工湖庫尾跌水工程設(shè)計方案探析
    傾聽雁鳴湖畔的人大足音
    人工湖清晨
    快樂游公園
    av在线天堂中文字幕| 人人妻人人澡欧美一区二区 | 12—13女人毛片做爰片一| 香蕉久久夜色| 亚洲av片天天在线观看| 看免费av毛片| 精品久久久久久久人妻蜜臀av | 在线观看66精品国产| 国产伦一二天堂av在线观看| 波多野结衣巨乳人妻| 99国产精品99久久久久| 亚洲成人久久性| 亚洲片人在线观看| 69av精品久久久久久| 三级毛片av免费| 成人18禁高潮啪啪吃奶动态图| 亚洲av成人一区二区三| 不卡一级毛片| 岛国视频午夜一区免费看| 欧美成人性av电影在线观看| www国产在线视频色| 可以在线观看毛片的网站| 麻豆成人av在线观看| 久久国产乱子伦精品免费另类| 老汉色av国产亚洲站长工具| 亚洲人成电影观看| 国产精品一区二区在线不卡| 欧美精品亚洲一区二区| avwww免费| av在线天堂中文字幕| 欧美人与性动交α欧美精品济南到| 午夜久久久在线观看| 国产熟女xx| 亚洲成人久久性| 亚洲激情在线av| 黄频高清免费视频| 免费在线观看亚洲国产| 在线国产一区二区在线| 99久久久亚洲精品蜜臀av| 免费av毛片视频| 亚洲欧美一区二区三区黑人| 国产成人欧美在线观看| 超碰成人久久| 欧美国产日韩亚洲一区| 欧美绝顶高潮抽搐喷水| 国产野战对白在线观看| 露出奶头的视频| 亚洲色图 男人天堂 中文字幕| 久久久久国产一级毛片高清牌| 男男h啪啪无遮挡| 国产精品免费视频内射| 欧美av亚洲av综合av国产av| 淫秽高清视频在线观看| aaaaa片日本免费| 啦啦啦观看免费观看视频高清 | 麻豆国产av国片精品| 黄色丝袜av网址大全| 欧美日韩瑟瑟在线播放| 在线观看一区二区三区| 中亚洲国语对白在线视频| 国产单亲对白刺激| 亚洲自偷自拍图片 自拍| 97碰自拍视频| 国产精品九九99| 国产色视频综合| 日本三级黄在线观看| 欧美乱色亚洲激情| 免费少妇av软件| 日韩欧美免费精品| 精品一区二区三区四区五区乱码| 精品卡一卡二卡四卡免费| 国产精品综合久久久久久久免费 | 亚洲第一欧美日韩一区二区三区| 99国产精品一区二区三区| 韩国精品一区二区三区| 黄色视频,在线免费观看| 日韩有码中文字幕| 热re99久久国产66热| 纯流量卡能插随身wifi吗| 亚洲男人天堂网一区| 亚洲精品久久国产高清桃花| 亚洲情色 制服丝袜| 国产精品1区2区在线观看.| 99re在线观看精品视频| 女性生殖器流出的白浆| 黄频高清免费视频| 又大又爽又粗| 性少妇av在线| 中出人妻视频一区二区| 高潮久久久久久久久久久不卡| 看免费av毛片| 亚洲一区二区三区不卡视频| 国产99白浆流出| 午夜精品国产一区二区电影| 国产精品久久久av美女十八| 久久精品aⅴ一区二区三区四区| 一进一出好大好爽视频| 淫秽高清视频在线观看| 久久精品国产清高在天天线| av中文乱码字幕在线| 欧美激情久久久久久爽电影 | 亚洲av成人不卡在线观看播放网| 999久久久精品免费观看国产| 亚洲精品国产区一区二| www.自偷自拍.com| 国产成人系列免费观看| videosex国产| а√天堂www在线а√下载| 亚洲五月天丁香| 如日韩欧美国产精品一区二区三区| 国产精品美女特级片免费视频播放器 | 在线免费观看的www视频| 国产精品精品国产色婷婷| av网站免费在线观看视频| 欧美日韩亚洲国产一区二区在线观看| 亚洲人成网站在线播放欧美日韩| 色综合亚洲欧美另类图片| 日韩大码丰满熟妇| 丁香六月欧美| 国产亚洲精品第一综合不卡| 大型av网站在线播放| 成熟少妇高潮喷水视频| 久久精品91蜜桃| 88av欧美| 久久香蕉激情| 一级毛片高清免费大全| 婷婷精品国产亚洲av在线| 国产精品自产拍在线观看55亚洲| 免费不卡黄色视频| 亚洲av五月六月丁香网| 国产一区二区三区综合在线观看| 国产成人av教育| 9191精品国产免费久久| 在线视频色国产色| 亚洲国产精品成人综合色| √禁漫天堂资源中文www| 又紧又爽又黄一区二区| 久久久国产精品麻豆| 99久久精品国产亚洲精品| 好看av亚洲va欧美ⅴa在| 久久人人精品亚洲av| 日本免费a在线| 自拍欧美九色日韩亚洲蝌蚪91| 日韩精品中文字幕看吧| 丰满的人妻完整版| 国产精华一区二区三区| 可以在线观看毛片的网站| 日韩欧美国产一区二区入口| 欧美日韩乱码在线| 最近最新中文字幕大全电影3 | 精品国产美女av久久久久小说| 每晚都被弄得嗷嗷叫到高潮| 日韩 欧美 亚洲 中文字幕| 成在线人永久免费视频| 十八禁人妻一区二区| 多毛熟女@视频| 9191精品国产免费久久| 成人国产一区最新在线观看| 亚洲第一青青草原| 狠狠狠狠99中文字幕| 午夜精品久久久久久毛片777| 91成人精品电影| 日韩 欧美 亚洲 中文字幕| 亚洲精品粉嫩美女一区| 日韩高清综合在线| 欧美成人免费av一区二区三区| 可以在线观看毛片的网站| 一级a爱视频在线免费观看| 18禁观看日本| 国产男靠女视频免费网站| 欧美午夜高清在线| 成人亚洲精品一区在线观看| 母亲3免费完整高清在线观看| 久久久国产成人免费| 日韩欧美在线二视频| 亚洲中文日韩欧美视频| 午夜两性在线视频| 中文字幕人成人乱码亚洲影| 成人亚洲精品一区在线观看| 好男人在线观看高清免费视频 | av片东京热男人的天堂| 日日爽夜夜爽网站| 变态另类成人亚洲欧美熟女 | 两人在一起打扑克的视频| 亚洲色图综合在线观看| 999久久久精品免费观看国产| 国产精品影院久久| 一本大道久久a久久精品| 国产亚洲精品久久久久久毛片| 岛国视频午夜一区免费看| 香蕉国产在线看| 国产人伦9x9x在线观看| 美女国产高潮福利片在线看| 成人三级黄色视频| av有码第一页| 亚洲色图 男人天堂 中文字幕| 我的亚洲天堂| 怎么达到女性高潮| 亚洲av电影不卡..在线观看| 桃色一区二区三区在线观看| 两人在一起打扑克的视频| 免费不卡黄色视频| 免费不卡黄色视频| 国产亚洲欧美98| 亚洲成a人片在线一区二区| 精品国产一区二区久久| 国产成人av教育| 亚洲欧美激情综合另类| 国产区一区二久久| 国产精品亚洲av一区麻豆| www国产在线视频色| 精品欧美国产一区二区三| 亚洲电影在线观看av| 亚洲欧美日韩另类电影网站| 日韩成人在线观看一区二区三区| 亚洲精品粉嫩美女一区| 成人手机av| 制服人妻中文乱码| 精品福利观看| 亚洲国产毛片av蜜桃av| 大陆偷拍与自拍| 日韩有码中文字幕| 亚洲精品av麻豆狂野| 亚洲精品中文字幕在线视频| 中文字幕人妻丝袜一区二区| 高清毛片免费观看视频网站| 丰满人妻熟妇乱又伦精品不卡| 欧美不卡视频在线免费观看 | e午夜精品久久久久久久| 欧美乱色亚洲激情| 国产高清激情床上av| 久久天躁狠狠躁夜夜2o2o| 美女扒开内裤让男人捅视频| 韩国精品一区二区三区| 给我免费播放毛片高清在线观看| 成年版毛片免费区| 黑丝袜美女国产一区| 丰满人妻熟妇乱又伦精品不卡| www.自偷自拍.com| 禁无遮挡网站| 777久久人妻少妇嫩草av网站| 亚洲精品美女久久久久99蜜臀| 女人被狂操c到高潮| 亚洲成人免费电影在线观看| 日韩av在线大香蕉| 久久精品亚洲精品国产色婷小说| 久久久水蜜桃国产精品网| 精品欧美一区二区三区在线| 成人av一区二区三区在线看| 久久久久国内视频| 天天一区二区日本电影三级 | 91九色精品人成在线观看| 亚洲视频免费观看视频| 欧美日韩亚洲国产一区二区在线观看| 女人高潮潮喷娇喘18禁视频| 在线观看免费日韩欧美大片| 亚洲av五月六月丁香网| 少妇裸体淫交视频免费看高清 | www.www免费av| 国产亚洲精品综合一区在线观看 | 亚洲成人精品中文字幕电影| 一本久久中文字幕| 国语自产精品视频在线第100页| 亚洲无线在线观看| 久久青草综合色| 国产真人三级小视频在线观看| 99香蕉大伊视频| 成人国产综合亚洲| 免费女性裸体啪啪无遮挡网站| 黄片播放在线免费| 一进一出好大好爽视频| 欧美另类亚洲清纯唯美| 亚洲国产毛片av蜜桃av| 性少妇av在线| 大型黄色视频在线免费观看| 国产亚洲av嫩草精品影院| 在线十欧美十亚洲十日本专区| 日韩三级视频一区二区三区| 久久午夜综合久久蜜桃| 在线观看免费视频网站a站| 三级毛片av免费| 午夜两性在线视频| 亚洲国产精品成人综合色| 热re99久久国产66热| 国产精品免费一区二区三区在线| 久久香蕉激情| 欧美乱码精品一区二区三区| 国产成人影院久久av| 色综合欧美亚洲国产小说| av免费在线观看网站| 1024香蕉在线观看| 男人操女人黄网站| 91九色精品人成在线观看| 精品人妻在线不人妻| 精品国产一区二区久久| 黄色成人免费大全| 在线观看66精品国产| 精品久久久久久久毛片微露脸| 婷婷丁香在线五月| 国产精品一区二区精品视频观看| 国产蜜桃级精品一区二区三区| 99国产精品一区二区蜜桃av| 亚洲男人的天堂狠狠| 国产精品99久久99久久久不卡| 国产私拍福利视频在线观看| 亚洲精品久久国产高清桃花| 亚洲 欧美一区二区三区| 成人国产一区最新在线观看| 叶爱在线成人免费视频播放| 欧美一级毛片孕妇| 亚洲午夜精品一区,二区,三区| 91麻豆精品激情在线观看国产| 欧美成狂野欧美在线观看| 黄色毛片三级朝国网站| 久久国产乱子伦精品免费另类| 国产精品久久久久久人妻精品电影| 中文字幕最新亚洲高清| 国产一区二区三区视频了| 淫秽高清视频在线观看| 亚洲avbb在线观看| 亚洲精品在线美女| 在线观看午夜福利视频| 一级a爱视频在线免费观看| 国产伦人伦偷精品视频| 精品久久久精品久久久| 国产精品爽爽va在线观看网站 | 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av成人不卡在线观看播放网| 一进一出抽搐动态| 大型黄色视频在线免费观看| 日本一区二区免费在线视频| 久久久久久人人人人人| 精品熟女少妇八av免费久了| 国产精品综合久久久久久久免费 | 欧美日韩瑟瑟在线播放| 亚洲精品中文字幕一二三四区| 9191精品国产免费久久| 日韩视频一区二区在线观看| 黑人操中国人逼视频| 亚洲熟妇熟女久久| 成人欧美大片| 国产午夜精品久久久久久| 九色国产91popny在线| 午夜福利一区二区在线看| 亚洲第一欧美日韩一区二区三区| 国产成人系列免费观看| 99久久综合精品五月天人人| 97人妻精品一区二区三区麻豆 | 精品久久蜜臀av无| 亚洲国产精品999在线| 美女 人体艺术 gogo| 国产成人欧美| cao死你这个sao货| 免费观看精品视频网站| 国产99久久九九免费精品| 亚洲午夜理论影院| 午夜福利18| 亚洲国产欧美一区二区综合| 最近最新中文字幕大全电影3 | av天堂久久9| 午夜精品国产一区二区电影| 狠狠狠狠99中文字幕| 午夜免费鲁丝| 亚洲成人精品中文字幕电影| www.www免费av| 精品卡一卡二卡四卡免费| 男女做爰动态图高潮gif福利片 | 男男h啪啪无遮挡| 午夜福利一区二区在线看| 日韩有码中文字幕| 亚洲色图av天堂| 亚洲无线在线观看| 亚洲第一电影网av| 婷婷精品国产亚洲av在线| 亚洲精品久久国产高清桃花| 中文字幕高清在线视频| 欧美黄色片欧美黄色片| 啦啦啦观看免费观看视频高清 | 美国免费a级毛片| 精品久久久久久久毛片微露脸| 中文字幕高清在线视频| 俄罗斯特黄特色一大片| 国产成+人综合+亚洲专区| 国产精品精品国产色婷婷| 亚洲av片天天在线观看| 中文字幕另类日韩欧美亚洲嫩草| 欧美精品啪啪一区二区三区| 久久亚洲真实| 一边摸一边做爽爽视频免费| 午夜精品在线福利| 亚洲三区欧美一区| 一夜夜www| 国产精品一区二区免费欧美| 午夜影院日韩av| 在线观看日韩欧美| 少妇粗大呻吟视频| 美女 人体艺术 gogo| 黑人巨大精品欧美一区二区蜜桃| 亚洲成a人片在线一区二区| 人人妻人人澡欧美一区二区 | 日韩欧美一区视频在线观看| 国产高清videossex| 国产精品久久久久久精品电影 | 国产精品1区2区在线观看.| a级毛片在线看网站| 变态另类成人亚洲欧美熟女 | 大型黄色视频在线免费观看| 欧美色欧美亚洲另类二区 | 国产精品永久免费网站| 午夜免费观看网址| 怎么达到女性高潮| 90打野战视频偷拍视频| 熟妇人妻久久中文字幕3abv| 满18在线观看网站| 黑人巨大精品欧美一区二区mp4| 国产成人影院久久av| 精品国产超薄肉色丝袜足j| 好男人在线观看高清免费视频 | 91精品三级在线观看| 看免费av毛片| 免费高清在线观看日韩| 免费在线观看视频国产中文字幕亚洲| 啪啪无遮挡十八禁网站| 亚洲熟女毛片儿| 如日韩欧美国产精品一区二区三区| 悠悠久久av| 国产成人av教育| 999久久久国产精品视频| 午夜精品在线福利| 动漫黄色视频在线观看| 午夜精品久久久久久毛片777| 亚洲精品国产区一区二| 国产精品美女特级片免费视频播放器 | 精品电影一区二区在线| 很黄的视频免费| 婷婷精品国产亚洲av在线| 十分钟在线观看高清视频www| 亚洲欧美日韩无卡精品| 国产精品久久电影中文字幕| 国产精品免费一区二区三区在线| 黄片播放在线免费| 亚洲五月婷婷丁香| 亚洲人成网站在线播放欧美日韩| 黄片大片在线免费观看| 老司机靠b影院| 99精品久久久久人妻精品| 日韩免费av在线播放| 免费在线观看黄色视频的| 成人国产一区最新在线观看| 少妇的丰满在线观看| 美女午夜性视频免费| 日韩欧美在线二视频| 国产免费av片在线观看野外av| 天堂影院成人在线观看| 老汉色av国产亚洲站长工具| 男男h啪啪无遮挡| 久久中文字幕人妻熟女| 曰老女人黄片| 黄片大片在线免费观看| 国产亚洲av高清不卡| 给我免费播放毛片高清在线观看| 999久久久国产精品视频| 精品国产一区二区久久| 久久久国产欧美日韩av| 亚洲成国产人片在线观看| 亚洲国产看品久久| 亚洲第一电影网av| 男女下面进入的视频免费午夜 | 亚洲成国产人片在线观看| 国产xxxxx性猛交| 久久伊人香网站| 1024香蕉在线观看| 一级毛片精品| 欧美日韩亚洲综合一区二区三区_| 日韩国内少妇激情av| 一个人观看的视频www高清免费观看 | 精品国产一区二区三区四区第35| 日韩欧美国产在线观看| 亚洲中文av在线| 熟妇人妻久久中文字幕3abv| 亚洲精品久久国产高清桃花| 日韩免费av在线播放| 欧美黄色淫秽网站| 色综合亚洲欧美另类图片| 免费观看精品视频网站| 少妇裸体淫交视频免费看高清 | 法律面前人人平等表现在哪些方面| 亚洲激情在线av| 日日爽夜夜爽网站| 亚洲七黄色美女视频| 一二三四在线观看免费中文在| 欧美色视频一区免费| 久久欧美精品欧美久久欧美| 国产激情久久老熟女| 欧美老熟妇乱子伦牲交| 一级毛片高清免费大全| 99国产精品一区二区蜜桃av| 亚洲精品国产区一区二| 欧美黄色淫秽网站| 18禁国产床啪视频网站| 精品卡一卡二卡四卡免费| 99国产精品一区二区蜜桃av| 身体一侧抽搐| 国产单亲对白刺激| 在线免费观看的www视频| 国产精品野战在线观看| 亚洲一码二码三码区别大吗| 亚洲一区二区三区色噜噜| 国产97色在线日韩免费| 1024香蕉在线观看| 99久久99久久久精品蜜桃| 啪啪无遮挡十八禁网站| 一进一出抽搐gif免费好疼| 欧美黄色淫秽网站| 在线免费观看的www视频| 亚洲午夜理论影院| 视频在线观看一区二区三区| 自线自在国产av| 欧美成人免费av一区二区三区| 岛国视频午夜一区免费看| 欧美日本视频| 一夜夜www| 亚洲成人免费电影在线观看| 少妇裸体淫交视频免费看高清 | 精品人妻1区二区| 日本撒尿小便嘘嘘汇集6| 日本免费一区二区三区高清不卡 | 亚洲中文字幕一区二区三区有码在线看 | 亚洲一码二码三码区别大吗| 欧美性长视频在线观看| 国产亚洲精品久久久久5区| 在线观看舔阴道视频| 极品教师在线免费播放| 香蕉久久夜色| 免费观看人在逋| 国产伦一二天堂av在线观看| 免费观看人在逋| 日韩国内少妇激情av| 日日干狠狠操夜夜爽| 岛国在线观看网站| 亚洲黑人精品在线| av免费在线观看网站| 午夜福利在线观看吧| 成人特级黄色片久久久久久久| a在线观看视频网站| 亚洲精品在线观看二区| 伦理电影免费视频| or卡值多少钱| 国产亚洲精品久久久久久毛片| aaaaa片日本免费| 久久精品国产清高在天天线| 女人爽到高潮嗷嗷叫在线视频| 18禁裸乳无遮挡免费网站照片 | 精品国产一区二区三区四区第35| 一级毛片女人18水好多| cao死你这个sao货| 好男人电影高清在线观看| 国产激情久久老熟女| 亚洲av日韩精品久久久久久密| 国产精品98久久久久久宅男小说| 美女 人体艺术 gogo| 国产精品久久久人人做人人爽| 看片在线看免费视频| 最近最新免费中文字幕在线| 女人被躁到高潮嗷嗷叫费观| 亚洲男人的天堂狠狠| 国产一区二区在线av高清观看| 亚洲色图av天堂| 欧美乱妇无乱码| 精品久久久久久,| 国内精品久久久久精免费| 满18在线观看网站| 18禁观看日本| 国产片内射在线| 久久中文字幕人妻熟女| 久久久久久免费高清国产稀缺| 亚洲第一欧美日韩一区二区三区| АⅤ资源中文在线天堂| 亚洲色图av天堂| 精品国产一区二区久久| 18禁国产床啪视频网站| 欧美 亚洲 国产 日韩一| 69精品国产乱码久久久| 免费少妇av软件| 免费一级毛片在线播放高清视频 | ponron亚洲| svipshipincom国产片| 大陆偷拍与自拍| 欧美日韩黄片免| 亚洲精品国产一区二区精华液| 女生性感内裤真人,穿戴方法视频| 九色亚洲精品在线播放| e午夜精品久久久久久久| 日本精品一区二区三区蜜桃| 国产精品国产高清国产av| 操出白浆在线播放| 久久婷婷人人爽人人干人人爱 | 免费人成视频x8x8入口观看| 韩国精品一区二区三区| 两个人看的免费小视频| 亚洲av熟女| 日韩欧美国产一区二区入口| 欧美大码av| 精品免费久久久久久久清纯| 日韩免费av在线播放| 亚洲天堂国产精品一区在线| 黄频高清免费视频| 精品电影一区二区在线| 午夜久久久久精精品| 久久人妻熟女aⅴ| 18禁观看日本| 久久人妻av系列| 免费看美女性在线毛片视频| 国产一区二区三区视频了| 午夜成年电影在线免费观看|