On Modelling Discrete Geological Structures as Markov Random Fields

The purpose of this paper is to extend the locally based prediction methodology of BayMar to a global one by modelling discrete spatial structures as Markov random fields. BayMar uses one-dimensional Markov-properties for estimating spatial correlation and Bayesian updating for locally integrating prior and additional information. The methodology of this paper introduces a new estimator of the field parameters based on the maximum likelihood technique for one-dimensional Markov chains. This makes the estimator straightforward to calculate also when there is a large amount of missing observations, which often is the case in geological applications. We make simulations (both unconditional and conditional on the observed data) and maximum a posteriori predictions (restorations) of the non-observed data using Markov chain Monte Carlo methods, in the restoration case by employing simulated annealing. The described method gives satisfactory predictions, while more work is needed in order to simulate, since it appears to have a tendency to overestimate strong spatial dependence. It provides an important development compared to the BayMar-methodology by facilitating global predictions and improved use of sparse data.     

极大熵准则在降雨强度-历时-频率模型参数估计中的应用

根据频率与重现期的关系推导出三类常用降雨强度-历时-频率模型的无条件及条件概率分布(概率密度)函数及与模型相对应的约束。极大熵与极大似然准则产生一致估计,本文尝试基于极大熵准则建立降雨强度-历时-频率模型参数估计的优化模型,应用模拟退火算法求解该优化模型。根据比较发现,极大熵估计有时比常用的极大似然估计和最小二乘估计更精确。     

混合模型试验等效水深截断系统优化设计研究

进行了深海海洋平台混合模型等效水深截断系统优化设计试验研究,优化方法选用混合离散变量模拟退火法,目标函数为截断系统和原系统的静力特性相似程度。以一个转塔式系泊浮式生产系统为例进行了等效水深截断系统优化设计计算。     

Effective medium theory and multiple linear regression-based velocity estimation in syn-rift clastic sequence of the Krishna–Godavari basin,India

An inclusion model, based on the Kuster–Toksöz effective medium theory along with Gassmann theory, is tested to forward model velocities for fluid-saturated rocks. A simulated annealing algorithm, along with the inclusion model, effectively inverts measured compressional velocity (V_P) to achieve an effective pore aspect ratio at each depth in a depth variant manner, continuously along with depth. Early Cretaceous syn-rift clastic sediments at two different depth intervals from two wells [well A (2160–2274 m) and well B (5222–5303 m)], in the Krishna–Godavari basin, India, are used for this study. Shear velocity (V_S) estimated using modelled pore aspect ratio offers a high correlation coefficient (>0.95 for both the wells) with measured data. The modelled pore aspect ratio distribution suggests the decrease in pore aspect ratio for the deeper interval, mainly due to increased effective vertical stress. The pore aspect ratio analysis in relation to total porosity and volume of clay reveals that the clay volume has insignificant influence in shaping the pore geometry in the studied intervals. An approach based on multiple linear regression method effectively predicts velocity as a linear function of total porosity, the volume of clay and the modelled pore-space aspect ratio of the rock. We achieved a significant match between measured and predicted velocities. The correlation coefficients between measured and modelled velocities are considerably high (approximately 0.85 and 0.8, for V_P and V_S, respectively). This process indicates the possible influence of pore geometry along with total porosity and volume of clay on velocity.     

Study of 4H-SiC junction barrier Schottky diode using field guard ring termination

This paper reports that the 4H-SiC Schottky barrier diode, PiN diode and junction barrier Schottky diode terminated by field guard rings are designed, fabricated and characterised. The measurements for forward and reverse characteristics have been done, and by comparison with each other, it shows that junction barrier Schottky diode has a lower reverse current density than that of the Schottky barrier diode and a higher forward drop than that of the PiN diode. High-temperature annealing is presented in this paper as well to figure out an optimised processing. The barrier height of 0.79 eV is formed with Ti in this work, the forward drop for the Schottky diode is 2.1 V, with an ideality factor of 3.2, and junction barrier Schottky diode with blocking voltage higher than 400 V was achieved by using field guard ring termination.     

SA算法优化BP神经网络应用于GPS高程拟合

BP神经网络算法用于高程拟合有训练速度慢和易于陷入局部最小值的缺点,基于BP算法学习特点,模拟退火算法（SA）在局部极小处的概率突变性,本文作者有效结合BP和SA算法,提出一种SA优化BP神经网络算法的BPSA混合学习策略,并以实例验证了该算法的有效性。     

区域土地灌溉环境评价与决策分析

国土资源是人们赖以生存和发展的物质基础,土地灌溉环境作为国土资源调查的重要内容,其关键基础设施机井的空间配置对于农业生产具有重要的现实意义。本文以北京市通州区永乐店镇为研究区,基于国土资源调查的多源空间数据,提出了一种对现有土地灌溉环境进行定量评价和优化决策的方法。结果表明,研究区既有机井配置有待进一步改进,优化方案给出的183口新增机井合理分配在优化单元内,能够满足当地农民的土地灌溉需求,同时为土地灌溉环境中机井的更新配置提供了依据。     

复合形模拟退火算法及其在水泥土墙优化设计中的应用

将传统优化方法中的复合形法和智能优化方法中的模拟退火算法有机结合,得到一种新型的智能计算方法-复合形模拟退火算法。仿真分析表明它比单一的复合形法和模拟退火算法搜索性能更优。按照我国现行规范建立了水泥土墙优化设计数学模型,探讨了提出的复合形模拟退火算法在水泥土墙支护结构参数优化设计中的应用。给出了工程计算实例,其计算结果验证了方法的可行性和可靠性。还讨论了基坑内被动土压力区的加固对水泥土墙的嵌固深度和有效宽度的影响。     

遗传模拟退火算法在边坡稳定分析中的应用

利用遗传模拟退火算法结合瑞典圆弧法,寻找最危险滑裂面进行边坡稳定分析。对构造数学模型、确定计算步骤、选取参数等方面作了一定的研究。通过实例计算,结果令人满意。     

Solutions of the inherent problem of the equivalence in direct current resistivity and electromagnetic methods through global optimization and joint inversion by successive refinement of model space

The problem of equivalence in direct current (DC) resistivity and electromagnetic methods for a thin resistive and conducting layer is well‐known. Attempts have been made in the past to resolve this problem through joint inversion. However, equivalence still remains an unresolved problem. In the present study, an effort is made to reduce non‐uniqueness due to equivalence using global optimization and joint inversion by successive refinement of the model space. A number of solutions derived for DC resistivity data using very fast simulated annealing global inversion that fits the observations equally well, follow the equivalence principle and show a definite trend. For a thin conductive layer, the quotient between resistivity and thickness is constant, while for a resistive one, the product between these magnitudes is constant. Three approaches to obtain very fast simulated annealing solutions are tested. In the first one, layer resistivities and thicknesses are optimized in a linear domain. In the second, layer resistivities are optimized in the logarithmic domain and thicknesses in the linear domain. Lastly, both layer resistivities and thicknesses are optimized in the logarithmic domain. Only model data from the mean models, corresponding to very fast simulated annealing solutions obtained for approach three, always fit the observations. The mean model defined by multiple very fast simulated annealing solutions shows extremely large uncertainty (almost 100%) in the final solution after inversion of individual DC resistivity or electromagnetic (EM) data sets. Uncertainty associated with the intermediate resistive and conducting layers after global optimization and joint inversion is still large. In order to reduce the large uncertainty associated with the intermediate layer, global optimization is performed over several iterations by reducing and redefining the search limits of model parameters according to the uncertainty in the solution. The new minimum and maximum limits are obtained from the uncertainty in the previous iteration. Though the misfit error reduces in the solution after successive refinement of the model space in individual inversion, it is observed that the mean model drifts away from the actual model. However, successive refinement of the model space using global optimization and joint inversion reduces uncertainty to a very low level in 4–5 iterations. This approach works very well in resolving the problem of equivalence for resistive as well as for conducting layers. The efficacy of the approach has been demonstrated using DC resistivity and EM data, however, it can be applied to any geophysical data to solve the inherent ambiguities in the interpretations.