999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

水量調(diào)度對內(nèi)秦淮河水質(zhì)改善的效應(yīng)評估

2015-02-22 07:01:52秦文浩
水資源保護(hù) 2015年2期
關(guān)鍵詞:水質(zhì)效果

夏 琨, 王 華, 秦文浩, 錢 鈞

(1.河海大學(xué)環(huán)境學(xué)院,江蘇 南京 210098; 2.江蘇省秦淮河水利工程管理處,江蘇 南京 210098)

水量調(diào)度對內(nèi)秦淮河水質(zhì)改善的效應(yīng)評估

夏 琨1, 王 華1, 秦文浩1, 錢 鈞2

(1.河海大學(xué)環(huán)境學(xué)院,江蘇 南京 210098; 2.江蘇省秦淮河水利工程管理處,江蘇 南京 210098)

為了分析論證內(nèi)秦淮河現(xiàn)狀引調(diào)水的效果和存在問題,采用MIKE11模型對河網(wǎng)進(jìn)行概化,建立了一維河網(wǎng)數(shù)學(xué)模型。以枯水期為例,選取6個控制斷面的氨氮濃度為水質(zhì)評價指標(biāo),運用驗證后的模型對內(nèi)秦淮河現(xiàn)狀引調(diào)水的效果進(jìn)行模擬與計算。結(jié)果表明:現(xiàn)狀條件下,內(nèi)秦淮河水體污染較嚴(yán)重,水源地水質(zhì)和引水流量達(dá)不到引水沖污的要求,需改善水源地水質(zhì),加大引水流量;每條引調(diào)水線路均存在引水盲區(qū),引水沖污應(yīng)采取4條線路聯(lián)合引調(diào)水的方式。

內(nèi)秦淮河;水量調(diào)度;引水沖污;一維河網(wǎng)數(shù)學(xué)模型;水質(zhì)改善效應(yīng);效應(yīng)評估

隨著城市人口的持續(xù)增長、城市經(jīng)濟(jì)的快速發(fā)展和城市化進(jìn)程的不斷推進(jìn),城市內(nèi)河在人類活動的影響下,水環(huán)境問題日益嚴(yán)重。而為了控制水位和防洪排澇,較多的閘控設(shè)施導(dǎo)致內(nèi)河的水體循環(huán)不暢,更加劇了水質(zhì)的惡化。時下,見效快、效果明顯,在全國范圍內(nèi)應(yīng)用廣泛[1-2]的城市內(nèi)河污染治理手段是引水沖污,即通過引調(diào)鄰近水域的清水到受污染水體來沖刷稀釋水體中污染物的濃度,以達(dá)到快速改善水體水環(huán)境的目的。目前,關(guān)于通過引水沖污改善城市內(nèi)河水環(huán)境的研究較多,如盧士強等[3]利用MIKE11模型模擬研究了蘇州河水系綜合調(diào)水對黃浦江以及蘊藻浜水質(zhì)的影響;盧士強等[4]利用MIKE11建立了上海市平原感潮河網(wǎng)的水動力和水質(zhì)模型,模擬和分析了主要調(diào)水方案對水質(zhì)的影響;劉云華等[5]對深圳河灣水系的珠江口引水和大鵬灣引水方案進(jìn)行了研究;陳志和等[6]研究了中山市水環(huán)境引水調(diào)控措施與實施效果;李梓嘉等[7]利用MIKE11模型模擬分析了各種引水沖污工程調(diào)控措施對泗洪縣城區(qū)水系水量水質(zhì)的影響;童朝鋒等[8]構(gòu)建了一維河網(wǎng)水動力和水質(zhì)模型模擬了不同時期不同水源補給以及在不同閘控系統(tǒng)運行調(diào)控方式下外秦淮河水質(zhì)的相應(yīng)變化情況。上述研究的重點都在達(dá)到相應(yīng)水質(zhì)目標(biāo)所需的引水流量和引水后水體水質(zhì)改善效果的認(rèn)定上,而基于現(xiàn)狀各單線引調(diào)水方案下的水質(zhì)情況和水源地情況,對現(xiàn)狀引調(diào)水的改善效果和存在問題的研究尚不多見。

自2007年起,江蘇省秦淮河水利工程管理處開始不定期引長江水沖洗河道,內(nèi)秦淮河水質(zhì)有所改善,但水質(zhì)時好時壞。據(jù)對內(nèi)秦淮河各段水質(zhì)評價的結(jié)果,2012年內(nèi)秦淮河全年水質(zhì)總體狀況為劣V類標(biāo)準(zhǔn),遠(yuǎn)低于Ⅳ類水體的功能區(qū)劃要求。內(nèi)秦淮河的主要污染物是有機物,表現(xiàn)為典型的城市污水污染特征,主要超標(biāo)項目有TN、NH3-N、DO等,其中TN、NH3-N是超標(biāo)最為嚴(yán)重的污染指標(biāo),超標(biāo)率在90%以上。筆者針對內(nèi)秦淮河面臨的水環(huán)境問題,基于現(xiàn)狀水量水質(zhì)監(jiān)測資料和引調(diào)水方案,建立MIKE11一維河網(wǎng)數(shù)學(xué)模型,利用驗證后的模型模擬現(xiàn)狀引水對內(nèi)秦淮河的水質(zhì)改善情況。以2012年枯水期(12月、1月和2月)為例,選擇秦淮河珍珠橋、竺橋、大中橋、五龍橋、文德橋和倉巷橋6個水質(zhì)控制斷面,選取主要超標(biāo)因子NH3-N為評價指標(biāo),定量評估現(xiàn)狀引水調(diào)度方案對內(nèi)秦淮河的水質(zhì)改善效果,分析內(nèi)秦淮河現(xiàn)狀引調(diào)水存在的問題,并提出相關(guān)建議,以期為內(nèi)秦淮河和類似城市內(nèi)河水系引水沖污的優(yōu)化提供借鑒。

1 研究區(qū)概況

內(nèi)秦淮河位于江蘇省南京市城南、城東地區(qū),處于北緯32°0′34.16″~32°4′0.07″、東經(jīng)118°45′25.77″~118°49′42.73″之間,沿途流經(jīng)秦淮、建鄴、白下、玄武4區(qū)及鼓樓區(qū)的一小部分,主流、支流全長25.2 km,匯水面積24.7 km2。內(nèi)秦淮河是主城區(qū)秦淮河水系內(nèi)最重要的景觀河道,穿越國家5A級夫子廟旅游風(fēng)景區(qū),主要功能是承接主城區(qū)秦淮河水系的洪水與澇水排入秦淮河干流。內(nèi)秦淮河主流分為南段、中段、東段及北段,支流包括進(jìn)香河、珍珠河、九華山溝、香林寺溝、清溪河、玉帶河及明御河等。內(nèi)秦淮河的地形呈東北高西南低,河水流向為自北向南、自東向西[9]。內(nèi)秦淮河地理位置及水系見圖1。

圖1 內(nèi)秦淮河地理位置

2 模型構(gòu)建

充分利用內(nèi)秦淮河水量水質(zhì)監(jiān)測數(shù)據(jù)和污染源調(diào)查資料,對內(nèi)秦淮河中段、南段、北段、東段、珍珠河、香林寺溝、明御河、玉帶河、清溪河等12條河段進(jìn)行概化,基于MIKE11平臺建立一維河網(wǎng)數(shù)學(xué)模型,并對模型進(jìn)行驗證,以確定主要河道的水量模型參數(shù)和水質(zhì)降解參數(shù)。

2.1 水量水質(zhì)耦合模型

計算水量的微分方程是建立在質(zhì)量和動量守恒定律基礎(chǔ)上的圣維南方程組,以流量和水位為未知變量,并補充考慮漫灘和旁側(cè)入流的完全形式的圣維南方程組為

(1)

式中:Q為流量;x為沿水流方向空間坐標(biāo);BW為調(diào)蓄寬度,指包括灘地在內(nèi)的全部河寬;Z為水位;t為時間;q為單位河長旁側(cè)入流流量,入流為正,出流為負(fù);u為斷面平均流速;g為重力加速度;A為過水?dāng)嗝婷娣e;B為主流斷面寬度;n為糙率;R為水力半徑。

對式(1)以普賴斯曼(Preissmann)四點隱式差分格式將其離散,輔以連接條件,形成河道方程,以微段、河段、汊點三級聯(lián)解的方法求解。

河網(wǎng)對流傳輸移動問題的河道和河道汊點基本方程為

(3)

式中:ρ為水流輸送的物質(zhì)質(zhì)量濃度;Ex為縱向分散系數(shù);Sc為與輸送物質(zhì)質(zhì)量濃度有關(guān)的衰減項,Sc=KdAρ,Kd為衰減因子;S為外部的源匯項;Ω為河道汊點—節(jié)點的水面面積;j為節(jié)點編號;i為與節(jié)點j相連接的河道編號。

對時間項采用向前差分、對流項采用上風(fēng)格式、擴(kuò)散項采用中心差分格式求解[10]。

2.2 河網(wǎng)概化

內(nèi)秦淮河河網(wǎng)可以概化成8個節(jié)點、12條河道,河網(wǎng)共6個邊界,其中4個入流流量控制邊界,包括作為玄武湖引水口的武廟閘、太平門閘、銅芯管閘,作為外秦淮河引水的象房村三孔閘;2個出流水位控制邊界,即鐵窗欞泵站和西水關(guān)泵站。內(nèi)秦淮河河網(wǎng)概化如圖2所示。

圖2 內(nèi)秦淮河河網(wǎng)概化圖(單位:km)

2.3 模型參數(shù)選擇及邊界條件選取

2.3.1 模型參數(shù)選擇

圖3 2012年內(nèi)秦淮河各控制斷面流量驗證結(jié)果

內(nèi)秦淮河糙率取為0.022~0.025,基于水質(zhì)監(jiān)測數(shù)據(jù)計算得到NH3-N的降解系數(shù),鼓樓區(qū)取0.023~0.025,玄武區(qū)取0.032~0.034,白下區(qū)取0.033~0.035,秦淮區(qū)取0.057~0.058。

2.3.2 模型邊界條件選取

根據(jù)內(nèi)秦淮河模型的概化結(jié)果,選取6個邊界斷面:武廟閘、太平門閘、銅芯管閘、象房村三孔閘、鐵窗欞泵站、西水關(guān)泵站。模型的邊界條件為以上6個邊界斷面各時段的流量、水位和水質(zhì)實測值。

2.4 模型驗證

2.4.1 水量模型驗證

水量模型驗證中,選取2012年為基準(zhǔn)年。基于2012年的6個邊界斷面的引調(diào)水資料,選取浮橋、竺橋、大中橋、倉巷橋、五龍橋和文德橋6個控制斷面進(jìn)行驗證,圖3為2012年內(nèi)秦淮河各斷面流量驗證結(jié)果。結(jié)果顯示相對誤差均在15%以內(nèi),計算值和實測值基本吻合,表明模型可以用來開展內(nèi)秦淮河水量的模擬和預(yù)測。

2.4.2 水質(zhì)模型驗證

水質(zhì)模型驗證中,選取2012年為基準(zhǔn)年。基于水質(zhì)監(jiān)測數(shù)據(jù)計算得到NH3-N的降解系數(shù),對浮橋、竺橋、大中橋、倉巷橋和五龍橋6個控制斷面的NH3-N質(zhì)量濃度進(jìn)行驗證。圖4為2012年內(nèi)秦淮河各斷面NH3-N質(zhì)量濃度驗證結(jié)果。結(jié)果表明模擬值和實測值吻合較好,所建模型可用來進(jìn)行內(nèi)秦淮河水質(zhì)的計算和預(yù)測。

3 水量調(diào)度改善效果評估

3.1 現(xiàn)狀引水調(diào)度

圖4 2012年內(nèi)秦淮河各控制斷面NH3-N質(zhì)量濃度驗證結(jié)果

內(nèi)秦淮河的引水水源有北源玄武湖、南源外秦淮河(武定門閘上游)和東源月牙湖,2012年枯水期3個水源地和6個控制斷面的NH3-N的質(zhì)量濃度情況見表1。現(xiàn)狀引調(diào)水采取常流量自流方式通過4條線路引水沖洗內(nèi)秦淮河(具體線路見圖5),其中線路一開武廟閘引玄武湖水約1 m3/s,線路二開太平門閘引玄武湖水約0.5 m3/s,線路三開象房村三孔閘引外秦淮河水約2 m3/s,線路四開銅芯管閘引月牙湖水約0.5 m3/s,4條線路均從鐵窗欞泵站和西水關(guān)泵站出水。

表1 2012年內(nèi)秦淮河枯水期水源地和控制斷面的NH3-N質(zhì)量濃度監(jiān)測值 mg/L

圖5 2012年內(nèi)秦淮河現(xiàn)狀引水調(diào)度線路

3.2 效果分析

運用所建一維河網(wǎng)數(shù)學(xué)模型,根據(jù)引調(diào)水實際工況中武廟閘、太平門閘、象房村三孔閘和銅芯管閘的引水流量,分別模擬和計算了枯水期4條線路現(xiàn)狀引調(diào)水流量下NH3-N質(zhì)量濃度的變化情況,具體計算結(jié)果見圖6。

圖6 2012年內(nèi)秦淮河各線路現(xiàn)狀引水流量下NH3-N質(zhì)量濃度變化情況

NH3-N的Ⅳ類標(biāo)準(zhǔn)為1.5 mg/L,Ⅴ類標(biāo)準(zhǔn)為2.0 mg/L,根據(jù)圖6的結(jié)果對引調(diào)水的改善效應(yīng)進(jìn)行評估與分析,結(jié)果表明:①在枯水期水源地現(xiàn)狀引水水質(zhì)和流量的條件下,內(nèi)秦淮河流域的NH3-N質(zhì)量濃度仍處于劣Ⅴ類,運高于其水環(huán)境功能要求。主要原因有:內(nèi)秦淮河污染較為嚴(yán)重,入河污染負(fù)荷較大;現(xiàn)狀引調(diào)水流量較小,引清治污能力有限,已不能滿足現(xiàn)狀水質(zhì)改善的要求;水源地均是間接引長江水沖洗內(nèi)秦淮河,長江水中NH3-N等的質(zhì)量濃度本來就偏高,水源地的NH3-N指標(biāo)除玄武湖較好外,外秦淮河和月牙湖的均為劣Ⅴ類,達(dá)不到水環(huán)境功能要求。②4條引調(diào)水線路分別存在改善區(qū)域和引水盲區(qū)。其中線路一的引調(diào)水對浮橋和大中橋斷面的改善效果較為明顯,NH3-N質(zhì)量濃度分別降低19%和12%,引水盲區(qū)為竺橋和大中橋斷面所在的城東水系;線路二的引調(diào)水對竺橋斷面的改善效果較為明顯,NH3-N質(zhì)量濃度降低11%,引水盲區(qū)為浮橋斷面所在的珍珠河和內(nèi)秦淮河北段;線路三的引調(diào)水對倉巷橋和文德橋斷面所在的內(nèi)秦淮河中段和南段的改善效果較為明顯,NH3-N質(zhì)量濃度分別降低18%和21%,引水盲區(qū)為浮橋、竺橋、五龍橋和大中橋斷面所在的河段;線路四的引調(diào)水對五龍橋的改善效果較為明顯,NH3-N質(zhì)量濃度降低11%,引水盲區(qū)為浮橋斷面所在的珍珠河和內(nèi)秦淮河北段。③對倉巷橋和文德橋斷面所在的重污染區(qū),北源和東源沖洗路徑過長,引調(diào)水對接近水源地河段的沖洗效果較好,對文德橋和倉巷橋斷面的水質(zhì)改善效果則不明顯,而南源沖洗路徑短,引調(diào)水對倉巷橋和文德橋斷面水質(zhì)改善效果最好,NH3-N質(zhì)量濃度降低幅度可達(dá)20%左右。

4 結(jié)論與建議

a. 現(xiàn)狀引調(diào)水條件下,引水流量不能滿足需求。各線路引水沖污對內(nèi)秦淮河水質(zhì)的改善效果不明顯,內(nèi)秦淮河流域水質(zhì)仍為劣Ⅴ類,需加大引調(diào)水流量;4條引水線路分別存在其引水改善效果較好區(qū)域和引水盲區(qū),引水沖污應(yīng)采用4條線路聯(lián)合引調(diào)水的方式;基于現(xiàn)狀引調(diào)水的不合理之處,應(yīng)進(jìn)一步研究改善秦淮河流域的水量水質(zhì)優(yōu)化調(diào)控方案。

b. 引水水源地的水質(zhì)需改善。水源地的水質(zhì)對引水工程改善內(nèi)秦淮河水環(huán)境有著至關(guān)重要的作用。干凈、可持續(xù)利用的水源是實施引水工程的首要前提條件,故需加大對水源地水環(huán)境的治理。當(dāng)前內(nèi)秦淮河引水水源有北源玄武湖、南源外秦淮河和東源月牙湖,這些水源地均受到不同程度的污染,尤其月牙湖接納了周邊娛樂場所和私人養(yǎng)殖所產(chǎn)生的污水,使得本線路引水沖污效果不明顯。

c. 其他整治措施力度需加大。內(nèi)秦淮河流域水體污染較嚴(yán)重,入河污染負(fù)荷較大,截污是內(nèi)秦淮河整治的根本性措施,只有從源頭上減少和控制污染物入河,才能保證內(nèi)秦淮河的水質(zhì)有改善的可能。引水沖污工程能緩解內(nèi)秦淮河水質(zhì)惡化,但治標(biāo)不治本,且排出的污水會成為外秦淮河的負(fù)擔(dān),故在引調(diào)清水的同時要實施其他整治措施,如截污、疏浚、曝氣復(fù)氧、生態(tài)修復(fù)等,內(nèi)秦淮河水環(huán)境的整治效果才會更理想。

[1] 徐貴泉,褚君達(dá).上海市引清調(diào)水改善水環(huán)境探討 [J].水資源保護(hù),2001(3):26-30.(XU Guiquan,CHU Junda.Water environment improvement by clean water diversion in Shanghai City [J].Water Resources Protection,2001(3):26-30.(in Chinese))

[2] 王超,衛(wèi)臻,張磊,等.平原河網(wǎng)區(qū)調(diào)水改善水環(huán)境試驗研究 [J].河海大學(xué)學(xué)報:自然科學(xué)版,2005,33(2):136-138.(WANG Chao,WEI Zhen,ZHANG Lei,et al.Experimental study on improvement of water environment by water diversion in plain river networks [J].Journal of Hohai University:Natural Sciences,2005,33(2):136-138.(in Chinese))

[3] 盧士強,徐祖信,林衛(wèi)青.綜合調(diào)水對蘇州河周邊水系水質(zhì)影響的數(shù)值模擬 [J].環(huán)境污染與防治,2005,27(6):473-475.(LU Shiqiang,XU Zuxin,LIN Weiqing.Numerical simulation of effect of Suzhou Creek flow augmentation on water quality of its ambient rivers [J].Environmental Pollution & Control,2005,27(6):473-475.(in Chinese))[4] 盧士強,徐祖信,羅海林,等.上海市主要河流調(diào)水方案的水質(zhì)影響分析 [J].河海大學(xué)學(xué)報:自然科學(xué)版,2006,34(l):32-36.(LU Shiqiang,XU Zuxin,LUO Hailin,et al.Influence of water diversion schemes on water quality of Shanghai river networks [J].Journal of Hohai University:Natural Sciences,2006,34(l):32-36.(in Chinese))[5] 劉云華,董增川,李朝方,等.深圳河流水系水質(zhì)改善引調(diào)水工程 [J].水資源保護(hù),2008,24(3):31-34.(LIU Yunhua,DONG Zengchuan,LI Chaofang,et al.Water diversion project for the improvement of water quality in Shenzhen River-Bay system [J].Water Resources Protection,2008,24(3):31-34.(in Chinese))

[6] 陳志和,陳曉宏,杜建,等.河網(wǎng)地區(qū)水環(huán)境引水調(diào)控及其效果預(yù)測 [J].水資源保護(hù),2012,28(3):16-21.(CHEN Zhi he,CHEN Xiaohong,DU Jian,et al.Water diversion and regulation of water environment and effect forecast in river network region [J].Water Resources Protection,2012,28(3):16-21.(in Chinese))

[7] 李梓嘉,董增川,樊孔明,等.MIKE11模型在泗洪縣城城區(qū)河網(wǎng)引水沖污工程中的應(yīng)用 [J].水電能源科學(xué),2012,30(8):100-103.(LI Zijia,DONG Zengchuan,FAN Kongming,et al.Application of MIKE11 model in water diversion and flushing pollutant of urban river network in Sihong City [J].Water Resources and Power,2012,30(8):100-103.(in Chinese))

[8] 童朝鋒,岳亮亮,郝嘉凌,等.南京市外秦淮河水質(zhì)模擬及引調(diào)水效果 [J].水資源保護(hù),2012,28(6):49-54.(TONG Chaofeng,YUE Liangliang,HAO Jialing,et al.Water quality simulation and water diversion effect analysis of external Qinhuai River in Nanjing [J].Water Resources Protection,2012,28(6):49-54.(in Chinese))

[9] 南京市地方志編纂委員會.南京環(huán)境保護(hù)志[M].北京:中國環(huán)境科學(xué)出版社,1996.

[10] 李曦,逄勇.江陰城區(qū)水環(huán)境數(shù)學(xué)模型建立及應(yīng)用 [J].水資源保護(hù),2009,25(2):38-41.(LI Xi,PANG Yong.Application of water environment numerical model in Jiangyin City [J].Water Resources Protection,2009,25(2):38-41.(in Chinese))

Assessment on effects of water quality improvement in Inner Qinhuai River by water dispatching

XIA Kun1, WANG Hua1, QIN Wenhao1, QIAN Jun2

(1.CollegeofEnvironment,HohaiUniversity,Nanjing210098,China;2.WaterConservancyProjectManagementOfficeoftheQinhuaiRiverinJiangsuProvince,Nanjing210098,China)

In order to analyze the effects and existing problems of water diversion in Inner Qinhuai River, a river network was generalized by MIKE11 model, and a one-dimension mathematical model of the river network was set up. Taking the typical dry season as an example, NH3-N concentration of six control sections were selected as the assessment index of water quality, and a validated model was applied to simulate and calculate the effects of water diversion. The results show that the Inner Qinhuai River is seriously polluted, and the water quality and the flow of water diversion are not up to the requirements of flushing pollutants, and the water quality in water source area should be improved and the flow of water diversion should be increased. There are blind spots of water diversion in every water diversion line, and we should take 4 lines of diverting water jointly to flush pollutants.

Inner Qinhuai River; water dispatching; water diversion for flushing pollutants; one-dimension mathematical model of river networks; effects of water quality improvement; effects assessment

10.3880/j.issn.1004-6933.2015.02.015

江蘇省水利科技項目(2012015)

夏琨(1991—),男,碩士研究生,研究方向為水環(huán)境保護(hù)及數(shù)值模擬。E-mail: xk19910523@163.com

X522

A

1004-6933(2015)02-0074-05

2014-05-06 編輯:彭桃英)

猜你喜歡
水質(zhì)效果
水質(zhì)抽檢豈容造假
環(huán)境(2023年5期)2023-06-30 01:20:01
按摩效果確有理論依據(jù)
一月冬棚養(yǎng)蝦常見水質(zhì)渾濁,要如何解決?這9大原因及處理方法你要知曉
這條魚供不應(yīng)求!蝦蟹養(yǎng)殖戶、垂釣者的最愛,不用投喂,還能凈化水質(zhì)
迅速制造慢門虛化效果
圖像識別在水質(zhì)檢測中的應(yīng)用
電子制作(2018年14期)2018-08-21 01:38:16
抓住“瞬間性”效果
中華詩詞(2018年11期)2018-03-26 06:41:34
模擬百種唇妝效果
Coco薇(2016年8期)2016-10-09 02:11:50
濟(jì)下水庫徑流水質(zhì)和垂向水質(zhì)分析及評價
3D—DSA與3D—CTA成像在顱內(nèi)動脈瘤早期診斷中的應(yīng)用效果比較
主站蜘蛛池模板: 国产在线精品99一区不卡| 91热爆在线| 久久国产精品夜色| 久久青草视频| 67194成是人免费无码| 国产精品19p| 国产成年无码AⅤ片在线| 美女裸体18禁网站| 亚洲码一区二区三区| 91在线激情在线观看| 黄片在线永久| 亚洲最黄视频| 尤物在线观看乱码| 国产极品美女在线播放| 久久无码高潮喷水| 熟妇丰满人妻av无码区| 欧美在线伊人| 精品国产免费观看一区| 日本高清免费不卡视频| 国产精品无码一区二区桃花视频| 欧美啪啪视频免码| 国产新AV天堂| 国产熟女一级毛片| 日本人妻丰满熟妇区| 国产在线自乱拍播放| 97视频在线精品国自产拍| 国产主播喷水| 亚洲国产成人综合精品2020| 全午夜免费一级毛片| 久久精品66| 久久99精品久久久久纯品| 99在线观看免费视频| 国产尤物jk自慰制服喷水| 亚洲AV无码不卡无码| 久久亚洲高清国产| 手机精品福利在线观看| 中文字幕av无码不卡免费| 波多野结衣一区二区三区AV| 91午夜福利在线观看| 亚洲制服中文字幕一区二区| 婷婷中文在线| 第一区免费在线观看| 国产18在线播放| 欧美怡红院视频一区二区三区| 天天干伊人| 91欧美在线| 欧美精品一区在线看| 91精品aⅴ无码中文字字幕蜜桃| 久久久久无码国产精品不卡| 手机精品视频在线观看免费| 黄色网在线| 国产精品久久久久久影院| 免费人成网站在线高清| 亚洲无线一二三四区男男| 亚洲日本一本dvd高清| av无码久久精品| 成人韩免费网站| 日韩精品一区二区三区免费在线观看| 成人免费视频一区二区三区| 欧美另类第一页| 国产成人在线无码免费视频| 久久黄色免费电影| 一级看片免费视频| 黄色在线不卡| 亚洲色图欧美视频| 国产草草影院18成年视频| 性视频一区| 欧美另类图片视频无弹跳第一页| 国产精品网址你懂的| 99精品一区二区免费视频| 91福利国产成人精品导航| 国产高清毛片| 2020国产精品视频| 日韩中文字幕亚洲无线码| 最新亚洲人成网站在线观看| 中文字幕无码中文字幕有码在线| 精品视频福利| yjizz国产在线视频网| 国产精品私拍在线爆乳| 99精品免费欧美成人小视频| 在线va视频| 国产激情在线视频|