首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Spectral analysis of the deformation of a homogenous electric field caused by a long, thin inclined bed, which is of considerable importance in the exploration of ground water and minerals, is presented using the Hartley transform. The Hartley transform is an alternative and real replacement for the well‐known complex Fourier transform in the field of spectral analysis. The thickness of the bed and the inclination are given as functions of frequency by simple expressions. A theoretical example illustrates the method while the applicability is demonstrated by the field examples from the fractured crystalline basement complex in Burkina Faso, Africa and the Precambrian limestones of the Cuddapah basin, Andhra Pradesh, India. The results obtained by this method agree well with those of the drilling.  相似文献   

2.
提出直接在序率域内用Walsh变换实现引力场球谐综合的问题。给出球谐函数展开式的Walsh变换及快速算法,讨论了Walsh变换和Walsh-Fourier变换、Fourier变换之间的差异,分析了用地球重力场模型OSU81的位系数作出的Walsh变换和Fourier变换的结果。研究表明:Walsh变换与Walsh-Fourier变换、Fourier变换对应向量在数量方面的差值通常都小于士10~(-5);对于给定的阶数和飞行高度,3种方法求得的球谐综合值总是完全一致的;1°×1°等网格数据和Walsh函数形状相近。在重力场研究中Walsh级数会比Fourier级数收敛得更快;Walsh变换在计算速度、计算准确度、数据储存、收敛速度和方法简单方面都好于Fourier变换。  相似文献   

3.
Self-potential anomalies due to two-dimensional inclined sheet of finite depth extent are analysed from their horizontal and vertical derivatives via the Hilbert transform. The depths to the top and bottom of the sheet and the inclination are obtained by means of simple mathematical expressions. The method is applied to a theoretical example and to a field S.P. anomaly over a sulphide deposit in the Kalava fault zone of Cuddapah basin, India. The presence of random noise is analysed for various levels and its effect is discussed. This interpretation process can be automated for all practical purposes by simple programming.  相似文献   

4.
提出了一种新的发震构造型式-转换破型型。其特征是地震不发生在主干断裂和共轭断裂面上,而是在其两侧一定距离的岩块中;发震断裂不与其平行而是成一定的交角;震源机制解的三个力轴与造力应场相对应的三个力轴不一致,而是有规律地旋转了一定有角度。  相似文献   

5.
金星  陈超  张明宇 《地震学报》2000,22(5):509-516
将震源简化为有一定破裂长度的线状断层,定量研究了单侧破裂震源对地震动空间相关性的影响.依据波谱随空间坐标的展开式,得到了描述地震动场相关性的表达式.在远场假定下,给出了理论解析公式,并依据地震动Fourier谱随空间坐标的展开条件,得到了计算地震动场孔径的理论公式.最后,以一次地震为例,给出了孔径随震中距、方位角、频率变化的三维图象.   相似文献   

6.
采用稳定型双共轭梯度快速Fourier变换(BCGS-FFT)算法精确计算二维均匀介质中的积分方程.采用一种新的插值函数作为基函数和试探函数对积分方程进行弱化离散,离散后的积分方程采用稳定型双共轭梯度迭代方法进行求解,从而得到异常体内电场的分布.计算时采用快速Fourier变换技术将积分方程内Green函数与电场的乘积表示成褶积形式以加快计算速度.数值计算举例说明了算法的精确性和有效性.  相似文献   

7.
Hartley变换化极   总被引:1,自引:1,他引:0       下载免费PDF全文
骆遥 《地球物理学报》2013,56(9):3163-3172
本文采用一种新颖的积分变换Hartley变换,系统推导了其频率域化极方法,提出在实数空间利用Hartley变换进行磁异常化极.结合现有的低纬度化极方法,讨论了低纬度Hartley变换化极特性,实现低纬度化极.理论模型计算表明建立的Hartley变换化极方法准确、可靠,实际资料处理表明该化极方法具有实用性.  相似文献   

8.
Potential, potential field and potential‐field gradient data are supplemental to each other for resolving sources of interest in both exploration and solid Earth studies. We propose flexible high‐accuracy practical techniques to perform 3D and 2D integral transformations from potential field components to potential and from potential‐field gradient components to potential field components in the space domain using cubic B‐splines. The spline techniques are applicable to either uniform or non‐uniform rectangular grids for the 3D case, and applicable to either regular or irregular grids for the 2D case. The spline‐based indefinite integrations can be computed at any point in the computational domain. In our synthetic 3D gravity and magnetic transformation examples, we show that the spline techniques are substantially more accurate than the Fourier transform techniques, and demonstrate that harmonicity is confirmed substantially better for the spline method than the Fourier transform method and that spline‐based integration and differentiation are invertible. The cost of the increase in accuracy is an increase in computing time. Our real data examples of 3D transformations show that the spline‐based results agree substantially better or better with the observed data than do the Fourier‐based results. The spline techniques would therefore be very useful for data quality control through comparisons of the computed and observed components. If certain desired components of the potential field or gradient data are not measured, they can be obtained using the spline‐based transformations as alternatives to the Fourier transform techniques.  相似文献   

9.
Acoustic plane wave scattering at a vertical fault structure represents the simplest two-dimensional model of geophysical exploration that can be investigated by analytical techniques. The exact and complete solution, in the time domain, for the scattering of the pressure field of an acoustic plane wave normally incident on a vertical fault structure is determined adapting previous results given for the frequency domain. The wave form of the pressure field of the incident plane wave is expressed by a causal time function that decays exponentially with time at every point above the fault (z<0). The zero-order term of the scattered pressure field has been computed above the fault. This zero-order term consists of an inverse Fourier transform which reduces to a closed expression forx=0, and contains an integral of a Hankel function forx#0. The high frequency part of the inverse Fourier transform forx#0 is computed employing asymptotic expressions for the Hankel function. The integral of the asymptotic expression of the Hankel function reduces to: (i) a Fresnel integral which contains a plane wave term for |x||z|; and (ii) a stationary point plane wave term plus an upper limit term for |x|=O(|z|). For the latter case the plane wave term cancels, leaving a cylindrical wave emanated from the edge of the fault. The wave front is well defined in shape, in phase and in amplitude. The amplitude of the scattered field is discontinuous atx=0, presents a jump and is well defined for |x| small and is rather smooth for |x| large.  相似文献   

10.
11.
The gravity effect of an infinite horizontal trapezoidal prism is derived and its Fourier spectrum is analyzed so as to yield information about four parameters of the causative structure, namely the depths to the upper and lower surfaces, width of the upper surface, and the inclination of the sides. In order to test the applicability of the method, synthetic data are constructed by digitizing the theoretical gravity effect. Subsequently, the corresponding Discrete Fourier Transform (DFT) is obtained. The parameters evaluated from the DFT are observed to be sufficiently close to the chosen values.  相似文献   

12.
在2008年5月12日汶川MS8.0地震中,四川数字强震台网共获取了133组三分向加速度记录. 本文选取了一些不同断层距的台站所获取的强震动记录进行了处理和分析.在数据处理中,采用基于聚类经验模态分解(EEMD)提取信号时频特性的方法,有效获得了信号能量的时频分布,提取了中心频率、 Hilbert能量、最大振幅对应的时频等特性,并与傅里叶变换、小波变换进行了对比研究.研究结果表明, 对非线性的强震记录采用聚类经验模态分解(EEMD)能抑制经验模态分解(EMD)中存在的模态混叠问题;与傅里叶变换和小波变换相比发现, HHT边际谱在低频处幅值高于傅里叶谱;与小波变换受到所选取的母波强烈影响不同, HHT直接从强震记录中分离出固有模态函数(IMF),更能反映出原始数据的固有特性, Hilbert谱反映出大部分能量都集中在一定的时间和频率范围内,而小波谱的能量却在频率范围内分布较为广泛.因此,基于EEMD的HHT在客观性和分辨率方面都具有明显的优越性,能提取到更多强震加速度记录的时频特性.   相似文献   

13.
New formulae are presented for the calculation of the horizontal magnetic field due to a point electrode situated at the vertical boundary between two quarter-spaces of different electrical conductivities. Previously the solution was obtained using a double Fourier transform of the expression for the vertical component of the magnetic field. The new formulae are given in the form of single integrals.  相似文献   

14.
An efficient method of computing spectrum and cross-spectrum of large scale aero-magnetic field (or of any other two-dimensional field) has been developed and programmed for a digital computer. The method uses fast Fourier transform techniques. Briefly, the method is as follows: a digitized aeromagnetic map is divided into a number of rectangular blocks. Fourier transforms of these blocks are computed using a two-dimensional fast Fourier transform method. Finally, the amplitude of the Fourier transforms is averaged to give the desired spectrum. Computation of cross-spectrum follows the same lines. In fact, the same programme may be used to a compute the spectrum as well as cross-spectrum. The method has a number of computational advantages, in particular it reduces greatly computational time and storage requirements. The programme has been tested on synthetic data as well as on real aeromagnetic data. It took less than 30 seconds on an IBM 360/50 computer to compute the spectrum of an aeromagnetic map covering an area of approximately 4500 square miles.  相似文献   

15.
The formal equivalency between a Hankel transform and an Abel transform followed by a Fourier transform is used as a basis for the computation of synthetic seismograms and for performing plane-wave decomposition of both synthetic and field seismograms. The Abel-Fourier method performs better for near-offset or small ray-parameter, although at greater computational effort than calculations based upon the asymptotic Hankel transform.  相似文献   

16.
利用Hilbert变换计算重力归一化总梯度   总被引:12,自引:3,他引:9       下载免费PDF全文
针对提高重力勘探正反演解释的分辨率问题,提出利用Hilbert变换计算和研究重力归一化总梯度.文中从理论上详细地证明了方法的可行性,给出了适合计算机实现的计算方法.为探讨该方法对油气藏的分辨能力,在模型计算中,分别利用Fourier级数法、Fourier变换法和本文提出的Hilbert变换法计算模型的GH场值,发现当三度体储油球冠模型(模拟似三度贮油气藏背斜模型)油气藏部分(低密度体)厚度减小到低于球冠厚度的十分之一时,Fourier级数及变换法不能分辨出低密度体所产生的异常,而Hilbert变换法仍然可以清晰地识别,这说明用Hilbert变换法计算的GH场对异常识别的分辨率优于其他两种方法.  相似文献   

17.
The spectrum of a magnetic or a gravity anomaly due to a body of a given shape with either homogeneous magnetization or uniform density distribution can be expressed as a product of the Fourier transforms of the source geometry and the Green's function. The transform of the source geometry for any irregularly-shaped body can be accurately determined by representing the body as closely as possible by a number of prismatic bodies. The Green's function is not dependent upon the source geometry. So the analytical expression for its transform remains the same for all causative bodies. It is, therefore, not difficult to obtain the spectrum of an anomaly by multiplying the transform of the source geometry by that of the Green's function. Then the inverse of this spectrum, which yields the anomaly in the space domain, is calculated by using the Fast Fourier Transform algorithm. Many examples show the reliability and accuracy of the method for calculating potential field anomalies.  相似文献   

18.
北京八宝山断裂带主要由八宝山内支、外支断裂和黄庄-高丽营断裂组成。根据野外地质资料和室内研究,认为八宝山内支断裂是燕山期的逆掩断层,已不具活动性,而外支断裂是一条规模和断距都不大的正断层,只有黄庄-高丽营断裂才是活动断层,它错断了中更新统黄土。 为了获得更多的有关八宝山断裂带活动性的资料,我们采用Kanaori提出的一种判定断层活动性的新方法,通过扫描电镜观察了八宝山断裂带的6个断层泥样品和4个有关地层中的石英颗粒表面微结构,所得结果与野外宏观地质资料相吻合  相似文献   

19.
Various authors have recently studied the Walsh and Hadamard transforms for functions with one or two variables. Like the Fourier transform they may be obtained by using the Cooley-Tukey algorithm. They can be calculated all the faster because Walsh functions have only values equal to + 1 or ? 1. The complex number multiplications in the Fourier transform are replaced by additions and subtractions of real numbers. These transforms are more general than that of Fourier because the concept of sequency, i.e. the mean number of zeros per time unit, is substituted for that of frequency. These transforms already have numerous applications to signal transmission. We describe some in the field of seismic prospecting, in particular the compression of the information contained in seismic traces.  相似文献   

20.
A fast algorithm is presented for numerical evaluation of forward and inverse Radon transforms. The algorithm does not perform exact one-to-one mapping as the discrete Fourier transform but, due to the use of band-limited basis functions, it is robust and sufficiently accurate for seismic applications. By rewriting the transform as a convolution, a computational speed is obtained similar to the speed of the 2D fast Fourier transform.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号