• / 10
  • 下载费用:5 下载币  

地震叠前逆时偏移中的去噪与存储

关 键 词:
物探 地震资料解释 地震处理 反演
资源描述:
书书书第53卷第9期2010年9月地球物理学报CHINESE JOURNAL OF GEOPHYSICSVol.53,No.9Sep.,2010刘红伟,刘洪,邹振.地震叠前逆时偏移中的去噪与存储.地球物理学报,2010,53(9):2171~2180,DOI:10.3969/j.issn.00015733.2010.09.017LiuHW,LiuH,ZouZ.Theproblemsofdenoiseandstorageinseismicreversetimemigration.犆犺犻狀犲狊犲犑.犌犲狅狆犺狔狊.(inChinese),2010,53(9):2171~2180,DOI:10.3969/j.issn.00015733.2010.09.017地震叠前逆时偏移中的去噪与存储刘红伟1,2,刘洪1,邹振1,2,崔永福31中国科学院地质与地球物理研究所中国科学院油气资源研究重点实验室,北京1000292中国科学院研究生院,北京1000393中国石油塔里木油田分公司勘探开发研究院,新疆库尔勒841000摘要地震叠前逆时偏移是当前公认的地震成像的有效途径,然而它面临着计算量甚巨,低频成像噪音以及存储量大等问题,因此,业内科研工作者对其研究乐此不疲.借助GPU/CPU协同计算可以有效解决计算量的难点,笔者已在另文中阐述,本文着重探讨成像噪音抑制以及存储问题.文中分析了叠前逆时偏移产生成像噪音的机制,据此提出在叠前地震资料中先对数据进行相位与振幅校正,进而在成像后运用拉普拉斯算子滤波法消除成像噪音,从而有效去除成像所产生的低频噪音;针对存储量,采用随机边界,用计算换存储,并借助GPU实现,节省了GPU与CPU之间的数据通讯,数值实验结果表明,采用随机边界方法的逆时偏移结果与直接存储波场的方法得到的结果差别甚小.关键词低频成像噪音,拉普拉斯滤波器,随机边界,数据通讯,GPU实现DOI:10.3969/j.issn.00015733.2010.09.017中图分类号P631收稿日期20100127,20100904收修定稿基金项目国家科技重大专项(2008ZX05008006003)与国家重大科研装备研制项目(ZDYZ20082A08)联合资助.作者简介刘红伟,男,1982年生,中国科学院地质与地球物理研究所在站博士后,主要从事地震成像与高性能计算方面的研究.Email:liuhw@mail.iggcas.ac.cn犜犺犲狆狉狅犫犾犲犿狊狅犳犱犲狀狅犻狊犲犪狀犱狊狋狅狉犪犵犲犻狀狊犲犻狊犿犻犮狉犲狏犲狉狊犲狋犻犿犲犿犻犵狉犪狋犻狅狀LIUHongWei1,2,LIUHong1,ZOUZhen1,2,CUIYongFu31犓犲狔犔犪犫狅狉犪狋狅狉狔狅犳犘犲狋狉狅犾犲狌犿犚犲狊狅狌狉犮犲狊犚犲狊犲犪狉犮犺,犐狀狊狋犻狋狌狋犲狅犳犌犲狅犾狅犵狔犪狀犱犌犲狅狆犺狔狊犻犮狊,犆犺犻狀犲狊犲犃犮犪犱犲犿狔狅犳犛犮犻犲狀犮犲狊,犅犲犻犼犻狀犵100029,犆犺犻狀犪2犌狉犪犱狌犪狋犲犝狀犻狏犲狉狊犻狋狔,狋犺犲犆犺犻狀犲狊犲犃犮犪犱犲犿狔狅犳犛犮犻犲狀犮犲狊,犅犲犻犼犻狀犵100049,犆犺犻狀犪3犚犲狊犲犪狉犮犺犐狀狊狋犻狋狌狋犲狅犳犈狓狆犾狅狉犪狋犻狅狀犪狀犱犇犲狏犲犾狅狆犿犲狀狋,犜犪狉犻犿犗犻犾犳犻犲犾犱犆狅犿狆犪狀狔,犘犲狋狉狅犆犺犻狀犪,犓狅狉犾犪,犡犻狀犼犻犪狀犵841000,犆犺犻狀犪犃犫狊狋狉犪犮狋 Prestackreversetimemigration(RTM)isaveryusefultoolforseismicimaging,however,ithassomeproblemssuchashighlyintensivecomputationcost,imagingnoise,andmassymemorydemand,whichtheresearchersinthefieldhavebeenstrivingtosolve.TheproblemofintensivecomputationcostcouldbesolvedbyGPU/CPUcollaborativecomputingandhasbeendiscussedinthelastpaperandtheothertwoproblemswillbeemphasizedinthispaper.Atthebeginning,weanalyzethegenerationprincipleofimagingnoiseandproposethatthephaseandfrequencyspectraofseismicdatashouldbemodifiedbeforemigrationandaLaplacianfiltercouldbeusedtoremovethelowwavenumbernoiseefficiently.Aimingattheproblemofmassymemorydemand,weadopttherandomboundaryconditionwhichreducesthememorydemandbutsacrificethecomputationcost.TheimplementationcouldbeperformedonGPUandthecommunicationbetweenCPUandGPUcouldbesavedwhichcouldreducethecomputationcostin地球物理学报(ChineseJ.Geophys.)53卷anotherway.Thetestsonsyntheticdataexamplesillustratethatthedifferencebetweenthemigrationresultfromthemethodofrandomboundaryconditionandthatfromthetraditionalmethodofsavingthefootprintofthewavefieldcouldbeoverlooked.犓犲狔狑狅狉犱狊 Lowfrequencyimagingnoise,Laplacianfilter,Randomboundary,Datacommunication,GPUrealization1 引言逆时偏移的历史可以追溯到1978年Hemon[1]的文章.真正将逆时偏移引入地震勘探这一领域要归功于1983年几位科学家Baysal[2]、McMechan[3]、McMechan[4]和Loewenthal[5]等独立的科研工作,所有这些方法最初都是用做叠后偏移.近年来,由于计算机技术的飞速发展以及对复杂构造油气藏成像要求的提高,逆时偏移重新登上了历史的舞台.目前,逆时偏移技术已经从叠后走向叠前,从二维走向三维[6],从声波方程走向弹性波方程[7],从各向同性介质走向TTI(TiltedTransverselyIsotropic,倾斜横向各向同性)介质[8~11].国内方面逆时偏移的研究可追溯到20世纪90年代,这些研究包括有限元法[12~14]、谱元法[15]、有限差分法[16]以及混合法[17].另外,薛东川等[18]采用最大振幅成像条件替代最小走时成像条件;杜启振等[19]将叠前逆时偏移技术应用到横向各向同性介质弹性波多分量数据;刘红伟等[20]用GPU实现了高阶有限差分法逆时偏移,将逆时偏移算法加速了30倍.计算量大、低频成像噪音以及存储量大是制约叠前逆时偏移应用的三大瓶颈,计算量大的问题笔者已经在文献[20]中介绍(GPU/CPU协同计算高效实现),本文主要探讨低频成像噪音消除以及存储问题.低频成像噪音主要是由于强反射界面内反射波造成的,噪音的频率一般较低,一种简单的方法是高通滤波,但这种消除往往是不彻底的,并且会破坏波场的有用信息(如盐丘边界信息).低频成像噪音消除的研究从逆时偏移诞生之日起就已经开始,它也是制约逆时偏移应用的一大瓶颈.消除低频成像噪音主要有三种途径:第一,在波场传播过程中消除.Baysal等[21]采用无反射方程,可以消除垂直入射的内反射,但这种方法只能应用于叠后数据,因为无反射方程对于非垂直入射的波反射系数不为零;Loewenthal等[22]对Baysal的方法做出改进,采用大于波长长度的窗函数对模型慢度做平滑,从而消除内反射,这种方法的效果比Baysal的方法效果要好,但是消除了反射的同时也消除了有用的信息,比如棱柱波等信息,会影响逆时偏移成像质量;RobinFFletcher等[23]采用类似于吸收边界条件的方法,通过在波动方程中引入衰减项,在计算区域产生噪音的位置采用方向衰减的边界条件,压制内反射造成的成像噪音,问题在于,这种方法需要人工交互判断噪音产生的位置,另外衰减边界的引入会衰减掉有用信息.第二,在应用成像条件时去噪.KwangjinYoon等[24]利用噪音产生位置炮点与检波点波场传播方向一致性的特点,采用坡印廷矢量成像条件,在成像条件上乘一个方向衰减因子,去噪效果较好,问题在于需要额外计算和存储波场传播矢量,另外,在波场复杂区域,提取波场传播矢量比较困难;BrunoKaelin等[25]分析了成像噪音与检波点波场的相关性,提出在相关成像条件基础上除检波点的波场,从而达到去噪的目的,这种方法实现起来比较方便,但去噪效果不明显;FaqiLiu等[26]介绍了一种新的成像条件去噪,通过把波场分解成(波数域)上下行波和左右行波(二维情形),运用成像条件时只保留方向相反的波场分量成像结果,可以有效的去噪,问题在于波场分解增加了额外的计算量,对于三维波场需要分解为八个方向分量,另外这种分解可能会损失垂直反射层等有效信息;同年,AndréBulc珘ao等[27]提出了一种类似的波场分离去噪方法,原理与上面的方法类似,面临的问题也相似.第三,成像后滤波法去噪.AntoineGuitton等[28]指出传统的滤波方法改变成像结果相位与频谱信息的缺点,采用预测误差滤波器实现最小二乘滤波,在去除成像噪音的同时,波场反射信息得以保存;Zhang等[29]指出拉普拉斯算子滤波相当于成像波场角度域衰减,并且这种操作不需要输出角度域道集.实际操作中,根据不同的问题应该考虑选用不同的去噪方法,另外,应该尽量避免人工干预,避免增加过多的计算量[30].本文提出在叠前对数据进行相位与振幅校正,在成像后运用拉普拉斯算子滤波法消除成像噪音的方法.波场存储是叠前逆时偏移的另外一个瓶颈.因为炮点波场和检波点波场沿相反的时间方向递推,2712 9期刘红伟等:地震叠前逆时偏移中的去噪与存储而应用互相关成像条件时需要相同时刻的波场.这种矛盾迫使叠前逆时偏移过程中,必须将炮点正传波场或者检波点反传波场的传播历史保存下来.这种波场传播的历史往往需要巨大的存储空间,以500×1000×1000的三维空间网格为例,时间方向的步数对地震勘探而言,往往是10000的量级,要保存整个传播历史,单炮数据偏移大约需要20TB的存储量和通讯量,这是目前的计算机硬件无法承受的,必须采取有效的算法减少存储的需要.RueyChyuanShih等[31]提出了一种分段多步偏移方法,将偏移区域在深度方向按照构造特征划分为几个连续的空间,对每个空间依次偏移;ChristopherEHarris等[32]采用类似的方法将成像空间划分为规则的矩形区域,避免人工干预确定边界的复杂性,但这两种方法需要保存分块偏移子空间边界的波场信息,而且增加了I/O和计算时间,并且这两种方法只用在叠后数据上;HGuan等[33]将这种思路用于叠前偏移,而且对于构造相对简单的空间区域,可以选择Kirchhoff积分或者单程波算法来偏移,进一步减少计算量和存储的需求.这种分段多步偏移方法的缺点在于空间的划分要求炮点与检波点在某一个划分的空间内的传播方向单一,这种要求在速度模型未知的情况下是比较困难的,如果空间划分不够合理,无法得到理想的成像结果.WilliamWSymes等[34]提出设置“检查点”(checkpointing)的最优方法.第一次正传时通过设置某些时刻作为“检查点”将波场存储下来,对于非检查点时刻的波场,通过从特定的检查点重新计算这一时刻的波场,选择合适的检查点,可以降低“重算比”(recomputationratio),问题在于,检查点选得太少,重算比太高,检查点选得过多,所需要存储的临时波场越大,存储问题依然存在;RobertGClapp等[35]提出随机边界的思想,用计算换存储.本文采用随机边界的思路,借助GPU/CPU协同计算并行实现,优势在于避免了GPU和CPU之间的数据通讯,可以部分弥补增加的计算量.本文主要包含以下几个部分:首先分析了低频成像噪音的产生机制,提出在叠前对数据进行相位与振幅校正,在成像后运用拉普拉斯算子滤波法消除成像噪音的方法,通过理论模型的测试证明方法的有效性;然后在总结叠前逆时偏移解决波场存储的五种方法的基础上,提出借助GPU/CPU协同计算实现随机边界法逆时偏移,通过盐丘模型的测试证明了方法的可行性.2 去噪2.1 噪音的产生机制地震叠前逆时偏移低频成像噪音主要是由于强反射界面内反射波造成的.图1a是一个简单的层状模型,第一层厚度1000m,速度为3000m/s.在狓=1500m的位置放炮、3000m位置接收的一道地震记录走时为0.833s.图1b中的椭圆是将这一道地震记录采用单程波方法的偏移结果,椭圆的底部位置是反射界面的位置,将多道和多炮偏移结果相加时相干加强,而椭圆的其他位置相干消除,单程波偏移结果没有低频噪音,这是由于单程波方法限制了波场的传播方向,成像只发生在炮点与检波点波场传播方向相反的位置.与单程波方法不同,逆时偏移会产生低频噪音.根据互相关成像条件:image(狓,狕)=∑time犛(狓,狕,狋)犚(狓,狕,狋),(1)凡是满足成像条件的位置都可以成像,在本文的例子中,凡是炮点正传时间加上检波点反传时间等于地震道走时0.833s的位置都可以成像,如图1a中圆圈和方框表示的位置.实际上,在图中箭头线所示的射线整条路径上,成像条件都是满足的,因此会在整条射线路径上成像,如图1c所示,这就是叠前逆时偏移低频噪音的产生机制.2.2 拉普拉斯算子去噪拉普拉斯算子的波数域表示如下:Δ2→FFT-犽2狓+犽2狔+犽2()狕=-犽I2,(2)(2)式中的犽I是成像域波数矢量,如图2所示,有如下关系式[36]:犽I=犽R-犽S,(3)应用余弦定理可以求得:犽I2=犽R2+犽S2-2犽R2犽S2cos(π-2θ)=4ω2狏2cos2θ,(4)(4)式中θ是入射角,犽R与犽S分别是检波点波场与炮点波场的波数矢量.从(4)式不难看出:将拉普拉斯算子作用于逆时偏移成像结果相当于在角度道集上做角度衰减,并且这种方法并不需要输出角度道集,因为拉普拉斯算子是线性操作.图1a中圆圈位置成像噪音的入射角为90°,根据(4)式,这部分噪音可以被完全消除;方框位置的噪音的入射角度不为90°,只能被部分的消除,如图1d所示,笔者认为,3712地球物理学报(ChineseJ.Geophys.)53卷图1 逆时偏移低频成像噪音产生机制(a)是速度模型,圆圈位置指示了入射角为90°的噪音产生机制,方框位置指示了入射角不为90°的噪音机制;(b)是单程波偏移结果,无噪音;(c)是逆时偏移结果,含低频噪音;(d)是逆时偏移去噪后的结果.Fig.1 ThegenerationmechanismofthelowfrequencybandnoiseinRTM(a)Thevelocitymodel.Thecircledenotesthegenerationmechanismofnoisewheretheincidentangleequalsto90°andthesquaredenotesthegenerationmechanismofnoisewheretheincidentangledoesnotequalto90°;(b)Themigrationresultofonewaywavemethodwithoutimagingnoise;(c)ThemigrationresultofRTMwithseverelowwavenumberimagingnoise;(d)ThemigrationresultofRTMafterdenoisewithLaplacianfilter.这种噪音去除结果是合理的,因为方框位置所示的波场传播路径与棱柱波(Prismaticwave)完全一样,如果这种噪音被完全去除,那么棱柱波成像结果也会被去除,逆时偏移的优势便荡然无存.与棱柱波不同的是,这种噪音会在多道叠加时相消,而棱柱波会相干增强.2.3 叠前处理从(4)式不难看出,拉普拉斯去噪的同时会破坏成像信息的振幅和相位信息.分母的速度平方项比较容易处理,分子中的频率项可以在成像前数据上补偿.补偿方法有两种:一种直观的方法是将数据变换到频率域,乘上滤波因子,再变换回时间域.另一种方法是在时间域积分,利用以下关系:犉(∫犵(狋)d狋)=1iω犉(犵(狋)).(5)时间域的方法操作简单,增加的额外计算量可以忽略不计,更重要的是这种方法可以保持波场的局部性;频率域滤波的方法需要做两次Fourier变换,增加了额外的计算量,并且可能会破坏波场的局部性,这是逆时偏移所不希望的.因为逆时偏移是沿时间方向递推,破坏了波场的局部性意味着引入噪音.图2 成像域波数矢量计算原理Fig.2 Theprincipleofcalculatingtheimagingwavenumbervector图3a是一个标准的线性相位Ricker子波,图3d是它的频谱,图3b和3c分别是时间域一次积分和频率域滤波(只对带限频率操作)后的子波,图3e和3f是二者的频谱图,虽然二者频谱类似,但时间域的子波却不同,频率域滤波后的子波的局部性遭到破坏,对逆时偏移是有害的.2.4 模型测试下面分别对Marmousi模型和Sigsbee2a模型做去噪测试.2.4.1 Marmousi模型图4a是Marmousi偏移速度模型,模型中存在4712 9期刘红伟等:地震叠前逆时偏移中的去噪与存储图3 Ricker子波测试(a)原始Ricker子波;(b)一次积分校正后的子波;(c)滤波运算之后的子波;(d)、(e)和(f)分别是前三种子波对应的频谱.Fig.3 Testonrickerwavelet(a)TheoriginalRickerwavelet;(b)Thewaveletafterintegralmodification;(c)Thewaveletbeingfilteredinfrequencydomain;(d),(e)and(f)arethefrequencyspectralofthethreewavelets,respectively.图4 Marmousi模型去噪测试(a)速度模型;(b)去噪之前的逆时偏移结果;(c)去噪之后的逆时偏移结果.Fig.4 ThedenoisetestonMarmousimodel(a)thevelocitymodel;(b)theresultofRTMbeforedenoising;(c)theresultofRTMafterdenoising.大量的反射界面,逆时偏移之后会在整个剖面上存在低频噪音,如图4b所示,应用拉普拉斯算子去噪之后的结果如图4c,低频噪音得到有效压制,断层与界面成像清楚.2.4.2 Sigsbee2a模型图5a是Sigsbee2a偏移速度模型,与Marmousi模型不同,这个偏移速度模型对速度场做了平滑,只在盐丘边界处存在速度间断,相应的,低频噪音主要集中在盐丘边界强反射处,如图5b所示,这与上文中关于噪音产生机制的分析是一致的,图5c是应用拉普拉斯算子去噪之后的结果,低频噪音得到有效压制,盐丘边界与盐丘底部散射点成像清楚.3 存储从成像条件不难发现,逆时偏移成像需要炮点波场和检波点波场相同时刻的波场.然而与之矛盾的是炮点波场和检波点波场沿相反的时间方向递推,EricDussaud等[37]总结了解决这一矛盾的五种途径:(1)从狋=0开始计算每次成像时需要的与检波点相同时刻的炮点波场,这种方法不存在波场存储问题,但是计算量是犗(狀2),无法承受;(2)把炮点正传历史信息全部存下来,成像时读5712地球物理学报(ChineseJ.Geophys.)53卷取相应的信息,计算量是犗(狀),但存储整个传播历史是目前计算机硬件无法承受的;(3)每隔犽步存一个波场,中间结果利用插值得到,这种方法的计算量是犗(狀),存储变为原来的1/犽,这也是目前很多商业RTM软件使用的方法,缺点在于插值会损失精度,并且波场存储量依然巨大;(4)只在边界存储正传波场信息,反传时作为边界条件重构波场,这种方法增加了一倍的计算量,但存储量大幅度降低到计算机硬件可以承受的地步;(5)文献[34]中的设置“检查点”的最优方法.第一次正传时通过设置某些时刻作为“检查点”将波场存储下来,对于非检查点时刻的波场,通过从特定的检查点重新计算这一时刻的波场,选择合适的检查点,可以降低“重算比”.问题在于,检查点选得太少,重算比太高,检查点选得过多,所需要存储的临时波场越大,存储问题依然存在.本文采用文献[35]中的随机边界的思想,用计算换存储.并将这种方法借助GPU/CPU协同计算并行实现,避免了GPU和CPU之间的数据通讯,可以部分弥补增加的计算量.图5 Sigsbee2A模型去噪测试(a)偏移速度模型;(b)去噪之前的逆时偏移结果;(c)去噪之后的逆时偏移结果.Fig.5 ThedenoisetestonSigsbee2Amodel(a)Thevelocitymodel.(b)TheresultofRTMbeforedenoising.(c)TheresultofRTMafterdenoising.3.1 随机边界模型图6a是带随机边界的速度模型;将炮点置于模型的中央,图6(b,c,d)分别是狋=0.25s,0.75s和1.25s的波场快照;波场遇到随机边界产生漫反射,波前面在1.25s时已经变成随机噪音,这种随机噪音一方面在成像时不会形成连续的同相轴,另一方面波场反传时可以重构波场有用信息.图6e是波场从狋=1.25s反传回狋=0.25s时的波场快照;图6f是正反传波场狋=0.25s时刻的波场差异放大10000倍的结果,二者的差异可以忽略.图7a是初始时刻波场的空间波数谱;图7b是狋=1.25s时波场的空间波数谱;采用不同的随机边界参数,相同的震源子波信息做了20次试验,图7c是狋=1.25s时刻20次试验叠加结果的空间波数谱.结果表明,高波数噪音经过叠加后已经有效的去除,只有低波数(长波场)的信息依然存在,这与水波中的“衍射”现象类似.炮域叠前逆时偏移结果是将所有炮的偏移结果累加,这些高波数的噪音可以在累加中被消除.3.2 犌犘犝实现随机边界逆时偏移随机边界的引入用计算换存储,炮点正传波场与检波点波场可以同时沿时间反方向传播并成像,不需要存储炮点正传的历史信息.本文借助GPU/CPU协同计算实现随机边界叠前逆时偏移,图8a是存储传播历史的GPU/CPU协同计算逆时偏移流程,图8b是GPU/CPU协同计算随机边界实现逆时偏移的计算流程,从流程图上可以看出,炮点波场正传加反传虽然增加了一倍的计算量,但是可以有效地减少GPU和CPU之间的通讯,所以在计算时间上并不是增加两倍,实测的结果表明,二者的计算时间相差不多.以Sigsbee2a模型为例,流程一计算单炮的偏移时间为111s,而流程二的计算时间为136s,而不是166s,额外增加的计算量通过减少通讯部分得到弥补.6712 9期刘红伟等:地震叠前逆时偏移中的去噪与存储图6 随机边界模型测试(a)带随机边界的速度模型;(b)、(c)、(d)分别是狋=0.25s,0.75s和1.25s的波场快照;(e)波场反传回狋=0.25s时的波场快照;(f)正反传波场狋=0.25s时刻的波场差异放大10000倍的结果.Fig.6 Randomboundarymodeltest(a)Thevelocitymodelwithrandomboundaries;(b),(c)and(d)isthewavefieldsnapshotsat狋=0.25s,0.75sand1.25s;(e)Thesnapshotofbackwardpropagatingwavefieldat狋=0.25s;(f)The10000magnifieddifferencebetween(b)and(e).图7 空间波数谱(a)初始时刻波场的空间波数谱;(b)狋=1.25s时波场的空间波数谱;(c)狋=1.25s时刻20次试验叠加结果的空间波数谱.Fig.7 Wavenumberspectral(a)Thewavenumberspectralofthewavefieldattheoriginaltime;(b)Thewavenumberspectralofthewavefieldat狋=1.25s;(c)Thewavenumberspectralofthestackingwavefieldfrom20testsat狋=1.25s.3.3 模型测试图9a是二维盐丘速度模型,图9b是存储炮点传播历史的GPU/CPU协同计算逆时偏移去噪之后的结果;图9c是借助随机边界的GPU/CPU协同计算逆时偏移结果;图9d是采用与随机边界相同的计算流程,在边界处采用自由反射边界的GPU/CPU协同计算逆时偏移结果.图9(e,f,g)分别是图9(b,c,d)的局部放大结果.通过比较,图9(g)中的椭圆位置存在严重的噪音,是由于自由反射边界将连续的波场信息反射到成像区域造成的,而借助随机边界的GPU/CPU协同计算逆时偏移结果与存储传播历史的逆时偏移结果差别很小,证明了随机边界方法的有效性.7712地球物理学报(ChineseJ.Geophys.)53卷图8 逆时偏移流程(a)原始的GPU/CPU协同计算逆时偏移流程;(b)GPU/CPU协同计算随机边界实现逆时偏移的计算流程.Fig.8 TheflowchartsofRTM(a)TheoriginalflowchartofGPU/CPUcollaborativecomputingRTM;(b)TheflowchartofGPU/CPUcollaborativecomputingRTMusingrandomboundarycondition.4 结论叠前逆时偏移低频成像噪音主要是由于强反射边界的内反射造成的,本文提出的在叠前对数据进行相位和振幅校正,在成像后运用拉普拉斯算子滤波法消除成像噪音的方法可以完全消除入射角为90°的低频噪音,同时指出入射角小于90°的噪音不能被完全消除的合理性.借助GPU/CPU协同计算并行实现的随机边界逆时偏移方法,用计算换存储,并且节省了GPU与CPU之间的数据通讯,实验结果表明,采用随机边界方法的偏移结果与直接存储波场的方法得到的结果差别很小,可以用来解决叠前逆时偏移存储量大的问题.图9 二维盐丘模型测试(a)速度模型;(b)原始的GPU/CPU协同计算逆时偏移结果;(c)借助随机边界的GPU/CPU协同计算逆时偏移结果;(d)自由反射边界GPU/CPU协同计算逆时偏移结果;(e)、(f)、(g)分别是(b)、(c)、(d)的局部放大结果.Fig.9 Teston2Dsaltmodel(a)Thevelocitymodel;(b)TheresultofthetraditionalGPU/CPUcollaborativecomputingRTM;(c)TheresultofGPU/CPUcollaborativecomputingRTMusingrandomboundarycondition;(d)TheresultofGPU/CPUcollaborativecomputingRTMusingfreereflectionboundarycondition;(e),(f)and(g)isthepartialzoomoutof(b),(c)and(d),respectively.下一步的研究重点是开展GPU集群实现超大规模逆时偏移研究.致谢对SmartJV提供的Sigsbee模型数据,谨致谢意,感谢李幼铭研究员对本文的指导和帮助.感谢吉星达科技有限公司刘钦和佟小龙对本文的指导和帮助,感谢吉星达科技有限公司允许本文发表.参考文献(References)[1]HemonC.Equationsd′ondeetmodeles.犌犲狅狆犺狔狊犻犮犪犾犘狉狅狊狆犲犮狋犻狀犵,1978,26:790~821[2]EdipBaysal,DanDKosloff,JohnWCSherwood.Reversetimemigration.犌犲狅狆犺狔狊犻犮狊,1983,48(11):1514~1524[3]McMechanGA.Migrationbyextrapolationoftimedependent8712 9期刘红伟等:地震叠前逆时偏移中的去噪与存储boundaryvalues.犌犲狅狆犺狔狊犻犮犪犾犘狉狅狊狆犲犮狋犻狀犵,1983,31:413~420[4]DanLoewenthal,IrshadRMufti.Reversedtimemigrationinspatialfrequencydomain.犌犲狅狆犺狔狊犻犮狊,1983,48(5):627~635[5]WhitmoreND.Iterativedepthmigrationbybackwardtimepropagation.53rdAnnualInternationalMeeting,SEG,ExpandedAbstracts,1983.827~830[6]McMechanGA,ChangWF.3Dacousticprestackreversetimemigration.犌犲狅狆犺狔狊犻犮犪犾犘狉狅狊狆犲犮狋犻狀犵,1990,38(7):737~756[7]ChangWenFong,GeorgeA.McMechan.3Delasticprestack,reversetimedepthmigration.犌犲狅狆犺狔狊犻犮狊,1994,59:597~609[8]ZhangHouzhu(James),ZhangYu.Reversetimemigrationin3DheterogeneousTTImedia.78thAnnualInternationalMeeting,SEGExpandedAbstracts,2008,27.2196~2200[9]RobinPFletcher,DuXiang,PaulJ.FowlerReversetimemigrationintiltedtransverselyisotropic(TTI)media.犌犲狅狆犺狔狊犻犮狊,2009,74:WCA179[10]HuangTony,ZhangYu,ZhangHouzhu,etal.SubsaltimagingusingTTIreversetimemigration.犜犺犲犔犲犪犱犻狀犵犈犱犵犲,2009,28:448~452[11]ZhangYu,ZhangHouzhu.AstableTTIreversetimemigrationanditsimplementation.79狋犺犃狀狀狌犪犾犐狀狋犲狉狀犪狋犻狅狀犪犾犕犲犲狋犻狀犵,犛犈犌犈狓狆犪狀犱犲犱犃犫狊狋狉犪犮狋狊,2009,28:2794~2798[12]邓玉琼,戴霆范,郭宗汾等.弹性波叠前有限元反时偏移.石油物探,1990,29(3):22~33TengYC,DaiTF,JohnT.Kuo.Prestackfiniteelementreversetimemigrationforelasticwaves.犌犲狅狆犺狔狊犻犮犪犾犘狉狅狊狆犲犮狋犻狀犵犳狅狉犘犲狋狉狅犾犲狌犿(inChinese),1990,29(3):22~33[13]张美根,王妙月.各向异性弹性波有限元叠前逆时偏移.地球物理学报,2001,44(5):711~719ZhangMG,WangMY.Prestackfiniteelementreversetimemigrationforanisotropicelasticwaves.犆犺犻狀犲狊犲犑.犌犲狅狆犺狔狊.(inChinese),2001,44(5):711~719[14]底青云,王妙月.弹性波有限元逆时偏移技术研究.地球物理学报,1997,40(4):570~579DiQY,WangMY.Thestudyoffiniteelementreversetimemigrationforelasticwave.犆犺犻狀犲狊犲犑.犌犲狅狆犺狔狊.(犃犮狋犪犌犲狅狆犺狔狊犻犮犪犛犻狀犻犮犪)(inChinese),1997,40(4):570~579[15]王童奎,付兴深,朱德献等.谱元法叠前逆时偏移研究.地球物理学进展,2008,23(3):681~685WangTK,FuXS,ZhuDX,etal.Spectralelementmethodforprestackreversetimemigration.犘狉狅犵狉犲狊狊犻狀犌犲狅狆犺狔狊犻犮狊(inChinese),2008,23(3):681~685[16]陈可洋.基于高阶有限差分的波动方程叠前逆时偏移方法.石油物探,2009,48(5):475~478ChenKY.Highorderfinitedifferenceprestackreversetimemigration.犌犲狅狆犺狔狊犻犮犪犾犘狉狅狊狆犲犮狋犻狀犵犳狅狉犘犲狋狉狅犾犲狌犿(inChinese),2009,48(5):475~478[17]王山山,聂勋碧,邓玉琼.混合法双程无反射波动方程偏移.石油地球物理勘探,1993,28(5):537~542WangSS,NieXB,DengYQ.Hybridmigrationusingtwowaynonreflectingwaveequation.犗犻犾犌犲狅狆犺狔狊犻犮犪犾犘狉狅狊狆犲犮狋犻狀犵(inChinese),1993,28(5):537~542[18]薛东川,王尚旭.波动方程有限元叠前逆时偏移.石油地球物理勘探,2008,43(1):17~21XueDC,WangSX.Waveequationfiniteelementprestackreversetimemigration.犗犻犾犌犲狅狆犺狔狊犻犮犪犾犘狉狅狊狆犲犮狋犻狀犵(inChinese),2008,43(1):17~21[19]杜启振,秦童.横向各向同性介质弹性波多分量叠前逆时偏移.地球物理学报,2009,52(3):801~807DuQZ,QinT.Multicomponentprestackreversetimemigrationofelasticwavesintransverseisotropicmedium.犆犺犻狀犲狊犲犑.犌犲狅狆犺狔狊.(inChinese),2009,52(3):801~807[20]刘红伟,李博,刘洪等.地震叠前逆时偏移高阶有限差分算法及GPU实现.地球物理学报,2010,53(7):1725~1733LiuHW,LiB,LiuH,etal.ThealgorithmofhighorderfinitedifferenceprestackreversetimemigrationandGPUimplementation.犆犺犻狀犲狊犲犑.犌犲狅狆犺狔狊.(inChinese),2010,53(7):1725~1733[21]BaysalE,KosloffDD,SherwoodJWC.Atwowaynonreflectingwaveequation.犌犲狅狆犺狔狊犻犮狊,1984,49:132~141[22]LoewenthalD,StoffaPA,Fa
展开阅读全文
  石油文库所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
0条评论

还可以输入200字符

暂无评论,赶快抢占沙发吧。

关于本文
本文标题:地震叠前逆时偏移中的去噪与存储
链接地址:http://www.oilwenku.com/p-11018.html
关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服客服 - 联系我们
copyright@ 2016-2020 石油文库网站版权所有
经营许可证编号:川B2-20120048,ICP备案号:蜀ICP备11026253号-10号
收起
展开