利用PETSc和FEPG编制海啸数值模拟程序的研究

海啸的数值模拟是海啸研究的一个重要领域, 它对于帮助理解海啸的基本物理特性和预防减灾具有重要意义。 海啸数值模拟程序的编制是一项繁杂的工作, 该文介绍了利用两种功能强大的通用软件(PETSc和FEPG)来进行海啸数值模拟程序编写的研究。 PETSc和FEPG采用有限差分、 有限元、 有限体积等多种离散方法, 可以对数值问题给出稳定的求解。 该文采用了直角坐标系下的非线性浅水波方程作为海啸波传播的控制方程进行离散求解, 并将其结果与TUNAMI N1模型进行了比较, 表明这两种方法方便而且有效。     

地震波有限差分数值模拟框架下的起伏地表处理方法综述

本文查阅了近30年来国内外关于有限差分方法模拟近地表复杂介质方面的文献,分析了在有限差分法框架下复杂地表地形起伏问题的处理及实现方法.总结了地震波有限差分数值模拟方面的成果.复杂地表类型很多,对地震资料的影响程度也不尽相同,要想完全描述近地表影响因素是不可能的,而关键的问题就是自由边界条件方程的处理.所有处理及实现方法中不外两种:直接法(局部)和间接(全局)法.文中列出的方法,各有特点.我们的目的是寻找合适的方法解决适合的问题.     

考虑动水压力影响的单柱式桥墩地震反应分析

在Morison方程的基础上,用附加水质量法考虑动水压力对桥墩的影响,以单柱式桥墩为研究对象,以ABAQUS有限元软件为计算平台,建立了考虑桩-土动力相互作用的单柱式桥墩地震反应分析模型,考虑土体和桥墩混凝土的动力非线性特征,分析了地震动作用下动水压力对单柱式桥墩的墩顶相对墩底位移、加速度、剪力和弯矩反应的影响,并探讨了水位对单柱式桥墩地震反应特性的影响。结果表明:动水压力改变桥墩的地震反应特性,增大了桥墩顶部相对底部的位移、墩顶绝对加速度和墩底的内力,水位变化影响桥墩的地震反应特性。对于深水桥墩抗震设计计算,考虑动水压力效应、水位变化是有必要的。     

浅海内波垂直模式方程的一种解析解

给出了一些海区实测浮力频率剖面的数学模型.利用超几何函数对内波垂直模式方程进行了解析求解,得到内波模式方程的解析解和频散方程.利用得到的解析解给出了一种理想条件下内波前三阶简正模式的结构.     

A time-dependent numerical model of the mild-slope equation

On the basis of the previous studies, the simplest hyperbolic mild-slope equation has been gained and the linear time-dependent numerical model for the water wave propagation has been established combined with different boundary conditions. Through computing the effective surface displacement and transforming into the real transient wave motion, related wave factors will be calculated. Compared with Lin’s model, analysis shows that calculation stability of the present model is enhanced efficiently, because the truncation errors of this model are only contributed by the dissipation terms, but those of Lin’s model are induced by the convection terms, dissipation terms and source terms. The tests show that the present model succeeds the merit in Lin’s model and the computational program is simpler, the computational time is shorter, and the computational stability is enhanced efficiently. The present model has the capability of simulating transient wave motion by correctly predicting at the speed of wave propagation, which is important for the real-time forecast of the arrival time of surface waves generated in the deep sea. The model is validated against analytical solution for wave diffraction and experimental data for combined wave refraction and diffraction over a submerged elliptic shoal on a slope. Good agreements are obtained. The model can be applied to the theory research an d engineering applications about the wave propagation in a biggish area.     

1975年8月河南省特大暴雨雨强极值的重新估算

通过对云中气载液态水含量以及水汽和气载液态水之间的平衡方程的讨论,建立了改进的可降水量Ia公式。按照改进的公式,重新估算了1975年8月河南省特大暴雨7日20时雨强。结果表明,分别以遂平气象站记录的8 m/s和按物象估算的12 m/s地面风速为基础的改进的板桥水库附近8月7日20时雨强估算值Iae1=75.0 mm/h和Iae1=95.8 mm/h很接近实际的雨强极值Ioe=99.7 mm/h。     

海口永庄气候平均差值的订正分析

根据海口、永庄两站2006年1~12月对比观测气象资料进行分析,得出海拔高度和周围环境不同是形成海口与永庄两地气候差异的主要原因。利用两地气候平均差值与有关气象要素存在很好相关关系,应用回归法建立气候平均差值的订正方程。经检验,方程回归效果较好,因此,可利用气候平均差值订正方程将永庄站的气象资料订正到搬迁前等长序列,并求取各气候要素的各个时段历史平均值。     

L和C波段雷达干涉数据矿区地面沉降监测能力分析

雷达成像的波长、入射角、地面分辨率等参数严重影响着SAR差分干涉测量地面沉降的监测能力和精度,论文通过理论推导和矿区实际沉降差分干涉相位模拟,从监测到的最大沉降梯度和沉降量、保相能力、对微小沉降的敏感程度等方面对L和 C波段雷达干涉数据的矿区地面沉降监测能力进行分析;精化双轨D-InSAR数据处理的流程、方法和相应参数,使用ALOS PALSAR和ENVISAT ASAR数据获取济宁某矿区2009年12月到2010年02月期间更为精确的矿区地面沉降结果,并对沉降结果进行详细比较和系统分析.理论推导、相位模拟和真实数据实验都表明,相对于C波段的雷达干涉数据而言,L波段雷达干涉数据具有较强的保相能力,能够更好地降低失相干和相位不连续性的影响,更容易监测到沉降梯度和沉降量较大的矿区地面沉降,但对微小矿区地面沉降的敏感程度较低.     

黏弹TTI介质中旋转交错网格高阶有限差分数值模拟

以Carcione黏弹各向异性理论为基础,给出了适用于黏弹性具有任意倾斜对称轴横向各向同性介质(黏弹TTI介质)的二维三分量一阶速度-应力方程,采用旋转交错网格任意偶数阶精度有限差分格式求解该方程,并推导出了二维黏弹TTI介质完全匹配层(PML)吸收边界条件公式和相应的旋转交错网格任意偶数阶精度有限差分格式,实现了该类介质的地震波场数值模拟.数值模拟结果表明:该方法模拟精度高,边界吸收效果好,可以得到高精度的波场快照和合成记录;并且波场快照和合成记录能较好地反映地下介质的各向异性特征和黏弹性特征.     

一类测井响应方程存在的必要条件

本文将测井分析中所用的绝大多数响应方程归并为一类,给出了这类方程的一般形式。同时,根据一个基本事实导出了这类方程存在的必要条件,从而建立了检验和评价它们的理论标准,并已得到了实验验证。