主动源与被动源地震数据插值及联合数据成像
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
本文提出了两种情况下主动源数据和被动源数据的插值方法,并研究了两种数据在偏移成像中的互补效果.基于互相关法被动源数据重构原理,本文提出了结合多域迭代去噪技术的重构方法.提出了两种时间域主动源和被动源数据的插值方法,分别是共炮点域能量匹配插值和共检波点域最小平方匹配插值.然后对获得的主动源和被动源联合地震数据进行叠前深度偏移成像.在被动源活跃度不是很高的地区进行被动源地震勘探时,少量的主动源地震数据可以有效控制和补充被动源数据的成像效果.在稀疏炮点的主动源勘探中,有效利用被动源的信息能够在成像中增加更多的细节信息,提高成像质量.
In seismic acquisition,we may only acquire a limited amount of data or lose parts of data because of the environment,cost or efficiency.This will inevitably influence the quality of final imaging results.In data processing,we usually use interpolation to solve the problem of sparse data.The conventional interpolation methods are usually based on mathematic principles,and compensate data after calculation.The passive sources are actually located in the subsurface.The passive seismic waves transmit to the surface according to the same principles as active seismic waves.As a result,they contain abundant real information about the subsurface.We consider compensating active source data with passive source data using interpolation methods.In addition,in the areas suitable for passive seismic exploration,when there are not enough sources,the quality of reconstruction is low.So we also consider compensating passive source data with active source data using interpolation methods.As the SNR of passive retrieved data reconstructed by the conventional method is low,the quality of imaging is influenced seriously.Based on conventional methods,we propose amultidomain iterative denoising reconstruction method,which is based on the Curvelet threshold iterative denoising and multilevel median filter.As to the two cases mentioned above,and considering the matching degree of passive virtual wavelets and active wavelets,we suggest two time-domain interpolation methods using active and passive seismic data.They are the energy matching interpolation in the common source domain and least square matching interpolation in the common receiver domain.When the two kinds of wavelets match well,we adopt the energy matching method to calculate the interpolation factors.When the two kinds of wavelets do not match well,we adopt the least square matching interpolation and linear interpolation methods to calculate matching factors.Using the above reconstruction and interpolation methods,we analyze three cases and compare the results before and after pre-stack depth migration.For the ideal case,we adopt 1000 random noise sources,which are distributed uniformly both in axial and lateral directions.The recording time is 1200 s.After reconstruction,we conduct pre-stack depth migration with active and passive source data.The results show that in the ideal case,there is a large amount of effective reflection information.In the imaging results,the imaging result of passive-source seismic data is very close to that of active-source data.The complex structures and reservoirs are displayed clearly.For the case with inadequate passive sources,we adopt 400 locally distributed noise sources. We estimate the source region by comparing the causal and un-causal parts of reconstruction results.Then we conduct interpolation and migration.We can see from the results that the deep detailed information is clearer after interpolation.For the case of sparse data in active-source seismic exploration,we use a few noise sources located at depth.From the interpolation results,we can see that the active-source data and passive-source data match well after matching.The imaging results show that the complex structures and deep information become clearer after interpolation.The numerical results show that effective use of active-and passive-source data can greatly improve the quality of imaging.The multidomain iterative denoising reconstruction method can suppress noises effectively.By adjusting the threshold value and filter windows,we can extract effective signals gradually.In the interpolation processing of active-and passive-source data,energy matching interpolation in the common source domain and least square matching interpolation in the common receiver domain can obtain good results in different cases.In the area with a few passive sources,we can estimate the source region from the causal and un-causal part of data.Adding a few active data in the quiet area can compensate the imaging results effectively.In the sparse shot active seismic exploration,the use of passive data can add some detailed information.
引文
Artman B,Draganov D,Wapenaar K,et al.2004.Direct Migrationof passive seismic data.∥66th EAGE Conference and ExhibitionExtended Abstracts,Z-99.
    Bakulin A,Calvert R.2006.The virtual source method:Theoryand case study.Geophysics,71(4):SI139-SI150.
    Berkhout A J,Verschuur D J.2011.A scientific framework foractive and passive seismic imaging,with applications to blendeddata and micro-earthquake responses.Geophysical JournalInternational,184(2):777-792,doi:10.1111/j.1365-246X.2010.04855.x.
    Claerbout J F.1968.Synthesis of a layered medium from its acoustictransmission response.Geophysics,33(2):264-269.
    Doulgeris P,Mahdad A,Blacquiere G.2010.Separation of blendedimpulsive sources using an iterative approach.∥72nd EAGEConference and Exhibition Extended Abstracts,B004.
    Draganov D,Wapenaar K,Thorbecke J.2006.Seismic interferometry:Reconstructing the earth′s reflection response.Geophysics,71(4):SI61-SI70.
    Draganov D,Wapenaar K,Mulder W,et al.2007.Retrieval ofreflections from seismic background-noise measurements.Geophysical Research Letters,34:L04305,doi:10.1029/2006GL028735.
    Ferguson R J,Margrave G F.2005.Planned seismic imaging usingexplicit one-way operators.Geophysics,70(5):S101-S109.
    Groenestijn G V,Verschuur D J.2010.Estimation of primaries bysparse inversion from passive seismic data.Geophysics,75(4):SA61-SA69.
    Han L G,Tan C Q,LüQ T,et al.2013.Separation of multisource blended seismic acquisition data by iterative denoising.Chinese J.Geophys.(in Chinese),56(7):2402-2412,doi:10.6038/cjg20130726.
    Qi C,Chen Q F,Chen Y.2007.A new method for seismic imagingfrom ambient seismic noise.Progress in Geophysics(in Chinese),22(3):771-777.
    Ramirez A C,Weglein A B.2009.Greens theorem as acomprehensive framework for data reconstruction,regularization,wavefield separation,seismic interferometry,and wavelet estimation:A tutorial.Geophysics,74(6):W35-W62,doi:10.1190/1.3237118.
    Ruigrok E N,Vidal C A.2013.Body-wave Receiver-pair SeismicInterferometry.∥75th EAGE Conference and ExhibitionExtended Abstracts,Th-P15-02.
    Schuster G T,Rickett J.2001.Daylight imaging in V(x,y,z)media.∥Stanford Exploration Project,Report SEP-105,209-227.
    Schuster G T,Yu J,Sheng J,et al.2004.Interferometric/daylightseismic imaging.Geophysics Journal International,157(2):838-852.
    Sheng J.2001.Migrating multiples and primaries in CDP data bycross-correlation migration.∥71st SEG Annual meetingExpanded Abstracts,1297-1300.
    Snieder R,Miyazawa M,Slob E,et al.2009.A comparison ofstrategies for seismic interferometry.Surveys in Geophysics,30(4-5):503-523,doi:10.1007/s10712-009-9069-z.
    Stehly L,Campillo M,Shapiro N M.2006.A study of the seismicnoise from its long-range correlation properties.Journal ofGeophysical Research,111(B10306):1-12,doi:10.1029/2005JB004237.
    Tan C Q,Han L G,Zhang Y H,et al.2012.Separation of blended
    data by iterative denoising.∥74th EAGE Conference andExhibition Extended Abstracts,A045.
    Thorbecke J W,Draganov D.2011.Finite-difference modelingexperiments for seismic interferometry.Geophysics,76(6):H1-H18,doi:10.1190/GEO2010-0039.1.
    Wapenaar K.2003.Synthesis of an inhomogeneous medium from itsacoustic transmission response.Geophysics,68(5):1756-1759.
    Wapenaar K.2004.Retrieving the elastodynamic Green′s functionof an arbitrary inhomogeneous medium by cross correlations.Physical Review Letters,93:1-4,254301,doi:10.1103/PhysRevLett.93.254301.
    Wapenaar K,Thorbecke J,Draganov D.2004.Relations betweenreflection and transmission responses of 3-D inhomogeneousmedia.Geophysical Journal International,156(2):179-194.
    Wapenaar K,Fokkema J.2006.Green′s function representationsfor seismic interferometry.Geophysics,71(4):SI33-SI46,doi:10.1190/1.2213955.
    Wapenaar K,Neut J,Ruigrok E.2008.Passive seismic interferometryby multidimensional deconvolution.Geophysics,73(6):51-56.
    Xu Z,Juhlin C,Gudmundsson O,et al.2012.Reconstruction ofsubsurface structure from ambient seismic noise:an examplefrom Ketzin,Germany.Geophysical Journal International,189(2):1085-1102.
    Zheng Y C.2010.Retrieving the exact Green′s function by wavefieldcrosscorrelation.J.Acoust.Soc.Am.,127(3):EL93-EL98.
    Zheng Y C,He Y F,Fehler M C.2011.Crosscorrelation kernels inacoustic Green′s function retrieval by wavefield correlation forpoint sources on a plane and a sphere.Geophysical JournalInternational,184(2):853-859.
    韩立国,谭尘青,吕庆田等.2013.基于迭代去噪的多源地震混合采集数据分离.地球物理学报,56(7):2402-2412,doi:10.6038/cjg20130726.
    齐诚,陈棋福,陈颙.2007.利用背景噪声进行地震成像的新方法.地球物理学进展,22(3):771-777.

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心