(19)中华 人民共和国 国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111505364.1 (22)申请日 2021.12.10 (71)申请人 西南石油大 学 地址 610500 四川省成 都市新都区新都大 道8号 (72)发明人 郭晶晶　王思博　王海涛　 (74)专利代理 机构 成都方圆聿联专利代理事务 所(普通合伙) 51241 代理人 邓永红 (51)Int.Cl. G06F 30/28(2020.01) G06F 111/10(2020.01) G06F 113/08(2020.01) G06F 119/14(2020.01) (54)发明名称 一种复杂情况下的致密油藏压裂井试井模 拟方法 (57)摘要 本发明公开了一种复杂情况下的致密油藏 压裂井试井模拟方法， 其包括： 构建同时考虑了 压裂缝的变导流效应、 储层的应力敏感 效应及储 层的径向非均质性的， 能同时准确模拟井底压力 和不同裂缝部位产液状况的致密油藏压裂井试 井模型。 本发 明可为压敏性复合致密油藏变导流 裂缝的开采提供更全面的、 能综合考虑多种因素 影响的不稳定试井模型和模拟方法。 权利要求书5页 说明书15页 附图4页 CN 114201932 A 2022.03.18 CN 114201932 A 1.一种复杂情况 下的致密油藏压裂 井试井模拟方法， 其特 征在于， 包括： 构建同时考虑了压裂缝变导流效应、 储层应力敏感效应、 储层的径向非均质性的致密 油藏压裂井试井模型， 通过所述试井模型获得井底压力变化和裂缝不同部位产液状况的模 拟结果； 其中， 所述变导 流效应是指压裂缝渗透率随位置不同而发生变化的现象； 所述应力敏感效应是指油藏 渗透率随压力的变化而变化的现象； 所述径向非均质性是指储层的内区渗透率、 孔隙度和综合压缩系数与外区的这些参数 具有差异的现象。 2.根据权利要求1所述的模拟方法， 其特征在于， 所述致密油藏压裂井试井模型的构建 包含以下设定条件： 该压裂井为压敏性复合致密油藏的开采井， 该压敏性复合致密油藏的储层中存在应力 敏感效应， 且 储层中油藏的渗透率随压力的变化遵循指数式变化 规律； 该压敏性复合 致密油藏所在的渗 流场为等温渗流场； 该压裂井压裂过程中形成的压裂缝存在所述变导 流效应。 3.根据权利要求1所述的模拟方法， 其特征在于， 所述致密油藏压裂井试井模型的构建 包括： 建立压敏性复合致密油藏的线汇井渗流模型， 其包括： 给定考虑所述应力敏感效应的 线汇井模型定解条件； 求解所述线汇井渗 流模型， 获得其摄动解； 基于所述摄动解， 采用叠加原理， 获得由所述压裂井引起的所述压敏性复合致密油藏 的压力响应解； 基于压裂缝内流动特征， 建立考虑所述应力敏感效应和所述变导流效应的压裂缝渗流 模型； 求解所述压裂缝渗流模型， 获得描述压裂缝模型中渗流的积分方程， 即积分形式的渗 流计算模型； 采用单元离散法、 分布积分法对所述渗流的积分方程进行求解， 获得用来刻画井底压 力和流量密度的解矩阵。 4.根据权利要求3所述的模拟方法， 其特征在于， 所述构建还包括： 对所述线汇井渗流 模型和所述压裂缝渗流模型进行无因次化处理， 获得其对应的无因次模型， 其后的求解进 一步基于所述无因次模型进行。 5.根据权利要求4所述的模拟方法， 其特征在于， 所述无因次模型的求解通过Pedrosa 变换、 摄动变换、 Laplace变换、 贝塞尔函数 理论、 二重积分、 积分方程单元离散、 分部积分中 的一种或多种实现。 6.根据权利要求 4所述的模拟方法， 其特 征在于， 其还 包括： 通过高斯消元法对所述解矩阵进行求解， 并通过Stehfest数值反演法将Laplace域的 解反演到真实时间域， 获得真实时间域的解； 进一步的， 将上述真实时间域的井底压力解代入Pedrosa变换， 获得考虑应力敏感效应 的真实时间域内的无因次井底瞬态压力。权　利　要　求　书 1/5 页 2 CN 114201932 A 27.根据权利要求4所述的模拟方法， 其特征在于， 其中， 所述压力响应解的获得包括： 利 用压降叠加原理， 对所述摄动解沿缝长进 行积分， 和/或， 所述渗流计算模型的获得包括： 对 所得压裂缝渗流模型的无 因次模型进 行求解； 求解中对压力的二阶导数项关于横坐标进 行 二重积分； 其后将所 得积分式与所述压力响应解联立， 得到积分形式的所述 渗流计算模型。 8.根据权利要求3所述的模拟方法， 其特 征在于， 所述 渗流计算模型的求 解包括： 获得所述渗流计算模型的离 散化模型； 将所述离散化模型与离散化的总流量与流量密度的关系式进行联立， 得到所述解矩 阵； 其中， 所述离 散化包括采用单 元离散和分布积分的处 理过程。 9.根据权利要求1 ‑8中任一项所述的模拟方法， 其特征在于， 所述致密油藏压裂井试井 模型包括以下模型中的一项或多 项： 油藏中的线汇井渗 流模型， 其包括： 油藏内区渗 流微分方程： 油藏外区渗 流微分方程： 内边界条件： 外边界条件： p2(r,t)|r→∞＝0     (4)， 内、 外区衔接面 流量连续条件： 衔接面压力相等条件： 初始条件： p1|t＝0＝p2|t＝0＝pi         (7)， 其中， p1为内区油藏压力， p2为外区油藏压力， r为油藏中某点距线汇井的径向距离， γ 为油藏渗透率的应力敏感系数， pi为原始地层压力， φ1为内区油藏的孔隙度， φ2为外区油 藏的孔隙度， Ct1为内区油藏的综合压缩系数， Ct2为外区油藏的综合压缩系数， μ为原油粘 度， k1i为内区油藏的初始渗透率， k2i为外区油藏的初始渗透率， e为自然指数， t为时间， ξ为 无穷小量， 为线汇井的产 量， B为地层原油体积系数， h为油藏厚度， π为圆周率， rf为内、 外区衔接面径向距离； 上述线汇井渗 流模型的无因次模型：权　利　要　求　书 2/5 页 3 CN 114201932 A 3