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

    基于熵權(quán)云模型評價法的省級考核斷面水質(zhì)評價與對比

    2022-05-27 01:58:08馬聰聰劉登峰楊志勇潘保柱
    人民珠江 2022年5期
    關(guān)鍵詞:變權(quán)歐式水質(zhì)評價

    馬聰聰,劉登峰*,黃 強(qiáng),楊志勇,潘保柱,林 木

    (1.西安理工大學(xué) 省部共建西北旱區(qū)生態(tài)水利國家重點實驗室, 陜西 西安 710048;2.中國水利水電科學(xué)研究院,北京 100038;3.中央財經(jīng)大學(xué)統(tǒng)計與數(shù)學(xué)學(xué)院,北京 100081)

    水質(zhì)評價是合理開發(fā)、利用和保護(hù)水資源和生態(tài)環(huán)境的一項基本工作[1-2]。自從Jacobs提出水質(zhì)指數(shù)的概念和公式以來,廣大專家學(xué)者關(guān)于水質(zhì)評價的方法與思路不斷出現(xiàn),填補了水質(zhì)評價方面的空缺[3]。

    中國在河流水質(zhì)綜合評價的研究已有40多年的歷史,提出的評價模型及其方法有二三十種。由于水質(zhì)綜合評價的特殊性,不同方法對同一數(shù)據(jù)會有不同的結(jié)果,因此沒有一個被大家公認(rèn)的、具有可比性的數(shù)學(xué)模型[4]。方案評價問題是廣泛存在的一類問題,除了水質(zhì)評價,在水文模型參數(shù)優(yōu)化[5]、工程規(guī)劃設(shè)計、水量調(diào)配[6]、水庫調(diào)度方案優(yōu)選[7-8]等領(lǐng)域也廣泛存在,方法多樣。在水質(zhì)評價中,單因子評價法操作簡單、使用廣泛,但由于其污染程度最高的評價指標(biāo)作為最后的評價結(jié)果,無視了其他指標(biāo)對水體的影響作用,不能全面反映水體污染的情況[9-10];模糊綜合評價法需構(gòu)造隸屬函數(shù)矩陣和模糊關(guān)系矩陣,通過運算得出結(jié)論。其結(jié)果更加符合客觀的規(guī)律,具有較強(qiáng)的合理性,是到目前為止水質(zhì)評價中應(yīng)用最廣的方法。但由于其概念模糊,計算復(fù)雜,在多評價斷面、多評價指標(biāo)體系中存在過程冗長、結(jié)果容易出錯的問題,影響其使用與推廣[11-13]。綜合指數(shù)法是先確定污染指標(biāo)的權(quán)重,然后進(jìn)行評價的方法,陳仁杰等[14]建立了一種適用于集中式生活飲用水、地表水的改良綜合指數(shù)法,具有簡單、合理和可操作性強(qiáng)的優(yōu)點,可以應(yīng)用于其他城市的水質(zhì)評價;BP神經(jīng)網(wǎng)絡(luò)評價法無需確定權(quán)重,適用于地區(qū)間評價,具有精度高的優(yōu)點,但是其原理和計算過程較為復(fù)雜[15];集對分析法復(fù)合了決策和統(tǒng)計學(xué)理論,它精度高,更能反映水質(zhì)變化的實際情況,但其計算方法相對復(fù)雜[16];灰色關(guān)聯(lián)法打破了傳統(tǒng)的精確數(shù)學(xué)的束縛,因此,在水質(zhì)評價中的使用率極高,但計算過程十分繁瑣,計算起來比較復(fù)雜和困難[17];熵權(quán)云模型綜合評價法中,熵權(quán)法在計算權(quán)重的過程中靈活性不夠,云模型又可以實現(xiàn)定量數(shù)據(jù)和定性概念之間的轉(zhuǎn)換,兩者結(jié)合起來可以客觀地反映實際情況[18];張慶慶等[19]使用變權(quán)歐式距離模型評價錢塘江支流的水質(zhì)進(jìn)行綜合評價,驗證了變權(quán)歐式距離模型應(yīng)用于水質(zhì)綜合評價的科學(xué)性和合理性。新的評價方法不斷出現(xiàn),對評價方法的對比研究仍然是研究的重要內(nèi)容。

    單因子評價法是實際過程中最常用到的水質(zhì)評價方法。該方法可以幫助快速確定水體中的主要污染因子;變權(quán)歐式距離法全面考慮了各污染因子,而且強(qiáng)調(diào)了超標(biāo)污染物對水質(zhì)的影響,便于對水質(zhì)進(jìn)行因子分析,適用于進(jìn)一步對水質(zhì)因子的變化規(guī)律進(jìn)行分析;熵權(quán)云模型綜合評價法計算簡單高效,還能有效地削弱界限值對最終評價結(jié)果的影響,適用于對水質(zhì)的綜合評價中取得一個較易接受的結(jié)果。因此,以2020年陜西省渭河、漢江、丹江3條主要河流的水質(zhì)監(jiān)測斷面數(shù)據(jù),選擇熵權(quán)云模型綜合評價法、單因子評價法和變權(quán)歐式距離法,評價并對比水質(zhì)評價的結(jié)果,研究結(jié)果為河長制考核和生態(tài)保護(hù)的評價提供參考。

    1 數(shù)據(jù)與方法

    1.1 數(shù)據(jù)來源

    選取陜西省渭河、漢江、丹江考核斷面為評價對象,對2020年監(jiān)測數(shù)據(jù)進(jìn)行評價,監(jiān)測數(shù)據(jù)來自陜西省生態(tài)環(huán)境廳(http://sthjt.shaanxi.gov.cn/shjzl.html),水質(zhì)指標(biāo)有化學(xué)需氧量(COD)、氨氮(NH3-N)、溶解氧(DO)、總磷(TP)。水質(zhì)監(jiān)測斷面序號1—40號(其中,1—27號位于渭河干支流,28—38號位于漢江干支流,39—40號位于丹江干流),見表1、圖1。水質(zhì)監(jiān)測項目的等級根據(jù)GB 3838—2002《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》進(jìn)行劃分。因陜西省生態(tài)廳在漢江、丹江流域內(nèi)不考核溶解氧與總磷指標(biāo),因此后續(xù)分析中漢江、丹江流域的污染因子少于渭河流域的污染因子。

    表1 水質(zhì)監(jiān)測點分布信息

    圖1 研究區(qū)域及水質(zhì)監(jiān)測點分布

    1.2 熵權(quán)云模型綜合評價法

    熵權(quán)云模型綜合評價法的計算步驟如下:

    a)利用熵權(quán)法[20]計算各水質(zhì)指標(biāo)的權(quán)重矩陣W=wi。熵權(quán)法是在沒有專家權(quán)重的情況下,為了減少主觀因素的影響,根據(jù)各指標(biāo)的變異程度,利用信息熵計算各指標(biāo)的熵權(quán),通過評價值構(gòu)成的判斷矩陣來確定評價指標(biāo)的權(quán)重。設(shè)評價對象有m個(i=1,2,…,m),評價指標(biāo)有n個(j=1,2,…,n),得到原始矩陣A=(aij)m×n采用min-max標(biāo)準(zhǔn)化方法對指標(biāo)值進(jìn)行歸一化處理,得到指標(biāo)值的歸一化矩陣R=(rij)m×n。計算各個指標(biāo)的熵值。將第j項指標(biāo)的熵Ej定義為:

    (1)

    (2)

    式中,如果yij=0時,yijlnyij=0。

    計算各指標(biāo)權(quán)重Wj。

    (3)

    b)利用X條件云發(fā)生器求解水質(zhì)指標(biāo)的隸屬度。根據(jù)GB 3838—2002《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》確定水質(zhì)評價標(biāo)準(zhǔn)的等級標(biāo)注。計算不同等級的特征值以及不同指標(biāo)對應(yīng)的隸屬度。

    (4)

    (5)

    Heij=Enij×0.1

    (6)

    式中aijmax、aijmin——第i個指標(biāo)在第j個等級中的上下限;Exij、Enij——第i個指標(biāo)在第j個等級中的期望與熵;Heij——超熵。

    利用X條件云發(fā)生器,輸入各個標(biāo)準(zhǔn)值,重復(fù)計算1 000次,得出每個等級的隸屬度,最后輸出隸屬度矩陣Z。

    c)根據(jù)已經(jīng)求出的權(quán)重矩陣W和隸屬度矩陣Z,由矩陣相乘計算得到各個評價等級的確定度。根據(jù)最大隸屬度原則,對應(yīng)隸屬度最大值的水質(zhì)評價等級即為該水質(zhì)監(jiān)測斷面的評價級別。

    1.3 變權(quán)歐式距離法

    建立以各個污染因子為坐標(biāo)軸的n維空間,以水質(zhì)標(biāo)準(zhǔn)中的I級濃度限制為“原點”,然后結(jié)合污染因子的權(quán)重系數(shù)計算每個待評價單元以及水質(zhì)評價標(biāo)準(zhǔn)中其他各級別水質(zhì)與“原點”間的變權(quán)歐式距離,并將后者作為劃分每個待評價單元水質(zhì)級別的標(biāo)準(zhǔn)[21]。計算步驟如下。

    a)建立歐式空間。設(shè)實測樣本為A=(aji)n×m(i=1,2,…,m;j=1,2,…,n),水質(zhì)標(biāo)準(zhǔn)矩陣B=(bik)m×5(i=1,2,…,m;k=1,2,…,5)。將A中的每一行和B中的每一列均看作以n個污染因子為坐標(biāo)軸的歐式空間中的點。

    b)歐式空間坐標(biāo)變換。將水質(zhì)標(biāo)準(zhǔn)矩陣B的Ⅰ級濃度限值作為歐式空間的原點,則A和B的坐標(biāo)變?yōu)椋?/p>

    (7)

    (8)

    c)權(quán)重確定[22]。根據(jù)污染貢獻(xiàn)率計算各個監(jiān)測值的權(quán)重,即某監(jiān)測值超標(biāo)倍數(shù)越多,其權(quán)重(wji)越大:

    (9)

    對于像DO這樣濃度越高越好的因子,權(quán)重(wjλ)為:

    (10)

    d)變權(quán)歐式距離的確定[23]。相對等級距離(Djk)確定:

    (11)

    相對實測值距離(dj)確定:

    (12)

    e)評價等級(Rj)的確定。當(dāng)Djk≤dj≤Dj(k-1)時,第j個評價單元的具體等級為:

    (13)

    2 結(jié)果與分析

    2.1 單因子評價法

    對40個斷面2020年1—6月份的監(jiān)測結(jié)果進(jìn)行評價,1、2月基于單因子評價法的評價結(jié)果見圖2、3。其中:寶雞峽總干渠(斷面序號1)、小韋河(斷面序號7)、新河(斷面序號13)、皂河(斷面序號16)、臨河(斷面序號23)、漕運明渠(斷面序號22)、幸福渠(斷面序號19)、沋河(斷面序號25)等8條基本無生態(tài)水的支流干渠不考核溶解氧指標(biāo),漢江流域、丹江流域(斷面序號28—40)不考核溶解氧、總磷指標(biāo)。1月寶雞峽總干渠(斷面序號1)未監(jiān)測。2月臨河、漕運明渠、幸福渠、涇河(斷面序號15)、金錢河(斷面序號38)和丹江商洛出境(斷面序號40)等6個斷面受新冠肺炎疫情影響,未監(jiān)測;渭河寶雞出境(斷面序號5)、渭河咸陽出境(斷面序號17)、小韋河、皂河、沋河等5個自動站總磷未監(jiān)測。5月寶雞峽總干渠(斷面序號1)未監(jiān)測。表2為1—6月監(jiān)測斷面水質(zhì)類別單因子評價結(jié)果進(jìn)行統(tǒng)計。

    圖2 1月各監(jiān)測斷面單因子評價結(jié)果

    圖3 2月各監(jiān)測斷面單因子評價結(jié)果

    表2 監(jiān)測斷面水質(zhì)類別單因子評價結(jié)果統(tǒng)計

    單因子評價顯示,1月水質(zhì)較差,只有65%左右的水質(zhì)達(dá)到了III類標(biāo)準(zhǔn),化學(xué)需氧量、氨氮和總磷超標(biāo)情況比較嚴(yán)重,2—6月的水質(zhì)達(dá)標(biāo)率在75%以上,化學(xué)需氧量和氨氮的超標(biāo)情況較為嚴(yán)重。1—6月的污染出現(xiàn)在渭河支流的小韋河、新河、皂河、幸福渠、灞河、臨河、沋河、北洛河、沋河等9個斷面出現(xiàn)Ⅳ、Ⅴ類水質(zhì)的概率較大。

    2.2 熵權(quán)云模型綜合評價法

    a)用熵權(quán)法得到的1月指標(biāo)權(quán)重見表3。

    表3 1月各水質(zhì)指標(biāo)權(quán)重

    b)計算云模型特征參數(shù)及不同評價等級的正態(tài)隸屬度云圖,隸屬度云圖見圖4。

    a)化學(xué)需氧量

    b)氨氮

    c)溶解氧

    以1月渭河干流寶雞出境為例,計算指標(biāo)隸屬度,結(jié)果見表4。

    由表3的權(quán)重和表4的隸屬度,計算水質(zhì)隸屬度,并且根據(jù)最大隸屬度原則確定寶雞出境評價結(jié)果,評價結(jié)果見表5。表6為1—6月監(jiān)測斷面水質(zhì)類別熵權(quán)云模型綜合評價結(jié)果。

    表5 1月渭河干流寶雞出境斷面熵權(quán)云模型綜合評價結(jié)果

    表6 監(jiān)測斷面水質(zhì)類別的熵權(quán)云模型綜合評價結(jié)果

    由表6可以看出,1—6月的Ⅲ類水占比全為75%以上。由熵權(quán)重的計算結(jié)果可以看出,總磷和氨氮對綜合評價的影響較大,即其為主要污染物質(zhì)。1—6月在渭河支流的新河、皂河、臨河和沋河等4個斷面出現(xiàn)Ⅳ、Ⅴ類水質(zhì)的概率較大。

    2.3 變權(quán)歐式距離法

    應(yīng)用式(7)將污染因子進(jìn)行坐標(biāo)轉(zhuǎn)化,根據(jù)污染物的貢獻(xiàn)率,計算其權(quán)重(式(9)、(10))。對每個單元的水質(zhì)進(jìn)行評價(式(13)),進(jìn)而得到評價結(jié)果,監(jiān)測斷面水質(zhì)類別變權(quán)歐式距離法評價結(jié)果統(tǒng)計,見表7。

    表7 監(jiān)測斷面水質(zhì)類別變權(quán)歐式距離法評價結(jié)果

    由表7可以看出,1—6月的Ⅲ類水占比全為80%以上,甚至4、6月占比達(dá)到90%。1—6月在渭河支流的新河、皂河、幸福渠、臨河和沋河等5個斷面出現(xiàn)Ⅳ、Ⅴ類水質(zhì)的概率較大。表8為1月5個斷面各污染因子權(quán)重,表9為1月5個斷面水質(zhì)評價結(jié)果。

    表8 1月5個斷面各污染因子權(quán)重

    表9 1月5個斷面變權(quán)歐式距離法評價結(jié)果

    歐式距離超過1的評價單元,其水質(zhì)級別為Ⅴ類,級別中的小數(shù),代表了水質(zhì)結(jié)果偏離水質(zhì)標(biāo)準(zhǔn)的程度,如幸福渠的級別為4.732 5,表明水質(zhì)介于Ⅳ—Ⅴ類,但更偏向于Ⅴ類。

    2.4 討論

    圖5為1—6月評價結(jié)果對比。從圖中可以看出,3種評價方法結(jié)果略有差異。其中:熵權(quán)云模型綜合評價法和變權(quán)歐式距離法的評價結(jié)果更為相近,在1—6月的評價結(jié)果中,Ⅲ類及以上的評價結(jié)果占比相近,都在84%之上,表明2種方法對水質(zhì)的評價結(jié)果更為相近,驗證了2種方法用于水質(zhì)評價的合理性和科學(xué)性。熵權(quán)云模型綜合評價法通過熵權(quán)對實際數(shù)據(jù)進(jìn)行修正,得到了客觀的指標(biāo)權(quán)重,但是其權(quán)重計算過程更加依賴數(shù)據(jù)本身,忽略了污染因子本身的重要程度,無法像專家打分法更多地考慮主觀因素。變權(quán)歐式距離法在考慮各個污染物因子的基礎(chǔ)上,同時考慮了污染物的貢獻(xiàn)率,可以與單因子評價一起確定主要的污染因子。與其余2種方法相比,它既可得出評價斷面的水質(zhì)類別,還可得到評價斷面對不同級別水質(zhì)標(biāo)準(zhǔn)的距離。

    d)2月渭河流域

    g)3月渭河流域

    j)4月渭河流域

    m)5月渭河流域

    p)6月渭河流域

    漢江流域和丹江流域因其污染因子監(jiān)測較少,對已有的監(jiān)測結(jié)果進(jìn)行水質(zhì)分析可知,3種評價結(jié)果相似,且水質(zhì)等級均處于Ⅰ—Ⅱ類,水質(zhì)情況良好,因此主要分析渭河流域。

    渭河支流中的新河、皂河、幸福渠、臨河和沋河5條支流是污染較為嚴(yán)重的支流,主要位于西安和咸陽境內(nèi)。如1月的幸福渠,單因子評價法表明NH3-N濃度最高,根據(jù)單因子評價法的判定,導(dǎo)致了幸福渠的水質(zhì)等級為Ⅴ類,產(chǎn)生了“過保護(hù)”的評價結(jié)果;在熵權(quán)云模型綜合評價法中,通過對各個評價指標(biāo)的權(quán)重計算,綜合考慮了不同指標(biāo)對水質(zhì)等級的影響,雖然NH3-N的濃度超過Ⅴ類水的評價標(biāo)準(zhǔn)很多,但是DO和COD的濃度均在Ⅰ—Ⅱ類,熵權(quán)云模型綜合評價法權(quán)衡了DO和COD的較優(yōu)的表現(xiàn),使得水質(zhì)的整體評價結(jié)果為Ⅲ類;在變權(quán)歐式距離法中,以評價因子的實際污染貢獻(xiàn)率作為各個評價單元的權(quán)重,同時給出了各個評價單元的不同等級的劃分標(biāo)準(zhǔn),NH3-N的權(quán)重為0.448 2,明顯高于其他3項污染因子,因而NH3-N是幸福渠的主要污染因子,水質(zhì)級別為4.732 5,表明水質(zhì)處于Ⅳ~Ⅴ類,水質(zhì)評價結(jié)果為Ⅳ類,因主要污染因子的含量較高,在一定程度上影響了其他表現(xiàn)較好的污染因子,從而使得評價結(jié)果較差。

    渭河流域陜西段的污染范圍主要集中于西安和咸陽兩市,污染源主要是生活污水、農(nóng)業(yè)的面源污染和工業(yè)廢水排放。監(jiān)測結(jié)果可以表明,渭河氮磷污染較為嚴(yán)重,而渭河的氮磷污染主要來源于農(nóng)業(yè)的面源污染,見圖6。枯水期(1—2、11—12月)渭河流域陜西段的總流量較小,在污染物數(shù)量不變的情況下,污染物濃度要大一些,因此枯水期的水質(zhì)要次于平水期(3—5月)。豐水期(6—10月)的降水增大,導(dǎo)致地面徑流和地下徑流增加,從而使得更多的污染物進(jìn)入河道,污染物濃度因雨季稀釋的濃度小于因污染物進(jìn)入而增大的濃度,導(dǎo)致整體豐水期污染物濃度比平水期的污染物濃度大,比枯水期的污染物濃度小。根據(jù)圖1可知,污染較為嚴(yán)重的支流位于整個流域的中下游地區(qū),渭河干流兩岸的人口密集,導(dǎo)致了生活污水排放量較大,同時干流兩岸的國控、省控的重污染企業(yè)多,導(dǎo)致工業(yè)廢水排放量大。同時污染物隨著河流不斷累加,因此導(dǎo)致了上游水質(zhì)較好、中下游水質(zhì)較差的結(jié)果。

    渭河支流的新河、皂河、幸福渠、臨河和沋河等5個斷面出現(xiàn)Ⅳ、Ⅴ類水質(zhì)的概率較大,因此在表10中主要展示這5個斷面的不同時期的污染物貢獻(xiàn)比。

    a)豐水期

    b)平水期

    c)枯水期

    表10 不同斷面不同時期污染物貢獻(xiàn)比

    所用的評價指標(biāo)在計算權(quán)重的時候,只考慮了客觀權(quán)重計算方法,未考慮主觀因素,同時水質(zhì)的時空變化規(guī)律趨勢分析較為淺薄,在以后的研究中應(yīng)考慮主觀權(quán)重計算法,并在用綜合權(quán)重來評價水質(zhì)結(jié)果的同時也考慮水質(zhì)時空變化的原因。

    3 結(jié)論

    采用單因子評價法、熵權(quán)云模型綜合評價法、變權(quán)歐式距離法等,評價了陜西境內(nèi)渭河、漢江、丹江河段2020年1—6月的水質(zhì),主要結(jié)論如下。

    a)研究區(qū)域內(nèi)水質(zhì)狀況良好,70%以上的站點水質(zhì)達(dá)到Ⅲ類。污染較為嚴(yán)重的斷面主要出現(xiàn)在渭河支流,即新河、皂河、臨河、幸福渠和沋河等斷面。在這些斷面中,主要污染物指標(biāo)為化學(xué)需氧量和氨氮。

    b)3條主要河流的水質(zhì)評價結(jié)果中,熵權(quán)云模型綜合評價法評價結(jié)果低于單因子評價法,高于變權(quán)歐式距離法。尤其是小韋河入杏林監(jiān)測斷面,單因子評價法評價結(jié)果均為Ⅳ—Ⅴ類,而熵權(quán)云模型綜合評價法結(jié)果處于Ⅱ—Ⅲ類,變權(quán)歐式距離法的評價結(jié)果表明,2、5月處于Ⅳ—Ⅴ類,其余月份均為Ⅲ類。單因子評價法中用一個因子的污染指數(shù)判定水體的污染程度不夠嚴(yán)謹(jǐn),而熵權(quán)云模型綜合評價法與變權(quán)歐式距離法則削弱了各個水質(zhì)等級的界限對評價結(jié)果的影響,能較好地權(quán)衡表現(xiàn)較優(yōu)、較差2種污染因子對水體水質(zhì)的影響,尤其是對于熵權(quán)云模型綜合評價法而言,能更多地考慮不同污染因子的污染程度,評價結(jié)果偏向更為安全,因此,熵權(quán)云模型綜合評價法評價結(jié)果更為合理。

    猜你喜歡
    變權(quán)歐式水質(zhì)評價
    基于Creo軟件的石材歐式壁爐三維造型設(shè)計
    石材(2020年2期)2020-03-16 13:12:56
    一類特殊混合跳擴(kuò)散Black-Scholes模型的歐式回望期權(quán)定價
    歐式城堡——木炭與色彩的碰撞
    對我國小城鎮(zhèn)建設(shè)過程中歐式古典風(fēng)格建筑興起的思考
    變權(quán)空間權(quán)重構(gòu)造及空間效應(yīng)分析
    滻灞河水質(zhì)評價方法研究
    新疆生產(chǎn)建設(shè)兵團(tuán)城鎮(zhèn)化水平的變權(quán)組合預(yù)測
    基于概率統(tǒng)計和模糊綜合評價法的水質(zhì)評價模型及其應(yīng)用——以拉薩河水質(zhì)評價為例
    基于變權(quán)的物流資源公平分配方法
    基于SPAM的河流水質(zhì)評價模型
    99香蕉大伊视频| 亚洲精品国产色婷婷电影| 一本一本久久a久久精品综合妖精| 在线观看免费视频网站a站| 欧美 日韩 精品 国产| 女的被弄到高潮叫床怎么办| 欧美精品亚洲一区二区| 热re99久久精品国产66热6| 国产日韩一区二区三区精品不卡| 波多野结衣一区麻豆| 中文字幕人妻丝袜制服| 黑丝袜美女国产一区| 成人亚洲欧美一区二区av| 丝袜美腿诱惑在线| 久久天堂一区二区三区四区| 高清在线视频一区二区三区| 18禁观看日本| 丝袜在线中文字幕| 性色av一级| 九九爱精品视频在线观看| 精品亚洲乱码少妇综合久久| 一个人免费看片子| kizo精华| 九草在线视频观看| 免费观看人在逋| 一区在线观看完整版| 色综合欧美亚洲国产小说| 日韩成人av中文字幕在线观看| 天天躁日日躁夜夜躁夜夜| 国产亚洲一区二区精品| 免费少妇av软件| 国产成人精品福利久久| 又大又爽又粗| 99国产精品免费福利视频| av在线app专区| 秋霞伦理黄片| 啦啦啦 在线观看视频| 伊人亚洲综合成人网| 看非洲黑人一级黄片| 一区二区av电影网| 国产精品亚洲av一区麻豆 | 亚洲精品国产色婷婷电影| 亚洲婷婷狠狠爱综合网| 国产亚洲一区二区精品| 久久久国产欧美日韩av| 亚洲男人天堂网一区| 国产免费福利视频在线观看| 又大又黄又爽视频免费| 男女国产视频网站| 中文字幕av电影在线播放| 五月天丁香电影| videos熟女内射| 在线天堂中文资源库| 午夜日本视频在线| 亚洲综合色网址| 欧美日韩国产mv在线观看视频| 久久久久久免费高清国产稀缺| 国产免费又黄又爽又色| 亚洲精品美女久久av网站| 美女大奶头黄色视频| bbb黄色大片| 国产成人av激情在线播放| 午夜激情av网站| 在线亚洲精品国产二区图片欧美| 日韩视频在线欧美| 国产精品av久久久久免费| 高清黄色对白视频在线免费看| 日韩伦理黄色片| 啦啦啦在线观看免费高清www| 久久综合国产亚洲精品| 精品久久久精品久久久| 在线免费观看不下载黄p国产| 国产精品蜜桃在线观看| 美女午夜性视频免费| 成人免费观看视频高清| 欧美精品一区二区免费开放| 亚洲国产欧美在线一区| e午夜精品久久久久久久| 国产一区二区 视频在线| av不卡在线播放| 亚洲av综合色区一区| 成年动漫av网址| 日韩熟女老妇一区二区性免费视频| 成人午夜精彩视频在线观看| 亚洲成人免费av在线播放| 久久婷婷青草| 久久久久久久久久久免费av| 精品福利永久在线观看| 少妇被粗大猛烈的视频| 多毛熟女@视频| 咕卡用的链子| 美女视频免费永久观看网站| a级片在线免费高清观看视频| 国产女主播在线喷水免费视频网站| 人人妻人人澡人人看| 日韩大片免费观看网站| 久久久久视频综合| 亚洲 欧美一区二区三区| 9热在线视频观看99| 在线亚洲精品国产二区图片欧美| 国产又爽黄色视频| 国产伦理片在线播放av一区| 欧美日韩一区二区视频在线观看视频在线| 亚洲激情五月婷婷啪啪| 国产亚洲欧美精品永久| 国产精品熟女久久久久浪| 菩萨蛮人人尽说江南好唐韦庄| 丰满乱子伦码专区| 我要看黄色一级片免费的| 蜜桃国产av成人99| 欧美变态另类bdsm刘玥| 欧美变态另类bdsm刘玥| 免费观看人在逋| 精品免费久久久久久久清纯 | 飞空精品影院首页| 巨乳人妻的诱惑在线观看| 欧美激情极品国产一区二区三区| 国产免费福利视频在线观看| 七月丁香在线播放| 美女福利国产在线| 一级黄片播放器| 日本91视频免费播放| 久久精品国产亚洲av涩爱| 老鸭窝网址在线观看| 人人妻人人澡人人看| 国产麻豆69| 大码成人一级视频| 看免费av毛片| 午夜激情av网站| 国产精品香港三级国产av潘金莲 | 国产又色又爽无遮挡免| 精品免费久久久久久久清纯 | 欧美黄色片欧美黄色片| av网站免费在线观看视频| 免费不卡黄色视频| 亚洲精品av麻豆狂野| 一级a爱视频在线免费观看| 飞空精品影院首页| 日韩电影二区| 免费观看人在逋| 国产麻豆69| 亚洲 欧美一区二区三区| 亚洲三区欧美一区| av在线老鸭窝| 亚洲成人av在线免费| 又大又黄又爽视频免费| 黄片播放在线免费| 水蜜桃什么品种好| 美女午夜性视频免费| 国产成人精品久久久久久| 操美女的视频在线观看| 嫩草影院入口| 人人妻,人人澡人人爽秒播 | 久久久精品区二区三区| 一个人免费看片子| 久久久国产精品麻豆| 街头女战士在线观看网站| 欧美精品av麻豆av| 久久人人爽人人片av| 亚洲精品美女久久av网站| 亚洲免费av在线视频| 曰老女人黄片| 国产伦理片在线播放av一区| 日韩中文字幕视频在线看片| 黄网站色视频无遮挡免费观看| videosex国产| av一本久久久久| 亚洲国产成人一精品久久久| 国产福利在线免费观看视频| 看免费av毛片| 成人午夜精彩视频在线观看| 免费不卡黄色视频| 少妇人妻精品综合一区二区| 老司机影院毛片| 男女边吃奶边做爰视频| 免费黄网站久久成人精品| 男女免费视频国产| 国产成人免费无遮挡视频| 少妇精品久久久久久久| 制服诱惑二区| 欧美精品亚洲一区二区| 丝袜美腿诱惑在线| 黑人欧美特级aaaaaa片| 91成人精品电影| 纵有疾风起免费观看全集完整版| 午夜免费鲁丝| 青春草视频在线免费观看| 国产精品一国产av| 黑人欧美特级aaaaaa片| 国精品久久久久久国模美| 亚洲欧洲国产日韩| 亚洲精品av麻豆狂野| 国产无遮挡羞羞视频在线观看| 国产又爽黄色视频| 肉色欧美久久久久久久蜜桃| 一级片免费观看大全| 99热国产这里只有精品6| 亚洲国产精品一区三区| 最近2019中文字幕mv第一页| 少妇精品久久久久久久| 男女免费视频国产| 又粗又硬又长又爽又黄的视频| 女性被躁到高潮视频| 国产色婷婷99| 国产精品秋霞免费鲁丝片| 亚洲欧美精品综合一区二区三区| 天天添夜夜摸| 国产精品女同一区二区软件| 午夜福利乱码中文字幕| 一本—道久久a久久精品蜜桃钙片| 国产精品一二三区在线看| 香蕉丝袜av| 在线观看国产h片| 国产亚洲av片在线观看秒播厂| 午夜日本视频在线| 久久久久网色| 美女高潮到喷水免费观看| 中国三级夫妇交换| www.自偷自拍.com| 欧美精品一区二区免费开放| 最近2019中文字幕mv第一页| 国产成人精品久久久久久| 久久久久久久精品精品| 国产又色又爽无遮挡免| 一级片'在线观看视频| 久久久久久人人人人人| 亚洲国产欧美一区二区综合| 最近手机中文字幕大全| 大香蕉久久网| 青青草视频在线视频观看| tube8黄色片| 精品一品国产午夜福利视频| 制服丝袜香蕉在线| www.熟女人妻精品国产| 亚洲国产精品成人久久小说| 国产亚洲最大av| 777久久人妻少妇嫩草av网站| 国产色婷婷99| 黄频高清免费视频| 国产国语露脸激情在线看| 国产免费视频播放在线视频| 亚洲精品国产av蜜桃| 国产伦理片在线播放av一区| 国产精品一国产av| 久久韩国三级中文字幕| videos熟女内射| 满18在线观看网站| 黑人猛操日本美女一级片| 日韩一区二区三区影片| 狠狠婷婷综合久久久久久88av| 精品一区二区三区av网在线观看 | 99热全是精品| 一级爰片在线观看| av网站免费在线观看视频| 亚洲精品国产一区二区精华液| 亚洲欧美日韩另类电影网站| 欧美人与善性xxx| 国产熟女午夜一区二区三区| 精品久久久精品久久久| 男女床上黄色一级片免费看| 丰满乱子伦码专区| 巨乳人妻的诱惑在线观看| 99香蕉大伊视频| 国产精品99久久99久久久不卡 | 水蜜桃什么品种好| 男女下面插进去视频免费观看| 国产成人a∨麻豆精品| 黄片小视频在线播放| av网站在线播放免费| 亚洲伊人色综图| 亚洲熟女精品中文字幕| 国产乱来视频区| 亚洲第一区二区三区不卡| 久久国产亚洲av麻豆专区| 亚洲视频免费观看视频| 欧美国产精品一级二级三级| 美女中出高潮动态图| 国产一区二区在线观看av| 性色av一级| videosex国产| 国产激情久久老熟女| 极品少妇高潮喷水抽搐| 18在线观看网站| 久久99热这里只频精品6学生| 久久精品熟女亚洲av麻豆精品| 少妇人妻久久综合中文| 欧美成人精品欧美一级黄| 一级片免费观看大全| 亚洲自偷自拍图片 自拍| 欧美激情极品国产一区二区三区| 最近中文字幕高清免费大全6| 亚洲成人免费av在线播放| 日韩熟女老妇一区二区性免费视频| 国产精品久久久久久精品古装| 赤兔流量卡办理| 日本色播在线视频| 免费观看性生交大片5| 久久精品国产a三级三级三级| 国产精品av久久久久免费| xxx大片免费视频| 狂野欧美激情性bbbbbb| 精品人妻熟女毛片av久久网站| 18禁裸乳无遮挡动漫免费视频| 一边摸一边抽搐一进一出视频| 亚洲,欧美精品.| 久久久精品国产亚洲av高清涩受| 999久久久国产精品视频| 青草久久国产| 免费看av在线观看网站| 制服诱惑二区| 丝袜在线中文字幕| av免费观看日本| 国产日韩欧美视频二区| 免费高清在线观看视频在线观看| 国产日韩欧美视频二区| av在线观看视频网站免费| 国产精品一区二区精品视频观看| 欧美精品人与动牲交sv欧美| 老司机在亚洲福利影院| 亚洲av日韩在线播放| av电影中文网址| 久久久国产精品麻豆| 一区二区三区四区激情视频| 国产成人精品福利久久| 啦啦啦视频在线资源免费观看| 亚洲成人国产一区在线观看 | 观看av在线不卡| 国产成人a∨麻豆精品| 精品亚洲成a人片在线观看| 永久免费av网站大全| 国产视频首页在线观看| 国产爽快片一区二区三区| 狠狠婷婷综合久久久久久88av| 久久ye,这里只有精品| 亚洲av成人不卡在线观看播放网 | 一个人免费看片子| 香蕉国产在线看| 免费日韩欧美在线观看| 赤兔流量卡办理| 女人精品久久久久毛片| 国产精品一二三区在线看| 成年美女黄网站色视频大全免费| 久久久久久久久免费视频了| 久久久久久人妻| 永久免费av网站大全| 最近中文字幕高清免费大全6| 国产又爽黄色视频| 观看av在线不卡| 中文乱码字字幕精品一区二区三区| 在线精品无人区一区二区三| 亚洲成国产人片在线观看| 久久精品久久久久久久性| 免费观看av网站的网址| 高清欧美精品videossex| 免费高清在线观看日韩| 一级毛片我不卡| 99精国产麻豆久久婷婷| 亚洲欧美清纯卡通| 欧美xxⅹ黑人| 欧美乱码精品一区二区三区| 欧美激情 高清一区二区三区| 国产97色在线日韩免费| 免费看av在线观看网站| 国产女主播在线喷水免费视频网站| 男男h啪啪无遮挡| 亚洲男人天堂网一区| 精品午夜福利在线看| 午夜日韩欧美国产| 国产精品香港三级国产av潘金莲 | 99香蕉大伊视频| 一级a爱视频在线免费观看| 免费高清在线观看日韩| 制服人妻中文乱码| 免费少妇av软件| 操出白浆在线播放| 久久久久国产精品人妻一区二区| 纵有疾风起免费观看全集完整版| 一本久久精品| 捣出白浆h1v1| 伦理电影免费视频| 欧美另类一区| 中文字幕人妻丝袜制服| av一本久久久久| 午夜影院在线不卡| 极品人妻少妇av视频| 美女视频免费永久观看网站| 久久天躁狠狠躁夜夜2o2o | 久久国产亚洲av麻豆专区| 热99久久久久精品小说推荐| 久久久久精品久久久久真实原创| 777久久人妻少妇嫩草av网站| 国产免费现黄频在线看| 9色porny在线观看| 欧美av亚洲av综合av国产av | 亚洲精品国产区一区二| 精品国产露脸久久av麻豆| 国产精品欧美亚洲77777| 亚洲少妇的诱惑av| 在线观看免费日韩欧美大片| 精品国产国语对白av| 亚洲国产精品999| 可以免费在线观看a视频的电影网站 | 亚洲国产欧美日韩在线播放| 观看av在线不卡| 亚洲美女搞黄在线观看| 爱豆传媒免费全集在线观看| 一区二区日韩欧美中文字幕| 欧美日韩视频高清一区二区三区二| 精品人妻熟女毛片av久久网站| 一二三四中文在线观看免费高清| 亚洲成人av在线免费| 我要看黄色一级片免费的| 黑丝袜美女国产一区| 久久影院123| 精品一区二区三卡| 国产成人免费观看mmmm| www日本在线高清视频| 看十八女毛片水多多多| 99九九在线精品视频| 中文字幕另类日韩欧美亚洲嫩草| 午夜激情av网站| 欧美亚洲 丝袜 人妻 在线| 婷婷色av中文字幕| 精品一区二区三区四区五区乱码 | 蜜桃在线观看..| 丰满迷人的少妇在线观看| 91精品三级在线观看| 久久韩国三级中文字幕| 51午夜福利影视在线观看| 九色亚洲精品在线播放| 成人黄色视频免费在线看| 日韩制服丝袜自拍偷拍| 色吧在线观看| 亚洲精品久久成人aⅴ小说| 久久久久精品国产欧美久久久 | 女的被弄到高潮叫床怎么办| 色吧在线观看| 80岁老熟妇乱子伦牲交| 亚洲精品久久午夜乱码| 最近中文字幕2019免费版| 国产精品嫩草影院av在线观看| 精品福利永久在线观看| 亚洲av综合色区一区| 亚洲美女视频黄频| 美女中出高潮动态图| 中国国产av一级| 国产一区二区激情短视频 | 国精品久久久久久国模美| 亚洲成人av在线免费| 国产精品嫩草影院av在线观看| 亚洲精品视频女| 亚洲精品在线美女| 亚洲国产av影院在线观看| 国产免费又黄又爽又色| 久久久久人妻精品一区果冻| 久久久久久久久久久免费av| 亚洲精品日本国产第一区| 汤姆久久久久久久影院中文字幕| 午夜福利视频在线观看免费| 日韩大片免费观看网站| av视频免费观看在线观看| 午夜免费鲁丝| 国产不卡av网站在线观看| 一级毛片电影观看| 欧美变态另类bdsm刘玥| 国产亚洲午夜精品一区二区久久| 精品视频人人做人人爽| 国产国语露脸激情在线看| 久久久久久久大尺度免费视频| 中文欧美无线码| 免费看av在线观看网站| 国产精品无大码| 建设人人有责人人尽责人人享有的| 久久天躁狠狠躁夜夜2o2o | 激情视频va一区二区三区| 日韩av在线免费看完整版不卡| 男人爽女人下面视频在线观看| 亚洲av日韩在线播放| 中文字幕精品免费在线观看视频| 免费高清在线观看日韩| 色播在线永久视频| 大陆偷拍与自拍| av国产久精品久网站免费入址| 久久精品人人爽人人爽视色| 啦啦啦在线免费观看视频4| 国产在线视频一区二区| 久久精品久久精品一区二区三区| 国产伦理片在线播放av一区| 好男人视频免费观看在线| 人成视频在线观看免费观看| 如何舔出高潮| 在线观看三级黄色| 亚洲久久久国产精品| 美国免费a级毛片| 男女边摸边吃奶| 久久久久久久大尺度免费视频| 爱豆传媒免费全集在线观看| 亚洲男人天堂网一区| 丰满少妇做爰视频| 十八禁人妻一区二区| 亚洲图色成人| svipshipincom国产片| 看免费成人av毛片| 超碰成人久久| 国产 一区精品| 国产精品麻豆人妻色哟哟久久| 伦理电影免费视频| 亚洲人成77777在线视频| 国产精品99久久99久久久不卡 | 国产xxxxx性猛交| 久久热在线av| av.在线天堂| 亚洲三区欧美一区| 色吧在线观看| 午夜福利,免费看| 精品久久久精品久久久| 超色免费av| 极品少妇高潮喷水抽搐| 青春草亚洲视频在线观看| 黑丝袜美女国产一区| 侵犯人妻中文字幕一二三四区| 看非洲黑人一级黄片| 精品一区二区三卡| 啦啦啦在线免费观看视频4| 男女下面插进去视频免费观看| 亚洲国产欧美日韩在线播放| 午夜激情久久久久久久| 日韩,欧美,国产一区二区三区| 欧美人与性动交α欧美软件| 在线观看免费日韩欧美大片| 欧美av亚洲av综合av国产av | 丝袜美腿诱惑在线| 激情视频va一区二区三区| 成人国语在线视频| 中文字幕人妻丝袜一区二区 | 亚洲美女黄色视频免费看| 国产亚洲av片在线观看秒播厂| 成人18禁高潮啪啪吃奶动态图| 欧美日韩亚洲综合一区二区三区_| 麻豆精品久久久久久蜜桃| 国产精品秋霞免费鲁丝片| 国产精品久久久人人做人人爽| 久久99热这里只频精品6学生| 成人亚洲精品一区在线观看| 最近手机中文字幕大全| 亚洲一码二码三码区别大吗| 亚洲欧美精品综合一区二区三区| 久久99精品国语久久久| 性色av一级| 免费在线观看黄色视频的| 97在线人人人人妻| 成人三级做爰电影| h视频一区二区三区| 男女国产视频网站| 十八禁网站网址无遮挡| 日日爽夜夜爽网站| 女人久久www免费人成看片| 精品国产国语对白av| 男女下面插进去视频免费观看| 日韩中文字幕欧美一区二区 | 大陆偷拍与自拍| 老司机深夜福利视频在线观看 | 丝袜在线中文字幕| 精品久久久久久电影网| 天天影视国产精品| 日本猛色少妇xxxxx猛交久久| 宅男免费午夜| 国精品久久久久久国模美| 亚洲成人av在线免费| 成年人午夜在线观看视频| 亚洲av国产av综合av卡| 蜜桃国产av成人99| 性少妇av在线| 一区二区日韩欧美中文字幕| 精品国产国语对白av| 青草久久国产| 国产精品香港三级国产av潘金莲 | 女人高潮潮喷娇喘18禁视频| 人人澡人人妻人| 国产成人精品久久久久久| 久久久久久久国产电影| 免费av中文字幕在线| 叶爱在线成人免费视频播放| 国产一级毛片在线| 国产男女内射视频| 黑丝袜美女国产一区| 婷婷色av中文字幕| 国产亚洲午夜精品一区二区久久| 欧美日韩综合久久久久久| 激情视频va一区二区三区| 大片电影免费在线观看免费| 亚洲av日韩精品久久久久久密 | 亚洲av在线观看美女高潮| 国产一区二区 视频在线| 亚洲欧美一区二区三区黑人| 久久人人97超碰香蕉20202| 久久影院123| 国产福利在线免费观看视频| 天天躁夜夜躁狠狠躁躁| 久久久久精品国产欧美久久久 | 日韩电影二区| 99九九在线精品视频| 午夜av观看不卡| 国产av一区二区精品久久| 国语对白做爰xxxⅹ性视频网站| 波多野结衣一区麻豆| 欧美变态另类bdsm刘玥| 欧美人与性动交α欧美软件| 日韩精品有码人妻一区| 国产成人午夜福利电影在线观看| 99久国产av精品国产电影| 成人国产麻豆网| 高清视频免费观看一区二区| 777久久人妻少妇嫩草av网站| 最近最新中文字幕免费大全7|