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

    長江宜昌段溢油事故影響數學模型分析研究

    2018-04-03 07:41:52辛小康陳和春王繼保李光浩
    關鍵詞:江段豐水期溢油

    胡 旭,辛小康,尹 煒,陳和春,王繼保,李光浩

    (1.三峽大學 水利與環(huán)境學院,湖北 宜昌 443002;2.長江水資源保護科學研究所,武漢 430051)

    長江是我國主要的河流運輸通道,客貨運輸量大。交通運輸部發(fā)布《2016年交通運輸行業(yè)發(fā)展統(tǒng)計公報》顯示[1],內河運輸完成貨物量達35.72億噸,占水路總運輸貨物量的55.97%;內河港口貨物吞吐量完成47.76億噸,增長3.1%。隨著國家依托黃金水道推動長江經濟帶發(fā)展戰(zhàn)略的實施,內陸自貿實驗區(qū)相繼開放,貨運量在今后仍將繼續(xù)增加[2]。在同條件及預警技術管理水平不變前提下,這無疑增大了溢油事故風險發(fā)生的概率。眾所周知,溢油事故產生的石油類污染物會給整個河流水體造成巨大水環(huán)境污染,嚴重破壞水環(huán)境平衡,危害水生生物。相關研究表明,中華鱘在含油濃度為10 mg/L時,死亡時間僅為2天[3]。

    為了更好地應對內河溢油事故風險,有必要對內河溢油事故進行模擬研究。目前,工程界和學術界對溢油事故的模擬研究多采用數值模擬方法[4-5]。現有研究對象多數為入??诩昂Q蟮貐^(qū)。對于環(huán)境差異較大、水文條件特殊、河道情況復雜多變的內河,溢油事故研究較少[6-9]。本文采用MIKE21軟件,建立平面二維溢油模型,模擬葛洲壩至下游枝江江段豐水期和枯水期突發(fā)溢油事故,并分析不同溢油量條件下油膜的擴展、輸移范圍及其對下游造成的影響。

    1 模型計算原理

    1.1 水動力模塊

    MIKE21是由丹麥水力研究院研究開發(fā)的應用于近海海域、海灣、港區(qū)及河道等水環(huán)境模擬的系統(tǒng)。水動力模塊(HD)是其核心,可以模擬由于各種力的作用而產生的水位及水流變化。其控制方程是基于不可壓縮流體和Reynolds值均布的Navier-Stokes方程,并遵從Boussinesq假設和靜水壓力假設。二維非恒定淺水方程組如下:

    式中,t為時間,η為水面高程,h為總水深,g為重力加速度,ρ為水的密度,ρ0為(淡)水的參考密度,f=2Ωsin?為科氏力系數(Ω為旋轉角速率,?為地理緯度),Pa為大氣壓強,Sij為輻射應力張量,S和(us,vs)分別為點源的排放量和速度,uˉ和vˉ分別為沿x,y方向上的流速在深度上的平均值;(τsx,τsy)和(τbx,τby)為水面風應力張量和河床床面應力張量,Tij為側向應力。

    該模塊利用隱式交替方向ADI(Alternating di?rection implicit)技術對上述連續(xù)方程和動量方程進行離散,所得矩陣用追趕法求解,各項微分和主要系數均采用中心差分格式,Taylor級數展開的截斷誤差可達二至三階精度。

    1.2 溢油模塊

    MIKE21的Oil Spill模塊(OS模塊)是在HD模塊基礎之上建立的溢油模型,用于模擬溢油擴散運動規(guī)律和歸宿。溢油模塊采用“油粒子”對溢油量進行概化,采用拉格朗日方法描述油粒子的輸移軌跡,包括擴展、漂移、擴散等過程[9],其輸移方程見公式(4)-(7)。同時,在輸移過程中油粒子也會發(fā)生如蒸發(fā)、乳化、溶解等風化過程,該軟件通過計算油粒子質量損失來體現。

    (1)擴展:油膜擴展運動采用修正的Fay重力-粘力公式計算。公式如下:

    式中:Ao為油膜面積;t為時間;Kα為擴散系數;Vo為油膜體積;Ro為油膜半徑,h0為初始油膜厚度,取10 cm。

    (2)漂移:影響油粒子漂移速度的主要因素是水流和風作用力。漂移速度為:

    式中:UW為水面上10 m處的風速,Us為表面流速,cw為漂移系數,一般取值0.02~0.03。

    (3)擴散:由于單個粒子不能被分成幾片,因此擴散的過程被解釋為在隨機方向上的運動。對于二維的情況,可以將隨機走動的距離形式表示為一個時間步長α方向上的擴散距離。計算公式如下:

    式中:Sα為在α方向上的一個時間步長內可能擴散走動的距離,Dα為α方向上的擴散系數,R為-1到1的隨機數。Δtp為計算時間步長。

    2 研究區(qū)域概況

    本文模擬范圍為長江葛洲壩至下游枝江江段。上游起始斷面為葛洲壩大壩,下游終止斷面為枝江昌門溪,范圍全長80 km。根據《全國重要江河湖泊水功能區(qū)劃(2011-2030年)》,該江段包含長江宜昌飲用水源、工業(yè)用水區(qū)和長江宜昌中華鱘保護區(qū)兩個水功能區(qū)。根據水功能區(qū)管理要求,長江宜昌飲用水源、工業(yè)用水區(qū)水質目標為Ⅱ類,水質目標要求較高。模擬范圍如圖1。

    圖1 模擬范圍內各水功能區(qū)、監(jiān)測斷面、概化溢油事故點相對位置示意圖Fig.1 Relative position of each water function area,monitoring section,and generalized point of oil spill accident

    3 模型參數設定及計算方案

    3.1 網格劃分

    本文采用2015年實測水下地形數據。葛洲壩至枝江區(qū)間的長江干流,河道形狀相對規(guī)則,變化梯度不大。為了保證更高計算精度和收斂性,模型網格劃分時采用四邊形網格,網格尺寸為200 m×40 m;支流清江入匯口段,地形復雜,區(qū)域較小,采用三角形網格劃分,三角形網格控制面積小于1 000 m2。計算區(qū)域共剖分計算網格單元16 093個,并根據該實測水下地形進行插值,如圖2和圖3所示。

    圖2 葛洲壩-枝江區(qū)間水質數學模型計算網格Fig.2 Computing grid for water quality mathematical model between Gezhouba and Zhijiang section

    圖3 葛洲壩-枝江區(qū)間計算網格地形插值Fig.3 Terrain interpolation for computing grid between Gezhouba and Zhijiang section

    3.2 水動力模塊邊界條件及參數設定

    葛洲壩是三峽的反調節(jié)水庫,為日調節(jié)水庫,設計水文條件可參照三峽出庫流量。三峽水庫2003年才開始初期蓄水,2010年才開始175 m試驗性蓄水,水文情勢尚不穩(wěn)定。于是本文便采用三峽水庫近10年最枯、最豐月的月均流量作為設計水文條件。本模型枯水期上游邊界設計流量為三峽2月份平均出庫流量5 280 m3/s,參照枝江水位站水位—流量關系曲線,得出模型下游邊界設計水位為33.42 m。豐水期上游邊界設計流量為三峽8月份平均出庫流量21 357 m3/s,對應下游邊界設計水位為36.15 m。參照多年平均水溫,枯水期設定水溫10℃,汛期設定水溫20℃。

    設定濕水深度為0.100 m,淹沒深度為0.050 m,干水深度為0.005 m,根據辛小康[10]等對長江宜昌江段水動力模型的率定結果,曼寧系數為0.031,渦粘性系數的Smagorinsky系數為0.28,擴散系數D設定為1 m2/s。

    3.3 溢油模塊參數設定

    本文首先計算各個油粒子的位置變化和組成變化,然后統(tǒng)計各個網格上的油粒子個數和組分含量,模擬出油膜的濃度時空分布和組分變化。MIKE21中溢油模型將單個油粒子質量定義為揮發(fā)性輕組分的質量、不揮發(fā)性重組分的質量、蠟狀物的質量以及瀝青的質量。本文只考慮輕組分和重組分,兩者比例為3∶2,即1個1 000 kg油粒子,包含輕組分600 kg,重組分400 kg。由于石油類污染物易發(fā)生蒸發(fā)、溶解、乳化、沉淀、生物降解、光氧化等特點,本文參考姜衛(wèi)星[11]2007年對狹長型感潮河流溢油模型的研究,采用縱向擴散系數為0.25 m2/s、橫向為0.10 m2/s。風速和風向,可根據氣象資料獲得。

    枯水期葛洲壩下游江段石油類背景濃度采用宜昌黃陵廟斷面2010-2014年2月份實測平均值0.03 mg/L。豐水期采用該斷面2010-2014年8月份實測平均值0.01 mg/L。

    3.4 計算方案

    本次模擬設計溢油點位置為該區(qū)域航運交通要塞三江引航道和大江引航道的匯合點,坐標經緯度為(E111°16′4.8″,N30°42′10.8″),垂向水深為零。參考我國《船舶油污染事故等級》(JT/458-2001)[12],考慮最不利情形,即油船造成的水域油污染事故。根據油船大事故和重大事故的劃分依據,擬取設計總溢油量分別為5 t和10 t,并分別模擬枯水期和豐水期兩種不同水文條件下的溢油事故情形,共4種模擬工況,提取每種工況1 h,3 h,10 h后的溢油模擬計算結果進行對比分析。

    4 計算結果與分析

    4.1 枯水期

    枯水期瞬時溢油為10 t時,事故發(fā)生后1 h,3 h和10 h以后油膜位置和面積如圖4所示??梢钥闯觯斠缬?0 t后1 h,油膜尚未完全展開,油膜厚度≥0.3μm的范圍為600 m×300 m;油膜厚度大部分為100μm~0.1 mm的范圍為400 m×200 m。溢油后3 h,油膜厚度≥0.3μm的范圍為3600 m×300 m;油膜厚度為100μm~0.1 mm的范圍為3 200 m×280 m。溢油10 h以后,油類物質經過擴散,污染帶范圍進一步增加,但由于乳化、蒸發(fā)、風化、光解等綜合作用,油膜厚度有所減小,油膜厚度≥0.3μm的范圍可達到14 000 m×400 m;油膜厚度為100μm~0.1 mm的范圍為600 m×200 m。

    圖4 枯水期溢油10 t后1 h,3 h,10 h油膜厚度分布圖Fig.4 Oil film thickness distribution at 1 h,3 h and 10 h after 10 t spillage during low water period

    枯水期瞬時溢油5 t后1 h,3 h和10 h油膜演進情況見圖5??梢钥闯觯斠缬秃? h,油膜厚度≥0.3μm的范圍為400 m×240 m;油膜厚度大部分為100μm~0.1 mm,其范圍為300 m×200 m。溢油后3 h,油膜厚度≥0.3μm的范圍為2 200 m×300 m;油膜厚度為100μm~0.1 mm的范圍為2 000 m×260m。溢油10 h以后,油類物質經過擴散,并有少量附著在河岸,污染帶范圍進一步增加,油膜厚度有所減小,油膜厚度≥0.3μm的范圍可達到10 000 m×400 m,而油膜厚度在100μm~0.1 mm幾乎消失。

    圖5 枯水期溢油5 t后1 h,3 h,10 h油膜厚度分布圖Fig.5 Oil film thickness distribution at 1 h,3 h and 10 h after 5 t spillage during low water period

    4.2 豐水期

    豐水期瞬時溢油10 t時,事故發(fā)生后1 h,3 h和10 h以后油膜位置和面積如圖6所示??梢钥闯?,當溢油10 t后1 h,油膜厚度為100μm~1 mm,范圍為800 m×100 m。溢油后3 h,油膜面積達到4 400 m×200 m,油膜厚度超過100μm。溢油10 h以后,油類物質經過擴散,污染帶范圍進一步增加,油膜范圍達到10 000 m×400 m,同樣由于乳化、蒸發(fā)、風化、光解等綜合作用,油膜厚度減小,基本小于100μm。

    圖6 豐水期溢油10 t后1 h,3 h,10 h油膜厚度分布圖Fig.6 Oil film thickness distribution at 1 h,3 h and 10 h after 10 t spillage during high water period

    豐水期瞬時溢油5 t時,事故發(fā)生后1 h,3 h和10 h以后油膜位置和面積如圖7所示??梢钥闯?,當溢油5 t后1 h,油膜尚未完全展開,油膜厚度為100μm~1 mm,范圍800 m×100 m。溢油后3 h,油膜面積達2 600 m×200 m,油膜厚度超過100μm。溢油10 h以后,油類物質經過擴散,飄移,污染帶范圍進一步增加,油膜范圍達到8 000 m×300 m,而油膜厚度減小,基本小于100μm。

    圖7 豐水期溢油5 t后1 h,3 h,10 h油膜厚度分布圖Fig.7 Oil film thickness distribution at1 h,3 h and 10 h after 5 t spillage during high water period

    4.3 結果分析

    根據模擬結果可知,對于所有工況,漏油5 t以上均會產生明顯的油膜厚度≥100μm的區(qū)域,且持續(xù)時間都超過3 h。溢油1 h后油膜處于或剛完成溢油擴散的初始階段,油膜展開范圍不大,厚度均大于100μm。溢油3 h后,油膜范圍繼續(xù)擴大,油膜影響范圍均大于2 000 m×260 m,最大可達3 600 m×300 m。溢油10 h以后,油類物質經過擴散,污染帶范圍進一步增加,無論是枯水期還是豐水期,兩種溢油量油膜影響范圍均大于8 000 m×300 m,最大可達10 000 m×400 m,但由于物理、化學等綜合作用,油膜厚度減小較明顯,普遍小于100μm。在同等水力條件下,溢油量越多油膜相對更厚,但油膜擴散范圍的差異較小,這可能與河道較狹窄有關。

    對比枯水期和豐水期可知,枯水期由于水量偏小,流速小,油污范圍較小,持續(xù)的時間較長,而豐水期水量大,流速大,油膜掃過的范圍相對更大更遠,油膜厚度減小速率更大,持續(xù)的時間也相對較短??菟谟湍哌^的范圍可達到猇亭一帶,豐水期可達到枝城一帶。無論是枯水期還是豐水期,當瞬時溢油量為5 t及以上時,均會對中華鱘自然保護區(qū)產生明顯影響。

    5 結論

    我國內河流域環(huán)境普遍復雜多變,與海洋有較大差別,對于內河溢油事故模擬,需要考慮的因素眾多。本文根據實測地形資料借用MIKE21軟件對宜昌江段進行溢油事故模擬預測,考慮了不同水文條件,不同溢油量的情形,而模擬結果也表明這對油膜厚度、飄移范圍、持續(xù)時間的差別有較為顯著的影響,這對內河溢油應急對策具有指導意義,也為內河溢油事故對水環(huán)境影響評估提供了很好的參考價值。

    [1]中華人民共和國交通運輸部.交通運輸行業(yè)發(fā)展統(tǒng)計公報[R].北京:綜合規(guī)劃司,2016.

    [2]袁群.淺談長江航運油污染現狀、原因及其對策[J].水運管理,2004(5):22-23.

    [3]倪朝輝,翟良安.石油對魚類等水生生物的毒性[J].淡水漁業(yè),1997,27(6):38-40.

    [4]陳家興,杜娟,付金鋒.基于GIS的三峽庫區(qū)溢油模型研究[J].中國水運,2010,10(9):96-98.

    [5]高龍驤.內河船舶溢油擴散的建模與仿真[D].上海:復旦大學,2013.

    [6]趙琰鑫,王永桂,張萬順,等.河道溢油污染事故二維數值模型研究[J].人民長江,2012,43(15):81-84.

    [7]GARCIA M R,TOVAR F H.Computer modeling of oil spill trajectories with a high accuracy method[J].Spill Science&Technology Bulletin,1999,5(5/6):323-330.

    [8]巫麗俊.感潮江段事故溢油二維數值模擬[D].南京:河海大學,2006.

    [9]張帆,黃立文,鄧建,等.重慶主城區(qū)江段溢油模型及數值試驗研究[J].武漢理工大學學報(交通科學與工程版),2011,35(1):87-90.

    [10]辛小康,葉閩,尹煒.長江宜昌江段水污染事故的水庫調度措施研究[J].水電能源科學,2011,29(06):46-48,95.

    [11]姜衛(wèi)星.黃浦江溢油事故的數值模擬研究[D].上海:同濟大學,2007.

    [12]中華人民共和國交通部.船舶油污染事故等級:JT/458-2001[S].北京:中國標準出版社,2001.

    猜你喜歡
    江段豐水期溢油
    豐水期巢湖水體與沉積物重金屬污染及分配特征研究
    牛欄江流域浮游動物群落結構特征及驅動因子
    膠東沿海某地區(qū)生活飲用水總α、總β放射性水平調查分析
    松花江干流哈爾濱江段封凍
    中國水利(2020年23期)2020-12-24 06:24:40
    近岸溢油漂移擴散預測方法研究——以膠州灣溢油事件為例
    海洋通報(2020年2期)2020-09-04 09:22:48
    基于GF-1衛(wèi)星的海上溢油定量監(jiān)測——以青島溢油事故為例
    海洋通報(2020年2期)2020-09-04 09:22:46
    黑龍江遜克江段跑冰排
    科學導報(2020年27期)2020-05-13 14:23:43
    小城大愛——監(jiān)利江段“東方之星”號游輪翻沉事件愛心幫扶紀實
    學習月刊(2015年23期)2015-07-09 05:41:56
    滇池流域入湖河流水文形貌特征對豐水期大型底棲動物群落結構的影響
    對白茆沙水域溢油事故后修復治理的思考
    中國水利(2015年4期)2015-02-28 15:12:23
    中文字幕人妻熟女乱码| 午夜av观看不卡| 老鸭窝网址在线观看| 中文字幕人妻熟女乱码| 在线观看三级黄色| 日韩成人av中文字幕在线观看| 精品少妇久久久久久888优播| 亚洲欧美清纯卡通| 最近2019中文字幕mv第一页| 伊人久久国产一区二区| 久久人人97超碰香蕉20202| 一区二区三区激情视频| 香蕉精品网在线| 十八禁高潮呻吟视频| 精品视频人人做人人爽| 欧美日本中文国产一区发布| 成人毛片60女人毛片免费| 欧美人与性动交α欧美精品济南到 | 大片免费播放器 马上看| 卡戴珊不雅视频在线播放| 亚洲精品日本国产第一区| 精品亚洲成国产av| 中文字幕人妻丝袜制服| 日日撸夜夜添| 久久国内精品自在自线图片| 国产xxxxx性猛交| 欧美激情高清一区二区三区 | 91精品国产国语对白视频| 免费看av在线观看网站| 在线精品无人区一区二区三| 国产极品天堂在线| 人人妻人人添人人爽欧美一区卜| 亚洲国产成人一精品久久久| 国产免费现黄频在线看| 中国三级夫妇交换| 夜夜骑夜夜射夜夜干| 亚洲国产精品国产精品| 丝袜美腿诱惑在线| 欧美精品国产亚洲| 久久久久国产精品人妻一区二区| 国产精品香港三级国产av潘金莲 | videosex国产| 美女国产视频在线观看| 超碰成人久久| 久久精品国产亚洲av天美| 性高湖久久久久久久久免费观看| 免费在线观看黄色视频的| 精品少妇久久久久久888优播| 国产精品一二三区在线看| 18禁国产床啪视频网站| 欧美日韩精品成人综合77777| 日韩成人av中文字幕在线观看| 久久久精品区二区三区| 国产一区二区激情短视频 | 国产野战对白在线观看| 最近最新中文字幕大全免费视频 | 少妇被粗大猛烈的视频| 国产成人a∨麻豆精品| 在线亚洲精品国产二区图片欧美| 精品国产一区二区三区久久久樱花| 少妇人妻 视频| 91精品国产国语对白视频| 男人舔女人的私密视频| 日韩熟女老妇一区二区性免费视频| 观看av在线不卡| 9色porny在线观看| 激情视频va一区二区三区| 亚洲男人天堂网一区| 两性夫妻黄色片| 中国国产av一级| 久久人人爽人人片av| 亚洲欧美成人精品一区二区| 欧美日韩一级在线毛片| 日日摸夜夜添夜夜爱| 欧美成人午夜免费资源| 国产福利在线免费观看视频| 黑人欧美特级aaaaaa片| 999精品在线视频| 成年人免费黄色播放视频| 国产一区二区三区av在线| 亚洲精品国产色婷婷电影| 九色亚洲精品在线播放| 免费大片黄手机在线观看| 99香蕉大伊视频| 国精品久久久久久国模美| 999久久久国产精品视频| 99国产精品免费福利视频| 国产 一区精品| 校园人妻丝袜中文字幕| 精品亚洲成a人片在线观看| 亚洲欧美精品自产自拍| 亚洲欧洲精品一区二区精品久久久 | 黄片播放在线免费| 亚洲精品国产av蜜桃| 成人免费观看视频高清| 亚洲国产欧美在线一区| 麻豆精品久久久久久蜜桃| 1024香蕉在线观看| 哪个播放器可以免费观看大片| 免费久久久久久久精品成人欧美视频| 久久人人爽av亚洲精品天堂| 欧美 亚洲 国产 日韩一| 水蜜桃什么品种好| 精品亚洲成a人片在线观看| 热99久久久久精品小说推荐| 蜜桃在线观看..| 最新的欧美精品一区二区| 国产日韩一区二区三区精品不卡| 亚洲成人一二三区av| 亚洲精品美女久久久久99蜜臀 | 麻豆av在线久日| 亚洲国产精品一区三区| 午夜精品国产一区二区电影| 永久网站在线| 丰满迷人的少妇在线观看| 亚洲精品中文字幕在线视频| 成年人免费黄色播放视频| 久久久久精品久久久久真实原创| 一本久久精品| 尾随美女入室| 久久久久久伊人网av| 超碰成人久久| 精品少妇黑人巨大在线播放| 性少妇av在线| 少妇 在线观看| 亚洲第一青青草原| 最黄视频免费看| 精品国产国语对白av| 在线观看一区二区三区激情| 亚洲精品aⅴ在线观看| √禁漫天堂资源中文www| 一级a爱视频在线免费观看| 黄片小视频在线播放| 婷婷色综合大香蕉| 欧美国产精品va在线观看不卡| 午夜福利乱码中文字幕| 午夜av观看不卡| 精品卡一卡二卡四卡免费| 久久久国产欧美日韩av| 成人手机av| 黄色毛片三级朝国网站| 日韩视频在线欧美| 日本午夜av视频| 香蕉精品网在线| 久久女婷五月综合色啪小说| 黑人猛操日本美女一级片| 蜜桃国产av成人99| 欧美xxⅹ黑人| 天美传媒精品一区二区| 国产免费又黄又爽又色| 26uuu在线亚洲综合色| 午夜91福利影院| 麻豆精品久久久久久蜜桃| 亚洲国产精品999| 2022亚洲国产成人精品| 国产亚洲欧美精品永久| 国产精品免费视频内射| 亚洲精品美女久久久久99蜜臀 | 最黄视频免费看| 国语对白做爰xxxⅹ性视频网站| 男女无遮挡免费网站观看| 91国产中文字幕| 男男h啪啪无遮挡| 国产精品久久久av美女十八| 男人操女人黄网站| 久久久久精品性色| 亚洲经典国产精华液单| 成人影院久久| 天天躁狠狠躁夜夜躁狠狠躁| 18在线观看网站| 亚洲,一卡二卡三卡| 一级爰片在线观看| 婷婷色麻豆天堂久久| 91在线精品国自产拍蜜月| 老女人水多毛片| 久久久亚洲精品成人影院| 久久久久久久精品精品| 日韩一区二区视频免费看| 国产精品无大码| 9191精品国产免费久久| 亚洲国产精品999| 国产成人91sexporn| 丝袜脚勾引网站| 亚洲在久久综合| 日日撸夜夜添| 女人被躁到高潮嗷嗷叫费观| 欧美精品人与动牲交sv欧美| 高清不卡的av网站| 亚洲精品一二三| 国产一区二区在线观看av| 成年人午夜在线观看视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久国内精品自在自线图片| 亚洲久久久国产精品| 成人国产av品久久久| 99久久人妻综合| 国产精品 国内视频| 精品少妇久久久久久888优播| 国产精品麻豆人妻色哟哟久久| 只有这里有精品99| 2021少妇久久久久久久久久久| 蜜桃在线观看..| 观看av在线不卡| av电影中文网址| 久久婷婷青草| av福利片在线| 国产精品一二三区在线看| 叶爱在线成人免费视频播放| 黄色怎么调成土黄色| 少妇精品久久久久久久| 欧美bdsm另类| 亚洲五月色婷婷综合| 久久精品人人爽人人爽视色| 国产综合精华液| 男女午夜视频在线观看| 超色免费av| 丰满少妇做爰视频| 不卡视频在线观看欧美| 精品一区在线观看国产| 亚洲国产毛片av蜜桃av| 国产麻豆69| 国产精品99久久99久久久不卡 | 日韩不卡一区二区三区视频在线| 97精品久久久久久久久久精品| 久久精品熟女亚洲av麻豆精品| 美女大奶头黄色视频| 我的亚洲天堂| 乱人伦中国视频| 成人毛片60女人毛片免费| 啦啦啦在线免费观看视频4| 精品一区二区三卡| 精品久久久久久电影网| 亚洲精品中文字幕在线视频| 亚洲一区二区三区欧美精品| 老司机亚洲免费影院| 亚洲久久久国产精品| 爱豆传媒免费全集在线观看| 亚洲精品中文字幕在线视频| 国产成人精品久久二区二区91 | 中文字幕精品免费在线观看视频| 国产免费现黄频在线看| 91在线精品国自产拍蜜月| 日韩电影二区| 中文字幕色久视频| 女人高潮潮喷娇喘18禁视频| 99香蕉大伊视频| 狂野欧美激情性bbbbbb| 免费女性裸体啪啪无遮挡网站| 亚洲精品久久午夜乱码| 亚洲视频免费观看视频| 国产97色在线日韩免费| 天堂俺去俺来也www色官网| 欧美日韩一区二区视频在线观看视频在线| 亚洲成人手机| av视频免费观看在线观看| 午夜日本视频在线| 国产高清国产精品国产三级| 麻豆乱淫一区二区| 亚洲,欧美精品.| 两性夫妻黄色片| 嫩草影院入口| av在线老鸭窝| 99香蕉大伊视频| 人妻少妇偷人精品九色| 免费观看av网站的网址| 午夜老司机福利剧场| 捣出白浆h1v1| 成人漫画全彩无遮挡| 亚洲第一av免费看| 亚洲综合色惰| 亚洲av在线观看美女高潮| 亚洲人成电影观看| 成年人午夜在线观看视频| 久久国产精品大桥未久av| 国产成人av激情在线播放| 亚洲三级黄色毛片| 一本色道久久久久久精品综合| 国产欧美日韩一区二区三区在线| 如日韩欧美国产精品一区二区三区| 免费在线观看视频国产中文字幕亚洲 | 国产一区二区在线观看av| 超碰97精品在线观看| 纵有疾风起免费观看全集完整版| 夜夜骑夜夜射夜夜干| 日韩熟女老妇一区二区性免费视频| 丰满迷人的少妇在线观看| 精品视频人人做人人爽| 啦啦啦在线观看免费高清www| 成年av动漫网址| 黄片无遮挡物在线观看| 99热全是精品| freevideosex欧美| 亚洲精品第二区| 免费女性裸体啪啪无遮挡网站| 亚洲五月色婷婷综合| 在线看a的网站| 久久精品国产综合久久久| 久久久久精品人妻al黑| 久久人人爽av亚洲精品天堂| 久久av网站| 熟女少妇亚洲综合色aaa.| 午夜福利视频在线观看免费| 成人国产av品久久久| 80岁老熟妇乱子伦牲交| 久久久国产一区二区| 国产野战对白在线观看| 十分钟在线观看高清视频www| 久久久久久久国产电影| 婷婷色综合大香蕉| 中文字幕亚洲精品专区| 欧美人与性动交α欧美精品济南到 | 两个人免费观看高清视频| 亚洲av福利一区| 亚洲伊人色综图| 国产 一区精品| 亚洲三区欧美一区| 亚洲精品久久久久久婷婷小说| 国产精品嫩草影院av在线观看| 人体艺术视频欧美日本| 久久精品国产综合久久久| 美女高潮到喷水免费观看| 综合色丁香网| 老汉色av国产亚洲站长工具| 日韩成人av中文字幕在线观看| 国产精品不卡视频一区二区| 成人18禁高潮啪啪吃奶动态图| 一区二区日韩欧美中文字幕| 午夜久久久在线观看| av线在线观看网站| 亚洲美女视频黄频| 高清在线视频一区二区三区| 少妇人妻久久综合中文| 国产乱人偷精品视频| 精品国产乱码久久久久久小说| 亚洲欧美日韩另类电影网站| 欧美精品人与动牲交sv欧美| 国产精品免费视频内射| 国产黄色免费在线视频| 日韩 亚洲 欧美在线| 狠狠婷婷综合久久久久久88av| 午夜久久久在线观看| 老女人水多毛片| 视频区图区小说| 精品99又大又爽又粗少妇毛片| 国产一区二区三区av在线| 建设人人有责人人尽责人人享有的| 母亲3免费完整高清在线观看 | 国产白丝娇喘喷水9色精品| 综合色丁香网| 欧美激情极品国产一区二区三区| 国产成人av激情在线播放| 欧美变态另类bdsm刘玥| av天堂久久9| 考比视频在线观看| 在线观看国产h片| 少妇的逼水好多| √禁漫天堂资源中文www| 国产精品香港三级国产av潘金莲 | 桃花免费在线播放| 国产精品久久久av美女十八| 日韩av免费高清视频| 久久久a久久爽久久v久久| 亚洲伊人久久精品综合| 久久久久久久大尺度免费视频| 婷婷色综合大香蕉| av线在线观看网站| 各种免费的搞黄视频| 91成人精品电影| 少妇人妻久久综合中文| 国产精品麻豆人妻色哟哟久久| 少妇的丰满在线观看| 免费黄频网站在线观看国产| 中文精品一卡2卡3卡4更新| 国产精品麻豆人妻色哟哟久久| 中文字幕色久视频| 久久这里有精品视频免费| 国产精品 欧美亚洲| 考比视频在线观看| 激情五月婷婷亚洲| 亚洲成av片中文字幕在线观看 | 伊人久久国产一区二区| 两性夫妻黄色片| 成年人午夜在线观看视频| 边亲边吃奶的免费视频| 精品久久久精品久久久| 黄色配什么色好看| 一区二区三区四区激情视频| 老司机影院毛片| 美国免费a级毛片| 狠狠婷婷综合久久久久久88av| 超碰成人久久| 一级,二级,三级黄色视频| 中文乱码字字幕精品一区二区三区| 成人毛片60女人毛片免费| 中文乱码字字幕精品一区二区三区| 免费高清在线观看日韩| 久久午夜福利片| 大陆偷拍与自拍| 99久久精品国产国产毛片| 黄频高清免费视频| 黄色视频在线播放观看不卡| 国产黄频视频在线观看| 香蕉精品网在线| 人人妻人人澡人人看| 丝瓜视频免费看黄片| 色婷婷av一区二区三区视频| 少妇精品久久久久久久| 夫妻午夜视频| 深夜精品福利| 啦啦啦视频在线资源免费观看| 在线天堂中文资源库| 日韩av免费高清视频| 另类亚洲欧美激情| 国产亚洲最大av| 另类亚洲欧美激情| 我的亚洲天堂| 国产高清不卡午夜福利| 国产精品熟女久久久久浪| 老鸭窝网址在线观看| 一级a爱视频在线免费观看| 国产精品偷伦视频观看了| 亚洲成人一二三区av| 69精品国产乱码久久久| 久久久久久久国产电影| 少妇的丰满在线观看| 秋霞伦理黄片| 成人影院久久| 美女高潮到喷水免费观看| 99九九在线精品视频| 狠狠婷婷综合久久久久久88av| av女优亚洲男人天堂| 2022亚洲国产成人精品| 纯流量卡能插随身wifi吗| 国产一区二区在线观看av| 一级爰片在线观看| av天堂久久9| 国产一区二区三区综合在线观看| 咕卡用的链子| 成人毛片a级毛片在线播放| 一级a爱视频在线免费观看| 两性夫妻黄色片| 男人舔女人的私密视频| 精品午夜福利在线看| 国产精品国产av在线观看| 黑人欧美特级aaaaaa片| 欧美最新免费一区二区三区| 中文乱码字字幕精品一区二区三区| 欧美 亚洲 国产 日韩一| 国产女主播在线喷水免费视频网站| 最近的中文字幕免费完整| 国产精品一国产av| 亚洲精品一区蜜桃| www.精华液| a级片在线免费高清观看视频| 黄色怎么调成土黄色| 欧美最新免费一区二区三区| 下体分泌物呈黄色| 亚洲熟女精品中文字幕| 国产乱人偷精品视频| 欧美日韩亚洲国产一区二区在线观看 | 天堂8中文在线网| 精品酒店卫生间| 高清欧美精品videossex| 男女啪啪激烈高潮av片| 美女视频免费永久观看网站| 日韩制服骚丝袜av| 日韩一区二区视频免费看| 欧美日韩视频高清一区二区三区二| 色94色欧美一区二区| 色婷婷久久久亚洲欧美| 两性夫妻黄色片| 亚洲av在线观看美女高潮| 十八禁高潮呻吟视频| 啦啦啦在线观看免费高清www| 18禁国产床啪视频网站| 精品卡一卡二卡四卡免费| 永久免费av网站大全| 热re99久久精品国产66热6| 日韩精品有码人妻一区| 久久久a久久爽久久v久久| 有码 亚洲区| 人妻一区二区av| 精品国产国语对白av| 欧美精品人与动牲交sv欧美| 女人被躁到高潮嗷嗷叫费观| 大码成人一级视频| 欧美老熟妇乱子伦牲交| 9191精品国产免费久久| 亚洲国产精品一区二区三区在线| 国产精品国产三级国产专区5o| av福利片在线| 国产日韩欧美在线精品| 国产成人av激情在线播放| 国产av精品麻豆| 久久精品夜色国产| av在线观看视频网站免费| 精品99又大又爽又粗少妇毛片| xxx大片免费视频| 人人妻人人爽人人添夜夜欢视频| 国语对白做爰xxxⅹ性视频网站| 男女边摸边吃奶| 飞空精品影院首页| 欧美精品国产亚洲| 久久久久精品人妻al黑| 免费日韩欧美在线观看| 国产 一区精品| 少妇精品久久久久久久| 亚洲情色 制服丝袜| 欧美最新免费一区二区三区| 亚洲四区av| 免费人妻精品一区二区三区视频| 成年av动漫网址| 一级片'在线观看视频| 国产成人一区二区在线| 黄色毛片三级朝国网站| 在线看a的网站| 婷婷成人精品国产| 日韩中文字幕视频在线看片| 亚洲国产精品999| 人妻人人澡人人爽人人| 久久精品亚洲av国产电影网| 亚洲少妇的诱惑av| 精品一区在线观看国产| 亚洲人成77777在线视频| 一二三四在线观看免费中文在| 99香蕉大伊视频| 精品国产一区二区三区四区第35| 搡老乐熟女国产| 亚洲中文av在线| 下体分泌物呈黄色| 日韩伦理黄色片| 国产在线一区二区三区精| 一级片'在线观看视频| 赤兔流量卡办理| 久久国产亚洲av麻豆专区| 国产福利在线免费观看视频| 亚洲国产成人一精品久久久| 国产精品免费大片| 国产白丝娇喘喷水9色精品| 国产精品99久久99久久久不卡 | 毛片一级片免费看久久久久| 男女边摸边吃奶| 亚洲av在线观看美女高潮| 亚洲欧洲精品一区二区精品久久久 | 久久毛片免费看一区二区三区| 久久国产亚洲av麻豆专区| 久久久久久久久久久久大奶| 国产深夜福利视频在线观看| 国产高清国产精品国产三级| 国产有黄有色有爽视频| 免费看av在线观看网站| 亚洲综合色网址| 亚洲精品美女久久av网站| 曰老女人黄片| 国产精品秋霞免费鲁丝片| 久热这里只有精品99| 一二三四中文在线观看免费高清| 亚洲av欧美aⅴ国产| 亚洲精品国产av蜜桃| 国产精品 欧美亚洲| 国产一区二区三区综合在线观看| a 毛片基地| 日韩制服骚丝袜av| 亚洲,欧美精品.| 亚洲av免费高清在线观看| 男人操女人黄网站| 秋霞伦理黄片| 欧美日韩国产mv在线观看视频| 丝袜人妻中文字幕| 亚洲国产毛片av蜜桃av| 国精品久久久久久国模美| 天天躁夜夜躁狠狠久久av| 高清在线视频一区二区三区| 天堂8中文在线网| 国产女主播在线喷水免费视频网站| 久久久久久人妻| 丰满乱子伦码专区| 黑丝袜美女国产一区| 免费看av在线观看网站| 熟女少妇亚洲综合色aaa.| 90打野战视频偷拍视频| 一边摸一边做爽爽视频免费| 国产97色在线日韩免费| 日韩,欧美,国产一区二区三区| 曰老女人黄片| 极品人妻少妇av视频| 欧美人与善性xxx| 五月天丁香电影| 精品国产一区二区久久| 免费在线观看完整版高清| av不卡在线播放| 免费播放大片免费观看视频在线观看| 1024香蕉在线观看| www日本在线高清视频| 亚洲美女视频黄频| 人妻一区二区av| 亚洲精品国产av成人精品| 最黄视频免费看| kizo精华| 一本久久精品| 久久久亚洲精品成人影院| 亚洲成人一二三区av| av片东京热男人的天堂| 国产色婷婷99| 久久99热这里只频精品6学生| 交换朋友夫妻互换小说| av在线app专区| 亚洲 欧美一区二区三区| 亚洲精品日韩在线中文字幕| 久久久a久久爽久久v久久| 日日摸夜夜添夜夜爱| 亚洲精品,欧美精品| 成年女人在线观看亚洲视频| 一级毛片我不卡| 哪个播放器可以免费观看大片|