基于优化参数的地震P、S波震相到时自动拾取及质量评估

准确拾取P、S波震相到时是深入开展地震波研究工作的基础,本文改进了自动拾取参数优化函数算法和质量评估方案,引入了拾取到时优化方案,使用基于参数优化的频带-带宽拾取算法、AICD拾取算法和峰度拾取算法对腾冲地区7个宽频带地震台站记录的地震资料开展了地震P、S波到时自动拾取,对拾取结果进行了优化和质量判定.结果表明：经参数优化、拾取优化后,采用3种方法自动拾取的P、S波到时与人工拾取到时的时差在0.1 s内的记录占比分别达到74.66%、70.98%.这些参数值均优于算法改进前的同类参数,证明了优化方法的可靠性.     

基于改进粒子群算法的地震标量波方程反演

针对标准粒子群优化(PSO)算法存在易出现早熟而陷入局部最优以及进化后期收敛速度慢等缺陷,通过考虑粒子所处位置间相互作用,提出了一种改进的并行粒子群优化算法.由于引入粒子位置间的相互影响,减少了粒子搜索过程盲目性,因此能有效提高算法的收敛速度.数值试验表明,这种改进的粒子群算法适用于二维标量波方程的速度反演,且算法具有...     

基于改进差分进化算法的大地电磁测深和重力宽范围物性约束联合反演

宽范围物性约束技术容易实现、具有一定容错性,目前已在大地电磁测深（MT）和地震、MT和重力联合反演中实现,但该技术是结合模拟退火算法实现的.差分进化算法（DE）是一种全局优化算法,但该算法在地球物理联合反演领域应用较少.基于此,本文以双种群设置方案为框架改进了DE算法,并提出了基于改进DE算法的宽范围物性约束技术.MT和重力联合反演的模型试验表明：与传统的DE算法相比,改进的DE算法收敛速度更快,寻优能力更强;基于改进DE算法的宽范围物性约束技术可以促进不同岩石物性参数在一定"范围"内实现耦合,既可以利用岩石物性关联的导向作用,又可以发挥优化算法的寻优能力,进而降低地球物理联合反演对先验信息的要求;此外,该技术的实现也验证了宽范围物性约束思想在联合反演领域中的适用性,具有进一步推广至其他优化算法中的潜质.     

最大准则优化技术在贴体网格中的应用

贴体网格在地质数值模拟中具有广阔的应用前景,为解决贴体网格生成时边界离散问题,提出了最大长度准则和最大面积准则,把曲线逼近和曲面网格优化问题转化为数学优化问题,为求解该问题,提出了改进的单粒子优化算法.试验表明,最大长度准则和最大面积准则的优化效果好于常规方法;以改进的单粒子优化算法求解该问题时,计算效率是智能单粒子优化算法的30倍左右（节点量为200）,从而实现最大长度准则和最大面积准则在贴体网格生成中的应用.针对最大面积准则优化曲面网格不能控制网格步长的情况,提出了限定步长的网格优化算法,使网格步长合理化,并通过实例验证了该算法的有效性.研究成果提供了生成贴体网格时边界优化准则和求解方法,对今后复杂边界的贴体网格生成具有重要意义.     

基于子空间的二维大地电磁量子遗传反演法研究

量子遗传算法作为一种高效的优化算法,仍存在容易陷入局部极值的缺点.为提高算法的高效性,并探讨将算法应用于大地电磁二维反演的可行性和有效性,本文对算法进行了改进,并通过一维两层D型和四层HK型模型数值试验验证了改进的有效性.然后将改进后的算法引入二维大地电磁反演,在引入滑动子空间思想,同时只考虑最简化反演条件的前提下,对...     

面向目标自适应三维大地电磁正演模拟

本文将面向目标的自适应算法应用于三维大地电磁数值模拟.使用基于非结构网格的矢量有限单元法对起伏地表大地电磁正演模拟问题进行求解.使用利用垂向电流密度在物性界面上的连续性对后验误差进行估算的算法指导网格优化.由于全局自适应算法针对观测点优化网格的能力较差,本文通过求解正演问题的对偶问题计算后验误差的加权系数,并对相关加权系数进行改进,从而实现了面向目标的自适应算法.与传统基于结构化网格的电磁正演算法相比,采用非结构网格能够更好地拟合起伏地表和地下不规则异常体.由于使用了面向目标的自适应算法,本文能够使用更少的网格达到较高的计算精度.通过对比本文模拟结果与半空间响应和全局自适应算法计算结果,并通过对比使用改进前和改进后加权系数得到的网格剖分结果验证了本文算法的有效性.     

基于模拟植物生长算法的桁架结构截面与形状优化分析

针对桁架结构优化问题，给出基于改进的生长空间限定与并行搜索的模拟植物生长算法（GSL&PS-PGSA）优化模型，对典型平面桁架、空间桁架结构分别进行截面优化及形状优化分析，结果表明：改进GSL&PS-PGSA算法的优化效果显著，其优化所得结构经济高效；对于平面桁架截面优化，相比二级优化法、相对差商法及连续变量法，改进GSL&PS-PGSA所得结构较其减少了8%～15%的用钢量；对于空间桁架优化，截面优化所得结构总质量较序列两级算法轻约8.1%，较蚁群算法轻约3.6%，而形状优化所得结构总质量则分别较序列两级算法、蚁群算法轻约31.2%、22.4%。因此，改进的GSL&PS-PGSA算法对桁架结构优化问题具有良好的适用性。     

大地电磁阻尼粒子群优化反演法研究

粒子群优化算法（PSO）是模仿鸟群寻找食物的社会行为的一种全局最优化算法,在多维空间函数寻优、动态目标寻优等方面有着收敛速度快、解质量高且需要设置的参数较少等优点.本文在研究常规粒子群优化算法的基础上,对常规的粒子群算法进行了改进,提出了一种新的惯性权重ω参数振荡递减策略,加快了PSO算法的收敛速度,构造的新算法称为阻尼粒子群优化算法.在MATLAB 6.5 编程环境中对阻尼PSO算法进行了数值实验,并对大地电磁测深的理论模型和实测数据进行了反演试算,结果表明,阻尼PSO算法不依赖于初始模型、能够搜索到全局极值,不易陷入局部极值,是一种快速有效的地球物理反演方法.     

改进的模拟退火算法及其在叠前储层参数反演中的应用

快速模拟退火算法中,要求扰动后新的模型参数必须在取值范围内,对此,人们通常采用两种方法:(1)反复扰动直到新的模型参数满足要求;(2)当新的模型参数超过最大值时,参数等于最大值;新的模型参数小于最小值时,参数等于最小值.研究发现,方法1在计算过程中增加了寻找新模型参数的时间开销,计算效率明显降低;方法2当目标函数在边界或边界的小邻域存在多个极小值时,在有限的时间内,其陷入局部极小的概率大大增加,优化效果明显下降.本文在前人研究的基础上提出了改进的模拟退火算法,同时利用改进的算法反演了地层泊松比、纵波速度和密度.结果表明,改进的算法加快了模拟退火算法的收敛速度,提高了模拟退火算法的效率和精度.     

基于改进余弦组合窗的解耦方程弹性波逆时偏移

弹性波矢量波场逆时偏移可以综合利用纵横波场信息,对地下空间进行清晰成像,且对于成像介质没有角度限制,可以对复杂构造进行更清晰的成像.而弹性波逆时偏移中最重要的就是求解波动方程的算法,其直接影响成像的精度以及效率.本文引进电力系统谐波分析中常用的余弦组合窗函数,并通过一种新的优化算法得到了改进的余弦组合窗函数从而得到优化后的有限差分算子.并将此算子应用于解耦方程的矢量波场分离算法从而提高了成像精度.数值测试表明基于新算法的逆时偏移的成像精度和清晰度得到了明显的提高.     

A fuzzy logic based dynamic wave model inversion algorithm for canal regulation

A fuzzy logic based centralized control algorithm for irrigation canals is presented. Purpose of the algorithm is to control downstream discharge and water level of pools in the canal, by adjusting discharge release from the upstream end and gates settings. The algorithm is based on the dynamic wave model (Saint‐Venant equations) inversion in space, wherein the momentum equation is replaced by a fuzzy rule based model, while retaining the continuity equation in its complete form. The fuzzy rule based model is developed on fuzzification of a new mathematical model for wave velocity, the derivational details of which are given. The advantages of the fuzzy control algorithm, over other conventional control algorithms, are described. It is transparent and intuitive, and no linearizations of the governing equations are involved. Tuning of the algorithm and method of computation are explained. It is shown that the tuning is easy and the computations are straightforward. The algorithm provides stable, realistic and robust outputs. The disadvantage of the algorithm is reduced precision in its outputs due to the approximation inherent in the fuzzy logic. Feed back control logic is adopted to eliminate error caused by the system disturbances as well as error caused by the reduced precision in the outputs. The algorithm is tested by applying it to water level control problem in a fictitious canal with a single pool and also in a real canal with a series of pools. It is found that results obtained from the algorithm are comparable to those obtained from conventional control algorithms. Copyright © 2009 John Wiley & Sons, Ltd.     

Damage detection of base‐isolated buildings using multi‐input multi‐output subspace identification

A damage detection algorithm of structural health monitoring systems for base‐isolated buildings is proposed. The algorithm consists of the multiple‐input multiple‐output subspace identification method and the complex modal analysis. The algorithm is applicable to linear and non‐linear systems. The story stiffness and damping as damage indices of a shear structure are identified by the algorithm. The algorithm is further tuned for base‐isolated buildings considering their unique dynamic characteristics by simplifying the systems to single‐degree‐of‐freedom systems. The isolation layer and the superstructure of a base‐isolated building are treated as separate substructures as they are distinctly different in their dynamic properties. The effectiveness of the algorithm is evaluated through the numerical analysis and experiment. Finally, the algorithm is applied to the existing 7‐story base‐isolated building that is equipped with an Internet‐based monitoring system. Copyright © 2004 John Wiley & Sons, Ltd.     

A split operator approach to reactive transport with the forward particle tracking Eulerian Lagrangian localized adjoint method

A forward particle tracking Eulerian Lagrangian localized adjoint method (ELLAM) is applied to the multicomponent reactive transport problem using a split operator approach. Two split operator algorithms are compared, the Strang algorithm and the sequential non-iterative algorithm (SNIA). The reaction equations are integrated using a coupled predictor corrector algorithm with adaptive time stepping. Reaction time steps are adjusted at the inflow boundary to reflect the actual time of transport inside the solution domain.Results show that split operator ELLAM formulations are competitive with direct or fully coupled ELLAM solutions for reactive transport problems. The SNIA algorithm is more accurate than the Strang splitting algorithm when large time steps are used. The reaction algorithm employed dominates computational effort in runs with large time step sizes. To illustrate the use of the method in practical problems, the model is fitted to aerobic aniline degradation data from laboratory scale column experiments. Model inversion is achieved using non-linear regression with a shuffled complex evolution optimization algorithm and parameter uncertainty is assessed using a Bayesian uncertainty analysis procedure.     

基于BS-MUSIC的定向钻孔雷达三维成像算法

定向钻孔雷达由位于井中的一个发射天线和分布于圆环上的呈等角度分布的四个接收天线组成,能确定井旁地质目标的深度和方位.然而,目标方位确定和三维成像的算法问题一直是该技术面临的挑战.本文将空间谱估计中的BS-MUSIC算法应用于定向钻孔雷达目标方位识别算法中,在此基础上提出并实现了利用四个接收天线的信号对井周地质体进行三维空间成像的算法.使用合成数据对以上算法的效果进行验证,并和通常的MUSIC算法进行比较,结果表明,MUSIC算法和基于波束空间的MUSIC算法均可以在方位确定和三维成像方面取得非常好的效果,但计算速度方面,BS-MUSIC算法相比于MUSIC算法有了非常明显的提升.     

CT 成像的改进算法及其实验

Ｆｏｕｒｉｅｒ算法是求解Ｒａｄｏｎ变换的一种方法,优点是计算量小,速度快,但因其成像质量较差,所以在医疗ＣＴ方面的应用较少。本文叙述了一种改进的Ｆｏｕｒｉｅｒ算法,并对其作了成像实验,实验结果表明改进Ｆｏｕｒｉｅｒ 算法的成像质量不低于国产ＣＴ所采用的卷积反投影算法,但在速度上却有很大优势。     

BFA-CM最优化测井解释方法

最优化测井解释方法能充分利用各种测井资料及地质信息,可以有效地评价复杂岩性油气藏.优化算法的选择是最优化测井解释方法的关键,影响着测井解释结果的准确性.细菌觅食算法(BFA)是新兴的一种智能优化算法,具有较强的全局搜索能力,但在寻优后期收敛速度较慢.复合形算法(CM)局部搜索能力极强,将其与BFA算法相结合构成BFA-CM混合算法,既提高了搜索精度又提高了搜索效率.利用BFA-CM最优化测井解释方法对苏里格致密砂岩储层实际资料进行了处理,计算结果与岩心及薄片分析资料吻合得很好.     

基于模拟退火算法的设计反应谱标定方法

反应谱的标定是抗震设计的基础工作之一,模拟退火算法是基于模拟固体退火过程而提出的多参数优化组合方法。本文提出将模拟退火算法应用于设计反应谱的标定,并概述了模拟退火算法的基本原理及特点,介绍了运用MATLAB基于模拟退火算法对反应谱进行标定的过程。检验了将模拟退火算法应用于场地相关谱标定的可行性和合理性,给出了将其应用于工程场地地震安全性评价中设计反应谱标定实例。通过检验和实例分析可以看出,基于模拟退火算法的反应谱标定方法所给出的设计反应谱谱形真实地反映了原地震反应谱的特征,较客观的反映了场地相关反应谱的峰值和周期特征。     

PROGRESS IN THE DIRECT INTERPRETATION OF RESISTIVITY SOUNDINGS: AN ALGORITHM*

An algorithm is presented for the direct interpretation of resistivity sounding data. The algorithm is based on the method of successive reductions to lower boundary plane of the resistivity transform function. A novel aspect of the algorithm is that error limits are assigned to the initial values of the resistivity transform, and these error limits are carried through in all the subsequent computations. The width of the error range is then used as the basis for assigning weight factors in the final computation of thicknesses and resistivities of the layers. The errors in the resistivity transform derived from the solution given by the algorithm are usually not more than twice as large as those in the original data.     

Hybrid genetic algorithms in view of the evolution theories with application for the electrical sounding method

The classical genetic algorithm is a stochastic process which operates by natural selection. Although the algorithm may localize a point around the global minimum of the misfit function, it is not efficient at finding the precise solution. This paper suggests some hybrid genetic algorithms, derived from evolution theories, to overcome this problem. Firstly, sexual selection has been incorporated in the classical genetic algorithm to obtain a full representation of the Darwinist evolution concept. The simulation of sexual selection is performed by assigning a higher probability of surviving to some parameters that satisfy some algebraic relationships. This method is called the ‘marked constraints’ algorithm since it permits us to insert geological and geophysical constraints into the problem. The algorithm implementation is realized by progressively shrinking the parameter search space through successive generations. In this way, the genetic algorithm gains some degree of determinism. Secondly, since the evolution theory of Lamarck postulates that the acquired traits are passed on to the next generation, a hybrid use of the damped least‐squares method and the genetic algorithm is called Lamarckian inversion. Lamarckian inversion involves some improvement procedures that simulate the reduction of the misfit with the help of a derivative‐based method between two generations. Finally, although there is no correspondence in nature, Lamarckian and Darwinist evolution concepts are combined to strengthen the deterministic part of the solution algorithm. This is called the Lamarckian‐marked‐constraint algorithm. The merits and behaviours of the suggested algorithms are discussed using two examples. The first is a hypothetical example affected by a multiminima problem. The second examines the equivalence problem using vertical electrical sounding data.     

Partial eigensolution of damped structural systems by Arnoldi's method

An efficient numerical algorithm is developed to solve the quadratic eigenvalue problems arising in the dynamic analysis of damped structural systems. The algorithm can even be applied to structural systems with non-symmetric matrices. The algorithm is based on the use of Arnoldi's method to generate a Krylov subspace of trial vectors, which is then used to reduce a large eigenvalue problem to a much smaller one. The reduced eigenvalue problem is solved and the solutions are used to construct approximate solutions to the original large system. In the process, the algorithm takes full advantage of the sparseness and symmetry of the system matrices and requires no complex arithmetic, therefore, making it very economical for use in solving large problems. The numerical results from test examples are presented to demonstrate that a large fraction of the approximate solutions calculated are very accurate, indicating that the algorithm is highly effective for extracting a number of vibration modes for a large dynamic system, whether it is lightly or heavily damped.