一种针对非常规储层样品的相对渗透率曲线计算方法 | |
李曹雄![]() ![]() ![]() ![]() ![]() | |
2020-03-17 | |
Rights Holder | 中国科学院力学研究所 |
Abstract | 本发明提供一种针对非常规储层样品的相对渗透率曲线计算方法,包括如下步骤:使用岩石成像技术对真实岩样进行扫描和识别后得到相应图像,将图像中岩样的矿物颗粒类型划分为组分特征,将矿物颗粒中的粒内孔,矿物颗粒间的粒间孔信息划分为孔隙特征,将岩样中的其它特征划分为辅助特征;利用混合分形方法分别拟合每种特征的分形参数;建立立方体构成的表征单元体,对各特征按照分形重构的块体大小和个数顺序放进表征单元体,建立连通概率矩阵,通过计算即可得到该模型的相渗曲线图。本发明使用混合分形单元法提取分形特征,能够快速计算非常规岩石样样品相渗曲线。 |
Application Date | 2018-05-31 |
Application Number | CN201810546641.5 |
Patent Number | CN108956415B |
Claim | 1.一种针对非常规储层样品的相对渗透率曲线计算方法,其特征在于,包括如下步骤: 步骤100,使用岩石成像技术对真实岩样进行扫描和识别后得到相应图像,将图像中岩样的矿物颗粒类型划分为组分特征,将矿物颗粒中的粒内孔,矿物颗粒间的粒间孔信息划分为孔隙特征,将岩样中的其它特征划分为辅助特征; 步骤200,将组分特征根据矿物颗粒的类型不同而细分为A1、A2……An个对应的子组分特征;将孔隙特征对应所述子组分特征的不同而细分为B1、B2……Bn个对应的子孔隙特征;将辅助特征根据与不同所述子组分特征的关系而细分为C1、C2……Cn个对应的子辅助特征;利用混合分形方法分别拟合每种子组分特征的分形参数; 步骤300,建立立方体构成的表征单元体,将表征单元体整体作为子组分特征A1,再将A2、A3按照分形重构的块体大小和个数顺序放进表征单元体,A2、A3的块体所占的几何空间的组分信息替换掉A1的组分信息,再将子孔隙特征B1随机放进表征单元体,统计每个孔的壁面的最大组成成分,作为该孔隙的组成成分; 步骤400,设定模型为油驱水,油水界面张力30mN/m,油粘度1.39mPa.s,水粘度1.05mPa.s,岩心表观渗透率1000mD,孔隙度22%,接触角0°,建立连通概率矩阵,通过计算即可得到该模型的相渗曲线图; 计算所述相渗曲线的公式为: 其中,每级孔对应的等效孔径为λ1,λ2,…,λi,…,λm,每级孔对应的个数分别为N1,N2,…,Ni,…,Nm;在排序后的连通概率矩阵中使用连通孔隙概率算法得到第i级连通概率P1,i;模型为非润湿相驱替润湿相,在第i级的Ni各孔中,驱扫概率为P2,i,使每级孔中驱扫成功的孔隙个数为P2,i×Ni,驱扫概率P2,i通过真实岩心压汞曲线和孔隙成像结果综合分析给出;第i级孔导通概率为P3,i,krw(j)为润湿相相对渗透率,kro(j)为非润湿相相对渗透率,Srw(j)为润湿相饱和度,Ci为第i级孔隙形状因子;其中1≤j≤m且j为整数; 当驱替动力大于不等径孔隙连接处阻力时,孔隙导通(P3,i=1),否则孔隙不能导通(P3,i=0),影响P3,i的参数包括毛管力、润湿性、驱替压力、孔隙等效直径和界面张力参数。 |
Language | 中文 |
Classification | 发明授权 |
Status | 有效 |
Note | 授权 |
Country | 中国 |
Agency | 北京和信华成知识产权代理事务所(普通合伙) |
Document Type | 专利 |
Identifier | http://dspace.imech.ac.cn/handle/311007/82232 |
Collection | 流固耦合系统力学重点实验室 |
Recommended Citation GB/T 7714 | 李曹雄,林缅,江文滨,等. 一种针对非常规储层样品的相对渗透率曲线计算方法. CN108956415B[P]. 2020-03-17. |
Files in This Item: | Download All | |||||
File Name/Size | DocType | Version | Access | License | ||
CN201810546641.5.pdf(764KB) | 专利 | 开放获取 | CC BY-NC-SA | View Download |
Items in the repository are protected by copyright, with all rights reserved, unless otherwise indicated.
Edit Comment