[发明专利]含气饱和度确定方法及装置有效
申请号: | 201611263196.9 | 申请日: | 2016-12-30 |
公开(公告)号: | CN106842303B | 公开(公告)日: | 2019-07-05 |
发明(设计)人: | 谢伟彪;司兆伟;庄东志;徐风;张建林;陈晶莹;马明明;殷秋丽 | 申请(专利权)人: | 中国石油天然气股份有限公司 |
主分类号: | G01V1/28 | 分类号: | G01V1/28;G01V1/30 |
代理公司: | 北京三高永信知识产权代理有限责任公司 11138 | 代理人: | 董亚军 |
地址: | 100007 北京市*** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 饱和度 确定 方法 装置 | ||
1.一种含气饱和度确定方法,其特征在于,所述方法包括:
分别获取纵波和横波在至少两个岩样中传播的至少八个样本波速参数,所述至少两个岩样中的每个岩样对应至少两个含气饱和度,所述至少八个样本波速参数包括至少四个横波样本参数和至少四个纵波样本参数,每个岩样对应至少两个横波样本参数和至少两个纵波样本参数;
采用波速-时差转换公式对所述至少八个样本波速参数中的每个样本波速参数进行转换,得到至少八个样本时差,所述波速-时差转换公式为:
其中,DT为样本时差,v为样本波速参数;
根据所述至少四个纵波样本时差,采用模型参数计算公式得到所述模型参数,所述模型参数计算公式为:
其中,ΔDT为所述模型参数,DTC为预设岩样处于不同含水饱和度时对应的纵波时差,DTCsg=0为所述预设岩样为第一岩样时对应的第一纵波样本时差,DTCsg=1为所述预设岩样为第二岩样时对应的第二纵波样本时差,所述预设岩样为所述至少一个岩样中的任意一个岩样;
根据每个岩样对应的模型参数和每个岩样对应的含气饱和度,采用最小二乘法确定所述含气饱和度计算公式,所述含气饱和度计算公式为:
Sg∝ΔDT
其中,Sg为含气饱和度,ΔDT为模型参数;
分别获取横波和纵波在待确定气藏储层中的传播速度的参数,所述传播速度的参数包括横波速度参数和纵波速度参数;
根据所述横波速度参数以及预先确定的转换参数,分别确定第一纵波时差和第二纵波时差,所述第一纵波时差是纵波在含气饱和度为0的待确定气藏储层中传播时的纵波时差,所述第二纵波时差是纵波在含气饱和度为1的待确定气藏储层中传播时的纵波时差;
根据所述纵波速度参数、所述第一纵波时差、所述第二纵波时差和预先确定的模型参数,确定所述待确定气藏储层的含气饱和度;
所述根据所述纵波速度参数、所述第一纵波时差、所述第二纵波时差和预先确定的模型参数,确定所述待确定气藏储层的含气饱和度,包括:
将所述纵波速度参数、所述第一纵波时差、所述第二纵波时差和所述模型参数代入预先确定的含气饱和度计算公式,得到所述待确定气藏储层的含气饱和度。
2.根据权利要求1所述的方法,其特征在于,所述至少八个样本时差包括与所述至少四个横波样本参数对应的至少四个横波样本时差,以及与所述至少四个纵波样本参数对应的至少四个纵波样本时差;
在所述采用波速-时差转换公式对所述至少八个样本波速参数中的每个样本波速参数进行转换之后,所述方法还包括:
根据所述至少四个横波样本时差和所述至少四个纵波样本时差,采用最小二乘法确定所述转换参数。
3.根据权利要求1所述的方法,其特征在于,每个岩样包括含气饱和度为0时的第一岩样以及含气饱和度为1时的第二岩样,所述至少四个纵波样本参数包括纵波在所述第一岩样中传播的至少两个第一纵波样本参数,以及在所述第二岩样中传播的至少两个第二纵波样本参数,所述至少八个样本时差包括与所述至少四个纵波样本参数对应的至少四个纵波样本时差,所述至少四个纵波样本时差包括分别与所述至少两个第一纵波样本参数对应的至少两个第一纵波样本时差,以及与所述至少两个第二纵波样本参数对应的至少两个第二纵波样本时差。
4.根据权利要求3所述的方法,其特征在于,在所述分别获取纵波和横波在至少一个岩样中传播的至少八个样本波速参数之前,所述方法还包括:
获取所述至少一个岩样中每个岩样对应的含水饱和度;
根据每个岩样对应的含水饱和度,采用饱和度转换公式计算得到每个岩样的含气饱和度,所述饱和度转换公式为:
Sg=1-Sw
其中,Sg为含气饱和度,Sw为含水饱和度。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国石油天然气股份有限公司,未经中国石油天然气股份有限公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201611263196.9/1.html,转载请声明来源钻瓜专利网。