Flow time estimation with spatially variable hillslope velocity in ungauged basins

The automated spatial estimation of the hillslope runoff dynamics is used as a valuable tool for the estimation of the travel time distribution (flow time), a major factor for the hydrologic prediction in ungauged basins. In fact, while the flow time function is usually obtained by rescaling the flow paths with constant channel and hillslope velocities, in this work a spatially distributed kinematic component, as a function of terrain features and in particular slope and land use, is implemented and its influence on the hydrologic response is tested by means of the Width Function Instantaneous Unit Hydrograph (WFIUH) framework. Hillslope surface flow velocities are evaluated by applying different uniform flow formulas within an automated DEM-based (terrain analysis) algorithm. A comparison test of the performances of the Manning, Darcy, Maidment and Soil Conservation Service uniform flow equations is performed for several case studies in Italy pertaining to different climatic and geomorphic conditions. Results provide new insights for a better understanding of the flow time function also introducing a more parsimonious and physically-based calibration scheme of the WFIUH.     

基于双曲正切函数的土动力非线性本构模型

基于双曲正切函数,构造了一种新的土动力非线性本构模型。该模型从双曲正切函数出发,通过平移和缩放构造了土动力应力应变关系的骨架曲线和卸载、反向加载的滞回曲线。相对于其他模型而言,这种新的数学模型参数较少,物理概念明确,形式简单,更利于编程实现,加卸载应力应变关系曲线只需要记忆加卸载转折点处的应力和应变值。将该数学模型应用于场地地震反应分析中,并与常用的等效线性化方法的地震反应结果进行比较,结果初步证明了该本构模型的可行性。     

江西数字地震台网近震震级偏差与量规函数关系的研究

采用江西数字地震台网数字地震记录资料,选取自2007年10月正式运行至2009年12月所记录的138次ML≥1.5的地震事件,计算了各台站单台震级与台网平均震级的偏差、平均偏差和标准偏差,从震级偏差频次分布、量规函数、平均周期等因素对单台震级的影响上作了定量统计和分析。结果表明:①除九江、赣州、石城、都昌、高安5个台的震级与台网的平均震级偏差比较大外,其余台站均比较小;②目前所用量规函数在震中距小于40 km时偏小,大于330 km时偏大,给出了量规函数校正值,③平均周期小于0.17 s时,变得极不稳定,离散程度明显增大,这种情况应该是由于单台震中距超过300 km的样本数太少所导致的。     

用接收函数方法研究上海地震台阵下地壳结构

利用上海地震台阵16个台站记录的远震资料,采用接收函数线性反演方法,对台阵下的地壳速度结构进行研究,获得了研究区域内地壳厚度和地壳速度的分布特征。研究结果表明,研究区域Moho面深度约为33±2 km,Moho面深度基本不变,地幔顶部S波速度约4.4 km/s,地壳内没有发现明显的低速层。     

Stationary phase approximation in the ambient noise method revisited

The method of extracting Green's function between stations from cross correlation has proven to be effective theoretically and experimentally.It has been widely applied to surface wave tomography of the crust and upmost mantle.However,there are still controversies about why this method works.Snieder employed stationary phase approximation in evaluating contribution to cross correlation function from scatterers in the whole space,and concluded that it is the constructive interference of waves emitted by the scatterers near the receiver line that leads to the emergence of Green's function.His derivation demonstrates that cross correlation function is just the convolution of noise power spectrum and the Green's function.However,his derivation ignores influence from the two stationary points at infinities,therefore it may fail when attenuation is absent.In order to obtain accurate noise-correlation function due to scatters over the whole space,we compute the total contribution with numerical integration in polar coordinates.Our numerical computation of cross correlation function indicates that the incomplete stationary phase approximation introduces remarkable errors to the cross correlation function,in both amplitude and phase,when the frequency is low with reasonable quality factor Q.Our results argue that the distance between stations has to be beyond several wavelengths in order to reduce the influence of this inaccuracy on the applications of ambient noise method,and only the station pairs whose distances are above several (5) wavelengths can be used.     

Turbulent intensity and its similarity function over an Inner Mongolian grassland during spring

Macroscopic reason for the development of turbulences in fluid is the conjunction action of velocity and temperature shears. We use the 2008 spring flux observations in the surface layer over Inner Mongolian grassland, obtain the similarity functions of turbulent intensity under different stratification stability conditions, and calculate the clear day turbulent intensity in the spring, which is then analyzed and compared with the observed values of turbulent intensity. The results show that in the Inner Mo...     

台风“麦莎”移动过程中水汽输送特征

利用NCEP1°×1°的再分析资料,用大尺度水汽通量流函数和速度势以及水汽收支对台风“麦莎”在移动过程中水汽输送流函数和速度势进行了分析。结果表明：台风“麦沙”在北上过程中,主要有2支水汽通道,一支是从太平洋开始向西经赤道到达印度洋,在索马里转向,经过孟家加拉湾呈西南气流向西北方向输送;另一支来自副高南侧偏东气流。计算分析还表明,用无辐散风流场来定台风中心要比用总的流场更精确。     

改进的模拟退火遗传算法在地下水管理中的应用

对于高度非线性、非凸的地下水管理模型,传统优化方法难以找到全局最优解。本文采用模拟退火遗传算法求解地下水管理模型,并从三个方面对算法进行改进:引入小生境技术,采用自适应交叉和变异概率,在选择过程中采用最优保存策略,从而提高算法的全局寻优能力和收敛速度。采用惩罚函数法处理约束条件。用Fortran 90语言编制了计算程序,并通过Schaffer测试函数验证了该算法不仅具有强大的全局寻优能力和局部搜索能力,而且具有较快的收敛速度和较高的优化精度。将该算法应用到某研究区地下水管理中,取得了较好的效果。     

基于核最小二乘模型的矿产靶区预测

地质统计单元的含矿性与地质找矿证据之间存在复杂的非线性关系,建立这种复杂关系的多元非线性统计模型并预测矿产靶区,对矿产勘查具有重要指导意义。以核函数为理论工具,在核最小二乘原理基础上提出了矿产靶区预测的核最小二乘模型;在GDAL数字图像输入输出函数库和CLAPACK线性代数软件包基础上,用VC++语言开发了面向栅格数据的矿产靶区预测核最小二乘模型算法程序,并把模型应用于新疆阿勒泰地区的矿产靶区预测研究。在MapInfo中生成包含100×151个网格统计单元的栅格图层,把栅格化后的15种找矿证据图层转化成100×151×15的数字图像数据立方体,用自行开发的程序计算每个网格统计单元的核最小二乘判别得分。结果表明,网格统计单元判别得分的高值区与已知矿床（点）的空间分布基本一致。     

基于粒子群算法与混合罚函数法的有限元优化反演模型及应用

岩土工程优化反分析是一个典型的复杂非线性函数优化问题,采用全局优化算法是解决这个问题的理想途径。针对常规反演方法应用于岩土工程参数反演时搜索效率低的缺点,结合粒子群算法和遗传算法的特点,充分考虑二者的互补性,提出一种效率较高的全局优化算法,以测点的实测值与计算值建立一种新的评价函数,将多目标优化问题转化为单目标优化问题,用混合罚函数法将约束问题变为无约束问题,构建了一种新的目标函数,将有限元程序ABAQUS作为一个模块嵌入到优化算法程序中,编制了有限元优化反演分析程序。并给出了应用实例验证了该法的有效性和实用性,是一种可行的参数反演方法,可应用于实际工程中复杂岩土介质初始应力场反演、渗流场以及位移反分析