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

    黃河下游影響帶(河南段)水化學(xué)演化規(guī)律研究

    2021-12-17 00:46:36孫龍王莉莉曹文庚
    人民黃河 2021年12期
    關(guān)鍵詞:黃河水水化學(xué)黃河

    孫龍 王莉莉 曹文庚

    摘 要:黃河下游地區(qū)(河南段)是河南省重要的工業(yè)和糧食基地,該區(qū)域地下水化學(xué)演化受自然因素和人類活動的交互影響,而地下水與地表水的綜合研究有助于揭示地下水成因及其水文地球化學(xué)過程。通過在黃河下游影響帶(河南段)采集黃河水和地下水樣,分析黃河水和地下水的水化學(xué)特征演化規(guī)律及二者之間的聯(lián)系,結(jié)果表明:研究區(qū)黃河水pH值和TDS(總?cè)芙夤腆w)均表現(xiàn)出季節(jié)性差異,地下水pH值為7.30~8.10,總體上低于河水;研究區(qū)黃河水側(cè)滲轉(zhuǎn)化為地下水,水化學(xué)類型由SO4·Cl·HCO3-Na·Ca·Mg、SO4·HCO3·Cl-Na·Ca·Mg、HCO3·SO4·Cl-Na·Ca·Mg型演變?yōu)镠CO3·Ca-Mg·Na、HCO3-Mg·Ca·Na型;研究區(qū)黃河水和地下水均受巖石風(fēng)化和蒸發(fā)濃縮作用的影響,地下水中碳酸廣泛參與碳酸鹽礦物的溶解;研究區(qū)黃河水Cl-、SO2-4和NO-3質(zhì)量濃度較高可能與上游景區(qū)排污有關(guān),地下水個別取樣點陰離子Cl-質(zhì)量濃度較高可能是農(nóng)田施用含Cl-農(nóng)藥或化肥污染地下水所致。

    關(guān)鍵詞:水化學(xué)特征;黃河水;地下水;黃河下游

    中圖分類號:TV62;TV882.1 ? 文獻標志碼:A

    doi:10.3969/j.issn.1000-1379.2021.12.018

    引用格式:孫龍,王莉莉,曹文庚.黃河下游影響帶(河南段)水化學(xué)演化規(guī)律研究[J].人民黃河,2021,43(12):91-99.

    Abstract: The lower reaches of the Yellow River (Henan section) is an important industrial and grain base in Henan Province. The hydrochemical evolution of groundwater in this region is influenced by the interaction of natural factors and strong human activities. The comprehensive study of groundwater and surface water is helpful to reveal the genesis of groundwater and the hydrogeochemical process. Water samples of the Yellow River and underground water were collected in the influence zone of the lower Yellow River in Henan to study the hydrochemical characteristics evolution of the Yellow River and underground water and the relationship between them. The results show that a)the pH value and TDS of the Yellow River water in the study area show seasonal differences. The pH range of groundwater is 7.3~8.1, which is generally lower than that of the river water and generally consistent with the seasonal differences of the Yellow River water in the same period. b)With the lateral flow from the Yellow River to groundwater in the influence zone of the Yellow River, the hydrochemical types are changed from SO4·Cl·HCO3-Na·Ca·Mg, SO4·HCO3·Cl-Na·Ca·Mg and HCO3·SO4·Cl-Na·Ca·Mg to HCO3·Ca-Mg·Na and HCO3-Mg·Ca·Na. c)Both the Yellow River water and groundwater in the study area are affected by rock weathering and evaporation and concentration, and carbonic acid in the groundwater is widely involved in the dissolution of carbonate minerals. d)Cl-, SO2-4 and NO-3 of the Yellow River in the study area are relatively higher, which may be related to the discharge of the upstream scenic spots. The high level of anion Cl- in some points of groundwater may be caused by the pollution of groundwater by cl-containing pesticides or chemical fertilizers of farmland.

    Key words: hydrochemical characteristics; Yellow River water; groundwater; Lower Yellow River

    在自然系統(tǒng)的循環(huán)過程中,地下水與其接觸的巖石圈、生物圈和大氣圈進行著極其復(fù)雜的物質(zhì)、能量和信息交換,水化學(xué)特征在不斷變化[1]。在全球氣候變化的大背景下,人類活動往往在較短時間使在自然條件下經(jīng)過漫長地質(zhì)演化形成的穩(wěn)定地下水化學(xué)環(huán)境發(fā)生劇烈變化[2-4]。由于水體化學(xué)組分記錄了水體形成、循環(huán)和演化的歷史[5],因此地下水化學(xué)時空變異特征與演變規(guī)律可以揭示地下水與環(huán)境的相互作用機制[6]。

    自然條件變化和經(jīng)濟發(fā)展速度不同導(dǎo)致各地區(qū)地下水化學(xué)特征演化的主控因素不同:河南平原淺層地下水豐枯期水化學(xué)特征受氣候條件(降雨、蒸發(fā)等)、地下水循環(huán)條件、含水介質(zhì)以及人類活動等因素共同影響[7];黃河三角洲淺層地下水化學(xué)特征演化的驅(qū)動因素是黃河入海流路的變遷和海水入侵[8];城市化和農(nóng)村工業(yè)化過程中東莞地區(qū)地下水化學(xué)特征演化的主控因素為海水入侵、水巖相互作用和河流側(cè)向流動等自然過程以及工業(yè)廢水排放和農(nóng)業(yè)面源污染等人類活動[9];2001年生態(tài)輸水以來,內(nèi)蒙古西部額濟納綠洲淺層地下水化學(xué)特征與水位埋深有一定響應(yīng)關(guān)系[10];石家莊淺層地下水水質(zhì)演化與水文地質(zhì)條件、污染物排放及地下水開采等因素有關(guān)[11];河西走廊疏勒河流域地下水化學(xué)特征演化的主控因素是巖石風(fēng)化[12]。一般情況下,在人類活動影響較小的地區(qū)地下水化學(xué)特征主要受控于區(qū)域水文地質(zhì)條件、蒸散發(fā)等自然因素,在人口密度較大及人類活動較為頻繁的地區(qū)地下水化學(xué)特征主要受人類生活、工農(nóng)業(yè)排污等人為因素和蒸發(fā)濃縮、溶濾等自然因素共同影響,且人為因素占主導(dǎo)地位。

    黃河下游地區(qū)(河南段)是河南省重要的工業(yè)和糧食基地,供水格局以引黃河水為主、開采地下水為輔[13]。近20 a來隨著經(jīng)濟的發(fā)展和人類需求的增加,黃河中下游(河南段)人類活動增強,導(dǎo)致該區(qū)域地下水化學(xué)特征演化受到了自然因素和強烈人類活動的交互影響。然而,目前對黃河影響帶(河南段)的研究主要集中在地表水和地下水資源開發(fā)及其影響方面[14-15],對于變化條件下黃河下游影響帶(河南段)地下水化學(xué)特征、水化學(xué)類型演變規(guī)律的研究較少。地下水與地表水的綜合研究有助于揭示地下水成因及其水文地球化學(xué)過程[16],因此以黃河下游影響帶(河南段)為研究區(qū),運用水化學(xué)分析法、Piper三線圖、Gibbs圖和離子比例系數(shù)分析研究區(qū)內(nèi)黃河水和地下水的水化學(xué)特征及其演化規(guī)律,并在此基礎(chǔ)上分析其控制因素,以期為研究區(qū)水循環(huán)研究及地下水可持續(xù)開發(fā)提供參考。

    1 研究區(qū)概況

    研究區(qū)位于河南省北部,西起桃花峪,東至濮陽市臺前縣(黃河沿程長240 km左右),以現(xiàn)代黃河為軸,向兩側(cè)擴展數(shù)千米至數(shù)十千米——黃河影響帶范圍[14,17],涉及焦作、新鄉(xiāng)、濮陽、鄭州及開封共5市9個區(qū)縣,面積約1萬km2。地面高程40~85 m,黃河以北地勢向東北緩傾斜,黃河以南地勢向東南緩傾斜。研究區(qū)屬暖溫帶季風(fēng)氣候區(qū),四季分明,總的氣候特征為春季干旱風(fēng)沙多,夏熱多雨,秋季晴和日照長,冬季寒冷雪少,年內(nèi)降水分配不均,降水集中在6—9月。黃河下游平原廣泛分布新生代新近紀和第四紀松散堆積物,松散層中夾有較多各類砂層,這些砂層構(gòu)成研究區(qū)主要含水層,賦存有較豐富的地下水資源。淺層地下水流場以黃河為中心軸,呈放射狀流向東南和東北,地下水接受黃河側(cè)滲補給,且補給條件優(yōu)越,含水層富水性好,具有較大的開發(fā)潛力。含水層在垂向上呈由細至粗(粉土—粉砂—粉細砂—細中砂—中粗砂)的變化特點,砂層厚度為40~100 m,含砂比50%以上,個別地段可達80%[18]。開封—封丘以西為黃河沖積扇,以東為黃河泛濫平原,受黃河擺動時間不同影響,主流帶和泛流帶相間分布。沖積扇由灰、灰黑、黑灰色粉土,粉質(zhì)黏土與厚層粉細砂、細粉砂、局部中細砂組成,形成一較厚的具二元結(jié)構(gòu)的旋回層。研究區(qū)典型水文地質(zhì)剖面見圖1,其中S4、ZK34等為鉆孔編號。

    黃河下游,河床一般高出堤外平原3~5 m,最高達10 m,成為黃河下游沖積平原的分水嶺[15]。研究區(qū)淺層地下水主要接受大氣降水、河渠入滲、灌溉回滲及黃河側(cè)滲等補給,潛水蒸發(fā)和人工開采為主要排泄方式。天然條件下,研究區(qū)淺層地下水動態(tài)主要受水文、氣象因素制約,隨著距黃河水邊線距離的增加,淺層地下水受黃河影響由強漸弱,由強影響帶到弱影響帶,逐漸過渡到未影響帶,而氣象因素由次要地位漸變?yōu)橹鲗?dǎo)地位[18]。研究區(qū)農(nóng)業(yè)活動廣泛,主要種植玉米、小麥、花生等農(nóng)作物,部分區(qū)域有蘋果、桃子等經(jīng)濟果樹,農(nóng)業(yè)以引黃灌溉和開采淺層地下水井灌為主。此外,沿黃兩岸地區(qū)淺層地下水是沿黃河兩岸大中城市(河南段)的重要供水水源。

    2 數(shù)據(jù)和方法

    2.1 樣品采集

    結(jié)合研究區(qū)實際情況,在桃花峪至河南與山東省界的黃河河段區(qū)間垂直于黃河布設(shè)4個剖面,各剖面分別布設(shè)黃河水和地下水取樣點,其中武陟剖面取樣點編號為WZHH、WZ1、WZ2、WZ3、WZ4,原陽剖面取樣點編號為YYHH、YY1、YY2、YY3,封丘剖面取樣點編號為FQHH、FQ1、FQ2,濮陽剖面取樣點編號為PYHH、PY1、PY2、PY3、LK。2017年7月底(潛水豐水期)、2017年11月底(潛水平水期)和2018年3月底(潛水枯水期)分別采集地表水樣和地下水樣,每期采集黃河水樣4個、地下水樣13個,共計3期51個,每個樣品均采集平行樣。每期采樣在2~3 d內(nèi)完成,同時避開降雨天。

    2.2 樣品測試

    樣品測試在河南省地質(zhì)環(huán)境監(jiān)測院實驗測試中心進行,主要測試項目包括pH值、化學(xué)需氧量(COD)、溶解性總固體(TDS)、總硬度,以及Na+、K+、Ca2+、Mg2+、Cl-、SO2-4、NO-3、HCO-3、CO2-3含量等指標,其中Na+、K+含量采用火焰原子吸收分光光度法測定,Ca2+、Mg2+含量采用EDTA滴定法測定,Cl-、SO2-4、NO-3含量采用離子色譜儀測定,HCO-3和CO2-3含量采用酸堿滴定法測定。

    3 結(jié)果與分析

    黃河下游影響帶(河南段)黃河水、地下水主要化學(xué)指標統(tǒng)計結(jié)果見表1。

    3.1 pH值與TDS

    研究區(qū)所有水樣pH值均大于7.0小于8.5,呈弱堿性。黃河水的pH值為7.4~8.1(見圖2),平均值為7.85,與蘇小四等(7.49~8.18)[19]、陳靜生等(7.9~8.0)[20]及Zhang等(7.05~8.71)[21]的研究結(jié)果相近。其中2017年7月和2017年11月黃河水pH值在取樣點YYHH和WZHH有較大波動,而2018年3月黃河水pH值變化范圍較小,整體上黃河水2018年3月pH值>2017年11月pH值>2017年7月pH值,表現(xiàn)出顯著的季節(jié)性變化特征。地下水的pH值總體上低于河水的,為7.30~8.10,平均值為7.63,接近于黃河影響帶北岸新鄉(xiāng)段2008年10月地下水的pH值(7.34~7.73)[22],地下水pH值大體表現(xiàn)出與同期黃河水一致的變化趨勢,即2018年3月pH值>2017年11月pH值>2017年7月pH值,另外地下水采樣點pH值年內(nèi)波動范圍為0.1~0.4。

    研究區(qū)黃河水TDS為520.55~719.87 mg/L,均值為636.83 mg/L,各期黃河水的TDS沒有明顯波動,年際2017年7月TDS>2018年3月TDS>2017年11月TDS,季節(jié)性差異明顯。地下水TDS為438.80~1 250.40 mg/L,均值為646.66 mg/L,各采樣點TDS年內(nèi)變化不大。除了2018年3月的LK取樣點外,地下水TDS均低于1 000 mg/L,達到《地下水質(zhì)量標準》(GB/T 14848—2017)Ⅲ類水標準。

    3.2 水化學(xué)成分變化特征

    研究區(qū)黃河水標準無機離子平衡常數(shù)(NIBC)為8.64%~19.72%,而黃河所有河段河水NIBC為-23%~21%[23],說明本研究所得水質(zhì)分析結(jié)果可靠。研究區(qū)黃河水的主要陽離子質(zhì)量濃度Na+>Ca2+>Mg2+>K+(見圖2),優(yōu)勢陽離子為Na+,其質(zhì)量濃度變化范圍為82.13~145.00 mg/L、均值為111.99 mg/L;主要陰離子質(zhì)量濃度分幾種情況:2017年7月SO2-4>HCO-3>Cl->NO-3,2017年11月HCO-3和SO2-4質(zhì)量濃度相當(dāng),2018年3月HCO-3>SO2-4>Cl->NO-3,HCO-3和SO2-4質(zhì)量濃度各季節(jié)略有差異。黃河水各離子質(zhì)量濃度與TDS變化相對應(yīng),在研究河段沒有顯示出沿程增大的特征,同時河水各離子的標準差和變異系數(shù)均較小,表明黃河水中各離子質(zhì)量濃度變化相對較小。

    跟黃河水一樣,地下水各陽離子質(zhì)量濃度大小也為Na+>Ca2+>Mg2+>K+,優(yōu)勢陽離子為Na+和Ca2+,其質(zhì)量濃度分別為40.67~396.56 mg/L、30.26~129.86 mg/L;各陰離子質(zhì)量濃度大小為HCO-3>SO2-4>Cl->NO-3,個別取樣點Cl->SO2-4,HCO-3在陰離子中占絕對優(yōu)勢,變化范圍在338.66~762.14 mg/L,均值為507.36 mg/L。地下水中除了NO-3外,各離子變異系數(shù)均小于1,但HCO-3、SO2-4、Cl-、Na+和Ca2+的平均值和標準差都較大,說明他們在地下水中絕對含量較大,但相對含量差別不大。

    總體上,地下水中各離子質(zhì)量濃度的波動幅度大于河水的,且不同離子質(zhì)量濃度與河水的差別較大(見圖3,WZ、YY、FQ、PY分別表示武陟、原陽、封丘、濮陽斷面)。HCO-3表現(xiàn)為地下水質(zhì)量濃度遠高于黃河水,地下水中HCO-3質(zhì)量濃度是黃河水的1.47~4.54倍,且總體上與黃河距離越近其質(zhì)量濃度越低;Mg2+也表現(xiàn)出了相同特征,即ρ(Mg2+)地下水>ρ(Mg2+)河水(ρ表示質(zhì)量濃度);黃河水NO-3、SO2-4、Cl-和K+質(zhì)量濃度高于地下水的,其中NO-3質(zhì)量濃度差別最大,河水中NO-3質(zhì)量濃度(16.78~21.92 mg/L)遠大于地下水的(<0.01~5.70 mg/L),說明黃河水相對于地下水受到了更多人類活動的影響。

    3.3 水化學(xué)類型

    水化學(xué)Piper圖可直觀揭示水體中主要離子的空間關(guān)系,顯示研究區(qū)內(nèi)地下水化學(xué)成分的空間演變規(guī)律[24]。研究區(qū)黃河水和地下水Piper三線圖見圖4。

    從圖4可以看出,黃河水陰離子較為集中地落在了4區(qū),表明黃河水中強酸根SO2-4+Cl-質(zhì)量濃度大于弱酸根HCO-3,同時陽離子多數(shù)分布在2區(qū),說明堿金屬離子Na+質(zhì)量濃度大于堿土金屬離子Ca2++Mg2+,其中Na+占比46.61%~57.97%,因此黃河水主要水化學(xué)類型為SO4·Cl·HCO3-Na·Ca·Mg、SO4·HCO3·Cl-Na·Ca·Mg和HCO3·SO4·Cl-Na·Ca·Mg。地下水陰離子主要落在3區(qū),碳酸鹽硬度超過55%,陽離子主要落在了1區(qū),Ca2++Mg2+為主要陽離子,地下水化學(xué)類型主要為HCO3-Ca·Mg·Na和HCO3-Mg·Ca·Na。另外,LK取樣點陽離子落在2區(qū),水化學(xué)類型為HCO3-Na·Mg;PY3取樣點陰離子落在了4區(qū),水化學(xué)類型為HCO3·SO4·Cl-Ca·Mg。

    4 討 論

    4.1 地表水與地下水的關(guān)系

    (1)黃河水對地下水的補給。黃河水pH值與TDS均表現(xiàn)出明顯的季節(jié)性差異,這可能與不同季節(jié)各種來源的水對黃河水的補給貢獻不同有關(guān)[25]。研究區(qū)黃河水均為上中游來水,每年3月中旬開始黃河源頭氣溫回升至零度以上,上游冰雪融水[21]和凍土層融化水[25]大量補給黃河水,使得黃河水pH值相對較高。4月中旬以后冰雪融水補給貢獻減少、地下水補給貢獻增大,黃河中下游大量引黃灌區(qū)的灌溉回歸水也會補給黃河水,綜合作用下導(dǎo)致黃河水pH值表現(xiàn)為季節(jié)性變化。地下水pH值大體表現(xiàn)出與同期黃河水一致的變化趨勢,可能原因是,研究區(qū)屬于黃河懸河段,黃河側(cè)滲為研究區(qū)淺層地下水重要補給來源,其占地下水總補給量的11.23%[26]~12.86%[27],同時傍河水源地地下水大量開采等加速了地下水循環(huán)交替。此外,研究區(qū)地下水補給源的多樣性導(dǎo)致距黃河不同距離的地下水pH值波動幅度不同。

    研究區(qū)黃河水TDS均值為636.83 mg/L,高于蘇小四等[19]、何姜毅等[23]的研究結(jié)果,這與黃河流域近年來降水量減少以及農(nóng)田灌溉退水、工業(yè)廢水和生活污水排放等人為輸入有關(guān)。地下水TDS為438.80~1 250.40 mg/L,與距黃河3~9 km的原陽縣灘區(qū)地下水TDS(486.0~1 572.7 mg/L)[22]較為接近。研究區(qū)地下水各采樣點TDS年內(nèi)變化不大,主要原因是TDS受黃河側(cè)滲補給和地下水蒸發(fā)濃縮等多重影響。

    (2)水體轉(zhuǎn)化對應(yīng)的水化學(xué)類型演變趨勢。地下水主要離子的空間變化可以反映區(qū)域含水層的異質(zhì)性和連通性[28],揭示主導(dǎo)研究區(qū)地下水演化的機制[29]。黃河水側(cè)滲補給黃河影響帶的地下水,使得水體從地表水轉(zhuǎn)化為地下水,因此沿著地下水徑流方向距離黃河越近的點,其地下水與黃河水的水化學(xué)特征越相似。由圖4可以看出,大多數(shù)地下水陽離子聚集在左下方三角形中間位置,與黃河水陽離子所在區(qū)域極為接近;右下方三角形中大多數(shù)地下水陰離子落在了左邊重碳酸根的位置,而黃河水陰離子則落在地下水略右上方的三角形中間位置。黃河水轉(zhuǎn)為地下水,陽離子Ca2+和Mg2+占比略有增加、Na+占比減少,陰離子HCO-3占比明顯增加、Cl-和SO2-4占比減少,水化學(xué)類型由SO4·Cl·HCO3-Na·Ca·Mg、SO4·HCO3·Cl-Na·Ca·Mg和HCO3·SO4·Cl-Na·Ca·Mg型演變?yōu)镠CO3·Ca-Mg·Na和HCO3-Mg·Ca·Na型。從黃河水到地下水離子質(zhì)量濃度變化可以初步推斷,研究區(qū)地下水徑流過程中存在碳酸鹽礦物溶解[30],比如方解石和白云巖。此外,地下水PY3、LK取樣點與大多數(shù)地下水取樣點有不同的水化學(xué)特征,表明這2個地下水取樣點與黃河水之間的聯(lián)系較弱。

    4.2 水化學(xué)組分控制因素分析

    (1)水巖作用。根據(jù)γ(Na+)/γ(Na++Ca2+)或γ(Cl-)/γ(Cl-+HCO3-)的值(γ表示毫克當(dāng)量濃度)可以定性判斷區(qū)域巖石風(fēng)化、大氣降水、蒸發(fā)濃縮作用等對水體水化學(xué)的影響[31]。由研究區(qū)黃河水和地下水的Gibbs圖(見圖5)可知,黃河水γ(Na+)/γ(Na++Ca2+)值為0.57~0.71,位于巖石風(fēng)化帶和蒸發(fā)濃縮作用帶的過渡帶,且遠離降水控制區(qū),反映了研究區(qū)黃河水主要組分受巖石風(fēng)化和蒸發(fā)濃縮作用共同影響,這與張旺等[32]的研究結(jié)果基本一致。地下水γ(Na+)/γ(Na++Ca2+)值橫跨范圍較大,從巖石風(fēng)化帶到巖石風(fēng)化和蒸發(fā)濃縮作用的過渡帶再到蒸發(fā)濃縮作用帶均有分布,說明地下水主要組分從主要受巖石風(fēng)化作用到巖石風(fēng)化和蒸發(fā)濃縮二者共同作用,再到主要受蒸發(fā)濃縮作用,具體主控因素由地下水采樣點所處位置決定:最左側(cè)地下水PY3取樣點距離圖5中黃河水取樣點較遠,與黃河的實際距離也較遠,主要受巖石風(fēng)化作用;中間大部分地下水取樣點越接近黃河,與黃河水γ(Na+)/γ(Na++Ca2+)值越相近,說明隨著接近黃河,地下水位埋深變淺,巖石風(fēng)化作用逐漸減弱,地下水受巖石風(fēng)化和蒸發(fā)濃縮共同作用;地下水LK取樣點在黃河水取樣點右側(cè)的蒸發(fā)濃縮帶,與大多數(shù)地下水取樣點距離較遠,說明該點與大多數(shù)地下水取樣點水化學(xué)成因不同,該取樣點位于黃河右岸距黃河垂直距離為2.8 km的蘭考縣東壩頭鄉(xiāng),該區(qū)域地下水成因較為復(fù)雜,黃河水對當(dāng)?shù)厮瘜W(xué)影響力較弱。

    地下水中各組分含量之間的比例系數(shù)常常被用來研究某些水文地球化學(xué)問題[33],γ(HCO-3)/γ(Cl-)可以反映陰離子沿徑流途徑的水文地球化學(xué)過程[34],地下水成因系數(shù)γ(Na+)/γ(Cl-)是表征地下水中鈉離子富集程度的一個水文地球化學(xué)參數(shù)[1],可以揭示地下水系統(tǒng)中巖鹽的來源[35]。研究區(qū)地下水γ(HCO-3)/γ(Cl-)為2.15~9.28(見圖6),說明研究區(qū)地下水中存在方解石、白云石礦物溶解。地下水Na+和Cl-主要來源于含NaCl巖鹽的溶解[36],圖6(b)中地下水點據(jù)幾乎都位于1∶1直線以上,說明Na+濃度大多大于Cl-濃度,因此Na+除了來源于巖鹽溶解外還可能來源于與Ca2+、Mg2+發(fā)生陽離子交換作用[37]。如果地下水演化主要受蒸發(fā)作用的影響而不受水巖交互作用,則γ(Na+)/γ(Cl-)值保持恒定[38],圖6(b)顯示γ(Na+)/γ(Cl-)隨著Cl-濃度增大而增大,表明研究區(qū)地下水化學(xué)組分在地下徑流過程中受到巖石溶解和蒸發(fā)濃縮共同作用。此外,地下水PY3取樣點位于1∶1直線下,地下水LK取樣點位于1∶1直線上方,其Na+與Cl-的變化關(guān)系也揭示了二者與其他地下水取樣點的成因不同。

    地下水中HCO-3主要來源于碳酸鹽礦物溶解[39-40]。γ(Ca2++Mg2+)/γ(HCO-3)的值可以反映地下水碳酸鹽溶解的特征[41]。根據(jù)圖7(a),黃河水γ(Ca2++Mg2+)與γ(HCO-3)關(guān)系點據(jù)聚集于1∶1線上方,說明Ca2++Mg2+不能被HCO-3完全平衡,河水中主要離子來源除了碳酸鹽的溶解外可能還有其他來源;大部分地下水取樣點γ(Ca2++Mg2+)與γ(HCO-3)關(guān)系點據(jù)分布在1∶1線兩側(cè),說明地下水Ca2++Mg2+與HCO-3具有較好的正相關(guān)關(guān)系,地下水中Ca2+、Mg2+、HCO-3來源主要為碳酸鹽巖(方解石和白云石)的溶解。由圖7(b)可知,研究區(qū)黃河水點據(jù)位于硫酸參與碳酸鹽溶解區(qū)域,說明SO2-4參與了Ca2++Mg2+不能被HCO-3平衡部分的溶解,地下水中碳酸廣泛參與碳酸鹽礦物的溶解,個別取樣點碳酸和硫酸共同參與碳酸鹽礦物的溶解(如PY3)。此外,地下水LK取樣點存在石膏溶解。

    (2)人類活動的影響。水體受人類活動排放污染物影響的特征是河水富含Cl-、SO2-4和NO-3[42],研究區(qū)河水的Cl-、SO2-4和NO-3質(zhì)量濃度大于地下水的,說明研究區(qū)黃河水受人類活動的影響大于地下水。黃河從上游至下游沿程水樣水化學(xué)組成中陰離子質(zhì)量濃度HCO-3大于SO2-4的[19,23],本研究2017年7月黃河水樣均出現(xiàn)陰離子質(zhì)量濃度SO2-4大于HCO-3的結(jié)果,河水中SO2-4除來源于自然界蒸發(fā)鹽巖的溶解外,還可能來源于生活污水[32]、工業(yè)活動和大氣沉降等[43]。研究區(qū)上游小浪底附近旅游業(yè)發(fā)達,由此推測夏季旅游旺季,小浪底水庫受排放污染物影響導(dǎo)致庫區(qū)水體SO2-4質(zhì)量濃度增大,同時小浪底水庫汛期排水造成下游出現(xiàn)SO2-4質(zhì)量濃度大于HCO-3的。黃河水NO-3質(zhì)量濃度遠大于地下水說明黃河水相對于地下水受到了更多人類活動的影響。

    地下水個別取樣點在3期水樣中均出現(xiàn)陰離子Cl-質(zhì)量濃度大于SO2-4的情況,而附近其他取樣點SO2-4質(zhì)量濃度大于Cl-的,可能是取樣點附近農(nóng)田施用含Cl-農(nóng)藥或化肥污染地下水所致。

    5 結(jié) 論

    (1)研究區(qū)黃河水pH值為7.40~8.10,TDS為520.55~719.87 mg/L,pH值與TDS均表現(xiàn)出明顯的季節(jié)性差異,這可能與不同季節(jié)各種來源的水對黃河水的補給貢獻不同有關(guān)。地下水pH值總體上低于河水的,大體上與同期黃河水表現(xiàn)出一致的季節(jié)性差異,表明研究區(qū)黃河水對地下水的側(cè)滲補給水循環(huán)交替很快。

    (2)研究區(qū)沿著黃河水側(cè)滲補給地下水方向陽離子Ca2+和Mg2+比例增大、Na+比例減小,陰離子HCO-3比例增大、Cl-和SO42-比例減小,水化學(xué)類型由SO4·Cl·HCO3-Na·Ca·Mg、SO4·HCO3·Cl-Na·Ca·Mg、HCO3·SO4·Cl-Na·Ca·Mg型演變?yōu)镠CO3·Ca-Mg·Na、HCO3-Mg·Ca·Na型。

    (3)研究區(qū)黃河水和地下水均受巖石風(fēng)化和蒸發(fā)濃縮作用影響。地下水中Na+來源除NaCl的溶解外還有離子交換作用,Ca2+、Mg2+主要來源于碳酸鹽巖(方解石和白云石)溶解,個別取樣點存在石膏溶解。地下水中碳酸廣泛參與碳酸鹽礦物的溶解,個別取樣點碳酸和硫酸共同參與碳酸鹽礦物的溶解。

    (4)研究區(qū)黃河水Cl-、SO2-4和NO-3質(zhì)量濃度較高,可能與上游景區(qū)排污有關(guān)。地下水個別取樣點Cl-質(zhì)量濃度較高,可能是取樣點附近施用含Cl-農(nóng)藥或化肥污染地下水所致。

    (5)水化學(xué)成分分析結(jié)果可以作為判斷黃河水與地下水間聯(lián)系是否密切進而劃分黃河影響帶的重要依據(jù)。

    參考文獻:

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

    [2] 劉白薇,董少剛,唐仲華,等.內(nèi)蒙古土默川平原地下水化學(xué)季節(jié)性變化特征[J].工程勘察,2020,48(5):34-40.

    [3] CRONIN A A, BARTH J A C, ELLIOT T, et al. Recharge Velocity and Geochemical Evolution for the Permo-Triassic Sherwood Sandstone, Northern Ireland[J]. Journal of Hydrology, 2005, 315(1-4): 308-324.

    [4] HAMZAOUI-AZAZA F, KETATA M, BOUHLILA R, et al. Hydrogeochemical Characteristics and Assessment of Drinking Water Quality in Zeuss-Koutine Aquifer, Southeastern Tunisia[J]. Environmental Monitoring & Assessment, 2011(1-4): 283-298.

    [5] 王雨山,郭媛,周殷竹,等.基于水化學(xué)和同位素評價馬蓮河下游地下水補給河水的時空變化[J].干旱區(qū)地理,2020,43(2):290-298.

    [6] ADAMS S, TITUS R, PIETERSEN K, et al. Hydrochemical Characteristics of Aquifers Near Sutherland in the Western Karoo, South Africa[J]. Journal of Hydrology, 2001, 241(1-2): 91-103.

    [7] 張巖,董維紅,李滿洲,等.河南平原淺層地下水總?cè)芙夤腆w和水化學(xué)類型的分布特征[J].水文,2011,31(2):79-83.

    [8] 安樂生,趙全升,葉思源,等.黃河三角洲淺層地下水化學(xué)特征及形成作用[J].環(huán)境科學(xué),2012,33(2):370-378.

    [9] HUANG G X, SUN J C, ZHANG Y, et al. Impact of Anthropogenic and Natural Processes on the Evolution of Groundwater Chemistry in a Rapidly Urbanized Coastal Area, South China[J]. Science of the Total Environment, 2013, 463-464(1): 209-221.

    [10] 張學(xué)靜,王平,王田野,等.輸水條件下額濟納綠洲淺層地下水水化學(xué)特征與水位埋深關(guān)系[J].南水北調(diào)與水利科技,2019,17(6):86-94.

    [11] 單曉杰,何江濤,張小文,等.基于人為活動影響識別的區(qū)域地下水水質(zhì)演化預(yù)測:以石家莊地區(qū)為例[J].現(xiàn)代地質(zhì),2020,34(1):189-198.

    [12] 張百祖,宋天琪.疏勒河流域地下水化學(xué)特征及控制因素研究[J].人民黃河,2020,42(6):68-72,78.

    [13] 王福剛.同位素技術(shù)在黃河下游懸河段(河南段)水循環(huán)特征研究中的應(yīng)用[D].長春:吉林大學(xué),2006:1.

    [14] 趙云章,邵景力,閻振鵬,等.黃河下游影響帶地下水資源評價及可持續(xù)開發(fā)利用[M].北京:中國大地出版社,2002:76.

    [15] 孫躍.河南省沿黃地下水?dāng)?shù)值模擬及地下水資源評價[D].北京:中國地質(zhì)大學(xué),2011:14.

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

    [17] 趙云章,邵景力,閆震鵬,等.黃河水對兩側(cè)地下水補給范圍的初步研究[J].人民黃河,2003,25(1):3-5.

    [18] 趙云章,閆震鵬,焦紅軍,等.河南省沿黃城市后備地下水水源地勘查開發(fā)與保護[M].北京:地質(zhì)出版社,2013:51.

    [19] 蘇小四,吳曉芳,林學(xué)鈺,等.黃河水主要化學(xué)組分與δ13C的沿程變化特征[J].人民黃河,2006,28(5):29-31,80.

    [20] 陳靜生,王飛越,何大偉.黃河水質(zhì)地球化學(xué)[J].地學(xué)前緣,2006,13(1):58-73.

    [21] ZHANG Q Q, JIN Z D, ZHANG F, et al. Seasonal Variation in River Water Chemistry of the Middle Reaches of the Yellow River and Its Controlling Factors[J]. Journal of Geochemical Exploration, 2015, 156(3): 101-113.

    [22] ZHANG Y, LI F D, ZHAO G S, et al. An Attempt to Evaluate the Recharge Source and Extent Using Hydrogeochemistry and Stable Isotopes in North Henan Plain, China[J]. Environmental Monitoring and Assessment, 2014(8): 5185-5197.

    [23] 何姜毅,張東,趙志琦.黃河流域河水水化學(xué)組成的時間和空間變化特征[J].生態(tài)學(xué)雜志,2017,36(5):1390-1401.

    [24] 黃奇波,覃小群,劉朋雨,等.柳林泉域巖溶地下水區(qū)域演化規(guī)律及控制因素[J].環(huán)境科學(xué),2019,40(5):2132-2142.

    [25] 蘇小四,林學(xué)鈺,廖資生,等.黃河水δ18O、δD和3H的沿程變化特征及其影響因素研究[J].地球化學(xué),2003,32(4):349-357.

    [26] 曹劍峰,冶雪艷,王福剛,等.黃河懸河段地下水系統(tǒng)分析與模擬[M].鄭州:黃河水利出版社,2006:69.

    [27] 趙云章,朱中道,王繼華.河南省地下水資源與環(huán)境[M].北京:中國大地出版社,2004:176.

    [28] BELKHIRI L, MOUNI L, BOUDOUKHA A. Geochemical Evolution of Groundwater in an Alluvial Aquifer: Case of El Eulma Aquifer, East Algeria[J]. Journal of African Earth Ences, 2012, 66-67: 46-55.

    [29] STADLER S, SLTENFU J, HOLLNDERD H M, et al. Isotopic and Geochemical Indicators for Groundwater Flow and Multi-Component Mixing Near Disturbed Salt Anticlines[J]. Chemical Geology, 2012, 294-295: 226-242.

    [30] 柳鳳霞,史紫薇,錢會,等.銀川地區(qū)地下水水化學(xué)特征演化規(guī)律及水質(zhì)評價[J].環(huán)境化學(xué),2019,38(9):2055-2066.

    [31] 王亞平,王嵐,許春雪,等.長江水系水文地球化學(xué)特征及主要離子的化學(xué)成因[J].地質(zhì)通報,2010,29(增刊1):446-456.

    [32] 張旺,王殿武,雷坤,等.黃河中下游豐水期水化學(xué)特征及影響因素[J].水土保持研究,2020,27(1):380-386,393.

    [33] 沈照理.水文地球化學(xué)基礎(chǔ)[M].北京:地質(zhì)出版社,1999:40.

    [34] HAN D M, LIANG X, JIN M G. Hydrogeochemical Indicators of Groundwater Flow Systems in the Yangwu River Alluvial Fan, Xinzhou Basin, Shanxi, China[J]. Environmental Management, 2009, 44(2): 243-255.

    [35] KUMAR S K, RAMMOHAN V, SAHAYAM J D, et al. Assessment of Groundwater Quality and Hydrogeochemistry of Manimuktha River Basin, Tamil Nadu, India[J]. Environmental Monitoring & Assessment, 2009, 159(1-4): 341-351.

    [36] 劉兵,王賀,姜永海,等.基于水化學(xué)和氫氧同位素的東宮河流域不同水體轉(zhuǎn)化關(guān)系研究[J].環(huán)境科學(xué)研究,2020,33(9):1979-1990.

    [37] CHIDAMBARAM S, ANANDHAN P, PRASANNA M V, et al. Major Ion Chemistry and Identification of Hydrogeochemical Processes Controlling Groundwater in and Around Neyveli Lignite Mines, Tamil Nadu, South India[J]. Arabian Journal of Geosciences, 2013, 6(9): 3451-3467.

    [38] XIAO J, JIN Z D, WANG J, et al. Hydrochemical Characteristics, Controlling Factors and Solute Sources of Groundwater Within the Tarim River Basin in the Extreme Arid Region, NW Tibetan Plateau[J]. Quaternary International, 2015, 380-381: 237-246.

    [39] MONTETY V D, RADAKOVITCH O, VALLET-COULOMB C, et al. Origin of Groundwater Salinity and Hydrogeochemical Processes in a Confined Coastal Aquifer Case of the Rhone delta (Southern France)[J]. Applied Geochemistry, 2008, 23(8): 2337-2349.

    [40] 欒風(fēng)嬌,周金龍,賈瑞亮,等.新疆巴里坤-伊吾盆地地下水水化學(xué)特征及成因[J].環(huán)境化學(xué),2017,36(2):380-389.

    [41] 劉江濤,蔡五田,曹月婷,等.沁河沖洪積扇地下水水化學(xué)特征及成因分析[J].環(huán)境科學(xué),2018,39(12):5428-5439.

    [42] 胡春華,周文斌,夏思奇.鄱陽湖流域水化學(xué)主離子特征及其來源分析[J].環(huán)境化學(xué),2011,30(9):1620-1626.

    [43] 劉叢強,蔣穎魁,陶發(fā)祥,等.西南喀斯特流域碳酸鹽巖的硫酸侵蝕與碳循環(huán)[J].地球化學(xué),2008,37(4):404-414.

    【責(zé)任編輯 呂艷梅】

    猜你喜歡
    黃河水水化學(xué)黃河
    幸福湖泉域地下水水化學(xué)特征及其演化機制
    人民黃河(2025年1期)2025-01-16 00:00:00
    秦皇島北部地表水和地下水水化學(xué)特征及水質(zhì)評價
    多彩黃河
    金橋(2020年11期)2020-12-14 07:52:46
    黃河寧,天下平
    金橋(2020年11期)2020-12-14 07:52:42
    黃河水是膠體嗎
    『黃河』
    黃河知道我愛誰
    兩滴黃河水
    文學(xué)港(2019年5期)2019-05-24 14:19:42
    衡水市新近系館陶組地?zé)崴瘜W(xué)特征與分布規(guī)律研究
    筠連煤層氣井產(chǎn)出水化學(xué)特征及意義
    中國煤層氣(2015年5期)2015-08-22 03:26:32
    亚洲内射少妇av| 久久久久久九九精品二区国产| ponron亚洲| 97超视频在线观看视频| 18禁在线无遮挡免费观看视频| 免费黄频网站在线观看国产| 男女啪啪激烈高潮av片| 高清视频免费观看一区二区 | 午夜福利成人在线免费观看| av国产久精品久网站免费入址| 国产精品.久久久| 欧美性感艳星| 一级毛片 在线播放| 高清在线视频一区二区三区| 亚洲国产精品国产精品| 亚洲欧美一区二区三区黑人 | 街头女战士在线观看网站| 国产淫片久久久久久久久| 久久久精品欧美日韩精品| 日韩精品青青久久久久久| 日韩国内少妇激情av| 伊人久久国产一区二区| 亚洲乱码一区二区免费版| 国产午夜精品论理片| 精品久久久久久成人av| 亚洲四区av| 日本爱情动作片www.在线观看| 亚洲精华国产精华液的使用体验| 亚洲激情五月婷婷啪啪| 婷婷六月久久综合丁香| 国产视频内射| 婷婷色综合大香蕉| 亚洲精品乱码久久久v下载方式| 菩萨蛮人人尽说江南好唐韦庄| 麻豆国产97在线/欧美| 免费看光身美女| av国产久精品久网站免费入址| 久久久久网色| 欧美 日韩 精品 国产| 亚洲av不卡在线观看| 黄片无遮挡物在线观看| 国产精品国产三级国产av玫瑰| 内地一区二区视频在线| 精品久久久久久久久亚洲| 18禁动态无遮挡网站| 日韩av在线免费看完整版不卡| 亚洲成人一二三区av| eeuss影院久久| 成人亚洲精品一区在线观看 | 亚洲图色成人| 五月天丁香电影| 亚洲成色77777| 欧美成人a在线观看| 精品一区二区免费观看| 久久久久久久久久人人人人人人| 蜜桃亚洲精品一区二区三区| 亚洲精品第二区| 22中文网久久字幕| 免费av不卡在线播放| 欧美潮喷喷水| 免费黄频网站在线观看国产| 久久久久久久午夜电影| 97人妻精品一区二区三区麻豆| 白带黄色成豆腐渣| 天堂av国产一区二区熟女人妻| 精品99又大又爽又粗少妇毛片| 最近最新中文字幕免费大全7| 久久久精品免费免费高清| 校园人妻丝袜中文字幕| 久久精品夜夜夜夜夜久久蜜豆| 男女边摸边吃奶| 中文字幕免费在线视频6| 久久久久九九精品影院| 国产综合精华液| 亚洲电影在线观看av| 免费av观看视频| 亚洲av在线观看美女高潮| 色哟哟·www| 又大又黄又爽视频免费| 99久久人妻综合| 久久久国产一区二区| 日本与韩国留学比较| 欧美一级a爱片免费观看看| 婷婷色av中文字幕| 爱豆传媒免费全集在线观看| 可以在线观看毛片的网站| 久久久久久久大尺度免费视频| 男人舔女人下体高潮全视频| 国内少妇人妻偷人精品xxx网站| 欧美人与善性xxx| 久久久a久久爽久久v久久| 搡老乐熟女国产| 国产毛片a区久久久久| 久久久精品免费免费高清| 亚洲欧美日韩卡通动漫| 最近视频中文字幕2019在线8| 日日撸夜夜添| 亚洲欧美日韩无卡精品| 成人毛片60女人毛片免费| 蜜桃亚洲精品一区二区三区| 国产精品一区www在线观看| 亚洲精品成人久久久久久| 日本-黄色视频高清免费观看| 久久亚洲国产成人精品v| 舔av片在线| 有码 亚洲区| 99久久精品热视频| av免费观看日本| 天堂av国产一区二区熟女人妻| 丝瓜视频免费看黄片| 日本猛色少妇xxxxx猛交久久| 亚洲经典国产精华液单| 搞女人的毛片| 国产精品人妻久久久久久| 成人欧美大片| 中文字幕免费在线视频6| 99热全是精品| 久久久久久久午夜电影| 黑人高潮一二区| 看免费成人av毛片| 欧美人与善性xxx| 日韩欧美一区视频在线观看 | 久久午夜福利片| av在线天堂中文字幕| 一级片'在线观看视频| 全区人妻精品视频| 成人亚洲欧美一区二区av| 久久精品国产亚洲网站| 日韩 亚洲 欧美在线| 午夜精品在线福利| 简卡轻食公司| 日韩亚洲欧美综合| 男插女下体视频免费在线播放| 性插视频无遮挡在线免费观看| 色哟哟·www| 男女边吃奶边做爰视频| 人妻夜夜爽99麻豆av| 一级毛片aaaaaa免费看小| 久久久久久国产a免费观看| 国产成人一区二区在线| 男人爽女人下面视频在线观看| 欧美成人午夜免费资源| 91狼人影院| 91精品伊人久久大香线蕉| 成年av动漫网址| 成人亚洲欧美一区二区av| 免费播放大片免费观看视频在线观看| 精品久久久精品久久久| 亚洲av不卡在线观看| 亚洲最大成人av| 人人妻人人澡欧美一区二区| 亚洲精品日韩在线中文字幕| 欧美成人午夜免费资源| 免费观看性生交大片5| 永久网站在线| 国产精品不卡视频一区二区| 人人妻人人看人人澡| 国产欧美日韩精品一区二区| 亚洲丝袜综合中文字幕| 欧美成人一区二区免费高清观看| 国产成人a∨麻豆精品| 午夜精品在线福利| 午夜老司机福利剧场| 亚洲高清免费不卡视频| 欧美性猛交╳xxx乱大交人| 亚洲熟妇中文字幕五十中出| 国产欧美另类精品又又久久亚洲欧美| 日本免费在线观看一区| 午夜久久久久精精品| 久久久久久久久久久丰满| 国产在视频线精品| 亚洲婷婷狠狠爱综合网| 91在线精品国自产拍蜜月| 亚洲人与动物交配视频| 久久精品夜色国产| 男女边吃奶边做爰视频| 夫妻午夜视频| 色网站视频免费| 亚洲人与动物交配视频| 最近的中文字幕免费完整| 永久网站在线| 99久久精品一区二区三区| 中文乱码字字幕精品一区二区三区 | 国产伦精品一区二区三区四那| 亚洲av免费高清在线观看| 少妇的逼好多水| 插逼视频在线观看| 老女人水多毛片| 亚洲综合精品二区| 日本黄色片子视频| 干丝袜人妻中文字幕| 精品人妻视频免费看| 日韩一区二区视频免费看| 水蜜桃什么品种好| 九九久久精品国产亚洲av麻豆| 亚洲熟女精品中文字幕| 国产精品人妻久久久久久| www.av在线官网国产| 午夜精品国产一区二区电影 | 国产精品av视频在线免费观看| 街头女战士在线观看网站| 午夜免费男女啪啪视频观看| 国产片特级美女逼逼视频| 欧美日韩在线观看h| 欧美一级a爱片免费观看看| 美女被艹到高潮喷水动态| 直男gayav资源| 精品久久久久久成人av| 欧美高清性xxxxhd video| 久久久久久久亚洲中文字幕| 啦啦啦啦在线视频资源| 中文字幕免费在线视频6| 老女人水多毛片| 国产在线男女| 欧美激情久久久久久爽电影| 插逼视频在线观看| 国产在线一区二区三区精| 午夜精品在线福利| 天堂俺去俺来也www色官网 | 最近的中文字幕免费完整| 日日撸夜夜添| 亚洲一级一片aⅴ在线观看| 午夜亚洲福利在线播放| 国产亚洲av嫩草精品影院| 最近2019中文字幕mv第一页| 欧美性感艳星| 免费观看的影片在线观看| 又爽又黄无遮挡网站| 欧美性感艳星| 春色校园在线视频观看| 一级av片app| 18+在线观看网站| 亚洲在久久综合| 国产在视频线精品| 久久久精品欧美日韩精品| 好男人视频免费观看在线| 亚洲av福利一区| 街头女战士在线观看网站| 日本-黄色视频高清免费观看| 国产亚洲5aaaaa淫片| 日本黄大片高清| 五月伊人婷婷丁香| 97超视频在线观看视频| 99九九线精品视频在线观看视频| av卡一久久| 国产亚洲最大av| 亚洲欧美成人综合另类久久久| 22中文网久久字幕| 看非洲黑人一级黄片| 国产午夜福利久久久久久| 久久99精品国语久久久| 三级男女做爰猛烈吃奶摸视频| 男女边摸边吃奶| 免费av不卡在线播放| 国产亚洲5aaaaa淫片| 欧美xxⅹ黑人| 十八禁网站网址无遮挡 | 老师上课跳d突然被开到最大视频| 精品久久久久久久人妻蜜臀av| 一个人看的www免费观看视频| 欧美不卡视频在线免费观看| 九九在线视频观看精品| 大话2 男鬼变身卡| 国产精品.久久久| 亚洲在线自拍视频| 午夜免费观看性视频| 亚洲一区高清亚洲精品| 精品99又大又爽又粗少妇毛片| 国产黄色视频一区二区在线观看| 日本爱情动作片www.在线观看| 一边亲一边摸免费视频| 国产免费又黄又爽又色| 啦啦啦啦在线视频资源| 免费看av在线观看网站| 毛片女人毛片| 国产精品国产三级国产专区5o| 久久精品综合一区二区三区| 国产精品av视频在线免费观看| 两个人的视频大全免费| 亚洲精品一二三| 高清欧美精品videossex| ponron亚洲| 国产熟女欧美一区二区| h日本视频在线播放| 久久久久免费精品人妻一区二区| 别揉我奶头 嗯啊视频| 卡戴珊不雅视频在线播放| 亚洲国产精品成人综合色| 日韩制服骚丝袜av| 婷婷色麻豆天堂久久| 免费观看在线日韩| 亚洲欧美精品专区久久| 欧美区成人在线视频| 国产精品一区二区在线观看99 | 国产午夜精品久久久久久一区二区三区| 伦理电影大哥的女人| 亚洲精品亚洲一区二区| 黑人高潮一二区| 黄片无遮挡物在线观看| 三级国产精品欧美在线观看| 国产真实伦视频高清在线观看| 99久久精品一区二区三区| av天堂中文字幕网| 哪个播放器可以免费观看大片| 男女边摸边吃奶| 国产精品人妻久久久久久| 99视频精品全部免费 在线| 亚洲18禁久久av| 麻豆av噜噜一区二区三区| 色网站视频免费| 亚洲精品,欧美精品| 又粗又硬又长又爽又黄的视频| 国模一区二区三区四区视频| 波多野结衣巨乳人妻| 日韩不卡一区二区三区视频在线| 中文精品一卡2卡3卡4更新| 六月丁香七月| 亚洲精品中文字幕在线视频 | 日韩国内少妇激情av| av又黄又爽大尺度在线免费看| 天天躁日日操中文字幕| 亚洲欧美成人精品一区二区| 又黄又爽又刺激的免费视频.| 国产永久视频网站| 国产午夜精品久久久久久一区二区三区| 国产男人的电影天堂91| 亚洲高清免费不卡视频| 免费观看的影片在线观看| 别揉我奶头 嗯啊视频| 欧美激情久久久久久爽电影| 日本与韩国留学比较| 国产视频首页在线观看| 少妇的逼好多水| 精品一区二区三区视频在线| 国产成人a区在线观看| 欧美日韩综合久久久久久| 99视频精品全部免费 在线| 成人国产麻豆网| 在线观看美女被高潮喷水网站| 国产精品国产三级国产专区5o| 99久久精品国产国产毛片| 女的被弄到高潮叫床怎么办| 国产高清三级在线| 777米奇影视久久| 国产成人免费观看mmmm| 免费看a级黄色片| 国产黄色小视频在线观看| 久久精品夜色国产| 99热全是精品| 欧美成人精品欧美一级黄| 久久精品国产亚洲av天美| 特级一级黄色大片| 只有这里有精品99| 只有这里有精品99| 国产午夜精品久久久久久一区二区三区| 97人妻精品一区二区三区麻豆| 免费观看无遮挡的男女| 午夜亚洲福利在线播放| 高清欧美精品videossex| 插阴视频在线观看视频| 卡戴珊不雅视频在线播放| 嫩草影院精品99| 99热这里只有是精品在线观看| 欧美97在线视频| 韩国高清视频一区二区三区| 久久99热6这里只有精品| 亚洲精品一二三| 日本熟妇午夜| 久久精品久久久久久噜噜老黄| 亚洲国产av新网站| 国产女主播在线喷水免费视频网站 | 看免费成人av毛片| 精华霜和精华液先用哪个| 成人漫画全彩无遮挡| 亚洲av成人av| 国产成人aa在线观看| 大又大粗又爽又黄少妇毛片口| 国产永久视频网站| 久久99精品国语久久久| 国产精品av视频在线免费观看| 禁无遮挡网站| 日韩欧美一区视频在线观看 | 欧美+日韩+精品| 高清日韩中文字幕在线| 免费观看的影片在线观看| 欧美一级a爱片免费观看看| 亚洲av中文av极速乱| av免费观看日本| 国产国拍精品亚洲av在线观看| 一级av片app| 搞女人的毛片| 日韩一本色道免费dvd| 在线免费观看的www视频| a级毛片免费高清观看在线播放| 日韩成人av中文字幕在线观看| 亚洲精品久久午夜乱码| 国产免费一级a男人的天堂| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 老司机影院毛片| av专区在线播放| 在线观看av片永久免费下载| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品成人av观看孕妇| 久久国内精品自在自线图片| av线在线观看网站| 国产精品美女特级片免费视频播放器| 亚洲欧美成人精品一区二区| 国产色爽女视频免费观看| 国产乱人视频| 69人妻影院| 欧美日韩亚洲高清精品| 久久精品久久精品一区二区三区| 免费不卡的大黄色大毛片视频在线观看 | 又爽又黄a免费视频| 国产精品综合久久久久久久免费| 嫩草影院精品99| 国产久久久一区二区三区| 22中文网久久字幕| 偷拍熟女少妇极品色| 如何舔出高潮| 精品熟女少妇av免费看| av在线亚洲专区| 久久精品国产自在天天线| 一级毛片我不卡| 国产一区二区亚洲精品在线观看| 一级黄片播放器| 少妇高潮的动态图| 免费av毛片视频| 最近的中文字幕免费完整| 女人十人毛片免费观看3o分钟| 亚洲天堂国产精品一区在线| 午夜福利网站1000一区二区三区| 性色avwww在线观看| 国产v大片淫在线免费观看| 精品国产一区二区三区久久久樱花 | 黄色配什么色好看| 亚洲精品成人久久久久久| 美女内射精品一级片tv| 欧美精品国产亚洲| 在线免费十八禁| 在线观看人妻少妇| av免费在线看不卡| 舔av片在线| 欧美激情国产日韩精品一区| 久久久久精品久久久久真实原创| 91久久精品电影网| 噜噜噜噜噜久久久久久91| 色尼玛亚洲综合影院| 91午夜精品亚洲一区二区三区| 国内精品一区二区在线观看| 亚洲一级一片aⅴ在线观看| 少妇熟女欧美另类| 成人无遮挡网站| 久久这里有精品视频免费| 亚洲丝袜综合中文字幕| 国产精品麻豆人妻色哟哟久久 | 大片免费播放器 马上看| 搡老妇女老女人老熟妇| 乱码一卡2卡4卡精品| 国产精品福利在线免费观看| 国产亚洲精品久久久com| 亚洲乱码一区二区免费版| av在线天堂中文字幕| 国产综合懂色| 深夜a级毛片| 大香蕉97超碰在线| 午夜亚洲福利在线播放| 久久久久久久亚洲中文字幕| 日本色播在线视频| 两个人的视频大全免费| 日韩强制内射视频| 老女人水多毛片| 欧美日韩在线观看h| 亚洲精华国产精华液的使用体验| 亚洲人成网站高清观看| 日日啪夜夜撸| 日日干狠狠操夜夜爽| 国产综合懂色| 哪个播放器可以免费观看大片| 国产一区二区三区综合在线观看 | 国产91av在线免费观看| 欧美成人a在线观看| 真实男女啪啪啪动态图| 精品人妻视频免费看| av卡一久久| 国产精品精品国产色婷婷| 中文字幕人妻熟人妻熟丝袜美| 2022亚洲国产成人精品| 欧美xxⅹ黑人| 大香蕉久久网| 国内精品一区二区在线观看| 黑人高潮一二区| 国产一区二区三区av在线| 亚洲人成网站在线观看播放| 老女人水多毛片| 国产永久视频网站| 伊人久久国产一区二区| 国产精品久久久久久精品电影小说 | 国产成人精品福利久久| 国产白丝娇喘喷水9色精品| 日日干狠狠操夜夜爽| 黄色配什么色好看| 丝袜美腿在线中文| 三级男女做爰猛烈吃奶摸视频| 中国美白少妇内射xxxbb| 日韩不卡一区二区三区视频在线| 国产亚洲最大av| 简卡轻食公司| 日韩av不卡免费在线播放| 精品午夜福利在线看| 亚洲av成人精品一区久久| 亚洲av福利一区| 男插女下体视频免费在线播放| 国产精品1区2区在线观看.| 毛片一级片免费看久久久久| 久久亚洲国产成人精品v| av女优亚洲男人天堂| 精品酒店卫生间| 免费看不卡的av| 波多野结衣巨乳人妻| 亚洲欧洲国产日韩| 日韩精品有码人妻一区| 亚洲精品日本国产第一区| 少妇猛男粗大的猛烈进出视频 | 非洲黑人性xxxx精品又粗又长| 国产成人精品久久久久久| av.在线天堂| 18禁在线播放成人免费| 成人美女网站在线观看视频| 美女大奶头视频| 能在线免费观看的黄片| 大香蕉久久网| 久久草成人影院| 纵有疾风起免费观看全集完整版 | 免费黄色在线免费观看| 欧美日韩国产mv在线观看视频 | 最近中文字幕高清免费大全6| 91在线精品国自产拍蜜月| 国产午夜精品一二区理论片| 国产精品人妻久久久久久| av一本久久久久| 亚洲国产精品成人综合色| 中文资源天堂在线| 精品久久久久久久久av| 免费黄频网站在线观看国产| 精品久久久久久久久久久久久| 国产伦精品一区二区三区视频9| 99久国产av精品| 欧美日韩在线观看h| 99视频精品全部免费 在线| 久久99热这里只频精品6学生| 日韩av在线大香蕉| 看免费成人av毛片| 天美传媒精品一区二区| 91狼人影院| 狂野欧美激情性xxxx在线观看| 午夜免费观看性视频| 免费观看精品视频网站| 黑人高潮一二区| 国产成人aa在线观看| 亚洲国产成人一精品久久久| 啦啦啦啦在线视频资源| freevideosex欧美| 少妇丰满av| 亚洲久久久久久中文字幕| 看非洲黑人一级黄片| 久久久成人免费电影| av黄色大香蕉| 亚洲一级一片aⅴ在线观看| 亚洲国产成人一精品久久久| 亚洲在线观看片| 国产精品一及| 精品人妻一区二区三区麻豆| 成人综合一区亚洲| 亚洲一级一片aⅴ在线观看| 日韩欧美国产在线观看| 亚洲综合色惰| 老师上课跳d突然被开到最大视频| 久久久久久久久大av| 国产在视频线在精品| 免费观看a级毛片全部| 国产三级在线视频| 久久精品久久久久久久性| 日日啪夜夜爽| 亚洲三级黄色毛片| 在线免费观看的www视频| 国产亚洲av片在线观看秒播厂 | 在线免费观看的www视频| 丝袜喷水一区| 天堂av国产一区二区熟女人妻| av又黄又爽大尺度在线免费看| 欧美一区二区亚洲| 亚洲国产日韩欧美精品在线观看| 99九九线精品视频在线观看视频| 国产精品人妻久久久久久| 亚洲精品久久午夜乱码| 最近的中文字幕免费完整| 晚上一个人看的免费电影| 又黄又爽又刺激的免费视频.| 久久久久久久久久黄片| 午夜精品一区二区三区免费看| 免费电影在线观看免费观看| 国产 一区 欧美 日韩| 久久久久久久久久黄片| 国产淫语在线视频| 国产黄片美女视频| 久久午夜福利片| 高清午夜精品一区二区三区| 欧美xxⅹ黑人| 国产亚洲5aaaaa淫片| 国产精品人妻久久久影院| 成人一区二区视频在线观看| 色哟哟·www| 国产精品人妻久久久久久| 日韩成人av中文字幕在线观看|