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

反褶积干涉成像及其在随钻地震数据处理中的应用

关 键 词:
反褶积 干涉 成像 及其 地震 数据处理 中的 应用
资源描述:
书书书第25卷第3期2010年6月(页码:951~956)地球物理学进展PROGRESS IN GEOPHYSICSVol.25,No.3June2010黄伟传,葛洪魁,王宝善,等.反褶积干涉成像及其在随钻地震数据处理中的应用.地球物理学进展,2010,25(3):951~956,DOI:10.3969/j.issn.10042903.2010.03.032.HuangWC,GeHK,WangBS,犲狋犪犾.Deconvolutioninterferometryanditsapplicationtoseismicwhiledrillingdataprocessing.犘狉狅犵狉犲狊狊犻狀犌犲狅狆犺狔狊.(inChinese),2010,25(3):951~956,DOI:10.3969/j.issn.10042903.2010.03.032.反褶积干涉成像及其在随钻地震数据处理中的应用黄伟传1,葛洪魁1,王宝善1,王兆明2(1.中国地震局地球物理研究所,北京100081;2.中国石油勘探开发研究院,北京100083)摘要随钻地震是石油工程领域的一项新技术,是利用钻头破岩的连续随机振动信号进行钻前预测,通常需要利用钻头震源函数与地面记录互相关进行处理,而地下数千米震源函数的获取难度很大,成为制约随钻地震技术的一个关键问题.反褶积干涉成像是在互相关基础上发展起来的,可以利用随机噪音获得地震台站间的格林函数,对地下构造和速度结构成像,同时这种方法不依赖于震源函数信息.本文回顾了干涉成像的发展历程,分析了互相关干涉和反褶积干涉的原理、方法和特点,结合随钻地震数据处理需要,结合实例分析了反褶积干涉应用于随钻地震数据处理的前景,为发展新型随钻地震技术提供了重要参考信息.关键词反褶积,干涉成像,互相关,随钻地震,数据处理DOI:10.3969/j.issn.10042903.2010.03.032    中图分类号P631     文献标识码A     犇犲犮狅狀狏狅犾狌狋犻狅狀犻狀狋犲狉犳犲狉狅犿犲狋狉狔犪狀犱犻狋狊犪狆狆犾犻犮犪狋犻狅狀狋狅狊犲犻狊犿犻犮狑犺犻犾犲犱狉犻犾犾犻狀犵犱犪狋犪狆狉狅犮犲狊狊犻狀犵HUANGWeichuan1,GEHongkui1,WANGBaoshan1,WANGZhaoming2(1.犐狀狊狋犻狋狌狋犲狅犳犌犲狅狆犺狔狊犻犮狊,犆犺犻狀犪犈犪狉狋犺狇狌犪犽犲犃犱犿犻狀犻狊狋狉犪狋犻狅狀,犅犲犻犼犻狀犵100081,犆犺犻狀犪;2.犚犲狊犲犪狉犮犺犐狀狊狋犻狋狌狋犲狅犳犘犲狋狉狅犾犲狌犿犈狓狆犾狅狉犪狋犻狅狀牔犇犲狏犲犾狅狆犿犲狀狋,犘犲狋狉狅犆犺犻狀犪,犅犲犻犼犻狀犵100083,犆犺犻狀犪)犃犫狊狋狉犪犮狋 SeismicWhileDrilling(SWD)isanewtechnologyappliedinthefieldofpetroleum.Itusesthecontinuousrandomvibrationtomakebitforecastbycrosscorrelationbetweenthesourcefunctionofthousandsmetersundergroundandthewaverecordedonthesurface.ThekeyissueoftheSWDishowtoobtaintheseismicsourcefunctionatsuchadepthwhichblocksthewideapplicationofSWD.Thedeconvolutioninterferometry,developedbasedoncrosscorrelation,isanothertechnologythatcanbeusedinSWD.Itdoesn’tneedthesourcefunctionandcanderivetheGreen’sfunctionofstationpairsbytheambientnoise,andimagethestructureandvelocitypropertyofthesubsurface.Itcaneliminatetheeffectsofthesourceandgetgoodimageresult.Thispaperreviewsthedevelopmentsofthedeconvolutioninterferometry,andanalyzesthetheoryandcharactersofthecrosscorrelationanddeconvolutioninterferometry.WithanexampleofdeconvolutioninterferometryapplicationtotheSWDdataprocessing,itdemonstratesthegreatpotentialofdeconvolutioninterferometryinSWD,andprovidesimportantinformationfornewtechnologydevelopmentsofSWD.犓犲狔狑狅狉犱狊 deconvolution,interferometry,crosscorrelation,seismicwhiledrilling,dataprocessing收稿日期20090920;修回日期20100115.基金项目国家“863”计划项目“随钻地震技术研究”(2006AA06A108)和中央级公益性科研院所基本科研业务专项项目(DQJB08B15)联合资助.作者简介黄伟传,1969年生,博士,从事地震资料的处理和解释研究工作.(Email:sucun2003@163.com)地球物理学进展25卷0 引言随钻地震(SWDSeismicWhileDrilling)是石油工程领域中的一项先进技术,它利用钻头旋转破岩振动作震源,进行地层探测,从而实现钻前预测,对于安全高效钻井施工具有重要意义.由于钻头振动为连续随机信号,处理方法为利用由钻柱顶端传感器记录到的参考信号提取钻头震源函数,与地面测线上记录到的信号作互相关处理,获取从地层中传播的直达波和反射波信息(图1).图1 随钻地震互相关处理原理图Fig.1 Diagramshowingtheprincipalofthecrosscorrelaionintheseismicwhiledrilling根据互相关原理可知,源函数估计精度对互相关处理具有决定性影响.而钻头信号在沿钻柱的传播过程中受到钻柱滤波效应的强烈改造,加上来自地面和井壁的非稳态噪音,使得利用钻柱顶端振动信号恢复钻头源函数的难度很大,往往失真严重,成为发展SWD技术的制约因素.如能避开钻头源函数估计,进行钻前速度结构成像,是解决了SWD技术发展的一个关键问题.同时,不进行钻柱顶端振动测量,可大大简化现场施工,也有可能使许多没有进行参考信号测量的地面记录资料得以利用,有利于SWD技术的推广,对于发展SWD技术具有非常重要的意义.连续随机信号成像在地震观测数据分析中得到较好的发展,以噪音记录互相关为基础的干涉成像取得了重要的成果.近几年在互相关干涉处理基础上发展起了一种新型成像技术—反褶积干涉成像.反褶积干涉成像的一个显著优点是不需要震源函数信息,特别适合于随钻地震数据处理.本文系统介绍了干涉成像原理与方法,对比分析了反褶积干涉和互相关干涉特点,结合实例探讨了将反褶积干涉成像方法用于随钻地震资料处理的优点和发展前景,为发展新型随钻地震技术提供重要的参考信息.1 干涉成像技术的发展在1976年,Claerbout就通过分析一维模型的介质反射和透射特征,发现通过互相关处理可以得到地层反射响应特征[1],Wapenaar又根据格林函数的互易性理论,进行了严格的证明,地震干涉成像在任意介质中都是成立的.干涉测量一般是指通过研究信号对之间的干涉现象,从它们之间的差异中获得信息的方法.1987年Cole为了证明这个猜想,在斯坦福校园大约0.5km2的范围内,利用4096个检波器记录地震数据,尽管实验工作没有取得预期效果,但是地震学家发现了从远距离传来的高频地震波[2].2000年的斯坦福大学勘探项目中,Gerard利用互相关处理方法对环境噪声和人工震源数据都进行了处理,并引入了干涉成像概念,主要包括了互相关和偏移[3].现已研究表明,噪声中包含有大量有用的信息,通过对地面上两个地震台站记录的背景地震噪声信号在长时间段内进行互相关运算,可以提取格林函数,从而对地下深部进行地震成像研究.利用环境噪声干涉成像,进行地下构造成像是地球物理的新方法,国际上有多个研究小组从事干涉成像研究.科罗拉多矿业学院波现象研究中心Snieder和物理系地球内部成像中心Shapiro等人,在2004年[2]利用地震资料互相关处理,提取到台站间的格林函数,通过对几天的环境噪声进行互相关计算,得到宽频带的瑞雷波散射特征和群速度.在2005年,对美国加州USArray台阵中62个地震台站记录到的一个月的地震背景噪声进行互相关计算,得到的波形与台站间的格林函数仅有幅度的差异,并由此得到了短周期(7~18s)的面波频散图.法国地球内部物理实验室的Campillo等利用散射地震的Coda波互相关,得到低频的相干散射波场,并能够得到波的极化特性和群速度.Utah大学层析成像、模拟与偏移研究组的Schuster等,利用干涉理论对随钻地震数据进行了处理,它主要包括地震干涉和偏移,并对实际随钻数据进行了处理.2007259 3期黄伟传,等:反褶积干涉成像及其在随钻地震数据处理中的应用年Florent等[4]对火山噪声进行检测,并利用火山噪声重建了地震台站对间的瑞雷面波格林函数,并提取了剪切波速,建立了三维速度模型,与其它成像结果相吻合.Pierre等[5]在地震勘探中使用互相关处理技术,研究散射的影响和震源分布对格林函数收敛的影响.试验是在1平方公里的面积上布置了1600个检波器,网格大小是25×25m,采用主动震源,炮在两个方向上的移动距离都是12.5m.当利用直达波时,震源的位置控制着互相关提取格林函数的收敛,对于尾波来说,震源分布影响很小,同样也可提取格林函数;在各向均匀波场中,重建格林函数时,不仅可以提取不同偏移距的相位,也可以提取相对振幅.对地震干涉研究,在国内许多研究人员也做了大量的研究,取得非常好的成果.房立华等[6]利用华北地震科学台阵2007年1~4月的地震噪声数据,通过互相关提取了瑞利波的经验格林函数,并处理得到了瑞利波面波的群速度频散曲线,计算了多层的瑞利面波速度分布图,揭示了地壳内部横向变化.在2007年王宝善等[8]通过在云南一个多月的试验,利用尾波(coda)干涉方法对地下界面的速度变化进行检测,发现了速度变化的趋势和每天的周期变化.2009年王宝善等的个人通讯中,利用四川台网地震前后各50天的连续波形记录,得到很好的互相关格林函数,计算了面波的走时变化,利用面波成像的方法得到因汶川地震导致该地区的地震波速变化.齐诚等[10]认为通过对地面上两个地震台站记录的背景地震噪声信号在长时间段内进行互相关运算,提取格林函数,从而进行地震成像,获取对地下结构的认识,将对地震学的发展会有不可估量的影响.在地震干涉研究过程中,发现互相关干涉受到震源特性的重要影响,Snieder等在互相关干涉成像的基础上,分析了地震响应的褶积特征,提出了反褶积干涉成像理论,通过反褶积可消除震源函数的影响.Snieder等[11]在2006年利用不同楼层高度地震记录来计算大楼的脉冲响应,通过反褶积处理来恢复建筑物的弹性响应.Curtis等人[12]研究表明,在任意多维介质中可以通过反褶积来实现干涉成像,而且成像结果不受震源特征影响.在2008年,Wapenaar等[13]对噪音进行了多维反褶积数值模拟试验,表明多维反褶积干涉成像具有震源畸变校正的潜在能力.Ivan和Snieder在圣安蒂斯断层监测中,将反褶积干涉成像方法应用于随钻地震资料处理,在没有参考信号的条件下得到了很好成像效果[14~15].2 干涉成像理论地震干涉成像的方法主要是互相关处理,通过互相关处理,利用地面的两个地震记录,可以得到一个台站作为震源,另一台站接收时的格林函数.在互相关处理基础上,Snieder等又发展了反褶积干涉处理技术.三维的格林函数在接收点之间的关系是互相关的负时间导数[16~17],公式为:-dd狋犆disp犻犼(狉1,狉2,狋[])≈犌disp犻犼(狉1,狉2,狋)-犌disp犻犼(狉1,狉2,-狋).(1)这里的犌disp犻犼是单位力在犻方向狉1接收点对犼方向狉2点的位移响应,[]表示一段很长时间的平均,在处理中通常需要特别长时间的记录.在较窄的频带时,对时间的导数,就近似为他们的互相关函数.犮犻犼()狋=犌犻犼()狋. 狋>0(2)利用互相关处理对地下构造成像,结果稳定,信噪比高,获得了很好的效果.目前,在瑞利波和勒夫波成像的基础上,又展开了对小距离台站之间的数据处理,得到体波的反射特征,并利用体波进行高分辨率成像.反褶积干涉是依据格林函数的互易理论[14],利用任意两个接收点的地面地震资料进行反褶积处理,可以得到他们之间的格林函数.对于同一域内的两个点接收点犪、犫,在频率域内犪点的褶积公式为狌(狉犪,狊,ω)=犠(狊,ω)犌(狉犪,狊,ω),其中:犌为格林函数,犠为频率域源函数.犫点的褶积公式为狌(狉犫,狊,ω)=犠(狊,ω)犌(狉犫,狊,ω).a、b两点地震信息做互相关处理,那么它们的频率域公式为:犆犪犫=犠(狊)2犌(狉犪,狊)犌(狉犫,狊),(3)其中犠(狊)2为震源影响因子,是非常复杂的频率函数,随震源s的变化而变化,这里的星号表示是复数的共轭.从方程得到,互相关犆犪犫受到犠(狊)能量谱的影响.对接收点A、B点的记录信息狌(狉犫,狊)和狌(狉犪,狊)进行反褶积处理的频率域公式为:犇犪犫=狌(狉犪,狊)狌(狉犫,狊)=狌(狉犪,狊)狌(狉犫,狊)狘狌(狉犫,狊)狘2=犌(狉犪,狊)犌(狉犫,狊)狘犌(狉犫,狊)狘2.(4)通过反褶积处理去掉了震源函数特征对处理结359地球物理学进展25卷果的影响,在不知震源特征的条件下,可以对信号进行波场分离和速度成像.对于弹性介质中的反褶积干涉方法,就像在声波介质中一样,可以定义为下列的运算:犇(狆,狇)犃犅=犇(狆,狇)犃犅,犓=狌(狏,φ)(狆,犓)(狉犃,狊,ω)狌(狏,φ)(狇,犓)(狉犅,狊,ω)=犌(狏,φ)(狆,犓)(狉犃,狊,ω)犌(狏,φ)(狆,犓)(狉犃,狊,ω).(5)这里狌(狏,φ)(狆,犓)=犠犓犌(狏,φ)(狆,犓)是测量的响应在频率域,下标狆和狇表示一个特定的记录粒子速度的一个分量(上标狏表示场的离子速度大小).3 反褶积干涉成像在随钻地震中的应用地震干涉在随钻地震处理中得到了一定的应用,并发展虚拟炮集和自相关偏移新的成像方法,反褶积干涉方法也在实际数据处理中得到应用.Snieder等人对反褶积干涉方法进行了详细的研究,并在圣安蒂斯断层监测中,利用反褶积干涉成像对随钻地震资料进行了成像.圣安蒂斯断层监测作为美国透镜计划(EarthScope)计划的一部分,美国地震科学家在Parkfield断层带进行了科学钻探工作,位置在靠近1966年发生的帕克菲尔德地震震中,钻一口4000m的斜井,穿过断层带达到另一侧未受扰动的围岩.在钻井中进行取心、井中测量,了解断层带的成分、物理状态和力学性质,建立圣安德列斯断层深部观测站(SAFOD),进行长期连续的直接测量,开展地震物理机制研究.在钻井过程中同时进行了随钻监测,钻探井有主井和参考井两个井筒,检波器放在参考井中,共放置了32个检波器.主井由钻机钻进,作为震源.由于在这次的实验中没有记录参考信号,这样依赖于震源函数的互相关干涉成像就遇到了困难,反褶积干涉技术为地下界面成像提供了技术方法.Ivan和Snieder[15]将反褶积干涉用于SAFOD随钻地震成像,在处理中对所有信号进行了低通滤波,保留了最高频率到55hz,并进行了垂直叠加和反褶积干涉成像处理,对不同分量的地震资料进行了处理,得到了高分辨率成像成果.由于钻头激发的地震波是连续的,震源脉冲对于时间来说是均匀的.图2 是通过钻头地震噪声成像得到的结果(灰色)叠放在Sigbee模型上(a)利用虚炮集记录通过波动方程偏移得到的;(b)通过互相关处理得到的虚炮集记录(引自Ivan[15])Fig.2 Imagesobtanedfromdrillbitnoiseinterferometry.Images(ingray)aresuperposedontheSigbeemodel(a)Theimageobtainedfromshot-profilewave-equationmigrationofpseudoshotgathersgeneratedfromdeconvoluionintergerometry;(b)Theresultofmigratingcorrelation-basedinterferometricshotgathers.Redlinesintheimagesrepresentthereceiverarray(from[15])反褶积干涉成像能够更好地压缩震源函数,提高干涉成像对地下构造成像能力,同时还可以分辨出盐丘沉积反射体,具有更高的分辨率.图2是地震干涉互相关和反褶积处理的对比图,图2a是反褶积干涉处理的结果,可以看到反褶积的处理结果分辨率高,成像精度高;图2b互相关处理的结果受到的干扰影响严重,而且分辨率低.这主要是由于反褶积干涉处理能够很好地压制钻头震源函数,消除震源的影响,得到虚炮集波场特征明显,成像结果就不会出现严重的畸变.4 反褶积干涉成像在SWD数据处理中的应用展望地震干涉是最新的地球物理技术,具有非常广的应用前景.利用地震干涉可以进行地壳速度结构成像,利用无时无刻不存在的自然噪音对地震危险区进行动态监测,分析监测目的地区的应力状态、速度变化,有助于获取深部动态信息,开展地震预报新459 3期黄伟传,等:反褶积干涉成像及其在随钻地震数据处理中的应用途径的探索.通过地震干涉技术,可以通过实时监控矿井中的微震,对微震事件进行及时分析研究,预测矿井中的瓦斯突出、塌方等灾害,做到对矿井灾难的监测和提前预报.反褶积干涉成像新方法可以有效地消除震源的影响,提高干涉处理的分辨率,可以在强噪声背景下提取有效信号.在随钻地震处理中,将反褶积干涉成像技术引入SWD数据处理中,在一定程度上能解决随钻地震数据处理中钻井环境和震源性质对处理结果的严重影响,提高成像结果稳定性和分辨率.通过开发研究高分辨率的干涉成像技术,可以突破石油勘探原来地震成像空白区,如火成岩、碳酸盐岩的屏蔽区,对岩下构造进行研究,可以带来勘探上新发现和突破,对隐蔽油气藏、特殊油气藏进行定量的描述.在油田开发中可以实现对开发井实时成像、实时监督,根据实时资料对开发井做出合理调整.新技术的应用可改变现有的勘探开发模式,给石油工业带来强劲发展动力.参考文献(References):[1]JonFC.Synthesisofalayeredmediumfromitsacoustictransmissionresponse[J].Geophysics.1968,33(2):264~269.[2]ShapiroNM,CampilloM.EmergenceofbroadbandRayleighwavesfromcorrelationsoftheambientseismicnoise[J].GeophysicalResearchLetters.2004,31,L07614,doi:10.1029/2004GL019491.[3]KeesW,DeyanD,JohanR.Introductiontothesupplementonseismicinterferometry[J].Geophysics,2006,71(4):1~4.[4]FlorentB,NikolaiMS,MichelC,SlexandreN,ValérieF.3DsurfacewavetomographyofthePitondelaFournaisevolcanousingseismicnoisecorrelations[J].GeophysicalResearchLetters,2007,34,L02305,doi:10.1029/2006GL028585.[5]PierreG,PhilippeR,MichelC,犲狋犪犾.ConvergenceofthetwopointcorrelationfunctiontowardtheGreen’sfunctioninthecontextofaseismicprospectingdataset[J].Geophysics,2008,73(6):47~53.[6]房立华,吴建平,吕作勇.华北地区基于噪声的瑞利面波群速度层析成像[J].地球物理学报,2009,52(3):663~671.FangLH,WuJP,LuZY.RayleighwavegroupvelocitytomographyfromambientseismicnoiseinNorthChina[J].ChineseJ.Geophys.(inChinese),2009,52(3):663~671.[7]齐诚,赵大鹏,陈,等.首都圈地区地壳P波和S波三维速度结构及其与大地震的关系[J].地球物理学报,2006,49(3):805~815.QiC,ZhaoDP,ChenY,犲狋犪犾.3DPandSwavevelocitystructuresandtheirrelationshiptostrongearthquakesintheChinesecapitalregion[J].ChineseJ.Geophys.(inChinese),2006,49(3):805~815.[8]WangBS,ZhuP,ChenY,犲狋犪犾.Continuoussubsurfacevelocitymeasurementwithcodawaveinterferometry[J].JournalofGeophysicalResearch2007,113,B12312,doi:10.1029/2007JB005023.[9]刘启元,李昱,陈九辉,等.汶川Ms8.0地震:地壳上地幔S波速度结构的初步研究[J].地球物理学报,2009,52(2):309~319.LiuQY,LiY,ChenJH,犲狋犪犾.WenchuanMs8.0earthquake:preliminarystudyoftheSwavevelocitystructureofthecrustanduppermantle[J].ChineseJ.Geophys.(inChinese),2009,52(2):309~319.[10]齐诚,陈棋福,陈.利用背景噪声进行地震成像的新方法[J].地球物理学进展,2007,22(3):771~777.QiC,ChenQF,ChenY.Anewmethodforseismicimagingfromambientseismicnoise[J].ProgressInGephysics,2007,22(3):771~777.[11]ReolS,ErdalS.ExtractingtheBuildingResponseUsingSeismicInterferometry:TheoryandApplicationtotheMillikanLibraryinPasadena[J].California.BulletinoftheSeismologicalSocietyofAmerica,2005,96(2):586~598.[12]AndrewC,PeterG,HaruoS,犲狋犪犾.Seismicinterferometryturningnoiseintosignal[J].TheLeadingEdge,2006,25(9):1082~1092.[13]KeesW,JoostvanderNeut,ElmerR.Passiveseismicinterferometrybymultidimensionaldeconvolution[J].Geophysics,2008,73(6):51~56.[14]IvanVR.Interferometrybydecomvolution,Part1Theoryforacousticwavesandnumericalexamples[J].Geophysics,2008,73(3):115~128.[15]IvanVR.Interferometrybydecomvolution,Part2Theoryforelasticwavesandapplicationtodrillbitseismicimaging[J].Geophysics,2008,73(3):129~141.[16]KeesW.Green’sfunctionretrievalbycrosscorrelationincaseofonesidedillumination[J].GeophysicalResearchLetters,2006,33,L19304,doi:10.1029/2006GL027747.[17]KarimGS,PeterG,PhilippeR,犲狋犪犾.ExtractingtimedomainGreen’sfunctionestimatesfromambientseismicnoise[J].Geophysicalresearchletters,2005,32,L03310,doi:10.1029/2004GL021862.[18]SergeyF,LongJ.Timelapseimageregistrationusingthelocalsimilarityattribute[J].Geophysics,2009,74(2):7~11.[19]MichelC,AnneP.LongRangeCorrelationsintheDiffuseSeismicCoda.Science,2003,299,5606:547~549.[20]AndrewK,XunL.Methodologyfortomographicimagingaheadofminingusingtheshearerasaseismicsource[J].Geophysics,2009,74(2),1~8.[21]LangqiuFS,BerndM,DouglasRS.Measuringvelocitydispersionandattenuationintheexplorationseismicfrequencyband[J].Geophysics,2009,72(2):113~122.[22]NikolaiMS,MichelC,LaurentS,MichaelHR.HighResolutionSurfacewavetomographyfromAmbientSeismicNoise[J].Science,2005,307(5715):1615~1618.559地球物理学进展25卷[23]PhilippeR,KarimGS,PterG,犲狋犪犾.Pwavesfromcrosscorrelationofseismicnoise[J].GeophysicalResearchLetters,2005,32,L19303,doi:10.1029/2005GL023803.[24]吴庆举,李永华,张瑞青,等.用多道反褶积方法测定台站接收函数[J].地球物理学报,2007,50(3):791~796.WuQJ,LiYH,ZhangRQ,犲狋犪犾.Receiverfunctionestimatedbymultichanneldeconvolution[J].ChineseJ.Geophys.(inChinese),2007,50(3):791~796.[25]NicholasH,PeterG,CatherineAR,犲狋犪犾.PhaseVelocitiesfromseismicnoiseusingbeamformingandcrosscorrelationinCostaRicaandNicaragua[J].GeophysicalResearchLetters,2008,35,L19303,doi:10.1029/2008GL035387.[26]郝春月,边银菊.地震台阵、台阵地震学在中国的现状与发展[J].国际地震动态,2007,S345(9):36~42.HaoCY,BianYJ.TheStatusandDevelopmentofSeismicArrayandArraySeismologyinChina[J].RecentDevelopmentsinWorldSeismology.2007,S345(9):36~42.[27]李平,许厚泽,朱介寿,等.地震面波频散的层析成像[J].东北地震研究,1998,14(3):39~49.LiP,XuHZ,ZhuJS,犲狋犪犾.TomographyofRayleighWaveDispersion[J].SeismologicalResearchofNortheastChina.1998,14(3):39~49.[28]陈,朱日祥.设立“地下明灯研究计划”的建议[J].地球物理学进展,2005,20(5):485~489.ChenY,ZhuRX.ProposedProjectof“UmdergroundBrightLump”[J].AdvancesinEarthScience,2005,20(5):485~489.[29]李世愚,和雪松,张天中,等.地震学在减轻矿山地质灾害中的应用进展[J].国际地震动态,2006,328:1~9.LiSY,HeXS,ZhangTZ,犲狋犪犾.RecentDevelopmentsofMiningSeismologyforMiningHazards[J].RecentDevelopmentsinWorldSeismology,2006,328:1~9.[30]张碧星,鲁来玉.用频率波数法分析瑞利波频散曲线.工程地球物理学报,2005,2(4):245~255.ZhangBX,LuLY.InvestigationontheDispersionCurvesofRayleighWaveByFrequencyWavenumberAnalysisMethod[J].ChineseJournalofEngineeringGeophysic,2005,2(4):245~255.[31]SteckLK,PhillipsWS,MackeyK,犲狋犪犾.SeismictomographyofcrustalPandSacrossEurasia[J].GeophysJInt,2009,177:17781~92.[32]LauraAB,PeterG,DavidPK.Multichannelarraydiagnosisusingnoisecrosscorrelation[J].J.Acoust.Soc.Am,2008,124(4):203~209.[33]SebastianR,ChristineT.ArraySeismlolgy:MethodsandApplications[J].Geophysics,2002,40(3):1~26.[34]VanWK,LevshinAL.Surfacewavedispersionfromsmallverticalscatterers[J].GeoPhysicalResearchLetters,2004,31,L20602,doi:10.1029/2004GL021007.[35]ShuoM,GermanAP,GregoryCB.Testingcommunityvelocitymodelsforsoutherncaliforniausingtheambientseismicfield[J].BulletinoftheSeismologicalSocietyofAmerica,2008,98(6):2694~2714.[36]KuleshM,DialloMS,HolschneiderM,atel.Polarizationanalysisinthewaveletdomainbasedontheadaptivecovariancemethod[J].GeophysJInt,2007,68,doi:10.1111/j.1365246X.2007.03417.x.[37]BensenGD,RitzwollerMH,BarminMP,atel.Processingseismicambientnoisedatatoobtainreliablebroadbandsurfacewavedispersionmeasurements[J].GeophysJInt,2007,169:1239~1260.[38]RoelS,SlexandreG,HuubD.CodaWaveInterferometryforestimatingnonlinearbehaviorinseismicvelocity[J].Science,2002,295(22):2253~2255.[39]MehtaK,SniederR,GraizerV.Extractionofnearsurfacepropertiesforalossylayeredmediumusingthepropagatormatrix[J].Geophys.J.Int,2007,169,doi:10.1111/j.1365246X.2006.03303.x.659
展开阅读全文
  石油文库所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
0条评论

还可以输入200字符

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

关于本文
本文标题:反褶积干涉成像及其在随钻地震数据处理中的应用
链接地址:http://www.oilwenku.com/p-18893.html
关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服客服 - 联系我们
copyright@ 2016-2020 石油文库网站版权所有
经营许可证编号:川B2-20120048,ICP备案号:蜀ICP备11026253号-10号
收起
展开