• <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è)面源污染技術研究
    美女cb高潮喷水在线观看| 国产精品久久久久久久电影| 男女无遮挡免费网站观看| 在线 av 中文字幕| 好男人视频免费观看在线| 99热国产这里只有精品6| 91久久精品电影网| 一级爰片在线观看| 色5月婷婷丁香| 精品国产露脸久久av麻豆| 欧美三级亚洲精品| 精品一区二区免费观看| 成人综合一区亚洲| 成人综合一区亚洲| 成人亚洲欧美一区二区av| 免费观看av网站的网址| 成人国产麻豆网| 搡女人真爽免费视频火全软件| 亚洲自偷自拍三级| 精品国产国语对白av| 久久久久网色| 综合色丁香网| 国产一区二区三区av在线| 亚洲无线观看免费| 午夜免费男女啪啪视频观看| 男女无遮挡免费网站观看| 久久久久久人妻| 亚洲精品乱码久久久v下载方式| 在线 av 中文字幕| 国产精品久久久久成人av| 一区二区三区四区激情视频| 日本av手机在线免费观看| 亚洲精品国产成人久久av| 十分钟在线观看高清视频www | 97超碰精品成人国产| 亚洲精品一区蜜桃| 晚上一个人看的免费电影| 免费黄色在线免费观看| 美女cb高潮喷水在线观看| 成人无遮挡网站| 99九九在线精品视频 | 内射极品少妇av片p| 国产精品不卡视频一区二区| 一级毛片aaaaaa免费看小| 国产成人freesex在线| 99国产精品免费福利视频| 国产精品.久久久| 在线天堂最新版资源| 一级,二级,三级黄色视频| 亚洲av成人精品一二三区| 777米奇影视久久| 高清不卡的av网站| 91久久精品国产一区二区三区| 麻豆精品久久久久久蜜桃| 日韩伦理黄色片| 亚洲真实伦在线观看| 麻豆成人av视频| 国产高清三级在线| 99久国产av精品国产电影| 欧美另类一区| 亚洲熟女精品中文字幕| 成人毛片a级毛片在线播放| 国产av码专区亚洲av| av有码第一页| 国产av码专区亚洲av| 交换朋友夫妻互换小说| 美女大奶头黄色视频| 日韩av在线免费看完整版不卡| 天天躁夜夜躁狠狠久久av| 一级毛片黄色毛片免费观看视频| 午夜视频国产福利| 一级片'在线观看视频| 国产视频内射| 亚洲欧美成人综合另类久久久| 成人午夜精彩视频在线观看| 热99国产精品久久久久久7| 欧美日韩在线观看h| 97超碰精品成人国产| 亚洲美女搞黄在线观看| 高清不卡的av网站| 欧美精品国产亚洲| 你懂的网址亚洲精品在线观看| 久久久国产欧美日韩av| 大陆偷拍与自拍| 青春草亚洲视频在线观看| 69精品国产乱码久久久| 日日啪夜夜撸| 久久久久久久久久久久大奶| 22中文网久久字幕| 少妇丰满av| av又黄又爽大尺度在线免费看| √禁漫天堂资源中文www| 国产免费视频播放在线视频| 丝袜在线中文字幕| 国产精品99久久久久久久久| 观看av在线不卡| 亚洲国产精品一区三区| 亚洲精品第二区| 国产乱人偷精品视频| 国产在线男女| 国产精品一区二区三区四区免费观看| 日韩精品免费视频一区二区三区 | 亚洲第一av免费看| 亚洲欧美成人综合另类久久久| 久久精品国产亚洲av涩爱| 日韩精品免费视频一区二区三区 | 久久精品国产自在天天线| 久久久久视频综合| 99九九在线精品视频 | 亚洲av国产av综合av卡| 这个男人来自地球电影免费观看 | 两个人免费观看高清视频 | a 毛片基地| 亚洲情色 制服丝袜| 久久久久久久精品精品| 99久国产av精品国产电影| 国产精品一二三区在线看| 建设人人有责人人尽责人人享有的| 新久久久久国产一级毛片| 国产av精品麻豆| 久久99蜜桃精品久久| 国产视频首页在线观看| 晚上一个人看的免费电影| 国产在线免费精品| 少妇被粗大猛烈的视频| 日日爽夜夜爽网站| 精品亚洲乱码少妇综合久久| 纵有疾风起免费观看全集完整版| 美女xxoo啪啪120秒动态图| 日本黄色片子视频| 大陆偷拍与自拍| 日韩中文字幕视频在线看片| 国产美女午夜福利| 黄色配什么色好看| 欧美日韩精品成人综合77777| 这个男人来自地球电影免费观看 | 麻豆成人av视频| 亚洲成色77777| 汤姆久久久久久久影院中文字幕| 日韩伦理黄色片| 国产精品国产三级专区第一集| 秋霞伦理黄片| 午夜免费鲁丝| 日本与韩国留学比较| av天堂中文字幕网| 女人久久www免费人成看片| 精品少妇黑人巨大在线播放| 免费看光身美女| a级一级毛片免费在线观看| 久久鲁丝午夜福利片| 永久网站在线| 久久精品久久久久久噜噜老黄| 国产精品欧美亚洲77777| 丰满迷人的少妇在线观看| 老司机亚洲免费影院| 黄色欧美视频在线观看| 国产成人精品久久久久久| 视频区图区小说| 久久6这里有精品| 91成人精品电影| 在线精品无人区一区二区三| 日韩制服骚丝袜av| 久久久久网色| 国产又色又爽无遮挡免| 成人毛片a级毛片在线播放| 大陆偷拍与自拍| 国产精品99久久久久久久久| 少妇精品久久久久久久| 女的被弄到高潮叫床怎么办| 久久久久久久大尺度免费视频| 日韩一区二区三区影片| 一级黄片播放器| 狠狠精品人妻久久久久久综合| 午夜视频国产福利| 欧美日韩一区二区视频在线观看视频在线| av免费观看日本| 一级毛片久久久久久久久女| 插阴视频在线观看视频| 国产免费视频播放在线视频| av免费观看日本| 国产精品成人在线| 久久久久精品久久久久真实原创| 男的添女的下面高潮视频| 久久人妻熟女aⅴ| 日韩,欧美,国产一区二区三区| 久久鲁丝午夜福利片| 亚洲性久久影院| 一级毛片久久久久久久久女| 激情五月婷婷亚洲| tube8黄色片| 久久亚洲国产成人精品v| 日本-黄色视频高清免费观看| 91午夜精品亚洲一区二区三区| 国产一区二区在线观看日韩| 日韩强制内射视频| 国产片特级美女逼逼视频| 亚洲欧美中文字幕日韩二区| 嫩草影院入口| 亚洲国产欧美在线一区| 精品国产露脸久久av麻豆| 国产成人精品福利久久| 内射极品少妇av片p| 成人毛片a级毛片在线播放| 边亲边吃奶的免费视频| 久久人妻熟女aⅴ| 日韩一本色道免费dvd| 在线天堂最新版资源| 永久免费av网站大全| 看非洲黑人一级黄片| 亚洲精品国产色婷婷电影| 国产淫片久久久久久久久| 国产在线一区二区三区精| 99精国产麻豆久久婷婷| 又粗又硬又长又爽又黄的视频| 热re99久久精品国产66热6| 大陆偷拍与自拍| 国国产精品蜜臀av免费| 99九九在线精品视频 | 久久精品国产亚洲网站| 国产日韩欧美视频二区| 国产中年淑女户外野战色| 丝袜脚勾引网站| 亚洲精品456在线播放app| 欧美97在线视频| 一级毛片黄色毛片免费观看视频| 美女国产视频在线观看| 欧美日韩国产mv在线观看视频| 亚洲图色成人| 亚洲av日韩在线播放| 18禁裸乳无遮挡动漫免费视频| av一本久久久久| 一区二区三区免费毛片| 最近手机中文字幕大全| 亚洲国产精品成人久久小说| 18禁动态无遮挡网站| av线在线观看网站| 久久久a久久爽久久v久久| 久久久久久久久久久丰满| 亚洲一区二区三区欧美精品| 中国美白少妇内射xxxbb| 国产亚洲欧美精品永久| 欧美性感艳星| 人妻系列 视频| 97超碰精品成人国产| 最近中文字幕2019免费版| 人人妻人人澡人人看| 久久影院123| 国产成人aa在线观看| 男女无遮挡免费网站观看| 亚洲欧美日韩东京热| 久久这里有精品视频免费| 国产一区亚洲一区在线观看| 日本wwww免费看| 精品久久久久久久久亚洲| 久久久久久久久大av| 日韩大片免费观看网站| 精品久久久久久电影网| 精品一区在线观看国产| 视频中文字幕在线观看| av在线老鸭窝| 在线精品无人区一区二区三| 我的老师免费观看完整版| 亚洲欧美精品专区久久| 黑人高潮一二区| 一级a做视频免费观看| 国精品久久久久久国模美| 99热6这里只有精品| 极品少妇高潮喷水抽搐| 在线精品无人区一区二区三| 欧美精品一区二区大全| a级毛片免费高清观看在线播放| 一二三四中文在线观看免费高清| 如何舔出高潮| 亚洲人成网站在线观看播放| 性高湖久久久久久久久免费观看| 中文天堂在线官网| 国产 一区精品| 久久久亚洲精品成人影院| 国产欧美亚洲国产| 国产伦理片在线播放av一区| 中文资源天堂在线| 国产有黄有色有爽视频| 天堂中文最新版在线下载| 国产一区有黄有色的免费视频| h视频一区二区三区| 久久婷婷青草| 99视频精品全部免费 在线| 久热这里只有精品99| 亚洲成人av在线免费| 婷婷色综合大香蕉| 毛片一级片免费看久久久久| av黄色大香蕉| 99re6热这里在线精品视频| 九草在线视频观看| 天天躁夜夜躁狠狠久久av| 爱豆传媒免费全集在线观看| 亚洲av电影在线观看一区二区三区| 热re99久久国产66热| 久久精品国产鲁丝片午夜精品| 国产成人freesex在线| 精品国产一区二区久久| 亚洲av电影在线观看一区二区三区| 嫩草影院入口| 国产日韩欧美在线精品| 欧美日韩一区二区视频在线观看视频在线| 少妇丰满av| 高清欧美精品videossex| 亚洲国产欧美在线一区| 观看av在线不卡| 欧美日韩在线观看h| av福利片在线| 日韩精品有码人妻一区| 亚洲欧洲精品一区二区精品久久久 | 色吧在线观看| 亚洲一级一片aⅴ在线观看| 国产精品国产三级国产av玫瑰| 午夜老司机福利剧场| 最近2019中文字幕mv第一页| 日韩 亚洲 欧美在线| 偷拍熟女少妇极品色| 高清av免费在线| 日韩成人av中文字幕在线观看| 国产高清有码在线观看视频| 免费在线观看成人毛片| 亚洲精品第二区| 成人美女网站在线观看视频| 永久免费av网站大全| 亚洲精品乱久久久久久| 亚洲四区av| 国产日韩欧美亚洲二区| 91精品一卡2卡3卡4卡| 免费人妻精品一区二区三区视频| 免费播放大片免费观看视频在线观看| 亚洲av福利一区| 女性生殖器流出的白浆| 欧美+日韩+精品| 国产一区二区三区综合在线观看 | 大话2 男鬼变身卡| 亚洲中文av在线| 99精国产麻豆久久婷婷| 色婷婷av一区二区三区视频| 久久久久久伊人网av| 久久99热这里只频精品6学生| 午夜福利影视在线免费观看| 亚洲欧美日韩卡通动漫| 各种免费的搞黄视频| 免费观看无遮挡的男女| 丁香六月天网| 日韩成人伦理影院| 久久狼人影院| 一个人免费看片子| 99久久精品一区二区三区| 亚洲无线观看免费| 亚洲精品视频女| 欧美日韩综合久久久久久| 天天躁夜夜躁狠狠久久av| 人妻夜夜爽99麻豆av| 日韩av免费高清视频| 美女福利国产在线| 色婷婷久久久亚洲欧美| 性色av一级| 亚洲国产精品国产精品| 欧美老熟妇乱子伦牲交| .国产精品久久| 天天躁夜夜躁狠狠久久av| 亚洲欧美成人综合另类久久久| 午夜激情久久久久久久| 人人澡人人妻人| 51国产日韩欧美| 国产深夜福利视频在线观看| 国产日韩欧美亚洲二区| kizo精华| 成年av动漫网址| 亚洲图色成人| 观看免费一级毛片| 久久久亚洲精品成人影院| 欧美日韩精品成人综合77777| 日本wwww免费看| 国产日韩一区二区三区精品不卡 | 人妻系列 视频| 精品人妻一区二区三区麻豆| 另类亚洲欧美激情| a级毛片在线看网站| 欧美bdsm另类| 精品午夜福利在线看| 大片免费播放器 马上看| 国产精品三级大全| av卡一久久| 国产中年淑女户外野战色| 精品人妻一区二区三区麻豆| a级片在线免费高清观看视频| 五月伊人婷婷丁香| 中文字幕av电影在线播放| 美女中出高潮动态图| 国产精品一区二区性色av| av国产精品久久久久影院| 黑人巨大精品欧美一区二区蜜桃 | 亚洲丝袜综合中文字幕| 免费播放大片免费观看视频在线观看| 久久久久久人妻| 麻豆成人av视频| av福利片在线| 亚洲色图综合在线观看| 久久精品国产亚洲网站| 又爽又黄a免费视频| 男人添女人高潮全过程视频| 久久99精品国语久久久| 亚洲av电影在线观看一区二区三区| 午夜91福利影院| 国产精品国产av在线观看| 欧美日韩一区二区视频在线观看视频在线| 国国产精品蜜臀av免费| 国产精品一二三区在线看| 看十八女毛片水多多多| 少妇丰满av| 日韩亚洲欧美综合| 久久久国产精品麻豆| 我的女老师完整版在线观看| 国产美女午夜福利| 中文字幕制服av| 欧美老熟妇乱子伦牲交| www.av在线官网国产| 又大又黄又爽视频免费| 国产精品久久久久久精品电影小说| 日本91视频免费播放| 一区二区三区精品91| 国产视频首页在线观看| 久久人人爽人人片av| 久久国产精品大桥未久av | 久久精品国产亚洲网站| 免费观看在线日韩| 欧美日韩av久久| 日本黄色日本黄色录像| 狂野欧美激情性xxxx在线观看| 熟女av电影| 亚洲国产欧美在线一区| 人人妻人人添人人爽欧美一区卜| 天美传媒精品一区二区| 三上悠亚av全集在线观看 | 水蜜桃什么品种好| 国产淫语在线视频| 亚洲怡红院男人天堂| 亚洲成人一二三区av| 观看av在线不卡| 亚洲精品一区蜜桃| 亚洲av男天堂| 日韩电影二区| 久久午夜综合久久蜜桃| 精品亚洲成国产av| 街头女战士在线观看网站| 国产精品久久久久久av不卡| 国产欧美日韩一区二区三区在线 | 日本av手机在线免费观看| 久久久久久久国产电影| 不卡视频在线观看欧美| 国产毛片在线视频| 超碰97精品在线观看| 亚洲无线观看免费| 免费人妻精品一区二区三区视频| 男女国产视频网站| 久久久久网色| 精品一区二区三卡| 国内揄拍国产精品人妻在线| 狂野欧美白嫩少妇大欣赏| 99国产精品免费福利视频| 免费久久久久久久精品成人欧美视频 | 麻豆乱淫一区二区| 午夜福利网站1000一区二区三区| 激情五月婷婷亚洲| videos熟女内射| 亚洲国产精品一区二区三区在线| 综合色丁香网| 久热久热在线精品观看| 国产伦精品一区二区三区四那| 我的女老师完整版在线观看| 18禁在线播放成人免费| 韩国高清视频一区二区三区| 人体艺术视频欧美日本| 久久久久视频综合| 欧美性感艳星| 久久精品熟女亚洲av麻豆精品| 久热久热在线精品观看| 三级经典国产精品| 最近最新中文字幕免费大全7| 久久久久久久久久久丰满| 成人黄色视频免费在线看| 精品一区二区三区视频在线| 大又大粗又爽又黄少妇毛片口| 人人妻人人澡人人看| 国产精品国产三级国产专区5o| 国产白丝娇喘喷水9色精品| 黄色怎么调成土黄色| 免费看av在线观看网站| 久久午夜福利片| 一级av片app| 国产成人精品久久久久久| 国产av国产精品国产| 观看免费一级毛片| 精品国产一区二区三区久久久樱花| av在线播放精品| 夜夜爽夜夜爽视频| 一级二级三级毛片免费看| 精品亚洲成a人片在线观看| 日韩中文字幕视频在线看片| 日本91视频免费播放| 国产女主播在线喷水免费视频网站| 成人18禁高潮啪啪吃奶动态图 | 日产精品乱码卡一卡2卡三| 久久6这里有精品| 中文字幕制服av| 日韩av在线免费看完整版不卡| 亚洲国产精品一区三区| 中文字幕久久专区| 亚洲精品亚洲一区二区| 丰满少妇做爰视频| 久久亚洲国产成人精品v| 男女边吃奶边做爰视频| 极品教师在线视频| 在线观看av片永久免费下载| 国产综合精华液| 青春草国产在线视频| 亚洲第一区二区三区不卡| a级毛片免费高清观看在线播放| av黄色大香蕉| 国产av码专区亚洲av| 久久综合国产亚洲精品| 九九爱精品视频在线观看| 成年人免费黄色播放视频 | 少妇高潮的动态图| 午夜av观看不卡| 国产成人免费无遮挡视频| 亚洲国产精品999| 看非洲黑人一级黄片| 一级黄片播放器| 人体艺术视频欧美日本| 黄色毛片三级朝国网站 | 精品国产露脸久久av麻豆| 精品少妇久久久久久888优播| 国产日韩欧美在线精品| 日本vs欧美在线观看视频 | 大陆偷拍与自拍| 国产男人的电影天堂91| 99热这里只有是精品在线观看| 欧美激情国产日韩精品一区| 2022亚洲国产成人精品| 简卡轻食公司| 免费少妇av软件| 精品一区二区免费观看| 久久久久久人妻| 人妻夜夜爽99麻豆av| 曰老女人黄片| 美女大奶头黄色视频| 嘟嘟电影网在线观看| 一区在线观看完整版| 久久狼人影院| 人人妻人人添人人爽欧美一区卜| 久久久欧美国产精品| 亚洲欧洲精品一区二区精品久久久 | 成年av动漫网址| 久久国产精品大桥未久av | 99久国产av精品国产电影| 伦理电影大哥的女人| 久久精品久久久久久噜噜老黄| 亚洲精品色激情综合| 久久99热这里只频精品6学生| 久久精品国产鲁丝片午夜精品| 69精品国产乱码久久久| 精品少妇久久久久久888优播| 午夜视频国产福利| 在线观看免费日韩欧美大片 | 国产精品久久久久久久电影| 色婷婷av一区二区三区视频| 欧美xxxx性猛交bbbb| 人人妻人人爽人人添夜夜欢视频 | 丝袜在线中文字幕| 亚州av有码| 国产色爽女视频免费观看| 日日撸夜夜添| 丝袜在线中文字幕| 国产熟女午夜一区二区三区 | 日韩中字成人| av又黄又爽大尺度在线免费看| 亚洲国产精品国产精品| 两个人的视频大全免费| 婷婷色综合www| 欧美老熟妇乱子伦牲交| 日韩一本色道免费dvd| 日韩成人av中文字幕在线观看| 中文在线观看免费www的网站| 菩萨蛮人人尽说江南好唐韦庄| 国产精品久久久久久精品电影小说| 精品熟女少妇av免费看| 91精品伊人久久大香线蕉| 成人亚洲精品一区在线观看| 一级片'在线观看视频| av专区在线播放| 精品一区二区免费观看| 免费少妇av软件| 色婷婷av一区二区三区视频| 丁香六月天网| 黄片无遮挡物在线观看| 男人爽女人下面视频在线观看| 少妇高潮的动态图| 免费av中文字幕在线| 18+在线观看网站| 国产 一区精品| 精品99又大又爽又粗少妇毛片| 男的添女的下面高潮视频| 男人爽女人下面视频在线观看| av有码第一页| 精品久久久精品久久久| 日本黄色日本黄色录像| 六月丁香七月| 99re6热这里在线精品视频| 久久综合国产亚洲精品|