• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 21
    • 下载费用:30 金币  

    重庆时时彩37期开: 等效孔隙裂缝介质的弹性阻抗及广义流体因子分析方法.pdf

    关 键 词:
    等效 孔隙 裂缝 介质 弹性 阻抗 广义 流体 因子分析 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201710084219.8

    申请日:

    2017.02.16

    公开号:

    CN106772608A

    公开日:

    2017.05.31

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 实质审查的生效IPC(主分类):G01V 1/30申请日:20170216|||公开
    IPC分类号: G01V1/30; G01V1/28 主分类号: G01V1/30
    申请人: 甘肃省地震局
    发明人: 窦喜英; 王恩利; 张博; 张辉; 闫国亮; 陈启艳
    地址: 730000 甘肃省兰州市城关区渭源路街道东岗西路450号
    优先权:
    专利代理机构: 北京科亿知识产权代理事务所(普通合伙) 11350 代理人: 汤东凤
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201710084219.8

    授权公告号:

    |||

    法律状态公告日:

    2017.06.23|||2017.05.31

    法律状态类型:

    实质审查的生效|||公开

    摘要

    等效孔隙裂缝介质的弹性阻抗及广义流体因子分析方法,将Thomsen等效孔隙裂缝模型与线性滑动模型等价对比,得出Thomsen裂缝模型的广义各向异性参数,引入到HTI介质的纵波、转换横波弹性阻抗公式中,最终得到对应Thomsen裂缝介质的弹性阻抗公式,建立裂缝参数、岩性参数以及物性参数与弹性阻抗的直接函数关系,对弹性阻抗归一化后进行数值模拟。模拟结果表明,裂缝型油气藏弹性阻抗对裂缝填充物性质极其敏感,且在不同流体饱和条件下对基岩的横纵波速度平方比、基质孔隙度、介质裂缝密度的响应程度差异显著。

    权利要求书

    1.等效孔隙裂缝介质的弹性阻抗及广义流体因子分析方法,其特征在于,由如下步骤
    构成:
    步骤一、以低频孔隙裂缝模型为等效孔隙裂缝模型,将其与线性滑动模型等价比较后
    建立的新柔量表达式引入以柔量为参数的广义各向异性参数中,得到新广义各向异性参数
    关系式(1)~(4)
    <mrow> <msup> <mi>&epsiv;</mi> <mrow> <mo>(</mo> <mi>v</mi> <mo>)</mo> </mrow> </msup> <mo>=</mo> <mo>-</mo> <mfrac> <mn>8</mn> <mn>3</mn> </mfrac> <mi>e</mi> <mi>q</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
    <mrow> <msup> <mi>&delta;</mi> <mrow> <mo>(</mo> <mi>v</mi> <mo>)</mo> </mrow> </msup> <mo>=</mo> <mo>-</mo> <mfrac> <mn>8</mn> <mn>3</mn> </mfrac> <mi>e</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <mn>2</mn> <mi>g</mi> </mrow> <mrow> <mn>1</mn> <mo>-</mo> <mi>g</mi> </mrow> </mfrac> <mi>q</mi> <mo>+</mo> <mfrac> <mrow> <mn>4</mn> <mi>g</mi> </mrow> <mrow> <mn>3</mn> <mo>-</mo> <mn>2</mn> <mi>g</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
    <mrow> <msup> <mi>&gamma;</mi> <mrow> <mo>(</mo> <mi>v</mi> <mo>)</mo> </mrow> </msup> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <mn>8</mn> <mi>e</mi> </mrow> <mrow> <mn>3</mn> <mrow> <mo>(</mo> <mn>3</mn> <mo>-</mo> <mn>2</mn> <mi>g</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
    <mrow> <msup> <mi>&eta;</mi> <mrow> <mo>(</mo> <mi>v</mi> <mo>)</mo> </mrow> </msup> <mo>=</mo> <mfrac> <mrow> <mn>8</mn> <mi>e</mi> </mrow> <mn>3</mn> </mfrac> <mi>g</mi> <mo>&lsqb;</mo> <mfrac> <mn>4</mn> <mrow> <mn>3</mn> <mo>-</mo> <mn>2</mn> <mi>g</mi> </mrow> </mfrac> <mo>-</mo> <mfrac> <mn>1</mn> <mrow> <mn>1</mn> <mo>-</mo> <mi>g</mi> </mrow> </mfrac> <mi>q</mi> <mo>&rsqb;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
    其中,e为裂缝密度,g为横、纵波速度平方比,q为与流体相关的系数;
    步骤二、将新广义Thomsen各向异性参数关系式(1)~(4)代入陈天胜推导的HTI介质PP
    波弹性阻抗公式中,整理得到等效孔隙裂缝介质的PP波弹性阻抗公式如下:

    公式中的q与步骤一中相同,将q中的Kf替换为气体体积模量,即Kf=0.02GPa,变为气
    饱和裂缝介质的弹性阻抗,Kf替换为水体积模量,即Kf=2GPa,即为水饱和裂缝介质的弹性
    阻抗;
    步骤三、把关系式(1)~(4)代入HTI介质PSV波弹性阻抗公式中,整理得到等效孔隙裂
    缝介质的PSV波弹性阻抗公式,

    同理代入气和水的体积模量,能够得到气饱和与水饱和裂缝介质的PSV波弹性阻抗。
    2.根据权利要求1所述的方法,其特征在于,选用Hilterman的第1类AVO气砂岩的P波速
    度、S波速度和密度作为背景参数,修改该砂岩为基质孔隙度在0~10%范围内的等效孔隙
    裂缝介质,在裂缝密度e(模型①)、基质孔隙度φp(模型②)、横纵波速度平方比g(模型③)3
    方面进行参数变化的改造,并考虑水饱和与气饱和条件的影响(见表1),开展PP波和PSV波
    弹性阻抗的数值模拟;
    数值模拟的具体流程为:
    第一步:计算每种改造对应的归一化各向同性弹性阻抗;
    第二步:计算气饱和与水饱和等效孔隙裂缝介质的归一化弹性阻抗;
    第三步:将第二步中弹性阻抗结果与相应第一步中的各向同性弹性阻抗结果做差,得
    到相应流体饱和状态下裂缝介质的阻抗差;
    第四步:提取入射角θ=30°,方位角获取相同流体饱和状态下阻抗差的方
    位各向异性情况。
    3.根据权利要求1所述的方法,其特征在于,等效孔隙裂缝介质的广义流体因子构建方
    法为:
    为实现裂缝型介质的流体识别,尝试用裂缝介质的PP波和PSV波弹性阻抗替代
    Goodwayg构建的λρ和μρ、Russell构建的ρf、宁忠华构建的σHSFIF和σHSFIF1、王栋构建的FIFW和
    F1、以及许平构建的FIFP这几种常规流体因子公式中的纵波和横波阻抗,将裂缝密度、入射
    角和方位角参数引入其中,对应形成g_λρ、g_μρ、g_ρf、g_σHSFIF、g_σHSFIF1、g_FIFW、g_F1和g_
    FIFP广义流体识别因子,具体形式如下:








    说明书

    等效孔隙裂缝介质的弹性阻抗及广义流体因子分析方法

    技术领域

    本发明涉及裂缝储层弹性阻抗研究领域,特别涉及等效孔隙裂缝介质的广义流体
    因子分析方法研究。

    背景技术

    弹性阻抗是重要的叠前储层预测工具。自1999年Connolly提出各向同性介质纵波
    弹性阻抗的概念以来,弹性阻抗理论得到迅速发展,一系列具有明显物理意义的概念被相
    继推出,如VerWest的射线弹性阻抗(VEI),Whitcombe等提出的归一化弹性阻抗、扩展性弹
    性阻抗(EEI),马劲风提出的广义弹性阻抗(GEI)、反射率阻抗(RI)、Zoeppritz弹性阻抗
    (ZEI),Snatos提出的射线弹性阻抗(REI),苑春芳对Whitcombe的扩展弹性阻抗进行改进得
    到的新弹性阻抗。在性阻抗发展的过程中,为解决P-SV转换横波的层位标定和岩性反演问
    题,Landro等、Kenneth等先后提出了横波弹性阻抗(SEI),基于此概念,Ezequiel推导了以
    反射横波角度和入射纵波角度为参数的转换横波弹性阻抗(PSEI),Ma和Morozov基于
    Zoeppritz方程推导了转换横波弹性阻抗的精确表达形式(PSZEI)。

    与此同时,由于裂缝型油气藏逐渐成为重要的剩余油气资源和研究目标,各向异
    性介质的弹性阻抗成为研究热点,Martins于2002年首先提出了各向异性弹性阻抗概念,并
    在2006年推导了弱各向异性介质的纵波弹性阻抗公式,将弹性阻抗分为各向同性与各向异
    性两部分,分析它们之间的差异,强调了各向异性在弹性阻抗反演中所起的重要作用。
    Jilek基于弱各向异性转换横波反射系数公式,推导出复杂的P-SV波、P-SH波弹性阻抗公
    式。陈天胜推导出HTI介质的纵波弹性阻抗公式。刘前坤基于Cherepanov的HTI介质PS波反
    射系数公式,引入Thomsen等效各向异性参数,推导了HTI介质中PS波弹性阻抗公式,并进行
    了含有不同流体介质的弹性阻抗讨论。崔杰等基于简化的Jilek反射系数公式,重新推导了
    新的各向异性转换波弹性阻抗公式。王恩利基于Ruger的HTI介质PP波反射系数公式,引入
    裂缝密度参数化的各向异性参数,推导出包含裂缝参数的弹性阻抗公式。阻抗分析与流体
    识别密不可分。对AVO信息以及多波多分量地震资料的深入研究推动了储层流体识别技术
    的发展?;诟飨蛲越橹世砺?,国内外许多学者提出很多流体指示因子和识别方法,如
    Goodway提出了流体识别异常的LMR法、Hedllin提出了孔隙模量法、Russell提出了流体组
    分识别的Russell法、Dillon提出了直接油气指示(DHI)的波阻抗差分法、陈遵德 提出了利
    用纵、横波速度、密度以及一些其它地震属性实现岩性油气藏储层预测的方法?;菩鞯绿岢?br />了泊松比流体识别因子。宁忠华,贺振华提出了由波阻抗量纲的0次幂和2次幂的组合而成
    的具有较高灵敏度的流体识别因子。尹川研究了流体因子的敏感性,利用不同敏感程度的
    流体因子进行交会分析对流体的识别效果。王栋,贺振华考虑到波阻抗的高次量纲能将不
    同含流体砂岩的差异放大,构建了包含有波阻抗量纲4次方的高灵敏度流体识别因子,又于
    2009年对该流体因子进行修改,增加对流体敏感的λρ项,得到新流体识别因子,取得了良好
    的效果。许平提出了类似于王栋2008年给出的流体识别因子的新流体因子。Pei联合K/μ、λ/
    μ两个参数,提出了可以有效识别低孔渗储层中流体的新流体因子λK/μ2。郑静静综合流体
    因子ρf和泊松阻抗PI的优点,构建了新的流体因子,更好地提高了抗噪性能和流体检测能
    力。

    传统的裂缝储层弹性阻抗研究多基于静态等效介质理论,由于未充分考虑岩石孔
    隙和流体,距离地下介质的真实情况依然有较大差距,理论误差较大。本文将Thomsen等效
    多孔裂缝介质模型(低频条件下)的广义各向异性参数引入HTI介质的PP波、P-SV波弹性阻
    抗公式中,建立裂缝参数、基质孔隙度和弹性阻抗之间的显性关系式,分析裂缝密度、裂缝
    填充物、基岩的纵横波速度比、基质孔隙度对弹性阻抗的影响,以避免由各向异性参数向裂
    缝参数和基质孔隙度参数转换过程中引入的误差。此外,在继承前人研究成果的同时,尝试
    将流体因子引入到等效孔隙裂缝介质中,研究裂缝参数、基质孔隙度、入射角和方位角对流
    体因子的影响。

    发明内容

    为解决上述现有技术存在的问题,本发明的目的在于提供由裂缝参数、基质孔隙
    度、裂缝填充物等参数直接表示的等效孔隙裂缝介质的弹性阻抗及广义流体因子分析方
    法。本发明中,等效裂缝介质理论作为裂缝介质的一种近似理论,更加接近于地下的真实介
    质,那么对于该类介质的PP波、PSV波弹性阻抗以及流体识别因子的研究是非常有意义的。

    为达到上述目的,本发明的技术方案为:

    等效孔隙裂缝介质的弹性阻抗及广义流体因子分析方法,由如下步骤构成:

    步骤一、以低频孔隙裂缝模型为等效孔隙裂缝模型,将其与线性滑动模型等价比
    较后建立的新柔量表达式引入以柔量为参数的广义各向异性参数中,得到新广义各向异性
    参数关系式(1)~(4)





    其中,e为裂缝密度,g为横、纵波速度平方比,q为与流体相关的系数;

    步骤二、将等效孔隙裂缝介质的广义Thomsen各向异性参数关系式(1)~(4)代入
    陈天胜推导的HTI介质PP波弹性阻抗公式中,整理得到等效孔隙裂缝介质的PP波弹性阻抗
    公式如下:


    公式中的q与(1)式中相同,将q的Kf替换为气体体积模量,即Kf=0.02GPa,变为气
    饱和裂缝介质的弹性阻抗,Kf替换为水体积模量,即Kf=2GPa,即为水饱和裂缝介质的弹性
    阻抗;

    步骤三、把关系式(1)~(4)代入HTI介质PSV波弹性阻抗公式中,整理得到等效孔
    隙裂缝介质的PSV波弹性阻抗公式,


    通过代入气和水的体积模量,能够得到气饱和与水饱和裂缝介质的PSV波弹性阻
    抗。

    进一步的,选用Hilterman的第1类AVO气砂岩的P波速度、S波速度和密度作为背景
    参数,修改砂岩为基质孔隙度在0~10%范围内的等效孔隙裂缝介质,在裂缝密度e(模型
    ①)、基质孔隙度φp(模型②)、横纵波速度平方比g(模型③)3方面进行参数变化的改造,并
    考虑水饱和与气饱和条件的影响(见表1),开展PP波和PSV波弹性阻抗的数值模拟;

    进一步的,所述数值模拟的具体流程为:

    第一步:计算每种改造对应的归一化各向同性弹性阻抗;

    第二步:计算气饱和与水饱和等效孔隙裂缝介质的归一化弹性阻抗;

    第三步:将第二步中弹性阻抗结果与相应第一步中的各向同性弹性阻抗

    结果做差,得到相应流体饱和状态下裂缝介质的阻抗差;

    第四步:选择入射角θ=30°,方位角提取反映方位信息的阻抗差。

    进一步的,上述等效孔隙裂缝介质的广义流体因子构建方法为:

    为实现裂缝型介质的流体识别,尝试用裂缝介质的PP波和转PSV波的弹性阻抗替
    代Goodwayg构建的λρ和μρ、Russell构建的ρf、宁忠华构建的σHSFIF和σHSFIF1、王栋构建的FIFW
    和F1、以及许平构建的FIFP这几种常规流体因子公式中的纵波和横波阻抗,将裂缝密度、入
    射角和方位角参数引入其中,对应形成g_λρ、g_μρ、g_ρf、g_σHSFIF、g_σHSFIF1、g_FIFW、g_F1和
    g_FIFP广义流体识别因子,具体形式如下:









    相对于现有技术,本发明的有益效果为:

    本发明借助等效裂缝模型中“柔量参数”这个桥梁架起等效孔隙裂缝模型中裂缝
    参数、基质孔隙度与广义各向异性之间的关系,建立新的弹性阻抗与流体因子,并对其中的
    关键参数,包括裂缝密度、纵横波速度比和基质孔隙度进行数值模拟讨论,所得结论如下:

    1)通过直接分析裂缝密度、裂缝填充物、基岩的纵横波速度比、基质孔隙度对弹性
    阻抗的影响,避免了由各向异性参数向裂缝参数和基质孔隙度参数转换过程中引入的误
    差;

    2)裂缝密度、基质孔隙度或横纵波速度平方比参数变化引起的阻抗差变化,可以
    用于区分裂缝型储层的填充物性质,其中裂缝密度和基质孔隙度引起含气介质的阻抗差变
    化高于含水介质,横纵波速度平方比引起的变化则相反;

    3)将常规流体因子延伸的裂缝介质流体因子按流体识别灵敏程度归为3类,其中
    g-FIFW>g-σHSFIF≈g-σHSFIF1,在应用于实际数据处理中具有很好的适用性。

    等效裂缝介质理论作为裂缝介质的一种近似理论,更加接近于地下的真实介质,
    那么对于该类介质的PP波、PSV波弹性阻抗以及流体识别因子的研究是非常有意义的。

    附图说明

    图1气饱和与水饱和等效孔隙裂缝介质的PP波阻抗差。

    其中(a)、(b)、(c)为气饱和等效孔隙裂隙介质,(d)、(e)、(f)为水饱和等效孔隙裂
    隙介质。

    图2气饱和与水饱和等效孔隙裂缝介质的PSV波阻抗差。

    其中(a)、(b)、(c)为气饱和等效孔隙裂隙介质,(d)、(e)、(f)为水饱和等效孔隙裂
    隙介质。

    图3气饱和与水饱和等效孔隙裂缝介质的PP波阻抗差极坐标图。

    其中(a)、(b)、(c)为气饱和等效孔隙裂隙介质,(d)、(e)、(f)为水饱和等效孔隙裂
    隙介质。

    图4气饱和与水饱和等效孔隙裂缝介质的PSV波阻抗差极坐标图。

    其中(a)、(b)、(c)为气饱和等效孔隙裂隙介质,(d)、(e)、(f)为水饱和等效孔隙裂
    隙介质。

    图5不同流体因子的流体识别对比。

    图6气、水饱和裂缝介质流体因子之间的相对差异随入射角变化的曲线图。

    图7气、水饱和裂缝介质流体因子之间的相对差异随方位角变化的曲线图。

    图8测井曲线。

    图9广义流体因子与对应常规流体因子的流体识别效果对比。(灰色横线为干层/
    气层分界线)。

    具体实施方式

    下面结合附图和具体实施方式对本发明技术方案做进一步详细描述:

    等效孔隙裂缝介质的弹性阻抗及广义流体因子分析方法,由如下步骤构成:

    步骤一、以低频孔隙裂缝模型为等效孔隙裂缝模型,将其与线性滑动模型等价比
    较后建立的新柔量表达式引入以柔量为参数的广义各向异性参数中,得到新广义各向异性
    参数关系式(1)~(4)





    其中,e为裂缝密度,g为横、纵波速度平方比,q为与流体相关的系数;

    步骤二、将等效孔隙裂缝介质的广义Thomsen各向异性参数关系式(1)~(4)代入
    陈天胜推导的HTI介质PP波弹性阻抗公式中,整理得到等效孔隙裂缝介质的PP波弹性阻抗
    公式如下:


    公式中的q与(1)式中相同,将q的Kf替换为气体体积模量,即Kf=0.02GPa,变为
    气饱和裂缝介质的弹性阻抗,Kf替换为水体积模量,即Kf=2GPa,即为水饱和裂缝介质的弹
    性阻抗;

    步骤三、把关系式(1)~(4)代入HTI介质PSV波弹性阻抗公式中,整理得到等效孔
    隙裂缝介质的PSV波弹性阻抗公式,


    通过代入气和水的体积模量,得到气饱和与水饱和裂缝介质的PSV波弹性阻抗。

    进一步的,选用Hilterman的第1类AVO气砂岩的P波速度、S波速度和密度作为背景
    参数,修改砂岩为基质孔隙度在0~10%范围内的等效孔隙裂缝介质,在裂缝密度e(模型
    ①)、基质孔隙度φp(模型②)、横纵波速度平方比g(模型③)3方面进行参数变化的改造,并
    考虑水饱和与气饱和条件的影响(见表1),开展PP波和PSV波弹性阻抗的数值模拟;

    进一步的,所述数值模拟的具体流程为:

    第一步:计算每种改造对应的归一化各向同性弹性阻抗;

    第二步:计算气饱和与水饱和等效孔隙裂缝介质的归一化弹性阻抗;

    第三步:将第二步中弹性阻抗结果与相应第一步中的各向同性弹性阻抗结果做
    差,得到相应流体饱和状态下裂缝介质的阻抗差;

    第四步:选择入射角θ=30°,方位角φ∈(0,360°),得到包含方位信息的阻抗差;

    进一步的,上述等效孔隙裂缝介质的广义流体因子构建方法为:

    为实现裂缝型介质的流体识别,尝试用裂缝介质的PP波和转PSV波的弹性阻抗替
    代Goodwayg构建的λρ和μρ、Russell构建的ρf、宁忠华构建的σHSFIF和σHSFIF1、王栋构建的FIFW
    和F1、以及许平构建的FIFP这几种常规流体因子公式中的纵波和横波阻抗,将裂缝密度、入
    射角和方位角参数引入其中,对应形成g_λρ、g_μρ、g_ρf、g_σHSFIF、g_σHSFIF1、g_FIFW、g_F1和
    g_FIFP广义流体识别因子,具体形式如下:









    实验例:

    1等效裂缝介质理论

    根据等效介质理论,许多研究者提出了多种HTI介质模型,常用的主要有三类,孤
    立平行薄扁平裂缝模型(Hudson裂缝模型)、线性滑动模型、等效孔隙裂缝模型,其中等效孔
    隙模型中孔隙和裂缝连通,处于前两种模型的中间状态。它假定裂缝为硬币状,平行排列且
    稀疏分布于由均匀颗粒构成的多孔介质中。不同频带条件下流体影响因子不同,导致了各
    向异性参数值有差异。低频带内,裂缝与相邻孔隙间的流体保持局部压力均衡,即裂缝间的
    流体通过孔隙可以互相流通;高频带内,裂缝间的流体不能相互流通,即各个裂缝独立且相
    互没有影响。地震频带属于低频带,是我们通常研究最多的频带范围,这种低频限制下,如
    果基质孔隙度较低(φp<10%),那么如同模拟平行裂缝一样,可以把介质模型中的孔隙看
    作少量集中分布,即介质的基质孔隙可模拟成各向同性固体背景中的球体稀疏分布,本文
    把这种低频裂缝模型称为等效孔隙裂缝模型,将其与线性滑动模型等价比较后建立的新柔
    量表达式引入以柔量为参数的广义各向异性参数中,得到新广义各向异性参数关系式(1)
    ~(4)





    其中,e为裂缝密度,g为横、纵波速度平方比,q为与流体相关的系数。

    显然,由于孔隙的存在,使得广义各向异性参数变得较为复杂。公式中,γ(V)不受
    填充流体性质的影响,并与裂缝密度e呈线性关系;参数ε(V)、δ(V)和η(V)除与裂缝密度e和横、
    纵波速度平方比g有关之外,还与q系数中的基质孔隙度φp和裂缝填充物相关,而对q系数
    中的纵横比c/a不敏感。

    上述分析可以推论,裂缝填充物、裂缝密度、围岩的横纵波速度比、基质孔隙度是
    直接影响裂缝介质各向异性参数的因素,而作为裂缝型油气藏检测的重要指标,它们必将
    对各向异性弹性阻抗产生显著影响。本文将所推导的广义各向异性参数关系式(1)~(4)引
    入到裂缝介质弹性阻抗公式当中,进行数值模拟,探索对各向异性参数有影响的几种因素
    变化时,等效孔隙裂缝介质弹性阻抗的变化规律。

    2等效孔隙裂缝介质弹性阻抗

    2.1弹性阻抗

    将等效孔隙裂缝介质的广义Thomsen各向异性参数关系式(1)~(4)代入陈天胜推
    导的HTI介质PP波弹性阻抗公式中,整理得到等效孔隙裂缝介质的PP波弹性阻抗公式,


    公式中的q与(1)式中相同,将q的Kf替换为气体体积模量,即Kf=0.02GPa,变为
    气饱和裂缝介质的弹性阻抗,Kf替换为水体积模量,即Kf=2GPa,即为水饱和裂缝介质的弹
    性阻抗。

    同理,把关系式(1)~(4)代入HTI介质PSV波弹性阻抗公式中,整理得到等效孔隙
    裂缝介质的PSV波弹性阻抗公式,


    通过代入气和水的体积模量,得到气饱和与水饱和裂缝介质的PSV波弹性阻抗。

    2.2数值模拟

    选用Hilterman的第1类AVO气砂岩作为背景参数,修改砂岩为等效孔隙度在0~
    10%范围内的等效孔隙裂缝介质,在裂缝密度e(模型①)、基质孔隙度φp(模型②)、横纵波
    速度平方比g(模型③)3方面进行参数改造,并考虑水饱和与气饱和条件的影响(见表1),开
    展PP波和PSV波弹性阻抗的数值模拟,具体流程为:

    第一步:计算每种改造对应的归一化各向同性弹性阻抗;

    第二步:计算气饱和与水饱和等效孔隙裂缝介质的归一化弹性阻抗;

    第三步:将第二步中弹性阻抗结果与相应第一步中的各向同性弹性阻抗结果做
    差,得到

    相应流体饱和状态下裂缝介质的阻抗差。

    模拟结果的组织方式参考表1,所得PP波的阻抗差见图1,PSV波的阻抗差见图2。其
    中,(a)、(d)讨论模型①中裂缝密度为0.05、0.1和0.15时的响应;(b)、(e)讨论模型②中基
    质孔隙度φp为0.001、0.01和0.1时的响应;(c)、(f)讨论模型③中横纵波速度平方比g为
    0.25、0.36和0.49时的响应。每图中曲面从上到下对应每种参数从小到大而变化。

    表1双层模型介质参数设置


    为分析阻抗的方位各向异性,进一步抽取在入射角θ=30°,方位角时
    相应PP波和PSV波阻抗差绝对值于极坐标中显示为图3和图4,其中,(a)和(d)对应模型①的
    裂缝密度e参数变化响应,(b)和(e)对应模型②的等径孔隙度φp参数变化响应,(c)和(f)
    对应模型③的横纵波速度平方比g参数变化响应,实线、点划线和点线代表对应参数依次增
    大。

    综合分析图1~图4,总结如下:

    1)PP波、PSV波阻抗差随着入射角、方位角的变化呈现周期性的方位各向异性特征
    (图3与图4),但有所区别。区别之一:曲线、曲面形态。PP波阻抗差方位角在垂直裂缝面方向
    (0°和180°)的相邻区域内,阻抗差幅值最大,各向异性最强,方位角在平行裂缝面方向
    (90°、270°)的相邻区域内,阻抗差为0,即各向同性。而PS波阻抗差呈4片花瓣状,方位角在
    垂直裂缝面方向(0°和180°)和平行裂缝面方向(90°、270°)的相邻区域内,都具有方位各向
    异性差异;在平行裂缝面方向,同一参数变化情况下,气、水饱和状态下阻抗差幅值相等;在
    垂直裂缝面方向,阻抗差幅值存在显著的梯度变化,是PSV波阻抗差各向异性差异性的主要
    来源。区别之二:阻抗差数值分布。PP波阻抗差始终为负,且在单一象限内数值幅度随方位
    角增大、随入射角增大呈现稳定的单调递增特征,这表明HTI介质PP波阻抗在数值上小于各
    向同性阻抗。而PS波阻抗差随入射角、方位角变化会出现符号反转现象,换言之PSV波各向
    异性阻抗值不总是低于各向同性阻抗。

    2)裂缝密度是影响阻抗的最重要因素,无论是气饱和还是水饱和条件下,由裂缝
    密度引起的阻抗差变化幅度都要高于基质孔隙度和横、纵波速度平方比参数。另外,对比图
    1、2、3、4中图(a)、(d)可见,PP波、PSV波阻抗差均随裂缝密度增大而增大,其中,PP波为负
    向,PSV波为正向。对比阻抗差的幅值,可见气饱和状态均高于水饱和状态。

    3)横、纵波速度平方比对阻抗的影响程度略弱于裂缝密度。对比图1、2、3、4中图
    (c)、(f)可见,PP波、PSV波阻抗差均随速度平方比的增大而增大,其中,PP波为负向,PSV波
    为正向。对比阻抗差幅值的梯度变化可见,PP波阻抗差在气饱和状态下均弱于水饱和状态,
    而PSV波阻抗差的变化恰恰相反,在气饱和状态下强于水饱和状态。

    4)基质孔隙度是阻抗的另外一个影响因素,对比图1、2、3、4中图(b)、(e)可见,在
    孔隙介质中,随着基质孔隙度的增加,PP波、PSV波阻抗差均有增大的趋势,且气饱和条件下
    的阻抗差梯度变化更加显著。

    可见,阻抗差可明显区分填充物性质,识别介质含气或含水性。裂缝密度、基质孔
    隙度或横纵波速度平方比的变化对于区分裂缝型储层的填充物性质也具有很好的优势。

    3广义流体因子构建及敏感性分析

    以声波的纵波、横波阻抗组合而成的几种常见流体识别因子主要有Goodwayg构建
    的 λρ和μρ、Russell构建的ρf、宁忠华构建的σHSFIF和σHSFIF1、王栋构建的FIFW和F1、以及
    许平构建的FIFP。为实现裂缝型介质的流体识别,尝试用裂缝介质的纵波和转换横波的弹
    性阻抗替代这几种常见流体因子公式中的纵波和横波阻抗,将裂缝密度、入射角和方位角
    参数引入流体因子中,形成广义流体识别因子,形式如下:









    以Hilterman给出的第1类砂岩参数(见表2)作为双层介质的背景参数,将下层砂
    岩修改为等效孔隙裂缝介质,取裂缝密度e=0.1,进行不同流体饱和状态下介质的PP波和
    PSV波弹性阻抗计算,进而获得几种广义流体因子,探讨它们对裂缝介质流体识别的适用
    性。

    表2模型二介质背景参数


    抽取入射角θ=30°,方位角时的流体因子曲线如图5所示,图中增加
    g-EIp、g-EIs以示表明纵波、转换横波弹性阻抗作为流体指示因子效果很差,不提倡直接用
    于流体识别。g-μρ在水气两种饱和状态下差异不明显,而g-λρ差异较显著,这与宁忠华于
    2006年所述观点一致。g-F1突跳点太多且值很大,无法实现整体方位角范围的流体识别,不
    适于裂缝介质模型,其它广义流体因子气、水饱和差异较大。

    由上述可见,g-μρ、g-λρ、g-ρf、g-σHSFIF、g-σHSFIF1、g-FIFW和g-FIFP均适用于裂
    缝介质的流体识别,为比较它们灵敏度差异性,引入一个相对差异参数,即水、气饱和状 态
    下流体因子之差的绝对值除以两者之和的绝对值,相对差异值越大,说明其流体识别的灵
    敏度越高。设裂缝密度不变,固定方位角入射角取θ∈(0°,30°),
    得到这些广义流体因子的相对差异随入射角变化曲线如图6所示,固定入射角θ=[10°,
    20°,30°],方位角取则随方位角变化的曲线如图7所示。图6中,为各向同
    性面,不含流体,流体因子相对差异为零,和时相对差异曲线基本相近,只保留
    其一进行讨论即可。将图6中30°,45°,60°与图7综合分析,根据灵敏度差异程度,可将
    流体因子g-FIFW和g-FIFP归为一类,记为第1类灵敏度流体因子,g-ρf和g-σHSFIF归为一
    类,记为第2类灵敏度流体因子,g-λρ和g-σHSFIF1归为一类,记为第3类灵敏度流体因子;g-
    μρ流体因子相对差异值不超过0.05,灵敏度非常低,单独用以识别流体效果极差。

    虽然3类灵敏度流体因子均可明显区分介质中的流体性质,但顾及到在方位角和
    入射角两种范围内相对差异值的覆盖程度,申请人认为第1类灵敏度流体因子最佳,即g-
    FIFW区分流体的能力更强。为验证这一观点,下面就以实际数据为例分析这3类灵敏度流体
    因子的识别效果以及探讨最优灵敏度流体因子。

    4应用实例

    已知某工区的测井曲线如图8所示,深度在5046m以上的厚度为各向同性层,在
    5046m~5048m为干层,裂缝密度为0.03,5048m以下为含气层,裂缝密度为0.2。将g-FIFW,g-
    σHSFIF和g-σHSFIF1这3类流体因子应用于该测井数据处理中,并与相应常规流体因子进行
    对比,结果如图9所示,实线代表入射角和方位角均为30°时的广义流体因子,虚线代表常规
    流体因子。从图中可知,广义因子与常规因子的曲线形态基本相似,说明使用广义流体因子
    进行流体识别具有较好可信度;常规流体因子FIFW,σHSFIF和σHSFIF1基本无法区分干层和
    含气层,而广义流体因子从气层上界面(5048m)开始出现明显的变化;3类广义流体因子与
    对应常规流体因子对比,在气层,g-FIFW与对应常规流体因子的差异性最大,g-σHSFIF次
    之,g-σHSFIF1差异性最小,由此验证了推导公式的正确性。

    5结论

    裂缝流体研究,是世界剩余油气研究的重要领域。本申请借助等效裂缝模型中“柔
    量参数”这个桥梁架起等效孔隙裂缝模型中裂缝参数、基质孔隙度与广义各向异性之间的
    关系,建立新的弹性阻抗与流体因子,并对其中的关键参数,包括裂缝密度、纵横波速度比
    和基质孔隙度进行数值模拟讨论,所得结论如下:

    1)等效孔隙裂缝介质的纵波、转换横波弹性阻抗呈现明显的方位各向异性,且两
    者差异 明显;在垂直裂缝方向,两者都具有明显的方位性,而沿裂缝方向转换横波存在方
    位各向异性,纵波却不具有这一特点;

    2)裂缝密度、基质孔隙度或横纵波速度平方比参数变化引起的阻抗差变化,可以
    用于区分裂缝型储层的填充物性质,其中裂缝密度和基质孔隙度引起含气介质的阻抗差变
    化高于含水介质,横纵波速度平方比引起的变化则相反;

    3)流体因子源于纵、横波弹性阻抗,由裂缝介质的转换横波弹性阻抗单独形成的
    流体因子g-μρ不适于进行介质中的流体识别;

    4)将常规流体因子延伸的裂缝介质流体因子按流体识别灵敏程度归为3类,其中
    g-FIFW>g-σHSFIF≈g-σHSFIF1,在应用于实际数据处理中具有很好的适用性。

    等效裂缝介质理论作为裂缝介质的一种近似理论,更加接近于地下的真实介质,
    那么对于该类介质的PP波、PSV波弹性阻抗以及流体识别因子的研究是非常有意义的。

    以上所述,仅为本发明的具体实施方式,但本发明的?;し段Р⒉痪窒抻诖?,任何
    不经过创造性劳动想到的变化或替换,都应涵盖在本发明的?;し段е?。因此,本发明的
    ?;し段вΩ靡匀ɡ笫樗薅ǖ谋;し段?。

    关于本文
    本文标题:等效孔隙裂缝介质的弹性阻抗及广义流体因子分析方法.pdf
    链接地址://www.4mum.com.cn/p-6001500.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • 助赢彩票软件 手机 新浪双色球基本走势图 一全网最早原创36码特开 北京pk赛车龙虎技巧论坛 重庆时时全天计划单期 投注技巧规律码 今晚七星彩开奖结果是 北京时时是正规的 购彩票大厅 pc28怎么才能稳赚 北京飞艇计划赛车全天计划网页 黑龙江时时开奖号码0 手机赌博赢了几十万 dnf铁马有多稳 江苏快三2期计划软件 双色球投注单详解