[发明专利]一种射孔损伤室内综合评价方法在审
申请号: | 202110520310.6 | 申请日: | 2021-05-13 |
公开(公告)号: | CN113189122A | 公开(公告)日: | 2021-07-30 |
发明(设计)人: | 梁豪;吴木旺;裴健翔;郭书生;姜洪丰 | 申请(专利权)人: | 中海石油(中国)有限公司海南分公司;中海石油(中国)有限公司湛江分公司 |
主分类号: | G01N23/046 | 分类号: | G01N23/046;G01N15/08;G06F30/23 |
代理公司: | 北京中索知识产权代理有限公司 11640 | 代理人: | 唐亭 |
地址: | 570000 海南省海口*** | 国省代码: | 海南;46 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 损伤 室内 综合 评价 方法 | ||
1.一种射孔损伤室内综合评价方法,其特征在于,包括:
步骤S10、根据岩心扫描实验确定岩心的渗透率平面分布数据;
步骤S20、根据渗透率平面分布数据、Q型聚类分析法对岩心的原装带、过渡带、破碎带及压实带进行分带划分;
步骤S30、根据分带划分后岩心的原装带、过渡带、破碎带及压实带分别取原状带、射孔尖端破碎带、射孔尖端压实带的岩样;
步骤S40、分别对原状带、射孔尖端破碎带、射孔尖端压实带的岩样进行微CT扫描,并获得CT扫描结果、数字三维图像,再利用图像处理技术对岩样孔隙结构进行还原得到孔隙网络模型;
步骤S50、根据CT扫描结果建立数字岩心,然后将数字岩心的孔隙按孔径分布频率进行十等分,利用分割出10种大小级别的孔隙;
步骤S60、根据数字三维图像统计得到等效孔径的分布与孔隙数量的分布,绘制的数据散点,最后通过这些点的拟合直线的斜率取绝对值,得到岩心的孔隙度分形维数D;
步骤S70、基于孔隙网络模型统计连通孔道的实际长度LT和表观长度L的分布情况,绘制lgLT~lgλ的数据散点,最后通过这些点的拟合线斜率求绝对值,可以得到孔隙的迂曲度分形维数T;
步骤S80、统计计算10个孔隙区间内的区间中值孔隙直径R、孔隙度分形维数D、迂曲度分形维数T、形状因子A、以及孤立孔隙占比E,并计算得到孔隙结构综合因子ω;
步骤S90、最后根据孔隙结构综合因子ω评价岩心的射孔压实和破碎程度。
2.根据权利要求1所述的一种射孔损伤室内综合评价方法,其特征在于,所述步骤S10中的渗透率进行渗透率校正。
3.根据权利要求2所述的一种射孔损伤室内综合评价方法,其特征在于,所述渗透率校正的具体步骤为:
1)利用有限元方法计算得到岩样的固有渗透率kint;
2)根据实际流体的粘度、温度、压力状态,计算流体校正系数β:
3)计算岩心的视渗透率校正系数kr:
4)利用视渗透率校正系数kr,通过下式便可计算得到对应温度、压力和流体性质的表观渗透率k:
k=krkint
式中:k为表观渗透率;kr为视渗透率校正系数;kint为固有渗透率。
4.根据权利要求2所述的一种射孔损伤室内综合评价方法,其特征在于,所述流体校正系数β的计算公式为:
式中:T为温度;R为气体常数;ρ为密度;μ为动力粘度;λmean为平均孔隙直径。
5.根据权利要求2所述的一种射孔损伤室内综合评价方法,其特征在于,所述视渗透率校正系数kr的计算公式为:
kr=0.2379β2+4.4462β+1
式中:kr为视渗透率校正系数;β为流体校正系数。
6.根据权利要求1所述的一种射孔损伤室内综合评价方法,其特征在于,所述孔隙度分形维数D的计算公式为:
式中:D为孔隙度分形维数;λ为孔隙直径;λmax为最大孔隙直径;N为等效孔隙直径大于等于λ的孔隙个数。
7.根据权利要求1所述的一种射孔损伤室内综合评价方法,其特征在于,所述迂曲度分形维数T的计算公式为:
lgLT=(1-T)lgλ+TlgL
式中:T为迂曲度分形维数,λ为毛管直径,L为岩心长度,LT为直径λ的孔隙实际长度。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中海石油(中国)有限公司海南分公司;中海石油(中国)有限公司湛江分公司,未经中海石油(中国)有限公司海南分公司;中海石油(中国)有限公司湛江分公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202110520310.6/1.html,转载请声明来源钻瓜专利网。