地下核磁共振小尺寸线圈设计和实验

很多国家在进行隧道挖掘、地下矿体开采时,经常会遇到突水涌泥事故.核磁共振方法（Magnetic Resonance Sounding,MRS）作为目前唯一一种直接找水的地球物理技术,在探测诱发灾害的水体方面具有明显的优势.本文首先给出了多匝分离式线圈的MRS信号响应和核函数的计算公式,分析了线圈的匝数、边长和含水层位置对MRS信号响应的影响.分析表明,在限定空间内,利用多匝分离式线圈的设计方法可以提高小尺寸线圈的探测能力.同时,对高匝数线圈给探测系统带来的不利因素进行了分析和改进.最后,为了验证本研究所采用的方法,利用6 m边长的方形线圈和改进型MRS探测系统在中国温州一处在建隧道内的地下场地进行实地探测,探测结果与瞬变电磁方法结果对比表明,6 m边长的方形线圈可以对地下30 m内的含水层进行有效探测.

     

二维阵列线圈核磁共振地下水探测理论研究

核磁共振法(Magnetic Resonance Sounding, MRS)是一种直接探测地下水的地球物理方法,目前只能对水平层状的含水层进行一维测深,对于尺寸小于线圈直径的二维或三维含水构造成像时,其灵敏度和横向分辨率很低.本文从研究二维阵列线圈核磁共振地下水探测方式的可行性出发,推导了地面发射线圈产生的椭圆极化激发场和阵列接收线圈的核磁共振响应核函数的表达式,数值计算了二维阵列线圈测量方式的MRS信号初始振幅图像和探测灵敏度,通过与二维发射/接收同一线圈测量方式对比,得到结论如下:二维阵列线圈测量方式可直接确定地下含水构造的水平位置,但需更大的激发电流强度和接收灵敏度;当中浅层(2,以上三种分辨率均由仪器噪声水平所决定.通过实施增加接收线圈的数量和匝数、提高激发脉冲矩和接收灵敏度等改进措施,现有的仪器系统可用二维阵列线圈测量方式.本文提出的理论和方法,将促使二维阵列线圈核磁共振方法在地下暗河、喀斯特溶洞和堤坝渗漏等二维或三维含水构造快速探测方面得到应用.     

羌塘盆地冻土层结构核磁共振信号响应研究与实践

为了了解羌塘盆地冻土层的结构,利用核磁共振测深(MRS)技术对羌塘盆地的两个研究区进行了野外勘探。结合室内实验及野外实际资料,冻土层基本无法引起仪器核磁共振响应是划分冻土层和探测水合物的重要证据。在钻井附近进行的MRS方法对试验点进行测量并采集到了NMR信号,分析认为是天然气水合物处在亚稳定状态,分解产生水或者液态烃,形成水-水合物-气三相共存的情况,引起了NMR反应。实践表明可以将MRS方法运用于探测陆域冻土带分布以及冻土带,该方法对于冻土层的结构划分具有十分重要的科学意义和应用价值。     

核磁共振找水技术的研究现状与发展趋势

近20多年来,用核磁共振(Nuclear Magnetic Resonance,NMR)方法形成的一种直接非侵害性的探测地下水的地球物理新技术,与传统的地球物理探测地下水的方法相比具有高分辨力、高效率、信息量丰富和解的唯一性等优点,是一种很有发展前景的找水方法技术.我国的水资源短缺,对地下水资源的勘探、开发与利用十分重视,已将核磁共振找水技术研究列入国家十一五科技支撑计划.本文在广泛收集迄今为止的国内外大量资料的基础上,并根据作者近年来有关核磁共振找水技术的研究经历,综述了核磁共振找水技术的发展历史、现状和发展趋势,以推进我国核磁共振找水技术的发展.     

含水溶洞的地面核磁共振响应数值模拟研究

针对目前水平层状含水层的地面核磁共振响应特征已研究的比较透彻,而含水溶洞的响应特征却鲜有涉及的现状,本文给出含水溶洞地面核磁共振响应的计算方法.在文中,作者使用华军对双重Bessel函数无限积分的计算方法,实现了大回线源激发磁场的快速计算.再通过对地面核磁共振正演公式的离散化,实现了三维结构含水体的地面核磁共振正演计算...     

浅海底瞬变电磁法接收天线频率特性研究——带宽对浅层分辨的影响

运用电路理论建立接收线圈的等效电路;分析接收线圈的频率特性;探讨多匝小线圈的带宽对探测效果的影响,线圈的谐振频率较高,截止频率也相应较高,可以提高地下浅层的分辨率;设计一种高带宽多匝接收小线圈,通过对比高、低带宽接收线圈对同一断层进行的实测结果,阐明了接收线圈带宽在浅层瞬变电磁探测中的重要性.实验表明,高带宽的接收线圈会提高探测结果的浅层分辨率.     

基于参考线圈和变步长自适应的磁共振信号噪声压制方法

核磁共振地下水探测仪的灵敏度高,接收到的纳伏级磁共振探测信号极易受到强工频谐波噪声的干扰,导致信号特征参数提取的准确度降低,影响反演解释的水文地质参数结果.为了解决这一难题,基于相关抵消的原理,针对全波磁共振信号,设计带有参考线圈的90°移相自适应噪声抵消系统,理论计算了参考线圈相对于探测线圈的距离,提出变步长LMS算法进行噪声压制.仿真结果表明,在不同的信号强度及不同的信噪比下,当信号与工频谐波干扰频谱不重合时,采用设计的自适应噪声抵消系统和变步长算法,信噪比可以提高到5.94 dB以上,初始振幅、弛豫时间特征参数的拟合误差在2.8%以内;当信号与工频谐波干扰频谱重合时,采用双向自适应滤波算法,信噪比可以达到5dB以上,初始振幅、弛豫时间特征参数的拟合误差在10%以内,可以满足实际应用的要求;实测数据处理进一步证明了方法的有效性.     

基于参考线圈和变步长自适应的磁共振信号噪声压制方法

核磁共振地下水探测仪的灵敏度高,接收到的纳伏级磁共振探测信号极易受到强工频谐波噪声的干扰,导致信号特征参数提取的准确度降低,影响反演解释的水文地质参数结果.为了解决这一难题,基于相关抵消的原理,针对全波磁共振信号,设计带有参考线圈的90°移相自适应噪声抵消系统,理论计算了参考线圈相对于探测线圈的距离,提出变步长LMS算法进行噪声压制.仿真结果表明,在不同的信号强度及不同的信噪比下,当信号与工频谐波干扰频谱不重合时,采用设计的自适应噪声抵消系统和变步长算法,信噪比可以提高到5.94 dB以上,初始振幅、弛豫时间特征参数的拟合误差在2.8%以内;当信号与工频谐波干扰频谱重合时,采用双向自适应滤波算法,信噪比可以达到5dB以上,初始振幅、弛豫时间特征参数的拟合误差在10%以内,可以满足实际应用的要求;实测数据处理进一步证明了方法的有效性.     

基于椭圆极化的核磁共振找水理论研究(英文)

核磁共振找水是目前唯一能够直接探测地下水的地球物理方法,本文在weichman等人改进的核磁共振理论的基础上,应用经典的Chave算法对含有贝塞尔函数的积分核进行积分,求得了地下磁场和垂直激发场的空间分布。通过计算地下垂直激发场的椭圆极化率,验证了垂直激发场的椭圆极化效应,发现当地下电导率较大时椭圆极化效应将导致垂直激发场严重畸变;地下垂直激发场的正旋和逆旋分量显示了椭圆极化效应对核磁共振找水激发和接收上的不同影响,发射线圈和接收线圈间的延迟相位揭示了相位延迟效应的存在性及其相关性质。将以上理论应用到共圈模式下的核磁共振找水响应函数中,得到了地下氢核的扳倒角、地面核磁共振找水核函数和单一含水层模型的响应曲线,发现椭圆极化效应和相位延迟效应将导致核磁共振找水响应发生明显改变,因此对核磁共振找水理论和实际研究考虑椭圆极化效应是十分必要的。     

地下工程灾害水源的磁共振探测研究

地面磁共振探测(Surface Magnetic Resonance Sounding,SMRS)是近年来发展起来的一种地球物理新方法,这种在地面直接探测地下介质中氢核丰度的技术,不仅可以用于缺水地区的地下水资源勘查与评价,还可以在地下水引起的堤坝渗漏、滑坡、海水入侵等地质灾害水源的探测预警中发挥独特的作用.本文首次提出了地下磁共振探测方法(Underground Magnetic Resonance Sounding,UMRS),将SMRS方法引入到地下工程领域,实现隧道工程和煤矿开采等地下狭窄空间极端环境的探测.为应用UMRS方法,需要深入研究地下水超前探测理论、准全空间处理与反演方法、旋转多匝小线圈探测模式,强电磁干扰环境自适应噪声压制策略、以及复杂地质环境磁共振与瞬变电磁联合探测关键技术等难题.论文还简要介绍了超导磁探测技术和工程盾构及掘进实时探测等新技术在地下工程生产安全探测预警中的应用前景.     

基于双频磁共振的未知拉莫尔频率地下水探测方法

在利用地面磁共振方法进行地下水探测时,要求发射频率与拉莫尔（Larmor）频率相同.由于地磁场的不均匀性和噪声干扰,实际测量时无法获得准确的Larmor频率,导致发生偏共振现象.如果假设磁共振或频率偏量未知时,反演得到的含水量和弛豫时间（T₂^*）将存在较大误差.为了解决未知Larmor频率情况下的准确探测问题,本文在频率环测量方案的基础上提出了双频磁共振探测方法,即只需知道Larmor频率的范围区间,通过区间外的两个频率进行偏共振激发,利用频率差值的一半作为频率偏量计算核函数,即可消除未知Larmor频率的影响.通过仿真实验,证明了在电阻率大于100 Ωm时,双频磁共振和偏共振的核函数偏差较小.对于双层含水层模型,双频磁共振和偏共振的信号偏差小于3 nV,因此得到双频磁共振信号可以用偏共振核函数进行反演.通过对假设磁共振、假设偏共振和双频磁共振反演结果的对比,可以得到在未知频率偏量和设定频率偏量改变时,双频磁共振均能得到准确的反演结果.但是,当环境噪声增加和电阻率降低时,反演结果的准确度降低.最后,通过长春市烧锅镇采集数据的反演结果与已知钻探资料进行对比,验证了双频磁共振探测方法的有效性和准确性.

     

Numerical study of the variations of magnetic resonance signals caused by surface slope

When performing forward modelling and inversion of Magnetic Resonance Sounding (MRS) data, the water-content distribution is typically assumed to be horizontal (1D case). This assumption is fully justified because MRS is often used for characterizing continuous aquifers in a nearly flat environment. However, MRS can also be used in areas with sharp topographical variations. Following a review of the standard MRS equations when using a coincident transmitter/receiver loop, the mathematical terms potentially affected by tilting of the loop are discussed. We present the results of a numerical modelling exercise, studying a case where the surface is not horizontal and the loop cannot be considered to be parallel to the top of the aquifer. This shows that maximum variations in the MRS-signal amplitude are caused mainly by north- or south-dipping slopes. Slope effects depend on the loop size (a larger loop produces a larger error) especially in the presence of shallow water. With a geomagnetic-field inclination of 65° and a slope angle ≤ 10°, the topography causes a maximum variation in amplitude of less than 10%. Near magnetic poles and equator, the slope effect is lower and undetectable in most cases. It was found that within a 10% range of variation in the amplitude, errors introduced into inversions are within the typical uncertainty for MRS inversion and hence no topographic corrections are necessary. Thus, a significant effect from non-horizontal topography might be expected only when data uncertainty is lower than the slope effect (the slope effect is lower than equivalence when data quality is poor). Today, most field data sets are inverted using the modulus of the MRS signal, but some new developments consider the complex signal (both modulus and phase). However, inversion of complex MRS signals, which would provide a higher sensitivity to groundwater distribution, may be affected by slope effect. Thus, the slope orientation and dip angle should be accurately measured in the field when the phase of MRS signals is inverted too.     

Contribution of magnetic resonance sounding to aquifer characterization and recharge estimate in semiarid Niger

To improve the knowledge of the regionally important Continental Terminal 3 (CT3) aquifer in south-western Niger, fifteen magnetic resonance soundings (MRS) were carried out in December 2005 in the vicinity of wells and boreholes. The output MRS geophysical parameters, i.e. water content and decay constants versus depth, were compared to hydrogeological characteristics, i.e. water table depth, total porosity, specific yield and transmissivity estimated from direct measurements, pumping tests and transient groundwater modelling. The MRS-determined parameters were then used to estimate the rates of groundwater recharge.Contained in poorly consolidated Tertiary sandstones, the CT3 aquifer's water table has continuously risen by 4 m in total over the past four decades. Additionally, a significant portion of this increase has occurred in the past decade alone, with an annual rise now ranging between 0.1 and 0.3 m depending on the monitored well. Increase in groundwater recharge due to land clearance and deforestation explains this situation. According to previous estimations, the pre-clearing recharge ranged from 1 to 5 mm per year in 1950–60 s, while more recent recharge rates (1990s–2000s) range from 20 to 50 mm per year. These recharge values are directly affected by estimated aquifer specific yield value, while the spatial variation of rates of water table rise can be attributed to large scale hydrodynamic heterogeneities in the aquifer. However, few field measurements were available to confirm these assumptions.The main results of this study are: (1) The water table depth and aquifer transmissivity are estimated from MRS output parameters with an average accuracy of ± 10% and ± 9% respectively. (2) The MRS-determined water content is linked to both the total porosity and the specific yield of the aquifer, but no quantitative formulation can be proposed as yet. (3) Using the average MRS-determined water content over the investigated area, i.e. 13%, the groundwater recharge rates can be estimated to be ~ 2 mm per year in the 1950–1960s (pre-clearing period), and ~ 23 mm per year for the last decade. (4) The variations in specific yield and transmissivity cannot explain by themselves the spatial variability of the rise of the water table. (5) The ranges in transmissivity and water content obtained from MRS are more realistic than the groundwater modelling outputs. Therefore, MRS could be used to better constrain the aquifer parameters in groundwater modelling with a dense site network.Finally, this work illustrates how MRS can successfully improve characterisation and transient multi-year groundwater balance of commonly found sedimentary aquifers, particularly when integrated with well observations and pumping tests.     

均匀介质中MRS方法三维模型的核磁共振响应

MRS方法三维含水体的核磁共振正演计算涉及两个难点:①激发场表达式的快速计算,该表达式为包含双重Bessel函数乘积的积分式,将该表达式分为两个积分区间,第一积分区间应用贝塞尔函数的汉克尔函数表述式及其大宗量渐近特性,将其转化为Fourier正(余)弦变换,采用快速计算方法实现数值计算;另一积分区间的计算直接应用一般的数值积分算法即能获得较高的计算精度。②三维含水体的空间离散化,提供柱坐标和直角坐标中的两种离散方法,可对任意形状的三维含水体进行构建。在此基础上,开发出相应的正演程序,并进行几种常见的三维含水体的MRS正演,对其初始振幅曲线进行规律总结。结果表明,三维含水体的MRS振幅响应与层状含水层振幅响应类似,使用收发共圈的工作模式将难以辨别地下含水体的空间形态。     

基于谐波建模和自相关的磁共振信号消噪与提取方法研究

磁共振探测技术（Magnetic Resonance Sounding,MRS）以其无损、定量、直接等优势,被广泛应用于地下水调查、水文环境评价以及灾害水源预警等领域.在实际应用中,强工频谐波和随机噪声等严重影响MRS信号的质量,导致后续水文参数解释不准确.针对这一问题,提出谐波建模和自相关相结合的方法进行消噪以及信号特征参数提取.首先构建工频谐波模型,针对建模算法严重依赖工频基频精度的问题,采用自适应扫描方式搜索方案,大幅提高搜索准确度和速度;其次推导了MRS信号自相关表达式,提出了自相关参数提取的非线性拟合方法.仿真数据结果表明,建模消噪方法有效消除了工频谐波,信噪比平均提升了17.03 dB;自相关处理后,信噪比进一步提升了16.10 dB,初始振幅和弛豫时间参数提取结果的准确度比处理前分别提高了3.8倍和2.8倍.通过不同信噪比和弛豫时间的重复实验,得到当噪声水平小于200 nV和弛豫时间大于200 ms时,自相关参数提取具有较高的稳定性.最后,通过野外实测数据处理实验,进一步验证了联合消噪和参数提取方法的有效性.

     

基于双通道参考技术的磁共振地下水探测电磁噪声分布规律研究

磁共振地下水探测（MRS）技术能直接测量水中氢质子的拉莫尔进动,可表征地下水含水量和介质孔隙度等信息,因而获得了广泛应用.然而,MRS信号微弱,实际探测中电磁噪声强度大,若噪声未被有效去除,会导致后续水文参数解释不准确,采用"8"字形线圈或运用带参考线圈的噪声抵消等方法可抑制噪声干扰,但这些方法消噪效果的好坏取决于电磁噪声的空间分布.另外,在强噪声或多噪声源环境下应用MRS方法,选择恰当的探测位置对提高探测效率以及增加探测结果的有效性具有重要意义.因此,本文使用双通道噪声采集装置在实验前记录环境中的电磁噪声,通过双通道参考技术推导同一时刻不同测点的噪声强度,再得到测区电磁噪声分布结果,该噪声分布情况对MRS方法探测方式的选取具有指导性的作用.同时,可根据噪声分布情况为探测线圈选择最佳测点以避开噪声强度大的区域,进而提高采集信号的信噪比.最后,通过采集室外环境下的电磁噪声,验证本文提出的基于双通道参考技术的电磁噪声分布规律研究方法的有效性,为可靠高效地开展磁共振地下水探测提供了新的技术手段.     

2D magnetic resonance tomography applied to karstic conduit imaging

Karstic conduits play a crucial role for water supply in many parts of the world. However, the imaging of such targets is generally a difficult task for most geophysical methods. Magnetic Resonance Sounding (MRS) is a geophysical method designed for imaging of water bearing structures. Initially, MRS was developed for characterizing horizontally stratified aquifers. However, when applying a 1D MRS measuring setup to the imaging of 2D–3D targets, the size of which may be much smaller than the loop, the accuracy and the lateral resolution may not be sufficient. We have studied the possibility of simultaneously processing several MRS aligned along a profile to perform a Magnetic Resonance Tomography (MRT). This work emphasizes the gain of resolution for 2D–3D imagery of MRT versus the interpolation of 1D inversion results of MRS along the same profile. Numerical modelling results show that the MRT response is sensitive to the size and location of the 2D target in the subsurface. Sensitivity studies reveal that by using the coincident transmitting/receiving (TX/RX) setup and shifting the loop around the anomaly area, the depth, section and position of a single karstic conduit with a size smaller than the MRS loop size can be resolved. The accuracy of the results depends on the noise level and signal level, the latter parameter being linked to the depth and volume of the karstic conduit and the water content in the limestone matrix. It was shown that when applying MRT to the localization of 2D anomalies such as karstic conduits, the inclination of the geomagnetic field, the orientation of the MRT profile and the angle of crossover of the conduit by the MRT profile must be taken into account. Otherwise additional errors in interpretation should be expected. A 2D inversion scheme was developed and tested. Both numerical and experimental results confirm the efficiency of the developed approach.     

Direct and indirect methods for groundwater investigations: A case-study of MRS and VES in the southern part of Sweden

Magnetic Resonance Sounding (MRS) has been successfully tested for detecting groundwater in two areas in southern Sweden. Measurements of Schlumberger VES have been conducted in the same place as the MRS and the results are generally consistent. Low resistivity layers interpreted as clay are sometimes identified close to the surface. The MRS result in site 2 is a good example of signals penetrating through the clay and deeper aquifer still being detected. The MRS data suggest aquifers that are not only hosted in soft sediment materials (moraine, sand, and mixed materials), but hosted in basement rocks. Based on the MRS and borehole pumping test data, the results agree with yield, average water content and subsurface geological data.