基于QSSP的震源破裂过程反演方法

在Qt平台上,基于QSSP软件,利用C⁺⁺进行地震破裂过程反演方法研究。采用按阶段反演的方法,第1阶段反演采用热浴算法进行全局搜索,先按断层区域的构造背景等初始条件对参数范围进行划取,在该范围内随机给定一组初始值开始迭代,使波形初步拟合避开局部最优解;第2阶段反演使用拟牛顿法进行快速收敛,提高波形拟合程度,迭代至目标函数小于误差条件时停止,输出满足误差条件的待解模型参数。为避免模型参数出现病态问题,使用拉普拉斯方程建立平滑矩阵并引入平滑因子对断层模型进行平滑约束。使用棋盘模型验证该方法的稳定性和可靠性。最后,将全国13个台站的重力数据积分后对2013-04-20芦山7.0级地震的破裂过程进行反演,并与其他研究结论进行对比分析。     

An iterative pressure‐stabilized fractional step algorithm in saturated soil dynamics

This paper presents an iterative, incremental pressure‐stabilized fractional step algorithm for coupled hydro‐mechanical problems with mixed formulations of the displacement–pressure ( u –p) model in saturated soil dynamics that allows the use of finite elements with equal low order of interpolation approximation of u and p. In comparison with the original fractional step algorithm, the distinct features of the proposed algorithm lie in its enhanced stability owing to the introduction of both an iteration procedure and a finite increment calculus (FIC) process into the algorithm. The introduction of the iterative procedure makes the velocity term satisfy the momentum conservation equation in an implicit sense and allows much larger time step sizes to be used than those limited in existing explicit and semi‐implicit versions of the algorithm. The introduction of the FIC process removes the dependence of the stability of the proposed algorithm on the time step size, as a result it allows to using the incremental version of the algorithm and evades the minimum time step size requirement presented in the existing versions of the fractional step algorithm that restricts the application of the algorithm to saturated soil dynamics problems with high frequencies. Numerical experiments demonstrate the effectiveness and improved performance of the proposed iterative pressure‐stabilized fractional step algorithm. Copyright © 2009 John Wiley & Sons, Ltd.     

Multi-objective autocalibration of SWAT model for improved low flow performance for a small snowfed catchment

ABSTRACT

A reliable modelling framework needs to ensure that the model is simulating reality with limited uncertainty, thus enhancing its predictive ability. In the literature, hydrological model assessment using one or more metrics is reported to be inadequate when the river flow regime is required to be reproduced comprehensively. This research is aimed to: (a) calibrate the Soil and Water Assessment Tool (SWAT) based on the concept of multi-objective optimization by applying the Borg multi-objective evolutionary algorithm (MOEA); (b) apply hydrological signatures as objective functions; and (c) adopt a multi-metric approach for model evaluation. The SWAT model was coupled with a relatively newer and powerful Borg MOEA. The inclusion of hydrological signatures as objective functions along with the conventional statistical functions assisted in improving the performance for low flows by 135% in terms of volume efficiency and 65% for flow time series simulation.     

Estimation of Nash's IUH parameters using stochastic search algorithms

In rainfall–runoff studies, it is often necessary to change the duration of a given unit hydrograph. Nash's Instantaneous Unit Hydrograph (IUH) is an ideal method that eliminates the hydrograph duration. This paper presents the results of the application of search algorithms, namely a genetic algorithm and hill climbing, to develop the IUH that minimizes the error between the observed and generated hydrographs. Also the performance of these methods has been compared with that of the classical method used for estimation of IUH, namely the method of moments. The genetic algorithm is a popular search procedure for function optimization that applies the mechanics of natural genetics and natural selection to explore a given search space. Hill climbing is an optimization technique that belongs to the family of local search and algorithms can be used to solve problems that have many solutions, with some solutions better than others. The results obtained from both the genetic algorithm and hill climbing algorithm for estimation of Nash's IUH parameters were compared with the results obtained by the method of moments for storms from two river basins that are located in different climatic regions. It was found that both the genetic algorithm and hill climbing provided improved and consistent results. Copyright © 2008 John Wiley & Sons, Ltd.     

用微机进行Mie散射算法研究

本文介绍在东海微机上用伦茨算法计算Mie级数的系数,其计算程序已运行成功。计算结果表明此法可使误差传播显著减小,计算精度提高,降低了对微机内存客量的要求,同时,本文得到的计算结果与国外公开发表的计算结果达到同样精度。     

A modified land surface temperature split window retrieval algorithm and its applications over China

Due to the difficulties in correcting the influences of the atmosphere absorbability and the Earth surface emissivity diversification, the retrieval of LST (land surface temperature) from satellite data is a challenging task. In this paper, a modified Becker's split window LST inversion algorithm is developed for retrieving LST from the NOAA-16/17 AVHRR data. A new set of parameters for Becker's LST algorithm is proposed. The algorithm is developed from a surface brightness temperature dataset generated from the MODTRAN program, which uses a range of surface parameters and atmospheric quantities as inputs. The 10-day composites of the channels 4 and 5 brightness temperature data of NOAA-17 AVHRR (1-km resolution) are used to generate the clear-sky LST. As a validation of the algorithm, the retrieved LST is compared with MODIS LST of same period and area. The two LST products are found to be consistent, with the absolute difference being about 2.5 K for most areas. The NOAA retrieved LST is also compared with in-situ ground surface 0-cm measurements taken from 257 meteorological stations, which cover overall China area for the three periods of satellite observations. The comparison shows that the correlation between the retrieved LST and in-situ measurements is over 0.90 and the RMSE (root mean square error) is about 3.4 K.     

一类非线性投影方程解的扰动迭代算法及其收敛性

在Hilbert空间中讨论了一类非线性投影方程解的扰动迭代算法及其收敛性分析．     

基于替换的组合电路的等价性检验方法

介绍了基于替换的组合电路的等价性检验算法，利用待检验的两个电路的结构相似性来逐步约简电路，从而加速了验证过程。     

球面等积网格系统生成算法与软件模型研究

利用全球离散网格系统管理海量空间信息是一种受到普遍关注的新途径,具有等积特性的网格系统在大范围遥感影像数据处理方面具有重要的理论意义和应用前景。文中论述的算法在二十面体展开图上建立坐标系描述不同剖分产生的多分辨率网格,利用施奈德多面体投影映射到球面。针对算法特点设计具有网络模拟功能的软件模型。实验结果表明,该软件模型具有较好的扩展性。     

SISA完备性监测与WUL可靠性算法

详述了Galileo完备性概念,分析了空间信号有效参数完备性监测系统和监测机理,给出Galileo系统中描述SISA预报形式和算法模型.基于空间信号完备性参数的算法,优化计算最坏用户位置方法,建立了算法模型,得出相应结论.算法和结论初步验证适用于全球,为进一步研究Galileo系统完备性算法和估计用户保护水平奠定基础.