基于局部运动学属性参数的窄方位角三维多域插值算法研究与应用
|
摘要
基于局部运动学属性参数的三维多域插值算法在实际地震数据处理中具有广阔的应用前景。但这种算法通常仅适用于宽方位角三维地震数据,而海洋三维地震数据往往是窄方位角的,这就给传统算法的实现带来了一定的难度。为此,提出了一种适应于窄方位角三维地震数据的三维多域插值与规则化新算法,并设计了相应的实现流程。其策略是通过零偏移距三维频率-波数(F-K)域时间偏移与反偏移,结合零偏移距三维输出道成像方式的共反射面元成像算法(CRS-OIS),获得高信噪比的叠后模型道数据;再通过反叠加生成三维叠前CMP模型道道集,提取其运动学属性参数,基于这些运动学属性参数拟合原始三维叠前数据,从而获得信噪比向模型道逼近的优化后叠前数据体。南海深水区某三维窄方位角地震数据的试应用结果表明,新的算法完全适应窄方位角三维地震数据的特点,有效提高了数据的信噪比和规则化程度,证明了新算法的稳定性和可靠性。
The local kinematic attributes based multi-domain interpolation algorithm has huge potential in real applications.However,such kind of algorithm is always not applicable for 3Dnarrow-azimuth seismic data.In this paper,a novel multi-domain interpolation algorithm for 3Dnarrow-azimuth seismic data based on the local kinematic attributes is developed.The basic idea of this algorithm is to combine the zero-offset 3Dfrequency-wavenumber domain time migration/de-migration method and the zero-offset 3D-CRS-OIS method to achieve high signal-noise-ratio zero-offset post-stack model trace data.Then,a 3D model CMP gather can be generated with a smearing style de-stacking operator.The local kinematic attributes of3 D model CMP gathers will be extracted to guide for the optimization of the raw 3DCMP gathers.Thus a practical workflow is designed and applied to deep-water 3Dseismic data acquired in the South China Sea.The application result demonstrated the robustness and high efficiency of the proposed workflow and algorithm.
The local kinematic attributes based multi-domain interpolation algorithm has huge potential in real applications.However,such kind of algorithm is always not applicable for 3Dnarrow-azimuth seismic data.In this paper,a novel multi-domain interpolation algorithm for 3Dnarrow-azimuth seismic data based on the local kinematic attributes is developed.The basic idea of this algorithm is to combine the zero-offset 3Dfrequency-wavenumber domain time migration/de-migration method and the zero-offset 3D-CRS-OIS method to achieve high signal-noise-ratio zero-offset post-stack model trace data.Then,a 3D model CMP gather can be generated with a smearing style de-stacking operator.The local kinematic attributes of3 D model CMP gathers will be extracted to guide for the optimization of the raw 3DCMP gathers.Thus a practical workflow is designed and applied to deep-water 3Dseismic data acquired in the South China Sea.The application result demonstrated the robustness and high efficiency of the proposed workflow and algorithm.
引文
[1]Xu S,Zhang Y,LambaréG.Antileakage Fourier transform for seismic data regularization in higher dimensions[J].Geophysics,2010,75(6):WB113-WB120
[2]Naghizadeh M,Sacchi M.Multidimensional de-aliased Cadzow reconstruction of seismic records[J].Geophysics,2013,78(1):A1-A5
[3]Wang X,Wang H Z,Zhou D H,et al.Minimizing ellipsoidal norm seismic data interpolation with radon spectrum constraints[J].Expanded Abstracts of 73rd EAGE Annual Conference,2011,A046
[4]辛可锋,王华忠,王成礼,等.叠前地震数据的规则化[J].石油地球物理勘探,2002,37(4):311-317Xin K F,Wang H Z,Wang C L,et al.Regularization of pre-stack seismic data[J].Oil Geophysical Prospecting,2002,37(4):311-317
[5]刘玉金,李振春.局部平面波模型约束下的迭代加权最小二乘反演三维地震数据规则化[J].石油地球物理勘探,2012,47(3):418-424Liu Y J,Li Z C.Three dimensional data regularization by iteratively re-weighted least squares inversion based on the local plane wave model[J].Oil Geophysical Prospecting,2012,47(3):418-424
[6]杨锴,王毓玮,王华忠,等.基于偏移/反偏移的地震数据映射方法[J].石油地球物理勘探,2008,43(1):16-22Yang K,Wang Y W,Wang H Z,et al.Seismic data mapping based on migration/inverse-migration[J].Oil Geophysical Prospecting,2008,43(1):16-22
[7]Buzlukov V,Landa E.Imaging improvement by prestack signal enhancement[J].Geophysical Prospecting,2013,61(6):1150-1158
[8]Hoecht G,Ricarte P,Bergler S,et al.Operator-oriented CRS interpolation[J].Geophysical Prospecting,2009,57(6):957-981
[9]Yang K,Wang H Z,Ma Z T.An output imaging scheme of the common reflection surface stack[J].Journal of Seismic Exploration,2005,14(1):131-154
[10]Baykulov M,Gajewski D.Prestack seismic data enhancement with partial common-reflection-surface(CRS)stack[J].Geophysics,2009,74(3):V49-V58
[11]Yang K,Wang H Z,Dong L G,et al.An output imaging scheme of the common reflection surface stack:applications to real data[J].Expanded Abstracts of76th Annual Internat SEG Mtg,2006,2529-2533
[12]Yang K,Jiang X G.Application of an output imaging scheme of 3Dcommon reflection surface stack(3D CRS-OIS)to 3Dreal data acquired in Chinese mountainous area[J].Expanded Abstracts of 80th Annual Internat SEG Mtg,2010,3614-3618
[13]Ni Y,Yang K,Wang Y X.Noise suppression and data interpolation with 3DCRS-OIS[J].Expanded Abstracts of 82nd Annual Internat SEG Mtg,2012,2529-2533
[14]杨锴,倪瑶.基于稳相原理与克希霍夫统一成像理论对射线束类成像方法的一种新观点:以CRS叠加为例[C]∥CPS-SEG北京2014国际地球物理会议摘要.北京:石油出版社,2014:490-493Yang K,Ni Y.A new insight of beam-ray imaging method based on the stationary-phase principle and unified Kirchhoff imaging theory:for example of CRS stack[C]∥CPS-SEG Beijing 2014International Geophysical Conference.Beijing:Petroleum Press,2014:490-493
[15]童思友,杨锴,孔剑冰,等.利用CRS-OIS实现噪声压制与数据插值[J].同济大学学报(自然科学版).2010,38(5):753-757Tong S Y,Yang,K,Kong J B,et al.Noise suppression and data interpolation by using CRS-OIS[J].Journal of Tongji University(Natural Science),2010,38(5):753-757
[16]倪瑶,杨锴.基于GPU计算平台实现三维输出道方式的共反射面元(3D-CRS-OIS)叠加[J].石油地球物理勘探,2013,48(1):57-65Ni Y,Yang,K.3Dcommon reflection surface stack algorithm based GPU with the output imaging scheme(3D-CRS-OIS)[J].Oil Geophysical Prospecting,2013,48(1):57-65
[17]刘春成,刘畅,王大利,基于重构信号模型的去噪技术[J].铀矿地质,2011,27(3):173-179Liu C C,Liu C,Wang D L.A new technology for noise removal in the model-based signal reconstruction[J].Uranium Geology,2011,27(3):173-179
[18]Dubrulle A.Numerical methods for the migration of constant-offset sections in homogeneous and horizontally layered media[J].Geophysics,1983,48(9):1195-1203
[19]Yang K,Chen B S,Wang X J,et al,Handling conflicting dip problem in CRS stack via combination of CRS-OIS and migration/demigration[J].Geophysical Prospecting,2012,60(2):255-269
[2]Naghizadeh M,Sacchi M.Multidimensional de-aliased Cadzow reconstruction of seismic records[J].Geophysics,2013,78(1):A1-A5
[3]Wang X,Wang H Z,Zhou D H,et al.Minimizing ellipsoidal norm seismic data interpolation with radon spectrum constraints[J].Expanded Abstracts of 73rd EAGE Annual Conference,2011,A046
[4]辛可锋,王华忠,王成礼,等.叠前地震数据的规则化[J].石油地球物理勘探,2002,37(4):311-317Xin K F,Wang H Z,Wang C L,et al.Regularization of pre-stack seismic data[J].Oil Geophysical Prospecting,2002,37(4):311-317
[5]刘玉金,李振春.局部平面波模型约束下的迭代加权最小二乘反演三维地震数据规则化[J].石油地球物理勘探,2012,47(3):418-424Liu Y J,Li Z C.Three dimensional data regularization by iteratively re-weighted least squares inversion based on the local plane wave model[J].Oil Geophysical Prospecting,2012,47(3):418-424
[6]杨锴,王毓玮,王华忠,等.基于偏移/反偏移的地震数据映射方法[J].石油地球物理勘探,2008,43(1):16-22Yang K,Wang Y W,Wang H Z,et al.Seismic data mapping based on migration/inverse-migration[J].Oil Geophysical Prospecting,2008,43(1):16-22
[7]Buzlukov V,Landa E.Imaging improvement by prestack signal enhancement[J].Geophysical Prospecting,2013,61(6):1150-1158
[8]Hoecht G,Ricarte P,Bergler S,et al.Operator-oriented CRS interpolation[J].Geophysical Prospecting,2009,57(6):957-981
[9]Yang K,Wang H Z,Ma Z T.An output imaging scheme of the common reflection surface stack[J].Journal of Seismic Exploration,2005,14(1):131-154
[10]Baykulov M,Gajewski D.Prestack seismic data enhancement with partial common-reflection-surface(CRS)stack[J].Geophysics,2009,74(3):V49-V58
[11]Yang K,Wang H Z,Dong L G,et al.An output imaging scheme of the common reflection surface stack:applications to real data[J].Expanded Abstracts of76th Annual Internat SEG Mtg,2006,2529-2533
[12]Yang K,Jiang X G.Application of an output imaging scheme of 3Dcommon reflection surface stack(3D CRS-OIS)to 3Dreal data acquired in Chinese mountainous area[J].Expanded Abstracts of 80th Annual Internat SEG Mtg,2010,3614-3618
[13]Ni Y,Yang K,Wang Y X.Noise suppression and data interpolation with 3DCRS-OIS[J].Expanded Abstracts of 82nd Annual Internat SEG Mtg,2012,2529-2533
[14]杨锴,倪瑶.基于稳相原理与克希霍夫统一成像理论对射线束类成像方法的一种新观点:以CRS叠加为例[C]∥CPS-SEG北京2014国际地球物理会议摘要.北京:石油出版社,2014:490-493Yang K,Ni Y.A new insight of beam-ray imaging method based on the stationary-phase principle and unified Kirchhoff imaging theory:for example of CRS stack[C]∥CPS-SEG Beijing 2014International Geophysical Conference.Beijing:Petroleum Press,2014:490-493
[15]童思友,杨锴,孔剑冰,等.利用CRS-OIS实现噪声压制与数据插值[J].同济大学学报(自然科学版).2010,38(5):753-757Tong S Y,Yang,K,Kong J B,et al.Noise suppression and data interpolation by using CRS-OIS[J].Journal of Tongji University(Natural Science),2010,38(5):753-757
[16]倪瑶,杨锴.基于GPU计算平台实现三维输出道方式的共反射面元(3D-CRS-OIS)叠加[J].石油地球物理勘探,2013,48(1):57-65Ni Y,Yang,K.3Dcommon reflection surface stack algorithm based GPU with the output imaging scheme(3D-CRS-OIS)[J].Oil Geophysical Prospecting,2013,48(1):57-65
[17]刘春成,刘畅,王大利,基于重构信号模型的去噪技术[J].铀矿地质,2011,27(3):173-179Liu C C,Liu C,Wang D L.A new technology for noise removal in the model-based signal reconstruction[J].Uranium Geology,2011,27(3):173-179
[18]Dubrulle A.Numerical methods for the migration of constant-offset sections in homogeneous and horizontally layered media[J].Geophysics,1983,48(9):1195-1203
[19]Yang K,Chen B S,Wang X J,et al,Handling conflicting dip problem in CRS stack via combination of CRS-OIS and migration/demigration[J].Geophysical Prospecting,2012,60(2):255-269
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心