单程波方程偏移算子高阶辛格式

针对地震偏移算法中单程波算子中e指数项的近似，提出一种高阶辛格式算法. 与非辛格式的近似、一阶辛格式近似和两步法进行数值对比表明，文中提出的高阶辛格式具有更高的精度. 脉冲响应的计算也表明文中的近似方法具有降低频散效应和适应横向变速的能力. 用Marmousi模型验证了对复杂构造的成像的能力，实际资料计算表明本文方法适合于较大深度的稳定准确成像.     

局部角度域波传播步进算法研究

本文从非均匀介质中波动方程出发,提出了基于一般标架的相空间（局部角度域）波传播的步进算法. 该方法在构造单程波的步进算法时,在选择标架或正交基等方面有更大的自由度. 我们以不随频率及深度变化且具有变尺度特性的Gabor_Daubechies紧标架为例,给出了单程波传输算子的具体形式及相应的波场步进算法;详细讨论了基于Gabor_Daubechies标架的传输算子的高频渐近展开问题,得出了在高频、小传输步长条件下传输算子的近似解析表达式,并给出使用条件. 通过模型算例,比较了精确传输算子与高频近似传输算子用于非均匀介质中波传播的结果,说明在一定条件下由两者得出的波场几乎是相同的.     

共炮检距道集波动方程保幅叠前深度偏移方法

本文提出了一种基于双平方根算子的共炮检距道集波动方程保幅叠前深度偏移方法,将振幅误差补偿作为偏移的一部分与“运动学偏移”一起在偏移过程中实现.其基本内容包括：(1)从保幅的单平方根算子方程出发,推导出由双平方根算子定义的保幅单程波方程;(2)根据地震波摄动理论把速度场分裂为层内常速背景和变速扰动,分别在频率－波数域和频率－空间域求得波场深度延拓的偏移时移量及振幅校正系数,从而得到最终的DSR保幅波场延拓算子;(3)在高频假设条件下,把DSR保幅波场延拓公式中的积分运算进行稳相近似,得到保幅波场延拓的相移公式.理论分析和模型数值试验表明,该方法不但可以使散射能量聚焦、归位,提高成像精度;而且可以输出正确反映地下反射系数的振幅信息,为后续的地震属性分析（如AVO/AVA）提供更真实的地震信息.     

一种双向测距与时间同步系统的设计与分析

在研究双向单程伪距测量原理的基础上,设计了双向测距与时间同步（DRTS）终端系统总体构架,阐述了在系统中使用的技术,并搭建了基于DSP＋FPGA的双向测距与时间同步系统软硬件平台。实验结果表明,此系统的码速率为5MHz、中心频率为15MI-Iz时,测距和时间同步的分辨率可达0．15cm和5ps（@1S）,采用不同频率源时测距和时间同步的精度分别为1．038m和3．46ns,采用相同频率源时分别为0．28cm和9．43ps（参考频率稳定度1×10^-10／d量级）。与国外同类产品相比具有测量精度优势,但考虑通用性,此系统的硬件仍需进一步优化,软件上需要做到码速率可调。     

逆时偏移速度误差影响分析

逆时偏移作为一种先进的地震偏移成像方法,其成像结果的好坏取决于很多因素,其中最关键的是速度.速度越准确,成像效果就越好.但在实践当中,速度往往是未知的,只能通过速度反演等各种手段来估计近似速度.由此导致逆时偏移成像结果存在误差,从而降低后续地震解释的精度.那么速度误差对逆时偏移成像的影响到底有多大呢?实际当中一般估计出来的近似速度可能在模型每一个地方都不同,很难直接进行逆时偏移速度误差系统分析.为了简化,本文采用等效速度误差——速度模型整体平均误差,进行逆时偏移成像分析.逆时偏移成像采用互相关成像条件,低频噪声压制采用振幅补偿拉普拉斯滤波方法,源波场正推和接收波场逆推采用波动方程一阶应力.速度形式及交错网格有限差分方法.首先对比分析不同速度误差对逆时偏移成像结果的影响,并和单程波波动方程偏移进行对比分析;然后分析速度误差对逆时偏移和单程波波动方程偏移成像位置的影响;最后基于偏移速度分析方法估计得到的近似速度模型,进行逆时偏移和单程波波动方程偏移成像试验.结果 表明:在不同速度误差情况下,逆时偏移成像在同相轴连续性和能量聚焦等方面要好于单程波波动方程偏移;速度误差对逆时偏移和单程波波动方程偏移成像位置的影响程度相当.本文研究成果对逆时偏移在实践中的应用具有一定的参考价值.     

基于单程波偏移算子的地表相关多次波成像

在常规地震资料处理中,多次反射波被视为噪声并从地震数据中去除,以免在之后的地震资料解释中造成误解.而事实上,多次波也是地震信号,是照明波场的一部分,能够对地下构造成像的精度做出贡献.本文分析了多次波在传统单程波叠前深度偏移中产生构造假象的机制和表现,为实现基于单程波偏移算子的多次波成像,修改了单程波叠前深度偏移的边界条件,即将输入的震源波场用包含多次波的记录来替代,输入的记录波场用预测出的表层相关多次波来替代,实现了基于单程波偏移算子的地表相关多次波成像,并从理论上给出了其成像依据.通过基于二范式最小能量差原则求取的匹配因子,将多次波成像结果与一次波成像结果进行匹配叠加,应用多次波成像来弥补一次波成像的不足.简单模型验证了基于单程波偏移算子的多次波成像方法的有效性,最后对Sigsbee2B模型进行了一次波与多次波联合成像试算,盐边界高陡构造成像质量得到了明显改善.     

一阶弹性波方程数值模拟中的混合吸收边界条件

Liu和Sen（2010和2012）在地震波场数值模拟中提出一种混合吸收边界条件. 该方法具有计算量小、容易实现及吸收效果好等优点. 但现有的混合吸收边界条件是针对二阶位移-应力方程设计的,存在稳定性问题. 本文首先推导了两种速度-应力单程波方程：二阶Higdon单程波方程和一阶Higdon单程波方程. 进而提出基于一阶弹性波方程的混合吸收边界条件方法. 在内部区域和边界之间引入一个过渡区域,通过单程波与双程波方程平滑过渡来消除人工边界反射. 为了改善混合吸收边界条件的吸收效果和稳定性,我们采用了能同时吸收纵、横波反射的一阶单程波方程和与变量位置有关的加权系数. 为了验证混合吸收边界条件的有效性,将其与常规分裂完全匹配层（PML）方法进行了比较. 数值模拟结果表明,与PML边界条件相比,混合吸收边界条件在耗用更小计算时间和存储量的前提下,可以获得更好的吸收效果. 另外,本文提出的两种混合吸收边界条件中,混合一阶Higdon吸收边界条件具有更好的稳定性.     

基于波场垂向和水平方向外推的高陡构造偏移

常规的单程波波动方程偏移成像方法对大角度的高陡构造偏移成像存在内在的限制.根据波动方程在各个空间方向的数学特性和高陡构造反射地震波的传播特征,通过把地震波分解为垂向的上下行波、水平方向的前后行波和左右行波,提出基于波场垂向外推和水平方向外推相结合的单程波波动方程高陡构造偏移成像方法,即用波场垂向外推的单程波波动方程偏移成像方法解决中低角度平缓构造的偏移成像,用波场水平方向外推的单程波波动方程偏移成像方法解决中高角度陡倾构造的偏移成像.这种基于波场垂向和水平方向外推相结合的高陡构造偏移成像方法是常规单程波波动方程叠前深度偏移成像方法的补充和改进,它相对基于全波方程的逆时偏移具有计算效率上的优势.     

常用叠前深度偏移方法特点分析与实例对比

目前常用叠前深度偏移主要有四种:基于射线理论的Kirchhoff偏移和高斯束偏移以及基于波动方程理论的单程波偏移和逆时偏移.要想在实际应用中取得好的成像效果,必须选取合理、有效的叠前深度偏移算法.本文首先简单阐述了这四种常用叠前深度偏移方法的基本原理,对比分析了四种方法在陡倾角、速度依赖性、起伏地表、非规则数据、低信噪比数据等方面的优劣性以及计算效率和计算精度的差异,然后利用模型试算和实际资料处理分析说明这四种叠前深度偏移方法的成像效果,最后给出不同情况下叠前深度偏移方法的选取建议.     

嫦娥二号卫星X波段测控体制试验定位分析

2010-10-03～2010-10-04,在嫦娥二号任务实时观测阶段进行了X波段测控体制试验.从定位归算角度对试验期间的S波段数传信号(S1频点)测量时延(S1时延)和X波段单程测距差分(DOR)信号测量时延(DOR时延)进行了比对分析.初步结果表明,DOR时延的实测资料精度与链路计算理论估值基本相符.而且由于信号的...