数据自相关多次波偏移成像

在常规偏移方法中一般都需要压制地震数据中的多次波,仅利用一次波信息成像,把自由表面反射的多次波视为噪声,但是在多次波中也包含着地下结构信息,应该将其充分利用到成像中来.事实上,已经有不少成像方法试图利用多次波信息,但是大部分方法都需要对多次波进行预测.本文提出了基于傅里叶有限差分偏移算子的数据自相关偏移方法.在这种偏移方法中,对含有一次波和多次波的地震数据,分别进行下行和上行延拓,然后直接利用常规的互相关成像条件成像.由于波场延拓采用了傅里叶有限差分算子,其计算效率高,能够很好地对复杂介质中的地震数据进行延拓.在数值试验中,使用了一个含散射点的三层模型和Marmousi模型.合成数据测试结果表明,这种方法可以对更大范围的地下构造成像,比常规的只利用一次波的傅里叶有限差分法照明度更好,并且在浅层可以提供更高的分辨率.我们提出的数据自相关策略易于实现且避免了繁杂的多次波预测,这对于复杂地下构造成像可能有着重大意义.     

再论地震数据偏移成像

利用地震波正向传播方程对属于波形线性反演问题近似求解方法的地震数据偏移成像进行重新推导,得到了适合散射地震数据的散射偏移成像方法和适合反射地震数据的反射偏移成像方法.以地震波传播的散射理论为出发点,首先根据描述一次散射波正向传播的线性方程研究建立散射地震数据的偏移成像方法理论;利用高频近似对产生散射波场的地下速度扰动函数的空间变化进行近似,推导出地下反射率函数,再由散射波传播方程推导出基于反射率函数的反射波传播方程,然后根据描述一次反射波正向传播的线性方程研究建立反射地震数据的偏移成像方法理论.本文指出和修正了Claerbout偏移成像方法中的不足,提出的地震数据偏移成像方法是对当前偏移成像方法理论的完善,使反射地震数据偏移成像具有了更坚实的数学物理理论基础,得到的偏移成像结果相位正确、位置准确、分辨率提高.     

从瞬变电磁测深数据到平面电磁波场数据的等效转换

中心回线源装置的瞬变电磁场是一种涡流场,在地下主要以扩散形式传播.平面电磁波在地下介质中传播时具有反射、折射等电磁学特性,为了借鉴平面电磁波场测深成熟的解释方法,文中开展了从瞬变电磁涡流场测深数据到平面波场数据的快速等效转换的研究.通过大量的理论模型的正演计算、曲线对比,误差统计,对瞬变电磁场、平面波场在地下传播机制特性的分析,建立了一种由瞬变电磁测深视电阻率数据向平面波场测深视电阻率数据转换的时间-频率等效对应关系.利用这一关系式,把扩散场数据快速转换成等效平面波场测深数据. 为进一步进行瞬变电磁拟平面波成像解释打下理论基础.     

Functional data clustering using K-means and random projection with applications to climatological data

In this paper, an efficient pattern recognition method for functional data is introduced. The proposed method works based on reproducing kernel Hilbert space (RKHS), random projection and K-means algorithm. First, the infinite dimensional data are projected onto RKHS, then they are projected iteratively onto some spaces with increasing dimension via random projection. K-means algorithm is applied to the projected data, and its solution is used to start K-means on the projected data in the next spaces. We implement the proposed algorithm on some simulated and climatological datasets and compare the obtained results with those achieved by K-means clustering using a single random projection and classical K-means. The proposed algorithm presents better results based on mean square distance (MSD) and Rand index as we have expected. Furthermore, a new kernel based on a wavelet function is used that gives a suitable reconstruction of curves, and the results are satisfactory.     

地下流体数据处理实用技巧

介绍了在前兆地下流体数据预处理时常见到的情况、处理依据及处理实用技巧,讨论了妥当处理观测数据的方法和依据以便提高数据的可靠性和可用性,为地震预报人员和数据使用者提供可靠的研究数据.     

网络化数据中心NetDC

IRIS等于1996年提出的网络化数据中心NetDC为用户访问各地震数据中心数据提供了一个统一的入口和界面，给用户访问数据带来极大的便利。NetDC是用C语言编写，建立在各地震数据中心的数据管理系统之上的软件系统，它很好地实现了与各数据中心现有的数据管理系统的集成，使得每个数据中心都能为用户提供更好的数据服务。本文简要介绍了网络化数据中心NetDC的概念、数据请求方法和格式，以及数据请求的处理过程等，并探讨了在我国地震系统的应用前景。     

三维叠后地震数据可纵向显示的处理方法

在三维地震资料处理中,为使叠后相邻测线的地震数据在绘图纸上能够纵向同时显示,达到节约绘图纸、降低成本的目的,本文提出了一种利用地震资料处理模块生成可纵向显示的三维地震绘图数据的处理方法.该方法通过进行叠后三维数据体的分离、延长测线数据的记录长度、修改线号和数据体舍并等处理工作,产生可纵向显示的绘图数据体.经实际绘图显示,不同测线的地震剖面由原来的横向单张显示,转为了纵向的多张绘图.因此,大大节约了纸张,降低了地震资料处理的成本.     

Worth of secondary data compared to piezometric data for the probabilistic assessment of radionuclide migration

A common approach for the performance assessment of radionuclide migration from a nuclear waste repository is by means of Monte-Carlo techniques. Multiple realizations of the parameters controlling radionuclide transport are generated and each one of these realizations is used in a numerical model to provide a transport prediction. The statistical analysis of all transport predictions is then used in performance assessment. In order to reduce the uncertainty on the predictions is necessary to incorporate as much information as possible in the generation of the parameter fields. In this regard, this paper focuses in the impact that conditioning the transmissivity fields to geophysical data and/or piezometric head data has on convective transport predictions in a two-dimensional heterogeneous formation. The Walker Lake data based is used to produce a heterogeneous log-transmissivity field with distinct non-Gaussian characteristics and a secondary variable that represents some geophysical attribute. In addition, the piezometric head field resulting from the steady-state solution of the groundwater flow equation is computed. These three reference fields are sampled to mimic a sampling campaign. Then, a series of Monte-Carlo exercises using different combinations of sampled data shows the relative worth of secondary data with respect to piezometric head data for transport predictions. The analysis shows that secondary data allows to reproduce the main spatial patterns of the reference transmissivity field and improves the mass transport predictions with respect to the case in which only transmissivity data is used. However, a few piezometric head measurements could be equally effective for the characterization of transport predictions.     

Worth of secondary data compared to piezometric data for the probabilistic assessment of radionuclide migration

A common approach for the performance assessment of radionuclide migration from a nuclear waste repository is by means of Monte-Carlo techniques. Multiple realizations of the parameters controlling radionuclide transport are generated and each one of these realizations is used in a numerical model to provide a transport prediction. The statistical analysis of all transport predictions is then used in performance assessment. In order to reduce the uncertainty on the predictions is necessary to incorporate as much information as possible in the generation of the parameter fields. In this regard, this paper focuses in the impact that conditioning the transmissivity fields to geophysical data and/or piezometric head data has on convective transport predictions in a two-dimensional heterogeneous formation. The Walker Lake data based is used to produce a heterogeneous log-transmissivity field with distinct non-Gaussian characteristics and a secondary variable that represents some geophysical attribute. In addition, the piezometric head field resulting from the steady-state solution of the groundwater flow equation is computed. These three reference fields are sampled to mimic a sampling campaign. Then, a series of Monte-Carlo exercises using different combinations of sampled data shows the relative worth of secondary data with respect to piezometric head data for transport predictions. The analysis shows that secondary data allows to reproduce the main spatial patterns of the reference transmissivity field and improves the mass transport predictions with respect to the case in which only transmissivity data is used. However, a few piezometric head measurements could be equally effective for the characterization of transport predictions.     

压缩感知技术在地震数据重建中的应用

地震资料室内处理过程要求野外采集的地震资料越多越好, 而地震数据远距离快速传输又要求野外地震数据量越少越好． 为解决这一矛盾, 将基于曲波变换与压缩感知的数据重建技术引入到地震资料处理中, 对实际的野外不完整数据进行压缩重建． 结果表明, 曲波变换相对于傅里叶变换在数据压缩采样方法中占有一定的优势． 但是, 在对实际资料进行处理时, 首先要对资料中的面波进行处理, 同时还要在一定曲波基元尺寸的情况下, 考虑缺失道数量的影响． 最终, 得到的重建数据图像纹理清晰、 连接自然, 从而验证了该方法的实用性和有效性．      

稀疏时频分解方法的研究与运用（英文）

Gabor变换和S变换是常用的时频分析工具。根据测不准原理,它们的时频分解结果无法在时间域和频率域同时具有很高的分辨率。为了提高非平稳信号时频分解结果的分辨率,本文提出瞬时频率分布函数（IFDF）并利用它表达非平稳信号。当非平稳信号时频成分的分布满足测不准原理对信号可分辨的要求时,瞬时频率分布函数的支集和短时Fourier变换的小波脊支集是同一个集合。利用IFDF的该特征,本文提出一种迭代算法（Sparse-STFT）实现了信号的稀疏时频分解。该算法在每次迭代过程中利用残留信号的短时Fourier变换结果的脊支集更新信号的时频成分,每次迭代得到的时频成分的叠加结果即为最终的稀疏时频分解结果。文中的数值实验证明了Sparse-STFT可以有效地提高非平稳信号时频分解结果的分辨率。最后,本文将该方法应用于地震数据面波的压制中,取得了理想的处理结果。     

Using data assimilation method to calibrate a heterogeneous conductivity field conditioning on transient flow test data

A data assimilation method is developed to calibrate a heterogeneous hydraulic conductivity field conditioning on transient pumping test data. The ensemble Kalman filter (EnKF) approach is used to update model parameters such as hydraulic conductivity and model variables such as hydraulic head using available data. A synthetical two-dimensional flow case is used to assess the capability of the EnKF method to calibrate a heterogeneous conductivity field by assimilating transient flow data from observation wells under different hydraulic boundary conditions. The study results indicate that the EnKF method will significantly improve the estimation of the hydraulic conductivity field by assimilating continuous hydraulic head measurements and the hydraulic boundary condition will significantly affect the simulation results. For our cases, after a few data assimilation steps, the assimilated conductivity field with four Neumann boundaries matches the real field well while the assimilated conductivity field with mixed Dirichlet and Neumann boundaries does not. We found in our cases that the ensemble size should be 300 or larger for the numerical simulation. The number and the locations of the observation wells will significantly affect the hydraulic conductivity field calibration.     

A new inverse method and application to ocean data

A new method is proposed to inverse normalization data of hidden variables in a dynamical system by embedding a time series in multidimensional spaces and applying a normalization analysis to the conditional probability density of points in the reconstructed phase spaces. The method is robust in the application to Lorenz system and 4-dimensional R?ssler system by testing quantitatively and qualitatively the correlation coefficient between inverse data and original data in time domain and in frequency domain, respectively. By applying the method to analyzing the South China Sea data, the normalization data of wind speed is extracted from the sea surface temperature time series.     

A data assimilation approach to discharge estimation from space

River discharge is currently monitored by a diminishing network of gauges, which provide a spatially incomplete picture of global discharges. This study assimilated water level information derived from a fused satellite Synthetic Aperture Radar (SAR) image and digital terrain model (DTM) with simulations from a coupled hydrological and hydrodynamic model to estimate discharge in an un‐gauged basin scenario. Assimilating water level measurements led to a 79% reduction in ensemble discharge uncertainty over the coupled hydrological hydrodynamic model alone. Measurement bias was evident, but the method still provided a means of improving estimates of discharge for high flows. The study demonstrates the potential of currently available synthetic aperture radar imagery to reduce discharge uncertainty in un‐gauged basins when combined with model simulations in a data assimilation framework, where sufficient topographic data are available. The work is timely because in the near future the launch of satellite radar missions will lead to a significant increase in the volume of data available for space‐borne discharge estimation. Copyright © 2009 John Wiley & Sons, Ltd.     

Sensitivity of goodness-of-fit statistics to rainfall data rounding off

An analysis based on the L-moments theory suggests of adopting the generalized Pareto distribution to interpret daily rainfall depths recorded by the rain-gauge network of the Hydrological Survey of the Sardinia Region. Nevertheless, a big problem, not yet completely resolved, arises in the estimation of a left-censoring threshold able to assure a good fitting of rainfall data with the generalized Pareto distribution. In order to detect an optimal threshold, keeping the largest possible number of data, we chose to apply a “failure-to-reject” method based on goodness-of-fit tests, as it was proposed by Choulakian and Stephens [Choulakian, V., Stephens, M.A., 2001. Goodness-of-fit tests for the generalized Pareto distribution. Technometrics 43, 478-484]. Unfortunately, the application of the test, using percentage points provided by Choulakian and Stephens (2001), did not succeed in detecting a useful threshold value in most analyzed time series. A deeper analysis revealed that these failures are mainly due to the presence of large quantities of rounding off values among sample data, affecting the distribution of goodness-of-fit statistics and leading to significant departures from percentage points expected for continuous random variables. A procedure based on Monte Carlo simulations is thus proposed to overcome these problems.     

Applicability of some infiltration formulae to rangeland infiltrometer data

Three algebraic infiltration equations (Kostiakov's, Horton's and Philip's) were examined to determine which one would best fit infiltrometer data collected from a variety of mostly semi-arid rangeland plant communities from both Australia and the United States. Approximately 1,100 infiltrometer plots were included in the analysis. Results indicated that, in every instance, Horton's equation best fit the infiltrometer data. Variability of “point” measures of short-term infiltration rates were never satisfactorily accounted for by using either Kostiakov's or Philip's equation. Though Horton's equation provided a best fit to the overall infiltration data, R² values indicated a potential usefulness of this equation only under the certain conditions that were sampled in several rangeland plant communities in the Northern Territory, Australia. The equation could not be considered consistently useful under conditions sampled on rangelands in the United States.     

A new approach to the integration of accelerometer data

A direct approach to the integration of accelerometer data is presented which depends only on the fact that the final velocity and possibly, the final displacement are known. Using a dynamic programming formulation, the accelerometer record is corrected to account for these cases. The formulation is presented for the general case where velocity, displacement, as well as accelerometer data, are available. A numerical example is included.     

Fitting a stochastic partial differential equation to aquifer data

The steady state two dimensional groundwater flow equation with constant transmissivities was studied by Whittle in 1954 as a stochastic Laplace equation. He showed that the correlation function consisted of a modified Bessel function of the second kind, order 1, multiplied by its argument. This paper uses this pioneering work of Whittle to fit an aquifer head field to unequally spaced observations by maximum likelihood. Observational error is also included in the model. Both the isotropic and anisotropic cases are considered. The fitted field is then calculated on a two dimensional grid together with its standard deviation. The method is closely related to the use of two-dimensional splines for fitting surfaces to irregularly spaced observations.     

IGS服务数据下载软件的研制

在GPS数据解算和研究中,IGS服务数据的下载工作繁琐,尤其网络不正常时,下载更是费时费力.针对地震系统需要,作者研制了在UNIX环境下运行专用于下载IGS服务数据的软件.用户只需要提供台站名,起始、终止时间便可下载需要的文件,下载的同时能够实时显示网络运行情况及下载进度,并对网络故障做出处理.作者采用套接口API编写FTP客户端程序,采用GTK 多线程并发及互斥等方法使软件更为健壮.