• <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)大:最先感知“水溫”變化
    亚洲黑人精品在线| 国产在线一区二区三区精| 一个人免费在线观看的高清视频| 成人精品一区二区免费| 久久午夜综合久久蜜桃| 99在线人妻在线中文字幕 | 国产高清国产精品国产三级| 美女高潮喷水抽搐中文字幕| 亚洲情色 制服丝袜| 操出白浆在线播放| 人成视频在线观看免费观看| 国产伦理片在线播放av一区| 老汉色∧v一级毛片| 欧美 亚洲 国产 日韩一| 亚洲精品国产色婷婷电影| 久久久久久亚洲精品国产蜜桃av| 久久国产精品男人的天堂亚洲| 成人18禁在线播放| av天堂久久9| 国产成人啪精品午夜网站| www.精华液| 日本撒尿小便嘘嘘汇集6| 成人国语在线视频| 狂野欧美激情性xxxx| 亚洲人成电影观看| 热99re8久久精品国产| av福利片在线| 老司机福利观看| 欧美国产精品一级二级三级| 成人18禁在线播放| 天堂8中文在线网| 正在播放国产对白刺激| 91精品国产国语对白视频| 亚洲精品国产一区二区精华液| 久久久久久久久免费视频了| svipshipincom国产片| 午夜老司机福利片| 成人国产av品久久久| 肉色欧美久久久久久久蜜桃| 美女国产高潮福利片在线看| 十八禁网站免费在线| 在线观看www视频免费| 亚洲国产av新网站| 亚洲av国产av综合av卡| 老汉色av国产亚洲站长工具| 一区福利在线观看| 侵犯人妻中文字幕一二三四区| netflix在线观看网站| 中文字幕高清在线视频| 午夜视频精品福利| 一级片免费观看大全| 老熟女久久久| 我的亚洲天堂| 成人18禁在线播放| 日本av手机在线免费观看| 99国产精品99久久久久| 亚洲精品国产色婷婷电影| 亚洲五月婷婷丁香| 国产免费现黄频在线看| 国产欧美日韩精品亚洲av| 午夜福利一区二区在线看| 国产精品久久久av美女十八| 日韩视频在线欧美| 亚洲av成人不卡在线观看播放网| 黑人巨大精品欧美一区二区mp4| 精品少妇久久久久久888优播| 老司机福利观看| 亚洲国产av新网站| 国产色视频综合| 俄罗斯特黄特色一大片| 日韩 欧美 亚洲 中文字幕| 欧美成人免费av一区二区三区 | 极品教师在线免费播放| 国产免费av片在线观看野外av| av线在线观看网站| 一级毛片电影观看| 欧美人与性动交α欧美精品济南到| 亚洲第一av免费看| 亚洲欧美一区二区三区久久| 男女高潮啪啪啪动态图| 国产一区二区在线观看av| 国产麻豆69| 人人妻人人澡人人爽人人夜夜| 欧美激情 高清一区二区三区| 成人av一区二区三区在线看| 日日夜夜操网爽| 国产人伦9x9x在线观看| www.熟女人妻精品国产| 18在线观看网站| 一级片免费观看大全| 国产免费福利视频在线观看| 国产精品一区二区在线观看99| 国产一区二区三区综合在线观看| 亚洲久久久国产精品| 午夜福利,免费看| 搡老岳熟女国产| 99九九在线精品视频| 国产高清videossex| 亚洲av国产av综合av卡| 欧美黄色片欧美黄色片| 丁香六月天网| 中文亚洲av片在线观看爽 | 久久久国产精品麻豆| 免费少妇av软件| 久久久精品区二区三区| 美国免费a级毛片| 制服人妻中文乱码| 性少妇av在线| 高清毛片免费观看视频网站 | 69精品国产乱码久久久| 久久久国产精品麻豆| 国产精品久久久人人做人人爽| 欧美 日韩 精品 国产| 美女国产高潮福利片在线看| 久久久国产一区二区| 人成视频在线观看免费观看| 老司机靠b影院| 黄色 视频免费看| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品久久久人人做人人爽| 欧美午夜高清在线| 91精品三级在线观看| 亚洲熟妇熟女久久| 国产欧美日韩精品亚洲av| 欧美日韩中文字幕国产精品一区二区三区 | 国产成人免费观看mmmm| 狠狠狠狠99中文字幕| 人妻一区二区av| 少妇被粗大的猛进出69影院| 久久久欧美国产精品| 丰满饥渴人妻一区二区三| 亚洲自偷自拍图片 自拍| av电影中文网址| 一边摸一边抽搐一进一出视频| 成年人午夜在线观看视频| 欧美日韩av久久| 午夜两性在线视频| 欧美乱码精品一区二区三区| av电影中文网址| 亚洲av成人一区二区三| 王馨瑶露胸无遮挡在线观看| 日韩免费高清中文字幕av| 欧美人与性动交α欧美软件| 精品国产国语对白av| 欧美日韩一级在线毛片| 国产男女超爽视频在线观看| 另类亚洲欧美激情| 国产精品国产高清国产av | 黄色丝袜av网址大全| 99精品久久久久人妻精品| 巨乳人妻的诱惑在线观看| 国产精品一区二区精品视频观看| 高清视频免费观看一区二区| 国产无遮挡羞羞视频在线观看| 麻豆av在线久日| videos熟女内射| 亚洲成人免费av在线播放| 老熟女久久久| 老熟妇仑乱视频hdxx| 国产成人欧美| 自线自在国产av| 人人妻,人人澡人人爽秒播| 水蜜桃什么品种好| 亚洲一区中文字幕在线| 手机成人av网站| 国产老妇伦熟女老妇高清| 中文亚洲av片在线观看爽 | 日日夜夜操网爽| 91国产中文字幕| 亚洲成人免费电影在线观看| 大片电影免费在线观看免费| 亚洲国产欧美一区二区综合| 国产免费av片在线观看野外av| 久久精品熟女亚洲av麻豆精品| 国产xxxxx性猛交| 色在线成人网| tube8黄色片| 国产日韩欧美在线精品| 免费看a级黄色片| 99国产极品粉嫩在线观看| 免费日韩欧美在线观看| 久久 成人 亚洲| 色视频在线一区二区三区| 9色porny在线观看| 色尼玛亚洲综合影院| 人妻久久中文字幕网| 别揉我奶头~嗯~啊~动态视频| 视频在线观看一区二区三区| 国产淫语在线视频| 啦啦啦免费观看视频1| 两个人免费观看高清视频| 91精品三级在线观看| 日韩欧美免费精品| 免费看十八禁软件| 欧美亚洲 丝袜 人妻 在线| 亚洲人成伊人成综合网2020| 精品一区二区三区四区五区乱码| 久久免费观看电影| 欧美日本中文国产一区发布| 高清视频免费观看一区二区| tube8黄色片| 欧美在线黄色| 亚洲成人国产一区在线观看| www.自偷自拍.com| www.自偷自拍.com| 精品视频人人做人人爽| a级片在线免费高清观看视频| 美女高潮到喷水免费观看| 亚洲av欧美aⅴ国产| 国内毛片毛片毛片毛片毛片| 91av网站免费观看| 中文欧美无线码| 国产熟女午夜一区二区三区| 国产精品美女特级片免费视频播放器 | 精品亚洲乱码少妇综合久久| 丝袜在线中文字幕| 超碰成人久久| 一级毛片女人18水好多| 啦啦啦视频在线资源免费观看| 久久精品国产99精品国产亚洲性色 | 精品久久蜜臀av无| 国产亚洲午夜精品一区二区久久| 我的亚洲天堂| 亚洲va日本ⅴa欧美va伊人久久| 亚洲专区国产一区二区| 国产高清视频在线播放一区| 丁香六月天网| 高清视频免费观看一区二区| 亚洲一码二码三码区别大吗| 国产欧美日韩一区二区精品| 亚洲国产中文字幕在线视频| 国产国语露脸激情在线看| 操美女的视频在线观看| 男人操女人黄网站| 亚洲熟妇熟女久久| 在线看a的网站| 丝袜美腿诱惑在线| 岛国毛片在线播放| 五月开心婷婷网| 十八禁网站网址无遮挡| 久久国产精品大桥未久av| 久久人妻熟女aⅴ| 午夜福利视频精品| 老司机影院毛片| 国产精品免费一区二区三区在线 | 精品少妇一区二区三区视频日本电影| 在线观看免费视频网站a站| 亚洲av欧美aⅴ国产| 国产精品亚洲av一区麻豆| 国产日韩一区二区三区精品不卡| 久久国产精品影院| 欧美黑人精品巨大| 久久人妻福利社区极品人妻图片| 最黄视频免费看| 国产亚洲午夜精品一区二区久久| 亚洲av片天天在线观看| 国产在视频线精品| 亚洲欧美一区二区三区久久| 国产高清视频在线播放一区| www.999成人在线观看| 国产午夜精品久久久久久| 久久人妻熟女aⅴ| 久久午夜亚洲精品久久| 交换朋友夫妻互换小说| 欧美成狂野欧美在线观看| 亚洲人成伊人成综合网2020| 欧美成人免费av一区二区三区 | 国产av精品麻豆| 国产精品久久久久久人妻精品电影 | 国产一区二区在线观看av| 亚洲午夜理论影院| 国产精品亚洲一级av第二区| 丁香六月天网| 日韩熟女老妇一区二区性免费视频| 夜夜爽天天搞| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品国产色婷婷电影| 久久精品国产综合久久久| 久久久久久人人人人人| 精品视频人人做人人爽| 男女之事视频高清在线观看| 午夜老司机福利片| 2018国产大陆天天弄谢| 最近最新中文字幕大全电影3 | 亚洲一区中文字幕在线| 嫩草影视91久久| 精品一品国产午夜福利视频| 狠狠狠狠99中文字幕| 最近最新免费中文字幕在线| 18禁黄网站禁片午夜丰满| 中文字幕最新亚洲高清| 午夜免费鲁丝| 成年人午夜在线观看视频| 97在线人人人人妻| 涩涩av久久男人的天堂| av不卡在线播放| 亚洲熟女毛片儿| 日韩免费av在线播放| 免费看a级黄色片| 一二三四社区在线视频社区8| 51午夜福利影视在线观看| 国产欧美日韩一区二区精品| 高清av免费在线| 亚洲自偷自拍图片 自拍| 国产日韩欧美视频二区| 涩涩av久久男人的天堂| av有码第一页| 1024视频免费在线观看| 一本—道久久a久久精品蜜桃钙片| 午夜成年电影在线免费观看| 日日摸夜夜添夜夜添小说| 午夜两性在线视频| 成年动漫av网址| 俄罗斯特黄特色一大片| 国产成人av教育| 精品国产乱子伦一区二区三区| 手机成人av网站| 丝袜美腿诱惑在线| 丝袜人妻中文字幕| 亚洲国产av新网站| 18禁裸乳无遮挡动漫免费视频| 制服人妻中文乱码| 国产老妇伦熟女老妇高清| 人人妻,人人澡人人爽秒播| 亚洲成人国产一区在线观看| av网站免费在线观看视频| 99热网站在线观看| 亚洲欧美一区二区三区久久| 欧美成人午夜精品| 一进一出好大好爽视频| 狠狠狠狠99中文字幕| 亚洲成a人片在线一区二区| 成人国语在线视频| 黄色怎么调成土黄色| 亚洲欧美激情在线| 亚洲,欧美精品.| 午夜福利在线观看吧| 亚洲综合色网址| 一区二区三区国产精品乱码| 国产精品电影一区二区三区 | 亚洲一卡2卡3卡4卡5卡精品中文| av网站在线播放免费| 精品卡一卡二卡四卡免费| 超碰97精品在线观看| 久久精品aⅴ一区二区三区四区| 天天添夜夜摸| tocl精华| 久久人妻av系列| 国产精品98久久久久久宅男小说| 国产在线观看jvid| 大陆偷拍与自拍| 99久久国产精品久久久| 日日爽夜夜爽网站| 精品一区二区三区视频在线观看免费 | 真人做人爱边吃奶动态| 每晚都被弄得嗷嗷叫到高潮| 久久天堂一区二区三区四区| 又大又爽又粗| 男女高潮啪啪啪动态图| 久久精品亚洲熟妇少妇任你| 丁香六月天网| 久久中文字幕一级| 日本wwww免费看| 国产亚洲精品一区二区www | av免费在线观看网站| 在线永久观看黄色视频| 在线观看免费日韩欧美大片| 黄色a级毛片大全视频| 老熟妇仑乱视频hdxx| 亚洲精品国产精品久久久不卡| 一夜夜www| 日本一区二区免费在线视频| av国产精品久久久久影院| 中文亚洲av片在线观看爽 | av天堂在线播放| 国内毛片毛片毛片毛片毛片| 日韩大片免费观看网站| 一本综合久久免费| 波多野结衣av一区二区av| 国产伦人伦偷精品视频| 国产精品久久电影中文字幕 | 美女主播在线视频| 精品高清国产在线一区| 男女午夜视频在线观看| 国产三级黄色录像| 两人在一起打扑克的视频| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品 欧美亚洲| 天堂8中文在线网| 女性被躁到高潮视频| 伊人久久大香线蕉亚洲五| 80岁老熟妇乱子伦牲交| 国产精品一区二区在线观看99| 1024视频免费在线观看| 日韩大片免费观看网站| 亚洲精品久久午夜乱码| 欧美激情久久久久久爽电影 | 飞空精品影院首页| 韩国精品一区二区三区| 亚洲成人免费av在线播放| 亚洲欧美一区二区三区黑人| 欧美成狂野欧美在线观看| avwww免费| 老司机影院毛片| 欧美成狂野欧美在线观看| 精品少妇黑人巨大在线播放| 嫁个100分男人电影在线观看| 在线观看免费高清a一片| 18禁观看日本| 露出奶头的视频| 成年女人毛片免费观看观看9 | 国产亚洲欧美在线一区二区| 欧美黄色淫秽网站| 久久精品亚洲精品国产色婷小说| 亚洲成人免费av在线播放| 欧美精品人与动牲交sv欧美| 国产高清videossex| 日韩中文字幕欧美一区二区| 美女高潮到喷水免费观看| 免费不卡黄色视频| 国产精品99久久99久久久不卡| 国产91精品成人一区二区三区 | 精品国内亚洲2022精品成人 | 国产成人一区二区三区免费视频网站| 国产成人av激情在线播放| 少妇粗大呻吟视频| 欧美黄色片欧美黄色片| 久久九九热精品免费| 国产高清videossex| 老熟妇仑乱视频hdxx| 国产一区二区三区在线臀色熟女 | 日韩免费高清中文字幕av| 欧美av亚洲av综合av国产av| 51午夜福利影视在线观看| 首页视频小说图片口味搜索| av网站免费在线观看视频| 国产高清videossex| 少妇粗大呻吟视频| 最近最新中文字幕大全免费视频| 伦理电影免费视频| 国产高清视频在线播放一区| 久久人妻福利社区极品人妻图片| 热99re8久久精品国产| 欧美日韩国产mv在线观看视频| 午夜福利欧美成人| 亚洲精品自拍成人| 少妇猛男粗大的猛烈进出视频| 久久国产精品男人的天堂亚洲| 精品国产乱码久久久久久小说| 欧美激情久久久久久爽电影 | 久久 成人 亚洲| 成人精品一区二区免费| 妹子高潮喷水视频| 精品乱码久久久久久99久播| 亚洲天堂av无毛| 亚洲自偷自拍图片 自拍| 久久久精品免费免费高清| 欧美日韩成人在线一区二区| 国精品久久久久久国模美| 国产成人系列免费观看| 丝袜人妻中文字幕| 热re99久久国产66热| 成年动漫av网址| 午夜精品久久久久久毛片777| av不卡在线播放| 久久中文看片网| 日本vs欧美在线观看视频| 国产精品1区2区在线观看. | 我要看黄色一级片免费的| 亚洲国产看品久久| 色精品久久人妻99蜜桃| 啦啦啦视频在线资源免费观看| 亚洲精品在线观看二区| 亚洲黑人精品在线| 精品国产乱子伦一区二区三区| 人人妻,人人澡人人爽秒播| 丝袜人妻中文字幕| 国产免费视频播放在线视频| 国产淫语在线视频| av片东京热男人的天堂| 国产一区二区三区视频了| 久久免费观看电影| 香蕉丝袜av| 老司机福利观看| 亚洲精品在线观看二区| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲专区国产一区二区| 国产成人系列免费观看| 在线永久观看黄色视频| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美 亚洲 国产 日韩一| 精品国产一区二区三区四区第35| 免费人妻精品一区二区三区视频| 国产成人免费无遮挡视频| videos熟女内射| 欧美日韩成人在线一区二区| 亚洲欧美一区二区三区久久| 1024香蕉在线观看| 首页视频小说图片口味搜索| 亚洲欧美日韩另类电影网站| 成在线人永久免费视频| 我的亚洲天堂| avwww免费| 欧美日韩一级在线毛片| 岛国在线观看网站| 99精国产麻豆久久婷婷| 久久午夜综合久久蜜桃| 日韩成人在线观看一区二区三区| 香蕉久久夜色| 一级a爱视频在线免费观看| 欧美黑人精品巨大| 无人区码免费观看不卡 | 国产精品电影一区二区三区 | 久久久久久人人人人人| 欧美激情高清一区二区三区| 久久人妻av系列| 91成人精品电影| 叶爱在线成人免费视频播放| 国产一区二区 视频在线| 免费黄频网站在线观看国产| 欧美日韩成人在线一区二区| 亚洲专区字幕在线| 99国产精品免费福利视频| 母亲3免费完整高清在线观看| 我的亚洲天堂| 午夜福利欧美成人| 菩萨蛮人人尽说江南好唐韦庄| 久久久久网色| 成人免费观看视频高清| 免费av中文字幕在线| 好男人电影高清在线观看| 亚洲一区二区三区欧美精品| 午夜福利视频精品| a级片在线免费高清观看视频| 一进一出抽搐动态| 欧美成人午夜精品| 国产一区二区三区在线臀色熟女 | 美女高潮喷水抽搐中文字幕| 男女高潮啪啪啪动态图| 久久亚洲精品不卡| 日韩三级视频一区二区三区| 亚洲熟女精品中文字幕| 国产精品98久久久久久宅男小说| 国产精品久久久久久精品电影小说| 丝袜美腿诱惑在线| 午夜日韩欧美国产| 久久中文字幕人妻熟女| 99九九在线精品视频| 成在线人永久免费视频| 王馨瑶露胸无遮挡在线观看| 色在线成人网| 一本色道久久久久久精品综合| 欧美精品啪啪一区二区三区| 亚洲成人免费av在线播放| 国产一区二区在线观看av| 国产日韩欧美亚洲二区| 国产成人免费观看mmmm| 丁香六月天网| 亚洲综合色网址| 国产av一区二区精品久久| 最新美女视频免费是黄的| 一二三四社区在线视频社区8| 中文字幕av电影在线播放| 亚洲精品美女久久久久99蜜臀| 国产免费视频播放在线视频| 成年人免费黄色播放视频| 亚洲国产欧美网| 久久精品国产亚洲av高清一级| 啦啦啦视频在线资源免费观看| 国产精品.久久久| 一本久久精品| 亚洲精品中文字幕一二三四区 | 亚洲精品国产区一区二| 18在线观看网站| 亚洲va日本ⅴa欧美va伊人久久| 一夜夜www| 老司机在亚洲福利影院| 国产精品自产拍在线观看55亚洲 | 丝袜美足系列| 99精国产麻豆久久婷婷| 老司机深夜福利视频在线观看| 精品国产一区二区三区久久久樱花| 五月天丁香电影| 亚洲人成电影观看| 操美女的视频在线观看| 免费黄频网站在线观看国产| 国产成人欧美在线观看 | 国产精品香港三级国产av潘金莲| 高清av免费在线| 精品一区二区三卡| 成人18禁在线播放| 欧美黄色片欧美黄色片| 国产1区2区3区精品| 精品国产一区二区三区久久久樱花| 久久精品亚洲熟妇少妇任你| 国产精品久久久久久人妻精品电影 | 精品卡一卡二卡四卡免费| 90打野战视频偷拍视频| 一本久久精品| 视频在线观看一区二区三区| 成人特级黄色片久久久久久久 | 免费在线观看日本一区| 五月开心婷婷网| 国产精品98久久久久久宅男小说| 极品人妻少妇av视频| 大型黄色视频在线免费观看| 国产一区二区在线观看av| 日韩精品免费视频一区二区三区| 精品高清国产在线一区| 嫁个100分男人电影在线观看| 国产成人系列免费观看| 在线看a的网站| 侵犯人妻中文字幕一二三四区| 高清毛片免费观看视频网站 | 丝袜美腿诱惑在线|