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

    三角模糊數優(yōu)化的貝葉斯水質模型評價方法研究*
    ——以太湖竺山灣緩沖帶濕地水質評價為例

    2016-03-13 01:38:44呂美婷李春華王亞空
    環(huán)境污染與防治 2016年11期
    關鍵詞:山灣緩沖帶水質評價

    呂美婷 葉 春 李春華 王亞空 干 磊

    (中國環(huán)境科學研究院湖泊工程技術中心,環(huán)境基準與風險評估國家重點實驗室,北京 100012)

    水質評價是水資源保護的重要基礎性工作。目前,有關水質評價方法的研究很多,常見的有:單因子污染指數法[1];綜合污染指數法,如綜合水質標識指數法、霍頓水質指數法[2]、布朗水質指數法[3]、內梅羅水污染指數法[4]等;模糊數學評價法,如模糊綜合指數法[5]、模糊層次分析法[6-7]、模糊綜合評價分級保證率法[8]等;灰色系統(tǒng)評價法[9],如灰色聚類法[10]、區(qū)域灰色決策法等;人工神經網絡評價法,如反向傳播算法(BP)網絡模型法[11-12]、Hopfield網絡模型法[13-14]、徑向基函數神經網絡(RBFNNs)模型法[15]等。由于水文條件、氣候因素等客觀原因,以及監(jiān)測資料不足、精度不夠、研究人員認識限制等主觀原因,給水質評價工作帶來了困難。目前,有關水質評價不確定性的研究主要運用模糊數學、概率統(tǒng)計、灰色系統(tǒng)和未確知數學等方法[16-18]。貝葉斯水質模型是解決水質評價過程中不確定性問題的新思路[19-21]。它是以概率統(tǒng)計為基礎,但又與傳統(tǒng)統(tǒng)計學理論不同。貝葉斯水質模型將所有參數作為隨機變量,先于觀測信息確定先驗分布[22-24]。模糊數學是利用隸屬度思想,將定性描述轉化為定量描述,更好地解決不確定性強、難以量化的問題。三角模糊數是適用于數據資料較少或不夠精確的情況的模糊數學法[25-26],被廣泛應用于水環(huán)境健康、河流水質污染、水環(huán)境風險等評價研究中[27]。本研究將三角模糊數引入到傳統(tǒng)的貝葉斯水質模型中,建立三角模糊數優(yōu)化的貝葉斯水質模型,并將此方法應用于太湖竺山灣緩沖帶濕地。具體以實際監(jiān)測數據均值及標準差為基礎建立三角模糊數α截集,代入優(yōu)化的貝葉斯水質模型;通過計算隸屬度,確定水質級別;并將水質評價結果與綜合水質標識指數法以及模糊標識指數法評價結果進行比較,為太湖竺山灣緩沖帶濕地水環(huán)境管理和生態(tài)建設決策提供理論基礎。

    1 評價方法

    1.1 三角模糊數原理及運算法則

    (1)

    式中:μA為模糊變量的隸屬度;x為某一變量。

    為方便計算,通常采用α截集對三角模糊數進行處理。α為可信度,且α∈[0,1],不同可信度水平對應不同數據區(qū)間,一般取α為0.9[29]。α截集處理的計算方法如下:

    (2)

    (3)

    (4)

    (5)

    (6)

    1.2 貝葉斯水質模型

    1.2.1 傳統(tǒng)貝葉斯水質模型

    傳統(tǒng)貝葉斯水質模型[30]見式(7):

    (7)

    式中:P(Bi,j|Aj)為第j水質指標的實測值為Aj時,水質屬于i類水的概率;Bi,j表示第j項水質指標i類水質的標準值;Aj表示第j項水質指標的實測值;j為水質指標;i為水質級別,包括Ⅰ~Ⅴ,共5類;P(Bi,j)為水質屬于i類水的先驗概率,先于觀測信息而推測得出;P(Aj|Bi,j)為第j項水質指標水質級別為i時,其值恰好為Aj的條件概率。Bi,j和Aj的單位視具體情況而定。

    1.2.2 三角模糊數優(yōu)化的貝葉斯水質模型

    (8)

    根據分層的思想及幾何距離的概念計算先驗概率和條件概率,計算公式如下:

    (9)

    (10)

    1.3 水質指標權重的確定

    采用變異系數法確定各水質指標的權重[31]906,[32],計算公式如下:

    (11)

    (12)

    1.4 水質級別及綜合得分的計算

    得出單項指標屬于各水質級別的條件概率及其在綜合水質評價中的權重,采取加權求和的方法計算水質屬于各級別的綜合后驗概率,計算公式如下:

    (13)

    得出水質屬于各級別的綜合后驗概率后,根據式(14)計算水質綜合得分:

    (14)

    計算出水質綜合得分后,可進一步計算水質綜合得分區(qū)間[P1,P2]對水質分級標準區(qū)間[Q1,Q2]的隸屬度,計算公式如下[28]1702,[33]:

    (15)

    式中:∣∣為區(qū)間的幾何長度;∩為兩個區(qū)間的交集。

    水質分級標準區(qū)間見表1。

    表1 水質分級標準區(qū)間

    2 實例分析

    2.1 研究區(qū)概況及采樣點

    太湖是中國第三大淡水湖泊,因地理位置及經濟社會地位突出,其生態(tài)環(huán)境問題備受關注[34-35]。竺山灣位于太湖西北部,是太湖上游來水的主要匯入湖區(qū)[36]。因竺山灣具有半封閉性,以及接受河道污水匯入等原因,已成為北部湖區(qū)污染最嚴重,且惡化速率最快的湖區(qū)之一[37]。竺山灣緩沖帶內生活污水、工業(yè)廢水以及農業(yè)非點源污染是造成水質污染的主要原因之一。

    2013年5月,在竺山灣緩沖帶核心示范區(qū)內選取了13個典型濕地斑塊進行采樣,對溫度、pH、溶解氧(DO)、TN、TP、氨氮、高錳酸鹽指數(IMn)7個指標進行監(jiān)測分析。TN采用堿性過硫酸鉀消解/紫外分光光度法測定;TP采用鉬銻鈧分光光度法測定;氨氮采用納氏試劑分光光度法測定;IMn按照《水質 高錳酸鹽指數的測定》(GB 11892—89)測定;溫度、pH和DO采用HORIBA多功能水質參數儀現場測定。采樣點分布見圖1,采樣點基本涵蓋了示范區(qū)內不同類型匯水區(qū)的濕地斑塊。

    2.2 太湖竺山灣緩沖帶濕地水質評價

    2.2.1 模型參數的選取及三角模糊化處理

    根據竺山灣緩沖帶濕地水質實際情況以及現有實測水質數據,選取TN、TP、氨氮、IMn和DO作為評價因子。根據三角模糊數對實測水質數據進行處理,為方便計算,采用α截集處理將三角模糊數轉化為可信度水平下對應的區(qū)間數,α取0.9,經過轉化的水質實測數據如表2所示。

    2.2.2 先驗概率及后驗概率的確定

    根據《地表水環(huán)境質量標準》(GB 3838—2002)中TN、TP、氨氮、IMn、DO 5項水質指標的各級水質標準以及式(9)計算先驗概率,結果如表3所示。

    根據式(10)計算出條件概率,并將結果代入式(8),即可計算出5個單項水質指標對各水質級別的后驗概率,結果見表4。從表4可以分析出各個濕地斑塊的重污染因子。例如:采樣點1#的TN對Ⅳ類水的后驗概率最大,而氨氮、TP、IMn和DO依次對Ⅲ、Ⅱ、Ⅲ、Ⅰ類水的后驗概率最大,因此TN為采樣點1#的重污染因子。采樣點4#、6#、7#、12#、13#同時存在TN和COD嚴重污染(COD污染由IMn表征)。采樣點1#、5#和9#僅存在TN嚴重污染,采樣點2#、3#和10#的COD污染嚴重,采樣點8#和11#同時存在TN、氨氮、TP、COD污染。

    根據水質分級標準矩陣和變異系數法,計算得出TN、TP、氨氮、IMn、DO 5項指標的權重分別為0.198 6、0.204 6、0.264 8、0.198 0、0.134 0。將單項水質指標對各水質級別的后驗概率及其對應的權重

    圖1 太湖竺山灣緩沖帶濕地采樣點分布Fig.1 Sampling site distribution in Zhushan Bay buffer zone wetland of Taihu Lake

    Table 2 Water quality monitoring data of Zhushan Bay buffer zone wetland treated withα-cut set technology

    mg/L

    表3 單項水質指標屬于各水質級別的先驗概率

    表4 單項水質指標屬于各水質級別的后驗概率

    表4 單項水質指標屬于各水質級別的后驗概率(續(xù))

    表5 各采樣點水質屬于各水質級別的綜合后驗概率

    代入式(13),即可得出各采樣點水質對各水質級別的綜合后驗概率,結果見表5。從表5可以看出,各采樣點的水質級別并不具有唯一性,只是相對來說,對某一級別的綜合后驗概率最大,如采樣點1#,水質級別為Ⅲ類的可能性最大。

    2.2.3 不同水質評價方法的評價結果比較

    采用余勛等[31]906對各級別水質的綜合打分方法,對Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ類水依次賦值5、4、3、2、1,然后根據式(14),計算各采樣點水質綜合得分,得出水質綜合得分區(qū)間,并根據式(15)計算水質綜合得分區(qū)間對水質分級標準區(qū)間的隸屬度,確定水質級別。同時,采用綜合水質標識指數法[38]和模糊標識指數法[39]對研究區(qū)內13個采樣點的水質情況進行了評價,并對不同水質評價方法產生的結果進行了比較(見表6)。

    從表6可以看出,3種評價方法對同一水質的評價結果存在差異,僅少數采樣點的水質評價結果一致,表明水質評價工作中,評價模型對評價結果有一定的不確定性??偟膩碚f,三角模糊數優(yōu)化的貝葉斯水質模型法對竺山灣緩沖帶濕地水質的評價結果與綜合水質標識指數法、模糊標識指數法得出的評價結果相似度均為46%左右;而綜合水質標識指數法和模糊標識指數法的評價結果相似度較低。

    綜合水質標識指數法和模糊標識指數法的共同優(yōu)點在于,百分位表示參與評價的水質指標中劣于目標水質的指標數,千分位表示綜合水質級別劣于目標水質的等級;而且,兩者都可以對劣Ⅴ類水進行評價。然而,也分別存在以下缺點:綜合水質標識指數法忽略了水質評價過程中的各種不確定性因素,如測量結果及分析過程的不確定性、模型結構的不確定性等等,而將各項評價指標進行平均,未考慮各指標對水質級別的貢獻存在一定差異,因而在水質評價定級上有一定誤差;模糊標識指數法在傳統(tǒng)模糊層次分析法的基礎上有所改進,改善了傳統(tǒng)方法低估水體污染程度方面的不足。三角模糊數優(yōu)化的貝葉斯水質模型在水質評價過程中綜合考慮了指標監(jiān)測數據的不確定性以及模型結構的不確定性。首先,根據幾何距離概念對先驗概率進行分層確定;其次,在確定單項指標后驗概率過程中,將實測值的平均值和標準差進行三角模糊化處理后引入評價模型,考慮了參評指標在一定置信度水平下的波動范圍,大大提高了水質評價結果的真實性。

    表6 基于不同水質評價方法的竺山灣緩沖帶濕地水質評價結果1)

    注:1)綜合水質標識指數法和模糊標識指數法的水質綜合得分越高,則水質越差;相反,三角模糊數優(yōu)化的貝葉斯水質模型法的水質綜合得分越高,則水質越好。

    三角模糊數優(yōu)化的貝葉斯水質模型對太湖竺山灣濕地水質的評價結果顯示:13個采樣點水質級別評價結果為Ⅲ~Ⅳ類水,而根據《太湖流域水環(huán)境綜合治理總體方案》和《江蘇省太湖流域水環(huán)境綜合治理實施方案》對水功能區(qū)級別的劃分,竺山灣隸屬一級功能區(qū),目標水質級別為Ⅰ類。因此,竺山灣緩沖帶濕地內的13個濕地斑塊的水質均未達到目標水質標準。

    3 討論與結論

    3.1 討 論

    三角模糊數考慮了實測值的隨機誤差、采樣點具體位置引起的監(jiān)測結果差異等不確定性因素,尤其適合處理資料較少、精度不夠等情況下的水質信息。因此,將三角模糊數引入太湖竺山灣緩沖帶濕地的水質評價過程,使評價結果更為客觀合理。

    貝葉斯水質模型利用已知信息,統(tǒng)計推斷出多種可能性,并以最大可能性做出最終判斷,是一種針對小樣本的、可行的水質評價方法。在水質評價中則體現為根據水質綜合得分區(qū)間對水質分級標準區(qū)間的隸屬度,判定水質級別。傳統(tǒng)貝葉斯水質模型中,先驗概率均勻分布。本研究根據GB 3838—2002各指標水質分級標準,引入幾何距離的概念,分層計算先驗概率;同時以各因子加權求和的方法計算得出綜合后驗概率,全面考慮了各參評水質指標對水質評價定級的貢獻率,使得評價結果更加真實可靠。

    三角模糊數優(yōu)化的貝葉斯水質模型法綜合了三角模糊數和傳統(tǒng)貝葉斯水質模型的優(yōu)點,同時考慮了水質評價模型結構和參數的不確定性。得出的水質綜合得分及水質級別均是基于模糊數學思想的,帶有隸屬度信息的,能更好地反映竺山灣緩沖帶濕地水質真實情況。

    3.2 結 論

    竺山灣緩沖帶濕地內水體水質都為Ⅲ~Ⅳ類,未能符合Ⅰ類目標水質要求。從總體上看,竺山灣緩沖帶濕地水體的TN和COD污染較嚴重,TP、氨氮和DO并非重污染因子。因此,在后期的生態(tài)修復及水質改善研究工作中應主要針對削減氮和有機物污染開展。

    [1] 徐祖信.我國河流單因子水質標識指數評價方法研究[J].同濟大學學報(自然科學版),2005,33(3):321-325.

    [2] RAMOSJ A. Applications of TLS and related methods in the environmental sciences[J].Computational Statistics & Data Analysis,2007,52(2):1234-1267.

    [3] 周敏.飲用水水質指數分析及其在基層的應用研究[D].上海:復旦大學,2011.

    [4] 黃海東,張克峰.小城鎮(zhèn)水源水質評價方法選擇方案探討[J].水利科技與經濟,2010,16(7):736-738,740.

    [5] 楊靜.改進的模糊綜合評價法在水質評價中的應用[D].重慶:重慶大學,2014.

    [6] 徐兵兵,張妙仙,王肖肖.改進的模糊層次分析法在南苕溪臨安段水質評價中的應用[J].環(huán)境科學學報,2011,31(9):2066-2072.

    [7] 盧文喜,李迪,張蕾,等.基于層次分析法的模糊綜合評價在水質評價中的應用[J].節(jié)水灌溉,2011(3):43-46.

    [8] 張美華,陳宏.幾種模糊聚類法在環(huán)境質量綜合評價中的應用[J].重慶環(huán)境科學,1999,21(3):13-16.

    [9] 陸洲,夏秋穎,周琳,等.等斜率灰色聚類法在地面水環(huán)境質量評價中的應用[J].環(huán)境保護科學,2000,26(5):43-46.

    [10] 賀北方,王效宇,賀曉菊,等.基于灰色聚類決策的水質評價方法[J].鄭州大學學報(工學版),2002,23(1):10-13.

    [11] 李文娟.改進BP神經網絡在水質評價中的應用研究[D].重慶:重慶理工大學,2011.

    [12] 黃勝偉,董曼玲.自適應變步長BP神經網絡在水質評價中的應用[J].水利學報,2002,33(10):119-123.

    [13] 初海波,盧文喜,尹津航,等.BP網絡、Hopfield網絡在水質評價應用中的比較研究[J].中國農村水利水電,2011(10):70-72.

    [14] 盧文喜,初海波,王喜華,等.基于因子分析的Hopfield神經網絡在水質評價的應用[J].水土保持通報,2012,32(1):197-200.

    [15] 周豐,郭懷成,劉永,等.基于多元統(tǒng)計分析和RBFNNs的水質評價方法[J].環(huán)境科學學報,2007,27(5):846-853.

    [16] GRONEWOLD A D,STOW C A,VIJAYAVEL K A,et al.Differentiatingenterococcusconcentration spatial,temporal,and analytical variability in recreational waters[J].Water Research,2013,47(7):2141-2152.

    [17] GRONEWOLD A D,BORSUK M E.Improving water quality assessments through a hierarchical Bayesian analysis of variability[J].Environmental Science & Technology,2010,44(20):7858-7864.

    [18] 尹海龍,徐祖信.河流綜合水質評價方法比較研究[J].長江流域資源與環(huán)境,2008,17(5):729-733.

    [19] WELLEN C,ARHONDITSIS G B,LONG Tanya,et al.Quantifying the uncertainty of nonpoint source attribution in distributed water quality models:a Bayesian assessment of SWAT’s sediment export predictions[J].Journal of Hydrology,2014,519:3353-3368.

    [20] HADDAD K,EGODAWATTA P,RAHMAN A A.Uncertainty analysis of pollutant build-up modelling based on a Bayesian weighted least squares approach[J].Science of the Total Environment,2013,449:410-417.

    [21] 孫鵬程,陳吉寧.基于貝葉斯網絡的河流突發(fā)性水質污染事故風險評估[J].環(huán)境科學,2009,30(1):47-51.

    [22] 黃凱,張曉玲.貝葉斯方法在水環(huán)境系統(tǒng)不確定性分析中的應用述評[J].水電能源科學,2012,30(9):47-49.

    [23] 廖杰,汪嘉楊,丁晶.基于改進貝葉斯模型的四川省主要河流水質評價[J].四川師范大學學報(自然科學版),2009,32(4):518-521.

    [24] 郁章文.基于貝葉斯方法的三峽庫區(qū)水環(huán)境質量研究[D].濟南:山東大學,2008.

    [25] SANGAIAHA A K,THANGAVELU A K.An exploration of FMCDM approach for evaluating the outcome/success of GSD projects[J].Central European Journal of Engineering,2013,3(3):419-435.

    [26] HUANG Guangqiu,PAN Changbo.Introduction of probabilistic risk assessment approach to analyze the safety of MVS[J].Procedia Engineering,2011,26:2027-2031.

    [27] 趙鋒霞.基于不確定性參數的水環(huán)境健康風險評價方法及應用[D].大連:遼寧師范大學,2013.

    [28] 樊夢佳,袁興中,?;勰?,等.基于三角模糊數的河流沉積物中重金屬污染評價模型[J].環(huán)境科學學報,2010,30(8).

    [29] 李如忠,洪天求,金菊良.河流水質模糊風險評價模型研究[J].武漢理工大學學報,2007,29(2):43-46.

    [30] YANG Likun,ZHAO Xinhua,PENG Sen,et al.Integration of Bayesian analysis for eutrophication prediction and assessment in a landscape lake[J].Environmental Monitoring and Assessment,2014,187(1):1-19.

    [31] 余勛,梁婕,曾光明,等.基于三角模糊數的貝葉斯水質評價模型[J].環(huán)境科學學報,2013,33(3).

    [32] 屈衛(wèi)軍.基于變異系數法的灰色關聯決策模型在節(jié)水工程方案優(yōu)選中的應用[J].地下水,2010,32(5):50-51.

    [33] ?;勰?,袁興中,曾光明,等.基于區(qū)間數的河流水環(huán)境健康風險模糊綜合評價模型[J].環(huán)境科學學報,2009,29(7):1527-1533.

    [34] 許嫚,劉曉收,劉清河,等.太湖竺山灣湖泊緩沖帶內濕地浮游植物季節(jié)變化[J].海洋湖沼通報,2014(4):41-51.

    [35] 呂美婷,趙德華,楊棠武,等.湖泊重要入湖河流穩(wěn)定水像元提取技術[J].環(huán)境監(jiān)控與預警,2013,5(5):6-9.

    [36] 張明禮,楊浩,林振山,等.太湖竺山灣底泥中有害物質含量與環(huán)境污染評價[J].中國環(huán)境科學,2011,31(5):852-857.

    [37] 王雯雯,姜霞,王書航,等.太湖竺山灣污染底泥環(huán)保疏浚深度的推算[J].中國環(huán)境科學,2011,31(6):1013-1018.

    [38] 徐祖信.我國河流綜合水質標識指數評價方法研究[J].同濟大學學報(自然科學版),2005,33(4):482-488.

    [39] 王肖肖,張妙仙,徐兵兵.模糊標識指數與對應分析法在水質評價中的聯合應用[J].環(huán)境科學學報,2012,32(5):1227-1235.

    猜你喜歡
    山灣緩沖帶水質評價
    統(tǒng)編語文教材九(下)第五單元拓展閱讀
    河岸緩沖帶植被布局對氮流失的影響
    統(tǒng)編語文教材九(上)第五單元拓展閱讀
    滻灞河水質評價方法研究
    太湖流域流經不同類型緩沖帶入湖河流秋、冬季氮污染特征?
    湖泊科學(2016年6期)2016-11-24 09:31:28
    基于概率統(tǒng)計和模糊綜合評價法的水質評價模型及其應用——以拉薩河水質評價為例
    基于SPAM的河流水質評價模型
    臨澤縣紅山灣水庫工程壩基處理工程技術的應用
    坡耕地作物緩沖帶對徑流和氮流失的影響分析
    濱岸緩沖帶去除農業(yè)面源污染技術研究
    午夜免费鲁丝| 亚洲av不卡在线观看| 天堂网av新在线| 白带黄色成豆腐渣| 免费观看a级毛片全部| 嫩草影院精品99| 国内精品宾馆在线| 在线观看一区二区三区激情| 久久鲁丝午夜福利片| 成人毛片a级毛片在线播放| av免费在线看不卡| 91久久精品国产一区二区成人| 韩国高清视频一区二区三区| 亚洲av中文字字幕乱码综合| 婷婷色麻豆天堂久久| 日韩av免费高清视频| 午夜福利视频1000在线观看| 午夜精品一区二区三区免费看| 欧美日韩综合久久久久久| 在线观看一区二区三区| 午夜亚洲福利在线播放| 国产亚洲一区二区精品| 日日摸夜夜添夜夜添av毛片| 边亲边吃奶的免费视频| 亚洲内射少妇av| 日韩亚洲欧美综合| 亚洲三级黄色毛片| 边亲边吃奶的免费视频| 高清av免费在线| 国产亚洲精品久久久com| 涩涩av久久男人的天堂| 天堂中文最新版在线下载 | 日韩成人av中文字幕在线观看| 51国产日韩欧美| 免费黄频网站在线观看国产| 插逼视频在线观看| 少妇人妻一区二区三区视频| 日韩欧美精品v在线| 97人妻精品一区二区三区麻豆| 街头女战士在线观看网站| 草草在线视频免费看| 欧美少妇被猛烈插入视频| 亚洲精品国产成人久久av| 超碰av人人做人人爽久久| 国模一区二区三区四区视频| 丝袜脚勾引网站| 美女视频免费永久观看网站| 能在线免费看毛片的网站| 美女xxoo啪啪120秒动态图| videossex国产| 亚洲国产精品成人久久小说| 一级二级三级毛片免费看| 精品国产乱码久久久久久小说| www.av在线官网国产| 国产色爽女视频免费观看| 啦啦啦中文免费视频观看日本| 91在线精品国自产拍蜜月| 亚洲高清免费不卡视频| 亚洲aⅴ乱码一区二区在线播放| 国产视频首页在线观看| 国产亚洲91精品色在线| 欧美性猛交╳xxx乱大交人| 久久人人爽av亚洲精品天堂 | 高清日韩中文字幕在线| 免费观看a级毛片全部| 婷婷色综合大香蕉| 亚洲国产av新网站| 99久久中文字幕三级久久日本| 高清日韩中文字幕在线| 亚洲av免费在线观看| 激情五月婷婷亚洲| 国产视频首页在线观看| 亚洲精品影视一区二区三区av| 国产精品久久久久久精品电影小说 | 成人亚洲精品av一区二区| 色哟哟·www| 人人妻人人澡人人爽人人夜夜| 欧美变态另类bdsm刘玥| 亚洲精品一区蜜桃| 青春草视频在线免费观看| 亚洲av二区三区四区| 亚洲精品久久久久久婷婷小说| 高清日韩中文字幕在线| 激情 狠狠 欧美| 精品久久久精品久久久| 美女高潮的动态| 特大巨黑吊av在线直播| 80岁老熟妇乱子伦牲交| 伦精品一区二区三区| 国产爽快片一区二区三区| 亚洲国产色片| 一区二区三区四区激情视频| 美女xxoo啪啪120秒动态图| 少妇人妻久久综合中文| 久久久成人免费电影| 亚洲欧美日韩无卡精品| 在线观看三级黄色| 国产成人免费无遮挡视频| 天天躁日日操中文字幕| 亚洲最大成人手机在线| 午夜精品国产一区二区电影 | 男人添女人高潮全过程视频| 国产在线男女| 亚洲成色77777| 免费看不卡的av| 天堂网av新在线| 最近最新中文字幕大全电影3| 久久久久久久亚洲中文字幕| 蜜桃亚洲精品一区二区三区| 水蜜桃什么品种好| 日韩欧美精品免费久久| 精品久久久久久久末码| 欧美zozozo另类| 亚洲精品久久午夜乱码| 观看美女的网站| 99九九线精品视频在线观看视频| 亚洲欧美日韩卡通动漫| 天美传媒精品一区二区| 成人国产av品久久久| 国内少妇人妻偷人精品xxx网站| 免费观看在线日韩| 插逼视频在线观看| 成人免费观看视频高清| 欧美最新免费一区二区三区| 中文天堂在线官网| 中文欧美无线码| 精品一区二区免费观看| 国产精品人妻久久久影院| 国产亚洲最大av| 国产亚洲av嫩草精品影院| 国产精品久久久久久久久免| 亚洲国产成人一精品久久久| 美女被艹到高潮喷水动态| 肉色欧美久久久久久久蜜桃 | 建设人人有责人人尽责人人享有的 | av又黄又爽大尺度在线免费看| 国产一区二区亚洲精品在线观看| 黄色怎么调成土黄色| 日本欧美国产在线视频| 亚洲av免费高清在线观看| 亚洲精品成人av观看孕妇| 欧美成人午夜免费资源| kizo精华| av在线蜜桃| 亚洲激情五月婷婷啪啪| 亚洲av成人精品一二三区| 97在线视频观看| 最后的刺客免费高清国语| 国产v大片淫在线免费观看| 波多野结衣巨乳人妻| 亚洲av日韩在线播放| 国产 一区精品| 久久精品综合一区二区三区| 高清毛片免费看| 高清午夜精品一区二区三区| 久久综合国产亚洲精品| 免费黄频网站在线观看国产| 亚洲熟女精品中文字幕| 久久久久网色| h日本视频在线播放| 黑人高潮一二区| 久久精品国产亚洲网站| 免费观看av网站的网址| 啦啦啦中文免费视频观看日本| 久久精品国产亚洲av天美| 国产白丝娇喘喷水9色精品| 看非洲黑人一级黄片| 69av精品久久久久久| av线在线观看网站| av在线观看视频网站免费| 一二三四中文在线观看免费高清| 99久久精品一区二区三区| 亚洲成人精品中文字幕电影| 亚洲欧美中文字幕日韩二区| 特级一级黄色大片| 啦啦啦在线观看免费高清www| 日韩一区二区三区影片| 国产免费福利视频在线观看| 欧美激情在线99| 国产高清三级在线| 国产欧美亚洲国产| 搡女人真爽免费视频火全软件| 亚洲四区av| 18+在线观看网站| 精华霜和精华液先用哪个| 国模一区二区三区四区视频| 亚洲在久久综合| 各种免费的搞黄视频| 麻豆乱淫一区二区| 亚洲成人av在线免费| 久久久久久国产a免费观看| 国产一级毛片在线| 日韩欧美 国产精品| 国产精品精品国产色婷婷| 秋霞在线观看毛片| 2021少妇久久久久久久久久久| 亚洲欧美日韩无卡精品| 777米奇影视久久| 日日摸夜夜添夜夜爱| 久久久色成人| 国产欧美日韩一区二区三区在线 | 成人免费观看视频高清| 精品国产乱码久久久久久小说| av在线播放精品| 成人二区视频| 欧美激情国产日韩精品一区| 婷婷色麻豆天堂久久| av一本久久久久| 午夜激情福利司机影院| 美女被艹到高潮喷水动态| 成人毛片a级毛片在线播放| 日韩av免费高清视频| 欧美日韩亚洲高清精品| 我的老师免费观看完整版| 亚洲色图综合在线观看| 欧美97在线视频| 久久国产乱子免费精品| av女优亚洲男人天堂| 特级一级黄色大片| 国产亚洲最大av| videos熟女内射| 国产黄频视频在线观看| 免费观看的影片在线观看| 国产毛片在线视频| 亚洲国产欧美在线一区| 亚洲精品一区蜜桃| 91午夜精品亚洲一区二区三区| 亚洲最大成人av| 欧美3d第一页| 不卡视频在线观看欧美| 全区人妻精品视频| 欧美日韩视频高清一区二区三区二| 免费黄频网站在线观看国产| 国产精品av视频在线免费观看| 99热6这里只有精品| 99视频精品全部免费 在线| 免费不卡的大黄色大毛片视频在线观看| 亚洲最大成人手机在线| 久久久久久久久久人人人人人人| 亚洲精品乱久久久久久| 如何舔出高潮| 亚洲经典国产精华液单| 久久久久久伊人网av| 日本-黄色视频高清免费观看| eeuss影院久久| 国产av不卡久久| 一级毛片我不卡| 午夜精品一区二区三区免费看| 男女啪啪激烈高潮av片| 欧美日韩视频高清一区二区三区二| 99久久精品国产国产毛片| 秋霞伦理黄片| 国产一区二区三区av在线| 王馨瑶露胸无遮挡在线观看| 欧美一级a爱片免费观看看| 91久久精品电影网| 舔av片在线| 久久99热6这里只有精品| 久久久久国产网址| 神马国产精品三级电影在线观看| 69人妻影院| 免费看av在线观看网站| 97超碰精品成人国产| 亚洲精品久久久久久婷婷小说| 久久久精品欧美日韩精品| 韩国高清视频一区二区三区| 日韩国内少妇激情av| 久久久久精品久久久久真实原创| 嫩草影院入口| 日本一本二区三区精品| 插逼视频在线观看| 午夜激情久久久久久久| 国产精品久久久久久精品电影小说 | av在线播放精品| 丝瓜视频免费看黄片| 国产精品一区二区三区四区免费观看| 99热国产这里只有精品6| 成人高潮视频无遮挡免费网站| 免费电影在线观看免费观看| 免费黄网站久久成人精品| 久久久午夜欧美精品| 一区二区三区乱码不卡18| 国产精品99久久久久久久久| 偷拍熟女少妇极品色| 夜夜爽夜夜爽视频| 免费电影在线观看免费观看| 高清毛片免费看| 日韩欧美精品v在线| 亚洲人与动物交配视频| av.在线天堂| 欧美日韩视频高清一区二区三区二| 精品视频人人做人人爽| 在线观看三级黄色| 天堂网av新在线| 99热网站在线观看| 亚洲人成网站在线播| 91aial.com中文字幕在线观看| 简卡轻食公司| 中文天堂在线官网| 在线免费十八禁| 欧美日韩综合久久久久久| 日日摸夜夜添夜夜爱| 亚洲,欧美,日韩| 精品一区在线观看国产| 成人午夜精彩视频在线观看| 精品午夜福利在线看| 别揉我奶头 嗯啊视频| 国产亚洲5aaaaa淫片| 国产精品嫩草影院av在线观看| 中国美白少妇内射xxxbb| 久久久成人免费电影| 亚洲最大成人手机在线| 看黄色毛片网站| 精品熟女少妇av免费看| 大片电影免费在线观看免费| 麻豆成人午夜福利视频| tube8黄色片| 交换朋友夫妻互换小说| 少妇裸体淫交视频免费看高清| 在线观看人妻少妇| 免费看a级黄色片| 欧美成人a在线观看| 国产成人精品一,二区| 蜜臀久久99精品久久宅男| 视频中文字幕在线观看| 中文在线观看免费www的网站| 免费看a级黄色片| 免费观看的影片在线观看| 91精品一卡2卡3卡4卡| 成人黄色视频免费在线看| 国产高清有码在线观看视频| 老女人水多毛片| 午夜日本视频在线| 亚洲婷婷狠狠爱综合网| 久热这里只有精品99| av免费观看日本| 五月天丁香电影| 国产欧美亚洲国产| 精品国产乱码久久久久久小说| 一本久久精品| 欧美变态另类bdsm刘玥| 欧美bdsm另类| 大香蕉97超碰在线| 91狼人影院| 七月丁香在线播放| 亚洲精品日本国产第一区| 美女被艹到高潮喷水动态| 亚洲真实伦在线观看| 午夜免费鲁丝| 亚洲人成网站高清观看| 亚洲精品视频女| 好男人视频免费观看在线| 午夜精品一区二区三区免费看| 亚洲自拍偷在线| 国语对白做爰xxxⅹ性视频网站| 男人添女人高潮全过程视频| 亚洲一区二区三区欧美精品 | 国产高清三级在线| 婷婷色麻豆天堂久久| 校园人妻丝袜中文字幕| 成年av动漫网址| 在线观看av片永久免费下载| 亚洲最大成人手机在线| 大香蕉久久网| 天天躁日日操中文字幕| 黄色一级大片看看| 亚洲aⅴ乱码一区二区在线播放| 岛国毛片在线播放| 国产精品久久久久久av不卡| 免费av不卡在线播放| 久久6这里有精品| 在线精品无人区一区二区三 | 久热久热在线精品观看| 国产一区二区在线观看日韩| 久久精品国产a三级三级三级| 搡女人真爽免费视频火全软件| 亚洲不卡免费看| 男女下面进入的视频免费午夜| 在线亚洲精品国产二区图片欧美 | 大又大粗又爽又黄少妇毛片口| xxx大片免费视频| 成人综合一区亚洲| 久久精品综合一区二区三区| 91狼人影院| 视频中文字幕在线观看| 伦理电影大哥的女人| 日日摸夜夜添夜夜爱| 久久久久久久久久久免费av| 一级毛片我不卡| 免费观看av网站的网址| 国产精品一区www在线观看| 久久久亚洲精品成人影院| 国产精品99久久99久久久不卡 | 亚洲国产成人一精品久久久| 99热国产这里只有精品6| 欧美成人午夜免费资源| 一本色道久久久久久精品综合| tube8黄色片| av又黄又爽大尺度在线免费看| 精品久久久噜噜| 在线观看三级黄色| 午夜免费观看性视频| 久久ye,这里只有精品| 免费av观看视频| 一级毛片久久久久久久久女| 精品酒店卫生间| 一区二区三区免费毛片| 亚洲欧美一区二区三区黑人 | 国产在视频线精品| 亚洲精品aⅴ在线观看| 欧美xxxx性猛交bbbb| 久久99精品国语久久久| 91精品国产九色| 精品国产三级普通话版| 亚洲av成人精品一区久久| 欧美精品国产亚洲| 美女被艹到高潮喷水动态| 噜噜噜噜噜久久久久久91| 伦精品一区二区三区| 欧美激情在线99| 欧美日韩综合久久久久久| 国内精品宾馆在线| 婷婷色麻豆天堂久久| 亚洲精品中文字幕在线视频 | 久久久久久久久久人人人人人人| 日韩亚洲欧美综合| 日产精品乱码卡一卡2卡三| 热re99久久精品国产66热6| 五月伊人婷婷丁香| 中文在线观看免费www的网站| 啦啦啦在线观看免费高清www| 波多野结衣巨乳人妻| 一级毛片久久久久久久久女| 亚洲经典国产精华液单| 久久久色成人| 一区二区三区精品91| 国产精品伦人一区二区| 99久久精品热视频| 男女无遮挡免费网站观看| 亚洲精品国产色婷婷电影| 亚洲在久久综合| 亚洲婷婷狠狠爱综合网| 纵有疾风起免费观看全集完整版| 久久精品久久久久久久性| 26uuu在线亚洲综合色| 又爽又黄a免费视频| 久久人人爽人人爽人人片va| 国产亚洲一区二区精品| 亚洲怡红院男人天堂| 在线观看一区二区三区激情| 搡老乐熟女国产| 九草在线视频观看| 成人无遮挡网站| 全区人妻精品视频| 国产成人aa在线观看| 亚洲欧美日韩无卡精品| 别揉我奶头 嗯啊视频| 精品午夜福利在线看| 人人妻人人看人人澡| 69人妻影院| 99视频精品全部免费 在线| 国产淫片久久久久久久久| 久久99精品国语久久久| 久久精品国产a三级三级三级| 午夜精品一区二区三区免费看| 最近最新中文字幕免费大全7| 别揉我奶头 嗯啊视频| 亚洲av福利一区| av专区在线播放| av在线天堂中文字幕| 日韩欧美一区视频在线观看 | 亚洲国产欧美在线一区| 尾随美女入室| 性插视频无遮挡在线免费观看| 国产老妇女一区| 国产老妇伦熟女老妇高清| 亚洲精品久久久久久婷婷小说| 日本黄色片子视频| 日本wwww免费看| 男女无遮挡免费网站观看| 一级毛片久久久久久久久女| 插逼视频在线观看| 亚洲精品成人av观看孕妇| 少妇丰满av| 国产免费又黄又爽又色| 少妇裸体淫交视频免费看高清| 亚洲成人av在线免费| 乱系列少妇在线播放| 中文字幕av成人在线电影| 美女视频免费永久观看网站| 亚州av有码| 国产精品人妻久久久久久| 国产成人免费无遮挡视频| 日本黄大片高清| 美女高潮的动态| 亚洲av男天堂| 国产在线男女| 欧美zozozo另类| 三级国产精品欧美在线观看| 亚洲精品中文字幕在线视频 | 天堂中文最新版在线下载 | av网站免费在线观看视频| 啦啦啦啦在线视频资源| 高清日韩中文字幕在线| 亚洲欧美成人精品一区二区| 欧美成人一区二区免费高清观看| 亚洲精品自拍成人| av线在线观看网站| 久久久久久国产a免费观看| 免费大片黄手机在线观看| 禁无遮挡网站| 国产亚洲5aaaaa淫片| 伊人久久国产一区二区| 久久精品综合一区二区三区| 亚洲色图av天堂| 男人爽女人下面视频在线观看| 亚洲熟女精品中文字幕| 午夜精品一区二区三区免费看| 在线观看美女被高潮喷水网站| 国产av码专区亚洲av| 最近的中文字幕免费完整| 日韩不卡一区二区三区视频在线| 国产伦精品一区二区三区四那| 欧美zozozo另类| 边亲边吃奶的免费视频| 久久久久久久国产电影| 午夜激情久久久久久久| 免费大片黄手机在线观看| 久久久久网色| 日韩伦理黄色片| 成人欧美大片| 国产精品久久久久久av不卡| 亚洲精品乱码久久久久久按摩| 日本三级黄在线观看| 男女国产视频网站| 色网站视频免费| 免费在线观看成人毛片| 我的女老师完整版在线观看| 97超碰精品成人国产| 一级毛片我不卡| 国产黄频视频在线观看| 99热这里只有是精品50| 久久精品久久精品一区二区三区| 交换朋友夫妻互换小说| 欧美高清性xxxxhd video| 丰满乱子伦码专区| 九草在线视频观看| 国产成人aa在线观看| 国产精品av视频在线免费观看| 另类亚洲欧美激情| 插阴视频在线观看视频| 国产精品成人在线| 午夜精品一区二区三区免费看| 日本爱情动作片www.在线观看| 国产免费福利视频在线观看| 国产v大片淫在线免费观看| 国产av不卡久久| 99久久九九国产精品国产免费| 国产美女午夜福利| 亚洲国产av新网站| 尾随美女入室| 热99国产精品久久久久久7| av免费观看日本| 国产淫语在线视频| 91久久精品国产一区二区成人| 久久精品人妻少妇| 九九爱精品视频在线观看| 亚洲精品国产色婷婷电影| 精品国产露脸久久av麻豆| 久久99热6这里只有精品| 国产成年人精品一区二区| 国产成人午夜福利电影在线观看| 天美传媒精品一区二区| 黄色日韩在线| 亚洲欧美中文字幕日韩二区| 国产精品无大码| 熟妇人妻不卡中文字幕| 国产精品国产三级国产av玫瑰| 成人美女网站在线观看视频| 激情五月婷婷亚洲| 伦精品一区二区三区| 在线观看免费高清a一片| 色网站视频免费| 成年人午夜在线观看视频| 汤姆久久久久久久影院中文字幕| 成人鲁丝片一二三区免费| 看非洲黑人一级黄片| 别揉我奶头 嗯啊视频| 亚洲国产欧美在线一区| 免费av观看视频| 国产av国产精品国产| 国产老妇女一区| 国产老妇伦熟女老妇高清| 国产黄色免费在线视频| 国产精品蜜桃在线观看| 亚州av有码| 97超碰精品成人国产| 91久久精品国产一区二区成人| 天堂中文最新版在线下载 | 国产精品99久久久久久久久| 美女内射精品一级片tv| 永久免费av网站大全| 成年免费大片在线观看| 免费播放大片免费观看视频在线观看| 极品少妇高潮喷水抽搐| 欧美少妇被猛烈插入视频| 亚洲在久久综合| 亚洲欧美清纯卡通| freevideosex欧美| 久久精品国产亚洲av天美| 亚洲国产欧美在线一区| 成年版毛片免费区| 青春草国产在线视频| 又粗又硬又长又爽又黄的视频| 老司机影院毛片|