作者简介:王晓明(1964-),男,高级工程师,长期从事水工环地质勘查、评价、研究工作。E-mail: cnahwxm@163.com
地面核磁共振找水技术是一种非介入探测技术,通过识别水分子中氢原子信号的分布来推测水分子在地层中的分布,进而确定是否有地下水存在,并给出地下水位、富水程度和孔隙率等地层结构信息。以安徽某水库土质大坝为例,采用地面核磁共振方法对大坝渗漏进行勘查,取得初始振幅( E0)、纵向弛豫时间( T1)、横向弛豫时间(
Surface nuclear magnetic resonance(SNMR) technique is a non-invasive detection technology, which, through the recognition of the distribution of hydrogen atom signals in water molecules, infers the distribution of water molecules in the strata, and then determines whether there exists groundwater or not and provides some information concerning the water table, the water-rich degree and the porosity. In this paper, taking a soil dam of a reservoir in Anhui Province as an example, the authors used the SNMR to investigate the dam leakage, and obtained such physical parameters as initial amplitude ( E0), longitudinal relaxation time ( T1), and transverse relaxation time (
我国在20世纪50年代末至70年代修建了大量的平原和山区水库, 很大一部分属边勘测、边设计、边施工的水利工程, 限于材料和技术等条件, 水库坝体以土石坝居多, 经过几十年的运行加之管理粗放、维护较少, 水库坝体、坝基出现渗漏的现象较为普遍, 坝体的稳定性存在较大隐患, 进而威胁水库下游的生命与财产安全。据不完全统计, 全国3 100多座大中型水库中, 出现渗漏的坝体与坝基约占20%; 小型水库大坝也约有40%存在渗漏问题[1]。目前, 探测大坝渗漏点的技术方法有很多, 一般是通过地球物理资料反演间接推断坝体、坝基的含水信息, 再推测渗漏介质和渗漏途径[2, 3]。而地面核磁共振(Surface nuclear magnetic resonance, SNMR)方法是目前可以直接探测地下水的物探新方法, 通过对物理参数与地质结构对应关系的研究, 做出大坝渗漏位置、渗漏量的判断[4]。本次选取安徽省某水库作为地面核磁共振方法研究点, 以检验SNMR方法在大坝渗漏勘查中的应用效果。
地面核磁共振方法探测的目标是水中的氢质子。在地磁场B0的作用下, 氢质子以Larmor角频率ω L=γ p|B0|=2π fL( 其中γ p为质子旋磁比)进行旋进运动, 产生沿地磁场方向的磁化强度M0(见图1a)[5]。
世界范围内, 地磁场强度B0为赤道的25 000 nT到两极的65 000 nT, 即Larmor 频率范围为1.0~3.1kHz。当铺设在地面上的线圈内通入Larmor 频率fL 的交变电流I(t)=I0cos(2π fLt)脉冲时, 在地面以下产生高强度非均匀的激发磁场BT。激发磁场BT的能量由激发脉冲矩q所决定, q是激发电流脉冲的幅值I0和激发脉冲持续时间τ 的乘积(τ 一般为40 ms)。
激发磁场BT垂直于地磁场B0方向, 则对自旋的氢质子施加了一个磁矩, 使其旋转轴偏离地磁场B0方向(见图1b), 导致平行于B0方向的纵向磁化强度减弱, 而垂直于B0方向的横向磁化强度得到增强。当切断产生BT的电流脉冲时, 核自旋磁化强度的纵向分量和横向分量按指数规律(纵向分量增强, 横向分量衰减)变化至各自的原始状态(见图1c)。而横向分量衰减的弛豫现象产生一个微弱的, 但可检测到的感应宏观交变磁场, 用铺设在地面上的同一线圈采集随时间衰减的感应电动势, 即核磁共振(NMR)信号[4]。
应用一个脉冲矩q进行激发, 只能使地面以下某一深度区域内氢质子的自旋方向由平衡状态发生倾斜。从小到大选择合适的一系列脉冲矩q, 可由浅到深地实现探测地下各深度位置的含水层分布。在每次测量过程中, 记录所有按指数规律衰减的核磁共振响应信号, 通过叠加和滤波等数据处理方法后, 对信号进行参数提取, 获得NMR信号的初始振幅E0、衰减时间常数T
在每个NMR测点上是通过由小到大(100~9 000 A· ms)依次改变激发脉冲矩 (q=I0 · tp, I0、tp分别为激发电流脉冲的幅值和持续时间)参数, 观测和研究NMR信号的变化规律来了解由坝基浅到深各层的含水量[9]。通过分析坝基含水量的变化, 可分析出坝基泄露点位置。
本次选用法国IRIS公司研制的NUMIS系统开展研究, 该系统由直流电源(汽车电瓶2块)、DC/DC转换器、发射机、发射/接收天线、调谐单元、接收机、外接计算机及部分辅助连线构成。
研究水库属长江流域滁河水系, 设计洪水标准重现期30年一遇, 校核洪水标准重现期300年一遇, 水库大坝长350 m, 坝顶高程89.68~90.2 m, 顶宽4 m, 最大坝高12.7 m。蓄水位86.4 m, 设计洪水位87.48 m, 总库容154 万m3。水库枢纽工程由大坝、溢洪道和高、低放水涵组成。大坝坝体为土质坝, 北部坝基为岩体花岗二长岩, 南部坝基为砂岩。该水库1976年完成, 在运行中多次出现严重渗漏。2009年曾采用多头小直径深层搅拌喷灌浆建造防渗墙对水库大坝进行处理, 桩径为380 mm, 钻杆间距320 mm, 最小成墙厚度205 mm, 平均钻入深度为11.4 m, 取得了一定的防渗效果。
库区及周边地貌类型为浅丘, 地形呈南北高, 中间低, 标高约80~90 m。地层属华南地层大区下扬子地层区, 主要岩性有第四系上更新统下蜀组粘性土、新近系上新统黄岗组的玄武岩及古近系古新统的舜山集组的砂砾岩、泥砂岩、含钙砂岩, 并出露燕山期侵入的花岗二长岩岩体(图2)。
区内地下水类型划分为松散岩类孔隙水、碎屑岩类孔隙裂隙水、基岩裂隙水和玄武岩孔洞裂隙水四种类型。松散岩类孔隙水分布在大坝下及水库下游, 含水岩组主要由第四系上更新统粉质黏土组成。单井涌水量10~100 m3/d。碎屑岩类孔隙裂隙水分布在大坝西侧和水库南部, 含水岩组主要由古近系的砾岩、砂岩、泥岩、页岩组成, 单井涌水量< 10 m3/d。基岩裂隙水的含水岩组由燕山期侵入岩组组成, 分布在大坝北侧和水库北部, 地下水主要赋存于块状岩类的风化裂隙、构造裂隙、层间裂隙中, 还有一些与断层直接相关的脉状水。单井涌水量10~100 m3/d。玄武岩孔洞裂隙水的主要含水岩组由新近系上新统黄岗组玄武岩组成, 主要分布于大坝南侧。地下水主要赋存于玄武岩气孔状、杏仁状构造及垂直节理、水平“ 层理” 中, 单井涌水量100~1 000 m3/d。
首先确定测区激发频率, 利用NUMIS系统进行测量。利用磁力仪测量并算出该工区地磁场平均感应磁化强度B0=50 258.22 nT, 据此确定NUMIS系统发射频率为2 141.35 Hz。
在水库主坝体上布设一条SNMR观测剖面, 剖面上敷设3个SNMR测深点:H1、H2、H3(见图2)。平均信噪比较大(1.18), 每个点使用 “ 8” 字形线圈。根据工区地质背景和工作条件, 确定NUMIS系统的记录长度为240 ms, 脉冲的持续时间置成40 ms, 脉冲间歇时间为30 ms; 选用16个脉冲矩, 测量数值范围是在500~60 000 nV之间设置; 叠加次数为32~64次。
利用NUMIS系统自带程序和预处理系统对记录下NMR信号各参数进行数据处理和反演解释。表1为测量参数与反演所得水文地质参数的对应关系。
自由水
根据核磁共振测深各测点的反演解释结果, 可以得到各测点的含水量直方图(含
由H1点含水量直方图(图3)可见, 在深度2 m以内
图4显示, 在深度1~2 m范围内平均含水量为9%,
图5可见, 在深度1~4 m附近区域, 含水量平均等于10%,
根据3个核磁共振测点所得含水量直方图数据, 人工绘制含水量断面等值线图, 并与钻孔揭露的地层进行对比(图6)。
由图6可以得到如下结论:在埋深30 m以内, 含水量大于10%的区域有4处, 编号L1~L4, L1、L4现场确认为放水涵。L2是渗漏点, 推测为大坝放水涵下部断裂所致, 埋深6~9 m; L3为渗漏带, 埋深主要在10~25.5 m, 长度约为180 m, 基本对应了花岗岩侵入接触带和岩石风化带上。
通过钻孔和野外调查也揭露相应的泄露点和渗漏带(图7), 不含水的地方为水泥混合土, 钻孔取出的岩性较为完整, 室内试验的渗透系数在9× 10-7cm/s。说明地面核磁共振方法确定堤坝渗漏点(带)是有效的。
地面核磁共振方法探测土石堤坝的完好性(例如坝体、坝基渗漏), 是一种直接、准确、有效的方法, 已为现场调查与钻探资料所验证。坝体、坝基的渗漏与地层含水性、渗透性和孔隙性密切相关, 实际测量的E0、T1、
目前, 地面核磁共振方法的数据处理和反演解释还是一维的, 水文地质解释是建立在地表和地下介质水平、均质, 且各项同性的基础上, 研究成果还有一定的局限性; 而二、三维数据处理、反演和水文地质解释, 随着SNMR技术的改进而完善, 将进一步提高新方法的勘查精度。
致谢:感谢安徽省地勘局科技信息处的大力支持,感谢中国地质大学(武汉)地空学院李振宇教授的辛苦工作和提供的地面核磁共振勘查资料,感谢项目组无私提供完整报告及原始资料供研究参考。
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|