基于球坐标系下有限差分的地磁测深三维正演

为了计算全球尺度电磁感应的响应,本文介绍地磁测深频率域三维正演。正演算法采用球坐标系下的交错网格有限差分方法,从Maxwell方程的积分形式出发,采用PARDISO对离散后的方程组求解,避免了迭代求解的散度校正。为了验证本文结果的正确性和精度,与前人的有限元和有限差分方法进行了对比,一维层状模型的三维交错网格有限差分数值结果和解析解相对误差小于5%,双半球模型的计算结果与前人的计算结果完全吻合。三维"棋盘模型"计算表明磁场分量对异常体的大小和位置具有很好的分辨能力。     

利用归一化总梯度方法实现中梯激电剖面数据的解释

将归一化总梯度数据处理技术运用于中间梯度剖面视极化率观测数据的延拓处理,研制了一套实用的中梯激电数据解释软件系统.实验数据与野外数据的处理结果表明,在选择合理的计算参数基础上,归一化总梯度能给出激电异常较为可靠的解释.     

三维电性结构揭示的中国东北地区新生代火山深部起源

利用上地幔尺度的三维电性结构探讨了中国东北地区新生代火山的深部起源及彼此间联系.文章在该地区采集并获得了高质量的长周期大地电磁全阻抗张量数据,并采用高斯-牛顿法进行三维反演,得到了电阻率模型.该电阻率模型显示,中国东北地区高阻岩石圈内,存在多个局部低阻异常,它们在空间上和该区新生代火山对应;在上地幔深部,环松辽盆地分布三个低阻异常;软流圈中存在着脉状低阻异常,联系着岩石圈内和上地幔深部的低阻异常.岩石物理分析表明,这些低阻异常由流体熔融物质引起.俯冲到中国东北地区地幔转换带中的西太平洋板块发生持续脱水作用,模型则进一步揭示,这些水向上迁移到松辽盆地下方的转换带顶部并引起上地幔底部的部分熔融,在浮力作用下形成热物质上升流,融透软流圈后底辟于干枯岩石圈底部,最后刺穿较为薄弱或者活化的缝合带等构造边界,排出地表形成包括长白山火山、哈拉哈火山等整个中国东北地区新生代大规模火山事件.因此认为,中国东北地区的新生代火山"本是同根生",它们均起源于上地幔深部的环状低阻异常.     

基于可控源电磁法阻抗信息的有限内存拟牛顿法三维反演

本文研究了利用阻抗信息进行可控源电磁勘探有限内存拟牛顿法三维反演的技术。首先用理论模型来验证有限内存拟牛顿法反演的准确性和可行性。观测参数为复阻抗Z_xy分量,采用交错网格有限差分方法计算模型响应,反演采用有限内存拟牛顿法。数值模拟结果表明：有限内存拟牛顿法反演迭代速度较快,每4 min迭代一次,拟合差由146.00下降到1.78,收敛稳定;异常体的位置与理论模型吻合较好,有效地验证了有限内存拟牛顿法可控源三维反演的正确性。为了进一步验证该方法的实用性,将其应用到隐伏钼矿可控源电磁勘探工作中。工区的反演结果显示：在工区北西段深部存在高阻异常,其上为低阻异常。截取过钻孔的3号测线发现,-500~-100 m的位置表现为低阻,东侧存在向上涌起的高阻。推测此低阻为矿化蚀变带,延伸较深。该异常与钻探资料揭示的钼矿脉一致,证明了反演结果的准确性。因此,利用可控源观测得到的阻抗信息进行有限内存拟牛顿法三维反演,可以获得可靠的三维电阻率分布。     

音频大地电磁数据非线性共轭梯度3D反演在江西相山铀矿勘查中的应用

相山铀矿田是中国重要的铀矿生产基地。为研究相山铀矿田地球物理特征,进而寻找相山地区铀成矿的有利部位,笔者利用非线性共轭梯度3D反演方法进行反演计算。结果显示,该地区地下结构表现出明显的3层结构,推测第一层高阻异常是下白垩统鹅湖岭组,第二层低阻异常是下白垩统打鼓顶组,第三层高阻异常是中元古界基底变质岩。结合地质情况,根据测井曲线得到第一层与第二层之间的第一成矿界面电阻率约为2 000Ω·m,第二层与第三层之间的第二成矿界面电阻率约为500Ω·m。依据两成矿界面电阻率描绘出测区的成矿界面的形态特征呈现地区分布不均匀性,第一成矿界面西南部保存较好,测区南部被断层切割明显;第二成矿界面在中部保存较好,边缘缺失。根据成矿界面和断裂复合部位得到成矿有利区,第一成矿界面的有利区位于测区中部和南部,第二成矿界面的有利区位于测区边缘。     

油水双相饱和孔隙模型核磁特性理论研究

利用特征函数展开方法，给出了基于扩散效应的核磁共振Bloch控制方程的弛豫模式解表达式, 在此基础上分别给出了水饱和与油水双相饱和孔隙模型的核磁共振特性模拟结果. 结果表明，孔隙流体的核磁弛豫与孔隙大小、孔隙表面弛豫率、孔隙流体的扩散能力以及含油饱和度等有密切的关系. 对于水饱和孔隙，弛豫时间主要由孔隙大小控制. 当孔隙较大时，即使均匀大小孔隙，孔隙流体的弛豫也会表现为多指数弛豫. 而且最小模式弛豫时间与孔隙大小为非线性关系. 对于油水双相饱和孔隙，在孔隙较小时，含油饱和度对弛豫的影响主要表现在弛豫时间随含油饱和度的增加而线性减小；但在孔隙较大时，含油饱和度的增加对弛豫影响表现在两个方面，其一，孔隙水弛豫由多个弛豫模式控制逐渐转变为由最小弛豫模式控制；其二，孔隙水弛豫时间与含油饱和度表现为非线性关系. 对由实际岩芯抽象出的孔隙模型，采用本文获得的理论公式，在油水双相饱和时进行了正反演模拟. 计算结果与已有的实验结果较为一致.     

综合物探方法在应急抗旱找水中的应用

通过结合区域地质资料,采用重力、高密度电法、激电测深及音频大地电磁测深等多种地球物理勘探技术,计算机数值模拟反演计算等工作手段,为河南省鹤壁市气象干旱缺水区开展应急抗旱找水。多种物探方法的运用,对于应急抗旱找水打井的井位合理的布置、快速定井,保证成井后的出水率和地下水的合理利用提供了重要的依据,为综合物探方法的应用和类似应急找水工作提供了一定的借鉴。     

基于姿态参数的海洋可控源电磁数据一维反演

基于姿态信息的海洋可控源数据一维反演采用Occam方法进行,反演数据为不同收发距离观测到的电场振幅(MVO)。接收器和发射器姿态参数同时引入到反演中,正演过程引入发射源姿态参数,借助矢量分解将任意姿态的发射源分解为3个正交坐标轴分量源,分别计算接收点处电场响应,叠加得到总电场响应;根据接收器姿态信息,将计算坐标系中的电场矢量进行欧拉旋转获取任意姿态时测量到的电场响应。理论模型与实际数据反演结果表明,基于姿态信息进行海洋可控源勘探数据反演是可行的。     

层状导电介质中地面核磁共振响应特征理论研究

水平层状导电大地中层状含水层地面核磁共振响应的数值计算技术涉及到导电介质中回线源 磁场的计算以及地下含水层中质子磁矩在线圈中产生感应信号的体积积分. 文中采用直接数 值积分技术,对具有振荡核函数的Hankel变换进行积分,以求取回线源磁场的径向与垂向分 量,计算并研究了回线产生的径向与垂向磁场分量随空间位置的变化规律. 基于磁场空间分 布特点,利用不等间距空间剖分技术计算地面核磁共振的体积积分,模拟了不同模型的地面 核磁共振响应并讨论了其影响因素. 结果表明,结合能对任意层状导电介质中磁场进行稳定 快速计算的直接数值积分技术与不等距空间模型剖分技术,可正确模拟地面核磁共振响应. 导电性是产生地面核磁共振信号相位的先决条件,但影响响应振幅强度与相位的因素还有含 水层的埋深、厚度以及装置大小等.