(19)国家知识产权局 (12)发明 专利 (10)授权公告 号 (45)授权公告日 (21)申请 号 202111444127.9 (22)申请日 2021.11.30 (65)同一申请的已公布的文献号 申请公布号 CN 114169158 A (43)申请公布日 2022.03.11 (73)专利权人 中国石油大 学(北京) 地址 102249 北京市昌平区府学路18号 (72)发明人 田冷　柴晓龙　董鹏举　张瑞天　 王嘉新　黄灿　蒋丽丽　王恒力　 王泽川　 (74)专利代理 机构 北京三友知识产权代理有限 公司 11127 专利代理师 刘飞　许曼 (51)Int.Cl. G06F 30/20(2020.01)G06F 113/08(2020.01) G06F 119/14(2020.01) (56)对比文件 CN 111523696 A,2020.08.1 1 CN 113610292 A,2021.1 1.05 姬安召等.致密砂岩双重介质复合油藏偏心 井试井分析. 《石油化工 应用》 .2017,(第07期), 审查员 慈丽雁 (54)发明名称 一种致密砂岩气藏多层合采产能预测方法、 装置及设备 (57)摘要 本文提供了一种致密砂岩气藏多层合采产 能预测方法、 装置及设备， 所述方法包括： 获取致 密砂岩气藏中每个地层的地质参数； 根据每个地 层的地质参数， 结合预先建立的考虑井筒储集效 应和表皮效应的致密砂岩气藏单层产能渗流模 型， 计算得到每个地层在指定时间段内的致密砂 岩气藏单层产能渗流模型在拉式空间中的无因 次井底压力解； 根据每个地层在指定时间段内的 致密砂岩气藏单层产能渗流模型在拉式空间中 的无因次井底压力解， 以及地层叠加原理， 计算 获得致密砂岩气藏多层合采的无因 次产量， 本文 能够在考虑井筒储集效应和表皮效应的基础上， 准确全面快速的进行致密砂岩气藏多层合采产 能的预测。 权利要求书4页 说明书13页 附图3页 CN 114169158 B 2022.11.22 CN 114169158 B 1.一种致密 砂岩气藏多层合采产能预测方法， 其特 征在于， 所述方法包括： 获取致密 砂岩气藏中每 个地层的地质参数； 根据每个地层的地质参数， 结合预先建立的考虑井筒储集效应和表皮效应的致密砂岩 气藏单层产能渗流模型， 计算得到每个地层在指定时间段内的致密砂岩气藏单层产能渗流 模型在拉式空间中的无因次井底压力解； 根据每个地层在指定时间段内的致密砂岩气藏单层产能渗流模型在拉式空间中的无 因次井底压力解， 以及地层叠加原理， 计算获得致密 砂岩气藏多层合采的无因次产量； 其中， 所述地质参数至少包括渗流半径、 井筒气体粘度、 地层系数、 井筒储容能力和原 始地层压力； 其中， 所述预先建立的考虑井筒储集效应和表皮效应的致密砂岩气藏单层产能渗流模 型通过如下公式表示： 初始条件： GD(rD,0)＝GD， 内边界条件： 外边界条件： 其中， GD为开采时 的无因次压力， β 为地层的地层系数， w为地层的井筒储容能力， rD为无 因次渗流半径， tD为无因次 时间， ReD为无因次边界半径。 2.根据权利要求1所述的方法， 其特征在于， 所述根据每个地层的地质参数， 结合预先 建立的考虑井筒储集效应和表皮 效应的致密砂岩气藏单层产能渗流模型， 计算得到每个地 层在指定时间段内的致密砂岩气藏单层产能渗流模型在拉式空间中的无因次井底压力解， 包括： 针对每个地层进行如下步骤： 将每个地层的地质参数带入到预先建立的考虑井筒储集效应和表皮效应的致密砂岩 气藏单层产能渗流模型中， 并引入拉普拉斯变换函数， 得到在拉式空间下 的渗流模型以及 在拉式空间下的渗 流模型的压力通 解； 根据所述在拉式空间下的渗流模型的压力通解， 结合所述预先建立的致密砂岩气藏单 层产能渗流模型的内边界条件和外边界条件， 计算获得所述在拉式空间下的渗流模型的压 力通解的系数； 根据所述在拉式空间下的渗流模型的压力通解的系数， 得到所述地层在指定时间段内 的致密砂岩气藏单层产能渗 流模型在拉式空间中的无因次井底压力解。 3.根据权利要求2所述的方法， 其特征在于， 所述地层在指定时间段内的致密砂岩气藏 单层产能渗 流模型在拉式空间中的无因次井底压力解 通过如下公式表示：权　利　要　求　书 1/4 页 2 CN 114169158 B 2其中， 为在拉氏空间中单层开采的无因次井底压力解， β 为地层的地层系数， u为拉 氏变量， reD为无因次边界半径， GiD为单层开采的无因次原始地层压力， z＝ω/β， 为地层井 筒储容能力 ω与地层系数β 的比值， I0为第一类型的第一类零阶虚宗量的Bessel函数， I1为 第二类型的第一类零阶虚宗量的Bessel函数， K0为第一类型的第二类零阶虚宗量的Bessel 函数， K1为第二类型的第二类零阶虚宗量的Bes sel函数。 4.根据权利要求1所述的方法， 其特征在于， 根据每个地层在指定时间段内的致密砂岩 气藏单层产能渗流模型在拉式空间中的无因次井底压力解， 以及地层叠加原理， 计算获得 致密砂岩气藏多层合采的无因次产量， 包括： 根据所述致密砂岩气藏单层产能渗流模型和底层叠加 原理， 获得致密砂岩气藏多层合 采产能渗 流模型， 所述 致密砂岩气藏多层合采产能渗 流模型通过如下公式表示： 初始条件： pDj(rD,0)＝piDj， 考虑井筒储集效应 影响的内边界条件： 考虑表皮效应的内边界条件： 外边界条件： 其中， pD是无因次平均地层压力， pDj为第j层地层的无因次地层压力， tD为无因次时间， piDj为第j层地层的无因次原始地层压力， PwDj为第j层地层的无因次井底压力， rD为无因次 渗流半径， qsfDj为第j层地层的井底无因次流量， qscDj为第j层地层的地面无因次流量， wj为 第j层地层的井筒储容能力， βj为第j层地层的地层系数， Sj为第j层地层的表皮系数， ReD为 无因次边界半径， CDj为第j层地层的无因次井筒储集系数； 根据致密砂岩气藏多层合采产能渗流模型， 结合开采井井壁压力相同原理， 获得拉氏 空间中总的井底无因次流 量模型和拉氏空间总的地 面无因次流 量模型； 根据所述每个地层在指定时间段内的致密砂岩气藏单层产能渗流模型在拉式空间中 的无因次井底压力解、 拉氏空间中总的井底无 因次流量模型和拉氏空间总的地面无因次流 量模型， 计算获得在指定时间内致密砂岩气藏多层合采在拉氏空间中总的地面无因次流 量。 5.根据权利要求4所述的方法， 其特征在于， 在指定时间内致密砂岩气藏多层合采在拉 氏空间中总的地 面无因次流 量通过如下公式表示：权　利　要　求　书 2/4 页 3 CN 114169158 B 3