基于空间填充曲线排列码划分的并行缓冲区算法

缓冲区计算是地理信息系统空间分析的基本功能之一。在矢量缓冲区研究方面,以往大多聚焦缓冲区生成算法,这些优化方法在单机计算环境中针对大规模空间数量数据的计算效率提升是有限的,因而需要并行环境下的改进方法来进一步破解计算瓶颈。本文基于Map Reduce计算模型和分布式内存计算框架,提出了一种基于空间填充曲线排列码划分的并行缓冲区分析算法(SPBM),基于空间填充曲线编码进行数据排序和范围划分,对跨越网格的数据采用近似切分方法,在考虑负载均衡的前提下对任务并行分解,之后按照给定深度的"树状"方式进行结果合并。本文在单机和集群两种环境下利用实际道路网等数据进行了实验。同等环境下,相较于目前流行的GIS软件——QGIS和Post GIS计算性能提升明显,相较于现有其他并行优化方法效率提升超过50%。这样的优化分析方法对于GIS中其他矢量分析算法的并行算法也具有一定的借鉴意义。     

三层模式中风场对背风波的影响

为了研究风场对背风波的影响,针对边界层附近为弱稳定层结的背风波,建立了一个三维三层的理论模型和线性计算模式,分析了各层中风速和风向的变化对背风波特征的影响,揭示了气流过孤立山脉产生背风波的有利风场条件。结果表明：背风波的波长、振幅等特征对各层风速和风向的变化具有相当的敏感性,波长随着低、高层风速的增大而增大,随着中层风速的增大先减小后增大;振幅随着低、中层风速的增大先增大后减小,随着高层风速的增大而增大。此外,风速和上下层风向切变的增大均使背风波的形态逐渐由横波型转为辐散型,但是上下层风向的切变对背风波形态的影响比风速更为显著。     

含气介质临界点、流体和骨架弹性参数的数值计算方法

目前临界孔隙度及孔隙介质的研究主要为测试分析的实验方法,如何运用数值计算方法求取孔隙流体介质临界点和组分弹性参数一直是备受关注的课题.本文提出了求取孔隙介质临界点、孔隙内流体和骨架弹性参数的数值计算公式及方法,并结合含气样品测试数据实现了这种数值计算.文章首先给出了3个包含孔隙度的线性方程,通过有机组合每个线性方程中的两个系数得到求取有关弹性参数的数值计算公式.然后,详细阐述了计算的方法和步骤以及需要注意的问题,把室内含气砂岩样品实测的综合介质密度、纵波速度和横波速度作为数值计算的输入数据,求得了临界点和流体及骨架等相关弹性参数的具体数值.通过比较分析本文数学方法求得的计算数据与实验方法测得的实测数据,表明了该数值计算公式的正确性和数值实现方法的有效性.     

基于交错网格Fourier伪谱微分矩阵算子的地震波场模拟GPU加速方案

作为高精度波形反演或逆时偏移的重要组成部分,地震波数值模拟对计算速度和效率提出了更高要求.GPU通用计算技术的产生及其内在数据并行性,为高效地震波数值模拟应用和研究得以有效开展奠定了基础.本文借助交错网格的Fourier伪谱微分矩阵算子和GPU上高效矩阵乘法,实现了复杂介质中地震波模拟的高效算法.数值试验表明,优化后的GPU计算相比CPU单核计算在大规模二维地震波场计算中获得至少100x以上的加速比.这对我们快速分析目标反射层在地震剖面中同相轴位置,制定优化采集方案具有重要意义.     

金融危机下国家大规模投资政策模拟

构建了一个中国多区域可计算一般均衡模型(MReCGE-C),利用该模型分析了我国应对国际金融危机的大规模投资政策。MReCGE-C根据城乡属性和年龄结构对居民进行了分组,放松了要素流动的限制,使劳动力和资本可以同时跨区域跨部门流动,并设计了一种基于居民终生效用的区域均衡机制,通过区域差距变量可以调控和测度区域差距。最后利用MReCGE-C模型针对国家大规模投资计划下的3种典型投资方案进行了模拟。结果发现在等量新增投资下,针对中部地区投资方案的GDP增长率高于针对中西部和全国投资方案的GDP增长率,同时针对中部地区投资方案对于城乡就业的拉动也要优于其余两种投资方案,而且可以实现缩小区域差距的目标。综合比较,目前国家大规模投资的优势应该指向中部地区。     

浅埋隧道围岩应力场的计算复变函数求解法

对于受地表边界和地面荷载影响的浅埋隧道的围岩应力场,由于在数学处理上存在一定的困难,很难用解析解来进行分析,而通常采用边界元或有限元的数值方法来解答。为了求解浅埋隧道的应力场,采用边界配点来确定边界条件,同时用保角映射将一个含圆孔的半无限空间区域映射为圆环域,然后把这个区域内的解析函数展开成Laurent级数的形式,利用Muskhelishvili的复变函数理论和最小二乘法来确定解析函数的各项系数,从而求得浅埋隧道围岩压力的半数值、半解析解,最后通过算例给出了围岩应力的分布情况。计算结果表明,该方法计算精度高、计算量小,具有应用价值。     

不同频率循环荷载作用下岩石阻尼比试验与变形破坏机制细观分析

通过对取自海南昌江核电厂一期工程主厂区的黑云母花岗岩,在RMT－150B多功能全自动刚性岩石伺服试验机进行单轴不同频率的循环加卸载试验,研究花岗岩动应变-应力滞回曲线、动弹性模量和阻尼比同循环频率之间的规律,然后对破坏断裂断口进行细观电镜扫描试验,分析了不同频率下花岗岩断口细观形貌与结构特征。结果表明,花岗岩在循环荷载作用下的加卸应力-应变曲线并不重合,而是形成一个封闭的滞回环;当循环频率在0.01～1.0 Hz范围内变化时,随着频率增大,花岗岩滞回环面积、动弹模和阻尼比都增大;通过细观数字图像试验分析发现,循环荷载作用下,随着频率的增大,岩石断口细观裂纹总的数目、面积、分形维数均变大。文中给出不同频率循环荷载作用下岩石破裂的细观损伤力学机制,研究成果对海南昌江核电厂的地震反应分析和场地安全性评价有着重要参考价值。     

非均布列车荷载作用下软土路基的振动分析

基于实测资料,假设列车轮载经钢轨传递下来的荷载由5根枕木承担。根据分担比,获得了道渣层表面非均布荷载的表达式。将道渣层视为单相弹性介质,软土地基被看成是考虑水土耦合作用的饱和多孔介质。借助势函数,利用Helmholtz矢量分解法及Fourier变换技术分别对弹性土体和饱和半空间土体进行求解,得到非均布移动荷载作用下软土地基位移、应力及孔压响应在变换域内的精确解答。利用FFT算法得到了数值结果,详细分析了荷载分布形式、观察点位置、道渣性质以及软土地基渗透系数对动力响应的影响。研究结果表明：高速情况下动力响应与低速情况有很大不同;对于软土路基,应特别注意列车高速运行时路基浅层范围内产生的孔压响应。     

循环荷载下双曲线模型修正土体一维固结解

假定土体的应力-应变关系满足双曲线模型、渗透系数的降低与压缩系数的减小成正比以及初始有效应力沿深度均匀分布,推导出任意荷载下单层地基的非线性一维固结方程。结合具体的定解条件采用分离变量法求出土体有效应力的通解,并进一步求出超静孔隙水压力、按沉降定义的平均固结度Us以及按孔隙水压力定义的平均固结度Up的通解。根据通解求出梯形低频循环荷载作用下的解,并找出循环荷载作用下土体固结性状的影响因素。基于循环荷载作用下的解,采用FORTRAN语言编制了计算程序,计算并分析考虑土体本构关系非线性时的固结性状,以及反映土体本构特征的参数E和n对固结性状的影响。计算结果表明：①Us稍大于Up;②距排水平面越近,附加有效应力受加载方式的影响越明显,距排水平面越远则相反;③E和n对Us的影响很小,可以忽略不计;④E和n对Up和附加有效应力的影响较大且不容忽略。     

海面-地面系统中三维古地面反演计算的剥蚀算法研究

三维古地面反演计算是海面-地面系统研究中的一个基础研究。该文在前人研究的基础上提出了三维古地面反演计算的剥蚀算法,并与传统的加积算法做了比较分析,得出剥蚀算法比加积算法有对时间序列离散判断简单、方便计算等优点。同时,在GIS平台下实现了该算法,并以长江三角洲地区为试验区,进行了古地面虚拟反演计算,取得了较好的效果。