海洋可控源电磁法二维有限元正演及其反演.pdf

谨以此论文献给我的老师、同学、亲人和朋友们 刘颖 万方数据 万方数据 万方数据 海洋可控源电磁法二维有限元正演及反演 摘 要 进入 21 世纪，石油和天然气勘探的重点已由陆地向海洋转移。三维海洋地 震勘探一直是为海洋钻探提供海底地质构造特征的重要手段，但石油公司面对海 洋油气钻探的巨额风险，并不满足于仅获知海底地质构造的特征，还希望在钻探 前得知构造内流体的性质，以便提高钻探成功率。对岩性变化十分敏感的电阻率 一直是油层解释(如电阻率测井资料解释) 的重要标志性参数之一。海洋可控源电 磁法是勘探海底油气资源的一种地球物理方法，通过观测油气储层和低阻沉积围 岩之间的电阻率差异所产生的电磁异常，从而可以确定油气储层的分布范围。在 国外，海洋可控源电磁装备研制和方法研究已较为成熟，正反演研究成果较多， 而我国海洋可控源电磁探测技术的研究相对滞后。因此，我们需要加强对海洋可 控源电磁场的数值模拟和反演方法的研究工作。 本论文在总结前人研究成果的基础上，深入探讨了海洋可控源电磁一维和二 维正反演方法和关键技术，主要研究内容和取得的研究成果如下： (1) 从麦克斯韦方程组出发，利用矢量位函数，详细推导了直角坐标系下三 个方向上电偶源分别在层状各向同性和垂直各向异性介质中的电磁场解析表达 式。利用水平旋转角与倾角来分解任意取向的电偶源，求得三个方向上的电偶极 距，采用汉克尔变换计算含虚宗量贝塞尔函数的无穷积分，得到三个方向电偶源 产生的电磁场，经矢量叠加，即可求得任意取向电偶源的电磁场。采用积分法或 者剖分法计算有限长电偶源的电磁场，并分析了这两种方法的计算效率及精度。 在收发距较小的情况下，不能将有限长电偶源看作为电偶极子计算电磁场响应， 而应采用积分法来计算。用 Fortran 计算机语言编写了计算任意取向有限长电偶 源电磁响应的计算机程序，并验证了算法的正确性和精度。 (2) 以空气-海水-岩层为背景模型，分析了电磁场响应的振幅与海底岩层电阻 率之间的关系，提出了一种计算等效视电阻率的方法，以海洋油气模型为例证明 了等效视电阻率方法的有效性。 万方数据 (3) 从层状各向同性介质的电磁场解析表达式出发，详细推导了电磁场对各 层电导率的偏导数，采用高斯牛顿方法，实现了海洋可控源电磁场一维反演。 (4) 海洋可控源电磁场 2D 正演模拟中采用自适应有限元法分别求解关于总 场和异常场的边值问题。异常场法中，通过水平层状模型的一次场解析解和二次 场有限元数值解的矢量叠加得到电磁场响应。对海洋可控源电磁 2D 正演模拟中 的难点问题进行了详细的测试和分析：1）采用伪delta 函数对源进行离散，分析 了不同离散方法的精度；2 ）针对海洋可控源电磁场 2D 有限元模拟的线性方程 组，采用不同的求解器求解，确定了最优的求解器；3 ）有限元求解过程中，采 用非结构化三角单元网格模拟复杂地电模型，利用后验误差估计方法指导网格加 密细化，这样在同样网格数量下能够得