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

    南湖流域典型支流面源污染負(fù)荷時(shí)空變化分析與模擬

    2025-08-03 00:00:00劉康于鵬躍孫新平鄒高峰伍藝鄒韜黃穎費(fèi)講馳
    湖南農(nóng)業(yè)科學(xué) 2025年6期
    關(guān)鍵詞:磷素南湖監(jiān)測點(diǎn)

    中圖分類號:X52 文獻(xiàn)標(biāo)識碼:A 文章編號:1006-060X(2025)06-0052-07

    引用格式:,等.南湖流域典型支流面源污染負(fù)荷時(shí)空變化分析與模擬[J].湖南農(nóng)業(yè)科學(xué),2025(6):52-58

    DOI:10.16498/j.cnki.hnnykx.2025.006.009

    Analysis and Simulation of Spatio-Temporal Variations of Non-Point Source Pollution Loads in Typical Tributaries of the Nanhu Lake Basin

    LIUKang1,YU Peng-yue1,SUNXin-ping2,ZOUGao-feng3,WU Yi4,ZOUTao5, HUANG Ying1,F(xiàn)EI Jiang-chi1

    (1.HunanAcuualUsitC;urallteterfCt PRC; 3.AgricultureandRural Affirs BureauofNingxiang CityChangsha 410699,PRC; 4.Meitanba TounAgricultural Comprehensive Service Center of Ningxiang City, Changsha 410609, PRC; 5. Hunan Tobacco Company Hengyang Company Changning Branch, Changning 42150o, PRC)

    Abstract:Thisstudytookthebasinof WangjiaRiver,atypicaltributaryof teNanhuLakebasininYueyangCity,HunanProvince, as theresearcharea.Onthebasisofconducting the monitoringof phosphorus polutants inthe water environment,andfurther consideringtheinfenesofthehydrologicalenvironment,topographyandlandforms,ndlandusetypesofthebasi,tisstudy constructedandoptimizedamulti-salephosphorus source-sinktransportmodel integratingenvironmental informationandspatial remotesensingthrough MIKE21 software,andanalyzedandsimulatedthe spatio-temporalvariationsof phosphorus non-point sourcepolltionloads,soastoprovidereferencefortheprecisepreventionandcontrolofphosporuspollutioninwaterbodie.he monitoringresultsshowedthatthepeakperiodofphosphorusconcentrationfromvarious polutionsources in WangjiaRiver from September2O23toAugust2O24wasconcentratedinwinterandspring.Theaveragetotalphosphorusconcentrationinthewaterbody ofWangjiaRiverwas 0.78mg/L ,and the average total phosphorus concentration in the water body of the catchment area was 0.61 mg/L .Both exceeded the Class ΔV water quality concentration standard in the GB 3838—2O02 Surface Water Environmental Quality Standards.Intemsoftemporalvariation,theaveragetotalphosphorusconcentrationinthewaterbodiesofWangiaRiverandthe catchment area bothshowedthehighest levelsinspring andthelowest levelsinautumn.Intermsofspatial variation,theaverage totalphosphorusconcentrationinthewaterbodyofWangjiaRiverfromtheupperreaches tothelowerreachesshowedatrendoffirst increasingand thendecreasing.The average totalphosphorusconcentrationofthe waterbodyin the middleofthecatchmentarea was higherthanthatonboth sidesinspring,while theopposite trendoccurred inautumn.The average acuracyoftheoptimized model simulation increasedtoabout 70% ,andthe trend of thesimulation results was basically consistent with that of the measured results.

    Keywords: non-point source pollution; total phosphorus in the water body;MIKE 21;spatio-temporal variations;Nanhu Lake; Wangjia River

    位于湖南省岳陽市的南湖是洞庭湖的主要湖汊之一,其歷來是周邊漁民們的生存之本,也成為周邊大部分生活、生產(chǎn)活動的排污地[1-2]。近年來,在周邊區(qū)域工業(yè)化、城鎮(zhèn)化進(jìn)程加快的背景下,由于環(huán)保規(guī)劃、管理不夠到位等,南湖流域出現(xiàn)面源污染加重、生態(tài)功能退化等一系列問題[3]。岳陽市環(huán)境監(jiān)測中心數(shù)據(jù)表明,南湖水體中磷污染程度最高[46]。磷污染問題是我國水體污染防治領(lǐng)域的重要關(guān)注點(diǎn)之一[7-8],磷污染不僅會加速水體富營養(yǎng)化進(jìn)程,減少生物多樣性,也會影響到工農(nóng)業(yè)生產(chǎn)和居民生活。然而,目前磷污染物在南湖流域的遷移動態(tài)尚未得到很好的評估,因此,開展南湖流域磷污染程度趨勢變化的評價(jià)研究具有迫切的現(xiàn)實(shí)意義[69]。

    當(dāng)前主流的2種評價(jià)流域污染程度趨勢變化的方法是進(jìn)行實(shí)地勘探和構(gòu)建水質(zhì)模型?;趯?shí)地勘探的評價(jià)方法通常需要大量的觀測數(shù)據(jù),而觀測對時(shí)間和資源的要求較高。因此,通過構(gòu)建水質(zhì)模型對水體污染負(fù)荷進(jìn)行模擬分析成為國內(nèi)外學(xué)者的研究重點(diǎn)。如Wang等[1o使用水質(zhì)分析模擬程序(WaterQualityAnalysisSimulationProgram,WASP)研究污染物排放總量和地理分布特征;Venishetty等[11]使用土壤與水評估工具(SoilandWaterAssessmentTool,SWAT)對流域內(nèi)污染物負(fù)荷進(jìn)行評價(jià)分析。這幾種模型主要針對長時(shí)段的流域綜合模擬,但針對小流域以及排污口分布模擬的模塊應(yīng)用則較為單一。相比于上述幾種模型,MIKE模型在污染物傳輸過程模擬以及污染物空間分布輸出方面具有一定優(yōu)勢。如Lu等[12]利用MIKE耦合模型模擬流域內(nèi)沉積物和農(nóng)業(yè)面源污染物的分布情況;王思雅等[13]基于MIKEZERO2020軟件建立水文、水質(zhì)信息耦合模型,對郁江貴港河段水位、流量以及污染物在水體中的遷移、擴(kuò)散、衰減情況進(jìn)行模擬。這些研究通過MIKE模型在極大程度上對流域內(nèi)污染物傳輸過程進(jìn)行了模擬,但沒有綜合考慮地形地貌、土地利用類型等因素對模型精準(zhǔn)度造成的影響。

    因此,該研究以南湖流域典型支流王家河流域?yàn)檠芯繀^(qū)域,在開展水環(huán)境磷污染物監(jiān)測的基礎(chǔ)上,進(jìn)一步考慮流域水文環(huán)境、地形地貌、土地利用類型等因素的影響,通過MIKE21軟件構(gòu)建多尺度、環(huán)境信息與空間遙感融合的磷素源-匯傳輸模型并對模型進(jìn)行優(yōu)化,對磷素面源污染負(fù)荷的時(shí)空變化進(jìn)行分析與模擬,以期建立洞庭湖流域磷污染物精準(zhǔn)溯源和傳輸路徑高分辨解析范式,揭示水體磷素面源污染負(fù)荷的時(shí)空演變規(guī)律,為水體磷污染精準(zhǔn)防控提供依據(jù)。

    1 材料與方法

    1.1 南湖流域概況

    南湖所處的岳陽市屬中亞熱帶向北亞熱帶過渡的濕潤季風(fēng)性氣候。南湖因南津港大堤的修筑而與洞庭湖相分隔,僅在出口處有1個(gè)與洞庭湖相通的控制性閘口,成為典型的半封閉城市內(nèi)湖泊,湖水補(bǔ)給主要依賴湖面降水、集雨區(qū)徑流水和城市污水。南湖現(xiàn)有水面面積約為 13.78km2 ,湖岸線總長約為28.2km[14] ,且均屬于城市區(qū)域;南湖平均水深約為 3m 最大水深約為9m,集雨面積約為150km2[1,15]。 。

    1.2 試驗(yàn)設(shè)計(jì)

    1.2.1研究區(qū)域及監(jiān)測點(diǎn)布設(shè)課題組于2023年9月—2024年8月在湖南省岳陽市岳陽樓區(qū)南湖流域開展水體磷素污染情況監(jiān)測,以南湖流域的主要支流之一——王家河流域?yàn)檠芯繀^(qū)域,將監(jiān)測點(diǎn)主要布設(shè)在王家河和匯水區(qū)及其周邊區(qū)域。王家河位于南湖流域東北區(qū)域,河流由東北向西南匯入南湖;匯水區(qū)位于南湖流域東部區(qū)域,與王家河下游相連。如圖1所示,該研究根據(jù)王家河兩側(cè)和匯水區(qū)污染源類型及排污口情況共布設(shè)8個(gè)污染源監(jiān)測點(diǎn),分別為南津古渡溢流口(S1)、木橋市政生活污水排口(S2)王家橋河小區(qū)溢流口(S3)、錦繡河山小區(qū)溢流口(S4)、牛欄坡排口(S5)、王家橋河長山臨時(shí)泵站排口(S6)、吉祥園林塘溢流口(S7)、畔湖灣社區(qū)排污口(S8);根據(jù)王家河的長度和水文情況劃分采樣斷面、設(shè)置采樣垂線,從王家河的上游到下游加密布點(diǎn),沿途共布設(shè)25個(gè)排放路徑監(jiān)測點(diǎn),分別為W1\~W25;根據(jù)匯水區(qū)情況,在匯水區(qū)布設(shè)20個(gè)加密監(jiān)測點(diǎn),分別為H1\~H20。

    1.2.2樣品采集及分析測定該研究中的樣品采集、運(yùn)輸、保存及質(zhì)量控制方法按照HJ/T91—2002地表水和污水監(jiān)測技術(shù)規(guī)范、HJ495—2009水質(zhì)采樣方案設(shè)計(jì)技術(shù)規(guī)定和SL219—2013水環(huán)境監(jiān)測規(guī)范等執(zhí)行。該研究采用GB11893—89水質(zhì)總磷的測定鉬酸銨分光光度法測定水體總磷,將樣品檢測統(tǒng)一為3次重復(fù),采用Excel、Origin、MIKE21、Arcmap軟件進(jìn)行數(shù)據(jù)計(jì)算與統(tǒng)計(jì)分析。

    1.3 水質(zhì)模型模擬

    1.3.1 水動力模塊模擬 水動力模塊基于三向不可壓縮和Reynolds值均布的Navier-Stokes方程,并服從于Boussinesq假定和靜水壓力的假定。此即流體低速流動中密度的變化分析僅考慮溫度變化產(chǎn)生的影響,而不考慮壓強(qiáng)變化產(chǎn)生的影響。二維非恒定平面淺水流方程組如公式(1)~(3)所示[16-19]。

    圖1南湖流域概況及監(jiān)測點(diǎn)布設(shè)情況(該圖來源于91衛(wèi)圖助手軟件)

    式 Ψ(1)~Ψ(3) 中: t 為時(shí)間; η 為水位; h 為總水深; u,ν 分別為 x,y 方向上的速度分量; g 為重力加速度; ρ 為水的密度, ρ0 為相對密度; Pa 為大氣壓強(qiáng); Sxx 、 Sxy 、 Syx 、 Syy 為輻射應(yīng)力分量; Txx 、Txy,Tyy 為粘滯應(yīng)力; s 為源項(xiàng); us,νs 為源項(xiàng)水流速度。

    1.3.2水質(zhì)模塊模擬濃度對流擴(kuò)散方程可以反映水中模擬對象2個(gè)方面的輸移機(jī)理,即由平均水流速度引起的對流輸移機(jī)理和由濃度梯度引起的擴(kuò)散(彌散)輸移機(jī)理。污染物在自由表面的二維水體中的對流擴(kuò)散輸移過程如公式(4)所示[20]。

    式(4)中: C 是水體污染物濃度; Dx,Dy 分別是 x,y 方向上的擴(kuò)散系數(shù); F 是線性衰減系數(shù);s 為源匯項(xiàng),即污染物的源頭和匯合。

    1.3.3數(shù)據(jù)來源該研究所需的模型庫數(shù)據(jù)主要分為屬性數(shù)據(jù)和空間數(shù)據(jù)2個(gè)部分。屬性數(shù)據(jù)包括水質(zhì)數(shù)據(jù)和農(nóng)業(yè)管理數(shù)據(jù),其中,水質(zhì)數(shù)據(jù)為實(shí)際監(jiān)測數(shù)據(jù),農(nóng)業(yè)管理數(shù)據(jù)來源于《湖南統(tǒng)計(jì)年鑒2024》??臻g數(shù)據(jù)包括土地利用類型、數(shù)字高程模型(DigitalElevationModel,DEM)、坡度、坡向等參數(shù),其中,土地利用數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)平臺,分辨率為 1 000m×1 000m ,DEM、坡度、坡向數(shù)據(jù)來源于地理空間數(shù)據(jù)云平臺。

    該研究在王家河流域水文環(huán)境信息和監(jiān)測點(diǎn)實(shí)測數(shù)據(jù)的基礎(chǔ)上,通過MIKE21、Arcmap軟件構(gòu)建多尺度、環(huán)境信息與空間遙感融合的磷素源-匯傳輸水質(zhì)模型,對王家河及匯水區(qū)的水質(zhì)變化特征進(jìn)行系統(tǒng)模擬分析;并通過引入整合岳陽樓區(qū)土地利用類型、高程、坡度及坡向等多源數(shù)據(jù)對模型進(jìn)行優(yōu)化、提高模型精準(zhǔn)度,對模型模擬數(shù)據(jù)與實(shí)際監(jiān)測數(shù)據(jù)進(jìn)行對比分析。

    2 結(jié)果與分析

    2.1 王家河污染源定點(diǎn)磷素監(jiān)測情況

    如圖2所示,2023年9月—2024年8月S3、S4、S5、S6監(jiān)測點(diǎn)平均總磷濃度接近GB18918—2002城鎮(zhèn)污水處理廠污染物排放標(biāo)準(zhǔn)中一級標(biāo)準(zhǔn)的A標(biāo)準(zhǔn)( 1mg/L ),且總磷濃度最高值大多數(shù)出現(xiàn)在冬季(1一2月,12月的數(shù)據(jù)因未監(jiān)測而沒有在圖中呈現(xiàn))。S1、S2、S7、S8監(jiān)測點(diǎn)平均總磷濃度遠(yuǎn)低于該標(biāo)準(zhǔn)。8個(gè)監(jiān)測點(diǎn)的磷素濃度高峰期主要集中在冬季和春季,其中,S3、S6在冬季時(shí)平均總磷排放濃度最高,最高值為 1.3mg/L 。

    2.2王家河各季節(jié)磷素加密監(jiān)測情況

    該研究對王家河沿途進(jìn)行實(shí)際監(jiān)測,得到王家河W1\~W25各季節(jié)水體平均總磷濃度,并使用Arcmap軟件對結(jié)果進(jìn)行繪圖(圖3)。如圖3所示,王家河W1\~W25秋季平均總磷濃度為 0.62mg/L 冬季平均總磷濃度為 0.65mg/L ,春季平均總磷濃度為 1.00mg/L ,夏季平均總磷濃度為 0.84mg/L 。根據(jù)GB3838—2002地表水環(huán)境質(zhì)量標(biāo)準(zhǔn),王家河水質(zhì)總體而言超過V類水質(zhì)濃度標(biāo)準(zhǔn)( 0.4mg/L )。在時(shí)間變化上,王家河各監(jiān)測點(diǎn)春季、夏季的水體平均總磷濃度明顯高于秋季、冬季,具體表現(xiàn)為春季 gt; 夏季 gt; 冬季 gt; 秋季。在空間變化上,王家河從上游到下游的平均總磷濃度呈先升高后降低的變化趨勢:王家河上游(W1\~W9)因沒有排污口的布設(shè)而水體平均總磷濃度較低;中游(W10\~W20)的排污口相對集中,導(dǎo)致水體平均總磷濃度逐漸升高;下游(W21\~W25)的水體平均總磷濃度呈下降趨勢,其主要原因是河道的變寬以及各類水生植物的消納降解作用[21-22]

    圖2王家河8個(gè)污染源監(jiān)測點(diǎn)各月份磷素濃度
    圖3王家河各季節(jié)磷素加密監(jiān)測情況

    2.3匯水區(qū)各季節(jié)磷素加密監(jiān)測情況

    該研究對匯水區(qū)沿途進(jìn)行實(shí)際監(jiān)測,得到匯水區(qū)H1~H20 各季節(jié)水體平均總磷濃度,并使用Arcmap軟件對結(jié)果進(jìn)行繪圖(圖4。從圖4可以看出,匯水區(qū) H1~H20 水體秋季平均總磷濃度為 0.42mg/L ,冬季平均總磷濃度為 0.71mg/L ,春季平均總磷濃度為 0.74mg/L ,夏季平均總磷濃度為 0.57mg/L 。根據(jù)GB3838—2002地表水環(huán)境質(zhì)量標(biāo)準(zhǔn),匯水區(qū)水質(zhì)總體而言超過V類水質(zhì)濃度標(biāo)準(zhǔn)。在時(shí)間變化上,匯水區(qū)水體春季、冬季平均總磷濃度明顯高于夏季、秋季,具體表現(xiàn)為春季 gt; 冬季 gt; 夏季 gt; 秋季。在空間變化上,春季時(shí)匯水區(qū)中部(H8\~H15)的水體平均總磷濃度高于兩側(cè)(H1\~H7和H16\~H20),秋季時(shí)則呈現(xiàn)相反趨勢。

    2.4水質(zhì)模型的構(gòu)建、優(yōu)化與精準(zhǔn)度驗(yàn)證

    實(shí)測結(jié)果顯示,2023年9月—2024年8月王家河及匯水區(qū)總磷濃度范圍為 0.18~1.75mg/L ,參照GB3838—2002地表水環(huán)境質(zhì)量標(biāo)準(zhǔn),研究區(qū)域水體的平均總磷濃度超過V類水質(zhì)標(biāo)準(zhǔn)限值。該研究基于王家河流域水文環(huán)境信息及監(jiān)測點(diǎn)實(shí)測數(shù)據(jù),采用MIKE21軟件構(gòu)建水質(zhì)模型,依據(jù)王家河流域監(jiān)測斷面磷污染物濃度實(shí)測值設(shè)定模型初始水質(zhì)濃度,經(jīng)模型計(jì)算后得出王家河及匯水區(qū)磷素空間分布結(jié)果。為了進(jìn)一步提升模型精準(zhǔn)度,該研究通過引入整合岳陽樓區(qū)土地利用類型、高程、坡度及坡向等多源數(shù)據(jù),對王家河及匯水區(qū)的水質(zhì)變化特征模擬輸出結(jié)果進(jìn)行系統(tǒng)優(yōu)化。模型模擬結(jié)果顯示,研究區(qū)域總磷濃度范圍為 0.03~1.80mg/L ,與實(shí)際監(jiān)測結(jié)果較為相近,且其濃度變化趨勢與實(shí)際監(jiān)測結(jié)果基本一致,這說明模型模擬結(jié)果較好地揭示了王家河及匯水區(qū)磷素的空間分布特征。通過將模型模擬值與實(shí)測值進(jìn)行對比發(fā)現(xiàn),優(yōu)化后的模型對水體磷素含量的預(yù)測精準(zhǔn)度較優(yōu)化前有明顯改善,平均精準(zhǔn)度為 70% 左右(秋、冬、春、夏季 R2 分別為0.758、0.720、0.572、0.869)。上述結(jié)果表明,優(yōu)化后的模型模擬預(yù)測精準(zhǔn)度顯著提升,且模擬預(yù)測結(jié)果與實(shí)際監(jiān)測數(shù)據(jù)之間具有良好的一致性。

    圖4匯水區(qū)各季節(jié)磷素加密監(jiān)測情況

    3 討論

    研究結(jié)果表明,2023年9月—2024年8月南湖王家河8個(gè)污染源監(jiān)測點(diǎn)的磷素濃度高峰期主要集中在冬季和春季,其中,S3、S6在冬季時(shí)平均總磷排放濃度最高,最高值為 1.3mg/L 。王家河的水體平均總磷濃度為 0.78mg/L ,匯水區(qū)的水體平均總磷濃度為 0.61mg/L ,根據(jù)GB3838—2002地表水環(huán)境質(zhì)量標(biāo)準(zhǔn),其均超過 ΔV 類水質(zhì)濃度標(biāo)準(zhǔn)??梢?,王家河的水體平均總磷濃度高于匯水區(qū),這可能與支流對中心湖面有一定的保護(hù)作用有關(guān)[19]。在時(shí)間變化上,王家河和匯水區(qū)的水體平均總磷濃度變化雖有不同,但均表現(xiàn)為春季最高、秋季最低。這與當(dāng)?shù)鼐用裆钅J降募竟?jié)性差異有較大關(guān)系,由于人類活動強(qiáng)度、氣候條件等不確定性因素的影響,水體污染物濃度呈現(xiàn)出顯著的季節(jié)性波動特征[21-23];水產(chǎn)養(yǎng)殖、禽畜養(yǎng)殖污染等也是造成水體污染的重要原因,這與已有研究的結(jié)論一致[24-28]。在空間變化上,王家河從上游到下游(W1\~W25)的水體平均總磷濃度呈先升高后降低的趨勢;春季時(shí)匯水區(qū)中部(H8\~H15)的水體平均總磷濃度高于兩側(cè)(H1\~H7和H16\~H20),秋季時(shí)則呈現(xiàn)相反趨勢。總體來看,王家河與匯水區(qū)的水體平均總磷濃度與季節(jié)變化及居民活動密切相關(guān)。

    此外,該研究構(gòu)建并優(yōu)化多尺度、環(huán)境信息與空間遙感融合的磷素源-匯傳輸模型,對王家河流域水體總磷濃度進(jìn)行模擬預(yù)測,發(fā)現(xiàn)模擬預(yù)測結(jié)果與實(shí)際監(jiān)測結(jié)果趨勢較為一致,模型模擬預(yù)測的平均精準(zhǔn)度為 70% 左右。模型模擬精準(zhǔn)度總體較高的原因是污染源監(jiān)測數(shù)據(jù)、流域水文環(huán)境信息、土地利用類型、高程、坡度及坡向等與研究區(qū)域水質(zhì)變化密切相關(guān),王家河與匯水區(qū)周邊大多為城市用地和農(nóng)田,工業(yè)廢水、居民生活污水、農(nóng)業(yè)污水等大多排入南湖流域。值得注意的是,春季的模擬數(shù)據(jù)與實(shí)測數(shù)據(jù)之間存在較大偏差( R2=0.572 ),這一現(xiàn)象可能主要?dú)w因于季節(jié)性環(huán)境變化以及研究區(qū)域周邊居民生活模式的季節(jié)性差異。

    參考文獻(xiàn):

    [1]朱江,林小莉.湖泊濕地生態(tài)修復(fù)規(guī)劃研究:以岳陽南湖濕地生態(tài)修復(fù)為例[J].濕地科學(xué)與管理,2020,16(3):12-16.

    [2]宋麗,符海文.岳陽南湖水環(huán)境現(xiàn)狀及綜合治理措施探討[J].廣東化工,2013,40(14):136-137.

    [3]王麗婧,田澤斌,李瑩杰,等.洞庭湖近30年水環(huán)境演變態(tài)勢及影響因素研究[J].環(huán)境科學(xué)研究,2020,33(5):1140-1149.

    [4]彭嫻,朱丹丹,熊劍.洞庭湖水環(huán)境質(zhì)量狀況與污染來源影響分析[J].湖南理工學(xué)院學(xué)報(bào)(自然科學(xué)版),2023,36(3):50-56.

    [5]楊文,周腳根,張文釗.洞庭湖源頭流域典型塘庫水體氮磷組分、pH、DO及水溫晝夜變異特征[J].長江流域資源與環(huán)境,2023,32(2):384-393.

    [6]符哲,郭晶,黃代中,等.洞庭湖的富營養(yǎng)演變特征及影響因素[J].環(huán)境化學(xué),2022,41(8):2636-2645.

    [7]王博,劉大喜,羅朝,等.我國水域的磷污染及其治理[J].生態(tài)經(jīng)濟(jì),2024,40(9):5-8.

    [8]崔鍵,杜易,丁程成,等.中國湖泊水體磷的賦存形態(tài)及污染治理措施進(jìn)展[J].生態(tài)環(huán)境學(xué)報(bào),2022,31(3):621-633.

    [9]李勝男,李冀,何康,等.洞庭湖區(qū)滄浪河流域農(nóng)業(yè)面源污染現(xiàn)狀調(diào)查與分析[J].湖南農(nóng)業(yè)科學(xué),2018(5):56-60.

    [10]WANGFE,LIYN,YANGJ,etal.ApplicationofWASP modeland Gini coefficient in total mass control of water pollutants ∵ AcasestudyinXichengCanal,China[J].Desalinationand WaterTreatment,2016,57(7):2903-2916.

    [11]VENISHETTY V,PARAJULI PB.Assessment of BMPs byestimatinghydrologic andwaterqualityoutputsusingSWATinYazoo River Watershed[J].Agriculture,2022,12: 477.

    [12] LU X,DONGY,LIUQ,et al. Simulation on TNand TPdistributionofsedimentinLiaoheEstuaryNationalWetlandParkusingMIKE21-couplingmodel[J].Water,2023,15(15):2727.

    [13]王思雅,呂俊.基于MIKEZERO2020的郁江貴港河段一維水質(zhì)模型[J].廣西水利水電,2024(5):130-137.

    [14]南湖新區(qū)辦公室.南湖概況[J].(2025-04-27)[2025-05-07].https://yynanhu.yueyang.gov.cn/29503/29504/content_1609666.html.

    [15]覃鴻,楊偉,黃菊梅,等.氣象因子對岳陽市南湖水質(zhì)富營養(yǎng)化進(jìn)程的影響與關(guān)系研究[J].湖北農(nóng)業(yè)科學(xué),2021,60(S2):152-155,251.

    [16]王愿珠,葛長成.基于Mike21模型的鴨綠江支流排污口水質(zhì)影響分析[J].濕地科學(xué)與管理,2024,20(5):59-64.

    [17]耿玉瑩,錢秀紅,金炳琪,等.基于MIKE21的小清河感潮河段水質(zhì)數(shù)值模擬研究[J].人民黃河,2024,46(10):96-102.

    [18]高濤,黃偉杰,吳瓊,等.Mike21模型在河口大型綜合工程洪水影響評價(jià)方面的應(yīng)用[J].城市建設(shè)理論研究(電子版),2024(22):196-198.

    [19]WANGD,TAOT,ZHANG S,et al.Tributaries protect the waterqualityofcentral lake inanartificial spiderweb-structured watersystem,from the perspectiveof runoff pollution control[J]. JournalofCleanerProduction,2023,429:139566.

    [20]YEY,ZHANGJ,LIUH,etal.Predictive simulationstudyonthe effect of smalland medium river basin outfall treatment measuresonwaterqualityimprovement[J].Water,2023,15(13):2359.

    [21]韋必穎,成建梅,蘇曉煜,等.深圳市2015—2021年雨源型河流水質(zhì)時(shí)空變化及其對降雨的響應(yīng)[J].環(huán)境科學(xué),2024,45(2):780-791.

    [22]李昂,劉加強(qiáng),劉超,等.城市降雨及雨水徑流對河道水質(zhì)的影響[J].供水技術(shù),2022,16(3):49-53.

    [23]王登超,李發(fā)東,李曹樂,等.降雨對河流水質(zhì)影響及污染源解析:以鑒江茂名段為例[J].環(huán)境科學(xué),2025,46(4):2165-2178.

    [24]談毅奇,李向前,殷文潔,等.岳陽縣東洞庭湖自然保護(hù)區(qū)實(shí)驗(yàn)區(qū)水產(chǎn)養(yǎng)殖污染現(xiàn)狀及綜合治理對策研究[J].當(dāng)代水產(chǎn),2020,45(8):80-81.

    [25]鄭立國,王海萍,李靈,等.湖南省化肥農(nóng)藥使用的時(shí)空分布特征研究[J].綠色科技,2024,26(16):208-212.

    [26]王風(fēng)燕.畜牧養(yǎng)殖環(huán)境污染的現(xiàn)狀及治理[J].畜牧業(yè)環(huán)境,2023(2):35-36.

    [27]張瑞姿.水產(chǎn)養(yǎng)殖環(huán)境污染問題與治理措施[J].農(nóng)家參謀,2022(3):162-164.

    [28]蒙小俊,張建東,王秋利,等.漢江流域安康段面源污染特征分析[J].湖南農(nóng)業(yè)科學(xué),2023(10):35-42,49.

    (責(zé)任編輯:袁萍萍)

    猜你喜歡
    磷素南湖監(jiān)測點(diǎn)
    CLDAS土壤相對濕度產(chǎn)品在皖北地區(qū)適應(yīng)性評估及應(yīng)用
    2021年-2023年怒江州草地貪夜蛾成蟲和幼蟲發(fā)生分析
    南湖畫舫制作技藝
    氮用量和基追肥氮比例對稻茬烤煙氮、磷、鉀積累和產(chǎn)量、產(chǎn)值的影響
    山水筆記
    優(yōu)化營商環(huán)境的晉州實(shí)踐
    公民與法治(2025年7期)2025-08-05 00:00:00
    城市園林景觀設(shè)計(jì)對河道水污染的防治效果研究
    99久久人妻综合| 欧美日本视频| 日韩欧美 国产精品| 免费播放大片免费观看视频在线观看| 亚洲av电影在线观看一区二区三区 | 欧美精品国产亚洲| 好男人在线观看高清免费视频| 久久草成人影院| 日韩成人伦理影院| 国产一区二区三区av在线| 亚洲熟妇中文字幕五十中出| 成年人午夜在线观看视频 | 亚洲av中文av极速乱| 黑人高潮一二区| 亚洲国产精品成人综合色| 干丝袜人妻中文字幕| 国产精品人妻久久久久久| 男的添女的下面高潮视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 性插视频无遮挡在线免费观看| 秋霞在线观看毛片| 久久亚洲国产成人精品v| 亚洲美女搞黄在线观看| 在线a可以看的网站| 男女下面进入的视频免费午夜| 青青草视频在线视频观看| 午夜激情欧美在线| 精品一区二区三卡| 成年女人看的毛片在线观看| 亚洲18禁久久av| 在线a可以看的网站| av网站免费在线观看视频 | 高清毛片免费看| 十八禁网站网址无遮挡 | 国产在线一区二区三区精| 人妻一区二区av| 亚洲精品日韩在线中文字幕| 久久综合国产亚洲精品| 亚洲av不卡在线观看| 亚洲在久久综合| 久久99蜜桃精品久久| 最近视频中文字幕2019在线8| 国产黄色视频一区二区在线观看| 欧美日本视频| 亚洲性久久影院| 2022亚洲国产成人精品| 日本av手机在线免费观看| 成年免费大片在线观看| 国产探花极品一区二区| 婷婷色综合大香蕉| 六月丁香七月| videos熟女内射| 色网站视频免费| 菩萨蛮人人尽说江南好唐韦庄| 丝袜美腿在线中文| 日韩欧美国产在线观看| 两个人的视频大全免费| 欧美性猛交╳xxx乱大交人| 99热6这里只有精品| 国产午夜精品一二区理论片| 国产一区二区三区av在线| 国产黄色小视频在线观看| 又爽又黄无遮挡网站| 亚洲va在线va天堂va国产| 伊人久久国产一区二区| 肉色欧美久久久久久久蜜桃 | 国内精品宾馆在线| a级一级毛片免费在线观看| 国内少妇人妻偷人精品xxx网站| 精品久久久久久久久av| 麻豆精品久久久久久蜜桃| 亚洲欧美成人精品一区二区| eeuss影院久久| 丝袜美腿在线中文| videos熟女内射| 亚洲电影在线观看av| 亚洲一级一片aⅴ在线观看| 春色校园在线视频观看| a级一级毛片免费在线观看| 亚洲精品视频女| 丰满乱子伦码专区| 99久国产av精品| 国内精品一区二区在线观看| 亚洲一级一片aⅴ在线观看| 国内精品宾馆在线| 免费看不卡的av| 国产极品天堂在线| 亚洲欧美日韩卡通动漫| 真实男女啪啪啪动态图| 亚洲精品乱码久久久久久按摩| 搡老乐熟女国产| 色视频www国产| 亚洲,欧美,日韩| 亚洲成人中文字幕在线播放| 69av精品久久久久久| 亚洲精品乱码久久久v下载方式| 波野结衣二区三区在线| 国产精品三级大全| 如何舔出高潮| 亚洲av二区三区四区| 九九久久精品国产亚洲av麻豆| 欧美xxxx黑人xx丫x性爽| 午夜久久久久精精品| 在线天堂最新版资源| 国产欧美另类精品又又久久亚洲欧美| 久久人人爽人人片av| 尾随美女入室| 久久精品夜夜夜夜夜久久蜜豆| 久久午夜福利片| 久久精品国产鲁丝片午夜精品| 国产精品久久久久久久电影| 97在线视频观看| 国产视频内射| 欧美日韩综合久久久久久| 国产精品一区二区三区四区久久| 久久99热6这里只有精品| 麻豆成人av视频| 在线观看一区二区三区| 在现免费观看毛片| 亚洲在久久综合| 亚洲欧洲国产日韩| 午夜爱爱视频在线播放| 国内精品一区二区在线观看| 九九久久精品国产亚洲av麻豆| 毛片女人毛片| 麻豆成人av视频| 国内揄拍国产精品人妻在线| 亚洲成人久久爱视频| 91精品一卡2卡3卡4卡| 少妇人妻一区二区三区视频| 高清欧美精品videossex| 天堂av国产一区二区熟女人妻| 蜜臀久久99精品久久宅男| 久久97久久精品| 高清在线视频一区二区三区| 男女啪啪激烈高潮av片| 性插视频无遮挡在线免费观看| 丝袜喷水一区| av在线观看视频网站免费| 日日摸夜夜添夜夜爱| 亚洲精品乱码久久久v下载方式| 久久精品国产亚洲av涩爱| 精品久久久久久久久亚洲| 日韩av在线大香蕉| av黄色大香蕉| 午夜福利在线观看免费完整高清在| 国产欧美另类精品又又久久亚洲欧美| av女优亚洲男人天堂| kizo精华| 国产一区二区在线观看日韩| 午夜福利在线观看吧| 国产精品麻豆人妻色哟哟久久 | 少妇人妻精品综合一区二区| 日韩成人伦理影院| 亚洲精品成人久久久久久| 色播亚洲综合网| 中文在线观看免费www的网站| 久久久久网色| 欧美一级a爱片免费观看看| 日韩,欧美,国产一区二区三区| 日韩一本色道免费dvd| 2018国产大陆天天弄谢| 国产一区二区亚洲精品在线观看| xxx大片免费视频| 夫妻午夜视频| 日韩强制内射视频| 深爱激情五月婷婷| 日本av手机在线免费观看| 国产黄片视频在线免费观看| 日本av手机在线免费观看| 69av精品久久久久久| videossex国产| 国产老妇伦熟女老妇高清| 精品一区二区三区视频在线| 97精品久久久久久久久久精品| 成年免费大片在线观看| 看黄色毛片网站| 欧美区成人在线视频| 又大又黄又爽视频免费| 久久韩国三级中文字幕| 亚洲精品第二区| 久久这里有精品视频免费| 成人亚洲欧美一区二区av| 亚洲欧洲国产日韩| 一区二区三区免费毛片| 十八禁国产超污无遮挡网站| 国产精品1区2区在线观看.| av在线观看视频网站免费| 亚洲一级一片aⅴ在线观看| 国产黄片视频在线免费观看| 国产乱来视频区| 久久久久免费精品人妻一区二区| 精品一区二区三区人妻视频| 日本wwww免费看| 亚洲欧美日韩卡通动漫| 女的被弄到高潮叫床怎么办| 精品少妇黑人巨大在线播放| 人人妻人人澡欧美一区二区| 噜噜噜噜噜久久久久久91| 国内精品宾馆在线| 高清视频免费观看一区二区 | 直男gayav资源| av网站免费在线观看视频 | kizo精华| 国产亚洲91精品色在线| 亚洲精品一区蜜桃| 亚洲精品日韩av片在线观看| 亚洲国产精品成人久久小说| 人妻制服诱惑在线中文字幕| 久久久精品94久久精品| 极品少妇高潮喷水抽搐| 国产一区二区在线观看日韩| 一级黄片播放器| 三级男女做爰猛烈吃奶摸视频| 久久久国产一区二区| 免费看日本二区| 国内精品美女久久久久久| 免费观看性生交大片5| 成人美女网站在线观看视频| 精品一区二区三区人妻视频| 五月伊人婷婷丁香| 日本av手机在线免费观看| 久热久热在线精品观看| 亚洲最大成人中文| 精品久久久久久电影网| 亚洲av免费高清在线观看| 秋霞伦理黄片| 免费高清在线观看视频在线观看| 亚洲精品久久午夜乱码| 国产综合精华液| 欧美日韩视频高清一区二区三区二| 久久人人爽人人片av| av线在线观看网站| 在线a可以看的网站| 国产欧美另类精品又又久久亚洲欧美| 毛片一级片免费看久久久久| 免费黄频网站在线观看国产| 国产永久视频网站| 精品久久久久久久久av| 国精品久久久久久国模美| 亚洲av电影不卡..在线观看| h日本视频在线播放| 日韩伦理黄色片| 一区二区三区免费毛片| 成人美女网站在线观看视频| 成人毛片60女人毛片免费| 夫妻午夜视频| 麻豆精品久久久久久蜜桃| 麻豆成人av视频| 久久精品国产自在天天线| 性插视频无遮挡在线免费观看| 天堂影院成人在线观看| 日本-黄色视频高清免费观看| 国产亚洲午夜精品一区二区久久 | 乱系列少妇在线播放| 国产黄片美女视频| 国产黄片视频在线免费观看| 久久久久免费精品人妻一区二区| 丰满乱子伦码专区| 久久久a久久爽久久v久久| 深夜a级毛片| 午夜爱爱视频在线播放| 精品一区二区三区人妻视频| 天天躁日日操中文字幕| 国产成人freesex在线| 床上黄色一级片| 极品少妇高潮喷水抽搐| 国产伦精品一区二区三区视频9| 韩国av在线不卡| 免费观看在线日韩| 美女内射精品一级片tv| 99久久九九国产精品国产免费| 激情 狠狠 欧美| 日韩视频在线欧美| 国产欧美另类精品又又久久亚洲欧美| videos熟女内射| 美女主播在线视频| 亚洲一区高清亚洲精品| av在线观看视频网站免费| 久久久久网色| 国产免费福利视频在线观看| 婷婷色综合www| 岛国毛片在线播放| 亚洲精品久久午夜乱码| 特级一级黄色大片| 亚洲av一区综合| 可以在线观看毛片的网站| 国产精品一及| 精品人妻一区二区三区麻豆| 久久久成人免费电影| 男人狂女人下面高潮的视频| 精品一区二区三区人妻视频| or卡值多少钱| 国产精品嫩草影院av在线观看| 国国产精品蜜臀av免费| 一级毛片电影观看| 日韩欧美精品免费久久| 亚洲国产精品国产精品| 三级经典国产精品| 亚洲国产欧美在线一区| 男人狂女人下面高潮的视频| 欧美最新免费一区二区三区| 亚洲内射少妇av| 性色avwww在线观看| 91精品国产九色| av.在线天堂| 在线观看美女被高潮喷水网站| 亚洲在线观看片| av国产免费在线观看| 好男人在线观看高清免费视频| 亚洲在线观看片| 日韩av不卡免费在线播放| 成人漫画全彩无遮挡| 99久久精品国产国产毛片| 人妻制服诱惑在线中文字幕| 精品久久久久久电影网| 亚洲精品成人av观看孕妇| 91狼人影院| 人人妻人人澡欧美一区二区| 亚洲国产精品成人综合色| 看黄色毛片网站| 日韩三级伦理在线观看| 成人av在线播放网站| 三级经典国产精品| 五月玫瑰六月丁香| 一级av片app| 在线a可以看的网站| 最近最新中文字幕免费大全7| 免费看日本二区| 亚洲国产高清在线一区二区三| 嫩草影院精品99| 亚洲图色成人| 国产精品国产三级专区第一集| 卡戴珊不雅视频在线播放| 中文在线观看免费www的网站| 日本午夜av视频| 在线 av 中文字幕| 亚洲三级黄色毛片| 久久久久久伊人网av| 99热6这里只有精品| 永久网站在线| av在线天堂中文字幕| 国产成人精品一,二区| 亚洲乱码一区二区免费版| 青青草视频在线视频观看| av线在线观看网站| 2021天堂中文幕一二区在线观| 国产亚洲精品久久久com| 少妇高潮的动态图| 久久国产乱子免费精品| 欧美日韩国产mv在线观看视频 | 国产精品无大码| 日日干狠狠操夜夜爽| 欧美97在线视频| 在线观看一区二区三区| 国产亚洲精品久久久com| 中文资源天堂在线| 亚洲国产精品sss在线观看| 亚洲国产精品成人综合色| 联通29元200g的流量卡| 欧美bdsm另类| 国产精品一区www在线观看| 免费大片黄手机在线观看| 国产精品一区二区三区四区久久| 特大巨黑吊av在线直播| 中文天堂在线官网| 卡戴珊不雅视频在线播放| 99久久精品一区二区三区| 在线观看一区二区三区| 久久精品夜色国产| 欧美日韩综合久久久久久| 欧美性感艳星| 国产淫语在线视频| 婷婷色综合大香蕉| 国产精品一区www在线观看| 久久久久久久久久久免费av| 欧美性猛交╳xxx乱大交人| 高清毛片免费看| 日韩精品有码人妻一区| 亚洲精品乱码久久久v下载方式| 草草在线视频免费看| 久久精品国产亚洲av天美| 99久久精品热视频| 国产乱人偷精品视频| 亚洲av男天堂| 真实男女啪啪啪动态图| 日韩三级伦理在线观看| 国产成人精品一,二区| 中文精品一卡2卡3卡4更新| 青春草视频在线免费观看| 有码 亚洲区| 波多野结衣巨乳人妻| 国产高清国产精品国产三级 | 熟妇人妻不卡中文字幕| 久久久午夜欧美精品| 一级毛片我不卡| 成人综合一区亚洲| 91午夜精品亚洲一区二区三区| 激情 狠狠 欧美| 亚洲av在线观看美女高潮| freevideosex欧美| 大香蕉97超碰在线| 免费大片黄手机在线观看| 国产一区二区亚洲精品在线观看| 99久久人妻综合| 日本与韩国留学比较| 五月天丁香电影| .国产精品久久| 综合色av麻豆| 免费黄频网站在线观看国产| 精品人妻视频免费看| 亚洲成人久久爱视频| 美女脱内裤让男人舔精品视频| 天美传媒精品一区二区| 日韩欧美一区视频在线观看 | 久久97久久精品| 国产黄片美女视频| 一区二区三区乱码不卡18| 精品熟女少妇av免费看| 亚洲国产精品成人久久小说| 亚洲精品亚洲一区二区| 美女高潮的动态| 看非洲黑人一级黄片| 久久精品国产鲁丝片午夜精品| 亚洲一级一片aⅴ在线观看| 大又大粗又爽又黄少妇毛片口| 激情五月婷婷亚洲| 久久久久久久久久成人| 一二三四中文在线观看免费高清| 国产精品1区2区在线观看.| 高清欧美精品videossex| 一区二区三区高清视频在线| 久久人人爽人人片av| 深夜a级毛片| 国产亚洲一区二区精品| 在线观看av片永久免费下载| 亚洲精品影视一区二区三区av| 97超碰精品成人国产| 高清av免费在线| 国产探花在线观看一区二区| 黄片wwwwww| 久热久热在线精品观看| 亚洲欧美日韩无卡精品| 欧美3d第一页| 午夜日本视频在线| 美女cb高潮喷水在线观看| 黄色配什么色好看| 亚洲av电影在线观看一区二区三区 | 啦啦啦中文免费视频观看日本| 国产亚洲精品久久久com| 色综合色国产| 校园人妻丝袜中文字幕| 欧美高清成人免费视频www| 久久久国产一区二区| 国产大屁股一区二区在线视频| 免费播放大片免费观看视频在线观看| 亚洲成人一二三区av| 男的添女的下面高潮视频| 亚洲av成人精品一区久久| 中文精品一卡2卡3卡4更新| 又黄又爽又刺激的免费视频.| 亚洲精品日韩av片在线观看| 成人综合一区亚洲| av.在线天堂| 偷拍熟女少妇极品色| 日韩欧美国产在线观看| 免费高清在线观看视频在线观看| 69av精品久久久久久| 中文资源天堂在线| 一边亲一边摸免费视频| 久久久a久久爽久久v久久| 五月伊人婷婷丁香| 欧美丝袜亚洲另类| 亚洲av日韩在线播放| 午夜亚洲福利在线播放| 国产在线男女| 七月丁香在线播放| 免费黄网站久久成人精品| 国产乱人视频| 精品人妻熟女av久视频| 国内少妇人妻偷人精品xxx网站| 国产片特级美女逼逼视频| 亚洲真实伦在线观看| 亚洲美女搞黄在线观看| 2018国产大陆天天弄谢| 亚洲最大成人中文| 亚洲av免费高清在线观看| 人人妻人人澡欧美一区二区| 国产乱人偷精品视频| 有码 亚洲区| 日韩国内少妇激情av| 欧美成人一区二区免费高清观看| 日韩一本色道免费dvd| 亚洲国产成人一精品久久久| 国产男人的电影天堂91| 亚洲,欧美,日韩| 最后的刺客免费高清国语| 91午夜精品亚洲一区二区三区| 亚洲国产成人一精品久久久| 午夜激情久久久久久久| 欧美激情国产日韩精品一区| 亚洲欧美清纯卡通| 只有这里有精品99| 中文字幕人妻熟人妻熟丝袜美| 国国产精品蜜臀av免费| 高清视频免费观看一区二区 | 又爽又黄a免费视频| 成年女人在线观看亚洲视频 | freevideosex欧美| 毛片女人毛片| 高清午夜精品一区二区三区| av国产久精品久网站免费入址| 青春草视频在线免费观看| av在线亚洲专区| 爱豆传媒免费全集在线观看| 特大巨黑吊av在线直播| 街头女战士在线观看网站| 非洲黑人性xxxx精品又粗又长| 成人毛片a级毛片在线播放| 精品欧美国产一区二区三| 国产精品一区二区性色av| 18禁裸乳无遮挡免费网站照片| 精品国产一区二区三区久久久樱花 | 亚洲国产精品成人综合色| h日本视频在线播放| 一区二区三区免费毛片| 国产老妇女一区| 美女cb高潮喷水在线观看| 大话2 男鬼变身卡| 国产人妻一区二区三区在| 中文字幕av成人在线电影| 国产中年淑女户外野战色| 汤姆久久久久久久影院中文字幕 | 内射极品少妇av片p| 91在线精品国自产拍蜜月| 色综合色国产| 丰满少妇做爰视频| 人妻夜夜爽99麻豆av| 熟女人妻精品中文字幕| 久久99精品国语久久久| 久久草成人影院| 久久精品国产自在天天线| 成年av动漫网址| 亚洲伊人久久精品综合| 国产大屁股一区二区在线视频| 国产成年人精品一区二区| 欧美潮喷喷水| 舔av片在线| 国产成人一区二区在线| 嫩草影院新地址| 女人久久www免费人成看片| 精品久久久久久久末码| 久久久久久久久久人人人人人人| 一区二区三区乱码不卡18| 亚洲欧美一区二区三区黑人 | 天堂√8在线中文| 岛国毛片在线播放| 最近2019中文字幕mv第一页| 亚洲精品成人久久久久久| 亚洲图色成人| av又黄又爽大尺度在线免费看| 欧美+日韩+精品| 老司机影院成人| 男女边摸边吃奶| 成人鲁丝片一二三区免费| 日韩不卡一区二区三区视频在线| 看黄色毛片网站| 一级毛片我不卡| 国产精品无大码| 大又大粗又爽又黄少妇毛片口| 亚洲aⅴ乱码一区二区在线播放| 在线观看av片永久免费下载| 日日啪夜夜爽| 亚洲av电影在线观看一区二区三区 | 亚洲精品aⅴ在线观看| 内地一区二区视频在线| or卡值多少钱| 成年女人看的毛片在线观看| 欧美日韩亚洲高清精品| 久久久精品94久久精品| 国产av码专区亚洲av| 婷婷色av中文字幕| 久久久国产一区二区| 白带黄色成豆腐渣| 亚洲精品成人av观看孕妇| 日韩欧美三级三区| 国内精品一区二区在线观看| 亚洲人与动物交配视频| av免费在线看不卡| 日韩,欧美,国产一区二区三区| 国产成人aa在线观看| 国产成人免费观看mmmm| 不卡视频在线观看欧美| 亚洲欧美一区二区三区国产| 我的女老师完整版在线观看| 少妇熟女aⅴ在线视频| 欧美xxxx性猛交bbbb| 我的女老师完整版在线观看| 国产一区二区三区综合在线观看 | 国产精品久久久久久av不卡| 2022亚洲国产成人精品| 午夜精品一区二区三区免费看| 国产成人a∨麻豆精品| 2021天堂中文幕一二区在线观| 夜夜爽夜夜爽视频| 大片免费播放器 马上看| 亚洲自拍偷在线| 青春草视频在线免费观看| 国产精品久久视频播放| 深夜a级毛片| 精品久久久久久久人妻蜜臀av| 午夜免费激情av| 男人舔女人下体高潮全视频| 欧美日本视频| 在线免费十八禁|