长时效的正压原始方程能量完全守恒(拟)谱模式

遵循误差反演补偿新计算原理,对正压原始方程传统气象全球拟谱模式方案进行了改造,构造了正压原始方程能量完全守恒全球拟增模式新计算方案,解决了正压原始方程的（非线性）计算稳定性问题和能量守恒整体性质保持问题,改进了相应正压原始方程传统气象全球拟谱模式方案的计算效能。新方案的数值试验表明：在计算实践上,新方案在解决能量守恒问题的同时,可解决（非线性）计算稳定性问题,并在一定条件下可解决非线性计算收敛性问题。进一步的比较数值试验还表明：在计算实践上,新方案具有在提高相应传统气象方案的计算精度,减少其计算量的同时,延长其计算时效,解决其中一类特定“气候漂移”问题方面的效用。本工作原理也适用于斜压原始方程情形。     

An improvement of the mass flux convection parameterization scheme and its sensitivity tests for seasonal prediction over China

A modified cumulus parameterization scheme, suitable for use in a seasonal forecast model, is presented. This parameterization scheme is an improvement of the mass flux convection scheme developed by Gregory and Rowntree (1989; 1990). This convection scheme uses a “bulk” cloud model to present an ensemble of convective clouds, and aims to represent shallow, deep, and mid-level convection. At present,this convection scheme is employed in the NCC T63L20 model (National Climate Center, China Meteorological Administration). Simulation results with this scheme have revealed some deficiencies in the scheme,although to some extent, it improves the accuracy of the simulation. In order to alleviate the deficiencies and reflect the effect of cumulus convection in the actual atmosphere, the scheme is modified and improved.The improvements include (i) the full estimation of the effects of the large-scale convergence in the lower layer upon cumulus convection, (ii) the revision of the initial convective mass flux, and (iii) the regulation of convective-scale downdrafts. A comparison of the results obtained by using the original model and the modified one shows that the improvement and modification of the original convection scheme is successful in simulating the precipitation and general circulation field, because the modified scheme provides a good simulation of the main features of seasonal precipitation in China, and an analysis of the anomaly correlation eoetfieient between the simulation and the observations confirms the improved results.     

声波方程频率域高精度正演的17点格式及数值实现

频率域正演计算是频率域全波形反演的基础.传统的最优9点格式只具有二阶精度,不能满足高精度地震成像的需要.本文考虑两个四阶精度的格式,即经典的四阶9点格式和优化的17点格式.17点格式可将最小波长内所需网格点数减小到2.56.通过在简单模型和Overthrust模型上的数值实验,比较分析了三种格式的正演效果;简单模型数值实验显示了17点格式克服频散误差的能力优于四阶9点格式和最优9点格式;复杂模型数值实验则进一步承认了算法的可行性.     

基于贝叶斯框架的各向异性页岩储层岩石物理反演技术

页岩岩石物理建模旨在建立页岩矿物组分、微观结构、流体填充与岩石弹性参数的关系.对四川盆地龙马溪组页岩进行岩石物理建模研究,针对页岩黏土含量高、层间微裂缝发育等特点,利用Backus平均理论描述页岩黏土矿物弹性参数,利用Chapman理论计算与水平微裂缝有关的VTI各向异性,并利用Bond变换考虑地层倾角的影响.提出以黏土矿物纵、横波速度和孔隙纵横比为拟合参数进行岩石物理反演的方法,并引入贝叶斯框架减小反演的多解性.由已知的黏土矿物纵、横波速度和孔隙纵横比作为先验信息,并以测井纵、横波速度作为约束条件建立反演的目标函数,同时利用粒子群算法进行最优化搜索.计算结果表明,基于先验约束和粒子群算法的反演方法能够较准确地反演黏土矿物的弹性参数、孔隙形态参数以及裂缝密度等参数.计算得到的黏土纵、横波速度较高,并且在一定范围内变化,这可能与龙马溪组页岩的黏土矿物组分中具有较高弹性模量的伊利石含量较高有关,同时也与黏土定向排列等微观物性特征有关.反演得到的裂缝密度与纵波各向异性参数ε呈明显的正相关,而与横波各向异性参数γ相关性较小.另外,页岩各向异性参数与黏土垂向的纵横波速度有较强的相关性.     

利用密钥策略进行数据分离的算法

在分析现有的数据处理策略的基础上,将经过改进的基于二进制Goppa编码的McEliece密码体制与数据份额的产生相结合,提出了一种基于纠错码密码体制的数据分离算法(KPDA)。与经典秘密共享算法相比,该算法具有同等或更高的效率;而与同类的数据分离技术相比,基于编码的密码体制可以提供更好的安全保障。     

Non-point source pollution estimation in the Pingqiao River Basin,China, using a spatial hydrograph-separation approach

For efficient and targeted management, this study demonstrates a recently developed non-point source (NPS) pollution model for a year-long estimation in the Pingqiao River Basin (22.3 km²) in China. This simple but physically reasonable model estimates NPS export in terms of land use by reflecting spatial hydrological features and source runoff measurements under different land-use types. The NPS export was separately analysed by a distributed hydrological model, a spatial hydrograph-separation technique, and an empirical water quality sub-model. Simulation results suggest that 57 890 kg of total nitrogen (TN) and 1148 kg of total phosphorus (TP) were delivered. The results, validated with observed stream concentrations, show relative errors of 23.3% for TN and 47.4% for TP. Countermeasures for urban areas (5.3% of total area) were prioritized because of the high contribution rate to TN (14.1%) and TP (26.2%) which is caused by the high degree of runoff (8.5%) and pollution source.     

Numerical modeling of lock-exchange gravity/turbidity currents by a high-order upwinding combined compact difference scheme

This study presents two-dimensional direct numerical simulations for sediment-laden current with higher density propagating forward through a lighter ambient water.The incompressible NavierStokes equations including the buoyancy force for the density difference between the light and heavy fluids are solved by a finite difference scheme based on a structured mesh.The concentration transport equations are used to explore such rich transport phenomena as gravity and turbidity currents.Within the framework of an Upwinding Combined Compact finite Difference(UCCD)scheme,rigorous determination of weighting coefficients underlies the modified equation analysis and the minimization of the numerical modified wavenumber.This sixth-order UCCD scheme is implemented in a four-point grid stencil to approximate advection and diffusion terms in the concentration transport equations and the first-order derivative terms in the Navier-Stokes equations,which can greatly enhance convective stability and increase dispersive accuracy at the same time.The initial discontinuous concentration field is smoothed by solving a newly proposed Heaviside function to prevent numerical instabilities and unreasonable concentration values.A two-step projection method is then applied to obtain the velocity field.The numerical algorithm shows a satisfying ability to capture the generation,development,and dissipation of the Kelvin-Helmholz instabilities and turbulent billows at the interface between the current and the ambient fluid.The simulation results also are compared with the data in published literatures and good agreements are found to prove that the present numerical model can well reproduce the propagation,particle deposition,and mixing processes of lock-exchange gravity and turbidity currents.     

BDA方案在台风路径预报中的应用

利用PSU／NCAR中尺度非静力有限区域MM5及其伴随模式,以T106分析资料为背景场,设计两种台风Bogus方案对台风的初始场进行优化,并进行了数值模拟试验。对9608号台风个例的数值模拟试验研究表明,经过优化的台风初始场较好地改进了由于海洋上资料缺乏所造成的台风中心位置不准、台风环流偏弱和台风内部结构不完整等问题,提高了台风路径预报的准确率。通过试验对比发现,BDA方案优化的初始场更合理,其台风路径预报效果优于GFDL方案。     

白鹤滩水电站左岸坝肩开挖边坡稳定性分析

白鹤滩水电站左岸坝基边坡地质构造错综复杂,坡内多条陡倾断层和顺坡向错动带是控制边坡稳定性的主要因素。结合现场地质资料、监测数据及勘察认识,利用通用离散元软件UDEC建立左岸边坡开挖变形分析模型。首先在无支护开挖工况下研究主要结构面F17、LS331、LS3319上下盘岩体关键点应力及位移变化,并指出潜在失稳破坏区;同时采用预应力锚索单元在不同施工时序下对潜在不稳定块体进行加固计算,对比研究了不同方案下左岸边坡的变形破坏机制和整体稳定性。计算结果和现场监测表明：左岸边坡变形破坏与结构面密切相关,在开挖卸荷作用下主要表现为沿错动带的压剪破坏和陡倾断层的张拉破坏,采用预应力锚索加固处理能够有效提高边坡稳定性,并且在开挖后及时跟进支护可以有效抑制边坡剪切变形。研究结果对白鹤滩左岸坝基边坡后续工程加固和施工工序具有一定的参考价值。     

Well-balanced numerical modelling of non-uniform sediment transport in alluvial rivers

The last two decades have witnessed the development and application of well-balanced numerical models for shallow flows in natural rivers.However,until now there have been no such models for flows with non-uniform sediment transport.This paper presents a 1D well-balanced model to simulate flows and non-capacity transport of non-uniform sediment in alluvial rivers.The active layer formulation is adopted to resolve the change of bed sediment composition.In the framework of the finite volume Slope Llmiter Centred(SLIC) scheme,a surface gradient method is incorporated to attain well-balanced solutions to the governing equations.The proposed model is tested against typical cases with irregular topography,including the refilling of dredged trenches,aggradation due to sediment overloading and flood flow due to landslide dam failure.The agreement between the computed results and measured data is encouraging.Compared to a non-well-balanced model,the well-balanced model features improved performance in reproducing stage,velocity and bed deformation.It should find general applications for non-uniform sediment transport modelling in alluvial rivers,especially in mountain areas where the bed topography is mostly irregular.