有限区域细网格模式产品在江汉平原暴雨预报中的应用

利用湖北省提供的有限区域细网格模式数值预报产品及降水实况资料，用逐步回归法拟合出江汉平原等区域暴雨方程。方程的拟合效果比较理想，此方程尚有待台站业务考核及改进。     

岩石有限应变测量的莫尔圆图解网法

通常进行的应变测量需经过繁杂的数学计算或莫尔圆作图才能获得应变结果。本文介绍了笔者新设计的莫尔圆图解网及其使用方法。它可以快速准确地作出莫尔圆,进而求得应变结果,也可用于求解应力莫尔圆问题。     

琼西戈枕断裂带的有限应变分析及其意义

琼西戈枕断裂带是海南岛最重要的金矿带,断裂带构造岩及金矿体中广泛发育变形石英颗粒.这些变形石英颗粒的三维有限应变分析表明:糜棱岩型金矿床形成于总体伸展构造背景,与断裂带印支期的韧性拆离构造有关;破碎带蚀变岩型金矿床形成于总体挤压构造背景,与断裂带中印支晚-燕山早期的韧-脆性逆冲推覆构造有关;石英脉型金矿床形成于总体平移剪切状态,与断裂带中燕山中晚期的脆性平移剪切运动、块断滑脱运动有关.     

有限断层破裂时空过程反演结果误差的随机加权估计

针对非线性方法对有限断层破裂时空过程反演中存在的不唯一性问题,应用随机加权方法对反演结果进行误差估计．通过数值模型计算,定量揭示了滑动量、滑动角、破裂时间及上升时间在断层面上的分辨率值分布．应用随机加权方法对台湾海峡地震及缅中边境地震的误差估计结果表明,对应两次地震破裂中心位置上的参数分辨率值最高,结果最可靠,而在边缘子断层处滑动量与上升时间的分辨率值较低．      

基于BISQ模型双相各向同性介质中地震波数值模拟

基于BISQ机制,推导了双相各向同性介质中弹性波数值模拟方程,并借助交错网格有限差分方法进行了波场模拟。证实了三种波(快纵波、慢纵波、横波SV)的存在和波场特征;在合成地震记录中,观测列三种波的反射,由于各种波相互之间的转换,转换波在地震记录中显示也很清楚,这使得波场变得更为复杂化。     

孔隙水压力静力触探动态贯入过程的有限元模拟

孔隙水压力静力触探作为一种新型的土体勘测技术已被广泛应用,但是,其机理的研究尚存在很大的不足。笔者系统分析了孔压触探的贯入特征,提出了“孔隙水压力体积基准值”的概念,借此推导出用土体体积变形量计算超孔隙水压力的公式。在此基础上,用轴对称有限单元法模拟探头的动态贯入过程,揭示应力－应变的变化特征,寻找超孔隙水压力生成与消散的规律,对孔压触探的机理进行了研究。     

非均匀介质中交错网格高阶有限差分数值模拟

地震波场的数值模拟一直是地球物理学的一个重要的研究领域,而在数值正演模拟方法的研究中,计算精度和计算效率是评价该方法有效性及优越性的二个关键问题。这里从一阶速度—应力弹性波动方程出发,着重介绍如何构造离散化模型的网格,如何求解空间导数,如何选取边界条件等内容,从而更有效地提高数值计算的精度与计算效率。文中构造了不同类型的介质模型,并在交错网格中,利用高阶有限差分模拟非均匀介质的波场传播。模拟结果表明,该方法实现简单,具有很好地稳定性和较高的精度,能够直观、高效地反映出介质中波场的传播规律。     

井地直流电法三维数值模拟中若干问题研究

讨论了地下垂直线源分段计算和场叠加的方法,并实现了在套管上供直流电的三维数值模拟。讨论了大型容量矩阵的压缩存储方式,采用数组和结构体相结合的方法实现容量矩阵的一维链表式压缩存储。在求解超大型稀疏线形方程组时引入不完全Cholesky分解稳定化的双共轭梯度算法(ICBG),通过与均匀半空间垂直线源解析解的对比,证明了该算法是准确可靠的。     

洞庭湖泥沙淤积数值模拟模式

在确保沙量守恒的非恒定流、非均匀沙的二维水沙数学模型的基础上,通过洞庭湖水沙输运和河床变形计算实践,比较了两种冲泻质与床沙质的转化模式,提出了包括冲泻质在内的泥沙淤积模式和水沙动边界计算模式.从洞庭湖泥沙输运和河床变形计算的结果看,洞庭湖泥沙淤积数值模式具有较好的稳定性和较高的模拟精度,说明所建湖泊泥沙模型的良好性能和具有推广运用的前景.     

改进的传染型余震序列模型

将有限幂律公式（LPL）应用到传染型余震序列模型（ETAS）中，并对原ETAS模型加以改进。以台湾集集地震早期余震序列为例，对原ETAS模型与改进的ETAS模型进行了对比分析，结果表明改进的ETAS模型要优于原ETAS模型。