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

    地表水重金屬評價(jià)方法的改進(jìn)與運(yùn)用比較研究

    2021-07-19 09:44:38趙一蔚高良敏陳曉晴查甫更石闖張海強(qiáng)
    關(guān)鍵詞:污染方法研究

    趙一蔚 高良敏 陳曉晴 查甫更 石闖 張海強(qiáng)

    摘 要:針對重金屬污染指數(shù)(HPI)、重金屬評估指數(shù)(HEI)方法的部分局限性,對重金屬污染指數(shù)(HPI)進(jìn)行了改進(jìn),形成了改進(jìn)的地表水重金屬污染指數(shù)|r-HPI|方法,并分別應(yīng)用重金屬污染指數(shù)(HPI)、重金屬評估指數(shù)(HEI)、改進(jìn)的重金屬污染指數(shù)(|r-HPI|)對研究區(qū)域連續(xù)三年共1 164個(gè)地表水樣品統(tǒng)計(jì)分析,驗(yàn)證了|r-HPI|方法使用方便,具有靈活性,并且綜合考慮了不同重金屬之間的內(nèi)在聯(lián)系,有利于反映出水體重金屬污染的真實(shí)狀況。基于普通克里金插值的地理空間圖清楚地反映了|r-HPI|的年際變化。通過Kolmogorov-Smirnov檢驗(yàn)數(shù)據(jù)集的正態(tài)性,觀察到顯著差異(P <0.05)。|r-HPI|通過Bland-Altman分析并驗(yàn)證了其有效性,表明在地表水重金屬污染研究中的良好適用性。

    關(guān)鍵詞:地表水;重金屬評價(jià);修正重金屬污染指數(shù)(|r-HPI|);Bland-Altman分析

    中圖分類號: X824文獻(xiàn)標(biāo)志碼:A

    文章編號:1672-1098(2021)02-0045-06

    收稿日期:2020-11-08

    基金項(xiàng)目:安徽省重點(diǎn)研究與開發(fā)計(jì)劃基金資助項(xiàng)目(202004i07020012)

    作者簡介:趙一蔚(1998-),女,山東菏澤人,在讀碩士,研究方向:生態(tài)環(huán)境效應(yīng)、水污染控制、環(huán)境規(guī)劃管理與評價(jià)。

    Comparative Study on Improvement and Application of Evaluation Methods of Heavy Metals in Surface Water

    ZHAO Yiwei, GAO Liangmin, CHEN Xiaoqing, ZHA Fugeng SHI Chuang, ZHANG Haiqiang

    (School of Earth and Environment, Anhui University of Science and Technology, HuainanAnhui 232001,China)

    Abstract:To some of the limitations of the Heavy Metal Pollution Index (HPI) and Heavy Metal Evaluation Index (HEI) methods, HPI was specially improved in this article to form an improved surface water heavy metal pollution index |r-HPI| method. 1164 surface water samples in the study area for three consecutive years were analyzed statistically with HPI, HEI and |r-HPI|, verifying the convenience and flexibility of |r- HPI |method, which was in favor of reflecting the real situation of heavy metal pollution in water because of the internal relationship between different heavy metals comprehensively considered. The geospatial map based on ordinary Kriging interpolation reflected clearly the inter-annual changes of |r-HPI|. When the normality of the data set was tested by Kolmogorov-Smirnov, a significant difference was observed (P<0.05). |r-HPI| has been verified effective by Bland-Altman, showing its good applicability in the study of heavy metal pollution in surface water.

    Key words:surface water; Heavy metal assessment; modified heavy metal pollution index (|r-HPI|); Bland-Altman analysis

    在世界范圍內(nèi),水體重金屬的不可逆轉(zhuǎn)性和長期累積性對人類健康和生態(tài)安全的損害引起了廣泛關(guān)注[1-2]。重金屬污染指數(shù)[3](heavy metal pollution index,HPI)和重金屬評估指數(shù)[4-9](heavy metal evaluation index,HEI)被廣泛采用[10-13],但是也有不少弊端:HPI方法中的重要參數(shù)Qi與重金屬濃度、最高期望濃度的絕對差值有關(guān)。因此,使用HPI方法表示水體重金屬污染狀況存在一定誤差。HEI方法主要取決于允許濃度,未將水中重金屬濃度的標(biāo)準(zhǔn)限值考慮在內(nèi)。另外,無論是HPI還是HEI均是將不同種類的重金屬單獨(dú)劃分為一類,通過對不同類的指數(shù)求和獲得整體水體重金屬污染狀況,忽略了不同重金屬的潛在聯(lián)系,很有可能與實(shí)際情況產(chǎn)生偏差。

    針對現(xiàn)有HPI、HEI方法使用的局限性問題,本研究特提出通過使用數(shù)理統(tǒng)計(jì)方法修正HPI,旨在:①基于現(xiàn)有的重金屬污染評價(jià)指數(shù)評估淮河流域的水體重金屬,并改進(jìn)重金屬污染指數(shù)(|r-HPI|),使用多元統(tǒng)計(jì)技術(shù)并將其與GIS集成以進(jìn)行水體重金屬的空間表示;②通過Bland-Altman分析對HPI、HEI、|r-HPI|進(jìn)行有效性驗(yàn)證。

    1 資料與方法

    1.1 研究區(qū)域

    本研究選擇中國的淮河流域作為研究區(qū),該區(qū)域人口密集,地理位置特殊。研究區(qū)位于半濕潤與半干旱區(qū)域的交界地帶,是黃河流域和長江流域的過渡區(qū),地形主要以平原為主,西部及西南部、東北部是山區(qū)和丘陵,在行政劃分上,淮河流域橫跨河南、安徽、江蘇等省。

    1.2 樣品收集

    在收集研究區(qū)點(diǎn)位數(shù)據(jù)之前查閱了大量該地區(qū)的文獻(xiàn)資料,并進(jìn)行了實(shí)地勘察,本課題組通過水文年鑒收集到研究區(qū)的各個(gè)點(diǎn)位水質(zhì)數(shù)據(jù),并前往水文局及生態(tài)環(huán)境局,對已掌握的數(shù)據(jù)(站點(diǎn)實(shí)際位置等)進(jìn)行核實(shí)以及對未掌握或掌握不全的數(shù)據(jù)進(jìn)行溝通獲取。在此基礎(chǔ)上,本研究中數(shù)據(jù)點(diǎn)的選擇綜合考慮了淮河水系的一、二級支流的交界處,水系的較大支流匯入前的河口處,主要河流、水庫(響洪甸水庫、佛子嶺水庫、梅山水庫)的出入口,有水工建筑物(蚌埠閘、王家壩等)并受人工控制的河段。本研究共收集了2015~2017a連續(xù)三年36個(gè)數(shù)據(jù)點(diǎn)逐月水質(zhì)數(shù)據(jù),對數(shù)據(jù)初步整理后,共獲得1 164個(gè)可利用水體重金屬數(shù)據(jù)。

    1.3 水體重金屬污染評價(jià)方法

    重金屬污染指數(shù)HPI針對單個(gè)重金屬對水質(zhì)的相對重要性劃分權(quán)重,計(jì)算公式為

    ∑ni=1WiQi∑ni-1Wi(1)

    式中:Wi為單位權(quán)重,Qi是第i個(gè)參數(shù)的子指數(shù),在該研究中,選取6個(gè)重金屬指標(biāo)(Pb、Cu、As、Zn、Cd、Cr)。

    Qi=∑ni=1|Mi-Ii|Si-Ii×100(2)

    式中:Si代表第i個(gè)參數(shù)的標(biāo)準(zhǔn)值[14],Mi是監(jiān)測到的相應(yīng)參數(shù)的濃度。Ii是第i個(gè)參數(shù)的最高期望濃度[15],一般認(rèn)為低于50,表示污染較低;介于50~100,表示處于中度污染;大于100,表示污染程度高[16]。

    重金屬評估指數(shù)HEI是通過考慮每個(gè)指標(biāo)的監(jiān)測值與標(biāo)準(zhǔn)值之比的總和,一般認(rèn)為低于40,表示污染較低;介于40~80,表示處于中度污染;大于80,表示污染程度高[17]。計(jì)算公式為

    ∑ni=1MiSi(3)

    針對現(xiàn)有HPI、 HEI方法的使用局限性問題,本研究特提出通過使用數(shù)理統(tǒng)計(jì)方法修正HPI—|r-HPI|, |r-HPI|方法使用方便,具有靈活性,綜合考慮了不同重金屬之間的內(nèi)在聯(lián)系,有效反映污染的真實(shí)情況。

    xi,j=∑ni=1Mi-Ii|Si-Ii|×100(4)

    li,j=xi,jzi,j(5)

    pi=li,j×Zj(6)

    r-HPI=∑ni=1ki×pi(7)

    式中:xi,j為第j個(gè)重金屬指標(biāo)的第i個(gè)成分矩陣重因子載荷分布,zi,j為對應(yīng)的第j個(gè)重金屬指標(biāo)的第i個(gè)成份特征值。li,j是對應(yīng)的因子載荷權(quán)重,Zj為第j個(gè)重金屬指標(biāo)的標(biāo)準(zhǔn)化子指數(shù),pi為第j個(gè)重金屬參數(shù)的分配率,ki為相應(yīng)的方差貢獻(xiàn)率。

    首先對原始數(shù)據(jù)(xi,j)進(jìn)行標(biāo)準(zhǔn)化映射,然后將因子分析應(yīng)用于每個(gè)獨(dú)立簇的數(shù)據(jù)集,獲得了相應(yīng)子指數(shù)的因子載荷。使用(Keiser-Meyer-Olkin,KMO)準(zhǔn)則驗(yàn)證了數(shù)據(jù)的充分性, 與此同時(shí), Bartlett的球形度檢驗(yàn)表明了該數(shù)據(jù)集對因子分析的適用性。 將經(jīng)Varimax旋轉(zhuǎn)后成分矩陣系數(shù)與特征值進(jìn)行相除獲得對應(yīng)參數(shù)的因子載荷權(quán)重li,j,可以實(shí)現(xiàn)數(shù)據(jù)不相關(guān)化且正交。值得注意的是,累計(jì)方差百分比要高于80%。在獲得子得分后與對應(yīng)的貢獻(xiàn)率相乘獲得每個(gè)樣本的修正重金屬污染指數(shù)?;谏鲜霾襟E,本研究根據(jù)獲得的 |r-HPI|值將水樣劃分成3類:|r-HPI|<0.5時(shí), 重金屬污染程度低;0.5≤|r-HPI|<1時(shí),屬于中度污染;|r-HPI|≥1時(shí),屬于高度污染。

    2 結(jié)果與分析

    2.1 |r-HPI|

    在該項(xiàng)研究中, 使用6個(gè)重金屬參數(shù)計(jì)算了研

    究區(qū)36個(gè)點(diǎn)位連續(xù)三年的|r-HPI|值,以評估地表水的重金屬污染程度。2015a的樣品中|r-HPI|保持在0~0.8范圍,平均值為0.14,重金屬含量偏低,但在2016年,部分點(diǎn)位重金屬含量升高,尤其是25、36號點(diǎn)位,|r-HPI|值達(dá)到了0.5以上,該項(xiàng)研究中的地區(qū)擁有較高的工農(nóng)業(yè)生產(chǎn)水平,采礦、煤炭燃燒以及農(nóng)業(yè)耕作等都有可能是造成水體重金屬濃度增加的原因[18]。盡管仍在本研究所規(guī)定的地表水重金屬標(biāo)準(zhǔn)范圍內(nèi),但隨著重金屬濃度的增加,人體和水中動(dòng)植物的健康將受到嚴(yán)重威脅。2017a的整體點(diǎn)位水體重金屬水平良好,盡管在4、24等點(diǎn)位重金屬濃度較高,但|r-HPI|值仍符合本研究規(guī)定的水體重金屬良好范疇內(nèi)。

    2.2 |r-HPI|的空間分布

    利用克里金插值技術(shù)在GIS平臺上繪制了2015~2017a間|r-HPI|的地理空間變化圖,以觀察整個(gè)研究區(qū)域內(nèi)水體重金屬水平的年際變化和空間分布??死锝鸱ㄊ亲詈玫木€性無偏技術(shù)之一,它可以在稀疏的數(shù)據(jù)區(qū)域中制作出令人滿意的空間地圖,并給出地圖的隨機(jī)歧義性[19]。研究區(qū)連續(xù)三年的|r-HPI|空間分布變化如圖3所示,整體而言,三年來研究區(qū)|r-HPI|的空間分布趨勢高度相似,|r-HPI|熱點(diǎn)出現(xiàn)在研究區(qū)西部。

    2.3 |r-HPI|與HEI、HPI結(jié)果比較

    本研究中提出的|r-HPI|方法評估地表水中重金屬水平的一個(gè)優(yōu)勢是在確保評估指標(biāo)穩(wěn)定性的前提下考慮到各個(gè)參數(shù)的差異,盡管降低了任何低于理想極限的變量以及高于允許極限的變量對指標(biāo)的影響,但“特殊值”并不因此被忽略。此外,即使一個(gè)變量變化,修改后的|r-HPI|結(jié)果也提供了準(zhǔn)確而科學(xué)的評估。彌補(bǔ)了一些現(xiàn)有的評估方法的不足,因?yàn)樗鼈儾⑽纯紤]單因素評估。

    通過對該區(qū)域數(shù)據(jù)點(diǎn)進(jìn)行HEI、HPI及新提出|r-HPI|三種方法的評估,從表1中可看出,HPI在區(qū)域中的值并沒有顯著差異,均顯示低污染水平。HEI值在1、2、23、24、25、35號點(diǎn)位出現(xiàn)了高值,由此方法可以判斷對應(yīng)點(diǎn)位所在區(qū)域受到較大程度的污染,其余點(diǎn)位也受到一定程度的污染。所有點(diǎn)位的 r-HPI 均小于1,表明研究區(qū)處于重金屬低污染水平,與HPI的結(jié)果一致。進(jìn)一步發(fā)現(xiàn),當(dāng)各重金屬指標(biāo)中有兩個(gè)指標(biāo)超過《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)(GB3838-2002)》I級標(biāo)準(zhǔn)時(shí), r-HPI 范圍保持在0.4~0.58之間,只有一個(gè)指標(biāo)為II級標(biāo)準(zhǔn)時(shí), r-HPI 小于0.4,因此, r-HPI 值準(zhǔn)確反映出了重金屬各指標(biāo)的超標(biāo)情況。

    2.4 |r-HPI|有效性檢驗(yàn)

    Bland-Altman適用于連續(xù)變量的兩種測量方法進(jìn)行一致性比較,通過繪制散點(diǎn)圖量化測量方法之間的差異。由于本研究中測量方法的結(jié)果不同,故先對原始數(shù)據(jù)進(jìn)行了標(biāo)準(zhǔn)化處理,并使用Kolmogorov-Smirnov檢驗(yàn)數(shù)據(jù)的正態(tài)性,結(jié)果表明,無論是HPI、HEI還是本研究中提出的|r-HPI|方法,其對應(yīng)的K-S檢驗(yàn)的漸進(jìn)相關(guān)性均為0,顯著小于標(biāo)準(zhǔn)值(0.05),由此得出,以上三種水體重金屬評價(jià)方法之間的差異具有顯著的統(tǒng)計(jì)意義。在95%的一致限(limits of agreement, LoA)內(nèi)觀察到最大差異(見圖4),其偏差值為8.333×10-7(見圖4b)和-2.778×10-7(見圖4c)。表明使用|r-HPI|代替HPI的錯(cuò)誤概率僅為0.000 08%,代替HEI錯(cuò)誤概率僅為0.000 03%,可忽略不計(jì)。在HPI與HEI組合的情況下,錯(cuò)誤概率為0.000 1%,差值也分散在95%一致限之間,偏差值為1.111×10-6(見圖4a)。該分析表明,|r-HPI|可以替代現(xiàn)有方法,因?yàn)閨r-HPI|對任何定性參數(shù)均有用。

    3 結(jié)論

    運(yùn)用重金屬污染指數(shù)(HPI)、重金屬評估指數(shù)(HEI)、修正重金屬污染指數(shù)(|r-HPI|)分析研究區(qū)域連續(xù)三年共1 164個(gè)地表水樣品(點(diǎn)位數(shù)為36),表明|r-HPI|方法使用方便,具有靈活性。基于普通克里金插值的地理空間圖清楚地反映了|r-HPI|的年際變化,整體而言,三年來研究區(qū)|r-HPI|空間分布趨勢高度相似,|r-HPI|熱點(diǎn)出現(xiàn)在研究區(qū)西部。通過Kolmogorov-Smirnov檢驗(yàn)數(shù)據(jù)集的正態(tài)性,觀察到顯著差異(P<0.05)。|r-HPI|通過Bland-Altman分析并驗(yàn)證了其有效性,表明它在地表水重金屬污染研究中的良好適用性。

    參考文獻(xiàn):

    [1] 李世龍, 趙增鋒, 邱小琮, 等. 寧夏清水河流域重金屬分布特征及風(fēng)險(xiǎn)評價(jià)[J]. 灌溉排水學(xué)報(bào), 2020,39(7): 128-137.

    [2] 馬明真, 高揚(yáng), 宋賢威, 等. 鄱陽湖地區(qū)多尺度流域水體重金屬輸送特征及其污染風(fēng)險(xiǎn)評價(jià)[J]. 生態(tài)學(xué)報(bào), 2019,39(17): 6 404-6 415.

    [3] TIWARI A K,SINGH P K,SINGH A K,et al.Estimation of Heavy Metal Contamination in Groundwater and Development of a Heavy Metal Pollution Index by Using GIS Technique[J]. Bulletin of Environmental Contamination and Toxicology, 2016,96(4): 508-515.

    [4] DASH S, BORAH S S, KALAMDHAD A. A modified indexing approach for assessment of heavy metal contamination in Deepor Beel, India[J]. Ecological Indicators, 2019,106(11): 105444.1-105444.12.

    [5] GIRI S,SINGH A K.Assessment of Surface Water Quality Using Heavy Metal Pollution Index in Subarnarekha River,India[J]. Water Quality Exposure & Health, 2014,5(4): 173-182.

    [6] BHARDWAJ R, GUPTA A, GARG J K. Evaluation of heavy metal contamination using environmetrics and indexing approach for River Yamuna, Delhi stretch, India[J]. Water Science, 2017,31(1): 52-66.

    [7] KUMAR V, PARIHAR R D, SHARMA A, et al. Global evaluation of heavy metal content in surface water bodies: A meta-analysis using heavy metal pollution indices and multivariate statistical analyses[J]. Chemosphere, 2019,236: 124364.

    [8] HOSSAIN M, PATRA P K. Water pollution index – A new integrated approach to rank water quality[J]. Ecological Indicators, 2020,117: 106668.

    [9] SINGH R, VENKATESH A S, SYED T H, et al. Assessment of potentially toxic trace elements contamination in groundwater resources of the coal mining area of the Korba Coalfield, Central India[J]. Environmental Earth Sciences, 2017,76(16): 566.

    [10] BHUIYAN M A H, PARVEZ L, ISLAM M A, et al. Heavy metal pollution of coal mine-affected agricultural soils in the northern part of Bangladesh[J]. Journal of Hazardous Materials, 2010,173(1): 384-392.

    [11] HEROJEET R,RISHI M S,KISHORE N.Integrated approach of heavy metal pollution indices and complexity quantification using chemometric models in the Sirsa Basin,Nalagarh valley,Himachal Pradesh, India[J].Chinese Journal of Geochemistry,2015,34(4):1-14.

    [12] 楊剛, 伍鈞, 孫百曄, 等. 雅安市耕地土壤重金屬健康風(fēng)險(xiǎn)評價(jià)[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2010,29(S1): 74-79.

    [13] 郝守寧, 李成林, 郭永剛. 西藏尼洋河流域飲用水源地重金屬健康風(fēng)險(xiǎn)評價(jià)[J]. 環(huán)境科學(xué)與技術(shù), 2020,43(5): 154-163.

    [14] REZAEI A, HASSANI H, HASSANI S, et al. Evaluation of groundwater quality and heavy metal pollution indices in Bazman basin, southeastern Iran[J].Groundwater for Sustainable Development, 2019,9: 100245.

    [15] EDET A E, OFFIONG O E. Evaluation of water quality pollution indices for heavy metal contamination monitoring. A study case from Akpabuyo-Odukpani area, Lower Cross River Basin (southeastern Nigeria)[J]. GeoJournal, 2002,57(4): 295-304.

    [16] PHILIP L V. Frederic R. Siegel: Environmental Geochemistry of Potentially Toxic Metals[J]. Mineralium Deposita, 2003,38(2):258.

    [17] BHUIYAN M A H,PARVEZ L,ISLAM M A,et al.Heavy metal pollution of coal mine-affected agricultural soils in the northern part of Bangladesh[J].Journal of hazardous materials, 2010,173(1-3):384-392.

    [18] XU M, WANG R, YANG X, et al. Spatial distribution and ecological risk assessment of heavy metal pollution in surface sediments from shallow lakes in East China[J]. Journal of Geochemical Exploration, 2020,213: 106490.

    [19] GUISANDE C, RUEDA-QUECHO A J, RANGEL-SILVA F A, et al. SINENVAP: An algorithm that employs kriging to identify optimal spatial interpolation models in polygons[J]. Ecological Informatics, 2019,53(C): 100975.

    (責(zé)任編輯:李 麗)

    猜你喜歡
    污染方法研究
    FMS與YBT相關(guān)性的實(shí)證研究
    遼代千人邑研究述論
    視錯(cuò)覺在平面設(shè)計(jì)中的應(yīng)用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    EMA伺服控制系統(tǒng)研究
    堅(jiān)決打好污染防治攻堅(jiān)戰(zhàn)
    堅(jiān)決打好污染防治攻堅(jiān)戰(zhàn)
    可能是方法不對
    用對方法才能瘦
    Coco薇(2016年2期)2016-03-22 02:42:52
    四大方法 教你不再“坐以待病”!
    Coco薇(2015年1期)2015-08-13 02:47:34
    捕魚
    午夜精品国产一区二区电影| 国产亚洲精品久久久久5区| 如日韩欧美国产精品一区二区三区| 中出人妻视频一区二区| 香蕉久久夜色| 久久人人97超碰香蕉20202| 午夜福利,免费看| 久久精品国产亚洲av高清一级| 很黄的视频免费| 欧美大码av| 99国产精品一区二区三区| 乱人伦中国视频| 久久中文看片网| 母亲3免费完整高清在线观看| 女人被狂操c到高潮| 色综合婷婷激情| 夜夜看夜夜爽夜夜摸 | 天天添夜夜摸| 精品久久久久久电影网| 国产亚洲精品一区二区www| 天天躁夜夜躁狠狠躁躁| 久久久水蜜桃国产精品网| 嫁个100分男人电影在线观看| 久久伊人香网站| 免费在线观看视频国产中文字幕亚洲| 国产真人三级小视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 超碰97精品在线观看| av视频免费观看在线观看| 69av精品久久久久久| 十分钟在线观看高清视频www| 99国产精品99久久久久| 91字幕亚洲| 久久精品影院6| 欧美精品一区二区免费开放| 丁香六月欧美| 韩国精品一区二区三区| 久久久久久大精品| 欧美中文日本在线观看视频| 美女大奶头视频| videosex国产| 首页视频小说图片口味搜索| 美国免费a级毛片| 亚洲伊人色综图| 女性生殖器流出的白浆| 色婷婷av一区二区三区视频| 男男h啪啪无遮挡| 老熟妇仑乱视频hdxx| 天堂√8在线中文| 国产欧美日韩一区二区精品| 亚洲精品久久午夜乱码| 天天躁夜夜躁狠狠躁躁| 一个人免费在线观看的高清视频| 欧美黄色淫秽网站| 国产免费男女视频| 久久国产精品男人的天堂亚洲| 香蕉丝袜av| 日韩一卡2卡3卡4卡2021年| 麻豆一二三区av精品| 欧美一区二区精品小视频在线| 黄色视频不卡| 亚洲欧美激情综合另类| 免费日韩欧美在线观看| av视频免费观看在线观看| 久久天躁狠狠躁夜夜2o2o| 精品午夜福利视频在线观看一区| 在线十欧美十亚洲十日本专区| 欧美av亚洲av综合av国产av| 老鸭窝网址在线观看| 午夜久久久在线观看| 又黄又爽又免费观看的视频| 久久精品国产亚洲av香蕉五月| 啪啪无遮挡十八禁网站| 高潮久久久久久久久久久不卡| 在线观看www视频免费| 国产91精品成人一区二区三区| 啪啪无遮挡十八禁网站| 国产精品一区二区在线不卡| 精品一区二区三区视频在线观看免费 | a级毛片在线看网站| 丰满人妻熟妇乱又伦精品不卡| 亚洲五月色婷婷综合| 国产精品国产av在线观看| 亚洲一区二区三区欧美精品| 国产精品一区二区在线不卡| 嫩草影视91久久| 五月开心婷婷网| 69精品国产乱码久久久| 一级毛片高清免费大全| 婷婷精品国产亚洲av在线| 日韩视频一区二区在线观看| 一级片'在线观看视频| 水蜜桃什么品种好| 亚洲av成人av| e午夜精品久久久久久久| 咕卡用的链子| 国产xxxxx性猛交| 国产精品av久久久久免费| 亚洲成a人片在线一区二区| 亚洲自偷自拍图片 自拍| 亚洲第一av免费看| 国产成人精品久久二区二区91| 欧美另类亚洲清纯唯美| 国产无遮挡羞羞视频在线观看| 亚洲熟妇熟女久久| 最新在线观看一区二区三区| 国产视频一区二区在线看| 岛国视频午夜一区免费看| 后天国语完整版免费观看| 黄色毛片三级朝国网站| 啦啦啦免费观看视频1| 免费女性裸体啪啪无遮挡网站| 成人亚洲精品av一区二区 | 一进一出抽搐动态| 国产成+人综合+亚洲专区| 亚洲第一av免费看| 国产av又大| 免费一级毛片在线播放高清视频 | av在线天堂中文字幕 | 亚洲精品在线观看二区| 亚洲精华国产精华精| 欧美大码av| 12—13女人毛片做爰片一| 色综合站精品国产| 日本三级黄在线观看| 中文亚洲av片在线观看爽| 天堂动漫精品| 老司机午夜福利在线观看视频| 9热在线视频观看99| 丰满的人妻完整版| 欧美日韩视频精品一区| 一本大道久久a久久精品| 丁香欧美五月| 亚洲人成电影观看| 午夜91福利影院| 国产在线观看jvid| 免费人成视频x8x8入口观看| 免费观看精品视频网站| 免费看a级黄色片| 久久99一区二区三区| 国产片内射在线| 女人被躁到高潮嗷嗷叫费观| 成在线人永久免费视频| 国内久久婷婷六月综合欲色啪| 亚洲三区欧美一区| 黄色视频,在线免费观看| 99久久综合精品五月天人人| 精品少妇一区二区三区视频日本电影| 一二三四社区在线视频社区8| 免费在线观看黄色视频的| 午夜福利影视在线免费观看| 午夜福利,免费看| 久久久精品国产亚洲av高清涩受| 国产精品免费视频内射| 久久久精品欧美日韩精品| 成熟少妇高潮喷水视频| 国产无遮挡羞羞视频在线观看| 淫秽高清视频在线观看| 亚洲欧洲精品一区二区精品久久久| 老司机靠b影院| 成在线人永久免费视频| 午夜久久久在线观看| 超碰97精品在线观看| 成在线人永久免费视频| 黄片播放在线免费| 国产又色又爽无遮挡免费看| av福利片在线| 欧美精品亚洲一区二区| 免费高清在线观看日韩| 久久影院123| 精品电影一区二区在线| 色精品久久人妻99蜜桃| 一区二区三区国产精品乱码| 他把我摸到了高潮在线观看| 女人精品久久久久毛片| 激情视频va一区二区三区| 色综合站精品国产| 制服人妻中文乱码| 国产亚洲精品第一综合不卡| 麻豆国产av国片精品| 9191精品国产免费久久| 亚洲五月色婷婷综合| 亚洲中文av在线| av片东京热男人的天堂| 成人特级黄色片久久久久久久| 欧美日韩中文字幕国产精品一区二区三区 | 真人一进一出gif抽搐免费| 亚洲av片天天在线观看| 国产欧美日韩一区二区三| 成年人黄色毛片网站| 丝袜美足系列| 国产精品免费一区二区三区在线| 日本vs欧美在线观看视频| 成年版毛片免费区| 淫秽高清视频在线观看| 亚洲情色 制服丝袜| 亚洲精品在线观看二区| 亚洲色图 男人天堂 中文字幕| 伦理电影免费视频| 18禁观看日本| 99riav亚洲国产免费| 中文字幕高清在线视频| 丰满的人妻完整版| 老司机午夜福利在线观看视频| 久久精品91蜜桃| 天天添夜夜摸| av国产精品久久久久影院| 国产精品久久久av美女十八| 香蕉国产在线看| 又黄又粗又硬又大视频| 亚洲熟妇熟女久久| 国产精品一区二区在线不卡| 久久国产乱子伦精品免费另类| 大香蕉久久成人网| 不卡av一区二区三区| 免费av毛片视频| 亚洲欧美激情综合另类| 丝袜在线中文字幕| 久久精品国产亚洲av高清一级| 人妻丰满熟妇av一区二区三区| 色在线成人网| 香蕉国产在线看| 不卡av一区二区三区| 一二三四社区在线视频社区8| 精品午夜福利视频在线观看一区| 波多野结衣高清无吗| 母亲3免费完整高清在线观看| 麻豆一二三区av精品| 少妇被粗大的猛进出69影院| 中出人妻视频一区二区| 男人舔女人的私密视频| 嫁个100分男人电影在线观看| 国产伦一二天堂av在线观看| 亚洲成人国产一区在线观看| 欧美一级毛片孕妇| 婷婷精品国产亚洲av在线| 高清毛片免费观看视频网站 | 亚洲熟妇熟女久久| 欧美中文综合在线视频| 黑人巨大精品欧美一区二区蜜桃| 99riav亚洲国产免费| 国产成+人综合+亚洲专区| www.999成人在线观看| 久久草成人影院| 免费在线观看影片大全网站| 久久精品91无色码中文字幕| 老司机深夜福利视频在线观看| 天堂动漫精品| 夫妻午夜视频| 日本欧美视频一区| 很黄的视频免费| 精品国产亚洲在线| 老鸭窝网址在线观看| 久久精品国产综合久久久| 美女高潮喷水抽搐中文字幕| 丰满人妻熟妇乱又伦精品不卡| 欧美日本亚洲视频在线播放| 国产精品秋霞免费鲁丝片| 午夜两性在线视频| 日日干狠狠操夜夜爽| 国产精品野战在线观看 | 麻豆一二三区av精品| 国产欧美日韩一区二区三| 国产成人欧美在线观看| 女人高潮潮喷娇喘18禁视频| 日本一区二区免费在线视频| 亚洲一区二区三区色噜噜 | 成人永久免费在线观看视频| 久久久久久大精品| 日日摸夜夜添夜夜添小说| 亚洲精品久久成人aⅴ小说| 极品人妻少妇av视频| 不卡一级毛片| 高清毛片免费观看视频网站 | 99国产综合亚洲精品| 午夜免费鲁丝| 国产蜜桃级精品一区二区三区| e午夜精品久久久久久久| 每晚都被弄得嗷嗷叫到高潮| 久久午夜亚洲精品久久| 亚洲成人免费电影在线观看| 电影成人av| 欧美乱码精品一区二区三区| videosex国产| 91精品国产国语对白视频| 久久天堂一区二区三区四区| 男女高潮啪啪啪动态图| 婷婷六月久久综合丁香| 性欧美人与动物交配| 男人舔女人的私密视频| 操出白浆在线播放| 亚洲成人久久性| 国产成人影院久久av| 色老头精品视频在线观看| 亚洲中文av在线| 久久人人爽av亚洲精品天堂| 亚洲视频免费观看视频| 夫妻午夜视频| 日韩欧美一区视频在线观看| 欧美成人免费av一区二区三区| www.自偷自拍.com| 中文字幕人妻熟女乱码| 身体一侧抽搐| 午夜精品在线福利| 动漫黄色视频在线观看| 青草久久国产| 久久久久久免费高清国产稀缺| 日韩av在线大香蕉| 色婷婷久久久亚洲欧美| 天堂中文最新版在线下载| 在线观看免费高清a一片| 欧美黑人欧美精品刺激| 黄色a级毛片大全视频| 免费av毛片视频| 亚洲男人的天堂狠狠| 国产成人系列免费观看| 在线av久久热| 日本五十路高清| 日韩有码中文字幕| 久久久久久人人人人人| 麻豆国产av国片精品| 久久久精品欧美日韩精品| 国产精品电影一区二区三区| 欧美日本中文国产一区发布| 怎么达到女性高潮| 亚洲精品美女久久av网站| 久久午夜综合久久蜜桃| 少妇被粗大的猛进出69影院| 亚洲精品国产色婷婷电影| 国产不卡一卡二| 最新在线观看一区二区三区| 亚洲欧美日韩无卡精品| www国产在线视频色| 精品人妻1区二区| 人人妻人人澡人人看| 亚洲国产欧美日韩在线播放| 麻豆国产av国片精品| 久久久久久免费高清国产稀缺| 国产成人影院久久av| 亚洲专区国产一区二区| 免费看十八禁软件| 亚洲熟女毛片儿| 国产精品 国内视频| 国产成人精品久久二区二区免费| 免费高清在线观看日韩| 咕卡用的链子| 欧美日韩乱码在线| 又黄又粗又硬又大视频| 亚洲国产精品一区二区三区在线| 手机成人av网站| 18禁美女被吸乳视频| 欧美日韩国产mv在线观看视频| 91成年电影在线观看| 91大片在线观看| 80岁老熟妇乱子伦牲交| 国产高清国产精品国产三级| 在线免费观看的www视频| 天堂影院成人在线观看| 日韩欧美三级三区| 人人妻人人添人人爽欧美一区卜| 人人妻人人爽人人添夜夜欢视频| 一级毛片女人18水好多| 男女床上黄色一级片免费看| 搡老乐熟女国产| 国产欧美日韩一区二区三| 亚洲欧洲精品一区二区精品久久久| 国产av一区二区精品久久| 欧美中文综合在线视频| 国产一卡二卡三卡精品| 成人影院久久| 日本精品一区二区三区蜜桃| 日日摸夜夜添夜夜添小说| 亚洲精品国产精品久久久不卡| 一级毛片精品| 一个人免费在线观看的高清视频| 国产伦一二天堂av在线观看| 国产成人精品无人区| 日韩精品中文字幕看吧| 天堂中文最新版在线下载| 欧美性长视频在线观看| 啪啪无遮挡十八禁网站| 欧美成人午夜精品| 成年人黄色毛片网站| 999精品在线视频| 亚洲一区中文字幕在线| 99国产精品免费福利视频| 亚洲欧美精品综合久久99| 欧美亚洲日本最大视频资源| 99riav亚洲国产免费| 午夜激情av网站| 大陆偷拍与自拍| 91老司机精品| 免费在线观看影片大全网站| 天堂中文最新版在线下载| 中文字幕最新亚洲高清| 午夜福利一区二区在线看| 免费av中文字幕在线| 他把我摸到了高潮在线观看| 午夜免费鲁丝| 91成年电影在线观看| 国产成人精品久久二区二区免费| 少妇裸体淫交视频免费看高清 | 十八禁网站免费在线| 国产激情欧美一区二区| 日韩视频一区二区在线观看| 老司机亚洲免费影院| 精品福利永久在线观看| 日韩中文字幕欧美一区二区| 91国产中文字幕| 国产成人欧美| 精品久久久久久电影网| 国产真人三级小视频在线观看| 亚洲avbb在线观看| 黄片小视频在线播放| 久久人人爽av亚洲精品天堂| 久久久久久久午夜电影 | а√天堂www在线а√下载| 老熟妇乱子伦视频在线观看| 丝袜人妻中文字幕| 黄色 视频免费看| 免费在线观看黄色视频的| 人人妻,人人澡人人爽秒播| 亚洲成国产人片在线观看| 国产又色又爽无遮挡免费看| 欧美大码av| 一二三四社区在线视频社区8| 久久国产亚洲av麻豆专区| 黄色视频,在线免费观看| 亚洲va日本ⅴa欧美va伊人久久| 人成视频在线观看免费观看| 成人特级黄色片久久久久久久| 在线看a的网站| 超碰成人久久| 中亚洲国语对白在线视频| 丰满饥渴人妻一区二区三| 真人一进一出gif抽搐免费| 男人舔女人的私密视频| 老司机午夜福利在线观看视频| 啦啦啦免费观看视频1| 高清在线国产一区| 天堂√8在线中文| 中文字幕高清在线视频| a级片在线免费高清观看视频| 一级片'在线观看视频| 他把我摸到了高潮在线观看| 窝窝影院91人妻| 免费高清在线观看日韩| 9热在线视频观看99| 大香蕉久久成人网| 黑人巨大精品欧美一区二区蜜桃| 老司机福利观看| 宅男免费午夜| 国产成人精品久久二区二区免费| 高清黄色对白视频在线免费看| 亚洲精品成人av观看孕妇| 男人操女人黄网站| av天堂在线播放| 婷婷丁香在线五月| 99热只有精品国产| 国产亚洲精品综合一区在线观看 | 一区二区三区激情视频| 国产成人一区二区三区免费视频网站| 搡老岳熟女国产| 正在播放国产对白刺激| 国产精品日韩av在线免费观看 | 国产不卡一卡二| 国产aⅴ精品一区二区三区波| 亚洲精品粉嫩美女一区| 不卡av一区二区三区| 日韩高清综合在线| 欧美+亚洲+日韩+国产| 久久久国产成人精品二区 | 日日爽夜夜爽网站| 亚洲欧美一区二区三区黑人| 久久精品91蜜桃| 日本wwww免费看| 日本欧美视频一区| 国产成人av激情在线播放| 亚洲av电影在线进入| 欧美日韩瑟瑟在线播放| 午夜免费激情av| 亚洲全国av大片| 欧美中文日本在线观看视频| 国产精品免费一区二区三区在线| 真人做人爱边吃奶动态| 国产成人精品无人区| 在线观看www视频免费| 老司机午夜十八禁免费视频| 一区福利在线观看| 亚洲在线自拍视频| 日韩视频一区二区在线观看| 变态另类成人亚洲欧美熟女 | 丝袜美足系列| 久久久久久久久久久久大奶| 丝袜在线中文字幕| 精品久久蜜臀av无| 亚洲成av片中文字幕在线观看| 韩国av一区二区三区四区| 欧美在线一区亚洲| 久久久久久人人人人人| 一级片'在线观看视频| 中文字幕人妻丝袜一区二区| 国产伦一二天堂av在线观看| 久久久久九九精品影院| 亚洲精品在线观看二区| 黄色丝袜av网址大全| 久久久久久久午夜电影 | 成人黄色视频免费在线看| 91成人精品电影| 成人18禁在线播放| 成人亚洲精品一区在线观看| 国产精品99久久99久久久不卡| 国产欧美日韩一区二区精品| 一级a爱片免费观看的视频| 久久午夜综合久久蜜桃| 俄罗斯特黄特色一大片| 亚洲欧美激情在线| 老司机午夜十八禁免费视频| 欧美午夜高清在线| 在线观看舔阴道视频| www.www免费av| 黑人巨大精品欧美一区二区蜜桃| 色尼玛亚洲综合影院| 老汉色av国产亚洲站长工具| 免费日韩欧美在线观看| 69精品国产乱码久久久| 在线观看一区二区三区| 国产激情欧美一区二区| 国产精品久久久久久人妻精品电影| 女人被躁到高潮嗷嗷叫费观| 高清av免费在线| 国产午夜精品久久久久久| 久久精品亚洲熟妇少妇任你| 欧美色视频一区免费| 波多野结衣av一区二区av| 黑人欧美特级aaaaaa片| 黄频高清免费视频| 黄色 视频免费看| 自拍欧美九色日韩亚洲蝌蚪91| 男女床上黄色一级片免费看| 两性午夜刺激爽爽歪歪视频在线观看 | 深夜精品福利| 中文字幕精品免费在线观看视频| 青草久久国产| 国产又爽黄色视频| 国产精品电影一区二区三区| 啪啪无遮挡十八禁网站| 韩国精品一区二区三区| 国产高清国产精品国产三级| 欧美激情高清一区二区三区| 999精品在线视频| 欧美成人免费av一区二区三区| 老汉色av国产亚洲站长工具| 欧美不卡视频在线免费观看 | av欧美777| 国产精品野战在线观看 | 国产欧美日韩一区二区三区在线| 亚洲国产中文字幕在线视频| 90打野战视频偷拍视频| 亚洲欧美日韩另类电影网站| 日韩一卡2卡3卡4卡2021年| 他把我摸到了高潮在线观看| 国产亚洲精品一区二区www| 久久性视频一级片| 国产xxxxx性猛交| 可以在线观看毛片的网站| 欧美在线一区亚洲| 国产深夜福利视频在线观看| 最新美女视频免费是黄的| 少妇裸体淫交视频免费看高清 | 高清黄色对白视频在线免费看| 国产欧美日韩一区二区三区在线| 成人18禁在线播放| 亚洲九九香蕉| 精品一区二区三区视频在线观看免费 | 日韩国内少妇激情av| 久久香蕉激情| 国产精品香港三级国产av潘金莲| 99国产精品99久久久久| 国产成人影院久久av| 亚洲国产精品999在线| 午夜免费鲁丝| 俄罗斯特黄特色一大片| 欧美亚洲日本最大视频资源| 国产成人啪精品午夜网站| 视频区欧美日本亚洲| 国产精品亚洲一级av第二区| 中文字幕色久视频| 亚洲精品粉嫩美女一区| 亚洲欧美日韩另类电影网站| 亚洲成国产人片在线观看| 亚洲精品粉嫩美女一区| 变态另类成人亚洲欧美熟女 | 老司机亚洲免费影院| 操出白浆在线播放| 午夜免费观看网址| 亚洲va日本ⅴa欧美va伊人久久| 99热国产这里只有精品6| 91麻豆精品激情在线观看国产 | 十八禁人妻一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品野战在线观看 | 国产精品影院久久| 黄片播放在线免费| 久久久久久久久久久久大奶| 涩涩av久久男人的天堂| 久久青草综合色| 久久精品成人免费网站| 国产在线观看jvid| 男男h啪啪无遮挡| 久久九九热精品免费| 欧美黑人精品巨大| 国产伦一二天堂av在线观看| 日韩有码中文字幕| 欧美老熟妇乱子伦牲交|