张叔伦
- 作品数:25 被引量:148H指数:8
- 供职机构:大连理工大学数学科学学院更多>>
- 发文基金:国家自然科学基金地震科学联合基金中国石油天然气集团公司科技项目更多>>
- 相关领域:天文地球石油与天然气工程理学自动化与计算机技术更多>>
- 双平方根叠前深度偏移的广义高阶屏方法被引量:10
- 2001年
- 叠前深度偏移技术既可在炮点—接收点域实现,也可在共中心点—炮检距域实现。前人将裂步延拓算子 推广到共中心点—炮检距域相移法双平方根叠前偏移中。在波场向下延拓的每一步长内,仅通过一次时移量来校正常速相移延拓产生的误差,得出了简单、高效的裂步双平方根叠前深度偏移方法,但精度较低。本文基于波场延拓的非稳态相移公式,通过引入参考速度,并对双平方根项中的两个平方根项作泰勒级数展开,经过适当的数学推导,得出了共中心点—炮检距域波场延拓的双平方根非稳态相移新的高阶屏近似公式。该公式可直接 在共中心点—炮检距域高效地延拓叠前波场。通过增加屏的阶数,提高了剧烈横向变速条件 下叠前深度偏移的精度。它是裂步双平方根叠前深度偏移方法的推广。理论模型和实际资料 的试算结果表明,本文方法是有效而实用的。
- 孙沛勇张叔伦
- 关键词:地震勘探
- 按速度分区的裂步傅里叶叠前深度偏移方法被引量:1
- 2004年
- 本文提出了一种按速度构造对速度场进行分区的裂步傅里叶叠前深度偏移方法。裂步傅里叶法 (SSF)的理论基础是摄动项 ,在一个步长的波场延拓中 ,当速度反差较大时 ,摄动项的误差也随之加大。所以当横向速度剧烈变化时 ,要根据速度构造对速度场进行分区 ,在每一个速度区内 ,分别采用不同的参考速度 ,这样可以最大限度地减小参考速度和真实速度之间的反差 ,进而提高裂步法成像精度 ,同时可以最大限度地减少参考速度的个数 ,使计算量降到最小。实现步骤为 :1在每一个延拓步长内 ,首先求出某层的最大速度和最小速度 ,然后从两者之间确定一个合适的分界速度 ;2以分界速度为界 ,将该速度场分为若干高速区和低速区 ,各个区之间不必连续 ;3在每一速度区内分别采用不同的参考速度 ,应用裂步傅里叶方法进行波场延拓 ,即先在频率—波数域中做相位移 ,然后再回到频率—空间域逐点实现裂步法 ;4对分区延拓后得到的不同波场 ,按照实际对应的网格点进行合并 ,合并波场时应做适当的光滑处理 ,即可得到该层延拓后的完整波场 ;5重复以上步骤 ,直至完成所有速度分区的波场延拓及合并。理论分析和简单模型的叠后深度偏移试算及 Marmousi模型叠前深度偏移的试算结果表明了该方法的可行性。
- 张叔伦王昌龙赵景霞孙沛勇
- 关键词:叠前深度偏移波场延拓MARMOUSI模型成像精度步长
- 平面波最大能量叠前深度偏移被引量:8
- 2000年
- 本文提出了一种新的基于最大能量震源波场的平面波叠前深度偏移方法。该方法保留了平面波偏移方法的优点 ,即在成像过程中 ,通过同时延拓频率域稀疏采样的平面波震源和最大能量旅行时 ,求取震源波场中的最大能量部分 ,从而显著地提高了平面波叠前深度偏移的效率。利用震源波场中的最大能量部分成像 ,还保证了较高的成像精度。 Marmousi模型的试算结果表明 ,在复杂构造情况下 。
- 孙沛勇张叔伦
- 关键词:叠前深度偏移平面波地震勘探地震数据处理
- 最小二乘傅立叶有限差分偏移被引量:27
- 2008年
- 一般偏移算法是用反演算子通过解析方法求解.最小二乘偏移方法采用另一种思路,即采用数值方法,通过解一个线性离散反问题来索求解.这样我们试着寻找一个模型匹配地震数据并能表现出其某些特点来得到偏移图像.最小二乘法能减少偏移赝像,得到更精确的偏移效果.Kirchhoff算子在最小二乘偏移方法中应用较广,但需要较多的迭代次数,而且具有Kirchhoff偏移的缺点.本文把最小二乘方法运用到基于波长延拓的波动方程偏移方法中,为提高最小二乘偏移的效率,可采用效率较高的正传播算子和反传播算子.我们利用效率较高,能适应剧烈横向变速的傅立叶有限差分正传播和反传播算子来做叠后最小二乘偏移.数值实例表明,通过少数的共轭梯度法迭代,就能得到与真实模型差别很微小的偏移效果.对于傅式变换我们采用了数值软件FFTW,其变换速度比常规FFT算法一般要快六倍以上,进一步提高了效率.本文算法很容易在并行机上实现,这些特点在处理大型数据时大有裨益.
- 杨其强张叔伦
- 关键词:共轭梯度法偏移算子
- 相位移方法的两个新算法
- 1999年
- 文中提出了有关相位移方法的两个新算法。一个是相位移分层算法,无论是模型还是偏移都可以实现剧烈纵横向变速;另一个是共炮检距剖面的相位移算法,运用几何地震学的绕射点公式和平面波分析方法,推出了常速介质条件下,共炮检剖面的模型和偏移的相位移公式。两种新的相位移算法都保留了原相位移方法的优点,合成记录的列子,进一步验证了方法、原理的正确性。
- 张叔伦
- 关键词:平面波地震勘探
- 逆时偏移的方向导数算法被引量:1
- 1991年
- 本文从实用观点出发,在延时坐标下提出了一种新的逆时偏移算法。该算法运用方向导数,把二维标量波动方程变换成形式上与Claerbout的15°方程形式完全相同的方程,从而很容易在计算机上求解。该算法无倾角限制,计算区域与15°方程的要求相同,因而计算量也基本相当。
- 张叔伦陈静张关泉
- 关键词:逆时偏移方向导数
- 起伏地表下合成震源记录剩余偏移速度分析
- 2009年
- 提出一种起伏地表下新的剩余偏移速度分析方法.应用控制照明技术,克服了由于在起伏地表上合成的震源波场不是平面波震源波场而不能应用剩余速度校正公式的障碍,使得射线参数在偏移和速度分析过程中保持一致,且避免了因横向变速而导致的平面波震源波场在传播过程中的畸变.利用改进的合成震源记录叠前深度偏移,直接从起伏地表延拓波场和成像,偏移前无需对地震数据做静校正或重定基准面.合成数据试算表明该方法是有效和实用的.
- 王昌龙张叔伦
- 关键词:起伏地表
- 自适应空间分区裂步傅立叶偏移方法(英文)
- 2005年
- 本文提出了一种新的偏移方法——自适应空间分区裂步傅立叶(ASDSSF)偏移方法。该方法将剥层相位移方法的思想推广到裂步傅立叶偏移方法,使之当速度场出现强间断时也能精确而有效地成像。原理上ASDSSF偏移属于多参考慢度(MRS)偏移方法,本文的重点是,在不损失精度的同时选取比同类MRS偏移方法更少的参考慢度。我们根据全局速度函数的变化和误差控制参数来选择参考慢度,同一个参考慢度所对应的速度构成一个分区,每一个空间分区可以由几个空间上不连续的子分区组成,从而有效地减小了参考慢度的个数。每一延拓步的参考慢度以及参考慢度的个数和如何构建分区都是根据速度函数自动生成,因此更为合理。为了消除速度场强间断产生的人为噪音,设计了简单有效的f-k域的光滑滤波。我们对一个生成的二维叠前模型和SEG/EAEG盐丘模型进行了试算。
- 赵景霞张叔伦王昌龙倪逸
- 相位移的分层模型算法
- 1996年
- 相位移的分层模型算法是以常途相位移算法为基本组成部分,运用非齐次单程波声波方程,把每一层的层速度作为常速,分层计算波场,即把下一层波场作为上一层的震源函数计算上一层的波场,如此由深至浅向上延拓,直到地面为止。这种方法既克服了常速相移法因速度突变而产生的误差,又保持了常速相移法精度高的优点,非常适合于速度纵横向突变的情况。文中给出的合成记录说明了该方法的可行性。
- 张叔伦赵达壮
- 关键词:波场延拓
- 曲网络伪谱法二维声波模拟被引量:16
- 2003年
- 为更精确地模拟地质界面和波的传播 ,给出了一种曲界面波场模拟算法。提出了具体的计算方案 ,并改进了快速傅氏变换的计算。为使用方便 ,推导了稳定性条件。该方法利用多块映射和超限插值技术将直角坐标系下的曲界面变换成曲坐标系下的规则界面 ,在曲坐标系下利用伪谱法模拟波场。它能有效地消除因离散化曲界面产生的人为阶梯状界面所带来的不必要的绕射。直网格下的合成记录和曲网格下的合成记录的比较 。
- 赵景霞张叔伦孙沛勇
- 关键词:伪谱法波场模拟地质界面波声波