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

    松原市淺層地下水污染源識別及風險研究

    2021-11-29 06:02:12開偉萌梁秀娟肖長來孟凡傲楊偉飛
    中國農(nóng)村水利水電 2021年11期
    關(guān)鍵詞:松原市淺層污染源

    開偉萌,梁秀娟,肖長來,孟凡傲,楊偉飛,王 鴿

    (1.吉林大學地下水資源與環(huán)境教育部重點實驗室,長春130021;2.吉林大學水資源與水環(huán)境重點實驗室,長春130021;3.吉林大學油頁巖地下原位轉(zhuǎn)化與鉆采技術(shù)國家地方聯(lián)合工程實驗室,長春130021;4.吉林大學新能源與環(huán)境學院,長春130021)

    0 引 言

    近年來,隨著城市化和工業(yè)化進程的加快,我國地下水污染問題日趨嚴重。其成因可簡要歸納為原生水質(zhì)問題和人類活動導致的水質(zhì)問題[1]。對地下水水質(zhì)進行檢測和評價,識別其主要污染源、空間分布及污染風險,在地下水治理和開發(fā)利用中具有重要意義。

    以往的水質(zhì)評價方法有單因子指數(shù)法、內(nèi)梅洛指數(shù)法、模糊數(shù)學評價法、灰色聚類法等[2,3]。這些方法可以對當前水質(zhì)進行評價與分級,但難以識別地下水污染源,且忽略了地下水污染的空間特性。前人為確定地下水污染源做出了大量的研究。唐軍等[4]將哈斯圖解法與綜合評價法結(jié)合用于地下水優(yōu)控污染物識別;馮茜等[5]利用因子分析法識別了德惠市淺層地下水的污染源,并分析了污染源的空間分布特征;郭濤等[6]利用正定矩陣因子分解法識別了拉林河流域地下水的主要污染源。曹陽等[7]利用因子分析法和系統(tǒng)聚類法對水質(zhì)影響因素進行評價。綜合來看,因子分析法,不但可以從宏觀上描述水質(zhì)因子,識別出污染源,還可以利用因子得分將研究區(qū)各污染因子及綜合污染因子的空間分布情況在圖中直觀表達,是識別地下水污染源的理想方法。

    本文以松原市為研究對象,選取2017年11月58 口潛水監(jiān)測井的水質(zhì)檢測數(shù)據(jù)進行分析研究,首次采用因子分析法確定松原市淺層地下水污染因子,分析污染源分布特征,并結(jié)合USEPA 的非致癌風險模型找出污染因子中的高風險因子,為研究區(qū)有關(guān)部門治理和開發(fā)利用地下水提供依據(jù)。

    1 研究區(qū)概況及數(shù)據(jù)來源

    1.1 研究區(qū)概況

    松原市位于松嫩平原以南,東與長春市、四平市接壤,西與白城市接壤,屬溫帶大陸性氣候。區(qū)內(nèi)東部為由王府-伏龍泉砂礫石臺地和松拉河間地塊組成的高平原,其第四系孔隙潛水含水層為下更新統(tǒng)砂礫石含水巖組(Q1)和中更新統(tǒng)沖積砂礫石含水巖組(Q2);西部為低平原,其地勢低平開闊,向二松河谷方向緩傾,其第四系孔隙潛水含水層為上更新統(tǒng)沖湖積粉細砂、砂礫石含水巖組(Q3);東部高平原與西部低平原之間為河谷平原,其第四系孔隙潛水含水層為全新統(tǒng)沖湖積、湖沼積淤泥質(zhì)亞黏土含水巖組(Q4)[8],如圖1(a)所示。區(qū)內(nèi)淺層地下水的補給來源主要為降水入滲補給,在西部低平原區(qū),由于徑流較為緩慢,淺層地下水的排泄方式主要以蒸發(fā)和人工開采為主。

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

    選用松原市2017年11月58組淺層地下水水質(zhì)檢測數(shù)據(jù)進行分析,取樣分布,如圖1(b)所示。檢測的潛水含水層厚度在7~40 m 之間,按均勻分布的原則,每個取樣點取三瓶水樣,取樣完后立即用錫紙密封送往譜尼測試(吉林分公司)化驗,檢測結(jié)果真實可靠。利用Piper 圖法確定研究區(qū)地下水化學類型時選取鉀、鈣、鈉、鎂、碳酸根、重碳酸根、硫酸根、氯離子8項指標,利用因子分析法識別淺層地下水污染源時選取pH、鋇、鐵、鉬、鈉、鋅、氟化物、氯化物、硝酸鹽、硫酸鹽、溶解性總固體(TDS)、總硬度、耗氧量、氨氮等18項指標。

    圖1 研究區(qū)地形地貌及采樣分布圖Fig.1 Topography and sampling distribution map

    2 研究方法與步驟

    2.1 研究方法

    本研究采用spss 軟件對地下水各離子指標進行統(tǒng)計性分析;運用Aquachem軟件繪制Piper圖,以確定淺層地下水水化學類型。采用因子分析法識別研究區(qū)地下水中的污染源;將因子分析提取的污染因子代入非致癌物風險模型,計算風險指數(shù)。因子分析數(shù)學模型與非致癌物風險模型如下。

    因子分析數(shù)學模型[1]:

    式中:f1,f2,…,fk為因子變量,a11,a12,…,apk為因子xp與變量fk之間的相關(guān)系數(shù),也稱為載荷;εi為特殊因子,代表公因子以外的影響因素,實際分析時可忽略不計。

    非致癌物風險模型[9]:

    式中:RI為非致癌物所致的個人平均風險,無量綱;DI為經(jīng)飲水途徑所導致的單位日均暴露劑量,mg/(kg·d);RiDI為經(jīng)飲水途徑的參考劑量,mg/(kg·d);CI為水中污染物的濃度,mg/L;IR為攝入率,一般取值2.1 L/d;EF為暴露頻率,一般取值365 d/a;ED為暴露持續(xù)時間,a;BW為平均體重,一般取值70 kg;AT平均壽命,取70 a[9]。

    2.2 具體步驟

    2.2.1 數(shù)據(jù)標準化

    將選取的各項指標的數(shù)據(jù)進行標準化處理,其公式如下。

    式中:xi為第i個水樣的原始測量值;x′i為第i個水樣去除量綱后的結(jié)果;s為其標準差。

    2.2.2 數(shù)據(jù)檢驗

    為檢驗變量之間是否有足夠的相關(guān)性,選用KMO 檢驗和Bartlett球型檢測法,作為檢驗方法。檢驗結(jié)果如表1所示。

    表1 KMO 和巴特利特檢驗Tab.1 KMO and Bartlett tests

    檢測結(jié)果顯示,KMO 檢驗結(jié)果為0.698,大于0.6,巴特利特檢驗結(jié)果為0,小于0.05,檢驗結(jié)果適合做因子分析。

    2.2.3 因子提取

    提取其中幾個能反映原始變量絕大部分信息的公因子,其特征值應(yīng)大于1,為使公因子具有良好的代表性,其累計方差值應(yīng)大于70%。

    2.2.4 因子旋轉(zhuǎn)

    采用最大方差法對矩陣進行旋轉(zhuǎn),旋轉(zhuǎn)后的載荷矩陣應(yīng)使各個因子的載荷方差之和最大,以便對評價內(nèi)容進行求解與分析。

    2.2.5 計算因子得分

    因子得分是指所提取的公共因子f1,f2,…,fk在每個取樣點的具體得分,每個因子得分由以下因子得分公式給出。

    以上步驟借助spss軟件中的因子分析模塊實現(xiàn)。

    2.2.6 因子風險研究

    采用USEPA 的非致癌物風險模型[9]對因子分析法所提取的污染因子進行風險指數(shù)計算,以EPA 推薦的最大可接受風險度(1.0×10-6)[10]為參照依據(jù),研究分析各污染因子的風險程度。

    3 結(jié)果與分析

    3.1 區(qū)域水化學特征

    通過分析研究區(qū)淺層地下水Piper 圖,如圖2 所示,可以看出研究區(qū)水化學類型主要以Ca-HCO3型水,CaNa-HCO3型水為主,少部分地區(qū)為Ca-Cl-HCO3型水。由表2(研究區(qū)主要水化學指標統(tǒng)計表)可知,該區(qū)淺層地下水的pH均值為7.64,呈弱堿性。TDS 和總硬度的變異系數(shù)較大,說明在全區(qū)范圍內(nèi)礦化度相差較大,應(yīng)存在地下水富集的高污染地區(qū),研究區(qū)Cl-、NO2-、NO3-、NH4+、Fe、As、Zn 等離子的標準差和變異系數(shù)較大,說明部分地區(qū)相關(guān)離子的含量偏高,可能是潛在的污染源。

    表2 研究區(qū)主要水化學指標統(tǒng)計表Tab.2 Statistical table of main hydrochemical indexes in the study area

    圖2 淺層地下水Piper圖Fig.2 Piper diagram of shallow groundwater

    3.2 地下水污染源識別及空間分布情況

    在對選取的18 項指標進行標準化處理,KMO 檢驗和Bartlett 球型檢測后,提取5 個公因子,累計方差貢獻率為75.546%,如表3 所示。為使各公因子的典型代表變量更加突出,將因子荷載矩陣進行正交旋轉(zhuǎn),旋轉(zhuǎn)因子荷載矩陣,見表4。

    表3 總方差解釋Tab.3 Total variance interpretation

    表4 旋轉(zhuǎn)因子載荷矩陣Tab.4 Rotation factor load matrix

    將各公因子得分及綜合因子得分,與對應(yīng)的坐標結(jié)合,運用ArcGIS 中的克里金插值法得到淺層地下水各公因子與綜合污染因子得分的空間分布圖,如圖3 所示。圖中某區(qū)域單項因子得分較高,說明該地區(qū)某項離子污染程度較為嚴重;綜合因子得分較高,說明該地區(qū)整體離子濃度偏高,可能存在多項離子污染超標的現(xiàn)象。

    f1因子的方差貢獻率為38.151%,是區(qū)域淺層地下水的主要成分,f1因子的主要代表變量有溶解性總固體(TDS)、氯化物、總硬度、硝酸鹽氮、鈉離子、硫酸鹽,將f1命名為“主成分”因子。f1因子的得分在區(qū)域內(nèi)呈西高東低,得分較高的地區(qū)為乾安縣及長嶺縣西部區(qū)域,如圖3(a)所示。f1因子的得分與水文地質(zhì)條件、水動力條件以及人類活動有關(guān)。松原市西部地區(qū)為低平原區(qū),干旱半干旱氣候,地勢低洼,地下水水位較低,周邊的地下水與巖石長期發(fā)生溶濾作用攜帶大量離子向此處運移,在長期蒸發(fā)濃縮作用下,地下水中的離子濃度不斷增高,導致該地區(qū)f1因子得分較高。其次,生活污水和工業(yè)污水中的酸和鹽類物質(zhì)排入地表水后,由地表水下滲到地下水中,與巖層發(fā)生陽離子交替作用,使巖層中的鈣、鎂離子析出,這使得地下水的總硬度逐漸增高[11]。

    f2因子的方差貢獻率為13.036%,主要代表變量有亞硝酸鹽氮和氨氮。f2因子得分自西向東逐漸降低,其中長嶺縣及乾安縣西部及寧江區(qū)因子得分較高,如圖3(b)所示。從空間分布來看,f2因子得分較高的地區(qū)主要分布在西部低平原區(qū)及人口密集的寧江區(qū),這是由于亞硝酸鹽氮和氨氮主要和農(nóng)藥化肥的使用及污水的排放有關(guān)。在低平原地區(qū),由于地勢平坦開闊,農(nóng)田較為集中,單位面積農(nóng)藥化肥的使用量遠高于東部地區(qū),過量的氮肥進入土壤后隨雨水進入地下水,造成地下水中氮含量逐漸增高。在寧江區(qū),工廠排放的工業(yè)廢水及集中排放的生活污水存在大量的氮化物,處理未達標的污水被排入第二松花江河道內(nèi)后,在遷移時發(fā)生下滲和越流補給現(xiàn)象,使得市區(qū)附近的淺層地下水中氮含量逐漸升高。

    f3因子的方差貢獻率為10.128%,主要代表變量有鉬、鐵、和有機物。f3因子主要分布在乾安縣西部地區(qū),如圖3(c)所示。參照《地下水質(zhì)量標準》(GB/T14848-2017),鉬和三氯甲烷的離子濃度均在III類標準以內(nèi),屬正常范圍。地下水中鐵元素的含量與TDS 的含量呈正相關(guān),與Ph 的值呈負相關(guān),同時還與地下水中的氧化還原環(huán)境有直接關(guān)系,當?shù)叵滤幱谶€原環(huán)境時鐵元素更容易從土壤中析出。在西部低平原區(qū),鐵元素主要來源于土壤及整個含水層[12]。

    f4因子的方差貢獻率為7.724%,主要代表變量是氟化物。f4因子主要分布于長嶺縣和乾安縣兩地,其中長嶺縣因子得分最高,如圖3(d)所示。氟元素含量主要受水文地質(zhì)條件和原生地質(zhì)環(huán)境影響。長嶺、乾安兩地,屬低平原區(qū),為干旱半干旱氣候,其地形低洼,沉積物以沖湖積、風積為主,含水層滲透性差,地下水交替緩慢,蒸發(fā)排泄為其主要排泄方式,易于離子的富集。其次,周邊裸露的氟礦物(方解石、石膏、長石類礦物和螢石等)經(jīng)過長期的風化、淋溶、并借助水動力條件富集于此。發(fā)源于內(nèi)蒙古大興安嶺德福勒罕山北麓的霍林河,是運輸大興安嶺含氟巖石的主要載體,為該地區(qū)源源不斷的提供氟離子[13]。

    圖3 各公因子與綜合污染因子得分空間分布圖Fig.3 Spatial distribution of the scores of common factors and comprehensive pollution factors

    f5因子的方差貢獻率為6.507%,主要代表變量是砷元素。f5因子分布于松原市絕大部分地區(qū),如圖3(e)所示。參照《地下水質(zhì)量標準》(GB/T14848-2017),砷的離子濃度在III 類標準以內(nèi),屬正常范圍。原生地質(zhì)作用是地下水中砷的主要來源,含砷巖石經(jīng)風化和雨水沖蝕等過程釋放可溶態(tài)砷吸附于土壤中,又在雨水作用下經(jīng)過長時間的下滲進入地下水中。另外,地下水中的鐵元素會和固態(tài)砷產(chǎn)生鐵砷循環(huán)耦合作用,地下水中鐵元素含量越高,越容易促進砷的釋放[14]。

    通過公式f綜合=(38.151f1+13.036f2+10.128f3+7.724f4+6.057f5)/75.546 計算松原市淺層地下水污染綜合得分f。結(jié)果顯示,長嶺縣西部地區(qū)的太平川鎮(zhèn)、北正鎮(zhèn),以及乾安縣西部地區(qū)永字鎮(zhèn)、大布蘇鎮(zhèn)污染最為嚴重,長嶺縣縣區(qū)、扶余縣縣區(qū)、及前郭爾羅斯蒙古族自治縣哈拉毛都鎮(zhèn)、長龍鄉(xiāng)、東三家子鄉(xiāng)水質(zhì)較好,如圖3(f)所示。綜上所述,綜合污染因子得分較高地區(qū)主要為低平原區(qū),其主要受原生地質(zhì)條件、水文地質(zhì)條件及地下水動力條件共同作用;其次,人類生產(chǎn)生活所排放的生活污水、工業(yè)污水及過度施用的農(nóng)藥也是造成該地淺層地下水污染的原因。

    3.3 因子風險研究

    在污染源識別的基礎(chǔ)上,選取氟、鐵、硝酸鹽氮、氨氮4個污染因子進行風險指數(shù)計算。依據(jù)美國環(huán)保署公眾健康評估手冊[15],硝酸鹽氮、氨氮、氟、鐵的參考劑量如表5所示。將指標含量及各參數(shù)代入非致癌風險模型計算,結(jié)果如圖4所示。

    圖4 各因子風險指數(shù)Fig.4 Risk index of each factor

    表5 非致癌物參考劑量 mg/(kg·d)Tab.5 Reference dose of non-carcinogens

    以EPA 推薦的最大可接受風險度(1.0×10-6)[12]為參照依據(jù),在58 個點的檢測樣本中,氟化物有13 個點的風險指數(shù)高于最大可接受風險度,超標率達22.41%,最大超標倍數(shù)為2.48倍;硝酸鹽氮有6 個點的風險指數(shù)高于最大可接受風險度,超標率達10.34%,最大超標倍數(shù)為3.40 倍;鐵元素有1 個超標,超標率為1.72%,氨氮無超標。結(jié)果表明地下水中氟化物、硝酸鹽氮的風險指數(shù)超標率偏高,鐵元素和氨氮的風險指數(shù)超標率較低。綜上所述,氟化物為風險最大的污染因子,其次是硝酸鹽氮類,鐵和氨氮因子的風險性較小。

    4 結(jié) 論

    (1)松原市淺層地下水總體上呈弱堿性,部分地區(qū)的礦化度較高。地下水中主要陽離子有Ca+、Na+、Mg2+,主要陰離子有HCO3-、Cl-。區(qū)內(nèi)地下水水化學類型主要為:Ca-HCO3、CaNa-HCO3型水。

    (2)影響松原市淺層地下水水質(zhì)的主要因子有TDS 和硝酸鹽氮,其方差貢獻率達38.151%;其次是亞硝酸鹽氮、氨氮和氟化物;其方差貢獻率總和達20.76%;鐵和砷的方差貢獻率之和在17%左右。

    (3)污染因子風險研究表明,區(qū)內(nèi)氟化物是風險最大的污染因子,其次是硝酸鹽氮類,鐵和氨氮因子的風險較小。

    (4)區(qū)內(nèi)淺層地下水污染較嚴重的區(qū)域主要集中在乾安縣與長嶺縣西部地區(qū),其主要受原生地質(zhì)環(huán)境作用和離子的遷移-富集作用影響。另外污水排放及化肥農(nóng)藥的過量施用也對淺層地下水水質(zhì)造成污染。 □

    猜你喜歡
    松原市淺層污染源
    淺層換填技術(shù)在深厚軟土路基中的應(yīng)用
    持續(xù)推進固定污染源排污許可管理全覆蓋
    基于淺層曝氣原理的好氧顆粒污泥的快速培養(yǎng)
    基于污染源解析的空氣污染治理對策研究
    十二五”期間佳木斯市污染源排放狀況分析
    看不見的污染源——臭氧
    淺層地下水超采區(qū)劃分探究
    松原市備春耕堅持高起點高標準
    向大作家學寫作秘訣
    ——讀《曹文軒教你寫作》有感
    小學時代(2016年17期)2016-02-23 06:45:59
    松原市林下經(jīng)濟發(fā)展模式和措施探討
    久久久久精品性色| 国产乱人偷精品视频| 精品视频人人做人人爽| 永久网站在线| 久久这里只有精品19| 一级毛片电影观看| 美女国产高潮福利片在线看| 日日啪夜夜爽| 免费高清在线观看日韩| 国产国语露脸激情在线看| 精品一区在线观看国产| 亚洲国产精品成人久久小说| 国产爽快片一区二区三区| 18+在线观看网站| 满18在线观看网站| 自拍欧美九色日韩亚洲蝌蚪91| 好男人视频免费观看在线| 日韩一区二区三区影片| 亚洲人成网站在线观看播放| 日本-黄色视频高清免费观看| 一区二区av电影网| 国产免费一区二区三区四区乱码| 久久精品国产亚洲av天美| 1024香蕉在线观看| 男女啪啪激烈高潮av片| 欧美精品一区二区免费开放| 国产97色在线日韩免费| 啦啦啦中文免费视频观看日本| 久久久久国产精品人妻一区二区| 伊人久久大香线蕉亚洲五| 少妇精品久久久久久久| 岛国毛片在线播放| 久久精品人人爽人人爽视色| 久久精品亚洲av国产电影网| www.熟女人妻精品国产| 亚洲国产精品成人久久小说| 热99久久久久精品小说推荐| 丝袜人妻中文字幕| 国产男女内射视频| 国产精品99久久99久久久不卡 | 色吧在线观看| 国产免费又黄又爽又色| 精品一品国产午夜福利视频| 成年人免费黄色播放视频| 午夜福利在线观看免费完整高清在| 欧美亚洲 丝袜 人妻 在线| 亚洲天堂av无毛| 看免费av毛片| 在线观看国产h片| 青青草视频在线视频观看| 一本大道久久a久久精品| 少妇人妻精品综合一区二区| 九色亚洲精品在线播放| 亚洲av综合色区一区| 亚洲精品中文字幕在线视频| 亚洲五月色婷婷综合| 成人二区视频| 丝袜在线中文字幕| 一个人免费看片子| 一个人免费看片子| 免费看不卡的av| 尾随美女入室| 成人免费观看视频高清| 亚洲av电影在线进入| 日本免费在线观看一区| 亚洲国产精品一区二区三区在线| 两性夫妻黄色片| 久久ye,这里只有精品| 人妻少妇偷人精品九色| 亚洲国产看品久久| 亚洲av中文av极速乱| 欧美日韩亚洲国产一区二区在线观看 | 日韩 亚洲 欧美在线| 日韩视频在线欧美| 久久人妻熟女aⅴ| 国产97色在线日韩免费| 水蜜桃什么品种好| 欧美激情 高清一区二区三区| 亚洲人成电影观看| 久久热在线av| 久久久久久久大尺度免费视频| 精品少妇黑人巨大在线播放| 一区二区av电影网| 成年人午夜在线观看视频| 国产免费又黄又爽又色| 亚洲人成电影观看| 国产成人a∨麻豆精品| 精品99又大又爽又粗少妇毛片| 十八禁网站网址无遮挡| 亚洲熟女精品中文字幕| 色播在线永久视频| 国产淫语在线视频| 亚洲精品,欧美精品| 久久午夜综合久久蜜桃| 三上悠亚av全集在线观看| 国产一区二区三区综合在线观看| 日韩制服骚丝袜av| 亚洲四区av| 国产色婷婷99| 国产精品国产三级专区第一集| av天堂久久9| 麻豆乱淫一区二区| 女性生殖器流出的白浆| 可以免费在线观看a视频的电影网站 | 婷婷色综合www| 久久久国产欧美日韩av| 婷婷色麻豆天堂久久| 国产精品女同一区二区软件| 亚洲国产成人一精品久久久| av.在线天堂| 可以免费在线观看a视频的电影网站 | 一区二区三区激情视频| 亚洲国产日韩一区二区| 国产97色在线日韩免费| 最近的中文字幕免费完整| 考比视频在线观看| 免费播放大片免费观看视频在线观看| 我要看黄色一级片免费的| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲av中文av极速乱| 成人漫画全彩无遮挡| 亚洲精品在线美女| 在线天堂最新版资源| 精品久久蜜臀av无| 亚洲三区欧美一区| 亚洲av电影在线进入| 国产精品亚洲av一区麻豆 | 女人高潮潮喷娇喘18禁视频| 老司机影院毛片| 母亲3免费完整高清在线观看 | 精品亚洲乱码少妇综合久久| 国产成人精品久久久久久| 777米奇影视久久| 国产老妇伦熟女老妇高清| 国产精品香港三级国产av潘金莲 | 欧美激情高清一区二区三区 | 国产精品久久久久久精品电影小说| 欧美97在线视频| a级毛片在线看网站| 国产一区二区 视频在线| 欧美亚洲日本最大视频资源| 日韩熟女老妇一区二区性免费视频| 久久久久精品久久久久真实原创| 黄色视频在线播放观看不卡| 国产免费现黄频在线看| 国产精品三级大全| 丝瓜视频免费看黄片| 亚洲人成77777在线视频| 中国三级夫妇交换| 天堂俺去俺来也www色官网| 久久久久精品人妻al黑| 性少妇av在线| 人妻少妇偷人精品九色| 看免费成人av毛片| 精品午夜福利在线看| 黄色视频在线播放观看不卡| 在线精品无人区一区二区三| 成人国产麻豆网| 午夜福利影视在线免费观看| 一区福利在线观看| 在线观看免费日韩欧美大片| 日韩免费高清中文字幕av| 男人添女人高潮全过程视频| 国产精品一国产av| 免费在线观看视频国产中文字幕亚洲 | 日韩精品有码人妻一区| 久久久久精品性色| 在线精品无人区一区二区三| 免费播放大片免费观看视频在线观看| 亚洲欧美精品综合一区二区三区 | 18禁裸乳无遮挡动漫免费视频| 女性被躁到高潮视频| 国产成人午夜福利电影在线观看| 69精品国产乱码久久久| 男人爽女人下面视频在线观看| 80岁老熟妇乱子伦牲交| 波多野结衣一区麻豆| 午夜日韩欧美国产| 久热这里只有精品99| 日本av免费视频播放| 亚洲国产欧美网| 国产xxxxx性猛交| 美女大奶头黄色视频| 亚洲情色 制服丝袜| 免费黄网站久久成人精品| 大话2 男鬼变身卡| 久久热在线av| 国产日韩欧美亚洲二区| freevideosex欧美| 精品卡一卡二卡四卡免费| 欧美97在线视频| 香蕉国产在线看| 大片免费播放器 马上看| 日本欧美视频一区| 久久99精品国语久久久| 热re99久久国产66热| 亚洲国产毛片av蜜桃av| 亚洲天堂av无毛| 中文字幕精品免费在线观看视频| 在线精品无人区一区二区三| 国产精品免费大片| 日韩制服丝袜自拍偷拍| 亚洲av国产av综合av卡| 一区二区三区乱码不卡18| 午夜免费观看性视频| 一区在线观看完整版| 亚洲国产最新在线播放| 亚洲男人天堂网一区| 极品人妻少妇av视频| av女优亚洲男人天堂| 中文字幕人妻丝袜一区二区 | 人妻人人澡人人爽人人| 亚洲国产精品999| 五月伊人婷婷丁香| 国产日韩一区二区三区精品不卡| 国产高清国产精品国产三级| 只有这里有精品99| 中文字幕色久视频| 超碰成人久久| 国产在线一区二区三区精| 成人亚洲精品一区在线观看| 婷婷色综合大香蕉| 在线观看美女被高潮喷水网站| 2021少妇久久久久久久久久久| 国产不卡av网站在线观看| 日韩精品有码人妻一区| 国产亚洲一区二区精品| 国产精品av久久久久免费| 亚洲欧洲日产国产| 精品亚洲乱码少妇综合久久| 如日韩欧美国产精品一区二区三区| 人妻人人澡人人爽人人| 日本黄色日本黄色录像| 丰满乱子伦码专区| tube8黄色片| 如日韩欧美国产精品一区二区三区| 久久国内精品自在自线图片| 亚洲国产日韩一区二区| 制服人妻中文乱码| 精品一区二区三卡| 精品国产乱码久久久久久小说| 97人妻天天添夜夜摸| 精品酒店卫生间| 日本vs欧美在线观看视频| 国产精品久久久久久av不卡| 精品午夜福利在线看| 街头女战士在线观看网站| 制服人妻中文乱码| 国产无遮挡羞羞视频在线观看| 欧美激情 高清一区二区三区| 欧美精品一区二区免费开放| 久久午夜综合久久蜜桃| 97精品久久久久久久久久精品| 国产亚洲精品第一综合不卡| 看十八女毛片水多多多| 国产欧美日韩综合在线一区二区| 我的亚洲天堂| 亚洲精品在线美女| 人妻一区二区av| 亚洲av国产av综合av卡| 中文字幕最新亚洲高清| 国产 精品1| 欧美97在线视频| 狠狠婷婷综合久久久久久88av| 最近最新中文字幕免费大全7| 国产1区2区3区精品| 亚洲av在线观看美女高潮| 黄片无遮挡物在线观看| 最近2019中文字幕mv第一页| 免费高清在线观看日韩| 看免费av毛片| 亚洲久久久国产精品| 亚洲综合色网址| 精品99又大又爽又粗少妇毛片| 亚洲人成电影观看| 一二三四在线观看免费中文在| 免费高清在线观看日韩| 午夜福利视频精品| 亚洲精品美女久久久久99蜜臀 | 欧美精品人与动牲交sv欧美| 一区二区三区乱码不卡18| 日本黄色日本黄色录像| 在线观看美女被高潮喷水网站| 久久精品aⅴ一区二区三区四区 | 欧美 日韩 精品 国产| 日韩av在线免费看完整版不卡| 亚洲欧美清纯卡通| 国产成人精品在线电影| av国产精品久久久久影院| 一区二区日韩欧美中文字幕| xxx大片免费视频| 久久久久久人妻| 黄片无遮挡物在线观看| 韩国高清视频一区二区三区| 1024视频免费在线观看| 国产精品久久久久久精品电影小说| 国产高清不卡午夜福利| 9191精品国产免费久久| 女人久久www免费人成看片| 欧美国产精品一级二级三级| 午夜福利乱码中文字幕| 亚洲综合色惰| 中文字幕亚洲精品专区| 中文欧美无线码| 最近的中文字幕免费完整| 国产成人精品久久久久久| av国产久精品久网站免费入址| 亚洲成av片中文字幕在线观看 | 久久久久久久亚洲中文字幕| 国产精品秋霞免费鲁丝片| 性色av一级| av网站免费在线观看视频| 一级毛片我不卡| 女性生殖器流出的白浆| 黄色怎么调成土黄色| 男女无遮挡免费网站观看| 男人添女人高潮全过程视频| 久久久精品国产亚洲av高清涩受| 午夜精品国产一区二区电影| av免费观看日本| 在线免费观看不下载黄p国产| 日本-黄色视频高清免费观看| av国产久精品久网站免费入址| 99国产精品免费福利视频| 久久久精品94久久精品| 男人添女人高潮全过程视频| 欧美日韩亚洲国产一区二区在线观看 | 成年女人毛片免费观看观看9 | 一级毛片我不卡| 伦理电影大哥的女人| 日本av手机在线免费观看| 亚洲第一青青草原| 国产精品一区二区在线不卡| 少妇的逼水好多| 国产一区亚洲一区在线观看| 波多野结衣av一区二区av| 日本vs欧美在线观看视频| 欧美成人精品欧美一级黄| 亚洲精品久久久久久婷婷小说| 水蜜桃什么品种好| 男女边摸边吃奶| 老司机亚洲免费影院| 中文字幕人妻丝袜制服| 男男h啪啪无遮挡| 国产一区二区激情短视频 | 制服丝袜香蕉在线| 在线 av 中文字幕| 欧美精品人与动牲交sv欧美| 免费黄色在线免费观看| 天天躁夜夜躁狠狠久久av| 国产在视频线精品| 国产av精品麻豆| 91精品国产国语对白视频| 十八禁网站网址无遮挡| 日日摸夜夜添夜夜爱| 亚洲,欧美精品.| 搡老乐熟女国产| 秋霞伦理黄片| 成人影院久久| 亚洲久久久国产精品| 久久综合国产亚洲精品| 亚洲男人天堂网一区| 丰满迷人的少妇在线观看| 日韩精品免费视频一区二区三区| 欧美国产精品一级二级三级| 1024香蕉在线观看| 在线观看美女被高潮喷水网站| kizo精华| 国产精品女同一区二区软件| 成年人免费黄色播放视频| 啦啦啦中文免费视频观看日本| 伦精品一区二区三区| 精品一区在线观看国产| 热99久久久久精品小说推荐| 久久久久久久久免费视频了| 在线天堂中文资源库| 校园人妻丝袜中文字幕| 亚洲国产日韩一区二区| 午夜日韩欧美国产| av.在线天堂| av不卡在线播放| 老司机影院毛片| 午夜老司机福利剧场| 丝袜美腿诱惑在线| 亚洲综合精品二区| 免费观看在线日韩| 曰老女人黄片| 亚洲精品aⅴ在线观看| 亚洲av免费高清在线观看| 中文字幕另类日韩欧美亚洲嫩草| av不卡在线播放| av.在线天堂| 国产精品嫩草影院av在线观看| 日本免费在线观看一区| 下体分泌物呈黄色| 久久99精品国语久久久| 国产精品蜜桃在线观看| 欧美日韩视频精品一区| 成人二区视频| 最黄视频免费看| 国产欧美日韩一区二区三区在线| 性少妇av在线| 国产免费又黄又爽又色| 精品一区在线观看国产| 国产毛片在线视频| 卡戴珊不雅视频在线播放| 成年人免费黄色播放视频| 久久国产亚洲av麻豆专区| 伦理电影大哥的女人| 亚洲成人手机| 日本免费在线观看一区| 人妻少妇偷人精品九色| 91国产中文字幕| 激情视频va一区二区三区| 欧美精品一区二区免费开放| 自线自在国产av| 久久久久久伊人网av| 亚洲欧美色中文字幕在线| 一本—道久久a久久精品蜜桃钙片| 日本欧美视频一区| 亚洲欧洲精品一区二区精品久久久 | 男人舔女人的私密视频| 亚洲国产av新网站| 丝袜喷水一区| 国产成人91sexporn| 黑人巨大精品欧美一区二区蜜桃| 极品人妻少妇av视频| 天天躁狠狠躁夜夜躁狠狠躁| 97在线视频观看| 欧美黄色片欧美黄色片| 国产成人a∨麻豆精品| 丝袜美腿诱惑在线| 国产 精品1| 一区二区三区精品91| 成人二区视频| 亚洲综合色惰| 一区福利在线观看| 天堂8中文在线网| 久久亚洲国产成人精品v| 国产精品女同一区二区软件| 色94色欧美一区二区| videossex国产| 亚洲成av片中文字幕在线观看 | 精品视频人人做人人爽| 久久久久久久久免费视频了| 看免费成人av毛片| 我要看黄色一级片免费的| 美女大奶头黄色视频| 国产亚洲一区二区精品| 日本wwww免费看| 波多野结衣av一区二区av| 欧美中文综合在线视频| 欧美黄色片欧美黄色片| 王馨瑶露胸无遮挡在线观看| 欧美精品一区二区大全| 国产亚洲一区二区精品| 极品人妻少妇av视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲av免费高清在线观看| 亚洲精品自拍成人| 久久久久久久国产电影| 最近最新中文字幕免费大全7| 日韩一区二区视频免费看| 人人澡人人妻人| www.av在线官网国产| 精品国产国语对白av| 免费女性裸体啪啪无遮挡网站| 水蜜桃什么品种好| 美女主播在线视频| 国产精品久久久av美女十八| 亚洲av电影在线观看一区二区三区| 国产成人精品福利久久| 色婷婷久久久亚洲欧美| 两性夫妻黄色片| 久久精品国产鲁丝片午夜精品| 精品久久蜜臀av无| 欧美黄色片欧美黄色片| 咕卡用的链子| 亚洲精品自拍成人| 久久久久久人妻| 亚洲综合色网址| 亚洲国产精品999| 99久久综合免费| 人成视频在线观看免费观看| 美女视频免费永久观看网站| 亚洲av男天堂| 1024香蕉在线观看| 久久午夜福利片| 母亲3免费完整高清在线观看 | 日产精品乱码卡一卡2卡三| 国产爽快片一区二区三区| 18禁裸乳无遮挡动漫免费视频| 国产精品国产av在线观看| 国产1区2区3区精品| 寂寞人妻少妇视频99o| 少妇人妻精品综合一区二区| 99热国产这里只有精品6| 欧美日韩一区二区视频在线观看视频在线| 日韩一区二区三区影片| 性少妇av在线| 国产不卡av网站在线观看| av卡一久久| 亚洲国产精品一区二区三区在线| 国产精品久久久久久精品古装| 一级毛片 在线播放| 日韩电影二区| 黄片播放在线免费| 最近中文字幕2019免费版| 国产 精品1| 欧美日韩一级在线毛片| 91成人精品电影| 国产精品一区二区在线不卡| 欧美日韩综合久久久久久| a 毛片基地| 国产精品成人在线| 美女视频免费永久观看网站| 精品午夜福利在线看| 亚洲久久久国产精品| 又黄又粗又硬又大视频| 欧美精品一区二区免费开放| 久久精品国产a三级三级三级| 日韩三级伦理在线观看| 久久毛片免费看一区二区三区| 91精品伊人久久大香线蕉| 亚洲国产精品一区二区三区在线| 美女视频免费永久观看网站| 777米奇影视久久| 亚洲精品美女久久久久99蜜臀 | 久久午夜福利片| 亚洲激情五月婷婷啪啪| 伊人久久国产一区二区| 在线观看国产h片| 99国产综合亚洲精品| 久久久久国产精品人妻一区二区| 亚洲国产av新网站| 亚洲精华国产精华液的使用体验| 午夜老司机福利剧场| 亚洲人成电影观看| 在线观看人妻少妇| 我的亚洲天堂| 亚洲av日韩在线播放| 亚洲欧美日韩另类电影网站| 日本色播在线视频| 国产亚洲av片在线观看秒播厂| 一级片'在线观看视频| 亚洲精品美女久久久久99蜜臀 | 亚洲国产最新在线播放| 超色免费av| 伊人久久国产一区二区| 两个人看的免费小视频| 午夜福利在线观看免费完整高清在| 国产精品久久久久久精品古装| xxx大片免费视频| 亚洲国产日韩一区二区| 欧美日韩av久久| 久久久精品国产亚洲av高清涩受| 免费看不卡的av| 人妻人人澡人人爽人人| 国产一区二区三区av在线| 国产亚洲精品第一综合不卡| 综合色丁香网| 国产女主播在线喷水免费视频网站| 国产欧美日韩一区二区三区在线| 国产精品欧美亚洲77777| 成年人免费黄色播放视频| 9热在线视频观看99| 国产精品秋霞免费鲁丝片| 宅男免费午夜| 亚洲四区av| 日韩伦理黄色片| 午夜影院在线不卡| 中文字幕av电影在线播放| a级毛片在线看网站| 十八禁高潮呻吟视频| 9色porny在线观看| 亚洲av日韩在线播放| 精品99又大又爽又粗少妇毛片| av一本久久久久| 国产熟女欧美一区二区| 性色av一级| 国产亚洲欧美精品永久| 国产成人精品在线电影| 韩国高清视频一区二区三区| 中文精品一卡2卡3卡4更新| 午夜福利一区二区在线看| 日本爱情动作片www.在线观看| 大话2 男鬼变身卡| www.自偷自拍.com| 丰满迷人的少妇在线观看| 亚洲av日韩在线播放| 男人舔女人的私密视频| 男女边吃奶边做爰视频| 99国产综合亚洲精品| 999精品在线视频| 国产精品三级大全| 叶爱在线成人免费视频播放| 狂野欧美激情性bbbbbb| 性少妇av在线| 久久久久精品人妻al黑| 国产黄色视频一区二区在线观看| 成人国产av品久久久| 最近中文字幕2019免费版| 超碰97精品在线观看| www.精华液| 国产精品av久久久久免费| 久热久热在线精品观看| 有码 亚洲区| 丝袜美腿诱惑在线| 女的被弄到高潮叫床怎么办| 国产亚洲最大av| 丝袜美腿诱惑在线| 久久久久久人人人人人| 亚洲欧洲日产国产| 在线天堂最新版资源| 精品99又大又爽又粗少妇毛片|