长时效的正压原始方程能量完全守恒(拟)谱模式

遵循误差反演补偿新计算原理,对正压原始方程传统气象全球拟谱模式方案进行了改造,构造了正压原始方程能量完全守恒全球拟增模式新计算方案,解决了正压原始方程的（非线性）计算稳定性问题和能量守恒整体性质保持问题,改进了相应正压原始方程传统气象全球拟谱模式方案的计算效能。新方案的数值试验表明：在计算实践上,新方案在解决能量守恒问题的同时,可解决（非线性）计算稳定性问题,并在一定条件下可解决非线性计算收敛性问题。进一步的比较数值试验还表明：在计算实践上,新方案具有在提高相应传统气象方案的计算精度,减少其计算量的同时,延长其计算时效,解决其中一类特定“气候漂移”问题方面的效用。本工作原理也适用于斜压原始方程情形。     

甘肃省地震局高性能计算系统

高性能计算是现代科学研究、工程技术开发和大规模数据处理的关键支撑技术,在地震系统发挥着日益重要的作用。甘肃省地震局高性能计算集群系统由60个节点、采用40 GB InfiniBand网络互连组成,峰值计算速度达9.2万亿次/秒。作为地震系统重要的科技服务平台之一,为地震反演、数据分析、数值模拟等领域以及科研创新平台提供高性能计算环境。     

低管电压与低对比剂量在肝脏CT增强中的应用

目的:探讨低管电压(110kVp)与低对比剂量(1.0mL/kg)在肝脏CT增强扫描中的可行性。方法:116例行肝脏CT增强扫描,且BMI<25kg/m2的患者,随机分成A、B两组进行研究,每组58例。A组:管电压110kVp,对比剂碘帕醇(300mgI/mL),对比剂量(1.0mL/kg),采用正弦迭代重建(SAFIRE)。B组:管电压130kVp,对比剂碘帕醇(300mgI/mL),对比剂量(1.2mL/kg),采用滤波反投影重建算法(FBP)。测量动脉期腹主动脉CT值、肝门脉期门静脉CT值及增强扫描3期(动脉期、门脉期、延迟期)的肝右叶肝实质CT值;计算两组图像的肝动脉期对比噪声比(CNR)、信噪比(SNR)及总有效辐射剂量ED;并对两组图像进行评分,进行统计学分析。结果:A组和B组CNR、SNR差异无统计学意义(P>0.05);两组分别在增强扫描3期测得肝右叶肝实质CT值及动脉期腹主动脉CT值差异亦无统计学意义(P>0.05);肝门脉期A组的门静脉CT值低于B组统计学差异有意义(P<0.05)。A组的3期总有效辐射剂量(ED)(9.88±3.56)mSv低于B组(13.91±4.71)mSv,统计学差异有显著意义(P<0.001);A组比B组的对比剂量减少17%。结论:128层MSCT在肝脏增强CT扫描中,应用低管电压(110kVp)和低对比剂量(1.0mL/kg),联合迭代重建技术(SAFIRE),可以得到较好的图像质量,同时降低辐射剂量及人体碘摄入量。      

基于计算岩石物理方法的页岩储层弹性参数提取

发展了应用数值计算方法获取页岩储层的速度、各向异性参数的计算岩石物理系列方法.该系列方法包括了大尺度精细地质模型数值建模、计算网格尺度的地球物理建模和地震波数值模拟提取岩石物理弹性参数.本文方法利用储层的统计数据而不是具体岩心的测量数据,可获得储层岩石物理弹性参数的变化规律.相比于基于岩心测试的岩石物理方法,本文方法可精细考虑实际储层的非均匀特征,可得到岩心测试难以求取的与尺寸效应高度相关的弹性参数,也避免了求取弹性参数变化规律时获取不同地质特征岩心的困难.本文发展了计算岩石物理方法,为计算岩石物理面临的大尺度地质建模和计算能力限制问题提供了有效的解决方案.文中以胜利罗家的页岩储层为例,求得了储层TOC含量从3%到21%变化情况下储层的P波、S波速度以及各向异性参数变化规律.     

地震预测中的地电阻率数据处理方法

本文介绍了地电阻率时间序列数据的常用处理方法,即消除年变化法、无量纲法、相对均方差法和差分法共4类8种方法,包括“九五”以来对原方法的改进和新发展的方法,讨论了各方法的原理、数据处理过程、异常分辨能力、异常指标和异常物理机制及其在数据处理中存在的不足.结果显示:① 一般来说,用有效的数据处理方法才能分析、识别出原始数据曲线上的“弱变化”异常;② 文中8种方法的原理简明,异常物理机制清晰或较清晰,各方法的异常识别指标分别为明确、基本明确和定性的,定性的异常在震情研判中仅有参考意义;③ 消除年变化法和无量纲法通常用于识别地震中期、短临异常,而相对均方差法和差分法通常用于识别短临异常;④ 经数据处理得到的异常与原始曲线的“弱变化”异常相协调;⑤ 时间序列数据出现的异常并不等同与地震孕育、发生过程有直接联系的前兆异常,出现数据异常的台站附近不一定会发生显著地震.      

三维复杂山地多级次群推进迎风混合法多波型走时计算

三维复杂山地条件下的各种地震波型的走时计算技术,可以直接用于复杂山地区域地震波运动学特性的分析、地震数据采集观测系统的设计以及直接基于三维复杂地表的地震数据处理技术的研发.为了在三维复杂地表条件下准确、灵活且稳定地计算各种地震波型的走时,提出一种多级次群推进迎风混合法.该算法利用不等距迎风差分法简洁稳定地处理三维复杂地表及附近的局部走时计算问题,利用计算精度不错的迎风双线性插值法处理绝大部分均匀正方体网格中的局部走时计算问题,利用群推进法模拟三维复杂地表条件下地震波前的扩展问题,利用多级次算法处理各种类型的地震波的走时计算问题.算法分析和计算实例表明:新方法具有很好的计算精度与效率,且能灵活稳定地处理三维复杂地表复杂介质条件下的多波型走时计算问题.     

柴达木盆地东部地震地面运动放大效应

柴达木盆地是青藏高原东北部大型断陷山间盆地,该地区的流动观测记录了2008年11月10日发生于大柴旦附近的M_W6.3地震。和附近的基岩上的记录相比,盆地内部的记录显示出非常显著的地面运动放大效应,表现为峰值速度的增大、持续时间的延长,其呈现出长持续时间的后续震相。傅里叶频谱分析表明盆地内部显著的后续震相的频率和直达波相比较低,地面质点运动轨迹图显示后续震相为面波运动特征。为了解释地面运动的差异,构建二维模型,通过交错网格高阶有限差分方法计算了地震波在盆地内部的传播过程,结果显示盆地内部低速层的存在造成直达波的放大以及多次反射与转换,盆地边缘结构造成的波的相干叠加产生了强烈的次生面波,其低频、大振幅、长持续时间的特征是盆地内部地面运动放大的主要原因。     

“房桥合一”铁路客站抗震计算的快速高效多维虚拟激励法

针对"房桥合一"铁路客站节点数多、空间相互作用显著等特点,利用沿结构整体坐标系方向功率谱的表达式和合理谱强度因子的确定方法,提出能合理考虑地震激励输入角度变化的快速多维虚拟激励法,并以天津西站Ⅱ区作为工程案例,对利用快速多维虚拟激励法、时程响应分析法和已有虚拟激励法的计算结果进行比较。结果表明:三种算法的计算结果吻合较好,证明快速多维虚拟激励法的计算值是合理的;且与已有虚拟激励法相比、快速多维虚拟激励法的计算效率显著提高,更适合于"房桥合一"铁路客站的抗震计算。     

一种锥束CT动态Bowtie的蒙特卡罗模拟研究

在CT中,前置于患者的衰减器（"Bowtie滤波器"或"Bowtie"）被用来减少患者所受的辐射剂量和帮助满足探测器的动态范围,同时Bowtie也能提高CT图像质量。传统的静态滤波器对调制X射线能力有限,我们之前提出了一种能够有效调制X射线的动态Bowtie滤波器。本文使用蒙特卡罗模拟方法设计和研究提议动态Bowtie滤波器。首先,建立了一个仿真框架,包括一个目标试样和一个由反射靶生成的多能谱射线源。接着,一个设计试样被用来辅助设计针对于目标试样的动态Bowtie滤波器。最后,我们分别从没有Bowtie滤波器和带有动态Bowtie滤波器两种模式中,模拟获得两个试样的投影数据并进行分析研究。结果显示,在设计的动态Bowtie作用下,到达探测器的X射线动态范围明显减小、X射线光子通量一致性很好,并且动态Bowtie可以减小辐射剂量和提高CT图像质量。      

基于多相态KE-I关系的降水动能计算方法

降水动能描述了水凝物到达地表的动能大小,是计算土壤颗粒分离和土壤侵蚀等地表物理化学过程的重要参数之一。针对目前降水动能计算多考虑降雨特征,忽略其他相态降水的问题,本研究基于2018年全国25个雨滴谱仪站点的降水观测数据,结合速度(V)-直径(D)关系对水凝物进行相态划分,拟合了不同相态降水的降水动能-降水强度(KE_V-I)经验公式,总结出适用于雨、雪、雹的动能计算关系模型,并对各相态降水动能进行偏差分析和校正。结果表明：(1)降雨动能计算适用KE_V-I指数型关系公式,降雪和降雹动能计算适用幂函数型关系公式;(2)降水动能偏差与总降水中非液态占比呈正相关,仅以降雨KE_V-I经验公式估算所有相态降水的动能会显著高估降雪动能和低估降雹动能,对降水整体估算偏差在-0.05到0.31之间;(3)通过建立各相态的KE_V-I计算模型,降水整体估算误差得到有效降低,降雪偏差从7.08降至-0.05,降雹偏差从-0.29降至0.03。本文所提出的方法能够有效提高降水动能估算的准确度,为区域土壤水力侵蚀计算提供更好的支...