发展方程计算稳定性的比较分析

对线性与非线性发展方程的几种差分格式,以一维线性和非线性平流方程为例,对线性与非线性发展方程差分格式的计算稳定性进行了比较分析,揭示了差分格式结构和初值形式与非线性计算稳定性的关系。     

Dynamic and Numerical Study of Waves in the Tibetan Plateau Vortex简

In terms of its dynamics, The Tibetan Plateau Vortex （TPV） is assumed to be a vortex in the botmdary layer forced by diabatic heating and friction. In order to analyze the basic characteristics of waves in the vortex, the governing equations for the vortex were established in column coordinates with the balance of gradient wind. Based on this, the type of mixed waves and their dispersion characteristics were deduced by solving the linear model. Two numerical simulations with triple-nested domains--one idealized large-eddy simulation and one of a TPV that took place on 14 August 2006---were also carried out. The aim of the simulations was to validate the mixed wave deduced from the governing equations. The high-resolution model output data were analyzed and the results showed that the tangential flow field of the TPV in the form of center heating was cyclonic and convergent in the lower levels and anticyclonic and divergent in the upper levels. The simulations also showed that the vorticity of the vortex is uneven and might have shear flow along the radial direction. The changing vorticity causes the formation and spreading of vortex Rossby （VR） waves, and divergence will cause changes to the n~otion of the excitation and evolution of inertial gravity （IG） waves. Therefore, the vortex may contain what we call mixed ：inertial gravity-vortex Rossby （IG-VR） waves. It is suggested that some strongly developed TPVs should be studied in the future, because of their effects on weather in downstream areas.     

地震海啸的激发与传播

本文围绕海啸的激发和传播两个方面简述了有限矩形源产生的地表形变场、浅水波浪理论和Boussinesq方程等基本理论,讨论了震源参数对海啸传播的影响,比较分析了海啸的数值模拟方法,介绍了我国地震海啸的研究进展等。     

Receiver function estimated by maximum entropy deconvolution

Maximum entropy deconvolution is presented to estimate receiver function, with the maximum entropy as the rule to determine auto-correlation and cross-correlation functions. The Toeplitz equation and Levinson algorithm are used to calculate the iterative formula of error-predicting filter, and receiver function is then estimated. During extrapolation, reflective coefficient is always less than I, which keeps maximum entropy deconvolution stable. The maximum entropy of the data outside window increases the resolution of receiver function. Both synthetic and real seismograms show that maximum entropy deconvolution is an effective method to measure receiver function in time-domain.     

弹性介质中密度ρ与拉梅常数λ、μ的衍射层析成像方法研究

利用接收到的纵波信息，通过解纵波方程得到介质的密度和拉梅常数分布函数．首先，在频率波数域建立了描述入射纵波信号、接收纵波信号和介质三参数扰动量之间关系的积分方程．然后，建立了关于接收信号的各频率成份的方程组．经过正反傅氏变换求解此方程组即可得到介质三参数的分布函数．最后，就跨孔观测系统进行了数值模拟成像试验。结果表明方法有效．     

用于检测同相轴的Duffing型系统恢复力项的讨论

在用于检测同相轴的Duffing型系统中,非线性恢复力项的选取至关重要.针对具有不同非线性恢复力项的混沌系统,本文首先就非线性恢复力项的来源进行了说明,通过对哈密顿系统和耗散系统的仿真实验分析,以及理论推导,确定了势能—哈密顿多项式方程中可用于混沌振子检测系统的各项幂次和各项系数组合关系,即非线性恢复力项的选取规则:幂次最高项为奇次,其系数大于零.     

含CO2盐水流体包裹体摩尔体积和组分求解新方法

借鉴Parry(1986)方法的思路, 提出求解含CO₂盐水流体包裹体摩尔体积和组分(V_m-X)的新方法. 新方法以实测包裹体气-液相CO₂部分均一温度及均一方式(T_h,CO₂), 包裹体盐度(S)和包裹体完全均一温度(T_h)为原始数据, 构建了含X_CO2、X_NaCl、V_m及F(包裹体气-液相CO₂部分均一时CO₂相的充填度)四个未知量的四个关联方程. 通过解四个方程构成的方程组, 求取包裹体的V_m-X值. 前三个方程为X_CO2、X_NaCl和V_m的数学表达式, 它们只与T_h,CO₂、S和F相关, 其简化的形式可表示为: ; ; . 第四个方程为包裹体完全均一时X_CO2、X_NaCl、V_m和T_h间的热力学关系式, 简化形式为: . 解方程组要使用迭代求解法, 过程如下: 先给定F值代入前三个方程, 可分别求得X_CO2、X_NaCl和一个摩尔体积值V_m1, 然后把求得的X_CO2、X_NaCl代入方程f₄ 求出另一个摩尔体积值V_m2, 当V_m1＝V_m2, V_m1(V_m2)、F、X_CO2和X_NaCl即为整个方程组的解, 如符合地质意义, 即求得了包裹体的V_m-X值. 与Parry(1986)方法相比, 该方法更易于使用, 对X_CO2的求解也更精确. 新方法适用于求解CO₂气-液相部分均一时, 温度高于笼合物熔化温度, 且不含固相石盐的含CO₂盐水体系流体包裹体.     

基于环境星HJ1B数据反演大亚湾海表温度的三种算法

利用大亚湾2011年2月10日一景环境星HJ₁B热红外遥感影像,通过MODTRAN计算得到大气上、下行辐射和大气透过率,再运用经验公式计算星下辐亮度,分别用大气辐射传输方程模型法、覃志豪单窗算法和Jimenez-Munoz&Sobrino单通道算法反演大亚湾地区海表温度。将三种方法反演的结果与MODIS海表温度对比,其中覃志豪普适性单窗算法反演的海表温度与MODIS海表温度一致性最好,获得的精度较高,RTM算法与覃志豪单窗反演的结果平均误差在0.7K左右。     

基于“S”型函数的隧道开挖变形预测模型

基于不同开挖阶段围岩阶段变形曲线几何特征,提出一种基于"S"型函数的隧道开挖全过程变形预测模型,实现隧道开挖全过程中的围岩变形预测。对比分析了常用"S"型函数适应性,建立了基于Logistic模型的隧道围岩变形分析模型,模型精度控制在2%以内,并提出了模型参数辨识过程。利用Logistic隧道围岩变形模型,讨论了隧道开挖瞬时变形特征,提出了一种求解隧道瞬时开挖变形值的方法,计算出不同工程条件下瞬时变形值处于2. 192～43. 288 mm,瞬时位移占最终位移的6. 61%～18. 33%。以马尾山隧道为工程实例,建立了具体的工程隧道围岩变形模型,得到隧道开挖瞬时位移值处于1. 041～9. 475 mm,瞬时位移占总位移比值处于5. 52%～36. 44%。     

利用古成熟度方程初步恢复济阳坳陷晚三叠世印支运动剥蚀量

与新生界镜质组反射率(Ro)值比较,济阳坳陷晚古生界Ro值整体上出现了显著的异常,通过对济阳坳陷大量的晚古生界Ro异常数据点与新生界Ro数据点的分析、比较,建立了济阳坳陷新生界成熟度方程（RoH方程）和晚古生界古成熟度方程,利用这些古成熟度方程,对济阳坳陷Mz/Pz间断面剥蚀量进行了初步恢复,研究结果表明东营凹陷、沾化凹陷、车镇凹陷Pz/Mz间断面剥蚀量分别达31~41 km、20~41 km、31~41 km左右,平均在35 km左右,佐证了华北地台在印支运动开始前是个统一的大型沉积盆地,济阳坳陷部分晚古生界煤系地层在晚三叠世初期之前已经开始热解生烃。