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

    沙潁河安徽段水環(huán)境容量計(jì)算及動(dòng)態(tài)分析

    2019-09-10 19:32:39王瑩張幼寬梁修雨謝顯傳
    人民黃河 2019年1期
    關(guān)鍵詞:沙潁河

    王瑩 張幼寬 梁修雨 謝顯傳

    摘要:基于對(duì)沙潁河安徽段污染負(fù)荷的調(diào)查以及水質(zhì)和流量監(jiān)測(cè)數(shù)據(jù)的分析,采用MIKE11模型建立水流與水質(zhì)模型,采用線性規(guī)劃法,以入淮河干流Ⅲ類水為水質(zhì)目標(biāo),計(jì)算沙潁河安徽段水環(huán)境容量并分析其動(dòng)態(tài)特征。結(jié)果表明:沙潁河安徽段CODM。和氨氮年水環(huán)境容量分別為2.39萬(wàn)、0.76萬(wàn)t/a:受流量及上游水質(zhì)影響,COD和氨氮水環(huán)境容量月變化較大,8月水環(huán)境容量最大,分別為3 557、1363 t,3月水環(huán)境容量最小,分別為621、181 t;上游水質(zhì)對(duì)沙潁河安徽段水環(huán)境容量影響顯著。

    關(guān)鍵詞:水環(huán)境容量;水質(zhì)模型;線性規(guī)劃法;沙潁河

    中圖分類號(hào):X522

    文獻(xiàn)標(biāo)志碼:A

    doi:10. 3969/j .issn. 1000- 13 79.2019.01. 018

    水環(huán)境容量是指水體環(huán)境在規(guī)定的環(huán)境目標(biāo)下所能容納的污染負(fù)荷[1]。水環(huán)境容量與水體特征、水質(zhì)目標(biāo)及污染物特性有關(guān),同時(shí)還與污染物的排放方式及排放的時(shí)空分布密切相關(guān)[2],是容量總量控制的核心內(nèi)容之一[3]。水環(huán)境容量計(jì)算的基本思路:通過(guò)解析或數(shù)值等方法計(jì)算河流的水流過(guò)程及污染物遷移轉(zhuǎn)化規(guī)律,建立污染負(fù)荷排放與水質(zhì)響應(yīng)的關(guān)系,基于給定的環(huán)境目標(biāo)(如水質(zhì)目標(biāo))計(jì)算指定區(qū)域所能容納的最大污染負(fù)荷。對(duì)于特定的河流,在水質(zhì)目標(biāo)不變的情況下,由于水文特征、自凈能力和排污口排污方式的時(shí)空變異性,因此水環(huán)境容量是隨時(shí)間和空間變化的。流域?qū)嶋H的水環(huán)境容量需要針對(duì)實(shí)際背景條件計(jì)算,即計(jì)算動(dòng)態(tài)水環(huán)境容量[4]。水環(huán)境容量的計(jì)算方法主要包括公式法、模型試錯(cuò)法、系統(tǒng)最優(yōu)化法(如線性規(guī)劃法和隨機(jī)規(guī)劃法)、概率稀釋模型法和模糊數(shù)學(xué)法[5-6],其中系統(tǒng)最優(yōu)化法適用范圍廣,對(duì)邊界條件和設(shè)計(jì)流量的適應(yīng)能力較強(qiáng),是目前水環(huán)境容量計(jì)算的主流方法之一[7]。

    本研究以沙潁河流域安徽段為例,采用系統(tǒng)最優(yōu)化法中的線性規(guī)劃法計(jì)算水環(huán)境容量。整體思路:根據(jù)排污數(shù)據(jù)和水功能區(qū)劃確定排污口和控制斷面的位置:通過(guò)河流水流和水質(zhì)數(shù)值模擬,建立排污口與控制斷面的水質(zhì)響應(yīng)關(guān)系,計(jì)算實(shí)際本底濃度值和各排污口負(fù)荷貢獻(xiàn)度:建立線性最優(yōu)化模型,計(jì)算流域水環(huán)境容量并分析其動(dòng)態(tài)特征。

    1 研究區(qū)概況

    沙潁河是淮河最大的一級(jí)支流,全長(zhǎng)624 km,流域面積39 980 km2。沙潁河流域年平均氣溫約15℃,年平均降水量約700 mm/a。沙潁河流域跨越河南和安徽兩省,在安徽境內(nèi)長(zhǎng)208 km,流域面積4% km2,占沙潁河流域總面積的10%。沙潁河在安徽境內(nèi)被稱為潁河,其主要支流為汾泉河(見(jiàn)圖1)。雖然沙潁河流域僅占淮河流域總面積的1/7,但其污染負(fù)荷排放量占淮河流域的1/3。近年來(lái),隨著工業(yè)發(fā)展和人口增長(zhǎng),沙潁河水質(zhì)日趨惡化,污染嚴(yán)重,水質(zhì)劣于V類的水功能區(qū)占總數(shù)的60%[8]。為控制流域水污染,實(shí)現(xiàn)沙潁河人淮河水質(zhì)達(dá)到Ⅲ類水的目標(biāo),需要制定科學(xué)合理的污染控制方案,而河流的納污能力(水環(huán)境容量)是制定該方案的基礎(chǔ)。

    沙潁河安徽段的主要污染物來(lái)自城鎮(zhèn)生活廢水排放和工業(yè)點(diǎn)源污染。根據(jù)2012年安徽省環(huán)境統(tǒng)計(jì)數(shù)據(jù),該段共有排污口約166個(gè),其空間分布及排污負(fù)荷見(jiàn)圖2。

    為計(jì)算水環(huán)境容量,采用重心概化法將干流沿岸的排污口概化為6個(gè)排污口。

    2 數(shù)值模擬

    2.1 水流模擬

    河流水流過(guò)程是水質(zhì)模擬和水環(huán)境容量計(jì)算的基礎(chǔ),采用數(shù)值方法模擬水流過(guò)程。研究?jī)H考慮污染物沿河流流向的一維遷移過(guò)程[9],因此選用一維水動(dòng)力模型MIKE11的HD模塊模擬河流徑流量。MIKE11水動(dòng)力計(jì)算模型是基于一維非恒定流Saint -Venant方程組來(lái)模擬河流或河口的水流狀態(tài).Saint - Venant方程組為

    對(duì)式(1)采用Abbott -Ionescu六點(diǎn)隱式有限差分格式求解。剖分網(wǎng)格尺寸最大值為5 000 m,并在局部加密,共生成179個(gè)計(jì)算節(jié)點(diǎn)。時(shí)間步長(zhǎng)為5 min,總計(jì)算時(shí)段為2 a(2011-01-01-2012-12-31)。以界首閘和沈丘閘實(shí)測(cè)日流量為上游流量邊界,以潁上閘實(shí)測(cè)日河水位為下游水位邊界,初始河水位通過(guò)上下游觀測(cè)水位線性插值獲得。采用阜陽(yáng)閘和潁上閘的實(shí)測(cè)日流量與水位進(jìn)行模型率定驗(yàn)證,率定期為2011年,驗(yàn)證期為2012年。根據(jù)MIKE11模型按過(guò)水?dāng)嗝嬖O(shè)置糙率的方法,將斷面分為上、中、下三部分,率定參數(shù)為斷面上、中、下三部分的曼寧糙率n。研究區(qū)河道斷面變化較大,為有效校正模型[10].斷面三部分的初始n值均為0.02[11]?;谟^測(cè)流量和水位,通過(guò)調(diào)整參數(shù),獲得最佳擬合參數(shù)。

    2012年阜陽(yáng)閘和潁上閘流量與水位的模擬值與觀測(cè)值見(jiàn)圖3、圖4.其中:2012年阜陽(yáng)閘流量的相對(duì)誤差為11.4%(Nash系數(shù)為0.83),水位的相對(duì)誤差為5.3%(Nash系數(shù)為0.93);潁上閘流量的相對(duì)誤差為12.3%(Nash系數(shù)為0.75),水位的相對(duì)誤差為0.3%(Nash系數(shù)為0.97)。模擬結(jié)果基本符合實(shí)際水文情況。通過(guò)驗(yàn)證,獲得斷面上、中、下三部分的n值分別為0.020,0.022和0.025。驗(yàn)證的水流模型可用于河流COD和氨氮運(yùn)移模擬以及水環(huán)境容量計(jì)算。

    2.2 水質(zhì)模擬

    污染物在河流中遷移的控制方程為一維對(duì)流彌散方程,考慮一級(jí)降解,其方程為

    式(3)主要參數(shù)是彌敢度a和降解系數(shù)k,筆者參考沙潁河流域已有研究成果[12],選取合適的彌散度a,見(jiàn)表1。姚瑞珍[13]通過(guò)室內(nèi)試驗(yàn)方法確定了CODM。和氨氮在河流中的初始降解系數(shù)k分別為0.053 6和0.034 9,并基于實(shí)測(cè)楊湖斷面水質(zhì),通過(guò)敏感度分析及手動(dòng)調(diào)整參數(shù),獲得了最佳擬合參數(shù),見(jiàn)表1。2011年楊湖斷面水質(zhì)模擬值和實(shí)測(cè)值見(jiàn)圖5,其中COD模擬的相對(duì)誤差為28%,氨氮模擬的相對(duì)誤差為10%.水質(zhì)模型基本符合實(shí)際情況,可用于沙潁河安徽段的水環(huán)境容量計(jì)算。

    數(shù)值模型網(wǎng)格剖分方式及尺寸與水流模型相同,計(jì)算節(jié)點(diǎn)位置與水流交叉分布。模擬時(shí)間步長(zhǎng)與計(jì)算時(shí)段亦與水流模型設(shè)置相同。以往水環(huán)境容量計(jì)算以水質(zhì)目標(biāo)值代替上游來(lái)水水質(zhì)的真實(shí)值,忽略了上游輸入的污染負(fù)荷對(duì)計(jì)算河段環(huán)境容量的影響。研究區(qū)上游為豫皖省界斷面,該斷面水質(zhì)的變化對(duì)安徽段水環(huán)境容量有不可忽略的影響,因此以上游豫皖省界斷面每周觀測(cè)的水質(zhì)為模型上游邊界的濃度進(jìn)行計(jì)算。

    為獲取研究區(qū)河流的背景值(考慮上游污染負(fù)荷輸入).以上游豫皖省界斷面實(shí)測(cè)水質(zhì)為模型上游邊界,假設(shè)區(qū)內(nèi)排污量為0,計(jì)算2011-2012年河流污染物濃度。由于研究區(qū)實(shí)測(cè)水質(zhì)資料缺乏,無(wú)法給定符合實(shí)際的初始濃度,因此采用數(shù)值模擬方法估算初始條件影響范圍,以此消除初始條件對(duì)水質(zhì)模擬的影響。具體方法為,假設(shè)4組COD和氨氮的初始濃度值(C0)分別為1、5、10、15 mg/L。采用上述邊界條件及參數(shù),計(jì)算人淮口濃度時(shí)間序列(見(jiàn)圖6)。從圖6可以看出,當(dāng)COD和氨氮模擬時(shí)間為50 d時(shí),初始條件的影響可忽略(濃度差別小于0.03 mg/L)。因此選取第51 d的模擬濃度值為初始條件,作為河流背景濃度值的估算值,用于河流氨氮和COD的遷移轉(zhuǎn)化過(guò)程模擬。

    3 水環(huán)境容量計(jì)算

    針對(duì)水功能區(qū)劃及水質(zhì)目標(biāo)[8],劃定3個(gè)水環(huán)境容量計(jì)算單元及相應(yīng)的控制斷面(見(jiàn)圖2),控制斷面的水質(zhì)目標(biāo)均為Ⅲ類。以往的水環(huán)境容量研究大都采用單一的水文設(shè)計(jì)條件,而根據(jù)不同的水文條件計(jì)算出的結(jié)果差別較大,利用最枯流量計(jì)算的結(jié)果往往嚴(yán)重低于實(shí)際環(huán)境容量值[14]。為了解實(shí)際動(dòng)態(tài)水環(huán)境容量動(dòng)態(tài)特征,從2010-2014年日觀測(cè)數(shù)據(jù)中選取2012年作為水文設(shè)計(jì)條件,通過(guò)線性規(guī)劃模型及水質(zhì)模擬,計(jì)算出各排污口允許最大日排污量及河流日環(huán)境容量,累加可得排污口允許最大年排污量和河流月環(huán)境容量。

    由于對(duì)流彌散方程是線性的,具有可疊加性,因此可以把多個(gè)排污口對(duì)控制斷面水質(zhì)的綜合影響劃分為單個(gè)排污口對(duì)其影響的線性疊加[15]?;诏B加原理,李適宇等[16]提出貢獻(xiàn)度的概念,結(jié)合線性規(guī)劃法求解環(huán)境容量。該方法的思路:設(shè)某個(gè)排污口(i)單位時(shí)間排放單位負(fù)荷,而其他排污口無(wú)負(fù)荷量排放,背景濃度值為0,采用一維水質(zhì)模型計(jì)算各水質(zhì)控制斷面(j)的濃度值,該濃度值即為i排污口對(duì)各斷面的污染貢獻(xiàn)度(a);改變排污口,重復(fù)以上步驟,即可求出每個(gè)排污口的污染貢獻(xiàn)度:基于排污口的污染貢獻(xiàn)度,考慮河流背景濃度值,以控制斷面水質(zhì)目標(biāo)為約束,建立線性規(guī)劃模型的約束條件,確定目標(biāo)函數(shù),即總排污負(fù)荷最大(環(huán)境容量)。具體目標(biāo)函數(shù)為

    4 結(jié)果與討論

    4.1 排污口環(huán)境容量

    采用MIKE11模型計(jì)算出各排污口對(duì)各控制斷面COD和氨氮的貢獻(xiàn)度,年均日貢獻(xiàn)度見(jiàn)表2、表3(其中0表示排污口不在此斷面控制范圍內(nèi))??梢钥闯?,排污口1#、2#、3#對(duì)控制斷面I的貢獻(xiàn)度基本相等;僅排污口6#對(duì)控制斷面Ⅱ有貢獻(xiàn);所有排污口對(duì)控制斷面Ⅲ均有貢獻(xiàn),其中排污口5#貢獻(xiàn)度最大,排污口6#貢獻(xiàn)度最小。排污口離控制斷面越近,貢獻(xiàn)度越大,即排污口對(duì)控制斷面的影響隨距離的減小而增大。

    基于排污口對(duì)各斷面的貢獻(xiàn)度,利用線性規(guī)劃模型計(jì)算出2012年各排污口的環(huán)境容量,結(jié)果見(jiàn)圖7。沙潁河安徽段COD環(huán)境容量約2.39萬(wàn)t/a,氨氮環(huán)境容量約0.76萬(wàn)t/a。其中:排污口5#的COD和氨氮容量最大,分別為7 892、2 970 t/a,占總水環(huán)境容量的38.90/0和32.9%;排污口3#的COD和氨氮環(huán)境容量最小,分別為948、165 t/a,僅占總水環(huán)境容量的3.9%和2.2%。原因是:排污口5#主要受斷面Ⅲ控制,對(duì)斷面I和斷面Ⅱ的貢獻(xiàn)度為0,受到的限制較小;而排污口3#對(duì)斷面I和斷面Ⅲ的貢獻(xiàn)度較大,受到的限制較大。排污口4#、5#、6#的環(huán)境容量較大,而排污口1#、2#、3#的環(huán)境容量較小,原因是上游來(lái)水(省界斷面)COD和氨氮濃度較大,大部分時(shí)期超過(guò)Ⅲ類水質(zhì)標(biāo)準(zhǔn),余留的環(huán)境容量小,而隨水流過(guò)程中污染物的降解和稀釋,下游的環(huán)境容量逐漸增大。

    4.2 月環(huán)境容量變化分析

    2012年沙潁河安徽段月水環(huán)境容量及河流月徑流量見(jiàn)圖8。從圖8可知,COD和氨氮水環(huán)境容量年內(nèi)變幅較大,豐水期(7-10月)環(huán)境容量較大,枯水期(11月至翌年5月)環(huán)境容量較小。COD和氨氮環(huán)境容量均在8月最大,分別為3 557、1363 t;COD和氨氮環(huán)境容量在3月最小,分別為621、181 t;氨氮和COD環(huán)境容量的變化趨勢(shì)與徑流量不完全一致,1-7月徑流量逐漸減少,而COD和氨氮環(huán)境容量逐漸增大,直到8月隨徑流量增大而達(dá)到頂峰,說(shuō)明水環(huán)境容量受徑流量、上游來(lái)水水質(zhì)等因素的綜合影響。

    4.3 省界水質(zhì)對(duì)環(huán)境容量的影響

    為分析省界斷面水質(zhì)對(duì)沙潁河安徽段環(huán)境容量的影響,設(shè)定3種情景討論環(huán)境容量的變化。情景1為省界斷面COD和氨氮濃度為0(即無(wú)負(fù)荷補(bǔ)給)、背景CODM。和氨氮濃度為0,計(jì)算水環(huán)境容量(L0);情景2為省界斷面水質(zhì)為Ⅲ類水(水質(zhì)目標(biāo)),計(jì)算水環(huán)境容量(L);情景3為省界斷面COD氨氮濃度為實(shí)測(cè)值(實(shí)際情況),計(jì)算COD、氨氮環(huán)境容量(L)。

    3種情景月水環(huán)境容量見(jiàn)圖9。在無(wú)上游污染負(fù)荷輸入的情景1下,環(huán)境容量最大,COD、氨氮年容量分別為36.0萬(wàn)、9.1萬(wàn)t:上游污染負(fù)荷達(dá)標(biāo)排放情景下,環(huán)境容量較小,河段的環(huán)境容量主要來(lái)自CODM。和氨氮在河流中天然降解量.COD、氨氮年容量分別為17.5萬(wàn)、5.5萬(wàn)t;情景3,環(huán)境容量最小,COD、氨氮年容量分別為1.96萬(wàn)、0.62萬(wàn)t.水環(huán)境容量?jī)H約為情景2的9%。COD和氨氮的環(huán)境容量均為L(zhǎng)0>L>L,說(shuō)明上游水質(zhì)越差,水環(huán)境容量越小。

    2012年豫皖省界斷面水質(zhì)情況見(jiàn)圖10.由圖10可知,COD的超標(biāo)率在泉河省控?cái)嗝孢_(dá)到55%,在潁河省控?cái)嗝孢_(dá)到29%:氨氮的超標(biāo)率在泉河省控?cái)嗝孢_(dá)到18%.在潁河省控?cái)嗝孢_(dá)到46%,水質(zhì)達(dá)標(biāo)率較低。其中,COD超標(biāo)時(shí)間主要為1-6月,氨氮超標(biāo)時(shí)間主要為1-6月(潁河省控?cái)嗝妫┖?-9月(泉河省控?cái)嗝妫?。由此可?jiàn),水環(huán)境容量除受河流流量控制外,還受上游水質(zhì)的影響。

    5 結(jié)論

    基于河流流量及水質(zhì)觀測(cè)數(shù)據(jù),利用MIKE11水流及水質(zhì)模型,結(jié)合線性規(guī)劃法,計(jì)算沙潁河流域安徽段水環(huán)境容量,并分析水環(huán)境容量影響因素,結(jié)果表明:MIKE11模型對(duì)水文和水質(zhì)的模擬基本符合實(shí)際水文情況,模擬結(jié)果與實(shí)測(cè)數(shù)據(jù)基本一致;沙潁河安徽段COD№、氨氮年環(huán)境容量總和約為2.39萬(wàn)、0.76萬(wàn)t/a,COD和氨氮水環(huán)境容量年內(nèi)變幅較大,8月容量最大,分別為3 557、1 363 t,3月容量最小,分別為621、181 t:上游水質(zhì)對(duì)沙潁河安徽段水環(huán)境容量影響較大,實(shí)際水環(huán)境容量?jī)H為上游水質(zhì)達(dá)標(biāo)(Ⅲ類水)前提下理論水環(huán)境容量的9%。

    參考文獻(xiàn):

    [1] 張玉清,張?zhí)N華,張景霞,河流功能區(qū)水污染物容量總量控制的原理和方法[J].環(huán)境科學(xué),1998,19(增刊1):23-33.

    [2] 張永良,水環(huán)境容量基本概念的發(fā)展[J].環(huán)境科學(xué)研究,1992,5(3):59-61.

    [3]孟偉,張楠,張遠(yuǎn),等,流域水質(zhì)目標(biāo)管理技術(shù)研究(I):控制單元的總量控制技術(shù)[J].環(huán)境科學(xué)研究,2007,20 (4):1—8.

    [4]姜欣,許士國(guó),練建軍,等,北方河流動(dòng)態(tài)水環(huán)境容量分析與計(jì)算[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2013,29(4):409-414.

    [5]鄭孝宇,褚君達(dá),朱維斌,河網(wǎng)非穩(wěn)態(tài)水環(huán)境容量研究[J].水科學(xué)進(jìn)展,1997,8(1):25-31.

    [6]周孝德,郭瑾瓏,程文,等,水環(huán)境容量計(jì)算方法研究[J].西安理工大學(xué)學(xué)報(bào),1999,15(3):1-6.

    [7] 董飛,劉曉波,彭文啟,等,地表水水環(huán)境容量計(jì)算方法回顧與展望[J].水科學(xué)進(jìn)展,2014,25(3):451-463.

    [8]郝守寧,潁河流域水功能區(qū)面源污染控制研究[D].鄭州:華北水利水電大學(xué),2014:1-2.

    [9]逄勇,陸桂華,水環(huán)境容量計(jì)算理論及應(yīng)用[M].北京:科學(xué)出版社,2010:4.

    [10] 中國(guó)環(huán)境規(guī)劃院,全國(guó)水環(huán)境容量核定技術(shù)指南[R].北京:中國(guó)環(huán)境規(guī)劃院,2003:46-76.

    [11] 陳煉鋼,施勇,錢新,等,閘控河網(wǎng)水文一水動(dòng)力一水質(zhì)耦合數(shù)學(xué)模型:Ⅱ應(yīng)用[J].水科學(xué)進(jìn)展,2014,25(6):856-863.

    [12] 劉玉年,施勇,程緒水,等,準(zhǔn)河中游水量水質(zhì)聯(lián)合調(diào)度模型研究[J].水科學(xué)進(jìn)展,2009,20(2):177-183.

    [13] 姚瑞珍,人工控制的半封閉河流水環(huán)境容量測(cè)算研究:以老府河為實(shí)例[D].武漢:華中科技大學(xué),2006:34-43.

    [14] 李紅亮,王樹(shù)峰,不同設(shè)計(jì)水文條件下河北省水功能區(qū)納污能力研究[J].南水北調(diào)與水利科技,2010,8(3):68-70.

    [15] LIS Y, TOHRU M. Optimal Allocation of Waste Loadsin aRiver with Probabilistic Tributary Flow Under TransverseMixing[ J]. Water Environment Research, 1999, 71(2): 156-162.

    [16] 李適宇,李耀初,陳炳祿,等,分區(qū)達(dá)標(biāo)控制法求解海域環(huán)境容量[J].環(huán)境科學(xué),1999,20(4):97-100.

    猜你喜歡
    沙潁河
    在沙潁河畔
    沙潁河聽(tīng)水
    沙潁河周口段氮形態(tài)的時(shí)空變化及高濃度成因解析
    水資源與經(jīng)濟(jì)社會(huì)發(fā)展時(shí)空匹配特征量化分析
    人民黃河(2022年2期)2022-02-13 11:50:04
    閘壩控制下沙潁河沉積物重金屬污染特征與生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)
    沙潁河流域管理
    淺析“城市雙修”視角下的濱水公園設(shè)計(jì)——以界首沙潁河公園為例
    北京園林(2020年1期)2020-01-18 05:19:28
    沙潁河耿樓復(fù)線船閘通航水位分析
    河南省沙潁河岸線評(píng)價(jià)及布局
    沙潁河流域管理
    亚洲aⅴ乱码一区二区在线播放 | av在线天堂中文字幕| 亚洲黑人精品在线| 少妇的丰满在线观看| 国内精品久久久久精免费| 特级一级黄色大片| 女人高潮潮喷娇喘18禁视频| 变态另类丝袜制服| 男女视频在线观看网站免费 | 宅男免费午夜| 欧美在线一区亚洲| 久久久久国产一级毛片高清牌| 一区二区三区国产精品乱码| 午夜福利高清视频| 色综合站精品国产| 国产真实乱freesex| 国产aⅴ精品一区二区三区波| 亚洲色图av天堂| 日韩成人在线观看一区二区三区| 少妇裸体淫交视频免费看高清 | 国产视频一区二区在线看| 国产主播在线观看一区二区| 亚洲av熟女| 99热6这里只有精品| 亚洲 欧美一区二区三区| 小说图片视频综合网站| 在线免费观看的www视频| 亚洲国产精品成人综合色| 久久精品国产99精品国产亚洲性色| 成人精品一区二区免费| 好男人电影高清在线观看| 在线免费观看的www视频| 久久精品成人免费网站| 久久国产精品人妻蜜桃| 在线观看免费视频日本深夜| 丰满的人妻完整版| 成人手机av| 久久午夜综合久久蜜桃| 精品国产超薄肉色丝袜足j| 午夜免费观看网址| 久久久久九九精品影院| 久久久久久久久免费视频了| 成人特级黄色片久久久久久久| 久久99热这里只有精品18| 成人午夜高清在线视频| xxxwww97欧美| 亚洲av中文字字幕乱码综合| 亚洲真实伦在线观看| 视频区欧美日本亚洲| 1024手机看黄色片| 看黄色毛片网站| 男人舔女人的私密视频| 我要搜黄色片| 欧洲精品卡2卡3卡4卡5卡区| 国产乱人伦免费视频| 99re在线观看精品视频| 欧美在线一区亚洲| 丰满人妻熟妇乱又伦精品不卡| 激情在线观看视频在线高清| 亚洲美女视频黄频| 久久久久国内视频| 久久国产精品人妻蜜桃| 变态另类成人亚洲欧美熟女| 免费在线观看完整版高清| av超薄肉色丝袜交足视频| e午夜精品久久久久久久| 99精品在免费线老司机午夜| 亚洲国产高清在线一区二区三| 男女做爰动态图高潮gif福利片| 天天添夜夜摸| 99热6这里只有精品| 国产成人av激情在线播放| 日本黄大片高清| 亚洲国产看品久久| 精品久久久久久久末码| 成人18禁高潮啪啪吃奶动态图| 亚洲专区国产一区二区| 亚洲美女视频黄频| 啪啪无遮挡十八禁网站| 正在播放国产对白刺激| 亚洲中文字幕日韩| 长腿黑丝高跟| 99久久精品国产亚洲精品| 搡老熟女国产l中国老女人| 99国产极品粉嫩在线观看| 九九热线精品视视频播放| 久久九九热精品免费| 午夜久久久久精精品| 亚洲色图av天堂| 男人舔女人下体高潮全视频| 这个男人来自地球电影免费观看| 国产日本99.免费观看| 亚洲av中文字字幕乱码综合| 在线视频色国产色| 久久人妻福利社区极品人妻图片| 少妇的丰满在线观看| 超碰成人久久| 精品国产美女av久久久久小说| 18美女黄网站色大片免费观看| 国产久久久一区二区三区| 午夜激情av网站| 国产又色又爽无遮挡免费看| a级毛片在线看网站| 欧美国产日韩亚洲一区| 岛国在线观看网站| 欧美黑人欧美精品刺激| 欧美乱妇无乱码| 美女扒开内裤让男人捅视频| 国产午夜福利久久久久久| 岛国视频午夜一区免费看| 亚洲成av人片免费观看| 亚洲精品av麻豆狂野| 久久精品人妻少妇| 亚洲av成人av| 99久久国产精品久久久| 好看av亚洲va欧美ⅴa在| 亚洲成人精品中文字幕电影| 一边摸一边抽搐一进一小说| 久久久久久久久免费视频了| av天堂在线播放| 嫩草影院精品99| 99久久精品国产亚洲精品| 一本综合久久免费| 久久99热这里只有精品18| 日韩大码丰满熟妇| 麻豆一二三区av精品| 日韩精品青青久久久久久| 国产精品精品国产色婷婷| xxxwww97欧美| 亚洲最大成人中文| 久久中文看片网| 国产精华一区二区三区| 脱女人内裤的视频| 午夜福利欧美成人| 久久性视频一级片| 久久久久免费精品人妻一区二区| 欧美黑人精品巨大| 欧美成人一区二区免费高清观看 | 老司机在亚洲福利影院| 国产精品久久久久久精品电影| 麻豆av在线久日| 国产人伦9x9x在线观看| 日本黄大片高清| 国产三级黄色录像| 最新在线观看一区二区三区| 在线观看66精品国产| 午夜免费观看网址| 老司机午夜十八禁免费视频| 好看av亚洲va欧美ⅴa在| 精华霜和精华液先用哪个| 99国产精品一区二区蜜桃av| 久久久精品大字幕| 久久人人精品亚洲av| 成人午夜高清在线视频| 一本精品99久久精品77| 伊人久久大香线蕉亚洲五| 国产精品精品国产色婷婷| 99热6这里只有精品| 搞女人的毛片| 人人妻,人人澡人人爽秒播| or卡值多少钱| 天堂影院成人在线观看| 精品人妻1区二区| 国产黄色小视频在线观看| 亚洲avbb在线观看| 一进一出抽搐gif免费好疼| 国产精品精品国产色婷婷| 特大巨黑吊av在线直播| 人人妻人人看人人澡| 日韩免费av在线播放| 99久久久亚洲精品蜜臀av| 日本撒尿小便嘘嘘汇集6| 麻豆国产av国片精品| 狂野欧美激情性xxxx| 亚洲男人的天堂狠狠| 日日爽夜夜爽网站| 国产精品98久久久久久宅男小说| 国产精品影院久久| 国产精品香港三级国产av潘金莲| 国产av不卡久久| 丝袜美腿诱惑在线| 好男人在线观看高清免费视频| 免费看日本二区| 国产av麻豆久久久久久久| 在线观看66精品国产| 欧美av亚洲av综合av国产av| 哪里可以看免费的av片| 夜夜爽天天搞| 日本撒尿小便嘘嘘汇集6| 免费在线观看日本一区| 国产精品亚洲一级av第二区| 久久精品夜夜夜夜夜久久蜜豆 | 欧美3d第一页| 在线观看66精品国产| 亚洲精品美女久久久久99蜜臀| 国产区一区二久久| 亚洲成人免费电影在线观看| 久久久精品欧美日韩精品| 亚洲免费av在线视频| 日日夜夜操网爽| 夜夜看夜夜爽夜夜摸| 国产高清视频在线观看网站| 欧美中文日本在线观看视频| 亚洲中文字幕一区二区三区有码在线看 | 女同久久另类99精品国产91| 欧美性猛交黑人性爽| 老司机午夜福利在线观看视频| 手机成人av网站| 男女床上黄色一级片免费看| 可以免费在线观看a视频的电影网站| 久久久久国产精品人妻aⅴ院| 丰满人妻一区二区三区视频av | 亚洲中文av在线| 国内精品久久久久久久电影| 成人精品一区二区免费| 悠悠久久av| 巨乳人妻的诱惑在线观看| 国内精品久久久久精免费| 精品国产乱子伦一区二区三区| 免费在线观看视频国产中文字幕亚洲| 国产亚洲精品久久久久久毛片| 12—13女人毛片做爰片一| 欧美午夜高清在线| 在线观看午夜福利视频| 国产成人啪精品午夜网站| 一本大道久久a久久精品| 日本五十路高清| 91成年电影在线观看| 午夜老司机福利片| 看黄色毛片网站| 亚洲av成人不卡在线观看播放网| 国产av不卡久久| 成年人黄色毛片网站| 国产成人av激情在线播放| 免费在线观看影片大全网站| 亚洲av美国av| 香蕉国产在线看| 亚洲人成网站高清观看| 国产三级在线视频| bbb黄色大片| 精品福利观看| 亚洲精品中文字幕一二三四区| 国内精品久久久久精免费| 欧美在线黄色| 女警被强在线播放| 欧美三级亚洲精品| 中亚洲国语对白在线视频| 99久久99久久久精品蜜桃| 丰满的人妻完整版| 亚洲avbb在线观看| 婷婷精品国产亚洲av在线| 这个男人来自地球电影免费观看| 婷婷亚洲欧美| 男女之事视频高清在线观看| 色综合婷婷激情| 国产精品久久久久久精品电影| 99久久国产精品久久久| 国产一级毛片七仙女欲春2| 亚洲欧美一区二区三区黑人| 日韩欧美在线乱码| 国产一区二区激情短视频| 欧美国产日韩亚洲一区| 午夜福利欧美成人| 国产亚洲欧美98| 午夜两性在线视频| 51午夜福利影视在线观看| 亚洲在线自拍视频| 久久婷婷人人爽人人干人人爱| 18禁黄网站禁片午夜丰满| 激情在线观看视频在线高清| 国内毛片毛片毛片毛片毛片| 淫妇啪啪啪对白视频| 亚洲国产精品sss在线观看| 后天国语完整版免费观看| 少妇熟女aⅴ在线视频| 国产午夜精品久久久久久| 熟女少妇亚洲综合色aaa.| 好男人电影高清在线观看| 亚洲国产欧美一区二区综合| 亚洲自偷自拍图片 自拍| 五月玫瑰六月丁香| 国产高清videossex| 人人妻,人人澡人人爽秒播| 久久久久九九精品影院| 午夜激情av网站| 亚洲狠狠婷婷综合久久图片| 成人欧美大片| 宅男免费午夜| 级片在线观看| 国产一区二区在线观看日韩 | 一区福利在线观看| 午夜福利成人在线免费观看| 午夜福利18| 日本一区二区免费在线视频| 少妇粗大呻吟视频| 男人的好看免费观看在线视频 | 午夜精品在线福利| 国产亚洲av高清不卡| 制服丝袜大香蕉在线| 可以在线观看毛片的网站| 国产欧美日韩一区二区三| 免费看日本二区| 大型黄色视频在线免费观看| 99在线视频只有这里精品首页| 国产精品98久久久久久宅男小说| 99久久无色码亚洲精品果冻| 免费在线观看视频国产中文字幕亚洲| 亚洲人成77777在线视频| 国产99久久九九免费精品| 久久99热这里只有精品18| 亚洲欧美精品综合久久99| 国产av在哪里看| 操出白浆在线播放| 日本免费一区二区三区高清不卡| 99在线视频只有这里精品首页| 一边摸一边做爽爽视频免费| 色播亚洲综合网| 99国产精品一区二区三区| 色综合站精品国产| 中文字幕av在线有码专区| 长腿黑丝高跟| 啦啦啦韩国在线观看视频| 午夜免费观看网址| 中文在线观看免费www的网站 | 亚洲人成网站高清观看| 亚洲av美国av| 国内久久婷婷六月综合欲色啪| 少妇人妻一区二区三区视频| 可以免费在线观看a视频的电影网站| 神马国产精品三级电影在线观看 | 亚洲成人精品中文字幕电影| 操出白浆在线播放| 白带黄色成豆腐渣| 99在线人妻在线中文字幕| 欧美成人午夜精品| 白带黄色成豆腐渣| 成人三级黄色视频| 久久久国产精品麻豆| 大型av网站在线播放| av片东京热男人的天堂| 精品第一国产精品| 午夜免费成人在线视频| 免费搜索国产男女视频| 国内精品久久久久精免费| 亚洲自拍偷在线| 青草久久国产| 国产亚洲精品一区二区www| av欧美777| 国产精品一区二区三区四区免费观看 | 两人在一起打扑克的视频| 变态另类成人亚洲欧美熟女| 久久久久久国产a免费观看| 成人三级做爰电影| 国产乱人伦免费视频| 国产精品免费一区二区三区在线| 给我免费播放毛片高清在线观看| 日韩中文字幕欧美一区二区| 久久久久九九精品影院| 国产午夜福利久久久久久| 嫁个100分男人电影在线观看| av福利片在线观看| a在线观看视频网站| 国产午夜精品论理片| 美女高潮喷水抽搐中文字幕| 哪里可以看免费的av片| 在线免费观看的www视频| 国产成人影院久久av| 人妻丰满熟妇av一区二区三区| 午夜福利高清视频| 中亚洲国语对白在线视频| 亚洲中文日韩欧美视频| 久久这里只有精品中国| bbb黄色大片| 777久久人妻少妇嫩草av网站| 欧美日韩中文字幕国产精品一区二区三区| 午夜成年电影在线免费观看| 国产男靠女视频免费网站| 日韩高清综合在线| 日本熟妇午夜| 亚洲精品国产精品久久久不卡| 欧美日韩亚洲国产一区二区在线观看| av福利片在线观看| 777久久人妻少妇嫩草av网站| 麻豆国产97在线/欧美 | 一个人免费在线观看的高清视频| xxx96com| 久久精品国产亚洲av香蕉五月| 久久这里只有精品19| av免费在线观看网站| 国产一区二区三区在线臀色熟女| 怎么达到女性高潮| 在线永久观看黄色视频| 99久久综合精品五月天人人| 亚洲欧美精品综合一区二区三区| 国内毛片毛片毛片毛片毛片| 国产精品电影一区二区三区| 黄色女人牲交| 1024视频免费在线观看| 午夜精品一区二区三区免费看| 欧美又色又爽又黄视频| 国产精品久久久久久人妻精品电影| 一级毛片女人18水好多| 国产高清激情床上av| 精品无人区乱码1区二区| 久久香蕉精品热| 精品免费久久久久久久清纯| 欧美一级a爱片免费观看看 | 一卡2卡三卡四卡精品乱码亚洲| 久久久久久九九精品二区国产 | 精品不卡国产一区二区三区| 少妇裸体淫交视频免费看高清 | 久久欧美精品欧美久久欧美| 在线观看午夜福利视频| 老司机深夜福利视频在线观看| 久久久水蜜桃国产精品网| 国产单亲对白刺激| 亚洲国产精品999在线| 特大巨黑吊av在线直播| 丰满的人妻完整版| 国产探花在线观看一区二区| videosex国产| 亚洲专区字幕在线| av超薄肉色丝袜交足视频| 国产亚洲av高清不卡| 欧美乱妇无乱码| 国产爱豆传媒在线观看 | av在线播放免费不卡| 欧美另类亚洲清纯唯美| 在线免费观看的www视频| 久久久久久久久免费视频了| 久久久国产精品麻豆| 搡老熟女国产l中国老女人| 亚洲av第一区精品v没综合| 欧美日韩精品网址| 黄片大片在线免费观看| 亚洲一区二区三区色噜噜| 麻豆成人午夜福利视频| 亚洲国产日韩欧美精品在线观看 | 亚洲精品久久成人aⅴ小说| 亚洲欧洲精品一区二区精品久久久| 国产成年人精品一区二区| 一级毛片女人18水好多| 极品教师在线免费播放| 白带黄色成豆腐渣| 久久午夜亚洲精品久久| 又粗又爽又猛毛片免费看| aaaaa片日本免费| 18禁黄网站禁片免费观看直播| 亚洲avbb在线观看| 99精品欧美一区二区三区四区| 精品国内亚洲2022精品成人| 最近视频中文字幕2019在线8| 亚洲欧美日韩东京热| 最近最新免费中文字幕在线| 三级国产精品欧美在线观看 | 国产1区2区3区精品| 欧美日韩亚洲综合一区二区三区_| 国产精品国产高清国产av| 一个人观看的视频www高清免费观看 | 黑人操中国人逼视频| 黄频高清免费视频| 日本成人三级电影网站| 国产熟女午夜一区二区三区| 日本五十路高清| 欧美日韩亚洲国产一区二区在线观看| 99热只有精品国产| 国产亚洲精品综合一区在线观看 | 成人特级黄色片久久久久久久| av福利片在线| 久久久国产欧美日韩av| 国产精品 欧美亚洲| 亚洲男人天堂网一区| 成人欧美大片| 我要搜黄色片| 夜夜看夜夜爽夜夜摸| 热99re8久久精品国产| 中文字幕人成人乱码亚洲影| 欧美久久黑人一区二区| av在线天堂中文字幕| 国产精品一区二区精品视频观看| 欧美色视频一区免费| 人人妻人人澡欧美一区二区| 国产成人啪精品午夜网站| 亚洲精品美女久久久久99蜜臀| 欧美乱码精品一区二区三区| 男女视频在线观看网站免费 | 国产区一区二久久| 免费av毛片视频| 窝窝影院91人妻| 精品免费久久久久久久清纯| 精品一区二区三区视频在线观看免费| 免费在线观看亚洲国产| 久久99热这里只有精品18| 亚洲精品av麻豆狂野| 亚洲精品中文字幕一二三四区| 日本黄色视频三级网站网址| 免费在线观看亚洲国产| 国产精品久久久久久精品电影| 免费一级毛片在线播放高清视频| 欧美zozozo另类| 在线观看www视频免费| av片东京热男人的天堂| 18禁国产床啪视频网站| 亚洲av成人一区二区三| 麻豆av在线久日| 亚洲 国产 在线| 美女扒开内裤让男人捅视频| 19禁男女啪啪无遮挡网站| 天天一区二区日本电影三级| 亚洲色图 男人天堂 中文字幕| 久久国产乱子伦精品免费另类| 国产av一区在线观看免费| 老司机深夜福利视频在线观看| 久久精品成人免费网站| 国产精品影院久久| 久久久久精品国产欧美久久久| 亚洲av日韩精品久久久久久密| 亚洲av五月六月丁香网| 麻豆成人午夜福利视频| 国产不卡一卡二| 欧美一级a爱片免费观看看 | 亚洲人成网站在线播放欧美日韩| 精品无人区乱码1区二区| 久久精品国产亚洲av香蕉五月| 国产1区2区3区精品| 免费人成视频x8x8入口观看| 欧美一区二区精品小视频在线| 成人亚洲精品av一区二区| 亚洲熟妇熟女久久| 欧美不卡视频在线免费观看 | 午夜福利在线在线| 天堂影院成人在线观看| 丝袜美腿诱惑在线| 亚洲精品在线美女| 国产v大片淫在线免费观看| 欧美黄色片欧美黄色片| 欧美日韩国产亚洲二区| 男人的好看免费观看在线视频 | 一本一本综合久久| 久久久久久人人人人人| 免费在线观看亚洲国产| 久久这里只有精品19| 免费无遮挡裸体视频| 亚洲国产欧美人成| 在线看三级毛片| 在线播放国产精品三级| 观看免费一级毛片| 男人舔女人下体高潮全视频| 我要搜黄色片| 看免费av毛片| 国产精品久久久av美女十八| 91麻豆av在线| 亚洲精品久久成人aⅴ小说| 久久久久精品国产欧美久久久| av中文乱码字幕在线| 后天国语完整版免费观看| 亚洲 国产 在线| 欧美丝袜亚洲另类 | 在线观看www视频免费| 久久久精品欧美日韩精品| 熟妇人妻久久中文字幕3abv| 日韩欧美一区二区三区在线观看| 级片在线观看| 久久精品亚洲精品国产色婷小说| 19禁男女啪啪无遮挡网站| 国产片内射在线| 国产精品综合久久久久久久免费| 日本一区二区免费在线视频| 777久久人妻少妇嫩草av网站| 亚洲人成伊人成综合网2020| 国产精品久久久av美女十八| 天堂影院成人在线观看| АⅤ资源中文在线天堂| 国产97色在线日韩免费| 成熟少妇高潮喷水视频| 2021天堂中文幕一二区在线观| 欧美成人一区二区免费高清观看 | 国产亚洲精品第一综合不卡| 免费看十八禁软件| 在线观看免费视频日本深夜| 18禁美女被吸乳视频| 国产乱人伦免费视频| 国产精品精品国产色婷婷| 亚洲国产精品合色在线| 一边摸一边抽搐一进一小说| 亚洲专区字幕在线| 最近最新免费中文字幕在线| 欧美性长视频在线观看| 久久久久久久精品吃奶| 亚洲熟女毛片儿| 可以在线观看的亚洲视频| 日韩欧美国产在线观看| 在线观看舔阴道视频| 亚洲av成人av| 国产探花在线观看一区二区| 久久中文看片网| 9191精品国产免费久久| 黄频高清免费视频| 看免费av毛片| 伊人久久大香线蕉亚洲五| 久久久久久大精品| 午夜成年电影在线免费观看| 亚洲精品久久成人aⅴ小说| 国产精品一区二区免费欧美| 久久精品国产亚洲av香蕉五月| 欧美国产日韩亚洲一区| 久久99热这里只有精品18| 日日摸夜夜添夜夜添小说| 18禁美女被吸乳视频| 欧美日韩国产亚洲二区| 国产精品久久久久久人妻精品电影| 国产精品免费视频内射| 国产精品久久久久久人妻精品电影| 成人特级黄色片久久久久久久| bbb黄色大片|