提高测井曲线纵向分辨率的方法

常规处理方法往往漏掉一些有价值的薄油气储层。利用频率域反槽积方法、贝叶斯反褶积方法、反褶积与光滑滤波组合（ＤＦＳＮ）方法和分辨率匹配与非线性拟合方法对测井曲线进行了提高纵向分辨率处理，测井曲线的纵向分辨率可提高１～４倍，一些结果得到了证实，为准确地定量评价薄储居提供了较可靠的参数     

基于点云空洞修复和TPS模型的数学形态学机载LIDAR滤波

在分析现有滤波算法的基础上,结合数据驱动和模型驱动式算子各自的优点,提出基于点云空洞修复和TPS变形模型的数学形态学机载LIDAR点云滤波,该方法首先提取和修复由水域造成的大面积点云空洞,采用多尺度形态学开算子作用于修复的数据,得到近似裸露地表面;然后利用2D空间的TPS变形模型,以近似地表面为基础,插值原始点云,根据插值与原始点云高程的差值大小去识别地面点和非地面点。通过定量分析,验证该方法不仅有较高的滤波精度,而且也能较好的保留裸露地表的细节特征,同时该方法有助于辅助人工处理,提高数据处理的质量。     

基于非稳态匹配的角度域叠前道集去调谐方法

地层调谐效应改变了不同偏移距处反射波的干涉模式,使得远偏移距数据振幅和频率信息发生畸变,将会降低AVO/AVA分析及反演的置信度和分辨率.文中首先给出了调谐作用下反射波合成记录,通过正演模拟说明了调谐对叠前道集的影响.考虑到不同地层厚度下调谐效应的差异性,借助局域Lamoureux窗实现地震数据的自适应分解,利用角度数据之间的差异,构建了角度域地震数据的非平稳匹配目标函数,形成了非平稳匹配去调谐方法,实现了叠前道集的振幅和波形拉伸校正.本文方法在实际应用中取得了较好的效果,能够有效地改善大角度地震数据品质,为储层预测与流体识别奠定了数据基础.     

基于两步弥散算子的共反射面元射线束成像方法研究

基于克希霍夫统一成像理论,对共反射面元射线束成像方法进行了深入分析.通过这种分析可以得到一种新的共反射面元射线束类成像实现方式：基于两步弥散算子的共反射面元射线束类成像方法.当仅实施第一步弥散时,即可获得信噪比大幅提高同时更为规则的叠前道集.二维、三维的理论和实际数据算例证实了上述观点.

     

Uncertainty assessment of ephemeral gully identification,characteristics and topographic threshold when using aerial photographs in agricultural settings

Manual digitizing on aerial photographs is still commonly used for characterizing gully erosion over large areas. Even when automated detection procedures are implemented, manual digitizing is frequently being resorted to in order to constitute reference datasets used for training and validation. In both cases, manual digitizing entails some subjective decisions on behalf of the operator, which introduces uncertainty into the resulting datasets. To assess the magnitude of this uncertainty, 11 experienced operators were asked to digitize and classify ephemeral gullies (EGs) on cropland following a standardized methodology. The resulting 11 datasets were compared in terms of number, type and location of EGs. Furthermore, for EGs located on a well‐defined runoff flow concentration axis, the slope versus contributing area topographic thresholds required for initiating gully channels were assessed using four thresholding methods, and compared across the 11 datasets. The operators identified 259 different EGs. However, the number (52–139) and sum total length (8.9–23.7 km) of EGs varied widely across operators. Only 34% of the EGs were digitized by more than half of the operators, and 7% were identified by all. Identification of EGs located on a well‐defined flow concentration axis proved least subjective. The longer the EG and the more fields the EG crossed, the larger the number of operators that were able to identify it. EGs were also most easily identified when located in sugar beet fields as compared to other crops. EG classification and topographic threshold lines were also found to be strongly operator‐dependent. Quantile regression appeared to be one of the most robust thresholding methods. Operator subjectivity when digitizing EGs on orthophotographs introduces uncertainty that should be taken into account in future remote sensing‐based studies of EG erosion whenever they rely, in part or in full, on manual photograph interpretation. Copyright © 2014 John Wiley & Sons, Ltd.     

Multiridge Euler deconvolution

Potential field interpretation can be carried out using multiscale methods. This class of methods analyses a multiscale data set, which is built by upward continuation of the original data to a number of altitudes conveniently chosen. Euler deconvolution can be cast into this multiscale environment by analysing data along ridges of potential fields, e.g., at those points along lines across scales where the field or its horizontal or vertical derivative respectively is zero. Previous work has shown that Euler equations are notably simplified along any of these ridges. Since a given anomaly may generate one or more ridges we describe in this paper how Euler deconvolution may be used to jointly invert data along all of them, so performing a multiridge Euler deconvolution. The method enjoys the stable and high‐resolution properties of multiscale methods, due to the composite upward continuation/vertical differentiation filter used. Such a physically‐based field transformation can have a positive effect on reducing both high‐wavenumber noise and interference or regional field effects. Multiridge Euler deconvolution can also be applied to the modulus of an analytic signal, gravity/magnetic gradient tensor components or Hilbert transform components. The advantages of using multiridge Euler deconvolution compared to single ridge Euler deconvolution include improved solution clustering, increased number of solutions, improvement of accuracy of the results obtainable from some types of ridges and greater ease in the selection of ridges to invert. The multiscale approach is particularly well suited to deal with non‐ideal sources. In these cases, our strategy is to find the optimal combination of upward continuation altitude range and data differentiation order, such that the field could be sensed as approximately homogeneous and then characterized by a structural index close to an integer value. This allows us to estimate depths related to the top or the centre of the structure.     

Generalized diffraction‐stack migration and filtering of coherent noise

We reformulate the equation of reverse‐time migration so that it can be interpreted as summing data along a series of hyperbola‐like curves, each one representing a different type of event such as a reflection or multiple. This is a generalization of the familiar diffraction‐stack migration algorithm where the migration image at a point is computed by the sum of trace amplitudes along an appropriate hyperbola‐like curve. Instead of summing along the curve associated with the primary reflection, the sum is over all scattering events and so this method is named generalized diffraction‐stack migration. This formulation leads to filters that can be applied to the generalized diffraction‐stack migration operator to mitigate coherent migration artefacts due to, e.g., crosstalk and aliasing. Results with both synthetic and field data show that generalized diffraction‐stack migration images have fewer artefacts than those computed by the standard reverse‐time migration algorithm. The main drawback is that generalized diffraction‐stack migration is much more memory intensive and I/O limited than the standard reverse‐time migration method.     

无人机核线影像的稀疏匹配与稠密匹配

无人机影像转化为水平核线影像后,能够有效地减少同名点的搜索空间。在此基础上,本文使用SIFT算子进行了稀疏匹配,并用BP算法进行了稠密匹配。结果表明：①SIFT算子获取的同名点比较少,但是计算方法简单,同名点空间坐标精确,适用于大范围获取简要的空间三维信息;②BP算法计算复杂度高,可以获取地物大量的同名点,适用于小范围的地物三维重建。总体而言,两者各有优缺点,在实际的应用中可互补。     

基于VB的测量坐标正算与反算研究

坐标正反算在测量实际工作中经常用到,本文介绍了测量坐标正算与反算的数学模型,采用VB开发语言,设计了测量坐标正算与反算程序,为测绘工作中坐标转换提供了思路。     

空间数据仓库的认知过程

从理论上探讨了空间数据仓库的认知过程．主要是十一个层次的世界模型通过十个算子转换的过程，对十一个层次世界及十个算子概念进行了较详细地讨论，并用代数系统对其进行了定义。