基于规则格网的DEM插值实验

通过地形建模,将6个地形隶属函数按照101×101,112×112,126×126,144×144,168×168,201×201,257×257 7种格网密度生成6种局部地形单元的规则格网DEM;使用反距离加权(IDW)等7种插值算法,将前6种格网密度下的DEM插值成257×257规格;从原始257×257DEM中随机抽取检查点计算残差,并对残差中误差进行分析。通过分组插值实验,运用控制变量法、方差分析等方法研究地貌类型、采样密度和插值算法对DEM插值精度的影响。     

平行轮廓线自适应插值的三维表面重构算法

针对稀疏点轮廓线在三维表面重构中存在的不够光滑、匹配易出错和难以交互等问题,提出了基于轮廓线插值的三维表面重构方法,并通过实验进行验证。文章分析了三维表面重构的形式化表达方法和准则,提出了基于虚跨距解决稀疏点轮廓线三维重构的基本思想;基于同步前进算法和虚跨距,设计了自适应等比例插值的同步前进三角面镶嵌算法;针对轮廓线插值违背三角表面重构拓扑准则的规律进行了分析,提出了相应的修正算法,保证了三角网的拓扑一致性。试验结果表明:本算法所构建的三角网表面更为光滑,各个三角面的大小更为均匀,具有更好的鲁棒性和交互性,对稀疏和稠密坐标点的轮廓线均能构建出更为理想的三维表面模型。     

抗差估计的多波束测深数据内插方法

针对传统多波束测深数据网格化内插方法抗差性不足的问题,该文将抗差估计理论应用于网格化内插方法中。基于不确定度指标建立了距离不确定度联合加权内插模型,在该模型的权函数中引入了等价权,同时给出了等价权设计方案及迭代初值确定原则,以增强模型的抗差性;介绍了一种节点邻域参考点合理选取方法,给出了改进内插法实现的基本步骤,最后通过实测数据对反距离加权法、联合加权法及改进方法内插水深质量进行比较。结果表明,改进方法内插结果要优于另外两种方法,可应用于多波束测深数据网格化处理中。     

Dasymetric Spatiotemporal Interpolation

This research applies the principles of dasymetric mapping to spatiotemporal interpolation by extending the spatial concepts of zone and area to their temporal analogs of interval and duration, respectively. An example application of dasymetric spatiotemporal interpolation using crime event data is presented. Results indicate that dasymetric spatiotemporal interpolation significantly improves the accuracy of estimates over areal or duration weighting. In addition, even when dasymetric interpolation in either the spatial or temporal dimension is relatively weak, combining dasymetric estimation in both space and time dimensions simultaneously has the potential to amplify the accuracy of the overall dasymetric estimation.     

Hybrid N-order Lagrangian Interpolation Eulerian?Lagrangian Method for Salinity Calculation

The Eulerian?Lagrangian method (ELM) has been used by many ocean models as the solution of the advection equation, but the numerical error caused by interpolation imposes restriction on its accuracy. In the present study, hybrid N-order Lagrangian interpolation ELM (LiELM) is put forward in which the N-order Lagrangian interpolation is used at first, then the lower order Lagrangian interpolation is applied in the points where the interpolation results are abnormally higher or lower. The calculation results of a step-shaped salinity advection model are analyzed, which show that higher order (N=3?8) LiELM can reduce the mean numerical error of salinity calculation, but the numerical oscillation error is still significant. Even number order LiELM makes larger numerical oscillation error than its adjacent odd number order LiELM. Hybrid N-order LiELM can remove numerical oscillation, and it significantly reduces the mean numerical error when N is even and the current is in fixed direction, while it makes less effect on mean numerical error when N is odd or the current direction changes periodically. Hybrid odd number order LiELM makes less mean numerical error than its adjacent even number order LiELM when the current is in the fixed direction, while the mean numerical error decreases as N increases when the current direction changes periodically, so odd number of N may be better for application. Among various types of Hybrid N-order LiELM, the scheme reducing N-order directly to 1st-order may be the optimal for synthetic selection of accuracy and computational efficiency.     

基于非均质地层参数表征的颗粒流模型构建

非均质地层参数的获取与表征方法是浅层地质体改造的重要依据。基于非均质地层的地质建模,表征具有随机不确定性和模糊不确定性的土体表观参数和物理力学参数。利用 Weibull 分布统计描述局部地层信息,并应用Diamond-Square分形插值方法,结合差分盒子维计算非均质地层分形维数,进一步合理演算整体地层参数。采用PFC3D数值模拟软件,结合马尔科夫链蒙特卡洛方法采样数据,对非均质地层块石颗粒及土颗粒进行建模。通过非均质地层模拟结果,发现其破坏过程与应力应变曲线要比均质土体复杂得多,由于块石大颗粒的含量、分布特性对土样力学性质影响显著。本文建模方法有助于非均质地层参数表征研究。     

Brittleness index predictions from Lower Barnett Shale well-log data applying an optimized data matching algorithm at various sampling densities

The capability of accurately predicting mineralogical brittleness index(BI)from basic suites of well logs is desir-able as it provides a useful indicator of the fracability of tight formations.Measuring mineralogical components in rocks is expensive and time consuming.However,the basic well log curves are not well correlated with BI so correlation-based,machine-learning methods are not able to derive highly accurate BI predictions using such data.A correlation-free,optimized data-matching algorithm is configured to predict BI on a supervised basis from well log and core data available from two published wells in the Lower Barnett Shale Formation(Texas).This transparent open box(TOB)algorithm matches data records by calculating the sum of squared errors be-tween their variables and selecting the best matches as those with the minimum squared errors.It then applies optimizers to adjust weights applied to individual variable errors to minimize the root mean square error(RMSE)between calculated and predicted(BI).The prediction accuracy achieved by TOB using just five well logs(Gr,pb,Ns,Rs,Dt)to predict BI is dependent on the density of data records sampled.At a sampling density of about one sample per 0.5 ft BI is predicted with RMSE～0.056 and R2～0.790.At a sampling density of about one sample per 0.1 ft BI is predicted with RMSE～0.008 and R2～0.995.Adding a stratigraphic height index as an additional(sixth)input variable method improves BI prediction accuracy to RMSE～0.003 and R2～0.999 for the two wells with only 1 record in 10,000 yielding a BI prediction error of＞±0.1.The model has the potential to be applied in an unsupervised basis to predict BI from basic well log data in surrounding wells lacking mineralogical measure-ments but with similar lithofacies and burial histories.The method could also be extended to predict elastic rock properties in and seismic attributes from wells and seismic data to improve the precision of brittleness index and fracability mapping spatially.     

梅江河流域年均降雨量空间插值方法研究

分析了降雨空间分布的影响要素及降雨空间插值的主要方法,并在此基础上,利用地理信息系统(GIS)的空间分析功能,综合考虑高程、坡度、坡向等影响降雨空间分布的多种要素,应用协克里金方法,探讨流域降雨量的空间插值问题及其应用,并与反距离权重法、克里金方法等单要素方法进行了降雨空间插值的效果比较.结果表明:应用协克里金方法,在综合考虑高程、坡度、坡向等影响要素的条件下,对降雨进行空间插值,其结果具有比较好的精度.     

矿体三维数据模型及品位插值方法研究

传统的储量计算方法采用几何平均法计算矿体的体积,用部分化验数据的平均品位代替矿块的整体品位,其计算精度很难满足需要.以三维矿体模型为基础计算矿体体积能够提高体积计算精度,以空间插值技术对矿体任意空间位置的品位进行估值,有助于提高品位估值的精度.提出了一种适于矿体储量计算的混合数据模型（线框块段模型）,详细论述了根据剖面数据生成矿体线框模型进而在线框模型基础上生成块段模型的方法.介绍了距离反比法和克立格法的原理,利用距离反比法结合某金矿体数据进行了储量估算,结果表明以线框块段模型为基础,采用空间插值技术计算矿体储量具有一定的实用价值.     

海底地形数据处理与DEM生成

分析了海底DEM生成经常会遇见的问题,介绍了坐标变换、接边处理、插值加密等数据预处理方法,最终生成海底规则格网DEM。