基于免疫粒子群算法的马斯京根模型参数识别

目前在河道汇流计算中广泛应用的是马斯京根模型,该模型参数的优化求解是影响能否精确模拟实测水文过程的关键问题.以往的模型参数求解大多采用试错法、矩法等方法,计算过程繁琐,计算精度不高,影响到模型的模拟精度.针对此类问题,笔者将免疫原理与粒子群优化算法有机结合.提出了免疫粒子群优化新方法,有效解决了传统方法的计算结果精度不高的问题,并在马斯京根模型参数优化求解中得以应用,结果表明免疫粒子群算法的计算结果精度令人满意,为河道洪水演算方面研究提供了一种新的研究方法和研究模式.     

用美国暴雨水管理模型计算北京城区防洪排水

城市化的发展，迅速地改变着城镇的下垫面条件，城区排洪河道的大规模整治与改建，地下排水系统的管网化直接改变着原有雨洪径流的形成规律。根据这一特点，首先，简要介绍了美国暴雨水管理模型ＳＷＭＭ的计算及软件系统；其次，通过ＳＷＭＭ在北京太平湖及安定门两个试验小区的应用，认为模型在流域模拟、不同阶段水流运动结构等方面均较传统的计算方法合理，文中还对太平湖小区作了典型下垫面分析介绍，说明只要基本资料收集符合实     

岩土力学数值模拟中力学参数的确定方法

力学参数的选取是进行岩土力学数值模拟的关键，它直接决定了模拟结果的合理性。不同的力学模型需要选取不同的计算参数，经验值在参数的确定过程中具有一定程度的参考价值，几种变形参数和强度参数的现场尺度转化方法解决了室内土工实验参数的现场应用问题，提出了进行力学参数和计算模型合理确定的推断模型的概念，并对推断模型实现的思想和步骤进行了较为详细的讨论。     

基于GIS的水资源管理系统的开发和应用

以地区水资源总开发量最大为优化目标，以人工开发利用水资源量为决策变量，以地下水位为状态变量，将水资源系统的模拟模型作为约束条件集，嵌入到水资源系统的优化模型中，建立了水资源可持续利用量化模型。以水资源可持续利用量化模型为核心，在地理信息系统(GIS)支持下开发出城市水资源管理系统。该系统利用集成技术将多个应用软件的功能模块有机地组合在一起，以邯郸市为例完成了对城市水资源可持续利用的量化模拟、优化计算和科学管理。     

四川南江县窑厂坪滑坡地下水数值模拟

地下水流场的分布是评价滑坡体稳定性的关键因素。使用Visual MODFLOW地下水数值模拟软件,建立窑厂坪滑坡在天然、暴雨和排水条件下的三个模型。通过参数反演提高天然模型的精度。分析暴雨和排水模型表明,强降雨抬升滑坡区地下水位,并在部分区域形成雍高水位,而廊道地下排水工程降低地下水水位,工程效果明显。     

泥石流防治工程优化设计方案初探

本文介绍一种以可靠理论为指导的泥石流防治工程优化设计方法，文中首先提出了由于泥石流荷载的不确定性而以概率约束表示的泥石流防治工程参数优化设计的模型；然后建立了用离散变量优化与随机模拟相结合的一种求解方法；最后应用实例对该方法的原理和计算程序进行了说明。     

三种水文模型不确定性分析方法比较

模型不确定性研究是水文科学的重要课题。以尼泊尔Bagmati流域为案例,采用了马尔科夫链蒙托卡罗(Markov Chain Monte Carlo)、蒙托卡罗(Monte Carlo)和拉丁超立方体(Latin Hypercube)等三种方法,分析了水箱模型输出成果的不确定性,并将三种方法所获得参数不确定性进行了比较。另外,运用Meta-Gaussian模型计算了总体不确定性,在基于所采用的似然函数基础上,对由参数导致模型输出的不确定性和模型输出的总体不确定性进行了比较。结果显示,模型的不确定性比参数不确定性更为重要,同时也表明,尽管蒙托卡罗和拉丁超立方体两种模拟方法产生几乎相同的结果,但两者都与马尔科夫链蒙托卡罗方法有很大的不同。     

地下排水对滑坡稳定性影响动态研究

提出了考虑孔隙水压力的滑坡稳定性计算方法,对工程中常用的不平衡推力法进行改进。以某个水电工程滑坡为例,采用FLAC3D进行三维数值分析,对滑坡在地下隧洞排水过程中的孔隙水压、水位及流量等变化进行了模拟。利用FLAC3D计算结果,并使用FISH语言编程,得到改进后的不平衡推力法的所需数据,对滑坡的稳定性变化进行了动态研究。计算结果表明,上述方法很好地反映了地下排水过程中稳定性增长过程,同时表明常规的计算方法低估了地下排水在提高滑坡稳定性和降低剩余下滑力方面的作用。     

新安江产流模型与改进的BP汇流模型耦合应用

为提高新安江模型的汇流计算精度并减少经验因素对参数率定的影响,将新安江产流模型与改进的BP汇流模型相耦合,建立XBK(XAJ-BP-KNN)模型。该模型以前期模拟流量和新安江产流模型计算的产流量作为BP网络的输入,出口断面流量作为网络输出,拟合汇流的非线性关系,代替新安江模型的分水源、线性水库及河道马斯京根法的汇流计算;采用相似原理和K-最近邻算法,基于历史样本的模拟误差及相应影响要素对网络输出进行误差修正,实现了无前期实测流量的连续模拟;模型使用SCE-UA算法与遗传早停止LM算法相结合的全局优化方法进行参数优选。在呈村流域的验证表明XBK模型的模拟精度高于新安江模型,全局优化方法能找到最优参数,降低了模型的使用难度。     

降雨对滑坡的作用机理及动态模拟研究

从动态分析和数值模拟两方面讨论了滑坡过程中的地下水作用机理和地下水动力场的时空分布问题，建立了二维极限平衡状态下的临界水位高度模型和描述滑坡剖面二维不稳定流动问题的动力学模型；以临界水位高度模型对重庆土台镇滑坡稳定性分析的研究结果表明，该滑坡的防治措施宜首先考虑地表水防渗和地下排水方案，在滑坡位移监测的同时更应考虑位移、降雨量、地下水动态三者的结合；采用数值模型计算了长江三峡库区鸡扒子滑坡因降雨入渗而引起的地下水动力场的时空分布，研究结果表明，地下水水位和流速场的时空演化过程与滑坡发生的时空特点基本吻合，当降雨达到46h时，滑坡体中地下水流向与滑动面倾向基本一致，尤其是在滑体底部出现潜水流，此时的地下水流动对滑坡稳定性最为不利。     

Geochemistry of Thermal Waters of Continental Margin of Far East of Russia

Studied waters belong to warm(T=30-50℃),alkaline(pH=8.9-9.3),low mineralized(TDS235 mg/1)Na-HCO_3 or Na-SO_4-HCO_3 thermal waters with high content of SiO_2(up to 81 mg/l)and F(up to 3.9 mg/1),occur on modern volcano-tectonic rejuvenated areas of Eastern Sikhote-Alin orogenic belt.Low~3He concentration as well as N_2/O_2 and N_2/Ar ratios exclude influence of deep mantle fluid.New rare earth element data constrain our understandmg of water-rock interaction occurring in the water source region.Meteoric origin of waters is proved by stable isotope values varying from-71‰to-136.1‰and from-10.8‰to-18.8‰forδ~2U andδ~(18)O respectively.REE patterns reflect high pH,resultfing from water-rock interaction and oxidative conditions.Calculations of deep aquifer temperature using Na-K and quartz geothermometers show 116.8-131.1°C and 82.2-125.8℃respectively.Presence of deep faults both with abnormal thermal gradient(~45-50 K/km)define unique geochemical shape of thermal waters of Sikhote-Alin,area,where no present volcanic activity is registered.     

黄土路基沉降曲线特征分析

在分析郑西客运专线试验段路基沉降实测资料的基础上,提出沉降变形常规预测方法受观测时间及单次观测结果的影响比较大,不利于用来预测最终沉降。通过分析现场沉降变形与时间关系曲线,提出路基沉降的发展主要受固结的影响,其发展是一个从沉降增加到趋缓的过程,可以明显分为固结变形及流变变形,并用非饱和土等效固结理论加以解释,提出了一种利用s- 曲线特征预测最终沉降量的方法,经检验对比发现其预测结果受观测时间及单次观测结果的影响比较小,优于其它常规预测方法     

沉管隧道基槽爆破施工对既有堤岸稳定性影响的数值仿真分析

结合爆破数值模拟中爆破荷载的特点,分析了各种荷载施加方法的特点及爆破震动场模拟的要求,基于圣维南荷载等效原理,提出了爆破荷载的等效施加方法,即将作用在炮孔壁上的爆破荷载等效后,施加在同排炮孔中心线或者面上来模拟爆破荷载的作用。结果表明,等效施加方法在近区存在一定差异,但在中远区吻合得较好。利用该等效荷载施加方法,采用三维实体模型计算分析了某沉管隧道基槽爆破施工对既有堤岸稳定性影响,计算了周围不同位置、不同深度位置处的振动响应,并在此基础上,分析研究了减小爆破振动对既有结构物影响的减振措施,可为类似工程提供借鉴。     

岩石坚硬程度对核电厂地基-基础地震响应特征的影响分析

考虑核电厂地基-基础的动力相互作用,应用显式动力有限差分法分析了地震作用下极软岩、较软岩、坚硬岩上核电厂建筑结构基础的地震响应特征,比较了岩石坚硬程度对基础加速度反应谱的影响。研究表明：随着岩石坚硬程度的提高,核电厂建筑物结构基础的地震响应有增加的趋势;在周期轴上,基础处的加速度反应谱曲线会随着岩石坚硬程度的提高逐渐向短周期（高频段）方向移动。在高频段,建造于较坚硬岩石上基础结构的加速度反应谱值偏大;在中等频段,建造于较软岩石上基础结构的加速度反应谱值偏大;在低频段,岩石坚硬程度对加速度反应谱的影响不显著。     

Strains of scattering of near-field of a point source

Three dimensional scattering of near-field is studied for dilatation and rotation in the time domain. The perturbation method is applied to solve the equation of motion for the first order scattering from a weak inhomogeneity in an otherwise homogeneous medium. The inhomogeneity is assumed close enough to the point source so that the near-field intermediate wave is dominating over the far-field sphericalP andS pulses. The integral expressions are derived to relate dilatation and rotation of scattering to the radial fluctuations of velocities and density in the inhomogeneity. These integrals are solved to calculate the strains of scattering from (a part of) an inhomogeneous spherical shell of arbitrary curvature. Variable curvature may allow the shape of inhomogeneity volume element to change uniformly from spherical to rectangular. Rotation of scattering from a spherical shell is independent ofP wave velocity inhomogeneity. Dilatation of scattering does not involveS wave velocity inhomogeneity but its gradient. The back scattering results are obtained as a special case. Strains are computed numerically, for hypothetical models to study the effects of various parameters viz., velocity inhomogeneity, distance of source from inhomogeneity and from receiver, and thickness of inhomogeneity. The curvature of the spherical shell is varied to study the effects of the shape of inhomogeneous volume element on scattering.     

Structure of sediments of kaolinite

In an assembly of clay particles placed in a fluid, each particle is typically subjected to: (1) double-layer repulsive forces; (2) van der Waals attractive forces; and (3) contact mechanical forces. The study presented here concerns an approximate, quantitative analysis of clay suspensions, with considerations to the first two - the physico-chemical forces. Using recent theories to calculate the physico-chemical forces between two clay particles in an approximate model of an assembly, the equilibrium void ratio of a clay suspension is computed. The mechanical forces are ignored in the analysis. The results serve to verify the validity of physico-chemical theories employed and help interpret experimental data more fundamentally in terms of the system variables.     

土石坝心墙水力劈裂分析

简要分析了水力劈裂的发生条件、力学机理。基于水压楔劈机理,利用有限元方法建立了一种水力劈裂的发生判定方法。该方法假定心墙预先存在局部渗透弱面（裂缝）,通过将裂缝位置的单元材料改为裂缝软材料,考虑库水进入裂缝后对裂缝周围土体的作用,建立水力劈裂分析的平面有限元模型,确定裂缝端部垂直于裂缝面的正应力,进而依据该正应力判断水力劈裂发生的可能性,该方法同时可模拟水力劈裂的发展过程。     

土钉支护失稳分析

某基坑采用复合土钉墙支护,开挖至坑底后,基坑北侧发生失稳事故。本文就其事故发生原因作一探讨。     

地震次生地质灾害风险评估集对态势分析

5．12四川汶川8级强震引发了大量次生地质灾害,加剧了灾情并严重影响抗震救灾。地震次生地质灾害受多重因素制约,其孕育与致灾过程具有不确定性。利用集对分析原理与方法,对地震次生地质灾害进行同异反态势分析,拟定了风险评估的集对分析同一度、差异度、对立度等指标体系的构建原则与赋值标准,对承灾体系统的不确定性及其作用作了刻画与分析,建立了不同风险分区代表性的集对分析联系度表达式,为地震灾害研究提供了可资借鉴的新思路和技术方法。