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

    基于 WQImin 模型的貴州草海水質(zhì)評(píng)價(jià)分析

    2025-08-04 00:00:00萬穎杜明普郭艷敏劉坤云孫赫英王英才李瑞雯
    關(guān)鍵詞:草海入湖水質(zhì)評(píng)價(jià)

    近年來,隨著社會(huì)經(jīng)濟(jì)發(fā)展和人口的增長,大量工業(yè)、農(nóng)業(yè)廢水和生活污水進(jìn)入湖泊,導(dǎo)致湖泊水質(zhì)惡化,對(duì)人類健康、生物多樣性產(chǎn)生嚴(yán)重威脅[1-2].水質(zhì)評(píng)價(jià)是水環(huán)境治理的重要基礎(chǔ)之一,旨在通過監(jiān)測(cè)水體的物理、化學(xué)和生物指標(biāo),采用科學(xué)方法對(duì)水體污染程度進(jìn)行綜合評(píng)估.在此基礎(chǔ)上,深入分析污染時(shí)空變化特征,識(shí)別關(guān)鍵污染因子,并采取針對(duì)性的水污染防治措施,具有重要意義[3].

    當(dāng)前,水質(zhì)評(píng)價(jià)方法主要有單因子評(píng)價(jià)法、主成分分析法、水質(zhì)綜合評(píng)價(jià)指數(shù)法等[4-6].其中,水質(zhì)綜合評(píng)價(jià)指數(shù)(WQI)法是將各種水質(zhì)參數(shù)整合,得到一個(gè)可代表整體水質(zhì)的無量綱數(shù)值,具有充分利用水質(zhì)參數(shù)信息、綜合反映水質(zhì)狀況等優(yōu)點(diǎn),在水資源管理中發(fā)揮著重要的作用 [7-8] .然而, WQI 法需要考慮大量水質(zhì)參數(shù),數(shù)據(jù)處理過程較為復(fù)雜.最簡水質(zhì)綜合評(píng)價(jià)指數(shù) )在 WQI 的基礎(chǔ)上,選取有限的關(guān)鍵指標(biāo)對(duì)水質(zhì)進(jìn)行評(píng)價(jià),能夠使分析和監(jiān)測(cè)成本大幅降低[5.9].研究表明[9-10], WQI 與 WQImin 具有較強(qiáng)的線性相關(guān)性,因此可通過 快速得到WQI.值得注意的是, WQImin 的關(guān)鍵指標(biāo)會(huì)隨著不同區(qū)域的水體特征而發(fā)生變化.由于所選定關(guān)鍵指標(biāo)和權(quán)重的差異,模型的精確度也會(huì)存在一定的差別[10-11].因此, 模型的性能需要根據(jù)特定的研究對(duì)象進(jìn)行精準(zhǔn)評(píng)估,

    草海作為貴州最大的高原天然淡水湖泊,對(duì)當(dāng)?shù)刈匀缓蜕鐣?huì)環(huán)境有著重要影響.當(dāng)前,草海面臨水質(zhì)富營養(yǎng)化和水生植物衰退等問題.對(duì)草海水質(zhì)時(shí)空分布特征進(jìn)行分析評(píng)價(jià),制定針對(duì)性的水污染防治措施十分必要.采用基于多項(xiàng)水質(zhì)指標(biāo)的WQI法對(duì)草海水質(zhì)進(jìn)行評(píng)價(jià),將有助于更加全面地理解草海水質(zhì)變化規(guī)律[12].而WQI與 模型用于高原湖泊的水質(zhì)評(píng)價(jià)還未見報(bào)道.因此,本研究采用單因子評(píng)價(jià)法探討了草海水質(zhì)基本特性,運(yùn)用WQI法揭示了草海及入湖支流水質(zhì)時(shí)空分布規(guī)律.在此基礎(chǔ)上,結(jié)合草海流域特點(diǎn)、氣候條件構(gòu)建僅考慮關(guān)鍵參數(shù)的 模型,以期降低分析測(cè)試成本.

    1 材料與方法

    1.1 研究區(qū)域概況

    草海 (26°47~26°52N,104°09~104°20E) 位于貴州省威寧彝族回族苗族自治縣,是一個(gè)由水域、沼澤、草甸、林地構(gòu)成的完整、典型的高原湖泊濕地生態(tài)系統(tǒng),為國家級(jí)自然保護(hù)區(qū).草海屬亞熱帶高原季風(fēng)氣候,流域年平均氣溫 10.9°C ,年日照時(shí)數(shù)介于 1 374~1 634h 之間,年均降雨量介于 626.5~1124.1mm 之間,豐水期為 6~10 月份,平水期為5、11月份,枯水期為 1~4 月份與12月份.草海流域面積 96km2 ,水域面積 25km2 ,最大水深小于 5m ,屬于淺水湖泊.草海表層沉積物氧化還原電位平均值為 (-102.1±4.34)mV 底質(zhì)環(huán)境受到歷史積累因素的影響,湖區(qū)重金屬風(fēng)險(xiǎn)指數(shù)較高.湖區(qū)水生植物以沉水植物和浮葉植物為主,敞水區(qū)物種相對(duì)較少.匯入草海的河流大多數(shù)是發(fā)源于泉水的短小河溪,主要有大中河、東山河、白馬河和萬下河.

    1.2 樣品采集與檢測(cè)

    研究區(qū)域共布設(shè)17個(gè)采樣點(diǎn),其中湖區(qū)13個(gè)采樣點(diǎn),4條入湖支流各1個(gè)采樣點(diǎn)(如圖1所示).采樣時(shí)間為2022年1月至2023年3月,按照《水質(zhì)采樣技術(shù)指導(dǎo)》(HJ494-2009)每月采樣1次,水質(zhì)指標(biāo)為水溫(t)、pH 值、溶解氧(DO)、電導(dǎo)率、濁度( Tur )、高錳酸鹽指數(shù)( CODMn )、化學(xué)需氧量( CODCr )、總氮(TN)、硝酸鹽氮( NO3 -N)、亞硝酸鹽氮( )、氨氮( NH4-N) !總磷(TP)、懸浮物(SS)、葉綠素 等14個(gè).

    圖1研究區(qū)域及水質(zhì)監(jiān)測(cè)點(diǎn)位置示意圖Fig.1Location of the study area and water quality monitoring sites

    pH 值、溶解氧、電導(dǎo)率 采用SX751型pH/ORP/電

    導(dǎo)率/溶解氧儀測(cè)定,濁度采用WGZ-20O 濁度計(jì)測(cè)定, CODMn 采用《水質(zhì)高酸鹽指數(shù)的測(cè)定》(GB/T11892一89)測(cè)定, CODCr 采用重鉻酸鹽法(HJ828-2017)測(cè)定,TN采用堿性過硫酸鉀消解紫外分光光度法(HJ636-2012)測(cè)定, NH4-N 采用納氏試劑分光光度法(HJ535-2009)測(cè)定,TP采用鉬酸銨分光光度法(GB11893一89)測(cè)定, NO3-N,NO2-N 采用CIC-D100 離子色譜儀測(cè)定(HJ84-2016),SS采用重量法(GB11901-89)測(cè)定, Chla 采用分光光度法(HJ897-2017)測(cè)定.

    1.3 WQI分析方法

    本研究采用WQI法對(duì)水質(zhì)進(jìn)行評(píng)價(jià),計(jì)算公式如下:

    式中,WQI為水質(zhì)綜合指數(shù); Ci 為第 i 種水質(zhì)指標(biāo)的標(biāo)準(zhǔn)化得分; Pi 為第 i 種水質(zhì)指標(biāo)的權(quán)重, Pi 的最小值為1,最大值為4.基于地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)( (GB3838-2002) 以及相關(guān)文獻(xiàn),確定各指標(biāo)的權(quán)重[10,13-14].隨后將各指標(biāo)進(jìn)行標(biāo)準(zhǔn)化處理, Ci 計(jì)算公式如下:

    式中, Ti 為第 i 個(gè)指標(biāo)的實(shí)際監(jiān)測(cè)值, Si,k 和 Si,k+n 為第 i 個(gè)指標(biāo)在標(biāo)準(zhǔn)化過程中所處的第 k 和 k+n 級(jí)對(duì)應(yīng)的標(biāo)準(zhǔn)化分?jǐn)?shù), Ii,k 為第 i 個(gè)指標(biāo)的歸一化值, n 為標(biāo)準(zhǔn)化值相同的個(gè)數(shù).值得注意的是,本研究中僅 CODCr 的第1類和第 I 類的標(biāo)準(zhǔn)化值相同, n=2 ,其他情況下 n=1 采用多元線性逐步回歸法確定 模型的關(guān)鍵參數(shù).基于 值將水質(zhì)劃分為5類,(80,100]為優(yōu),(60,80]為良,(40,60]為中,(20,40]為差,[0,20]為非常差.采用決定系數(shù) (R2) 評(píng)估 模型的擬合度,不同模型的預(yù)測(cè)能力由百分比誤差 (PE) 和均方根誤差 RMSE )判斷,二者數(shù)值越小,模型準(zhǔn)確度越高[15].采用2022年1至12月份監(jiān)測(cè)數(shù)據(jù)為 模型訓(xùn)練集數(shù)據(jù),2023年1至3月份為測(cè)試集數(shù)據(jù).

    表1WQI中各水質(zhì)因子分?jǐn)?shù)及權(quán)重Tab.1ScoreandweightofvariouswaterqualityfactorsinWQI
    注: pH 值在 6~9 之間時(shí),賦分為100分,否則為0.

    2 結(jié)果與分析

    2.1 草海水質(zhì)特性分析

    2022 年草海水質(zhì)變化特性如圖2所示.2022年草海最高水溫為 23.18°C ,最低水溫為 2.4°C.pH 平均值為8.37,表明草海水體整體處于堿性狀態(tài).年平均溶解氧質(zhì)量濃度為 6.82mg/L ,除2022年4月份與8月份外,整體優(yōu)于Ⅱ類標(biāo)準(zhǔn).水體電導(dǎo)率在 1~7 月份較高, 8~12 月份較低,約為 250μS/cm ,而濁度在2月份最低,隨后至8月份呈現(xiàn)出增加的趨勢(shì),二者可能與夏季降雨量較大有關(guān). CODMn 年均值為 8.46mg/L ,總體上高于Ⅲ類水體標(biāo)準(zhǔn).值得注意的是, CODCr 年均值為 31.41mg/L ,劣于地表水Ⅳ類水平.湖區(qū)TN與 NH4-N 值波動(dòng)較大,其中8月份 ΔNH4-N 達(dá)到最低值.而 Chla 與TN、 ΔNH4-N 的變化趨勢(shì)相反,7月份質(zhì)量濃度最高為 20.54μg/L. 研究表明,不同湖泊水質(zhì)特性在豐水期、枯水期存在明顯差異,部分湖泊水質(zhì)在時(shí)間上差異不明顯,而部分湖泊降雨徑流污染,導(dǎo)致豐水期水體中氮、磷濃度較高[16-17].草海 TP 年均質(zhì)量濃度為0.02mg/L ,最大值為 0.09mg/L ,優(yōu)于 I 類水.與貴州紅楓湖、云南滇池等高原湖泊類似,草海呈現(xiàn)出TP質(zhì)量濃度較低,TN質(zhì)量濃度較高的特性[18].將草海 CODCr、CODMn 、氮、磷等主要單因子指標(biāo)按照豐水期-平水期-枯水期進(jìn)行分析,結(jié)果如圖3所示.可以看到,豐水期草海 NH4–N 值較低,而TP在平水期值較高. CODCr 與 CODMn 在豐水期相對(duì)較高,但顯著性分析表明兩項(xiàng)指標(biāo)在豐水期和平水期沒有明顯差異( Plt;0.05 ,說明降雨徑流污染對(duì)湖泊 CODCr F CODMn 具有一定貢獻(xiàn),但影響程度有限,

    圖22022年草?;舅|(zhì)指標(biāo)變化特征Fig.2Variation ofwater quality indicators in Caohai in 2022
    圖32022年不同水期草海主要水質(zhì)指標(biāo)分析結(jié)果

    2.2 草海及入湖支流WQI分析評(píng)價(jià)

    圖4為湖區(qū)不同點(diǎn)位的WQI評(píng)價(jià)時(shí)間分布差異.從時(shí)間上來看,15月的平均值均為中等,所有月份的水質(zhì)指數(shù)均值均大于40.其中,2022年1月份和9月份水質(zhì)評(píng)價(jià)結(jié)果最佳,評(píng)分分別為56.82和55.67,5月份水質(zhì)評(píng)價(jià)結(jié)果最差,評(píng)分僅為49.83.從2022年1月份至7月份,草海WQI呈現(xiàn)出逐漸下降的趨勢(shì),隨后呈現(xiàn)出升高趨勢(shì),在9月份達(dá)到最佳水平,之后呈降低趨勢(shì).

    圖4草海湖區(qū)WQI時(shí)間分布差異Fig.4 Difference in temporal distribution of WQI in Caohai areε

    如圖5所示,草海4條人湖支流WQI指數(shù)均值均大于50.2022年1月份和3月份WQI指數(shù)表現(xiàn)最佳,均值分別為67.17、68.55.而4至7月份期間,WQI降低至較低水平.人湖河流的水質(zhì)在時(shí)間上與草海水質(zhì)呈現(xiàn)出相同的變化趨勢(shì),這可能是由于雨季集中在4至7月份,降雨徑流攜帶污染物進(jìn)入河道或湖泊,使水質(zhì)變差.草海人湖河流流量較小,水體中污染物濃度極易受到降雨的影響,導(dǎo)致WQI值隨時(shí)間波動(dòng)較大[19].

    草海及入湖支流的WQI空間分布如圖6所示.WQI指數(shù)在空間上的變化較小,不同點(diǎn)位的WQI指數(shù)均大于50,評(píng)價(jià)為中等.其中,點(diǎn)位C3和C4的WQI值較高,分別為53.36 和 53.29.C8~C13 點(diǎn)位WQI指數(shù)較低,說明草海西部水質(zhì)總體上優(yōu)于東部水質(zhì).對(duì)4條入湖河流而言,WQI均值由高到低:Z3(東山河)、Z4(萬下河)、Z2(白馬河)、Z1(大中河).其中東山河與萬下河水質(zhì)較為接近,WQI分別為63.29、63.26,水質(zhì)評(píng)價(jià)為良好;大中河水質(zhì)最差,WQI為56.99,水質(zhì)評(píng)價(jià)為中等.研究表明,湖泊水質(zhì)情況與湖泊周邊城市化率、人為活動(dòng)強(qiáng)度、土地利用類型等因素存在一定的相關(guān)性[19].草海流域東北部區(qū)域城市化率較高,污染負(fù)荷輸出顯著高于西南部區(qū)域,導(dǎo)致草海水質(zhì)污染程度由西向東漸增,城市生活污水的匯人可能是草海富營養(yǎng)化的主要原因[20].值得注意的是,4條入湖河流的WQI均值顯著高于湖區(qū),說明入湖河流水質(zhì)優(yōu)于湖區(qū).如表2所示,單因子評(píng)價(jià)得到草海湖區(qū)2022與2023年水質(zhì)類別均為V類,超標(biāo)因子為 CODCr ;4條入湖河流中,大中河和白馬河水質(zhì)較差,2022與2023年水質(zhì)類別均為V類,人湖支流水質(zhì)優(yōu)于湖區(qū)水質(zhì),該結(jié)果與WQI評(píng)價(jià)結(jié)果一致.由于草海濕地分布30 余條入湖溝渠,雨季時(shí)地表徑流沖刷的污染物可能通過溝渠進(jìn)入湖泊.此外,草海水生植物近年來大面積退化,其殘?bào)w沉積于湖底腐爛消解產(chǎn)生的內(nèi)源污染也可能成為重要的污染源[21].

    圖54條入湖河流WQI時(shí)間分布差異
    圖6草海及入湖支流WQI的空間分布差異Fig.6Differencesin spatial distribution of WQIin Caohai and rivers entering the lake

    2.3 基于 的水質(zhì)分析評(píng)價(jià)

    以14個(gè)水質(zhì)指標(biāo)為自變量,以WQI為因變量進(jìn)行多元線性逐步回歸分析,如表3所示.可以看到,CODMn 對(duì)WQI的貢獻(xiàn)最大,能夠解釋W(xué)QI變異的 46.4%(?lt;0.05) .將 TN、SS、 NH4. -N納入模型時(shí),模型的顯著性依次增加,分別能解釋W(xué)QI變異的 59.8%.68%?76.4% .而溶解氧、濁度、 Chla 作為模型的第5、第6、第7參數(shù),模型性能進(jìn)一步提高 (R2=0.873) :

    根據(jù)多元線性逐步回歸的結(jié)果,選擇 CODMn,TN,SS,NH4-N 作為 模型的基本參數(shù),溶解氧、濁度、 Chla 作為補(bǔ)充參數(shù).各 WQImin 模型包含的指標(biāo)如表4所示,其在訓(xùn)練集上的表現(xiàn)如圖7所示.通過比較RMSE,PE 和 R2 可以發(fā)現(xiàn),隨著參數(shù)的增多, 模型的準(zhǔn)確性逐步提高.其中,當(dāng) WQImin2 增加溶解氧指標(biāo)時(shí),模型擬合系數(shù)由0.775增加到了0.846,且 的決定系數(shù)均大于0.8,均表現(xiàn)出良好的預(yù)測(cè)性能.其中, 表現(xiàn)最佳,RMSE和 PE 值最低,分別為 1.72、2、44% ,決定系數(shù)達(dá)到0.895.此外,按照豐水期-平水期-枯水期分別構(gòu)建了不同水期的 WQImin 模型,如附錄表 S1~S2 所示.結(jié)果表明,豐水期 模型中關(guān)鍵參數(shù)為 CODMn 、 NO2-N 、濁度、SS、 NH4-N ,此時(shí)豐水期 模型 R2 為0.744;在平水期, 模型中關(guān)鍵參數(shù)為 CODMn 、溫度、 NH4-N 、濁度、SS,當(dāng)考慮以上參數(shù)時(shí),模型 R2 為0.900.枯水期 模型中關(guān)鍵參數(shù)為 CODMn 、 TN 、 Chla 、SS、 ΔNH4-N ,當(dāng)考慮以上參數(shù)時(shí),模型 R2 為0.882.

    Tab.2The average concentration of main water quality indicators in Caohai and rivers from 2022 to 2023
    表22022-2023年草海及入湖河流主要水質(zhì)指標(biāo)平均濃度表3多元線性逐步回歸分析
    Tab. 3 Multiplelinearstepwiseregressionanalysis表4不同 模型采用的指標(biāo)Tab.4Basic parameters of different models

    基于2023年1至3月份數(shù)據(jù)對(duì) WQImin 模型進(jìn)行了驗(yàn)證,如圖8所示.可以看到,4個(gè)模型在測(cè)試集上均表現(xiàn)出良好的擬合效果,其中 擬合效果最好,具有最高的決定系數(shù)( R2=0.788) , RMSE=2.68 ,PE=4.66% .值得注意的是, 模型在訓(xùn)練集上決定系數(shù) R2=0.846 ,當(dāng)繼續(xù)增加水質(zhì)指標(biāo)時(shí),相關(guān)性盡管有一定程度的增加,但增加幅度較小;而在測(cè)試集上,增加水質(zhì)指標(biāo), WQImin3 與 WQImin4 的準(zhǔn)確性略微下降.因此,綜合考慮模型精度和監(jiān)測(cè)成本后,可將 作為草海最佳的水質(zhì)預(yù)測(cè)模型.

    Fig.7Results of fitting the WQImin model to the WQI model based on the training set
    圖8基于測(cè)試集的 WQImin 模型與WQI模型的擬合結(jié)果Fig.8 Fitting results of WQIain model to WQI model based on test set

    3討論

    研究表明,評(píng)價(jià)對(duì)象水質(zhì)特征直接決定 WQImin 模型的關(guān)鍵因子[14].陳麗華等[12]采用TLI和PCA法對(duì)草海水質(zhì)進(jìn)行分析評(píng)價(jià)時(shí)發(fā)現(xiàn), CODMn 和 TN可能是影響草海水質(zhì)的主要原因.2013至2022年,草海CODMn,NH4-N,TN 顯著升高,水質(zhì)評(píng)價(jià)結(jié)果主要為V類,超標(biāo)因子為 CODMn 、 TN[22] .在本研究中, CODMn 、TN 在草海 模型中具有最高的解釋度.其中, CODMn 為首要參數(shù),是因?yàn)椴莺Kw中 CODMn 較高.草海沉積物中TN含量較高,甚至高于部分典型的富營養(yǎng)化湖泊[18].因此,TN作為反映草海水質(zhì)狀況的重要指標(biāo),有必要將其納人 模型[23].在WQI模型構(gòu)建中,選取的SS權(quán)重 Pi 要小于 ΔNH4-N 和溶解氧(表1),但在 模型中SS具有更高的解釋度,說明懸浮物顆粒在草海水質(zhì)評(píng)價(jià)中起到重要作用.溶解氧對(duì)水生生物生長等復(fù)雜的生化過程具有顯著影響,能夠反映出水體自凈能力,是水生態(tài)系統(tǒng)的敏感性指標(biāo),在草海 WQImin 模型得到了考慮[24].TP作為湖泊營養(yǎng)狀況評(píng)價(jià)的重要指標(biāo),并未納入草海 模型,是由于草海TP質(zhì)量濃度較低,對(duì)評(píng)價(jià)結(jié)果解釋力較弱[25].

    WQI法不僅能定性地評(píng)價(jià)水質(zhì)隨時(shí)間的變化,也能夠在空間上對(duì)不同水體斷面或人湖河流進(jìn)行單獨(dú)評(píng)價(jià)[26].WQI法對(duì)數(shù)據(jù)有較高的要求,需要盡可能包含全面的監(jiān)測(cè)指標(biāo),導(dǎo)致監(jiān)測(cè)成本增加. 通過篩選能夠準(zhǔn)確解釋水質(zhì)變化的關(guān)鍵指標(biāo),降低了水質(zhì)監(jiān)測(cè)的復(fù)雜性,在一定程度上彌補(bǔ)了WQI法的不足.本研究構(gòu)建的草海 WQImin2 模型中,所選取的5個(gè)指標(biāo)均易于監(jiān)測(cè),有利于當(dāng)?shù)厮|(zhì)分析評(píng)價(jià)工作的開展.值得注意的是,水質(zhì)評(píng)估中存在不確定性,主要受采樣時(shí)間和方式的影響[15].此外,標(biāo)準(zhǔn)化處理在WQI分析中同樣具有重要影響,部分指標(biāo)可以依據(jù)國家地表水水質(zhì)標(biāo)準(zhǔn) (GB3838-2002) 進(jìn)行評(píng)估,部分指標(biāo)還需要進(jìn)一步研究,以消除權(quán)重和歸一化參考范圍的影響[27].總體而言, 模型的構(gòu)建和驗(yàn)證過程具有很強(qiáng)的推廣性.在后續(xù)草海生態(tài)環(huán)境治理工作中,應(yīng)重點(diǎn)關(guān)注 CODMn 、 TN 、SS、 NH4 -N、DO等5項(xiàng)指標(biāo),同時(shí)可將濁度、 Chla 作為輔助參考指標(biāo).在污染源方面,除入湖河流外,環(huán)草海溝渠中輸入的污染物也需引起關(guān)注,并重視湖泊底泥產(chǎn)生的內(nèi)源污染.

    4結(jié)論

    本研究采用WQI法對(duì)草海水質(zhì)進(jìn)行綜合評(píng)價(jià),篩選得到關(guān)鍵指標(biāo)并構(gòu)建 模型,得到以下結(jié)論:(1)草海水體中TP質(zhì)量濃度較低,TN質(zhì)量濃度較高.此外,水體中 CODCr 質(zhì)量濃度長期高于地表水N類水平.(2)WQI評(píng)價(jià)表明草海和入湖河流水質(zhì)呈現(xiàn)出季節(jié)性變化,總體上評(píng)價(jià)為“中等\"水平,其中8至9月份水質(zhì)評(píng)價(jià)結(jié)果最優(yōu).而草海4條主要入湖支流的WQI評(píng)價(jià)結(jié)果優(yōu)于草海湖區(qū),表明入湖河流對(duì)草海水質(zhì)影響較小.(3)研究提出草海 模型,關(guān)鍵水質(zhì)指標(biāo)包括 CODMn 、TN、SS、 NH4 -N、溶解氧等 模型在水質(zhì)評(píng)價(jià)方面表現(xiàn)出優(yōu)異的性能.

    附錄見電子版(DOI:10.16366/j.cnki.1000-2367.2024.08.30.0001).

    參考文獻(xiàn)

    [1]GIRIS.Waterqualityprospetiveintwentyfirstcenturystatusofwaterqualityinmajorriverbasins,ontemporarystrategiesadediments:a review[J].Environmental Pollution,2021,271:116332.

    [2]宋高飛,朱宇軒,米武娟,等.武漢市湖泊浮游動(dòng)物群落特征及其影響因素[J].河南師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2022,50(3):135-142. SONGGF,ZHUYX,MWJ,etalCharacteristicsofzooplanktoncommunityandtheirinfluencingfactorsinlakesofWuhanJoual ofHenan Normal University (Natural Science Edition),2022,50(3):135-142.

    [3]HUANGYK,WANG XY,XIANG WJ,etal.Forward-lokingroadmapsforlong-termcontinuous waterqualitymonitoring:otle necks,innovationsndprospects incriticalrevieJEnvironmentalSienceamp;Technology,O22,56(9)534-354.

    [4]YANT,HENSLZHOUANIdicesand modelsofsurface waterqualityaessment:Reviewand perspectivesJ]EnvirotalPollution,2022,308:119611.

    [5] WUT,WANG SR,SUBL,etalUnderstandingthe waterqualitychangeof theYilongLakebasedoncomprehensive asessmentmethOceanology and Limnology,2020,42(2) :166-174.

    [7]UDDIN MG,NASH S,RAHMANA,etal.Acomprehensive methodfor improvementof waterquality index(WQI) modelsforcoastal water quality assessment[J].Water Research,2022,219:118532.

    [8]WUZS,LAIXJ,IKYWaterqualityassessentofriversinLkeChaouBasin(Chia)using waterqualityidexJEologicalIndi cators,2021,121:107021.

    [9]朱長軍,趙方星,李步東,等.基于主成分分析及 WQImin 的大黑汀水庫水質(zhì)評(píng)價(jià)[J].河南師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2021,49(3):52-58. ZHU CJ,ZHAOF X,LIBD,et al.Water quality evaluationof Daheiting reservoir based on principalcomponent analysis and WQImin [J]. Journal of Henan Normal University(Natural Science Edition),2021,49(3) :52-58.

    [10]程鵬,孫明東,郝韶楠.基于最簡水質(zhì)綜合評(píng)價(jià)指數(shù)的官廳水庫上游河流水質(zhì)評(píng)價(jià)[J].生態(tài)環(huán)境學(xué)報(bào),2023,32(2):372-380. CHENG P,SUNMD,HAOSN.WaterqualityassessmentofupstreamriversofGuantingReservoirbasedonthe simplest watequality index[J].Ecology and Environmental Sciences,2023,32(2) :372-380.

    [11]楊列,閆霄珂,劉艷麗,等.基于 ΔWQImin 模型的武漢市金銀湖水質(zhì)時(shí)空特性[J].長江科學(xué)院院報(bào),2022,39(4):49-55. YANG L,YANX K,LIUYL,etal.Analyzingspatio-temporalcharacteristicsof waterqualityinJinyinLakeusing waterqualitycrucial index( WQImin )[J].Journal of Yangtze River Scientific Research Institute,2022,39(4) :49-55.

    [12]陳麗華,喻記新,李麗,等.基于TLI和PCA的貴州草海水質(zhì)狀況評(píng)價(jià)[J].長江科學(xué)院院報(bào),2018,35(8):40-45. CHENLH,YUJX,LIL,etalWaterqualityassessmentofCaohai Lakein Guizhou basedontrophic state indexand principalcomo nent analysis[J].Journal of Yangtze River Scientific Research Institute,2O18,35(8) :40-45.

    [13]HOU WSUNSH,WANGMQ,etalAsessng waterqualityofivetypicalreservoirsilowerreachesofYellwRiver,ChnaUsinga water quality index method[J].Ecological Indicators,2o16,61:309-316.

    [14]唐琦,劉兵,王璞,等.改進(jìn)WQI在川中丘陵地區(qū)典型流域水質(zhì)評(píng)價(jià)中的應(yīng)用:以瓊江流域上游段為例[J].環(huán)境工程技術(shù)學(xué)報(bào),2022, 12(2):615-623. TANGQ,LIUG,etalApplicationofimprovedWQmodelinwaterqualitysessmntoftyalwatershedintehillf central Sichuan Province:a C se study in the upper reaches of Qiongjiang River Basin[J].Journal of Environmental Engineering Technolo gy,2022,12(2):615-623.

    [15]WU ZS,WANGXL,CHENYW,etal.Asessngriverwatequalityusing waterqualityindex inLakeTaiuBasi,ChinaJ]Sciencef The Total Environment,2018,6l2:914-922.

    [16]楊水化,彭正洪,焦洪贊,等.城市富營養(yǎng)化湖泊的外源污染負(fù)荷與貢獻(xiàn)解析;以武漢市后官湖為例[J].湖泊科學(xué),2020,32(4):941-951. YANG S H,PENG Z H,JIAOH Z,et alExternal pollution sourceload andcontributionofurban eutrophiclakes:TakingLake Houguan hu of Wuhan as an example[J].Journal of Lake Sciences,2020,32(4) :941-951.

    [17]劉發(fā)根,李梅,郭玉銀.鄱陽湖水質(zhì)時(shí)空變化及受水位影響的定量分析[J].水文,2014,34(4):37-43. LIUFG,LIM,GUOYY.Spatial-temporalvariationofwaterqualityandwaterlevelefectonwaterqualityinPoyangLakeJJounal of China Hydrology,2014,34(4) :37-43.

    [18]周茜,彭秋蕭,張才巧,等.貴州草海沉積物中養(yǎng)分的空間分布與評(píng)價(jià)[J].西南大學(xué)學(xué)報(bào)(自然科學(xué)版),2021,43(11):33-39. ZHOU X,PENGQX,ZHANGCQ,etal.patialdistributionand evaluationofnutrientsinthesedimentsofCahaiinGuzhouJJour nal of Southwest University(Natural Science Edition),2021,43(11) :33-39.

    [19]于超,儲(chǔ)金宇,白曉華,等.洱海入湖河流彌苴河下游氮磷季節(jié)性變化特征及主要影響因素[J].生態(tài)學(xué)報(bào),2011,31(23):119-126. YU C,CHUJY,BAIX H,etalSeasonal variationofntrogenand phosphorus in MijuRiverandLakeErhaiandinfluencingfactorsJ]. Acta Ecologica Sinica,2011,31(23):119-126.

    [20]周晨,喻理飛,蔡國俊,等.草海高原濕地湖泊水質(zhì)時(shí)空變化及水質(zhì)分區(qū)研究[J].水生態(tài)學(xué)雜志,2016,37(1):24-30. ZHOUC,YULF,CAIGJ,etalSpatial-temporalvariationofwaterqualityiCaohaiLakeJ]JurnalofHydroecology,6,7: 24-30.

    [21]ZHAOB,HU YS,YUHR,et alA methodforresearching the eutrophication andN/Ploads of plateau lakes:Lugu Lake asacase[J]. Science of The Total Environment,2023,876:162747.

    [22]徐林,郭艷敏,劉迪,等.2013—2022年貴州草海水質(zhì)變化與影響因素分析[J].環(huán)境保護(hù)科學(xué),2024,50(5):127-133. XUL,GUOYM,IUD,etal.TmporalandinfluencingfactorsofCaohai waterqualityfrom213to02J.EnvironmentalProtection Science,2024,50(5):127-133.

    [23]田瑜,趙海麗,常艷文,等.基于PCA- ?WQImin 模型的湯河水庫水質(zhì)評(píng)價(jià)[J].四川環(huán)境,2023,42(3):123-129. TIAN Y,ZHAO H L,CHANG Y W,et al.Water quality assessment of Tanghe Reservoir based on PCA and WQImin model[J].Sichuan Environment,2023,42(3):123-129.

    [24]QIJH,YNGLY,IUEF.Aolistcframeworkofwaterqualityevaluationusing waterqualityindex(WQIintheYieRiverChna) [J].Environmental Science and Pollution Research,2022,29(53):80937-80951.

    [25]ALITURANKOCERM,SEVGILIHParameters selectionforwaterqualityindexintheasessmntoftheenvironmentalimpactsof land-based trout farms[J].Ecological Indicators,2014,36:672-681.

    [26]VASISTHAPGANGULYWaterqualityaessmentofnaturallakesanditsimportanceanoverviewJ]aterialsTodroceed ings,2020,32:544-552.

    [27]NONG X Z,SHAODG,ZHONGH,etalEvaluationof waterquaityinthe South-to-North WaterDiversionProjetofChinausing the water quality index(WQI) method[J].Water Research,2020,178:115781.

    Evaluation of water quality of Guizhou Caohai based on the model

    Wan Yingl,Du Mingpul,Guo Yanmin',Liu Kunyun2,Sun Heying1,Wang Yingcai',Li Ruiwen'

    1.EcologyandEnvioentMonitoingndSitifiResearchCente,ChangngBasinEcologandEnvioentdmiistratioMistryf Ecologyand Environment,Wuhan430olo,China;2.Guizhou Ecological EnvironmentMonitoringCenter,Guiyang 550oo,China)

    Abstract:Caohai,thelargest plateau lake in Guizhou Province,has garnered significant attention due to recent water qualitydeterioration.This study assessed water qualityof Caohai using the water quality index(WQI)method,analyzing month-to-month monitoringdata spanning January 2022 to March 2O23from Caohai lakeand rivers entering the lake.Results indicatedalowtotalphosphorus(TP)contentbuthightotalnitrogen(TN)levels inthelake,withrainfallrunoffcontributing relatively minimally to nitrogen and phosphorus pollution.The average WQI value across the study area was 54.49 ,reflecting a \"medium\" waterquality evaluationgrade.Temporall,the best water qualityevaluation results were presentedfrom August to September.Spatiall,the western region exhibited superior waterqualitycompared totheeastern area,and inlet rivers showed better evaluations thanthelake,which exerted leser influence on lake water quality.Employing stepwise multiple linearregression,5 key indicators,including permanganate index( CODMn ),total nitrogen(TN),suspended solids(SS),ammonia nitrogen( NH4 -N),dissolved oxygen(DO),were selected from 14 water quality indicators,establishing the minimum comprehensive water quality index( )model. The model demonstrated robust predictive capabilities in both the training and test sets,showing an excellent fit (R2=0.846 ) with low root mean square error (RMSE=2.09 ) and a percentage error ratio of (20 3.45% in the training set.The resultsof thestudyare intended to providea decision-making basis for the ecological management of the Caohai.

    Keywords: plateau lake; Caohai; water quality assessment; minimum comprehensive water quality index( WQImin )spatial and temporal distribution characteristics

    表 S1豐水期多元線性逐步回歸分析
    表S2平水期多元線性逐步回歸分析
    表 S3枯水期多元線性逐步回歸分析

    猜你喜歡
    草海入湖水質(zhì)評(píng)價(jià)
    基于Python的因子分析在水質(zhì)自動(dòng)監(jiān)測(cè)數(shù)據(jù)評(píng)價(jià)中的應(yīng)用
    撫仙湖流域水環(huán)境問題與保護(hù)治理研究
    滇池流域水環(huán)境問題與保護(hù)治理對(duì)策研究
    異龍湖流域入河排污口排查與整治研究
    久久九九热精品免费| 亚洲人成网站高清观看| 免费看日本二区| 最新在线观看一区二区三区| 亚洲五月婷婷丁香| 黄片大片在线免费观看| 性色avwww在线观看| 99国产精品一区二区蜜桃av| 综合色av麻豆| 五月玫瑰六月丁香| 久久久成人免费电影| 男女下面进入的视频免费午夜| 长腿黑丝高跟| 国产精品 国内视频| 99久国产av精品| 熟妇人妻久久中文字幕3abv| 精品久久久久久,| 噜噜噜噜噜久久久久久91| 一本久久中文字幕| 国产中年淑女户外野战色| 一夜夜www| 欧美日韩乱码在线| 亚洲av电影在线进入| 最近视频中文字幕2019在线8| 国产精品一区二区三区四区久久| 男女下面进入的视频免费午夜| 日本五十路高清| 国产三级中文精品| 女人高潮潮喷娇喘18禁视频| 国产三级中文精品| 此物有八面人人有两片| 97碰自拍视频| 好男人电影高清在线观看| 日韩人妻高清精品专区| 波多野结衣高清无吗| 99精品欧美一区二区三区四区| 日韩欧美精品免费久久 | 在线播放国产精品三级| 婷婷精品国产亚洲av| 日本与韩国留学比较| 亚洲在线自拍视频| 国产aⅴ精品一区二区三区波| 国产视频内射| 亚洲精品影视一区二区三区av| 日本免费a在线| 免费大片18禁| 一夜夜www| 夜夜爽天天搞| 女警被强在线播放| 国产美女午夜福利| 国产日本99.免费观看| e午夜精品久久久久久久| 国产精品亚洲av一区麻豆| 欧美中文日本在线观看视频| 一卡2卡三卡四卡精品乱码亚洲| 久久国产精品人妻蜜桃| 成人国产综合亚洲| 国产精品亚洲一级av第二区| 欧洲精品卡2卡3卡4卡5卡区| 91在线精品国自产拍蜜月 | 日本一二三区视频观看| 免费在线观看亚洲国产| 精品人妻一区二区三区麻豆 | 中文字幕高清在线视频| 国产精品亚洲一级av第二区| 丁香欧美五月| 久久亚洲精品不卡| 淫秽高清视频在线观看| 亚洲av第一区精品v没综合| 亚洲最大成人手机在线| 免费在线观看成人毛片| 丰满人妻熟妇乱又伦精品不卡| 国产av一区在线观看免费| 一本精品99久久精品77| 日韩欧美在线二视频| 首页视频小说图片口味搜索| 男插女下体视频免费在线播放| 欧美日韩福利视频一区二区| 国产成人aa在线观看| 久久天躁狠狠躁夜夜2o2o| 日韩亚洲欧美综合| 精品欧美国产一区二区三| 久久久色成人| 久久久久性生活片| 一卡2卡三卡四卡精品乱码亚洲| 黄色丝袜av网址大全| 欧美3d第一页| 88av欧美| 亚洲18禁久久av| 国产精品日韩av在线免费观看| 美女大奶头视频| avwww免费| 午夜福利在线观看免费完整高清在 | 热99在线观看视频| 97超视频在线观看视频| 久久香蕉国产精品| 90打野战视频偷拍视频| 国产av不卡久久| 欧美在线一区亚洲| 美女cb高潮喷水在线观看| 成熟少妇高潮喷水视频| 又黄又粗又硬又大视频| 在线免费观看的www视频| 真人一进一出gif抽搐免费| 国产激情欧美一区二区| 在线视频色国产色| 高清毛片免费观看视频网站| 久久精品国产综合久久久| 最好的美女福利视频网| 成人无遮挡网站| 淫妇啪啪啪对白视频| 99久久久亚洲精品蜜臀av| 亚洲av日韩精品久久久久久密| 97超级碰碰碰精品色视频在线观看| 国产成+人综合+亚洲专区| 3wmmmm亚洲av在线观看| 午夜福利在线观看吧| 国产 一区 欧美 日韩| 人妻夜夜爽99麻豆av| 欧美日韩中文字幕国产精品一区二区三区| 亚洲成a人片在线一区二区| 老司机深夜福利视频在线观看| 午夜福利欧美成人| 99在线视频只有这里精品首页| 一级毛片高清免费大全| 搡女人真爽免费视频火全软件 | 熟妇人妻久久中文字幕3abv| 久久久久久久久久黄片| 久久久国产成人精品二区| 国产黄片美女视频| www.999成人在线观看| 国产成人aa在线观看| 国产成+人综合+亚洲专区| 最新中文字幕久久久久| 国产真实乱freesex| 亚洲国产日韩欧美精品在线观看 | 真实男女啪啪啪动态图| 国产野战对白在线观看| 男人的好看免费观看在线视频| 窝窝影院91人妻| 亚洲精品在线观看二区| 99国产精品一区二区蜜桃av| 精品国内亚洲2022精品成人| 久久人人精品亚洲av| 日韩av在线大香蕉| 看免费av毛片| 国产成人系列免费观看| 国产真实乱freesex| 搡老熟女国产l中国老女人| avwww免费| 欧美乱妇无乱码| 757午夜福利合集在线观看| 久久精品国产清高在天天线| 中文字幕精品亚洲无线码一区| 成人永久免费在线观看视频| 男女床上黄色一级片免费看| 少妇的丰满在线观看| 又紧又爽又黄一区二区| 成人特级av手机在线观看| 又爽又黄无遮挡网站| 国产精品亚洲一级av第二区| 久久久久久国产a免费观看| 久久草成人影院| 麻豆成人av在线观看| 亚洲熟妇熟女久久| 精品99又大又爽又粗少妇毛片 | 精品福利观看| 成人特级av手机在线观看| 亚洲av不卡在线观看| 日韩中文字幕欧美一区二区| 亚洲av电影在线进入| 欧美中文综合在线视频| 中文字幕av成人在线电影| 日本a在线网址| 亚洲第一电影网av| 国产午夜精品久久久久久一区二区三区 | 国产成人欧美在线观看| 欧美日韩黄片免| av在线天堂中文字幕| 日本在线视频免费播放| 亚洲美女黄片视频| 欧美色视频一区免费| 特级一级黄色大片| 国产精品三级大全| 亚洲不卡免费看| 一区二区三区高清视频在线| 老鸭窝网址在线观看| 窝窝影院91人妻| 国产真人三级小视频在线观看| 亚洲av免费高清在线观看| 国产成+人综合+亚洲专区| 欧美午夜高清在线| 女人被狂操c到高潮| 日本成人三级电影网站| 有码 亚洲区| 国产成人欧美在线观看| 日韩免费av在线播放| 69av精品久久久久久| 亚洲欧美日韩高清在线视频| 色在线成人网| av女优亚洲男人天堂| 国产免费男女视频| 一二三四社区在线视频社区8| 久久草成人影院| 日本黄大片高清| а√天堂www在线а√下载| 国产成人av激情在线播放| 亚洲美女黄片视频| bbb黄色大片| 一个人看的www免费观看视频| 国产高清视频在线观看网站| 内地一区二区视频在线| 国产真实伦视频高清在线观看 | 久久久久久久精品吃奶| 免费av观看视频| 久久精品国产清高在天天线| 又爽又黄无遮挡网站| 国产一级毛片七仙女欲春2| 日本a在线网址| 一级毛片女人18水好多| 亚洲美女视频黄频| 国产亚洲欧美98| 九九在线视频观看精品| 久久6这里有精品| 成人欧美大片| 欧美黄色片欧美黄色片| 亚洲久久久久久中文字幕| 日韩人妻高清精品专区| 午夜精品在线福利| 变态另类成人亚洲欧美熟女| 岛国在线免费视频观看| 亚洲国产日韩欧美精品在线观看 | 中文字幕av在线有码专区| 久久这里只有精品中国| 国产精品美女特级片免费视频播放器| 夜夜看夜夜爽夜夜摸| 色综合站精品国产| 日韩大尺度精品在线看网址| 日韩欧美一区二区三区在线观看| 成熟少妇高潮喷水视频| 日韩成人在线观看一区二区三区| 欧美xxxx黑人xx丫x性爽| 亚洲精品乱码久久久v下载方式 | 亚洲av免费在线观看| 欧美性猛交黑人性爽| 日韩有码中文字幕| 脱女人内裤的视频| 在线观看66精品国产| 男女做爰动态图高潮gif福利片| 欧美成人性av电影在线观看| 欧美色欧美亚洲另类二区| 欧美高清成人免费视频www| 国产精品综合久久久久久久免费| 日本精品一区二区三区蜜桃| 久久精品91无色码中文字幕| 精品久久久久久成人av| 欧美日韩瑟瑟在线播放| 无遮挡黄片免费观看| 国产精品日韩av在线免费观看| 美女cb高潮喷水在线观看| 亚洲久久久久久中文字幕| 欧美精品啪啪一区二区三区| 亚洲成人精品中文字幕电影| 成人特级av手机在线观看| 两人在一起打扑克的视频| 免费av不卡在线播放| 乱人视频在线观看| 色综合婷婷激情| 国产成年人精品一区二区| 亚洲久久久久久中文字幕| av中文乱码字幕在线| 青草久久国产| 日韩国内少妇激情av| 亚洲av日韩精品久久久久久密| 最近最新中文字幕大全电影3| 国产成人系列免费观看| 18禁黄网站禁片午夜丰满| 亚洲人成网站在线播放欧美日韩| 国产精品久久久久久人妻精品电影| 人人妻人人看人人澡| 男女床上黄色一级片免费看| svipshipincom国产片| 美女大奶头视频| 美女 人体艺术 gogo| 日韩中文字幕欧美一区二区| 亚洲av第一区精品v没综合| 国产激情偷乱视频一区二区| 黄色女人牲交| 特大巨黑吊av在线直播| 欧美av亚洲av综合av国产av| xxxwww97欧美| 欧美又色又爽又黄视频| 少妇人妻精品综合一区二区 | 欧美日韩瑟瑟在线播放| 国产精品爽爽va在线观看网站| 淫妇啪啪啪对白视频| 每晚都被弄得嗷嗷叫到高潮| 免费高清视频大片| 国产精品日韩av在线免费观看| 成人午夜高清在线视频| 成人18禁在线播放| 少妇熟女aⅴ在线视频| 国产成人系列免费观看| 欧美zozozo另类| 在线看三级毛片| 精品久久久久久久人妻蜜臀av| 一区二区三区高清视频在线| 宅男免费午夜| 精品熟女少妇八av免费久了| 国产精品一区二区免费欧美| 身体一侧抽搐| 午夜福利免费观看在线| e午夜精品久久久久久久| 在线十欧美十亚洲十日本专区| 久久久成人免费电影| 黄色片一级片一级黄色片| 九九在线视频观看精品| 午夜老司机福利剧场| 亚洲av五月六月丁香网| 岛国在线免费视频观看| 日韩精品中文字幕看吧| 亚洲av一区综合| 18+在线观看网站| 亚洲熟妇熟女久久| 日日摸夜夜添夜夜添小说| 国产高清视频在线观看网站| 国产色婷婷99| 99在线人妻在线中文字幕| 久久久久久久亚洲中文字幕 | 99国产精品一区二区蜜桃av| 国产精品国产高清国产av| 国内毛片毛片毛片毛片毛片| 搡老熟女国产l中国老女人| 成人特级av手机在线观看| 国产激情欧美一区二区| 国产一区二区亚洲精品在线观看| 精品久久久久久久人妻蜜臀av| 少妇的逼好多水| 黑人欧美特级aaaaaa片| 欧美黑人欧美精品刺激| 高潮久久久久久久久久久不卡| 久久中文看片网| 亚洲国产精品久久男人天堂| 婷婷六月久久综合丁香| 丰满人妻熟妇乱又伦精品不卡| 久久亚洲精品不卡| 久久久久精品国产欧美久久久| 热99在线观看视频| 国产在视频线在精品| 久久久成人免费电影| 久久伊人香网站| 国产精品香港三级国产av潘金莲| 日本在线视频免费播放| 免费电影在线观看免费观看| 黑人欧美特级aaaaaa片| 国产国拍精品亚洲av在线观看 | 婷婷六月久久综合丁香| 欧美3d第一页| 老熟妇仑乱视频hdxx| 国产黄色小视频在线观看| 欧美中文综合在线视频| 老汉色∧v一级毛片| 欧美日韩中文字幕国产精品一区二区三区| 夜夜夜夜夜久久久久| 搞女人的毛片| 男女做爰动态图高潮gif福利片| 亚洲欧美日韩高清专用| 麻豆一二三区av精品| 日韩 欧美 亚洲 中文字幕| 给我免费播放毛片高清在线观看| 国产伦一二天堂av在线观看| 久久精品夜夜夜夜夜久久蜜豆| 久久久色成人| 老司机深夜福利视频在线观看| 丰满人妻一区二区三区视频av | 国产精品99久久久久久久久| 亚洲国产欧洲综合997久久,| 又爽又黄无遮挡网站| 日本免费一区二区三区高清不卡| 久99久视频精品免费| 日本免费一区二区三区高清不卡| 欧美色欧美亚洲另类二区| 日本与韩国留学比较| 国产97色在线日韩免费| 97人妻精品一区二区三区麻豆| 欧美中文日本在线观看视频| 国产不卡一卡二| 亚洲男人的天堂狠狠| 男女做爰动态图高潮gif福利片| 丰满人妻熟妇乱又伦精品不卡| 精品国产美女av久久久久小说| 国产高清videossex| 高清日韩中文字幕在线| 亚洲在线自拍视频| 欧美一区二区国产精品久久精品| 亚洲成av人片在线播放无| 老汉色av国产亚洲站长工具| 色av中文字幕| 一级黄色大片毛片| 他把我摸到了高潮在线观看| 国产亚洲av嫩草精品影院| 国产激情欧美一区二区| 在线视频色国产色| 女警被强在线播放| 欧美不卡视频在线免费观看| 亚洲国产欧美人成| svipshipincom国产片| 亚洲真实伦在线观看| 中国美女看黄片| 无限看片的www在线观看| АⅤ资源中文在线天堂| 欧美国产日韩亚洲一区| 人人妻,人人澡人人爽秒播| 免费看十八禁软件| 免费看光身美女| 亚洲无线在线观看| 国产高清三级在线| 午夜福利欧美成人| aaaaa片日本免费| 日本黄大片高清| 欧美xxxx黑人xx丫x性爽| 久久精品国产亚洲av涩爱 | 叶爱在线成人免费视频播放| 99久久九九国产精品国产免费| 久久天躁狠狠躁夜夜2o2o| 亚洲国产高清在线一区二区三| 一本一本综合久久| 91字幕亚洲| 国产97色在线日韩免费| 搡女人真爽免费视频火全软件 | 久久这里只有精品中国| 精品一区二区三区av网在线观看| 成年女人永久免费观看视频| 窝窝影院91人妻| 欧美黄色片欧美黄色片| 日本撒尿小便嘘嘘汇集6| 两性午夜刺激爽爽歪歪视频在线观看| 国产免费av片在线观看野外av| 精品无人区乱码1区二区| 69av精品久久久久久| 小蜜桃在线观看免费完整版高清| 国产精品久久久久久精品电影| 欧美中文日本在线观看视频| 日韩人妻高清精品专区| 久久久色成人| av福利片在线观看| 成人午夜高清在线视频| 亚洲无线观看免费| 色综合站精品国产| 久久精品亚洲精品国产色婷小说| 色精品久久人妻99蜜桃| 亚洲成人久久性| 97超级碰碰碰精品色视频在线观看| 18禁黄网站禁片免费观看直播| 女同久久另类99精品国产91| 日韩亚洲欧美综合| 国产野战对白在线观看| 亚洲第一电影网av| 成年女人永久免费观看视频| 亚洲最大成人中文| 国产亚洲精品久久久com| 久久精品人妻少妇| 久久草成人影院| 欧美午夜高清在线| 小蜜桃在线观看免费完整版高清| а√天堂www在线а√下载| 精品不卡国产一区二区三区| 一本一本综合久久| 午夜影院日韩av| 日本三级黄在线观看| 波多野结衣高清作品| 亚洲国产日韩欧美精品在线观看 | 女同久久另类99精品国产91| 可以在线观看的亚洲视频| 色在线成人网| 亚洲avbb在线观看| 日本精品一区二区三区蜜桃| 91九色精品人成在线观看| 国产精品久久久久久人妻精品电影| 日韩中文字幕欧美一区二区| 国内精品久久久久久久电影| 国产伦精品一区二区三区视频9 | 亚洲,欧美精品.| 中文字幕高清在线视频| 国产精品久久久久久人妻精品电影| 国产99白浆流出| 国产视频内射| 香蕉av资源在线| av专区在线播放| 性色av乱码一区二区三区2| 小说图片视频综合网站| 少妇人妻精品综合一区二区 | 一二三四社区在线视频社区8| 在线观看舔阴道视频| 亚洲国产精品成人综合色| 精品一区二区三区人妻视频| 又爽又黄无遮挡网站| 宅男免费午夜| 在线观看舔阴道视频| 久久天躁狠狠躁夜夜2o2o| 一进一出抽搐gif免费好疼| 亚洲一区高清亚洲精品| 久久久久久久久大av| 亚洲av中文字字幕乱码综合| 亚洲aⅴ乱码一区二区在线播放| 国内精品一区二区在线观看| 最近视频中文字幕2019在线8| 啦啦啦韩国在线观看视频| 亚洲,欧美精品.| 亚洲天堂国产精品一区在线| 五月伊人婷婷丁香| 欧美3d第一页| 久久精品人妻少妇| 久久九九热精品免费| 国产老妇女一区| 免费搜索国产男女视频| h日本视频在线播放| 母亲3免费完整高清在线观看| 国产亚洲精品综合一区在线观看| 日韩欧美三级三区| 日本一本二区三区精品| 18禁国产床啪视频网站| 国产高清videossex| 久久久久久久久大av| 最新美女视频免费是黄的| 久久中文看片网| 国产精品一区二区免费欧美| 88av欧美| 啦啦啦观看免费观看视频高清| 老汉色av国产亚洲站长工具| 啦啦啦韩国在线观看视频| 亚洲精品粉嫩美女一区| 日日摸夜夜添夜夜添小说| 久久精品综合一区二区三区| 丁香六月欧美| 露出奶头的视频| 韩国av一区二区三区四区| 露出奶头的视频| 老司机在亚洲福利影院| 三级毛片av免费| 丰满的人妻完整版| 91九色精品人成在线观看| 一本久久中文字幕| 在线视频色国产色| 脱女人内裤的视频| 一个人观看的视频www高清免费观看| 99久久久亚洲精品蜜臀av| 精品99又大又爽又粗少妇毛片 | 国产精品免费一区二区三区在线| 国产伦在线观看视频一区| 毛片女人毛片| 成人三级黄色视频| 国产精品电影一区二区三区| 夜夜爽天天搞| 人人妻,人人澡人人爽秒播| 国产综合懂色| 男女之事视频高清在线观看| 每晚都被弄得嗷嗷叫到高潮| 国产成人啪精品午夜网站| 淫秽高清视频在线观看| 久久香蕉精品热| 99精品久久久久人妻精品| 国产精品影院久久| 久久久久性生活片| 亚洲乱码一区二区免费版| 99久国产av精品| 香蕉av资源在线| 国产精品一区二区免费欧美| 国产高清视频在线观看网站| 啦啦啦韩国在线观看视频| 午夜福利免费观看在线| 欧美乱色亚洲激情| 夜夜爽天天搞| 一区二区三区激情视频| 精品欧美国产一区二区三| 91字幕亚洲| 亚洲精品456在线播放app | 在线观看舔阴道视频| 成人午夜高清在线视频| 在线观看免费视频日本深夜| 麻豆国产av国片精品| 午夜免费成人在线视频| 亚洲欧美精品综合久久99| 一进一出好大好爽视频| 99精品在免费线老司机午夜| av在线天堂中文字幕| 日韩精品中文字幕看吧| 美女cb高潮喷水在线观看| 日韩人妻高清精品专区| 狂野欧美激情性xxxx| 岛国视频午夜一区免费看| 午夜日韩欧美国产| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产高清在线一区二区三| 熟女人妻精品中文字幕| 成人三级黄色视频| 级片在线观看| 日韩精品中文字幕看吧| 久久亚洲真实| 精品电影一区二区在线| 亚洲国产中文字幕在线视频| 国产成人福利小说| 成人高潮视频无遮挡免费网站| 青草久久国产| 亚洲自拍偷在线| 成熟少妇高潮喷水视频| 哪里可以看免费的av片| 搞女人的毛片| 毛片女人毛片| 亚洲精品日韩av片在线观看 | 日本熟妇午夜| 禁无遮挡网站| 亚洲成人精品中文字幕电影| 国产av一区在线观看免费| 国产av不卡久久|