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

    基于主成分聚類分析的地下水化學(xué)分類方法

    2016-06-15 01:58:55周金龍賈瑞亮曾妍妍
    地下水 2016年1期
    關(guān)鍵詞:分類法水化學(xué)組分

    孟 奇,周金龍,賈瑞亮,曾妍妍

    (新疆農(nóng)業(yè)大學(xué)水利與土木工程學(xué)院,新疆 烏魯木齊 830052)

    ?

    基于主成分聚類分析的地下水化學(xué)分類方法

    孟奇,周金龍,賈瑞亮,曾妍妍

    (新疆農(nóng)業(yè)大學(xué)水利與土木工程學(xué)院,新疆 烏魯木齊 830052)

    [摘要]根據(jù)新疆巴里坤縣漢水泉地區(qū)20組承壓水水樣的K+、Na+、Ca2+、Mg2+、SO42-、Cl-和HCO3-等7個(gè)化學(xué)指標(biāo)檢測(cè)數(shù)據(jù),采用舒卡列夫分類法進(jìn)行地下水水化學(xué)類型分類。結(jié)果顯示,部分同類型之間組分含量差異較大而不同類之間組分含量差異較小,分析認(rèn)為是舒卡列夫分類法劃分界限人為性太強(qiáng)導(dǎo)致的。因此,我們采用主成分聚類分析法代替舒卡列夫分類法,對(duì)水化學(xué)類型進(jìn)行分類。本文借助SPSS軟件進(jìn)行主成分分析,提取2個(gè)主成分,將主成分得分作為指標(biāo)進(jìn)行聚類分析,最終將20組水樣劃分為6類。主成分聚類法通過(guò)消除變量之間的相關(guān)性使結(jié)果更加客觀合理。

    [關(guān)鍵詞]主成分分析;聚類分析;地下水化學(xué)分類

    正確合理地劃分地下水化學(xué)類型有助于分析地下水化學(xué)特征及其成因。目前地下水水化學(xué)分類方法有多種,其中舒卡列夫分類法是應(yīng)用較為廣泛的方法之一,但是此方法存在一些不足。

    漢水泉地區(qū)位于新疆維吾爾自治區(qū)巴里坤哈薩克自治縣北部,北部為大哈甫提克山、呼洪得雷山,南部為白依山。研究區(qū)內(nèi)海拔460~1 700 m,地形較平坦開(kāi)闊,地形坡降約為40‰~90‰,面積為1 482 km2[1]。該地區(qū)基本無(wú)人類工程活動(dòng),人類對(duì)地下水的擾動(dòng)很少。

    本文采用主成分聚類分析,將漢水泉地區(qū)20組承壓水水樣劃分為6類,通過(guò)與舒卡列夫分類法比對(duì),發(fā)現(xiàn)該方法分類結(jié)果更加客觀準(zhǔn)確,對(duì)于其他研究區(qū)的地下水水化學(xué)分類也有一定參考價(jià)值。

    1研究區(qū)水文地質(zhì)條件概況

    漢水泉地區(qū)內(nèi)無(wú)常年性地表水流,研究區(qū)內(nèi)及附近分布零星的泉點(diǎn),以及在融雪季節(jié)和夏季暴雨過(guò)后在溝谷中可形成暫時(shí)性的地表水流。研究區(qū)第四系、新近系形成了許多層面裂隙、構(gòu)造斷裂和孔隙等空隙,為地下水的儲(chǔ)存、循環(huán)和運(yùn)移提供了良好的場(chǎng)所,賦存基巖裂隙水、第四系松散巖類孔隙潛水和新近系碎屑巖類孔隙裂隙承壓水。研究區(qū)內(nèi)地下水類型為多層結(jié)構(gòu),即上覆第四系松散巖類孔隙潛水,下伏碎屑巖類孔隙裂隙承壓水及基巖裂隙水。

    研究區(qū)地下水補(bǔ)給來(lái)源為:山區(qū)為基巖裂隙水補(bǔ)給,平原區(qū)主要為大氣降水入滲、側(cè)向徑流補(bǔ)給。山區(qū)地下水主要以基巖裂隙泉水和側(cè)向徑流為主要排泄方式。平原區(qū)以蒸發(fā)、植被蒸騰、自流井溢出為主要排泄方式。

    2 地下水化學(xué)分類方法對(duì)比

    2.1舒卡列夫分類法

    舒卡列夫分類方法是常用的地下水化學(xué)分類方法[2]。趙麗等(2015)通過(guò)對(duì)研究區(qū)部分表層土壤和地下水取樣進(jìn)行測(cè)試分析,結(jié)合舒卡列夫分類法對(duì)研究區(qū)地下水水化學(xué)類型進(jìn)行分類,對(duì)地下水進(jìn)行了質(zhì)量綜合評(píng)價(jià)并分析了污染原因,結(jié)果表明,研究區(qū)地下水質(zhì)量不容樂(lè)觀[3]。尹恒等(2012)對(duì)德陽(yáng)市城市規(guī)劃區(qū)地下水水樣水化學(xué)分析結(jié)果,在采用舒卡列夫分類法水化學(xué)分類的基礎(chǔ)上,對(duì)研究區(qū)地下水水化學(xué)特征以及空間分布進(jìn)行了探究[4]。陸海燕等(2014)在已有資料的基礎(chǔ)上運(yùn)用舒卡列夫分類法對(duì)通州-龍旺莊地區(qū)地下水水化學(xué)類型進(jìn)行劃分,對(duì)水質(zhì)演化規(guī)律進(jìn)行了探究并分析了水質(zhì)惡化的原因[5]。郭春艷等(2014)運(yùn)用舒卡列夫法在對(duì)唐山市平原區(qū)地下水進(jìn)行了基本分類的基礎(chǔ)上,通過(guò)對(duì)研究區(qū)化學(xué)組分的濃度和分布規(guī)律的研究,更進(jìn)一步探究了研究區(qū)地下水環(huán)境的影響因素,并為該研究區(qū)地下水資源保護(hù)及開(kāi)發(fā)利用提供了依據(jù)[6]。

    舒卡列夫分類法采用地下水6種主要離子Na+(K+合并于Na+中)、Ca2+、Mg2+、SO42-、Cl-和HCO3-(CO32-合并于HCO3-中),以及TDS劃分地下水化學(xué)類型。規(guī)定含量大于25%毫克當(dāng)量的陰離子與陽(yáng)離子進(jìn)行組合,將地下水化學(xué)類型劃分為49類。周金龍等(2003)認(rèn)為:該方法依據(jù)25%毫克當(dāng)量的劃分帶有一定人為性,會(huì)對(duì)分類結(jié)果造成很大影響,導(dǎo)致分類不盡合理[7]。表1是采用舒卡列夫分類法對(duì)漢水泉地區(qū)20組承壓水地下水水樣的分類結(jié)果。

    表1  漢水泉地區(qū)地下水舒卡列夫法分類結(jié)果

    2~5類地下水化學(xué)組分均值和變異系數(shù)見(jiàn)表2。由表2可以看出:變異系數(shù)≤0.25占28.57%,變異系數(shù)0.25~0.50占60.71%,變異系數(shù)0.50~0.75占7.15%,變異系數(shù)0.75~1.00占3.57%;平均變異系數(shù)為0.33。除4類外,同一類型的地下水化學(xué)組分含量變異系數(shù)較大(一般>0.25)。由表2還可以看出:2類和3類、4類和5類地下水化學(xué)組分差異較小。由此說(shuō)明,舒卡列夫法分類結(jié)果在很大程度上達(dá)不到同類組分含量差異較小、不同類組分含量差異較大的分類目的[7]。

    表2  舒卡列夫分類法不同類型各組分平均值和變異系數(shù)一覽表

    注:由于HK24-2與HK14-1單獨(dú)分為一類(1類和6類),故未列入本表;斜線左側(cè)為組分的平均含量(毫克當(dāng)量/升),斜線右側(cè)為組分含量的變異系數(shù)。

    2.2主成分聚類分類法

    2.2.1主成分分析

    主成分分析是將原來(lái)眾多具有相關(guān)性指標(biāo),重新組合成新的無(wú)相關(guān)的綜合指標(biāo)從而替代原來(lái)的指標(biāo)。選取漢水泉地區(qū)20組承壓水7個(gè)化學(xué)指標(biāo)(K+、Na+、Ca2+、Mg2+、SO42-、Cl-和HCO3-)進(jìn)行分析,選用數(shù)據(jù)單位為毫克當(dāng)量/升。通過(guò)運(yùn)用SPSS軟件因子分析模塊對(duì)數(shù)據(jù)進(jìn)行主成分分析,為聚類分析做準(zhǔn)備。

    主成分分析過(guò)程如下:

    1)數(shù)據(jù)相關(guān)性

    根據(jù)數(shù)據(jù)生成的相關(guān)矩陣進(jìn)行檢驗(yàn),一般認(rèn)為相關(guān)系數(shù)多數(shù)大于0.3,則認(rèn)為適合做主成分分析[8]。由表3知大多相關(guān)系數(shù)大于0.3,相關(guān)性較強(qiáng),則可以做主成分分析。KMO(Kaiser-Meyer-Olkin)檢驗(yàn)統(tǒng)計(jì)量是用于比較變量間簡(jiǎn)單相關(guān)系數(shù)和偏相關(guān)系數(shù)的指標(biāo),一般認(rèn)為低于0.5是不適合做主成分分析[8]。本文計(jì)算KMO值為0.514,說(shuō)明原始指標(biāo)中存在共線性,選用的數(shù)據(jù)適合做主成分分析。以上兩種檢驗(yàn)都判斷該組數(shù)據(jù)可以做主成分分析。

    表3 相關(guān)系數(shù)矩陣

    2)主成分的選取

    主成分選取采用累積貢獻(xiàn)率法,大多以累計(jì)方差貢獻(xiàn)率超過(guò)85%為標(biāo)準(zhǔn),但是在本文中,通過(guò)生成結(jié)果(表4)發(fā)現(xiàn)第三主成分與原始因子相關(guān)性不強(qiáng),故提取前兩個(gè)主成分,累積貢獻(xiàn)率為80.37%(表5),即前兩個(gè)主成分代表了原始成分的80.37%的信息量,其中第一主成分累積貢獻(xiàn)率為62.13%。

    表4  因子荷載矩陣

    表5  主成分特征值及方差貢獻(xiàn)率

    通過(guò)本文數(shù)據(jù)分析,得出因子荷載矩陣(表4)。根據(jù)數(shù)據(jù)知,主成分1代替了原始數(shù)據(jù)中的K+、Na+、Ca2+、Cl-和SO42-,其中Mg2+和HCO3-載荷相對(duì)較低,即主成分1對(duì)其解釋度較低。主成分2代替了Mg2+和HCO3-。

    3)主成分得分

    主成分得分的計(jì)算方法:將表4中第i列的每個(gè)元素分別除以第i個(gè)特征根的平方根,得到第i個(gè)主成分系數(shù)(表6)。

    表6  主成分系數(shù)

    根據(jù)表6得出前兩個(gè)主成分y1和y2線性組合為:

    y1=0.36x1+0.45x2+0.45x3+0.10x4+0.44x5+0.43x6-0.28x7;

    y2=-0.07x1+0.13x2+0.06x3+0.75x4+0.24x5-0.18x6+0.57x7

    其中:x1,x2,x3,x4,x5,x6, x7為對(duì)原始變量標(biāo)準(zhǔn)化后的變量。

    原始變量標(biāo)準(zhǔn)化采用Z-Score標(biāo)準(zhǔn)化,然后通過(guò)上式計(jì)算出主成分得分(表7)[9,10]。

    表7  主成分得分

    圖1 主成分聚類分析樹(shù)狀圖

    2.2.2 聚類分析

    用主成分得分替代原始數(shù)據(jù)進(jìn)行樣品聚類。通過(guò)SPSS聚類模塊,對(duì)20組數(shù)據(jù)的進(jìn)行系統(tǒng)聚類分析。

    聚類方法選用Ward法(離差平方和法)[11-13],必須選用歐式距離[14],將主成分分析中提取的兩個(gè)主成分得分作為聚類的新指標(biāo)。根據(jù)聚類分析的結(jié)果(圖1),將20個(gè)水樣點(diǎn)劃分為6類(表8)。(其中樣點(diǎn)HK14-1與HK22-1各自聚為一類,這兩個(gè)樣點(diǎn)的TDS分別為6139.8 mg/L與6414.3 mg/L,為其余樣點(diǎn)TDS的1.4~6.5倍)。

    根據(jù)上述分類計(jì)算各類離子含量平均值及變異系數(shù)(表9)。1、2、3和6類地下水化學(xué)組分均值和變異系數(shù)見(jiàn)表9。由表8可以看出:變異系數(shù)≤0.25占71.43%,變異系數(shù)0.25~0.50占28.57%,變異系數(shù)0.50~0.75占0.00%,變異系數(shù)0.75~1.00占0.00%;平均變異系數(shù)為0.21。同一類型的地下水化學(xué)組分含量變異系數(shù)明顯小于舒卡列夫分類法劃分的同一類型的地下水化學(xué)組分含量變異系數(shù)。由表9還可以看出:各類地下水化學(xué)組分含量的差異明顯大于舒卡列夫分類法劃分的同一類型的地下水化學(xué)組分含量差異。

    對(duì)比表9和表2的不同類型各組分平均值和變異系數(shù)發(fā)現(xiàn),主成分聚類分析法對(duì)地下水化學(xué)類型的劃分比較客觀合理。

    表8  主成分聚類分析結(jié)果

    表9 主成分聚類分析法不同類型各組分平均值和變異系數(shù)一覽表

    注:同表2。

    3 結(jié)語(yǔ)

    比較主成分聚類分析法和舒卡列夫分類法對(duì)漢水泉地區(qū)承壓水水化學(xué)的分類結(jié)果,得出以下結(jié)論:

    (1)舒卡列夫法分類法劃分地下水化學(xué)類型的界限人為性強(qiáng),在很大程度上達(dá)不到同類組分含量差異較小、不同類組分含量差異較大的分類目的。

    (2)主成分聚類法通過(guò)消除了變量之間的相關(guān)性使結(jié)果更加客觀合理,同一類型的地下水化學(xué)組分含量變異系數(shù)明顯小于舒卡列夫分類法劃分的同一類型的地下水化學(xué)組分含量變異系數(shù),不同類型地下水化學(xué)組分含量的差異明顯大于舒卡列夫分類法劃分的不同類型的地下水化學(xué)組分含量差異。

    參考文獻(xiàn)

    [1] 吳霞,吳津蓉,李巧,等.新疆漢水泉地區(qū)地下水環(huán)境背景值計(jì)算[J].人民黃河.2015,(1):83-86,90.

    [2] 張人權(quán),梁杏,靳孟貴,等.水文地質(zhì)學(xué)基礎(chǔ)(第六版)[M].北京:地質(zhì)出版社.2011:67.

    [3] 趙麗,王同桂,鄧春光,等.重慶市西部地區(qū)淺層地下水化學(xué)特征及污染現(xiàn)狀分析[J].地下水.2015,(1):70-72.

    [4] 尹恒,吳勇,高東東,等.德陽(yáng)市城市規(guī)劃區(qū)地下水化學(xué)特征分析[J].地下水.2012,(1):30-35.

    [5] 陸海燕,楊帆,辛寶東,等.北京大興迭隆起隱伏巖溶水化學(xué)特征及水質(zhì)變化原因[J].南水北調(diào)與水利科技.2014, (5):155-159.

    [6] 郭春艷,馬震,張兆吉,等.唐山市平原區(qū)淺層地下水環(huán)境特征研究[J].南水北調(diào)與水利科技.2014,12(4).

    [7]周金龍,董新光,蘭衛(wèi)松.分步聚類分析方法劃分地下水水化學(xué)類型[J].新疆農(nóng)業(yè)大學(xué)學(xué)報(bào).2003,26(2):72-75.

    [8] 管于.實(shí)用多元統(tǒng)計(jì)分析[M].杭州:浙江大學(xué)出版社.2011:188-213.

    [9] 何曉群.多元統(tǒng)計(jì)分析(第三版)[M].北京:中國(guó)人民大學(xué)出版社.2012:129-135.

    [10] 王曉鵬,曹廣超.基于多元統(tǒng)計(jì)和GIS的環(huán)境質(zhì)量評(píng)價(jià)研究[M].北京:科學(xué)出版社.2013:19.

    [11] 呂晉,鄔紅娟,林濟(jì)東,等.主成分及聚類分析在水生態(tài)系統(tǒng)區(qū)劃中的應(yīng)用[J].武漢大學(xué)學(xué)報(bào):理學(xué)版.2005,51(4): 461-466.

    [12] 蔣群,許光泉,梁修雨.主成分和聚類分析應(yīng)用于淮南礦區(qū)地下水水質(zhì)評(píng)價(jià)[J].能源環(huán)境保護(hù).2007,21(2):51-53.

    [13] 魯守剛,鄧婭敏,張美雁,等. 基于多元統(tǒng)計(jì)分析的高砷地下水水化學(xué)特征研究——以內(nèi)蒙古杭錦后旗為例[J].中國(guó)農(nóng)村水利水電.2014,(12):179-182.

    [14] 胡雷芳.五種常用系統(tǒng)聚類分析方法及其比較[J].浙江統(tǒng)計(jì).2007,(4):11-13.

    Groundwater Chemical Classification Method Based on Principal Component Cluster Analysis

    MENGQi,ZHOU Jin-long,JIA Rui-liang,ZENG Yan-yan

    (College of Water Conservancy and Civil Engineering, Xinjiang Agricultural University, Urumqi 830052, Xinjiang)

    Abstract:According to the determination result of seven chemical indexes including K+, Na+, Ca2+, Mg2+, SO42-, Cl- and HCO3- of 20 groups of confined groundwater samples in Hanshuiquan region of Balikun County in Xinjiang , Shkalev classification is applied to classification of hydrochemical types. Results show that the difference among the same type of samples is relative greater and the difference among different type of samples is relative smaller which is probably due to certain artificial factors of Shkalev classification which can lead to the inaccuracy of classification results. Therefore, principal component cluster analysis is applied to classification of hydrochemical types instead of Shkalev classification. Using principal component analysis with SPSS software,two principal components were extracted in our study. The 20 groups of samples were divided into six categories according to the principal components scores which are regarded as the metric of cluster analysis. The principal component cluster analysis eliminates the correlation between variables and makes the results more objective and reasonable.

    Key words:Principal component analysis;cluster analysis and groundwater chemical classification

    [收稿日期]2015-07-31

    [基金項(xiàng)目]中國(guó)地質(zhì)調(diào)查局“西北地區(qū)主要城市地下水污染調(diào)查評(píng)價(jià)”(1212011220982)和新疆自治區(qū)水文學(xué)及水資源重點(diǎn)學(xué)科基金(xjswszyzdxk20101202)

    [作者簡(jiǎn)介]孟奇(1989-),男,江蘇徐州人,在讀碩士研究生,主攻方向:地下水利用與保護(hù)。[通訊作者]周金龍(1964-),男,浙江龍游人,教授,主要從事灌區(qū)土壤水—鹽運(yùn)移監(jiān)測(cè)與模擬、干旱區(qū)地下水利用與保護(hù)等教學(xué)與科研工作。

    [中圖分類號(hào)]P641.13

    [文獻(xiàn)標(biāo)識(shí)碼]A

    [文章編號(hào)]1004-1184(2016)01-0005-04

    猜你喜歡
    分類法水化學(xué)組分
    幸福湖泉域地下水水化學(xué)特征及其演化機(jī)制
    人民黃河(2025年1期)2025-01-16 00:00:00
    組分分發(fā)管理系統(tǒng)在天然氣計(jì)量的應(yīng)用
    一種難溶難熔未知組分板材的定性分析
    秦皇島北部地表水和地下水水化學(xué)特征及水質(zhì)評(píng)價(jià)
    分類法在高中化學(xué)中的應(yīng)用
    黑順片不同組分對(duì)正常小鼠的急性毒性
    中成藥(2018年8期)2018-08-29 01:28:26
    金雀花中黃酮苷類組分鑒定及2種成分測(cè)定
    中成藥(2018年2期)2018-05-09 07:20:09
    K 近鄰分類法在巖屑數(shù)字圖像巖性分析中的應(yīng)用
    錄井工程(2017年1期)2017-07-31 17:44:42
    衡水市新近系館陶組地?zé)崴瘜W(xué)特征與分布規(guī)律研究
    筠連煤層氣井產(chǎn)出水化學(xué)特征及意義
    内射极品少妇av片p| 99久久精品一区二区三区| 国产欧美亚洲国产| 亚洲欧洲日产国产| 久久久久久久久久人人人人人人| 亚洲四区av| 黑人猛操日本美女一级片| 免费少妇av软件| 狂野欧美白嫩少妇大欣赏| 男女国产视频网站| 日本av免费视频播放| av.在线天堂| 精品一品国产午夜福利视频| 超碰97精品在线观看| 男人爽女人下面视频在线观看| 最近中文字幕高清免费大全6| 日韩精品有码人妻一区| 两个人的视频大全免费| 九草在线视频观看| 少妇被粗大的猛进出69影院 | 熟女人妻精品中文字幕| 熟女人妻精品中文字幕| 91久久精品国产一区二区成人| 亚洲精品国产成人久久av| 亚洲精品日韩av片在线观看| 亚洲精品亚洲一区二区| 蜜桃在线观看..| 在线免费观看不下载黄p国产| 精品久久久精品久久久| 精品一区在线观看国产| 久久久久久伊人网av| 亚洲一区二区三区欧美精品| 99视频精品全部免费 在线| 国产91av在线免费观看| 国产精品久久久久久精品电影小说| 91精品伊人久久大香线蕉| 国产精品一区二区在线不卡| 在线观看三级黄色| 男女国产视频网站| 亚洲欧洲日产国产| 免费观看a级毛片全部| h日本视频在线播放| 看十八女毛片水多多多| 精品久久久噜噜| 秋霞在线观看毛片| 热99国产精品久久久久久7| 婷婷色综合大香蕉| 永久网站在线| 国产成人一区二区在线| 蜜桃久久精品国产亚洲av| 在线 av 中文字幕| 午夜免费鲁丝| 桃花免费在线播放| 国产黄频视频在线观看| 伦理电影大哥的女人| 精品人妻一区二区三区麻豆| 黑人巨大精品欧美一区二区蜜桃 | 久久国产乱子免费精品| av在线app专区| 免费观看的影片在线观看| 久久久久久久精品精品| 视频中文字幕在线观看| 久久99蜜桃精品久久| 日韩,欧美,国产一区二区三区| 精品99又大又爽又粗少妇毛片| 下体分泌物呈黄色| 一级毛片久久久久久久久女| 色吧在线观看| 久久久久久久久大av| 男女边摸边吃奶| 国产成人精品久久久久久| 亚洲成人av在线免费| h视频一区二区三区| 久久久精品免费免费高清| 一级爰片在线观看| 自线自在国产av| 日韩 亚洲 欧美在线| 欧美少妇被猛烈插入视频| 亚洲av成人精品一二三区| 如日韩欧美国产精品一区二区三区 | 妹子高潮喷水视频| 亚州av有码| 久久久久久久久大av| a级毛片在线看网站| 99精国产麻豆久久婷婷| 高清毛片免费看| 国产精品免费大片| 搡女人真爽免费视频火全软件| 国产午夜精品一二区理论片| 一级毛片我不卡| a级毛片免费高清观看在线播放| 国产精品秋霞免费鲁丝片| 少妇人妻一区二区三区视频| h日本视频在线播放| 国产精品人妻久久久久久| 美女视频免费永久观看网站| 七月丁香在线播放| 美女主播在线视频| 日韩不卡一区二区三区视频在线| 久久久久精品久久久久真实原创| 国产精品无大码| 国产精品国产三级专区第一集| 91午夜精品亚洲一区二区三区| 亚洲欧美精品专区久久| 欧美亚洲 丝袜 人妻 在线| 天天操日日干夜夜撸| 18+在线观看网站| 久久久久久久久久成人| 黑丝袜美女国产一区| 色94色欧美一区二区| 国产精品.久久久| 久久这里有精品视频免费| 成人美女网站在线观看视频| 久久婷婷青草| 国产日韩欧美视频二区| 国产精品久久久久久av不卡| 国产一区亚洲一区在线观看| 亚洲av成人精品一二三区| 成人毛片60女人毛片免费| 亚洲国产日韩一区二区| 免费人成在线观看视频色| 日韩在线高清观看一区二区三区| 久久久久久久国产电影| 久久精品久久久久久久性| 欧美丝袜亚洲另类| a 毛片基地| 国产精品久久久久成人av| 少妇猛男粗大的猛烈进出视频| 国产精品人妻久久久影院| 2022亚洲国产成人精品| 午夜福利在线观看免费完整高清在| 国产伦理片在线播放av一区| 国产欧美日韩一区二区三区在线 | 男女啪啪激烈高潮av片| www.av在线官网国产| 色94色欧美一区二区| av有码第一页| 欧美激情极品国产一区二区三区 | 成人国产麻豆网| 最新的欧美精品一区二区| 中文天堂在线官网| 夫妻性生交免费视频一级片| 黄色日韩在线| 亚洲av日韩在线播放| 亚洲成人av在线免费| 亚洲国产精品一区二区三区在线| 亚洲欧美一区二区三区黑人 | kizo精华| 欧美老熟妇乱子伦牲交| 欧美一级a爱片免费观看看| 久久久a久久爽久久v久久| 国产色爽女视频免费观看| 国产午夜精品一二区理论片| 插逼视频在线观看| 成人黄色视频免费在线看| 女人精品久久久久毛片| 久久人人爽av亚洲精品天堂| 校园人妻丝袜中文字幕| 91在线精品国自产拍蜜月| 少妇猛男粗大的猛烈进出视频| www.色视频.com| 免费观看av网站的网址| 一个人免费看片子| 性色avwww在线观看| 国产精品偷伦视频观看了| 国产淫语在线视频| 女性被躁到高潮视频| 夫妻午夜视频| 97在线视频观看| 自拍偷自拍亚洲精品老妇| 岛国毛片在线播放| 国产高清三级在线| 久久99热6这里只有精品| 天堂俺去俺来也www色官网| 亚洲精品中文字幕在线视频 | av国产久精品久网站免费入址| 在线观看免费高清a一片| 久久人人爽人人片av| 99久久精品一区二区三区| 我的老师免费观看完整版| 一本色道久久久久久精品综合| av卡一久久| 秋霞在线观看毛片| 久久久国产欧美日韩av| 国产伦理片在线播放av一区| 亚洲欧美一区二区三区黑人 | 国产真实伦视频高清在线观看| 国产一区亚洲一区在线观看| 久久久久久久精品精品| 在现免费观看毛片| 亚洲国产欧美日韩在线播放 | 中文欧美无线码| 欧美精品国产亚洲| 噜噜噜噜噜久久久久久91| 精品少妇内射三级| 午夜激情福利司机影院| 午夜精品国产一区二区电影| 日韩一区二区视频免费看| 插逼视频在线观看| 免费播放大片免费观看视频在线观看| 成人漫画全彩无遮挡| 观看美女的网站| 久久国内精品自在自线图片| 国产精品久久久久久av不卡| 亚洲av免费高清在线观看| 美女cb高潮喷水在线观看| 三级国产精品片| 日韩免费高清中文字幕av| 五月天丁香电影| 国产欧美日韩一区二区三区在线 | 少妇猛男粗大的猛烈进出视频| 色视频在线一区二区三区| 91在线精品国自产拍蜜月| 欧美日韩国产mv在线观看视频| 色婷婷久久久亚洲欧美| 中文字幕精品免费在线观看视频 | 97超碰精品成人国产| 你懂的网址亚洲精品在线观看| 人人妻人人澡人人爽人人夜夜| 久久久国产精品麻豆| 99精国产麻豆久久婷婷| 寂寞人妻少妇视频99o| 免费人妻精品一区二区三区视频| 国产精品无大码| 啦啦啦在线观看免费高清www| 一级毛片久久久久久久久女| 久久久久久久精品精品| 啦啦啦视频在线资源免费观看| 国产亚洲精品久久久com| 狂野欧美激情性bbbbbb| 久久久久久人妻| 日本欧美视频一区| 91精品一卡2卡3卡4卡| 国产精品一区二区三区四区免费观看| 夜夜看夜夜爽夜夜摸| 两个人的视频大全免费| 国产熟女午夜一区二区三区 | 久久99精品国语久久久| 国产成人精品久久久久久| 久久久午夜欧美精品| 精品久久久久久电影网| 噜噜噜噜噜久久久久久91| 午夜福利影视在线免费观看| 一本大道久久a久久精品| 国产精品人妻久久久影院| 黑人猛操日本美女一级片| a级毛片免费高清观看在线播放| 一级a做视频免费观看| a级毛片免费高清观看在线播放| 日日爽夜夜爽网站| 秋霞在线观看毛片| 久久久久久久久久久丰满| 国产欧美日韩一区二区三区在线 | 久久久久视频综合| 人人妻人人澡人人看| 国产精品国产三级国产av玫瑰| 2022亚洲国产成人精品| 蜜臀久久99精品久久宅男| 最黄视频免费看| 国产精品无大码| 国产精品无大码| 交换朋友夫妻互换小说| av播播在线观看一区| 啦啦啦中文免费视频观看日本| 一级毛片久久久久久久久女| 日韩强制内射视频| 最近中文字幕高清免费大全6| 91精品伊人久久大香线蕉| 国产老妇伦熟女老妇高清| 国产精品蜜桃在线观看| 老女人水多毛片| 久久婷婷青草| 亚洲av二区三区四区| av福利片在线观看| 伊人亚洲综合成人网| 国产欧美日韩综合在线一区二区 | 国精品久久久久久国模美| 一本一本综合久久| 在线观看人妻少妇| 乱系列少妇在线播放| av在线app专区| 久久久欧美国产精品| 国产黄片美女视频| 欧美三级亚洲精品| av在线观看视频网站免费| 七月丁香在线播放| a级毛片免费高清观看在线播放| 99久久精品国产国产毛片| 男人狂女人下面高潮的视频| 亚洲不卡免费看| a级一级毛片免费在线观看| 国产男女超爽视频在线观看| 亚洲精品视频女| 菩萨蛮人人尽说江南好唐韦庄| 国产极品粉嫩免费观看在线 | 免费av不卡在线播放| 国产成人精品久久久久久| 交换朋友夫妻互换小说| 久久久亚洲精品成人影院| 成人国产av品久久久| 成人国产av品久久久| 久久精品久久久久久噜噜老黄| 色婷婷av一区二区三区视频| 综合色丁香网| 在线免费观看不下载黄p国产| 欧美精品高潮呻吟av久久| 国产精品一区二区在线观看99| 国产精品一区二区在线观看99| 全区人妻精品视频| 日本av免费视频播放| 夫妻午夜视频| 黄色怎么调成土黄色| 大香蕉97超碰在线| 国产老妇伦熟女老妇高清| 久久久久久久亚洲中文字幕| 一个人看视频在线观看www免费| 91久久精品电影网| 一区在线观看完整版| 久久久精品94久久精品| 国产在视频线精品| 日日啪夜夜爽| 亚洲成色77777| 国产成人91sexporn| 美女福利国产在线| 精品人妻熟女av久视频| 欧美成人午夜免费资源| 伦理电影大哥的女人| 2022亚洲国产成人精品| 欧美另类一区| 国产老妇伦熟女老妇高清| videos熟女内射| 九草在线视频观看| 男的添女的下面高潮视频| 99久久精品热视频| 亚洲国产最新在线播放| 日日啪夜夜爽| 久久国产精品大桥未久av | 美女内射精品一级片tv| 亚洲精品国产av蜜桃| 国产欧美亚洲国产| 免费观看av网站的网址| 亚洲内射少妇av| 一二三四中文在线观看免费高清| 99久久人妻综合| 永久免费av网站大全| 久久影院123| 91精品国产九色| 人妻制服诱惑在线中文字幕| 日韩中字成人| 久久6这里有精品| av在线老鸭窝| 亚洲国产精品一区二区三区在线| 日韩三级伦理在线观看| 亚洲欧美一区二区三区国产| 成年女人在线观看亚洲视频| 国产精品无大码| 久久韩国三级中文字幕| 一个人免费看片子| 国产免费一级a男人的天堂| 国产精品国产三级国产av玫瑰| 国产免费视频播放在线视频| 中文在线观看免费www的网站| 久久鲁丝午夜福利片| 美女国产视频在线观看| 我要看黄色一级片免费的| 国产精品久久久久久精品古装| 美女中出高潮动态图| 久久人妻熟女aⅴ| 我的女老师完整版在线观看| 2018国产大陆天天弄谢| 国产成人a∨麻豆精品| 国产91av在线免费观看| 日韩,欧美,国产一区二区三区| 草草在线视频免费看| 我要看黄色一级片免费的| 成年女人在线观看亚洲视频| 这个男人来自地球电影免费观看 | 国产又色又爽无遮挡免| 黄片无遮挡物在线观看| av国产精品久久久久影院| 熟妇人妻不卡中文字幕| 国产av码专区亚洲av| 亚洲一级一片aⅴ在线观看| av女优亚洲男人天堂| av专区在线播放| 亚洲精品亚洲一区二区| 我要看日韩黄色一级片| 免费看av在线观看网站| 波野结衣二区三区在线| 我的女老师完整版在线观看| 国产片特级美女逼逼视频| 两个人免费观看高清视频 | 97精品久久久久久久久久精品| 日韩,欧美,国产一区二区三区| 丰满人妻一区二区三区视频av| 亚洲欧美中文字幕日韩二区| 97超碰精品成人国产| 国产极品粉嫩免费观看在线 | 亚洲国产欧美在线一区| 99热全是精品| 性高湖久久久久久久久免费观看| 国产国拍精品亚洲av在线观看| 美女视频免费永久观看网站| 不卡视频在线观看欧美| 九九久久精品国产亚洲av麻豆| 我要看黄色一级片免费的| 国产又色又爽无遮挡免| 婷婷色综合大香蕉| 我要看黄色一级片免费的| 91成人精品电影| 久久99一区二区三区| 人人澡人人妻人| 国产爽快片一区二区三区| 亚洲欧美一区二区三区黑人 | 一级片'在线观看视频| 日本黄色片子视频| 欧美日韩一区二区视频在线观看视频在线| 亚洲情色 制服丝袜| 在线观看人妻少妇| 国产精品一区二区性色av| 久久久久久久精品精品| a级一级毛片免费在线观看| 深夜a级毛片| 亚洲国产最新在线播放| 赤兔流量卡办理| 99热网站在线观看| 久久99热6这里只有精品| 欧美最新免费一区二区三区| 精品国产一区二区久久| 欧美xxxx性猛交bbbb| 欧美日本中文国产一区发布| 亚洲av福利一区| 在线亚洲精品国产二区图片欧美 | 水蜜桃什么品种好| av专区在线播放| 少妇人妻 视频| 丰满少妇做爰视频| 国产欧美日韩精品一区二区| 99热网站在线观看| 久久人人爽人人片av| 老女人水多毛片| 久久热精品热| 天堂8中文在线网| 极品教师在线视频| 黑丝袜美女国产一区| 欧美精品国产亚洲| 国产精品一区www在线观看| 国产在线一区二区三区精| 人妻夜夜爽99麻豆av| 欧美日韩视频高清一区二区三区二| videos熟女内射| 国产精品一区www在线观看| 欧美 日韩 精品 国产| 国产免费福利视频在线观看| 久久97久久精品| 国产精品人妻久久久久久| 日韩三级伦理在线观看| 亚洲av不卡在线观看| 午夜福利在线观看免费完整高清在| 91久久精品国产一区二区三区| av在线老鸭窝| 国产熟女午夜一区二区三区 | 国产免费一区二区三区四区乱码| 久热这里只有精品99| 九九久久精品国产亚洲av麻豆| 少妇熟女欧美另类| 三级经典国产精品| 国产精品免费大片| 精品亚洲乱码少妇综合久久| 国产永久视频网站| 国产av国产精品国产| 成人午夜精彩视频在线观看| 亚洲国产av新网站| 国产黄色视频一区二区在线观看| 最近中文字幕2019免费版| 亚洲国产精品一区三区| 涩涩av久久男人的天堂| 蜜臀久久99精品久久宅男| 美女内射精品一级片tv| 中国美白少妇内射xxxbb| 最近中文字幕高清免费大全6| 新久久久久国产一级毛片| 男人舔奶头视频| 久久午夜福利片| 精品少妇黑人巨大在线播放| 亚洲自偷自拍三级| 丝瓜视频免费看黄片| 欧美亚洲 丝袜 人妻 在线| 99久久精品国产国产毛片| 国产在线免费精品| 多毛熟女@视频| 免费黄频网站在线观看国产| 日本色播在线视频| 久久精品夜色国产| 国产av码专区亚洲av| 啦啦啦中文免费视频观看日本| 午夜视频国产福利| 亚洲精品乱久久久久久| 18禁动态无遮挡网站| 这个男人来自地球电影免费观看 | 国产 一区精品| 国产精品一区www在线观看| 九九在线视频观看精品| 久久影院123| tube8黄色片| 日韩中字成人| 久久精品国产亚洲网站| 欧美日韩一区二区视频在线观看视频在线| 久久久国产一区二区| 免费看日本二区| 国产亚洲一区二区精品| 亚洲欧美日韩东京热| 久久久国产一区二区| 久久精品夜色国产| h日本视频在线播放| 大片电影免费在线观看免费| 老熟女久久久| 看非洲黑人一级黄片| 妹子高潮喷水视频| 大话2 男鬼变身卡| 日韩免费高清中文字幕av| 伊人久久精品亚洲午夜| 久久人人爽人人爽人人片va| 老熟女久久久| 国产精品一区二区性色av| 综合色丁香网| 亚洲人成网站在线播| 男女免费视频国产| 青春草亚洲视频在线观看| 日本91视频免费播放| 看非洲黑人一级黄片| 色婷婷久久久亚洲欧美| 久久鲁丝午夜福利片| 日韩强制内射视频| 亚洲成人一二三区av| 亚洲欧美精品自产自拍| 毛片一级片免费看久久久久| 亚洲精华国产精华液的使用体验| 欧美精品国产亚洲| 99热这里只有是精品在线观看| 亚洲av欧美aⅴ国产| 春色校园在线视频观看| 亚洲无线观看免费| 国产精品三级大全| 午夜精品国产一区二区电影| 午夜av观看不卡| 女性被躁到高潮视频| 最近中文字幕2019免费版| 搡老乐熟女国产| 日韩欧美一区视频在线观看 | 午夜91福利影院| 18+在线观看网站| 一级a做视频免费观看| 久久久久久久久久久丰满| 国产精品成人在线| 久久久久精品性色| 天天操日日干夜夜撸| 高清不卡的av网站| 欧美激情极品国产一区二区三区 | 亚洲精品中文字幕在线视频 | 久久热精品热| 男人爽女人下面视频在线观看| 国产av码专区亚洲av| 精品国产国语对白av| 欧美+日韩+精品| √禁漫天堂资源中文www| 精品亚洲成a人片在线观看| 久久精品国产自在天天线| 国产成人aa在线观看| 亚洲欧美精品自产自拍| 日日摸夜夜添夜夜爱| 国产在线免费精品| 国产高清国产精品国产三级| 赤兔流量卡办理| 亚洲av男天堂| 亚洲美女搞黄在线观看| 国产亚洲5aaaaa淫片| 一级黄片播放器| 边亲边吃奶的免费视频| 亚洲欧美中文字幕日韩二区| 欧美日韩综合久久久久久| 久久人人爽av亚洲精品天堂| a级片在线免费高清观看视频| 少妇人妻久久综合中文| 香蕉精品网在线| 精华霜和精华液先用哪个| 午夜激情福利司机影院| 亚洲欧美精品自产自拍| 日日摸夜夜添夜夜爱| 午夜视频国产福利| 最近最新中文字幕免费大全7| av线在线观看网站| 欧美日韩视频精品一区| 99热全是精品| 国产精品福利在线免费观看| 成年人午夜在线观看视频| 日本黄色日本黄色录像| 久久6这里有精品| 精品一区二区三卡| 日韩制服骚丝袜av| 欧美成人午夜免费资源| 日韩av在线免费看完整版不卡| 丰满少妇做爰视频| 精品国产一区二区三区久久久樱花| 亚洲国产精品一区三区| 免费少妇av软件| 日韩av免费高清视频| 在线观看www视频免费| 我的老师免费观看完整版| 亚洲va在线va天堂va国产| av.在线天堂| 卡戴珊不雅视频在线播放| 夫妻性生交免费视频一级片| 亚洲综合精品二区| 国产在线视频一区二区|