基于C＋＋类的快速实现测量数据数值计算的方法研究

采用面向对象的程序设计思想,基于C＋＋语言实现了矩阵类,避免了用二维数组来处理矩阵运算时遇到的大量繁琐问题。同时,在矩阵类的基础上实现了线性方程组类,对几种常用的求解线性方程组的算法进行了封装,便于进行频繁的调用。     

基于速度估计的改进型线性变换有限差分偏移在探地雷达中的应用

为了解决传统线性变换有限差分偏移算法对变速介质的高敏感性问题,对加入了速度变量的标量波动方程进行线性变换,在横向非均匀介质条件下,详细推导了包含混合因子项的线性变换差分偏移方程,实现了横向非均匀介质中陡倾角反射界面的正确归位,同时改进型算法可以灵活地选择混合因子项,这样就增强了原算法对复杂地质问题的适应能力,在此基础上...     

不规则区域面积的矩形四等分割计算法

如何对图纸上的区域面积进行计算,尤其是曲边不规则区域的面积计算,是一项基础工作。本文在详细分析了现有区域面积计算法的基础上,提出了一种应用于图纸上不规则区域的面积计算方法。应用本文提出的不规则区域面积的矩形四等分割计算法逐次将目标区域所在的矩形四等分割,直到满足精度阀值要求为止,将区域内所有有效值相加,得出区域总面积。...     

A statistical assessment of the uncertainty in a 3-D geological framework model

Three-dimensional framework models are the state of the art to present geologists’ understanding of a region in a form that can be used to support planning and decision making. However, there is little information on the uncertainty of such framework models. This paper reports an experiment in which five geologists each produced a framework model of a single region in the east of England. Each modeller was provided with a unique set of borehole observations from which to make their model. Each set was made by withholding five unique validation boreholes from the set of all available boreholes. The models could then be compared with the validation observations. There was no significant between-modeller source of variation in framework model error. There was no evidence of systematic bias in the modelled depth for any unit, and a statistically significant but small tendency for the mean error to increase with depth below the surface. The confidence interval for the predicted height of a surface at a point ranged from ±5.6 m to ±6.4 m. There was some evidence that the variance of the model error increased with depth, but no evidence that it differed between modellers or varied with the number of close-neighbouring boreholes or distance to the outcrop. These results are specific to the area that has been modelled, with relatively simple geology, and reflect the relatively dense set of boreholes available for modelling. The method should be applied under a range of conditions to derive more general conclusions.     

The time-dependent motion of a floating elastic or rigid body in two dimensions

The linear time-dependent motion of a floating elastic or rigid body, subject to some initial displacement, which subsequently evolves freely is considered. The solution is derived by a Fourier transform and by the generalized eigenfunction method. Compared to other solutions, such as the Cummins method, the present solution requires neither time-stepping nor high-frequency calculations. A series of new identities for the frequency-domain problem are also presented. The Fourier transform solution allows an approximate solution to be calculated by an expansion over the complex resonances known as the singularity expansion method. Simple expressions for the singularity expansion method approximation are given. The method is illustrated with a series of numerical calculations.     

Monitoring cross-border trails using airborne digital multispectral imagery and interactive image analysis techniques

The objective of this study was to evaluate image-based procedures for monitoring cross-border foot trails in the US – Mexico border zone in eastern San Diego County using airborne remote sensing techniques. Specifically, digital multi-spectral and multi-temporal imagery from an airborne digital multi-spectral imaging system, digital image processing, and visual image analysis techniques were explored in the context of detecting and delineating new trail features and updating trail GIS layers. Three trail updating approaches: map-to-image (M-I) overlay, map and image-to-image (M/I-I) differencing, map and image-to-image (M/I-I) swiping and two types of spectral transform, PCA and NDVI, were tested and compared. The M-I overlay was found to be the most reliable trail updating approach. The optimal image enhancement method for the M-I overlay approach varied with vegetation structure. PCA imagery yielded better results than NDVI imagery in a highly disturbed area and NDVI imagery performed better in a densely vegetated area. The M/I-I swiping approach was useful for distinguishing misregistered extant trails from new trail features.     

基于高精度地层倾角线性反演的地震全体积拉平方法

地震数据的拉平分析是三维地震资料解释中的常见方法,一般通过拾取某一标准层位并将其校正到一个基准面的方式拉平地震数据,这样会导致距离标准层越远的层位拉平效果越差.近几年发展起来的三维地震数据全体积拉平方法是避免这种拉平缺陷的有效途径.本文提出一种基于高精度地层倾角的全体积自动拉平方法,该方法将一种自适应加权的向量滤波法获得高精度的局部地层倾角作为观测数据,根据拉平变换对每个数据样点产生的垂向偏移量与地层倾角之间的近似线性关系,建立最小二乘反演机制.在优化反演模型时,引入由梯度结构张量法得到的误差控制模板,抑制断层等复杂地质构造对拉平效果的影响;此外,增加模型参数的平滑度量项,避免拉平后的数据出现较大的畸变,并采用模型重新参数化方法,加快反演算法的收敛速度.三维合成数据和实际地震资料的试算结果表明这种三维全体积拉平方法是有效和可行的.     

On the prediction error variance of three common spatial interpolation schemes

Three forms of linear interpolation are routinely implemented in geographical information science, by interpolating between measurements made at the endpoints of a line, the vertices of a triangle, and the vertices of a rectangle (bilinear interpolation). Assuming the linear form of interpolation to be correct, we study the propagation of error when measurement error variances and covariances are known for the samples at the vertices of these geometric objects. We derive prediction error variances associated with interpolated values at generic points in the above objects, as well as expected (average) prediction error variances over random locations in these objects. We also place all the three variants of linear interpolation mentioned above within a geostatistical framework, and illustrate that they can be seen as particular cases of Universal Kriging (UK). We demonstrate that different definitions of measurement error in UK lead to different UK variants that, for particular expected profiles or surfaces (drift models), yield weights and predictions identical with the interpolation methods considered above, but produce fundamentally different (yet equally plausible from a pure data standpoint) prediction error variances.     

Implicit structural inversion of gravity data using linear programming,a validation study*

In this study, a regional scale gravity data set has been inverted to infer the structure (topography) of the top of the basement underlying sub‐horizontal strata. We apply our method to this real data set for further proof of concept, validation and benchmarking against results from an earlier forward modelling done elsewhere. Our aim is to carry out implicit structural inversion, i.e., to obtain a geologically reasonable model, without specifically solving for structure. The 2.5D volume of interest is parametrized with homogeneous horizontal prisms and a two‐lithology medium is assumed. A possible regional linear trend and a general floating reference are also inverted for. Using a gridded parametrization, linear programming is used to minimize the L₁ ‐norm of the data misfit, relative to a floating reference. Given a known density contrast between the lithologies, an inversion using linear programming has the intrinsic advantage that a relatively sharp image of the sub‐surface is retrieved instead of a smooth one. The model recovered is almost bi‐modal and its general features seem to be robust with respect to several parametrization scenarios investigated. The floating reference and a linear trend in the data were also retrieved simultaneously. The inversion results, indicating two depressions in the basement, are robust and agree with those obtained earlier based upon detailed 2D forward modelling using many narrow, near‐vertical prisms.