-
摘要:
矿渣在沟谷的堆积形态、压占比例影响其在特定降雨概率下的稳定性。文中以豫西某金矿2条主沟16个渣堆为研究对象,提出在不同降雨概率雨力条件下单个渣堆稳定性计算方法。即近10年来最大观测雨强、25年一遇(P=4%)雨强、50年一遇雨强(P=2%)、100年一遇雨强(P=1%)条件下的启动流速(Uc)以及泄洪流速(Us),令稳定系数Fs=Uc/Us,计算出渣堆的稳定系数,研究成果对于不同特定雨力工况下渣堆的危险性评价及分类防治有一定意义。
Abstract:The slag accumulation form and the proportion of slag in the valley affect its stability under a specific rainfall probability. This paper takes 16 slag piles in two main trenches of a gold mine in western Henan as the research object and proposes a calculation method for the stability of a single slag pile under different rainfall probability conditions. That is, the start-up under the conditions of maximum observed rain intensity in recent years,rain intensity once in 25 years (P=4%), rain intensity once in 50 years (P=2%), and rain intensity once in 100 years (P=1%). The flow rate and the flow rate of the flood discharge, the ratio of which is calculated to analyze the stability of the slag pile. The research results have certain significance for the risk assessment of specific slag piles and the classification and prevention of slag piles.
-
Key words:
- debris flow /
- mine debris /
- flood discharge flow rate /
- startup flow rate /
- flood peak
-
海南岛地下水资源丰富,相对地表水,地下水具有不易受污染的优点。在地表水资源匮乏的季节或地区,地下水资源具有不可替代的补充作用。近年来,随着海南旅游业、房地产、农业和工业的快速发展,加之填海造田的影响,使地下水位下降、水质污染和海水入侵成为海南地下水资源安全最大最紧迫的威胁。因此,掌握海南岛的地下水化学成分信息,揭示其经历的水文地球化学过程,在保障海南岛城乡生活用水、农业用水和工业用水,支撑海南国家生态文明建设试验区和自由贸易岛、国际旅游岛,以及社会经济发展等方面发挥重要作用。
水文地球化学主要研究地下水化学成分的形成及其在地下水中的迁移演化过程[1-2]。地下水化学成分受控于当地地下水所处的围岩岩性、地表径流、植被、气候及人类活动状况,特别是海岸带地下水还可能受海水入侵影响。地下水在运动过程中,会与含水层介质中的化学成分和矿物组分发生一系列地球化学反应,这些过程决定了地下水化学成分的组成和演化过程[3-4]。因此,地下水中的化学组分能够反映影响水文地球化学演化的主要水-岩作用[5],揭示地下水的成因类型。
前人[6]研究了海南岛南部沿岸地下水水化学要素随大、小潮的波动和海水入侵特征,而对沿岸地下水的水文地球化学过程、水岩交互作用、成矿趋势及地下水的成因并未探讨。本文以海南南部沿海黎安港Z1、三亚湾Z2和板桥Z3三口监测井大、小潮时段的地下水为研究对象,利用Piper图、矿物相平衡分析、离子比例系数、Gibbs分析等,揭示海南南部沿海地下水的水文地球化学过程、水岩交互作用及成矿趋势和地下水的成因,为当地地下水的开发利用、污染防治等提供理论指导。
1. 研究区概况
本次研究的3口地下水监测井(Z1、Z2和Z3)分别位于海南南部的陵水县黎安港、三亚湾和东方市板桥镇的潮滩高潮线以上(表 1;图 1)。海南南部地下水资源按含水岩性及地下水动力特征可划为第四系松散岩类孔隙潜水、新近系松散岩类孔隙承压水、岩溶水及基岩裂隙水[7]。第四系松散岩类孔隙潜水分布于研究区山前斜地、山间谷地及滨海堆积平原。其中,滨海堆积平原一般平行于海岸,宽1~2 km,含水层厚度一般为5~15 m,含水层岩性主要为含贝壳中细砂、含砾亚砂土、中粗砂、砂砾石,富水性中等—丰富,地下水动态不稳定,随季节变化的趋势较明显,水位年动态变幅一般1~4 m,水位埋深多小于2 m。Z1监测井揭露的地层为第四系海陆交互相,Z2和Z3监测井揭露的地层为第四系海陆交互相和新近系上新统望楼港组。Z1、Z2和Z3监测的皆为潜水含水层,其中,Z1监测井含水层最厚,达12.35 m,岩性主要为砾砂、中砂、细砂等,Z2监测井含水层厚5.5 m,岩性主要为粗砂、细砂、粘土质砂等,Z3监测井含水层厚4.0 m,含水层岩性为砾砂、粗砂(图 2)。所在地层为第四系海陆交互相,其地下水为滨海平原第四系松散岩类孔隙潜水。孔隙潜水含水岩组岩性松散,透水性强,补给条件好,以大气降雨的垂直渗入补给为主。另外,埋藏于河流两侧的含水层与河水有密切联系,河水与潜水存在互补关系。研究区地下水总的流向是由山前和河谷上游向平原及沿海地区径流。除蒸发排泄外,主要是以水平径流方式排泄入海或地形低洼处及沙堤边缘排出地表,其次是在谷地上游排入裂隙溶洞含水层。滨海平原区孔隙潜水,矿化度一般为0.05~0.3 g/L,阴离子以Cl-和HCO3-为主,阳离子主要为Na+、Na+和Ca2+、Na+和Mg2+,水化学类型为Cl·HCO3-Na型、Cl·HCO3-Na·Ca型和Cl·HCO3-Na·Mg型。
表 1. 监测井地理位置和井深信息Table 1. Location and depth information of the monitoring wells监测井编号 东经 北纬 点位 井深/m Z1 110°03′59.51″ 18°24′39.87″ 陵水县黎安港 21.0 Z2 109°26′30.02″ 18°17′14.07″ 三亚市三亚湾 41.0 Z3 108°41′10.94″ 18°45′05.46″ 东方市板桥镇 42.0 2. 样品采集与分析
2014年6月在大潮、小潮时段分别对Z1、Z2和Z3监测井进行了一次连续26 h的地下水取样工作,地下水样每隔1 h取一次。Z1、Z2和Z3监测井大、小潮水样各26件,共156件。Z1监测井水样的取样深度为12.0~12.5 m,Z2监测井水样的取样深度为7.0~7.5 m,Z3监测井水样的取样深度为6.0~6.5 m。大潮时段Z1、Z2和Z3监测井的地下水静水位埋深分别为4.65 m、2.30 m和2.63 m;小潮时段Z1、Z2和Z3监测井的地下水静水位埋深分别为5.67 m、2.45 m和2.65 m。Z1监测井的地下水水位受潮汐影响较大,Z2和Z3监测井的地下水水位受潮汐影响较小。取样前取样瓶均采用蒸馏水清洗,再用所采水样润洗3次,水样用0.45 μm的滤膜过滤后冷藏。样品测试由海南省地质测试研究中心完成,测试内容包括Na+、K+、Mg2+、Ca2+、Cl-、SO42-、HCO3-、Br-、SiO2、TDS、总硬度、电导率和pH。Na+、K+、Mg2+、Ca2+、SO42-、SiO2使用电感耦合等离子体光谱仪IRIS Intrepid Ⅱ XSP,检测限分别为0.1 mg/L,0.1 mg/L,0.02 mg/L,0.05 mg/L,4 mg/L和0.1 mg/L。Cl-和HCO3-使用容量法,检测限分别为1.0 mg/L和1.5 mg/L,Br-使用比色法,其检测限为0.1 mg/L,总硬度和TDS分别使用容量法和重量法,pH使用酸度计,电导率采用电导率检测仪DDS-307A。所有测试依照地下水质检验方法DZ/T 0064—93标准执行。为方便讨论,本文以监测井名+st或+nt分别代表对应监测井大潮或小潮时段所取水样,如Z1st表示大潮时段Z1监测井所取水样。对所取的156个水样进行阴阳离子平衡检查,结果表明,大潮时段QZ1和QZ3各有一个水样的阴阳离子平衡误差(E)大于可接受误差5%,因此在数据分析中不考虑这2个水样。
3. 分析结果
3.1 地下水化学组分和参数特征
海南南部沿海3口监测井大、小潮时段的地下水化学组分和参数统计结果见表 2。Z1st和Z3st的主要阳离子平均浓度的顺序为Na+>Mg2+>K+>Ca2+,Z1nt和Z3nt的主要阳离子平均浓度顺序为Na+>Mg2+>Ca2+>K+,阴离子平均浓度顺序均为Cl->SO42->HCO3->Br-。这2口监测井中地下水的阳离子构成受大、小潮影响。Z2st和Z2nt的主要阳离子平均浓度顺序为Na+>Ca2+>K+>Mg2+,阴离子平均浓度顺序为HCO3->Cl->SO42-> Br-,此监测井中地下水的组分不受大、小潮影响。从图 3得知,在大、小潮时段内,所有监测井地下水中的化学组分和参数的灰色箱体四分位错开,但主体含量差异不大。
表 2. 大、小潮时段地下水化学参数统计Table 2. Chemical parameters of groundwater during spring tide and neap tide编号 样品数 统计项 化学成分/(mg·L-1) 化学参数 Na+ K+ Mg2+ Ca2+ Cl- SO42- HCO3- Br- 可溶性SiO2 TDS/(mg·L-1) 总硬度/(mg·L-1) 电导率/(mS·cm-1) pH Z1st 25 平均值 11412 430 1422 408 20124 2518 125.8 0.72 10.47 36444 6876 61.96 8.41 最小值 10700 406 1330 384 19100 2390 118 0.6 8.8 34400 6469 58.5 8.34 最大值 12200 458 1520 430 21200 2660 128 1.0 12.8 38600 7327 65.6 8.67 标准差 390.85 14.64 47.41 13.2 623.35 76.03 2.92 0.07 0.89 1153.65 219.89 1.96 0.06 变异系数 0.03 0.03 0.03 0.03 0.03 0.03 0.02 0.1 0.08 0.03 0.03 0.03 0.01 偏态系数 -0.1 -0.08 -0.06 0.08 0.03 -0.07 -0.94 3.12 0.72 -0.13 -0.03 -0.11 3.08 峰态系数 -0.6 -0.79 -0.45 -0.89 -0.96 -0.9 0 12.8 1.44 -0.80 -0.49 -0.84 13.36 Z1nt 26 平均值 11254 401 1298 413 18839 2594 125.1 0.4 3.69 34908 6357 59.42 8.36 最小值 10900 385 1230 390 18400 2400 123 0.4 2.8 34000 6050 57.9 8.32 最大值 11900 426 1380 623 19800 2780 128 0.4 6.8 36800 6756 62.6 8.38 标准差 264.17 10.5 34.33 43.59 386.86 88.82 2.52 0 0.94 710.45 165.44 1.20 0.02 变异系数 0.02 0.03 0.03 0.11 0.02 0.03 0.02 0 0.26 0.02 0.03 0.02 0 偏态系数 0.66 0.69 0.34 4.8 1.04 0.15 0.33 -1.06 1.96 1.13 0.46 1.14 -1.18 峰态系数 0.53 0.74 0.87 23.88 1.08 0.18 -2.06 -2.17 4.08 1.63 0.78 1.63 0.61 Z2st 26 平均值 80.32 10.69 8.63 58.38 96.78 37.88 227.2 0.34 19.19 426.85 181.81 0.726 7.92 最小值 68.7 9.8 8.3 49.8 83.3 35 197 L 17.9 404 160 0.688 7.88 最大值 129 11.2 9 61.5 160 49.9 242 0.7 22.3 527 190 0.897 8.00 标准差 15.41 0.36 0.14 2.93 20.17 3.74 10.63 0.22 1.07 32.61 7.46 0.05 0.03 变异系数 0.19 0.03 0.02 0.05 0.21 0.1 0.05 0.63 0.06 0.08 0.04 0.08 0 偏态系数 2.61 -0.92 0.12 -1.88 2.64 2.71 -1.87 1.02 1.98 2.55 -1.89 2.52 0.97 峰态系数 6.48 0.46 1.07 3.66 6.55 7.03 3.77 -0.25 4.22 6.15 3.73 6.03 0.34 Z2nt 26 平均值 78.4 10.89 7.85 56.85 94.97 41.78 226.4 0.33 19.46 423.65 174.31 0.722 7.78 最小值 67.8 10.3 7.3 42.2 80.5 39.3 202 L 18.9 385 136 0.656 7.62 最大值 119 11.4 8.4 63.6 142 49.3 237 0.4 21.6 501 193 0.854 8.13 标准差 10.59 0.27 0.22 5.42 12.71 1.97 7.26 0.11 0.55 21.84 13.73 0.04 0.09 变异系数 0.14 0.03 0.03 0.1 0.13 0.05 0.03 0.34 0.03 0.05 0.08 0.05 0.01 偏态系数 2.63 -0.84 0.33 -1.28 2.18 2.2 -1.31 -1.53 2.49 1.76 -1.31 1.78 2.24 峰态系数 8.42 0.66 1.66 1.33 6.81 7.95 4.1 1.26 8.46 5.80 1.66 5.87 9.56 Z3st 25 平均值 10588 451 1196 397 17820 2427 167.6 0.5 11.42 33040 5920 56.18 7.99 最小值 10200 432 1150 382 17000 2320 158 0.5 9.4 31900 5690 54.2 7.86 最大值 11600 492 1290 432 19100 2580 168 0.5 13.5 35600 6386 60.5 8.39 标准差 290.57 12.7 32.13 10.09 411.3 66.63 2 0 0.84 777.82 156.18 1.30 0.10 变异系数 0.03 0.03 0.03 0.03 0.02 0.03 0.01 0 0.07 0.02 0.03 0.02 0.01 偏态系数 1.9 1.49 1.12 1.88 1.08 0.58 -5 - 0.20 1.64 1.17 1.64 3.06 峰态系数 5.28 3.42 1.72 5.25 3.13 -0.07 25 - 1.20 3.81 2.07 3.97 13.18 Z3nt 26 平均值 10735 391 1266 416 18223 2445 157.6 0.29 5.22 33458 6249 56.97 8.10 最小值 10200 362 1210 398 17700 2370 153 L 4.7 31300 5975 53.3 7.92 最大值 11400 419 1360 438 19100 2540 163 0.3 5.7 35400 6701 60.3 8.23 标准差 308.47 12.8 39.82 11.29 417.91 53.46 1.96 0.06 0.22 952.12 191.76 1.65 0.09 变异系数 0.03 0.03 0.03 0.03 0.02 0.02 0.01 0.2 0.04 0.03 0.03 0.03 0.01 偏态系数 0.33 0.07 0.64 0.23 0.56 0.2 -0.79 -5.1 -0.26 0.09 0.51 0.08 -0.50 峰态系数 -0.34 0.09 -0.07 -0.62 -0.46 -1.09 4.18 26 0.42 0.18 -0.17 0.02 -0.29 标准海水[9] 10800 392 1290 410 19400 2700 142 67.3 注:“L”表示低于检测限;“-”表示不存在;“st”和“nt”分别代表大潮、小潮 变异系数(Cv)是样本标准差与均值的比值,能反映样本的离散程度,用来表征水化学参数随时空变化的稳定程度[8]。通常认为,变异系数0≤Cv≤10%时为弱变异,10%≤Cv≤100%为中等变异,Cv>100%为强变异。从表 2可知,3口监测井地下水化学参数的变异系数基本小于10%,为弱变异;Br-和可溶性SiO2偶尔为中等变异。化学参数在大、小潮时段较稳定。
4. 讨论
4.1 地下水分类
Robinove[10]根据TDS将自然水体划分为淡水、微咸水、中度咸水、极度咸水和盐水5类(表 3)。对大、小潮时段3口监测井地下水的TDS统计分析结果显示,Z1中地下水的TDS均大于35000 mg/L(图 4),属于盐水;Z3中的地下水小于35000 mg/L,大于10000 mg/L,为极度咸水;而Z2中的地下水落在淡水分界线以下,说明Z2井中水体为淡水。
表 3. 水质分类Table 3. Water quality classification分类 淡水 微咸水 中度咸水 极度咸水 盐水 TDS范围/(mg·L-1) 0~1000 1000~3000 3000~10000 10000~35000 >35000 4.2 地下水水化学类型
传统水文地球化学分析常用的手段有水化学类型分析法、离子和离子比值分析法。常用的水化学类型分析为Piper图示法。通过Piper图可以确定地下水的水化学类型,阴阳离子的变化规律,还可以估测地下水的来源及形成作用[11]。如Hamed[12]利用Piper三线图法确定Gafsa-Sidi Boubaker地区不同含水层地下水水化学类型;宋献方等[13]利用三线图分析潮白河流域地下水的形成作用及来源。根据阴阳离子的毫克当量百分数可以统计各监测井大小潮时段水样的水化学类型(表 4),并绘制出Piper三线图(图 5)。从表 4可以看出各个水样的水化学类型:大、小潮时段Z1和Z3监测井地下水水化学类型为Cl-Na型;Z2监测井地下水为HCO3·Cl-Na·Ca型。从Piper图(图 5)可以看到监测井的水文地球化学特征:3口监测井无论大潮还是小潮时段水样数据均落于2个区域。Z1和Z3监测井水样分布于海水占优势区域,表现为海水的水化学特征,说明这2个监测井区域已经完全被海水咸化;Z2监测井水样落在无优势型区域,即溶液中无明显占优势的阴阳离子,各离子含量相差不大,但是阳离子Na+稍高,占主导地位,偶有Ca2+显优势,阴离子HCO3-较高,偶尔Cl-稍高,阳离子Na+比阴离子Cl-的含量略小。由此推断,碳酸盐岩风化作用对该监测井区域地下水有主要影响,导致地下水中HCO3-占主导地位;由于经历阳离子交替吸附作用,使地下水中的Na+含量增加,Ca2+减少;岩盐的溶解或人为活动的影响,引起地下水中的Cl-溶度增加;3口监测井中的地下水与受大气降水影响的地表水体存在水力联系。
表 4. 大、小潮时段地下水水化学类型Table 4. Hydrochemical types of the groundwater during spring tide and neap tide样品编号 水化学类型 样品编号 水化学类型 Z1st Cl-Na Z1nt Cl-Na Z2st HCO3·Cl-Na·Ca Z2nt HCO3·Cl-Na·Ca Z3st Cl-Na Z3nt Cl-Na 4.3 地下水中主要离子间关系
Cl-具有很强的迁移能力,在低矿化水中较保守,难以形成难溶的矿物、也难以被生物积累,不被胶体吸附。监测井中地下水的6种离子K+、Na+、Ca2+、Mg2+、SO42-、HCO3-和Cl-(以meq/L计)的关系见图 6。各图横纵坐标均用对数坐标,斜线表示相应离子理论上的海水稀释线,即沿线的离子和氯离子的比值,又叫海水比值线[14],一般为标准海水的理论稀释线(TSDL),是用标准海水的Cl-浓度作为海水端元。假如某个具体的离子浓度在海水稀释线以下,则表示该离子被吸收或损耗,或由于岩盐或其他含氯化合物的溶解所致,相反则表示该离子富集[15]。
从图 6可知,Z2监测井所有样点的K+、Na+、Ca2+、SO42-、HCO3-值均位于海水比值线以上,且Ca2+、HCO3-偏离程度较大,而Mg2+值(图 6-c)基本落在海水比值线上或跳动,表明Z2监测井中的K+、Na+、Ca2+、SO42-、HCO3-出现富集。其主要原因可能是含钠矿物、含钾矿物、石膏等的风化溶解。如钠长石(Na2Al2Si6O16):
Na2Al2Si6O16+2CO2+3H2O→2HCO3−+2Na++H4Al2Si2O9+4SiO2 除Z3nt样点的Mg2+值(图 6-c)落在海水比值线以下外,Z1和Z3其余样点的K+、Na+、Ca2+、Mg2+、SO42-、HCO3-值均落在海水比值线上或其附近。这说明Z1和Z3监测井的地下水主要受海水影响,且Z3监测井小潮时段Mg2+被吸收或损耗,如通过阳离子交替吸附作用,Mg2+置换了Z3监测井围岩岩土颗粒所吸附的部分K+、Na+,使水体中的Mg2+减少,或者Mg2+与CO32-结合形成难溶于水的MgCO3。
4.4 地下水矿物相平衡分析
地下水矿物相的动态平衡是水体与围岩相互作用的结果。水体与围岩的相互作用是研究地球化学演化、地下水水化学特征、环境污染物质时空变化的基本理论[16]。饱和指数(Saturation Indices,SI)是水岩作用研究的基础,也是确定矿物在水体中是否达到饱和的重要依据[17]。Nath等[18]利用饱和指数建立的水文地球化学模型分析了地下水中砷的来源和迁移。饱和指数也被用于判定地下热水的一些物理化学性质[19-20]。建立饱和指数模型判定新安江库水浸泡坝址混凝土新、老龄期岩样不同时段的水溶液与固相介质之间的反应状态,并进一步评价混凝土岩样在浸泡过程中的溶蚀情况。
根据饱和指数SI可以判断和确定水与岩石、矿物之间处于何种反应状态,当SI大于0,说明矿物在水中已达到饱和状态,表明该矿物会从该溶液中沉淀出来;当SI小于0,该矿物为反应性矿物,在该溶液中能继续溶解直到达到平衡状态;当SI等于0,矿物溶解与沉淀达到平衡。一般视-0.5<SI<+0.5为平衡态,小于-0.5为溶解态,大于0.5为饱和态[2]。
对研究区大、小潮时段3口监测井Z1、Z2和Z3中的地下水样品运用水文地球化学软件Aquachem进行饱和指数模型分析,得到水体中硬石膏、文石、方解石、白云石、石膏、岩盐、石英、滑石8种矿物的饱和指数。
从表 5和图 7得知,大、小潮时段监测井Z1地下水中的文石、方解石、白云石、滑石的饱和指数大于0.5,说明Z1井中文石、方解石、白云石、滑石4种矿物在地下水中已经达到饱和状态,这4种矿物有发生沉淀的趋势,其中滑石的饱和指数达到8以上,很容易被析出,其次为白云石。石英的饱和指数均大于-0.5,小于0.5,说明其溶解和沉淀达到动态平衡。硬石膏、石膏和岩盐的饱和指数均小于-0.5,说明硬石膏、石膏、岩盐在Z1井溶解还未达到饱和,可继续在该水体中溶解。
表 5. 大、小潮时段地下水的矿物饱和指数统计结果Table 5. Statistics of mineral saturation index of groundwater during spring tide and neap tide编号 统计项 硬石膏 文石 方解石 白云石 石膏 岩盐 石英 滑石 Z1st 平均值 -0.89 1.01 1.16 3.23 -0.69 -2.47 0.27 8.93 最小值 -0.92 0.96 1.10 2.42 -0.73 -2.52 0.20 5.78 最大值 -0.77 1.13 1.27 3.51 -0.57 -2.42 0.36 10.71 标准差 0.03 0.04 0.04 0.18 0.03 0.03 0.03 0.77 反应状态 溶解态 饱和态 饱和态 饱和态 溶解态 溶解态 平衡态 饱和态 Z1nt 平均值 -0.86 0.98 1.13 3.17 -0.67 -2.50 -0.19 6.75 最小值 -0.91 0.93 1.07 3.07 -0.72 -2.53 -0.30 6.32 最大值 -0.67 1.15 1.29 3.33 -0.47 -2.46 0.08 7.83 标准差 0.04 0.04 0.04 0.05 0.04 0.02 0.10 0.40 反应状态 溶解态 饱和态 饱和态 饱和态 溶解态 溶解态 平衡态 饱和态 Z2st 平均值 -2.31 0.45 0.59 0.74 -2.10 -6.71 0.43 1.54 最小值 -2.34 0.31 0.45 0.52 -2.13 -6.82 0.40 1.26 最大值 -2.27 0.55 0.69 0.91 -2.06 -6.28 0.50 2.11 标准差 0.02 0.05 0.05 0.08 0.02 0.14 0.02 0.21 反应状态 溶解态 平衡态 饱和态 饱和态 溶解态 溶解态 平衡态 饱和态 Z2nt 平均值 -2.27 0.30 0.45 0.42 -2.07 -6.72 0.43 0.66 最小值 -2.40 0.16 0.30 0.10 -2.19 -6.84 0.42 -0.31 最大值 -2.23 0.54 0.68 0.98 -2.03 -6.37 0.47 2.64 标准差 0.04 0.07 0.07 0.15 0.04 0.10 0.01 0.53 反应状态 溶解态 平衡态 平衡态 平衡态 溶解态 溶解态 平衡态 饱和态 Z3st 平均值 -0.89 0.52 0.66 2.21 -0.69 -2.55 0.32 6.48 最小值 -0.91 0.40 0.54 1.97 -0.71 -2.58 0.23 5.86 最大值 -0.85 1.10 1.24 3.37 -0.65 -2.48 0.39 8.83 标准差 0.01 0.13 0.13 0.26 0.01 0.02 0.03 0.55 反应状态 溶解态 饱和态 饱和态 饱和态 溶解态 溶解态 平衡态 饱和态 Z3nt 平均值 -0.87 0.59 0.73 2.35 -0.68 -2.54 -0.03 5.83 最小值 -0.89 0.45 0.59 2.08 -0.69 -2.57 -0.08 4.88 最大值 -0.85 0.69 0.83 2.56 -0.65 -2.49 0.02 6.55 标准差 0.01 0.07 0.07 0.13 0.01 0.02 0.02 0.45 反应状态 溶解态 饱和态 饱和态 饱和态 溶解态 溶解态 平衡态 饱和态 Z2井地下水中矿物的溶解与沉淀趋势如表 5和图 7所示。大、小潮时段表现有差异,其中大潮时段,滑石、白云石、方解石的饱和指数大于0.5,说明大潮时段Z2井中滑石、白云石、方解石3种矿物在地下水中已经达到饱和状态,这3种矿物有发生沉淀的趋势,其中滑石容易析出;石英和文石的饱和指数介于-0.5~0.5之间,其溶解和沉淀达到动态平衡。硬石膏、石膏、岩盐的饱和指数小于-0.5,说明硬石膏、石膏、岩盐在Z2井溶解还未达到饱和度,可继续在该水体中溶解,其中岩盐饱和指数小于-6,最容易溶解。小潮时段滑石的饱和指数多为大于0.5,偶有小于0.5,但大于-0.5,说明小潮时段滑石多表现为沉淀趋势,偶达平衡态。文石、方解石、白云石和石英的饱和指数介于-0.5~0.5之间,说明这4种矿物相在小潮时段在地下水中的溶解与沉淀达到动态平衡。硬石膏、石膏、岩盐的饱和指数与大潮时段表现一致,小于-0.5,说明硬石膏、石膏、岩盐在Z2井溶解还未达到饱和,可继续在该水体中溶解,其中岩盐饱和指数小于-6,最容易溶解。
大、小潮时段Z3井地下水中矿物的溶解与沉淀趋势相同,具体如表 5和图 7所示。滑石、白云石、方解石和文石的饱和指数大于0.5,说明Z3井中滑石、白云石、方解石和文石4种矿物在地下水中已经达到饱和状态,这4种矿物有发生沉淀的趋势,其中滑石最容易析出。石英的饱和指数介于-0.5~0.5之间,其溶解和沉淀达到动态平衡。硬石膏、石膏、岩盐的饱和指数小于-0.5,说明硬石膏、石膏、岩盐在Z3井溶解还未达到饱和,硫酸盐和岩盐的溶解是正在发生并持续的过程,其中岩盐最易溶解。
Z1、Z2和Z3无论在大潮时段还是小潮时段,其地下水的滑石均处于饱和态,石英处于平衡态,硬石膏、石膏和岩盐均处于溶解态。Z1、Z3监测井中的文石、方解石和白云石在大潮和小潮时段均处于饱和态;Z2监测井中的文石处于平衡态,但是方解石和白云石在大潮时段处于饱和态,而小潮时段处于平衡态(表 5)。结合以上分析,推断硬石膏、石膏和岩盐溶解是海南南部沿海地下水中的主要离子来源,且大、小潮影响Z2中的地下水所含钙镁碳酸盐的水文地球化学行为。
4.5 地下水成因
4.5.1 γNa/γCl分析
不同条件或成因环境下形成的地下水,部分离子在浓度上的比值有明显差异,可以利用这些比值判定地下水化学组分的成因类型。如Vengosh等[21]利用Na/Cl值和Br/Cl值判断土耳其西部地下热水的来源。γNa/γCl值被用来分析地下水的成因及经历的水文地球化学作用[22]。
标准海水的γNa/γCl因子平均值为0.85。海相沉积水在地质历史过程中,如果水中的Na+与地层中的交换性Ca2+产生阳离子交换或经浓缩至析出石盐阶段,则Na+含量会下降,γNa/γCl值将小于0.85[23-24];若是含岩盐地层溶滤形成的地下水,其地下水成因系数γNa/γCl接近于1;若γNa/γCl大于1时,则受含钠硅酸盐(如前文提到的钠长石)风化溶解影响[25-26]或人为活动排放的生活污水所致[27-28]。
图 8显示,大、小潮时段3口监测井中地下水的γNa/γCl值Z2大于Z3和Z1,低矿化度水具有较高的γNa/γCl系数,高矿化度水具有较低的γNa/γCl系数[29-30]。但是3口井均大于标准海水的0.85,因此地下水中的Na+不与地层中的交换性Ca2+产生阳离子交换;Z1和Z3监测井中地下水的γNa/γCl值接近于0.85,其来源与海水紧密相关;Z2的γNa/γCl值大于1,说明地下水中存在过量的Na,该监测井远离化粪池,因此,该监测井中的地下水更可能受含钠硅酸盐风化溶解影响。同时,3口井中地下水的γNa/γCl值与1相差不大,说明地下水的化学组分与岩盐溶滤作用有关,这与前文的饱和指数推断结果一致。
4.5.2 吉布斯(Gibbs)分析
通过对世界大量雨水、河水、湖水及海水化学组分的分析,Gibbs将天然水体组分的控制因素划分为3种类型:降雨控制型、岩石风化控制型和蒸发-浓缩控制型[31]。并利用矿化度(TDS)与阳离子Na+/(Na++Ca2+)和阴离子Cl-/(Cl-+HCO3-)之间的内在关系图,即吉布斯(Gibbs)图划分出3种控制因素下天然水体对应的特征区域。Gibbs图的对数纵坐标表示TDS的质量浓度,横坐标为普通坐标,表示阳离子Na+/(Na++Ca2+)或阴离子Cl-/(Cl-+HCO3-)的质量浓度比。其后被广泛应用于评价形成地下水水化学成分的水文地球化学作用,分析水化学离子的成因类型[32]。如王亚平等[33]通过吉布斯(Gibbs)图分析了长江水系水文地球化学特征及成因。
从图 9可以看出,无论在大潮时段还是在小潮时段,Z2监测井中地下水的Cl-/(Cl-+HCO3-)值介于0.30~0.5之间,Na+/(Na++Ca2+)值介于0.5~0.58之间,水化学环境主要受岩石矿物风化溶解作用影响;Z1和Z3监测井中地下水的Cl-/(Cl-+HCO3-)和Na+/(Na++Ca2+)值均大于0.9,水样落在蒸发-浓缩控制型区域,其主要受高矿化度的海水影响。
5. 结论
(1) 依据Robinove水质分类标准,监测井Z1、Z2和Z3中地下水分别为盐水、淡水和极度咸水。根据舒卡列夫分类和Piper三线图分析结果,Z1、Z3中地下水的水化学类型为Cl-Na型,Z2中的为HCO3·Cl-Na·Ca型。
(2) 在大、小潮时段,Z1、Z2和Z3中的地下水的滑石均处于饱和态,石英处于平衡态,硬石膏、石膏和岩盐均处于溶解态。Z1、Z3中的文石、方解石和白云石均处于饱和态;Z2中的文石处于平衡态,但方解石和白云石在大潮时段处于饱和态,小潮时段处于平衡态。结合以上分析,推断硬石膏、石膏和岩盐溶解是海南南部沿海地下水中的主要离子来源。大、小潮影响Z2的地下水中钙镁碳酸盐的水文地球化学行为。
(3)γNa/γCl和Gibbs水化学成因分析表明,Z1和Z3的地下水来源与受蒸发浓缩控制型的高矿化度的海水紧密相关,Z2的水化学环境主要受含钠硅酸盐矿物的岩石风化溶解作用影响。3口井中地下水的化学组分还与岩盐溶滤作用有关。
致谢
衷心感谢海南省地质调查研究院对野外数据采集提供的帮助,海南省地质测试研究中心提供良好的实验条件使文稿获得了精确的实验数据,及审稿专家在文稿修改过程中提出的宝贵意见。
-
表 1 豫西某金矿渣堆体积及压占沟谷比例统计表
Table 1. Statistical table of volume and proportion of the slag in a gold slag pile in western Henan
矿渣 ZD1 ZD2 ZD3 ZD4 ZD5 ZD6 ZD7 ZD8 体积/(104 m3) 1.2 0.09 0.12 0.07 0.14 0.5 0.05 1.4 压占沟谷比例/% 50 43 83 69 72 71 74 48 矿渣 ZD9 ZD10 ZD11 ZD12 ZD13 ZD14 ZD15 ZD16 体积/(104 m3) 2 2.5 1.2 0.29 0.4 0.3 0.25 1.54 压占沟谷比例/% 74 55 77 71 60 67 49 88 表 2 渣堆泄洪流速Us计算参数测量结果表
Table 2. The measurement result of calculation parameters of flood discharge velocity of the slag pile
渣堆编号 ZD1 ZD2 ZD3 ZD4 ZD5 ZD6 ZD7 ZD8 F/km2 0.176 0.176 0.27 0.27 0.469 0.469 0.616 0.176 L/km 0.47 0.47 0.61 0.61 0.72 0.79 0.79 1.01 I/‰ 462 462 418 418 387 373 373 314 h/m 2 2 2.5 2.5 2.5 2 2 2 b/m 1.5 3.5 2.5 3.35 3.05 11 11.5 1.5 渣堆编号 ZD9 ZD10 ZD11 ZD12 ZD13 ZD14 ZD15 ZD16 F/km2 0.85 1.68 1.82 1.82 0.232 0.283 0.431 0.511 L/km 1.52 1.52 1.66 0.52 0.66 0.84 0.97 1.13 I/‰ 283 283 269 514 456 400 374 332 h/m 2.5 2 2 2 2.5 2.5 2 2 b/m 2.7 9.5 26.5 18.5 4.35 2 4.5 4 表 3 渣堆启动流速(Uc)计算参数测量结果表
Table 3. The measurement result of calculation parameters of startup flow rate of the slag pile
渣堆编号 ZD1 ZD2 ZD3 ZD4 ZD5 ZD6 ZD7 ZD8 α/(°) 14 14 9 10.84 14.16 10.45 6.76 6 φ/(°) 34 32 20 22 41 23 18 33 渣堆编号 ZD9 ZD10 ZD11 ZD12 ZD13 ZD14 ZD15 ZD16 α/(°) 7 4 4 12.84 12.95 8.38 11.14 6 φ/(°) 30 34 29 30 33 24 17 29 表 4 不同雨力工况下渣堆稳定性计算结果表
Table 4. The calculation result of slag pile stability under different rain conditions
近10 a最大值计算
结果(49 mm/h)渣堆编号 ZD1 ZD2 ZD3 ZD4 ZD5 ZD6 ZD7 ZD8 Fs 0.90 2.11 1.40 1.85 0.90 2.76 2.26 0.51 渣堆编号 ZD9 ZD10 ZD11 ZD12 ZD13 ZD14 ZD15 ZD16 Fs 0.72 1.88 1.33 2.64 1.01 1.25 0.95 0.68 25 a一遇计算
结果(60 mm/h)渣堆编号 ZD1 ZD2 ZD3 ZD4 ZD5 ZD6 ZD7 ZD8 Fs 0.72 1.70 1.12 1.49 0.72 2.22 1.81 0.41 渣堆编号 ZD9 ZD10 ZD11 ZD12 ZD13 ZD14 ZD15 ZD16 Fs 0.58 1.51 1.07 2.12 0.81 1.00 0.76 0.54 50 a一遇计算
结果(80 mm/h)渣堆编号 ZD1 ZD2 ZD3 ZD4 ZD5 ZD6 ZD7 ZD8 Fs 0.53 1.25 0.83 1.09 0.53 1.63 1.33 0.30 渣堆编号 ZD9 ZD10 ZD11 ZD12 ZD13 ZD14 ZD15 ZD16 Fs 0.42 1.11 0.79 1.56 0.60 0.74 0.56 0.40 100 a一遇计算
结果(90 mm/h)渣堆编号 ZD1 ZD2 ZD3 ZD4 ZD5 ZD6 ZD7 ZD8 Fs 0.47 1.10 0.73 0.96 0.47 1.44 1.18 0.26 渣堆编号 ZD9 ZD10 ZD11 ZD12 ZD13 ZD14 ZD15 ZD16 Fs 0.37 0.98 0.69 1.38 0.53 0.65 0.49 0.35 表 5 不同雨力下渣堆危险性以及稳定性系数
Table 5. Ranking table of slag pile stability under different rain conditions
渣堆
编号ZD
8ZD
16ZD
9ZD
1ZD
5ZD
15ZD
13ZD
14ZD
11ZD
3ZD
4ZD
10ZD
2ZD
7ZD
12ZD
649 mm/h Fs 0.51 0.68 0.72 0.90 0.90 0.95 1.01 1.25 1.33 1.40 1.85 1.88 2.11 2.26 2.64 2.76 60 mm/h Fs 0.41 0.54 0.58 0.72 0.72 0.76 0.81 1.00 1.07 1.12 1.49 1.51 1.70 1.81 2.12 2.22 80 mm/h Fs 0.30 0.40 0.42 0.53 0.53 0.56 0.60 0.74 0.79 0.83 1.09 1.11 1.25 1.33 1.56 1.63 90 mm/h Fs 0.26 0.35 0.37 0.47 0.47 0.49 0.53 0.65 0.69 0.73 0.96 0.98 1.10 1.18 1.38 1.44 表 6 不同雨力工况下渣堆稳定性修正计算结果表
Table 6. The calculation result of slag pile stability correction under different rain conditions
近10年最大观测雨强
修正结果(49 mm/h)渣堆编号 ZD1 ZD2 ZD3 ZD4 ZD5 ZD6 ZD7 ZD8 Fs 0.90 1.41 1.01 1.33 0.70 2.26 1.80 0.43 渣堆编号 ZD9 ZD10 ZD11 ZD12 ZD13 ZD14 ZD15 ZD16 Fs 0.62 1.48 1.12 2.64 1.01 1.25 0.95 0.62 25年一遇雨强修正
结果(60 mm/h)渣堆编号 ZD1 ZD2 ZD3 ZD4 ZD5 ZD6 ZD7 ZD8 Fs 0.72 1.13 0.81 1.06 0.56 1.70 1.44 0.34 渣堆编号 ZD9 ZD10 ZD11 ZD12 ZD13 ZD14 ZD15 ZD16 Fs 0.47 1.22 0.87 2.12 0.81 0.84 0.63 0.45 50年一遇雨强修正
结果(80 mm/h)渣堆编号 ZD1 ZD2 ZD3 ZD4 ZD5 ZD6 ZD7 ZD8 Fs 0.53 0.83 0.59 0.78 0.41 1.24 1.05 0.25 渣堆编号 ZD9 ZD10 ZD11 ZD12 ZD13 ZD14 ZD15 ZD16 Fs 0.35 0.88 0.62 1.56 0.60 0.62 0.46 0.33 100年一遇雨强修正
结果(90 mm/h)渣堆编号 ZD1 ZD2 ZD3 ZD4 ZD5 ZD6 ZD7 ZD8 Fs 0.47 0.73 0.52 0.68 0.36 1.10 0.93 0.22 渣堆编号 ZD9 ZD10 ZD11 ZD12 ZD13 ZD14 ZD15 ZD16 Fs 0.31 0.78 0.55 1.38 0.53 0.54 0.41 0.29 表 7 修正后不同雨力下渣堆危险性以及稳定性系数
Table 7. Ranking table of slag pile stability under different rain conditions after correction
49 mm/h 渣堆编号 ZD
8ZD
9ZD
16ZD
5ZD
1ZD
15ZD
3ZD
13ZD
11ZD
14ZD
4ZD
2ZD
10ZD
7ZD
6ZD
12Fs 0.43 0.62 0.62 0.70 0.90 0.95 1.01 1.01 1.12 1.25 1.33 1.41 1.48 1.80 2.26 2.64 60 mm/h 渣堆编号 ZD
8ZD
16ZD
9ZD
5ZD
15ZD
1ZD
3ZD
13ZD
14ZD
11ZD
4ZD
2ZD
10ZD
7ZD
6ZD
12Fs 0.34 0.45 0.47 0.56 0.63 0.72 0.81 0.81 0.84 0.87 1.06 1.13 1.22 1.44 1.70 2.12 80 mm/h 渣堆编号 ZD
8ZD
16ZD
9ZD
5ZD
15ZD
1ZD
3ZD
13ZD
11ZD
14ZD
4ZD
2ZD
10ZD
7ZD
6ZD
12Fs 0.25 0.33 0.35 0.41 0.46 0.53 0.59 0.60 0.62 0.62 0.78 0.83 0.88 1.05 1.24 1.56 90 mm/h 渣堆编号 ZD
8ZD
16ZD
9ZD
5ZD
15ZD
1ZD
3ZD
13ZD
14ZD
11ZD
4ZD
2ZD
10ZD
7ZD
6ZD
12Fs 0.22 0.29 0.31 0.36 0.41 0.47 0.52 0.53 0.54 0.55 0.68 0.73 0.78 0.93 1.10 1.38 -
[1] 徐友宁,陈华清,杨敏,等. 采矿废渣颗粒粒径对矿渣型泥石流起动的控制作用—以小秦岭金矿区为例[J]. 地质通报,2015,34(11):1993 − 2000. [XU Youning,CHEN Huaqing,YANG Min,et al. Controlling role of particle sizes of mining waste residues in the initiation of mine debris flow:A case study of the Xiaoqinling gold mining area[J]. Geological Bulletin of China,2015,34(11):1993 − 2000. (in Chinese with English abstract)
[2] 徐友宁,何芳,张江华,等. 矿山泥石流特点及其防灾减灾对策[J]. 山地学报,2010,28(4):463 − 469. [XU Youning,HE Fang,ZHANG Jianghua,et al. Characteristics of mine debris flow and its disaster prevention and mitigation countermeasures[J]. Journal of Mountain Science,2010,28(4):463 − 469. (in Chinese with English abstract)
[3] 徐友宁,何芳,陈华清. 西北地区矿山泥石流及分布特征[J]. 山地学报,2007,25(6):729 − 736. [XU Youning,HE Fang,CHEN Huaqing. Mine debris flow and its distribution in Northwestern China[J]. Journal of Mountain Science,2007,25(6):729 − 736. (in Chinese with English abstract)
[4] 陈华清. 小秦岭金矿区“7·23”泥石流形成特征及其启示[C]//第八届海峡两岸山地灾害与环境保育学术研讨会论文集. 中国水土保持学会、台湾中华水土保持学会、中国水土保持学会泥石流滑坡防治专业委员会, 2011.
CHEN Huaqing. Formation characteristics and enlightenment of "7·23" debris flow in Xiaoqinling gold mining area[C]//Proceedings of the 8th Cross-Strait Mountain Disaster and Environmental Conservation Symposium. Chinese Society of Soil and Water Conservation, Chinese Society of Soil and Water Conservation in Taiwan, Professional Committee of Debris Flow and Landslide Prevention of Chinese Society of Soil and Water Conservation, 2011. (in Chinese with English abstract)
[5] 徐友宁,陈华清,张江华,等. 小秦岭金矿区 7·23蒿岔峪泥石流成灾模式及启示[J]. 地质通报,2015,34(11):2001 − 2008. [XU Youning,CHEN Huaqing,ZHANG Jianghua,et al. Hazard mode of the 7·23 mine debris flow in the Xiaoqinling gold mine area[J]. Geological Bulletin of China,2015,34(11):2001 − 2008. (in Chinese with English abstract) doi: 10.3969/j.issn.1671-2552.2015.11.002
[6] 徐友宁. 影响小秦岭金矿区矿山泥石流形成的物源特征分析[C]// 第八届海峡两岸山地灾害与环境保育学术研讨会论文集. 中国水土保持学会、台湾中华水土保持学会、中国水土保持学会泥石流滑坡防治专业委员会, 2011
XU Youning. Analysis of Provenance Characteristics Affecting the Formation of Mine Debris Flow in Xiaoqinling Gold Mine Area[C]//Proceedings of the 8th Cross-Strait Mountain Disaster and Environmental Conservation Symposium. Chinese Society of Soil and Water Conservation, Chinese Society of Soil and Water Conservation in Taiwan, Professional Committee of Debris Flow and Landslide Prevention of Chinese Society of Soil and Water Conservation, 2011. (in Chinese with English abstract)
[7] 何芳,徐友宁,乔冈,等. 中国矿山地质灾害分布特征[J]. 地质通报,2012,31(增刊 1):476 − 485. [HE Fang,XU Youning,QIAO Gang,et al. Distribution characteristics of mine geological hazards in China[J]. Geological Bulletin of China,2012,31(Sup 1):476 − 485. (in Chinese with English abstract)
[8] 邓龙胜,范文,熊炜,等. 矿渣型泥石流发育特征及危险性评价[J]. 工程地质学报,2009,17(3):415 − 420. [DENG Longsheng,FAN Wen,XIONG Wei,et al. Development features and risk of inducing slag debris flow at Daxicha gully[J]. Journal of Engineering Geology,2009,17(3):415 − 420. (in Chinese with English abstract) doi: 10.3969/j.issn.1004-9665.2009.03.022
[9] 李荣,李义天,王迎春. 非均匀沙起动规律研究[J]. 泥沙研究,1999(1):27 − 32. [LI Rong,LI Yitian,WANG Yingchun. Study of non-uniform sand starting law[J]. Journal of Sediment research,1999(1):27 − 32. (in Chinese with English abstract) doi: 10.16239/j.cnki.0468-155x.1999.01.005
[10] 陈嫒儿, 谢鉴衡. 非均匀沙起动规律初探[J]. 武汉水利电力学院学报,1988,21(3):28 − 37. [CHEN Yuaner, XIE Jianheng. Preliminary research on the laws of threshold motion for non-uniform sediments[J]. Journal of Wuhan University of Hydraulic and Electric Engineering,1988,21(3):28 − 37. (in Chinese with English abstract)
[11] 谢鉴衡. 河流泥沙动力学[M]. 北京: 中国水利水电出版社, 1989: 63 − 84
XIE Jianheng. River sediment dynamics[M]. Beijing: China Water Conservancy and Hydropower Press, 1989: 63 − 84. (in Chinese)
[12] 秦荣昱. 不均匀砂的启动规律[J]. 泥沙研究, 1980(复刊号): 83 − 91
QIN Rongyu. Startup rule of uneven sand[J]. Sediment Research, 1980(Reissue number): 83 − 91. (in Chinese)
[13] 彭润译, 吕秀贞. 长江寸滩站卵石推移质输沙规律[J]. 水利学报, 1990, 21(1): 38 − 43
PENG Runyi, LYU Xiuzhen. The law of bedding and sediment transport of pebble at Cuntan station of the Yangtze River[J]. Journal of Hydraulic Engineering. 1990, 21(1): 38 − 43.(in Chinese)
[14] 林玫玲,简文彬,胡海瑞,等. 基于离散元的矿渣泥石流运动过程研究[J]. 中国地质灾害与防治学报,2017,28(2):10 − 14. [LIN Meiling,JIAN Wenbin,HU Hairui,et al. The movement process of slag debris flow based on the discrete element method[J]. The Chinese Journal of Geological Hazard and Control,2017,28(2):10 − 14. (in Chinese with English abstract) doi: 10.16031/j.cnki.issn.1003-8035.2017.02.02
[15] 李建林,张洪云,李振林. 矿山泥石流沟谷形态的分形分维[J]. 中国地质灾害与防治学报,2012,23(1):1 − 5. [LI Jianlin,ZHANG Hongyun,LI Zhenlin. Fractal analysis on morphology of mine debris flow gully[J]. The Chinese Journal of Geological Hazard and Control,2012,23(1):1 − 5. (in Chinese with English abstract) doi: 10.3969/j.issn.1003-8035.2012.01.001
[16] 唐亚明,武立,冯凡,等. 泥石流风险减缓措施及经济决策—以山西吉县城北沟为例[J]. 西北地质,2021,54(4):227 − 238. [TANG Yaming,WU Li,FENG Fan,et al. Risk mitigation measures and economic decisions on debris flow:Taking Beigou of Jixian County,Shanxi Province as an example[J]. Northwestern Geology,2021,54(4):227 − 238. (in Chinese with English abstract) doi: 10.19751/j.cnki.61-1149/p.2021.04.018
[17] 杨敏, 徐友宁. 小秦岭金矿区矿渣泥石流成因机理及防治对策[M]. 北京: 冶金工业出版社, 2021
YANG Min, XU Youning. Genesis mechanism and prevention countermeasures of slag and debris flow in Xiaoqinling gold mine area[M]. Beijing: Metallurgical Industry Press, 2021. (in Chinese)
[18] 徐友宁,陈社斌,李育敬,等. 陕西潼关金矿区泥石流潜势度评价[J]. 水文地质工程地质,2006,33(2):89 − 92. [XU Youning,CHEN Shebin,LI Yujing,et al. Potentiality degree assessment of slag mudslide geo-hazard on the gold mine area in Tongguan, Shaanxi Province[J]. Hydrogeology & Engineering Geology,2006,33(2):89 − 92. (in Chinese with English abstract) doi: 10.3969/j.issn.1000-3665.2006.02.021
[19] 高东光. 桥涵水文[M]. 北京: 人民交通出版社, 2005: 76 − 87
GAO Dongguang. Hydrology and hydraulics for bridge engineering[M]. Beijing: China Communications Press, 2005: 76 − 87. (in Chinese)
[20] 燕建设,王铭生,杨建朝,等. 豫西马超营断裂带的构造演化及其与金等成矿的关系[J]. 中国区域地质,2000,19(2):166 − 171. [YAN Jianshe,WANG Mingsheng,YANG Jianchao,et al. Tectonic evolution of the Machaoying fault zone in western Henan and its relationship with Au polymetallic mineralization[J]. The Regional Geology of China,2000,19(2):166 − 171. (in Chinese with English abstract)
[21] 中国科学院水利部成都山地灾害与环境研究所. 中国泥石流[M]. 北京: 商务印书馆, 2000
Institute of Mountain Hazards and Environment, CAS. China debris flow[M]. Beijing: The Ommercial Press, 2000. (in Chinese)
[22] 常士骠, 张苏民. 工程地质手册(第五版)[M]. 北京: 中国建筑工业出版社, 2018: 691
CHANG Shibiao, ZHANG Sumin. Geological engineering handbook(The fifth edition)[M]. Beijing: China Architecture & Building Press, 2018: 691. (in Chinese)
-