亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于MIKE21的大型湖庫型飲用水水源保護(hù)區(qū)劃分研究*
        ——以花亭湖水庫為例

        2021-02-01 02:20:08言文杰王曉輝石先陽
        環(huán)境污染與防治 2021年1期
        關(guān)鍵詞:水質(zhì)模型

        言文杰 王曉輝 石先陽#

        (1.安徽大學(xué)資源與環(huán)境工程學(xué)院,安徽 合肥 230601;2.安徽省環(huán)境科學(xué)研究院,安徽 合肥 230061)

        飲用水水源地保護(hù)是關(guān)系國計(jì)民生的重大問題,飲用水水源地環(huán)境管理受到世界各國的廣泛關(guān)注[1-6]。由于人類活動(dòng)加劇和不合理開發(fā)等因素,我國約有70%以上的水源地已遭受不同程度的污染[7]。湖庫型飲用水水源地在城鎮(zhèn)飲水安全工程中占有重要地位[8],花亭湖作為典型的大型湖庫型飲用水水源地,是太湖縣近58萬人民的主要供水水源,其水質(zhì)保護(hù)受到高度重視,而飲用水水源保護(hù)區(qū)劃分在水源地水質(zhì)保護(hù)中起到重要的基礎(chǔ)性作用。根據(jù)新修訂的《飲用水水源保護(hù)區(qū)劃分技術(shù)規(guī)范》(HJ 338—2018),大江大河和湖庫水源地皆鼓勵(lì)采用數(shù)值模型法進(jìn)行保護(hù)區(qū)劃分。目前,諸多學(xué)者基于一維、二維水環(huán)境數(shù)值模擬模型在飲用水水源保護(hù)區(qū)劃分方面做了大量的實(shí)踐與研究工作,如:李云禎等[9]使用《飲用水水源保護(hù)區(qū)劃分技術(shù)規(guī)范》(HJ/T 338—2007)中提出的排放口在兩個(gè)邊界中間的連續(xù)點(diǎn)源穩(wěn)態(tài)排放二維水質(zhì)模型,根據(jù)模擬點(diǎn)源排出的污染物濃度衰減到《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB 3838—2002)Ⅱ類、Ⅲ類水質(zhì)的距離劃分了飲用水水源一、二級保護(hù)區(qū);焦軍麗等[10]利用MIKE模型系列中的一維水質(zhì)模型MIKE11對萊壩飲用水水源地月波—涼姜溝斷面進(jìn)行了水質(zhì)沿程變化規(guī)律的模擬,并結(jié)合類比經(jīng)驗(yàn)法確定了水源保護(hù)區(qū)的長度范圍;周訓(xùn)華等[11]將一維點(diǎn)源稀釋自凈模型與完全混合式模型相結(jié)合,發(fā)現(xiàn)了赤田水庫上游點(diǎn)源排污對水庫的影響較小,只限于局部區(qū)域,并據(jù)此制定了保護(hù)區(qū)的劃分方案;汪洋等[12]基于GMS軟件構(gòu)建了三維非穩(wěn)定流地下水?dāng)?shù)值模型,對鄒城市某巖溶地下水源保護(hù)區(qū)進(jìn)行了劃分。但是,過往數(shù)值模型法劃分水源保護(hù)區(qū)研究多圍繞河流及地下水型水源地,且以一維水動(dòng)力水質(zhì)模型為主,基于二維水動(dòng)力-水質(zhì)耦合模型的湖庫型水源保護(hù)區(qū)劃分研究仍鮮有應(yīng)用。

        本研究選取安慶市太湖縣花亭湖水庫作為研究對象,建立二維湖庫水動(dòng)力-水質(zhì)耦合模型MIKE21,設(shè)置枯水期、平水期、豐水期3種不同情況,選取水庫代表性水質(zhì)指標(biāo)TN、TP、DO進(jìn)行模型的率定與驗(yàn)證,從而確定對水庫水質(zhì)造成較大波動(dòng)的時(shí)期(豐水期),根據(jù)該時(shí)期內(nèi)水質(zhì)指標(biāo)濃度分布以及水源地水質(zhì)保護(hù)的實(shí)際需要,結(jié)合數(shù)值模擬結(jié)果與類比經(jīng)驗(yàn)法進(jìn)行飲用水水源保護(hù)區(qū)劃分,使劃分結(jié)果具有科學(xué)性與合理性,以期為湖庫型水源地水域保護(hù)區(qū)劃分提供科學(xué)依據(jù)和借鑒案例。

        1 區(qū)域概況

        花亭湖水庫位于安徽省安慶市太湖縣境內(nèi),地處大別山南麓、長江北岸,介于北緯30°9′~30°46′、東經(jīng)115°45′~116°30′?;ㄍず饔蛩募痉置?,具有北亞熱帶季風(fēng)氣候特征。冬季受西北內(nèi)陸氣流控制,氣溫較低;夏季受東南海洋氣流影響,炎熱潮濕;年平均氣溫16.4 ℃,年平均日照1 938 h,無霜期249 d,年平均降雨量1 368.4 mm,常年主導(dǎo)風(fēng)向?yàn)槲鞅憋L(fēng)與東南風(fēng),年平均風(fēng)速為3.2 m/s。

        花亭湖水庫距太湖縣城4 km,是太湖縣重要的地表水資源和城鎮(zhèn)水源地?;ㄍず畮熳鳛榘矐c市唯一的大型水庫,流域面積1 870 km2,多年平均入庫流量46.3 m3/s,多年平均來水量14.6億m3,總庫容23.66億m3,水庫正常蓄水位88.00 m,設(shè)計(jì)洪水位95.21 m?;ㄍず畮熳鳛橥掏滦秃?,其來流主要有寺前河、羅溪河、樹堰河、青石河與潘家河,常年平均入庫流量分別為:12.0、10.0、8.0、7.6、6.0 m3/s,花亭湖水庫通過大壩向下級河流泄水。水庫取水口位于大壩上游480 m南岸,主要供給太湖縣一水廠和二水廠,取水規(guī)模分別為1.5萬、4.5萬m3/d。

        2 模型搭建與應(yīng)用

        MIKE模型系列軟件是由丹麥水力學(xué)研究所開發(fā),在國內(nèi)許多河流中已經(jīng)廣泛應(yīng)用[13-16]。MIKE21主要用于模擬湖庫、河口及海灣水流、波浪、泥沙及水環(huán)境等。該模型主要包括水動(dòng)力、水質(zhì)、波浪與泥沙等模塊,可進(jìn)行復(fù)雜二維平面水流計(jì)算、洪水淹沒計(jì)算、水質(zhì)模擬預(yù)報(bào)和水環(huán)境治理規(guī)劃等。MIKE21中地形網(wǎng)格具有自動(dòng)剖分功能,既可采用不規(guī)則三角網(wǎng),也可采用規(guī)則多邊形進(jìn)行網(wǎng)格劃分。

        2.1 模型原理

        模型遵循N-S方程,并服從Boussinesq假設(shè),即流體低速流動(dòng)中,密度的變化僅考慮溫度的影響,不考慮由壓強(qiáng)所引起的影響,同時(shí)服從靜水壓假設(shè)。

        2.2 地形與網(wǎng)格

        以花亭湖水庫為研究對象,選取潘家河與青石河交匯處斷面(記為潘家河/情石河)至花亭湖大壩斷面為模擬范圍,水域東西長13.2 km,南北寬10.2 km。采用三角網(wǎng)格對模擬區(qū)域進(jìn)行網(wǎng)格劃分,對于岸線邊界網(wǎng)格適當(dāng)加密以提高模擬精度,網(wǎng)格共計(jì)8 381個(gè)。

        研究區(qū)原始高程數(shù)據(jù)來源于衛(wèi)星,分辨率為50 m×50 m。本研究在高程質(zhì)量檢查無誤基礎(chǔ)上,對花亭湖水庫模擬范圍水下地形采用加權(quán)反離法插值到每個(gè)網(wǎng)格節(jié)點(diǎn),得到花亭湖水庫飲用水水源地的數(shù)字高程模型。

        2.3 邊界條件

        本研究選取上游入湖支流來水為入流邊界條件,水庫下泄流量為出流條件,其余岸線為固壁邊界條件。模型水文與氣象邊界數(shù)據(jù)來源于安慶市水利局網(wǎng)站2018年逐日監(jiān)測數(shù)據(jù)(因降雨集中于1—10月,因此只關(guān)注該時(shí)段數(shù)據(jù))。圖1反映了花亭湖水庫降雨量與水位日變化。

        蒸發(fā)數(shù)據(jù)來源于中國氣象數(shù)據(jù)網(wǎng)月尺度蒸發(fā)量,結(jié)果見表1?;ㄍず畮?018年典型枯水期(1月)蒸發(fā)量20.5 mm,平水期(4月)蒸發(fā)量50.8 mm,豐水期(6月)蒸發(fā)量102.1 mm。

        本文所構(gòu)建的水質(zhì)模型選用浮游植物碳(PC)、浮游植物氮(PN)、浮游植物磷(PP)、葉綠素(Chl-a)、浮游動(dòng)物碳(ZC)、風(fēng)化碳(DC)、風(fēng)化氮(DN)、風(fēng)化磷(DP)、無機(jī)氮(IN)、無機(jī)磷(IP)、溶解氧(DO)作為狀態(tài)變量,TN、TP作為衍生變量。其中,Chl-a、DO、IP取值結(jié)合地方水功能區(qū)劃目標(biāo)與不同水期來水監(jiān)測數(shù)據(jù)確定,剩余的狀態(tài)變量則通過MIKE21推薦的轉(zhuǎn)換方法或者推薦的常數(shù)值來確定?;ㄍず饔蚋髦Я鱽硭克|(zhì)見表2,除流量外,表中其他取不同水期平均值。

        圖1 花亭湖水庫降雨量與水位日變化Fig.1 Diurnal process of rainfall and water level in Huating Lake Reservoir

        表1 2018年花亭湖流域月蒸發(fā)量

        表2 花亭湖流域來水水量水質(zhì)

        2.4 模型驗(yàn)證

        2.4.1 水動(dòng)力驗(yàn)證

        結(jié)合2018年1、4、6月水動(dòng)力過程模擬結(jié)果與現(xiàn)場監(jiān)測數(shù)據(jù)進(jìn)行豐水期、平水期、枯水期驗(yàn)證,主要率定和驗(yàn)證的指標(biāo)有渦黏系數(shù)、曼寧系數(shù)、風(fēng)應(yīng)力拖曳系數(shù)。

        模型水位監(jiān)測點(diǎn)位布置見圖2??紤]到與湖泊水位變化趨勢相近,選擇庫區(qū)監(jiān)測點(diǎn)位L1(花亭湖大壩上游500 m)實(shí)測的水位變化過程進(jìn)行模型的參數(shù)率定、驗(yàn)證工作。率定結(jié)果顯示,花亭湖水庫渦黏系數(shù)0.27 m2/s、曼寧系數(shù)32 m1/3/s、風(fēng)應(yīng)力拖曳系數(shù)0.001 255。

        圖2 花亭湖水庫監(jiān)測點(diǎn)位分布Fig.2 Monitoring point sketch distribution of Huating Lake Reservoir

        圖3至圖5分別為1、4、6月花亭湖水庫水位實(shí)測值與模擬值對比結(jié)果。驗(yàn)證結(jié)果顯示,1、4、6月模擬值與實(shí)測值相對誤差均在10%之內(nèi)。對比曲線趨勢和指標(biāo)誤差統(tǒng)計(jì)可以看出,模型模擬精度較好,能真實(shí)反映花亭湖水庫水動(dòng)力變化過程。水動(dòng)力模擬結(jié)果具備與水質(zhì)模塊進(jìn)行耦合的條件。

        圖3 1月水位模擬值與實(shí)測值Fig.3 Simulated and measured values of water level in January

        2.4.2 水質(zhì)驗(yàn)證

        參數(shù)率定是水質(zhì)模擬過程中的重要一環(huán)。在前人對花亭湖水庫進(jìn)行水質(zhì)模型研究的基礎(chǔ)上,結(jié)合2018年1、4、6月監(jiān)測點(diǎn)位L2(花亭湖大壩上游2 500 m)TN、TP以及DO水質(zhì)監(jiān)測結(jié)果,選取12個(gè)水質(zhì)敏感性指標(biāo),以默認(rèn)推薦指標(biāo)值為基值,通過分析花亭湖水庫水源地相關(guān)資料,變動(dòng)1個(gè)指標(biāo)而其他指標(biāo)均保持不變,指標(biāo)取值上下浮動(dòng)50%進(jìn)行模擬計(jì)算對比。

        圖4 4月水位模擬值與實(shí)測值Fig.4 Simulated and measured values of water level in April

        圖5 6月水位模擬與實(shí)測結(jié)果Fig.5 Simulated and measured values of water level in June

        經(jīng)過不斷率定與調(diào)試,各個(gè)水期水質(zhì)的模擬值與實(shí)測值對比見表3。整體來看,各水質(zhì)指標(biāo)的相對誤差均在20%以內(nèi),說明各個(gè)水期模型對污染物濃度的模擬效果較好,符合應(yīng)用要求。

        花亭湖水庫相關(guān)參數(shù)最終率定結(jié)果見表4。

        2.5 模擬分析

        結(jié)合2018年1、4、6月監(jiān)測與模擬結(jié)果可知,受上游山區(qū)來水和非點(diǎn)源的影響,花亭湖水庫汛期水質(zhì)指標(biāo)波動(dòng)較大,污染風(fēng)險(xiǎn)高于非汛期。因此,選擇2018年6月模擬結(jié)果開展污染物污染空間分布與保護(hù)區(qū)劃分工作。圖6至圖8分別為2018年6月庫區(qū)內(nèi)TN、TP與DO空間分布。

        3 保護(hù)區(qū)劃分

        根據(jù)HJ 338—2018要求,飲用水水源一級保護(hù)區(qū)水質(zhì)不得低于GB 3838—2002Ⅱ類水質(zhì)要求,二級保護(hù)區(qū)水質(zhì)不得低于Ⅲ類水質(zhì)要求。

        由MIKE21模擬的花亭湖水庫6月TN、TP濃度分布可知,TN、TP形成的污染物濃度帶分別在距離花亭湖大壩850、720 m的范圍內(nèi)達(dá)到了GB 3838—2002Ⅱ類水質(zhì)要求(TN≤0.5 mg/L,TP≤0.025 mg/L),而剩余水庫水域基本上均可以滿足Ⅲ類水質(zhì)要求(TN≤1.0 mg/L,TP≤0.05 mg/L),DO則在全部區(qū)域內(nèi)幾乎均能滿足Ⅱ類水質(zhì)要求(DO≥6 mg/L)。

        3.1 一級保護(hù)區(qū)劃分

        根據(jù)MIKE21模擬結(jié)果可知,若按照類比經(jīng)驗(yàn)法進(jìn)行劃分,取水口上游240~500 m的區(qū)域水質(zhì)將無法全部滿足GB 3838—2002Ⅱ類水質(zhì)要求。

        因此水域一級保護(hù)區(qū)宜根據(jù)MIKE21模型模擬的結(jié)果進(jìn)行劃分,范圍為花亭湖水庫大壩至上游720 m,寬度為庫區(qū)正常水位淹沒區(qū)域。一級保護(hù)區(qū)陸域根據(jù)HJ 338—2018類比經(jīng)驗(yàn)法要求,將水域邊界向沿岸陸域縱深延伸200 m,但不超過流域分水嶺和花亭湖大壩的陸地區(qū)域劃分為陸域一級保護(hù)區(qū)。

        3.2 二級保護(hù)區(qū)劃分

        根據(jù)MIKE21模擬結(jié)果,水域一級保護(hù)區(qū)以外全部水域的TN、TP、DO濃度基本上均能滿足GB 3838—2002Ⅲ類水質(zhì)要求,如果將一級保護(hù)區(qū)水域外的全部水體劃分為水域二級保護(hù)區(qū),則庫區(qū)內(nèi)劃分的水域面積將達(dá)到52.35 km2,遠(yuǎn)超出了飲用水水源地保護(hù)實(shí)際需求,將對山區(qū)居民出行(現(xiàn)有航道)和庫區(qū)生態(tài)旅游發(fā)展產(chǎn)生較大影響。因此,二級保護(hù)區(qū)宜選用類比經(jīng)驗(yàn)法進(jìn)行劃分,將一級保護(hù)區(qū)外徑向距離2 000 m,寬度為庫區(qū)正常水位淹沒的區(qū)域設(shè)為水域二級保護(hù)區(qū);將陸域一級保護(hù)區(qū)外徑向3 000 m的區(qū)域劃分為陸域二級保護(hù)區(qū),但不超過相應(yīng)的流域分水嶺。

        表3 監(jiān)測點(diǎn)位L2水質(zhì)模擬值與實(shí)測值對比

        表4 參數(shù)率定結(jié)果

        圖6 TN空間分布Fig.6 Spatial distribution of TN

        3.3 準(zhǔn)保護(hù)區(qū)劃分

        HJ 338—2018規(guī)定,按照湖庫流域范圍、污染源分布及對飲用水水源水質(zhì)的影響程度,二級保護(hù)區(qū)以外的匯水區(qū)域可以設(shè)定為準(zhǔn)保護(hù)區(qū)。本研究結(jié)合花亭湖水庫飲用水水源地實(shí)際情況以及保護(hù)需要,將一、二級保護(hù)區(qū)外的匯水區(qū)域劃定為準(zhǔn)保護(hù)區(qū)。

        3.4 保護(hù)區(qū)定界

        圖7 TP空間分布Fig.7 Spatial distribution of TP

        圖8 DO空間分布Fig.8 Spatial distribution of DO

        依據(jù)飲用水源保護(hù)區(qū)初步劃分范圍,再根據(jù)現(xiàn)場調(diào)研和資料搜集,并結(jié)合花亭湖水庫飲用水水源地的地形、地表、地物的特點(diǎn),最終確定各級保護(hù)區(qū)的界限及面積(見表5)。劃分一級保護(hù)區(qū)0.76 km2,二級保護(hù)區(qū)19.96 km2,準(zhǔn)保護(hù)區(qū)1 858.25 km2。

        使用ArcGIS軟件制作花亭湖水庫水源一、二級保護(hù)區(qū)及準(zhǔn)保護(hù)區(qū)劃分圖,結(jié)果見圖9至圖10。

        4 劃分合理性討論

        一級保護(hù)區(qū)的劃分過程中利用數(shù)值模型準(zhǔn)確識(shí)別了水質(zhì)指標(biāo)濃度分布,確保了一級保護(hù)區(qū)水域滿足GB 3838—2002 Ⅱ類水質(zhì)要求。該劃分方法劃分的范圍較類比經(jīng)驗(yàn)法劃分范圍有所縮小,但考慮到劃分區(qū)域內(nèi)無點(diǎn)源排放口,周邊僅有極少的原住居民區(qū)且無景區(qū),面源污染對該區(qū)域水質(zhì)的影響不大。因此,使用數(shù)值模型法劃分的一級保護(hù)區(qū)水域可以達(dá)到防范水質(zhì)污染、保障飲水安全的目的。同時(shí),將距離取水口較近的旅游大道劃入陸域一級保護(hù)區(qū)范圍并嚴(yán)格管理,可預(yù)防交通運(yùn)輸移動(dòng)源的環(huán)境風(fēng)險(xiǎn)。

        表5 花亭湖水庫水源保護(hù)區(qū)范圍與面積

        圖9 花亭湖水庫飲用水水源一、二級保護(hù)區(qū)劃分Fig.9 Delineation of Huating Lake Reservoir drinking water source’s key and secondary protection area

        圖10 花亭湖水庫飲用水水源準(zhǔn)保護(hù)區(qū)劃分Fig.10 Delineation of Huating Lake Reservoir drinking water source’s quasi protection area

        二級保護(hù)區(qū)陸域范圍包含了約12.36 km2的林地、茶園、經(jīng)濟(jì)果樹林、水田等易造成山區(qū)徑流面源污染的區(qū)域,可有效監(jiān)管、防護(hù)雨季由于雨水徑流挾帶腐殖質(zhì)等營養(yǎng)物質(zhì)進(jìn)入水庫而造成的庫區(qū)水質(zhì)污染。同時(shí),將污水統(tǒng)一收集、處理工作還不是很完善的禪心島、情人島、休閑山莊、西風(fēng)禪寺景區(qū)以及較為分散的原住居民區(qū)納入二級保護(hù)區(qū),按照二級保護(hù)區(qū)要求對其進(jìn)行相應(yīng)管控,達(dá)到了水源有效保護(hù)的目的。二級保護(hù)區(qū)的劃分也充分考慮到了社會(huì)經(jīng)濟(jì)因素,劃分的陸域二級保護(hù)區(qū)中僅有14.21 km2屬于花亭湖風(fēng)景區(qū)范圍(占風(fēng)景區(qū)總面積的4%),同時(shí)將風(fēng)景區(qū)中的核心景區(qū)二祖禪堂景區(qū)、佛圖寺景區(qū)、大壩景區(qū)、龍山景區(qū)排除在二級保護(hù)區(qū)外,以劃入準(zhǔn)保護(hù)區(qū)的形式進(jìn)行保護(hù),為當(dāng)?shù)厣鷳B(tài)旅游發(fā)展預(yù)留了足夠的空間,同時(shí)保障了當(dāng)?shù)鼐用竦氖杖?。二級保護(hù)區(qū)內(nèi)有約2.93 km的航線,大大方便了山區(qū)居民水路出行。

        可見,結(jié)合數(shù)值模型法與類比經(jīng)驗(yàn)法劃分水源保護(hù)區(qū),可兼顧技術(shù)與經(jīng)濟(jì)可行性操作原則,具有一定的科學(xué)性和合理性?;ㄍず畮祜嬘盟幢Wo(hù)區(qū)的劃分兼顧了流域水質(zhì)污染風(fēng)險(xiǎn)控制與景區(qū)生態(tài)旅游發(fā)展遠(yuǎn)期規(guī)劃需求,達(dá)到了水源保護(hù)與經(jīng)濟(jì)發(fā)展的雙贏目的。

        5 結(jié) 論

        (1) 花亭湖水庫二維水動(dòng)力-水質(zhì)耦合模型MIKE21在飲用水水源保護(hù)區(qū)劃分工作中具有較強(qiáng)的適用性,各項(xiàng)水質(zhì)指標(biāo)的模擬精度滿足模型應(yīng)用要求。

        (2) 結(jié)合數(shù)值模型法與類比經(jīng)驗(yàn)法,劃分花亭湖水庫飲用水水源一級保護(hù)區(qū)面積0.76 km2,二級保護(hù)區(qū)面積19.96 km2,準(zhǔn)保護(hù)區(qū)面積1 858.25 km2。本次劃分兼顧了流域水質(zhì)污染風(fēng)險(xiǎn)控制與景區(qū)生態(tài)旅游發(fā)展遠(yuǎn)期規(guī)劃需求,達(dá)到了水源保護(hù)與經(jīng)濟(jì)發(fā)展的雙贏目的。

        猜你喜歡
        水質(zhì)模型
        一半模型
        水質(zhì)抽檢豈容造假
        環(huán)境(2023年5期)2023-06-30 01:20:01
        重要模型『一線三等角』
        重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
        一月冬棚養(yǎng)蝦常見水質(zhì)渾濁,要如何解決?這9大原因及處理方法你要知曉
        這條魚供不應(yīng)求!蝦蟹養(yǎng)殖戶、垂釣者的最愛,不用投喂,還能凈化水質(zhì)
        圖像識(shí)別在水質(zhì)檢測中的應(yīng)用
        電子制作(2018年14期)2018-08-21 01:38:16
        3D打印中的模型分割與打包
        濟(jì)下水庫徑流水質(zhì)和垂向水質(zhì)分析及評價(jià)
        FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
        中文字幕亚洲中文第一| 新婚人妻不戴套国产精品| 国产一区二区三区四色av| 中文字幕网伦射乱中文| 亚洲欧美日韩综合久久| 日本高清一区二区不卡视频| 熟女白浆精品一区二区| 国产精品亚洲在钱视频| 亚洲天堂精品成人影院| 午夜不卡无码中文字幕影院| 亚洲欧美一区二区三区 | 人妻丰满av无码中文字幕| 96精品免费视频大全| 国产成人av三级三级三级在线 | 久久国产精品一区av瑜伽| 中文字幕亚洲乱码成熟女1区| 久久精品成人无码观看不卡| 欧美孕妇xxxx做受欧美88| 亚洲精品中文字幕不卡在线| 免费看黄片视频在线观看| 少妇高潮av久久久久久| 精品无码国产自产野外拍在线| 久久与欧美视频| 国产精品亚洲一二三区| 国产精品国产三级第一集| 久久天天躁狠狠躁夜夜不卡| 婷婷四房色播| 欧美巨大xxxx做受中文字幕| 久久dvd| av网站入口在线免费观看| 亚洲一区二区三区亚洲| 国产一区二区三区在线男友| 摸进她的内裤里疯狂揉她动图视频 | 国产无套粉嫩白浆在线观看| 好爽…又高潮了毛片免费看| 国产高清丝袜美腿视频在线观看| 亚洲精品一区二区三区麻豆| 美女网站免费观看视频| 欧美性大战久久久久久久| 完整在线视频免费黄片| 亚洲综合伊人久久综合|