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

    水庫水溫模擬預(yù)測中常見問題探討

    2016-09-05 04:57:23李克鋒梁瑞峰脫友才鄧云李嘉
    環(huán)境影響評價 2016年3期
    關(guān)鍵詞:水溫氣溫水庫

    李克鋒,梁瑞峰,脫友才,鄧云,李嘉

    (四川大學(xué)水力學(xué)與山區(qū)河流開發(fā)保護(hù)國家重點實驗室,四川成都 610065)

    ?

    水庫水溫模擬預(yù)測中常見問題探討

    李克鋒,梁瑞峰,脫友才,鄧云,李嘉

    (四川大學(xué)水力學(xué)與山區(qū)河流開發(fā)保護(hù)國家重點實驗室,四川成都610065)

    針對影響水庫水溫模擬效果的常見問題,提出以水溫-氣溫相關(guān)法確定缺乏資料地區(qū)的現(xiàn)狀水溫,對水庫地形校正效果開展定量分析,討論了氣溫直減率對水氣界面長波輻射、蒸發(fā)、傳導(dǎo)的影響程度,對比了部分氣象經(jīng)驗公式的熱通量計算差異,以二維模型為例定性分析了零方程與雙方程紊流模式、靜水壓力假定與完整垂向動量方程的差異,提出了可通過熱量收支平衡、逐月徑流量與庫容關(guān)系等方法定量檢驗?zāi)M結(jié)果的合理性。

    水庫水溫;地形;氣象;數(shù)學(xué)模型

    水電資源作為一種清潔能源一直是國際社會能源行業(yè)的重點發(fā)展領(lǐng)域,但水電站的興建也極大地改變了河道水溫的自然變化過程,特別是在水能資源豐富的河流,巨型電站或者多個梯級水庫聯(lián)合運行時,對河道水溫的影響范圍更寬,對區(qū)域性的環(huán)境影響更大。水溫影響逐漸成為水電開發(fā)過程中環(huán)境影響評價的重要內(nèi)容[1-3]。

    目前,國內(nèi)外水庫水溫影響研究一般從數(shù)值模擬、室內(nèi)試驗、原位觀測3個方面開展,其中數(shù)值模擬方法具有物理機(jī)制清晰、考慮因素全面、可有效結(jié)合理論分析與實驗研究成果等優(yōu)勢,越來越多地應(yīng)用于新建項目的環(huán)境影響評價中,其模擬精度直接關(guān)系到對工程環(huán)境影響的正確認(rèn)識。但水溫的模擬由于受模型選擇、參數(shù)取值、離散方法、邊界條件、初始條件等多種因素的影響,在實際應(yīng)用中如果處理不當(dāng)可能影響模擬效果。本文擬對水庫水溫模擬中常見的主要問題展開分析和討論,以期為相關(guān)研究提供參考。

    1 現(xiàn)狀水溫的確定

    現(xiàn)狀水溫作為數(shù)值模擬的輸入條件和水溫影響程度的對比基礎(chǔ),是模擬的重要基礎(chǔ)資料。在平原河流或者具有布設(shè)水文儀器條件的河段,水溫作為常規(guī)水文監(jiān)測項目,其數(shù)據(jù)獲取尚不至于制約模擬預(yù)測。隨著水電開發(fā)逐漸向高海拔寒冷地區(qū)拓展,研究河段缺乏水溫觀測資料,難以確定模擬所需的現(xiàn)狀水溫,但通過氣象相關(guān)、鄰近流域監(jiān)測數(shù)據(jù)類比等方法,仍可得到合理的基礎(chǔ)水溫數(shù)據(jù)[4]。

    氣象條件是影響天然水溫變化的重要因素。河段沿程有氣象監(jiān)測數(shù)據(jù)且周邊氣候相近的流域存在水溫監(jiān)測數(shù)據(jù)時,可采用水溫-氣溫相關(guān)方法確定該區(qū)域的現(xiàn)狀水溫。以西藏自治區(qū)境內(nèi)河流為例,怒江上游僅嘉玉橋、龍普、扎那水文站曾短期監(jiān)測水溫,現(xiàn)有數(shù)據(jù)難以滿足上游規(guī)劃梯級電站現(xiàn)狀水溫確定的需要,但可根據(jù)西藏境內(nèi)其他流域的氣溫、水溫關(guān)系進(jìn)行推算。

    水溫-氣溫相關(guān)關(guān)系可表示為:

    T水=a×T氣+b

    (1)

    式中,a為水溫隨氣溫的變化率;b為某處水溫受上游來流水溫的影響程度,℃。

    西藏自治區(qū)雅魯藏布江、年楚河、拉薩河、尼洋河、怒江等干支流的水溫-氣溫相關(guān)系數(shù)見表1。從相關(guān)性來看,各水文站水溫均與氣溫呈高度相關(guān)關(guān)系(相關(guān)系數(shù)>0.8);水溫隨氣溫的變化率平均為0.757 9;各干支流相關(guān)關(guān)系中的參數(shù)b均為下游大于上游,說明各河流下游水溫受上游影響程度存在差異性。因此,可采用式(1)計算怒江干支流的天然水溫。

    此外,缺乏水溫觀測資料地區(qū)的現(xiàn)狀水溫的確定,也可直接移用鄰近流域流量相近站點的監(jiān)測資料,對一些監(jiān)測年限較短的資料,還需要與氣象數(shù)據(jù)建立相關(guān)關(guān)系后進(jìn)一步插補延長。

    2 地形的概化與校正

    由地形概化得到的網(wǎng)格是水庫庫容、面積等幾何尺寸在數(shù)學(xué)模型中的直接體現(xiàn),因而地形精度一定程度上決定了預(yù)測精度。從大斷面、地形圖等概化后的水庫庫容曲線受庫區(qū)形狀不規(guī)則、測量精度、提取方法等影響,與實際庫容、面積曲線存在差異,因而在計算中應(yīng)對概化庫容、面積曲線做進(jìn)一步校正,消除地形誤差帶來的水庫調(diào)節(jié)能力變化。

    浯溪口水利樞紐工程是江西省景德鎮(zhèn)市昌江干流中游一座以防洪為主,兼顧供水、發(fā)電等的綜合利用工程。該工程正常蓄水位為56 m,水庫總庫容為4.274×108m3,具有不完全年調(diào)節(jié)性能。圖1和圖2顯示了浯溪口水庫地形校正對水庫水溫的影響。浯溪口庫區(qū)地形系根據(jù)分辨率較低的地形圖得到,經(jīng)與設(shè)計院提供的高程-面積-容積曲線對比,從地形圖概化得到的正常蓄水位水面面積為17.66 km2,與實際水面面積誤差接近30%。對此地形不做校正時,采用立面二維水溫模型計算得到的平水年水庫下泄水溫最大升幅達(dá)6.0℃;經(jīng)對地形開展復(fù)核,根據(jù)設(shè)計方提供的復(fù)核后地形校正模型計算采用的庫容面積曲線(圖1),正常蓄水位水面面積校正為24.70 km2,下泄水溫的最大升幅減小至3.3℃(圖2)。復(fù)核前后水面面積差異帶來水面熱通量計算存在較大差異,特別是水庫蓄水后期的高水位運行,面積差異更趨放大,顯著改變了水庫的冬季蓄熱能力,導(dǎo)致不合理的下泄水溫大幅升高現(xiàn)象。

    圖1 地形校正前后浯溪口水庫水面面積Fig.1 Surface area of Wuxikou Reservoir before and after the topographical correction

    圖2 地形修正前后浯溪口水庫下泄水溫Fig.2 Discharge water temperature of Wuxikou Reservoir before and after the topographical correction

    3 氣象資料與計算公式的選擇

    影響水庫水溫模擬的氣象要素一般包括太陽輻射、氣溫、云量、風(fēng)速風(fēng)向、相對濕度等。

    太陽輻射和氣溫是影響水體熱量收支的最重要因素。河谷深切的山區(qū)水庫應(yīng)充分考慮兩側(cè)山體遮蔽對輻射到達(dá)水面時段的影響,特別是南北走向的河流。因地形遮蔽作用,水面接收的太陽直射輻射最大可減少25%左右[5]。

    氣象站點的選擇也間接影響到水面熱通量的計算。金沙江下游的溪洛渡庫區(qū)有金陽、永善等氣象站,金陽站位于庫區(qū)中段,在地理位置上更有代表性,但金陽氣象站海拔1 450 m,與溪洛渡正常蓄水位600 m相差太大,而氣溫具有隨海拔上升而降低的特性,直接采用金陽站氣溫計算水氣界面熱通量將會帶來較大誤差。

    圖3為金陽氣象站氣溫修正前后對溪洛渡庫區(qū)河段水氣傳導(dǎo)、蒸發(fā)、長波輻射的影響程度。在采用該地區(qū)山區(qū)氣溫直減率5℃/km[6]將1 450 m高程的氣溫修正至600 m高程后,該區(qū)域年均氣溫比下游永善站高出3.2℃,河段水氣界面熱通量年均增加達(dá)73.7 W/m2,也從側(cè)面說明氣象因素是金沙江溪洛渡段增溫率(0.43℃/100 km)高于向家壩段增溫率(0.07℃/100 km)的原因之一。

    圖3 金沙江下游溪洛渡庫區(qū)河段氣溫修正對水氣界面熱量的影響Fig.3 Impact of air temperature correction on heat flux in Xiluodu Reservoir

    此外,水面熱通量經(jīng)驗公式的選取也對計算存在較大影響,本文以大氣長波輻射為例進(jìn)行說明。大氣中的水珠、蒸汽、塵埃等可以吸收來自太陽的短波輻射等,使其自身增溫,然后又通過長波輻射形式輻射到水面,其強(qiáng)度取決于氣溫的高低。大氣長波輻射絕大部分可被水體吸收,約有3%被水面反射到大氣中。大氣長波輻射波長范圍為4~120 μm,輻射強(qiáng)度一般取決于氣溫和云量。

    大氣長波輻射計算公式為:

    φan=(1-γa)·σ·εa(273+Ta)4

    (2)

    式中,Ta為水面上2 m處的氣溫;γa為長波反射率,取0.03;σ為Stefan-Boltzman常數(shù),5.67×10-8(W/m2·K4);εa為大氣發(fā)射率。

    大氣發(fā)射率的常用計算公式為:

    (3)

    (4)

    (5)

    式中,C為云層覆蓋率;Kε與云層高度有關(guān),均值為0.17;ea為水面上空氣蒸發(fā)壓力;Aa為經(jīng)驗系數(shù),取值范圍為0.5~0.7,常取值0.6;Ab為經(jīng)驗系數(shù),取值范圍為0.031~0.076 mmHg-0.5,常取值0.031;eair為空氣中的水蒸氣分壓力,mmHg。

    空氣中的水蒸氣分壓力可用下式表示:

    (6)

    式中,Td為露點溫度。

    仍以浯溪口水庫庫區(qū)氣象條件來說明公式(3)、(4)、(5)的差異(圖4)。從圖4中大氣長波輻射在一年內(nèi)的變化可發(fā)現(xiàn),采用不同的大氣發(fā)射率經(jīng)驗公式,冬季的差異小于夏季,但最小也相差36.7 W/m2;夏季差異較大,6月差異高達(dá)97.3 W/m2。

    短波輻射、水體長波輻射、蒸發(fā)通量、傳導(dǎo)通量在選擇不同的計算公式時,也存在與大氣長波輻射計算類似的問題,在預(yù)測中宜根據(jù)研究區(qū)域的氣候特征對氣象公式作適用性驗證。在商業(yè)軟件中,由于熱通量計算模式相對固定,尤其應(yīng)注意公式的適用性。

    圖4 浯溪口庫區(qū)采用不同大氣發(fā)射率公式的水面長波輻射通量Fig.4 Water surface longwave radiation flux under different air emission rate formulae in Wuxikou Reservoir

    4 垂向動量交換的模擬

    應(yīng)用于水庫水溫預(yù)測的數(shù)學(xué)模型主要從空間維度上分為一維、二維、三維??v向一維模型可用于水深較淺的水庫,垂向一維模型用于深且短的水庫;二維模型主要有深度平均的平面模型和寬度平均的立面模型;三維數(shù)學(xué)模型則能充分反映水庫水溫的變化細(xì)節(jié),特別是泄流孔口附近水流及水溫的三維特征[7]。由于三維水溫計算工作量巨大,目前主要用于近壩區(qū)電站進(jìn)口處的模擬,而一維模型適用性有限,因此主要就立面二維水溫模型作簡要對比。

    水庫的流場分布與水庫的溫度分布具有強(qiáng)烈的相互作用,準(zhǔn)確模擬水庫溫度場的關(guān)鍵在于精確模擬浮力流流場,因此立面二維模型之間的差異主要是紊流模式和垂向動量方程的差異。

    根據(jù)紊流模型采用的微分輸運方程的個數(shù),可以分為零方程、單方程、雙方程3類。美國陸軍工程兵團(tuán)的CE-QUAL-W2模型推薦采用W2零方程模式計算渦粘系數(shù):

    (7)

    式中,νt為垂向渦流系數(shù);κ為范卡門常數(shù);lm為混合長度;u為縱向流速;τwy為因風(fēng)而產(chǎn)生的橫向剪應(yīng)力;k為波數(shù);τytributary為因支流入流而產(chǎn)生的橫向剪應(yīng)力;ρ為水體密度;C為常數(shù);Δzmax為垂向網(wǎng)格間距的最大值;Ri為理查森數(shù)。

    W2模式通過Ri考慮了垂向密度梯度帶來的紊動粘性變化,混合長則由網(wǎng)格垂向最大間距確定,這就使紊動粘性系數(shù)νt的確定依賴于網(wǎng)格劃分。W2可通過原位觀測率定相對合理的混合長,但率定難免將水溫實測中的不確定因素帶入其中,影響模型的適用性。

    單方程模式增加紊動動能的精確輸運方程,其難點仍在于長度比尺的確定。雙方程模型解決了單方程模型長度比尺的確定難題,但在使用中如采用完全的不可壓縮假定,垂向溫差引起的密度流不能通過浮力項影響渦粘系數(shù),也將導(dǎo)致浮力流模擬偏差。

    立面二維模型垂向動量方程的差異主要在于是否采用完整的垂向動量方程。完整的垂向動量方程和靜水壓力假定動量方程分別為:

    (8)

    (9)

    式中,w為垂向流速;p為壓強(qiáng);T為水溫;ρa為參考密度;B為河寬;g為重力加速度;β為等壓膨脹系數(shù);νe為分子粘性系數(shù)ν與紊動渦粘系數(shù)νt之和。

    與完整的垂向動量方程相比,靜水壓力假定動量方程忽略了垂向流速的時變項、對流項、擴(kuò)散項,使得模型主要適用于流速的各向梯度變化不大的流動模擬,而當(dāng)模擬存在明顯上浮或者下潛現(xiàn)象的溫差異重流時,模擬結(jié)果存在偏差。

    5 預(yù)測結(jié)果合理性分析

    水庫水溫模擬預(yù)測的結(jié)果是否合理,需要通過對熱量平衡、延遲時間等進(jìn)行分析來確定。分析結(jié)果如表2、表3所示。

    表2列出了金沙江下游烏東德水庫在平水年單獨運行時的熱量收支情況。從表2可知,烏東德水庫在6—9月入庫、出庫熱量占全年熱量收支的65%以上,水庫3—5月入庫、出庫熱量的全年占比雖小,但出庫熱量總體小于入庫熱量,反映了在此期間水庫下泄低溫水現(xiàn)象的存在。

    表2 金沙江下游烏東德水庫平水年的熱量收支情況

    表3 土溪口水庫典型平水年逐月徑流量與庫容的替換關(guān)系

    對于調(diào)節(jié)性能較強(qiáng)的水庫,其下泄水流的低溫水降幅可通過來流到達(dá)壩前的延遲時間做輔助分析。表3為四川宣漢縣渠江支流前河土溪口水庫的逐月徑流量與庫容的替換關(guān)系,替換次數(shù)與延遲時間均采用當(dāng)月水位下的全庫容計算,即假定水庫整個斷面參與流動。表3中可見水庫來流量年內(nèi)分布極不均勻,典型平水年在3—5月來流量僅占全年徑流量的25%,庫中水體難以充分替換,11月至次年4月,當(dāng)月來流均不能完全替換當(dāng)月庫內(nèi)水體;水庫調(diào)節(jié)引起的低溫來流延遲出庫是較大低溫水降幅的主要成因。

    6 結(jié)語

    水溫模擬受模型選擇、參數(shù)取值、離散方法、邊界條件、初始條件等多種因素的影響,在實際應(yīng)用中如處理不當(dāng),可能影響模擬效果。本文對水庫水溫模擬預(yù)測中常見的主要問題開展了討論,對直接影響模擬過程的現(xiàn)狀水溫的確定、地形的概化與校正、氣象資料的處理與計算公式的選擇、垂向動量交換的模擬進(jìn)行了分析,并提出了具體方法與處理思路,對采用熱量收支平衡、逐月徑流量與庫容關(guān)系等方法分析模擬結(jié)果的合理性進(jìn)行了闡述。

    由于數(shù)值模擬包括前后處理、離散方法與過程、計算控制、驗證分析等多個環(huán)節(jié),文中所涉僅為少數(shù)常見問題,實際處理中需要在更多方面加以關(guān)注,如垂向一維計算中綜合擴(kuò)散系數(shù)的確定、河流水溫計算中熱源項處理、二維計算中水庫泄流孔口喇叭口的概化、泄流孔口及周邊最小網(wǎng)格尺寸的確定、多維嵌套中的界面?zhèn)鬟f、分層取水過程中的動邊界處理、流場與溫度場的匹配驗證等。預(yù)測結(jié)果的合理性分析在水庫水溫預(yù)測中也至關(guān)重要。環(huán)評中可結(jié)合水庫入流和出流的熱量平衡分析、分時段的庫水交換次數(shù)分析等方法,總體上判斷預(yù)測結(jié)果的合理性。此外,還可將預(yù)測結(jié)果與具有可類比性的已建水庫水溫觀測資料進(jìn)行對比,來分析預(yù)測結(jié)果的合理性。

    [1]宋策, 周孝德, 辛向文. 龍羊峽水庫水溫結(jié)構(gòu)演變及其對下游河道水溫影響[J]. 水科學(xué)進(jìn)展, 2011, 22(3): 421-428.

    [2]李鐘順, 陳永燦, 劉昭偉, 等. 密云水庫水溫分布特征[J]. 清華大學(xué)學(xué)報(自然科學(xué)版), 2012, 52(6): 798-803.

    [3]脫友才, 劉志國, 鄧云, 等. 豐滿水庫水溫的原型觀測及分析[J]. 水科學(xué)進(jìn)展, 2014, 25(5): 739-746.

    [4]李克鋒, 郝紅升, 莊春義, 等. 利用氣象因子估算天然河道水溫的新公式[J]. 四川大學(xué)學(xué)報(工程科學(xué)版), 2006, 38(1): 1-4.

    [5]張艷林, 程國棟, 李新, 等. 山區(qū)太陽輻射對水熱過程影響的敏感性分析[J]. 冰川凍土, 2012, 34(3): 650-659.

    [6]翁篤鳴, 孫治安. 我國山地氣溫直減率的初步研究[J]. 地理研究, 1984, 3(2): 24-34.

    [7]馬方凱, 江春波. 三峽水庫近壩區(qū)三維流場及溫度場的數(shù)值模擬[J]. 水利水電科技進(jìn)展, 2007, 27(3): 17-20.

    Discussion on Common Problems in the Prediction of Reservoir Water Temperature Simulation

    LI Ke-feng, LIANG Rui-feng, TUO You-cai, DENG Yun, LI Jia

    (State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu 610065, China)

    Concerning the common problems affecting prediction results of reservoir water temperature, the method for the determination of natural water temperature in ungauged basins was put forward based on water-air temperature correlation and quantitative analysis for reservoir topographic correction was conducted. The impact of air temperature lapse rate on long-wave radiation, evaporation and conduction across water-air interface was discussed. The difference of heat flux calculated by some empirical meteorological formulas was compared. Qualitative analysis was conducted by comparing the zero equation turbulence model and double equation turbulence model, as well as investigating the difference between the hydrostatic symposium and complete z-momentum equation. A method for quantitative rationality verification of simulation results was put forward through the analysis on heat balance and correlation between reservoir capacity and monthly runoff.

    reservoir water temperature; topography; meteorology; mathematical model

    2015-11-22

    國家自然科學(xué)基金(51379136;50979063;51279114)

    李克鋒(1965—),男,湖北石首人,研究員,博士,主要研究方向為生態(tài)與環(huán)境水力學(xué),E-mail:kefengli@scu.edu.cn

    10.14068/j.ceia.2016.03.015

    X820.3;TV697.2

    A

    2095-6444(2016)03-0057-05

    猜你喜歡
    水溫氣溫水庫
    氣溫的神奇“魔法”
    基于FY-3D和FY-4A的氣溫時空融合
    中型水庫的工程建設(shè)與管理探討
    出山店水庫
    白沙水庫
    基于PLC的水溫控制系統(tǒng)設(shè)計
    電子制作(2019年7期)2019-04-25 13:18:10
    基于DS18B20水溫控制系統(tǒng)設(shè)計
    電子制作(2018年17期)2018-09-28 01:56:38
    出山店水庫
    與氣溫成反比的東西
    大聯(lián)大:最先感知“水溫”變化
    99riav亚洲国产免费| 精品第一国产精品| 国产亚洲欧美在线一区二区| 大型黄色视频在线免费观看| 午夜免费鲁丝| 亚洲国产精品一区二区三区在线| 好男人电影高清在线观看| 黄色怎么调成土黄色| 一区在线观看完整版| 国产欧美亚洲国产| 丝瓜视频免费看黄片| 怎么达到女性高潮| 国产精品av久久久久免费| 精品亚洲成a人片在线观看| 伦理电影免费视频| 满18在线观看网站| 日韩精品免费视频一区二区三区| 美女视频免费永久观看网站| 国产高清视频在线播放一区| 岛国在线观看网站| 激情在线观看视频在线高清 | bbb黄色大片| 91国产中文字幕| 人成视频在线观看免费观看| 亚洲av日韩精品久久久久久密| 日韩欧美三级三区| 99re6热这里在线精品视频| 80岁老熟妇乱子伦牲交| 久久狼人影院| 精品视频人人做人人爽| 国产精品香港三级国产av潘金莲| 日本欧美视频一区| 日韩精品免费视频一区二区三区| 国产高清激情床上av| 欧美黑人精品巨大| 久热这里只有精品99| 国产91精品成人一区二区三区 | 成人国语在线视频| 国产成人啪精品午夜网站| 免费在线观看日本一区| 老司机靠b影院| 国产欧美日韩一区二区三| 久久中文字幕一级| 免费不卡黄色视频| 99久久99久久久精品蜜桃| 叶爱在线成人免费视频播放| 午夜老司机福利片| 脱女人内裤的视频| 美女视频免费永久观看网站| 亚洲精品美女久久久久99蜜臀| 免费女性裸体啪啪无遮挡网站| 国产熟女午夜一区二区三区| 免费在线观看日本一区| 91字幕亚洲| 美女高潮喷水抽搐中文字幕| 久久精品亚洲av国产电影网| 我的亚洲天堂| 妹子高潮喷水视频| 亚洲精品美女久久久久99蜜臀| 超色免费av| 欧美另类亚洲清纯唯美| 欧美性长视频在线观看| 亚洲av成人不卡在线观看播放网| 90打野战视频偷拍视频| 视频区图区小说| 怎么达到女性高潮| 999久久久国产精品视频| kizo精华| 国产欧美日韩一区二区精品| xxxhd国产人妻xxx| 老汉色∧v一级毛片| 他把我摸到了高潮在线观看 | 女同久久另类99精品国产91| 亚洲精品久久成人aⅴ小说| 亚洲性夜色夜夜综合| 午夜精品久久久久久毛片777| 制服人妻中文乱码| 国产精品久久久久久精品古装| 国产免费av片在线观看野外av| 天堂中文最新版在线下载| 十分钟在线观看高清视频www| 天堂8中文在线网| 精品国产国语对白av| 人人妻人人爽人人添夜夜欢视频| 新久久久久国产一级毛片| 国产精品久久久av美女十八| avwww免费| 亚洲精品美女久久久久99蜜臀| 99精国产麻豆久久婷婷| 在线天堂中文资源库| 国产在线观看jvid| 国产成人免费观看mmmm| 性少妇av在线| 国产成人欧美在线观看 | 中亚洲国语对白在线视频| 中国美女看黄片| 国内毛片毛片毛片毛片毛片| 亚洲七黄色美女视频| 电影成人av| 成人国产一区最新在线观看| 妹子高潮喷水视频| 乱人伦中国视频| 悠悠久久av| 久久久久久免费高清国产稀缺| 女性被躁到高潮视频| 麻豆av在线久日| 亚洲成av片中文字幕在线观看| 日韩免费av在线播放| 十八禁高潮呻吟视频| 亚洲av成人不卡在线观看播放网| 亚洲av成人不卡在线观看播放网| 在线观看免费视频网站a站| 国产又色又爽无遮挡免费看| 国产淫语在线视频| 精品乱码久久久久久99久播| 亚洲男人天堂网一区| 美女高潮喷水抽搐中文字幕| 人人妻人人澡人人爽人人夜夜| 成年人免费黄色播放视频| 精品人妻熟女毛片av久久网站| 99国产综合亚洲精品| 天天添夜夜摸| 在线观看免费日韩欧美大片| 免费在线观看日本一区| 热re99久久精品国产66热6| 国产精品国产av在线观看| 深夜精品福利| 黄网站色视频无遮挡免费观看| 在线观看www视频免费| 亚洲视频免费观看视频| 久久亚洲真实| 最新美女视频免费是黄的| 美国免费a级毛片| videos熟女内射| 国产欧美日韩精品亚洲av| 精品福利永久在线观看| 老司机福利观看| 我要看黄色一级片免费的| 18禁观看日本| 中文字幕人妻丝袜制服| 精品久久久久久电影网| 久久精品亚洲精品国产色婷小说| 丝瓜视频免费看黄片| 国产91精品成人一区二区三区 | 久久青草综合色| 飞空精品影院首页| 99国产精品免费福利视频| 国产熟女午夜一区二区三区| 国产99久久九九免费精品| 999久久久精品免费观看国产| 久久精品熟女亚洲av麻豆精品| 国产伦理片在线播放av一区| 国产av又大| 国产单亲对白刺激| 91九色精品人成在线观看| 国产精品久久久久久精品古装| 国产高清视频在线播放一区| 中文字幕人妻丝袜一区二区| 如日韩欧美国产精品一区二区三区| 亚洲国产av新网站| 黑人操中国人逼视频| 欧美国产精品va在线观看不卡| 亚洲色图av天堂| 免费看a级黄色片| 国产亚洲欧美在线一区二区| 麻豆av在线久日| 一区在线观看完整版| 精品第一国产精品| 亚洲欧美日韩另类电影网站| 欧美老熟妇乱子伦牲交| 这个男人来自地球电影免费观看| 国产成人免费无遮挡视频| 亚洲欧美激情在线| 90打野战视频偷拍视频| 色婷婷av一区二区三区视频| a级毛片黄视频| 精品亚洲成a人片在线观看| 久久精品国产亚洲av高清一级| 国产精品秋霞免费鲁丝片| 在线观看免费高清a一片| 妹子高潮喷水视频| 欧美午夜高清在线| 无人区码免费观看不卡 | 精品乱码久久久久久99久播| 在线观看免费日韩欧美大片| 啦啦啦在线免费观看视频4| 又紧又爽又黄一区二区| 美女高潮到喷水免费观看| 99精国产麻豆久久婷婷| 视频区图区小说| 日日夜夜操网爽| 自线自在国产av| 国产在线视频一区二区| 久久久精品国产亚洲av高清涩受| 一区二区三区国产精品乱码| 亚洲欧美色中文字幕在线| 精品久久久久久久毛片微露脸| 十分钟在线观看高清视频www| 精品少妇内射三级| 日韩欧美国产一区二区入口| 最近最新中文字幕大全免费视频| 成人18禁在线播放| 亚洲av成人不卡在线观看播放网| 日韩视频一区二区在线观看| 老司机福利观看| 国产精品久久久久成人av| 久久亚洲精品不卡| bbb黄色大片| 久久精品亚洲熟妇少妇任你| 精品亚洲乱码少妇综合久久| 美女高潮到喷水免费观看| 亚洲精品美女久久久久99蜜臀| 另类精品久久| 久久精品熟女亚洲av麻豆精品| 国产午夜精品久久久久久| 国产亚洲精品久久久久5区| 视频区图区小说| 精品国产一区二区久久| 黑人欧美特级aaaaaa片| 王馨瑶露胸无遮挡在线观看| 久久久久久久大尺度免费视频| 国产精品久久久久久精品古装| 国产日韩欧美亚洲二区| 亚洲熟妇熟女久久| 国产伦人伦偷精品视频| 欧美日韩成人在线一区二区| 国产精品免费大片| 窝窝影院91人妻| 一本大道久久a久久精品| 乱人伦中国视频| 欧美在线黄色| 极品少妇高潮喷水抽搐| 日本精品一区二区三区蜜桃| 久久人妻熟女aⅴ| 亚洲五月婷婷丁香| 在线播放国产精品三级| 一级毛片精品| 黄色毛片三级朝国网站| 亚洲精品在线美女| 亚洲精品在线观看二区| 考比视频在线观看| 久久久久久久久久久久大奶| 日本一区二区免费在线视频| 无人区码免费观看不卡 | 国产色视频综合| 亚洲成a人片在线一区二区| 啦啦啦中文免费视频观看日本| 国产精品国产高清国产av | 国产免费视频播放在线视频| 亚洲精品成人av观看孕妇| 久久婷婷成人综合色麻豆| 亚洲一区中文字幕在线| 老熟女久久久| 在线永久观看黄色视频| 12—13女人毛片做爰片一| 美女国产高潮福利片在线看| 老司机在亚洲福利影院| 免费在线观看影片大全网站| av天堂在线播放| 国产精品.久久久| 啪啪无遮挡十八禁网站| 电影成人av| 久久久国产一区二区| 欧美久久黑人一区二区| 日日摸夜夜添夜夜添小说| 欧美 日韩 精品 国产| 国产野战对白在线观看| 美女高潮喷水抽搐中文字幕| 大香蕉久久成人网| 日韩一卡2卡3卡4卡2021年| 夫妻午夜视频| 国产精品国产高清国产av | 精品亚洲成国产av| 国产成人系列免费观看| 丝袜人妻中文字幕| 9色porny在线观看| 无遮挡黄片免费观看| 国产精品成人在线| 亚洲精品av麻豆狂野| 免费少妇av软件| 久久免费观看电影| 女人爽到高潮嗷嗷叫在线视频| 国产精品麻豆人妻色哟哟久久| 91国产中文字幕| 久久精品国产99精品国产亚洲性色 | 极品教师在线免费播放| 精品欧美一区二区三区在线| 黑人巨大精品欧美一区二区蜜桃| 欧美在线一区亚洲| 久久99热这里只频精品6学生| 嫩草影视91久久| 亚洲av国产av综合av卡| 精品卡一卡二卡四卡免费| 一区二区日韩欧美中文字幕| 别揉我奶头~嗯~啊~动态视频| 国产精品熟女久久久久浪| 国产主播在线观看一区二区| 亚洲伊人久久精品综合| 成年版毛片免费区| 俄罗斯特黄特色一大片| 国产亚洲午夜精品一区二区久久| 精品亚洲成a人片在线观看| 亚洲黑人精品在线| 丰满少妇做爰视频| 国产高清视频在线播放一区| 男女边摸边吃奶| 在线播放国产精品三级| 国产欧美日韩一区二区三区在线| 久久久久久久久久久久大奶| 精品国产乱子伦一区二区三区| 精品一区二区三卡| 两性夫妻黄色片| 亚洲自偷自拍图片 自拍| 国产高清激情床上av| 国内毛片毛片毛片毛片毛片| 国产av国产精品国产| 免费在线观看影片大全网站| 国产福利在线免费观看视频| 国产高清国产精品国产三级| 国产高清激情床上av| 久久久久国内视频| 成年版毛片免费区| 日本精品一区二区三区蜜桃| 国产亚洲午夜精品一区二区久久| 999久久久国产精品视频| 啪啪无遮挡十八禁网站| av国产精品久久久久影院| 精品人妻在线不人妻| 亚洲国产av新网站| 少妇 在线观看| 国产伦人伦偷精品视频| 久久99热这里只频精品6学生| 国产不卡av网站在线观看| 国产极品粉嫩免费观看在线| 欧美 日韩 精品 国产| 黑丝袜美女国产一区| 精品第一国产精品| 一个人免费看片子| 一区二区av电影网| 九色亚洲精品在线播放| 性色av乱码一区二区三区2| 色尼玛亚洲综合影院| 极品少妇高潮喷水抽搐| 视频在线观看一区二区三区| av不卡在线播放| 日韩欧美一区视频在线观看| 久久久久视频综合| 热re99久久精品国产66热6| 国产精品久久久久久人妻精品电影 | 纯流量卡能插随身wifi吗| 露出奶头的视频| 国产一区二区 视频在线| 日韩制服丝袜自拍偷拍| 最近最新中文字幕大全免费视频| 日韩免费高清中文字幕av| 亚洲精品国产精品久久久不卡| 国产97色在线日韩免费| 婷婷成人精品国产| 亚洲av第一区精品v没综合| 99久久99久久久精品蜜桃| avwww免费| 成人国语在线视频| 亚洲熟女毛片儿| 久久中文字幕人妻熟女| 成人影院久久| 日韩精品免费视频一区二区三区| 成年人黄色毛片网站| 下体分泌物呈黄色| 精品国产一区二区三区久久久樱花| 人人澡人人妻人| 久久久久久久国产电影| 日韩大码丰满熟妇| 亚洲成人手机| 国产在线观看jvid| 国产97色在线日韩免费| 在线观看www视频免费| 中亚洲国语对白在线视频| 亚洲精品av麻豆狂野| 久久av网站| 18禁黄网站禁片午夜丰满| 色老头精品视频在线观看| 日韩视频在线欧美| 欧美黑人精品巨大| 乱人伦中国视频| 国产亚洲欧美在线一区二区| av超薄肉色丝袜交足视频| 99国产精品99久久久久| 中文字幕制服av| 国产在线免费精品| 日韩大码丰满熟妇| 国产精品久久久人人做人人爽| 亚洲成av片中文字幕在线观看| 黄频高清免费视频| www.精华液| 五月天丁香电影| 十八禁高潮呻吟视频| 男女下面插进去视频免费观看| 久久性视频一级片| 久久人妻av系列| 国产视频一区二区在线看| 亚洲精品国产精品久久久不卡| 日本一区二区免费在线视频| 美女福利国产在线| 国产淫语在线视频| 免费看十八禁软件| 亚洲国产欧美网| av线在线观看网站| 久久久水蜜桃国产精品网| 精品一区二区三区四区五区乱码| 久久99热这里只频精品6学生| 国产一区二区三区视频了| 欧美激情高清一区二区三区| 女警被强在线播放| 欧美性长视频在线观看| 国产亚洲精品久久久久5区| 精品第一国产精品| 亚洲伊人久久精品综合| 在线观看人妻少妇| 欧美日韩福利视频一区二区| 每晚都被弄得嗷嗷叫到高潮| 久久人妻福利社区极品人妻图片| 脱女人内裤的视频| 在线天堂中文资源库| 久久久久久亚洲精品国产蜜桃av| 午夜福利视频在线观看免费| 午夜精品久久久久久毛片777| 美女视频免费永久观看网站| 十八禁人妻一区二区| 在线观看66精品国产| 性高湖久久久久久久久免费观看| 国产av国产精品国产| 成年动漫av网址| av网站在线播放免费| 另类精品久久| 午夜精品久久久久久毛片777| 国产高清视频在线播放一区| 色94色欧美一区二区| 人妻久久中文字幕网| 一区二区三区乱码不卡18| 在线观看一区二区三区激情| 熟女少妇亚洲综合色aaa.| 丰满人妻熟妇乱又伦精品不卡| 巨乳人妻的诱惑在线观看| 九色亚洲精品在线播放| 精品乱码久久久久久99久播| 一本色道久久久久久精品综合| 久久久久久人人人人人| 中文字幕人妻丝袜制服| 啦啦啦中文免费视频观看日本| videosex国产| 国产av国产精品国产| 亚洲av日韩精品久久久久久密| 男女高潮啪啪啪动态图| 十八禁网站免费在线| 日日爽夜夜爽网站| 亚洲精品一二三| 天堂8中文在线网| 欧美黑人精品巨大| 亚洲欧美日韩高清在线视频 | av超薄肉色丝袜交足视频| av电影中文网址| 69av精品久久久久久 | 人成视频在线观看免费观看| 一边摸一边做爽爽视频免费| 国产成人精品无人区| 国产亚洲av高清不卡| 丝袜美腿诱惑在线| 手机成人av网站| 亚洲性夜色夜夜综合| 老司机深夜福利视频在线观看| 99国产精品免费福利视频| 午夜日韩欧美国产| 免费观看av网站的网址| 一本综合久久免费| 久久久久久久久久久久大奶| a级毛片黄视频| 最新的欧美精品一区二区| 夜夜骑夜夜射夜夜干| 久久久久久久精品吃奶| 午夜激情av网站| 国产免费现黄频在线看| 91精品国产国语对白视频| 99在线人妻在线中文字幕 | 日本a在线网址| 91老司机精品| 一进一出好大好爽视频| 老司机午夜福利在线观看视频 | 18禁美女被吸乳视频| 亚洲av成人一区二区三| 纵有疾风起免费观看全集完整版| 老司机亚洲免费影院| 男女边摸边吃奶| 手机成人av网站| 如日韩欧美国产精品一区二区三区| 成人精品一区二区免费| 成在线人永久免费视频| 国产麻豆69| 日日摸夜夜添夜夜添小说| 国产av一区二区精品久久| 天堂8中文在线网| 99国产精品99久久久久| 亚洲成av片中文字幕在线观看| 桃花免费在线播放| 日韩精品免费视频一区二区三区| 国产男女超爽视频在线观看| 一本大道久久a久久精品| 日本欧美视频一区| svipshipincom国产片| av欧美777| 99久久国产精品久久久| 亚洲第一av免费看| 老司机午夜十八禁免费视频| av在线播放免费不卡| 日本撒尿小便嘘嘘汇集6| 国产日韩欧美在线精品| 国产麻豆69| 国产高清激情床上av| 亚洲精品自拍成人| 伊人久久大香线蕉亚洲五| 亚洲人成电影免费在线| 好男人电影高清在线观看| 人人妻人人澡人人看| 日韩免费av在线播放| 久久久国产成人免费| 12—13女人毛片做爰片一| 国产av又大| 蜜桃国产av成人99| 成年人黄色毛片网站| 中文亚洲av片在线观看爽 | 变态另类成人亚洲欧美熟女 | 精品福利观看| 亚洲国产欧美一区二区综合| 狠狠精品人妻久久久久久综合| 久9热在线精品视频| 丁香六月天网| 日本vs欧美在线观看视频| 五月天丁香电影| 午夜福利,免费看| 亚洲精品国产一区二区精华液| 午夜91福利影院| 亚洲欧美精品综合一区二区三区| 久久亚洲真实| 制服人妻中文乱码| 亚洲久久久国产精品| 久久精品91无色码中文字幕| 少妇 在线观看| 免费人妻精品一区二区三区视频| 18禁观看日本| 母亲3免费完整高清在线观看| 狠狠精品人妻久久久久久综合| 黄片大片在线免费观看| 好男人电影高清在线观看| 老司机影院毛片| 国产精品98久久久久久宅男小说| 国产成人精品无人区| 九色亚洲精品在线播放| 亚洲精品中文字幕一二三四区 | 91精品三级在线观看| 丰满迷人的少妇在线观看| 久久国产精品大桥未久av| 最黄视频免费看| 91成人精品电影| 十八禁人妻一区二区| 久久香蕉激情| 日本欧美视频一区| 欧美久久黑人一区二区| 久久天堂一区二区三区四区| 国产在线精品亚洲第一网站| 亚洲美女黄片视频| 日韩一卡2卡3卡4卡2021年| 热re99久久精品国产66热6| 日韩欧美免费精品| 色综合婷婷激情| 午夜福利在线免费观看网站| 一区二区三区国产精品乱码| 久久国产精品大桥未久av| 麻豆乱淫一区二区| 精品国产一区二区三区久久久樱花| 757午夜福利合集在线观看| 国产精品av久久久久免费| 午夜免费鲁丝| 国产精品.久久久| 中文欧美无线码| av又黄又爽大尺度在线免费看| 国产日韩欧美在线精品| 夜夜夜夜夜久久久久| 波多野结衣av一区二区av| 99精品久久久久人妻精品| 久久精品亚洲精品国产色婷小说| 黄色视频不卡| www.精华液| 涩涩av久久男人的天堂| 欧美日韩亚洲综合一区二区三区_| 99国产精品99久久久久| 国产一卡二卡三卡精品| 免费av中文字幕在线| 亚洲av成人不卡在线观看播放网| 两人在一起打扑克的视频| 久久ye,这里只有精品| 18禁观看日本| 欧美亚洲 丝袜 人妻 在线| 亚洲专区字幕在线| a级片在线免费高清观看视频| 欧美成狂野欧美在线观看| 国产精品免费大片| 久久精品国产亚洲av香蕉五月 | 久久精品91无色码中文字幕| 亚洲欧美日韩另类电影网站| 男女免费视频国产| 国产精品秋霞免费鲁丝片| 日本黄色视频三级网站网址 | 一边摸一边抽搐一进一出视频| 亚洲中文av在线| 久久人妻av系列| 国产一区二区激情短视频|