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

    長白山源頭區(qū)地下水質評價及監(jiān)測指標優(yōu)化

    2023-10-26 04:03:28谷志琪卞建民馬麗欣孫曉慶阮冬梅
    中國環(huán)境科學 2023年10期
    關鍵詞:正態(tài)水質評價水質

    谷志琪,卞建民,王 宇*,馬麗欣,孫曉慶,阮冬梅

    長白山源頭區(qū)地下水質評價及監(jiān)測指標優(yōu)化

    谷志琪1,卞建民1,王 宇1*,馬麗欣2,孫曉慶1,阮冬梅1

    (1.吉林大學新能源與環(huán)境學院,地下水資源與環(huán)境教育部重點實驗室,吉林 長春 130021;2.吉林省生態(tài)環(huán)境監(jiān)測中心,吉林 長春 130021)

    以長白山源頭區(qū)-安圖縣作為研究區(qū),分析地下水化學特征及其形成過程,基于組合熵權的云模型評估地下水質量,并利用隨機森林聯(lián)合多元線性回歸構建水質指標優(yōu)化模型,確定地下水源保障監(jiān)測的關鍵指標.結果表明,研究區(qū)地下水水化學類型主要為HCO3-Mg·Na·Ca型,主要受巖石風化作用影響;地下水質量評價等級為Ⅰ~Ⅲ類的樣品數(shù)達到74%;最佳指標優(yōu)化模型的2和RMSE值分別為0.6333和0.726,優(yōu)化指標為F-、Na+、TDS、Cl-,作為該區(qū)地下水監(jiān)測的關鍵指標,能夠有效減少監(jiān)測費用,并為強化水源地安全保障提供科學依據(jù).

    長白山優(yōu)質水源;云模型;隨機森林;水質評價;指標優(yōu)化;地下水

    長白山天池水資源豐富,是該地區(qū)的重要水源[1],隨著經(jīng)濟的發(fā)展,人們對優(yōu)質水源的需求愈加迫切,對水資源的開發(fā)利用程度也逐漸增加,而水資源大量開發(fā)過程中容易對水源地水質造成影響[2].因此,加強優(yōu)質水源地保護,探究區(qū)域地下水水質現(xiàn)狀,構建合理的地下水監(jiān)測體系等工作亟待開展,對于保障優(yōu)質水資源的可持續(xù)利用,支撐社會經(jīng)濟和生態(tài)環(huán)境的協(xié)調發(fā)展具有重要意義.

    目前我國針對優(yōu)質水源地已經(jīng)開展了水質適宜性分析[3]、構建水源評價指標體系[4]和決策供水水源模型[5]等一系列工作,但是部分水源地水質狀況仍不明確,沒有建立可靠的監(jiān)測保障體系.地下水質量評價涉及多個指標的不同決策,大量使用的Nemerow指數(shù)法[6]、模糊綜合評價法[7]、支持向量機[8]、主成分分析[9]和投影尋蹤法[10]等方法由于適用條件及評估標準的不同,當指標監(jiān)測值在水質分類標準閾值附近時具有一定的模糊性和隨機性[11],難以做出準確判斷進而導致評價結果的不確定性.針對這種不確定性,基于概率統(tǒng)計和模糊數(shù)學的正態(tài)云模型在定性概念與定量數(shù)值相互轉換過程中具有優(yōu)越性[12],還可以考慮評價標準的非線性和評價參數(shù)的多樣性,從而確保評價結果更加可靠.

    基于不同類別的大量地下水監(jiān)測數(shù)據(jù)信息基本可以反映出當前環(huán)境下的水質狀況,但不同的監(jiān)測指標與水環(huán)境之間存在不匹配問題,長期處于符合標準范圍內的非重要指標容易減弱重要指標對水環(huán)境的影響[13],使水質評價結果出現(xiàn)一定的偏差,難以表征真實的水環(huán)境狀況,還會造成時間和人力物力上的浪費.目前對水環(huán)境監(jiān)測指標進行優(yōu)化的研究往往結合水質評價結果進行主觀篩選[14],依據(jù)不同指標間的相關性[15]以及構建水質指數(shù)模型[16]優(yōu)化關鍵指標等.隨機森林方法不但數(shù)據(jù)處理速度快,較其他同類算法準確率也更高,已經(jīng)在地下水潛力評估[17]、地下水水質評價[18]、地下水污染預測[19]等領域進行了應用[20],而根據(jù)隨機森林方法定量評價監(jiān)測指標對于水質類別的重要性,構建指標優(yōu)化模型的研究鮮見報道.因此,本文選取長白山源頭區(qū)安圖縣作為研究對象,進行地下水調查和樣品采集測試,分析研究區(qū)地下水水化學特征及其形成過程,在此基礎上結合熵權理論運用正態(tài)云模型評價地下水水質,驗證該方法可靠性的同時將其評價結果與隨機森林方法聯(lián)用,構建地下水監(jiān)測指標的優(yōu)化模型,篩選影響該地區(qū)地下水水質的關鍵指標,為了解區(qū)域地下水水質狀況和保護優(yōu)質水源地提供理論依據(jù).

    1 數(shù)據(jù)與方法

    1.1 研究區(qū)概況

    長白山地處我國東北部(東經(jīng)127°~132°,北緯39°~46°),是中國濕潤地區(qū)最大的火山臺地山脈,面積約19600km2,屬于東亞季風氣候區(qū),境內多年平均氣溫3.54℃,年平均降水量在 700~1400mm之間[21].區(qū)內總體地勢東南高西北低,地質構造復雜,以天池為中心形成了眾多環(huán)狀斷裂及放射狀斷裂,各種構造斷裂相互穿插、交切,為地下水提供了大量的循環(huán)及儲存空間[22].松花江、鴨綠江、圖們江三大江皆發(fā)源于長白山頂端天池,在流經(jīng)途中分叉形成眾多細小分散支流,形成流域面積超過20km2.

    以長白山北側源頭區(qū)—安圖縣(127°48¢~ 129°08¢E,北緯42°01¢~43°24¢N)作為研究區(qū),該區(qū)域總面積約為7400km2,總人口約22萬,植被覆蓋率達87%以上,總體地勢南高北低,東高西低,呈階梯狀臺地展布,多年平均氣溫約為3.8℃,多年平均降水量約為752mm.區(qū)內地下水主要為玄武巖孔洞裂隙水,主要賦存于軍艦山組玄武巖地層中,其氣孔率一般為20%~30%,該含水層的補給來源包括大氣降水和凝結水[23],通過在山前谷地富集,排泄成泉,成為地表徑流的補給源.地表徑流中二道白河直接發(fā)源于長白山天池,與頭道白河、三道白河等均屬于松花江水系.

    圖1 研究區(qū)位置及采樣點分布示意

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

    利用本課題組在研究區(qū)采集的35個點位的地下水樣品測試數(shù)據(jù),采樣點的布設結合了研究區(qū)地質及水文地質條件,并沿地下水流動路徑布置(圖1).樣品的采集、運輸和儲存工作嚴格遵守《地下水環(huán)境監(jiān)測技術規(guī)范》(HJ/T 164-2004),測試工作由吉林大學測試科學實驗中心承擔.檢測指標14項,測試方法見表1.

    表1 地下水水化學組分測定方法

    1.3 研究方法

    1.3.1 正態(tài)云理論 正態(tài)分布是概率論中最重要的分布之一,根據(jù)正態(tài)分布曲線延伸出了鐘形隸屬函數(shù),在其基礎上發(fā)展了用于處理定性概念與定量描述的不確定轉換模型——正態(tài)云模型[24],其不確定性使用(、、)等數(shù)字特征進行確定,其計算方程如下:

    (1):期望,云滴在論域空間分布的期望,代表云滴的定性概念.

    式中:min、max分別為評價因子隸屬于某一水質等級的標準上下限值.

    (2)En:熵,反映了代表定性概念云滴的離散程度和在論域空間可被概念接受的云滴的取值范圍.

    (3)He:超熵,反映了云的離散程度.

    式中:為常數(shù),根據(jù)系列的分布情況調整[25].耦合熵權法的正態(tài)云模型的流程如圖2所示.

    圖2 耦合熵權的正態(tài)云模型流程

    1.3.2 隨機森林方法 隨機森林(RF)利用裝袋重抽樣和特征隨機性來構建每棵單獨的決策樹,生成不相關的森林,其對每個決策結果進行投票表決的最終結果比任何單獨的樹更準確[26].該方法的一個突出特點是可以通過計算Mean Decreased Gini Coefficient(平均下降Gini指數(shù),MDG)來評價每個因素的相對重要性及其對預測輸出的影響[27].MDG是數(shù)據(jù)集中隨機索引被錯誤劃分的概率,用于指示節(jié)點的純度,MDG值越大,該指標的重要性就越大,其計算公式如下:

    式中:(i,j)表示第個變量的第個指數(shù);表示回歸樹的節(jié)點編號.當節(jié)點分成兩個子節(jié)點時,兩個子節(jié)點處的基尼指數(shù)將趨于最小化,并將用作節(jié)點處基尼指數(shù)的最終值:

    當節(jié)點被拆分為兩個子節(jié)點時,其子節(jié)點的基尼指數(shù)之和低于其父節(jié)點的基尼系數(shù)之和,并且減少的絕對值計算如下:

    式中:Gini(l)是左側子節(jié)點的Gini指數(shù); Gini(r)是右側子節(jié)點的基尼指數(shù);l和r分別是左側和右側節(jié)點的樣本大小.

    2 結果與討論

    2.1 地下水化學特征及其形成過程

    根據(jù)研究區(qū)地下水樣品指標檢出率情況,結合《地下水質量標準》(GB/T14848—2017)[28]選取14個水質指標進行描述性統(tǒng)計分析,并繪制Piper三線圖判別水化學類型.結果表明,研究區(qū)地下水普遍呈中性,pH值變化范圍為6.16~8;地下水以淡軟水類型為主,TDS變化范圍為96~2029.68mg/L, TH變化范圍為24.09~895.41mg/L,分別有約14.29%和11.43%的樣品超過了地下水Ⅲ類水閾值;其他指標除F-和TFe外,均未超過地下水Ⅲ類水閾值,其中F-在所有指標中超標率最高,為22.86%.經(jīng)分析,該地區(qū)F-含量較高的原因是受到地下水深循環(huán)影響,火山活動形成的大型斷裂中高溫熱水溶解花崗巖中的氟元素形成高氟水,補給到淺層地下水中使其中F-含量相對較高,與前人的研究結論吻合[29].TFe的超標率為14.28%,主要受到長白山噴發(fā)的火成巖鐵元素本底值高的影響.pH值的變異系數(shù)非常小,表明區(qū)內地下水pH值總體上分布均勻,除K+、F-和H2SiO3以外其余指標變異系數(shù)均大于1,在研究區(qū)范圍內空間分布存在一定的差異,TFe的變異系數(shù)最大,其形成受火山地質作用影響是其空間差異性的主要原因.研究區(qū)地下水主要陽離子是Na+和Ca2+,主要陰離子是HCO3-,地下水水化學類型主要為HCO3- Mg·Na·Ca型和HCO3-Na·Mg·Ca型(圖3).

    圖3 研究區(qū)地下水piper三線圖

    圖4 研究區(qū)地下水Gibbs圖

    基于水化學數(shù)據(jù)繪制Gibbs圖,根據(jù)各點在 Gibbs圖中的位置,分析溶解于水中的化學成分的主要來源,判斷地下水的水化學成因.如圖4所示,大部分水樣點均落在Gibbs圖中間的巖石風化控制區(qū),少部分水樣點分布在蒸發(fā)濃縮作用控制區(qū),而在大氣降雨作用控制區(qū)沒有水樣點的分布.表明地下水化學成分的形成以巖石風化作用的影響為主,蒸發(fā)濃縮作用的影響為次,而大氣降雨作用幾乎沒有影響.個別水樣點落在Gibbs圖外,這表明還可能受到陽離子交換作用的影響.

    2.2 地下水水質評價

    考慮到地下水質量標準中K+、Ca2+、Mg2+、HCO3-、H2SiO3和pH值等指標沒有具體類別對應的標準閾值,故選擇Na+、SO42-、Cl-、TDS、TH、TFe、F-、和NO3-這8個指標,通過正態(tài)云模型理論耦合熵權法評價安圖縣地下水水質,了解地下水環(huán)境質量.

    表2 各指標權重計算結果

    圖5 水質評價結果對比圖

    不同指標的熵權計算結果如表2所示,在確定各指標權重的基礎上,結合云模型參數(shù)得出水質綜合評價結果,并與內梅羅綜合計算結果進行比較(圖5).由圖5可知:研究區(qū)內地下水水質較好,大部分均為Ⅰ~Ⅲ類,基于組合熵權的正態(tài)云模型評價結果中等于或好于Ⅲ類水的樣品有26個,內梅羅評價結果中等于或好于Ⅲ類水的樣品有22個,兩種方法對于區(qū)域水質的綜合評價結果具有較好的一致性.內梅羅綜合評價結果中絕大部分樣品點屬于Ⅰ類或Ⅴ等邊緣水質類別,占比超過全部樣品的66%,而基于組合熵權的正態(tài)云模型評價結果中各級別樣品數(shù)量較為平均,不同類別樣品占比均小于35%.

    這是因為內梅羅綜合評價根據(jù)計算F值進行分類,評價結果易受極大值影響,忽略水質較為適中的類別,而在大量云滴基礎上建立的正態(tài)云模型減少了極大值對評價結果的影響程度,同時與熵權法耦合,降低了評價過程中的模糊性,使得評價結果更能反映出地下水環(huán)境的真實情況.云模型結合不同權重方法應用于評價某地區(qū)的地下水水質的相關研究均得出了相似結論[30-31],也證明了本文結論的可靠性.

    2.3 地下水監(jiān)測指標優(yōu)化

    基于組合熵權的正態(tài)云模型的地下水質量評價結果,通過隨機森林定量評價水質指標對水質類別的特征重要性,結果表明,F-是對研究區(qū)地下水水質影響最大的指標,MDG值為3.08,Na+、TH、TDS、Cl-與SO42-的重要性差別較小,均在2.05~2.3范圍內,TFe和NO3-的重要性相對較低(圖6(a)),因此,可選擇F-、Na+、TH、TDS、Cl-與SO42-等指標構建優(yōu)化指標模型.

    表3 優(yōu)化模型指標選擇及分析結果

    同時,為確定構建優(yōu)化模型的最佳指標數(shù)量,本文采取十折交叉驗證法,計算不同指標數(shù)量下水質類別的擬合誤差(圖6(b)).結果表明,水質擬合誤差隨著選取指標數(shù)量的增加先增加,后減小,當選擇4個指標時,水質擬合的誤差最小,之后在選取指標數(shù)量增加的情況下,擬合誤差也趨于穩(wěn)定,因此,選擇4個指標構建優(yōu)化指標模型.此外,由于F-的相對重要性排序顯著高于其他指標(圖6(a)),因此,以F-作為優(yōu)化模型的基礎指標,分別加入Na+、TH、TDS、Cl-與SO42-構建10個優(yōu)化指標模型,通過多元線性回歸分析對地下水質量進行擬合,比較模擬結果與組合熵權的云模型評價結果間的擬合優(yōu)度(R2和RMSE值),判斷優(yōu)化模型的可行性.

    各模型選擇的指標及模型的2、RMSE值如表3所示.結果表明,優(yōu)化模型3、4和5的2值顯著高于其他模型,分別為0.6312、0.6333、0.6344,且RMSE值也低于其他模型,分別為0.7393、0.7260、0.7371,其中,模型4的模擬誤差較模型5少1.53%,且擬合度較模型5僅低0.17%.另外,從圖7中可以看出,優(yōu)化模型4、5對Ⅱ、Ⅲ類水擬合效果均十分準確,而模型4對Ⅰ、Ⅳ類水的模擬效果較模型5更好,且由相似的模型7、8也可證明,選擇Cl-的模型較選擇SO42-的模型對Ⅰ、Ⅳ類水的模擬效果更準確.

    由2.1節(jié)可知,F-、TDS是研究區(qū)內地下水超標率最高的兩個指標,可代表研究區(qū)內水質較差區(qū)域的主要監(jiān)測指標.Na+是地下水中陽離子的最主要的成分,與Cl-、SO42-共同作為研究區(qū)內未超標組分,可以表征高質量的地下水,通過前文分析可知,選擇Cl-的模型較選擇SO42-的模型擬合效果更好.因此,包含F(xiàn)-、Na+、TDS、Cl-這4個關鍵水質指標的優(yōu)化模型4可作為研究區(qū)監(jiān)測指標的最佳優(yōu)化模型,以上優(yōu)化指標可作為研究區(qū)地下水水質監(jiān)測及評價主要參考指標.

    3 結論

    3.1 研究區(qū)地下水普遍呈中性,以軟淡水類型為主,水化學類型主要為HCO3-Mg·Na·Ca型,超標組分主要包括F-、TFe、TDS和TH,地下水中主要離子來源于含鹽礦物的溶解,水化學特征主要受巖石風化作用控制.

    3.2 區(qū)內約43%的地下水采樣點水質為Ⅰ、Ⅱ類水,基于組合熵權的正態(tài)云模型評價結果與內梅羅綜合評價法對區(qū)域水質的整體判斷具有較好的一致性,而組合熵權的云模型對水質類別的劃分更為細致.

    3.3 基于隨機森林計算的相對重要性排序為F-、Na+、TH、TDS、Cl-、SO42-、TFe、NO3-,構建優(yōu)化模型的最佳指標數(shù)量為4個,優(yōu)化的關鍵指標為F-、Na+、TDS、Cl-,可作為研究區(qū)地下水水質監(jiān)測及評價主要參考指標,使得減少監(jiān)測費用的同時能夠有效保障優(yōu)質水源地安全.

    [1] 呂 琳.吉林省礦泉水區(qū)域與非礦泉水區(qū)域環(huán)境中礦物質元素含量分析 [D]. 長春:吉林大學, 2015. Lu L. Research on the environmental mineral elements contents of the mineral water area and the non-mineral water area in Jilin Provinoe [D]. Changchun: Jilin University, 2015.

    [2] 汪明武,周天龍,葉 暉,等.基于聯(lián)系云的地下水水質可拓評價模型 [J]. 中國環(huán)境科學, 2018,38(8):3035-3041. Wang M W, Zhou T L, Ye H, et al. A novel extension evaluation model of groundwater quality based on connection cloud model [J]. China Environmental Science, 2018,38(8):3035-3041.

    [3] 朱 琴,左麗明,彭錦添,等.河北省近岸海域海水淡化取水水質適宜性分析 [J]. 海洋開發(fā)與管理, 2017,34(7):60-66. Zhu Q, Zuo L M, Peng J T, et al. Analysis on optimal water quality sources for seawater desalination in Hebei Coastal Waters [J]. Ocean Development and Management, 2017,34(7):60-66.

    [4] 沈 強,胡 俊,胡菊香.浙江省綠色水源地評價體系研究 [J]. 環(huán)境工程, 2014,32(1):847-851,859. Shen Q, Hu J, Hu J X. Study on assessment system for green water-source in Zhejiang province [J]. Environmental Engineering, 2014,32(1):847-851,859.

    [5] 于 冰,梁國華,何 斌,等.城市供水系統(tǒng)多水源聯(lián)合調度模型及應用 [J]. 水科學進展, 2015,26(6):874-884. Yu B, Liang G H, He B, et al. Modeling of joint operation for urban water-supply system with multi-water sources and its application [J]. Advances in Water Science, 2015,26(6):874-884.

    [6] Yang Q, Zhang J, Hou Z, et al. Shallow groundwater quality assessment: Use of the improved Nemerow pollution index, wavelet transform and neural networks [J]. Journal of Hydroinformatics, 2017, 19(5):784-794.

    [7] Wang A F, Yang X T, Gu X B. The risk assessment of rockburst intensity in the highway tunnel based on the variable fuzzy sets theory [J]. Scientific reports, 2023,13(1):4755.

    [8] Gu Z Q, Bian J M, Wu J J, et al. Effects of anthropogenic activities on hydrochemical characteristics of ground water of Da'an irrigation area in Western of Jilin Province [J]. Environ Sci Pollut Res Int, 2021, 29(14):20479-20495.

    [9] 劉 瀟,薛 瑩,紀毓鵬,等.基于主成分分析法的黃河口及其鄰近水域水質評價[J]. 中國環(huán)境科學, 2015,35(10):3187-3192. Liu X, Xue Y, Ji Y P, et al. An assessment of water quality in the Yellow River estuary and its adjacent waters based on principal component analysis [J]. China Environmental Science, 2015,35(10): 3187-3192.

    [10] Jia Z, Bian J M, Wang Y, et al. Assessment and validation of groundwater vulnerability to nitrate in porous aquifers based on a DRASTIC method modified by projection pursuit dynamic clustering model [J]. Journal of Contaminant Hydrology, 2019,226:103522.

    [11] 高玉琴,賴麗娟,姚 敏,等.基于正態(tài)云-模糊可變耦合模型的水環(huán)境質量評價 [J]. 水資源與水工程學報, 2018,29(5):1-7. Gao Y Q, Lai L J, Yao M, et al. Water environment quality assessment based on normal cloud-fuzzy variable coupling model [J]. Journal of Water Resources and Water Engineering, 2018,29(5):1-7.

    [12] 李德毅,孟海軍,史雪梅.隸屬云和隸屬云發(fā)生器 [J]. 計算機研究與發(fā)展, 1995:15-20. Li D Y, Meng H J, Shi X M. Membership clouds and membership cloud generators [J]. Journal of Computer Research and Development, 1995:15-20.

    [13] 姜厚竹.松花江流域省界緩沖區(qū)水質監(jiān)測指標與斷面優(yōu)化 [D]. 哈爾濱:東北林業(yè)大學, 2017. Jiang H Z. Water quality monitoring index and section optimization of provincial buffer zone in Songhua River basin [D]. Harbin: Northeast Forestry University, 2017.

    [14] 曲茉莉.黑龍江水系干流監(jiān)測水質指標與斷面優(yōu)化研究 [D]. 哈爾濱:哈爾濱工業(yè)大學, 2012. Qu M L. Study on monitoring indexes and monitoring sections optimization of the main stream of Heilongjiang River [D]. Harbin: Harbin Institute of Technology, 2012.

    [15] 侯佳均.北京市平谷區(qū)地下水污染監(jiān)控指標優(yōu)化研究[D]. 成都:成都理工大學, 2020. Hou J J. Study on optimization of groundwater pollution monitoring index in Pinggu District of Beijing [D]. Chengdu: Chengdu University of Technolog, 2020.

    [16] Naghibi Seyed A K, Daneshi Alireza. Application of support vector machine, random forest, and genetic algorithm optimized random forest models in groundwater potential mapping [J]. Water Resources Management, 2017,31(9):2761-2775.

    [17] Naghibi Seyed A K, Daneshi Alireza. Application of support vector machine, random forest, and genetic algorithm optimized random forest models in groundwater potential mapping [J]. Water Resources Management, 2017,31(9):2761-2775.

    [18] 許力文,卞建民,孫曉慶,等.灌區(qū)退水對區(qū)域地下水質影響與健康風險評估 [J/OL]. 中國環(huán)境科學:1-10[2023-02-22].DOI:10.19674/ j.cnki.issn1000-6923.20221117.017. Xu L W, Bian J M, Sun X Q, et al. The influence and health risk assessment of groundwater quality in irrigated area [J/OL]. China Environmental Science: 1-10[2023-02-22].DOI:10.19674/j.cnki. issn1000-6923.20221117.017.

    [19] He S, W J H, Wang D, et al. Predictive modeling of groundwater nitrate pollution and evaluating its main impact factors using random forest [J]. Chemosphere, 2022,290.

    [20] 杜尚海,古成科,張文靜.隨機森林理論及其在水文地質領域的研究進展 [J]. 中國環(huán)境科學, 2022,42(9):4285-4295. Du S H, Gu C K, Zhang W J. A review on the progresses in random forests theory and its applications in hydrogeology [J]. China Environmental Science, 2022,42(9):4285-4295.

    [21] Xu T T, Zheng J Q, Han S J, et al. Responses of soil nitrogen transformation to long-term nitrogen fertilization and precipitation changes in a broad-leaved Korean pine forest in Changbai Mountains, China [J]. The journal of applied ecology, 2018,29(9):2797-2807.

    [22] 張文卿,王文鳳,劉淑芹,等.長白山礦泉水補給徑流與排泄關系 [J]. 河海大學學報(自然科學版), 2019,47(2):108-113. Zhang W Q, Wang W F, Liu S Q, et al. Relationship of recharge runoff and drainage for the mineral water in the Changbai Mountain [J] Journal of Hohai University (Natural Sciences), 2019,47(2):108-113.

    [23] 馬于曦.安圖縣礦泉水水源涵養(yǎng)能力評估與生態(tài)紅線劃定 [D]. 長春:吉林大學, 2021. Ma Y X. Assessment of water conservation capacity and delineation of ecological red line of mineral water in Antu County [D]. Changchun: Jilin University, 2021.

    [24] Liu, Q Y, Wang M W, Zhou T L, et al. A connection cloud model coupled with extenics for water eutrophication evaluation [J]. Earth Science Informatics, 2019,12(4):659-669.

    [25] 雷麗萍.基于云模型的水質評價方法優(yōu)化研究 [D]. 成都:南交通大學, 2019. Lei L P. Research on optimization of water quality assessment method based on cloud model [D]. Chengdu: Southwest Jiaotong University, 2019.

    [26] 楊明悅,毛獻忠.基于變量重要性評分-隨機森林的溶解氧預測模型——以深圳灣為例 [J]. 中國環(huán)境科學, 2022,42(8):3876-3881. Yang M Y, Mao X Z. Dissolved oxygen prediction model based on variable importance measures and random forest: A case study of Shenzhen Bay [J]. China Environmental Science, 2022,42(8):3876- 3381.

    [27] Cutler D R, Edwards Jr T C, Beard K H, et al. Random forests for classification in ecology [J]. Ecology, 2007,88(11):2783-2792.

    [28] GB/T 14848-2017 地下水質量標準 [S]. GB/T 14848-2017 Groundwater quality standard [S].

    [29] 曹玉和,齊佳偉,熊紹禮.吉林省氟中毒病區(qū)水文質地特征及防氟改水對策 [J]. 中國地質, 2010:690-695. Cao Y H, Qi J W, Xiong S L. Hydrogeological characteristics of endemic fluorine disease areas of Jilin Province and water source project countermeasures for fluorine prevention [J] Geology in China, 2010:690-695.

    [30] 侯炳江.基于組合賦權的水質綜合評價云模型及其應用 [J]. 水電能源科學, 2016,34(8):24-27. Hou B J. Water quality evaluation based on combination weight- normal cloud model and its application [J]. Water Resources and Power, 2016,34(8):24-27.

    [31] 趙夢婷,黃顯峰,金國裕,等.基于改進云模型的區(qū)域水質評價研究[J]. 安徽農業(yè)大學學報, 2020,47(5):778-783.Zhao M T, Huang X F, Jin G Y, et al. Study on regional water quality evaluation based on improved cloud model [J]. Journal of Anhui Agricultural University, 2020,47(5):778-783.

    Groundwater quality assessment and index optimization of water quality monitoring in the water source area of Changbai Mountain.

    GU Zhi-qi1, BIAN Jian-min1, WANG Yu1*, MA Li-xin2, SUN Xiao-qing1, RUAN Dong-mei1

    (1.Key Laboratory of Groundwater Resources and Environment, Ministry of Education, College of New Energy and Environment, Jilin University, Changchun 130021, China;2.Jilin Ecological Environment Monitoring Center, Changchun 130021, China)., 2023,43(10):5257~5264

    This study chose Antu County, the water source area of Changbai Mountain as the study area. the hydrochemical characteristics and formation mechanism of groundwater in the study area was analyzed. The cloud model based on entropy weight was used to evaluate the groundwater quality. Furthermore, an optimization model of water quality index was constructed by coupling random forest and stepwise multiple linear regression analysis to determine the key indicators of groundwater source security monitoring. The results showed that the primary water chemistry type of groundwater in the study area was HCO3-Mg·Na·Ca type, which was mainly controlled by the rock weathering and dissolution effects. About 74% of groundwater samples were classified as Class Ⅰ-Ⅲ. The2and RMSE values of the best index optimization model were 0.6333 and 0.726, respectively. F-, Na+, TDS and Cl-were identified as the key indicators of groundwater quality monitoring in the study area. This optimized water quality index can effectively reduce monitoring costs and provide scientific basis for guaranteeing the safety of water sources.

    high-quality water source in Changbai Mountain;cloud model;random forest;water quality evaluation;index optimization;groundwater

    X824

    A

    1000-6923(2023)10-5257-08

    2023-02-24

    吉林省環(huán)保廳環(huán)境保護科研項目(吉環(huán)科字第2022-10號);吉林省科技廳重大科技專項(20230303007SF)

    * 責任作者, 副教授, chair19881224@126.com

    谷志琪(1997-),男,吉林長春人,吉林大學博士研究生,主要從事地下水資源評價與承載力研究.發(fā)表論文2篇.lh2568@163.com.

    谷志琪,卞建民,王 宇,等.長白山源頭區(qū)地下水質評價及監(jiān)測指標優(yōu)化 [J]. 中國環(huán)境科學, 2023,43(10):5257-5264.

    Gu Z Q, Bian J M, Wang Y, et al. Groundwater quality assessment and index optimization of water quality monitoring in the water source area of Changbai Mountain [J]. China Environmental Science, 2023,43(10):5257-5264.

    猜你喜歡
    正態(tài)水質評價水質
    水質抽檢豈容造假
    環(huán)境(2023年5期)2023-06-30 01:20:01
    一月冬棚養(yǎng)蝦常見水質渾濁,要如何解決?這9大原因及處理方法你要知曉
    滻灞河水質評價方法研究
    雙冪變換下正態(tài)線性回歸模型參數(shù)的假設檢驗
    基于概率統(tǒng)計和模糊綜合評價法的水質評價模型及其應用——以拉薩河水質評價為例
    基于SPAM的河流水質評價模型
    基于泛正態(tài)阻抗云的諧波發(fā)射水平估計
    半?yún)?shù)EV模型二階段估計的漸近正態(tài)性
    基于正態(tài)變換的貸款組合定價模型構建及實證
    水質總磷測定存在的問題初探
    河南科技(2014年23期)2014-02-27 14:19:07
    麻豆av在线久日| aaaaa片日本免费| 无人区码免费观看不卡| 欧美日韩亚洲国产一区二区在线观看| 久久热在线av| 真人做人爱边吃奶动态| 亚洲精品中文字幕在线视频| 欧美在线黄色| 亚洲熟妇熟女久久| 国产熟女午夜一区二区三区| 国产99久久九九免费精品| 97人妻天天添夜夜摸| 国产亚洲精品第一综合不卡| 欧美激情极品国产一区二区三区| tocl精华| 黄色怎么调成土黄色| 在线观看免费视频日本深夜| 精品国产一区二区三区四区第35| 国产三级黄色录像| 欧美乱妇无乱码| 午夜成年电影在线免费观看| 精品福利观看| 一区二区日韩欧美中文字幕| 老鸭窝网址在线观看| 桃色一区二区三区在线观看| 成人永久免费在线观看视频| 嫩草影视91久久| 久久中文看片网| 老鸭窝网址在线观看| 欧美一级毛片孕妇| 欧美日韩瑟瑟在线播放| 村上凉子中文字幕在线| 免费女性裸体啪啪无遮挡网站| 国产国语露脸激情在线看| 黄片播放在线免费| 一区二区三区激情视频| 亚洲第一青青草原| 一边摸一边抽搐一进一小说| 中文字幕最新亚洲高清| 大型黄色视频在线免费观看| 欧美激情 高清一区二区三区| 啦啦啦 在线观看视频| 亚洲精品中文字幕一二三四区| 午夜免费激情av| av福利片在线| 两性夫妻黄色片| 国产精华一区二区三区| www.www免费av| www.www免费av| 男人的好看免费观看在线视频 | 黄频高清免费视频| 国产亚洲精品久久久久5区| 国产av精品麻豆| 免费高清在线观看日韩| 后天国语完整版免费观看| 久久久久精品国产欧美久久久| 女性被躁到高潮视频| 久久久国产欧美日韩av| 51午夜福利影视在线观看| 夜夜夜夜夜久久久久| 日日爽夜夜爽网站| 日韩精品青青久久久久久| 女生性感内裤真人,穿戴方法视频| 一区二区三区精品91| 男女高潮啪啪啪动态图| 国产精品免费视频内射| 亚洲成人国产一区在线观看| 成年人黄色毛片网站| av天堂久久9| 九色亚洲精品在线播放| 18禁国产床啪视频网站| 欧美乱妇无乱码| 淫妇啪啪啪对白视频| 丰满饥渴人妻一区二区三| 亚洲激情在线av| 午夜福利在线观看吧| 国产成人啪精品午夜网站| 757午夜福利合集在线观看| 国产成人系列免费观看| 国产精品 国内视频| 欧美成人免费av一区二区三区| av电影中文网址| 悠悠久久av| 黑人巨大精品欧美一区二区蜜桃| 成人亚洲精品一区在线观看| 免费看十八禁软件| 成人亚洲精品av一区二区 | 一级片免费观看大全| 亚洲精品国产一区二区精华液| 午夜91福利影院| 久久久久久免费高清国产稀缺| 91九色精品人成在线观看| 一边摸一边抽搐一进一出视频| 亚洲精品在线美女| 亚洲中文字幕日韩| 波多野结衣av一区二区av| 日本精品一区二区三区蜜桃| 久久久久久久午夜电影 | 自拍欧美九色日韩亚洲蝌蚪91| 午夜成年电影在线免费观看| 我的亚洲天堂| 国产精品野战在线观看 | 90打野战视频偷拍视频| 一夜夜www| 中亚洲国语对白在线视频| 高清毛片免费观看视频网站 | 中文字幕av电影在线播放| 国产精品国产高清国产av| 国产一区在线观看成人免费| 亚洲五月婷婷丁香| 在线播放国产精品三级| 麻豆久久精品国产亚洲av | 亚洲精品国产一区二区精华液| 精品福利观看| 亚洲久久久国产精品| 久久九九热精品免费| 国产精品一区二区免费欧美| 日本免费一区二区三区高清不卡 | 国产国语露脸激情在线看| 一区二区三区国产精品乱码| 天天影视国产精品| 久久精品影院6| 交换朋友夫妻互换小说| 男女下面插进去视频免费观看| 亚洲九九香蕉| 久久久国产精品麻豆| 99国产极品粉嫩在线观看| 乱人伦中国视频| 婷婷六月久久综合丁香| 亚洲熟妇中文字幕五十中出 | 窝窝影院91人妻| 无限看片的www在线观看| 少妇 在线观看| 女人被躁到高潮嗷嗷叫费观| 一区福利在线观看| 99国产综合亚洲精品| 国产av又大| www.熟女人妻精品国产| 狂野欧美激情性xxxx| 男女下面插进去视频免费观看| 欧美丝袜亚洲另类 | 免费在线观看亚洲国产| 国产高清videossex| 午夜免费鲁丝| 窝窝影院91人妻| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲av第一区精品v没综合| 久久精品国产亚洲av高清一级| 中文字幕人妻熟女乱码| 国产亚洲欧美98| 国产精品永久免费网站| 一级毛片女人18水好多| 99久久99久久久精品蜜桃| 国产亚洲av高清不卡| 最好的美女福利视频网| 欧美日韩中文字幕国产精品一区二区三区 | 在线观看午夜福利视频| 看黄色毛片网站| 亚洲欧美日韩无卡精品| 中文欧美无线码| 淫妇啪啪啪对白视频| 人人妻,人人澡人人爽秒播| 久久中文字幕人妻熟女| 亚洲情色 制服丝袜| 成人国产一区最新在线观看| 99riav亚洲国产免费| 国产成人免费无遮挡视频| 国产在线观看jvid| 国产av又大| 又黄又粗又硬又大视频| 日韩欧美一区二区三区在线观看| 老司机午夜十八禁免费视频| 亚洲成人免费电影在线观看| 午夜亚洲福利在线播放| 50天的宝宝边吃奶边哭怎么回事| 国产又色又爽无遮挡免费看| 中文欧美无线码| 好看av亚洲va欧美ⅴa在| 欧美黑人精品巨大| 亚洲人成伊人成综合网2020| 国产视频一区二区在线看| 久久久国产精品麻豆| 精品久久久久久成人av| 久久欧美精品欧美久久欧美| 男人的好看免费观看在线视频 | 视频在线观看一区二区三区| 欧美日韩国产mv在线观看视频| 国产精品久久久人人做人人爽| 99香蕉大伊视频| 国产精品久久久久成人av| 国产蜜桃级精品一区二区三区| 国产av又大| 超碰成人久久| 久久久久久久午夜电影 | 19禁男女啪啪无遮挡网站| 国产亚洲欧美98| 在线观看免费视频日本深夜| 日日干狠狠操夜夜爽| 精品第一国产精品| 人人妻,人人澡人人爽秒播| 一级,二级,三级黄色视频| 免费一级毛片在线播放高清视频 | 啪啪无遮挡十八禁网站| 亚洲av片天天在线观看| 美女高潮喷水抽搐中文字幕| 中文字幕色久视频| 999久久久国产精品视频| 亚洲黑人精品在线| 午夜日韩欧美国产| 亚洲,欧美精品.| 操美女的视频在线观看| 欧美亚洲日本最大视频资源| 免费搜索国产男女视频| 亚洲色图综合在线观看| 另类亚洲欧美激情| 97超级碰碰碰精品色视频在线观看| 一边摸一边做爽爽视频免费| 国产伦人伦偷精品视频| 亚洲在线自拍视频| 美女午夜性视频免费| 视频在线观看一区二区三区| 看免费av毛片| 亚洲成人精品中文字幕电影 | 在线观看免费日韩欧美大片| 午夜激情av网站| 久久狼人影院| 精品一区二区三区四区五区乱码| 成年版毛片免费区| 国产精品九九99| 黑人猛操日本美女一级片| 欧美国产精品va在线观看不卡| 十八禁人妻一区二区| 亚洲色图综合在线观看| 黄片小视频在线播放| 别揉我奶头~嗯~啊~动态视频| 久久精品国产综合久久久| 神马国产精品三级电影在线观看 | 啦啦啦在线免费观看视频4| 成人三级黄色视频| 老鸭窝网址在线观看| 精品国产国语对白av| 日韩精品青青久久久久久| 国产黄色免费在线视频| 精品国产乱子伦一区二区三区| 成人三级黄色视频| 操出白浆在线播放| 国产伦一二天堂av在线观看| 亚洲av成人av| 极品人妻少妇av视频| 午夜两性在线视频| 中文字幕另类日韩欧美亚洲嫩草| 免费av中文字幕在线| 欧美日本亚洲视频在线播放| 夜夜爽天天搞| 色婷婷久久久亚洲欧美| av视频免费观看在线观看| 中国美女看黄片| 久99久视频精品免费| 日本三级黄在线观看| 国产精品电影一区二区三区| 女性被躁到高潮视频| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲国产看品久久| 日韩有码中文字幕| 韩国av一区二区三区四区| 午夜福利,免费看| 国产单亲对白刺激| 天堂俺去俺来也www色官网| 一夜夜www| 91大片在线观看| 99re在线观看精品视频| 午夜精品在线福利| 天堂俺去俺来也www色官网| 亚洲av熟女| 91大片在线观看| 亚洲国产精品sss在线观看 | 国产伦一二天堂av在线观看| 久久人人精品亚洲av| 中亚洲国语对白在线视频| 最新美女视频免费是黄的| 精品国产亚洲在线| videosex国产| 看片在线看免费视频| 成人三级黄色视频| 女人被狂操c到高潮| 免费在线观看视频国产中文字幕亚洲| 琪琪午夜伦伦电影理论片6080| 高清黄色对白视频在线免费看| 久久热在线av| 婷婷六月久久综合丁香| 黄色片一级片一级黄色片| 久久香蕉精品热| 久久久国产精品麻豆| 夜夜躁狠狠躁天天躁| 亚洲av成人不卡在线观看播放网| 在线观看日韩欧美| 国产成人精品久久二区二区免费| 欧美黑人欧美精品刺激| 法律面前人人平等表现在哪些方面| 国产成人欧美| 男人的好看免费观看在线视频 | 嫩草影院精品99| 国产在线精品亚洲第一网站| 涩涩av久久男人的天堂| 视频区图区小说| 黄色a级毛片大全视频| 国产男靠女视频免费网站| 精品人妻在线不人妻| 69精品国产乱码久久久| 亚洲精品在线美女| 国产一区二区在线av高清观看| 男女高潮啪啪啪动态图| 国产成人欧美| a级毛片在线看网站| 国产免费男女视频| 国产精品九九99| 午夜免费鲁丝| 黄色女人牲交| 男女高潮啪啪啪动态图| 日韩精品免费视频一区二区三区| 老司机午夜十八禁免费视频| 欧美+亚洲+日韩+国产| 18美女黄网站色大片免费观看| 无限看片的www在线观看| 巨乳人妻的诱惑在线观看| 性少妇av在线| 波多野结衣高清无吗| 嫁个100分男人电影在线观看| 国产精品1区2区在线观看.| 视频区欧美日本亚洲| 日韩大尺度精品在线看网址 | 别揉我奶头~嗯~啊~动态视频| 757午夜福利合集在线观看| 国产精品二区激情视频| 久久久久国产一级毛片高清牌| 免费在线观看影片大全网站| 在线观看日韩欧美| 美女国产高潮福利片在线看| 久久久久久久精品吃奶| 色老头精品视频在线观看| 99在线人妻在线中文字幕| 久久热在线av| 在线视频色国产色| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲成人精品中文字幕电影 | 在线观看免费高清a一片| 国产一区在线观看成人免费| 亚洲男人的天堂狠狠| 亚洲国产看品久久| 国产精品美女特级片免费视频播放器 | 热99re8久久精品国产| 超碰97精品在线观看| 免费高清在线观看日韩| 午夜福利影视在线免费观看| 国产91精品成人一区二区三区| 久久久国产一区二区| 精品一区二区三卡| 一进一出抽搐动态| 嫩草影视91久久| 电影成人av| a级毛片黄视频| cao死你这个sao货| netflix在线观看网站| 免费高清在线观看日韩| 国产精品 欧美亚洲| 免费在线观看影片大全网站| 国产精品偷伦视频观看了| 咕卡用的链子| 99热只有精品国产| 免费女性裸体啪啪无遮挡网站| 精品国内亚洲2022精品成人| 国产一区二区激情短视频| 国产1区2区3区精品| 成人免费观看视频高清| 黄色毛片三级朝国网站| 国产精品免费一区二区三区在线| 精品国产乱子伦一区二区三区| www.999成人在线观看| 一进一出抽搐动态| 免费看a级黄色片| 日韩一卡2卡3卡4卡2021年| 韩国精品一区二区三区| 欧美色视频一区免费| 一级a爱视频在线免费观看| 一级片免费观看大全| 久久国产精品影院| 少妇裸体淫交视频免费看高清 | 中文字幕最新亚洲高清| 91麻豆精品激情在线观看国产 | 日本黄色视频三级网站网址| 18禁裸乳无遮挡免费网站照片 | 亚洲精品国产精品久久久不卡| 9色porny在线观看| 久久久久久免费高清国产稀缺| tocl精华| 精品电影一区二区在线| 欧美乱妇无乱码| 国产精品免费视频内射| 国产成年人精品一区二区 | 99re在线观看精品视频| 精品熟女少妇八av免费久了| 婷婷精品国产亚洲av在线| 久久天躁狠狠躁夜夜2o2o| 国产成+人综合+亚洲专区| 好男人电影高清在线观看| 亚洲性夜色夜夜综合| 狂野欧美激情性xxxx| 精品少妇一区二区三区视频日本电影| 国产又色又爽无遮挡免费看| 十八禁人妻一区二区| 50天的宝宝边吃奶边哭怎么回事| 悠悠久久av| 国产一区二区三区在线臀色熟女 | 国产成+人综合+亚洲专区| 天堂动漫精品| 欧美日韩视频精品一区| 欧美日韩av久久| 宅男免费午夜| 在线av久久热| 9热在线视频观看99| 久久久久久久午夜电影 | 国产精品一区二区精品视频观看| 久久中文看片网| 精品日产1卡2卡| 亚洲精品在线美女| 成人黄色视频免费在线看| 在线观看免费视频日本深夜| 中文字幕人妻熟女乱码| 国产男靠女视频免费网站| 久久香蕉激情| 最好的美女福利视频网| av福利片在线| 波多野结衣av一区二区av| 精品久久久久久成人av| 男人舔女人下体高潮全视频| av天堂久久9| 99在线人妻在线中文字幕| 国产在线精品亚洲第一网站| 国产成人精品在线电影| 黑人欧美特级aaaaaa片| 久久中文看片网| 丁香六月欧美| 午夜免费成人在线视频| 欧美中文综合在线视频| 亚洲精品美女久久久久99蜜臀| 久久国产精品人妻蜜桃| 国产三级黄色录像| 久久精品亚洲精品国产色婷小说| 欧美日韩瑟瑟在线播放| 国产一区二区三区在线臀色熟女 | 一a级毛片在线观看| 男女高潮啪啪啪动态图| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品亚洲av一区麻豆| 欧美黑人精品巨大| 可以免费在线观看a视频的电影网站| 午夜福利欧美成人| av天堂久久9| 在线永久观看黄色视频| 一级毛片女人18水好多| 黄色女人牲交| 亚洲精品久久午夜乱码| 香蕉丝袜av| 9191精品国产免费久久| 大型av网站在线播放| 一区二区三区激情视频| 中文欧美无线码| 国产精品影院久久| 黄色片一级片一级黄色片| 丁香欧美五月| 亚洲精品国产色婷婷电影| 日韩精品免费视频一区二区三区| 高潮久久久久久久久久久不卡| 长腿黑丝高跟| 亚洲精品av麻豆狂野| av在线天堂中文字幕 | 国产一区在线观看成人免费| 99久久综合精品五月天人人| 亚洲成人国产一区在线观看| 亚洲国产欧美网| 精品国产一区二区三区四区第35| 淫妇啪啪啪对白视频| 国产精品自产拍在线观看55亚洲| 女人爽到高潮嗷嗷叫在线视频| 亚洲aⅴ乱码一区二区在线播放 | 欧美不卡视频在线免费观看 | 久久人妻福利社区极品人妻图片| 搡老乐熟女国产| 精品一区二区三卡| 国产极品粉嫩免费观看在线| 悠悠久久av| 人成视频在线观看免费观看| 久久精品亚洲熟妇少妇任你| 精品乱码久久久久久99久播| 免费久久久久久久精品成人欧美视频| 老司机靠b影院| 中文欧美无线码| 97人妻天天添夜夜摸| 精品免费久久久久久久清纯| 免费观看精品视频网站| 欧美老熟妇乱子伦牲交| 国产国语露脸激情在线看| 在线十欧美十亚洲十日本专区| 成在线人永久免费视频| 欧美日韩一级在线毛片| 亚洲精品在线美女| 亚洲国产精品一区二区三区在线| 亚洲午夜理论影院| 亚洲av电影在线进入| a级毛片黄视频| 一区在线观看完整版| 三上悠亚av全集在线观看| 国产单亲对白刺激| 亚洲,欧美精品.| 看免费av毛片| 国产成人精品久久二区二区91| 亚洲av片天天在线观看| 欧美激情高清一区二区三区| 黑人巨大精品欧美一区二区mp4| 久久久久九九精品影院| 欧美日本亚洲视频在线播放| 成人影院久久| 国产一卡二卡三卡精品| av电影中文网址| 久久久久国内视频| 在线观看66精品国产| 久久精品成人免费网站| 亚洲欧洲精品一区二区精品久久久| 他把我摸到了高潮在线观看| 国产色视频综合| 久久久久久免费高清国产稀缺| 国内毛片毛片毛片毛片毛片| 99久久99久久久精品蜜桃| 国产高清激情床上av| 中文字幕色久视频| 丝袜在线中文字幕| 99精品久久久久人妻精品| 丝袜在线中文字幕| 这个男人来自地球电影免费观看| 国产精品免费视频内射| 国产三级黄色录像| 国产1区2区3区精品| 在线视频色国产色| 国产区一区二久久| 欧美日韩黄片免| 久久午夜亚洲精品久久| 动漫黄色视频在线观看| 母亲3免费完整高清在线观看| 一边摸一边抽搐一进一出视频| 国产一区二区三区在线臀色熟女 | 999久久久国产精品视频| 熟女少妇亚洲综合色aaa.| 免费看a级黄色片| 精品第一国产精品| 亚洲欧美一区二区三区黑人| 又黄又爽又免费观看的视频| 女警被强在线播放| 一区二区日韩欧美中文字幕| 亚洲avbb在线观看| 五月开心婷婷网| 少妇 在线观看| tocl精华| 一区福利在线观看| 99久久99久久久精品蜜桃| 久久中文字幕一级| 国产成人欧美| √禁漫天堂资源中文www| 亚洲国产看品久久| 欧美大码av| 十八禁网站免费在线| 国内久久婷婷六月综合欲色啪| 久热这里只有精品99| 91在线观看av| 一区二区三区精品91| 97碰自拍视频| xxx96com| 亚洲国产精品一区二区三区在线| 黄网站色视频无遮挡免费观看| 大香蕉久久成人网| 国产成人欧美| 91麻豆av在线| 亚洲成人国产一区在线观看| 可以在线观看毛片的网站| 琪琪午夜伦伦电影理论片6080| 精品少妇一区二区三区视频日本电影| 成人国语在线视频| 亚洲国产精品999在线| 啦啦啦在线免费观看视频4| 久久天躁狠狠躁夜夜2o2o| 丝袜在线中文字幕| 黄色片一级片一级黄色片| 亚洲精品国产区一区二| 老鸭窝网址在线观看| 亚洲 欧美 日韩 在线 免费| av福利片在线| 欧美在线黄色| 99久久99久久久精品蜜桃| 欧美成狂野欧美在线观看| 国产三级在线视频| 99久久综合精品五月天人人| 久久久久久大精品| 麻豆一二三区av精品| 久久精品成人免费网站| 丁香欧美五月| 国产日韩一区二区三区精品不卡| 国产伦人伦偷精品视频| 天天躁夜夜躁狠狠躁躁| 亚洲欧美精品综合久久99| 国产亚洲欧美在线一区二区| 麻豆成人av在线观看| 欧美成狂野欧美在线观看| 婷婷六月久久综合丁香| 高清欧美精品videossex| 啪啪无遮挡十八禁网站| 国产精品综合久久久久久久免费 |