均匀介质中的F-K反偏移:基本概念、基本公式及其在非均匀介质中的应用

从经典的常速度F-K偏移成像理论公式出发,通过由垂直波数到角频率映射的途径首先建立了常速度F-K反偏移的基本理论公式和基本实现算法,然后再将其用于解决非均匀介质中的反偏移问题。与同样条件下的Kirchhoff型反偏移理论相比,在建立常速度F-K反偏移理论时没有引入任何近似。因此,所提出的是一个在均匀介质中严格精确的反偏移理论。对于一般条件下的非均匀介质,利用了在现代反射地震偏移理论研究中常用的局部化处理方法和相移加插值(PSPI)偏移的基本思想。具体地讲,在处理非均匀介质中的反偏移问题时假设偏移场的形成完全由局部薄板(层)中的速度结构决定。此外,还假设偏移场是速度的连续函数。因此,可以在局部薄板(层)中利用关于垂直坐标的Fourier变换和均匀介质中的频散关系,以及一组均匀参考速度来近似地构造出在不均匀速度模型中任意一个给定深度上的反偏移场。在算法上,在均匀介质中的F-K反偏移是一个一步过程,而在非均匀介质中基于PSPI的F-K反偏移是一个多步递归过程。为了实现反偏移必须要从速度模型的最大深度开始逐层地进行插值和反向相移计算,直至到达地表为止。     

基于空间特征的二维矢量地图可逆水印算法原理

可逆水印技术能够在精确提取水印信息的同时,完全无损恢复原始载体数据,更适用于矢量地图版权保护.提出一种基于空间特征的二维矢量地图可逆水印算法.该算法在保持矢量地图特征点相对位置的前提下,根据密钥对特征点位置进行非线性置乱,避免非授权用户对矢量地图高精度数据的使用,并在置乱后的特征,占'中嵌入水印信息,能够在提取水印信息的同时无损恢复原始矢量数据.实验结果表明该水印方案严格可逆,具有较好的隐蔽性.     

2.5维地震波数值模拟评述：声波模型

本文的目的是对基于声波模型的2.5维地震波数值模拟工作进行评述,以便能够找出其存在的问题和解决这些问题的可能途径,根据定义,2.5维问题是三维问题中的一种特殊情况,其特点是:(1)介质参数沿走向保持为常数;(2)场源具有球对称性;(3)场源和接收点均位于垂直于走向的直线上,与三维数值模拟问题不同,2.5维数值模拟问题分为两部分:(1)在垂直于走向的平面内用数值方法解相应的微分方程,这在实质上是二维问题;(2)采用积分变换或其他方法处理来自于计算平面外的影响,这实际上是将一个特殊的三维问题转化成为了无限多个(在离散情况下是有限多个)二维问题的叠加,与二维模型相比,2.5维模型能得到计算平面内的精确地震波振幅信息,鉴于声波模型是反射地震偏移成像理论和应用研究中的基本数据模型,所以对2.5维声波数值模拟的研究具有重要的意义,根据对计算平面外传播效应的处理方式可以将到目前为止提出的2.5维声波数值模拟方法分为四类:(1)几何射线法;(2)滤波校正法;(3)Fourier变换法;(4)近似波动方程法.其中,几何射线法具有直观、快速的特点,但是在焦散区内失效,滤波校正法只在均匀介质条件下严格成立,在一般条件下只是一种精度难以估计的近似,Fourier变换法是一种经典方法,其研究程度已经相当深入,该方法的基本思想是通过沿走向的Fourier变换将2.5维问题转化为有限多个二维问题,从而,对反变换的数值实现直接影响到该方法的精度和效率,近似波动方程法的宗旨是针对2.5维波动问题建立专门的波动方程,与Fourier变换法相比,近似波动方程法等同于一个二维数值模拟,因此可以大大地降低计算量,但是,根据笔者所掌握的资料,到目前为止提出的几个近似波动方程不是具有很大的振幅误差,就是难以进行数值计算,因此,有必要对近似波动方程的形式进行进一步的研究.     

几种单程波动方程解法的对比

在均匀介质条件下, 比较了ω-k域内Stolt算法和Gazdag相移法、t-x域内部分分式算法、混合域内SSF算法和FFD算法的计算流程, 零炮间距叠后时间偏移的单位脉冲响应、脉冲振幅变化、中心地震道数据的子波振幅变化和脉冲波场频谱特征等体现了波场动力学特性的信息。分析表明, 在均匀介质条件下, ω-k域方法在脉冲响应、振幅变化和频谱特征方面能得到较好的结果, 但该方法的子波分辨率相对其他方法较低, 无法应用于横向速度变化介质问题, 所以复杂地质问题应采用FFD法进行处理。同时以Stolt算法为例具体分析了子波振幅、子波频谱在偏移过程中的变化。     

基于完全三叉树堆排序的波前扩展有限差分地震波走时快速算法

波前扩展有限差分地震波走时算法具有物理意义明确、因果稳定性强的特点,但每次波前扩展都要寻找波前面上的最小走时点。当计算网格点数较多,特别是涉及到三维走时计算时,寻找波前面上的最小走时点是一项十分耗时的工作。研究发现,波前扩展有限差分地震波走时算法的波前点具有两个突出特点:①波前点更新十分频繁,通常每次取出波前最小走时点后都要插入若干新的波前点;②新计算出的波前点的走时通常比较大。数据结构中的二叉树堆排序方法可以提高寻找波前面上最小走时点的效率,根据特点①,在原始二叉树堆排序方法的基础上,优化了插入新波前点和移除波前面上最小走时点的流程,实际计算结果表明,与原始的二叉树堆排序方法相比,改进后的二叉树堆排序方法可以提高大约20%的计算效率。根据特点②,将原始的二叉树堆排序方法推广到多叉树,实际计算结果表明,完全三叉树堆排序方法优于二叉树和四叉树堆排序方法,可以再提高5%的计算效率。     

二维起伏地表直流电场插值法数值模拟

本文提出了一种采用插值法实现边界条件的有限差分数值模拟新方法, 地表函数通过三次样条插值能很好的拟合真实地表形状。其地表边界条件的实现思想是在地表法线方向上寻找一个与地表点位于同一等位线上的点, 该点通过周围点插值来建立离散方程。采用不规则四边形距离加权插值法使插值形式统一化。在水平地表情况下做的数值解与解析解的对比分析表明: 该方法的计算精度适合要求, 对山谷模型的实例计算也得到了视电阻率与地形成镜像对称的结果。     

三维直流电场标量拟解析近似法数值模拟精度分析与评价

拟解析近似方法是一种可以处理强散射或大扰动问题的求解积分方程的近似方法,该方法可以避免用传统数值方法解决问题时,所遇到的大型矩阵或大型代数方程组问题。根据孙建国在文章[14]中的求解异常电场积分方程的标量拟解析近似理论公式,对标量拟解析近似进行研究。对均匀场中的异常球体标量拟解析近似解进行计算,并与理论解进行了对比,验证了其精度,证明拟解析近似方法求解直流电场积分方程是可行的。通过对标量拟解析近似方法在直流电场数值模拟的研究,为三维直流电场快速正反演模拟打下了基础。     

P-SV波和SS波反射系数之间的近似关系

本文通过对Aki&Richards给出的P SV波反射系数的一阶近似 ,给出了一个新的描述P SV波和SS波反射系数之间关系的近似公式。在 4个不同类型的含油气砂岩模型上的定量计算表明 :新的近似关系能够反映SS波阻抗和P SV波反射系数之间的内在关系 ,而且形式简单、近似效果好、适用范围广 ,适合于SS波阻抗的定量反演和理论研究 ;如果再结合PP波AVO信息 ,还可以简单、准确地进行多参数反演。另外 ,在 2 0°以内 ,利用该近似关系 ,纯SS波零偏移距剖面和P SV转换波剖面之间可以相互转换。     

一种利用TM图像自动提取洪积扇的方法

在分析洪积扇与其他地物光谱特征差异的基础上,针对TM图像的红光波段与近红外波段的比值能增大洪积扇与其他地物间光谱差异,以及地形阴影在蓝、绿光波段亮度值降低速率差异较大的特征,基于比值运算和差值运算,构建了洪积扇自动提取模型.利用该模型,首先结合阈值算法将洪积扇从其他地物中提取出来,然后采用数学形态学膨胀和腐蚀算法进行空...     

标量拟解析近似方法模拟立方体的异常电场

拟解析近似方法是一种求解积分方程的一种近似方法,它可以处理强散射或者大扰动的电磁散射问题,在计算过程中避免了传统微分数值方法解决问题时所遇到的大型矩阵或大型代数方程组的求解。孙建国[5]将其引入直流电场的积分方程中,并给出了求解异常电场积分方程的标量拟解析近似公式。在以前的研究中,已经验证了均匀场中异常球体的拟解析近似解的精度,这里对均匀场中的立方体异常体进行数值模拟,得到了直流电场中异常立方体模型的标量拟解析近似解。由于复杂地电模型可以用立方体的组合进行模拟,因此对立方体异常电场拟解析近似解的研究,为三维直流电场中复杂地电模型的快速正反演模拟打下了基础。