申请/专利权人：中国科学技术大学

申请日：2021-03-26

公开（公告）日：2024-03-29

公开（公告）号：CN113065246B

主分类号：G06F30/20

分类号：G06F30/20;G06F30/28;E21B49/00;G06F113/08;G06F119/08;G06F119/14

优先权：

专利状态码：有效-授权

法律状态：2024.03.29#授权;2021.07.20#实质审查的生效;2021.07.02#公开

摘要：本发明提供了一种基于射孔测试联作系统的地层参数反演方法，包括以下步骤：1根据射孔测试联作技术获得的一组井筒压力p和井筒温度T数据；2将所述压力p和所述温度T数据转换为井筒流量q数据；3利用图版拟合试井解释方法，将现场压力数据与理论曲线进行拟合匹配，反演获得地层参数，所述地层参数包括渗透率、井筒存储常数和表皮因子。本发明解决了常规试井无法获得准确非等温地层参数的问题，适用于非等温油气藏实际生产过程中的地层参数反演计算，大大缩短了试油周期，降低了试油成本。

主权项：1.一种基于射孔测试联作系统的地层参数反演方法，包括以下步骤：1根据射孔测试联作技术获得的一组井筒压力p和井筒温度T数据；2将所述压力p和所述温度T数据转换为井筒流量q数据；3利用图版拟合试井解释方法，将现场压力数据与理论曲线进行拟合匹配，反演获得地层参数，所述地层参数包括渗透率、井筒存储常数和表皮因子；步骤2中，所述压力p和所述温度T数据转换为井筒流量q数据的过程由以下推出：首先，SRK状态方程： 当为单组分气体时，Tr＝TTc， 当为多组分气体时，则必须考虑各组分气体之间的相互作用，此时的 其中，p为井筒压力，单位为Pa；T为井筒温度，单位为K；V为气体体积，单位为m3；R为气体常数，R＝8.314Jmol·K；ω为偏心因子；Tc代表临界温度，单位为K；Pc代表临界压力，单位为Pa；xi和xj为i组分和j组分的摩尔数，kij为i组分和j组分之间的二元交互系数，ai和aj则为假设各自为单组分气体通过方程2计算得到的；同理，bi也是假设各自为单组分气体通过方程3’计算得到；气体状态方程：pV＝zRT3其次，联合公式1和4，得到气体偏差因子z方程：z3-z2+Am-Bm-Bm2z-AmBm＝04 根据测量得到的所述压力p和所述温度T数据，根据方程5-7得到气体偏差因子z，再次，地层中气体流动的连续性方程： 其中，r为地层计算半径，单位为m；φ为孔隙度；μ为流体的粘度，单位为Pa·s；k为渗透率，单位为m2；Ct为总压缩系数，Ct＝Cg+Cr，Cg为气体压缩系数，Cr为岩石压缩系数，Cg和Cr单位均为Pa-1，然后，考虑表皮效应： 其中，pwf表示井底压力，单位为Pa；Skin表示表皮因子；t为时间，单位为s；r为地层计算半径，单位为m；rw为井筒半径，单位为m；然后，内边界条件： 其中，C为井筒存储常数，单位为m3Pa；外边界条件：pr→∞,t＝010初始条件：pr,t＝0＝pi11然后，将控制方程8-1214改写为拟压力形式13，拟压力定义： 其中，m为拟压力定义；从地层流入井筒的气体产量： 其中，h为地层厚度，单位为m，然后，采用如下的无量纲定义：无量纲时间无量纲拟压力无量纲半径无量纲井筒存储常数无量纲流量然后，将控制方程13-14改写无量纲方程组： 然后，对方程组15进行Laplace变换得到方程组16： 然后，求解方程组16，得到公式17： 然后，对Laplace空间上的气体产量进行Laplace数值反演，真实空间气体产量qD，由下述公式18-19Stehfest数值反演获得： 气体体积产量： 然后，根据10、14和15联合得到： C为井筒储存常数：C＝CgVw20Vw是井筒体积，单位为m3；Cg为气体压缩系数，单位为Pa-1： 然后，根据公式20-23，将测量得到的所述压力p和所述温度T数据，转换为井筒流量qt： 气体体积Vt：

全文数据：

权利要求：

百度查询： 中国科学技术大学 一种基于射孔测试联作系统的地层参数反演方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD