基于二维三次卷积插值算法的辛几何射线追踪

射线追踪是地震波走时层析成像的基础,射线空间位置的准确性及射线走时的精度决定了层析成像的可靠性.本文根据哈密尔顿系统可以有效提高程函方程解稳定性的特性,采用辛几何算法（SAM-Symplectic Algorithm Method）及二维三次卷积插值技术进行地震波射线追踪.由于采用了SAM算法,保证了地震波波前精度,提高了射线空间位置的准确性.数值模拟结果表明SAM既能保证哈密尔顿系统的稳定性又具有运算速度快的特点,提高了射线追踪的计算精度.     

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

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

基于SVD和BRD的二维核磁共振测井正则化反演算法研究

二维核磁共振能从"弛豫-扩散"两个维度上展现流体性质,在识别稠油储层方面具有理论优势,是当前核磁共振测井技术的研究热点.本文深入研究二维核磁共振测井原理,系统分析CPMG-DE脉冲序列测量扩散系数与弛豫时间的方法,结合核磁共振二维谱数理模型,提出一种基于SVD和BRD的正则化反演算法.该算法通过SVD压缩数据,采用带非负约束的Tikhonov正则化方法求解流体"弛豫-扩散"分布,并基于BRD算法迭代确定最佳正则化因子.模拟实验与数值分析表明,该算法无需先验信息、运算效率高、相对误差小,在原始数据信噪比低至50时,仍可有效获取流体(T₂,D)二维分布.在二维核磁共振测井数据实时解释应用中,该方法较传统反演算法(如TSVD)具有较大优势.同时,在自主研发的核磁共振测井仪测量CuSO₄溶液(T₂,D)分布的实验显示,本文设计算法对弛豫时间和扩散系数的反演误差分别仅为2%和4%,较TSVD算法有较大改善.     

北京地区主要活断层活动参数的重力反演

采用1991~1999年首都圈重力测网的重力测量资料,根据地质调查结果确定断层初始运动模型,利用稳健—贝叶斯最小二乘方法反演北京地区4条活断层的活动参数(走滑、倾滑、引张)。反演结果表明:在1991~1999年间,北京地区活断层的运动性质均以走滑和垂直升降为主,兼有伸展特征,但运动速率都不大,约为2 mm/a。北京地区北西走向断层反演的运动特征为左旋走滑正断层,北东走向断层反演的运动特征为右旋走滑正断层。     

看得见的“沧海桑田”

沧海桑田,一直都用以形容世间万物的极大变迁,这种变化往往需要数千,甚至数万年的漫长岁月。可是今日当你立足于黄河三角洲时,沧海桑田的变化、海陆分界的变迁,一览无余。黄河三角洲是以黄河入海河道为中心、向着大海辐射的巨大扇形。     

圆柱体的拟合与质量检测

为检测空间圆柱体的位置和形态变化,利用无反射电子全站仪观测圆柱体外表面上均匀分布的若干点位,将不同测站上观测的点位坐标通过公共点归算到同一坐标系下,根据观测点的三维坐标,选择合适的方法拟合出圆柱体的空间位置参数和形状参数。文章介绍了一种基于最小二乘原理的空间圆柱体拟合方法,详细阐述了圆柱体质量检测的原理与计算过程。文中算例验证了该圆柱体拟合方法的正确性与可靠性,以及质量检测方法的可行性。     

总体最小二乘平差方法在GPS高程拟合中的应用研究

采用最小二乘(LS)进行GPS高程拟合参数估计未考虑系数矩阵误差,尝试采用总体最小二乘(TLS)平差方法进行参数估计。利用本文提出的基于TLS平差的粗差探测方法进行粗差剔除的基础上,对TLS平差方法在GPS高程拟合中求解的参数及其精度进行了分析,通过与LS的对比表明,混合总体最小二乘的拟合结果最为合理。     

北斗非GEO卫星信号码并行捕获

2012年12月,北斗二代的接口控制文件正式公布。北斗B1I频点的非GEO卫星上调制20bit的NH码,因此每1ms信号可能存在比特跳变。借鉴传统GPS并行码相位捕获算法,本文研究一种基于补零FFT的并行码相位算法。该算法能够实现1-20ms相干积分,利用非相干累加提高信噪比,并采用一种对信号强度和采样率不敏感的判决变量。仿真表明,通过相干积分和非相干累加,该算法使用400ms数据可以捕获24dBHz的弱信号。     

数字化形变仪器校准

依据最新形变观测规范要求,结合数字化形变仪器日常维护的经验,分析校准工作中常见故障与校准结果数据表象特征的关联关系及注意事项,提出"二次校准"法和"快、准、简"3字校准要诀,建立仪器校准新流程,以供形变仪器校准参考。     

泉州盆地地震效应特征的一维等效线性和二维非线性分析的比较

采用等效线性动粘弹性模型描述土的动力非线性特性,基于一维等效线性波传法,对泉州盆地地震效应进行了分析;同时,采用修正Martin-Seed-Davidenkov动粘弹塑性模型描述土的动力非线性特性,对泉州盆地非线性地震效应进行了大尺度二维精细化有限元分析,研究了地形地貌和土层横向不均匀性对地震效应的影响。将两种分析结果进行对比,结果表明：①随着基岩输入地震动强度增大,地表峰值加速度PGA放大效应总体呈现减小趋势,中震与小震、大震与小震的地表PGA放大系数之比依次为0.83~0.99、0.72~0.97;②该盆地Ⅲ类场地处,基岩、地表起伏不大,且土层横向分布较均匀,两种方法计算得到的地震效应特征类似;基岩或地表起伏剧烈、土层横向分布明显不均匀的Ⅱ类场地上,二维非线性分析给出的地表PGA放大系数明显大于一维等效线性结果,两种方法得到的地表加速度反应谱及PGA随土层深度的变化特征存在显著差异,二维非线性分析给出的地表加速度反应谱大多呈现双峰甚至多峰现象,且PGA在土层特定深度处存在聚集效应,使PGA随土层深度的变化呈现非单调性。