裂缝性孔隙介质波场模拟与频变AVO储层参数反演
|
摘要
随着油气勘探开发程度的不断深入，裂缝性油气藏将在世界范围内的油气勘探与开发中占据越来越重要的地位。裂缝性储层介质的等效介质模型建立了地震数据与储层特征的联系。传统的等效介质模型未考虑岩石孔隙的存在，也无法对实际裂缝性储层中地震波的频散和衰减给出合理解释，利用基于传统的等效介质模型所建立的裂缝性储层地震波传播理论能够获得裂缝方位和裂缝密度信息，但无法获得裂缝的几何结构信息（如裂缝尺寸、裂缝开度）。然而，裂缝几何结构信息对裂缝性储层的渗透性具有重要影响。因此，为了实现裂缝性储层的精细描述，应采用更为完善的裂缝性储层介质等效介质模型来研究裂缝性储层的相关问题，进而力求通过所建立的地震波传播理论获得诸如裂缝几何结构信息等更多的储层信息，从而实现储层的定性评价以及储层参数的定量反演。
裂缝性储层介质实际上是一种包含裂缝、岩石孔隙以及流体的含流体裂缝性孔隙介质，裂缝性孔隙介质的等效介质模型能够更精确地描述实际裂缝性储层介质的性质，通过相应的等效介质模型研究裂缝性孔隙介质中的波动现象，认清地震波传播规律和响应特征可为裂缝性储层预测、流体识别以及储层参数反演提供新的技术手段。本文分析了几种常见的裂缝性孔隙介质等效介质模型的特点，Chapman(2003)模型能够对地震波频散和衰减给出合理的解释，在描述裂缝性储层介质性质方面较其他模型更具有优势，而且，目前已有的实际应用研究表明，利用基于该模型所发展起来的储层参数反演方法能够从野外观测数据中反演出裂缝尺寸、裂缝密度和裂缝方位，该模型具有良好的应用前景。因而，本文对基于Chapman模型的裂缝性孔隙介质的相关问题进行了研究，主要包括裂缝性孔隙介质地震波传播特性分析、波场数值模拟、频变AVO分析以及储层参数反演。
为了充分认识地震波在裂缝性储层中的频散、衰减、方位各向异性等特性，本文开展了基于Chapman模型的裂缝性孔隙介质中地震波传播特性的研究。首先，从HTI介质的频率域弹性波方程出发，通过假定各种波的平面谐波解，推导基于Chapman模型的三维HTI介质中的Christoffel方程。然后，通过求解频散关系方程确定复速度，并给出相速度和逆品质因子的表达式。在此基础上，研究入射角、渗透率、裂缝尺寸、裂缝密度和流体类型对频散和衰减的影响，并对频散和衰减的方位特性及其影响因素进行分析，同时，从横波分裂的频率相关性和Thomsen参数的频率相关性两个角度来分析裂缝性储层介质的频率相关各向异性特征。研究结果表明：拟P波和拟SV波在裂缝性孔隙介质中传播时，将诱导流体在裂缝与孔隙之间发生流动，从而发生速度频散和衰减，拟SH波由于不挤压裂缝，不会诱导流体流动，从而不发生速度频散和衰减；由于裂缝的定向排列，速度频散和衰减将表现出方位各向异性，进而地震各向异性和横波分裂将表现出频率相关性；入射角、裂缝密度、渗透率、裂缝尺寸和流体类型对于地震波频散和衰减具有重要影响；入射角、裂缝密度、裂缝尺寸和流体类型也将影响相速度和衰减的方位特性；同时，与频率相关的各向异性和横波分裂也将受到裂缝密度、裂缝尺寸和流体类型的影响。
为了从波场正演角度揭示裂缝性孔隙介质中由地震波诱导的流体中观流动引起的地震波频散、衰减以及与频率相关的各向异性，研究地震波在裂缝性孔隙介质中的传播机理及其波场响应特征，本文提出了基于Chapman模型的裂缝性孔隙介质波场数值模拟的时间域方法。根据Chapman模型有效刚度参数是与频率相关的复数值这一特点，采用Zener粘弹性模型匹配这些有效刚度参数，进而通过求解与裂缝性孔隙介质等效的粘弹各向异性介质的波动方程进行波场数值模拟。数值模拟算例表明，所提出的数值模拟方法正确地模拟了裂缝性孔隙介质中地震波场的传播特征，从波场快照、波形图、合成地震记录中可以明显观测到地震波的频散和衰减特性。
地震波反射问题是地震勘探研究中的最基本也是最重要的问题之一，是反射波地震勘探的核心问题。本文研究了裂缝性孔隙介质中波的反射问题。以Chapman模型描述裂缝性孔隙介质，根据粘弹各向异性介质与基于Chapman模型的HTI介质在频率域本构方程的一致性，将粘弹各向异性介质中的波反射、透射系数的计算方法与Chapman模型相结合，给出了基于Chapman模型的HTI介质的地震波精确反射系数的计算方法。同时，采用前人提出的衰减各向异性介质反射波近似反射系数的计算公式对基于Chapman模型的HTI介质的反射系数进行了计算，两种方法在近似公式适用的范围内具有很好的一致性，表明了两种方法的有效性。在此基础上，就重要的储层参数对PP波频变反射系数的影响进行了分析，并进一步分析了反射PP波的频变AVO和频变AVOZ特性，研究表明：PP波频变反射系数对于裂缝密度、裂缝尺寸和孔隙度的变化具有良好的敏感性；PP波的频变AVO特性与上、下层介质的波阻抗关系有关，通过三个模型算例的计算，并参照Ren等(2009)对斑块饱和含气储层的分类方式，按照PP波反射系数的频率相关性将储层分为低频暗点储层、相移储层和低频亮点储层；PP波反射系数及其频变性质随着方位角的变化而变化。
对于裂缝性储层的研究，裂缝尺寸、裂缝密度、孔隙度以及裂缝方位都是储层识别和评价的关键参数，储层参数反演对于储层的定性预测和定量描述至关重要。本文根据PP波频变反射系数对裂缝密度、裂缝尺寸、孔隙度变化具有良好敏感性这一特点，开展了利用遗传算法通过裂缝性孔隙介质的频变AVO、频变AVOZ进行储层参数反演的研究。在已知裂缝方位的情况下，通过频变AVO数据反演了裂缝密度、裂缝尺寸和孔隙度；而在裂缝方位未知的情况下，通过频变AVOZ数据对裂缝方位、裂缝密度、裂缝尺寸和孔隙度进行了同时反演；针对裂缝性薄储层，通过频变AVO数据实现薄层厚度、裂缝密度、裂缝尺寸和孔隙度的同时反演。反演方法具有一定的稳定性和抗噪性，为裂缝性储层参数反演提供了新的技术手段。
With the deepening of the degree of oil and gas exploration and development,fractured reservoirs will become more and more important in the worldwide oil andgas exploration and development. The relationship between seismic data and reservoircharacteristics is build with the equivalent medium model for fractured reservoirsmedia. Traditional equivalent medium model do not consider the existence of thepores of the rock, cannot give a reasonable explanation for seismic wave dispersionand attenuation in actual fractured reservoirs, the information about fracture azimuthand fracture density can be got with the use of seismic wave propagation theory infractured reservoirs which is established based on traditional equivalent mediummodel, but the information about fracture geometry (such as fracture size, fractureaperture) can't be got. However, the information about fracture geometry has animportant effect on reservoir permeability. Therefore，in order to achieve sophisticateddescription of fractured reservoirs, the more perfect equivalent medium model offractured reservoirs media should be used to research issues related to fracturedreservoirs, and then by using the established seismic wave propagation theory toobtain more reservoir information such as the information about fracture geometry, soas to realize the qualitative evaluation of the reservoir and reservoir parametersinversion.
In fact, fractured reservoir is a kind of fractured porous media containing fractures,pore space of rocks and fluid. The equivalent medium models for fractured porousmedia can describe the property of the actual fractured reservoirs media moreaccurately. To research wave phenomena and recognize seismic wave propagation law and the response characteristics in fractured porous media through the correspondingequivalent medium model can provide new technical ideas for fractured reservoirprediction, fluid identification and reservoir parameters inversion. This paper analyzesthe characteristics of several common equivalent medium models for fractured porousmedia, the equivalent medium model of Chapman (2003) can give a reasonableexplanation for seismic wave dispersion and attenuation in fractured reservoirs, It hasmore advantages than other models on the description of the property of fracturedreservoirs media, moreover, the existing investigation for some practical applicationsshow that fracture size, fracture density and fracture azimuth can be inverted fromfield observation data with the use of the reservoir parameter inversion methodadvanced based on this model, this model has good application prospect. Therefore,this paper researches the issues related to fractured porous media based on Chapmanmodel, such as the characteristics analysis of seismic wave propagation in fracturedporous media, numerical simulation of wave field, frequency-dependent amplitudeversus offset analysis and reservoir parameter inversion.
In order to fully understand seismic wave dispersion, attenuation, azimuthanisotropy and other characteristics in fractured reservoirs, this paper carry outseismic wave propagation characteristics research in fractured porous media based onChapman model. First, begin with frequency domain wave equation for HTI media,the Christoffel equation in the three-dimension HTI media based on Chapman modelis derived by assuming plane harmonic solution of various wave. Then, throughsolving the dispersion relation equation to determine the complex velocity and givethe expressions of phase velocity and quality factor. On this basis, researching theeffect of incidence angle, permeability, fracture size, fracture density and fluid type onseismic wave dispersion and attenuation, and analyzing the azimuth characteristics ofdispersion and attenuation and its influencing factors, meanwhile, through frequency-dependent shear wave splitting and the frequency-dependent property of Thomsenparameters to analyze frequency-dependent anisotropy characteristics of fracturedreservoirs media. The results show that when quasi P wave or quasi SV wave througha fractured porous rock, It will induce fluid flow occurs between fractures and pores result in velocity dispersion and attenuation, as to quasi SH wave, It will not squeezethe fractures and induce fluid flow, thus no velocity dispersion and attenuation.Because of the fractures are aligned in space, then velocity dispersion and attenuationwill show azimuth anisotropy, and then seismic anisotropy and shear wave splittingwill show frequency-dependent. Incidence angle, fracture density, permeability,fracture size and fluid type have important effect on seismic wave dispersion andattenuation. Incidence angle, fracture density, fracture size and fluid type have alsoeffect on the azimuth characteristics of phase velocity and attenuation. Meanwhile,frequency-dependent anisotropy and shear wave splitting will also be affected byfracture density, fracture size and fluid type.
In order to reveal seismic wave dispersion, attenuation and frequency-dependentanisotropy caused by wave-induced mesoscopic fluid flow in fractured porous mediafrom the angle of the wave field forward modeling, and then to study the mechanismof seismic wave propagation and seismic response characteristics in fractured porousmedia, this paper advances a time-domain numerical simulation method of wave fieldin fractured porous media based on Chapman model. According to the characteristicsthat the effective stiffness coefficients of Chapman model are complex andfrequency-dependent, by using Zener models to fit the stiffness coefficients, and thencarrying out numerical simulation of wave field by solving the wave equation inviscoelastic anisotropic media which is equivalent to fractured porous media. Thenumerical simulation examples show that the advanced numerical simulation methodsimulates correctly seismic wave propagation characteristics in fractured porousmedia, seismic wave dispersion and attenuation characteristics can be observedclearly from the wave field snapshots, waveform and synthetic seismogram.
Seismic wave reflection problem is one of the most basic and most important issues,It is the core issue of the reflection survey. This paper studies reflection problem ofwaves in fractured porous media. By using the Chapman model to describe thefractured porous media, according to the consistency between constitutive equationfor viscoelastic anisotropic media in frequency domain and the HTI media based onChapman model, combining the calculation method of reflection and transmission coefficients of waves in viscoelastic anisotropic media and the Chapman model topresent the calculation method of exact reflection coefficients of seismic waves inHTI media based on Chapman model. Meanwhile, by using the approximate formulasof the reflection coefficients in attenuative anisotropic media to calculate thereflection coefficients in the HTI media based on the Chapman model, two methodshave good consistency within the scope of the approximate formula applicable, whichshow that two methods are both effective. On this basis, analyzing the effects ofimportant reservoir parameters on PP-wave frequency-dependent reflectioncoefficient and further analyzing the frequency-dependent AVO characteristics and thefrequency-dependent AVOZ characteristics of reflected PP-wave, the researches showthat the PP-wave frequency-dependent reflection coefficient is sensitive to variationsof fracture density, fracture size and porosity; frequency-dependent AVOcharacteristics of reflected PP-wave have some connection with the relation of waveimpedance between the overlying formation and underlying formation. Through thecalculations of three model examples, and referring to the classification method forpatchy-saturated gas reservoirs advanced by Ren et al (2009), according tofrequency-dependent property of PP-wave reflection coefficient, classifying fracturedreservoirs into three classification, they are low-frequency dim-out reservoirs,phase-shift reservoirs and low-frequency bright-spot reservoirs. PP-wave reflectioncoefficient and its frequency-dependent characteristics vary with changes in azimuth.
As to the studies on fractured reservoirs, fracture size, fracture density, porosityand fracture azimuth are all the key parameters in reservoir identification andevaluation, and the reservoir parameters inversion plays the significant role onqualitative prediction and quantitative description of fractured reservoirs. On the basisof the characteristics that PP-wave frequency-dependent reflection coefficient hasgood sensitivity to the variations of fracture density, fracture size and porosity, thispaper adopts genetic algorithm to carry out reservoir parameter inversion by using thefrequency-dependent AVO or frequency-dependent AVOZ in fractured porous media.In the case of known fracture azimuth, fracture density, fracture size and porosity areinverted by frequency-dependent AVO data. While in the case of unknown fracture azimuth, the fracture azimuth, fracture density, fracture size and porosity are invertedsimultaneously. For the thin fractured reservoir, the thickness of thin lay, fracturedensity, fracture size and porosity are inverted simultaneously. The advanced reservoirparameters inversion method in this paper has a certain stability and noise immunity,provides new techniques for reservoir parameters inversion in fractured reservoirs.
裂缝性储层介质实际上是一种包含裂缝、岩石孔隙以及流体的含流体裂缝性孔隙介质，裂缝性孔隙介质的等效介质模型能够更精确地描述实际裂缝性储层介质的性质，通过相应的等效介质模型研究裂缝性孔隙介质中的波动现象，认清地震波传播规律和响应特征可为裂缝性储层预测、流体识别以及储层参数反演提供新的技术手段。本文分析了几种常见的裂缝性孔隙介质等效介质模型的特点，Chapman(2003)模型能够对地震波频散和衰减给出合理的解释，在描述裂缝性储层介质性质方面较其他模型更具有优势，而且，目前已有的实际应用研究表明，利用基于该模型所发展起来的储层参数反演方法能够从野外观测数据中反演出裂缝尺寸、裂缝密度和裂缝方位，该模型具有良好的应用前景。因而，本文对基于Chapman模型的裂缝性孔隙介质的相关问题进行了研究，主要包括裂缝性孔隙介质地震波传播特性分析、波场数值模拟、频变AVO分析以及储层参数反演。
为了充分认识地震波在裂缝性储层中的频散、衰减、方位各向异性等特性，本文开展了基于Chapman模型的裂缝性孔隙介质中地震波传播特性的研究。首先，从HTI介质的频率域弹性波方程出发，通过假定各种波的平面谐波解，推导基于Chapman模型的三维HTI介质中的Christoffel方程。然后，通过求解频散关系方程确定复速度，并给出相速度和逆品质因子的表达式。在此基础上，研究入射角、渗透率、裂缝尺寸、裂缝密度和流体类型对频散和衰减的影响，并对频散和衰减的方位特性及其影响因素进行分析，同时，从横波分裂的频率相关性和Thomsen参数的频率相关性两个角度来分析裂缝性储层介质的频率相关各向异性特征。研究结果表明：拟P波和拟SV波在裂缝性孔隙介质中传播时，将诱导流体在裂缝与孔隙之间发生流动，从而发生速度频散和衰减，拟SH波由于不挤压裂缝，不会诱导流体流动，从而不发生速度频散和衰减；由于裂缝的定向排列，速度频散和衰减将表现出方位各向异性，进而地震各向异性和横波分裂将表现出频率相关性；入射角、裂缝密度、渗透率、裂缝尺寸和流体类型对于地震波频散和衰减具有重要影响；入射角、裂缝密度、裂缝尺寸和流体类型也将影响相速度和衰减的方位特性；同时，与频率相关的各向异性和横波分裂也将受到裂缝密度、裂缝尺寸和流体类型的影响。
为了从波场正演角度揭示裂缝性孔隙介质中由地震波诱导的流体中观流动引起的地震波频散、衰减以及与频率相关的各向异性，研究地震波在裂缝性孔隙介质中的传播机理及其波场响应特征，本文提出了基于Chapman模型的裂缝性孔隙介质波场数值模拟的时间域方法。根据Chapman模型有效刚度参数是与频率相关的复数值这一特点，采用Zener粘弹性模型匹配这些有效刚度参数，进而通过求解与裂缝性孔隙介质等效的粘弹各向异性介质的波动方程进行波场数值模拟。数值模拟算例表明，所提出的数值模拟方法正确地模拟了裂缝性孔隙介质中地震波场的传播特征，从波场快照、波形图、合成地震记录中可以明显观测到地震波的频散和衰减特性。
地震波反射问题是地震勘探研究中的最基本也是最重要的问题之一，是反射波地震勘探的核心问题。本文研究了裂缝性孔隙介质中波的反射问题。以Chapman模型描述裂缝性孔隙介质，根据粘弹各向异性介质与基于Chapman模型的HTI介质在频率域本构方程的一致性，将粘弹各向异性介质中的波反射、透射系数的计算方法与Chapman模型相结合，给出了基于Chapman模型的HTI介质的地震波精确反射系数的计算方法。同时，采用前人提出的衰减各向异性介质反射波近似反射系数的计算公式对基于Chapman模型的HTI介质的反射系数进行了计算，两种方法在近似公式适用的范围内具有很好的一致性，表明了两种方法的有效性。在此基础上，就重要的储层参数对PP波频变反射系数的影响进行了分析，并进一步分析了反射PP波的频变AVO和频变AVOZ特性，研究表明：PP波频变反射系数对于裂缝密度、裂缝尺寸和孔隙度的变化具有良好的敏感性；PP波的频变AVO特性与上、下层介质的波阻抗关系有关，通过三个模型算例的计算，并参照Ren等(2009)对斑块饱和含气储层的分类方式，按照PP波反射系数的频率相关性将储层分为低频暗点储层、相移储层和低频亮点储层；PP波反射系数及其频变性质随着方位角的变化而变化。
对于裂缝性储层的研究，裂缝尺寸、裂缝密度、孔隙度以及裂缝方位都是储层识别和评价的关键参数，储层参数反演对于储层的定性预测和定量描述至关重要。本文根据PP波频变反射系数对裂缝密度、裂缝尺寸、孔隙度变化具有良好敏感性这一特点，开展了利用遗传算法通过裂缝性孔隙介质的频变AVO、频变AVOZ进行储层参数反演的研究。在已知裂缝方位的情况下，通过频变AVO数据反演了裂缝密度、裂缝尺寸和孔隙度；而在裂缝方位未知的情况下，通过频变AVOZ数据对裂缝方位、裂缝密度、裂缝尺寸和孔隙度进行了同时反演；针对裂缝性薄储层，通过频变AVO数据实现薄层厚度、裂缝密度、裂缝尺寸和孔隙度的同时反演。反演方法具有一定的稳定性和抗噪性，为裂缝性储层参数反演提供了新的技术手段。
With the deepening of the degree of oil and gas exploration and development,fractured reservoirs will become more and more important in the worldwide oil andgas exploration and development. The relationship between seismic data and reservoircharacteristics is build with the equivalent medium model for fractured reservoirsmedia. Traditional equivalent medium model do not consider the existence of thepores of the rock, cannot give a reasonable explanation for seismic wave dispersionand attenuation in actual fractured reservoirs, the information about fracture azimuthand fracture density can be got with the use of seismic wave propagation theory infractured reservoirs which is established based on traditional equivalent mediummodel, but the information about fracture geometry (such as fracture size, fractureaperture) can't be got. However, the information about fracture geometry has animportant effect on reservoir permeability. Therefore，in order to achieve sophisticateddescription of fractured reservoirs, the more perfect equivalent medium model offractured reservoirs media should be used to research issues related to fracturedreservoirs, and then by using the established seismic wave propagation theory toobtain more reservoir information such as the information about fracture geometry, soas to realize the qualitative evaluation of the reservoir and reservoir parametersinversion.
In fact, fractured reservoir is a kind of fractured porous media containing fractures,pore space of rocks and fluid. The equivalent medium models for fractured porousmedia can describe the property of the actual fractured reservoirs media moreaccurately. To research wave phenomena and recognize seismic wave propagation law and the response characteristics in fractured porous media through the correspondingequivalent medium model can provide new technical ideas for fractured reservoirprediction, fluid identification and reservoir parameters inversion. This paper analyzesthe characteristics of several common equivalent medium models for fractured porousmedia, the equivalent medium model of Chapman (2003) can give a reasonableexplanation for seismic wave dispersion and attenuation in fractured reservoirs, It hasmore advantages than other models on the description of the property of fracturedreservoirs media, moreover, the existing investigation for some practical applicationsshow that fracture size, fracture density and fracture azimuth can be inverted fromfield observation data with the use of the reservoir parameter inversion methodadvanced based on this model, this model has good application prospect. Therefore,this paper researches the issues related to fractured porous media based on Chapmanmodel, such as the characteristics analysis of seismic wave propagation in fracturedporous media, numerical simulation of wave field, frequency-dependent amplitudeversus offset analysis and reservoir parameter inversion.
In order to fully understand seismic wave dispersion, attenuation, azimuthanisotropy and other characteristics in fractured reservoirs, this paper carry outseismic wave propagation characteristics research in fractured porous media based onChapman model. First, begin with frequency domain wave equation for HTI media,the Christoffel equation in the three-dimension HTI media based on Chapman modelis derived by assuming plane harmonic solution of various wave. Then, throughsolving the dispersion relation equation to determine the complex velocity and givethe expressions of phase velocity and quality factor. On this basis, researching theeffect of incidence angle, permeability, fracture size, fracture density and fluid type onseismic wave dispersion and attenuation, and analyzing the azimuth characteristics ofdispersion and attenuation and its influencing factors, meanwhile, through frequency-dependent shear wave splitting and the frequency-dependent property of Thomsenparameters to analyze frequency-dependent anisotropy characteristics of fracturedreservoirs media. The results show that when quasi P wave or quasi SV wave througha fractured porous rock, It will induce fluid flow occurs between fractures and pores result in velocity dispersion and attenuation, as to quasi SH wave, It will not squeezethe fractures and induce fluid flow, thus no velocity dispersion and attenuation.Because of the fractures are aligned in space, then velocity dispersion and attenuationwill show azimuth anisotropy, and then seismic anisotropy and shear wave splittingwill show frequency-dependent. Incidence angle, fracture density, permeability,fracture size and fluid type have important effect on seismic wave dispersion andattenuation. Incidence angle, fracture density, fracture size and fluid type have alsoeffect on the azimuth characteristics of phase velocity and attenuation. Meanwhile,frequency-dependent anisotropy and shear wave splitting will also be affected byfracture density, fracture size and fluid type.
In order to reveal seismic wave dispersion, attenuation and frequency-dependentanisotropy caused by wave-induced mesoscopic fluid flow in fractured porous mediafrom the angle of the wave field forward modeling, and then to study the mechanismof seismic wave propagation and seismic response characteristics in fractured porousmedia, this paper advances a time-domain numerical simulation method of wave fieldin fractured porous media based on Chapman model. According to the characteristicsthat the effective stiffness coefficients of Chapman model are complex andfrequency-dependent, by using Zener models to fit the stiffness coefficients, and thencarrying out numerical simulation of wave field by solving the wave equation inviscoelastic anisotropic media which is equivalent to fractured porous media. Thenumerical simulation examples show that the advanced numerical simulation methodsimulates correctly seismic wave propagation characteristics in fractured porousmedia, seismic wave dispersion and attenuation characteristics can be observedclearly from the wave field snapshots, waveform and synthetic seismogram.
Seismic wave reflection problem is one of the most basic and most important issues,It is the core issue of the reflection survey. This paper studies reflection problem ofwaves in fractured porous media. By using the Chapman model to describe thefractured porous media, according to the consistency between constitutive equationfor viscoelastic anisotropic media in frequency domain and the HTI media based onChapman model, combining the calculation method of reflection and transmission coefficients of waves in viscoelastic anisotropic media and the Chapman model topresent the calculation method of exact reflection coefficients of seismic waves inHTI media based on Chapman model. Meanwhile, by using the approximate formulasof the reflection coefficients in attenuative anisotropic media to calculate thereflection coefficients in the HTI media based on the Chapman model, two methodshave good consistency within the scope of the approximate formula applicable, whichshow that two methods are both effective. On this basis, analyzing the effects ofimportant reservoir parameters on PP-wave frequency-dependent reflectioncoefficient and further analyzing the frequency-dependent AVO characteristics and thefrequency-dependent AVOZ characteristics of reflected PP-wave, the researches showthat the PP-wave frequency-dependent reflection coefficient is sensitive to variationsof fracture density, fracture size and porosity; frequency-dependent AVOcharacteristics of reflected PP-wave have some connection with the relation of waveimpedance between the overlying formation and underlying formation. Through thecalculations of three model examples, and referring to the classification method forpatchy-saturated gas reservoirs advanced by Ren et al (2009), according tofrequency-dependent property of PP-wave reflection coefficient, classifying fracturedreservoirs into three classification, they are low-frequency dim-out reservoirs,phase-shift reservoirs and low-frequency bright-spot reservoirs. PP-wave reflectioncoefficient and its frequency-dependent characteristics vary with changes in azimuth.
As to the studies on fractured reservoirs, fracture size, fracture density, porosityand fracture azimuth are all the key parameters in reservoir identification andevaluation, and the reservoir parameters inversion plays the significant role onqualitative prediction and quantitative description of fractured reservoirs. On the basisof the characteristics that PP-wave frequency-dependent reflection coefficient hasgood sensitivity to the variations of fracture density, fracture size and porosity, thispaper adopts genetic algorithm to carry out reservoir parameter inversion by using thefrequency-dependent AVO or frequency-dependent AVOZ in fractured porous media.In the case of known fracture azimuth, fracture density, fracture size and porosity areinverted by frequency-dependent AVO data. While in the case of unknown fracture azimuth, the fracture azimuth, fracture density, fracture size and porosity are invertedsimultaneously. For the thin fractured reservoir, the thickness of thin lay, fracturedensity, fracture size and porosity are inverted simultaneously. The advanced reservoirparameters inversion method in this paper has a certain stability and noise immunity,provides new techniques for reservoir parameters inversion in fractured reservoirs.
引文
Agersborg R, Jakobsen M, Ruud B O, et al. Effects of pore fluid pressure on the seismic responseof a fractured carbonate reservoir[J]. Studia Geophysica et Geodaetica,2007,51(1):89–118.
Al-Anboori A. Anisotropy, focal mechanisms and state of stress in an oilfield: passive seismicmonitoring in Oman[D]. University of Leeds.
Al-Harrasi O H, Kendall J M, Chapman M. Fracture characterization using frequency-dependentshear wave anisotropy analysis of microseismic data[J]. Geophysical Journal International,2011,185:1059–1070.
Ali A, Jakobsen M. Seismic characterization of reservoirs with multiple fracture sets usingvelocity and attenuation anisotropy data[J]. Journal of Applied Geophysics,2011,75:590-602.
Ali A, Jakobsen M. Anisotropic permeability in fractured reservoirs from frequency-dependentseismic Amplitude Versus Angle and Azimuth data[J]. Geophysical Prospecting,2014,62:293–314.
Auld B A. Acoustic Fields and Waves in Solids[M]. Malabar, FL: Krieger Publishing,1990.
Baird A F, Kendall J-M, Angus D A. Frequency-dependent seismic anisotropy due to fractures:Fluid flow versus scattering[J]. Geophysics,2013,78(2): WA111-WA122.
Ba J, Nie J X, Cao H, et al. Mesoscopic fluid flow simulation in double-porosity rocks[J].Geophys. Res. Lett.,2008,35, L04303.
Batzle M L, Han D H, Hofmann R. Fluid mobility and frequency-dependent seismic velocity—Direct measurements[J]. Geophysics,2006,71(1): N1–N9.
Behura J, Tsvankin I. Reflection coefficients in attenuative anisotropic media[J]. Geophysics,2009,74(5): WB193—WB202.
Beretta M M, Bernasconi G, Drufuca G. AVO and AVA inversion for fractured reservoircharacterization[J]. Geophysics,2002,67(1):300–306.
Berryman J G. Seismic waves in rocks with fluids and fractures[J]. Geophysical JournalInternational,2007,171:954–974.
Borcherdt R D. Viscoelastic waves in layered media[M]. Cambridge University Press,2009.
Brajanovski M, Gurevich B, Schoenberg M. A model for P-wave attenuation and dispersion in aporous medium permeated by aligned fractures[J]. Geophysical Journal International,2005,163:372–384.
Brajanovski M, Müller T M, Parra J O. A model for strong attenuation and dispersion of seismicP-waves in a partially saturated fractured reservoir[J]. Science China-Physics, Mechanics andAstronomy,2010,53:1383-1387.
Brown R J S, Korringa J. On the dependence of the elastic properties of a porous rock on thecompressibility of the pore fluid[J]. Geophysics,1975,40:608–616.
Carcione J M. Wave propogation in anisotropic linear viscoelastic media: theory and simulatedwavefields[J]. Geophys. J. Int.,1990,101:739-750.
Carcione J M. Constitutive model and wave equations for linear, viscoelastic, anisotropic media[J].Geophysics,1995,60(2):537-548.
Carcione J M. Wave propagation in anisotropic, saturated porous media: Plane wave theory andnumerical simulation[J]. Journal of the Acoustical Society of America,1996,99:2655–2666.
Carcione J M. Reflection and transmission of qP-qS plane waves at a plane boundary betweenviscoelastic transversely isotropic media[J]. Geophysics,1997,129:669–680.
Carcione J M. Staggered mesh for the anisotropic and viscoelastic wave equation[J]. Geophysics,1999,64(6):1863-1866.
Carcione J M. A model for seismic velocity and attenuation in petroleum source rocks[J].Geophysics,2000,65:1080-1092.
Carcione J M. Wave fields in real media: wave propagation in anisotropic, anelastic and porousmedia[M]. Pergamon,2001a.
Carcione J M. AVO effects of a hydrocarbon source-rock layer[J]. Geophysics,2001b,66(2):283-293.
Carcione J M. Wave fields in real media: Theory and numerical simulation of wave propagation inanisotropic, anelastic, porous and electromagnetic media,2nd ed.: Elsevier Scientific Publ. Co.,Inc.,2007.
Carcione, J M., and Gurevich, B.,2011, Differential form and numerical implementation of Biot'sporoelasticity equations with squirt dissipation, Geophysics,76, N55-N64.
Carcione J M, Helle H B. Numerical solution of the poroviscoelastic wave equation on a staggeredmesh[J]. J. Comput. Phys.,1999,154:520–527.
Carcione J M, Herman G C, ten Kroode A P E. Seismic modeling[J]. Geophysics,2002,67(4):1304-1325.
Carcione J M, Kosloff D, Kosloff R. Wave propagation in a linear viscoelastic medium[J]. Geophy.J. Int.,1988,95:597-611.
Carcione J M, Picotti S, Gei D, Rossi G. Physics and seismic modeling for monitoring CO2storage[J]. Pure and Applied Geophysics,2006,163:175–207.
Cerjan C, Kosloff D, Kosloff R, et al. A nonreflecting boundary condition for discrete acoustic andelastic wave equation[J]. Geophysics,1985,50(4):705-708.
Chapman M. Frequency dependent anisotropy due to mesoscale fractures in the presence ofequant porosity[J]. Geophysical Prospecting,2003,51:369–379.
Chapman M. Modeling the effect of multiple fracture sets of mesoscale fractures in porous rock onfrequency-dependent anisotropy[J]. Geophysics,2009,74(6): D97–D103.
Chapman M, Liu E, Li X Y. The influence of fluid-sensitive dispersion and attenuation on AVOanalysis[J]. Geophysics,2006,167,89-105.
Chapman M, Maultzsch S, Liu E, et al. The effect of fluid saturation in an anisotropic multi-scaleequant porosity model[J]. Journal of Applied Geophysics,2003,54:191-202.
Chapman M, Zatsepin S V, Crampin S. Derivation of a microstructural poroelastic model[J].Geophysical Journal International,2002,151:427–451.
Du Q Z, Wang X M, Ba J, et al. An equivalent medium model for wave simulation in fracturedporous rocks[J]. Geophysical Prospecting,2012,60(5):940-956.
Duxbury A, White D, Samson C, et al. Fracture mapping using seismic amplitude variation withoffset and azimuth analysis at the Weyburn CO2storage site[J]. Geophysics,2012,77(6):N17-N28.
Dvorkin J, Mavko G. Modeling attenuation in reservoir and non-reservoir rock[J]. The LeadingEdge,2006,25:194–197.
Eshelby J D. The determination of the elastic field of an ellipsoidal inclusion, and relatedproblems[J]. Proceedings of the Royal Society of London,1957(A241):376-396.
Galvin R J, Gurevich B. Scattering of a longitudinal wave by a circular crack in a fluid-saturatedporous medium[J]. International Journal of Solids and Structures,2007,44:7389–7398.
Galvin R J, Gurevich B. Effective properties of a poroelastic medium containing a distribution ofaligned cracks[J]. Journal of Geophysical Research—Solid Earth,2009,114, B07305.
Giroux B. Performance of convolutional perfectly matched layers for pseudospectral time domainporoviscoelastic schemes[J]. Computers&Geosciences,2012,45:149-160.
Gurevich B. Elastic properties of saturated porous rocks with aligned fractures[J]. Journal ofApplied Geophysics,2003,54:203-218.
Gurevich B, Brajanovski R, Galvin R J, et al. P-wave dispersion and attenuation in fractured andporous reservoirs—Poroelasticity approach[J]. Geophysical Prospecting,2009,57(2):225-237.
Hall S A, Kendall J M. Fracture characterization at Valhall: Application of P-wave amplitudevariation with offset and azimuth (AVOA) analysis to a3D ocean-bottom data set[J].Geophysics,2003,68(4):1150–1160.
Hudson J A. Wave speeds and attenuation of elastic waves in material containing cracks[J].Geophys. J. R. Astron. Soc.,1981,64:133–50.
Hudson J A, Liu E, Crampin S. The mechanical properties of materials with interconnected cracksand pores[J]. Geophysical Journal International,1996,124:105–112.
Hudson J A, Pointer T, Liu E. Effective-medium theories for fluid-saturated materials with alignedcracks[J]. Geophysical Prospecting,2001,49:509–522.
Hu Y, McMechan G. Sensitivity of three-component3D finite-difference elastic seismic modelingto inclusion parameters in HTI and TTI media with high inclusion density[J]. Geophysics,2010,75(2): T49-T61.
Jakobsen M. The interacting inclusion model of wave-induced fluid flow[J]. Geophysical JournalInternational,2004,158:1168–1176.
Jakobsen M, Chapman M. Unified theory of global flow and squirt flow in cracked porousmedia[J]. Geophysics,2009,74(2): WA65–WA76.
Jakobsen M, Johansen T A, McCann C. The acoustic signature of fluid flow in complex porousmedia[J]. Journal of Applied Geophysics,2003,54:219–246.
Jilek P. Converted PS-wave reflection coefficients in weakly anisotropic media[J]. Pure andApplied Geophysics,2002,159:1527-1562.
Johansen T A, Jakobsen M, Agersborg R, et al. The P-P and P-S response of a fractured reservoir.74th Annual International Meeting, SEG, Expanded Abstracts,2004,199–202.
Johnson D L. Theory of frequency dependent acoustics in patchy saturated porous media[J].Journal of the Acoustical Society of America,2001,110:682–694.
Klin P, Priolo E, Seriani G. Numerical simulation of seismic wave propagation in realistic3-Dgeo-models with a Fourier pseudo-spectral method[J]. Geophysical Journal International,2010,183:905-922.
Kosloff D, Baysal E. Forward modeling by the Fourier method[J]. Geophysics,1982,47:1402–1412.
Lambert M-A, Saenger E H, Quintal B, et al. Numerical simulation of ambient seismic wavefieldmodification caused by pore-fluid effects in an oil reservoir[J]. Geophysics,2013,78(1):T41–T52.
Liu C, Liu Y W, Feng X, et al. Constructing the convex quadratic function for the evaluation ofcrack density of HTI media using P-and converted waves[J]. J. Geophys. Eng.,2012,9:729–736.
Liu E, Queen J H, Li X Y, et al. Observation and analysis of frequency-dependent anisotropy froma multicomponent VSP at Bluebell-Altamont field[J]. Journal of Applied Geophysics,2003,54:319–333.
Liu J, Ba J, Ma J W, et al. An analysis of seismic attenuation in random porous media[J].Science China—Physics, Mechanics&Astronomy,2010,53:628–637.
Liu L F, Cao S Y, Wang L. Poroelastic analysis of frequency-dependent amplitude-versus-offsetvariations[J]. Geophysics,2011,76(3): C31-C40.
Liu X, Greenhalgh S, Zhou B. Transient solution for poro-viscoacoustic wave propagation indouble porosity media and its limitations[J]. Geophysical Journal International,2009,178:375–393.
Liu X, Greenhalgh S, Zhou B. Approximating the wave moduli of double porosity media at lowfrequencies by a single Zener or Kelvin-Voigt element[J]. Geophysical Journal International,2010,181:391–398.
Maultzsch S, Chapman M, Liu E, et al. Modelling frequency-dependent seismic anisotropy influid-saturated rock with aligned fractures: Implication of fracture size estimation fromanisotropic measurements[J]. Geophysical Prospecting,2003,51:381–392.
Maultzsch S, Chapman M, Liu E, et al. Modelling and analysis of attenuation anisotropy inmulti-azimuth VSP data from the Clair field[J]. Geophysical Prospecting,2007,55:627–642.
Mavko G, Jizba D. Estimating grain-scale fluid effects on velocity dispersion in rocks[J].Geophysics,1991,56:1940–1949.
Mavko G, Mukerji T, Dvorkin J. The rock physics handbook[M].2rd ed. Cambridge UniversityPress,2009.
Muir F, Dellinger J, Etgen J. et al. Modeling elastic fields across irregular boundaries[J].Geophysics,1992,57:1189-1193.
Müller T M, Gurevich B. One-dimensional random patchy saturation model for velocity andattenuation in porous rocks[J]. Geophysics,2004,69:1166–1172.
Müller T M, Gurevich B, Lebedev M. Seismic wave attenuation and dispersion resulting fromwave-induced flow in porous rocks—A review[J]. Geophysics,2010,75(5):75A147-75A164.
Mürphy W F I, Winkler K W, Kleinberg R L. Acoustic relaxation in sedimentary rocks:Dependence on grain contacts and fluid saturation[J]. Geophysics,1986,51:757–766.
Nelson R A. Geological analysis of naturally fractured reservoirs[M].2rd ed. Gulf Publ. Co.,2001.
Parra J O. Poroelastic model to relate seismic wave attenuation and dispersion to permeabilityanisotropy [J]. Geophysics,2000,65(1):202-210.
Picotti S, Carcione J M, Gei D, et al. Seismic modeling to monitor CO2geological storage: TheAtzbach-Schwanenstadt gas field, J. Geophys. Res.,2012,117, B06103.
Picotti S, Carcione J M, Rubino J G, et al. P-wave seismic attenuation by slow-wave diffusion:Numerical experiments in partially saturated rocks[J]. Geophysics,2007,72(4): N11-N21.
Picotti S, Carcione J M, Rubino J G, et al. A visoelastic representation of wave attenuation inporous media[J]. Computers&Geosciences,2010,36:44–53.
Pride S R, Berryman J G, Harris J M. Seismic attenuation due to wave induced flow[J]. Journal ofGeophysical Research,2004,109(B1), B01201.
Pyrak-Nolte L J, Nolte D D. Frequency dependence of fracture stiffness[J]. Geophysical ResearchLetters,1992,19:325–328.
Quintal B, Schmalholz S M, Podladchikov Y Y. Low-frequency reflections from a thin layer withhigh attenuation caused by interlayer flow[J]. Geophysics,2009,74(1): N14–N22.
Quintal B, Steeb H, Frehner M, et al. Quasi-static finite element modeling of seismic attenuationand dispersion due to wave-induced fluid flow in poroelastic media[J]. Journal of GeophysicalResearch Solid Earth,2011,116, B01201.
Ren H T, Goloshubin G, Hilterman F J. Poroelastic analysis of amplitude-versus-frequencyvariations[J]. Geophysics,2009,74(6): N41-N48.
Reshef M, Kosloff D, Edwards M, et al. Three-dimensional acoustic modeling by the Fouriermethod[J]. Geophysics,1988,53:1175-1183.
Robertsson J O A, Blance J O, Symes W W. Viscoelastic finite-difference modeling[J].Geophysics,1994,59:1444-1456.
Rubino J G, Ravazzoli C L, Santos J E. Equivalent viscoelastic solids for heterogeneousfluid-saturated porous rocks[J]. Geophysics,2009,74(1): N1–N13.
Rubino J G, Velis D R, Sacchi M D. Numerical analysis of wave-induced fluid flow effects onseismic data: Application to monitoring of CO2storage at the Sleipner field[J]. J. Geophys. Res.,2011,116, B03306.
Rüger A. P-wave reflection coefficient for transversely isotropic models with vertical andhorizontal axis of symmetry[J]. Geophysics,1997,62:713-722.
Rüger A. Variation of P-wave reflectivity with offset and azimuth in anisotropic media[J].Geophysics,1998,63:935–947.
Ruud B O, Hestholm S. Modeling seismic waves in orthorombic, viscoelastic media by finitedifferences.75th Annual International Meeting, SEG, Expanded Abstracts,2005,1771–1774.
Santos J E, Rubino J G, Ravazzoli C L. A numerical upscaling procedure to estimate effectiveplane wave and shear moduli in heterogeneous fluid-saturated poroelastic media[J]. ComputerMethods in Applied Mechanics and Engineering,2009,198:2067–2077.
Santos J E, Rubino J G, Ravazzoli C L. Modeling mesoscopic attenuation in a highlyheterogeneous Biot’s medium employing an equivalent viscoelastic model.78th AnnualInternational Meeting, SEG, Expanded Abstracts,2008,2112–2116.
Schoenberg M. Elastic wave behavior across linear slip interfaces[J]. J. Acoust. Soc. Am.,1980,68:1516–1521.
Schoenberg M, Sayers C M. Seismic anisotropy of fractured rock[J]. Geophysics,1995,60(1):204-211.
Sil S, Sen M K, Gurevich B. Analysis of fluid substitution in a porous and fractured medium[J].Geophysics,2011,76(3): WA157-WA166.
Thomsen L. Weak elastic anisotropic[J]. Geophysics,1986,51(10):1954-1966.
Thomsen L. Elastic anisotropy due to aligned cracks in porous rock[J]. Geophysical Prospecting,1995,43:805-829.
Tonn R. The determination of the seismic quality Q from VSP data: A comparison of differentcomputational methods[J]. Geophysical Prospecting,1991,39:1-27.
Tsvankin I. Reflection moveout and parameter estimation for horizontal transverse isotropy[J].Geophysics,62(2):614-629.
Van Der Kolk C M, Guest W S, Potters, J H H M. The3D shear experiment over the Natih field inOman: The effect of fracture-filling fluids on shear propagation[J]. Geophysical Prospecting,2001,49:179-197.
Vavrycuk V, Psencik I. PP-wave reflection coefficients in weakly anisotropic elastic media[J].Geophysics,1998,63:2129-2141.
White J E. Computed seismic speeds and attenuation in rocks with partial gas saturation[J].Geophysics,1975,40:224–232.
White J E, Mikhaylova N G, Lyakhovitskiy F M. Low-frequency seismic waves in fluid-saturatedlayered rocks[J]. Physics of the Solid Earth,1976,11:654–659.
Winkler K. Dispersion analysis of velocity and attenuation in Berea sandstone[J]. Journal ofGeophysical Research,1985,90:6793–6800.
Witte D, Richards P G. Contribution to the pseudo-spectral method for computing syntheticseismograms. Expanded Abstracts of57th SEG Annual Int. Mtg,1987,517-519.
Zhang S Q, Han L G, Liu C C, et al. Inverting reservoir parameters in a two-phase fracturedmedium with a niche genetic algorithm[J]. Applied Geophysics,2012,9(4):440-450.
巴晶，曹宏，姚逢昌等.利用错格虚谱法模拟礁、滩相储层地震波场[J].石油地球物理勘探,2010,45(2):177-184.
巴晶，卢明辉，胡彬等.黏弹双相介质中的松弛骨架模型[J].地球物理学报,2008,51(5):1527-1537.
杜启振.各向异性黏弹性介质伪谱法波场模拟[J]，物理学报，2004,53(12):4428-4434.
杜启振，孔丽云，韩世春.裂缝诱导各向异性双孔隙介质波场传播特征[J].地球物理学报，2009,52(4):1049-1058.
杜启振，刘莲莲，孙晶波.各向异性粘弹性孔隙介质地震波场伪谱法正演模拟[J].物理学报，2007，56(10):6143-6149.
郭智奇.粘弹各向异性介质波场模拟与储层信息研究[D].长春：吉林大学，2008.
郭智奇，刘财，杨宝俊等.粘弹各向异性介质中地震波场模拟与特征[J].地球物理学进展，2007,22(3):804-810.
孔丽云，王一博，杨慧珠.裂缝诱导HTI双孔隙介质中的裂缝参数分析[J].地球物理学报，2012,55(1):189-196.
雷英杰，张善文，李续武等. Matlab遗传算法工具箱及应用[M].西安：西安电子科技大学出版社，2005.
李敏强，寇纪淞，林丹等.遗传算法的基本理论与应用[M].北京：科学出版社，2002.
刘财，郭智奇，杨宝俊等.粘弹各向异性介质中波的反射与透射问题分析[J].地球物理学报，2007,50(4):1216-1224.
刘财，兰慧田，郭智奇等．基于改进BISQ机制的双相HTI介质波传播伪谱法模拟与特征分析[J]．地球物理学报，2013，56（10）：3461-3473.
刘财，刘宇巍，冯晅等.基于方位相交的纵波AVA数据运用SVD反演HIT介质裂缝密度[J].吉林大学学报（地球科学版），2013,43(5)：1655-1662.
刘炯,马坚伟,杨慧珠等.缝洞型储层错格伪谱法模拟[J].石油地球物理勘探,2008,43(6):723-727.
刘洋，李承楚.双相各向异性介质中弹性波传播伪谱法数值模拟研究.地震学报，2000,22(2):132-138.
聂建新,杨顶辉,杨慧珠.基于非饱和多孔隙介质BISQ模型的储层参数反演[J].地球物理学报,2004,47(6):1101-1105.
石耀霖，金文.面波频散反演地球内部构造的遗传算法[J].地球物理学报，1995，38(2):189-198.
轩义华，何樵登，孟庆生等.基于BISQ机制的双相EDA介质的波场分析[J].石油地球物理勘探，2006,41(5):550-556.
张广智,陈怀震,印兴耀等.基于各向异性AVO的裂缝弹性参数叠前反演方法[J].吉林大学学报（地球科学版），2012,42(3):845-851.
张军舵.孔隙流体介质波动方程正演模拟与参数反演[D].北京：中国石油大学，2008.
曾联波，柯式镇，刘洋.低渗透油气[M].北京：石油工业出版社，2010.
周明，孙树栋.遗传算法原理及应用[M].北京：国防工业出版社，1999.
朱培民,王家映,於文辉等.用纵波AVO数据反演储层裂隙密度参数[J].石油物探，2001,40(2)：1-12.
张文生，何樵登.二维横向各向同性介质的伪谱法正演模拟[J].石油地球物理勘探,1998,33(3):310-319.
张显文，王德利，王者江等.基于BISQ机制三维双相正交裂隙各向异性介质衰减及频散方位特性研究[J].地球物理学报，2010,53(10):2452-2459.
张显文.储层地震预测技术与储层地球物理参数反演[D].长春：吉林大学，2010.
赵志新,徐纪人,堀内茂木.错格实数傅立叶变换微分算子及其在非均匀介质波动传播研究中的应用[J].地球物理学报,2003,46(2):234-240.
Al-Anboori A. Anisotropy, focal mechanisms and state of stress in an oilfield: passive seismicmonitoring in Oman[D]. University of Leeds.
Al-Harrasi O H, Kendall J M, Chapman M. Fracture characterization using frequency-dependentshear wave anisotropy analysis of microseismic data[J]. Geophysical Journal International,2011,185:1059–1070.
Ali A, Jakobsen M. Seismic characterization of reservoirs with multiple fracture sets usingvelocity and attenuation anisotropy data[J]. Journal of Applied Geophysics,2011,75:590-602.
Ali A, Jakobsen M. Anisotropic permeability in fractured reservoirs from frequency-dependentseismic Amplitude Versus Angle and Azimuth data[J]. Geophysical Prospecting,2014,62:293–314.
Auld B A. Acoustic Fields and Waves in Solids[M]. Malabar, FL: Krieger Publishing,1990.
Baird A F, Kendall J-M, Angus D A. Frequency-dependent seismic anisotropy due to fractures:Fluid flow versus scattering[J]. Geophysics,2013,78(2): WA111-WA122.
Ba J, Nie J X, Cao H, et al. Mesoscopic fluid flow simulation in double-porosity rocks[J].Geophys. Res. Lett.,2008,35, L04303.
Batzle M L, Han D H, Hofmann R. Fluid mobility and frequency-dependent seismic velocity—Direct measurements[J]. Geophysics,2006,71(1): N1–N9.
Behura J, Tsvankin I. Reflection coefficients in attenuative anisotropic media[J]. Geophysics,2009,74(5): WB193—WB202.
Beretta M M, Bernasconi G, Drufuca G. AVO and AVA inversion for fractured reservoircharacterization[J]. Geophysics,2002,67(1):300–306.
Berryman J G. Seismic waves in rocks with fluids and fractures[J]. Geophysical JournalInternational,2007,171:954–974.
Borcherdt R D. Viscoelastic waves in layered media[M]. Cambridge University Press,2009.
Brajanovski M, Gurevich B, Schoenberg M. A model for P-wave attenuation and dispersion in aporous medium permeated by aligned fractures[J]. Geophysical Journal International,2005,163:372–384.
Brajanovski M, Müller T M, Parra J O. A model for strong attenuation and dispersion of seismicP-waves in a partially saturated fractured reservoir[J]. Science China-Physics, Mechanics andAstronomy,2010,53:1383-1387.
Brown R J S, Korringa J. On the dependence of the elastic properties of a porous rock on thecompressibility of the pore fluid[J]. Geophysics,1975,40:608–616.
Carcione J M. Wave propogation in anisotropic linear viscoelastic media: theory and simulatedwavefields[J]. Geophys. J. Int.,1990,101:739-750.
Carcione J M. Constitutive model and wave equations for linear, viscoelastic, anisotropic media[J].Geophysics,1995,60(2):537-548.
Carcione J M. Wave propagation in anisotropic, saturated porous media: Plane wave theory andnumerical simulation[J]. Journal of the Acoustical Society of America,1996,99:2655–2666.
Carcione J M. Reflection and transmission of qP-qS plane waves at a plane boundary betweenviscoelastic transversely isotropic media[J]. Geophysics,1997,129:669–680.
Carcione J M. Staggered mesh for the anisotropic and viscoelastic wave equation[J]. Geophysics,1999,64(6):1863-1866.
Carcione J M. A model for seismic velocity and attenuation in petroleum source rocks[J].Geophysics,2000,65:1080-1092.
Carcione J M. Wave fields in real media: wave propagation in anisotropic, anelastic and porousmedia[M]. Pergamon,2001a.
Carcione J M. AVO effects of a hydrocarbon source-rock layer[J]. Geophysics,2001b,66(2):283-293.
Carcione J M. Wave fields in real media: Theory and numerical simulation of wave propagation inanisotropic, anelastic, porous and electromagnetic media,2nd ed.: Elsevier Scientific Publ. Co.,Inc.,2007.
Carcione, J M., and Gurevich, B.,2011, Differential form and numerical implementation of Biot'sporoelasticity equations with squirt dissipation, Geophysics,76, N55-N64.
Carcione J M, Helle H B. Numerical solution of the poroviscoelastic wave equation on a staggeredmesh[J]. J. Comput. Phys.,1999,154:520–527.
Carcione J M, Herman G C, ten Kroode A P E. Seismic modeling[J]. Geophysics,2002,67(4):1304-1325.
Carcione J M, Kosloff D, Kosloff R. Wave propagation in a linear viscoelastic medium[J]. Geophy.J. Int.,1988,95:597-611.
Carcione J M, Picotti S, Gei D, Rossi G. Physics and seismic modeling for monitoring CO2storage[J]. Pure and Applied Geophysics,2006,163:175–207.
Cerjan C, Kosloff D, Kosloff R, et al. A nonreflecting boundary condition for discrete acoustic andelastic wave equation[J]. Geophysics,1985,50(4):705-708.
Chapman M. Frequency dependent anisotropy due to mesoscale fractures in the presence ofequant porosity[J]. Geophysical Prospecting,2003,51:369–379.
Chapman M. Modeling the effect of multiple fracture sets of mesoscale fractures in porous rock onfrequency-dependent anisotropy[J]. Geophysics,2009,74(6): D97–D103.
Chapman M, Liu E, Li X Y. The influence of fluid-sensitive dispersion and attenuation on AVOanalysis[J]. Geophysics,2006,167,89-105.
Chapman M, Maultzsch S, Liu E, et al. The effect of fluid saturation in an anisotropic multi-scaleequant porosity model[J]. Journal of Applied Geophysics,2003,54:191-202.
Chapman M, Zatsepin S V, Crampin S. Derivation of a microstructural poroelastic model[J].Geophysical Journal International,2002,151:427–451.
Du Q Z, Wang X M, Ba J, et al. An equivalent medium model for wave simulation in fracturedporous rocks[J]. Geophysical Prospecting,2012,60(5):940-956.
Duxbury A, White D, Samson C, et al. Fracture mapping using seismic amplitude variation withoffset and azimuth analysis at the Weyburn CO2storage site[J]. Geophysics,2012,77(6):N17-N28.
Dvorkin J, Mavko G. Modeling attenuation in reservoir and non-reservoir rock[J]. The LeadingEdge,2006,25:194–197.
Eshelby J D. The determination of the elastic field of an ellipsoidal inclusion, and relatedproblems[J]. Proceedings of the Royal Society of London,1957(A241):376-396.
Galvin R J, Gurevich B. Scattering of a longitudinal wave by a circular crack in a fluid-saturatedporous medium[J]. International Journal of Solids and Structures,2007,44:7389–7398.
Galvin R J, Gurevich B. Effective properties of a poroelastic medium containing a distribution ofaligned cracks[J]. Journal of Geophysical Research—Solid Earth,2009,114, B07305.
Giroux B. Performance of convolutional perfectly matched layers for pseudospectral time domainporoviscoelastic schemes[J]. Computers&Geosciences,2012,45:149-160.
Gurevich B. Elastic properties of saturated porous rocks with aligned fractures[J]. Journal ofApplied Geophysics,2003,54:203-218.
Gurevich B, Brajanovski R, Galvin R J, et al. P-wave dispersion and attenuation in fractured andporous reservoirs—Poroelasticity approach[J]. Geophysical Prospecting,2009,57(2):225-237.
Hall S A, Kendall J M. Fracture characterization at Valhall: Application of P-wave amplitudevariation with offset and azimuth (AVOA) analysis to a3D ocean-bottom data set[J].Geophysics,2003,68(4):1150–1160.
Hudson J A. Wave speeds and attenuation of elastic waves in material containing cracks[J].Geophys. J. R. Astron. Soc.,1981,64:133–50.
Hudson J A, Liu E, Crampin S. The mechanical properties of materials with interconnected cracksand pores[J]. Geophysical Journal International,1996,124:105–112.
Hudson J A, Pointer T, Liu E. Effective-medium theories for fluid-saturated materials with alignedcracks[J]. Geophysical Prospecting,2001,49:509–522.
Hu Y, McMechan G. Sensitivity of three-component3D finite-difference elastic seismic modelingto inclusion parameters in HTI and TTI media with high inclusion density[J]. Geophysics,2010,75(2): T49-T61.
Jakobsen M. The interacting inclusion model of wave-induced fluid flow[J]. Geophysical JournalInternational,2004,158:1168–1176.
Jakobsen M, Chapman M. Unified theory of global flow and squirt flow in cracked porousmedia[J]. Geophysics,2009,74(2): WA65–WA76.
Jakobsen M, Johansen T A, McCann C. The acoustic signature of fluid flow in complex porousmedia[J]. Journal of Applied Geophysics,2003,54:219–246.
Jilek P. Converted PS-wave reflection coefficients in weakly anisotropic media[J]. Pure andApplied Geophysics,2002,159:1527-1562.
Johansen T A, Jakobsen M, Agersborg R, et al. The P-P and P-S response of a fractured reservoir.74th Annual International Meeting, SEG, Expanded Abstracts,2004,199–202.
Johnson D L. Theory of frequency dependent acoustics in patchy saturated porous media[J].Journal of the Acoustical Society of America,2001,110:682–694.
Klin P, Priolo E, Seriani G. Numerical simulation of seismic wave propagation in realistic3-Dgeo-models with a Fourier pseudo-spectral method[J]. Geophysical Journal International,2010,183:905-922.
Kosloff D, Baysal E. Forward modeling by the Fourier method[J]. Geophysics,1982,47:1402–1412.
Lambert M-A, Saenger E H, Quintal B, et al. Numerical simulation of ambient seismic wavefieldmodification caused by pore-fluid effects in an oil reservoir[J]. Geophysics,2013,78(1):T41–T52.
Liu C, Liu Y W, Feng X, et al. Constructing the convex quadratic function for the evaluation ofcrack density of HTI media using P-and converted waves[J]. J. Geophys. Eng.,2012,9:729–736.
Liu E, Queen J H, Li X Y, et al. Observation and analysis of frequency-dependent anisotropy froma multicomponent VSP at Bluebell-Altamont field[J]. Journal of Applied Geophysics,2003,54:319–333.
Liu J, Ba J, Ma J W, et al. An analysis of seismic attenuation in random porous media[J].Science China—Physics, Mechanics&Astronomy,2010,53:628–637.
Liu L F, Cao S Y, Wang L. Poroelastic analysis of frequency-dependent amplitude-versus-offsetvariations[J]. Geophysics,2011,76(3): C31-C40.
Liu X, Greenhalgh S, Zhou B. Transient solution for poro-viscoacoustic wave propagation indouble porosity media and its limitations[J]. Geophysical Journal International,2009,178:375–393.
Liu X, Greenhalgh S, Zhou B. Approximating the wave moduli of double porosity media at lowfrequencies by a single Zener or Kelvin-Voigt element[J]. Geophysical Journal International,2010,181:391–398.
Maultzsch S, Chapman M, Liu E, et al. Modelling frequency-dependent seismic anisotropy influid-saturated rock with aligned fractures: Implication of fracture size estimation fromanisotropic measurements[J]. Geophysical Prospecting,2003,51:381–392.
Maultzsch S, Chapman M, Liu E, et al. Modelling and analysis of attenuation anisotropy inmulti-azimuth VSP data from the Clair field[J]. Geophysical Prospecting,2007,55:627–642.
Mavko G, Jizba D. Estimating grain-scale fluid effects on velocity dispersion in rocks[J].Geophysics,1991,56:1940–1949.
Mavko G, Mukerji T, Dvorkin J. The rock physics handbook[M].2rd ed. Cambridge UniversityPress,2009.
Muir F, Dellinger J, Etgen J. et al. Modeling elastic fields across irregular boundaries[J].Geophysics,1992,57:1189-1193.
Müller T M, Gurevich B. One-dimensional random patchy saturation model for velocity andattenuation in porous rocks[J]. Geophysics,2004,69:1166–1172.
Müller T M, Gurevich B, Lebedev M. Seismic wave attenuation and dispersion resulting fromwave-induced flow in porous rocks—A review[J]. Geophysics,2010,75(5):75A147-75A164.
Mürphy W F I, Winkler K W, Kleinberg R L. Acoustic relaxation in sedimentary rocks:Dependence on grain contacts and fluid saturation[J]. Geophysics,1986,51:757–766.
Nelson R A. Geological analysis of naturally fractured reservoirs[M].2rd ed. Gulf Publ. Co.,2001.
Parra J O. Poroelastic model to relate seismic wave attenuation and dispersion to permeabilityanisotropy [J]. Geophysics,2000,65(1):202-210.
Picotti S, Carcione J M, Gei D, et al. Seismic modeling to monitor CO2geological storage: TheAtzbach-Schwanenstadt gas field, J. Geophys. Res.,2012,117, B06103.
Picotti S, Carcione J M, Rubino J G, et al. P-wave seismic attenuation by slow-wave diffusion:Numerical experiments in partially saturated rocks[J]. Geophysics,2007,72(4): N11-N21.
Picotti S, Carcione J M, Rubino J G, et al. A visoelastic representation of wave attenuation inporous media[J]. Computers&Geosciences,2010,36:44–53.
Pride S R, Berryman J G, Harris J M. Seismic attenuation due to wave induced flow[J]. Journal ofGeophysical Research,2004,109(B1), B01201.
Pyrak-Nolte L J, Nolte D D. Frequency dependence of fracture stiffness[J]. Geophysical ResearchLetters,1992,19:325–328.
Quintal B, Schmalholz S M, Podladchikov Y Y. Low-frequency reflections from a thin layer withhigh attenuation caused by interlayer flow[J]. Geophysics,2009,74(1): N14–N22.
Quintal B, Steeb H, Frehner M, et al. Quasi-static finite element modeling of seismic attenuationand dispersion due to wave-induced fluid flow in poroelastic media[J]. Journal of GeophysicalResearch Solid Earth,2011,116, B01201.
Ren H T, Goloshubin G, Hilterman F J. Poroelastic analysis of amplitude-versus-frequencyvariations[J]. Geophysics,2009,74(6): N41-N48.
Reshef M, Kosloff D, Edwards M, et al. Three-dimensional acoustic modeling by the Fouriermethod[J]. Geophysics,1988,53:1175-1183.
Robertsson J O A, Blance J O, Symes W W. Viscoelastic finite-difference modeling[J].Geophysics,1994,59:1444-1456.
Rubino J G, Ravazzoli C L, Santos J E. Equivalent viscoelastic solids for heterogeneousfluid-saturated porous rocks[J]. Geophysics,2009,74(1): N1–N13.
Rubino J G, Velis D R, Sacchi M D. Numerical analysis of wave-induced fluid flow effects onseismic data: Application to monitoring of CO2storage at the Sleipner field[J]. J. Geophys. Res.,2011,116, B03306.
Rüger A. P-wave reflection coefficient for transversely isotropic models with vertical andhorizontal axis of symmetry[J]. Geophysics,1997,62:713-722.
Rüger A. Variation of P-wave reflectivity with offset and azimuth in anisotropic media[J].Geophysics,1998,63:935–947.
Ruud B O, Hestholm S. Modeling seismic waves in orthorombic, viscoelastic media by finitedifferences.75th Annual International Meeting, SEG, Expanded Abstracts,2005,1771–1774.
Santos J E, Rubino J G, Ravazzoli C L. A numerical upscaling procedure to estimate effectiveplane wave and shear moduli in heterogeneous fluid-saturated poroelastic media[J]. ComputerMethods in Applied Mechanics and Engineering,2009,198:2067–2077.
Santos J E, Rubino J G, Ravazzoli C L. Modeling mesoscopic attenuation in a highlyheterogeneous Biot’s medium employing an equivalent viscoelastic model.78th AnnualInternational Meeting, SEG, Expanded Abstracts,2008,2112–2116.
Schoenberg M. Elastic wave behavior across linear slip interfaces[J]. J. Acoust. Soc. Am.,1980,68:1516–1521.
Schoenberg M, Sayers C M. Seismic anisotropy of fractured rock[J]. Geophysics,1995,60(1):204-211.
Sil S, Sen M K, Gurevich B. Analysis of fluid substitution in a porous and fractured medium[J].Geophysics,2011,76(3): WA157-WA166.
Thomsen L. Weak elastic anisotropic[J]. Geophysics,1986,51(10):1954-1966.
Thomsen L. Elastic anisotropy due to aligned cracks in porous rock[J]. Geophysical Prospecting,1995,43:805-829.
Tonn R. The determination of the seismic quality Q from VSP data: A comparison of differentcomputational methods[J]. Geophysical Prospecting,1991,39:1-27.
Tsvankin I. Reflection moveout and parameter estimation for horizontal transverse isotropy[J].Geophysics,62(2):614-629.
Van Der Kolk C M, Guest W S, Potters, J H H M. The3D shear experiment over the Natih field inOman: The effect of fracture-filling fluids on shear propagation[J]. Geophysical Prospecting,2001,49:179-197.
Vavrycuk V, Psencik I. PP-wave reflection coefficients in weakly anisotropic elastic media[J].Geophysics,1998,63:2129-2141.
White J E. Computed seismic speeds and attenuation in rocks with partial gas saturation[J].Geophysics,1975,40:224–232.
White J E, Mikhaylova N G, Lyakhovitskiy F M. Low-frequency seismic waves in fluid-saturatedlayered rocks[J]. Physics of the Solid Earth,1976,11:654–659.
Winkler K. Dispersion analysis of velocity and attenuation in Berea sandstone[J]. Journal ofGeophysical Research,1985,90:6793–6800.
Witte D, Richards P G. Contribution to the pseudo-spectral method for computing syntheticseismograms. Expanded Abstracts of57th SEG Annual Int. Mtg,1987,517-519.
Zhang S Q, Han L G, Liu C C, et al. Inverting reservoir parameters in a two-phase fracturedmedium with a niche genetic algorithm[J]. Applied Geophysics,2012,9(4):440-450.
巴晶，曹宏，姚逢昌等.利用错格虚谱法模拟礁、滩相储层地震波场[J].石油地球物理勘探,2010,45(2):177-184.
巴晶，卢明辉，胡彬等.黏弹双相介质中的松弛骨架模型[J].地球物理学报,2008,51(5):1527-1537.
杜启振.各向异性黏弹性介质伪谱法波场模拟[J]，物理学报，2004,53(12):4428-4434.
杜启振，孔丽云，韩世春.裂缝诱导各向异性双孔隙介质波场传播特征[J].地球物理学报，2009,52(4):1049-1058.
杜启振，刘莲莲，孙晶波.各向异性粘弹性孔隙介质地震波场伪谱法正演模拟[J].物理学报，2007，56(10):6143-6149.
郭智奇.粘弹各向异性介质波场模拟与储层信息研究[D].长春：吉林大学，2008.
郭智奇，刘财，杨宝俊等.粘弹各向异性介质中地震波场模拟与特征[J].地球物理学进展，2007,22(3):804-810.
孔丽云，王一博，杨慧珠.裂缝诱导HTI双孔隙介质中的裂缝参数分析[J].地球物理学报，2012,55(1):189-196.
雷英杰，张善文，李续武等. Matlab遗传算法工具箱及应用[M].西安：西安电子科技大学出版社，2005.
李敏强，寇纪淞，林丹等.遗传算法的基本理论与应用[M].北京：科学出版社，2002.
刘财，郭智奇，杨宝俊等.粘弹各向异性介质中波的反射与透射问题分析[J].地球物理学报，2007,50(4):1216-1224.
刘财，兰慧田，郭智奇等．基于改进BISQ机制的双相HTI介质波传播伪谱法模拟与特征分析[J]．地球物理学报，2013，56（10）：3461-3473.
刘财，刘宇巍，冯晅等.基于方位相交的纵波AVA数据运用SVD反演HIT介质裂缝密度[J].吉林大学学报（地球科学版），2013,43(5)：1655-1662.
刘炯,马坚伟,杨慧珠等.缝洞型储层错格伪谱法模拟[J].石油地球物理勘探,2008,43(6):723-727.
刘洋，李承楚.双相各向异性介质中弹性波传播伪谱法数值模拟研究.地震学报，2000,22(2):132-138.
聂建新,杨顶辉,杨慧珠.基于非饱和多孔隙介质BISQ模型的储层参数反演[J].地球物理学报,2004,47(6):1101-1105.
石耀霖，金文.面波频散反演地球内部构造的遗传算法[J].地球物理学报，1995，38(2):189-198.
轩义华，何樵登，孟庆生等.基于BISQ机制的双相EDA介质的波场分析[J].石油地球物理勘探，2006,41(5):550-556.
张广智,陈怀震,印兴耀等.基于各向异性AVO的裂缝弹性参数叠前反演方法[J].吉林大学学报（地球科学版），2012,42(3):845-851.
张军舵.孔隙流体介质波动方程正演模拟与参数反演[D].北京：中国石油大学，2008.
曾联波，柯式镇，刘洋.低渗透油气[M].北京：石油工业出版社，2010.
周明，孙树栋.遗传算法原理及应用[M].北京：国防工业出版社，1999.
朱培民,王家映,於文辉等.用纵波AVO数据反演储层裂隙密度参数[J].石油物探，2001,40(2)：1-12.
张文生，何樵登.二维横向各向同性介质的伪谱法正演模拟[J].石油地球物理勘探,1998,33(3):310-319.
张显文，王德利，王者江等.基于BISQ机制三维双相正交裂隙各向异性介质衰减及频散方位特性研究[J].地球物理学报，2010,53(10):2452-2459.
张显文.储层地震预测技术与储层地球物理参数反演[D].长春：吉林大学，2010.
赵志新,徐纪人,堀内茂木.错格实数傅立叶变换微分算子及其在非均匀介质波动传播研究中的应用[J].地球物理学报,2003,46(2):234-240.