[发明专利]一种流体识别方法及系统有效
申请号: | 201410634719.0 | 申请日: | 2014-11-12 |
公开(公告)号: | CN104316966A | 公开(公告)日: | 2015-01-28 |
发明(设计)人: | 吴国忱;印兴耀;郭俊超;宗兆云 | 申请(专利权)人: | 中国石油大学(华东) |
主分类号: | G01V1/30 | 分类号: | G01V1/30;G06F19/00 |
代理公司: | 北京集佳知识产权代理有限公司 11227 | 代理人: | 王宝筠 |
地址: | 266580 山东省*** | 国省代码: | 山东;37 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 流体 识别 方法 系统 | ||
1.一种流体识别方法,其特征在于,该方法包括:
将经过保衰减信息处理的叠前CMP道集的数据转化为叠前CRP道集的数据;
根据所述叠前CRP道集提取至少3个角度的叠加数据;
根据所述叠加数据计算粘弹性流体因子、剪切模量和密度;
利用所述粘弹性流体因子、所述剪切模量和所述密度对储层流体进行识别。
2.根据权利要求1所述的流体识别方法,其特征在于,所述计算粘弹性流体因子、剪切模量和密度包括:
建立粘弹性阻抗模型;
根据所述粘弹性阻抗反演方法进行衰减粘弹性流体因子反演,得到所述粘弹性流体因子、所述剪切模量和所述密度。
3.根据权利要求2所述的流体识别方法,其特征在于,所述建立粘弹性阻抗模型包括:
建立粘弹性介质的纵波反射系数方程;
利用弱粘弹性和相似介质假设,得到所述纵波反射系数近似表达式;
利用所述纵波反射系数近似表达式建立计算粘弹性流体因子、剪切模量和体积密度的方程。
4.根据权利要求3所述的流体识别方法,其特征在于,所述建立粘弹性介质的纵波反射系数方程具体包括:
确定粘弹性流体因子表达式;
确定非弹性介质纵波反射系数表达式;
将所述粘弹性流体因子表达式带入到所述非弹性介质纵波反射系数表达式中,得到用所述粘弹性流体因子表示的纵波反射系数方程。
5.根据权利要求2所述的流体识别方法,其特征在于,所述根据粘弹性阻抗反演方法进行衰减粘弹性流体因子反演,得到所述粘弹性流体因子、所述剪切模量和所述密度包括:
对所述粘弹性阻抗模型进行线性变换;
获取所述至少3个角度对应的粘弹性阻抗数据;
根据所述粘弹性阻抗数据和测井曲线,计算得到所述粘弹性流体因子、所述剪切模量和所述密度。
6.根据权利要求5所述的流体识别方法,其特征在于,所述获取至少3个角度对应的粘弹性阻抗数据包括:
提取角度子波;
进行粘弹性阻抗反演,得到至少3个角度对应的粘弹性阻抗数据。
7.根据权利要求6所述的流体识别方法,其特征在于,所述提取角度子波的提取条件为:
E=∑(si-di)2
其中,E为误差能量,s=r*w为合成地震记录,r是反射系数,w是提取的子波,d为地震记录,选择最佳的子波,同时使误差能量最小。
8.根据权利要求7所述的流体识别方法,其特征在于,所述进行粘弹性阻抗反演,得到至少3个角度对应的粘弹性阻抗数据具体为:
用稀疏脉冲反演方法反演各角度的粘弹性阻抗数据体,其目标优化函数为:
F=Lp(r(θ))+λLq(s-d)+α-1(ΔZtrend)
其中,r(θ)为角度反射系数序列,ΔZ为与粘弹性阻抗趋势的差序列,d为地震道序列,s为合成地震道序列,λ为残差权重因子,α为趋势权重因子,p和q为L模因子。
9.根据权利要求1所述的流体识别方法,其特征在于,所述经过保衰减信息处理的叠前CMP道集的数据具体为包含了受粘弹性影响的叠前CMP道集数据。
10.根据权利要求1所述的流体识别方法,其特征在于,所述根据叠前CRP道集提取至少3个角度的叠加数据具体为根据所述叠前CRP道集提取小角度、中角度和大角度的叠加数据。
11.一种流体识别系统,其特征在于,包括:
转化单元:用于将经过保衰减信息处理的叠前CMP道集的数据转化为叠前CRP道集的数据;
第一提取单元:用于根据所述叠前CRP道集提取至少3个角度的叠加数据;
第一计算单元:用于根据所述叠加数据计算粘弹性流体因子、剪切模量和密度;
识别单元:用于利用所述粘弹性流体因子、所述剪切模量和所述密度对储层流体进行识别。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国石油大学(华东),未经中国石油大学(华东)许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201410634719.0/1.html,转载请声明来源钻瓜专利网。
- 上一篇:一种页岩气储层的岩石模量计算方法
- 下一篇:电缆护套橡胶及其制备方法