c语言雷克子波函数表达式,雷克子波频率研究

引用本文

张莹莹, 王云专, 石颖, 等. 2017. 雷克子波频率研究. 地球物理学进展, 2017, 32(5): 2162-2167, doi: 10.6038/pg20170542

ZHANG Ying-ying, WANG Yun-zhuan, SHI Ying, et al. 2017. Frequencies of the Ricker wavelet. Progress in Geophysics (in Chinese), 32(5): 2162-2167, doi: 10.6038/pg20170542

雷克子波频率研究

张莹莹1, 王云专1, 石颖1,2, 柯璇1

1. 东北石油大学地球科学学院, 大庆 163318

2. 黑龙江省普通高校科技创新团队"断层变形、封闭性及与流体运移", 大庆 163318

收稿日期：2017-02-10; 修回日期：2017-08-21.

基金项目：国家自然科学基金项目(41574117，41474118)、黑龙江省杰出青年科学基金项目(JC2016006)、大连理工大学海岸和近海工程国家重点实验室开放基金项目(LP1509)、黑龙江省自然科学基金项目(D2015011)和东北石油大学研究生创新科研项目(YJSCX2016-002NEPU)联合资助.

作者简介：张莹莹, 女, 1989年生, 在读硕士研究生, 主要从事地震资料处理方面的工作.(E-mail:569429532@qq.com)

通讯作者：王云专, 女, 1965年生, 硕士, 教授, 主要从事地震资料处理方面的工作.(E-mail:wyzmg@163.com)

摘要：雷克子波理论上是斯托克斯微分方程的解，它考虑到了牛顿黏度的影响，并适用于地震波通过黏弹性均匀介质的情况.本文定义了雷克子波的时域宽度及其频带宽度，并且分析了Lambert W函数的功能，采用功率谱估算平均频率，进而确定中心频率.文中证明了平均频率的标准差不是雷克子波的半频带宽度.此外，建立了计算理论频率(中心频率和半频带宽度)的数学方法，获得了理论频率与通过数值测量所得的频率(平均频率和标准差)之间的关系，并推导出了理论频率与测量频率关于雷克子波峰值频率的解析表达式.

关键词：雷克子波    Lambert W函数    功率谱    中心频率    标准差

Frequencies of the Ricker wavelet

ZHANG Ying-ying1, WANG Yun-zhuan1, SHI Ying1,2, KE Xuan1

1. Earth Science College of Northeast Petroleum University, Daqing 163318, China

2. Science and Technology Innovation Team on Fault Deformation, Sealing and Fluid Migration, Daqing 163318, China

Abstract:

The Ricker wavelet is theoretically a solution of the Stokes differential equation, which takes into account the effect of Newtonian viscosity, and is applicable to seismic waves propagated through viscoelastic homogeneous media. In this paper, we defined the time-domain breadth and the frequency-domain bandwidth of the Ricker wavelet and developed quantities analytically in terms of the Lambert W function. We determined that the central frequency, the geometric center of the frequency band, is close to the mean frequency statistically evaluated using the power spectrum, rather than the amplitude spectrum used in some of the published literature. We also proved that the standard deviation from the mean frequency is not, as suggested by the literature, the half-bandwidth of the frequency spectrum of the Ricker wavelet. Moreover, we established mathematically the relationships between the theoretical frequencies (the central frequency and the half-bandwidth) and the numerical measurements (the mean frequency and its standard deviation) and produced each of these frequency quantities analytically in terms of the peak frequency of the Ricker wavelet.

Key words:

Ricker wavelet

Lambert W function

power spectrum

frequency center

standard deviation

0 引言

雷克子波的波形理论上是斯托克斯微分方程的解，在数学上，雷克子波是一个高斯函数二阶导数(

频带的几何中心被称为中心频率.如果推导出雷克子波中心频率的解析表达式，就能够找到中心频率和平均频率之间的关系.后者是一个统计量，可以从地震数据的离散傅里叶谱得到.在已发表的文献中，通常使用功率谱(n≈2.15874时，n次谱估算平均频率可近似为中心频率.有的文献(

在本文中，我们分析了雷克子波的特点，描述了频率谱的关键参数.首先定义雷克子波在时间域的子波宽度和在频率域的频带宽度.这两个量分别是根据时间域波峰和频率域波峰的一半所定义的.根据以上二者的定