基于蒙特卡罗法的三源锥束CT散射分析

在医用X-CT中,散射现象对重建图像的质量有重要的影响。多源锥束CT是医用CT的一种发展趋势,可成倍提高时间分辨率。但与单源CT相比,多源CT额外的光源会导致散射现象的加重,因此在实际应用中需要着重考虑这一问题。本文在开源工具EGSnrc的基础上,开发了一套基于蒙特卡罗法的CT仿真软件SimProj,用该软件获得了贴近...     

Shallow landslides in deeply weathered slates of the Sierra de Collcerola (Catalonian Coastal Range, Spain)

Stability conditions in an area located NW of Barcelona (Spain) are discussed. Here, several mass movements were observed, mainly affecting weathered Paleozoic slates. Many of these failures involved slopes cut along recent infrastructures: debris flows, wedge and plane failures, generally surficial, occurred more frequently. After a detailed geological and geomorphologic survey, geomechanic characterization was carried out, according to RMR and SMR classifications. This rating gave a prediction of slope behaviour, in fairly good agreement with the real observed one.

Stability numerical analysis was carried out for the main cut slopes, based upon the Limit Equilibrium Method. First of all, the deterministic factor of safety was computed using the mean values of parameters. After that, a simulation technique based upon the Monte Carlo Method was applied in order to obtain factor of safety distributions. The probability of failure was estimated as P(F<1).

Finally, results from deterministic and probabilistic approaches were compared. The effectiveness of different possible remedial measures was highlighted by means of a sensitivity analysis, which showed that the more important parameters in the study area are the geometrical ones (height, slope and failure plane angles). The final technical solutions adopted are briefly outlined.     

基于岩体精细化描述的围岩分类及力学参数概率分布特征分析

针对目前水电站地下厂房工程中不同围岩分类方法存在评价结果不一致、围岩力学参数存在室内试验值与实际情况不吻合的现象,现推荐采用岩体精细化描述体系对围岩岩体结构进行定量化评价。将常用围岩分类方法(RMR、Q、RMi、GSI、BQ、HC)评价指标予以归纳分组,并通过各组内不同指标对比分析获得围岩分类方法中的基础评价指标。以大岗山水电站主厂房某区段为分析对象,采取现场岩体精细化地质素描与后期数据挖掘、拟合相结合方法,并依据评价指标间的关联关系,获得了基础、非基础评价指标的分布概型及对应参数,实现对该段围岩岩体精细化描述认知;基于精细化描述结果,应用Monte Carlo法生成符合各评价指标分布概型的大量随机数,而后参照各分类方法评价思路与评分流程,得到评价指标在各分类方法对应的大量随机评分值,通过归纳统计获得不同围岩分类方法评价结果的分布概型;基于各围岩分类方法评价结果与力学参数值之间的关联关系实现对力学参数概率特征分析。该分析方法与思路可为类似工程围岩质量及力学参数的精确确定提供一定借鉴,并可为实现围岩支护极限状态设计提供必要的原始参数支持。     

基于Monte Carlo方法的由GPS观测计算地应变率的误差分析

在利用GPS观测计算地应变率中,应变率的误差分布在地球动力学研究中非常重要. 为解决这一问题,文中提出基于Monte Carlo技术,利用GPS观测的速度计算地应变率进行误差分析的方法;以青藏高原地区为例,利用带有误差的GPS观测数据进行了大量、独立的应变率重复计算,然后运用统计方法给出了应变率分量及其误差的分布. 结果显示青藏高原的喜马拉雅及中东部地区地应变率的东西向、南北向及剪应变率分量的误差相对较大,但整个青藏高原及附近地区应变率结果稳定可靠;研究结果进一步表明,由GPS观测计算应变率的误差主要决定于GPS观测速度的误差, 在变形小的区域（如中国东部地区）,目前GPS的观测精度还不能完全满足地应变率计算的要求.     

双三次数值模式与理想Lamb波、龙卷扰动对比模拟分析

通过"双三次曲面拟合—时间步长积分—双三次曲面拟合……",实现成为一种双三次曲面拟合新算法数值模式。它对模式大气作非线性(空间二阶可导)描述,是用数值分析(三次样条/双三次曲面)与计算方法近似求解大气运动非线性偏微分方程组。用Z坐标系高分辨率双三次数值模式,采用原始大气运动方程组和采用欧拉时间积分方案,设置理想扰动,对大气声波Lamb波做数值模拟;采用大气不可压假定,除去声波,对大气重力波龙卷扰动做数值模拟。对比模拟发现:理想Lamb波以声速在水平方向传播,扰动能量频散快,仅能在原地持续数秒钟;理想龙卷扰动则在原地可持续约1 min,且波动内圈(外圈)自下至上始终处于风场辐合(辐散)状态,具冷心结构和气压驻波特征,可因水汽进入扰动与借助凝结潜热释放而得到发展。     

提取山脊线和山谷线的一种新方法

在研究了现有的仅从山脊线和山谷线的几何特性或物理特性的单一方面设计的提取山脊线和山谷线的算法后,提出了一种基于地形表面流水分析与等高线几何分析相结合的提取山脊线和山谷线的方法。该方法把等高线几何分析的方法与地形表面流水模拟分析的方法有机地结合起来,能够克服各自所具有的弊端。实验结果表明,用本文方法所提取的山脊线和山谷线与实际地形相符合。     

基于局部采样MCMC方法的时移探地雷达反演

马尔科夫链蒙特卡洛方法（MCMC）是一种启发式的全局寻优算法,可以用来解决概率反演的问题.基于MCMC方法的反演不依赖于准确的初始模型,可以引入任意复杂的先验信息,通过对先验概率密度函数的采样来获得大量的后验概率分布样本,在寻找最优解的过程中可以跳出局部最优得到全局最优解.MCMC方法由于计算量巨大,应用难度较高,在地球物理反演中的应用尚处于起步阶段.作者将MCMC方法应用到时移探地雷达（GPR）反演中,并结合双差法对时移的目标区域进行准确反演.在双差法时移反演的过程中,作者对目标区域使用拓展的Metropolis算法进行局部采样,减小了计算量的同时,进一步提高了目标区域的反演精度.最后对合成的数据进行了反演测试,对目标区域的反演误差进行了分析,说明了基于局部采样的MCMC反演方法的有效性,目标区域反演结果特征明显,反演效果好.

     

天文定位中几何精度衰减因子最小值分析

天文定位是一种重要的导航定位方法,被广泛应用于大地天文测量、天文航海等领域。该方法中观测恒星的选择会影响最终的定位精度,目前缺少针对同时测定经纬度天文定位算法中最优选星问题的研究。随着观测仪器自动化水平的提高,观测数据的获取变得更加高效,这就要求研究最优的选星方案以达到最高的定位精度。本文借鉴卫星导航中几何精度衰减因子GDOP的概念,研究了天顶距法中恒星的数量以及分布对定位精度的影响,最后通过仿真试验和实测数据验证得到结论:在天顶距观测误差的统计特性一定时,GDOP能够用来描述恒星的分布对定位结果影响的优劣,且观测的恒星方位角均匀分布时定位误差最小。考虑到不同高度的恒星天顶距大气折射改正残差不同,在实际测量中应尽量采用等天顶距且方位角均匀分布的恒星。     

线性矩法估计参数的保证修正值系数B的推求

在设计洪水的抽样误差估计中,采用保证修正值系数B方法估计均方差简单易用。应用线性矩法估计频率分布曲线的参数,有较好的无偏性和有效性,具有一定的应用前景。目前线性矩法的B值诺模图尚未确定,给用线性矩法确定参数时估计抽样误差带来了难度与不便。采用统计试验,分析了离差系数、偏态系数和样本长度对B值的影响,结果表明线性矩法估计参数时仍可采用保证修正值系数B值诺模图方法来估计抽样误差,并推求了不同偏态系数、设计频率情形下的B值,制成了线性矩法的B值诺模图备用。     

蒙特卡罗法在天然气水合物资源量计算中的应用

天然气水合物甲烷资源量的估算是天然气水合物研究中的热点问题。运用体积法计算含水合物沉积物中的甲烷资源量时,存在参数赋值不确定的问题,从而影响了计算结果的可信度。采用蒙特卡罗法,通过计算样本的频率可以较好地评价和描述计算结果的信度,弥补体积法的不足。对采用蒙特卡罗法估算甲烷资源量的原理进行了分析和探讨,并以麦肯齐三角洲和南海海域水合物为例,计算了水合物赋存区的甲烷资源量。