• <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新网站| 丰满迷人的少妇在线观看| a 毛片基地| 亚洲一级一片aⅴ在线观看| 亚洲熟女精品中文字幕| 亚洲欧美精品自产自拍| 亚洲久久久国产精品| 日韩制服丝袜自拍偷拍| 女性被躁到高潮视频| 亚洲 欧美一区二区三区| 亚洲欧洲日产国产| 精品一区二区三区四区五区乱码 | 男人舔女人的私密视频| 国产高清三级在线| 97超碰精品成人国产| 大香蕉久久网| 婷婷色麻豆天堂久久| 国产一区二区三区av在线| 国产一区二区三区av在线| 精品久久蜜臀av无| 一区二区三区四区激情视频| av电影中文网址| 性色av一级| 亚洲伊人色综图| 亚洲第一av免费看| 国产精品人妻久久久久久| 宅男免费午夜| av有码第一页| 波多野结衣一区麻豆| 国产成人免费观看mmmm| 亚洲色图综合在线观看| 亚洲五月色婷婷综合| 两个人免费观看高清视频| 色网站视频免费| 精品久久久久久电影网| xxxhd国产人妻xxx| 精品99又大又爽又粗少妇毛片| 午夜福利乱码中文字幕| 亚洲av电影在线进入| 最新中文字幕久久久久| 男女边吃奶边做爰视频| 春色校园在线视频观看| 久久久久久人妻| 精品人妻在线不人妻| videosex国产| 亚洲精品久久久久久婷婷小说| 午夜av观看不卡| 草草在线视频免费看| 亚洲国产欧美在线一区| 国产精品欧美亚洲77777| 国产亚洲午夜精品一区二区久久| 色婷婷av一区二区三区视频| 精品午夜福利在线看| 九色成人免费人妻av| 有码 亚洲区| 亚洲欧洲日产国产| videossex国产| 一二三四在线观看免费中文在 | 亚洲四区av| 最近手机中文字幕大全| 9热在线视频观看99| 午夜福利影视在线免费观看| 亚洲欧美日韩卡通动漫| 欧美精品av麻豆av| 中文字幕人妻熟女乱码| 18在线观看网站| 亚洲欧美精品自产自拍| 精品国产一区二区三区四区第35| 91精品伊人久久大香线蕉| 久久久久人妻精品一区果冻| 国产精品欧美亚洲77777| 少妇的丰满在线观看| 中文字幕人妻熟女乱码| 日韩中字成人| 69精品国产乱码久久久| 免费观看在线日韩| 18在线观看网站| 精品午夜福利在线看| 久久久久网色| 丝袜喷水一区| 91午夜精品亚洲一区二区三区| 十分钟在线观看高清视频www| 咕卡用的链子| 日韩制服丝袜自拍偷拍| 女性生殖器流出的白浆| 一级片'在线观看视频| 在线观看免费高清a一片| 亚洲,一卡二卡三卡| 香蕉精品网在线| 亚洲精品久久成人aⅴ小说| 国产白丝娇喘喷水9色精品| 亚洲国产精品专区欧美| 亚洲成人手机| 日韩一区二区三区影片| 捣出白浆h1v1| 亚洲色图 男人天堂 中文字幕 | 色婷婷久久久亚洲欧美| 免费黄频网站在线观看国产| 国产白丝娇喘喷水9色精品| 国产黄频视频在线观看| 18禁裸乳无遮挡动漫免费视频| 中文字幕av电影在线播放| 晚上一个人看的免费电影| 热re99久久精品国产66热6| 婷婷色麻豆天堂久久| 国产精品偷伦视频观看了| 69精品国产乱码久久久| 免费播放大片免费观看视频在线观看| 黄片播放在线免费| 亚洲欧美清纯卡通| 又黄又粗又硬又大视频| 99热6这里只有精品| 九色亚洲精品在线播放| 黄色 视频免费看| 乱人伦中国视频| 久久女婷五月综合色啪小说| 中文字幕另类日韩欧美亚洲嫩草| 国产精品 国内视频| 国产一区二区在线观看日韩| 啦啦啦在线观看免费高清www| 黑人欧美特级aaaaaa片| 久久精品国产综合久久久 | 欧美成人精品欧美一级黄| 亚洲av在线观看美女高潮| 亚洲第一av免费看| 亚洲国产日韩一区二区| 满18在线观看网站| 亚洲av中文av极速乱| 69精品国产乱码久久久| 国产麻豆69| 日韩大片免费观看网站| 亚洲欧洲日产国产| 中国三级夫妇交换| 多毛熟女@视频| 婷婷色综合www| www日本在线高清视频| 咕卡用的链子| 99九九在线精品视频| 亚洲成人一二三区av| 成人午夜精彩视频在线观看| 精品一区二区三区四区五区乱码 | 国产无遮挡羞羞视频在线观看| 国产精品国产av在线观看| 日本黄大片高清| 亚洲美女视频黄频| 99九九在线精品视频| a级毛色黄片| av有码第一页| 18禁裸乳无遮挡动漫免费视频| 99九九在线精品视频| 久久人人97超碰香蕉20202| 人妻 亚洲 视频| 国产精品一国产av| 亚洲精品av麻豆狂野| 亚洲av综合色区一区| 国产成人午夜福利电影在线观看| 男女免费视频国产| 九色成人免费人妻av| 日本欧美国产在线视频| 热99久久久久精品小说推荐| 汤姆久久久久久久影院中文字幕| 黑人猛操日本美女一级片| a 毛片基地| 日韩制服丝袜自拍偷拍| 亚洲一码二码三码区别大吗| 街头女战士在线观看网站| 少妇被粗大猛烈的视频| 欧美老熟妇乱子伦牲交| 黄色一级大片看看| 久久国产精品大桥未久av| 亚洲精品久久久久久婷婷小说| 七月丁香在线播放| 亚洲精品av麻豆狂野| 最近最新中文字幕免费大全7| 久久久久国产网址| 成人亚洲欧美一区二区av| 精品亚洲成a人片在线观看| 大香蕉久久成人网| 又大又黄又爽视频免费| 欧美日韩一区二区视频在线观看视频在线| 丝袜人妻中文字幕| 波多野结衣一区麻豆| 69精品国产乱码久久久| 免费久久久久久久精品成人欧美视频 | 人人妻人人添人人爽欧美一区卜| 成人漫画全彩无遮挡| 亚洲 欧美一区二区三区| 欧美国产精品va在线观看不卡| av国产久精品久网站免费入址| 日韩av免费高清视频| 激情五月婷婷亚洲| 久久精品久久精品一区二区三区| 五月天丁香电影| 午夜福利视频在线观看免费| 看十八女毛片水多多多| 色婷婷av一区二区三区视频| 久久精品久久久久久久性| 黄色怎么调成土黄色| 国产乱人偷精品视频| 国产不卡av网站在线观看| 日韩熟女老妇一区二区性免费视频| 最黄视频免费看| 毛片一级片免费看久久久久| 亚洲人与动物交配视频| 日本-黄色视频高清免费观看| 日韩中文字幕视频在线看片| 人妻少妇偷人精品九色| 少妇的逼水好多| 国产永久视频网站| 丰满少妇做爰视频| 免费黄频网站在线观看国产| 亚洲精品自拍成人| 久久精品久久精品一区二区三区| 亚洲国产看品久久| 国产一区二区三区av在线| 国产精品女同一区二区软件| 亚洲精品中文字幕在线视频| 国产精品 国内视频| www.熟女人妻精品国产 | 精品亚洲成a人片在线观看| 亚洲av国产av综合av卡| 亚洲av福利一区| 老司机影院成人| 欧美97在线视频| 亚洲第一区二区三区不卡| 日本vs欧美在线观看视频| 亚洲综合精品二区| 欧美激情 高清一区二区三区| 亚洲av日韩在线播放| 丰满乱子伦码专区| 亚洲欧美一区二区三区国产| 成年女人在线观看亚洲视频| 男人操女人黄网站| 国产精品欧美亚洲77777| 国产在线一区二区三区精| 最近中文字幕2019免费版| 夫妻午夜视频| 免费黄网站久久成人精品| 午夜视频国产福利| 国产 一区精品| 9热在线视频观看99| 18禁国产床啪视频网站| 国产欧美日韩综合在线一区二区| 久久久久精品性色| 久久久久国产网址| 久久久国产一区二区| 内地一区二区视频在线| 亚洲熟女精品中文字幕| 国产免费一区二区三区四区乱码| 国产精品人妻久久久久久| 欧美bdsm另类| 国产片内射在线| 午夜福利视频在线观看免费| 亚洲第一区二区三区不卡| 国产一区亚洲一区在线观看| 久久97久久精品| freevideosex欧美| 欧美另类一区| 免费日韩欧美在线观看| 18禁观看日本| 人妻人人澡人人爽人人| 国产乱人偷精品视频| 伦精品一区二区三区| 免费观看在线日韩| 国产欧美日韩综合在线一区二区| 欧美国产精品一级二级三级| 人人妻人人添人人爽欧美一区卜| 如何舔出高潮| h视频一区二区三区| 99re6热这里在线精品视频| 视频在线观看一区二区三区| 久久ye,这里只有精品| 久久国内精品自在自线图片| 天天躁夜夜躁狠狠躁躁| 久久久精品免费免费高清| 最近中文字幕2019免费版| 国产一区二区在线观看日韩| 一本色道久久久久久精品综合| 五月玫瑰六月丁香| 最近的中文字幕免费完整| 欧美日韩一区二区视频在线观看视频在线| 色哟哟·www| 99热这里只有是精品在线观看| 美女中出高潮动态图| 26uuu在线亚洲综合色| videossex国产| 午夜福利乱码中文字幕| 七月丁香在线播放| 亚洲精品456在线播放app| 亚洲精品自拍成人| 国产一区二区在线观看av| a级毛色黄片| 国产精品一区二区在线观看99| 午夜免费鲁丝| 国产精品久久久av美女十八| 蜜臀久久99精品久久宅男| 国产精品久久久久久精品电影小说| 国产伦理片在线播放av一区| 日本vs欧美在线观看视频| 国产黄色视频一区二区在线观看| 深夜精品福利| 爱豆传媒免费全集在线观看| 精品少妇内射三级| 国产精品不卡视频一区二区| 一本—道久久a久久精品蜜桃钙片| 九九爱精品视频在线观看| 天堂中文最新版在线下载| 久久午夜福利片| 久久久久久久亚洲中文字幕| 成人国语在线视频| 丝袜在线中文字幕| 国产精品一国产av| 王馨瑶露胸无遮挡在线观看| 国产精品.久久久| 精品久久蜜臀av无| 国产一级毛片在线| 亚洲天堂av无毛| 熟女人妻精品中文字幕| 人人妻人人澡人人爽人人夜夜| 久久精品久久精品一区二区三区| 国产乱人偷精品视频| 国产精品国产三级专区第一集| 国产片特级美女逼逼视频| 久久久久久久久久成人| 一边亲一边摸免费视频| 国产精品一二三区在线看| 日韩人妻精品一区2区三区| 综合色丁香网| 亚洲精品久久成人aⅴ小说| 欧美日韩视频精品一区| 国产在视频线精品| kizo精华| 久久这里只有精品19| 一本大道久久a久久精品| 中国美白少妇内射xxxbb| 亚洲国产av影院在线观看| 欧美精品亚洲一区二区| 精品国产露脸久久av麻豆| 中文字幕最新亚洲高清| 十八禁高潮呻吟视频| 亚洲精品日韩在线中文字幕| 一级毛片我不卡| 亚洲国产精品999| 中国三级夫妇交换| 热re99久久精品国产66热6| av.在线天堂| 黑人欧美特级aaaaaa片| 免费黄频网站在线观看国产| 精品午夜福利在线看| 国产乱来视频区| 国产免费现黄频在线看| 国产精品欧美亚洲77777| 侵犯人妻中文字幕一二三四区| 赤兔流量卡办理| 日本wwww免费看| 97人妻天天添夜夜摸| 90打野战视频偷拍视频| 2018国产大陆天天弄谢| 内地一区二区视频在线| 两性夫妻黄色片 | 中国三级夫妇交换| 亚洲精品av麻豆狂野| 精品人妻在线不人妻| 免费在线观看黄色视频的| 中文字幕亚洲精品专区| 国产精品国产三级专区第一集| 最近手机中文字幕大全| 免费观看在线日韩| 亚洲欧美清纯卡通| 久久久久国产网址| 丰满迷人的少妇在线观看| 欧美精品高潮呻吟av久久| 国产黄频视频在线观看| 国产精品久久久久久精品电影小说| 高清视频免费观看一区二区| 国国产精品蜜臀av免费| 人妻 亚洲 视频| 各种免费的搞黄视频| 999精品在线视频| 国产成人精品无人区| 美女福利国产在线| 秋霞伦理黄片| av免费在线看不卡| 亚洲av综合色区一区| 永久网站在线| 99九九在线精品视频| 亚洲国产日韩一区二区| videossex国产| 亚洲精品久久成人aⅴ小说| 免费av中文字幕在线| 成年动漫av网址| 18禁裸乳无遮挡动漫免费视频| 在线免费观看不下载黄p国产| 亚洲精品久久久久久婷婷小说| 国语对白做爰xxxⅹ性视频网站| 在线精品无人区一区二区三| 欧美精品av麻豆av| 国产日韩欧美在线精品| 亚洲一码二码三码区别大吗| 秋霞在线观看毛片| 国产成人精品无人区| 久久午夜福利片| 亚洲精品美女久久久久99蜜臀 | 女性被躁到高潮视频| 色94色欧美一区二区| 成人影院久久| av黄色大香蕉| 香蕉国产在线看| 午夜福利影视在线免费观看| 亚洲欧美精品自产自拍| 一区二区三区乱码不卡18| 亚洲成人av在线免费| xxxhd国产人妻xxx| 最近中文字幕2019免费版| av不卡在线播放| 亚洲成色77777| 在现免费观看毛片| 国产探花极品一区二区| 中文字幕最新亚洲高清| 欧美+日韩+精品| 9色porny在线观看| 国产一级毛片在线| 亚洲av成人精品一二三区| 秋霞在线观看毛片| 伊人亚洲综合成人网| 美女内射精品一级片tv| 肉色欧美久久久久久久蜜桃| 亚洲婷婷狠狠爱综合网| 熟女av电影| 国精品久久久久久国模美| 精品少妇黑人巨大在线播放| 91在线精品国自产拍蜜月| 熟女av电影| 国国产精品蜜臀av免费| 巨乳人妻的诱惑在线观看| 免费久久久久久久精品成人欧美视频 | 亚洲美女搞黄在线观看| 国产一区二区三区av在线| 人妻少妇偷人精品九色| 超色免费av| 国产黄色视频一区二区在线观看| 边亲边吃奶的免费视频| 麻豆精品久久久久久蜜桃| 国产精品 国内视频| 一区二区三区乱码不卡18| 午夜福利网站1000一区二区三区| 日日撸夜夜添| 日韩一区二区视频免费看| 在线观看国产h片| 久久久久久久久久成人| 校园人妻丝袜中文字幕| 在线精品无人区一区二区三| 亚洲国产最新在线播放| 亚洲av电影在线进入| 欧美国产精品va在线观看不卡| 美女脱内裤让男人舔精品视频| 啦啦啦视频在线资源免费观看| 亚洲中文av在线| 日日摸夜夜添夜夜爱| 久久久久精品人妻al黑| 成人无遮挡网站| 不卡视频在线观看欧美| 国产成人91sexporn| 性色av一级| 天天躁夜夜躁狠狠久久av| 美女内射精品一级片tv| 美女大奶头黄色视频| 欧美成人午夜精品| 免费人妻精品一区二区三区视频| 在线观看www视频免费| 最黄视频免费看| 18禁动态无遮挡网站| 亚洲美女黄色视频免费看| 欧美成人精品欧美一级黄| 五月伊人婷婷丁香| 97在线人人人人妻| 女人精品久久久久毛片| h视频一区二区三区| 国产成人精品婷婷| 精品亚洲成国产av| 久久午夜综合久久蜜桃| av黄色大香蕉| 狂野欧美激情性xxxx在线观看| 赤兔流量卡办理| 亚洲,一卡二卡三卡| 伦精品一区二区三区| 国产精品人妻久久久影院| 欧美激情 高清一区二区三区| 亚洲精品国产色婷婷电影| 丁香六月天网| 国产欧美另类精品又又久久亚洲欧美| 在线观看一区二区三区激情| 亚洲欧洲精品一区二区精品久久久 | 高清黄色对白视频在线免费看| 欧美 亚洲 国产 日韩一| 日本免费在线观看一区| 国产黄色视频一区二区在线观看| 丝袜人妻中文字幕| 少妇 在线观看| 人人妻人人澡人人看| 熟妇人妻不卡中文字幕| 观看av在线不卡| 国产色爽女视频免费观看| av免费在线看不卡| 熟妇人妻不卡中文字幕| 多毛熟女@视频| 精品一品国产午夜福利视频| 天美传媒精品一区二区| 精品国产露脸久久av麻豆| 大片免费播放器 马上看| 天天操日日干夜夜撸| 亚洲国产色片| 男人操女人黄网站| 内地一区二区视频在线| 99香蕉大伊视频| 日本欧美视频一区| 一级片'在线观看视频| 日韩在线高清观看一区二区三区| 人人妻人人添人人爽欧美一区卜| 国产午夜精品一二区理论片| 国产免费视频播放在线视频| 夫妻性生交免费视频一级片| 亚洲精品久久午夜乱码| av有码第一页| 亚洲精品av麻豆狂野| 亚洲精品aⅴ在线观看| 下体分泌物呈黄色| 欧美人与性动交α欧美精品济南到 | √禁漫天堂资源中文www| 亚洲精品美女久久av网站| av卡一久久| 高清黄色对白视频在线免费看| 日韩 亚洲 欧美在线| 天天操日日干夜夜撸| 三上悠亚av全集在线观看| 国产极品粉嫩免费观看在线| 成年女人在线观看亚洲视频| 精品一区在线观看国产| 亚洲成色77777| 中文字幕精品免费在线观看视频 | 在现免费观看毛片| 国产福利在线免费观看视频| 国产成人欧美| 极品人妻少妇av视频| 男人添女人高潮全过程视频| 久久久精品94久久精品| 免费在线观看完整版高清| 咕卡用的链子| 街头女战士在线观看网站| 亚洲欧美色中文字幕在线| 日本爱情动作片www.在线观看| 亚洲一码二码三码区别大吗| 狠狠精品人妻久久久久久综合| 日韩av免费高清视频| 国产精品.久久久| 久久久久国产精品人妻一区二区| 久久久久久伊人网av| 久久亚洲国产成人精品v| 高清欧美精品videossex| 波多野结衣一区麻豆| 美女中出高潮动态图| 久久韩国三级中文字幕| 丝袜人妻中文字幕| 国产高清三级在线| 久久99热这里只频精品6学生| 边亲边吃奶的免费视频| av线在线观看网站| 欧美性感艳星| 亚洲第一区二区三区不卡| 晚上一个人看的免费电影| 亚洲欧洲日产国产| 老司机影院毛片| 欧美成人午夜精品| 日韩中文字幕视频在线看片| 少妇被粗大的猛进出69影院 | 男人爽女人下面视频在线观看| 最近最新中文字幕大全免费视频 | 精品人妻一区二区三区麻豆| 99久久中文字幕三级久久日本| 午夜影院在线不卡| 国产亚洲午夜精品一区二区久久| 赤兔流量卡办理| 毛片一级片免费看久久久久| 男女免费视频国产| 少妇熟女欧美另类| 免费观看a级毛片全部| 狂野欧美激情性bbbbbb| 制服诱惑二区| 国产日韩欧美亚洲二区| 亚洲精品色激情综合| 精品一品国产午夜福利视频| 午夜免费鲁丝| 在线精品无人区一区二区三| 亚洲色图综合在线观看| 一级片免费观看大全| 国产亚洲一区二区精品| 免费播放大片免费观看视频在线观看| 亚洲美女搞黄在线观看| 国产精品女同一区二区软件| 在线看a的网站| 麻豆精品久久久久久蜜桃| 亚洲综合色网址| 日韩成人伦理影院| 国产片特级美女逼逼视频| 大香蕉久久成人网| 丰满饥渴人妻一区二区三| 一区二区av电影网| 午夜影院在线不卡| 日日摸夜夜添夜夜爱| 在线看a的网站| 熟女人妻精品中文字幕| 十分钟在线观看高清视频www| 人妻 亚洲 视频| 国国产精品蜜臀av免费| 精品熟女少妇av免费看|