Disjunctive Kriging with Hard and Imprecise Data

This paper presents a methodology for assessing local probability distributions by disjunctive kriging when the available data set contains some imprecise measurements, like noisy or soft information or interval constraints. The basic idea consists in replacing the set of imprecise data by a set of pseudohard data simulated from their posterior distribution; an iterative algorithm based on the Gibbs sampler is proposed to achieve such a simulation step. The whole procedure is repeated many times and the final result is the average of the disjunctive kriging estimates computed from each simulated data set. Being data-independent, the kriging weights need to be calculated only once, which enables fast computing. The simulation procedure requires encoding each datum as a pre-posterior distribution and assuming a Markov property to allow the updating of pre-posterior distributions into posterior ones. Although it suffers some imperfections, disjunctive kriging turns out to be a much more flexible approach than conditional expectation, because of the vast class of models that allows its computation, namely isofactorial models.     

洪家渡水电站左坝肩高边坡稳定及加固措施研究

贵州洪家渡水电站左坝肩高边坡开挖高度310 m,采用弹塑性有限元数值模拟方法对左坝肩高边坡施工期岩体应力应变状态进行分析研究,并结合现场变形观测资料对计算结果进行验证,从而对左坝肩高边坡的稳定性进行评价并提出相应的加固措施。     

对下扬子与华南边界结合带东延问题的地球物理探讨

以江山-绍兴、铅山-宜春等断裂带作为华南与扬子块体之间的边界结合带，这种认识目前已基本为大家所接受．但是，这条边界结合带向东延伸入海之后的位置与去向一直是多年来研究的热点．本文基于黄、东海研究区的地球物理数据(空间重力数据、布格重力数据和地震层析成像结果)，利用方向导数等处理方法，对研究区的地球物理数据进行处理，并对研究区的地球物理场进行了分析，划分出不同的区块．结合研究区的磁力与地质资料，利用各种成图、成像技术，形成一系列分析图件．在此基础上对华南与扬子块体之间的边界结合带进行了追踪．研究结果认为：该结合带的位置有可能比传统认识中自长江口至大黑山群岛的位置更向南一些，在杭州湾-长崎、对马海峡一线，并呈现向北略微凸出的弓形．边界结合带在深度上属于深大断裂带，一直可以追索到上地幔的顶部．与中朝与扬子块体之间的边界结合带相比较，本条结合带的踪迹不是十分清晰．表明加里东运动之后，该结合带的运动明显减弱．     

科学计算可视化技术的研究与应用

本文介绍了科学计算可视化技术的发展和应用情况，论述了这一新技术的特点及其在油气藏描述，工程设计ＣＡＤ等学科的应用和发展趋势。阐述了目前国际国内计算机图形学领域具体的研究方向和方法及图形软件、图形硬件的发展概况。     

Recovering magnetic susceptibility from electromagnetic data over a one-dimensional earth

While the inversion of electromagnetic data to recover electrical conductivity has received much attention, the inversion of those data to recover magnetic susceptibility has not been fully studied. In this paper we invert frequency-domain electromagnetic (EM) data from a horizontal coplanar system to recover a 1-D distribution of magnetic susceptibility under the assumption that the electrical conductivity is known. The inversion is carried out by dividing the earth into layers of constant susceptibility and minimizing an objective function of the susceptibility subject to fitting the data. An adjoint Green's function solution is used in the calculation of sensitivities, and it is apparent that the sensitivity problem is driven by three sources. One of the sources is the scaled electric field in the layer of interest, and the other two, related to effective magnetic charges, are located at the upper and lower boundaries of the layer. These charges give rise to a frequency-independent term in the sensitivities. Because different frequencies penetrate to different depths in the earth, the EM data contain inherent information about the depth distribution of susceptibility. This contrasts with static field measurements, which can be reproduced by a surface layer of magnetization. We illustrate the effectiveness of the inversion algorithm on synthetic and field data and show also the importance of knowing the background conductivity. In practical circumstances, where there is no a priori information about conductivity distribution, a simultaneous inversion of EM data to recover both electrical conductivity and susceptibility will be required.     

Field vibration mitigation by honeycomb WIB for pile foundations of a high-speed train viaduct

Environmental vibrations from recent high-speed trains are becoming a special concern in the civil and environmental engineering field since they can give detrimental effects to residents, sensitive equipments and high-tech production facilities in the vicinity of train tacks. Herein, aiming at the vibration mitigation for a specific high-tech industrial area, the low-frequency vibrations from a train viaduct are targeted over an anticipated range. A theoretically designed innovative countermeasure, called honeycomb wave impeding barrier (WIB) for a wave impeding barrier, is introduced and its effects are investigated by computer simulation. The present WIB is based on the wave dispersion phenomenon that can modulate the incoming wavelengths into the shorter wavelengths, creating an apparent wave cut-off characteristic in the wave field across WIB installation area. The shorter wavelengths are further impeded due to the impedance ratio of the WIB walls and in-fill materials and absorbed by the in-fills more. The three-dimensional FEM simulation demonstrates a dramatic reduction effect that is difficult to achieve by conventional measures.     

运用改进系统建模法对南海气象数据的建模研究

在系统建模理论的基础上,运用改进的动态数据建模方法,对南海气象数据中的温度进行建模并验证了模型的适用性.根据模型推导出格林函数、逆函数和自协方差函数等,并讨论了南海气象数据中温度模型的稳定性、可逆性和合理性.对系统的频率特性和谱函数进行分析讨论,并给出建模过程中的一些图像.根据模型的适用性检验发现,对所研究的气象数据而言,ARMA(4, 3)模型是最合适的,具有平稳可逆性.所有的建模和分析过程在MATLAB上实现.实验结果表明这种建模方案简便易行,能够快速准确地确定系统的合理模型.     

航空物探平剖图自动绘制的可视化处理

简要介绍了航空物探平面剖面图的自动绘制方法,给出了绘制中对一些问题的处理方法;重点介绍了彩色渐变填充着色、彩色渐变渲染着色等几种可视化处理方法,并提出了着色中渐变度、透明度概念及应用方法。     

变形监测基准点的稳定性检验

基准点的稳定性评价，是变形观测数据处理时不可忽视的重要内容。该文对基准点构网和不构网2种情况时的数据处理方法进行了探讨。     

A New Model for Quantifying Anisotropic Scale Invariance and for Decomposition of Mixing Patterns

A new power–law function has been derived to represent the relationship between area of the set consisting of wave numbers with spectral energy density above S (A(>S)) on the two-dimensional frequency plane and S. The power–law relation holds if the field concerned possessing isotropic scale invariance or generalized scaling invariance involves rotational and ratio-scale changing transforms. The equation is valid for dealing with common exploration geophysical and geochemical fields encountered in mineral exploration and environmental assessment. This power–law function not only provides a new model for characterizing anisotropic scaling invariance for generalized scaling field, for example, estimating the power exponent of power spectrum of generalized scale invariance measure in frequency domain, but also forms a theoretical base for the S–A filtering technique developed for decomposing a mixing field into components on the basis of distinct scaling properties in the frequency domain. It is demonstrated that the method has potential to become a general technique for image processing and pattern recognition.