Hydrochemical Zonality and Geochemical Modeling of Deep-Lying Pore Water in Taiyuan City
-
摘要: 太原市深层孔隙水具有明显的水化学分带性, 具体表现为由山前到盆地依次分布硫酸-重碳酸型水、重碳酸-硫酸型水、重碳酸型水, 且各类地下水均大体在南北向上呈条带状展布, 这与补给水的水化学状况密切相关.利用地球化学模拟软件PHREEQC建立一系列地下水混合模型对深层孔隙水的水化学形成过程进行模拟, 结果显示: 盆地北部的深层孔隙水受到北部边山岩溶水、盆地北部浅层孔隙水、汾河水的补给, 其中北部边山岩溶水是最主要的补给源; 盆地西部的深层孔隙水由西部边山岩溶水与盆地西部浅层孔隙水混合而成; 盆地南部的深层孔隙水则由盆地北部与西部的深层水混合而成.混合作用是控制区域水化学状况的最重要的因素.Abstract: Hydrochemical zonality was delineated for deep-lying pore water in Taiyuan. SO4-HCO3 water, HCO3-SO4 water and HCO3 water occur sequentially from mountain area to basin area, in north to south zonal distribution, which is highly correlative with the hydrochemistry of supply water. To understand the forming process of deep-lying pore water hydrochemistry at Taiyuan, a series of groundwater mixing models were set up using PHREEQC. The results indicate that the deep pore water of the northern basin is recharged by karst water from the northern mountain area, shallow pore water from the northern basin, and Fenhe River water. Karst water from the northern mountain area is the most important source of recharge. Deep pore water in the western basin is recharged by karst water from the western mountain area and shallow pore water from the western basin. Hydrochemistry of deep pore water in the southern basin is affected by the mixing of deep pore water from the western basin with that from the northern basin. According to above analysis, the mixing process is the most important in forming the hydrochemistry of deep pore water at Taiyuan.
-
Key words:
- Taiyuan /
- deep-lying pore water /
- hydrochemical zonality /
- geochemical modeling /
- mixing process
-
分带性是自然界的普遍规律, 如植被分布、地貌形态、气候特征、水资源量及大气降水量变化等, 常在地域上呈现明显的分带现象(Guillen and Palanques, 1997; Sarkar and Guha, 1997; Mallik et al., 2001).在各类地质现象中, 分带性更为常见, 地质构造、岩浆活动乃至矿床、风化壳、地球化学元素等地质因素的分布规律时常体现出分带特点(Naito et al., 1995; Robion et al., 1997; Wogelius et al., 1997).地下水的水文地球化学特征作为地下水系统的重要属性, 其空间变化也往往表现出极强的分带性.太原市深层孔隙水的水化学特征受地形地貌、地层岩性、地质构造以及水文气象等多种因素的影响, 由边山向盆地中心呈现水平分带性, 便是一个极好的实例.
太原市为华北能源重工业城市, 地下水在工农业生产与城市生活供水方面起着举足轻重的作用.区内各类型地下水均有分布, 总体上山区以岩溶水与裂隙水为主, 盆地区以孔隙水为主.限于篇幅, 本文仅研究区内深层孔隙水(含水层埋藏深度在50~200 m之间) 的水文地球化学特征及其水化学分带现象, 并对深层孔隙水水化学组分的形成过程进行地球化学模拟.
1. 区域水文地质条件
太原市地处晋中断陷盆地北端, 北、东、西三面环山, 地形北高南低.盆地东部与西部分别属太行山系和吕梁山系, 北部则为棋子山, 标高1 300~1 780 m.汾河自兰村峡谷进入盆地后向南形成较为开阔的冲积平原, 标高771~810 m.整体地形从山区向盆地呈阶梯状下降, 但两侧不对称, 该地形特征对地下水径流具有明显的控制作用(图 1).
研究区内的前第四纪地层主要出露于盆地两侧及北部的山区地带.寒武系、奥陶系出露于山西汾河河谷, 并广泛分布于东西山的北部; 石炭系则分布于东西山区的中北部及西边山地带; 二叠系零星出露于东山与晋祠西北山区.
区内广泛分布第四纪地层, 成因类型俱全.下更新统(Q1) 自丘陵区至冲积平原大面积分布, 为一套河湖相杂色粘土亚粘土, 含砾石及泥灰岩地层; 中更新统(Q2) 分布也极广泛, 在丘陵区为黄土状土(Q2eol), 盆地区为冲积物(Q2al)和洪积物(Q2pl); 上更新统(Q3) 广泛分布于丘陵、平原区, 丘陵区以亚砂土为主, 山前地带为洪积黄土, 汾河两侧为冲积型砂砾卵石及亚砂土, 南部为砂、亚砂及亚粘土; 全新统(Q4) 分布于盆地边缘及河床地带, 由砂砾、砂、亚砂土及少量亚粘土组成地表岩性分布, 多具双层或多层结构.
松散岩类分布面积广、厚度大、孔隙发育, 形成了良好的天然储水构造, 含水层受外围山区岩溶裂隙水的侧向补给以及大气降水与汾河的渗漏补给, 水量非常丰富, 其主要排泄形式为人工开采及沿地下径流排泄于下游地区.按照含水层埋藏深度的差异, 可将区内的松散岩类孔隙水分为浅层孔隙水与深层孔隙水.其中浅层孔隙水埋深小于50 m, 而深层孔隙水埋深在50~200 m之间, 是本次研究的主要对象.
2. 深层孔隙水的水化学特征及其空间分带性
研究采用的深层孔隙水水样共25件, 均采集于2002年5月.采样点位置与水化学分析结果分别见图 1与表 1.由表 1可知, 各水样的pH值介于7.1~7.8之间, 为中性-略偏碱性水.地下水中的主要阴离子为HCO3-、SO42-, 主要阳离子为Ca2+、Mg2+, 部分水样的Cl-与Na+含量也很高.水样的TDS值变化非常大, 最小值(S09) 为252.2 mg/L, 最大值(S19) 达2 847.3 mg/L, 这表明研究区深层孔隙水的水化学特征差异较大.
表 1 深层孔隙水水化学特征由研究区深层孔隙水的TDS等值线图(图 2) 可知, 盆地北部孔隙水的TDS一般小于0.5 g/L, 向南逐渐变大, 在盆地西南部已大于1 g/L, 局部地区高达2.8 g/L.地下水的HCO3-、SO42-、Ca2+等离子的浓度(图略) 也表现出与TDS相似的空间分布特征, 在区内均由北向南渐渐升高.此外, 深层孔隙水中Cl-、Mg2+、Na+等离子浓度的变化趋势也与TDS的变化呈现出很强的一致性(图略).
边山的岩溶裂隙水在运移过程中侧向补给盆地孔隙水, 致使盆地内的地下水水化学类型在补给区水化学类型的基础上, 随着地形、地貌以及水文地质条件的变化, 呈现有规律的水平分带, 从山前到盆地依次有硫酸-重碳酸型、重碳酸-硫酸型、重碳酸型等水化学类型(图 3).
在盆地的北部山区及东西山北部, 碳酸盐岩大面积出露, 岩溶水水化学类型为HCO3-Ca·Mg型, 受岩溶水补给的影响, 在研究区北部的上兰镇、向阳镇、阳曲县一带, 深层孔隙水的水化学类型为重碳酸型(图 3).在汾河两侧, 由于汾河水库在放水期的大量渗透补给, 地下水水质被淡化, 致使重碳酸型水沿汾河两侧呈条带状展布(图 3).但沿汾河由南向北选取水化学类型为重碳酸型的地下水水样, 做各主要离子浓度变化曲线.结果显示: 地下水的TDS与HCO3-、SO42-、Cl-、Ca2+、Mg2+、Na+等离子含量大都呈现增加趋势(图 4).
东西山中部岩性变为碎屑岩类, 在煤系地层中含硫矿物的溶解作用下, 加以矿坑排水及其他人类活动的影响, 边山岩溶裂隙水的水化学类型为HCO3·SO4-Ca·Mg型与SO4·HCO3-Ca·Mg型.因而在盆地东西两侧的山前一带, 深层孔隙水受岩溶裂隙水补给, 水化学类型为硫酸-重碳酸型, 向盆地中心则逐渐过渡为重碳酸-硫酸型(图 3).
重碳酸-硫酸-氯化物型水只在晋阳湖与盆地东南侧的西温庄一带有小块分布(图 3), 这与该区是农业生产与生活区, 地下水水质受人类活动影响比较大有关.晋阳湖北侧还有一小块劣质水区域, 其水化学类型为硫酸-氯化物型(图 3).这种劣质水的形成除与该区地下水径流缓慢、蒸发浓缩作用较强有关外, 人类向晋阳湖中排污, 以致地下水水质恶化, 也是一个重要的因素.
3. 深层孔隙水水化学组成的地球化学模拟
前面的分析表明: 研究区深层孔隙水的水化学类型及其分带性主要受补给源地下水水质的控制.为深刻认识深层孔隙水水化学特征的形成过程, 笔者利用地球化学模拟软件PHREEQC来建立一系列地下水混合模型, 进而对不同区域的深层孔隙水的形成原因进行分析.在对比不同来源补给水混合结果与深层孔隙水的水化学特征时, 笔者采用硬石膏、文石、方解石、白云石、石膏5种矿物的饱和指数为对比指标, 这是由于在地下水由山区向盆地的径流及混合过程中, 碳酸盐矿物与硫酸盐矿物的溶解和沉淀是最重要的水文地球化学过程之一.建模所用的北部山区岩溶水、西部山区岩溶水、盆地北部浅层孔隙水、盆地西部浅层孔隙水及汾河水的水化学组成资料分别采用各区所采集水样的平均值(水样均采集于2002年5月), 详见表 2所示.
表 2 岩溶水、浅层孔隙水、汾河水的水化学组成对盆地北部的重碳酸型深层孔隙水(水样为S1—S11) 而言, 北部山区的岩溶水及汾河渗漏水是它的2个可能来源, 此外, 盆地北部的浅层孔隙水(含水层埋深小于50 m) 也会向下越流补给深层水.由此建立北部重碳酸型深层水的混合模型, 当岩溶水: 汾河水: 浅层孔隙水的比例为70%: 20%: 10%时, 最终混合溶液的各主要矿物(硬石膏、文石、方解石、白云石、石膏) 的饱和指数与深层孔隙水各主要矿物的平均饱和指数十分接近(表 3), 因此可以推断: 北部山区的岩溶水是盆地北部重碳酸型深层水的主要补给来源, 而汾河水的渗漏与浅层孔隙水的越流补给则是深层水水化学组分形成的次要原因.
表 3 混合溶液及深层孔隙水中硬石膏、文石、方解石、白云石、石膏的饱和指数Table Supplementary Table Saturation indexes of anhydrite, aragonite, calcite, dolomite and gypsum in the mixed solution and deep pore water盆地西部深层孔隙水的水化学类型为硫酸-重碳酸型、重碳酸-硫酸型、硫酸-氯化物型(代表性水样为S12、S15、S19、S23、S24), 推断西部边山岩溶水与盆地西部浅层孔隙水的混合作用为其水化学组分形成的主要原因.同样建立混合模型, 当岩溶水: 浅层孔隙水的比例为50%: 50%时, 最终混合溶液的硬石膏、文石、方解石、白云石、石膏等矿物的饱和指数与深层孔隙水中上述矿物的平均饱和指数最为接近(表 3), 可以认为西部边山岩溶水与盆地西部浅层孔隙水对该区域深层孔隙水化学组分的形成起到了同样重要的作用.
盆地南部深层孔隙水的水化学类型为重碳酸型(代表性水样为S16、S20、S25), 但与盆地北部的重碳酸型水相比, 其各主要离子含量均有不同程度的增加.结合盆地深层水的流线图(图 1) 来判断, 南部深层水应由盆地北部与西部的深层水混合而成, 混合模型表明, 当北部深层水: 西部深层水的混合比例为50%: 50%时, 最终混合溶液的硬石膏、文石、方解石、白云石、石膏等矿物的饱和指数与深层孔隙水中上述矿物的平均饱和指数最为接近(表 3).可以得出结论, 盆地南部是区内深层水水位最低的区域之一, 该处汇集了盆地北部与西部的深层水, 其水化学组分的形成过程最为复杂.
由于缺乏东部边山岩溶水与盆地东部浅层孔隙水的水化学资料, 故本文未分析盆地东部深层水的混合情况.但盆地北部、西部、南部深层水混合过程的地球化学模拟结果是非常令人满意的, 这表明混合作用是控制区域水化学状况的最重要的因素.
4. 结论
(1) 太原市深层孔隙水的水化学类型在补给区水化学类型的基础上, 随着地形、地貌、水文地质条件的变化, 呈现有规律的水平分带, 从山前到盆地依次分布有硫酸-重碳酸型水、重碳酸-硫酸型水、重碳酸型水等. (2) 混合作用是控制区域水化学状况的最重要的因素.盆地北部的深层孔隙水受到北部边山岩溶水、盆地北部浅层孔隙水、汾河水的补给, 其中北部边山岩溶水是最主要的补给源.而盆地西部的深层孔隙水由西部边山岩溶水与盆地西部浅层孔隙水混合而成.盆地南部的深层孔隙水则由盆地北部与西部的深层水混合而成.
-
表 1 深层孔隙水水化学特征
表 2 岩溶水、浅层孔隙水、汾河水的水化学组成
表 3 混合溶液及深层孔隙水中硬石膏、文石、方解石、白云石、石膏的饱和指数
Table 3. Saturation indexes of anhydrite, aragonite, calcite, dolomite and gypsum in the mixed solution and deep pore water
-
[1] Guillen, J., Palanques, A., 1997. A shoreface zonation in the Ebro Delta based on grain size distribution. Journal of Coastal Research, 13: 867 -878. https://www.jstor.org/stable/4298679 [2] Mallik, A.U., Lamb, E.G., Rasid, H., 2001. Vegetation zonation among the microhabitats in a lacustrine environment: Analysis and application of belowground species trait patterns. Ecological Engineering, 18: 135 -146. doi: 10.1016/S0925-8574(01)00069-6 [3] Naito, K., Fukahori, Y., Peiming, H., et al., 1995. Oxygen and carbon isotope zonations of wall rocks around the Kamioka Pb-Zn skarn deposits, central Japan: Application to prospecting. Journal of Geochemical Exploration, 54: 199 -211. doi: 10.1016/0375-6742(95)00044-5 [4] Robion, P., Kissel, C., Lamotte, D.F., et al., 1997. Magnetic mineralogy and metamorphic zonation in the Ardennes massif(France-Belgium). Tectonophysics, 271: 231-248. doi: 10.1016/S0040-1951(96)00268-5 [5] Sarkar, A., Guha, A.K., 1997. Pleistocene paleoclimatic zonation in northern Indian Ocean as revealed from Globorotalia Menardii abundance. Indian Journal of Marine Sciences, 26: 84 -87. http://www.infona.pl/resource/bwmeta1.element.elsevier-aca4fe88-7601-326f-9a90-9efb34b05a3b/tab/summary [6] Wogelius, R. A., Fraser, D.G., Wall, R. T., et al., 1997. Trace element and isotopic zonation in vein calcite from the Mendip Hills, UK, with spatial-process correlation analysis. Geochimica et Cosmochimica Acta, 61: 2037-2051. doi: 10.1016/S0016-7037(97)00065-3 -