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

    滹沱河沖洪積扇地下水硝酸鹽含量的空間分布特征及影響因素

    2016-07-13 09:37:31喬肖翠江秋楓朱媛媛王山軍鄭丙輝中國環(huán)境科學(xué)研究院環(huán)境基準(zhǔn)與風(fēng)險評估國家重點(diǎn)實(shí)驗(yàn)室國家環(huán)境保護(hù)飲用水水源地保護(hù)重點(diǎn)實(shí)驗(yàn)室北京0002北京科技大學(xué)土木與環(huán)境工程學(xué)院環(huán)境工程系北京00083
    關(guān)鍵詞:洪積扇滹沱河點(diǎn)位

    劉 琰,喬肖翠,江秋楓,昌 盛,朱媛媛,2,王山軍,鄭丙輝*(.中國環(huán)境科學(xué)研究院環(huán)境基準(zhǔn)與風(fēng)險評估國家重點(diǎn)實(shí)驗(yàn)室,國家環(huán)境保護(hù)飲用水水源地保護(hù)重點(diǎn)實(shí)驗(yàn)室,北京0002;2.北京科技大學(xué)土木與環(huán)境工程學(xué)院環(huán)境工程系,北京00083)

    ?

    滹沱河沖洪積扇地下水硝酸鹽含量的空間分布特征及影響因素

    劉琰1,喬肖翠1,江秋楓1,昌盛1,朱媛媛1,2,王山軍1,鄭丙輝1*
    (1.中國環(huán)境科學(xué)研究院環(huán)境基準(zhǔn)與風(fēng)險評估國家重點(diǎn)實(shí)驗(yàn)室,國家環(huán)境保護(hù)飲用水水源地保護(hù)重點(diǎn)實(shí)驗(yàn)室,北京100012;2.北京科技大學(xué)土木與環(huán)境工程學(xué)院環(huán)境工程系,北京100083)

    摘要:為了解華北平原地下水中硝酸鹽的分布狀況及影響因素,對石家莊市滹沱河沖洪積扇四個地下水單元的地下水樣品進(jìn)行采集與分析。結(jié)果表明:在115個樣品中,64.35%的樣品水質(zhì)符合《地下水質(zhì)量標(biāo)準(zhǔn)》(GB 14848—1993)Ⅴ類標(biāo)準(zhǔn),地下水質(zhì)量總體較差,主要的污染指標(biāo)為硝酸鹽和總硬度,超標(biāo)點(diǎn)位主要分布在西北部兩個地下水單元。硝酸鹽的濃度范圍為0.24~440 mg·L-1,其平均濃度與最高濃度分別為GB/T 14848—1993中硝酸鹽Ⅲ類標(biāo)準(zhǔn)值的3倍與22倍,屬于嚴(yán)重污染;不同地下水單元中硝酸鹽濃度差異較大,在研究區(qū)內(nèi)基本呈現(xiàn)從西北到東南、從淺層到深層逐漸降低趨勢。除水文地質(zhì)條件對地下水中硝酸鹽含量有重要影響外,污染源類型及其分布也起到了重要作用。通過相關(guān)性分析,初步認(rèn)為研究區(qū)西北部兩個地下水單元中硝酸鹽污染主要來源于點(diǎn)源類型的生活污水及畜禽養(yǎng)殖廢水,東南部的兩個地下水單元的硝酸鹽污染主要來源于面源類型的農(nóng)業(yè)施肥。

    關(guān)鍵詞:硝酸鹽;地下水;污染特征;空間分布;滹沱河沖洪積扇

    劉琰,喬肖翠,江秋楓,等.滹沱河沖洪積扇地下水硝酸鹽含量的空間分布特征及影響因素[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2016, 35(5):947-954.

    LIU Yan, QIAO Xiao-cui, JIANG Qiu-feng, et al. Spatial distribution and influencing factors of nitrate content in groundwater of alluvial-pluvial fan of Hutuo River[J]. Journal of Agro-Environment Science, 2016, 35(5):947-954.

    華北平原是我國重要的經(jīng)濟(jì)區(qū),人口近1.3億,地下水是其最主要供水水源,占總供水量的69.5%[1-2]。滹沱河沖洪積扇地處華北平原西部,處于山前補(bǔ)給區(qū)的第一個過流區(qū)[3-4],也是南水北調(diào)中線工程地下水調(diào)蓄、地下水回灌的重點(diǎn)區(qū)域[5-6],其地下水質(zhì)量狀況對于整個華北平原中部地區(qū)有著至關(guān)重要的意義。早在20世紀(jì)地下水硝態(tài)氮污染已經(jīng)成為國際上普遍關(guān)注的問題[7],特別是對于以地下水作為主要飲用水源的國家和地區(qū)來說,硝態(tài)氮污染地下水對人體健康的潛在威脅更為嚴(yán)峻。由于近些年來受人類活動的影響,滹沱河沖洪積扇平原區(qū)地下水硝酸鹽含量日益增加:1959年其平均含量為2.35 mg·L-1,1978—1985年期間大部分區(qū)域成為上升區(qū),到1992年,該區(qū)NO-3含量為40~50 mg·L-1的水點(diǎn)已占31.2%,超標(biāo)率為35.5%,其中淺層地下水中氮濃度1986—1991年的平均增長速率為1.85 mg·L-1·a-1[8]。針對該地區(qū)開展地下水硝酸鹽污染特征研究,對于區(qū)域地下水資源評價、地下水污染評價等工作都具有十分重要的參考價值。

    目前關(guān)于滹沱河平原地區(qū)地下水的研究較為廣泛。王翠華等[9]開展了滹沱河平原地下水礦化度和硬度變化分析,結(jié)果表明礦化度和總硬度與地下水位埋深密切相關(guān);葉浩等[10]對滹沱河石家莊段地下水進(jìn)行了污染評價,結(jié)果表明該區(qū)域地下水已經(jīng)受到污染,主要的污染成分為硫酸根、硝酸根、亞硝酸根、化學(xué)需氧量(COD);李亞松等[3]研究了滹沱河沖積平原淺層地下水的有機(jī)污染狀況,結(jié)果表明區(qū)域內(nèi)淺層地下水已經(jīng)受到有機(jī)污染威脅,主要檢出組分為鹵代烴類和氯代苯類;王金哲等[11]對滹沱河流域平原區(qū)人類活動對淺層地下水干擾程度進(jìn)行了量化研究,結(jié)果表明干擾程度在不同時代都表現(xiàn)出從西部向東部逐漸減小的趨勢。以硝酸鹽形式存在的溶解態(tài)氮是地下水中最常見的污染物,許多研究表明硝酸鹽氮含量不斷增加是地下水水質(zhì)惡化的主要表現(xiàn)之一,因此國內(nèi)外都把硝酸鹽氮作為評價地下水水質(zhì)的主要指標(biāo)[12-13]。從目前滹沱河沖洪積平原已開展的研究看,雖然針對該地區(qū)地下水污染特征的相關(guān)研究較多,但是通過在研究區(qū)內(nèi)劃分不同地下水單元而系統(tǒng)研究硝酸鹽分區(qū)污染特征及影響因素的研究還相對較少。

    因此,基于滹沱河沖洪積扇在華北平原所處的重要位置,以及其地下水質(zhì)量狀況對華北平原地下水的重要影響,并結(jié)合該區(qū)域地下水的研究現(xiàn)狀,本研究選取滹沱河沖洪積扇作為研究區(qū),系統(tǒng)分析了地下水中硝酸鹽的分區(qū)污染特征,并對影響硝酸鹽分布的因素及污染來源進(jìn)行了探討,以期為華北平原地下水污染防治工作提供依據(jù)。

    1 材料與方法

    1.1研究區(qū)概況

    研究區(qū)位于太行山前崗南水庫和黃壁莊水庫之間的河谷平原和滹沱河沖洪積扇的頂部及中部地帶,是典型的山前沖洪積平原,為明顯分帶性、邊界條件清楚的典型扇型蓄水構(gòu)造。該區(qū)地勢西高東低,向東傾斜,屬溫帶大陸性季風(fēng)氣候區(qū),根據(jù)石家莊1955—2001年共47年的年降水與平均氣溫?cái)?shù)據(jù),該區(qū)多年平均氣溫13.3℃,年平均降水量534 mm。該區(qū)地下水系統(tǒng)為滹沱河沖洪積扇孔隙淡水系統(tǒng),包氣帶巖性以砂性土為主,降水入滲系數(shù)0.2~0.3。含水層為第四系松散沉積物,巖性以砂石、卵石、粗砂、中細(xì)砂為主,含水層垂向連續(xù)性好,從山前向東部,含水層介質(zhì)顆粒由粗變細(xì),總厚度10~40 m。地下水總體流向是西北向東南,地下水主要補(bǔ)給方式為大氣降水、河渠水入滲和灌溉回歸等,主要排泄方式為人工開采。自20世紀(jì)80年代以來,石家莊地區(qū)農(nóng)業(yè)生產(chǎn)的發(fā)展、人口的增長以及工礦企業(yè)不斷增多,對水資源的需求量越來越大,導(dǎo)致該地區(qū)地下水水位埋深逐年降低,且在局部地區(qū)形成了水位降落漏斗[5]。滹沱河貫穿整個研究區(qū)域,滹沱河河床及河漫灘地帶是較為理想的地下調(diào)蓄滲漏補(bǔ)給源區(qū)[14]。石家莊滹沱河沖洪積扇區(qū)域地下水資源豐富,是華北平原地區(qū)主要的飲用水水源和農(nóng)業(yè)灌溉水源。由于研究區(qū)內(nèi)不同區(qū)域地質(zhì)條件、地下水埋深及含水層介質(zhì)巖性存在較大差異,本研究將該區(qū)域劃分為四個地下水單元,分別為崗黃水庫之間河谷裂隙孔隙水單元(G1)、滹沱河沖洪積扇扇頂部孔隙水單元(G2)、滹沱河沖洪積扇扇中部淺層孔隙水單元(G3)和滹沱河沖洪積扇扇中部深層孔隙水單元(G4)。

    1.2樣品采集

    根據(jù)研究區(qū)內(nèi)地下水單元的劃分,在各單元內(nèi)使用網(wǎng)格法并遵循代表性和均勻性的原則布設(shè)了115個地下水采樣點(diǎn)。采樣點(diǎn)分布示意圖和各單元采樣點(diǎn)數(shù)量分別見圖1和表1。G3和G4所涵蓋的水平區(qū)域完全相同,只是垂向深度不同。

    于2014年10月對研究區(qū)內(nèi)115個點(diǎn)位進(jìn)行水質(zhì)樣品采集。對于地下水樣品,在采集前先清洗井孔,排除井內(nèi)滯水,當(dāng)水質(zhì)分析儀參數(shù)達(dá)到穩(wěn)定后開始采樣,并利用便攜式水質(zhì)分析儀對水體pH、溫度、溶解氧(DO)、電導(dǎo)率(EC)、氧化還原電位(Eh)進(jìn)行現(xiàn)場參數(shù)測定。樣品采集后放置于恒溫4℃的冷藏箱中低溫保存,并于2日內(nèi)送至位于石家莊的中國地質(zhì)科學(xué)院水文地質(zhì)環(huán)境地質(zhì)研究所進(jìn)行測試。

    圖1 滹沱河沖洪積扇地下水采樣點(diǎn)分布圖Figure 1 Location map for groundwater sampling sites in alluvial-pluvial fan of Hutuo River

    G1為崗黃水庫之間河谷裂隙孔隙水單元;G2為滹沱河沖洪積扇扇頂部孔隙水單元;G3為滹沱河沖洪積扇扇中部淺層孔隙水單元;G4為滹沱河沖洪積扇扇中部深層孔隙水單元。下同

    G1-Fissure and pore water unit in the valley between Gangnan Reservoir and Huangbizhuang Reservoir; G2-Pore water unit in the top alluvial-pluvial fan of Hutuo River; G3-Shallow pore water unit in the middle alluvial-pluvial fan of Hutuo River; G4-Deep pore water unit in the middle alluvial-pluvial fan of Hutuo River. The same below

    表1 不同地下水單元采樣點(diǎn)數(shù)量的分布情況Table 1 Distribution of sampling sites in different groundwater units

    1.3測試方法

    采用YSI6600(YSI公司,美國)對pH、溫度、DO、EC、Eh等指標(biāo)進(jìn)行現(xiàn)場測定??傆捕龋ㄒ訡aCO3計(jì))采用地下水質(zhì)檢驗(yàn)方法(DZ/T 0064—1993)中EDTA滴定法測試,HCO-3、CO2-3采用滴定法測定,硫酸鹽、氯離子采用離子色譜法檢測,硝酸鹽(以N計(jì))和亞硝酸鹽(以N計(jì))采用分光光度法檢測,氨氮(以N計(jì))采用納氏試劑比色法測定。

    1.4評價標(biāo)準(zhǔn)

    采用《地下水質(zhì)量標(biāo)準(zhǔn)》(GB/T 14848—1993)對地下水樣品進(jìn)行水質(zhì)評價;采用《生活飲用水衛(wèi)生標(biāo)準(zhǔn)》(GB 5749—2006)對硝酸鹽、亞硝酸鹽和總硬度的達(dá)標(biāo)情況進(jìn)行評價。

    GB 14848—1993和GB 5749—2006中關(guān)于硝酸鹽、亞硝酸鹽和總硬度的標(biāo)準(zhǔn)值見表2。

    1.5數(shù)據(jù)處理

    將各點(diǎn)位的坐標(biāo)信息導(dǎo)入MAPGIS軟件,繪制采樣點(diǎn)分布圖為后期數(shù)據(jù)分析提供依據(jù),并將硝酸根濃度數(shù)據(jù)導(dǎo)入ARCGIS軟件進(jìn)行分析,繪制研究區(qū)內(nèi)硝酸鹽空間分布圖。在進(jìn)行硝酸鹽污染分析過程中使用Origin軟件繪制箱型圖,進(jìn)一步對硝酸鹽污染特征進(jìn)行分析。此外,通過IBM SPSS Statistics軟件對各指標(biāo)與硝酸根的相關(guān)性進(jìn)行分析。

    表2 硝酸鹽、亞硝酸鹽和總硬度在地下水和飲用水標(biāo)準(zhǔn)中的標(biāo)準(zhǔn)值(mg·L-1)Table 2 Values of nitrogen and total hardness in groundwater and drinking water standards(mg·L-1)

    2 結(jié)果與討論

    2.1水質(zhì)基本狀況

    根據(jù)pH值、總硬度、硫酸鹽、氯離子、硝酸鹽、亞硝酸鹽和氨氮7項(xiàng)指標(biāo)的測定結(jié)果,采用GB/T 14848—1993中單項(xiàng)組分評價方法對研究區(qū)的地下水質(zhì)量進(jìn)行評價。結(jié)果表明,115個地下水樣品中,水質(zhì)為Ⅱ和Ⅲ類的點(diǎn)位個數(shù)分別為9個和22個,共占全部點(diǎn)位的26.96%;74個點(diǎn)位的水質(zhì)為Ⅴ類,占全部點(diǎn)位的64.35%。由此可見,研究區(qū)地下水質(zhì)量總體較差,主要的污染指標(biāo)為硝酸鹽和總硬度。研究區(qū)四個地下水單元的水質(zhì)評價結(jié)果見表3。G1和G2的水質(zhì)較差,均以Ⅴ類水為主,兩個地下水單元Ⅴ類水點(diǎn)位數(shù)量分別占該地下水單元總點(diǎn)位數(shù)的90.9%與93.48%;沖洪積扇扇中部孔隙水單元的水質(zhì)相對較好,其中G3點(diǎn)位中Ⅲ類水居多,占該地下水單元淺層孔隙水總點(diǎn)位的44.44%,Ⅴ類水?dāng)?shù)量占該地下水單元淺層孔隙水總點(diǎn)位的30.56%;G4的水質(zhì)最好,多為Ⅱ類和Ⅲ類水,分別占該地下水單元總點(diǎn)位數(shù)的63.64%和36.36%。這表明地下水水位越深,受污染程度相對越小。

    表3 不同單元地下水水質(zhì)評價結(jié)果Table 3 Classification of water quality in different groundwater units

    此外,對所采地下水進(jìn)行的基本化學(xué)指標(biāo)測試(表4)結(jié)果表明:研究區(qū)淺層地下水的水溫為13.1~ 18.2℃,平均值為15.4℃;DO范圍0.03~13.76 mg·L-1,平均值為7.42 mg·L-1。研究區(qū)深層地下水水溫為12.7~17.7℃,平均值為15.17℃;DO范圍1.38~8.03 mg·L-1,平均值為3.99 mg·L-1。一般來說,在沒有受海水入侵的含水層,EC值可作為地下水污染程度的指標(biāo),未受污染的地下水EC值一般小于3×10-2S·m-1[15]。研究區(qū)淺層地下水EC范圍為(3.67~31)×10-2S·m-1,平均值為9.55×10-2S·m-1,而深層地下水EC范圍是(3.08~10.20)×10-2S·m-1,平均值為6.30×10-2S·m-1。由此可見,研究區(qū)內(nèi)不論深層地下水還是淺層地下水均受到了嚴(yán)重的污染,而相對于淺層地下水來說,由于深層地下水采樣井深度較大,受污染程度相對較輕。淺層地下水Eh的范圍為-87.0~253 mV,平均值為138.8 mV,深層地下水Eh范圍是-36.4~241.0 mV,平均值為154.3 mV,整個研究區(qū)的地下水都處于弱還原性環(huán)境。

    表4 研究區(qū)內(nèi)地下水基本化學(xué)參數(shù)測試值Table 4 Basic chemical properties of groundwater in researched area

    2.2研究區(qū)硝酸鹽的總體含量特征

    由測試結(jié)果可知采樣點(diǎn)地下水硝酸鹽的濃度范圍為0.24~440 mg·L-1,平均濃度為63.09 mg·L-1,是GB/T 14848—1993中硝酸鹽Ⅲ類標(biāo)準(zhǔn)值的3倍以上,也遠(yuǎn)高于GB 5749—2006中硝酸鹽的標(biāo)準(zhǔn)限值要求。115個點(diǎn)位中,硝酸鹽濃度超過GB/T 14848—1993中Ⅲ類標(biāo)準(zhǔn)值要求的有82個,占全部點(diǎn)位的71.30%。由此可見,研究區(qū)地下水硝酸鹽污染非常嚴(yán)重。在所有點(diǎn)位中,有5個點(diǎn)位氨氮有檢出(檢出限為0.04 mg·L-1),其中位于G2的A29和A58氨氮濃度分別為3.00 mg·L-1和2.65 mg·L-1,按照GB/T 14848—1993判定為Ⅴ類水,其他點(diǎn)位能夠滿足Ⅲ類水標(biāo)準(zhǔn)。有42個點(diǎn)位亞硝酸鹽有檢出(檢出限為0.002 mg·L-1),檢出濃度范圍為0.002~33.20 mg·L-1,其中16個點(diǎn)位的亞硝酸鹽濃度不能滿足Ⅲ類水標(biāo)準(zhǔn)值要求,亞硝酸鹽檢出率高且濃度較高的點(diǎn)位多位于G1。依據(jù)各點(diǎn)位的硝酸鹽含量,通過插值法得到硝酸鹽在整個研究區(qū)的濃度分布情況(圖2),顯示G1和G2的南部區(qū)域硝酸鹽濃度整體較高??傮w而言,硝酸鹽為研究區(qū)地下水中主要的氮存在形式,是主要的污染指標(biāo)。

    圖2 研究區(qū)內(nèi)硝酸鹽空間分布特征Figure 2 Spatial distribution characteristics of nitrate in researched area

    2.3不同地下水單元硝酸鹽含量水平及影響因素

    各地下水單元硝酸鹽濃度的箱型圖見圖3。四個地下水單元中硝酸鹽的總體濃度水平由高到低依次為G1、G2、G3和G4。其中,G1硝酸鹽的濃度范圍為21.68~440 mg·L-1,平均濃度為128.17 mg·L-1,標(biāo)準(zhǔn)差為119.94 mg·L-1;G2硝酸鹽的濃度范圍為8.56~161.8 mg·L-1,平均濃度為73.01 mg·L-1,標(biāo)準(zhǔn)差為34.33 mg·L-1;G3硝酸鹽的濃度范圍為2.04~119 mg· L-1,平均濃度為28.62 mg·L-1,標(biāo)準(zhǔn)差為26.3 mg·L-1;G4硝酸鹽的濃度范圍為0.24~10.6 mg·L-1,平均濃度為4.2 mg·L-1,標(biāo)準(zhǔn)差為2.72 mg·L-1。研究區(qū)硝酸鹽濃度整體上呈現(xiàn)由西北到東南逐漸降低的趨勢。此外,G3、G4同處滹沱河沖洪積扇扇中部孔隙水單元,而G3淺層孔隙水中硝酸鹽濃度與G4深層孔隙水相比較高,此現(xiàn)象表明,地下水埋深越淺,硝態(tài)氮含量越高,超標(biāo)率和嚴(yán)重超標(biāo)率越高。Lotse等[16]的研究也得到了類似結(jié)果,這也意味著,地下水硝態(tài)氮很可能來源于上部土層的淋溶作用,淺層地下水最先受到污染,然后逐漸向深層滲漏。

    圖3 研究區(qū)內(nèi)不同地下水單元硝酸鹽含量水平Figure 3 Concentrations of nitrate in different groundwater units

    為識別影響硝酸鹽在不同地下水單元中分布的因素,對地下水單元含水層介質(zhì)的巖性及地下水埋深進(jìn)行了分析。G1含水層以各種片麻巖及大理巖裂隙、溶蝕裂隙和第四系松散沖積物為主,含水層厚度為10~20 m,水位埋藏較淺。G2含水層巖性為礫卵石和砂礫卵石,導(dǎo)水性和富水性好,滲透系數(shù)一般為300~400 m·d-1,防污性能較差,含水層厚度為20~35 m。G3 和G4含水層巖性顆粒變細(xì),粘土夾層增多,滲透系數(shù)為100~200 m·d-1,與G1和G2相比,滲透性降低,對垂向水流向下遷移的水中硝酸根離子起到了阻滯作用,同時這樣的環(huán)境(含水量高、氧化還原電勢低、高氮量)十分有利于反硝化作用的進(jìn)行而使氮元素?fù)p失。G3和G4的地下水巖性相同,但淺層孔隙水的含水層厚度為25~60 m,深層孔隙水的含水層厚度則為100~200 m,這使得深層孔隙水更不易被污染[17]。由此可見,含水層巖性及地下水埋深對硝酸鹽在不同地下水單元中的分布均具有重要影響。

    根據(jù)現(xiàn)場污染源調(diào)查情況看,在研究區(qū)西北部兩個地下水單元區(qū)域內(nèi)污水處理廠數(shù)量較少,而東南地區(qū)地下水單元區(qū)域內(nèi)污水處理廠數(shù)量相對較多。此外,畜禽養(yǎng)殖場在西北的兩個地下水單元中分布較多,而在東南地下水單元區(qū)域內(nèi)分布較少,基本符合研究區(qū)西北部污水處理廠較少且畜禽養(yǎng)殖場較多,而東南部污水處理廠相對較多而畜禽養(yǎng)殖場較少的趨勢。這一現(xiàn)象與硝酸鹽濃度在研究區(qū)內(nèi)的變化趨勢呈一定相關(guān)性,初步認(rèn)為除水文地質(zhì)條件的差異外,污水處理廠分布的情況及畜禽養(yǎng)殖場的數(shù)量是引起硝酸鹽濃度分布差異性的重要原因[18]。

    2.4地下水中硝酸鹽的來源分析

    流域氮素流失量與降雨徑流過程、土壤性質(zhì)、植被、管理措施、農(nóng)業(yè)生產(chǎn)方式(如化肥類型、施肥方式)等因素密切相關(guān)[17],硝酸鹽污染主要是由于氮的天然化合物溶解度較大,同時,人類的活動也起主要影響作用。人類活動會改變地下水的地球化學(xué)環(huán)境并相應(yīng)地改變其中化學(xué)元素的特性,如大量抽取地下水不但使得包氣帶加厚,而且改變了地下水的天然循環(huán)方式,不同程度地影響地下水的再循環(huán),而且人類活動排放的各種廢棄物也在一定程度上改變了地下水原來的氧化還原環(huán)境。這兩者都不同程度地改變了地下水的水文地球化學(xué)環(huán)境[19]。

    硝酸鹽是自然存在的離子,是氮循環(huán)的一部分,自然界存在的硝酸鹽在地下水中通常很少,一般來自人類活動,如污水灌溉、氮肥過度施用、人畜糞便及垃圾堆放等。在氧化性很強(qiáng)的地下水中,硝酸鹽會呈穩(wěn)定的溶解態(tài)氮形式,且會隨著地下水運(yùn)動遷移[13]。Rahmati等[20]對伊朗西部庫爾德斯坦的地下水水源地硝酸鹽污染程度進(jìn)行的研究表明,所有被污染的井水都位于農(nóng)業(yè)灌溉區(qū),且向灌溉區(qū)施加氮肥的量與當(dāng)?shù)氐叵滤邢跛猁}濃度密切相關(guān)。Ka?arogˇlu等[21]研究發(fā)現(xiàn),土耳其埃斯基謝希爾地區(qū)將城市污水處理廠出水以井灌形式進(jìn)行地下水回灌,從而導(dǎo)致回灌區(qū)地下水產(chǎn)生硝酸鹽污染問題。除了人類活動所帶來的氮污染外,硝酸鹽在地下水中的遷移與水文地質(zhì)條件也有很大的關(guān)系。Arauzo等[22]通過主成分分析表明,地下水硝酸鹽污染的主要原因是地下水的固有脆弱性(由水文地質(zhì)條件與土壤介質(zhì)共同決定),土地利用方式以及地下水流向都會對硝酸鹽的遷移產(chǎn)生影響。Chen等[23]通過對黃河三角洲地下水硝酸鹽的濃度進(jìn)行檢測分析得出,該研究區(qū)硝酸鹽污染主要發(fā)生在淺層地下水且濃度與該地區(qū)的地貌及土地利用方式呈空間相關(guān)性。

    一般來說,由城市污水和垃圾滲漏造成的污染,其硝酸根和鹽污染有較好的一致性。鹽污染包括溶解性總固體、硬度、硫酸鹽、氯化物等。而農(nóng)業(yè)施肥(不包括農(nóng)家肥)所帶來的硝酸鹽污染,水中硝酸根與鹽濃度則不一定呈較好的相關(guān)性。在環(huán)境學(xué)研究中,氯離子和硝酸根離子是兩種十分重要的離子。硝酸根離子由于其本身在環(huán)境中的行為而受到環(huán)境化學(xué)家的重視,而氯離子除了污染環(huán)境、導(dǎo)致土壤鹽漬化之外,由于其在地下水運(yùn)動過程中不受物理、化學(xué)、微生物的影響,還經(jīng)常作為指示地表水和地下水混合和反硝化作用的指示劑[24]。此外,總硬度也是生活飲用水水質(zhì)的一個重要指標(biāo),是由各種可溶性多價金屬離子形成的,主要是指Ca2+和Mg2+。地下水的總硬度主要來自地層,工業(yè)廢水和污染物的滲入也可導(dǎo)致總硬度升高。分析研究區(qū)內(nèi)硝酸根離子和氯離子及總硬度濃度的相互關(guān)系可以在一定程度上判斷出硝酸根離子的來源,本研究分析了不同地下水單元與硝酸鹽顯著相關(guān)的指標(biāo)及相關(guān)系數(shù),分析結(jié)果如表5所示。G1中,EC、Ca2+、Cl-、總硬度和溶解性總固體量值與硝酸鹽氮的濃度的相關(guān)系數(shù)分別為0.938、0.910、0.905、0.906、0.903,表明相關(guān)性較好,與G3和G4相比,G2各指標(biāo)與硝酸根相關(guān)系數(shù)較高,而G4鹽離子與硝酸根濃度的相關(guān)性并不顯著。這一現(xiàn)象也表明研究區(qū)內(nèi)G1的硝酸鹽污染可能來源于城市污水與垃圾滲濾液,而G3和G4的硝酸鹽則可能存在不同的來源。

    針對氯離子進(jìn)行分析可知,與G3和G4相比,在G1和G2中,硝酸鹽和氯離子相關(guān)性較好;G3中硝酸鹽氮與氯離子相關(guān)系數(shù)較低,相關(guān)性較差;而G4中雖然硝酸鹽氮與氯離子相關(guān)系數(shù)相對較高,但此結(jié)果是由于點(diǎn)數(shù)較少導(dǎo)致的,總的來說此單元中這兩個指標(biāo)相關(guān)性也較差。針對總硬度指標(biāo)進(jìn)行分析可知,硝酸鹽氮與總硬度相關(guān)性較好,但在G3和G4中相關(guān)系數(shù)則較低。一般認(rèn)為城市地下水硬度升高主要是由于污水滲漏造成的,其形成的化學(xué)機(jī)理是在污水入滲過程中鈣鎂易溶鹽和難溶鹽的溶解以及陽離子的交換作用,污水中有機(jī)物的分解是導(dǎo)致難溶鹽巖溶解的動力[25]。通過相關(guān)性分析可知,氯離子及總硬度和硝酸鹽氮相關(guān)性的好壞基本保持一致,都是西北兩個地下水單元的相關(guān)性較高,而東南地下水單元的相關(guān)性較差。由此可以初步推測,西北部兩個地下水單元的硝酸鹽污染主要來源于生活污水及畜禽養(yǎng)殖廢水,屬于點(diǎn)源污染;而東南滹沱河沖洪積扇地下水單元中硝酸鹽濃度較低且與氯離子及總硬度相關(guān)性較低,初步判斷該區(qū)硝酸鹽污染主要來源于農(nóng)業(yè)施肥,為面源污染。

    表5 不同地下水單元中各地下水化學(xué)參數(shù)與硝酸鹽的相關(guān)性分析Table 5 Correlation coefficients between nitrate and chemical properties in different groundwater units

    3 結(jié)論

    (1)研究區(qū)內(nèi)水質(zhì)為Ⅱ和Ⅲ類的點(diǎn)位占全部點(diǎn)位的26.96%;水質(zhì)為Ⅴ類的點(diǎn)位占全部點(diǎn)位的64.35%,地下水質(zhì)量總體較差,主要的污染指標(biāo)為硝酸鹽和總硬度。其中硝酸鹽的平均濃度與最高濃度分別為GB/T 14848—1993中硝酸鹽Ⅲ類標(biāo)準(zhǔn)值的3倍與22倍,屬于嚴(yán)重污染。

    (2)四個地下水單元中硝酸鹽濃度差異較大,硝酸鹽的平均濃度基本符合從西北到東南、從淺層到深層逐漸降低的趨勢。除當(dāng)?shù)氐乃牡刭|(zhì)條件外,污水處理廠及畜禽養(yǎng)殖場的分布也是影響硝酸鹽濃度的主要因素之一。

    (3)研究區(qū)西北部兩個地下水單元中硝酸鹽主要來源于點(diǎn)源類型的生活污水及畜禽養(yǎng)殖廢水,東南部兩個地下水單元的硝酸鹽主要來源于面源類型的農(nóng)業(yè)施肥。

    參考文獻(xiàn):

    [1]吳愛民,李長青,徐彥澤,等.華北平原地下水可持續(xù)利用的主要問題及對策建議[J].南水北調(diào)與水利科技, 2010, 8(6):110-113, 128. WU Ai-min, LI Chang-qing, XU Yan-ze, et al. Key issues influencing sustainable groundwater utilization and its countermeasures in North China Plain[J]. South-to-North Water Transfers and Water Science & Technology, 2010, 8(6):110-113, 128.

    [2]張兆吉,費(fèi)宇紅,郭春艷,等.華北平原區(qū)域地下水污染評價[J].吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 2012, 42(5):1456-1461. ZHANG Zhao-ji, FEI Yu-hong, GUO Chun-yan, et al. Reginoal groundwater contamination assessment in the North China Plain[J]. Journal of Jilin University(Earth Science Edition), 2012, 42(5):1456-1461.

    [3]李亞松,張兆吉,費(fèi)宇紅,等.河北省滹沱河沖積平原地下水質(zhì)量及污染特征研究[J].地球?qū)W報(bào), 2014, 35(2):169-176. LI Ya-song, ZHANG Zhao-ji, FEI Yu-hong, et al. Groundwater quality and contamination characteristics in the Hutuo River Plain area, Hebei Province[J]. Acta Geoscientica Sinica, 2014, 35(2):169-176.

    [4]衛(wèi)文,陳宗宇,趙紅梅,等.河北平原第四系承壓水4He與14C測年對比[J].吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 2011, 41(4):1144-1150.WEI Wen, CHEN Zong-yu, ZHAO Hong-mei, et al. Comparison of4He and14C dating of groundwater from quaternary confined aquifers in Hebei Plain[J]. Journal of Jilin University(Earth Science Edition), 2011, 41(4):1144-1150.

    [5]趙紅梅,趙華,毛洪亮,等.華北平原滹沱河沖洪積扇第四紀(jì)地層劃分[J].地層學(xué)雜志, 2014, 38(2):137-146. ZHAO Hong-mei, ZHAO Hua, MAO Hong-liang, et al. Quaternary stratigraphy division of Hutouhe alluvial fan deposits in the North China Plain[J]. Journal of Stratigraphy, 2014, 38(2):137-146.

    [6]劉少玉,靳盛海,韓雙平,等.華北東部平原淺層高礦化弱滲透地下水的開發(fā)技術(shù)示范研究[J].吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 2010, 40 (1):114-120. LIU Shao-yu, JIN Sheng-hai, HAN Shuang-ping, et al. Development and utilization mode of high salinity shallow groundwater from Low permeability aquifer in Eastern North China Plain[J]. Journal of Jilin University(Earth Science Edition), 2010, 40(1):114-120.

    [7]Fenech C, Rock L, Nolan K, et al. The potential for a suite of isotope and chemical markers to differentiate sources of nitrate contamination:A review[J]. Water Research, 2012, 46:2023-2041.

    [8]畢二平,李政紅.石家莊市地下水中氮污染分析[J].水文地質(zhì)工程地質(zhì), 2001, 9(2):31-34. BI Er-ping, LI Zheng-hong. Analysis of nitrates pollution in groundwater in Shijiazhuang City[J]. Hydrogeology & Engineering Geology, 2001, 9(2):31-34.

    [9]王翠華,孟素花,費(fèi)宇紅.滹沱河平原地下水礦化度和硬度變化分析[J].海河水利, 2009, 3(6):30-32. WANG Cui-hua, MENG Su-hua, FEI Yu-hong. Analysis on the change of ground water salinity and hardness in Hutuohe Plain[J]. Haihe Water Resources, 2009, 3(6):30-32.

    [10]葉浩,劉長禮,姜建梅,等.滹沱河石家莊段地下水污染風(fēng)險評價[J].地質(zhì)通報(bào), 2008, 27(7):1065-1070. YE Hao, LIU Chang-li, JIANG Jian-mei, et al. Risk assessments of groundwater pollution in the Shijiazhuang reach of the Hutuo River, China[J]. Geological Bulletin of China, 2008, 27(7):1065-1070.

    [11]王金哲,張光輝,聶振龍,等.滹沱河流域平原區(qū)人類活動對淺層地下水干擾程度量化研究[J].水土保持通報(bào), 2010, 30(2):65-69. WANG Jin-zhe, ZHANG Guang-hui, NIE Zhen-long, et al. Disturbing degree of mankind activities to groundwater in the Hutuo River Valley area[J]. Bulletin of Soil and Water Conservation, 2010, 30(2):65-69.

    [12]王開然,郭芳,姜光輝,等.15N和18O在桂林巖溶水氮污染源示蹤中的應(yīng)用[J].中國環(huán)境科學(xué), 2014, 34(9):2223-2230. WANG Kai-ran, GUO Fang, JIANG Guang-hui, et al. Application of15N and18O to nitrogen pollution source in karst water in Eastern Guilin [J]. China Environmental Science, 2014, 34(9):2223-2230.

    [13]劉興權(quán),許晶玉,江麗華,等.山東省種植區(qū)地下水硝酸鹽污染空間變異及分布規(guī)律研究[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2010, 29(6):1172-1179. LIU Xing-quan, XU Jing-yu, JIANG Li-hua, et al. Spatial variability and distribution pattern of groundwater nitrate pollution in farming regions of Shandong Province, China[J]. Journal of Agro-Environment Science, 2010, 29(6):1172-1179.

    [14]史入宇,崔亞莉,趙婕,等.滹沱河地區(qū)地下水適宜水位研究[J].水文地質(zhì)工程地質(zhì), 2013, 40(2):36-41. SHI Ru-yu, CUI Ya-li, ZHAO Jie, et al. A study of the suitable groundwater level of the Hutuo River area[J]. Hydrogeolgy & Eegineering Geology, 2013, 40(2):36-41.

    [15]劉中業(yè),徐建國,祁曉凡,等.地下水電導(dǎo)率與礦化度相關(guān)關(guān)系分析——以魯北平原為例[J].山東國土資源, 2013, 29(9):57-60, 64. LIU Zhong-ye, XU Jian-guo, QI Xiao-fan, et al. Analysis on groundwater conductivity and salinity correlation:Setting Lubei Plain in Shandong Province as an example[J]. Land and Resources in Shandong Province, 2013, 29(9):57-60, 64.

    [16]Lotse E G, Simmons K E, Baker D E. A field study of macropore flow under saturated conditions using a bromide tracer[J]. Journal of Soil and Water Conservation, 1991, 46(5):376-380.

    [17]陳丹,張冰,曾逸凡,等.基于SWAT模型的青山湖流域氮污染時空分布特征研究[J].中國環(huán)境科學(xué), 2015, 35(4):1216-1222. CHEN Dan, ZHANG Bing, ZENG Yi-fan, et al. Analysis of temporal and spatial characteristics of nitrogen pollution in Qingshan Lake basin based on SWAT model[J]. China Environmental Science, 2015, 35(4):1216-1222.

    [18]李亞松,張兆吉,費(fèi)宇紅,等.河北省滹沱河沖積平原地下水質(zhì)量及污染特征研究[J].地球?qū)W報(bào), 2014, 35(5):169-176. LI Ya-song, ZHANG Zhao-ji, FEI Yu-hong, et al. Groundwater quality and contamination characteristics in the Hutuo River Plain Area, Hebei Province[J]. Acta Geoscientica Sinica, 2014, 35(5):169-176.

    [19]劉君,陳宗宇.利用穩(wěn)定同位素追蹤石家莊市地下水中硝酸鹽來源[J].環(huán)境科學(xué), 2009, 30(6):1602-1607. LIU Jun, CHEN Zong-yu. Using stable isotope to trace the sources of nitrate in groundwater in Shijiazhuang[J]. Environmental Science, 2009, 30(6):1602-1607.

    [20]Rahmati O, Samani A N, Mahmoodi N, et al. Assessment of the contribution of N-fertilizers to nitrate pollution of groundwater in western Iran(case study:Ghorveh-Dehgelan aquifer)[J]. Water Qual Expo Health, 2015, 7:143-151.

    [21]Ka?arogˇluF, Günay G. Groundwater nitrate pollution in an alluvium aquifer, Eskis Eehir rurban area and its vicinity, Turkey[J]. Environmental Geology, 1997, 31(3/4):178-184.

    [22]Arauzo M, Martínez-Bastida J J. Environmental factors affecting diffuse nitrate pollution in the major aquifers of central Spain:Groundwater vulnerability vs. groundwater pollution[J]. Environ Earth Sci, 2015, 73(12):8271-8286.

    [23]Chen J Y, Taniguchi M, Liu G Q, et al. Nitrate pollution of groundwater in the Yellow River delta, China[J]. Hydrogeology Journal, 2007, 15 (8):1605-1614.

    [24]Liu C Q, Li S L, Lang Y C, et al. Using δ15N- and δ18O-values to identify nitrate sources in karst groundwater, Guiyang, Southwest China [J]. Environmental Science & Technology, 2006, 40(22):6928-6933.

    [25]王欣寶,谷明旭,王新友.石家莊市城區(qū)地下水硬度升高機(jī)理分析[J].安徽農(nóng)業(yè)科學(xué), 2013, 41(35):13709-13711. WANG Xin-bao, GU Ming-xu, WANG Xin-you. Mechanism of increasing groundwater hardness in urban district of Shijiazhuang[J]. Journal of Anhui Agri Sci, 2013, 41(35):13709-13711.

    中圖分類號:X522

    文獻(xiàn)標(biāo)志碼:A

    文章編號:1672-2043(2016)05-0947-08

    doi:10.11654/jaes.2016.05.019

    收稿日期:2015-11-05

    基金項(xiàng)目:環(huán)保公益性行業(yè)科研專項(xiàng)(201409029)

    作者簡介:劉琰(1978—),博士,研究員,主要從事流域水環(huán)境質(zhì)量評價及標(biāo)準(zhǔn)研究。E-mail:liuyan@craes.org.cn

    *通信作者:鄭丙輝E-mail:zhengbh@craes.org.cn

    Spatial distribution and influencing factors of nitrate content in groundwater of alluvial-pluvial fan of Hutuo River

    LIU Yan1, QIAO Xiao-cui1, JIANG Qiu-feng1, CHANG Sheng1, ZHU Yuan-yuan1,2, WANG Shan-jun1, ZHENG Bing-hui1*
    (1.State Key Laboratory of Environmental Criteria and Risk Assessment, State Environmental Protection Key Laboratory of Drinking Water Source Protection, Chinese Research Academy of Environmental Sciences, Beijing 100012, China; 2.Department of Environmental Engineering, School of Civil and Environmen tal Engineering, University of Science and Technology Beijing, Beijing 100083, China)

    Abstract:In order to identify the contamination profiles and spatial distribution of nitrate in the North China Plain, 115 groundwater samples,collected from Hutuo River pluvial fan in Shijiazhuang City, Hebei Province, were analyzed according to the parameters listed in the Quality Standard for Ground Water(GB 14848—1993). Results indicated that the water quality was generally poor. The 64.35%of total samples only met the Class V standards. Nitrate and hardness were the main failed indicators. The sampling sites with poor water quality were mainly located in the two groundwater units in the northwest area. The nitrate concentrations ranged from 0.24 to 440 mg·L-1, with the maximum and the average concentrations being 22 and 3 times higher than the ClassⅢstandard value, respectively. The nitrate concentrations spatially decreased from the north-west to the south-east and from shallow to deep stratum of the research area. Hydrogeological conditions had an important effect on the distribution of nitrate in groundwater. In addition, the type and distribution of pollution sources also had great contributions. Correlation analysis between chloridion and nitrate, and between total hardness and nitrate indicated that nitrate in the two units of the north-west area came from the discharges of domestic wastewater and livestock and poultry breeding wastewater from the point sources, while nitrate in the other two units located in the south-east area was mainly resulted from the agricultural non-point sources.

    Keywords:nitrate; groundwater; pollution profile; spatial distribution; alluvial-pluvial fan of Hutuo River

    猜你喜歡
    洪積扇滹沱河點(diǎn)位
    1年期國債收益率跌至1%
    銀行家(2025年1期)2025-02-08 00:00:00
    西藏拉薩河流域中下游洪積扇植被的物種組成與多樣性特征
    機(jī)器人快速示教方法及示教點(diǎn)位姿變換的研究
    青海托拉海沖洪積扇區(qū)水文地質(zhì)及水質(zhì)分析
    滹沱河
    公民與法治(2020年6期)2020-05-30 12:43:52
    垂直面內(nèi)建立基線的特殊點(diǎn)位高程測量法
    山前洪積扇形成的學(xué)生分組實(shí)驗(yàn)及改進(jìn)措施
    基于模糊數(shù)學(xué)的水資源價值評估及應(yīng)用:以滹沱河流域河北段為例
    卡普斯朗-臺勒丘克河洪積扇二維地下水流數(shù)值模擬
    山西滹沱河山區(qū)濕地生態(tài)系統(tǒng)健康評價
    国产精品三级大全| 18禁裸乳无遮挡免费网站照片| 一个人看视频在线观看www免费| 老司机福利观看| 91麻豆av在线| 最新中文字幕久久久久| 国内揄拍国产精品人妻在线| 国产精品美女特级片免费视频播放器| 国产91精品成人一区二区三区| 亚洲精品成人久久久久久| 久久这里只有精品中国| 99久久精品热视频| 久久亚洲精品不卡| 日韩欧美一区二区三区在线观看| 亚洲精品在线观看二区| 人人妻人人澡欧美一区二区| 日韩欧美精品v在线| 亚洲成人中文字幕在线播放| 亚洲国产精品sss在线观看| 五月伊人婷婷丁香| 久久亚洲精品不卡| 我的女老师完整版在线观看| 国产伦精品一区二区三区视频9| 男女做爰动态图高潮gif福利片| 九九久久精品国产亚洲av麻豆| 国产精品免费一区二区三区在线| 人人妻人人澡欧美一区二区| 欧美bdsm另类| 午夜福利欧美成人| 欧美日本亚洲视频在线播放| 日日啪夜夜撸| 国产男靠女视频免费网站| av中文乱码字幕在线| 国产精品伦人一区二区| 露出奶头的视频| 久久人人爽人人爽人人片va| 琪琪午夜伦伦电影理论片6080| 精品免费久久久久久久清纯| 精品久久久久久久末码| 免费搜索国产男女视频| 亚洲最大成人av| 午夜亚洲福利在线播放| 69av精品久久久久久| 夜夜爽天天搞| 国内精品宾馆在线| 能在线免费观看的黄片| 在线看三级毛片| 亚洲国产色片| 日本一二三区视频观看| 一本一本综合久久| 亚洲av免费高清在线观看| 欧美色视频一区免费| 男人和女人高潮做爰伦理| 日韩高清综合在线| 在现免费观看毛片| 国产欧美日韩一区二区精品| 一进一出抽搐动态| xxxwww97欧美| 亚洲人成网站在线播| 国内久久婷婷六月综合欲色啪| 老师上课跳d突然被开到最大视频| 级片在线观看| 日韩大尺度精品在线看网址| 亚洲午夜理论影院| 欧美人与善性xxx| 国产亚洲精品久久久久久毛片| 色5月婷婷丁香| 国内久久婷婷六月综合欲色啪| 一本精品99久久精品77| av黄色大香蕉| 亚洲av成人精品一区久久| 日韩欧美在线乱码| 黄色配什么色好看| 日本黄色视频三级网站网址| 欧美+亚洲+日韩+国产| 神马国产精品三级电影在线观看| 久久欧美精品欧美久久欧美| 美女免费视频网站| 香蕉av资源在线| 亚洲国产欧美人成| 一夜夜www| 午夜免费激情av| 亚洲国产色片| 成熟少妇高潮喷水视频| 国产成人a区在线观看| 赤兔流量卡办理| 日韩在线高清观看一区二区三区 | 日本黄大片高清| 日本免费一区二区三区高清不卡| 婷婷色综合大香蕉| 99在线人妻在线中文字幕| 一夜夜www| 欧美色视频一区免费| 午夜福利成人在线免费观看| 淫妇啪啪啪对白视频| netflix在线观看网站| 三级男女做爰猛烈吃奶摸视频| 悠悠久久av| 在线观看av片永久免费下载| 美女大奶头视频| 天堂动漫精品| 国产一区二区三区在线臀色熟女| 日本 欧美在线| 1024手机看黄色片| 欧美中文日本在线观看视频| 日韩欧美国产在线观看| 在线国产一区二区在线| 又紧又爽又黄一区二区| 国产高潮美女av| 美女免费视频网站| 久久精品国产清高在天天线| 久久精品综合一区二区三区| 久久亚洲真实| 少妇人妻精品综合一区二区 | 久久热精品热| 九九久久精品国产亚洲av麻豆| 亚洲精品国产成人久久av| 欧美成人a在线观看| 日本一二三区视频观看| 成年女人毛片免费观看观看9| 亚洲av成人精品一区久久| 成人国产综合亚洲| 日韩欧美国产一区二区入口| 精品一区二区免费观看| 乱人视频在线观看| 欧美成人一区二区免费高清观看| 国产69精品久久久久777片| 久久久精品大字幕| 成人国产综合亚洲| av在线蜜桃| 午夜福利欧美成人| 十八禁网站免费在线| 亚洲狠狠婷婷综合久久图片| 少妇猛男粗大的猛烈进出视频 | 欧美最黄视频在线播放免费| eeuss影院久久| 久久这里只有精品中国| 国产伦人伦偷精品视频| 日本欧美国产在线视频| 一卡2卡三卡四卡精品乱码亚洲| av在线亚洲专区| 又紧又爽又黄一区二区| 午夜福利在线在线| 在线观看一区二区三区| 变态另类丝袜制服| 成人午夜高清在线视频| 午夜精品在线福利| 欧美成人免费av一区二区三区| 天堂av国产一区二区熟女人妻| 女同久久另类99精品国产91| 一边摸一边抽搐一进一小说| 乱系列少妇在线播放| 国产精品伦人一区二区| 中文字幕久久专区| 成人高潮视频无遮挡免费网站| 国产高潮美女av| 99热网站在线观看| 狠狠狠狠99中文字幕| 可以在线观看的亚洲视频| 中文字幕免费在线视频6| 国产av麻豆久久久久久久| 亚洲色图av天堂| av天堂在线播放| 国产私拍福利视频在线观看| 欧美一区二区亚洲| 别揉我奶头~嗯~啊~动态视频| 免费黄网站久久成人精品| 国产精品人妻久久久久久| 亚洲中文字幕日韩| 欧美日韩精品成人综合77777| 91久久精品国产一区二区三区| 国内精品宾馆在线| 亚洲不卡免费看| netflix在线观看网站| 国产精品一区二区性色av| 少妇被粗大猛烈的视频| 在线观看美女被高潮喷水网站| 国产人妻一区二区三区在| 国产色爽女视频免费观看| 日本黄色视频三级网站网址| 日本色播在线视频| 九九热线精品视视频播放| 国产亚洲av嫩草精品影院| 欧美极品一区二区三区四区| 久久久久久久久中文| 在线播放无遮挡| 日韩精品青青久久久久久| 长腿黑丝高跟| 日本一二三区视频观看| 亚洲久久久久久中文字幕| 尾随美女入室| 欧美三级亚洲精品| 精品一区二区免费观看| 啪啪无遮挡十八禁网站| 成人二区视频| 18+在线观看网站| 男人舔奶头视频| 欧美一区二区国产精品久久精品| 1024手机看黄色片| 午夜精品一区二区三区免费看| 国产视频一区二区在线看| 成人国产麻豆网| 一区二区三区四区激情视频 | 热99在线观看视频| 制服丝袜大香蕉在线| 日韩欧美精品免费久久| 级片在线观看| 日本成人三级电影网站| 深夜a级毛片| 69人妻影院| 国产精品嫩草影院av在线观看 | 又紧又爽又黄一区二区| 草草在线视频免费看| 成人特级av手机在线观看| 久99久视频精品免费| 亚洲三级黄色毛片| 国产爱豆传媒在线观看| 天堂动漫精品| 男人的好看免费观看在线视频| 一本精品99久久精品77| 在线a可以看的网站| 男女下面进入的视频免费午夜| 欧美一区二区亚洲| 成年免费大片在线观看| 国产精品电影一区二区三区| 22中文网久久字幕| 国产成人福利小说| 在线免费十八禁| 国产综合懂色| 亚洲经典国产精华液单| 欧美日韩精品成人综合77777| 久久久久久久午夜电影| 别揉我奶头~嗯~啊~动态视频| 久久久精品大字幕| 国产免费一级a男人的天堂| 国产一区二区激情短视频| 成人av在线播放网站| 真人做人爱边吃奶动态| 国产精品综合久久久久久久免费| 亚洲熟妇中文字幕五十中出| 国产精品久久视频播放| 精品人妻一区二区三区麻豆 | 最新在线观看一区二区三区| 99久久精品一区二区三区| 国产一区二区三区视频了| 99久久中文字幕三级久久日本| 日韩高清综合在线| 十八禁国产超污无遮挡网站| 白带黄色成豆腐渣| 精品一区二区免费观看| 最后的刺客免费高清国语| 99久久成人亚洲精品观看| 变态另类成人亚洲欧美熟女| 精品午夜福利视频在线观看一区| 免费av观看视频| 久久久久久九九精品二区国产| 少妇的逼好多水| 国产精品一区www在线观看 | 黄色女人牲交| 久久久久国内视频| 99久久成人亚洲精品观看| 欧美区成人在线视频| 国产免费av片在线观看野外av| 亚洲精品久久国产高清桃花| 女的被弄到高潮叫床怎么办 | 狠狠狠狠99中文字幕| 老司机福利观看| 国产高清不卡午夜福利| 网址你懂的国产日韩在线| 午夜福利18| 亚洲av一区综合| 99久国产av精品| 国产精品乱码一区二三区的特点| 美女免费视频网站| 丝袜美腿在线中文| 中文在线观看免费www的网站| 欧美日韩黄片免| 99在线人妻在线中文字幕| 国内精品宾馆在线| 日日摸夜夜添夜夜添小说| 国产乱人视频| 淫秽高清视频在线观看| 非洲黑人性xxxx精品又粗又长| 日韩在线高清观看一区二区三区 | 超碰av人人做人人爽久久| 婷婷丁香在线五月| 婷婷亚洲欧美| 久久久久久久久大av| 国内精品宾馆在线| 欧美潮喷喷水| 人人妻人人澡欧美一区二区| 少妇裸体淫交视频免费看高清| 国产精品伦人一区二区| av天堂中文字幕网| 我的女老师完整版在线观看| 国产精品美女特级片免费视频播放器| 国产精品人妻久久久久久| 少妇人妻一区二区三区视频| 免费人成视频x8x8入口观看| 日本在线视频免费播放| 熟女电影av网| 久久久久久国产a免费观看| 中出人妻视频一区二区| 免费看光身美女| 色综合亚洲欧美另类图片| 日韩亚洲欧美综合| 亚洲人成网站高清观看| 成人av一区二区三区在线看| 丰满乱子伦码专区| 国产v大片淫在线免费观看| 欧美不卡视频在线免费观看| 亚洲一级一片aⅴ在线观看| 别揉我奶头~嗯~啊~动态视频| 中文资源天堂在线| 国产一区二区在线观看日韩| videossex国产| 日韩大尺度精品在线看网址| 亚洲乱码一区二区免费版| 男女视频在线观看网站免费| 老女人水多毛片| 久久久久久九九精品二区国产| 天堂动漫精品| 久久精品国产亚洲网站| 国产乱人伦免费视频| 精品久久国产蜜桃| 久久久久性生活片| 日本熟妇午夜| 日本精品一区二区三区蜜桃| 三级毛片av免费| h日本视频在线播放| 天堂动漫精品| 久久99热6这里只有精品| 国产乱人伦免费视频| 亚洲av中文av极速乱 | 天天一区二区日本电影三级| 一级a爱片免费观看的视频| 国产亚洲av嫩草精品影院| 亚洲天堂国产精品一区在线| 看十八女毛片水多多多| 最新中文字幕久久久久| 国产精品一区www在线观看 | 69av精品久久久久久| 在线播放国产精品三级| av在线蜜桃| 听说在线观看完整版免费高清| 99久国产av精品| 国产一级毛片七仙女欲春2| 九九热线精品视视频播放| 婷婷丁香在线五月| 两性午夜刺激爽爽歪歪视频在线观看| 精品久久久久久久久亚洲 | 免费看光身美女| 国产精品国产三级国产av玫瑰| 97热精品久久久久久| 中文字幕久久专区| 成人毛片a级毛片在线播放| 亚洲av第一区精品v没综合| 日本黄色片子视频| 国内精品一区二区在线观看| 丰满人妻一区二区三区视频av| 能在线免费观看的黄片| 色5月婷婷丁香| 国产不卡一卡二| 亚洲aⅴ乱码一区二区在线播放| 天堂影院成人在线观看| 一区二区三区激情视频| 久久久久久国产a免费观看| 午夜福利成人在线免费观看| 欧美三级亚洲精品| 深爱激情五月婷婷| 精品人妻一区二区三区麻豆 | 免费一级毛片在线播放高清视频| 亚洲 国产 在线| 日韩欧美免费精品| 国产av一区在线观看免费| 午夜福利在线观看免费完整高清在 | 又粗又爽又猛毛片免费看| 深夜精品福利| 国产精品日韩av在线免费观看| 性欧美人与动物交配| 免费无遮挡裸体视频| 亚洲综合色惰| 国产中年淑女户外野战色| .国产精品久久| 国产老妇女一区| 女生性感内裤真人,穿戴方法视频| 午夜福利成人在线免费观看| 91久久精品电影网| 亚洲av熟女| 啦啦啦韩国在线观看视频| 午夜福利18| 男人和女人高潮做爰伦理| 不卡视频在线观看欧美| 亚洲国产欧洲综合997久久,| 性欧美人与动物交配| 亚洲中文字幕一区二区三区有码在线看| 日日干狠狠操夜夜爽| 大又大粗又爽又黄少妇毛片口| 一区福利在线观看| 精品人妻熟女av久视频| 婷婷亚洲欧美| 欧美性猛交黑人性爽| 九九爱精品视频在线观看| 亚洲 国产 在线| 日韩欧美精品免费久久| 久久人妻av系列| 日韩一本色道免费dvd| 最近中文字幕高清免费大全6 | 麻豆久久精品国产亚洲av| 久久久久久久精品吃奶| 亚洲av中文字字幕乱码综合| 丰满的人妻完整版| 亚洲成a人片在线一区二区| 国产精品1区2区在线观看.| av在线天堂中文字幕| 国产色婷婷99| www日本黄色视频网| 国产人妻一区二区三区在| 亚洲av二区三区四区| 制服丝袜大香蕉在线| 国产伦精品一区二区三区四那| 国产蜜桃级精品一区二区三区| 欧美色欧美亚洲另类二区| 在线a可以看的网站| 国产精品精品国产色婷婷| 天堂网av新在线| 一区福利在线观看| 无人区码免费观看不卡| 九九在线视频观看精品| 欧美日韩乱码在线| 免费一级毛片在线播放高清视频| avwww免费| 精品无人区乱码1区二区| 欧美区成人在线视频| 欧美成人性av电影在线观看| 天堂av国产一区二区熟女人妻| 他把我摸到了高潮在线观看| 最新在线观看一区二区三区| 999久久久精品免费观看国产| 成人午夜高清在线视频| 午夜日韩欧美国产| 长腿黑丝高跟| 尤物成人国产欧美一区二区三区| 国产男靠女视频免费网站| 日韩亚洲欧美综合| 九九在线视频观看精品| 欧美成人性av电影在线观看| 午夜视频国产福利| 我的女老师完整版在线观看| 男人狂女人下面高潮的视频| 12—13女人毛片做爰片一| 国模一区二区三区四区视频| aaaaa片日本免费| 大又大粗又爽又黄少妇毛片口| 男人舔女人下体高潮全视频| 午夜精品久久久久久毛片777| 国产精品一区二区性色av| 高清毛片免费观看视频网站| 国产精品一区二区三区四区免费观看 | 国内少妇人妻偷人精品xxx网站| 最后的刺客免费高清国语| 成年版毛片免费区| 天美传媒精品一区二区| 大又大粗又爽又黄少妇毛片口| 国产91精品成人一区二区三区| 人妻久久中文字幕网| 亚洲国产精品sss在线观看| 人妻少妇偷人精品九色| 999久久久精品免费观看国产| 一区二区三区免费毛片| 亚洲最大成人中文| 黄色视频,在线免费观看| 精品久久久久久久人妻蜜臀av| 久久精品影院6| 成人美女网站在线观看视频| 国产私拍福利视频在线观看| 九九爱精品视频在线观看| 黄片wwwwww| 麻豆精品久久久久久蜜桃| 日本免费一区二区三区高清不卡| 男人和女人高潮做爰伦理| 国产视频一区二区在线看| 一个人观看的视频www高清免费观看| 黄色配什么色好看| 国产精品精品国产色婷婷| 午夜视频国产福利| 亚洲国产欧美人成| 日韩大尺度精品在线看网址| 久久久久久久久久黄片| 观看免费一级毛片| 日本 av在线| 又爽又黄无遮挡网站| 永久网站在线| 久久九九热精品免费| 日本黄色视频三级网站网址| 亚洲va在线va天堂va国产| 国产精品亚洲一级av第二区| 欧美xxxx黑人xx丫x性爽| 国语自产精品视频在线第100页| 国产高清有码在线观看视频| 人妻制服诱惑在线中文字幕| 俄罗斯特黄特色一大片| 91在线观看av| 成人特级av手机在线观看| 天堂av国产一区二区熟女人妻| 男女边吃奶边做爰视频| 日韩一区二区视频免费看| 国产高清视频在线观看网站| 久久人妻av系列| 亚洲精品色激情综合| 黄色丝袜av网址大全| 日韩欧美一区二区三区在线观看| 老司机福利观看| 亚洲精品成人久久久久久| 精品99又大又爽又粗少妇毛片 | 国产精品一区www在线观看 | 国产成人影院久久av| 波野结衣二区三区在线| 久久6这里有精品| 国产白丝娇喘喷水9色精品| 免费人成视频x8x8入口观看| 看十八女毛片水多多多| 丝袜美腿在线中文| 精品福利观看| 精品久久久久久久末码| 亚洲经典国产精华液单| 男女视频在线观看网站免费| 精品福利观看| 毛片一级片免费看久久久久 | 黄色女人牲交| 又爽又黄a免费视频| 色综合亚洲欧美另类图片| 蜜桃久久精品国产亚洲av| 久久久成人免费电影| 欧美国产日韩亚洲一区| 国产伦精品一区二区三区视频9| 99久久成人亚洲精品观看| 久久婷婷人人爽人人干人人爱| 极品教师在线视频| 深夜a级毛片| 极品教师在线视频| 日本 av在线| 免费黄网站久久成人精品| 他把我摸到了高潮在线观看| 又爽又黄a免费视频| 免费不卡的大黄色大毛片视频在线观看 | 欧美黑人欧美精品刺激| 中文在线观看免费www的网站| 久久中文看片网| 国产av一区在线观看免费| 亚洲美女视频黄频| 蜜桃久久精品国产亚洲av| 国产免费av片在线观看野外av| 久久久久精品国产欧美久久久| 丰满人妻一区二区三区视频av| 日本一本二区三区精品| 日本精品一区二区三区蜜桃| 窝窝影院91人妻| or卡值多少钱| 欧美激情国产日韩精品一区| 观看美女的网站| 黄色一级大片看看| 天天躁日日操中文字幕| 国语自产精品视频在线第100页| 97热精品久久久久久| 国内精品美女久久久久久| 国产精华一区二区三区| av福利片在线观看| 免费搜索国产男女视频| 国产高清不卡午夜福利| 美女cb高潮喷水在线观看| 欧美一区二区精品小视频在线| xxxwww97欧美| 久久国产乱子免费精品| 亚洲图色成人| 亚洲18禁久久av| 亚洲三级黄色毛片| 久久久精品欧美日韩精品| 伊人久久精品亚洲午夜| 中亚洲国语对白在线视频| 国产老妇女一区| av女优亚洲男人天堂| 此物有八面人人有两片| 干丝袜人妻中文字幕| 日本 av在线| 久久久精品大字幕| 国产精品三级大全| 伦精品一区二区三区| 色精品久久人妻99蜜桃| 亚洲一区二区三区色噜噜| 在线免费观看不下载黄p国产 | 欧美3d第一页| 午夜a级毛片| 国产aⅴ精品一区二区三区波| 丰满人妻一区二区三区视频av| 成人欧美大片| 琪琪午夜伦伦电影理论片6080| 国产亚洲精品av在线| 欧美zozozo另类| 国产黄片美女视频| 亚洲最大成人手机在线| 国产高清视频在线播放一区| 热99re8久久精品国产| 久久国产乱子免费精品| 高清在线国产一区| 国产av不卡久久| 亚洲,欧美,日韩| 欧美zozozo另类| 熟女电影av网| 午夜激情欧美在线| 18禁黄网站禁片免费观看直播| 精品人妻一区二区三区麻豆 | 97人妻精品一区二区三区麻豆| 日本三级黄在线观看| 国产高清视频在线观看网站|