Geodetic data inversion based on Bayesian formulation with direct and indirect prior information     

Mitsuhiro Matsu'ura  Akemi Noda  Yukitoshi Fukahata 《Geophysical Journal International》2007,171(3):1342-1351

  相似文献   

A model for estimating the value of sampling programs and the optimal number of samples for contaminated soil   总被引：1，自引：0，他引：1  

Pär-Erik Back 《Environmental Geology》2007,52(3):573-585

A model is presented for estimating the value of information of sampling programs for contaminated soil. The purpose is to calculate the optimal number of samples when the objective is to estimate the mean concentration. A Bayesian risk–cost–benefit decision analysis framework is applied and the approach is design-based. The model explicitly includes sample uncertainty at a complexity level that can be applied to practical contaminated land problems with limited amount of data. Prior information about the contamination level is modelled by probability density functions. The value of information is expressed in monetary terms. The most cost-effective sampling program is the one with the highest expected net value. The model was applied to a contaminated scrap yard in Göteborg, Sweden, contaminated by metals. The optimal number of samples was determined to be in the range of 16–18 for a remediation unit of 100 m². Sensitivity analysis indicates that the perspective of the decision-maker is important, and that the cost of failure and the future land use are the most important factors to consider. The model can also be applied for other sampling problems, for example, sampling and testing of wastes to meet landfill waste acceptance procedures.  相似文献   

Time independent seismic hazard analysis in Alborz and surrounding area   总被引：1，自引：0，他引：1  

Mohammad Ashtari Jafari 《Natural Hazards》2007,42(1):237-252

The Bayesian probability estimation seems to have efficiencies that make it suitable for calculating different parameters of seismicity. Generally this method is able to combine prior information on seismicity while at the same time including statistical uncertainty associated with the estimation of the parameters used to quantify seismicity, in addition to the probabilistic uncertainties associated with the inherent randomness of earthquake occurrence. In this article a time-independent Bayesian approach, which yields the probability that a certain cut-off magnitude will be exceeded at certain time intervals is examined for the region of Alborz, Iran, in order to consider the following consequences for the city of Tehran. This area is located within the Alpine-Himalayan active mountain belt. Many active faults affect the Alborz, most of which are parallel to the range and accommodate the present day oblique convergence across it. Tehran, the capital of Iran, with millions of inhabitants is located near the foothills of the southern Central Alborz. This region has been affected several times by historical and recent earthquakes that confirm the importance of seismic hazard assessment through it. As the first step in this study an updated earthquake catalog is compiled for the Alborz. Then, by assuming a Poisson distribution for the number of earthquakes which occur at a certain time interval, the probabilistic earthquake occurrence is computed by the Bayesian approach. The highest probabilities are found for zone AA and the lowest probabilities for zones KD and CA, meanwhile the overall probability is high.  相似文献   

On the Bayesian analysis of the earthquake hazard in the North-East Indian peninsula     

Imtiyaz A. Parvez 《Natural Hazards》2007,40(2):397-412

The Bayesian extreme-value distribution of earthquake occurrences has been used to estimate the seismic hazard in 12 seismogenic zones of the North-East Indian peninsula. The Bayesian approach has been used very efficiently to combine the prior information on seismicity obtained from geological data with historical observations in many seismogenic zones of the world. The basic parameters to obtain the prior estimate of seismicity are the seismic moment, slip rate, earthquake recurrence rate and magnitude. These estimates are then updated in terms of Bayes’ theorem and historical evaluations of seismicity associated with each zone. From the Bayesian analysis of extreme earthquake occurrences for North-East Indian peninsula, it is found that for T = 5 years, the probability of occurrences of magnitude (M _w = 5.0–5.5) is greater than 0.9 for all zones. For M _w = 6.0, four zones namely Z1 (Central Himalayas), Z5 (Indo-Burma border), Z7 (Burmese arc) and Z8 (Burma region) exhibit high probabilities. Lower probability is shown by some zones namely␣Z4, Z12, and rest of the zones Z2, Z3, Z6, Z9, Z10 and Z11 show moderate probabilities.  相似文献   

Archaeomagnetic secular variation in the UK during the past 4000 years and its application to archaeomagnetic dating     

I. Zananiri  C.M. Batt  D.H. Tarling 《Physics of the Earth and Planetary Interiors》2007,160(2):97-107

This paper examines the limitations and deficiencies of the current British archaeomagnetic calibration curve and applies several mathematical approaches in an attempt to produce an improved secular variation curve for the UK for use in archaeomagnetic dating. The dataset compiled is the most complete available in the UK, incorporating published results, PhD theses and unpublished laboratory reports. It comprises 620 archaeomagnetic (directional) data and 238 direct observations of the geomagnetic field, and includes all relevant information available about the site, the archaeomagnetic direction and the archaeological age. A thorough examination of the data was performed to assess their quality and reliability. Various techniques were employed in order to use the data to construct a secular variation (SV) record: moving window with averaging and median, as well as Bayesian statistical modelling. The SV reference curve obtained for the past 4000 years is very similar to that from France, most differences occurring during the early medieval period (or Dark Ages). Two examples of dating of archaeological structures, medieval and pre-Roman, are presented based on the new SV curve for the UK and the implications for archaeomagnetic dating are discussed.  相似文献   

Maximum-likelihood ambiguity resolution based on Bayesian principle   总被引：2，自引：1，他引：2  

J. Zhu  X. Ding  Y. Chen 《Journal of Geodesy》2001,75(4):175-187

 Based on the Bayesian principle and the fact that GPS carrier-phase ambiguities are integers, the posterior distribution of the ambiguities and the position parameters is derived. This is then used to derive the maximum posterior likelihood solution of the ambiguities. The accuracy of the integer ambiguity solution and the position parameters is also studied according to the posterior distribution. It is found that the accuracy of the integer solution depends not only on the variance of the corresponding float ambiguity solution but also on its values. Received: 27 July 1999 / Accepted: 22 November 2000  相似文献   

Risk assessment for the Yucca Mountain high-level nuclear waste repository site: Estimation of volcanic disruption     

Chih-Hsiang Ho 《Mathematical Geology》1992,24(4):347-364

In this article, we model the volcanism near the proposed nuclear waste repository at Yucca Mountain, Nevada, U.S.A. by estimating the instantaneous recurrence rate using a nonhomogeneous Poisson process with Weibull intensity and by using a homogeneous Poisson process to predict future eruptions. We then quantify the probability that any single eruption is disruptive in terms of a (prior) probability distribution, since not every eruption would result in disruption of the repository. Bayesian analysis is performed to evaluate the volcanic risk. Based on the Quaternary data, a 90% confidence interval for the instantaneous recurrence rate near the Yucca Mountain site is (1.85×10^–6/yr, 1.26×10^–5/yr). Also, using-these confidence bounds, the corresponding 90% confidence interval for the risk (probability of at least one disruptive eruption) for an isolation time of 10⁴ years is (1.0×10^–3, 6.7×10^–3), if it is assumed that the intensity remains constant during the projected time frame.  相似文献   

基于改进全变差正则化的GPR多尺度全波形双参数同步反演          下载免费PDF全文

王珣  冯德山  王向宇 《地球物理学报》1954,63(12):4485-4501

针对探地雷达（GPR）双参数全波形反演中电导率反演精度差、双参数存在串扰现象、反演计算量大、易陷入局部极值等问题.作者将具有多参数调节功能的L-BFGS算法引入到GPR时间域全波形反演中,它避免了对Hessian矩阵的直接存储与精确求解,减小了存储量和计算量.结合参数调节因子的选取,有效减小了同步反演时介电常数与电导率的串扰影响,在不降低介电常数反演精度的前提下,提高电导率参数的反演精度.通过在反演目标函数中加载改进全变差正则化方法,提高了反演的稳定性,使目标体边缘轮廓更加清晰.首先以简单模型为例,对比了单尺度反演与多尺度串行反演策略的优劣,说明多尺度串行反演有利于逐步搜索全局最优解;而开展参数调节因子的选取实验,说明合适的参数调节因子可以有效改善介质电导率的反演精度;测试了不同正则化的反演效果,表明改进全变差正则化能提高反演稳定性,显著降低模型重构误差.最后,分别对含噪合成数据和实测数据进行了反演测试,说明本文提出的多尺度、双参数反演具有较强的鲁棒性,能提供更丰富的信息约束,重构图像界面清晰、反演效果好.  相似文献   

均方误差意义下的正则化参数二次优化方法     

林东方  朱建军  付海强  张兵 《测绘学报》2020,49(4):443-451

Tikhonov正则化法是大地测量中应用最为广泛的病态问题解算方法之一。影响正则化法解算效果的重要因素是正则化参数,然而,最优正则化参数的确定一直是正则化解算的难题,如L曲线法确定的正则化参数具有稳定性好、可靠性高的优点,但存在过度平滑问题,导致正则化法对模型参数估值精度改善较小。本文从均方误差角度分析了正则化参数对模型参数估计质量的影响。基于奇异值分解技术,提出了由模型参数投影值分块计算均方误差的方法,避免了均方误差迭代计算,并基于均方误差最小准则给出了正则化参数优化方法,实现了对L曲线正则化参数的优化。数值模拟试验与PolInSAR植被高反演试验结果表明,正则化参数优化方法有效改善了正则化法解算效果,提高了模型参数估计精度。  相似文献   

Multiobjective analysis of chaotic dynamic systems with sparse learning machines     

Abedalrazq F. Khalil  Mac McKeeMariush Kemblowski  Tirusew AsefaLuis Bastidas 《Advances in water resources》2006

Sparse learning machines provide a viable framework for modeling chaotic time-series systems. A powerful state-space reconstruction methodology using both support vector machines (SVM) and relevance vector machines (RVM) within a multiobjective optimization framework is presented in this paper. The utility and practicality of the proposed approaches have been demonstrated on the time series of the Great Salt Lake (GSL) biweekly volumes from 1848 to 2004. A comparison of the two methods is made based on their predictive power and robustness. The reconstruction of the dynamics of the Great Salt Lake volume time series is attained using the most relevant feature subset of the training data. In this paper, efforts are also made to assess the uncertainty and robustness of the machines in learning and forecasting as a function of model structure, model parameters, and bootstrapping samples. The resulting model will normally have a structure, including parameterization, that suits the information content of the available data, and can be used to develop time series forecasts for multiple lead times ranging from two weeks to several months.  相似文献