深层地震

第17卷 第4期 地 球 物 理 学 进 展 Vol.17 No.4 2002年12月(564~574) PROGRESS IN GEOPHYSICS Dec.2002深层地震勘探的地震波传播理论研究前景田春志1 刘 洪2(大庆石油学院1,安达151400;中国科学院地质与地球物理研究所2,北京100101)[摘 要] 深层地震勘探为地震波传播理论研究提出了新的挑战和机遇.深层地震勘探的主要难点是上覆层的影响甚大,使后续的处理有隔靴挠痒之感,必须应用波场延拓消除上覆层影响.深层波速的高速性和横向不均匀性决定了大角散射和弹性波处理方法的重要性.本文具体评述了深层地震勘探的主要方法对策,深入探讨了波场延拓的李群方法和弹性反演的某些问题,目的在于为深化深层地震勘探提供新的研究手段和方法.[关键词] 深层地震勘探;波场延拓;李群方法;弹性反演[中图分类号] P618 [文献标识码] A [文章编号] 1004 2903(2002)04 0564 110 引 言深层天然气是今后天然气勘探的重要领域[1],是实现并完成天然气储量任务的重要保障,深层天然气勘探技术的深入研究对松辽盆地油气勘探与开发具有重要的现实意义.深层特殊岩性储层的形成是松辽盆地由裂谷盆地向凹陷盆地发展过程中形成的产物,因地质年代相对久远,地层岩性多变,加之多期不同的构造运动的叠合,使特殊岩性储层的构造格局十分复杂,其中登娄库组储层具有高地温场、强成岩作用、弱反射系数的特征,侏罗系具有多期火山活动,多种类型储层,构造起伏剧烈,反射连续性差的特点,在利用目前地震资料进行地质解释研究时,往往存在着构造解释多解性强,构造、岩性圈闭识别精度低的问题,针对深层的构造复杂,砂砾岩、火山岩储层具强非均匀特征,地震波能量弱、分辨能力低、信噪比低特点,在地震资料的采集、处理、解释等方面仍有许多理论问题尚需深入研究.针对地震波能量弱,可以考虑扩大面元叠加[2].但由于走时曲线非双曲线性质[3],直接叠加会损失信号.一种可能的考虑是,通过剩余动校正,将非双曲线校正为双曲线再叠加,但剩余动校正不能校正焦散走时曲线,剩余动校正扩大面元叠加有可能损失复杂构造引起的绕射波,从而影响复杂构造断点的成像.信噪比低特点,需要考虑消除规则干扰和随机噪声的方法.规则干扰包括面波、多次波[4].随机噪声主要表现为野值、高斯型随机噪声、非高斯型随机噪声.随机噪声的消除办法,主要靠叠加.束叠加是一个重要探索方向[5],由此引发对边缘小波、曲线小波的探索[6,7].分辨能力低主要考虑粘弹性损耗对波形的补偿问题,还需要考虑速度频散造成的迭加损失分辨率问题[8,9].为了考虑强非均匀储层特征的刻画,需要考虑非水平地层剥蚀不整合面的散射,这需要考虑对称轴非垂直的横向各向同性介质的反射和透射[10 12].为了考虑强非均匀储层特征的刻画,还需要考虑近源沉积造成的相带快速变化,即地震波的振幅有较强的横向变化.[收稿日期] 2002 05 20; [修回日期] 2002 09 30.[基金项目] 中国科学院知识创新工程重大项目(KZCXI y01),国家自然科学基金和大庆石油管理局项目(49894190)资助.[作者简介] 田春志,1963年3月生,男,黑龙江省庆安县人,1985年长春地质学院获学士学位,现为大庆石油学院博士研究生,从事油气勘探开发研究.波动方程基准面通过波动方程波场深度延拓,可将炮点和接收点都延拓到地下某一深度的平面上.这样可以校正上覆层的透镜效应和对走时的影响,使目标层的多值走时曲线能展开成单值曲线并接近双曲线,振幅更逼真.同时可大幅度减少后续处理和成象的工作量.波动方程基准面方法从根本上消除了上覆层速度横向不均匀的影响,使依赖近双曲线的分频叠加和去噪技术有可能用于深层资料.因此,为开展深层高分辨率三维地震勘探方法研究,需要进行理论研究,其核心关键技术是波动方程基准面和叠前反演技术,其基础是深化地震波理论研究.本文评述了通过李群方法研究深化地震波理论研究的研究前景.1 深层地震勘探方法国内外发展现状美国科罗拉多矿业学院地震波场中心(CWP)近年着重研究复杂构造的速度 深度模型建立、地震波场的正演模拟、叠前深度偏移和反演方法,在以层析成像技术为基础的速度分析方面取得了重要进展.斯坦福大学地球物理项目组主要集中于叠前去噪和波动方程叠前深度偏移,分别在信号保真的噪音压制预测方法和方位角域叠前深度偏移取得重要研究进展.法国石油研究院重点是利用反射波层析成像方法确定速度模型,在无需层位拾取的速度分析方面取得了进展.德国卡尔斯鲁厄大学地球物理研究所提出了共反射面元叠加方法,可大幅度提高地震资料信噪比和分辨率,但用于深层尚有焦散问题待解决.对于裂隙性储层的地球物理预测,美国麻省理工学院地球物理系基于各向异性介质的地震波传播理论,研究了利用叠前信息识别裂隙储层发育条件的理论方法,在振幅、频率、相位随入射角(或炮检距)变化的储层条件分析方面取得进展,并在委内瑞拉油气田勘探中取得效果,这为利用叠前信息识别和预测裂隙储层的研究提供了可借鉴的经验.国内新疆、大港、胜利等油田对火山岩油气藏研究积累了一些资料和经验[10],但在深层地震资料处理和解释方面,尤其在火山岩储层预测、储层描述等方面尚未形成成套的勘探方法和技术.大庆油田在深层方面的研究,盆地内部侏罗系和火山岩分布、层序划分、深层油气成藏机制、地球物理识别等方面取得了重要研究进展,并在实际勘探中取得了重大突破.但因深层地质结构复杂、沉积类型多样、储集介质类型多、气藏控制因素复杂,尚需深入开展储层定量、半定量描述的地质与地球物理综合方法研究.通过 七五!、 八五!、 九五!连续攻关,我国东部主力油田的中浅层高分辨地震勘探已取得丰硕成果,形成了一套行之有效的高分辨率、高信噪比地震勘探采集、处理、解释方法及流程,其关键技术是分频速度估计、分频叠加、反褶积、剩余静校正.将但其应用与深层有焦散问题待解决.若能与波动方程叠前深度延拓技术相结合,则可为更复杂的深层储层建模提供了前期技术保障.在此基础上进行精细储层地质建模,需要大力发展断层检测技术和地层分析.尽管相干数据体技术用于断层检测已引起了广泛关注[13 15],但其用于地层分析效果并不明显[16,17].随着对分辨率、处理速度、地质信息含量的要求越来越高,需要更合理有效的断层检测和地层分析技术[18,19].最近,在国际信号处理领域正在研究噪声为非高斯分布的边缘检测和信号分析技术,提出了独立分量分析[20,21]、Hilbert Huang 变换[22],在石油勘探信号处理领域已取得一定进展.其科学依据是,Nature 杂志1996年曾报告了一种图像处理中的稀疏分量分析方法.其目的是如何用最少的数据表示一个图像,所得到的稀疏分量主要∀565∀4期 田春志,等:深层地震勘探的地震波传播理论研究前景以稀疏线段为主,引起了图像处理、计算机视觉、人工神经网络分析等领域的广泛重视.这个方向的进展,表明稀疏分量分析与小波分析和 p 变换有密切关系,因此人们提出了第二代小波如S 小波、边缘小波(Ridgelet)、楔形小波(Wedgelets)、曲线小波(Curvelet)[6,7]、束形成技术[5]、适合非高斯分布噪声的保边缘去噪技术[18].将上述有关技术与窗式Hilbert 变换结合,在断层检测和地层分析中得到成功应用[18,19].中国科学院地质与地球物理研究所在火山岩喷发期次与储层特征提取等方面提出了小波Hilbert 变换与此原理类似[23 25].同时,从非平稳统计分析发展的独立分量分析,也显示了分解不同特征信号的强大功能,1999年系统总结了Hilbert Huang 变换,已导致了不同特征信号的快速分解算法.在地震波传播理论研究中应用李群和辛群是目前中国科学院地质与地球物理研究所探索的一个方向[26 36].这项工作的背景是用哈密顿体系描述波传播的研究.目前,对应三种等价的波传播动力学体系存在三种离散化算法:牛顿体系 有限差分法;拉格朗日体系 有限元方法;哈密顿体系 辛几何算法.对于前两者,已经有了很多研究,后者的研究国内外刚刚开始.单程波场随深度z 的延拓过程,实际上是波场在一个以z 为参数的单参数李群作用下的变换过程,保持波场演化的李群变换性质的离散化算法为李群算法.辛几何算法是李群算法的一种.研究辛算法和李群算法的最终目的是使地震波偏移成像具有保真性、稳定性和快捷性.他们的主要进展是:(1)发展了基于e 指数的二阶对角Pade 近似辛格式的深度域波场延拓算法,证明了其精度优于两步法的数值特性,为实现波动方程叠前深度偏移提供了计算方法[38,39];(2)在时间域的波传播辛几何算法研究中,导出了波场计算显式、隐式、蛙跳式三类辛格式,在均匀介质和Marmousi 模型上就计算效率、精度、哈密顿量性质进行了比较,说明了辛算法的优越性[35 37];(3)给出了逆时延拓辛几何算法+射线方法的深度偏移方法,计算的青海柴达木盆地北缘地质模型,为合理进行地震采集设计提供了理论依据;(4)在深度域的波场延拓研究中,结合谱法和直接解法的优点,提出了大型稀疏矩阵求逆LU 分解的混合方法,在精度和速度上能更好满足三维隐式李群算法求逆运算的需要.1999年12月用4~6个节点微机集群,在松辽盆地四条二维深大剖面上实现了波动方程叠前深度偏移实际资料处理,不仅证明了波场延拓辛几何算法的数值稳定性和精度,而且揭示了微机集群方案的高效性和低成本特点[38,39].在此基础上,2000年该所与大庆油田联合装配了24节点和32节点微机集群并行机,经过对所研制的辛几何算法处理能力比测,32节点机较SP2(24个CPU)机的运算速度快4倍[40].针对大庆东部深层陆上三维观测系统,完成了炮集、面炮、分裂AMO 等波动方程叠前深度偏移方法效果试验[26,27].表明了波动方程叠前深度延拓方法取得了成功,从而可用于波动方程基准面、波动方程叠前反演方法的研究.为提高分辨率,高精度叠前波场延拓还需考虑粘弹性影响.宋守根推导出考虑速度频散的波动方程深度延拓方法[8,9],该方法以小波理论为基础通过在延拓中利用速度频散和振幅反褶积,实现等效的空变和时变反褶积,并且更符合粘弹性介质波传播规律.同时求出的反褶积和频散参数可做为提取吸收和衰减参数的依据.2 运用李群方法深化地震波理论研究2.1 李群方法概述李群方法在物理学和应用科学中已取得进展.在理论声波和地震波领域,与李群表达密∀566∀ 地 球 物 理 学 进 展 17卷切相关的保辛保结构算法、Weyl Heisenbeger 框架下的步进算法的研究,近年来都取得迅速发展.在光学和雷达方面,关于合成孔径雷达、多普勒扫频、射线光学、高斯束方法与李群、李代数、Weyl Heisenbeger 框架之间的关系也已有了理论的阐述.将水平层状介质下传播矩阵研究的结果,推广到横向非均匀介质,这对旨在于淡化对先验性速度模型要求的条件下,研究多次波成象、多次波消除、面波计算、裂缝带成象等方面的研究开展,均有重要理论意义和实用价值[41 44].在水平层状介质中,频率域的弹性波表达遵循Hamilton Jacobi 方程.其求解已有成熟的传播矩阵方法和正反向传播方法,其在数学上对应于初射法和LU 分解法.在某些情况下,如求面波频散函数时,由于倏逝波的存在,传播矩阵方法会存在数值不稳定性.通过本征向量矩阵元素组成的反对称矩阵,构成了高效准确的递推关系,面波频散函数也得到了较精确的计算[4].从李群角度看,传播矩阵算法的递推关系与李代数有密切的关系.多方向分裂在三维声波单程波算子的隐格式构造和元胞自动机研究中都起到重要作用.将这些研究推广到弹性波三维深度延拓中,借助于李群方法的成果,将起到巨大的推进作用.在水平层状介质条件下,弹性波的深度延拓可归结为如下Hamilton Jacobi 方程f ~z=A f ~.(1)其中f ~为位移应力向量,A 由广义虎克定律、牛顿运动方程导出,此处称之为Hamilton Jac obi 矩阵.为便于理解,给出两个具体实例[45].在SH 波情况下A =0 (z )-1(z )k 2- !2,(2)在P SV 波情况下A =k-1-k ∀(z )[∀(z )+2 (z )]-100[∀(z )+2 (z )]-14 (z )[∀(z )+ (z )][∀(z )+2 (z )]00k ∀(z )[∀(z )+2 (z )]-1-!2(z )-k,(3)通过矩阵的指数函数式求解(1)式:f ~(z +#z )=e A #z f ~(z ),(4)其中e A #z 为传播矩阵.在传播矩阵的求解中,A 的特征分解或对角化十分重要,有A =MEM -1.(5)M 为合成矩阵,其物理意义是将特征波合成为位移应力矢量,其逆为分解矩阵,即将位移应力矢量用特征波展开.在横向非均匀介质中,如果弹性参数横向变化不大,仿照声波方程时的分裂步相移方法[46]-A 0(0,0,z ))+(6)∀567∀4期 田春志,等:深层地震勘探的地震波传播理论研究前景其中e (A (x ,y ,0,0,z)-A 0(0,0,z))#z可在空间域假定波为垂直入射完成,e A 0x ,y ,z#z在波数域分为各个模式正反向部分波后完成.但如果需要更高精度,则需要研究声波方程单程波的各种近似方法[4649]与李群方法[50,51]相结合.(2)、(3)式给出的A 有一个共同的性质tr A =0,(7)从而有 det(exp (A #z ))=exp (tr A #z )=1,(8)说明exp(A #z )#SL (n),A #sl (n),SL (n)为行列式为1的n ∃n 矩阵组成的特殊线性群,sl (n)迹为零的n ∃n 矩阵组成的特殊线性代数.SL (n )是一种李群,sl (n)是其李代数.李群方法的应用包括:(1)波算子的多方向分裂;(2)李代数与递推计算;(3)算子解耦的李代数方法.2.2 波算子的多方向分裂Hamilton Jacobi 矩阵特征分解是求传播矩阵的重要方法.在推广到横向非均匀介质时,要求解单程波方程.已导出了一系列显格式和隐格式方法.单程波算子的多方向分裂是导出隐格式的重要手段.根据元胞自动机研究,这和波的粒子表示是等价的.多方向分裂公式为f (k x ,k y )=12∃%∃g(k x cos (%)+k y sin (%),%)d %.(9)(9)式可通过Radon 变换的反投影定理证明g(u ,%)=1∃ p (u,%) u *P f 1u,(10)p (u,%)=%&-&f (u cos (%)-v sin (%),u sin (%)+v cos (%))d v .(11)若,f (k x ,k y )=f 1(k 2x +k 2y ),则(10)和(11)式还可简化,例如f 1(k 2)=k 2a 2-k2,(12)p (u,%)=p (u)=i ∃2a 2u (a 2-u 2)32,(13)g(u,%)=g (u)=%&-&i ∃2a 2v (a 2-u 2)3dvu -v =12a 2v(a 2-v 2)32+%&0∃at (1+t 2)32dt u -iat .(14)2.3 李代数与递推计算面波频散函数是分解矩阵的反对称子矩阵经正向传播和反向传播递推计算得到.当考虑大角度入射时,例如深部有高速层,由于倏逝波的存在,用传播矩阵方法计算面波频散函数会存在数值不稳定性.但如果将正向传播和反向传播同时计算,则有可能得到较高精度[4].设K L,1为从1层面到L 层面的传播矩阵(即exp (%zlz1A d z )),Y L +1为L 层面以下分解矩阵M 的子矩阵构成的一个二次反对称子矩阵,有∀568∀ 地 球 物 理 学 进 展 17卷Y L+1=M T L+1,2M L+1,4-M TL+1,4M L+1,2,(15)D L ,1=K T L,1Y L+1K L,1.(16)则面波频散函数是矩阵D L,1的一个元素.设D L,L =K TL,L-1Y L+1K L,L-1,(7)则利用 D L,i =K Ti+1,i D L,i+1K i +1,i ,(18)从i =L -1到i =1递推,可求得D L ,1.注意到D L,i 与D L ,i +1是线性依赖的关系,其线性变换系数可以解析推出,从而可提高计算精度.在地震成象领域,有很多变换是由正反变换序列实施组合而成,如变换到零炮检距(MZO)定义为非零炮检距偏移与零炮检距反偏移的组合.又如,在有限孔径地震波成象中,畸变函数是由孔径函数经正向传播和反向传播形成的.剩余偏移是某一速度反偏移和和另一速度正偏移的组合而成.利用李群方法将同时计算的正向传播和反向传播加以简化是十分重要的.从李群方法的角度看,递推关系与李代数有密切关系.当层厚趋于零时,(18)式就可化为新的Ha milton Jacobi 方程,其新的李代数为dd z(exp(d z A T )exp(d z A ))|d z=0.(19)用几何方法已导出MZO 和剩余偏移的李代数[53 55],孔径畸变函数的李代数也可导出.2.4 算子解耦的李代数方法深度延拓算子或传播矩阵属于一个行列式为1的特殊线性群.在传播区内,深度延拓算子或传播矩阵用辛几何算法是合适的.在传播区外,辛几何算法是不稳定的[52].在单平方根算子情形,传播区边缘的算子逼近非常困难,用非辛的显格式,算子较长;切比雪夫、台劳展开也往往不能保持精度.在物理上,在传播区的边缘,意味着波沿水平方向传播,出现焦散,通常用WKBJ 或Maslov 方法处理.李群方法的基本思想是,乘积分解后李代数不变.行列式为1的特殊线性群应表示为基本特殊线性群的乘积,基本特殊线性群应表示为李代数基的指数映射.这样一来就从保辛的思想进一步深入到保李代数.保李代数的想法不仅适应于辛群,而且适应更广泛的特殊线性群.保李代数比保辛要求更高.设G 是一个李群,g 是其相应的d 维李代数.我们选择一个g 的基{V 1,V 2,V 3,∋,V d }.这样对每个充分靠近I 的元素Y #G ,可以表示为Y =exp (&1Y 1)exp (&2V 2)exp (&3V 3)∋exp (&d V d ).(20)这里exp:g (G 是指数函数映射.在数值计算上,经常考虑g gl (n )G GL (n),此时,e xp为通常的矩阵指数函数.给定B #g ,可唯一地表示为B =)di=1∋i V i .(21)有可能将exp(t B )写成U(t)=exp (t B )=exp (g 1(t)V 1)exp (g 2(t)V 2)exp (g 3(t)V 3)∋exp (g d (t)V d ).(22)除某些简单的低维李代数外,g i (t)的形式一般比较复杂.Celledini and Damtp [50]研究了用多项式(i (t)近似表示g i (t)的方法,即求满足exp (t B ∗exp ((1(t)V 1)exp ((2(t)V 2)exp ((3(t)V 3)∋exp ((d (t)V d )的(i (t).并给出了李代数的结构常数求(i (t)的方法.∀569∀4期 田春志,等:深层地震勘探的地震波传播理论研究前景Lorenz 群是一个6维群SO(3,1),它是由所有满足AJA T =J 的4∃4矩阵A 组成,其中j =diag(1,1,1,-1).其李代数为Lorenz 代数so (3,1),它是由所有满足BJ +JB T =O 的4∃4矩阵B 组成:B =∋1∋2∋3-∋10∋4∋5-∋2-∋40∋6∋3∋5∋6.(23)李代数基为{V 1,V 2,V 3,V 4,V 5,V 6}=0100-10000000000,00100000-1000000,00000100-100000,0001000000001000,0000000100000100,0000000000010010.(24)其交换算子的表为[V 1,V 2]=-V 3,[V 1,V 3]=V 2,[V 1,V 4]=-V 5,[V 1,V 5]=-V 4,[V 1,V 6]=O,[V 2,V 3]=-V 1,[V 2,V 4]=-V 6,[V 2,V 5]=O,[V 2,V 6]=V 4,[V 3,V 4]=O,[V 3,V 5]=-V 6,[V 3,V 6]=V 5,[V 4,V 5]=V 1,[V 4,V 6]=V 2,[V 5,V 6]=V 3.利用交换算子的表中的结构常数,利用Celledini and Iserles A [51]的方法可以导出(1(t)=∋1t +12(∋2∋3-∋4∋5)t 2,(2(t)=∋2t -12(∋1∋3+∋4∋6)t 2,(3(t)=∋3t +12(∋1∋2-∋5∋6)t 2,(4(t)=∋4t -12(∋1∋5+∋2∋6)t 2,(5(t)=∋5t +12(∋1∋4-∋3∋6)t 2,(6(t)=∋6t +12(∋1∋53-∋2∋4)t 2.(25)此理论对单程波算子的多方向分裂、DMO 和MZO 理论中指数算子的构造十分有益[54 56].可见弹性波的研究与李群方法密切相关,加强李群方法的研究有重要意义.3 弹性波反演方法考虑一个各向同性的大地,x 表示大地内的一点,而!为时间变量, 0(x )、∀(x )、 0(x )∀570∀ 地 球 物 理 学 进 展 17卷分别为背景的密度和拉梅常数,有(x )= 0(x )+) 0(x ),∀(x )=∀0(x )+)∀0(x ), (x )= 0(x )+)0(x ).(26)u i (x ,!,x s )=u i 0(x ,!,x s )+)u i (x ,!,x s ).(27)为了避免问题所固有的病态,重要的是识别独立参数,根据物理上的证据I p (x )= (x )((x )=(x )[∀(x )+2 (x )],(28)I s (x )= (x )∋(x )=(x ) (x ).(29)纵波和横波波阻抗是独立参量.弹性参量和它们的换算关系是∀(x )=1 x(I 2p (x )-2I 2s (x )),(30) (x )=1 x I 2s (x ),(31))∀(x )=-((2(x )-2∋2(x ))) (x )+2((x ))I p -4∋(x ))I s (x ),(32)) (x )=-∋2(x )) (x )+2∋(x ))I s (x ).(33)由于介质扰动产生的频率域散射场可表为[56])u =A 0)I p +B 0)I s +C 0) =(A 0 B 0 C 0))I p)I s).(34)其中:A i0)I p =%d x A i 0(x r ,x s ,!,x ))I p (x ),(35)B i 0)I p =%d x B i(x r ,x s ,!,x ))I s (x ),(36)C i 0)I p=%d x C i 0(x r ,x s ,!,x )) (x ),(37)A i(x r ,!,x s |x )=-2((x ) ∗i j x j(x r ,!|x )∗mm(x s ,!|x ),(38)B i (x r ,!,x s |x )=-4∋(x ) ∗i j x m (x r ,!|x )∗jm (x s ,!|x )- ∗i j x j(x r ,!|x )∗mm(x s ,!|x ),(39)C i (x r ,!,x s |x )=-!2∗i j (x r ,!|x )∗j (x s ,!)+((2(x )-2∋2(x )) ∗i j x j (x r ,!|x )∗m m (x s ,!|x )+2∋2(x ) ∗ij x m(x r ,!|x )∗jm(x s ,!|x )(40)对应的弹性波反传公式为)u ,(41)∀571∀4期 田春志,等:深层地震勘探的地震波传播理论研究前景对应的反演公式为)I p)I s)=A T 0B T 0C T(A 0 B 0 C 0))I p )I s )=A T 0A 0)I pB T 0B 0)I sC T0C 0).(42)这是波阻抗反演在深度域和弹性波的推广[57].在高分辨率和高保真的要求下,(42)式是多解的,需要引进先验模型约束.先验模型的计算机生成的思路是,在构造模式、地层时空配置及叠置关系约束下,结合地质知识、测井解释、地震资料建立储层地层厚度,分布范围及几何学形态的统计特征和储层发育的统计规律,借助于随机建模方法和虚拟井技术,用以快速地建立储层模型.(42)的反演仅能反映的储层的弹性参数信息,而用以描述储层参数主要是包括孔隙度、渗透率、含油气性.为利用地震资料描述或预测储层特征参数或借助于正演模拟技术进行储层条件地震波场的影响特征研究,需要根据岩石物性理论和技术搞清储层参数与弹性参数张量之间转换关系,利用分析力学原理和张量反演方法,建立储层参数与弹性参数张量的转换关系,进而利用测井资料、岩石物性测量数据,建立储层参数与弹性参数之间的统计规律,并进行利用弹性参数预测储层条件及含油气性精度与可靠度的定量分析研究.参考文献[1] 张颖,余云英,陈立康.东部油区深层勘探中的地震处理技术[J].勘探家,1999,4(2):50~52.[2] H +ocht G,de Bazelaire E,Maj er P,Hubral P.Seismic andoptics :hyperbola and curvature[J].J.Appl.Geophys.,1999,42:261~281.[3] 高亮,李幼铭,陈旭荣,杨孔庆.地震射线辛几何算法初探[J].地球物理学报,2000,43(3):402~410.[4] 李幼铭,束沛镒.层状介质中地震面波频散函数和体波广义反射系数的计算[J].地球物理学报,1982,25(2):130~139.[5] Hu T,Whi te R E.Robus t multiple suppression usi ng adaptive beamforming [J ].Geophysical Prospec ting,1998,46:227~248.[6] Donoho D L.Wedgelets:nearly mini max es ti mation of edges[J].Ann.Statist,1999,27:859~897.[7] Chen S,Donoho D L,Saunders M A.Atomic Dec omposi ti onby Basis Purs uit[J].SIAM p.,1998,20(1):33~61.[8] Shougen Song,Rongfeng Zhang,Taseusz J Ulrych.Highres ol ution wavefield reconstruction i n s ei smic mi gration for viscoelas tic medi a [J ].Journal of Seis mic Exploration ,2000,9:1~18.[9] 王志章,石占中.现代油藏描述技术[M].北京:石油工业出版社,1999,163~176.[10] Bakulin A,Grechka V,Tsvankin I.Estimati on of frac tureparameters from reflec tion seis mic data Part I:HTI model due to a single fracture set[J].Geophysics,2000,65(6):1788~1802.[11] Bakulin A,Grechka V,Tsvankin I.Estimati on of frac tureparameters fro m reflection s eis mic data Part II:Fratured models with orthorhombic symmetry [J].Geophysics,2000,65(6):1803~1817.[12] Bakulin A,Grechka V,Tsvankin I.Estimati on of frac tureparameters from reflecti on seis mic data Part III:Fratured models with monoclinic symmetry [J ].Geophys ics,2000,65(6):1788~1802.[13] Mike Bahoric h,Steve Farmer,3D seismie di scontinui ty forfaults and stratigraphic features the coherence cube[J].The Leading Edge,1995,14(10):[14] Kurt J Marfurt,Ronald M Scheet,J ohn A Sharp,Mark GHarper.Suppress of the acquisition footprint for s eis mic se quence attribute mapping[J].Geophysics,1998,63(3):1024~1035.[15] Kurt J Marfurt,RLynn kirlin,Steven L Farmer,M ichael S∀572∀ 地 球 物 理 学 进 展 17卷Bahorich.3 D s eis mic attributes usi ng a semblance basedcoherency al gori thm[J].Geophysics,1998,63(4):1150~1165.[16] 李玲等.用地震相干数据体进行断层自动解释[J].石油地球物理勘探,1998,33(增刊1):105~111. [17] 刘传虎.地震相干分析技术在裂缝油气藏预测中的应用[J].石油地球物理勘探,2001,36(2):238~244.[18] Y Luo,W G Higgs,W S Kowalik.Edge detecti on ands tratigraphic analysis using3D seis mic data[A].Soc.Expl.Geophys.69th Ann Internat.Mtg Abs tract[C].1996,324~325.[19] N Carter,L Lines.Fault i maging using edge detecti on andcoherency meas ures on Hibernia3 D seismic data[J].TheLeading Edge,2001,20(1):64~69.[20] Bruno A,Ols haus en,David J Field.Emergence of si mplecell receptive field properties by learning a s pars e code fornatural i maging[J].Nature,1996,381:607~609. [21] A Hyvarinen.Fas t and robus t fixe-point algori thms forindependent component analysis[J].IEEE Trans.O n Neural Networks,1999,10(3):626~634.[22] N E Huang,Z Shen,S R Long,et al.The empirical modedecomposi ti on and the Hilbert s pectrum for nonlinear andnon stati onary time s eries analysis[J].Proc R Soc Lond,1998,454:903~995.[23] Gao Jing huai,Dong Xiao long,WangWen bing,Li Youmi ng,Pan Cun huan.Instantaneous Parameters extracti onvia wavelet transform[J].IEEE Tran.On Geoscience andRemote Sensing,,1999,37:867~870.[24] 高静怀,等.小波变换与信号瞬时特征分析[J].地球物理学报,1997,40(6):821~832.[25] 王西文,等.地震信号瞬时特征在小波域分频提取的方法和应用[J].石油地球物理勘探,2000,35(4):452~460,478.[26] 刘礼农,刘洪,李幼铭.波动方程辛几何算法三维叠前深度偏移流程及应用[J].中国科学院研究生院学报,2002,19(2):134~139.[27] 陈树民,刘洪,李幼铭.适于窄线三维地震资料条件的面炮方法[J].地球物理学进展,2001,16(3):85~94.[28] 罗明秋,刘洪,李幼铭.适用于波场成像的谱法LU分解[J].地球物理学报(待刊).[29] 刘洪,李幼铭.勘探地球物理学家协会第70届年会简况[J].地球物理学进展,2001,16(1):110~113. [30] 高红伟,李幼铭,刘洪.图法及其在Toeplitz矩阵分解中的应用[J].地球物理学进展,2001,16(4):35~42.[31] 罗明秋,刘洪,李幼铭.地震波传播的哈密顿表述及辛几何算法[J].地球物理学报,2001,44(1):120~127.[32] 罗明秋,刘洪,李幼铭.基于螺旋线上谱因式分解的地震波场隐式辛算法[J].地球物理学报,2001,44(3):379~388.[33] 李幼铭,高亮,刘洪,杨孔庆,曹禹.地震波传播的Hamilton描述及模拟[A].见:反射地震学论文集[C].同济大学海洋地质与地球物理系编,上海:同济大学出版社,2000,116~120.[34] 杨孔庆,曹禹,李幼铭.二维波场的几何量子化和单向传播[A].反射地震学论文集[C].同济大学海洋地质与地球物理系编,上海:同济大学出版社,2000,276~279.[35] Luo Mingqiu,Gao Hongwei,Liu Hong,Li Youming.Acomparison on the LU decomposition of Laplacian matri xwith spectral and graphic methods[A].In:Soc.Expl.Geophys.70th Ann Internat.Mtg Abstract[C].2000,2349~2352.[36] Luo M i ngqiu,Li Youming,Li u Hong.The Sympleetic Geometric Deseri ption and Algortthm of Sei smic Wave Propagati on[A].In:Soc.Expl.Geophys.69th Ann Internat.Mtg Abstract[C].1999,1825~1828.[37] 杨辉,刘洪,李幼铭.一种单程波方程算子的辛格式描述[A].见:反射地震学论文集[C].同济大学海洋地质与地球物理系编,上海:同济大学出版社,298~302,2000.[38] Yang Hui,Li u Hong,Li Youming,Wang Chunmei,BiMin.The Application of Symplec tic Depth Prestack Migration in Daqing XJWZ Area[A].In:the70th Ann,SEG,meeting[C].2000,882~885.[39] 杨辉,高亮,刘洪,李幼铭,范兴才.微机群并行实现M armousi模型叠前深度偏移[J].地球物理学进展,2001,16(3):58~75.[40] Zhang Jianfeng,Li Youming.Numerical simulation of elastic wave propagation in inhomogeneous medi a[J].WaveMotion,1997,25:109~125.[41] A B Weglein,K H Matson,D J Foster,P M Carvalho,DCorri gan,S A Shaw.Imaging and i nversi on at depth without a velocity model:theory,concepts and initial evaluation[A].In:Soc Expl Geophys.70th Ann Internat.MtgAbs tract[C].2000,1016~1019.[42] Vers chuur D J,Berkhout A J,Wapenaar C P A.Adapti vesurface related multi ple elimination[J].Geophysics,1992,57:1166~1177.∀573∀4期 田春志,等:深层地震勘探的地震波传播理论研究前景。