专利名称：一种电磁数据地形改正方法
技术领域：
本发明涉及地球物理勘探方法，是一种电磁勘探数据的数据处理技术中电磁数据地形改正方法。
背景技术：
现在的电磁方法勘探，包括大地电磁法(MT)、音频大地电磁法(AMT)、连续大地电磁剖面法(CEMP)勘探主要集中在复杂山前带等地区，在复杂地形条件下采集的电磁数据会产生一定的畸变，这些资料如果不进行地形改正，就会影响后续的处理结果。地形对电磁资料的影响问题自20世纪80年代开始就得到国内外学者的关注，Redding、Winnamaker等人用有限元法在二维MT地形模拟方面做了很多工作。到20世纪90年代末期，Chouteau、王绪本、晋光文、张翔等人也分别研究了二维MT的地形影响，并公 开了相应的改正方法，如比值法，位场延拓法等。自从Chouteau公开二维MT地形改正方法以来,很多文献公开了比值法均勻半空间带地形的模拟数据进行地形改正，如下式(I)和(2)构成了传统的比值法地形改正公式P c(f, x) = P obs(f, x) P o/P T(f x)(I)
φΛ ^χ) = Ψο^( ^χ)-Ψτ( ^χ)+π/4(2)其中P e(f，X)分别为地形改正后的视电阻率与相位，P Obs(f X)、U/，x)分别为观测的视电阻率与相位，PT(f，Χ)、Α(/，χ)分别为带地形均匀半空间模型正演的视电阻率与相位，P ^选择为测线地表的统计电阻率值。传统的比值法对地形影响的压制有一定的作用，但采用比值法进行地形影响改正后的视电阻率分布相对于水平地形模型的结果仍有差异，其残余部分与地形呈镜像关系。并且这些方法大都集中在室内简单地形的理论模拟，很少见到对大数据量的实测数据进行地形改正的报道。由于模拟地形时需要消耗大量计算机的内存资源，从而限制了处理的数据点数，不能满足现在长电磁剖面的实测数据要求。传统的比值法进行电磁数据地形改正后效果不好。
发明内容
本发明目的是提供一种与地形存不在镜像关系，适合并行机群使用的电磁数据地形改正方法。本发明具体方法步骤如下I)采集实测电磁数据，对数据的视电阻率和相位做平均；步骤I)所述的电磁数据包括大地电磁数据(MT)、音频大地电磁数据(AMT)、连续剖面法大地电磁数据(CEMP)。
步骤I)所述的平均是对所有实测电磁数据的视电阻率和相位做几何和算术平均。2)对平均视电阻率和相位曲线进行反演，将平均实测曲线的频率-视电阻率关系变为深度-电阻率的关系，获取改正区域的背景层状电阻率模型，正演获得视电阻率Play(f，x)和相位~(·^);步骤2)所述的反演采用一维反演方法，正演为一维正演的方法。步骤2)所述的反演包括广义逆、奇异值分解、马奎特、奥克姆、模拟退火、遗传算法、量子退火、演化算法。
3)根据测点相对高程数据进行横纵向网格剖分，由获取的改正区域的背景层状电阻率模型，在每个测点的深度方向上进行电阻率线性插值，从而获得随地形起伏的深度-电阻率模型，对该模型按照频率进行并行二维正演，获得带地形的层状模型二维正演视电阻率和相位<。(/，1)和<。(/，4 ;步骤3)所述的测点相对高程数据指的是，以第一个测点的高程数据为基准零点，其他测点高程逐一与第一个测点的高程数据想减，在这之后得到的所有测点的高程数据。步骤3)所述的并行二维正演，指的是利用MPI并行函数库，按照频率进行并行设计，执行并行二维正演。步骤3)所述的二维正演包括有限元法、有限差分法。4)利用以下公式(3)和⑷求取一层状模型的视电阻率P e(f，x)和一层状模型的相位Μ/，χ).
权利要求
1.一种电磁数据地形改正方法，特点是采用以下步骤实现 1)采集实测电磁数据，对数据的视电阻率和相位做平均； 2)对平均视电阻率和相位曲线进行反演，将平均实测曲线的频率-视电阻率关系变为深度-电阻率的关系，获取改正区域的背景层状电阻率模型，正演获得视电阻率Play(f，x)和相位 (/，x)； 3)根据测点相对高程数据进行横纵向网格剖分，由获取的改正区域的背景层状电阻率模型，在每个测点的深度方向上进行电阻率线性插值，从而获得随地形起伏的深度-电阻率模型，对该模型按照频率进行并行二维正演，获得带地形的层状模型二维正演视电阻率和相位<。(/，4和<。(/，4 ； 4)利用以下公式(3)和(4)求取一层状模型的视电阻率Pjf，x)和一层状模型的相位呎(/，x);
2.根据权利要求I所述的方法，特点是步骤I)所述的电磁数据包括大地电磁数据、音频大地电磁数据、连续剖面法大地电磁数据。
3.根据权利要求I所述的方法，特点是步骤I)所述的平均是对所有实测电磁数据的视电阻率和相位做几何和算术平均。
4.根据权利要求I所述的方法，特点是步骤2)所述的反演采用一维反演方法，正演为一维正演的方法。
5.根据权利要求I所述的方法，特点是步骤2)所述的反演包括广义逆、奇异值分解、马奎特、奥克姆、模拟退火、遗传算法、量子退火、演化算法。
6.根据权利要求I所述的方法，特点是步骤3)所述的测点相对高程数据指的是，以第一个测点的高程数据为基准零点，其他测点高程逐一与第一个测点的高程数据想减，在这之后得到的所有测点的高程数据。
7.根据权利要求I所述的方法，特点是步骤3)所述的并行二维正演，指的是利用MPI并行函数库，按照频率进行并行设计，执行并行二维正演。
8.根据权利要求I所述的方法，特点是步骤3)所述的二维正演包括有限元法、有限差分法。
9.根据权利要求I所述的方法，特点是步骤4)中的延拓高度hi为步骤3)中所述的各个测点对应的相对高程。
10.根据权利要求I所述的方法，特点是步骤4)中的公式(3)、(4)和(5)的计算包含在频率并行正演之中，在每个频率正演计算完毕之后就进行这三个公式的计算。
11.根据权利要求I所述的方法，特点是步骤5)所述的上下延拓改正在上延时，式(5)中第一项取负号，下延时则取正号。
12.根据权利要求I所述的方法，特点是步骤5)所述的测点相对高程为负时，进行上延，当测点相对高程为正时，进行下延。
全文摘要
一种电磁数据地形改正方法，实测电磁数据，对视电阻率和相位做平均，将频率-视电阻率关系变为深度-电阻率的关系，获得随地形起伏的模型，按照频率进行并行二维正演，获得带地形的层状模型二维正演视电阻率和相位，求取一层状模型的视电阻率，根据大地电磁场一维阻抗递推公式和测点相对高程，对进行上下延拓改正。本发明采用与地形相关的层状模型以及上下延拓技术，可以压制应用传统方法后视电阻率和相位剖面与地形存在镜像关系的问题，改善了地形改正的效果；同时通过采用并行处理方法进行电磁数据的改正，可以提高地形改正的能力与处理效率。
文档编号G01V3/38GK102778699SQ20111012460
公开日2012年11月14日 申请日期2011年5月13日 优先权日2011年5月13日
发明者何展翔, 胡祖志 申请人:中国石油天然气集团公司, 中国石油集团东方地球物理勘探有限责任公司