大气污染对太行山中部地区地表风的影响

为了研究大气污染对太行山中部地区地表风的影响,我们对阳泉、榆社(高山站),石家庄、邢台(平原站)4个站点1966~2005年间的能见度、近地面温度、近地面风速数据进行了统计计算与趋势分析。结果显示:在平原站能见度相对山坡站下降更加明显的背景下,平原站的近地面温度、近地面风速、850hPa风速都呈下降趋势;而山坡站的近地面温度、近地面风速呈上升趋势。分析表明:(1)由于气溶胶的辐射效应与冷却效应,抑制了垂直通量的上下传输,致使平原站下午的近地面气温呈下降趋势,平原站和高山站的地表风速呈相反的变化趋势。(2)平原站850hPa (与高山站高度相近)风速呈现下降趋势,印证了高山站的近地面风速增加是气溶胶的辐射效应减弱了垂直能量交换造成的。      

西南天山隆起时代的河床砂岩屑磷灰石裂变径迹证据

现代的天山山脉是在古生代造山基础上,于新生代强烈抬升而形成.其新生代造山和隆升过程,造就了现今的天山地貌格局.本文选取西南天山作为研究区域,采用河床砂岩屑裂变径迹测年分析,从统计角度限定西南天山的隆升-剥露过程.样品采集于特克斯河支流阿克雅孜河、夏特河、木扎河以及特克斯河干流的沉积河床.磷灰石裂变径迹测试和统计分析表明,存在代表源区热史演化不同阶段的年龄峰值.尽管不同样品的年龄众数分布有少许差别,颗粒年龄众数的去褶积分析获得了西南天山山体新生代冷却的三个基本一致的阶段：6~8 Ma,12~19 Ma以及32~40 Ma.结合山脉隆起的地质地貌模型,无论是整体抬升或掀斜抬升,以及压扭性背景的花状挤出抬升,根据磷灰石裂变径迹封闭温度推断的抬升量与现今天山高度基本相当的事实,都可以确认西南天山山体是6~8 Ma以来形成的.天山这三期快速抬升冷却事件与青藏高原及其周边的主要隆升时期有较好的对应,证明了天山隆升和印度-欧亚板块碰撞远程效应的关系.另外,6~8 Ma的冷却事件与沉积地层学研究揭示的6 Ma左右的气候显著变化相互印证,显示了研究区域山脉隆升和气候变化之间存在的密切关系.     

基于非牛顿指数描述的非达西渗流一维固结分析

基于非牛顿指数描述的非达西渗流定律,同时考虑地基内部竖向附加应力随深度线性变化以及变荷载的影响,建立了一维固结控制方程并应用有限差分法进行数值求解,同时对不同参数单级加荷下的固结性状进行分析。结果表明：基于非达西渗流比达西渗流下固结速率要慢,且渗流模型中非牛顿指数越大,土层的固结速率越慢;土层厚度越厚,固结速率越慢,因此,传统固结理论中室内土样与地基土层之间的相似关系不再成立;作用于土层的平均附加应力越大,土层的固结速率越快;在单面排水情况下,附加应力分布对土层固结速率有较大影响;相反,双面排水条件下土层固结速率与附加应力的分布是无关的;荷载的加荷速率越快,则土层的固结速率越快。最后,讨论了达西渗流计算固结变形的适用范围。     

基于逐减随机震源采样法的频率域二维黏滞声波方程全波形反演

全波形反演方法利用叠前地震波场的运动学和动力学信息重建地下速度结构,具有揭示复杂地质背景下构造与岩性细节信息的潜力。然而,巨大的计算量是阻碍其发展的一个瓶颈问题。为此,研究者们提出了震源编码技术来减少计算量,但是此方法在模型更新过程中会引进随机串扰噪声,降低反演结果准确性。所以,在保证计算精度的情况下,本文提出了采用逐减随机震源采样的方法来高效计算全波形反演问题。笔者将此方法应用于频率域二维黏滞声波波动方程全波形反演,开始了在频率域进行随机震源采样类方法的研究,计算过程中共使用了依次增大的8个频率段;并应用Overthrust模型来验证此类随机震源采样法的正确性。实验结果表明：基于逐减随机震源采样法的反演结果与实际Overthrust模型的拟合误差为0.065 65,而应用基于全部震源的全波形反演方法得到的反演结果与实际Overthrust模型的拟合误差为0.064 64,两者差别不大;但计算用时由740 min减少到291.2 min,即计算效率提高了2.54倍。为了更好地确定方法的有效性,将其应用于Marmousi模型进行试算。模型试算结果表明：基于逐减随机震源和基于全部震源得到的反演结果与实际Marmousi模型的拟合误差分别为0.080 12和0.078 97,相差不大;但计算用时由1 218.9 min减少到274.4 min,计算效率提高了4.44倍。综上,在保证反演精度的情况下,基于逐减随机震源采样法的频率域全波形反演方法大大减少了计算量,具有不可替代的计算优势,并且没有引进随机串扰噪声。     

椭球面三角形外接大地圆圆心测地坐标的另一种解法

为了在测地坐标系中确定外接大地圆圆心,采用了解非线性方程组的牛顿迭代法,解法简便,易于理解.数据验证表明,由该数值法求解的结果与由解析公式获得的结果完全一致.由此亦表明,一些繁难的大地测量计算也可改用数值法来进行.     

GPS测码伪距绝对定位的几种算法

GPS定位方程是非线性的,一般处理方法是按泰勒级数展开取至一次项进行线性化,再利用最小二乘原理求解,如果所取观测站坐标的初始值具有较大的偏差,略去二次微小量的模型误差,对解算结果将产生不能忽略的影响。本文研究分析了GPS测码伪距绝对定位的传统算法,并提出了一种通过求差法将GPS绝对定位的非线性观测方程转化成线性方程直接求解测站坐标的新算法,通过实例计算表明该方法计算简单,不需要测站的初始坐标信息,不需要求导计算和迭代计算,对于提高GPS测码伪距单点定位的解算速度和精度具有重要的意义。     

基于插值理论的GRACE卫星精密轨道内插的研究

利用GRACE卫星高采样率精密轨道,进行了拉格朗日插值、牛顿插值和线性插值方法的研究,并对插值结果进行了比较和分析,探讨了插值方法的插值精度及其应用于低轨卫星轨道内插的适用性。     

一种求解贝叶斯模型平均的新方法

贝叶斯模型平均（Bayesian model averaging,BMA）是最近提出的一种用于多模式集合预报的统计方法.进行贝叶斯模型平均需要准确估算模型集合中每个竞争模型的权重与方差,经常采用的方法是期望最大化（Expectation-Maximization,EM）方法与马尔可夫链蒙特卡洛（Markov Chain Monte Carlo,MCMC）方法,两种方法各有优劣.本文首先对BMA的（对数）似然函数进行改进使之无需BMA权重之和为1的显式约束,并利用一种有限记忆的拟牛顿优化算法（LBFGS-B）对其进行极大化,由此提出了一种求解贝叶斯模型平均的新方法（BMA-BFGS）.采用三个陆面模式进行的土壤湿度多模式数值模拟试验表明：在计算精度方面,BMA-BFGS的精度与MCMC方法几乎一致,优于EM算法;在计算耗时性方面,BMA-BFGS的计算耗时与EM算法相当,远小于MCMC方法.     

柴木铁路沼泽化冻土区热管冷却半径的观测研究

柴达尔-木里铁路沿线高温沼泽化多年冻土广泛发育,约20km的冻土路基采用了热管降温措施．借鉴青藏铁路的成功经验,热管的纵向和横向间距均设计为3m,并针对铁路东南-西北的不利走向,在阳坡布设了双排热管．为验证热管措施在多年冻土湿地路段设计参数的合理性,2008年底在柴达尔-木里铁路中段布设了热管冷却半径观测试验场,该观测...     

日本福岛核危机回顾

日本福岛核电站位于日本福岛工业区双叶郡,包括福岛第一核电站(6台机组)、福岛第二核电站(4台机组),均使用混合氧化物燃料。福岛第一核电站1号机组于1971年3月投入商业运行,第二核电站1号机组于1982年4月投入商业运行。福岛核电站的核反应堆都是单循环沸水堆,只有一条冷却回路,蒸汽直接从堆芯中产生,推动汽轮机。福岛核电站一号机组已经服役了40年,此前出现了许多老化的迹象。这一机组计划延寿20年,正式退役要到2031年。