伪普法零偏记录.docx

  1. 1、本文档共8页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:19108035856(电话支持时间:9:00-19:00)。
2011.11.1---2011.11.8 本周通过阅读文献及相关书籍主要学习了离散傅里叶变换的定义和二维离散傅里叶变换的基本原理,同时也学习了波动方程的有限差分方法,为之后学习伪普法合成零偏移距地震记录打下理论基础。 伪谱法也称傅立叶变换法,它的两个核心部分是:(1)借助快速傅立叶变换(FFT)计算空间偏导数;(2)采用二阶中心差分递归计算时间偏导数。主要优点是我们不仅可以得到记录的时间相位信息,而且还可以得到精确的地震波动力学信息一振幅,在计算波场空间精确导数时比其它方法需要更少的点数1.离散傅里叶变换(DFT) 傅里叶变换(FT)在地震数据处理中有着重要的作用。通过傅里叶变换把时空域的信号转换到频率或波数域,可以简化计算,减少计算工作量,由于计算离散傅里叶变换的有效快速算法的存在,离散傅里叶变换在地震数字信号处理的算法中起着核心作用。下面简要介绍DFT。 周期序列实际上只有有限个序列值有意义,因而它的离散傅里叶级数表示式也适用于有限长序列,这就得到有限长序列的离散傅里叶变换。实际上,可以把长度为N的有限长序列x(n)看成周期为N的周期序列的一个周期,这样利用离散傅里叶级数计算周期序列的一个周期,也就是计算了有限长序列。设x(n)为有限长序列,点数为N,即x(n)只在n=0~N-1有值,其他n时,x(n)=0。我们把它看成周期为N的周期序列x(n)的一个周期,而把x(n)看成x(n)的以N为周期的周期延拓,即表示成 Tf/Z)=£睥+洒)(1-2)同理,对频率域的周期序列X(k)也可看成是对有限长序列x(k)的周期延拓,而有限长序列x(k)看成周期序列X(k)的主值序列。因而我们可以定义有限长序列的离散傅里叶变换为: 正变换X( 正变换 X(k)=DFT[X(n)]=^TXn?W nkN 反变换X(n)=IDFT[X(k)]」云X(k)?W-nk0<n<N-1k=0 其中WN=ej为傅里叶变换权函数。通过上面的算式可以看出,获得每一个k的计算量都大得惊人,其时间复杂性为0(N2)因而在实时性要求较强的系统中,几乎没有实用价值。考察权函数W,其具有 周期性:WX(k+N)=WXk 对称性:W(N-X)k=(WXk)*其中(WXk)*为WXk的共轭复数这样,我们就能以权函数W这两个固有性质为基础,对X(n)序列作一些适当的变换,将一个较大的N分解成若干个较小的N,使整个傅里叶变换的计算过程变成一系列的迭代过程,就能大大减少运算次数,提高运算效率。 另外,值得注意的是x(n)和X(k)是一个有限长序列的离散傅里叶变换对。已知其中的一个序列,就能唯一地确定另一个序列。 二维离散傅里叶变换 u=0,1…MT,v=0,1…N-1x=0,1…M-1,y=0,1…N-1二维DFT定义为,设f u=0,1…MT,v=0,1…N-1 x=0,1…M-1,y=0,1…N-1 F(u,v)=劈2N1f(x,y)?Wux?Wvyx=0y=0 反变换 f(x,y)=勿无1F(u,v)?W-ux?W-vy MNu=0v=0 其中M,N为输入阵列点数,W=e-j3,W=e-j?分别为为x,y方向变换的权函数。二维DFT的直接运算比一维DFT要复杂得多,然而我们可以依据二维FT的可分离性,将其转化成一维FT,再运用FFT算法对其求解。 考察F(u,v)=勿2'f(x,y)?Wux?Wvy—^2VWux^2f(X,y)^WNvyx=0y=0x=0y=0 令gv)=21f(X,y)?",F(X,v)为f(X,y)在y方向N点的DFT,则y=0 F(",v)=2F",v)."为F3们在X方向M5时 于是,一个MXN阵列的二维DFT可转化成一次y方向N点的一维DFT和一次x方向M点的一维DFT,而这两次一维FT均可使用FFT算法。二维离散傅里叶反变换(IDFT)存在相同的可分离性,即可将MXN点频域阵列的IDFT转化成一次v方向N点的一维IDFT和一次u方向M点的一维IDFT,且先变换v方向或先变换u方向不影响结果。 弹性波动方程有限差分方法 有限差分方法基本思想是把连续的定解区域用有限个离散点构成的网格来代替,这些离散点称作网格的节点;把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似;把原方程和定解条件中的微商用差商来近似,积分用积分和来近似,于是原微分方程和定解条件就近似地代之以代数方程组,即有限差分方程组,解此方程组就可以得到原问题在离散点上的近似解。 有限差分法的主要内容包括:如何根据问题的特点将定解区域作网格剖分;如何把原微分方程离散化为差分方程组以及如何解此代数方程组。此外为了保证计算过程的可行和计算结果的正确,还需从理论上分析差分方程组的性态,包括解的唯一性、存在性和差分格式的相容性、收敛性和稳定性。 3.1声波方程有限差分格式的

您可能关注的文档

文档评论(0)

zl55662210 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

认证主体周**

相关文档

相关课程推荐