采空区的物探勘查方法

采空区给人类生产生活带来危害,采空区的探测已成为重要研究课题.文章介绍了采空区的地球物理特征,对采空区探测的各种物探方法从重、电、震、放、综合物探等方面进行分类和介绍,有以往使用的老方法,也介绍了采空区探测中使用的新方法.     

环境与工程地球物理的发展趋势

2006年6月在武汉召开了第二届环境与工程地球物理国际会议.会议围绕近年来兴起的近地表探测问题展开,涉及方法原理、仪器采集、资料处理等各个环节,反映出国内外环境与工程探测领域的若干前沿热点和发展动态.结合本次会议,可以看出该学科主要趋势为：应用领域不断扩宽,方法技术灵活多样,重视基础理论研究,开发新型（特种）仪器,通过高精度测量、延时测量、先进的数据处理技术和信息技术等手段改善探测质量.我国对环境与工程地球物理有很大的需要,应根据我国实际情况,吸收国外先进经验,积极开展自主创新.     

皖南云岭金矿地球物理找矿模型

皖南云岭金矿位于长江中下游多金属成矿带与华南成矿带之间,属于江南造山带上的一个热液脉型矿床。基于云岭金矿区域地质背景,在该区开展了岩石露头测定、岩石标本物性测定和地球物理学相关研究。运用地面高精度磁测法,识别本区主要控矿构造为NE向断裂带。根据激发极化法圈定高极化率异常区,并结合相应区域电阻率(可控源音频大地电磁法)和磁异常分布特征建立4种地球物理模型:低阻-高极化模型、高阻-高极化模型、高极化-磁异常过渡带模型和高极化-高磁异常模型。这4种地球物理模型建立在该区岩石物性与矿床成因基础之上,分别对应不同含矿模型:低阻-高极化模型对应接触带上硫铁矿成矿区;高阻-高极化模型对应含矿硅化带;高极化-磁异常过渡带模型对应围岩蚀变带;高极化-高磁异常模型对应岩浆侵入区。依据以上模型,可圈定NE向断裂带、断裂带西侧高阻硅化带及东侧蚀变岩体为找矿有利靶区。     

内蒙狼山成矿带地球物理场特征

内蒙狼山区域深部构造非常发育、复杂,岩浆岩具有多期侵入且明显受到深部断裂构造控制的分布特征.目前对于该地区深部断裂构造、区域地质结构特征、成矿远景区的划分,还没有详细深入的研究.结合已有的地质资料,对狼山成矿带主要构造区进行了地球物理场分析,研究了狼山成矿带6条断裂构造、岩浆岩分布及区域深部地质结构特征,圈定了4处成矿远景区,为深入认识狼山成矿带成矿环境地质结构特征和成矿潜力提供了依据.      

滇西龙陵地区地壳电性结构及其对大瑞铁路地质选线影响研究

高黎贡山地处印度板块与欧亚板块碰撞缝合带附近的横断山脉南段,是大理一瑞丽铁路(大瑞线)的必经之地,地形起伏大、构造复杂、活动性强,高黎贡山隧道作为全线控制性工程之一,其地质选线的最大困难就是对隧道深部构造环境的了解,特别是缺少对与地热、地震等联系紧密的深部地质构造的认识.为此,本文以大地电磁方法为手段,以高黎贡山隧道为主要研究对象,通过对滇西龙陵地区高黎贡山隧道越岭段两条大地电磁剖面数据的处理解释对研究区的地壳电性结构特征进行了勘探研究.结合区域地质构造特征与主要工程地质问题之间关系的分析,根据隧道主要断层地质条件设计了三维垂直断层模型,利用三维有限元开展正演模拟研究发现,测点点距、位置与横向分辨率密切相关,点距越密,分辨率越高,测点位于断层在地表投影位置能有效提高分辨率.采用大地电磁阻抗张量分解技术对两条剖面上各测点的二维偏离度和电性走向进行了计算和分析,对剖面视电阻率和阻抗相位数据进行了二维NLCG联合反演研究,揭示了沿剖面的腾冲地块、龙陵一瑞丽断裂带及保山地块10 km深度的电性结构特征及相互关系.结果表明:剖面CD电性结构呈现区域构造的三分性,腾冲地块电性结构成层性较好,保山地块成层性较差,两者均以中高阻电性特征为主,中间夹龙陵一瑞丽断裂带,电性结构反映从3 km深度以下存在几乎近于直立延伸的低阻带,推测为班公湖一怒江缝合带滇西段丁青一怒江缝合带的反映;剖面AB共划分了6条与工程密切相关的深部隐伏断裂,结合地震地质、地表地质及龙陵地震深部背景研究,推测F7-3断裂为1975龙陵7.3级地震断裂;从地表黄草坝断裂开始向下延伸,有一条发育最大深度约为4 km的低阻通道,推测为地热断裂深循环通道,其与黄草坝断裂共同控制研究区地下热水的补给、径流和排泄条件,在高黎贡山隧道线位位置形成了一个相对低温通道,为隧道方案成立的关键工程地质条件.勘探结果表明:滇西龙陵地区地壳电性结构有效的反映了高黎贡山隧道深部隐伏断裂和地热断裂深循环通道等深部构造特征,为大瑞线隧道工程地质选线提供了深部地质背景依据.     

青藏高原岩石圈三维电性结构

本文报道通过综合大地电磁调查数据研究青藏高原岩石圈三维电阻率模型的初步成果.大地电磁法调查区域已经覆盖了高原大部分面积,为全区三维电阻率成像研究打下了可靠的基础.对多个测区大地电磁数据进行精细的同化处理和反演成像,取得了青藏高原可靠的岩石圈三维电阻率结构图像.成像的区域为28°N—35°N,80°E—104°E.三维反演计算时采用的网格尺寸为20km×20km,垂直方向不等间距剖分为26层.结果表明,青藏高原现今岩石圈电阻率扰动主要反映印度克拉通对亚欧大陆板块俯冲引起的热流体运动和大陆碰撞和拆离产生的构造.在岩石圈地幔,察隅地块、喜马拉雅地块和拉萨地块东部联成统一的高电阻率地块,它们反映了向北东俯冲的印度克拉通.雅鲁藏布江、班公—怒江和金沙江缝合带都有明显的低电阻率异常,表明岩石圈深处有热流体活动.雅鲁藏布江、班公—怒江和金沙江缝合带都有明显的低电阻率异常,也表明它们的岩石圈还有流体活动.青藏高原东部的低阻区沿100°E向地幔下方扩大,反映了金沙江断裂带有切穿岩石圈的趋势.地幔电阻率平面扰动的模式显示,青藏高原东西部的地体碰撞拼合形式和方向是不同的.在青藏高原西部,羌塘、拉萨和喜马拉雅等地体从北到南碰撞拼合.在青藏高原东部,羌塘—拉萨、察隅、印支、雅安和扬子等地体多方向拆离拼合,在地壳造成不正交的拆离带和压扭构造系.从高阻-低阻区的分布看,东部的地体拼合有地幔的根源,今后还会进一步发展.察隅地块岩石圈对青藏高原东部的楔入,使其北部和东部地块的岩石圈发生拆离撕裂,也造成热流体上涌的低电阻率异常.     

青藏高原东北部多尺度重力场及其地球动力学意义

由于受到SN和EW双向挤压作用的综合影响,处于年轻高原与古老地块交接区的青藏高原东北部岩石圈强烈变形,构造活动十分活跃.为了整体性地了解青藏高原东北部重力场和深部动力学机制,本文基于EGM2008全球重力场模型数据,利用小波多尺度分析获得不同尺度的重力异常信息;同时反演了研究区的地壳厚度;通过构建穿越龙门山造山带和西秦岭造山带两条剖面的岩石圈密度结构模型,分析了地壳上地幔内不同介质的分布特征.研究结果表明,青藏高原东北部岩石圈显示十分复杂的塑性体的特征,其重力异常走向多以EW或SSE为主,反映了高原岩石圈物质向东运移的趋势;地壳厚度由西向东逐渐减薄,边缘造山带深部并没有发现"山根"痕迹,结合该地区低重力的特征,推测西秦岭—松潘构造结岩石圈发生过大规模地幔流底侵作用;地幔流上涌的动力可能来源于印度板块向欧亚大陆板块俯冲,激发了地幔流体侧向移动,在扬子地台和华北地台附近受到坚硬岩石圈的阻挡而被迫上移,并因此造成龙门山与西秦岭的隆升.     

基于MPI的一维大地电磁并行计算研究

在大地电磁中运用并行处理技术来减少计算时间,提高运算效率.结合一维广义逆矩阵法反演的计算特点,详细分析了串行程序并行化方法,主要采用主从并行模式、分频并行计算的并行方案.在Linux系统上使用Fortran和MPICH2相结合的开发工具编写了并行程序,通过理论模型和实测数据对实现的并行程序进行试算,并与串行程序对比验证了该算法的可行性、正确性,为二维、三维的正反演并行计算提供了研究基础.     

导线环垂直激发场高精度数值模拟

导线环垂直激发场是导线环在地下产生的激发场中垂直于地磁场的分量,在地面核磁共振找水方法中,垂直激发场在地下的分布直接决定地面核磁共振响应,探讨影响垂直激发场地下分布的诸因素进而研究这些因素对地面核磁共振响应的作用,对核磁共振找水正确野外施工具有一定的指导意义,有助于采集更高质量的野外数据.利用经典Chave算法对含有贝塞尔函数的强振荡积分核进行了数值积分试验,将其计算结果和对应的解析解进行了对比,试验结果显示该算法用于强震荡函数无穷积分的计算具有良好的收敛性、较高的精度和较快的计算速度,采用Chave算法对含有贝塞尔函数的强振荡函数进行积分是可行的.在此基础上,进一步计算并数值模拟了导线环垂直激发场及其相位,分别从地磁倾角、磁偏角、地下介质导电性、深度等方面探讨了影响垂直激发场及其相位分布的主要因素.