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

    平原河網(wǎng)地表水空間分布特征

    2021-06-17 11:56:39祝新明李莉王翡吳偉超
    安徽農(nóng)業(yè)科學(xué) 2021年2期
    關(guān)鍵詞:河網(wǎng)嘉興市入境

    祝新明 李莉 王翡 吳偉超

    摘要 以布設(shè)在嘉興市主要水系32個(gè)水質(zhì)自動(dòng)監(jiān)測(cè)站2019年水質(zhì)監(jiān)測(cè)數(shù)據(jù)為基礎(chǔ),綜合運(yùn)用主成分分析和聚類分析等多元統(tǒng)計(jì)分析方法,對(duì)嘉興市平原河網(wǎng)地表水的空間分布特征進(jìn)行分析。結(jié)果表明,主成分分析表明嘉興市平原河網(wǎng)地表水主要受總磷、氨氮、總氮等指標(biāo)的影響,pH、溶解氧、電導(dǎo)率、高錳酸鹽指數(shù)等指標(biāo)的影響較小。聚類分析表明嘉興市平原河網(wǎng)水系主要分為3類,第1類以瀾溪塘為代表的河道,水質(zhì)較好;第2類為南部排入錢塘江、東部排入黃浦江、中部環(huán)嘉興市區(qū)河道,水質(zhì)次之;第3類為西部和北部入境河道,水質(zhì)最差;聚類分析結(jié)果與主成分綜合評(píng)判結(jié)果基本一致。

    關(guān)鍵詞 地表水;空間分布;主成分分析;聚類分析;嘉興市

    中圖分類號(hào) X832? 文獻(xiàn)標(biāo)識(shí)碼 A

    文章編號(hào) 0517-6611(2021)02-0038-04

    doi:10.3969/j.issn.0517-6611.2021.02.012

    開(kāi)放科學(xué)(資源服務(wù))標(biāo)識(shí)碼(OSID):

    Characteristics of Spatial Distribution of Surface Water in Plain River Networks—Taking Jiaxing City as an Example

    ZHU Xinming,LI Li,WANG Fei et al (Jiaxing Ecoenvironmental Monitoring Center of Zhejiang Province,Jiaxing,Zhejiang 314000)

    Abstract Based on the water quality monitoring data of 32 automatic water quality monitoring stations deployed in the main water systems of Jiaxing City in 2019,the spatial distribution characteristics of surface water in plain river networks of Jiaxing City were discussed by using the methods of principal component analysis (PCA) and cluster analysis (CA). PCA showed that surface water in plain river networks of Jiaxing City was mainly affected by total phosphorus, ammonia nitrogen and total nitrogen. pH, dissolved oxygen, conductivity and permanganate index had little influence. CA showed that the plain river networks of Jiaxing City were mainly divided into three clusters. The first cluster was the river represented by Lanxitang River, with good water quality;the second cluster was discharged into Qiantang River in the south, discharged into Huangpu River in the East and around Jiaxing Rivers in the middle,followed by water quality;the third cluster was imported river in the west and the north, with the worst water quality. The results of CA were basically consistent with the results of PCA evaluation.

    Key words Surface water;Spatial distribution;Principal component analysis;Cluster analysis;Jiaxing City

    近幾十年經(jīng)濟(jì)飛速發(fā)展和城市化造成的水體污染,嚴(yán)重制約了高質(zhì)量發(fā)展和美麗中國(guó)生態(tài)目標(biāo)的實(shí)現(xiàn)。影響河流水體水質(zhì)的因素是多方面的,通過(guò)對(duì)河流的持續(xù)監(jiān)測(cè)和評(píng)價(jià),摸清其空間分布特征,有利于管理部門分區(qū)分類開(kāi)展水環(huán)境管理[1-2]。聚類分析、主成分分析、方差分析、判別分析等多元統(tǒng)計(jì)方法,已廣泛應(yīng)用于識(shí)別不同水體時(shí)空分布特征和污染源[3-6]。黃金良等[7]借助主成分分析、聚類分析和一元線性回歸等方法,識(shí)別九龍江流域水質(zhì)的時(shí)空分布和影響因素;李典寶等[8]利用聚類和判別分析,研究了上海河流秋季水質(zhì)在空間分布的相識(shí)性、差異性,并表征了空間差異性的顯著水質(zhì)指標(biāo);唐玉蘭等[9]綜合運(yùn)用方差分析、聚類分析和多維多尺度分析方法,對(duì)渾河流域沈撫段水質(zhì)時(shí)空變化特征進(jìn)行了分析。

    嘉興市是典型的江南水鄉(xiāng),自2013年實(shí)施“五水共治”、剿滅劣V類水體、碧水行動(dòng)等一系列舉措,水環(huán)境質(zhì)量得到明顯改善,但根據(jù)2018年嘉興市生態(tài)環(huán)境局公布的地表水監(jiān)測(cè)斷面結(jié)果,IV~V類斷面占比依然過(guò)半[10-11]。因此筆者以分布在嘉興平原河網(wǎng)上的32個(gè)自動(dòng)監(jiān)測(cè)斷面監(jiān)測(cè)結(jié)果為基礎(chǔ),采用主成分分析和聚類分析等多元統(tǒng)計(jì)方法,對(duì)平原河網(wǎng)開(kāi)展空間特征分析,解析主要污染源,為地表水治理提供一定的科學(xué)依據(jù)。

    1 材料與方法

    1.1 研究區(qū)域概況

    嘉興市位于太湖流域東南部、杭嘉湖平原腹地,市境介于30°21′~30°01′N、120°18′~121°16′E,陸域東西長(zhǎng)92 km,南北寬76 km。市內(nèi)河道交錯(cuò),湖蕩眾多,河道總長(zhǎng)1.38萬(wàn)km,骨干河道57條,河道密度200 m/km2[12]。市域水系總體屬太湖流域,京杭大運(yùn)河縱貫其間,其余各水系如長(zhǎng)山河、海鹽塘、嘉善塘、平湖塘等河道與之相連,形成網(wǎng)狀水系。

    根據(jù)水質(zhì)研究的需要,確定了嘉興市域范圍內(nèi)32個(gè)地表水自動(dòng)監(jiān)測(cè)斷面,其中國(guó)家環(huán)境質(zhì)量監(jiān)測(cè)網(wǎng)站點(diǎn)13個(gè)、省級(jí)水站19個(gè),除湘家蕩外均設(shè)在市域縣區(qū)交界斷面上,具體如圖1所示。

    1.2 數(shù)據(jù)來(lái)源

    地表水自動(dòng)監(jiān)測(cè)站可實(shí)時(shí)連續(xù)監(jiān)測(cè),及時(shí)反映流域內(nèi)主要斷面水體水質(zhì)狀況的優(yōu)勢(shì)。該研究采用2019年32個(gè)地表水水質(zhì)監(jiān)測(cè)斷面的月均值數(shù)據(jù),監(jiān)測(cè)指標(biāo)包含水溫、pH、溶解氧、濁度、電導(dǎo)率、高錳酸鹽指數(shù)、氨氮、總磷、總氮共9個(gè)指標(biāo),數(shù)據(jù)描述見(jiàn)表1。

    1.3 評(píng)價(jià)方法

    1.3.1 主成分分析法。主成分分析法通過(guò)降維處理技術(shù),客觀確定各個(gè)指標(biāo)的權(quán)重,用少數(shù)指標(biāo)代替原來(lái)指標(biāo),使評(píng)價(jià)結(jié)果更加客觀。主要步驟如下:對(duì)相關(guān)原始數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化變換,求解相關(guān)系數(shù)矩陣,計(jì)算特征值和特征向量,解釋方差,確定主成分,對(duì)各主成分進(jìn)行評(píng)價(jià)[13-14]。

    1.3.2 聚類分析法。

    聚類分析是按照某些數(shù)量特征將觀察對(duì)象進(jìn)行空間分類的一種統(tǒng)計(jì)方法[15]。層次聚類分析又是其中一種常用水質(zhì)空間分析的方法,其實(shí)質(zhì)是根據(jù)變量或觀察對(duì)象的相似程度,以逐次聚合的方法,將最相似的對(duì)象結(jié)合在一起,直到聚成一類[16]。

    該研究在數(shù)據(jù)標(biāo)準(zhǔn)化的基礎(chǔ)上,采用常用的歐氏距離平方和離差平方和法開(kāi)展各監(jiān)測(cè)斷面的空間相似性分析[17]。應(yīng)用SPSS 22.0軟件對(duì)地表水自動(dòng)站監(jiān)測(cè)結(jié)果開(kāi)展主成分分析和聚類分析。

    2 結(jié)果與分析

    2.1 基于主成分分析法的水質(zhì)質(zhì)量評(píng)價(jià)

    在主成分分析法前,對(duì)數(shù)據(jù)的KMO檢驗(yàn)結(jié)果為0.579(>0.5),顯著性為0(Sig.<0.05),基本適合做主成分分析。對(duì)各水質(zhì)自動(dòng)監(jiān)測(cè)站監(jiān)測(cè)結(jié)果開(kāi)展主成分分析,結(jié)果見(jiàn)表2。以累計(jì)方差貢獻(xiàn)率大于85%的原則,提取前4個(gè)主成分因子(累計(jì)貢獻(xiàn)率達(dá)87.9%),可以代表原有監(jiān)測(cè)項(xiàng)目的絕大部分信息。

    第1主成分的方差貢獻(xiàn)率為40.5%,其中總氮、氨氮、總磷、溶解氧所占載荷因子較大,與第1主成分的相關(guān)系數(shù)的絕對(duì)值均超過(guò)0.6,其中氨氮、總磷、總氮呈正相關(guān),溶解氧呈負(fù)相關(guān),表明監(jiān)測(cè)點(diǎn)位所在河流營(yíng)養(yǎng)鹽濃度較高,而有機(jī)污染物的增加消耗了水中的溶解氧,反映人類活動(dòng)(生活污水排放、農(nóng)業(yè)面源污染、水產(chǎn)養(yǎng)殖等污染)的影響[18]。第2主成分的方差貢獻(xiàn)率為21.8%,其中電導(dǎo)率、高錳酸鹽指數(shù)的所占載荷因子較大,與第2主成分的相關(guān)系數(shù)分別為0.671和-0.667,表明上游引錢塘江水沖刷河道和南排工程開(kāi)閘放水引起的海水倒灌影響了所在水體的離子水平和河流的有機(jī)污染水平,其中電導(dǎo)率呈正相關(guān),高錳酸鹽指數(shù)呈負(fù)相關(guān)。第3主成分的方差貢獻(xiàn)率為13.4%,其中電導(dǎo)率所占的載荷因子最大,與第3主成分的相關(guān)系數(shù)為-0.574,反映監(jiān)測(cè)點(diǎn)位所在水體的離子水平。第4主成分的方差貢獻(xiàn)率為12.2%,其總磷所占的載荷因子最大,與第4主成分的相關(guān)系數(shù)為0.562,反映了水體中磷的污染水平。

    根據(jù)主成分載荷值除以對(duì)應(yīng)特征值的平方根,并與相應(yīng)的變量進(jìn)行線性組合可得到對(duì)應(yīng)各主成分的評(píng)價(jià)得分。各主成分評(píng)價(jià)得分與相應(yīng)特征值占提取主成分特征值之和的比重乘積[19],求和后即為各斷面的綜合得分,結(jié)果見(jiàn)表3。

    綜合評(píng)價(jià)得分反映了各斷面所在河流的污染水平,得分越高說(shuō)明污染越嚴(yán)重。由表3可知,嘉興市水質(zhì)大體呈出境水質(zhì)好于入境水質(zhì)。從入境看,入境斷面大麻渡口、渡船橋和王江涇綜合得分較高,表明上游京杭運(yùn)河杭州塘、鹽官下河和京杭運(yùn)河蘇州塘的入境水體污染水平較高,而晚村、新塍港、民主水文站綜合得分較低,表明湖州入境長(zhǎng)山河,蘇州入境的新塍塘和蘆墟塘水體污染水平較低。從出境看,楓南大橋、南臺(tái)頭閘一號(hào)橋綜合得分較高,表明其所在嘉善塘和海鹽塘出境水質(zhì)較差,而烏鎮(zhèn)北、池家浜水文站、清涼大橋和鹽官排澇樞紐綜合得分較低,表明其對(duì)應(yīng)的瀾溪塘、俞匯塘、清涼港和鹽官下河出境水質(zhì)得到改善。

    2.2 基于聚類分析地表水監(jiān)測(cè)斷面空間分布特征

    選擇pH、溶解氧、電導(dǎo)率、高錳酸鹽指數(shù)、氨氮、總磷、總氮等項(xiàng)目開(kāi)展層次聚類分析,分析前對(duì)數(shù)據(jù)標(biāo)準(zhǔn)化并進(jìn)行正態(tài)分布檢驗(yàn),結(jié)果峰度為-1.047~3.759,偏度為-0.691~1.546,除個(gè)別指標(biāo)外,數(shù)據(jù)大體呈正態(tài)分布。

    基于水質(zhì)指標(biāo)的監(jiān)測(cè)斷面的聚類分析結(jié)果見(jiàn)圖2。依據(jù)監(jiān)測(cè)斷面水體理化性質(zhì)的相似性,監(jiān)測(cè)斷面聚類分析結(jié)果在空間上可劃分為3類。第1類包括斜橋、大麻渡口、王店南梅、楓南大橋、王江涇、杭申公路橋、渡船橋、聯(lián)合橋8個(gè)監(jiān)測(cè)斷面,對(duì)應(yīng)鹽平塘、京杭運(yùn)河杭州塘入境、長(zhǎng)水塘、楓涇塘、京杭運(yùn)河蘇州塘入境、鹽官下河入境、上塘河入境、長(zhǎng)山河中段,水質(zhì)污染水平較高,主要污染物為耗氧型有機(jī)污染物。

    第2類包括尤甪村等20個(gè)斷面,其中尤甪村、南臺(tái)頭閘一號(hào)橋和長(zhǎng)山河大橋、長(zhǎng)山閘一號(hào)橋分別對(duì)應(yīng)為排入錢塘江的海鹽塘和長(zhǎng)山河;荒田浜、青陽(yáng)匯、小新村和紅旗塘大壩分別對(duì)應(yīng)為以上海塘、廣陳塘為骨干排入黃浦江河道和經(jīng)紅旗塘入浦的河道,其中長(zhǎng)山河和海鹽塘開(kāi)閘后受海水倒灌影響導(dǎo)致電導(dǎo)率水平較高;善西、大云、天凝、楊廟大橋、新塍大通、新文橋?qū)?yīng)為嘉善塘、大云市河、紅旗塘入口、三店塘、運(yùn)河嘉興段等環(huán)市區(qū)河道;鹽官排澇樞紐、上塘河排澇閘對(duì)應(yīng)為排入錢塘江的鹽官下河和上塘河;清涼大橋、烏鎮(zhèn)、新塍港、晚村對(duì)應(yīng)為清涼港、雙林塘、新塍塘和長(zhǎng)山河入境。第2類河道污染水平相較第1類有明顯降低,表明水體自入境后經(jīng)過(guò)各河道凈化后,消耗了一部分污染物。

    第3類包括民主水文站、池家浜水文站、湘家蕩、烏鎮(zhèn)北4個(gè)監(jiān)測(cè)斷面,表明蘆墟塘入境、俞匯塘、三店塘支流和瀾溪塘的水體總體污染水平較低;其中氨氮、總磷和總氮平均濃度低于其余各組,且水體中溶解氧濃度較高。

    聚類分析結(jié)果與主成分分析所得結(jié)果基本一致。嘉興市平原河網(wǎng)水質(zhì)大體分為3類,第1類為瀾溪塘、俞匯塘、蘆墟塘和湘家蕩,水質(zhì)最優(yōu);第2類為東部入黃浦江河道、南部排入錢塘江河道、中部環(huán)嘉興市區(qū)各河道,水質(zhì)次之;第3類為北部入境的京杭運(yùn)河蘇州塘(引長(zhǎng)濟(jì)太后因太湖水位抬升,王江涇出境改為入境)、西部入境的京杭運(yùn)河杭州塘、上塘河和鹽官下河,水質(zhì)較差。

    隨著“五水共治”“河長(zhǎng)制”“污水零直排區(qū)”創(chuàng)建等措施的推進(jìn),目前工業(yè)企業(yè)污水均已入網(wǎng)排放,嘉興平原河網(wǎng)各水系中西部和北部入境的京杭運(yùn)河杭州塘、鹽官下河、上塘河、京杭運(yùn)河蘇州塘均為受上游來(lái)水較差影響。市區(qū)內(nèi)主要河道主要受城鎮(zhèn)和農(nóng)村生活污水、農(nóng)業(yè)種植業(yè)、養(yǎng)殖業(yè)和航道污染等影響。河道兩側(cè)分布各式居民區(qū),部分尚未實(shí)現(xiàn)污水管網(wǎng)全覆蓋,或雖經(jīng)納管,但雨污分流不徹底,導(dǎo)致雨天生活污水外溢情況的發(fā)生。沿河大部農(nóng)村生活污水雖經(jīng)微動(dòng)力等終端處理設(shè)施集中處理,但處理后大多排入平原河網(wǎng)各水系。嘉興素有魚(yú)米之鄉(xiāng)的美譽(yù),根據(jù)2019年嘉興統(tǒng)計(jì)年鑒,現(xiàn)有耕地面積204 897 hm2[20],基本沿河分布,由于化肥利用率不高,大部分肥料隨著雨水沖刷,進(jìn)入河道,造成磷氮污染。畜禽養(yǎng)殖和水產(chǎn)養(yǎng)殖也易導(dǎo)致水體富營(yíng)養(yǎng)化污染。嘉興地勢(shì)平坦,平均海拔只有3.7? m,水體流動(dòng)性差,污染物易累積,水體自凈能力較差。河網(wǎng)密布,大多河道屬“杭申線”或“湖嘉申線”航道,航運(yùn)繁忙,船只來(lái)往航行攪動(dòng)底泥,引起水體渾濁導(dǎo)致水質(zhì)較差。城市化導(dǎo)致的水系變化劇

    烈[21],部分水系水流不暢容易導(dǎo)致徑流排泄能力不足。

    3 結(jié)論

    對(duì)2019年分布在嘉興市平原河網(wǎng)的32個(gè)水站水質(zhì)監(jiān)測(cè)結(jié)果開(kāi)展主成分分析表明,嘉興平原河網(wǎng)水質(zhì)主要污染物為氨氮、總磷、總氮,其中pH、溶解氧、高錳酸鹽指數(shù)和電導(dǎo)率的影響較小。

    聚類分析結(jié)果表明嘉興市平原河網(wǎng)水系主要可分為3類,其中瀾溪塘等河道水質(zhì)最優(yōu),水體污染較小,南部排入錢塘江、東部排入浦河道和中部環(huán)市區(qū)河道水質(zhì)次之,西部和北部等入境河道水質(zhì)最差。西部、北部等主要入境河道水質(zhì)較差,主要受上游來(lái)水水質(zhì)較差影響。市域河網(wǎng)各水系主要污染來(lái)源于城鎮(zhèn)和農(nóng)村生活污染源、農(nóng)業(yè)面源污染、養(yǎng)殖業(yè)、航道污染等。

    參考文獻(xiàn)

    [1]

    紅梅,劉文治,張全發(fā).多元統(tǒng)計(jì)方法在金水河水質(zhì)時(shí)空變化分析中的應(yīng)用[J].資源科學(xué),2009,31(3):429-434.

    [2] 陳帥,莫彩芬,李艷薔,等.洪湖水質(zhì)時(shí)空特征及污染驅(qū)動(dòng)力分析[J].環(huán)境污染與防治,2019,41(4):421-425.

    [3] 何萬(wàn)謙,黃金良.澳門半島近岸海域水質(zhì)時(shí)空變異分析[J].環(huán)境科學(xué),2010,31(3):606-610.

    [4] 汪紅軍,顏昌龍,李嗣新,等.洪湖水質(zhì)空間特異性及主導(dǎo)因子分析[J].中國(guó)環(huán)境監(jiān)測(cè),2012,28(3):72-75.

    [5] 錢天鳴.杭州市錢塘江干支流水質(zhì)多元統(tǒng)計(jì)分析[J].中國(guó)環(huán)境監(jiān)測(cè),2015,31(2):74-77.

    [6] 岳程鵬,李興.基于聚類分析的烏梁素海水質(zhì)因子研究[J].人民黃河,2019,41(5):67-70,76.

    [7] 黃金良,黃亞玲,李青生,等.流域水質(zhì)時(shí)空分布特征及其影響因素初析[J].環(huán)境科學(xué),2012,33(4):1098-1107.

    [8] 李典寶,張瑋,王麗卿,等.2013年上海市河流秋季水質(zhì)空間分布特征[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2015,31(1):50-58.

    [9] 唐玉蘭,項(xiàng)瑩雪,馬甜甜,等.基于多元統(tǒng)計(jì)分析方法的渾河流域沈撫段水質(zhì)時(shí)空特征[J].安全與環(huán)境學(xué)報(bào),2018,18(5):2008-2012.

    [10] 嘉興市生態(tài)環(huán)境局.嘉興市2013年環(huán)境狀況公報(bào)[EB/OL].(2014-06-06)[2020-01-15].http://www.jiaxing.gov.cn/art/2014/6/6/art_1555291_26762230.html.

    [11] 嘉興市生態(tài)環(huán)境局.嘉興市2018年環(huán)境狀況公報(bào)[EB/OL].(2019-06-18) [2020-01-15].http://www.jiaxing.gov.cn/art/2019/6/18/art_1555291_34702773.html.

    [12] 葉永強(qiáng).嘉興年鑒(2017)[M]. 北京:方志出版社,2017.

    [13] 郝海,蹤家峰.系統(tǒng)分析與評(píng)價(jià)方法[M].北京:經(jīng)濟(jì)科學(xué)出版社,2007.

    [14] 朱琳,王雅南,韓美,等.武水河水質(zhì)時(shí)空分布特征及污染成因的解析[J].環(huán)境科學(xué)學(xué)報(bào),2018,38(6):2150-2156.

    [15] 王昱,盧世國(guó),劉娟娟,等.春季枯水期黑河水體理化性質(zhì)的空間分布特征[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2019,35(4):433-441.

    [16] 周豐,郭懷成,黃凱,等.基于多元統(tǒng)計(jì)方法的河流水質(zhì)空間分析[J].水科學(xué)進(jìn)展,2007,18(4):544-551.

    [17] 陳永娟,龐樹(shù)江,耿潤(rùn)哲,等.北運(yùn)河水系主要污染物通量特征研究[J].環(huán)境科學(xué)學(xué)報(bào),2015,35(7):2167-2176.

    [18] 馬小雪,王臘春,廖玲玲.溫瑞塘河流域水體污染時(shí)空分異特征及污染源識(shí)別[J].環(huán)境科學(xué),2015,36(1):64-71.

    [19] 吉祝美,方里,張俊,等.主成分分析法在SPSS軟件中的操作及在河流水質(zhì)評(píng)價(jià)中的應(yīng)用[J].環(huán)境研究與監(jiān)測(cè),2012,25(4):68-73.

    [20] 嘉興市統(tǒng)計(jì)局.嘉興統(tǒng)計(jì)年鑒-2019[M]. 北京:中國(guó)統(tǒng)計(jì)出版社,2019.

    [21] 鄧曉軍,許有鵬,韓龍飛,等.城市化背景下嘉興市河流水系的時(shí)空變化[J].地理學(xué)報(bào),2016,71(1):75-85.

    猜你喜歡
    河網(wǎng)嘉興市入境
    基于小世界網(wǎng)絡(luò)的海河流域河網(wǎng)結(jié)構(gòu)及功能響應(yīng)
    嘉興市第一醫(yī)院
    嘉興市第一醫(yī)院
    嘉興市第一醫(yī)院
    基于PSR模型的上海地區(qū)河網(wǎng)脆弱性探討
    中國(guó)第三個(gè)國(guó)家級(jí)入境再制造示范區(qū)通過(guò)驗(yàn)收
    嘉興市:多措并舉推進(jìn)環(huán)境治理上臺(tái)階
    不同引水水源對(duì)平原河網(wǎng)影響分析
    MIKE11模型在城市河網(wǎng)生態(tài)調(diào)水工程中的應(yīng)用
    非法入境
    99久国产av精品| 人妻制服诱惑在线中文字幕| 91午夜精品亚洲一区二区三区| 精品久久久久久电影网 | 亚洲av成人精品一区久久| 国产精品一区二区三区四区免费观看| 亚洲欧美精品专区久久| 久久人人爽人人爽人人片va| 国产老妇伦熟女老妇高清| 国产精品永久免费网站| 亚洲欧美日韩东京热| 1000部很黄的大片| 亚洲一级一片aⅴ在线观看| 观看美女的网站| 日韩一区二区三区影片| 日韩欧美在线乱码| 天天躁夜夜躁狠狠久久av| 国产欧美另类精品又又久久亚洲欧美| 六月丁香七月| 久久精品国产亚洲av天美| 久久鲁丝午夜福利片| 免费看a级黄色片| 国产亚洲午夜精品一区二区久久 | 久久6这里有精品| 日本黄色片子视频| 熟妇人妻久久中文字幕3abv| 国产一区二区三区av在线| 中文字幕熟女人妻在线| 少妇熟女欧美另类| 91狼人影院| 国产女主播在线喷水免费视频网站 | 国产成人免费观看mmmm| 亚洲综合精品二区| 久久人妻av系列| .国产精品久久| 久久精品久久久久久噜噜老黄 | 熟女电影av网| 日韩欧美三级三区| 日本欧美国产在线视频| 成人一区二区视频在线观看| 亚洲第一区二区三区不卡| 色综合亚洲欧美另类图片| 亚洲aⅴ乱码一区二区在线播放| 啦啦啦啦在线视频资源| 91狼人影院| 人人妻人人澡欧美一区二区| 少妇裸体淫交视频免费看高清| 国产中年淑女户外野战色| 亚洲天堂国产精品一区在线| 国产av码专区亚洲av| 精品一区二区免费观看| 精品久久久久久久末码| 亚洲图色成人| 国产黄片美女视频| 高清日韩中文字幕在线| 最近中文字幕2019免费版| eeuss影院久久| 一夜夜www| 如何舔出高潮| 网址你懂的国产日韩在线| 91av网一区二区| 一个人免费在线观看电影| 久久久国产成人精品二区| 内射极品少妇av片p| 国产精品福利在线免费观看| 超碰av人人做人人爽久久| 桃色一区二区三区在线观看| 特级一级黄色大片| 亚洲成人中文字幕在线播放| 毛片一级片免费看久久久久| 国产精品蜜桃在线观看| 内射极品少妇av片p| 日本色播在线视频| 午夜福利在线观看吧| 午夜爱爱视频在线播放| 丝袜喷水一区| .国产精品久久| 午夜福利高清视频| 国产高清国产精品国产三级 | 日韩中字成人| 嫩草影院新地址| 婷婷色av中文字幕| 国产白丝娇喘喷水9色精品| 午夜a级毛片| 男女视频在线观看网站免费| 中文在线观看免费www的网站| 一级黄色大片毛片| 亚洲国产欧洲综合997久久,| 好男人在线观看高清免费视频| 午夜精品一区二区三区免费看| 永久免费av网站大全| 九色成人免费人妻av| 51国产日韩欧美| 看黄色毛片网站| 汤姆久久久久久久影院中文字幕 | 欧美成人一区二区免费高清观看| 午夜精品一区二区三区免费看| 精品一区二区三区人妻视频| 噜噜噜噜噜久久久久久91| 少妇人妻精品综合一区二区| 一个人看视频在线观看www免费| 汤姆久久久久久久影院中文字幕 | 男女视频在线观看网站免费| 久久精品国产亚洲av天美| 亚洲成人av在线免费| 九九久久精品国产亚洲av麻豆| 亚洲国产最新在线播放| 午夜精品一区二区三区免费看| 1000部很黄的大片| 美女xxoo啪啪120秒动态图| 成人鲁丝片一二三区免费| 精品少妇黑人巨大在线播放 | 日本一本二区三区精品| 久久欧美精品欧美久久欧美| 狠狠狠狠99中文字幕| 亚洲av电影不卡..在线观看| www.av在线官网国产| 在线观看av片永久免费下载| 97在线视频观看| 国产黄色视频一区二区在线观看 | 噜噜噜噜噜久久久久久91| 在线免费观看不下载黄p国产| 亚洲丝袜综合中文字幕| 日韩av在线免费看完整版不卡| 成年女人永久免费观看视频| 成年免费大片在线观看| 大香蕉久久网| 日日撸夜夜添| 国产精品嫩草影院av在线观看| 欧美一级a爱片免费观看看| 我的老师免费观看完整版| 亚洲熟妇中文字幕五十中出| 秋霞在线观看毛片| 日本熟妇午夜| 亚洲怡红院男人天堂| 男女视频在线观看网站免费| 亚洲在线观看片| 麻豆久久精品国产亚洲av| 国产成年人精品一区二区| 国语对白做爰xxxⅹ性视频网站| 亚洲三级黄色毛片| 日本与韩国留学比较| 日本爱情动作片www.在线观看| 精品久久久久久成人av| 久久精品91蜜桃| 国产黄a三级三级三级人| 最近的中文字幕免费完整| 中国美白少妇内射xxxbb| 五月伊人婷婷丁香| 国产片特级美女逼逼视频| 中国美白少妇内射xxxbb| 少妇被粗大猛烈的视频| 99在线视频只有这里精品首页| 久久久久久久久中文| 国产午夜精品久久久久久一区二区三区| 嫩草影院新地址| av黄色大香蕉| 婷婷色麻豆天堂久久 | 最近中文字幕2019免费版| av又黄又爽大尺度在线免费看 | 你懂的网址亚洲精品在线观看 | 欧美+日韩+精品| 看非洲黑人一级黄片| 国产精品久久视频播放| 国产在线男女| 久久亚洲精品不卡| 久久精品夜色国产| 精品不卡国产一区二区三区| 欧美一区二区国产精品久久精品| 亚洲不卡免费看| 国产免费福利视频在线观看| 视频中文字幕在线观看| 国产乱人偷精品视频| 国产免费福利视频在线观看| 国产黄色小视频在线观看| 久久精品国产亚洲av天美| 亚洲av日韩在线播放| 亚洲最大成人av| 中文字幕亚洲精品专区| 大话2 男鬼变身卡| 亚洲av一区综合| 国产在线一区二区三区精 | 久久韩国三级中文字幕| 欧美变态另类bdsm刘玥| av免费观看日本| 色哟哟·www| 久久这里有精品视频免费| 亚洲在久久综合| 欧美日韩在线观看h| 麻豆成人av视频| 亚洲成色77777| 村上凉子中文字幕在线| 91久久精品国产一区二区三区| 国产亚洲精品久久久com| 嘟嘟电影网在线观看| 国产精品人妻久久久久久| 国产乱人视频| 18禁在线无遮挡免费观看视频| 免费一级毛片在线播放高清视频| 精品人妻视频免费看| 亚洲高清免费不卡视频| 日韩欧美在线乱码| 亚洲精品色激情综合| 日韩一区二区三区影片| 国产午夜福利久久久久久| 一级毛片久久久久久久久女| 成人一区二区视频在线观看| 亚洲伊人久久精品综合 | 乱码一卡2卡4卡精品| 男女啪啪激烈高潮av片| 成人无遮挡网站| 麻豆一二三区av精品| 国内揄拍国产精品人妻在线| 精品欧美国产一区二区三| 色5月婷婷丁香| 久久久久久伊人网av| 免费在线观看成人毛片| 身体一侧抽搐| 一区二区三区免费毛片| 久久久久久久久中文| 婷婷六月久久综合丁香| 亚洲国产精品成人久久小说| 夫妻性生交免费视频一级片| 国产黄色小视频在线观看| 精品久久久噜噜| 日本色播在线视频| av女优亚洲男人天堂| 黄色配什么色好看| 国产视频首页在线观看| 夜夜爽夜夜爽视频| 精品久久久久久电影网 | 久久婷婷人人爽人人干人人爱| 最近的中文字幕免费完整| av天堂中文字幕网| 我要搜黄色片| 国产欧美日韩精品一区二区| 久久久成人免费电影| 麻豆精品久久久久久蜜桃| 亚洲成色77777| 国产精品一区二区三区四区久久| 嫩草影院入口| 亚洲精品乱码久久久v下载方式| 亚洲第一区二区三区不卡| 免费大片18禁| 国产探花在线观看一区二区| 免费播放大片免费观看视频在线观看 | 高清在线视频一区二区三区 | 乱系列少妇在线播放| 白带黄色成豆腐渣| 亚洲av电影在线观看一区二区三区 | 18禁裸乳无遮挡免费网站照片| 中国美白少妇内射xxxbb| 美女黄网站色视频| 特大巨黑吊av在线直播| 男人舔奶头视频| 最后的刺客免费高清国语| 91精品一卡2卡3卡4卡| 久久久久久久久大av| 国产精品一二三区在线看| 夜夜看夜夜爽夜夜摸| 偷拍熟女少妇极品色| 国产 一区 欧美 日韩| 亚洲精品,欧美精品| 亚洲最大成人手机在线| 人妻制服诱惑在线中文字幕| 国产视频首页在线观看| 又黄又爽又刺激的免费视频.| 欧美日韩精品成人综合77777| 亚洲成人久久爱视频| 国产精品一区二区在线观看99 | av在线天堂中文字幕| 国产精品麻豆人妻色哟哟久久 | 欧美一区二区亚洲| 看片在线看免费视频| 亚洲美女视频黄频| 晚上一个人看的免费电影| 男女啪啪激烈高潮av片| 1024手机看黄色片| 在线a可以看的网站| 日韩强制内射视频| 日韩视频在线欧美| 高清日韩中文字幕在线| 99久久人妻综合| 久久久午夜欧美精品| 亚洲最大成人手机在线| 欧美区成人在线视频| 久久综合国产亚洲精品| kizo精华| 高清av免费在线| 亚洲美女搞黄在线观看| 看非洲黑人一级黄片| 秋霞在线观看毛片| 亚洲av日韩在线播放| 嫩草影院入口| 岛国毛片在线播放| or卡值多少钱| 久久这里有精品视频免费| 美女黄网站色视频| av黄色大香蕉| 国产精品av视频在线免费观看| a级毛色黄片| 午夜福利网站1000一区二区三区| 免费观看的影片在线观看| 国产精品一区二区三区四区免费观看| 日韩成人av中文字幕在线观看| 啦啦啦韩国在线观看视频| 日韩制服骚丝袜av| 九九爱精品视频在线观看| 91av网一区二区| 99热全是精品| 久久久午夜欧美精品| 1024手机看黄色片| 永久网站在线| 午夜福利视频1000在线观看| 国产精品久久久久久久电影| 国语对白做爰xxxⅹ性视频网站| 色综合色国产| 日日摸夜夜添夜夜爱| 国产伦一二天堂av在线观看| 欧美性猛交黑人性爽| 欧美97在线视频| 亚洲av.av天堂| 欧美成人免费av一区二区三区| 中文字幕av成人在线电影| 春色校园在线视频观看| 美女脱内裤让男人舔精品视频| 一二三四中文在线观看免费高清| 九草在线视频观看| 99久久无色码亚洲精品果冻| 国产精品国产高清国产av| 18禁在线播放成人免费| 91狼人影院| 国产亚洲最大av| 成人午夜精彩视频在线观看| 欧美一区二区国产精品久久精品| 波多野结衣巨乳人妻| 九色成人免费人妻av| 又黄又爽又刺激的免费视频.| 草草在线视频免费看| 超碰97精品在线观看| 日韩人妻高清精品专区| 成人av在线播放网站| 在线天堂最新版资源| 欧美xxxx黑人xx丫x性爽| 亚洲欧美精品专区久久| a级毛色黄片| 男女啪啪激烈高潮av片| 三级毛片av免费| 人妻夜夜爽99麻豆av| 91久久精品国产一区二区成人| 97超碰精品成人国产| 久久精品综合一区二区三区| 看片在线看免费视频| 国产国拍精品亚洲av在线观看| 国产 一区 欧美 日韩| 国产精品av视频在线免费观看| 成人亚洲精品av一区二区| 哪个播放器可以免费观看大片| 国产麻豆成人av免费视频| 国产免费一级a男人的天堂| 日韩欧美精品免费久久| 精品国产露脸久久av麻豆 | 亚洲精品国产av成人精品| 欧美最新免费一区二区三区| 国产精品久久视频播放| 夜夜看夜夜爽夜夜摸| 久久久久九九精品影院| 99视频精品全部免费 在线| 亚洲高清免费不卡视频| 黄色日韩在线| 日韩欧美精品免费久久| 你懂的网址亚洲精品在线观看 | 2022亚洲国产成人精品| 中文字幕精品亚洲无线码一区| 日韩av不卡免费在线播放| 国产淫片久久久久久久久| 国产免费一级a男人的天堂| 日韩亚洲欧美综合| 男的添女的下面高潮视频| 国产精品美女特级片免费视频播放器| 精品人妻偷拍中文字幕| 亚洲欧美一区二区三区国产| 欧美一区二区精品小视频在线| 久久久国产成人免费| 大话2 男鬼变身卡| 国产成人aa在线观看| 中文字幕av成人在线电影| 51国产日韩欧美| 日韩av不卡免费在线播放| 亚洲成人久久爱视频| ponron亚洲| 午夜爱爱视频在线播放| 亚洲激情五月婷婷啪啪| a级毛色黄片| 国产一区二区三区av在线| 男人狂女人下面高潮的视频| 日韩中字成人| 成年av动漫网址| 久久精品久久久久久噜噜老黄 | 久久精品影院6| 中文字幕久久专区| 国产黄a三级三级三级人| 成人一区二区视频在线观看| a级毛片免费高清观看在线播放| 亚洲av中文字字幕乱码综合| 乱人视频在线观看| 人妻制服诱惑在线中文字幕| 国产精品久久视频播放| 你懂的网址亚洲精品在线观看 | 日韩欧美精品v在线| 亚洲精品456在线播放app| 男人舔女人下体高潮全视频| 午夜久久久久精精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久午夜福利片| 三级男女做爰猛烈吃奶摸视频| videossex国产| 人人妻人人澡人人爽人人夜夜 | 午夜精品在线福利| 久久久亚洲精品成人影院| 麻豆精品久久久久久蜜桃| 黄片无遮挡物在线观看| 极品教师在线视频| 精品久久久久久电影网 | 淫秽高清视频在线观看| 伦精品一区二区三区| 国产精品精品国产色婷婷| 亚洲在久久综合| 一个人观看的视频www高清免费观看| 亚洲色图av天堂| 欧美xxxx黑人xx丫x性爽| 国产精品久久久久久久电影| 亚洲国产精品成人久久小说| 欧美一区二区精品小视频在线| videossex国产| 亚洲在久久综合| 亚洲av电影在线观看一区二区三区 | 青春草国产在线视频| 自拍偷自拍亚洲精品老妇| 婷婷色av中文字幕| 亚洲怡红院男人天堂| 老女人水多毛片| 国产成年人精品一区二区| 免费人成在线观看视频色| 男女国产视频网站| 国产伦一二天堂av在线观看| 六月丁香七月| 国产高潮美女av| 99热全是精品| 3wmmmm亚洲av在线观看| 六月丁香七月| 日本免费在线观看一区| 午夜爱爱视频在线播放| 亚洲乱码一区二区免费版| 少妇人妻精品综合一区二区| 1000部很黄的大片| 国产亚洲精品久久久com| 搞女人的毛片| 亚洲无线观看免费| 在线观看一区二区三区| 中文乱码字字幕精品一区二区三区 | 色播亚洲综合网| 尾随美女入室| 国产精品三级大全| 又黄又爽又刺激的免费视频.| 少妇人妻一区二区三区视频| 国产成人aa在线观看| 欧美日韩一区二区视频在线观看视频在线 | 狂野欧美白嫩少妇大欣赏| 国产一级毛片七仙女欲春2| 国模一区二区三区四区视频| 欧美不卡视频在线免费观看| 免费观看人在逋| 色5月婷婷丁香| 激情 狠狠 欧美| 久久精品熟女亚洲av麻豆精品 | 97人妻精品一区二区三区麻豆| 色视频www国产| 成人国产麻豆网| 成人美女网站在线观看视频| 美女xxoo啪啪120秒动态图| 国产精品国产高清国产av| 日本一二三区视频观看| 汤姆久久久久久久影院中文字幕 | 日韩一区二区视频免费看| 青青草视频在线视频观看| 一个人看的www免费观看视频| 精品欧美国产一区二区三| 婷婷色麻豆天堂久久 | 亚洲内射少妇av| 日韩成人伦理影院| 美女内射精品一级片tv| 永久免费av网站大全| 女人被狂操c到高潮| 国产成人91sexporn| 成人一区二区视频在线观看| 最近中文字幕高清免费大全6| av福利片在线观看| 欧美xxxx性猛交bbbb| 久久精品久久久久久噜噜老黄 | 亚洲av一区综合| 亚洲成人av在线免费| 国产精品美女特级片免费视频播放器| 日韩中字成人| 欧美一级a爱片免费观看看| 欧美性感艳星| 国产熟女欧美一区二区| 日韩一区二区三区影片| 久久99热这里只有精品18| 免费人成在线观看视频色| 搡老妇女老女人老熟妇| 精品国产一区二区三区久久久樱花 | 九九热线精品视视频播放| 国产毛片a区久久久久| 午夜福利视频1000在线观看| 午夜日本视频在线| 亚洲精品乱码久久久v下载方式| 天堂√8在线中文| 亚洲精华国产精华液的使用体验| 丝袜美腿在线中文| 国产精品一区二区性色av| 麻豆av噜噜一区二区三区| 亚洲国产精品专区欧美| 波多野结衣巨乳人妻| 99久久无色码亚洲精品果冻| 日韩人妻高清精品专区| 免费黄网站久久成人精品| 亚洲国产欧美在线一区| 久久久久久久久久成人| 国语对白做爰xxxⅹ性视频网站| 男人舔奶头视频| 久久久久久久久久久丰满| 国产av不卡久久| 伊人久久精品亚洲午夜| 国产亚洲av嫩草精品影院| 亚洲真实伦在线观看| 欧美成人精品欧美一级黄| 日韩三级伦理在线观看| 国产精品99久久久久久久久| 国产在视频线精品| 久久99蜜桃精品久久| 天堂影院成人在线观看| 激情 狠狠 欧美| 在线观看66精品国产| 熟女人妻精品中文字幕| 欧美bdsm另类| 亚洲丝袜综合中文字幕| 卡戴珊不雅视频在线播放| 身体一侧抽搐| 亚洲欧美精品综合久久99| 国产精品一区二区三区四区久久| 免费看光身美女| 国产片特级美女逼逼视频| 少妇丰满av| 中文字幕久久专区| 两个人视频免费观看高清| 国内精品美女久久久久久| 特大巨黑吊av在线直播| 成年女人永久免费观看视频| 国产精品一区二区在线观看99 | 精品久久久久久成人av| 有码 亚洲区| 少妇高潮的动态图| 亚洲欧美一区二区三区国产| 亚洲伊人久久精品综合 | 日日撸夜夜添| 色5月婷婷丁香| 一本一本综合久久| www日本黄色视频网| 国产成人精品一,二区| 日韩成人av中文字幕在线观看| 日韩,欧美,国产一区二区三区 | 日韩在线高清观看一区二区三区| 精品国产露脸久久av麻豆 | 亚洲欧洲国产日韩| 欧美一级a爱片免费观看看| 美女xxoo啪啪120秒动态图| 国产精品一区二区在线观看99 | 国产亚洲精品久久久com| 波多野结衣高清无吗| 99热这里只有精品一区| 高清av免费在线| 亚洲性久久影院| 欧美xxxx黑人xx丫x性爽| 国产爱豆传媒在线观看| 美女高潮的动态| 99久久中文字幕三级久久日本| 少妇的逼水好多| 国产免费视频播放在线视频 | 草草在线视频免费看| 成人美女网站在线观看视频| 淫秽高清视频在线观看| 三级经典国产精品| 99在线视频只有这里精品首页| 91av网一区二区| 亚洲av成人精品一区久久| 日韩强制内射视频| 国产黄片美女视频| 国产美女午夜福利| 男人舔女人下体高潮全视频| 成人综合一区亚洲| 超碰av人人做人人爽久久| 美女脱内裤让男人舔精品视频| 欧美xxxx性猛交bbbb| 久久精品夜夜夜夜夜久久蜜豆| 乱码一卡2卡4卡精品| 久久精品影院6| 国产精品人妻久久久久久| 久久综合国产亚洲精品| 伊人久久精品亚洲午夜| 久久午夜福利片| 天天躁夜夜躁狠狠久久av| videos熟女内射| 超碰av人人做人人爽久久| 人妻夜夜爽99麻豆av|