Analysis of electrical resistivity and acoustic wave velocity characteristics in fault structures and their combined application in the detection process
-
摘要:
查明断裂构造的具体位置对于工程场址的选择和稳定性评价具有重要意义。高密度电阻率法和声波测井是广泛运用于探测断裂构造的地球物理方法。高密度电阻率法具有对低阻体敏感且探测范围较大的特点，但高密度电阻率法的分辨率有限，难以准确圈定地质体的边界位置；声波测井具有准确识别地层的特点，但其横向探测范围非常有限。为了克服上述问题，结合2种方法在断裂构造探测中各自的优势以提高圈定地质体的精度，文章首先通过物性试验证明岩石破裂前后电阻率与声波波速变化的相关性，随后用一实例分析断裂构造的电阻率与声波波速特征以及两者在探测过程中的联合运用，最后借助声波测井数据对高密度电阻率法反演结果进行深度校正，得到断裂构造在地下空间的三维分布位置。结果表明：岩石破裂后的电阻率与声波波速较未破裂前均偏低；断裂构造的破碎带在高密度电阻率法反演剖面上以低阻异常区域存在，在声波波速曲线上破碎带区域的波速值会有明显的骤降，2种方法的组合探测模式可为断裂构造具体位置的圈定提供更好的指导。钻探结果验证了借助声波测井数据校正后的高密度电阻率法反演结果比未校正之前更加符合实际地质情况。
Abstract:Finding out the specific location of the fault structure is of great significance for the engineering site selection and stability evaluation. Electrical resistivity tomography and acoustic logging are geophysical methods widely used in detecting fracture structures. Although electrical resistivity tomography is sensitive to low resistance bodies and has a large detection range, it is difficult to accurately delineate the boundary position of geological bodies because of its limited resolution. Acoustic logging can identify strata accurately, but its lateral detection range is limited. This study combined the advantages of the two methods in fault structure detection to improve the precision of geological body delineation. Firstly, the correlation between the resistivity of rock before and after fracture and the change of acoustic wave velocity was proved by physical property experiments. Then the case analysis of the resistivity and acoustic wave velocity characteristics of the fault structure and their combined application in the detection process was conducted. Finally, the depth correction of the inversion result of the electrical resistivity tomography was carried out using the acoustic logging data; and the three-dimensional distribution position of the fault structure in the underground space was obtained. The
results show that the resistivity and acoustic wave velocity of the rock after fracture are lower than those before fracture. The fracture zone of the fault structure exists in the inversion section of the electrical resistivity tomography with a low resistivity anomaly region. On the wave velocity curve, the wave velocity value in the fracture zone decreases significantly. The combined detection mode of the two methods could provide a guidance for the delineation of the specific location of the fracture structure. The drilling results also verified that the inversion results of the electrical resistivity tomography corrected by acoustic logging data are more consistent with the actual geological conditions than those without correction.
-
涌水突泥灾害是威胁隧道施工安全最严重的地质灾害之一。开展隧址区涌水突泥致灾构造、富水特征和地下水循环过程研究是隧道涌水突泥灾害早期识别和涌水突泥灾害防范措施制定的基础[1-2]。断层类致灾构造是隧道涌水突泥致灾构造的主要类型[3] ，是涌水突泥灾害孕灾机制研究的重点之一。断裂的储水能力、物源补给能力是决定灾害威胁大小的重要因素，其中断裂的力学性质[4]、两盘的岩性组合[5]、断裂破碎带内物质组成[6]等方面会影响断裂带的导水、充水能力，物源补给能力与补给量、补给强度和持续性相关，因而断裂带构造裂隙水的补给、径流和排泄特征与断裂带致灾地质构造的规模和破坏性息息相关[7]，是涌水突泥灾害孕灾机制研究的重点之一。
川藏铁路西藏波密段位于高山峡谷区，山高谷深，平均高差超过2000 m。海拔4000 m以上的高山区有冰川覆盖。区内发育多条大型河流，水量十分充沛。河流受高位发育的冰川融水和大气降水补给，水力梯度大，水流湍急。该段铁路处于嘉黎—察隅断裂带影响范围内，多条河流与断裂发育位置重合，存在地表水与断裂带空间沟通并形成高压涌水突泥致灾构造的可能性。但由于该区山体陡峭、植被茂盛、地层岩性多样、构造条件复杂，缺乏对断裂带的导水或充水能力、地表水与地下水之间的关联机制、地下水循环过程等方面的认识，难以直接获取证据来开展涌水突泥灾害风险评价。
针对以上问题，本文采用多期测流和对大气降水、冰川、地表水和地下水水化学、同位素特征以及水文地球化学过程的分析，开展了西藏波密冰川覆盖区大型河流的流量变化特征和地表水与构造裂隙水的转化关系研究，对铁路隧道高压涌水突泥灾害的早期识别和工程预案设计具有一定指导作用。
1. 研究区概况
1.1 自然地理
西藏波密县位于西藏东南部高山峡谷地貌区。高山区海拔4000～6000 m，最高点海拔6643 m；河谷区海拔2000～2300 m，相对高差2000～4000 m[8]，地形高差大、岸坡陡峻。研究区属西藏东部温带山地湿润季风气候，年均气温8.7 ℃；1月气温最低，平均为−0.2 ℃；最高温度一般出现在7月，平均温度为16.4 ℃[9]。研究区降水充沛，多年平均降水量1276 mm，5—9月为雨季，占全年降水量的84.5%。
研究区整体属于雅鲁藏布江最大支流帕隆藏布流域。帕隆藏布发源于西藏八宿县阿扎贡拉山，自南东向北西穿过研究区，在林芝地区墨脱县甘登乡附近汇入雅鲁藏布江，年均流量约为3.12×1011 m3。隧址位于帕隆藏布右岸，共分布13条NE向大型河流和多条较小规模的河流[10]，自沟顶冰川覆盖区至沟口高差一般大于2000 m，平均坡度260‰，水量大，水流湍急（图1）。
![]()
图 1 研究区地质地貌与采样点位置图Figure 1. Geological geomorphology and location of the samplingof study area帕隆藏布两岸海拔4000 m以上的高山区有冰川覆盖，面积约324.72 km2。受热辐射和降水量变化影响[11-12]，近30年来冰川面积减少了58.2%，整体处于退缩状态[9, 11-12]，冰川融水量较大且逐年增加。冰川融水径流汇入沟中，为河流的重要补给项。
1.2 地质背景
研究区位于喜马拉雅东构造结东北部，一级构造单元属冈底斯—喜马拉雅造山系。主要断裂为NW向的嘉黎—察隅断裂，属全新世活动断裂[13]。小型断裂为张性和压扭性NE向断裂，其中古乡—随弄断裂、比通曲断裂、鸭容曲—索通断裂和赛隆卡曲断裂分别与古乡沟、比通曲、龙冲曲和赛隆卡曲的发育位置重合[14]。
研究区出露的基岩包括中—元古代念青唐古拉岩群片麻岩夹大理岩、燕山期侵入岩、变质岩、泥盆纪、石炭纪碎屑岩、碳酸盐岩和变质岩。出露的松散岩土体主要为更新世冰积碎石、全新世冲洪积碎石类土、泥石流堆积物等，主要分布于帕隆藏布河谷、冲沟低海拔处（图1）。
区内地下水类型分为松散岩类孔隙水和基岩裂隙水两大类。松散岩类孔隙水分布在帕隆藏布河谷、冲沟的河漫滩、阶地中，呈带状分布。基岩裂隙水分为风化裂隙水和构造裂隙水。风化裂隙水分布于花岗岩、片麻岩、石英岩等硬质岩的风化层；构造裂隙水呈带状或脉状分布于断裂带中，其补径排条件和赋存特征是本文讨论的主要内容。
2. 研究方法与数据获取
在研究区各河流近沟口位置选择水流较平稳、断面形状较规则的位置进行流量测定。共完成了3期测流工作：在2019年8月22—25日和2020年9月8—10日分别开展了雨季流量测定，在2020年4月22—25日开展了旱季流量测定，见图2。
![]()
图 2 研究区主要地表水流量统计Figure 2. Statistics of the main surface water flow in the study area在坎戈曲流域实施水文地质勘察钻孔1眼，孔深200 m，目的是对断裂带地下水赋存和径流特征进行研究。在该钻孔进行3个落程的抽水试验，获取了水文地质参数。
2019年8月和2020年4月，采集地表水样品各13件，采集当季降水样品各1件（2019年8月为降雨样品，2020年4月为降雪样品），采集冰川样品1件，采集泉点样品3件，结合水文地质勘察钻孔，采集埋深170 m地下水样品1件（图1）。采集样品时，使用意大利HANNA公司生产的多参数水质分析仪现场测试水样的温度、pH值、氧化还原电位（Eh）、溶解性总固体（TDS）等参数。对采集的样品进行全分析、δ2H、δ18O、氚同位素（3H）分析，利用14C同位素测定钻孔ZK01地下水样品的年龄。用于阳离子和微量元素分析的水样加入HNO3，使pH降到2以下，其余样品不加保护液。所有样品装瓶后均不留空气，拧紧瓶盖后用封口胶封闭送检。
水样全分析在自然资源部地下水矿泉水及环境监测中心完成。阳离子和阴离子测试分别采用ICP-AES（iCAP 6300, Thermo）和DX-120（Dionex）型离子色谱仪，执行DZ/T0064-1993和GB8538-2016标准，测试结果阴阳离子平衡相对误差小于±5%。δ2H、δ18O、3H测试在自然资源部地下水科学与工程重点实验室完成。δ2H和δ18O的测试利用波长扫描-光腔衰荡光谱仪（Picarro L2130i），测试精度分别为±1‰和±0.1‰。3H分析由超低本底液体闪烁谱仪（1220 Quantulus）计数测定，测试精度为±0.5TU。14C测年在美国Beta实验室利用加速器质谱（AMS）完成，测量精度为0.123%。
3. 结果与讨论
3.1 主要河流流量特征
在实际测定流量的17条主要河流中，2019年雨季河流流量为0.17×104～55.67×104 m3/d，平均值为17.96×104 m3/d，离散系数0.98，流量最大的是赛隆卡曲，最小的是优易沟。有5条河水流量超过30×104 m3/d，分别为扎隆沟、古乡沟、龙冲曲、茶隆隆巴曲和赛隆卡曲。
2020年雨季河流流量为0.10×104～40.68×104 m3/d，最大的为古乡沟，最小的为巴卡沟，平均值为15.03×104 m3/d，离散系数0.87，略小于2019年夏季。有4处流量超过30×104 m3/d，分别为古乡沟、龙冲曲、茶隆隆巴曲和赛隆卡曲。另外，比通曲流量超过了20×104 m3/d。
2020年旱季河流流量为0.10×104～23.86×104 m3/d，最大的为茶隆隆巴曲，最小为坎戈曲；平均值为5.72×104 m3/d，仅为2019年和2020年雨季的1/3左右。其中茶隆隆巴曲（23.86×104 m3/d）、赛隆卡曲（22.31×104 m3/d）和龙冲曲（10.36×104 m3/d）的河水流量超过10×104 m3/d，古乡沟（6.68×104 m3/d）、角隆曲（5.63×104 m3/d）和比通曲（4.30×104 m3/d）流量也较大。
对比连续两个雨季和一个旱季各河流流量变化，可以发现研究区河流流量在雨季和旱季差别明显。雨季河流流量与旱季河流流量之比为1.31～12.20。其中嘎朗沟、瓦普沟、达打隆巴等河流的流域面积和流量较小，但丰枯期流量差别较大；面积和流量较大的河流（如比通曲、龙冲曲、赛隆卡曲等与断裂复合发育的河流）的流量波动相对较小（图2），说明这些流域内的地质条件具有较强的地表水流量调蓄能力。
2019年8月和2020年9月的流量平均相差1.32倍。当日照、气温和降水量等气象条件相似时，地表水体得到的补给总量差别不大，说明气象条件对河流流量具有较强的影响。
3.2 水化学特征
地表水、地下水、冰川和降水的水化学指标统计分析见表1。雨、旱两季河流和地下水的pH值均为7.33～8.05，基本为中性偏碱性水；大气降水略偏酸性；冰川的pH值为7.45，基本呈中性。地表水各指标的变异系数范围为0.02～0.90，同一指标在2019年8月和2020年4月的变异系数相差不大，说明研究区河流的水化学指标随季节整体等幅变动。地下水各点中，既包括浅层地下水，又包括钻孔揭露的深层断裂带水，变异系数明显大于地表水，说明不同深度地下水循环过程存在差异。2019年雨季降水（雪）主要离子浓度小于2020年旱季，2019年8月和2020年4月溶解性总固体（TDS）分别为10.63 mg/L和30.07 mg/L，Ca2+浓度分别为1.07 mg/L和4.15 mg/L，
${\rm{SO}}_4^{2-} $ 浓度分别为1.41 mg/L和4.33 mg/L。受溶滤作用长期影响，冰川样品各指标浓度明显高于降水浓度。表 1 水样主要水化学指标值Table 1. Concentration statistics of the main hydrochemical indexes of water samples样品分类 项目 pH TDS K+ Na+ Ca2+ Mg2+ Cl− ${\rm{SO}}_4^{2-} $ ${\rm{HCO}}_3^- $ /（mg·L−1） 地表水—
2020年
4月最大值 7.99 551.00 5.16 2.84 147.40 23.07 1.75 353.50 158.50 最小值 7.33 84.00 1.29 0.85 21.95 1.88 1.40 15.19 53.84 平均值 7.57 248.53 2.65 1.53 62.65 7.79 1.69 118.35 85.86 标准差 0.18 143.92 1.08 0.58 37.50 5.21 0.13 102.28 27.91 中位数 7.57 210.00 2.24 1.36 46.16 6.65 1.75 90.23 83.75 变异系数 0.02 0.58 0.41 0.38 0.60 0.67 0.08 0.86 0.33 地表水—
2019年
8月最大值 8.05 196.17 4.39 1.30 56.20 8.69 0.20 80.31 173.20 最小值 7.48 37.09 0.40 0.47 8.62 0.98 0.05 4.29 27.83 平均值 7.83 96.18 1.64 0.72 26.20 3.46 0.09 25.15 64.26 标准差 0.19 55.89 1.15 0.22 16.10 2.26 0.04 22.53 41.16 中位数 7.90 83.82 1.45 0.69 22.30 3.06 0.08 21.67 52.57 变异系数 0.02 0.58 0.70 0.31 0.61 0.65 0.45 0.90 0.64 地下水 最大值 7.69 1325.00 6.28 3.30 340.00 23.04 2.79 848.60 191.40 最小值 7.50 176.00 4.69 1.02 48.23 3.51 1.75 32.03 122.60 平均值 7.58 612.50 5.62 1.94 158.11 14.59 2.18 332.11 142.25 标准差 0.08 462.82 0.58 0.86 118.45 8.55 0.45 336.14 28.48 中位数 7.57 474.50 5.76 1.73 122.11 15.90 2.10 223.91 127.50 变异系数 0.01 0.76 0.10 0.44 0.75 0.59 0.21 1.01 0.20 降水及
冰川最大值 7.45 67.37 0.55 5.24 12.57 1.59 1.75 8.50 53.78 最小值 6.38 10.63 0.39 0.10 1.07 0.08 1.40 1.41 6.52 平均值 6.80 36.02 0.49 1.85 5.93 1.06 1.52 4.75 31.07 标准差 0.47 23.54 0.07 2.40 4.86 0.69 0.16 2.91 19.34 中位数 6.56 30.07 0.54 0.21 4.15 1.51 1.40 4.33 32.90 变异系数 0.07 0.65 0.15 1.30 0.82 0.65 0.11 0.61 0.62 对比4种分类下各指标值的平均浓度，除Na++K+和Cl−外，各种离子浓度和TDS（溶解性总固体）都表现出地下水>地表水（2020年4月）>地表水（2019年8月）>降水及冰川。地表水在雨季与旱季离子指标存在差异的原因有两点：其一是雨季降水量较大，产生了稀释作用；其二是雨季地下水与含水介质之间反应时间有限[15]。
河水中主要阳离子为Ca2+和Mg2+，主要阴离子为
${\rm{HCO}}_3^- $ 和${\rm{SO}}_4^{2-} $ 。与河水相似，地下水的水化学类型主要为HCO3·SO4—Ca型，说明地下水与河水存在关联关系。大气降水的水化学类型为HCO3·SO4—Ca（2019年8月）和HCO3—Ca·Mg（2020年4月）型，冰川的水化学类型为HCO3—Ca型（图3）。大气降水及冰川水化学类型与地表水和地下水水化学组分存在差异，地下水和河水主要组分受岩石风化溶滤作用影响（图4），水化学组分主要来自于水—岩相互作用[16-17]。![]()
图 4 研究区河水和地下水吉布斯图Figure 4. Gibbs diagram of river water and groundwater in the study area地下水和地表水中主要阴离子为
${\rm{SO}}_4^{2-} $ 和${\rm{HCO}}_3^- $ ，而降水（雪、雨水）、冰川等补给源中的${\rm{SO}}_4^{2-} $ 和${\rm{HCO}}_3^- $ 浓度不高，说明在河水和地下水径流过程中，${\rm{SO}}_4^{2-} $ 和${\rm{HCO}}_3^- $ 发生了富集。有研究表明，青藏高原冰川径流中的${\rm{HCO}}_3^- $ 主要来自方解石的风化[18]，${\rm{SO}}_4^{2-} $ 主要来自黄铁矿氧化或石膏溶解[15,19]。图5（a）表明研究区内大理岩、灰岩的方解石溶解贡献了水中的${\rm{HCO}}_3^- $ ，图5（b）说明石膏溶解贡献了水中的${\rm{SO}}_4^{2-} $ 。研究区${\rm{SO}}_4^{2-} $ 的浓度远大于Fe，因此铁的硫化矿物的氧化作用对${\rm{SO}}_4^{2-} $ 的贡献小于石膏溶解。![]()
图 5 Ca2+/Na+和${\rm{HCO}}_3^- $ /Na+（a）与Ca2+和${\rm{SO}}_4^{2-} $ （b）毫克当量浓度关系图Figure 5. Milligram equivalent concentration diagram of (a) Ca2+/Na+ vs${\rm{HCO}}_3^- $ /Na and (b) Ca2+ vs${\rm{SO}}_4^{2-} $ 地表水和地下水中铁和锰的浓度较高，质量浓度范围分别为0.014～9.400 mg/L和0.006～0.343 mg/L，平均值分别为0.719 mg/L和0.032 mg/L。降水样品中铁和锰的浓度低。降水样品中铁低于0.001 mg/L，2019年8月降水中锰浓度为0.003 mg/L，2020年4月降水中未检出锰。在排除人类活动影响的前提下，地下水和河水中的铁和锰主要来源于水岩作用。沿嘉黎—察隅断裂带发育挤压破碎带，带内见糜棱岩、断层角砾岩、断层泥，含褐铁矿和锰质矿物[20]。通常还原环境有利于沉积物中Fe和Mn氧化物的溶解释放[21]，易溶的高价Fe和Mn氧化物会还原成溶解性更强的离子态Fe2+和Mn2+，提升水中铁和锰的浓度[22-23]；酸性环境更利于离子态Fe和Mn的富集[24]。通过图6可以看出，地表水和地下水分布在Fe2+和Fe（OH）3的界线附近，略偏向Fe2+区域；Mn主要以Mn2+的形式存在。
![]()
图 6 研究区水体Eh-pH值关系图Figure 6. Relationship between Eh and pH in the water samples in the study area研究区各条河水中铁和锰浓度存在较大差异，见图7。河水中Fe浓度较高的有古乡沟（9.400 mg/L）、来曲（1.822 mg/L）、龙冲曲（2.058 mg/L）和比通曲（1.381mg/L）和茶隆隆巴曲（0.896 mg/L），明显大于其他河水。河水中Mn浓度较大的有古乡沟（0.343 mg/L）、比通曲（0.083 mg/L）和龙冲曲（0.151 mg/L）。古乡沟、比通曲和龙冲曲等与断裂复合发育的地表水体的Fe和Mn质量浓度较高。
由此可以推断，大气降水和冰川融水补给地表水后，一部分转化为断裂带地下水，在还原环境下铁锰矿物溶解进入水中，使得地下水中的铁和锰浓度升高。当断裂带与地表水流复合发育，断裂带地下水循环时间长，与含水介质充分相互作用，铁锰浓度高；非断裂影响地区含水介质中铁锰浓度低。地下水最终排泄进入各地表水流中，并与地表水发生混合，汇入帕隆藏布。
3.3 同位素特征
研究区水体氢氧同位素分析测试结果见图8。所有水点均位于全球大气降水线（GMWL）和青藏高原东部大气降水线附近（青藏高原东部LMWL）[25]，说明降水是研究区各类型水体的主要补给来源[26-27]。地下水与2020年4月地表水点高度重合，说明地下水与地表水水力联系紧密。一般来说，在其他环境条件一定时，气温越低，δD和δ18O的值越小[28]。冰川形成时期气温低于现在的气温，图8中冰川点的δD、δ18O值相比于降水点更小。地表水和地下水是由冰川融水和大气降水共同补给形成的，δD和δ18O值位于在二者之间。
2020年4月气温低于2019年8月，但图8出现2020年4月各点的δD和δ18O的值大于2019年8月的异常现象，这可能是由于河水补给源比例不同造成的。为证明这一推测，选择龙冲曲作为典型流域，利用δD和δ18O计算雨季和旱季大气降水和冰川融水对地表水补给量的比重（图9），可以看出2019年8月冰川补给量占比为58.3%，2020年4月占比为38.5%。河水中冰川和大气降水的补给比例不同，造成不同时期地表水中δD和δ18O的值不同。由此可以看出，雨季冰川融水补给比例大于旱季。
利用氚同位素和14C同位素对地下水循环速率进行半定量分析（表2）。YGA10、YGA28和YGA34泉水样品的氚含量为3.9～6.2TU，根据Clark等[29]的氚方法定年解释，YGA10、YGA28泉水接受现代地下水补给，地下水年龄5～10 a；泉水YGA34接受现代地下水补给并混有少量1952年以前补给的地下水。3个泉点地下水循环时间短，更新速率较快，循环深度浅。
表 2 样品年代半定量对应表Table 2. Half quantitive corresponding table between content and age样品类型及编号 氚含量/TU 半定量年龄* 分类 泉-YGA10 6.2±0.9 5～10 a 浅循环地下水 泉-YGA28 4.1±0.8 5～10 a 浅循环地下水 泉-YGA34 3.9±0.8 1952年前补给与
5～10 a补给混合水浅循环地下水 钻孔样品 2.9±0.7 4640 a 中深循环地下水 地表水-YGA01 5.9±0.9 − 河水与地下水混合 地表水-YGB17 7.0±0.9 − 河水与地下水混合 地表水-YGB30 8.5±1.2 − 河水与地下水混合 注：*钻孔样品为14C测定的表现年龄。 在钻孔ZK01中约132 m埋深处采集了裂隙地下水，测得氚含量为2.9 TU，说明地下水的补给时间早于1952年。14C同位素测试得出该水样的表征年龄为4640 a。相比于YGA10、YGA28和YGA34泉水样品，钻孔中裂隙地下水的循环深度较大，更新速率较慢，属于中深循环的地下水。
YGA01（古乡沟）、YGB17（比通曲）和YGB30（龙冲曲）河水样品的氚含量为5.9～8.5 TU，高于泉点氚含量，但低于现代降雨中的氚含量（一般为101～102 TU[29]），说明河水中存在混入地下水的现象。对比地表水、泉和钻孔采集水样的年龄，可以推断，与断裂复合发育的地表水中有地下水的混合，其中以现代地下水为主。
3.4 冰川融水-大气降水-地下水-地表水转化过程
基于冰川融水、大气降水、河水和地下水等不同水体水化学和同位素组分的差异，分析研究区多种水体的转化过程。研究发现，西藏波密冰川覆盖区大型河流的补给受气象条件控制，径流过程中河水与地下水的转化过程受与河流复合发育的断裂控制。
区内河流发源于帕隆藏布右岸海拔4000 m以上的冰川覆盖区，河水和地下水水体δD和δ18O值均落在区域大气降水线附近，显示其主要接受冰川融水和大气降水补给。在雨季，河水的δD和δ18O的值偏向冰川，而在旱季河水的δD和δ18O的值偏向大气降水，说明河水补给来源的比例在旱季和雨季存在差异。旱季大气降水补给量占比大，雨季冰川融水补给量占优。这可能是由于旱季属全年气温较低的时期，冰川融水量小。因而大气降水占总补给量的比例超过60%。雨季集中在气温较高的时期，冰川融水量激增，对河水补给的贡献超过大气降水，约占58.3%。
河水通常在高海拔处接受补给，因而河水和地下水相对排泄基准面（帕隆藏布）具有较高的水头。河水沿各沟谷自北西向南东方向径流，断裂带地下水在河谷附近或以泉的形式、或以线状排泄的方式排泄，最终汇入帕隆藏布。
河流自高位补给到汇入帕隆藏布的径流过程中，存在2种不同的水体转化过程类型（图10）。
![]()
图 10 不同类型水体转化过程模式图Figure 10. Transformation process model of different types of water第一种类型是古乡沟、龙冲曲、比通曲和赛隆卡曲等与断裂复合发育的河流。研究区内所有河流2019年雨季与2020年旱季流量之比为1.71～11.15，2020年雨季与2020年旱季流量之比为1.31～12.20。而与断裂复合发育的河流2019年雨季与2020年旱季流量之比为2.50～7.28，2020年雨季与2020年旱季流量之比为1.61～6.09，明显小于全区河流的同类流量比值。说明以上几条河流的流量较为稳定。其原因为：在气温、降水量一定的前提下，断裂带地下水对河水持续补给，减小了河水流量在雨季和旱季的变化幅度。
与断裂复合发育的河流流域内岩体完整性差，利于地表水入渗。地下水接受河水持续的补给，方向基本与地形坡度保持一致。地下水与含水介质发生水岩作用较为充分，TDS值升高，方解石的风化和石膏的溶解使水中
${\rm{SO}}_4^{2-} $ 和${\rm{HCO}}_3^- $ 富集，水化学类型由补给源的HCO3—Ca（Mg）型逐渐变为HCO3·SO4—Ca型。裂隙含水介质的还原环境也使地下水中铁、锰浓度升高。大部分地下水循环深度小，年龄为5～10 a，更新速率快，与河流存在密切的水力联系和相互转化。少部分地下水沿断裂向深部径流，地下水年龄在4000 a以上。即便在断裂带渗透性较好的条件下，参与中深循环的地下水更新速率也较慢，水岩作用较充分，TDS值明显高于河水和浅层循环地下水。地下水与地表水混合后汇入帕隆藏布。
第二种类型是非断裂影响区的河流。基岩风化裂隙水接受河水渗漏补给。风化裂隙发育深度有限，地下水对河水的补给量和水量调节能力小于第一种类型，河水在雨季和旱季流量比值大。较弱的还原环境不利于地下水中铁、锰溶解，地下水中铁、锰浓度低。地下水年龄5～10 a，循环深度浅。临近排泄基准面的河水水化学类型为HCO3—Ca型和HCO3·SO4—Ca型，说明混入了少量浅循环的地下水。
从工程建设防灾减灾的角度考虑，第一种类型的河流，由于水量大且较为稳定，水头高，断裂带地下水与地表水存在较好的水力联系，因而在施工过程中产生高压涌水突泥灾害的风险较大。在施工过程中，应充分予以关注，充分考虑复杂地质构造的致灾风险。
4. 结论
（1）西藏波密冰川覆盖区河水主要接受冰川融水和大气降水补给。雨季河水的 δ 18O和 δ D值小于旱季，说明河水雨季和旱季的补给源结构不同。在雨季冰川融水为主要补给源，旱季以大气降水补给为主。
（2）河水成分及来源与断裂有密切关系。根据河水与断裂空间展布关系，河流流域内水体径流过程可分为2种类型。其一是古乡沟、比通曲、龙冲曲等与断裂带复合发育的河流，可补给断裂带裂隙地下水。浅层循环断裂带水年龄5～10 a，循环速率快；中深层断裂带水年龄超过4000 a，循环速率慢，水岩作用较充分。由于持续接受地下水补给，年内河水流量波动较小。其二是非断裂影响区内的河流，河水可补给风化裂隙地下水。地下水循环快，与河水交换较频繁，水岩作用程度弱。年内流量波动较大。河水与地下水混合后，排泄至帕隆藏布。
（3）研究区河水补给区海拔高，河水流量大，且与断裂带地下水存在较好的水力联系，在地下工程建设过程中，与断裂带复合发育的古乡沟、比通曲、龙冲曲存在高压涌水突泥灾害风险，应予以高度关注。
-
![]()
图 3 L1测线高密度电阻率法反演剖面图
Figure 3. Profile of electrical resistivity tomography inversion for L1 survey
![]()
图 4 ZK-1声波波速曲线与钻孔资料图
Figure 4. Acoustic velocity curve of ZK-1 compared with drilling data
![]()
图 5 断层破碎带电阻率和声波波速曲线对比图
Figure 5. Comparison of resistivity and acoustic velocity curves of a fractured fault zone
![]()
图 6 断层破碎带处不同垂直滤波比系数的电阻率曲线图
Figure 6. Resistivity with different vertical filtering coefficient values at the fractured fault zone
![]()
图 7 断层破碎带处不同垂直滤波比系数的电阻率曲线和声波波速曲线对比图
Figure 7. Comparison of resistivity with different vertical filtering coefficient values in the fractured fault zone and the acoustic wave velocity
![]()
图 8 系数为1.6的L1测线高密度电阻率反演剖面图
Figure 8. Profile of electrical resistivity tomography inversion for L1 survey with a coefficient of 1.6
表 1 岩石破裂前后的电阻率与声波波速变化特征
Table 1 Characteristics of resistivity and acoustic wave velocity changes in the rocks before and after hydraulic fracturing
岩石
样品岩样
直径
/mm岩样
长度
/mm破裂前 破裂后 电阻率值
/（Ω·m）声波波速值
/（km·s−1）电阻率值
/（Ω·m）声波波速值
/（km·s−1）1 50.50 110.40 1261 6.42 86 5.16 2 50.63 111.54 1892 5.36 139 4.28 3 50.60 110.38 1182 5.21 171 4.72 4 50.47 110.41 1291 5.02 210 4.68 
下载: 导出CSV
表 2 KV与岩体完整程度的对应关系[26]
Table 2 Correspondence between KV and the integrity degree of rock mass
KV >0.75 >0.55～0.75 >0.35～0.55 0.15～0.35 ＜0.15 岩体完整程度 完整 较完整 较破碎 破碎 极破碎
下载: 导出CSV
表 3 声波测井结果
Table 3 Results of acoustic well logging
孔深/m 层厚/m 岩性 声波波速值/（km·s−1） KV值 破碎程度 极小值 极大值 平均值 7.00 10.40 3.40 砾岩 2.00 3.05 2.54 0.26 破碎 19.80 9.40 石英砂岩 3.11 5.00 3.83 0.61 较完整 31.20 11.40 闪长玢岩 1.98 2.48 2.07 0.15 极破碎
下载: 导出CSV
-
[1] FAVREAU P,WOLF S. Theoretical and numerical stress analysis at edges of interacting faults:application to quasi-static fault propagation modelling[J]. Geophysical Journal International,2009,179(2):733 − 750. DOI: 10.1111/j.1365-246X.2009.04314.x
[2] PERRIN C,WALDHAUSER F,SCHOLZ C H. The shear deformation zone and the smoothing of faults with displacement[J]. Journal of Geophysical Research:Solid Earth,2021,126(5):e2020JB020447. DOI: 10.1029/2020JB020447
[3] 周保,孙皓,魏赛拉加,等. 冷湖赛什腾山天文台工程地质选址评价[J]. 水文地质工程地质,2022,49(4):165 − 174. [ZHOU Bao,SUN Hao,WEI Sailajia,et al. Site selection evaluation of engineering geology of the Lenghu Saishiteng Mountain Observatory[J]. Hydrogeology & Engineering Geology,2022,49(4):165 − 174. (in Chinese with English abstract)] ZHOU Bao, SUN Hao, WEI Sailajia, et al . Site selection evaluation of engineering geology of the Lenghu Saishiteng Mountain Observatory[J]. Hydrogeology & Engineering Geology,2022 ,49 （4 ）:165 −174 . （in Chinese with English abstract）[4] 刘明辉,薛建,王者江,等. 工程场地隐伏断裂的探测与地震活动性评价[J]. 物探与化探,2018,42(4):839 − 845. [LIU Minghui,XUE Jian,WANG Zhejiang,et al. The detection of buried faults in engineering sites and the evaluation of seismic activity[J]. Geophysical and Geochemical Exploration,2018,42(4):839 − 845. (in Chinese with English abstract)] LIU Minghui, XUE Jian, WANG Zhejiang, et al . The detection of buried faults in engineering sites and the evaluation of seismic activity[J]. Geophysical and Geochemical Exploration,2018 ,42 （4 ）:839 −845 . （in Chinese with English abstract）[5] 李坚,邓宏科,张家德,等. 可控源音频大地电磁勘探在大瑞铁路高黎贡山隧道地质选线中的应用[J]. 水文地质工程地质,2009,36(2):72 − 76. [LI Jian,DENG Hongke,ZHANG Jiade,et al. Application of geological routing about CSAMT exploration in Gaoligong Mountain tunnel of Dali-Ruili Railway[J]. Hydrogeology & Engineering Geology,2009,36(2):72 − 76. (in Chinese with English abstract)] LI Jian, DENG Hongke, ZHANG Jiade, et al . Application of geological routing about CSAMT exploration in Gaoligong Mountain tunnel of Dali-Ruili Railway[J]. Hydrogeology & Engineering Geology,2009 ,36 （2 ）:72 −76 . （in Chinese with English abstract）[6] 杨志华,郭长宝,吴瑞安,等. 青藏高原巴塘断裂带地震滑坡危险性预测研究[J]. 水文地质工程地质,2021,48(5):91 − 101. [YANG Zhihua,GUO Changbao,WU Ruian,et al. Predicting seismic landslide hazard in the Batang fault zone of the Qinghai-Tibet Plateau[J]. Hydrogeology & Engineering Geology,2021,48(5):91 − 101. (in Chinese with English abstract)] YANG Zhihua, GUO Changbao, WU Ruian, et al . Predicting seismic landslide hazard in the Batang fault zone of the Qinghai-Tibet Plateau[J]. Hydrogeology & Engineering Geology,2021 ,48 （5 ）:91 −101 . （in Chinese with English abstract）[7] SANA H,TABORIK P,VALENTA J,et al. Detecting active faults in intramountain basins using electrical resistivity tomography:A focus on Kashmir Basin,NW Himalaya[J]. Journal of Applied Geophysics,2021,192:104395. DOI: 10.1016/j.jappgeo.2021.104395
[8] 汪佳蓓,陈青云,胡怡,等. 深部煤层小断层断距识别精度研究[J]. 地球物理学进展,2022,37(3):1159 − 1169. [WANG Jiabei,CHEN Qingyun,HU Yi,et al. Study on the identification accuracy of small fault displacement in deep coal bed[J]. Progress in Geophysics,2022,37(3):1159 − 1169. (in Chinese with English abstract)] WANG Jiabei, CHEN Qingyun, HU Yi, et al . Study on the identification accuracy of small fault displacement in deep coal bed[J]. Progress in Geophysics,2022 ,37 （3 ）:1159 −1169 . （in Chinese with English abstract）[9] 付光明,黄进调,刘阳,等. 高密度电阻率法与CSAMT法在江西会昌县坝背地热勘查中的综合探测[J]. 中国地质,2019,46(4):927 − 936. [FU Guangming,HUANG Jintiao,LIU Yang,et al. Multi-electrode resistivity method and CSAMT method in geothermal exploration of Babei Area in Huichang County,Jiangxi Province[J]. Geology in China,2019,46(4):927 − 936. (in Chinese with English abstract)] FU Guangming, HUANG Jintiao, LIU Yang, et al . Multi-electrode resistivity method and CSAMT method in geothermal exploration of Babei Area in Huichang County, Jiangxi Province[J]. Geology in China,2019 ,46 （4 ）:927 −936 . （in Chinese with English abstract）[10] SEMINSKY K Z,ZARIPOV R M,OLENCHENKO V V. Interpretation of shallow electrical resistivity images of faults:tectonophysical approach[J]. Russian Geology and Geophysics,2016,57(9):1349 − 1358. DOI: 10.1016/j.rgg.2016.08.020
[11] 李华,王东辉,张伟,等. 地球物理探测技术在成都市浅表地质结构调查中的应用研究[J]. 中国地质,2022,49(5):1438 − 1457. [LI Hua,WANG Donghui,ZHANG Wei,et al. Application research of geophysical exploration technology in the investigation of shallow geological structure in Chengdu[J]. Geology in China,2022,49(5):1438 − 1457. (in Chinese with English abstract)] LI Hua, WANG Donghui, ZHANG Wei, et al . Application research of geophysical exploration technology in the investigation of shallow geological structure in Chengdu[J]. Geology in China,2022 ,49 （5 ）:1438 −1457 . （in Chinese with English abstract）[12] 周越,曾昭发,唐海燕,等. 公路勘察中滑坡体的地球物理特征与分析——以张榆线公路勘察为例[J]. 吉林大学学报(地球科学版),2021,51(2):638 − 644. [ZHOU Yue,ZENG Zhaofa,TANG Haiyan,et al. Geophysical Characteristics of Landslide Body in Highway Reconnaisance:A Case Study in Highway Prospecting of Zhangyu Line. Journal of Jilin University(Earth Science Edition),2021,51(2):638 − 644. (in Chinese with English abstract)] ZHOU Yue, ZENG Zhaofa, TANG Haiyan, et al. Geophysical Characteristics of Landslide Body in Highway Reconnaisance: A Case Study in Highway Prospecting of Zhangyu Line. Journal of Jilin University（Earth Science Edition）, 2021, 51（2）: 638 − 644. （in Chinese with English abstract）
[13] 方艺翔,李卓,范光亚,等. 监测资料、压水试验与综合物探法在某心墙坝渗漏识别中的应用研究[J]. 水利水电技术(中英文),2022,53(2):87 − 97. [FANG Yixiang,LI Zhuo,FAN Guangya,et al. Application of monitoring data,water pressure test and comprehensive geophysical method to identification of leakage of a core dam[J]. Water Resources and Hydropower Engineering,2022,53(2):87 − 97. (in Chinese with English abstract)] FANG Yixiang, LI Zhuo, FAN Guangya, et al . Application of monitoring data, water pressure test and comprehensive geophysical method to identification of leakage of a core dam[J]. Water Resources and Hydropower Engineering,2022 ,53 （2 ）:87 −97 . （in Chinese with English abstract）[14] 丁超,解阳波,张家豪,等. 三维高密度电法揭示SAGD地面窜漏通道[J]. 吉林大学学报(地球科学版),2022,52(6):2021 − 2033. [DING Chao,XIE Yangbo,ZHANG Jiahao,et al. Imaging the ground leakage channel of SAGD based on three-dimensional electrical resistivity tomography[J]. Journal of Jilin University (Earth Science Edition),2022,52(6):2021 − 2033. (in Chinese with English abstract)] DING Chao, XIE Yangbo, ZHANG Jiahao, et al . Imaging the ground leakage channel of SAGD based on three-dimensional electrical resistivity tomography[J]. Journal of Jilin University （Earth Science Edition）,2022 ,52 （6 ）:2021 −2033 . （in Chinese with English abstract）[15] 陈松,陈长敬,罗士新,等. 广州南沙厚覆盖区近地表地层结构分析:基于二维地球物理多方法探测结果[J]. 地质通报,2023,42(1):168 − 179. [CHEN Song,CHEN Changjing,LUO Shixin,et al. Near surface stratigraphic structure analysis in Nansha of Guangzhou thick overburden area:Based on multi method results of 2D geophysics exploration[J]. Geological Bulletin of China,2023,42(1):168 − 179. (in Chinese with English abstract)] CHEN Song, CHEN Changjing, LUO Shixin, et al . Near surface stratigraphic structure analysis in Nansha of Guangzhou thick overburden area: Based on multi method results of 2D geophysics exploration[J]. Geological Bulletin of China,2023 ,42 （1 ）:168 −179 . （in Chinese with English abstract）[16] 贾龙,雷明堂,程小杰. 基于井中超声波成像的岩溶特征高精度探测和评价[J]. 地质通报,2022,41(增刊1):453 − 460. [JIA Long,LEI Mingtang,CHENG Xiaojie. High precision detection and evaluation of karst features based on borehole ultrasonic imaging[J]. Geological Bulletin of China,2022,41(Sup 1):453 − 460. (in Chinese with English abstract)] JIA Long, LEI Mingtang, CHENG Xiaojie. High precision detection and evaluation of karst features based on borehole ultrasonic imaging[J]. Geological Bulletin of China, 2022, 41（Sup 1）: 453 − 460. （in Chinese with English abstract）
[17] 胡刚,何正勤,李娜,等. 跨安宁河断裂带浅孔综合地球物理测井成果分析[J]. 地震学报,2016,38(5):684 − 692. [HU Gang,HE Zhengqin,LI Na,et al. Analyses on integrated geophysical logging results in shallow hole across the Anninghe fault zone[J]. Acta Seismologica Sinica,2016,38(5):684 − 692. (in Chinese with English abstract)] HU Gang, HE Zhengqin, LI Na, et al . Analyses on integrated geophysical logging results in shallow hole across the Anninghe fault zone[J]. Acta Seismologica Sinica,2016 ,38 （5 ）:684 −692 . （in Chinese with English abstract）[18] 王新杰,杨进平,王营超,等. 声波测试技术在地下油库选址勘查中的应用[J]. 工程地球物理学报,2013,10(5):648 − 651. [WANG Xinjie,YANG Jinping,WANG Yingchao,et al. The application of acoustical testing technique in site selection of oil cellar[J]. Chinese Journal of Engineering Geophysics,2013,10(5):648 − 651. (in Chinese with English abstract)] WANG Xinjie, YANG Jinping, WANG Yingchao, et al . The application of acoustical testing technique in site selection of oil cellar[J]. Chinese Journal of Engineering Geophysics,2013 ,10 （5 ）:648 −651 . （in Chinese with English abstract）[19] 朱飞飞. 地井联合物探技术在岩溶注浆检测中的应用[J]. 工程地球物理学报,2022,19(4):450 − 458. [ZHU Feifei. Application of geophysical prospecting technology combined with ground and well in Karst grouting detection[J]. Chinese Journal of Engineering Geophysics,2022,19(4):450 − 458. (in Chinese with English abstract)] ZHU Feifei . Application of geophysical prospecting technology combined with ground and well in Karst grouting detection[J]. Chinese Journal of Engineering Geophysics,2022 ,19 （4 ）:450 −458 . （in Chinese with English abstract）[20] 易强,李望明,景营利,等. 高密度电法数据处理结合钻孔波速测试在土洞发育区的应用实例[J]. 物探化探计算技术,2021,43(1):104 − 107. [YI Qiang,LI Wangming,JING Yingli,et al. The application of high-density resistivity method data processing combined with borehole wave velocity measurement in soil hole development area[J]. Computing Techniques for Geophysical and Geochemical Exploration,2021,43(1):104 − 107. (in Chinese with English abstract)] YI Qiang, LI Wangming, JING Yingli, et al . The application of high-density resistivity method data processing combined with borehole wave velocity measurement in soil hole development area[J]. Computing Techniques for Geophysical and Geochemical Exploration,2021 ,43 （1 ）:104 −107 . （in Chinese with English abstract）[21] 陈清. 高密度电法与跨孔层析成像在岩溶探测中的 应用分析[J]. 工程地球物理学报,2017,14(4):427 − 434. [CHEN Qing. Application of high density electrical method and cross-hole tomography imaging to Karst detection[J]. Chinese Journal of Engineering Geophysics,2017,14(4):427 − 434. (in Chinese with English abstract)] CHEN Qing . Application of high density electrical method and cross-hole tomography imaging to Karst detection[J]. Chinese Journal of Engineering Geophysics,2017 ,14 （4 ）:427 −434 . （in Chinese with English abstract）[22] MENG Fangsong,ZHANG Gang,QI Yaping,et al. Application of combined electrical resistivity tomography and seismic reflection method to explore hidden active faults in Pingwu,Sichuan,China[J]. Open Geosciences,2020,12(1):174 − 189. DOI: 10.1515/geo-2020-0040
[23] 潘剑伟,占嘉诚,洪涛,等. 地面核磁共振方法和高密度电阻率法联合找水[J]. 地质科技情报,2018,37(3):253 − 262. [PAN Jianwei,ZHAN Jiacheng,HONG Tao,et al. Combined use of surface nuclear magnetic resonance and electrical resistivity imaging in detecting groundwater[J]. Geological Science and Technology Information,2018,37(3):253 − 262. (in Chinese with English abstract)] PAN Jianwei, ZHAN Jiacheng, HONG Tao, et al . Combined use of surface nuclear magnetic resonance and electrical resistivity imaging in detecting groundwater[J]. Geological Science and Technology Information,2018 ,37 （3 ）:253 −262 . （in Chinese with English abstract）[24] 楚泽涵,黄隆基,高杰,等. 地球物理测井方法与原理[M]. 北京:石油工业出版社,2007:225 − 227. [CHU Zehan,HUANG Longji,GAO Jie,et al. Geophysical logging methods and principles[M]. Beijing:Petrpleum Industry Press,2007:225 − 227. (in Chinese)] CHU Zehan, HUANG Longji, GAO Jie, et al. Geophysical logging methods and principles[M]. Beijing: Petrpleum Industry Press, 2007: 225 − 227. （in Chinese）
[25] 谭礼洪,张国强,谭忠健,等. 利用阵列声波测井资料评价变质岩储层有效性[J]. 石油地球物理勘探,2022,57(6):1464 − 1472. [TAN Lihong,ZHANG Guoqiang,TAN Zhongjian,et al. Evaluation of metamorphic rock reservoir effectiveness by array acoustic logging data[J]. Oil Geophysical Prospecting,2022,57(6):1464 − 1472. (in Chinese with English abstract)] TAN Lihong, ZHANG Guoqiang, TAN Zhongjian, et al . Evaluation of metamorphic rock reservoir effectiveness by array acoustic logging data[J]. Oil Geophysical Prospecting,2022 ,57 （6 ）:1464 −1472 . （in Chinese with English abstract）[26] 中华人民共和国住房和城乡建设部. 工程岩体分级标准:GB/T 50218—2014[S]. 北京:中国计划出版社,2015. [Ministry of Housing and Urban-Rural Development of the People’s Republic of China. Standard for engineering classification of rock mass:GB/T 50218—2014[S]. Beijing:China Planning Press,2015. (in Chinese)] Ministry of Housing and Urban-Rural Development of the People’s Republic of China. Standard for engineering classification of rock mass: GB/T 50218—2014[S]. Beijing: China Planning Press, 2015. （in Chinese）
-
期刊类型引用(3)
1. 邵杰，滕超，陈喜庆，杨欣杰，曹军，朱宁，肖登，吕菲. 藏东川西交通廊道波密至林芝段重大工程水文地质条件及问题研究. 河北工程大学学报(自然科学版). 2023(02): 105-112 . 
百度学术
2. 张春潮，李向全，马剑飞，付昌昌，白占学，余启明. 喜马拉雅东构造结地热资源赋存特征与开发利用潜力. 地质学报. 2023(08): 2728-2741 . 百度学术
3. 张田田，杨为民，万飞鹏. 浑河断裂带地质灾害发育特征及其成因机制. 吉林大学学报(地球科学版). 2022(01): 149-161 . 百度学术
其他类型引用(2)
计量
- 文章访问数: 111
- HTML全文浏览量: 13
- PDF下载量: 36
- 被引次数: 5
邮件订阅
RSS