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)用效果比較
主站蜘蛛池模板: 久久精品波多野结衣| www欧美在线观看| 免费网站成人亚洲| 日韩中文字幕免费在线观看| 国产91熟女高潮一区二区| 中文字幕第4页| 人人看人人鲁狠狠高清| 国产91在线|日本| 亚洲天堂网在线播放| 91在线高清视频| 就去色综合| 亚洲精品国产自在现线最新| 欧美精品导航| 国产综合精品日本亚洲777| 国产一区二区色淫影院| 欧美日韩中文国产| 国模视频一区二区| 国产丝袜一区二区三区视频免下载| 天堂在线视频精品| 最新国语自产精品视频在| 国产在线拍偷自揄拍精品| 狠狠色丁香婷婷| 免费va国产在线观看| 伊人五月丁香综合AⅤ| 日韩精品成人网页视频在线 | 3344在线观看无码| 免费毛片全部不收费的| 国产在线高清一级毛片| 国产亚洲高清视频| 国产福利微拍精品一区二区| 国产乱子伦精品视频| 亚洲国产成人超福利久久精品| 亚洲区一区| 国产性猛交XXXX免费看| 日本手机在线视频| 国产精品漂亮美女在线观看| 波多野结衣一区二区三区四区 | 老司国产精品视频91| 国产人在线成免费视频| 欧美在线中文字幕| 午夜毛片免费观看视频 | 三上悠亚精品二区在线观看| 最新亚洲av女人的天堂| 国产一区二区三区在线观看视频 | 欧美精品1区2区| 黄色福利在线| 国产精品蜜芽在线观看| 国产精品jizz在线观看软件| 亚洲欧美一区二区三区图片| 99久久精品视香蕉蕉| 19国产精品麻豆免费观看| 久久久国产精品免费视频| 亚洲男人的天堂在线观看| 国产精品亚洲专区一区| 欧美色伊人| 久久99国产精品成人欧美| 91亚洲视频下载| 中文字幕有乳无码| 国产无码精品在线| 91黄色在线观看| 波多野结衣视频网站| 欧美日韩久久综合| 欧美亚洲国产视频| 欧美成人第一页| 丝袜亚洲综合| 成人免费黄色小视频| 亚洲第一成年人网站| 美女无遮挡拍拍拍免费视频| 免费人成网站在线观看欧美| 久久国产亚洲欧美日韩精品| 91久久夜色精品国产网站| 国产一级毛片在线| 国产精品毛片一区| 五月婷婷导航| 国内精品视频| 国产第一页免费浮力影院| 无码AV动漫| 一本大道视频精品人妻| 国产成人一区免费观看| 亚洲人成网线在线播放va| 亚洲精品片911| 日本高清在线看免费观看|