• <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ù)值模擬
    女人久久www免费人成看片| 欧美亚洲 丝袜 人妻 在线| 欧美日韩精品成人综合77777| 色吧在线观看| 一区二区三区免费毛片| 超碰97精品在线观看| 国产亚洲av片在线观看秒播厂| 亚洲无线观看免费| 日韩欧美 国产精品| 国产伦精品一区二区三区视频9| 日本猛色少妇xxxxx猛交久久| 免费大片18禁| 午夜激情福利司机影院| 一级毛片aaaaaa免费看小| 网址你懂的国产日韩在线| 免费看a级黄色片| 最近最新中文字幕免费大全7| 校园人妻丝袜中文字幕| 成年免费大片在线观看| 国产精品国产av在线观看| 插逼视频在线观看| 亚洲国产av新网站| 亚洲aⅴ乱码一区二区在线播放| 禁无遮挡网站| 老师上课跳d突然被开到最大视频| 不卡视频在线观看欧美| 在线a可以看的网站| 日本色播在线视频| 校园人妻丝袜中文字幕| 麻豆精品久久久久久蜜桃| 欧美一级a爱片免费观看看| 成人免费观看视频高清| 婷婷色综合大香蕉| 观看美女的网站| 国产黄色视频一区二区在线观看| 国产高清三级在线| 国产有黄有色有爽视频| 我的女老师完整版在线观看| 丰满少妇做爰视频| 亚洲精品一区蜜桃| 99热国产这里只有精品6| 久久精品国产自在天天线| 亚洲精品国产色婷婷电影| 嫩草影院入口| 国产爽快片一区二区三区| 日韩成人伦理影院| 午夜日本视频在线| 我要看日韩黄色一级片| 三级国产精品片| 日韩一区二区视频免费看| 在线a可以看的网站| 黄片wwwwww| 欧美性猛交╳xxx乱大交人| 欧美xxxx性猛交bbbb| 在线观看一区二区三区激情| 91久久精品电影网| 久久人人爽人人爽人人片va| 精品一区二区三区视频在线| 最近最新中文字幕大全电影3| 欧美日韩精品成人综合77777| 青春草亚洲视频在线观看| 好男人在线观看高清免费视频| 观看免费一级毛片| 中文天堂在线官网| 亚洲最大成人av| 亚洲欧美日韩另类电影网站 | 亚洲天堂国产精品一区在线| av.在线天堂| 菩萨蛮人人尽说江南好唐韦庄| 日日摸夜夜添夜夜爱| 国产精品久久久久久久久免| 大香蕉久久网| 精品一区二区免费观看| 日本一二三区视频观看| 国产探花在线观看一区二区| 久久热精品热| 99九九线精品视频在线观看视频| 成年免费大片在线观看| 日本熟妇午夜| 欧美bdsm另类| 久久韩国三级中文字幕| 亚洲人成网站在线播| 在线精品无人区一区二区三 | 欧美性猛交╳xxx乱大交人| 狂野欧美白嫩少妇大欣赏| 成人欧美大片| 日韩av在线免费看完整版不卡| 七月丁香在线播放| 春色校园在线视频观看| 国产亚洲av片在线观看秒播厂| 国产精品一区二区在线观看99| 97人妻精品一区二区三区麻豆| 亚洲av一区综合| 男人狂女人下面高潮的视频| 麻豆成人av视频| 国产亚洲5aaaaa淫片| 亚洲精品国产成人久久av| 国产黄色视频一区二区在线观看| 尾随美女入室| 亚洲精品成人久久久久久| 观看美女的网站| 国产精品久久久久久久久免| 亚洲经典国产精华液单| 99久久中文字幕三级久久日本| av国产精品久久久久影院| 久久精品久久久久久久性| 中文资源天堂在线| 日韩欧美 国产精品| 亚洲内射少妇av| 91久久精品国产一区二区成人| 国产69精品久久久久777片| 亚洲国产最新在线播放| 别揉我奶头 嗯啊视频| 老司机影院毛片| 亚洲av.av天堂| 欧美性猛交╳xxx乱大交人| 欧美xxxx黑人xx丫x性爽| 蜜桃久久精品国产亚洲av| 91狼人影院| 成人黄色视频免费在线看| 97在线视频观看| 一本一本综合久久| 久久久久久久国产电影| 国产一区二区三区综合在线观看 | 久久99热这里只有精品18| 国产高清有码在线观看视频| 欧美日韩亚洲高清精品| 欧美最新免费一区二区三区| 久久99热这里只有精品18| 最近手机中文字幕大全| 精品久久久噜噜| 亚洲精品aⅴ在线观看| 午夜激情久久久久久久| 在现免费观看毛片| 国产在线男女| 2021天堂中文幕一二区在线观| 久久ye,这里只有精品| 日韩 亚洲 欧美在线| 欧美老熟妇乱子伦牲交| 欧美精品人与动牲交sv欧美| 草草在线视频免费看| 一级毛片我不卡| 国产v大片淫在线免费观看| 夜夜爽夜夜爽视频| 特级一级黄色大片| 自拍偷自拍亚洲精品老妇| 少妇丰满av| 日韩强制内射视频| 免费黄色在线免费观看| 欧美人与善性xxx| 国产精品无大码| 亚洲精品国产成人久久av| 国产69精品久久久久777片| 国产欧美日韩精品一区二区| 日韩伦理黄色片| 综合色丁香网| 精品人妻一区二区三区麻豆| 国产熟女欧美一区二区| 婷婷色麻豆天堂久久| 免费大片18禁| 蜜桃久久精品国产亚洲av| 夜夜看夜夜爽夜夜摸| 精品久久久久久电影网| 国产精品爽爽va在线观看网站| 午夜亚洲福利在线播放| 久久久久久久精品精品| 国产成人freesex在线| 黄色一级大片看看| av天堂中文字幕网| av.在线天堂| 熟女电影av网| 色网站视频免费| 男的添女的下面高潮视频| 日韩一本色道免费dvd| 婷婷色麻豆天堂久久| 黄片无遮挡物在线观看| 最近中文字幕2019免费版| 人妻夜夜爽99麻豆av| 欧美精品一区二区大全| 欧美人与善性xxx| 看黄色毛片网站| 一级毛片aaaaaa免费看小| 97在线人人人人妻| tube8黄色片| 国产成人a区在线观看| 日韩在线高清观看一区二区三区| 亚洲精品乱码久久久久久按摩| 婷婷色综合大香蕉| 偷拍熟女少妇极品色| 免费看日本二区| 伊人久久精品亚洲午夜| 亚洲精品色激情综合| 久久久久性生活片| 女人久久www免费人成看片| 各种免费的搞黄视频| 2021少妇久久久久久久久久久| 涩涩av久久男人的天堂| 国内少妇人妻偷人精品xxx网站| 99久久精品热视频| 免费在线观看成人毛片| 岛国毛片在线播放| 亚州av有码| 97精品久久久久久久久久精品| 一级毛片aaaaaa免费看小| 2022亚洲国产成人精品| 69av精品久久久久久| 99热这里只有是精品50| 五月伊人婷婷丁香| 大片免费播放器 马上看| 国产淫语在线视频| 日本色播在线视频| 久久久色成人| 久久久久久久久久成人| 国产精品国产三级国产av玫瑰| 亚洲精品日本国产第一区| 观看免费一级毛片| 精华霜和精华液先用哪个| 一个人观看的视频www高清免费观看| 九九爱精品视频在线观看| 欧美日韩综合久久久久久| 国产一区亚洲一区在线观看| 久久久久国产精品人妻一区二区| 女人久久www免费人成看片| 丰满人妻一区二区三区视频av| 国模一区二区三区四区视频| 国产精品麻豆人妻色哟哟久久| av福利片在线观看| 男女边吃奶边做爰视频| 欧美 日韩 精品 国产| 伦理电影大哥的女人| 自拍欧美九色日韩亚洲蝌蚪91 | 精品国产露脸久久av麻豆| 亚洲综合色惰| 中文乱码字字幕精品一区二区三区| 在线播放无遮挡| 纵有疾风起免费观看全集完整版| 久久久久久久精品精品| 又黄又爽又刺激的免费视频.| 只有这里有精品99| 美女cb高潮喷水在线观看| 麻豆成人av视频| 日韩精品有码人妻一区| 日本一二三区视频观看| h日本视频在线播放| xxx大片免费视频| 99热6这里只有精品| 色5月婷婷丁香| 搡女人真爽免费视频火全软件| 五月天丁香电影| 国产伦理片在线播放av一区| 大码成人一级视频| 久久久久久久国产电影| 日韩av不卡免费在线播放| 在线观看国产h片| 国产在线一区二区三区精| 亚洲精品aⅴ在线观看| 亚洲,一卡二卡三卡| 日韩视频在线欧美| 成人国产av品久久久| 亚洲av电影在线观看一区二区三区 | 精品视频人人做人人爽| 亚洲av欧美aⅴ国产| 一区二区三区四区激情视频| 欧美97在线视频| 日韩不卡一区二区三区视频在线| 国产美女午夜福利| 天堂中文最新版在线下载 | 中文字幕亚洲精品专区| 日日啪夜夜撸| 一本一本综合久久| 99久久九九国产精品国产免费| 爱豆传媒免费全集在线观看| 亚洲人成网站在线播| 欧美精品一区二区大全| 国产精品女同一区二区软件| 另类亚洲欧美激情| av免费在线看不卡| 极品教师在线视频| 精品熟女少妇av免费看| 欧美性猛交╳xxx乱大交人| 国产乱人偷精品视频| 久久精品人妻少妇| 老司机影院成人| 18+在线观看网站| 免费av毛片视频| 国产片特级美女逼逼视频| 老师上课跳d突然被开到最大视频| h日本视频在线播放| 精品久久国产蜜桃| 成年av动漫网址| 国产欧美日韩精品一区二区| 成人黄色视频免费在线看| 国产精品一及| 亚洲精品中文字幕在线视频 | 国产黄频视频在线观看| 久久ye,这里只有精品| a级毛色黄片| 人妻制服诱惑在线中文字幕| 七月丁香在线播放| 国产美女午夜福利| h日本视频在线播放| 国产 一区精品| 免费观看性生交大片5| 国产综合懂色| 亚洲欧美精品专区久久| 狂野欧美白嫩少妇大欣赏| 成人亚洲精品一区在线观看 | 免费看不卡的av| 日韩欧美精品v在线| 亚洲aⅴ乱码一区二区在线播放| 超碰97精品在线观看| av福利片在线观看| 欧美一区二区亚洲| 中文字幕久久专区| 欧美人与善性xxx| 亚洲av福利一区| 国产中年淑女户外野战色| 亚洲国产精品999| 亚洲av成人精品一二三区| 婷婷色综合www| 精品少妇黑人巨大在线播放| 一级毛片电影观看| 蜜桃久久精品国产亚洲av| 欧美精品人与动牲交sv欧美| 久久午夜福利片| 国产精品不卡视频一区二区| av又黄又爽大尺度在线免费看| 国产精品一二三区在线看| 欧美区成人在线视频| 高清日韩中文字幕在线| 久久久久国产精品人妻一区二区| 日本免费在线观看一区| 搡老乐熟女国产| 最近最新中文字幕大全电影3| 成人亚洲欧美一区二区av| 久久韩国三级中文字幕| 在线观看人妻少妇| 三级国产精品欧美在线观看| 国产精品女同一区二区软件| 久久久久久久久久成人| 免费不卡的大黄色大毛片视频在线观看| 97精品久久久久久久久久精品| 久久女婷五月综合色啪小说 | 91久久精品国产一区二区成人| av专区在线播放| 国产免费一级a男人的天堂| 18禁裸乳无遮挡动漫免费视频 | 你懂的网址亚洲精品在线观看| 黄色怎么调成土黄色| 大又大粗又爽又黄少妇毛片口| 亚洲第一区二区三区不卡| 久久韩国三级中文字幕| 在线观看一区二区三区激情| 搞女人的毛片| 久久久精品欧美日韩精品| av在线老鸭窝| 看非洲黑人一级黄片| 人妻少妇偷人精品九色| 老司机影院毛片| 18禁在线播放成人免费| 成年av动漫网址| 成人毛片a级毛片在线播放| 亚洲欧美日韩卡通动漫| 大又大粗又爽又黄少妇毛片口| 日本与韩国留学比较| 亚洲aⅴ乱码一区二区在线播放| 久久精品熟女亚洲av麻豆精品| 少妇人妻一区二区三区视频| 久久人人爽人人片av| 国产 一区精品| 毛片女人毛片| 免费看不卡的av| 少妇人妻久久综合中文| 国产v大片淫在线免费观看| 久久久午夜欧美精品| 欧美xxxx黑人xx丫x性爽| 亚洲性久久影院| 在线观看三级黄色| 在线观看人妻少妇| 搡女人真爽免费视频火全软件| 一级av片app| 亚洲激情五月婷婷啪啪| 国产精品熟女久久久久浪| 搞女人的毛片| av在线天堂中文字幕| 天堂网av新在线| 久久久成人免费电影| 免费看av在线观看网站| 久久99热6这里只有精品| h日本视频在线播放| 人妻夜夜爽99麻豆av| 永久免费av网站大全| 亚洲国产精品专区欧美| 亚洲国产最新在线播放| 又黄又爽又刺激的免费视频.| 亚洲欧美精品自产自拍| 中文精品一卡2卡3卡4更新| 狂野欧美白嫩少妇大欣赏| 新久久久久国产一级毛片| 亚洲成人av在线免费| 天天躁夜夜躁狠狠久久av| 亚洲精品国产色婷婷电影| 亚洲精品乱码久久久v下载方式| 身体一侧抽搐| 国产成人福利小说| 亚洲av在线观看美女高潮| 交换朋友夫妻互换小说| 嫩草影院新地址| 亚洲国产最新在线播放| 在线亚洲精品国产二区图片欧美 | 干丝袜人妻中文字幕| 午夜视频国产福利| 日本与韩国留学比较| 亚洲国产精品成人综合色| 久久久久久久大尺度免费视频| 久久久久久久国产电影| 久久久亚洲精品成人影院| 超碰av人人做人人爽久久| 最近中文字幕高清免费大全6| 国产淫片久久久久久久久| 国产爽快片一区二区三区| 精品国产露脸久久av麻豆| 七月丁香在线播放| 精品一区在线观看国产| 国产一级毛片在线| 国产一区二区三区综合在线观看 | 人体艺术视频欧美日本| 国产黄片美女视频| 日韩av在线免费看完整版不卡| 国产亚洲91精品色在线| 国产高清有码在线观看视频| 亚洲人成网站高清观看| 欧美bdsm另类| 嫩草影院入口| 在线免费十八禁| 亚洲精品成人av观看孕妇| 一本久久精品| 亚洲成人一二三区av| 国产精品国产av在线观看| 日产精品乱码卡一卡2卡三| 我的老师免费观看完整版| 97超视频在线观看视频| 欧美成人a在线观看| 国产视频内射| 久久久久国产网址| 久久久成人免费电影| 色播亚洲综合网| 亚洲av福利一区| 麻豆成人午夜福利视频| 欧美日韩亚洲高清精品| 国产午夜精品一二区理论片| 国产成年人精品一区二区| 亚洲国产精品成人久久小说| 精品亚洲乱码少妇综合久久| 高清午夜精品一区二区三区| 国产精品一及| 亚洲va在线va天堂va国产| 日韩欧美一区视频在线观看 | 欧美激情国产日韩精品一区| 黄色一级大片看看| 五月玫瑰六月丁香| 国内精品宾馆在线| 老司机影院成人| 真实男女啪啪啪动态图| 一区二区三区精品91| 欧美另类一区| 国产日韩欧美在线精品| 七月丁香在线播放| 黄色视频在线播放观看不卡| 午夜免费观看性视频| 亚洲成人av在线免费| 久久久久久国产a免费观看| 亚洲三级黄色毛片| 男女国产视频网站| 久久人人爽av亚洲精品天堂 | 日韩亚洲欧美综合| 国产淫语在线视频| 在线观看三级黄色| 在线亚洲精品国产二区图片欧美 | 国产永久视频网站| 午夜免费鲁丝| 99久久精品一区二区三区| 亚洲伊人久久精品综合| 久久精品国产亚洲网站| 汤姆久久久久久久影院中文字幕| 日韩 亚洲 欧美在线| 一个人看的www免费观看视频| 性插视频无遮挡在线免费观看| 国产精品一区二区三区四区免费观看| 亚洲综合色惰| 简卡轻食公司| 中文天堂在线官网| 国内少妇人妻偷人精品xxx网站| 丝袜美腿在线中文| 免费看不卡的av| 激情 狠狠 欧美| 国产亚洲av片在线观看秒播厂| 成人毛片a级毛片在线播放| 丰满少妇做爰视频| 久久久久国产网址| 观看美女的网站| 乱系列少妇在线播放| 夜夜看夜夜爽夜夜摸| 亚洲欧美清纯卡通| 欧美区成人在线视频| 亚洲婷婷狠狠爱综合网| kizo精华| 国产 一区精品| 国产美女午夜福利| 99久久人妻综合| 亚洲成人精品中文字幕电影| 国产91av在线免费观看| 搡女人真爽免费视频火全软件| 国产日韩欧美亚洲二区| 99热网站在线观看| 高清毛片免费看| 精品国产乱码久久久久久小说| 日本爱情动作片www.在线观看| 国产成人精品福利久久| 丰满人妻一区二区三区视频av| 在线观看人妻少妇| 成年版毛片免费区| xxx大片免费视频| 国产探花在线观看一区二区| 99久久九九国产精品国产免费| av一本久久久久| 亚洲人成网站在线观看播放| 中国国产av一级| 欧美+日韩+精品| 午夜福利视频精品| 成人一区二区视频在线观看| 99re6热这里在线精品视频| 18禁裸乳无遮挡动漫免费视频 | 日韩欧美精品免费久久| www.av在线官网国产| 欧美性感艳星| 亚洲,欧美,日韩| 噜噜噜噜噜久久久久久91| 久久久精品免费免费高清| 中国国产av一级| 高清视频免费观看一区二区| 久久久精品94久久精品| 欧美高清成人免费视频www| 国产精品爽爽va在线观看网站| 亚洲人成网站高清观看| 久久ye,这里只有精品| 在线播放无遮挡| 69av精品久久久久久| 久久久久国产网址| 夫妻午夜视频| 国产精品一二三区在线看| 夜夜爽夜夜爽视频| 亚洲精品成人久久久久久| 高清毛片免费看| 成人国产麻豆网| 亚洲精品成人av观看孕妇| av在线观看视频网站免费| 最近手机中文字幕大全| 国产免费一级a男人的天堂| 亚洲精品成人av观看孕妇| 精品久久久久久久久亚洲| 中文字幕亚洲精品专区| 夫妻性生交免费视频一级片| 性插视频无遮挡在线免费观看| 国产亚洲av片在线观看秒播厂| 国产精品秋霞免费鲁丝片| 极品教师在线视频| www.色视频.com| 丝瓜视频免费看黄片| 啦啦啦在线观看免费高清www| 成人午夜精彩视频在线观看| 久久久久久久精品精品| 九九在线视频观看精品| 久久久久精品性色| 全区人妻精品视频| 精品少妇久久久久久888优播| 精品酒店卫生间| 亚洲va在线va天堂va国产| 交换朋友夫妻互换小说| 下体分泌物呈黄色| 又大又黄又爽视频免费| 亚洲天堂国产精品一区在线| 嫩草影院精品99| 伊人久久精品亚洲午夜| 听说在线观看完整版免费高清| 搞女人的毛片| 美女cb高潮喷水在线观看| 成年女人看的毛片在线观看| 99热6这里只有精品| av网站免费在线观看视频| 久久人人爽人人片av| 少妇被粗大猛烈的视频| 欧美成人a在线观看| 亚洲第一区二区三区不卡| 国产中年淑女户外野战色| 国产黄色免费在线视频| 国产精品久久久久久精品电影小说 | 国产人妻一区二区三区在| 高清午夜精品一区二区三区| 91在线精品国自产拍蜜月| 日韩成人伦理影院| 99久久九九国产精品国产免费| 综合色丁香网| 联通29元200g的流量卡| 一级毛片 在线播放| 在线免费观看不下载黄p国产| 大陆偷拍与自拍| 在线观看免费高清a一片| 日韩在线高清观看一区二区三区| 亚洲av.av天堂| 91精品伊人久久大香线蕉| 国产午夜精品一二区理论片| 精品一区二区免费观看| 国产精品久久久久久精品电影| 晚上一个人看的免费电影|