[发明专利]一种低渗气藏压裂水平井产出剖面反演解释方法有效
申请号: | 202010643544.5 | 申请日: | 2020-07-07 |
公开(公告)号: | CN112302606B | 公开(公告)日: | 2021-08-24 |
发明(设计)人: | 罗红文;李海涛;李颖;蒋贝贝;崔小江;刘畅;朱晓萍;邹顺良;李保吉 | 申请(专利权)人: | 西南石油大学 |
主分类号: | E21B43/26 | 分类号: | E21B43/26;E21B49/00;G06F30/20;G06Q50/02 |
代理公司: | 重庆嘉禾共聚知识产权代理事务所(普通合伙) 50220 | 代理人: | 吴迪 |
地址: | 610500 四*** | 国省代码: | 四川;51 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 低渗气藏压裂 水平 产出 剖面 反演 解释 方法 | ||
1.一种低渗气藏压裂水平井产出剖面反演解释方法,其特征在于,包括以下步骤:
S1、根据目标压裂水平井的实测温度剖面识别和定位有效裂缝,初始化反演系统温度并根据实测温度剖面在各裂缝位置处温度降的温度降ΔT初始化裂缝半长分布并初始化Markov链长为LN;
S2、当前反演系统温度下的第n次反演迭代时,将当前裂缝半长分布代入温度正演预测模型,模拟计算目标压裂水平井生产时的温度剖面然后采用误差方程计算当前拟合评价目标函数
S3、通过扰动方程对当前裂缝半长分布产生随机扰动,生成新的裂缝半长分布并再通过误差方程计算新的拟合评价目标函数
S4、将新的拟合评价目标函数和当前拟合评价目标函数进行比较;若则将步骤S3中新的裂缝半长分布作为当前解;
若则通过概率方程计算新的裂缝半长分布的接受概率,然后在[0,1)区间生成一个均匀分布的随机数λ,若随机数λ小于接受概率,则将新的裂缝半长分布作为当前解,否则将作为当前解;
S5、在当前反演系统温度下转步骤S2,重复迭代LN次,直至生成一个长度为LN的Markov链;
S6、判断当前拟合评价目标函数是否满足收敛方程,若当前拟合评价目标函数满足所述的收敛方程,则输出裂缝半长分布当前反演解;否则根据退火方程计算新的反演系统温度并重新初始化Markov链长LN,然后转步骤S2,直至拟合评价目标函数满足所述的收敛方程,然后输出裂缝半长分布当前反演解;
S7、验证裂缝半长分布当前反演解对应的模拟产量是否收敛于目标压裂水平井的实测产量,若模拟产量收敛于实测产量,则反演完成;否则转步骤S1,直至反演完成,输出裂缝半长分布反演解;
S8、将获得的裂缝半长分布反演解输入温度正演预测模型,计算出目标压裂水平井的各级裂缝流量和产出剖面反演解释结果。
2.根据权利要求1所述的一种低渗气藏压裂水平井产出剖面反演解释方法,其特征在于,所述温度正演预测模型包括储层渗流模型、储层热学模型、裂缝渗流模型、裂缝热学模型、井筒流动模型、井筒热学模型。
3.根据权利要求1所述的一种低渗气藏压裂水平井产出剖面反演解释方法,其特征在于,所述步骤S2和S3中的误差方程如下:
式中,为拟合评价目标函数,为实测的温度剖面,为模拟计算的温度剖面,上标T表示向量转置。
4.根据权利要求3所述的一种低渗气藏压裂水平井产出剖面反演解释方法,其特征在于,所述步骤S3中的扰动方程如下:
式中,ξ为扰动系数,rand(0,1)表示生成0-1的随机数,为当前裂缝半长分布;为新的裂缝半长分布。
5.根据权利要求1所述的一种低渗气藏压裂水平井产出剖面反演解释方法,其特征在于,所述步骤S4中的概率方程如下:
式中,为新生成的裂缝半长分布的接受概率,min{}表示取最小值,表示生成的建议分布,表示生成的建议分布。
6.根据权利要求1所述的一种低渗气藏压裂水平井产出剖面反演解释方法,其特征在于,所述步骤S6中的收敛方程如下:
式中,ε为反演误差精度;为拟合评价目标函数。
7.根据权利要求6所述的一种低渗气藏压裂水平井产出剖面反演解释方法,其特征在于,所述步骤S6中的退火方程如下:
式中,α为温度衰减率,k=0,1,2,…,n;为当前反演系统温度;为新的反演系统温度。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于西南石油大学,未经西南石油大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202010643544.5/1.html,转载请声明来源钻瓜专利网。