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

    奎屯河流域地下水水化學(xué)特征及成因分析

    2024-06-12 11:56:55江宇威李巧陶洪飛馬合木江·艾合買提
    人民長江 2024年5期
    關(guān)鍵詞:水化學(xué)承壓水平原

    江宇威 李巧 陶洪飛 馬合木江·艾合買提

    摘要:了解地下水水化學(xué)特征及其成因,有助于評(píng)估地下水水質(zhì)、識(shí)別污染源、制定水資源管理策略,保障人類健康和環(huán)境?;?017年和2023年在新疆奎屯河流域相同點(diǎn)位采集的共78組地下水水樣,采用Piper三線圖、Gibbs圖、飽和指數(shù)(SI)和離子比值法分析該區(qū)域地下水特征及成因。結(jié)果表明:① 奎屯河流域2023年地下水水化學(xué)特征類型相較于2017年趨于集中,陽離子和陰離子分別向Na型和Cl型轉(zhuǎn)化。② 山前洪積傾斜礫質(zhì)平原和沖洪積細(xì)土平原主要受巖石風(fēng)化作用控制,沖湖積細(xì)土平原主要受蒸發(fā)濃縮作用控制且有增強(qiáng)的趨勢。中下游細(xì)土平原的離子組分主要來源于蒸發(fā)巖溶解且有增強(qiáng)的趨勢,山前洪積傾斜礫質(zhì)平原的離子組分主要源于碳酸鹽和硅酸鹽礦物溶解。下游沖湖積細(xì)土平原受陽離子交替吸附作用較強(qiáng),導(dǎo)致Na+濃度上升??秃恿饔虻叵滤x子在沖洪積細(xì)土平原北部和沖湖積細(xì)土平原出現(xiàn)富集,部分指標(biāo)已超過中國飲用水標(biāo)準(zhǔn),且已有10%的水樣不適合灌溉。研究成果可為奎屯河流域地下水合理利用以及防治提供參考。

    關(guān) 鍵 詞:地下水;水化學(xué)特征;空間分布;成因分析;飽和指數(shù);奎屯河流域

    中圖法分類號(hào):X523;P641.12

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

    DOI:10.16232/j.cnki.1001-4179.2024.05.009

    0 引 言

    地下水化學(xué)成分形成是多種水文地球化學(xué)過程相互作用的結(jié)果,包括巖石的溶解沉淀、地下水和巖石土壤之間的離子交換、人類活動(dòng)和人為輸入污染物等[1。西北干旱半干旱地區(qū)水資源匱乏,地下水是當(dāng)?shù)孛癖娚詈徒?jīng)濟(jì)發(fā)展的重要支撐,其水化學(xué)特征對(duì)當(dāng)?shù)氐霓r(nóng)業(yè)灌溉以及日常用水的影響重大。研究地下水的化學(xué)特征及組成,有助于了解地下水的起源、補(bǔ)給機(jī)制和水質(zhì)演變規(guī)律,可以更有效地進(jìn)行水資源規(guī)劃、管理和可持續(xù)利用[2

    對(duì)新疆奎屯河流域地下水水化學(xué)組成及其成因的探討具有特殊意義。首先,奎屯河流域內(nèi)的烏蘇市經(jīng)濟(jì)較發(fā)達(dá)、人口密集且產(chǎn)業(yè)集中,農(nóng)業(yè)作為三大經(jīng)濟(jì)支撐之一,其用水主要來源于地下水。因此,對(duì)該地區(qū)地下水水化學(xué)組成及其成因的研究,有助于判別地下水灌溉適宜性以及種植作物類別,具有明確的實(shí)用價(jià)值。其次,奎屯河流域?qū)儆诘湫偷母呱榈叵滤植紖^(qū)域,砷元素通過直接飲用或作物吸收間接進(jìn)入人體,危害人體健康。而地下水中鐵錳氧化物還原、SO2-4的還原和硫化礦物沉淀等均會(huì)對(duì)砷含量產(chǎn)生影響[3。因此,研究地下水水化學(xué)組成及其成因,可以為研究區(qū)后續(xù)探究砷氟環(huán)境形成機(jī)制提供一定參考。

    前人對(duì)奎屯河流域地下水水化學(xué)特征及成因研究較多,江軍等[4在研究砷遷移規(guī)律時(shí),分析了2017年水化學(xué)特征及礦物飽和指數(shù);劉亞楠等[5在研究地下水及沉積物化學(xué)組分特征對(duì)氟釋放的影響時(shí),分析了2017年地下水化學(xué)成因機(jī)制。張福初等[6研究較為系統(tǒng),分析了2019年以前全流域不同地貌下的水化學(xué)特征及成因。以上學(xué)者對(duì)該區(qū)域的研究未涉及時(shí)間變化特征分析,且研究時(shí)間較早。近年來氣候變暖,水資源總量增加和氣溫升高等因素會(huì)影響淺層水化學(xué)特征及成因。因此,針對(duì)相同點(diǎn)位、不同時(shí)期水化學(xué)組成變化的研究可以精確分析地下水水化學(xué)特征變化的程度及可能原因,對(duì)水資源合理利用具有參考意義。

    本文在2017年采集39組地下水樣的基礎(chǔ)上,于2023年在相同點(diǎn)位重新采集。采用Piper三線圖、Gibbs圖、礦物指數(shù)、離子比值法和氯堿指數(shù)分析水化學(xué)特征、空間分布和水化學(xué)成因,并對(duì)比分析其變化,以期為后續(xù)的砷氟成因規(guī)律探討提供基礎(chǔ)參考,為研究區(qū)地下水的開發(fā)利用提供科學(xué)依據(jù)。

    1 研究區(qū)概況

    奎屯河流域坐落于新疆天山北麓,準(zhǔn)噶爾盆地西南部。地處83°22′E~85°47′E,43°30′N~47°04′N。東西寬度約160 km,南北長度約240 km,總面積達(dá)2.83萬km2。氣候類型為溫帶大陸性干旱氣候,具有少雨、高蒸發(fā)和大溫差等特點(diǎn),年均降水量為162.6 mm,年均蒸發(fā)量為1 208.3 mm。

    研究區(qū)高程范圍在247~1 210 m之間,地勢有南高北低、東高西低的特點(diǎn)。隨著地勢由南向北降低,地貌逐漸從洪積傾斜礫質(zhì)平原轉(zhuǎn)變?yōu)闆_洪積細(xì)土平原(圖1),最后演變?yōu)闆_積湖積細(xì)土平原。烏伊公路以南區(qū)域?yàn)樯角昂榉e傾斜礫質(zhì)平原,巖性為卵礫石或砂礫石,地下水類型為單層潛水且潛水埋深大于20 m。作為研究區(qū)地下水的主要補(bǔ)給區(qū),該區(qū)域主要接收山區(qū)地下水的側(cè)向補(bǔ)給,以及河道水、洪水、渠系水滲漏。烏伊公路以北至一二三團(tuán)區(qū)域?yàn)闆_洪積細(xì)土平原,巖性為砂礫石、砂和黏性土,地下水類型為多層潛水-承壓水[7,潛水埋深一般小于5 m。該區(qū)域主要接收山前洪積傾斜礫質(zhì)平原的側(cè)向補(bǔ)給,以及灌區(qū)內(nèi)的滲漏補(bǔ)給。一二三團(tuán)以北為沖湖積細(xì)土平原,巖性為粉細(xì)砂以及黏土,地下水類型為多層潛水-承壓水,潛水埋深一般小于10 m,主要接收沖洪積細(xì)土平原的側(cè)向補(bǔ)給。由于潛水位埋藏淺且蒸發(fā)蒸騰作用強(qiáng),導(dǎo)致該區(qū)域?yàn)檠芯繀^(qū)地下水主要排泄區(qū)8。

    2 材料與方法

    2.1 水樣采集與分析

    本次研究使用的水樣為2017年6月于奎屯河流域已有機(jī)井采集的39組水樣和2023年7月于原位點(diǎn)采集的39組水樣,采樣點(diǎn)分布見圖1。水樣的采集、保存及送樣過程嚴(yán)格按照HJ/T 164-2004《地下水環(huán)境監(jiān)測技術(shù)規(guī)范》[9的要求執(zhí)行。樣品的測試由新疆地礦局第二水文地質(zhì)大隊(duì)化驗(yàn)室完成,測試項(xiàng)目包括可溶性鹽濃度(EC)、氧化還原電位(Eh)、pH、Na+、K+、Ca2+、Mg2+、總?cè)芙夤腆w(TDS)、總硬度(TH)、Cl-、SO2-4、HCO-3、NO-3、F-、As。運(yùn)用陰陽離子平衡法計(jì)算可知,所有水樣數(shù)據(jù)的電荷平衡絕對(duì)誤差均小于5%,為可靠數(shù)據(jù)。

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

    研究區(qū)沖洪積細(xì)土平原和沖湖積細(xì)土平原的地下水類型按埋深劃分為烏伊公路以南的地下水均為潛水[7、烏伊公路以北井深小于100 m為淺層承壓水、大于100 m為深層承壓水。研究過程中使用的散點(diǎn)圖基于Origin 2018繪制,空間分布圖基于ArcMAP 10.2繪制,Gibbs圖基于Excel 2010繪制。

    3 結(jié)果與討論

    3.1 地下水水化學(xué)特征

    3.1.1 地下水水化學(xué)組分統(tǒng)計(jì)分析

    2017年水化學(xué)組分統(tǒng)計(jì)結(jié)果表明(表1),pH范圍為7.04~9.00,承壓水整體為弱堿性至堿性。溶解性總固體(TDS)范圍為114.48~11 632.96 mg/L,其中淡水(TDS<1 000 mg/L)和微咸水(1 000 mg/L≤TDS<3 000 mg/L)分別占比74%和15%??傆捕龋═H)范圍為18.10~3 639.51 mg/L,其中軟水(75 mg/L450 mg/L)分別占比33.30%和30.76%。水的硬度主要由鈣離子(Ca2+)和鎂離子(Mg2+)含量決定,離子組分Na+、Ca2+、Mg2+、Cl-、SO2-4的變異系數(shù)較大,存在較強(qiáng)的空間異質(zhì)性。主要陽離子平均濃度排序?yàn)镹a+>Ca2+>Mg2+>K+,主要陰離子平均濃度排序?yàn)镾O2-4>Cl->HCO-3。

    2023年水化學(xué)組分統(tǒng)計(jì)結(jié)果表明(表2),pH范圍為7.18~9.21,極值上升但均值明顯下降。淡水和微咸水分別占比69.2%和17.4%,微咸水占比呈現(xiàn)增長趨勢。軟水和極硬水分別占比28.20%和30.76%,極硬水占比未發(fā)生變化,部分軟水向微硬水轉(zhuǎn)化。陰陽離子平均濃度排序未出現(xiàn)變化。相較2017年,淺層承壓水主要水化學(xué)組分明顯增大,按平均濃度的增長幅度排序?yàn)镃l-(137.47%)>SO2-4(67.72%)>Na+(37.09%)>TDS(26.11%)>HCO-3(17.96%)>Ca2+(13.92%)>TH(10.14%)>Mg2+(-6.31%)>K+(-46.87%)。深層承壓水主要水化學(xué)組分明顯降低,按平均濃度的增長幅度排序?yàn)镃a2+(3.55%)>HCO-3(2.64%)>K+(-8.37%)>Cl-(-12.63%)>Na+(-19.02%)>TDS(-21.91%)>SO2-4(-25.39%)>TH(-27.37%)>Mg2+(-55.52%)。

    3.1.2 水化學(xué)特征

    Piper三線圖[10是根據(jù)地下水樣品的離子組成和濃度比例進(jìn)行分類,通過圖中不同區(qū)域的分布確定地下水水化學(xué)類型(圖2)。 圖2中的菱形區(qū)域被劃分為5,6,7,8,9區(qū),依次對(duì)應(yīng)Ca-HCO3型、Ca-SO4型、Na-Cl或Na-SO4型、Na-HCO3型和混合型。三角形區(qū)域被劃分為A、B、C、D、E、F、G區(qū),依次對(duì)應(yīng)Ca型、無主導(dǎo)型、Mg型、Na型、HCO3型、SO4型和Cl型。

    從水化學(xué)類型看,2017年淺層承壓水集中分布在7區(qū)和9區(qū),表明水化學(xué)類型主要為Na-Cl型、混合型或Na-SO4型。潛水分布在5區(qū),水化學(xué)類型為Ca-HCO3型。深層承壓水在5,7,9區(qū)均有分布,水化學(xué)類型為Ca-HCO3型、Na-Cl型或Na-SO4型、混合型。2023年水化學(xué)類型分布和2017年類似,但出現(xiàn)5區(qū)向7區(qū)和9區(qū)轉(zhuǎn)化的趨勢,即部分Ca-HCO3型向Na-Cl型或Na-SO4型、混合型轉(zhuǎn)化。

    從陽離子看,2017年淺層承壓水陽離子集中分布在B和D區(qū),表明優(yōu)勢陽離子為無主導(dǎo)型和Na型。潛水分布在A區(qū),表明優(yōu)勢陽離子為Ca型。深層承壓水陽離子集中分布在A,B,D區(qū),表明優(yōu)勢陽離子為Ca型、無主導(dǎo)型和Na型。2023年淺層承壓水陽離子集中分布在D區(qū),表明5 a間Ca型向Na型轉(zhuǎn)化。潛水和深層承壓水的陽離子分布變化不大。

    從陰離子看,2017年淺層承壓水和深層承壓水陰離子均集中分布在B,E,F(xiàn)區(qū),表明優(yōu)勢陰離子為無主導(dǎo)型、HCO3型和SO4型。潛水分布在E區(qū),優(yōu)勢陰離子為HCO3型。2023年淺層承壓水出現(xiàn)B區(qū)向F區(qū)和G區(qū)遷移的趨勢,表明無主導(dǎo)型在向SO4型和Cl型轉(zhuǎn)化。潛水和深層承壓水陰離子均出現(xiàn)E區(qū)向B區(qū)和G區(qū)遷移,表明部分HCO3型向無主導(dǎo)型和Cl型轉(zhuǎn)化。

    3.2 地下水中主要離子空間分布

    基于ArcGIS軟件中的克里金插值法繪制研究區(qū)地下水優(yōu)勢離子空間分布圖(圖3)。除Ca2+外,其余指標(biāo)根據(jù)GB/T 14848-2017《地下水質(zhì)量標(biāo)準(zhǔn)》規(guī)定分為5類,其中,Ⅰ~Ⅲ類水適合飲用,Ⅵ~Ⅴ類水不適合飲用。根據(jù)2023年的數(shù)據(jù)可以得出,TDS、SO2-4和Cl-的高值區(qū)集中于一二六團(tuán)南部、一二三團(tuán)、車排子鎮(zhèn)南部、一二五團(tuán)北部、一二九團(tuán)北部和甘河子鎮(zhèn)。在一二七團(tuán)L9樣本點(diǎn)和一二六團(tuán)L2樣本點(diǎn)出現(xiàn)多項(xiàng)指標(biāo)異常值,TDS分別為7 262.00 mg/L和6 204.00 mg/L,Cl-濃度分別為2 269.71 mg/L和1 808.67 mg/L,Na+濃度分別為1 734.41 mg/L和1 421.51 mg/L,SO2-4濃度分別為2 605.96 mg/L和2 029.03 mg/L。對(duì)比2017年和2023年兩期數(shù)據(jù),包括西湖鎮(zhèn)、八十四戶鄉(xiāng)、西大溝鎮(zhèn)和九間樓鄉(xiāng)在內(nèi)的上游地區(qū),主要離子濃度維持在較低水平,未出現(xiàn)明顯變化。中下游地區(qū)的一二三團(tuán)、古爾圖鎮(zhèn)和一二六團(tuán)變化較明顯且規(guī)律相似。TDS、SO2-4、Cl-三項(xiàng)指標(biāo)呈明顯增長趨勢,大面積的Ⅲ類水向Ⅴ類水轉(zhuǎn)化,Ca2+呈現(xiàn)增長趨勢但未超標(biāo)。Na+濃度增長幅度相對(duì)較緩,下游右岸大面積的Ⅰ類水向Ⅳ類水轉(zhuǎn)化,古爾圖鎮(zhèn)的Ⅴ類水呈減少趨勢,表明研究區(qū)中下游地區(qū)的地下水已不適合飲用。

    3.3 地下水水化學(xué)特征成因分析

    3.3.1 水化學(xué)作用分析

    Gibbs圖是根據(jù)c(Na+)/c(Na++Ca2+)和c(Cl-)/c(Cl-+HCO-3)的離子濃度比與TDS之間的比值關(guān)系繪制散點(diǎn)圖,同時(shí)劃分了蒸發(fā)濃縮、巖石風(fēng)化和降雨控制3個(gè)區(qū)域,從點(diǎn)位的分布位置判斷樣本點(diǎn)水化學(xué)特征形成的主要機(jī)制[11

    對(duì)2017年和2023年的地下水樣本進(jìn)行分析(圖4),研究區(qū)地下水樣大部分分布在巖石風(fēng)化型區(qū)域,表明水巖相互作用在地下水水質(zhì)形成的水文地球化學(xué)過程中發(fā)揮著主導(dǎo)作用。部分點(diǎn)位落在3個(gè)區(qū)域外,可能受到離子交換作用或人類活動(dòng)的影響[12。山前洪積傾斜礫質(zhì)平原的地下水徑流條件較好,溶濾作用較顯著。沖湖積細(xì)土平原因地形平緩和水交替循環(huán)條件較差,地下水的溶濾作用相對(duì)較弱,導(dǎo)致蒸發(fā)濃縮型點(diǎn)位集中分布在沖湖積細(xì)土平原,巖石風(fēng)化型集中分布在山前洪積傾斜礫質(zhì)平原。對(duì)比兩期數(shù)據(jù)發(fā)現(xiàn),淺層承壓水中c(Na+)/c(Na++Ca2+)平均增大0.018,變化不顯著;但TDS增大26.11%,導(dǎo)致點(diǎn)位向蒸發(fā)濃縮區(qū)域集中。c(Cl-)/c(Cl-+HCO-3)平均增大0.372,是導(dǎo)致點(diǎn)位向區(qū)域外遷移的主要原因。對(duì)比兩期Cl-和HCO-3的平均濃度發(fā)現(xiàn),Cl-平均濃度增長137.47%,導(dǎo)致c(Cl-)/c(Cl-+HCO-3)增大。因此出現(xiàn)集中在巖石風(fēng)化型區(qū)域的樣本點(diǎn)開始向蒸發(fā)濃縮區(qū)域以及區(qū)域外遷移。深層承壓水中c(Na+)/c(Na++Ca2+)平均增大0.058,c(Cl-)/c(Cl-+HCO-3)平均增大0.062,變化不顯著。

    3.3.2 研究區(qū)礦物飽和指數(shù)

    為了評(píng)估地下水化學(xué)形成過程中水巖相互作用的影響,可以使用PHREEQC軟件計(jì)算主要礦物的飽和指數(shù)(SI)來判定它們的溶解程度[13。SI的大小表示礦物的溶解程度,正值表示過飽和,零值表示飽和,負(fù)值表示欠飽和[14。

    計(jì)算研究區(qū)2017年和2023共78組水樣的礦物飽和指數(shù)(SI)(圖5),隨著TDS增大,文石和方解石的飽和指數(shù)增長不顯著,TDS處于0~1 000 mg/L時(shí),白云石和石膏的飽和指數(shù)增長明顯,表明白云石和石膏的溶解促進(jìn)了TDS濃度增大。TDS濃度大于1 000 mg/L后,礦物飽和指數(shù)變化不明顯、接近飽和,且點(diǎn)位集中于沖湖積細(xì)土平原,表明TDS濃度增大受到蒸發(fā)濃縮作用的影響。一二六團(tuán)作為離子濃度高值區(qū),其中L2采樣點(diǎn)的文石、方解石和白云石在兩個(gè)時(shí)間點(diǎn)均處于過飽和狀態(tài),但2023年較2017年,飽和指數(shù)呈降低趨勢。L7采樣點(diǎn)2017年方解石和白云石處于過飽和狀態(tài),2023年在此基礎(chǔ)上出現(xiàn)文石過飽和現(xiàn)象,且飽和指數(shù)呈上升趨勢。此外,甘河子鎮(zhèn)的K23采樣點(diǎn)在2017年僅有方解石過飽和,2023年文石、方解石和白云石均處于過飽和狀態(tài),呈上升趨勢。

    3.3.3 主要離子來源

    離子比值法(圖6)是通過比較不同離子毫克相當(dāng)濃度(mEq,符號(hào)為E)間的關(guān)系,判斷離子主要來源。E(Na+)/E(Cl-[15可以判斷Na+是否來源于鹽巖和硅酸鹽巖溶解。研究區(qū)大部分地下水樣品落在y=x線附近且偏下方,說明Na+濃度略高于Cl-濃度,巖鹽溶解是水樣中Na+的主要來源,少部分Na+來源于硅酸鹽溶解或陽離子交替吸附作用。

    E(HCO-3+SO2-4)/E(Mg2++Ca2+)可以判斷該區(qū)水化學(xué)是否由碳酸鹽巖或者石膏溶解控制[16。研究區(qū)上游區(qū)域地下水樣品落在y=x線附近,說明Mg2+和Ca2+來源主要是碳酸鹽和硅酸鹽溶解。下游區(qū)域的地下水樣品落在y=x線的下方,說明Mg2+和Ca2+的來源是硅酸鹽和蒸發(fā)巖礦物的溶解,下游區(qū)域以蒸發(fā)濃縮作用為主導(dǎo),說明蒸發(fā)巖礦物的溶解為主要來源。對(duì)比2017年和2023年兩期數(shù)據(jù),個(gè)別落于y=x線上方的點(diǎn)位向y=x線及下方遷移,表明蒸發(fā)巖礦物的溶解提供的Mg2+和Ca2+增多。采用E(Ca2++Mg2+)/E(HCO-3)可以進(jìn)一步判斷Mg2+和Ca2+的來源是否為碳酸鹽溶解,研究區(qū)點(diǎn)位分布在y=x線附近及下方,表明碳酸鹽溶解是Mg2+和Ca2+的主要來源。2023年較2017年,y=x線下方點(diǎn)位分布增多,表明碳酸鹽溶解呈現(xiàn)增長趨勢。

    E(Ca2+)/E(Mg2+) 可以分析Mg2+和Ca2+來源與白云石和方解石溶解的關(guān)系17。研究區(qū)2017年樣本點(diǎn)分布較均勻,在各個(gè)區(qū)域都有點(diǎn)位分布,表明白云石、方解石和石膏對(duì)Mg2+和Ca2+的貢獻(xiàn)相當(dāng)。2023年樣本點(diǎn)集中分布于E(Ca2+)/E(Mg2+)>2的區(qū)域,表明石膏溶解是Ca2+的主要來源,且石膏溶解的貢獻(xiàn)率呈增長趨勢。為了進(jìn)一步驗(yàn)證Ca2+的來源,采用E(Ca2+)/E(HCO-3)判斷Ca2+是否由石膏溶解產(chǎn)生,研究區(qū)點(diǎn)位分布在y=x線附近及下方,對(duì)比E(Ca2++Mg2+)/E(HCO-3)的關(guān)系可以得出碳酸鹽溶解主要由石膏產(chǎn)生。2023年較2017年,y=x線下方點(diǎn)位分布增多,表明石膏溶解呈現(xiàn)增長趨勢。

    E(SO2-4+Cl-)/E(HCO-3)的關(guān)系可以判斷SO2-4和Cl-的主要來源[18。研究區(qū)地下水樣品落在y=x線上方,表明蒸發(fā)巖溶解是SO2-4和Cl-的主要來源。

    E(Ca2++Mg2+-HCO-3)/E(SO2-4-Na++Cl-)可以進(jìn)一步分析地下水中SO2-4的來源。研究區(qū)大部分樣本點(diǎn)分布在y=x線附近,小部分分布在y=x線下方,表明除石膏溶解外還有其他硫酸鹽的溶解提供SO2-4

    E(Ca2++Na+-Cl-)/E(SO2-4)可以分析SO2-4和石膏-芒硝溶解的關(guān)系[19。研究區(qū)大部分樣本點(diǎn)分布在y=x線附近,表明石膏和芒硝溶解提供了SO2-4。個(gè)別點(diǎn)位分布在y=x線上方,表明SO2-4濃度可能受到了人類活動(dòng)的影響。

    3.3.4 陽離子交替吸附作用

    氯堿指數(shù)(CAI)是用于評(píng)估陽離子交替吸附強(qiáng)度的指標(biāo)[20。通過觀察CAI1與CAI2的絕對(duì)值變化判斷陽離子交換的強(qiáng)度。絕對(duì)值越大,陽離子交換作用越顯著,即陽離子交替吸附的可能性更高[21。

    根據(jù)CAI1與CAI2關(guān)系圖(圖7)得出,研究區(qū)樣本點(diǎn)集中分布在CAI1和CAI2均小于0的區(qū)域,絕對(duì)值較小,表明研究區(qū)主要發(fā)生Ca2+、Mg2+置換Na+、K+的反應(yīng),但不顯著。在CAI1和CAI2均大于0的區(qū)域,大部分點(diǎn)位位于沖洪積細(xì)土平原,包括位于山前洪積傾斜礫質(zhì)平原的潛水點(diǎn)位,表明該區(qū)域主要發(fā)生Na+、K+置換Ca2+、Mg2+的反應(yīng),與上游區(qū)域Na+濃度長期維持低值的規(guī)律一致。在CAI1和CAI2均小于0區(qū)域,存在個(gè)別絕對(duì)值較高的點(diǎn)位,且集中分布在沖湖積細(xì)土平原。

    3.4 灌溉水質(zhì)評(píng)價(jià)

    USSL圖[22以EC為橫坐標(biāo),鈉吸附比(SAR)為縱坐標(biāo)繪制,EC指示鹽害程度,SAR指示堿害程度,位于C1S1、C2S1和C3S1區(qū)域的地下水適宜灌溉。Wilcox圖以EC為橫坐標(biāo),可溶性鈉百分比(SSP)為縱坐標(biāo)繪制,SSP指示水的軟硬程度,位于水質(zhì)極好、水質(zhì)良好和水質(zhì)適宜區(qū)域的地下水適宜灌溉。對(duì)USSL圖和Wilcox圖分區(qū),可以判斷地下水對(duì)灌溉的適宜程度。USSL圖結(jié)果表明(圖8),研究區(qū)大部分地下水集中在S1區(qū)域,受堿害程度低。EC在C1~C3均勻分布,表明高鹽危害仍在可控范圍。有14%的水樣不適合灌溉,且在不適合灌溉的水樣中,2023淺層承壓水占比50%。Wilcox圖結(jié)果表明(圖9),研究區(qū)有18%的水樣不適合灌溉,且在不適合灌溉的水樣中,2023年淺層承壓水占比38.9%。

    4 討 論

    沖洪積細(xì)土平原北部和沖湖積細(xì)土平原是離子濃度高值區(qū),呈現(xiàn)濃度上升趨勢。原因是地下水流向?yàn)橛赡舷虮?,中下游?xì)土平原屬于地下水排泄區(qū)域,地下水流速緩慢導(dǎo)致離子在北部的沖湖積細(xì)土平原出現(xiàn)累積。

    淺層承壓水由蒸發(fā)濃縮作用和巖石風(fēng)化作用共同控制。Gibbs圖結(jié)果表明,2017年采樣點(diǎn)位分布在蒸發(fā)濃縮區(qū)域和區(qū)域外,結(jié)合氯堿指數(shù)發(fā)現(xiàn),區(qū)域外點(diǎn)位受到陽離子交替吸附作用,Ca2+置換出Na+導(dǎo)致c(Na+)/c(Na++Ca2+)增大。2023年由于TDS增大,部分點(diǎn)位向蒸發(fā)濃縮區(qū)域遷移。但在c(Cl-)/c(Cl-+HCO-3)與TDS關(guān)系圖中,Cl-濃度增幅較大導(dǎo)致c(Cl-)/c(Cl-+HCO-3)增大,個(gè)別點(diǎn)位移動(dòng)至區(qū)域外。深層承壓水主要由巖石風(fēng)化作用控制,且2023年較2017年無明顯變化。點(diǎn)位均分布在巖石風(fēng)化型和區(qū)域外,區(qū)域外的點(diǎn)位受陽離子交替吸附作用影響。2023年較2017年,各類礦物的平均飽和指數(shù)均呈增長趨勢,導(dǎo)致TDS升高,淺層承壓水中巖鹽(NaCl)的平均飽和指數(shù)增幅(0.51)最大,一定程度上解釋了淺層承壓水中Cl-濃度大幅增長的現(xiàn)象。

    綜合上述分析,研究區(qū)中下游離子濃度增大主要是蒸發(fā)濃縮和巖石風(fēng)化共同作用的結(jié)果。但Cl-和SO2-4濃度增長幅度遠(yuǎn)大于其他陽離子,表明離子濃度增長可能受到人類活動(dòng)的影響。

    5 結(jié) 論

    本文基于2017年和2023年灌溉期共78組原位點(diǎn)水樣,通過Piper三線圖分析地下水化學(xué)特征,采用Gibbs圖、飽和指數(shù)、氯堿指數(shù)和離子比值法分析其成因。結(jié)果表明:

    (1)2017年陽離子以Na和Ca型為主,陰離子以SO2-4和無主導(dǎo)型為主。2023年陽離子以Na型為主,陰離子以Cl型為主。2023年較2017年,水化學(xué)類型趨于集中。

    (2)研究區(qū)山前洪積傾斜礫質(zhì)平原和沖洪積細(xì)土平原南部受巖石風(fēng)化作用控制,沖洪積細(xì)土平原北部和沖湖積細(xì)土平原整體受蒸發(fā)和巖石風(fēng)化濃縮共同控制。2023年較2017年,中下游細(xì)土平原淺層承壓水受到的蒸發(fā)濃縮作用增強(qiáng),深層承壓水沒有明顯變化。

    (3)研究區(qū)山前洪積傾斜礫質(zhì)平原的離子組分主要源于碳酸鹽和硅酸鹽礦物溶解。中下游細(xì)土平原的離子組分主要來源于蒸發(fā)巖溶解且有增強(qiáng)的趨勢。下游沖湖積細(xì)土平原受陽離子交替吸附作用較強(qiáng)。

    (4)研究區(qū)以兵團(tuán)一二六團(tuán)為主的下游地區(qū),地下水多項(xiàng)指標(biāo)超過Ⅲ類水級(jí)別,不適合飲用。從灌溉適宜性看,已有超過10%的水樣水適合灌溉。在不適合灌溉的水樣中,2023年的淺層承壓水占比最大。

    參考文獻(xiàn):

    [1]王玉蛟,賈亮亮,楊志光,等.淺層地下水水化學(xué)特征及適宜性評(píng)價(jià)[J].排灌機(jī)械工程學(xué)報(bào),2023,41(4):363-368.

    [2]胡爽,武剛,吳彬,等.阜康市平原區(qū)近20年地下水化學(xué)演變及成因[J].環(huán)境科學(xué)學(xué)報(bào),2023,43(5):79-88.

    [3]袁翰卿,李巧,陶洪飛,等.新疆奎屯河流域地下水砷富集因素[J].環(huán)境化學(xué),2020,39(2):524-530.

    [4]江軍,鮮虎勝,李巧,等.奎屯河流域地下水地球化學(xué)特征及其對(duì)砷運(yùn)移的影響[J].環(huán)境化學(xué),2021,40(6):1775-1786.

    [5]劉亞楠,李巧,宿彥鵬,等.奎屯河流域地下水及沉積物化學(xué)組分特征對(duì)氟釋放的影響[J].長江科學(xué)院院報(bào),2023,40(10):59-65.

    [6]張福初,吳彬,高凡,等.奎屯河流域平原區(qū)地下水水化學(xué)特征及成因分析[J].環(huán)境科學(xué)研究,2021,34(7):1663-1671.

    [7]葛國際.烏蘇市平原區(qū)水文地質(zhì)條件簡析[J].地下水,2016,38(2):176-178.

    [8]沈蕊芯,呂樹萍,杜明亮,等.新疆奎屯河流域平原區(qū)地下水資源量演變情勢[J].干旱區(qū)研究,2020,37(4):839-846.

    [9]國家環(huán)境保護(hù)總局.地下水環(huán)境監(jiān)測技術(shù)規(guī)范:HJ/T 164-2004[S].北京:中國環(huán)境科學(xué)出版社,2004.

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

    [11]趙艷麗,王敏,李常鎖.山東省淄河流域地下水化學(xué)污染源解析[J].人民長江,2022,53(12):37-43,49.

    [12]王繼華,于福榮,田鵬州,等.南太行山前沖積平原潛水水化學(xué)特征及控制因素分析[J].干旱區(qū)資源與環(huán)境,2023,37(5):130-138.

    [13]韓朝輝,王郅睿,田輝,等.漢中盆地地下水水化學(xué)特征及其成因研究[J].西北地質(zhì),2023,56(4):263-273.

    [14]魯涵,曾妍妍,周金龍,等.巴楚縣淺層地下水中氟的分布特征及影響因素分析[J].環(huán)境化學(xué),2021,40(11):3455-3463.

    [15]TALIB A M,TANG Z,SHAHAB A,et al.Hydrogeochemical characterization and suitability assessment of groundwater:a case study in Central Sindh,Pakistan[J].International Journal of Environmental Research and Public Health,2019,16(5):886-896.

    [16]蘇東,龔緒龍,楊磊,等.常州市地下水化學(xué)特征與成因分析[J].地質(zhì)論評(píng),2023,69(3):1039-1049.

    [17]劉元晴,周樂,呂琳,等.河北省順平縣地下水化學(xué)特征及其成因分析[J].環(huán)境科學(xué),2023,44(5):2601-2612.

    [18]付舒,毛龍富,劉宏,等.阿子營黃龍泉域地下水化學(xué)特征及演化[J].地球與環(huán)境,2023,51(2):153-163.

    [19]葛婷婷.克里雅河流域平原區(qū)潛水水質(zhì)評(píng)價(jià)及水質(zhì)成因分析[D].烏魯木齊:新疆農(nóng)業(yè)大學(xué),2022.

    [20]鄒嘉文.黑龍港典型壓采區(qū)地下水水化學(xué)變化特征及影響因素研究[D].邯鄲:河北工程大學(xué),2022.

    [21]王婧,高世昌,方媛,等.清水河平原地下水水化學(xué)特征及成因分析[J].人民黃河,2022,44(增2):153-156.

    [22]JIE M,SONG C,DIANDIAN D,et al.Hydrochemical characteristics and groundwater quality appraisal for irrigation uses in the Lan-gan region,Northern Anhui Province,East China[J].Water Science & Technology:Water Supply,2022,22(1):686-696.

    (編輯:劉 媛)

    Chemical characteristics and genetic analysis of groundwater in Kuitun River BasinJIANG Yuwei1,2,LI Qiao1,2,TAO Hongfei1,2,MAHEMUJIANG Aihemaiti1,2

    (1.College of Water Conservancy and Civil Engineering,Xinjiang Agricultural University,Urumqi 830052 China;2.Xinjiang Key Laboratory of Water Conservancy Engineering Safety and Water Disaster Prevention,Urumq 830052,China)

    Abstract:The analysis of groundwater chemical characteristics and their origins is essential for evaluating groundwater quality,identifying pollution sources,developing water resource management strategies,and safeguarding human health and the environment.Based on a total of 78 groundwater samples collected at the same sites in the Kuitun River Basin in 2017 and 2023,the characteristics and causes of groundwater in this area were analyzed by Piper trilinear diagram,Gibbs diagram,saturation index (SI)method and ion ratio method.The findings indicated that:① There was a converging trend in the types of groundwater chemical characteristics in the Kuitun River Basin in 2023 compared to 2017,with cations transforming into Na type and anions into Cl type.② Rock weathering primarily influenced the piedmont diluvial inclined gravel plain and alluvial fine soil plain,while evaporation and concentration mainly affect the alluvial lacustrine fine soil plain with a tendency to increase.The ion components in the middle and lower reaches of the fine soil plain were predominantly derived from evaporative karst decomposition with an increasing trend,whereas those in the fore mountain diluvial inclined gravel plain originate mainly from carbonate and silicate mineral dissolution.Strong adsorption by alternating cations has led to increased Na+concentration in downstream fine soil plains.Groundwater quality in northern alluvial fine soil plains as well as those around alluvial lacustrine fine soil plain within Kuitun River basin has exceeded China′s drinking water standard,and about 10% of the water samples have surpass irrigation water standards in future.These research findings can serve as a reference for rational utilization and management of groundwater resources within Kuitun River Basin.

    Key words:groundwater;hydrochemical characteristics;spatial distribution;genetic analysis;saturation index;Kuitun River Basin

    猜你喜歡
    水化學(xué)承壓水平原
    那一片平原
    黃河之聲(2022年6期)2022-08-26 06:48:50
    地鐵深基坑承壓水控制研究
    深層承壓水污染途徑及防治研究
    秦皇島北部地表水和地下水水化學(xué)特征及水質(zhì)評(píng)價(jià)
    平原的草
    浪起山走
    滇池(2019年1期)2019-02-14 02:36:54
    衡水市新近系館陶組地?zé)崴瘜W(xué)特征與分布規(guī)律研究
    筠連煤層氣井產(chǎn)出水化學(xué)特征及意義
    中國煤層氣(2015年5期)2015-08-22 03:26:32
    紅平原上的“小肉山”
    一種微生態(tài)制劑對(duì)鯽魚養(yǎng)殖水化學(xué)指標(biāo)的影響
    黑人巨大精品欧美一区二区蜜桃| 香蕉丝袜av| 久久久久久久久久久免费av| av免费观看日本| 高清av免费在线| 国精品久久久久久国模美| 国产精品久久久久成人av| 国产一级毛片在线| 国产成人一区二区在线| 人成视频在线观看免费观看| 久久久久国产一级毛片高清牌| 天美传媒精品一区二区| 丰满少妇做爰视频| 国产日韩欧美视频二区| 久久人人97超碰香蕉20202| 男人添女人高潮全过程视频| 少妇精品久久久久久久| av在线app专区| 丁香六月天网| 一个人免费看片子| 美女午夜性视频免费| 国产精品一区二区在线观看99| 中文欧美无线码| 国产激情久久老熟女| 天天躁日日躁夜夜躁夜夜| 三上悠亚av全集在线观看| 午夜福利,免费看| 母亲3免费完整高清在线观看| 最近最新中文字幕免费大全7| 国产免费视频播放在线视频| 天天躁夜夜躁狠狠躁躁| 在线精品无人区一区二区三| 免费观看a级毛片全部| av天堂久久9| 丝袜脚勾引网站| 国产成人系列免费观看| 男女免费视频国产| 一级毛片 在线播放| 男男h啪啪无遮挡| 国产成人欧美| 啦啦啦 在线观看视频| 欧美日韩成人在线一区二区| 精品第一国产精品| 如何舔出高潮| 日韩熟女老妇一区二区性免费视频| 制服诱惑二区| 啦啦啦视频在线资源免费观看| 国产精品一二三区在线看| 赤兔流量卡办理| 欧美精品一区二区大全| 男女高潮啪啪啪动态图| 少妇 在线观看| 国产精品香港三级国产av潘金莲 | 波多野结衣av一区二区av| 搡老岳熟女国产| 夫妻性生交免费视频一级片| 亚洲欧美成人综合另类久久久| 美女扒开内裤让男人捅视频| 国产精品免费大片| 不卡视频在线观看欧美| 国产免费现黄频在线看| 精品国产露脸久久av麻豆| 亚洲国产日韩一区二区| 新久久久久国产一级毛片| www.自偷自拍.com| 无遮挡黄片免费观看| 999精品在线视频| 18禁动态无遮挡网站| www.自偷自拍.com| 伦理电影大哥的女人| 久久99精品国语久久久| 少妇被粗大猛烈的视频| 亚洲伊人色综图| 黑人巨大精品欧美一区二区蜜桃| 女人高潮潮喷娇喘18禁视频| 欧美在线黄色| 黄色一级大片看看| 精品国产露脸久久av麻豆| 亚洲精品日韩在线中文字幕| 国产黄色免费在线视频| 亚洲美女黄色视频免费看| 一区在线观看完整版| 国产成人免费观看mmmm| 国产亚洲午夜精品一区二区久久| www.熟女人妻精品国产| 两个人看的免费小视频| 亚洲精品国产av成人精品| 日本91视频免费播放| 老司机影院毛片| 99热网站在线观看| 老司机亚洲免费影院| 国产精品蜜桃在线观看| 久热这里只有精品99| 91aial.com中文字幕在线观看| 青草久久国产| 十分钟在线观看高清视频www| 亚洲情色 制服丝袜| 婷婷色综合www| 女人被躁到高潮嗷嗷叫费观| 日韩欧美一区视频在线观看| 尾随美女入室| 亚洲精品美女久久av网站| 国产精品一区二区在线观看99| 亚洲av欧美aⅴ国产| 免费观看a级毛片全部| 久久久国产一区二区| 欧美日韩视频高清一区二区三区二| 国产精品一国产av| 国产高清国产精品国产三级| 中文字幕制服av| av有码第一页| 国产97色在线日韩免费| 五月天丁香电影| 欧美黑人精品巨大| 亚洲精品视频女| 成人漫画全彩无遮挡| 丝袜在线中文字幕| 狠狠婷婷综合久久久久久88av| xxx大片免费视频| av女优亚洲男人天堂| 少妇被粗大的猛进出69影院| 香蕉国产在线看| 久久99热这里只频精品6学生| 妹子高潮喷水视频| 国产精品av久久久久免费| 18禁裸乳无遮挡动漫免费视频| 亚洲四区av| 精品人妻在线不人妻| 亚洲一区二区三区欧美精品| 国产1区2区3区精品| 国产日韩欧美亚洲二区| 母亲3免费完整高清在线观看| 欧美少妇被猛烈插入视频| 亚洲熟女毛片儿| 亚洲久久久国产精品| 日本午夜av视频| 丝袜美足系列| h视频一区二区三区| 亚洲一区二区三区欧美精品| 丝袜美足系列| av视频免费观看在线观看| 亚洲一区二区三区欧美精品| 精品卡一卡二卡四卡免费| av在线老鸭窝| 国产极品天堂在线| 9191精品国产免费久久| 亚洲人成网站在线观看播放| 青草久久国产| 大片免费播放器 马上看| 丁香六月天网| 久久精品亚洲av国产电影网| 亚洲欧美一区二区三区久久| 另类亚洲欧美激情| 又大又黄又爽视频免费| 国产精品久久久久久精品电影小说| 国产精品久久久人人做人人爽| 大香蕉久久成人网| 黄色毛片三级朝国网站| 精品视频人人做人人爽| 久久99精品国语久久久| 一级毛片我不卡| 国产淫语在线视频| 亚洲美女黄色视频免费看| 丝袜美足系列| a级片在线免费高清观看视频| 观看av在线不卡| 亚洲欧美一区二区三区久久| 亚洲成人国产一区在线观看 | 男人操女人黄网站| 五月天丁香电影| 宅男免费午夜| 一区二区av电影网| 久热爱精品视频在线9| 高清视频免费观看一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久久久久久免费视频了| 中文字幕人妻丝袜制服| 亚洲av日韩精品久久久久久密 | 99热国产这里只有精品6| 欧美日韩视频精品一区| 少妇被粗大的猛进出69影院| 在线亚洲精品国产二区图片欧美| 久久狼人影院| 大香蕉久久网| 欧美日韩成人在线一区二区| 久久久精品94久久精品| 亚洲综合精品二区| 日韩制服丝袜自拍偷拍| 1024香蕉在线观看| 中文天堂在线官网| 亚洲成国产人片在线观看| 中文字幕av电影在线播放| 97精品久久久久久久久久精品| 日韩视频在线欧美| 国产av精品麻豆| 日韩伦理黄色片| 少妇猛男粗大的猛烈进出视频| 在线观看www视频免费| kizo精华| 国产精品 国内视频| 大码成人一级视频| 丝袜美足系列| 亚洲国产精品国产精品| 欧美最新免费一区二区三区| 久久热在线av| 搡老岳熟女国产| 91国产中文字幕| 嫩草影院入口| 亚洲国产精品999| 99精品久久久久人妻精品| 国产99久久九九免费精品| 国产日韩欧美视频二区| h视频一区二区三区| 最近最新中文字幕大全免费视频 | 十八禁人妻一区二区| 一级毛片电影观看| 水蜜桃什么品种好| 久久久久国产精品人妻一区二区| 亚洲av福利一区| 欧美国产精品一级二级三级| 91精品国产国语对白视频| 久久久久国产精品人妻一区二区| 国产xxxxx性猛交| 精品免费久久久久久久清纯 | 不卡视频在线观看欧美| 一级毛片电影观看| 国产在线一区二区三区精| 电影成人av| 91精品伊人久久大香线蕉| 久久婷婷青草| 亚洲美女视频黄频| 蜜桃国产av成人99| a级毛片在线看网站| 热re99久久国产66热| 考比视频在线观看| 国产亚洲最大av| 美女视频免费永久观看网站| 男女边摸边吃奶| 大片免费播放器 马上看| 另类亚洲欧美激情| 自拍欧美九色日韩亚洲蝌蚪91| 飞空精品影院首页| 欧美亚洲 丝袜 人妻 在线| av国产久精品久网站免费入址| 国产午夜精品一二区理论片| 亚洲av福利一区| 亚洲国产成人一精品久久久| 免费人妻精品一区二区三区视频| 久久精品国产综合久久久| 国产一卡二卡三卡精品 | 妹子高潮喷水视频| 日韩一本色道免费dvd| 久久国产亚洲av麻豆专区| 成年动漫av网址| 午夜老司机福利片| 亚洲欧美一区二区三区国产| 你懂的网址亚洲精品在线观看| 丰满迷人的少妇在线观看| 欧美av亚洲av综合av国产av | 国产日韩欧美亚洲二区| 久久精品久久久久久噜噜老黄| 最近2019中文字幕mv第一页| 亚洲七黄色美女视频| 中文字幕av电影在线播放| 成年动漫av网址| 菩萨蛮人人尽说江南好唐韦庄| 秋霞在线观看毛片| 亚洲一码二码三码区别大吗| netflix在线观看网站| 亚洲国产精品999| 大码成人一级视频| 国产亚洲av片在线观看秒播厂| 大香蕉久久成人网| 免费黄色在线免费观看| 男的添女的下面高潮视频| 亚洲国产最新在线播放| av有码第一页| 最近中文字幕高清免费大全6| 日韩电影二区| 各种免费的搞黄视频| 国产欧美日韩一区二区三区在线| 国产黄频视频在线观看| 最近手机中文字幕大全| 久久久国产精品麻豆| 亚洲av电影在线观看一区二区三区| 国产野战对白在线观看| 亚洲欧洲国产日韩| 国产成人欧美| 亚洲美女黄色视频免费看| a级毛片黄视频| 国产男女超爽视频在线观看| 久久久久久久久免费视频了| 国产在线一区二区三区精| 国产成人精品久久二区二区91 | 在线观看www视频免费| 黄色一级大片看看| 亚洲精品,欧美精品| 秋霞伦理黄片| 亚洲精品自拍成人| 亚洲婷婷狠狠爱综合网| 久久婷婷青草| 亚洲熟女精品中文字幕| 人妻人人澡人人爽人人| 赤兔流量卡办理| 亚洲国产av新网站| 日本欧美国产在线视频| 99re6热这里在线精品视频| 哪个播放器可以免费观看大片| 久久韩国三级中文字幕| 美女中出高潮动态图| 成人黄色视频免费在线看| 丁香六月欧美| 国产深夜福利视频在线观看| 制服丝袜香蕉在线| 中国国产av一级| 男人舔女人的私密视频| 国产av码专区亚洲av| 国产男人的电影天堂91| 国产乱来视频区| 精品酒店卫生间| 中国三级夫妇交换| 国产人伦9x9x在线观看| 成人18禁高潮啪啪吃奶动态图| 亚洲五月色婷婷综合| 人妻 亚洲 视频| 久久99精品国语久久久| 久久97久久精品| 丁香六月天网| 日韩av不卡免费在线播放| 啦啦啦中文免费视频观看日本| 午夜福利视频在线观看免费| 在线观看免费午夜福利视频| 欧美成人精品欧美一级黄| 最新的欧美精品一区二区| 亚洲欧美中文字幕日韩二区| 天天操日日干夜夜撸| 欧美人与性动交α欧美软件| 狠狠婷婷综合久久久久久88av| 人人澡人人妻人| 青春草国产在线视频| 卡戴珊不雅视频在线播放| 欧美亚洲日本最大视频资源| 高清视频免费观看一区二区| 91国产中文字幕| 色播在线永久视频| 免费黄频网站在线观看国产| 日韩人妻精品一区2区三区| 色视频在线一区二区三区| 久久久国产欧美日韩av| 亚洲av欧美aⅴ国产| 看免费av毛片| 女的被弄到高潮叫床怎么办| 99re6热这里在线精品视频| 欧美精品人与动牲交sv欧美| www.熟女人妻精品国产| 超碰97精品在线观看| 免费高清在线观看日韩| 国产精品人妻久久久影院| 最近2019中文字幕mv第一页| 亚洲精品美女久久av网站| 99热网站在线观看| 黑人巨大精品欧美一区二区蜜桃| 日本wwww免费看| 免费黄网站久久成人精品| 美女主播在线视频| 亚洲国产成人一精品久久久| bbb黄色大片| 狠狠精品人妻久久久久久综合| 日本av免费视频播放| 激情视频va一区二区三区| 亚洲av国产av综合av卡| 建设人人有责人人尽责人人享有的| 日韩av在线免费看完整版不卡| 欧美 亚洲 国产 日韩一| 一级毛片黄色毛片免费观看视频| 午夜日韩欧美国产| 黑丝袜美女国产一区| 一二三四中文在线观看免费高清| 国产精品一区二区在线观看99| 中文字幕高清在线视频| 国产成人精品在线电影| 国产在视频线精品| 欧美xxⅹ黑人| 午夜激情av网站| 一级片免费观看大全| 国产精品99久久99久久久不卡 | 99久久99久久久精品蜜桃| 久久久国产精品麻豆| 精品少妇内射三级| 亚洲人成77777在线视频| 亚洲婷婷狠狠爱综合网| 欧美激情高清一区二区三区 | 丰满乱子伦码专区| 成年动漫av网址| 欧美中文综合在线视频| 1024香蕉在线观看| 人人妻人人添人人爽欧美一区卜| 国产一区二区激情短视频 | 激情五月婷婷亚洲| 赤兔流量卡办理| 中文字幕最新亚洲高清| 精品少妇黑人巨大在线播放| 美女福利国产在线| 叶爱在线成人免费视频播放| 国产精品久久久久成人av| 大话2 男鬼变身卡| 国产精品免费视频内射| 国产女主播在线喷水免费视频网站| 在现免费观看毛片| 亚洲激情五月婷婷啪啪| 免费看av在线观看网站| 欧美日韩福利视频一区二区| 亚洲av在线观看美女高潮| 精品一区二区三区四区五区乱码 | avwww免费| 美女扒开内裤让男人捅视频| 亚洲精品一二三| 人人澡人人妻人| 国产在线视频一区二区| 好男人视频免费观看在线| 国产精品av久久久久免费| 免费在线观看视频国产中文字幕亚洲 | 亚洲精品乱久久久久久| 日韩电影二区| 成人18禁高潮啪啪吃奶动态图| 欧美人与善性xxx| 我的亚洲天堂| 国产免费现黄频在线看| 免费在线观看视频国产中文字幕亚洲 | 丝瓜视频免费看黄片| 精品一区二区免费观看| 一个人免费看片子| 多毛熟女@视频| 欧美激情 高清一区二区三区| 黄片播放在线免费| 久久久久久久久久久免费av| 老司机深夜福利视频在线观看 | 亚洲av电影在线进入| 国产在视频线精品| 母亲3免费完整高清在线观看| 一级黄片播放器| 91aial.com中文字幕在线观看| 在线亚洲精品国产二区图片欧美| 七月丁香在线播放| 一区福利在线观看| 青青草视频在线视频观看| 国产一区二区激情短视频 | 欧美久久黑人一区二区| 啦啦啦 在线观看视频| 亚洲婷婷狠狠爱综合网| 久久99一区二区三区| 水蜜桃什么品种好| 成人黄色视频免费在线看| 丝瓜视频免费看黄片| 少妇的丰满在线观看| 精品一区二区三卡| 久久热在线av| 90打野战视频偷拍视频| av不卡在线播放| 纵有疾风起免费观看全集完整版| 韩国高清视频一区二区三区| 亚洲av成人精品一二三区| 波多野结衣av一区二区av| 亚洲欧洲日产国产| 精品少妇黑人巨大在线播放| 一级毛片我不卡| 国产成人精品久久二区二区91 | 大话2 男鬼变身卡| 亚洲欧美激情在线| 少妇被粗大猛烈的视频| 午夜影院在线不卡| 亚洲av电影在线观看一区二区三区| 美女福利国产在线| 国产亚洲欧美精品永久| 18禁动态无遮挡网站| 最近最新中文字幕大全免费视频 | 美女扒开内裤让男人捅视频| 亚洲少妇的诱惑av| 女人精品久久久久毛片| 赤兔流量卡办理| 最近中文字幕高清免费大全6| 免费高清在线观看日韩| 国产精品一区二区精品视频观看| 亚洲一卡2卡3卡4卡5卡精品中文| 精品亚洲成国产av| 国产精品欧美亚洲77777| 久久精品国产亚洲av涩爱| 亚洲国产最新在线播放| 色精品久久人妻99蜜桃| 999久久久国产精品视频| 日韩中文字幕视频在线看片| 欧美激情 高清一区二区三区| 卡戴珊不雅视频在线播放| 国产成人精品久久久久久| 777久久人妻少妇嫩草av网站| 欧美日韩视频精品一区| 国产亚洲av片在线观看秒播厂| 亚洲欧洲日产国产| 看免费成人av毛片| 在现免费观看毛片| 香蕉国产在线看| 极品少妇高潮喷水抽搐| 男女边吃奶边做爰视频| 成年人免费黄色播放视频| 一区二区三区乱码不卡18| 美女大奶头黄色视频| 亚洲国产欧美日韩在线播放| 十分钟在线观看高清视频www| 国产午夜精品一二区理论片| 国产成人a∨麻豆精品| 午夜福利网站1000一区二区三区| 亚洲精品在线美女| e午夜精品久久久久久久| 一区二区日韩欧美中文字幕| 黄色毛片三级朝国网站| 交换朋友夫妻互换小说| 久久久久久久大尺度免费视频| 少妇人妻久久综合中文| 久久毛片免费看一区二区三区| 国产1区2区3区精品| 国产成人精品无人区| 青春草国产在线视频| 精品久久久久久电影网| 欧美黑人精品巨大| 亚洲色图综合在线观看| 亚洲欧美精品自产自拍| 免费看av在线观看网站| 可以免费在线观看a视频的电影网站 | 少妇人妻久久综合中文| 国产精品国产av在线观看| 毛片一级片免费看久久久久| 日韩人妻精品一区2区三区| 久久久久久久国产电影| 欧美中文综合在线视频| 国产精品.久久久| 久久久久精品性色| 99热国产这里只有精品6| 男女无遮挡免费网站观看| xxxhd国产人妻xxx| 大话2 男鬼变身卡| 女人被躁到高潮嗷嗷叫费观| 国产亚洲欧美精品永久| 黄色一级大片看看| 亚洲国产精品成人久久小说| 欧美日韩福利视频一区二区| 国产精品香港三级国产av潘金莲 | 狠狠精品人妻久久久久久综合| 爱豆传媒免费全集在线观看| 国产精品二区激情视频| 亚洲欧美清纯卡通| 日本av免费视频播放| 国产成人系列免费观看| 精品少妇久久久久久888优播| 国产成人精品久久久久久| 国产成人欧美| 高清av免费在线| 日韩,欧美,国产一区二区三区| 在线 av 中文字幕| 久久精品亚洲av国产电影网| av在线老鸭窝| 日本爱情动作片www.在线观看| 国产av精品麻豆| 欧美国产精品va在线观看不卡| 亚洲精品aⅴ在线观看| 亚洲精品成人av观看孕妇| 街头女战士在线观看网站| 成年av动漫网址| 日韩大片免费观看网站| 精品国产一区二区三区四区第35| 亚洲av福利一区| 日韩一区二区视频免费看| √禁漫天堂资源中文www| 亚洲第一青青草原| 久久99热这里只频精品6学生| 不卡视频在线观看欧美| 久久精品久久精品一区二区三区| e午夜精品久久久久久久| 日本91视频免费播放| 国产片内射在线| 国产男女超爽视频在线观看| 一区二区三区激情视频| 狂野欧美激情性xxxx| 久久鲁丝午夜福利片| 亚洲国产精品成人久久小说| 日韩不卡一区二区三区视频在线| 国产av一区二区精品久久| 欧美日韩国产mv在线观看视频| 亚洲一区二区三区欧美精品| 丁香六月天网| 久久99热这里只频精品6学生| 男的添女的下面高潮视频| www.自偷自拍.com| 在线观看国产h片| 精品免费久久久久久久清纯 | 久久精品国产亚洲av涩爱| 一边摸一边抽搐一进一出视频| 少妇人妻 视频| 久久亚洲国产成人精品v| 热99国产精品久久久久久7| 免费在线观看黄色视频的| 亚洲人成网站在线观看播放| 亚洲综合精品二区| 国产一区二区三区av在线| 成年动漫av网址| 亚洲国产最新在线播放| 国产精品一国产av| 久久婷婷青草| 亚洲国产最新在线播放| 亚洲熟女精品中文字幕| 久热这里只有精品99| netflix在线观看网站| 亚洲美女黄色视频免费看| 精品一品国产午夜福利视频| 欧美日韩国产mv在线观看视频| 男人操女人黄网站|