[发明专利]一种直流电阻率小波伽辽金三维正演方法

说明书

本发明属于地球物理勘探领域，具体地而言为一种直流电阻率小波伽辽金三维正演方法，该方法包括：对地电模型计算区域进行不均匀网格剖分为多个网格单元，得到网格单元编号，节点编号和坐标；读取设计的地电模型参数，包括背景层参数、网格单元编号、节点编号和坐标；计算设计的背景层的地下背景介质的一次场电位值；对所有网格单元进行循环计算单元系数矩阵，接着总体合成所有网格单元的系数矩阵；根据一次场电位值计算直流电阻率法小波伽辽金线性方程组的右端项；加载本质边界条件，求解直流电阻率法小波伽辽金线性方程组，得到各个节点的二次场电位，进而得到总场电位值。解决现有的方法存在的计算速度慢，计算精度低的问题。

技术领域

本发明属于地球物理勘探领域，具体地而言为一种直流电阻率小波伽辽金三维正演方法。

背景技术

随着经济建设的快速持续发展，矿产资源的短缺形势日益凸显，当前传统矿产勘探开发技术的理论研究不足和探测深度以及精度上的局限性以及目标体越隐蔽、复杂，导致油气勘探开发难度相应越来越大。因此，开展矿产资源勘探技术研究，提高勘探的精度，对于提高矿产资源勘探开发效益，具有重要的意义。

众所周知，直流电阻率法作为地球物理勘探中最古老的勘探方法之一，早已广泛应用于矿产资源勘探，水文地质调查，工程勘察等领域。随着勘探复杂度的加大，如何提高勘探的精度已然成为国内外研究的热点问题。正演是反演的基础，因此，研究高精度的电阻率三维正演对实现高效的三维反演是具有重要的理论意义和实际价值。目前为止，常用的直流电阻率数值模拟方法主要有积分方程法(IEM),有限差分法(FDM)，有限单元法(FEM)等。

直流电阻率法三维正演中存在计算速度慢，计算精度低等问题。

发明内容

本发明所要解决的技术问题在于提供一种直流电阻率小波伽辽金三维正演方法，解决现有的方法存在的计算速度慢，计算精度低的问题。

本发明是这样实现的，

一种直流电阻率小波伽辽金三维正演方法，该方法包括：

对地电模型计算区域进行不均匀网格剖分为多个网格单元，得到网格单元编号，节点编号和坐标；

读取设计的地电模型参数，包括背景层参数、网格单元编号、节点编号和坐标；

计算设计的背景层的地下背景介质的一次场电位值；

对所有网格单元进行循环计算单元系数矩阵，接着总体合成所有网格单元的系数矩阵；

根据一次场电位值计算直流电阻率法小波伽辽金线性方程组的右端项；

加载本质边界条件，求解直流电阻率法小波伽辽金线性方程组，得到各个节点的二次场电位，进而得到总场电位值。

进一步地：

所述直流电阻率法小波伽辽金线性方程组为：

K_su_s＝b

其中：

K_s为系数矩阵，u_s为求解域中各个节点待求的场值，b为右端项，均为2-term连接系数。

进一步地：

所述直流电阻率法小波伽辽金线性方程组通过下列方法得到：

采用狄利克雷边界条件的二次场电位所满足的边值问题表达式：