基于Excel VBA实现箱图剔除法确定化探异常下限

尽管目前地球化学异常下限的计算方法较多,但不同程度上存在计算过程复杂、应用范围较窄、相关专业软件普及程度低等问题。针对这些问题,笔者提出将EDA技术与经典统计概率法相结合的方法,并应用Excel下的程序开发工具进行VBA编程,在Excel中实现异常下限的快速计算。该方法具有计算简便、异常下限稳健以及应用范围广等优点。依据方法原理,该方法仅适用于来自单一母体的数据样本,来源母体可通过分形图给出初步判定方法。     

饱和多孔介质三维时域黏弹性人工边界与动力反应分析的显式有限元法

本文基于Biot的饱和多孔介质本构方程,考察具有辐射阻尼的外行球面波,推导了饱和多孔介质三维黏弹性人工边界的法向和切向边界方程;在已有的饱和多孔介质二维显式有限元数值计算方法基础上,提出该理论的三维方法,并开发了实现该三维方法的有限元程序.算例表明饱和多孔介质三维时域黏弹性人工边界与动力反应分析的显式有限元法具有较好的精度和稳定性.     

基于非分裂CFS-PML边界条件的频散介质GPR时域有限元模拟

通过数值模拟研究探地雷达（GPR）高频电磁波在频散介质中的传播规律,对提高实测资料的解释精度具有重要意义.复频移完全匹配层边界条件（CFS-PML）以其优越的吸收特性被广泛用于一阶电磁波动方程的GPR时域有限差分数值模拟中,其实现过程大都涉及电磁场的卷积计算,辅助变量较多,降低计算效率.为此,本文从复拉伸坐标系下的Debye频散介质电磁波动方程出发,通过合理构造辅助微分方程,推导了二阶Debye频散介质电磁波动方程的非分裂CFS-PML边界条件实现公式,避免了电磁波场的分裂和卷积计算.在此基础上,利用Galerkin法和Newmark-β差分法推导了基于非分裂CFS-PML边界条件的GPR有限元方程及其时域差分离散格式.两个GPR模型的模拟结果表明：本文提出的基于辅助微分方程的非分裂CFS-PML边界条件实现方法可有效地吸收大角度入射的低频虚假反射波,提高模拟精度;相比于非频散介质,高频电磁波在频散介质中传播衰减更强、子波持续时间增大、分辨率和传播速度降低、直达波和反射波的主频更小,分析结果有助于提高实测GPR资料的解译精度.

     

基于紧凑存储的任意四边形网格有限元法地震波场数值模拟

采用有限元法精确模拟复杂介质条件下的地震波场,一般采用三角单元,但在节点数相同的情况下,三角单元的计算精度不如矩形单元高,采用三角单元模拟复杂界面时,编制确定结构刚度矩阵非零元素位置的程序也较为麻烦。采用矩形单元离散含有倾斜或起伏界面的地质模型时,无法避免绕射噪声,若加密网格又会增加计算量。为此,本文基于任意四边形单元模拟声波的传播,在倾斜或起伏界面条件下,可以有效避免因离散引起的“阶梯状”界面,在不增加计算量以及内存占用的前提下,有效地消除离散绕射噪声。采用对角的集中质量矩阵代替一致质量矩阵,避免矩阵的求逆运算,从而提高显式有限元法的计算效率;对结构刚度矩阵采用紧凑存储格式,每一行需要存储的元素最多为5个,同时零元素不参与运算,既减少内存的占用,又极大地提高计算效率。     

1911年红海湾6级地震前珠江三角洲及其邻区的区域构造应力场

本文包括：（１）通过拟合区内断层错动方式确定区域边界上作用的远场构造应力;（２）根据发震前震中所在断层段要进入塑性状态,确定１９１１年红海湾地震前区内构造应力场．两个应力场均采用正交设计法与弹塑性有限元方法快捷地得到．结果表明：该区远场最大主压应力σ３的方向为Ｎ６８°Ｗ,σ３＝－２００ＭＰａ,最小主压应力σ１＝－８０ＭＰａ;算得的红海湾地震前区内应力场的特点是除滨海断裂的红海湾段已进入屈服状态外,其余断裂均处于弹性状态．     

北京平原第四纪下限问题研究

本文首先回顾了北京平原区第四纪下限的研究历史,指出迄今为止在这一领域存在的问题。根据我国目前的研究进展和北京地区研究现状,按照全国地层委员会的相关标准和要求,确定了划分第四纪下限的5大原则。北京平原区第四系沉积特征复杂,根据多年资料和研究结果,我们提出了划分该界限的方法。广大的北京平原区,不具备完整的第四系剖面露头,生物地层方法划分存在困难;古地麟、热释光及孢粉地质测试为准确定位这一界限提供了可能。对平原区的钻孔而言,笔者发现以古地麟测试结合孢粉分析、岩性特征、地球物理测井以及地球化学测试（碳酸钙分析）等综合手段确定这一界限的证据更有说服力,新5孔应用效果较好。     

EFG法在土体固结中的应用

EFGM作为一种新的数值计算方法，具有只需节点信息而无段单元的特性，故在解固结方程方面有很大的优势。在计算分析中，此法容易构造固结方程的EFGM刚度矩阵和处理不同边界条件。对单面排水等条件的计算结果表明，EFGM在解决固结变形问题上，精度较高，处理边界准确。     

小型主动型氢原子钟的核心——蓝宝石谐振腔的仿真研究

该文对小型主动型氢原子钟的核心——蓝宝石谐振腔进行了研究,利用采用有限元数值分析方法的软件HFSS V10(High Frequency Simulation Software)对日本NICT(National Institute of Informa-tionand Communication Technology)和上海天文台研制的两种不同的蓝宝石谐振腔进行仿真分析。对HFSS软件中TE011振动模式的识别进行了讨论,并得出了识别TE011模的简便方法。对上述两种腔的腔频随外金属腔及填充蓝宝石几何尺寸和温度的变化进行了计算。比较了两种腔的有载Q值,对判断蓝宝石谐振腔是否符合脉泽自激振荡条件的腔的填充因子η′和S参数进行了分析计算。     

滤管两端均不在含水层层面的承压不完整井近似计算方法研究

不完整井在实际工程中应用广泛,且存在某些承压不完整井,其滤管顶、底均不在含水层层面,不能直接按常规的方法进行计算。可根据流线图中的水平流线对渗流场分层后进行简化计算,但分层位置的确定非常困难。因此,提出了简便可行的渗流场分层方法,包括假设上下部渗流场渗径相等的等渗径改进方法、假设上下部渗流场中井滤管长度与其含水层厚度比值相等的等比例改进方法和以等比例方法为基础的等效阻力法,并以此为基础对计算简便、精度较高的沙金煊公式进行完善,得到了适用于滤管顶、底均不在含水层层面的承压不完整井近似计算方法。经过算例验证和与有限元比较,该方法计算精度较高,实用性较强。     

软土固结中孔压消散的有限元方法

在有凶中实现土体固结时孔压消散有两种方法,并对一实例作了分析,比较了它们的异同。