n—D Radon变换与波动方程偏移

本文把Ｒａｄｏｎ变换公式推广到任意ｎ维的情况。同时结合ｎ维Ｒａｄｏｎ变换和摄动理论提出了一种既能用于地面资料又能用于ＶＳＰ资料的偏移方法。     

一类广义Radon变换及反投影

本文提出了一类广义Ｒａｄｏｎ变换及反投影公式。这种变换形式比常规Ｒａｄｏｎ变换形式具有更广泛的意义及应用。数值计算结果表明，反投影恢复性较好。     

广义Radon变换的正交展开和反演

本文具体给出了广义Ｒａｄｏｎ变换在函数空间Ｌ＾２（Ｒ＾ｎ,Ｗｎ）上的正交函数展开形式,由此得到了广义Ｒａｄｏｎ变换的反演公式,从而推广了Ｌｏｕｉｓ和Ｄａｖｉｓｏｎ的工作。     

τ—P域中地震模型解释的若干问题

本文讨论了τ－Ｐ变换及其在地震方法中的应用，此变换已得到进一步发展，并用于研究压制多次波理论。τ－Ｐ变换时，由Ｐｅｔｒａｓｈｅｎ射线理论产生理论地震图，并通过Ｒａｄｏｎ变换实现τ－Ｐ变换，τ－Ｐ变换与决定检波器排列长度的近地表层厚度、Ｐ波速度以及Ｖ，／Ｖ，密切相关。由于海水层的均匀性，故水中声波速度最低，并且无需作静校正。τ－Ｐ变换适用于海洋地震探测。这种方法不仅可以用于深地震方法，而且也可用于浅     

数据严重缺少的Radon变换的反演

在许多实际问题中，由于客观条件的限制经常会遇到数据不全的Ｒａｄｏｎ变换的反演问题．在地球物理学的一些领域中更会遇到数据严重缺少的Ｒａｄｏｎ变换的反演问题──有限个点源及有限个方向投影问题．本文讨论此类反演问题，并给出具体的数值计算方法．     

分区连续函数的Radon变换的高精度反演法

Ｒａｄｏｎ变换及其反演是ＣＴ技术的数学基础，常用的反演算法虽然在函数比较光滑的情况下都能获得较满意的结果，但当函数有较强的奇性时精度就比较差．对函数在每一子区域中是光滑，但在不同于区域交界处有间断的情况，本文给出了有高精度的反演算法．     

沿圆曲线的Radon变换数值解

研究具有紧支集且在支集内连续的二元函数沿上半圆曲线的Radon变换反演问题。基于对投影函数的Fourier变换，反演问题可以归结为具有弱奇性及震荡核的Abel积分方程的求解。我们证明了当圆曲线中心及半径在一定范围内变化时，在已知沿上半圆曲线的Radon变换情况下，这个积分方程的解具有唯一性，并给出了消除Abel积分方程弱奇性的数值方法。在考虑投影数据噪声的情况下，给出了多次加权改善系数矩阵条件数稳定的数值方法，并通过数值模拟验证所提出方法的有效性。      

台湾地区地震活动的时间多分形特征研究

利用板块构造、活断层分布、地壳运动这三个方面的最新资料将台湾地区划分为东、西两个地震区和六个地震带。对各地震带11a来的地震目录计算了广义维数Dq和多标度分形谱f(α)，然后绘制成Dq—q曲线和多标度分形谱f(α)—α以及分维时程曲线。综合分析后发现：(1)地震活动的时间结构是多分形而非单一分形结构；(2)台湾地区在大震前后的1～2a内，大震所在的地震带内地震时间结构的Dq—q曲线和多标度分形谱f(α)—α大都出现异常形态；(3)功时程曲线分析可以为强震的时间预报提供有效的客观依据。     

陇东黄土塬地区中生界石油砂岩储层预测

在鄂尔多斯盆地陇东黄土源地区地震解释中,采用广义S变换地震层序处理解释技术、三叠系顶部古地貌解释技术,以及各种储层横向预测技术,对比分析了应用效果。研究表明,用广义S变换地层层序处理解释技术可识别出约6m的层序界面,进而确认了陇东地区长8段砂体的展布方向,扩大了陇东地区西峰油田的勘探范围。     

大华北10次中强震前后广义时空分维的变化及其在中短期预报中的应用

本文应用高安秀树提出的广义分维方法分析计算了１９７０年以来大华北地区１０次Ｍｓ≥５．８级地震前后震中周围小震活动广义时（Ｄｔ）空（Ｄｓ）分维的变化，结果表明中强震前普遍检测到Ｄｔ，Ｄｓ的低值异常。具有实用意义的是至少有一条曲线可以统一地取Ｄｔ，Ｄｓ低于震前数年一倍均方差线，而另外一条曲线较同步地低于均值线达半年之久作为异常判别标准，中强震一般于异常结束后３－９个月发生，其中个别地震于异常之中发生。     

THE QUANTITATIVE INTERPRETATION OF DIPOLE SOUNDINGS BY MEANS OF THE RESISTIVITY TRANSFORM FUNCTION*

In this paper a theorem is demonstrated which allows—after the introduction of a suitable dipole kernel function or dipole resistivity transform function—to write the apparent resistivity function as an Hankel transformable integral expression. As a practical application of the theorem a procedure of quantitative interpretation of dipole soundings is suggested in which the dipole resistivity transform function obtained after inversion of the original dipole apparent resistivity data is used to control the goodness of the set of layering parameters which have been derived with our previous method of transformation of dipole sounding curves into equivalent Schlumberger diagrams.     

Intensity measure conversion of fragility curves

In seismic risk assessment of structures, fragility functions are the probabilistic characterization of vulnerability at the component and/or structural level, expressing the probability of failure as a function of a ground motion intensity measure (IM). Fragility curves, in general, are structure- and site-specific, thus a comparison of fragility curves, then of vulnerability, is not straightforward across multiple structures. Also, it could be the case that hazard at a site of interest is not available for the IM originally considered in the fragility assessment. These situations require to convert fragility curves from an original IM to a target one. The present study addresses a hazard-consistent probabilistic framework for converting spectral acceleration-based IMs from an original IM to a target IM at a given site. In particular, three conversion cases, under different assumptions on the explanatory power of the involved IMs with respect to structural failure, are discussed: (a) a vector-valued IM consisting of the original and target IMs, magnitude, and source-to-site distance; (b) a vector-valued IM consisting of the original and target IMs; and (c) the original (scalar) IM only, assuming that structural response, given the IM, is statistically independent of the other ground motion variables. In this framework, the original fragility functions are characterized using the state-of-the-art methods in performance-based earthquake engineering, then the fragility curves as a function of the target IM are evaluated through applications of the probability calculus rules, ensuring consistency with the seismic hazard at the site of interest. The conversion strategy is illustrated through the applications to three-, six-, and nine-story Italian code-conforming reinforced concrete buildings designed for a high-hazard site in Italy. The study shows that, in most of the cases, the converted fragility curves have agreement with the reference curves directly developed in terms of the target IM. Cases in which least agreement was found are likely due to the models used to obtain the terms required by the conversion equations.     

DEFINITION AND APPLICATION OF THE FREQUENCY -EFFECT TRANSFORM FUNCTION TO THE INTERPRETATION OF IP SOUNDINGS*

This paper deals with a new method of quantitative interpretation of induced polarization soundings in the frequency-domain. From the general expression of the apparent frequency-effect for soundings carried out on a multi-layered earth the application of Hankel's inversion theorem allows to introduce a new function, called here the “frequency-effect transform”. The new interpretation method consists of two steps: 1) the inversion of field data to obtain the frequency-effect transform graph and 2) the analysis of this graph to derive the layering parameters. The first step is performed by means of a slightly revised version of a simple numerical procedure, previously suggested by the author for the inversion of d.c. resistivity sounding data. The second step is carried out by a complete curve-matching procedure, applied directly on the transform graph. This implies suitable master curves, whose preparation doesn't meet all the mathematical difficulties which are present when preparing master curves of the apparent frequency-effect function.     

基于拉普拉斯变换构建校正窄频速度记录的传递函数

受窄频地震仪平坦响应范围影响,窄频速度记录存在低频成分失真问题,导致地震记录可用范围受限。针对此问题,本文推导基于拉普拉斯变换和双线性变换的传递函数,实现由窄频地震记录向宽频地震记录的校正,并以日本Hi-net速度记录为例进行验证,将校正后的速度记录与同台KiK-net加速度积分所得的速度记录予以对比。结果显示,原始速度记录在低频处存在失真,而校正后的波形与KiK-net加速度积分速度记录波形一致,这表明改进的传递函数能有效地解决原速度记录中的低频成分失真问题,有效地拓宽了低频可使用范围,而且相较于Nakata校正速度记录方法,以本文给出的传递函数校正的速度记录在振幅和波形等方面的精度更高。     

EMD Method Applied to Identification of Logging Sequence Strata

In this work, we compare Fourier transform, wavelet transform, and empirical mode decomposition (EMD), and point out that EMD method decomposes complex signal into a series of component functions through curves of local mean value. Each of Intrinsic Mode Functions (IMFs — component functions) contains all the information on the original signal. Therefore, it is more suitable for the interface identification of logging sequence strata.Well logging data reflect rich geological information and belong to non-linear and non-stationary signals and EMD method can deal with non-stationary and non-linear signals very well. By selecting sensitive parameters combination that reflects the regional geological structure and lithology, the combined parameter can be decomposed through EMD method to study the correlation and the physical meaning of each intrinsic mode function. Meanwhile, it identifies the stratigraphy and cycle sequence perfectly and provides an effective signal treatment method for sequence interface.     

THE DETERMINATION OF FILTER COEFFICIENTS FOR THE COMPUTATION OF STANDARD CURVES FOR DIPOLE RESISTIVITY SOUNDING OVER LAYERED EARTH BY LINEAR DIGITAL FILTERING *

The technique of linear digital filtering developed for the computation of standard curves for conventional resistivity and electromagnetic depth soundings is applied to the determination of filter coefficients for the computation of dipole curves from the resistivity transform function by convolution. In designing the filter function from which the coefficients are derived, a sampling interval shorter than the one used in the earlier work on resistivity sounding is found to be necessary. The performance of the filter sets is tested and found to be highly accurate. The method is also simple and very fast in application.     

Rayleigh波频散曲线“交叉”及多模式耦合作用研究

Rayleigh波可以用来反演近地表结构,在工程物探、石油物探、地球内部结构探测中均有重要意义.数值计算得到的含低速层的层状介质对应的Rayleigh波频散曲线会出现看似“交叉”的现象,但是对于这种现象目前还没有进行系统的研究.事实上可以验证,有些看似交叉的频散曲线实际上不相交.改变低速层的厚度和横波速度发现低速层越明显（即低速层速度越低或层厚越厚）频散曲线越不容易相交.凡友华等在2007年提出频散曲线对应着四种基本模式,在频散曲线发生“交叉”现象的区域实际上存在两个以上模式的频散曲线.本文主要研究了存在R模和S₂模的区域内频散曲线的“交叉”现象.首先利用竖直本征振动曲线研究R模和S2模Rayleigh 波的振动特点,发现R模对应的本征振动主要集中在地表,随着深度变化能量快速衰减,S₂模对应的本征振动主要集中在第2层.研究“交叉点”附近频散点对应的本征振动曲线发现这一区域有些Rayleigh波同时具有R模和S₂模的振动特点,对应着一种耦合模式.通过对实例的研究发现,在“交叉点”附近,若两条频散曲线不发生交叉,则每条曲线对应的模式会发生R模和S₂模之间经由耦合模式的转变,本文称这种现象为两种模式发生耦合;若两条频散曲线相交,则同一条频散曲线上的Rayleigh波模式几乎相同,只是在离交点很近的区域会存在一些耦合模式,本文称此时两种模式不发生耦合.本文研究结果主要供Rayleigh波对低速层结构的反演研究参考.     

Contourlet based seismic reflection data non-local noise suppression

In this paper, we propose a non-local, transform domain noise suppression framework to improve the quality of seismic reflection data. The original non-local means (NLM) algorithm measures similarities in the data domain and we generalize it in the nonsubsampled contourlet transform (NSCT) domain. NSCT gives a multiscale, multiresolution and anisotropy representation of the noisy input. The redundancy information in NSCT subbands can be utilized to enhance the structures in the original seismic data. Like the wavelet transform, NSCT coefficients in each subband follow the generalized Gaussian distribution and the parameters can be estimated using appropriate techniques. These parameters are used to construct our proposed NSCT domain filtering algorithm. Applications for synthetic and real seismic data of the proposed algorithm demonstrate its effectiveness on seismic data random noise suppression.     

曲面上的位场转换（二维和三维情况）

在文献[1]中我们利用复变函数论的方法只解决了二维的曲面上的位场转换问题.在本文中,我们利用单层位势理论将它化为一个奇异积分方程求解的问题,并且提出了用奇点分离法来解此奇异积分方程,这对二维和三维情况都适用.通过理论模型的计算,结果的精度是令人满意的.     

小波变换在裂缝识别和裂缝密度评价中的研究(英文)

通过常规测井曲线小波变换来评价裂缝参数是一个较新的研究课题。本文针对火山岩裂缝储层,对常规测井曲线做coif5、bior4.4和db5基小波的小波变换,通过变换后的分解信号与成像测井裂缝密度的对比研究,寻找出与裂缝密度相匹配的基小波,再将常规测井曲线在这个基小波下的分解信号与曲线变化率法相结合建立裂缝指示曲线来识别裂缝发育段,最后通过裂缝指示曲线与成像测井的裂缝密度之间的量化关系来比较准确的评价裂缝密度。将裂缝指示曲线法在松辽盆地南部进行了应用,通过裂缝指示曲线求取的裂缝密度与成像测井的裂缝密度之间的相对误差为0.53。