交错网格中基于波数域插值的波场分离方法研究（英文）

应用多分量地震资料进行成像时通常需要先做波场分离,然后再对分离的波型进行成像。其中,波场分离可以在空间域或波数域实现。然而,由于用交错网格有限差分进行弹性波场数值模拟时,用来进行波数域波场分离的质点振动速度分量定义在不同网格节点上,本文提出了利用波数域插值方法来估算同一网格节点所需质点振动速度值;进而给出了先进行波数域插值后进行波场分离的波数域保幅波场分离方案。数值实验结果表明波数域插值方法具有较高的插值精度且保幅波场分离方法具有较好的保幅性,将本文方法进一步应用于弹性波逆时偏移可以获得保幅性较好的成像结果且对存在一定程度速度误差情况具有较好的适应性。     

不同截止高度角多模GNSS组合单点定位性能分析

讨论多模GNSS时空统一方法,建立多模单点定位数学模型。采用8个测站的实测数据,对比单GPS与GPS/GLONASS/BDS/Galileo四大GNSS组合多模系统在截止高度角分别为15°、30°、45°时的单点定位性能。结果表明,15°时,多模较单GPS定位精度有一定改善,但效果不明显| 30°、45°时,单GPS不能实现实时定位,多模则能提供稳定可靠的实时定位结果。在北斗卫星导航系统可见卫星数较多的亚太区域,95%以上历元可用,水平方向偏差RMS<8 m,高程方向偏差RMS<20 m,对城市峡谷和密林中导航具有一定的参考价值。     

台风莫拉克登陆期间的中尺度波动特征分析

利用波作用理论对台风莫拉克登陆期间的降水进行诊断分析。结果表明:波作用密度异常能在一定程度上指示暴雨雨区发展移动,其异常值的空间分布能够反映雨区上空动力场和热力场的典型垂直结构特征。台风内中尺度波动与暴雨落区宏观上具有一定联系。为了详细研究台风内部的波动特征,利用WRF(Weather Research and Forecasting model)模式模拟的高分辨率资料对台风登陆过程中波动特征进行分析,低波数波动的传播主导强对流的非对称分布,2波在登陆初期对对流分布起着关键作用,中尺度波动中同时存在涡旋罗斯贝波以及重力惯性波的特征,对登陆期间涡旋混合的现象有重要作用。     

基于理论地震图方法的波导结构对 地下核爆炸Lg波传播特性研究

基于频率-波数域算法的理论地震波形图方法, 可以数值模拟频率达到10 Hz、 震中距达1000 km的区域理论地震波形图. 该算法适用于计算大量分层地壳结构中激发的导波Lg波. 本文在前苏联东哈萨克斯坦地下核试验场至我国乌鲁木齐台站间的地球介质速度模型中, 引入速度梯度结构、 速度扰动分布的薄叠加层结构、 降低Q值结构以及速度扰动与Q值变化的综合结构来模拟实际地壳波导结构, 较好地模拟出东哈萨克斯坦地下核爆炸地震在乌鲁木齐台站记录的宽频带地震波形图, 模拟出完整的Lg波序列, 该序列符合Lg波能量分布特征, 且能够解释Lg波波尾的特征. 结果表明, Lg波的形状和峰值结构均依赖于地壳的不同波导结构.      

六种遥感海表温度产品的比对分析

海表温度(sea surface temperature, SST)产品是研究全球海洋大气系统的重要数据源, 在海洋相关领域的研究和应用方面具有重要价值。对2007年1和7月的六种不同SST产品(AVHRR-only, AMSR+AVHRR, NCODA, RSS, RTG-HR和OSTIA)在南大洋阿古拉斯回流(Agulhas Return Current, ARC)与绕极环流交汇区域的产品特性进行了比对统计分析, 包括SST分析、SST梯度分析、统计参数分析和波数谱分析。分析结果表明产品之间SST时空变化分布的整体趋势一致, RTG-HR在空间分布上过于平滑, OSTIA解析的大洋锋面最弱, RSS包含噪声较多, NCODA相对其他产品存在较大偏差。发现在AVHRR数据的覆盖率较好的情况下与以红外数据构建的AVHRR-only相比, AMSR数据并不能提供更多的SST信息, 反而会降低产品的空间解析能力。     

基于卫星海面高度资料的太平洋海区涡旋能谱分析

在空间均匀和各向同性的假定下,地转湍流理论认为涡旋动能随尺度分布遵循k–n定律(k为波数)。但实际海洋受边界、地形、层结等,具有明显的非均匀、各向异性特征。鉴于此,我们基于近30a卫星高度计资料,分别计算了热带、副热带、中高纬度等海洋涡旋强度不同区域的海面高度异常(SSHA)波数谱,进而利用线性回归拟合方法估算出中尺度波段上SSHA波数谱的斜率,并与地转湍流理论预测进行了对比。研究结果显示:SSHA波数谱从赤道到中高纬度逐渐变陡,其斜率由–4减到–5,基本符合赤道线性波动理论和准地转湍流理论的预测。SSHA波数谱斜率存在纬向与经向差异,例如在赤道地区,纬向谱比相应的经向谱陡;而在南极绕极流区域,经向谱斜率大于纬向谱斜率。SSHA波数谱斜率的各向异性表明海洋中尺度运动受β效应影响,具有明显的经向和纬向差异。以上结果表明,海洋中尺度运动介于准二维和三维之间,不能用一个全球普适的湍流理论模型来描述。     

Chebyshev低通滤波器在位场数据波数域导数换算中的应用

导数换算作为常用的位场数据处理手段具有重要的物理意义。由于实际测量误差和噪声的影响,波数域导数换算存在计算过程不稳定、求导结果精度差的问题。为了降低噪声干扰,提高导数计算精度,提出了在波数域常规导数算子基础上附加Chebyshev低通滤波器的任意阶导数换算方法。通过分析该滤波器的滤波特性,并结合位场异常径向平均功率谱曲线特征,确定了Chebyshev低通滤波法的滤波参数,减少了人为因素对求导结果的干扰。二度体和三度体模型试验表明,与其他3种方法(常规FFT(fast Fourier transformation)法、向上延拓法和ISVD(integrated second vertical derivative)法)相比,该方法计算的导数与理论值的均方根误差最小,求导结果精度较高。在虎林盆地布格重力异常数据处理中,利用该方法计算的垂向二阶导数受噪声干扰小,结果可靠性较高,依此划分出的8条大断裂和11条小断裂在以往研究中均得到证实。     

二维Wiener滤波器及在物探中应用的一些问题

本文结合二维随机场函数z（x,y）及其统计特征、平稳性和各态遍历性的概念,导出了二维Wiener滤波器的两个波数响应H_opt（u,v）和 _opt（u,v）,并详细地讨论了导出它们的条件。为了正确地应用Wiener滤波技术,我们提供了检验物探观测数据的平稳性和各态遍历性的简单而适用的方法。最后还指出,当应用Wiener滤波的不相关技术时,将在区域异常的谱函数中引入一个非线性畸变因子,这个区域异常是从具有局部干扰的观测数据中提取的。于是,所提取的区域异常必然会产生某种程度的失真。     

上海地震台阵响应图及其计算方法

地震台阵的响应特性是衡量台阵对地震监测能力的重要标志。根据台阵响应特性，可分析台阵子台布局的合理性和波数分辨率。本文采用简单延时相加方法，并用不同频率下的波数k为上海地震台阵计算并绘制二维与三维的波数响应图。结果表明，该台阵响应图具有显著的主瓣峰，而且边瓣少且峰值低，具有较高的波数分辨率。