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

    內(nèi)梅羅指數(shù)法在環(huán)境質(zhì)量評價中的適用性與修正原則

    2017-12-25 00:58:52韓術(shù)鑫王利紅趙長盛
    關(guān)鍵詞:內(nèi)梅羅環(huán)境要素適用性

    韓術(shù)鑫,王利紅,趙長盛

    內(nèi)梅羅指數(shù)法在環(huán)境質(zhì)量評價中的適用性與修正原則

    韓術(shù)鑫,王利紅,趙長盛

    (山東省分析測試中心,濟(jì)南 250014)

    為研究內(nèi)梅羅指數(shù)法在環(huán)境要素質(zhì)量評價中的適用性及確定其修正原則,將方法的計(jì)算公式進(jìn)行了數(shù)學(xué)推導(dǎo)。結(jié)果表明:典型內(nèi)梅羅指數(shù)法的計(jì)算公式可變換為線性表達(dá)式,其斜率和截距分別為最大值(Fj,max)和算術(shù)平均值(Fj,ave)的權(quán)重系數(shù);隨著Fj,max與Fj,ave比值的增大,F(xiàn)j,max權(quán)重增大,F(xiàn)j,ave權(quán)重減小。參評污染因子數(shù)量(n)對方法的評分結(jié)果具有重要影響,典型內(nèi)梅羅指數(shù)法更適用于區(qū)域中目的性和指向性較明確或總體污染程度較輕的環(huán)境要素質(zhì)量評價。為提高方法的適用性,確定了2種修正內(nèi)梅羅指數(shù)法的劃分原則,Ⅰ型和Ⅱ型修正內(nèi)梅羅指數(shù)法分別以n和Fj,max/Fj,ave等于5作為劃分結(jié)點(diǎn),其中:Ⅰ型修正內(nèi)梅羅指數(shù)法傾向于參評因子不多、評價目的性和指向性較強(qiáng)的環(huán)境要素質(zhì)量評價;Ⅱ型修正內(nèi)梅羅指數(shù)法是對Ⅰ型修正內(nèi)梅羅指數(shù)法的拓展和延伸,更適用于多因子綜合性質(zhì)量評價。由內(nèi)梅羅指數(shù)法線性表達(dá)式的數(shù)學(xué)內(nèi)涵可知,在n或Fj,max/Fj,ave≤5時,方法可不作修正,此時并不影響評價結(jié)果的客觀性和準(zhǔn)確性。

    內(nèi)梅羅指數(shù)法;環(huán)境要素;方法適用性;方法修正原則;污染因子權(quán)重

    內(nèi)梅羅指數(shù)法由單因子指數(shù)法發(fā)展而來,是當(dāng)前國內(nèi)外進(jìn)行綜合污染指數(shù)計(jì)算最常用的方法之一。該方法是一種兼顧極值或突出最大值的計(jì)權(quán)型多因子環(huán)境質(zhì)量指數(shù)[1],能夠較全面地反映各種污染物的協(xié)同作用,從而較準(zhǔn)確地評價各類環(huán)境要素的污染程度。該方法最早由美國雪城大學(xué)(Syracuse University)內(nèi)梅羅(N.L.Nemerow)教授在其所著的《河流污染科學(xué)分析》一書中提出。自20世紀(jì)70年代以來,被我國一些環(huán)??蒲腥藛T應(yīng)用于某些地區(qū)的水污染評價中[2]。發(fā)展至今,該方法在近岸海域海水[2-5]、河流[6-10]、湖泊[11-12]、濕地水環(huán)境[13-14]、飲用水源地水環(huán)境[15-16]、地下水[17-26]等水質(zhì)類環(huán)境質(zhì)量評價,以及河道沉積物[27]、海洋沉積物[28-30]、土壤[1,31-38]等環(huán)境要素中重金屬、有機(jī)氯、多氯聯(lián)苯、多環(huán)芳烴等高毒高危害性污染物的環(huán)境健康風(fēng)險(xiǎn)評價中均作為重要的評價方法廣泛應(yīng)用,倍受科研工作者的青睞。

    內(nèi)梅羅指數(shù)法具有計(jì)算過程簡潔、物理概念清晰、可操作性強(qiáng)等優(yōu)點(diǎn)[18],但在實(shí)際應(yīng)用過程中出現(xiàn)了一些問題,導(dǎo)致評價結(jié)果有時不能準(zhǔn)確真實(shí)地反映所評價的環(huán)境要素質(zhì)量。主要體現(xiàn)在如下幾個方面:(1)評價時只考慮了單因子污染指數(shù)算術(shù)平均值和最大值,且最大值權(quán)重過高[8,16,18-25];(2)未考慮參評污染因子的毒性、危害性在評價的環(huán)境要素中的權(quán)重[9,14,15,25-26,32,37];(3)未揭示參評污染因子種類和數(shù)量對評價值的影響。針對以上問題,許多科研人員通過各種途徑和方法對內(nèi)梅羅指數(shù)進(jìn)行改進(jìn)或修正,取得了一定的成效,但并未從根本上講清楚該方法應(yīng)用于環(huán)境要素的適用性問題以及進(jìn)行修正的原則劃分。本文通過對典型內(nèi)梅羅指數(shù)法計(jì)算公式的數(shù)學(xué)推導(dǎo),探討了該方法的數(shù)學(xué)原理,并以此為依據(jù)深入分析了其在環(huán)境要素質(zhì)量評價中的適用性。另外,劃分了該方法修正與否的基本原則,為科研人員和環(huán)保決策部門開展環(huán)境質(zhì)量評價提供實(shí)用性較強(qiáng)的操作方法。

    1 材料與方法

    1.1 典型內(nèi)梅羅指數(shù)法

    典型內(nèi)梅羅指數(shù)法包含單因子污染指數(shù)的算術(shù)平均值和最大值,單因子污染指數(shù)法和內(nèi)梅羅指數(shù)法的計(jì)算公式分別如下:

    單因子指數(shù)法:Fij=Ci/Sij

    式中:Fj為按某環(huán)境要素質(zhì)量標(biāo)準(zhǔn)中第j種標(biāo)準(zhǔn)計(jì)算得出的評分值;Fj,max為在某環(huán)境要素質(zhì)量標(biāo)準(zhǔn)第j種標(biāo)準(zhǔn)下Fij的最大值,F(xiàn)j,max=max{Fij};Fj,ave為在某環(huán)境要素質(zhì)量標(biāo)準(zhǔn)第j種標(biāo)準(zhǔn)下Fij的算術(shù)平均值,F(xiàn)j,ave=環(huán)境要素質(zhì)量標(biāo)準(zhǔn)選取的污染因子所對應(yīng)的質(zhì)量類別,j=1,2,…,m;Ci為第 i種污染因子的實(shí)測濃度;Sij為第i種污染因子在某環(huán)境要素質(zhì)量標(biāo)準(zhǔn)中第j種標(biāo)準(zhǔn)下的標(biāo)準(zhǔn)值。

    1.2 基于方法修正引入的不同環(huán)境要素中污染因子權(quán)重的計(jì)算

    引入污染因子權(quán)重是對典型內(nèi)梅羅指數(shù)法進(jìn)行修正的重大突破,因環(huán)境要素的特點(diǎn)不同,污染因子權(quán)重的計(jì)算方式略有不同。

    1.2.1 水質(zhì)類污染因子權(quán)重的計(jì)算

    由于各類污染物對人體與環(huán)境的毒性與危害程度存在較大的差異性,一般地,某種污染因子的危害性與該污染因子的排放標(biāo)準(zhǔn)基本上呈反比例的內(nèi)在關(guān)系[39]。具體計(jì)算過程如下:首先以環(huán)境要素質(zhì)量標(biāo)準(zhǔn)為依據(jù),確定評價標(biāo)準(zhǔn)級別j;其次,將j類標(biāo)準(zhǔn)等級下各種參評污染因子的排放標(biāo)準(zhǔn)Si按由小到大順序排列 S1,S2,…,Sn,將其最大值 Smax與 Si比較,令 Ri表示第i種污染因子的相關(guān)性比值,可得:

    式中:Wi為第i種污染因子的權(quán)重值,(Ri=Smax/Si),

    1.2.2 沉積物和土壤中污染因子權(quán)重的計(jì)算

    開展地表水、地下水等環(huán)境要素的質(zhì)量評價時,一般參評的污染因子較多且類別分散,評價結(jié)果更傾向于綜合性;而針對土壤和沉積物的質(zhì)量評價主要針對特定的某類污染因子(如重金屬、多環(huán)芳烴等),具有更明顯的指向性。因此,其污染因子的權(quán)重計(jì)算方式差異較大。關(guān)于重金屬污染因子權(quán)重的計(jì)算,簡要介紹幾種方式:

    第一種權(quán)重和類別依據(jù)Swaine提出的環(huán)境影響程度劃分[40],把 Hg、Pb、Cd、As、Cr、Se幾種毒性較高的元素歸為Ⅰ類,Mn、Mo、V、Be、Th、U、Ni、Zn、Cu 幾種毒性較低的元素歸為Ⅱ類,Ⅰ類和Ⅱ類分別設(shè)權(quán)重為3 和 2[37]。其他重金屬元素如 Co、Ba、Ra、Sb、Sn、Ti等劃歸Ⅲ類。

    第二種權(quán)重的計(jì)算通過徐爭啟等[41]應(yīng)用潛在生態(tài)危害指數(shù)法(RI),依據(jù)Hakanson的計(jì)算原則并結(jié)合陳靜生[42]的計(jì)算方法,得到12種重金屬元素的毒性系數(shù)(表1)。

    第三種依據(jù)林麗欽[43]研究的基于毒理學(xué)安全評價數(shù)據(jù)推算出重金屬毒性系數(shù),17種重金屬的規(guī)正化權(quán)重見表2。

    第四種計(jì)算方法是李雪梅等[44]基于改進(jìn)的層次分析法(AHP法),以重金屬污染物在糧食中的限量值為依據(jù),確定了7種重金屬污染因子的權(quán)重系數(shù)(表 3)。

    另外,鄒喬等[32]通過多元線性回歸的方式計(jì)算出16種多環(huán)芳烴濃度權(quán)重,具體見表4。

    2 結(jié)果與討論

    2.1 典型內(nèi)梅羅指數(shù)法計(jì)算公式的線性推導(dǎo)及方法適用性分析

    典型內(nèi)梅羅指數(shù)法的計(jì)算公式等于單因子污染指數(shù)的最大值和算術(shù)平均值平方和一半的幾何平均,在原公式的框架內(nèi),對參數(shù) Fj、Fj,max和 Fj,ave的相關(guān)性分析較為困難。因此,本文對3個參數(shù)進(jìn)行了技術(shù)性見圖1。將Y與X分別用上述變換式帶入公式后可得:Fj=aFj,max+bFj,ave。通過線性推導(dǎo)可知內(nèi)梅羅指數(shù)值實(shí)際上是單因子污染指數(shù)最大值和算術(shù)平均值加權(quán)后的和,斜率a和截距b分別代表的是Fj,max和Fj,ave的權(quán)重系數(shù)。

    圖 1 Fj/Fj,ave 與 Fj,max/Fj,ave 的線性關(guān)系圖Figure 1 Linear relation between Fj/Fj,aveand Fj,max/Fj,ave

    表1 基于RI法確定的12種重金屬毒性系數(shù)Table 1 The toxicity coefficients of 12 heavy metals based on RI method

    表2 基于毒理學(xué)確定的17種重金屬規(guī)正化權(quán)重Table 2 The toxicity coefficients of 17 heavy metals based on toxicology

    表3 基于改進(jìn)AHP法確定的7種重金屬污染因子權(quán)重系數(shù)Table 3 The toxicity coefficients of 7 heavy metals based on revised AHP method

    表4 基于多元線性回歸確定的16種多環(huán)芳烴污染因子權(quán)重系數(shù)Table 4 The toxicity coefficients of 16 PAHs based on multiple linear regression

    表5 典型內(nèi)梅羅指數(shù)法線性推導(dǎo)相關(guān)參數(shù)一覽表Table 5 Parameters of linear derivation based on typical NI method

    表5是X取不同值時斜率a、截距b和線性相關(guān)性r的相應(yīng)數(shù)值。由表5可知:隨著X的增加,F(xiàn)j,max權(quán)重系數(shù)a逐漸增大,F(xiàn)j,ave權(quán)重系數(shù)b逐漸減小,權(quán)重系數(shù)之和不斷減小,權(quán)重系數(shù)之比逐漸增大。這說明單因子污染指數(shù)的最大值與算術(shù)平均值之比越大,則Fj,max的權(quán)重系數(shù)就越大,F(xiàn)j,ave的權(quán)重系數(shù)越小,F(xiàn)j,max在Fj中的作用越突出。更為重要的是,隨著Fj,max與Fj,ave比值的增大,權(quán)重系數(shù)之和與理論權(quán)重系數(shù)之和1.00的差距不斷擴(kuò)大,顯示Fj在某類環(huán)境要素質(zhì)量評價的代表性逐步下降,因?yàn)榇嬖谝豁?xiàng)權(quán)重系數(shù)為[1-(a+b)]的未知參數(shù)沒有參與 Fj的計(jì)算。根據(jù) Fj,ave與Fj,max之間的關(guān)系可知:[Fj,max/n]≤Fj,ave≤Fj,max(n指參評的污染因子的數(shù)量),則 1≤[Fj,max/Fj,ave]≤n,即 1≤X≤n。這說明參評的污染因子數(shù)量n實(shí)際上對典型內(nèi)梅羅指數(shù)法的評分值Fj存在著重要影響。n越大,則出現(xiàn) Fj,max權(quán)重系數(shù)過大,F(xiàn)j,ave權(quán)重系數(shù)過小的風(fēng)險(xiǎn)越高。由于未知參數(shù)的權(quán)重系數(shù)隨著X的增加而不斷增大,其對Fj的貢獻(xiàn)越來越不可忽視,由表5可知:當(dāng)X≥20時,未知參數(shù)的權(quán)重系數(shù)甚至超過了Fj,ave的權(quán)重系數(shù),所以僅以 Fj,max和 Fj,ave兩項(xiàng)參數(shù)參與 Fj的計(jì)算將導(dǎo)致最終的評價結(jié)果可信度降低,適用性下降。

    典型內(nèi)梅羅指數(shù)法的線性推導(dǎo)從數(shù)學(xué)基本原理的角度反映了其內(nèi)部存在的問題實(shí)際上是受多種因素共同影響的,除了 Fj,max和 Fj,ave外,n 的作用往往被評價者所忽略?;谝陨蠁栴}的討論,為了達(dá)到準(zhǔn)確客觀的評價結(jié)果,有必要對典型內(nèi)梅羅指數(shù)法在環(huán)境要素中的適用性進(jìn)行嚴(yán)謹(jǐn)?shù)姆治?。為便于進(jìn)行分析,將[1-(a+b)]稱為未知項(xiàng)權(quán)重系數(shù)c。由表5可分析,導(dǎo)致典型內(nèi)梅羅指數(shù)法出現(xiàn)適用性問題的根本原因是c值對應(yīng)的未知參數(shù)未參與到Fj的計(jì)算中,只要解決了這個問題即可從根本上消除內(nèi)梅羅指數(shù)法的不足。由于準(zhǔn)確推導(dǎo)未知參數(shù)的數(shù)學(xué)表達(dá)式不具可操作性,本研究從相反的角度著手,即盡可能地降低c值,使得未知參數(shù)在整個Fj的計(jì)算過程中的貢獻(xiàn)值可忽略不計(jì),從而確保最終評價結(jié)果客觀、可信。

    由表 5 可知:當(dāng) n=5、10、20、30 時,c值分別為0.041、0.096、0.153、0.182。所以,要維持典型內(nèi)梅羅指數(shù)法的適應(yīng)性可從2個途徑進(jìn)行構(gòu)建:第一,針對某類環(huán)境要素質(zhì)量評價的目的和意義,合理選擇參評因子種類和數(shù)量,選擇有典型代表性的幾項(xiàng)污染因子。例如針對土壤中重金屬的環(huán)境風(fēng)險(xiǎn)評價,一般選擇Hg、Pb、Cd、As、Ni、Zn、Cr、Cu 這 8 種重金屬,為了提高方法適用性可進(jìn)一步根據(jù)重金屬毒性高低篩選出Hg、Pb、Cd、As、Cr這 5 項(xiàng)作為參評污染因子。其他的環(huán)境要素可依此類推。因此,從該角度可知,典型內(nèi)梅羅指數(shù)法更適用于目的性和指向性較明確的環(huán)境要素質(zhì)量評價;第二,不限制參評污染因子數(shù)量,從控制Fj,max與 Fj,ave的比值著手,二者的比值越小,方法的適用性越高。對于未受污染或污染極重的環(huán)境要素評價,二者的比值雖然較小但評價結(jié)果無法分級,是沒有實(shí)際意義的。因此,從反向推論可得典型內(nèi)梅羅指數(shù)法適用于某個區(qū)域中污染程度較輕的環(huán)境要素的質(zhì)量評價,這與一些學(xué)者的研究結(jié)論相吻合[13]。

    2.2 典型內(nèi)梅羅指數(shù)法修正原則的劃分

    在實(shí)際環(huán)境要素質(zhì)量評價過程中,出現(xiàn)了典型內(nèi)梅羅指數(shù)法適用性不強(qiáng)的情況,將導(dǎo)致評價結(jié)果不夠客觀公正,說服力不足。許多科研工作者通過多種途徑對該方法進(jìn)行改進(jìn)或修正,提高了方法適用性,獲得了較為真實(shí)客觀可信的評價結(jié)果。但典型內(nèi)梅羅指數(shù)法應(yīng)該采用怎樣的修正原則,缺乏系統(tǒng)地、科學(xué)地分析。為此,我們以該法數(shù)學(xué)推導(dǎo)得到的線性方程式為依據(jù),結(jié)合方法適應(yīng)性的基本原則,詳細(xì)探討典型內(nèi)梅羅指數(shù)法的修正原則。

    修正的根本目的是擴(kuò)大方法的適用性,提高評價結(jié)果的準(zhǔn)確性。從上文可知,典型內(nèi)梅羅指數(shù)法存在2個主要的適用性限制。首先,基于第一種適用性分析,從參評污染因子數(shù)量的角度確定如下修正原則,并命名為Ⅰ型修正內(nèi)梅羅指數(shù)法(表6)。從表5可知,有5個參評的污染因子時,出現(xiàn)的極端情況是X=5,此時對應(yīng)的 Fj,max與 Fj,ave的權(quán)重之和為 0.959,權(quán)重之比2.2,算術(shù)平均值的權(quán)重大于0.3,最大值的作用未出現(xiàn)過分突出的情況,且未知項(xiàng)的權(quán)重很低,對總評分的貢獻(xiàn)可以忽略,典型內(nèi)梅羅指數(shù)法可以不作修正即可保證評價結(jié)果的可靠性。因此,為了既突出最大值的作用又使算術(shù)平均值在總評分值中的貢獻(xiàn)適當(dāng),選定n等于5作為修正劃分的結(jié)點(diǎn)[17]。以此為界限劃分出條件(1)和條件(2)。在條件(2)中再以X等于5為結(jié)點(diǎn),劃分成條件(2-1)和(2-2)。值得注意的是,當(dāng)出現(xiàn)條件(2-2)時,為了提高 Fj,ave對 Fj的貢獻(xiàn),防止Fj,max作用過分突出,以及控制未知項(xiàng)的權(quán)重c,污染因子中最有代表性的5項(xiàng)重新計(jì)算出F′j,ave。其篩選過程其實(shí)是另一種形式的加權(quán)過程,篩選原則可以根據(jù)不同污染因子的毒性高低(表1~表4)或者人體對污染因子的敏感程度進(jìn)行。該修正方式使得Fj,max與F′j,ave的比值重新拉回到不修正的條件范圍內(nèi),從而滿足了典型內(nèi)梅羅指數(shù)法的適用性原則。Ⅰ型修正內(nèi)梅羅指數(shù)法的修正策略更加傾向于參評因子不多、評價目的性和指向性較強(qiáng)的環(huán)境要素質(zhì)量評價,在上述情況下,更容易從全部參評因子中選擇出能代表此次評價結(jié)果的污染因子。例如,環(huán)境要素中某些特定類別的污染物(重金屬、多環(huán)芳烴等)的環(huán)境污染風(fēng)險(xiǎn)評價可以用該方法予以修正。

    表6 Ⅰ型修正內(nèi)梅羅指數(shù)法的修正原則Table 6 Modifications ofⅠ-RNI method

    基于上述分析,注重環(huán)境要素多方位多層次多參評因子的綜合性質(zhì)量評價不適用于Ⅰ型修正內(nèi)梅羅指數(shù)法。為此,諸多研究人員打破典型內(nèi)梅羅指數(shù)法的固有計(jì)算方式,通過增加污染因子權(quán)重的方式,對Fj,max降權(quán)和(或)對 Fj,ave加權(quán)的方式,得到Ⅱ型修正內(nèi)梅羅指數(shù)法(表7)。這種修正方式相當(dāng)于在原來方法的基礎(chǔ)上引入1項(xiàng)或幾項(xiàng)新的參數(shù),參與Fj的計(jì)算。由表7可知:Ⅱ型修正內(nèi)梅羅指數(shù)法不限制參評污染因子數(shù)量,主要以 Fj,max/Fj,ave作為修正的必要條件。如上所述,當(dāng)該比值不大于5時,典型內(nèi)梅羅指數(shù)法具有適用性,故不作修正;比值大于5時,相關(guān)學(xué)者采用了從4種不同的修正方式對方法進(jìn)行了修正。其中針對Fj,max降權(quán)的修正方式又包含4類修正表達(dá)式,特別是修正表達(dá)式①在實(shí)際應(yīng)用中最為常見。降權(quán)的總體思路是引入?yún)⒃u污染因子中權(quán)重最大的1項(xiàng)或者幾項(xiàng)的評分值與Fj,max進(jìn)行算術(shù)平均,使其權(quán)重降低一半,另一半權(quán)重賦予污染因子權(quán)重較高的評分值,避免了 Fj,max過高引起的 Fj評分失真的情況。將Fj,ave加權(quán)有2種修正表達(dá)式:一種由算術(shù)平均值修正為加權(quán)平均值,另一種是Ⅰ型修正內(nèi)梅羅指數(shù)法采用的選擇有代表性的5項(xiàng)污染因子參與Fj,ave的計(jì)算,2種不同形式的加權(quán)均增加了污染因子權(quán)重對Fj的貢獻(xiàn),使結(jié)果更加全面客觀地反映環(huán)境要素的質(zhì)量狀況。將Fj,max降權(quán)和Fj,ave加權(quán)結(jié)合起來則雙向增加了污染因子權(quán)重對Fj的影響。以上3種修正方式都基于內(nèi)梅羅指數(shù)法的計(jì)算公式框架不變的前提下進(jìn)行。楊琳[2]則突破了原先的計(jì)算模型,將污染因子權(quán)重最大值對應(yīng)的 Fj,w直接作為第三項(xiàng)參數(shù),F(xiàn)j為 Fj,max、Fj,ave和 Fj,w3項(xiàng)的幾何平均值??傮w而言,Ⅱ型修正內(nèi)梅羅指數(shù)法是對典型內(nèi)梅羅指數(shù)法的大膽突破和創(chuàng)新,較大地?cái)U(kuò)展了方法的適用性,提高了評價結(jié)果的客觀性和準(zhǔn)確性。對于參評污染因子眾多的環(huán)境要素的綜合性質(zhì)量評價,修正意義尤其重大。同時,該法也是對Ⅰ型修正內(nèi)梅羅指數(shù)法的拓展和延伸,對于評價指向性和目的性較強(qiáng)的環(huán)境要素質(zhì)量評價同樣適用。

    表7 Ⅱ型修正內(nèi)梅羅指數(shù)法的修正原則Table 7 Modifications ofⅡ-RNI method

    3 結(jié)論

    (1)經(jīng)數(shù)學(xué)推導(dǎo),典型內(nèi)梅羅指數(shù)法計(jì)算公式可變換為線性表達(dá)式:Fj=aFj,max+bFj,ave。其中,斜率 a 和截距b分別代表單因子污染指數(shù)最大值(Fj,max)和算術(shù)平均值(Fj,ave)的權(quán)重系數(shù)。a 和 b 的大小受 Fj,max與Fj,ave比值的影響顯著,隨著比值的增大,F(xiàn)j,max權(quán)重增大,F(xiàn)j,ave權(quán)重減小。參評污染因子數(shù)量(n)對典型內(nèi)梅羅指數(shù)法的評分結(jié)果存在著重要影響,n越大,典型內(nèi)梅羅指數(shù)法對環(huán)境要素質(zhì)量評價不適用的風(fēng)險(xiǎn)越高。

    (2)典型內(nèi)梅羅指數(shù)法更適用于區(qū)域中目的性和指向性較明確或者總體污染程度較輕的環(huán)境要素質(zhì)量評價。

    (3)為提高典型內(nèi)梅羅指數(shù)法的適用性和評價結(jié)果的準(zhǔn)確性,確定了2種修正內(nèi)梅羅指數(shù)法的劃分原則。Ⅰ型修正內(nèi)梅羅指數(shù)法以n等于5作為劃分結(jié)點(diǎn),設(shè)置了1種修正方式,該法更加傾向于參評因子不多、評價目的性和指向性較強(qiáng)的環(huán)境要素質(zhì)量評價,如環(huán)境要素中某類特定類別的污染物(重金屬、多環(huán)芳烴等)的環(huán)境污染風(fēng)險(xiǎn)評價;Ⅱ型修正內(nèi)梅羅指數(shù)法是對Ⅰ型修正內(nèi)梅羅指數(shù)法的拓展和延伸,以Fj,max與 Fj,ave比值等于 5 為劃分結(jié)點(diǎn),修正方式包括Fj,max降權(quán)、Fj,ave加權(quán)、Fj,max降權(quán)同時 Fj,ave加權(quán)以及改變計(jì)算模型等4種,對多因子綜合性質(zhì)量評價更具適用性。

    (4)在 n≤5 或 Fj,max/Fj,ave≤5 時,典型內(nèi)梅羅指數(shù)法可不修正,不影響評價結(jié)果的客觀性和準(zhǔn)確性。

    [1]陳京都,戴其根,許學(xué)宏,等.江蘇省典型區(qū)農(nóng)田土壤及小麥中重金屬含量與評價[J].生態(tài)學(xué)報(bào),2012,32(11):3487-3496.CHEN Jing-du,DAI Qi-gen,XU Xue-hong,et al.Heavy metal contents and evaluation of farmland soil and wheat in typical area of Jiangsu Province[J].Acta Ecologica Sinica,2012,32(11):3487-3496.

    [2]楊 琳.應(yīng)用權(quán)重修正內(nèi)梅羅污染指數(shù)法對陸源入海排污口鄰近海域水環(huán)境質(zhì)量綜合評價[J].福建水產(chǎn),2010(2):50-54.YANG Lin.A comprehensive evaluation on coastal waters near pollutant discharge outlet with weighting-factor modified Nemerow index method[J].Journal of Fujian Fisheries,2010(2):50-54.

    [3]丁 菁.內(nèi)梅羅污染指數(shù)法在排污口鄰近海域水環(huán)境質(zhì)量評價中的應(yīng)用[J].福建水產(chǎn),2016(1):1-4.DING Jing.An application of the Nemerow pollution index in the quality assessment of the seawater environment in the vicinity of terrestrial drains[J].Journal of Fujian Fisheries,2016(1):1-4.

    [4]羅萬次,蘇 搏,劉 熊,等.北侖河口附近海域冬季海洋環(huán)境質(zhì)量評價[J].廣西科學(xué)院學(xué)報(bào),2014,30(2):107-111,119.LUO Wan-ci,SU Bo,LIU Xiong,et al.Winter marine environment quality assessment around the Beilun estuary[J].Journal of Guangxi A-cademy of Sciences,2014,30(2):107-111,119.

    [5]宋永剛,吳金浩,邵澤偉,等.遼東灣近岸表層海水重金屬污染分析與評價[J].漁業(yè)科學(xué)進(jìn)展,2016,37(3):14-19.SONG Yong-gang,WU Jin-hao,SHAO Ze-wei,et al.Evaluation of heavy metal pollution in the offshore surface seawater of the Liaodong Bay[J].Progress in Fishery Sciences,2016,37(3):14-19.

    [6]Yi Q,Dou X D,Huang Q R,et al.Pollution characteristics of Pb,Zn,As,Cd in the Bijiang River[J].Procedia Environmental Sciences,2012,13(13):43-52.

    [7]梁德華,蔣火華.河流水質(zhì)綜合評價方法的統(tǒng)一和改進(jìn)[J].中國環(huán)境監(jiān)測,2002,18(2):63-66.LIANG De-hua,JIANG Huo-hua.Unifying and improving the comprehensive assessment methods of river water quality[J].Environmental Monitoring in China,2002,18(2):63-66.

    [8]蔡 曄,林怡雯,李月娥,等.利用改進(jìn)的內(nèi)梅羅指數(shù)法模型評價蘇州市內(nèi)外城河水質(zhì)[J].化學(xué)分析計(jì)量,2015,24(2):84-87.CAI Ye,LIN Yi-wen,LI Yue-e,et al.Evaluation of water quality of city river in Suzhou by improved Nemerow index method[J].Chemical Analysis and Meterage,2015,24(2):84-87.

    [9]李小麗,黎小東,敖天其.改進(jìn)內(nèi)梅羅指數(shù)法在西充河水質(zhì)評價中的應(yīng)用[J].人民黃河,2016,38(8):65-68.LI Xiao-li,LI Xiao-dong,AO Tian-qi.Application of improved Nemerow index method used in water quality assessment[J].Yellow River,2016,38(8):65-68.

    [10]鄭利杰,高紅杰,宋永會,等.我國典型城市地表水水質(zhì)綜合評價與分析[J].環(huán)境工程技術(shù)學(xué)報(bào),2016,6(3):252-258.ZHENG Li-jie,GAO Hong-jie,SONG Yong-hui,et al.Comprehensive evaluation and analysis of surface water quality for typical cities of Chi-na[J].Journal of Environmental Engineering Technology,2016,6(3):252-258.

    [11]周笑白,張寧紅,張 詠,等.太湖水質(zhì)與水生生物健康的關(guān)聯(lián)性初探[J].環(huán)境科學(xué),2014,35(1):271-278.ZHOU Xiao-bai,ZHANG Ning-hong,ZHANG Yong,et al.Preliminary study on the relationship between the water quality and the aquatic biological health status of Taihu Lake[J].Environmental Science,2014,35(1):271-278.

    [12]蔡曉強(qiáng),高強(qiáng)立.基于內(nèi)梅羅指數(shù)法評價開封市湖泊水質(zhì)污染特征[J].吉林農(nóng)業(yè),2016(11):82-83.CAI Xiao-qiang,GAO Qiang-li.Pollution characteristics of water quality in Kaifeng based on Nemerow index method[J].Jilin Agriculture,2016(11):82-83.

    [13]孫清展,臧淑英,張囡囡,等.基于蒙特卡羅方法的扎龍濕地水環(huán)境質(zhì)量評價[J].濕地科學(xué),2013,11(1):75-81.SUN Qing-zhan,ZANG Shu-ying,ZHANG Nan-nan,et al.Evaluation of water environment quality of Zhalong wetlands based on Monte Carlo method[J].Wetland Science,2013,11(1):75-81.

    [14]呂俊芳,王書錦,劉云根,等.羅時江河口濕地水質(zhì)季節(jié)變化特征及評價[J].環(huán)境工程,2016(S1):845-849,877.Lü Jun-fang,WANG Shu-jin,LIU Yun-gen,et al.Seasonal variation characteristics of water quality and evaluation at the estuarine wetland of Luoshijiang River[J].Environmental Engineering,2016(S1):845-849,877.

    [15]何增輝.修正內(nèi)梅羅污染指數(shù)法在水源地環(huán)境質(zhì)量評價中的應(yīng)用[J].廣東化工,2011,38(7):141-143.HE Zeng-hui.Application of improved Nemerow exponent method for evaluation of environmental quality of waterhead area[J].Guangdong Chemical Engineering,2011,38(7):141-143.

    [16]丁雪卿.改進(jìn)的內(nèi)梅羅污染指數(shù)法在集中式飲用水源地環(huán)境質(zhì)量評價中的應(yīng)用[J].四川環(huán)境,2010,29(2):47-51.DING Xue-qing.Application of improved Nemerow index method for evaluation of environmental quality of waterhead area of public drinking water supply[J].Sichuan Environment,2010,29(2):47-51.

    [17]寇文杰,林 健,陳忠榮,等.內(nèi)梅羅指數(shù)法在水質(zhì)評價中存在問題及修正[J].南水北調(diào)與水利科技,2012,10(4):39-41,47.KOU Wen-jie,LIN Jian,CHEN Zhong-rong,et al.Existing problems and modifications of using Nemerow index method in water quality assessment[J].South-to-North Water Diversion and Water Science&Technology,2012,10(4):39-41,47.

    [18]谷朝君,潘 穎,潘明杰.內(nèi)梅羅指數(shù)法在地下水水質(zhì)評價中的應(yīng)用及存在問題[J].環(huán)境保護(hù)科學(xué),2002,28(1):45-47.GU Zhao-jun,PAN Ying,PAN Ming-jie.The application and existed problems of Nemero index in groundwater quality evaluation[J].Environmental Protection Science,2002,28(1):45-47.

    [19]馬成有,曹劍鋒,姜紀(jì)沂,等.改進(jìn)的內(nèi)梅羅污染指數(shù)法及其應(yīng)用——以磐石市地下水環(huán)境質(zhì)量評價為例[J].水資源保護(hù),2006,22(4):53-55.MA Cheng-you,CAO Jian-feng,JIANG Ji-yi,et al.Application of improved Nemerow index method to evaluation of groundwater quality in Panshi City[J].Water Resources Protection,2006,22(4):53-55.

    [20]關(guān)云鵬.利用內(nèi)梅羅指數(shù)法模型評價地下水水質(zhì)的探討[J].山西水利科技,2012(1):81-84.GUAN Yun-peng.A discussion on evaluating groundwater quality by Nemerow index method[J].Shanxi Hydrotechnics,2012(1):81-84.

    [21]張安昌,莊付磊,申銅飛,等.改進(jìn)內(nèi)梅羅指數(shù)法在聊城市地下水水質(zhì)評價中的應(yīng)用[J].治淮,2016(12):81-83.ZHANG An-chang,ZHUANG Fu-lei,SHEN Tong-fei,et al.Application of improved Nemerow index method to evaluation of groundwater quality in Liaocheng City[J].Zhihuai,2016(12):81-83.

    [22]于福榮,盧文喜,卞玉梅,等.改進(jìn)的尼梅羅污染指數(shù)法在黃龍工業(yè)園水質(zhì)評價中的應(yīng)用[J].世界地質(zhì),2008,27(1):59-62.YU Fu-rong,LU Wen-xi,BIAN Yu-mei,et al.Application of improved Nemerow index method in evaluation of water quality in Huanglong industrial park catchment[J].Global Geology,2008,27(1):59-62.

    [23]李亞松,張兆吉,費(fèi)宇紅,等.內(nèi)梅羅指數(shù)評價法的修正及其應(yīng)用[J].水資源保護(hù),2009,25(6):48-50.LI Ya-song,ZHANG Zhao-ji,FEI Yu-hong,et al.Improvement of Nemerow index method and its application[J].Water Resources Protection,2009,25(6):48-50.

    [24]Jie C,Qing L,Hui Q.Application of improved Nemerow index method based on entropy weight for groundwater quality evaluation[J].International Journal of Environmental Sciences,2012,2(3):1284-1290.

    [25]李玲玲.內(nèi)梅羅指數(shù)法的修正及應(yīng)用[J].地下水,2015,37(5):101-102,107.LI Ling-ling.Improvement and application of Nemerow index method[J].Groundwater,2015,37(5):101-102,107.

    [26]閆欣榮.修正的內(nèi)梅羅指數(shù)法及其在城市地下飲用水源地水質(zhì)評價中的應(yīng)用[J].地下水,2010,32(1):6-7.YAN Xin-rong.Modified Nemerow exponent and application of the water quality evaluation in urban underground drinking water source area[J].Groundwater,2010,32(1):6-7.

    [27]章熙鋒,唐家良,薛 菲,等.5·12汶川大地震后涪江綿陽段沉積物重金屬污染評價[J].環(huán)境保護(hù)科學(xué),2015,41(5):142-146.ZHANG Xi-feng,TANG Jia-liang,XUE Fei,et al.Evaluation of heavy metal pollution in the sediments from Mianyang section of Fu River after the earthquake in Wenchuan on May 12[J].Environmental Protection Science,2015,41(5):142-146.

    [28]鄭懿珉,高茂生,劉 森,等.萊州灣表層沉積物重金屬分布特征及生態(tài)環(huán)境評價[J].海洋環(huán)境科學(xué),2015,34(3):354-360.ZHENG Yi-min,GAO Mao-sheng,LIU Sen,et al.Distribution patterns and ecological assessment on heavy metals in the surface sediments of Laizhou Bay[J].Marine Environmental Science,2015,34(3):354-360.

    [29]Brady J P,Ayoko G A,Martens W N,et al.Enrichment,distribution and sources of heavy metals in the sediments of Deception Bay,Queensland,Australia[J].Marine Pollution Bulletin,2014,81(1):248-255.

    [30]Liu S,Zhang Y,Bi S,et al.Heavy metals distribution and environmental quality assessment for sediments off the southern coast of the Shandong Peninsula,China[J].Marine Pollution Bulletin,2015,100(1):483-488.

    [31]Kowalska J,Mazurek R,G?siorek M,et al.Soil pollution indices conditioned by medieval metallurgical activity:A case study from Krakow(Poland)[J].Environmental Pollution,2016,218:1023-1036.

    [32]鄒 喬,王 瑤,杜顯元,等.內(nèi)梅羅指數(shù)法在石油開采區(qū)土壤PAHs污染分級評價中的應(yīng)用 [J].安徽農(nóng)業(yè)科學(xué),2011,39(12):7350-7353.ZOU Qiao,WANG Yao,DU Xian-yuan,et al.Application of Nemerow index method in grade evaluation of PAHs pollution in the oilfield soils[J].Journal of Anhui Agricultural Science,2011,39(12):7350-7353.

    [33]安 婧,宮曉雙,陳宏偉,等.沈撫灌區(qū)農(nóng)田土壤重金屬污染時空變化特征及生態(tài)健康風(fēng)險(xiǎn)評價[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2016,35(1):37-44.AN Jing,GONG Xiao-shuang,CHEN Hong-wei,et al.Temporal and spatial characteristics and health risk assessments of heavy metal pollution in soils of Shenfu irrigation area[J].Journal of Agro-Environment Science,2016,35(1):37-44.

    [34]劉亞納,朱書法,魏學(xué)鋒,等.河南洛陽市不同功能區(qū)土壤重金屬污染特征及評價[J].環(huán)境科學(xué),2016,37(6):2322-2328.LIU Ya-na,ZHU Shu-fa,WEI Xue-feng,et al.Assessment and pollution characteristics of heavy metals in soil of different functional areas in Luoyang[J].Environmental Science,2016,37(6):2322-2328.

    [35]楊玉敏,師學(xué)義,張 琛.基于內(nèi)梅羅指數(shù)法的復(fù)墾村莊土壤重金屬污染評價及空間分布[J].水土保持研究,2016,23(4):338-343.YANG Yu-min,SHI Xue-yi,ZHANG Chen.Spatial distribution and evaluation of heavy metal pollution of reclaiming village based on Nemerow integrated pollution index method[J].Research of Soil and Water Conservation,2016,23(4):338-343.

    [36]張金婷,孫 華.內(nèi)梅羅指數(shù)法和模糊綜合評價法在土壤重金屬污染評價應(yīng)用中的差異分析 [J].環(huán)境監(jiān)測管理與技術(shù),2016,28(4):27-31.ZHANG Jin-ting,SUN Hua.Differences of Nemerow index method and fuzzy comprehensive evaluation method in evaluation heavy metal pollution in soil[J].The Administration and Technique of Environmental Monitoring,2016,28(4):27-31.

    [37]李小曼,徐夢潔,劉 勤,等.基于內(nèi)梅羅指數(shù)法及其改進(jìn)方法的小尺度區(qū)域土壤重金屬污染評價:以蘇南地區(qū)為例[J].江蘇農(nóng)業(yè)科學(xué),2017,45(3):241-245.LI Xiao-man,XU Meng-jie,LIU Qin,et al.Evaluation of heavy metal pollution of soil in small scale region based on typical and improved Nemerow index method in the case of South of Jiangsu[J].Journal of Jiangsu Agricultural Science,2017,45(3):241-245.

    [38]林 勇.濟(jì)南市可耕土壤有機(jī)氯、多氯聯(lián)苯污染分析[D].濟(jì)南:山東大學(xué),2011:40-41.LIN Yong.Pollution analysis of OCPs and PCBs from soil of Ji′nan[D].Ji′nan:Shandong University,2011:40-41.

    [39]楊賢智,李景錕,廖延梅.環(huán)境管理學(xué)[M].北京:高等教育出版社,1990:185.YANG Xian-zhi,LI Jing-kun,LIAO Yan-mei.Environmental management[M].Beijing:Higher Education Press,1990:185.

    [40]Swaine D J.Why trace elements are important?[J].Fuel Processing Technology,2000,s65-66(1):21-33.

    [41]徐爭啟,倪師軍,庹先國,等.潛在生態(tài)危害指數(shù)法評價中重金屬毒性系數(shù)計(jì)算[J].環(huán)境科學(xué)與技術(shù),2008,31(2):112-115.XU Zheng-qi,NI Shi-jun,TUO Xian-guo,et al.Calculation of heavy metals′toxicity coefficient in the evaluation of potential ecological risk index[J].Environmental Science&Technology,2008,31(2):112-115.

    [42]陳靜生,王 忠,劉玉機(jī).水體金屬污染潛在危害:應(yīng)用沉積學(xué)方法評價[J].環(huán)境科技(遼寧),1989,9(1):16-25.CHEN Jing-sheng,WANG Zhong,LIU Yu-ji.Potential risk of heavy metal pollution in water:Assessment by sedimentology measure[J].Environmental Science and Technology(Liaoning),1989,9(1):16-25.

    [43]林麗欽.應(yīng)用毒理學(xué)安全評價數(shù)據(jù)推算重金屬毒性系數(shù)的探討[C]//2009重金屬污染監(jiān)測、風(fēng)險(xiǎn)評價及修復(fù)技術(shù)高級研討會論文集.北京:中華環(huán)保聯(lián)合會能源環(huán)境專業(yè)委員會,2009:57-60.LIN Li-qin.Discussion on calculation of heavy metals’toxicity coefficient based on toxicological safety evaluation data[C]//Proceedings of advanced seminar on 2009 heavy metal pollution monitoring,risk assessment and remediation technology.Beijing:Energy Resources and Environment Professional Council of All-China Environment Federation,2009:57-60.

    [44]李雪梅,王祖?zhèn)?湯顯強(qiáng),等.重金屬污染因子權(quán)重的確定及其在土壤環(huán)境質(zhì)量評價中的應(yīng)用[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2007,26(6):2281-2286.LI Xue-mei,WANG Zu-wei,TANG Xian-qiang,et al.Determining weights of heavy metal contaminations and its application to soil environmental quality assessment[J].Journal of Agro-Environment Science,2007,26(6):2281-2286.

    Applicability and modifications of the Nemerow index method in evaluating environmental quality

    HAN Shu-xin,WANG Li-hong,ZHAO Chang-sheng
    (Shandong Analysis and Test Center,Ji′nan 250014,China)

    In order to study the applicability and modifications of the Nemerow index(NI)method used for evaluations of environmental quality,mathematical derivations were performed based on the original formula.The results showed that the formula of the NI method could be transformed into a linear equation.The slope and intercept of the equation could be the weight coefficients of the maximum value(Fj,max)and arithmetic mean value(Fj,ave),respectively.With increases in the ratio of Fj,maxto Fj,ave,the weight of Fj,maxincreased while that of Fj,avedecreased.The number of pollution factors(n)participating in the evaluation had a significant influence on the results.The NI method was found to be more suitable for evaluation of environmental factors with clear purposes and directions,or those associated with low pollution in a region.In order to improve the applicability of the NI method,some modifications were proposed for two revised Nemerow index(RNI)methods.For theⅠ-RNI method,n=5 was used as the partition node,and the ratio of Fj,maxto Fj,avewas equal to five for theⅡ-RNI method.TheⅠ-RNI method was used to evaluate factors with clear purposes and directions,and those associated with less pollution;theⅡ-RNI method,an extension of theⅠ-RNI method,was more suitable for multifactor comprehensive evaluation.Based on the mathematical connotation of the linear equation for the NI method,when n or the ratio of Fj,maxto Fj,aveis not more than 5,the original NI method can be used without affecting the objectivity and accuracy of the evaluation results.

    Nemerow index method;environmental factors;applicability of methods;modifications of methods;weights of pollution factors

    X820.2

    A

    1672-2043(2017)10-2153-08

    10.11654/jaes.2017-0519

    韓術(shù)鑫,王利紅,趙長盛,等.內(nèi)梅羅指數(shù)法在環(huán)境質(zhì)量評價中的適用性與修正原則[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2017,36(10):2153-2160.

    HAN Shu-xin,WANG Li-hong,ZHAO Chang-sheng.Applicability and modifications of the Nemerow index method in evaluating environmental quality[J].Journal of Agro-Environment Science,2017,36(10):2153-2160.

    2017-04-10 錄用日期:2017-06-15

    韓術(shù)鑫(1984—),男,山東諸城人,助理研究員,主要研究方向?yàn)榄h(huán)境監(jiān)測與環(huán)境質(zhì)量評價。E-mail:hsxaa369@qq.com

    山東省重點(diǎn)研發(fā)計(jì)劃(2015GSF120010)

    Project supported:The Key Technology Research and Development Program of Shandong Province,China(2015GSF120010)

    猜你喜歡
    內(nèi)梅羅環(huán)境要素適用性
    環(huán)境要素與住區(qū)戶外活動影響機(jī)制研究
    住區(qū)(2023年4期)2023-12-01 08:05:10
    長江經(jīng)濟(jì)帶水質(zhì)監(jiān)測斷面嘉興段評價分析
    綠色科技(2022年10期)2022-06-23 10:14:26
    強(qiáng)調(diào)簡潔和適用性 MICHI by Rotel X5/X3合并功放
    環(huán)保技術(shù)在土木工程領(lǐng)域的適用性探討
    小議階級分析理論的適用性
    一種新內(nèi)梅羅指數(shù)法的修正方法及應(yīng)用
    內(nèi)梅羅污染指數(shù)法和單因子評價法在水質(zhì)評價中的應(yīng)用
    新一代車上轉(zhuǎn)換裝置在廠礦企業(yè)鐵路運(yùn)輸中的適用性研究
    湛江東海島潮間帶表層沉積物粒度的分布及與環(huán)境要素的相關(guān)性
    用內(nèi)梅羅指數(shù)法評價武陟縣境內(nèi)沁河水質(zhì)
    夜夜看夜夜爽夜夜摸| 国产成人精品无人区| 日韩中文字幕欧美一区二区| 成年人黄色毛片网站| 好男人电影高清在线观看| 日韩欧美国产在线观看| 日本精品一区二区三区蜜桃| 又黄又爽又免费观看的视频| 热99re8久久精品国产| 99精品欧美一区二区三区四区| 在线十欧美十亚洲十日本专区| 夜夜夜夜夜久久久久| 中亚洲国语对白在线视频| 每晚都被弄得嗷嗷叫到高潮| 制服人妻中文乱码| 91成年电影在线观看| 不卡一级毛片| 最近最新中文字幕大全电影3 | 男女午夜视频在线观看| 国产激情偷乱视频一区二区| 男人的好看免费观看在线视频 | 欧美乱码精品一区二区三区| 夜夜夜夜夜久久久久| 国产精品av久久久久免费| 999久久久国产精品视频| 亚洲国产欧美网| 精品久久蜜臀av无| 亚洲国产精品成人综合色| 欧美黑人精品巨大| 亚洲一区二区三区不卡视频| 十分钟在线观看高清视频www| 欧美激情 高清一区二区三区| 变态另类成人亚洲欧美熟女| 国产成人啪精品午夜网站| 成在线人永久免费视频| 午夜视频精品福利| 欧美日本视频| 亚洲精品国产一区二区精华液| 成人午夜高清在线视频 | 国产欧美日韩精品亚洲av| 色婷婷久久久亚洲欧美| 两人在一起打扑克的视频| 色播亚洲综合网| 天堂√8在线中文| 成人三级黄色视频| 国产一区二区三区视频了| netflix在线观看网站| 很黄的视频免费| 国产欧美日韩一区二区三| 亚洲午夜精品一区,二区,三区| 最新在线观看一区二区三区| 成人免费观看视频高清| 97碰自拍视频| 精品国产乱子伦一区二区三区| 国产国语露脸激情在线看| 欧美激情极品国产一区二区三区| 日韩精品青青久久久久久| 午夜久久久久精精品| 最近最新中文字幕大全免费视频| 女人高潮潮喷娇喘18禁视频| 高潮久久久久久久久久久不卡| 亚洲人成电影免费在线| 亚洲av日韩精品久久久久久密| 老司机靠b影院| 日日夜夜操网爽| 日韩高清综合在线| 啦啦啦免费观看视频1| 国产一级毛片七仙女欲春2 | 一二三四在线观看免费中文在| 成年女人毛片免费观看观看9| 黄频高清免费视频| 少妇被粗大的猛进出69影院| 久久香蕉精品热| 欧美大码av| 51午夜福利影视在线观看| 国产又爽黄色视频| 亚洲自偷自拍图片 自拍| 色综合欧美亚洲国产小说| 露出奶头的视频| 国产一区二区三区视频了| 夜夜躁狠狠躁天天躁| 90打野战视频偷拍视频| 亚洲中文字幕一区二区三区有码在线看 | 波多野结衣av一区二区av| 最近最新中文字幕大全电影3 | 国产1区2区3区精品| 制服诱惑二区| www.999成人在线观看| 成人国产一区最新在线观看| 老熟妇乱子伦视频在线观看| 最近最新中文字幕大全电影3 | 亚洲免费av在线视频| 少妇裸体淫交视频免费看高清 | √禁漫天堂资源中文www| 黄片大片在线免费观看| 久久久久久国产a免费观看| 97人妻精品一区二区三区麻豆 | 此物有八面人人有两片| 欧美乱色亚洲激情| 午夜久久久久精精品| 亚洲色图 男人天堂 中文字幕| 嫩草影院精品99| 久久婷婷人人爽人人干人人爱| 亚洲国产欧美日韩在线播放| 国产免费男女视频| 国产99久久九九免费精品| 最新在线观看一区二区三区| 免费在线观看黄色视频的| 国产成人精品久久二区二区91| 88av欧美| 黄色女人牲交| 午夜视频精品福利| 亚洲av五月六月丁香网| 一级毛片高清免费大全| 国产精品一区二区免费欧美| 亚洲精品国产一区二区精华液| 精品一区二区三区av网在线观看| 中文资源天堂在线| 少妇 在线观看| 日韩成人在线观看一区二区三区| 国产成年人精品一区二区| 久久人妻福利社区极品人妻图片| 国产精华一区二区三区| 午夜福利成人在线免费观看| 91在线观看av| 999久久久精品免费观看国产| 18禁国产床啪视频网站| 99热6这里只有精品| 欧美乱色亚洲激情| 久久精品夜夜夜夜夜久久蜜豆 | 大香蕉久久成人网| 女人爽到高潮嗷嗷叫在线视频| 啦啦啦免费观看视频1| 亚洲精华国产精华精| 91大片在线观看| 亚洲中文字幕一区二区三区有码在线看 | 亚洲人成网站在线播放欧美日韩| 老汉色av国产亚洲站长工具| 精品日产1卡2卡| 亚洲五月色婷婷综合| 亚洲精品色激情综合| 日本一区二区免费在线视频| 亚洲av五月六月丁香网| 亚洲精品美女久久av网站| 无人区码免费观看不卡| 亚洲午夜精品一区,二区,三区| 宅男免费午夜| 在线av久久热| 精品国产亚洲在线| 精品免费久久久久久久清纯| 日韩欧美国产在线观看| 99久久无色码亚洲精品果冻| 久久精品国产亚洲av高清一级| 国产真实乱freesex| 法律面前人人平等表现在哪些方面| 中文字幕精品免费在线观看视频| 久久天堂一区二区三区四区| 一级作爱视频免费观看| 国产伦在线观看视频一区| 美女大奶头视频| 大型av网站在线播放| 精品久久久久久成人av| 久久久久久免费高清国产稀缺| 天堂√8在线中文| 一进一出抽搐动态| 少妇熟女aⅴ在线视频| 国产精品香港三级国产av潘金莲| 制服诱惑二区| 老司机深夜福利视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 又大又爽又粗| 亚洲中文日韩欧美视频| 久久国产精品男人的天堂亚洲| 色播亚洲综合网| 一本精品99久久精品77| 91国产中文字幕| 50天的宝宝边吃奶边哭怎么回事| 国产野战对白在线观看| 亚洲av熟女| 一本精品99久久精品77| 午夜免费成人在线视频| 最近最新免费中文字幕在线| 丁香六月欧美| av欧美777| 国产成人影院久久av| 我的亚洲天堂| 丁香欧美五月| 成年人黄色毛片网站| 波多野结衣巨乳人妻| 村上凉子中文字幕在线| 在线观看午夜福利视频| 国产成人精品久久二区二区免费| av欧美777| 香蕉国产在线看| 真人做人爱边吃奶动态| 午夜福利一区二区在线看| 欧美乱妇无乱码| 精品卡一卡二卡四卡免费| 美女免费视频网站| 久久久国产成人免费| 亚洲国产日韩欧美精品在线观看 | 亚洲精品久久国产高清桃花| 午夜日韩欧美国产| 国产私拍福利视频在线观看| 99久久精品国产亚洲精品| 亚洲国产精品合色在线| 99在线人妻在线中文字幕| 最新在线观看一区二区三区| 不卡一级毛片| 欧美zozozo另类| 亚洲人成网站高清观看| 亚洲国产毛片av蜜桃av| 1024手机看黄色片| 婷婷亚洲欧美| 免费高清在线观看日韩| 级片在线观看| 亚洲人成伊人成综合网2020| 婷婷精品国产亚洲av| 一二三四社区在线视频社区8| e午夜精品久久久久久久| 欧美日韩精品网址| 精品久久久久久久毛片微露脸| 一卡2卡三卡四卡精品乱码亚洲| 熟女电影av网| 香蕉丝袜av| 国产高清videossex| 亚洲 国产 在线| 51午夜福利影视在线观看| 天堂影院成人在线观看| 国产精品久久久av美女十八| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美精品亚洲一区二区| 国产精品精品国产色婷婷| 日本在线视频免费播放| 女人高潮潮喷娇喘18禁视频| 俺也久久电影网| 免费观看精品视频网站| 免费在线观看影片大全网站| 亚洲av熟女| 久久青草综合色| 欧美在线黄色| 三级毛片av免费| 一二三四在线观看免费中文在| 亚洲av日韩精品久久久久久密| 亚洲国产中文字幕在线视频| 久久人人精品亚洲av| 久久久久国内视频| 老司机深夜福利视频在线观看| 亚洲成av片中文字幕在线观看| 国产单亲对白刺激| 亚洲 国产 在线| 精品久久久久久久末码| 老司机深夜福利视频在线观看| 亚洲国产欧洲综合997久久, | av欧美777| 国产精品美女特级片免费视频播放器 | 亚洲成人久久爱视频| 欧美日韩精品网址| 男女之事视频高清在线观看| 久久国产亚洲av麻豆专区| 午夜a级毛片| 国产成人一区二区三区免费视频网站| 国产亚洲精品久久久久久毛片| 黄色视频,在线免费观看| 精品第一国产精品| 午夜a级毛片| 非洲黑人性xxxx精品又粗又长| 女同久久另类99精品国产91| 国产成人精品无人区| 一级毛片高清免费大全| 亚洲国产欧洲综合997久久, | 日本精品一区二区三区蜜桃| 在线播放国产精品三级| 精品国产一区二区三区四区第35| 免费在线观看亚洲国产| 欧美亚洲日本最大视频资源| 18禁裸乳无遮挡免费网站照片 | 黄色视频不卡| 最新美女视频免费是黄的| 波多野结衣av一区二区av| 久久久国产成人免费| 国产精品1区2区在线观看.| 成在线人永久免费视频| 真人做人爱边吃奶动态| 长腿黑丝高跟| 观看免费一级毛片| 少妇被粗大的猛进出69影院| 欧美日韩亚洲国产一区二区在线观看| 久久久国产欧美日韩av| 免费在线观看成人毛片| 男人舔女人的私密视频| 一本一本综合久久| 麻豆成人av在线观看| 性色av乱码一区二区三区2| 丁香六月欧美| 哪里可以看免费的av片| 久久久久久久久免费视频了| 国产av一区二区精品久久| 成年免费大片在线观看| 神马国产精品三级电影在线观看 | 国产精品一区二区免费欧美| 婷婷六月久久综合丁香| 午夜福利成人在线免费观看| 搡老岳熟女国产| 久久久精品国产亚洲av高清涩受| 在线观看www视频免费| av视频在线观看入口| 亚洲午夜理论影院| 欧美日韩中文字幕国产精品一区二区三区| 精品久久久久久久久久免费视频| 国产主播在线观看一区二区| 国产亚洲欧美精品永久| 亚洲国产毛片av蜜桃av| 首页视频小说图片口味搜索| 久久久久亚洲av毛片大全| 1024手机看黄色片| 国产一区二区三区视频了| a级毛片在线看网站| 女人爽到高潮嗷嗷叫在线视频| 国产av不卡久久| 欧美大码av| 国产高清激情床上av| 很黄的视频免费| 老司机靠b影院| 一本久久中文字幕| 国产亚洲精品久久久久久毛片| 成人午夜高清在线视频 | 成人午夜高清在线视频 | 亚洲国产精品成人综合色| 少妇的丰满在线观看| 999久久久国产精品视频| 国产亚洲av嫩草精品影院| 色老头精品视频在线观看| 极品教师在线免费播放| 欧美乱码精品一区二区三区| 成在线人永久免费视频| netflix在线观看网站| 欧美激情极品国产一区二区三区| 国产精品电影一区二区三区| 哪里可以看免费的av片| 啦啦啦韩国在线观看视频| 18禁国产床啪视频网站| 国产又黄又爽又无遮挡在线| 日韩欧美一区二区三区在线观看| 国产单亲对白刺激| 国产精品久久久av美女十八| 中出人妻视频一区二区| www日本黄色视频网| 欧美丝袜亚洲另类 | 欧美乱色亚洲激情| 制服丝袜大香蕉在线| 色播在线永久视频| 成在线人永久免费视频| 香蕉久久夜色| 亚洲九九香蕉| 一本综合久久免费| 国产伦人伦偷精品视频| 亚洲最大成人中文| 人成视频在线观看免费观看| 老熟妇乱子伦视频在线观看| 99国产精品一区二区蜜桃av| 精品福利观看| 99国产综合亚洲精品| 亚洲av片天天在线观看| 黄色视频,在线免费观看| 国产成人精品久久二区二区免费| 精品久久久久久久毛片微露脸| 国产1区2区3区精品| 日韩三级视频一区二区三区| 99精品在免费线老司机午夜| 啦啦啦 在线观看视频| 亚洲欧美日韩高清在线视频| 变态另类丝袜制服| 成人精品一区二区免费| 精品久久久久久久毛片微露脸| 一级毛片高清免费大全| 日韩视频一区二区在线观看| 久久精品国产综合久久久| 久久国产亚洲av麻豆专区| 黄色视频不卡| 色综合婷婷激情| 成人亚洲精品一区在线观看| 亚洲av成人一区二区三| 精品国产乱码久久久久久男人| 国产熟女xx| or卡值多少钱| 欧美乱色亚洲激情| 亚洲av熟女| 曰老女人黄片| av超薄肉色丝袜交足视频| 亚洲午夜精品一区,二区,三区| 满18在线观看网站| 一级a爱视频在线免费观看| av视频在线观看入口| 日韩精品中文字幕看吧| a级毛片在线看网站| 国产精品1区2区在线观看.| 亚洲五月天丁香| 国产精品1区2区在线观看.| 中亚洲国语对白在线视频| 久久国产精品影院| 久久午夜亚洲精品久久| 亚洲专区中文字幕在线| 两人在一起打扑克的视频| 久久久国产欧美日韩av| 欧美日韩亚洲国产一区二区在线观看| 少妇被粗大的猛进出69影院| 国产精品一区二区免费欧美| 国产黄a三级三级三级人| 久久 成人 亚洲| 亚洲精品国产精品久久久不卡| 久久 成人 亚洲| 90打野战视频偷拍视频| 一本大道久久a久久精品| 亚洲七黄色美女视频| 亚洲va日本ⅴa欧美va伊人久久| 国产高清有码在线观看视频 | 欧美日韩黄片免| 午夜精品久久久久久毛片777| 午夜福利高清视频| 黄色片一级片一级黄色片| 久久欧美精品欧美久久欧美| 国产单亲对白刺激| 久久亚洲精品不卡| 黑人欧美特级aaaaaa片| 变态另类丝袜制服| 韩国av一区二区三区四区| 欧美日韩乱码在线| av免费在线观看网站| 午夜免费激情av| 黄色 视频免费看| 黄色片一级片一级黄色片| 亚洲真实伦在线观看| 免费看日本二区| avwww免费| 正在播放国产对白刺激| 国产av不卡久久| 亚洲av片天天在线观看| 欧美成人性av电影在线观看| 精品久久久久久久久久久久久 | 久久久国产成人精品二区| 18禁黄网站禁片免费观看直播| 美女扒开内裤让男人捅视频| 久久天堂一区二区三区四区| 国产激情欧美一区二区| 极品教师在线免费播放| 老熟妇乱子伦视频在线观看| 国产精华一区二区三区| 嫩草影视91久久| 国产精品亚洲av一区麻豆| 男人舔女人下体高潮全视频| 91大片在线观看| АⅤ资源中文在线天堂| 亚洲国产日韩欧美精品在线观看 | 日本黄色视频三级网站网址| 一进一出抽搐动态| 丁香欧美五月| 中文字幕最新亚洲高清| 国产三级黄色录像| 99国产精品一区二区三区| 男女视频在线观看网站免费 | 亚洲熟妇中文字幕五十中出| 亚洲av片天天在线观看| 欧美zozozo另类| 午夜两性在线视频| 久99久视频精品免费| 久久人人精品亚洲av| 亚洲精品国产一区二区精华液| 免费电影在线观看免费观看| 精品久久久久久久人妻蜜臀av| 久久中文字幕一级| 亚洲专区国产一区二区| 99精品久久久久人妻精品| 国产三级在线视频| 日本在线视频免费播放| 在线天堂中文资源库| 搡老熟女国产l中国老女人| 69av精品久久久久久| 精品不卡国产一区二区三区| 日本一区二区免费在线视频| 巨乳人妻的诱惑在线观看| 级片在线观看| 视频区欧美日本亚洲| 无限看片的www在线观看| 国产精品久久视频播放| 日韩欧美 国产精品| 久久久久久亚洲精品国产蜜桃av| 日韩欧美一区二区三区在线观看| 一级毛片精品| 最近最新中文字幕大全免费视频| 精品国产亚洲在线| av片东京热男人的天堂| 国产免费av片在线观看野外av| 免费观看精品视频网站| 性色av乱码一区二区三区2| 日本成人三级电影网站| 无人区码免费观看不卡| 国产精品香港三级国产av潘金莲| 亚洲欧美日韩高清在线视频| 制服丝袜大香蕉在线| 女同久久另类99精品国产91| 无人区码免费观看不卡| 日本 av在线| 看免费av毛片| 制服丝袜大香蕉在线| 国产精品亚洲一级av第二区| 无人区码免费观看不卡| 亚洲片人在线观看| 啦啦啦韩国在线观看视频| 国产亚洲精品久久久久5区| 一边摸一边做爽爽视频免费| x7x7x7水蜜桃| 黄色女人牲交| 在线国产一区二区在线| 精品久久久久久久末码| 欧美+亚洲+日韩+国产| 校园春色视频在线观看| 一二三四社区在线视频社区8| 国产伦在线观看视频一区| 亚洲精品久久国产高清桃花| 波多野结衣高清无吗| 午夜激情福利司机影院| 精品国产一区二区三区四区第35| 成人亚洲精品av一区二区| 美女午夜性视频免费| 国产精品 欧美亚洲| 亚洲熟妇熟女久久| 丁香六月欧美| 免费高清在线观看日韩| 超碰成人久久| 好男人在线观看高清免费视频 | 久久国产精品人妻蜜桃| 国产精品av久久久久免费| 激情在线观看视频在线高清| 亚洲人成电影免费在线| avwww免费| 丁香六月欧美| 日本 欧美在线| 日本免费一区二区三区高清不卡| www.精华液| 中文字幕人妻丝袜一区二区| 又黄又粗又硬又大视频| 禁无遮挡网站| 搡老岳熟女国产| 免费在线观看日本一区| 两个人视频免费观看高清| 黑人巨大精品欧美一区二区mp4| 国产日本99.免费观看| 宅男免费午夜| 搞女人的毛片| 麻豆成人午夜福利视频| 午夜福利欧美成人| 真人做人爱边吃奶动态| 一边摸一边抽搐一进一小说| 婷婷精品国产亚洲av在线| 成人18禁高潮啪啪吃奶动态图| 久久久精品欧美日韩精品| 变态另类成人亚洲欧美熟女| 97超级碰碰碰精品色视频在线观看| 国产精品,欧美在线| 久久久久免费精品人妻一区二区 | 国产精品一区二区精品视频观看| 99热6这里只有精品| 宅男免费午夜| 日韩国内少妇激情av| 亚洲色图 男人天堂 中文字幕| 搡老岳熟女国产| 中文字幕av电影在线播放| 最近最新免费中文字幕在线| 亚洲精品一区av在线观看| 亚洲在线自拍视频| 欧美在线一区亚洲| 国产精品综合久久久久久久免费| 国产精品永久免费网站| 99国产综合亚洲精品| 欧美日韩黄片免| 熟女少妇亚洲综合色aaa.| a级毛片a级免费在线| 亚洲 欧美 日韩 在线 免费| 91麻豆av在线| videosex国产| 久久精品91蜜桃| 亚洲av电影不卡..在线观看| 视频区欧美日本亚洲| 午夜两性在线视频| 亚洲欧美精品综合久久99| 国产精品日韩av在线免费观看| 黄色成人免费大全| 成人国语在线视频| 国产午夜精品久久久久久| 又紧又爽又黄一区二区| 国产单亲对白刺激| av有码第一页| 亚洲人成网站在线播放欧美日韩| 成年女人毛片免费观看观看9| 19禁男女啪啪无遮挡网站| 麻豆成人av在线观看| 可以免费在线观看a视频的电影网站| 亚洲精品粉嫩美女一区| 亚洲欧美精品综合一区二区三区| 十分钟在线观看高清视频www| 久久久水蜜桃国产精品网| 日本成人三级电影网站| 久久草成人影院| 女生性感内裤真人,穿戴方法视频| 少妇 在线观看| www国产在线视频色| 亚洲av日韩精品久久久久久密| 日韩中文字幕欧美一区二区| 午夜福利成人在线免费观看| 国产一区二区三区在线臀色熟女| 国产一级毛片七仙女欲春2 | 国产午夜精品久久久久久| 女同久久另类99精品国产91| 亚洲欧美精品综合久久99|