利用矩阵分解理论的双程波方程叠前深度偏移方法

叠前深度偏移理论及方法一直是地震数据成像中研究的热点问题.业界对单程波叠前深度偏移方法和逆时深度偏移开展了深入的研究,但对双程波方程波场深度延拓理论及成像方法的研究还鲜有报道.本文以地表记录的波场值为基础,利用单程波传播算子估计波场对深度的偏导数,为在深度域求解双程波方程提供充分的边界条件,并提出利用矩阵分解理论实现双程波方程的波场深度外推.通过对强速度变化介质中传播波场的计算,与传统的单程波偏移方法相比,本文提出的偏移方法计算的波场与常规有限差分技术计算的波场相一致,证明了本方法计算的准确性.通过对SEAM模型的成像,在相同的成像参数下,与传统的单程波偏移算法和逆时深度偏移算法方法相比,本文提出的偏移方法能够提供更少的虚假成像和更清晰的成像结果.本文所提偏移算法具有深度偏移和双程波偏移的双重特色,推动和发展了双程波叠前深度偏移的理论和实践.     

长江河口北槽弯道环流的涡度研究

基于ADCP走航观测得到长江河口北槽弯道附近3个横向断面(AD3、AD5和AD6)的流速资料,采用涡度方法,本文计算、分析了弯道环流与混合在垂直横向上的时空分布、影响因素及其重要性。3个横向断面上均存在由不规则界面分开的二层结构的横向环流。半拉格朗日余流的计算结果显示:(1)小潮期间,AD3、AD5和AD6断面呈现表层向海、底层向陆的纵向环状半拉格朗日余流;大潮期间,呈现表、底层均向海的纵向半拉格朗日余流;(2)小潮期间,AD3断面呈现表层向北导堤、底层向南导堤的横向环状半拉格朗日余流;大潮期间,AD3断面中间区域呈现表层向北导堤、底层向南导堤的、而断面两端区域则呈现表层向南导堤、底层向北导堤的横向环状半拉格朗日余流;(3)小、大潮期间,横向断面AD5和AD6均呈现表层向北导堤、底层向南导堤的横向环状半拉格朗日余流;(4)"纵向半拉格朗日余流"在–0.2—0.7m/s;横向半拉格朗日余流"在–0.15—0.2m/s;(5)纵向半拉格朗日余流在横向上有明显变化。对弯道环流的进一步分析表明:(1)斜压梯度、内部摩擦致混合和底部摩擦致混合这三项各自的纵向分量是驱动纵向环流形成的主要因素,"纵向动量的横向重新分布项"次之,离心力和地转的影响可忽略;(2)横向斜压梯度和内部摩擦致混合项是驱动横向环流形成的主要因素,离心力、地转和底部摩擦致混合次之;(3)横向环流可能通过"纵向动量的横向重新分布项"减弱纵向动量,从而可能减弱纵向环流。     

质心法雷达回波外推的原理及应用

介绍了多普勒雷达原始数据的格式,并通过计算机对两个相邻时次的雷达回波进行自动外推,实现了雷达回波外推的自动化和客观化.在2004年5月～2005年9月的应用中,这种方法对回波的移向、移速外推结果和实况基本吻合.     

空间目标激光测距轨道预报研究

为了实现地面卫星激光观测站对空间目标的轨道预报,设计了空间目标激光测距轨道预报软件。首先,利用空间目标的两行根数结合轨道外推算法,计算出空间目标在地心惯性坐标系下的坐标;其次,通过空间目标、地心和观测站三者之间的坐标系转换,计算出空间目标相对于站心坐标系下的位置;再次,根据空间目标相对地面激光观测站的方位角和俯仰角,画出空间目标的测站预报图;最后,以北斗卫星的测站预报实例对该软件进行验证,所得出的测站轨道预报图与同类商业软件的图结果一致。结果表明,该软件能够为地面观测站捕获空间目标,并同步跟踪其轨道运行提供测站轨道预报服务。     

一个浅水波方程的Semi—Lagrangian—半隐式格式试验方案

在正形投影坐标系下,半拉格朗日法与半隐式法相结合积分浅水波方程导出的模式,用一小时的时间步长,是无条件稳定。我们作了一个60小时的实例时间积分,计算稳定。半隐式积分格式导致求解一个赫姆霍茨方程,为比较其高低精度的优劣,我们提出在空间上用二阶精度和四阶精度解赫姆霍茨方程的超松驰法,同时,作了有无地形的试验。应该注意列,半拉格朗日法含有大量的插值,因而,在每一步时间步长里此积分格式比通常的欧拉法含有更多的计算,但是每一时间步长的计算量的增加却从使用很大的时间步长得到了更多的补偿。     

基于稀疏离散τ-p变换的叠后地震道外推

这里研究的是一种基于稀疏离散τ-p变换的地震道缺失重建方法。该方法根据叠后零偏移距剖面,在局部时窗内可以看作是一系列线性同相轴的组合,使用稀疏离散τ-p和预条件双共轭梯度算法进行地震道外推,使空间方向的缺失道得到恢复。τ-p变换由于信息不足,有限的孔径和离散等因素,导致在τ-p域的结果存在假象、不准确,而在τ-p域进行稀疏对于重建缺失的信息是非常有用和必要的。同时,还可使用迭代的预条件双共轭梯度算法进行计算。     

测井约束下的神经网络地震储层参数反演

针对目前大多数地震反演方法基于褶积模型及相应假设条件的现状,本文提出了在假定储层参数与地震资料之间存在某种非线性映射的情况下,以测井作为约束条件,利用Ｅｌｍａｎ局部递归神经网络法,建立地震数据到地层参数的映射模型,将地震数据直接由已知井点向外递推地层速度、孔隙度和泥质含量的方法。文中给出了方法原理、处理流程及应用实例。     

道真县干旱强度指数方差分析周期动态计算程序设计及预报

针对方差分析周期叠加外推法在长期天气预报中应用较多,而周期寻找和周期元值提纯相当繁琐的问题,设计了一个小程序,对道真干旱强度指数进行周期分析,在道真夏旱预报中应用效果较好.此方法从传统的手工计算转向计算机运算,不但提高了计算速度,而且保证了计算结果的准确率.计算程序不受样本长度限制,随时可以增加样本长度,补充最新资料进行周期元值提纯,提高拟合度.     

利用中国区域电离层数据拟合 Klobuchar参数

对于单频G PS用户而言,电离层延迟是最重要的误差来源之一。G PS系统使用Klobuchar模型对电离层延迟进行改正,其改正数从根据经验模型和历史数据得到的370组常数选取。在缺少全球观测数据的情况下,仅利用中国区域观测数据拟合电离层模型参数,存在区域外精度下降,参数超限等问题。针对这些问题,提出了利用经验模型外推进行参数拟合的方法。相比广播Klobuchar参数,在太阳活动高峰年RMS误差减小接近7 TECU ,正常年减小1.5 T EC U ,平静年精度也略有提升。2001年至2012年的拟合参数不存在参数超限现象,且量化后或者使用预报的电离层数据进行拟合,精度下降很小,可以用于预报全球电离层模型参数。