频率波数域预测和减去法压制多次波

研究了压制与海面有关的地震勘探多次反射波的两种方法。两种方法均以波动方程为基础。当地震子波已知时，可在频率波数域中通过子波反褶积方法消除与海水表面有关的多次波；在未知子波的情况下，可依据消除后波场能量最小为准则，求出一滤波因子，用它去标定估计的多次波，然后从地震记录中减去预测的多次波，达到压制它们的目的。为了降低维数，x、y的二维空间坐标用柱坐标的径向距离代替，其波数通过汉克尔变换求得。     

运动荷载附近有限层厚软土地基的振动研究

基于Biot多孔弹性介质的波动理论,研究了运动荷载附近软土地基的振动问题。假设一条形均布荷载作用在地基表面,则该模型可视为平面应变问题进行分析。通过引入4个势函数和Helmholtz原理,并利用Fourier变换及逆变换技术,获得了运动荷载作用下软土地基的应力、位移和孔隙水压力的解答。利用离散Fourier逆变换得到数值计算结果,分析了荷载速度,频率以及软土的渗透系数及多孔弹性参数对运动荷载作用下地表竖向位移及土体中任一点孔隙水压力分布的影响。     

一种加权抛物Radon变换地震波场重建方法

基于带限正反抛物Radon变换的基本原理和计算方法,笔者给出了一种加权抛物Radon变换（PRT）算法,该方法可应用于叠前地震波场恢复和重建,如缺失偏移距地震数据的插值外推和反假频重采样,其权系数完全由数据域来确定,也就是说在迭代计算过程中,该权系数也在变化。对于均匀采样的地震数据,笔者给出了近偏移距外推和波场的反假频重采样的计算结果,试算表明笔者给出的加权PRT方法可有效地应用于地震波场重建中,并具有算法稳健、计算精度高和效率高的特点。     

多层横观各向同性地基轴对称固结的传递矩阵解

从横观各向同性地基轴对称Biot固结的基本方程出发,通过关于t的Laplace变换和关于r的Hankel变换,得到关于z的一阶常微分方程组。然后,对变换域内的基本未知量进行线性化处理,建立了变换域内的基本状态变量在z = 0处和任意深度处z的显式关系。利用传递矩阵法,结合层间连续性条件和边界条件,得到了多层横观各向同性地基的Biot固结轴对称问题的解答。该解答能避免随着层数增加而需要求解大型方程组的困难,明显地提高了计算效率。     

考虑热-水-力耦合效应多孔弹性地基的动力响应

通过对Biot波动方程的修正,得到考虑热-水-力学耦合效应的多孔弹性介质动力响应的控制方程,研究了简谐均布荷载作用下地基土体的热-水-力耦合动力响应问题。利用Fourier变换技术,得到地基中的应力、位移和孔隙水压力积分形式的解答。利用Fourier逆变换得到数值结果,分析了热-水-力学耦合条件下地基土体中温度增量、应力、位移和孔隙水压力响应的分布,并讨论了热源输入的影响, 结果表明：应力、位移和孔隙水压力随 的增大而有一定的减小。     

广义S变换时频域滤波在MT数据处理中的应用

短时傅里叶变换是建立在稳态信号基础之上,它仅能提供信号的频域信息,对信号的时间分辨能力差。这影响了它在大地电磁测深数据处理中的应用效果。S变换是一种优于短时傅立叶变换的时频分析方法,能够提供信号时-频域信息。利用S变换对大地电磁测深数据进行时频分析,有助于实现大地电磁测深数据噪声的时频-域滤波,从而提高大地电、磁分量数据的频谱分析精度。从广义S变换理论出发,分析了各类波形噪声的时-频域特征及其对大地电磁测深数据的影响。针对大地电磁测深数据处理特点,利用广义S变换得到时频谱,采用时频比值和门槛值方法,研究适合压制电磁噪声的时频滤波器和滤波方法。对实际大地电磁测深数据的处理结果表明这个方法提高了阻抗张量的估算质量,验证了该方法的有效性。     

矩形运动荷载作用下软土地基的三维振动研究

通过引入势函数,并利用Helmholtz原理和Fourier变换技术,研究了运动荷载作用下有限层厚软土地基的振动,考虑了矩形分布荷载作用下振动的三维特性,使得分析更符合工程实际,给出了运动荷载作用下饱和黏弹性地基三维振动的应力、位移和孔隙水压力响应的积分形式解答。利用Fourier数值逆变换进行算例分析,讨论了荷载速度对位移及孔隙水压力分布的影响。结果表明,位移幅值随荷载速度的增加而增大,荷载不同位置处孔隙水压力的分布有很大差异。     

Fast sequential indicator simulation: Beyond reproduction of indicator variograms

Sequential Indicator Simulation (SIS), although widely used, is relatively slow, and requires tedious inference of a large number of indicator variogram models. SIS is designed only to estimate class proportions and to reproduce indicator variogram models; the statistics of the continuous attribute being simulated,z-histogram and variogram, may be poorly reproduced. Several implementations of the SIS algorithm are proposed resulting in better reproduction of statistics yet with better CPU performance.     

应用改进的KL变换提高叠后地震记录信噪比

常规KL变换叠后去噪仅能加强水平方向同相轴,对倾斜或弯曲同相轴处理效果较差,对能量较强的相干噪声去除效果不佳,且在处理大量数据时计算成本很高,在实际生产中难以广泛应用。本文针对常规KL变换的缺陷分别进行了三点改进:使用倾角扫描叠加KL变换,可在有效压制地震剖面噪声的同时,较好保持倾斜同相轴和弯曲同相轴;使用本文提出的时空变-倾角KL变换能去除能量较强的相干噪声;使用数据分块技术可减少运算量。模型测试和实际资料处理效果表明:改进后的KL变换适应性强,去噪效果大大改善。      

On the Optimal Estimation of the Cumulative Distribution Function in Presence of Spatial Dependence

For spatial analysis of data, the optimal prediction of the cumulative distribution function (c.d.f.) is made difficult by the presence of spatial dependence, particularly when data are clustered. Several methods have been proposed in the literature to correct for the redundancy which is present when the observed values are correlated. Most of them are based on an unequal weighting of the values for the computation of the c.d.f. The methodology is extended by defining weights which depend both on the spatial correlation and on the threshold value. If a parametric multivariate distribution is assumed, optimal weights can be obtained conditionally to the choice of this distribution, and a nonparametric indicator estimate of the c.d.f. can be defined a posteriori. The method can be extended in order to provide kernel smoothed estimates of the c.d.f., by replacing the indicator functions by other functions taken here as complementary differentiable c.d.f., whose parameters are defined in order to preserve the unbiasedness and minimum variance property of the estimate, conditionally to the choice of the a priori distribution. A simulation study is conducted in order to emphasize the properties of the various methods. It is found that, even for bounded and positively skewed distributions, a simplified version of the method seems to perform quite well without significant loss of efficiency, and provides smooth estimates of the c.d.f. which can be differentiated in order to obtain density estimates. These estimates do not depend highly on the choice of the parameters for the a priori distribution, and yet they respect major characteristics of the a priori c.d.f., such as lower and upper bounds, which are mainly dictated by the nature of the variable under study. A practical case study is conducted on cobalt concentrations for various rock types in the Swiss Jura, illustrating the applicability of the methodology for real datasets.     

Aftershock sequences in the mid-ocean ridge environment: an analysis using hydroacoustic data

Hydroacoustic data from autonomous arrays and the U.S. Navy's Sound Surveillance System (SOSUS) provide an opportunity to examine the temporal and spatial properties of seismicity along portions of the slow-spreading Mid-Atlantic Ridge (MAR), intermediate-spreading Juan de Fuca Ridge (JdFR) and fast-spreading East Pacific Rise (EPR). Aftershock and foreshock events are selected from the hydroacoustic earthquake catalog using single-link cluster (SLC) analysis, with a combined space–time metric. In the regions examined, hydroacoustic data improve the completeness level of the earthquake catalog by 1.5–2.0 orders of magnitude, allowing the decay constant, p, of the modified Omori law (MOL) to be determined for individual sequences. A non-parametric goodness-of-fit test indicates six of the seven sequences examined are described well by a MOL model. The p-values obtained for individual ridge and transform sequences using hydroacoustic data are larger than that previously estimated from the analysis of a stacked sequence generated from teleseismic data. For three sequences along the Siqueiros, Discovery and western Blanco Transforms, p-values are estimated to be 0.94–1.29. The spatial distribution of aftershocks suggests that the mainshock rupture is constrained by intra-transform spreading centers at these locations. An aftershock sequence following a 7.1M_s thrust event near the northern edge of the Easter Microplate exhibits p=1.02±0.11. Within the sequence, aftershocks are located to the north of a large topographic ridge, which may represent the surface expression of the shallow-dipping fault that ruptured during the mainshock. Two aftershock sequences near 24°25′N and 16°35′N on the MAR exhibit higher p-values, 1.74±0.23 and 2.37±1.65, although the latter estimate is not well constrained because of the small number of aftershocks. Larger p-values along the ridge crest might reflect a hotter thermal regime in this setting. Additional monitoring, however, will be needed to determine if p-value differences between the ridge and transform sequences are robust. A 1999 sequence on the Endeavour segment of the JdFR, which has been correlated with changes in the hydrothermal system, is described poorly by the MOL model. The failure of the MOL model, the anomalously large number of earthquakes within the sequence and absence of a clearly dominant mainshock are inconsistent with aftershock activity and the simple tectonic origin that has been proposed previously for this sequence.     

Predictive GIS-Based Model of Rockfall Activity in Mountain Cliffs

Rockfall susceptibility has been analysed in mountain cliffs of the Cantabrian Range, North Spain. The main aim of this analysis has been to build a predictive model of rockfall activity from a low number of environmental and geological variables. The rockfall activity has been quantified in a GIS. The cartographic information used shows the spatial distribution of all the recent talus screes as well as their associated source areas in the rock-slopes. The area relation At/Ar (recent talus scree polygon/source basins) in the rock slopes has been used as the rockfall activity indicator. This relation has been validated in 50 pilot rock-slopes and compared with the relation number of recent rock fragments/source basin, obtained from field work. The environmental factors causing rockfall depend on the rock slope situation, and these are: altitude and sun radiation on the rock cliff. The geological factors considered are: lithology, relative position of the main discontinuities with respect to the topographic surface and two morphologic parameters: the roughness and slope gradient. A logistic regression analysis has been applied to a population of 442 limestone and quartzite rock cliffs. The dependent variable is the rockfall activity indicator, which allows the definition of two classes of rock cliff units: low and high activity. The independent variables are altitude, sun radiation (equinox radiation, summer solstice radiation, winter solstice radiation), slope roughness, slope gradient,anisotropy and lithology. Results suggest that it is possible tobuild a valid cartographic predictive model for rockfall activity in mountain rock cliffs from a limited number of easily obtainable variables. The method is especially applicable in massive rock slopes or in regions with uniform rock mass characteristics.     

Quantitative analysis of anomalous sources and geochemical signatures in the Walker Lake quadrangle of Nevada and California

A new quantitative transportation model for geochemical elements in stream sediment samples is proposed and used to estimate mineral sources of anomalous concentrations and to adjust both location and concentration for the transportation effect. Subsequent to transportation analysis, data on 14 geochemical elements in the Walker Lake quadrangle are analyzed for geochemical signatures for three types of mineralization environments: epithermal gold-silver deposits, porphyry-skarn copper deposits, and a mixed occurrence. Geochemical signatures are analyzed using a multidimensional scaling method referred to as the p_ijk model, which is a factor-like model generalized to include nonlinear associations, as well as linear, among triple variables (elements). A result of this analysis is that Pb is the best indicator for epithermal gold-silver mineralization, Mo and Zn for mixed mineralization, and Cu for porphyry-skarn mineralization.     

基于多模式分离的S变换与小波变换提取面波频散曲线对比分析

为了对比分析小波变换和S变换在计算瑞雷波频散曲线时的精度,通过设计两层介质和六层介质模型,采用交错网格高阶有限差分方法对两个模型进行正演模拟,利用F-K法分离出时间—空间域不同模态的瑞雷波记录,再采用小波变换和S变换计算出不同模态的瑞雷波频散曲线,分别与理论值进行对比分析,最后利用两种方法对实际资料进行处理。结果表明:利用小波变换计算得到的频散曲线在低频段误差相对较大,S变换计算得到的频散曲线精度总体高于小波变换。     

Analytic signals for multivariate data

Hilbert transformation is a standard tool in univariate signal-processing. It leaves the information content unaffected for, apart from a change of sign, the iterated transform reproduces the original data. As we know from Fourier transforms, such alternative representations of the same information, notwithstanding their theoretical equivalence with the data, can assist powerfully with extracting and interpreting that information. Although the extension to multivariate data is not so obvious for Hilbert as for Fourier transforms, Nabighian gave a treatment of the bivariale situation in 1984. Fueter, some 50 years earlier, had worked on an analogous extension problem, seeking to generalize complex function theory. On comparing these two developments we learn that, although Nabighian's transforms fit naturally into Fueter's theory, they are only one among many alternative possibilities. This paper presents a general theory, of higher dimensional Hilbert transforms and analytic signals, applicable to data of all dimensions less than eight. The change-of-field-direction fillers used in geophysical data processing are shown to arise as special situations.     

τ-p变换的两种实现方法

τ-p变换是一种经典的函数投影变换方法,在地震数据处理中有着广泛的应用,可以根据地震数据射线参数的差异,实现信噪分离、地震道插值、平面波分解等。但是,由于τ-p变换的精度和分辨率受到相应数学反问题的限制,在保证较高变换精度的前提下,τ-p变换计算速度的加快,以及相应滤波算子的设计等都值得研究。本文提出了基于径向道变换的τ-p变换方法和基于斜率分解的τ-p变换方法,阐明了两种方法的基本原理并与常规τ-p变换方法进行对比。基于径向道变换的τ-p变换方法利用快速Fourier变换和径向道变换,能有效减少τ-p变换的耗时;基于斜率分解的τ-p变换方法运用斜率分解算法,能进行高分辨率τ-p变换,并且提供冗余变换维度（τ-x-p域）,使滤波器算子的设计更灵活。数值实验结果表明,两种实现方法分别在计算速度和重构精度上优于传统τ-p变换算法。通过面波噪声压制的实际数据处理和比较,证明本文提出的两种τ-p变换实现方法可以为实际处理提供更加有效和灵活的实施方案。     

基于脊波域的低信噪比地震资料处理技术

在小波分析的基础上发展起来的脊波变换方法,能更好地处理含线状变化特征的信号。针对地震资料中的同相轴信息,尝试利用脊波变换方法对其进行处理,提高剖面信噪比,突出同相轴信息。在对某工区实际地震资料的处理中发现,处理后的地震剖面同相轴品质及连续性有了明显改善,信噪比增强,分辨率相应提高,体现了该方法相对常规小波分析方法的优越性。      

Geochemical markers based on concentration ratios of PAH in oils and oil-polluted areas

The paper presents data on the possibility of using the proportions of concentrations of polycyclic aromatic hydrocarbons (PAH) as indicators of the pollution sources with oil hydrocarbons. Approaches are suggested to estimate the efficiency of these indicator ratios, and the efficiency of the currently used ratios is evaluated. Multivariate data analysis is applied to demonstrate how the ratios depend on the degree of transformation of the petroleum products. The data presented in the paper are utilized to more accurately determine the indicator ratios of oil-bearing samples with regard for the contamination age.     

Undiscovered Petroleum Accumulation Mapping Using Model-Based Stochastic Simulation

Stochastic simulation has been proven to be a useful tool for revealing uncertainties in petroleum exploration and exploitation. The application to petroleum resource assessment would result in predicted potential accumulations with geographic locations, a desirable feature for improving both resource management and exploration efficiency. The associated uncertainties with the prediction provide information useful for exploration risk analysis. This attempt has been encumbered by two typical technical difficulties: biased observation data and lack of information with respect to the undiscovered accumulation locations. In this paper we propose a model-based simulation approach, in which models are used to perform unbiased parameter estimation from biased data and to facilitate the location of undiscovered petroleum accumulations based on reasoning of available geological and geophysical observations. The Fourier transform algorithm is chosen for the simulation because the spatial correlation and location-specific features can be studied separately from different data sources and integrated in the simulation in a frequency domain. The proposed approach is illustrated by an example from the Rainbow petroleum play in the West Canadian Sedimentary Basin. In the application example, a pre-1994 exploration history data set was used as input, and the predictions are then checked against the locations of post-1993 exploratory drilling results. The comparison of the predictions from the proposed approach and the traditional conditional simulation shows that the model-based approach captures the essentials of geological controls on the spatial distribution of petroleum accumulation, thus improving the projections of undiscovered petroleum accumulations.     

广义S变换时频分析在广域电磁法数据处理中的应用

广义S变换是一种可逆的局部时频分析方法,是短时傅里叶变换和小波变换的延伸和推广。利用广义S变换对广域电磁时间序列数据进行分析,除了可以划定噪声区域便于实施滤波去噪之外,还可以提高数据频谱分析精度。针对广域电磁数据特点,基于广义S变换提出了适合广域电磁法的广义S变换处理流程,进行了信号的仿真实验。同时,也对广域电磁法实测数据进行了处理,高信噪比数据的处理结果证实了广义S变换的有效性和准确性。强干扰数据处理的结果表明,获取的信号质量比基于傅里叶变换的传统方法有明显提高。处理效果表明,提出的方法适用于广域电磁法数据处理。