Hartley变换在位场（重、磁）异常导数计算中的应用

导数计算是位场（重磁）数据处理中必不可少的技术手段,现今大多采用Fourier变换来进行。Hartley变换是在Fourier变换基础上定义的一种实数域运算,比Fourier变换更加对称,所需的运算量更少。笔者推导出基于Hartley变换的重磁异常导数计算公式,计算结果与理论值之间误差小于5%,通过理论模型证明Hartley变换可替代Fourier变换进行位场异常的导数计算,且受噪音干扰较小。将Hartley变换用于位场边界识别滤波器的计算,获得了清晰的断裂分布。     

运用Hartley变换模拟地震波在正交各向异性介质中的传播

在地震传播理论中,地震波场的正演数值模拟一直是研究的热点。而在正演数值模拟方法的研究中,计算精度和计算效率是评价此方法的有效性及优越性的二个关键问题。伪谱法在计算精度与计算效率方面优越性十分明显,常用的有Fourier变换法和Hartley变换法;虽然Hartley变换法在求取导数时较Fourier变换法复杂,但由于Fourier变换法的计算同时涉及复数的实部与虚部,在计算速度和占用内存方面不如Hartley变换法。这里详细地阐述了利用Hartley变换求解正交各向异性介质波动方程的数值方法,并模拟了多种正交各向异性介质模型,对比分析了地震波在各向同性介质与正交各向异性介质中的传播差异。结论表明,该方法能正确、高效、直观地反应地震波在正交各向异性介质中的传播规律。     

地震信号广义时频分析及其数值实现

基于分数阶Fourier变换的信号时频分析称为广义时频分析,将其引入地震信号处理中,是对现有基于传统变换时频分析方法的广义化拓展。这里着重数值实现,探索基于分数阶Fourier变换的时频分析,并应用了其在物探技术中新思路和比较的优势。在分数阶Fourier变换基本定义的基础上,数值实现离散分数阶Fourier变换;研究三类广义时频分析方法及其数值实现;针对地震信号,以伪Wigner分布为例,展示广义时频分析效果。经研究表明:实现的数值算法正确有效;地震信号新的广义时频分布能量更加集中,时频分析的聚焦度和信息特征更加清晰。由此可见,广义时频分析为地震信息特征分析、预测和提取提供新的视角和工具,值得进一步研究和推广。     

全场冰漂流测量技术研究

结合国内先进的航海雷达测冰技术,介绍了全场冰漂流测量技术的互相关分析原理及其Fourier实现方法, 以及冰漂流场中部分错误矢量的修正方法.利用全场冰漂流测量技术, 通过对一些实测的雷达海冰图像进行互相关分析处理, 获得的二维海冰漂流场的矢量分布图, 充分验证了全场冰漂流测量技术的实用性和可靠性.     

基于二维小波变换的混沌图像加密算法

为了提高图像信息的安全性,利用混沌序列的特性,提出了一种新的基于混沌序列的图像加密技术。在对图像进行二维小波变换后,以混沌序列为基础,利用图像置乱变换技术,实现基于频域的数字图像加密算法。该算法具有较高的安全性,实验结果令人满意。     

基于小波变换的遥感影像压缩编码

基于小波变换的图像压缩方法能在高压缩比的前提下保持好的重建图像质量,它的多分辨率特性与人类视觉系统非常相近,小波变换的编码技术被广泛地应用于遥感影像的压缩中。介绍了小波理论的相关知识,着重对嵌入零树小波、等级树分割、小波数据的形态学表示算法等方法进行了详细的论述,并对图像编码技术的发展方向进行了展望。     

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

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

基于小波变换的岩石声发射信号互相关分析及时延估计

声发射全波形采集为基于波形分析的声源定位方法提供了可能。声发射源定位方法中最常用的是时差定位方法,因此,时差获取精确与否直接影响定位的精度。借助于小波变换和互相关技术研究声发射信号的时延估计,首先利用小波分析技术对信号进行分解,确定有意义信号的频带宽度,再提取分解后的相应频带的小波系数,利用互相关技术计算出该频带信号到达各传感器的时差,进而利用时差定位法反演声发射源,实现了声发射源更精确的定位。经声发射监测验证结果表明,该时延估计方法能有效地减小声发射源定位的误差,提高声发射源的定位精度,为准确反演声发射源提供了一种有效途径。     

基于全相位多次互相关的Geiger声发射源定位方法

直接互相关在岩石声发射源定位中,会出现多次峰值与波形密集的问题,在求解过程中会出现无解的情况且定位精度不足,为了解决这些问题,提出一种基于全相位快速傅里叶变换（FFT）分析及多次互相关的岩石声发射源定位的算法。通过对50 mm×100 mm×50 mm的花岗岩进行断铅试验,利用小波变换对互相关分析的频谱进行预处理,同时发挥全相位FFT的相位不变性的优点,找出波形峰值,得到时间差,最后利用最小二乘法得到定位初始点,再通过Geiger迭代,得到最终定位坐标。在选取8个测点的基础上,经声发射监测试验结果表明：该方法比单独直接互相关更容易找到波峰,且波形清晰,初始点的选择有效解决了Geiger迭代法对初始值要求严格的问题,通过全相位分析及多次互相关计算,然后Geiger程序迭代,得到的定位结果比直接互相关精度更高,平均绝对误差比美国生产PCI–2型声发射仪器Geiger定位的结果低 2.40 mm,比基于多次互相关和Geiger算法的声发射源定位结果低3.31 mm,比全相位相位差法声发射测试时间延迟定位结果低1.74 mm 。     

大地电磁测深数据的时频分析

大地电磁测深(MT)数据表现出非线性、非平稳、非最小相位特征,不符合以Fourier变化为基础的传统谱分析的基本要求.为寻求更好的谱分析方法,把最新发展的Hilbert-Huang变换(HHT)引入到大地电磁测深数据的时频分析当中来,针对湖北宜万线实测的MT信号,与短时Fourier变换(STFT)、Winger-Ville分布(WVD)和小波变换(WT)进行了时频分析比较.比较分析表明:HHT方法克服了其它一些方法的缺陷,完全取消了窗函数的作用,其结果不受核函数影响与时频测不准原理限制,具有完全的局部时频特性和更好的时频聚集性.分析认为HHT是一种全新而更优越的分析与处理大地电磁测深信号的时频方法,能更好地提取信号的时频本质特征,用于指导生产勘测实践.     

Hartley spectral analysis of self-potential anomalies caused by a 2-D horizontal circular cylinder

The real spectral analysis of SP anomalies due to 2-D horizontal circular cylinder is carried out using Hartley transform. The Hartley transform is an alternate means of realizing spectral analysis in real domain unlike the Fourier spectral analysis. The Hartley transform yields a straightforward interpretation of SP anomalies caused by horizontal circular cylinder wherein all the parameters are derived independently as a function of frequency. A theoretical example illustrates the procedure. The effect of random noise, up to 10% of white Gaussian noise, on the interpretation scheme was also studied and found to be of negligible. The field example of the ??Sulleymonkey?? anomaly in the Ergoni copper district, Turkey exemplifies the applicability of the proposed method.     

实信号快速傅立叶变换在地震信号处理中的应用

地震信号的数据量非常大,如果简单利用复序列快速傅立叶变换算法对地震信号进行处理,将会成倍增加数据量,增加计算量,因此针对地震信号这一实数序列,讨论实信号的快速傅立叶变换 。文章提出了3种实序列的快速傅立叶变换算法,对运算量的减少进行了定量的分析,并对实际地震数据进行了计算,从分析及计算结果可以看出,这3种方法可以明显减少运算量,同时节省了计算机的存储空间。     

关于瞬变电磁法2.5维正演中的几个问题

这里讨论了瞬变电磁法2.5维正演模拟中的两个积分变换:傅里叶逆变换和拉普拉斯逆变换。针对发收距为零的中心回线方式瞬变电磁法,提出了在傅氏域中考察傅氏变换函数随波数的变化规律,进而根据曲线首尾支渐近线来划定波数覆盖的范围,然后同解析解对比,确定出最少个数的傅氏域波数。另外,介绍了只需对较少的拉氏变换变量作纯实数运算的拉普拉斯数值反演计算方法。通过对均匀半空间表面上垂直磁偶极子源形成的瞬变电磁场进行正演模拟,结果表明,关于傅氏域波数的范围划定原则及个数选取方案是合理的、拉氏逆变换算法也是切实可行的。     

Spectral-type simulation of spatially correlated fracture set properties

Fracture set properties such as orientation, spacing, trace length, and waviness tend to be spatially correlated. These properties can be efficiently simulated by spectral analysis procedures that take advantage of the computational speed of the fast Fourier transform. The covariance function of each property to be simulated is obtained from the variogram function estimated from mapped fracture set data and is typically referenced to the mean vector of the set. Simulation procedures for normally and exponentially distributed data involve generating uncorrelated Fourier coefficients that are assigned proper variance according to the spectral density, which is the Fourier transform of the covariance function. These coefficients are then reverse Fourier transformed to produce simulated set properties that have the desired variance and variogram function.     

时频分析在石油地球物理勘探中的应用

时频分析方法经历了富氏变换、短时富利叶变换、小波变换三个阶段。随 改进,其在地学中的应用及研究范围也在不断拓宽,从对有关地质信号的能量补偿、滤波、压缩等处理,发展到通过信号的时频特征对地层的变化特征、沉积旋回性质等方面的研究。主要阐述了时频分布的概念,方法的改进及其在石油地球物理勘探中的应用和前景。     

地震波混合阶褶积算法模拟

基于正反傅立叶变换,提出了地震波模拟的混合阶褶积算法。该方法原理简单、易于实现,结合了有限差分与伪谱法的优点,具有较高的精度和计算效率,适用于地震波场的正演计算,同时还给出了二维地震波场的理论计算实例。计算结果表明,此算法模拟结果正确、精度高、速度快、能适应较为复杂地质模型,并且易于推广到各向异性介质中去。     

基于分数阶傅里叶变换的探地雷达子波提取算法

针对探地雷达子波提取问题,提出了一种基于分数阶傅里叶变换的子波提取算法。假设探地雷达地下介质模型为线性系统并是时变的,并且子波在传播过程中不变,把子波和地层信息描述为一个时间的积分形式。分数阶傅里叶变换是酉变换,具有保范性,把子波的提取问题转化为在分数阶傅里叶变换域上的最优滤波问题。对分数阶傅里叶变换的阶数p在(－2,2]范围以1/300为步长进行了搜索,在分数阶傅里叶变换域上找到最优解,最后经过分数阶傅里叶变换反变换得到子波估计的时域信号。仿真实验首先验证了分数阶傅里叶变换最优滤波算法的正确性,然后利用探地雷达实验数据提取了子波,验证了该算法在提取子波中的有效性。