一种用于确定剩余油饱和度的碳氧比值计算方法
【技术领域】
[0001] 本发明涉及石油天然气开发领域，具体涉及一种针对碳氧比测井的碳氧比值的方 法。
【背景技术】
[0002] 碳氧比能谱测井通过记录中子与地层元素原子核作用产生的非弹散射及俘获伽 马射线，能够在未知矿化度及矿化度变化较大的地层条件下定量评价储层剩余油饱和度， 在油田开发中后期的剩余油分布度动态监测、提高采收率方面发挥了重要作用。目前，通过 处理测量非弹及俘获伽马能谱获取碳氧比值的主要方法是能窗法。
[0003] 单能伽马射线的探测器响应能谱，主要由全能峰及康普顿平台等部分组成。利用 能窗法计算碳氧比值时，碳能窗伽马计数会受其他高能伽马射线康普顿平台计数的影响， 降低了复杂岩性及低孔隙度储层条件下，碳氧比测井响应灵敏度及饱和度计算精度。
【发明内容】
[0004] 针对上述能窗法计算碳氧比值相对含油饱和度变化灵敏度低的问题，本发明提出 了一种用于确定剩余油饱和度的碳氧比值计算方法，该方法能够降低干扰伽马对碳氧比值 影响，明显提高低孔隙度条件下碳氧比值对含油饱和度的响应灵敏度。
[0005] 其技术解决方案包括：
[0006] -种用于确定剩余油饱和度的碳氧比值计算方法，依次包括以下步骤：
[0007] a利用高斯函数及线性函数组合模型分别对测量非弹谱及标准伽马能谱的碳、氧 能峰进行非线性拟合，使用Levenberg-Marquardt方法计算得到对应的特征系数，采用的拟 合模型如下所示：
[0008]
(1)
[0009] C2 = dXx+e (2)
[0010] C = Ci+C2 (3)
[0011] 式中：Ci为拟合能峰的高斯部分计数;C2为拟合能峰的线性部分计数;a，b，c分别为 高斯部分的高度、中心位置及标准方差拟合系数;d，e分别为线性部分的一次系数及常系 数;C为测量伽马能谱的能道计数;X为道址对应能量，单位为MeV;
[0012] 结合式（1)、（2)、（3)计算得到测量非弹谱及标准伽马能谱的碳、氧能峰高斯拟合 系数￡lc、￡l〇、Else及&so ;
[0013] b利用公式(4)、（5)计算碳、氧两种元素产生非弹伽马射线计数，
[0014] ⑷
[0015] (5)
[0016] 式中：c_ntsS测量谱中碳元素产生非弹伽马射线计数，0_nts为测量谱中氧元素 产生非弹伽马射线计数，&、s。分别为标准化后的碳、氧元素非弹标准伽马能谱总计数；
[0017] C根据公式(6)计算得到地层碳氧比值，
[0018]
(6')_。.
[0019] 上述计算方法用于测量地层含油饱和度。
[0020] 本发明首先利用基于高斯模型及线性模型的组合模型拟合实测非弹伽马能谱数 据的碳、氧能峰，获取特征峰拟合系数，特征系数主要包括测量谱碳、氧能峰的的高斯部分 高度、中心位置及标准方差；以及线性组成部分的一次系数及常系数；
[0021] 然后利用相同的非线性拟合模型，拟合碳、氧标准伽马能谱对应能峰，获取其对应 拟合系数，拟合系数同样包括氧能峰的的高斯部分高度、中心位置及标准方差；以及线性组 成部分的一次系数及常系数；
[0022] 利用测量非弹能谱碳峰高斯组成部分的高度与标准谱碳峰高斯组成部分的高度 的比值乘以碳元素标准谱伽马能谱总计数，计算得到碳元素产生非弹伽马计数；
[0023] 利用测量非弹能谱氧峰高斯组成部分的高度与标准谱氧峰高斯组成部分的高度 的比值乘以氧元素标准谱伽马能谱总计数，计算得到氧元素产生非弹伽马计数；
[0024] 最后，利用计算碳元素产生非弹伽马计数除以氧元素非弹伽马计数，生成碳氧比 值，用于计算地层含油饱和度。
[0025]相比现有技术，本发明基于高斯模型及线性模型组合的计算方法，并结合碳、氧元 素标准伽马能谱，降低了其他元素产生伽马射线对计算碳氧比值的影响，提高了碳氧比值 相对含油饱和度响应灵敏度。
【附图说明】
[0026] 下面结合附图对本发明做进一步说明：
[0027] 图1为本发明计算方法确定地层碳氧比值的具体流程图；
[0028] 图2为计算碳氧比过程中采用的碳、氧元素标准伽马能谱图；
[0029] 图3为孔隙度为10%条件下本发明方法与常规能窗法计算碳氧比值对比图；
[0030] 图4为孔隙度为15%条件下本发明方法与常规能窗法计算碳氧比值对比图。
【具体实施方式】
[0031] 下面结合具体实施例对本发明做详细说明。
[0032] 对于碳氧比测井仪器使用的闪烁晶体探测器，具有特定能量的入射伽马射线能量 沉积产生的可见光子数成高斯分布，其平均值为对应入射伽马光子能量;测量谱中某能量 范围内若没有能峰贡献，可以利用线性模型近似反应康普顿及本底贡献。
[0033] 实施例1:
[0034] 本发明，一种碳氧比值的计算方法，如图1所示，具体包括以下步骤：
[0035] 步骤1、利用高斯函数与线性函数组合模型对测量非弹谱的碳、氧特征能峰进行拟 合，使用Levenberg-Marquardt方法计算得到相关参数，采用拟合模型如下所示：
[0036]
[0037]
[0038]
[0039] 式中：Ci为拟合能峰的高斯部分;C2为拟合能峰的线性部分;a，b，c分别为高斯部分 的高度、中心位置及标准方差;d，e分别为线性部分的一次系数及常系数;C为测量伽马能谱 道计数;X为道址对应能量，单位为MeV，拟合完成后得到特征峰非线性拟合系数。
[0040] 利用上述相同方法，对碳、氧标准伽马能谱如图2所示，进行非线性拟合，得到非线 性拟合系数，并计算碳、氧元素标准谱伽马计数。
[0041] 采用上述方法分别计算得到仪器测量谱及标准伽马能谱的碳、氧能峰高斯拟合系 数ac、a。、as。及a s。，测量谱是由单元素标准谱及本底贡献组成，测量谱的碳、氧能峰的高斯组 成部分越大，说明碳、氧元素标准谱所占贡献越大，利用如下公式(4)(5)可以计算碳、氧两 种元素产生非弹伽马射线总计数：
[0042] (4)
[0043] (5)
[0044] 式中：C_nts为测量谱中碳元素产生非弹伽马射线，0__为测量谱中氧元素产生 非弹伽马射线，&、S。分别为标准化后的碳、氧元素非弹标准伽马能谱总计数。进而可以计算 得到地层碳氧比值：
[0045]
(6)
[0046] 利用上述具体实施方案，处理孔隙度为10%及15%条件下，含油饱和度分别为 0%、20%、40%、60%、80%与100%地层的非弹伽马能谱，计算得到碳氧比值;并与常规能 窗法进行对比，得到碳氧比值与含油饱和度关系如图3、4所示。可以看出孔隙度固定时，随 着含油饱和度的增加，两种方法计算碳氧比值都基本呈线性增大。相对常规能窗法，本发明 方法计算碳氧比值变化更加剧烈，对含油饱和度变化响应更加敏感。
[0047] 利用公式(7)、（8)分别计算两种低孔隙度条件下，碳氧比动态变化范围及灵敏度 如表1所示。
[0050]表1碳氧比动态范围及灵敏度
[0048] (7)
[0049]
[0052]与现有能窗法相比，新方法计算碳氧比值具有更大动态变化范围及更高灵敏度； 在孔隙度为10%时，新方法计算碳氧比值的动态变化范围约为常规能窗法计算值的2.5倍； 因此新方法计算碳氧比值在确定低孔隙度储层剩余油饱和度方面具有优势。
【主权项】
1. 一种应用于确定剩余油饱和度的碳氧比值计算方法，其特征在于，依次包括以下步 骤： a利用高斯函数及线性函数组合模型分别对测量非弹谱及标准伽马能谱的碳、氧能峰 进行非线性拟合，使用Levenberg-Marquardt方法计算得到对应的特征系数，采用的拟合模 型如下所示：〇2 = d X x+e (2) C = Ci+C2 (3) 式中：Cl为拟合能峰的高斯部分计数;C2为拟合能峰的线性部分计数;a，b，c分别为高斯 组成部分的高度、中心位置及标准方差拟合系数;d，e分别为线性组成部分的一次系数及常 系数;C为测量伽马能谱的能道计数;X为道址对应能量，单位为MeV; 结合式（1)、（2)、（3)计算得到测量非弹谱及标准伽马能谱的碳、氧能峰高斯拟合系数 Be、&〇、Bsc_S_￡lso ; b利用公式(4)、（5)计算碳、氧两种元素产生非弹伽马射线计数，式中：(：__为测量谱中碳元素产生非弹伽马射线计数，0__为测量谱中氧元素产生 非弹伽马射线计数，&、S。分别为标准化后的碳、氧元素非弹标准伽马能谱总计数； c根据公式(6)计算得到地层碳氧比值，2. 根据权利要求1所述的碳氧比值计算方法，其特征在于:所述计算方法用于测量地层 含油饱和度。
【专利摘要】本发明公开了一种应用于确定剩余油饱和度的碳氧比值计算方法，首先利用高斯与线性双模型拟合实测非弹伽马能谱数据的碳、氧能峰，获取特征峰拟合系数；然后利用相同的拟合模型，拟合碳、氧标准伽马能谱对应能峰；利用测量谱碳峰高斯组成部分的高度与标准谱碳峰高斯组成部分的高度的比值乘以碳元素标准谱总计数，得碳元素产生非弹伽马计数；利用测量谱氧峰高斯组成部分的高度与标准谱氧峰高斯组成部分的高度的比值乘以氧元素标准谱总计数，得氧元素产生非弹伽马计数；最后二者相除生成碳氧比值。该方法降低了其他元素产生伽马射线对碳氧比值的影响，提高了碳氧比值相对含油饱和度的响应灵敏度。
【IPC分类】E21B49/00
【公开号】CN105673006
【申请号】CN201610008112
【发明人】刘军涛, 张锋, 张泉滢, 吴赫, 田立立, 王新光
【申请人】中国石油大学(华东)
【公开日】2016年6月15日
【申请日】2016年1月7日