| 一种颗粒样品两级渗透率的获取方法 | |
江文滨 ; 林缅; 周羁; 曹高辉; 姬莉莉
| |
| 2022-02-01 | |
| Rights Holder | 中国科学院力学研究所 |
| Abstract | 本发明实施例涉及一种颗粒样品两级渗透率的获取方法，针对现有的颗粒样品渗透率测试中发现的部分样品的压力衰减曲线具有快速和慢速两个明显不同阶段的特点，假定页岩内部具有两级渗流通道和空间，可分别提取两级空间的孔隙度和渗透率，所提取的孔、渗参数正演计算的结果较现有方法与实测值更为吻合，有助于深入认识页岩气储集空间孔、渗特征，为储层评价提供更全面和准确的参数。 |
| Application Date | 2019-09-30 |
| Application Number | CN201910942485.9 |
| Patent Number | CN110579433B |
| Claim | 1.一种颗粒样品两级渗透率的获取方法，其特征在于，所述方法包括： 获取压力随时间变化的曲线P(t)、基准腔体积Vb、放入样品前样品腔体积Vs、样品总体体积Vsb、颗粒平均半径Ra、测试温度T、测试阀开启前基准腔的气体密度ρ1和样品腔的气体密度ρ0； 基于氦气状态方程、压力衰减曲线P(t)最终稳定值P2和温度T，计算对应的气体密度ρ2、样品孔隙体积VP、孔隙度φ、颗粒样品外部空间和内部孔隙空间的体积比Kc和测试阀未开启时颗粒外部空间的气体平均密度ρc0； 确定气体开始进入样品对应的压力P3，从所述压力随时间变化的曲线P(t)中确定压力P3对应的时刻t3； 选择所有大于所述时刻t3的测点，基于氦气状态方程计算气体密度ρ(t)和质量分数FR(t)，对ln(FR(t))中所有大于所述时刻t3的测点和时间t进行线性拟合，得到斜率s1，根据所述斜率s1计算渗透率k0； 确定渗透率k′，设定i的初始值为4，以第i个测点为界，将压力衰减曲线分为两段：第i个测点及以前的测点记为A段，压力衰减曲线为PA(t)；第i个测点及以后的测点记为B段，压力衰减曲线为PB(t)； 基于氦气状态方程、第i个测点的压力Pi和温度T，计算对应的气体密度ρi，计算A段的孔隙体积VpA、孔隙度φA(i)、颗粒样品外部空间和内部孔隙空间的体积比KcA和测试阀未开启时颗粒外部空间的气体平均密度ρc0A，输入(kA(i)，D2A(i))＝f(PA(t)，KcA，ρc0A，ρ0，Ra，φA(i)，k0，k′)求得对应的渗透率最优解kA(i)和A段最小偏差平方和D2A(i)； 计算B段的孔隙体积VpB、孔隙度φB(i)、颗粒样品外部空间和内部孔隙空间的体积比KcB和测试阀未开启时颗粒外部空间的气体平均密度ρc0B，输入(kB(i)，D2B(i))＝f(PB(t)，KcB，ρc0B，ρ0，Ra，φB(i)，k0，k′)求得对应的渗透率最优解kB(i)和B段最小偏差平方和D2B(i)； 确定kA(i)和kB(i)为两级渗流通道的渗透率值。 |
| Language | 中文 |
| Classification | 发明授权 |
| Status | 有效 |
| Note | 授权 |
| Country | 中国 |
| Agency | 北京和信华成知识产权代理事务所(普通合伙) |
| Document Type | 专利 |
| Identifier | http://dspace.imech.ac.cn/handle/311007/90035 |
| Collection | 流固耦合系统力学重点实验室 |
| Recommended Citation GB/T 7714 | 江文滨,林缅,周羁,等. 一种颗粒样品两级渗透率的获取方法. CN110579433B[P]. 2022-02-01. |
| Files in This Item: | Download All | |||||
| File Name/Size | DocType | Version | Access | License | ||
| 20220201_0C_CN_0.pdf（689KB） | 专利 | 开放获取 | CC BY-NC-SA | View Download | ||
Items in the repository are protected by copyright, with all rights reserved, unless otherwise indicated.
Edit Comment