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

    基于模糊Borda法的組合評價在林芝水質(zhì)評價中的應(yīng)用

    2020-12-21 01:55:32段小龍宗永臣黃德才郝凱越李遠(yuǎn)威
    水力發(fā)電 2020年9期
    關(guān)鍵詞:林芝水質(zhì)評價水源地

    段小龍,宗永臣,2,黃德才,郝凱越,李遠(yuǎn)威

    (1.西藏農(nóng)牧學(xué)院水利土木工程學(xué)院,西藏 林芝 860000;2.西藏林芝高山森林生態(tài)系統(tǒng)國家野外科學(xué)觀測研究站,西藏 林芝 860000)

    林芝地區(qū)位于西藏東南部,雅魯藏布江中下游,有“西藏江南”之美譽(yù),現(xiàn)為西藏經(jīng)濟(jì)社會最為發(fā)達(dá)的區(qū)域。由于特殊的地理與氣候原因,林芝市擁有豐富的水資源;同時,該地區(qū)生態(tài)環(huán)境脆弱,水源易受到污染[1-2];水資源是支撐當(dāng)?shù)卣麄€生態(tài)系統(tǒng)十分重要的物質(zhì)基礎(chǔ),一旦遭受污染,其經(jīng)濟(jì)社會發(fā)展必然受到制約。因此,亟需開展飲用水源水質(zhì)評價,為水源保護(hù)和水質(zhì)安全提供依據(jù)。水質(zhì)評價的實質(zhì)是選取一些水質(zhì)影響因子,運用數(shù)學(xué)模型對水質(zhì)現(xiàn)狀的優(yōu)劣進(jìn)行定性或定量地評價[3-5]。國內(nèi)外在水質(zhì)評價領(lǐng)域的研究成果較多。其中,具有代表性的方法有:模糊綜合評價法、主成分分析法、綜合污染指數(shù)、內(nèi)梅羅水質(zhì)指數(shù)等[6-11]。學(xué)者們對西藏地區(qū)水質(zhì)評價展開了富有成果的研究:Qu等通過分析西藏的11條河流水質(zhì)現(xiàn)狀,結(jié)果表明由于自然過程和采礦作業(yè)等人為活動,一些地區(qū)水體中的As、Cd等幾種重金屬元素含量高于世界衛(wèi)生組織的飲用水指南[12];Zhang等對藏東北的水質(zhì)進(jìn)行了調(diào)查,結(jié)果表明所有采樣點的水質(zhì)均為Ⅰ級,氨氮是主要污染物[13];Zong對林芝市農(nóng)村地區(qū)飲用水源質(zhì)量進(jìn)行的評估結(jié)果表明,9個村的水質(zhì)均滿足飲用水要求[14]。上述水質(zhì)研究成果對開展林芝市水源地水質(zhì)評價研究提供了參考。

    1 材料與方法

    1.1 研究區(qū)概況

    林芝市平均海拔3 100 m,處于北緯26°52~30°40′N,東經(jīng)92°09′~98°47′E之間,以有世界上最深的大峽谷著稱于世;并有世界上落差最大的垂直地貌分布,典型的高山峽谷地貌,有異常豐富的野生動植物資源。該地區(qū)屬于溫帶季風(fēng)氣候,降水充沛,大部分地區(qū)年平均降水量500~1 000 mm;陽光充足,年平均溫度8.7℃。全境水網(wǎng)密集,有冰川、冰湖及多條河流,城鎮(zhèn)供水基本來自山脈峽谷的河流。轄有林芝、米林、工布江達(dá)等7個縣55個鄉(xiāng)鎮(zhèn),約有23.1萬人口。

    1.2 水樣采集與測定

    本次研究從林芝市7個縣級以上水源地采集水樣(分別是S1市區(qū)、S2工布江達(dá)縣、S3米林縣、S4朗縣、S5波密縣、S6墨脫縣、S7察隅縣),監(jiān)測頻率為2019年雨季旱季各1次,分析數(shù)值為采樣點各水質(zhì)指標(biāo)實測平均值。根據(jù)林芝市往年水質(zhì)數(shù)據(jù)和研究需要,選取檢測項目有:化學(xué)需氧量、總磷、總氮、銅、鎘等11項指標(biāo)。水樣采集和保存依照有關(guān)國家標(biāo)準(zhǔn)的要求進(jìn)行[15-16]。常規(guī)水質(zhì)指標(biāo)現(xiàn)場使用便攜式多參數(shù)水質(zhì)分析儀檢測;其他指標(biāo)在實驗室按照《水和廢水監(jiān)測分析方法(第四版)》[17]中的方法進(jìn)行檢測分析。為保證數(shù)據(jù)可靠性,采集過程中設(shè)置空白樣和平行樣;實驗過程中用加標(biāo)回收等措施來控制實驗質(zhì)量,且數(shù)據(jù)的相對標(biāo)準(zhǔn)偏和加標(biāo)回收率均在合理范圍內(nèi)。

    1.3 水質(zhì)評價方法

    由于水體多種污染因子之間關(guān)系復(fù)雜,并有明顯的隨機(jī)性和模糊性,用單一的方法進(jìn)行評價很難解決各污染因子與水質(zhì)等級間的復(fù)雜關(guān)系[18]。鑒此,本文選取用基于主成分分析法和綜合污染指數(shù)法相組合的模糊Borda法對林芝市水源水質(zhì)進(jìn)行評價。

    1.3.1綜合污染指數(shù)法

    綜合污染指數(shù)法是對各種污染物的相對污染指數(shù)進(jìn)行歸納和統(tǒng)計得到一個代表水質(zhì)的綜合污染指數(shù)[19]。該方法能對水質(zhì)做定量分析,可以反映水體的污染程度和主要污染因子;同時,在進(jìn)行水質(zhì)綜合評價時,可以將各水源地的污染程度排序[20]。結(jié)合相關(guān)的研究,綜合污染指數(shù)與水質(zhì)污染程度對應(yīng)見表1[21]。本論文中采用該方法中算術(shù)平均法對水質(zhì)進(jìn)行綜合評價,其表達(dá)式為

    (1)

    pi=ci/si

    (2)

    式中,P為綜合污染指數(shù);Pi為第i種污染因子的單因子指數(shù);n為污染因子總數(shù);Ci為污染因子實測濃度;si為第i種污染因子的評價標(biāo)準(zhǔn)值,參見GB3838—2002《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》中的III類水標(biāo)準(zhǔn)值。

    1.3.2主成分分析法

    主成分分析法是一種多維因子納入同一個系統(tǒng)中進(jìn)行定量化研究的分析方法。其將較多影響因子轉(zhuǎn)化為獨立且關(guān)鍵的綜合因子,并利用方差較大的綜合因子反映原多個因子所包含的主要信息[22]。若前幾個主成分包含的信息量足夠分析所需時,可以忽略后面的主成分,從而使得水質(zhì)評價中各種數(shù)據(jù)的分析效率得以提高。主成分分析法基本分析步驟如下:①對原始數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化以消除量綱的影響;②利用標(biāo)準(zhǔn)化相關(guān)系數(shù)矩陣對林芝市水源地水質(zhì)影響因子進(jìn)行主成分分析,得到特征值、累計貢獻(xiàn)率;③進(jìn)一步求出主成分表達(dá)式,選取主成分;④根據(jù)新得出的主成分,綜合評價評價。

    1.3.3基于模糊 Borda法的組合評價方法

    模糊Barda法可以綜合多種評價方法的不同結(jié)果,該方法既考慮不同方法下排序名次的差異,又考慮相應(yīng)評價方法下各項目的得分值,能更好地利用已有的評價信息,從而使得評價結(jié)果具有較高的合理性和優(yōu)越性[23]。組合評價法的基本步驟如下。

    (1)用綜合污染指數(shù)法和主成分分析法分別對7個水源地進(jìn)行評價。

    (2)采用Spearman等級相關(guān)系數(shù)對上述評價的結(jié)果進(jìn)行一致性檢驗,具體計算為

    (3)

    式中,di為評價結(jié)果排序的等級差;n為指標(biāo)數(shù)量;r取值區(qū)間為[0,1],且當(dāng)取值越接近于1,其相關(guān)性越好。

    (3)運用基于模糊 Borda的組合評價方法進(jìn)行二次評價,組合計算過程詳見文獻(xiàn)[23]。

    (4)最后采用Kendall一致性系數(shù)對組合評價結(jié)果與單一評價結(jié)果的密切性進(jìn)行檢驗,具體計算如下

    W=s/[(1/12)k2(N3-N)]

    (4)

    (5)

    式中,W為一致性系數(shù);K為方法數(shù)量;N為采樣點數(shù)量;s為每個被評對象的離差平方和;Ri每個采樣點的等級之和。W取值區(qū)間為[0,1],取值越接近于1,其相關(guān)性越好。

    2 結(jié)果與分析

    2.1 綜合污染指數(shù)法評價結(jié)果

    根據(jù)各水源地實測的水質(zhì)指標(biāo)濃度,利用綜合污染指數(shù)法對所得到的數(shù)據(jù)進(jìn)行分析與處理,對林芝水源地水環(huán)境質(zhì)量現(xiàn)狀進(jìn)行評價,計算結(jié)果見表1。

    表1 綜合污染指數(shù)法分類標(biāo)準(zhǔn)及評價結(jié)果

    根據(jù)綜合污染指數(shù)法的評價結(jié)果可知,7個縣級以上水源地水質(zhì)的評價排序(由優(yōu)到劣):巴宜區(qū)、墨脫縣、波密縣朗縣、工布江達(dá)、察隅縣、米林縣。水源地水質(zhì)均在清潔和尚清潔范圍內(nèi),綜合污染指數(shù)在0.195~0.257之間。其中,巴宜區(qū)綜合污染指數(shù)最低(0.195),受污染程度最輕,可能與該地屬市級水源地,市政府對水源地保護(hù)區(qū)內(nèi)的居民和農(nóng)牧業(yè)等廢水排放管理相對更嚴(yán)格、故水體受到的外來污染以及人們干擾較少有關(guān);米林縣最高(0.257),受污染相對最嚴(yán)重,可能是米林縣水源地處在旅游區(qū)且水環(huán)境容量較小導(dǎo)致水體自凈能力較弱,加上生活污水和牧業(yè)廢水等點源和面源污染物不能及時被稀釋,造成水質(zhì)相對最差;巴宜區(qū)、墨脫縣、波密縣水源受到污染很少,多數(shù)項目未檢出,個別項目檢出但在標(biāo)準(zhǔn)內(nèi),水質(zhì)屬于清潔(P<0.2);朗縣、工布江達(dá)、察隅縣和米林縣水源地水質(zhì)屬于清潔受到污染較少,檢出值在標(biāo)準(zhǔn)內(nèi),個別項目接近超標(biāo),水質(zhì)尚屬清潔(P=0.2~0.4)??傮w來說,各水源地的綜合污染指數(shù)為均低于0.4,表明參評指標(biāo)的達(dá)標(biāo)情況較好,林芝天然飲用水水質(zhì)較為潔凈,可以作為飲用水水源使用。這結(jié)論這與靳會姣[24]和羅丹[25]的研究結(jié)論相吻合。

    2.2 主成分分析法評價結(jié)果

    利用主成分分析法對林芝7個縣級以上水源地進(jìn)行水質(zhì)評價,具體計算過程在SPSS軟件中進(jìn)行。為了進(jìn)一步明確使每個主成分的所代表意義,通過最大方差法結(jié)合主成分的特征值進(jìn)行因子旋轉(zhuǎn),旋轉(zhuǎn)后的因子載荷量,見表2。

    表2 旋轉(zhuǎn)后的因子載荷量

    從表2的各主成分的貢獻(xiàn)率可知,前4個主成分已經(jīng)能夠反映水質(zhì)狀況,第一主成分到第四主成分的累積貢獻(xiàn)率達(dá)到92.843%,第一主成分貢獻(xiàn)率32.638%攜帶的水質(zhì)信息最多,說明林芝市縣級以上飲用水源地水質(zhì)主要受第一主成分影響。根據(jù)計算出的載荷因子來看,氨氮、總磷、總氮等指標(biāo)與第一主成分相關(guān)系數(shù)絕對值都在0.75以上,鉛、鎘離子的載荷也達(dá)到了約0.82和0.65,表明林芝市飲用水源地氮、磷營養(yǎng)鹽污染和重金屬離子污染情況較為突出。這可能與大部分水源地保護(hù)區(qū)大多在有農(nóng)牧民生活、耕牧的山谷之中,且部分水源地保護(hù)區(qū)內(nèi)有民宅、工礦業(yè)有關(guān)。第二、三、四主成分中化學(xué)需氧量和硝態(tài)氮、銅和鋅離子等載荷較高,進(jìn)一步反映了水質(zhì)受到氮、磷營養(yǎng)鹽污染和重金屬污染情況明顯。

    根據(jù)上表計算出各水源地4個主成分得分F1~F4(見圖1),并根據(jù)4個主成分的貢獻(xiàn)率為權(quán)重最終得出綜合主成分得分F(見圖2)。如圖1所示,各主成分得分在空間上呈明顯的變化,察隅和米林水源地的F1與F2的值均較高,說明這兩個水源地的水質(zhì)污染物濃度均較高;其原因可能是農(nóng)牧業(yè)面源污染物滲流有有關(guān)。墨脫縣、波密縣和巴宜區(qū)的F1、F2、F3的值都較低,說明這3各水源地水質(zhì)污染物濃度都較低,表明水質(zhì)受污染程度較小。綜合得分越高,說明水樣中污染物濃度越高,水質(zhì)越差。根據(jù)主成分綜合得分可以得出,7個縣級以上水源地水質(zhì)的評價排序(由優(yōu)到劣):墨脫縣、波密縣、巴宜區(qū)、工布江達(dá)縣、朗縣、察隅縣、米林縣。各水源地得分在-1.32~1.89之間。從綜合得分的數(shù)值特征來看,墨脫縣、巴宜區(qū)、波密縣和工布江達(dá)縣得分值均為負(fù)數(shù),水質(zhì)狀況較好,然后依次是朗縣、察隅縣和米林縣得分值均為正數(shù),水質(zhì)狀況相對較差。其中,墨脫水質(zhì)最好,得分為-1.32;而米林縣水質(zhì)狀況相對是相對最差,得分為1.89,二者相差3.21。上述結(jié)論與綜合污染指數(shù)法的結(jié)論基本一致。

    圖1 主成分得分在水源地的分布

    圖2 林芝市各水源地水質(zhì)主成分法綜合評價結(jié)果

    2.3 基于模糊 Borda法的組合評價結(jié)果

    兩種評價模型的得分值和排序顯示,評價結(jié)果具有一定的差異,各水源地水質(zhì)通過綜合污染指數(shù)法所得分值較為接近,主成分分析法的得分值變化較大;且巴宜區(qū)在兩種評價方法中結(jié)果有明顯不同。所以,需要用Spearman等級相關(guān)系數(shù)進(jìn)行事前檢驗,經(jīng)過計算相關(guān)系數(shù)為0.857。由此說明,這兩種評價方法具有高度正相關(guān),可以將兩種方法進(jìn)行組合,形成組合評價法。

    上面兩種單一評價模型通過了一致性檢驗,但是水源地水質(zhì)排序有一定的差異。為了消除這種差異性同時克服單一評價方法的不足,將兩種評價方法用模糊Borda 法形成一種新的組合評價法,通過上文的分析步驟得到評價得分和排序,為了便于比較分析將三種方法的評價結(jié)果及排序整合在一起,見表4。從組合評價結(jié)果可以看出,墨脫縣、巴宜區(qū)和波密縣的水源地在組合排名結(jié)果靠前,朗縣、工布江達(dá)縣和察隅縣排序較為靠后,而米林縣排名在最后。組合的模糊Borda法結(jié)果排序與兩種單一評價模型結(jié)果排序差別較小(排序上下浮動一名),但是只從數(shù)據(jù)上無法判斷組合評價結(jié)果是否與原始評價結(jié)果密切相關(guān),故對組合評價模型的評價結(jié)果進(jìn)行檢驗。

    表4 單一評價模型與組合模型的得分及排序

    在實際應(yīng)用中,采用組合評價方法需要滿足原始方法的評價結(jié)果和組合方法評價的結(jié)果密切相關(guān)。為滿足上述條件需要采用Kendall一致性系數(shù)進(jìn)行事后檢驗。計算得出一致性系數(shù)為0.94,由此認(rèn)為兩種方法與基于模糊Borda法組合評價法的結(jié)果具有很好的一致性。即,組合評價法有效性較高。通過單一評價模型與組合模型的得分及排序的評價結(jié)果排序?qū)φ漳軌蚩闯觯?種方法的結(jié)果排序基本一致。這說明基于模糊Borda法的組合評價模型在水質(zhì)評價中是可行的。結(jié)合水源地周邊的工礦業(yè)、農(nóng)牧業(yè)及居民區(qū)等污染源的分布看,組合評價結(jié)果與水源地自然環(huán)境現(xiàn)狀相符。該方法汲取了兩種單一評價方法各自的優(yōu)勢,同時考慮了指標(biāo)間的相互關(guān)系和數(shù)據(jù)的離散特征,使得評價結(jié)果更為全面,具有明顯的合理性。

    3 結(jié) 論

    本文對林芝市7個縣級以上水源地水體的化學(xué)需氧量、氨氮、鎘、鉛等11項指標(biāo)為對象,應(yīng)用基于主成分分析法和綜合污染指數(shù)法的模糊Borda組合評價法對其進(jìn)行了水質(zhì)評價,并通過了事前與事后檢驗,結(jié)果表明:林芝市7個縣級以上水源地全年的綜合污染指數(shù)在0.195~0.257之間,綜合污染指數(shù)得分表明水質(zhì)較好水源地有巴宜區(qū)和墨脫,水質(zhì)較差的為察隅和米林縣,各水源地的各項水質(zhì)檢測值均符合國家地表水水質(zhì)標(biāo)準(zhǔn),總體上林芝水源地水質(zhì)較為潔凈。由主成分分析法篩選出4個因子來看,重金屬離子和氮磷營養(yǎng)鹽在林芝水源地中起主導(dǎo)作用,對水質(zhì)影響最大,如需要應(yīng)加大監(jiān)測頻率和采取相應(yīng)治污措施,保障人們的飲水安全。模糊Borda組合評價法結(jié)果與兩種單一評價法結(jié)果排序基本一致,且評價結(jié)果更為全面,故建議采用多種評價方法相組合的方法來開展水質(zhì)評價工作。

    猜你喜歡
    林芝水質(zhì)評價水源地
    林芝市雷電風(fēng)險區(qū)劃
    林芝桃花香
    滻灞河水質(zhì)評價方法研究
    基于概率統(tǒng)計和模糊綜合評價法的水質(zhì)評價模型及其應(yīng)用——以拉薩河水質(zhì)評價為例
    基于因果分析圖法的飲用水源地保護(hù)探討
    基于SPAM的河流水質(zhì)評價模型
    A new species of the genus Amolops (Amphibia: Ranidae) from southeastern Tibet, China
    南水北調(diào)中線水源地土壤侵蝕經(jīng)濟(jì)損失估算
    吉林省供水水源地問題探討
    礦山疏干排水對地下水源地影響的數(shù)值模擬
    最近中文字幕2019免费版| 18禁黄网站禁片午夜丰满| 香蕉国产在线看| 各种免费的搞黄视频| 亚洲欧洲日产国产| 国产成人免费无遮挡视频| 高清不卡的av网站| 97精品久久久久久久久久精品| 久久久久久久大尺度免费视频| 天天影视国产精品| 岛国毛片在线播放| 这个男人来自地球电影免费观看| 99热网站在线观看| 亚洲欧美日韩高清在线视频 | 男女之事视频高清在线观看 | 啦啦啦视频在线资源免费观看| 欧美 亚洲 国产 日韩一| 欧美国产精品va在线观看不卡| 久热爱精品视频在线9| 黑人巨大精品欧美一区二区蜜桃| 国产不卡av网站在线观看| 大陆偷拍与自拍| 国产精品国产三级专区第一集| 一区二区av电影网| 久久久久久久久免费视频了| 国产精品免费视频内射| 晚上一个人看的免费电影| 亚洲三区欧美一区| 免费日韩欧美在线观看| 99国产精品99久久久久| 午夜免费鲁丝| 两个人看的免费小视频| 久久午夜综合久久蜜桃| 精品一区二区三区四区五区乱码 | 各种免费的搞黄视频| 国产一区二区三区av在线| 青草久久国产| 欧美日韩视频精品一区| 国产成人啪精品午夜网站| 中文字幕最新亚洲高清| 曰老女人黄片| 国产精品久久久久久精品古装| 国产在线观看jvid| 国产精品.久久久| 欧美激情高清一区二区三区| 别揉我奶头~嗯~啊~动态视频 | 欧美日韩视频精品一区| 欧美大码av| 不卡av一区二区三区| 亚洲国产欧美日韩在线播放| 亚洲精品一二三| 免费女性裸体啪啪无遮挡网站| 一边摸一边做爽爽视频免费| 丝瓜视频免费看黄片| 国产主播在线观看一区二区 | 精品一区在线观看国产| 国产成人免费无遮挡视频| 亚洲五月色婷婷综合| 久久精品国产a三级三级三级| 免费在线观看完整版高清| 99国产综合亚洲精品| 国产色视频综合| 黑人猛操日本美女一级片| 丝袜美足系列| 久久人人97超碰香蕉20202| 国产精品久久久人人做人人爽| videosex国产| 欧美日韩福利视频一区二区| 精品少妇久久久久久888优播| 国产精品久久久久久精品电影小说| 精品亚洲乱码少妇综合久久| 成人亚洲欧美一区二区av| 青草久久国产| 国产成人啪精品午夜网站| 青青草视频在线视频观看| 99九九在线精品视频| 高清欧美精品videossex| 亚洲精品国产av成人精品| 国产精品久久久久久精品电影小说| 日韩 亚洲 欧美在线| 咕卡用的链子| 亚洲国产欧美网| 女人爽到高潮嗷嗷叫在线视频| 日日爽夜夜爽网站| 中文字幕av电影在线播放| 大香蕉久久成人网| 汤姆久久久久久久影院中文字幕| 麻豆av在线久日| 这个男人来自地球电影免费观看| 高清黄色对白视频在线免费看| 国产成人精品无人区| 1024视频免费在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 黄色片一级片一级黄色片| 人人妻人人澡人人看| 欧美另类一区| 高清欧美精品videossex| 另类精品久久| 国产91精品成人一区二区三区 | 欧美在线黄色| 国产亚洲精品久久久久5区| 国产精品一国产av| www.熟女人妻精品国产| 亚洲国产最新在线播放| 嫩草影视91久久| av视频免费观看在线观看| 亚洲七黄色美女视频| 中文欧美无线码| 美女扒开内裤让男人捅视频| 国产精品久久久久成人av| 伊人久久大香线蕉亚洲五| 尾随美女入室| 制服诱惑二区| 黄片播放在线免费| 亚洲欧美激情在线| 欧美变态另类bdsm刘玥| 夫妻性生交免费视频一级片| 国产成人一区二区三区免费视频网站 | 女人精品久久久久毛片| 中文字幕av电影在线播放| 999精品在线视频| 涩涩av久久男人的天堂| 深夜精品福利| 成年人黄色毛片网站| 制服人妻中文乱码| 丝袜美足系列| 欧美av亚洲av综合av国产av| 大香蕉久久网| bbb黄色大片| 久久人人爽av亚洲精品天堂| 国产亚洲av高清不卡| 日韩电影二区| 国产成人精品久久二区二区91| 国产精品一国产av| 成人18禁高潮啪啪吃奶动态图| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲熟女毛片儿| 亚洲精品久久成人aⅴ小说| 狠狠婷婷综合久久久久久88av| 视频区图区小说| 看免费av毛片| av线在线观看网站| 亚洲av电影在线进入| 国产欧美日韩综合在线一区二区| 丰满迷人的少妇在线观看| 久久人人爽人人片av| 日韩欧美一区视频在线观看| 欧美亚洲 丝袜 人妻 在线| www.熟女人妻精品国产| 精品福利观看| 两个人免费观看高清视频| 满18在线观看网站| 午夜久久久在线观看| 久久精品亚洲熟妇少妇任你| 精品亚洲成a人片在线观看| 午夜福利在线免费观看网站| 午夜老司机福利片| av国产久精品久网站免费入址| 日韩制服丝袜自拍偷拍| 人人妻人人澡人人爽人人夜夜| 午夜av观看不卡| 久久午夜综合久久蜜桃| 亚洲国产欧美日韩在线播放| 国产高清不卡午夜福利| 日本vs欧美在线观看视频| 亚洲一区中文字幕在线| 在线精品无人区一区二区三| 欧美日韩亚洲国产一区二区在线观看 | 亚洲欧美色中文字幕在线| 麻豆av在线久日| 中文字幕av电影在线播放| 搡老岳熟女国产| 18在线观看网站| 黄色毛片三级朝国网站| 日本91视频免费播放| 精品免费久久久久久久清纯 | 五月天丁香电影| 国产日韩欧美在线精品| 永久免费av网站大全| 国产极品粉嫩免费观看在线| 精品高清国产在线一区| 免费少妇av软件| 久久99一区二区三区| 啦啦啦在线观看免费高清www| 9热在线视频观看99| 国产精品香港三级国产av潘金莲 | 在现免费观看毛片| 久久毛片免费看一区二区三区| 久久人妻福利社区极品人妻图片 | 精品亚洲成国产av| 熟女av电影| 精品国产乱码久久久久久男人| 色精品久久人妻99蜜桃| 日韩制服骚丝袜av| 国产高清国产精品国产三级| 久久精品国产综合久久久| 亚洲欧美日韩另类电影网站| 亚洲成人免费av在线播放| 另类亚洲欧美激情| 国产精品.久久久| 久久精品成人免费网站| 精品国产一区二区久久| 午夜视频精品福利| 女人高潮潮喷娇喘18禁视频| 日本91视频免费播放| 国产精品.久久久| 亚洲国产欧美一区二区综合| 日韩大片免费观看网站| 久久人人爽人人片av| 99国产精品免费福利视频| 欧美日韩一级在线毛片| 久久精品亚洲熟妇少妇任你| 深夜精品福利| 后天国语完整版免费观看| 校园人妻丝袜中文字幕| 亚洲少妇的诱惑av| 久久免费观看电影| 亚洲中文日韩欧美视频| 欧美在线黄色| 麻豆国产av国片精品| 日日摸夜夜添夜夜爱| 丁香六月天网| 亚洲av日韩在线播放| 国产片内射在线| 啦啦啦啦在线视频资源| 亚洲成av片中文字幕在线观看| 国产国语露脸激情在线看| 欧美少妇被猛烈插入视频| 久久久久久免费高清国产稀缺| 99国产精品免费福利视频| 国产麻豆69| 婷婷色综合www| 久热这里只有精品99| 欧美精品一区二区免费开放| 性高湖久久久久久久久免费观看| 各种免费的搞黄视频| 亚洲欧美精品综合一区二区三区| 在线天堂中文资源库| bbb黄色大片| 一本久久精品| 亚洲成人国产一区在线观看 | 午夜日韩欧美国产| 在线观看免费视频网站a站| 婷婷丁香在线五月| 日韩免费高清中文字幕av| 久久亚洲精品不卡| 99精品久久久久人妻精品| 99久久99久久久精品蜜桃| 波多野结衣av一区二区av| 国产福利在线免费观看视频| 亚洲国产精品一区三区| 纯流量卡能插随身wifi吗| 免费高清在线观看日韩| 久久久国产一区二区| 一区二区三区激情视频| 成在线人永久免费视频| 久久久久视频综合| 免费一级毛片在线播放高清视频 | 伦理电影免费视频| 男人爽女人下面视频在线观看| 国产野战对白在线观看| 国产精品国产三级专区第一集| 真人做人爱边吃奶动态| 久久毛片免费看一区二区三区| 麻豆国产av国片精品| 欧美性长视频在线观看| 黄色一级大片看看| 中文精品一卡2卡3卡4更新| 麻豆乱淫一区二区| 日韩熟女老妇一区二区性免费视频| 中文字幕色久视频| 中文字幕亚洲精品专区| 亚洲人成77777在线视频| 欧美日韩亚洲国产一区二区在线观看 | 夫妻午夜视频| 欧美+亚洲+日韩+国产| 性色av乱码一区二区三区2| 国产人伦9x9x在线观看| 日本91视频免费播放| 这个男人来自地球电影免费观看| 亚洲欧美清纯卡通| 波多野结衣av一区二区av| 国产精品成人在线| 久久久久网色| www.av在线官网国产| 色综合欧美亚洲国产小说| 国产一级毛片在线| 纵有疾风起免费观看全集完整版| 午夜免费鲁丝| 天天影视国产精品| 欧美日韩一级在线毛片| 国产黄色视频一区二区在线观看| 宅男免费午夜| 人人妻人人澡人人看| 午夜激情av网站| 在线天堂中文资源库| 色综合欧美亚洲国产小说| 另类精品久久| 交换朋友夫妻互换小说| 高清不卡的av网站| 少妇被粗大的猛进出69影院| 久久久精品国产亚洲av高清涩受| 国产亚洲av高清不卡| 岛国毛片在线播放| 又粗又硬又长又爽又黄的视频| 久久久精品区二区三区| 国产有黄有色有爽视频| 亚洲,一卡二卡三卡| 99久久人妻综合| 99香蕉大伊视频| 日本欧美国产在线视频| 日本黄色日本黄色录像| 国产精品一区二区在线不卡| 免费看不卡的av| 少妇的丰满在线观看| 亚洲伊人色综图| 久9热在线精品视频| 午夜两性在线视频| 国产精品偷伦视频观看了| 97人妻天天添夜夜摸| 肉色欧美久久久久久久蜜桃| 亚洲国产最新在线播放| 久久中文字幕一级| xxxhd国产人妻xxx| 国产一区二区三区综合在线观看| 精品少妇久久久久久888优播| 欧美激情 高清一区二区三区| 黑人猛操日本美女一级片| 久久狼人影院| 爱豆传媒免费全集在线观看| 欧美av亚洲av综合av国产av| 91麻豆av在线| 男人添女人高潮全过程视频| 亚洲伊人久久精品综合| 亚洲伊人色综图| 人人妻人人澡人人看| 亚洲国产成人一精品久久久| 欧美日韩福利视频一区二区| 伊人久久大香线蕉亚洲五| 最黄视频免费看| 欧美亚洲日本最大视频资源| 99香蕉大伊视频| 免费在线观看影片大全网站 | 91精品国产国语对白视频| 欧美老熟妇乱子伦牲交| 中国国产av一级| 亚洲视频免费观看视频| 制服人妻中文乱码| 亚洲,欧美,日韩| 久久国产精品人妻蜜桃| 99国产精品免费福利视频| 久久综合国产亚洲精品| 亚洲国产日韩一区二区| 一区二区三区四区激情视频| 看十八女毛片水多多多| 美女福利国产在线| 国产av精品麻豆| 精品卡一卡二卡四卡免费| videosex国产| 国语对白做爰xxxⅹ性视频网站| 成人亚洲欧美一区二区av| 9191精品国产免费久久| 少妇裸体淫交视频免费看高清 | 亚洲精品第二区| 中国美女看黄片| 一级a爱视频在线免费观看| 欧美精品人与动牲交sv欧美| 久久人人爽人人片av| 国产黄色免费在线视频| 亚洲欧美日韩另类电影网站| 中文欧美无线码| 妹子高潮喷水视频| 99热国产这里只有精品6| 女人被躁到高潮嗷嗷叫费观| 亚洲人成电影免费在线| 十八禁高潮呻吟视频| 久久精品亚洲av国产电影网| 老汉色∧v一级毛片| 后天国语完整版免费观看| 欧美亚洲日本最大视频资源| 视频区图区小说| 老熟女久久久| 国产一区亚洲一区在线观看| 男女无遮挡免费网站观看| 国产男人的电影天堂91| h视频一区二区三区| 欧美国产精品一级二级三级| 久久久精品免费免费高清| 激情五月婷婷亚洲| 亚洲第一青青草原| 国产99久久九九免费精品| 国产欧美日韩精品亚洲av| 国精品久久久久久国模美| 性高湖久久久久久久久免费观看| 大片免费播放器 马上看| 色网站视频免费| 一本大道久久a久久精品| 亚洲七黄色美女视频| 亚洲一区中文字幕在线| 又大又爽又粗| 欧美亚洲日本最大视频资源| 少妇人妻 视频| 精品人妻一区二区三区麻豆| av又黄又爽大尺度在线免费看| 美女中出高潮动态图| 高清av免费在线| 一区福利在线观看| 欧美久久黑人一区二区| 亚洲精品美女久久av网站| 少妇猛男粗大的猛烈进出视频| 色94色欧美一区二区| 人妻人人澡人人爽人人| 欧美在线一区亚洲| 一区二区三区激情视频| 女警被强在线播放| 国产成人系列免费观看| svipshipincom国产片| 亚洲成人国产一区在线观看 | 中文字幕人妻丝袜制服| 男人爽女人下面视频在线观看| 一级,二级,三级黄色视频| 建设人人有责人人尽责人人享有的| 久久精品亚洲av国产电影网| 91字幕亚洲| 99热网站在线观看| 久久狼人影院| 欧美久久黑人一区二区| 国产免费现黄频在线看| 老司机亚洲免费影院| 久久免费观看电影| 国产成人欧美| 国产成人免费观看mmmm| 日韩 欧美 亚洲 中文字幕| svipshipincom国产片| 免费在线观看完整版高清| 国精品久久久久久国模美| 天天添夜夜摸| 午夜激情av网站| 中国美女看黄片| 国产一级毛片在线| 最新在线观看一区二区三区 | 免费不卡黄色视频| 色精品久久人妻99蜜桃| 欧美日韩av久久| 免费人妻精品一区二区三区视频| 飞空精品影院首页| 国产精品亚洲av一区麻豆| 美女午夜性视频免费| 国产精品久久久久成人av| 欧美日韩亚洲综合一区二区三区_| 一本大道久久a久久精品| 久久久久久久久久久久大奶| 在线观看www视频免费| 精品国产乱码久久久久久小说| 女人高潮潮喷娇喘18禁视频| √禁漫天堂资源中文www| 欧美激情 高清一区二区三区| 久久精品久久精品一区二区三区| 视频区欧美日本亚洲| 免费在线观看视频国产中文字幕亚洲 | 最新的欧美精品一区二区| 超碰97精品在线观看| 久久精品久久久久久噜噜老黄| 日本91视频免费播放| 国产精品 欧美亚洲| 尾随美女入室| 校园人妻丝袜中文字幕| 香蕉丝袜av| 9热在线视频观看99| 久久九九热精品免费| 老司机靠b影院| 人人妻人人澡人人看| 在线观看一区二区三区激情| 欧美大码av| 视频区欧美日本亚洲| 在线 av 中文字幕| 两性夫妻黄色片| www.自偷自拍.com| 2018国产大陆天天弄谢| 人人妻人人爽人人添夜夜欢视频| 丰满少妇做爰视频| 中国美女看黄片| 午夜精品国产一区二区电影| 一本色道久久久久久精品综合| 在线观看免费视频网站a站| 色视频在线一区二区三区| 99精国产麻豆久久婷婷| 日韩 亚洲 欧美在线| 别揉我奶头~嗯~啊~动态视频 | 永久免费av网站大全| 久久久久久免费高清国产稀缺| 黄色片一级片一级黄色片| 精品一区二区三卡| 看免费av毛片| 国产精品 国内视频| 91精品伊人久久大香线蕉| 男女免费视频国产| a 毛片基地| 麻豆av在线久日| 精品福利观看| av网站免费在线观看视频| 成年人午夜在线观看视频| 亚洲一码二码三码区别大吗| 中文字幕精品免费在线观看视频| 亚洲国产日韩一区二区| 精品高清国产在线一区| 国产熟女欧美一区二区| 超碰成人久久| 51午夜福利影视在线观看| 色播在线永久视频| 黄色一级大片看看| 亚洲,一卡二卡三卡| 久久鲁丝午夜福利片| 黄色怎么调成土黄色| 亚洲精品美女久久久久99蜜臀 | 少妇被粗大的猛进出69影院| 亚洲欧洲日产国产| 少妇被粗大的猛进出69影院| 精品一区在线观看国产| 熟女av电影| 精品少妇久久久久久888优播| 亚洲精品美女久久久久99蜜臀 | 亚洲国产欧美网| 十八禁高潮呻吟视频| 亚洲国产av影院在线观看| 午夜视频精品福利| 精品国产乱码久久久久久小说| 亚洲专区中文字幕在线| 日韩制服骚丝袜av| 搡老岳熟女国产| 一级毛片电影观看| 尾随美女入室| 久久综合国产亚洲精品| 亚洲精品美女久久久久99蜜臀 | 久久天堂一区二区三区四区| 激情五月婷婷亚洲| 一本一本久久a久久精品综合妖精| 在线观看国产h片| 一本一本久久a久久精品综合妖精| 黑人欧美特级aaaaaa片| 免费在线观看完整版高清| 免费在线观看日本一区| 一本一本久久a久久精品综合妖精| xxxhd国产人妻xxx| videosex国产| 久久精品亚洲熟妇少妇任你| cao死你这个sao货| 老司机在亚洲福利影院| 亚洲av欧美aⅴ国产| 亚洲精品第二区| 91国产中文字幕| 欧美日韩亚洲高清精品| 日韩制服骚丝袜av| 国产精品麻豆人妻色哟哟久久| 国产亚洲欧美精品永久| 亚洲人成网站在线观看播放| 久久综合国产亚洲精品| 黄色视频不卡| 久久精品国产综合久久久| 肉色欧美久久久久久久蜜桃| 男人爽女人下面视频在线观看| 亚洲激情五月婷婷啪啪| 超碰97精品在线观看| bbb黄色大片| 波多野结衣av一区二区av| 午夜久久久在线观看| 亚洲,欧美精品.| 久久 成人 亚洲| 男女之事视频高清在线观看 | a级毛片在线看网站| av线在线观看网站| 9色porny在线观看| 你懂的网址亚洲精品在线观看| 999精品在线视频| 亚洲国产中文字幕在线视频| 免费高清在线观看日韩| 亚洲精品在线美女| 亚洲精品国产一区二区精华液| 人人妻人人澡人人爽人人夜夜| 亚洲三区欧美一区| 中文精品一卡2卡3卡4更新| 美女高潮到喷水免费观看| 一级黄片播放器| 久久久久视频综合| 国产精品九九99| 欧美成人午夜精品| 18禁裸乳无遮挡动漫免费视频| 成年人免费黄色播放视频| 国产精品免费视频内射| 国产精品亚洲av一区麻豆| 久久鲁丝午夜福利片| 久久久久国产精品人妻一区二区| 一本—道久久a久久精品蜜桃钙片| 91精品三级在线观看| 水蜜桃什么品种好| 最新在线观看一区二区三区 | 国产日韩一区二区三区精品不卡| 久久人妻福利社区极品人妻图片 | 成人三级做爰电影| 女性生殖器流出的白浆| 一区二区三区激情视频| 熟女少妇亚洲综合色aaa.| 国产高清不卡午夜福利| 美女视频免费永久观看网站| 一级毛片电影观看| 男女之事视频高清在线观看 | 天堂中文最新版在线下载| 99热全是精品| 亚洲天堂av无毛| 午夜免费成人在线视频| 亚洲国产欧美在线一区| 国产午夜精品一二区理论片| 日日爽夜夜爽网站| 最近最新中文字幕大全免费视频 | 视频在线观看一区二区三区|