简介：标量CSAMT只适合一维及测量方向与构造方向垂直的二维情况,对于复杂的三维地电结构,CSAMT需采用张量测量。本文试图采用矢量有限元法实现三维张量CSAMT的正演模拟。为了验证算法的正确性,本文在层状介质中计算了三维CSAMT远区的电场,磁场及阻抗张量,并且与层状介质中的理论解进行了比较,接着还模拟了均匀半空间中含有三维异常体的模型,并且分析了四个阻抗张量、视电阻率及阻抗相位的响应特征。得出如下结论：采用矢量有限元法来模拟三维张量CSAMT,其电磁场及阻抗张量的实虚部计算精度都比较高,并且该方法本身满足电场法向不连续,不用进行散度校正。

简介：角度域共成像点道集(ADCIGS)是偏移速度分析和振幅随角度变化分析(AVA)的基础数据。传统Kirchhoff叠前深度偏移(KPSDM)按偏移距组织数据,能方便的输出偏移距域共成像点道集(ODCIGS),其高效的角度道集输出是有挑战的。本文提出基于旅行时梯度场的KPSDM角道集输出方法。其核心步骤为:(1)利用任意介质中的动态规划法旅行时计算方法提供炮点和检波点的旅行时场;(2)根据旅行时场的梯度方向计算反射张角;(3)在偏移过程中抽取ADCIGS。由于本文旅行时计算方法没有射线阴影区,也没有对速度光滑性的要求,其角度道集输出在阴影区比传统射线追踪更有优势。基于该角度道集输出方法,本文发展了一种适合大规模三维地震数据的KPSDM及角道集输出的并行实现方案。其基本思想是:(1)按照炮数据来组织输入数据;(2)旅行时场的输入与单炮覆盖范围相联系以节省内存;(3)多炮数据间采用MPI并行处理,单炮深度切片之间采用OpenMp并行处理,可进一步提高内存利用率和并行力度。数值试验结果证明本文角度道集生成方法的优越性和本文实现方案的有效性。

简介：被动源地震数据包含丰富的低频信息,本文有效地提取并利用这些信息对缺失低频的主动源地震数据进行低频重构,提出了基于多正弦窗的被动源多窗谱重构方法,并给出了相应的多震源多道重构公式。与常规互相关法和常规反褶积法重构的被动源记录相比,该方法能重构出更为准确的相对振幅信息。通过分析被动源数据重构前后的频谱特性,发现被动源的低频特性在重构和去噪处理后能更明显的体现出来。并提出了一种用被动源数据重构主动源低频信息的方法,即在功率谱上进行匹配,并在频域进行补偿和平滑。最后进行了数值算例的验证,对低频重构后的数据进行了叠前深度偏移处理。能量匹配方法能够用被动源的低频信息有效地重构主动源缺失的低频信息,低频重构后的记录在偏移成像中能体现更多的细节信息和深部构造。

简介：标量CSAMT只适合一维及测量方向与构造方向垂直的二维情况,对于复杂的三维地电结构,CSAMT需采用张量测量。本文试图采用矢量有限元法实现三维张量CSAMT的正演模拟。为了验证算法的正确性,本文在层状介质中计算了三维CSAMT远区的电场,磁场及阻抗张量,并且与层状介质中的理论解进行了比较,接着还模拟了均匀半空间中含有三维异常体的模型,并且分析了四个阻抗张量、视电阻率及阻抗相位的响应特征。得出如下结论:采用矢量有限元法来模拟三维张量CSAMT,其电磁场及阻抗张量的实虚部计算精度都比较高,并且该方法本身满足电场法向不连续,不用进行散度校正。

简介：坐标变换法通过将物理空间的曲网格映射为计算空间的矩形网格,将起伏地表转化为水平地表,同时将物理空间的波动方程转化为计算空间的波动方程,在计算空间完成数值模拟,坐标变换的方法对处理起伏自由边界具有较好的适应性和应用效果。本文在传统坐标变换方法的基础上,根据计算区域速度差异采用不同的网格大小和采样时间步长,提出了一种基于时空双变网格的起伏地表坐标变换正演模拟方法。在编程实现算法的基础上,通过典型模型波场模拟试算结果分析可知：（1）变网格方法与常规方法波场模拟误差在0.5%左右;（2）变网格方法计算效率视不同的变网格区域面积及变网格大小可提高几倍量级,在本文模型和计算参数下提高约5倍。（3）在满足模拟精度及频散条件要求下,变网格方法较全局细网格算法能显著节约计算内存。为此,针对起伏地表数值模拟,本文方法具有较高的模拟计算精度和一定的适应性。

简介：无网格法形函数构造不依赖预定义的单元，具有计算精度高、处理复杂模型便利等优点。本文介绍了无单元Galerkin法（EFGM）、点插值法（PIM）与径向基点插值法（RPIM）三种全域弱式无网格法的近似原理及特点；以二维泊松方程为例研究了支持域无量纲尺寸、场节点与背景网格设置对无网格法计算精度的影响。将RPIM与EFGM应用于频率域线源二维正演，给出了RPIM形状参数的推荐值；分析了均匀介质模型大地电磁（MT）二维正演无网格法边界条件直接加载与罚函数法加载的精度差异，结合PIM与RPIM边界条件加载便利及EFGM计算复杂模型精度高的优势，提出了EFG—PIM及EFG．RPIM耦合算法，数值计算结果验证了耦合算法的有效性。研究发现：无网格法及其耦合方法适用于电磁法数值模拟；支持域无量纲尺寸取1．0时无网格法精度与效率高，场节点与背景网格重合时计算效果佳；泊松方程求解PIM及RPIM精度较EFGM低，计算均匀介质MT响应精度较EFGM高；RPIM改善了PIM计算涉及的奇异性问题，对应支持域无量纲尺寸选择空间大。

简介：天然气水合物沉积层在常规叠后地震剖面上，一般表现为似海底反射（BSR）特征，常作为水合物识别的重要标志，但由于地震数据分辨率的限制与多解性的存在，类似的反射特征不一定是水合物的表现。本文先从南中国海神狐海域地质条件分析人手，综述研究区域天然气水合物发育的构造、沉积与运移环境；其次以加权方程与Zoeppritz公式的正演模拟为基础，讨论天然气水合物与BSR特征的关系；然后选取经过三口已钻获天然气水合物站位的地震测线，提取角度域共成像道集用于BSR的AVA特征研究；最后依据叠前地震弹性参数同步反演思路，对水合物重点赋存区域开展精细预测研究，获取了研究区域含水合物沉积层稳定带的分布特征与饱和度数。该方法可作为采用多道地震数据估算天然气水合物含量的一种有效技术。