Shuffled Complex Evolution model calibrating algorithm: enhancing its robustness and efficiency

Shuffled Complex Evolution—University of Arizona (SCE‐UA) has been used extensively and proved to be a robust and efficient global optimization method for the calibration of conceptual models. In this paper, two enhancements to the SCE‐UA algorithm are proposed, one to improve its exploration and another to improve its exploitation of the search space. A strategically located initial population is used to improve the exploration capability and a modification to the downhill simplex search method enhances its exploitation capability. This enhanced version of SCE‐UA is tested, first on a suite of test functions and then on a conceptual rainfall‐runoff model using synthetically generated runoff values. It is observed that the strategically located initial population drastically reduces the number of failures and the modified simplex search also leads to a significant reduction in the number of function evaluations to reach the global optimum, when compared with the original SCE‐UA. Thus, the two enhancements significantly improve the robustness and efficiency of the SCE‐UA model calibrating algorithm. Copyright © 2008 John Wiley & Sons, Ltd.     

Influence of plasticity on the seismic soil–micropiles–structure interaction

This paper includes an analysis of the influence of soil plasticity on the seismic response of micropiles. Analysis is carried out using a global three-dimensional modeling in the time domain. The soil behavior is described using the non-associated Mohr–Coulomb criterion. Both the micropiles and the superstructure are modeled as three-dimensional beam elements. Proper boundary conditions are used to ensure waves transmission through the lateral boundaries of the soil mass. Analyses are first conducted for harmonic loadings and then for real earthquake records. They show that plasticity could have a significant influence on the seismic response of the soil–micropiles–structure systems. This influence depends on the amplitude of the seismic loading and the dominant frequencies of both the input motion and the soil–piles–structure system.     

The influence of dimensionality on simulations of mass recovery from nonuniform dense non-aqueous phase liquid (DNAPL) source zones

The influence of model dimensionality on predictions of mass recovery from dense non-aqueous phase liquid (DNAPL) source zones in nonuniform permeability fields was investigated using a modified version of the modular three-dimensional transport simulator (MT3DMS). Thirty-two initial two- (2D) and three-dimensional (3D) tetrachloroethene–DNAPL source zone architectures, taken from a recent modeling study, were used as initial conditions for this analysis. Commonly employed source zone metrics were analyzed to determine differences between 2D and 3D predictions: (i) down-gradient flux-averaged contaminant concentration, (ii) reductions in contaminant mass flux through a down-gradient boundary, (iii) source zone ganglia-to-pool (GTP) ratio, and (iv) time required to achieve a remediation objective. 3D flux-averaged contaminant concentrations were approximately 3.5 times lower than concentrations simulated in 2D. This difference was attributed to dilution of the contaminant concentrations down gradient of the source zone. Contaminant flux reduction predictions for a given mass recovery were generally 5% higher in 3D simulations than in 2D simulations. The GTP ratio declined over time as mass was recovered in both 2D and 3D simulations. Although the source longevity (i.e., time required to achieve 99.99% mass recovery) differed between individual 2D and 3D realizations, the mean source longevity for the 2D and 3D simulation ensembles was within 2%. 2D simulations tended to over-predict the time required to achieve lower mass recovery levels (e.g. 50% mass recovery) due to a smaller contaminated area exposed to uncontaminated water. These findings suggest that ensemble averages of 2D numerical simulations of DNAPL migration, entrapment, dissolution, and mass recovery in statistically homogenous, nonuniform media may provide reasonable approximations to average behavior obtained using simulations conducted in fully three-dimensional domains.     

岩石爆破的相对损伤与损伤累积计算

为了便于工程应用,利用损伤变量定义,确定既真实反映力学作用过程,又易于获取和进行客观度量损伤变量的一种方法。通过提出的相对损伤概念,用相对损伤变量表示含损伤材料在载荷作用下损伤扩展程度和不同部位的相对损伤程度,使衡量受力作用程度的损伤变量确定更易实现。同时,提出了一种基于裂纹扩展的损伤累积计算方法,并用实例计算说明了其合理可行。     

松辽盆地徐家围子断陷芳深9井区火山岩储层裂缝预测

本文在松辽盆地北部徐家围子芳深9井区构造地质、石油地质及构造应力场分析的基础上,对早白垩世嫩江组末期、晚白垩世明水组末期和现今三个构造时期在营城组火山岩地层中的构造应力场进行了三维有限元数值模拟,进而分析了三维构造应力场在裂缝预测分析中的应用,研究结果与构造地质、石油地质分析结论相一致,对油气勘探开发具有指导意义。      

利用FFT技术计算大地水准面高若干问题研究

本文首先从普通二维 FFT算法的基本定义和要求出发 ,深入分析和比较了大地水准面二维平面和二维球面 FFT算法的特点和差异 ,找出了影响二维球面 FFT计算精度的主要误差源 ,并给出了二维球面 FFT计算公式的改进形式。在此基础上 ,本文详细分析讨论了参考场选取、积分球冠半径确定以及核函数改化对大地水准面计算结果的作用和影响。最后 ,本文利用全国 5′× 5′实测重力异常对我国海陆大地水准面进行了试算 ,并将计算结果分别同陆上GPS水准和海上卫星测高观测值作了比较。     

圆形建筑物倾斜观测的新方法

基于全站仪和光电测距仪,本文提出了圆形建筑物倾斜观测的新方法,并推导出了相应的计算公式。     

颅颌面外科与三维重建技术

目的 探讨电子束ＣＴ三维重建技术在颅颌面外科的适应症和应用价值。方法 采用美国Ｉｍａｔｒｏｎ公司的电子束ＣＴ（ｅｌｅｃｔｒｏｎ ｂｅａｍ ＣＴ,以下简称ＥＢＣＴ）Ｃ－１５０,对７６例严重颅颌面病人实行薄层ＣＴ容积扫描。将所获ＣＴ 象经数字接口传至加拿大ＩＳＧ公司生产的Ａｌｌｅｇｒｏ工作站进行三维重建。结果 ＥＢＣＴ成像技术能立体的、详尽和精确的显示机体组织三维解剖结构极其相互关系。其再现畸形或病体     

Nonlinear and 3D effects in leaky infragravity waves

In this paper, infragravity (IG) waves, forced by normally and obliquely incident wave groups, are studied using the quasi-3D (Q3D) nearshore circulation model SHORECIRC [Van Dongeren, A.R., I.A. Svendsen, 1997b. Quasi 3-D modeling of nearshore hydrodynamics. Research report CACR-97-04. Center for Applied Coastal Research, University of Delaware, Newark, 243 pp.], which includes the Q3D effects. The governing equations that form the basis of the model, as well as the numerical model and the boundary conditions, are described. The model is applied to the case of leaky IG waves. It is shown that the Q3D terms have a significant effect on the cross-shore variation of the surface elevation envelope, especially around the breakpoint and in the inner surf zone. The effect of wave groupiness on the temporal and spatial variation of all Q3D terms is shown after which their contribution to the momentum equations is analyzed. This reveals that only those Q3D coefficients, which appear in combination with the largest horizontal velocity shears make a significant contribution to the momentum equations. As a result of the calculation of the Q3D coefficients, the IG wave velocity profiles can be determined. This shows that in the surf zone, the velocity profiles exhibit a large curvature and time variation in the cross-shore direction, and a small — but essential — depth variation in the longshore direction.     

基于动力学震源模型的三维沉积盆地直下型断层地震动模拟

沉积盆地与近断层地震共同作用会增加地震破坏的风险水平,尤其是盆地下方直下型断层发震情况。采用动力学震源模型刻画断层破裂发震过程,开展沉积盆地直下型断层谱元法地震动模拟研究,探讨不同断层面初始剪应力和成核区位置下三维沉积盆地地表响应规律。研究结果表明,断层面应力降对盆地地表地震动的影响显著,在断层面强度一定的情况下,随着初始剪应力的增大,即应力降增大,盆地地表峰值响应增大,原因在于应力降的改变影响了断层破裂释放能量,进而引起断层破裂速度改变,最终导致盆地地表响应发生变化;改变断层面成核区位置会对盆地内部地震动分布规律产生影响,当成核区位置从断层中间向断层左侧移动时,盆地左侧地震动逐渐减小,而右侧地震动逐渐增大,最终表现为盆地右侧地震动显著高于盆地左侧,原因在于改变成核区位置后,导致近断层地震动的方向性效应发生变化。