复杂区低信噪比地震资料处理方法研究
|
摘要
本文针对目前复杂区低信噪比地震资料处理面临的主要瓶颈技术问题,以提高叠前资料的信噪比、补偿有效信号能量和保证复杂构造精确成像为目的,主要通过对叠前散射噪音去除、多次波衰减、透射损失补偿和消除非双曲时差影响的高精度速度分析以及考虑各向异性和地形起伏影响的偏移成像等处理方法的研究和探讨。经过研究提出了基于信号子空间分解的空间变换叠前散射噪音衰减方法;频率域自适应抛物线拉东变换多次波衰减;基于波动理论地震信号透射损失补偿;高密度逐点双谱速度分析和基于起伏地表的各向异性叠前时间偏移五项处理技术方法,归纳总结出一套适合复杂区低信噪比地震资料处理的技术思路和流程。
通过对这些方法的原理进行详细的阐述、理论模型测试并对大庆外围盆地三江绥滨二维和方正复杂区三维实际地震数据进行处理,最终形成了可以进行实际资料处理的技术方法和流程,其实际处理结果与目前常规的处理方法结果进行对比:可以明显地看出采用本文所研究的方法处理结果在资料的信噪比、有效信号的能量、保真度以及成像精度等方面都明显优于目前常规处理方法的结果。说明了该方法的先进性、有效性和适用性。为复杂区低信噪比地震资料处理提供了行之有效的方法和可借鉴的思路和技术流程。
In recent years, with the advance of exploration technology and the continuous rising level of exploration, the easily found structural traps with simple geological conditions are becoming less and less and the complex oil and gas reservoirs in the complex areas are becoming more and more explored. These complex exploration objects present new challenges to the oil and gas geophysical technology with seismic exploration as the main method. Especially in the area of low S/N ration seismic data processing there are still many bottleneck techniques hard to be solved. Firstly, the drastically changed surface relief and the complicated surface velocities caused by the effects of complex surface geological conditions, together with the complexity of surface heterogeneity, lead to the original seismic data acquired with low S/N ratio and serious interference and in the mean time, the surface waves, the refraction waves and especially the comparatively evolutional scattering noises present many difficulties for raising the S/N ratio. Secondly, because of the effects of complex underground seismic and geological conditions, the data structures from complex areas are extremely complicated with developed fractures, complicated wave fields and very changed lateral velocities. Furthermore, the igneous rocks are comparatively developed and widely distributed. Due to the presence of igneous rocks, on the one hand, the igneous rocks have very strong energy and the strong reflection energy has very good isolation for the underlying effective reflection. On the other hand, the velocity is high and has very big wave impedance difference with surrounding rocks. Thus, the strong multiple waves are produced under the reflection interface. The existence of these factors makes the already complex wave fields much complicated and the accurate imaging of structures difficult.
Based on the actual S/N ratio data from complex areas in Fangzheng Fault of Sanjiang and Yishu Graben of Peripheral Basin of Daqing Oilfield, the study in this paper aims at the methods for suppressing the pre-stack noises, especially the strong scattering and multiple-wave interference to improve the S/N ratio of actual data and ensure the accurate imaging of complex structures. Through method argumentation of principles, theoretical model design, software programming, as well as testing of actual data and final practical application of actual data and based on the results of previous researches, this paper puts forward some new innovative technical methods and understanding and has obtained the following five technological achievements suitable for low S/N ratio seismic data processing in complex areas.①Pre-stack scattering noise attenuation technique based on the spatial conversion of signal sub-space decomposition.②Multiple-wave removal technique with self-adaptive parabolic Radon transformation in the frequency domain.③Compensation technique for seismic signal transmission loss based on the wave theory.④High-density and point by point dual spectral velocity analysis technique.⑤Pre-stack time migration processing technique based on the anisotropy of surface relief. After making the systematic summary and conclusion as well as conducting the actual data processing, the technological process framework suitable for low S/N ratio data processing in complex areas has been obtained, which provides effective methods, ideas for reference and technological process for low S/N ratio data processing in complex areas. The work done mainly focuses on the following areas.
1. By the analysis of formation mechanism and forward modeling for scattering noises, the characteristics and formation reasons of scattering noises have been obtained. The display forms of scattering noises in original records are generally linear. However, when the scattering anomaly is buried very deep, the similar hyperbolic characteristic, which has relatively strong energy and wide frequency distribution, will be displayed. It is very difficult to be removed by current conventional methods. As to the characteristics of scattering noises, this paper introduces the pre-stack scattering noise attenuation method based on the spatial conversion of signal sub-space decomposition. The thinking of this method is: to use the characteristic that the waveform of scattering noises has a certain similarity in the time-space domain for simulating and predicting the scattering noises and getting the estimation of scattering wave field and then remove the scattering wave filed from the seismic records and realize the attenuation of seismic scattering noises. This paper conducts detailed analysis and argumentation of this method on basic principles, theoretical model design and actual data processing and finally puts forward the pre-stack scattering noise attenuation technique based on the space conversion of the signal sub-space decomposition. With the help of corresponding developed software, this technique has been loaded to CGG processing system on PC clusters and can be mass-produced. It is a new method for pre-stack scattering noise suppression.
2. Based on the review of previous Radon transformation technique for the removal of multiple waves, this paper also focuses on the characteristics and existing problems of Radon transformation in the time domain. Through the analysis of basic formulas and physical meaning of Radon transformation as well as multiple tests for actual multiple-wave separation with the theoretical model data, it is concluded that although PRT in the time domain has a clear physical and intuitive meaning, the unavoidable shortcomings of Radon transformation in the time domain still exist when conducting actual programming. After having found the drawbacks of Radon transformation in the time domain, this paper presents a new method that uses the parabolic Radon transformation to remove multiple waves in the frequency domain. Through different implementations of Radon transformation method in the frequency domain, inverse transformation least-square algorithm and forward transformation least-square algorithm, as well as comparison of results obtained by these two methods, it is concluded that the forward transformation least-square algorithm can effectively focus the energy and be convenient to distinguish and remove multiple waves and also prevent the presence of alias in the q-ωdomain. In the meanwhile, in order to further enhance the effects of multiple-wave removal and improve the accuracy of algorithm, this paper also focuses on discussing the application factors such as aliasing problems in the Radon domain, solution of double coefficient matrix equation, selection of white noise coefficient, removal of multiple waves in the Radon domain and self-adaptive filtering, etc. and makes a conclusion with a certain guiding significance. And then, after making a summary and conclusion, this paper puts forward the concrete realization process and procedures with the self-adaptive parabolic Radon transformation to remove multiple waves in the frequency domain. And with the characteristics of ultra low S/N ratio data from Suibin Depression in Sanjiang Basin in consideration, the Radon transformation method in the frequency domain was used to process the actual data and the results are very good.
3. This paper also emphasizes the necessity of compensation for transmission loss, which is one factor that cannot be ignored for achieving high-fidelity amplitude compensation. Especially when the objective layer is covered by a large number of strong reflection interfaces or there is the presence of volcanic rocks with strong energy, the energy of effective reflection objective layer is very weak because of the severe transmission loss, which decreases the interpretation accuracy of final results. This paper points out the existing problems on the current conventional compensation method for transmission loss by analyzing the previous research results. The previous studies focused on the transmission loss of seismic waves in the random medium and only have theoretical significance. But the actual medium in sedimentation process is not at random and has cyclic characteristic. Thus, the transmission coefficient calculation in the actual non-random medium must be studied. In order to solve this problem, this paper puts forward the compensation technical method for seismic signal transmission loss based on the wave theory. Firstly, use the geometric theory and wave theory to calculate respectively the transmission coefficient for the thick layer model and thin layer model, describe the difference of physical nature between the two methods and then point out that the wave theory should be used to calculate the transmission coefficient of actual thin-layered medium. Secondly, establish the transmission compensation method based on the wave theory, use the modeling technique to realize the quantitative analysis of transmission loss mechanism, use the wave theory to calculate the compensation factor of transmission loss, point out the processing procedures for realizing post-stack transmission compensation with the velocity spectrum and seismic traces as the input data and then achieve the true restoration of seismic wave energy. Based on a large number of researches on theoretical methods and after making the summary and conclusion, this paper gives the concrete principles and technological process for post-stack transmission compensation. This compensation method was used to process the actual data from Suibin Depression in Sanjiang Basin and the result is very remarkable. The continuity on same phase axis in many parts of transmission-compensated sections is improved after the amplitude of reflection wave being restored and the basement as well as the form of deep strata is clear. Therefore, the post-stack transmission compensation is useful.
4. This part is about the conventional velocity analysis method for the data from complex areas based on the time-offset hyperbolic curve of seismic reflection wave travel. Because of the approximation of hyperbolic hypothesis, the conventional velocity analysis cannot meet the requirement for establishing subsurface accurate velocity model with the data from complex areas. This paper puts forward the high-density and point by point dual spectral velocity analysis technique on the basis of previous research results. By introducing principles of high-density and point by point dual spectral velocity analysis, dynamic correction parameters to affect all the offsets, intelligent sorting of the pairν,η, geological filtering of high-density time differential correction parameter field (ν,η), unrelated dynamic correction parameters (dtn,τ0), interpolation and filtering ofνandηfields, this paper presents the theoretical formula of high-order non-hyperbolic time difference and puts forward the processing procedures used for high-density and point by point velocity analysis. This method uses the bi-directional scanning of dual-time parameters to obtain accurate velocity parameters. The actual processing of data acquired from Fangzheng Fault in Daqing Peripheral Basin demonstrates that this method can effectively solve the problem that the large-offset non-hyperbolic time difference correction gather couldn’t be leveled. Meanwhile, the point by point velocity analysis and calculation effectively solve the problem of lateral change of velocity and thus the accurate velocity fields of underground complex structures can be obtained.
5. In order to ensure the accuracy of migration imaging in complex areas, this paper analyzes the main problems faced by migration imaging in complex areas and puts forward the pre-stack time migration processing technique based on the anisotropy of surface relief. On the basis of imaging principles of surface common-reflection point, this paper uses the equivalent offset method for the data from the terrain with surface relief to eliminate the effect of surface relief on the structural morphologies and then formulates in detail the anisotropic pre-stack time migration processing technique. With the introduction of important concepts of anisotropy, phase velocity and ray velocity, NMO velocity in anisotropic medium, three dimensional description of NMO velocity in the anisotropic medium and principles of anisotropic pre-stack time migration method, this paper focuses on the anisotropic pre-stack time migration velocity analysis, high-accuracy calculation of velocity field, calculation of anisotropic parameters in VTI medium and the implementation process of anisotropic pre-stack time migration, and finally gets high-accuracy velocity modeling method and technical process for three dimensional anisotropic pre-stack time migration,ηfield method and process for obtaining anisotropic parameters in VTI medium with surface seismic data, three dimensional anisotropic pre-stack time migration processing method and procedures. It is proved by actual example that the surface–based anisotropic time migration processing technique can eliminate the effects of anisotropy on seismic wave travel path and velocity as well as the effects of surface relief on seismic data imaging and thus improves the accuracy of migration imaging. At the same time, because of the elimination of every anisotropic effect, it also solves the problem that the geophone offset on the CRP gather is over-corrected and in the meanwhile increases the effective multiplicity used for pre-stack inversion gather. And thus, the S/N ratio of migration result can be improved and it is also favorable for AVO analysis and pre-stack elastic impedance inversion.
At the end of this paper, the five techniques used for processing seismic data from complex areas, obtained from this research, are systematically summarized and concluded and a set of technological processes, suitable for processing low S/N ratio data from complex surfaces and complex structures, is given. And then, the actual data of 100 km2 from Fangzheng Fault in Daqing Peripheral Basin was processed with the above techniques. By comparing the overall effect of the final results, the new processed results with the previous 2D and 3D results, analyzing quantitatively the S/N ratios of new and old data and comparing the new interpreted results with the old one, it is demonstrated that not only the S/N ration but also the resolution of data are improved significantly, especially the structural morphologies have very big changes. After comparing with the exploration well data, the new-processed results match perfectly with the well data, which proves that the processing methods in this paper are advanced and practical. This paper also provides a set of new techniques for processing the low S/N ratio data from complex areas, which can be very reliable technical support for accurate structure interpretation, inversion prediction of reservoirs and for speeding up the seismic exploration of complex surfaces and complex structures.
通过对这些方法的原理进行详细的阐述、理论模型测试并对大庆外围盆地三江绥滨二维和方正复杂区三维实际地震数据进行处理,最终形成了可以进行实际资料处理的技术方法和流程,其实际处理结果与目前常规的处理方法结果进行对比:可以明显地看出采用本文所研究的方法处理结果在资料的信噪比、有效信号的能量、保真度以及成像精度等方面都明显优于目前常规处理方法的结果。说明了该方法的先进性、有效性和适用性。为复杂区低信噪比地震资料处理提供了行之有效的方法和可借鉴的思路和技术流程。
In recent years, with the advance of exploration technology and the continuous rising level of exploration, the easily found structural traps with simple geological conditions are becoming less and less and the complex oil and gas reservoirs in the complex areas are becoming more and more explored. These complex exploration objects present new challenges to the oil and gas geophysical technology with seismic exploration as the main method. Especially in the area of low S/N ration seismic data processing there are still many bottleneck techniques hard to be solved. Firstly, the drastically changed surface relief and the complicated surface velocities caused by the effects of complex surface geological conditions, together with the complexity of surface heterogeneity, lead to the original seismic data acquired with low S/N ratio and serious interference and in the mean time, the surface waves, the refraction waves and especially the comparatively evolutional scattering noises present many difficulties for raising the S/N ratio. Secondly, because of the effects of complex underground seismic and geological conditions, the data structures from complex areas are extremely complicated with developed fractures, complicated wave fields and very changed lateral velocities. Furthermore, the igneous rocks are comparatively developed and widely distributed. Due to the presence of igneous rocks, on the one hand, the igneous rocks have very strong energy and the strong reflection energy has very good isolation for the underlying effective reflection. On the other hand, the velocity is high and has very big wave impedance difference with surrounding rocks. Thus, the strong multiple waves are produced under the reflection interface. The existence of these factors makes the already complex wave fields much complicated and the accurate imaging of structures difficult.
Based on the actual S/N ratio data from complex areas in Fangzheng Fault of Sanjiang and Yishu Graben of Peripheral Basin of Daqing Oilfield, the study in this paper aims at the methods for suppressing the pre-stack noises, especially the strong scattering and multiple-wave interference to improve the S/N ratio of actual data and ensure the accurate imaging of complex structures. Through method argumentation of principles, theoretical model design, software programming, as well as testing of actual data and final practical application of actual data and based on the results of previous researches, this paper puts forward some new innovative technical methods and understanding and has obtained the following five technological achievements suitable for low S/N ratio seismic data processing in complex areas.①Pre-stack scattering noise attenuation technique based on the spatial conversion of signal sub-space decomposition.②Multiple-wave removal technique with self-adaptive parabolic Radon transformation in the frequency domain.③Compensation technique for seismic signal transmission loss based on the wave theory.④High-density and point by point dual spectral velocity analysis technique.⑤Pre-stack time migration processing technique based on the anisotropy of surface relief. After making the systematic summary and conclusion as well as conducting the actual data processing, the technological process framework suitable for low S/N ratio data processing in complex areas has been obtained, which provides effective methods, ideas for reference and technological process for low S/N ratio data processing in complex areas. The work done mainly focuses on the following areas.
1. By the analysis of formation mechanism and forward modeling for scattering noises, the characteristics and formation reasons of scattering noises have been obtained. The display forms of scattering noises in original records are generally linear. However, when the scattering anomaly is buried very deep, the similar hyperbolic characteristic, which has relatively strong energy and wide frequency distribution, will be displayed. It is very difficult to be removed by current conventional methods. As to the characteristics of scattering noises, this paper introduces the pre-stack scattering noise attenuation method based on the spatial conversion of signal sub-space decomposition. The thinking of this method is: to use the characteristic that the waveform of scattering noises has a certain similarity in the time-space domain for simulating and predicting the scattering noises and getting the estimation of scattering wave field and then remove the scattering wave filed from the seismic records and realize the attenuation of seismic scattering noises. This paper conducts detailed analysis and argumentation of this method on basic principles, theoretical model design and actual data processing and finally puts forward the pre-stack scattering noise attenuation technique based on the space conversion of the signal sub-space decomposition. With the help of corresponding developed software, this technique has been loaded to CGG processing system on PC clusters and can be mass-produced. It is a new method for pre-stack scattering noise suppression.
2. Based on the review of previous Radon transformation technique for the removal of multiple waves, this paper also focuses on the characteristics and existing problems of Radon transformation in the time domain. Through the analysis of basic formulas and physical meaning of Radon transformation as well as multiple tests for actual multiple-wave separation with the theoretical model data, it is concluded that although PRT in the time domain has a clear physical and intuitive meaning, the unavoidable shortcomings of Radon transformation in the time domain still exist when conducting actual programming. After having found the drawbacks of Radon transformation in the time domain, this paper presents a new method that uses the parabolic Radon transformation to remove multiple waves in the frequency domain. Through different implementations of Radon transformation method in the frequency domain, inverse transformation least-square algorithm and forward transformation least-square algorithm, as well as comparison of results obtained by these two methods, it is concluded that the forward transformation least-square algorithm can effectively focus the energy and be convenient to distinguish and remove multiple waves and also prevent the presence of alias in the q-ωdomain. In the meanwhile, in order to further enhance the effects of multiple-wave removal and improve the accuracy of algorithm, this paper also focuses on discussing the application factors such as aliasing problems in the Radon domain, solution of double coefficient matrix equation, selection of white noise coefficient, removal of multiple waves in the Radon domain and self-adaptive filtering, etc. and makes a conclusion with a certain guiding significance. And then, after making a summary and conclusion, this paper puts forward the concrete realization process and procedures with the self-adaptive parabolic Radon transformation to remove multiple waves in the frequency domain. And with the characteristics of ultra low S/N ratio data from Suibin Depression in Sanjiang Basin in consideration, the Radon transformation method in the frequency domain was used to process the actual data and the results are very good.
3. This paper also emphasizes the necessity of compensation for transmission loss, which is one factor that cannot be ignored for achieving high-fidelity amplitude compensation. Especially when the objective layer is covered by a large number of strong reflection interfaces or there is the presence of volcanic rocks with strong energy, the energy of effective reflection objective layer is very weak because of the severe transmission loss, which decreases the interpretation accuracy of final results. This paper points out the existing problems on the current conventional compensation method for transmission loss by analyzing the previous research results. The previous studies focused on the transmission loss of seismic waves in the random medium and only have theoretical significance. But the actual medium in sedimentation process is not at random and has cyclic characteristic. Thus, the transmission coefficient calculation in the actual non-random medium must be studied. In order to solve this problem, this paper puts forward the compensation technical method for seismic signal transmission loss based on the wave theory. Firstly, use the geometric theory and wave theory to calculate respectively the transmission coefficient for the thick layer model and thin layer model, describe the difference of physical nature between the two methods and then point out that the wave theory should be used to calculate the transmission coefficient of actual thin-layered medium. Secondly, establish the transmission compensation method based on the wave theory, use the modeling technique to realize the quantitative analysis of transmission loss mechanism, use the wave theory to calculate the compensation factor of transmission loss, point out the processing procedures for realizing post-stack transmission compensation with the velocity spectrum and seismic traces as the input data and then achieve the true restoration of seismic wave energy. Based on a large number of researches on theoretical methods and after making the summary and conclusion, this paper gives the concrete principles and technological process for post-stack transmission compensation. This compensation method was used to process the actual data from Suibin Depression in Sanjiang Basin and the result is very remarkable. The continuity on same phase axis in many parts of transmission-compensated sections is improved after the amplitude of reflection wave being restored and the basement as well as the form of deep strata is clear. Therefore, the post-stack transmission compensation is useful.
4. This part is about the conventional velocity analysis method for the data from complex areas based on the time-offset hyperbolic curve of seismic reflection wave travel. Because of the approximation of hyperbolic hypothesis, the conventional velocity analysis cannot meet the requirement for establishing subsurface accurate velocity model with the data from complex areas. This paper puts forward the high-density and point by point dual spectral velocity analysis technique on the basis of previous research results. By introducing principles of high-density and point by point dual spectral velocity analysis, dynamic correction parameters to affect all the offsets, intelligent sorting of the pairν,η, geological filtering of high-density time differential correction parameter field (ν,η), unrelated dynamic correction parameters (dtn,τ0), interpolation and filtering ofνandηfields, this paper presents the theoretical formula of high-order non-hyperbolic time difference and puts forward the processing procedures used for high-density and point by point velocity analysis. This method uses the bi-directional scanning of dual-time parameters to obtain accurate velocity parameters. The actual processing of data acquired from Fangzheng Fault in Daqing Peripheral Basin demonstrates that this method can effectively solve the problem that the large-offset non-hyperbolic time difference correction gather couldn’t be leveled. Meanwhile, the point by point velocity analysis and calculation effectively solve the problem of lateral change of velocity and thus the accurate velocity fields of underground complex structures can be obtained.
5. In order to ensure the accuracy of migration imaging in complex areas, this paper analyzes the main problems faced by migration imaging in complex areas and puts forward the pre-stack time migration processing technique based on the anisotropy of surface relief. On the basis of imaging principles of surface common-reflection point, this paper uses the equivalent offset method for the data from the terrain with surface relief to eliminate the effect of surface relief on the structural morphologies and then formulates in detail the anisotropic pre-stack time migration processing technique. With the introduction of important concepts of anisotropy, phase velocity and ray velocity, NMO velocity in anisotropic medium, three dimensional description of NMO velocity in the anisotropic medium and principles of anisotropic pre-stack time migration method, this paper focuses on the anisotropic pre-stack time migration velocity analysis, high-accuracy calculation of velocity field, calculation of anisotropic parameters in VTI medium and the implementation process of anisotropic pre-stack time migration, and finally gets high-accuracy velocity modeling method and technical process for three dimensional anisotropic pre-stack time migration,ηfield method and process for obtaining anisotropic parameters in VTI medium with surface seismic data, three dimensional anisotropic pre-stack time migration processing method and procedures. It is proved by actual example that the surface–based anisotropic time migration processing technique can eliminate the effects of anisotropy on seismic wave travel path and velocity as well as the effects of surface relief on seismic data imaging and thus improves the accuracy of migration imaging. At the same time, because of the elimination of every anisotropic effect, it also solves the problem that the geophone offset on the CRP gather is over-corrected and in the meanwhile increases the effective multiplicity used for pre-stack inversion gather. And thus, the S/N ratio of migration result can be improved and it is also favorable for AVO analysis and pre-stack elastic impedance inversion.
At the end of this paper, the five techniques used for processing seismic data from complex areas, obtained from this research, are systematically summarized and concluded and a set of technological processes, suitable for processing low S/N ratio data from complex surfaces and complex structures, is given. And then, the actual data of 100 km2 from Fangzheng Fault in Daqing Peripheral Basin was processed with the above techniques. By comparing the overall effect of the final results, the new processed results with the previous 2D and 3D results, analyzing quantitatively the S/N ratios of new and old data and comparing the new interpreted results with the old one, it is demonstrated that not only the S/N ration but also the resolution of data are improved significantly, especially the structural morphologies have very big changes. After comparing with the exploration well data, the new-processed results match perfectly with the well data, which proves that the processing methods in this paper are advanced and practical. This paper also provides a set of new techniques for processing the low S/N ratio data from complex areas, which can be very reliable technical support for accurate structure interpretation, inversion prediction of reservoirs and for speeding up the seismic exploration of complex surfaces and complex structures.
引文
1.吴春红,张智勇,胡斌.低信噪比地震资料处理方法[J].中国科技信息, 2005,(10).
2.熊翥.复杂地区地震数据处理思路.北京:石油工业出版社,2002.
3.阎世信等.山地地球物理勘探技术.北京:石油工业出版社,2000.
4.熊翥.我国物探技术的进步与展望[J].石油地球物理勘探,2003,38(4):447-459.
5.张永刚.油气地球物理技术新进展—第二季度3届SEG年会论文概要.石油工业出版社,2004.
6. SEG 71、72届年会论文摘要,2003.
7.李世雄,刘家琦.小波变换和反演数学基础.地质出版社,2005.
8. SEG 73届年会论文摘要[C],2004.
9. SEG 77届年会论文摘要[C],2008.
10. SEG 76届年会论文摘要[C],2007.
11.牛滨华.波动方程压制多次波的技术方法.地学前缘,2002,9(2).
12. Perroud H , Hubral P , Hocht G,et al. Migrating around in circles - Part III. The Leading Edge ,1997 ,16(8):875-883.
13. Hubral P, Hoecht G, Jaeger, R. Seismic illumination. The Leading Edge, 1999,(11):1268-1271.
14. Schleicher J, Tygle M, Hubral P. Parabolic and hyperbolic paraxial two-point traveltimes in 3D media. Geophys. Prosp. , 1993,41:495-514.
15. Hubral P. Computing true amplitude reflections in a laterally inhomogeneous earth. Geophysics,1983,48:1051-1062.
16. SEG 75届年会论文摘要[C],2006.
17. SEG 74届年会论文摘要[C],2005.
18.熊翥.地震数据处理应用技术.石油工业出版社,2008.
19.杨勤勇.油气地球物理技术发展新动向.勘探地球物理进展,2007,30(2).
20.王喜双.中油股份公司物探技术现状及发展趋势.中国石油勘探,2006.4.12.
21.熊翥.21世纪初中期油气地球物理技术展望.石油工业出版社,2006.
22.吴如山.地震波散射:理论与应用[J].地球物理学进展,1989,4(1):1-23.
23. Wu R S,Aki K.弹性非均匀体引起的弹性波的散射特征[J].胡中平译.石油物探译丛,1987,(2):43-58.
24.吴如山,安艺敬一.地震波的散射与衰减[M].李裕澈,卢寿德译.北京:地震出版社,1993.1-50.
25.李灿苹,刘学伟,王祥春,杨丽.地震波的散射理论和散射特征及其应用[J].勘探地球物理进展,2005:81-88.
26.徐基祥,杨旭明,胡英,张才.一种在地震数据处理中衰减瑞雷波散射噪声的方法,申请号,(200810110902.5).
27.杨旭明.近地表地震波散射的研究与应用.成都理工大学博士学位论文,2002.
28.张文坡,陈湛文.线性化空间变换叠前与叠后地震记录去噪方法[J].地球物理学进展,2006,21(1).
29.陈湛文,刘勇.基于空间变换的F-X地震去噪新方法[J].大庆石油地质发,2002,21(3).
30.俞寿朋,蔡希玲,苏永昌.用地震信号多项式拟合提高叠加剖面信噪比[J].石油地球物理勘探,1988.
31. Hampson D. Inverse velocity stacking and multiple elimination. Journal of the Canadian Society of Exploration Geophysicists,1986,22,44-55.
32. Foster J. and Mosher C. Suppression of multiple reflections using the Radon transform. Geophysics, 1992,57,386-395.
33. Sacchi M.D.and Ulrych T.J. Improving resolution of Radon operators using a model reweighed least-squares procedure. Journal of Seismic Exploration,1995,4,315-328.
34. Zhou B. and Greenhalgh S.A.. Multiple suppression by 2D filtering in the parabolicτ?p domain: a wave-equation-based method. Geophysical prospecting, 1996,44,375-401.
35. Kabir M.M.N. and Marfurt K.J. Toward true amplitude multiple removal. The Leading Edge,1999,18,66-73.
36. Landa E., Belfer I. and Keydar S. Multiple attenuation in the parabolicτ?p domain using wavefront characteristics of multiple generating primaries. Geophysics, 1999,64,1806-1815.
37.沈操,牛滨华,余钦范.拉东变换的MATLAB实现.物探化探计算技术,2000,22(4):346-350.
38.朱生旺等.用抛物线Radon变换稀疏解分离和压制多次波[J].石油地球物理勘探,2002,37(2):110-115.
39. Wang Y.. Multiple attenuation: coping with the spatial truncation effect in the Radon transform domain. Geophysical prospecting,2003,51,75-87.
40.张军华,吕宁,田连玉.地震资料去噪方法、技术综合评述[J].地球物理学进展,2005,20(4):1083-1091.
41.张军华等.抛物线拉东变换消除多次波的应用要素分析[J].石油地球物理勘探,2004,39(4):398-405.
42.张军华等.频率域PRT方法及其在多次波去噪中的应用.六大勘探领域会议,中国青岛,2005年9月.
43. Marfurt K.J., Schneider R.V. and Mueller M.C.. Pitfalls of using conventional and discrete Radon transforms on poorly sampled data. Geophysics, 1996,61,1467-1482.
44.牛滨华等.多项式Radon变换[J].地球物理学报,2001,44(2):263-271.
45. Turner G. Aliasing in the tau-p transform and the remove of spatial aliased coherent noise. Geophysics,1990,55,1496-1503.
46.徐士良.Fortran常用算法程序集.清华大学出版社,1992.
47.李振春,张军华.地震资料数字处理讲义.石油大学出版社,2003.
48.张军华等.三江盆地低信噪比资料的多次波处理[J].石油地球物理勘探,2006,41(5):573-579.
49. O’Doherty R FandAnstay N A.Reflection on amplitudes.eophysics Prosp. 1971,430-439.
50. Sharpiro S A and peter Hubral, Elastic waves in fineiy latered sediments medium and generlinsed O’Doherty-Anstay formulas. Geophysics. 1996, 61(5):1282-1300.
51. Martin.T.widmeir, Sargei.A ,Peter Hubral , AVO correction for scalar waves in the case a thinly layered reflector overburden. Geophysics,1996,61(2):520-528.
52.徐仲达.透射损失校正在AVO处理中的应用[J].石油物探,2002,1(4):406-409.
53.王兆湖等.叠后地震信号透射损失的补偿[J].石油地球物理勘探,2008,43(2):217-221.
54. Tsvankin L , Thomsen L Nonhyperbolic reflection moveout inanisotropic media[J].Geophys -ics,1994,59(8):1290-1304.
55. Alkhalifah T,Tsvankin I.Velocityanalysis for transversely isotropic media[J].Geophysics , 1995 ,60 (5) :1550– 1566.
56. Siliqi R ,D. Le Meur , Gamar F ,et al. High2density moveout parameter fields V andη[C] . 73rd SEG International Exposition &Annual Meeting - Dallas ,Texas ,26 - 31 October ,2003.
57. Siliqi,R. and Bousquie,N., Anelliptic time processing based on a shifted hyperbola approach, 2000,70 th Ann.Internat.Mtg:Soc of Expl.Geophys.,2245-2248
58. Castle,R.J., Atheory of normal moveout :Geophysics, 1994,59.983-999.
59. Siliqi,R., Technological leap in time processing focuses the data throughout anisotropic media:First Break,2001,19,no.11,612-618.
60. Le Meur,D., Siliqi,R., High-density moveout parameter field V and ?. Part two:Simultaneous geostatistical filtering:2003,73rd Ann.Internat.Mtg., Soc of Expl.Geophys
61. De Bazelaire,E. and Viallix,J.R., Normal moveout in focus:Geophys.Prosp.,Eur. Assn.Geosci.Eng.,1994,42,477-618.
62. Adler,F.,and Brandwood,S., Robust estimation of dense 3D stacking velocities from automated picking:69th Ann.Internat.Mtg.,SEG Expanded Abstracts.1999.
63. Alkhalifah,T.,Velocity analysis using nonhyperbolic moveout in transversely isotropic media:Geophysics,1997,62,1839-1854.
64. Anderson,J.E.,and Tsvankin,I.,Dip-moveout processing by Fourier transform in anisotropic media[J].Geophysics,1997,62:1260-1269 .
65. Siliqi,R., Le Meur,D., High-density moveout parameter field V and ?. Part one:Simultaneous automatic picking:2003,73rd Ann.Internat.Mtg., SEG Expanded Abstracts.
66. Le Meur,D.and Hermann,F., Pre-conditioning of the densely sampled stacking velocity field., 2002, 64th EAGE Annual Meeting.
67. Matheron G., Principles of geostatistics. 1963,Economic Geol.58,1246-1290.
68.马在田.地震成像技术.北京:石油工业出版社,1898.
69. Z.dongetal.各向异性介质中的三维叠前偏移. Geophysics, 1993.
70.刘彦强,董敏煜.对各向同性到各向异性带来的几个问题的初探.石油地球物理勘探,1992,27(1):29-44.
71. Tariq Alkhalifah等.横向各向异性介质的速度分析. Geophysics,1995.
72. Tariq Alkhalifah.横向各向异性介质中应用非双曲线时距曲线进行速度分析.GRI通讯, 1996增刊.
73. Thomsen,L.,Weak elastic anisotropy [J].Geophysics,1986,51:1954-1966
74.刘怀山等译.地震资料分析-地震资料处理、反演和解释[M].北京:石油工业出版社,2006.
75. Berryman,J.G.,Long-wave elastic anisotropy in transversely isotropic media[J].Geophysics, 1979,44:896-917.
76. Tsvankin,I.,Seismic Signatures and Analysis of Reflection Data in Anisotropic Media,Handbook of Geophysical Exploration seismic Exploration [M].Volume 29.
77. IlyaTsvankin.方位角各向异性介质中群速度和动校正速度的解析表达. GRI通讯,1996增刊.
78.张翠兰.地层各向异性研究综述.山地地震勘探,1995.
79.陈宝书等.地震资料叠前偏移处理技术应用.中国海上油气(地质),2001,15(5):361-364.
80.张中杰.地震各向异性研究进展[J].地球物理学进展, 2002.
81.马在田.高阶有限差分偏移.石油地球物理勘探,1982,17:6-15.
82.徐亦鸣,黄中玉,刘路佳.各向异性介质纵波速度分析[J].石油物探,2004,43(5)438-528.
83.杜启振,李辉.各向异性介质非双曲时差速度分析.油气地球物理,2005.
2.熊翥.复杂地区地震数据处理思路.北京:石油工业出版社,2002.
3.阎世信等.山地地球物理勘探技术.北京:石油工业出版社,2000.
4.熊翥.我国物探技术的进步与展望[J].石油地球物理勘探,2003,38(4):447-459.
5.张永刚.油气地球物理技术新进展—第二季度3届SEG年会论文概要.石油工业出版社,2004.
6. SEG 71、72届年会论文摘要,2003.
7.李世雄,刘家琦.小波变换和反演数学基础.地质出版社,2005.
8. SEG 73届年会论文摘要[C],2004.
9. SEG 77届年会论文摘要[C],2008.
10. SEG 76届年会论文摘要[C],2007.
11.牛滨华.波动方程压制多次波的技术方法.地学前缘,2002,9(2).
12. Perroud H , Hubral P , Hocht G,et al. Migrating around in circles - Part III. The Leading Edge ,1997 ,16(8):875-883.
13. Hubral P, Hoecht G, Jaeger, R. Seismic illumination. The Leading Edge, 1999,(11):1268-1271.
14. Schleicher J, Tygle M, Hubral P. Parabolic and hyperbolic paraxial two-point traveltimes in 3D media. Geophys. Prosp. , 1993,41:495-514.
15. Hubral P. Computing true amplitude reflections in a laterally inhomogeneous earth. Geophysics,1983,48:1051-1062.
16. SEG 75届年会论文摘要[C],2006.
17. SEG 74届年会论文摘要[C],2005.
18.熊翥.地震数据处理应用技术.石油工业出版社,2008.
19.杨勤勇.油气地球物理技术发展新动向.勘探地球物理进展,2007,30(2).
20.王喜双.中油股份公司物探技术现状及发展趋势.中国石油勘探,2006.4.12.
21.熊翥.21世纪初中期油气地球物理技术展望.石油工业出版社,2006.
22.吴如山.地震波散射:理论与应用[J].地球物理学进展,1989,4(1):1-23.
23. Wu R S,Aki K.弹性非均匀体引起的弹性波的散射特征[J].胡中平译.石油物探译丛,1987,(2):43-58.
24.吴如山,安艺敬一.地震波的散射与衰减[M].李裕澈,卢寿德译.北京:地震出版社,1993.1-50.
25.李灿苹,刘学伟,王祥春,杨丽.地震波的散射理论和散射特征及其应用[J].勘探地球物理进展,2005:81-88.
26.徐基祥,杨旭明,胡英,张才.一种在地震数据处理中衰减瑞雷波散射噪声的方法,申请号,(200810110902.5).
27.杨旭明.近地表地震波散射的研究与应用.成都理工大学博士学位论文,2002.
28.张文坡,陈湛文.线性化空间变换叠前与叠后地震记录去噪方法[J].地球物理学进展,2006,21(1).
29.陈湛文,刘勇.基于空间变换的F-X地震去噪新方法[J].大庆石油地质发,2002,21(3).
30.俞寿朋,蔡希玲,苏永昌.用地震信号多项式拟合提高叠加剖面信噪比[J].石油地球物理勘探,1988.
31. Hampson D. Inverse velocity stacking and multiple elimination. Journal of the Canadian Society of Exploration Geophysicists,1986,22,44-55.
32. Foster J. and Mosher C. Suppression of multiple reflections using the Radon transform. Geophysics, 1992,57,386-395.
33. Sacchi M.D.and Ulrych T.J. Improving resolution of Radon operators using a model reweighed least-squares procedure. Journal of Seismic Exploration,1995,4,315-328.
34. Zhou B. and Greenhalgh S.A.. Multiple suppression by 2D filtering in the parabolicτ?p domain: a wave-equation-based method. Geophysical prospecting, 1996,44,375-401.
35. Kabir M.M.N. and Marfurt K.J. Toward true amplitude multiple removal. The Leading Edge,1999,18,66-73.
36. Landa E., Belfer I. and Keydar S. Multiple attenuation in the parabolicτ?p domain using wavefront characteristics of multiple generating primaries. Geophysics, 1999,64,1806-1815.
37.沈操,牛滨华,余钦范.拉东变换的MATLAB实现.物探化探计算技术,2000,22(4):346-350.
38.朱生旺等.用抛物线Radon变换稀疏解分离和压制多次波[J].石油地球物理勘探,2002,37(2):110-115.
39. Wang Y.. Multiple attenuation: coping with the spatial truncation effect in the Radon transform domain. Geophysical prospecting,2003,51,75-87.
40.张军华,吕宁,田连玉.地震资料去噪方法、技术综合评述[J].地球物理学进展,2005,20(4):1083-1091.
41.张军华等.抛物线拉东变换消除多次波的应用要素分析[J].石油地球物理勘探,2004,39(4):398-405.
42.张军华等.频率域PRT方法及其在多次波去噪中的应用.六大勘探领域会议,中国青岛,2005年9月.
43. Marfurt K.J., Schneider R.V. and Mueller M.C.. Pitfalls of using conventional and discrete Radon transforms on poorly sampled data. Geophysics, 1996,61,1467-1482.
44.牛滨华等.多项式Radon变换[J].地球物理学报,2001,44(2):263-271.
45. Turner G. Aliasing in the tau-p transform and the remove of spatial aliased coherent noise. Geophysics,1990,55,1496-1503.
46.徐士良.Fortran常用算法程序集.清华大学出版社,1992.
47.李振春,张军华.地震资料数字处理讲义.石油大学出版社,2003.
48.张军华等.三江盆地低信噪比资料的多次波处理[J].石油地球物理勘探,2006,41(5):573-579.
49. O’Doherty R FandAnstay N A.Reflection on amplitudes.eophysics Prosp. 1971,430-439.
50. Sharpiro S A and peter Hubral, Elastic waves in fineiy latered sediments medium and generlinsed O’Doherty-Anstay formulas. Geophysics. 1996, 61(5):1282-1300.
51. Martin.T.widmeir, Sargei.A ,Peter Hubral , AVO correction for scalar waves in the case a thinly layered reflector overburden. Geophysics,1996,61(2):520-528.
52.徐仲达.透射损失校正在AVO处理中的应用[J].石油物探,2002,1(4):406-409.
53.王兆湖等.叠后地震信号透射损失的补偿[J].石油地球物理勘探,2008,43(2):217-221.
54. Tsvankin L , Thomsen L Nonhyperbolic reflection moveout inanisotropic media[J].Geophys -ics,1994,59(8):1290-1304.
55. Alkhalifah T,Tsvankin I.Velocityanalysis for transversely isotropic media[J].Geophysics , 1995 ,60 (5) :1550– 1566.
56. Siliqi R ,D. Le Meur , Gamar F ,et al. High2density moveout parameter fields V andη[C] . 73rd SEG International Exposition &Annual Meeting - Dallas ,Texas ,26 - 31 October ,2003.
57. Siliqi,R. and Bousquie,N., Anelliptic time processing based on a shifted hyperbola approach, 2000,70 th Ann.Internat.Mtg:Soc of Expl.Geophys.,2245-2248
58. Castle,R.J., Atheory of normal moveout :Geophysics, 1994,59.983-999.
59. Siliqi,R., Technological leap in time processing focuses the data throughout anisotropic media:First Break,2001,19,no.11,612-618.
60. Le Meur,D., Siliqi,R., High-density moveout parameter field V and ?. Part two:Simultaneous geostatistical filtering:2003,73rd Ann.Internat.Mtg., Soc of Expl.Geophys
61. De Bazelaire,E. and Viallix,J.R., Normal moveout in focus:Geophys.Prosp.,Eur. Assn.Geosci.Eng.,1994,42,477-618.
62. Adler,F.,and Brandwood,S., Robust estimation of dense 3D stacking velocities from automated picking:69th Ann.Internat.Mtg.,SEG Expanded Abstracts.1999.
63. Alkhalifah,T.,Velocity analysis using nonhyperbolic moveout in transversely isotropic media:Geophysics,1997,62,1839-1854.
64. Anderson,J.E.,and Tsvankin,I.,Dip-moveout processing by Fourier transform in anisotropic media[J].Geophysics,1997,62:1260-1269 .
65. Siliqi,R., Le Meur,D., High-density moveout parameter field V and ?. Part one:Simultaneous automatic picking:2003,73rd Ann.Internat.Mtg., SEG Expanded Abstracts.
66. Le Meur,D.and Hermann,F., Pre-conditioning of the densely sampled stacking velocity field., 2002, 64th EAGE Annual Meeting.
67. Matheron G., Principles of geostatistics. 1963,Economic Geol.58,1246-1290.
68.马在田.地震成像技术.北京:石油工业出版社,1898.
69. Z.dongetal.各向异性介质中的三维叠前偏移. Geophysics, 1993.
70.刘彦强,董敏煜.对各向同性到各向异性带来的几个问题的初探.石油地球物理勘探,1992,27(1):29-44.
71. Tariq Alkhalifah等.横向各向异性介质的速度分析. Geophysics,1995.
72. Tariq Alkhalifah.横向各向异性介质中应用非双曲线时距曲线进行速度分析.GRI通讯, 1996增刊.
73. Thomsen,L.,Weak elastic anisotropy [J].Geophysics,1986,51:1954-1966
74.刘怀山等译.地震资料分析-地震资料处理、反演和解释[M].北京:石油工业出版社,2006.
75. Berryman,J.G.,Long-wave elastic anisotropy in transversely isotropic media[J].Geophysics, 1979,44:896-917.
76. Tsvankin,I.,Seismic Signatures and Analysis of Reflection Data in Anisotropic Media,Handbook of Geophysical Exploration seismic Exploration [M].Volume 29.
77. IlyaTsvankin.方位角各向异性介质中群速度和动校正速度的解析表达. GRI通讯,1996增刊.
78.张翠兰.地层各向异性研究综述.山地地震勘探,1995.
79.陈宝书等.地震资料叠前偏移处理技术应用.中国海上油气(地质),2001,15(5):361-364.
80.张中杰.地震各向异性研究进展[J].地球物理学进展, 2002.
81.马在田.高阶有限差分偏移.石油地球物理勘探,1982,17:6-15.
82.徐亦鸣,黄中玉,刘路佳.各向异性介质纵波速度分析[J].石油物探,2004,43(5)438-528.
83.杜启振,李辉.各向异性介质非双曲时差速度分析.油气地球物理,2005.
© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号 地址：北京市海淀区学院路29号 邮编：100083 电话：办公室：(+86 10)66554848；文献借阅、咨询服务、科技查新：66554700 |