一阶弹性波交错网格时间高阶差分格式及稳定性分析

弹性波模拟或逆时偏移时,对空间偏导数采用高阶差分格式可提高计算精度,但这种算法的稳定性条件过于严格,要求差分离散的时间步长必须足够小以确保算法稳定。在常规空间高阶差分格式的基础上,将速度(应力)对时间的高阶导数转化为不同精度的应力(速度)对空间的差分,得到了一种新的基于交错网格的时间高阶、空间高阶差分格式。通过对交错网格时间高阶差分格式稳定性的分析,认为该算法的稳定性条件较常规算法宽松,在弹性波场的求解过程中可以采用更大的时间步长。     

覆盖层影响下典型地-井模型瞬变电磁法正演

地-井瞬变电磁法是用于地质找矿的有效方法之一,开展该方法的三维正演研究能对资料处理与解释提供帮助。在时间域有限差分算法的基础上,建立了三维地质模型和覆盖层模型,模拟了该模型下的地-井瞬变电磁响应;引入了参数EA,分析了方形低阻体和覆盖层的响应特征以及覆盖层影响因素。结果显示：低阻体位于发射场源下方时,其响应曲线呈现“双高一低”的极值特征,高峰值分别对应低阻体上下边界,低峰值对应低阻体中心;低阻体位于发射场源旁边时,其响应曲线呈现单极值特征,该特征与覆盖层响应特征相同。当覆盖层电阻率较低时,“趋肤效应”作用使覆盖层的响应强度强于低阻体;但当发射场源离钻孔足够远时,覆盖层影响可以忽略,低阻体响应强度更与覆盖层-低阻体电阻率的比值有关。孔中发射能有效降低覆盖层影响。     

非接触极化曲线视反应电位差区分矿体矿物成分及确定其边界的实验研究

本文用非接触极化曲线视反应电位差区分矿体矿物成分和确定矿体边界的基本原理，着重阐述矿物非接触极化曲线视反应电位差特征的实验研究方法和结果，并指出非接触极化曲线视反应电位差区分矿体矿物成分和确定矿体边界是有可能的。     

基于有限差分原理的预应力锚索抗滑桩改进计算方法

针对现有预应力锚索抗滑桩计算模型及其相应计算理论存在的问题,基于有限差分法的原理,提出了改进的锚索桩计算模型并进行了详细的理论推导。根据实际的施工和受力过程,该模型包括预应力施加和滑坡推力作用两个阶段,通过滑坡推力分布形式的改进和桩锚协调变形方程的修正来真实反映预应力锚索抗滑桩的主动式受力特点。同时利用Matlab编制成相应的计算程序进行算例验证,并同目前普遍采用的计算理论进行对比分析。结果表明,该改进计算方法及其相应的计算程序充分考虑了预应力施加这一锚索桩特有的过程所带来的影响,且在计算中不必将锚索拉力等效为在桩顶处施加的剪力和弯矩。另外还可对桩前滑面以上存在任意厚度滑体和任意复杂地基条件的情况进行合理地设计计算,相对于目前已有的计算理论,该改进方法更加合理、可靠。     

含初始缺陷地铁隧道的抗震动力学行为研究

针对含初始缺陷的既有地铁隧道受邻近穿越隧道施工影响的首都地铁隧道工程背景,基于有限差分方法,通过对既有地铁车站结构受邻近隧道穿越后的抗震性能研究,揭示其抗震动力响应机理。结果表明,在邻近隧道穿越影响下,含初始缺陷的既有隧道各薄弱控制点在地震作用下位移和应力动力时程的波动趋势具有相似性,说明衬砌结构是整体运动的,其位移、应力较穿越前均有所增加,边墙和仰拱在穿越后抗震性能影响最大,剪切破坏主要以车站结构周边土体为主,新线穿越后其剪切破坏分布范围有所增大,应进行注浆加固处理,与自身结构刚度相匹配,提高抗震能力。     

Summertime Mixing Heights At Vienna, Austria, Estimated from Vertical Soundings and by a Numerical Model

Mixing heights calculated by the Danish OML meteorological pre-processor are compared to those diagnosed from radio- and tether-sonde vertical potential temperature profiles. All methods give reliable estimates of noon mixing heights deduced from radiosoundings, especially when the boundary layer is fully convective. Differences are larger during convective conditions without a well-defined capping inversion in the radiosonde potential temperature profile or when OML calculates a mechanical mixing height. The OML model is also able to calculate the daily course of the mixing height as expected. The tethersonde-derived mixing heights are especially valuable during the morning rise of the elevated inversion. Modifications to all three methods to improve mixing-height predictions are discussed.     

我国氨的排放量和时空分布

根据我国的实际情况讨论了各种氨排放源的排放因子，计算了各种氨源的贡献，以动物的贡献最大，占52%。1993 年我国氨的排放量为12 Mt，排放量最大的是河南、山东和四川，排放密度最大的是上海。在北京、广西、广东、湖南、江西和山东测定了大气氨的浓度，北京大气氨浓度有明显的季节变化，春季平均为16 μg/m3，夏季最高为41μg/m3，冬季只有3.2 μg/m3。一般夜间氨的浓度高于白天。测定了氨的垂直分布，并计算了氨浓度的标高。     

The Neutral, Barotropic Planetary Boundary Layer, Capped by a Low-Level Inversion

The presence of a low-level, capping inversion layer will affect the height and structure of the planetary boundary layer (PBL). Results from models of varying levels of sophistication, including analytical, turbulent kinetic energy (TKE), second-order closure (SOC), large-eddy simulation (LES) and direct numerical simulation (DNS) models, are used to investigate this influence for the neutral, barotropic PBL. Predicted and observed profiles of stress and geostrophic departure components, and integral measures, such as the parameters of Rossby-number similarity theory, are compared for the KONTUR, Marine Stratocumulus, JASIN, Leipzig, Pre-Wangara and Upavon field experiments.Analytical models of the equilibrium value of inversion height z_i, which depend on the surface friction velocity u_*, and both the Coriolis parameter f and the free-flow Brunt-Väisälä frequency N, are found to give reasonable estimates of the PBL height. They also indicate that only the KONTUR and Marine Stratocumulus experiments were strongly influenced by N. More quantitative comparisons would require larger, more comprehensive datasets. The effects of the presence of a capping inversion on the profile structure were found to be insignificant for h_* = |f|z_i/u_* > 0.15.The simple analytical model performed quite well over all values of h_*; it predicted the profiles of the longitudinal stress component (in the direction of the surface stress) better than the lateral component. The more advanced models performed well for small values of h_* (for flow over the sea), but systematically underestimated the cross-isobaric angle for flow over land. These models predicted the profiles of the lateral stress component better than the longitudinal component. The profiles of the analytical model agreed with those of the advanced models when the constant eddy viscosity of the outer layer was increased.Agreement with DNS was achieved by increasing the eddyviscosity of the analytical model by a factor of 5.Zilitinkevich and Esau(2002, Boundary-Layer Meteorology, 104, 371–379)suggest that the neutral, barotropic values of A and B of Rossby-numbersimilarity theory are not universal constants, but depend on the ratio N/|f|. The dependence for A and B is calculated using the analytical model and TKE models. Over the sea (h_* 0.1; N/|f| 100, where we have used the Zilitinkevich-Esau relation to convert between h_* and N/|f|) there is agreement between the model predictions and observations; however over land where the equilibrium boundary-layer height is greater (h_* 0.35; N/|f| 10) the inconsistency between the advanced model predictions (TKE, SOC, LES, and DNS) and observations, as noted previously by Hess and Garratt, still exists. We attribute this disagreement to violations of the strict assumptions of steady, horizontally homogeneous, neutral, barotropic conditions implicit in the observations. At small values of z_i and a strongly stable background stratification (h_* 0.04; N/|f| 1000) both the TKE and analytical models predict that A and B depend significantly on h_*, however observations are unavailable to confirm these predictions. Zilitinkevich and Esau call this case the `long-lived near-neutral PBL', and state that it is found in cold weather at high latitudes.     

四川自动与人工观测气温的差异分析

利用四川135个站自动与人工第2年平行观测气温资料,就自动与人工观测的气温的差异、引起差异的原因进行了分析。结果表明:气温各项目的自动观测比人工观测平均偏高,平均差值基本在0.2℃以内;气温、日最高、日最低的自动与人工观测比较,差值在±0.2℃之间的分别占58.71%、51.58%、62.68%,自动与人工观测值一致的分别占15.07%、12.16%、14.78%;自动与人工观测气温之间的对比差值存在明显的日变化,无明显的季节性和地域性差异;两种观测体制在观测时间上、观测方式上、测温传感器安装位置上,以及感应元件和测温原理的不同、仪器误差、热滞效应,都会造成观测结果出现差异。     

两次秸秆焚烧导致辽宁重污染过程的气象条件对比分析

利用环境监测站大气污染物数据、地面自动气象站观测资料、L波段加密探空资料和0.125°×0.125°的EC再分析资料,结合MODIS遥感火点监测和HYSPLIT4后向轨迹模拟结果,对比分析了2015年11月8日和2016年11月5日的两次由于东北地区秸秆焚烧导致辽宁重污染天气过程的大气边界层特征、气象扩散条件和大气污染物输送来源等。结果表明:两次过程地面PM_(2.5)浓度均出现快速上升和下降,其中2015年11月8日重污染过程的污染强度较2016年11月5日强,且持续时间更长。2015年11月8日重污染过程的混合层高度较低,其上层的中性层结转变为逆温层结,抑制混合层高度的发展。同时低层冷平流不断侵入到暖平流下方,使得大气层结稳定性增强,维持时间较2016年11月5日重污染过程更长,低层下沉运动和黑龙江西南部、吉林西部污染物的远距离输送增强使得辽宁地面污染物浓度快速累积。而2016年11月5日重污染天气过程主要受深厚冷空气影响,东北地区西部污染物的区域输送和地面风场辐合是地面污染物浓度快速上升的主要原因。