Efficient uncertainty quantification techniques in inverse problems for Richards’ equation using coarse-scale simulation models

This paper concerns efficient uncertainty quantification techniques in inverse problems for Richards’ equation which use coarse-scale simulation models. We consider the problem of determining saturated hydraulic conductivity fields conditioned to some integrated response. We use a stochastic parameterization of the saturated hydraulic conductivity and sample using Markov chain Monte Carlo methods (MCMC). The main advantage of the method presented in this paper is the use of multiscale methods within an MCMC method based on Langevin diffusion. Additionally, we discuss techniques to combine multiscale methods with stochastic solution techniques, specifically sparse grid collocation methods. We show that the proposed algorithms dramatically reduce the computational cost associated with traditional Langevin MCMC methods while providing similar sampling performance.     

�й���½�ٶȳ�ģ�ͽ��������о�

???й?????????????????????????????????5?????壬????????????????С???????÷?????С?????????????????????3????????????????о??й?????????????????????????????????????????????????????С???????÷????????壨??????????????????壩????????????????????|?????????????????????????????????????????????????С???????÷?????й????1°×1°?????????????????????????????????С??     

基于方差分量估计的拟合推估及其在GIS误差纠正的应用

拟合推估解算必须首先求得信号向量的方差协方差矩阵,该协方差矩阵一般通过选定的协方差函数,并通过已测点数据进行拟合得到。显然观测噪声的先验方差协方差阵与拟合得到的随机信号的方差协方差矩阵必须相互协调,即观测噪声向量和信号向量的权矩阵所对应的方差因子应该一致,否则将对固定效应和随机效应参数的估计带来系统性的影响。应用方差分量估计来协调拟合推估模型中观测噪声和信号向量的随机模型,并分别从极大似然估计、MINQUE估计、赫尔默特方差分量估计三方面构建了拟合推估模型的方差分量解,最后利用新提出的理论与方法,对一幅实际的扫描地形图进行误差纠正,结果表明基于方差分量估计的拟合推估法能够提高扫描地形图的精度。     

多波束测深数据处理的抗差最小二乘配置迭代解法

针对利用最小二乘配置处理多波束测深数据,存在二次曲面数学模型通常无法精确表征海底地形的整体变化趋势以及观测数据存在粗差或异常点时,常规方法给出的协方差函数不能精确表征其统计特性的问题,本文提出了一种抗差最小二乘配置迭代解法。该方法首先进行协方差函数和观测值方差阵初始化,以多面函数拟合趋势项,然后应用等价权抗差估计并通过迭代计算,最终给出稳健的协方差函数参数解及最小二乘配置解。利用本文提出的方法及传统的方法处理实测的多波束测深数据,试验结果表明,相比于传统的方法,本文提出的方法能够较好地表征海底地形的整体变化趋势,一定程度上克服了多波束测深数据中粗差或异常点的影响。相比于传统的抗差方法,本文方法更为有效地识别出测深数据中异常点,推估效果较好,具有稳健性。     

湖南地区似大地水准面模型再精化研究

湖南地区2007版似大地水准面模型存在数据源质量相对不高、局部区域存在粗差等问题。通过引入EGM2008地球重力场模型、高分辨率DEM数据、高精度GPS/水准数据等新数据源,采用最小二乘配置法与格林基函数样条插值算法,在原有模型的基础上得到1′×1′分辨率的再精化新模型。精度检测结果表明,新模型的内、外符合精度相较2007模型提升80.7%和22.9%,有效降低汉寿等地区的残差,对湖南省地区的复杂地形具有更好的适用性。     

Computation of spherical harmonic coefficients and their error estimates using least-squares collocation

 Equations expressing the covariances between spherical harmonic coefficients and linear functionals applied on the anomalous gravity potential, T, are derived. The functionals are the evaluation functionals, and those associated with first- and second-order derivatives of T. These equations form the basis for the prediction of spherical harmonic coefficients using least-squares collocation (LSC). The equations were implemented in the GRAVSOFT program GEOCOL. Initially, tests using EGM96 were performed using global and regional sets of geoid heights, gravity anomalies and second-order vertical gravity gradients at ground level and at altitude. The global tests confirm that coefficients may be estimated consistently using LSC while the error estimates are much too large for the lower-order coefficients. The validity of an error estimate calculated using LSC with an isotropic covariance function is based on a hypothesis that the coefficients of a specific degree all belong to the same normal distribution. However, the coefficients of lower degree do not fulfil this, and this seems to be the reason for the too-pessimistic error estimates. In order to test this the coefficients of EGM96 were perturbed, so that the pertubations for a specific degree all belonged to a normal distribution with the variance equal to the mean error variance of the coefficients. The pertubations were used to generate residual geoid heights, gravity anomalies and second-order vertical gravity gradients. These data were then used to calculate estimates of the perturbed coefficients as well as error estimates of the quantities, which now have a very good agreement with the errors computed from the simulated observed minus calculated coefficients. Tests with regionally distributed data showed that long-wavelength information is lost, but also that it seems to be recovered for specific coefficients depending on where the data are located. Received: 3 February 2000 / Accepted: 23 October 2000     

不同饵料搭配对日本对虾人工育苗存活率的影响

本文着眼于生产需求，以螺旋藻粉+虾片+日配饵料0号为Z     

区域水资源优化配置多目标模型求解的Matlab与Excel集成实现

区域水资源优化配置模型是决策变量多、结构复杂的非线性多目标模型,简洁快速的求解方法是模型用于实践的关键。利用excel和matlab的特点,通过excellink工具建立excel和matlab的数据通道,实现两者之间的通信,把excel作为数据输入端,调用matlab的函数求解多目标规划问题,实现了对区域水资源配置多目标模型的求解.以郑州市水资源优化配置为例进行了实践应用,应用结果表明,求解方法具有操作简便、易于实现、通用性强的特点。     

The Adaptive Wavelet Collocation Method and Its Application in Front Simulation

The adaptive wavelet collocation method(AWCM)is a variable grid technology for solving partial differential equations(PDEs)with high singularities.Based on interpolating wavelets,the AWCM adapts the grid so that a higher resolution is automatically attributed to domain regions with high singularities. Accuracy problems with the AWCM have been reported in the literature,and in this paper problems of effciency with the AWCM are discussed in detail through a simple one-dimensional(1D)nonlinear advection equation whose analytic solution is easily obtained.A simple and effcient implementation of the AWCM is investigated.Through studying the maximum errors at the moment of frontogenesis of the 1D nonlinear advection equation with different initial values and a comparison with the finite difference method(FDM) on a uniform grid,the AWCM shows good potential for modeling the front effciently.The AWCM is also applied to a two-dimensional(2D)unbalanced frontogenesis model in its first attempt at numerical simulation of a meteorological front.Some important characteristics about the model are revealed by the new scheme.