基于场延拓的海洋可控源电磁正演模拟及各向异性特征识别

海洋可控源电磁法近年来被广泛应用于海洋油气勘探之中,凭借对高阻含油储层的有效识别大大降低干井率,节约钻探成本。受沉积环境的影响,海底地层常常表现出宏观各向异性,而对海洋各向异性模拟算法进行研究可以更有效地进行海洋电磁数据解释。本文从频率域麦克斯韦方程出发,基于标量势分解理论,通过标量势函数与场的连续性条件,在波数域中将电磁场分别向海底深部和海水中延拓,并在海底耦合到发射源上完成电磁场的求解,进而利用汉克尔变换得到了空间域电磁场。通过对海底介质一维各向异性模型电磁响应特征的分析研究发现,覆盖层各向异性对电磁响应有明显的影响;然而,海底高阻层仅垂向电阻率和厚度对海洋电磁响应产生明显影响。通过极性图的方式可以对海底介质各向异性特征进行有效识别。     

各向异性介质中的浅海海洋可控源电磁响应特征

由于受到空气波的影响,浅海海洋可控源电磁数据对海底储油层的反映较弱,如何对浅海数据进行处理和解释一直是海洋电磁理论研究的热点。随着近年来海洋电磁理论的不断完善,浅海数据已经可以被较好地处理与反演,但是其解释水平仍然受基本理论研究不足的制约。针对这一现状,本文开展了海底电各向异性对浅海数据影响规律的研究。具体方法为:利用欧拉旋转建立不同的海底电性各向异性模型,然后采用交错网格有限差分法计算浅海海洋可控源电磁响应,最后通过分析同线情况下电场E_x分量的振幅和相位曲线特征以及海底电场及电流密度分布规律,分析各向异性对浅海海洋可控源电磁响应影响的物理机制,并讨论浅海各向异性情况下海洋电磁对高阻储油层的识别能力。得出的结论为各向异性介质中的浅海海洋电磁响应特征与深海有较大区别,在进行数据的处理、反演和解释时应区别于深海情况。     

时间域航空电磁数据加权横向约束反演

传统的单点一维反演技术处理时间域航空电磁数据得到的电阻率和层厚度容易呈现横向不连续性,导致层界面不光滑。即使相邻测点的反演结果也会出现突变现象,与实际地质情况不符,给解释工作带来困难。笔者针对时间域航空电磁数据空间采样密集、在沉积岩地区工作时地下电性结构横向连续性较好的特点,研究了时间域航空电磁数据参数加权横向约束反演方法。反演时同时考虑数据拟合、相邻测点电阻率和层厚度横向约束以及深度横向约束,利用参数加权来调整对各层模型参数和深度横向光滑度的约束强度。通过对理论和实测数据反演并与传统的单点一维反演结果进行对比,验证了拟二维横向加权约束反演是处理时间域航空电磁数据的有效方法。横向约束反演结果电阻率和层厚度横向光滑连续,层界面清晰。引入参数加权实现对不同参数横向光滑度和连续性的制约,可进一步提高反演解的稳定性。加权因子对反演结果有较大影响,在实际应用中,应该综合考虑地质资料选择合适的加权因子。     

起伏地表频域/时域航空电磁系统三维正演模拟研究

由于航空电磁系统具有工作频率低、时间延迟短等特点,地形对航空电磁响应有很大影响,忽略地形影响会给航空电磁数据解释造成很大误差.本文将基于非结构化网格的矢量有限元法应用于模拟起伏地表条件下频域/时域(FD/TD)三维航空电磁系统响应.该方法由于采用非结构网格,与传统的结构化网格电磁正演算法相比,能更好地拟合地形和地下不规则异常体,提高对不规则地形和地下介质航空电磁响应的计算精度.通过将计算结果与半空间模型的半解析解及已发表的结果进行对比,检验了本文算法的精度.通过对典型山峰和山谷地形航空电磁响应分析对比,总结了地形对航空电磁响应的影响特征.研究结果对航空电磁地形效应的识别和校正具有指导意义.     

时间域航空电磁扩散特征和成像深度研究

为了研究时间域航空电磁系统扩散特征,基于连续性边界条件将电磁场向下延拓得到地下介质中各点处的频率域响应,通过傅里叶变换将其变换到时间域.计算阶跃波形发射时垂直磁偶极子和水平磁偶极子在均匀半空间介质和两层介质中产生的电磁场,通过电流密度矢量图和随时间变化的电流密度等值线展示电磁场在地下介质中的扩散过程.电磁场在地下的扩散受电阻率结构的影响.在良导地区电磁场扩散慢、衰减快,而在高阻地区电磁场扩散快、衰减慢.垂直磁偶极子在地下产生的感应电流形成一个电流环随时间向下、向外扩散.水平磁偶极子在地下产生的感应电流形成两个互相叠加的电流环.时间域电磁场扩散代表着真正物理意义上的电磁扩散,因为它展示了地下电磁场空间分布随时间的变化和场强的衰减.基于对感应电流环的研究,我们探究了扩散深度和时间域航空电磁数据成像深度的关系.时间域电磁场扩散规律的研究不仅加深了对电磁扩散特征的理解,更有助于提高航空电磁数据解释水平.     

任意各向异性介质中三维可控源音频大地电磁正演模拟

在一些地层层理发育的地区,地下介质存在显著的电各向异性,此时基于各向同性模型解释含各向异性效应的可控源音频大地电磁（CSAMT）测深观测数据会导致错误的结果.本文通过引入3×3的对称正定张量表征电导率各向异性,采用非结构四面体网格和矢量有限元方法离散电场满足的矢量Helmholtz方程,并将电磁场源等效为系列电偶极子,实现任意各向异性介质中CSAMT高效数值模拟.本文首先通过层状各向异性模型检验三维有限元算法的精度和有效性,进一步建立三维地电模型研究异常体各向异性和围岩各向异性对CSAMT响应的影响,最后使用视电阻率极性图来识别各向异性电导率主轴方向.数值模拟结果表明,各向异性电导率对CSAMT视电阻率幅值及分布规律都有很大影响,视电阻率极性图能够很好地识别各向异性主轴方向.     

基于谱元法的频率域三维海洋可控源电磁正演模拟

高精度、快速有效的正演模拟算法是三维电磁正反演的前提.为了提高海洋电磁三维数值模拟的精度和效率,本文提出利用基于Gauss-Lobatto-Chebyshev（GLC）基函数的谱元法进行海洋可控源三维电磁正演模拟.谱元法结合有限元法和谱方法的优点.我们通过应用伽辽金加权残差法离散二次电场矢量亥姆赫兹方程,在单元内选择混合阶GLC多项式的张量积作为高阶矢量插值基函数,在求解大型稀疏线性方程组时利用直接求解器进行快速求解,从而实现了三维海洋可控源电磁快速高精度正演模拟.一维和三维模型正演结果验证了本文算法的有效性和准确性.典型模型的数值结果表明谱元法是一种有效的三维海洋可控源电磁正演数值方法,能在稀疏网格剖分情况下获得精确的海洋电磁正演模拟响应.     

基于小波变换的重力压缩正演和多尺度反演研究

解决多解性和计算效率等问题是重力反演中的主要难点。针对计算效率问题,本文基于小波压缩理论,提出一种改进的重力矩阵压缩正演方法。该方法通过引入一种新的灵敏度矩阵排序规则,降低灵敏度峰值和波动、减少高频信息,在保证精度的同时改善小波变换后灵敏度矩阵的稀疏性、提高压缩比,从而降低内存需求,提高正演速度。正演模拟分析结果表明,本文提出的正演矩阵压缩方法在保证正演精度的前提下压缩比可达300以上。其次,本文将基于小波变换的多尺度反演应用到三维重力反演中,通过将重力反演问题分解成不同尺度的子问题,充分利用不同尺度的数据信息来改善重力反演的多解性和稳定性问题。最后,我们分别对理论数据和实测数据进行常规聚焦反演和多尺度反演以证明多尺度反演在重力反演中的有效性。     

起伏海底地形时间域海洋电磁三维自适应正演模拟

本文基于自适应非结构有限元算法实现海洋电磁起伏海底地形三维正演模拟.通过采用隐式后推欧拉时间离散技术,保证在较大的时间步长条件下获得正确结果.为获得多时间道海洋电磁正演模拟的有效网格,我们采用基于法向电流连续的后验误差估计的自适应方法和网格融合技术;同时为了控制网格数量和保证正演模拟稳定性,我们还在网格融合过程中应用了随机网格挑选技术.对于方程组求解我们使用MUMPUS直接求解器.当时间步长不变时,只需对系数矩阵进行一次分解,大大提高计算速度.将本文计算结果与半空间模型解析解进行对比,验证了本文算法精度.针对海洋电磁半拖曳式和双船拖曳式工作方式,我们通过典型模型的模拟计算,研究海底地形影响及海底高阻层识别特征.     

ElasticSearch分布式搜索引擎在地名地址检索中的应用

针对传统全文搜索引擎Lucene无法满足海量地名地址数据高速分布式检索的问题,本文通过对分布式搜索引擎ElasticSearch的集群架构、索引分片与复制集、地理要素索引存储、数据检索及位置聚合等关键点进行详细研究,结合地名地址库数据结构的特点,提出了一种面向Web电子地图开发的分布式地名地址搜索服务构架,实现了地名地址模糊搜索、专题搜索、视野内搜索、周边搜索等功能。该架构具备可伸缩性、容错性、支持大规模并发访问的特点。