弱各向异性介质转换波AVO分析与流体识别方法研究
|
摘要
应用振幅随偏移距变化(AVO)等叠前地震信息预测含油气储层已得到了广泛的实践,在国内外油气地球物理勘探与开发领域中取得了大量成功实例,但由于影响地震波振幅变化的因素非常复杂,常规的AVO分析也存在着多解性问题,解决这一多解性的有效途径是利用多波多分量地震勘探信息。纵波和P-SV转换波的联合AVO分析及应用可以大大减少多解性,这点已得到了地球物理工作者的普遍认同；同时,综合二者的AVO属性还可以获得更多更准确的储层岩石物理参数,因此,转换波AVO技术的研究与完善是储层地震勘探亟待解决的课题之一。
科学合理地构建地质-地球物理模型也是储层预测和流体识别的关键问题。从南海深水区等油气藏实际数据分析来看,基于双相介质和各向异性介质的储层建模比常规的各向同性介质储层模型更加科学合理。为了更好地利用叠前各种地震信息来精确预测储层和识别流体,人们必须首先研究各向异性等复杂介质中的纵波及转换波的AVO的理论及分析技术,并据此来获取储层的岩性及空间展布信息等。本文从实效性就弱各向异性介质中转换波的反射系数公式进行了系统研究,为了提高计算效率推导了简化公式,建立了各种属性参数计算方法,同时推导了弱各向异性介质中P-SV波弹性阻抗,构建了各向异性介质中广义流体识别因子,得出了一系列有意义的研究成果。
本文首先给出了各向同性介质中转换波的反射系数线性近似公式,系统分析了基于该公式的AVO分析算子,认为AVO分析的不同近似式对模型都有不同程度的依赖性。根据四种典型的含气含水层实例分析,P-SV波的AVO响应特征较PP波要复杂,不能同时区分四类含气含水层,尤其是对第一类和第二类基本没有识别效果,对第三类的识别效果较好,说明了转换波的AVO对储层的流体成分不是很敏感。而通过归一化后的反射系数AVA响应曲线呈单调增加或减小同时还放大了异常,对四类AVO类型的含气含水层识别效果都有所提高,对前三类,归一化后的反射系数随着入射角的增加呈单调递增现象,对第四类,反射系数是呈单调递减的,但是反射系数均为正值。
最为流行的各向异性介质中转换波反射系数近似公式是Vavrycuk(1999)给出的,但该公式在AVO分析中的适用性不强,因此Petr Jilek(2002)给出了PS波反射系数近似公式的另一种表达形式,适用于任意半空间各向异性介质,由于假设了界面两侧纵横波速度和密度的弱差异性和任意半空间的弱各向异性,所以称此公式为弱差异弱各向异性近似公式,它是关于Thomsen各向异性参数、入射角和方位角的函数,并适用于以下介质组合：各向同性、VTI、HTI和正交各向异性,垂直对称平面的任意方位角适用于HTI介质和正交介质。但是Petr Jilek的表达式比较复杂,因此我们对正交各向异性介质中P-SV波反射系数公式进行了简化,对各向异性参数进行了转换,将其退化到VTI介质和HTI介质两种常用的典型各向异性介质中,结合各向同性中转换波反射系数的近似方式,尝试将反射系数表达成sinθ奇函数的形式,在保证精度的基础上分别保留到sinθ的一次幂、三次幂、五次幂和七次幂,本文将其称之为一阶、三阶、五阶、七阶新近似公式,均为各向同性反射系数与各向异性反射系数扰动项加和的形式,通过不同的理论模型对比了几种表达形式的精度,并做了相应的误差分析,同时针对不同模型做了反射系数的影响因素分析,得出了相应的规律。
在得出新近似公式之后,本文首次尝试将归一化反射系数思想用于VTI和HTI介质的P-SV波反射系数三阶新近似公式,利用各向同性情况下四类经典AVO含气含水砂岩模型的具体参数值,在此基础上增加了Thomsen各向异性参数,计算并绘制含气含水砂岩AVA的响应特征,获得了与各向同性介质情况下类似的结果。仿照各向同性条件下的P-SV波AVO截距梯度属性分析方法,以VTI介质为例讨论了P-SV波AVO截距梯度属性,计算结果表明,各向异性介质转换波截距梯度的正负规律跟各向同性介质情况下的规律是类似的,由各向异性参数引起的区别在于截距和梯度取值的大小。不同储层模型的交会图的特征是不同的,但是它们的背景趋势可以拟合成一条直线,由直线的斜率可以求取平均横纵波速度比。在进行AVO属性交会图分析之后可以根据需要得到横波速度剖面,密度剖面和剪切模量剖面,为AVO定性分析提供了有利的途径,在横纵波速度比k一定的情况下对各向异性参数可以进行定量分析,利用该公式也可以进行参数反演。
各向同性储层模型弹性阻抗公式也具有多种形式,并且大部分已经应用于叠前实际数据处理中。而对转换波弹性阻抗以及各向异性介质中弹性阻抗的研究及成果都比较少,在研究过程中需借鉴各向同性和各向异性介质中纵波弹性阻抗的方法和理论。因此,利用本文给出的VTI和HTI介质转换波反射系数近似公式推导出了最新的转换波弹性阻抗公式,通过不同模型的正演模拟,对公式精度进行了误差分析,并结合whitcombe的归一化思想,形成了弱各向异性介质中扩展的转换波弹性阻抗公式,最后对影响弹性阻抗公式的主要因素及规律做了总结。通过弹性阻抗的交会分析,也表明了在储层识别方面弹性阻抗比声波阻抗更灵敏。
为了从实际地震资料中提取更多地反映储层流体的信息,人们提出了许多流体因子,地震纵波、横波和转换波的综合利用为此也提供了一条有效的途径。这些方法都是使用线性Zoeppritz方程来提取物性参数,例如P波和S波阻抗、速度、泊松比等,同时也可以获得一些利于解释的参数,如弹性模量、拉梅常数以及一些流体因子。通常这些流体因子都可以较好地识别烃类区域。从流体因子建立的岩石物理基础Biot-Gassmann方程出发,分析了填充不同类型的流体时岩石对速度和波阻抗的影响,在此基础上列举了常见的阻抗类流体识别因子,使用四类典型的砂岩模型做数值试验,总结了各个流体因子的对含气含水砂岩的识别能力及敏感性。为了将各向异性参数、入射角和方位角引入到流体识别因子公式中,本文尝试把弱各向异性介质中PP波和P-SV波弹性波阻抗公式代入到已有的流体因子公式中,构成新的广义流体因子公式,这里已有的流体识别因子选择的是高灵敏度流体识别因子Fdanj和灵敏度较弱的μρani识别公式,便于进行灵敏度比较,分析了流体因子Fdani和μρani对三种岩性组合模型的识别效果,同时也分析了Thomsen各向异性参数的影响,得出了相应规律。
转换波信息是纵波的补充,同时AVO、弹性阻抗和流体识别都是储层预测的关键技术,在加入各向异性因素之后,可以提高它们对储层预测和流体识别的精度。
Application of AVO and other pre-stack seismic data to predict oil and gas reservoirs has been widely used, which has made a lot of success stories in domestic and international oil and gas geophysical exploration and development. But due to the complex factors of AVO changes, conventional AVO analysis also exists multiple solutions, and multi-wave and multi-component seismic exploration information is an effective way to solve multiple solutions. The combination using of converted wave would greatly reduce multiple solutions, which is generally accepted by geophysicists. At the same time, comprehensive AVO attributes can get more accurate physical parameters of reservoir, therefore, the research and improvement of converted wave AVO technology is an urgent subject of reservoir seismic exploration.
Building the geological-geophysical model scientifically and reasonably is the key issue for reservoir prediction and fluid identification. From the actual oil and gas reservior data analysis of the South China Sea deep water, the reservoir modeling based on two-phase media and anisotropic media is more scientific and reasonable than the conventional isotropic media. In order to make better use of pre-stack seismic information to accurately predict reservoir and identify fluid, we must first study PP wave and P-SV wave AVO in anisotropy media, according to which we can get lithology information and spatial distribution of reservoir. In this paper, the weak anisotropic media converted waves linear approximation formula is simplified, and AVO attribute analysis is done, while the weak anisotropic media P-SV wave elastic impedance is derived, the generalized anisotropic medium fluid identification factors are constructed, which are very meaningful.
The paper first gives the advantage of converted wave AVO analysis, and summarizes converted wave reflection coefficient linear approximate formula in isotropic medium, recognizing that even with the same formula, accuracies of different models are not the same, which means that approximate formulas have some dependence on the model. For the four typical gas aquifer, P-SV wave AVA response characteristics are more complex than PP wave, and cannot distinguish the four types of gas-bearing aquifer, particularly did not identify first and second base, has better identification of the third category, indicating that the converted wave AVO is not very sensitive to the fluid composition of the reservoir. Reflection coefficient AVA response curve after normalized is monotonic increase or decrease, but also enlarge the exception, recognition effect for four types of AVO gas aquifer has increased, for former three categories, normalized reflection coefficient decreases with the increasing angle, and for the fourth category, the reflection coefficient is monotonically decreasing, but the reflection coefficient is positive.
The most common converted wave reflection coefficient approximate formula in anisotropic media is given in 1999 by Vavrycuk, but the applicability of formula in the actual of AVO analysis is not strong. Therefore. Petr Jilek (2002) gives another PS wave reflection coefficient approximate formula expression for arbitrary anisotropic dielectric half-space, the derivation of the formula is based on the following assumptions:weak contrasts of elastic medium parameters across the interface and weak anisotropy in half-space, so this formula is called weak contrasts and weak anisotropy approximations, this formula is the function of the Thomsen anisotropic parameters, incident angle and azimuth. This resulting expressions can be applied for any combination of isotropic. VTI. HTI and orthorhombic halfspaces.. But it is complex, so we simplified orthotropic medium P-SV wave reflection coefficient formula, converted the anisotropy parameters, and degradated it to the VTI and HTI media, which are two typical anisotropic media.we commonly used it in combination with converted wave reflection coefficient in the isotropic approximation method, try to express the reflection coefficient as odd functions of sinθ.on the base of accuracy, we reserved it to one power.three power, five power and seven power, the paper called it a first order approximation, the third order, fifth-order, seven-order new approximate formula. Both are the combined-items of isotropic reflection coefficient and anisotropic disturbed item, through different models we compared the accuracy of several forms of expression, and made the corresponding error analysis, at the same time for different models the impact of factors of the reflection coefficient obtained corresponding laws.
After obtaining the new approximate formula, we firstly attempt to normalize the P-SV wave new third-order approximation reflection coefficient in VTI and HTI media. Based on using of four classic case of isotropic AVO gas sandstone model parameters, adding to the anisotropy factors-Thomsen anisotropic parameters of the assignment, we draw gas/water sandstone response characteristics of AVA to get the effect of identification, obtained the similar results under the isotropic medium. Similar to the isotropic P-SV wave AVO gradient and intercept method, an example of P-SV AVO gradient attribute intercept is discussed in the VTI medium wave, and we obtain anisotropic converted wave positive and negative gradient case law with the law of isotropic media is similar, while the difference caused by the anisotropic parameter is the size of the intercept and gradient values. Forms of cross plots of different models are different, but their background trend can fit a straight line. The slope of the line can strike the average k value. After the research on AVO attributes analysis by new equation,we can get S wave velocity profile, density profile and shear modulus profile, which can provide a favorable way for qualitative analysis. When attribution parameters are constant, quantitative analysis and inversion of anisotropic parameters can be carried out.
Many scholars have proposed different types of elastic impedance formula, the P-wave elastic impedance inversion of pre-stack data has become more sophisticated, while converted wave elastic impedance and elastic impedance in anisotropic media research is less than P-wave, so it needs to draw on the research based on the isotropic elastic impedance method and theory. Based on this, P-SV wave elastic impedance formulations in VTI and HTI were proposed by approximation formula.
Then we present shale/gas sand of models to test the accuracy of P-SV wave elastic impedance.We normalized it with whitcombe's ideas, and formed converted wave elastic impedance formula in weak anisotropic media, and elastic impedance of the final formula of the main factors and rules are also summarized. The intersection of elastic impedance analysis also shows that recognition of reservoir types in terms of elastic impedance is more sensitive than the acoustic impedance.
In order to predict reservoir by using AVO and multi-components seismic data, some fluid factors and fluid identification methods have been studied, seismic wave, shear wave and converted wave of comprehensive utilization of this also provides an effective way. These methods are used to extract linear Zoeppritz equation parameters, such as P wave and S wave impedance, velocity, Poisson ratio, but also can get some help to explain the parameters, such as elastic modulus. Lame constant, and some fluid factor. Generally, these fluids can better identify hydrocarbon region. Starting from the physical basis of Biot-Gassmann equations, we analyze the impact of filling the rock with different types of fluids on velocity and impedance, and in this basis list the common impedance fluid identification factors, which are used to identify four typical sandstones. Here recognition ability and sensitivity of the various factors on gas/water sandstone are analyzed. In order to introduce the anisotropic parameters and azimuthal angle into the equation of fluid factor, we insert the PP and P-SV wave elastic impedance into the fluid factor. The high sensitive fluid identification Fdam and low sensitive fluid identificationμpam are compared, which is easy to compare the sensitivity analysis of the fluid factor Fdam andμρam. Though three lithological model identification, we also analyze the Thomsen anisotropic parameters'effect and obtain the corresponding laws.
Converted wave is the supplement for P-wave information, and at the same time AVO. elastic impedance and fluid identification are the key technology,joining the anisotropy factors for reservoir prediction, which can improve the accuracy of their reservoir interpretation, so our study is very meaningflul.
科学合理地构建地质-地球物理模型也是储层预测和流体识别的关键问题。从南海深水区等油气藏实际数据分析来看,基于双相介质和各向异性介质的储层建模比常规的各向同性介质储层模型更加科学合理。为了更好地利用叠前各种地震信息来精确预测储层和识别流体,人们必须首先研究各向异性等复杂介质中的纵波及转换波的AVO的理论及分析技术,并据此来获取储层的岩性及空间展布信息等。本文从实效性就弱各向异性介质中转换波的反射系数公式进行了系统研究,为了提高计算效率推导了简化公式,建立了各种属性参数计算方法,同时推导了弱各向异性介质中P-SV波弹性阻抗,构建了各向异性介质中广义流体识别因子,得出了一系列有意义的研究成果。
本文首先给出了各向同性介质中转换波的反射系数线性近似公式,系统分析了基于该公式的AVO分析算子,认为AVO分析的不同近似式对模型都有不同程度的依赖性。根据四种典型的含气含水层实例分析,P-SV波的AVO响应特征较PP波要复杂,不能同时区分四类含气含水层,尤其是对第一类和第二类基本没有识别效果,对第三类的识别效果较好,说明了转换波的AVO对储层的流体成分不是很敏感。而通过归一化后的反射系数AVA响应曲线呈单调增加或减小同时还放大了异常,对四类AVO类型的含气含水层识别效果都有所提高,对前三类,归一化后的反射系数随着入射角的增加呈单调递增现象,对第四类,反射系数是呈单调递减的,但是反射系数均为正值。
最为流行的各向异性介质中转换波反射系数近似公式是Vavrycuk(1999)给出的,但该公式在AVO分析中的适用性不强,因此Petr Jilek(2002)给出了PS波反射系数近似公式的另一种表达形式,适用于任意半空间各向异性介质,由于假设了界面两侧纵横波速度和密度的弱差异性和任意半空间的弱各向异性,所以称此公式为弱差异弱各向异性近似公式,它是关于Thomsen各向异性参数、入射角和方位角的函数,并适用于以下介质组合：各向同性、VTI、HTI和正交各向异性,垂直对称平面的任意方位角适用于HTI介质和正交介质。但是Petr Jilek的表达式比较复杂,因此我们对正交各向异性介质中P-SV波反射系数公式进行了简化,对各向异性参数进行了转换,将其退化到VTI介质和HTI介质两种常用的典型各向异性介质中,结合各向同性中转换波反射系数的近似方式,尝试将反射系数表达成sinθ奇函数的形式,在保证精度的基础上分别保留到sinθ的一次幂、三次幂、五次幂和七次幂,本文将其称之为一阶、三阶、五阶、七阶新近似公式,均为各向同性反射系数与各向异性反射系数扰动项加和的形式,通过不同的理论模型对比了几种表达形式的精度,并做了相应的误差分析,同时针对不同模型做了反射系数的影响因素分析,得出了相应的规律。
在得出新近似公式之后,本文首次尝试将归一化反射系数思想用于VTI和HTI介质的P-SV波反射系数三阶新近似公式,利用各向同性情况下四类经典AVO含气含水砂岩模型的具体参数值,在此基础上增加了Thomsen各向异性参数,计算并绘制含气含水砂岩AVA的响应特征,获得了与各向同性介质情况下类似的结果。仿照各向同性条件下的P-SV波AVO截距梯度属性分析方法,以VTI介质为例讨论了P-SV波AVO截距梯度属性,计算结果表明,各向异性介质转换波截距梯度的正负规律跟各向同性介质情况下的规律是类似的,由各向异性参数引起的区别在于截距和梯度取值的大小。不同储层模型的交会图的特征是不同的,但是它们的背景趋势可以拟合成一条直线,由直线的斜率可以求取平均横纵波速度比。在进行AVO属性交会图分析之后可以根据需要得到横波速度剖面,密度剖面和剪切模量剖面,为AVO定性分析提供了有利的途径,在横纵波速度比k一定的情况下对各向异性参数可以进行定量分析,利用该公式也可以进行参数反演。
各向同性储层模型弹性阻抗公式也具有多种形式,并且大部分已经应用于叠前实际数据处理中。而对转换波弹性阻抗以及各向异性介质中弹性阻抗的研究及成果都比较少,在研究过程中需借鉴各向同性和各向异性介质中纵波弹性阻抗的方法和理论。因此,利用本文给出的VTI和HTI介质转换波反射系数近似公式推导出了最新的转换波弹性阻抗公式,通过不同模型的正演模拟,对公式精度进行了误差分析,并结合whitcombe的归一化思想,形成了弱各向异性介质中扩展的转换波弹性阻抗公式,最后对影响弹性阻抗公式的主要因素及规律做了总结。通过弹性阻抗的交会分析,也表明了在储层识别方面弹性阻抗比声波阻抗更灵敏。
为了从实际地震资料中提取更多地反映储层流体的信息,人们提出了许多流体因子,地震纵波、横波和转换波的综合利用为此也提供了一条有效的途径。这些方法都是使用线性Zoeppritz方程来提取物性参数,例如P波和S波阻抗、速度、泊松比等,同时也可以获得一些利于解释的参数,如弹性模量、拉梅常数以及一些流体因子。通常这些流体因子都可以较好地识别烃类区域。从流体因子建立的岩石物理基础Biot-Gassmann方程出发,分析了填充不同类型的流体时岩石对速度和波阻抗的影响,在此基础上列举了常见的阻抗类流体识别因子,使用四类典型的砂岩模型做数值试验,总结了各个流体因子的对含气含水砂岩的识别能力及敏感性。为了将各向异性参数、入射角和方位角引入到流体识别因子公式中,本文尝试把弱各向异性介质中PP波和P-SV波弹性波阻抗公式代入到已有的流体因子公式中,构成新的广义流体因子公式,这里已有的流体识别因子选择的是高灵敏度流体识别因子Fdanj和灵敏度较弱的μρani识别公式,便于进行灵敏度比较,分析了流体因子Fdani和μρani对三种岩性组合模型的识别效果,同时也分析了Thomsen各向异性参数的影响,得出了相应规律。
转换波信息是纵波的补充,同时AVO、弹性阻抗和流体识别都是储层预测的关键技术,在加入各向异性因素之后,可以提高它们对储层预测和流体识别的精度。
Application of AVO and other pre-stack seismic data to predict oil and gas reservoirs has been widely used, which has made a lot of success stories in domestic and international oil and gas geophysical exploration and development. But due to the complex factors of AVO changes, conventional AVO analysis also exists multiple solutions, and multi-wave and multi-component seismic exploration information is an effective way to solve multiple solutions. The combination using of converted wave would greatly reduce multiple solutions, which is generally accepted by geophysicists. At the same time, comprehensive AVO attributes can get more accurate physical parameters of reservoir, therefore, the research and improvement of converted wave AVO technology is an urgent subject of reservoir seismic exploration.
Building the geological-geophysical model scientifically and reasonably is the key issue for reservoir prediction and fluid identification. From the actual oil and gas reservior data analysis of the South China Sea deep water, the reservoir modeling based on two-phase media and anisotropic media is more scientific and reasonable than the conventional isotropic media. In order to make better use of pre-stack seismic information to accurately predict reservoir and identify fluid, we must first study PP wave and P-SV wave AVO in anisotropy media, according to which we can get lithology information and spatial distribution of reservoir. In this paper, the weak anisotropic media converted waves linear approximation formula is simplified, and AVO attribute analysis is done, while the weak anisotropic media P-SV wave elastic impedance is derived, the generalized anisotropic medium fluid identification factors are constructed, which are very meaningful.
The paper first gives the advantage of converted wave AVO analysis, and summarizes converted wave reflection coefficient linear approximate formula in isotropic medium, recognizing that even with the same formula, accuracies of different models are not the same, which means that approximate formulas have some dependence on the model. For the four typical gas aquifer, P-SV wave AVA response characteristics are more complex than PP wave, and cannot distinguish the four types of gas-bearing aquifer, particularly did not identify first and second base, has better identification of the third category, indicating that the converted wave AVO is not very sensitive to the fluid composition of the reservoir. Reflection coefficient AVA response curve after normalized is monotonic increase or decrease, but also enlarge the exception, recognition effect for four types of AVO gas aquifer has increased, for former three categories, normalized reflection coefficient decreases with the increasing angle, and for the fourth category, the reflection coefficient is monotonically decreasing, but the reflection coefficient is positive.
The most common converted wave reflection coefficient approximate formula in anisotropic media is given in 1999 by Vavrycuk, but the applicability of formula in the actual of AVO analysis is not strong. Therefore. Petr Jilek (2002) gives another PS wave reflection coefficient approximate formula expression for arbitrary anisotropic dielectric half-space, the derivation of the formula is based on the following assumptions:weak contrasts of elastic medium parameters across the interface and weak anisotropy in half-space, so this formula is called weak contrasts and weak anisotropy approximations, this formula is the function of the Thomsen anisotropic parameters, incident angle and azimuth. This resulting expressions can be applied for any combination of isotropic. VTI. HTI and orthorhombic halfspaces.. But it is complex, so we simplified orthotropic medium P-SV wave reflection coefficient formula, converted the anisotropy parameters, and degradated it to the VTI and HTI media, which are two typical anisotropic media.we commonly used it in combination with converted wave reflection coefficient in the isotropic approximation method, try to express the reflection coefficient as odd functions of sinθ.on the base of accuracy, we reserved it to one power.three power, five power and seven power, the paper called it a first order approximation, the third order, fifth-order, seven-order new approximate formula. Both are the combined-items of isotropic reflection coefficient and anisotropic disturbed item, through different models we compared the accuracy of several forms of expression, and made the corresponding error analysis, at the same time for different models the impact of factors of the reflection coefficient obtained corresponding laws.
After obtaining the new approximate formula, we firstly attempt to normalize the P-SV wave new third-order approximation reflection coefficient in VTI and HTI media. Based on using of four classic case of isotropic AVO gas sandstone model parameters, adding to the anisotropy factors-Thomsen anisotropic parameters of the assignment, we draw gas/water sandstone response characteristics of AVA to get the effect of identification, obtained the similar results under the isotropic medium. Similar to the isotropic P-SV wave AVO gradient and intercept method, an example of P-SV AVO gradient attribute intercept is discussed in the VTI medium wave, and we obtain anisotropic converted wave positive and negative gradient case law with the law of isotropic media is similar, while the difference caused by the anisotropic parameter is the size of the intercept and gradient values. Forms of cross plots of different models are different, but their background trend can fit a straight line. The slope of the line can strike the average k value. After the research on AVO attributes analysis by new equation,we can get S wave velocity profile, density profile and shear modulus profile, which can provide a favorable way for qualitative analysis. When attribution parameters are constant, quantitative analysis and inversion of anisotropic parameters can be carried out.
Many scholars have proposed different types of elastic impedance formula, the P-wave elastic impedance inversion of pre-stack data has become more sophisticated, while converted wave elastic impedance and elastic impedance in anisotropic media research is less than P-wave, so it needs to draw on the research based on the isotropic elastic impedance method and theory. Based on this, P-SV wave elastic impedance formulations in VTI and HTI were proposed by approximation formula.
Then we present shale/gas sand of models to test the accuracy of P-SV wave elastic impedance.We normalized it with whitcombe's ideas, and formed converted wave elastic impedance formula in weak anisotropic media, and elastic impedance of the final formula of the main factors and rules are also summarized. The intersection of elastic impedance analysis also shows that recognition of reservoir types in terms of elastic impedance is more sensitive than the acoustic impedance.
In order to predict reservoir by using AVO and multi-components seismic data, some fluid factors and fluid identification methods have been studied, seismic wave, shear wave and converted wave of comprehensive utilization of this also provides an effective way. These methods are used to extract linear Zoeppritz equation parameters, such as P wave and S wave impedance, velocity, Poisson ratio, but also can get some help to explain the parameters, such as elastic modulus. Lame constant, and some fluid factor. Generally, these fluids can better identify hydrocarbon region. Starting from the physical basis of Biot-Gassmann equations, we analyze the impact of filling the rock with different types of fluids on velocity and impedance, and in this basis list the common impedance fluid identification factors, which are used to identify four typical sandstones. Here recognition ability and sensitivity of the various factors on gas/water sandstone are analyzed. In order to introduce the anisotropic parameters and azimuthal angle into the equation of fluid factor, we insert the PP and P-SV wave elastic impedance into the fluid factor. The high sensitive fluid identification Fdam and low sensitive fluid identificationμpam are compared, which is easy to compare the sensitivity analysis of the fluid factor Fdam andμρam. Though three lithological model identification, we also analyze the Thomsen anisotropic parameters'effect and obtain the corresponding laws.
Converted wave is the supplement for P-wave information, and at the same time AVO. elastic impedance and fluid identification are the key technology,joining the anisotropy factors for reservoir prediction, which can improve the accuracy of their reservoir interpretation, so our study is very meaningflul.
引文
[1]Gaiser J., Moldoveanu N., Macbeth.et al. Multicomponent technology:the players,problems,application,and trends[J]. The Leading Edge,2001,20(9): 974-977.
[2]Crampin S., Anisotropy in exploration seismics[J].,First Break,1984b, Vol.23, 19-21.
[3]Crampin,S. Evaluation of anisotropy by shear-wave splitting[J]. Geophysics, 1985.50:142-152.
[4]Crampin S. Suggestios for a consistent terminology for seismic anisotropy [J]. Geophysical Prospecting.1989,37(7):753-770.
[5]Stewart,R.R.,&Gaiser J.E. Application and interpretation of converted waves[J]. SEG Continuing Education Short Course.
[6]黄中玉.多分量地震勘探的机遇和挑战[J].石油物探,2001,40(2)：131-137.
[7]石玉梅,姚逢昌,曹宏.多波多分量天然气勘探技术的进展[J].勘探地球物理进展,2003,26(3)：172-177.
[8]陈天胜,刘洋,魏修成.纵波和转换波联合AVO反演方法研究[J].中国石油人学学报(自然科学版),2006,30(1)：33-37.
[9]Feng Hong. Hydrocarbon indicators derived from AVO analysis[D]. Ph.D. thesis,University of Calgary,2009.
[10]Ramos A.C.B, Castagna J. P, Useful approximations for converted-wave AVO[J]. Geophysics,2001,66(6):1721-1734.
[11]Hudson J A. Overall properties of a cracked solid [J]. Math Proc Camb phil Soc.1980,88:371-384.
[12]Xu X, Tsvankin I. Azimuthal AVO analysis with anisotropic spreading correction: A synthetic study [J]. The Leading Edge,2006,25(11),1336-1342.
[13]Xu X, Tsvankin I. A case study of azimuthal AVO analysis with anisotropic spreading correction[J]. The Leading Edge,2007,26(12),1552-1561.
[14]Connoly P. Elastic impedance[J]. The Leading Eage,1999,18(4):438-452.
[15]曹孟起,王九栓,邵林海.叠前弹性波阻抗反演技术及应用[J].石油地球物理勘探,2006,4(3)：323-326.
[16]王保丽,印兴耀,张繁昌.弹性阻抗反演及应用研究[J].地球物理学进展,2005,20(1)：89-92.
[17]倪逸.弹性阻抗计算的一种新方法[J].石油地球物理勘探,2003,38(2)：147-155.
[18]王保丽.叠前弹性阻抗反演理论与应用[D].东营：石油大学(华东),2005.
[19]邹文.基于地震资料的流体识别技术[D].成都：成都理工大学,2008.
[20]Stewart R. R., Gaiser J. E., Brown R. J.,and Lawton D. C. Converted-wave seismic exploration:Applications[J].Geophysics.2003.68,40-57.
[21]Larsen J.,Margrave G. and Lu, H. X. AVO analysis by simultaneous P-P and P-S weighted stacking applied to 3-C-3-D seismic data[J].69th Ann. Internat. Mtg,. Soc. of Expl. Geophys., Expanded Abstracts.1999,721-724.
[22]Kelly M. Skidmore C.and Cotton R.. P-P and P-S angle stack inversionfJ].69th Ann. Internat. Mtg,. Soc. of Expl. Geophys., Expanded Abstracts.2000.222-223.
[23]Margrave G. F.. Stewart R. R. and Larsen J. A. Joint PP and PS seismic inversion[J].The Leading Edge,2001,20,1048-1052.
[24]Veire H..and Landro M. Jiont inversion of PP-and PS-seismic data[J].71st Ann. Internat. Mtg.. Soc. of Expl. Geophys., Expanded Abstracts.2001.861-864.
[25]Engelmark F. Using converted shear waves to image reservoirs with low-inpedance contrast[J]. The Leading Edge.2000,19:600-603.
[26]Gonzalez E. F.,Mukerji. T. and Mavko G.,Facies identification using P-to-P and P-to-S AVO attributes [J].70th Ann. Internat. Mtg, Soc of Expl Geophys. Expanded Abstracts,2000:98-101.
[27]Jin S.,Cambois G. and Vuillermoz C. Shear-wave velocity and density estimation from PS-wave AVO analysis:Application to OBS dataset from the North Sea[J].Geophysics,2000,65:1446-1454.
[28]Wu Y. Estimation of gas saturation using P-to-S converted waves s[J].70th Ann. Internat. Mtg. Soc of Expl Geophys. Expanded Abstracts,2000:158-161.
[29]Zhu F., Gibson R. L.. Atkins J.,and Yuh S. H.. Distingrishing fizz gas from commercial gas reservoirs using multicomponent seismic data[J]. The Leading Edge.19,1238-1245.
[30]Bortfeld R. Approximation to the reflection and transmission coefficients of plane longitudinal and transverse waves [J]. Geophysical Prospecting,1961,9:485-502.
[31]Aki K.,Richards P. G.. Quantitative Seismology:Theory and Methods[M]. San Francisco:Freenman and Coopration.1980,146-177.
[32]郑晓东.Zoeppritz方程的近似及应用[J].石油地球物理勘探,1991,26(2)：129-144
[33]郑晓东.AVO理论和方法的一些新进展[J].石油地球物理勘探,1992,27(3)： 305-317.
[34]周竹生.P-SV波和SH波的AVO分析[J].石油地球物理勘探,1993,28(4)：430-438
[35]杨绍国,周熙襄.Zoeppritz方程的级数表达式及近似[J].石油地球物理勘探,1994,29(4)：399-412.
[36]甘永忠.用P-SV波叠加剖面合成横波层速度剖面[J].石油地球物理勘探1995.30(1)：56-61.
[37]Donati M., Martin N W. A comparison of approximations for the converted wave reflection[R]. Crewes Research Reports. Calgary,1998,159-179.
[38]Xu Y., Bancroft J C. Joint AVO analysis of PP and PS data[R]. Crewes Research Reports 34.1997,1-44.
[39]Xu Y.,AVO developments applied to Blackfoot 3C-2D broadband line[D]. Calgary:University of Calgary,1999.
[40]Wang Y H., Approximations to Zoeppritz equations and their use in AVO analysis [J].Geophysics,1999,64(6):1920-1927.
[41]Ezequiel F G, Tapan M, Gary M. Facies classification using P-to-P and P-to-S AVO attributes[J]. Expanded Abstracts of 70th Annual Internat SEG Mtg,2000.98-101
[42]Vant A, Brown R. J, Approximation the P-S reflection coefficient for small incidence angles[J]. C rewes Research Reports 13,2001,159-179.
[43]孙鹏远,孙建国,卢秀丽.P-SV波AVO分析[J].石油地球物理勘探,2003(02)：131-135
[44]孙鹏远,孙建国,卢秀丽.P-SV波AVO方法研究进展[J].地球物理学进展,2003(04)：602-607.
[45]孙鹏远.多属性AVO分析及弹性参数反演方法研究[D]. 吉林大学,2004.
[46]孙鹏远,孙建国,卢秀丽.P-SV波在提取SS波零偏移距反射率中的应用[J].石油物探,2004(02)：103-105
[47]孙鹏远,孙建国,卢秀丽.P-SV转换波AVO加权叠加反演[J].石油地球物理勘探,2004(03)：271-274
[48]孙鹏远,孙建国,卢秀丽.P-SV波和SS波反射系数之间的近似关系[J].地球学报,2004(06)：677-682.
[49]孙鹏远.AVO技术新进展[J].勘探地球物理进展,2005(06)：432-438.
[50]孙鹏远,孙建国,李彦鹏.P-SV波AVO响应特征分析[J].石油地球物理勘探,2005(04)：417-422.
[51]孙鹏远,孙建国,卢秀丽.P-SV波反射系数近似及其AVO属性特征[J].地球学报,2006(01)：85-89.
[52]陈天胜,刘洋,魏修成,几种P-SV转换波反射系数近似公式的比较[J].石油物探,2005,44(1)：1-3.
[53]唐旭东,甘利灯,李凌高,转换波反射系数近似公式及精度分析[J].石油物探,2008,47(2)：150-152.
[54]王者江.基于BISQ机制的三维双相正交介质正演模拟及传播特性研究[D].长春：吉林大学,2008.
[55]王恩利.基于各向异性的复合介质弹性波场模拟与特征分析[D].长春：吉林大学,2008.
[56]刘前坤.方位各向异性介质AVO及弹性波阻抗研究[D].长春：吉林大学,2008.
[57]张立勤,彭苏萍,李国发,王璞.方位AVO技术检测储层各向异性的方法和实践[J].天然气工业.2005,(10),38-40.
[58]Alford, R. M.. Shear data in the presence of azimuthal anisotropy[C]:56th Annual International Meeting. SEG, Expanded Abstracts.1986,476-479.
[59]Thomsen. L., Weak elastic anisotropy[J].Geophysics,1986.51,1954-1966.
[60]Banik N.C. An effective anisotropy parameter in transversely isotropic media[J]. Geophysics,1987,52:1654-1664.
[61]Thomsen, L., Elastic anisotropy due to aligned cracks in porous rock:Geophys. Prosp.,1995.43.805-829.
[62]Ruger. A.. Reflection coefficients and azimuthal AVO analysis in anisotropic media[D]. Ph.D. thesis. Colorado School of Mines.1996.
[63]Ruger, A.. P-wave reflection coefficients for transversely isotropic models with vertical and horizontal axis of symmetry[J]. Geophysics,1997,62:713-722;
[64]Ruger, A., Using AVO for fracture detection:analytic basis and practical solutions[J]:The Leading Edge.1997,16(10):1429-1438.
[65]Ruger, A., Analytic insight into shear-wave AVO for fractured reservoirs[C]: Proceedings of 7IWSA, special SEG volume on seismic anisotropy, Internat. Workshop on Seismic Anisotropy,1998.
[66]Ruger, A., Variation of P-wave reflectivity with offset and azimuth in anisotropic media[J]:Geophysics,1998,63,935-947.
[67]Ruger, A., Reflection coefficients and azimuthal AVO analysis in anisotropic media[C]:SEG,2001.
[68]Vavrycuk,V, and I. Psencik, PP-wave reflection coefficients in weakly elastic media[J]:Geophysics,1998,63,2129-2141
[69]Vavry C uk, V.,Weak-contrast reflection/transmission coefficients in weakly anisotropic elastic media. P-wave incidence[J]:Geophysical Journal International, 1999,138.553-562
[70]Petr Jilek. Modeling and inversion of converted-wave reflection coefficients in anisotropic media:A tool for quantitative AVO analysis,Ph.D. thesis,Colorado School of Mines,2002.
[71]Petr Jilek.Converted PS-wave Reflection Coefficients in Weakly Anisotropic Media[J]. Pure and Applied Geophysics.2002,159:1527-1562.
[72]Tsvankin,I.Seismic signatures and analysis of reflection data in anisotropic media[M]. Pergamon,2001.
[73]Jyoti Behura and Ilya Tsvankin. Small-angle AVO response of PS-waves in tilted TImedia. Geophysics,2006;,71(5):C69-C79
[74]Maxim V.Cherepanov and Tatyana N. Nefedkina,2004, Analytic description of PS wave reflection in weakly anisotropic media:SEG Technical Program Expanded Abstracts 23(1):191-194
[75]Liu Qiankun,Han Liguo.Wang ENli.and Shan Gangyi,2008,Reflection coefficients of P-SV waves in weak anisotropic media:Applied Geophysics, 5(1):18-23
[76]刘前坤,韩立国,廉玉广等,基于等效Thomsen参数的P-SV波AVO属性研究[J].西北地震学报,2008,30(3)：214-220.
[77]Cui Jie.Han Liguo, et al. P-SV wave elastic impedance and fluid identification factor in weakly anisotropic media[J].Applied Geophysics,2010,Vol.7,No.2, 135-142.
[78]甘利灯,赵邦六,杜文辉,等.弹性阻抗在岩性与流体预测中的潜力分析[J].石油物探,2005.44(5)：504-506.。
[79]曹谊,娄晓东.弹性波阻抗反演[J].油气地球物理,2003,1(3)：29-33.
[80]Whitcombe D N. Elastic impedance normalization[J]. Geophysics,2002,67(1): 60-62.
[81]Whitcombe D N.,et al. Extended elastic impedance for fluid and lithology prediction[J]. Geophysics,2002,67(1):63-67.
[82]Bruce Verwest.,et al. Elastic impedance inversion[J]. Expanded Abstracts of 70th SEG Mtg.2000:1580-1582.
[83]Santos L T., Tygel M., Ramos A C B. Reflection Impedance[J]. Expanded Abstracts of 70th SEG Mtg,2002.
[84]Santos L T., Tygel M. Impedance-type approximations of the P-P elastic reflection coefficient:Modeling and AVO inversion[J].Geophysics,2004. 69(2):592-598.
[85]Gray.D., Elastic inversion for Lame parameters [J].70th Ann. Internat. Mtg. Soc of Expl Geophys. Expanded Abstracts.2002:213-216.
[86]王保丽,印兴耀,张繁昌.基于Gray近似的弹性波阻抗方程及反演[J].石油地球物理勘探,2007.42(4)：435-439.
[87]Wang B L.,Yin X Y.,Lame parameters inversion based on elastic impedance and its application [J]. Applied Geophysics.2006.3 (3):174-178.
[88]王保丽,印兴耀,张繁昌等.基于Fatti近似的弹性波阻抗方程及反演[J].地球物理学进展,2008.23(1)：192-197.
[89]马劲风.地震勘探中广义弹性阻抗的正反演[J].地球物理学报,2003.46(3)：118-124.
[90]Ma J F., Morozov I B. A fluid detection study from Zoeppritz elastic impedance [A].76th SEG Annual Meeting Expanded Abstracts.. AVO, New Orleans. Louisiana. USA.2006:284-288.
[91]Ma J F., Morozov I B. Ray-path elastic impedance[J]. CSEG National Convention. May 2004:10-12.
[92]Kenneth Duffaut.Martin landro.et al. Shear wave elastic impedance [J].The leading Edge.2000:1223-1229
[93]Gonzalez E. F,Tapan M,et al. Near and far PS elastic impedace for discriminating fizz water from commercial gas[J]. The leading Edge,2003, 22(10):1012-1015
[94]Gonzalez E. F,Tapan M,et al. Far offset P-to-S elastic impedance for lithology and partial gas saturation(fizz water)identification:Application with well logs. Soc.Expl.Geophys.73rd Ann. Int. Mtg.Expanded Abstracts.2003,183-186.
[95]Gonzalez E. F. Physical and quantitative interpretation of seismic attributes for rocks and fluids identification[D]. Ph.D. thesis,Stanford university.2006.
[96]Jinfeng Ma,Igor B.Morozov. The exact elastic impedance for P-SV wave[J]. SEG Expanded Abstracts,2007:288-292.
[97]Igor B.Exact elastic P/SV impedance [J]. Geophysics. Vol.75.No.2:C7-C13.
[98]Martins.J.L.,An approach for elastic impedance in weakly anisotropic media[J].Expanded Abstracts of 72th Annual Meeting.2002:185-188.
[99]Martins,J.L., Elastic impedance in weakly anisotropic media[J],Geophysics,2006, Vol.71,No.3:D73-D83.
[100]苑书金,于长青.各向异性介质中的弹性阻抗及其反演[J].地球物理学进展,2006,21(2)：520-523.
[101]陈天胜,魏修成,刘洋.一种新的各向异性弹性阻抗近似公式[J].石油物探,2006,45(6)：563-569.
[102]李爱山.基于EI的岩石物性参数提取方法研究与应用[D].东营：石油大学(华东),2007.
[103]李爱山,印兴耀,张繁昌等.VTI介质中的弹性阻抗与参数提取[J].地球物理学进展,2008,23(6)：1878-1885.
[104]Castagna J P.,and Swan H. W., Principles of AVO crossplotting[J]. The Leading Edge,1998.16,337-342.
[105]Castagna J P.,Swan H. W.,and Foster D.J. Framework for AVO gradient and intercept interpretation [J]. Geophysics,1998,63,948-956.
[106]王栋,贺振华,黄德济.含气含水砂岩的高灵敏度流体识别因子的研究[J].长江大学学报(自然科学版),2008,5(3)：45-47.
[107]Smith G C.,and Gidlow P M. Weighted stacking for rock property estimation and detection of gas[J]. Geophysics prospecting,1987,35:993-1014.
[108]Castagna J P.,Batzle M L.,and Eastwood R L. Relationship between compressional and shear-wave velocities in clastic silicate rocks[J]. Geophysics,1985,50:571-581.
[109]Fatti J L., Vail P J., Smith G C.et al. Detection of gas in sandstone reservoirs using AVO analysis:A 3-D seismic case history using the geostack technique[J].Geophysics,1994,59:1362-1376.
[110]Castagna John P, Steven W Smith. Comparison of AVO indicators:A modeling study[J].Geophysics,1994,59(12):1849-1855.
[111]Goodway W,Chen T,Downton J. Improved AVO fluid detection and lithology discrimination using Lame petrophysical parameters from P and S inversion[J]. Expanded Abstracts of 67th Annual Internat SEG Mtg,1997,183-186.
[112]Goodway W. AVO and Lame constants for rock parameterization and fluid detection[J]. Recorder,26:39-60.
[113]Gray D., Goodway W.,and Chen T.,Bridging the gap:Using AVO to detect changes infundamental elastic constants[J]. SEG meeting abstracts,1999:852-855.
[114]Hedlin K.,Pore space modulus and extraction using AVO[J].70th Ann. Internat. Mtg., Soc. Expl.Geophys..Expanded Abstracts.2000:170-173.
[115]Batzle M.L..Hafmann R.. Han D H.,et al. Fluids and frequency dependent seismic velocities of rocks[J].The Leading Edge,2001(2):168-171.
[116]Batzle M.L.. Han D H..Hafmann R.Optimal hydrocarbon indicators[J]. SEG Expanded Abstracts.2001.vol20:p.1697.
[117]Russell B H. Hedlin K.Hilterman F J.et al. Fluid property discrimination with AVO:a Biot-Gassmann perspective[J]. Geophysics,2003,68(1):29-39
[118]Dillon L.Schwedersky G.Guilherme V,et al.A multiscale DHI elastic attributes evaluation[J]. The Leading Edge.2003.22(10):1024-1029.
[119]宁忠华,贺振华,黄德济.基于地震资料的高灵敏度流体识别因子[J].石油物探,2006,45(3)：239-241
[120]尹川,顾汉明.流体指示因子的敏感性分析[J].工程地球物理学报,2008,5(1)：85-88.
[121]李景叶,陈小宏.基于地震资料的储层流体识别[J].石油学报,2008,29(2)：235-238.
[122]王栋,贺振华,黄德济.新流体识别因子的构建与应用分析[J].石油物探,2009,48(2)：141-145.
[123]王栋.AVO分析与流体识别[D].成都：成都理工大学,2009.
[124]Pei F G., Zou C C, He T. Fluid sensitivity study of elastic parameters in low-medium porosity and permeability reservoir rocks[J]. Applied Geophysics,2010,7(1):1-9.
[125]Ostrander W.J. Plane waves reflection coefficient for gas sands at normal angles of incidence[J].Geophysics.1984.49(10):1637-1648.
[126]殷八斤等.AVO技术的理论与实践[M].北京：石油工业出版社,2002.
[127]Hilterman,F J.Seismic amplitude interpretation.Seismic Amplitude Interpretation[M]:Soc.Of Expl.Geophys,2001.
[128]安艺敬一,理查德.定量地震学：理论与方法[J].北京：地震出版社,1980.
[129]李邗,何诚.AVO处理解释和岩性识别[J].物探化探计算技术,2003,25(2)：145-150.
[130]何樵登.地震勘探原理和方法.北京：地质出版社,1986
[131]陆基孟.地震勘探原理.北京：石油工业出版社,1993.
[132]Knott C.G. On the reflection and refraction of elastic wave with seismological application. Phil Mag 5th ser,1899,48:64-97.
[133]Zoeppritz, K.,Erdbebenwellen VIII B, On the reflection and penetration of seismic waves through unstable layers:Goettinger Nachr.,pages 66-84.
[134]程冰洁,张玉芬,AVO简化方程的物理意义及其在油气识别中的应用[J]：物探化探计算技术,2003,25(I),26-30.
[135]王立明,提高多波速度分析精度的研究[D],西安：长安大学,2005.
[136]傅旦丹,朱宏彰,何汉漪,一种确定转换点的算法及其应用,石油物探[J],2002,41(1)：26-30.
[137]R.H.Tatham,M.D.Mccormack,Multicomponent Seismology in Petroleum Exploration[J]. Geophysics,2002,Vol62,NO.4:PP115-119;
[138]何兵寿,张会星,共转换点速度分析[J],石油地球物理勘探,2005,40(5)：510-514：
[139]毕丽飞,李振春等,转换波处理中共转换点的定位[J],勘探地球物理进展,2005,28(5)：353-356；
[140]崔杰.转换波高分辨率速度分析方法研究[D].长春：吉林大学,2008.
[141]郭向宇,凌云,魏修成,P-SV转换波共转换点的几种计算方法及实际应用[J],2002,41(2),141-148；
[142]苑春方,彭苏萍等,水平界面上P-SV转换波转换点的精确解[J],地球物理学报,2005,48(5)：1179-1184；
[143]Rutherford S. R, Williams R. H.,.Amplitude-versus-offset variations in gas sands[J]. Geophysics,1989,54(6):680-688.
[144]王云专,朱伦,王开燕等.AVO截距与梯度交会图分析及应用[J].大庆石油学院学报,2007,31(1)：12-14.
[145]王卫红,姜在兴,潘仁芳.AVO交会分析及其应用[J].西安石油学院学报：自然科学版,2003.18(2)：5-8.
[146]张显文.储层地震预测技术与储层地球物理参数反演[D].长春：吉林大学,2010.
[147]何樵登,韩立国,朱建伟等.地震波理论.长春：吉林大学出版社,2004.
[148]Crampin S. Suggestios for a consistent terminology for seismic anisotropy [J]. Geophysical Prospecting.1989,37(7):753-770.
[149]Hudson J A. Overall properties of a cracked solid[J]. Math Proc Camb phil Soc,1980,88:371-384.
[150]王德利,何樵登.裂隙型单斜介质中弹性系数的计算及波的传播特性研究[J].吉林大学学报(地球科学版),2002,32(2)：91-96.
[151]Thomsen, L. Reflection seismology over azimuthally anisotropic media: Geophysics,1988,53.304-313.
[152]徐善辉.TTI介质PS波小角度AVO分析[D].长春：吉林大学.2009.
[153]Mukerji, T., A. Jorstad, G. Mavko, and J. R. Granli, Near and far offset impedances:Seismic attributes for identifying lithofacies and pore fluids[J] Geophysical Research Letters,1998.25.4557-560.
[154]Martins. J. L., A second-order approach for P-wave elastic impedance technology[J].Studia Geophysica et Geodaetica,2003.47,545-564.
[155]Wang, Y.,Seismic amplitude inversion in reflection tomography[M]. Pergamon Press,2003.
[156]Rowbotham. P., D. Marion. R. Eden, P. Williamson. P. Lamy. and P.Swaby, The implications of anisotropy for seismic impedance inversion[J].First Break. 2003,21.53-57.
[157]Biot, M. A., Theory of propagation of elastic waves in a fluid saturated porous solid. I. Low frequency range[J]. J. Acoust. Soc. Am.,1956a.28:168-178;
[158]Biot, M. A., Theory of propagation of elastic waves in a fluid saturated porous solid. Ⅱ. Higher frequency range[J]. J. Acoust. Soc. Am.,1956b,28:179-191:
[159]Gassmann, F. Uber die Elastizitat poroser Medien. Vieteljahrsschrift der Naturforschenden Gesellchaft in Zurich[J].1951,96:1-23.
[160]王炳章,岩石物理学基本准则[J].石油物探译丛,2001,4：1-20；
[161]郭智奇.粘弹各向异性介质波场模拟与储层信息研究[D].长春：吉林大学,2008.
[2]Crampin S., Anisotropy in exploration seismics[J].,First Break,1984b, Vol.23, 19-21.
[3]Crampin,S. Evaluation of anisotropy by shear-wave splitting[J]. Geophysics, 1985.50:142-152.
[4]Crampin S. Suggestios for a consistent terminology for seismic anisotropy [J]. Geophysical Prospecting.1989,37(7):753-770.
[5]Stewart,R.R.,&Gaiser J.E. Application and interpretation of converted waves[J]. SEG Continuing Education Short Course.
[6]黄中玉.多分量地震勘探的机遇和挑战[J].石油物探,2001,40(2)：131-137.
[7]石玉梅,姚逢昌,曹宏.多波多分量天然气勘探技术的进展[J].勘探地球物理进展,2003,26(3)：172-177.
[8]陈天胜,刘洋,魏修成.纵波和转换波联合AVO反演方法研究[J].中国石油人学学报(自然科学版),2006,30(1)：33-37.
[9]Feng Hong. Hydrocarbon indicators derived from AVO analysis[D]. Ph.D. thesis,University of Calgary,2009.
[10]Ramos A.C.B, Castagna J. P, Useful approximations for converted-wave AVO[J]. Geophysics,2001,66(6):1721-1734.
[11]Hudson J A. Overall properties of a cracked solid [J]. Math Proc Camb phil Soc.1980,88:371-384.
[12]Xu X, Tsvankin I. Azimuthal AVO analysis with anisotropic spreading correction: A synthetic study [J]. The Leading Edge,2006,25(11),1336-1342.
[13]Xu X, Tsvankin I. A case study of azimuthal AVO analysis with anisotropic spreading correction[J]. The Leading Edge,2007,26(12),1552-1561.
[14]Connoly P. Elastic impedance[J]. The Leading Eage,1999,18(4):438-452.
[15]曹孟起,王九栓,邵林海.叠前弹性波阻抗反演技术及应用[J].石油地球物理勘探,2006,4(3)：323-326.
[16]王保丽,印兴耀,张繁昌.弹性阻抗反演及应用研究[J].地球物理学进展,2005,20(1)：89-92.
[17]倪逸.弹性阻抗计算的一种新方法[J].石油地球物理勘探,2003,38(2)：147-155.
[18]王保丽.叠前弹性阻抗反演理论与应用[D].东营：石油大学(华东),2005.
[19]邹文.基于地震资料的流体识别技术[D].成都：成都理工大学,2008.
[20]Stewart R. R., Gaiser J. E., Brown R. J.,and Lawton D. C. Converted-wave seismic exploration:Applications[J].Geophysics.2003.68,40-57.
[21]Larsen J.,Margrave G. and Lu, H. X. AVO analysis by simultaneous P-P and P-S weighted stacking applied to 3-C-3-D seismic data[J].69th Ann. Internat. Mtg,. Soc. of Expl. Geophys., Expanded Abstracts.1999,721-724.
[22]Kelly M. Skidmore C.and Cotton R.. P-P and P-S angle stack inversionfJ].69th Ann. Internat. Mtg,. Soc. of Expl. Geophys., Expanded Abstracts.2000.222-223.
[23]Margrave G. F.. Stewart R. R. and Larsen J. A. Joint PP and PS seismic inversion[J].The Leading Edge,2001,20,1048-1052.
[24]Veire H..and Landro M. Jiont inversion of PP-and PS-seismic data[J].71st Ann. Internat. Mtg.. Soc. of Expl. Geophys., Expanded Abstracts.2001.861-864.
[25]Engelmark F. Using converted shear waves to image reservoirs with low-inpedance contrast[J]. The Leading Edge.2000,19:600-603.
[26]Gonzalez E. F.,Mukerji. T. and Mavko G.,Facies identification using P-to-P and P-to-S AVO attributes [J].70th Ann. Internat. Mtg, Soc of Expl Geophys. Expanded Abstracts,2000:98-101.
[27]Jin S.,Cambois G. and Vuillermoz C. Shear-wave velocity and density estimation from PS-wave AVO analysis:Application to OBS dataset from the North Sea[J].Geophysics,2000,65:1446-1454.
[28]Wu Y. Estimation of gas saturation using P-to-S converted waves s[J].70th Ann. Internat. Mtg. Soc of Expl Geophys. Expanded Abstracts,2000:158-161.
[29]Zhu F., Gibson R. L.. Atkins J.,and Yuh S. H.. Distingrishing fizz gas from commercial gas reservoirs using multicomponent seismic data[J]. The Leading Edge.19,1238-1245.
[30]Bortfeld R. Approximation to the reflection and transmission coefficients of plane longitudinal and transverse waves [J]. Geophysical Prospecting,1961,9:485-502.
[31]Aki K.,Richards P. G.. Quantitative Seismology:Theory and Methods[M]. San Francisco:Freenman and Coopration.1980,146-177.
[32]郑晓东.Zoeppritz方程的近似及应用[J].石油地球物理勘探,1991,26(2)：129-144
[33]郑晓东.AVO理论和方法的一些新进展[J].石油地球物理勘探,1992,27(3)： 305-317.
[34]周竹生.P-SV波和SH波的AVO分析[J].石油地球物理勘探,1993,28(4)：430-438
[35]杨绍国,周熙襄.Zoeppritz方程的级数表达式及近似[J].石油地球物理勘探,1994,29(4)：399-412.
[36]甘永忠.用P-SV波叠加剖面合成横波层速度剖面[J].石油地球物理勘探1995.30(1)：56-61.
[37]Donati M., Martin N W. A comparison of approximations for the converted wave reflection[R]. Crewes Research Reports. Calgary,1998,159-179.
[38]Xu Y., Bancroft J C. Joint AVO analysis of PP and PS data[R]. Crewes Research Reports 34.1997,1-44.
[39]Xu Y.,AVO developments applied to Blackfoot 3C-2D broadband line[D]. Calgary:University of Calgary,1999.
[40]Wang Y H., Approximations to Zoeppritz equations and their use in AVO analysis [J].Geophysics,1999,64(6):1920-1927.
[41]Ezequiel F G, Tapan M, Gary M. Facies classification using P-to-P and P-to-S AVO attributes[J]. Expanded Abstracts of 70th Annual Internat SEG Mtg,2000.98-101
[42]Vant A, Brown R. J, Approximation the P-S reflection coefficient for small incidence angles[J]. C rewes Research Reports 13,2001,159-179.
[43]孙鹏远,孙建国,卢秀丽.P-SV波AVO分析[J].石油地球物理勘探,2003(02)：131-135
[44]孙鹏远,孙建国,卢秀丽.P-SV波AVO方法研究进展[J].地球物理学进展,2003(04)：602-607.
[45]孙鹏远.多属性AVO分析及弹性参数反演方法研究[D]. 吉林大学,2004.
[46]孙鹏远,孙建国,卢秀丽.P-SV波在提取SS波零偏移距反射率中的应用[J].石油物探,2004(02)：103-105
[47]孙鹏远,孙建国,卢秀丽.P-SV转换波AVO加权叠加反演[J].石油地球物理勘探,2004(03)：271-274
[48]孙鹏远,孙建国,卢秀丽.P-SV波和SS波反射系数之间的近似关系[J].地球学报,2004(06)：677-682.
[49]孙鹏远.AVO技术新进展[J].勘探地球物理进展,2005(06)：432-438.
[50]孙鹏远,孙建国,李彦鹏.P-SV波AVO响应特征分析[J].石油地球物理勘探,2005(04)：417-422.
[51]孙鹏远,孙建国,卢秀丽.P-SV波反射系数近似及其AVO属性特征[J].地球学报,2006(01)：85-89.
[52]陈天胜,刘洋,魏修成,几种P-SV转换波反射系数近似公式的比较[J].石油物探,2005,44(1)：1-3.
[53]唐旭东,甘利灯,李凌高,转换波反射系数近似公式及精度分析[J].石油物探,2008,47(2)：150-152.
[54]王者江.基于BISQ机制的三维双相正交介质正演模拟及传播特性研究[D].长春：吉林大学,2008.
[55]王恩利.基于各向异性的复合介质弹性波场模拟与特征分析[D].长春：吉林大学,2008.
[56]刘前坤.方位各向异性介质AVO及弹性波阻抗研究[D].长春：吉林大学,2008.
[57]张立勤,彭苏萍,李国发,王璞.方位AVO技术检测储层各向异性的方法和实践[J].天然气工业.2005,(10),38-40.
[58]Alford, R. M.. Shear data in the presence of azimuthal anisotropy[C]:56th Annual International Meeting. SEG, Expanded Abstracts.1986,476-479.
[59]Thomsen. L., Weak elastic anisotropy[J].Geophysics,1986.51,1954-1966.
[60]Banik N.C. An effective anisotropy parameter in transversely isotropic media[J]. Geophysics,1987,52:1654-1664.
[61]Thomsen, L., Elastic anisotropy due to aligned cracks in porous rock:Geophys. Prosp.,1995.43.805-829.
[62]Ruger. A.. Reflection coefficients and azimuthal AVO analysis in anisotropic media[D]. Ph.D. thesis. Colorado School of Mines.1996.
[63]Ruger, A.. P-wave reflection coefficients for transversely isotropic models with vertical and horizontal axis of symmetry[J]. Geophysics,1997,62:713-722;
[64]Ruger, A., Using AVO for fracture detection:analytic basis and practical solutions[J]:The Leading Edge.1997,16(10):1429-1438.
[65]Ruger, A., Analytic insight into shear-wave AVO for fractured reservoirs[C]: Proceedings of 7IWSA, special SEG volume on seismic anisotropy, Internat. Workshop on Seismic Anisotropy,1998.
[66]Ruger, A., Variation of P-wave reflectivity with offset and azimuth in anisotropic media[J]:Geophysics,1998,63,935-947.
[67]Ruger, A., Reflection coefficients and azimuthal AVO analysis in anisotropic media[C]:SEG,2001.
[68]Vavrycuk,V, and I. Psencik, PP-wave reflection coefficients in weakly elastic media[J]:Geophysics,1998,63,2129-2141
[69]Vavry C uk, V.,Weak-contrast reflection/transmission coefficients in weakly anisotropic elastic media. P-wave incidence[J]:Geophysical Journal International, 1999,138.553-562
[70]Petr Jilek. Modeling and inversion of converted-wave reflection coefficients in anisotropic media:A tool for quantitative AVO analysis,Ph.D. thesis,Colorado School of Mines,2002.
[71]Petr Jilek.Converted PS-wave Reflection Coefficients in Weakly Anisotropic Media[J]. Pure and Applied Geophysics.2002,159:1527-1562.
[72]Tsvankin,I.Seismic signatures and analysis of reflection data in anisotropic media[M]. Pergamon,2001.
[73]Jyoti Behura and Ilya Tsvankin. Small-angle AVO response of PS-waves in tilted TImedia. Geophysics,2006;,71(5):C69-C79
[74]Maxim V.Cherepanov and Tatyana N. Nefedkina,2004, Analytic description of PS wave reflection in weakly anisotropic media:SEG Technical Program Expanded Abstracts 23(1):191-194
[75]Liu Qiankun,Han Liguo.Wang ENli.and Shan Gangyi,2008,Reflection coefficients of P-SV waves in weak anisotropic media:Applied Geophysics, 5(1):18-23
[76]刘前坤,韩立国,廉玉广等,基于等效Thomsen参数的P-SV波AVO属性研究[J].西北地震学报,2008,30(3)：214-220.
[77]Cui Jie.Han Liguo, et al. P-SV wave elastic impedance and fluid identification factor in weakly anisotropic media[J].Applied Geophysics,2010,Vol.7,No.2, 135-142.
[78]甘利灯,赵邦六,杜文辉,等.弹性阻抗在岩性与流体预测中的潜力分析[J].石油物探,2005.44(5)：504-506.。
[79]曹谊,娄晓东.弹性波阻抗反演[J].油气地球物理,2003,1(3)：29-33.
[80]Whitcombe D N. Elastic impedance normalization[J]. Geophysics,2002,67(1): 60-62.
[81]Whitcombe D N.,et al. Extended elastic impedance for fluid and lithology prediction[J]. Geophysics,2002,67(1):63-67.
[82]Bruce Verwest.,et al. Elastic impedance inversion[J]. Expanded Abstracts of 70th SEG Mtg.2000:1580-1582.
[83]Santos L T., Tygel M., Ramos A C B. Reflection Impedance[J]. Expanded Abstracts of 70th SEG Mtg,2002.
[84]Santos L T., Tygel M. Impedance-type approximations of the P-P elastic reflection coefficient:Modeling and AVO inversion[J].Geophysics,2004. 69(2):592-598.
[85]Gray.D., Elastic inversion for Lame parameters [J].70th Ann. Internat. Mtg. Soc of Expl Geophys. Expanded Abstracts.2002:213-216.
[86]王保丽,印兴耀,张繁昌.基于Gray近似的弹性波阻抗方程及反演[J].石油地球物理勘探,2007.42(4)：435-439.
[87]Wang B L.,Yin X Y.,Lame parameters inversion based on elastic impedance and its application [J]. Applied Geophysics.2006.3 (3):174-178.
[88]王保丽,印兴耀,张繁昌等.基于Fatti近似的弹性波阻抗方程及反演[J].地球物理学进展,2008.23(1)：192-197.
[89]马劲风.地震勘探中广义弹性阻抗的正反演[J].地球物理学报,2003.46(3)：118-124.
[90]Ma J F., Morozov I B. A fluid detection study from Zoeppritz elastic impedance [A].76th SEG Annual Meeting Expanded Abstracts.. AVO, New Orleans. Louisiana. USA.2006:284-288.
[91]Ma J F., Morozov I B. Ray-path elastic impedance[J]. CSEG National Convention. May 2004:10-12.
[92]Kenneth Duffaut.Martin landro.et al. Shear wave elastic impedance [J].The leading Edge.2000:1223-1229
[93]Gonzalez E. F,Tapan M,et al. Near and far PS elastic impedace for discriminating fizz water from commercial gas[J]. The leading Edge,2003, 22(10):1012-1015
[94]Gonzalez E. F,Tapan M,et al. Far offset P-to-S elastic impedance for lithology and partial gas saturation(fizz water)identification:Application with well logs. Soc.Expl.Geophys.73rd Ann. Int. Mtg.Expanded Abstracts.2003,183-186.
[95]Gonzalez E. F. Physical and quantitative interpretation of seismic attributes for rocks and fluids identification[D]. Ph.D. thesis,Stanford university.2006.
[96]Jinfeng Ma,Igor B.Morozov. The exact elastic impedance for P-SV wave[J]. SEG Expanded Abstracts,2007:288-292.
[97]Igor B.Exact elastic P/SV impedance [J]. Geophysics. Vol.75.No.2:C7-C13.
[98]Martins.J.L.,An approach for elastic impedance in weakly anisotropic media[J].Expanded Abstracts of 72th Annual Meeting.2002:185-188.
[99]Martins,J.L., Elastic impedance in weakly anisotropic media[J],Geophysics,2006, Vol.71,No.3:D73-D83.
[100]苑书金,于长青.各向异性介质中的弹性阻抗及其反演[J].地球物理学进展,2006,21(2)：520-523.
[101]陈天胜,魏修成,刘洋.一种新的各向异性弹性阻抗近似公式[J].石油物探,2006,45(6)：563-569.
[102]李爱山.基于EI的岩石物性参数提取方法研究与应用[D].东营：石油大学(华东),2007.
[103]李爱山,印兴耀,张繁昌等.VTI介质中的弹性阻抗与参数提取[J].地球物理学进展,2008,23(6)：1878-1885.
[104]Castagna J P.,and Swan H. W., Principles of AVO crossplotting[J]. The Leading Edge,1998.16,337-342.
[105]Castagna J P.,Swan H. W.,and Foster D.J. Framework for AVO gradient and intercept interpretation [J]. Geophysics,1998,63,948-956.
[106]王栋,贺振华,黄德济.含气含水砂岩的高灵敏度流体识别因子的研究[J].长江大学学报(自然科学版),2008,5(3)：45-47.
[107]Smith G C.,and Gidlow P M. Weighted stacking for rock property estimation and detection of gas[J]. Geophysics prospecting,1987,35:993-1014.
[108]Castagna J P.,Batzle M L.,and Eastwood R L. Relationship between compressional and shear-wave velocities in clastic silicate rocks[J]. Geophysics,1985,50:571-581.
[109]Fatti J L., Vail P J., Smith G C.et al. Detection of gas in sandstone reservoirs using AVO analysis:A 3-D seismic case history using the geostack technique[J].Geophysics,1994,59:1362-1376.
[110]Castagna John P, Steven W Smith. Comparison of AVO indicators:A modeling study[J].Geophysics,1994,59(12):1849-1855.
[111]Goodway W,Chen T,Downton J. Improved AVO fluid detection and lithology discrimination using Lame petrophysical parameters from P and S inversion[J]. Expanded Abstracts of 67th Annual Internat SEG Mtg,1997,183-186.
[112]Goodway W. AVO and Lame constants for rock parameterization and fluid detection[J]. Recorder,26:39-60.
[113]Gray D., Goodway W.,and Chen T.,Bridging the gap:Using AVO to detect changes infundamental elastic constants[J]. SEG meeting abstracts,1999:852-855.
[114]Hedlin K.,Pore space modulus and extraction using AVO[J].70th Ann. Internat. Mtg., Soc. Expl.Geophys..Expanded Abstracts.2000:170-173.
[115]Batzle M.L..Hafmann R.. Han D H.,et al. Fluids and frequency dependent seismic velocities of rocks[J].The Leading Edge,2001(2):168-171.
[116]Batzle M.L.. Han D H..Hafmann R.Optimal hydrocarbon indicators[J]. SEG Expanded Abstracts.2001.vol20:p.1697.
[117]Russell B H. Hedlin K.Hilterman F J.et al. Fluid property discrimination with AVO:a Biot-Gassmann perspective[J]. Geophysics,2003,68(1):29-39
[118]Dillon L.Schwedersky G.Guilherme V,et al.A multiscale DHI elastic attributes evaluation[J]. The Leading Edge.2003.22(10):1024-1029.
[119]宁忠华,贺振华,黄德济.基于地震资料的高灵敏度流体识别因子[J].石油物探,2006,45(3)：239-241
[120]尹川,顾汉明.流体指示因子的敏感性分析[J].工程地球物理学报,2008,5(1)：85-88.
[121]李景叶,陈小宏.基于地震资料的储层流体识别[J].石油学报,2008,29(2)：235-238.
[122]王栋,贺振华,黄德济.新流体识别因子的构建与应用分析[J].石油物探,2009,48(2)：141-145.
[123]王栋.AVO分析与流体识别[D].成都：成都理工大学,2009.
[124]Pei F G., Zou C C, He T. Fluid sensitivity study of elastic parameters in low-medium porosity and permeability reservoir rocks[J]. Applied Geophysics,2010,7(1):1-9.
[125]Ostrander W.J. Plane waves reflection coefficient for gas sands at normal angles of incidence[J].Geophysics.1984.49(10):1637-1648.
[126]殷八斤等.AVO技术的理论与实践[M].北京：石油工业出版社,2002.
[127]Hilterman,F J.Seismic amplitude interpretation.Seismic Amplitude Interpretation[M]:Soc.Of Expl.Geophys,2001.
[128]安艺敬一,理查德.定量地震学：理论与方法[J].北京：地震出版社,1980.
[129]李邗,何诚.AVO处理解释和岩性识别[J].物探化探计算技术,2003,25(2)：145-150.
[130]何樵登.地震勘探原理和方法.北京：地质出版社,1986
[131]陆基孟.地震勘探原理.北京：石油工业出版社,1993.
[132]Knott C.G. On the reflection and refraction of elastic wave with seismological application. Phil Mag 5th ser,1899,48:64-97.
[133]Zoeppritz, K.,Erdbebenwellen VIII B, On the reflection and penetration of seismic waves through unstable layers:Goettinger Nachr.,pages 66-84.
[134]程冰洁,张玉芬,AVO简化方程的物理意义及其在油气识别中的应用[J]：物探化探计算技术,2003,25(I),26-30.
[135]王立明,提高多波速度分析精度的研究[D],西安：长安大学,2005.
[136]傅旦丹,朱宏彰,何汉漪,一种确定转换点的算法及其应用,石油物探[J],2002,41(1)：26-30.
[137]R.H.Tatham,M.D.Mccormack,Multicomponent Seismology in Petroleum Exploration[J]. Geophysics,2002,Vol62,NO.4:PP115-119;
[138]何兵寿,张会星,共转换点速度分析[J],石油地球物理勘探,2005,40(5)：510-514：
[139]毕丽飞,李振春等,转换波处理中共转换点的定位[J],勘探地球物理进展,2005,28(5)：353-356；
[140]崔杰.转换波高分辨率速度分析方法研究[D].长春：吉林大学,2008.
[141]郭向宇,凌云,魏修成,P-SV转换波共转换点的几种计算方法及实际应用[J],2002,41(2),141-148；
[142]苑春方,彭苏萍等,水平界面上P-SV转换波转换点的精确解[J],地球物理学报,2005,48(5)：1179-1184；
[143]Rutherford S. R, Williams R. H.,.Amplitude-versus-offset variations in gas sands[J]. Geophysics,1989,54(6):680-688.
[144]王云专,朱伦,王开燕等.AVO截距与梯度交会图分析及应用[J].大庆石油学院学报,2007,31(1)：12-14.
[145]王卫红,姜在兴,潘仁芳.AVO交会分析及其应用[J].西安石油学院学报：自然科学版,2003.18(2)：5-8.
[146]张显文.储层地震预测技术与储层地球物理参数反演[D].长春：吉林大学,2010.
[147]何樵登,韩立国,朱建伟等.地震波理论.长春：吉林大学出版社,2004.
[148]Crampin S. Suggestios for a consistent terminology for seismic anisotropy [J]. Geophysical Prospecting.1989,37(7):753-770.
[149]Hudson J A. Overall properties of a cracked solid[J]. Math Proc Camb phil Soc,1980,88:371-384.
[150]王德利,何樵登.裂隙型单斜介质中弹性系数的计算及波的传播特性研究[J].吉林大学学报(地球科学版),2002,32(2)：91-96.
[151]Thomsen, L. Reflection seismology over azimuthally anisotropic media: Geophysics,1988,53.304-313.
[152]徐善辉.TTI介质PS波小角度AVO分析[D].长春：吉林大学.2009.
[153]Mukerji, T., A. Jorstad, G. Mavko, and J. R. Granli, Near and far offset impedances:Seismic attributes for identifying lithofacies and pore fluids[J] Geophysical Research Letters,1998.25.4557-560.
[154]Martins. J. L., A second-order approach for P-wave elastic impedance technology[J].Studia Geophysica et Geodaetica,2003.47,545-564.
[155]Wang, Y.,Seismic amplitude inversion in reflection tomography[M]. Pergamon Press,2003.
[156]Rowbotham. P., D. Marion. R. Eden, P. Williamson. P. Lamy. and P.Swaby, The implications of anisotropy for seismic impedance inversion[J].First Break. 2003,21.53-57.
[157]Biot, M. A., Theory of propagation of elastic waves in a fluid saturated porous solid. I. Low frequency range[J]. J. Acoust. Soc. Am.,1956a.28:168-178;
[158]Biot, M. A., Theory of propagation of elastic waves in a fluid saturated porous solid. Ⅱ. Higher frequency range[J]. J. Acoust. Soc. Am.,1956b,28:179-191:
[159]Gassmann, F. Uber die Elastizitat poroser Medien. Vieteljahrsschrift der Naturforschenden Gesellchaft in Zurich[J].1951,96:1-23.
[160]王炳章,岩石物理学基本准则[J].石油物探译丛,2001,4：1-20；
[161]郭智奇.粘弹各向异性介质波场模拟与储层信息研究[D].长春：吉林大学,2008.
© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号 地址：北京市海淀区学院路29号 邮编：100083 电话：办公室：(+86 10)66554848；文献借阅、咨询服务、科技查新：66554700 |