漂航状态矢量修正海流测量数据方法研究

针对舰船在漂航状态(停航不抛锚)时,海流计实测数据不准问题,基于天气海况条件良好的试验海区,分析海流计数据误差原因,提出流速矢量修正海流计实测数据的方法,将修正后结果与《T-D值表》潮流推算值对比发现:13:20时刻的流向误差为0.6%,流速误差为3%;14:20时刻流向误差为1.6%,流速误差为0%,修正效果理想,误差较低,所以此方法对于漂航状态海流测量数据的修正是有效的;由此可知在浅海、近岸海域,漂航状态的舰船在海面风较小时,海流计的测量数据经流速矢量修正后精度较高,可用性较强,可以解决漂航状态海流数据不准问题。     

MCSEM电磁场能流密度分布特征研究

采用海洋可控源电磁方法进行海底油气和水合物的勘查,需事先在海底布放接收机阵列,而后用一艘携带深拖缆的科考船将大功率电磁发射机牵引至海底,沿预设路线慢速走航,并将大功率电磁场导入海底.海底的人工电磁场中,空气波和直达波对提取勘探目标体的异常信息不做贡献.只有当折射波经过被测目标体时,产生电磁感应,海底的接收机对此感应信号进行数据采集,经数据处理,揭示被测目标体的几何大小和物理属性.异常体的有用信息与导入海底电磁场的能流密度(坡印廷矢量)有着密切的关系,以致于后者决定了前者的有效探测范围、有效带宽以及尽可能达到的最佳测量精度.能流密度在海底的分布情况十分复杂,但在均匀全空间中可较清楚地认识其分布特征,可为复杂地质条件的情况作理论铺垫.以此推广,在水平层状介质的地电条件下,建立能流密度的数理模型,利用海洋生产作业中的典型物性与几何参数,对海底以下有覆盖层的高阻体(可模拟油气或天然气水合物)进行正演,建立电磁场能流密度分布特征与数据采集技术的有机联系,为海洋可控源电磁发射机和接收机的优化改进提供理论依据.     

可控源电磁场三维自适应矢量有限元正演模拟

本文实现了可控源电磁(CSEM)场三维自适应矢量有限元正演算法,该算法采用非结构四面体单元进行三维网格剖分,能够真实模拟地形起伏和复杂电性异常体.采用一次场和二次场分离的方式计算电磁场响应,能够有效解决有限元模拟中的源点奇异性,提高场源附近电磁场数值精度,其中一次场利用CSEM一维正演算法解析求得,二次场采用矢量有限元方法求得.并利用基于后验误差估计的自适应网格细化算法指导网格细化,以减少人为设计网格导致的误差.通过一维和三维模型的数值模拟,验证了本文算法的有效性:一维模型有限元解与解析解吻合得很好,电磁场振幅相对误差在1%左右,相位差整体小于1°;三维模型有限元解与有限体积解吻合得也很好.模拟了一个含三维倾斜板状异常体的可控源电磁响应,表明了本文算法模拟复杂地电结构电磁场的能力和有效性.     

电导率任意各向异性海洋可控源电磁三维矢量有限元数值模拟

现有海洋可控源电磁三维数值模拟方法大多基于电导率各向同性介质理论,不能模拟海底地层电导率各向异性的实际情况.本文给出了电导率各向异性三维介质中电性源海洋可控源电磁二次电场的边值问题以及相应的变分问题,采用长方体单元对研究区域剖分,将场分量定义在剖分单元的边上,利用矢量有限单元法求解变分问题,实现了电导率任意各向异性海洋可控源电磁三维矢量有限元数值模拟.这个新的正演方法可以计算电导率任意各向异性三维地电模型的海洋可控源电磁响应,基于二次场矢量有限元法直接求解电磁场,避免了传统有限元方法可能遇到的伪解问题和难于处理电场法向分量不连续的问题,提高了数值模拟计算精度.一维电导率各向异性模型电磁场数值解与解析解吻合得相当好,无论在源附近还是远离源处相对误差均不超过1%.电导率各向异性二维模型的计算结果与已有文献采用的非结构有限元模拟结果十分吻合.三维地电模型数值模拟结果显示,电导率各向异性张量电导率主轴分量和欧拉角对不同装置海洋可控源电磁响应均有着明显的影响.     

频率域电磁法三维有限元正演线性方程组迭代算法

在三维频率域电磁法的正演模拟方法中,有限元方法具有计算精度高、适应性强的优点,近年来来得到了越来越多的关注.在正演过程中,主要的计算量集中在求解由偏微分方程组离散得到的线性方程组上,因此求解线性方程组关系着正演计算速度以及模拟精度.由于由有限元方法离散得到的复系数线性方程组条件数非常大,使用常规的迭代法和预条件很难收敛.目前大多数的研究工作采用直接解法,需要大量的计算机内存,限制了可求解问题的规模.本文研究了线性方程组的迭代解法,通过将复系数线性方程组转化为其实对称形式,构造分块对角预条件.在应用预条件的过程中,需要求解两个较小的实数方程,通过辅助空间解法求解.本文的算法适用于可控源电磁法和大地电磁法,对一系列的数值算例的模拟结果证明了迭代算法的效率,结果表明迭代算法可以在小于20次迭代内收敛,同时迭代次数与模型电阻率、问题规模和频率无关.     

利用免费DEM数据及制图软件绘制跨断层场地图

涉及跨断层测量的很多专业报告都需要以图片形式展示观测场地测网布设、地形地貌等基础信息。由于测网改造及城市建设,很多观测场地的场地图已不能反映场地的真实面貌,有必要绘制新的场地图。在不进行实地测量和没有航测设备的情况下,利用免费的DEM数据及多个制图软件可绘制出满足展示需求的场地图。     

基于波场分离的弹性波逆时偏移

尽管叠前逆时偏移成像精度高,但仅针对单一纵波的成像也可能形成地下介质成像盲区,由于基于弹性波方程的逆时偏移成像可形成多波模式的成像数据,因此弹性波逆时偏移成像可提供更为丰富的地下构造信息.本文依据各向同性介质的一阶速度-应力方程组构建震源和检波点矢量波场,再利用Helmholtz分解提取纯纵波和纯横波波场,使用震源归一化的互相关成像条件获得纯波成像,避免了直接使用坐标分量成像而引起的纵横波串扰问题.针对转换波成像的极性反转问题,文中提出一种共炮域极性校正方法.为有效节约存储成本,也提出一种适用于弹性波逆时偏移的震源波场逆时重建方法,在震源波场正传过程中,仅保存PML边界内若干层的速度分量波场,进而逆时重建出所有分量的震源波场.本文分别对地堑模型和Marmousi2模型进行了弹性波逆时偏移成像测试,结果表明:所提出的共炮域极性校正方法正确有效,基于波场分离的弹性波逆时偏移成像的纯波数据能够对复杂地下构造准确成像.     

面向目标自适应海洋可控源电磁三维矢量有限元正演

本文实现了一种面向目标自适应海洋可控源电磁三维矢量有限元方法.为满足三维复杂电性结构模拟的需求,网格剖分采用非结构化六面体.在组装刚度矩阵之后,形成的大型复数线性方程组分解为等价的实数形式,利用带预条件的广义最小残差法进行求解.在获得微分方程的解之后,为提高解的准确性,通过面向目标的自适应误差估计来指示网格细化,重点加密能使观测点数值模拟精度提高的网格.对于大规模三维数据,为了使模型空间的并行计算达到均衡负载的效果,我们使用METIS函数库来进行网格计算任务量的划分.最后,通过对比一维解析解与三维自适应矢量有限元计算结果,验证了程序的正确性;通过自适应过程中误差指示子的分布,验证了面向目标自适应的有效性;通过对三维复杂模型进行均衡负载下的并行计算,测试了程序的可扩展性.     

三台山森林公园建设对宿迁台地磁干扰分析

在三台山森林公园建设过程中,宿迁地震台地磁观测受到严重干扰。分析了对地磁场造成干扰的主要因素,通过对比干扰前后的地磁场观测数据,运用矢量合成的规则,计算出干扰磁场的7个要素,再使用IAGA推荐的理论公式对建筑中铁磁性材料造成的干扰进行计算,得到铁磁性材料产生的干扰数据情况。