预优共轭梯度法及井资料约束的磁化强度成像

磁化强度成像反演方法可以得到地下介质的磁化强度分布,可以反演任意形状磁性体,比模型参数反演占优势。但因为该方法计算量和存储量较大,分辨率低,"趋肤效应"严重,导致其应用受到限制。共轭梯度法迭代计算的基本单位是列向量,减少了计算时间和存储空间。再将核矩阵作用于与单元体深度有关的预优矩阵,可以提高反演分辨率并减弱"趋肤效应",并且对于有钻孔的区域,进行井资料约束,克服了反演问题的多解性。该理论模型模拟及青海省尕林格矿区地面高精度磁测磁化强度成像反演说明该方法具有良好的应用效果。

磁化强度矢量反演方程及二维模型正反演研究

推导了磁化强度矢量层析成像方程 ,并与磁化强度标量层析成像方程进行了对比 .使该矢量层析成像方程既适用于三维的也适用于二维 ,既适用于使用磁场垂直分量资料 ,也适用于使用磁场总强度资料 .本文采用改进的高斯 -赛德尔迭代求解磁化强度矢量层析成像方程 ,在求解方程中引进了与深度有关的权系数 .并对二维模型开展了正反演研究 .当模型层数为二层 ,且每层 5 1个柱体时 ,采用零初始模型就能获得较好的磁化强度垂直分量和水平分量反演结果 ,除了异常体边部外 ,磁化方向比较可靠 .当模型层数为五层 ,且每层 5 1个柱体时 ,采用零初始模型不能得到较好结果 ,改用接近背景值的均匀初始模型 ,反演的磁化强度垂直分量和水平分量大致接近真实 。

二维磁异常的粒子群快速成像方法及其应用

通过将地下介质划分为一定规模的网格单元,并构造一个与磁性体中心网格位置有关的指数函数模拟地下介质磁性分布特征,利用粒子群算法反演该函数参数,根据求得的函数表达式,推算各网格单元磁化强度,由此实现磁化强度的快速成像。理论模型试验显示:该方法不仅能将一个求解20×10网格规模的物性分布问题简化为求解指数函数的10个参数值的问题;而且,当在不同模型地质体磁异常中加入20%的随机噪声时,也能取得理想的效果,且需要调整的参数少,算法通用性强。利用该算法处理老挝爬立山铁矿区某磁测剖面,得到厚度近140m的强磁性体;钻孔证实,实际钻遇矿体累积最大厚度达130.6m,且矿体中心与成像中心相对应:由此证明该方法应用于实际资料处理与解释的可行性和有效性。

根据磁异常确定磁性体上顶埋深的几种反演方法——以尕林格矿区磁测资料解释为例

根据磁异常确定磁性体上顶埋深是磁法勘探中寻找磁性物体及磁性矿床的重要问题之一。笔者总结了4种用磁异常确定磁性体上顶埋深的反演方法,它们是切线法、欧拉法、功率谱法、磁化强度成像法。本文首先将这4种方法用于板状体模型的反演,反演深度与板状体实际深度一致。然后,将它们应用于青海省尕林格矿区地面高精度磁测资料的反演解释,以勘查被广厚的第四纪砂砾沉积覆盖的磁铁矿,取得较好找矿效果。