井间地震纵横波分离方法研究及应用

StudyonPandSwaveseparationmethodincross-wellseismicanditsapplicationFENGYuping1，DIAORui2，XUWeixiu1，YANGJing1，SONGJianguo3（1.ShengliBranchCompany，GeophysicalCorporation，SinopecOilfieldServiceCorporation，Dongying257086，China；2.GeophysicsResearchInstituteofShengliOilfieldBranchCompany，Sinopec，Dongying257000，China；3.SchoolofGeophysics，ChinaUniversityofPetroleum，Qingdao266555，China）Abstract：Three-componentseismicacquisitiontechnologyusedinVSPandcross-wellseismiccanrecordwealthyseismicwavefield.ThePandSwavecoupletogether,whichbecomesthedifficultyandemphasisofdataprocessing.Throughanalyzingboreholeseismicacquisitionsystem,combinedwithpolarityrelationofPandSwaveindifferentcomponents,wecanassociatepolarityofPandSwavewithapparentvelocityanddirectionofseismicwavereachingthegeophonearray.Accordingtothisrelation,weobtainedthedistributionrelationofPandSwaveinτ-pdomainbasedonthelinearRadontransform（LRT）.ConsideringtheconvergenceofLRT,wedesignedapolarizationfilteringmethodusingthehigh-resolutionRadontransform,andthenaccomplishedtheseparationofPandSwave.Theprocessingofsimulateddataoftheoreticalmodelandpracticalcross-wellmodelseparatedPandSwavefieldsuccessfully,whichisusefulforindependentimagingofPandSwave.Keywords：cross-wellseismic；Radontransform；PandSwave；polarizationfiltering；wavefieldseparation第25卷第5期2013年10月岩性油气藏LITHOLOGICRESERVOIRSVol.25No.5Oct.20131，2，1，1，3（1.中国石化石油工程地球物理有限公司胜利分公司，山东东营257086；2.中国石化胜利油田分公司物探研究院，山东东营257000；3.中国石油大学（华东）地球科学与技术学院，山东青岛266555）：以VSP和井间地震勘探为主的井孔地震资料采集需采用三分量地震采集技术，该技术能够记录下信息丰富的地震波场，但是纵波波场和横波波场耦合在一起，成为多分量地震资料处理的重点和难点。通过对井孔地震资料采集观测系统的分析和研究，结合纵波波场和横波波场在不同分量上的关系，可以得到纵波和横波的极性与视速度，以及地震波到达检波器接收排列方向之间的关系。据此在基于Cauchy分布的高分辨率Radon变换的基础上，得到纵波波场和横波波场的分布关系，并在τ（截距时间）-p（视慢度或斜率）域中利用极化滤波方法，实现纵波波场与横波波场的有效分离。通过正演模拟数据试算和实际井间地震资料处理效果分析，认为该方法能够有效地分离纵波波场和横波波场，为后续纵横波的独立成像奠定了基础。：井间地震；Radon变换；纵横波；极化滤波；波场分离：P631.4：A：1673-8926（2013）05-0094-06：2013-03-06；：2013-05-20：“863”“、”（：2011AA060303）。：（1987-），，，。：（257086）。E-mail：[e-mail]2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年0，。，［1-2］。、，，［3-4］。［5］Radon，。［6］，———。［7］，，。［8］，。，［9-12］。［t（）-x（）］，，［13］。Radon，τ（）-p（）。Radon，τ-p，，［14］。τ-p。，τ-p，（P）（S），。τ-p，，，。，CauchyRadon，，τ-p，。，。1CauchyRadonRadon，，pτ。，VSP、，Radon［15］。Bayes，，，Cauchy，，J=μiΣln（1+mk/b）+‖d-Lm‖（1）：mkRadonm；μ（）b（）；d；L。m，LTLm-LTd+Qm=0（2）：Q。Q，，：Qpii=2np×μb+mimax（3）：npp；mimaxnpmi。，，，CauchyRadon。1（a），100，10m，2ms。1（b）Radon。，，。1（c）CauchyRadon。，1（b），CauchyRadon，Radon。22冯玉苹等：井间地震纵横波分离方法研究及应用95岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期2Radon。w（）u（），，。，，［16-17］。，。，uw，（SVSH），，。，，，。2。，（b）Radon（a）（c）CauchyRadon1RadonFig.1ThecontrastofRadontransformresultsofforwardmodelingdata2（P）（SV、SH）Fig.2ThepolarizationdirectionchartofPwaveandconvertedSwave10012141618100.400.20/s20115110115100.400.20/s“”20115110115100.400.20/sSVSSHGPOPSVGSOH0θSVZP“”，。，，3，uw。，962013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年（b）（S）3Fig.3TheupcomingPandSwavepropagatetogeophone（a）（P）4Fig.4ThepolarizationrelationshipofupcomingPandSwavePwu×w-PSu+w-u-w+u-w-u+w+u×w+uSwu2，。，u，w；，u，w。，，uw。，，2，。，uw；，uw。，，uw。，4。，uw，uw。Radon，t-xRadon，τ-p，t-xτ-p。CauchyRadon，τ-p。，Radon，τ-p，，。，，Radon，uw。，，，p（x，t）=sign［up（x，t）］*up2（x，t）+wp2（x，t）姨s（x，t）=sign［us（x，t）］*us2（x，t）+ws2（x，t）姨姨（4）：p；upwpuw；s；uswsuw；x；t。，，（x，y，z）（u，v，w），（SH），。，Radon，（SV）。3uw，CauchyRadon，τ-p，。5（a）5（b）uw，。Radon。CauchyRadon，τ-p，，Radon，。5（c）5（d），uw。冯玉苹等：井间地震纵横波分离方法研究及应用97岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期岩性油气藏第25卷第5期6Fig.6ThePandSwaveseparationresultsofcross-wellseismicdata（c）（a）u（b）w（d）5Fig.5ThePandSwavefieldseparationresultsofmodeldata（d）（a）u（b）w（c）10001400180022002600200300400100014001800220026002003004001000140018002200260020030040010001400180022002600200300400/ms145800015720001686000100200145800015720001686000100200145800015720001686000100200145800015720001686000100200/ms4，Radon。τ-p，Radon，uw，。，，。6（a）6（b）uw，，。Radon。［6（c）6（d）］，，，。5（1）CauchyRadon，，，τ-p。，，。（2）CauchyRadon，τ-p，，，，。（3），Radon982013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年2013年（本文编辑：杨琦）uw，Radon2，t-x，。：本研究得到了中国石油大学（华东）地球科学与技术学院莫言工作室提供的大力帮助，以及中石化石油工程地球物理有限公司胜利分公司地震勘探研究所的鼎力支持，在此一并致谢！：［1］，，.［J］.，2007，19（2）：117-120.［2］WongJ.CrossholeseismicimagingforsulfideorebodydelineationnearSudbury，Ontario，Canada［J］.Geophysics，2000，65（6）：1900-1907.［3］，．［J］．，1995，34（4）：10-13.［4］PrattGR，ShippRM.Seismicwaveforminversioninthefrequencydomain，Part2：Faultdelineationinsedimentsusingcr