基于ICCP算法的重力辅助惯性导航

迭代最近等值线算法(ICCP)是一种重要的匹配导航算法,文中首先介绍ICCP算法的基本原理,随后在0.2′×0.2′重力异常数据库的基础上,利用ICCP算法进行仿真计算得到最佳匹配位置。最后为了验证匹配位置是否可用于修正惯导误差,提出将匹配位置误差作为观测量,用卡尔曼滤波对惯导系统误差进行最优估计。由最后的仿真结果可以看出,ICCP算法可有效抑制惯导纬度误差的增长,且最大纬度误差不超过2,′以匹配位置误差作为观测量可以用来估计惯导方位误差角。     

顾及DEM误差自相关的坡度计算模型精度分析

基于DEM的坡度计算,其误差来源于DEM误差、DEM结构和坡度计算模型。在顾及DEM误差自相关的前提下,对四种DEM坡度计算模型进行了分析和评价。研究表明,三阶不带权差分能给出较高的坡度计算精度;在局部窗口中,格网点数量越多,坡度计算越准确;等权比不等权的坡度计算模型更准确;DEM误差自相关结构形式对坡度计算无影响。进一步的理论分析和试验分析还表明：DEM误差自相关性的存在,不仅能够改善地形分析的精度,也能改善DEM自身精度。     

基于规则库和网格算法的土地利用现状图自动数字注记

图形自动注记可极大地提高制图的工作效率。针对土地利用现状图数字注记的特点,本文从地图认知角度出发,建立注记自动配置的规则库,提出面状要素自动注记的一种新的算法——网格法。该算法最大的优点是容易实现、效率高,可方便地解决注记压盖冲突等问题。通过开发一个程序实例,对图斑要素进行了自动数字注记,实验表明95％以上的图形注记效果都较为理想。     

Aspherical correction for free vibration of elastically anisotropic solid by means of theXYZ algorithm

By applying the perturbation theory to theXYZ algorithm (a kind of variational method), the difference f in free vibration frequencies between sphere and ellipsoid was approximated as , where i and _i (i = x,y andz) (i=x, y andz) are aspherical coefficients and asphericities of the ellipsoid, respectively. We developed an analytic method to compute the aspherical coefficients719-4 by using theXYZ algorithm. A numerical example was given for an ellipsoidal olivine, and an attempt was made to estimate the asphericities of the specimen by a least-squares method, based on the relationship between frequency shift and asphericity.     

基于改进PSO-SVM的噪声稳健建筑物变形监测方法

对高层、超高层建筑物进行实时,高精度的变形监测对提前预防安全隐患,保证人民生命财产安全具有重要意义。建筑物变形作为一种典型的随机性和微弱性过程,噪声等误差的存在会影响从中提取有用的变形信息。针对该问题,提出一种改进粒子群(Particle Swarm Optimization,PSO)算法优化支持向量机(Support Vector Machine,SVM)的噪声稳健建筑物变形监测方法,利用改进PSO算法的全局搜索能力对SVM的核参数进行优化,提升预测精度的同时增强算法的噪声稳健性。基于实测数据的试验结果表明,相对于传统交叉验证SVM和小波方法,所提方法可以获得更高的变形预测精度,并且在低信噪比条件下优势更加明显。     

Estimation for partially observed Markov processes

Many stochastic process models for environmental data sets assume a process of relatively simple structure which is in some sense partially observed. That is, there is an underlying process (X_n, n 0) or (X_t, t 0) for which the parameters are of interest and physically meaningful, and an observable process (Y_n, n 0) or (Y_t, t 0) which depends on the X process but not otherwise on those parameters. Examples are wide ranging: the Y process may be the X process with missing observations; the Y process may be the X process observed with a noise component; the X process might constitute a random environment for the Y process, as with hidden Markov models; the Y process might be a lower dimensional function or reduction of the X process. In principle, maximum likelihood estimation for the X process parameters can be carried out by some form of the EM algorithm applied to the Y process data. In the paper we review some current methods for exact and approximate maximum likelihood estimation. We illustrate some of the issues by considering how to estimate the parameters of a stochastic Nash cascade model for runoff. In the case of k reservoirs, the outputs of these reservoirs form a k dimensional vector Markov process, of which only the kth coordinate process is observed, usually at a discrete sample of time points.     

Kriging with imprecise (fuzzy) variograms. I: Theory

Imprecise variogram parameters are modeled with fuzzy set theory. The fit of a variogram model to experimental variograms is often subjective. The accuracy of the fit is modeled with imprecise variogram parameters. Measurement data often are insufficient to create good experimental variograms. In this case, prior knowledge and experience can contribute to determination of the variogram model parameters. A methodology for kriging with imprecise variogram parameters is developed. Both kriged values and estimation variances are calculated as fuzzy numbers and characterized by their membership functions. Besides estimation variance, the membership functions are used to create another uncertainty measure. This measure depends on both homogeneity and configuration of the data.     

Reservoir description using dynamic parameterisation selection with a combined stochastic and gradient search

There is a correspondence between flow in a reservoir and large scale permeability trends. This correspondence can be derived by constraining reservoir models using observed production data. One of the challenges in deriving the permeability distribution of a field using production data involves determination of the scale of resolution of the permeability. The Adaptive Multiscale Estimation (AME) seeks to overcome the problems related to choosing the resolution of the permeability field by a dynamic parameterisation selection. The standard AME uses a gradient algorithm in solving several optimisation problems with increasing permeability resolution. This paper presents a hybrid algorithm which combines a gradient search and a stochastic algorithm to improve the robustness of the dynamic parameterisation selection. At low dimension, we use the stochastic algorithm to generate several optimised models. We use information from all these produced models to find new optimal refinements, and start out new optimisations with several unequally suggested parameterisations. At higher dimensions we change to a gradient-type optimiser, where the initial solution is chosen from the ensemble of models suggested by the stochastic algorithm. The selection is based on a predefined criterion. We demonstrate the robustness of the hybrid algorithm on sample synthetic cases, which most of them were considered insolvable using the standard AME algorithm.     

Temporal behavior of the background seismicity rate in central Japan, 1998 to mid-2003

We studied the temporal behavior of the background shallow seismicity rate in 700 circular areas across inland Japan. To search for and test the significance of the possible rate changes in background seismicity, we developed an efficient computational method that applies the space–time ETAS model proposed by Ogata in 1998 to the areas. Also, we conducted Monte Carlo tests using a simulated catalog to validate the model we applied. Our first finding was that the activation anomalies were found so frequently that the constant background seismicity hypothesis may not be appropriate and/or the triggered event model with constraints on the parameters may not adequately describe the observed seismicity. However, quiescence occasionally occurs merely by chance. Another outcome of our study was that we could automatically find several anomalous background seismicity rate changes associated with the occurrence of large earthquakes. Very significant seismic activation was found before the M6.1 Mt. Iwate earthquake of 1998. Also, possible seismic quiescence was found in an area 150 km southwest of the focal region of the M7.3 Western Tottori earthquake of 2000. The seismicity rate in the area recovered after the mainshock.     

T -distributed Random Fields: A Parametric Model for Heavy-tailedWell-log Data1

Histograms of observations from spatial phenomena are often found to be more heavy-tailed than Gaussian distributions, which makes the Gaussian random field model unsuited. A T-distributed random field model with heavy-tailed marginal probability density functions is defined. The model is a generalization of the familiar Student-T distribution, and it may be given a Bayesian interpretation. The increased variability appears cross-realizations, contrary to in-realizations, since all realizations are Gaussian-like with varying variance between realizations. The T-distributed random field model is analytically tractable and the conditional model is developed, which provides algorithms for conditional simulation and prediction, so-called T-kriging. The model compares favourably with most previously defined random field models. The Gaussian random field model appears as a special, limiting case of the T-distributed random field model. The model is particularly useful whenever multiple, sparsely sampled realizations of the random field are available, and is clearly favourable to the Gaussian model in this case. The properties of the T-distributed random field model is demonstrated on well log observations from the Gullfaks field in the North Sea. The predictions correspond to traditional kriging predictions, while the associated prediction variances are more representative, as they are layer specific and include uncertainty caused by using variance estimates.