基于一种改进凸集投影方法的地震数据同时插值和去噪

基于稀疏反演的地震插值方法是一种重要的插值方法, 然而大多数这类方法只针对无噪声数据或者高信噪比数据插值.实际上, 地震数据含有各种噪声, 使得插值问题变得更加困难.凸集投影方法是一种高效的插值算法, 但是对于含噪声数据的插值效果不理想, 针对含噪声数据提出的加权凸集投影方法能够实现同时插值和去噪, 但是除了最小阈值需要认真选取外, 增加一个权重因子来实现去噪功能.本文由迭代阈值算法推导出加权凸集投影方法, 证明其是解无约束优化问题的一种方法, 加权因子可以看作拟合误差项的系数.本文还提出了一种改进的凸集投影方法, 与原始凸集投影方法相比该方法不需要增加任何计算量, 只要通过阈值的选择来进行插值和去噪.数值模拟证明了该算法的计算效率, 并且对含噪声数据能够实现较好的插值效果; 先插值后去噪的结果证明了同时去噪和插值算法的可靠性和稳定性.     

Hydraulic-Head Formulation for Density-Dependent Flow and Transport

Density-dependent flow and transport solutions for coastal saltwater intrusion investigations, analyses of fluid injection into deep brines, and studies of convective fingering and instabilities of denser fluids moving through less dense fluids typically formulate the groundwater flow equation in terms of pressure or equivalent freshwater head. A formulation of the flow equation in terms of hydraulic head is presented here as an alternative. The hydraulic-head formulation can facilitate adaptation of existing constant-density groundwater flow codes to include density-driven flow by avoiding the need to convert between freshwater head and hydraulic head within the code and by incorporating density-dependent terms as a compartmentalized “correction” to constant-density calculations already performed by the code. The hydraulic-head formulation also accommodates complexities such as unconfined groundwater flow and Newton-Raphson solution schemes more readily than the freshwater-head formulation. Simulation results are presented for four example problems solved using an implementation of the hydraulic-head formulation in MODFLOW.     

主动源与被动源地震数据插值及联合数据成像

本文提出了两种情况下主动源数据和被动源数据的插值方法,并研究了两种数据在偏移成像中的互补效果.基于互相关法被动源数据重构原理,本文提出了结合多域迭代去噪技术的重构方法.提出了两种时间域主动源和被动源数据的插值方法,分别是共炮点域能量匹配插值和共检波点域最小平方匹配插值.然后对获得的主动源和被动源联合地震数据进行叠前深度偏移成像.在被动源活跃度不是很高的地区进行被动源地震勘探时,少量的主动源地震数据可以有效控制和补充被动源数据的成像效果.在稀疏炮点的主动源勘探中,有效利用被动源的信息能够在成像中增加更多的细节信息,提高成像质量.     

岩石圈磁场数据网格化方法研究

选取最小曲率、克里格、改进Shepard、反距离加权和径向基函数等5种网格化数学模型,对小江断裂地磁总强度加密区岩石圈磁场数据进行数据网格化,采用均方根预测误差和插值数据残差均方根等评价指标对网格化结果进行评价,结果表明,克里格插值与反距离加权插值法的精度最高。进一步比较克里格插值与反距离加权插值法的网格化图形质量,结果显示克里格插值网格化过程中兼顾了数据的平滑性和各实测点与待估点之间的空间位置关系,避免了系统误差,得出克里格插值更适用于岩石圈磁场数据网格化的结论。     

地震数据Bregman整形迭代插值方法

人工地震方法由于受到野外观测系统和经济因素等的限制,采集的数据在空间方向总是不规则分布.但是,许多地震数据处理技术的应用（如：多次波衰减,偏移和时移地震）都基于空间规则分布条件下的地震数据体.因此,数据插值技术是地震数据处理流程中关键环节之一.失败的插值方法往往会引入虚假信息,给后续处理环节带来严重的影响.迭代插值方法是目前广泛应用的地震数据重建思路,但是常规的迭代插值方法往往很难保证插值精度,并且迭代收敛速度较慢,尤其存在随机噪声的情况下,插值地震道与原始地震道之间存在较大的信噪比差异.因此开发快速的、有效的迭代数据插值方法具有重要的工业价值.本文将地震数据插值归纳为数学基追踪问题,在压缩感知理论框架下,提出新的非线性Bregman整形迭代算法来求解约束最小化问题,同时在迭代过程中提出两种匹配的迭代控制准则,通过有效的稀疏变换对缺失数据进行重建.通过理论模型和实际数据测试本文方法,并且与常规迭代插值算法进行比较,结果表明Bregman整形迭代插值方法能够更加有效地恢复含有随机噪声的缺失地震信息.

     

Analysis of radial basis function interpolation approach

The radial basis function （RBF） interpolation approach proposed by Freedman is used to solve inverse problems encountered in well-logging and other petrophysical issues. The approach is to predict petrophysical properties in the laboratory on the basis of physical rock datasets, which include the formation factor, viscosity, permeability, and molecular composition. However, this approach does not consider the effect of spatial distribution of the calibration data on the interpolation result. This study proposes a new RBF interpolation approach based on the Freedman＇s RBF interpolation approach, by which the unit basis functions are uniformly populated in the space domain. The inverse results of the two approaches are comparatively analyzed by using our datasets. We determine that although the interpolation effects of the two approaches are equivalent, the new approach is more flexible and beneficial for reducing the number of basis functions when the database is large, resulting in simplification of the interpolation function expression. However, the predicted results of the central data are not sufficiently satisfied when the data clusters are far apart.     

地震同震形变场失相干恢复方法研究

在合成孔径雷达干涉测量中,干涉相干是局部干涉条纹质量的重要评价标准,但是由于各种失相干源的影响,会出现失相干现象,而这可能会丢失重要的形变数据。本文首先利用ENVISAT ASAR的两景影像获得了2009年意大利拉奎拉M_W6.3地震的带有失相干现象的同震形变场,对数据预处理后再使用移动窗口克里金法对拉奎拉地震的同震形变场进行插值处理,恢复了形变图的失相干区域,之后进一步通过正反演获得同震形变场,并与上一步的插值结果进行对比分析。结果显示,插值结果与反演结果在极震区吻合得很好,验证了移动窗口克里金插值法恢复失相干的可靠性,说明插值法可以成为实现地震同震形变场失相干恢复的一种途径。      

一种边缘保持的地震数据插值方法

在地震数据处理中,地震数据插值方法常常用来解决地震数据空间采样率低和不规则的问题.本文提出了一种基于边缘保持滤波器的地震数据插值方法.在该方法中,对于一个1D信号,逐点滑动一个处理窗口,将信号分成多个信号片段.对于某一个待恢复的缺失采样点,存在多个包含(或邻近)该采样点的信号片段可以用来估计它.采用多项式来拟合这些信号片段,并选择拟合误差最小的信号片段估计此缺失采样点,达到边缘保持的目的.对于2D地震信号,先沿不同方向扫描抽取1D信号,然后采用上述1D边缘保持插值算法分别进行处理,得到沿不同方向的插值结果.对于任一待插值采样点,选取对应拟合误差最小的方向的插值结果作为最后输出的2D数据的插值结果.理论模型和实际资料的处理结果表明,所提方法具有保边缘、抗假频及能够进行不规则数据重建等特点,既能有效的实现不规则地震数据的重建,又能很好的保持原有数据的边缘特征.     

地震信号内插与噪音剔除（一）

信号内插与噪音剔除是相互有着有机联系的两个方面。首先讨论地震信号的内插方法,提出了一种“最佳内插算子”,从而导出了另一种“检噪算子”,它可以有效地把干扰波从记录中识别出来,并将其“剔除”出去。进而讨论空间域数据的内插问题--道内插问题。在一定的条件下,空间域数据内插与时间域内插是完全等效的,因此道内插也可以采用与时间域相同的“最佳内插算子”。地震记录上的干扰波其绝大部分表现为空间域的脉冲式干扰,所以用“检噪算子”在空间域作褶积,便能识别干扰,从而对干扰加以剔除。本文解决了一种迭代的剔除方法,取得了好的效果。     

Bilinear surface smoothing for spatial interpolation with optional incorporation of an explanatory variable. Part 2: Application to synthesized and rainfall data

ABSTRACT

The non-parametric mathematical framework of bilinear surface smoothing (BSS) methodology provides flexible means for spatial (two dimensional) interpolation of variables. As presented in a companion paper, interpolation is accomplished by means of fitting consecutive bilinear surface into a regression model with known break points and adjustable smoothing terms defined by means of angles formed by those bilinear surface. Additionally, the second version of the methodology (BSSE) incorporates, in an objective manner, the influence of an explanatory variable available at a considerably denser dataset. In the present study, both versions are explored and illustrated using both synthesized and real world (hydrological) data, and practical aspects of their application are discussed. Also, comparison and validation against the results of commonly used spatial interpolation methods (inverse distance weighted, spline, ordinary kriging and ordinary cokriging) are performed in the context of the real world application. In every case, the method’s efficiency to perform interpolation between data points that are interrelated in a complicated manner was confirmed. Especially during the validation procedure presented in the real world case study, BSSE yielded very good results, outperforming those of the other interpolation methods. Given the simplicity of the approach, the proposed mathematical framework’s overall performance is quite satisfactory, indicating its applicability for diverse tasks of scientific and engineering hydrology and beyond.

Editor Z. W. Kundzewicz; Associate editor A. Carsteanu     

A new approach to spatial data interpolation using higher-order statistics

Interpolation techniques for spatial data have been applied frequently in various fields of geosciences. Although most conventional interpolation methods assume that it is sufficient to use first- and second-order statistics to characterize random fields, researchers have now realized that these methods cannot always provide reliable interpolation results, since geological and environmental phenomena tend to be very complex, presenting non-Gaussian distribution and/or non-linear inter-variable relationship. This paper proposes a new approach to the interpolation of spatial data, which can be applied with great flexibility. Suitable cross-variable higher-order spatial statistics are developed to measure the spatial relationship between the random variable at an unsampled location and those in its neighbourhood. Given the computed cross-variable higher-order spatial statistics, the conditional probability density function is approximated via polynomial expansions, which is then utilized to determine the interpolated value at the unsampled location as an expectation. In addition, the uncertainty associated with the interpolation is quantified by constructing prediction intervals of interpolated values. The proposed method is applied to a mineral deposit dataset, and the results demonstrate that it outperforms kriging methods in uncertainty quantification. The introduction of the cross-variable higher-order spatial statistics noticeably improves the quality of the interpolation since it enriches the information that can be extracted from the observed data, and this benefit is substantial when working with data that are sparse or have non-trivial dependence structures.     

Modelling the Velocity Field in a Regular Grid in the Area of Poland on the Basis of the Velocities of European Permanent Stations

The paper presents the results of testing the various methods of permanent stations’ velocity residua interpolation in a regular grid, which constitutes a continuous model of the velocity field in the territory of Poland. Three packages of software were used in the research from the point of view of interpolation: GMT (The Generic Mapping Tools), Surfer and ArcGIS. The following methods were tested in the softwares: the Nearest Neighbor, Triangulation (TIN), Spline Interpolation, Surface, Inverse Distance to a Power, Minimum Curvature and Kriging. The presented research used the absolute velocities’ values expressed in the ITRF2005 reference frame and the intraplate velocities related to the NUVEL model of over 300 permanent reference stations of the EPN and ASG-EUPOS networks covering the area of Europe. Interpolation for the area of Poland was done using data from the whole area of Europe to make the results at the borders of the interpolation area reliable. As a result of this research, an optimum method of such data interpolation was developed. All the mentioned methods were tested for being local or global, for the possibility to compute errors of the interpolated values, for explicitness and fidelity of the interpolation functions or the smoothing mode. In the authors’ opinion, the best data interpolation method is Kriging with the linear semivariogram model run in the Surfer programme because it allows for the computation of errors in the interpolated values and it is a global method (it distorts the results in the least way). Alternately, it is acceptable to use the Minimum Curvature method. Empirical analysis of the interpolation results obtained by means of the two methods showed that the results are identical. The tests were conducted using the intraplate velocities of the European sites. Statistics in the form of computing the minimum, maximum and mean values of the interpolated North and East components of the velocity residuum were prepared for all the tested methods, and each of the resulting continuous velocity fields was visualized by means of the GMT programme. The interpolated components of the velocities and their residua are presented in the form of tables and bar diagrams.     

三维地震数据离散光滑插值的共轭梯度法

针对三维地震数据插值,提出采用Laplacian算子进行光滑约束的插值方法,并借鉴Mallet研究的离散光滑插值思路,采用预条件共轭梯度法,直接生成网格节点上的值,从而回避寻求满足插值方程的函数. 为了实现其中Laplacian算子的快速求逆,文中引入Claerbout螺旋坐标系谱因式分解理论. 在螺旋坐标系下,Laplacian算子的表示矩阵具有Toeplitz结构,其快速求逆可由谱法LU分解实现. 基于二维离散光滑插值,文中还给出共轭梯度法与NMO相结合的沿时间切片逐层处理的离散光滑插值流程. 最后,应用该方法对模型数据和实际三维地震数据进行了处理.     

三维地震数据离散光滑插值的共轭梯度法

针对三维地震数据插值,提出采用Laplacian算子进行光滑约束的插值方法,并借鉴Mallet研究的离散光滑插值思路,采用预条件共轭梯度法,直接生成网格节点上的值,从而回避寻求满足插值方程的函数. 为了实现其中Laplacian算子的快速求逆,文中引入Claerbout螺旋坐标系谱因式分解理论. 在螺旋坐标系下,Laplacian算子的表示矩阵具有Toeplitz结构,其快速求逆可由谱法LU分解实现. 基于二维离散光滑插值,文中还给出共轭梯度法与NMO相结合的沿时间切片逐层处理的离散光滑插值流程. 最后,应用该方法对模型数据和实际三维地震数据进行了处理.     

Simultaneous Denoising and Interpolation of Seismic Data via the Deep Learning Method

Utilizing data from controlled seismic sources to image the subsurface structures and invert the physical properties of the subsurface media is a major effort in exploration geophysics. Dense seismic records with high signal-to-noise ratio (SNR) and high fidelity helps in producing high quality imaging results. Therefore, seismic data denoising and missing traces reconstruction are significant for seismic data processing. Traditional denoising and interpolation methods rarely occasioned rely on noise level estimations, thus requiring heavy manual work to deal with records and the selection of optimal parameters. We propose a simultaneous denoising and interpolation method based on deep learning. For noisy records with missing traces, we adopt an iterative alternating optimization strategy and separate the objective function of the data restoring problem into two sub-problems. The seismic records can be reconstructed by solving a least-square problem and applying a set of pre-trained denoising models alternatively and iteratively.We demonstrate this method with synthetic and field data.     

Discrete stability of the Differential System Method evaluated with geostatistical techniques

The Differential System Method (DSM) permits identification of the physical parameters of finite-difference groundwater flow models in a confined aquifer when piezometric head and source terms are known at each point of the finite-difference lattice for at least two independent flow situations for which the hydraulic gradients are not parallel. Since piezometric head data are usually few and sparse, interpolation of the measured data onto a regular grid can be performed with geostatistical techniques. We apply kriging to the sparse data of a synthetic aquifer to evaluate the stability of the DSM with respect to uncorrelated measurement errors and interpolation errors. The numerical results show that the DSM is stable.     

Discrete stability of the Differential System Method evaluated with geostatistical techniques

The Differential System Method (DSM) permits identification of the physical parameters of finite-difference groundwater flow models in a confined aquifer when piezometric head and source terms are known at each point of the finite-difference lattice for at least two independent flow situations for which the hydraulic gradients are not parallel. Since piezometric head data are usually few and sparse, interpolation of the measured data onto a regular grid can be performed with geostatistical techniques. We apply kriging to the sparse data of a synthetic aquifer to evaluate the stability of the DSM with respect to uncorrelated measurement errors and interpolation errors. The numerical results show that the DSM is stable.     

分形插值法在地电场干扰数据处理中的应用

当地电场和地电阻率同场地观测时,地电场观测会受到地电阻率观测的供电干扰,这类干扰时间短、干扰形态和出现时间固定,影响了地电场观测数据的正常变化形态,给数据分析和地震科学研究造成困难.为解决这一干扰问题,本文在比较分形插值方法与传统插值方法优劣的基础上,采用分形插值方法对受干扰的地电场观测数据进行重建,以提高信号重建的精...     

Cokriging for enhanced spatial interpolation of rainfall in two Australian catchments

Rainfall data in continuous space provide an essential input for most hydrological and water resources planning studies. Spatial distribution of rainfall is usually estimated using ground‐based point rainfall data from sparsely positioned rain‐gauge stations in a rain‐gauge network. Kriging has become a widely used interpolation method to estimate the spatial distribution of climate variables including rainfall. The objective of this study is to evaluate three geostatistical (ordinary kriging [OK], ordinary cokriging [OCK], kriging with an external drift [KED]), and two deterministic (inverse distance weighting, radial basis function) interpolation methods for enhanced spatial interpolation of monthly rainfall in the Middle Yarra River catchment and the Ovens River catchment in Victoria, Australia. Historical rainfall records from existing rain‐gauge stations of the catchments during 1980–2012 period are used for the analysis. A digital elevation model of each catchment is used as the supplementary information in addition to rainfall for the OCK and kriging with an external drift methods. The prediction performance of the adopted interpolation methods is assessed through cross‐validation. Results indicate that the geostatistical methods outperform the deterministic methods for spatial interpolation of rainfall. Results also indicate that among the geostatistical methods, the OCK method is found to be the best interpolator for estimating spatial rainfall distribution in both the catchments with the lowest prediction error between the observed and estimated monthly rainfall. Thus, this study demonstrates that the use of elevation as an auxiliary variable in addition to rainfall data in the geostatistical framework can significantly enhance the estimation of rainfall over a catchment.