Statistical inverse analysis based on genetic algorithm and principal component analysis: Applications to excavation problems and pressuremeter tests

This study concerns the identification of constitutive models from geotechnical measurements by inverse analysis. Soil parameters are identified from measured horizontal displacements of sheet pile walls and from a measured pressuremeter curve. An optimization method based on a genetic algorithm (GA) and a principal component analysis (PCA), developed and tested on synthetic data in a previous paper, is applied. These applications show that the conclusions deduced from synthetic problems can be extrapolated to real problems. The GA is a robust optimization method that is able to deal with the non‐uniqueness of the solution in identifying a set of solutions for a given uncertainty on the measurements. This set is then characterized by a PCA that gives a first‐order approximation of the solution as an ellipsoid. When the solution set is not too curved in the research space, this ellipsoid characterizes the soil properties considering the measured data and the tolerate margins for the response of the numerical model. Besides, optimizations from different measurements provide solution sets with a common area in the research space. This intersection gives a more relevant and accurate identification of parameters. Finally, we show that these identified parameters permit to reproduce geotechnical measurements not used in the identification process. Copyright © 2009 John Wiley & Sons, Ltd.     

Enhancement of backtracking search algorithm for identifying soil parameters

In this paper, an enhanced backtracking search algorithm (so-called MBSA-LS) for parameter identification is proposed with two modifications: (a) modifying the mutation of original backtracking search algorithm (BSA) considering the contribution of current best individual for accelerating convergence speed and (b) novelly incorporating an efficient differential evolution (DE) as local search for improving the quality of population. The proposed MBSA-LS is first validated with better performance than the original BSA and some other typical state-of-the-art optimization algorithms on a benchmark of soil parameter identification in terms of effectiveness, efficiency, and robustness. Then, the efficiency of the MBSA-LS is further illustrated by two representative cases: identifying soil parameters from both laboratory tests and field measurements. All comparisons demonstrate that the proposed MBSA-LS algorithm can give accurate results in a short time. Finally, to conveniently solve the problems of parameter identification, a practical tool ErosOpt for parameter identification is developed by integrating the proposed MBSA-LS and some other efficient algorithms for readers to conduct the parameter identification using optimisation algorithms.     

Single‐and multi‐objective genetic algorithm optimization for identifying soil parameters

This paper discusses the quality of the procedure employed in identifying soil parameters by inverse analysis. This procedure includes a FEM‐simulation for which two constitutive models—a linear elastic perfectly plastic Mohr–Coulomb model and a strain‐hardening elasto‐plastic model—are successively considered. Two kinds of optimization algorithms have been used: a deterministic simplex method and a stochastic genetic method. The soil data come from the results of two pressuremeter tests, complemented by triaxial and resonant column testing. First, the inverse analysis has been performed separately on each pressuremeter test. The genetic method presents the advantage of providing a collection of satisfactory solutions, among which a geotechnical engineer has to choose the optimal one based on his scientific background and/or additional analyses based on further experimental test results. This advantage is enhanced when all the constitutive parameters sensitive to the considered problem have to be identified without restrictions in the search space. Second, the experimental values of the two pressuremeter tests have been processed simultaneously, so that the inverse analysis becomes a multi‐objective optimization problem. The genetic method allows the user to choose the most suitable parameter set according to the Pareto frontier and to guarantee the coherence between the tests. The sets of optimized parameters obtained from inverse analyses are then used to calculate the response of a spread footing, which is part of a predictive benchmark. The numerical results with respect to both the constitutive models and the inverse analysis procedure are discussed. Copyright © 2011 John Wiley & Sons, Ltd.     

Investigation into MOGA for identifying parameters of a critical-state-based sand model and parameters correlation by factor analysis

Adding refinement and accuracy to constitutive models of soil results in the introduction of complexities along with more model parameters. These parameters (such as hardening-/softening-, dilatancy-/contractancy-related parameters and critical state parameters) are usually not easily obtained in a straightforward way. How to identify these key parameters and estimate their correlations of advanced soil models is a particular issue for geotechnical engineering. This paper was aimed to investigate multi-objective genetic algorithms for identifying parameters of advanced sand models based on standard laboratory tests, followed by the correlation analysis of parameters. A critical-state-based sand model has been developed to simulate three triaxial compression tests performed on loose and dense Hostun sand. Two widely used genetic algorithms with two initialisation methods are examined. The performance of the two genetic algorithms is assessed by comparing their simulation performance using optimal parameters, the convergence speed and the distribution of solutions on the Pareto front. The optimal parameters can then be classified into two factors by their correlation relationship.     

Statistical inverse analysis based on genetic algorithm and principal component analysis: Method and developments using synthetic data

This study concerns the identification of parameters of soil constitutive models from geotechnical measurements by inverse analysis. To deal with the non‐uniqueness of the solution, the inverse analysis is based on a genetic algorithm (GA) optimization process. For a given uncertainty on the measurements, the GA identifies a set of solutions. A statistical method based on a principal component analysis (PCA) is, then, proposed to evaluate the representativeness of this set. It is shown that this representativeness is controlled by the GA population size for which an optimal value can be defined. The PCA also gives a first‐order approximation of the solution set of the inverse problem as an ellipsoid. These developments are first made on a synthetic excavation problem and on a pressuremeter test. Some experimental applications are, then, studied in a companion paper, to show the reliability of the method. Copyright © 2009 John Wiley & Sons, Ltd.     

Selection of sand models and identification of parameters using an enhanced genetic algorithm

Numerous constitutive models of granular soils have been developed during the last few decades. As a consequence, how to select an appropriate model with the necessary features based on conventional tests and with an easy way of identifying parameters for geotechnical applications has become a major issue. This paper aims to discuss the selection of sand models and parameters identification by using genetic algorithm. A real‐coded genetic algorithm is enhanced for the optimization with high efficiency. Models with gradually varying features (elastic‐perfectly plastic modelling, nonlinear stress–strain hardening, critical state concept and two‐surface concept) are selected from numerous sand models as examples for optimization. Conventional triaxial tests on Hostun sand are selected as the objectives in the optimization. Four key points are then discussed in turn: (i) which features are necessary to be accounted for in constitutive modelling of sand; (ii) which type of tests (drained and/or undrained) should be selected for an optimal identification of parameters; (iii) what is the minimum number of tests that should be selected for parameter identification; and (iv) what is the suitable and least strain level of objective tests to obtain reliable and reasonable parameters. Finally, a useful guide, based on all comparisons, is provided at the end of the discussion. Copyright © 2015 John Wiley & Sons, Ltd.     

挡土墙库仑土压力的遗传算法求解分析

在对破裂面上滑动土体静力极限平衡分析的基础上,建立了基于优化方法求解无黏性土、黏性土库仑土压力的自变量取值区间和目标函数模型,并采用遗传进化方法进行了实例求解分析。研究结果表明,遗传算法在计算挡土墙库仑主动土压力的过程中,收敛速度快、用时短,并具有较高的计算精度。算例1中5组无黏性土挡土墙的主动土压力的计算结果与经典库仑解析解非常接近,平均误差为1.748 %,平均进化代数为15代。算例2中8组黏性土挡土墙的主动土压力计算结果与文献的解答非常吻合,平均误差仅为0.017 %,平均进化代数为17.125代。遗传算法具有良好的适应性和强大的搜索性能,非常适合求解岩土工程优化问题。     

Soil parameter identification using a genetic algorithm

This paper is dedicated to the identification of constitutive parameters of the Mohr–Coulomb constitutive model from in situ geotechnical measurements. A pressuremeter curve and the horizontal displacements of a sheet pile wall retaining an excavation are successively used as measurements. Two kinds of optimization algorithms are used to minimize the error function, the first one based on a gradient method and the second one based on a genetic algorithm. The efficiency of each algorithm related to the error function topology is discussed. Finally, it is shown that the use of a genetic algorithm to identify the soil parameters seems particularly suitable when the topology of the error function is complex. Copyright © 2007 John Wiley & Sons, Ltd.     

带有横向约束的全局优化波阻抗反演方法及应用

地震波阻抗反演是一种基于地震数据得到波阻抗参数的非线性优化问题。不依赖于目标函数梯度信息的全局优化算法是求解地震波阻抗反演问题的有效方法。然而这类方法采用逐道反演的策略,忽略了相邻地震道的空间相关性,导致反演结果的横向连续性差。针对该问题,提出了一种融入旁道最优解的模型空间初始化方法约束波阻抗反演的搜索空间范围,以改善反演结果的横向连续性,进而提出了一种带有横向约束的多组变异差分进化地震波阻抗反演方法。合成地震记录算例表明,该方法相比于传统方法不仅具有更快的收敛速度,而且反演结果具有更好的横向连续性。此外,该方法被应用于胜利油田某区块的储层波阻抗参数反演,反演结果与测井资料有很好的一致性,且有效地刻画了储层砂岩厚度。     

基于INPGA的地下水污染治理多目标优化管理模型:Ⅰ——理论方法与算例验证

基于小生境技术的Pareto遗传算法(NPGA)是一种求解多目标问题的智能搜索方法,适用于优化多种非线性、不连续等复杂多目标问题.但该算法存在局部早熟收敛和收敛速度慢两个不足,在求解Pareto前沿上效果不佳.本文在NPGA的基础上,提出了改进NPGA方法(INPGA),通过Pareto解集过滤器、精英个体保留策略、邻...     

Evaluating the performance of an explicit dynamic relaxation technique in analyzing non-linear geotechnical engineering problems

Explicit dynamic relaxation is an efficient tool that has been used to solve problems involving highly non-linear differential equations. The key feature of this method is the ability to use explicit dynamic algorithms in solving static problems. Few attempts have been made to date to apply this technique in conventional geotechnical engineering. In this study, an algorithm that incorporates the application of a stiffness dependent time step scheme is proposed. The algorithm has been successfully used to solve 2D and 3D non-linear geotechnical engineering problems. To calibrate the developed algorithm, numerical simulations have been conducted for a strip and square footings supported by Mohr–Coulomb material. Performance of four different types of brick elements used in collapse load calculation is examined in terms of convergence speed and accuracy. In addition, the role of employing adaptive time steps in reducing the number of iterations needed for convergence is also evaluated.     

基于改进麻雀搜索算法的瑞利波频散曲线反演

非线性优化算法在给定的参数搜索范围内对最优解进行全局搜索,在全局搜索方面具有先天的优势,具有一定的跳出局部极值的能力。本文将一种新兴的非线性优化算法——麻雀搜索算法引入瑞利波频散曲线反演问题,针对频散曲线反演问题瑞利波频散曲线反演问题多参数、多局部极值的特点,引入自适应t分布对算法进行改进。三种理论模型的反演实验数据表明,改进的麻雀搜索算法与传统麻雀搜索算法相比具有更好的反演精度和稳定性,同时具有较好的抗随机噪声的能力。与粒子群算法和差分进化算法两种较成熟的非线性优化算法进行对比,改进的麻雀搜索算法较好地平衡了迭代前期的全局搜索和迭代后期的局部搜索,取得了与粒子群算法和差分进化算法相比更好的效果。     

SBPT测定饱和黏土不排水强度的数值分析

自钻式旁压试验（SBPT）因其扰动小、测试深度大、可以获得应力-应变、超孔隙水压力-时间等数据,在确定地基土性参数和地基承载力上有广阔的应用前景。然而由于用以解释SBPT的柱孔扩张理论（Gibson解）所采用的平面应变假设与实际旁压腔几何特征存在差异,导致试验所确定的黏土不排水剪切强度su与其他原位试验或室内试验结果存在差别。针对旁压腔几何尺寸及应变区间的选择对确定su的影响,基于修正剑桥模型,采用低渗透系数控制加载过程中不排水条件,利用有限元法模拟SBPT,建议了不同应力历史下确定su的应变区间,并给出考虑几何尺寸影响时相应应变区间上su的修正系数。     

Identifying parameters controlling soil delayed behaviour from laboratory and in situ pressuremeter testing

The aim of this paper is to present a methodology for identifying the soil parameters controlling the delayed behaviour from laboratory and in situ pressuremeter tests by using an elasto‐viscoplastic model (EVP‐MCC) based on Perzyna's overstress theory and on the elasto‐plastic Modified Cam Clay model. The influence of both the model parameters and the soil permeability was studied under the loading condition of pressuremeter tests by coupling the proposed model equations with Biot's consolidation theory. On the basis of the parametric study, a methodology for identifying model parameters and soil permeability by inverse analysis from three levels of constant strain rate pressuremeter tests was then proposed and applied on tests performed on natural Saint‐Herblain clay. The methodology was validated by comparing the optimized values of soil parameters and the values of the same parameters obtained from laboratory test results, and also by using the identified parameters to simulate other tests on the same samples. The analysis of the drainage condition and the strain rate effect during a pressuremeter test demonstrated the coupled influence of consolidation and viscous effects on the test results. The numerical results also showed that the inverse analysis procedure could successfully determine the parameters controlling the time‐dependent soil behaviour. Copyright © 2008 John Wiley & Sons, Ltd.     

PSO-RBFNN模型及其在岩土工程非线性时间序列预测中的应用

岩土工程受力变形演化是一个典型的非线性问题,其演化的高度非线性和复杂性,很难用简单的力学、数学模型描述,但可用粒子群优化径向基神经网络对岩土工程应力、位移非线性时间序列进行动态实时预测。网络径向基层的单元数通过均值聚类法确定后,所有其它参数：中心位置、形状参数、网络权值,均通过粒子群优化算法在全局空间优化确定。工程实例应用表明,该模型预测结果准确、精度高,有良好的应用前景。     

水文地质参数反演的Hooke-Jeeves粒子群混合算法

水文地质参数寻优结果的好坏会直接影响到地下水数值模拟的精度,而参数寻优结果很大程度上取决于寻优方法的选择。粒子群算法是一种基于群智能的随机全局寻优方法,算法的缺陷是后期搜索效率低劣。基于随机寻优算法的混合策略,引入有效的约束处理手段和粒子群算法惯性因子的动态非线性调整技术,有机融合粒子群算法与Hooke-Jeeves方法,提出一种适用于水文地质参数反演的HJPSO混合算法。应用研究表明,HJPSO混合算法在参数反演计算中求解精度高、收敛速度快、寻优性能强,是一种值得推广的水文地质参数识别方法。     

Simulation-based calibration of geotechnical parameters using parallel hybrid moving boundary particle swarm optimization

Simulation-based optimization methods have been recently proposed for calibrating geotechnical models from laboratory and field tests. In these methods, geotechnical parameters are identified by matching model predictions to experimental data, i.e. by minimizing an objective function that measures the difference between the two. Expensive computational models, such as finite difference or finite element models are often required to simulate laboratory or field geotechnical tests. In such cases, simulation-based optimization might prove demanding since every evaluation of the objective function requires a new model simulation until the optimum set of parameter values is achieved. This paper introduces a novel simulation-based “hybrid moving boundary particle swarm optimization” (hmPSO) algorithm that enables calibration of geotechnical models from laboratory or field data. The hmPSO has proven effective in searching for model parameter values and, unlike other optimization methods, does not require information about the gradient of the objective function. Serial and parallel implementations of hmPSO have been validated in this work against a number of benchmarks, including numerical tests, and a challenging geotechnical problem consisting of the calibration of a water infiltration model for unsaturated soils. The latter application demonstrates the potential of hmPSO for interpreting laboratory and field tests as well as a tool for general back-analysis of geotechnical case studies.     

多维连续参数的随机试验蚁群算法优化反演分析

针对一般蚁群算法难以求解优化反演中的多维连续参数优化问题,借鉴进化思想,提出随机试验蚁群算法求解多维连续参数全局优化问题的方法。该方法主要运用随机试验求出每个参数的任一水平对目标函数的影响度,以此影响度为蚂蚁选择下一个参数水平的能见度,从而确定参数水平被选中的概率。同时考虑到当参数被划分的水平很多时目标函数收敛较慢的缺点,提出用方差分析法来确定参数的水平对目标函数的影响是否显著,将对目标函数影响不显著的参数水平其信息素设置为很小的值,从而加速目标函数的收敛。通过对一个岩土工程优化实例,表明了该算法的有效性。     

Application of Seagull Optimization Algorithm Log Interpretation to Shale Gas Reservoir of Well H in Sichuan Basin Yuxi Block

Due to the diversity of mineral types in shale gas reservoirs, it is difficult to establish reservoir parameter volume model by conventional log interpretation methods. The optimization log interpretation method can evaluate complex lithology reservoirs effectively, and the key is optimization algorithm. With the newly proposed seagull optimization algorithm method, we calculate the mineral and physical parameters of shale gas reservoir in Well H of Yuxi block, Sichuan Basin, and compare with the genetic algorithm and the genetic algorithm-complex hybrid algorithm. It shows that calculation results of seagull optimization algorithm optimization log interpretation match well with core analysis data, and calculation error is small, calculation speed is fast. Seagull optimization algorithm also makes up for the shortcomings of premature convergence and easy to fall into local optimization of genetic algorithm, the need for secondary optimization and slow search speed of genetic-complex hybrid algorithm. It provides a reference for the application of seagull optimization algorithm in other shale gas reservoirs regions.     

利用混沌粒子群算法确定河流水质模型参数

将混沌寻优思想引入到粒子群优化算法中,提出了混沌粒子群算法,这种方法利用混沌运动的随机性、遍历性和规律性等特性对当前粒子群体中的粒子进行混沌寻优。通过这种处理使得粒子群体的进化速度加快,从而改善了粒子群优化算法摆脱局部极值点的能力,提高了算法的收敛速度和精度。并将混沌粒子群算法应用于求解分析瞬时投放示踪剂情况下的一维河流水团示踪试验数据以及确定河流水质参数的函数优化问题,结果表明,混沌粒子群算法的收敛性能明显优于粒子群优化算法。