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

    衡水地區(qū)深層地下水水化學(xué)特征及其演化過程

    2018-07-03 11:49:46張振國何江濤
    現(xiàn)代地質(zhì) 2018年3期
    關(guān)鍵詞:特征研究

    張振國,何江濤,王 磊,彭 聰

    (1.中國地質(zhì)大學(xué)(北京) 水資源與環(huán)境學(xué)院,北京 100083;2.江蘇省地質(zhì)調(diào)查研究院,江蘇 南京 210018)

    0 引 言

    衡水地區(qū)地處河北平原中部,水資源匱乏,人均水資源量只有148 m3,是河北省最缺水的地區(qū)之一[1]。由于地表水資源短缺,地下水是居民生活用水和工農(nóng)業(yè)用水的主要來源。20世紀(jì)70年代以來,隨著社會(huì)經(jīng)濟(jì)的發(fā)展,大量超采深層淡水導(dǎo)致區(qū)域深層地下水位逐年下降,形成了以衡水市區(qū)為中心的區(qū)域降落漏斗,改變了初始的水動(dòng)力場(chǎng)和水化學(xué)場(chǎng),并導(dǎo)致了地面沉降、咸水下移、深層地下水質(zhì)惡化等一系列環(huán)境地質(zhì)問題,嚴(yán)重阻礙了當(dāng)?shù)亟?jīng)濟(jì)社會(huì)的發(fā)展,給居民生活帶來了諸多不便。

    近年來已有學(xué)者對(duì)超采情況下衡水地區(qū)的地下水問題進(jìn)行了研究[2-5]。程雙虎[3]基于地下水均衡理論建立了地下水均衡模型,明確了深層地下水開采量主要來源于側(cè)向徑流補(bǔ)給、越流補(bǔ)給和黏性土壓密釋水。高業(yè)新[4]通過抽水試驗(yàn)結(jié)合水化學(xué)數(shù)據(jù)分析發(fā)現(xiàn)淺層咸水對(duì)深層淡水的越流補(bǔ)給主要發(fā)生在第一含水層組和第二含水層組之間,第三、第四、第五含水層組開采量主要來源于儲(chǔ)存量和側(cè)向徑流補(bǔ)給。宋海波[5]通過空間分析、水質(zhì)分析和數(shù)理統(tǒng)計(jì)等方法研究了地下水動(dòng)力場(chǎng)及咸淡水界面的演化,發(fā)現(xiàn)衡水市咸淡水界面以年均0.4 m的速度下移。這些研究的關(guān)注點(diǎn)主要集中在水質(zhì)和水量上,對(duì)開采條件下水化學(xué)場(chǎng)的演變關(guān)注較少。咸水下移主要導(dǎo)致與其相鄰的第二含水層組中地下水TDS及Cl-、Na+濃度升高,對(duì)第三、第四、第五含水層組影響較小。作為主開采層的第三含水層水化學(xué)特征較大規(guī)模開采前出現(xiàn)了明顯的變化,其演變機(jī)制尚不明了。本文通過分析區(qū)域水文地球化學(xué)特征,分析了研究區(qū)以第三含水層組為代表的深層地下水在過量開采條件下的水化學(xué)演化特征,探討了水化學(xué)組分的來源及形成機(jī)制。

    1 研究區(qū)概況

    1.1 自然地理環(huán)境

    研究區(qū)位于華北平原中部衡水市及其周邊地區(qū),屬于溫帶半濕潤大陸性季風(fēng)氣候區(qū),四季分明,多年平均降水量約522 mm,主要集中在汛期的6—9月份[5]。研究區(qū)內(nèi)地勢(shì)平坦,東部屬于山前沖洪積傾斜平原末端,中部和西部屬于沖積湖積平原。滏陽河是流經(jīng)該地區(qū)的最大河流,自西南向東北流經(jīng)研究區(qū)南部及東部,滹沱河自西向東流經(jīng)研究區(qū)北部。

    1.2 水文地質(zhì)條件

    本區(qū)屬河北平原水文地質(zhì)區(qū),第四系堆積厚度一般為400~600 m,巖層沉積顆粒由西向東逐漸變細(xì)[6]。第四系松散巖層構(gòu)成的孔隙含水巖組為主要含水巖組,由上而下可劃分為4個(gè)以一定厚度黏土層相隔的含水層組[7],通常稱第一含水層組賦存的地下水為淺層地下水,第二、第三、第四含水層組賦存的地下水為深層地下水。由于大陸鹽化和海侵作用,華北平原中東部地區(qū)淺部普遍分布一層咸水體,研究區(qū)在安平—寧晉一線以西為全淡水區(qū),以東為有咸水區(qū)。

    研究區(qū)地下水的天然流向整體為自西南流向東北,北部則由西向東[5]。淺層地下水接受大氣降水和地表水入滲補(bǔ)給,通過蒸發(fā)、側(cè)向徑流等方式排泄。深層地下水接受緩慢側(cè)向徑流補(bǔ)給,通過徑流和托頂越流排泄[4]。20世紀(jì)70年代以來,大量超采深層地下水形成了區(qū)域地下水降落漏斗,改變了初始的地下水流場(chǎng)。深層地下水除西北部流向自西向東外,其他大部分區(qū)域流向漏斗中心區(qū)(圖1)[5]。同時(shí),在抽采活動(dòng)強(qiáng)烈的地區(qū),深層含水層產(chǎn)生了明顯的水平徑流,造成其水動(dòng)力條件好于淺層含水層[1]。補(bǔ)給來源與含水層組的埋深有關(guān),第二含水層組主要來自凈儲(chǔ)量以及越流補(bǔ)給,其他含水層組主要來源于側(cè)向徑流補(bǔ)給[4]。排泄方式主要為人工開采排泄。

    圖1 研究區(qū)水文地質(zhì)簡(jiǎn)圖及采樣點(diǎn)分布(等水頭線據(jù)張兆吉等[30]修改)Fig.1 Hydrogeological sketch map of the study area and distribution of groundwater sampling sites

    2 材料與方法

    2.1 樣品采集與處理

    為充分反映研究區(qū)深層地下水水化學(xué)特征及其隨時(shí)間和空間的演化規(guī)律,本文使用了2個(gè)時(shí)期的地下水監(jiān)測(cè)數(shù)據(jù):歷史時(shí)期(1974—1978年)和現(xiàn)階段(2007—2008年)數(shù)據(jù)。現(xiàn)階段數(shù)據(jù)來源于全國地下水污染調(diào)查評(píng)價(jià)綜合研究項(xiàng)目,共124組,采樣井主要為民用井,井深在170~350 m之間。采樣時(shí)現(xiàn)場(chǎng)測(cè)定水溫、氧化還原電位(Eh)、酸堿度(pH)、電導(dǎo)率和溶解氧等現(xiàn)場(chǎng)指標(biāo)。實(shí)驗(yàn)室測(cè)定總硬度、Ca2+、Mg2+、Na+、K+、Cl-、SO42-、HCO3-、CO32-等指標(biāo),溶解性總固體(TDS)為計(jì)算得出。歷史時(shí)期數(shù)據(jù)來源于衡水地區(qū)地下水質(zhì)年鑒[8],采樣井主要為勘探孔,井深在150~350 m之間,共48組。根據(jù)采樣井深度和水文地質(zhì)條件可知,兩個(gè)時(shí)期的地下水樣品均采自第三含水層組,采樣點(diǎn)位置如圖1所示。

    2.2 研究方法

    本文通過對(duì)水質(zhì)分析結(jié)果進(jìn)行數(shù)理統(tǒng)計(jì)并繪制Piper三線圖,分析研究區(qū)深層地下水的水化學(xué)特征。利用SPSS19統(tǒng)計(jì)分析軟件Q型聚類分析模塊對(duì)研究區(qū)進(jìn)行水化學(xué)特征分區(qū),結(jié)合Gibbs圖、氯堿指數(shù)、離子相關(guān)關(guān)系方法分析各水化學(xué)分區(qū)的地下水水化學(xué)特征及其形成機(jī)制;對(duì)比不同時(shí)期水動(dòng)力場(chǎng)和水化學(xué)場(chǎng)的響應(yīng)關(guān)系,探討研究區(qū)地下水水化學(xué)演化規(guī)律,分析人類活動(dòng)對(duì)地下水水化學(xué)場(chǎng)的影響。

    3 結(jié)果分析與討論

    3.1 水化學(xué)分區(qū)及特征

    3.1.1 水化學(xué)分區(qū)

    圖2 研究區(qū)水化學(xué)分區(qū)圖(等水頭線據(jù)張兆吉和費(fèi)宇紅[30],2009)Fig.2 Hydrochemical fields of the study area

    圖3 現(xiàn)階段各區(qū)深層地下水主要指標(biāo)含量特征Fig.3 Characteristics of the major indicative contents of deep groundwater at present

    不同區(qū)域水文地質(zhì)條件等因素不同導(dǎo)致地下水水化學(xué)特征具有一定的叢集性[9-10]。本文利用SPSS19統(tǒng)計(jì)分析軟件[11],根據(jù)pH、Ca2+、Mg2+、Na++K+、Cl-、SO42-和HCO3-共7項(xiàng)指標(biāo)對(duì)現(xiàn)階段124組數(shù)據(jù)進(jìn)行了Q型聚類分析。樣品分為3類,各類樣品數(shù)分別為23個(gè)、51個(gè)和50個(gè)。按分類情況將研究區(qū)劃分為3個(gè)水化學(xué)分區(qū)(圖2)。其中,Ⅰ區(qū)位于研究區(qū)東部,靠近太行山前沖積傾斜平原末端,為補(bǔ)給區(qū);Ⅲ區(qū)位于冀棗衡降落漏斗中心區(qū),為排泄區(qū);Ⅱ區(qū)位于Ⅰ區(qū)和Ⅲ區(qū)之間的饒陽—深州—寧晉一帶,在Ⅲ區(qū)以東武強(qiáng)—武邑—冀州一帶也有分布,為徑流區(qū)。

    3.1.2 各區(qū)水化學(xué)特征

    圖4 歷史時(shí)期各區(qū)深層地下水主要指標(biāo)箱型圖Fig.4 Box plot of the major indicative contents of deep groundwater in the past

    圖5 深層地下水歷史時(shí)期(a)和現(xiàn)階段(b)Piper三線圖Fig.5 Piper diagrams of historic (a) and present (b) deep groundwater

    現(xiàn)階段各區(qū)TDS值表現(xiàn)為Ⅲ區(qū)>Ⅱ區(qū)>Ⅰ區(qū),這符合地下水徑流路徑越長,溶濾作用越充分,TDS越高的特征。陽離子濃度由Ⅰ區(qū)到Ⅲ區(qū)均呈上升趨勢(shì),具體表現(xiàn)為Na+升高速率高于Ca2+、Mg2+升高速率。陰離子濃度由Ⅰ區(qū)到Ⅲ區(qū)Cl-呈上升趨勢(shì),HCO3-呈下降趨勢(shì),SO42-濃度由Ⅰ區(qū)到Ⅱ區(qū)上升,Ⅱ到Ⅲ區(qū)無明顯變化(圖3)。歷史時(shí)期TDS表現(xiàn)為Ⅲ區(qū)略高于Ⅱ區(qū)。陽離子濃度由Ⅱ區(qū)到Ⅲ區(qū)Na+呈上升趨勢(shì),兩個(gè)區(qū)Ca2+、Mg2+濃度均較低,以Na+為主。陰離子濃度由Ⅱ區(qū)到Ⅲ區(qū)HCO3-和Cl-呈微弱的上升趨勢(shì),SO42-略有下降趨勢(shì)(圖4)。結(jié)合Piper圖可知,現(xiàn)階段各區(qū)地下水水化學(xué)類型具有明顯的叢集性,其中Ⅰ區(qū)水化學(xué)類型主要是HCO3-Na、HCO3·SO4-Na和HCO3-Ca,Ⅱ區(qū)水化學(xué)類型主要是SO4·Cl-Na和Cl·SO4-Na,Ⅲ區(qū)水化學(xué)類型主要是Cl-Na。由Ⅰ區(qū)到Ⅲ區(qū)地下水水化學(xué)類型基本符合由補(bǔ)給區(qū)到排泄區(qū)的水化學(xué)演化規(guī)律。歷史時(shí)期Ⅱ區(qū)和Ⅲ區(qū)地下水水化學(xué)特征沒有明顯的分帶性,主要是因?yàn)樵诮德渎┒樊a(chǎn)生前Ⅱ區(qū)和Ⅲ區(qū)深層地下水由西、西南向東、東北緩慢側(cè)向徑流,水交替較弱,沒有明顯的補(bǔ)排關(guān)系。其中Ⅱ區(qū)水化學(xué)類型主要是Cl·HCO3-Na和SO4·Cl-Na,Ⅲ區(qū)水化學(xué)類型主要是Cl·HCO3-Na和Cl-Na(圖5)。

    3.2 人類活動(dòng)影響下水化學(xué)場(chǎng)時(shí)空演化特征

    3.2.1 水化學(xué)演化特征

    對(duì)比不同時(shí)期研究區(qū)地下水水化學(xué)特征及演化規(guī)律可知,兩個(gè)時(shí)期Ⅱ、Ⅲ區(qū)陽離子均以Na+為主。現(xiàn)階段陰離子Ⅱ區(qū)以SO42-占優(yōu)勢(shì),Ⅲ區(qū)以Cl-占優(yōu)勢(shì),歷史時(shí)期均為Cl-、SO42-占優(yōu)勢(shì)?,F(xiàn)階段陽離子濃度由Ⅱ區(qū)到Ⅲ區(qū)均呈上升趨勢(shì),歷史時(shí)期Na+呈上升趨勢(shì),Ca2+、Mg2+呈下降趨勢(shì)?,F(xiàn)階段陰離子濃度由Ⅱ區(qū)到Ⅲ區(qū)Cl-呈上升趨勢(shì),HCO3-呈下降趨勢(shì),SO42-由Ⅰ區(qū)到Ⅱ區(qū)上升,Ⅱ到Ⅲ區(qū)無明顯變化。歷史時(shí)期Ⅱ區(qū)到Ⅲ區(qū)HCO3-和Cl-呈微弱上升趨勢(shì),SO42-略有下降。TDS在兩個(gè)時(shí)期由Ⅱ區(qū)到Ⅲ區(qū)均呈上升趨勢(shì),現(xiàn)階段升高速率高于歷史時(shí)期。TDS、Na++K+和Cl-濃度由Ⅰ區(qū)到Ⅲ區(qū)呈上升趨勢(shì),且現(xiàn)階段上升速率高于歷史時(shí)期,表明地下水降落漏斗的產(chǎn)生和發(fā)展加強(qiáng)了地下水對(duì)圍巖的溶濾作用。歷史時(shí)期HCO3-濃度由Ⅱ區(qū)到Ⅲ區(qū)略有上升,而現(xiàn)階段呈下降趨勢(shì),可能是漏斗區(qū)水頭下降、壓力降低、脫碳酸作用發(fā)生所導(dǎo)致。兩個(gè)時(shí)期地下水水化學(xué)類型的時(shí)空演變特征為(圖5):歷史時(shí)期Ⅱ區(qū)到Ⅲ區(qū)地下水化學(xué)類型沒有明顯分帶性,主要為Cl·HCO3-Na、SO4·Cl-Na和Cl-Na。現(xiàn)階段則出現(xiàn)了明顯分帶性,Ⅱ區(qū)轉(zhuǎn)化為SO4·Cl-Na和Cl·SO4-Na,Ⅲ區(qū)轉(zhuǎn)化為Cl-Na。

    3.2.2 人類活動(dòng)影響

    冀棗衡漏斗形成于1968年[12-13],隨采水量持續(xù)增加,漏斗發(fā)展很快,在1980—2003年的24年間,衡水市區(qū)及南部的巨鹿、南宮等地深層地下水位下降了約45 m,到2005年衡水市區(qū)地下水位埋深達(dá)到87.5 m(圖6)。在降落漏斗產(chǎn)生前研究區(qū)深層地下水由西、西南向東、東北緩慢徑流,地下水化學(xué)類型主要為HCO3·Cl-Na和Cl·HCO3-Na。隨降落漏斗形成和發(fā)展,水力坡度變大,水交替變強(qiáng),加劇了地下水對(duì)圍巖的溶濾作用。地下水對(duì)圍巖的溶濾可分為兩個(gè)階段,沿徑流路徑初期圍巖中鹽巖的易溶組分快速進(jìn)入地下水中導(dǎo)致地下水TDS升高,隨溶濾作用的進(jìn)行,易溶組分溶濾完全,難溶的碳酸鹽巖等開始溶解,這時(shí)地下水TDS開始降低。研究區(qū)深層地下水14C測(cè)年(圖7)顯示,由Ⅰ區(qū)到Ⅲ區(qū)地下水年齡逐漸增大,且地下水TDS由Ⅰ區(qū)到Ⅲ區(qū)逐漸升高,表明其處于溶濾作用的第一階段。同時(shí),漏斗區(qū)水頭下降,壓力降低,導(dǎo)致了脫碳酸作用的發(fā)生,使Ⅰ區(qū)到Ⅲ區(qū)HCO3-濃度逐漸降低。地下水降落漏斗的產(chǎn)生導(dǎo)致地下水環(huán)境發(fā)生了變化,新的補(bǔ)徑排關(guān)系及其強(qiáng)度使研究區(qū)原有的水化學(xué)場(chǎng)發(fā)生了改變,Ⅱ區(qū)地下水水化學(xué)類型由Cl·HCO3-Na、HCO3·Cl-Na和SO4·Cl-Na轉(zhuǎn)化為SO4·Cl-Na和Cl·SO4-Na,Ⅲ區(qū)地下水水化學(xué)類型由Cl·HCO3-Na和Cl-Na轉(zhuǎn)化為Cl-Na(圖6)。

    圖6 研究區(qū)深層地下水埋深及1980—2003年水位變差圖(據(jù)張兆吉等修改[30])Fig.6 Potentiometric surface and depth variation of the Hengshui deep groundwater (1980—2003)

    圖7 石家莊—衡水2005年地下水14C年齡剖面圖(據(jù)張兆吉等[30]修改)Fig.7 Shijiazhuang-Hengshui 14C age profile (modified from Zhang et al.[30])

    圖8 研究區(qū)深層地下水Gibbs圖Fig.8 Gibbs diagrams of the Hengshui deep groundwater

    3.3 水化學(xué)組分來源及形成作用

    地下水水化學(xué)組分的形成受氣候條件、水文地質(zhì)條件、地下水補(bǔ)給來源等多種因素共同影響,其主要形成作用有溶濾作用、陽離子交替吸附作用、氧化還原作用、混合作用、濃縮作用和脫碳酸作用[14-18]。根據(jù)地下水中各離子之間的關(guān)系,可以解釋水化學(xué)組分的來源[19-25]。將兩個(gè)時(shí)期的水化學(xué)數(shù)據(jù)投繪在Gibbs圖[26](圖8)上,可以看出各時(shí)期地下水TDS主要集中在400~1 000 mg/L,屬中等含量,表明該含水層組不同時(shí)期地下水水化學(xué)組分受蒸發(fā)濃縮和大氣降水影響均較小。各區(qū)地下水ρ(Cl-)/ρ(Cl-+HCO3-)比值具有較好的集群性,ρ(Na+)/ρ(Na++Ca2+)比值則沒有明顯的規(guī)律性,這與陽離子交替吸附作用的影響有關(guān)?,F(xiàn)階段Ⅰ區(qū)地下水ρ(Cl-)/ρ(Cl-+HCO3-)比值均小于0.2,Ⅱ區(qū)主要分布在0.2~0.8之間,Ⅲ區(qū)均大于0.6,表明由Ⅰ區(qū)到Ⅲ區(qū)受溶濾作用影響Cl-占比逐漸升高。歷史時(shí)期Ⅱ區(qū)和Ⅲ區(qū)地下水ρ(Cl-)/ρ(Cl-+HCO3-)比值規(guī)律性不明顯。

    3.3.1 溶濾作用

    研究區(qū)深層地下水Gibbs圖顯示巖石風(fēng)化是控制該含水層組地下水水化學(xué)組分形成的主要因素之一。通過地下水中HCO3-/Na+和Ca2+/Na+摩爾比值可以進(jìn)一步判斷其受不同礦物風(fēng)化影響的程度。當(dāng)HCO3-/Na+、Ca2+/Na+的摩爾比值較高時(shí),說明地下水水化學(xué)成分受碳酸鹽礦物溶解的影響較大;反之,受巖鹽溶解影響較大[27]。從研究區(qū)深層地下水HCO3-/Na+-Ca2+/Na+關(guān)系圖(圖9)中可以看出,整體上現(xiàn)階段地下水中HCO3-/Na+、Ca2+/Na+的摩爾比值由Ⅰ區(qū)到Ⅲ區(qū)逐漸降低,歷史時(shí)期則沒有明顯變化,表明現(xiàn)階段深層地下水在降落漏斗影響下由Ⅰ區(qū)到Ⅲ區(qū)對(duì)碳酸鹽礦物的溶濾逐漸減弱,對(duì)巖鹽礦物的溶濾逐漸增強(qiáng),反映在地下水組分濃度關(guān)系上由Ⅰ區(qū)到Ⅲ區(qū)Na+、Cl-升高速率遠(yuǎn)高于Ca2+、Mg2+。

    圖9 研究區(qū)深層地下水HCO3-/Na+、Ca2+/Na+關(guān)系圖Fig.9 Hydrochemical relationships of HCO3- vs.Na+ and Ca2+ vs.Na+ of the Hengshui deep groundwater

    圖10 研究區(qū)深層地下水氯堿指數(shù)與TDS關(guān)系圖Fig.10 Chlor-alkali index of the Hengshui deep groundwater

    3.3.2 陽離子交替吸附作用

    陽離子交替吸附作用的強(qiáng)度和交換方向受巖石粒度、離子性質(zhì)和離子濃度等因素影響[18]。本文引入氯堿指數(shù)(CAI1和CAI2)[28-29]進(jìn)一步分析研究區(qū)地下水陽離子交替吸附作用發(fā)生的方向和強(qiáng)度(圖10),除Ⅲ區(qū)少數(shù)樣品外,三個(gè)區(qū)CAI1和CAI2均小于0,且表現(xiàn)為Ⅲ區(qū)>Ⅱ區(qū)>Ⅰ區(qū),表明研究區(qū)地下水中Ca2+、Mg2+與圍巖中Na+、K+發(fā)生了交換,由Ⅰ區(qū)到Ⅲ區(qū)交換強(qiáng)度逐漸減弱,Ⅲ區(qū)少數(shù)樣品CAI1和CAI2大于0,表明Ⅲ區(qū)部分地區(qū)地下水中Na+、K+與圍巖中Ca2+、Mg2+發(fā)生了反向陽離子交換,反映在地下水組分關(guān)系上表現(xiàn)為Ⅱ區(qū)到Ⅲ區(qū)Na+質(zhì)量濃度的升高速率低于Ⅰ區(qū)到Ⅱ區(qū),同時(shí)Na+摩爾濃度大于Cl-摩爾濃度,表明除巖鹽溶解外Na+還有其他來源,如鈉長石(NaAlSi3O8)。

    (1)

    (2)

    3.3.3 脫碳酸作用

    研究區(qū)由Ⅰ區(qū)到Ⅲ區(qū)隨地下水降落漏斗加深,地下水中壓力逐漸降低,且pH逐漸升高,由碳酸平衡可知地下水中CO2含量也應(yīng)逐漸降低。碳酸鹽類在水中的溶解度取決于水中所溶解CO2的量,當(dāng)溫度升高或壓力減小時(shí),水中溶解的CO2減少,HCO3-與Ca2+、Mg2+等結(jié)合產(chǎn)生碳酸鹽沉淀[18]。離子含量變化表現(xiàn)為:歷史時(shí)期Ⅱ區(qū)到Ⅲ區(qū)HCO3-質(zhì)量濃度升高,而現(xiàn)階段Ⅰ區(qū)到Ⅲ區(qū)HCO3-質(zhì)量濃度逐漸降低,Ca2+、Mg2+質(zhì)量濃度略有升高,同時(shí)SO42-質(zhì)量濃度升高(圖4),說明現(xiàn)階段深層地下水在降落漏斗影響下Ⅱ區(qū)和Ⅲ區(qū)發(fā)生了脫碳酸作用,消耗了HCO3-、Ca2+和Mg2+。Ca2+、Mg2+質(zhì)量濃度的升高表明地下水混入了其他來源,結(jié)合SO42-質(zhì)量濃度的升高,推測(cè)地下水中發(fā)生了石膏的溶解。

    4 結(jié) 論

    (1)地下水降落漏斗的產(chǎn)生及發(fā)展導(dǎo)致研究區(qū)水動(dòng)力場(chǎng)發(fā)生了改變,形成了新的補(bǔ)徑排關(guān)系,水化學(xué)場(chǎng)也隨之改變,出現(xiàn)了分帶性。由補(bǔ)給區(qū)到排泄區(qū),TDS逐漸變大,水化學(xué)類型由HCO3-Na、HCO3·SO4-Na型逐步演變?yōu)镾O4·Cl-Na和Cl·SO4-Na型,排泄區(qū)為Cl-Na型。

    (2)超采深層地下水對(duì)研究區(qū)深層地下水水化學(xué)演化影響較大,隨地下水降落漏斗形成和發(fā)展,水力坡度變大,地下水滲透速度升高,地下水對(duì)圍巖的溶濾加劇,導(dǎo)致溶解組分含量整體升高。漏斗區(qū)水頭下降,壓力降低,地下水發(fā)生脫碳酸作用。在漏斗形成后,Ⅱ區(qū)地下水水化學(xué)類型由Cl·HCO3-Na、HCO3·Cl-Na和SO4·Cl-Na型轉(zhuǎn)化為SO4·Cl-Na和Cl·SO4-Na型,Ⅲ區(qū)地下水水化學(xué)類型由Cl·HCO3-Na和Cl-Na型轉(zhuǎn)化為Cl-Na型。

    (3)研究區(qū)深層地下水水化學(xué)組分主要受巖石風(fēng)化、陽離子交替吸附和脫碳酸作用控制。由補(bǔ)給區(qū)到排泄區(qū)表現(xiàn)為:地下水對(duì)碳酸鹽礦物的溶濾逐漸減弱,對(duì)巖鹽礦物的溶濾逐漸增強(qiáng);陽離子交換強(qiáng)度逐漸減弱,在排泄區(qū)部分地區(qū)地下水中Na+、K+與圍巖中Ca2+、Mg2+發(fā)生了反向陽離子交換;在徑流區(qū)發(fā)生了脫碳酸作用和石膏的溶解。

    參考文獻(xiàn):

    [1] 李秀瑩.衡水市深層地下水水環(huán)境狀況分析[J].地下水,2011,33(6):29-30.

    [2] 李莎, 成建梅, 李敏敏, 等. 人類開采活動(dòng)影響下的衡水地區(qū)地下水水質(zhì)特征及演化[J]. 南水北調(diào)與水利科技,2016,14(3): 55-61, 100.

    [3] 程雙虎, 李明良, 許健, 等. 衡水市區(qū)深層地下水資源構(gòu)成分析[J]. 地下水, 2016,35(5): 32-33.

    [4] 高業(yè)新. 華北平原典型地區(qū)大規(guī)模開采條件下不同層位含水組地下水互動(dòng)關(guān)系研究[D]. 北京: 中國地質(zhì)科學(xué)院, 2008: 100-119.

    [5] 宋海波. 華北平原典型區(qū)地下水質(zhì)及咸淡水界面變化研究[D]. 北京: 中國地質(zhì)科學(xué)院, 2006: 1-78.

    [6] 張宗祜,李烈榮,中國地下水資源/河北卷[M].北京:中國地圖出版社, 2005: 5-12.

    [7] 陳望和. 河北地下水[M]. 北京: 地震出版社, 1999: 5-25.

    [8] 河北省水文地質(zhì)工程地質(zhì)觀測(cè)總站.河北省衡水地區(qū)地下水質(zhì)年鑒(1967—1980)[R].石家莊:河北省水文地質(zhì)工程地質(zhì)觀測(cè)總站,1981: 1-176.

    [9] 王凱. 晉祠泉域巖溶地下水水化學(xué)特征及水文地球化學(xué)模擬[D]. 太原:太原理工大學(xué), 2016: 47-66.

    [10] 竇妍. 鄂爾多斯盆地北部白堊系地下水水文地球化學(xué)演化及循環(huán)規(guī)律研究[D]. 西安:長安大學(xué), 2010: 64-78.

    [11] 鄧維斌. SPSS 19 (中文版) 統(tǒng)計(jì)分析實(shí)用教程[M]. 北京:電子工業(yè)出版社,2012: 203-217.

    [12] 張宗祜. 華北平原地下水環(huán)境演化[M]. 北京:地質(zhì)出版社, 2000: 36-46.

    [13] 河北省地質(zhì)局第九地質(zhì)大隊(duì).衡水地區(qū)“冀棗衡地下水位下降漏斗”的發(fā)展情況及其防治意見[R].衡水:河北省地質(zhì)局第九地質(zhì)大隊(duì),1976: 1-10.

    [14] PARLANGE J Y. Concepts and models in groundwater hydrology[J]. American Scientist,1971(4):494-495.

    [15] SCHUH W M, KLINKEBIEL D L, GARDNER J C, et al. Tracer and nitrate movement to groundwater in the Northern Great Plains[J]. Journal of Environmental Quality, 1997, 26(5): 1335-1347.

    [16] GAOFENG Z, YONGHONG S, CHUNLIN H, et al. Hydrogeochemical processes in the groundwater environment of Heihe River Basin, northwest China[J]. Environmental Earth Sciences, 2010, 60(1): 139-153.

    [17] 郭清海, 王焰新. 典型新生代斷陷盆地內(nèi)孔隙地下水地球化學(xué)過程及其模擬:以山西太原盆地為例[J]. 地學(xué)前緣,2014,25(4):83-90.

    [18] 李學(xué)禮. 水文地球化學(xué)[M]. 北京:原子能出版社, 1988: 148-156.

    [19] 蔣萬軍, 趙丹, 王廣才, 等. 新疆吐—哈盆地地下水水文地球化學(xué)特征及形成作用[J]. 現(xiàn)代地質(zhì), 2016, 30(4): 825-833.

    [20] 蔡月梅, 王文祥, 張明江, 等. 吐魯番盆地地表水和地下水主要化學(xué)成分的形成與演化[J]. 現(xiàn)代地質(zhì),2016,30(3):680-687.

    [21] 郭華明, 倪萍, 賈永鋒, 等. 內(nèi)蒙古河套盆地地表水-淺層地下水化學(xué)特征及成因[J]. 現(xiàn)代地質(zhì), 2015,29(2):229-237.

    [22] 吳初, 武雄, 錢程, 等. 內(nèi)蒙古杭錦旗氣田區(qū)地下水化學(xué)特征及其形成機(jī)制[J]. 現(xiàn)代地質(zhì),2017,31(3):629-636.

    [23] 錢程, 武雄. 鹽池內(nèi)流區(qū)地下水水化學(xué)特征及其形成作用[J]. 干旱區(qū)資源與環(huán)境, 2016, 30(3):169-175.

    [24] SAMANI S, MOGHADDAM A A. Hydrogeochemical characteristics and origin of salinity in Ajabshir aquifer, East Azerbaijan, Iran[J]. Quarterly Journal of Engineering Geology & Hydrogeology, 2015, 48(3):2014-2070.

    [25] WU Y, GIBSON C E. Mechanisms controlling the water chemistry of small lakes in Northern Ireland[J]. Water Research, 1996, 30(1): 178-182.

    [26] GIBBS R J. Mechanisms controlling world water chemistry[J]. Science, 1970, 170: 1088-1090.

    [27] 張志強(qiáng). 成都新機(jī)場(chǎng)地下水水化學(xué)特征及形成機(jī)制研究[D]. 成都:成都理工大學(xué), 2016: 39-41.

    [28] LI P, WU J, QIAN H. Assessment of groundwater quality for irrigation purposes and identification of hydrogeochemical evolution mechanisms in Pengyang County, China[J]. Environmental Earth Sciences, 2013, 69(7):2211-2225.

    [29] HOUATMIA F, AZOUZI R, CHAREF A, et al. Assessment of groundwater quality for irrigation and drinking purposes and identification of hydrogeochemical mechanisms evolution in Northeastern, Tunisia[J]. Environmental Earth Sciences, 2016, 75(9):746.

    [30] 張兆吉.華北平原地下水可持續(xù)利用圖集[M].北京:中國地圖出版社, 2009: 92-141.

    猜你喜歡
    特征研究
    抓住特征巧觀察
    FMS與YBT相關(guān)性的實(shí)證研究
    2020年國內(nèi)翻譯研究述評(píng)
    遼代千人邑研究述論
    新型冠狀病毒及其流行病學(xué)特征認(rèn)識(shí)
    視錯(cuò)覺在平面設(shè)計(jì)中的應(yīng)用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    如何表達(dá)“特征”
    不忠誠的四個(gè)特征
    EMA伺服控制系統(tǒng)研究
    抓住特征巧觀察
    成年动漫av网址| av在线老鸭窝| 精品久久久久久电影网| 街头女战士在线观看网站| 日韩av在线免费看完整版不卡| 国产福利在线免费观看视频| 波多野结衣一区麻豆| 久久这里只有精品19| 丰满乱子伦码专区| 王馨瑶露胸无遮挡在线观看| 亚洲av在线观看美女高潮| 国产伦理片在线播放av一区| 汤姆久久久久久久影院中文字幕| www.精华液| 一本大道久久a久久精品| 热re99久久精品国产66热6| 久久人人97超碰香蕉20202| 亚洲,欧美精品.| 国产福利在线免费观看视频| 高清av免费在线| 亚洲成人手机| 丝瓜视频免费看黄片| av视频免费观看在线观看| 中文字幕人妻丝袜制服| 亚洲综合色惰| av国产久精品久网站免费入址| 亚洲av综合色区一区| 久久这里有精品视频免费| 国产男女超爽视频在线观看| 久久这里只有精品19| 看免费av毛片| 久热久热在线精品观看| 少妇的逼水好多| 亚洲欧美精品自产自拍| 一区在线观看完整版| 国产成人免费观看mmmm| 人妻人人澡人人爽人人| 国产日韩欧美在线精品| 99国产精品免费福利视频| a级毛片在线看网站| 久久久久国产一级毛片高清牌| 亚洲国产欧美网| 亚洲欧美清纯卡通| 免费黄网站久久成人精品| 精品一区二区三卡| 亚洲成av片中文字幕在线观看 | 国产97色在线日韩免费| 国产免费视频播放在线视频| 久久精品国产亚洲av天美| 亚洲国产日韩一区二区| 久久久久精品久久久久真实原创| 91精品伊人久久大香线蕉| 18禁动态无遮挡网站| 高清视频免费观看一区二区| 欧美人与性动交α欧美精品济南到 | 一级片'在线观看视频| 精品少妇内射三级| 国产亚洲精品第一综合不卡| 永久网站在线| 欧美中文综合在线视频| 日韩一区二区三区影片| 91精品三级在线观看| 只有这里有精品99| 久久97久久精品| 中文字幕亚洲精品专区| 国产男女超爽视频在线观看| 在线观看三级黄色| 国产精品99久久99久久久不卡 | 男的添女的下面高潮视频| 亚洲综合色网址| 久久国内精品自在自线图片| 成年女人在线观看亚洲视频| 亚洲精品久久久久久婷婷小说| 国产精品免费大片| 一级a爱视频在线免费观看| 超碰97精品在线观看| 在线 av 中文字幕| xxxhd国产人妻xxx| 青春草亚洲视频在线观看| 国产精品欧美亚洲77777| 99久久中文字幕三级久久日本| 一区二区三区精品91| 日本免费在线观看一区| 久久国产精品大桥未久av| 国产成人精品福利久久| 高清视频免费观看一区二区| 不卡视频在线观看欧美| 高清视频免费观看一区二区| 亚洲成av片中文字幕在线观看 | 久久99蜜桃精品久久| 桃花免费在线播放| 不卡av一区二区三区| 天天影视国产精品| 男女啪啪激烈高潮av片| 免费观看性生交大片5| 免费看av在线观看网站| 国产高清不卡午夜福利| 亚洲综合色网址| 欧美日韩av久久| 久久婷婷青草| 日本色播在线视频| 一区二区三区乱码不卡18| 91成人精品电影| 欧美日韩一区二区视频在线观看视频在线| 97在线视频观看| a 毛片基地| 成人午夜精彩视频在线观看| 亚洲熟女精品中文字幕| 国产免费一区二区三区四区乱码| 久久婷婷青草| av免费观看日本| 久久国产亚洲av麻豆专区| 免费黄网站久久成人精品| 国产男女超爽视频在线观看| √禁漫天堂资源中文www| 国产精品免费视频内射| 日韩,欧美,国产一区二区三区| 欧美日韩精品成人综合77777| 中文字幕av电影在线播放| 天天躁夜夜躁狠狠躁躁| 人体艺术视频欧美日本| 国产精品欧美亚洲77777| 国产精品国产三级国产专区5o| 亚洲,欧美精品.| 69精品国产乱码久久久| 欧美人与性动交α欧美软件| 国产不卡av网站在线观看| 91国产中文字幕| 日日爽夜夜爽网站| 国产精品久久久av美女十八| 最新的欧美精品一区二区| 两个人看的免费小视频| 精品99又大又爽又粗少妇毛片| 成年动漫av网址| 91久久精品国产一区二区三区| 一级,二级,三级黄色视频| 叶爱在线成人免费视频播放| 人妻一区二区av| 99久久人妻综合| 乱人伦中国视频| 涩涩av久久男人的天堂| 中文天堂在线官网| 满18在线观看网站| 91午夜精品亚洲一区二区三区| 母亲3免费完整高清在线观看 | 免费播放大片免费观看视频在线观看| 99九九在线精品视频| 日本91视频免费播放| 欧美精品一区二区大全| 青青草视频在线视频观看| 亚洲精品国产av蜜桃| 色网站视频免费| 国产亚洲欧美精品永久| 免费在线观看黄色视频的| 丝袜美足系列| av视频免费观看在线观看| 制服诱惑二区| 久久韩国三级中文字幕| 最近的中文字幕免费完整| 宅男免费午夜| 亚洲一区二区三区欧美精品| 美女xxoo啪啪120秒动态图| 久久久久久久大尺度免费视频| 波野结衣二区三区在线| 男女国产视频网站| 亚洲精品在线美女| 亚洲av日韩在线播放| 在线观看国产h片| 下体分泌物呈黄色| 国产精品.久久久| 精品午夜福利在线看| 秋霞在线观看毛片| 伦理电影大哥的女人| 中文字幕色久视频| 黄色 视频免费看| 最近中文字幕2019免费版| 国产精品三级大全| 涩涩av久久男人的天堂| 亚洲精品视频女| 少妇猛男粗大的猛烈进出视频| av在线老鸭窝| 久久精品久久久久久噜噜老黄| 色婷婷久久久亚洲欧美| 亚洲中文av在线| 大码成人一级视频| 亚洲精品国产av蜜桃| 亚洲一码二码三码区别大吗| 久久青草综合色| 老女人水多毛片| 亚洲精品自拍成人| 国产成人免费无遮挡视频| 国产日韩欧美在线精品| 日韩不卡一区二区三区视频在线| 日韩一区二区三区影片| 亚洲av日韩在线播放| 校园人妻丝袜中文字幕| 在线观看人妻少妇| 99香蕉大伊视频| 高清欧美精品videossex| 久久久欧美国产精品| 国产国语露脸激情在线看| 在线 av 中文字幕| 老司机影院毛片| √禁漫天堂资源中文www| 免费黄网站久久成人精品| av网站免费在线观看视频| 啦啦啦中文免费视频观看日本| 久久久久久久亚洲中文字幕| 99香蕉大伊视频| 亚洲一区中文字幕在线| 国产男女超爽视频在线观看| 久久久久久久精品精品| 日韩欧美一区视频在线观看| 七月丁香在线播放| 999久久久国产精品视频| 午夜福利在线观看免费完整高清在| 国产视频首页在线观看| 美女午夜性视频免费| 欧美bdsm另类| 亚洲欧美中文字幕日韩二区| 色网站视频免费| 中国国产av一级| 男女午夜视频在线观看| 免费在线观看黄色视频的| 国产成人aa在线观看| 看十八女毛片水多多多| 国产成人av激情在线播放| 巨乳人妻的诱惑在线观看| av免费观看日本| 亚洲欧美中文字幕日韩二区| 黑丝袜美女国产一区| 亚洲国产精品一区三区| 亚洲精品久久成人aⅴ小说| 欧美bdsm另类| 另类精品久久| a级毛片在线看网站| 久久午夜福利片| 亚洲欧美成人精品一区二区| 一级片'在线观看视频| 18禁观看日本| 国产精品偷伦视频观看了| 人妻人人澡人人爽人人| videosex国产| 成人影院久久| 丝袜喷水一区| 免费观看a级毛片全部| 老司机影院成人| 免费大片黄手机在线观看| 色哟哟·www| 国产在线一区二区三区精| 一区二区三区精品91| 国产精品嫩草影院av在线观看| 久久鲁丝午夜福利片| 黑人猛操日本美女一级片| 精品人妻一区二区三区麻豆| 99香蕉大伊视频| 亚洲欧美精品综合一区二区三区 | 9色porny在线观看| 久久精品aⅴ一区二区三区四区 | 欧美成人午夜免费资源| 亚洲熟女精品中文字幕| 亚洲欧美精品综合一区二区三区 | 黄色毛片三级朝国网站| 在线精品无人区一区二区三| 在线观看国产h片| 免费少妇av软件| 亚洲美女视频黄频| 老汉色∧v一级毛片| 成人毛片a级毛片在线播放| 王馨瑶露胸无遮挡在线观看| 国产精品.久久久| 亚洲情色 制服丝袜| 日韩电影二区| 亚洲成国产人片在线观看| 丝袜在线中文字幕| 国产一区二区三区av在线| 9热在线视频观看99| 夫妻性生交免费视频一级片| 日日爽夜夜爽网站| 亚洲精品美女久久av网站| 国产精品欧美亚洲77777| 免费黄频网站在线观看国产| av片东京热男人的天堂| 亚洲美女搞黄在线观看| 精品国产乱码久久久久久小说| 一级a爱视频在线免费观看| 免费观看在线日韩| 日韩欧美精品免费久久| 成人亚洲精品一区在线观看| 在线观看美女被高潮喷水网站| 精品视频人人做人人爽| 爱豆传媒免费全集在线观看| 丁香六月天网| 又大又黄又爽视频免费| 18禁动态无遮挡网站| 汤姆久久久久久久影院中文字幕| 少妇人妻精品综合一区二区| 91午夜精品亚洲一区二区三区| av不卡在线播放| 性高湖久久久久久久久免费观看| 最近中文字幕高清免费大全6| 国产av一区二区精品久久| 各种免费的搞黄视频| 日韩制服丝袜自拍偷拍| 少妇人妻 视频| 性色av一级| 视频在线观看一区二区三区| 水蜜桃什么品种好| 精品人妻偷拍中文字幕| 狂野欧美激情性bbbbbb| 一级片免费观看大全| 亚洲欧美清纯卡通| 免费观看av网站的网址| 一本久久精品| 精品国产乱码久久久久久男人| av卡一久久| 亚洲久久久国产精品| 丰满乱子伦码专区| 久久久欧美国产精品| 日韩av免费高清视频| 秋霞伦理黄片| 欧美中文综合在线视频| 人妻 亚洲 视频| 国产精品人妻久久久影院| 国产在线一区二区三区精| 欧美bdsm另类| 久久人人97超碰香蕉20202| 2022亚洲国产成人精品| 免费黄色在线免费观看| 中文天堂在线官网| 一级a爱视频在线免费观看| 亚洲人成77777在线视频| 五月伊人婷婷丁香| 国产探花极品一区二区| 亚洲精品,欧美精品| 性色avwww在线观看| www.av在线官网国产| 国产精品免费大片| 日韩一区二区三区影片| 欧美激情极品国产一区二区三区| 青草久久国产| 最近中文字幕高清免费大全6| 亚洲伊人色综图| 国产女主播在线喷水免费视频网站| 国产精品 欧美亚洲| 欧美成人午夜精品| 老司机影院毛片| 又粗又硬又长又爽又黄的视频| 一二三四在线观看免费中文在| 蜜桃在线观看..| 日本vs欧美在线观看视频| av女优亚洲男人天堂| 亚洲男人天堂网一区| 国产精品免费大片| 亚洲国产精品国产精品| 久久精品人人爽人人爽视色| 国产视频首页在线观看| 日本91视频免费播放| 国产一区二区三区av在线| 一本久久精品| 丝袜喷水一区| 天天躁夜夜躁狠狠久久av| 久久精品久久精品一区二区三区| 亚洲图色成人| 男女免费视频国产| 亚洲欧洲精品一区二区精品久久久 | 亚洲欧美成人精品一区二区| 又粗又硬又长又爽又黄的视频| 美女视频免费永久观看网站| 午夜福利在线观看免费完整高清在| 日日摸夜夜添夜夜爱| 久久久a久久爽久久v久久| 精品福利永久在线观看| 午夜激情久久久久久久| 国产亚洲av片在线观看秒播厂| 国产精品 国内视频| videossex国产| 久久鲁丝午夜福利片| av天堂久久9| av线在线观看网站| 午夜av观看不卡| 久久久久网色| 国产白丝娇喘喷水9色精品| 亚洲欧洲日产国产| 999精品在线视频| 国产精品av久久久久免费| 国产亚洲最大av| 国产欧美日韩一区二区三区在线| 成年美女黄网站色视频大全免费| 观看美女的网站| 男女无遮挡免费网站观看| 亚洲欧洲国产日韩| 欧美成人午夜精品| 制服丝袜香蕉在线| 少妇猛男粗大的猛烈进出视频| 国产成人午夜福利电影在线观看| 极品人妻少妇av视频| 精品国产超薄肉色丝袜足j| 肉色欧美久久久久久久蜜桃| 久久97久久精品| 国产不卡av网站在线观看| 一本—道久久a久久精品蜜桃钙片| 成人二区视频| 交换朋友夫妻互换小说| 午夜老司机福利剧场| 99国产精品免费福利视频| 亚洲 欧美一区二区三区| 精品酒店卫生间| 王馨瑶露胸无遮挡在线观看| 久久精品国产亚洲av高清一级| 下体分泌物呈黄色| 91国产中文字幕| 香蕉精品网在线| 国产 一区精品| 女人被躁到高潮嗷嗷叫费观| 免费播放大片免费观看视频在线观看| 国产在线一区二区三区精| 日韩大片免费观看网站| 日韩欧美一区视频在线观看| 人人澡人人妻人| 精品酒店卫生间| 咕卡用的链子| 成人毛片60女人毛片免费| 一级a爱视频在线免费观看| 久久99精品国语久久久| 女人被躁到高潮嗷嗷叫费观| 女人精品久久久久毛片| 久久婷婷青草| 欧美日韩视频高清一区二区三区二| 欧美激情极品国产一区二区三区| 亚洲综合精品二区| 美女大奶头黄色视频| 韩国av在线不卡| av免费观看日本| 99久久人妻综合| 日韩三级伦理在线观看| 热99久久久久精品小说推荐| 一区二区三区四区激情视频| 美女脱内裤让男人舔精品视频| 女性被躁到高潮视频| av网站免费在线观看视频| 青春草国产在线视频| 免费观看a级毛片全部| 亚洲国产精品国产精品| 日韩av免费高清视频| 宅男免费午夜| 毛片一级片免费看久久久久| 在线看a的网站| 在线免费观看不下载黄p国产| 亚洲熟女精品中文字幕| videos熟女内射| 桃花免费在线播放| 欧美变态另类bdsm刘玥| 亚洲 欧美一区二区三区| 久久久久久久久久人人人人人人| 黑丝袜美女国产一区| 国产免费一区二区三区四区乱码| 国产女主播在线喷水免费视频网站| 99热全是精品| 亚洲av成人精品一二三区| 女人高潮潮喷娇喘18禁视频| 夫妻午夜视频| 国产欧美亚洲国产| 婷婷色综合www| 美女大奶头黄色视频| 亚洲成人一二三区av| 亚洲欧洲国产日韩| 美女中出高潮动态图| 精品第一国产精品| 国产成人精品在线电影| a 毛片基地| 亚洲色图 男人天堂 中文字幕| av线在线观看网站| 国产亚洲午夜精品一区二区久久| 国产精品二区激情视频| 亚洲综合精品二区| 国产极品天堂在线| 日韩制服丝袜自拍偷拍| www.精华液| 久久久久精品人妻al黑| 国产精品久久久久成人av| a级毛片在线看网站| 如何舔出高潮| 国产精品亚洲av一区麻豆 | 国产成人精品久久二区二区91 | 精品卡一卡二卡四卡免费| 大话2 男鬼变身卡| 激情视频va一区二区三区| 视频区图区小说| 午夜影院在线不卡| 天堂8中文在线网| 午夜福利在线免费观看网站| 国产黄色视频一区二区在线观看| 观看av在线不卡| 亚洲欧美色中文字幕在线| 国产亚洲最大av| 女性生殖器流出的白浆| 久久久久精品性色| 69精品国产乱码久久久| 国产免费又黄又爽又色| 国产精品女同一区二区软件| 久久国产精品大桥未久av| 欧美精品亚洲一区二区| 26uuu在线亚洲综合色| 国产xxxxx性猛交| 成人漫画全彩无遮挡| 国产高清不卡午夜福利| 亚洲欧美日韩另类电影网站| 国产国语露脸激情在线看| 国产淫语在线视频| av在线播放精品| 国产免费福利视频在线观看| 999精品在线视频| 一本色道久久久久久精品综合| 另类亚洲欧美激情| 国产成人a∨麻豆精品| 交换朋友夫妻互换小说| 亚洲精品国产av蜜桃| 欧美日韩亚洲国产一区二区在线观看 | av国产精品久久久久影院| 日本-黄色视频高清免费观看| 日韩一区二区三区影片| 亚洲av.av天堂| 1024香蕉在线观看| 咕卡用的链子| 久久99精品国语久久久| 久久久国产欧美日韩av| 久久影院123| 亚洲av电影在线进入| a 毛片基地| 久久久久精品久久久久真实原创| 久久这里有精品视频免费| 国产1区2区3区精品| 国产福利在线免费观看视频| 亚洲国产成人一精品久久久| 欧美中文综合在线视频| 亚洲经典国产精华液单| 侵犯人妻中文字幕一二三四区| 又大又黄又爽视频免费| 80岁老熟妇乱子伦牲交| 亚洲国产欧美网| 亚洲av国产av综合av卡| 国产亚洲精品第一综合不卡| 两个人看的免费小视频| 成年美女黄网站色视频大全免费| 婷婷成人精品国产| 自拍欧美九色日韩亚洲蝌蚪91| 中文字幕人妻丝袜一区二区 | 国产日韩欧美在线精品| 久久精品久久精品一区二区三区| 永久免费av网站大全| 久久久久久久国产电影| 欧美少妇被猛烈插入视频| 午夜福利在线免费观看网站| 人妻系列 视频| 久久久久久伊人网av| 精品人妻在线不人妻| 王馨瑶露胸无遮挡在线观看| 免费人妻精品一区二区三区视频| 精品亚洲成国产av| 黄色怎么调成土黄色| 制服丝袜香蕉在线| 三上悠亚av全集在线观看| 欧美另类一区| 国产男女超爽视频在线观看| 激情视频va一区二区三区| 国产一区亚洲一区在线观看| 免费大片黄手机在线观看| 亚洲欧美一区二区三区国产| 99热网站在线观看| 亚洲久久久国产精品| 美女国产高潮福利片在线看| 精品国产一区二区久久| 久久毛片免费看一区二区三区| 一级毛片我不卡| 午夜av观看不卡| 一区福利在线观看| av卡一久久| 大香蕉久久成人网| 欧美97在线视频| 久久 成人 亚洲| 丝袜人妻中文字幕| 久久久国产欧美日韩av| 啦啦啦中文免费视频观看日本| 一级片免费观看大全| 新久久久久国产一级毛片| 人人妻人人添人人爽欧美一区卜| 中文天堂在线官网| 亚洲少妇的诱惑av| 日本爱情动作片www.在线观看| 老鸭窝网址在线观看| 乱人伦中国视频| 国产精品女同一区二区软件| 纵有疾风起免费观看全集完整版| 熟妇人妻不卡中文字幕| 91久久精品国产一区二区三区| 狠狠婷婷综合久久久久久88av| 麻豆乱淫一区二区| 看非洲黑人一级黄片| 少妇精品久久久久久久| 亚洲精品自拍成人| 成人18禁高潮啪啪吃奶动态图| 看免费成人av毛片| av电影中文网址| 母亲3免费完整高清在线观看 | 两个人免费观看高清视频| 丰满少妇做爰视频| 婷婷色麻豆天堂久久| 国产又爽黄色视频| 亚洲婷婷狠狠爱综合网| 街头女战士在线观看网站| 亚洲精品久久成人aⅴ小说| 一区二区三区乱码不卡18| 日韩中字成人| 另类亚洲欧美激情|