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

    越南某電廠溫排水?dāng)?shù)值模擬

    2017-09-22 09:56:09李紹武陳漢寶
    水道港口 2017年4期
    關(guān)鍵詞:港池防波堤取水口

    李紹武,韓 飛,陳漢寶

    (1.天津大學(xué) 水利工程仿真與安全國家重點實驗室,天津 300072;2.交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究所 港口水工建筑技術(shù)國家工程實驗室 工程泥沙交通行業(yè)重點實驗室,天津 300456)

    越南某電廠溫排水?dāng)?shù)值模擬

    李紹武1,韓 飛1,陳漢寶2

    (1.天津大學(xué) 水利工程仿真與安全國家重點實驗室,天津 300072;2.交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究所 港口水工建筑技術(shù)國家工程實驗室 工程泥沙交通行業(yè)重點實驗室,天津 300456)

    為了解越南某電廠新建一期三期工程對已投入運(yùn)行的二期工程取水口溫升的影響,采用MIKE 21中的潮流和溫度輸運(yùn)模塊進(jìn)行了潮流和溫度場的數(shù)值模擬。選取冬季和夏季兩種氣象背景,模擬了典型大潮、小潮情況下一、二、三期工程同時運(yùn)行時的溫度場擴(kuò)散過程,得出不同工況下的溫升包絡(luò)圖。模擬結(jié)果表明:一期和三期工程投入運(yùn)營后,二期工程取水口溫升不超過1℃,符合使用條件。針對工程中的透水防波堤,文章提出了杜比公式與MIKE模型相結(jié)合的模擬方法,并根據(jù)計算結(jié)果分析了防波堤的透水性對溫度場擴(kuò)散的影響。結(jié)果表明,防波堤透水性對溫度場擴(kuò)散影響顯著。

    MIKE21;溫排水;透水防波堤;杜比公式

    數(shù)值模擬方法是一種工程中廣泛采用的較為成熟的研究溫排水運(yùn)動過程的常規(guī)方法[1-5],既有效且經(jīng)濟(jì),許多開源程序或商業(yè)軟件中都有淺水和溫排模塊,但無法直接處理透水堤情況。一般斜坡式防波堤由于拋石孔隙率較大,難以阻擋潮流的穿透,而透過堤身的溫排水對港內(nèi)水體的溫度影響相當(dāng)顯著。MIKE軟件是在國內(nèi)外應(yīng)用十分普遍的水動力數(shù)值模擬軟件。本文提出一種將杜比公式與MIKE21模型相結(jié)合的方法,實現(xiàn)了透水防波堤溫度場輸運(yùn)過程的模擬。

    1 模型理論

    1.1淺水波水動力模塊

    質(zhì)量守恒方程

    (1)

    (2)

    (3)

    1.2溫度模塊

    MIKE21 FM的溫度模塊的控制方程為二維熱輸運(yùn)方程[6]

    (4)

    2 工程概況

    圖1 電廠布置及水文測點Fig.1 Plane layout of the power plant and location of observation stations for tide and current

    越南某燃煤電廠工程位于越南南部。電廠包括:一期2×600 MW 二期2×600 MW 和三期2×1 000 MW,其中二期工程已投入運(yùn)行,一期與三期工程正在建設(shè)。電廠循環(huán)水系統(tǒng)采用海水直流循環(huán)冷卻系統(tǒng),排水口設(shè)置在電廠防波堤外側(cè),取水口設(shè)置在港池內(nèi)部,采用港池明渠取水方案,一、三期工程共用取水明渠(圖1)。

    廠址附近海域的潮汐類型屬于不規(guī)則全日潮,全日分潮占主導(dǎo),其間夾雜不規(guī)則半日潮,每月有6~7 d呈半日潮特征。海平面變化在雨季較小,符合海平面的月際變化規(guī)律,最高潮位一般出現(xiàn)在10月,最低潮位出現(xiàn)在8月。由實測流速資料分析可知,該海域流場為從東北到西南方向的單向流,漲潮時流速較大,落潮時流速較小。根據(jù)1981~1991年的水溫統(tǒng)計資料,該地區(qū)屬于熱帶氣候,年內(nèi)各月的平均水溫變化不大,在26.4~28.0℃,多年平均水溫為27.4℃,最高水溫為32.4℃,最低水溫為22.3℃ 。年內(nèi)各月平均氣溫在25.2~28.4℃,多年平均氣溫為26.7℃,最高氣溫為37.2℃,最低氣溫為16.7℃ 。

    3 數(shù)值模型的建立及計算結(jié)果討論

    3.1計算區(qū)域和網(wǎng)格劃分

    2-a計算區(qū)域地形 2-b計算區(qū)域網(wǎng)格剖分圖2 計算區(qū)域地形及網(wǎng)格劃分Fig.2 Numerical domain and meshes

    根據(jù)電廠周邊區(qū)域的岸線特點,結(jié)合電廠取排水工程布置及地形情況,確定模型計算區(qū)域(圖2),面積約為234.9 km2。采用非結(jié)構(gòu)化三角形網(wǎng)格,節(jié)點數(shù)9 811,網(wǎng)格單元數(shù)19 110,并對取水口、排水口及進(jìn)港航道進(jìn)行網(wǎng)格加密,分辨率最小約12 m。

    3.2模型的初始和邊界條件

    潮流場計算的初始條件采用零啟動條件,即全場流速給定為零。岸邊界采用滑移陸地邊界?;顒舆吔绮捎酶蓾窬W(wǎng)格法,干網(wǎng)格、半淹沒網(wǎng)格和濕水網(wǎng)格截斷水深分別為0.005 m、0.05 m、0.1 m,當(dāng)水深小于干水深時視為干網(wǎng)格,水深大于濕水深時視為濕網(wǎng)格,水深處于干水深與濕水深之間時不計算動量方程,只計算連續(xù)方程。水邊界采用潮位和流速過程。

    溫度場計算的初始條件為環(huán)境水體溫度,冬季取30.8℃,夏季取32.4℃。陸邊界采用絕熱條件,對于開邊界,水流流出,邊界溫度由計算給出,水流流入,邊界溫度為環(huán)境水體溫度。排水口溫度與取水口溫度關(guān)聯(lián),為取水溫度加排水溫差,一期工程循環(huán)水流量為61 m3s,二、三期工程循環(huán)水流量均為56 m3s,取、排水溫差為7℃ 。

    3.3模型參數(shù)設(shè)置

    水流和溫排模型參數(shù)主要包括底部糙率和熱交換系數(shù)。海床底部糙率根據(jù)實測水文資料,通過模型調(diào)試確定。本次計算取曼寧系數(shù)為27~50 m13s。水體與大氣的熱量交換主要包括潛熱通量、可感熱通量、短波輻射、長波輻射4個基本物理過程,取決于赤緯角、與太陽的距離、緯度、地外輻射、風(fēng)速、蒸發(fā)量、云量、氣溫、相對濕度等[7]。模型中指定電廠所在時區(qū)為東七區(qū)??紤]電廠最不利條件,夏季和冬季氣溫分別取月平均氣溫37.2℃、35.3℃,相對濕度分別取為88%和83%,其他相關(guān)系數(shù)通過模型率定確定。

    表1 模擬工況Tab.1 Cases of simulation

    3.4模擬方案

    本次模擬根據(jù)二期工程單獨運(yùn)行和一、二、三期工程同時運(yùn)行、冬季和夏季、典型大潮和典型小潮等不同情況共設(shè)置8種不同計算工況(表1),其中二期工程單獨運(yùn)行工況用來驗證模型,一、二、三期工程同時運(yùn)行工況用來預(yù)測一、三期工程投入運(yùn)行之后對二期取水口溫升的影響。

    圖3 WL1大潮潮位驗證圖 圖4 WL1小潮潮位驗證圖 圖5 CS1大潮流速驗證圖Fig.3 Verification of water surface elevation for spring tide at WL1 Fig.4 Verification of water surface elevation for neap tide at WL1 Fig.5 Verification of velocity for spring tide at CS1

    圖6 CS1小潮流速驗證圖 圖7 CS2大潮流速驗證圖 圖8 CS2小潮流速驗證圖Fig.6 Verification of velocity for neap tide at CS1 Fig.7 Verification of velocity for spring tide at CS2 Fig.8 Verification of velocity for neap tide at CS2

    3.5模型驗證

    3.5.1 潮位和流速驗證

    潮位驗證點為WL1,流速驗證點為CS1、CS2(圖1)。通過模型調(diào)試,水位及流速模擬結(jié)果與實測結(jié)果對比分別如圖3~8所示,實測與模擬值吻合良好。

    3.5.2 溫升

    電廠產(chǎn)生的熱水源源不斷地由排水口排入環(huán)境水體,并隨著潮流向周圍輸移,造成周圍水體溫度升高。在相同的動力-熱力強(qiáng)迫條件下,有溫排水與無溫排水模式計算的溫度場之差為增溫場。根據(jù)穩(wěn)定后的增溫場做出的最大溫升包絡(luò)圖,可以得出模型區(qū)域范圍內(nèi)的任何位置可能出現(xiàn)的最大溫升值。8種情況得出的最大溫升包絡(luò)圖如圖9~16所示。

    由圖9~12可以看出,溫排水輸移的情況與模型區(qū)域流場情況相對應(yīng),由于流場以單向流為主,所以溫度場為單方向輸移。防波堤在溫排水的輸移過程中起到較好的阻隔作用,在夏季情況下,較少的溫排水進(jìn)入港池內(nèi)部從而對取水口水溫產(chǎn)生影響,在冬季情況下,從口門處進(jìn)入港池的溫排水較夏季增加,但0.5℃等溫升線仍未到達(dá)二期工程取水口位置,對取水口水溫的影響仍然較小。二期工程的數(shù)值模擬結(jié)果與實際情況基本吻合。

    圖9 大潮夏季二期溫升包絡(luò)圖 圖10 小潮夏季二期溫升包絡(luò)圖Fig.9 Temperature rising envelope of spring tide in summer for PhaseⅡ Fig.10 Temperature rising envelope of neap tide in summer for PhaseⅡ

    圖11 大潮冬季二期溫升包絡(luò)圖 圖12 小潮冬季二期溫升包絡(luò)圖Fig.11 Temperature rising envelope of spring tide in winter for PhaseⅡ Fig.12 Temperature rising envelope of neap tide in winter for PhaseⅡ

    圖13 大潮夏季一二三期溫升包絡(luò)圖 圖14 小潮夏季一二三期溫升包絡(luò)圖Fig.13 Temperature rising envelope of spring tide in summer for PhaseⅠ,Ⅱ and Ⅲ Fig.14 Temperature rising envelope of neap tide in summer for PhaseⅠ,Ⅱ and Ⅲ

    圖15 大潮冬季一二三期溫升包絡(luò)圖 圖16 小潮冬季一二三期溫升包絡(luò)圖Fig.15 Temperature rising envelope of spring tide in winter for PhaseⅠ,Ⅱ and Ⅲ Fig.16 Temperature rising envelope of neap tide in winter for PhaseⅠ,Ⅱ and Ⅲ

    由圖13~16可以看出,當(dāng)一二三期工程同時投入運(yùn)營之后,由于單位時間內(nèi)排放到環(huán)境水體的溫排水大量增加,造成較多的溫排水通過口門進(jìn)入港池內(nèi)部,在夏季情況下,港池內(nèi)部溫升均超過0.5℃,部分區(qū)域達(dá)到1℃,在冬季情況下,港池內(nèi)部大部分溫升超過了1℃,但由于二期工程取水口位于港池最內(nèi)側(cè),1℃等溫升線仍未到達(dá)二期工程取水口位置。

    由圖9~16可以看出,由于典型大、小潮的平均流速比較接近,溫度場的輸移與潮型關(guān)系不大,且由于冬季水溫相對較低,有利于溫排水的輸移,導(dǎo)致較多的溫排水進(jìn)入港池,使港池內(nèi)部溫升增大。

    表2 不同工況下二期工程取水口溫升特征值Tab.2 Numerical results of water temperature rise at intakes of Phase Ⅱ for each case

    不同工況下二期工程取水口的溫升值如表2所示。由表2可以看出,在一二三期工程同時運(yùn)行的情況下,取水口的最大平均溫升為0.7℃,最大溫升為0.9℃,均小于1℃,符合使用要求。

    3.6透水防波堤的模擬

    斜坡式拋石堤有一定透水性,港外溫排水在港內(nèi)外水位差影響下可以向港內(nèi)滲透,從而引起港內(nèi)取水口處水溫升高,影響機(jī)組冷卻效果。MIKE21 FM水流模型中沒有直接考慮結(jié)構(gòu)透水性的方法。為此,提出以下近似處理方法。

    根據(jù)杜比公式[8],滲流速度可表示為

    (5)

    圖17 簡化滲透模型示意圖Fig.17 Sketch of permeable model

    圖17為水平地基上的防波堤示意圖,防波堤外側(cè)液體通過邊界AB滲入堤內(nèi),在堤內(nèi)形成自由表面光(浸潤面)AC,C點為逸出點,ABCD區(qū)域為滲流區(qū)。防波堤滲流采用“分段法”進(jìn)行計算,A′B′GC為外側(cè)段(將上游楔形體ABE用矩形體AEB′A′取代),內(nèi)側(cè)段為GCD。

    等效矩形體的寬度

    (6)

    由杜比公式可以推導(dǎo)得出水平不透水層上防波堤的滲流計算公式

    防波堤外側(cè)

    (7)

    防波堤內(nèi)側(cè)

    (8)

    式中:m1,m2分別為防波堤外側(cè)和內(nèi)側(cè)的邊坡;H1,H2分別為防波堤外側(cè)和內(nèi)側(cè)的水位;hk為逸出點高度;q為滲流單寬流量;k為滲透系數(shù)。

    通過式(6)、式(7)、式(8)三式聯(lián)立即可求出防波堤的單寬滲流量和逸出點高度。

    首先,在模型中將透水堤按不透水考慮,得到堤內(nèi)外的水位和水溫的變化過程;再假設(shè)滲流不引起堤內(nèi)外水位變化,利用式(6)、式(7)計算出單寬滲流量的變化過程。以單寬滲流量作為流量源,以防波堤內(nèi)外測的水溫差作為熱源,加入到模型中進(jìn)行二次計算,即可求出防波堤的透水性對溫度場擴(kuò)散的影響。

    由于防波堤地基和堤身填土材料的差異,滲透系數(shù)也有所不同[9]。本文選取靠近取水口1 100 m長的兩段防波堤(圖1)按透水堤處理,分別取滲流系數(shù)為0、0.2、0.4和 0.6 cms共4種情況進(jìn)行計算,得到不同透水性對二期工程取水口溫升的影響(表3)。

    表3 不同滲透系數(shù)下二期工程取水口溫升特征值Tab.3 Numerical results of water temperature rise at intakes of Phase Ⅱ for different permeability coefficient

    從表3可以看出,防波堤的滲透性對取水口的溫升有顯著影響。拋石堤的堤心石一般為10~100 kg的塊石,其透水系數(shù)目前尚無較準(zhǔn)確的取值方法。本研究參照卵石滲透系數(shù)取0.2 cms。

    由上述可知,對于取水口處于溫排口下游的情況,若防波堤結(jié)構(gòu)選用斜坡式拋石堤結(jié)構(gòu),對于兼做隔熱和防波雙重作用的斜坡式拋石堤,應(yīng)考慮采取適當(dāng)?shù)姆罎B措施,以免影響取水口取水的水質(zhì)。

    4 結(jié)論

    本文采用MIKE21FM模型對越南某電廠的潮流場和溫度場進(jìn)行了數(shù)值模擬。計算結(jié)果表明,在增加一、三期工程之后,二期工程取水口最大溫升值均未達(dá)到1℃,符合工程要求。在考慮防波堤的滲透性時,本文提出了杜比公式與MIKE21FM模型相結(jié)合的方法,并對不同滲透系數(shù)下防波堤透水性對溫排水過程的影響進(jìn)行了模擬,結(jié)果表明,防波堤的透水性對港池內(nèi)部水體溫升的影響顯著,并且隨著滲透系數(shù)的增大而增大,在進(jìn)行與溫排水相關(guān)的工程建設(shè)時,建筑物的透水性應(yīng)該引起高度重視。

    [1]王慶改, 戴文楠, 趙曉宏, 等. 基于Mike21FM的來賓電廠擴(kuò)建工程溫排水?dāng)?shù)值模擬研究[J]. 環(huán)境科學(xué)研究, 2009(3):332-336. WANG Q G, DAI W N, ZHAO X H, et al. Numerical simulation research was done on hot water discharged by Laibin power plant impacting the surrounding waters environment based on Mike21FM[J]. Research of Environment Sciences,2009(3):332-336.

    [2]丁躍平, 郭遠(yuǎn)明, 李鐵軍, 等. 三門灣春季溫排水增溫效應(yīng)數(shù)值模擬研究[J]. 海洋科學(xué), 2015(4):101-107. DING Y P, GUO Y M, LI T J, et al. The numerical simulation studies of the effect of the increasing temperature on warmly discharged water in Sanmen bay[J]. Marine Sciences, 2015(4):101-107.

    [3]趙懿珺, 王靖宇, 陳小莉. 鐵山灣內(nèi)電廠溫排水累積熱影響數(shù)值模擬[J]. 水利水電科技進(jìn)展, 2015(2):47-51. ZHAO Y J, WANG J Y, CHEN X L. Numerical simulation of accumulated effect of thermal discharge from power plants in Tieshan Bay[J]. Advances in Science and Technology of Water Resources,2015(2):47-51.

    [4]周巧菊. 大亞灣海域溫排水三維數(shù)值模擬[J]. 海洋湖沼通報, 2007(4):37-46. ZHOU Q J. A Three-dimensional Numercal Simulation of Thermal Discharge in Daya Bay[J].Transactions of Oceanology and Limnology,2007(4): 37-46.

    [5]陳雪峰,劉松.長興島某電廠溫排水?dāng)?shù)值模擬研究[J].水道港口,2016,37(3):269-274. CHEN X F,LIU S.Numerical simulation for the plant of Changxing Island in Dalian[J].Journal of Waterway and Harbor,2016,37(3):269-274.

    [6]Rodi W. Turbulence models and their applications in hydraulics[M].United States: IHAR publication,1984.

    [7]Iqbal M. An Introduction to solar Radiation [M].New York: Academic Press, 1983.

    [8]吳持恭.水力學(xué): 第四版[M].北京:高等教育出版社,2007.

    [9]王向峰.馬鑾海堤開口改造防滲技術(shù)研究[D] .成都:西南交通大學(xué), 2009.

    Numerical simulation study on the thermal water discharge of one power plant in Vietnam

    LIShao-wu1,HANFei1,CHENHan-bao2

    (1.StateKeyLaboratoryofHydraulicEngineeringSimulationandSafety,TianjinUniversity,Tianjin300072,China; 2.TianjinResearchInstituteforWaterTransportEngineering,NationalEngineeringLaboratoryforPortHydraulicConstructionTechnology,KeyLaboratoryofEngineeringSediment,MinistryofTransport,Tianjin300456,China)

    In order to understand the effect of the newly-built Phase I and III project on the Phase II project of one power plant of Vietnam, a hydrodynamic and transport model of MIKE 21 were adopted to simulate tidal flow field and water temperature field. The thermal water diffusion process in the synchronizing operating case of Phase I, II and III, for both the case of spring tide and neap tide, summer and winter, were simulated. The simulation results show that the rise of water temperature around the water intakes in Phase II project does not exceed 1℃ after the Phase I and II projects putting into operation, which meets the requirements of operation. A combination method of Dupuit formula and MIKE 21 model was proposed to simulate the permeable breakwater, based on which, the effect of the permeability of the rubble mound breakwater on thermal water diffusion was evaluated. The results show that the permeability of breakwater has a great effect on thermal water diffusion.

    MIKE 21; thermal discharge; permeable breakwater; Dupuit formula

    TM 621;O 242.1

    :A

    :1005-8443(2017)04-0361-07

    2016-09-30;

    :2016-10-27

    李紹武(1962-),男,山東萊州人,教授,主要從事海岸動力學(xué)及岸灘演變、海岸建筑物波浪數(shù)值模型研究。

    Biography: LI Shao-wu (1962-) , male, professor.

    猜你喜歡
    港池防波堤取水口
    水資源取水口數(shù)字化管理技術(shù)探討
    治淮(2022年3期)2023-01-03 11:05:22
    寧夏地表水一級取水口評價與調(diào)整
    中國水利(2022年1期)2023-01-03 09:16:29
    基于Boussinesq波浪模型的港池波浪數(shù)值模擬與泊穩(wěn)分析
    海洋通報(2021年3期)2021-08-14 02:20:50
    黃河中游干流取水口設(shè)計討論
    陜西水利(2021年10期)2021-01-18 08:32:47
    基于數(shù)值水槽的半封閉港池水體交換能力研究
    水道港口(2020年3期)2020-07-28 06:49:26
    寬肩臺防波堤穩(wěn)定性數(shù)值模擬方法研究
    海洋通報(2020年6期)2020-03-19 02:10:34
    關(guān)于浮式防波堤消能效果及透射系數(shù)的研究
    提升環(huán)抱式港池水交換能力研究進(jìn)展
    水道港口(2018年5期)2018-12-04 08:18:44
    深降幅水電站進(jìn)水口分層取水口設(shè)計
    頂升平臺在強(qiáng)涌浪海域深水防波堤地基處理中的應(yīng)用
    免费在线观看视频国产中文字幕亚洲| 欧美不卡视频在线免费观看 | 一本久久中文字幕| 国产精品久久久久久亚洲av鲁大| 国产熟女午夜一区二区三区| 99re在线观看精品视频| 丝袜美足系列| 制服诱惑二区| 丰满的人妻完整版| 又黄又爽又免费观看的视频| 久久久精品国产亚洲av高清涩受| 亚洲av成人不卡在线观看播放网| 成年女人毛片免费观看观看9| 女生性感内裤真人,穿戴方法视频| 午夜精品国产一区二区电影| 99国产精品99久久久久| 日韩高清综合在线| 成人三级黄色视频| 精品国产国语对白av| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲男人天堂网一区| 精品少妇一区二区三区视频日本电影| 精品国内亚洲2022精品成人| 亚洲精华国产精华精| 久久天堂一区二区三区四区| 精品福利观看| 一本久久中文字幕| 美女大奶头视频| 亚洲一区中文字幕在线| 欧美中文日本在线观看视频| 亚洲第一av免费看| 丁香六月欧美| 99在线视频只有这里精品首页| 亚洲精品一区av在线观看| 欧美日本视频| 一本久久中文字幕| 国产黄a三级三级三级人| 每晚都被弄得嗷嗷叫到高潮| 亚洲国产中文字幕在线视频| 好男人在线观看高清免费视频 | 欧美日本视频| 三级毛片av免费| 99国产精品一区二区蜜桃av| 神马国产精品三级电影在线观看 | 亚洲国产看品久久| 色尼玛亚洲综合影院| 狠狠狠狠99中文字幕| 久久狼人影院| 久热爱精品视频在线9| 99久久精品国产亚洲精品| 国产一卡二卡三卡精品| ponron亚洲| 国产欧美日韩综合在线一区二区| 欧美午夜高清在线| 人人澡人人妻人| 琪琪午夜伦伦电影理论片6080| 麻豆久久精品国产亚洲av| 国产av一区在线观看免费| 亚洲精品中文字幕一二三四区| 国产一区二区三区视频了| 国产高清videossex| 久久久久国内视频| 国产午夜精品久久久久久| 自线自在国产av| 九色国产91popny在线| 成年人黄色毛片网站| 精品久久久久久,| 亚洲欧美激情在线| 成人18禁高潮啪啪吃奶动态图| 日本 欧美在线| 亚洲激情在线av| 国产精品国产高清国产av| 色综合欧美亚洲国产小说| 欧美黑人精品巨大| 日韩欧美国产在线观看| 亚洲七黄色美女视频| www日本在线高清视频| 自拍欧美九色日韩亚洲蝌蚪91| 大型av网站在线播放| 久久久国产欧美日韩av| 人人妻,人人澡人人爽秒播| 久久精品国产99精品国产亚洲性色 | 给我免费播放毛片高清在线观看| av中文乱码字幕在线| 美女国产高潮福利片在线看| 国产一区二区三区视频了| 日本三级黄在线观看| 国产精品日韩av在线免费观看 | 日日干狠狠操夜夜爽| 老汉色av国产亚洲站长工具| 亚洲欧美精品综合一区二区三区| 亚洲在线自拍视频| 母亲3免费完整高清在线观看| 丰满的人妻完整版| 激情视频va一区二区三区| 国产精品野战在线观看| 久久精品国产亚洲av高清一级| 亚洲av日韩精品久久久久久密| 亚洲伊人色综图| 日日夜夜操网爽| 婷婷丁香在线五月| 男女之事视频高清在线观看| 无遮挡黄片免费观看| 欧美成人一区二区免费高清观看 | 侵犯人妻中文字幕一二三四区| 国产精品 欧美亚洲| 男人的好看免费观看在线视频 | 长腿黑丝高跟| 国产免费男女视频| 精品国产乱码久久久久久男人| 国产午夜福利久久久久久| 91在线观看av| 亚洲成国产人片在线观看| 搡老妇女老女人老熟妇| av中文乱码字幕在线| 精品高清国产在线一区| 精品无人区乱码1区二区| 嫩草影院精品99| 一进一出抽搐gif免费好疼| 搞女人的毛片| 黄片播放在线免费| 国产欧美日韩一区二区精品| 在线观看午夜福利视频| 成人亚洲精品一区在线观看| 欧美av亚洲av综合av国产av| 精品国产国语对白av| 法律面前人人平等表现在哪些方面| 日本一区二区免费在线视频| 国产精品久久久av美女十八| 中文亚洲av片在线观看爽| 一级a爱视频在线免费观看| 色婷婷久久久亚洲欧美| 淫妇啪啪啪对白视频| 波多野结衣一区麻豆| 18美女黄网站色大片免费观看| 美女大奶头视频| 999久久久精品免费观看国产| 精品国内亚洲2022精品成人| 日韩欧美在线二视频| 啪啪无遮挡十八禁网站| 90打野战视频偷拍视频| 每晚都被弄得嗷嗷叫到高潮| 欧美亚洲日本最大视频资源| 国产乱人伦免费视频| 久久精品亚洲熟妇少妇任你| 人人澡人人妻人| 正在播放国产对白刺激| 免费搜索国产男女视频| 成人免费观看视频高清| av视频免费观看在线观看| 国产野战对白在线观看| 91字幕亚洲| 色老头精品视频在线观看| 搡老熟女国产l中国老女人| 亚洲国产欧美网| 成人精品一区二区免费| 国内精品久久久久精免费| 最近最新中文字幕大全电影3 | 一区二区三区高清视频在线| 91精品国产国语对白视频| 精品第一国产精品| 老司机福利观看| 久久久精品国产亚洲av高清涩受| 亚洲精品av麻豆狂野| 欧美色视频一区免费| 亚洲精华国产精华精| 亚洲男人天堂网一区| 国产精品日韩av在线免费观看 | 成人三级做爰电影| 一进一出抽搐动态| 久久伊人香网站| 人人澡人人妻人| 久久久久久亚洲精品国产蜜桃av| 69精品国产乱码久久久| 欧美色视频一区免费| 国产色视频综合| 18美女黄网站色大片免费观看| 欧美日韩中文字幕国产精品一区二区三区 | 又黄又粗又硬又大视频| 亚洲欧美日韩高清在线视频| 国产aⅴ精品一区二区三区波| 最新在线观看一区二区三区| 一边摸一边抽搐一进一小说| 亚洲av日韩精品久久久久久密| 久久精品国产清高在天天线| 两个人免费观看高清视频| 午夜福利一区二区在线看| 此物有八面人人有两片| 久久精品91蜜桃| 亚洲午夜精品一区,二区,三区| 老熟妇乱子伦视频在线观看| 久久精品亚洲熟妇少妇任你| 中文字幕精品免费在线观看视频| 不卡av一区二区三区| 欧美精品啪啪一区二区三区| 两个人看的免费小视频| 国内精品久久久久精免费| 99国产精品99久久久久| 久久久久久久久中文| 国产精品野战在线观看| 中文字幕高清在线视频| 欧美最黄视频在线播放免费| 国产在线精品亚洲第一网站| 脱女人内裤的视频| 咕卡用的链子| 亚洲国产欧美网| 中文字幕高清在线视频| 国产1区2区3区精品| 免费在线观看日本一区| 淫妇啪啪啪对白视频| 日韩欧美国产在线观看| www.熟女人妻精品国产| 久久午夜综合久久蜜桃| 国产日韩一区二区三区精品不卡| 两个人看的免费小视频| 一级毛片女人18水好多| 可以在线观看毛片的网站| 好看av亚洲va欧美ⅴa在| 多毛熟女@视频| 后天国语完整版免费观看| 黄色视频,在线免费观看| 亚洲自拍偷在线| 女人被躁到高潮嗷嗷叫费观| 免费无遮挡裸体视频| 久久精品亚洲熟妇少妇任你| 亚洲精品国产色婷婷电影| 正在播放国产对白刺激| 欧美最黄视频在线播放免费| 一区二区三区国产精品乱码| 超碰成人久久| 老熟妇仑乱视频hdxx| 88av欧美| 免费看十八禁软件| 麻豆国产av国片精品| 亚洲免费av在线视频| 一级毛片精品| 激情在线观看视频在线高清| 欧美av亚洲av综合av国产av| 亚洲精品美女久久久久99蜜臀| 亚洲精品国产精品久久久不卡| 亚洲九九香蕉| 欧美在线一区亚洲| 母亲3免费完整高清在线观看| 一本久久中文字幕| 久久人妻熟女aⅴ| 最好的美女福利视频网| 亚洲美女黄片视频| 一区二区日韩欧美中文字幕| 国产蜜桃级精品一区二区三区| 老汉色av国产亚洲站长工具| 国产av精品麻豆| а√天堂www在线а√下载| a级毛片在线看网站| 欧美成人一区二区免费高清观看 | 国产三级在线视频| 亚洲自拍偷在线| 日韩欧美三级三区| 国产午夜精品久久久久久| 午夜福利免费观看在线| 中文字幕另类日韩欧美亚洲嫩草| 一个人免费在线观看的高清视频| 欧美日韩乱码在线| 老鸭窝网址在线观看| 别揉我奶头~嗯~啊~动态视频| 黄色片一级片一级黄色片| 大码成人一级视频| 黄片播放在线免费| 老司机深夜福利视频在线观看| 手机成人av网站| 精品一区二区三区视频在线观看免费| 伊人久久大香线蕉亚洲五| 精品一品国产午夜福利视频| 深夜精品福利| 丝袜美腿诱惑在线| 一区二区三区精品91| 免费在线观看黄色视频的| 国产单亲对白刺激| 精品久久久久久久久久免费视频| 两性午夜刺激爽爽歪歪视频在线观看 | 9191精品国产免费久久| 激情在线观看视频在线高清| 日本黄色视频三级网站网址| 欧美亚洲日本最大视频资源| 最好的美女福利视频网| 一a级毛片在线观看| 免费高清视频大片| 视频区欧美日本亚洲| 日韩高清综合在线| 国内精品久久久久久久电影| 婷婷六月久久综合丁香| 久久久国产精品麻豆| 青草久久国产| 精品欧美一区二区三区在线| 国产成人免费无遮挡视频| 九色亚洲精品在线播放| 999久久久精品免费观看国产| 男女做爰动态图高潮gif福利片 | 露出奶头的视频| 这个男人来自地球电影免费观看| 国产麻豆成人av免费视频| 国产精品 欧美亚洲| 中文字幕高清在线视频| 亚洲三区欧美一区| 国产真人三级小视频在线观看| 国产欧美日韩精品亚洲av| 免费人成视频x8x8入口观看| 免费在线观看亚洲国产| 亚洲欧洲精品一区二区精品久久久| 香蕉久久夜色| 一进一出抽搐gif免费好疼| 久久精品aⅴ一区二区三区四区| or卡值多少钱| 成年人黄色毛片网站| 免费久久久久久久精品成人欧美视频| 热re99久久国产66热| 97人妻天天添夜夜摸| 欧美日本亚洲视频在线播放| 精品乱码久久久久久99久播| 亚洲精品中文字幕一二三四区| 三级毛片av免费| 久久人妻av系列| 久久精品亚洲熟妇少妇任你| 男人舔女人的私密视频| 黄片小视频在线播放| 高清黄色对白视频在线免费看| 国产国语露脸激情在线看| 黑人欧美特级aaaaaa片| 麻豆成人av在线观看| 久久国产精品男人的天堂亚洲| 亚洲av熟女| 久久久久国产一级毛片高清牌| 丁香六月欧美| 精品国产亚洲在线| 在线观看一区二区三区| 国产精品 国内视频| 久久国产精品人妻蜜桃| 亚洲成人久久性| 在线免费观看的www视频| 如日韩欧美国产精品一区二区三区| svipshipincom国产片| 亚洲一区二区三区不卡视频| 成人国产一区最新在线观看| 国产免费av片在线观看野外av| 久久人妻福利社区极品人妻图片| 国产精华一区二区三区| 亚洲精品国产一区二区精华液| 国产免费男女视频| 亚洲成人精品中文字幕电影| av在线天堂中文字幕| 亚洲一卡2卡3卡4卡5卡精品中文| 免费观看精品视频网站| 精品一区二区三区视频在线观看免费| 亚洲成人国产一区在线观看| 久久久久国内视频| 亚洲精品国产色婷婷电影| 老熟妇乱子伦视频在线观看| 一级a爱片免费观看的视频| 国产av精品麻豆| 欧美乱妇无乱码| 我的亚洲天堂| 国产高清视频在线播放一区| 欧美日本中文国产一区发布| av福利片在线| 日日爽夜夜爽网站| 日韩精品青青久久久久久| 两人在一起打扑克的视频| 精品欧美一区二区三区在线| 无限看片的www在线观看| 最好的美女福利视频网| 亚洲人成77777在线视频| 精品一区二区三区四区五区乱码| 亚洲色图综合在线观看| 身体一侧抽搐| 十八禁网站免费在线| 亚洲五月婷婷丁香| 免费观看精品视频网站| 丝袜在线中文字幕| av欧美777| 国产av精品麻豆| 老司机午夜十八禁免费视频| 国产精品电影一区二区三区| 身体一侧抽搐| 亚洲精品美女久久av网站| 国产色视频综合| 我的亚洲天堂| 麻豆久久精品国产亚洲av| АⅤ资源中文在线天堂| 亚洲av第一区精品v没综合| 亚洲午夜精品一区,二区,三区| 国产又爽黄色视频| 日韩av在线大香蕉| 高清在线国产一区| 亚洲欧美激情在线| АⅤ资源中文在线天堂| 成年人黄色毛片网站| 午夜免费鲁丝| 黄色 视频免费看| 欧美日韩中文字幕国产精品一区二区三区 | 久久国产精品男人的天堂亚洲| 亚洲av第一区精品v没综合| 熟女少妇亚洲综合色aaa.| 777久久人妻少妇嫩草av网站| 国产蜜桃级精品一区二区三区| 国产成人影院久久av| 九色国产91popny在线| 国产精品影院久久| 欧美激情高清一区二区三区| www国产在线视频色| 老司机在亚洲福利影院| 久久天堂一区二区三区四区| 母亲3免费完整高清在线观看| 亚洲精品国产一区二区精华液| 亚洲专区国产一区二区| 他把我摸到了高潮在线观看| 日韩大码丰满熟妇| 可以在线观看毛片的网站| 午夜视频精品福利| 久久久久久国产a免费观看| 久久久精品国产亚洲av高清涩受| 男男h啪啪无遮挡| 免费少妇av软件| 一级,二级,三级黄色视频| 91大片在线观看| 一进一出抽搐gif免费好疼| 极品教师在线免费播放| 99久久99久久久精品蜜桃| 女人精品久久久久毛片| 精品卡一卡二卡四卡免费| 村上凉子中文字幕在线| or卡值多少钱| 国产97色在线日韩免费| 亚洲最大成人中文| 每晚都被弄得嗷嗷叫到高潮| 欧美老熟妇乱子伦牲交| 这个男人来自地球电影免费观看| 国产精品一区二区在线不卡| 美国免费a级毛片| 国产高清视频在线播放一区| 丰满的人妻完整版| 国产亚洲欧美在线一区二区| 国产精品亚洲av一区麻豆| 激情视频va一区二区三区| 天天添夜夜摸| 欧美一级a爱片免费观看看 | 国产精品一区二区三区四区久久 | 男人舔女人下体高潮全视频| 国产高清有码在线观看视频 | 国产精品永久免费网站| 精品福利观看| 国产成年人精品一区二区| 国产国语露脸激情在线看| 女人爽到高潮嗷嗷叫在线视频| 久久国产精品男人的天堂亚洲| 18禁裸乳无遮挡免费网站照片 | 免费观看人在逋| 日本三级黄在线观看| 老司机靠b影院| 免费观看人在逋| 成人18禁在线播放| 久久久久久免费高清国产稀缺| 精品久久久久久久毛片微露脸| 91麻豆av在线| 精品不卡国产一区二区三区| 国产欧美日韩精品亚洲av| 无人区码免费观看不卡| 色婷婷久久久亚洲欧美| 99精品欧美一区二区三区四区| 精品卡一卡二卡四卡免费| 18禁观看日本| 亚洲av电影在线进入| 99久久综合精品五月天人人| 777久久人妻少妇嫩草av网站| 美女高潮喷水抽搐中文字幕| av片东京热男人的天堂| 亚洲午夜理论影院| 巨乳人妻的诱惑在线观看| 国产色视频综合| 香蕉久久夜色| 宅男免费午夜| 欧美色欧美亚洲另类二区 | 十八禁网站免费在线| 99久久精品国产亚洲精品| 亚洲午夜精品一区,二区,三区| 丰满的人妻完整版| 热99re8久久精品国产| 又黄又爽又免费观看的视频| 黄片小视频在线播放| 非洲黑人性xxxx精品又粗又长| 搞女人的毛片| 深夜精品福利| 久久久国产成人精品二区| 亚洲三区欧美一区| 日韩欧美免费精品| 91老司机精品| 日韩 欧美 亚洲 中文字幕| 男人操女人黄网站| 超碰成人久久| 免费在线观看完整版高清| 在线免费观看的www视频| 黄色片一级片一级黄色片| 又黄又爽又免费观看的视频| 久久精品影院6| 色av中文字幕| 色哟哟哟哟哟哟| 黄频高清免费视频| 91大片在线观看| 最好的美女福利视频网| 日韩高清综合在线| 久久青草综合色| 欧美在线黄色| 两性夫妻黄色片| 成人国产综合亚洲| 美女高潮到喷水免费观看| 国产激情久久老熟女| 国产高清视频在线播放一区| 老鸭窝网址在线观看| 亚洲av成人一区二区三| 日韩欧美一区二区三区在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美激情久久久久久爽电影 | 少妇的丰满在线观看| 在线观看66精品国产| 亚洲无线在线观看| 两个人看的免费小视频| 亚洲精品一区av在线观看| 男人操女人黄网站| 欧美日韩瑟瑟在线播放| 精品久久久久久久毛片微露脸| 天堂影院成人在线观看| 97碰自拍视频| 日韩高清综合在线| 国产伦一二天堂av在线观看| 国产精品一区二区三区四区久久 | 国产精品,欧美在线| 亚洲一区高清亚洲精品| 精品国产亚洲在线| 两个人视频免费观看高清| 精品国产国语对白av| 亚洲熟妇中文字幕五十中出| 中出人妻视频一区二区| 国产精品综合久久久久久久免费 | 人人妻人人澡人人看| 在线观看免费视频日本深夜| 成在线人永久免费视频| 亚洲欧美日韩高清在线视频| 1024香蕉在线观看| 美女扒开内裤让男人捅视频| 亚洲中文字幕一区二区三区有码在线看 | 国产精品免费视频内射| 国产av精品麻豆| 一级毛片女人18水好多| 少妇被粗大的猛进出69影院| 女人被狂操c到高潮| 亚洲一区二区三区色噜噜| 一级片免费观看大全| 欧美乱码精品一区二区三区| 88av欧美| 国产精品一区二区在线不卡| 啦啦啦免费观看视频1| 97超级碰碰碰精品色视频在线观看| 日本免费一区二区三区高清不卡 | 欧美日韩精品网址| 在线天堂中文资源库| 极品人妻少妇av视频| 亚洲成人久久性| 久久精品91无色码中文字幕| www日本在线高清视频| 嫩草影院精品99| 窝窝影院91人妻| 日本精品一区二区三区蜜桃| 欧洲精品卡2卡3卡4卡5卡区| 久久久国产欧美日韩av| 免费在线观看亚洲国产| 在线观看www视频免费| 久久精品成人免费网站| 亚洲一区中文字幕在线| 高清黄色对白视频在线免费看| 国产高清视频在线播放一区| 高清黄色对白视频在线免费看| 香蕉久久夜色| 亚洲狠狠婷婷综合久久图片| 9191精品国产免费久久| 久久久久久大精品| 黄色毛片三级朝国网站| 91精品三级在线观看| 最新美女视频免费是黄的| 久久亚洲精品不卡| 色综合欧美亚洲国产小说| 亚洲国产精品sss在线观看| 久久这里只有精品19| www.自偷自拍.com| 女警被强在线播放| 天天躁夜夜躁狠狠躁躁| 99精品久久久久人妻精品| 色哟哟哟哟哟哟| 亚洲成国产人片在线观看| 激情在线观看视频在线高清| 久久人妻熟女aⅴ| 免费在线观看视频国产中文字幕亚洲| 久久久久久大精品| 99久久99久久久精品蜜桃| 亚洲精品美女久久av网站| 看免费av毛片| 免费高清视频大片| 亚洲色图综合在线观看| 99在线人妻在线中文字幕| 亚洲精品av麻豆狂野| 久久香蕉国产精品| 亚洲在线自拍视频| 午夜福利,免费看| 国产在线精品亚洲第一网站| 美女高潮喷水抽搐中文字幕| xxx96com| АⅤ资源中文在线天堂| 国内精品久久久久精免费|