一种裂缝型储层流体的识别方法及装置制造方法
【专利摘要】本申请公开了一种裂缝型储层流体的识别方法及装置，该识别方法包括以下步骤：将三维地震纵波资料进行处理，生成多方位的道集数据；利用所生成的道集数据进行反演得到储层的各向同性梯度Biso和各向异性梯度Bani，并计算各向异性梯度和各向同性梯度的比值利用裂缝横向宽度与纵向高度的比值d、储层流体的第二拉梅常数μ'以及背景介质拉梅常数λ,μ对进行修正，获得流体因子Ifluid；利用该流体因子Ifluid对储层流体进行识别。识别装置相应的包括：道集数据生成单元、第一级流体因子计算单元、第二级流体因子获取单元、流体因子获取单元以及流体识别单元。通过本申请中的识别方法及装置能减小表征储层流体异常特性的流体因子的相对误差，准确识别裂缝型储层流体的类型。
【专利说明】一种裂缝型储层流体的识别方法及装置
【技术领域】
[0001] 本申请涉及石油勘探领域，尤其涉及一种裂缝型储层流体的识别方法及装置。
【背景技术】
[0002] 裂缝型油气藏的产量占全世界石油天然气总产量的一半以上，在国内的油气勘探 中，裂缝型油气藏是一个重要的勘探领域。碳酸盐岩、泥岩、砂岩、火山岩和变质岩中均有裂 缝发育。裂缝既是油、气、水等地下流体的运移通道，也是流体的储集空间，裂缝能极大地提 升储层的渗透率。因此，识别裂缝中所贮存的流体，对于寻找油气储层具有重大意义。
[0003] 利用地震资料研究岩石裂缝所贮存流体的类型和特性，是根据岩石物理等效介质 理论，将与储层流体有关的异常特性表征为流体因子，再依托流体因子实现储层流体的类 型判识。流体识别的质量主要取决于两方面，一是反演计算的弹性参数是否可靠；二是构 建的流体因子对流体类型是否敏感且稳定。储层内发育裂缝时对外表现为各向异性特征， 所以用各向异性参数来构建流体因子。所谓各向异性是指地震波的物理性质（振幅、速度、 频率等）随观测方位的变化而变化。对于垂直排列的平行裂缝，可以将其等效为HTI介质 (Horizontal Transverse Isotropy,即具有水平对称轴的横向各向同性介质）模型来研 究。若入射角为9,测线方位角为（KRuger认为HTI介质分界面处的纵波反射系数是各 向异性参数、入射角和方位角的函数：
[0004]
【权利要求】
1. 一种裂缝型储层流体的识别方法，其特征在于，包括以下步骤： S1，将获取的整个研究区域的三维地震纵波资料进行各向异性偏移处理，生成道集数 据，所述道集数据包括三维地震纵波射线的方位角、入射角和反射振幅的信息； 52, 利用所生成的道集数据中反射振幅R与入射角Θ和方位角φ之间的关系，并考虑 大于30度的入射角Θ对反射振幅R的影响，进行反演得到储层的各向同性梯度B is°和各 向异性梯度Bani，并计算各向异性梯度Bani和各向同性梯度B is°的比值，将所述比值记为第 一级流体因子/JL ; 53, 利用裂缝横向宽度与纵向高度的比值d对所述第一级流体因子/^,进行修正，获 得第二级流体因子，所述裂缝横向宽度与纵向高度的比值d是通过运用背景介质裂缝 结构的各向异性反演计算所得到的； 54, 利用所获得的第二级流体因子对储层流体类型进行区分。
2. 根据权利要求1所述的裂缝型储层流体的识别方法，其特征在于，所述利用所生成 的道集数据中反射振幅R与入射角Θ和方位角φ之间的关系，并考虑入射角Θ大于30度 时对反射振幅R的影响，进行反演得到储层流体的各向同性梯度B is°和各向异性梯度Bani, 包括： 利用反射振幅R与入射角Θ和方位角φ之间的下述关系式进行AVAZ反演，得到储层 的各向同性梯度Bis°、各向异性梯度Bani以及其它四个参数Rtl, Cis°，Canil和Cani2 : R( θ , φ) = R0+ (Biso+Banicos2 Φ ') sin2 θ + (Ciso+Canilcos4 Φ ' +Cani2Sin2 Φ ' cos2 Φ ') sin2 θ tan2 θ。
3. 根据权利要求2所述的裂缝型储层流体的识别方法，其特征在于，所述AVAZ反演包 括： 将所述反射振幅R与入射角Θ和方位角φ之间的关系式简化为： R( θ，φ) = △&+88"。+〇8-+0(^<^(^1+卩(：-2，其中六、8、(：、04和？是入射角和方位角 的组合： A=I B = sin2 Θ C = sin2 Θ cos2 Φ ' D = sin2 Θ tan2 θ E = sin2 Θ tan2 Θ cos4 Φ ' F = sin2 θ tan2 θ sin2 Φ ' cos2 Φ ' 将上述公式构建超定方程组=Mnx6X6xi = Rnxi 其中η为覆盖次数，
通过求解所述超定方程组，得到储层流体的各向同性梯度Bis°和各向异性梯度Bani。
4. 根据权利要求1-3任一项所述的裂缝型储层流体的识别方法，其特征在于，所述步 骤S2还包括：在获得所述第二级流体因子/g,后，利用储层流体的第二拉梅常数μ'以及 背景介质的拉梅常数λ，μ对第二级流体因子/iL进行修正，获得流体因子Ifluid，所述背景 介质的拉梅常数λ，μ是通过利用所述道集数据中的反射振幅信息反演计算得到的； 所述步骤S3为利用流体因子Ifluid对储层流体进行区分。
5. 根据权利要求4所述的裂缝型储层流体的识别方法，其特征在于，所述利用流体因 子Ifluid对储层流体进行区分包括： 将所获得的流体因子Ifluid与预设的阈值进行对比来判断储层流体的类型， 或， 将所获得的流体因子Ifluid与现有的流体因子谱库中的数据进行对比来判断储层流体 的类型。
6. -种裂缝型储层流体的识别装置，其特征在于，包括： 道集数据生成单元，所述道集数据生成单元用于将获取的整个研究区域的三维地震纵 波资料进行各向异性偏移处理，生成道集数据，所述道集数据包括三维地震纵波射线的方 位角、入射角和反射振幅的信息； 第一级流体因子计算单元，所述第一级流体因子计算单元用于利用所生成的道集数 据，通过对方位角和入射角进行反演得到储层流体的各向同性梯度Bis°和各向异性梯度 Bani，并计算各向异性梯度Bani和各向同性梯度Bis°的比值，将所述比值记为第一级流体因子 第二级流体因子获取单元，所述第二级流体因子获取单元用于利用裂缝横向宽度与纵 向高度的比值d对所述第一级流体因子巧^进行修正，获得第二级流体因子/^，裂缝横向 宽度与纵向高度的比值d是通过运用背景介质裂缝结构的各向异性反演计算得到的； 流体识别单元，所述流体识别单元用于利用所获得的第二级流体因子对储层流体 进行识别。
7. 根据权利要求6所述的裂缝型储层流体的识别装置，其特征在于，所述第一级流体 因子计算单元包括第一子单元，所述第一子单元用于利用反射振幅R与入射角Θ和方位角 Φ之间的下述关系式进行AVAZ反演，得到储层流体的各向同性梯度Bis°和各向异性梯度 Bani以及其它四个参数R。，Cis°，Canil和Cani2 : R( θ , φ) = R0+ (Biso+Banicos2 Φ ') sin2 θ + (Ciso+Canilcos4 Φ ' +Cani2Sin2 Φ ' cos2 Φ ') sin2 θ tan2 θ。
8. 根据权利要求7所述的裂缝型储层流体的识别装置，其特征在于，所述第一子单元 包括： 关系式简化单元，所述关系式简化单元用于将所述反射振幅R与入射角Θ和方位角φ 之间的关系式简化为： R( θ，φ) = △&+88"。+〇8-+0(^<^(^1+卩(：-2，其中六、8、(：、04和？是入射角和方位角 的组合： A=I B = sin2 θ C = sin2 Θ cos2 Φ ' D = sin2 θ tan2 θ E = sin2 θ tan2 θ cos4 Φ ' F = sin2 θ tan2 θ sin2 Φ ' cos2 Φ ' 超定方程组构建单元，所述超定方程组构建单元用于将关系式简化单元简化后的关系 式构建超定方程组：Mnx6X6x i = Rnxi 其中η为覆盖次数，
超定方程组求解单元，所述超定方程组求解单元用于求解所述超定方程组，得到储层 流体的各向同性梯度Bis°和各向异性梯度Bani。
9. 根据权利要求6-8任一项所述的裂缝型储层流体的识别装置，其特征在于，所述第 二级流体因子获取单元包括流体因子获取单元，所述流体因子获取单元用于利用储层流体 的第二拉梅常数μ'以及背景介质的拉梅常数λ，μ对第二级流体因子/^^进行修正，获 得流体因子Ifluid，所述背景介质的拉梅常数λ，μ是通过利用所述道集数据中的反射振幅 反演计算得到的； 所述流体识别单元用于利用流体因子Ifluid对储层流体进行区分。
10. 根据权利要求9所述的裂缝型储层流体的识别装置，其特征在于，所述流体识别单 元包括第二子单元或第三子单元，所述第二子单元用于将所获得的流体因子、^与预设的 阈值进行对比来判断储层流体的类型，所述第三子单元用于将所获得的流体因子I fluid与现 有的流体因子谱库中的数据进行对比来判断储层流体的类型。
【文档编号】G01V1/30GK104375182SQ201410645053
【公开日】2015年2月25日 申请日期:2014年11月10日 优先权日:2014年11月10日
【发明者】谢春辉, 雍学善, 杨午阳, 周春雷, 王洪求 申请人:中国石油天然气股份有限公司