基于MPI并行算法的电阻率法多种装置数据的三维联合反演

由于地表电阻率法受到浅地表局部异常体的干扰，反演精度受到影响，井中装置数据资料参与反演虽然可以减小浅地表局部异常体的干扰，但是由于钻井位置的局限性，数据量得不到保障，也会导致反演精度降低.为此，本文开发了一套结合地表、地-井、井-地和井-井装置数据的三维联合反演算法.首先，利用有限差分法实现正演模拟，采用非线性共轭梯度法（NLCG）恢复电阻率结构；其次，调用Message Passing Interface（MPI）函数库加速正演模拟和灵敏度矩阵运算，当开辟12个进程时，反演程序获得最大加速比4.51；最后，通过合成数据和实测数据算例证明该反演算法的有效性，也证实了多种装置组合数据体反演结果明显优于单一地表装置数据体反演结果.     

三维张量CSAMT响应特征研究

张量CSAMT(可控源音频大地电磁法)能完整地描述三维的地电信息,更适应复杂地质条件下的地质勘探。采用将总场分解为一次场和二次场的策略,用差分近似微分形成大型线性方程组,通过解方程组得到张量CSAMT三维正演结果。通过计算低阻和高阻目标体的响应,发现卡尼亚视电阻率和视阻抗相位对低阻体表现为低值异常,对高阻体表现为高值异常,视阻抗相位达到极值时的频率比视电阻率达到极值时的频率要高。张量CSAMT对低阻体比高阻体更灵敏,在同样条件下,高阻体的异常区位置相比于低阻体会上移。某工区的勘察案例表明,根据卡尼亚视电阻率异常能定性地判断地下异常体的存在,为解决地质类问题提供了依据。本研究为张量CSAMT三维反演提供了基础。     

铜陵矿集区深部找矿工作的一些思考

铜陵矿集区深部结构特征为该区大规模成岩成矿作用奠定了良好的基础,浅表已发现的矿床只是矿集区成矿系统的一部分。除此之外,深部还应有多个可以找寻的矿床(矿体)的富集空间。与成矿相关的花岗岩研究表明,铜陵矿集区可能存在与宁芜火山盆地同期的成岩成矿作用,是铜陵矿集区找矿的新方向。铜陵矿集区奥陶系灰岩中铜金矿化作用的初步研究表明,除传统矿床类型外,前泥盆纪地层中的成矿作用应该在今后的找矿工作中得到关注。覆盖区成矿预测初步研究表明,铜陵矿集区向斜核部及褶皱倾伏端等地段具有与已知典型矿床相似的成矿条件。铜陵矿集区的进一步找矿工作应该在上述地段运用矿田构造系统分析、深部地球物理探测和工程验证等方法逐步展开。     

复杂地区油气地球物理勘探技术集成

随着我国油气勘探程度的不断提高,勘探难度不断加大.全国油气资源战略选区调查与评价项目历经六年,通过对深水海域、西部复杂山地、西藏高原、南方碳酸盐岩、火山岩覆盖区等几类典型地球物理勘探久攻不克的复杂地区开展地震、重磁电、综合地球物理勘探联合攻关,以及天然地震层析成像攻关实验,取得了长足进步,直接带动获得了一批有价值的油气...     

倾子资料三维共轭梯度反演研究

在对倾子响应和共轭梯度算法深入分析的基础上,我们实现了倾子资料三维共轭梯度反演算法.基于倾子资料的三维共轭梯度反演研究,探讨了利用倾子资料进行三维反演定量解释的方法.通过对理论模型合成数据 进行反演试算,验证了所实现的倾子资料三维共轭梯度反演算法的有效性和稳定性.该反演算法可用于对大地电磁测深和地磁测深(地震地磁台站进...     

利用远震接收函数揭示的喜马拉雅东构造结台阵下方地壳结构及其动力学意义

在喜马拉雅造山带的东缘,雅鲁藏布江缝合带在这里发生急剧转折,南迦巴瓦变质体快速隆起,然而关于东构造结的形成机制一直未有定论.利用围绕南迦巴瓦峰的48个宽频带地震台站记录的远震数据提取P波接收函数,采用改进的H-κ叠加方法和共转换点叠加方法综合研究了东构造结的地壳厚度、波速比分布和地壳结构特征.结果表明:研究区平均地壳厚度为64.03 km,大部分台站介于60.48~66.55 km范围;平均波速比为1.728,主要集中范围为1.696~1.742.东构造结地壳厚度横向变化剧烈,构造结西端和北端厚而中间薄,东构造结核部Moho面呈现上隆的构造形态,东西向上隆幅度约为6~7 km,南北向的上隆超过9~10 km.东构造结核部地壳上隆减薄可能由高密度、高波速的岩石圈撕裂残片拆沉到上地幔软流圈后重力失衡所致.平均波速比超过1.8的高值异常展布于东构造结的两侧,推测为环东构造结的壳内部分熔融体.东构造结地壳上隆减薄和壳内部分熔融的存在很可能均与幔源热物质的上涌有关,而软流圈地幔的上涌则可能由印度板片的撕裂引起.     

班公—怒江构造带的电性结构特征——大地电磁探测结果

对青藏高原过班公—怒江构造带的三条大地电磁剖面进行探测,获得班公—怒江构造带及其邻区的电性结构模型,研究了班公—怒江构造带的深部结构与构造特征.研究结果表明：构造带及其两侧上地壳内广泛分布不连续高阻体,反映了岩浆岩的空间分布特征,表明构造带南北两侧岩浆的活动规律可能存在较大差别.研究区内的冈底斯及羌塘地体的中、下地壳普遍发育高导层,反映了印度大陆碰撞、俯冲过程的效应与痕迹,而高导层之下的高阻块体则可能是向北俯冲、冷的、刚性的印度大陆地壳.羌塘地体的电性结构模型可以分为南北两个区段,南羌塘块体的壳内高导层与班公—怒江构造带对印度板块俯冲的阻挡作用有关;而北羌塘块体壳内高导层与亚洲大陆对印度板块向北俯冲的“阻挡”与向南“对冲”有关.印度板块向北的俯冲与挤入,受到班公—怒江构造带及亚洲板块的阻挡,可能没有越过班公—怒江构造带,并在班公—怒江构造带附近向下插入软流圈,导致幔源物质上涌,形成壳、幔热交换与物质交换的通道和规模巨大、延伸至上地幔的高导体.班公—怒江构造带的电性结构证明了该构造带是一组产状陡立、巨型的超壳深断裂带.     

华北中部岩石圈电性结构——应县—商河剖面大地电磁测深研究

20 0 1年, 沿着山西应县到山东商河, 重新布置大地电磁测深剖面进行研究.采用现代先进的大地电磁数据处理技术和快速松弛二维反演方法获得该剖面二维电性结构模型, 从而充分展示了华北地区岩石圈电性结构的特点.从电性特征上讲, 华北岩石圈以太行山前断裂为界划分为东、西两区, 东区为低阻区, 西区为高阻区.在东区, 上地壳电性结构基本与华北裂谷系的隆、坳构造格局相对应, 岩石圈的电导最高达3× 104 S, 远远大于强烈活动的安第斯山岩浆弧区和西藏高原岩石圈的电导.这里, 在构造连接部位的地壳中有不连续的高导体存在, 电导率大约0.1~ 0.8S/m.在西区, 太行山和恒山的岩石圈为高阻块体, 表现出稳定大陆区岩石圈导电性结构的特点.但恒山高阻块体之下发现一组向西缓倾的高导层, 其电导率为0.0 4~0.2 5S/m, 顶面在2 0km深处, 底面深度大约40km.      

被动源激发极化法的野外试验结果

对被动源激发极化法场源的基本性质、激发极化效应的提取方法进行了讨论,以野外试验结果为例研究了矿体与异常的对应规律,对该方法的有效性和可行性进行了说明。     

西藏泽当岩体含铬超镁铁岩体的深部构造特征与找矿前景探讨

在中国最大的铬铁矿罗布莎周边地区有着数目众多的含铬超镁铁岩体;分布广阔,具有很好的找矿前景,选择泽当岩体进行了重、磁、大地电磁法的剖面探测,结果表明:泽当岩体被近南北向断裂切割成多块,其西端岩体向北推移的距离远于东段的岩体,总体上呈北西西向展布;岩体磁异常由蛇纹石化产生的磁铁矿和岩体生成时剩余磁化引起,磁异常指出了岩体浅部的范围;航磁较大范围的异常表明存在较深部的磁性体;重力异常说明岩体局部存在高密度体,且与磁异常相吻合,地表出露的矿点多分布在此范围内;岩体表层为高阻薄层,在剖面中段高阻体厚度最大处达1.5 km,范围不大.在薄高阻体下为薄层低阻层.电阻率剖面两侧以高阻体为主,中段有多个低阻体,且与深部5～25km深度上的低阻体相连,中地壳范围的低阻体与幔源物质上涌在地壳范围留下的熔融体有关,与地壳上部的超镁铁岩体有关.综合各种地球物理场的特征认为在磁异常分布区内具有局部重力异常、低阻体异常较多的地段是泽当岩体最有利于进一步详细探测找矿的部位.