多尺度地震资料正反演方法研究
|
摘要
油藏地球物理是地球物理技术面向油气藏开发生产综合应用的新型学科,油藏地球物理的内涵包括油藏描述与油藏管理,利用地球物理方法进行油藏圈定和描述,并在油藏开发开采过程中监测油藏变化。油藏地球物理依托于高精度三维地震、井中地震、时延地震、多波多分量地震等多种地震技术的发展,其中高精度三维地震与其它技术的综合应用是油藏地球物理的核心技术。开展多尺度地震资料正反演方法研究,提高多尺度地震资料协同解决油藏地质问题的能力;基于流变学理论研究高精度粘弹性波动方程正演模拟方法,准确认识地质体在不同观测系统下的多尺度地震响应特征;利用贝叶斯反演理论建立多尺度地震资料联合反演方法,减小地震反演多解性、提高反演分辨率,充分发挥地震资料为油田开发服务的作用。
地震波在实际介质中的固有传播衰减特征很复杂,导致在理论上出现了多种描述粘弹性介质的衰减模型,但大部分用于粘弹性介质波动方程时间域正演模拟的衰减模型对其常Q特征都存在较大的近似。基于流变学理论研究,开展了标准线性体、Kelvin等多种衰减模型对粘弹性介质波动方程常Q正演模拟的适应性研究,分析了能够用于精确描述粘弹性介质中地震波衰减规律的多种广义流变学模型特征,并证明了广义流变模型之间具有等价性。分析了不同阶次广义Maxwell流变模型对粘弹性介质常Q衰减特征的拟合程度,完善了利用广义流变模型开展时间域粘弹性介质波动方程常Q正演模拟的基础,并将模型特征与常用的Futterman等衰减模型特征进行比较。对多种衰减模型开展了一维平面波正演模拟分析,解释了由于模型不同所引起的地震波传播衰减特征差异,指出了广义流变模型在波动方程时间域正演模拟中的优势。
将广义流变模型引入到时间域粘弹性介质波动方程常Q正演模拟过程中,并分别针对广义Maxwell流变模型和广义标准线性体两种模型特征进行分析,对前者考虑了两种非弹性函数定义方式来建立对应的波动方程,分别采用弹性模量和松弛模量来描述基于广义标准线性体模型的波动方程,重新定义记忆变量使得通过两种广义流变模型所建立的波动方程表达式得到统一。完整的给出了基于两种模量表示的各种广义流变模型所对应的位移-应力、位移、位移-速度-应力和速度-应力方程,指出了早期波动方程由于采用松弛模量进行表示所引起的误解。基于五阶广义Maxwell流变模型建立了高精度粘弹性介质高阶交错网格有限差分波动方程常Q正演模拟表达式,采用PML边界条件有效的消除人工边界反射,讨论了模拟过程中的稳定性条件和震源问题,将有限差分正演模拟结果与二维均匀粘弹性介质中的解析解进行对比,分析了高精度粘弹性介质波动方程有限差分数值模拟方法的精度。
采用高精度粘弹性介质常Q波动方程正演模拟方法开展模型响应特征分析认识到地震波在粘弹性介质中具有较强的衰减频散特征,随着Q值的减小和传播距离的增大,地震波产生振幅衰减、相位改变、主频降低等现象。针对不同的观测系统开展正演模拟,通过波场、共炮集及偏移成像剖面特征进行对比分析可知,地面地震、VSP和井间地震三种多尺度地震资料之间既有联系也有差异,联系在于多尺度地震资料是相同地下地质体在相应观测方式下的地震波响应特征;差异在于资料所包含的波场信息丰富程度不一样、反映地质体信息的范围不一样、分辨地质体细节特征的能力也不一样。加入近地表强衰减层的影响表明,多尺度地震资料之间存在一定的相位差异,进一步认识了实际多尺度地震资料之间既具有一致性又存在一定差异的特征。
基于贝叶斯反演理论,通过似然函数有机的将地面地震、VSP和井间地震三种多尺度地震资料进行整合,建立起基于修正柯西分布的多尺度地震资料联合反演方法体系,并完善了联合反演流程。利用一维模型验证了联合反演算法的正确性,加噪后通过卡方分布测试了联合反演参数的合理选择范围,并采用二维多尺度正演模拟资料验证了联合反演方法的有效性及对多尺度资料的适应性。对联合反演方法进行退化建立了时间域单一尺度地震资料反演方法;针对连井井间地震剖面的特征重点强调了高精度深时转换、井地时深关系一致性及加强横向约束等改进策略,完善了深度域井间地震资料的反演方法和流程。
将多尺度地震资料反演方法对胜利油田垦71区块开展了应用研究,利用深度域井间地震反演方法获得连井井间地震波阻抗反演剖面,从该反演结果中提取的虚拟井能够有效的匹配地面地震资料反射特征。虚拟井约束下的井地联合反演只能够提高地面地震反演精度,而基于贝叶斯理论框架建立的多尺度地震资料联合反演方法则能够有效的提高地面地震资料反演结果的分辨率,实现多尺度地震资料特征的融合。
Reservoir geophysics is a new discipline that geophysics technique towards the integrated application of oilfield development and production, the reservoir description and reservoir management is included. The geophysics technique is adapted to determin the reservoir, descript the reservoir and monitor the change in the production. It is relying on the development of 3-D high-precision seismic technique, borehole seismic technique, time-lapse seismic technique and multi-wave/multi-component technique, the integrated research of high-precision 3-D seismic and other methods is the key technology of reservoir geophysics. Research the modeling and inversion method of multiscale seismic data is improving the capacity of seismic data to resolve the geology problem. Develop the high precision seismic modeling method of wave-equation in viscoelastic media, determine the mulitscale seismic response characteristics of geology with different observation system. Build up the multiscale seismic data joint inversion method based on the Bayes theory, overcome the multiple solution of inversion problem and improve the resolution of seismic data inversion, then service for the oilfield development with seismic data.
The nature attenuation law of seismic propagating in actual media is very complicated. many attenuation model of attenuation model is present in theory. The approximation to Q is obvious to most models in seismic wave-equation modeling at viscoelastic media in time domain. Based on the rheology theory, the adaptability of SLS (Standard Linear Solid), Kelvin and other attenuation model is researched in viscoelastic seismic modeling. The generalized rheology models that can be use to accurately describe the seismic wave propagating and attenuation characteristics have been analysied, the equivalence of generalize Maxwell model and generalized SLS model is proofed. The fitting degree to the constant Q of different order generalized Maxwell model is analysized and compared with the Futterman model and other models, the foundation is perfecing to modeling the constant Q characteristics in viscoelastic media. One dimension plane wave modeling is done with these attenuation models, the difference of seismic wave propagation and attenuation characteristic can be inteprated with these results, the advantage of generalized rheology model is clear.
The generalized rheology model is used in wave-equation seismic modeling at viscoelastic media in time domain. For the generalized Maxwell model, two kinds of unelastic function is adapted to derive the wave-equation. The unelastic modulus and relaxed modulus have been used to build up the wave-equation for generalized SLS model. The unified expression between these two kinds of generalized rheology model after the memory variable is redefined. The expression of displacement equation, the displacement-stress equation, the displacement-velocity-stress equation and the velocity-stress equation is presented according different modulus form.With the equation expansion or variable transform, the commom viscoelastic wave equation in articles can be get, the disadvantage of early equation expressed by relaxed modulus is clear. The high-order staggered grid finite difference scheme of high-precision viscoelastic wave equation based on five-order generalized Maxwell model is present. It is the high-precision wave-equation modeling method to constant Q in viscoelastic media. The PML boundary condition is used to absorb the artificial boundary reflection, the stability condition and source problem is discussed, the high-order finite difference scheme is accurate enough according to the comparison with the analytical resolution of 2D homogeneous viscoelastic media.
According the high-precision viscoelastic modeling, the attenuation and dispersion characterstic of seismic wave is obvious in viscoelastic media. With the Q decrease and the propagation distance increase, the amplitude is attenuated. the phase is lagged. the dominant frequency is decreased. Modeling with different observation system and compare with the wave field, common shot gather and migrated section characteristic, the relationship and difference among multiscale seismic data can be recognized. The relationship is multiscale seimic data is the seismic wave reflection response of the same geological body under different observation way. The difference is represented at the wave field information richness, the geological range and the detail characteristic in the multiscale seismic data. When the near-surface strong attenuation layer is considered, the phase difference is obvious between surface seismic data and others, the relationship and difference among multiscale seismic data can be inteprated.
Based on the Bayesian theorem, the new likelihood fundction is integrated with surface, vsp and crosswell seismic data, the multiscale seismic data joint inversion framework is built with modified Cauchy distribution and the inversion flow is improved. The one-dimension model testing shows that the method is right, the chi-square distribution is used to select the parameter range of the joint inversion with gauss noise, the joint inversion flow is valid to the two-dimension multiscale seismic data modeling according to the wave equation. For the singlescale seismic data, the crosswell seismic data inversion workflow in depth domain is improved. The high-precision interpolation function was adapted in depth-time transformation, the reflection characteristic consistence between surface seismic data and crosswell seismic data is emphasized, the transverse constraint is strengthened according to the characteristic of multi-crosswell seismic data.
The multiscale sesmic data inversion method is tested in Ken71 area of Shengli oilfield. The reasonable inversion result of multi-crosswell seismic data is get, the pseudo well exracted from it can mathch the characteristics of surface seismic data well. The constraint inversion with pseudo well can improve the precision of result. The multiscale seismic joint inversion method can reduce the multiplicity of solution, improve the inversion resolution, the multiscale data can be matched with each other according to the joint inversion result.
地震波在实际介质中的固有传播衰减特征很复杂,导致在理论上出现了多种描述粘弹性介质的衰减模型,但大部分用于粘弹性介质波动方程时间域正演模拟的衰减模型对其常Q特征都存在较大的近似。基于流变学理论研究,开展了标准线性体、Kelvin等多种衰减模型对粘弹性介质波动方程常Q正演模拟的适应性研究,分析了能够用于精确描述粘弹性介质中地震波衰减规律的多种广义流变学模型特征,并证明了广义流变模型之间具有等价性。分析了不同阶次广义Maxwell流变模型对粘弹性介质常Q衰减特征的拟合程度,完善了利用广义流变模型开展时间域粘弹性介质波动方程常Q正演模拟的基础,并将模型特征与常用的Futterman等衰减模型特征进行比较。对多种衰减模型开展了一维平面波正演模拟分析,解释了由于模型不同所引起的地震波传播衰减特征差异,指出了广义流变模型在波动方程时间域正演模拟中的优势。
将广义流变模型引入到时间域粘弹性介质波动方程常Q正演模拟过程中,并分别针对广义Maxwell流变模型和广义标准线性体两种模型特征进行分析,对前者考虑了两种非弹性函数定义方式来建立对应的波动方程,分别采用弹性模量和松弛模量来描述基于广义标准线性体模型的波动方程,重新定义记忆变量使得通过两种广义流变模型所建立的波动方程表达式得到统一。完整的给出了基于两种模量表示的各种广义流变模型所对应的位移-应力、位移、位移-速度-应力和速度-应力方程,指出了早期波动方程由于采用松弛模量进行表示所引起的误解。基于五阶广义Maxwell流变模型建立了高精度粘弹性介质高阶交错网格有限差分波动方程常Q正演模拟表达式,采用PML边界条件有效的消除人工边界反射,讨论了模拟过程中的稳定性条件和震源问题,将有限差分正演模拟结果与二维均匀粘弹性介质中的解析解进行对比,分析了高精度粘弹性介质波动方程有限差分数值模拟方法的精度。
采用高精度粘弹性介质常Q波动方程正演模拟方法开展模型响应特征分析认识到地震波在粘弹性介质中具有较强的衰减频散特征,随着Q值的减小和传播距离的增大,地震波产生振幅衰减、相位改变、主频降低等现象。针对不同的观测系统开展正演模拟,通过波场、共炮集及偏移成像剖面特征进行对比分析可知,地面地震、VSP和井间地震三种多尺度地震资料之间既有联系也有差异,联系在于多尺度地震资料是相同地下地质体在相应观测方式下的地震波响应特征;差异在于资料所包含的波场信息丰富程度不一样、反映地质体信息的范围不一样、分辨地质体细节特征的能力也不一样。加入近地表强衰减层的影响表明,多尺度地震资料之间存在一定的相位差异,进一步认识了实际多尺度地震资料之间既具有一致性又存在一定差异的特征。
基于贝叶斯反演理论,通过似然函数有机的将地面地震、VSP和井间地震三种多尺度地震资料进行整合,建立起基于修正柯西分布的多尺度地震资料联合反演方法体系,并完善了联合反演流程。利用一维模型验证了联合反演算法的正确性,加噪后通过卡方分布测试了联合反演参数的合理选择范围,并采用二维多尺度正演模拟资料验证了联合反演方法的有效性及对多尺度资料的适应性。对联合反演方法进行退化建立了时间域单一尺度地震资料反演方法;针对连井井间地震剖面的特征重点强调了高精度深时转换、井地时深关系一致性及加强横向约束等改进策略,完善了深度域井间地震资料的反演方法和流程。
将多尺度地震资料反演方法对胜利油田垦71区块开展了应用研究,利用深度域井间地震反演方法获得连井井间地震波阻抗反演剖面,从该反演结果中提取的虚拟井能够有效的匹配地面地震资料反射特征。虚拟井约束下的井地联合反演只能够提高地面地震反演精度,而基于贝叶斯理论框架建立的多尺度地震资料联合反演方法则能够有效的提高地面地震资料反演结果的分辨率,实现多尺度地震资料特征的融合。
Reservoir geophysics is a new discipline that geophysics technique towards the integrated application of oilfield development and production, the reservoir description and reservoir management is included. The geophysics technique is adapted to determin the reservoir, descript the reservoir and monitor the change in the production. It is relying on the development of 3-D high-precision seismic technique, borehole seismic technique, time-lapse seismic technique and multi-wave/multi-component technique, the integrated research of high-precision 3-D seismic and other methods is the key technology of reservoir geophysics. Research the modeling and inversion method of multiscale seismic data is improving the capacity of seismic data to resolve the geology problem. Develop the high precision seismic modeling method of wave-equation in viscoelastic media, determine the mulitscale seismic response characteristics of geology with different observation system. Build up the multiscale seismic data joint inversion method based on the Bayes theory, overcome the multiple solution of inversion problem and improve the resolution of seismic data inversion, then service for the oilfield development with seismic data.
The nature attenuation law of seismic propagating in actual media is very complicated. many attenuation model of attenuation model is present in theory. The approximation to Q is obvious to most models in seismic wave-equation modeling at viscoelastic media in time domain. Based on the rheology theory, the adaptability of SLS (Standard Linear Solid), Kelvin and other attenuation model is researched in viscoelastic seismic modeling. The generalized rheology models that can be use to accurately describe the seismic wave propagating and attenuation characteristics have been analysied, the equivalence of generalize Maxwell model and generalized SLS model is proofed. The fitting degree to the constant Q of different order generalized Maxwell model is analysized and compared with the Futterman model and other models, the foundation is perfecing to modeling the constant Q characteristics in viscoelastic media. One dimension plane wave modeling is done with these attenuation models, the difference of seismic wave propagation and attenuation characteristic can be inteprated with these results, the advantage of generalized rheology model is clear.
The generalized rheology model is used in wave-equation seismic modeling at viscoelastic media in time domain. For the generalized Maxwell model, two kinds of unelastic function is adapted to derive the wave-equation. The unelastic modulus and relaxed modulus have been used to build up the wave-equation for generalized SLS model. The unified expression between these two kinds of generalized rheology model after the memory variable is redefined. The expression of displacement equation, the displacement-stress equation, the displacement-velocity-stress equation and the velocity-stress equation is presented according different modulus form.With the equation expansion or variable transform, the commom viscoelastic wave equation in articles can be get, the disadvantage of early equation expressed by relaxed modulus is clear. The high-order staggered grid finite difference scheme of high-precision viscoelastic wave equation based on five-order generalized Maxwell model is present. It is the high-precision wave-equation modeling method to constant Q in viscoelastic media. The PML boundary condition is used to absorb the artificial boundary reflection, the stability condition and source problem is discussed, the high-order finite difference scheme is accurate enough according to the comparison with the analytical resolution of 2D homogeneous viscoelastic media.
According the high-precision viscoelastic modeling, the attenuation and dispersion characterstic of seismic wave is obvious in viscoelastic media. With the Q decrease and the propagation distance increase, the amplitude is attenuated. the phase is lagged. the dominant frequency is decreased. Modeling with different observation system and compare with the wave field, common shot gather and migrated section characteristic, the relationship and difference among multiscale seismic data can be recognized. The relationship is multiscale seimic data is the seismic wave reflection response of the same geological body under different observation way. The difference is represented at the wave field information richness, the geological range and the detail characteristic in the multiscale seismic data. When the near-surface strong attenuation layer is considered, the phase difference is obvious between surface seismic data and others, the relationship and difference among multiscale seismic data can be inteprated.
Based on the Bayesian theorem, the new likelihood fundction is integrated with surface, vsp and crosswell seismic data, the multiscale seismic data joint inversion framework is built with modified Cauchy distribution and the inversion flow is improved. The one-dimension model testing shows that the method is right, the chi-square distribution is used to select the parameter range of the joint inversion with gauss noise, the joint inversion flow is valid to the two-dimension multiscale seismic data modeling according to the wave equation. For the singlescale seismic data, the crosswell seismic data inversion workflow in depth domain is improved. The high-precision interpolation function was adapted in depth-time transformation, the reflection characteristic consistence between surface seismic data and crosswell seismic data is emphasized, the transverse constraint is strengthened according to the characteristic of multi-crosswell seismic data.
The multiscale sesmic data inversion method is tested in Ken71 area of Shengli oilfield. The reasonable inversion result of multi-crosswell seismic data is get, the pseudo well exracted from it can mathch the characteristics of surface seismic data well. The constraint inversion with pseudo well can improve the precision of result. The multiscale seismic joint inversion method can reduce the multiplicity of solution, improve the inversion resolution, the multiscale data can be matched with each other according to the joint inversion result.
引文
[1] Pennington W D. Reservoir geophysics[J]. Geophysics. 2001, 66(1): 25~30.
[2]李阳.胜利油田油藏综合地球物理研究现状及展望[J].石油勘探与开发. 2004, 31(3).
[3]王喜双,甘利灯,易维启,等.油藏地球物理技术进展[J].石油地球物理勘探. 2006, 41(5): 606~613.
[4]杨文采.评地球物理反演的发展趋向[J].地学前缘. 2002, 9(4): 389~396.
[5]牟永光,裴正林.三维复杂介质地震数值模拟[M].北京:石油工业出版社, 2005.
[6] Carcione J M, Herman G C, Ten K A. Seismic modeling[J]. Geophysics. 2002, 67(4): 1304-1325.
[7] Alterman Z, Karal F C. Propagation of elastic waves in layered media by finite-difference methods[J]. Bulletin of the Seismological Society of America. 1968, 58(1): 367~398.
[8] Alford R M, Kelly K R, Boore D M. accuracy of finite-difference modeling of the acoustic wave equation[J]. Geophysics. 1974, 39(4): 834～842.
[9] Kelly K R, Ward R W, Treitel S, et al. Synthetic seismograms: a finite-difference approach[J]. 1976, 41(1).
[10]刘洋.裂缝性油气藏多波多分量地震勘探方法研究[D].北京:石油大学, 1998.
[11] Madariaga R A. Dynamics of an expanding circular fault[J]. Bulletin of the Seismological Society of America. 1976, 66(3): 639~666.
[12] Virieux J. P-SV wave propagation in heterogeneous media: Velocity-stress finite-difference method[J]. Geophysics. 1986, 51(4): 889~901.
[13] Virieux J. SH-wave propagation in heterogeneous media: Velocity-stress finite-difference method[J]. Geophysics. 1984, 49(11): 1933~1942.
[14] Levander A R. Fourth-order finite-difference P-SV seismograms[J]. Geophysics. 1988, 53(11): 1425~1436.
[15] Bayliss A, Jordan K E, Lemesurier B J, et al. A fourth-order accurate finite-difference scheme for the computation of elastic waves[J]. Bulletin of the Seismological Society ofAmerica. 1986, 76(4): 1115~1132.
[16] Dablain M A. The application of high-order differencing to the scalar wave equation[J]. 1986, 51(1): 54~66.
[17]董良国,马在田,曹景忠,等.一阶弹性波方程交错网格高阶差分解法[J].地球物理学报. 2000, 43(3): 411~419.
[18] Graves R W. Simulating seismic wave propagation in 3D elastic media using staggered-grid finite differences[J]. Bulletin of the Seismological Society of America. 1996, 86(4): 1091~1106.
[19] Reynolds A C. Boundary conditions for the numerical solution of wave propagation problems[J]. Geophysics. 1978, 43(6): 1099~1110.
[20] Clayton R O, Engquist B J. Absorbing boundary conditions for acoustic and elastic wave equations[J]. Bulletin of the Seismological Society of America. 1977, 67(6): 1529~1540.
[21] Cerjan C, Kosloff D, Kosloff R, et al. A nonreflecting boundary condition for discrete acoustic and elastic wave equations[J]. Geophysics. 1985, 50(4): 705~708.
[22] Berenger J P. A perfectly matched layer for the absorption of electromagnetic waves[J]. J. Comput. Phys. 1994, 114(2): 185~200.
[23] Collino F, Tsogka C. Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media[J]. Geophysics. 2001, 66(1): 294~307.
[24] Festa G, Nielsen S. PML Absorbing Boundaries[J]. Bulletin of the Seismological Society of America. 2003, 93(2): 891~903.
[25]蔡其新,何佩军,秦广胜,等.有限差分数值模拟的最小频散算法及其应用[J].石油地球物理勘探. 2003, 38(3): 247~251,262.
[26]董良国,李培明.地震波传播数值模拟中的频散问题[J].天然气工业. 2004, 24(6): 53~56.
[27]吴国忱,王华忠.波场模拟中的数值频散分析与校正策略[J].地球物理学进展. 2005, 20(1): 58~65.
[28] Ricker N.粘弹性介质中的地震波[M].北京:地质出版社, 1981.
[29] Born W T. The attenuation constant of earth materials[J]. Geophysics. 1941, 6(2): 132-148.
[30] Mcdonal F J, Angona F A, Mills R L, et al. Attenuation of shear and compressional waves in Pierre shale[J]. Geophysics. 1958, 23(3): 421~439.
[31] Wuenschel P C. Dispersive body waves - an experimental study[J]. Geophysics. 1965, 30(4): 539-551.
[32] Liu H, Anderson D L, Kanamori H. Velocity dispersion due to anelasticity; implications for seismology and mantle composition[J]. Geophysical Journal International. 1976, 47(1): 41~58.
[33] Futterman W I. Dispersive body waves[J]. Journal of geophysical research. 1962, 67(13): 5279-5291.
[34] Kjartansson E. Attenuation of seismic waves in rocks and applications in energy exploration[D]. Stanford Univeristy, 1980.
[35] Aki K, Richards P G. Quantitative seismology : theory and methods[M]. Freeman, 1980.
[36] Christensen R M. Theory of Viscoelasticity: An Introduction[M]. New York: Academic Press, 1982.
[37] Day S M, Minster J B. Numerical simulation of attenuated wavefields using a Pade approximant method[J]. Geophysical Journal of the Royal Astronomical Society. 1984, 78(1): 105-118.
[38] Emmerich H, Korn M. Incorporation of attenuation into time-domain computations of seismic wave fields[J]. Geophysics. 1987, 52(9): 1252~1264.
[39] Carcione J M, Kosloff D, Kosloff R. Viscoacoustic wave propagation simulation in the earth[J]. Geophysics. 1988, 53(6): 769~777.
[40] Carcione J M. Wave propagation simulation in a linear viscoelastic medium[J]. Geophysical Journal International. 1988, 95(3): 597~611.
[41] Carcione J M. Seismic modeling in viscoelastic media[J]. Geophysics. 1993, 58(1): 110~120.
[42] Tal-ezer A H, Carcione J M, Kosloff D. An accurate and efficient scheme for wavepropagation in linear viscoelastic media[J]. Geophysics. 1990, 55(10): 1366-1379.
[43] Robertsson J O, Blanch J O, Symes W W. Viscoelastic finite-difference modeling[J]. Geophysics. 1994, 59(9): 1444~1456.
[44] Blanch J O, Robertsson J O, Symes W W. Modeling of a constant Q: Methodology and algorithm for an efficient and optimally inexpensive viscoelastic technique[J]. Geophysics. 1995, 60(1): 176~184.
[45] Xu T, Mc G A. Efficient 3-D viscoelastic modeling with application to near-surface land seismic data[J]. Geophysics. 1998, 63(2): 601~612.
[46] Stekl I, Pratt R G. Accurate viscoelastic modeling by frequency-domain finite differences using rotated operators[J]. Geophysics. 1998, 63(5): 1779~1794.
[47] Bohlen T. Parallel 3-D viscoelastic finite difference seismic modelling[J]. Computers & Geosciences. 2002, 28(8): 887-899.
[48] Russo G, Zollo A. A constant Q technique for the numerical simulation of attenuation of seismic body waves[J]. Geophysics. 2003, 68(5): 1744~1748.
[49] Kristek J, Moczo P. Seismic-Wave Propagation in Viscoelastic Media with Material Discontinuities: A 3D Fourth-Order Staggered-Grid Finite-Difference Modeling[J]. Bulletin of the Seismological Society of America. 2003, 93(5): 2273~2280.
[50] Moczo P, Kristek J. On the rheological models used for time-domain methods of seismic wave propagation[J]. Geophysical Research Letter. 2005, 32: L01306.
[51]范祯祥,郑仙种.地震波数值模拟与偏移成像[M].郑州:河南科学技术出版社, 1994.
[52]毕玉英,杨宝俊.二维粘弹介质中完全波场正演模拟[J].石油地球物理勘探. 1995, 30(3): 351~362.
[53]崔建军,何继善.粘弹性波动方程正演和偏移[J].中南工业大学学报(自然科学版). 2001, 32(5): 441~444.
[54]杜启振,杨慧珠.线性黏弹性各向异性介质速度频散和衰减特征研究[J].物理学报. 2002, 51(9): 2101~2108.
[55]杜启振,刘莲莲,孙晶波.各向异性粘弹性孔隙介质地震波场伪谱法正演模拟[J].物理学报. 2007, 56(10): 6143~6149.
[56]苑春方. Kelvin-Voigt均匀黏弹性介质中传播的地震波[J].中国科学D辑. 2005, 35(10): 957~962.
[57]单启铜.粘弹性波动方程正演模拟与参数反演[D].东营:中国石油大学(华东), 2007.
[58]单启铜,乐友喜. PML边界条件下二维粘弹性介质波场模拟[J]. 2007, 46(2): 126~130,206.
[59]王德利,雍运动,韩立国,等.三维粘弹介质地震波场有限差分并行模拟[J].西北地震学报. 2007, 29(1): 30~34.
[60] Day S M, Bradley C R. Memory-Efficient Simulation of Anelastic Wave Propagation[J]. Bulletin of the Seismological Society of America. 2001, 91(3): 520-531.
[61] Day S M. Efficient simulation of constant Q using coarse-grained memory variables[J]. Bulletin of the Seismological Society of America. 1998, 88(4): 1051-1062.
[62] Liu P, Archuleta R J. Efficient Modeling of Q for 3D Numerical Simulation of Wave Propagation[J]. Bulletin of the Seismological Society of America. 2006, 96(4A): 1352~1358.
[63] Moczo P, Bystricky E, Kristek J, et al. Hybrid modeling of P-SV seismic motion at inhomogeneous viscoelastic topographic structures[J]. Bulletin of the Seismological Society of America. 1997, 87(5): 1305~1323.
[64] Krebes E S, Goode G Q. A standard finite-difference scheme for the time-domain computation of anelastic wavefields[J]. Geophysics. 1994, 59(2): 290~296.
[65] Emmerich H. PSV-wave propagation in a medium with local heterogeneities: a hybrid formulation and its application[J]. Geophysical Journal International. 1992, 109(1): 54~64.
[66] Moczo P, Lucka M, Kristek J, et al. 3D Displacement finite differences and a combined memory optimization[J]. Bulletin of the Seismological Society of America. 1999, 89(1): 69-79.
[67] Hestholm S, Ruud B. 3D free-boundary conditions for coordinate-transform finite-difference seismic modelling.[J]. Geophysical Prospecting. 2002, 50(5): 463~474.
[68] Hestholm S. Composite memory variable velocity-stress viscoelastic modelling.[J]. Geophysical Journal International. 2002, 148(1): 153~162.
[69] Ruud B, Hestholm S. 2D surface topography boundary conditions in seismic wave modelling.[J]. Geophysical Prospecting. 2001, 49(4): 445~460.
[70] Hestholm S O, Ruud B. 2D finite-difference viscoelastic wave modelling including surface topography.[J]. Geophysical Prospecting. 2000, 48(2): 341~373.
[71] Hestholm S. Three-dimensional finite difference viscoelastic wave modelling including surface topography.[J]. Geophysical Journal International. 1999, 139(3): 852~878.
[72] Carcione J M. Wave fields in real media: Wave propagation in anisotropic, anelastic and porous media[M]. Pergamon Press, 2001.
[73]郭建. VSP技术应用现状及发展趋势[J].勘探地球物理进展. 2004, 27(1): 1~8.
[74] Wyatt K D. Synthetic vertical seismic profile[J]. Geophysics. 1981, 46(6): 880~891.
[75] Young T K. Computer modeling of vertical seismic profiling [J]. Geophysical Prospecting. 1984, 32(5): 851~870.
[76] Mcmechan G A. Synthetic finite-offset vertical seismic profiles for laterally varying media[J]. Geophysics. 1985, 50(4): 627~636.
[77]肖春燕,庄东海,金耀祖.垂直地震剖面的有限元正演模拟[J].江汉石油学院学报. 1994, 16(S1): 81~86.
[78]吴清岭,张平,马振轶.弹性波方程的VSP合成记录[J].大庆石油地质与开发. 1998, 17(4): 37~38,40.
[79]涂齐催. VSP射线追踪模拟与处理方法研究及应用[D].中国海洋大学, 2006.
[80]王成礼,韩文功,王延光,等. VSP弹性波方程正演模拟与波场分析[J].油气地球物理. 2006, 4(4): 7~14.
[81]崔杰,韩立国,单刚义.含衰减的零偏移距VSP正演模拟及Q值的提取[J].吉林大学学报(地球科学版). 2006, 36(S2): 125~129.
[82]曹辉.井间地震技术发展现状[J].勘探地球物理进展. 2002, 25(6): 6~10.
[83]陈世军,刘洪,周建宇,等.井间地震技术的现状与展望[J].地球物理学进展. 2003, 18(3): 524~529.
[84]周建宇.井间地震研究与应用[D].北京:中国科学院研究生院(兰州地质研究所), 2002.
[85]周建宇,李建平,张永辉.井间地震应用效果分析[J].地球物理学进展. 2003, 18(3): 567~575.
[86]曹辉,郭全仕,唐金良,等.井间地震资料特点分析[J].勘探地球物理进展. 2006, 29(5): 312~317.
[87]杜世通.用有限单元法作井间地震观测数据的模拟和偏移[J].石油地球物理勘探. 1990, 25(3): 27~36,138.
[88]缪林昌.井间地震Hartley变换法正演模拟与偏移[J].石油物探. 1994, 33(4): 70~77.
[89]吴律.层析基础及其在井间地震中的应用[M].北京:石油工业出版社, 1997.
[90]张汝杰,贺振华,王理.井间地震高斯射线束正演方法[J].物探化探计算技术. 1997, 19(2): 128~137.
[91]赵连锋.井间地震波速与衰减联合层析成像方法研究[D].成都:成都理工大学, 2002.
[92]李庆忠,王建花.井间地震勘探的误区及出路[J].石油地球物理勘探. 2004, 39(1): 1~11.
[93]马德堂.弹性波场数值模拟及井间地震初至波旅行时层析成像[D].长安大学, 2005.
[94]张文波.井间地震交错网格高阶差分数值模拟及逆时偏移成像研究[D].长安大学, 2005.
[95]窦玉坛,郭常升,毛中华,等.井间地震波动方程正演数值模拟及其应用[J].地球物理学进展. 2007, 22(6): 1867~1872.
[96]窦玉坛.基于正演方法的井间地震响应特征研究[D].北京:中国科学院研究生院(海洋研究所), 2007.
[97]宋常瑜,裴正林.井间地震粘弹性波场特征的数值模拟研究[J].石油物探. 2006, 45(5): 508~513.
[98]李万万,裴正林.井间地震弹性波传播特征数值模拟[J].物探与化探. 2008, 32(2): 207~211.
[99]左建军.井间地震粘弹性数值模拟及效果分析[J].油气地球物理. 2008, 6(1): 25~27.
[100] Cooke D A, Schneider W A. Generalized linear inversion of reflection seismic data[J]. Geophysics. 1983, 48(6): 665~676.
[101]张永刚.地震波阻抗反演技术的现状和发展[J].石油物探. 2002, 41(4): 385~390.
[102]何惺华.深度域地震资料若干问题初探[J].石油物探. 2004, 43(4): 353~358.
[103]吕铁良,王永刚,谢万学,等.稀疏脉冲反演技术在井间地震反演中的应用[J].石油物探. 2007, 46(1): 58~63.
[104]孟宪军,王延光,孙振涛,等.井间地震资料时间域波阻抗反演研究[J].石油学报. 2005, 26(1): 47~49,54.
[105]许丽,刘磊,吴长双.应用井间地震资料建立砂泥岩薄互层地质模型[J].石油地球物理勘探. 2006, 41(1): 58~61.
[106]沈财余,江洁,赵华,等.测井约束地震反演解决地质问题能力的探讨[J].石油地球物理勘探. 2002, 37(4): 372-376.
[107]李庆忠.论地震约束反演的策略[J].石油地球物理勘探. 1998, 33(4): 423~438.
[108] Alekseev A S. Quantitative statement and general properties of solutions of cooperative inverse problems (integral geophysics)[J]. SEG Technical Program Expanded Abstracts. 1991, 10(1): 1021~1022.
[109] Vozoff K, Jupp D L B. Effective search for a buried layer: An approach to experimental design in geophysics[J]. Exploration. Geophysics. 1977, 8(1): 6~15.
[110]敬荣中.地球物理非线性联合反演方法研究[D].长沙:中南大学, 2002.
[111]李爽.大地测量联合反演的模式及算法研究[D].武汉:武汉大学, 2005.
[112]杨辉,戴世坤,宋海斌,等.综合地球物理联合反演综述[J].地球物理学进展. 2002, 17(2): 262-271.
[113]敬荣中,鲍光淑,陈绍裘.地球物理联合反演研究综述[J].地球物理学进展. 2003, 18(3): 535-540.
[114]于鹏,王家林,吴健生,等.重力与地震资料的模拟退火约束联合反演[J].地球物理学报. 2007, 50(2): 529-538.
[115]陈洁,温宁,陈邦彦.重磁电震联合反演研究进展与展望[J].地球物理学进展. 2007, 22(5): 1427-1438.
[116] Bosch M, Mcgaughey J. Joint inversion of gravity and magnetic data under lithologic constraints[J]. The Leading Edge. 2001, 20(8): 877~881.
[117] Laurence R L. Cooperative inversion of geophysical data[J]. Geophysics. 1988, 53(1): 8-20.
[118] Colombo D, De S M. Geophysical modeling via simultaneous joint inversion of seismic, gravity, and electromagnetic data: Application to prestack depth imaging[J]. The Leading Edge. 2007, 26(3): 326~331.
[119] Hu W, Abubakar A, Habashy T. Joint inversion algorithm for electromagnetic and seismic data[J]. SEG Technical Program Expanded Abstracts. 2007, 26(1): 1745~1749.
[120] Mackie R, Watts M D, Rodi W. Joint 3D inversion of marine CSEM and MT data[J]. SEG Technical Program Expanded Abstracts. 2007, 26(1): 574~578.
[121] Heincke B, Jegen M. Hobbs R. Joint inversion of MT, gravity and seismic data applied to sub~basalt imaging[J]. SEG Technical Program Expanded Abstracts. 2006, 25(1): 784~789.
[122] Maxwell A M. Joint inversion of TEM and distorted MT soundings; some effective practical considerations[J]. Geophysics. 1996, 61(1): 56~65.
[123] Helene H V, Martin L. Simultaneous inversion of PP and PS seismic data[J]. Geophysics. 2006, 71(3).
[124] Fu L Y. Joint inversion of seismic data for acoustic impedance[J]. Geophysics. 2004, 69(4): 994~1004.
[125] Tarantola A. Inverse Problem Theory and Methods for Model Parameter Estimation[M]. philadephia: Society for Industrial and Applied Mathematics, 2005.
[126] Ulrych T J, Sacchi M D, Woodbury A. A Bayes tour of inversion: A tutorial[J]. Geophysics. 2001, 66(1): 55-69.
[127] Youzwishen C F. Non-linear sparse and blocky constraint for seismic inverse problems[D]. University of Alberta, 2001.
[2]李阳.胜利油田油藏综合地球物理研究现状及展望[J].石油勘探与开发. 2004, 31(3).
[3]王喜双,甘利灯,易维启,等.油藏地球物理技术进展[J].石油地球物理勘探. 2006, 41(5): 606~613.
[4]杨文采.评地球物理反演的发展趋向[J].地学前缘. 2002, 9(4): 389~396.
[5]牟永光,裴正林.三维复杂介质地震数值模拟[M].北京:石油工业出版社, 2005.
[6] Carcione J M, Herman G C, Ten K A. Seismic modeling[J]. Geophysics. 2002, 67(4): 1304-1325.
[7] Alterman Z, Karal F C. Propagation of elastic waves in layered media by finite-difference methods[J]. Bulletin of the Seismological Society of America. 1968, 58(1): 367~398.
[8] Alford R M, Kelly K R, Boore D M. accuracy of finite-difference modeling of the acoustic wave equation[J]. Geophysics. 1974, 39(4): 834～842.
[9] Kelly K R, Ward R W, Treitel S, et al. Synthetic seismograms: a finite-difference approach[J]. 1976, 41(1).
[10]刘洋.裂缝性油气藏多波多分量地震勘探方法研究[D].北京:石油大学, 1998.
[11] Madariaga R A. Dynamics of an expanding circular fault[J]. Bulletin of the Seismological Society of America. 1976, 66(3): 639~666.
[12] Virieux J. P-SV wave propagation in heterogeneous media: Velocity-stress finite-difference method[J]. Geophysics. 1986, 51(4): 889~901.
[13] Virieux J. SH-wave propagation in heterogeneous media: Velocity-stress finite-difference method[J]. Geophysics. 1984, 49(11): 1933~1942.
[14] Levander A R. Fourth-order finite-difference P-SV seismograms[J]. Geophysics. 1988, 53(11): 1425~1436.
[15] Bayliss A, Jordan K E, Lemesurier B J, et al. A fourth-order accurate finite-difference scheme for the computation of elastic waves[J]. Bulletin of the Seismological Society ofAmerica. 1986, 76(4): 1115~1132.
[16] Dablain M A. The application of high-order differencing to the scalar wave equation[J]. 1986, 51(1): 54~66.
[17]董良国,马在田,曹景忠,等.一阶弹性波方程交错网格高阶差分解法[J].地球物理学报. 2000, 43(3): 411~419.
[18] Graves R W. Simulating seismic wave propagation in 3D elastic media using staggered-grid finite differences[J]. Bulletin of the Seismological Society of America. 1996, 86(4): 1091~1106.
[19] Reynolds A C. Boundary conditions for the numerical solution of wave propagation problems[J]. Geophysics. 1978, 43(6): 1099~1110.
[20] Clayton R O, Engquist B J. Absorbing boundary conditions for acoustic and elastic wave equations[J]. Bulletin of the Seismological Society of America. 1977, 67(6): 1529~1540.
[21] Cerjan C, Kosloff D, Kosloff R, et al. A nonreflecting boundary condition for discrete acoustic and elastic wave equations[J]. Geophysics. 1985, 50(4): 705~708.
[22] Berenger J P. A perfectly matched layer for the absorption of electromagnetic waves[J]. J. Comput. Phys. 1994, 114(2): 185~200.
[23] Collino F, Tsogka C. Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media[J]. Geophysics. 2001, 66(1): 294~307.
[24] Festa G, Nielsen S. PML Absorbing Boundaries[J]. Bulletin of the Seismological Society of America. 2003, 93(2): 891~903.
[25]蔡其新,何佩军,秦广胜,等.有限差分数值模拟的最小频散算法及其应用[J].石油地球物理勘探. 2003, 38(3): 247~251,262.
[26]董良国,李培明.地震波传播数值模拟中的频散问题[J].天然气工业. 2004, 24(6): 53~56.
[27]吴国忱,王华忠.波场模拟中的数值频散分析与校正策略[J].地球物理学进展. 2005, 20(1): 58~65.
[28] Ricker N.粘弹性介质中的地震波[M].北京:地质出版社, 1981.
[29] Born W T. The attenuation constant of earth materials[J]. Geophysics. 1941, 6(2): 132-148.
[30] Mcdonal F J, Angona F A, Mills R L, et al. Attenuation of shear and compressional waves in Pierre shale[J]. Geophysics. 1958, 23(3): 421~439.
[31] Wuenschel P C. Dispersive body waves - an experimental study[J]. Geophysics. 1965, 30(4): 539-551.
[32] Liu H, Anderson D L, Kanamori H. Velocity dispersion due to anelasticity; implications for seismology and mantle composition[J]. Geophysical Journal International. 1976, 47(1): 41~58.
[33] Futterman W I. Dispersive body waves[J]. Journal of geophysical research. 1962, 67(13): 5279-5291.
[34] Kjartansson E. Attenuation of seismic waves in rocks and applications in energy exploration[D]. Stanford Univeristy, 1980.
[35] Aki K, Richards P G. Quantitative seismology : theory and methods[M]. Freeman, 1980.
[36] Christensen R M. Theory of Viscoelasticity: An Introduction[M]. New York: Academic Press, 1982.
[37] Day S M, Minster J B. Numerical simulation of attenuated wavefields using a Pade approximant method[J]. Geophysical Journal of the Royal Astronomical Society. 1984, 78(1): 105-118.
[38] Emmerich H, Korn M. Incorporation of attenuation into time-domain computations of seismic wave fields[J]. Geophysics. 1987, 52(9): 1252~1264.
[39] Carcione J M, Kosloff D, Kosloff R. Viscoacoustic wave propagation simulation in the earth[J]. Geophysics. 1988, 53(6): 769~777.
[40] Carcione J M. Wave propagation simulation in a linear viscoelastic medium[J]. Geophysical Journal International. 1988, 95(3): 597~611.
[41] Carcione J M. Seismic modeling in viscoelastic media[J]. Geophysics. 1993, 58(1): 110~120.
[42] Tal-ezer A H, Carcione J M, Kosloff D. An accurate and efficient scheme for wavepropagation in linear viscoelastic media[J]. Geophysics. 1990, 55(10): 1366-1379.
[43] Robertsson J O, Blanch J O, Symes W W. Viscoelastic finite-difference modeling[J]. Geophysics. 1994, 59(9): 1444~1456.
[44] Blanch J O, Robertsson J O, Symes W W. Modeling of a constant Q: Methodology and algorithm for an efficient and optimally inexpensive viscoelastic technique[J]. Geophysics. 1995, 60(1): 176~184.
[45] Xu T, Mc G A. Efficient 3-D viscoelastic modeling with application to near-surface land seismic data[J]. Geophysics. 1998, 63(2): 601~612.
[46] Stekl I, Pratt R G. Accurate viscoelastic modeling by frequency-domain finite differences using rotated operators[J]. Geophysics. 1998, 63(5): 1779~1794.
[47] Bohlen T. Parallel 3-D viscoelastic finite difference seismic modelling[J]. Computers & Geosciences. 2002, 28(8): 887-899.
[48] Russo G, Zollo A. A constant Q technique for the numerical simulation of attenuation of seismic body waves[J]. Geophysics. 2003, 68(5): 1744~1748.
[49] Kristek J, Moczo P. Seismic-Wave Propagation in Viscoelastic Media with Material Discontinuities: A 3D Fourth-Order Staggered-Grid Finite-Difference Modeling[J]. Bulletin of the Seismological Society of America. 2003, 93(5): 2273~2280.
[50] Moczo P, Kristek J. On the rheological models used for time-domain methods of seismic wave propagation[J]. Geophysical Research Letter. 2005, 32: L01306.
[51]范祯祥,郑仙种.地震波数值模拟与偏移成像[M].郑州:河南科学技术出版社, 1994.
[52]毕玉英,杨宝俊.二维粘弹介质中完全波场正演模拟[J].石油地球物理勘探. 1995, 30(3): 351~362.
[53]崔建军,何继善.粘弹性波动方程正演和偏移[J].中南工业大学学报(自然科学版). 2001, 32(5): 441~444.
[54]杜启振,杨慧珠.线性黏弹性各向异性介质速度频散和衰减特征研究[J].物理学报. 2002, 51(9): 2101~2108.
[55]杜启振,刘莲莲,孙晶波.各向异性粘弹性孔隙介质地震波场伪谱法正演模拟[J].物理学报. 2007, 56(10): 6143~6149.
[56]苑春方. Kelvin-Voigt均匀黏弹性介质中传播的地震波[J].中国科学D辑. 2005, 35(10): 957~962.
[57]单启铜.粘弹性波动方程正演模拟与参数反演[D].东营:中国石油大学(华东), 2007.
[58]单启铜,乐友喜. PML边界条件下二维粘弹性介质波场模拟[J]. 2007, 46(2): 126~130,206.
[59]王德利,雍运动,韩立国,等.三维粘弹介质地震波场有限差分并行模拟[J].西北地震学报. 2007, 29(1): 30~34.
[60] Day S M, Bradley C R. Memory-Efficient Simulation of Anelastic Wave Propagation[J]. Bulletin of the Seismological Society of America. 2001, 91(3): 520-531.
[61] Day S M. Efficient simulation of constant Q using coarse-grained memory variables[J]. Bulletin of the Seismological Society of America. 1998, 88(4): 1051-1062.
[62] Liu P, Archuleta R J. Efficient Modeling of Q for 3D Numerical Simulation of Wave Propagation[J]. Bulletin of the Seismological Society of America. 2006, 96(4A): 1352~1358.
[63] Moczo P, Bystricky E, Kristek J, et al. Hybrid modeling of P-SV seismic motion at inhomogeneous viscoelastic topographic structures[J]. Bulletin of the Seismological Society of America. 1997, 87(5): 1305~1323.
[64] Krebes E S, Goode G Q. A standard finite-difference scheme for the time-domain computation of anelastic wavefields[J]. Geophysics. 1994, 59(2): 290~296.
[65] Emmerich H. PSV-wave propagation in a medium with local heterogeneities: a hybrid formulation and its application[J]. Geophysical Journal International. 1992, 109(1): 54~64.
[66] Moczo P, Lucka M, Kristek J, et al. 3D Displacement finite differences and a combined memory optimization[J]. Bulletin of the Seismological Society of America. 1999, 89(1): 69-79.
[67] Hestholm S, Ruud B. 3D free-boundary conditions for coordinate-transform finite-difference seismic modelling.[J]. Geophysical Prospecting. 2002, 50(5): 463~474.
[68] Hestholm S. Composite memory variable velocity-stress viscoelastic modelling.[J]. Geophysical Journal International. 2002, 148(1): 153~162.
[69] Ruud B, Hestholm S. 2D surface topography boundary conditions in seismic wave modelling.[J]. Geophysical Prospecting. 2001, 49(4): 445~460.
[70] Hestholm S O, Ruud B. 2D finite-difference viscoelastic wave modelling including surface topography.[J]. Geophysical Prospecting. 2000, 48(2): 341~373.
[71] Hestholm S. Three-dimensional finite difference viscoelastic wave modelling including surface topography.[J]. Geophysical Journal International. 1999, 139(3): 852~878.
[72] Carcione J M. Wave fields in real media: Wave propagation in anisotropic, anelastic and porous media[M]. Pergamon Press, 2001.
[73]郭建. VSP技术应用现状及发展趋势[J].勘探地球物理进展. 2004, 27(1): 1~8.
[74] Wyatt K D. Synthetic vertical seismic profile[J]. Geophysics. 1981, 46(6): 880~891.
[75] Young T K. Computer modeling of vertical seismic profiling [J]. Geophysical Prospecting. 1984, 32(5): 851~870.
[76] Mcmechan G A. Synthetic finite-offset vertical seismic profiles for laterally varying media[J]. Geophysics. 1985, 50(4): 627~636.
[77]肖春燕,庄东海,金耀祖.垂直地震剖面的有限元正演模拟[J].江汉石油学院学报. 1994, 16(S1): 81~86.
[78]吴清岭,张平,马振轶.弹性波方程的VSP合成记录[J].大庆石油地质与开发. 1998, 17(4): 37~38,40.
[79]涂齐催. VSP射线追踪模拟与处理方法研究及应用[D].中国海洋大学, 2006.
[80]王成礼,韩文功,王延光,等. VSP弹性波方程正演模拟与波场分析[J].油气地球物理. 2006, 4(4): 7~14.
[81]崔杰,韩立国,单刚义.含衰减的零偏移距VSP正演模拟及Q值的提取[J].吉林大学学报(地球科学版). 2006, 36(S2): 125~129.
[82]曹辉.井间地震技术发展现状[J].勘探地球物理进展. 2002, 25(6): 6~10.
[83]陈世军,刘洪,周建宇,等.井间地震技术的现状与展望[J].地球物理学进展. 2003, 18(3): 524~529.
[84]周建宇.井间地震研究与应用[D].北京:中国科学院研究生院(兰州地质研究所), 2002.
[85]周建宇,李建平,张永辉.井间地震应用效果分析[J].地球物理学进展. 2003, 18(3): 567~575.
[86]曹辉,郭全仕,唐金良,等.井间地震资料特点分析[J].勘探地球物理进展. 2006, 29(5): 312~317.
[87]杜世通.用有限单元法作井间地震观测数据的模拟和偏移[J].石油地球物理勘探. 1990, 25(3): 27~36,138.
[88]缪林昌.井间地震Hartley变换法正演模拟与偏移[J].石油物探. 1994, 33(4): 70~77.
[89]吴律.层析基础及其在井间地震中的应用[M].北京:石油工业出版社, 1997.
[90]张汝杰,贺振华,王理.井间地震高斯射线束正演方法[J].物探化探计算技术. 1997, 19(2): 128~137.
[91]赵连锋.井间地震波速与衰减联合层析成像方法研究[D].成都:成都理工大学, 2002.
[92]李庆忠,王建花.井间地震勘探的误区及出路[J].石油地球物理勘探. 2004, 39(1): 1~11.
[93]马德堂.弹性波场数值模拟及井间地震初至波旅行时层析成像[D].长安大学, 2005.
[94]张文波.井间地震交错网格高阶差分数值模拟及逆时偏移成像研究[D].长安大学, 2005.
[95]窦玉坛,郭常升,毛中华,等.井间地震波动方程正演数值模拟及其应用[J].地球物理学进展. 2007, 22(6): 1867~1872.
[96]窦玉坛.基于正演方法的井间地震响应特征研究[D].北京:中国科学院研究生院(海洋研究所), 2007.
[97]宋常瑜,裴正林.井间地震粘弹性波场特征的数值模拟研究[J].石油物探. 2006, 45(5): 508~513.
[98]李万万,裴正林.井间地震弹性波传播特征数值模拟[J].物探与化探. 2008, 32(2): 207~211.
[99]左建军.井间地震粘弹性数值模拟及效果分析[J].油气地球物理. 2008, 6(1): 25~27.
[100] Cooke D A, Schneider W A. Generalized linear inversion of reflection seismic data[J]. Geophysics. 1983, 48(6): 665~676.
[101]张永刚.地震波阻抗反演技术的现状和发展[J].石油物探. 2002, 41(4): 385~390.
[102]何惺华.深度域地震资料若干问题初探[J].石油物探. 2004, 43(4): 353~358.
[103]吕铁良,王永刚,谢万学,等.稀疏脉冲反演技术在井间地震反演中的应用[J].石油物探. 2007, 46(1): 58~63.
[104]孟宪军,王延光,孙振涛,等.井间地震资料时间域波阻抗反演研究[J].石油学报. 2005, 26(1): 47~49,54.
[105]许丽,刘磊,吴长双.应用井间地震资料建立砂泥岩薄互层地质模型[J].石油地球物理勘探. 2006, 41(1): 58~61.
[106]沈财余,江洁,赵华,等.测井约束地震反演解决地质问题能力的探讨[J].石油地球物理勘探. 2002, 37(4): 372-376.
[107]李庆忠.论地震约束反演的策略[J].石油地球物理勘探. 1998, 33(4): 423~438.
[108] Alekseev A S. Quantitative statement and general properties of solutions of cooperative inverse problems (integral geophysics)[J]. SEG Technical Program Expanded Abstracts. 1991, 10(1): 1021~1022.
[109] Vozoff K, Jupp D L B. Effective search for a buried layer: An approach to experimental design in geophysics[J]. Exploration. Geophysics. 1977, 8(1): 6~15.
[110]敬荣中.地球物理非线性联合反演方法研究[D].长沙:中南大学, 2002.
[111]李爽.大地测量联合反演的模式及算法研究[D].武汉:武汉大学, 2005.
[112]杨辉,戴世坤,宋海斌,等.综合地球物理联合反演综述[J].地球物理学进展. 2002, 17(2): 262-271.
[113]敬荣中,鲍光淑,陈绍裘.地球物理联合反演研究综述[J].地球物理学进展. 2003, 18(3): 535-540.
[114]于鹏,王家林,吴健生,等.重力与地震资料的模拟退火约束联合反演[J].地球物理学报. 2007, 50(2): 529-538.
[115]陈洁,温宁,陈邦彦.重磁电震联合反演研究进展与展望[J].地球物理学进展. 2007, 22(5): 1427-1438.
[116] Bosch M, Mcgaughey J. Joint inversion of gravity and magnetic data under lithologic constraints[J]. The Leading Edge. 2001, 20(8): 877~881.
[117] Laurence R L. Cooperative inversion of geophysical data[J]. Geophysics. 1988, 53(1): 8-20.
[118] Colombo D, De S M. Geophysical modeling via simultaneous joint inversion of seismic, gravity, and electromagnetic data: Application to prestack depth imaging[J]. The Leading Edge. 2007, 26(3): 326~331.
[119] Hu W, Abubakar A, Habashy T. Joint inversion algorithm for electromagnetic and seismic data[J]. SEG Technical Program Expanded Abstracts. 2007, 26(1): 1745~1749.
[120] Mackie R, Watts M D, Rodi W. Joint 3D inversion of marine CSEM and MT data[J]. SEG Technical Program Expanded Abstracts. 2007, 26(1): 574~578.
[121] Heincke B, Jegen M. Hobbs R. Joint inversion of MT, gravity and seismic data applied to sub~basalt imaging[J]. SEG Technical Program Expanded Abstracts. 2006, 25(1): 784~789.
[122] Maxwell A M. Joint inversion of TEM and distorted MT soundings; some effective practical considerations[J]. Geophysics. 1996, 61(1): 56~65.
[123] Helene H V, Martin L. Simultaneous inversion of PP and PS seismic data[J]. Geophysics. 2006, 71(3).
[124] Fu L Y. Joint inversion of seismic data for acoustic impedance[J]. Geophysics. 2004, 69(4): 994~1004.
[125] Tarantola A. Inverse Problem Theory and Methods for Model Parameter Estimation[M]. philadephia: Society for Industrial and Applied Mathematics, 2005.
[126] Ulrych T J, Sacchi M D, Woodbury A. A Bayes tour of inversion: A tutorial[J]. Geophysics. 2001, 66(1): 55-69.
[127] Youzwishen C F. Non-linear sparse and blocky constraint for seismic inverse problems[D]. University of Alberta, 2001.
© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号 地址：北京市海淀区学院路29号 邮编：100083 电话：办公室：(+86 10)66554848；文献借阅、咨询服务、科技查新：66554700 |