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

    贛江中游淺層地下水水化學(xué)特征演化及成因分析

    2024-07-04 10:49:58胡健范超盧志剛
    人民長江 2024年6期
    關(guān)鍵詞:贛江水化學(xué)陽離子

    胡健 范超 盧志剛

    收稿日期:2023-11-27;接受日期:2024-01-26

    基金項目:國家自然科學(xué)基金地區(qū)項目(52364018);江西省地質(zhì)局第七地質(zhì)大隊井開區(qū)1∶10 000環(huán)境水文地質(zhì)調(diào)查項目

    作者簡介:

    胡? 健,男,講師,碩士,主要從事工程地質(zhì)水文地質(zhì)方面的研究。E-mail:974244132@qq.com

    Editorial Office of Yangtze River. This is an open access article under the CC BY-NC-ND 4.0 license.

    文章編號:1001-4179(2024) 06-0037-08

    引用本文:胡健,范超,盧志剛.

    贛江中游淺層地下水水化學(xué)特征演化及成因分析

    [J].人民長江,2024,55(6):37-44.

    摘要:贛江是江西省主要水源,研究其地下水水化學(xué)特征及成因機理對流域水質(zhì)保護至關(guān)重要。選取贛江中游吉安市井開區(qū)為研究區(qū),于2019年豐-平-枯水期各采集研究區(qū)淺層地下水樣45組,運用數(shù)理統(tǒng)計、Durov圖、Gibbs圖、離子比例系數(shù)、PHREEQC反向模擬等方法,研究了吉安市井開區(qū)淺層地下水的水化學(xué)分布特征和組成類型,討論了其主控因素、主成分離子來源,并開展了典型徑流路徑水文地球化學(xué)模擬。結(jié)果表明:研究區(qū)淺層地下水屬低礦化度淡水,主要陰陽離子為HCO-3、Ca2+,個別樣本中NH+4、NO-3、SO2-4的異常值均與人類活動有關(guān),水化學(xué)類型主要有HCO3-Ca型、HCO3-Ca·Mg型、HCO3·SO4-Ca型3種;研究區(qū)地下水化學(xué)形成的主控因素為巖石風(fēng)化作用;離子比值關(guān)系表明水中組分主要源于碳酸鹽巖礦物的溶解作用,其中Na+、K+主要源于鹽巖的溶解,SO2-4源于少量蒸發(fā)巖礦物(石膏)的溶解,Ca2+、Mg2+源于方解石、白云石和硅酸鹽礦物的溶解;水文地球化學(xué)反向模擬表明水流路徑均發(fā)生白云石、石膏、巖鹽的溶解作用,結(jié)合氯堿指數(shù),地下水流存在一定的逆向陽離子交換作用。研究成果可為贛江中游地區(qū)地下水開發(fā)利用與保護提供理論依據(jù)。

    關(guān)? 鍵? 詞:淺層地下水; 水化學(xué)特征; 離子來源; PHREEQC; 贛江中游

    中圖法分類號: P641.3

    文獻標志碼: A

    DOI:10.16232/j.cnki.1001-4179.2024.06.006

    0? 引 言

    地下水資源是人類生活和工業(yè)生產(chǎn)的重要資源,淺層地下水的利用最為廣泛,其有埋藏廣、易開采、水質(zhì)較好等優(yōu)點。然而隨著社會發(fā)展,用水量增加,人類活動影響地下水水質(zhì)變化,進一步威脅生態(tài)環(huán)境和人類健康[1-2]。贛江是長江的支流,是江西省內(nèi)最大河流,流域面積大,水系分支多,上中下游地區(qū)地下水與河流互相補給關(guān)系密切,地下水動態(tài)直接影響贛江流域水質(zhì)安全。吳通航等對鄱陽湖流域贛江下游地區(qū)進行了地下水化學(xué)特征分析,并開展了水質(zhì)及人類健康風(fēng)險評價,發(fā)現(xiàn)地下水對嬰兒存在嚴重潛在非致癌風(fēng)險[3]。劉小真等分析了贛江上、中、下游河段氮磷營養(yǎng)鹽濃度變化和分布特征,發(fā)現(xiàn)沿岸城市如豐城、吉安市等氮磷營養(yǎng)鹽濃度明顯高于其他點位,2011年中游的禾水、遂川江NH3-N和NO3-N濃度超標[4]。因此,研究贛江中游地區(qū)地下水化學(xué)特征、成因規(guī)律及演化機理意義重大。

    近年來,關(guān)于水文地球化學(xué)特征的研究持續(xù)進行,研究方法眾多,如:圖解法、離子比例系數(shù)法、空間分析、水文地球化學(xué)模擬等方法。王潤濤等運用Piper、Gibbs圖等多種方法對南淝河下游淺層地下水的水化學(xué)特征進行了研究,揭示了巖石風(fēng)化溶解為水化學(xué)主控作用,碳酸鹽礦物溶解有較大影響[5]。鄒嘉文等運用飽和指數(shù)、氯堿指數(shù)等方法研究了邯鄲市黑龍港平原潛水水化學(xué)特征及成因,發(fā)現(xiàn)水質(zhì)水位變化與南水北調(diào)中線工程通水密切相關(guān),蒸發(fā)巖和碳酸鹽巖礦物的溶解/沉淀作用為離子濃度變化的主因[6]。李捷等分析了北京市涼水河流域地下水和河水的氘氧關(guān)系,結(jié)合水化學(xué)特征,表明山前沖洪積扇河水入滲為研究區(qū)地下水主要補給源[7]。唐春雷等運用Durov圖、氫氧同位素、Sr同位素等方法分析了古堆泉域巖溶地下水化學(xué)特征和徑流特征,發(fā)現(xiàn)古堆泉口水樣變化原因是三泉水庫補給,而三泉水庫的變化是受引黃水的影響[8]。

    吉安市井開區(qū)位于贛江流域中游段,谷地平原地形,水資源豐富,城市化和工業(yè)加速發(fā)展,地下水資源不合理開采導(dǎo)致水質(zhì)受到一定影響。贛江流域吉安水文站數(shù)據(jù)指出,2008~2018年該河段水質(zhì)污染指標氨氮-總磷高度顯著上升,受人類活動的影響顯著[9]。而區(qū)域地下水水質(zhì)與地表水、大氣降水存在密切聯(lián)系,因此對其淺層地下水化學(xué)特征分析和成因展開研究工作尤為重要。本文以贛江中游吉安市井開區(qū)淺層地下水為研究對象,于2019年豐-平-枯水期各采集45組水樣品,查明研究區(qū)水文地質(zhì)條件,對水樣品數(shù)據(jù)進行整理分析,基于Durov圖、Gibbs圖、離子比例系數(shù)、PHREEQC反向模擬法等方法,開展了區(qū)域水化學(xué)特征、演化成因研究及其與環(huán)境相互作用機理分析討論工作,以期為今后贛江中游吉安市井開區(qū)地下水資源開發(fā)利用保護和水質(zhì)評價工作奠定理論基礎(chǔ)。

    1? 研究區(qū)概況

    研究區(qū)為贛江中游吉安市井開區(qū),位于吉安市東南部,面積為42.99 km2,規(guī)劃范圍西起105國道,東至繞城高速公路,北起廬陵大道以北200 m,吉州窯文化保護區(qū)南側(cè),南至南環(huán)路。區(qū)內(nèi)屬平原-低丘陵地貌,海拔由南向北逐漸降低。區(qū)內(nèi)水系以贛江為主流,禾水為支流,亞熱帶季風(fēng)濕潤氣候,四季分明,年均氣溫為17.5~18.6 ℃,年降水量為1 566.9 mm。區(qū)內(nèi)主要發(fā)育第四系地層,按其從老到新地層為石炭系梓山組下段灰、灰白色厚層狀含礫不等礫石石英砂巖、粗粒石英砂巖泥質(zhì)粉砂巖,二疊系茅口組硅質(zhì)巖、硅質(zhì)灰?guī)r,二疊系小江邊組鈣質(zhì)泥巖、鈣質(zhì)泥巖夾灰?guī)r透鏡體及薄層灰?guī)r,白堊系上統(tǒng)南雄組泥質(zhì)粉砂巖、粉砂巖,第四系上部為粉質(zhì)黏土,下部為礫質(zhì)黏土,中下部細砂及礫砂。區(qū)內(nèi)地下水類型主要有4種:第四系松散巖類孔隙水,富水性中等(單井涌水量142.34~687.69 m3/d)、弱富水(小于100 m3/d);紅層碎屑巖類孔隙裂隙水,水量貧乏(小于100 m3/d);碳酸巖類裂隙溶洞水,泉流量小于1 L/s;基巖裂隙水,泉流量小于0.1 L/s。區(qū)域地下水水位埋藏較淺,埋深1.0~6.8 m,以大氣降水、地表水補給為主,由南往北徑流,以泉、人工開采或向贛江、禾水河排泄。

    2? 樣品采集與分析方法

    2.1? 樣品采集與測定

    于2019年4月(豐水期)、9月(平水期)、11月(枯水期)開展本次研究的樣品采集工作,各取45組淺層地下水樣(民井和鉆井),井深8~25 m,采樣點位分布見圖1。取樣時,按規(guī)范操作需待取水樣潤洗3次,完成后,將容器瓶密封并貼好標簽。現(xiàn)場檢測水位和溫度,其余指標K+、Na+、Ca2+、Mg2+、Cl-、SO2-4、HCO-3、CO2-3、NO-3、TDS、pH等由江西省地質(zhì)礦產(chǎn)勘查開發(fā)局九一二實驗室完成。

    2.2? 數(shù)據(jù)分析方法

    采用SPSS 22.0軟件對研究區(qū)地下水檢測數(shù)據(jù)進行描述性統(tǒng)計分析,運用Origin 12b繪制Durov圖、Gibbs圖、離子比例系數(shù)圖?;贓XCEL實現(xiàn)舒卡列夫分類,分析研究區(qū)淺層地下水水化學(xué)類型和離子含量特征,進一步討論其水化學(xué)組分形成的主要水文地球化學(xué)作用及離子來源。運用PHREEQC軟件對選取的典型水流路徑進行反向水文地球化學(xué)模擬,分析水化學(xué)組分形成的演化過程。

    3? 結(jié)果與討論

    3.1? 研究區(qū)淺層地下水水化學(xué)特征

    3.1.1? 水化學(xué)組成特征

    對所采集的水樣檢測數(shù)據(jù)整理并統(tǒng)計結(jié)果,主要離子相關(guān)參數(shù)見表1。

    (1) 由表1可知,研究區(qū)淺層地下水pH范圍為5.33~8.21,平均值為7.15,總體呈中性—弱堿性,豐-平-枯水期數(shù)值略有變化。TDS值范圍為3~566 mg/L,平均值為203.17 mg/L,屬低礦化度水,豐平水期略高于枯水期。研究區(qū)主要陰陽離子濃度在豐-平-枯水期呈現(xiàn)差異,陽離子含量總體為Ca2+>Mg2+>Na+>K+,其中Ca2+范圍值為4.18~113.90 mg/L,豐-平-枯均值分別為20.30,27.83,49.04 mg/L;Mg2+豐-平-枯均值分別為5.00,7.36,6.68 mg/L;HCO-3為優(yōu)勢陰離子,范圍值為3.13~358.25 mg/L,

    豐-平-枯均值分別為151.75,159.07,133.22 mg/L;

    其次為SO2-4,范圍值為0~431.09 mg/L,豐-平-枯均值分別為23.51,14.26,26.81 mg/L。出現(xiàn)上述現(xiàn)象的原因是研究區(qū)下伏二疊系碳酸鹽類巖層的溶解使得Ca2+、Mg2+、HCO-3含量較高。

    (2) 研究區(qū)地下水水位埋深較淺,水化學(xué)組分必定受人為活動的干擾,少數(shù)水樣品檢測出較高數(shù)值的NH+4、NO-3、SO2-4。如,豐水期NO-3濃度范圍值為0~69.65 mg/L,均值達7.96 mg/L,主要有4處異常點位,推測與人類生產(chǎn)活動、生活垃圾不合理堆放有關(guān)。硫酸鹽礦物的溶解或硫化物的氧化是SO2-4的重要來源,點號DXS30枯水期SO2-4出現(xiàn)最高值431.09 mg/L,可能與周邊工廠廢棄物堆放存在較大影響,但豐水期、平水期值分別為98.92 mg/L和2.11 mg/L,這表明該區(qū)域地下水系統(tǒng)循環(huán)快、水質(zhì)不斷更新,變化幅度大。點號DXS32平水期NH+4出現(xiàn)最高值24.51 mg/L,但豐-枯水期均未檢出,說明豐水期受到一定的人為污染,而之后水質(zhì)得到更新恢復(fù)。

    3.1.2? 水化學(xué)類型特征

    如圖2所示,地下水樣集中分布在某一區(qū)域,陽離子主要靠近Ca軸,表明研究區(qū)淺層地下水陽離子以Ca2+、Mg2+為主,其毫克當量百分數(shù)平均值為74.07%和21.44%。陰離子數(shù)值點主要靠近三角形的HCO3端,說明陰離子主要為HCO-3,其毫克當量百分數(shù)平均值達78.08%。豐-平水期水樣分布特征較為相似,主要分布在Ca軸、三角形的HCO3端,反映出可能為碳酸鹽巖礦物的風(fēng)化溶解形成主要離子。枯水期水樣分布較豐平水期分散,主要分布在三角形的Ca端、HCO3端和陰離子中間過渡區(qū)域,說明人類活動使部分水化學(xué)類型由HCO3-Ca型轉(zhuǎn)為HCO3·Cl-Ca型、HCO3·SO4-Ca型。

    從中間的正方形比較Ca2+、Na++K+和HCO-3、Cl-的百分數(shù)含量,絕大多數(shù)水樣均落在右上部區(qū)域(HCO-3>Cl-,Ca2+>Na++K+),包括豐-平水期所有水樣,枯水期3/4的水樣,說明了研究區(qū)淺層地下水主導(dǎo)離子為HCO-3和Ca2+??菟趨^(qū)域外水樣均位于左上部區(qū)域(Cl->HCO-3,Ca2+>Na++K+),說明存在一定程度人為污染導(dǎo)致Cl-含量增加。

    從Cl-、HCO-3的相對含量與pH關(guān)系分析,Cl-相對含量增加,pH明顯有降低趨勢,HCO-3含量增大,則pH增大,枯水期水樣具有最小的pH。從Ca2+、Na++K+與TDS的關(guān)系分析,Ca2+相對含量增大,TDS呈下降趨勢,表明影響研究區(qū)地下水TDS變化的主要離子為陰離子。

    對水中離子進行舒卡列夫分類,結(jié)果見表2,研究區(qū)淺層地下水一共有10種水化學(xué)類型,主要有HCO3-Ca型、HCO3-Ca·Mg型、HCO3·SO4-Ca型3類,其占比分別為55.56%,23.70%,8.15%,其余水化學(xué)類型均低于5%。

    3.2? 研究區(qū)淺層地下水化學(xué)成因分析

    3.2.1? 水化學(xué)組分主要控制因素

    地下水水化學(xué)組分隨時間和空間的改變而動態(tài)變化,導(dǎo)致水化學(xué)特征各異和成因復(fù)雜。Gibbs圖是根據(jù)c(Na+)/c(Na++Ca2+)和c(Cl-)/c(Cl-+HCO-3)的離子濃度比與TDS之間的比值關(guān)系繪制散點圖,同時劃分了蒸發(fā)濃縮、巖石風(fēng)化和降雨控制3個區(qū)域,從點位的分布位置判斷樣本點水化學(xué)特征形成的主要機制[10-12]。研究區(qū)豐-平-枯水期地下水中c(Na+)/c(Na++Ca2+)和c(Cl-)/c(Cl-+HCO-3)與TDS之間的比值關(guān)系Gibbs圖見圖3。由圖3可知:研究區(qū)豐-平-枯水期淺層地下水中離子含量c(Na+)/c(Na++Ca2+)比值均小于0.5,所有水樣品數(shù)據(jù)點集中分布在圖的中部左側(cè),表明地下水離子主要受到巖石風(fēng)化作用的控制,而蒸發(fā)濃縮和降雨作用微弱。然而c(Cl-)/c(Cl-+HCO-3)比值大部分小于0.5,分布集中,僅枯水期少部分水樣比值大于0.5,TDS中等偏低,位于圖右下方,表明局部區(qū)域地下水中離子受到降水作用的影響,也可能是該時段人類活動造成Cl-離子的輸入。

    通常運用M(Mg2+)/M(Na+)、M(HCO-3)/M(Na+)和M(Ca2+)/M(Na+)的摩爾濃度比值來分析離子來源于碳酸鹽巖、硅酸鹽巖和蒸發(fā)巖的風(fēng)化溶解[13]。研究區(qū)豐-平-枯水期地下水中M(Mg2+)/M(Na+) 、M(HCO-3)/M(Na+) 和M(Ca2+)/M(Na+)關(guān)系(不包括Na+未檢出的樣品數(shù)據(jù))比值圖見圖4。據(jù)圖4所示,研究區(qū)水樣品數(shù)據(jù)點大部分靠近在碳酸鹽礦物的端元,說明該區(qū)域地下水離子組分主要由碳酸鹽巖礦物風(fēng)化溶解作用控制,與下伏含水層巖性為二疊系硅質(zhì)巖、硅質(zhì)灰?guī)r相對應(yīng)??菟诓糠炙畼悠窋?shù)據(jù)投影點靠近硅酸鹽巖礦物端元,也說明研究區(qū)地下水離子也受到硅酸鹽礦物風(fēng)化作用的影響。

    3.2.2? 水化學(xué)組分主要離子來源

    為進一步分析研究區(qū)淺層地下水水化學(xué)成因及變化規(guī)律,運用離子比例系數(shù)法分析地下水的演化規(guī)律、識別水化學(xué)組分來源,即通過比較不同離子毫克當量濃度(單位為mEq,用符號γ表示)間的關(guān)系,判斷離子主要來源。

    (1) 利用γ(Na++K+)/γ(Cl-)比值來判定地下水中Na+、K+的來源。若γ(Na++K+)/γ(Cl-)接近于1,說明主要來自于鹽巖溶解;若大于1,則說明Na+、K+除鹽巖外還來自于硅酸鹽巖溶解,或者存在陽離子交換作用[14-16]。如圖5(a)所示,樣品點集中于左下角,濃度較低,部分樣品點位于1∶1線附近,表明地下水中Na+、K+主要來源于鹽巖溶解,部分數(shù)據(jù)投影點偏離1∶1線并處于線下方,表明人為污染改變了Cl-離子含量。

    (2) 利用γ(Ca2++Mg2+)/γ(HCO-3)和γ(Ca2++Mg2+)/γ(HCO-3+SO2-4)比值來判定地下水中Ca2+、Mg2+、HCO-3、SO2-4的來源。若γ(Ca2++Mg2+)/γ(HCO-3+SO2-4)比值大于1,表明存在硅酸鹽巖的溶解;而小于1,反映出水化學(xué)過程伴隨了離子交換吸附作用[17]。研究區(qū)水樣點總體分布在1∶1線附近,說明既有碳酸鹽礦物的溶解也有硫酸鹽礦物(石膏)的溶解,尤其是枯水期水樣點最靠近1∶1;而豐-平水期水樣點大部分位于1∶1線下方,則表明此時段研究區(qū)淺層地下水中的Ca2+、Mg2+離子來自碳酸鹽巖、硫酸鹽巖礦物的溶解,且同時陽離子交替吸附作用不可忽略,或者存在其他陽離子與HCO-3、SO2-4相平衡。γ(Ca2++Mg2+)/γ(HCO-3)比值圖與其圖5(c)相似,而隨HCO-3含量的增加,水樣品數(shù)據(jù)點分布整體輕微向上移動,表明石膏的溶解產(chǎn)生的SO2-4、Ca2+較少。

    (3) 利用γ(Ca2+)/γ(Mg2+)來判定主要溶解的碳酸鹽礦物是白云石或方解石。若比值為1,表明為白云石的溶解,隨著數(shù)值的增大,方解石的溶解增加;而比值大于2,表明硅酸鹽礦物的溶解對地下水中Ca2+、Mg2+ 來源有顯著貢獻[18]。如圖5(d)所示,研究區(qū)水樣品數(shù)據(jù)點多數(shù)分布于2:1線外,尤其枯水期,表明水化學(xué)離子組分Ca2+、Mg2+不僅源于白云石、方解石的溶解,硅酸鹽巖礦物的溶解也有較大貢獻。

    (4) 利用γ(Cl-+SO2-4)/γ(HCO-3)比值來判定水化學(xué)組分主要蒸發(fā)巖或碳酸鹽巖礦物的溶解。若比值大于1,表明蒸發(fā)巖的溶解作用大;若小于1,則表明碳酸鹽巖的溶解作用大[19-20]。如圖5(e)所示,研究區(qū)豐-平-枯水期水樣點絕大多數(shù)位于1∶1線的下方,意味著碳酸鹽巖礦物的溶解作用對水中各離子的形成有較大影響。

    (5) 利用γ(Ca2++Mg2+-HCO-3-SO2-4)/γ(Na++K+-Cl-)比值來判定水化學(xué)組分受到陽離子交替吸附作用的影響大?。?1-22]。如圖5(f)所示,研究區(qū)枯-平水期部分水樣品數(shù)據(jù)點分布于-1∶1線上,表明該時期地下水中存在陽離子交換作用,而豐水期水樣品數(shù)據(jù)點的分布偏離-1∶1線,表明該時期陽離子交換作用不明顯。

    (6) 利用氯-堿指數(shù)CAI-Ⅰ和CAI-Ⅱ判定陽離子交換作用的程度和方向。若CAI-Ⅰ、CAI-Ⅱ均為負值,表示巖土表面的Na+、K+代換了地下水中的Ca2+、Mg2+;若均為正值,則表示Ca2+、Mg2+交換了

    Na+、K+[23-24]。如圖6所示,研究區(qū)枯-平水期大部

    水樣點CAI-Ⅰ和CAI-Ⅱ為正值,占樣品總數(shù)的75%,說明含水層礦物表面的Ca2+、Mg2+取代了地下

    水中Na+、K+,即存在逆向陽離子交換作用。研究區(qū)

    淺層地下水水流交替迅速,含有豐富的二氧化碳,加速了

    溶濾作用,Ca2+、Mg2+逐漸變成主要的陽離子。

    3.3? 反向水文地球化學(xué)模擬

    PHREEQC是用于模擬水文地球化學(xué)反應(yīng)的軟件,計算礦物飽和指數(shù)(SI)和電荷平衡系數(shù)為前置工

    作[25]。若礦物飽和指數(shù)SI>0,則表明該礦物在水中

    達到飽和;若SI<0,則表明該礦物在水中不飽和[26]。結(jié)合水文地質(zhì)條件、水文地質(zhì)調(diào)查,選取兩個剖面的兩條水流路徑進行模擬,分別為路徑1(DXS14→DXS06)、路徑2(DXS39→DXS35),根據(jù)路徑水樣點化學(xué)組分(表3),從上游到下游地下水各化學(xué)組分隨水流含量升高,路徑2總體略高于路徑1。計算各模擬路徑起始點電荷平衡系數(shù)和礦物飽和指數(shù)如表4所示,各模擬點電荷平衡系數(shù)幾乎為0,表明溶液基本處于電荷平衡狀態(tài);路徑1中硬石膏、文石、方解石、白云石、石膏、巖鹽、鉀鹽處于不飽和狀態(tài);而路徑2中文石、方解石、白云石處于飽和及平衡狀態(tài),表明地下水在流動過程中存在礦物的溶解或沉淀。

    根據(jù)模擬路徑點的礦物飽和指數(shù),結(jié)合淺層地下水的賦存條件,第四紀地層、紅層碎屑巖及碳酸鹽地層中的常見礦物,結(jié)合前人研究成果[27-28],選取白云石、方解石、石膏、巖鹽、鉀鹽、鉀長石、鈉長石、鈣長石、CO2(g)以及Na+、Mg2+、Ca2+作為可能礦物相。運用

    PHREEQC建立反向水文地球化學(xué)模型,調(diào)試不確定系數(shù)為0.05,運行模擬出路徑1、2發(fā)生水巖相互作用

    的礦物轉(zhuǎn)移質(zhì)量(表5)。結(jié)果表明:路徑1中白云石、石膏、巖鹽、鉀長石、CO2(g)的溶解,溶解量最大為巖鹽的1.02 mmol/L,方解石發(fā)生沉淀,遷出量為0.611 mmol/L,水巖相互作用較弱;路徑2中存在白云石、方

    解石、石膏、巖鹽、CO2(g)的溶解,少量鉀長石遷出,摩

    爾轉(zhuǎn)移量為0.027 mmol/L,鈉長石的溶解、沉淀與鉀

    長石的狀態(tài)高度相關(guān)。兩條路徑均發(fā)生巖石中的

    Ca2+置換地下水中的Na+的逆向陽離子交換作用,與前文氯堿指數(shù)分析得出結(jié)論相同,物質(zhì)轉(zhuǎn)移量與溶解沉淀量相當,表明陽離子交換作用為重要的水文地球化學(xué)作用過程。

    4? 結(jié) 論

    (1) 研究區(qū)淺層地下水pH均值為7.15,總體呈中性—弱堿性,TDS均值為203.17 mg/L,屬低礦化度水,HCO-3、Ca2+為優(yōu)勢陰陽離子,少數(shù)水樣品檢測出較高數(shù)值的NH+4、SO2-4,豐-平-枯水期差異較大,表明區(qū)域地下水循環(huán)更新較快,恢復(fù)較快。

    (2) 根據(jù)Durov圖分析,地下水陰陽離子以HCO-3、Ca2+、Mg2+為主,可能為碳酸鹽巖礦物的風(fēng)化溶解作用形成;從Ca2+、Na++K+與TDS關(guān)系分析,影響研究區(qū)地下水TDS變化的主要離子為陰離子;研究區(qū)淺層水化學(xué)類型有10種,主要為HCO3-Ca型,HCO3-Ca·Mg型,HCO3·SO4-Ca型3種。

    (3) 根據(jù)Gibbs圖、主要離子濃度比值,表明研究區(qū)淺層地下水水文地球化學(xué)作用為巖石風(fēng)化作用。碳酸鹽巖礦物的溶解對研究區(qū)淺層地下水水化學(xué)組分的形成起決定作用,同時硅酸鹽礦物也存在一定影響。其中Na+、K+主要來自鹽巖的溶解,SO2-4源于少量蒸發(fā)巖礦物的溶解,據(jù)75%水樣品氯堿指數(shù)均為正值,表明發(fā)生了逆向陽離子交換作用。

    (4) 根據(jù)兩條水流路徑的反向水文地球化學(xué)模擬結(jié)果,表明白云石、石膏、巖鹽發(fā)生溶解作用,方解石狀態(tài)存在差異,存在一定的逆向陽離子交換作用影響Ca2+、Na+的濃度。

    綜上所述,研究區(qū)淺層地下水總體水質(zhì)較好,水化學(xué)特征受氣候、巖性、人類活動等影響,不同時期略有差異。地下水補徑排條件好,水化學(xué)組分主要受巖石風(fēng)化溶解作用控制,各離子來源不一,水巖作用中各礦物溶解沉淀交互,陽離子交換作用不可忽視。

    參考文獻:[1]? 黃旭娟.環(huán)鄱陽湖區(qū)淺層地下水化學(xué)和同位素特征及成因分析[D].武漢:中國地質(zhì)大學(xué),2017.

    [2]? 王慧瑋,郭小嬌,張千千,等.濾沱河流域地下水水化學(xué)特征演化及成因分析[J].環(huán)境化學(xué),2021,40(12):3838-3845.

    [3]? 吳通航,劉海燕,張衛(wèi)民,等.鄱陽湖流域贛江下游水化學(xué)特征及人類健康風(fēng)險評價[J].現(xiàn)代地質(zhì),2022,36(2):427-438.

    [4]? 劉小真,楊子依,曹炳偉,等.贛江主要河流豐水期-枯水期氮磷營養(yǎng)鹽分布特征[J].中國環(huán)境監(jiān)測,2023,39(2):21-32.

    [5]? 王潤濤,錢家忠,王宿軒,等.南淝河下游淺層地下水水化學(xué)特征及成因分析[J/OL].現(xiàn)代地質(zhì),1-10[2024-04-08].

    https:∥doi.org/10.19657/j.geoscience.1000-8527.2022.028.

    [6]? 鄒嘉文,劉飛,張靖坤.南水北調(diào)典型受水區(qū)淺層地下水水化學(xué)特征及成因[J].中國環(huán)境科學(xué),2022,42(5):2260-2268.

    [7]? 李捷,姜穎,劉玉蓮,等.涼水河流域地下水水化學(xué)特征和時空變化規(guī)律[J].中國環(huán)境科學(xué),2022,42(4):1847-1853.

    [8]? 唐春雷,申豪勇,趙春紅,等.古堆泉域巖溶地下水水化學(xué)特征及成因[J].環(huán)境科學(xué),2023,44(9):4874-4883.

    [9]? 候林麗,鄒武,徐鵬,等.贛江吉安河段水質(zhì)變化及其趨勢分析[J].人民長江,2020,51(11):58-63.

    [10]GIBBS R J.Mechanisms controlling world water chemistry[J].Science,1970,170(3962):1088-1090.

    [11]LIANG H T,WANG W,LI J L et al.Hydrochemical characteristics and health risk assessment of groundwater in Dingbian county of the Chinese Loess Plateau,northwest China[J].Environmental Earth Sciences,2022,81(11):1-16.

    [12]XIAO Y,ZHANG J H,LONG A H,et al.Hydrochemical characteristics and formation mechanism of quaternary groundwater in Baoshan Basin,Western Yunnan,China[J].Water,2023,15(15):2736.

    [13]吳璽,安永會,魏世博,等.黑河下游鼎新谷地淺層地下水水化學(xué)特征及演化規(guī)律[J].干旱區(qū)資源與環(huán)境,2021,35(9):103-109.

    [14]鄭慧銘,張豐,史啟朋,等.煤炭資源型城市淺層地下水化學(xué)特征及成因分析[J].人民長江,2024,55(2):51-56.

    [15]吳銘杰,李勝偉,王東輝,等.岷江流域含石膏地層地下水水化學(xué)差異性研究[J].人民長江,2023,54(7):14-21.

    [16]蔡小虎,康叢軒.錢塘江下游濱江地區(qū)地下水水化學(xué)特征與水質(zhì)評價[J].人民長江,2023,54(6):27-33.

    [17]劉松韜,張東,李玉紅,等.伊洛河流域河水來源及水化學(xué)組成控制因素[J].環(huán)境科學(xué),2020,41(3):1184-1196.

    [18]李瀟瀚,張翼龍,王瑞,等.呼和浩特盆地地下水化學(xué)特征及成因[J].南水北調(diào)與水利科技,2018,16(4):136-145.

    [19]ZHANG L Z,LIANG X J,XIAO C L et al.Hydrochemical characteristics and the impact of human activities on groundwater in a semi-arid plain:a case study of western Jilin Province,Northeast China[J].Environmental Science and Pollution Research,2023,30(2):110204.

    [20]李會亞,馮起,陳朋娟,等.民勤綠洲灌區(qū)地下水水化學(xué)特征及其演化驅(qū)動機理[J].干旱區(qū)研究,2017,34(4):733-740.

    [21]ZHANG J W,LIANG X,JIN M G,et al.Identifying the groundwater flow systems in a condensed rivet-network interlluve between the Han River and Yangtze River (China) using hydrogeochemi-cal indicators[J].Hydrogeology Journal,2019,27(7):2415-2430.

    [22]ZHAO Q,LUO Z J,DU J G et al.Hydrogeochemical characteristics and genesis of geothermal water in Xiaoyangkou,Nantong,China[J].Environmental Earth Sciences,2023,82:456.

    [23]吳嘉鈴,王瑩,胡倩,等.雷州半島地下水水化學(xué)特征及成因分析[J].安全與環(huán)境工程,2022,29(1):145-154.

    [24]胡昱欣,周瑞靜,宋煒,等.北京市豐臺區(qū)永定河以東淺層地下水水化學(xué)演變規(guī)律及成因分析[J].環(huán)境科學(xué),2024,45(5):2651-2664.

    [25]萬朔陽,吳勇,李仁海,等.什邡市城市規(guī)劃區(qū)地下水化學(xué)特征及演化模擬[J].科學(xué)技術(shù)與工程,2019,19(26):94-100.

    [26]趙江濤,周金龍,梁川,等.新疆焉耆盆地平原區(qū)地下水反向水文地球化學(xué)模擬[J].干旱區(qū)資源與環(huán)境,2017,31(10):65-70.

    [27]余東,周金龍,張杰,等.新疆喀什地區(qū)地下水鐵錳水文地球化學(xué)及演化規(guī)律[J].環(huán)境科學(xué)學(xué)報,2021,41(6):2169-2181.

    [28]陳盟,吳勇,高東東,等.廣漢市平原區(qū)淺層地下水化學(xué)演化及其控制因素[J].吉林大學(xué)學(xué)報(地球科學(xué)版),2016,46(3):831-843.

    (編輯:劉 媛)

    Evolution of groundwater hydrochemical characteristics and origin analysis in middle reaches of Ganjiang River

    HU Jian1,2,F(xiàn)AN Chao3,LU Zhigang1,2

    (1.School of Resources and Civil Engineering,Gannan University of Science and Technology,Ganzhou 341000,China;

    2.Ganzhou Key Laboratory of Mine Geological Disaster Prevention and Control and Ecological Restoration,Gannan University of Science and Technology,Ganzhou 341000,China;

    3.The Seventh Geological Brigade of Jiangxi Bureau of Geology,Ganzhou 341000,China)

    Abstract:

    Ganjiang River is the main water source of Jiangxi Province,studying its groundwater hydrochemical characteristics and genesis mechanism is crucial for the protection of the water quality in the basin.This study selected the Jingkai area of Ji′an City in the middle reaches of the Ganjiang River as the research area and collected 45 sets of shallow groundwater samples during the flood season,normal season,and dry season in 2019.Then,methods such as mathematical statistics,Durov diagrams,Gibbs diagrams,ion ratio coefficients,and PHREEQC reverse simulation were used to study the hydrochemical distribution characteristics and composition types of shallow groundwater in the region.The main controlling factors and ion sources of the main components were discussed,and hydrogeochemical simulations of typical runoff paths were carried out.The results showed that:the shallow groundwater in the study area belongs to low salinity freshwater,with the main anions and cations being HCO-3 and Ca2+.The abnormal values of NH+4,NO-3,and SO2-4 in some samples are related to human activities.The main chemical types of groundwater in the study area are HCO3-Ca,HCO3-Ca·Mg,and HCO3·SO4-Ca.The main controlling factor of the formation of groundwater chemistry in the study area is rock weathering.From the ion ratio relationship,it can be concluded that the main source of water components is the dissolution of carbonate minerals,among them,Na+ and K+ mainly come from the dissolution of salt rocks,SO2-4 comes from the dissolution of a small amount of evaporite minerals (gypsum),Ca2+ and Mg2+ come from the dissolution of calcite,dolomite,and silicate minerals.The reverse simulation results of hydrogeochemistry indicated that the dissolution of dolomite,gypsum,and rock salt occurs along the water flow path.Combined with the chlor alkali index,there is a certain degree of reverse cation exchange in groundwater flow.The results can provide some theoretical basis for the development,utilization,and protection of groundwater in the middle reaches of the Ganjiang River.

    Key words:

    shallow groundwater; hydrochemical characteristics; ion source; PHREEQC; middle reaches of Ganjiang River

    猜你喜歡
    贛江水化學(xué)陽離子
    幸福湖泉域地下水水化學(xué)特征及其演化機制
    人民黃河(2025年1期)2025-01-16 00:00:00
    千里贛江第一城
    心聲歌刊(2023年5期)2023-10-10 13:31:32
    贛江奔流
    心聲歌刊(2023年3期)2023-07-16 02:11:52
    秦皇島北部地表水和地下水水化學(xué)特征及水質(zhì)評價
    贛江頌
    心聲歌刊(2019年5期)2020-01-19 01:52:52
    烷基胺插層蒙脫土的陽離子交換容量研究
    雨落贛江 〔外一首〕
    中國詩歌(2017年12期)2017-11-15 04:04:05
    衡水市新近系館陶組地?zé)崴瘜W(xué)特征與分布規(guī)律研究
    筠連煤層氣井產(chǎn)出水化學(xué)特征及意義
    中國煤層氣(2015年5期)2015-08-22 03:26:32
    陽離子Gemini表面活性劑的應(yīng)用研究進展
    大又大粗又爽又黄少妇毛片口| 99久久精品热视频| 亚洲电影在线观看av| 国产高潮美女av| av在线观看视频网站免费| 精品国产乱码久久久久久小说| 久久久久久久久大av| 久久久精品94久久精品| 欧美成人一区二区免费高清观看| 亚洲av综合色区一区| 精品亚洲乱码少妇综合久久| 九九爱精品视频在线观看| 久久精品国产鲁丝片午夜精品| 亚洲精品456在线播放app| 国产精品久久久久久久久免| 久久久精品免费免费高清| 成人高潮视频无遮挡免费网站| 一级毛片黄色毛片免费观看视频| 国产一区二区三区av在线| 亚洲一级一片aⅴ在线观看| 国产视频首页在线观看| 亚洲av免费高清在线观看| 插阴视频在线观看视频| 我要看黄色一级片免费的| 久久精品夜色国产| 国产乱来视频区| 精品一区二区三卡| 九草在线视频观看| 两个人的视频大全免费| 日韩伦理黄色片| 色婷婷av一区二区三区视频| 国产精品人妻久久久影院| 一级毛片久久久久久久久女| 欧美精品亚洲一区二区| 国产精品久久久久久精品古装| 国产精品伦人一区二区| 精品亚洲乱码少妇综合久久| 国产精品.久久久| 久久久久久九九精品二区国产| 日韩伦理黄色片| 精品久久国产蜜桃| 国模一区二区三区四区视频| 韩国av在线不卡| 精品亚洲乱码少妇综合久久| 国产白丝娇喘喷水9色精品| 亚洲av国产av综合av卡| 亚洲伊人久久精品综合| 搡老乐熟女国产| 久久综合国产亚洲精品| 美女xxoo啪啪120秒动态图| 久久青草综合色| av女优亚洲男人天堂| 久久精品国产亚洲av天美| 欧美老熟妇乱子伦牲交| 国产av精品麻豆| 街头女战士在线观看网站| 日日撸夜夜添| av一本久久久久| 激情五月婷婷亚洲| 97精品久久久久久久久久精品| 晚上一个人看的免费电影| 欧美丝袜亚洲另类| 亚洲精品456在线播放app| 青春草视频在线免费观看| 久久女婷五月综合色啪小说| 亚洲经典国产精华液单| 乱码一卡2卡4卡精品| 日日撸夜夜添| 亚洲av日韩在线播放| 精品一品国产午夜福利视频| 精品久久久噜噜| 国产高清三级在线| 在线观看人妻少妇| 高清在线视频一区二区三区| 男人狂女人下面高潮的视频| 亚洲欧美日韩另类电影网站 | 亚洲欧美成人综合另类久久久| 欧美日韩在线观看h| 精品人妻偷拍中文字幕| 亚洲精品aⅴ在线观看| 人体艺术视频欧美日本| 99精国产麻豆久久婷婷| 美女脱内裤让男人舔精品视频| 边亲边吃奶的免费视频| 又粗又硬又长又爽又黄的视频| 久久6这里有精品| 午夜激情福利司机影院| 日韩一区二区视频免费看| 夫妻性生交免费视频一级片| 亚洲av日韩在线播放| 久久99精品国语久久久| 亚洲真实伦在线观看| 日本av免费视频播放| 久久久亚洲精品成人影院| 久久久久久九九精品二区国产| 国产 一区精品| 久久热精品热| 两个人的视频大全免费| 亚洲av中文av极速乱| 观看av在线不卡| 色婷婷av一区二区三区视频| 一级片'在线观看视频| 成人国产av品久久久| 久久精品夜色国产| av不卡在线播放| 大香蕉97超碰在线| 一区二区三区免费毛片| 性色avwww在线观看| 国产精品久久久久成人av| 99久国产av精品国产电影| 国产精品久久久久成人av| 亚洲精品国产av成人精品| 国产精品嫩草影院av在线观看| 久久久色成人| 免费在线观看成人毛片| 又黄又爽又刺激的免费视频.| 国产精品一区二区在线观看99| 精华霜和精华液先用哪个| 一级二级三级毛片免费看| 精品人妻熟女av久视频| 久久久精品免费免费高清| 亚洲欧美一区二区三区黑人 | 成人亚洲精品一区在线观看 | 亚洲精品国产av成人精品| 亚洲第一av免费看| 国产亚洲5aaaaa淫片| 国产精品一区二区在线不卡| 中文字幕免费在线视频6| 成人亚洲欧美一区二区av| 亚州av有码| 女性被躁到高潮视频| 人妻 亚洲 视频| 婷婷色麻豆天堂久久| 日韩av免费高清视频| 99久国产av精品国产电影| 新久久久久国产一级毛片| 久久这里有精品视频免费| 亚洲在久久综合| 街头女战士在线观看网站| 欧美国产精品一级二级三级 | 国产成人aa在线观看| 在线观看一区二区三区激情| 肉色欧美久久久久久久蜜桃| 久久久久久伊人网av| 舔av片在线| 午夜激情福利司机影院| av专区在线播放| 久久精品国产亚洲网站| 中国美白少妇内射xxxbb| 这个男人来自地球电影免费观看 | h视频一区二区三区| 国产高清不卡午夜福利| 国产午夜精品久久久久久一区二区三区| 久久久精品94久久精品| 欧美3d第一页| 国产欧美亚洲国产| 最黄视频免费看| 日韩成人av中文字幕在线观看| 青春草国产在线视频| 赤兔流量卡办理| 亚洲精品乱码久久久久久按摩| 国产亚洲午夜精品一区二区久久| 国产人妻一区二区三区在| 夜夜爽夜夜爽视频| 纵有疾风起免费观看全集完整版| 国产高清国产精品国产三级 | av在线app专区| 一个人看的www免费观看视频| 国内揄拍国产精品人妻在线| 99九九线精品视频在线观看视频| a级毛色黄片| 国产成人一区二区在线| 美女主播在线视频| 嫩草影院新地址| 嘟嘟电影网在线观看| 欧美激情极品国产一区二区三区 | 国产精品福利在线免费观看| 精品少妇久久久久久888优播| 国产精品久久久久久精品古装| 国产成人freesex在线| 亚洲欧美日韩无卡精品| 免费大片18禁| 欧美xxxx性猛交bbbb| 成人一区二区视频在线观看| 男的添女的下面高潮视频| 在线观看三级黄色| 亚洲人与动物交配视频| 久久亚洲国产成人精品v| 欧美亚洲 丝袜 人妻 在线| 男女边吃奶边做爰视频| 美女高潮的动态| 在线播放无遮挡| 男女边摸边吃奶| 成人特级av手机在线观看| 51国产日韩欧美| h视频一区二区三区| 在线观看av片永久免费下载| 黑人猛操日本美女一级片| 国产成人免费无遮挡视频| 最近的中文字幕免费完整| 国产亚洲一区二区精品| 99久久人妻综合| a级一级毛片免费在线观看| 九草在线视频观看| 99视频精品全部免费 在线| av.在线天堂| 一区二区三区精品91| 观看av在线不卡| 欧美精品国产亚洲| 激情 狠狠 欧美| 日韩欧美精品免费久久| 三级国产精品片| 少妇精品久久久久久久| 91精品一卡2卡3卡4卡| 欧美区成人在线视频| 日本色播在线视频| 18禁在线播放成人免费| 国产乱来视频区| 国产毛片在线视频| av在线观看视频网站免费| 天美传媒精品一区二区| 久久久精品免费免费高清| 少妇丰满av| 中文乱码字字幕精品一区二区三区| 一个人看视频在线观看www免费| 99久久中文字幕三级久久日本| 久久久久性生活片| 寂寞人妻少妇视频99o| 国产精品女同一区二区软件| 美女中出高潮动态图| av免费观看日本| 日本欧美国产在线视频| 大香蕉久久网| 国产精品爽爽va在线观看网站| 亚洲人与动物交配视频| 女性生殖器流出的白浆| 国产精品三级大全| 成年免费大片在线观看| 精品熟女少妇av免费看| 久久97久久精品| 男人添女人高潮全过程视频| 九九久久精品国产亚洲av麻豆| 简卡轻食公司| 制服丝袜香蕉在线| 王馨瑶露胸无遮挡在线观看| 18+在线观看网站| 亚洲综合色惰| 久热久热在线精品观看| 亚洲激情五月婷婷啪啪| 亚洲国产色片| 国产熟女欧美一区二区| 久久ye,这里只有精品| tube8黄色片| 最黄视频免费看| 午夜激情福利司机影院| 国产免费一级a男人的天堂| 少妇猛男粗大的猛烈进出视频| 黄色欧美视频在线观看| 亚洲欧美一区二区三区国产| 亚洲电影在线观看av| av在线播放精品| 99久久精品一区二区三区| 99热网站在线观看| 国产91av在线免费观看| 五月伊人婷婷丁香| 乱系列少妇在线播放| 黄色一级大片看看| 久久99蜜桃精品久久| 寂寞人妻少妇视频99o| 我的女老师完整版在线观看| 中文字幕精品免费在线观看视频 | 国产日韩欧美亚洲二区| 欧美区成人在线视频| 一边亲一边摸免费视频| 99热6这里只有精品| 亚洲久久久国产精品| 久久久午夜欧美精品| 亚洲精品国产色婷婷电影| 一级黄片播放器| 亚洲最大成人中文| av国产免费在线观看| av国产久精品久网站免费入址| 乱码一卡2卡4卡精品| 欧美日本视频| 欧美bdsm另类| 在线观看免费高清a一片| 国产v大片淫在线免费观看| 欧美日本视频| 亚洲激情五月婷婷啪啪| 亚洲av电影在线观看一区二区三区| 干丝袜人妻中文字幕| 夫妻性生交免费视频一级片| 视频中文字幕在线观看| 欧美日本视频| 激情 狠狠 欧美| 亚洲第一av免费看| 日韩中文字幕视频在线看片 | 美女国产视频在线观看| 直男gayav资源| 亚洲天堂av无毛| 乱系列少妇在线播放| 看十八女毛片水多多多| 一级片'在线观看视频| 涩涩av久久男人的天堂| 国产精品av视频在线免费观看| 麻豆乱淫一区二区| 精品人妻偷拍中文字幕| 麻豆成人午夜福利视频| 边亲边吃奶的免费视频| 午夜老司机福利剧场| 色综合色国产| 嫩草影院新地址| 尾随美女入室| 搡老乐熟女国产| 秋霞伦理黄片| 18禁裸乳无遮挡动漫免费视频| 蜜桃亚洲精品一区二区三区| 亚洲一区二区三区欧美精品| 亚洲精品第二区| 日本黄大片高清| 一区二区三区乱码不卡18| 国产在线视频一区二区| 啦啦啦中文免费视频观看日本| 欧美高清成人免费视频www| 一级片'在线观看视频| 涩涩av久久男人的天堂| 丰满人妻一区二区三区视频av| 中文字幕制服av| 亚洲第一区二区三区不卡| av又黄又爽大尺度在线免费看| 亚洲欧美精品专区久久| 亚洲,一卡二卡三卡| 99热这里只有精品一区| 美女脱内裤让男人舔精品视频| 美女国产视频在线观看| 麻豆乱淫一区二区| 亚洲久久久国产精品| 免费不卡的大黄色大毛片视频在线观看| 国产黄色视频一区二区在线观看| 黄色一级大片看看| 黑人高潮一二区| 日本一二三区视频观看| 精品亚洲乱码少妇综合久久| av天堂中文字幕网| 亚洲美女视频黄频| 国产精品熟女久久久久浪| 欧美日本视频| 国产伦精品一区二区三区视频9| 美女cb高潮喷水在线观看| 国产av一区二区精品久久 | 国产精品国产三级国产av玫瑰| 乱系列少妇在线播放| 黄色一级大片看看| 一级a做视频免费观看| 国产无遮挡羞羞视频在线观看| 色网站视频免费| 舔av片在线| 涩涩av久久男人的天堂| 蜜臀久久99精品久久宅男| 免费观看无遮挡的男女| 日韩视频在线欧美| 26uuu在线亚洲综合色| 一本一本综合久久| 高清午夜精品一区二区三区| 一区二区三区乱码不卡18| 欧美区成人在线视频| 大片免费播放器 马上看| 国产黄片美女视频| av国产久精品久网站免费入址| 老司机影院毛片| 欧美成人a在线观看| 伦理电影大哥的女人| 亚洲欧美日韩卡通动漫| 久久久久久久久久久免费av| av卡一久久| 三级经典国产精品| 久久午夜福利片| 大陆偷拍与自拍| 91精品伊人久久大香线蕉| 国产在线免费精品| 亚洲高清免费不卡视频| 黄色日韩在线| 国产免费一区二区三区四区乱码| 精品国产露脸久久av麻豆| 人人妻人人添人人爽欧美一区卜 | 天美传媒精品一区二区| 欧美高清性xxxxhd video| 国产高潮美女av| 亚洲精品久久午夜乱码| 国产免费又黄又爽又色| 国产精品成人在线| 免费观看的影片在线观看| 亚洲欧美成人精品一区二区| 日日撸夜夜添| 久久精品夜色国产| 老司机影院毛片| 久久久精品94久久精品| 亚洲精品乱码久久久v下载方式| 欧美变态另类bdsm刘玥| 成人亚洲欧美一区二区av| 91精品伊人久久大香线蕉| 亚洲,欧美,日韩| 久久影院123| 精品一区二区三卡| 黑人高潮一二区| 九九爱精品视频在线观看| 国产真实伦视频高清在线观看| 少妇熟女欧美另类| 国内少妇人妻偷人精品xxx网站| 精品久久久久久久末码| 久久av网站| 五月玫瑰六月丁香| 涩涩av久久男人的天堂| 少妇人妻 视频| 中文字幕亚洲精品专区| 国产一区二区三区av在线| 国产黄色视频一区二区在线观看| 精品国产露脸久久av麻豆| 天天躁日日操中文字幕| 人妻一区二区av| 国产欧美另类精品又又久久亚洲欧美| 国产精品国产三级国产专区5o| 在线天堂最新版资源| 啦啦啦视频在线资源免费观看| 蜜桃在线观看..| 免费久久久久久久精品成人欧美视频 | 中文字幕av成人在线电影| 亚洲四区av| 精品国产一区二区三区久久久樱花 | 男女边吃奶边做爰视频| 一级爰片在线观看| 亚洲伊人久久精品综合| 熟女电影av网| 国产精品欧美亚洲77777| 亚洲美女搞黄在线观看| 一区在线观看完整版| 午夜激情福利司机影院| 五月开心婷婷网| 欧美zozozo另类| 亚洲久久久国产精品| 久久久精品94久久精品| 国产爽快片一区二区三区| 国产黄色免费在线视频| 中文字幕精品免费在线观看视频 | 久久久亚洲精品成人影院| 中文字幕精品免费在线观看视频 | 亚洲怡红院男人天堂| 超碰97精品在线观看| 久久久久久九九精品二区国产| 夜夜骑夜夜射夜夜干| 国产白丝娇喘喷水9色精品| 在线观看国产h片| 精品视频人人做人人爽| 日韩制服骚丝袜av| 亚洲精品456在线播放app| 国产亚洲5aaaaa淫片| 国产精品一区www在线观看| 黑丝袜美女国产一区| 一二三四中文在线观看免费高清| 欧美+日韩+精品| 国产高清国产精品国产三级 | 国产亚洲午夜精品一区二区久久| 女人久久www免费人成看片| 啦啦啦中文免费视频观看日本| 网址你懂的国产日韩在线| 18禁裸乳无遮挡免费网站照片| 伦理电影大哥的女人| 国语对白做爰xxxⅹ性视频网站| 欧美三级亚洲精品| 亚洲欧美日韩另类电影网站 | 亚洲精品久久午夜乱码| 色网站视频免费| 国产片特级美女逼逼视频| 免费人妻精品一区二区三区视频| 蜜桃久久精品国产亚洲av| 各种免费的搞黄视频| 久久久精品免费免费高清| 成人影院久久| 色综合色国产| 色婷婷久久久亚洲欧美| 在线观看免费日韩欧美大片 | 欧美精品亚洲一区二区| 亚洲美女黄色视频免费看| 精品亚洲乱码少妇综合久久| 亚洲内射少妇av| 麻豆成人av视频| 免费观看无遮挡的男女| 高清黄色对白视频在线免费看 | 亚洲性久久影院| 深爱激情五月婷婷| 亚洲成色77777| 老司机影院成人| 亚洲av中文av极速乱| 中文字幕亚洲精品专区| 免费播放大片免费观看视频在线观看| 三级国产精品片| 成人高潮视频无遮挡免费网站| 日韩av不卡免费在线播放| 最后的刺客免费高清国语| 夫妻性生交免费视频一级片| 国产黄片视频在线免费观看| 尾随美女入室| 久久久国产一区二区| 亚洲性久久影院| 国产成人a∨麻豆精品| 人体艺术视频欧美日本| 少妇 在线观看| 亚洲欧美成人综合另类久久久| 国产色爽女视频免费观看| 色哟哟·www| 久久亚洲国产成人精品v| 卡戴珊不雅视频在线播放| 亚洲精品aⅴ在线观看| 干丝袜人妻中文字幕| 国产成人91sexporn| 青春草视频在线免费观看| 午夜免费鲁丝| 夜夜看夜夜爽夜夜摸| 麻豆乱淫一区二区| 久久99热这里只有精品18| 秋霞伦理黄片| 日日摸夜夜添夜夜爱| 日本黄大片高清| 久久久久性生活片| 国产成人91sexporn| 免费看不卡的av| 亚洲aⅴ乱码一区二区在线播放| 我的老师免费观看完整版| 国产69精品久久久久777片| 精品久久久久久电影网| 小蜜桃在线观看免费完整版高清| 久久6这里有精品| 大话2 男鬼变身卡| 日本黄色片子视频| 日韩三级伦理在线观看| 欧美人与善性xxx| 久久国内精品自在自线图片| 亚洲精品久久久久久婷婷小说| 丝瓜视频免费看黄片| 国产精品欧美亚洲77777| www.av在线官网国产| 亚洲国产色片| 日韩伦理黄色片| 久久久精品免费免费高清| 国产成人免费无遮挡视频| 亚洲av日韩在线播放| 国产中年淑女户外野战色| 国产精品成人在线| 国产精品蜜桃在线观看| 亚洲欧洲国产日韩| 久久 成人 亚洲| 一级毛片aaaaaa免费看小| 精品久久久久久久久av| 国产伦理片在线播放av一区| 亚洲欧美日韩卡通动漫| 久久人妻熟女aⅴ| 亚洲人与动物交配视频| av播播在线观看一区| 中文字幕免费在线视频6| av免费在线看不卡| 人人妻人人看人人澡| 一级毛片久久久久久久久女| 伊人久久精品亚洲午夜| 王馨瑶露胸无遮挡在线观看| 国产av码专区亚洲av| 婷婷色综合www| 久久久亚洲精品成人影院| 性色av一级| 亚洲国产精品专区欧美| 国产 一区 欧美 日韩| 欧美精品国产亚洲| 国产精品99久久久久久久久| 国产成人免费无遮挡视频| 肉色欧美久久久久久久蜜桃| 简卡轻食公司| 蜜臀久久99精品久久宅男| 国产久久久一区二区三区| 亚洲精品乱久久久久久| 国产精品人妻久久久久久| 人人妻人人添人人爽欧美一区卜 | 精品人妻偷拍中文字幕| 欧美丝袜亚洲另类| 久久99精品国语久久久| 直男gayav资源| 国产国拍精品亚洲av在线观看| 午夜视频国产福利| 久久人人爽av亚洲精品天堂 | 亚洲,一卡二卡三卡| 夫妻性生交免费视频一级片| 日韩免费高清中文字幕av| 人妻夜夜爽99麻豆av| 婷婷色av中文字幕| www.色视频.com| 日韩一区二区视频免费看| 亚洲精品久久午夜乱码| 亚洲国产日韩一区二区| 午夜老司机福利剧场| 1000部很黄的大片| 亚洲精品日本国产第一区| 国产乱人视频| 色视频www国产| 日产精品乱码卡一卡2卡三| 欧美3d第一页| 久久久久久久久大av| 联通29元200g的流量卡| 18禁在线播放成人免费| 亚洲性久久影院| 一级毛片我不卡| 亚洲精品aⅴ在线观看| 国产乱来视频区| 日韩不卡一区二区三区视频在线| av在线老鸭窝| 国产亚洲最大av| 国产欧美亚洲国产| 大香蕉97超碰在线|