基于WA、ANN和水文频率分析法相结合的中长期水文预报模型的研究

针对传统中长期水文预报方法模拟预测结果精度低、未考虑水文不确定性因素的影响等问题,本文将小波分析(WA),人工神经网络(ANN)和水文频率分析法联合使用,建立了不确定性中长期水文预报模型:即在应用WA揭示水文序列变化特性的基础上,将原序列分为主序列和随机序列两部分,然后利用ANN对主序列进行模拟预测,对随机序列进行水文频率分析,最后将两部分结果叠加作为最终预测值.将该模型用于黄河河口地区作中长期水文预报,并与传统方法作对比,进行模型验证.结果显示:该模型能同时揭示序列的时、频结构和变化特性;预报值结果精度高;且合格率高;能定量分析和描述水文不确定性因素对预报结果的影响,可得到不同频率对应水文序列的模拟预测值.因此该模型的预报结果更加合理有效,对实际生产应用更具有指导意义.     

新安江模型参数多目标优化研究

水文模型的参数优化率定一直以来是水文预报领域的重要研究内容,当水文模型的结构确定后,水文模型参数的选择对水文模型整体性能和水文预报结果的好坏有着至关重要的影响.针对传统水文模型参数优选采用单一目标不能充分全面挖掘水文观测资料中蕴含的水文特征信息的缺陷,本文以新安江三水源模型为例,尝试采用多目标优化算法优化率定水文模型,算例应用分析表明,通过合理的选择目标函数的种类和数目,采用多目标进化算法优化率定模型参数,可以获得相对于单目标率定模型参数更优的结果.进一步,研究工作针对模型参数优化的结果进行分析,可以明显看出模型参数优化中存在“异参同效”现象,为后续模型参数不确定性分析等相关研究工作的开展做好了铺垫.     

基于分布式水文模拟的三峡区间洪水预报(Ⅰ)——模型构建及验证

长江三峡区间因暴雨形成的洪水峰高量大,对三峡水库的防洪安全和运行调度的影响很大.本论文依据三峡地区的地形地貌特征,采用基于GIS的机理性分布式水文模型,来模拟三峡区间入库洪水,以尽量减少洪水预报中的不确定性.利用近期建成的78个自动雨量站网监测的小时降雨信息作为模型的输入,对模型参数进行了率定和验证,结果表明:大多数洪水过程的模拟精度较好,但也有的模拟结果较差,其中降雨信息缺失是洪水预报不确定性的主要来源.     

贝叶斯概率水文预报系统在中长期径流预报中的应用

采用贝叶斯概率水文预报理论制订水电站水库中长期径流预报模型,以概率分布的形式定量地描述水文预报的不确定度,探索概率水文预报理论及其应用价值。采用气象因子灰关联预报模型处理输入因子的不确定度,将实时气象信息和历史水文资料有效结合,突破传统确定性预报方法在信息利用和样本学习方面的局限性,以提高水文预报的精确度。以丰满水电厂水库为例对所建模型进行检验,模拟计算结果表明,该模型与确定性径流预报方法相比,不仅有利于决策人员定量考虑不确定性,而且在期望意义上提高了径流预报精度,具有较高的应用价值。     

基于贝叶斯模型平均的径流模拟及不确定性分析

水文模型是模拟水循环过程重要手段,依靠单个模型进行模拟往往存在很大的不确定性,使通过多模型进行组合模拟成为必然趋势。选取3个集总式水文模型应用贝叶斯模型平均(BMA)进行流域月径流量的多模型模拟,采用期望最大化算法推求BMA分布参数以得到BMA均值模拟序列和90%不确定性区间。以武烈河实测数据为例进行分析,结果表明:BMA方法既能通过均值模拟提供更高精度的模拟效果,还可通过不确定性置信区间定量评价模型结构不确定性,为径流模拟提供丰富信息。     

水文模型参数区域化研究进展

参数区域化方法是解决资料缺乏地区水文模拟和预报的有效手段,主要包括回归法、空间邻近法和属性相似法三类方法,可将有资料流域的水文模型参数移用到资料缺乏流域。首先回顾了区域化方法的基本原理和应用方法,并分析了三类主要区域化方法的适用性。从流域特征因子、水文模型及参数、不确定性探讨三个方面综述了区域化方法的研究进展。分析发现,当前区域化方法缺乏完善的理论基础,流域特征因子选择存在主观性,水文模型及参数的适用性方面研究不足。最后展望了未来的研究重点：（1）多维度适用性比较;（2）水文过程和参数的空间分布规律;（3）参数的尺度问题;（4）参数区域化的不确定性问题。     

缺资料流域水文模型参数区域化研究进展

缺资料流域由于缺乏历史径流资料无法进行水文模型参数率定,因此模型参数识别具有很大的难度和不确定性。目前国内外学者对缺资料流域水文模型参数识别一般采用区域化方法,即通过某种途径,利用有资料流域的模型参数推求缺资料流域的模型参数,从而对缺资料流域进行预报。文章总结分析了缺资料流域水文模型参数区域化方法中的参数移植法和回归法的研究进展,对存在的方法选择问题、尺度问题等进行了讨论,并指出在参数不确定性、尺度转换以及多种信息源利用等方面还有待于进一步研究。     

自适应神经模糊推理系统(ANFIS)在水文模型综合中的应用

由于目前已有很多比较成熟的流域水文模型,因此我们可以选用几个流域水文模型进行并行运算,来同时模拟流域降雨—径流关系。在相同的降雨输入情况下,不同模型得到的模拟流量必然会有所不同,模型效率系数和模拟精度也会不同。因此,如何将不同模型的模拟结果进行综合以进一步提高流量模拟精度是一个关键问题。本文选用自适应神经模糊推理系统(ANFIS)作为水文模型综合平台,以牧马河流域为试验区域,对两个并行运算水文模型(三水源新安江模型和总径流响应模型)的结果进行综合处理,得到了更稳健的流量模拟结果,大大提高了模型效率和模拟精度。该方法值得在实践中借鉴。     

河道洪水实时概率预报模型与应用

通过数据同化方法合理地将实时水文观测数据融入到洪水预报模型中,可提高洪水预报模型的实时性和精确度。选取沿程断面流量、水位和糙率系数作为代表水流状态的基本粒子,以监测断面实测水位数据作为观测信息,建立了基于粒子滤波数据同化算法的河道洪水实时概率预报模型。模型应用于黄河中下游河道洪水预报计算的结果表明,采用粒子滤波方法同化观测水位后,不仅可以直接校正水位,同时也可以有效地校正流量和糙率,为未来时刻模型预报计算提供更准确的水流初始条件和糙率取值区间,进而有效地提高模型预报结果的精度,给出合理的概率预报区间。不同预报期的预报结果表明,随着预报期的增长,同化效果减弱,模型预报结果的精度会有所降低,水位概率预报结果受粒子间糙率不同的影响不确定性增加,而流量概率预报结果受给定模型边界条件的影响不确定性降低。所提出模型可以有效同化真实水位观测数据,适合应用于实际的洪水预报工作中。     

大尺度水文模型参数不确定性分析的挑战与综合研究框架

水文模型是对自然界复杂水文现象与过程的一种综合近似描述,在水旱灾害防治、水资源管理与开发利用等方面应用广泛。本文分析了大尺度水文模型应用的难点,总结了参数不确定性研究的主要进展,介绍了参数不确定性分析框架“敏感性分析—参数优化—参数区域化”(SOR)的基本概念、重要性与应用情况。论文基于已有认识,建议在水文建模优化过程中引入更全面的参数不确定性分析SOR框架,并加强新一代分布式水文模型与更加成熟的水文气象数据观测系统的开发,以减少来自模型结构与模型驱动数据的不确定性,提高全球变化背景下大尺度水文模型水循环过程模拟和预测的准确性。     

Determination of site-specific soil-water characteristic curve from a limited number of test data— A Bayesian perspective

Determining soilewater characteristic curve(SWCC) at a site is an essential step for implementing unsaturated soil mechanics in geotechnical engineering practice, which can be measured directly through various in-situ and/or laboratory tests. Such direct measurements are, however, costly and timeconsuming due to high standards for equipment and procedural control and limits in testing apparatus. As a result, only a limited number of data points(e.g., volumetric water content vs. matric suction)on SWCC at some values of matric suction are obtained in practice. How to use a limited number of data points to estimate the site-specific SWCC and to quantify the uncertainty(or degrees-of-belief) in the estimated SWCC remains a challenging task. This paper proposes a Bayesian approach to determine a site-specific SWCC based on a limited number of test data and prior knowledge(e.g., engineering experience and judgment). The proposed Bayesian approach quantifies the degrees-of-belief on the estimated SWCC according to site-specific test data and prior knowledge, and simultaneously selects a suitable SWCC model from a number of candidates based on the probability logic. To address computational issues involved in Bayesian analyses, Markov Chain Monte Carlo Simulation(MCMCS), specifically Metropolis-Hastings(M-H) algorithm, is used to solve the posterior distribution of SWCC model parameters, and Gaussian copula is applied to evaluating model evidence based on MCMCS samples for selecting the most probable SWCC model from a pool of candidates. This removes one key limitation of the M-H algorithm, making it feasible in Bayesian model selection problems. The proposed approach is illustrated using real data in Unsaturated Soil Database(UNSODA) developed by U.S. Department of Agriculture. It is shown that the proposed approach properly estimates the SWCC based on a limited number of site-specific test data and prior knowledge, and reflects the degrees-of-belief on the estimated SWCC in a rational and quantitative manner.     

Determining the initial stage loss content of tectonic coal gas desorption via the Bayesian probability method

Mine gas disaster prediction and prevention are based on gas content measurement, which results in initial stage loss when determining coal gas desorption contents in engineering applications. We propose a Bayesian probability statistical method in the coal gas desorption model on the basis of constrained prior information. First, we use a self-made coal sample gas desorption device to test initial stage gas desorption data of tectonic coal and undeformed coal. Second, we calculate the initial stage loss of different coal samples with the power exponential function parameters by using Bayesian probability statistics and least squares estimation. Results show that Bayesian probability statistics and least squares estimation can be used to obtain regression and desorption coefficients, thereby illustrating the Bayesian estimation method’s validity and reliability. Given that the Bayesian probability method can apply prior information to constrain the model’s posterior parameters, it provides results that are statistically significant in the initial stage loss of coal gas desorption by connecting observation data and prior information.

     

基于贝叶斯网络的地震相分类

为了解决传统多地震属性的地震相分类方法中"难以引入先验信息用以指导分类,难以给出地震相分类结果可靠程度的定量估计,且各分类参数的权值较难确定"这3个问题,提出了一种新的基于贝叶斯网络的地震相分类方法.该分类方法有效地融合了先验信息和训练样本的分布特征,对提取的多种地震属性进行智能分析,以概率推理的方式得到各地震相类别的概率值,并根据概率分布估计分类结果的可靠程度.详述了贝叶斯网络用于地震相分类的原理与方法,并结合理论地震数据,验证了该方法的可行性和正确性.      

Landslide susceptibility zonation method based on C5.0 decision tree and K-means cluster algorithms to improve the efficiency of risk management

Machine learning algorithms are an important measure with which to perform landslide susceptibility assessments,but most studies use GIS-based classification methods to conduct susceptibility zonation.This study presents a machine learning approach based on the C5.0 decision tree(DT)model and the K-means cluster algorithm to produce a regional landslide susceptibility map.Yanchang County,a typical landslide-prone area located in northwestern China,was taken as the area of interest to introduce the proposed application procedure.A landslide inventory containing 82 landslides was prepared and subse-quently randomly partitioned into two subsets:training data(70％landslide pixels)and validation data(30％landslide pixels).Fourteen landslide influencing factors were considered in the input dataset and were used to calculate the landslide occurrence probability based on the C5.0 decision tree model.Susceptibility zonation was implemented according to the cut-off values calculated by the K-means clus-ter algorithm.The validation results of the model performance analysis showed that the AUC(area under the receiver operating characteristic(ROC)curve)of the proposed model was the highest,reaching 0.88,compared with traditional models(support vector machine(SVM)=0.85,Bayesian network(BN)=0.81,frequency ratio(FR)=0.75,weight of evidence(WOE)=0.76).The landslide frequency ratio and fre-quency density of the high susceptibility zones were 6.76/km2 and 0.88/km2,respectively,which were much higher than those of the low susceptibility zones.The top 20％interval of landslide occurrence probability contained 89％of the historical landslides but only accounted for 10.3％of the total area.Our results indicate that the distribution of high susceptibility zones was more focused without contain-ing more"stable"pixels.Therefore,the obtained susceptibility map is suitable for application to landslide risk management practices.     

Seismic hazard assessment in TERESA test areas based on a Bayesian technique

A method based on Bayesian techniques has been applied to evaluate the seismic hazard in the two test areas selected by the participants in the ESC/SC8-TERESA project: Sannio-Matese in Italy and the northern Rhine region (BGN). A prior site occurrence model (prior SOM) is obtain from a seismicity distribution modeled in wide seismic sources. The posterior occurrence model (posterior SOM) is calculated after a Bayesian correction which, basically, recovers the spatial information of the epicenter distribution and considers attenuation and location errors, not using source zones. The uncertainties of the occurrence probabilities are evaluated in both models.The results are displayed in terms of probability and variation coefficient contour maps for a chosen intensity level, and with plots of mean return period versus intensity in selected test sites, including the 90% probability intervals.It turns out that the posterior SOM gives a better resolution in the probability estimate, decreasing its uncertainty, especially in low seismic activity regions.     

Stochastic prediction in geomorphology using Bayesian inference models

Bayesian inference modeling may be applied to empirical stochastic prediction in geomorphology where outcomes of geomorphic processes can be expressed by probability density functions. Natural variations in process outputs are accommodated by the probability model. Uncertainty in the values of model parameters is reduced by considering statistically independent prior information on long-term, parameter behavior. Formal combination of model and parameter information yields a Bayesian probability distribution that accounts for parameter uncertainty, but not for model uncertainty or systematic error which is ignored herein. Prior information is determined by ordinary objective or subjective methods of geomorphic investigation. Examples involving simple stochastic models are given, as applied to the prediction of shifts in river courses, alpine rock avalanches, and fluctuating river bed levels. Bayesian inference models may be applied spatially and temporally as well as to functions of a random variable. They provide technically superior forecasts, for a given shortterm data set, to those of extrapolation or stochastic simulation models. In applications the contribution of the field geomorphologist is of fundamental quantitative importance.     

突发性水污染事件溯源方法

为快速准确地求解突发性水污染溯源问题,在微分进化与蒙特卡罗基础上提出了一种新的溯源方法。该方法将溯源问题视为贝叶斯估计问题,推导出污染源强度、位置和排放时刻等未知参数的后验概率密度函数;结合微分进化和蒙特卡罗模拟方法对后验概率分布进行采样,进而估计出这些未知参数,确定污染源项。通过算例与贝叶斯-蒙特卡罗方法进行对比,结果表明:该方法可使迭代次数有效缩减3/4,污染源强度、位置和排放时刻的平均相对误差分别减少1.23%、2.23%和4.15%,均值误差分别降低0.39%、0.83%和1.49%,其稳定性和可靠性明显高于贝叶斯-蒙特卡罗方法,能较好地识别突发性水污染源,为解决突发水污染事件中的追踪溯源难点问题提供了新的思路和方法。     

基于逻辑回归的四川青川县区域滑坡灾害预警模型

四川省青川县滑坡灾害群发,点多面广,区域滑坡灾害预警是有效防灾减灾的重要手段,预警模型是成功预警的核心。由于研究区滑坡诱发机理复杂、调查监测大数据及分析方法不足等原因,传统区域地质灾害预警模型存在预警精度有限、精细化不足等问题。文章在青川县地质灾害调查监测和降水监测成果集成整理与数据清洗基础上,构建了青川县区域滑坡灾害训练样本集,样本集包括地质环境、降雨等27个输入特征属性和1个输出特征属性,涵盖了青川县近9年（2010—2018年）全部样本,数量达1 826个（其中,正样本613个,负样本1 213个）。基于逻辑回归算法,对样本集进行5折交叉验证学习训练,采用贝叶斯优化算法进行模型优化,采用精确度、ROC曲线和AUC值等指标校验模型准确度和模型泛化能力。其中,ROC曲线也称为“受试者工作特征”曲线;AUC值表示ROC曲线下的面积。校验结果显示,基于逻辑回归算法的模型训练结果准确率和泛化能力均较好（准确率94.3%,AUC为0.980）。开展区域滑坡实际预警时,按训练样本特征属性格式,输入研究区各预警单元27个特征属性,调用预先学习训练好的模型,输出滑坡灾害发生概率,根据输出概率分段确定滑坡灾害预警等级。当输出概率P≥40%且P<60%时,发布黄色预警;当输出概率P≥60%且P<80%时,发布橙色预警;当输出概率P≥80%时,发布红色预警。     

Impact on loss/risk assessments of inter-model variability in vulnerability analysis

Fragility curves (FCs) constitute an emerging tool for the seismic risk assessment of all elements at risk. They express the probability of a structure being damaged beyond a specific damage state for a given seismic input motion parameter, incorporating the most important sources of uncertainties, that is, seismic demand, capacity and definition of damage states. Nevertheless, the implementation of FCs in loss/risk assessments introduces other important sources of uncertainty, related to the usually limited knowledge about the elements at risk (e.g., inventory, typology). In this paper, within a Bayesian framework, it is developed a general methodology to combine into a single model (Bayesian combined model, BCM) the information provided by multiple FC models, weighting them according to their credibility/applicability, and independent past data. This combination enables to efficiently capture inter-model variability (IMV) and to propagate it into risk/loss assessments, allowing the treatment of a large spectrum of vulnerability-related uncertainties, usually neglected. As case study, FCs for shallow tunnels in alluvial deposits, when subjected to transversal seismic loading, are developed with two conventional procedures, based on a quasi-static numerical approach. Noteworthy, loss/risk assessments resulting from such conventional methods show significant unexpected differences. Conventional fragilities are then combined in a Bayesian framework, in which also probability values are treated as random variables, characterized by their probability density functions. The results show that BCM efficiently projects the whole variability of input models into risk/loss estimations. This demonstrates that BCM is a suitable framework to treat IMV in vulnerability assessments, in a straightforward and explicit manner.     

Bayesian model selection for complex geological structures using polynomial chaos proxy

Different interpretation of sedimentary environments lead to “scenario uncertainty” where the prior reservoir model has a high level of discrete uncertainty. In a real field application, the scenario uncertainty has a considerable effect on flow response uncertainty and makes the uncertainty quantification problem highly nonlinear. We use clustering methods to address the scenario uncertainty. Our approach to cluster analysis is based on the posterior probabilities of models, known as “Bayesian model selection.” Accordingly, we integrate overall possible parameters in each scenario with respect to their corresponding priors to give the measure of how well a model is supported by observations. We propose a cluster-based reduced terms polynomial chaos proxy to efficiently estimate the posterior probability density function under each cluster and calculate the posterior probability of each model. We demonstrate that the convergence rate of the reduced terms polynomial chaos proxy is significantly improved under each cluster comparing to the non-clustered case. We apply the proposed cluster-based polynomial chaos proxy framework to study the plausibility of three training images based on different geological interpretation of the second layer of synthetic Stanford VI reservoir. We demonstrate that the proposed workflow can be efficiently used to calculate the posterior probability of each scenario and also sample from the posterior facies models within each scenario.