基于矢量有限元的带地形大地电磁三维反演研究

研究了基于矢量有限元方法的大地电磁带地形三维反演算法并开发了三维反演计算程序代码.在大地电磁场正演数值模拟方面,采用并行直接稀疏求解器PARDISO且无需进行散度校正的快速正演方案,对典型地形模型,在中等规模计算条件下,与双共轭梯度法(BICG)计算结果比较,发现PARDISO比BICG快10倍以上;通过理论模型试算,并与前人的有限元法计算结果对比,验证了带地形三维正演计算程序的正确性.在反演方面,本研究基于共轭梯度方法编写了大地电磁带地形三维反演代码,为了避免直接求取雅可比矩阵,将反演中的雅可比矩阵计算问题转为求解两次“拟正演”问题,进而将PARDISO的快速正演方案应用于“拟正演”问题的求解,以提高反演计算效率.利用开发的反演算法对多个带地形地电模型的合成数据进行了三维反演,反演结果能很好地重现理论模型的电性结构,验证了本文开发的三维反演算法的正确性和可靠性.最后,利用该算法反演了某矿区大地电磁实测数据,反演得到的三维电性结构清晰地反映了研究区的地电特征,将反演结果与该区已有地质资料结合进行解释,应用效果明显,进一步验证了本文算法的有效性.     

空间波数混合域磁异常场积分解三维数值模拟

本文提出一种空间波数混合域磁异常场三维数值模拟方法.该方法利用磁位三维空间域积分为卷积的特点,沿水平方向进行二维傅里叶变换,把空间域磁位满足的三维积分问题转化为不同波数之间相互独立的垂向一维积分问题.保留垂向为空间域,优势之一在于便于浅层单元剖分可适当加密,随着深度增加,单元剖分适当稀疏,可以准确模拟任意复杂地形和磁性体的磁异常,兼顾了计算精度与计算效率;优势之二在于一维积分垂向可离散为多个单元积分之和,每个单元采用二次形函数表征磁化强度,可得出单元积分的解析表达式,计算精度高、效率高.该方法充分利用一维形函数积分的高效和高精度、快速傅里叶变换的高效性及算法高度并行性,实现了磁异常场高效、高精度的数值模拟.设计棱柱体模型,将模型解析解与空间波数混合域法的数值解对比,结果表明该方法计算精度高、效率高.设计了组合棱柱体复杂模型,对比分析了标准FFT扩边法与Gauss-FFT法的计算精度与计算效率,总结了标准FFT的扩边系数选取策略.针对任意复杂地形条件下的磁异常模拟问题,本文提出一种适用于起伏地形条件下的磁异常场快速计算方法,并对其有效性进行了验证.     

三维地震属性技术在湖区的应用效果

微山湖湖区表层地震地质条件十分复杂,芦苇丛生,三维地震勘探野外数据采集采用自行研制的钻井平台加套管护壁的成孔工艺,沼泽检波器加上加长尾椎,提高了检波器的耦合效果,获取了高信噪比的三维地震数据。资料解释应用了三维可视化、相干数据体、地震属性技术进行三维构造精细化解释,明显提高了地震解释识别地下小地质构造及地震异常体的能力,它完善和发展了三维高分辨率地震勘探技术,扩大了地震勘探的应用范围,为解决煤矿开采阶段所遇到的构造、煤层等地质问题提供了一种综合应用的解释方法,在湖区的勘探开发中,取得了明显的地震地质效果。     

维波场延拓复杂表层模型校正方法及应用

如何正确地消除复杂地表对地震波场的影响，提高地下构造成像的质量一直是中国西部复杂地区地震勘探中存在的难题.本文在三维复杂表层速度模型层析反演\[1\]的基础上，研究了关于复杂地表的静校正问题，提出用三维波动方程在炮集上对地震波场进行正、反向延拓，消除复杂地表对波场的影响，实现三维复杂表层模型校正.理论和实际应用证明，该方法已超越常规静校正的含义，属时变校正方法.用本方法处理复杂地表问题，不但能消除表层对不同深度反射波产生的不同时差影响，提高叠加剖面质量，而且能使校正后的地震波场保持波动特征不发生畸变，可为建立正确的深层速度模型和波动方程叠前深度偏移奠定良好的基础.     

Gravity and magnetic data inversion for 3D topography of the Moho discontinuity in the northern Red Sea area, Egypt

The main goal of our study is to investigate 3D topography of the Moho boundary for the area of the northern Red Sea including Gulf of Suez and Gulf of Aqaba. For potential field data inversion we apply a new method of local corrections. The method is efficient and does not require trial-and-error forward modeling. To separate sources of gravity and magnetic field in depth, a method is suggested, based on upward and downward continuation. Both new methods are applied to isolate the contribution of the Moho interface to the total field and to find its 3D topography. At the first stage, we separate near-surface and deeper sources. According to the obtained field of shallow sources a model of the horizontal layer above the depth of 7 km is suggested, which includes a density interface between light sediments and crystalline basement. Its depressions and uplifts correspond to known geological structures. At the next stage, we isolate the effect of very deep sources (below 100 km) and sources outside the area of investigation. After subtracting this field from the total effect of deeper sources, we obtain the contribution of the Moho interface. We make inversion separately for the area of rifts (Red Sea, Gulf of Suez and Gulf of Aqaba) and for the rest of the area. In the rift area we look for the upper boundary of low-density, heated anomalous upper mantle. In the rest of the area the field is satisfied by means of topography for the interface between lower crust and normal upper mantle. Both algorithms are applied also to the magnetic field. The magnetic model of the Moho boundary is in agreement with the gravitational one.     

基于速度-状态摩擦本构定律的三维DDA方法

原有三维非连续变形分析(DDA)方法采用常摩擦系数的Mohr-Coulomb定律作为切向破坏准则,然而当描述更大尺度构造块体的运动与变形时,常摩擦系数不再适用.速度-状态摩擦本构定律能够定量描述地震周期各阶段断层面剪应力变化,解释发震断层行为.本文将速度-状态摩擦定律与三维DDA方法相结合,首先推导了计算摩擦系数的实用公式,随后通过滑动-保持-滑动实验与速度步进实验算例对改进的三维DDA方法进行了验证.结果表明,应用速度-状态摩擦本构定律的三维DDA方法能够比较准确地模拟静摩擦的时间依赖性与动摩擦的速度依赖性,解决了将三维DDA方法在地学中应用的基本问题.     

三维动态地形中瀑布的设计与实现

目的文中了说明三维地形的模拟是进行室外可视化系统开发的一个必不可少的工作.方法文中提出了一个利用粒子系统的方法来模拟三维地形中瀑布的运动和汹涌直下.结果通过粒子作用范围及浓度的变化模拟瀑布的运动,介绍了粒子的表示和粒子系统的一般原理.结论给出了一种基于粒子系统的快速绘制三维地形中瀑布的方法.     

海洋可控源电磁三维非结构矢量有限元数值模拟

本文实现了海洋可控源电磁三维矢量有限元数值模拟.由于采用非结构四面体单元进行三维网格剖分,该方法可以模拟复杂电性异常体和海底地形.一维模型的数值模拟结果表明,电场实、虚部均与解析解吻合得相当好,计算误差基本小于1%.二维模型的计算结果与已有的二维自适应非结构有限元模拟结果吻合很好.带地形模型的数值模拟结果显示,海底地形对电场影响很大,有可能掩盖海底油气藏产生的异常.     

三维三分量CSAMT法有限元正演模拟研究初探

首先从麦克斯韦方程出发,用伽里金方法推导了三维三分量CSAMT法的有限元方程.在研究过程中,认识到加入散度条件的必要性,在公式中强加了散度条件,提高了解的完备性.其次将成功应用于二维线源频率域电磁法有限元模拟中的两种技术推广到三维中,一是边界条件统一采用一阶吸收边界,使线源产生的电磁波在边界上按波的传播规律被吸收,以降低平面波假设造成的影响;二是总体系数矩阵的存储,用两个二维数组分别记录总体系数矩阵的非零元素及其在总体结点编号中所处的位置,使总体系数矩阵的存储量达到最小的同时,物理意义明确,迭代求解时迅速简便.最后用均匀半空间模型进行了验证.     

三维泊松方程数值模拟的多重网格方法

本文简要介绍了多重网格方法的基本思想和原理,然后应用多重网格(MG)方法求解三维泊松方程,网格尺度从17×17×17逐次增加至257×257×257,并与不完全Chelesky共轭梯度法(ICCG)、Gauss直接解法进行比较,结果表明,MG方法计算速度明显优于ICCG、Gauss方法,对于129×129×129网格的三维数值模拟费时43s,比ICCG法快7倍,而对于257×257×257超大型网格的三维数值模拟也仅需412s.     

三维分块倾斜界面的反演及其应用

介绍了利用反射波走时反演介质层速度和包含直立断层的三维界面的计算方法．对各层界面利用多个倾斜或平界面方程来描述．给出了数值计算实例．结果表明，反演解与真模型很接近，说明计算方法是有效的．利用该方法处理了唐山地震区的ＰｍＰ反射波资料，获得了该区域莫霍界面的三维分块形态．唐山附近莫霍界面错断与该区域１９７６年唐山地震及一系列余震的发生有密切的联系．     

基于快速推进迎风双线性插值法的三维地震波走时计算

三维地震波走时计算技术是三维地震反演、层析成像、偏移成像等诸多地震数据处理技术中非常重要的正演计算工具.为了获得精度高且兼顾效率的三维走时计算方法:首先,在常规双线性插值公式推导过程中,充分利用平面波双线性假设的结论,获得了二元极小值超越方程的解析解,进而推导出了准确的局部走时计算公式,同时构造性地证明了该计算公式满足地震波的传播规律和Eikonal方程;其次,引入迎风差分的基本思想,提出迎风双线性插值的局部走时计算策略,该计算策略能简化算法、提高效率且保证无条件稳定性;然后,将上述计算公式和迎风双线性插值策略与常规快速推进法中的窄带技术结合,获得了一种新的基于快速推进迎风双线性插值法的三维地震波走时计算方法;最后,通过精度和效率分析检验了新算法的精度、效率和正确性,并通过计算实例验证了算法在面对复杂介质时的稳定性和有效性.     

考虑视觉传达效果的震后区域三维图像虚拟重建

传统二维震害图像方法对震后区域进行研究时,由于其拍摄角度具有局限性,震后区域图像的视觉效果不理想。提出基于三维激光扫描技术的震后区域三维虚拟重建方法,采用三维激光扫描仪测量震后区域,获取该区域的点云数据,采用Cyclone软件合并点云数据后,得到震后区域拼接后的整体点云图,将该图点云数据进行去体外孤点、去噪声点以及点云取样等处理后实施封装,在封装的点云数据上采用Sketch模型实施贴图操作,实现视觉传达效果理想的震后区域三维图形的虚拟重建。实验证明,所提方法对震后区域的三维图像虚拟重建结果精度高、视觉效果好。     

直流电阻率三维正演的代数多重网格方法

多重网格方法在求解由偏微分方程的边值问题离散所得线性系统时,具有非常高的计算效率.但常用的几何多重网格法在处理带跃变系数的偏微分方程时存在一定缺陷,限制了其应用.本文应用代数多重网格（AMG）方法求解三维直流电阻率法正演模拟形成的有限差分线性方程组,通过求解二次场的方法消除了总场中由点电源导致的奇异性,从而获得快速、精确的三维电阻率数值模拟.对两个存在大的电性差异的模型进行了模拟计算,以验证代数多重网格法的收敛效率.计算结果表明,与不完全Cholesky共轭梯度（ICCG）方法相比,代数多重网格方法具有更高的计算效率及稳定性.而且,随着三维网格节点数的增加,代数多重网格方法计算的高效性更加明显.     

用积分方程法研究源与勘探区之间的三维体对CSAMT观测曲线的影响

在可控源音频大地电磁法野外作业中,源和勘探区间的距离可达几公里,为了了解源和勘探区间的异常体对勘探区内异常响应的影响,我们用三维压缩积分方程法做了数值模拟研究.首先对勘探区目标体进行了数值模拟,发现在高频时,观测到的异常中心位于目标体的正上方,随着频率降低,出现异常中心向远离源的方向略有移动的现象,所以对三维异常体最好用3D软件来解释.然后,对源和勘探区间存在三维异常体的情况进行了数值模拟与分析讨论.模拟结果表明只有当三维异常体达到较大的规模时,才会对目标体上方的观测曲线造成影响,否则其电阻率的变化及埋深的变化对观测曲线的影响较小,可以忽略.当异常体在源方向有延伸时,观测曲线受到的干扰最大,沿垂直源布设的方向延伸时引起的干扰中等,垂直地面向下延伸引起的干扰最小.     

基于高斯牛顿法的频率域可控源电磁三维反演研究

三维反演解释是电磁法勘探发展的重要趋势,而如何提高三维反演的可靠性、稳定性和计算效率是算法开发者们目前的研究重点.本文实现了一种频率域可控源电磁(CSEM)三维反演算法.其中正演基于拟态有限体积法离散化,利用直接矩阵分解技术来求解大型线性系统方程,不仅准确、稳定,而且特别有利于含有大量发射场源位置的CSEM勘探情况;对目标函数的最优化采用高斯牛顿法(GN),具有近似二次的收敛性;使用预条件共轭梯度法(PCG)求解每次GN迭代所得到的法方程,避免了显式求解和存储灵敏度矩阵,减小了计算量.以上这些方法的结合应用,使得本文的三维反演算法准确、稳定且高效.通过陆地和海洋CSEM勘探场景中的典型理论模型的反演测试,验证了本文算法的有效性.     

多频可控源电磁法三维有理函数Krylov子空间模型降阶正演算法研究

本文采用有理函数Krylov子空间模型降阶算法实现了同时求解多频可控源电磁法三维正演响应的快速计算.首先采用基于Yee氏交错网格的拟态有限体积法实现控制方程的空间离散,将任意频率的电场响应表示为关于频率参数的传递函数.采用有理函数Krylov子空间算法求解该传递函数.针对构建m维有理函数Krylov子空间需要求解m次（几十到上百）关于有理函数极点和离散控制方程系数矩阵的线性方程组的问题,本文提出采用单个重复极点的有理函数Krylov子空间模型降阶算法,结合直接法求解器PARDISO,采用Gram-Schmidt方法,只需要1次系数矩阵分解和m次矩阵回代即可实现有理函数Krylov子空间的构建,极大地减少了计算量.针对最优化有理函数极点选取问题,本文根据传递函数的有理函数Krylov子空间投影算法的误差分析理论,引入关于单个重复极点的收敛率函数,通过求解有理函数的最大收敛率直接给出最优化的单个重复极点公式.最终实现了不同发射频率的可控源电磁法三维正演响应的快速计算.分别计算了典型层状模型多发射频率的CSAMT和海洋CSEM的正演响应,通过与解析解的对比验证了本文算法在多发射频率正演的计算精度和计算效率;并通过一个三维海洋CSEM勘探设计最优化发射频率和接收区域选取的例子进一步说明本文算法的优点.     

基于图像压缩的震损结构三维模型快速重建方法

针对基于数字图像的大尺度三维结构模型建模时间过长的问题,提出了利用图像压缩算法对震损结构三维模型快速重建的方法.首先,拍摄得到结构的原始图像;其次,利用主成分分析算法压缩图像;最后,基于处理后的图像对结构三维模型进行重建.为验证提出方法的有效性,对一个混凝土试块,受损的剪力墙试验模型和实际单体建筑进行图像采集,利用论文...     

Two-step interface and velocity inversion—— Study o e of the Tangshan seismic region

ｎｔｒｏｄｕｃｔｉｏｎＩｎｒｅｃｅｎｔｔｅｎｙｅａｒｓ,ｒｅｓｅａｒｃｈｅｓｏｎｔｈｅ３Ｄｃｒｕｓｔａｌｓｔｒｕｃｔｕｒｅａｎｄｔｅｃｔｏｎｉｃｓｈａｖｅｂｅｅｎｐａｉｄｍｏｒｅａｎｄｍｏｒｅａｔｅｎｔｉｏｎ,ｓｔｕｄｉｅｓｏｆｒｅｌｅｖａｎｔｍｅ...     

3D3C陆地反射PS转换波CMP成像影响因素分析

三维三分量(3D3C)陆地反射PS转换波共中心点(CMP)叠加成像方法,虽然抽道集简单,但是对实际资料处理结果往往不理想.尤其当反射界面为三维倾斜界面时,其成像质量较差.本文提出有三个主要因素影响其成像质量:第一,转换点离散.运用实例计算得出,转换点离散度随着纵横波速度比、偏移距和界面倾角的增大而增大.相同界面倾角,不同测线方位的转换点离散度不同,视倾角的绝对值越大离散度也越大;第二,道集内静校正量差异增大.CMP道集中,由于转换点离散使得转换点横向跨度较大,经倾斜界面反射转换的S波出射到近地表地层时的角度差异也较大,导致静校突出;第三,加大动校叠加复杂性.三维倾斜界面PS波CMP道集近炮检距时距方程可表示为双曲形式,但是曲线的顶点位置和动校速度同时随测线方位变化,使得CMP道集同相轴很难校平,动校叠加过程很复杂.