气候变化对沈阳冬季供暖期的影响

利用1961-2010年沈阳市的年平均气温和供暖期及供暖期前一个月、供暖期后一个月的逐日平均气温资料,采用线性趋势分析方法及度日法,分析近50 a沈阳市年平均气温和供暖期平均气温变化的总体趋势,并对供暖初日、供暖终日和供暖期长度及供暖强度的变化特征进行分析。结果表明：近50 a沈阳市年平均气温与供暖期平均气温均呈上升趋势,近50 a年平均气温上升明显,通过了显著性检验;气候变暖使得供暖初日推迟;供暖终日提前;供暖期长度缩短;气候变暖使实际供暖强度逐渐减弱。     

一种基于统一引力场表示理论的外部扰动场赋值建模方法

提出一种新的外部扰动场赋值的建模方法。该方法将场赋值问题纳入统一引力场表示理论的框架,得到的模式具有结构简单、奇性弱的核,并顾及了地形效应,适合低空引力速算。     

合肥市近50年最高气温变化特征

利用合肥市1953～2005年逐日最高气温资料,采用滑动t检验等统计方法,分析了合肥市平均最高气温、高温日数和高温出现初终日的变化特征.结果表明:合肥市近50年来平均最高气温呈"上升-下降-上升"的趋势,并且平均最高气温在1968和1988年附近发生了突变.合肥市高温日数也呈"增多-减少-增多"的趋势.但年累计高温日数没有很明显的突变;高温日5～9月都有出现,其中7月出现的次数最多,占总日数的49.7%,8月次之.占40.5%.合肥市高温持续的时间经历了由长变短再变长的过程,20世纪80年初期最短.     

总体最小二乘法在初轨计算中的应用

在分析了UVM2方法条件方程特点的基础上,将总体最小二乘法(TLS)引入轨道计算中,避免了原算法中对活力公式的线性化处理．仿真和实测计算结果表明,将TLS应用于UVM2是有效的．     

论初轨计算的最佳精度及二重解

从最优估计的角度,研究短弧测角资料初轨计算的最佳精度及准最佳精度算法．按照统计理论,对一般的参数估计问题,其估计误差的方差存在Cram’e4-Rao下界,这就是初轨理论的最佳精度;通过大量的模拟计算表明,文［1］中提出参考矢量法,其精度已充分接近C·R下界,且模拟计算轨道统计精度与观测误差成正比,可以认为是一种初轨计算的准最佳精度算法;对于小偏心率情况,发现一个重要的事实：初轨计算有二重解,更准确地说,在M0的左右两个半平面内有真假两个解,偏心率愈小,平近点角M0愈接近±π／2,给出的假解的概率愈大．针对这一情况,提出一种算法,可以给出2个解,其中必定有一个是真解,它有较高的精度.     

基于S变换和瞬时旅行时属性的自动初至拾取方法研究

地震波走时拾取是地震资料前期处理过程中的一个重要步骤,同时也是后期层析成像方法中计算速度的重要依据。本文提出了一种新的自动初至拾取的方法,利用地震信号的瞬时旅行时属性进行地震波自动初至拾取,瞬时旅行时属性的计算包括两个独立但相关联的部分,首先通过S变换将地震信号转换到时频域,然后利用映射算子将时频域表征映像回时间域得到瞬时旅行时,从而实现自动拾取初至。理论模型和实际地震资料处理的结果表明,该方法可以有效地拾取地震波初至时间。     

山西近50a初霜冻的气候变化特征及其周期分析

初霜冻13的变化对农业生产的布局和耕作模式具有一定的影响。利用山西62个气象观测站点1961－2010年的逐年初霜冻日及逐日最低地温资料,应用统计学方法分析了山西初霜冻日的变化特征,结果表明：1）山西近50a平均初霜冻日的时间分布总体为推后趋势,且具有明显的年代际特征,20世纪70年代、90年代初霜冻日提早显著,20世纪60年代、80年代和21世纪初霜冻日明显推后。2）山西近50a平均初霜冻日的空间分布为“5级阶梯”型分布,从9月中旬到11月上旬,自北向南相继出现初霜冻。3）M—K突变检验表明,山西近50a平均初霜冻日在1989年发生了一次显著的气候突变。4）山西近50a初霜冻日变化趋势的分布具有明显的区域特征,大部分地区初霜冻日有推后的趋势,但南部和中南部的部分地区初霜冻日有提前的趋势。5）Morlet小波分析结果表明,山西的初霜冻、轻微初霜冻、中度初霜冻和重度初霜冻发生的主周期分别为5a、2a、5a、6a;目前山西正处在初霜冻和轻微初霜冻的推后期、中度初霜冻和重度初霜冻的提前期。     

解二维抛物型方程初边值问题的交替方向隐式方法

详细地研究了各种情形的求解二维抛物型方程初边值问题的交替方向隐式差分法。这种方法可以将二维隐式方法归结为求解三对角线性方程组。与一维情况类似,可以继续利用追赶 法求解。它具有运算速度快,存储量小,无条件稳定等优点。该方法是求解二维物型方程的有效方法,必将得到更广泛的应用。     

超定大地边值问题的变分解及相关理论

借助于变分法中求泛函极小化的方法提出求解超定大地边值问题的新方法. 首先依据变分原理引入了超定边值问题的变分解; 其次从理论上证明了变分解的存在唯一性以及最佳逼近性质; 然后讨论了变分解的一般算法, 并在球边界情况下给出了变分解的表达式; 最后利用EGM96重力场模型给出了一个算例, 计算结果表明变分解法能够有效地提高处理重力数据的精度. 总之, 该方法有利于各种类型数据的有效使用, 并且具有严格的数学基础, 能将不同类型的观测数据纳入统一的模型中进行研究.     

基于单频走时敏感度核函数的三维波动方程初至层析

地震波走时层析是主要的近地表速度建模技术.对于变速模型,尤其是复杂近地表模型,波动方程能够更加精确地描述地震波的一阶绕射效应,理论上反演精度高于射线(束)类层析.然而,波动方程的数值求解的计算量较射线理论增加了几个数量级,严重制约了波动方程走时层析技术的三维实用化.为解决波动方程走时层析中遇到的计算和存储问题,本文通过引入随机边界条件以及离散Fourier积分,提出了基于单频梯度反演策略的实用化三维波动方程初至层析方法.本文方法的优势在于：在计算和存储方面,相比于基于波场重构算法(需要三次波传播)的梯度计算策略,本文提出的单频梯度计算方法对内存需求极低(相对于波动方程正演,仅需要增加两个单频波场以及一个成像体),且计算量至少减少1/3(仅需要两次波传播且无需处理数值吸收边界).此外,随机边界的引入使得正演算法大幅简化因而更适用于GPU(Graphics Processing Unit,图形处理器)加速.在理论层面,由于影响地震波一阶散射效应的主要区域集中在连接炮检的中心射线邻域的第一菲涅尔带内,随机边界条件的引入降低了随机边界内产生的散射波的空间相关性,使得第一菲涅尔带内的单频梯度相干加强;同时,地震观测系统的多次覆盖特性有助于第一菲涅尔带外的高波数振荡相互干涉.数值实验表明,光滑处理后的单频梯度仍然保留了非常可观的低波数成分,有助于实现背景速度反演.此外,三维复杂模型反演结果验证了本文方法能够较为准确地反演近地表速度结构.相对于传统的波动方程走时层析实现方案,本文方法在显著降低内存需求的同时进一步提高了计算效率,有望将波动方程初至层析技术推向三维实用化.