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

    薊縣盤山地區(qū)淺層地下水水質(zhì)分布特征及影響因素識(shí)別

    2016-02-06 02:59:51曹陽(yáng)李明明楊耀棟申月芳馬晗宇
    華北地質(zhì) 2016年4期
    關(guān)鍵詞:盤山薊縣豐水期

    曹陽(yáng),李明明,楊耀棟,申月芳,馬晗宇

    (1.天津市地質(zhì)礦產(chǎn)測(cè)試中心,天津 300191;2.國(guó)土資源部天津礦產(chǎn)資源監(jiān)督檢測(cè)中心,天津 300191;3.天津市地質(zhì)調(diào)查研究院,天津 300191)

    薊縣盤山地區(qū)淺層地下水水質(zhì)分布特征及影響因素識(shí)別

    曹陽(yáng)1,2,李明明3,楊耀棟1,2,申月芳1,2,馬晗宇1,2

    (1.天津市地質(zhì)礦產(chǎn)測(cè)試中心,天津 300191;2.國(guó)土資源部天津礦產(chǎn)資源監(jiān)督檢測(cè)中心,天津 300191;3.天津市地質(zhì)調(diào)查研究院,天津 300191)

    針對(duì)天津市薊縣盤山地區(qū)官莊鎮(zhèn)開展了豐水期(2013年7~9月)與枯水期(2014年3~4月)地下水水質(zhì)取樣與分析。利用多元統(tǒng)計(jì)法研究了該區(qū)地下水總硬度、、、Cl-、TDS、Ca2+、Mg2+、Na+和CODMn等9項(xiàng)指標(biāo)的季節(jié)變化。采用因子分析法對(duì)地下水水質(zhì)影響因素進(jìn)行了識(shí)別,并通過(guò)一元線性回歸分析識(shí)別地下水水質(zhì)因子與研究區(qū)人口密度、耕地比例、禽類養(yǎng)殖、農(nóng)家院數(shù)量和工業(yè)產(chǎn)值等之間的相關(guān)性,識(shí)別出地下水污染來(lái)源。結(jié)果顯示:豐水期水質(zhì)劣于枯水期水質(zhì),溶濾作用和生活污染是研究區(qū)水質(zhì)的主要影響因素,工業(yè)污染和商業(yè)服務(wù)均對(duì)該區(qū)地下水造成明顯影響。

    地下水水質(zhì);一元線性回歸;因子分析法;污染識(shí)別;薊縣盤山

    地下水質(zhì)的分布特征反映了地下水形成演化特征和水質(zhì)變化情況,是自然因素和人為活動(dòng)對(duì)地下水共同影響的直接表現(xiàn)。地下水由于其埋藏條件的特殊性,形成演化受含水層介質(zhì)特性、化學(xué)組分、水動(dòng)力條件以及人為因素等影響,是一個(gè)相對(duì)復(fù)雜的過(guò)程[1-2]。隨著計(jì)算機(jī)技術(shù)和統(tǒng)計(jì)學(xué)的發(fā)展,因子分析法、主成分分析法和聚類分析法被廣泛應(yīng)用于水質(zhì)成因分析、污染來(lái)源及影響因素分析研究中[3-5]。開展相關(guān)研究不僅對(duì)揭示地下水化學(xué)特征和水化學(xué)形成作用具有重要的學(xué)術(shù)意義,同時(shí)對(duì)地下水水資源保護(hù)和可持續(xù)開發(fā)利用具有重要的實(shí)際價(jià)值[6-9]。

    天津市薊縣盤山地區(qū),蘊(yùn)藏著豐富的飲用天然礦泉水資源,分布著多家礦泉水生產(chǎn)企業(yè)。該區(qū)淺層地下水賦存在第四系松散沉積層和花崗巖全風(fēng)化裂隙含水層中,含水層底板埋深普遍小于60m,是礦泉水主要補(bǔ)給來(lái)源之一。近年來(lái)人們一方面對(duì)礦泉水的需求量增加,另一方面對(duì)礦泉水水源地的人為干擾也逐漸增大。目前針對(duì)區(qū)內(nèi)礦泉水開發(fā)利用與水質(zhì)現(xiàn)狀,開展了多項(xiàng)專題研究工作[10-13],但研究集中在礦泉水成因及資源開發(fā)利用上,地下水水質(zhì)污染調(diào)查和成因分析也以定性研究為主,忽略了淺層地下水水質(zhì)季節(jié)性變化與人為活動(dòng)的干擾。本文通過(guò)系統(tǒng)地采集研究區(qū)不同季節(jié)淺層地下水水質(zhì)樣品,利用多元統(tǒng)計(jì)法和一元線性回歸分析,對(duì)薊縣盤山淺層地下水水質(zhì)時(shí)空變化進(jìn)行研究,定量地研究了研究區(qū)水質(zhì)影響因素,識(shí)別出淺層地下水污染來(lái)源,可為該區(qū)地下水水質(zhì)監(jiān)測(cè)管理與礦泉水污染控制提供依據(jù)。

    1 研究區(qū)概況

    天津市薊縣盤山地區(qū)位于薊縣西部,區(qū)內(nèi)自然條件優(yōu)越、資源物產(chǎn)豐富,盛產(chǎn)優(yōu)質(zhì)低鈉、低礦化度偏硅酸型礦泉水。工作區(qū)位于燕山山脈與華北平原交接地帶,地勢(shì)北高南低,主要地貌單元為侵蝕-剝蝕中低山丘陵和沖洪積扇(圖1)。北部屬盤山的一部分,侵蝕切割強(qiáng)烈,山勢(shì)陡峭,溝谷狹窄,多呈“V”字型。山區(qū)與平原交接過(guò)渡帶,發(fā)育坡殘積洪積裙,寬度約1~2km,第四系厚0~55m,花崗巖全風(fēng)化層平均厚度30m,巖性以砂質(zhì)粘土、亞砂土、粉細(xì)砂為主,夾卵礫石透鏡體,沖洪積裙上發(fā)育有一系列沖溝,寬度由幾米到幾十米不等,切割深度1~5m,沖溝兩側(cè)多為陡坎。研究區(qū)降水量主要集中在6~9月,達(dá)到全年降水總量的80%以上,多年平均降水量為627.2mm。該區(qū)淺層地下水水位動(dòng)態(tài)主要受降水影響,地下水豐水期一般出現(xiàn)在7~9月,枯水期出現(xiàn)在3~5月。

    圖1 研究區(qū)位置及采樣點(diǎn)分布Fig.1 Location of the research area and distribution of the sampling points

    2 采樣與分析方法

    2.1 采樣與分析

    本次研究對(duì)薊縣盤山地區(qū)淺層地下水進(jìn)行了豐水期和枯水期采樣與分析,時(shí)間分別為2013年7~9月和2014年3~4月?,F(xiàn)場(chǎng)取樣采用2 L塑料桶,取井、孔中水樣時(shí),先抽水10分鐘左右再進(jìn)行采樣,采樣前用去離子水清洗,待采水樣潤(rùn)洗3遍。采樣后樣品置于4℃采樣箱保存,24小時(shí)內(nèi)送國(guó)土資源部天津礦產(chǎn)資源監(jiān)督檢測(cè)中心進(jìn)行檢測(cè)。檢測(cè)的水質(zhì)指標(biāo)為總硬度、Cl-、TDS(溶解性總固體)、Ca2+、Mg2+、Na+和CODMn(高錳酸鹽指數(shù))共9項(xiàng)。其中總硬度、、Ca2+、Mg2+、Na+采用ICP等離子體發(fā)射光譜法進(jìn)行測(cè)定采用紫外分光光度法進(jìn)行測(cè)定,Cl-采用硝酸銀容量法進(jìn)行測(cè)定,TDS采用重量法進(jìn)行測(cè)定、CODMn采用高錳酸鉀法進(jìn)行測(cè)定。分析過(guò)程中進(jìn)行質(zhì)量控制,每一批次樣品進(jìn)行全流程空白分析2個(gè),每10個(gè)樣品進(jìn)行加標(biāo)回收率樣品測(cè)試2個(gè),每一批次樣品選擇10%的樣品進(jìn)行平行樣分析,保證了測(cè)試結(jié)果準(zhǔn)確可靠。

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

    采用單因子標(biāo)準(zhǔn)指數(shù)法評(píng)價(jià)水質(zhì)現(xiàn)狀,利用S= Ci/Cs計(jì)算標(biāo)準(zhǔn)指數(shù),其中Ci為調(diào)查點(diǎn)水質(zhì)單因子實(shí)測(cè)值,Cs為地下水質(zhì)量標(biāo)準(zhǔn)中三類水限值,單位均為mg/L。

    多元統(tǒng)計(jì)方法能夠從統(tǒng)計(jì)數(shù)據(jù)中分析各水質(zhì)點(diǎn)潛在相關(guān)信息,被廣泛應(yīng)用于水質(zhì)時(shí)空變異性及污染源識(shí)別的研究中。本次研究利用因子分析法將原始變量和指標(biāo)化為少數(shù)幾個(gè)代表原變量信息的綜合指標(biāo)因子,達(dá)到降維和提取主要影響因素的目的,采用KMO和Bartlett球形檢驗(yàn)數(shù)據(jù)運(yùn)用因子分析的適用性。統(tǒng)計(jì)分析過(guò)程主要采用spss17.0完成。

    利用線性回歸分析方法研究地下水水質(zhì)影響因素與區(qū)內(nèi)人口、土地利用情況和工業(yè)產(chǎn)值的相關(guān)性,識(shí)別研究區(qū)地下水水質(zhì)主要影響因素。

    3 結(jié)果與分析

    3.1 地下水水質(zhì)及季節(jié)性變化特征

    研究區(qū)淺層地下水取水段埋深均小于60m,是該區(qū)自用民井的主要開采層位,根據(jù)《中華人民共和國(guó)地下水質(zhì)量標(biāo)準(zhǔn)》(T14848-93)對(duì)區(qū)內(nèi)46個(gè)采樣點(diǎn)指標(biāo)進(jìn)行評(píng)價(jià),總硬度、NO3-N、TDS、CODMn四項(xiàng)指標(biāo)存在超過(guò)地下水Ⅲ類標(biāo)準(zhǔn)的情況,其中總硬度和NO3-N超標(biāo)比例比較高,分別為19.6%和32.6%。超標(biāo)項(xiàng)目中,利用地下水質(zhì)量標(biāo)準(zhǔn)對(duì)總硬度進(jìn)行評(píng)價(jià),符合Ⅳ類水標(biāo)準(zhǔn)的調(diào)查點(diǎn)占8.7%,符合Ⅴ類水標(biāo)準(zhǔn)的調(diào)查點(diǎn)占10.9%,達(dá)到Ⅲ類及以上地下水標(biāo)準(zhǔn)的調(diào)查點(diǎn)占比為80.4%。若對(duì)NO3-N利用地下水質(zhì)量標(biāo)準(zhǔn)進(jìn)行評(píng)價(jià),符合Ⅳ類水標(biāo)準(zhǔn)的調(diào)查點(diǎn)占8.7%,符合Ⅴ類水標(biāo)準(zhǔn)的調(diào)查點(diǎn)占23.9%,達(dá)到Ⅲ類及以上地下水標(biāo)準(zhǔn)的調(diào)查點(diǎn)占比為67.4%。TDS、CODMn超標(biāo)率均為3.4%,為地下水Ⅳ類水。

    分別對(duì)研究區(qū)豐水期29個(gè)淺層水水樣和枯水期17個(gè)淺層水水樣的9項(xiàng)指標(biāo)進(jìn)行了統(tǒng)計(jì)描述,9項(xiàng)指標(biāo)平均濃度值均為豐水期大于枯水期;其中總硬度、NO3-N、Cl-、TDS、Ca2+、Mg2+和Na+等7項(xiàng)指標(biāo)的最小值和最大值均出現(xiàn)在豐水期;的最小值與最大值出現(xiàn)在枯水期,CODMn最小值出現(xiàn)在枯水期,最大值出現(xiàn)在豐水期。根據(jù)地下水質(zhì)量標(biāo)準(zhǔn)(GB/T-14848-93)超標(biāo)的項(xiàng)目主要有總硬度、NO3-N、TDS和CODMn,其中總硬度在豐水期超標(biāo)率為20.7%,最大標(biāo)準(zhǔn)指數(shù)為2.1,枯水期超標(biāo)率為17.6%,最大標(biāo)準(zhǔn)指數(shù)為1.4 ;NO3-N豐水期超標(biāo)率為34.5%,最大標(biāo)準(zhǔn)指數(shù)為3.7,枯水期超標(biāo)率為29.4%,最大標(biāo)準(zhǔn)指數(shù)為3.6 ;TDS和CODMn均在豐水期單個(gè)點(diǎn)存在超標(biāo)現(xiàn)象,超標(biāo)率為3.4%,最大標(biāo)準(zhǔn)指數(shù)均為1.2 。水質(zhì)總體超標(biāo)情況為豐水期超標(biāo)率大于枯水期(圖2)。

    表1 研究區(qū)采樣點(diǎn)水質(zhì)數(shù)據(jù)結(jié)果Tab.1 The results of water quality data of sampling points in the study area

    3.2 地下水主要污染因子識(shí)別

    對(duì)調(diào)查點(diǎn)水質(zhì)的9項(xiàng)指標(biāo)進(jìn)行數(shù)據(jù)檢驗(yàn),得到KMO(Kaiser-Meyer-Olkin檢驗(yàn)統(tǒng)計(jì)量)值為0.7,Bartlett球形檢驗(yàn)顯著性水平sig為0.00,低于0.01 的顯著性水平,說(shuō)明各因素具有一定相關(guān)性,適合運(yùn)用因子分析法進(jìn)行分析。利用主成分法計(jì)算載荷矩陣,最終選取3個(gè)因子作為主要影響因素,累計(jì)方差貢獻(xiàn)率達(dá)到89.4%,可以解釋研究區(qū)地下水約90%的影響因素(表2)。

    經(jīng)分析提取影響研究區(qū)水質(zhì)的主因子F1,總方差貢獻(xiàn)率為57.57%,綜合主成分貢獻(xiàn)率為64.4%,是該區(qū)水質(zhì)的主要影響因子,因子載荷較高的水質(zhì)因子主要為總硬度、Ca2+、TDS、Mg2+、NO3-N和Cl-;第二主因子F2,總方差貢獻(xiàn)率為18.74%,綜合主成分貢獻(xiàn)率20.96%,因子載荷較高的水質(zhì)因子主要為、Na+;第三主因子F3,總方差貢獻(xiàn)率為13.09%,綜合主成分貢獻(xiàn)率為14.64%,載荷因子主要為CODMn。

    3.3 地下水主要污染源識(shí)別

    利用一元回歸分析對(duì)水質(zhì)因子及特征進(jìn)行分析,研究各水質(zhì)因子與各村莊的人口密度、耕地比例、畜禽養(yǎng)殖種類和農(nóng)家院數(shù)量以及工業(yè)產(chǎn)值之間的相關(guān)性,相關(guān)系數(shù)結(jié)果見(jiàn)表3。結(jié)果顯示研究區(qū)各水質(zhì)參數(shù)與人口密度相關(guān)性明顯,其中人口密度與NO3-N具有強(qiáng)烈正相關(guān),與TDS具有強(qiáng)烈負(fù)相關(guān),同時(shí)人口密度還影響到總硬度、CODMn、Cl-和Na+,影響因子眾多;畜禽養(yǎng)殖和農(nóng)家院對(duì)水質(zhì)各因子影響不明顯,畜禽養(yǎng)殖與Na+有一定的正相關(guān)性,農(nóng)家院數(shù)量與CODMn呈正相關(guān)性,與Cl-呈負(fù)相關(guān);工業(yè)產(chǎn)值與各因子相關(guān)性最大,其中影響最為明顯的是,與工業(yè)產(chǎn)值相關(guān)性達(dá)到0.76,其次為Na+和總硬度,相關(guān)系數(shù)分別為0.604 和0.516 。

    圖2 淺層地下水水質(zhì)季節(jié)分布特征Fig.2 Seasonal variation of the water quality paramters in the shallow groundwater

    表2 因子載荷矩陣Tab.2 Factor loadingmatrix

    表3 水質(zhì)參數(shù)與影響因素的相關(guān)系數(shù)Tab.3 The correlation coefficient of the water quality parameters and influencing factors

    4 討論

    研究區(qū)各離子季節(jié)性變化顯示,總硬度、NO3-N、Cl-、TDS、Ca2+、Mg2+和Na+最大值與最小值均出現(xiàn)在豐水期,該區(qū)豐水期降水量相對(duì)集中,隨著降水入滲補(bǔ)給量的增加,各離子垂向淋溶補(bǔ)給量有所增加,導(dǎo)致地下水中各離子含量上升,同時(shí)該7項(xiàng)離子在空間分布上存在一定的差異性,部分地區(qū)離子富集,隨著降水量增大濃度逐漸升高;局部地區(qū)上覆地層離子較貧乏,隨著降水量增大濃度逐漸減小,影響因素主要呈面狀分布。最大值和最小值均出現(xiàn)在枯水期,說(shuō)明該因素在區(qū)域上受到降水量影響較小,與地下水徑流關(guān)系不明顯,主要受到其他因素的影響;CODMn最小值出現(xiàn)在枯水期,最大值出現(xiàn)在豐水期,表明CODMn不僅受到季節(jié)性變化的影響,同時(shí)還受到面狀分布趨勢(shì)的影響,在枯水期地表污染物垂向補(bǔ)給量減小,CODMn出現(xiàn)最小值,豐水期污染物隨著雨水遷移到地下水中,導(dǎo)致地下水中CODMn升高。

    回歸分析顯示區(qū)內(nèi)總硬度、NO3-N、和Cl-與人口密度具有正相關(guān)性,NO3-N、與人口密度呈強(qiáng)相關(guān),說(shuō)明研究區(qū)內(nèi)各村落人口數(shù)量和人類活動(dòng)對(duì)地下水水質(zhì)造成十分明顯的影響,尤其是生活污染源氨氮的排放導(dǎo)致地下水中硝酸鹽含量顯著升高,從研究區(qū)相同村落調(diào)查點(diǎn)NO3-N分布規(guī)律來(lái)看(表4),村內(nèi)地下水NO3-N標(biāo)準(zhǔn)指數(shù)遠(yuǎn)遠(yuǎn)高于村落周邊。其次工業(yè)污染源對(duì)地下水水質(zhì)的影響較為明顯,通過(guò)2013—2014年對(duì)工作區(qū)31處地表水的取樣分析,根據(jù)《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB3838-2002)進(jìn)行單因子指標(biāo)評(píng)價(jià),發(fā)現(xiàn)當(dāng)?shù)氐乇硭傮w水質(zhì)狀況較差,達(dá)到或優(yōu)于地表水III類水標(biāo)準(zhǔn)的僅占25.8%,超過(guò)一半的地表水為劣V類水。

    表4 研究區(qū)同村鄰近調(diào)查井中硝酸鹽標(biāo)準(zhǔn)指數(shù)差異Tab.4 The differences of nitrate standard survey index in study area near the village wells

    當(dāng)?shù)刂饕绊懙乇硭|(zhì)量的指標(biāo)為總氮,僅此一項(xiàng)就使51.6%的地表水成為劣V類水。此外,CODMn、CODCr(化學(xué)需氧量)、氨氮和總磷也是主要的影響指標(biāo),本次研究中,4處食品廠附近的坑塘均是劣V類地表水,其總氮、TDS等指標(biāo)濃度較高,其中總氮最高可達(dá)23.05mg/L,坑塘底泥中總氮含量也偏高,屬于典型的食品工業(yè)廢水污染。當(dāng)?shù)刂饕称饭I(yè)排污坑污水中Cl-含量都超過(guò)1 000mg/L,而在地下水環(huán)境中,Cl-是保守離子,性質(zhì)穩(wěn)定,一般不參與水文地球化學(xué)反應(yīng)。西后子峪和北小屯淺層地下水中Cl-含量偏高可以作為地下水受到食品工業(yè)污水影響的直接證據(jù)。同時(shí),食品工業(yè)污水中同樣含有大量的Na+,這些Na+大大增加了離子交換作用的強(qiáng)度,使更多的Ca2+、Mg2+進(jìn)入地下水中,導(dǎo)致了硬度的升高。同時(shí),據(jù)前人研究,區(qū)內(nèi)部分地下水井水中的等組分受到區(qū)內(nèi)工業(yè)尤其是食品廠等企業(yè)排污影響,濃度存在不同程度的逐年增高的趨勢(shì)[8],從側(cè)面反映了工業(yè)對(duì)濃度變化的影響。雖然工作區(qū)淺層地下水中硝酸鹽主要來(lái)自于生活污染源,但溶解性總固體含量和離子組分異常證明了食品工業(yè)污水對(duì)當(dāng)?shù)氐牡叵滤|(zhì)同樣造成了污染;農(nóng)家院數(shù)量與CODMn呈正相關(guān)性,說(shuō)明該區(qū)農(nóng)家院對(duì)地下水中有機(jī)物含量和還原性無(wú)機(jī)物含量造成一定的影響。

    因此研究區(qū)主要的影響因素為以總硬度、Ca2+、TDS、Mg2+、NO3-N和Cl-為載荷的自然因子和生活污染因子,主要成因是研究區(qū)長(zhǎng)期的地下水溶濾作用和生活垃圾及污水排放。其次是以、Na+為主要荷載的工業(yè)污染,主要來(lái)源于工業(yè)廢水排放和地下水自然成因。最后以CODMn為載荷的商業(yè)服務(wù)因素也是造成研究區(qū)地下水水質(zhì)變差的主要因素之一。研究區(qū)內(nèi)耕地和畜禽養(yǎng)殖對(duì)地下水影響不明顯。

    5 結(jié)論

    (1)研究區(qū)總體水質(zhì)指標(biāo)滿足地下水Ⅲ類標(biāo)準(zhǔn),但水質(zhì)普遍較差,地下水中總硬度、NO3-N、TDS、CODMn四項(xiàng)指標(biāo)存在超標(biāo)情況,其中總硬度和NO3-N超標(biāo)比例分別為19.6%和32.6%,NO3-N最大超標(biāo)倍數(shù)為3.7 。研究區(qū)地下水水質(zhì)總體表現(xiàn)為豐水期水質(zhì)劣于枯水期水質(zhì),季節(jié)變化明顯。

    (2)將影響研究區(qū)地下水水質(zhì)的主要因素概化為以總硬度、Ca2+、TDS、Mg2+、NO3-N和Cl-為荷載的自然因子和生活污染因子,綜合主成分貢獻(xiàn)率為64.4%;以、Na+為載荷的工業(yè)污染因子F2,總方差貢獻(xiàn)率為18.74%;以CODMn為載荷的商業(yè)服務(wù)污染因子,總方差貢獻(xiàn)率為13.09%。

    (3)研究區(qū)地下水水質(zhì)總體為自然成因,同時(shí)受到人類活動(dòng)影響,其中總硬度和硝酸鹽受生活垃圾及污水排放影響,與人口密度密切相關(guān),工業(yè)污水排放以及商業(yè)服務(wù)行業(yè)的發(fā)展均對(duì)該區(qū)地下水造成明顯影響。

    [1]滕彥國(guó),左銳,王金生,等.區(qū)域地下水演化的地球化學(xué)研究進(jìn)展[J].水科學(xué)進(jìn)展,2010,21(1):127-136.

    [2]曹陽(yáng),滕彥國(guó),劉昀竺.寧夏吳忠市金積水源地地下水水質(zhì)影響因素的多元統(tǒng)計(jì)分析[J].吉林大學(xué)學(xué)報(bào)(地球科學(xué)版),2013,43(1):235-244.

    [3]Shrestha S,Kazama F.Assessment of Surface Water Quality Usingmultivariate Statistical Techniques:A Case Study of the Fuji River Basin,Japan[J].Environmentalmodelling and Software,2007,22 (4):464-475.

    [4]Zhou F,Huang G,Guo H,et al.Spatio-temporal Patterns and Source Apportionment of Coastal Water Pollution in Eastern Hong Kong[J].Water Research,2007,41 (15):3429-3439.

    [5]張妍,尚金城,于相毅.主成分-聚類復(fù)合模型在水環(huán)境管理中的應(yīng)用:以松花江吉林段為例[J].水科學(xué)進(jìn)展,2005,16(4):592-595.

    [6]章光新,鄧偉,何巖,等.中國(guó)東北松嫩平原地下水水化學(xué)特征與演變規(guī)律[J].水科學(xué)進(jìn)展,2006,17(1):20-28.

    [7]王文祥,王瑞久,李文鵬,等.塔里木盆地河水氫氧同位素與水化學(xué)特征分析[J].水文地質(zhì)工程地質(zhì),2013,40(4):29-35.

    [8]張婷,劉靜玲,王雪梅.白洋淀水質(zhì)時(shí)空變化及影響因子評(píng)價(jià)與分析[J].環(huán)境科學(xué)學(xué)報(bào),2010,30(2):261-267.

    [9]黃金良,黃亞玲,李青生,等.流域水質(zhì)時(shí)空分布特征及其影響因素初析[J].環(huán)境科學(xué),2012,33(4):1098-1107.

    [10]李學(xué)美,趙子軍,崔亞莉,等.盤山山前地下水與礦泉水污染成因分析[J].水文地質(zhì)工程地質(zhì),2008(1):112-116.

    [11]劉兆茂,沈鍵,郭晶.薊縣盤山花崗巖區(qū)天然礦泉水形成機(jī)理分析[J].勘察科學(xué)技術(shù),2004(6):49-53.

    [12]王作友,張偉,李希武.天津市薊縣盤山礦泉水形成機(jī)理、賦存規(guī)律及開發(fā)利用研究報(bào)告[R].天津:天津市薊縣水資源開發(fā)技術(shù)研究所,1993:1-10.

    [13]馬晗宇,楊耀棟,曹陽(yáng),等.天津市薊縣盤山地區(qū)礦泉水水源地水質(zhì)評(píng)價(jià)研究成果報(bào)告[R].天津:天津市地質(zhì)礦產(chǎn)測(cè)試中心,2015:5-12.

    Analysis of variation and its influence factors of the shallow groundwater quality in the Panshan area of Jixian,Tianjin

    CAO Yang1,2,LIming-ming3,YANG Yao-dong1,2,SHEN Yue-fang1,2,MAHan-yu1,2
    (1.Tianjin Geology andmineral Testing Center,Tianjin 300191 ,China;2.Tianjin Supervision and Inspection Center ofmineral Resources,Ministry of Land and Resources,Tianjin 300191 ,China;3.Tianjin Institute ofGeologicalSurvey,Tianjin 300191 ,China)

    Water samples are collected in flood season(July,August and September)and average season(March and April)in the Panshan area of Jixian,Tianjin city in 2013.Multivariate statistics analysis are used to study the seasonal variation of the nine indicators of groundwater including of total hardness,,,Cl-,TDS,Ca2+,Mg2+, Na+and CODMn.Influencing factors of groundwater quality variation are identified by factor analysis.The unitary linearity regressionmethod are used to explore the correlation of groundwater quality and the population density, proportion of cultivated land,the number of livestock breeding,the number of farmyards and industrial production.This article shows the source of groundwater pollution.The results show that water quality in flood season is worse than that in average season.Leaching effect and domestic pollution aremain influencing factors. Industrial pollution and commercial services have significant impact to the groundwater.

    groundwater quality,unitary linearity regressionmethod,factor analysis,identification of pollution source,the Panshan area of Jixian

    P641

    A

    1672-4135(2016)04-0305-06

    2016-04-29

    天津市礦產(chǎn)資源補(bǔ)償費(fèi)項(xiàng)目“天津市薊縣盤山地區(qū)礦泉水水源地水質(zhì)評(píng)價(jià)研究(國(guó)土房任[2013]15號(hào))”

    曹陽(yáng)(1984-),女,工程師,2013年畢業(yè)于北京師范大學(xué)地下水科學(xué)與工程專業(yè),主要從事水工環(huán)地質(zhì)研究工作,E-mail:caoy84@163.com。

    猜你喜歡
    盤山薊縣豐水期
    豐水期巢湖水體與沉積物重金屬污染及分配特征研究
    牛欄江流域浮游動(dòng)物群落結(jié)構(gòu)特征及驅(qū)動(dòng)因子
    盤山抗日根據(jù)地舊址
    求知(2023年10期)2023-10-26 01:00:12
    膠東沿海某地區(qū)生活飲用水總α、總β放射性水平調(diào)查分析
    磐安縣盤山中心成人(社區(qū))學(xué)校
    ——培訓(xùn)掠影
    磐安縣盤山中心成人(社區(qū))學(xué)校
    薊縣風(fēng)采
    游薊縣梨木臺(tái)
    ——紀(jì)念上山下鄉(xiāng)48周年
    中華魂(2017年3期)2017-11-22 02:01:55
    聚盤山
    黃河之聲(2016年24期)2016-02-03 09:01:52
    區(qū)域重力調(diào)查在薊縣-滄州地區(qū)的應(yīng)用及認(rèn)識(shí)
    美女主播在线视频| 国产在视频线精品| 大片免费播放器 马上看| 超碰av人人做人人爽久久| 精品一区二区免费观看| 午夜福利高清视频| 亚洲国产精品999| 国产乱人偷精品视频| 亚洲成人一二三区av| 亚洲欧美一区二区三区国产| 18禁在线无遮挡免费观看视频| 国产亚洲最大av| 男女国产视频网站| 舔av片在线| 少妇被粗大猛烈的视频| 亚洲性久久影院| 婷婷色av中文字幕| 亚洲精品成人av观看孕妇| 中文乱码字字幕精品一区二区三区| 极品教师在线视频| 黄色配什么色好看| 国产精品偷伦视频观看了| 精品久久久久久久久亚洲| 欧美日韩在线观看h| 国产午夜福利久久久久久| 男人狂女人下面高潮的视频| 国产国拍精品亚洲av在线观看| 在线天堂最新版资源| 欧美激情国产日韩精品一区| 岛国毛片在线播放| 国产成人精品福利久久| 成年女人在线观看亚洲视频 | 男女边吃奶边做爰视频| 国产成人精品久久久久久| 欧美老熟妇乱子伦牲交| 人妻少妇偷人精品九色| 精品久久久久久久久av| 纵有疾风起免费观看全集完整版| 久热这里只有精品99| 插阴视频在线观看视频| 成人亚洲精品一区在线观看 | a级毛色黄片| 99久久精品热视频| 国产毛片a区久久久久| 中文资源天堂在线| 亚洲av二区三区四区| 搞女人的毛片| 最近中文字幕2019免费版| 精品一区二区三卡| 一区二区三区免费毛片| 亚洲在久久综合| 日韩欧美精品免费久久| 男女啪啪激烈高潮av片| 日韩,欧美,国产一区二区三区| 国产欧美日韩精品一区二区| 久久久a久久爽久久v久久| av黄色大香蕉| 国产精品麻豆人妻色哟哟久久| 久久这里有精品视频免费| 久久久精品免费免费高清| 别揉我奶头 嗯啊视频| 国产高清有码在线观看视频| 蜜臀久久99精品久久宅男| av网站免费在线观看视频| 蜜臀久久99精品久久宅男| 人妻夜夜爽99麻豆av| 国产精品一区二区性色av| 最近中文字幕高清免费大全6| 久久99热这里只频精品6学生| 99热这里只有是精品在线观看| 免费看不卡的av| 国内少妇人妻偷人精品xxx网站| 各种免费的搞黄视频| av在线天堂中文字幕| 色播亚洲综合网| 亚洲精品国产成人久久av| 国产亚洲91精品色在线| 亚洲精品日韩av片在线观看| 国产精品国产三级专区第一集| 在线观看av片永久免费下载| 中文字幕av成人在线电影| 一区二区三区精品91| 亚洲欧洲国产日韩| 免费不卡的大黄色大毛片视频在线观看| 五月伊人婷婷丁香| 三级国产精品欧美在线观看| 国产毛片a区久久久久| 你懂的网址亚洲精品在线观看| 国产av码专区亚洲av| 最近中文字幕2019免费版| av.在线天堂| 精品熟女少妇av免费看| 汤姆久久久久久久影院中文字幕| 一本色道久久久久久精品综合| 日韩欧美精品免费久久| 美女被艹到高潮喷水动态| 两个人的视频大全免费| 舔av片在线| a级毛色黄片| 人妻一区二区av| 国产午夜精品久久久久久一区二区三区| 激情 狠狠 欧美| 下体分泌物呈黄色| 亚洲精品一区蜜桃| 成年女人在线观看亚洲视频 | 国产探花极品一区二区| 久久久欧美国产精品| 99久久精品一区二区三区| 亚洲丝袜综合中文字幕| 青春草视频在线免费观看| 亚洲欧美成人精品一区二区| 波野结衣二区三区在线| 国内少妇人妻偷人精品xxx网站| 日韩中字成人| 国产精品三级大全| 亚洲精品久久午夜乱码| 欧美一区二区亚洲| 搡老乐熟女国产| 自拍欧美九色日韩亚洲蝌蚪91 | 日本欧美国产在线视频| 男人和女人高潮做爰伦理| 韩国高清视频一区二区三区| 嘟嘟电影网在线观看| 可以在线观看毛片的网站| 女人久久www免费人成看片| 亚洲第一区二区三区不卡| 久热这里只有精品99| 可以在线观看毛片的网站| 乱码一卡2卡4卡精品| 国产综合懂色| 亚洲精品一二三| 女人久久www免费人成看片| 在线看a的网站| 亚洲天堂av无毛| 中文精品一卡2卡3卡4更新| av在线天堂中文字幕| 五月天丁香电影| 国产91av在线免费观看| 中文字幕制服av| 国产69精品久久久久777片| 国产精品一区www在线观看| 久久精品夜色国产| 亚洲精品一区蜜桃| 国产探花在线观看一区二区| 亚洲欧美成人精品一区二区| 老司机影院毛片| 国产精品成人在线| 亚洲精品乱久久久久久| 最近手机中文字幕大全| 在线观看人妻少妇| 99精国产麻豆久久婷婷| 久久精品国产自在天天线| 激情 狠狠 欧美| 91精品伊人久久大香线蕉| 日韩欧美精品免费久久| 一个人观看的视频www高清免费观看| 我要看日韩黄色一级片| 色视频在线一区二区三区| 欧美激情在线99| 男人和女人高潮做爰伦理| 亚洲国产成人一精品久久久| 国产伦精品一区二区三区四那| 99视频精品全部免费 在线| 亚洲成人精品中文字幕电影| 精品久久国产蜜桃| 日韩免费高清中文字幕av| 亚洲精品国产av蜜桃| 丝袜喷水一区| 免费少妇av软件| 免费看光身美女| 日本免费在线观看一区| 成人高潮视频无遮挡免费网站| 亚洲性久久影院| 熟女电影av网| 高清午夜精品一区二区三区| 成人欧美大片| videos熟女内射| 伊人久久国产一区二区| 人妻系列 视频| av在线天堂中文字幕| 日本-黄色视频高清免费观看| 久久久精品免费免费高清| 亚洲精品国产av蜜桃| 亚洲国产精品成人久久小说| 欧美激情在线99| 国产男女内射视频| 777米奇影视久久| 中文精品一卡2卡3卡4更新| 各种免费的搞黄视频| 亚洲自偷自拍三级| 国产高清有码在线观看视频| 久久精品久久精品一区二区三区| 欧美一区二区亚洲| 亚洲精品久久午夜乱码| 18禁裸乳无遮挡免费网站照片| 人妻少妇偷人精品九色| 熟妇人妻不卡中文字幕| 18+在线观看网站| 永久免费av网站大全| 美女国产视频在线观看| 欧美激情在线99| 国产亚洲一区二区精品| 亚洲av中文字字幕乱码综合| 日韩欧美 国产精品| 一本色道久久久久久精品综合| 美女内射精品一级片tv| 成人鲁丝片一二三区免费| 观看免费一级毛片| 日韩视频在线欧美| 99久久九九国产精品国产免费| 边亲边吃奶的免费视频| 国产精品人妻久久久影院| 欧美日韩视频精品一区| 99热网站在线观看| 久久精品综合一区二区三区| 亚洲av电影在线观看一区二区三区 | 亚洲电影在线观看av| 久久热精品热| 精品久久久久久久久亚洲| 一级二级三级毛片免费看| 欧美xxxx性猛交bbbb| 一级毛片aaaaaa免费看小| 麻豆久久精品国产亚洲av| 国产一区亚洲一区在线观看| 最近手机中文字幕大全| 尤物成人国产欧美一区二区三区| 18+在线观看网站| 高清视频免费观看一区二区| av福利片在线观看| 97人妻精品一区二区三区麻豆| xxx大片免费视频| 国产成人a∨麻豆精品| 特大巨黑吊av在线直播| 国产一区二区三区综合在线观看 | 国产欧美日韩精品一区二区| 国产精品久久久久久久电影| 午夜福利网站1000一区二区三区| 午夜福利在线在线| 成人国产麻豆网| 观看美女的网站| 韩国高清视频一区二区三区| 麻豆精品久久久久久蜜桃| 国产精品久久久久久久久免| 国产探花极品一区二区| 在线观看一区二区三区| 欧美xxxx黑人xx丫x性爽| 精品亚洲乱码少妇综合久久| 蜜臀久久99精品久久宅男| 亚洲欧洲日产国产| 18禁动态无遮挡网站| 久久久久久久久久人人人人人人| 七月丁香在线播放| av国产久精品久网站免费入址| 韩国高清视频一区二区三区| 日韩视频在线欧美| 人妻一区二区av| 大陆偷拍与自拍| 最新中文字幕久久久久| 狂野欧美激情性xxxx在线观看| 亚洲国产色片| 在线观看一区二区三区激情| 日韩欧美 国产精品| 欧美少妇被猛烈插入视频| 水蜜桃什么品种好| 成人美女网站在线观看视频| 亚洲,欧美,日韩| 免费观看av网站的网址| 精品人妻偷拍中文字幕| 国产精品av视频在线免费观看| 夫妻午夜视频| 三级国产精品片| 国产成人一区二区在线| 久久久久精品久久久久真实原创| 好男人视频免费观看在线| 久久久a久久爽久久v久久| 寂寞人妻少妇视频99o| 国产伦精品一区二区三区视频9| 激情 狠狠 欧美| av免费观看日本| 精品午夜福利在线看| 搡老乐熟女国产| 一本一本综合久久| 久久精品国产a三级三级三级| 在线精品无人区一区二区三 | 久久亚洲国产成人精品v| 免费高清在线观看视频在线观看| 国产精品女同一区二区软件| 成人亚洲精品av一区二区| 久久鲁丝午夜福利片| 在线观看一区二区三区| 亚洲av一区综合| 久久久久精品性色| 国产在线一区二区三区精| 黄色一级大片看看| 七月丁香在线播放| 大香蕉97超碰在线| 日本一本二区三区精品| 免费观看av网站的网址| 九色成人免费人妻av| 亚洲美女搞黄在线观看| 精品久久久久久久久亚洲| 51国产日韩欧美| 欧美三级亚洲精品| 久久99热这里只频精品6学生| av在线app专区| 在线a可以看的网站| 最新中文字幕久久久久| 青春草亚洲视频在线观看| 国产 精品1| 免费少妇av软件| 制服丝袜香蕉在线| av又黄又爽大尺度在线免费看| 精品一区二区三区视频在线| 午夜免费观看性视频| 国产午夜精品久久久久久一区二区三区| 免费观看无遮挡的男女| av线在线观看网站| 亚洲经典国产精华液单| 国产精品福利在线免费观看| 午夜激情久久久久久久| 日韩欧美精品免费久久| 丝袜喷水一区| 老司机影院毛片| 亚州av有码| 久久人人爽人人片av| 在线看a的网站| 色视频在线一区二区三区| av在线app专区| 久久久久久九九精品二区国产| 国产黄a三级三级三级人| 性色avwww在线观看| 韩国av在线不卡| av在线亚洲专区| 国产精品偷伦视频观看了| 午夜免费鲁丝| 国产精品久久久久久精品古装| 人人妻人人爽人人添夜夜欢视频 | 99久久精品热视频| 激情五月婷婷亚洲| 亚洲人与动物交配视频| 美女被艹到高潮喷水动态| 久久久亚洲精品成人影院| 成人特级av手机在线观看| 久久久久网色| 国产免费又黄又爽又色| 只有这里有精品99| 日韩欧美精品免费久久| 真实男女啪啪啪动态图| 一级毛片黄色毛片免费观看视频| 亚洲欧美日韩无卡精品| 亚洲精品国产色婷婷电影| 成人鲁丝片一二三区免费| 女人久久www免费人成看片| 亚洲精品成人久久久久久| 能在线免费看毛片的网站| 尾随美女入室| 欧美日韩视频高清一区二区三区二| 老师上课跳d突然被开到最大视频| 日韩一本色道免费dvd| a级毛片免费高清观看在线播放| 国产乱来视频区| 日本一本二区三区精品| 国产女主播在线喷水免费视频网站| 国产av不卡久久| 国产黄色视频一区二区在线观看| 91精品一卡2卡3卡4卡| 精品人妻视频免费看| 精品久久久久久久末码| 99久久精品一区二区三区| 国产成人精品婷婷| 观看免费一级毛片| 777米奇影视久久| 99热这里只有是精品50| 久久久国产一区二区| 成人美女网站在线观看视频| 99热国产这里只有精品6| 国产人妻一区二区三区在| 精品酒店卫生间| 九色成人免费人妻av| 日韩一区二区三区影片| 99热这里只有是精品50| 久久热精品热| 18禁在线播放成人免费| 香蕉精品网在线| 国产男女超爽视频在线观看| 少妇的逼好多水| 成人亚洲精品av一区二区| 亚洲一区二区三区欧美精品 | 成人欧美大片| 精品人妻偷拍中文字幕| 各种免费的搞黄视频| 久久久久网色| 国产精品人妻久久久久久| 黄片无遮挡物在线观看| 深爱激情五月婷婷| 欧美高清性xxxxhd video| 亚洲精品自拍成人| 午夜老司机福利剧场| av福利片在线观看| 91久久精品国产一区二区成人| 成年女人在线观看亚洲视频 | 在线 av 中文字幕| 可以在线观看毛片的网站| 男人舔奶头视频| 国产精品精品国产色婷婷| 亚洲人成网站在线播| 国产免费一级a男人的天堂| 在线观看美女被高潮喷水网站| 亚洲av在线观看美女高潮| 男女边摸边吃奶| 日韩欧美精品免费久久| 免费看日本二区| 亚洲av日韩在线播放| 在线观看美女被高潮喷水网站| 黄色欧美视频在线观看| 99视频精品全部免费 在线| 激情 狠狠 欧美| 99热这里只有是精品50| 精品人妻视频免费看| 国产淫片久久久久久久久| 亚洲性久久影院| 少妇人妻久久综合中文| 国产黄色视频一区二区在线观看| 久久精品国产亚洲av涩爱| 成人国产av品久久久| 国产黄色免费在线视频| 成人鲁丝片一二三区免费| 哪个播放器可以免费观看大片| 美女高潮的动态| 永久免费av网站大全| 成人无遮挡网站| 日韩成人伦理影院| 99热网站在线观看| 最近2019中文字幕mv第一页| 免费看光身美女| 国产成人福利小说| 99热这里只有精品一区| 亚洲精品aⅴ在线观看| 啦啦啦中文免费视频观看日本| 一级爰片在线观看| 少妇人妻 视频| 亚洲欧美日韩另类电影网站 | 欧美成人午夜免费资源| 大陆偷拍与自拍| 97超视频在线观看视频| a级毛色黄片| 国产精品.久久久| 欧美另类一区| 别揉我奶头 嗯啊视频| 天天躁夜夜躁狠狠久久av| 亚洲aⅴ乱码一区二区在线播放| 美女被艹到高潮喷水动态| 免费看光身美女| 久久久精品欧美日韩精品| 一级爰片在线观看| 日韩中字成人| 亚洲欧美一区二区三区国产| 看十八女毛片水多多多| 亚洲国产欧美人成| 亚洲天堂国产精品一区在线| 成人黄色视频免费在线看| 国产精品熟女久久久久浪| 亚洲国产日韩一区二区| 久久久精品94久久精品| 真实男女啪啪啪动态图| 少妇人妻 视频| 亚洲精品中文字幕在线视频 | 久久久亚洲精品成人影院| 久久综合国产亚洲精品| 亚洲av成人精品一区久久| 少妇丰满av| 久久精品夜色国产| 人妻夜夜爽99麻豆av| 日日摸夜夜添夜夜爱| 欧美+日韩+精品| a级毛片免费高清观看在线播放| 新久久久久国产一级毛片| 黄片无遮挡物在线观看| 一级毛片我不卡| 成人欧美大片| 免费av毛片视频| 免费看av在线观看网站| av网站免费在线观看视频| 王馨瑶露胸无遮挡在线观看| 尾随美女入室| 国产成人福利小说| 秋霞伦理黄片| 2018国产大陆天天弄谢| 久久久久久久久久久丰满| 久久久久精品久久久久真实原创| 少妇熟女欧美另类| 国产高清不卡午夜福利| 舔av片在线| 成人综合一区亚洲| 日韩精品有码人妻一区| 99热全是精品| 国产免费福利视频在线观看| 91在线精品国自产拍蜜月| 午夜免费男女啪啪视频观看| 久久人人爽人人爽人人片va| 日韩在线高清观看一区二区三区| 国产成人a∨麻豆精品| 美女脱内裤让男人舔精品视频| 美女内射精品一级片tv| 网址你懂的国产日韩在线| 久久97久久精品| 人妻 亚洲 视频| 久久久久久伊人网av| 国产亚洲一区二区精品| 亚洲欧美一区二区三区黑人 | 日日摸夜夜添夜夜爱| 天天一区二区日本电影三级| 亚洲精品国产av成人精品| 亚洲精品亚洲一区二区| 成年av动漫网址| 成人国产av品久久久| 欧美成人精品欧美一级黄| 在线看a的网站| 男插女下体视频免费在线播放| 亚洲经典国产精华液单| 一二三四中文在线观看免费高清| 永久网站在线| 一级av片app| 免费电影在线观看免费观看| 国产欧美日韩精品一区二区| 国产黄a三级三级三级人| 精品视频人人做人人爽| 久久人人爽人人爽人人片va| 久久精品国产鲁丝片午夜精品| 国产精品无大码| 在线免费十八禁| 国产精品一区二区三区四区免费观看| 亚洲av电影在线观看一区二区三区 | 99热这里只有是精品50| 日韩不卡一区二区三区视频在线| 女人被狂操c到高潮| 亚洲欧美精品专区久久| 亚洲精品一区蜜桃| 国产视频首页在线观看| 成人亚洲欧美一区二区av| 在线免费观看不下载黄p国产| 日日啪夜夜撸| 日本黄色片子视频| 自拍偷自拍亚洲精品老妇| 一级黄片播放器| 国产永久视频网站| 在现免费观看毛片| 寂寞人妻少妇视频99o| 国产精品久久久久久av不卡| 国产伦理片在线播放av一区| 久久久久久久午夜电影| 日产精品乱码卡一卡2卡三| 性色av一级| 国产精品福利在线免费观看| 日韩一区二区三区影片| 日韩一本色道免费dvd| 最近的中文字幕免费完整| 日韩欧美 国产精品| 亚洲欧美一区二区三区国产| 日日啪夜夜撸| 久久99热6这里只有精品| 神马国产精品三级电影在线观看| 免费黄色在线免费观看| 中文天堂在线官网| 热re99久久精品国产66热6| 草草在线视频免费看| 国产又色又爽无遮挡免| 一本久久精品| 亚洲自拍偷在线| 又爽又黄a免费视频| 自拍欧美九色日韩亚洲蝌蚪91 | 婷婷色综合大香蕉| 国产精品久久久久久久久免| 2021天堂中文幕一二区在线观| 韩国av在线不卡| 国产午夜精品久久久久久一区二区三区| 久久久久性生活片| 国产黄色免费在线视频| 日韩一区二区三区影片| 2021天堂中文幕一二区在线观| 国产69精品久久久久777片| 成人黄色视频免费在线看| 亚洲av不卡在线观看| 各种免费的搞黄视频| 亚洲国产日韩一区二区| 国产精品99久久99久久久不卡 | 97超碰精品成人国产| 丰满乱子伦码专区| 久久这里有精品视频免费| 一级a做视频免费观看| 成人午夜精彩视频在线观看| 欧美性猛交╳xxx乱大交人| 精品少妇黑人巨大在线播放| 色播亚洲综合网| 在线a可以看的网站| 亚洲精品国产av成人精品| 国产高潮美女av| 99热全是精品| 色吧在线观看| 一个人看视频在线观看www免费| 亚洲国产色片| www.色视频.com| 国产精品99久久久久久久久| 国产片特级美女逼逼视频| 亚洲美女视频黄频| 91精品一卡2卡3卡4卡| 99久国产av精品国产电影| 亚洲av国产av综合av卡| 特大巨黑吊av在线直播| 亚洲av男天堂| 99久久精品一区二区三区| 亚洲精品国产色婷婷电影| 日韩在线高清观看一区二区三区| 一级毛片久久久久久久久女| 免费看光身美女|