复杂地质条件地震多波照明及地震采集方法研究
|
摘要
随着国民经济的快速发展,国内对能源的需求也越来越紧张。为了加大对能源的储备,地震勘探的重心从东部平原地区向中西部山地、高原地区转移。在这些地区,蕴藏着丰富的石油天然气,相继发现了克拉玛依油田、塔里木油田、吐哈油田等一些大的能源基地。但是在这些地区,沙漠、山地、高原等复杂的地形构成了不同于东部地区的复杂地质结构,因此,基于共中心点(CMP)的野外采集方式已经不适合复杂的地质条件。对于新的探区,野外地震勘探投资大,风险高,为了提高经济效益,迫切需要一种技术手段能够对新探区的采集参数进行模拟和研究分析,通过参数论证和应用研究,优化野外采集的观测系统,提高对勘探目标的勘探精度,为寻找新的隐伏油气藏提供理论指导和技术支持。
限差分算子,并且进一步研究了其稳定性、收敛性和频散关系。在地震勘探中,人工激发的地震波在岩石中传播,地震波遇到介质会发生反射、透射等波的运动形式,通过这种介质的机械振动研究岩石的物理特性。岩在野外地震勘探初期,地震采集方式可以通过模拟地质构造进行参数论证,对地震覆盖次数、偏移距、道距等采集参数进行室内模拟,因此,通过模拟获得的结论进行勘探目标的验证,进而修改采集参数,提高勘探效益。常规的建模技术是建立在水平层状介质的基础上,这种建模方式简单,对我国东部平原地区比较适合,但是对中西部高原地区由于复杂的地质结构进行模拟变得非常困难。本文在研究常规建模的基础上引入了点、断、面的二维建模方法,通过模拟西部复杂的逆掩断层等地质结构取得了较好效果。
波动方程正演方法不同于野外地震勘探,由于计算机硬件的限制,因此造成了人为的边界问题,同时由于是室内模拟也产生了对震源问题的研究。为了对波的传播这种连续问题离散化,国内外学者都在研究过程中获得了许多成果。在边界问题的研究中先后推导出了像完全吸收边界等一些效果明显的边界解决方法。通过对野外震源特性的研究分析,先后研究出了许多子波函数表达式,结合地震勘探的特点,其中一种能够模拟野外震源的雷克子波在波场模拟中得到了广泛运用,这种子波具有分辨率高、频带宽的特点。
在研究地震波传播的波场特征时,正演模拟是一种不可缺少的手段。地震波数值模拟主要分为几何射线法和波动方程法两大类。几何射线法是在高频近似的前提下模拟波的传播过程的,它主要反映波的运动学特征,而波动方程法在模拟波的传播过程时不但能提供波的运动学特征还能反映波在介质中的动力学参数,因此,波动方程法正演模拟得到了更广泛的应用。随着计算机软件、硬件的飞速发展,数值模拟技术在各个领域得到了迅速应用。在数值模拟技术发展中研究出了许多波动方程模拟方法,主要有有限元法、伪谱法和有限差分法,其中有限差分法编程简单,能够模拟各种地质结构,在地震波数值模拟中得到了广泛应用。由于对连续问题离散化,从简单的声波方程出发推导出了波动方程的有石圈中普遍存在各向异性介质即横向各向同性介质,从广义的胡克定律出发,根据描述弹性介质运动的微分方程建立了有限差分公式,以有限差分算子为基础进行地震照明模拟,通过对纵波和横波(本文主要指转换横波SV波)的模拟照明,得出了有意义的结论,对提高野外勘探效益有一定的指导意义。
地震照明能够模拟野外地震勘探采集方式,优化采集观测系统的设计,降低投资风险,提高经济效益,因此,经常作为一种评价野外采集的方法手段。
地震勘探采集参数(激发震源的药量、排列的接收方式、道距、偏移距等)以及地下地质构造是影响地震照明的因素。
最先用于对地震照明研究的方法是射线追踪,这种方法在基于水平层状结构的地质模型中能够很好的反映地震波在介质中的照明度。随着地震勘探的目标区越来越复杂,由于射线追踪法本身的缺陷已经不能适应复杂的地质模型,经常在模拟过程中得出不合理的结论。由于波动方程能够适应复杂的地质条件,能够提供更多的波场信息,有助于研究介质的物性参数。通过建立复杂的地质模型,模拟地下隐伏的含油气蕴藏环境,对不同的地质模型通过变换多种采集参数得到不同的照明结果进行分析,可以对勘探目标的采集方式进行优化,提高勘探精度,使在基于共反射点的观测系统的设计上进一步提高采集参数的优化,真正实现基于勘探目标的野外采集方式,为提高经济效益提供理论指导和技术支持。
地震资料去噪技术是地震勘探过程的一个重要环节,如何把采集的资料经过处理达到真正反映地下地质结构的成果,去噪技术发挥了重要作用,特别是在复杂的地质条件影响下,解决长波长、短波长校正量问题也是一个重要研究课题。
通过对复杂地质条件下的地质构造进行地震照明和地震资料处理,得到了较好的成果,对进一步研究复杂地质构造积累了经验。
Domestic demand for energy is more and more tensely along with the rapid development of the national economy.In order to increase the energy reserves, the strategic center of the sei-smic exploration has been moved from the eastern plain areas to the Midwest mountain and pl-ateau regions. Some big oilfields have been successively found in these areas, rich in oil and gas. such as Kelamayi oilfield、Talim oilfield and Tuha oil field. But in these areas, it is different from the eastern regions with the complex geological structure. Therefore, the acquisition methods based on CMP are not suitable for complicated geological conditions. Seismic exploration has large investment risk for the new exploraion areas.In order to improve the economic benefit,some new exploration techologies will be used to study and analysis the collection parameters in desperate need, which will optimize the field collection observation system and improve the accuracy of exploration targets and provide theoretical guidance and technical support for searching the new concealed reservoirs.
In the field seismic exploration early.it can be simulated for the geological structure parameters of argumentation,such as seismic folds、offset and recevers indoor simulation.Therefore, it can modify the gathering parameters and increase the exploration benefit by the conclusion obtained in exploration simulation.It is simple to establish the layered media on the basis of the conventional modeling technology and it is suitable for the eastern plain area but not for the Midwest plateau due to the complicated geological structure.In this thesis, it will be introduced a new two-dimensional modeling method by studying the conventional modeling.its basic elements are dot、broken and surface. The new modeling method is successful in the sim-ulation of the overthrust fault in the western complex geological structure.
Because the computer hardware restrictions and the man-made boundary issue.the wave equation modeling algorithm is different from the wild seismic exploration. At the same time,we must study the source in the indoor simulation. In order to discrete the continuous problem, scholars at home and abroad have received many achievements in the process of studying the wave propagation. There are some good effect in the studying boundary problem.such as fully absorting boundary. Many wavelet function have been successively developed by the research and analysis of the source characteristics. One of the famous wavelets is Ricker wavelet, and it has been widely used in the wave field sim-ulation with the high resolution wade frequency characteristics.
The forward modeling is a useful tool in the study of the wave field characteristic. There are two methods to simulate the seismic wave propagation.geometric ray method and wave equation method. Geometric ray method simulate the seismic wave propagation on the base of the high frequency approximation, it mainly reflects the wave of kinematic features, and wave equation method can not only provide the wave of kinematic features but also reflects the wave of the kinetic parameters. Therefore,the wave equation method has been applied extensiveely in the forward modeling. Along with the rapid development of the computer software and hardware, the numerical simulation technology has been used in many fields. In the wave equation simulation technology,there are three numerical simulation methods, including finite element method、the method of pseudo-spectral and finite difference method. Because of the simple and simulating various geological structure,the finite difference method has been used widely. In the process of studing the wave propagation with the discretization problem, the finite difference operators has been deduced, at the same time, its stability、convergence and frequency dispersion relationnship have been studied too. In the seismic exploration, the physical properties of the rocks are studied by the artificial mechanical vibration of the medium when the seismic wave propagats. The anisotropic medium namely transverse isotropic media exists widely in the lithosphere. The finite difference formula,describing the movement of the elastic medium, has been established from the broad Hooker law.and the seismic illumination simulation with finite difference operators began laterly.ln this text,it is concluded that it is very helpful to improve the field exploration by simulating the illumination of the compress waves and shear waves (it refers to convert sv waves here).
The seismic illumination can simulate the field seismic exploration collecting methods, optimize the designation of the system, so it can reduce the risk of investment and improve the economic benefit. Therefore, it are used as a kind of evaluation methods of field collection usually.
The seismic illumination are influenced by the seismic exploration gathering parameters.such as the source、the arrangement of receiving mode、the distance between the recevers and the subsurface geologic structure.
Firstly,the ray tracing was used in the seismic research method for the illumination.and this kind of method works well in the horizontal layer structure of geological model, but does bad in the complicated geological model due to the defects of ray tracing law itself. Because of adapting to the complicated geological model, the wave equation method can provide more information to study of the medium physical parameter in the wave field. By establishing the complicated geological model and underground concealed deposits of oil and gas with differrent geological environment, we can obtain the conclusion of the illumination, optimizing the acquisition method, increasing explorat-ion precision.So we can establish the optimizated acquisition system based on exploration targets truly,with which to improve the economic benefit with theoretical guidance and technical support.
The seismic denoising technology is important in the process of the seismic exploration, it plays an important role to obtain the true imaging underground geologic structure, especially in the complicated geological conditions,the influence of the long wavelength and the short wavelength statics problem is important in the process of the seismic data.
In this thesis,I have gotten the achievements in studing the seismic illumination and processing the seismic data under the complex geological conditions and the experience to study of the complex ge-ological region.
限差分算子,并且进一步研究了其稳定性、收敛性和频散关系。在地震勘探中,人工激发的地震波在岩石中传播,地震波遇到介质会发生反射、透射等波的运动形式,通过这种介质的机械振动研究岩石的物理特性。岩在野外地震勘探初期,地震采集方式可以通过模拟地质构造进行参数论证,对地震覆盖次数、偏移距、道距等采集参数进行室内模拟,因此,通过模拟获得的结论进行勘探目标的验证,进而修改采集参数,提高勘探效益。常规的建模技术是建立在水平层状介质的基础上,这种建模方式简单,对我国东部平原地区比较适合,但是对中西部高原地区由于复杂的地质结构进行模拟变得非常困难。本文在研究常规建模的基础上引入了点、断、面的二维建模方法,通过模拟西部复杂的逆掩断层等地质结构取得了较好效果。
波动方程正演方法不同于野外地震勘探,由于计算机硬件的限制,因此造成了人为的边界问题,同时由于是室内模拟也产生了对震源问题的研究。为了对波的传播这种连续问题离散化,国内外学者都在研究过程中获得了许多成果。在边界问题的研究中先后推导出了像完全吸收边界等一些效果明显的边界解决方法。通过对野外震源特性的研究分析,先后研究出了许多子波函数表达式,结合地震勘探的特点,其中一种能够模拟野外震源的雷克子波在波场模拟中得到了广泛运用,这种子波具有分辨率高、频带宽的特点。
在研究地震波传播的波场特征时,正演模拟是一种不可缺少的手段。地震波数值模拟主要分为几何射线法和波动方程法两大类。几何射线法是在高频近似的前提下模拟波的传播过程的,它主要反映波的运动学特征,而波动方程法在模拟波的传播过程时不但能提供波的运动学特征还能反映波在介质中的动力学参数,因此,波动方程法正演模拟得到了更广泛的应用。随着计算机软件、硬件的飞速发展,数值模拟技术在各个领域得到了迅速应用。在数值模拟技术发展中研究出了许多波动方程模拟方法,主要有有限元法、伪谱法和有限差分法,其中有限差分法编程简单,能够模拟各种地质结构,在地震波数值模拟中得到了广泛应用。由于对连续问题离散化,从简单的声波方程出发推导出了波动方程的有石圈中普遍存在各向异性介质即横向各向同性介质,从广义的胡克定律出发,根据描述弹性介质运动的微分方程建立了有限差分公式,以有限差分算子为基础进行地震照明模拟,通过对纵波和横波(本文主要指转换横波SV波)的模拟照明,得出了有意义的结论,对提高野外勘探效益有一定的指导意义。
地震照明能够模拟野外地震勘探采集方式,优化采集观测系统的设计,降低投资风险,提高经济效益,因此,经常作为一种评价野外采集的方法手段。
地震勘探采集参数(激发震源的药量、排列的接收方式、道距、偏移距等)以及地下地质构造是影响地震照明的因素。
最先用于对地震照明研究的方法是射线追踪,这种方法在基于水平层状结构的地质模型中能够很好的反映地震波在介质中的照明度。随着地震勘探的目标区越来越复杂,由于射线追踪法本身的缺陷已经不能适应复杂的地质模型,经常在模拟过程中得出不合理的结论。由于波动方程能够适应复杂的地质条件,能够提供更多的波场信息,有助于研究介质的物性参数。通过建立复杂的地质模型,模拟地下隐伏的含油气蕴藏环境,对不同的地质模型通过变换多种采集参数得到不同的照明结果进行分析,可以对勘探目标的采集方式进行优化,提高勘探精度,使在基于共反射点的观测系统的设计上进一步提高采集参数的优化,真正实现基于勘探目标的野外采集方式,为提高经济效益提供理论指导和技术支持。
地震资料去噪技术是地震勘探过程的一个重要环节,如何把采集的资料经过处理达到真正反映地下地质结构的成果,去噪技术发挥了重要作用,特别是在复杂的地质条件影响下,解决长波长、短波长校正量问题也是一个重要研究课题。
通过对复杂地质条件下的地质构造进行地震照明和地震资料处理,得到了较好的成果,对进一步研究复杂地质构造积累了经验。
Domestic demand for energy is more and more tensely along with the rapid development of the national economy.In order to increase the energy reserves, the strategic center of the sei-smic exploration has been moved from the eastern plain areas to the Midwest mountain and pl-ateau regions. Some big oilfields have been successively found in these areas, rich in oil and gas. such as Kelamayi oilfield、Talim oilfield and Tuha oil field. But in these areas, it is different from the eastern regions with the complex geological structure. Therefore, the acquisition methods based on CMP are not suitable for complicated geological conditions. Seismic exploration has large investment risk for the new exploraion areas.In order to improve the economic benefit,some new exploration techologies will be used to study and analysis the collection parameters in desperate need, which will optimize the field collection observation system and improve the accuracy of exploration targets and provide theoretical guidance and technical support for searching the new concealed reservoirs.
In the field seismic exploration early.it can be simulated for the geological structure parameters of argumentation,such as seismic folds、offset and recevers indoor simulation.Therefore, it can modify the gathering parameters and increase the exploration benefit by the conclusion obtained in exploration simulation.It is simple to establish the layered media on the basis of the conventional modeling technology and it is suitable for the eastern plain area but not for the Midwest plateau due to the complicated geological structure.In this thesis, it will be introduced a new two-dimensional modeling method by studying the conventional modeling.its basic elements are dot、broken and surface. The new modeling method is successful in the sim-ulation of the overthrust fault in the western complex geological structure.
Because the computer hardware restrictions and the man-made boundary issue.the wave equation modeling algorithm is different from the wild seismic exploration. At the same time,we must study the source in the indoor simulation. In order to discrete the continuous problem, scholars at home and abroad have received many achievements in the process of studying the wave propagation. There are some good effect in the studying boundary problem.such as fully absorting boundary. Many wavelet function have been successively developed by the research and analysis of the source characteristics. One of the famous wavelets is Ricker wavelet, and it has been widely used in the wave field sim-ulation with the high resolution wade frequency characteristics.
The forward modeling is a useful tool in the study of the wave field characteristic. There are two methods to simulate the seismic wave propagation.geometric ray method and wave equation method. Geometric ray method simulate the seismic wave propagation on the base of the high frequency approximation, it mainly reflects the wave of kinematic features, and wave equation method can not only provide the wave of kinematic features but also reflects the wave of the kinetic parameters. Therefore,the wave equation method has been applied extensiveely in the forward modeling. Along with the rapid development of the computer software and hardware, the numerical simulation technology has been used in many fields. In the wave equation simulation technology,there are three numerical simulation methods, including finite element method、the method of pseudo-spectral and finite difference method. Because of the simple and simulating various geological structure,the finite difference method has been used widely. In the process of studing the wave propagation with the discretization problem, the finite difference operators has been deduced, at the same time, its stability、convergence and frequency dispersion relationnship have been studied too. In the seismic exploration, the physical properties of the rocks are studied by the artificial mechanical vibration of the medium when the seismic wave propagats. The anisotropic medium namely transverse isotropic media exists widely in the lithosphere. The finite difference formula,describing the movement of the elastic medium, has been established from the broad Hooker law.and the seismic illumination simulation with finite difference operators began laterly.ln this text,it is concluded that it is very helpful to improve the field exploration by simulating the illumination of the compress waves and shear waves (it refers to convert sv waves here).
The seismic illumination can simulate the field seismic exploration collecting methods, optimize the designation of the system, so it can reduce the risk of investment and improve the economic benefit. Therefore, it are used as a kind of evaluation methods of field collection usually.
The seismic illumination are influenced by the seismic exploration gathering parameters.such as the source、the arrangement of receiving mode、the distance between the recevers and the subsurface geologic structure.
Firstly,the ray tracing was used in the seismic research method for the illumination.and this kind of method works well in the horizontal layer structure of geological model, but does bad in the complicated geological model due to the defects of ray tracing law itself. Because of adapting to the complicated geological model, the wave equation method can provide more information to study of the medium physical parameter in the wave field. By establishing the complicated geological model and underground concealed deposits of oil and gas with differrent geological environment, we can obtain the conclusion of the illumination, optimizing the acquisition method, increasing explorat-ion precision.So we can establish the optimizated acquisition system based on exploration targets truly,with which to improve the economic benefit with theoretical guidance and technical support.
The seismic denoising technology is important in the process of the seismic exploration, it plays an important role to obtain the true imaging underground geologic structure, especially in the complicated geological conditions,the influence of the long wavelength and the short wavelength statics problem is important in the process of the seismic data.
In this thesis,I have gotten the achievements in studing the seismic illumination and processing the seismic data under the complex geological conditions and the experience to study of the complex ge-ological region.
引文
[1]Alford.R.M,Kelly.K.R,and Boore.D.M,1974,Accuracy of finite-difference modeling of the acoustic wave equation:Geophysics,39:834-842
[2]Alterman Z,Karal F C. Propagation of seismic wave in layered media by finite difference methods[J].BSSA,1968,58(1):367-398
[3]Aminzadeh F.,Brac J.,Kunz T.3-D Salt and Overthrust Models, SEG/EAGE 3-D Modeling Series No.1.Soc.Explor.Geophysicists,Tulsa,1997
[4]Bear G,Lu C,Lu R.The construction of subsurface illumination and amplitude maps via ray tracing.The leading Edge,2000,19(7):726-731
[5]Berenger J P.A perfectly matched layer for the absorption of elect romagnetic waves.Journal of Computational Physics,1994,114(2):185-200
[6]Berenger J P.Three-dimensional perfectly matched layer for the absorption of electromagnetic waves.Journal of Computational Physics,1996,127(2):363-379
[7]Boor.D.M.1970.Finite-difference solution to the equations of elastic wave propagation with application to Love waves over dipping interfaces:Ph.D thesis.M.I.T
[8]B Fornberg.The pseudospectral method:Acurate representation of interfaces in elastic wave calculations. Geophysics.1988,53(5):625-637
[9]Chen K H.Seismic modeling:A finite element approach[J].SEG/China Geophysical Society Metting,Beijing,China,1981.
[10]Clayton R.Engquist B.Absorbing boundary conditions for acoustic and elastic wave equations [J].The Bulletin of the Seismological Society of America,1977,67:1529-1540
[11]Clough R W.The finite elementmethod in plane stress analysis[A].1960, In:Proceedings of 2nd ASCE Conference on Electronic Computation Pittsburgh, ASCE.
[12]Collino F,Tsogka C. Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media.Geo-physics,2001,66(1):294-307
[13]David Muerdter,Diamond Geoscience Research Corporation,Redmond, Washington,U.S.DAVIS RATCLIFF,Diamond Geophysical Service Corporation Houston.Texas.U.S.Understanding subsalt illumination through ray-trace modeling.Part 1:Simple 2-D salt models. The Leading Edge.2001,1(1):578-586
[14]David Muerdter等著,娄晓东译,孙夕平,赵该善校.通过射线追踪模拟理解盐下照射[M].石油物探译丛,2001,10(5)：59-70
[15]Falk J, Tessmer E,Gajewski D. Efficient finite-difference modeling of seismic waves using locally adjustable time steps [J]. Geophysical Prospecting,1998,46:603-616.
[16]Gelius L.J.Lecomte I.and Tabti H.Analysis of the resolution function in seismic prestack depth imaging.Geophys.Prosp.2002.50:505-515
[17]Greg J.Schurter.Yayat Supriatna,Angke Nuraeni.and Supriyono.Bp In donesia Carl J.Regone and Nurul Kabir.Bp EPTG.Houston.A 3D finitedifference modeling study of seismic imaging challenges in Bintuni Bay,Irian Jaya Barat.2009,1008-1018
[18]Gustavo Catao Alves*,Andre Bulcao, Djalma Manoel Soares Filho, Carlos Eduardo Theodoro,Luiz Alberto Santos.Marcos Antonio Gallotti Guima raes,Petrobras.Target oriented illumination analysis using wave equation [J]. SEG Las Vegas Annual Meeting.2008,163-167
[19]Hastings F.Schneider J B.Broschat S L.Application of the perfectly matchedlayer(PML) absorbing boundary condition to elastic wave pro pagation[J] Journal of the Acoustic Society of America,1996,100(5):3061-3069
[20]Hedstrom G W.Nonreflecting boundary conditions for nonlinear hype rbolic systerms.J Comput Phys.1979.30:222-237
[21]Higdon R L. Absorbing boundary conditions for elastic wave[J]. Geophysics.1991.56(2):231-241
[22]Hoffmann J.Illumination.resolution and image quality of PP- and PS- waves for survey planning.The leading Edge.2001.20(9):1008-1012
[23]Igel H, Weber M. P_SV wave propagation in the Earth s mantle using finite differences:Application to heterogeneous lower most mantle structure [J].G eophys,Res.Lett,1966,23:415-418
[24]Jastram C,Tessemer E.Elastic modeling on a grid with vertically varying spacing[J].Geophysical prospecting,1994,42:357-370
[25]K.R.KELLY*,R.W.WARD.SVEN TREITEL*.and R.M.ALFORD*.Synthetic Seismograms:A F-inite-Difference Approach,1976.41 (1):2-27
[26]Kneib G, Kerner C. Accurate and efficient seismic modeling in random media[J].G eophysics.1993,58(4):576-588
[27]Luo,M.,J.Cao.X.B.Xie.and R.S.Wu. Comparison of illumination analysis using one-way and full-wave propagators:74th Annual Internat. Mtg.,SEG.Expanded abstracts.2004.67-70
[28]Ostmo.S.,and Maao. F., Seismic illumination from finite-difference modeling,65th Conference&Technical Exhibition, EAGE, Expanded Abstracts.2003:146
[29]Peng C,Toksoz M N.An optimal absorbing boundary condition for elastic wave modeling,Geophysics,1995,60(1):296-301
[30]Ristow D.and Ruth T. Fourier finite-difference migration[J]. Geophysics,1994,59(12):1882-1893
[31]Robert C,Bjorn Eengquist.Absorbing boundary conditions for acoustic and ela-stic wave equations[J].Bull. Seism. Soc. Am.,1977,67(6):1529-1540
[32]Robert C,Bjorn Eengquist.Absorbing boundary conditions for wave equation migration[J].Geophysics,1980,45(3):895-904
[33]Roynods A C.Boundary conditions for the numerical solution of wave propagation problems.Geophysics,1978,43(6):1099-1110
[34]Strang.G.and Fix.G.J,1973,An analysis of the finite-element method:Englewood Cliffs,N.J.Prentice-Hall,Inc
[35]S A Orszag. Comparison of pseudo spectral and spectral approximation. Stu Appl Math:1972,51:253-259
[36]Schneider W A,Winbow G A.Efficient and accurate modeling of 3-D seismic il-lumination[J].Expanded Abstracts of 69th Annual Internation SEG Meeting,1999.633-636
[37]Steve Campbell,Bill Pramik,Bill Cafarelli. Comparative ray-based Illumination analysis[J].SEG Expanded abstracts.2002:41-44
[38]V. Bucha*,Department of Geophysics, Charles University, Prague, Czech Republic. Ray tracing computations in the smoothed SEG/EAGE Salt Model[J]. SEG Int'l Exposition and 74th Annual Meeting* Denver,Colorado*,2004
[39]W.E.A.Rietveld*,A.J.Berkhout*,C.P.A.Wapenaar*.Optinum seismic illumination of hydrocarbon reservoirs[J]. SEG Expanded abstracts,1992,57(10):1334-1345
[40]White.J.E,1965,Seismic waves:New York,McGrawHill Book Co.,Inc
[41]Xiaobi Xie, Shengwen Jin, Rushan Wu.Wave equation based illumination analysis. SEG Int'l Exposition and 74th Annual Meeting,2004
[42]Xiao-Bi Xie and Ru-Shan Wu,IGPP,University of California,Santa Cr uz,Michael Fehler and Lianjie Huang, Los Alamos National Labor atory Seismic resolution and illumination: A wave-equation-based analysis[J].SEG/Houston Annual Meeting.2005,1862-1867
[43]Zeng Y Q.He J Q,Liu Q H.The application of the perfectly matched layer in numerical modeling of wave propagation in poroelastic media[J].Geophysics.2001.66(4):1258-1266
[44]白焕新,刘洋.油气田与煤田地球物理勘探：地震照明技术.中国地球物理,2008：56-57
[45]白杰,彭德丽,蔡锡伟.基于射线追踪正演照明的观测系统分析研究[J].勘探地球物理学进展,2010,3(2)：87-93
[46]白振华.准噶尔盆地乌夏地区火山岩发育特征及其与油气成藏关系研究[J].中国西部科技，,2010,10(1)：31-34
[47]陈可洋.地震波数值模拟中差分近似的各向异性分析[J].石油物探,2010,49(1)：19-22
[48]陈可洋.起伏地表地震波正演数值模拟[J].高原地震,2010，,22(2)：39-46
[49]陈可洋，杨微.优化的三维地震波旁轴近似吸收边界条件[J].勘探地球物理进展,2009,32(3)：179-183
[50]陈生昌,毕亚杰,林晨.地震观测系统评价的波动方程方向照明方法[J].石油地球物理勘探,2009,44(6)：656-663
[51]陈生昌,马在田,吴如山.波动方程双程地下方向照明分析[J].同济大学学报,2007,35(5)：681-685
[52]董良国.弹性波数值模拟中的吸收边界条件[J].石油地球物理勘探.1999,34(1)：45-57
[53]董良国,吴晓丰,唐海忠,陈生昌,杨锴,刘守伟.逆掩推覆构造的地震波照明与观测系统优化[J].石油物探,2006,45(1)：40-48
[54]冯凯.三维地震观测系统最优化设计的方法研究.成都理工大学博士学位论文，2006
[55]冯伟，王华忠,吴如山,马在田.面向目标控制照明的合成波源偏移[J].石油物探,2004,43(3)：223-228
[56]冯英杰,杨长春,吴萍.地震波有限差分模拟综述[J].地球物理学进展，,2007,22(2)：487-491
[57]郭树祥.地震采集观测系统设计分析与探讨[J].油气地球物理,2009，7(4)：1-5
[58]国九英,周兴元,杨慧珠.三维f-x,y域随机噪音衰减[J].石油地球物理勘探,2005,302(2)：207-216
[59]韩刚.基于照明分析的观测系统优化.中国石油大学硕士学位论文,2009
[60]何光明,贺振华,黄得济,王翠华.非线性层次静校正在川西地区资料处理中的应用[J].石油物探,2006,45(1)：88-93
[61]何樵登,熊纬纲.地震勘探.地质出版社.1991.
[62]胡勐乾,邓志辉,陆远忠.有限元数值模拟方法在华北地区地震地质研究中的应用进展[J].地震地质.2010.32(1)：162-173
[63]胡天跃.地震资料叠前去噪技术的现状与未来[J].地球物理学进展,2002,17(2)：218-223
[64]栗宝鹃.基于射线理论的照明分析.中国石油大学硕士研究生学位论文.2004
[65]李佩,陈生昌,常鉴,徐国庆.基于波动理论的地震观测系统设计[J].2010,44(1)：203-208
[66]李庆忠.《走向精确勘探的道路高分辨率地震勘探系统工程剖析》.北京：石油工业出版社,1993
[67]李素华,卢齐军,王云专,范兴才,杨利伟.松辽盆地火山岩发育区数值模拟参数设计[J].石油地球物理勘探,2009,44(3)：319-326
[68]李万万.基于波动方程模拟的地震采集参数论证方法及应用研究.中国地质大学博士学位论文,2008
[69]李文杰,魏修成,宁俊瑞.一种新的弹性波数值模拟吸收边界条件[J].石油地球物理勘探,2009,44(4)：501-510
[70]李信富,李小凡,张美根.地震波数值模拟方法研究综述[J].防灾减灾工程学报,2007,27(2)：241-248
[71]李莹莹.中国石油发展现状、问题与前景分析[J].石油天然气,2010,32(12)：16-19
[72]刘传虎.模型正演技术在油气藏描述中的应用[J].石油地球物理勘探,1993,28(4)：477-484
[73]刘全新,杨宝莹.面向油气勘探需求加速物探技术发展[J].石油科技论坛,2008(5)：4-7
[74]刘洋,刘财,王典,李勤学,冯喧,王建民.时变中值滤波技术在地震随机噪声衰减中的应用[J].石油地球物理勘探，,2008，,43(3)：327-334
[75]刘洋,魏修成.复杂构造中地震波传播数值模拟[J].新疆石油地质，,2008,29(1)：12-14
[76]刘振武,撒利明,董世泰,邓志文,徐光成.中国石油物探技术现状与发展方向.石油勘探与开发,2010,37(1)：1-10
[77]刘振武,撒利明,张听,董世泰,甘利灯.中国石油开发地震技术应用现状和未来发展建议.石油学报,，2009,，30(5)：711-717
[78]刘志敏,包吉山.高阶有限差分地震波场正演[J].石油地球物理勘探,1989,24(1)：1-5
[79]吕公河,尹成,周星合.基于采集目标的地震照明度的精确模拟[J].石油地球物理勘探,2006,42(3)：258-261
[80]吕公河,徐锦玺,张在武,卢湘鹏,何京国.地震采集技术.2008,，78届SEG年会
[81]陆基孟.地震勘探原理[M].山东东营：石油大学出版社,1993
[82]罗大清,宋炜,吴律,柏森林.加权方向校正及插值预测吸收边界条件的三维波动方程叠前数值模拟[J].石油地球物理勘探，1999,34(3)：275-284
[83]贾晓峰，胡天跃,王润秋.无单元法用于地震波波动方程模拟与成像[J].地球物理学进展，2006.21(1)：11-17
[84]蒋先艺.基于二维与三维复杂结构模型正演的地震数据采集设计方法研究.成都理工大学博士论文,2004
[85]马中良,曾溅辉,张善文,王永诗,王洪玉,刘惠民.砂岩透镜体油运移过程模拟及成藏主控因素分析[J].岩性油气藏,2008,20(1).
[86]马永旺,詹正彬,顾汉明.伪谱法地震波传播数值模拟[J].石油物探,2001,40(2)：42-48
[87]牛滨华,吕景贵,孙春岩,黄新武,章明昱,沈操,吴亚东.叠前面波干扰压制方法的研究与应用[J].现代地质,2001,15(3)：326-332
[88]潘文锋,贾随良,梁春生,王小山,皇甫煊.灵活二维观测系统设计及其应用[J].石油物探,2003,42(1)：39-41
[89]裴正林,牟永光.地震波传播的数值模拟[J].地球物理学进展,2004,19(4)：933-941
[90]彭真明,梁国胜,王纬,樊振现,陈艳玲.砂、泥岩互层透镜体模型正演及分析[J].物探化探计算技术,2001,23(2)：101-104
[91]皮红梅.双程波动方程数值模拟和照明分析方法研究.吉林大学博士学位论文.2009.
[92]浦锡锋,张海波,田宙.岩石中强爆炸动力学过程数值模拟[J].岩石力学与工程学报,2005,24(1)：5007-5012
[93]钱光萍,康家光,王紫捐.基于模型的地震采集参数分析及应用研究[J].物探化探计算技术.2001,23(2)：109-114
[94]秦顺亭.当前中国西部新区物探工作的几点思考[J].中国西部石油地质,2005,1(2)：124-127
[95]单刚义.转换波处理与解释在储层预测中的应用.吉林大学硕士学位论文，2007.
[96]单刚义,韩立国,张丽华,董世学.压电式检波器在高分辨率地震勘探中的实验研究.石油物探，,2009,44(1)：91-95
[97]单刚义,韩立国,张丽华，王艳玲.转换波在储层预测中的研究分析.地球物理学进展.2010，,25(1)：282-287
[98]单联瑜,刘洪,匡斌,王华忠,沈财余,许学平.基于胜利典型地质模型的波动方程地震波照明分析研究[J].石油地球物理勘探,2009,44(1)：1-8
[99]孙书荣,凡友华.地震波数值模拟技术发展现状[J].油气地球物理,2009,7(1)：18-23
[100]孙伟家,符力耘，碗学检,基于相位编码的地震照明度分析[J].石油地球物理勘探，2007，42(5)：539-545
[101]谭绍泉.深层油气地震勘探数据采集与处理方法[M].北京：地质出版社.2004
[102]王春燕.高阶交错网格有限差分地震波场计算.成都理工大学硕士学位论文，2007.
[103]王德志,贾烈明,王帮助,尹军杰,刘学伟.基于勘探目标的观测系统设计[J].石油地球物理勘探,2006,41(5)：498-503
[104]王海波,刘兵,于建华,王长江.正演模拟在大平房火山岩发育区地震信息采集中的应用[J].新疆石油地质,2010,31(2)：200-202
[105]王建伟,宋书君,王新征,赵密福,万从礼,牛富玲.东营凹陷牛庄洼陷砂岩透镜体的成藏机理[J].石油学报,2007,28(5)：39-45
[106]王雪秋.《利用有限积分法变换一有限差分方法模拟地震波场》.吉林大学硕士学位论文,2002.
[107]王月,冯连勇,牛燕.世界石油物探市场及物探技术发展趋势分析[J].石油地球物理勘探,2008,43(6)：736-7443
[108]熊翥.山前冲断带物探技术改进的思路[J].勘探技术,2005,2：28-33
[109]徐义,张剑锋.地震波数值模拟的非规则网格PML吸收边界[J].地球物理学报,2008,51(5)：1520-1526
[110]徐果明,卫山,高尔根,林庆忠,蒋先艺,罗开云.二维复杂介质的块状建模及射线追踪[J].石油地球物理勘探,，2001,36(2)：213-220
[111]徐怀大.寻找非构造油气藏的新思路[J].勘探家,1996,1(1)：43-48
[112]尹成,吕公河,田继东,徐锦玺,尚应军.基于地球物理目标参数的三维观测系统优化设计[J].石油物探,2006,45(1)：74-78
[113]尹吴海,甄文胜,魏国显,张立军.三维地震采集观测系统设计技术-以柴达木盆地西部地区为例[J].勘探地球物理进展,2007,30(2)：111-118
[114]张付生.焉耆盆地山前推覆体区地震勘探技术研究.中国地质大学工程硕士学位论文,2008
[115]张海燕,李庆忠.几种常用解析子波的特性分析[J].石油地球物理勘探,2007,42(6)：651-660
[116]张加海,王敬绵,张光德,安峰,董长安.炸药震源激发参数与场地质点振动响应关系探讨[J].油气地质与采收率,2010,17(2)：57-61
[117]张金海,王卫民,赵连锋等.傅里叶有限差分法三维波动方程正演模拟[J].地球物理学报,2007,50(6)：1854-1862
[118]张军华.吕宁,田连玉,陆文志，钟磊.地震资料去噪方法技术综合评述[J].地球物理学进展，2006，21(2)：546-553
[119]张孟，李亚林，李大军,宋红玲.川东高陡构造石灰岩出露区的地震采集技术及其应用效 果[J].天然气工业,2007,(增刊A)：82-85
[120]张晓斌,王勤耕,陈江力,李彪,雷迎春.西部地区复杂山地地震采集技术[J].天然气工业,2007,(A)：76-78
[121]张政威,陈亚强,张付生.焉耆盆地西部逆冲推覆带二维地震采集方法研究[J].河南石油,2005,19(3)：14-16
[122]赵爱华,张美根,丁志峰.横向各向同性介质中地震波走时模拟[J].地球物理学报,2006,49(6)：1762-1769
[123]赵殿栋,郭建,王咸彬,魏福吉,吴长祥,贾洪顺,丁士平,李鼎民,董良国.基于模型面向目标的观测系统优化设计技术[J].中国西部石油地质.2006,2(2)：119-122
[124]赵海波,王秀明,王东,陈浩.完全匹配层吸收边界在孔隙介质弹性波模拟中的应用[J].地球物理学报,2007,50(2)：581-591
[125]赵明月,李桐林,刘希芳,林君.长春科技大学学报[J].1998,28(4)：448-452
[126]中石油勘探与生产分公司.中国中西部前陆盆地冲断带油气勘探文集[M].北京：石油工业出版社,2002：99-112
[127]周星合.地质模型创建及地震采集设计中地震波照明度模拟与分析研究,西南石油大学硕士论文，2006
[128]朱号锋,靳平,肖卫国.硬岩中地下爆炸震源函数的数值模拟[J].爆炸与冲击,2009.,29(6)：648-653
[2]Alterman Z,Karal F C. Propagation of seismic wave in layered media by finite difference methods[J].BSSA,1968,58(1):367-398
[3]Aminzadeh F.,Brac J.,Kunz T.3-D Salt and Overthrust Models, SEG/EAGE 3-D Modeling Series No.1.Soc.Explor.Geophysicists,Tulsa,1997
[4]Bear G,Lu C,Lu R.The construction of subsurface illumination and amplitude maps via ray tracing.The leading Edge,2000,19(7):726-731
[5]Berenger J P.A perfectly matched layer for the absorption of elect romagnetic waves.Journal of Computational Physics,1994,114(2):185-200
[6]Berenger J P.Three-dimensional perfectly matched layer for the absorption of electromagnetic waves.Journal of Computational Physics,1996,127(2):363-379
[7]Boor.D.M.1970.Finite-difference solution to the equations of elastic wave propagation with application to Love waves over dipping interfaces:Ph.D thesis.M.I.T
[8]B Fornberg.The pseudospectral method:Acurate representation of interfaces in elastic wave calculations. Geophysics.1988,53(5):625-637
[9]Chen K H.Seismic modeling:A finite element approach[J].SEG/China Geophysical Society Metting,Beijing,China,1981.
[10]Clayton R.Engquist B.Absorbing boundary conditions for acoustic and elastic wave equations [J].The Bulletin of the Seismological Society of America,1977,67:1529-1540
[11]Clough R W.The finite elementmethod in plane stress analysis[A].1960, In:Proceedings of 2nd ASCE Conference on Electronic Computation Pittsburgh, ASCE.
[12]Collino F,Tsogka C. Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media.Geo-physics,2001,66(1):294-307
[13]David Muerdter,Diamond Geoscience Research Corporation,Redmond, Washington,U.S.DAVIS RATCLIFF,Diamond Geophysical Service Corporation Houston.Texas.U.S.Understanding subsalt illumination through ray-trace modeling.Part 1:Simple 2-D salt models. The Leading Edge.2001,1(1):578-586
[14]David Muerdter等著,娄晓东译,孙夕平,赵该善校.通过射线追踪模拟理解盐下照射[M].石油物探译丛,2001,10(5)：59-70
[15]Falk J, Tessmer E,Gajewski D. Efficient finite-difference modeling of seismic waves using locally adjustable time steps [J]. Geophysical Prospecting,1998,46:603-616.
[16]Gelius L.J.Lecomte I.and Tabti H.Analysis of the resolution function in seismic prestack depth imaging.Geophys.Prosp.2002.50:505-515
[17]Greg J.Schurter.Yayat Supriatna,Angke Nuraeni.and Supriyono.Bp In donesia Carl J.Regone and Nurul Kabir.Bp EPTG.Houston.A 3D finitedifference modeling study of seismic imaging challenges in Bintuni Bay,Irian Jaya Barat.2009,1008-1018
[18]Gustavo Catao Alves*,Andre Bulcao, Djalma Manoel Soares Filho, Carlos Eduardo Theodoro,Luiz Alberto Santos.Marcos Antonio Gallotti Guima raes,Petrobras.Target oriented illumination analysis using wave equation [J]. SEG Las Vegas Annual Meeting.2008,163-167
[19]Hastings F.Schneider J B.Broschat S L.Application of the perfectly matchedlayer(PML) absorbing boundary condition to elastic wave pro pagation[J] Journal of the Acoustic Society of America,1996,100(5):3061-3069
[20]Hedstrom G W.Nonreflecting boundary conditions for nonlinear hype rbolic systerms.J Comput Phys.1979.30:222-237
[21]Higdon R L. Absorbing boundary conditions for elastic wave[J]. Geophysics.1991.56(2):231-241
[22]Hoffmann J.Illumination.resolution and image quality of PP- and PS- waves for survey planning.The leading Edge.2001.20(9):1008-1012
[23]Igel H, Weber M. P_SV wave propagation in the Earth s mantle using finite differences:Application to heterogeneous lower most mantle structure [J].G eophys,Res.Lett,1966,23:415-418
[24]Jastram C,Tessemer E.Elastic modeling on a grid with vertically varying spacing[J].Geophysical prospecting,1994,42:357-370
[25]K.R.KELLY*,R.W.WARD.SVEN TREITEL*.and R.M.ALFORD*.Synthetic Seismograms:A F-inite-Difference Approach,1976.41 (1):2-27
[26]Kneib G, Kerner C. Accurate and efficient seismic modeling in random media[J].G eophysics.1993,58(4):576-588
[27]Luo,M.,J.Cao.X.B.Xie.and R.S.Wu. Comparison of illumination analysis using one-way and full-wave propagators:74th Annual Internat. Mtg.,SEG.Expanded abstracts.2004.67-70
[28]Ostmo.S.,and Maao. F., Seismic illumination from finite-difference modeling,65th Conference&Technical Exhibition, EAGE, Expanded Abstracts.2003:146
[29]Peng C,Toksoz M N.An optimal absorbing boundary condition for elastic wave modeling,Geophysics,1995,60(1):296-301
[30]Ristow D.and Ruth T. Fourier finite-difference migration[J]. Geophysics,1994,59(12):1882-1893
[31]Robert C,Bjorn Eengquist.Absorbing boundary conditions for acoustic and ela-stic wave equations[J].Bull. Seism. Soc. Am.,1977,67(6):1529-1540
[32]Robert C,Bjorn Eengquist.Absorbing boundary conditions for wave equation migration[J].Geophysics,1980,45(3):895-904
[33]Roynods A C.Boundary conditions for the numerical solution of wave propagation problems.Geophysics,1978,43(6):1099-1110
[34]Strang.G.and Fix.G.J,1973,An analysis of the finite-element method:Englewood Cliffs,N.J.Prentice-Hall,Inc
[35]S A Orszag. Comparison of pseudo spectral and spectral approximation. Stu Appl Math:1972,51:253-259
[36]Schneider W A,Winbow G A.Efficient and accurate modeling of 3-D seismic il-lumination[J].Expanded Abstracts of 69th Annual Internation SEG Meeting,1999.633-636
[37]Steve Campbell,Bill Pramik,Bill Cafarelli. Comparative ray-based Illumination analysis[J].SEG Expanded abstracts.2002:41-44
[38]V. Bucha*,Department of Geophysics, Charles University, Prague, Czech Republic. Ray tracing computations in the smoothed SEG/EAGE Salt Model[J]. SEG Int'l Exposition and 74th Annual Meeting* Denver,Colorado*,2004
[39]W.E.A.Rietveld*,A.J.Berkhout*,C.P.A.Wapenaar*.Optinum seismic illumination of hydrocarbon reservoirs[J]. SEG Expanded abstracts,1992,57(10):1334-1345
[40]White.J.E,1965,Seismic waves:New York,McGrawHill Book Co.,Inc
[41]Xiaobi Xie, Shengwen Jin, Rushan Wu.Wave equation based illumination analysis. SEG Int'l Exposition and 74th Annual Meeting,2004
[42]Xiao-Bi Xie and Ru-Shan Wu,IGPP,University of California,Santa Cr uz,Michael Fehler and Lianjie Huang, Los Alamos National Labor atory Seismic resolution and illumination: A wave-equation-based analysis[J].SEG/Houston Annual Meeting.2005,1862-1867
[43]Zeng Y Q.He J Q,Liu Q H.The application of the perfectly matched layer in numerical modeling of wave propagation in poroelastic media[J].Geophysics.2001.66(4):1258-1266
[44]白焕新,刘洋.油气田与煤田地球物理勘探：地震照明技术.中国地球物理,2008：56-57
[45]白杰,彭德丽,蔡锡伟.基于射线追踪正演照明的观测系统分析研究[J].勘探地球物理学进展,2010,3(2)：87-93
[46]白振华.准噶尔盆地乌夏地区火山岩发育特征及其与油气成藏关系研究[J].中国西部科技，,2010,10(1)：31-34
[47]陈可洋.地震波数值模拟中差分近似的各向异性分析[J].石油物探,2010,49(1)：19-22
[48]陈可洋.起伏地表地震波正演数值模拟[J].高原地震,2010，,22(2)：39-46
[49]陈可洋，杨微.优化的三维地震波旁轴近似吸收边界条件[J].勘探地球物理进展,2009,32(3)：179-183
[50]陈生昌,毕亚杰,林晨.地震观测系统评价的波动方程方向照明方法[J].石油地球物理勘探,2009,44(6)：656-663
[51]陈生昌,马在田,吴如山.波动方程双程地下方向照明分析[J].同济大学学报,2007,35(5)：681-685
[52]董良国.弹性波数值模拟中的吸收边界条件[J].石油地球物理勘探.1999,34(1)：45-57
[53]董良国,吴晓丰,唐海忠,陈生昌,杨锴,刘守伟.逆掩推覆构造的地震波照明与观测系统优化[J].石油物探,2006,45(1)：40-48
[54]冯凯.三维地震观测系统最优化设计的方法研究.成都理工大学博士学位论文，2006
[55]冯伟，王华忠,吴如山,马在田.面向目标控制照明的合成波源偏移[J].石油物探,2004,43(3)：223-228
[56]冯英杰,杨长春,吴萍.地震波有限差分模拟综述[J].地球物理学进展，,2007,22(2)：487-491
[57]郭树祥.地震采集观测系统设计分析与探讨[J].油气地球物理,2009，7(4)：1-5
[58]国九英,周兴元,杨慧珠.三维f-x,y域随机噪音衰减[J].石油地球物理勘探,2005,302(2)：207-216
[59]韩刚.基于照明分析的观测系统优化.中国石油大学硕士学位论文,2009
[60]何光明,贺振华,黄得济,王翠华.非线性层次静校正在川西地区资料处理中的应用[J].石油物探,2006,45(1)：88-93
[61]何樵登,熊纬纲.地震勘探.地质出版社.1991.
[62]胡勐乾,邓志辉,陆远忠.有限元数值模拟方法在华北地区地震地质研究中的应用进展[J].地震地质.2010.32(1)：162-173
[63]胡天跃.地震资料叠前去噪技术的现状与未来[J].地球物理学进展,2002,17(2)：218-223
[64]栗宝鹃.基于射线理论的照明分析.中国石油大学硕士研究生学位论文.2004
[65]李佩,陈生昌,常鉴,徐国庆.基于波动理论的地震观测系统设计[J].2010,44(1)：203-208
[66]李庆忠.《走向精确勘探的道路高分辨率地震勘探系统工程剖析》.北京：石油工业出版社,1993
[67]李素华,卢齐军,王云专,范兴才,杨利伟.松辽盆地火山岩发育区数值模拟参数设计[J].石油地球物理勘探,2009,44(3)：319-326
[68]李万万.基于波动方程模拟的地震采集参数论证方法及应用研究.中国地质大学博士学位论文,2008
[69]李文杰,魏修成,宁俊瑞.一种新的弹性波数值模拟吸收边界条件[J].石油地球物理勘探,2009,44(4)：501-510
[70]李信富,李小凡,张美根.地震波数值模拟方法研究综述[J].防灾减灾工程学报,2007,27(2)：241-248
[71]李莹莹.中国石油发展现状、问题与前景分析[J].石油天然气,2010,32(12)：16-19
[72]刘传虎.模型正演技术在油气藏描述中的应用[J].石油地球物理勘探,1993,28(4)：477-484
[73]刘全新,杨宝莹.面向油气勘探需求加速物探技术发展[J].石油科技论坛,2008(5)：4-7
[74]刘洋,刘财,王典,李勤学,冯喧,王建民.时变中值滤波技术在地震随机噪声衰减中的应用[J].石油地球物理勘探，,2008，,43(3)：327-334
[75]刘洋,魏修成.复杂构造中地震波传播数值模拟[J].新疆石油地质，,2008,29(1)：12-14
[76]刘振武,撒利明,董世泰,邓志文,徐光成.中国石油物探技术现状与发展方向.石油勘探与开发,2010,37(1)：1-10
[77]刘振武,撒利明,张听,董世泰,甘利灯.中国石油开发地震技术应用现状和未来发展建议.石油学报,，2009,，30(5)：711-717
[78]刘志敏,包吉山.高阶有限差分地震波场正演[J].石油地球物理勘探,1989,24(1)：1-5
[79]吕公河,尹成,周星合.基于采集目标的地震照明度的精确模拟[J].石油地球物理勘探,2006,42(3)：258-261
[80]吕公河,徐锦玺,张在武,卢湘鹏,何京国.地震采集技术.2008,，78届SEG年会
[81]陆基孟.地震勘探原理[M].山东东营：石油大学出版社,1993
[82]罗大清,宋炜,吴律,柏森林.加权方向校正及插值预测吸收边界条件的三维波动方程叠前数值模拟[J].石油地球物理勘探，1999,34(3)：275-284
[83]贾晓峰，胡天跃,王润秋.无单元法用于地震波波动方程模拟与成像[J].地球物理学进展，2006.21(1)：11-17
[84]蒋先艺.基于二维与三维复杂结构模型正演的地震数据采集设计方法研究.成都理工大学博士论文,2004
[85]马中良,曾溅辉,张善文,王永诗,王洪玉,刘惠民.砂岩透镜体油运移过程模拟及成藏主控因素分析[J].岩性油气藏,2008,20(1).
[86]马永旺,詹正彬,顾汉明.伪谱法地震波传播数值模拟[J].石油物探,2001,40(2)：42-48
[87]牛滨华,吕景贵,孙春岩,黄新武,章明昱,沈操,吴亚东.叠前面波干扰压制方法的研究与应用[J].现代地质,2001,15(3)：326-332
[88]潘文锋,贾随良,梁春生,王小山,皇甫煊.灵活二维观测系统设计及其应用[J].石油物探,2003,42(1)：39-41
[89]裴正林,牟永光.地震波传播的数值模拟[J].地球物理学进展,2004,19(4)：933-941
[90]彭真明,梁国胜,王纬,樊振现,陈艳玲.砂、泥岩互层透镜体模型正演及分析[J].物探化探计算技术,2001,23(2)：101-104
[91]皮红梅.双程波动方程数值模拟和照明分析方法研究.吉林大学博士学位论文.2009.
[92]浦锡锋,张海波,田宙.岩石中强爆炸动力学过程数值模拟[J].岩石力学与工程学报,2005,24(1)：5007-5012
[93]钱光萍,康家光,王紫捐.基于模型的地震采集参数分析及应用研究[J].物探化探计算技术.2001,23(2)：109-114
[94]秦顺亭.当前中国西部新区物探工作的几点思考[J].中国西部石油地质,2005,1(2)：124-127
[95]单刚义.转换波处理与解释在储层预测中的应用.吉林大学硕士学位论文，2007.
[96]单刚义,韩立国,张丽华,董世学.压电式检波器在高分辨率地震勘探中的实验研究.石油物探，,2009,44(1)：91-95
[97]单刚义,韩立国,张丽华，王艳玲.转换波在储层预测中的研究分析.地球物理学进展.2010，,25(1)：282-287
[98]单联瑜,刘洪,匡斌,王华忠,沈财余,许学平.基于胜利典型地质模型的波动方程地震波照明分析研究[J].石油地球物理勘探,2009,44(1)：1-8
[99]孙书荣,凡友华.地震波数值模拟技术发展现状[J].油气地球物理,2009,7(1)：18-23
[100]孙伟家,符力耘，碗学检,基于相位编码的地震照明度分析[J].石油地球物理勘探，2007，42(5)：539-545
[101]谭绍泉.深层油气地震勘探数据采集与处理方法[M].北京：地质出版社.2004
[102]王春燕.高阶交错网格有限差分地震波场计算.成都理工大学硕士学位论文，2007.
[103]王德志,贾烈明,王帮助,尹军杰,刘学伟.基于勘探目标的观测系统设计[J].石油地球物理勘探,2006,41(5)：498-503
[104]王海波,刘兵,于建华,王长江.正演模拟在大平房火山岩发育区地震信息采集中的应用[J].新疆石油地质,2010,31(2)：200-202
[105]王建伟,宋书君,王新征,赵密福,万从礼,牛富玲.东营凹陷牛庄洼陷砂岩透镜体的成藏机理[J].石油学报,2007,28(5)：39-45
[106]王雪秋.《利用有限积分法变换一有限差分方法模拟地震波场》.吉林大学硕士学位论文,2002.
[107]王月,冯连勇,牛燕.世界石油物探市场及物探技术发展趋势分析[J].石油地球物理勘探,2008,43(6)：736-7443
[108]熊翥.山前冲断带物探技术改进的思路[J].勘探技术,2005,2：28-33
[109]徐义,张剑锋.地震波数值模拟的非规则网格PML吸收边界[J].地球物理学报,2008,51(5)：1520-1526
[110]徐果明,卫山,高尔根,林庆忠,蒋先艺,罗开云.二维复杂介质的块状建模及射线追踪[J].石油地球物理勘探,，2001,36(2)：213-220
[111]徐怀大.寻找非构造油气藏的新思路[J].勘探家,1996,1(1)：43-48
[112]尹成,吕公河,田继东,徐锦玺,尚应军.基于地球物理目标参数的三维观测系统优化设计[J].石油物探,2006,45(1)：74-78
[113]尹吴海,甄文胜,魏国显,张立军.三维地震采集观测系统设计技术-以柴达木盆地西部地区为例[J].勘探地球物理进展,2007,30(2)：111-118
[114]张付生.焉耆盆地山前推覆体区地震勘探技术研究.中国地质大学工程硕士学位论文,2008
[115]张海燕,李庆忠.几种常用解析子波的特性分析[J].石油地球物理勘探,2007,42(6)：651-660
[116]张加海,王敬绵,张光德,安峰,董长安.炸药震源激发参数与场地质点振动响应关系探讨[J].油气地质与采收率,2010,17(2)：57-61
[117]张金海,王卫民,赵连锋等.傅里叶有限差分法三维波动方程正演模拟[J].地球物理学报,2007,50(6)：1854-1862
[118]张军华.吕宁,田连玉,陆文志，钟磊.地震资料去噪方法技术综合评述[J].地球物理学进展，2006，21(2)：546-553
[119]张孟，李亚林，李大军,宋红玲.川东高陡构造石灰岩出露区的地震采集技术及其应用效 果[J].天然气工业,2007,(增刊A)：82-85
[120]张晓斌,王勤耕,陈江力,李彪,雷迎春.西部地区复杂山地地震采集技术[J].天然气工业,2007,(A)：76-78
[121]张政威,陈亚强,张付生.焉耆盆地西部逆冲推覆带二维地震采集方法研究[J].河南石油,2005,19(3)：14-16
[122]赵爱华,张美根,丁志峰.横向各向同性介质中地震波走时模拟[J].地球物理学报,2006,49(6)：1762-1769
[123]赵殿栋,郭建,王咸彬,魏福吉,吴长祥,贾洪顺,丁士平,李鼎民,董良国.基于模型面向目标的观测系统优化设计技术[J].中国西部石油地质.2006,2(2)：119-122
[124]赵海波,王秀明,王东,陈浩.完全匹配层吸收边界在孔隙介质弹性波模拟中的应用[J].地球物理学报,2007,50(2)：581-591
[125]赵明月,李桐林,刘希芳,林君.长春科技大学学报[J].1998,28(4)：448-452
[126]中石油勘探与生产分公司.中国中西部前陆盆地冲断带油气勘探文集[M].北京：石油工业出版社,2002：99-112
[127]周星合.地质模型创建及地震采集设计中地震波照明度模拟与分析研究,西南石油大学硕士论文，2006
[128]朱号锋,靳平,肖卫国.硬岩中地下爆炸震源函数的数值模拟[J].爆炸与冲击,2009.,29(6)：648-653