Identification of Produced Water and Characteristics of Hydrochemistry and Stable Hydrogen−Oxygen Isotopes of Contaminated Groundwater
-
摘要:
识别地下水污染来源、认识受该类污染源污染的地下水化学特征是地下水污染防治工作的重中之重。产出水作为石油、天然气工业的废水,具有组分复杂、危害性大的特点。针对受产出水污染的地下水研究较少,受污染地下水的特征以及识别该污染源的方法尚不明确的问题,笔者以延安某地下水污染场地为研究区,利用水文地球化学和氢氧稳定同位素的方法探讨受产出水污染的地下水的水化学和同位素特征,并通过对比地下水和油层水的钠氯系数、氯镁系数、脱硫系数和碳酸盐平衡系数对产出水进行识别。研究结果表明,该区域受产出水污染的地下水表现为高TDS和贫化的氢氧稳定同位素特征;其水化学类型以Cl−Na型、Cl−Mg·Ca·Na型为主,且随着受产出水影响程度降低,地下水由Cl−Na型转化为Cl−Mg·Ca·Na型,再到HCO3·SO4−Na·Ca·Mg型;离子比例关系较正常地下水混乱,无线性规律;受产出水污染的地下水的钠氯系数、氯镁系数、脱硫系数和碳酸盐平衡系数大小均在长6油层水的范围内,表明判断油气成藏条件的相关参数可以用来识别产出水污染。该研究探讨了受产出水污染的地下水的水化学特征和氢氧稳定同位素特征,提出对比地下水和油层水的相关参数来识别产出水污染的方法,对产出水污染场地的识别、认识、调查、监测、和修复具有重要意义。
Abstract:Identifying the source of pollution in groundwater and understanding the hydrochemical characteristics of contaminated groundwater by such pollution are very important for pollution prevention of groundwater. As the waste water of petroleum and natural gas industry, the produced water has the characteristics of complex components and great harmfulness. In view of the problems that there was less research on the contaminated groundwater by produced water, and the characteristics of the contaminated groundwater and the method to identify the pollution source were still unclear, the paper toke a polluted groundwater site in Yan’an as the research area, applied the methods of hydrogeochemistry and stable hydrogen−oxygen isotopes to describe the hydrochemical and isotopic characteristics of the groundwater contaminated by produced water, and compared the sodium−chloride coefficients, magnesium−chloride coefficients, desulfurization coefficients and carbonate balance coefficients of groundwater and reservoir water to identify the produced water. The results showed that the contaminated groundwater by produced water in this area was characterized by high TDS and depleted stable hydrogen−oxygen isotope; The chemical types were mainly Cl−Na type and Cl−Mg·Ca·Na type, and with the decrease of the influence of produced water, the groundwater changes from Cl−Na type to Cl−Mg·Ca·Na type and then to HCO3·SO4−Na·Ca·Mg type. The relationship of ion proportions was more chaotic than that of normal groundwater and had no linear law. The sodium−chloride coefficients, magnesium−chloride coefficients, desulfurization coefficients and carbonate balance coefficients of polluted groundwater by produced water were all within the range of Chang 6 reservoir water, indicating that the relevant parameters for judging the conditions of oil and gas accumulation could be used to identify produced water. This study described the hydrochemical and stable hydrogen−oxygen isotopic characteristics of contaminated groundwater by produced water, and proposed a method to identify produced water by comparing the relevant parameters of groundwater and reservoir water, which were of great significance to the identification, recognition, investigation, monitoring and repair of polluted sites by produced water.
-
青藏高原的形成与演化经历了多个洋盆的闭合以及陆陆碰撞过程,由此形成了高原上多个近EW向延伸的构造缝合带,将青藏高原划分为多个次级地块(Zhu et al., 2011; Kapp et al., 2019)。其中,班公湖–怒江缝合带(以下简称班–怒带)横亘于青藏高原中部,是中生代班公湖–怒江特提斯洋(以下简称班-怒洋)构造演化的残迹,其EW向延伸达2 000 km以上,构成了拉萨地块与羌塘地块之间的地质界线(图1a)。
图 1 青藏高原构造单元划分图(a)、南羌塘地块中—晚侏罗世侵入岩展布图(b)、卡易错地区地质简图(c)1.第四系;2.上三叠统日干配错组;3.古新统—始新统牛堡组;4.中—晚侏罗世花岗岩;5. 构造岩快;6. 断层;7. 角度不整合;8. 湖泊;9. 采样点;10. 锆石U-Pb年龄(本文);11. 锆石U-Pb年龄(引用);JSSZ.金沙江缝合带;LSSZ.龙木措–双湖缝合带;BNSZ.班公湖–怒江缝合带;IYZSZ.印度–雅鲁藏布缝合带;LT.拉萨地块;SQ.南羌塘地块;NQ.北羌塘地块;年龄数据引自Li等(2014)、Liu等(2014)、Wu等(2016)、Sun等(2020)、Yang等(2021)Figure 1. (a) Tectonic subdivision of the Tibetan Plateau, (b) distribution of Middle-Late Jurassic intrusive rocks within Southern Qiangtang block and (c) simplified geological map of the Kayico area, Tibet中生代以来,受班–怒洋俯冲闭合过程的影响,在班–怒带上及其两侧发育了大规模的火山岩浆活动,这些多期次岩浆岩记录了洋盆开合过程的信息,是反演区域构造–岩浆演化的关键,为揭示班–怒洋俯冲闭合过程提供了重要约束(李永飞等, 2005; Zhu et al., 2011, 2016; Pan et al., 2012; Wu et al., 2016, 2019a, 2019b; 刘海永等, 2019; 吴浩等, 2020)。然而,尽管国内外学者先后在青藏高原中部开展了大量的地质研究工作,但是关于区域上中生代多期次、多样性岩浆活动的成因与深部动力学过程尚存有较大的争议,是众多学者关注的热点(Kapp et al., 2007; Pan et al., 2012; Zhu et al., 2016; Fan et al., 2017)。近年来,越来越多的岩浆岩研究资料表明羌塘地块的南缘出露着大规模的中—晚侏罗世岩浆岩(图1b),并显示复杂的地球化学组成,仅在中酸性侵入岩中就先后识别了普通钙碱性I型花岗岩、高分异型花岗岩、富Na埃达克岩与富K埃达克岩等岩石类型(Li et al., 2014; Wu et al., 2018)。然而,羌塘地块南缘中—晚侏罗世多样性岩石类型之间究竟有何成因联系、形成于何种构造环境、反映了怎样的深部动力学过程尚不明确,亟待进一步研究。不仅如此,尽管关于班-怒洋的闭合时限仍存有争议(Kapp et al., 2007; Wu et al., 2019b; Fan et al., 2021),但是当前的研究普遍认为青藏高原中部中-晚侏罗世岩浆岩形成于洋壳俯冲背景,是班–怒洋洋壳俯冲消减引发的弧型岩浆活动(Li et al., 2014; Wu et al., 2016, 2018)。因此,查明羌塘地块南缘中—晚侏罗世岩浆作用过程,不仅对认识班–怒洋构造演化史具有重要的指示意义,同时对探讨俯冲带弧型岩浆起源与演化、壳幔物质循环与交换过程同样具有重要的约束。
藏北卡易错地区发育着大规模的晚侏罗世酸性侵入岩,为揭示班-怒洋俯冲过程、探讨俯冲带复杂的构造–岩浆活动提供了理想的研究对象(图1b)。本次在详细野外地质调查工作的基础上,对藏北卡易错地区出露的花岗闪长岩体进行了系统的岩石学、地质年代学、地球化学和同位素的研究工作,以此确定花岗闪长岩的形成时代与岩石成因,进一步对比区域岩浆岩研究资料,共同约束区域构造–岩浆过程,为揭示班–怒洋俯冲过程、探讨青藏高原早期形成与演化史提供新的约束。
1. 地质背景与岩石学特征
研究表明青藏高原自中生代以来经历了多个地块的闭合、碰撞过程,由此形成了青藏高原上近EW向延伸的多个构造缝合带,并将青藏高原从北至南划分为可可西里–松潘–甘孜地块、北羌塘地块、南羌塘地块、拉萨地块和喜马拉雅地块等多个次级地块(Zhu et al., 2011; Kapp et al., 2019)(图1a)。其中,班-怒带夹持于南羌塘地块和拉萨地块之间,大量的岩浆岩在南羌塘地块的南缘出露着大规模的中-晚侏罗世岩浆岩,该期岩浆岩以中酸性侵入岩为主,具有持续时间长、地球化学组分多样的特征,为探讨青藏高原早期形成与演化过程提供了理想的研究对象。
本次研究区卡易错地区位于日土县NE方向约45 km,大地构造位置处于班-怒带以北、南羌塘地块的南缘。区内构造–地层格架近NW–SE向展布,其研究区西南部主要以上三叠统日干配错组(T3r)灰岩夹砂岩和古新统—始新统牛堡组(E1-2n)砂、砾岩为主;而研究区东北部主要以酸性侵入岩为主。此外,受构造作用的影响,区内出露着大量的灰岩和砾岩的构造岩块(图1c)。前人已经对区内花岗岩体进行了初步的年代学和地球化学的研究工作,研究认为卡易错岩体形成于中—晚侏罗世(168~160 Ma),其地球化学组成指示岩石类型以高分异型花岗岩为主,是古老的变火成岩地壳部分熔融并经历广泛结晶分异作用的产物(Li et al., 2014; Liu et al., 2014)。
本次研究的花岗闪长岩(E 80°6′25″;N 33°28′20″)呈岩株状侵入于构造岩块之中(图2a),出露规模长约为10 m、宽约为3 m,岩石整体呈灰黑色,块状构造,中粗粒花岗结构(图2b),矿物组成以长石、角闪石、石英为主,粒度在0.5~2 mm之间,副矿物有锆石、磷灰石等(图2c、图2d)。
2. 分析方法
本次锆石U-Pb、全岩主微量地球化学与锆石Lu-Hf同位素测试分析工作均在武汉上谱分析科技有限责任公司完成。锆石U-Pb同位素定年和微量元素含量利用LA-ICP-MS同时分析完成,GeolasPro激光剥蚀系统由COMPexPro 102 ArF 193 nm准分子激光器和MicroLas光学系统组成,ICP-MS型号为Agilent 7700e,详细的仪器参数和分析流程见Zong等(2017)和李艳广等(2023)。分析数据的离线处理(包括对样品和空白信号的选择、仪器灵敏度漂移校正、元素含量及U-Pb同位素比值和年龄计算)采用软件ICPMSDataCal(Liu et al., 2008, 2010)完成。锆石样品的U-Pb年龄谐和图绘制和年龄加权平均计算采用Isoplot/Ex_ver3(Ludwig, 2003)完成。全岩主量元素含量利用日本理学PrimusⅡ X射线荧光光谱仪(XRF)分析完成,微量元素含量利用Agilent 7700e ICP-MS分析完成。原位微区锆石Lu-Hf同位素比值测试利用激光剥蚀多接收杯等离子体质谱(LA-MC-ICP-MS)完成。激光剥蚀系统为Geolas HD(Coherent,德国),MC-ICP-MS为Neptune Plus(Thermo Fisher Scientific,德国)。
3. 分析结果
3.1 锆石U-Pb定年与Lu-Hf同位素
花岗闪长岩锆石整体呈无色透明的长柱状,晶形完好,粒径为200~300 μm,长宽比为2∶1~3∶1(图3)。本次共对20颗锆石进行测试分析工作(表1),20颗锆石测点206Pb/238U年龄均集中在162~154 Ma之间。在谐和图上(图3),所有测点都落在谐和线上或附近区域,获得锆石206Pb/238U年龄加权平均值为(158.4±1.8)Ma(MSWD=0.15),这与前人在卡易错花岗岩体中获得的168~160 Ma的年龄信息基本一致,代表了花岗闪长岩的形成时代。
表 1 卡易错花岗闪长岩LA-ICP-MS锆石U-Pb定年分析结果Table 1. LA-ICP-MS U-Pb dating results for zircons of Kayico granodiorites点号 同位素比值(1σ) 年龄比值(Ma) 207Pb/206Pb 1σ 207Pb/235U 1σ 206Pb/238U 1σ 207Pb/206Pb 1σ 207Pb/235U 1σ 206Pb/238U 1σ 18T12-01 0.04954 0.00219 0.17017 0.00781 0.02491 0.00066 173 59 160 7 159 4 18T12-02 0.04929 0.00266 0.16824 0.00927 0.02475 0.00067 162 78 158 8 158 4 18T12-03 0.04965 0.00191 0.16944 0.00690 0.02475 0.00064 179 50 159 6 158 4 18T12-04 0.04989 0.00233 0.16974 0.00824 0.02467 0.00064 190 66 159 7 157 4 18T12-05 0.04939 0.00438 0.17064 0.01519 0.02505 0.00071 166 147 160 13 159 4 18T12-06 0.05092 0.00237 0.16985 0.00810 0.02419 0.00066 237 61 159 7 154 4 18T12-07 0.04885 0.00211 0.16623 0.00740 0.02467 0.00067 141 56 156 6 157 4 18T12-08 0.04955 0.00218 0.17037 0.00777 0.02493 0.00066 174 59 160 7 159 4 18T12-09 0.04888 0.00299 0.16655 0.01042 0.02471 0.00064 142 95 156 9 157 4 18T12-10 0.05012 0.00238 0.17289 0.00852 0.02501 0.00065 201 67 162 7 159 4 18T12-11 0.04908 0.00453 0.16867 0.01561 0.02492 0.00071 152 154 158 14 159 4 18T12-12 0.04889 0.00190 0.16671 0.00681 0.02472 0.00065 143 50 157 6 157 4 18T12-13 0.04940 0.00280 0.17026 0.00980 0.02499 0.00069 167 83 160 9 159 4 18T12-14 0.04970 0.00173 0.16964 0.00637 0.02475 0.00064 181 43 159 6 158 4 18T12-15 0.04905 0.00219 0.16888 0.00775 0.02496 0.00068 150 59 158 7 159 4 18T12-16 0.04950 0.00233 0.17051 0.00831 0.02498 0.00066 172 66 160 7 159 4 18T12-17 0.05017 0.00176 0.17246 0.00651 0.02492 0.00064 203 44 162 6 159 4 18T12-18 0.04892 0.00463 0.16803 0.01594 0.02490 0.00071 144 159 158 14 159 4 18T12-19 0.05000 0.00226 0.17212 0.00808 0.02496 0.00066 195 62 161 7 159 4 18T12-20 0.04862 0.00207 0.17018 0.00751 0.02538 0.00068 130 56 160 7 162 4 此外,对8颗获得谐和年龄的锆石进行原位Lu-Hf同位素测试(表2),8颗锆石测点的176Yb/177Hf在
0.012865 ~0.037412 之间,176Lu/177Hf在0.000427 ~0.001145 之间,表明所测锆石放射成因Hf的积累很少(吴福元等, 2007),测定的176Hf/177Hf(0.282602 ~0.282669 )值可以用来代替锆石的初始176Hf/177Hf值。8个测点Hf同位素地壳模式年龄TDMC为1090 ~1222 Ma,对应的εHf(t)值变化范围在−2.66~−0.27之间。表 2 卡易错花岗闪长岩锆石Lu-Hf同位素组成Table 2. Lu-Hf isotopes of zircons from the Kayico granodiorites.点号 年龄 (Ma) 176Hf/177Hf 1σ 176Lu/177Hf 1σ 176Yb/177Hf 1σ εHf(0) 1σ εHf(t) 1σ TDM1 TDMC fLu/Hf 18T12-01 159 0.282608 0.000008 0.000562 0.000003 0.017391 0.000117 −5.81 0.59 −2.38 0.60 901 1207 −0.98 18T12-02 158 0.282602 0.000009 0.000996 0.000005 0.031211 0.000271 −6.03 0.61 −2.66 0.62 921 1222 −0.97 18T12-03 158 0.282609 0.000010 0.000842 0.000026 0.026596 0.000831 −5.77 0.61 −2.39 0.62 907 1207 −0.97 18T12-04 154 0.282618 0.000010 0.000867 0.000013 0.027469 0.000376 −5.46 0.62 −2.17 0.63 895 1191 −0.97 18T12-05 159 0.282641 0.000009 0.000427 0.000002 0.012865 0.000067 −4.65 0.61 −1.20 0.62 853 1142 −0.99 18T12-06 157 0.282643 0.000009 0.000608 0.000005 0.019402 0.000131 −4.57 0.60 −1.18 0.61 854 1139 −0.98 18T12-07 159 0.282631 0.000008 0.000532 0.000001 0.016431 0.000077 −4.98 0.59 −1.55 0.60 868 1161 −0.98 18T12-08 159 0.282669 0.000009 0.001145 0.000043 0.037412 0.001346 −3.64 0.60 −0.27 0.61 829 1090 −0.97 3.2 地球化学
本次共采集4件花岗闪长岩样品进行全岩主微量元素地球化学分析工作,分析结果见表3。4件样品的地球化学组成相对均一,其SiO2含量为62.6%~65.2%,Al2O3含量为15.9%~16.6%,TiO2含量为0.68%~0.81%,TFe2O3含量为4.53%~5.37%,MgO含量为2.06%~2.34%,全碱(Na2O+K2O)含量为5.40%~5.56%,Na2O/K2O值为0.94~1.10,Mg#为50~52。在岩石类型判别图解中(图4a),样品均在花岗闪长岩区域;在K2O-SiO2图解中(图4b),样品显示中钾–高钾钙碱性的特征;在A/NK-A/CNK图解中(图3c),样品整体显示弱过铝质的特征(A/CNK=1.02~1.05)。以上主量元素地球化学组成表明样品整体显示弱过铝质中钾-高钾钙碱性花岗闪长岩的特征。
表 3 卡易错花岗闪长岩全岩主量(%)和微量(10−6)元素分析结果Table 3. Whole-rock major (%) and trace (10−6) element contents of Kayico granodiorites元素 T12h1 T12h2 T12h3 T12h4 元素 T12h1 T12h2 T12h3 T12h4 SiO2 63.6 65.2 64.5 62.6 Zr 248 230 244 245 TiO2 0.78 0.68 0.71 0.81 Nb 13.8 13.1 13.7 14.1 Al2O3 16.1 15.9 16.2 16.6 Sn 3.58 5.26 4.55 4.30 TFe2O3 5.37 4.53 4.82 5.12 Cs 17.5 21.1 14.6 20.0 MnO 0.08 0.06 0.07 0.07 Ba 457 383 401 428 MgO 2.34 2.06 2.07 2.33 La 29.8 32.7 32.1 30.2 CaO 4.48 4.44 4.38 4.83 Ce 66.7 68.5 63.0 61.6 Na2O 2.70 2.74 2.83 2.74 Pr 6.97 7.45 7.02 6.90 K2O 2.86 2.78 2.57 2.77 Nd 26.37 28.1 27.3 25.4 P2O5 0.15 0.13 0.14 0.15 Sm 5.26 5.68 5.75 5.24 LOI 0.99 1.11 1.30 1.30 Eu 1.11 1.05 1.14 1.16 SUM 99.5 99.7 99.5 99.3 Gd 5.08 4.65 4.55 4.66 Li 56.7 49.2 63.6 54.6 Tb 0.83 0.83 0.78 0.74 Be 1.97 2.31 2.18 2.07 Dy 4.95 5.01 4.54 4.49 Sc 14.9 12.7 13.8 14.8 Ho 1.00 1.01 0.91 0.92 V 74.0 61.3 67.1 74.1 Er 2.83 2.87 2.72 2.60 Cr 40.0 32.8 37.6 37.7 Tm 0.45 0.43 0.40 0.41 Co 15.4 11.5 12.9 12.8 Yb 2.70 2.78 2.46 2.42 Ni 16.6 14.3 15.2 15.5 Lu 0.41 0.44 0.37 0.37 Cu 104 45.7 44.1 45.7 Hf 6.64 6.24 6.79 6.29 Zn 56.6 43.1 45.7 51.6 Ta 1.02 1.10 1.05 0.98 Ga 19.2 18.8 19.1 19.5 Tl 1.33 1.40 1.13 1.34 Rb 170 185 169 178 Pb 8.40 9.57 7.76 10.2 Sr 188 181 194 203 Th 14.4 16.9 15.5 12.4 Y 28.2 28.0 26.8 24.8 U 2.18 2.88 2.52 2.23 图 4 卡易错花岗岩TAS图解(Middlemost, 1994)(a)、K2O-SiO2图解(Le Maitre et al., 1989; Rickwood, 1989)(b)和A/NK-A/CNK图解(c)(Shand, 1943)Figure 4. (a) TAS classification diagram, (b) K2O vs. SiO2 diagram and (c) A/NK vs. A/CNK diagram of Kayico granitic rocks在球粒陨石标准化稀土元素配分曲线中,样品呈轻稀土元素富集、重稀土元素亏损的右倾模式([La/Yb]N=7.90~9.34),同时具有不同程度的Eu负异常(Eu/Eu*=0.62~0.72)(图5a)。在原始地幔标准化蛛网图中,样品具有明显Nb、Ta等高场强元素以及Ba、Sr、Eu等大离子亲石元素的亏损(图5b)。
图 5 岩石球粒陨石标准化稀土元素模式图(a)和原始地幔标准化微量元素蛛网图(b)(标准化值引自Sun 等1989)Figure 5. (a) Chondrite-normalized rare earth element and (b) primitive-mantle–normalized multi-element patterns4. 讨论
4.1 中-晚侏罗世岩浆活动
近年来,国内外学者已经对青藏高原中部出露的岩浆岩开展了系统的研究工作,报道了大量的年代学和地球化学数据(Li et al., 2014; Liu et al., 2014; Wu et al., 2016; 刘海永等, 2022)。越来越多的研究表明在南羌塘地块南缘的卡易错、材玛、青草山、改则、荣玛、高保约一带存在一期大规模的中—晚侏罗世岩浆活动(约165~150 Ma),这些岩浆岩整体呈带状近平行于班–怒带出露,EW向延伸近上千千米,岩石类型整体以中酸性侵入岩为主,具有持续时间长、分布范围广的特征(图1b)。前人的研究在卡易错花岗岩体中获得了168~160 Ma的锆石U-Pb年龄信息(Liu et al., 2014; Li et al., 2014),本次锆石U-Pb定年工作在卡易错花岗闪长岩中获得了158 Ma的年龄,表明区内花岗岩形成于中—晚侏罗世之交,与南羌塘地块上大规模发育的中酸性侵入岩形成时代相近,应该是区域上同一期构造–岩浆活动的产物。
此外,Wu等(2018)对南羌塘地块上发育的中—晚侏罗世侵入岩进行了系统的地球化学对比研究,研究发现该期侵入岩整体显示复杂的地球化学特征,根据岩石的主、微量元素含量和同位素地球化学特征,可以划分为普通钙碱性I型花岗岩、高分异花岗岩、富Na埃达克岩与富K埃达克岩等不同的岩石类型,反应了区域上花岗质岩石复杂的岩浆源区与成岩过程。而根据现有地球化学资料,卡易错花岗岩可以划分为普通钙碱性I型花岗岩和高分异型花岗岩两类(Li et al., 2014; Liu et al., 2014),然而二者之间有何成因联系尚不明确。
4.2 构造背景
大量的岩浆岩研究资料表明青藏高原中部中生代岩浆作用整体沿班–怒带展布,主要发育在混杂带上及其两侧的南羌塘地块南缘和拉萨地块北缘,显示与与班–怒洋密切的时空联系,系统的研究工作已经初步建立中生代多期次岩浆作用与班–怒洋俯冲闭合过程之间的成因联系(Zhu et al., 2011, 2016; Kapp et al., 2019; Wu et al., 2019a, 2019b)。尽管对于班–怒洋的闭合时限一直存在争论,然而现有的争议均认为班–怒洋闭合与拉萨–羌塘地块陆陆碰撞应不早于白垩纪,而青藏高原中部侏罗纪构造演化主要受班–怒洋俯冲消减作用的影响(Liu et al., 2022)。此外,班–怒带上蛇绿混杂岩系统的年代学和地球化学的研究同样表明班–怒洋洋盆在中—晚侏罗世仍存有一定的规模(范建军等, 2019; 李志军等, 2019; 唐跃等, 2021)。
卡易错花岗闪长岩显著的Nb、Ta等元素亏损(图5),显示与俯冲成因的弧型岩浆岩相似的地球化学特征。近年来,不同学者对南羌塘地块上发育的中—晚侏罗世岩浆岩开展了大量的研究工作,陆续取得了众多的研究进展(Liu et al., 2014; Li et al., 2016; Wu et al., 2016)。其中,中酸性侵入岩整体显示陆缘弧的地球化学特征,而洋壳熔融成因的埃达克质花岗岩的识别则为深部俯冲洋壳的存在提供了最直接的岩石学证据(Li et al., 2016)。此外,区域上还陆续报道了中—晚侏罗世俯冲成因的钙碱性弧型安山岩与OIB型辉绿岩(李小波等, 2015; 董宇超等, 2016; Li et al., 2016),如此复杂多样的岩石组合反映了俯冲带上复杂的源区物质组成与循环过程。不仅如此,最近的研究提出班–怒洋的初始俯冲起始于晚三叠世末—早侏罗世(Qian et al., 2020; Liu et al., 2022),进一步表明班–怒洋中—晚侏罗世应处于洋壳俯冲背景。综上所述,南羌塘地块南缘出露的中—晚侏罗世岩浆岩应该形成于洋壳俯冲背景,是班-怒洋北向俯冲至南羌塘地块之下引发的弧型岩浆活动。
4.3 岩石类型与成因
花岗岩的类型划分与岩石成因长期以来一直以来是众多地质学家关注的热点研究问题,其中Chappell等(1974, 1992)根据花岗岩的岩浆源区物质组成和成岩构造环境的差异将花岗岩划分为I、S、M、A型4类,该分类方式已经被广泛运用于花岗岩成因与演化的研究工作(王亮等, 2022; 孙巍等, 2022)。Li等(2014)和Liu等(2014)对区内花岗岩岩石开展了年代学和地球化学的研究工作,识别了普通钙碱性I型和高分异型两类花岗岩。笔者在新获得的卡易错花岗闪长岩研究资料的基础上,进一步收集整理了前人已报道的卡易错花岗岩体的数据资料,以此准确约束卡易错地区花岗质岩石的岩石类型与成因。
卡易错花岗岩整体具有低的Zr+Ce+Nb+Y含量,然而根据其(Na2O+K2O)/CaO和FeOt/MgO值的不同,区内岩石可以划分为高分异与未分异两组(图6a、图6b )。不同于前人报道的高分异型花岗岩,本次采集的花岗闪长岩样品显示低的(Na2O+K2O)/CaO和FeOt/MgO值,具有未分异花岗岩的特征。在P2O5-SiO2和Th-Rb图解中,花岗闪长岩样品均显示出I型花岗岩的演化趋势(图6c、图6d)。结合花岗闪长岩弱过铝质的特征(图4c),卡易错花岗闪长岩应该属于普通钙碱性I型花岗岩。I型花岗岩一般认为起源于变火成岩下地壳部分熔融或者幔源物质对变沉积岩下地壳的改造(Petford et al., 1996; Chappell et al., 2001; Li et al., 2007),研究区内乃至整个南羌塘地块南缘均未发现大规模幔源岩浆活动成因的基性岩,结合花岗闪长岩相对均一的锆石原位Hf同位素组成(εHf(t)=−2.66~−0.27),卡易错花岗闪长岩难以解释为幔源物质与变沉积岩下地壳熔体混合的产物。同时,花岗闪长岩具有高的CaO/Na2O以及低的Al2O3/TiO2和Rb/Ba、Rb/Sr值,进一步指示着其起源于玄武质火成岩下地壳的部分熔融(图7a、 图7b)。此外,花岗闪长岩中不同程度的Eu、Sr、Ba等元素的负异常一般认为是成岩过程中存在着长石类矿物的结晶分异,而Nb、Ta等元素的亏损则和金红石/榍石的结晶分离相关。综上所述,卡易错花岗闪长岩应该是南羌塘地块之下古老的变火成岩下地壳熔体经历一定结晶分异作用的产物。
图 6 卡易错花岗岩(Na2O/K2O)/CaO-Zr+Ce+Nb+Y图解(a)、FeOt/MgO-Zr+Ce+Nb+Y图解(b)(Whalen et al., 1987)、P2O5-SiO2图解(c)和Th-Rb图解(d)(Li et al., 2007)Figure 6. Geochemical classification diagrams of Kayico granitic rocks: (a)(Na2O/K2O)/CaO vs. Zr+Ce+Nb+Y diagram, (b) FeOt/MgO vs. Zr+Ce+Nb+Y diagram, (c) P2O5 vs. SiO2diagram, (d) Th vs. Rb diagram of Kayico granitic rocks图 7 卡易错花岗岩CaO/Na2O-Al2O3/TiO2图解(a)、Rb/Ba-Rb/Sr图解(b)(Sylvester, 1998)、 Rb/Sr-Sr图解(c)、Ba-Sr图解(d)(Rollinson, 1993)Amp.闪石;Bi.黑云母;Kfs.钾长石;Pl.斜长石Figure 7. (a) CaO/Na2O vs. Al2O3/TiO2 diagram, (b) Rb/Ba vs. Rb/Sr diagram, (c) Rb/Sr vs. Sr diagram, (d) Ba vs. Sr diagram of Kayico granitic rocks卡易错高分异型花岗岩具有与钙碱性I型花岗岩形成时代一致,在空间上密切共生,并相似的锆石Hf同位素特征,表明二者具有相似的岩浆源区物质组成(Li et al., 2014; Liu et al., 2014)。然而,高分异型花岗岩显示强烈的亏损Eu、Sr、Ba等元素(图5),指示着岩浆在侵位过程中经历了强烈的结晶分异作用。在Rb/Sr-Sr和Ba-Sr图解中,卡易错两类花岗岩之间显示出明显的长石类矿物分离结晶趋势(图7c、图7d),进一步表明卡易错不同类型花岗岩地球化学的差异是后期岩浆侵位过程中经历不同演化过程的产物。近年来,晶粥体模型(MUSH)被广泛应用于解释同期共生花岗闪长岩与高分异型花岗岩之间的成因联系(Bachmann et al., 2004; Hildreth, 2004),该模型提出花岗质岩浆在浅层岩浆房中发生矿物结晶时,形成一种晶体与液体共存的晶粥体,其外围的物质结晶形成富含斑晶的花岗质岩石,而中心的残余岩浆则形成高分异的高硅花岗岩(Wu et al., 2017)。晶粥体模型无疑为卡易错地区空间上相伴生、时代上相一致的普通钙碱性I型花岗岩和高分异型花岗岩提供了合理的成因解释。因此,研究认为,在中—晚侏罗世班–怒洋持续北向俯冲过程中,底侵的幔源玄武质岩浆诱发南羌塘地块之下古老的变火成岩下地壳物质发生重熔,形成的熔体上升侵位、并在浅层岩浆房内发生显著低压结晶分异作用,其外围的晶粥体与内部的熔体分别冷凝形成了卡易错地区普通钙碱性I型和高分异型两类花岗质岩石(图8)。
5. 结论
(1)锆石U-Pb年龄表明,卡易错花岗闪长岩形成于158 Ma,与前人在卡易错岩体中获得的168~160 Ma的年龄相一致,是南羌塘地块南缘中—晚侏罗世构造-岩浆活动的产物。
(2)卡易错花岗闪长岩具有钙碱性I型花岗岩的地球化学特征,以及相对富集的锆石Hf同位素组成,是班-怒洋中-晚侏罗世北向俯冲背景下南羌塘地块古老的变火成岩下地壳熔融的产物。
(3)结合现有研究资料,晶粥体模式为卡易错地区钙碱性I型与高分异型两类花岗岩提供了合理的成因解释,二者是壳源熔体在浅层岩浆房内经历结晶分异后不同端元冷凝的产物。
致谢:野外工作得到西藏自治区地质调查院刘海永博士和吉林大学罗安波博士等人的帮助,审稿专家对稿件提出的宝贵意见对论文质量提高至关重要,在此一并致以衷心的感谢。
-
表 1 主要离子和TDS浓度统计表(mg/L)
Table 1 The concentrations of major ions and TDS (mg/L)
样品 时间 Na+ Ca2+ Mg2+ K+ Cl− HCO3− SO42− TDS Z1 2018 2 445 1 467 103 7.78 6 795 84.6 10.9 10 871 2019 1 150 627 37 6.63 2 853 85.4 158.2 4 875 Z2 2018 1 308 702 69 6.11 3 540 201 7.99 5 734 2019 669 312 26 8.04 1 530 293 40.8 2 733 2020 164 106 15.4 3.86 177 511 8.8 731 Z3 2018 1 368 567 73 5.47 3 391 153 7.7 5 488 2019 1 484 668 49 15.91 3 721 159 83.1 6 100 2020 1 460 616 55.8 4.26 3 250 149 1.3 5 462 Z4 2019 262 171 146 5.47 639 500 104.8 1 579 2020 286 207 166 3.31 778 665 94 1 867 Z5 2018 304 99 27 1.53 512 238 96.6 1 180 2019 329 112 30 1.98 542 287 89.1 1 247 Z6 2018 258 18 67 21.48 155 513 196.7 1 003 2019 303 104 90 25.53 195 824 228.6 1 358 2020 307 111 89 21.8 272 814 237 1 445 Z7 2018 241 69 32 2.42 236 370 131.3 912 2019 150 45 29 1.62 78 427 75.6 592 2020 167 51.9 28.4 1.80 110 437 83.2 661 H1 2018 59.5 33.2 18.8 4.66 35 196 59.6 309 2019 135.1 52.4 39.1 10.89 107.2 305.1 133.3 631 2020 144 45.1 34.1 9.99 131 271 138 638 H2 2018 69.3 29.4 18.6 3.54 45.6 196 64.2 329 2019 157.7 41 39.8 10.61 119.2 311.2 145.7 670 2020 149 43.4 33.7 8.63 121 295 143 646 注:Z代表地下水,H代表地表水。 表 2 2018~2020年样品 δD(‰)和 δ18O(‰)值
Table 2 The values of δD (‰) and δ18O (‰) from 2018 to 2020
采样时间 2018年 2019年 2020年 样品 δD(‰) δ18O (‰) δD (‰) δ18O (‰) δD (‰) δ18O (‰) Z1 −84.08 −11.32 −74.78 −9.65 − − Z2 −89.64 −11.74 −80.18 −10.67 −70.65 −9.19 Z3 −84.28 −11.53 −88.33 −12.03 − − Z4 − − −65.86 −8.55 −65.39 −8.60 Z5 −72.78 −9.47 −71.28 −9.84 − − Z6 −64 −8.72 −61.34 −8.43 −63.40 −8.66 Z7 −62.32 −8 −60.81 −8.21 −61.60 −7.69 H1 −53.71 −7.38 −56.08 −7.26 −55.36 −7.30 H2 −51.52 −7.39 −55.99 −6.98 −52.36 −7.13 表 3 理想混合模型计算结果
Table 3 The results of ideal mixing
2018~2019年计算结果(mg/L) 混合分数 Na+ Ca2+ Mg2+ K+ Cl− HCO3− SO42- NO3− 以Na+ 计,x=0.5403 669.01 343.21 48.31 6.91 1 678.80 236.87 67.22 8.77 以Cl− 计,x=0.5835 617.92 314.52 46.66 6.98 1 529.99 239.74 71.96 9.23 实测值 669 312 26 8.04 1 530 293 40.8 1.86 2019~2020年计算结果(mg/L) 混合分数 Na+ Ca2+ Mg2+ K+ Cl− HCO3− SO42- NO3− 以Na+ 计,x=0.9289 164.0 57.43 30.38 7.63 197.26 269.22 112.16 12.77 以Cl− 计,x=0.9430 156.33 53.56 30.44 7.62 177.03 268.86 113.24 12.93 实测值 164 106 15.4 3.86 177 511 8.8 0.14 -
焦艳军, 王广才, 崔霖峰, 等. 济源盆地地表水和地下水的水化学及氢、氧同位素特征[J]. 环境化学, 2014, 33(06): 962-968 doi: 10.7524/j.issn.0254-6108.2014.06.023 JIAO Yanjun, WANG Guangcai, CUI Linfeng, et al. Characteristics of hydrochemistry and stable hydrogen, oxygen isotopes in surface water and groundwater in Jiyuan Basin[J]. Environmental Chemistry, 2014, 33(06): 962-968. doi: 10.7524/j.issn.0254-6108.2014.06.023
梁晓伟, 牛小兵, 李卫成, 等. 鄂尔多斯盆地油田水化学特征及地质意义[J]. 成都理工大学学报(自然科学版), 2012, 39(05): 502-508 LIANG Xiaowei, NIU Xiaobing, LI Weicheng, et al. Chemical character of oil-field water in Ordos Basin and geological significance[J]. Journal of Chengdu University of Technology(Science & Technology Edition), 2012, 39(05): 502-508.
楼章华, 苏一哲, 朱蓉, 等. 四川盆地新场构造带上三叠统须家河组二段地层水化学动态特征及其成因[J]. 石油与天然气地质, 2021, 42(04): 841-851 doi: 10.11743/ogg20210406 LOU Zhanghua, SU Yizhe, ZHU Rong, et al. Dynamic chemical characteristics and origin of formation water in the second member of Xujiahe Formation, Xinchang structural belt, Sichuan Basin[J]. Oil & Gas Geology, 2021, 42(04): 841-851. doi: 10.11743/ogg20210406
张茂省, 孙传尧, 校培喜, 等. 延安市宝塔区地质灾害详细调查示范[J]. 西北地质, 2007,40(02): 29-55 doi: 10.3969/j.issn.1009-6248.2007.02.002 ZHANG Maosheng, SUN Chuanyao, XIAO Peixi, et al. A Demonstration Project for Detailed Geo-hazard Survey in Baota District, Yan’ an[J]. Northwestern Geology, 2007,40(02): 29-55. doi: 10.3969/j.issn.1009-6248.2007.02.002
曾溅辉, 吴琼, 杨海军, 等. 塔里木盆地塔中地区地层水化学特征及其石油地质意义[J]. 石油与天然气地质, 2008(02): 223-229 doi: 10.3321/j.issn:0253-9985.2008.02.011 ZENG Jianhui, WU Qiong, YANG Haijun, et al. Chemical characteristics of formation water in Tazhong area of the Tarim Basin and their petroleum geological significance[J]. Oil & Gas Geology, 2008(02): 223-229. doi: 10.3321/j.issn:0253-9985.2008.02.011
战常武. 曲塘次凹油田地层水特征及其石油地质意义[J]. 承德石油高等专科学校学报, 2020, 22(04): 21-23+34 doi: 10.3969/j.issn.1008-9446.2020.04.006 ZHAN Changwu. Characteristics of Formation Water and Its Petroleum Geological Significance in Qutang Sub Sag[J]. Journal of Chengde Petroleum College, 2020, 22(04): 21-23+34. doi: 10.3969/j.issn.1008-9446.2020.04.006
张景涛, 史浙明, 王广才, 等. 柴达木盆地大柴旦地区地下水水化学特征及演化规律[J]. 地学前缘, 2021, 28(04): 194-205 doi: 10.13745/j.esf.sf.2020.6.40 ZHANG Jingtao, SHI Zheming, WANG Guangcai, et al. Hydrochemical characteristics and evolution of groundwater in the Dachaidan area, Qaidam Basin[J]. Earth Science Frontiers, 2021, 28(04): 194-205. doi: 10.13745/j.esf.sf.2020.6.40
张俊, 尹立河, 顾小凡, 等. 同位素水化学指示的新疆孔雀河流域地下水与地表水关系[J]. 西北地质, 2021, 54(01): 185-195 doi: 10.19751/j.cnki.61-1149/p.2021.01.016 ZHANG Jun, YIN Lihe, GU Xiaofan, et al. Study on the Relationship Between Groundwater and Surface Water in Xinjiang Kongque River Basin Using Isotopes and Hydrochemistry method[J]. Northwestern Geology, 2021, 54(01): 185-195. doi: 10.19751/j.cnki.61-1149/p.2021.01.016
张治波, 刘腾, 李丽荣, 等. 鄂尔多斯盆地CD区块长6地层水化学性质及其地质意义[J]. 矿物岩石, 2017, 37(03): 61-68 doi: 10.19719/j.cnki.1001-6872.2017.03.009 ZHANG Zhibo, LIU Teng, LI Lirong, et al. Chemical Characteristics and Geological Significance of Chang 6 Formation Water in CD Area of Ordos Basin[J]. Journal of Mineralogy and Petrology, 2017, 37(03): 61-68. doi: 10.19719/j.cnki.1001-6872.2017.03.009
周新平, 邓秀芹, 李士祥, 等. 鄂尔多斯盆地延长组下组合地层水特征及其油气地质意义[J]. 岩性油气藏, 2021, 33(01): 109-120 doi: 10.12108/yxyqc.20210111 ZHOU Xinping, DENG Xiuqin, LI Shixiang, et al. Characteristics of formation water and its geological significance of lower combination of Yanchang Formation in Ordos Basin[J]. Lithologic Reservoirs, 2021, 33(01): 109-120. doi: 10.12108/yxyqc.20210111
周训, 胡伏生, 何江涛, 等. 地下水科学概论[M]. 北京: 地质出版社, 2014 ZHOU Xun, HU Fusheng, HE Jiangtao, et al. Introduction to Groundwater Science[M]. Beijing: Geological Publishing House, 2014.
Al-Ghouti M A, Al-Kaabi M A, Ashfaq M Y, et al. Produced water characteristics, treatment and reuse: A review[J]. Journal of Water Process Engineering, 2019, 28: 222-239. doi: 10.1016/j.jwpe.2019.02.001
Amakiri K T, Canon A R, Molinari M, et al. Review of oilfield produced water treatment technologies[J]. Chemosphere, 2022, 298: 134064. doi: 10.1016/j.chemosphere.2022.134064
Atekwana E A, Seeger E J. Carbonate and carbon isotopic evolution of groundwater contaminated by produced water brine with hydrocarbons[J]. Applied Geochemistry, 2015, 63: 105-115. doi: 10.1016/j.apgeochem.2015.08.001
Benko K L, Drewes J E. Produced Water in the Western United States: Geographical Distribution, Occurrence, and Composition[J]. Environmental engineering science, 2008, 25(2): 239–246. doi: 10.1089/ees.2007.0026
Chen X L, Sheng Y Z, Wang G C, et al. Microbial compositional and functional traits of BTEX and salinity co-contaminated shallow groundwater by produced water[J]. Water Research, 2022, 215.
Craig H. Isotopic Variations in Meteoric Waters[J]. Science (New York, N. Y. ), 1961, 133(3465): 1702-1703. doi: 10.1126/science.133.3465.1702
Dolan F C, Cath T Y, Hogue T S. Assessing the feasibility of using produced water for irrigation in Colorado[J]. Science of the Total Environment, 2018, 640-641: 619-628. doi: 10.1016/j.scitotenv.2018.05.200
Fakhru’l-Razi A, Pendashteh A, Abdullah L C, et al. Review of technologies for oil and gas produced water treatment[J]. Journal of Hazardous Materials, 2009, 170(2): 530-551.
Huang X J, Wang G C, Liang X Y, et al. Hydrochemical and Stable Isotope (δD and δ18O) Characteristics of Groundwater and Hydrogeochemical Processes in the Ningtiaota Coalfield, Northwest China[J]. Mine Water and the Environment, 2018, 37(1): 119-136. doi: 10.1007/s10230-017-0477-x
Jiang W B, Xu X S, Hall R, et al. Characterization of produced water and surrounding surface water in the Permian Basin, the United States[J]. Journal of Hazardous Materials, 2022, 430.
McDevitt B, McLaughlin M C, Blotevogel J, et al. Oil & gas produced water retention ponds as potential passive treatment for radium removal and beneficial reuse[J]. Environmental science. Processes & impacts, 2021, 23(3): 501-518.
Miller H, Dias K, Hare H, et al. Reusing oil and gas produced water for agricultural irrigation: Effects on soil health and the soil microbiome[J]. Science of the Total Environment, 2020, 722.
Ozgun H, Ersahin M E, Erdem S, et al. Effects of the pre-treatment alternatives on the treatment of oil-gas field produced water by nanofiltration and reverse osmosis membranes[J]. Journal of Chemical Technology & Biotechnology, 2013, 88(8): 1576-1583.
Patterson L A, Konschnik K E, Wiseman H, et al. Unconventional oil and gas spills: risks, mitigation priorities, and state reporting requirements[J]. Environmental science & technology, 2017, 51(5): 2563–2573.
Sanchez-Rosario R, Hildenbrand Z L. Produced Water Treatment and Valorization: A Techno-Economical Review[J]. Energies, 2022, 15(13): 4319.
Shariq L. Health Risks Associated With Arsenic and Cadmium Uptake in Wheat Grain Irrigated With Simulated Hydraulic Fracturing Flowback Water[J]. Journal of Environmental Health, 2019, 81(6): E1-E9.
Shih J S, Saiers J, Anisfeld S, et al. Characterization and analysis of liquid waste from Marcellus Shale gas development[J]. Abstracts of Papers of the Amercian Chemical Society, 2015, 49(16): 9557-9565.
Shores A, Laituri M. The state of produced water generation and risk for groundwater contamination in Weld County, Colorado[J]. Environmental science and pollution research international, 2018, 25(30): 30390-30400. doi: 10.1007/s11356-018-2810-8
Warner N R, Jackson R B, Darrah T H, et al. Geochemical evidence for possible natural migration of Marcellus Formation brine to shallow aquifers in Pennsylvania. [J]. Proceedings of the National Academy of Sciences of the United States of America, 2012, 109(30): 11961-11966. doi: 10.1073/pnas.1121181109
Xu B Y, Wang G C. Surface water and groundwater contaminations and the resultant hydrochemical evolution in the Yongxiu area, west of Poyang Lake, China[J]. Environmental Earth Sciences, 2016, 75(3): 1-16.
Zheng Z X, Zhang H D, Chen Z Y, et al. Hydrogeochemical and Isotopic Indicators of Hydraulic Fracturing Flowback Fluids in Shallow Groundwater and Stream Water, derived from Dameigou Shale Gas Extraction in the Northern Qaidam Basin. [J]. Environmental science & technology, 2017, 51(11): 5889-5898.
-
期刊类型引用(9)
1. 解铭威,周慧荻,陈耸,王向荣. 银川市生态系统服务价值评估及多情景模拟. 水土保持研究. 2025(01): 294-304 . 百度学术
2. 洪桦,张渝. 基于土地利用变化的生态系统服务价值演变研究——以福州滨海新城为例. 农业与技术. 2024(04): 107-112 . 百度学术
3. 龚芯磊 ,陈鸿申 ,胡世敏 ,李正荣 ,聂坤 ,梁龙 . 特定模型与GIS结合的生态系统服务功能评价及分区研究——以黔北芙蓉江流域绥正道地区为例. 贵州地质. 2024(03): 328-338 . 百度学术
4. 陈彦珍,张丽华,王亚军,宋学云,曾乐,郝爱华,贺甜田. 不同枸杞品种的区域差异性研究. 宁夏农林科技. 2024(07): 46-52 . 百度学术
5. 张永生,李建国,赵广臣. 农田回归用水引入鸣翠湖后水质改善的数值模拟研究. 水电能源科学. 2024(11): 51-54 . 百度学术
6. 朱青青,刘超,沈艳,马红彬,谭松伟,王国会,李燕,李千飞,李国强. 宁夏罗山草地生态系统服务价值的地形效应. 应用生态学报. 2024(12): 3267-3274 . 百度学术
7. 欧阳渊,刘洪,张景华,唐发伟,张腾蛟,黄勇,黄瀚霄,李富,陈敏华,宋雯洁. 西南山区生态地质调查技术方法研究. 西北地质. 2023(04): 218-242 . 本站查看
8. 江山,石绍山,郭常来,冯雨林,孙家全,孙秀波,周丽. 大凌河流域1998~2019年NDVI时空变化及其对气温和降水的响应. 西北地质. 2023(04): 254-262 . 本站查看
9. 高媛,李谦. 湿地生态系统服务价值评估研究现状及展望. 绿色科技. 2023(18): 62-68 . 百度学术
其他类型引用(1)