通过常规测井获取气藏克努森数的方法
【技术领域】
[0001] 本发明涉及石油勘探与开发，特别是涉及一种通过常规测井获取气藏克努森数的 方法。
【背景技术】
[0002] 在气藏流动特性研宄中，杨建和康毅力等在研宄气藏微观结构及渗流特征、许艳 芝和郑杰等在微尺度相似流动特性研宄、朱维耀等在纳微米级孔隙气体流动数学模型及压 裂直井产能分析以及李治平、李智锋和李勇明等页岩气藏纳米孔隙气体渗流特征分析等的 研宄都讨论了克努森数的应用。
[0003] 但如何获得克努森数，或者如何取得实际测井的应用等尚无进一步的研宄。而该 方法除应用于测井外，还能为油藏工程的其它方面提供有益的技术保障。但查阅国内外所 有相关的文献，目如未开展相关的研宄。
【发明内容】
[0004] 本发明的目的是为了克服上述【背景技术】的不足，提供一种通过常规测井获取气藏 克努森数的方法，为流体流动状态的判断提供了一种既切实可行又方便易行的思路。
[0005] 本发明提供的一种通过常规测井获取气藏克努森数的方法，包括如下步骤：步骤 一、读取常规测井数据；步骤二、测试并计算本地区地层参数；步骤三、利用步骤一与步骤 二获得的数据计算一个深度点的克努森数Kn;步骤四、一个深度点计算完成后，回到步骤 一，直到整个深度段的克努森数计算完成。
[0006] 在上述技术方案中，所述步骤一中，所述常规测井数据包括深探测电阻率RD和地 层孔隙度小。
[0007] 在上述技术方案中，所述步骤一中，所述地层孔隙度（i)通过包括密度、声波和中 子在内的孔隙度测井曲线计算而成。
[0008] 在上述技术方案中，所述步骤二中，所述本地区地层参数包括：阿尔奇胶结指数 m，阿尔奇饱和度指数n，阿尔奇常数a，地层水电阻率Rw，绝对温度T，压力P，汞的表面张力 〇，采在空气中的接触角Q，地层参数cl、c2和c3,Boltzmann常数kB以及气体分子直径 d〇
[0009] 在上述技术方案中，所述步骤三中，一个深度点的克努森数Kn的计算公式为：
【主权项】
1. 一种通过常规测井获取气藏克努森数的方法，其特征在于：包括如下步骤： 步骤一、读取常规测井数据； 步骤二、测试并计算本地区地层参数； 步骤三、利用步骤一与步骤二获得的数据计算一个深度点的克努森数Kn; 步骤四、一个深度点计算完成后，回到步骤一，直到整个深度段的克努森数计算完成。
2. 根据权利要求1所述通过常规测井获取气藏克努森数的方法，其特征在于：所述步 骤一中，所述常规测井数据包括深探测电阻率RD和地层孔隙度小。
3. 根据权利要求2所述通过常规测井获取气藏克努森数的方法，其特征在于：所述步 骤一中，所述地层孔隙度巾通过包括密度、声波和中子在内的孔隙度测井曲线计算而成。
4. 根据权利要求3所述通过常规测井获取气藏克努森数的方法，其特征在于：所述步 骤二中，所述本地区地层参数包括：阿尔奇胶结指数m，阿尔奇饱和度指数n，阿尔奇常数a， 地层水电阻率Rw，绝对温度T，压力P，汞的表面张力〇,汞在空气中的接触角0，地层参数 cl、c2和c3,Boltzmann常数kB以及气体分子直径d。
5. 根据权利要求4所述通过常规测井获取气藏克努森数的方法，其特征在于：所述步 骤三中，一个深度点的克努森数Kn的计算公式为：
式中： Kn-克努森数， RD-深探测电阻率， 小一地层孔隙度，m-阿尔奇胶结指数， n-阿尔奇饱和度指数， a-阿尔奇常数， Rw-地层水电阻率， T一绝对温度， P-压力， 0 -汞的表面张力， 9 一汞在空气中的接触角，cl、c2和c3 -地层参数， kB-Boltzmann常数， d-气体分子直径。
6. 根据权利要求5所述通过常规测井获取气藏克努森数的方法，其特征在于：所述步 骤三中，所述地层参数cl为2. 5?3. 2、c2为1. 11?1. 42、c3为3. 2?4. 7。
7. 根据权利要求6所述通过常规测井获取气藏克努森数的方法，其特征在于：所述步 骤三中，所述地层参数cl为2.8、c2为1.25、c3为3. 7。
8. 根据权利要求5至7中任一项所述通过常规测井获取气藏克努森数的方法，其特征 在于：所述步骤三中，所述采的表面张力〇为480dyne/cm，所述采在空气中的接触角0为 140。。
9. 根据权利要求5至7中任一项所述通过常规测井获取气藏克努森数的方法，其特征 在于：所述步骤三中，Boltzmann常数心为1. 3805*1(T23J/K。
10. 根据权利要求5至7中任一项所述通过常规测井获取气藏克努森数的方法，其特 征在于：所述步骤三中，所述气体分子为天然气分子，天然气分子直径didXKT1%，31 = 3.14159。
【专利摘要】本发明公开了一种通过常规测井获取气藏克努森数的方法，包括如下步骤：步骤一、读取常规测井数据；步骤二、测试并计算本地区地层参数；步骤三、利用步骤一与步骤二获得的数据计算一个深度点的克努森数Kn；步骤四、一个深度点计算完成后，回到步骤一，直到整个深度段的克努森数计算完成。本发明为流体流动状态的判断提供了一种既切实可行又方便易行的思路，可以广泛应用于石油勘探与开发领域。
【IPC分类】G06F19-00
【公开号】CN104573324
【申请号】CN201410776683
【发明人】张超谟, 李军, 张占松, 王晓畅, 张泽宇, 胡瑶, 张冲
【申请人】长江大学
【公开日】2015年4月29日
【申请日】2014年12月15日